Abstract

For earth-rock dams influenced by rainstorms, seepage status monitoring is very important and provides the basis for the safe and effective operation of earth-rock dams. The most influential factors concerning the seepage of earth-rock dams are the reservoir water level, precipitation, temperature, and timeliness, and the influence of the reservoir water level and precipitation on the seepage of an earth-rock dam exhibits hysteretic effects. The reservoir water level of an earth-rock dam abruptly increases and may exceed the historically highest water level, therein causing new deformations of the earth-rock dam or even plastic deformation. Thus, the permeability coefficient for parts of an earth-rock dam changes, and we present the exceeded water level factor. Considering the complexity of the seepage monitoring of earth-rock dams, based on the hysteretic reservoir water level and precipitation, temperature, timeliness, and the exceeded water level factor, a statistical model based on an explicit function and an artificial wavelet neural network model based on an implicit function are established. Based on these two models, an integrated monitoring model based on maximum entropy theory is established. At the end of this paper, three monitoring models are used for the seepage monitoring of a measuring point of an earth-rock dam influenced by rainstorms, and the results show that the three monitoring models obtain satisfactory predication accuracy.

1. Introduction

Due to the characteristics of low cost, fine environmental adaptability, and lower construction difficulty, earth-rock dams have been widely used and rapidly developed in the world, which make up over 80 percent of all dams [1]. During the construction and running period of earth-rock dams, safety problems such as seepage [2], cracks [3], and landslide [4] may happen. Seepage has considerable influence on earth-rock dams and often increases from small range to large range, which may cause dam settlement, collapse, and concentrated leakage passage for earth-rock dams. The structural damage may be a single form or multiple forms of damage in one part or different parts of the earth-rock dams. For the earth-rock dams in coastal areas suffer from rainstorms, the rapidly increased reservoir water lever and large amount of rain within a short time may cause threatened structural problems. Therefore, it is of great significance to study the seepage status of earth-rock dams influenced by rainstorms. Because earth-rock dams and the surrounding environment are rather complex and fickle, the potential seepage diseases are difficult to find out. Seepage monitoring analysis [5] of earth-rock dams could help to judge the existence of seepage damage and grasp the running status of earth-rock dams, which provides basis for the safety running of earth-rock dams.

To monitor and analyze the seepage status of earth-rock dams influenced by rainstorms accurately and timely, seepage monitoring models should be built up to help to find out seepage diseases conveniently and ensure the stable operation of the dams [5, 6]. The influence factors concerning earth-rock dam seepage are the reservoir water level, precipitation, temperature, timeliness, and so forth. In fact, the influence of the reservoir water level and precipitation on earth-rock dam seepage exhibits hysteretic effect. Moreover, rainstorms may lead the reservoir water level to exceed the historical highest water level, which causes new deformation or even plastic deformation of earth-rock dams. Then, the permeability property of earth-rock material would change, which influences the seepage state of earth-rock dams. The exceeded reservoir water level factor is put forward to consider the abrupt increased reservoir water level.

Therefore, the seepage statistical model considering the hysteretic effect of reservoir water level and precipitation and the exceeded water level factor is established. Meanwhile, the influence factors on seepage monitoring indexes are rather complex, which make it difficult to simulate with explicit function. Artificial wavelet neural network [79] is the implicit function, which has the advantage to explain complex relationship. Thus, earth-rock dam seepage monitoring model based on artificial wavelet neural network is established. Based on the two monitoring models and maximum entropy theory [1012], an integrated seepage monitoring model is set up to optimize the earth-rock dam seepage monitoring further. Finally, the three seepage monitoring models are applied to analyze the seepage status of an earth-rock dam influenced by rainstorms. The results show the three models with fine precision successfully used in earth-rock dam seepage monitoring, which provide technical support for seepage monitoring of other earth-rock dams.

2. Earth-Rock Dam Seepage Statistical Model considering the Hysteretic Effect of Reservoir Water Level and Precipitation and the Exceeded Water Level Factor

The most influential factors concerning earth-rock dam seepage include reservoir water level, precipitation, temperature, and timeliness. Given that the effect of reservoir water level and precipitation on seepage has the hysteretic effect, in the traditional seepage statistical model, the previous reservoir water level and precipitation are categorized based on averages over a number of previous days [5], for example, the previous two days, the previous five days, and the previous ten days. Practice has proven that the influences of reservoir water level and precipitation on seepage rise in the first stage and then decrease, which are presented as normal distribution [5]. The normal distribution curve is used to simulate the hysteretic effect, and the hysteretic days and influence days of reservoir water level and precipitation influence on earth-rock dam seepage are hard to determine. Considering the efficiency of general calculation methods is low, quantum genetic algorithm [13, 14] is used to calculate the hysteretic days and influence days to obtain the optimal earth-rock dam seepage statistical model. Moreover, due to the influence of rainstorms, the reservoir water level may exceed the historical highest water level. As a result, the abrupt increased reservoir water level is difficult to simulate. Hence, the exceeded water level factor is added.

2.1. Reservoir Water Level Component

For the seepage index, we take the piezometric tube level as an example. The piezometric tube level is hysteretic influenced by the reservoir water level, and the seepage index at the time is given as follows:where is the retardation time, , , and and are the reservoir water level at the corresponding time.

Equation (1) reflects the hysteretic relationship between the reservoir water level and the piezometric tube level, and the piezometric tube level at the time is continuously influenced by the previous reservoir water level. Because the hysteretic time is difficult to determine, the previous reservoir water level is often categorized based on averages over a number of previous days, for example, the previous two days, the previous five days, and the previous ten days. However, these factors are fuzzy and are unable to accurately reflect the hysteretic effect of reservoir water levels. Moreover, the influence of the reservoir water level on the piezometric tube level may be from the previous two days, the previous ten days, or longer.

Suppose that the seepage index is influenced by the reservoir water level from the previous days , and the equivalent water level is expressed as follows:where , is the weight of the th water level’s influence on the equivalent water level and is the function reflecting the hysteretic influence of the reservoir water level on the seepage index.

Considering the characteristic of the weight vector , the following equation is then obtained:

Therefore, the equivalent water level is obtained as follows:

Numerous studies have shown that presents a normal distribution [5]. Combined with the characteristics of normal distributions, the hysteretic days and influenced days are set as and . The process that describes the influence of the reservoir water level on the seepage index is shown in Figure 1, and the hysteretic influence function is given as follows:

If a water storage starts at the time , then, at the fixed time , the following equation is obtained:

For the piezometric tube level at a specific measuring point during a certain period of time, the hysteretic days and influenced days are constant.

in (5) can be treated as a function reflecting the influence of the previous water level on the equivalent water level . Suppose that the equivalent reservoir water level at the time is and that the influenced weight function . Then, the following equation is obtained:where is the undetermined hysteretic reservoir water level days, is the undetermined influenced normal distribution standard deviation, and is the reservoir water level at the time .

Quantum genetic algorithm is used to obtain the hysteretic days and influenced distributed parameter. Generally, the continuous integration is replaced by a discrete integration, and an integration range with the value of 3-4 times meets the accuracy requirements.

2.2. Precipitation Component

To reasonably consider the hysteretic effects of precipitation on earth-rock dam seepage, a lognormal distribution function is used to depict the precipitation effect on the earth-rock dam seepage. Specifically, the hysteretic influencing function of precipitation is given as follows: where is the hysteretic days of precipitation and is the influencing distribution parameter.

Suppose that the observation day is and that the starting day is . The influence range can be 1-2 months, and then,

Similarly, an exponential transformation is adopted for the influence of precipitation infiltration. Suppose that the equivalent precipitation at the time is . The hysteretic influenced function of precipitation meets the following condition:where is the precipitation at the time of , is the infiltration transformation index, , and the remaining symbols have the same meaning as in (8).

Therefore, the precipitation component is expressed as follows: where is the precipitation component, is the regression coefficient, is the undetermined hysteretic days of precipitation, is the undetermined distributed parameter of precipitation, is the precipitation at the time of , and is the equivalent precipitation.

Similarly, quantum genetic algorithm is used to optimally obtain and .

2.3. Temperature Component

Although seepage has no direct relation with the water temperature, air temperature, and batholith temperature, the temperature variation changes the viscosity of water and the batholith crack opening: as a result, the seepage status may change. The temperature component of seepage is mainly caused by the temperature variations in the dam foundation and the batholith on both sides of the dam, and the temperature variation exhibits annual periodic variations or semiannual periodic variations. However, sufficient measured temperature data have not been provided. Therefore, a temperature component consisting of sine and cosine functions with periods of a year and half a year is given as follows:where corresponds to the period of a year and half a year, , and and are the regression coefficient.

2.4. Timeliness Component

During operation of the earth-rock dam, the structure of soil grains may change. Meanwhile, the gradually produced deposit in front of the dam forms a natural blanket. The influence of these factors on the seepage exhibits the timeliness process and can be modeled using the following equation:where and are the regression coefficients of the timeliness component and is the cumulative number of days from the starting day to the measured day dividing 100.

2.5. Exceeded Eater Level Factor

For earth-rock dams influenced by rainstorms, the reservoir water level may exceed the previous historical highest water level. The reservoir water level is unable to accurately reflect abrupt increases in the reservoir water level. If the reservoir water level for a certain day exceeds the previous historical highest water level, the exceeded water level factor is used. The factor relates to the excess of the reservoir water level, the reservoir water level on the day of operation, and the rate of change of the water level. The greater the excess reservoir water level is, the greater the influence on the earth-rock dam seepage will be. Therefore, the exceeded water level factor is given as follows:where is the excess reservoir water level, is the product of and the water level on the measuring day, and is the rate of change over the previous days.

2.6. Earth-Rock Dam Seepage Statistical Model

By studying the earth-rock dam seepage statistical model considering the hysteretic effect of precipitation and reservoir water level and the exceeded water level factor, combined with the temperature component and timeliness component, the earth-rock dam seepage statistical model is obtained as follows:where is a constant term and the other parameters have the same meaning as mentioned above.

By selecting multiple correlation coefficients or the residual standard deviation as the objective function, the quantum genetic algorithm is used to obtain the optimal hysteretic parameters and the coefficient in (15). The calculation procedure is shown in Figure 2.

3. Earth-Rock Dam Seepage Monitoring Model Based on Artificial Wavelet Neural Network

The relationship between environment factors and earth-rock dam seepage is rather complex, and it is difficult for an explicit function to simulate this complex relationship. Artificial neural network algorithm is a type of implicit function with strong nonlinear fitting ability and adaptability. The wavelet transformation has the ability to better analyze the local details of the measured data and reflects the characteristics of the measured data. Therefore, combined with wavelet analysis and the artificial neural network algorithm, an earth-rock dam seepage monitoring model based on the artificial wavelet neural network is established.

3.1. Wavelet Theory

Wavelet analysis [15, 16] is a type of multiresolution data analysis method with the ability to analyze any details of an object. Through wavelet analysis, a signal is resolved into different frequency bands. Supposing that the wavelet function transforms with a magnitude of , the inner product between the original signal and the wavelet function using a scaling of is obtained as follows:

Through the transformation of the wavelet basis function and the analysis of the local features of the signal, the local characteristics of the seepage data are obtained, which help to reflect the multiscale change law of the measured data.

3.2. Earth-Rock Dam Seepage Monitoring Model Based on Artificial Wavelet Neural Network

Because it fuses the artificial neural network algorithm and wavelet theory, artificial wavelet neural networks have the advantage of providing multiscale analysis and implicit function. Therefore, the seepage monitoring model based on the artificial wavelet neural network more objectively reflects the seepage status. The transfer function of the nodes in the hidden layer of the neural network is the wavelet basis function. The weight value and threshold value are adjusted through error back propagation. Figure 3 shows the topological structure of the artificial wavelet neural network.

In Figure 3, are the input factors, which correspond to the reservoir water level component, precipitation component, temperature component, timeliness component, and the exceeded water level factor. is the output value, that is, the seepage index, and and are the weight values.

The data series of input factors are , and the output in the hidden layer iswhere is the output value of the th node in the hidden layer, is the weight between the input layer and the hidden layer, is the wavelet basis function, is the shift factor of the wavelet basis function, and is the scaling factor of the wavelet basis function.

As the excitation function in the network, the selection of wavelet function is very important to the fitting and predicated results. Since Morlet wavelet has the characteristic of calculation stability, small error, and fine robustness on error interference, a Morlet wavelet [17] with the function image shown in Figure 4 is used in the artificial wavelet neural network, and its specific formula is given as follows:

Then, the output of the seepage index is where is the exporting seepage index and the other parameters are similar to those in (17) and (18) and Figure 4.

The gradient descent method is used to modify the weights in the artificial wavelet neural network, through which the output values gradually approach the expectation until the termination condition is satisfied. The concrete steps are given as follows:(1)Parameter initialization: the scaling factor , the shift factor , and the weights and are initialized randomly, and the learning rate and the momentum coefficient are selected appropriately.(2)Classification of original samples: the original samples are classified as either the training set or the test set. The training set is used to train the network, and the test set is used to test the network.(3)Comparison of output values: the test set is imported into the network, and the mean squared error is obtained by subtracting the predicted output from the expected output as follows:where is the expected output and is the sample size.(4)Modification of weight: based on the error , the gradient descent method is used to modify the weight, scaling factor, and shift factor to make the predicted output approximate the expected output.(5)Convergence judgment: when the computation is converged, the computation is halted. Otherwise, the computation should return to Step (3), and the network parameters should be adjusted as follows:where , , and are the changes in the network parameters in the iteration.In addition, where , , and are the network parameters in the th iteration.

The weight and valve values are modified until the terminal condition is satisfied. The finally obtained parameters are substituted into the artificial wavelet neural network and combined with the import of reservoir water level component, precipitation component, temperature component, timeliness component, and the exceeded water level factor; the seepage index series are then exported. Finally, the monitoring model based on the artificial wavelet neural network is obtained and can be used to monitor earth-rock dam seepage.

4. The Integrated Earth-Rock Dam Seepage Monitoring Model Based on Maximum Entropy Theory

The two above-mentioned earth-rock dam seepage monitoring models mentioned above have their own strengths and weaknesses. Through the consideration of the hysteretic effect of seepage and the exceeded water level factor, the statistical model reflects the seepage status of the earth-rock dam to a certain extent. However, a reasonable explicit function that reflects the complex influencing factors of seepage is difficult to select. Artificial wavelet neural network helps to address this difficulty and is quite capable of determining the local details of measured data. Simultaneously, the artificial wavelet neural network suffers from artificial factors. Therefore, based on the application of maximum entropy theory, to obtain the advantages of the two above-mentioned monitoring models mentioned, the integrated earth-rock dam seepage model is established. This model is of great significance to the monitoring of the seepage status of earth-rock dams.

Maximum entropy theory originates from classical thermodynamics, which is an important theory in the frontier of modern physics. Entropy is applied to information theory, which has been successfully applied to different subject areas. Specifically, the expression of entropy is given as follows:where is the probability that the signal appears in a signal source and is the magnitude of the entropy, which measures the uncertainty of the system status.

Under the given condition, to solve an ill-posed problem, one probability distribution among all possible probability distributions is found to have the maximum value of entropy.

For the seepage index of earth-rock dams, the optimization problem in maximum entropy theory is given as follows:

In addition, the constraint condition is

When the entropy is maximized, the following equation holds:where is the entropy of the seepage index, is the probability when the value of the seepage index is , is a function, such as the first-order central moment and second-order central moment, of , is the mean value of , are Lagrange multiplier, corresponds to the earth-rock dam seepage statistical model, and corresponds to the earth-rock dam seepage monitoring model based on the artificial wavelet neural network.

The seepage index of earth-rock dams can be regarded as a discrete series continuously changing time, which is used in the two monitoring models to obtain predicted values and the difference between the predicted value and the measured value. With the obtained results regarding the constrained information of the predictor, maximum entropy theory is applied to solve the constraint information problem to improve the precision of the model predication. The concrete calculation steps are given as follows.

(1) Calculation of Model Eigenvalues. To obtain information about the earth-rock dam seepage statistical model and the earth-rock dam seepage monitoring model based on the artificial wavelet neural network, eigenvalues of the computation results of the models are obtained. Suppose that measured values are provided and the measured data set is , and the calculated values of the seepage index based on the two monitoring models are . The variance of the calculated values based on the two models is

(2) Solution of the Probability Density Function Based on Maximum Entropy. The earth-rock dam seepage index is a discrete random variable, and the following equation based on maximum entropy theory is given:

The constraint condition is

Substituting (26) into (29),

Combining (30) and (31), the Lagrange multipliers and are obtained, and the probability function of the th measured value is

(3) Prediction of Seepage Index. For the seepage index in the upcoming days, the predicated results of the th () seepage index based on the two monitoring models are , and the probability of the th predicated seepage index is

By integrating the probability, the desired value, which is the predicated value on the th day, is obtained. Then, repeating Step (3), the predicated series of the seepage monitoring values of the earth-rock dam is obtained.

5. Case Study

5.1. Project Profile

To verify the three models applied to the monitoring of the earth-rock dam seepage status, we take a reservoir located in Zhejiang, China, for analysis. The reservoir is an integrated large reservoir, which has the functions of flood protection, water supply, irrigation, and electricity generation. It controls drainage area of 258 square kilometers and has a total storage of 114 million cubic meters. The normal water level and the maximum flood level are 115.127 m and 123.527 m. The project consisted of diversion dam, normal spillway, emergency spillway, flood diversion sluice, water supply tunnel, and power station. The reservoir dam is Clay Core Wall Sand Dam, with the top elevation of dam, top elevation of wave wall, maximum dam height, width of dam crest, and length of dam crest being 124.527 m, 125.727 m, 36.9 m, 6 m, and 560 m, respectively. The upstream slope and downstream slope of the dam are 1 : 2.0–2.80 and 1 : 1.85–2.0.

Every year, there is abundant rain around the reservoir. Sometimes the reservoir suffered from rainstorm. In the period between 1 January 2012 and 31 December 2013, the reservoir suffered from rainstorm due to the influence of Typhoon Haikui which happened in August 2012 and Typhoon Fitow which happened in October 2013, and the maximum precipitation in one day is 137.2 mm and 183.0 mm. The corresponding increased reservoir water levels are 8.43 m and 7.89 m. After the rainstorm, the reservoir water level reduces gradually.

Therefore, we take the measuring point I in the cross section of 0–60.0 of the earth-rock dam for analysis, and the measured interval of measuring point I is 7 days. By setting the period between 1 July 2012 and 31 October 2013 as the modeling period, the three monitoring models are used to fit the piezometric tube level of the measuring point and predicate the seepage status of the measuring point between 1 November 2013 and 31 December 2013. Figures 5 and 6 show the monitoring data for precipitation and reservoir water level between 1 January 2012 and 31 December 2013, and their measured interval is 1 day.

5.2. The Application of Three Models in Monitoring the Seepage Status
5.2.1. Statistical Model considering the Hysteretic Effect of Reservoir Eater Level and Precipitation and the Factor of the Exceeded Water Level

Based on (15) and on the fitness function with the multiple correlation coefficients , a quantum genetic algorithm with a population of 40 is used to optimize the hysteretic parameters of the reservoir water level and precipitation. Figure 7 illustrates the iterative curve of the quantum genetic algorithm, and the optimal multiple correlation coefficient is 0.966. Table 1 shows the statistical metrics of the statistical model.

Through the solution obtained using the optimization algorithm, the hysteretic parameters , , , and are 19 days, 16 days, 40 days, and 40 days, respectively. The equivalent reservoir water level and precipitation based on (7) and (11) are shown in Figures 8 and 9. Figure 10(a) illustrates the comparison between the fitting values and the measured values, and Figure 10(b) shows the scatter plot of measured versus fitted with regression line. As shown in Figures 10(a) and 10(b), we can clearly see how model behaves in lower and higher regions. Table 2 shows the coefficients of the model.

For the prediction series between 1 November 2013 and 31 December 2013, the predicted values and measured values of piezometric level are shown in Table 3.

5.2.2. Seepage Monitoring Model Based on Artificial Wavelet Neural Network

The input of the monitoring model based on the artificial wavelet neural network includes the equivalent reservoir water level, equivalent precipitation, temperature, timeliness, and the exceeded water level factor, and the output is the measuring point I series. According to (15), 13 factors are input factors. By setting the number of network nodes of the hidden layer as 10, the topological structure of the network is 13-10-1.

Since the measured interval of measuring point I is 7 days, we select 7 days, 14 days, and 21 days as the sampling time to have the comparative analysis. Table 4 shows the mean squared error for the three cases, and the artificial wavelet neural network model with the sampling time of 7 days has the greater precision. Therefore, we select all the measured values of measuring point I to have analysis.

Figure 11 shows the mean squared error variability for the 30 epochs in the network training process using artificial wavelet neural network, and the minimum mean squared error is 0.426. Figure 12(a) illustrates the comparison between the fitted values of the monitoring model based on the artificial wavelet neural network and measured values, and Figure 12(b) shows the scatter plot of measured versus fitted with regression line. As shown in Figures 12(a) and 12(b), we can clearly see how model behaves in lower and higher regions. Table 5 shows the predicated values and measured values of piezometric tube level.

In addition, Table 6 shows the coefficients of artificial wavelet neural network model. Through the comparison between the coefficients of the statistical model and the artificial wavelet neural network, we can see the statistical model has a higher precision as a whole.

5.2.3. The Integrated Seepage Monitoring Model Based on Maximum Entropy Theory

Based on maximum entropy theory, the Lagrangian multipliers , , and are obtained through the solution of (28)–(32). corresponds to the statistical model and corresponds to the monitoring model based on the artificial wavelet neural network.

After applying the probability function, the predicated values in the predication period are obtained. Table 7 shows the predicated values of the integrated model based on maximum entropy theory.

To compare the prediction precision of the three models, we show the statistical metrics of the predicated values of the three models in Table 8, and we find that the integrated monitoring model based on maximum entropy theory fuses the advantages of the two models and thus obtains greater prediction accuracy.

6. Conclusions

This work studied the three monitoring models applied to the seepage status of earth-rock dams influenced by rainstorms. The main content of this paper was as follows:(1)Based on the research on the hysteresis effect of reservoir water levels and precipitation, the reservoir water level and precipitation were equivalently processed to obtain equivalent reservoir water levels and precipitation. To easily simulate the rapidly increased reservoir water lever, the exceeded water level factor was introduced. Combined with the components of temperature and timeliness, a statistical model was established.(2)By analyzing the characteristics of wavelet theory and neural networks, a monitoring model that fuses both wavelet theory and neural network was established. The monitoring model benefits from a strong nonlinearity mapping ability and thus can be used to analyze the detailed characteristics of seepage monitoring data.(3)Considering the deficiency of a single monitoring model, an integrated monitoring model based on maximum entropy theory was established and found to improve the predication accuracy of seepage status determination for earth-rock dams.(4)Through the application of the three monitoring models to a measuring point of an earth-rock dam influenced by rainstorms, we found that the three models have acceptable precision in fitting and predicating the seepage status, thereby providing technological support for the seepage monitoring of similar earth-rock dams.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 51579085, 41323001, 51139001, and 51279052), Project Funded by China State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering (Grant no. 20145028312), Jiangsu Province “333 High-Level Personnel Training Project” (Grant no. 2016-B1307101), Open Fund of Key Laboratory of Earth-Rock Dam Failure Mechanism and Safety Control Techniques, Ministry of Water Resources (Grant no. YK914022), and Huai’an Water Conservancy Academician Workstation.