Simulation of Underground Reservoir Stability of Pumped Storage Power Station Based on Fluid-Structure Coupling

Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its construction environment is more complex than that of a traditional reservoir.In particular,the stability of the rock strata in the underground reservoirs is affected by the seepage pressure and rock stress,which presents some challenges in achieving engineering safety and stability.Using the advantages of the numerical simulation method in dealing deal with nonlinear problems in engineering stability,in this study,the stability of the underground reservoir of the Shidangshan(SDS)pumped storage power station was numerically calculated and quantitatively analyzed based on fluid-structure coupling theory,providing an important reference for the safe operation and management of the underground reservoir.First,using the COMSOL software,a suitablemechanicalmodel was created in accordance with the geological structure and project characteristics of the underground reservoir.Next,the characteristics of the stress field,displacement field,and seepage field after excavation of the underground reservoir were simulated in light of the seepage effect of groundwater on the nearby rock of the underground reservoir.Finally,based on the construction specifications and Molar-Coulomb criterion,a thorough evaluation of the stability of the underground reservoir was performed through simulation of the filling and discharge conditions and anti-seepage strengthening measures.The findings demonstrate that the numerical simulation results have a certain level of reliability and are in accordance with the stress measured in the project area.The underground reservoir excavation resulted in a maximum displacement value of the rock mass around the caverns of 3.56 mm in a typical section,and the safety coefficient of the parts,as determined using the Molar-Coulomb criterion,was higher than 1,indicating that the project as a whole is in a stable state.

Role of land-atmosphere coupling in persistent extreme climate events in eastern China in summer 2022

2022年暖季,中国东部地区遭受持续性高温,少雨和土壤干旱的复合极端事件.特征分析指出,在研究时段内,中国东部地区的气温,降水和土壤湿度呈现明显的季节内变化和南北差异。由1940-2022年的气候态可知,长江流域和东南地区的土壤含水充足,蒸散主要受限于陆面有效能量.然而,潜在机制研究指出,2022年土壤湿度对蒸散的限制作用在上述区域异常偏强,土壤湿度与气候要素之间的强反馈可能在2022年复合极端事件的演变和持续中发挥了关键作用。

Implicit scheme for integrating constitutive model of unsaturated soils with coupling hydraulic and mechanical behavior

A constitutive model of unsaturated soils with coupling capillary hystere- sis and skeleton deformation is developed and implemented in a fully coupled transient hydro-mechanical finite-element model （computer code U-DYSAC2）. The obtained re- sults are compared with experimental results, showing that the proposed constitutive model can simulate the main mechanical and hydraulic behavior of unsaturated soils in a unified framework. The non-lineaxity of the soil-water characteristic relation is treated in a similar way of elastoplasticity. Two constitutive relations axe integrated by an implicit return-mapping scheme similar to that developed for saturated soils. A consistent tan- gential modulus is derived to preserve the asymptotic rate of the quadratic convergence of Newton＇s iteration. Combined with the integration of the constitutive model, a complete finite-element formulation of coupling hydro-mechanical problems for unsaturated soils is presented. A number of practical problems with different given initial and boundary conditions are analyzed to illustrate the performance and capabilities of the finite-element model.

Fully fluid-solid coupling dynamic model for seismic response of underground structures in saturated soils

The seismic response characteristics of underground structures in saturated soils are investigated.A fully fluid-solid coupling dynamic model is developed and implemented into ABAQUS with a user-defined element to simulate the dynamic behavior of saturated soils.The accuracy of the model is validated using a classic example in literature.The performance of the model is verified by its application on simulating the seismic response characteristics of a subway station built in saturated soils.The merits of the model are demonstrated by comparing the difference of the seismic response of an underground structure in saturated soils between using the fully coupling model and a single-phase medium model.The study finds that the fully coupling model developed herein can simulate the dynamic response characteristics of the underground structures in saturated soils with high accuracy.The seismic response of the underground structure tends to be underestimated by using the single-phase medium model compared with using the fully coupling model,which provides a weaker confining action to the underground structure.

Numerical simulation and experimental validation of moisture-heat coupling for saturated frozen soils

In seasonally frozen regions,freezing-and-thawing action is the main cause responsible for the destruction of canals,which is closely linked to the temperature gradient and water transport.To investigate the behaviour of soils under freezing-and-thawing actions,many numerical models have been established that consider the important coupling of moisture transport and temperature evolution;but they contain excessive parameters,some of which are rather difficult to determine.Based on the well-known Harlan's theory,a simple moisture-heat coupling model was recently proposed to quantify the coupled moisture-heat transport performance of soils in terms of the central temperature and porosity.The mathematical module of COMSOL Multiphysics was further employed to solve the governing equations numerically.To validate our model,a thorough experimental scheme was carried out in our lab.The measured temperature distribution was found to be consistent with the predicted results.

Formulation of thermo-hydro-mechanical coupling behavior of unsaturated soils based on hybrid mixture theory

Thermo-Hydro-Mechanical （THM） coupling pro- cesses in unsaturated soils are very important in both theoretical researches and engineering applications. A coupled formulation based on hybrid mixture theory is derived to model the THM coupling behavior of unsaturated soils. The free-energy and dissipative functions for different phases are derived from Taylor＇s series expansions. Constitutive relations for THM coupled behaviors of unsaturated soils, which include deformation, entropy change, fluid flow, heat conduction, and dynamic compatibility conditions on the interfaces, are then established. The number of field equations is shown to be equal to the number of unknown variables; thus, a closure of this coupling problem is established. In addition to modifications of the physical conservation equations with coupling effect terms, the constitutive equations, which consider the coupling between elastoplastic deformation of the soil skeleton, fluid flow, and heat transfer, are also derived.

Analysis of Coupling between Soil and Water Conservation and Economic-social Development

The coupling relation exists in water and soil conser-vation and economic-social development. The article analyses the relation of soil and water conservation and economic-social development stages as well as the coupling analytical method. Then calculates the expecting income by dispersing Markov decision and calculates the correlation coefficient and the re-lationship degree. The article obtains the relationship of soil and water conservation investments and all kinds of incomes. Finally, it analyzes the important meaning in socio-economic development of water and soil conservation.

Coupling numerical simulation with remotely sensed information for the study of frozen soil dynamics

The acquisition of spatial-temporal information of frozen soil is fundamental for the study of frozen soil dynamics and its feedback to climate change in cold regions.With advancement of remote sensing and better understanding of frozen soil dynamics,discrimination of freeze and thaw status of surface soil based on passive microwave remote sensing and numerical simulation of frozen soil processes under water and heat transfer principles provides valuable means for regional and global frozen soil dynamic monitoring and systematic spatial-temporal responses to global change.However,as an important data source of frozen soil processes,remotely sensed information has not yet been fully utilized in the numerical simulation of frozen soil processes.Although great progress has been made in remote sensing and frozen soil physics,yet few frozen soil research has been done on the application of remotely sensed information in association with the numerical model for frozen soil process studies.In the present study,a distributed numerical model for frozen soil dynamic studies based on coupled water-heat transferring theory in association with remotely sensed frozen soil datasets was developed.In order to reduce the uncertainty of the simulation,the remotely sensed frozen soil information was used to monitor and modify relevant parameters in the process of model simulation.The remotely sensed information and numerically simulated spatial-temporal frozen soil processes were validated by in-situ field observations in cold regions near the town of Naqu on the East-Central Tibetan Plateau.The results suggest that the overall accuracy of the algorithm for discriminating freeze and thaw status of surface soil based on passive microwave remote sensing was more than 95%.These results provided an accurate initial freeze and thaw status of surface soil for coupling and calibrating the numerical model of this study.The numerically simulated frozen soil processes demonstrated good performance of the distributed numerical model based on the coupled water-heat transferring theory.The relatively larger uncertainties of the numerical model were found in alternating periods between freezing and thawing of surface soil.The average accuracy increased by about 5%after integrating remotely sensed information on the surface soil.The simulation accuracy was significantly improved,especially in transition periods between freezing and thawing of the surface soil.

Simulation of Moving Bed Erosion Based on the Weakly Compressible Smoothed Particle Hydrodynamics-Discrete Element Coupling Method

A complex interface exists between waterflow and solid particles during hydraulic soil erosion.In this study,the particle discrete element method(DEM)has been used to simulate the hydraulic erosion of a granular soil under moving bed conditions and surrounding terrain changes.Moreover,the weakly compressible smoothed particle hydrodynamics(WCSPH)approach has been exploited to simulate the instability process of the free surfacefluid and its propagation characteristics at the solid–liquid interface.The influence of a suspended medium on the waterflow dynamics has been characterized using the mixed viscosity concept accounting for the solid–liquid mixed particle volume ratio.Numerical simulations of wall-jet scouring and reservoir sedimentflushing on a mobile bed were performed and validated with experiments.The results show that the proposed WCSPH–DEM coupling model is highly suitable for determining parameters,such as the local maximum scour depth,the scour pit width,and the sand bed profile.The effects on the hydraulic erosion process of two important para-meters of the mixed viscosity coefficient(initial solid volume concentration and initial viscosity coefficient)are also discussed to a certain extent in this study.

Upper Bound Solution of Soil Slope Stability under Coupling Effect of Rainfall and Earthquake

Based on the limit analysis upper bound method,a new model of soil slope collapse has been proposed which consists of two rigid block zones and a plastic shear zone.Soil slope was induced failure by coupling effect of rainfall and earthquake,and these blocks were also incorporated horizontal earthquake force and vertical gravitate.The velocities and forces were analyzed in three blocks,and the expression of velocity discontinuities was obtained by the principle of incompressibility.The external force work for the blocks,the internal energy of the plastic shear zone and the velocity discontinuous were solved.The present stability ratios are compared to the prevenient research,which shows the superiority of the mechanism and rationality of the analysis.The critical height of the soil slope can provide theoretical basis for slope support and design.

Simplification and Experimental Verification for Temperature and Humidity Field Coupling Model of Conservatory Soil

Based on the surface energy balance model which is widely used abroad, a temperature and humidity field coupling model of conservatory soil without crop vegetation in full illumination was established. Considering the relatively closed environment in conservatory, weak solar radiation and little surface evaporation of soil, the daily variation of water content in different soil layers may be neglected, then the temperature and humidity field coupling model was simplified to a one-dimensional thermal diffusion model. The simplified model and the temperature and humidity field coupling model adopted the same computational method of soil physical parameters and discrete format of heat diffusion differential equations, and were applied to the continuous simulation of temperature field in conservatory soil without crop vegetation in full illumination. Through the comparison between simulation results and experimental data, the precision of the simplified model was verified. The typical rule of soil heat flux variation in a 24 h cycle was also obtained.

WAVE-SOIL-PIPE COUPLING EFFECT ON SUBMARINE PIPELINE ON-BOTTOM STABILITY

Wave-soil-pipe coupling effect on the untrenched pipeline stability on sands is for the first time investigated experimentally. Tests are conducted in the U-shaped water tunnel, which generates an oscillatory how, simulating the water particle movements with periodically changing direction under the wave action. Characteristic times and phases during the instability process are revealed. Linear relationship between Froude number and non-dimensional pipe weight is obtained. Effects of initial embedment and loading history are observed. Test results between the wavesoil-pipe interaction and pipe-soil interaction under cyclic mechanical loading are compared. The mechanism is briefly discussed. For applying in the practical design, more extensive and systematic investigations are needed.

Stability analysis of unsaturated soil slope during rainfall infiltration using coupled liquid-gas-solid three-phase model

Generally, most soil slope failures are induced by rainfall infiltration, a process that involves interactions between the liquid phase, gas phase,and solid skeleton in an unsaturated soil slope. In this study, a loosely coupled liquid-gas-solid three-phase model, linking two numerical codes,TOUGH2/EOS3, which is used for water-air two-phase flow analysis, and FLAC^(3D), which is used for mechanical analysis, was established. The model was validated through a documented water drainage experiment over a sandy column and a comparison of the results with measured data and simulated results from other researchers. The proposed model was used to investigate the features of water-air two-phase flow and stress fields in an unsaturated soil slope during rainfall infiltration. The slope stability analysis was then performed based on the simulated water-air two-phase seepage and stress fields on a given slip surface. The results show that the safety factor for the given slip surface decreases first, then increases, and later decreases until the rainfall stops. Subsequently, a sudden rise occurs. After that, the safety factor decreases continually and reaches its lowest value, and then increases slowly to a steady value. The lowest value does not occur when the rainfall stops, indicating a delayed effect of the safety factor. The variations of the safety factor for the given slip surface are therefore caused by a combination of pore-air pressure, matric suction, normal stress, and net normal stress.

Regional Features and Seasonality of Land–Atmosphere Coupling over Eastern China

Land-atmosphere coupling is a key process of the climate system, and various coupling mechanisms have been proposed before based on observational and numerical analyses. The impact of soil moisture（SM） on evapotranspiration（ET） and further surface temperature（ST） is an important aspect of such coupling. Using ERA-Interim data and CLM4.0 offline simulation results, this study further explores the relationships between SM/ST and ET to better understand the complex nature of the land-atmosphere coupling（i.e., spatial and seasonal variations） in eastern China, a typical monsoon area. It is found that two diagnostics of land-atmosphere coupling（i.e., SM-ET correlation and ST-ET correlation） are highly dependent on the climatology of SM and ST. By combining the SM-ET and ST-ET relationships, two ＂hot spots＂ of land-atmosphere coupling over eastern China are identified： Southwest China and North China. In Southwest China, ST is relatively high throughout the year, but SM is lowest in spring, resulting in a strong coupling in spring. However, in North China, SM is relatively low throughout the year, but ST is highest in summer, which leads to the strongest coupling in summer. Our results emphasize the dependence of land-atmosphere coupling on the seasonal evolution of climatic conditions and have implications for future studies related to land surface feedbacks.

FEM-DEM coupled modeling of cone penetration tests in lunar soil

Cone penetration test(CPT)is an appropriate technique for quickly determining the geotechnical properties of lunar soil,which is valuable for in situ lunar exploration.Utilizing a typical coupling method recently developed by the authors,a finite element method(FEM)-discrete element method(DEM)coupled model of CPTs is obtained.A series of CPTs in lunar soil are simulated to qualitatively reveal the flow of particles and the development of resistance throughout the penetration process.In addition,the effects of major factors,such as penetration velocity,penetration depth,cone tip angle,and the low gravity on the Moon surface are investigated.

Coupled Transfer of Water and Heat in Red Soil: Experiment and Numerical Modelling

Coupled transfer of soil water and heat in closed columns of homogeneous red soil was studied under laboratory conditions. A coupled model was constructed using soil physical theory, empirical equations and experimental data to predict the coupled transfer. The results show that transport of soil water was affected by temperature gradient, and the largest net water transport was found in the soil column with initial water content of 0.148 m3 m-3. At the same time, temperature changes with the transport of soil water was in a nonlinear shape as heat parameters were function of water content, and the changes of temperature were positively correlated with the net amount of water transported. Numerical modelling results show that the predicted values of temperature distribution were close to the observed values, while the predicted values of water content exhibited limited deviation at both ends of the soil column due to the slight temperature changes at both ends. It was indicated that the model proposed here was applicable.

Different discretization method used in coupled water and heat transport mode for soil under freezing conditions

A coupled water and heat transport mode is established based on the Richards equation to study water flow and heat transport in soil during freezing process. Both the finite difference and finite element method are used in the discretization, respectively. Two different computer programs are written and used to simulate an indoor unidirectional frozen test. The freezing depth, freezing rate and temperature variation are compared among lab tests, finite difference calculation simulation and finite element calculation simulation. Result shows that： the finite difference method has a better performance in freezing depth simulation while the finite element method has a better performance in numerical stability in one-dimensional freezing simulation.

Hydro-mechanical coupling mechanism on joint of clay core-wall and concrete cut-off wall

The joint of clay core-wall and concrete cut-off wall is one of the weakest parts in high earth and rockftll dams.A kind of highly plastic clay is always fixed on the joint to fit the large shear deformation between clay core-wall and concrete cut-offwall,so the hydro-mechanical coupling mechanisms on the joint under high stress,high hydraulic gradient,and large shear deformation are of great importance for the evaluation of dam safety.The hydro-mechanical coupling characteristics of the joint of the highly plastic clay and the concrete cut-off wall in a high earth and rockfill dam in China were studied by using a newly designed soil-structure contact erosion apparatus.The experimental results indicate that：1） Shear failure on the joint is due to the hydro-mechanical coupling effect of stress and seepage failure.The seepage failure will induce the final shear failure when the ratio of deviatoric stress to confining pressure is within 1.0-1.2; 2） A negative exponential permeability empirical model for the joint denoted by a newly defined principal stress function,which considers the coupling effect of confining pressure and axial pressure on the permeability,is established based on hydro-mechanical coupling experiments.3） The variation of the settlement before and after seepage failure is very different.The settlement before seepage failure changes very slowly,while it increases significantly after the seepage failure.4） The stress-strain relationship is of a strain softening type.5） Flow along the joint still follows Darcian flow rule.The results will provide an important theoretical basis for the further evaluation on the safety of the high earth and rockfill dam.

Assimilation of ASAR Data with a Hydrologic and Semi-empirical Backscattering Coupled Model to Estimate Soil Moisture

The most promising approach for studying soil moisture is the assimilation of observation data and computational modeling. However, there is much uncertainty in the assimilation process, which affects the assimilation results. This research developed a one-dimensional soil moisture assimilation scheme based on the Ensemble Kalman Filter (EnKF) and Genetic Algorithm (GA). A two-dimensional hydrologic model-Distributed Hydrology-Soil-Vegetation Model (DHSVM) was coupled with a semi-empirical backscattering model (Oh). The Advanced Synthetic Aperture Radar (ASAR) data were assimilated with this coupled model and the field observation data were used to validate this scheme in the soil moisture assimilation experiment. In order to improve the assimilation results, a cost function was set up based on the distance between the simulated backscattering coefficient from the coupled model and the observed backscattering coefficient from ASAR. The EnKF and GA were used to re-initialize and re-parameterize the simulation process, respectively. The assimilation results were compared with the free-run simulations from hydrologic model and the field observation data. The results obtained indicate that this assimilation scheme is practical and it can improve the accuracy of soil moisture estimation significantly.

Factors controlling deep-profile soil organic carbon and water storage following Robinia pseudoacacia afforestation of the Loess Plateau in China

Background:Afforestation is a common and effective approach used for the restoration of degraded ecosystems worldwide.In China,Robinia pseudoacacia(RP)is among the main non-native tree species and has been widely planted in revegetation of the Loess Plateau.However,owing to uncertainties regarding soil water consumption and carbon sequestration,it is necessary to assess the suitability and sustainability of R.pseudoacacia in restoration.In this study,we aimed to analyse the dynamic effects of R.pseudoacacia forest on soil carbon storage(SCS)and soil water storage(SWS).Specifically,we investigated the association between soil water content(SWC)and soil organic carbon(SOC)and underlying factors in the 0-500-cm profile of a 10-to 50-year-old chronosequence.Results:The results obtained indicated that the dynamics of SWS and SCS on this time scale could be divided into an initial reduction phase(the initial 20 years after afforestation)and subsequent recovery(20-50 years after afforestation).Compared with in the abandoned land(AL),the net accumulation of SCS in R.pseudoacacia forest was 2.51 Mg·ha^(-1)at 50 years after afforestation,whereas there was a 398.76-mm deficit in SWS.Additionally,the natural succession of R.pseudoacacia forest has contributed to the continuous change in stand structure(e.g.vegetation cover(VC),understory vegetation coverage(UVC),and litter biomass(LB)).Conclusions:These findings indicate that vegetation restoration increases carbon sequestration while causing soil water deficit.Furthermore,stand density(SD)was established to make a predominant contribution to the dynamics of SWS and SCS via its effects in altering vegetation,soil,and litter characteristics.Therefore,high-density plantation forests in the Loess Plateau area should be appropriately thinned to reduce the density of forest stands on the basis of soil erosion control and wind and sand fixation,so as to increase carbon sink with lower water consumption,thus realizing the synergistic development of soil carbon sequestration and water connotation.