Updating the models and uncertainty of mechanical parameters for rock tunnels using Bayesian inference

Rock mechanical parameters and their uncertainties are critical to rock stability analysis,engineering design,and safe construction in rock mechanics and engineering.The back analysis is widely adopted in rock engineering to determine the mechanical parameters of the surrounding rock mass,but this does not consider the uncertainty.This problem is addressed here by the proposed approach by developing a system of Bayesian inferences for updating mechanical parameters and their statistical properties using monitored field data,then integrating the monitored data,prior knowledge of geotechnical parameters,and a mechanical model of a rock tunnel using Markov chain Monte Carlo(MCMC)simulation.The proposed approach is illustrated by a circular tunnel with an analytical solution,which was then applied to an experimental tunnel in Goupitan Hydropower Station,China.The mechanical properties and strength parameters of the surrounding rock mass were modeled as random variables.The displacement was predicted with the aid of the parameters updated by Bayesian inferences and agreed closely with monitored displacements.It indicates that Bayesian inferences combined the monitored data into the tunnel model to update its parameters dynamically.Further study indicated that the performance of Bayesian inferences is improved greatly by regularly supplementing field monitoring data.Bayesian inference is a significant and new approach for determining the mechanical parameters of the surrounding rock mass in a tunnel model and contributes to safe construction in rock engineering.

Study on the new shell bolting and shotcrete support and its application in soft rock tunnels

A new kind of tunnel support was put forward on the basis of the anchor spraysupport principle. The mechanics of the new three-dimensional steel bar shotcrete liningsupport was studied and the structure's internal forces were analyzed. The model experiment was done relying on the industrial test. The conclusion of numerical calculationsproved that the ANSYS program is reasonable and creditable. It was compared to otherkinds of support that are commonly used in soft rock tunnels. The technique and economiccontrasts of the typical tunnel with support three-dimensional steel bar were completed.

The role of polyurethane foam compressible layer in the mechanical behaviour of multi-layer yielding supports for deep soft rock tunnels

The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not been thoroughly investigated.To fill this gap,large-scale model tests were conducted in this study.The synergistic load-bearing mechanics were analyzed using the convergenceconfinement method.Two types of multi-layer yielding supports with different thicknesses(2.5 cm,3.75 cm and 5 cm)of PU compressible layers were investigated respectively.Digital image correlation(DIC)analysis and acoustic emission(AE)techniques were used for detecting the deformation fields and damage evolution of the multi-layer yielding supports in real-time.Results indicated that the loaddisplacement relationship of the multi-layer yielding supports could be divided into the crack initiation,crack propagation,strain-hardening,and failure stages.Compared with those of the stiff support,the toughness,deformability and ultimate load of the yielding supports were increased by an average of 225%,61%and 32%,respectively.Additionally,the PU compressible layer is positioned between two primary linings to allow the yielding support to have greater mechanical properties.The analysis of the synergistic bearing effect suggested that the thickness of PU compressible layer and its location significantly affect the mechanical properties of the yielding supports.The use of yielding supports with a compressible layer positioned between the primary and secondary linings is recommended to mitigate the effects of high geo-stress in squeezing tunnels.

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.

Time-dependent squeezing deformation mechanism of tunnels in layered soft-rock stratum under high geo-stress

Large squeezing deformation of layered soft rock tunnel under high geo-stress has a significant time-dependent deformation behavior.In this paper,we studied the deformation mechanism during the construction period of deep-buried softrock tunnel by means of a combination of field observations and a numerical method.First,a new classification criterion for large deformations based on the power exponent variation law between the deformation and the strength-stress ratio is proposed.Then,the initial damage tensor reflecting the bedding plane(joint)distribution and an equivalent damage evolution equation derived from the viscoplastic strain are introduced based on the geometric research method,i.e.,a new rheological damage model(RDL model)of layered soft rock is established consisting of elastic,viscous,viscoelastic,viscoplastic and plastic elements.A field test was conducted on the Maoxian tunnel in Sichuan province,southwestern China,which is in broken phyllite(layered soft rock)under high geo-stress.The tunnel has experienced large deformation due to serious squeezing pressure,thus we adopted double primary support method to overcome the supporting structure failure problems.The rheological parameters of phyllite in the Maoxian tunnel were recognized by using SA-PSO optimization,and the RDL model does a good job in describing the time-dependent deformation behavior of a layered soft-rock tunnel under high geo-stress.Thus,the RDL model was used to investigate the supporting effect and bearing mechanism of the double primary support method.Compared with the single primary support method,the surrounding rock pressure,secondary lining force,surrounding rock deformation,and the depth of the damage to the rock mass was reduced by 40%-60%after the double primary support method was used.

Zonal disintegration phenomenon in enclosing rock mass surrounding deep tunnels——Elasto-plastic analysis of stress field of enclosing rock mass

The zonal disintegration phenomenon （ZDP） is a typical phenomenon in deep block rock masses. In order to investigate the mechanism of ZDP, an improved non-linear Hock-Brown strength criterion and a bi-linear constitutive model of rock mass were used to analyze the elasto-plastic stress field of the enclosing rock mass around a deep round tunnel. The radius of the plastic region and stress of the enclosing rock mass were obtained by introducing dimensionless parameters of radial distance. The results show that tunneling in deep rock mass causes a maximum stress zone to appear in the vicinity of the boundary of the elastic and the plastic zone in the surrounding rock mass. Under the compression of a large tangential force and a small radial force, the rock mass in the maximum stress zone was in an approximate uniaxial loading state, which could lead to a split failure in the rock mass.

Zonal disintegration phenomenon in rock mass surrounding deep tunnels

Zonal disintegration is a typical static phenomenon of deep rock masses. It has been defined as alternating regions of fractured and relatively intact rock mass that appear around or in front of the working stope during excavation of a deep tunnel. Zonal disintegration phenomenon was successfully demonstrated in the laboratory with 3D tests on analogous gypsum models, two circular cracked zones were observed in the test. The linear Mohr-Coulomb yield criterion was used with a constitutive model that showed linear softening and ideal residual plastic to analyze the elasto-plastic field of the enclosing rock mass around a deep tunnel. The results show that tunneling causes a maximum stress zone to appear between an elastic and plastic zone in the surrounding rock. The zonal disintegration phenomenon is analyzed by considering the stress-strain state of the rock mass in the vicinity of the maximum stress zone. Creep instability failure of the rock due to the development of the plastic zone, and transfer of the maximum stress zone into the rock mass, are the cause of zonal disintegration. An analytical criterion for the critical depth at which zonal disintegration can occur is derived. This depth depends mainly on the character and stress concentration coefficient of the rock mass.

Analysis of mechanical behavior of soft rocks and stability control in deep tunnels

Due to the weakness in mechanical properties of chlorite schist and the high in situ stress in Jinping II hydropower station, the rock mass surrounding the diversion tunnels located in chlorite schist was observed with extremely large deformations. This may significantly increase the risk of tunnel instability during excavation. In order to assess the stability of the diversion tunnels laboratory tests were carried out in association with the petrophysical properties, mechanical behaviors and waterlweakening properties of chlorite schist. The continuous deformation of surrounding rock mass, the destruction of the support structure and a large-scale collapse induced by the weak chlorite schist and high in situ stress were analyzed. The distributions of compressive deformation in the excavation zone with large deformations were also studied. In this regard, two reinforcement schemes for the excavation of diversion tunnel bottom section were proposed accordingly. This study could offer theoretical basis for deed tunnel construction in similar geological condition~

3D DEM simulation of hard rock fracture in deep tunnel excavation induced by changes in principal stress magnitude and orientation

To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.

Model test study of frost heaving pressures in tunnels excavated in fractured rock mass in cold regions

Based on the similarity theory, taking the horseshoe, city-gate and round linings as examples, the value and distribution regularities of normal frost heaving pressures （hereinafter as frost heaving pressures） in tunnels excavated in fractured rock mass in cold regions under different constraints and freezing depths were studied by a test model. It was found that the larger the frozen depth, the larger the frost heaving pressure, and the stronger the top constraint, the larger the frost heaving pressure. For the horseshoe lining and city-gate lining, the top constraint has a greater effect on the frost heaving pressures on the arch and the inverted arch. For the round lining, the influences of the top constraint on the frost heaving pressure in all linings are almost the same. The frost heaving pressure is maximum on the city-gate lining and minimal on the round lining. The largest frost heaving pressure all occur near the foot of the inverted arch for the three kinds of lining. Thus, the test data basically coincide with the observed in situ data.

Blasting induced dynamic stress concentration and failure characteristics of deep-buried rock tunnel

In this study,the dynamic stress concentration factors(DSCF)around a straight-wall arch tunnel(SWAT)were solved analytically utilizing the complex variable function methods and Duhamel’s integral.The effects of wavelength,incident angle,and blasting rising time on the DSCF distribution were analyzed.Theoretical results pointed out dynamic disturbances resulting in compressive stress concentration in the vertical direction and tensile stress in the incident direction.As the wavelength and rising time increased,there was a tendency for the amplitude of stress concentration to initially rise and then converge.Moreover,a series of 3D FEM models were established to evaluate the effect of different initial stress states on the dynamic failure of the tunnel surrounding rock.The results indicated that the failure of the surrounding rock was significantly influenced by the direction of the static maximum principal stress and the direction of the dynamic disturbance.Under the coupling of static and blasting loading,damage around the tunnel was more prone to occur in the dynamic and static stress concentration coincidence zone.Finally,the damage modes of rock tunnel under static stress and blasting disturbance from different directions were summarized and a proposed support system was presented.The results reveal the mechanisms of deep-buried rock tunnel destruction and dynamically triggered rockburst.

Design Theory of Full Face Rock Tunnel Boring Machine Transition Cutter Edge Angle and Its Application

At present, the inner cutters of a full face rock tunnel boring machine (TBM) and transition cutter edge angles are designed on the basis of indentation test or linear grooving test. The inner and outer edge angles of disc cutters are characterized as symmetric to each other with respect to the cutter edge plane. This design has some practical defects, such as severe eccentric wear and tipping, etc. In this paper, the current design theory of disc cutter edge angle is analyzed, and the characteristics of the rock-breaking movement of disc cutters are studied. The researching results show that the rotational motion of disc cutters with the cutterhead gives rise to the difference between the interactions of inner rock and outer rock with the contact area of disc cutters, with shearing and extrusion on the inner rock and attrition on the outer rock. The wear of disc cutters at the contact area is unbalanced, among which the wear in the largest normal stress area is most apparent. Therefore, a three-dimensional model theory of rock breaking and an edge angle design theory of transition disc cutter are proposed to overcome the flaws of the currently used TBM cutter heads, such as short life span, camber wearing, tipping. And a corresponding equation is established. With reference to a specific construction case, the edge angle of the transition disc cutter has been designed based on the theory. The application of TBM in some practical project proves that the theory has obvious advantages in enhancing disc cutter life, decreasing replacement frequency, and making economic benefits. The proposed research provides a theoretical basis for the design of TBM three-dimensional disc cutters whose rock-breaking operation time can be effectively increased.

Wear Analysis of Disc Cutters of Full Face Rock Tunnel Boring Machine

Wear is a major factor of disc cutters’ failure. No current theory offers a standard for the prediction of disc cutter wear yet. In the field the wear prediction method commonly used is based on the excavation length of tunnel boring machine（TBM） to predict the disc cutter wear and its wear law, considering the location number of each disc cutter on the cutterhead（radius for installation）; in theory, there is a prediction method of using arc wear coefficient. However, the preceding two methods have their own errors, with their accuracy being 40% or so and largely relying on the technicians’ experience. Therefore, radial wear coefficient, axial wear coefficient and trajectory wear coefficient are defined on the basis of the operating characteristics of TBM. With reference to the installation and characteristics of disc cutters, those coefficients are modified according to penetration, which gives rise to the presentation of comprehensive axial wear coefficient, comprehensive radial wear coefficient and comprehensive trajectory wear coefficient. Calculation and determination of wear coefficients are made with consideration of data from a segment of TBM project（excavation length 173 m）. The resulting wear coefficient values, after modification, are adopted to predict the disc cutter wear in the follow-up segment of the TBM project（excavation length of 5621 m）. The prediction results show that the disc cutter wear predicted with comprehensive radial wear coefficient and comprehensive trajectory wear coefficient are not only accurate（accuracy 16.12%） but also highly congruous, whereas there is a larger deviation in the prediction with comprehensive axial wear coefficient（accuracy 41%, which is in agreement with the prediction of disc cutters’ life in the field）. This paper puts forth a new method concerning prediction of life span and wear of TBM disc cutters as well as timing for replacing disc cutters.

Tunnel failure in hard rock with multiple weak planes due to excavation unloading of in-situ stress

Natural geological structures in rock(e.g.,joints,weakness planes,defects)play a vital role in the stability of tunnels and underground operations during construction.We investigated the failure characteristics of a deep circular tunnel in a rock mass with multiple weakness planes using a 2D combined finite element method/discrete element method(FEM/DEM).Conventional triaxial compression tests were performed on typical hard rock(marble)specimens under a range of confinement stress conditions to validate the rationale and accuracy of the proposed numerical approach.Parametric analysis was subsequently conducted to investigate the influence of inclination angle,and length on the crack propagation behavior,failure mode,energy evolution,and displacement distribution of the surrounding rock.The results show that the inclination angle strongly affects tunnel stability,and the failure intensity and damage range increase with increasing inclination angle and then decrease.The dynamic disasters are more likely with increasing weak plane length.Shearing and sliding along multiple weak planes are also consistently accompanied by kinetic energy fluctuations and surges after unloading,which implies a potentially violent dynamic response around a deeply-buried tunnel.Interactions between slabbing and shearing near the excavation boundaries are also discussed.The results presented here provide important insight into deep tunnel failure in hard rock influenced by both unloading disturbance and tectonic activation.

A novel image-based approach for interactive characterization of rock fracture spacing in a tunnel face

This paper presents a novel integrated method for interactive characterization of fracture spacing in rock tunnel sections.The main procedure includes four steps:(1)Automatic extraction of fracture traces,(2)digitization of trace maps,(3)disconnection and grouping of traces,and(4)interactive measurement of fracture set spacing,total spacing,and surface rock quality designation(S-RQD)value.To evaluate the performance of the proposed method,sample images were obtained by employing a photogrammetrybased scheme in tunnel faces.Experiments were then conducted to determine the optimal parameter values(i.e.distance threshold,angle threshold,and number of fracture trace grouping)for characterizing rock fracture spacing.By applying the identified optimal parameters involved in the model,the proposed method could lead to excellent qualitative results to a new tunnel face.To perform a quantitative analysis,three methods(i.e.field,straightening,and the proposed method)were employed in the same study and comparisons were made.The proposed method agrees well with the field measurement in terms of the maximum and average values of measured spacing distribution.Overall,the proposed method has reasonably good accuracy and interactive advantage for estimating the ultimate fracture spacing and S-RQD.It can be a possible extension of existing methods for fracture spacing characterization for two-dimensional(2D)rock tunnel faces.

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.

Numerical modeling of large deformation and nonlinear frictional contact of excavation boundary of deep soft rock tunnel

Roadways excavated in soft rocks at great depth are difficult to be maintained due to large deformation of surrounding rocks, which greatly influences the safety and efficiency of deep resources exploitation. During the excavation process of a deep soft rock tunnel, the rock wall may be compacted due to large deformation. In this paper, the technique to address this problem by a two-dimensional (2D) finite element software, large deformation engineering analyses software (LDEAS 1.0), is provided. By using the Lagrange multiplier method, the kinematic constraint of non-penetrating condition and static constraint of Coulomb friction are introduced to the governing equations in the form of incremental displacement. The numerical example demonstrates the efficiency of this technology. Deformations of a transportation tunnel in inclined soft rock strata at the depth of 1 000 m in Qishan coal mine and a tunnel excavated to three different depths are analyzed by two models, i.e. the additive decomposition model and polar decomposition model. It can be found that the deformation of the transportation tunnel is asymmetrical due to the inclination of rock strata. For extremely soft rock, large deformation can converge only for the additive decomposition model. The deformation of surrounding rocks increases with the increase in the tunnel depth for both models. At the same depth, the deformation calculated by the additive decomposition model is smaller than that by the polar decomposition model.

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.

The three-dimension model for the rock-breaking mechanism of disc cutter andanalysis of rock-breaking forces

To study the rock deformation with three- dimensional model under rolling forces of disc cutter, by car- rying out the circular-grooving test with disc cutter rolling around on the rock, the rock mechanical behavior under rolling disc cutter is studied, the mechanical model of disc cutter rolling around the groove is established, and the the- ory of single-point and double-angle variables is proposed. Based on this theory, the physics equations and geometric equations of rock mechanical behavior under disc cutters of tunnel boring machine （TBM） are studied, and then the bal- ance equations of interactive forces between disc cutter and rock are established. Accordingly, formulas about normal force, rolling force and side force of a disc cutter are de- rived, and their validity is studied by tests. Therefore, a new method and theory is proposed to study rock- breaking mech- anism of disc cutters.

Tunnelling performance prediction of cantilever boring machine in sedimentary hard-rock tunnel using deep belief network

Evaluating the adaptability of cantilever boring machine(CBM) through in-depth excavation and analysis of tunnel excavation data and rock mass parameters is the premise of mechanical design and efficient excavation in the field of underground space engineering.This paper presented a case study of tunnelling performance prediction method of CBM in sedimentary hard-rock tunnel of Karst landform type by using tunneling data and surrounding rock parameters.The uniaxial compressive strength(UCS),rock integrity factor(Kv),basic quality index([BQ]),rock quality index RQD,brazilian tensile strength(BTS) and brittleness index(BI) were introduced to construct a performance prediction database based on the hard-rock tunnel of Guiyang Metro Line 1 and Line 3,and then established the performance prediction model of cantilever boring machine.Then the deep belief network(DBN) was introduced into the performance prediction model,and the reliability of performance prediction model was verified by combining with engineering data.The study showed that the influence degree of surrounding rock parameters on the tunneling performance of the cantilever boring machine is UCS > [BQ] > BTS >RQD > Kv > BI.The performance prediction model shows that the instantaneous cutting rate(ICR) has a good correlation with the surrounding rock parameters,and the predicting model accuracy is related to the reliability of construction data.The prediction of limestone and dolomite sections of Line 3 based on the DBN performance prediction model shows that the measured ICR and predicted ICR is consistent and the built performance prediction model is reliable.The research results have theoretical reference significance for the applicability analysis and mechanical selection of cantilever boring machine for hard rock tunnel.