Inverse Problems: Theory and Application to Science and Engineering 2015View this Special Issue
Research Article | Open Access
Identification of Hydraulic Conductivity in Aquifer for Coupled FEM and Adaptive Genetic Algorithm
The hydraulic conductivity of a natural rock mass was difficult to determine because of the complex structure and the significant influence of uncertain factors. In this paper, hydraulic conductivity was adopted to conduct an inversion analysis according to the measurement of head materials by combining the finite element method with the adaptive genetic algorithm. The results showed that the maximum relative error of the measuring and computation groundwater levels at the measuring points was 5.3%, and the average head error was 1.41%; the effective hydraulic conductivity of intensively weathered layer, moderately weathered layer, and fresh bedrock layer in riverbed formation tended to decline gradually; the effective permeability coefficient in direction was the minimum in the same aquifer. Therefore, the established hydraulic conductivity inversion analysis method was effective.
The seepage characteristics of a rock mass closely were related to the permeability of fissure structural planes. In process of inversion analysis, the seepage parameter was a wide value range mass or the calculation result exhibited discreteness. The main reasons were the complex composition of the structural plane, with the calculation result easily depending on the initial value of the parameter and falling into the local optimum. Thus, it was difficult to determine the value of seepage parameters with the inversion methods. In the 1960s, a field survey method was adopted to determine seepage parameters by snow. Although the method was simple and easy to implement, the required cost was high, and the measuring result had poor representativeness . In a later study, Kitamura, Chavent, Yeh, and others analysed the hydraulic conductivity inversely through an analytical or semianalytical solution [2–4]. Nevertheless, given that the assumed conditions always differed from actual conditions, the result was difficult to apply. After the 1980s, with a thorough study of the seepage model of rock mass, a back analysis method had been conducted to identify the seepage parameters of the discrete fracture medium [5, 6]. In the mid- and late 1990s, the inversion theory of seepage coefficients was developed rapidly [7–9]. In particular, with the rapid development of an optimal algorithm, the inversion method of seepage coefficients had significantly developed greatly and acquired numerous achievements [10–12]. Although the progress had been made on the seepage coefficient inversion analysis of rock mass, it was still difficult to identify that the seepage parameters of rock mass rigorously, accurately, rapidly, and economically had been proposed.
In this study, geological survey, water pressure test, and groundwater level were tested in the field, and other analyses were conducted in relation to hydraulic engineering. These field tests were conducted based on the combined influences of different actual seepage factors. Given that the observation period was short, problems such as the sudden change in drilling water data and high discreteness of the water pressure test data occurred periodically, thus making the determination of the seepage parameter of rock mass difficult. The finite element method combined with an optimal algorithm was adopted to analyze the seepage parameters of different areas in a rock mass inversely through the iteration method according to the measured hydraulic heads.
2. Basic Principle of the Inversion Method
2.1. Mathematical Model
The concept of hydraulic conductivity dates back to the middle nineteenth century. According to the current research results, three related concepts (effective hydraulic conductivity, equivalent conductivity, and interpreted conductivity) are defined according to the different conditions and assumptions . During the last 40 years, an enormous amount of work has been devoted to research and identify representative hydraulic conductivity. Actually, hydraulic conductivity () is arguably one of the most variable parameters in the Earth sciences. Laboratory measurements of span more than 12 orders of magnitude. Even in seemingly homogeneous aquifers, measured values of may range over some orders of magnitude. This realization spurred a large amount of work on the problem of finding the effect of heterogeneity on flow and transport.
Up to now, most theoretical results are related to effective conductivity, and their application to real world scenarios relies on ergodic assumptions. Fortunately, a number of results are available suggesting that conventional hydraulic test interpretations yield (interpreted) hydraulic conductivity values that can be closely linked to equivalent and/or effective hydraulic conductivities. Actually, the hydraulic conductivity is usually generalized as the ensemble averages of flux and head gradient in engineering projects, that is, effective hydraulic conductivities. Effective hydraulic conductivities will be analysed in this paper.
Equivalent continuum medium seepage model was widely applied in many actual engineering projects owing to the simplicity of model. The permeability of rock mass was considered as the effective hydraulic conductivities of a certain area. The three-dimensional seepage control equation was expressed as follows : where was the unit storage volume; was the hydraulic head; , , and were the hydraulic conductivity of rock mass in three directions. When conducting steady seepage analysis, the right side of the equation was equal to 0. In this paper, three-dimensional steady seepage analysis was used.
2.2. Objective Function
The effective hydraulic conductivity was analysed as the inversion parameters in this paper. Firstly, the water heads were calculated through three-dimensional finite element numerical analysis. Then, the error of measurement and computation head were calculated on the base of the measured groundwater levels of the natural seepage field. Finally, the effective hydraulic conductivity values were identified through the optimization of the adaptive genetic algorithm when the error was minimised. The objective function was given as follows:where , , and were the effective hydraulic conductivity to be used for the inversion analysis; was the computation head at measuring point ; was the measuring head at measuring point ; and was the total number of measuring points. The relative value of the head was adopted to make the objective function a dimensionless numerical function to avoid other problems during optimisation and to determine whether the algorithm was converged .
2.3. Boundary Conditions
According to the established field model and actual engineering practice, the boundary conditions were as follows:(1)The Class-I boundary referred to the measured groundwater level. The borehole was drilled on the cross section while the groundwater level was measured. The hydraulic heads of the leftmost and rightmost sides defined the Class-I head boundary, such that the hydraulic head was known.(2)The leftmost and rightmost sides of the model measuring the saturated surface referred to the Class-II (known discharge) boundary. The seepage discharge of these boundary surfaces was zero.
2.4. Adaptive Genetic Algorithms
In traditional genetic algorithm, crossover rate and mutation rate are the key factors influencing the behavior and performance of genetic algorithm. However, these two factors are predetermined in traditional genetic algorithm, and and should be determined through repeated tests specific to different optimization problems and tend to be trapped by local optimal solution. To solve the problem effectively, adaptive genetic algorithm in which and change automatically with the adaptability is presented in this paper; expressions are as follows:where refers to the maximum individual fitness in each generation; refers to average fitness of each generation; refers to the larger one between two crossover individuals; refers to the fitness of mutant individual. In general, the initial value of , , , and in the expressions above can be 0.9, 0.6, 0.1, and 0.001, respectively.
If the individual fitness in a population tends to be uniform or local optimum, the algorithm is able to increase the values of crossover rate and mutation rate automatically, and the optimal solution can be obtained avoiding the algorithm from falling into the local optimum solution. Meanwhile, individual with fitness higher than the average fitness corresponds to lower and , so that the next generation can be protected by the solution. Individual with fitness lower than average corresponds to higher crossover rate and mutation rate; hence the solution can be eliminated. Therefore, the adaptive algorithm constructed is able to guarantee the convergence of genetic algorithm while maintaining the population diversity.
2.5. Inversion Analysis Procedures
Since adaptive genetic algorithm was very effective to solve the multiobjective optimization problem, it was adopted to identify the effective hydraulic conductivity in this paper. The inversion analysis was divided into the following procedures:(1)Coding: adaptive genetic algorithm was adopted to transform , , and into coded strings.(2)Generation of initial group: initial groups were generated at random within the preset parameter scope.(3)Computation of objective function value: the initial solvable group was encoded into the finite element analysis program of seepage and the objective function value was then computed.(4)Evaluation of fitness: the result of the finite element analysis program for seepage was used to evaluate the individual fitness and to compute the average fitness of groups.(5)Convergence determination: the judging basis was considered when the difference of the maximum fitted value and the average fitted value of group was less than the preset error, and the maximum evolution generation was considered as the auxiliary judgment.(6)Genetic evolution operation: optional, cross, and mutation operations were performed, and a new generation of solvable groups was formed according to the operational approach of adaptive genetic algorithm.(7)When a new solvable group was obtained, Step (3) was repeated until the objective function value satisfies the error precision.(8)Optimised objective function value and optimised inversion parameters were obtained.
The procedure for the inversion of the hydraulic conductivity was as shown in Figure 1.
3. Engineering Situation
A hydropower station was situated in the main stream of the upper reach of Lancang River in the territory of Diqing Tibetan Autonomous Prefecture, Yunnan. Hydropower was mainly used to generate power. The normal groundwater level of the reservoir was 1,906.0 m, the downstream groundwater level was 1,814.86 m, and the normal storage capacity was 265,000,000 m3. The dam area was a cross valley with an asymmetric “V” shape. The natural grade of the left bank slope ranged from 40° to 50°. The rock mass of bank slope was not completely covered by the weathered layer but was partially covered by the quaternary waste slope gravelly soil with a thickness of approximately 7 m. The natural grade of the right bank below an altitude of 1,910 m ranges from 60° to 80°, whereas that above an altitude of 1,910 m ranged from 40° to 45°. The strata distributed in the dam area were mainly foliated sandy slate, argillaceous slate, and a spot of tuff. The groundwater under the dam area was divided into fissure water and pore phreatic water according to the mode of occurrence and occurrence medium. Fissure water occurred among rock fissures, whereas pore phreatic water occurred in the quaternary unconsolidated formation. Groundwater was mainly supplemented by atmospheric precipitation as well as melted mountain ice and snow. Such groundwater discharges into Lancang River. According to an adit exploration in the dam area, the groundwater in the rock mass of both banks was deeply buried underground. The hydraulic gradient of the left groundwater was approximately 0.266, whereas that of the right groundwater was approximately 0.568.
The riverbed covering layer was alluvial sand and gravel cobble, with a thickness of 19.8 m to 25 m. The rock roof altitude was approximately 1,782 m. The rock of the riverbed was mainly sandy and argillaceous slate. The average depth of the strongly weathered rock on the riverbed ranges from 5 m to 6 m. The vertical depth (below the rock plane) of the weakly weathered belt ranged from 6 m to 15 m, and the fresh rock was distributed at approximately 15 m to 46 m below the surface. The hydrogeological section map in the study area was shown in Figure 2.
4. Back Analyses of the Hydraulic Conductivity of Different Strata
4.1. Variation Range of Initial Hydraulic Conductivity
According to the geological plane of the dam area, the profiles of cross sections, and the geological material of the project, a natural seepage field model of the dam area was established. In the model, the permeability of the different strata significantly differed from top to bottom. The dam foundation area was divided into four layers: the alluvial, strongly weathered, weakly weathered, and fresh rock layers, from top to bottom. The dam shoulders on the left and right banks were considered to have varying permeability. The natural seepage field model of the dam area was shown in Figure 3.
The equivalent continuous medium seepage model was adopted on the analysis, and the effective hydraulic conductivity values were taken as the parameters for inversion. Considering the number of multiple layers, the degree of parameter inversion was reduced to lessen the computational load. The effective hydraulic conductivity of the riverbed alluvial layer could be determined based on engineering experience. The riverbed alluvial layer was formed through the alleviation of soil, gravel, sand, and so forth, and its hydraulic conductivity was 1 × 10−3 m/s according to relative engineering experience. The effective hydraulic conductivity values of the strongly and weakly weathered layers of the riverbed as well as of the left and right bank rocks were analysed in this paper. According to a similar project, the ranges of the initial value of the parameters were given in Table 1.
4.2. Finite Element Model
The dam axis was -axis, whereas the right bank direction was positive. -direction was perpendicular to -direction, and the water flow direction was positive. The vertical direction was -axis, whereas the positive direction was vertically upward (Figure 4).
The origin of the model was situated at -coordinate, which was 1,600 m below the intersection of - and -axes. The model measures 910 m × 820 m × 400 m. Based on the finite element subdivision, the model comprises 32,486 nodes and 34,320 hexahedron elements (Figure 4).
4.3. Computation Results and Analysis
The optimised inversion parameter combination corresponding to the minimum error of the measuring and computation heads was obtained through optimal computation. The measuring and computation heads at all measuring points were calculated for comparison (Table 2). In Table 2, the difference (186.5 m) between the maximum head (1,997 m) and the minimum head (1,810.5 m) was selected as the basis for the comparison between the measured and computed groundwater levels. The relative error refers to the ratio of the absolute error of the computation and measuring heads to the base.
Table 2 showed that the error of the measuring and computation heads of different measuring points was relatively small, and the average error of the measuring points in the dam area was only 1.41%. The computed groundwater level of the boreholes was higher than the measured value when the error was at the minimum value. The maximum error point was at ZK213, with a relative error of 5.3%. These results indicated that the inversion analysis model could reflect the permeability characteristic of the model area equivalently and macroscopically. Nevertheless, the position of the computation phreatic line remained relatively high.
For the validation of computation accuracy, visual graphic comparisons were established for the axis cross section, where intensive measuring points were located. The cross section on the upstream was 30 m away from the dam (Figures 5 and 6 as well as Tables 3 and 4).
Table 3 and Figure 5 showed that the error of the computation and measuring head values of the measuring points of the axis cross section was relatively small. The minimum error was 0.2 m at ZK227, with a relative error of 0.1%. The maximum head error was at ZK204, with an error value of 4.6 m and relative error of 2.49%. The average error of the eight measuring points on the axis cross section was only 0.90%, and the inversion analysis error was relatively small.
Table 4 and Figure 6 showed that the error of the computation and measuring head values at the measuring points on the cross section of the upstream 30 m away from the dam was relatively small. The minimum error value was 0.25 m at ZK225, with a relative error of 0.1%. The maximum error value was 4.86 m at ZK221, with a relative error of 2.63%. The average error of the seven measuring points on the cross section of the upstream 30 m away from the dam was only 1.11%, and the inversion analysis error was relatively small. The measuring and computation head error charts of the two cross sections indicated that the farther from the river channel the left or right bank was, the higher the error obtained.
The seepage field within the computation area could be obtained when the objective function was minimised. The seepage field of the axis cross section ( section), the cross section of the downstream 205 m away from the dam ( section), the 0 m to 70 m longitudinal section ( section) of the left and right banks, and the m longitudinal section ( section) were used to illustrate the condition of the seepage field (see Figures 7 and 8 for the location diagram and Figures 9–12 for the different positions of the seepage field).
The optimised effective hydraulic conductivity combinations of strata were obtained when the errors of the measuring and computation heads were minimized (Table 5).
The inversion analysis result showed the following. (1) The effective hydraulic conductivity values of the strongly weathered, weakly weathered, and fresh rock layers in the riverbed stratum gradually became small and increasingly weak. In the same stratum, -direction effective hydraulic conductivity value was minimised. (2) The effective hydraulic conductivity values of the right bank rock were less than that of the left bank rock, which indicates that the overall integrity of the former was better than that of the latter. This finding was associated with the fact that the hydraulic gradient of the right bank was higher than that of the left bank.
(1)To overcome the limitation of genetic algorithm, adaptive genetic algorithm was put forward. Combining with finite element method, effective hydraulic conductivity of different aquifers in natural dam foundation was identified. This algorithm overcame the limitations of the traditional optimisation algorithms, in which the inversion result significantly depended on the initial value of and and fell into the local optimum when the inversion analysis was performed.(2)According to the inversion analysis result of this study, the maximum relative error of the measuring and computation groundwater levels at the measuring points was 5.3%, and the average head error was 1.41%. The results showed that the established hydraulic conductivity inversion analysis method was effective. Meanwhile, the inversion results showed the effective hydraulic conductivity of intensively weathered layer, moderately weathered layer, and fresh bedrock layer in riverbed formation tended to decline gradually, which reflected that the permeability of aquifer became weaker and weaker, and in the same aquifer, the effective permeability coefficient in direction was the minimum.(3)The result showed that effective hydraulic conductivity of right bank was smaller than that of left bank, which indicated that the rock on left bank was more complete than those on left bank in whole, which coincided well with the condition that the average hydraulic gradient of right bank was larger than that of left bank.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This research was financially supported by the National Program on Key Basic Research Project of China (973 Program no. 2013CB036003), National Natural Science Foundation of China (NSFC Grant no. 50579092 sponsored), and the Scientific Research Foundation (2010JK599) by the Education Department of Shaanxi Province.
- D. T. Snow, “Rock fracture spacings, openings, and porosities,” Journal of the Soil Mechanics and Foundations Division, vol. 94, no. 1, pp. 73–92, 1968.
- S. Kitamura and S. Nakagiri, “Identifiability of spatially-varying and constant parameters in distributed systems of parabolic type,” The SIAM Journal on Control and Optimization, vol. 15, no. 5, pp. 785–802, 1977.
- G. Chavent, “Identification of distributed parameter system: about the output least square method, its implementation and identification,” in Identification and System Parameter Estimation, R. Isermann, Ed., vol. 1, pp. 85–97, Pergamon Press, New York, NY, USA, 1979.
- W. W.-G. Yeh and N.-Z. Sun, “An extended identifiability in aquifer parameter identification and optimal pumping test design,” Water Resources Research, vol. 20, no. 12, pp. 1837–1847, 1984.
- Y. W. Tsang and C. F. Tsang, “Channel model of flow through fractured media,” Water Resources Research, vol. 23, no. 3, pp. 467–479, 1987.
- Y. M. Zhu, L. J. Zhang, and Y. Wu, “Back analysis of coefficient tensor of permeability for fractured rock mass,” Chinese Journal of Rock Mechanics and Engineering, vol. 16, no. 1, pp. 461–470, 1997.
- M. Oda, “An equivalent model for coupled stress and fluid flow analysis in jointed rock masses,” Water Resources Research, vol. 22, no. 13, pp. 1845–1856, 1986.
- B. F. Zhu, “A new method for the back analysis of seepage problem,” Journal of Hydraulic Engineering, no. 9, pp. 42–46, 1994.
- M. Bai, F. Meng, D. Elsworth, Y. Abousleiman, and J.-C. Roegiers, “Numerical modelling of coupled flow and deformation in fractured rock specimens,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 23, no. 2, pp. 141–160, 1999.
- Z. F. Zhou, “Inverse analysis of parameters for groundwater movement in fissured double media,” Journal of Hydrodynamics, vol. 18, no. 6, pp. 742–747, 2003.
- X. He, S. J. Li, and Y. X. Liu, “Identification of hydraulic conductivity of rock mass in dam foundation based on genetic neural network,” Chinese Journal of Rock Mechanics and Engineering, vol. 23, no. 5, pp. 751–757, 2004.
- X.-S. Liu, C.-B. Zhou, and L.-J. Zhang, “Application of neural network based on simulated annealing Gauss-Newton algorithm to seepage back analysis,” Rock and Soil Mechanics, vol. 26, no. 3, pp. 404–406, 2005.
- X. Sanchez-Vila, A. Guadagnini, and J. Carrera, “Representative hydraulic conductivities in saturated groundwater flow,” Reviews of Geophysics, vol. 44, no. 3, 2006.
- J. R. Chai, S. Y. Li, and Y. Q. Wu, “Multi-level fracture network model for coupled seepage and stress fields in rock mass,” Communications in Numerical Methods in Engineering, vol. 20, no. 1, pp. 63–74, 2004.
- X. H. Deng, “Study on seepage coefficient Identification based on adaptive genetic algorithms,” Journal of Northwest A & F University (Natural Science Edition), vol. 37, no. 10, pp. 221–224, 2009.
Copyright © 2015 Xianghui Deng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.