Bamboo plants are an essential component of tropical ecosystems,yet their vulnerability to climate extremes,such as drought,is poorly understood due to limited knowledge of their hydraulic properties.Cephalostachyum p...Bamboo plants are an essential component of tropical ecosystems,yet their vulnerability to climate extremes,such as drought,is poorly understood due to limited knowledge of their hydraulic properties.Cephalostachyum pergracile,a commonly used tropical bamboo species,exhibited a substantially higher mortality rate than other co-occurring bamboos during a severe drought event in 2019,but the underlying mechanisms remain unclear.This study investigated the leaf and stem hydraulic traits related to drought responses,including leaf-stem embolism resistance(P50leaf;P50stem) estimated using optical and X-ray microtomography methods,leaf pressure-volume and water-releasing curves.Additionally,we investigated the seasonal water potentials,native embolism level(PLC) and xylem water source using stable isotope.We found that C.pergracile exhibited strong resistance to embolism,showing low P50leaf,P50stem,and turgor loss point,despite its rapid leaf water loss.Interestingly,its leaves displayed greater resistance to embolism than its stem,suggesting a lack of effective hydraulic vulnerability segmentation(HVS) to protect the stem from excessive xylem tension.During the dry season,approximately 49% of the water was absorbed from the upper 20-cm-deep soil layer.Consequently,significant diurnal variation in leaf water potentials and an increase in midday PLC from 5.87±2.33% in the wet season to 12.87±4.09%in the dry season were observed.In summary,this study demonstrated that the rapid leaf water loss,high reliance on surface water,and a lack of effective HVS in C.pergracile accelerated water depletion and increased xylem embolism even in the typical dry season,which may explain its high mortality rate during extreme drought events in 2019.展开更多
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.展开更多
The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determinati...The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.展开更多
The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a...The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.展开更多
The spatiotemporal distributions of microbes in soil by different methods could affect the efficacy of the microbes to reduce the soil hydraulic conductivity.In this study,the specimens of bio-mediated sands were prep...The spatiotemporal distributions of microbes in soil by different methods could affect the efficacy of the microbes to reduce the soil hydraulic conductivity.In this study,the specimens of bio-mediated sands were prepared using three different methods,i.e.injecting,mixing,and pouring a given microbial so-lution onto compacted sand specimens.The hydraulic conductivity was measured by constant-head tests,while any soil microstructural changes due to addition of the microbes were observed by scan-ning electron microscope(SEM)and mercury intrusion porosimetry(MIP)tests.The amount of dextran concentration produced by microbes in each type of specimen was quantified by a refractometer.Results show that dextran production increased exponentially after 5-7 d of microbial settling with the supply of culture medium.The injection and mixing methods resulted in a similar amount and uniform dis-tribution of dextran in the specimens.The pouring method,however,produced a nonuniform distri-bution,with a higher concentration near the specimen surface.As the supply of culture medium discontinued,the dextran content near the surface produced by the pouring method decreased dramatically due to high competition for nutrients with foreign colonies.Average dextran concentration was negatively and correlated with hydraulic conductivity of bio-mediated soils exponentially,due to the clogging of large soil pores by dextran.The hydraulic conductivity of the injection and mixing cases did not change significantly when the supply of culture medium was absent.展开更多
Deep and ultra-deep reservoirs have gradually become the primary focus of hydrocarbon exploration as a result of a series of significant discoveries in deep hydrocarbon exploration worldwide.These reservoirs present u...Deep and ultra-deep reservoirs have gradually become the primary focus of hydrocarbon exploration as a result of a series of significant discoveries in deep hydrocarbon exploration worldwide.These reservoirs present unique challenges due to their deep burial depth(4500-8882 m),low matrix permeability,complex crustal stress conditions,high temperature and pressure(HTHP,150-200℃,105-155 MPa),coupled with high salinity of formation water.Consequently,the costs associated with their exploitation and development are exceptionally high.In deep and ultra-deep reservoirs,hydraulic fracturing is commonly used to achieve high and stable production.During hydraulic fracturing,a substantial volume of fluid is injected into the reservoir.However,statistical analysis reveals that the flowback rate is typically less than 30%,leaving the majority of the fluid trapped within the reservoir.Therefore,hydraulic fracturing in deep reservoirs not only enhances the reservoir permeability by creating artificial fractures but also damages reservoirs due to the fracturing fluids involved.The challenging“three-high”environment of a deep reservoir,characterized by high temperature,high pressure,and high salinity,exacerbates conventional forms of damage,including water sensitivity,retention of fracturing fluids,rock creep,and proppant breakage.In addition,specific damage mechanisms come into play,such as fracturing fluid decomposition at elevated temperatures and proppant diagenetic reactions at HTHP conditions.Presently,the foremost concern in deep oil and gas development lies in effectively assessing the damage inflicted on these reservoirs by hydraulic fracturing,comprehending the underlying mechanisms,and selecting appropriate solutions.It's noteworthy that the majority of existing studies on reservoir damage primarily focus on conventional reservoirs,with limited attention given to deep reservoirs and a lack of systematic summaries.In light of this,our approach entails initially summarizing the current knowledge pertaining to the types of fracturing fluids employed in deep and ultra-deep reservoirs.Subsequently,we delve into a systematic examination of the damage processes and mechanisms caused by fracturing fluids within the context of hydraulic fracturing in deep reservoirs,taking into account the unique reservoir characteristics of high temperature,high pressure,and high in-situ stress.In addition,we provide an overview of research progress related to high-temperature deep reservoir fracturing fluid and the damage of aqueous fracturing fluids to rock matrix,both artificial and natural fractures,and sand-packed fractures.We conclude by offering a summary of current research advancements and future directions,which hold significant potential for facilitating the efficient development of deep oil and gas reservoirs while effectively mitigating reservoir damage.展开更多
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.展开更多
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.展开更多
With an extension of the geological entropy concept in porous media,the approach called directional entrogram is applied to link hydraulic behavior to the anisotropy of the 3D fracture networks.A metric called directi...With an extension of the geological entropy concept in porous media,the approach called directional entrogram is applied to link hydraulic behavior to the anisotropy of the 3D fracture networks.A metric called directional entropic scale is used to measure the anisotropy of spatial order in different directions.Compared with the traditional connectivity indexes based on the statistics of fracture geometry,the directional entropic scale is capable to quantify the anisotropy of connectivity and hydraulic conductivity in heterogeneous 3D fracture networks.According to the numerical analysis of directional entrogram and fluid flow in a number of the 3D fracture networks,the hydraulic conductivities and entropic scales in different directions both increase with spatial order(i.e.,trace length decreasing and spacing increasing)and are independent of the dip angle.As a result,the nonlinear correlation between the hydraulic conductivities and entropic scales from different directions can be unified as quadratic polynomial function,which can shed light on the anisotropic effect of spatial order and global entropy on the heterogeneous hydraulic behaviors.展开更多
Deepsea mining has been proposed since the 1960s to alleviate the lack of resources on land.Vertical hydraulic transport of collected ores from the seabed to the sea surface is considered the most promising method for...Deepsea mining has been proposed since the 1960s to alleviate the lack of resources on land.Vertical hydraulic transport of collected ores from the seabed to the sea surface is considered the most promising method for industrial applications.In the present study,an indoor model test of the vertical hydraulic transport of particles was conducted.A noncontact optical method has been proposed to measure the local characteristics of the particles inside a vertical pipe,including the local concentration and particle velocity.The hydraulic gradient of ore transport was evaluated with various particle size distributions,particle densities,feeding concentrations and mixture flow velocities.During transport,the local concentration is larger than the feeding concentration,whereas the particle velocity is less than the mixture velocity.The qualitative effects of the local concentration and local fluid velocity on the particle velocity and slip velocity were investigated.The local fluid velocity contributes significantly to particle velocity and slip velocity,whereas the effect of the local concentration is marginal.A higher feeding concentration and mixture flow velocity result in an increased hydraulic gradient.The effect of the particle size gradation is slight,whereas the particle density plays a crucial role in the transport.展开更多
A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the dr...A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.展开更多
Perforation is a pivotal technique employed to establish main flow channels within the reservoir formation at the outset of hydraulic fracturing operations.Optimizing perforation designs is critical for augmenting the...Perforation is a pivotal technique employed to establish main flow channels within the reservoir formation at the outset of hydraulic fracturing operations.Optimizing perforation designs is critical for augmenting the efficacy of hydraulic fracturing and boosting oil or gas production.In this study,we employ a hybrid finite-discrete element method,known as the continuous–discontinuous element method(CDEM),to simulate the initiation of post-perforation hydraulic fractures and to derive enhanced design parameters.The model incorporates the four most prevalent perforation geometries,as delineated in an engineering technical report.Real-world perforations deviate from the ideal cylindrical shape,exhibiting variable cross-sectional profiles that typically manifest as an initial constriction followed by an expansion,a feature consistent across all four perforation types.Our simulations take into account variations in perforation hole geometries,cross-sectional diameters,and perforation lengths.The findings show that perforations generated by the 39g DP3 HMX perforating bullet yield the lowest breakdown pressure,which inversely correlates with increases in sectional diameter and perforation length.Moreover,this study reveals the relationship between breakdown pressure and fracture degree,providing valuable insights for engineers and designers to refine perforation strategies.展开更多
Stomatal regulation is critical for mangroves to survive in the hyper-saline intertidal zone where water stress is severe and water availability is highly fluctuant.However,very little is known about the stomatal sens...Stomatal regulation is critical for mangroves to survive in the hyper-saline intertidal zone where water stress is severe and water availability is highly fluctuant.However,very little is known about the stomatal sensitivity to vapour pressure deficit(VPD)in mangroves,and its co-ordination with stomatal morphology and leaf hydraulic traits.We measured the stomatal response to a step increase in VPD in situ,stomatal anatomy,leaf hydraulic vulnerability and pressure-volume traits in nine true mangrove species of five families and collected the data of genome size.We aimed to answer two questions:(1)Does stomatal morphology influence stomatal dynamics in response to a high VPD in mangroves?with a consideration of possible influence of genome size on stomatal morphology;and(2)do leaf hydraulic traits influence stomatal sensitivity to VPD in mangroves?We found that the stomata of mangrove plants were highly sensitive to a step rise in VPD and the stomatal responses were directly affected by stomatal anatomy and hydraulic traits.Smaller,denser stomata was correlated with faster stomatal closure at high VPD across the species of Rhizophoraceae,and stomata size negatively and vein density positively correlated with genome size.Less negative leaf osmotic pressure at the full turgor(πo)was related to higher operating steady-state stomatal conductance(gs);and a higher leaf capacitance(Cleaf)and more embolism resistant leaf xylem were associated with slower stomatal responses to an increase in VPD.In addition,stomatal responsiveness to VPD was indirectly affected by leaf morphological traits,which were affected by site salinity and consequently leaf water status.Our results demonstrate that mangroves display a unique relationship between genome size,stomatal size and vein packing,and that stomatal responsiveness to VPD is regulated by leaf hydraulic traits and stomatal morphology.Our work provides a quantitative framework to better understand of stomatal regulation in mangroves in an environment with high salinity and dynamic water availability.展开更多
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.展开更多
The existence of high-density bedding planes is a typical characteristic of shale oil reservoirs.Understanding the behavior of hydraulic fracturing in high-density laminated rocks is significant for promoting shale oi...The existence of high-density bedding planes is a typical characteristic of shale oil reservoirs.Understanding the behavior of hydraulic fracturing in high-density laminated rocks is significant for promoting shale oil production.In this study,a hydraulic fracturing model considering tensile failure and frictional slip of the bedding planes is established within the framework of the unified pipe-interface element method(UP-IEM).The model developed for simulating the interaction between the hydraulic fracture and the bedding plane is validated by comparison with experimental results.The hydraulic fracturing patterns in sealed and unsealed bedding planes are compared.Additionally,the effects of differential stress,bedding plane permeability,spacing,and the friction coefficient of the bedding plane are investigated.The results showed that a single main fracture crossing the bedding planes is more likely to form in sealed bedding planes under high differential stress.The decrease in bedding plane permeability and the increase in the friction coefficient also promote the fracture propagating perpendicular to the bedding planes.Shale with high-density bedding planes has a poorer fracturing effect than that with low-density bedding planes,as the hydraulic fracture is prone to initiate and propagate along the bedding planes.Moreover,higher injection pressure is needed to maintain fracture propagation along the bedding.An increase in bedding density will lead to a smaller fracturing area.Fracturing fluid seepage into the bedding planes slows shale fracturing.It is recommended that increasing the injection flow rate,selecting alternative fracturing fluids,and employing multi-well/multi-cluster fracturing may be efficient methods to improve energy production in shale oil reservoirs.展开更多
In groundwater hydrology,aquitard heterogeneity is often less considered compared to aquifers,despite its significant impact on groundwater hydraulics and groundwater resources evaluation.A semi-analytical solution is...In groundwater hydrology,aquitard heterogeneity is often less considered compared to aquifers,despite its significant impact on groundwater hydraulics and groundwater resources evaluation.A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity(K)and specific storage(S_(s)).The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or S_(s).Two scenarios of analytical investigations are explored:one is the presence of aquitard interlayers with distinct K or S_(s) values,a common field-scale occurrence;another is an exponentially depth-decaying aquitard S_(s),a regional-scale phenomenon supported by statistical analysis.Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion;a high-S_(s) interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion.Locations of low-K or high-S_(s) interlayers also significantly impact well hydraulics and groundwater budget.In the context of an exponentially depth-decaying aquitard S_(s),a larger decay exponent can enhance aquifer drawdown.When using current models with a vertically homogeneous aquitard,half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard S_(s) should be used to calculate aquitard depletion and unconfined aquifer leakage.展开更多
Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the...Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the concept of large-scale stimulation by fracture network,balanced propagation and effective support of fracture network in fracturing design and developed the extreme massive hydraulic fracturing technique for deep coalbed methane(CBM)horizontal wells.This technique involves massive injection with high pumping rate+high-intensity proppant injection+perforation with equal apertures and limited flow+temporary plugging and diverting fractures+slick water with integrated variable viscosity+graded proppants with multiple sizes.The technique was applied in the pioneering test of a multi-stage fracturing horizontal well in deep CBM of Linxing Block,eastern margin of the Ordos Basin.The injection flow rate is 18 m^(3)/min,proppant intensity is 2.1 m^(3)/m,and fracturing fluid intensity is 16.5 m^(3)/m.After fracturing,a complex fracture network was formed,with an average fracture length of 205 m.The stimulated reservoir volume was 1987×10^(4)m^(3),and the peak gas production rate reached 6.0×10^(4)m^(3)/d,which achieved efficient development of deep CBM.展开更多
Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale...Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.展开更多
Hydraulic jetting is a form of fracturing that involves using a high-pressure jet of water to create fractures in the reservoir rock with a nozzle serving as the central component of the hydraulic sandblasting perfora...Hydraulic jetting is a form of fracturing that involves using a high-pressure jet of water to create fractures in the reservoir rock with a nozzle serving as the central component of the hydraulic sandblasting perforation tool.In this study,the flow behavior of the nozzle is simulated numerically in the framework of a SST k-ωturbulence model.The results show that the nozzle structure can significantly influence the jet performance and related cavitation effect.Through orthogonal experiments,the nozzle geometric parameters are optimized,and the following configuration is found accordingly:contraction angle 20°,contraction segment length 6 mm,cylindrical segment diameter 6 mm,cylindrical segment length 12 mm,spread segment length 10 mm,and spread angle 55°.展开更多
The efficient development and exploitation of shale oil depends on long-distance horizontal wells. As the degreeof cleaning of the wellbore plays a key role in these processes, in this study, this problem is investiga...The efficient development and exploitation of shale oil depends on long-distance horizontal wells. As the degreeof cleaning of the wellbore plays a key role in these processes, in this study, this problem is investigated experimentallyby focusing on the dimensionless cuttings bed height. A method is proposed to calculate the horizontalwellhydraulic extension taking into account the influence of the wellbore cleaning degree on the wellborepressure distribution and assess the effect of a variety of factors such as the bottom hole pressure, the circulatingpressure drop, the drilling pump performance and the formation properties. The analysis shows that the hydraulicextension of horizontal wells decreases with an increase in the cuttings bed height, and the higher the displacementof drilling fluid, the faster the hydraulic extension declines. The annular pressure drop of the horizontalsection increases with the increase of the cuttings bed height, resulting in a higher bottom-hole pressure. Severalarguments are provided to guide the safe drilling of shale oil horizontal wells and overcome the limits of currenttechnological approaches.展开更多
基金supported by the National Natural Science Foundation of China(Nos:32071735,32371576,32350410420,41861144016,and 31570406)CAS‘Light of West China’Program+2 种基金The 14th Five-Year Plan of the Xishuangbanna Tropical Botanical Garden,Chinese Academy of Sciences(E3ZKFF1K,E3ZKFF2B)Yunnan Provincial Science and Technology Department(2018HB068)Yunnan Revitalization Talents Support Plan(YNWR-QNBJ-2019177).
文摘Bamboo plants are an essential component of tropical ecosystems,yet their vulnerability to climate extremes,such as drought,is poorly understood due to limited knowledge of their hydraulic properties.Cephalostachyum pergracile,a commonly used tropical bamboo species,exhibited a substantially higher mortality rate than other co-occurring bamboos during a severe drought event in 2019,but the underlying mechanisms remain unclear.This study investigated the leaf and stem hydraulic traits related to drought responses,including leaf-stem embolism resistance(P50leaf;P50stem) estimated using optical and X-ray microtomography methods,leaf pressure-volume and water-releasing curves.Additionally,we investigated the seasonal water potentials,native embolism level(PLC) and xylem water source using stable isotope.We found that C.pergracile exhibited strong resistance to embolism,showing low P50leaf,P50stem,and turgor loss point,despite its rapid leaf water loss.Interestingly,its leaves displayed greater resistance to embolism than its stem,suggesting a lack of effective hydraulic vulnerability segmentation(HVS) to protect the stem from excessive xylem tension.During the dry season,approximately 49% of the water was absorbed from the upper 20-cm-deep soil layer.Consequently,significant diurnal variation in leaf water potentials and an increase in midday PLC from 5.87±2.33% in the wet season to 12.87±4.09%in the dry season were observed.In summary,this study demonstrated that the rapid leaf water loss,high reliance on surface water,and a lack of effective HVS in C.pergracile accelerated water depletion and increased xylem embolism even in the typical dry season,which may explain its high mortality rate during extreme drought events in 2019.
基金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.
基金supported by the National Natural Science Foundation of China(Grant No.42174118)a research grant(Grant No.ZDJ 2020-7)from the National Institute of Natural Hazards,Ministry of Emergency Management of China.
文摘The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.
基金Funding support from Heilongjiang"Open Competition"project(Grant No.DQYT2022-JS-758)is greatly acknowledgedFinancial support from the National Natural Science Foundation of China(Grant Nos.52304025 and 52174025)is acknowledged+1 种基金supports from Northeast Petroleum University and Guangdong Basic and Applied Basic Research Foundationsupport from the Heilongjiang Touyan Innovation Team Program.
文摘The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.
基金The first author(V.Kamchoom)acknowledges the grant(Grant No.FRB66065/0258-RE-KRIS/FF66/53)from King Mongkut’s Insti-tute of Technology Ladkrabang(KMITL)and National Science,Research and Innovation Fund(NSRF)the grant under Climate Change and Climate Variability Research in Monsoon Asia(CMON3)from the National Research Council of Thailand(NRCT)(Grant No.N10A650844)the National Natural Science Foundation of China(NSFC).
文摘The spatiotemporal distributions of microbes in soil by different methods could affect the efficacy of the microbes to reduce the soil hydraulic conductivity.In this study,the specimens of bio-mediated sands were prepared using three different methods,i.e.injecting,mixing,and pouring a given microbial so-lution onto compacted sand specimens.The hydraulic conductivity was measured by constant-head tests,while any soil microstructural changes due to addition of the microbes were observed by scan-ning electron microscope(SEM)and mercury intrusion porosimetry(MIP)tests.The amount of dextran concentration produced by microbes in each type of specimen was quantified by a refractometer.Results show that dextran production increased exponentially after 5-7 d of microbial settling with the supply of culture medium.The injection and mixing methods resulted in a similar amount and uniform dis-tribution of dextran in the specimens.The pouring method,however,produced a nonuniform distri-bution,with a higher concentration near the specimen surface.As the supply of culture medium discontinued,the dextran content near the surface produced by the pouring method decreased dramatically due to high competition for nutrients with foreign colonies.Average dextran concentration was negatively and correlated with hydraulic conductivity of bio-mediated soils exponentially,due to the clogging of large soil pores by dextran.The hydraulic conductivity of the injection and mixing cases did not change significantly when the supply of culture medium was absent.
基金Dao-Bing Wang was supported by the Beijing Natural Science Foundation Project(No.3222030)the National Natural Science Foundation of China(No.52274002)+1 种基金the PetroChina Science and Technology Innovation Foundation Project(No.2021DQ02-0201)Fu-Jian Zhou was supported by the National Natural Science Foundation of China(No.52174045).
文摘Deep and ultra-deep reservoirs have gradually become the primary focus of hydrocarbon exploration as a result of a series of significant discoveries in deep hydrocarbon exploration worldwide.These reservoirs present unique challenges due to their deep burial depth(4500-8882 m),low matrix permeability,complex crustal stress conditions,high temperature and pressure(HTHP,150-200℃,105-155 MPa),coupled with high salinity of formation water.Consequently,the costs associated with their exploitation and development are exceptionally high.In deep and ultra-deep reservoirs,hydraulic fracturing is commonly used to achieve high and stable production.During hydraulic fracturing,a substantial volume of fluid is injected into the reservoir.However,statistical analysis reveals that the flowback rate is typically less than 30%,leaving the majority of the fluid trapped within the reservoir.Therefore,hydraulic fracturing in deep reservoirs not only enhances the reservoir permeability by creating artificial fractures but also damages reservoirs due to the fracturing fluids involved.The challenging“three-high”environment of a deep reservoir,characterized by high temperature,high pressure,and high salinity,exacerbates conventional forms of damage,including water sensitivity,retention of fracturing fluids,rock creep,and proppant breakage.In addition,specific damage mechanisms come into play,such as fracturing fluid decomposition at elevated temperatures and proppant diagenetic reactions at HTHP conditions.Presently,the foremost concern in deep oil and gas development lies in effectively assessing the damage inflicted on these reservoirs by hydraulic fracturing,comprehending the underlying mechanisms,and selecting appropriate solutions.It's noteworthy that the majority of existing studies on reservoir damage primarily focus on conventional reservoirs,with limited attention given to deep reservoirs and a lack of systematic summaries.In light of this,our approach entails initially summarizing the current knowledge pertaining to the types of fracturing fluids employed in deep and ultra-deep reservoirs.Subsequently,we delve into a systematic examination of the damage processes and mechanisms caused by fracturing fluids within the context of hydraulic fracturing in deep reservoirs,taking into account the unique reservoir characteristics of high temperature,high pressure,and high in-situ stress.In addition,we provide an overview of research progress related to high-temperature deep reservoir fracturing fluid and the damage of aqueous fracturing fluids to rock matrix,both artificial and natural fractures,and sand-packed fractures.We conclude by offering a summary of current research advancements and future directions,which hold significant potential for facilitating the efficient development of deep oil and gas reservoirs while effectively mitigating reservoir damage.
基金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.
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.42077243,52209148,and 52079062).
文摘With an extension of the geological entropy concept in porous media,the approach called directional entrogram is applied to link hydraulic behavior to the anisotropy of the 3D fracture networks.A metric called directional entropic scale is used to measure the anisotropy of spatial order in different directions.Compared with the traditional connectivity indexes based on the statistics of fracture geometry,the directional entropic scale is capable to quantify the anisotropy of connectivity and hydraulic conductivity in heterogeneous 3D fracture networks.According to the numerical analysis of directional entrogram and fluid flow in a number of the 3D fracture networks,the hydraulic conductivities and entropic scales in different directions both increase with spatial order(i.e.,trace length decreasing and spacing increasing)and are independent of the dip angle.As a result,the nonlinear correlation between the hydraulic conductivities and entropic scales from different directions can be unified as quadratic polynomial function,which can shed light on the anisotropic effect of spatial order and global entropy on the heterogeneous hydraulic behaviors.
基金financially supported by the Hainan Provincial Joint Project of Sanya Yazhou Bay Science and Technology City(Grant No.520LH052)the National Natural Science Foundation of China(Grant No.51909164).
文摘Deepsea mining has been proposed since the 1960s to alleviate the lack of resources on land.Vertical hydraulic transport of collected ores from the seabed to the sea surface is considered the most promising method for industrial applications.In the present study,an indoor model test of the vertical hydraulic transport of particles was conducted.A noncontact optical method has been proposed to measure the local characteristics of the particles inside a vertical pipe,including the local concentration and particle velocity.The hydraulic gradient of ore transport was evaluated with various particle size distributions,particle densities,feeding concentrations and mixture flow velocities.During transport,the local concentration is larger than the feeding concentration,whereas the particle velocity is less than the mixture velocity.The qualitative effects of the local concentration and local fluid velocity on the particle velocity and slip velocity were investigated.The local fluid velocity contributes significantly to particle velocity and slip velocity,whereas the effect of the local concentration is marginal.A higher feeding concentration and mixture flow velocity result in an increased hydraulic gradient.The effect of the particle size gradation is slight,whereas the particle density plays a crucial role in the transport.
基金Project supported by the National Natural Science Foundation of China(No.42202314)。
文摘A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.
基金support from the National Natural Science Foundation of China(Grant Nos.52178324,12102059)the China Postdoctoral Science Foundation(Grant No.2023M743604)+1 种基金the Beijing Natural Science Foundation(Grant No.3212027),the National Key R&D Program of China(Grant No.2023YFC3007203)the 2019 Foreign Experts Plan of Hebei Province.
文摘Perforation is a pivotal technique employed to establish main flow channels within the reservoir formation at the outset of hydraulic fracturing operations.Optimizing perforation designs is critical for augmenting the efficacy of hydraulic fracturing and boosting oil or gas production.In this study,we employ a hybrid finite-discrete element method,known as the continuous–discontinuous element method(CDEM),to simulate the initiation of post-perforation hydraulic fractures and to derive enhanced design parameters.The model incorporates the four most prevalent perforation geometries,as delineated in an engineering technical report.Real-world perforations deviate from the ideal cylindrical shape,exhibiting variable cross-sectional profiles that typically manifest as an initial constriction followed by an expansion,a feature consistent across all four perforation types.Our simulations take into account variations in perforation hole geometries,cross-sectional diameters,and perforation lengths.The findings show that perforations generated by the 39g DP3 HMX perforating bullet yield the lowest breakdown pressure,which inversely correlates with increases in sectional diameter and perforation length.Moreover,this study reveals the relationship between breakdown pressure and fracture degree,providing valuable insights for engineers and designers to refine perforation strategies.
基金financially supported by a grant of the National Natural Science Foundation of China(31670406)the Bagui Fellow scholarship(C33600992001)of Guangxi Zhuang Autonomous Region to KFC.
文摘Stomatal regulation is critical for mangroves to survive in the hyper-saline intertidal zone where water stress is severe and water availability is highly fluctuant.However,very little is known about the stomatal sensitivity to vapour pressure deficit(VPD)in mangroves,and its co-ordination with stomatal morphology and leaf hydraulic traits.We measured the stomatal response to a step increase in VPD in situ,stomatal anatomy,leaf hydraulic vulnerability and pressure-volume traits in nine true mangrove species of five families and collected the data of genome size.We aimed to answer two questions:(1)Does stomatal morphology influence stomatal dynamics in response to a high VPD in mangroves?with a consideration of possible influence of genome size on stomatal morphology;and(2)do leaf hydraulic traits influence stomatal sensitivity to VPD in mangroves?We found that the stomata of mangrove plants were highly sensitive to a step rise in VPD and the stomatal responses were directly affected by stomatal anatomy and hydraulic traits.Smaller,denser stomata was correlated with faster stomatal closure at high VPD across the species of Rhizophoraceae,and stomata size negatively and vein density positively correlated with genome size.Less negative leaf osmotic pressure at the full turgor(πo)was related to higher operating steady-state stomatal conductance(gs);and a higher leaf capacitance(Cleaf)and more embolism resistant leaf xylem were associated with slower stomatal responses to an increase in VPD.In addition,stomatal responsiveness to VPD was indirectly affected by leaf morphological traits,which were affected by site salinity and consequently leaf water status.Our results demonstrate that mangroves display a unique relationship between genome size,stomatal size and vein packing,and that stomatal responsiveness to VPD is regulated by leaf hydraulic traits and stomatal morphology.Our work provides a quantitative framework to better understand of stomatal regulation in mangroves in an environment with high salinity and dynamic water availability.
基金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.
基金The authors wish to acknowledge the financial support from Key Laboratory of Deep Earth Science and Engineering(Sichuan University),Ministry of Education(DESE202202,H.Y)State Energy Center for Shale Oil Research and Development(33550000-22-ZC0613-0365,H.Y)+2 种基金National Natural Science Foundation of China(42307209,X.Y)China Postdoctoral Science Foundation(2022M712425,X.Y)Shanghai Pujiang Program(2022PJD076,X.Y).
文摘The existence of high-density bedding planes is a typical characteristic of shale oil reservoirs.Understanding the behavior of hydraulic fracturing in high-density laminated rocks is significant for promoting shale oil production.In this study,a hydraulic fracturing model considering tensile failure and frictional slip of the bedding planes is established within the framework of the unified pipe-interface element method(UP-IEM).The model developed for simulating the interaction between the hydraulic fracture and the bedding plane is validated by comparison with experimental results.The hydraulic fracturing patterns in sealed and unsealed bedding planes are compared.Additionally,the effects of differential stress,bedding plane permeability,spacing,and the friction coefficient of the bedding plane are investigated.The results showed that a single main fracture crossing the bedding planes is more likely to form in sealed bedding planes under high differential stress.The decrease in bedding plane permeability and the increase in the friction coefficient also promote the fracture propagating perpendicular to the bedding planes.Shale with high-density bedding planes has a poorer fracturing effect than that with low-density bedding planes,as the hydraulic fracture is prone to initiate and propagate along the bedding planes.Moreover,higher injection pressure is needed to maintain fracture propagation along the bedding.An increase in bedding density will lead to a smaller fracturing area.Fracturing fluid seepage into the bedding planes slows shale fracturing.It is recommended that increasing the injection flow rate,selecting alternative fracturing fluids,and employing multi-well/multi-cluster fracturing may be efficient methods to improve energy production in shale oil reservoirs.
基金financially supported by the National Key Research and Development Program of China(Grant No.2019YFC1804301)the National Science Fourdation of China(Grant No.42272279,41902244)partial support from a Discovery Grant awarded by the Natural Sciences and Engineering Research Council of Canada(NSERC)。
文摘In groundwater hydrology,aquitard heterogeneity is often less considered compared to aquifers,despite its significant impact on groundwater hydraulics and groundwater resources evaluation.A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity(K)and specific storage(S_(s)).The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or S_(s).Two scenarios of analytical investigations are explored:one is the presence of aquitard interlayers with distinct K or S_(s) values,a common field-scale occurrence;another is an exponentially depth-decaying aquitard S_(s),a regional-scale phenomenon supported by statistical analysis.Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion;a high-S_(s) interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion.Locations of low-K or high-S_(s) interlayers also significantly impact well hydraulics and groundwater budget.In the context of an exponentially depth-decaying aquitard S_(s),a larger decay exponent can enhance aquifer drawdown.When using current models with a vertically homogeneous aquitard,half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard S_(s) should be used to calculate aquitard depletion and unconfined aquifer leakage.
基金Supported by the National Natural Science Foundation of China Project(52274014)Comprehensive Scientific Research Project of China National Offshore Oil Corporation(KJZH-2023-2303)。
文摘Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the concept of large-scale stimulation by fracture network,balanced propagation and effective support of fracture network in fracturing design and developed the extreme massive hydraulic fracturing technique for deep coalbed methane(CBM)horizontal wells.This technique involves massive injection with high pumping rate+high-intensity proppant injection+perforation with equal apertures and limited flow+temporary plugging and diverting fractures+slick water with integrated variable viscosity+graded proppants with multiple sizes.The technique was applied in the pioneering test of a multi-stage fracturing horizontal well in deep CBM of Linxing Block,eastern margin of the Ordos Basin.The injection flow rate is 18 m^(3)/min,proppant intensity is 2.1 m^(3)/m,and fracturing fluid intensity is 16.5 m^(3)/m.After fracturing,a complex fracture network was formed,with an average fracture length of 205 m.The stimulated reservoir volume was 1987×10^(4)m^(3),and the peak gas production rate reached 6.0×10^(4)m^(3)/d,which achieved efficient development of deep CBM.
基金Supported by the National Natural Science Foundation of China(52374043)Key Program of the National Natural Science Foundation of China(52234003).
文摘Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.
基金The authors gratefully acknowledge the financial support by the Marine Economy Development Foundation of Guangdong Province“Technical Support for Stimulation and Testing of Gas Hydrate Reservoirs”(GDNRC[2022]44).
文摘Hydraulic jetting is a form of fracturing that involves using a high-pressure jet of water to create fractures in the reservoir rock with a nozzle serving as the central component of the hydraulic sandblasting perforation tool.In this study,the flow behavior of the nozzle is simulated numerically in the framework of a SST k-ωturbulence model.The results show that the nozzle structure can significantly influence the jet performance and related cavitation effect.Through orthogonal experiments,the nozzle geometric parameters are optimized,and the following configuration is found accordingly:contraction angle 20°,contraction segment length 6 mm,cylindrical segment diameter 6 mm,cylindrical segment length 12 mm,spread segment length 10 mm,and spread angle 55°.
基金supported by CNPC Key Core Technology Research Projects (2022ZG06)project funded by China Postdoctoral Science Foundation (2021M693508)Basic research and strategic reserve technology research fund project of institutes directly under CNPC.
文摘The efficient development and exploitation of shale oil depends on long-distance horizontal wells. As the degreeof cleaning of the wellbore plays a key role in these processes, in this study, this problem is investigated experimentallyby focusing on the dimensionless cuttings bed height. A method is proposed to calculate the horizontalwellhydraulic extension taking into account the influence of the wellbore cleaning degree on the wellborepressure distribution and assess the effect of a variety of factors such as the bottom hole pressure, the circulatingpressure drop, the drilling pump performance and the formation properties. The analysis shows that the hydraulicextension of horizontal wells decreases with an increase in the cuttings bed height, and the higher the displacementof drilling fluid, the faster the hydraulic extension declines. The annular pressure drop of the horizontalsection increases with the increase of the cuttings bed height, resulting in a higher bottom-hole pressure. Severalarguments are provided to guide the safe drilling of shale oil horizontal wells and overcome the limits of currenttechnological approaches.