A 2D stability analysis of the rock surrounding underground liquified natural gas storage cavern based on COMSOL Multiphysics

Underground liquified natural gas(LNG)storage is essential in guaranteeing national energy strategic reserves,and its construction is being accelerated.The stability of surrounding rock of underground LNG storage caverns under stress-low temperature coupling effect is the key factor determining the feasibility of LNG storage.First,a mathematical model used for controlling the stress-low temperature coupling and the processes of rock damage evolution is given,followed by a 2-D numerical execution process of the mathematical model mentioned above described based on Comsol Multiphysics and Matlab code.Finally,a series of 2-D simulations are performed to study the influence of LNG storage cavern layout,burial depth,temperature and internal pressure on the stability of surrounding rocks of these underground storage caverns.The results indicate that all the factors mentioned above affect the evolution of deformation and plastic zone of surrounding rocks.The research results contribute to the engineering design of underground LNG storage caverns.

Prediction model of surrounding rock deformation in doublecontinuous-arch tunnel based on the ABC-WNN

Purpose–The wavelet neural network(WNN)has the drawbacks of slow convergence speed and easy falling into local optima in data prediction.Although the artificial bee colony(ABC)algorithm has strong global optimization ability and fast convergence speed,it also has the drawbacks of slow speed while finding the optimal solution and weak optimization ability in the later stage.Design/methodology/approach–This article uses an ABC algorithm to optimize the WNN and establishes an ABC-WNN analysis model.Based on the example of the Jinan Yuhan underground tunnel project,the deformation of the surrounding rock of the double-arch tunnel crossing the fault fracture zone is predicted and analyzed,and the analysis results are compared with the actual detection amount.Findings–The comparison results show that the predicted values of the ABC-WNN model have a high degree of fitting with the actual engineering data,with a maximum relative error of only 4.73%.On this basis,the results show that the statistical features of ABC-WNN are the lowest,with the errors at 0.566 and 0.573,compared with the single back propagation(BP)neural network model and WNN model.Therefore,it can be derived that the ABC-WNN model has higher prediction accuracy,better computational stability and faster convergence speed for deformation.Originality/value–This article uses firstly the ABC-WNN for the deformation analysis of double-arch tunnels.This attempt laid the foundation for artificial intelligence prediction in deformation analysis of multiarch tunnels and small clearance tunnels.It can provide a new and effective way for deformation prediction in similar projects.

Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

Research on the stability analysis and design of soil tunnel surrounding rock

The paper first analyzes the failure mechanism and mode of tunnel according to model experiments and mechanical calculation and then discusses the deficiency of taking the limit value of displacement around the tunnel and the size of the plastic zone of surrounding rock as the criterion of stability. So the writers put forward the idea that the safety factor of surrounding rock calculated through strength reduction FEM(finit element method) should be regarded as the criterion of stability,which has strict mechanical basis and unified standard and would not be influenced by other factors. The paper also studies the safety factors of tunnel surrounding rock (safety factors of shear and tension failure) and lining and some methods of designing and calculating tunnels. At last,the writers take the loess tunnel for instance and show the design and calculation results of two-lane railway tunnel.

Numerical Simulation of Surrounding Rock Deformation and Grouting Reinforcement of Cross-Fault Tunnel under Different Excavation Methods

Tunnel construction is susceptible to accidents such as loosening, deformation, collapse, and water inrush, especiallyunder complex geological conditions like dense fault areas. These accidents can cause instability and damageto the tunnel. As a result, it is essential to conduct research on tunnel construction and grouting reinforcementtechnology in fault fracture zones to address these issues and ensure the safety of tunnel excavation projects. Thisstudy utilized the Xianglushan cross-fault tunnel to conduct a comprehensive analysis on the construction, support,and reinforcement of a tunnel crossing a fault fracture zone using the three-dimensional finite element numericalmethod. The study yielded the following research conclusions: The excavation conditions of the cross-fault tunnelarray were analyzed to determine the optimal construction method for excavation while controlling deformationand stress in the surrounding rock. The middle partition method (CD method) was found to be the most suitable.Additionally, the effects of advanced reinforcement grouting on the cross-fault fracture zone tunnel were studied,and the optimal combination of grouting reinforcement range (140°) and grouting thickness (1m) was determined.The stress and deformation data obtained fromon-site monitoring of the surrounding rock was slightly lower thanthe numerical simulation results. However, the change trend of both sets of data was found to be consistent. Theseresearch findings provide technical analysis and data support for the construction and design of cross-fault tunnels.

Heterogeneous information phase space reconstruction and stability prediction of filling body–surrounding rock combination

Traditional research believes that the filling body can effectively control stress concentration while ignoring the problems of unknown stability and the complex and changeable stress distribution of the filling body–surrounding rock combination under high-stress conditions.Current monitoring data processing methods cannot fully consider the complexity of monitoring objects,the diversity of monitoring methods,and the dynamics of monitoring data.To solve this problem,this paper proposes a phase space reconstruction and stability prediction method to process heterogeneous information of backfill–surrounding rock combinations.The three-dimensional monitoring system of a large-area filling body–surrounding rock combination in Longshou Mine was constructed by using drilling stress,multipoint displacement meter,and inclinometer.Varied information,such as the stress and displacement of the filling body–surrounding rock combination,was continuously obtained.Combined with the average mutual information method and the false nearest neighbor point method,the phase space of the heterogeneous information of the filling body–surrounding rock combination was then constructed.In this paper,the distance between the phase point and its nearest point was used as the index evaluation distance to evaluate the stability of the filling body–surrounding rock combination.The evaluated distances(ED)revealed a high sensitivity to the stability of the filling body–surrounding rock combination.The new method was then applied to calculate the time series of historically ED for 12 measuring points located at Longshou Mine.The moments of mutation in these time series were at least 3 months ahead of the roadway return dates.In the ED prediction experiments,the autoregressive integrated moving average model showed a higher prediction accuracy than the deep learning models(long short-term memory and Transformer).Furthermore,the root-mean-square error distribution of the prediction results peaked at 0.26,thus outperforming the no-prediction method in 70%of the cases.

Face stability analysis of circular tunnels in layered rock masses using the upper bound theorem

An analysis of tunnel face stability generally assumes a single homogeneous rock mass.However,most rock tunnel projects are excavated in stratified rock masses.This paper presents a two-dimensional(2D)analytical model for estimating the face stability of a rock tunnel in the presence of rock mass stratification.The model uses the kinematical limit analysis approach combined with the block calculation technique.A virtual support force is applied to the tunnel face,and then solved using an optimization method based on the upper limit theorem of limit analysis and the nonlinear Hoek-Brown yield criterion.Several design charts are provided to analyze the effects of rock layer thickness on tunnel face stability,tunnel diameter,the arrangement sequence of weak and strong rock layers,and the variation in rock layer parameters at different positions.The results indicate that the thickness of the rock layer,tunnel diameter,and arrangement sequence of weak and strong rock layers significantly affect the tunnel face stability.Variations in the parameters of the lower layer of the tunnel face have a greater effect on tunnel stability than those of the upper layer.

Deformation mechanism and collapse treatment of the rock surrounding a shallow tunnel based on on-site monitoring

When tunnels are constructed at shallow depths in areas with poor geological conditions,such as portal sections,valleys and hillsides in regions with granitic bedrock,considerable excavation-induced deformation of the surrounding rock may occur,potentially resulting in tunnel collapses.The main reason for these problems is the lack of understanding of the deformation mechanism and evolution of the soft granitic rock surrounding the tunnel and the adoption of inappropriate construction technology and methods.This article analyzes the deformation mechanism of the rock surrounding a shallow tunnel based on in situ monitoring data as a case study and suggests that certain measures should be taken to effectively control the deformation of the surrounding rock and to minimize the potential for tunnel collapse.The results show that the deformation of the granitic soil surrounding the tunnel can be divided into three stages:the rapid deformation stage,the slow deformation stage and the stabilization stage.Appropriate construction methods should be carefully selected to ensure safety during tunnel excavation in the first stage.To avoid secondary disasters caused by tunnel collapses,three treatment measures may be implemented as part of safety management:enhancing the monitoring of the surrounding rock deformation,adjusting the construction methods and optimizing the support systems.In particular,accurate monitoring data and timely information feedback play a vital role in tunnel construction.Therefore,engineers with considerable engineering experience and professional knowledge are needed to analyze the monitoring data and make accurate predictions of tunnel deformation to ensure that reasonable measures are taken in the process of shallow tunnel excavation.

Numerical simulation on deformation character of surrounding rock masses of Changjiashan tunnel through the gob of coalmine

Based on the construction project of the Changjiashan tunnel of the freeway,the variety rule of surrounding rock masses of the tunnel through the gob of coalmine wasstudied by using of finite element methed(FEM).The status of the stress and strain,thevariety of the plastic area were simulated in the whole rock mass before and after thetunnel was excavated.The characters of stress and deformation of surrounding rockmasses were analyzed when the tunnel was built.It concluded from the numerical simula-tion that the influence on the tunneling is great when the tunnel passing through the gob ofcoalmine is excavated,and the relative measures should be taken.

Engineering Geological Classification of Surrounding Rock Stability in Coal Mine Tunnels and Its Engineering Applications

The stability type and spatial distribution of surrounding rock are an important basis for the layout of hae tunuels and the selection of support patterns in a mine design. Based on a lot of investigations and testing studies,a new engineering geological classification scheme of surrounding rock stability is put forward,which is easy to be applied,and reliable verified by examples. According to the classification system,the spatial divisions of surrounding rock stability can be delineated in an exploration progamme, providing relevant engineering geological informations for mine designers to prevent them from making the uurealistic tunnel layout and support

Stability mechanism and control of the pumpable supports in longwall recovery room

The load-bearing performance(LBP)of pumpable supports(PPS)is crucial for the stability of longwall pre-driven recovery room(PRR)surrounding rock.However,the unbalanced bearing coefficient(UBC)of the PPS(undertaking unequal load along the mining direction)has not been investigated.A mechanical model of the PRR was established,considering the main roof cantilever beam structure,to derive an assessment formula for the load,the failure criteria,and the UBC of the PPS.Subsequently,the generation mechanisms,and influencing factors of the UBC were revealed.Global sensitivity analysis shows that the main roof hanging length(l_(2))and the spacing between the PPS(r)significantly impact the UBC.A novel design of the PPS and the coupling control technology were proposed and applied to reduce the UBC of the PPS in the adjacent longwall PRR.Monitor results showed no failure of the PPS at the test site,with the UBC(ζ)reduced to 1.1 consistent with the design value(1.15)basically,fully utilizing the collaborative LBP of the PPS.Finally,the maximum roof-to-floor convergence of the PRR was 234 mm,effectively controlling the stability of the surrounding rock of the PRR and ensuring the mining equipment recovery.

Study on creep deformation and energy development of underground surrounding rock under four‐dimensional support

There is an urgent need to develop optimal solutions for deformation control of deep high‐stress roadways,one of the critical problems in underground engineering.The previously proposed four‐dimensional support(hereinafter 4D support),as a new support technology,can set the roadway surrounding rock under three‐dimensional pressure in the new balanced structure,and prevent instability of surrounding rock in underground engineering.However,the influence of roadway depth and creep deformation on the surrounding rock supported by 4D support is still unknown.This study investigated the influence of roadway depth and creep deformation time on the instability of surrounding rock by analyzing the energy development.The elastic strain energy was analyzed using the program redeveloped in FLAC3D.The numerical simulation results indicate that the combined support mode of 4D roof supports and conventional side supports is highly applicable to the stability control of surrounding rock with a roadway depth exceeding 520 m.With the increase of roadway depth,4D support can effectively restrain the area and depth of plastic deformation in the surrounding rock.Further,4D support limits the accumulation range and rate of elastic strain energy as the creep deformation time increases.4D support can effectively reduce the plastic deformation of roadway surrounding rock and maintain the stability for a long deformation period of 6 months.As confirmed by in situ monitoring results,4D support is more effective for the long‐term stability control of surrounding rock than conventional support.

Stability classification model of mine-lane surrounding rock based on distance discriminant analysis method

Based on the principle of Mahalanobis distance discriminant analysis (DDA) theory, a stability classification model for mine-lane surrounding rock was established, including six indexes of discriminant factors that reflect the engineering quality of surrounding rock: lane depth below surface, span of lane, ratio of directly top layer thickness to coal thickness, uniaxial comprehensive strength of surrounding rock, development degree coefficient of surrounding rock joint and range of broken surrounding rock zone. A DDA model was obtained through training 15 practical measuring samples. The re-substitution method was introduced to verify the stability of DDA model and the ratio of mis-discrimination is zero. The DDA model was used to discriminate 3 new samples and the results are identical with actual rock kind. Compared with the artificial neural network method and support vector mechanic method, the results show that this model has high prediction accuracy and can be used in practical engineering.

Upper bound solution of surrounding rock pressure of shallow tunnel under nonlinear failure criterion

The nonlinear Baker failure criterion is introduced into the upper-bound limit analysis to examine the face stability of a shallow tunnel. The tunnel face under the ultimate limit state is analyzed from the perspective of energy balance. The work rates of external forces and internal energy dissipation are calculated. An analytical solution of necessary face pressures is derived. The optimal upper-bound solution of the face pressures is obtained by optimization. The results show that the three dimensionless parameters A, T, n of nonlinear Baker failure criterion have different effects on the necessary face pressures and the pattern failure mechanisms ahead of tunnel face. A is the most important one;n takes the second place, and T is the least one. The computed necessary face pressures are nonlinearly increasing when A is reduced. Combined with the actual monitoring data of Taxia tunnel, the calculation results in this paper is verified. It is suggested that the tunnel face supports should be strengthened timely in soft rocks to prevent the occurrence of face collapse.

Propagation characteristics of vibration waves induced in surrounding rock by tunneling blasting

The effect of blasting vibration waves on surrounding rock and supporting structures is an important field in underground engineering. In this paper, the separation variable method is used to solve the displacement potential function for the propagation of the blasting vibration waves. In the axis coordinate system, the particle motion and stress change with axial distance, radial distance and time is obtained in surrounding rock. The peak particle velocity law in surrounding rock under different blast loads and surrounding rock parameters is discussed.In addition, the particle vibration characteristics in the surrounding rock are studied using numerical simulations method. The results shows that the peak particle velocity in surrounding rock appears negative exponent attenuation with the increase of axial distance, but it appears positive and negative fluctuations in radial direction. This phenomenon is a new discovery and it has been rarely investigated before. Moreover, the peak particle velocity attenuates more quickly and intensely in the near blasting field,which means that the supporting structure in a shorter distance away from the heading face is vulnerable to the impact of blasting vibration. Theattenuation of blasting vibration velocity is closely related to charge length, blasting load amplitude,attenuation index and rock elastic modulus. The numerical simulation accomplishes the same results and then demonstrates the validity of theoretical results.

Influence of underground water seepage flow on surrounding rock deformation of multi-arch tunnel

Based on a typical multi-arch tunnel in a freeway, the fast Lagrangian analysis of continua in 3 dimensions(FLAC3D) was used to calculate the surrounding rock deformation of the tunnel under which the effect of underground water seepage flow was taken into account or not. The distribution of displacement field around the multi-arch tunnel, which is influenced by the seepage field, was gained. The result indicates that the settlement values of the vault derived from coupling analysis are bigger when considering the seepage flow effect than that not considering. Through the contrast of arch subsidence quantities calculated by two kinds of computation situations, and the comparison between the calculated and measured value of tunnel vault settlement, it is found that the calculated value(5.7-6.0 mm) derived from considering the seepage effect is more close to the measured value(5.8-6.8 mm). Therefore, it is quite necessary to consider the seepage flow effect of the underground water in aquiferous stratum for multi-arch tunnel design.

Study on calculation of rock pressure for ultra-shallow tunnel in poor surrounding rock and its tunneling procedure

A computational method of rock pressure applied to an ultra-shallow tunnel is presented by key block theory, and its mathematical formula is proposed according to a mechanical tunnel model with super-shallow depth. Theoretical analysis shows that the tunnel is subject to asymmetric rock pressure due to oblique topography. The rock pressure applied to the tunnel crown and sidewall is closely related to the surrounding rock bulk density, tunnel size, depth and angle of oblique ground slope. The rock pressure applied to the tunnel crown is much greater than that to the sidewalls, and the load applied to the left side-wall is also greater than that to the right sidewall. Mean-while, the safety of the lining for an ultra-shallow tunnel in strata with inclined surface is affected by rock pressure and tunnel support parameters. Steel pipe grouting from ground surface is used to consolidate the unfavorable surrounding rock before tunnel excavation, and the reinforcing scope is proposed according to the analysis of the asymmetric load induced by tunnel excavation in weak rock with inclined ground surface. The tunneling procedure of bench cut method with pipe roof protection is still discussed and carried out in this paper according to the special geological condition. The method and tunneling procedure have been successfully utilized to design and drive a real expressway tunnel. The practice in building the super-shallow tunnel＆amp;nbsp;has proved the feasibility of the calculation method and tunneling procedure presented in this paper.

Development of a new type of foam concrete and its application on stability analysis of large-span soft rock tunnel

The long-term stability of large-span soft rock tunnel is influenced greatly by the creep effect of surrounding rock.The development of a new type of foam concrete which has the property of high compressibility and low ductility was introduced.And it was made as filling material of reserved deformation layer between the first lining and the second lining used in large-span soft rock tunnel.The effect of the new type of foam concrete was simulated as filling material of reserved deformation layer using numerical simulation.Through the comparison with the common large-span soft rock tunnel,the vault settlement and surrounding convergence are reduced by about 61% and 45%,respectively,after creep of 100 a.And in the second lining,the plastic zone reduces apparently and the maximum equivalent plastic strain decreases relatively.So,it can be found that the application of the new type of foam concrete as the filling material of reserved deformation layer can relieve the excessive force in second lining induced by rock creep,reduce its deformation and improve the stability of tunnel.

Automatic identification of rock discontinuity and stability analysis of tunnel rock blocks using terrestrial laser scanning

Local geometric information and discontinuity features are key aspects of the analysis of the evolution and failure mechanisms of unstable rock blocks in rock tunnels.This study demonstrates the integration of terrestrial laser scanning(TLS)with distinct element method for rock mass characterization and stability analysis in tunnels.TLS records detailed geometric information of the surrounding rock mass by scanning and collecting the positions of millions of rock surface points without contact.By conducting a fuzzy K-means method,a discontinuity automatic identification algorithm was developed,and a method for obtaining the geometric parameters of discontinuities was proposed.This method permits the user to visually identify each discontinuity and acquire its spatial distribution features(e.g.occurrences,spac-ings,trace lengths)in great detail.Compared with hand mapping in conventional geotechnical surveys,the geometric information of discontinuities obtained by this approach is more accurate and the iden-tification is more efficient.Then,a discrete fracture network with the same statistical characteristics as the actual discontinuities was generated with the distinct element method,and a representative nu-merical model of the jointed surrounding rock mass was established.By means of numerical simulation,potential unstable rock blocks were assessed,and failure mechanisms were analyzed.This method was applied to detection and assessment of unstable rock blocks in the spillway and sand flushing tunnel of the Hongshiyan hydropower project after a collapse.The results show that the noncontact detection of blocks was more labor-saving with lower safety risks compared with manual surveys,and the stability assessment was more reliable since the numerical model built by this method was more consistent with the distribution characteristics of actual joints.This study can provide a reference for geological survey and unstable rock block hazard mitigation in tunnels subjected to complex geology and active rockfalls.

Forward prediction for tunnel geology and classification of surrounding rock based on seismic wave velocity layered tomography

Excavation under complex geological conditions requires effective and accurate geological forward-prospecting to detect the unfavorable geological structure and estimate the classification of surround-ing rock in front of the tunnel face.In this work,a forward-prediction method for tunnel geology and classification of surrounding rock is developed based on seismic wave velocity layered tomography.In particular,for the problem of strong multi-solution of wave velocity inversion caused by few ray paths in the narrow space of the tunnel,a layered inversion based on regularization is proposed.By reducing the inversion area of each iteration step and applying straight-line interface assumption,the convergence and accuracy of wave velocity inversion are effectively improved.Furthermore,a surrounding rock classification network based on autoencoder is constructed.The mapping relationship between wave velocity and classification of surrounding rock is established with density,Poisson’s ratio and elastic modulus as links.Two numerical examples with geological conditions similar to that in the field tunnel and a field case study in an urban subway tunnel verify the potential of the proposed method for practical application.