1. Introduction
Water plays a significant role in the survival and development of human society [
1,
2,
3,
4,
5,
6,
7]. Lakes, reservoirs, and rivers provide an important source of drinking water for people. Water quality is largely influenced by natural processes—for example, precipitation rate, weathering processes, and soil erosion, as well as anthropogenic influences, including industrial, urban, and agricultural activities [
8,
9,
10]. In recent years, the degradation of water quality has been regarded as one of the most serious global ecological issues, especially in developing countries [
11]. Therefore, the effective assessment of water quality is critical to ecological sustainability.
Water quality evaluation is a comprehensively multi-index decision-making process. The transformation of qualitative and quantitative concepts needs to be addressed, and the qualitative level should be quantified with different indices. A great deal of methods and methodologies have been proposed to determine water quality in the past few decades, including multivariate statistical models, fuzzy set-based models [
12], and hybrid models [
13,
14]. Multivariate statistical models for randomness have included cluster analysis (CA), factor analysis (FA), and principal component analysis (PCA) [
15,
16,
17]. The above methods take fuzziness and randomness less into account in the evaluation process. A new evaluation method is proposed based on the cloud model theory [
18] (, which considers randomness with fuzziness and is effective as a cognitive technology for water quality assessment. The cloud model makes it possible to travel from qualitative concept and quantitative characteristics. This theory has been widely applied to a great number of areas [
19,
20]. For example, in hydrology, Deng et al. [
21] proposed a hybrid forecasting method for water quality time series based on a cloud model and fuzzy forecasting. Wang et al. [
22] introduced the normal cloud model to evaluate the status of eutrophication.
Related studies have been carried out in lakes, reservoirs, and rivers [
20,
23,
24], but few studies have evaluated the water quality of springs. The spring is an important source of drinking water and one of the most famous tourist attractions in Jinan, China, where springs have attracted tourists from all over the world. In recent years, with the increase in the number of tourists, the water quality of the springs is worsening. There is a need to propose a comprehensive evaluation system to assess the water quality of springs. In the present study, the data of major springs are collected for 10 parameters at 10 different sites in Jinan. It is assumed that the evaluation of spring water quality has ambiguity and randomness. Uncertainties of variables and membership are graphically modeled via probability distributions and the membership function, respectively. A multi-index evaluation system of water quality for springs in Jinan is proposed by using normal cloud model theory.
2. The Normal Cloud Model
The cloud model theory is constructed on the idea of fuzzy sets theory [
25] and probability concepts [
26], proposed by professor Li [
27]. It is a new graphic model, which can synthetically picture the randomness and fuzziness of concepts, and completes the uncertain transformation between qualitative concepts and their quantitative instantiations [
28].
Let
be the universe of discourse and
a qualitative concept connected with
. If the quantitative number
x belongs to
, and
x is a random instantiation of the concept
, the certainty grade
for
x belonging to
is a random variable which has a stable tendency. The parameter
can be estimated as below:
Mathematically, the distribution of
, can be determined by three numerical characteristics: expectation parameter
, entropy parameter
and hyper entropy parameter
[
29,
30,
31]. According to
Figure 1,
is the expectancy of the cloud drop in the universe of discourse and the best characterization of a qualitative concept.
, the entropy of
, reflects the variation range of a cloud drop in the distribution, which is determined by the fuzziness and randomness of the concept. The degree of dispersion of the cloud drops is recognized by
, the entropy of
. In addition,
can measure the uncertainty of the membership degree and control the aggregation of cloud drops. For example,
Figure 1 illustrates a normal cloud model, in which the three numerical characteristics are
,
, and
.
The distribution of is defined as a cloud, and each x is called a cloud drop, which is a tool as a quantitative meaning to describe a qualitative concept. The production process of cloud drops represents the uncertainty of conversion between qualitative concept and quantitative values. The certainty grade of a cloud drop reflected the fuzziness and randomness of a concept.
If a number
is a random realization of the concept
and satisfies
and
, the certainty grade of
belonging to concept
, satisfies
then the distribution of
on the universal set
is called the one-dimension normal cloud, and each
is defined as a normal cloud drop.
The cloud generator is a cloud generation algorithm that is implemented in software. Cloud generators can achieve the mapping relationship between qualitative concepts and quantitative characteristics, and it is feasible to simulate the human cognition process using cloud generators by computer [
32]. There are two common types of generators: the normal cloud generator and the backward cloud generator.
Specifically, the normal cloud generators make it easier to realize the transformation from the qualitative concept to the quantitative characteristic. In contrast, the backward cloud generator can achieve the transformation from quantitative characteristic to qualitative concept. Using generators, as many cloud drops can be produced as are required. The normal cloud generator algorithm can be expressed in the algorithm as follows [
20]:
Input: three numerical characteristics of the cloud model (, , ) and the required number of cloud drops ;
Output: cloud drops and their certainty degree .
The algorithm steps can be presented as follows:
- (1)
Calculate expectation and the standard deviation ;
- (2)
Generate a random number : . The is a normally random number of which the expectation is En and standard deviation is He.
- (3)
Variable is a quantitative value of a qualitative concept—generate a normally random number with expectation and standard deviation :;
- (4)
Calculate certainty degree: ;
- (5)
Generate a cloud drop with the certainty degree and the normally random number ;
Repeat steps (1) to (5) until n required cloud drops are generated.
3. Study Area
Jinan is located in Shandong province, China, which is in the north of Mount Tai and south of the Yellow River, at latitude 36°32′–36°51′ N, longitude 116°49′–117°14′ E. There are more than 600 Springs, and the most famous of these springs are the four big springs: Baotu Springs, Heihu Springs, Pearl Springs, and Wulong Springs. As a result, Jinan is known as the “City of Springs”. The formation of springs is closely related to geological and hydrogeological conditions.
The geographical location and the geological conditions of springs in Jinan are shown in
Figure 2. Jinan is located at the conjunction of a mountain and a flood plain, which is northwest of Tai mountain and southeast of Lubei plain. Therefore, groundwater flow in this area is from southeast to northwest, which is controlled by geological layers and aquifers. Specifically, the main geological layers in this field are the Cambrian (Є), Ordovician (O), Quaternary (Q), Archaean (Ar), and the igneous rock (γ) layers. The main active faults are the Dongwu faults, Qianfoshan faults, Licheng faults, and Ganggou faults. Moreover, the springs are recharged by precipitation in the areas, with the Cambrian and Ordovician layers and leakage from the upper Quaternary aquifer. There are many reservoirs and rivers in this area, such as Jinxiuchuan reservoir, Wohushan reservoir, the Yellow River, and the Yufu river, the catchment area of which are up to 7000 m
2. Restricted by igneous rocks in the northwest part of the Springs field, groundwater flows upward as soon as the flows approach the north. Later, upward flows finally discharge to the ground, then the famous springs are formed. The Baotu springs, the Heihu springs, the Pearl springs, and Wulong Lake are the most famous springs in Jinan (
Figure 3). The area of the Baotu spring, the Pearl spring and the Wulong Lake are about 10,500 m
2, 1240 m
2, 54,000 m
2, respectively. However, with the increase of water demand in industry and agriculture, as well as the development of tourism, there is a declining trend in flow rates, and the pollution of springs is becoming more and more serious. To fully understand the pollution situation of springs field in Jinan, a brief assessment of water quality of springs in Jinan is included in this paper.
4. Water Quality Assessment of Springs
The flow chart of water quality assessment, based on a normal cloud model, is seen in
Figure 4. First, the evaluation criteria of water quality assessment are determined, then we divide them into five levels. Second, based on MATLAB software, three numerical characteristics are calculated, and cloud drops are built in a graphic way. Third, a comprehensive weight calculating algorithm is used to assign weights. Finally, reading the measured data and weights, the certainties of each spring are determined by a normal cloud generator. Then the level of water quality of springs can be determined.
4.1. Assessment Index System
Water quality is significantly affected by rainfall, human activities, and industrial impacts, as well as soil and geology. Water quality is mainly reflected in physicochemical characteristics, ecotoxicological characteristics, and the quantification of metal concentration [
22]. The evaluation criteria are classified into five levels by referring to the Standard of Natural Mineral Water for Drinking [
33] and the Hygienic Standard of Water for Drinking [
34] in China, as listed in
Table 1. Water whose level is I or II is safe for drinking, but water belonging to other levels must be comprehensively purified before drinking and irrigation.
To make the results correspond to the practical situation of springs, 10 parameters—total hardness, Pb, Cd, As, Hg, chloride, fluoride, sulphate, nitrate, and nitrite are adopted as evaluation factors. The evaluation factors and their classification standards are detailed in
Table 2.
4.2. Numerical Characteristics
Three numerical characteristics of the normal cloud model were determined according to the specific evaluation criteria (
Table 2) [
19,
22].
For a bilateral boundary of the form
, three numerical characteristics (
,
,
) can be determined as suggested in Equations (3)–(5):
Apart from that, if a variable has only a single boundary, like
or
, its default boundary parameter and
can be determined by the value of the upper and lower bounds:
In detail,
and
are the minimum and maximum values that can be accepted by a qualitative concept. For instance, for level III of total hardness,
and
, so the corresponding numerical characteristics can be calculated as
and
. Moreover,
is a constant (
k < 0.5), which controls the granularity of a cloud drop, and
can be adjusted appropriately according to the practical situation. In the present study,
k is assumed as 0.1 to balance the variation and robustness of water quality assessment. In that situation, the cloud parameters would be calculated by Equation (2), and all parameters are determined as in
Table 3.
4.3. Assessment Index Weight
Considering that the weight of each index is not absolute but relative, a comprehensive weight calculating algorithm is introduced to determine weights of indices. In particular, the analytic hierarchy process method (AHP) [
35] can fully draw on the experience of experts, and it has been widely used to weight factors based on multiple criteria. Therefore, this paper employed the AHP method for the subjective weight calculation. The projection pursuit classifications (PP) [
36,
37,
38] are adopted to calculate the objective weight, to reduce human factors and balance the subjective uncertainty and randomness of AHP. The PP method is a statistical method for analyzing and dealing with high-dimensional data. It can project high-dimensional data to low-dimensional subspaces, and the characteristics of structure of high-dimensional data can be studied based on analyzing the scattered structure of low-dimensional projected data. Furthermore, the relevancy between evaluation factors and water quality could be determined quantitatively by the comprehensive weight calculating algorithm. The proposed algorithm is a practical method, which could reconcile subjectivity and objectivity in the weight calculating.
The basic calculating steps are as follows.
To calculate the subjective weight by the AHP approach, a hierarchy structure model is constructed on the basis of an assessment index system, and is then scored by experts. Next is setting up a judgement matrix, ranking evaluation indices of water quality and completing the consistency checking, and calculating the total weights
[
39]. Finally, the relative importance of 10 indices is identified from high to low as Hg, As, Cd, Pb, Sulphate, Nitrate, Chloride, Fluoride, Nitrite, and Total Hardness (
Table 4). The analytic hierarchy process is used to calculate the subjective weight as the relative importance of the indicator, and thus the sum of the subjective weights may not be 1.
To calculate the objective weight, in view of the different evaluation units, the normalization method is used to make indices dimensionless. After normalization, the variation tendency of the primitive value is maintained. The indices are normalized using Equation (8):
In the formula, is the i index of the j samples for the water quality, i = 1, 2, …, n; j = 1, 2, …, m.
Data were projected from high-dimension to low-dimension by linear projection. Let
be the
m dimensional unit direction vector—then one-dimensional projection eigenvalue
(
i = 1, 2, …,
n) can be calculated by Equation (9):
After that, the objective function should be constructed as
In the formulas,
Sz is the standard deviation of the projection eigenvalue, which reflects the distances between different classes.
Dz is the local density of the projection value.
is the mean value of
;
R is the local density radius, and
k is assumed as
m;
is the distance between
and
; and
;
is a unit speed function, which can be calculated as
When Equation (10) achieve the maximum value, the projection direction can best reflect the data characteristics:
To solve this nonlinear optimization problem, a real-coded accelerating genetic algorithm (RAGA) was used. Considering
and
, the final weight
of each level was determined with the mean of
:
According to the principle of minimum entropy [
40], the hybrid AHP-PP weights are given as the formula below:
4.4. Generation of Normal Cloud Model
The normal cloud model of water quality of springs can be set up by the normal cloud generator. The certainty grade of each level of water quality assessment is calculated by the normal cloud generator.
Then, combining with the hybrid AHP-PP weights, the integrated certainty degree
M can be calculated as the following:
In the formula, is the certainty grade of an assessment level, and is the weight of index .
5. Results and Discussion
Water quality assessment of springs in Jinan is of great uncertainty, and randomness and fuzziness are the two most important uncertainties. In other words, accurate intensions and exact boundaries of the extensions are hard obtain in the evaluation process. It is essential to introduce the cloud model to discuss the randomness and fuzziness in water quality assessment.
In the present study, ten springs in Jinan were chosen as the study area. Some engineers and scientists were employed by the environmental protection department of Jinan to collect water quality data. The water quality data collected by the environmental protection department of Jinan are given in
Table 5 [
41,
42].
After completing the weight calculation, the distributions of certainty degrees of five grades are calculated and illustrated in
Figure 5. In
Figure 5, the abscissa represents the actual value of evaluation index, and the ordinate represents the certainty grade. Therefore, the more the cloud drops, the better the evaluation effect.
Based on the cloud model method, the process can be carried out as follows:
- (1)
First, we determined the evaluation criteria of water quality assessment, and divided it into five levels. Then the numerical characteristics can be calculated by Equations (3)–(7).
- (2)
In view of the comprehensive weight calculating algorithm, the weights of indices are given in
Table 4.
- (3)
Based on the obtained normal cloud model method and the weights, the certainty degree can be discerned by the normal cloud generator.
- (4)
The final evaluation results of ten springs are presented in
Table 6.
The reliability of the normal cloud model-based method was assessed by comparison with the single factor index (SFI) method [
43], the Nemerow index (NI) method [
44], the variable fuzzy sets (VFS) method [
24], and the artificial neural networks (ANN) method. The results of these methods are compared in
Table 7.
In the SFI method, the measured data of a spring are compared with water quality standards; then, the highest pollution grade of water quality was selected as the grade of this spring. In the NI method, water quality is assessed by calculation of comprehensive pollution index. In the variable fuzzy sets model, a comprehensive evaluation vector is established by performing an operation between the set of fuzzy weights and the fuzzy relationship matrix. Lastly, an artificial neural network is established with three layers in the ANN model, which represent independent and dependent variables of the system, respectively. Specifically, the input layer has ten neurons and the output layer has five neurons in the present study.
In
Table 6, each spring can be generally classified by five levels from low to high. The final levels of Baotu Spring, Heihu Spring, Baishi Spring, Jiunv Spring, Pearl Spring, Jiujiao Spring, Begonia Garden, and Wulong Lake are I; most of Wuyou Spring and Suyu Spring are level II. It can be concluded from
Table 7 that the calculated results obtained by various methods are basically the same, particularly for Wuyou Spring, whose final grade is II. The 90% (9/10) final levels are the same as those of the variable fuzzy sets method and the artificial neural networks method, and the 80% (8/10) results are exactly same as those of the Nemerow index method. The 20% (2/10) results are exactly same as those of the single factor index method. The data in
Table 6 show high consistency of the results between the proposed normal cloud model and other four methods, considering that the SFI method puts emphasis on the most polluting indicator, which is one-sided. In brief, the normal cloud method is highly representative in these methods.
It is worth noting that the cloud model not only achieves transformation between qualitative concepts and quantitative characteristics, but also provides more details about the water quality levels. For example, the total hardness of Baotu Spring was 285 mg/L, which should belong to grade II according to
Table 2. In addition, the certainty degrees of total hardness obtained by the normal cloud generator are
,
,
,
and
. Therefore, the certainty degree by quantitative analysis is
, and the total hardness of Baotu Spring more likely belongs to grade II, likely belongs to grade III, and almost impossibly belongs to grades I, IV, and V. The conclusion is consistent with qualitative analysis. Furthermore, the level of Pearl Spring is more likely to be I than that of Jiunv Spring and Jiujiao Spring, because the certainty degree for level I of Pearl Spring (0.9266) is higher than that of Jiunv Spring (0.7485) and Jiujiao Spring (0.8975). In other words, the results based on the normal cloud model not only reflect the water quality level accurately, but also further judge the pollution degree of water quality for different springs on the same level.
Among all the springs in Jinan, Baotu Springs are the most visited. Baotu Spring, Wuyou Spring, and Suyu Spring hosted an average of more than 20,000 visitors a day in the tourist season, up to a maximum of 50,000 during the Spring Festival and Lantern Festival. Moreover, there are some new toilets and a new subway being built near the Wuyou Spring and the Suyu Spring. The most likely causes of pollution for these springs are anthropogenic contamination and the recharge of groundwater during the construction of subways. Therefore, the evaluation results of the proposed method are basically consistent with the actual situation. In general, the water quality levels of springs in Jinan are mostly grade I, which is suitable for drinking.
6. Conclusions
A multi-index evaluation system is proposed to assess the water quality of springs based on normal cloud theory in the present study. Required cloud drops were generated based on three numerical characteristics calculated by cloud generator algorithm. The weights of indices are rationally assigned by a comprehensive weight calculating algorithm, which could reconcile subjectivity and objectivity in the weight calculating. Finally, the certainty of each spring belonging to each level can be calculated by a cloud generator, and the levels of water quality of springs are determined.
The present method is applied to the water quality assessment of springs in Jinan. The results are compared with those from other four methods: the SFI method, the NI method, the VFS method and the ANN method. The final levels of Baotu Spring, Heihu Spring, Baishi Spring, Jiunv Spring, Pearl Spring, Jiujiao Spring, Begonia Garden, and Wulong Lake are all level I. The final levels of Wuyou Spring and Suyu Spring are level II. The calculated results obtained by various methods are almost the same, particularly for Wuyou Spring. In other words, the results of the proposed normal cloud model are highly consistent with that of the other four methods. Moreover, the present method can not only reflect the water quality level accurately, but also further judge the pollution degree of water quality for different springs by comparing the certainty of the same level.
The normal cloud model provides a practical guide for the water quality evaluation of springs. In the present study, the one-dimensional normal cloud model is introduced into the water assessment of springs. In view of the successful application of the present model, a two-dimensional and a multi-dimensional cloud modeling strategy are worth exploring with regard to water assessment or other further applications.