Experimental research on the effect caused by effective stress change in injection production capacity of carbonate reservoir

Abstract

The effective stress of the reservoir rock in the underground gas storage (UGS) changes with the multi cycle injection and production, which changes the physical characteristics of the reservoir and affects the injection production capacity of the UGS. Taking carbonate fracture-pore type gas storage as the research object, combined with the actual production of UGS, establish a set of experimental evaluation methods based on production characteristics; Core damage evaluator, overburden pressure porosimeter and scanning electron microscope are used to carry out experimental research on the impact of effective stress change on injection production capacity, quantitative analysis of reservoir internal characteristics, multi cycle injection and production and effect of effective stress action time on injection-production capacity. The results show, the experimental method established based on the actual production characteristics of the UGS, simulation of multi cycle injection and production by forward and reverse displacement, simulation of effective stress change under constant confining pressure and variable internal pressure, the experimental evaluation of the injection production capacity of the simulation multi cycle injection and production of the UGS is realized; The change of effective stress makes the reservoir seepage characteristics worse and the reservoir space smaller, injection-production coefficient (Szc) is positively correlated with the internal characteristic value of the reservoir, the larger the internal characteristic value of the reservoir, the higher Szc, and the better the injection production capacity; The injection production capacity of the UGS reservoir deteriorates with the increase of the injection production cycle, in the 4th injection production cycle, the injection production capacity of the UGS changes from moderately strong to moderately weak; The damage of effective stress to the injection production capacity of UGS mainly occurs in the initial stage of injection production, and the damage tends to be stable with the progress of injection production. The conclusion is that the change of effective stress in carbonate fracture-pore type gas storage makes the injection-production capacity of the reservoir worse, it is proposed to carry out reservoir reconstruction before the UGS is put into operation, and expand the capacity of the gas storage after the completion of the 3rd injection and production cycle, which can effectively offset the impact of the decline of injection production capacity of the UGS and play the role of reservoir protection, it provides a constructive idea for improving the operation efficiency of the UGS.

1 Introduction

In recent years, with the rapid growth of natural gas consumption, the contradiction between insufficient underground gas storage (UGS) capacity with safe and stable supply of natural gas has further intensified. As the “ballast stone” to ensure the safe and stable operation of the natural gas energy market, UGS plays an indispensable role in the whole natural gas industry chain (Ma et al., 2022). According to international experience, once the external dependence of natural gas exceeds 30%, the working gas volume of UGS needs to exceed 12% of the consumption, and China’s external dependence on natural gas will climb to 44.9% in 2021 (Ding et al., 2022). However, the gas storage and peak shaving facilities are far from meeting the aforementioned standards. Compared with other countries, China’s UGS construction is still at the primary level: the working gas volume of UGS is only 3%–4% of the national natural gas consumption, and the international average is 12%–15% (Yang et al., 2021). To improve the gas storage and peak shaving capacity of UGS, we should not only pay attention to the construction of the number of UGS, but also strengthen the excavation of the potential of the built UGS to ensure the efficient operation of the UGS (Jiang et al., 2021). UGS must meet the requirements of “gas injection, storage and production”, and must complete gas injection and production in a short time (generally, gas injection is completed once in 8 months and gas production is completed once in 4 months), forming an injection production cycle. The operation of the UGS is repeated between the upper and lower limit pressures, and the effective stress borne by the reservoir rock changes periodically. This change will inevitably lead to changes in the reservoir seepage characteristics (Quinn and MacDonald, 1992; Castelletto et al., 2013), which will affect the injection production capacity and storage capacity of the UGS, and is not conducive to the efficient operation of UGS. Therefore, it is of great significance to clarify the impact of periodic stress changes on the injection production capacity of UGS reservoirs, put forward injection and production measures for UGS, scientific formulation of injection-production plan.

At present, many important results and understandings have been obtained in the experimental study of the influence of stress sensitivity on the percolation characteristics of oil and gas reservoirs, but less is involved in the multi cycle injection and production of carbonate fracture-pore UGS. The lower the permeability of sandstone reservoir, the stronger the stress sensitivity (Wang et al., 2021; Xu et al., 2022); Strong stress sensitive rocks have obvious fracture characteristics, while weak stress sensitive rocks have porosity characteristics (Cuss et al., 2003; Engelder, 2014); The pressure-sensitive experiment of ultra-high pressure gas reservoir is carried out, and the binomial productivity equation considering the influence of pressure-sensitive effect is established (Engelder and Fischer, 1994); Carry out rock permeability and porosity change experiments under effective stress loading, describe the relationship between porosity, permeability and effective stress by establishing functions, spatial models, etc (Oughanem et al., 2013; Rahmati et al., 2013; Moosavi et al., 2014; Selvadurai and Glowacki, 2016). A method for evaluating the stress sensitivity is established, which mainly includes using the constant in the relationship between regression permeability and effective stress index to express the stress sensitivity (Bazin et al., 2008), stress sensitivity coefficient method (You et al., 2013), and a series of evaluation methods such as the industry standard flow test evaluation method of reservoir sensitivity (SY/T 5358-2010).

Multi cycle injection and production of sandstone gas storage will reduce the permeability of reservoir rocks, but with the increase of injection and production rounds, the permeability decline will decrease. The lower the permeability, the greater the irreversible permeability loss rate, and the stronger the periodic stress sensitivity (Bakhtiari et al., 2021; Wang et al., 2023). The operation of the UGS has the characteristics of multi cycle, strong injection and production, and alternating load, it is easy to induce reservoir damage such as stress sensitivity, particle migration and sand production. Through the optimization design of injection production parameters and storage capacity parameters, the safe and stable operation of the UGS can be realized (Dehghan et al., 2016; Verga, 2018; Honari and Seyedi Hosseininia, 2021). Multi cycle injection and production of UGS, adopting different injection and production measures, can not only meet production needs, but also alleviate reservoir damage caused by particle migration and improve the operation efficiency of UGS (Arfaee and Sola, 2014; Xinping et al., 2020; Subbiah et al., 2021).

In summary, the research on stress sensitivity at home and abroad mainly focuses on the seepage characteristics and evaluation methods of oil and gas reservoirs, mainly serving the oil and gas productivity. The research on multi cycle injection and production stress sensitivity of UGS is relatively few, and most of the UGS involved are sandstone reservoirs, and the research on the impact of effective stress changes of carbonate fracture-pore type gas storage on injection production capacity is rarely involved, and lack of corresponding experimental methods and evaluation indicators. For carbonate rock fracture-pore type gas storgae, natural gas is stored and transported in the fracture pore, and the periodic change of effective stress will inevitably lead to the change of reservoir seepage characteristics, which will affect the injection production capacity and storage capacity of the UGS. However, at present, no research has been carried out on the change of injection production capacity caused by the periodic rise and fall of effective stress under the conditions of multi cycle injection and production for carbonate fracture-pore type gas storage.

During the operation of the UGS, the periodic change of effective stress will change the physical properties of the reservoir, especially the permeability and porosity of the reservoir. The author selected carbonate fracture-pore reservoir rock samples of X gas storage, established an experimental method based on the injection production characteristics of the X UGS, and gave the indexes for evaluating the injection production capacity of the UGS; Carry out experiments on the influence of internal characteristic values of the reservoir, multi cycle injection and production and considering the effective stress action time on the injection production capacity of the reservoir. By detecting the changes of permeability and porosity of the rock sample during the experiment, and observing the changes of fracture wall characteristics of the rock sample before and after the experiment with a scanning electron microscope, analyze the impact of effective stress changes on the injection production capacity of carbonate fracture-pore type gas storage, and reveal the mechanism of action. Based on the understanding of experimental research, provide reservoir protection measures and suggestions, it is expected to provide experimental technical support for optimizing the injection-production system of carbonate fracture-pore gas reservoirs and improving the operation efficiency of UGS.

2 Experimental samples and equipments

2.1 Experimental samples

The experimental rock sample is from the core in the injection production wells of X UGS, and the reservoir is carbonate fracture-pore type. According to the analysis of logging and coring data, the porosity of the reservoir is 0.39%–16.96%, and the permeability is 0.001–223.69 mD. The core is processed into a column sample with a diameter of 2.5 cm and a length of 5 cm (Figure 1), testing porosity and permeability, select core samples similar to reservoir physical properties as experimental rock samples.

FIGURE 1

2.2 Experimental equipments

Permeability is measured by core damage evaluation instrument, porosity is measured by overburden porosity-permeability instrument, and microstructure is observed by scanning electron microscope. The core damage evaluation instrument mainly consists of injection system, temperature control system, pressure control system, back pressure system, metering system, data acquisition and control system, the injection system can adjust the pressure by 0–50 MPa, the temperature control system room temperature ∼150°C, the maximum confining pressure of the pressure control system is 90 MPa, and the accuracy of the metering system is 0.001 g (Figure 2A). The overburden porosity-permeability instrument is mainly composed of pressure regulating system, sampling system, automatic control system, data acquisition and analysis system, etc., the porosity measurement range: 0.01%–40%, the permeability measurement range: 0.00005–15,000 mD, the pressure sensor: 0–250 psi, 0–10,000 psi (Figure 2B).

FIGURE 2

3 Experimental methods

The operation of UGS is a process in which the overlying formation pressure remains unchanged, the reservoir pore pressure increases during gas injection, and the reservoir pore pressure decreases during gas production, that is the operation is repeated between the upper and lower pressure of the UGS. According to the principle of similarity and similarity, it is converted into laboratory experiment, using the effective stress change mode of constant confining pressure and variable pore pressure, using confining pressure to simulate the overlying pressure of the UGS, the experimental injection pressure simulates the pressure in the pores of the gas reservoir, and the injection pressure changes to simulate the change of the operating pressure during the injection production of the gas reservoir.

3.1 Experimental parameters

For the pertinence and accuracy of experimental research, the experimental temperature 82 C (X UGS reservoir temperature), confining pressure 71 MPa (X UGS reservoir overlying formation pressure). Conduct triaxial rock mechanics experiments on experimental rock samples, under the experimental confining pressure of 71 MPa, the average compressive strength is 314.75 MPa, which means that under the experimental conditions of adding confining pressure, the experimental rock samples undergo elastic deformation when the pressure is less than 314.75 MPa; when the pressure reaches 314.75 MPa, plastic deformation occurs to breakage (Figure 3). The critical pressure for displacement deviation of the experimental rock sample is 314.75 MPa. Injection pressure is 13.20 MPa–28.06 MPa (X UGS operation of upper and lower limit pressure). The injection pressure is refined and decomposed in the experiment, the injection pressure with 3 MPa as the pressure interval from the lower limit 13 MPa to the upper limit 31 MPa, which are 13 MPa, 16 MPa, 19 MPa, 22 MPa, 25 MPa, 28 MPa, and 31 MPa. The effective stress is the difference between the confining pressure and the injection pressure, which are 58 MPa, 55 MPa, 52 MPa, 49 MPa, 46 MPa, 43 MPa, and 40 MPa (Table 1).

FIGURE 3

The experimental medium is 8% potassium chloride solution, in order to eliminate the effect of velocity sensitivity on the experimental results, the experimental displacement flow should be determined. The specific experimental steps are as follows: 1)Put the core sample into the core holder, increase the confining pressure, and let it stand for 4 h to eliminate the stress sensitivity; 2)Keep the confining pressure always higher than the inlet pressure 2 MPa, and test the permeability under the flow rates of 0.10 cm³/min, 0.25 cm³/min, 0.5 cm³/min, 0.75 cm³/min, 1.00 cm³/min, 1.5 cm³/min, and 2.0 cm³/min. When the permeability change rate is greater than 20%, it indicates that velocity sensitivity has occurred, and the flow corresponding previous point is the displacement flow. The permeability change rate can be expressed using Eq. (1).Where, D_w is the permeability change rate under different flow rates; K_i is the permeability of core under different flow rates, respectively (mD); K_w is the initial permeability of core, respectively (mD).

Test the displacement flow of different experimental samples, which shows that the flow rate when the permeability change rate is greater than 20% increases with the increase of initial permeability, but when the initial permeability of the sample is within a certain range, the flow rate when the permeability change rate is greater than 20% is a constant value (Table 2). The reservoir itself contains fine particles, particles migrate under the scouring of high-speed fluid, which may block the pore throat or accompany fluid discharge from the reservoir (Muecke, 1979; Gruesbeck and Collins, 1982; Civan, 2015; Li et al., 2017; Hou et al., 2015). When the initial permeability of the sample is within a certain range, its corresponding pore structure is consistent, and the fluid exhibits the same properties when passing through the pore. Therefore, different displacement flow rates are used for different experimental samples.

3.2 Experimental process

During the change of effective stress, test the permeability and porosity changes at different effective stress points. The experimental steps are as follows: 1)Put the core samples into the core holder of the overburden porosity-permeability instrument, test the porosity and permeability, select two core samples with similar porosity and permeability as a group, and select the cores representing the physical properties of different reservoirs as different groups; 2)A core sample of the same group is put into the core holder of the Core damage evaluation instrument, linearly load the confining pressure at a loading rate of 0.69 MPa/min, load the confining pressure to 71 MPa and keep it constant, and stabilize for 4 h; the displacement flow is used to pressurize to the injection pressure, maintain the pressure for 15 min, measure the permeability of this point; At the same time, another core sample of the same group is put into the core holder of the overburden porosity-permeability instrument, and the effective stress corresponding to the injection pressure is added as the confining pressure, maintain the pressure for 15 min, measure the pore of this point; 3)According to the injection pressure interval, increase the injection pressure in turn to the maximum, and test the permeability of core samples at different injection pressure points; The effective stress corresponding to the injection pressure is taken as the confining pressure to test the porosity of core samples at different effective stress points, this process simulates the gas injection process of UGS; 4)Change the displacement direction, according to the injection pressure interval, and then start from the maximum injection pressure, reduce the injection pressure to the minimum injection pressure by the displacement flow, and test the permeability of core samples at different injection pressure points; The effective stress corresponding to the injection pressure is taken as the confining pressure to test the porosity of core samples at different effective stress points, this process simulates the gas production process of UGS; 5)Repeat steps (2)–(4) to complete the experiment until to complete multi cycle injection and production; 6)Replace the core sample and repeat steps (2)–(5) until the experiment is completed. Finally, based on the change characteristics of permeability and porosity under different effective stress, define the effect of effective stress change on the injection production capacity of the reservoir.

3.3 Injection production capacity characterization

There is a correlation between porosity, permeability with effective stress. The pore volume and permeability of rocks decrease with the increase of the cycle of pressure reciprocating changes (Ranjith et al., 2013; Song et al., 2022). The natural gas in the UGS is stored and migrated in the fractures and pores of the reservoir, with the injection production of the UGS, the internal characteristics of the reservoir are also changing (Dake, 1983; Gobran et al., 1987; Economides and Nolte, 1989; Cosentino, 2001; Lucia et al., 2003; Jaeger et al., 2009; Lake, 2012; Ahmed, 2018) (Figure 4). In order to clarify the influence of the change of the internal characteristics of the reservoir on the injection production capacity of the reservoir, use the product of porosity and permeability to characterize the internal characteristic value of the reservoir. The ratio of the internal characteristic value of the reservoir at different injection and production stages of the UGS to the internal characteristic value of the reservoir before the construction UGS, called injection-production coefficient, which can be expressed using Eq. 2. Before the construction of the UGS, the maximum injection-production coefficient is 1, As the UGS is built and operated, the injection production capacity gradually changes. According to the operation status of the UGS, the injection-production coefficient is divided into four ranges, analyze the change of the injection-production coefficient with the periodic injection production of the UGS, indicates the injection production capacity of the UGS (Table 3).Where, S_zc is the injection-production coefficient of UGS at different injection-production stages, dimensionless; K_n is the permeability of UGS at different stages of injection and production, respectively (mD); Q_n is the porosity of UGS in different injection production stages, respectively (%); K_i is the initial permeability of UGS, respectively (mD); Q_i is the initial porosity of UGS, respectively (%).

FIGURE 4

4 Experimental result

4.1 Influence of reservoir internal characteristic values on reservoir injection production capacity

Select three groups of samples representing small, medium and large internal characteristic values of the UGS reservoir, namely, 1, 2 and 3 groups, to conduct single cycle injection and production, and test the change of injection-production coefficient under different effective stress (Figure 5).

FIGURE 5

For the same sample, under the same effective stress, Szc in the gas injection stage is greater than that in the gas production stage, which indicates that gas injection and production will damage the reservoir’s injection production capacity and is irreversible. Different samples have different Szc after injection and production, Group 3 (0.65) > Group 2 (0.55) > Group 1 (0.48), which shows that Szc is positively with the internal characteristic value of the reservoir, that is the larger the internal characteristic value of the reservoir after injection production, the higher Szc is, indicating that the higher the permeability and porosity of the reservoir, the better the injection production capacity; after the completion of injection and production, the variation value of Szc (Eq. 3) is different, Group 1 (0.08) > Group 2 (0.07) > Group 3 (0.05), indicating that the variation value of Szc is inversely correlated with the internal characteristic value of the reservoir, which shows that the greater the internal characteristic value of the reservoir, the smaller the change value of Szc, indicating that the higher the initial permeability and porosity of the reservoir, the less damage to the injection production capacity of the reservoir caused by gas injection and production. Therefore, reservoir reconstruction after completion of the injection production wells in UGS can improve the permeability and porosity of the reservoir and play a role in reservoir protection.Where, △S_zc is the variation value of injection-production coefficient; S_zc^’ is the injection-production coefficient in gas injection stage; S_zc^’’ is the injection production coefficient in gas production stage.

4.2 Influence of multi cycle injection and production on reservoir injection production capacity

Select 2nd group of internal characteristic values of representative UGS reservoir for multi cycle injection and production, and test the change of Szc under different effective stresses (Figure 6).

FIGURE 6

In the same injection production cycle, during gas injection reservoir pore pressure increases, while the overlying formation pressure remains unchanged, the increased pore pressure offsets part of the overlying formation pressure, and the effective stress acting on the reservoir decreases, making the reservoir internal fractures and pores space of the reservoir partially recovered, which is shown as Szc increases with the decrease of the effective stress; Contrary to gas injection during gas production, the reservoir pore pressure decreases, while the overlying formation pressure remains unchanged, and the effective stress acting on the reservoir increases, making the internal fractures and pores space compression of the reservoir smaller, which is shown as Szc decreases with the increase of effective stress. At the same effective stress point, Szc during gas injection is greater than Szc during gas production, and the difference increases with the increase of effective stress, indicating that the change of effective stress is irreversible to the damage of reservoir injection production capacity. In different injection production cycles, Szc decreases with the increase of injection production cycle, from 0.75 at the beginning of the first injection production cycle to 0.43 at the end of the fifth injection production cycle, but the decrease amplitude becomes smaller. After the third injection production cycle, the decrease value of injection-production coefficient tends to be stable, indicating that the damage to the reservoir injection production ability becomes weaker and stable due to the change of effective stress with the increase of injection production cycle. The injection-production coefficient of the 1st∼3rd injection production cycle is 0.75–0.51, and the injection production capacity of the reservoir is moderately strong; 0.49–0.42 after the fourth injection production cycle, the injection production capacity of the reservoir is moderately weak. Therefore, in order to maintain the injection production capacity of the UGS reservoir, the UGS capacity will be expanded after the completion of the 3rd injection and production cycle.

4.3 Effect of effective stress time on injection production capacity

Select 3 samples representing the internal characteristic value of the gas reservoir, the effective stress is constant, With the increase of the effective stress acting time, real time test of injection-production coefficient changes under different effective stresses (Figure 7).

FIGURE 7

The overburden pressure of the gas reservoir is 71 MPa, and the UGS operates repeatedly at the upper and lower pressure of 28.06 MPa–13.20 MPa (the effective stress is 42.94 MPa–57.8 MPa), so the effective stress is 42.94 MPa, 57.8 MPa and the average value is 50.37 MPa. Under the same effective stress, Szc shows three stages of rapidly decrease, stable and further decrease with the increase of action time. The rapid decrease of Szc occurred within 1.5 h, with a decrease range of 2%–5%; Szc occurs stably in 1.5 h–3.0 h; Szc decreased slightly after 3.0 h. For different effective stresses, the smaller the effective stress, the greater the Szc, under the same action time, the greater the effective stress, the smaller the Szc. It shows that the damage of effective stress to the injection production capacity of UGS reservoir occurs at the beginning of injection production, and the damage tends to be stable with the injection production.

5 Analysis of factors affecting injection production capacity

The change of effective stress will cause rock deformation, which mainly includes the deformation of fracture pore space, the deformation of skeleton particles and the deformation of skeleton particle arrangement (Gilbert et al., 2013; Liu et al., 2022). The deformation of skeleton particles belongs to elastic deformation, the effective stress can usually recover after unloading; The deformation due to the change of particle arrangement of the skeleton belongs to plastic deformation, and it is difficult to recover the effective stress after unloading (Sui et al., 2011).

5.1 Internal characteristic value of reservoir

According to the scanning electron microscope pictures of the experimental samples, it can be observed that the pore morphology in the sample is mainly circular pore, and at the same time, there are long and narrow natural fractures (Figure 8). Carbonate fracture-pore type reservoir rock is hard in compactness, and the space for skeleton particle deformation is very small, therefore, the effective stress first causes the deformation of the fracture pore space, with the increase of effective stress, skeleton particle deformation and skeleton particle arrangement change successively (Figure 9). When the effective stress is between 40 MPa and 49 MPa (injection production pressure 31 MPa–22 MPa), Szc is larger and changes rapidly, the reason is that the reservoir rock with the increase of the effective pressure has obvious compaction, and the porosity and permeability decrease rapidly, at this stage, Szc changes are mainly caused by the deformation of the fracture porosity space. When the effective stress changes between 58 MPa and 49 MPa (injection production pressure 13 MPa–22 MPa), Szc is relatively small and basically unchanged, indicating that the injection production capacity is relatively weak and stable, the reason is that the space available for compression in the rock sample is reduced, and the decrease range of porosity and permeability is reduced, at this stage, Szc changes are mainly caused by the deformation of skeleton particles and the deformation caused by the change of skeleton particle arrangement.

FIGURE 8

FIGURE 9

5.2 Multi cycle injection and production

Multi cycle injection and production reduces the injection production capacity of the UGS, but with the increase of injection and production cycle, the injection production capacity tends to be stable, indicating that the fractures and pores of the reservoir rock are reduced to a certain extent due to multi cycle effective stress, but still have injection production capacity. The scanning electron microscope results show that there are dense and a large number of microconvex bodies on the fracture and pore wall of the sample, when the fracture and pore are closed to a certain extent, these dense and large number of microconvex bodies contact the other wall of the fracture and pore, forming a supporting role for the fracture and pore (Figures 10, 11).

FIGURE 10

FIGURE 11

5.3 Effective stress action time

Figure 4 shows, Szc can be divided into three stages of rapid decrease, stability and redecrease with the increase of effective stress acting time. At the beginning, fractures and pores are closed to a certain extent under the effect of effective stress, resulting in a rapid reduction of Szc. After the effective stress is applied for a period of time, Szc is relatively stable due to the support of the microconvex body on the fracture and pore wall (Figure 12A). As the wall microconvex bodies, especially those in point contact, are prone to stress concentration, as the effective stress action time continues to increase, the microconvex body will break (Figure 12B), and the fractures and pores will be closed to a certain extent again, eventually causing Szc to decrease again (Figure 12C).

FIGURE 12

6 Conclusions and suggestions

References

• 1

  AhmedT. (2018). Reservoir engineering handbook. Cambridge, MA: Gulf professional publishing.

  • Google Scholar

• 2

  ArfaeeM. I. R.SolaB. S. (2014). Investigating the effect of fracture–matrix interaction in underground gas storage process at condensate naturally fractured reservoirs. J. Nat. Gas Sci. Eng.19, 161–174. 10.1016/j.jngse.2014.05.007

  • CrossRef
  • Google Scholar

• 3

  BakhtiariM.ShadS.ZivarD.RazaghiN. (2021). Coupled hydromechanical analysis of underground gas storage at Sarajeh field, Qom formation, Iran. J. Nat. Gas. Sci. Eng.92, 103996. 10.1016/j.jngse.2021.103996

  • CrossRef
  • Google Scholar

• 4

  BazinB.BekriS.VizikaO.HerzhaftB.AubryE. (2008). “Fracturing in tight-gas reservoirs: Application of SCAL methods to investigate formation-damage mechanisms,” in SPE International Symposium and Exhibition on Formation Damage Control (OnePetro). 10.2118/112460-MS

  • CrossRef
  • Google Scholar

• 5

  CastellettoN.FerronatoM.GambolatiG.JannaC.MarzoratiD.TeatiniP. (2013). Can natural fluid pore pressure Be safely exceeded in storing gas underground?Eng. Geol.153, 35–44. 10.1016/j.enggeo.2012.11.008

  • CrossRef
  • Google Scholar

• 6

  CivanF. (2015). Reservoir formation damage. Waltham, MA: Gulf professional publishing.

  • Google Scholar

• 7

  CosentinoL. (2001). Integrated reservoir studies. Paris, France: Institut Français du Pétrole Publications.

  • Google Scholar

• 8

  CussR. J.RutterE. H.HollowayR. F. (2003). Experimental observations of the mechanics of borehole failure in porous sandstone. Int. J. Rock Mech. Min. Sci.40 (5), 747–761. 10.1016/S1365-1609(03)00068-6

  • CrossRef
  • Google Scholar

• 9

  DakeL. P. (1983). Fundamentals of reservoir engineering. Amsterdam, The Netherlands: Elsevier.

  • Google Scholar

• 10

  DehghanA. N.GoshtasbiK.AhangariK.JinY. (2016). Mechanism of fracture initiation and propagation using a tri-axial hydraulic fracturing test system in naturally fractured reservoirs. Eur. J. Environ. Civ. Eng.20, 560–585. 10.1080/19648189.2015.1056384

  • CrossRef
  • Google Scholar

• 11

  DingG.DingY.LiY.TangL.WuZ.WanyanQ.et al (2022). Prospects of underground gas storage in China under the strategy of carbon neutrality. Oil Gas Storage Transp.41 (1), 1–7. 10.6047/j.issn.1000-8241.2022.01.001

  • CrossRef
  • Google Scholar

• 12

  EconomidesM. J.NolteK. G. (1989). Reservoir stimulation, Vol. 2. Englewood Cliffs, NJ: Prentice Hall, 6–10.

  • Google Scholar

• 13

  EngelderT.FischerM. P. (1994). Influence of poroelastic behavior on the magnitude of minimum horizontal stress, S h in overpressured parts of sedimentary basins. Geology22 (10), 949–952. 10.1130/0091-7613(1994)022<0949:IOPBOT>2.3.CO;2

  • CrossRef
  • Google Scholar

• 14

  EngelderT. (2014). Stress regimes in the lithospher, Vol. 151. New Jersey: Princeton University Press.

  • Google Scholar

• 15

  GilbertG.BackéR.PuspitasariZ.PallikathekathilJ.ManeyB.DewhurstD. (2013). Model. geomechanics gas storage a case study Iona gas field. Aust. J. Greenh. Gas. Control13, 138–148. 10.1016/j.ijggc.2012.12.009

  • CrossRef
  • Google Scholar

• 16

  GobranB. D.BrighamW. E.RameyH. J. (1987). Absolute permeability as a function of confining pressure, pore pressure, and temperature. SPE Form. Eval.2 (01), 77–84. 10.2118/10156-PA

  • CrossRef
  • Google Scholar

• 17

  GruesbeckC.CollinsR. E. (1982). Entrainment and deposition of fine particles in porous media. Soc. Petroleum Eng. J.22 (06), 847–856. 10.2118/8430-pa

  • CrossRef
  • Google Scholar

• 18

  HonariS.Seyedi HosseininiaE. (2021). Particulate modeling of sand production using coupled DEM-LBM. Energies14 (4), 906. 10.3390/en14040906

  • CrossRef
  • Google Scholar

• 19

  HouD.ZhangH.LiG. (2015). Study of pressure propagation characteristics to determine stress sensitivity of low-permeability gas reservoir. Chem. Technol. Fuels Oils51, 181–189. 10.1007/s10553-015-0591-z

  • CrossRef
  • Google Scholar

• 20

  JaegerJ. C.CookN. G.ZimmermanR. (2009). Fundamentals of rock mechanics. Malden, MA: Blackwell Publishing.

  • Google Scholar

• 21

  JiangT.WangJ.WangZ.ZhouD.LiY.WangS.et al (2021). A novel protein encoded by circMAPK1 inhibits progression of gastric cancer by suppressing activation of MAPK signaling. Nat. Gas. Ind.41 (9), 66–74. 10.1186/s12943-021-01358-y

  • CrossRef
  • Google Scholar

• 22

  LakeL. (Editor) (2012). Reservoir characterization (Orlando, Florida: Academic Press, Inc.).

  • Google Scholar

• 23

  LiX. P.CaoL. N.LuoC.ZhangB.ZhangJ. Q.TanX. H. (2017). Characteristics of transient production rate performance of horizontal well in fractured tight gas reservoirs with stress-sensitivity effect. J. Petroleum Sci. Eng.158, 92–106. 10.1016/j.petrol.2017.08.041

  • CrossRef
  • Google Scholar

• 24

  LiuJ.MeiL.DingW.XuK.YangH.LiuY. (2022). Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs. GSA Bull.135, 678–688. 10.1130/B36358.1

  • CrossRef
  • Google Scholar

• 25

  LuciaF. J.KeransC.JenningsJ. W. (2003). Carbonate reservoir characterization. J. Petroleum Technol.55 (06), 70–72. 10.2118/82071-ms

  • CrossRef
  • Google Scholar

• 26

  MaX.ZhengD.WeiG.DingG.ZhengS. (2022). Development directions of major scientific theories and technologies for underground gas storage. Nat. Gas. Ind.42 (5), 93–99. 10.3787/j.issn.1000-0976.2022.05.010

  • CrossRef
  • Google Scholar

• 27

  MoosaviS. A.GoshtasbiK.KazemzadehE.BakhtiariH. A.EsfahaniM. R.ValiJ. (2014). Relationship between porosity and permeability with stress using pore volume compressibility characteristic of reservoir rocks. Arabian J. Geosciences7 (1), 231–239. 10.1007/s12517-012-0760-x

  • CrossRef
  • Google Scholar

• 28

  MueckeT. W. (1979). Formation fines and factors controlling their movement in porous media. J. petroleum Technol.31 (02), 144–150. 10.2118/7007-pa

  • CrossRef
  • Google Scholar

• 29

  OughanemR.YoussefS.BazinB.MaireE.VizikaO. (2013). “Pore-scale to core-scale aspects of capillary desaturation curves using CT-scan imaging,” in IOR 2013-17th European Symposium on Improved Oil Recovery (European Association of Geoscientists & Engineers), 342. 10.3997/2214-4609.20142615

  • CrossRef
  • Google Scholar

• 30

  QuinnD. F.MacDonaldJ. A. (1992). Natural gas storage. Carbon30 (7), 1097–1103. 10.1016/0008-6223(92)90141-I

  • CrossRef
  • Google Scholar

• 31

  RahmatiH.JafarpourM.AzadbakhtS.NouriA.VaziriH.ChanD.et al (2013). Review of sand production prediction models. J. Petroleum Eng.2013, 1–16. 10.1155/2013/864981

  • CrossRef
  • Google Scholar

• 32

  RanjithP. G.PereraM. S. A.PereraW. K. G.WuB.ChoiS. K. (2013). Effective parameters for sand production in unconsolidated formations: An experimental study. J. Petroleum Sci. Eng.105, 34–42. 10.1016/j.petrol.2013.03.023

  • CrossRef
  • Google Scholar

• 33

  SelvaduraiA. P. S.GlowackiA. (2017). Stress-induced permeability alterations in an argillaceous limestone. Rock Mech. Rock Eng.50 (5), 1079–1096. 10.1007/s00603-016-1153-3

  • CrossRef
  • Google Scholar

• 34

  SongR.ZhangP.TianX.HuangF.LiZ.LiuJ. (2022). Study on critical drawdown pressure of sanding for wellbore of underground gas storage in a depleted gas reservoir. Energies15 (16), 5913. 10.3390/en15165913

  • CrossRef
  • Google Scholar

• 35

  SubbiahS. K.SamsuriA.Mohamad-HusseinA.JaafarM. Z.ChenY. R.KumarR. R. (2021). Root cause of sand production and methodologies for prediction. Petroleum7 (3), 263–271. 10.1016/j.petlm.2020.09.007

  • CrossRef
  • Google Scholar

• 36

  SuiY.LinT.LiuX.DongC.ChengW.ZhangG.et al (2011). Numerical simulation for the effect of cyclic stress on microstructure and mechanical properties of reservoir rocks during operations of gas storage. J. China Univ. Pet. Nat. Sci.44 (2), 78–84. 10.1016/j.proeng.2011.05.040

  • CrossRef
  • Google Scholar

• 37

  VergaF. (2018). What’s conventional and what’s special in a reservoir study for underground gas storage. Energies11 (5), 1245. 10.3390/en11051245

  • CrossRef
  • Google Scholar

• 38

  WangJ.FengX.WanyanQ.ZhaoQ.PeiZ.ZieG.et al (2021). Hysteresis effect of three-phase fluids in the high-intensity injection–production process of sandstone underground gas storages. Energy242, 123058. 10.1016/j.energy.2021.123058

  • CrossRef
  • Google Scholar

• 39

  WangQ-y.ShangL.XinC-y.SunY-c.GaoG-l.WangM.et al (2023). Application of the in-situ stress testing technology for the design of operating pressure of underground gas storage reservoir. Front. Earth Sci.10, 1002676. 10.3389/feart.2022.1002676

  • CrossRef
  • Google Scholar

• 40

  XinpingG. A. O.JunliangP. E. N. G.HuifenH. A. N.YuchaoZ. H. O. U.LimingL. I.YanF. E. N. G. (2020). Experimental research on the damage caused by particle migration in Xiang guo Si gas storage. Drill. Prod. Technol.43 (3), 46. 10.3969/J.ISSN.1006-768X.2020.03.14

  • CrossRef
  • Google Scholar

• 41

  XuD.ChenS.ChenJ.XueJ.YangH. (2022). Study on the imbibition damage mechanisms of fracturing fluid for the whole fracturing process in a tight sandstone gas reservoir. Energies15 (12), 4463. 10.3390/en15124463

  • CrossRef
  • Google Scholar

• 42

  YangY.ChenJ.WangL.ZhangH.NiJ.ZhangY.et al (2021). Analysis on the development prospects of China's gas storage business and the development path of operation models. PETROLEUM NEW ENERGY33 (3), 11–14. 10.3969/j.issn.2097-0021.2021.03.003

  • CrossRef
  • Google Scholar

• 43

  YouL.XueK.KangY.LiaoY.KongL. (2013). Pore structure and limit pressure of gas slippage effect in tight sandstone. Sci. World J.2013, 1–7. 10.1155/2013/572140

  • CrossRef
  • Google Scholar