A Simpliﬁed Method for Leakage Estimation of Clay Core Dams with Different Groundwater Levels

: Clay core dams are widely applied in reservoir construction, regulating water resource and provide electric power. Leakage is a common problem in reservoir construction, and the leakage amount, which not only affects the economic beneﬁts of the project, but also relates to the safety of the dam body, is difﬁcult to estimate. According to Darcy’s law and stable seepage theory, an analytical method can be proposed to calculate the leakage of the clay core dam to gain the seepage ﬂux in a short time. By making some reasonable assumptions, we propose formulae for seepage calculation in different conditions of the position of the groundwater levels, below or above the reservoir bottom. Both sets of formulae contain two parts of leakage calculation, i.e., leakage from the reservoir bottom and leakage from the dam body. By using the proposed analytical method, the leakage of clay core dams can be estimated considering the inﬂuence of the groundwater level. To prove the rationality of the analytical method, a simple numerical model can be established using Geo-studio 2020 to calculate the seepage ﬂux of the clay core dam, where relative errors between numerical solutions and analytical solutions are less than 10%. To verify the feasibility in engineering applications, the proposed method was applied to calculate the seepage of a clay core dam in Sichuan, China, which was also calculated using numerical methods by establishing a three-dimensional model. The results show the rationality of the analytical method, which can strike a balance between precision and efﬁciency.


Introduction
A clay core dam is a special kind of earth rock dam which has a long history and is widely used in hydraulic engineering to regulate water resources and provide electric power [1,2]. The dam body is made of local materials, with low transportation costs, a simple structure, and a long service life. Different from a homogeneous dam, a clay core dam is equipped with a special impervious core in the dam body, which uses sand and stone with high permeability as the dam shell, and the impervious core adopts cohesive soil with good impervious performance. In the normal operation of the clay core dam, leakage problems endanger the hydraulic engineering in many aspects [3][4][5][6][7]. Some serious leakage problems even result in dam breakage due to a variety of mechanisms such as fracturing and collapse [8]. The Mostiště embankment dam, the first compacted rockfill dam with a relatively thin inclined impervious core in the Czech Republic, was troubled by seepage problems 36 years after its completion. The defects could not be determined due to the lack of a drainage and seepage observation system. The leakage problem worsened, and at the end of 2004, the reservoir level was significantly lowered, threatening the life and property safety of downstream residents [9]. Internal erosion in embankment dams, which leads to ageing effects and further causes dam failures, accidents, and deterioration of dams, is still not a completely understood phenomenon. Mattsson et al. [10] developed a transition characteristics of seepage in different regions. Numerical solutions were used to compare with the analytical solutions, showing that analytical and numerical solutions are consistent in case of minor soil capillarity [46]. Rezk [47] developed an analytical solution to investigate the effect of core with different permeability levels on seepage flow and phreatic line position. Fakhari and Ghanbari [48] studied many numerical models of clay core dams to develop new approximate equations for seepage calculation and gained relatively high precision. Additionally, seepage through cores may be also obtained using velocity hodographs [45]. These studies have helped in obtaining seepage flux through the dam body, but ignore the leakage of water from the reservoir bottom and periphery, which are important parts of the total leakage. Additionally, the influence of different groundwater levels on the leakage amount is not considered.
In order to make the engineering design process quicker and more targeted, analytical methods should be applied to calculate the leakage amount of the clay core dam. By simplifying the calculation process, engineers can gain results in a short time to assess the engineering safety and economy. In this research, by making some reasonable assumptions, an analytical method is proposed for the seepage calculation of the clay core dam (including the dam body and the reservoir bottom). According to Darcy's law and stable seepage theory, two sets of formulae are proposed for seepage calculation of the clay core dam in different conditions of the position of the groundwater levels. A simple numerical model was built to calculate the seepage flux of the clay core dam using Geo-studio 2020 software to verify the rationality of the formulae. At last, the simplified method was used for leakage calculation of a clay core dam in Sichuan, China, which demonstrates its feasibility in engineering applications. Figure 1 is the cross-section of a simplified clay core dam. Some basic assumptions can be made for the calculation of seepage:

Basic Assumptions
(1) The position of the groundwater level remains constant. The focus of the following analytical method is to estimate the leakage amount in a steady state, and an unchanged groundwater level is necessary.
(2) There is no centralized leakage channel in the reservoir area or dam body. It is a common engineering problem that there exists a centralized leakage channel, through which a large amount of leakage water may flow away. The focus of this problem is how to find the leakage channel and plug it effectively, rather than calculate the seepage flow of the centralized leakage channel. Therefore, the centralized leakage channel is not considered in this paper.
(3) The rock-soil body is homogeneous and isotropic. In practical engineering, the geological conditions and project design are more complex and detailed. The rock-soil body is often inhomogeneous and anisotropic, and it is difficult to homogenize the dam filling. For example, some dam site areas are layered with inclination of layers, and probably result in anisotropy. It is unrealistic and complex to obtain the analytical solutions considering the geological conditions in detail. Therefore, such an assumption is needed in practical engineering to gain hydraulic conductivities of different zones for leakage estimation.
(4) Rock, soil, and water are incompressible, and the pore size and porosity of soil mass remain unchanged.
(5) The lag of the soil moisture characteristic curve and its correlation with volume deformation should be ignored. Studies on the seepage of earths dam usually focus on the gravitational water, which is a very complex problem, and it needs to be simplified in the analytical calculation. The structure of rock-soil body, the shape and size of particles, the porosity, and the saturation of water all affect the permeability. Therefore, in order to facilitate the calculation, assumptions (3)-(5) need to be made. Many engineering examples show that the calculation results can meet the accuracy requirements of the project if they are assumed to be within the allowable range. In the seepage field, if only a local area does not meet the assumptions of (3)-(5), the results are still reliable [43].

Simplified Calculation Method
According to different groundwater levels, i.e., below or above the reservoir bottom, the simplified calculation can be mainly divided into two situations.

Groundwater Level below the Reservoir Bottom
(1) The groundwater level is deep. When the groundwater level is deep, the rock-soil body below the reservoir is unsaturated. Since the hydraulic conductivity of the rockfill area is much greater than that of the core, the role of the rockfill area in dam anti-seepage is usually ignored. The leakage amount of the clay core dam is mainly associated with the hydraulic conductivity of the core and the reservoir bottom.
For leakage from the reservoir bottom, there are two main periods after the reservoir starts filling. Firstly, the leakage water from the reservoir gradually wets the rock-soil body to the groundwater level. With the increase in the range of the wet rock-soil body, the leakage amount will decrease gradually. This period can be called the unsteady free leakage stage (Figure 2a). Secondly, the leakage amount will remain constant. There is no direct hydraulic association between the leakage water and the groundwater; as a result, the leakage amount will not change. This period can be called the steady free leakage stage ( Figure 2b). Studies on the seepage of earths dam usually focus on the gravitational water, which is a very complex problem, and it needs to be simplified in the analytical calculation. The structure of rock-soil body, the shape and size of particles, the porosity, and the saturation of water all affect the permeability. Therefore, in order to facilitate the calculation, assumptions (3)-(5) need to be made. Many engineering examples show that the calculation results can meet the accuracy requirements of the project if they are assumed to be within the allowable range. In the seepage field, if only a local area does not meet the assumptions of (3)-(5), the results are still reliable [43].

Simplified Calculation Method
According to different groundwater levels, i.e., below or above the reservoir bottom, the simplified calculation can be mainly divided into two situations.

Groundwater Level below the Reservoir Bottom
(1) The groundwater level is deep. When the groundwater level is deep, the rock-soil body below the reservoir is unsaturated. Since the hydraulic conductivity of the rockfill area is much greater than that of the core, the role of the rockfill area in dam anti-seepage is usually ignored. The leakage amount of the clay core dam is mainly associated with the hydraulic conductivity of the core and the reservoir bottom.
For leakage from the reservoir bottom, there are two main periods after the reservoir starts filling. Firstly, the leakage water from the reservoir gradually wets the rock-soil body to the groundwater level. With the increase in the range of the wet rock-soil body, the leakage amount will decrease gradually. This period can be called the unsteady free leakage stage (Figure 2a). Secondly, the leakage amount will remain constant. There is no direct hydraulic association between the leakage water and the groundwater; as a result, the leakage amount will not change. This period can be called the steady free leakage stage (Figure 2b).
In the steady leakage stage, the leakage amount relates to the size of the reservoir and the hydraulic conductivity of the rock-soil body below the reservoir. In Figure 2, the slopes on both sides of the reservoir are 1:m 1 and 1:m 2 . The seepage flux can be calculated as follows [43]: where q s represents the seepage flux per unit length of reservoir, m 2 /s; k 0 represents the hydraulic conductivity of the rock-soil body close to the reservoir, m/s; B represents the width of the pool level in the cross-section of the reservoir, m; H is height of the pool level, m; and A is the coefficient depending on the value of B/H, m 1 , and m 2 . The value of A can be determined empirically by Figure 3 [43]. In the steady leakage stage, the leakage amount relates to the size of the reservoir and the hydraulic conductivity of the rock-soil body below the reservoir. In Figure 2, the slopes on both sides of the reservoir are 1:m1 and 1:m2. The seepage flux can be calculated as follows [43]: where s q represents the seepage flux per unit length of reservoir, m 2 /s;  The total leakage from the reservoir bottom s Q can be written as follows: where l is the length upstream of the dam body, m. If the shape of the reservoir valley section substantially varies, data from multiple sections should be taken for averaging to The total leakage from the reservoir bottom Q s can be written as follows: Water 2022, 14, 1961 6 of 23 where l is the length upstream of the dam body, m. If the shape of the reservoir valley section substantially varies, data from multiple sections should be taken for averaging to obtain the total leakage from the reservoir bottom. For leakage from the core, the core can be divided into three sections, i.e., the upstream section (blue section), the middle section (green section), and the downstream section (gray section), as is shown in Figure 4. β 1 is the angle of the upstream slope of the core, and β 2 is the angle of the downstream slope of the core. h 0 is the height of the seepage face at the intersection of the phreatic line and downstream slope.
The total leakage from the reservoir bottom s Q can be written as follows: where l is the length upstream of the dam body, m. If the shape of the reservoir valley section substantially varies, data from multiple sections should be taken for averaging to obtain the total leakage from the reservoir bottom. For leakage from the core, the core can be divided into three sections, i.e., the upstream section (blue section), the middle section (green section), and the downstream section (gray section), as is shown in Figure 4. 1  is the angle of the upstream slope of the core, and 2  is the angle of the downstream slope of the core.  In the upstream section, the triangular region can be transformed into an equivalent rectangular region for the convenience of seepage flux calculation. As is shown in Figure  5, a rectangular coordinate system can be established. In the upstream section, the triangular region can be transformed into an equivalent rectangular region for the convenience of seepage flux calculation. As is shown in Figure 5, a rectangular coordinate system can be established. The basic assumption of Dupuit's equation is that the pressure at each point on section conforms to the hydrostatic pressure distribution-that is, the piezometric h along the vertical are equal, and the isohead line is a vertical line. This situation is approximate when the seepage flow is slowly changing. Many experimental studies s that the seepage movement in the middle section of the earth dam is a slowly chan seepage movement, so the application of Dupuit's equation is also approximate [43] cording to the analysis of Dupuit's equation and relevant test data, Dupuit's equation sufficient accuracy.
The width of the equivalent rectangular region depends on the seepage resistan the upstream section, which is associated with the values of 1  and H . According t solution of fluid mechanics, a semi-empirical formula can be presented: The basic assumption of Dupuit's equation is that the pressure at each point on the section conforms to the hydrostatic pressure distribution-that is, the piezometric heads along the vertical are equal, and the isohead line is a vertical line. This situation is only approximate when the seepage flow is slowly changing. Many experimental studies show that the seepage movement in the middle section of the earth dam is a slowly changing seepage movement, so the application of Dupuit's equation is also approximate [43]. According to the analysis of Dupuit's equation and relevant test data, Dupuit's equation has sufficient accuracy. The width of the equivalent rectangular region depends on the seepage resistance of the upstream section, which is associated with the values of β 1 and H. According to the solution of fluid mechanics, a semi-empirical formula can be presented: where L 1 denotes the width of the equivalent rectangular region, m. According to Dupuit s postulates, the seepage flux per unit through the upstream section and the middle section q c1 can be written as [43]: where H is the pool level, m; δ 2 is the bottom width of the middle section and the downstream section, m; and k c is the hydraulic conductivity of the core, m/s. The discharge per unit through the downstream section q c2 can be written as: According to the continuity equation, the seepage flux through different sections should remain constant. Hence, the value of q c1 is equal to q c2 . Combining Equations (4) and (5), the following equation can be obtained as: The solution of the above equation is: When h 0 = (L 1 + δ 2 ) tan β 2 + (L 1 + δ 2 ) 2 tan 2 β 2 − H 2 , the value of h 0 is greater than the pool level H, which is unreasonable. As a result, the value of h 0 should be Therefore, we substitute the value of h 0 into Equation (5), and the seepage flux per unit through the core q c can be obtained as: The longitudinal section of the core can be simplified as a trapezoid, as is shown in Figure 6. The gray part of the longitudinal section is below the pool level. B 1 and B 2 represent the length, m, of the top and bottom of the gray part, respectively. α 1 and α 2 represent the included angle between the left and right sides of the core and the horizontal plane, respectively. The core can be divided into three parts for convenience of calculation.
Water 2022, 14,1961 8 of 26 h is greater than the pool level H , which is unreasonable. As a result, the value of 0 h should be Therefore, we substitute the value of 0 h into Equation (5), and the seepage flux per unit through the core c q can be obtained as: The longitudinal section of the core can be simplified as a trapezoid, as is shown in Figure 6. The gray part of the longitudinal section is below the pool level.     The seepage flux through the middle rectangular part Q c2 can be written as: To calculate the discharge through the left part Q c1 , a rectangular coordinate system, as shown in Figure 6, can be established.
Considering the seepage flux over the length dx, the seepage flux dQ c1 can be written as: Integrating Equation (10) with x from 0 to H cot α 1 , the seepage flux Q c1 can be obtained as: where Similarly, the seepage flux through the right part Q c3 can be written as: The total leakage from the core Q c can be written as: The total leakage Q of the clay core dam is combined by two parts, i.e., leakage from the reservoir bottom and leakage from the core, and it can be written as follows: (2) The groundwater level is shallow. In this situation, leakage from the reservoir bottom continuously replenishes groundwater, and the area below the reservoir bottom is saturated. There exists hydraulic association between the leakage water and the groundwater; as a result, the leakage amount relates not only to the permeability of the rock-soil body below the reservoir, but also to the position of the groundwater level. As is shown in Figure 7, leakage water flows into the rock-soil body vertically, and then flows towards the groundwater.
Water 2022, 14,1961 9 of 26 The total leakage from the core c Q can be written as: The total leakage Q of the clay core dam is combined by two parts, i.e., leakage from the reservoir bottom and leakage from the core, and it can be written as follows: (2) The groundwater level is shallow. In this situation, leakage from the reservoir bottom continuously replenishes groundwater, and the area below the reservoir bottom is saturated. There exists hydraulic association between the leakage water and the groundwater; as a result, the leakage amount relates not only to the permeability of the rock-soil body below the reservoir, but also to the position of the groundwater level. As is shown in Figure 7, leakage water flows into the rock-soil body vertically, and then flows towards the groundwater. For leakage from the reservoir bottom, a rectangular system, as shown in Figure 8, can be established. Suppose that the rock-soil body below the reservoir has n layers. l i and k i are the thickness and the hydraulic conductivity of each layer, respectively. ∆ϕ i is the hydraulic head drop of each layer. For leakage from the reservoir bottom, a rectangular system, as shown in Figure 8, can be established. Suppose that the rock-soil body below the reservoir has n layers.  Based on Darcy's law, the discharge through the reservoir si q can be written as [31]: Based on Darcy's law, the discharge through the reservoir q si can be written as [31]: where B 0 represents the reservoir bottom length, m, in the cross-section. Therefore, ∆ϕ i can be written as: It should be noted that k i = 0. Then, the total drop can be written as [35]: Replacing the above formation system with a continuous variation of hydraulic permeability, Equation (17) can be written as: By integrating Equation (18), the head drop ϕ(y 0 ) at depth y 0 can be written as: As a result, if the hydraulic head drop is known, the seepage flux q s (x) can be easily obtained as follows: The leakage from the bottom of the reservoir can be written as: For leakage from the reservoir periphery, rectangular systems, as shown in Figure 9, can be established. According to Equation (20), the seepage flux through the periphery of the reservoir on the left side 1 s q (shown in Figure 9a) can be obtained as follows: Similarly, the seepage flux through the periphery of the reservoir on the right side 2 s q (shown in Figure 9b) can be obtained as follows: For leakage from the core, because the downstream water level is still below the reservoir bottom, the calculation of the seepage flux ci Q can refer to Equations (9), (11), and (12).
The total leakage from the core c Q can be written as: The total leakage Q of the clay core dam can be expressed as follows: (3) The above two situations can be distinguished according to the following standards.
Assuming that the water table is ' L far from the earth surface, the seepage flux  According to Equation (20), the seepage flux through the periphery of the reservoir on the left side q s1 (shown in Figure 9a) can be obtained as follows: Similarly, the seepage flux through the periphery of the reservoir on the right side q s2 (shown in Figure 9b) can be obtained as follows: The total leakage from the reservoir bottom Q s can be written as: For leakage from the core, because the downstream water level is still below the reservoir bottom, the calculation of the seepage flux Q ci can refer to Equations (9), (11), and (12).
The total leakage from the core Q c can be written as: The total leakage Q of the clay core dam can be expressed as follows: (3) The above two situations can be distinguished according to the following standards.
Assuming that the water table is L far from the earth surface, the seepage flux Q s from the reservoir bottom can be written as follows, according to Equation (2): Additionally, according to Equations (21)- (24), the seepage flux Q s from the reservoir bottom can be written as follows: Combining Equations (27) and (28), the value of L can be obtained.

Groundwater Level above the Reservoir Bottom
In this case, the rock-soil body below the groundwater level is saturated and the downstream water level is higher than the reservoir bottom. There exists hydraulic association between the reservoir leakage and the groundwater. The calculation of the seepage flux mainly considers the hydraulic conductivity of the rock-soil body and the core, as well as the position of the downstream water level.
For leakage from the reservoir bottom, Figure 10 is the cross-section of the reservoir for simplified calculation.

Groundwater Level above the Reservoir Bottom
In this case, the rock-soil body below the groundwater level is saturated and the downstream water level is higher than the reservoir bottom. There exists hydraulic association between the reservoir leakage and the groundwater. The calculation of the seepage flux mainly considers the hydraulic conductivity of the rock-soil body and the core, as well as the position of the downstream water level.
For leakage from the reservoir bottom, Figure 10 is the cross-section of the reservoir for simplified calculation.
where 0 L represents the depth of the saturated zone, m, and 1 H represents the distance between the groundwater level and the reservoir bottom, m. Similar to situation in which the groundwater level is close to reservoir bottom, as is shown in Figure 11, the seepage flux through both sides of the periphery of the reservoir  Referring to Equation (21), the seepage flux through the reservoir bottom q s can be obtained as: where L 0 represents the depth of the saturated zone, m, and H 1 represents the distance between the groundwater level and the reservoir bottom, m. Similar to situation in which the groundwater level is close to reservoir bottom, as is shown in Figure 11, the seepage flux through both sides of the periphery of the reservoir q s1 (Figure 11a), q s2 (Figure 11b) can be obtained as:  Therefore, the total leakage from the reservoir Q s can be written as: Q s = (q s + q s1 + q s2 )l (32) For leakage from the core, like the previous segmentation method, the simplified diagram of the calculation is shown in Figure 12. H 2 is the downstream water level. Therefore, the total leakage from the reservoir s Q can be written as: 12 () s s s s Q q q q l = + + (32) For leakage from the core, like the previous segmentation method, the simplified diagram of the calculation is shown in Figure 12.  To calculate the seepage flux through the core, the core can be divided into two parts to calculate the seepage flow separately, i.e., the part of core below the downstream water level and the part of core above the downstream water level, as is shown in Figure 13.  To calculate the seepage flux through the core, the core can be divided into two parts to calculate the seepage flow separately, i.e., the part of core below the downstream water level and the part of core above the downstream water level, as is shown in Figure 13. Therefore, the total leakage from the reservoir s Q can be written as: 1 2 ( ) s s s s Q q q q l = + + For leakage from the core, like the previous segmentation method, the simpl agram of the calculation is shown in Figure 12. To calculate the seepage flux through the core, the core can be divided into tw to calculate the seepage flow separately, i.e., the part of core below the downstream level and the part of core above the downstream water level, as is shown in Figur   Figure 13. Calculation of seepage flux through the core. In the part of core below the downstream water level, the seepage flux can be calculated as seepage under pressure. To calculate the seepage flux through this part, a rectangular coordinate system, as shown in Figure 13, can be established. Based on Darcy's law, the discharge dq cd through this part can be written as: By integrating Equation (33) with y from 0 to H 2 , the seepage flux q cd through this part can be obtained as follows: In the part of core above the downstream water level, the calculation of seepage flux is similar to Equation (8), which can be written as: In Equation (35), L 2 and δ 4 respectively represent: The longitudinal section of the core can be simplified as a trapezoid, as is shown in Figure 14. The gray part of the longitudinal section is between the upstream water level and the downstream water level, and the blue part is below the downstream level. The core can be divided into six parts for convenience of calculation.
By integrating Equation (33) with y from 0 to 2 H , the seepage flux cd q through this part can be obtained as follows: In the part of core above the downstream water level, the calculation of seepage flux is similar to Equation (8), which can be written as: The longitudinal section of the core can be simplified as a trapezoid, as is shown in Figure 14. The gray part of the longitudinal section is between the upstream water level and the downstream water level, and the blue part is below the downstream level. The core can be divided into six parts for convenience of calculation. The seepage flux through the middle rectangular blue part Q cd2 can be written as: To calculate the seepage flux through the left blue part Q cd1 , a rectangular coordinate system, as shown in Figure 14, can be established.
Considering the seepage flux over the length dx, the seepage flux dQ cd1 can be written as: Integrating Equation (39) with x from 0 to H 2 cot α, the seepage flux Q cd1 can be obtained as: where Similarly, the seepage flux through the right part Q cd3 can be written as: To calculate the seepage flux through the gray part, the solutions to Q cu1 , Q cu2 , and Q cu3 , referring to Equations (9), (11), and (12), can be obtained as follows: where E = (L 2 + δ 4 ) tan β 2 and H 3 = H − H 2 . The total leakage from the core Q c can be written as: The total leakage Q of the clay core dam can be expressed as follows:

Case Verification
To prove the rationality of the proposed formulae, a simple numerical model of a clay core dam was established using Geo-studio 2020 SEEP/W [49], as is shown in Figure 15. Figure 15a is the cross-section of the core. The width of the bottom of the core δ 1 is 40 m, the angle of the upstream slope of the core β 1 is 77.5 • , and the angle of the downstream slope of the core β 2 is 78.1 • . The hydraulic conductivity of the core k c is 1 × 10 −9 m/s. Figure 15b is the cross-section of the reservoir. The width of the bottom of the reservoir B 0 is 59 m, the slope on the left side of the reservoir is 1:1.3, and the slope on the right side of the reservoir is 1:1.2. The rock-soil body can be divided into four layers under the reservoir, and the thicknesses and hydraulic conductivities are listed in Table 1.  Table 1.    For leakage from the reservoir bottom, suppose the pool level is 50 m; when the groundwater level is below the reservoir bottom, the value of L can be solved as 11.36 m. As a result, when L = 5 m, Equation (24) should be used, and when L = 40 m, Equation (2) should be used. When the groundwater level is above the reservoir bottom, Equation (32) should be used.
For leakage from the core, suppose the pool level is 50 m; when the downstream water level is below the reservoir bottom, Equation (13) should be used; when the downstream water level is above the reservoir bottom, Equation (45) should be used. By applying the proposed analytical method and the Geo-studio 2020 SEEP/W analysis, the seepage flux can be calculated.
The calculation results are shown in Tables 2 and 3. Figures 16 and 17 are the streamline directions analyzed by Geo-studio 2020, and different areas in different colors represent different equipotential regions.     Based on the above results, the analytical results are consistent with the software modeling results in calculating the discharge of both the reservoir bottom and the core, where the relative errors are less than 10%. By this simple case, it is verified that the proposed method is reliable.

Engineering Application
The Jinsha River is the name for the upper reaches of the Yangtze River. It originates from the southern foot of the Bayankala Mountains in Yushu County, Qinghai Province. Its largest tributary is the Yalong River. The total length of the river is about 1570 km, and the drainage area is about 136,000 km 2 , accounting for about 27.1% of the catchment area of the Jinsha River (above Yibin). The average annual discharge of the estuary is 1930 m 3 /s, and the annual runoff is 60.9 billion m 3 . The Lianghekou hydropower station is located on the main stream of the Yalong River in Yajiang Country, Sichuan Province, China (shown in Figure 18).
Water 2022, 14,1961  Based on the above results, the analytical results are consistent with modeling results in calculating the discharge of both the reservoir bottom where the relative errors are less than 10%. By this simple case, it is verifie posed method is reliable.

Engineering Application
The Jinsha River is the name for the upper reaches of the Yangtze Rive from the southern foot of the Bayankala Mountains in Yushu County, Qin Its largest tributary is the Yalong River. The total length of the river is abou the drainage area is about 136,000 km 2 , accounting for about 27.1% of the c of the Jinsha River (above Yibin). The average annual discharge of the estua and the annual runoff is 60.9 billion m 3 . The Lianghekou hydropower statio the main stream of the Yalong River in Yajiang Country, Sichuan Province, in Figure 18). The power station adopts a clay core dam to store water and generate e a maximum dam height of 295 meters. The normal pool level of the reservo a.s.l., and the storage is more than 10 billion m 3 . The dam crest elevation is 2 and the excavation elevation of the core bottom is 2580.00 m a.s.l. The to core is 6 m, the top elevation is 2874.00 m a.s.l., and the upstream and dow ratio of the core are both 1:0.2. The bottom width of the core is 153 m a elevation is 2583.00 m a.s.l. The cross-sections of the dam body and the The power station adopts a clay core dam to store water and generate electricity, with a maximum dam height of 295 meters. The normal pool level of the reservoir is 2865.00 m a.s.l., and the storage is more than 10 billion m 3 . The dam crest elevation is 2875.00 m a.s.l., and the excavation elevation of the core bottom is 2580.00 m a.s.l. The top width of the core is 6 m, the top elevation is 2874.00 m a.s.l., and the upstream and downstream slope ratio of the core are both 1:0.2. The bottom width of the core is 153 m and the bottom elevation is 2583.00 m a.s.l. The cross-sections of the dam body and the reservoir are shown in Figures 19 and 20, and the region in magenta is the range of the core. In Figures 19 and 20, different areas of the dam body have been marked. Circled Roman numerals indicate different rock grades and "f+ numerals" indicate different faults. bedrock fissure water, and the lower part is fissure-confined water. Pore phreatic water is related to the distribution of overburden and is locally exposed. Shallow weathering unloading fissure water is exposed in the weathering unloading zone of rock mass. The phreatic water is generally a one-sided aquifer with a relatively stable groundwater level. Fissure-confined water is mainly exposed in siltstone and fine sandstone, which is a banded water-bearing body. The river valley in the dam area is deep and the bank slope is steep. Generally, the groundwater supplies the river water. Based on engineering geological and hydrogeological data, the hydraulic conductivities of different zones of the dam body are shown in Table 4, and the hydraulic conductivities of different layers of the rock-soil body under the reservoir bottom are shown in Table 5.    bedrock fissure water, and the lower part is fissure-confined water. Pore phreatic water is related to the distribution of overburden and is locally exposed. Shallow weathering unloading fissure water is exposed in the weathering unloading zone of rock mass. The phreatic water is generally a one-sided aquifer with a relatively stable groundwater level. Fissure-confined water is mainly exposed in siltstone and fine sandstone, which is a banded water-bearing body. The river valley in the dam area is deep and the bank slope is steep. Generally, the groundwater supplies the river water. Based on engineering geological and hydrogeological data, the hydraulic conductivities of different zones of the dam body are shown in Table 4, and the hydraulic conductivities of different layers of the rock-soil body under the reservoir bottom are shown in Table 5.    Based on the occurrence conditions of groundwater and the characteristics of waterbearing medium, lithology, and its combination, the types of groundwater are divided into loose accumulation layer pore phreatic water and bedrock fissure water. The bedrock fissure water is shallow weathering unloading bedrock fissure water and fissure-confined water. The spatial form of groundwater in the dam site area is pore phreatic water and bedrock fissure water, and the lower part is fissure-confined water. Pore phreatic water is related to the distribution of overburden and is locally exposed. Shallow weathering unloading fissure water is exposed in the weathering unloading zone of rock mass. The phreatic water is generally a one-sided aquifer with a relatively stable groundwater level. Fissure-confined water is mainly exposed in siltstone and fine sandstone, which is a banded water-bearing body. The river valley in the dam area is deep and the bank slope is steep. Generally, the groundwater supplies the river water. Based on engineering geological and hydrogeological data, the hydraulic conductivities of different zones of the dam body are shown in Table 4, and the hydraulic conductivities of different layers of the rock-soil body under the reservoir bottom are shown in Table 5.

Analytical Solution
The normal pool level is 2865.00 m a.s.l., and the downstream water level is 2610 m a.s.l., indicating that Equations (32) and (45) should be used to calculate the seepage flux of the reservoir and the dam body, respectively. The values of different parameters in this condition are listed in Table 6.

Numerical Solution
In addition to the proposed analytical method, an FEM numerical method was built to model the seepage field of the clay core dam. Figure 21 is the three-dimensional model of the reservoir area. The numerical calculation adopts the self-developed three-dimensional stable seepage finite element calculation and analysis program, whose reliability has been proven in the seepage calculation of dozens of water conservancy projects at home and abroad. The zone below the water level of the upstream reservoir and the location of the groundwater level have a Dirichlet boundary. The bottom and surroundings of the model are considered as an aquifuge, with a Neumann boundary. The top of the model is a free surface boundary. Other calculation parameters are the same as in Section 3.1.
Water 2022, 14, 1961 22 of 26 has been proven in the seepage calculation of dozens of water conservancy projects at home and abroad. The zone below the water level of the upstream reservoir and the location of the groundwater level have a Dirichlet boundary. The bottom and surroundings of the model are considered as an aquifuge, with a Neumann boundary. The top of the model is a free surface boundary. Other calculation parameters are the same as in Section 3.1. Figure 21. A 3D FEM model of the clay core dam.
The potential distribution of the cross-section of the dam body is shown in Figure 22, and the potential distribution of the cross-section of the reservoir is shown in Figure 23.  The potential distribution of the cross-section of the dam body is shown in Figure 22, and the potential distribution of the cross-section of the reservoir is shown in Figure 23. tion of the groundwater level have a Dirichlet boundary. The bottom and surroundings of the model are considered as an aquifuge, with a Neumann boundary. The top of the model is a free surface boundary. Other calculation parameters are the same as in Section 3.1. The potential distribution of the cross-section of the dam body is shown in Figure 22, and the potential distribution of the cross-section of the reservoir is shown in Figure 23.  The results gained by the analytical method are consistent with those obtained by the numerical method, where the relative errors are less than 10%, indicating that the simplified method is feasible in engineering applications. The leakage amount from the core is far less than that from the reservoir bottom with the condition of a normal pool level, which reflects the good anti-seepage performance of the core. To reduce the leakage amount in case of a determined hydrogeological condition, more targeted engineering measures should be taken to control the seepage flux of the reservoir bottom, such as laying face slab for retaining water, geomembranes, etc. The calculation results of both numerical and analytical methods for the leakage of the clay core dam are presented in Table 7. The results gained by the analytical method are consistent with those obtained by the numerical method, where the relative errors are less than 10%, indicating that the simplified method is feasible in engineering applications. The leakage amount from the core is far less than that from the reservoir bottom with the condition of a normal pool level, which reflects the good anti-seepage performance of the core. To reduce the leakage amount in case of a determined hydrogeological condition, more targeted engineering measures should be taken to control the seepage flux of the reservoir bottom, such as laying face slab for retaining water, geomembranes, etc.

Conclusions
To obtain the leakage amount of the clay core dam in a short time, we proposed a simplified analytical method. The proposed method is composed of two sets of formulae according to different conditions of the position of the groundwater level, i.e., below or above the reservoir bottom. Both formulae can be used to calculate the seepage flux of the dam body and the reservoir bottom, which are the main leakage zones of a clay core dam. To prove the rationality of the proposed method, a numerical model was built using Geostudio 2020 SEEP/W to calculate the discharge of the clay core dam, where the relative errors between the results are no more than 10%, indicating the rationality of the analytical method. In leakage estimation of the clay core dam of the Lianghekou power station, the results obtained by the proposed analytical method and the numerical method are close, and the relative errors are no more than 10%, showing the feasibility of the proposed method in engineering applications.
In engineering applications, the proposed method can help obtain the leakage The calculation results of both numerical and analytical methods for the leakage of the clay core dam are presented in Table 7.

Conclusions
To obtain the leakage amount of the clay core dam in a short time, we proposed a simplified analytical method. The proposed method is composed of two sets of formulae according to different conditions of the position of the groundwater level, i.e., below or above the reservoir bottom. Both formulae can be used to calculate the seepage flux of the dam body and the reservoir bottom, which are the main leakage zones of a clay core dam. To prove the rationality of the proposed method, a numerical model was built using Geo-studio 2020 SEEP/W to calculate the discharge of the clay core dam, where the relative errors between the results are no more than 10%, indicating the rationality of the analytical method. In leakage estimation of the clay core dam of the Lianghekou power station, the results obtained by the proposed analytical method and the numerical method are close, and the relative errors are no more than 10%, showing the feasibility of the proposed method in engineering applications.
In engineering applications, the proposed method can help obtain the leakage amount and evaluate the rationality of the clay core dam quickly in the design state without the tedious and time-consuming modeling process of numerical methods. It is of vital importance to estimate the leakage amount of clay core dams before the engineering is commenced, which not only relates to the safe operation of the project, but is also related to the economic benefits of the project. The proposed analytical method can offer a solution for leakage amount estimation by striking a balance between precision and efficiency.