Numerical Simulation Three-Dimensional Nonlinear Seepage in a Pumped-Storage Power Station : Case Study

Due to high water pressure in the concrete reinforced hydraulic tunnels, surrounding rocks are confronted with nonlinear seepage problem in the pumped storage power station. In this study, to conduct nonlinear seepage numerical simulation, a nonlinear seepage numerical model combining the Forchheimer nonlinear flow theory, the discrete variational inequality formulation of Signorini’s type and an adaptive penalized Heaviside function is established. This numerical seepage model is employed to the seepage analysis of the hydraulic tunnel surrounding rocks in the Yangjiang pumped-storage power station, which is the highest water pressure tunnel under construction in China. Moreover, the permeability of the surrounding rocks under high water pressure is determined by high pressure packer test and its approximate analytical model. It is shown that the flow in the surrounding rocks is particularly prone to become nonlinear as a result of the high flow velocities and hydraulic gradients in the nearby of the seepage-control measures and the high permeability fault. The nonlinear flow theory generates smaller flow rate than the Darcy flow theory. With the increase of nonlinear flow, this observation would become more remarkable.


Introduction
Electricity demand in China has risen speedily and reached an unprecedented level due to the rapid economic growth and modern development [1].In China, primary electricity sources include the fossil fuel, hydropower, nuclear, and other renewable energy sources (e.g., wind, solar, biomass, wave and geothermal energy), which have generated total power capacity 5649.6 TWh in 2015 [2,3].Considering the economic, technical and environmental benefits of the hydropower, most countries give priority to its development.China has the richest water resource in the world, with maximum exploitable hydropower resources of 694 GW and currently economic exploitable potential 380-400 GW [4].However, due to the unbalanced hydropower distribution and economic development in China, the provinces in Southeast China, including Guangdong, Shanghai, Jiangsu, Zhejiang, Anhui, Fujian, Jiangxi and Shandong, accounted for more than 50% of electricity consumption in China, while Southwest China has the most fruitful hydropower resources in this country, which includes four provinces: Sichuan, Yunnan, Tibet and Guizhou.Pumped-storage power station offers a technically and economically feasible solution to the problem of the hydropower situation unbalance, which was used as early as 1890 in Italy and Switzerland, and promoted in the developed country [5].For instance, in the USA, Japan and the EU, the installed capacity of pumped-storage power reached 2.14%, 8.70% and 3.35% of their total installed capacity by 2010, respectively [6][7][8].
As shown in Figure 1, a typical pumped-storage power station usually contain a lower and an upper reservoir, a powerhouse and hydraulic tunnels.During off-peak electricity demand hours, the pumped-storage power station stores electricity by moving water from a lower to an upper reservoir [9][10][11].Electrical energy is converted to potential energy and stored in the upper elevation.And during periods of peak hours, the stored water is released back through the turbines and converted back to electricity like a conventional hydropower station.In the abovementioned process, the amount of energy stored depends on the height difference between the two reservoirs and the total volume of water stored.As part of its energy transition strategy, China has accelerated the construction of pumped-storage power station since 2010 [12].In the latest five year plan, China's government has set targets for pumped hydropower station capacity of 30 GW by 2015 to 70 GW by 2020, accounting for 3-5% of the total installed generation capacities in the country.As listed in Table 1, more than 10 pumped-storage power stations had been completed in the recent five years, such as: Bailianhe (1200 MW, Hubei), Baoquan (1200MW, Henan), Xianyou (1200 MW, Fujian), Qingyuan (1200 MW, Gangdong), and Xianju (1500 MW, Zhejiang) [13][14][15].Moreover, more than 20 pumped-storage power stations are under construction in other 10 provinces in China.
As shown in Figure 1, a typical pumped-storage power station usually contain a lower and an upper reservoir, a powerhouse and hydraulic tunnels.During off-peak electricity demand hours, the pumped-storage power station stores electricity by moving water from a lower to an upper reservoir [9][10][11].Electrical energy is converted to potential energy and stored in the upper elevation.And during periods of peak hours, the stored water is released back through the turbines and converted back to electricity like a conventional hydropower station.In the abovementioned process, the amount of energy stored depends on the height difference between the two reservoirs and the total volume of water stored.As part of its energy transition strategy, China has accelerated the construction of pumped-storage power station since 2010 [12].In the latest five year plan, China's government has set targets for pumped hydropower station capacity of 30 GW by 2015 to 70 GW by 2020, accounting for 3-5% of the total installed generation capacities in the country.As listed in Table 1, more than 10 pumped-storage power stations had been completed in the recent five years, such as: Bailianhe (1200 MW, Hubei), Baoquan (1200MW, Henan), Xianyou (1200 MW, Fujian), Qingyuan (1200 MW, Gangdong), and Xianju (1500 MW, Zhejiang) [13][14][15].Moreover, more than 20 pumped-storage power stations are under construction in other 10 provinces in China.Seepage control is one of the key technical issues in deep buried hydraulic tunnel and large-scale powerhouse of a pumped-storage power station.In those abovementioned pumped-storage power stations, concrete reinforced hydraulic tunnels were usually adopted, which are subjected to a maximum hydrostatic pressure of more than 2.0 MPa on the inner surface.The water pressure will inevitably induce tensile cracks in the concrete linings, and Table 2 lists the amount of leakage out of concrete reinforced hydraulic tunnels in the pumped-storage power stations in China [16][17][18][19], U.S. [20] and Norway [21], which results in enormouspower loss day by day [22].Consequently, there is an urgent need in understanding the nonlinear groundwater flow behaviors through porous and fractured rock, which is commonly encountered in various  Seepage control is one of the key technical issues in deep buried hydraulic tunnel and large-scale powerhouse of a pumped-storage power station.In those abovementioned pumped-storage power stations, concrete reinforced hydraulic tunnels were usually adopted, which are subjected to a maximum hydrostatic pressure of more than 2.0 MPa on the inner surface.The water pressure will inevitably induce tensile cracks in the concrete linings, and Table 2 lists the amount of leakage out of concrete reinforced hydraulic tunnels in the pumped-storage power stations in China [16][17][18][19], U.S. [20] and Norway [21], which results in enormous power loss day by day [22].Consequently, there is an urgent need in understanding the nonlinear groundwater flow behaviors through porous and fractured rock, which is commonly encountered in various engineering applications, such as underground tunneling, geothermal energy extraction, hydrocarbon production, and hazardous wastes isolation [23][24][25].For instance, non-Darcy flow behavior in the near-well region was simulated with a nonlinear finite element model by Zhang and Xing [23].A numerical model of three-dimensional unstable and nonlinear seepage was established, and nonlinear seepage movement characteristics of fluid in underground coal gasification (UCG) were studied by Yang [24].Samanta et al. conducted direct numerical simulations of the fully turbulent flow through a porous square duct to study the effect of the permeable wall on the secondary cross-stream flow [26].Li et al. presented a systematic performance assessment of the seepage control system designed for the underground caverns by finite element numerical modeling, with particular concerns on quantitative determination of the excavation-induced variation in hydraulic conductivity of the surrounding rock masses [27].However, under the high water pressure environment, nonlinear flow behaviors would occur in the surrounding rock of hydraulic tunnel, which is rarely reported by previous literatures.Seepage flow modeling is a very important step before system design, simulation and optimization.Given the fact that the high water pressure in the concrete reinforced manifolds, surrounding rocks are confronted with seepage stability problem in the pumped storage power station.In this study, the equivalent continuum finite element method (FEM) is employed in this study.A nonlinear flow numerical model combining the Forchheimer nonlinear flow theory, the discrete variational inequality formulation of Signorini's type and an adaptive penalized Heaviside function is established, which is employed to the seepage analysis of the hydraulic tunnel and powerhouse surrounding rocks in the Yangjiang pumped-storage power station.Moreover, the permeability of the rock mass under high water pressure is determined by high pressure packer test and its approximate analytical model.

Governing Equations
Given the fact that the high water pressure, nonlinear flow would occur in the concrete reinforced hydraulic tunnels, and this type of flow can be adequately described by the Forchheimer equation (Forchheimer 1901) [23]: where ∇h is the pressure gradient (or pressure drop), µ the dynamic viscosity of the fluid, ρ w is the fluid density, v is the average velocity, k is the intrinsicpermeability, and β is the nonlinear coefficient dependent on the properties of the fractured rock.Because of inconsistency in definitions and thus in critical values, no widely accepted criterion for nonlinear flow in porous media is available.The Forchheimer number F 0 is defined as the ratio of nonlinear (the second term in Equation ( 1)) to linear (the first term in Equation ( 1)) pressure losses in the Forchheimer's law, and it represents the ratio of the pressure gradient required to overcome inertial forces to that of viscous forces.The Forchheimer number F 0 is adopted as criterion for nonlinear flow in this study because of the clear physical meaning of variables involved.
Assuming the fractured rock is non-deformable and the fluid has a constant effective compressibility, substituting Equation ( 2) into (1) yields.
As demonstrated in Figure 2, seepage flow through rock domain Ω = Ω w ∪ Ω d is then governed by the following equation of continuity.
Energies 2019, 12, x FOR PEER REVIEW 4 of 16 inertial forces to that of viscous forces.The Forchheimer number F0is adopted as criterion for nonlinear flow in this studybecause of the clearphysical meaning of variables involved.
Assuming the fractured rock is non-deformable and the fluid has a constant effective compressibility, substituting Equation ( 2) into (1) yields.
As demonstrated in Figure 2, seepage flow through rock domain Ω=Ωw∪Ωd is then governed by the following equation of continuity.
Figure2.Sketch of seepage flow through an underground cavern.
Starting from Equations ( 3) and ( 4) and following the discrete variational inequality formulation of Signorini's type proposed by Chen et al. [28], one obtains the governing equation for nonlinear flow in fractured rock equation: is the permeability of the fractured rock.
Equation ( 5) is subjected to the following boundary conditions:


The water head boundary condition where  is the prescribed water head onΓφ, is the total water head, zthe vertical coordinate, pthe pore water pressure,γwthe unit weightof water.


The flux boundary condition Starting from Equations ( 3) and ( 4) and following the discrete variational inequality formulation of Signorini's type proposed by Chen et al. [28], one obtains the governing equation for nonlinear flow in fractured rock equation: in which K = µ k ρ w g is the permeability of the fractured rock.Equation ( 5) is subjected to the following boundary conditions:

•
The water head boundary condition where φ is the prescribed water head on Γ ϕ , φ = z + p/γ w is the total water head, z the vertical coordinate, p the pore water pressure, γ w the unit weight of water.

•
The flux boundary condition where q is the prescribed flux on Γ q , and n the outward unit normal vector to the boundary.For an impermeable boundary, q = 0.

•
The boundary condition of Signorini's type on the seepage surface where Γ s is the potential seepage boundary.Obviously, on section GHIJ, ϕ = z and q n ≤ 0; while on sections AFED and GKJ, ϕ < z and q n = 0. ϕ = z and q n = 0 are satisfied at seepage points G and J.

•
The boundary condition on the phreatic surface where Γ f ≡ {(x, y, z|ϕ = z)} is the phreatic surface, an interface between Ω w and Ω d .

Solution Procedure and Finite Element Numerical Methods
The mathematical statement for a discrete version of the iterative formulation is given as follows: find a vector VI , the following inequality holds: with where l is the iterative step, j the total number of nodal points in the finite element mesh, B the geometrical matrix of the finite element model, and H λ the adaptive penalized Heaviside function introduced to evade numerical instability and mesh dependency, given by Chen et al. [28].
where λ 1 and λ 2 are two parameters associated with each element.The parameter λ 1 is defined as the vertical distance between the lowest integration point and the lowest node, and λ 2 the vertical distance between the highest integration point and the highest node, in the FEM model concerned.The symbol ζ is introduced to scale values of parameters λ 1 and λ 2 , so that stronger convergence criteria could be ensured for highly nonlinear problems with coarse meshes.Moreover, the introduction of ζ in the suggested range does not significantly influence the behavior of the penalized Heaviside function, but results in very good numerical convergence [28].
The discretized PDE formulation given in Equation (10) was implemented in a FEM code, THYME3D, initially developed for coupled deformation/multiphase flow/thermal transport analysis in geological porous media [29].
In the iterative process, the following convergence criteria is used: in which φ l and φ l+1 denote the water head at the l and l + 1 step, respectively.The symbols ε 1 and ε 2 denote the user-specified error tolerances and take the values of ε 1 = 10 −3 and ε 2 = 10 −4 in this study.

General Description and Geological Condition
Yangjiang pumped-storage power station is under construction in Yangjiang County, Guangdong Province, China, which includes an upper reservoir, a lower reservoir and hydropower system with installed capacity and annual output were 2400 MW and 3432GWh, respectively.The upper reservoir is located in Efengling provincial nature reserve area, with effective catchment area of 7.54 km 2 and normal pool level 773.7 m.The lower reservoir is located in Litian River, with effective catchment area of 15.94 km 2 and normal pool level 103.7 m.As illustrated in Figure 3, the diversion tunnel, 7.4 m in diameter, is composed of an upper horizontal section, a middle horizontal section, a lower horizontal section and two inclined sections that connect the horizontal ones.At the lower horizontal section, the hydraulic tunnel is branched into three smaller tunnels of 3.0 m in diameter, which carry water to operate three reversible turbine/generator assemblies installed in the hydropower.The tunnels are concrete reinforced, and at the lower horizontal section, they are subjected to a maximum hydrostatic pressure of 7.99 MPa and an extra surge pressure up to 3.0 MPa on the inner surface in the condition of general operating. and in which l  and 1 l   denote the water head at the l and l+1 step, respectively.The symbols ε1 and ε2 denote the user-specified error tolerances and take the values of ε1=10 −3 and ε2=10 −4 in this study.

General Description and Geological Condition
Yangjiang pumped-storage power station is under construction in Yangjiang County, Guangdong Province, China, which includes an upper reservoir, a lower reservoir and hydropower system with installed capacity and annual output were 2400 MW and 3432GWh, respectively.The upper reservoir is located in Efengling provincial nature reserve area, with effective catchment area of 7.54 km 2 and normal pool level 773.7 m.The lower reservoir is located in Litian River, with effective catchment area of 15.94 km 2 and normal pool level 103.7 m.As illustrated in Figure 3, the diversion tunnel, 7.4 m in diameter, is composed of an upper horizontal section, a middle horizontal section, a lower horizontal section and two inclined sections that connect the horizontal ones.At the lower horizontal section, the hydraulic tunnel is branched into three smaller tunnels of 3.0 m in diameter, which carry water to operate three reversible turbine/generator assemblies installed in the hydropower.The tunnels are concrete reinforced, and at the lower horizontal section, they are subjected to a maximum hydrostatic pressure of 7.99 MPa and an extra surge pressure up to 3.0 MPa on the inner surface in the condition of general operating.The bedrocks in the study area mainly contain sedimentary rocks of the Cretaceous Lumuwan group (K1Lm) and Indosinian granites (ηγ5 1-3 ), with K1Lm unconformably overlying ηγ5 1-3 .The Cretaceous Lumuwan group (K1Lm) is abundantly distributed in the lower reservoir area, and it consists of sandstone and conglomerate.The Indosinian granites (ηγ5 1-3 ) are mostly distributed in the upper reservoir area, and consist of chloritization granite, biotite granite and medium-grained granite.As demonstrated in Figure 4, the main structures at the site consist of five sub-vertically The bedrocks in the study area mainly contain sedimentary rocks of the Cretaceous Lumuwan group (K 1Lm ) and Indosinian granites (ηγ 5 1-3 ), with K 1Lm unconformably overlying ηγ 5 1-3 .
The Cretaceous Lumuwan group (K 1Lm ) is abundantly distributed in the lower reservoir area, and it consists of sandstone and conglomerate.The Indosinian granites (ηγ 5 1-3 ) are mostly distributed in the upper reservoir area, and consist of chloritization granite, biotite granite and medium-grained granite.As demonstrated in Figure 4, the main structures at the site consist of five sub-vertically oriented faults, f721, f717, f715, f708 and f745, with f721 striking towards N70

Proposed Seepage-Control Measures
Figure 5 illustrates seepage-control measures for the hydropower in the Yangjiang pumped storage power station, which is composed of high pressure consolidation grouting, impervious

Proposed Seepage-Control Measures
Figure 5 illustrates seepage-control measures for the hydropower in the Yangjiang pumped storage power station, which is composed of high pressure consolidation grouting, impervious curtain, galleries and drainage holes.To form a complete impervious system in the branch of hydraulic tunnels, the high pressure consolidation grouting is arranged around the hydraulic tunnels, with 6.0 m in depth, 230 m in length and 1.3 times as the static head high in grouting pressure.The "V" shaped curtain grouting galleries (A1) is arranged right above the consolidation grouting, and the curtain extends 2.0 m and 3.0 m respectively from the upstream and downstream hydraulic tunnels.The proposed space between the consolidation grouting and curtain grouting holes is 2 m, with 76mm in the diameter.At the same time, the drainage galleries are deployed in the end of the hydraulic tunnels (A2) and in the surrounding rock of the hydropower (A3~A6) to form a drainage system.Drainage holes are vertically deployed in the surrounding rock immediately behind the grouting curtain, and the drainage galleries close to the hydropower are connected by the vertical drainage holes.The drainage holes are deployed with 5 m in spacing and 120.0 m in depth, and the size of the cross-sections of the holes is assumed to be 12.0 cm in diameter.
Energies 2019, 12, x FOR PEER REVIEW 8 of 16 curtain, galleries and drainage holes.To form a complete impervious system in the branch of hydraulic tunnels, the high pressure consolidation grouting is arranged around the hydraulic tunnels, with 6.0 m in depth, 230 m in length and 1.3 times as the static head high in grouting pressure.The "V" shaped curtain grouting galleries (A1) is arranged right above the consolidation grouting, and the curtain extends 2.0 m and 3.0 m respectively from the upstream and downstream hydraulic tunnels.The proposed space between the consolidation grouting and curtain grouting holes is 2 m, with 76mm in the diameter.At the same time, the drainage galleries are deployed in the end of the hydraulic tunnels (A2) and in the surrounding rock of the hydropower (A3~A6) to form a drainage system.Drainage holes are vertically deployed in the surrounding rock immediately behind the grouting curtain, and the drainage galleries close to the hydropower are connected by the vertical drainage holes.The drainage holes are deployed with 5 m in spacing and 120.0 m in depth, and the size of the cross-sections of the holes is assumed to be 12.0 cm in diameter.

The Finite Element Model
To simulate the nonlinear seepage filed and assess the performance of the seepage-control measures in the surrounding rock of the hydraulic tunnels and hydropower, a 3-D finite element (FE) mesh, with 455,110 brick elements and 316,464 nodes in total, was generated, as shown in Figure 6.The size of the mesh is 460 m along the pumping or generating hydraulic tunnels water flow direction, and 360 m in another horizontal direction.According to Chen et al. [28], for the proposed parabolic PDE and the discrete version, the accuracy of numerical simulation is dependent on the mesh density.Therefore, to improve the calculation accuracy of the seepage field around the hydraulic tunnels and water outflow, the grid around the hydraulic tunnels was more intensive.Simultaneously, the topographic and geological features at the study area, and the seepage-control measures containing high pressure consolidation grouting, impervious curtain, galleries and drainage holes.

The Finite Element Model
To simulate the nonlinear seepage filed and assess the performance of the seepage-control measures in the surrounding rock of the hydraulic tunnels and hydropower, a 3-D finite element (FE) mesh, with 455,110 brick elements and 316,464 nodes in total, was generated, as shown in Figure 6.The size of the mesh is 460 m along the pumping or generating hydraulic tunnels water flow direction, and 360 m in another horizontal direction.According to Chen et al. [28], for the proposed parabolic PDE and the discrete version, the accuracy of numerical simulation is dependent on the mesh density.Therefore, to improve the calculation accuracy of the seepage field around the hydraulic tunnels and water outflow, the grid around the hydraulic tunnels was more intensive.Simultaneously, the topographic and geological features at the study area, and the seepage-control measures containing high pressure consolidation grouting, impervious curtain, galleries and drainage holes.

The Calculation Parameters
To obtain the permeability of the surrounding rock under high water pressure, high pressure packer test (HPPT) was conducted of the boreholes in the study area, as plotted in the Figure 4. Thereafter, on the basis of Izbash's nonlinear empirical equation, an approximate analytical model was developed especially for interpreting data from these HPPTs.
Assuming the volumetric flow rate (q) injected into the infinitesimal section to be linearly proportional to the volumetric total flow rate (Q) injected into the test interval (L) yields Further, assuming radial velocity vρof the test interval yields (17) whereρ is the radial distance from the borehole axis.
On the plane vertically intersecting the center of the test interval, the Izbash's law is similarly written as r r r where Pr is the water pressure at the distance r from the borehole, vr the tangential velocity, and λ and m is empirical coefficient in the Izbash's law.For laminar flow (m=1), the Izbash's law reduces to the Darcy's law, and for m=2, Equation ( 18) represents a fully turbulent flow.In transitional flow condition, the value of m ranges from 1 to 2. Beyond the radius of influence R0, the pressure is assumed to vanish r lim 0

The Calculation Parameters
To obtain the permeability of the surrounding rock under high water pressure, high pressure packer test (HPPT) was conducted of the boreholes in the study area, as plotted in the Figure 4. Thereafter, on the basis of Izbash's nonlinear empirical equation, an approximate analytical model was developed especially for interpreting data from these HPPTs.
Assuming the volumetric flow rate (q) injected into the infinitesimal section to be linearly proportional to the volumetric total flow rate (Q) injected into the test interval (L) yields Further, assuming radial velocity v ρ of the test interval yields where ρ is the radial distance from the borehole axis.
On the plane vertically intersecting the center of the test interval, the Izbash's law is similarly written as where P r is the water pressure at the distance r from the borehole, v r the tangential velocity, and λ and m is empirical coefficient in the Izbash's law.For laminar flow (m = 1), the Izbash's law reduces to the Darcy's law, and for m = 2, Equation ( 18) represents a fully turbulent flow.In transitional flow condition, the value of m ranges from 1 to 2. Beyond the radius of influence R 0 , the pressure is assumed to vanish lim Substituting Equation ( 17) into ( 18) and integrating Equation ( 18) yields The analytical solution for the 2D radial flow is then derived as Rearranging Equation ( 21) yields the expression for permeability k: Based on the above mentioned in situ tests and approximate analytical model, permeability of the surrounding rock is listed in Table 4.Moreover, permeability of high pressure consolidation grouting and impervious curtain is 3.68 × 10 −6 and 5.0 × 10 −6 cm/s, respectively.

The Boundary Conditions
The back analysis showed that when the water head on the upper reservoir stakes a normal water level of 773.7 m, the objective function was minimized over all the boreholes with available groundwater level measurements.Under this condition, the upstream of the FEM model water level is 770 m, and the corresponding water level in the downstream is 208.5 m.The lateral boundary at the mountain side and the base boundary of the model was assumed to be impermeable.The ground and dam surfaces above the upstream and downstream water levels, the surfaces of the hydraulic tunnels, the drainage galleries and the drainage holes, are all taken as the potential seepage boundaries satisfying the Signorini's complementary condition.

Performance Assessment of the Seepage-Control Measures
Figure 7 shows the distribution of phreatic surface and pressure head at cross-section I-I, with its location shown in Figure 6c.One observes from Figure 7 that the pressure head around the surrounding rocks of the hydraulic tunnels is rather high due to the low-permeability high pressure consolidation grouting around the tunnels.Moreover, the closer to the branch of the hydraulic tunnels, more intensive distribution of the pressure head.Both abovementioned results showed that the high pressure consolidation grouting plays a very good anti-seepage effect.Figure 8 shows the distribution of phreatic surface and pressure head at cross-section II-II, with its location shown in Figure 6c.It demonstrates that the impervious curtain at the upstream of the branch, and the drainage facilities and holes has a dramatic impact on the seepage flow and can sharply depress the free surface, resulting in a remarkable fall of external water pressure above the hydraulic tunnels.These phenomenaare also observed in Figure 9, where water pressure decreases to about 200 m after the seepage-control measures around the hydraulic tunnels.The simulation results show that the suggested seepage-control measures is effective in lowering the groundwater level and reducing the pore water pressure in the surrounding rocks around the hydraulic tunnels, and the possible concentrated flow through Faults f721, f745, f708, f715 and f717 is effectively avoided.These phenomenaare also observed in Figure 9, where water pressure decreases to about 200 m after the seepage-control measures around the hydraulic tunnels.The simulation results show that the suggested seepage-control measures is effective in lowering the groundwater level and reducing the pore water pressure in the surrounding rocks around the hydraulic tunnels, and the possible concentrated flow through Faults f721, f745, f708, f715 and f717 is effectively avoided.These phenomenaare also observed in Figure 9, where water pressure decreases to about 200 m after the seepage-control measures around the hydraulic tunnels.The simulation results show that the suggested seepage-control measures is effective in lowering the groundwater level and reducing the pore water pressure in the surrounding rocks around the hydraulic tunnels, and the possible concentrated flow through Faults f721, f745, f708, f715 and f717 is effectively avoided.

Nonlinear Seepage Characteristic in the Surrounding Rocks of the High Water Pressure Hydraulic Tunnels
To illustrate the nonlinear flow characteristic in the surrounding rocks of the hydraulic tunnels, Figures 10 and 11 plot the distribution of the Forchheimer number F0 at cross-section of the A2 and A3 drainage galleries, respectively.One observes from Figure 10 that the nonlinear seepage mainly occurs in the areas nearby faults and seepage-control measures.The maximum values were 0.7 and 1 at the bottom of drainage holes and upstream of A3 drainage galleries, where the nonlinear flow is the most strongly.That is because the hydraulic gradient in the surrounding rocks near the seepage-control measures are so large, which would lead nonlinear flow easily.In addition, special attention must be focused on the fault f708 and f717, where the high water pressure is firstly occurs at the f708 and rapidly extend to f717.Then, a potential seepage channel will growth, which may result the leakage in the surrounding rocks of the hydraulic tunnels.In order to demonstrate the nonlinear effect on water inflow and outflow, the flow rate at the branch of the hydraulic tunnels, A1-A6 drainage galleries and drainage holes are calculated by the linear and nonlinear seepage theories.

Nonlinear Seepage Characteristic in the Surrounding Rocks of the High Water Pressure Hydraulic Tunnels
To illustrate the nonlinear flow characteristic in the surrounding rocks of the hydraulic tunnels, Figures 10 and 11 plot the distribution of the Forchheimer number F 0 at cross-section of the A2 and A3 drainage galleries, respectively.One observes from Figure 10 that the nonlinear seepage mainly occurs in the areas nearby faults and seepage-control measures.The maximum values were 0.7 and 1 at the bottom of drainage holes and upstream of A3 drainage galleries, where the nonlinear flow is the most strongly.That is because the hydraulic gradient in the surrounding rocks near the seepage-control measures are so large, which would lead nonlinear flow easily.In addition, special attention must be focused on the fault f708 and f717, where the high water pressure is firstly occurs at the f708 and rapidly extend to f717.Then, a potential seepage channel will growth, which may result the leakage in the surrounding rocks of the hydraulic tunnels.In order to demonstrate the nonlinear effect on water inflow and outflow, the flow rate at the branch of the hydraulic tunnels, A1-A6 drainage galleries and drainage holes are calculated by the linear and nonlinear seepage theories.
As listed in Table 5, the leakage is calculated by nonlinear flow theory are significantly less than the linear flow, such as leakage are 32.9 and 40.5 L/s at the branch of the hydraulic tunnels, and corresponding power loss are 2189 and 2694 kW•h, respectively.Moreover, the difference is getting bigger as the nonlinear flow stronger, and this is because the calculated flow velocity by nonlinear flow theory is smaller than the linear theory.As listed in Table 5, the leakage is calculated by nonlinear flow theory are significantly less than the linear flow, such as leakage are 32.9 and 40.5 L/s at the branch of the hydraulic tunnels, and corresponding power loss are 2189 and 2694 kW•h, respectively.Moreover, the difference is getting  As listed in Table 5, the leakage is calculated by nonlinear flow theory are significantly less than the linear flow, such as leakage are 32.9 and 40.5 L/s at the branch of the hydraulic tunnels, and corresponding power loss are 2189 and 2694 kW•h, respectively.Moreover, the difference is getting

Conclusions
A nonlinear seepage numerical model is built in this study for numerical simulation the nonlinear seepage in a pumped-storage power station.This model contains the Forchheimer nonlinear flow theory, the discrete variational inequality formulation of Signorini's type and an adaptive penalized Heaviside function, which is employed to demonstrate the performance assessment of the seepage-control measures and nonlinear flow characteristic in the surrounding rocks of the hydraulic tunnels.The major conclusions of this case study are summarized as follows: (1) The simulation results show that the suggested seepage-control measures is effective in lowering the groundwater level and reducing the pore water pressure in the surrounding rocks around the hydraulic tunnels, and the possible concentrated flow through Faults f721, f745, f708, f715 and f717 is effectively avoided.
(2) Flow in the surrounding rocks is particularly prone to become nonlinear as a result of the high flow velocities and hydraulic gradients in the nearby of the seepage-control measures and the high permeability fault.The flow rate derived from nonlinear flow theory is smaller than the Darcy flow theory, what is more remarkable as the increase of nonlinear flow.Therefore, it is necessary to adopt nonlinear flow theory in the seepage analysis in the condition of high water pressure.
(3) With the development of the pumped storage power stations in China, the impervious design and safety control of high head hydropower projects are the key technical problems that need to solved urgently.Preliminary work on the nonlinear seepage analysis of surrounding rock has carried out in this study, more efforts are needed in the theory, analysis method and control measures of nonlinear seepage in geotechnical engineering.

Figure 2 .
Figure 2. Sketch of seepage flow through an underground cavern.

Figure 3 .
Figure 3.Layout of the diversion tunnel in the Yangjiang pumped-storage power station.

Figure 3 .
Figure 3. Layout of the diversion tunnel in the Yangjiang pumped-storage power station.

Figure 4 .
Figure 4.Layout of the geological map in the Yangjiang pumped-storage power station.

Figure 4 .
Figure 4. Layout of the geological map in the Yangjiang pumped-storage power station.

Figure 5 .
Figure 5. Illustration of seepage-control measures for the hydropower in the Yangjiang pumped storage power station.

Figure 5 .
Figure 5. Illustration of seepage-control measures for the hydropower in the Yangjiang pumped storage power station.

Figure 7 .Figure 7 .
Figure 7. Phreatic surfaces and pressure distribution in the surrounding rocks at cross-section I-I.

Figure 7 .Figure 8 .
Figure 7. Phreatic surfaces and pressure distribution in the surrounding rocks at cross-section I-I.

Figure 8 .Figure 9 .
Figure 8. Phreatic surfaces and pressure distribution in the surrounding rocks at cross-section II-II.

Figure 9 .
Figure 9. Phreatic surfaces and pressure distribution in the surrounding rocks at cross-section III-III.

Figure 10 .
Figure 10.The F0 contours at cross-section of the A2 drainage galleries.

Figure 11 .
Figure 11.The F0 contours at cross-section of the A3 drainage galleries.

Figure 10 .Figure 10 .
Figure 10.The F 0 contours at cross-section of the A2 drainage galleries.

Figure 11 .
Figure 11.The F0 contours at cross-section of the A3 drainage galleries.

Figure 11 .
Figure 11.The F 0 contours at cross-section of the A3 drainage galleries.

Table 1 .
Major pumped-storage power stations completed in China.

Table 1 .
Major pumped-storage power stations completed in China.

Table 2 .
Leakage events occurred in the pumped-storage power stations.

Table 3 .
Key faults information selected in the study area.

Table 3 .
Key faults information selected in the study area.

Table 4 .
Permeability of the strata and faults at the study area (unit: cm/s).

Table 5 .
Comparison of the flow rate for linear and nonlinear flow theories (L/min).