Micromechanism and mathematical model of stress sensitivity in tight reservoirs of binary granular medium

Research on reservoir rock stress sensitivity has traditionally focused on unary granular structures,neglecting the binary nature of real reservoirs,especially tight reservoirs.Understanding the stresssensitive behavior and mathematical characterization of binary granular media remains a challenging task.In this study,we conducted online-NMR experiments to investigate the permeability and porosity evolution as well as stress-sensitive control mechanisms in tight sandy conglomerate samples.The results revealed stress sensitivity coefficients between 0.042 and 0.098 and permeability damage rates ranging from 65.6%to 90.9%,with an average pore compression coefficient of 0.0168—0.0208 MPa 1.Pore-scale compression occurred in three stages:filling,compression,and compaction,with matrix pores playing a dominant role in pore compression.The stress sensitivity of binary granular media was found to be influenced by the support structure and particle properties.High stress sensitivity was associated with small fine particle size,high fines content,high uniformity coefficient of particle size,high plastic deformation,and low Young's modulus.Matrix-supported samples exhibited a high irreversible permeability damage rate(average=74.2%)and stress sensitivity coefficients(average=0.089),with pore spaces more slit-like.In contrast,grain-supported samples showed low stress sensitivity coefficients(average=0.021)at high stress stages.Based on the experiments,we developed a mathematical model for stress sensitivity in binary granular media,considering binary granular properties and nested interactions using Hertz contact deformation and Poiseuille theory.By describing the change in activity content of fines under stress,we characterized the non-stationary state of compressive deformation in the binary granular structure and classified the reservoir into three categories.The model was applied for production prediction using actual data from the Mahu reservoir in China,showing that the energy retention rates of support-dominated,fill-dominated,and matrix-controlled reservoirs should be higher than 70.1%,88%,and 90.2%,respectively.

A continuous and long-term in-situ stress measuring method based on fiber optic. Part I: Theory of inverse differential strain analysis

A method for in-situ stress measurement via fiber optics was proposed. The method utilizes the relationship between rock mass elastic parameters and in-situ stress. The approach offers the advantage of long-term stress measurements with high spatial resolution and frequency, significantly enhancing the ability to measure in-situ stress. The sensing casing, spirally wrapped with fiber optic, is cemented into the formation to establish a formation sensing nerve. Injecting fluid into the casing generates strain disturbance, establishing the relationship between rock mass properties and treatment pressure.Moreover, an optimization algorithm is established to invert the elastic parameters of formation via fiber optic strains. In the first part of this paper series, we established the theoretical basis for the inverse differential strain analysis method for in-situ stress measurement, which was subsequently verified using an analytical model. This paper is the fundamental basis for the inverse differential strain analysis method.

Support design method for deep soft-rock tunnels in non-hydrostatic high in-situ stress field

Due to the long-term plate tectonic movements in southwestern China,the in-situ stress field in deep formations is complex.When passing through deep soft-rock mass under non-hydrostatic high in-situ stress field,tunnels will suffer serious asymmetric deformation.There is no available support design method for tunnels under such a situation in existing studies to clarify the support time and support stiffness.This study first analyzed the mechanical behavior of tunnels in non-hydrostatic in-situ stress field and derived the theoretical equations of the ground squeezing curve(GSC)and ground loosening curve(GLC).Then,based on the convergence confinement theory,the support design method of deep soft-rock tunnels under non-hydrostatic high in-situ stress field was established considering both squeezing and loosening pressures.In addition,this method can provide the clear support time and support stiffness of the second layer of initial support.The proposed design method was applied to the Wanhe tunnel of the China-Laos railway in China.Monitoring data indicated that the optimal support scheme had a good effect on controlling the tunnel deformation in non-hydrostatic high in-situ stress field.Field applications showed that the secondary lining could be constructed properly.

Performance of water-coupled charge blasting under different in-situ stresses

Water-coupled charge blasting is a promising technique to efficiently break rock masses.In this study,numerical models of double boreholes with water-coupled charge are established using LS-DYNA and are calibrated by the tests of rock masses subjected to explosion loads to examine its performance.The crack levels of rock mass induced by water-coupled charge blasting and air-coupled charge blasting are first compared.It is found that water-coupled charge blasting is more appropriate to fracture deep rock mass than air-coupled charge blasting.In addition,the effects of rock properties,water-coupled charge coefficients,and borehole connection angles on the performance of water-coupled charge blasting are investigated.The results show that rock properties and water-coupled charge coefficients can greatly influence the crack and fragmentation levels of rock mass induced by water-coupled charge blasting under uniform and non-uniform in-situ stresses.However,changing borehole-connection angles can only affect crack and fragmentation levels of rock mass under non-uniform in-situ stresses but barely affect those under uniform in-situ stresses.A formula is finally proposed by considering the above-mentioned factors to provide the design suggestion of water-coupled charge blasting to fracture rock mass with different in-situ stresses.

Influences of clean fracturing fluid viscosity and horizontal in-situ stress difference on hydraulic fracture propagation and morphology in coal seam

The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal samples for experiments, and clean fracturing fluid samples were prepared using CTAB surfactant. A series of hydraulic fracturing tests were conducted with an in-house developed triaxial hydraulic fracturing simulator and the fracturing process was monitored with an acoustic emission instrument to analyze the influences of fracturing fluid viscosity and horizontal in-situ stress difference on coal fracture propagation. The results show that the number of branched fractures decreased, the fracture pattern became simpler, the fractures width increased obviously, and the distribution of AE event points was concentrated with the increase of the fracturing fluid viscosity or the horizontal in-situ stress difference. The acoustic emission energy decreases with the increase of fracturing fluid viscosity and increases with the increase of horizontal in situ stress difference. The low viscosity clean fracturing fluid has strong elasticity and is easy to be compressed into the tip of fractures, resulting in complex fractures. The high viscosity clean fracturing fluids are the opposite. Our experimental results provide a reference and scientific basis for the design and optimization of field hydraulic fracturing parameters.

Impact of well placement and flow rate on production efficiency and stress field in the fractured geothermal reservoirs

Geothermal energy has gained wide attention as a renewable alternative for mitigating greenhouse gas emissions.The advancements in enhanced geothermal system technology have enabled the exploitation of previously inaccessible geothermal resources.However,the extraction of geothermal energy from deep reservoirs poses many challenges due to high‐temperature and high‐geostress conditions.These factors can significantly impact the surrounding rock and its fracture formation.A comprehensive understanding of the thermal–hydraulic–mechanical(THM)coupling effect is crucial to the safe and efficient exploitation of geothermal resources.This study presented a THM coupling numerical model for the geothermal reservoir of the Yangbajing geothermal system.This proposed model investigated the geothermal exploitation performance and the stress distribution within the reservoir under various combinations of geothermal wells and mass flow rates.The geothermal system performance was evaluated by the criteria of outlet temperature and geothermal productivity.The results indicate that the longer distance between wells can increase the outlet temperature of production wells and improve extraction efficiency in the short term.In contrast,the shorter distance between wells can reduce the heat exchange area and thus mitigate the impact on the reservoir stress.A larger mass flow rate is conducive to the production capacity enhancement of the geothermal system and,in turn causes a wider range of stress disturbance.These findings provide valuable insights into the optimization of geothermal energy extraction while considering reservoir safety and long‐term sustainability.This study deepens the understanding of the THM coupling effects in geothermal systems and provides an efficient and environmentally friendly strategy for a geothermal energy system.

In-situ stress of coal reservoirs in the Zhengzhuang area of the southern Qinshui Basin and its effects on coalbed methane development

In-situ stress is a critical factor influencing the permeability of coal reservoirs and the production capacity of coalbed methane(CBM)wells.Accurate prediction of in-situ stress and investigation of its influence on coal reservoir permeability and production capacity are significant for CBM development.This study investigated the CBM development zone in the Zhengzhuang area of the Qinshui Basin.According to the low mechanical strength of coal reservoirs,this study derived a calculation model of the in-situ stress of coal reservoirs based on the multi-loop hydraulic fracturing method and analyzed the impacts of initial fractures on the calculated results.Moreover,by combining the data such as the in-situ stress,permeability,and drainage and recovery data of CBM wells,this study revealed the spatial distribution patterns of the current in-situ stress of the coal reservoirs and discussed the impacts of the insitu stress on the permeability and production capacity.The results are as follows.(1)Under given fracturing pressure,longer initial fractures are associated with higher calculated maximum horizontal principal stress values.Therefore,ignoring the effects of the initial fractures will cause the calculated values of the in-situ stress to be less than the actual values.(2)As the burial depth increases,the fracturing pressure,closure pressure,and the maximum and minimum horizontal principal stress of the coal reservoirs in the Zhengzhuang area constantly increase.The average gradients of the maximum and minimum horizontal principal stress are 3.17 MPa/100 m and 2.05 MPa/100 m,respectively.(3)Coal reservoir permeability is significantly controlled by the magnitude and state of the current in-situ stress.The coal reservoir permeability decreases exponentially with an increase in the effective principal stress.Moreover,a low lateral pressure coefficient(less than 1)is associated with minor horizontal compressive effects and high coal reservoir permeability.(4)Under similar conditions,such as resource endowments,CBM well capacity is higher in primary structural coal regions with moderate paleotectonic stress modification,low current in-situ stress,and lateral pressure coefficient of less than 1.

Numerical analysis on mechanical difference of sandstone under in-situ stress,pore pressure preserved environment at depth

Deep in-situ rock mechanics considers the influence of the in-situ environment on mechanical properties,differentiating it from traditional rock mechanics.To investigate the effect of in-situ stress,pore pressure preserved environment on the mechanical difference of sandstone,four tests are numerically modeled by COMSOL:conventional triaxial test,conventional pore pressure test,in-situ stress restoration and reconstruction test,and in-situ pore pressure-preserved test(not yet realized in the laboratory).The in-situ stress restoration parameter is introduced to characterize the recovery effect of in-situ stress on elastic modulus and heterogeneous distribution of sandstone at different depths.A random function and nonuniform pore pressure coefficient are employed to describe the non-uniform distribution of pore pressure in the in-situ environment.Numerical results are compared with existing experimental data to validate the models and calibrate the numerical parameters.By extracting mechanical parameters from numerical cores,the stress-strain curves of the four tests under different depths,in-situ stress and pore pressure are compared.The influence of non-uniform pore pressure coefficient and depth on the peak strength of sandstone is analyzed.The results show a strong linear relationship between the in-situ stress restoration parameter and depth,effectively characterizing the enhanced effect of stress restoration and reconstruction methods on the elastic modulus of conventional cores at different depths.The in-situ pore pressurepreserved test exhibits lower peak stress and peak strain compared to the other three tests,and sandstone subjected to non-uniform pore pressure is more prone to plastic damage and failure.Moreover,the influence of non-uniform pore pressure on peak strength gradually diminished with increasing depth.

Stress sensitivity of carbonate gas reservoirs and its microscopic mechanism

In order to evaluate the stress sensitivity of carbonate reservoirs,a series of rock stress sensitivity tests were carried out under in-situ formation temperature and stress condition.Based on the calibration of capillary pressure curve,the variable fractal dimension was introduced to establish the conversion formula between relaxation time and pore size.By using the nuclear magnetic resonance(NMR)method,the pore volume loss caused by stress sensitivity within different scales of pore throat was quantitatively analyzed,and the microscopic mechanism of stress sensitivity of carbonate gas reservoirs was clarified.The results show that fractures can significantly affect the stress sensitivity of carbonate reservoirs.With the increase of initial permeability,the stress sensitivity coefficient decreases and then increases for porous reservoirs,but increases monotonously for fractured-porous reservoirs.The pore volume loss caused by stress sensitivity mainly occurs for mesopores(0.02–0.50μm),contributing more than 50%of the total volume loss.Single high-angle fracture contributes 9.6%of the stress sensitivity and 15.7%of the irreversible damage.The microscopic mechanism of the stress sensitivity of carbonate gas reservoirs can be concluded as fracture closure,elastic contraction of pores and plastic deformation of rock skeleton.

Influences of lithology on in-situ stress field in low permeability reservoirs in Bonan Oilfield, Bohai Bay Basin, China

The differences of rock mechanical properties were analyzed based on triaxial compression test in low permeability reservoirs of the Bonan Oilfield. Through the analysis of reservoir mechanics, the influence mechanisms of different mechanical properties of rocks on reservoir in-situ stress were studied. By means of stress ellipse and finite element simulation, the influence rules of different mechanical properties of rocks on in-situ stress field were discussed. For the low permeability reservoirs of the Bonan Oilfield, the coarser rock has a larger Young’s modulus value and a lower Poisson’s ratio. The rock mechanical parameters and stress-strain relationship of sandstone facies and mudstone facies are different. Different rocks have different mechanical properties, which cause extra stress at the lithological contact interface, and the existence of extra stress affects the reservoir in-situ stress. Without considering the influence of structural features on the in-situ stress field, the reservoir in-situ stress is controlled by the magnitude of extra stress and the angle between lithological contact surface and boundary stress.

Analysis on method for effective in-situ stress measurement in hot dry rock reservoir

With the rapid increase of energy demand and the increasingly highlighted environmental problems, clean, safe and widely distributed geothermal resources have become a hot spot for renewable resources development. The state of in-situ stress is a major control parameter for multiple links including well location, fracture inspiration and reservoir assessment, so how to determine the accurate state of in-situ stress in the deep thermal reservoir becomes a core problem drawing widely attention and urgent to be solved. Based on features of hot dry rock reservoir in terms of temperature and pressure and the comparison analysis, this article proposes the method of Anelastic Strain Recovery(ASR) as an effective method for determining the state of in-situ stress in the area with HDR resources distributed and explains the availability of ASR method by application examples.

Comparison between double caliper,imaging logs,and array sonic log for determining the in-situ stress direction:A case study from the ultra-deep fractured tight sandstone reservoirs,the Cretaceous Bashijiqike Formation in Keshen8 region of Kuqa depress

The tight sandstone in the Tarim Basin has the characteristics of large burial depth and development of nature fractures due to concentrated in-situ stress. Identifying the present-day in-situ stress orientation is important in hydrocarbon exploration and development, but also a key scientific question in understanding naturally fractured reservoirs. This paper presents a case study where we integrate various methods using wireline and image-log data, to identify present-day in-situ stress direction of ultra-deep fractured tight sandstone reservoirs, in the Kuqa depression. We discuss the formation mechanism of the elliptical borehole, compares the advantages and applicable conditions of the double caliper method,resistivity image logs and array sonic logs method. The well borehole diameter is measured orthogonally,then the ellipse is fitted, and the in-situ stress orientation is identified by the azimuth of the short-axis borehole, but it fails in the borehole expansion section, the fracture development section and the borehole collapse section. The micro-resistivity image logs method reveals the borehole breakouts azimuth, and also the strike of induced fractures, which are used to determine the orientation of in-situ stress. However, under water-based mud conditions, it’s hard to distinguish natural fractures from induced fractures by image logs. Under oil-based mud conditions, the induced fractures are difficult to identify due to the compromised image quality. As for the sonic log, shear waves will split when passing through an anisotropic formation, shear waves will split during propagation, and the azimuth of fast shear waves is consistent with the orientation of in-situ stress. However, it is usually affected by the anisotropy caused by the excessively fast rotation of the well log tools, so that the azimuth of fast shear wave cannot effectively reflect the orientation of the in-situ stress. Based on comprehensive assessment and comparison, in this paper we propose a method integrating various logging data to identify the orientation of in-situ stress. Among various types of logging data, the breakouts azimuth identified by image logs is proved to be the most credible in identifying the orientation of in-situ stress, while using the direction of induced fractures under water-based mud conditions is also viable. However, the azimuth of the fast shear wave is consistent with the orientation of maximum in-situ stress only when the rotation speed of the logging tool is low. The caliper method can be used as a reference for verifying the other two methods. Using this integrated method to study the orientation of in-situ stress in the Keshen8 trap, the results show that faults are an important factor affecting the direction of in-situ stress, while multi-level faults will produce superimposed effects that cause the current direction of in-situ stress to change.

Development status and research recommendations for thermal extraction technology in deep hot dry rock reservoirs

Based on a comprehensive review of domestic and foreign literature, this article discusses the technical difficulties and development status of enhanced geothermal system(EGS) concerning the thermal energy extraction of deep hot dry rock(HDR) reservoirs and proposes suggestions for the research focus of numerical simulation of HDR reservoir stimulation. Additionally, it summarizes the existing methods and mainstream working fluids for HDR reservoir stimulation. The article emphasizes the significance of factors such as well location, production well depth, artificial fracture orientation, and complexity in optimizing the thermal production efficiency of the EGS. Furthermore, this article delves into a detailed discussion on the influence of fracture spacing, fracture permeability,fracture length, fluid injection rate, and injected fluid temperature on the performance of the EGS. In light of the thermo-hydro-mechanical coupling challenges associated with high-temperature reservoirs, it is suggested that future research efforts should focus on investigating the impact of thermo-induced stresses on the stability of the artificial fracture network within the EGS during long-term(>30 years) circulation of hot and cold fluids.

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.

Rate transient analysis methods for water-producing gas wells in tight reservoirs with mobile water

Tight gas reservoirs with mobile water exhibit multi-phase flow and high stress sensitivity.Accurately analyzing the reservoir and well parameters using conventional single-phase rate transient analysis methods proves challenging.This study introduces novel rate transient analysis methods incorporating evaluation processes based on the conventional flowing material balance method and the Blasingame type-curve method to examine fractured gas wells producing water.By positing a gas-water two-phase equivalent homogenous phase that considers characteristics of mobile water,gas,and high stress sensitivity,the conventional single-phase rate transient analysis methods can be applied by integrating the phase's characteristics and defining the phase's normalized parameters and material balance pseudotime.The rate transient analysis methods based on the equivalent homogenous phase can be used to quantitatively assess the parameters of wells and gas reservoirs,such as original gas-in-place,fracture half-length,reservoir permeability,and well drainage radius.This facilitates the analysis of production dynamics of fractured wells and well-controlled areas,subsequently aiding in locating residual gas and guiding the configuration of well patterns.The specific evaluation processes are detailed.Additionally,a numerical simulation mechanism model was constructed to verify the reliability of the developed methods.The methods introduced have been successfully implemented in field water-producing gas wells within tight gas reservoirs containing mobile water.

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.

Numerical study on gas production via a horizontal well from hydrate reservoirs with different slope angles in the South China Sea

It is important to study the effect of hydrate production on the physical and mechanical properties of low-permeability clayey–silty reservoirs for the largescale exploitation of hydrate reservoirs in the South China Sea.In this study,a multiphysical-field coupling model,combined with actual exploration drilling data and the mechanical experimental data of hydrate cores in the laboratory,was established to investigate the physical and mechanical properties of low-permeability reservoirs with different slope angles during 5-year hydrate production by the depressurization method via a horizontal well.The result shows that the permeability of reservoirs severely affects gas production rate,and the maximum gas production amount of a 20-m-long horizontal well can reach186.8 m3/day during the 5-year hydrate production.Reservoirs with smaller slope angles show higher gas production rates.The depressurization propagation and hydrate dissociation mainly develop along the direction parallel to the slope.Besides,the mean effective stress of reservoirs is concentrated in the near-wellbore area with the on-going hydrate production,and gradually decreases with the increase of the slope angle.Different from the effective stress distribution law,the total reservoir settlement amount first decreases and then increases with the increase of the slope angle.The maximum settlement of reservoirs with a 0°slope angle is up to 3.4 m,and the displacement in the near-wellbore area is as high as2.2 m after 5 years of hydrate production.It is concluded that the pore pressure drop region of low-permeability reservoirs in the South China Sea is limited,and various slope angles further lead to differences in effective stress and strain of reservoirs during hydrate production,resulting in severe uneven settlement of reservoirs.

Influence of the Earth Stress on the Development of Gas Reservoirs

Gas reservoirs are located kilometers deep beneath the earth's surface under great earth stresses, including the overburden stress and the horizontal stress. After a well is drilled, the stress condition around the well bore will be changed. During the development, a pressure funnel forms around the hole, with the rock stress redistributed. In this paper, the influence of the earth stress on the gas reservoir development, including the output, the period of the steady output, the recovery and ratio, is researched thoroughly with the theory of reservoir seepage dynamics and clarified with the calculation method. The research shows that the earth stress produces impacts on the development of a gas reservoir when the stress is great.

Stress initialization methods for dynamic numerical simulation of rock mass with high in-situ stress

In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.

Time-sensitivity of the Kaiser effect of acoustic emission in limestone and its application to measurements of in-situ stress

Measuring in-situ stress by using the Kaiser effect in rocks has such advantages as timeefficiency, low cost and little limitation, but the precision of the method is dependent on rock properties and delay time of the measurement. In this paper, experiments on the Kaiser effect in limestones were performed, and it was found that the limestones had good ability to retain a memory of their recent stress history and high time-sensitivity. The longer the experiment was delayed from the extraction of the stone, the larger the Felicity ratio was. As the Felicity ratio approached l, significant Kaiser effect was observed. In-situ stress should be determined by the limestone measurements when the delay time was 40-120 days. Finally, the in-situ stress in a limestone formation could be successfully measured in practice.