Influence of the Variability of Compressed Air Temperature on Selected Parameters of the Deformation-Stress State of the Rock Mass Around a CAES Salt Cavern

The article presents the results of a numerical simulation of the deformation-stress state in the rock mass around a salt cavern which is a part of a CAES installation (Compressed Air Energy Storage). The model is based on the parameters of the Huntorf power plant installation. The influence of temperature and salt-creep speed on the stability of the storage cavern was determined on the basis of the three different stress criteria and the effort of the rock mass in three points of the cavern at different time intervals. The analysis includes two creep speeds, which represent two different types of salt. The solutions showed that the influence of temperature on the deformation-stress state around the CAES cavern is of importance when considering the stress state at a distance of less than 60 m from the cavern axis (at cavern diameter 30–35 m). With an increase in cavern diameter, it is possible that the impact range will be proportionately larger, but each case requires individual modeling that includes the shape of the cavern and the cavern working cycle.


Introduction
The stochastic process of the operation of power plants using renewable energy sources leads to the necessity to develop technologies that allow energy storage when there is no large demand for electricity. Therefore, to increase the use of renewable energy sources in the total energy balance, the most advantageous method turns out to be technologies which use a combination of both conventional fuels and renewable energy sources for the production of electricity.
CAES technology (Compressed Air Energy Storage) uses underground salt caverns to store compressed air. At times of increased demand for electricity in the electrical grid, a gas power plant is started up and the air stored in the cavern is supplied to the gas combustion process in the turbine.
There are two major industrial CAES installations in the world. One is the Huntorf power plant in Germany, which has been working since 1978 and originally had electrical power of 290 MW (currently 310 MW). The other is the McIntosh power plant in the United States (110 MW), operating since 1991.
The advantage of a power plant connected to a CAES system is that all electricity produced in the turbine can be transferred directly to the electrical network and the possibility of storage means that the surplus of energy produced does not have to be transferred to the network at the time of its generation.
When analyzing the work of a salt cavern (the storage part of a CAES installation), it should be noted that it operates in a relatively small pressure range compared to natural gas storage caverns (in the Huntorf power plant, it is in the range 5-7 MPa).
p max = g f raq ·h cem (1) where: • Cooling of the rock mass during leaching which is caused by the fact that the temperature of the water used for leaching is lower than the temperature of the rock mass; • Thermo-hydrodynamic changes in the air injected into the cavern; • Heat exchange between the stored medium, the rock mass and the pipes used for exploitation; • Heat propagation in the rock mass due to thermal conductivity; • Heat convection in the cavern.
However, in practice it is difficult to directly measure the temperature in the cavern during its operation or in the neighboring rock mass (theoretically, it is possible to install sensors inside the cavern, but due to difficulties in connecting to a recording device, this is not done). The phenomena occurring in the cavern can be determined by measuring pressure, temperature and flow capacity in the wellhead.
It is also very difficult to measure the heat transfer coefficients and precisely determine the cavern surface through which heat is exchanged with rock mass. Therefore, as an alternative method of obtaining the necessary information, it is proposed to use a computer model that will properly simulate the necessary parameters during the operation of the salt cavern.
Using the different models of creep law that are discussed in the literature, numerical calculations were performed of the deformation-stress state around a CAES salt cavern, and thermal influence was taken into account [5,6].
However, during laboratory tests of rock salt samples from Polish locations, it was noticed that determining the creep law coefficients for even such a simple law is a relatively difficult issue, even on the basis of tests that are performed for a specific single cavern. It is practically impossible to define a formula with a dozen or so coefficients. Therefore, it is recommended to use a simple Norton formula to describe the salt-creep phenomenon [7]. Therefore, the numerical creep law calculations in the model implemented Norton's law; these were performed to verify the extent of the influence of temperature changes on the deformation-stress state around the storage cavern. These results were compared [8], which showed that the influence of temperature on the deformation-stress state is important at a distance of less than 50 m from the cavern axis (for a cavern diameter of 30-35 m). At a greater distance, the models begin to converge [8].

Thermo-Mechanical Properties of Salt
The location of storage facilities in salt caverns is associated with finding optimal solutions and ensuring long-term stability of the rock mass and the tightness of the storage structure, while maintaining economic profitability. From the geomechanical point of view, the design consists of determining the depth of the cavern foundation and the maximum and minimum working pressure that can be achieved; the geometrical dimensions take the stored medium into account. Design calculations are carried out with the use of specialized computer programs that analyze the state of stresses in the rock mass. Unfortunately, in Energies 2021, 14, 6197 3 of 28 practice the results of these calculations are different than the actual values, as a result of the discrepancy between the parameters determined in laboratory tests and the actual behavior of the rock mass. The influence of temperature changes the heterogeneity of the salt deposit; the complexity of the equations used are also important [7].

Thermodynamic Model
The flow and energy of air mass in a CAES installation should be considered a complex process. The application of simple laws of thermodynamics regarding temperature and pressure changes in the cavern allows the behavior of the storage cavern to be predicted with an acceptable level of accuracy. To obtain an analytical solution in this part of the installation (salt cavern) and to determine changes in temperature and pressure, it is suggested to adopt the following assumptions [9][10][11]: • For thermodynamic problems, the cavern is considered as one volume ( Figure 1); • The compressibility of rock salt is very low compared to the compressibility of gas; therefore, it is assumed that the cavern volume is constant throughout the simulation (the change in cavern volume is much smaller than the total volume, therefore its influence on temperature and air pressure is negligible) [12]. When formulating the mass and energy balance equations in the first stage of the analysis, a constant volume V for the cavern is assumed. In the second stage of the analysis, in which the cavity stability is tested, the volumetric convergence is calculated to ensure that the applied thermo-mechanical loads meet the design requirements and stability is maintained; • Air flow during injection and withdraw is turbulent, so the air temperature and pressure gradients are very small [13,14]. It is assumed that air temperature, pressure, and air density are uniform throughout the cavern volume; • Injection and download air characteristics (temperature, enthalpy, and mass flow rate) are controlled by the CAES installation operator. Therefore, it can be assumed that the air properties are known during the injection and download processes to and from the cavern; • Air is considered an ideal gas-this assumption is true for the pressure and temperature inside the CAES storage cavern [10].

Thermo-Mechanical Properties of Salt
The location of storage facilities in salt caverns is associated with finding optimal solutions and ensuring long-term stability of the rock mass and the tightness of the storage structure, while maintaining economic profitability. From the geomechanical point of view, the design consists of determining the depth of the cavern foundation and the maximum and minimum working pressure that can be achieved; the geometrical dimensions take the stored medium into account. Design calculations are carried out with the use of specialized computer programs that analyze the state of stresses in the rock mass. Unfortunately, in practice the results of these calculations are different than the actual values, as a result of the discrepancy between the parameters determined in laboratory tests and the actual behavior of the rock mass. The influence of temperature changes the heterogeneity of the salt deposit; the complexity of the equations used are also important [7].

Thermodynamic Model
The flow and energy of air mass in a CAES installation should be considered a complex process. The application of simple laws of thermodynamics regarding temperature and pressure changes in the cavern allows the behavior of the storage cavern to be predicted with an acceptable level of accuracy. To obtain an analytical solution in this part of the installation (salt cavern) and to determine changes in temperature and pressure, it is suggested to adopt the following assumptions [9][10][11]: • For thermodynamic problems, the cavern is considered as one volume ( Figure 1); Figure 1. Considered parameters during the injection and download processes from a cavern [5].
• The compressibility of rock salt is very low compared to the compressibility of gas; therefore, it is assumed that the cavern volume is constant throughout the simulation (the change in cavern volume is much smaller than the total volume, therefore its influence on temperature and air pressure is negligible) [12]. When formulating the mass and energy balance equations in the first stage of the analysis, a constant volume V for the cavern is assumed. In the second stage of the analysis, in which the cavity stability is tested, the volumetric convergence is calculated to ensure that the applied thermo-mechanical loads meet the design requirements and stability is maintained; • Air flow during injection and withdraw is turbulent, so the air temperature and pressure gradients are very small [13,14]. It is assumed that air temperature, pressure, Figure 1. Considered parameters during the injection and download processes from a cavern [5].
Taking into account the above assumptions as well as the general principles of mass and energy balance, a mathematical formulation of the injection process and the air download from the storage cavern is obtained [5].

Cavern Operating Parameters
When performing thermodynamic calculations for the cavern, it should be assumed that the pressure in the reservoir is equal to the current pressure behind the compressor because it increases during injection due to the constant volume. A constant stream of air Energies 2021, 14, 6197 4 of 28 mass is maintained by regulating the speed of the compressor. During the calculations, a so-called "air cushion", i.e., the amount of air remaining in the cavern, is used to ensure the minimum pressure.
However, the fact that the pressure in the cavern before injection is equal to the pressure after expansion (emptying) does not mean that the temperature is also the same [15]. From the thermodynamic point of view, the storage cavern can be considered a reservoir from which there is an outflow through the nozzle due to the stability of the turbine parameters, assuming that an adiabatic process occurs.
It should be clearly emphasized that the purpose of the performed calculations was not to optimize the efficiency of the CAES installation but to maintain the safety of the storage part (the tightness of the cavern and the borehole) and the stability of the cavern during the cyclic operation of the storage facility. Therefore, the simulations focused only on the storage part (cavern). For the numerical calculations, the most extreme working scenario was selected to test the stability of the rock mass around the cavern. Therefore, a situation in which the efficiency of the installation would be reduced was not analyzed; it is expected that-since the stability of the cavern was maintained in the most extreme conditions-lower loading of the CAES installation will have a positive effect on the stability of the rock mass around the cavern.

Balance of Mass and Energy in the Cavern for Steady-State Conditions
Assuming a constant cavern volume V, the equation of the mass flow balance during injection and download takes the form: where: V-cavern volume (m 3 ); . ρ air -air density change in the cavern (kg/m 3 s); . m inj , . m dow -mass flow during injection/download (kg/s). The mass flow rates depend on the technical requirements of the installation (e.g., maximum compressor capacity or maximum air supply to the gas turbine) as well as the energy demand from the electrical grid.
The operator of the CAES installation can adjust the properties of the air inside the cavern to the desired range by controlling only these two factors.
Thus, the energy balance for the cavern takes the form: where: ρ air -air density in the cavern (kg/m 3 ); u-internal energy (J); h dow , h inj -air enthalpy during download and injection (kJ/kg); U-heat transfer coefficient between the cavern wall and the rock mass W m 2 ·K ; A cav -the area between the cavern and the rock mass (m 2 ); T cav , T rock -air temperature in the cavern and rock mass temperature (K). The last term on the right side of Equation (2) represents the heat exchange between the cavern and the salt rock mass.
If air is considered an ideal gas, the following equations can be used to calculate the pair pressure, internal energy u, and air enthalpy h. .
where: T inj -temperature of the injected air (K). The heat transfer coefficient U and the A cav area are values that cannot be easily determined [9,10], therefore the relationship U A c V should be replaced by the effective heat transfer coefficient [10], which is a function of the air mass flow rate during injection and the downloading of compressed air from the cavern: where: a 1 , a 2 , a 3 -parameters obtained during the calibration model with data from in-situ measurements [5].

PVT Equation for the State of Gas in a Cavern
In classical numerical models, air is usually treated as an ideal gas (Clapeyron equation): where: p-pressure (MPa); T-temperature (K); v-molar volume m 3 mol ; n-number of moles of gas (-). The general equation of state for the real gas is used in the calculations: where: Z-compressibility factor (-). The Z coefficient of compressibility is a function of the gas state in question, depending on the composition, temperature, and pressure. The model uses the Redlich-Kwong and Peng-Robinson real gas equations of state as well as other simplified PVT correlations.

Heat Exchange between Cavern and Rock Mass (Macroscopic Balance)
The isochoric cavern model (with the assumption that the variability of the stored air volume can be limited) can be described by the equation: where: Energies 2021, 14, 6197 6 of 28 κ-air adiabatic exponent (-); T A cav -temperature on the surface A cav (in the rock mass) (K).
The dynamic air temperature model can be described by the equation: where: m-air mass (kg); C rock -specific heat of the rock mass J kg·K [16]. However, the heat exchange between the compressed air inside the cavern and the rock mass can be described by the equation: (13) where: α-heat transfer coefficient W m 2 ·K ; λ rock -thermal conductivity coefficient in the rock mass W m·K ; dS-surface element; n-a unit vector normal (which points outward) to the surface of a cavern or gas-filled well (-); Ω-any fragment of the surface of the reservoir (-). The shape of the salt cavern cannot be described as an ideal sphere or cylinder because its side surface, due to its irregular shape, is much larger than that of a sphere with an analogous volume, hence the theoretical coefficient α should be treated as an approximation. In practice, inaccurate knowledge of the heat exchange has a relatively small effect on the cavern temperature. When using high coefficients (above 20 W m 2 ·K ), the heat inflow from the salt rock is limited by its thermal conductivity.

Temperature Changes in the Salt Rock Mass-Thermal Conductivity
Thermal conductivity in the rock mass is realized by the equation that compares the heat change in a given volume to the heat conduction balance through the side surface of this volume: ∂ ∂t Ω ρ rock C rock T rock dV = Γ λ rock ∇T rock ·ndS (14) where: ρ rock -rock mass density (kg/m 3 ); T-the surface of Ω element (m 2 ). This equation can determine the temperature distribution in the area on the edge of which the temperature or its gradient is known. Considering the rock mass in the vicinity of the cavern, the range of the model is so large that that the thermodynamic changes taking place in the cavern do not affect the conditions on the edge of the model. It can be assumed that the temperature at the edge of the considered model (in the rock mass at a great distance from the cavern side) does not depend on time but only on depth, according to the geothermal degree: where: h-depth below the ground surface (m) [17].

The Constitutive Equation of Rock Salt
A characteristic feature of rheological media, including rock salt, is the influence of loading time on the amount of deformation ( Figure 2). When salt samples are compressed, an increase in strain is observed at a given stress level as a function of the loading time.
to the geothermal degree: where: h-depth below the ground surface (m) [17].

The Constitutive Equation of Rock Salt
A characteristic feature of rheological media, including rock salt, is the influence of loading time on the amount of deformation ( Figure 2). When salt samples are compressed, an increase in strain is observed at a given stress level as a function of the loading time. Temperature, humidity, and grain size are also additional factors which may affect the deformation of the tested salt samples. Due to the rheological properties of the salt, even with constant loading of the samples, the deformation increases with time [7].
An important phenomenon in laboratory tests of salt loading is the so-called creep and salt relaxation (a decrease in stresses in a state of deformation that is constant over time).
In the description of the creep phenomenon, four phases should be distinguished ( Figure 3): 1. The initial deformation phase occurs so quickly that time is usually omitted in its description; 2. Primary creep, when creep speed decreases; 3. Stationary creep, characterized by a constant increase in strain over time. This phase has the greatest impact on the strain-stress state of the rock mass around the cavern; 4. Final phase. Under the influence of stresses, creep accelerates rapidly, which leads to material destruction [7]. Temperature, humidity, and grain size are also additional factors which may affect the deformation of the tested salt samples. Due to the rheological properties of the salt, even with constant loading of the samples, the deformation increases with time [7].
An important phenomenon in laboratory tests of salt loading is the so-called creep and salt relaxation (a decrease in stresses in a state of deformation that is constant over time).
In the description of the creep phenomenon, four phases should be distinguished ( Figure 3):

1.
The initial deformation phase occurs so quickly that time is usually omitted in its description; 2.
Stationary creep, characterized by a constant increase in strain over time. This phase has the greatest impact on the strain-stress state of the rock mass around the cavern; 4.
Final phase. Under the influence of stresses, creep accelerates rapidly, which leads to material destruction [7].
The creep speed in elastic-viscous media depends not only on their deformation-stress state but above all on the state of knowledge about them; this causes great difficulties in describing the properties of salt. In practice, a large number of simplifications are used [7].
To describe the properties of salt, it is assumed that it is an isotropic medium, i.e., the properties of the medium are independent of the direction in which they are measured. In addition, it is assumed that the change in the deformation state is a function of the current independent variables and depends on stress, strain, temperature and time [7].
Another important assumption is the application of the superposition principle for deformation; this principle determines the total increase in deformation as the sum of plastic, elastic and viscous deformations. The description of elastic deformation is based on Hooke's law, and the influence of temperature is also considered. Although the influence of plastic deformations on the rock mass in the vicinity of the storage cavern is smaller than that of elastic deformation, the theory of their description is essential because theories of plasticity are built on this basis. The description of plastic deformation includes the theory of plastic flow and the theory of plasticity. Viscous deformations require the assumption of zero creep dilatancy [7].  The creep speed in elastic-viscous media depends not only on their deformationstress state but above all on the state of knowledge about them; this causes great difficulties in describing the properties of salt. In practice, a large number of simplifications are used [7].
To describe the properties of salt, it is assumed that it is an isotropic medium, i.e., the properties of the medium are independent of the direction in which they are measured. In addition, it is assumed that the change in the deformation state is a function of the current independent variables and depends on stress, strain, temperature and time [7].
Another important assumption is the application of the superposition principle for deformation; this principle determines the total increase in deformation as the sum of plastic, elastic and viscous deformations. The description of elastic deformation is based on Hooke's law, and the influence of temperature is also considered. Although the influence of plastic deformations on the rock mass in the vicinity of the storage cavern is smaller than that of elastic deformation, the theory of their description is essential because theories of plasticity are built on this basis. The description of plastic deformation includes the theory of plastic flow and the theory of plasticity. Viscous deformations require the assumption of zero creep dilatancy [7].
The rheological properties of the medium are represented by the constitutive equation given in the form of two equations. The first describes the change in volume and relates the stress and strain axioms. The second is the deviator equation that describes the law of change of form [7].
where: Pσ, Qε-differential operators; Dσ-stress deviator (MPa); Dε-deformation deviator (‰). In special cases, the deviator equation may assume a linear form with constant coefficients. Then, the components of the deviators of both stresses and strains are related by the following relationship [7]: where: The rheological properties of the medium are represented by the constitutive equation given in the form of two equations. The first describes the change in volume and relates the stress and strain axioms. The second is the deviator equation that describes the law of change of form [7].
where: P σ , Q ε -differential operators; D σ -stress deviator (MPa); D ε -deformation deviator (‰). In special cases, the deviator equation may assume a linear form with constant coefficients. Then, the components of the deviators of both stresses and strains are related by the following relationship [7]: where: Various models are used to describe the properties of a rheological medium, e.g., Kelvin, Maxwell, Poynting-Thomson, and Prager. However, the rheological model of rock salt is actually much more complicated. The coefficients of the constitutive equation are not constant parameters as they depend on the stresses and side factors. The complexity of the models used to describe salt is mainly from taking into account various types of salt-creep deformations. The final form of the rheological equation is primarily determined by the manner in which the properties of a complex medium are presented as a set of properties of simple media [7]. As was mentioned in the introduction, during the laboratory tests of rock salts it was noticed that determining the creep law coefficients for even such a simple law is a relatively difficult issue, therefore it is recommended to use a simple Norton's formula to describe the salt-creep phenomenon [7].

Cavern Stability Parameters
Geomechanical problems in terms of the stability of the storage cavern as a result of its exploitation should be considered in five categories [18]:

Local Stability of the Cavern
There are several criteria to be met by the cavern in order not to damage the salt rock on the side walls and the ceiling:

•
Lack of tensile stresses. This criterion results from the very low tear strength of the salt, which means that in the evaluation of the stress state it is assumed to be zero. If this criterion is met, an appropriate vault is created during the leaching of the cavern. • Long-term durability is not exceeded. This is one of the most discussed criteria in geomechanics. In the case of rock salt, it is unacceptable for the rock mass to enter the third phase of creep because this leads to rapid destruction. In order to determine the permissible stress state, strength hypotheses and rock mass effort are introduced as these determine the value of the analyzed stress state in relation to the short-term strength. The allowable values of the effort are determined subjectively depending on the required level of security.

•
The permissible deformations are not exceeded. Taking into account the ambiguity of the stress and deformation states and the variable pressure inside the working storage cavern, it is necessary to determine the permissible effective deformation state. This criterion determines the allowable increase in effective deformation during pressure reduction and at minimum pressure, as well as the allowable rate of effective deformations.

Global Cavern Stability
The assessment of this criterion consists of examining the extent of the impact of the cavern on the rock mass. This criterion is determined on the basis of the range of stresses and effective deformations as well as the effort of the rock mass, assuming that the range of significant changes in these quantities determines the extent of the impact of the cavern.
It should be noted that the theoretical extent of the impact of the storage cavern on the rock mass depends primarily on the constitutive equation adopted for the description of the rock mass. With high rheological properties, the range should be greater.

The Tightness of the Cavern
The tightness criterion is related to the maximum storage pressure in the cavern and is determined on the basis of microfracturing tests carried out in the well during its drilling. On the basis of these tests, the so-called "resting pressure of fracturing", which is established after fracturing the rock mass and is close to the primary pressure in the rock mass, increases proportionally to the depth. The proportionality factor is called the microfracture gradient.
The stress distribution around the cavern when the stored medium (air) is at high pressure differs from the primary pressure of the rock mass. The migration of the stored medium is therefore possible at a pressure lower than that which results from the microfracture gradient. When determining the maximum storage pressure, an appropriate safety factor is adopted.
The issue of the tightness of the cavern is also related to its maximum allowable height. If it is too high, it is not possible to empty the brine from the cavern due to cracking of the Energies 2021, 14, 6197 10 of 28 rock mass in the ceiling. The permissible maximum height of the cavern is determined by the formula [7]: where: H cav -active height of the cavern, including the "neck"; H cem -cementation depth (m); φ s -resting microfracture gradient (MPa/m); x-brine flow resistance per 1 m of pipe (-); γ brine -specific gravity of brine (kg/m 3 ). This condition is essential in the case of shallow-set storage caverns.

Cavern Capacity Changes during Storage Operations
The loss of the cavern's active capacity due to convergence is one of the most important issues that arise during the exploitation of storage caverns in rock salt deposits. This phenomenon is inevitable, and its intensity depends on the depth of the foundation, rheological properties, as well as the operating pressure, in particular the minimum pressure.
At constant pressure, the rate of convergence stabilizes and is a value that can be approximated by the formula [1]: where: p ∞ (h center )-primary rock mass pressure at the depth of the cavern center (MPa); A 1 , m 1 , Q 1 , B 1 -coefficients; η-cavern slenderness (-); σ 0 = 1 MPa. Since the pressure in the cavern is not constant during the operation of the storage installation, the rate of convergence changes over time. Generally, when describing the storage operation scenario, it can be divided into four stages, and it can be assumed that:

•
The cavern remains at pressure p min for duration ∆t min ; • The cavern is filled to pressure p max for duration ∆t inj ; • The cavern remains filled to pressure pmax for duration ∆t max ; • The cavern is emptied to pressure p min for duration ∆t dow .
Then, the average annual rate of convergence can be described by the formula [1]:

Protection of the Site Surface
Surface settlement in the case of storage caverns is a marginal issue due to the large foundation depths of the caverns and the distance between them. Despite this, regular leveling measurements are carried out; however, as practice shows, the possible land subsidence is much smaller than in the case of the classic exploitation of rock salt using the borehole method [1,18].

Assumptions for the Assessment of the Influence of Temperature Variability on the Deformation-Stress State of Salt Rock
Based on the parameters of the original installation in Huntorf, certain design assumptions were made as a result of technical and economic conditions that are mainly related to the strength of production pipes, casing pipes and cementation under significant temperature fluctuations. The parameters of the Huntorf power plant indicate that the well productivity exceeds 1 million m 3 /h, which requires the use of 24 diameter production pipes. An increase in the depth of the cavern foundations would significantly affect the cost of piping and compressors and would increase the flow resistance; therefore, it would reduce the overall efficiency of the power plant.

Analysis of the Injection and Download Cycles of the Cavern
Taking as a reference the parameters of the Huntorf power plant [2], which has a power of 290 MW, a cavern volume of 310,000 m 3 , and which can be emptied in 2 h, the active volume of cavern V for generating power P should be [4]: where: V-volume in thousand cubic meters; t-operating time in hours; P-power in MW.
The number of caverns n to be leached to reach the required volume V: This depends on two factors: the deposit thickness, which limits the cavern volume [4]: where M is the thickness of the salt deposit; limitations resulting from the possibility of emptying the cavern in the required time: The number of caverns, n, is determined by the smaller of these values. Based on the above equations, it is possible to estimate the number of caverns that should be made in order to obtain the required power in a deposit of a certain thickness. Assuming the thickness of the rock salt deposit is approx. 200 m, two caverns (cavern height of 150 m) should be leached to achieve a capacity similar to that of the Huntorf power plant (300 MW) and to deliver energy for 2-3 h. However, if the energy delivery time were extended to 8 h, the number of caverns needed would increase to 4.

Preferred Range of Operating Pressure
The range of pressures at which a storage cavern operates depends primarily on the depth of its foundation, the geological structure of the deposit, and the stored medium, as well as the operating scenario. For this reason, the parameters of the CAES installation in Huntorf were adopted for numerical simulations. The maximum operating pressure will be close to the maximum operating pressure of the system, i.e., 66 bar. The minimum pressure (the pressure to which the caverns in the Huntorf are emptied) will be 46 bar. On this basis, the operating pressure amplitude in the simulated cavern was 20 bar.

Temperature Variation Range
For the simulation it was decided to assume that the volume of the cavern was close to the volume of the existing caverns in Huntorf, i.e., 150,000 m 3 : the cavern ceiling was at a depth of 650 m; the bottom of the cavern was at a depth of 800 m. The assumed working scenario was continuous operation of the CAES storage facility. The following parameters were adopted for the simulation:

•
Daily cyclic emptying and filling-1 daily cycle consists of two hours of emptying the cavern (download); 6 h of resting the cavern at the minimum pressure, 8 h of filling (injection); and 8 h of resting the cavern at the maximum pressure.

•
With these assumptions, a thermodynamic simulation of CAES cavern operation was carried out. On this basis, the temperature distribution in the cavern was obtained. Then, the temperature distribution in the cavern was included in the further geomechanical calculations using a modified version of the GEOSOLK program, which considers the thermo-mechanical effects, including the temperature distribution around the storage cavern.

•
For proper simulation in the GEOSOLK program, it was assumed that after the leaching process was completed, the cooling of the rock mass would stabilize the temperature in the cavern at 36.0 • C.

•
For the thermodynamic simulations, the minimum temperature in the cavern (in the first emptying to 46 bar) was 12.28 • C. After 365 daily operating cycles, it dropped to the minimum value: 11.19 • C.

•
The maximum temperature in the cavern after filling to the maximum pressure was 41.54 • C at the beginning of the simulation (after the first filling); it stabilized at 40.23 • C after filling and the simulation was complete after 365 daily work cycles.
The preliminary calculations showed that the maximum range of temperature changes (the difference between the temperature in the cavern after emptying and the temperature after filling) did not exceed 30 • C.

Numerical Model of CAES Salt Cavern Operation
The aim of the simulations was to determine the impact of changes in the temperature of the stored compressed air on the strain-stress state of the rock mass directly in its vicinity. The simulations included the range of temperature changes and load changes that corresponds to the operated CAES cavern. The temperature distribution around the cavern during operation was also made.
Numerical simulations of the stress-strain state of the rock mass during the exploitation of the CAES cavern were performed using the GEOSOLK software. Initially, the calculations were performed with the assumption that the original temperature in the salt deposit did not change over time. Then, the calculations were performed with the assumption that the temperature did change over time, and the temperature distribution in the rock mass was caused by the variable air temperature in the cavern and the conditions on the edges of the model: (1) Constant primary temperature on the vertical edge of the model; (2) The heat flux resulting from the geothermal gradient on the horizontal edge.
To compare the results of both simulations, the influence of temperature changes on the stress-strain state of the rock mass was assessed.
The simulations were performed for two sets of parameter values of the adopted "creep law". The scenario of a cyclical CAES cavern operating for 365 days was analyzed.
During the assessment of cavern stability, the stresses and effective deformations were used [8]. An analysis of the distribution of the rock mass effort around the cavern in the selected model elements at the edge of the cavern was also carried out.
Since no single simple criterion has been created so far to define the permissible stress and strain limit values in the cavern vicinity, several criteria describing the effort of the rock mass in the vicinity of the cavern are instead used for this. In this research, three criteria were used (described in Section 5) [4].
The effort values obtained according to these criteria were compared with the values obtained using the formulas developed by Lux, Hou and Düsterloh [18]. Figure 4 shows a diagram of the working CAES cavern modeling procedure, taking into account the changes in air pressure and temperature during download from and injection into the storage. The models and criteria described below were adopted for the thermo-mechanical analysis of the deformation-stress state. The effort values obtained according to these criteria were compared with the values obtained using the formulas developed by Lux, Hou and Düsterloh [18]. Figure 4 shows a diagram of the working CAES cavern modeling procedure, taking into account the changes in air pressure and temperature during download from and injection into the storage. The models and criteria described below were adopted for the thermo-mechanical analysis of the deformation-stress state.

Strength Criteria Adopted for Salt
The strength criteria commonly used for rock salts occurring in Poland were adopted for the numerical calculations. The GEOSOLK software's strength criteria file contains seven parameters. The values of these parameters were as follows: Rc = 20.

The Condition of Viscoelasticity-the Creep Law
In the adopted numerical model, Norton's creep law, taking into account elastic strains and viscous strains, was used to describe the viscoelasticity condition [7].
where: -total increase of strains; -increase of elastic strains;

Strength Criteria Adopted for Salt
The strength criteria commonly used for rock salts occurring in Poland were adopted for the numerical calculations. The GEOSOLK software's strength criteria file contains seven parameters. The values of these parameters were as follows: R c = 20.

The Condition of Viscoelasticity-the Creep Law
In the adopted numerical model, Norton's creep law, taking into account elastic strains and viscous strains, was used to describe the viscoelasticity condition [7].
where: dε-total increase of strains; dε el -increase of elastic strains; dε ν -increase of viscous strains. To compare salts with different rheological properties, the coefficients that describe the viscous strains part of the equation were modified in the numerical model.
Additionally, the adopted assumption of deformation superposition allows for a separate description of each type of deformation. The description of elastic strains is based on Hooke's law and the assumption of salt isotropy, so we can write [7]: where: ε el ij -components of the strain tensor; σ ij -components of the stress tensor;  The speed of the creep tensor components is expressed as [7]: where: However, the unknown in the above equation is the speed of effective deformation related to the creep dε ν e f dt . The creep law is a formula that defines the speed of effective creep strains as a function of effective stresses, effective strains, time and temperature.
Based on the results of laboratory tests and literature studies, it was decided to use Norton's creep law while creating the numerical models [7].
where: Q-activation energy (J/mol); R-universal gas constant J mol·K ; T-temperature (K); A, n-constant coefficients. To be able to compare the values obtained in the two models that depend on the salt-creep speeds, two values of the A coefficient were adopted so that the laws used for calculations significantly differed in terms of the creep speed.
For the calculations, the Q/R ratio = 6500 K was assumed, which means that at 50 • C the salt-creep rate will be almost eight times faster than at 20 • C. For this value of the Q/R coefficient, two values of the power coefficient were adopted, n = 4 and n = 5, because this range probably includes the actual value of this coefficient for the salt rock mass. The values of Norton's creep law coefficients adopted in the models are summarized in Table 1.

Thermal Properties of Salt
It was assumed that the primary temperature of the salt rock mass changes linearly [19]: where: α rock -temperature gradient; β-constant ratio; H-depth (below ground level). The value of the temperature gradient varies according to the salt structure under consideration; therefore, for the purposes of computer simulations, the average values of the parameters related to the thermal properties of salts are summarized in Table 2. On the basis of the adopted temperature gradient, it was determined that the range of temperature changes in the rock mass in the interval 0-2000 m was 15-75 • C.

The Adopted Calculation Scenario and Model Parameters
The scenario used for the simulations was based on a working installation using CAES (Compressed Air Energy Storage) technology for energy storage in compressed air [4,20].
Due to the fact that there are very rapid pressure changes in an operating CAES installation (pressure drop by 2 MPa in two hours), one can expect equally rapid changes in temperature inside the cavern that should be taken into account in the calculations [2,21]. The adopted calculation scenario consisted of 365 cycles of cavern daily operation. One cycle includes two hours of downloading the compressed air from the cavern; 6 h of pause (the cavern stays under the minimum pressure); 8 h of injection of compressed air into the cavern; and 8 h of pause at maximum pressure. In the numerical calculations, the maximum pressure in the cavern is 6.6 MPa when the cavern is full, and the minimum pressure is 4.6 MPa when the cavern is empty [3,[22][23][24].
The selection of the most representative points at which the effort was assessed was based on the representativeness assessment performed in [1], which is based on analysis of the correlation of effective stress values at selected points. The compared solution variants showed that the relevant results are strongly correlated [1].
The convergence values calculated with GEOSOLK software in the previous simulations were verified with the calculations carried out in FLAC 2D software [1]. The performed calculations showed a high similarity of obtained results.
To determine the extent of the impact of temperature changes in the cavern on the rock mass, the preliminary studies were carried out for two models [8]: 1-With a constant distribution of rock mass around the storage cavern; 2-Taking into account the impact of changes in the temperature of the rock mass around the storage cavern.
The model grid consisted of 660 elements. The horizontal range of the cavern (measured from the cavern center) was 250 m, including 17.84 m in the cavern radius. The vertical range (measured from the cavern center) was also 250 m [8].
The obtained results showed the greatest differences directly on the cavern wall and in its very close vicinity. On the basis of the obtained results, it was concluded that the influence of temperature on the deformation-stress case in the vicinity of CAES caverns is important when considering the stress state at a distance of less than 50 m from the axis of the storage cavity (when considering storage caverns with diameters of 30-35 m) [8].
Therefore, the same cavern model was used in the calculations (radius 17.84 m and height 150 m) for the 650-800 m interval below ground level, which resulted in a cavern volume close to 150,000 m 3 . A pipe with an external diameter of 21 (internal diameter 20 ) was adopted for the cavern exploitation. However, compared to the preliminary model [8], the number of mesh elements was increased to 996. The number of elements was optimized so as to obtain satisfactory computer operation because no large differences in the obtained results were noticed in the simulations performed for other meshes (increasing the number of elements in the mesh did not cause large differences in the obtained results).
Tables 3-5 summarize the most important parameters of the cavern model on which the simulations were carried out.   Figure 5 shows the CAES cavern mesh, with selected elements marked at points on the ceiling, side wall and bottom (measured from the wall of the cavern into the rock mass). It should be noted that the shape of the cavern used for the simulation is different in the upper wall and in the lower wall; this is dictated by the SMRI (Solution Mining Research Institute) guidelines regarding the relationship between the shape of the dome and the diameter and height of the cavern.
The results presented below take into account the calculations for the two models: fast-(FCL) and slow-creep law (SCL), taking into account the impact of changes in the temperature of the rock mass around the storage cavern. Figure 6 compares the stresses in the rock mass with the applied creep law.  The results presented below take into account the calculations for the two models: fast-(FCL) and slow-creep law (SCL), taking into account the impact of changes in the temperature of the rock mass around the storage cavern. Figure 6 compares the stresses in the rock mass with the applied creep law. When comparing both the models used in the simulations (Figure 7), it was noticed that the applied creep law affects the values obtained up to a distance of 50 m from the cavern axis, regardless of the considered cavern point in a given case (ceiling, side wall, bottom). This translates into the real range of cavern impact differences, depending on the When comparing both the models used in the simulations (Figure 7), it was noticed that the applied creep law affects the values obtained up to a distance of 50 m from the cavern axis, regardless of the considered cavern point in a given case (ceiling, side wall, bottom). This translates into the real range of cavern impact differences, depending on the type of salt (fast-creeping or slow-creeping) up to a maximum of 23-24 m from the cavern wall (the considered model has a radius of 17.84 m). It should be noted that in slow-creeping salts the intermediate stresses in the rock mass in the immediate vicinity of the cavern are about 3 MPa higher than in fast-creeping salts. When analyzing the results obtained for the deformation in model elements, it was noticed that the situation is opposite to the situation for stresses, i.e., in the fast-creeping model we deal with greater deformation in the immediate vicinity of the cavern than in the slow-creeping model. However, similar to the stresses, at a distance of approx. 23 m from the cavern axis, the strains begin to assume similar values. This means that in the case of fast-creeping salts, we deal with the greater phenomenon of stress relaxation.
It is worth noting that when analyzing the values of deformations in the model elements of the ceiling and the bottom, at a distance of 70 m from the cavern wall these deformations stabilize at the level of 0.00-0.01‰, which means that they are practically negligible. However, deformations of 0.01‰ in the side wall elements do not appear until 170 m from the cavern wall.

Strength Criteria and Rock Mass Effort
Researchers from Germany use the most precise criteria for assessing the stability of storage caverns. It is true that there are no unambiguously applicable limit values regarding the stress and deformation states in the cavern vicinity, but the three detailed criteria that were formulated and assessed by Lux et al. [18] were used to evaluate the stability of the cavern.
The first criterion concerns the maximum rock mass effort [18,25,26] and assumes When analyzing the results obtained for the deformation in model elements, it was noticed that the situation is opposite to the situation for stresses, i.e., in the fast-creeping model we deal with greater deformation in the immediate vicinity of the cavern than in the slow-creeping model. However, similar to the stresses, at a distance of approx. 23 m from the cavern axis, the strains begin to assume similar values. This means that in the case of fast-creeping salts, we deal with the greater phenomenon of stress relaxation.
It is worth noting that when analyzing the values of deformations in the model elements of the ceiling and the bottom, at a distance of 70 m from the cavern wall these deformations stabilize at the level of 0.00-0.01‰, which means that they are practically negligible. However, deformations of 0.01‰ in the side wall elements do not appear until 170 m from the cavern wall.

Strength Criteria and Rock Mass Effort
Researchers from Germany use the most precise criteria for assessing the stability of storage caverns. It is true that there are no unambiguously applicable limit values regarding the stress and deformation states in the cavern vicinity, but the three detailed criteria that were formulated and assessed by Lux et al. [18] were used to evaluate the stability of the cavern.
The first criterion concerns the maximum rock mass effort [18,25,26] and assumes that the maximum rock mass effort W ef at the edge of the cavern is [7]: σ max e f (σ m , θ, σ 3 ) = 1 cos θ + π 6 + asin θ + π 6 exp(bσ 3) σ TC e f (σ m ) where: a, b-constant coefficients (-); σ TC e f (σ m )-strength in conventional triaxial compression tests (MPa). At stress σ 3 , which is equal to the minimum gas pressure in the cavern, W ef cannot exceed 50% as this may result in loss of cavern stability.
The second criterion [18,25,26] concerns the effective stresses inside the interventricular pillar at the depth of the maximum stress at the wall. The effort in the distance L 2 − D from the cavern axis (where L is distance between cavern axes, D is cavern diameter) may not exceed 30%.
The third criterion [18,25,26] concerns the maximum effective deformation at the cavern wall; it is assumed that the deformation increment during the pressure drop and the period when the cavern is under minimum pressure may not exceed 3%.
However, the precision of the short-term strength description, which depends on the degree of complication of the strength hypothesis, applies to typical salts, for which there is an extensive database of measurements and experiments. Therefore, when we have a limited amount of data, calculations using less complicated formulas turn out to be better. At the same time, it should be noted that the problem of long-term strength is still unsolved, hence it is recommended to assume a subjectively acceptable stress effort degree of the rock mass around a particular cavern within 30-40% of the short-term strength.
The basic parameters used in this study to analyze the stability of a cavern in a working CAES system were the effective stresses and strains. In addition to these two values, analysis of the rock mass around the cavern was also carried out in three model elements at the cavern walls: in the ceiling, on the side wall and at the bottom (respectively, elements 816, 504, 192 of the cavern mesh shown in Figure 5). The scenario of cyclic cavern operation lasting 365 days, as was described in Section 4, was analyzed. The effective stress hypothesis obtained by applying the formulas developed by Lux et al. [18] was used to assess the effort; three strength hypotheses are commonly used for these purposes: Maximum principal stress hypothesis [7]: [27] with parameters: R c = 20.6 (MPa) a = 17.55 (-) b = 0.692 (-) Burzynski's hypothesis [7]: where: σ m = 1 3 (σ 11 + σ 22 + σ 33 ) (MPa) with parameters: R c = 20.6 (MPa) α = 0.025 (-) Huber's hypothesis [7]: In this case, based on the hypothesis of maximum principal stresses, the values of effective stresses can be analyzed: with parameters as above [7].
The highest values of the stresses are achieved immediately after emptying the cavern to the minimum pressure ( Table 6). The intermediate values of the stresses are achieved after 6 h of cavern rest under minimum pressure and immediately after injection to the maximum pressure (Tables 7 and 8) The lowest values of effective stresses were observed after the maximum pressure in the cavern stabilized (after the cavern had rested for 8 h at maximum pressure) ( Table 9). Diagrams of effective stresses for the two creep laws (365-day cycles of cavern operation) are shown in Figure 8. Table 6. Strength criteria for the two salt-creep models. Status at the end of a given downloading cycle (after emptying the cavern to the minimum pressure).  It should be noted that the differences in the values of effective stresses in extreme cases of cavern operation do not exceed 1 MPa.

Time (Days) Pressure in Cavern (MPa) Model Type
When the cavern is continuously operated, the stresses in the considered elements decrease and, after the simulation of 365 daily cycles, their values are 3-4 MPa lower than the values determined after the first operation cycle. This was the case with both the fastand slow salt-creep factors. The effective stresses at each considered are lower when using the fast-creep law (FCL) than with the slow-creep law (SCL). This is due to the greater value of the n-factor used in this law, which means greater ability to reduce extreme stress. This situation occurs both when there is minimum and maximum pressure in the cavern; in both cases, the effective stress in the side wall element drops below 10 MPa.
In the case of the slow-creep law, the effective deformations in the elements on the side wall and on the bottom increase by a maximum of 0.4‰ during a year of operation of the cavern. With the use of fast salt-creep, these deformations significantly increased compared to the first model and reached values from 0.5-0.6‰ for the considered element in the ceiling to 1.3‰ for the element in the bottom (Table 10 compares the values listed in Tables 6-9). This situation was expected. The diagrams of effective strains obtained for the two creep laws (365-day cycles of cavern operation) are shown in Figure 9.  The stresses and effective strains are the most important in terms of cavern stability. Comparing the values obtained in the cavern wall elements against the stress values presented in the diagrams in Section 4, different ranges of cavern impact can be observed depending on which creep law was applied. It can be noticed that a greater reduction of stresses (fast-creep law) causes a greater increase in stresses in the rock mass, i.e., the range  The stresses and effective strains are the most important in terms of cavern stability. Comparing the values obtained in the cavern wall elements against the stress values presented in the diagrams in Section 4, different ranges of cavern impact can be observed depending on which creep law was applied. It can be noticed that a greater reduction of Figure 9. Comparison of the effective strains obtained for the slow-(SCL) and fast-creep law (FCL) models, with temperature changes in the rock mass for the considered elements on the cavern contour. The stresses and effective strains are the most important in terms of cavern stability. Comparing the values obtained in the cavern wall elements against the stress values presented in the diagrams in Section 4, different ranges of cavern impact can be observed depending on which creep law was applied. It can be noticed that a greater reduction of stresses (fast-creep law) causes a greater increase in stresses in the rock mass, i.e., the range of cavern impact depends primarily on the type of salt (fast-or slow-creeping).
Different stress values cause differences in the obtained values of the rock mass effort coefficients in the considered elements at the wall of the cavern. The highest effort values are achieved by the W ef coefficient: in the element on the side wall, at the stage immediately after emptying the cavern on the first day of operation, it is equal to 50% (for the slow-creep law) (Appendix A Table A1). However, in no case does its value exceed the limits allowed by Hunsche and Lux et al.: 50% [18,25,26]. It should be emphasized that the W ef criterion determines the maximum effort that can be achieved in the rock mass, hence it should be assumed that the local stability of a CAES cavern operating under the assumed conditions will be maintained.
When analyzing the remaining strength criteria, it can be noticed that the highest values of the parameters describing the remaining hypotheses, W H' , W σ max 1 , W B , occur on the cavern side wall in stages immediately after emptying the cavern to the minimum pressure and after the cavern rested under the minimum pressure for 6 h (Appendix A  Tables A1 and A2). These values are similar in both stages, which may indicate a lack of stress relaxation in the rock mass in such a short period of time. With the passage of time of the cavern operation, the coefficients describing the effort of the rock mass begin to decrease, which indicates the relaxation of stresses in the rock mass despite the continuous cyclic operation of the cavern between the maximum and minimum pressure. It should be emphasized that the obtained values of these coefficients are so small that there is no risk of not maintaining cavern stability due to the use of insufficient minimum pressure.
Depending on the stage in question in individual work cycles, comparison of the obtained values of effective stresses and deformations (Appendix A Tables A1-A4) allows the conclusion that the deformation-stress state of the rock mass around the storage cavern undergoes some stabilization during operation of the cavern. This tendency is clearly noticeable for the simulations the with slow-creep law.
The annual increments of the strains presented in Table 10 are definitely lower than the permissible value defined by Lux et al.'s criteria, i.e., 3% [18,24,25]. It should also be noted that for the considered models, the stresses in the vicinity of the cavern wall drop below 12 MPa, especially in the side wall, which may lead to salt rock penetration.

Temperature Distribution in the Rock Mass Around CAES Cavern
During the operation of the CAES installation, rapid temperature changes occur in the cavern due to the injection of compressed air at a given temperature (35 • C) and the download of large amounts of stored air at a temperature of 15.2 • C for a very short period of time (2 h). Very high downloading and injection efficiency causes large temperature amplitudes.
As a result of the numerical simulations, temperature distributions in the rock mass around the CAES storage cavern were calculated for the adopted scenarios. Figure 10 shows the maps of temperature distribution obtained at selected moments of the simulation, i.e., on the 1st, 30th, 90th, 180th and 365th day of operation. Temperature fluctuations in the rock mass over the course of one day were about 1 °C. This is caused by the high specific heat of the rock mass and very fast temperature changes in the cavern. The range of temperature changes observed on these maps did not exceed 50 m from the cavern axis.
These maps show that at the beginning of the CAES cavern operation, the average temperature in the rock mass directly on the side wall increases; then, as the operation proceeds, the average temperature in the rock mass at all analyzed points drops by 2-3 °C compared to the value at the beginning of the operation time.

Conclusions
The performed numerical calculations show that the salt-creep process influences the stress variability in a range of up to 50 m from the cavern axis (approx. 23-24 m from the cavern wall in the considered model), regardless of the cavern point (ceiling, side wall, bottom). It should be noted that in slow-creeping salts, the stresses in the rock mass in the immediate vicinity of the cavern are greater than in fast-creeping salts.
The deformations obtained for fast-creeping salts show greater values in the immediate vicinity of the cavern than those obtained for slow-creeping salts. However, at a distance of approx. 23 m from the cavern axis, they have comparable values. This means that in the case of fast-creeping salts, we deal with a greater phenomenon of stress relaxation.
The highest effort values are achieved by the Wef coefficient, but in no case does its value exceed the permissible value according to Hunsche and Lux et al.'s criteria, i.e., 50% [18,25,26].
The highest values of the parameters described by hypotheses WH', , and WB occur on the cavern face in stages immediately after the cavern is emptied to the minimum pressure and after it is rested under the minimum pressure. During the exploitation of the CAES cavern, the coefficients describing the effort of the rock mass begin to decrease, Temperature fluctuations in the rock mass over the course of one day were about 1 • C. This is caused by the high specific heat of the rock mass and very fast temperature changes in the cavern. The range of temperature changes observed on these maps did not exceed 50 m from the cavern axis.
These maps show that at the beginning of the CAES cavern operation, the average temperature in the rock mass directly on the side wall increases; then, as the operation proceeds, the average temperature in the rock mass at all analyzed points drops by 2-3 • C compared to the value at the beginning of the operation time.

Conclusions
The performed numerical calculations show that the salt-creep process influences the stress variability in a range of up to 50 m from the cavern axis (approx. 23-24 m from the cavern wall in the considered model), regardless of the cavern point (ceiling, side wall, bottom). It should be noted that in slow-creeping salts, the stresses in the rock mass in the immediate vicinity of the cavern are greater than in fast-creeping salts.
The deformations obtained for fast-creeping salts show greater values in the immediate vicinity of the cavern than those obtained for slow-creeping salts. However, at a distance of approx. 23 m from the cavern axis, they have comparable values. This means that in the case of fast-creeping salts, we deal with a greater phenomenon of stress relaxation.
The highest effort values are achieved by the W ef coefficient, but in no case does its value exceed the permissible value according to Hunsche and Lux et al.'s criteria, i.e., 50% [18,25,26].
The highest values of the parameters described by hypotheses W H' , W σ max 1 , and W B occur on the cavern face in stages immediately after the cavern is emptied to the minimum pressure and after it is rested under the minimum pressure. During the exploitation of the CAES cavern, the coefficients describing the effort of the rock mass begin to decrease, which indicates the relaxation of stresses in the rock mass despite the continuous cyclic operation of the cavern between the maximum and minimum pressure. It should be noted here that the values of these coefficients, which were obtained as a result of the simulation, meet the adopted criteria, thus confirming that the temperature variation in the scope of the compressed air storage operation significantly affects the deformation-stress state of the cavern side, but it does not cause a loss of its stability. It should be expected that with many years of exploitation of the cavern, the deformation and stress state of the rock mass around the storage cavern will stabilize.
Annual deformation increments are definitely lower than the permissible value defined by Lux et al., i.e., 3% [18].
It should be emphasized that the CAES installation works when the demand for electricity in the network increases. The scenario adopted for the calculations was intended to verify the stability of the cavern in the most extreme load situation, when cyclic filling and emptying of the cavern occurs at very short intervals. The simulations show that it is possible to use CAES caverns even in the most extreme mode of operation. At the same time, it should be remembered that each case of a cavern located in a salt deposit should be considered individually, therefore it is necessary to conduct research (including numerical simulations) for each given case.
Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicity available due to company's policy.