<p>Operational asymmetry is a characteristic feature of underground gas storage systems—particularly in salt caverns—where injection phases are typically longer and steadier than withdrawal phases, which are shorter and more abrupt in response to fluctuating energy demand. Although large-scale hydrogen storage in salt caverns remains under development, similar operational patterns are expected once such facilities reach commercial maturity. Understanding how this temporal imbalance influences the mechanical behavior of the cavern is therefore important for assessing long-term performance and supporting the development of robust operational strategies for future hydrogen storage systems. This study introduces a dimensionless parameter, <InlineEquation ID="IEq1"> <EquationSource Format="TEX">$\chi$</EquationSource> <EquationSource Format="MATHML"><math> <mi>χ</mi> </math></EquationSource> </InlineEquation>, to quantify the temporal asymmetry of cyclic operation. <InlineEquation ID="IEq2"> <EquationSource Format="TEX">$\chi$</EquationSource> <EquationSource Format="MATHML"><math> <mi>χ</mi> </math></EquationSource> </InlineEquation> accounts not only for the residence times at high and low pressure but also for the different durations required to reach and leave these states. It is defined as the ratio between the effective durations of the high-pressure and low-pressure stages, including both the holding and ramp phases. The numerator represents the total effective time under high-pressure conditions, which includes half of the filling ramp (<InlineEquation ID="IEq3"> <EquationSource Format="TEX">$0.5\tau _{\textrm{fill}}$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>0.5</mn> <msub> <mi>τ</mi> <mtext>fill</mtext> </msub> </mrow> </math></EquationSource> </InlineEquation>), the full high-pressure holding phase (<InlineEquation ID="IEq4"> <EquationSource Format="TEX">$\tau _{\textrm{WGP}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>τ</mi> <mtext>WGP</mtext> </msub> </math></EquationSource> </InlineEquation>), and half of the withdrawal ramp (<InlineEquation ID="IEq5"> <EquationSource Format="TEX">$0.5\tau _{\textrm{withd}}$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>0.5</mn> <msub> <mi>τ</mi> <mtext>withd</mtext> </msub> </mrow> </math></EquationSource> </InlineEquation>). The denominator represents the complementary effective time under low-pressure conditions, consisting of half of the withdrawal ramp (<InlineEquation ID="IEq6"> <EquationSource Format="TEX">$0.5\tau _{\textrm{withd}}$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>0.5</mn> <msub> <mi>τ</mi> <mtext>withd</mtext> </msub> </mrow> </math></EquationSource> </InlineEquation>), the full low-pressure holding phase (<InlineEquation ID="IEq7"> <EquationSource Format="TEX">$\tau _{\textrm{CGP}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>τ</mi> <mtext>CGP</mtext> </msub> </math></EquationSource> </InlineEquation>), and half of the filling ramp (<InlineEquation ID="IEq8"> <EquationSource Format="TEX">$0.5\tau _{\textrm{fill}}$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>0.5</mn> <msub> <mi>τ</mi> <mtext>fill</mtext> </msub> </mrow> </math></EquationSource> </InlineEquation>). This formulation expresses the relative exposure of the rock mass to stabilizing (high-stress) versus creep-promoting (low-stress) conditions. A viscoplastic numerical model is developed to reproduce the time-dependent deformation of salt under cyclic pressure loading. Two measurable indicators are analyzed: the vertical displacement at the cavern roof and the volumetric shrinkage associated with viscoplastic mechanical closure. Results show that temporal asymmetry exerts a dominant control on the long-term mechanical response. When <InlineEquation ID="IEq9"> <EquationSource Format="TEX">$\chi >1$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>χ</mi> <mo>></mo> <mn>1</mn> </mrow> </math></EquationSource> </InlineEquation>—indicating longer effective residence under high pressure—the cavern evolves toward a steady-state equilibrium; when <InlineEquation ID="IEq10"> <EquationSource Format="TEX">$\chi <1$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>χ</mi> <mo><</mo> <mn>1</mn> </mrow> </math></EquationSource> </InlineEquation>, viscoplastic deformation accelerates, promoting mechanical cavern closure. Comparison with operational data from three underground natural gas storage facilities in the United States is consistent with the <InlineEquation ID="IEq11"> <EquationSource Format="TEX">$\chi$</EquationSource> <EquationSource Format="MATHML"><math> <mi>χ</mi> </math></EquationSource> </InlineEquation> ranges associated with the stability trends identified numerically. The proposed parameter <InlineEquation ID="IEq12"> <EquationSource Format="TEX">$\chi$</EquationSource> <EquationSource Format="MATHML"><math> <mi>χ</mi> </math></EquationSource> </InlineEquation> thus provides a physically grounded and practical first-order indicator for interpreting long-term mechanical behavior in future hydrogen storage caverns, within the scope of a viscoplastic mechanical framework that does not explicitly account for damage, fracture, or permeability evolution.</p>

Viscoplastic analysis of salt caverns under time-asymmetric pressure cycling: definition and validation of a dimensionless operational rule with real case data

期刊论文
出版模式：

开放获取
发表日期：

2026年6月19日

查看全文

F. Figueiras,
M. A. Soage-Quintans,
J. París,
I. Colominas,
L. Cueto-Felgueroso

摘要

Operational asymmetry is a characteristic feature of underground gas storage systems—particularly in salt caverns—where injection phases are typically longer and steadier than withdrawal phases, which are shorter and more abrupt in response to fluctuating energy demand. Although large-scale hydrogen storage in salt caverns remains under development, similar operational patterns are expected once such facilities reach commercial maturity. Understanding how this temporal imbalance influences the mechanical behavior of the cavern is therefore important for assessing long-term performance and supporting the development of robust operational strategies for future hydrogen storage systems. This study introduces a dimensionless parameter, $\chi$ $χ$ , to quantify the temporal asymmetry of cyclic operation. $\chi$ $χ$ accounts not only for the residence times at high and low pressure but also for the different durations required to reach and leave these states. It is defined as the ratio between the effective durations of the high-pressure and low-pressure stages, including both the holding and ramp phases. The numerator represents the total effective time under high-pressure conditions, which includes half of the filling ramp ( $0.5\tau _{\textrm{fill}}$ $0.5 τ_{fill}$ ), the full high-pressure holding phase ( $\tau _{\textrm{WGP}}$ $τ_{WGP}$ ), and half of the withdrawal ramp ( $0.5\tau _{\textrm{withd}}$ $0.5 τ_{withd}$ ). The denominator represents the complementary effective time under low-pressure conditions, consisting of half of the withdrawal ramp ( $0.5\tau _{\textrm{withd}}$ $0.5 τ_{withd}$ ), the full low-pressure holding phase ( $\tau _{\textrm{CGP}}$ $τ_{CGP}$ ), and half of the filling ramp ( $0.5\tau _{\textrm{fill}}$ $0.5 τ_{fill}$ ). This formulation expresses the relative exposure of the rock mass to stabilizing (high-stress) versus creep-promoting (low-stress) conditions. A viscoplastic numerical model is developed to reproduce the time-dependent deformation of salt under cyclic pressure loading. Two measurable indicators are analyzed: the vertical displacement at the cavern roof and the volumetric shrinkage associated with viscoplastic mechanical closure. Results show that temporal asymmetry exerts a dominant control on the long-term mechanical response. When $\chi >1$ $χ > 1$ —indicating longer effective residence under high pressure—the cavern evolves toward a steady-state equilibrium; when $\chi <1$ $χ < 1$ , viscoplastic deformation accelerates, promoting mechanical cavern closure. Comparison with operational data from three underground natural gas storage facilities in the United States is consistent with the $\chi$ $χ$ ranges associated with the stability trends identified numerically. The proposed parameter $\chi$ $χ$ thus provides a physically grounded and practical first-order indicator for interpreting long-term mechanical behavior in future hydrogen storage caverns, within the scope of a viscoplastic mechanical framework that does not explicitly account for damage, fracture, or permeability evolution.