Seepage simulation of high concrete-faced rockfill dams based on generalized equivalent continuum model

This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.

Monitoring models for base flow effect and daily variation of dam seepage elements considering time lag effect

Affected by external environmental factors and evolution of dam performance, dam seepage behavior shows nonlinear time-varying characteristics. In this study, to predict and evaluate the long-term development trend and short-term fluctuation of the dam seepage behavior, two monitoring models were developed, one for the base flow effect and one for daily variation of dam seepage elements. In the first model, to avoid the influence of the time lag effect on the evaluation of seepage variation with the time effect component of seepage elements, the base values of the seepage element and the reservoir water level were extracted using the wavelet multi-resolution analysis method, and the time effect component was separated by the established base flow effect monitoring model. For the development of the daily variation monitoring model for dam seepage elements, all the previous factors, of which the measured time series prior to the dam seepage element monitoring time may have certain influence on the monitored results, were considered. Those factors that were positively correlated with the analyzed seepage element were initially considered to be the support vector machine(SVM) model input factors, and then the SVM kernel function-based sensitivity analysis was performed to optimize the input factor set and establish the optimized daily variation SVM model. The efficiency and rationality of the two models were verified by case studies of the water level of two piezometric tubes buried under the slope of a concrete gravity dam.Sensitivity analysis of the optimized SVM model shows that the influences of the daily variation of the upstream reservoir water level and rainfall on the daily variation of piezometric tube water level are processes subject to normal distribution.

Seepage safety monitoring model for an earth rock dam under influence of high-impact typhoons based on particle swarm optimization algorithm

Extreme hydrological events induced by typhoons in reservoir areas have presented severe challenges to the safe operation of hydraulic structures. Based on analysis of the seepage characteristics of an earth rock dam, a novel seepage safety monitoring model was constructed in this study. The nonlinear influence processes of the antecedent reservoir water level and rainfall were assumed to follow normal distributions. The particle swarm optimization （PSO） algorithm was used to optimize the model parameters so as to raise the fitting accuracy. In addition, a mutation factor was introduced to simulate the sudden increase in the piezometric level induced by short-duration heavy rainfall and the possible historical extreme reservoir water level during a typhoon. In order to verify the efficacy of this model, the earth rock dam of the Siminghu Reservoir was used as an example. The piezometric level at the SW1-2 measuring point during Typhoon Fitow in 2013 was fitted with the present model, and a corresponding theoretical expression was established. Comparison of fitting results of the piezometric level obtained from the present statistical model and traditional statistical model with monitored values during the typhoon shows that the present model has a higher fitting accuracy and can simulate the uprush feature of the seepage pressure during the typhoon perfectly.

Three-Dimensional Tidal Model and Its Application to Numerical Simulation of Water Quality in Coastal Waters

The turbulence mechanism plays an important part in the mixing process and momentum transfer of turbulence. A three-dimensional Prandtl mixing length tidal model has been developed to simulate tidal flows and water quality. The eddy viscosities and diffusivities are computed from the Prandtl mixing length model. In order to model the water quality of an estuary or coastal area many interdependent processes need to be simulated. These may be conveniently separated into three main groups: transport and mixing processes, biochemical interaction of water quality variables and the utilization and re-cycling of nutrients by living matter. The model simulates full oxygen and nutrient balance, primary productivity and the transport, reaction mechanism and fate of pollutants over tidal time-scales. The model is applied to numerical simulation of tidal flows and water quality in Dalian Bay. The model has been calibrated against a limited data set of historical water quality observations and in general demonstrates excellent agreement with all available data.

The three-dimension model for the rock-breaking mechanism of disc cutter andanalysis of rock-breaking forces

To study the rock deformation with three- dimensional model under rolling forces of disc cutter, by car- rying out the circular-grooving test with disc cutter rolling around on the rock, the rock mechanical behavior under rolling disc cutter is studied, the mechanical model of disc cutter rolling around the groove is established, and the the- ory of single-point and double-angle variables is proposed. Based on this theory, the physics equations and geometric equations of rock mechanical behavior under disc cutters of tunnel boring machine （TBM） are studied, and then the bal- ance equations of interactive forces between disc cutter and rock are established. Accordingly, formulas about normal force, rolling force and side force of a disc cutter are de- rived, and their validity is studied by tests. Therefore, a new method and theory is proposed to study rock- breaking mech- anism of disc cutters.

Groundwater flow through fractured rocks and seepage control in geotechnical engineering: Theories and practices

Groundwater flow through fractured rocks has been recognized as an important issue in many geotechnical engineering practices.Several key aspects of fundamental mechanisms,numerical modeling and engineering applications of flow in fractured rocks are discussed.First,the microscopic mechanisms of fluid flow in fractured rocks,especially under the complex conditions of non-Darcian flow,multiphase flow,rock dissolution,and particle transport,have been revealed through a com-bined effort of visualized experiments and theoretical analysis.Then,laboratory and field methods of characterizing hydraulic properties(e.g.intrinsic permeability,inertial permeability,and unsaturated flow parameters)of fractured rocks in different flow regimes have been proposed.Subsequently,high-performance numerical simulation approaches for large-scale modeling of groundwater flow in frac-tured rocks and aquifers have been developed.Numerical procedures for optimization design of seepage control systems in various settings have also been proposed.Mechanisms of coupled hydro-mechanical processes and control of flow-induced deformation have been discussed.Finally,three case studies are presented to illustrate the applications of the improved theoretical understanding,characterization methods,modeling approaches,and seepage and deformation control strategies to geotechnical engi-neering projects.

A Well Productivity Model for Multi-Layered Marine and Continental Transitional Reservoirs with Complex Fracture Networks

Using the typical characteristics of multi-layered marine and continental transitional gas reservoirs as a basis,a model is developed to predict the related well production rate.This model relies on the fractal theory of tortuous capillary bundles and can take into account multiple gas flow mechanisms at the micrometer and nanometer scales,as well as the flow characteristics in different types of thin layers(tight sandstone gas,shale gas,and coalbed gas).Moreover,a source-sink function concept and a pressure drop superposition principle are utilized to introduce a coupled flow model in the reservoir.A semi-analytical solution for the production rate is obtained using a matrix iteration method.A specific well is selected for fitting dynamic production data,and the calculation results show that the tight sandstone has the highest gas production per unit thickness compared with the other types of reservoirs.Moreover,desorption and diffusion of coalbed gas and shale gas can significantly contribute to gas production,and the daily production of these two gases decreases rapidly with decreasing reservoir pressure.Interestingly,the gas production from fractures exhibits an approximately U-shaped distribution,indicating the need to optimize the spacing between clusters during hydraulic fracturing to reduce the area of overlapping fracture control.The coal matrix water saturation significantly affects the coalbed gas production,with higher water saturation leading to lower production.

Seepage-stress coupling constitutive model of anisotropic soft rock

To provide a seepage-stress coupling constitutive model that can directly describe the seepage-stress coupling relationship, a series of one-dimensional seepage-stress coupling tests on two kinds of soft rock (argillaceous siltstone and brown mudstone) were performed by using an MTS-815.02 tri-axial rock mechanics test system, with which the stress-strain curves according to the seepage variation were obtained. Based on the experimental results and by employing Hooke's law, the formulation of the coefficient of strain-dependent permeability was presented and introduced to establish a coupling model. In addition, the mathematical expression and the incremental formulation for coupling model were advanced, in which five parameters that can be respectively determined by using the experimental results were included. The calculated results show that the proposed coupling model is capable of simulating the stress-strain relationship with considering the seepage-stress coupling in the nonlinear elastic stage of two kinds of soft rock.

Coupling model for assessing anti-seepage behavior of curtain under dam foundation

The seepage under a dam foundation is mainly controlled by the performance of the curtain.Its anti-seepage behavior may be weakened by the long term physic-chemical actions from groundwater.According to seepage hydraulics and geochemistry theory,a coupling model for assessing the behavior of the curtain under a dam foundation is set up,which consists of seepage module,solute transport module,geochemistry module and curtain erosion module,solved by FEM.A case study was carried out.The result shows that the curtain efficiency is weakened all the time.Aqueous calcium from the curtain is always in dissolution during the stress period for simulation,which leads to the increasing amount in groundwater reaching 846.35-865.312 g/m3.Within the domain,reaction extent differs in different parts of the curtain.The dissolution of Ca(OH)2 accounts to 877.884 g/m3 near the bottom and is much higher than that of the other parts.The erosion is much more serious near the bottom of the curtain than the other parts,which is the same to the upstream and downstream.Calcium dissolution is mainly controlled by hydraulic condition and dispersion,and it varies in a non-linear way within the domain.

Applying seepage modeling to improve sediment yield predictions in contour ridge systems

Contour ridge systems may lead to seepage that could result in serious soil erosion. Modeling soil erosion under seepage conditions in a contour ridge system has been overlooked in most current soil erosion models. To address the importance of seepage in soil erosion modeling, a total of 23 treatments with 3 factors, row grade, field slope and ridge height, in 5 gradients were arranged in an orthogonal rotatable central composite design. The second-order polynomial regression model for predicting the sediment yield was improved by using the measured or predicted seepage discharge as an input factor, which increased the coefficient of determination(R^2) from 0.743 to 0.915 or 0.893. The improved regression models combined with the measured seepage discharge had a lower P(0.007) compared to those combined with the predicted seepage discharge(P=0.016). With the measured seepage discharge incorporated, some significant(P<0.050) effects and interactions of influential factors on sediment yield were detected, including the row grade and its interactions with the field slope, ridge height and seepage discharge, the quadratic terms of the field slope and its interactions with the row grade and seepage discharge. In the regression model with the predicted seepage discharge as an influencing factor, only the interaction between row grade and seepage discharge significantly affected the sediment yield. The regression model incorporated with predicted seepage discharge may be expressed simply and can be used effectively when measured seepage discharge data are not available.

Research on the temperature field of a partially freezing sand barrier with groundwater seepage

To study the distribution characteristics and variation regularity of the temperature field during the process of seepage freezing,a simulated-freezing test with seepage of Xuzhou sand was completed by using a model test developed in-house equipment.By means of three group freezing tests with different seepage velocities,we discovered the phenomenon of the asymmetry of the temperature field under the influence of seepage.The temperature upstream was obviously higher than that downstream.The temperature gradient upstream was also steeper than that downstream.With a higher seepage velocity,the asymmetry of the temperature field is more pronounced.The asymmetry for the interface temperature profile is more strongly manifest than for the main surface temperature profile.The cryogenic barrier section is somewhat"heartshaped".With the increasing velocity of the seepage flow,the cooling rate of the soil decreases.It takes much time to reach the equilibrium state of the soil mass.In our study,seepage flow velocities of 0 m/d,7.5 m/d,and 15 m/d showed the soilcooling rate of 4.35°C/h,4.96°C/h,and 1.72°C/h,respectively.

Test Research of Seepage Monitoring Based on Distributed Optical Fiber

In this paper,a theoretical model of temperature and velocity of a fiber is derived.A test model simulating seepage indoor is designed.The optical fiber heating temperatures under different compaction degree and seepage velocities are measured through applying AC voltage on the optical fiber.The analyzing results show that the optical fiber heating temperature and seepage velocity are related in quadratic function.The quantitative relations of optical fiber temperature and seepage velocity under different soil types and compaction degree are fitted.Analysis on how the compaction degree influences the relation of optical fiber temperature and seepage velocity shows that with the increase of compaction degree,optical fiber heating temperature will gradually decline.The influence of soil type on fiber heating temperature is very complex.In practice,according to the characteristics of the soil,determining quantitative relationship and implementing quantitative monitoring of the seepage velocity are needed.

Three-dimensional frequency-domain full waveform inversion based on the nearly-analytic discrete method

The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.

Seepage Mechanism and Transient Pressure Analysis of Shale Gas

The current research of nonlinear seepage theory of shale-gas reservoir is still in its infancy. According to the characteristics of shale gas in adsorption-desorption, diffusion, slippage and seepage during accumulation, migration and production, a mathematical model of unstable seepage in dual-porosity sealed shale-gas reservoir was developed while considering Knudsen diffusion, slip-flow effect and Langmuir desorption effect. By solving the model utilizing the Stehfest numerical inversion and computer programming in Laplace space, several typical curves of bottomhole pressure were obtained. In this paper, we discussed the effects of several parameters on the pressure dynamics, i.e. storativity ratio, Langmuir volume, Langmuir pressure, adsorption-desorption, tangential momentum accommodation coefficient, flow coefficient, boundary. The results show that the desorbed gas extends the time for fluid to flow from matrix system to fracture system;the changes of Langmuir volume and Langmuir pressure associated with desorption and adsorption effect are the internal causes of the storativity ratio change;when the tangential momentum accommodation coefficient decreases, the time for pressure wave to spread to the border reduces;interporosity flow coefficient determines the occurrence time of the transition stage;boundary range restricts the time for pressure wave to spread to the border.

An Approach for Quantifying the Influence of Seepage Dissolution on Seismic Performance of Concrete Dams

Many concrete dams seriously suffer from long-term seepage dissolution,and the induced mechanical property deterioration of concrete may significantly affect the structural performance,especially the seismic safety.An approach is presented in this paper to quantify the influence of seepage dissolution on seismic performance of concrete dams.To connect laboratory test with numerical simulation,dissolution tests are conducted for concrete specimens and using the cumulative relative leached calcium as an aging index,a deterioration model is established to predict the mechanical property of leached concrete in the first step.A coupled seepage-calcium dissolutionmigrationmodel containing two calculation modes is proposed to simulate the spatially non-uniformdeterioration of concrete dams.Based on the simulated state of a roller compacted concrete dam subjected to 100 years of seepage dissolution,seismic responses of the damare subsequently analyzed.During which the nonlinear cracking of concrete,the radiation damping of the far-field foundation is considered.Research results show that seepage dissolution will seriously weaken the seismic safety of concrete dams because of the dissolution-induced decrease of effective thickness of the dam body.The upstream surface,dam toe and gallery wall suffer from a large degree of dissolution,whereas it is minimal and basically the same inside the dam body,at a degree of 0.19%within 100 years.The horizontal displacements of dam crest under the design static load and fortification against earthquake increase by 6.9%and 21.9%,respectively,and the dissolution-induced seismic cracking leads to the failure of dam anti-seepage system.This study can provide engineers with a reference basis for reinforcement decision of old concrete dams.

3D Finite Element Analysis of TBM Water Diversion Tunnel Segment Coupled with Seepage Field

In most studies of tunnel boring machine(TBM)tunnelling, the groundwater pressure was not considered, or was simplified and exerted on the boundary of lining structure. Meanwhile, the leakage, which mainly occurs in the segment joints, was often ignored in the relevant studies of TBM tunnelling. Additionally, the geological models in these studies were simplified to different extents, and mostly were simplified as homogenous bodies. Considering the deficiencies above, a 3D refined model of the surrounding rock of a tunnel is firstly established using NURBS-TIN-BRe P hybrid data structure in this paper. Then the seepage field of the surrounding rock considering the leakage in the segment joints is simulated. Finally, the stability of TBM water diversion tunnel is studied coupled with the seepage simulation, to analyze the stress-strain conditions, the axial force and the bending moment of tunnel segment considering the leakage in the segment joints. The results illustrate that the maximum radial displacement, the minimum principal stress, the maximum principal stress and the axial force of segment lining considering the seepage effect are all larger than those disregarding the seepage effect.

ATheoretical and Applied Study of Seepage under Coupling Between Seepage Field and Stress Field

In civil engineering, more and more geological hazards are due to ignoring the interaction between seepage field and stress field(such as the water gushing in tunnel and other underground engineering). Faced this problem, the article has given a mathematical model on coupling between seepage field and stress field, and carried out numerical simulation with FEM (finite element method). Finally, the numerical simulation of coupling between fractured groundwater seepage field and fractured water bearing media stress field on the longest tunnel in China shows that this method is successful. At the same time, the prediction of water gushing yield in this tunnels construction is given.

DYNAMICS OF GAS SEEPAGE AND ITS APPLICATIONS

Dynamics of gas seepage as a borderline subject of geosciences mainly studies the flow and distribution of gas in coalseams or gas-bearing strata. In this paper new dynamic models for coal gas flow are developed.Using in-situ measured parameters of coal gas dynamics, the new models are tested with three existing dynamic models in the world. The results show that the new models approach the reality more cIose than the other three models.In addition, the other relations or indices helped to evaluate gas flow in coalseam are proposed.

Establishment of Unstable Flow Model and Well Testing Analysis for Viscoelastic Polymer Flooding

At present, the polymer solution is usually assumed to be Newtonian fluid or pseudoplastic fluid, and its elasticity is not considered on the study of polymer flooding well testing model. A large number of experiments have shown that polymer solutions have viscoelasticity, and disregarding the elasticity will cause certain errors in the analysis of polymer solution seepage law. Based on the percolation theory, this paper describes the polymer flooding mechanism from the two aspects of viscous effect and elastic effect, the mathematical model of oil water two-phase three components unsteady flow in viscoelastic polymer flooding was established, and solved by finite difference method, and the well-test curve was drawn to analyze the rule of well test curve in polymer flooding. The results show that, the degree of upward warping in the radial flow section of the pressure recovery curve when considering polymer elasticity is greater than the curve which not considering polymer elasticity. The relaxation time, power-law index, polymer injection concentration mainly affect the radial flow stage of the well testing curve. The relaxation time, power-law index, polymer injection concentration and other polymer flooding parameters mainly affect the radial flow stage of the well testing curve. The larger the polymer flooding parameters, the greater the degree of upwarping of the radial flow derivative curve. This model has important reference significance for well-testing research in polymer flooding oilfields.