Stability Of Caverns Research Articles

Comparative Study of Temperature and Pressure Variation Patterns in Hydrogen and Natural Gas Storage in Salt Cavern

Clarifying the distribution of temperature and pressure in the wellbore and cavern during hydrogen injection and extraction is crucial for quantitatively assessing cavern stability and wellbore integrity. This paper establishes an integrated flow and heat transfer model for the cavern and wellbore during hydrogen injection and withdrawal, analyzing the variations in temperature and pressure in both the wellbore and the cavern. The temperature and pressure parameters of hydrogen and natural gas within the chamber and wellbore were compared. The specific conclusions are as follows. (1) Under identical injection and withdrawal conditions, the temperature of hydrogen in the chamber was 10 °C higher than that of natural gas, and 16 °C higher in the wellbore. The pressure of hydrogen in the chamber was 2.9 MPa greater than that of natural gas, and 2.6 MPa higher in the wellbore. (2) A comparative analysis was conducted on the impact of surrounding rock’s horizontal and numerical distance on temperature during hydrogen and natural gas injection processes. As the distance from the cavity increases, from 5 to 15 m, the temperature fluctuation in the surrounding rock diminishes progressively, with the temperature effect in the hydrogen storage chamber extending to at least 10 m. (3) The influence of rock thermal conductivity parameters on temperature during the processes of hydrogen injection and natural gas extraction is also compared. The better the thermal conductivity, the deeper the thermal effects penetrate the rock layers, with the specific heat capacity having the most significant impact.

Effect of a gunite lining on the stability of salt caverns in bedded formations for hydrogen storage

In this paper, a novel methodology to improve the structural behaviour and durability of bedded salt caverns for hydrogen storage is proposed. The suggested technology consists in applying a gunite lining over the whole cavern surface by means of a pneumatic air blowing. This continuous flow projected at high speed onto the cavern surface produces a self-compacted gunite facing. The effect of adding a gunite lining on the structural performance of a bedded salt cavern is analysed. The nonlinear finite element method is then applied to solve the inherent nonlinear solid mechanics problem. It is shown that reinforcing salt caverns with a gunite lining improves the cavern behaviour, specially in the case of caverns located within a horizontal salt stratum confined between other strata. When dealing with these salt formations, the stratum geometry limits the cavern height. In general, wider caverns that exhibit larger displacements are required to achieve a huge storage capacity. Moreover, these caverns can be located at a considerable depth, where the temperature has an important role on the creep deformation process. Therefore, the proposition of stabilization techniques is required to ensure both their structural integrity and the ground subsidence. This technology allows the safe exploitation of salt caverns located at bedded-type stratified horizontal salt formations.

A multi-step calibration strategy for reliable parameter determination of salt rock mechanics constitutive models

The storage of renewable hydrogen in salt caverns requires fast injection and production rates to cope with the imbalance between energy production and consumption. This raises concerns about the mechanical stability of salt caverns under such operational conditions. The use of appropriate constitutive models for salt mechanics is an important step in investigating this issue, therefore many constitutive models with several parameters have been presented in the literature. However, a robust calibration strategy to reliably determine which model and parameter set represents the given rock, based on stress–strain data sets, remains an unsolved challenge. In this paper, for the first time in the community, we present a multi-step strategy to determine a single parameter set based on many deformation data sets for salt rocks. Towards this end, we first develop a comprehensive constitutive model able to capture all relevant nonlinear deformation physics of transient, reverse, and steady-state creep. The determination of the single set of representative material parameters is then achieved by framing the calibration process as an optimization problem, for which the global Particle Swarm Optimization algorithm is employed. To allow for dynamic data integration, a multi-step calibration strategy is developed for a situation where experiments are included one at a time, as they become available. Additionally, due to the existing mild heterogeneity in the experimental rock samples, our optimization strategy is made flexible to allow for this slight heterogeneity. The devised optimization strategy, based on the multi-physics comprehensive constitutive modeling framework, results in a single set of representative material properties of all the deformation data sets. As a rigorous mathematical analysis for the presented method and the lack of relevant experimental data sets, we consider a wide range of synthetic experimental data sets, inspired by the existing sparse relevant data in the literature. The results of our performance analyses show that the proposed calibration strategy is robust. Moreover, the results become increasingly more accurate as more data sets become available.

The influence of steep faults at various avoidance distances on the stability of hard rock caverns: Physical model experiments and DEM simulations

The influence of avoidance distances from steep faults on the stability of hard rock caverns is pivotal in deep engineering, yet research in this area remains scant. This study constructed 10 physical models of granite caverns, incorporating artificial cement-filled faults at various avoidance distances ranging from 0 to 2.4 cavern diameters. Biaxial compression experiments were conducted, utilizing acoustic emission and digital image correlation techniques for monitoring. Coupled with discrete element method simulations, the analysis explored the influence mechanisms of fault-cavern distance on the failure modes of hard rock caverns. The results reveal a positive correlation between the peak strength of specimens and fault-cavern distance. As the distance approaches 0.6 cavern diameters, the fault significantly influences the displacement field, crack number, distribution, and crack initiation location within the specimens. Furthermore, even at 1.2 cavern diameters, the fault still exhibits a notable influence on the stress field, external load transfer paths, cavern failure modes, and the size distribution of spalling debris. These findings can guide the selection of sites for deep hard rock caverns, ensuring a minimum fault-cavern distance of 1.2 cavern diameters to mitigate the risks posed by steep faults to such caverns.

Research of interlayer dip angle effect on stability of salt cavern energy and carbon storages in bedded salt rock

In order to clean and decarbonize energy, large-scale underground energy storage in salt caverns can be utilized. The construction and safety evaluations for these caverns are crucial, especially for salt caverns with significant interlayer dip angles. Aimed at the bedded salt rocks for energy and carbon storage, this study focuses on the impact of different interlayer dip angles on the stability of underground energy and carbon storage caverns in bedded salt rock formations in China. Four salt caverns with varying interlayer dip angles of 0–30° were established for stability comparison. Numerical simulations reveal that caverns with high interlayer dip angles demonstrate low volume shrinkage, general deformation resistance and small equivalent strain, while exhibiting a larger volume of plastic zones compared to caverns with lower interlayer dip angles. It is important to consider the deformation of the cavern wall rock in the down-dip direction for caverns with high interlayer dip angles. The presence of inclined interlayers affects the symmetry of the contour distribution of deformation, safety factor and equivalent strain. Furthermore, recommendations are provided for the construction of salt caverns in salt rock formations with high interlayer dip angles.

True triaxial modeling test of high-sidewall underground caverns subjected to dynamic disturbances

Seismicity resulting from the near- or in-field fault activation significantly affects the stability of large-scale underground caverns that are operating under high-stress conditions. A comprehensive scientific assessment of the operational safety of such caverns requires an in-depth understanding of the response characteristics of the rock mass subjected to dynamic disturbances. To address this issue, we conducted true triaxial modeling tests and dynamic numerical simulations on large underground caverns to investigate the impact of static stress levels, dynamic load parameters, and input directions on the response characteristics of the surrounding rock mass. The findings reveal that: (1) When subjected to identical incident stress waves and static loads, the surrounding rock mass exhibits the greatest stress response during horizontal incidence. When the incident direction is fixed, the mechanical response is more pronounced at the cavern wall parallel to the direction of dynamic loading. (2) A high initial static stress level specifically enhances the impact of dynamic loading. (3) The response of the surrounding rock mass is directly linked to the amplitude of the incident stress wave. High amplitude results in tensile damage in regions experiencing tensile stress concentration under static loading and shear damage in regions experiencing compressive stress concentration. These results have significant implications for the evaluation and prevention of dynamic disasters in the surrounding rock of underground caverns experiencing dynamic disturbances.

Long-term deformation of rock salt under creep–fatigue stress loading paths: Modeling and prediction

Rock salt, due to its water solubility, low permeability, high plasticity, and damage self-healing ability, is one of the best candidate rock types for underground energy storage. Utilizing salt caves to construct compressed air energy storage (CAES) facilities can effectively enhance the utilization of renewable energy. Due to the need for peak shaving, the surrounding rock of the salt cavern will undergo discontinuous cyclic loading with varying gas injection rates and pressures, namely, alternating creep–fatigue loading. Considering the actual peak-shaving cycle of a CAES plant, long-term creep–fatigue tests of rock salt with different loading cycles and stress levels were conducted. The results indicate that in long-term creep–fatigue tests for rock salt, the lower the loading stress rate is, the greater the deformation of the rock salt. The variation in the stress limit has a greater effect on creep than on fatigue loading and unloading. The deformation rate of rock salt is influenced by alterations in the stress state. Based on the test results, according to the Norton creep model, a new creep–fatigue constitutive model for rock salt was established by defining a state variable that characterizes the level of rock hardening and introducing unloading as well as crack factors. This model can accurately describe the impact of historical loading and unloading processes on the viscoplastic mechanical characteristics of rock salt. The rock salt creep–fatigue test results were used to verify the constitutive model. A comparison of the fitting curve of the different stress loading paths with the test curve reveals good consistency, indicating that the model comprehensively considers the effects of time, load, and state on rock salt creep–fatigue, effectively describing the viscoplastic deformation characteristics of rock salt under different stress paths. These research findings provide important guidance for ensuring the stability of salt caverns used for CAES.

Mechanical Behavior of Sediment-Type High-Impurity Salt Cavern Gas Storage during Long-Term Operation

With the development of salt cavern gas storage technology, the construction of large-scale salt cavern gas storage using sediment voids is expected to solve the problems of low effective volume formation rate and poor construction economy of high-impurity salt mines. At present, there are few studies on the long-term operational mechanical behavior of salt cavern gas storage under the influence of sediment accumulation. The present paper studies the influence of sediment height, particle gradation, and operating pressure on the stability of salt caverns by constructing a coupling model of sediment particle discontinuous medium and surrounding rock continuous medium. The continuous–discontinuous coupling algorithm is suitable for analyzing the influence of sediment height and particle gradation on the creep shrinkage of salt caverns. The increase in sediment height slows down the creep shrinkage of the cavern bottom, which strengthens the restraining effect on the surrounding rock of the cavern. As a result, the position of the maximum displacement of the surrounding rock moves to the upper part of the cavern. The sediment particle gradation has little effect on the cavern volume shrinkage rate. The greater the coarse particle content, the smaller the cavern volume shrinkage rate. The greater the operating pressure, the more conducive to maintaining the stability of the cavern. This situation slows down the upward movement of the sediment accumulation and increases the gas storage space in the upper part of the cavern. The obtained results can provide a reference for evaluating the long-term operational stability of sediment-type high-impurity salt cavern gas storage.

Probabilistic evaluation method for the stability of large underground cavern considering the uncertainty of rock mass mechanical parameters: A case study of Baihetan underground powerhouse project

For deep large underground cavern projects under complex geological conditions, the determination of rock mass mechanical parameters and stability analysis is fraught with a great deal of uncertainty. This paper proposes a set of methods for estimating rock mass mechanical parameters and probabilistic stability assessment suitable for construction characteristics of large underground caverns. On the one hand, according to the Hoek-Brown criterion and Monte Carlo simulation, a dynamic estimation method for rock mass mechanical parameters is proposed by using Rock Mass Rating (RMR), the uniaxial compressive strength (UCS) of intact rock, and the material constant (mb) as input parameters. On the other hand, considering the uncertainty of rock mass mechanical parameters, a probabilistic evaluation method applicable to the unloading response characteristics of rock mass (including displacement and fracturing zones) around the large cavern is put forward, combining the point estimation method and numerical simulation. The proposed methods were applied to an underground powerhouse, with an excavation span of 34 m and height of 88.7 m, at the Baihetan hydropower station in the southwest of China. Based on extensive field investigations and tests, the probability distributions of key rock mass mechanical parameters were obtained. Furthermore, the high sensitivity of input parameters in the H-B criterion for estimating rock mass mechanical parameters were revealed. Through dynamic simulation, the probability distributions of the displacement and fracturing zones in surrounding rock during the layer-by-layer excavation of the cavern were presented. Comprehensive on-site tests showed that the simulated results were in basic agreement with the actual unloading response behavior during the cavern excavation, validating the correctness and applicability of the proposed methods. The study provides a relatively comprehensive approach for estimating rock mass mechanical parameters and stability evaluation in similar underground engineering projects, and also has guiding significance for predicting and preventing engineering rock mass hazards.

Analysis and Experimental Study on the Stability of Large-Span Caverns’ Surrounding Rock Based on the Progressive Collapse Mechanism

The collapse failure of rock surrounding caverns involves a progressive collapse process. Based on the nonlinear Hoek–Brown failure criterion and the upper limit theorem, the whole process curve of the progressive collapse of the surrounding rock of a large-span cavern is outlined in this paper. The progressive collapse process of the surrounding rock of the large-span cavern is experimentally studied using an independently developed visualized large-span-cavern geomechanical model test device with variable angles. The results show that, through theoretical calculation and model tests, the surrounding rock at the top of the large-span cavern undergoes three collapses. Under the condition of rock mass and the shape of the cavern, the larger the span of the cavern, the more times the surrounding rock collapses; with the increase in surrounding rock pressure, the first collapse occurs in the middle part of the arch roof. When the overlying load reaches a certain level, the arch foot becomes the weakest part, and the rock undergoes shear failure along the arch foot, gradually extending upwards, accompanied by multiple collapses, forming a progressive collapse process. The theoretical calculation results of this paper are basically consistent with the scope of the model test, and the research results can provide a basis for the construction and support design of the large-span cavern.

Study on significance of pillar dimensions on stability of caverns in rocks with high overburden

The present study aims to examine the influence of pillar widths on the stability of caverns. Case study considering two caverns viz. Powerhouse Cavern (PHC) and Transformerhall Cavern (TC), in a major hydro-electric project in the eastern Himalayas is considered. 2D and 3D numerical analysis was carried out for w/B ratios 1.5 and 3.0 respectively. Primary aspects like major principal stress and development of plastic zone were investigated for the two pillar widths. An optimum pillar width was observed that resulted in reduced stress acting along the cavern periphery, a better stress distribution, and no overlap of plastic zones between the caverns. Further, the optimum pillar width resulted in a better stress-redistribution with progress of excavation and the in-situ stress became constant at an earlier stage of excavation. Observations from comparative analysis revealed that a pillar width nearly equal to the largest dimension or twice the width of the larger of the caverns in the group resulted in a better stability and hence can be considered as the optimum width. Furthermore, the analysis suggests that along the pillar width, maximum stress was observed at mid-height, and it is more in the vicinity of the face of the caverns.

Microfabrics impact on rock salt geomechanical diversity – an experimental investigation

Geomechanical studies of rock salt serve as the foundation for designing and ensuring the long-term stability of storage caverns. Research conducted in recent years has shown that rock salt exhibit highly variable mechanical parameters. Even macroscopically similar rock salt samples from the same formation have very different mechanical parameters in standard strength tests. In this paper, the effect of microfabrics on the mechanical behaviour of rock salt was studied using a uniaxial compression device with simultaneous observation under a microscope. The tests were carried out on rock salt samples, previously tested on standard uniaxial compression (UCS). Tested samples showed different strength values despite macroscopic similarity and similar anhydrite content. The tests revealed that the microfabrics, such as fluid inclusions and anhydrite crystals distributed along the halite grain boundaries, play a significant role in the strength of a rock salt sample. Four types of grain boundaries were recognized in both group of samples: (1) dominated by anhydrite crystals with a small number of fluid inclusions (FIs), (2) with FIs prevailing over anhydrite crystals, (3) only with FIs, and (4) without anhydrite crystals or FIs. The experimental uniaxial compression tests in microscale (mUCS) showed that the damage process is initiated by FIs activity at the halite grain boundaries and consists of three main stages. The samples dominated by anhydrite crystals at halite grain boundaries were characterised by the highest mUCS. Conversely, if FIs prevail at the grain boundaries the lowest mUCS values were registered. Therefore, anhydrite crystals located along halite grain boundaries influenced the mechanical behaviour and deformation process. On the contrary, anhydrite crystals dispersed within halite grains had no impact on failure and mechanical parameters. The experiment showed that apart from impurity content, the distribution of these impurities within halite grains should be examined during geomechanical evaluation of rock salt samples. Thus, the overall strength of rock salt samples determined in standard test results from the local variations of petrological features in micro- and macro scale.

Plastic zone and surrounding strain analysis in a lined cavern for compressed air energy storage subjected to cyclic high internal gas pressure and temperature

In the context of compressed air energy storage (CAES) technology, the strategic excavation of hard rock caverns is pivotal for overcoming geological challenges and advancing CAES deployment. This study adopts a circular profile for the cavern and investigates the impact of variables such as cavern diameter and depth. Considering the influence of cyclic internal gas pressure, temperature fluctuations, and potential leakage on the cavern, a thermo-hydro-mechanical coupled numerical model is developed using Comsol software. This model, grounded in classic elastoplastic constitutive relations, aligns with the conditions of a proposed project. It facilitates the computation of the plastic zone and surrounding strain for both the cavern lining and adjacent rock, post-excavation and throughout a typical ‘charging-storage-discharging-storage’ operational cycle. Through detailed analysis, the study delineates the variations and characteristics of the plastic zone and surrounding strain, thereby providing a solid foundation for a comprehensive examination of cavern stability, and the design of linings and sealing mechanisms.

A case study on the stability of a big underground powerhouse cavern cut by an interlayer shear zone in the China Baihetan hydropower plant

AbstractThe big underground powerhouse cavern of the China Baihetan hydropower plant is 438 m long, 34 m wide, and 88.7 m high. It is cut by a weak interlayer shear zone and its high sidewall poses a huge stability problem. This paper reports our successful solution of this problem through numerical simulations and a replacement‐tunnel scheme in the detailed design stage and close site monitoring in the excavation stage. Particularly, in the detail design stage, mechanical parameters of the shear zone were carefully determined through laboratory experiments and site tests. Then, deformation of the surrounding rocks and the shear zone under high in situ stress conditions was predicted using 3 Dimensional Distinct Element Code (3DEC). Subsequently, a replacement‐tunnel scheme was proposed for the treatment on the shear zone to prevent severe unloading relaxation of surrounding rocks. In the construction period, excavation responses were closely monitored on deformations of surrounding rocks and the shear zone. The effect of local cracking in the replacement tunnels on sidewall stability was evaluated using the strength reduction method. These monitoring results were compared with the predicted numerical simulation in the detailed design stage. It is found that the shear zone greatly modified the deformation mode of the cavern surrounding rocks. Without any treatment, rock mass deformation on the downstream sidewall was larger than 125 mm and the shearing deformation of the shear zone was 60–70 mm. These preset replacement tunnels can reduce not only the unloading and relaxation of rock masses but also the maximum shearing deformation of the shear zone by 10–20 mm. The predictions by numerical simulation were in good agreement with the monitoring results. The proposed tunnel‐replacement scheme can not only restrain the shear zone deformation but also enhance the safety of surrounding rocks and concrete tunnels. This design procedure offers a good reference for interaction between a big underground cavern and a weak layer zone in the future.

Pressure monitoring and deformation analysis of a brine-filled salt cavern - A case study of Jianghan, China

Deformation stability plays a crucial role in the construction and operation of hydrocarbon storage salt caverns. Monitoring and analyzing wellhead pressure data under closed well conditions provide valuable insights into the evaluation of cavern stability. In this study, an analysis of five factors influencing internal brine pressure was conducted, which revealed that salt creep and brine thermal expansion contribute to an increase in internal brine pressure, while brine penetration, salt dissolution and cavern compressibility lead to a decrease. The atmospheric pressure, surface temperature, and earth tides are the other factors influencing wellhead pressure fluctuations. Through the analysis of wellhead pressure data from brine-filled salt cavern under closed conditions, the creep parameters of salt formations were calculated using the Norton-Hoff creep model. To further understand the behavior of brine-filled salt cavern, numerical simulations were conducted to analyze the brine temperature distribution and deformation stability. The volumetric average temperature of the brine exhibits a nonlinear increase over time, characterized by a rapid rise in the initial stage followed by a slower rate of change. Finally, the brine temperature stabilizes and converges toward the formation temperature. From deformation analysis, it was observed that the closed condition is more favorable for the maintenance of cavern stability in comparison to the open condition. The subsidence of the cavern roof is particularly sensitive and requires careful attention, especially for flat-roof salt caverns. Under closed conditions, it is crucial to control the brine pressure decrease rate and prevent excessively high brine pressure to ensure the mechanical and penetration stability of the brine-filled salt cavern. These findings provide valuable insights into the factors influencing internal brine pressure, the determination of creep parameters, the distribution of brine temperature, and the deformation stability of brine-filled salt caverns.

Investigation of the Tunnel Water Inflow Prediction Method Based on the MODFLOW-DRAIN Module

The precise prediction of water inflow in tunnels is a key aspect of civil engineering. It is directly related to the progress of construction, the stability of caverns and the safety of construction personnel. Currently widely used calculation methods include the empirical formula method, the analytical method and the numerical simulation method. In situations where the geological conditions are complex and the boundary conditions are irregular, numerical simulation methods have clear advantages. However, there are still discussions about the rationality, accuracy and applicability of the definition of tunnel boundary conditions in numerical simulation methods. Based on the Groundwater Model System (GMS10.7) software, we investigate the feasibility of using the DRAIN module to describe tunnel boundaries and calculate water inflow. By comparing it with traditional empirical formulas, the feasibility of the method is verified, and the error is about 8%. In addition, this method is applied to predict the water inflow in the construction of an underground pumping station under complex geological conditions and the prediction is compared with the measured water inflow, and the error is about 20%. This article explains the physical meaning and the value method of the relevant parameters of the DRAIN module, providing basic support for the prediction of water inflow in tunnels.

Fracture characteristics and failure analysis of underground caverns in Longyou Grottoes induced by humidity change

The stability of underground caverns in Longyou Grottoes is greatly influenced by humidity change since they are mainly composed of argillaceous siltstone, which is highly sensitive to moisture. In this work, a coupled humidity diffusion-fracture model based on the finite-discrete element method (FDEM) is employed to study the fracture characteristics and failure mechanism of underground caverns No.4 and 5 in Longyou Grottoes caused by humidity change. Firstly, the input parameters are calibrated with uniaxial compression tests on dry and saturated siltstone samples. Then, the crack patterns of a disk rock specimen under drying and wetting conditions are captured, which are compared with experimental results to validate the coupled model and calibrated parameters. Subsequently, numerical models of caverns No.4 and 5 in the Longyou Grottoes are established based on the actual geological conditions. The failure processes of the surrounding rock under moisture evaporation and absorption are simulated, and the stress evolution caused by humidity change is analyzed. Results indicate that the caverns primarily experience tensile failure under the evaporation condition with a distribution of tensile cracks perpendicular to the cavern boundary. However, shear failure occurs in the surrounding rock during moisture absorption, and the cracks gradually propagate into the deep rock with humidity diffusion. Additionally, several factors influencing cavern failure are discussed, and the number of cracks is quantified. Finally, it is found that the HDZ of the caverns is significantly reduced by applying a waterproof layer. The findings in this study have practical significance for predicting and improving the long-term stability of underground caverns in Longyou Grottoes under complex climatic conditions.

Stability of the horizontal salt cavern used for different energy storage under varying geological conditions

Horizontal salt cavern underground energy storage (UES) is a key focus for future energy storage facility development in China. The country is actively advancing the implementation of salt cavern UES in various sectors. This research aims to investigate the impact of different gas frequencies on the stability of horizontal salt cavern UES under varying geological conditions. A creep model was developed to analyze salt rock behavior during gas cycling operations. Fatigue tests were conducted on salt rock, varying the frequencies and amplitudes. The experimental findings determined the model parameters and validated its accuracy. Geological-mechanical models were created for different energy storage areas: Yunying, Jintan, and Chuzhou. The creep model was implemented to analyze the stability of salt cavern UES under three scenarios: compressed air energy storage (high frequency), natural gas storage (moderate frequency), and helium storage (low frequency). The construction of horizontal salt cavern UES was optimized considering economic and stability factors associated with the geological conditions. Results indicated that shallow salt mines are suitable for compressed air energy storage, middle-depth salt mines are better for natural gas storage, and deep salt mines are appropriate for helium storage. This research provides valuable insights into optimizing the construction of horizontal salt cavern UES in China, considering different frequencies and geological conditions.

Numerical Simulation Study on the Mechanics and Pore Characteristics of Tectonically Deformed Coal under Multi-Level and Multi-Cycle Loading and Unloading Conditions

Horizontal well cavern completion and stress release is considered a potential technique for efficient development of coalbed methane in tectonically deformed coal (TDC). Pulsating loading and unloading is a key technique for the controlled expansion of caverns and broader stress release within the reservoir. However, current understanding of the mechanical characteristics and pore network structure evolution of TDC under cyclic loading and unloading conditions is still limited. This paper employs numerical simulation methods to study the mechanical behavior and damage characteristics of TDC under cyclic loading and unloading. After obtaining a set of micromechanical parameters reflecting the behavior of TDC samples under triaxial compression in high-stress states, the effects of different stress gradients and cyclic amplitudes on the stress–strain curve, porosity changes, and crack propagation in TDC samples were analyzed. The study results indicate that under various cyclic loading and unloading conditions, the mechanical response characteristics of TDC samples are broadly similar, primarily divided into compression, slow expansion, and accelerated expansion phases. Under low unloading level conditions, the volume expansion of TDC samples is minimal. Also, at the same unloading level, the strain increment decreases with an increasing number of cycles. Correspondingly, under these conditions, the porosity and microcrack expansion in TDC are less than in high-stress gradient scenarios. Under the same unloading level but different amplitudes, the volume expansion rate at 50% unloading amplitude is higher than at 1 MPa unloading amplitude for TDC, with an increased number of crack expansions. Therefore, under cyclic loading conditions, the sensitivity of crack propagation within TDC samples to amplitude is greater than that to unloading level. Under actual pulsating excitation conditions, a low-amplitude, low-stress gradient pulsation method should be used to maintain the stability of horizontal well caverns, and gradually increase the cyclic amplitude to achieve the efficient extraction of coalbed methane in TDC reservoirs. The findings of this study can serve as an important reference for optimizing process parameters in cyclic pulsating stress release engineering for TDC.

Optimization of construction parameters for salt cavern underground energy storage using horizontal multi-stage leaching method under different geological conditions

In the thinly bedded rock salt in China, the horizontal multi-stage leaching method can solve the problem of constructing an enormous salt cavern in a limited-thickness salt layer. This study aims to investigate how the final morphology of the salt cavern is affected by critical parameters of the horizontal multi-stage leaching method in different geological conditions, such as the injection rate of freshwater, tubing withdrawal distance, and the presence of an air cushion. These parameters will directly affect the shape of the salt cavern and thus influence the stability of the underground energy storage during operation. We developed an experimental device to simulate salt cavern construction using the horizontal multi-stage leaching method. By varying the parameters in the physical simulation experiments, horizontal salt cavern models with different morphologies and sizes were obtained. The experiments showed that the injection rate of freshwater significantly affects the salt cavern height, length, and maximum width. The tubing withdrawal distance mainly affects the shape of the roof of the salt cavern. Introducing an air cushion to protect the salt cavern roof during the construction of a salt cavern results in a “flat-top” morphology. According to these salt cavern models, three-dimensional geomechanical models for different geological conditions are established by using reverse engineering technology. After analyzing the stability of the salt cavern using these models and combining them with economic factors to recommend values for the critical parameters of horizontal multi-stage leaching method in different geological conditions are provided. This study provides reference value and a theoretical basis for the construction of horizontal salt cavern storage in China.