Evaluating the stability and volumetric flowback rate of proppant packs in hydraulic fractures using the lattice Boltzmann-discrete element coupling method

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.

Estimation of the anisotropy of hydraulic conductivity through 3D fracture networks using the directional geological entropy

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.

Extreme massive hydraulic fracturing in deep coalbed methane horizontal wells:A case study of the Linxing Block,eastern Ordos Basin,NW China

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.

Stress tensor determination by modified hydraulic tests on pre-existing fractures:Method and stress constraints

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.

Influence of Blade Number on the Performance of Hydraulic Turbines in the Transition Stage

To analyze the effect of blade number on the performance of hydraulic turbines during the transient stage in which theflow rate is not constant,six hydraulic turbines with different blade numbers are considered.The instantaneous hydraulic performance of the turbine and the pressure pulsation acting on the impeller are investigated numerically by using the ANSYS CFX software.The ensuing results are compared with the outcomes of experimental tests.It is shown that thefluctuation range of the pressure coefficient increases with time,but the corresponding range for the transient hydraulic efficiency decreases gradually when theflow velocity transits to larger values.During the transition to smallflow velocity,thefluctuation range of the pressure coefficient gradually decreases as time passes,but the correspondingfluctuation range of its transient hydraulic efficiency gradually becomes larger.Thefluctuation range in the Z9 case is small during the transition.The main frequency of transient hydraulic efficiency pulsation is equal to the blade frequency.At the main frequency,Z7 has the largest amplitude of the hydraulic efficiency pulsation,Z10 has the smallest amplitude,and the difference between Z7 and Z9 is limited.As the number of blades grows,the pressure pulsation during the transition process gradually decreases,but the pressure pulsation of Z10 at the volute tongue is larger.In the steady state,Z9 has the highest efficiency and in the transient stage,the pressure coefficientfluctuation range is small.Accordingly,for the hydraulic turbine Z9,the performance is optimal.

Hydraulic fracturing behaviors of shale under coupled stress and temperature conditions simulating different burial depths

Fracture propagation in shale under in situ conditions is a critical but poorly understood mechanical process in hydraulic fracturing for deep shale gas reservoirs. To address this, hydraulic fracturing experiments were conducted on hollow double-wing crack specimens of the Longmaxi shale under conditions simulating the ground surface(confining pressure σ_(cp)=0, room temperature(Tr)) and at depths of 1600 m(σ_(cp)=40 MPa, Ti=70 ℃) and 3300 m(σ_(cp)=80 MPa, high temperature Ti=110 ℃) in the study area.High in situ stress was found to significantly increase fracture toughness through constrained microcracking and particle frictional bridging mechanisms. Increasing the temperature enhances rather than weakens the fracture resistance because it increases the grain debonding length, which dissipates more plastic energy and enlarges grains to close microdefects and generate compressive stress to inhibit microcracking. Interestingly, the fracture toughness anisotropy in the shale was found to be nearly constant across burial depths, despite reported variations with increasing confining pressure. Heated water was not found to be as important as the in situ environment in influencing shale fracture. These findings emphasize the need to test the fracture toughness of deep shales under coupled in situ stress and temperature conditions rather than focusing on either in situ stress or temperature alone.

Pumping-induced Well Hydraulics and Groundwater Budget in a Leaky Aquifer System with Vertical Heterogeneity in Aquitard Hydraulic Properties

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.

Unraveling the hydraulic properties of loess for landslide prediction:A study on variations in loess landslides in Lanzhou,Dingxi,and Tianshui,China

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.

A phase-field model for simulating the propagation behavior of mixed-mode cracks during the hydraulic fracturing process in fractured reservoirs

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.

Spatiotemporal variations of sand hydraulic conductivity by microbial application methods

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.

Effect mechanism of seepage force on the hydraulic fracture propagation

The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a notable influence on rock deformation and failure,its effect on hydraulic fracture(HF)propagation remains ambiguous.Therefore,in this study,we improved a traditional fluid–solid coupling method by incorporating the role of seepage force during the fracturing fluid seepage,using the discrete element method.First,we validated the simulation results of the improved method by comparing them with an analytical solution of the seepage force and published experimental results.Next,we conducted numerical simulations in both homogeneous and heterogeneous sandstone formations to investigate the influence of seepage force on HF propagation.Our results indicate that fluid viscosity has a greater impact on the magnitude and extent of seepage force compared to injection rate,and that lower viscosity and injection rate correspond to shorter hydraulic fracture lengths.Furthermore,seepage force influences the direction of HF propagation,causing HFs to deflect towards the side of the reservoir with weaker cementation and higher permeability.

Reshaping force for deformed casing repairing with hydraulic rolling reshaper and its influencing factors

Hydraulic rolling reshaper is an advanced reshaping tool to solve the problem of casing deformation,which has been widely used in recent years.When it is used for well repair operation,the reshaping force provided by ground devices is generally determined by experience.However,too large reshaping force may destroy the deformed casing,and too small reshaping force may also prolong the construction period and affect the repairing effect.In this paper,based on Hertz contact theory and elastic-plastic theory,combined with the process parameters of shaping,and considering the structural characteristics of the deformed casing and reshaper,we propose a mathematical model for calculating the reshaping force required for repairing deformed casing by hydraulic rolling reshaper.Meanwhile,the finite element model and numerical method of hydraulic rolling reshaper repairing deformed casing are established by using the finite element method,and the reliability of the mathematical model is verified by several examples.On this basis,the control variable method is used to investigate the influence of each parameter on the reshaping force,and the influence degree of each parameter is explored by orthogonal simulation test and Pearson correlation analysis.The research results not only provide an important theoretical basis for the prediction of reshaping force in on-site construction,but also provide a reference for the subsequent improvement of the shaping process.

Numerical Analysis of Perforation during Hydraulic Fracture Initiation Based on Continuous-Discontinuous Element Method

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.

Study of a Hydraulic Jump in an Asymmetric Trapezoidal Channel with Different Sluice Gates

In this study,the main properties of the hydraulic jump in an asymmetric trapezoidal flume are analyzed experimentally,including the so-called sequent depths,characteristic lengths,and efficiency.In particular,an asymmetric trapezoidal flume with a length of 7 m and a width of 0.304 m is considered,with the bottom of the flume transversely inclined at an angle of m=0.296 and vertical lateral sides.The corresponding inflow Froude number is allowed to range in the interval(1.40<F1<6.11).The properties of this jump are compared to those of hydraulic jumps in channels with other types of cross-sections.A relationship for calculating hydraulic jump efficiency is proposed for the considered flume.For F1>5,the hydraulic jump is found to be more effective than that occurring in triangular and symmetric trapezoidal channels.Also,when■mes>8 and■>5,the hydraulic jump in the asymmetrical trapezoidal channel downstream of a parallelogram sluice gate is completely formed as opposed to the situation where a triangular sluice is considered.

Numerical analysis of hydraulic fracture propagation in deep shale reservoir with different injection strategies

Deep shale reservoirs are characterized by elevated breakdown pressures,diminished fracture complexity,and reduced modified volumes compared to medium and shallow reservoirs.Therefore,it is urgent to investigate particular injection strategies that can optimize breakdown pressure and fracturing efficiency to address the increasing demands for deep shale reservoir stimulation.In this study,the efficiency of various stimulation strategies,including multi-cluster simultaneous fracturing,modified alternating fracturing,alternating shut-in fracturing,and cyclic alternating fracturing,was evaluated.Subsequently,the sensitivity of factors such as the cycle index,shut-in time,cluster spacing,and horizontal permeability was investigated.Additionally,the flow distribution effect within the wellbore was discussed.The results indicate that relative to multi-cluster simultaneous fracturing,modified alternating fracturing exhibits reduced susceptibility to the stress shadow effect,which results in earlier breakdown,extended hydraulic fracture lengths,and more consistent propagation despite an increase in breakdown pressure.The alternating shut-in fracturing benefits the increase of fracture length,which is closely related to the shut-in time.Furthermore,cyclic alternating fracturing markedly lowers breakdown pressure and contributes to uniform fracture propagation,in which the cycle count plays an important role.Modified alternating fracturing demonstrates insensitivity to variations in cluster spacing,whereas horizontal permeability is a critical factor affecting fracture length.The wellbore effect restrains the accumulation of pressure and flow near the perforation,delaying the initiation of hydraulic fractures.The simulation results can provide valuable numerical insights for optimizing injection strategies for deep shale hydraulic fracturing.

Vertical distribution characteristics of soil saturated hydraulic conductivity under different land use patterns in the Mu Us sandy land

Background,aim,and scope Soil saturated hydraulic conductivity(K_(s))is a key parameter in the hydrological cycle of soil;however,we have very limited understanding of K_(s) characteristics and the factors that inf luence this key parameter in the Mu Us sandy land(MUSL).Quantifying the impact of changes in land use in the Mu Us sandy land on K_(s) will provide a key foundation for understanding the regional water cycle,but will also provide a scientific basis for the governance of the MUSL.Materials and methods In this study,we determined K_(s) and the basic physical and chemical properties of soil(i.e.,organic matter,bulk density,and soil particle composition)within the first 100 cm layer of four different land use patterns(farmland,tree,shrub,and grassland)in the MUSL.The vertical variation of K_(s) and the factors that influence this key parameter were analyzed and a transfer function for estimating K_(s) was established based on a multiple stepwise regression model.Results The K_(s) of farmland,tree,and shrub increased gradually with soil depth while that of grassland remained unchanged.The K_(s) of the four patterns of land use were moderately variable;mean K_(s)values were ranked as follows:grassland(1.38 mm·min^(-1))<tree(1.76 mm·min^(-1))<farmland(1.82 mm·min^(-1))<shrub(3.30 mm·min^(-1)).The correlation between K_(s) and organic matter,bulk density,and soil particle composition,varied across different land use patterns.A multiple stepwise regression model showed that silt,coarse sand,bulk density,and organic matter,were key predictive factors for the K_(s) of farmland,tree,shrub,and grassland,in the MUSL.Discussion The vertical distribution trend for K_(s) in farmland is known to be predominantly influenced by cultivation,fertilization,and other factors.The general aim is to improve the water-holding capacity of shallow soil on farmland(0-30 cm in depth)to conserve water and nutrients;research has shown that the K_(s) of farmland increases with soil depth.The root growth of tree and shrub in sandy land exerts mechanical force on the soil due to biophysical processes involving rhizospheres,thus leading to a significant change in K_(s).We found that shallow high-density fine roots increased the volume of soil pores and eliminated large pores,thus resulting in a reduction in shallow K_(s).Therefore,the K_(s) of tree and shrub increased with soil depth.Analysis also showed that the K_(s) of grassland did not change significantly and exhibited the lowest mean value when compared to other land use patterns.This finding was predominantly due to the shallow root system of grasslands and because this land use pattern is not subject to human activities such as cultivation and fertilization;consequently,there was no significant change in K_(s) with depth;grassland also had the lowest mean K_(s).We also established a transfer function for K_(s) for different land use patterns in the MUSL.However,the predictive factors for K_(s) in different land use patterns are known to be affected by soil cultivation methods,vegetation restoration modes,the distribution of soil moisture,and other factors,thus resulting in key differences.Therefore,when using the transfer function to predict K_(s) in other areas,it will be necessary to perform parameter calibration and further verification.Conclusions In the MUSL,the K_(s) of farmland,tree,and shrub gradually increased with soil depth;however,the K_(s) of grassland showed no significant variation in terms of vertical distribution.The mean K_(s) values of different land use patterns were ranked as follows:shrub>farmland>tree>grassland;all land use patterns showed moderate levels of variability.The K_(s) for different land use patterns exhibited differing degrees of correlation with soil physical and chemical properties;of these,clay,silt,sand,bulk density,and organic matter,were identified as important variables for predicting K_(s) in farmland,tree,shrub,and grassland,respectively.Recommendations and perspectives In this study,we used a stepwise multiple regression model to establish a transfer function prediction model for K_(s) for different land use patterns;this model possessed high estimation accuracy.The ability to predict K_(s) in the MUSL is very important in terms of the conservation of water and nutrients.

Geohydraulic Investigation of Aquifer Parameters in Abak, Southern Nigeria

A geophysical investigation was carried out to characterize aquifer parameters and assess the groundwater condition in Abak, southern Nigeria. Vertical electrical sounding (VES) was carried out using Schlumberger configuration in 10 locations within the area. Both manual and computer interpretation of the resistivity data reveals three to four geoelectric units (laterite topsoil, medium-grained sand, coarse-grained sand, and sandy clay sand) which agrees with the lithologic log from existing boreholes. The aquifer layer was identified along the third formation with resistivity values ranging from 1239 - 5719 Ωm and aquifer depth ranging from 30.2 - 54.8 m. The aquifer thickness ranged from 24.2 - 43.7 m. Hydraulic conductivity ranged from 0.1206 m/day - 0.5026 m/day with an average value of 0.29403 m/day. Formation factor ranged from 14.55 - 16.64. Porosity ranged from 15.98% - 22.40%, with an average value of 19.64%. The aquiferous zone falls within the medium/coarse-grained sand. The aquifer is shallow, unconfined, and prolific with a little overburden formation of 13.5 m. The area shows good prospects for groundwater development.

Effects of gravel on the water absorption characteristics and hydraulic parameters of stony soil

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.

Extended finite element-based cohesive zone method for modeling simultaneous hydraulic fracture height growth in layered reservoirs

In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.

补肾固胎方联合黄体酮注射液治疗肾虚型早期先兆流产疗效及对妊娠激素水平和Th1/Th2细胞因子的影响

目的 观察补肾固胎方联合黄体酮注射液治疗肾虚型早期先兆流产疗效及对妊娠激素水平和Th1/Th2细胞因子的影响。方法 采用随机数字表法将符合入选标准的120例肾虚型早期先兆流产患者分为观察组和对照组。对照组予黄体酮注射液治疗,观察组同联合补肾固胎方治疗。连续治疗2周后,比较两组中医证候积分(妊娠期阴道流血、小腹坠胀或疼痛、腰酸痛、膝关节酸软)、妊娠激素水平[雌二醇(Estradiol, E_(2))、人绒毛膜促性腺激素(Human chorionic gonadotropin, β-HCG)、孕酮(Progesterone, P)]、Th1/Th2细胞因子[白细胞介素2(interleukin-2,IL-2)、白细胞介素4(interleukin-4,IL-4)、γ干扰素(Interferon-γ,INF-γ)]改善情况,统计两组治疗总有效率、妊娠结局及不良反应发生率。结果 两组治疗后妊娠期阴道流血、小腹坠胀或疼痛、腰酸痛、膝关节酸软积分,IL-2、INF-γ水平均低于治疗前(P<0.05),且观察组治疗后以上中医证候积分和IL-2、INF-γ水平均低于对照组(P<0.05);两组治疗后IL-4水平均高于治疗前(P<0.05),且观察组治疗后IL-4水平高于对照组(P<0.05);观察组总有效率为(91.7%,55/60),高于对照组(75.0%,45/60)(P<0.05);观察组足月分娩率为(85.0%,51/60),高于对照组(56.7%,34/60)(P<0.05);观察组不良反应发生率为3.3%,低于对照组(25.0%)(P<0.05)。结论 补肾固胎方联合黄体酮注射液能够有效调节肾虚型早期先兆流产患者妊娠激素水平和Th1/Th2细胞因子水平,纠正免疫功能紊乱,降低患者流产率,提高足月妊娠率,改善妊娠结局。