Numerical investigation of geostress influence on the grouting reinforcement effectiveness of tunnel surrounding rock mass in fault fracture zones

Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.

Creep constitutive model considering nonlinear creep degradation of fractured rock

Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fractured at 30° and 45° angles. The experimental results indicate that the steady creep strain rates of intact and fractured rock present an exponential increase trend with the increase of stress level. A nonlinear creep model is developed based on the experimental results, in which the initial damage caused by fracture together with the damage caused by constant load have been taken into consideration. The fitting analysis results indicated that the model proposed is more accurate at identifying the full creep regions in fractured granite, especially the accelerated stage of creep deformation. The least-square fit error of the proposed creep model is significantly lower than that of Nishihara model by almost an order of magnitude. An analysis of the effects of elastic modulus, viscosity coefficient, and damage factors on fractured rock strain rate and creep strain is conducted. If no consideration is given to the effects of the damage, the proposed nonlinear creep model can degenerate into to the classical Nishihara model.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Borehole stability in naturally fractured rocks with drilling mud intrusion and associated fracture strength weakening:A coupled DFN-DEM approach

Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.

Multiple damage zones around hydraulic fractures generated by high-frequency pulsating hydraulic fracturing

Pulsating hydraulic fracturing(PHF)is a promising fracturing method and can generate a dynamic periodic pressure.The periodic pressure can induce fatigue failure of rocks and decrease initiation pressure of fracture.If the frequency of periodic pressure exceeds 10 Hz,the distribution of pressure along the main fracture will be heterogeneous,which is much different from the one induced by the common fracturing method.In this study,the impact of this special spatial feature of pressure on hydraulic fracture is mainly investigated.A coupled numerical simulation model is first proposed and verified through experimental and theoretical solutions.The mechanism of secondary fracture initiation around the main fracture is then discovered.In addition,sensitivity studies are conducted to find out the application potential of this new method.The results show that(1)this coupled numerical simulation model is accurate.Through comparison with experimental and theoretical data,the average error of this coupled model is less than 1.01%.(2)Even if a reservoir has no natural fracture,this heterogeneous distribution pressure can also cause many secondary fractures around the main fracture.(3)The mechanism of secondary fracture initiation is that this heterogeneous distribution pressure causes tensile stress at many locations along the main fracture.(4)Through adjusting the stimulation parameters,the stimulation efficiency can be improved.The average and amplitude of pressure can increase possibility of secondary fracture initiation.The frequency of this periodic pressure can increase number of secondary fractures.Even 6 secondary fractures along a 100 m-length main fracture can be generated.(5)The influence magnitudes of stimulation parameters are larger than ones of geomechanical properties,therefore,this new fracturing method has a wide application potential.

Dynamic Identification and Evolution of Urban-suburban-rural Transition Zones Based on the Blender of Natural and Humanistic Factors: A Case Study of Chengdu, China

Urban-suburban-rural(U-S-R)zones exhibit distinctive transitional characteristics in interaction between human and nature.U-S-R transition zones(U-S-RTZ)are also highlighting the function diversity and landscape heterogeneity across territorial spaces.As a super megacity in western China,Chengdu’s rapid urbanization has driven the evolution of U-S-R spaces,resulting in a sequential structure.To promote the high-quality spatial development of urban-rural region in a structured and efficient manner,it is essential to con-duct a scientific examination of the multidimensional interconnection within the U-S-RTZ framework.By proposing a novel identifica-tion method of U-S-RTZ and taking Chengdu,China as a case study,grounded in a blender of natural and humanistic factors,this study quantitatively delineated and explored the spatial evolutions of U-S-RTZ and stated the optimization orientation and sustainable devel-opment strategies of the production-living-ecological spaces along the U-S-R gradients.The results show that:1)it is suitable for the quantitative analysis of U-S-RTZ by established three-dimensional identification system in this study.2)In 1990-2020,the urban-sub-urban transition zones(U-STZ)in Chengdu have continuously undergone a substantial increase,and the scale of the suburban-rural transition zones(S-RTZ)has continued to expand slightly,while the space of rural-ecological transition zones(R-ETZ)has noticeably compressed.3)The landuse dynamics within U-S-RTZ has gradually increased in 1990-2020.The main direction of landuse transition was from farmland to construction land or woodlands,with the expansion of construction land being the most significant.4)R-ETZ primarily focus on ecological functions,and there is a trade-off relationship between the production-ecological function within the S-RTZ,and in the U-STZ,production-living composite functions are prioritized.This study emphasizes the importance of elastic planning and precise governance within the U-S-RTZ in a rapid urbanization region,particularly highlighting the role of suburbs as landscape corridors and service hubs in urban-rural integration.It elucidates to the practical implications for achieving high-quality development of integrated U-S-R territorial spaces.

Mechanisms of fracture propagation from multi-cluster using a phase field based HMD coupling model in fractured reservoir

Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recognized as significant factors affecting hydraulic fracture(HF)propagation,the combined influence of these factors remains poorly understood.To address this knowledge gap,a novel coupled hydromechanical-damage(HMD)model based on the phase field method is developed to investigate the propagation of multi-cluster HFs in fractured reservoirs.The comprehensive energy functional and control functions are established,while incorporating dynamic fluid distribution between multiple perforation clusters and refined changes in rock mechanical parameters during hydraulic fracturing.The HMD coupled multi-cluster HF propagation model investigates various scenarios,including single HF and single NF,reservoir heterogeneity,single HF and NF clusters,and multi-cluster HFs with NF clusters.The results show that the HMD coupling model can accurately capture the impact of approach angle(θ),stress difference and cementation strength on the interaction of HF and NF.The criterion of the open and cross zones is not fixed.The NF angle(a)is not a decisive parameter to discriminate the interaction.According to the relationship between approach angle(θ)and NF angle(a),the contact relationship of HF can be divided into three categories(θ=a,θ<a,andθ>a).The connected NF can increase the complexity of HF by inducing it to form branch fracture,resulting in a fractal dimension of HF as high as2.1280 at angles of±45°.Inter-fracture interference from the heel to the toe of HF shows the phenomenon of no,strong and weak interference.Interestingly,under the influence of NFs,distant HFs from the injection can become dominant fractures.However,as a gradually increases,inter-fracture stress interference becomes the primary factor influencing HF propagation,gradually superseding the dominance of NF induced fractures.

Formation mechanism of fault accommodation zones under combined stress in graben basin:Implications from geomechanical modeling

A fault accommodation zone is a type of structure that is defined as regulating displacement and strain between faults structure.Increasing numbers of fault accommodation zones are being identified in graben basins,indicating the potential exploration target and petroleum accumulation areas.This study aims to analyze the formation mechanism and development of fault accommodation zones under combined stress by a numerical simulation method considering geomechanical modeling.Using three-dimensional(3-D)seismic interpretation and fractal dimension method,exampled with the Dongxin fault zone,the fault activity and fault combination pattern were conducted to quantitatively characterize the activity difference in fault accommodation zones.Combined with mechanical experiment test,a geomehcanical model was established for fault accommodation zones in a graben basin.Integrating the paleostress numerical simulations and structural physical simulation experiment,the developmental characteristics and genetic mechanism of fault accommodation zones were summarized.Influenced by multi movements and combined stresses,three significant tectonic evolution stages of the Dongxing Fault Zone(DXFZ)were distinguished:During the E_(s)^(3)sedimentary period,the large difference in the stress,strain,and rupture distribution in various faults were significant,and this stage was the key generation period for the prototype of the DXFZ,including the FAZ between large-scale faults.During the E_(s)^(2)sedimentary period,the EW-trending symmetric with opposite dipping normal faults and the NE-SW trending faults with large scale were furtherly developed.The junction area of two secondary normal faults were prone to be ruptured,performing significant period for inheriting and developing characteristics of fault accommodation zones.During the Es1 sedimentary period,the high-order faults in the DXFZ exhibited the obvious fault depressions and strike-slip activity,and the fault accommodation zones were furtherly inherited and developed.This stage was the molded and formative period of the FAZ,the low-order faults,and the depression in the DXFZ.

Twenty-Three Years of Insect Pollinator-Dependent Crop Studies in Agro-Ecological Zones of Cameroon (2011-2020)

Cultivation of pollinator-dependent crops has expanded globally, increasing our reliance on insect pollination. This essential ecosystem service is provided by a wide range of managed and wild pollinators whose abundance and diversity are thought to be in decline, threatening sustainable food production. In Cameroon, several studies on pollinator-dependent crops carried out in different agro-ecological zones (AEZ) have been published in national and international journals, in order to present the importance and impact of flowering insects on fruit and seed yields of plant species. We proposed to produce a review article highlighting the different flowering insects and their importance for different plants according to AEZ, without however focusing on the quality of the journal (predator or non-predator) and how the different insects were identified (scientific names given in the publications). Thus, from 1997 to 2020, we collected 116 published papers from which only 26 were kept for this review. The results show that Hymenoptera, including the Apidae, followed by Megachilidae, are the most excellent pollinators of plant species in Cameroon, and they are present in different agro-ecological zones. The majority of publications focused on bees, particularly the honeybee Apis mellifera.

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.

Disparities in Prevalence and Determinants of Hypertension amongst Bamiléké Adults Residing in Two Different Agroecological Zones of Cameroon

Background: Cardiovascular diseases such as hypertension (HTN) are one of the main causes of death in Cameroon. This study aimed at assessing prevalence disparities and determinants of hypertension amongst Bamiléké adults residing in two different agroecological zones of Cameroon. Methods: A cross-sectional and descriptive survey was conducted among Bamiléké population living in the Highlands zone (Western region) and in the Monomodal Rainforest zone (Littoral region) of Cameroon from August 2016 to August 2017. Participants (962) were aged at least 20 years old. Data on sociodemographic, hemodynamic, anthropometric, and biochemical parameters and lifestyle of the participants were collected. Results: Results obtained revealed that 34.2% were hypertensive and those residing in the highland zone were more affected than those living in the monomodal rainforest zone (44.5% vs 22.9%). The different subtypes of HTN (Isolated systolic hypertension (14.1%), isolated diastolic hypertension (7.2%) and Systo-diastolic hypertension (23.3%)) were also more prevalent in the Highlands Zone. The most prevalent stage of HTN was pre-HTN (31.5%). However, people living in the monomodal rainforest zone were more affected by pre-HTN compared to Bamiléké living in the highland zone (33.6% vs. 29.6%). Results also showed that high consumption (≥ 3 times/week) of carbohydrate- and fat-rich foods, ageing, obesity, and marital status were associated with high blood pressure in both agroecological zones. Besides, secondary education (OR = 0.68;95% CI: 0.42 - 0.99) in the Highlands Zone and high (≥3 times/week) vegetable consumption (OR = 0.66;95% CI: 0.44 - 0.98) in the Monomodal Rainforest Zone had a protective effect on elevated blood pressure of population. Conclusion: There is a disparity in the prevalence of hypertension and some of its determinants among Bamiléké adults residing in different agroecological zones. This work highlights the need to advocate for local and ethno-cultural health policies to prevent, diagnose and manage hypertension.

A new elastoplastic model for bolt-grouted fractured rock

Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Identifying flow regime in the aquifer of fractured rock system in Germi Chai Basin,Iran

Considering the importance of fractured rock aquifers in the hydrogeologic process,this research aimed to analyze the flow regime,internal degree of karstification,and estimate storage volume in fractured rock aquifers of the Germi Chai Basin in northwest Iran,which is attributed to its active tectonics,erosion,and the lithological diversity.Given the geological setting,the hypothesis is that this basin is characterized by a high degree of karstification and diffuse or intermediate flow regime leading to variation in discharge flow rate.The hydrodynamic and hadrochemical analysis was conducted on 9 well distributed springs across the basin from 2019 to 2020.The maximum flow rate in most of the springs appeared in the early wet season despite their different levels of fluctuations on the monthly discharge time series.Analyzing the spring recession curve form revealed an aquifer containing multiple micro-regimes withαrecession coefficients and a degree of karstification ranging between 0.001 to 0.06 and 0.55 to 2.61,respectively.These findings indicated a dominant diffuse and intermediate flow system resulting from the development of a high density of fractures in this area.The electrical conductivity of the spring changes inversely proportional to the change in flow discharge,indicating the reasonable hydrological response of the aquifer to rainfall events.Hydrograph analysis revealed that the delay time of spring discharge after rainfall events mostly varies between 10 to 30 days.The total dynamic storage volume of the spring for a given period(2019-2020)was estimated to be approximately 1324 million cubic meters reflecting the long-term drainage potential and high perdurability of dynamic storage.Estimating the maximum and minimum ratio revealed that the springs recharging system in Germi Chai Basin comes under the slow aquifers category.This finding provides valuable insight into the hydrogeological properties of fractured rock aquifers contributing to effective water management strategy.

Evaluation of excavation damaged zones(EDZs)in Horonobe Underground Research Laboratory(URL)

Excavation of underground caverns,such as mountain tunnels and energy-storage caverns,may cause the damages to the surrounding rock as a result of the stress redistribution.In this influenced zone,new cracks and discontinuities are created or propagate in the rock mass.Therefore,it is effective to measure and evaluate the acoustic emission(AE)events generated by the rocks,which is a small elastic vibration,and permeability change.The authors have developed a long-term measurement device that incorporates an optical AE(O-AE)sensor,an optical pore pressure sensor,and an optical temperature sensor in a single multi-optical measurement probe(MOP).Japan Atomic Energy Agency has been conducting R&D activities to enhance the reliability of high-level radioactive waste(HLW)deep geological disposal technology.In a high-level radioactive disposal project,one of the challenges is the development of methods for long-term monitoring of rock mass behavior.Therefore,in January 2014,the long-term measurements of the hydro-mechanical behavior of the rock mass were launched using the developed MOP in the vicinity of 350 m below the surface at the Horonobe Underground Research Center.The measurement results show that AEs occur frequently up to 1.5 m from the wall during excavation.In addition,hydraulic conductivity increased by 2e4 orders of magnitude.Elastoplastic analysis revealed that the hydraulic behavior of the rock mass affected the pore pressure fluctuations and caused micro-fractures.Based on this,a conceptual model is developed to represent the excavation damaged zone(EDZ),which contributes to the safe geological disposal of radioactive waste.

Exploring pore-scale production characteristics of oil shale after CO_(2) huff‘n’puff in fractured shale with varied permeability

Recent studies have indicated that the injection of carbon dioxide(CO_(2))can lead to increased oil recovery in fractured shale reservoirs following natural depletion.Despite advancements in understanding mass exchange processes in subsurface formations,there remains a knowledge gap concerning the disparities in these processes between the matrix and fractures at the pore scale in formations with varying permeability.This study aims to experimentally investigate the CO_(2) diffusion behaviors and in situ oil recovery through a CO_(2) huff‘n’puff process in the Jimsar shale oil reservoir.To achieve this,we designed three matrix-fracture models with different permeabilities(0.074 mD,0.170 mD,and 0.466 mD)and experimented at 30 MPa and 91℃.The oil concentration in both the matrix and fracture was monitored using a low-field nuclear magnetic resonance(LF-NMR)technique to quantify in situ oil recovery and elucidate mass-exchange behaviors.The results showed that after three cycles of CO_(2) huff‘n’puff,the total recovery degree increased from 30.28%to 34.95%as the matrix permeability of the core samples increased from 0.074 to 0.466 mD,indicating a positive correlation between CO_(2) extraction efficiency and matrix permeability.Under similar fracture conditions,the increase in matrix permeability further promoted CO_(2) extraction efficiency during CO_(2) huff‘n’puff.Specifically,the increase in matrix permeability of the core had the greatest effect on the extraction of the first-cycle injection in large pores,which increased from 16.42%to 36.64%.The findings from our research provide valuable insights into the CO_(2) huff‘n’puff effects in different pore sizes following fracturing under varying permeability conditions,shedding light on the mechanisms of CO_(2)-enhanced oil recovery in fractured shale oil reservoirs.

Experimental Investigation on Condensate Revaporization During Gas Injection Development in Fractured Gas Condensate Reservoirs

The gas field in the Bohai Bay Basin is a fractured metamorphic buried-hill reservoir with dual-media characteristics. The retrograde vaporization mechanism observed in this type of gas condensate reservoir differs significantly from that observed in sand gas condensate reservoirs. However, studies on improving the recovery of fractured gas condensate reservoirs are limited;thus, the impact of retrograde vaporization on condensate within fractured metamorphic buried-hill reservoirs remains unclear. To address this gap, a series of gas injection experiments are conducted in pressure-volume-temperature(PVT) cells and long-cores to investigate the retrograde vaporization effect of condensate using different gas injection media in fractured gas condensate reservoirs. We analyze the variation in condensate volume, gas-to-oil ratio, and condensate recovery during gas injection and examine the influence of various gas injection media(CO_(2), N_(2), and dry gas) under different reservoir properties and varying gas injection times. The results demonstrate that the exchange of components between injected gas and condensate significantly influences condensate retrograde vaporization in the formation. Compared with dry gas injection and N_(2) injection,CO_(2) injection exhibits a superior retrograde vaporization effect. At a CO_(2) injection volume of 1 PV, the percentage shrinkage volume of condensate is 13.82%. Additionally, at the maximum retrograde condensation pressure, CO_(2) injection can increase the recovery of condensate by 22.4%. However, the condensate recovery is notably lower in fractured gas condensate reservoirs than in homogeneous reservoirs, owing to the creation of dominant gas channeling by fractures, which leads to decreased condensate recovery. Regarding gas injection timing, the effect of gas injection at reservoir pressure on improving condensate recovery is superior to that of gas injection at the maximum retrograde condensation pressure. This research provides valuable guidance for designing gas injection development plans and dynamic tracking adjustments for fractured gas condensate reservoirs.

Evaluation of Well Spacing for Primary Development of Fractured Horizontal Wells in Tight Sandstone Gas Reservoirs

Methods for horizontal well spacing calculation in tight gas reservoirs are still adversely affected by the complexity of related control factors,such as strong reservoir heterogeneity and seepage mechanisms.In this study,the stress sensitivity and threshold pressure gradient of various types of reservoirs are quantitatively evaluated through reservoir seepage experiments.On the basis of these experiments,a numerical simulation model(based on the special seepage mechanism)and an inverse dynamic reserve algorithm(with different equivalent drainage areas)were developed.The well spacing ranges of Classes I,II,and III wells in the Q gas field are determined to be 802–1,000,600–662,and 285–400 m,respectively,with their average ranges as 901,631,and 342.5 m,respectively.By considering both the pairs of parallel well groups and series well groups as examples,the reliability of the calculation results is verified.It is shown that the combination of the two models can reduce errors and provide accurate results.

Characteristics of transfer zones under the influence of preexisting faults and regional stress transformation:Wenchang A subsag,Zhujiang River Mouth Basin,northern South China Sea

A transfer zone in rift basins preserves important information on regional tectonic evolution and plays significant roles in hydrocarbon accumulation.Based on the systematic analysis of 3D seismic data and hydrocarbon accumulation conditions,the geometry,kinematics,and reservoir control of a large synthetic overlapping transfer zone in the south of the Wenchang A subsag in the Zhujiang(Pearl)River Mouth basin were investigated.Results indicate that the development and evolution of the transfer zone was controlled by the interaction between pre-existing faults and regional stress transformation.The intense rifting of the main faults of the transfer zone controlled the development of source rocks and faultcontrolled slope break paleogeomorphology.The strike-slip overprint since the Oligocene is conducive to the formation of a large-scale fault-anticline trap,and the secondary faults in the transfer zone contribute to the hydrocarbon transportation.The conjugate intersection area of the NE-and NW-trending faults offers more opportunity for hydrocarbon migration and accumulation.

A multi-mechanism numerical simulation model for CO_(2)-EOR and storage in fractured shale oil reservoirs

Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.