Displacement Solution of Salt Cavern with Shear Dilatation Behavior Based on Hoek-Brown Strength Criterion

In this study, an analytical solution of stress, strain, and displacement, in the elastic and plastic zone is proposed. *e solution is derived on the basis of ideal elastoplastic mechanical model of spherical salt cavern with shear dilatation behavior, by adopting Hoek-Brown (H-B) criterion. *e solution obtains not only in small and large strain stage but also in creep stage. *e proposed solution is validated, by comparison of the obtained results with numerical results in FLAC3D.*e results indicate that the result obtained adopting the H-B criterion is closer to that one obtained adopting the Mohr-Coulomb (M-C). *e H-B criterion is more applicable for the salt cavern construction as it considers the structural characteristics of the rock salt formation.*e displacement difference obtained by two different methods decreases with the increase of GSI or running pressure, but it increases with the enlarged angle of dilation. *e influence of different assumptions of elastic strain of plastic zone on displacements is more significant under large strain conditions. *e influence of the angle of dilation on displacements is more obvious when the elastic strain of plastic zone is given to stationary values, and the influence degree increases with the enlarged angle of dilation. Under the same conditions, the creep displacement decreases with the increase of GSI, and both the creep displacement and the effect degree enhance with the enlarged dilation angle. *e proposed solutions can be used in the stability analysis of surrounding rock in the construction and operation of salt cavern storage.


Introduction
A study on the analytical solution to stress and deformation of surrounding rock of simplified mechanical model is carried out based on the stress and strain redistribution.
e redistribution is caused by the salt cavern storage formation in underground caverns, which is formed in the salt cavern building with water solution in the deep salt formations.And the rock mass strength at the cavity wall is reduced during salt cavern leaching.In order to obtain reasonable deformation of surrounding rock, many scholars have carried out the analysis of stress and strain of surrounding rock based on the simplified mathematical model of underground caverns.Park and Kim [1] proposed a closed-form solution for the displacement of the plastic zone of a circular tunnel, by using the nonassociated flow rules combined with three definition methods for elastic strains of plastic zone.Sun et al. [2] proposed an approximate analytical solution for the circular tunnel with joints and analyzed the influence relationship of the rock mass parameters based on the H-B strength criterion and the assumption of the stress field under axisymmetrical and asymmetric stress conditions.Jimenez et al. [3] proposed a new equivalent method for symmetric circular tunnel based on the H-B criterion by assuming that the stress field in the plastic zone is linear.Fraldi and Guarracino [4] predicted the possibility of collapse accidents in circular tunnel excavation and proposed a simple and accurate solution based on the H-B, and they proposed elastoplastic analytical solutions of ground rock of tunnel based on convergence-confinement method.Rojat et al. [5] established a new dimensionless form of the generalized H-B failure criterion, completed the convergence calculations for tunnel in-depth by considering the problem of edge effects, and proposed new exact formulations for the case of an associated H-B plastic potential.Wang et al. [6] determined the H-B criterion parameters of the anhydrite intact rock and rock mass and performed a numerical simulation in FLAC3D.In order to analyze the stability of the deep cavern group, they proposed the element safety factor method based on the generalized H-B and stress distribution.Lu et al. [7] studied the plane strain elastic-plastic problem for the circular tunnel in H-B media subjected to nonhydrostatic stress, proposed the analytical optimization method for nonlinear equations.A solution of the elastoplastic interface is presented for the plane strain problem of a pressurized circular tunnel in H-B media subjected to nonhydrostatic stress at infinity.Llamas et al. [8] proposed a mini-CAES concept where the cavity is shallower than current CAES, applying the H-B criterion for the calculation of the subsidence profile.Other scholars have conducted a lot of research on the plane problem of circular tunnel [9][10][11][12].Hou et al. [13] proposed semianalytical and seminumerical solutions of the axisymmetric circular tunnel under ideal elastoplastic conditions by adopting the Levy-Mises constitutive relation and generalized H-B failure criterion.Cai et al. [14] analyzed the elastoplastic solution of the circular chamber under uniform stress field by the dimensionless H-B criterion and proposed that displacement can be effectively controlled by improving rock mass quality for soft rock chamber.Zhang et al. [15], by taking comprehensive influences of intermediate principal stress, brittle softening, dilation characteristic, Young's modulus, and elastic strains in the plastic zone into account, established a new solution for the displacement and characteristic curve of the plastic zone of the surrounding rock.Moreover, Liu et al. [16] analyzed the salt rock deformation under triaxial compressive stress condition, indicated the large compressive deformation properties of salt rock under higher confining pressure, and analyzed the stress-strain characteristics of the salt rock by using the logarithmic strain.Using large deformation (logarithmic strain), the problem of cavity expansion has been analyzed [17][18][19][20] in consideration of the shear dilation and the strain softening.e long-term behavior of salt causes the convergence of the salt cavity, so it is necessary to consider the stability and shrinkage of the cavity of long-term operation.Zhang and Xiong [21] proposed an analytical solution for the displacement of elasto-viscoplastic ground around the circular tunnel, which can not only reflect the effects of dilatancy and plastic softening on displacements but also take account into the creep property of the ground.Sun [22] discussed the research on rock rheology and the engineering applications.Yuan et al. [23] carried out viscoelastic-plasticity analysis of the rheological problem of deep soft rock tunnel taking the Drucker-Prager strength criterion as the plastic yielding conditions, and expanding effect of plastic zone is considered.Hou et al. [24] proposed an iterative algorithm based on the creep solution of constant passive support force for axisymmetric round well and used the iterative algorithm to calculate the creep deformation of surrounding rock.Wen et al. [25] and Cao et al. [26] determined an analytical solution for the rheological deformation of the surrounding rock by the generalized H-B criterion combined with the Nishihara model and considered the influence of the volumetric dilatancy on the solution.
e theoretical solutions for the deep spherical salt cavern are related to the complicated mathematical deduction; however, there is only a few theoretical research studies on the deep spherical cavern, especially the analytical solution of the rheological and the large deformation of the surrounding rock in the salt cavern storage project, and M-C and Drucker-Prager criterion has been used.Because the H-B criterion introduced the influence on rock mass underground construction, which can reflect the disturbance effect on salt cavern leaching process, it is more suitable for the practical engineering.From this point of view, this paper studied the elastoplastic (including two forms, i.e., logarithmic strain and engineering strain) and viscoelasticplastic displacement of surrounding rock in the salt cavern storage by introducing the H-B criterion combined with the equilibrium equation.is study is expected to provide theoretical reference for the design and construction of underground salt cavern.

Simplified Model for Salt Spherical Cavern Storage
In this study, the salt cavern storage is simplified into a deep spherical cavern.e distribution of the far field stress tends to be a approximate homogeneous hydrostatic pressure field where the stress value is p 0 .e rock mass of the surrounding rock is homogeneous and isotropic.e radius of the salt cavern is designed as R 0 , and the running pressure p i homogeneously acts on the cavern wall.Gravity and side pressure coefficient are not considered in this study.e coordinate origin of the salt cavern mechanics model is set in the center of the spherical cavern to meet the spherical symmetry.Initial in situ stress and running pressure are considered in this study; the salt rock conforms to the hypothesis of the ideal elastoplastic model, and the compressive stresses are considered positive.During the cavern building and operating process, there exist the elastic and plastic zones in the surrounding rock of the salt cavity.e region away from the salt cavern (R p ⟶ ∞) is in an elastic state, while the near-zone of the cavern (R 0 ⟶ R p ) enters into a plastic state.R p is designed as plastic radius.Figure 1 is the elastoplastic analysis model for surrounding rock of spherical salt cavern adits.And the surrounding rock satisfies continuity conditions in the elastic-plastic interface.

Strength Criterion
e generalized H-B criterion which is widely used can be written as 2 Advances in Civil Engineering where σ ci is the uniaxial compressive strength of the intact rock (MPa) and m b , s, and a are the petrographic constants of the rock mass.e disturbance coefficient D of the rock mass is introduced into the H-B criterion and a new method for parameter (m b , s, and a) value selection is proposed [27].
e parameters are based on the geological strength index (GSI).m b , s, and a are given as where GSI is the geological strength index, which varies from 10 to 100, it will be taken as 10 for a poor rock mass.D is the disturbance coefficient of the rock mass, which varies from 0 for an undisturbed in situ rock mass to 1 for a very disturbed rock mass.m i is the petrographic parameter reflecting the degree of stiffness of rock mass, which varies from 0 to 25 for the rock mass.s and a are petrographic parameters of the rock mass, which varies from 0 to 1, and the parameter a � 1/ 2 for most intact rock.

Elastoplastic Analytical Solutions
In the elastic zone, where R p ⟶ ∞, the boundary condition can be written as σ| r�R p � σ R p .Based on Lame's formula, compressive stresses are reckoned to be positive, and the stress in the elastic zone can be expressed as Based on the elasticity and without considering the effect of the initial in situ stress p 0 , strains can be written as In the elastic region, at elastic-plastic interface, r � R p , and the elastic strains can be expressed as where G � (E/(2(1 + μ))).e displacement caused by the salt cavern building with water solution can be obtained as E l a s t i c z o n e P l a s t i c z o n e By substituting the radial and tangential stresses in plastic zone with r � R p into equation ( 1), the nonlinear H-B strength criterion can be expressed as After combining equations ( 3) and ( 9) and considering the boundary condition σ p r � p i at r � R 0 , we can obtain the radial stress in the plastic zone: Combining equation ( 8) and considering the stress is continuous at the interface, r � R p , we can obtain that Combining equations ( 9) and ( 11) and considering the boundary condition Combining nonlinear equation ( 12) and using iteration method, the radial stress in the plastic zone can be determined.And, the tangential stress can be obtained based on equation (9).en, with equation ( 10), the plastic radius R p is ) As R p ≥ R 0 , the stress at the elastic-plastic interface should satisfy σ R p ≥ p i .And we inferred that the plastic radius is related to the in situ stress p 0 , the running pressure p i , and the H-B criterion parameters σ ci , m i , s, and a.

Deformation Analysis in Plastic
Zone.When the surrounding rock of salt cavern enters into plastic state, in order to determine the displacement in the plastic zone, it is necessary to establish displacement differential equations for surrounding rock using the geometric equation.erefore, it is necessary to determine the relationship between the strain components of the plastic zone where the total strain consists of elastic and plastic strain.
Many scholars have put forward different definitions for the elastic strain in the plastic zone.Some consider the elastic strain is very small, and it does not change.Hence, the elastic strain in the plastic zone is equal to the elastic strain at the elastic-plastic interface in the elastic zone.Others consider the stress-strain relationship conforms to the thick-walled cylinder theory and the generalized Hooke's law.e elastic strain in the plastic zone is a variable which changes with the position of the surrounding rock [15].In this study, the elastic strain in the plastic zone is discussed under the above two cases.
At the same time, the relationship between the plastic strain components is established using the flow rule.Based on the nonassociated flow law, plastic strain increments can be written as where h � (1 − sin ψ)/1 + sin ψ and ψ is the angle of dilation, which is obtained by triaxial compressive test.
After combining equations ( 14) and ( 15), we can obtain that Assuming the elastic strain in the plastic zone is a constant value, which is equal to the elastic strain at r � R p in the elastic zone and after combining equations ( 6) and ( 16), we can obtain that where Assuming the elastic strain in the plastic zone is a variable and abiding by Hooke's law, after combining equations ( 5) and ( 16), we can obtain that where To simplify the problem, we make the H-B parameter a � 0.5.Equation ( 10) can be written as where After combining equations ( 9) and ( 19), the tangential stress can be written as By substituting equations ( 19) and ( 20) into the function f 2 (r), we can obtain that 4

Advances in Civil Engineering
where In the infinitesimal deformation case, relationships of strains and displacement are as follows: By substituting equation ( 22) into equation ( 17), differential equation for the radial displacement u(r) can be written as h(du/dr) + 2(u/r) � f 1 (r).
e total displacement in the plastic zone can be obtained as en, considering continuous condition at the elasticplastic interface, we can obtain that 23) is caused by building, and it is obtained without considering the variation of the elastic strain in plastic zone.
In the same way, by substituting equation ( 22) into equation (18), differential equation for the radial displacement u(r) can be written as h(du/dr) + 2(u/r) � f 2 (r).e total displacement in the plastic zone is en, considering continuous condition at the elasticplastic interface, we can obtain that e displacement u p v (r) in equation (24), which is caused by building, is obtained with assuming the elastic strain in the plastic zone abiding by Hooke's law.
Table 1 shows displacements, which are obtained based on perfect elastoplastic constructive model, using H-B criterion and considering different definitions for elastic strain in the plastic zone.
In the case of ψ � 0 °(no plastic volume change), the displacements in the plastic zone are discussed.
(a) Assuming that the elastic strain in the plastic zone is a constant value, then f 1 (r) � 0, C 1 � (1/2)δR 3 p ; the displacement in the plastic zone is which is the same with that in the elastic zone, i.e., equation (7).(b) Assuming the elastic strain in the plastic zone is a variable, which abides by Hooke's law, we can obtain that e displacement in the plastic zone is Advances in Civil Engineering

Large Strain Solutions
According to the finite deformation theory, it is considered that the strain in the plastic zone is distributed in a logarithmic form.e logarithmic strain (true strain) can reflect the sum of the relative strain increments at various stages of total deformation, and it is closer to the actual accumulation process of deformation and the degree of large deformation.is study introduced the logarithmic strain ε �  r r 0 dl/l � ln(r/r 0 ) to represent large deformation of salt cavern [16,17].
e logarithmic strain indicates that the initial position r 0 gradually changes to r which will undergo several intermediate states during the deformation process of the plastic zone; the variables r and r 0 cannot be treated as constants during the calculation.en, the displacement of the salt cavern in plastic zone is derived adopting the large strain.
For large strain problem, the geometric equation can be written as Combing equations ( 17) and ( 29), we can obtain that where and u R p denotes the displacement in the elastic zone at elastic-plastic interface.e integration of equation (30) leads to the following expression for the displacement.
For the case of large deformation in plastic zone, which is caused by cavern building, equation (31) can fully express the displacement, which is obtained without consideration of the variation of elastic strain in the plastic zone.
In the same way, combing equations ( 18) and (39), we can obtain that e elastic strain in the plastic zone is a constant value, which equals to the elastic strain at the elasticplastic interface r � R 0 in the elastic zone i � 2 e elastic strain in the plastic zone is a variable, which abides by Hooke's law Advances in Civil Engineering e integration of equation (32) leads to the following form for displacement For the case of large deformation in plastic zone, which is caused by cavern building, equation (33) can fully express the displacement, which is obtained assuming the elastic strain and the total stress in the plastic zone abiding by Hooke's law.By the method of element changing in integral, equation (33) can be written as en, the numerical solution of equation ( 34) can be obtained by Romberg integration.
When ψ � 0 °(no plastic volume change), assuming the elastic strain in the plastic zone is a constant value, the displacement in the plastic zone can be obtained: Assuming the elastic strain in the plastic zone is a variable, the displacement in the plastic zone can be written as (36)

Viscoelastic-Plastic Deformation Analysis
e period of salt cavern building is long, and the working state is periodically injecting or producing circulated process.Creep displacement will take place in the surrounding rock in elastic and plastic zones, when the salt cavern kept operation at a reasonable inner pressure.It is necessary to investigate creep displacements, without considering the creep failure.As shown in Figure 2, the classical Nishihara model, which can fully reflect elasticity, viscoelasticity, and viscoplasticity in the transient creep and steady creep process [20,21,28], can be used to analyze the creep deformation of the surrounding rock of salt cavern.
e 1D creep constitutive equations for the classical Nishihara model are where E 0 stands for the elastic modulus, E 1 denotes the viscoelastic modulus, and η 1 and η 2 are the viscosity coefficients.e 1D creep constitutive equations show that the creep strain is related to time and the differential stress.

Viscoelastic Deformation Analysis.
e Nishihara model degenerates to a three-parameter model when the stress of the surrounding rock of salt cavern is lower than the yield stress σ s ; operators of the Nishihara model are listed in Table 2.

Advances in Civil Engineering
Taking Laplace inverse transform on equation (37), the creep strain in viscoelastic zone can be written as e analytic equation of displacement, i.e., equation ( 7) can be directly converted to Laplace space form as Substituting the Nishihara model's operators, which are presented in Table 2, into equation (39) and then taking Laplace inverse transform on equation (39), we can obtain the viscoelastic displacement in the viscoelastic zone R p ⟶ ∞, which is caused by caverns building, as the following expression: where G 0 � E/2(1 + μ), G 1 is viscoelastic shear modulus which can be determined experimentally, and η 1 is viscoelastic shear coefficient which can be determined experimentally.

Viscoplastic Deformation Analysis.
It is complicated to solve the displacement of surrounding rock of salt cavern in plastic zone R 0 ≤ r ≤ R p .e total strain of the surrounding rock is ε � ε e + ε p + ε c , which includes elastic strain ε e , plastic strain ε p , and creep strain ε c .In order to determine the total strain of plastic zone, assuming neither the elastic strain nor the creep strain causes the volume change, we can obtain that ε e r + 2ε e θ � 0, ε c r + 2ε c θ � 0.

􏼨 (41)
Based on equation (15), the strain increments relationship can be written as where I 0 � ε e θ + ε c θ .In order to simplify the problem, assuming the elastic strain in the viscoplastic zone is continuous at the viscoelastic-viscoplastic interface, based on equation ( 6), the tangential elastic strain in the viscoplastic zone ε e θ can be obtained.Hence, the elastic strain ε e θ in the viscoplastic zone is a constant independent of both the time t and the distance away from the center.e creep strain in the viscoplastic zone is obtained using viscoelasticity analysis method.When σ ≥ σ s , the 3D creep strain in the plastic zone can be written as e yield stress is defined as the difference between the tangential stress and the radial stress at elastic-plastic interface.Combining equations ( 19) and ( 20), we can obtain that σ s � A 1 + A 2 ln(R p /R 0 ).Combining equations ( 4) and ( 38) and combining equations ( 19), (20), and ( 43 Assuming that the creep strain in the viscoplastic zone, which can be approximated by the average of the creep strain in the viscoplastic zone at the elastoplastic interface and the creep strain at the cavern wall calculated based on the differential stress, is an r-independent value [23,26] and substituting ε e θ and equations (36) into (34), we can obtain that e parameter I 0 depends only on the time t and does not depend on the position of the surrounding rock r.Referring to the solving approach in the above chapter and combing equations ( 22) and ( 42), the first-order nonhomogeneous constant coefficient differential equation for u(r) can be written as h(du/dr) + 2(u/r) � 2(1 − h)I 0 .e viscoplastic displacement in the plastic zone can be obtained as where p .When t � 0, equation ( 45) is the same as equation (23), which can describe the displacement in the plastic zone without consideration of the creep strain.
When ψ � 0 °(no plastic volume change), the equation ( 45) is the same as equation ( 40), which can describe the displacement in the viscoelastic zone, and we can obtain that (46)

Case Study and Discussion
Taking the salt cavern gas storage being excavated in China as an example, the structural parameters and mechanical parameters are given in Table 3. Table 4 Shows the mechanical parameters of salt rock which are involved in the model.
e literature investigation results show that GSI and m i of the H-B criterion have a similar significant influence on the stress and displacement of the surrounding rock.By taking the different values 60 and 80 of GSI into account, a study has been conducted to evaluate their influences on the displacement.Table 5 shows the other parameters for calculation.In order to validate the result in this paper, as shown in Figure 3, the displacements, which are obtained adopting two different definitions of elastic strain of the plastic zone, and stresses are contrasted; it is not surprising that the results agree with previous outcomes [29].
In order to validate the analytical solutions, numerical computations are carried out for the salt cavity storage.Under the same conditions, the displacement of the surrounding rock can be derived based on the H-B and M-C criterions with the aid of the FLAC3D.In order to establish the model of the salt cavern under uniform stress field, set stress boundary conditions in FLAC3D 5.0.As shown in Figure 4, principal stress and displacement, which are obtained under the same conditions, are extracted from FLAC3D by setting GSI � 60 and using the H-B criterion in FLAC3D.Figure 5 shows the comparison between numerical solutions and analytical solutions.In Figure 5, the results show that the distribution of the principal stress and the displacement is fundamentally identical with the analytical solution.e results indicate that displacements in the plastic zone gradually decrease with the distance away from the salt cavity; the difference of the displacements that obtained under different conditions gradually declined.But the difference between displacements on the plastic zone still remains under different conditions.Compared with numerical solutions at the cavity wall, the analytical one is smaller.No matter the analytical solutions or the numerical solutions, the displacement for M-C rock mass lies in the range of results when GSI � 80 to GSI � 60. e results close to that based on the M-C criterion can be obtained by appropriate readjustment of the H-B criterion parameter GSI, but the H-B criterion, which introduced the structural characteristics of the rock mass, has better applicability than the M-C criterion.

Effect of Large Strain and Small Strain: eory to Results.
e displacement in the plastic zone is plotted in Figure 6 for cases of H-B parameter GSI � 60 and GSI � 80.It can be seen from Figure 6 that the range of the plastic zone reduced (the plastic radius decreased from 44.79 m to 35.91 m) and the displacement has decreased, with the H-B parameter GSI increasing from 60 to 80. e difference of the displacement, which caused no matter by two different strain theories or by two different elastic strain definitions, is greater, when GSI � 60.When the elastic strain of plastic zone takes on stationary values, the displacement obtained with small strain theory is larger than large strain theory.But when the elastic strain abides by Hooke's law, the displacement obtained with small strain theory is smaller than large strain theory.In large strain stage, for cases of GSI � 60 and GSI � 80, the difference of displacement at the cavity wall that is caused by two definitions for elastic strains, is 0.030 m and 0.008 m, respectively.However, in small strain stage, for cases of GSI � 60 and GSI � 80, the difference of displacement is 0.008 m and 0.003 m, respectively.It indicates that the difference of displacements has a little difference for good rock mass; while it has great difference for poor rock mass, such as broken rock mass, fractured rock mass, or rock formation with natural fractures.In large strain stage, for two cases of different assumptions about the elastic strains, the difference of displacement at the cavity wall, which is caused by two different values of GSI, is 0.035 m and 0.057 m, respectively.In small strain stage, the difference of displacements for different assumptions about the elastic strains is 0.044 m and 0.039 m, respectively.It indicates that the influence of different hypothesis of the elastic strains on the displacement is more obvious, as considering large strain.
Figure 7 shows the displacements, which are obtained by taking two different elastic strain definitions in the plastic zone and two different strain theories into account, in the plastic zone for cases of ψ � 0, ψ � 0.25φ, and ψ � 0.5φ.It can be seen from Figure 7 that the displacement at the same point increases with the enlarged dilation angle.e displacements intersect with each other at one point (displacement at elastoplastic interface) with the distance away from the salt cavity, which indicate that the influence of dilation angle on the displacements gradually reduces.When ψ � 0 (no plastic volume change), for cases of different elastic strain definitions, the displacements with reference to small strain are almost the same.It is important to note that when the displacement is obtained Advances in Civil Engineering based on large strain theory, the various definitions of elastic strain in the plastic zone affect the displacements differently.Note that when the elastic strain takes a constant value, the displacements with reference to large deformation theory agree with the displacements with reference to small deformation theory; however, when the elastic strain conforms to Hooke's law, the displacement which is obtained based on large deformation theory is greater than the displacement obtained based on small deformation theory.When ψ ≠ 0 (plastic volume increase), whether ψ � 0.25φ or ψ � 0.5φ, the displacements obtained based on large strain are smaller than that one obtained based on small strain when the elastic strain of plastic zone is constant.However, the displacements obtained based on large strain are greater than those obtained based on small strain when the elastic strain abides by Hooke's law.e above results relate to the solutions obtained by using the semianalytical method of Romberg numerical integration.Figure 8 shows relationships between the dilation angles and the relative increment of the displacements at the cavity wall.It can be seen from Figure 8 that the relative increment increases with the dilation angle.e dilation angle has the most significant influence on the displacement with reference to small strain, while it has the least influence on the displacement with reference to large strain when the elastic strain of plastic zone is constant.
Figure 9 illustrates relationships between the operation pressure and the displacement at the cavity wall.It can be seen from Figure 9, the displacement decreases with the increase of the operating pressure, and the maximum difference of the displacement is at operation pressure 2 MPa. e displacements at the cavity wall for different cases are 0.36 m, 0.39 m, 0.42 m, 0.55 m, and 0.74 m, respectively.ere is difference between displacements obtained based on different strain theories.When the elastic strain in the plastic zone is constant, the displacement with reference to large strain theory is smaller than the one with reference to small strain theory.However, there is a contrary relationship between the displacements, which is closer to the numerical results, when the elastic strain conforms to Hooke's law.erefore, appropriate control of inner pressure can reduce the displacements at the cavity wall effectively, and it can also reduce the impact of the two different elastic strain definitions and the two different strain theories.6.
Take the point r � 30 m and r � 35 m as the observation points; it can be seen from Figure 10 that the rheological displacement equation derived in this result is correct.e displacements of the surrounding rock decrease with the increase of GSI; the higher the parameter GSI, the lower the displacement.e rheological displacement of the surrounding  Advances in Civil Engineering rock can be well controlled and can enter the second creep stage (steady-state creep stage) faster.e displacements at the cavern wall still increase linearly in the second creep stage, and the increasing rate is related to the viscosity coefficient μ 2 .
Figure 11 shows the relationship between time and creep displacements at the cavern wall for cases of ψ � 0, ψ � 0.25φ, and ψ � 0.5φ.As can be seen from Figure 11, the initial value of the creep displacement gradually decreases with the reduced angle of dilation at the cavern wall.e total creep displacement increased with the enlarged dilation angle within 200h.e dilation angle increases from 0 to 0.5φ, the creep displacements are 0.1 m, 0.13 m, and 0.16 m, respectively, and the relative values of displacement (u(t � 200) − u(t � 0)/u(t � 0)) × 100% are 1.3%, 9.9%, and 16.2%, respectively.Both the variation gradient of the creep displacement rate and the creep rate in the steady-state creep stage increase with the enlarged dilation angle.Also, the region of the primary creep stage becomes larger with the enlarged dilation angle.Hence, the angle of dilation has great influence on the creep displacement.Advances in Civil Engineering equations, although the constitutive relation and the plastic flow rule are consistent in solving process.e results of the spherical model are more applicable to the natural conditions of the salt cavern storage.(2) e displacements of plastic zone are analyzed by taking two different elastic strain definitions of plastic zone, large strain, and small strain into account.e influence of H-B parameter GSI, operating pressure, and dilation angle on the four analytical solutions is studied.e maximum difference of the displacements in the plastic zone is at the cavern wall.e difference decreased gradually with the distance away from the cavern wall.e displacements of the surrounding rock can be controlled effectively with the increase of the H-B parameter GSI or the operating pressure, while it increases with the enlarged dilation angle.When the elastic strain of plastic zone is constant, the displacements obtained based on large strain are smaller than that obtained based on small strain under the same conditions, while there is a contrary relationship when the elastic strain conforms to Hooke's law.e influence of different elastic strain definitions on the displacements that obtained based on large strain is more significant.When the elastic strain of plastic zone is constant, the influence of dilation angle on the displacements is most obvious,

Figure 1 :
Figure 1: e ideal elastoplastic model of spherical salt cavern.(a) e elastic-plastic model for spherical salt cavern.(b) Material behavior model of salt rock.
), creep strains ε c θp     r�R p and ε c θp     r�R 0 can be obtained, respectively.
e creep curves of surrounding rock of the salt cavern have been obtained at r � 30 m and r � 35 m, by taking different values of GSI and taking the creep model parameters of surrounding rock of the salt cavern as shown in Table

Figure 7 :( 1 )
Figure 7: Relationships between plastic zone displacements and dilation angle (different elastic strain definitions and different strain theories).

Figure 8 :Figure 11 :
Figure 8: Relationships between displacements at the cavity wall and dilation angle (different elastic strain definitions and different strain theories).

Table 1 :
Summary of the displacements in the plastic zone.Geometric equations ε r � (du/dr) ε θ � (u/r) Displacement in the plastic zone in a compact form r

Table 3 :
Parameters for salt cavern storage.

Table 4 :
Parameters of salt rock.