Abstract
A mathematical model is established to predict the salt cavern development during leaching in high-insoluble salt formations. The salt-brine mass transfer rate is introduced, and the effects of the insoluble sediments on the development of the cavern are included. Considering the salt mass conservation in the cavern, the couple equations of the cavern shape, brine concentration and brine velocity are derived. According to the falling and accumulating rules of the insoluble particles, the governing equations of the insoluble sediments are deduced. A computer program using VC++ language is developed to obtain the numerical solution of these equations. To verify the proposed model, the leaching processes of two salt caverns of Jintan underground gas storage are simulated by the program, using the actual geological and technological parameters. The same simulation is performed by the current mainstream leaching software in China. The simulation results of the two programs are compared with the available field data. It shows that the proposed software is more accurate on the shape prediction of the cavern bottom and roof, which demonstrates the reliability and applicability of the model.
Similar content being viewed by others
Introduction
The usage of natural gas resources strongly fluctuates according to the season, while the gas well production is stable1,2. To balance the mismatch in gas supply and demand, the underground gas storage (UGS) has been used effectively as peak shaving means for nearly a century3,4,5. As one of the UGSs, salt cavern gas storage is constructed by cavern leaching in the sedimentary formation of salt rock, which formed from the drying up of enclosed salt lakes and seas6. Using one or more drilled wells, fresh water is injected into the salt formation to dissolve the salt deposits, brines are pumped out, and a cavern is developed gradually7. A salt cavern can be leached to hundreds of thousands of cubic meters, and it is an excellent choice for gas or oil storage, due to the good creep property, self-healing characteristics and extremely low permeability rate of salt rock8,9,10,11,12. There are 74 salt cavern UGSs in the world and the working gas volume is about 1.62 × 1010 m3, according to the research of Ozarslan13.
During the construction of salt cavern UGSs, the shapes of the caverns are the most important concern for the safety and stability of the cavern, which are always the primary considered issues14,15. According to the literatures16,17,18, an ellipsoidal shape is the most optimal structure for cavern stability. To ensure the salt cavern grows into the ideal shape, researchers have conducted many studies on the design and control of cavern shape during leaching. Duire and Jessen19 studied the solution mechanism of salt rock by laboratory experiments. They proposed an experience equation of the solution rate of salt rock. Nolen et al.20 designed a new digital program to simulate the cavern leaching process, which has been used successfully for years in Germany. The program took time, produced brine concentration and increased cavern volume as calculating parameters. Meanwhile, a certain amount of insolubles is considered. Reda and Russo21 conducted three salt cavity leaching experiments in the laboratory. They kept the cavity pressurized to actual cavern pressure and recorded the brine salinity and salt-wall recession. The experimental results are in good agreement with the numerical predictions. O’Hern et al.22 provided an experimental data record of measuring the leaching rates of salt walls. They found the dividing line between upper and lower regions of brines, roughly above and below the fresh water injection point. Li et al.23 set up a model to simulate the repair process of irregularly shaped salt cavern by re-leaching under gas. They conducted indoor experiments and field tests, the simulation results coincide well with the experimental results and field data, while the simulation of insolubles accumulation was not mentioned.
However, most of these researches only considered the case of cavern form control in pure or low- insoluble salt rock. For their target salt caverns are mostly constructed in salt dome or thick salt formations, where the content of water-insoluble substances is very low and the development of the cavern is less affected. However, in many places, the salt formations contain significant amounts of water-insoluble substances, such as mudstone, sandstone, gypsum and limestone24,25. Especially in China, where most salt formations have interbedded water-insoluble interlayers, the salt beds are thin and the content of insoluble substances are higher than 20%26. Take an example, cavern JT86 (Fig. 1) is one of the completed salt caverns in Jintan UGS of Jiangsu province, China. There are many insoluble interlayers in the salt formations, and the average insolubles content of the salt formations around the cavern is about 21.6%. During leaching, the insoluble substances detached from the salt wall, and accumulated on the cavern bottom as loose insoluble particles with larger volume. The volume content of the insoluble sediments in this cavern reached about 35%, and the cavern bottom rose about 56 meters. Actually, the bottom of the cavern boundary is actually the surface of the insoluble sediments. For most Chinese salt formations, the content of insoluble substances is too high to ignore, the accumulation of the insolubles has greatly limit the lower cavern growth and decreased cavern capacity27.
Little research has been performed on the construction of salt cavern in these high-insoluble salt formations. Li et al.27 proposed a prediction model of the insolubles accumulation shape during leaching in salt formations with high content of insolubles. However, the model is two-dimensional and the insolubles accumulation is not coupled with the development of cavern shape. Zheng et al.28 discussed the water-soluble mechanism of thick interlayers in salt cavern. They proposed that the collapse of interlayers would connect the salt deposits above and below, and consequently enlarge cavern volume. A prediction model of the collapse of the interlayers was presented, whereas the solution mechanism of salt rock was not mentioned. The most commonly used cavern leaching software in China, “WinUbro”, has no completed insolubles calculation module29. All of the simulated caverns have a flat bottom, which does not match the reality that the cavern bottom shapes are variable30. The deviation of the cavern bottom results in the inconformity of total cavern shape, volume and development. Thus the salt cavern has to be detected by sonar and corrective measures need to be taken after each leaching stage, which brings great schedule and financial pressure.
For the further development of the construction of salt cavern gas storage in China, this paper proposes a prediction model of the single-well cavern leaching process in high-content insolubles salt formations. The salt-brine mass transfer rate is introduced to be related to the brine concentration, salt-brine contact angle and temperature. Based on the law of conservation of mass and reasonable assumptions, the calculation equations of cavern radius, brine concentration, flow rate, cavern volume and insolubles height are deduced. According to the particle packing theory, the re-distribution equation of the insoluble sediments is derived. Using the finite difference method, a VC++ computer program is developed and the numerical solution of the model can be obtained. Two salt caverns in Jintan UGS are selected and the actual leaching processes are simulated by the developed program and by the current software WinUbro. The comparison of the simulation results and the field data demonstrates the validity and accuracy of the proposed model.
Methods
In the cavern leaching process, two concentric tubings are placed into the leaching cavern through a casing pipe (Fig. 2). The casing pipe is the flow passage of blanket material, usually diesel oil, which is used to protect the upper salt from being dissolved. The inner and outer tubings are used for water injection and brine discharge. The depth of the inner and outer tubings, the injection flow rate, and the oil-brine interface depth are adjusted in each leaching stage to obtain an ideal cavern shape. The fresh water can be injected through inner or outer leaching tubing, when the leaching process can be distinguished into “direct leaching mode” and “diverse leaching mode”. In this section, the equations of the developing cavern boundary of the two modes are derived respectively, and the control equations of the insolubles sediments surface are deduced.
Mathematical model of cavern boundary development
Direct leaching mode
Assume that the inner tubing orifice is placed at depth H in, outer tubing is at H out, and oil-brine interface is at H oil (Fig. 2). Fresh water is injected through inner tubing in the direct leaching mode, the flow rate is Q. In most cases, the brine concentration of the circulation zone between H in and H oil can be assumed as uniform during the direct leaching period22,23. For the water injection tube is at the bottom of the circulation zone, the injected fresh water is far lighter than the brines and it will rise. While the high saturated brines generated by the dissolution of cavern top and cavern wall are heavy and they will sink. The up flow of fresh water and the down flow of brines will be fully blended, thus the concentration of the brines in the circulation zone can be assumed as C in this mode. Since there are no fresh water flowed into the dead zone under H in, the brines in the dead zone will keep saturated.
Taking the whole circulation zone as a control volume, the mass of salt is conserved,
where, ρ is the desity of salt rock, V is the volume of the circulation zone, u is the content of insolubles in the salt formation, V d is the volume of the dead zone under H in, C s is the concentration of the saturated brines in the dead zone, Q is the flow of injected water, C is the concentration of the brines in the circulation zone, t is time.
The first term of Eq. (1) is the mass of the dissolved salt, the second term is the up flow of the dead zone due to the accumulation of insoluble sediments, and the third term is the mass of salt contained in the brine discharge through outer tubing.
\(\frac{\partial {V}_{d}}{\partial t}\) is the volume change of the dead zone, causing by the accumulation of the insolubles contained in the salt formations. It can be calculated by,
where, f is the expansion coefficient of insoluble substances in brine.
The volume of the circulation zone V can be calculated by,
where, H oil is the depth of oil-gas interface, H in is the depth of inner leaching tubing, R is the cavern radius, h is height, θ is circumferential angle.
The cavern radius R is related to the dissolution rate of cavern wall,
where, ω is the dissolution rate of the salt cavern wall, α is the angle between the normal direction of cavern wall and vertical upward direction, or salt-brine contact angle for short.
ω is actually the mass transfer rate between salt and brine, which is mainly related to the brine concentration, salt-brine contact angle, and temperature of brines31. According to the research of Xiao et al.32, ω can be calculated by an empirical formula. In addition, if the salt wall is covered by insoluble sediments, the salt is not exposed to the brines, and then the dissolution will stop.
where, T is the temperature of the brine, e is the natural logarithm, ω 0 is the dissolution rate of pure salt in fresh water, which can be easily measured in laboratory.
α can be calculated by,
Diverse leaching mode
During the diverse leaching period, the fresh water is injected through the higher outer tubing, while the brine is discharged through the lower inner tubing (Fig. 3). The circulation zone is divided into two parts, upper circulation zone above H out and lower circulation zone under H out. The distribution of the brine concentration in the upper circulation zone is just like that in the circulation zone of the direct leaching mode. Taking the upper circulation zone as a control volume, and the mass conservation equation of this part can be written as,
where, V can be calculated by,
where, H out is the depth of outer leaching tubing orifice.
The first term of Eq. (7) is the mass of the dissolved salt in the upper circulation zone, and the second term is the mass of salt contained in the brines discharged into the lower circulation zone. The term on the right side of Eq. (7) is the salt mass change of the control volume. Compared with Eq. (1), the up flow of the saturated brines in the dead zone due to the accumulation of insoluble sediments is discharged from the lower inner tubing, thus the related term does not appear in Eq. (7).
Due to the injection of fresh water, the volume of the brines in the upper circulation zone expands and the extra brines flow down to the lower circulation zone. The brines in the lower zone are pushed downward, dissolving salts along the way, and finally discharged through the inner tubing. The lower circulation zone can be divided into numerous horizontal micro sections, the heights of which are dh. Due to the horizontal stratification of brines, it can be assumed that the brines concentrations in each micro section are uniform23. The mass conservation equation of the micro sections can be expressed as,
where, h is depth as well as the code name of the micro section whose depth is h, C h,t is the brine concentration of section h, μ h,t is the vertical flow rate of section h, V h,t is the volume of section h.
V h,t can be written as,
Since the brine cannot be compressed, the brine flow μ h,t in Eq. (9) can be calculated by the volume conservation equation
Mathematical model of insoluble sediments surface
During cavern leaching, with salt dissolved into brines, the insoluble substances contained in the salt lose support. They detach from the cavern wall and accumulate on the cavern bottom. The salt cavern development is closely related to the insoluble sediments in the cavern. For the accumulation of the insolubles will cover the lower cavern wall and protect them from being dissolved. Furthermore, the bottom of a complete cavern is actually the surface of the insoluble sediments. In this section, the accumulation of the insoluble substances will be discussed.
Falling from the cavern wall, an insoluble particle will finally settle within a distance to the projection of its starting point on the cavern bottom. According to the research of Li et al.27, this distance is basically proportional to the height of the falling path, as shown in Eq. (12).
where, R fall is the maximum distance from the position of a fallen insoluble particle to the projection of its starting point on the cavern bottom, H fall is the height of the falling path of the insoluble particle, D is the scattering coefficient between R fall and H fall, the value of D can be obtained by laboratory experiments.
In other words, the insoluble substances, fallen from point (R, Θ, H), will accumulate within a radius around the projection of (R, Θ, H) on the cavern bottom (insoluble surface), as shown in Fig. 4(a). We can assume a set P r,θ,z to represent a number of cavern wall points, the insolubles from where may fall onto point (r, θ, z) on the insoluble sediments surface. Every point in the set will satisfy the following equation,
where, P r,θ,z is the set of the cavern wall points, the insolubles from where may fall onto (r, θ, z).
Assuming an even distribution of the insolubles in the settle circle, the height change of insoluble sediments surface can be calculated by,
where, U R,Θ,H is the insoluble content of the micro-element on point (R, Θ, H), \(\tfrac{\partial {V}_{R,\Theta ,H}}{\partial t}\) is the volume variation of the micro-element on point (R, Θ, H), ∈ is the expansion coefficient of the insoluble substances in brine, which has considered the effects of broken expansion and compaction.
According to the particle packing theory, the tangential angle of the insolubles surface should be always less than the repose angle of insoluble particles. To simplify the calculation, we only consider the tangential angle limit on the radial direction and circumferential direction, which can be expressed as,
where, α ins is the angle of repose of the insoluble particles in brines.
For all of the points on the insolubles surface which do not meet Eq. (15), the height of the insoluble sediments should be adjusted. For example, in the radial direction, the height of the insoluble sediments on point (γ, θ) is z γ,θ , it is too high and the unstable part fall down, as shown in Fig. 4(b). The height of the insolubles on this point decreases, while the insoluble heights of neighboring lower points increase. The volume conservation equation of the insoluble sediments can be expressed as,
where, z γ,θ is the insolubles height on point (γ, θ), z r,θ is the insolubles height on point (r, θ), r 0 and r 1 respectively are the radii of the leftmost and rightmost points affected by the fallen insolubles from (γ, θ).
The tangential angle of the re-distributed insolubles will be equal to α ins, thus
On point (r 0, θ), the height of insolubles does not change after redistribution,
The redistribution of the insoluble sediments is similar in the circumferential direction, and we can easily obtain the adjustment height in the circumferential directions. Thus, considering the redistribution, Eq. (14) should be re-written as,
where, dz r,θ and \({\rm{d}}{z}_{r,\theta }^{\prime} \) respectively are the height changes of the insolubles in the radial and circumferential directions at point (r, θ).
Software implementation for numerical solution
Eqs (1–19) provides the theoretical formulae for the two modes of the cavern leaching process. Given the technological parameters and the initial values of the control variables, the developing cavern boundary, brine concentration and insoluble sediments surface can be calculated. Based on the equations, a computer program “Single-well Salt Cavern Leaching Simulation V1.0”, “SSCLS” for short, is developed to solve the numerical solution of the model. Using finite difference method, SSCLS is written in the VC++ language, with independent operation interface as shown in Fig. 5. There are five function modules in SSCLS: formation information and sonar data import, technical parameters input, flow/concentration field calculation, cavern boundary calculation, and insoluble sediments surface calculation.
The cavern is divided into numerous micro-sectors, with different serial number, cylindrical coordinate, volume, brine flow velocity, brine concentration, salt-brine contact angle, and insolubles content. All of these sector information are stored and calculated in array in every time step, until the end of the leaching process. In the cavern system, a grid of the insoluble sediments surface is set up. The height of the insolubles on each grid is calculated and amended after the change of cavern shape. The calculation flow chart of the whole simulation process is shown in Fig. 6.
Data availability
The data generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Results
Since the whole construction period is quite long (3–7 years), we cannot verify the proposed model by simulation-construction-comparison. Second best, we can simulate the leaching process of a complete salt cavern, using the practical geological and technological parameters, and then compare the simulation results with the actual cavern shape to verify the model. Two complete salt caverns of Jintan UGS in Jiangsu province, a regular shaped cavern JT52 and an irregular shaped cavern JT103, are selected for simulation. The same simulation will be performed by WinUbro, the most commonly used cavern leaching software in China33. The simulation results are shown in this section.
Simulation results of the regular salt cavern JT52
Cavern JT52 is one of the earliest constructed salt caverns of Jintan UGS in Jiangsu, China. The buried depth range of the salt formations is between −1140 m and −990 m, the average salt purity is about 78.8%. There are about 7 insoluble interlayers in the formation, the major components of the interlayers are mudstone and gypsum. Cavern JT52 was built between 2003 and 2010, the volume of the complete cavern is about 186,675 m3 and the depth range is between −1089.5 m and −1015.4 m. The whole cavern construction process has 7 direct leaching stages with different technological parameters, including inner/outer tubing depth, oil-brine interface depth, stage time and water injection rate, as shown in Table 1.
The practical leaching parameters and geological data are imported into SSCLS and WinUbro, where the seven leaching stages are simulated. The comparison of the simulation results and the actual cavern shape is shown in Fig. 7(a). Since cavern JT52 is roughly axisymmetric, only the vertical section in the north-south direction is compared. Since there is no sonar data of the first leaching stage, only cavern shape comparison of stages 2–7 are shown.
Simulation results of the irregular salt cavern JT103
Cavern JT103 is a lately completed salt cavern in Jintan UGS, started in 2012 and completed in 2016. The buried depth range of the salt formations is from −1050 m to −880 m, with 12 insoluble mudstone interlayers. The average salt purity is about 84.5%. The volume of the complete cavern is about 181275 m3 and the depth range is between −1013.4 m and −923.4 m. The whole cavern construction process has 1 direct leaching stages and 4 diverse leaching stages, with different technological parameters (Table 2).
The leaching processes are simulated by SSCLS and WinUbro at the same time, using the actual geological and technological parameters. The simulated cavern shapes of the two programs and the actual cavern shapes (average radii) are drawn in Fig. 7(b). Since there is no sonar data of the 1st and 4th leaching stage, only cavern shape comparison of stages 2, 3, 5 are shown.
Discussion
Comparison of the simulation results and the field data of JT52
As shown in Fig. 7(a), the lateral dimensions of the simulated caverns by using the two simulators (SSCLS and WinUbro) are both fairly accurate. The maximum radius errors of the two simulations are less than 2 meters (Fig. 8(a)), accounting for about 5% of the maximum radius, which can meet the requirements of actual engineering. However, the simulations are quite different on the cavern top and bottom. The cavern bottom shapes of the simulation of SSCLS are closer to the actual shapes, the errors of the lowest depth are about 1–3 meters (Fig. 8(a)), which is a pretty high accuracy. Since WinUbro has no function module about the insolubles accumulation, the cavern bottoms are handled as flat, which have a large deviation from the actual shapes. There are 8–11 meters errors between the lowest depth of WinUbro’s simulations and the actual caverns. As for the top part of the cavern, the SSCLS simulated shapes basically coincide with the actual shapes, while the WinUbro’s simulated cavern has an inflection point after stage 4 (Fig. 7(a)). The reason is that, the control points of the cavern boundary of SSCLS move in the normal direction, which is the same as that of the actual cavern wall. The cavern control points of WinUbro moves in the radial direction, resulting in the rough of the cavern top boundary. Due to the perfectly match on the cavern top and bottom, the volume error of the simulation of SSCLS (about 2400 m3) is much less than that of WinUbro (about 23,000 m3), as shown in Fig. 8(b).
After stage 4, a protuberance appeared in the middle of the cavern bottom of the SSCLS simulated cavern (Fig. 7(a)), for the insolubles in the top center of the cavern fell and accumulate on the bottom center. However, the actual cavern bottom is pretty smooth in the middle part. The reason is that, fresh water was injected through the lower inner tubing near the insoluble surface, the insoluble sediments were re-distributed due to the injection flow. For example, the inner tubing was placed at −1099.1 m at the 4th stage, when the middle point of the insoluble surface was at −1098.8 m. The injection flow smoothen the cavern bottom, which is not considered in the proposed model and SSCLS at this stage.
In addition, the actual cavern shape has a “neck” at depth −1044 m in the 4th stage, which probably results from the effect of the insoluble interlayers. The SSCLS simulated cavern shape is pretty smooth at −1044 m, for the proposed model only considered the simulation of cavern shape change due to salt dissolution. The dissolution, soften and collapse of the insoluble interlayers is not included. The effect of these insoluble interlayers on the cavern development might be discussed in the next step of our work.
Comparison of the simulation results and the field data of JT103
As shown in Fig. 7(b), the SSCLS simulated cavern shapes basically coincide with the actual shapes of each stage. The radius errors are less than 2 meters (Fig. 9(a)), about 5% of the maximum cavern radius. The errors of the lowest cavern depth are less than 1.6 m in stage 2 and 5, about 1.8% of the cavern height. However, in stage 3, the lowest cavern depth error is about 9 meters (Fig. 9(a)). The simulated cavern bottom is lower and steeper than the actual cavern bottom. The inner tubing was placed at −1036.3 m at the 3rd stage, buried under the middle point of the insoluble surface, which rise from −1031.8 m to −1022.6 m. The brine discharge current around the inner tubing washed the insolubles from both sides to the center part. As mentioned in the comparison of cavern JT52, the impact of the flow field on the insoluble distribution was not considered in SSCLS, resulting in the difference between the simulation and actual cavern bottom in stage 3. The effect of the injection flow on the redistribution of insoluble sediments might be discussed in our next step work. However, the effect does not change the total volume of the cavern. The effective volume of the simulation results of SSCLS still has a high accuracy, the final volume error is less than 3000 m3 (Fig. 9(b)), about 1.6% of the total cavern volume.
The WinUbro simulated shapes have large deviation in the lower part of the cavern. In the second stage, the actual cavern is roughly obconical, the flat bottom processing mode of WinUbro results in large errors in maximum radius, volume and lowest point of the cavern. The error accumulates in the following stages. In the final stage, the radius error is 10 m, the volume error is about 19,000 m3, and the lowest depth error is 9.5 m (Fig. 9).
In general, the proposed model can accurately predict the cavern development and insolubles accumulation during leaching for gas storage within engineering precision limitation. Compared with WinUbro, the simulation results of SSCLS are more consistent with the actual shape on cavern bottom and top.
Conclusions
-
(I)
A mathematical model is proposed to predict the cavern development during cavern leaching for gas storage. The salt-brine mass transfer rate is introduced to be related to the brine concentration, salt-brine contact angle and temperature. The salt mass conservation is considered, and the couple equations of the cavern shape, brine concentration and flow rate of the direct and diverse leaching modes are deduced.
-
(II)
According to the falling and accumulating rules of the insoluble substances in brines, the calculation equations of the surface of the insoluble sediments are introduced for the first time. The redistribution of the insoluble sediments is considered according to the particle packing theory. The amendatory control equations of the insoluble sediments surface are derived.
-
(III)
A VC++ software SSCLS is developed to obtain the numerical solution of the proposed model. Using the actual geological data and technological parameters, the leaching processes of two salt caverns of Jintan UGS are simulated by SSCLS and by WinUbro, the most commonly used cavern leaching software in China. The simulation results of SSCLS coincide well with the actual cavern shape, and are more accurate on the shape prediction of the cavern bottom and top. It shows that the proposed model is reliable for the simulation of cavern leaching in China.
References
Kazeem, A. L., Mathilda, I. O., Stella, I. E., Saka, M. & Ayodeji, T. A. Underground storage as a solution for stranded associated gas in oil fields. J. Petrol. Sci. Eng. 150, 366–375 (2017).
Arfaee, M. I. R. & Sola, B. S. 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 (2014).
Danel, R., Otte, L., Vancura, V. & Řepka, M. Monitoring and balance of gas flow in underground gas storage. Pr. Earth Pl. Sci. 6, 485–491 (2013).
Shi, X. et al. Influences of filling abandoned salt caverns with alkali wastes on surface subsidence. Environ. Earth Sci. 73, 6939–6950 (2015).
Lawal, K. A., Ovuru, M. I., Eyitayo, S. I., Matemilola, S. & Adeniyi, A. T. Underground storage as a solution for stranded associated gas in oil fields. J. Petrol. Sci. Eng. 150, 366–375 (2016).
Li, Y. P., Liu, W., Yang, C. H. & Daemen, J. J. K. Experimental investigation of mechanical behavior of bedded rock salt containing inclined interlayer. Int. J. Rock Mech. Min. 69, 39–49 (2014).
Liu, W. et al. Investigation on the permeability characteristics of bedded salt rocks and the tightness of natural gas caverns in such formations. J. Nat. Gas sci. Eng. 35, 468–82 (2016).
Yuan, G. J. et al. Review of underground gas storage in the bedded salt deposit in china. (Society of Petroleum Engineers, 2006).
Ma, H. L. et al. Stability evaluation of the underground gas storage in rock salts based on new partitions of the surrounding rock. Environ. Earth Sci. 73, 6911–6925 (2015).
Chen, J., Ren, S., Yang, C. H., Jiang, D. Y. & Li, L. Self-healing characteristics of damaged rock salt under different healing conditions. Materials 6, 3438–3450 (2013).
Liu, W. et al. Permeability characteristics of mudstone cap rock and interlayers in bedded salt formations and tightness assessment for underground gas storage caverns. Eng. Geol. 193(2), 212–223 (2015).
Guo, Y. T., Yang, C. H. & Mao, H. J. Mechanical properties of Jintan mine rock salt under complex stress paths. Int. J.Rock Mech. Min. 56(12), 54–61 (2012).
Ozarslan, A. Large-scale hydrogen energy storage in salt caverns. Int. J. Hydrogen. Energ. 37(19), 14265–14277 (2012).
Wang, T. T. et al. Dynamic response of underground gas storage salt cavern under seismic loads. Tunn. Undergr. Sp. Tech. 43, 241–252 (2014).
Wang, T. T. et al. A new shape design method of salt cavern used as underground gas storage. Appl. Energ. 104, 50–61 (2013).
Charnavel, Y., Leca, D. & Poulain, F. Advanced geometrical modelling of salt dissolution during cavern leaching-Illustration with a case study. (Solution Mining Research Institute Meeting, 1999).
Cristescu, N. D. & Paraschiv, I. The optimal shape of rectangular-like caverns. Int. J.Rock Mech. Min. Geomech. Abstr. 32(4), 285–300 (1995).
Staudtmeister, K. & Rokahr, R. B. Rock mechanical design of storage caverns for natural gas in rock salt mass. Int. J.Rock Mech. Min. 34(3-4), 300 (1997).
Durie, R. W. & Jessen, F. W. Mechanism of the dissolution of salt in the formation of underground salt cavities. SPE J. 4(2), 183–190 (1964).
Nolen, J. S., von Hantlemann, O., Meister, S., Kleinitz, W. & Heiblinger, J. Numerical Simulation of the Solution Mining Process. (SPE European Spring Meeting, 1974).
Reda, D. C. & Russo, A. J. Experimental Studies of salt-cavity leaching by freshwater injection. SPE Prod. Eng. 1(1), 82–86 (1986).
Webb, S. W., O’Hern, T. J. & Hartenberger, J. D. SPR salt wall leaching experiments in lab-scale vessel- data report. (Office of Scientific & Technical Information Technical Reports, 2010).
Li, J. L. et al. Repair of irregularly shaped salt cavern gas storage by re-leaching under gas blanket. J. Nat. Gas Sci. Eng. 45, 848–859 (2017).
Wang, T. T. et al. Failure analysis of overhanging blocks in the walls of a gas storage salt cavern: a case study. Rock Mech. Rock Eng. 50(1), 1–13 (2016).
Xu, F. et al. Effect of confining pressure on the mechanical properties of thermally treated sandstone. Curr. Sci. India 112(6), 1101–1106 (2017).
Wang, T. T. et al. Failure analysis of thick interlayer from leaching of bedded salt caverns. Int. J. Rock Mech. Min. 73, 175–183 (2015).
Li, J. L. et al. A prediction model of the accumulation shape of insoluble sediments during the leaching of salt cavern for gas storage. J. Nat. Gas Sci. Eng. 33, 792–802 (2016).
Zheng, Y. et al. Solution mining technology of enlarging space for thick-sandwich salt cavern storage. Petrol. Explor. Dev. 44(1), 139–145 (2017).
Wang, W. Q. et al. Effect of cavity displacement of salt cavern gas storage on brine displacing concentration and cavity form (in Chinese). Oil Gas Stor. Transport. 34(2), 175–179 (2015).
Charnavel, Y. & Lubin, N. Insoluble deposit in salt cavern - test case. (Proc of Solution Mining Research Institute Fall Meeting, 2002).
Ma, H. L., Chen, F., Yang, C. H., Shi, X. L. & Zhang, C. Experimental research on dissolution rate of salt rock in depth (in Chinese). Min. Res. Dev. 30(5), 9–13 (2010).
Xiao, C. F., Yang, Y. K., Wu, G., Li, C. P. & Qiu, X. D. The dissolution characteristics of rock salt and mass transmission process (in Chinese). J. Chongqing. Univ. 16(2), 51–57 (1993).
Kunstman, A. & Urbanczyk, K. Application of winubro software to modelling of carven development in trona deposit. (Solution Mining Research Institute Technical Conference, 2009).
Acknowledgements
The authors wish to acknowledge the financial supports of National Natural Science Foundation of China (Grant No. 51774266; 41502296; 51404241), National Natural Science Foundation of China Innovative Research Team (Grant No. 51621006), Youth Innovation Promotion Association CAS (Grant No. 2016296), and Natural Science Foundation for Innovation Group of Hubei Province, China (Grant No. 2016CFA014). The authors are sincerely grateful to Jaak J Daemen, Mackay School of Earth Sciences and Engineering, University of Nevada, for his thoughtful proofreading of this paper.
Author information
Authors and Affiliations
Contributions
J.L. wrote the main manuscript text, and is responsible for the mathematical model and the software implementation; X.S. and T.W. have made contributions to the conception and design of the work. C.Y., Y.L., and H.M. have made contributions to the simulation analysis and comparison. All authors discussed the results and critically reviewed the manuscript.
Corresponding authors
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Li, J., Shi, X., Yang, C. et al. Mathematical model of salt cavern leaching for gas storage in high-insoluble salt formations. Sci Rep 8, 372 (2018). https://doi.org/10.1038/s41598-017-18546-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-18546-w
This article is cited by
-
Geomechanical simulation of energy storage in salt formations
Scientific Reports (2021)
-
Stability evaluation of underground gas storage salt caverns with micro-leakage interlayer in bedded rock salt of Jintan, China
Acta Geotechnica (2020)
-
Investigation on the Deformation and Strength Characteristics of Rock Salt Under Different Confining Pressures
Geotechnical and Geological Engineering (2020)
-
Prediction of pressure and temperature changes in the salt caverns of Tuz Golu underground natural gas storage site while withdrawing or injecting natural gas by numerical simulations
Arabian Journal of Geosciences (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.