Accurate and reliable river flow information is critical to planning and management for sustainable water resources utilization. Most of engineering activities related to hydrologic designs, flood, drought, reservoirs...Accurate and reliable river flow information is critical to planning and management for sustainable water resources utilization. Most of engineering activities related to hydrologic designs, flood, drought, reservoirs and their operations are heavily dependent on the river flow information derived from river rating curve. The rating curve for a given river section is normally developed from a set of direct stage-discharge measurements for different periods. This involves considerable labour, risk and resources, and presupposes a complex and extensive measuring survey. Extrapolating the rating curve beyond the measured range, as common in many cases, is fraught with errors and uncertainties, due to the complex hydraulic behaviour of the surface water profile in transition from section, channel, downstream and flood plain controls which are often poorly understood with direct measurements. Hydraulic modeling has recently emerged as one of the more promising methods to efficiently develop accurate rating curves for a river section with simple or complex hydraulic structures and conditions. This paper explores the use of a Hydraulic Engineering Center-River Analysis System (HEC-RAS) model to review and develop river rating curves for three hydrometric stations on two rivers in Kwale, coastal Kenya. The HEC-RAS models were set up based on topographical (cross section and longitudinal) survey data for the reaches and engineering drawings for the hydraulic structures commonly used as section controls for flow measurement. The model was calibrated under unsteady state conditions against measured stage-discharge data which were captured using a Velocity Current Meter (Valeport) and an Acoustic Doppler Current Profiler (ADCP) for both low and high flow. The rating curves were extracted from model results and the uncertainty associated with each rating curve analyzed. The results obtained by the HEC-RAS model were satisfactory and deemed acceptable for predicting discharge across the stage range at each river section.展开更多
This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fr...This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.展开更多
Loess has distinctive characteristics,leading to frequent landslide disasters and posing serious threats to the lives and properties of local re sidents.The involvement of water repre sents a critical factor in induci...Loess has distinctive characteristics,leading to frequent landslide disasters and posing serious threats to the lives and properties of local re sidents.The involvement of water repre sents a critical factor in inducing loess landslides.This study focuses on three neighboring cities sequentially situated on the Loess Plateau along the direction of aeolian deposition of loess,namely Lanzhou,Dingxi,and Tianshui,which are densely populated and prone to landslide disasters.The variations in hydraulic properties,including water retention capacity and permeability,are investigated through Soil Water Characteristic Curve(SWCC)test and hydraulic conductivity test.The experimental findings revealed that Tianshui loess exhibited the highest water retention capacity,followed by Dingxi loess,while Lanzhou loess demonstrated the lowest water retention capacity.Contrastingly,the results for the saturated permeability coefficient were found to be the opposite:Tianshui loess showed the lowest permeability,whereas Lanzhou loess displayed the highest permeability.These results are supported and analyzed by scanning electron microscopy(SEM)observation.In addition,the water retention capacity is mathematically expressed using the van Genuchten model and extended to predict unsaturated hydraulic properties of loess.The experimental results exhibit a strong accordance with one another and align with the regional distribution patterns of disasters.展开更多
Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operation...Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.展开更多
The eastern foothills of the Helan Mountains in China are a typical mountainous region of soil and gravel,where gravel could affect the water movement process in the soil.This study focused on the effects of different...The eastern foothills of the Helan Mountains in China are a typical mountainous region of soil and gravel,where gravel could affect the water movement process in the soil.This study focused on the effects of different gravel contents on the water absorption characteristics and hydraulic parameters of stony soil.The stony soil samples were collected from the eastern foothills of the Helan Mountains in April 2023 and used as the experimental materials to conduct a one-dimensional horizontal soil column absorption experiment.Six experimental groups with gravel contents of 0%,10%,20%,30%,40%,and 50%were established to determine the saturated hydraulic conductivity(K_(s)),saturated water content(θ_(s)),initial water content(θ_(i)),and retention water content(θ_(r)),and explore the changes in the wetting front depth and cumulative absorption volume during the absorption experiment.The Philip model was used to fit the soil absorption process and determine the soil water absorption rate.Then the length of the characteristic wetting front depth,shape coefficient,empirical parameter,inverse intake suction and soil water suction were derived from the van Genuchten model.Finally,the hydraulic parameters mentioned above were used to fit the soil water characteristic curves,unsaturated hydraulic conductivity(K_(θ))and specific water capacity(C(h)).The results showed that the wetting front depth and cumulative absorption volume of each treatment gradually decreased with increasing gravel content.Compared with control check treatment with gravel content of 0%,soil water absorption rates in the treatments with gravel contents of 10%,20%,30%,40%,and 50%decreased by 11.47%,17.97%,25.24%,29.83%,and 42.45%,respectively.As the gravel content increased,inverse intake suction gradually increased,and shape coefficient,K_(s),θ_(s),andθ_(r)gradually decreased.For the same soil water content,soil water suction and K_(θ)gradually decreased with increasing gravel content.At the same soil water suction,C(h)decreased with increasing gravel content,and the water use efficiency worsened.Overall,the water holding capacity,hydraulic conductivity,and water use efficiency of stony soil in the eastern foothills of the Helan Mountains decreased with increasing gravel content.This study could provide data support for improving soil water use efficiency in the eastern foothills of the Helan Mountains and other similar rocky mountainous areas.展开更多
In conjunction with the working characteristics of the high-clearance wheeled sprayer and the benefits of the closed hydraulic system,a series of reasonable working parameters should be established,and a hydraulic sys...In conjunction with the working characteristics of the high-clearance wheeled sprayer and the benefits of the closed hydraulic system,a series of reasonable working parameters should be established,and a hydraulic system that fulfills the requisite specifications should be designed.The AMESim software model is employed to construct a closed hydraulic transmission system,and the simulation analysis is then performed according to the data of hydraulic components.According to analysis results,the prototype can be optimized and upgraded,and a verification test is further carried out.The test results demonstrate that the designed closed hydraulic transmission system meets the actual working requirements of the high-clearance wheeled sprayer and provides a stable experimental platform for intelligent control of agricultural machinery.展开更多
A jet mill bit(JMB)is proposed to increase the drilling efficiency and safety of horizontal wells,which has the hydraulic characteristics of depressurization and cuttings cleaning.This paper fills the gap in the hydra...A jet mill bit(JMB)is proposed to increase the drilling efficiency and safety of horizontal wells,which has the hydraulic characteristics of depressurization and cuttings cleaning.This paper fills the gap in the hydraulic study of the JMB by focusing on the hydraulic modeling and optimization of the JMB and considering these two hydraulic characteristics.First,the hydraulic depressurization model and the hydraulic cuttings cleaning model of the JMB are developed respectively.In the models,the pressure ratio and efficiency are chosen as the evaluation parameters of the depressurization capacity of the JMB,and the jet hydraulic power and jet impact force are chosen as the evaluation parameters of cuttings cleaning capacity of the JMB.Second,based on the hydraulic models,the effects of model parameters[friction loss coefficient,target inclination angle,rate of penetration(ROP),flow ratio,and well depth]on the hydraulic performance of the JMB are investigated.The results show that an increase in the friction loss coefficient and target inclination angle cause a significant reduction in the hydraulic depressurization capacity,and the effect of ROP is negligible.The flow ratio is positively related to the hydraulic cuttings cleaning capacity,and the well depth determines the maximum hydraulic cuttings cleaning capacity.Finally,by combining the hydraulic depressurization model and hydraulic cuttings cleaning model,an optimization method of JMB hydraulics is proposed to simultaneously maximize the jet depressurization capacity and the cuttings cleaning capacity.According to the drilling parameters given,the optimal values of the drilling fluid flow rate,backward nozzle diameter,forward nozzle diameter,and throat diameter can be determined.Moreover,a case study is conducted to verify the effectiveness of the optimization method.展开更多
Unsaturated expansive soil is widely distributed in China and has complex engineering properties.This paper proposes the unified hydraulic effect shear strength theory of unsaturated expansive soil based on the effect...Unsaturated expansive soil is widely distributed in China and has complex engineering properties.This paper proposes the unified hydraulic effect shear strength theory of unsaturated expansive soil based on the effective stress principle,swelling force principle,and soil–water characteristics.Considering the viscoelasticity and structural damage of unsaturated expansive soil during loading,a fractional hardening–damage model of unsaturated expansive soil was established.The model parameters were established on the basis of the proposed calculation method of shear strength and the triaxial shear experiment on unsaturated expansive soil.The proposed model was verified by the experimental data and a traditional damage model.The proposed model can satisfactorily describe the entire process of the strain-hardening law of unsaturated expansive soil.Finally,by investigating the damage variables of the proposed model,it was found that:(a)when the values of confining pressure and matric suction are close,the coupling of confining pressure and matric suction contributes more to the shear strength;(b)there is a damage threshold for unsaturated expansive soil,and is mainly reflected by strength criterion of infinitesimal body;(c)the strain hardening law of unsaturated expansive soil is mainly reflected by fractional derivative operator.展开更多
Hydropower gains increasing importance as a steerable and controllable power source in a renewable energy mix and deregulated markets. Although hydropower produces fossil-free energy, it has a significant impact on th...Hydropower gains increasing importance as a steerable and controllable power source in a renewable energy mix and deregulated markets. Although hydropower produces fossil-free energy, it has a significant impact on the local environment. This review investigates the effects of flow alterations by hydropower on the downstream river system and the possibilities to integrate these effects into hydraulic modeling. The results show that various effects of flow regulation on the ecosystem, but also social and economic effects on related communities were observed in the last decades. The application of hydraulic models for investigations of ecological effects is common. Especially hydraulic effects and effects on fish were extensively modeled with the help of hydraulic 1D- and 2D-simulations. Current applications to investigate social and economic effects integrated into hydraulic modeling are meanwhile limited. Approaches to realizing this integration are presented. Further research on the economic valuation of ecosystems and integration of social and economic effects to hydraulic models is necessary to develop holistic tools to support decision-making on sustainable hydropower.展开更多
Hydraulic servo system plays an important role in industrial fields due to the advantages of high response,small size-to-power ratio and large driving force.However,inherent nonlinear behaviors and modeling uncertaint...Hydraulic servo system plays an important role in industrial fields due to the advantages of high response,small size-to-power ratio and large driving force.However,inherent nonlinear behaviors and modeling uncertainties are the main obstacles for hydraulic servo system to achieve high tracking perfor-mance.To deal with these difficulties,this paper presents a backstepping sliding mode controller to improve the dynamic tracking performance and anti-interfer-ence ability.For this purpose,the nonlinear dynamic model is firstly established,where the nonlinear behaviors and modeling uncertainties are lumped as one term.Then,the extended state observer is introduced to estimate the lumped distur-bance.The system stability is proved by using the Lyapunov stability theorem.Finally,comparative simulation and experimental are conducted on a hydraulic servo system platform to verify the efficiency of the proposed control scheme.展开更多
Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient applicat...Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.展开更多
The previous sensitivity analysis researches are not accurate enough and also have the limited reference value, because those mathematical models are relatively simple and the change of the load and the initial displa...The previous sensitivity analysis researches are not accurate enough and also have the limited reference value, because those mathematical models are relatively simple and the change of the load and the initial displacement changes of the piston are ignored, even experiment verification is not conducted. Therefore, in view of deficiencies above, a nonlinear mathematical model is established in this paper, including dynamic characteristics of servo valve, nonlinear characteristics of pressure-flow, initial displacement of servo cylinder piston and friction nonlinearity. The transfer function block diagram is built for the hydraulic drive unit closed loop position control, as well as the state equations. Through deriving the time-varying coefficient items matrix and time-varying free items matrix of sensitivity equations respectively, the expression of sensitivity equations based on the nonlinear mathematical model are obtained. According to structure parameters of hydraulic drive unit, working parameters, fluid transmission characteristics and measured friction-velocity curves, the simulation analysis of hydraulic drive unit is completed on the MATLAB/Simulink simulation platform with the displacement step 2 mm, 5 mm and 10 mm, respectively. The simulation results indicate that the developed nonlinear mathematical model is sufficient by comparing the characteristic curves of experimental step response and simulation step response under different constant load. Then, the sensitivity function time-history curves of seventeen parameters are obtained, basing on each state vector time-history curve of step response characteristic. The maximum value of displacement variation percentage and the sum of displacement variation absolute values in the sampling time are both taken as sensitivity indexes. The sensitivity indexes values above are calculated and shown visually in histograms under different working conditions, and change rules are analyzed. Then the sensitivity indexes values of four measurable parameters, such as supply pressure, proportional gain, initial position of servo cylinder piston and load force, are verified experimentally on test platform of hydraulic drive unit, and the experimental research shows that the sensitivity analysis results obtained through simulation are approximate to the test results. This research indicates each parameter sensitivity characteristics of hydraulic drive unit, the performance-affected main parameters and secondary parameters are got under different working conditions, which will provide the theoretical foundation for the control compensation and structure optimization of hydraulic drive unit.展开更多
This paper presents an improved understanding of coupled hydro-thermo-mechanical(HTM) hydraulic fracturing of quasi-brittle rock using the bonded particle model(BPM) within the discrete element method(DEM). BPM has be...This paper presents an improved understanding of coupled hydro-thermo-mechanical(HTM) hydraulic fracturing of quasi-brittle rock using the bonded particle model(BPM) within the discrete element method(DEM). BPM has been recently extended by the authors to account for coupled convective econductive heat flow and transport, and to enable full hydro-thermal fluidesolid coupled modeling.The application of the work is on enhanced geothermal systems(EGSs), and hydraulic fracturing of hot dry rock(HDR) is studied in terms of the impact of temperature difference between rock and a flowing fracturing fluid. Micro-mechanical investigation of temperature and fracturing fluid effects on hydraulic fracturing damage in rocks is presented. It was found that fracture is shorter with pronounced secondary microcracking along the main fracture for the case when the convectiveeconductive thermal heat exchange is considered. First, the convection heat exchange during low-viscosity fluid infiltration in permeable rock around the wellbore causes significant rock cooling, where a finger-like fluid infiltration was observed. Second, fluid infiltration inhibits pressure rise during pumping and delays fracture initiation and propagation. Additionally, thermal damage occurs in the whole area around the wellbore due to rock cooling and cold fluid infiltration. The size of a damaged area around the wellbore increases with decreasing fluid dynamic viscosity. Fluid and rock compressibility ratio was found to have significant effect on the fracture propagation velocity.展开更多
Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presen...Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.展开更多
Instead of establishing mathematical hydraulic system models from physical laws usually done with the problems of complex modelling processes, low reliability and practicality caused by large uncertainties, a novel mo...Instead of establishing mathematical hydraulic system models from physical laws usually done with the problems of complex modelling processes, low reliability and practicality caused by large uncertainties, a novel modelling method for a highly nonlinear system of a hydraulic excavator is presented. Based on the data collected in the excavator's arms driving experiments, a data-based excavator dynamic model using Simplified Refined Instrumental Variable (SRIV) identification and estimation algorithms is established. The validity of the proposed data-based model is indirectly demonstrated by the performance of computer simulation and the.real machine motion control exoeriments.展开更多
Modeling and digital simulation is an effective method to analyze the dynamiccharacteristics of hydraulic system. It is difficult to determine some performance parameters inthe hydraulic system by means of currently u...Modeling and digital simulation is an effective method to analyze the dynamiccharacteristics of hydraulic system. It is difficult to determine some performance parameters inthe hydraulic system by means of currently used modeling methods. The 'gray-box' modeling method forlarge-scale hydraulic system is introduced. The principle of the method, the submodels of somecomponents and the parameters identification of components or subsystem are researched.展开更多
Water resources management usually requires that hydraulic, ecological, and hydrological models be linked. The Hy- drologic Engineering Center River Analysis System (HEC-RAS) hydraulic model and the Hydrologic Enginee...Water resources management usually requires that hydraulic, ecological, and hydrological models be linked. The Hy- drologic Engineering Center River Analysis System (HEC-RAS) hydraulic model and the Hydrologic Engineering Center Geospatial River Analysis System (HEC-GEORAS), imitates flow and water profiles in the Neka river basin’s downstream flood plain. Hydrograph phases studied during the flood seasons of 1986-1999 and from 2002-2004 were used to calibrate and verify the hydraulic model respectively. Simulations of peak flood stages and hydrographs’ evaluations are congruent with studies and observations, with the former showing mean square errors between 4.8 - 10 cm. HECRAS calculations and forecast flood water levels. Nash-Sutcliffe effectiveness (CR3) is more than 0.92 along with elevated levels of water which were created with some effectiveness (CR5) of 0.94 for the validation period. The coupled two models show good performance in the water level modeling.展开更多
This paper presents a three-dimensional fully hydro-mechanical coupled distinct element study on fault reactivation and induced seismicity due to hydraulic fracturing injection and subsequent backflow process,based on...This paper presents a three-dimensional fully hydro-mechanical coupled distinct element study on fault reactivation and induced seismicity due to hydraulic fracturing injection and subsequent backflow process,based on the geological data in Horn River Basin,Northeast British Columbia,Canada.The modeling results indicate that the maximum magnitude of seismic events appears at the fracturing stage.The increment of fluid volume in the fault determines the cumulative moment and maximum fault slippage,both of which are essentially proportional to the fluid volume.After backflow starts,the fluid near the joint intersection keeps flowing into the critically stressed fault,rather than backflows to the wellbore.Although fault slippage is affected by the changes of both pore pressure and ambient rock stress,their contributions are different at fracturing and backflow stages.At fracturing stage,pore pressure change shows a dominant effect on induced fault slippage.While at backflow stage,because the fault plane is under a critical stress state,any minor disturbance would trigger a fault slippage.The energy analysis indicates that aseismic deformation takes up a majority of the total deformation energy during hydraulic fracturing.A common regularity is found in both fracturing-and backflow-induced seismicity that the cumulative moment and maximum fault slippage are nearly proportional to the injected fluid volume.This study shows some novel insights into interpreting fracturing-and backflowinduced seismicity,and provides useful information for controlling and mitigating seismic hazards due to hydraulic fracturing.展开更多
Modeling unsaturated flow in fractured rocks is essential in various subsurface engineering applications,but it remains a great challenge due to the difficulties in determining the unsaturated hydraulic properties of ...Modeling unsaturated flow in fractured rocks is essential in various subsurface engineering applications,but it remains a great challenge due to the difficulties in determining the unsaturated hydraulic properties of rocks that contain various scales of fractures.It is generally accepted that the van Genuchten(VG)model can be applied to fractured rocks,provided that the hydraulic parameters could be representatively determined.In this study,scaling relationships between the VG parameters(a and n)and hydraulic conductivity(K)across 8 orders of magnitude,from 10^(-10)m/s to 10^(-2)m/s,were proposed by statistical analysis of data obtained from 1416 soil samples.The correlations were then generalized to predict the upper bounds of VG parameters for fractured rocks from the K data that could be obtained more easily under field conditions,and were validated against a limited set of data from cores,fractures and fractured rocks available in the literature.The upper bound estimates significantly narrow the ranges of VG parameters,and the representative values of a and n for fractured rocks at the field scale can then be determined with confidence by inverse modeling using groundwater observations in saturated zones.The proposed methodology was applied to saturated-unsaturated flow modeling in the right-bank slope at the Baihetan dam site with a continuum approach,showing that most of the flow behaviors in fractured rocks in this complex hydrogeological condition could be properly reproduced.The proposed method overcomes difficulties in suction measurement in fractured rocks with strong heterogeneity,and provides a feasible way for modeling of saturated-unsaturated flow in fractured rocks with acceptable engineering accuracy.展开更多
文摘Accurate and reliable river flow information is critical to planning and management for sustainable water resources utilization. Most of engineering activities related to hydrologic designs, flood, drought, reservoirs and their operations are heavily dependent on the river flow information derived from river rating curve. The rating curve for a given river section is normally developed from a set of direct stage-discharge measurements for different periods. This involves considerable labour, risk and resources, and presupposes a complex and extensive measuring survey. Extrapolating the rating curve beyond the measured range, as common in many cases, is fraught with errors and uncertainties, due to the complex hydraulic behaviour of the surface water profile in transition from section, channel, downstream and flood plain controls which are often poorly understood with direct measurements. Hydraulic modeling has recently emerged as one of the more promising methods to efficiently develop accurate rating curves for a river section with simple or complex hydraulic structures and conditions. This paper explores the use of a Hydraulic Engineering Center-River Analysis System (HEC-RAS) model to review and develop river rating curves for three hydrometric stations on two rivers in Kwale, coastal Kenya. The HEC-RAS models were set up based on topographical (cross section and longitudinal) survey data for the reaches and engineering drawings for the hydraulic structures commonly used as section controls for flow measurement. The model was calibrated under unsteady state conditions against measured stage-discharge data which were captured using a Velocity Current Meter (Valeport) and an Acoustic Doppler Current Profiler (ADCP) for both low and high flow. The rating curves were extracted from model results and the uncertainty associated with each rating curve analyzed. The results obtained by the HEC-RAS model were satisfactory and deemed acceptable for predicting discharge across the stage range at each river section.
文摘This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.
基金the financial support for the research presented in this paper from National Natural Science Foundation of China(42201142,42067066,51778590)。
文摘Loess has distinctive characteristics,leading to frequent landslide disasters and posing serious threats to the lives and properties of local re sidents.The involvement of water repre sents a critical factor in inducing loess landslides.This study focuses on three neighboring cities sequentially situated on the Loess Plateau along the direction of aeolian deposition of loess,namely Lanzhou,Dingxi,and Tianshui,which are densely populated and prone to landslide disasters.The variations in hydraulic properties,including water retention capacity and permeability,are investigated through Soil Water Characteristic Curve(SWCC)test and hydraulic conductivity test.The experimental findings revealed that Tianshui loess exhibited the highest water retention capacity,followed by Dingxi loess,while Lanzhou loess demonstrated the lowest water retention capacity.Contrastingly,the results for the saturated permeability coefficient were found to be the opposite:Tianshui loess showed the lowest permeability,whereas Lanzhou loess displayed the highest permeability.These results are supported and analyzed by scanning electron microscopy(SEM)observation.In addition,the water retention capacity is mathematically expressed using the van Genuchten model and extended to predict unsaturated hydraulic properties of loess.The experimental results exhibit a strong accordance with one another and align with the regional distribution patterns of disasters.
基金funded by the joint fund of the National Key Research and Development Program of China(No.2021YFC2902101)National Natural Science Foundation of China(Grant No.52374084)+1 种基金Open Foundation of National Energy shale gas R&D(experiment) center(2022-KFKT-12)the 111 Project(B17009)。
文摘Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.
基金funded by the National Natural Science Foundation of China(32360321)the Natural Science Foundation of Ningxia Hui Autonomous Region,China(2023AAC03046,2023AAC02018)the Ningxia Key Research and Development Project(2021BEG02011).
文摘The eastern foothills of the Helan Mountains in China are a typical mountainous region of soil and gravel,where gravel could affect the water movement process in the soil.This study focused on the effects of different gravel contents on the water absorption characteristics and hydraulic parameters of stony soil.The stony soil samples were collected from the eastern foothills of the Helan Mountains in April 2023 and used as the experimental materials to conduct a one-dimensional horizontal soil column absorption experiment.Six experimental groups with gravel contents of 0%,10%,20%,30%,40%,and 50%were established to determine the saturated hydraulic conductivity(K_(s)),saturated water content(θ_(s)),initial water content(θ_(i)),and retention water content(θ_(r)),and explore the changes in the wetting front depth and cumulative absorption volume during the absorption experiment.The Philip model was used to fit the soil absorption process and determine the soil water absorption rate.Then the length of the characteristic wetting front depth,shape coefficient,empirical parameter,inverse intake suction and soil water suction were derived from the van Genuchten model.Finally,the hydraulic parameters mentioned above were used to fit the soil water characteristic curves,unsaturated hydraulic conductivity(K_(θ))and specific water capacity(C(h)).The results showed that the wetting front depth and cumulative absorption volume of each treatment gradually decreased with increasing gravel content.Compared with control check treatment with gravel content of 0%,soil water absorption rates in the treatments with gravel contents of 10%,20%,30%,40%,and 50%decreased by 11.47%,17.97%,25.24%,29.83%,and 42.45%,respectively.As the gravel content increased,inverse intake suction gradually increased,and shape coefficient,K_(s),θ_(s),andθ_(r)gradually decreased.For the same soil water content,soil water suction and K_(θ)gradually decreased with increasing gravel content.At the same soil water suction,C(h)decreased with increasing gravel content,and the water use efficiency worsened.Overall,the water holding capacity,hydraulic conductivity,and water use efficiency of stony soil in the eastern foothills of the Helan Mountains decreased with increasing gravel content.This study could provide data support for improving soil water use efficiency in the eastern foothills of the Helan Mountains and other similar rocky mountainous areas.
基金Supported by 2023 Xinjiang Uygur Autonomous Region R&D and Promotion and Application of Key Technologies of CNC Sprayer for Seed Corn(2023NC010).
文摘In conjunction with the working characteristics of the high-clearance wheeled sprayer and the benefits of the closed hydraulic system,a series of reasonable working parameters should be established,and a hydraulic system that fulfills the requisite specifications should be designed.The AMESim software model is employed to construct a closed hydraulic transmission system,and the simulation analysis is then performed according to the data of hydraulic components.According to analysis results,the prototype can be optimized and upgraded,and a verification test is further carried out.The test results demonstrate that the designed closed hydraulic transmission system meets the actual working requirements of the high-clearance wheeled sprayer and provides a stable experimental platform for intelligent control of agricultural machinery.
基金financially supported by Youth Project of Natural Science Basic Research Program of Shaanxi Province(Grant number:2023-JC-QN-0538)Scientifical Research Program for Youth Innovation Team Construction of Shaanxi Provincial Department of Education(Grant number:21JP054,22JP032)+1 种基金National Natural Science Foundation of China(Grant numbers:52174012,51804322,51821092,51774301,U1762214)Open Fund(PLC 20210404)of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Chengdu University of Technology)。
文摘A jet mill bit(JMB)is proposed to increase the drilling efficiency and safety of horizontal wells,which has the hydraulic characteristics of depressurization and cuttings cleaning.This paper fills the gap in the hydraulic study of the JMB by focusing on the hydraulic modeling and optimization of the JMB and considering these two hydraulic characteristics.First,the hydraulic depressurization model and the hydraulic cuttings cleaning model of the JMB are developed respectively.In the models,the pressure ratio and efficiency are chosen as the evaluation parameters of the depressurization capacity of the JMB,and the jet hydraulic power and jet impact force are chosen as the evaluation parameters of cuttings cleaning capacity of the JMB.Second,based on the hydraulic models,the effects of model parameters[friction loss coefficient,target inclination angle,rate of penetration(ROP),flow ratio,and well depth]on the hydraulic performance of the JMB are investigated.The results show that an increase in the friction loss coefficient and target inclination angle cause a significant reduction in the hydraulic depressurization capacity,and the effect of ROP is negligible.The flow ratio is positively related to the hydraulic cuttings cleaning capacity,and the well depth determines the maximum hydraulic cuttings cleaning capacity.Finally,by combining the hydraulic depressurization model and hydraulic cuttings cleaning model,an optimization method of JMB hydraulics is proposed to simultaneously maximize the jet depressurization capacity and the cuttings cleaning capacity.According to the drilling parameters given,the optimal values of the drilling fluid flow rate,backward nozzle diameter,forward nozzle diameter,and throat diameter can be determined.Moreover,a case study is conducted to verify the effectiveness of the optimization method.
基金financially supported by Sichuan Huaxi Group Co.,ltd.(No.HXKX2019/015,No.HXKX2019/019,No.HXKX2018/030)the Open Fund of Sichuan Provincial Engineering Research Center of City Solid Waste Energy and Building Materials Conversion and Utilization Technology(No.GF2022ZC009)the Open Fund of Sichuan Engineering Research Center for Mechanical Properties and Engineering Technology of Unsaturated Soils(No.SC-FBHT2022-04)。
文摘Unsaturated expansive soil is widely distributed in China and has complex engineering properties.This paper proposes the unified hydraulic effect shear strength theory of unsaturated expansive soil based on the effective stress principle,swelling force principle,and soil–water characteristics.Considering the viscoelasticity and structural damage of unsaturated expansive soil during loading,a fractional hardening–damage model of unsaturated expansive soil was established.The model parameters were established on the basis of the proposed calculation method of shear strength and the triaxial shear experiment on unsaturated expansive soil.The proposed model was verified by the experimental data and a traditional damage model.The proposed model can satisfactorily describe the entire process of the strain-hardening law of unsaturated expansive soil.Finally,by investigating the damage variables of the proposed model,it was found that:(a)when the values of confining pressure and matric suction are close,the coupling of confining pressure and matric suction contributes more to the shear strength;(b)there is a damage threshold for unsaturated expansive soil,and is mainly reflected by strength criterion of infinitesimal body;(c)the strain hardening law of unsaturated expansive soil is mainly reflected by fractional derivative operator.
文摘Hydropower gains increasing importance as a steerable and controllable power source in a renewable energy mix and deregulated markets. Although hydropower produces fossil-free energy, it has a significant impact on the local environment. This review investigates the effects of flow alterations by hydropower on the downstream river system and the possibilities to integrate these effects into hydraulic modeling. The results show that various effects of flow regulation on the ecosystem, but also social and economic effects on related communities were observed in the last decades. The application of hydraulic models for investigations of ecological effects is common. Especially hydraulic effects and effects on fish were extensively modeled with the help of hydraulic 1D- and 2D-simulations. Current applications to investigate social and economic effects integrated into hydraulic modeling are meanwhile limited. Approaches to realizing this integration are presented. Further research on the economic valuation of ecosystems and integration of social and economic effects to hydraulic models is necessary to develop holistic tools to support decision-making on sustainable hydropower.
基金Thework issupportedby the Key Scienceand Technology Programof Henan Province(Grant No.222102220104)the Science and Technology Key Project Foundation of Henan Provincial Education Department(Grant No.23A460014)the High Level Talent Foundation of Henan University of Technology(Grant No.2020BS043).
文摘Hydraulic servo system plays an important role in industrial fields due to the advantages of high response,small size-to-power ratio and large driving force.However,inherent nonlinear behaviors and modeling uncertainties are the main obstacles for hydraulic servo system to achieve high tracking perfor-mance.To deal with these difficulties,this paper presents a backstepping sliding mode controller to improve the dynamic tracking performance and anti-interfer-ence ability.For this purpose,the nonlinear dynamic model is firstly established,where the nonlinear behaviors and modeling uncertainties are lumped as one term.Then,the extended state observer is introduced to estimate the lumped distur-bance.The system stability is proved by using the Lyapunov stability theorem.Finally,comparative simulation and experimental are conducted on a hydraulic servo system platform to verify the efficiency of the proposed control scheme.
文摘Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.
基金Supported by National Key Basic Research Program of China(973 Program,Grant No.2014CB046405)Hebei Provincial Applied Basic Research Program(Grant No.12962147D)National Natural Science Foundation of China(Grant No.51375423)
文摘The previous sensitivity analysis researches are not accurate enough and also have the limited reference value, because those mathematical models are relatively simple and the change of the load and the initial displacement changes of the piston are ignored, even experiment verification is not conducted. Therefore, in view of deficiencies above, a nonlinear mathematical model is established in this paper, including dynamic characteristics of servo valve, nonlinear characteristics of pressure-flow, initial displacement of servo cylinder piston and friction nonlinearity. The transfer function block diagram is built for the hydraulic drive unit closed loop position control, as well as the state equations. Through deriving the time-varying coefficient items matrix and time-varying free items matrix of sensitivity equations respectively, the expression of sensitivity equations based on the nonlinear mathematical model are obtained. According to structure parameters of hydraulic drive unit, working parameters, fluid transmission characteristics and measured friction-velocity curves, the simulation analysis of hydraulic drive unit is completed on the MATLAB/Simulink simulation platform with the displacement step 2 mm, 5 mm and 10 mm, respectively. The simulation results indicate that the developed nonlinear mathematical model is sufficient by comparing the characteristic curves of experimental step response and simulation step response under different constant load. Then, the sensitivity function time-history curves of seventeen parameters are obtained, basing on each state vector time-history curve of step response characteristic. The maximum value of displacement variation percentage and the sum of displacement variation absolute values in the sampling time are both taken as sensitivity indexes. The sensitivity indexes values above are calculated and shown visually in histograms under different working conditions, and change rules are analyzed. Then the sensitivity indexes values of four measurable parameters, such as supply pressure, proportional gain, initial position of servo cylinder piston and load force, are verified experimentally on test platform of hydraulic drive unit, and the experimental research shows that the sensitivity analysis results obtained through simulation are approximate to the test results. This research indicates each parameter sensitivity characteristics of hydraulic drive unit, the performance-affected main parameters and secondary parameters are got under different working conditions, which will provide the theoretical foundation for the control compensation and structure optimization of hydraulic drive unit.
基金Financial support provided by the U.S. Department of Energy under DOE Grant No. DE-FE0002760
文摘This paper presents an improved understanding of coupled hydro-thermo-mechanical(HTM) hydraulic fracturing of quasi-brittle rock using the bonded particle model(BPM) within the discrete element method(DEM). BPM has been recently extended by the authors to account for coupled convective econductive heat flow and transport, and to enable full hydro-thermal fluidesolid coupled modeling.The application of the work is on enhanced geothermal systems(EGSs), and hydraulic fracturing of hot dry rock(HDR) is studied in terms of the impact of temperature difference between rock and a flowing fracturing fluid. Micro-mechanical investigation of temperature and fracturing fluid effects on hydraulic fracturing damage in rocks is presented. It was found that fracture is shorter with pronounced secondary microcracking along the main fracture for the case when the convectiveeconductive thermal heat exchange is considered. First, the convection heat exchange during low-viscosity fluid infiltration in permeable rock around the wellbore causes significant rock cooling, where a finger-like fluid infiltration was observed. Second, fluid infiltration inhibits pressure rise during pumping and delays fracture initiation and propagation. Additionally, thermal damage occurs in the whole area around the wellbore due to rock cooling and cold fluid infiltration. The size of a damaged area around the wellbore increases with decreasing fluid dynamic viscosity. Fluid and rock compressibility ratio was found to have significant effect on the fracture propagation velocity.
基金supported by the Research Directorate of the University of Cuenca(DIUC)
文摘Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.
文摘Instead of establishing mathematical hydraulic system models from physical laws usually done with the problems of complex modelling processes, low reliability and practicality caused by large uncertainties, a novel modelling method for a highly nonlinear system of a hydraulic excavator is presented. Based on the data collected in the excavator's arms driving experiments, a data-based excavator dynamic model using Simplified Refined Instrumental Variable (SRIV) identification and estimation algorithms is established. The validity of the proposed data-based model is indirectly demonstrated by the performance of computer simulation and the.real machine motion control exoeriments.
基金This project is supported by Provincial Natural Science Foundation of Shanxi, China(No.l9991062).
文摘Modeling and digital simulation is an effective method to analyze the dynamiccharacteristics of hydraulic system. It is difficult to determine some performance parameters inthe hydraulic system by means of currently used modeling methods. The 'gray-box' modeling method forlarge-scale hydraulic system is introduced. The principle of the method, the submodels of somecomponents and the parameters identification of components or subsystem are researched.
文摘Water resources management usually requires that hydraulic, ecological, and hydrological models be linked. The Hy- drologic Engineering Center River Analysis System (HEC-RAS) hydraulic model and the Hydrologic Engineering Center Geospatial River Analysis System (HEC-GEORAS), imitates flow and water profiles in the Neka river basin’s downstream flood plain. Hydrograph phases studied during the flood seasons of 1986-1999 and from 2002-2004 were used to calibrate and verify the hydraulic model respectively. Simulations of peak flood stages and hydrographs’ evaluations are congruent with studies and observations, with the former showing mean square errors between 4.8 - 10 cm. HECRAS calculations and forecast flood water levels. Nash-Sutcliffe effectiveness (CR3) is more than 0.92 along with elevated levels of water which were created with some effectiveness (CR5) of 0.94 for the validation period. The coupled two models show good performance in the water level modeling.
基金supported by the Key Innovation Team Program of Innovation Talents Promotion Plan by Ministry of Science and Technology of China(Grant No.2016RA4059)National Natural Science Foundation of China(Grant Nos.41672268 and 41772286)。
文摘This paper presents a three-dimensional fully hydro-mechanical coupled distinct element study on fault reactivation and induced seismicity due to hydraulic fracturing injection and subsequent backflow process,based on the geological data in Horn River Basin,Northeast British Columbia,Canada.The modeling results indicate that the maximum magnitude of seismic events appears at the fracturing stage.The increment of fluid volume in the fault determines the cumulative moment and maximum fault slippage,both of which are essentially proportional to the fluid volume.After backflow starts,the fluid near the joint intersection keeps flowing into the critically stressed fault,rather than backflows to the wellbore.Although fault slippage is affected by the changes of both pore pressure and ambient rock stress,their contributions are different at fracturing and backflow stages.At fracturing stage,pore pressure change shows a dominant effect on induced fault slippage.While at backflow stage,because the fault plane is under a critical stress state,any minor disturbance would trigger a fault slippage.The energy analysis indicates that aseismic deformation takes up a majority of the total deformation energy during hydraulic fracturing.A common regularity is found in both fracturing-and backflow-induced seismicity that the cumulative moment and maximum fault slippage are nearly proportional to the injected fluid volume.This study shows some novel insights into interpreting fracturing-and backflowinduced seismicity,and provides useful information for controlling and mitigating seismic hazards due to hydraulic fracturing.
基金financial supports from the National Natural Science Foundation of China(Grant Nos.51925906 and 51988101)the National Key R&D Program of China(Grant No.2018YFC0407001)。
文摘Modeling unsaturated flow in fractured rocks is essential in various subsurface engineering applications,but it remains a great challenge due to the difficulties in determining the unsaturated hydraulic properties of rocks that contain various scales of fractures.It is generally accepted that the van Genuchten(VG)model can be applied to fractured rocks,provided that the hydraulic parameters could be representatively determined.In this study,scaling relationships between the VG parameters(a and n)and hydraulic conductivity(K)across 8 orders of magnitude,from 10^(-10)m/s to 10^(-2)m/s,were proposed by statistical analysis of data obtained from 1416 soil samples.The correlations were then generalized to predict the upper bounds of VG parameters for fractured rocks from the K data that could be obtained more easily under field conditions,and were validated against a limited set of data from cores,fractures and fractured rocks available in the literature.The upper bound estimates significantly narrow the ranges of VG parameters,and the representative values of a and n for fractured rocks at the field scale can then be determined with confidence by inverse modeling using groundwater observations in saturated zones.The proposed methodology was applied to saturated-unsaturated flow modeling in the right-bank slope at the Baihetan dam site with a continuum approach,showing that most of the flow behaviors in fractured rocks in this complex hydrogeological condition could be properly reproduced.The proposed method overcomes difficulties in suction measurement in fractured rocks with strong heterogeneity,and provides a feasible way for modeling of saturated-unsaturated flow in fractured rocks with acceptable engineering accuracy.