In this study, novel methods such as wavelet–artificial neural network hybrid models and artificial neural network models were used to predict seepage from the Zonouz earthen dam. The dataset consisted of 972 piezometric data points. Statistical fitting methods such as root mean squared error, determination coefficient, scatter plots, and data distribution diagrams were used to evaluate the results. The findings indicated that the wavelet–artificial neural network hybrid model was more accurate than the artificial neural network model. Specifically, during training, the wavelet–artificial neural network hybrid model had determination coefficients and root mean squared errors of 0.820, 0.680, 743.39, and 792.52, while the artificial neural network model had 0.700, 0.600, 426.39, and 131.45. Similarly, during validation, the wavelet–artificial neural network hybrid model had determination coefficients and root mean squared errors of 0.700, 0.600, 426.39, and 131.45, while the artificial neural network model had 0.823, 0.680, 743.39, and 792.52. Therefore, the wavelet–artificial neural network hybrid model can be proposed as a precise method for predicting seepage in earthen dams and is more accurate than the artificial neural network model. This study highlights the importance of preventing dam failures and using advanced modeling techniques for better predictions and preventive measures.

  • This study focuses on predicting seepage from the Zonouz earthen dam using novel methods.

  • Statistical fitting methods, including root mean square error, determination coefficient, scatter plots, and data distribution diagrams, were employed.

  • This study emphasizes the significance of preventing dam failures and underscores the use of advanced modeling techniques.

  • In this study, the hybrid method of wavelet–artificial neural network has been used to predict the seepage phenomenon.

Building earthen dams is essential in countries with specific climatic conditions and a need for storing water. These dams act as surface water retention structures and control floods, providing the opportunity for increased use of river water. Prior to construction, studies are conducted, but accurately predicting the hydraulic behavior of the dam body and surrounding geological formations is not always possible. Therefore, seepage after construction is almost certain. The intensity of seepage is acceptable as long as dam safety is not endangered, but assessing the risks associated with seepage and infiltration due to their inherent complexity is crucial for safety. Many dams have seepage from geological formations at the construction site, downstream of the dam or the dam body itself, which can lead to adverse consequences such as economic issues, high hydraulic gradients leading to erosion or boiling, and increased pore pressure that reduces effective stress. Therefore, addressing seepage from the downstream and dam body is important during study phases, implementation operations, and post-construction. Precisely calculating the amount of seepage flow and investigating methods to control or reduce it is essential to prevent potential life and financial hazards. Research on dams and water seepage is crucial due to the growing need for water resources and the importance of optimal resource use. Given the water shortage in Iran and the effects of climate change, research in this area can greatly assist in planning water resources, agriculture, drinking water supply, and economic development. Ersayin (2006) employed a feed-forward neural network (FF-NN) with a sigmoid activation function and one hidden layer to model seepage in earthen dams. The findings indicated that the model's accuracy did not improve and, in fact, the network became unstable when increasing the number of hidden layers. Miao et al. (2012) employed a genetic algorithm (GA) and a novel neural network model to predict seepage from the ‘Diaoli’ earthen dam in China. The findings indicated that the Genetic Algorithm and Levenberg-Marquardt (GA-LM) model had the ability to accurately predict seepage from the dam and outperformed conventional neural networks in extrapolation and interpolation. Nourani et al. (2012) analyzed the piezometric head in the core of an earthen dam in Sattarkhan, Iran, using artificial neural networks (ANNs). The study showed a strong agreement between predicted and measured values with an R2 of 0.798 for the single ANN model. For the integrated ANN models of feed-forward backpropagation and radial basis function (RBF), the R2 was found to be 0.87 and 0.67, respectively. Poorkarimi et al. (2013) introduced a new approach for estimating seepage flow from the foundation and body of an earthen dam based on data mining techniques. Soil permeability of the dam was estimated through the use of SEEP/W software. The results demonstrated that ANNs performed well as a tool for detecting patterns in data and accurately predicting seepage flow from the Fileh-Khase earthen dam's foundation and body. Kamanbedast & Delvari (2013) conducted a study on the Maroun earthen dam located in the Maroun River in Northern Bahaman. The study aimed to investigate the behavior of the dam's soil and seepage phenomenon using ANSYS software and comparing results with Geo-Studio software. Researchers found that the value of pore water pressure is influenced by several factors, including material stiffness, moisture, soil permeability, and applied loads. Pore water pressure can also reduce the shear strength of the soil, which may cause soil particles to move and result in piping and erosion of finer particles from the dam's body if the rate of reduction in pore water pressure caused by leakage is greater than the resistance of the soil particles. Jafari (2014) employed ANNs as a data mining method to predict seepage from the earthen dam body of Sattarkhan. The dataset consisted of 1,684 piezometric data points, which were divided into training and validation sets at an 80:20 ratio. Using appropriate statistical parameters, the neural network was effectively trained and showed a high capability for predicting seepage. Parkam Shadbad et al. (2023) conducted a study on seepage from the body of the Sattarkhan dam in Ahar, using the finite difference method through Flac3D software. The results revealed that the ratio of permeability between the shell and core had a significant effect on the seepage behavior of the dam, and reducing it led to a more uniform and homogeneous dam behavior, which was observed by examining the trend of changes in pore water pressure in the shell and core while moving from upstream to downstream. With a constant ratio, the seepage pattern in the dam remained unchanged, and only the time of occurrence of generated pore water pressure changed based on the permeability values of the materials. Ranković et al. (2014) developed a FNN model using an improved resilient backpropagation algorithm to predict the piezometric water level in dams. The measured data were compared with the results obtained from both the FNN models and multiple linear regression models, which are commonly used for analyzing the structural behavior of dams. These models were developed and tested on experimental data collected over a period of 9 years. The study findings suggest that FNN models can serve as powerful and essential tools for evaluating dams. Zhang et al. (2020) employed backpropagation neural network (BPNN) theories in conjunction with a GA to model and predict water overflow at dams. They compared the prediction results of the BPNN-GA model with those of a statistically guided model that used supervised values, using a sample dam in China. The findings revealed that the improved model enhanced nonlinear visualization and generalization capabilities, accurately predicting future water overflow with an acceptable level of precision. Using the established criterion, the safety status of the dam during flood season was evaluated. Ishfaque et al. (2022) developed a deep learning model to predict the leakage rate of Tarbela Dam in Pakistan, which is the world's second largest earth-filled dam. The dataset consisted of hydroclimatological, geophysical, and engineering features related to water inflow to the dam from 2014 to 2020. The data were organized as a time series, and recurrent neural network (RNN) and long short-term memory (LSTM) algorithms were applied for time series neural networks. The RNN-LSTM model showed an average mean squared error (MSE) of 0.12 and a performance score of 0.9451, with the lowest errors and high accuracy, leading to the best predicted results for dam leakage. Beiranvand & Rajaee (2022) conducted a review in which they revealed that machine learning models (37.53%), neural networks (27.63%), and hybrid models (21.05%) are the most common techniques for predicting infiltration. Among single models, ANNs, support vector regression, random forests (RF), and FF-NN have been used more frequently than other models. In addition, 81.25% of the hybrid models utilized neural network models, and 31.25% of the hybrid models used GAs. They reviewed a total of 46 research articles from 2005 to 2022. Duong & Tran (2023) developed ANN models for primary prediction of soil sediment velocity and resistance based on fiber type, soil type, fiber length, and fiber content. They utilized an extensive dataset of input parameters to develop the ANN models and compared their performance to that of regression models. In addition, sensitivity analyses were performed using the Garson algorithm and connection weight method to rank the effective input parameters in the ANN formulas. The results demonstrated that the ANN models had excellent prediction performance, with high coefficient of determination (DC) values (R2 = 0.995 for sediment velocity and R2 = 0.998 for resistance), indicating their potential for predicting these issues. Emami et al. (2019) analyzed water infiltration in the earthen dam of Shahid Kazemi in Boukan, Iran, using RBF and FF-NN models from ANNs. The results revealed good agreement between predicted and observed values, and it was concluded that the RBF model trained by the Levenberg–Marquardt algorithm with four hidden layers had high potential for predicting infiltration. In addition, the statistical parameter values of R2 and root mean squared error (RMSE) were 0.81 and 33.12, respectively. Parsaie et al. (2021) utilized various soft computing models, including multilayer perceptron neural network (MLPNN), support vector machine, multiple adaptive regression splines (MARS), genetic programming, M5 algorithm, and group method of data handling, to predict piezometric pressure in the dam core and infiltration flow from the earthen dam body. The dataset collected from the absolute instrument mechanism during the past 94 months at the Shahid Kazemi Dam in Boukan were used for this purpose. The results showed that all models used for predicting piezometric pressure had an acceptable level of accuracy, with average error indices of R2 = 0.957 and RMSE = 0.806 in the training phase and R2 = 0.949 and RMSE = 0.932 in the testing phase. Moreover, the performance of all models, except for M5 and MARS algorithms, in predicting infiltration flow was almost similar. However, the MARS method exhibited the best performance, while the M5 algorithm had the weakest performance. Rehamnia et al. (2021) developed an efficient data-information artificial intelligence paradigm that combined an extended Kalman filter with a feed-forward ANN (EKF-ANN) to accurately estimate daily infiltration flow through earthen dams at the Fontaine Gazelles dam in Algeria. The study examined three powerful machine learning methods, including MLPNNs, RBF neural networks (RBF-NN), and RF, to assess the capability of EKF-ANN in predicting infiltration flow. The results showed that the EKF-ANN paradigm outperformed MLP, RF, and RBF-NN. In addition, the leverage approach was utilized to report the range of application of the presented models. Zhang et al. (2021) introduced a supervised infiltration flow monitoring model for concrete dams by considering delays in input factors and proposed an effective identification method for the delay process. The results indicated that the improved harmony search global optimization (HGWO) has stronger global optimization capability, and the HGWO-XGBoost model performs well for predicting infiltration flow in concrete dams. Moreover, compared to the traditional trial-and-error method, the proposed delay process calculation method in this study provides better effectiveness, which has significant value for monitoring and controlling infiltration flow in concrete dams. Markovic et al. (2021) developed a new method for predicting permeation pressure in incompatible cells using multiple sequential ANNs to address the issue of an incompatible cell with a small available dataset. The results showed that this approach achieved high accuracy in predicting values, which were more precise than those obtained using only one ANN for prediction. Chen et al. (2020) developed a framework for data mining and monitoring to control the safety of dam leakage. The study first investigated the effective factors in dam leakage, followed by training a kernel extreme learning machine (KELM) to predict dam leakage. Finally, an innovative global sensitivity analysis was provided to evaluate the relative importance of each input variable based on KELM. The simulation results from the case study demonstrated that KELM provides an acceptable prediction of leakage flow. Furthermore, water level fluctuations and precipitation were found to have a significant impact on the size of leakage. The sensitivity analysis offers a useful qualitative criterion for dam leakage, which has high value for monitoring and safety operations of dams. Ersayin (2006) studied seepage in a body of earthfill dam by ANNs. Seepage is investigated since both in the dam's body and under the foundation, it adversely affects dam's stability. This study specifically investigated seepage in the dam's body. Seepage in the dam's body follows a phreatic line. To understand the degree of seepage, it is necessary to measure the level of the phreatic line. This measurement is called piezometric measurement. Piezometric datasets that were collected from the Jeziorsko earthfill dam in Poland were used for training and testing the developed ANN model. Sazzad & Islam (2019) have studied about different seepage control methods in earthen dams, and the objective of that paper is to perform a comprehensive study of the incorporation of different seepage control measures to an earthen dam using a finite element method (FEM). SEEP/W, an FEM-based software, has been used for modeling and analysis of different seepage control measures. From the numerical analysis, it is observed that the use of rock toe combined with a horizontal filter is more beneficiary than they are used alone. The length of a horizontal blanket filter is a controlling factor to reduce the pore water pressure rather than its thickness. Attia et al. (2021) studied seepage through earthen dams with internal cut-off, Seepage through an earthen dam with internal cut-off is experimentally studied in the laboratory of Irrigation Engineering and Hydraulics Department, Faculty of Engineering, Alexandria University, Egypt, on a Hele-Shaw model. The effect of the relative penetration depth (d/D) and the relative distance of cut-off (X/Dw) on both the relative drop in the total head upstream (ΔH/H) and relative length of filter (Lf/H) is studied. The pervious layer thickness beneath the dam is also studied to investigate its effect on the filter length and total seepage discharge. The optimum position of the cut-off also is obtained from analysis of the experimental results. Experimental results are drawn in dimensionless form. Salem et al. (2019), in their research, studied seepage passing through the body of the earthen dam with/without an internal core experimentally and numerically, It is concluded that by increasing the side slopes (H:V), a noticeable drop of seepage line is observed with increasing downstream slope stability; by decreasing core permeability, increasing core thickness, base core thickness, and core penetration, the quantity of seepage and exit gradient decreased the upstream slope stability and increased the downstream slope stability. To explore the multifrequency mode of the stock, this paper proposes an adaptive wavelet transform model (AWTM). AWTM integrates the advantages of XGBoost algorithm, wavelet transform, LSTM, and adaptive layer in feature selection, time–frequency decomposition, data prediction, and dynamic weighting (Liu et al. 2020). Zhu et al. (2023) investigated the combined method of variational mode decomposition (VMD)-wavelet packet denoising method and the improved temporal convolutional network model to predict seepage in dams. Consequently, this paper combines the advantages of VMD decomposition and wavelet packet analysis. First, the contaminated signals are decomposed into several components by VMD, and the components with more noise are denoised by the wavelet packet threshold. Finally, the signal pieces are reconstructed to obtain the denoised signal.

Accurate estimation of seepage from earthen dams has always been an important challenge in the design of these huge structures. The use of instrumentation may provide an accurate estimate of this phenomenon to some extent, but the upcoming problems such as the failure of the instruments due to time, expenditure of money and manpower for reading, and so on, greatly reduce the efficiency of this method. It is common to use analytical solution methods proposed by different researchers to evaluate the amount of seepage from the body of earthen dams located on an impermeable bed due to their ease of use. However, in these analytical methods, assumptions are used to simplify the construction of equations, which may lead to big errors. Therefore, according to the stated content, the aim of this study is to provide a combined wavelet–ANN model to more accurately predict the amount of seepage from the body of earthen dams and eliminate the above problems, which is actually an accurate smart method to determine the amount of seepage in the body of nonhomogenous earthen dams. So, it has been tried to predict the phenomenon of seepage based on the data of the precise instruments of a specific dam and using the above method, and this method, the problems in the discussion of accurate estimation of seepage that existed in the previous methods have been solved. In this research, the capabilities of the ANN model as a black box model, as well as the wavelet soft calculations, which is a suitable tool for the preprocessing of unstable data, have been discussed; using the combined wavelet–ANN models, the amount of percolation in earthen dams has been predicted, and finally the results have been compared with the ANN model. In general, this method can be used in the future to determine the exact amount of seepage in porous media such as the body of earthen dams. One of the innovations of the present research is the use of the wavelet transform, unlike Fourier transformation (FT), which decomposes the signal or information series into sine and cosine functions and their harmonics; the signal is depicted on a group of functions called wavelets, which are derived from the mother wave. In this research, artificial intelligence and soft computing tools such as ANNs and wavelet transform are used to predict the piezometric pressure. Also, the combined neural wavelet method is able to provide a more accurate prediction of a phenomenon with uncertainty such as seepage in the porous medium, so the range and scope of using this method in engineering problems can be very wide. Hence, at the beginning of the work, the input data are extracted and then normalization and preprocessing of the data are done. Then, prediction is done in two cases, with wavelet transformation and without wavelet transformation. In the case of wavelet transformation, the data are divided into sub-data; these analyzed data are entered into the ANN model and prediction data are extracted. In the case without wavelet transformation, the data without decomposition is entered into the mentioned models and the output is obtained. In the next step, the results obtained from the different models of this research are validated. Regionally and globally, the phenomenon of leakage in all hydraulic structures that can be eroded is a natural phenomenon. How to predict leakage with advanced models is important, which is discussed in this study. The problems that are more visible in the past research are one-dimensional in terms of the research method, while in the present research, the combined method of ANNs and wavelet transform has been discussed, which is mentioned previously. The neural wavelet model is an advanced method in machine learning, which is a combination of wavelet transform and neural networks. One of the reasons the combined model was used for this research is its high efficiency in predicting the amount of leakage in earthen dams; previously where only one type of model was used, the accuracy level of the model was compared to the known models. It was low and not reliable, but in the combined model of this research, which is based on the neural wavelet and considered a hybrid model, it can practically be used to predict the amount of leakage in the body of all earthen dams.

The Zonouz dam, a reservoir dam with a height of 59 m from the riverbed and 60 m from the foundation, was constructed using 560,000 m3 of earthwork. The dam has a useful reservoir volume of 8.5 million m3 and is located on the Zonouz River, approximately 24 km the northeast of Marand County (Figure 1). It is an earthfill dam with a clayey silt impermeable core. The dam's crest elevation is set at 1,845 m above sea level. The annual river discharge is 11 million m3, and the dam stands 59 m above the riverbed and 60 m above the foundation. The embankment volume is about 0.55 million m3, and the total reservoir capacity is 8.5 million m3. The dam features a side free overflow spillway with a discharge capacity of 373 m3/s. The diversion system consists of a honeycomb-shaped gallery with internal dimensions of 8.1 × 5.1 m in the upper section and 1.3 × 3.2 m in the lower section, spanning a length of 263 m. An earthfill cofferdam with a height of 12 m is used for diversion, providing a capacity of 21 m3/s. The inlet channel can handle a flow rate of 4.1 m3/s, while the underflow channel has a capacity of 4.2 m3/s. The total cost of the project implementation, based on 1,380 exchange rates, amounts to 4.8 million USD (The Zonouz Reservoir Dam Phase I Study, 1377).
Figure 1

Project area location and access roads to the Zonouz dam (Sani et al. 2019).

Figure 1

Project area location and access roads to the Zonouz dam (Sani et al. 2019).

Close modal

Piezometers position in the dam

Three instrumented sections have been designated along the length of the dam to facilitate stability control based on Figure 2, and the chosen sections for study are 82 + 0, 110 + 0, and 150 + 0. Piezometers have been installed at different elevations in each of these sections. In cross-section 1 (82 + 0), piezometers 101–109 have been embedded. Piezometers 102 and 104 are positioned downstream of the dam, while piezometers 101 and 103 are located upstream. Piezometers 105 and 108 are situated upstream of the clayey silt core, while piezometers 107 and 109 are placed downstream of the core. In addition, piezometer 106 is installed on the dam axis. Moving to cross-section 2 (110 + 0), piezometers 201–213 have been installed. Piezometers 202 and 204 are positioned below the downstream injection gallery, while piezometers 201 and 203 are situated upstream of the gallery. Piezometers 205, 206, and 209 are located above the upstream clayey silt core, whereas piezometers 207, 208, and 211 are placed below the core. Piezometer 210 is situated on the dam axis, piezometer 212 is positioned upstream, and piezometer 213 is positioned downstream of the dam. In cross-section 3 (150 + 0), piezometers 301–305 have been embedded. Piezometer 301 is located above the upstream clayey silt core, piezometer 302 is installed on the dam axis, and piezometer 303 is placed below the downstream clayey silt core. Furthermore, piezometer 304 is positioned upstream, while piezometer 305 is placed downstream of the dam.
Figure 2

Various cross-sections of the Zonouz dam and the location of piezometers.

Figure 2

Various cross-sections of the Zonouz dam and the location of piezometers.

Close modal

The list and specifications of the installed piezometers in the Zonouz dam are provided in Table 1.

Table 1

List and specifications of the installed vibrating piezometers in the Zonouz dam (taken from the 2016 Behavior Monitoring Report, Regional Water Archive of East Azerbaijan Province)

RowTool IDInstalled sectionPositionInstallation levelInstallation dateDistance from axis (m)
101 Injection gallery 1,819.199 83.10.07 u.s. = +0.96 
102 Injection gallery 1,819.99 83.10.07 d.s. = −1 
103 Injection gallery 1,829.996 83.10.07 u.s. = +1.12 
104 Injection gallery 1,829.996 83.10.07 d.s. = −1.11 
105 Clay core 1,840 82.06.15 u.s. = +10 
106 Clay core 1,840 82.06.15 
107 Clay core 1,840 82.06.15 d.s. = −10 
108 Clay core 1,860 83.03.07 u.s. = +7 
109 Clay core 1,860 83.03.07 d.s. = −3 
10 201 Injection gallery 1,816 83.03.07 u.s. = +0.95 
11 202 Injection gallery 1,816 83.03.07 d.s. = −0.97 
12 203 Injection gallery 1,826 83.03.07 u.s. = +1.07 
13 204 Injection gallery 1,826 83.03.07 d.s. = −1.10 
14 205 Clay core 1,834 82.05.14 u.s. = +15 
15 206 Clay core 1,834 82.05.14 d.s. = −5 
16 207 Clay core 1,834 82.05.14 u.s. = +5 
17 208 Clay core 1,834 82.05.14 d.s. = −15 
18 209 Clay core 1,852 82.08.13 u.s. = +10 
19 210 Clay core 1,852 82.08.13 
20 211 Clay core 1,852 82.08.13 d.s. = −10 
21 212 Clay core 1,866 83.04.06 u.s. = +7 
22 213 Clay core 1,866 83.04.06 d.s. = −3 
23 301 Clay core 1,852 82.08.18 u.s. = +10 
24 302 Clay core 1,852 82.08.18 
25 303 Clay core 1,852 82.08.18 d.s. = −10 
26 304 Clay core 1,866 83.04.06 u.s. = +7 
27 305 Clay core 1,866 83.04.06 d.s. = −3 
RowTool IDInstalled sectionPositionInstallation levelInstallation dateDistance from axis (m)
101 Injection gallery 1,819.199 83.10.07 u.s. = +0.96 
102 Injection gallery 1,819.99 83.10.07 d.s. = −1 
103 Injection gallery 1,829.996 83.10.07 u.s. = +1.12 
104 Injection gallery 1,829.996 83.10.07 d.s. = −1.11 
105 Clay core 1,840 82.06.15 u.s. = +10 
106 Clay core 1,840 82.06.15 
107 Clay core 1,840 82.06.15 d.s. = −10 
108 Clay core 1,860 83.03.07 u.s. = +7 
109 Clay core 1,860 83.03.07 d.s. = −3 
10 201 Injection gallery 1,816 83.03.07 u.s. = +0.95 
11 202 Injection gallery 1,816 83.03.07 d.s. = −0.97 
12 203 Injection gallery 1,826 83.03.07 u.s. = +1.07 
13 204 Injection gallery 1,826 83.03.07 d.s. = −1.10 
14 205 Clay core 1,834 82.05.14 u.s. = +15 
15 206 Clay core 1,834 82.05.14 d.s. = −5 
16 207 Clay core 1,834 82.05.14 u.s. = +5 
17 208 Clay core 1,834 82.05.14 d.s. = −15 
18 209 Clay core 1,852 82.08.13 u.s. = +10 
19 210 Clay core 1,852 82.08.13 
20 211 Clay core 1,852 82.08.13 d.s. = −10 
21 212 Clay core 1,866 83.04.06 u.s. = +7 
22 213 Clay core 1,866 83.04.06 d.s. = −3 
23 301 Clay core 1,852 82.08.18 u.s. = +10 
24 302 Clay core 1,852 82.08.18 
25 303 Clay core 1,852 82.08.18 d.s. = −10 
26 304 Clay core 1,866 83.04.06 u.s. = +7 
27 305 Clay core 1,866 83.04.06 d.s. = −3 

u.s., upstream; d.s., downstream.

Wavelet transform and governing equations

The wavelet transform, known for its ability to analyze signals at multiple resolutions, has been widely used in signal processing and is recognized as a valuable tool. In the 1800s, Joseph Fourier, a mathematician from France, showed that any periodic function could be represented as a sum of infinite sinusoidal and exponential functions. The advantages of the neural wavelet model are as follows: (1) The ability to analyze signals in both time and frequency domains; this allows important signal features to be extracted in both domains. (2) The ability to process high-dimensional data due to the use of wavelet transform that reduces the dimensions of the data; this makes calculations simpler and more efficient. (3) High adaptability due to the layered structure of the neural network; this feature makes the model easily adaptable to a variety of problems. (4) High interpretability due to the clear and understandable structure of the neural network; this feature helps analysts to better understand the results. Approximately 150 years later, in 1965, a new algorithm called the fast FT emerged in the field of computer calculations. The FT decomposes a signal into a sum of infinite complex exponential functions, each with different frequencies. As per the definition, the continuous-time FT of a signal x(t) is obtained through Equation (1). Discrete wavelet transform (DWT) can be used to decompose the signal payment so that the signals that correspond to sets of high and low frequencies are separated. One of the advantages of the wavelet is that this transformation uses filter banks to reconstruct the level of multifrequency time resolution.
(1)
In this equation, the variable t symbolizes time while f stands for frequency. The equation signifies the FT of a signal. The time-limited version of the FT is quite similar to the original, with the main contrast being that in the time-limited version, and the signal is segmented into small enough parts to assume they are stationary. This segmentation is facilitated by employing a time window function, represented as w. The length of the window is chosen to be the minimum required length for the separated segments to be considered valid and stationary parts of the signal. Therefore, the time-limited FT of the signal, using the time window ψ, is defined by Equation (2).
(2)
The continuous wavelet transform was introduced as a substitute for the short-time FT in time–frequency analysis to address its accuracy constraints. In wavelet analysis, akin to the short-time FT, the signal of interest is convolved with a specific function (wavelet) tailored to it. Similarly, the wavelet transform is executed on distinct time intervals of the signal. Nonetheless, two primary distinctions exist between the wavelet transform and the short-time FT, as outlined by Chen and Liu in 1998. (a) In wavelet transform, there is no FT performed on the windowed signal, and hence individual peaks corresponding to sinusoids, or in other words, negative frequencies are not computed. (b) In wavelet transform, the width of the window varies alongside the changing frequency components, which is considered the most significant characteristic of wavelet transform. Consequently, the continuous wavelet transform is defined by Equation (3) (Addison et al. 2006).
(3)
where and s represent the translation and scale parameters, respectively. The concept of translation is exactly analogous to the time shift concept in the time–frequency short-time FT, determining the amount of window displacement and clearly incorporating temporal information into the transform. However, unlike the time–frequency short-time FT, direct frequency parameters are not present in wavelet transform. Instead, we have the scale parameter, which is inversely related to frequency. In other words, s = 1/f. In the recent equation, ψ is the window function, commonly referred to as the mother wavelet. The term ‘wavelet’ signifies a small wave.
The continuous wavelet transform can be expressed as the inner product of the signal and a basis function, given by Equation (4) (Addison et al. 2006).
(4)
where
(5)
Based on the provided definition in Equation (5), which is expressed as an inner product, the wavelet transform can be understood as a measurement of similarity between the signal and the basic functions (wavelets). The wavelet transform can be interpreted as an assessment of the resemblance between the signal and the basic functions (wavelets). In any given space, there are multiple sets of basic functions, and among them, the orthogonal basis functions are particularly significant. This is because orthogonal basis functions possess favorable properties and greatly facilitate the computation of transform coefficients. Consequently, by leveraging the orthogonality property of the basic functions, the transform coefficients are calculated according to Equation (6) (Addison et al. 2006).
(6)
By having these coefficients, it is possible to reconstruct the function according to Equation (7) (Addison et al. 2006).
(7)

Two famous mother wavelets with explicit mathematical relationships are presented (Addison et al. 2006).

  • (A) Mexican Hat wavelet, which is the second derivative of the Gaussian function (Equation (8)):
    (8)
  • (B) Morlet wavelet, which consists of a Gaussian term and a sinusoidal term (Equation (9)):
    (9)

Artificial neural network

ANNs are one of the oldest methods used in data mining. In recent years, the application of ANNs has greatly extended in various engineering fields. A review of the technical literature reveals that ANNs have been successfully utilized for predicting the load-bearing capacity of piles, modeling soil behavior, site characterization, retaining wall design, settlement analysis, slope stability, tunnel and cavern design, groundwater flow, soil permeability, soil compaction, soil swelling, and soil classification.

Different technical and engineering disciplines utilize a range of neural network types such as Perceptron, Hamming, and Hopfield networks. The Perceptron network falls under the category of feed-forward networks, and various training algorithms have been developed for these networks. In this research, the backpropagation learning algorithm has been employed to train the specified networks. The performance index, known as the MSE, is the operational controller in this algorithm and is minimized during supervised learning using the following training data pairs (Equation (10)):
(10)
where pi represents the input vector to the network and ti represents the desired output (response) of the network for input pi. After applying input p(k) (kth pattern) to the network, the error signal at the output neuron j in the output layer at time k or iteration k is obtained from Equation (11).
(11)
The momentary squared error value for the j-the neuron in the output layer is denoted as . Consequently, the instantaneous sum of squared errors is obtained by summing for all neurons in the output layer. It should be noted that only the neurons in the output layer are considered observable neurons. Therefore, the network's performance is expressed by the performance index in Equation (12).
(12)
where SL is the number of neurons in the output layer. In practice, the performance index in Equation (13) is often used.
(13)
The parameters of the network, which consist of weight matrices and biases, are adjusted to minimize the sum of squared instantaneous errors. In other words, these parameters, including the weight matrices and biases, are tuned to achieve the minimum MSE across a set of Q′ patterns (number of training data pairs), as described by Equation (14).
(14)
Finally, the weight matrices and bias vectors of the network are adjusted using Equations (15) and (16):
(15)
(16)

As the learning rate α decreases, the number of required trial-and-error iterations to reach the desired performance index increases, but the problem-solving accuracy improves.

The wavelet–ANN and ANN have similarities and differences in their methodologies and applications. Both wavelet–ANNs and ANNs are widely used for prediction in various fields, such as hydraulic phenomena, meteorology, climatology, and hydrology.

The similarities are as follows:

  • Both wavelet–ANNs and ANNs are used for prediction in diverse areas.

  • They involve the processing of data to make predictions or classifications.

  • Both networks can be trained using data to improve their predictive capabilities.

The differences are as follows:

  • wavelet–ANNs specifically incorporate wavelet functions as activation functions, allowing for time-dependent spectral analysis and time–frequency decomposition of time series data.

  • ANNs have a broader scope and can encompass various architectures beyond wavelet-based approaches.

To use data mining methods such as neural networks and fuzzy models, having an appropriate dataset is one of the basic requirements, and for each dataset, three features are necessary.

  • 1. Reliability: This means being real and accurate.

  • 2. Sufficient data: Having enough data based on the dimensions and complexities of the problem.

  • 3. Comprehensive coverage: In an audit environment, reliability of data means that data are applicable for audit purpose and are sufficiently complete and accurate.

Based on the relationship and proportion of changes in the water level behind the dam and changes in the pressure of the cavity water, faulty piezometers are identified and the resulting data are collected. Accordingly, the collected dataset includes 972 piezometric data, which contain information such as the coordinates of the location of the piezometer, the date of piezometer readings, and the water level behind the dam. To develop the model, the date of piezometer readings was combined with the water level behind the dam, and the input ratio of the time interval was considered as in Equation (17). Piezometers are geotechnical sensors used to measure pore water pressure in the ground, soil, earth/rock fills, foundations, and concrete structures. They are considered one of the most important geotechnical monitoring instruments and are used in large civil engineering structures such as tall buildings, dams, and underground tunnels.
(17)
Recent investigations have revealed that the manner in which a dataset is partitioned into training and testing subsets can significantly influence the results. Since computational models, like all empirical models, cannot extend beyond the training data range, it is crucial to divide the data in a manner that, while random, is similar in terms of parameters such as mean, standard deviation, minimum, and maximum. This is essential for normalizing the test data based on the characteristics of the training data during model simulation. To achieve this, several random combinations were chosen to separate the training and testing sets, resulting in two sets that were statistically compatible with each other. In machine learning, the dataset is divided into training, validation, and test sets. The test set is not utilized during the training process. The validation dataset, which is a subset of data reserved for training the model, is used to assess the model performance while fine-tuning model hyperparameters. The evaluation becomes more biased as the model's performance on the validation dataset is integrated into the model configuration. Data are typically divided into training and testing sets in two ways (80:20 and 70:30), and the choice of method depends on the number of data and inputs. In this study, 80% of the data set (778 data) was used for training and 20% of the data set (194 data) was used for testing the proposed model. To ensure the accuracy and speed of the model, normalization is performed to equalize the data value for the model. Data normalization is usually performed in the range of 1 and 0. In this study, Equation (18) is used to normalize all input data before applying it to the model.
(18)

The text explains that Table 2 presents the normalized value, actual value, and maximum and minimum values of the data for input and output variables used in the ANN model. In addition, the table includes statistical parameters related to the two sets of training and testing data. The purpose of this presentation is to provide a clear understanding of the data used in the ANN model. Data normalization is a technique used to normalize the range of features in a dataset so that no feature is dominated by the others. This technique is used to improve the performance of the ANN model. The input and output variables are preprocessed using normalization or standardization techniques before training the ANN model. These techniques are used to improve the stability and performance of the ANN model.

Table 2

Statistical parameters related to the training and testing data are presented

ParameterStatisticTraining dataTest data
X (m) Max 15 10 
Mean 0.420 0.526 
Min −15 −10 
Std Dev 6.865 6.964 
Y (m) Max 1,866 10 
Mean 1,838.050 0.526 
Min 1,816 −10 
Std Dev 15.110 6.937 
Z (m) Max 110 150 
Mean 98.339 147.113 
Min 82 110 
Std Dev 13.812 10.377 
Data rate Max 
Mean −0.014 0.080 
Min −6.85 −6.85 
Std Dev 3.587 3.599 
Elevation Max 1,887.6 1,887.6 
Mean 1,873.958 1,874.592 
Min 1,855.95 1,855.95 
Std Dev 8.011 8.020 
P (kPa) Max 409.947 241.117 
Mean 172.545 109.278 
Min 9.535 10.565 
Std Dev 92.804 71.448 
ParameterStatisticTraining dataTest data
X (m) Max 15 10 
Mean 0.420 0.526 
Min −15 −10 
Std Dev 6.865 6.964 
Y (m) Max 1,866 10 
Mean 1,838.050 0.526 
Min 1,816 −10 
Std Dev 15.110 6.937 
Z (m) Max 110 150 
Mean 98.339 147.113 
Min 82 110 
Std Dev 13.812 10.377 
Data rate Max 
Mean −0.014 0.080 
Min −6.85 −6.85 
Std Dev 3.587 3.599 
Elevation Max 1,887.6 1,887.6 
Mean 1,873.958 1,874.592 
Min 1,855.95 1,855.95 
Std Dev 8.011 8.020 
P (kPa) Max 409.947 241.117 
Mean 172.545 109.278 
Min 9.535 10.565 
Std Dev 92.804 71.448 

The aim of this study is to predict the piezometric pressure in the Zonouz dam using ANN models and a combined wavelet–ANN model. The modeling system involves examining different models with different memories. In this modeling system, the best input pattern is determined first, followed by comparing the best pattern of each model with the other models to determine the most accurate model. Two important graphical criteria used in this study are.

  • (a) a scaled linear plot between observed and computed data for verification and validation stages and

  • (b) a scatter plot between observed and computed data for verification and validation stages.

In addition, the performance of the models can be examined using five numerical criteria, including the DC, MSE, RMSE, mean absolute error (MAE), and standard error of prediction (SEP). Among these criteria, the DC has a wide range of applications in evaluating predictive models.

This criterion is very sensitive to extreme values and can lead to high values for ideal predictions, while it is not very sensitive to differences between observed and computed data in general. Due to this limitation, this criterion may evaluate the model as a good predictive model in many cases, even if the model is not suitable. Therefore, this criterion cannot be used alone as a suitable criterion for evaluating hydrological models, and when using it, other statistics and criteria such as absolute error or mean error should also be used simultaneously (Legates & McCabe 1999). In this regard, the RMSE criterion, which is obtained by taking the square root of a positive number, can be used as a suitable criterion for measuring the level of prediction error. The range of RMSE varies from zero for accurate predictions to large positive values for cases where the difference between computed and observed data is high and subject to oscillation. In addition, high values of DC (with an ideal state of one) and small values of MAE, SEP, and RMSE indicate high model performance. Therefore, in this study, to evaluate the accuracy of the models, more than two statistics, including the DC and the root MSE, will be used with Equations (19) and (20).
(19)
(20)

In the above relationships, the value of N represents the number of observed, are computed data, are predicted values, and are mean of observed data.

In complex problems where the relationship between variables is unknown, ANNs are a powerful tool for predicting phenomena. To clarify the effect of different parameters and achieve a method for estimating piezometric pressure, an ANN with five input parameters, including the coordinates of the location of piezometers (X, Y, Z), the ratio of time interval and the level of water behind the dam, has been created using different input patterns. The model is shown as the following functions:
ANNs are a broader class of neural networks that encompass various architectures, including wavelet–ANNs. ANNs, including deep learning neural networks, face challenges such as local minimum problems, instability, and limited precision. The methodology flowchart of the combined model is shown in Figure 3.
Figure 3

Methodology scheme of the wavelet–ANNs and ANNs.

Figure 3

Methodology scheme of the wavelet–ANNs and ANNs.

Close modal

The most common method for finding the optimal combination of weights in FF-NN is the backpropagation algorithm, which works based on reducing the first derivative. Since the backpropagation algorithm uses the minimization of the first derivative to determine the connection weights of the network, if the initial weight guesses are not suitable, the optimization process may get stuck in a local minimum. Therefore, in this study, the training phase has been repeated several times with different specifications for the layer or hidden layers to find the optimal network. The results of the seven input patterns are shown in Table 3.

Table 3

The results of the coefficient and the errors of the ANN model

 
 

The results indicate that the input pattern four with three neurons in the input layer, eight neurons in the hidden layer, and only one output neuron representing piezometric pressure had the best performance among other input patterns (shown in bold in Table 3). The results also indicate that the proposed ANN model has been properly trained and can be used to predict seepage.

In ANN training, emphasis is placed on the scale with the greatest impact. Subsequently, to enhance model accuracy, preprocessed data utilizing DWT is utilized as input for the neural network. This approach ensures that key features are represented across both the detailed and coarse scales of the time series, enabling the neural network to assign distinct weights to each decomposed subseries.

Therefore, high weights are assigned to signal-valued levels, which in turn increase the accuracy of the model. In signal decomposition using wavelets, all hidden information in the signal is extracted, and basic data preprocessing methods such as normalization of raw data have no effect on the decomposition results. In this study, the wavelet packet decomposition was used to decompose the signal, and four common mother wavelets, including db2, Harr, Coif1, and sym3, were used for discrete wavelet decomposition (Nourani et al. 2009).
(21)
In this relationship, L is the analysis level, N is the number of time series data, and int[.] is the correct part of the function. In the time domain of this research with a third-order decomposition, the result of wavelet decomposition with the mother wavelet db2 is shown in Figure 4.
Figure 4

Under the partial and approximate subseries obtained from wavelet decomposition with the mother wavelet db2, the time series of the Zonouz dam reservoir (third order) is aligned.

Figure 4

Under the partial and approximate subseries obtained from wavelet decomposition with the mother wavelet db2, the time series of the Zonouz dam reservoir (third order) is aligned.

Close modal

A summary of the modeling results using this method can be seen in Table 4. As seen in the table, the best structure in the input pattern 1 is obtained with six neurons in the hidden layer and 130 iterations using the mother wavelet coif1 (shown in bold in Table 4). The comparative results of two prediction methods are given in Table 5.

Table 4

The results of the DC and errors obtained from the wavelet–ANN model are summarized

 
 
Table 5

The comparative results of the DC and errors obtained from prediction with two models are presented

ModelType of modelDC
RMSE
TrainingValidationTrainingValidation
ANN Nonlinear 0.68 0.60 52.792 45.131 
Wavelet–ANN Combination 0.82 0.70 39.743 39.426 
ModelType of modelDC
RMSE
TrainingValidationTrainingValidation
ANN Nonlinear 0.68 0.60 52.792 45.131 
Wavelet–ANN Combination 0.82 0.70 39.743 39.426 

The findings suggest that the wavelet–ANN prediction model outperformed the alternative model. The superior performance of the wavelet–ANN model can be attributed to the spatial and temporal variations in delays across stations. In this integrated model, apart from preprocessing and segmenting data into subseries, a proportional weight is assigned to each subseries before modeling. The enhanced results compared to the ANN model highlight the beneficial impact of wavelet transform on modeling and scaling down processes. The clearer proof of the higher accuracy of the combined wavelet–ANN model compared to the ANN model can also be observed from the graphs in Figure 5.
Figure 5

Optimal evaluation of the results of the two models based on success evaluation criteria: (a) RMSE and (b) DC.

Figure 5

Optimal evaluation of the results of the two models based on success evaluation criteria: (a) RMSE and (b) DC.

Close modal

As can be seen in the results of MSE and explanation coefficient for the two methods of ANN and wavelet–ANN, the result of test and training in both statistical communities shows the appropriate and high efficiency of the wavelet method. From Figure 5, it can be seen that the combined wavelet–ANN method has an extremely high accuracy compared to the ANN model. The observational versus computational figures and the scatter of observational and computational data for the optimal state of both the wavelet–ANN and ANN models are shown in Figures (6)(8).

As can be seen in Figures 6 and 7, the dispersion of observational and computational data for the neural wavelet method is shown in comparison with the observational trend, which is the result of the high efficiency of the neural wavelet method. This comparison based on seepage pressure has been observed and calculated in both models.
Figure 6

Scatter of observational and computational data around the median in wavelet–ANNs model.

Figure 6

Scatter of observational and computational data around the median in wavelet–ANNs model.

Close modal
Figure 8 shows the trend of pressure dispersion in two observational and computational modes for wavelet–ANN and ANN models in a monthly period, which can be easily concluded that the efficiency of this method is very high. It is higher and more accurate than other methods, and even with a difference of 5%, it has been placed with the observed trend.
Figure 7

Scatter of observational and computational data around the semi-meter in the most optimal state in ANN.

Figure 7

Scatter of observational and computational data around the semi-meter in the most optimal state in ANN.

Close modal
Figure 8

The scatter plot of the comparative observational-computational values over time in the optimal state of the models.

Figure 8

The scatter plot of the comparative observational-computational values over time in the optimal state of the models.

Close modal

In this research, an attempt was made to simulate piezometric pressure in the body of the Zonouz dam using ANN and wavelet–ANN models. The simulation was performed based on the dataset obtained from the piezometers installed in the body of the dam. In this paper, the first step was to ensure the accuracy of the piezometers and the recorded data, and then to use the data. The accuracy of the piezometers was investigated to make sure that the data collected were reliable. The model carried out in this research was based on a combined wavelet–ANN model, in which the amount of seepage in the body of earthen dams can be accurately predicted to a great extent, and this method is considered a fundamental method.

The use of these methods for model training and evaluation based on appropriate statistical parameters led to significant results:

  • 1. Randomly dividing the data for model training and testing showed that the 80–20%, which is recommended in many references, is the optimal state.

  • 2. Due to the multiscale structure of the wavelets, the used hybrid model has the ability to extract local and partial features, and in general, it provides a better understanding of the data structure in problems where there is uncertainty; in this research, it provided a more accurate prediction.

  • 3. Using appropriate statistical parameters to evaluate the presented models can significantly help in identifying their accuracy and capability. Each of the statistical parameters used evaluates a specific criterion of the model. Based on such comprehensive and complete validation, it can be stated that the above models have good usability in predicting piezometric pressure in earthen dams.

  • 4. Using the presented models and replacing them with the available precise tools, it is possible to witness a reduction in the costs associated with instrumentation in the body of earthen dams, a reduction in the significant costs of human resources for continuous reading of piezometers, and predicting seepage for the future.

  • 5. Using the presented models and assigning more weight to the scale that has the most effect during the training of the ANN model, preprocessed data by the DWT has been considered as the input to the neural network. In this way, all the main features are visible in the fine and coarse time series scales, and the neural network assigns a specific weight to each of the decomposition subseries, leading to an improvement in the results of the wavelet–ANN model compared to ANN model.

  • 6. According to the modeling done, it was found that the proposed model is well trained and has a high capability in predicting the seepage phenomenon, and this method can be introduced as a basic method for predicting the seepage phenomenon in porous media.

  • 7. Due to the difference in geometric and physical specifications in earthen dams, in this study, such specifications were not considered as input to only investigate the exact amount of seepage through the combined intelligent model in a dam. Therefore, the results of this study can be used for all earthen dams.

For future research, this method can be introduced as a basic method for predicting the seepage phenomenon in porous media.

Data cannot be made publicly available; readers should contact the corresponding author for details.

The authors declare there is no conflict.

Addison
P. S.
,
Morvidone
M.
,
Watson
J. N.
&
Clifton
D.
2006
Wavelet transform reassignment and the use of low-oscillation complex wavelets
.
Mechanical Systems and Signal Processing
20
(
6
),
1429
1443
.
Attia
M.
,
Abdel Razek
M.
&
Salam
A. A.
2021
Seepage through earthen dams with internal cut off
.
Geotechnical and Geological Engineering
39
(
8
),
5767
5774
.
Chen
W. Y. C.
&
Liu
Z.
1998
Parameter augmentation for basic hypergeometric series, I
. In:
Bruce E. Sagan & Richard P. Stanley (Eds)
.
Mathematical Essays in Honor of Gian-Carlo Rota
.
Progress in Mathematics Vol. 161, Springer, Berlin
, pp.
111
129
.
Emami
S.
,
Choopan
Y.
&
Parsa
J.
2019
Dam seepage prediction using RBF and GFF models of artificial neural network; Case study: Boukan Shahid Kazemi's Dam
.
Journal of Rehabilitation in Civil Engineering
7
(
3
),
15
32
.
Ersayin
D.
2006
Studying Seepage in a Body of Earth-Fill Dam by Artificial Neural Networks. Doctoral Dissertation, İzmir Institute of Technology
.
Ishfaque
M.
,
Dai
Q.
,
Haq
N. U.
,
Jadoon
K.
,
Shahzad
S. M.
&
Janjuhah
H. T.
2022
Use of recurrent neural network with long short-term memory for seepage prediction at Tarbela dam, KP, Pakistan
.
Energies
15
(
9
),
3123
.
Jafari
M. H.
2014
Prediction of Seepage From Earthen Dams Using Data Mining Techniques
.
Master's Thesis
,
Islamic Azad University of Semnan
,
Iran
.
Kamanbedast
A.
&
Delvari
A.
2013
Analysis of earthen dam: Seepage and stability using Ansys and Geo-Studio software
.
World Applied Sciences Journal
17
(
9
),
1087
1094
.
Liu
X.
,
Liu
H.
,
Guo
Q.
&
Zhang
C.
2020
Adaptive wavelet transform model for time series data prediction
.
Soft Computing
24
,
5877
5884
.
Markovic
M.
,
Brankovic
J. M.
,
Stosovic
M. A.
,
Zivkovic
S.
&
Brankovic
B.
2021
A new method for pore pressure prediction on malfunctioning cells using artificial neural networks
.
Water Resources Management
35
,
979
992
.
Miao
X. Y.
,
Chu
J. K.
,
Qiao
J.
&
Zhang
L. H.
2012
Predicting seepage of earthen dams using neural network and genetic algorithm
.
Advanced Materials Research
403
,
3081
3085
.
Nourani
V.
,
Komasi
M.
&
Mano
A.
2009
A multivariate ANN-wavelet approach for rainfall–runoff modeling
.
Water Resources Management
23
(
1
),
2877
2894
.
Nourani
V.
,
Sharghi
E.
&
Aminfar
M. H.
2012
Integrated ANN model for earthfill dams seepage analysis: Sattarkhan Dam in Iran
.
Artificial Intelligence Research
1
(
2
),
1
22
.
Parkam Shadbad
R.
,
Fard Moradinia
S.
&
Alizadeh Mjdi
A.
2023
Evaluation of permeability ratio of embankment dam body materials with clay core using three-dimensional seepage behavior model (Case study: Setarehkhan Ahraz Dam)
.
Journal of Civil and Environmental Engineering
.
doi: 10.22034/ceej.2023.54880.2214
.
Parsaie
A.
,
Haghiabi
A. H.
,
Latif
S. D.
&
Tripathi
R. P.
2021
Predictive modelling of piezometric head and seepage discharge in earthen dam using soft computational models
.
Environmental Science and Pollution Research
28
(
43
),
60842
60856
.
Poorkarimi
S.
,
Maghsoodian
S.
,
MolaAbasi
H.
&
Kordnaeij
A.
2013
Seepage evaluation of an earthen dam using group method of data handling (GMDH) type neural network: A case study
.
Scientific Research and Essays Journal
8
(
3
),
120
127
.
Ranković
V.
,
Novaković
A.
,
Grujović
N.
,
Divac
D.
&
Milivojević
N.
2014
Predicting piezometric water level in dams via artificial neural networks
.
Neural Computing and Applications
24
,
1115
1121
.
Salem
M. N.
,
Eldeeb
H. M.
&
Nofal
S. A.
2019
Analysis of seepage through earthen dams with internal core
.
International Journal of Engineering Research
8
,
768
777
.
Sani
H.
,
Roushangar
K.
&
Ghasempour
R.
2019
Comparative study of the performance of finite element method and evolutionary model in seepage discharge predicting from the body of an earthen dam
.
Journal of Civil Infrastructure Researches
4
(
2
),
1
15
.
doi:10.22091/cer.2018.2349.1094
.
Sazzad
M. M.
&
Islam
M. M.
2019
A comprehensive study of different types of seepage control measures for earthen dam using FEM
.
Journal of Civil and Construction Engineering
5
(
1
),
24
37
.
Zhang
X.
,
Chen
X.
&
Li
J.
2020
Improving dam seepage prediction using back-propagation neural network and genetic algorithm
.
Mathematical Problems in Engineering
2020
,
1
8
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).