Hoek-Brown Failure Criterion-Based Creep Constitutive Model and BP Neural Network Parameter Inversion for Soft Surrounding Rock Mass of Tunnels

: This paper summarizes the main factors affecting the large deformation of soft rock tunnels, including the lithology combination, weathering effect, and underground water status, by reviewing the typical cases of largely-deformed soft rock tunnels. The engineering geological properties of the rock mass were quantiﬁed using the rock mass block index ( RBI ) and the absolute weathering index ( AWI ) to calculate the geological strength index ( GSI ). Then, the long-term strength σ r and the elastic modulus E 0 of the rock mass were calculated according to the Hoek–Brown failure criterion and substituted into the creep constitutive model based on the Nashihara model. Finally, the creep parameters of the surrounding rock mass of the Ganbao tunnel were inverted and validated by integrating the on-site monitoring and BP neural network. The inversion results were consistent with the measured convergence during monitoring and satisﬁed the engineering requirements of accuracy. The method proposed in this paper can be used to invert the geological parameters of the surrounding rock mass for a certain point, which can provide important mechanical parameters for the design and construction of tunnels, and ensure the stability of the surrounding rock mass during the period of construction and the safety of the lining structure during operation.


Introduction
Tunnel projects in western China often encounter soft rocks with well-developed bedding, such as carbonaceous phyllite, sericite phyllite, schist, carbonaceous slate, sandy slate, and carbonaceous shale. Under high in situ stress, the laminated soft surrounding rock masses are vulnerable to large and rapid deformation and local destruction [1][2][3]. In these conditions, tunnels frequently suffer from large deformation hazards.
Laminated rock mass, frequently seen in engineering practice and presenting oriented grouped bedding, has more heterogeneous mechanical properties than normal rock mass. Many researchers have developed constitutive models for laminated surrounding rock masses. For instance, Jia et al. [4] applied the microscopic element method with the constitutive model based on damage mechanics and statistical theory to the simulation of rock tunnel stability using the finite element method (FEM). Li et al. [5] developed a three-dimensional creep constitutive model for transversely isotropic rock mass, based on the Burgers viscoelastic model. Li et al. [6] proposed three basic creep patterns of shale and a general methodology for developing the anisotropic creep model. These constitutive models have been used to probe the deformation and failure characteristics of laminated rock masses with varied dip angles. However, research on the constitutive theory of laminated surrounding rock masses in tunnels is fairly limited and cannot provide the theoretical guidance for the engineering design of actual tunnels. In terms of the mechanisms of laminated rock mass deformation, previous studies have mainly investigated patterns of stability [7], deformation [8][9][10], mechanical behavior [11], and mechanical properties [12]. The main triggers for large tunnel deformation through highlydipped laminated soft rock under high in situ stress are high structural stress, unfavorable rock occurrences, and low rock strength [13]. On-site tunnel monitoring has shown that the deformation-failure zone of the surrounding rock mass of laminated soft rock tunnels is concentrated along the direction perpendicular to the rock bedding, instead of the direction of the maximum principal stress; the local deformation of laminated soft rock tunnels is affected by the topography, rock mass structure, and in situ stress [14]. Four mechanisms have been proposed to explain large tunnel deformation, namely, along-bedding sliding, flexural deformation, toppling deformation, and plastic extrusion. For each mechanism, a specific anchoring support plan should be developed [15]. Numerical simulation has been performed to comprehensively analyze the stress environment, the deviator stress distribution [16][17][18][19], and the plastic zone distribution [20,21] of the surrounding rock mass during tunnel excavation. Such simulations, incorporating the rock characteristics of deterioration, stabilization, and accelerated creep [22], reveal the characteristics of the surrounding rock mass failure with varied surrounding rock mass strengths, buried depths, and roof strengths [23], which provide guiding values to support the design of the analogous tunnels.
Due to the complexity of geology, geotechnical materials are typically characterized by discontinuity and heterogeneity, and the rheological parameters of the actual surrounding rock mass of the tunnel are often hard to measure. There are two main methods to obtain creep parameters: one is to calculate creep parameters through laboratory testing of rock combined with a creep constitutive model, but there is a large deviation between the laboratory results for rock and on-site rock mass parameters. Second, to obtain the mechanical parameters of the surrounding rock mass through displacement back analysis based on monitoring data [24]. Wenzheng Cao [25] proposed a novel back analysis program based on a BP neural network, which can realize automatic correction and the adjustment of parameters and adapt to most tunnel projects. Qingdong Wu [26] introduced a support vector machine (SVM) and an artificial neural network (ANN) to predict tunnel surrounding rock mass displacement, and compared and analyzed the results of the two methods. Xianghui Deng [27] established a tunnel risk assessment model by combining a fuzzy method with a BP neural network based on historical data from 50 tunnels. Numerical simulation has been widely combined with back-propagation (BP) neural network analysis to calculate the physical and mechanical parameters of surrounding rock mass, which has achieved good application performance [28][29][30].
This paper proposed a solution for obtaining the creep parameters of the surrounding rock mass in highway tunnels after excavation. First, the long-term strength of the engineering rock mass was obtained using the Hoek-Brown failure criterion, and a creep constitutive model for the rock mass was developed. Then, back-propagation analysis was performed using the on-site monitoring data to obtain the initial values of the creep parameters of the tunnel's surrounding rock mass. Subsequently, training samples were generated via orthogonal experiments. Then, a BP neural network and numerical simulation were integrated to inversely calculate the tunnel's surrounding rock mass creep parameters that matched the reality.

Factors Affecting the Deformation of Soft Rock Tunnels under High In Situ Stresses
The engineering data from some typical largely-deformed soft rock tunnels in China are summarized in Table 1, including the buried depth, formation lithology, in-situ stress, deformation magnitude, and characteristics of deformation and failure. Surrounding rock masses that are found with large deformations are mostly soft rock with a highly developed bedding texture, such as shale, slate, and phyllite ( Table 1). The bedding plane profoundly affects the strength and deformation characteristics of tunnel-surrounding rock mass and is considered the key factor determining the mechanical behavior and deformation-failure characteristics of anisotropic rock mass. In general, the compressional deformation of the surrounding rock mass of tunnels in laminated soft rock occurs mainly at the two sides along the normal direction of the weakest bedding plane. Moreover, the smaller the bedding spacing is, the worse the integrity and self-stabilizing ability of the surrounding rock mass are, which leads to increased proneness to local large deformation and thus a relatively large compressional load on the support structure and highly-deformed lining and surrounding rock mass in local areas. The factors affecting the large deformation of the surrounding rock mass can be summarized as follows: Based on lithology, the surrounding rock mass can be grouped into two types, namely, plastic and brittle rock. The former consists of soft rock and is typically characterized by low mechanical strength, softening and swelling on wetting, and low slaking resistance, which are all unfavorable for the load-bearing and stability of the tunnel's surrounding rock mass. On the contrary, the latter is of hard rock and generally has a high load-bearing capacity, low softening and swelling tendencies on wetting, and high slaking resistance, which are favorable for the load-bearing and stability of the tunnel's surrounding rock mass. The lithology combination can be generally summarized into four cases, namely, the consistent lithology, hard interbeds (in much softer surrounding rock mass), soft interbeds (in much harder surrounding rock mass), and soft-hard alternating lithology. Different lithology combinations are associated with varied deformation strengths. A soft interbed in the surrounding rock mass is equivalent to a weak part in the surrounding rock mass. In most cases, deformation-failure occurs first in the soft interbed during tunnel excavation. If a hard interbed exists in the surrounding rock mass, it will often suppress the rock deformation, due to its more rigid mechanical properties than those of the surrounding rock mass. Therefore, the deformation magnitude is often relatively small at the hard interbed.
The weathering effect is a common geological phenomenon in nature, which can weaken the cementation between rock particles, form damage cracks, reduce surface roughness, and worsen the physical and mechanical properties of rock. Rapid weathering of soft rock in tunnels will reduce the stability of the surrounding rock mass and affect the safety of the tunnel construction. After the tunnel excavation, the weak layered surrounding rock mass will form internal cracks due to the unloading effect; then, an effective seepage channels are formed, groundwater seeps through these channels, reducing the effective stresses and thus the shear strength along discontinuities and therefore the strength of the rock mass. Therefore, the presence of underground water makes the rock mass more prone to large deformation. The varied underground water seepage pathways or inconsistent thickness of the broken rock zone in the surrounding rock mass after excavation may result in local large deformation of the surrounding rock mass.
Due to the complex geological conditions outlined above, obvious asymmetry is observed for the deformation occurring at the two sides of the central axis, crown, and inverted arch, and local deformation is a common phenomenon for the practice of tunnel engineering. Therefore, it is necessary to quantitatively analyze the engineering geology of each part of the tunnel's surrounding rock mass and accurately invert the creep parameters at each local point from the on-site deformation monitoring data, so as to guide the design and construction of tunnels.

Rock Strength and Creep Constitutive Model Based on the Hoek-Brown Failure Criterion
The classic rock mass failure criterion considers rock as a continuous homogeneous medium, which greatly simplifies the reality. A constitutive model for the rock mass with multiple structural planes was developed via superimposition of the constitutive equations of the rock mass with a single joint, which cannot accurately capture the complex mechanical characteristics of the tiny fissures and multiple groups of structural planes in the rock mass [31,32]. Hence, the semi-empirical semi-theoretical approach was used to characterize the laminated rock mass in engineering practice.
Hoek and E. T. Brown derived the correlation expression among the limiting principal stresses for the failure of rock (mass) (namely, the Hoek-Brown failure criterion) via a trial-and-error process, based on the statistical analysis of massive data from rock triaxial tests and field testing of rock mass, and the Griffith theory. Then, E. Hoek further proposed the generalized Hoek-Brown empirical failure criterion. The Hoek-Brown empirical failure criterion can reflect the effects of various factors of the rock mass-such as the rock strength, the number of structural planes, and the in situ stress-on the rock mass strength, and overcome the disadvantages of the conventional theoretical equation. The generalized Hoek-Brown failure criterion can be expressed as below: where σ r is the maximum principal stress when the rock fails, namely, the rock strength; σ 3 is the minimum principal stress when the rock fails; σ c is the uniaxial compressive strength of the intact rock; m b , s, and a are all empirical parameters; m b represents the hardness of the rock, with a value of 0.0000001-25 (0.0000001 for severely-disturbed rock and 25 for intact hard rock); s represents the broken degree of the rock mass, with a value of 0-1 (zero for completely broken rock and one for intact rock); a is related to the rock mass quality. For rock mass with good quality, its strength characteristics are mainly controlled by the rock particle strength, due to the tight packing of rock particles. The restricted Hoek-Brown empirical failure criterion is more applicable to this case and a = 0.5. On the contrary, for rock mass with poor quality, the packing of fragments in the rock mass is loosened by shearing or weathering, which results in loss of the tensile strength of the rock mass (zero cohesion). Under such circumstances, if there is no confinement, the rock mass will collapse and a should be assigned to other values. The expressions of the parameters in Equation (1) are listed below: where D is the disturbance factor (zero for the undisturbed rock mass, and one for completely disturbed rock mass). m i indicates the values of constants for intact rock, which can be determined by laboratory testing (uniaxial compression test and conventional triaxial compression test). Based on laboratory data and engineering experience, E. Hoek et al. [33,34] developed a comprehensive and detailed m i value table that was able to cover a variety of rocks; the m i index of intact phyllite in this paper refers to this m i value table.
Hoek and E. J. Brown further derived the estimations of the relevant mechanical parameters of the rock mass [35], based on Equation (1): Equation (1) shows that the GSI value is the key to determining the mechanical parameters of the Hoek-Brown failure criterion for the rock mass. However, E. Hoek proposed only a general range for the texture type and weathering conditions of the rock mass and offered no quantitative calculation method.
Rock quality designation (RQD) is the most commonly used index to describe the structural characteristics of rock mass, which is defined as the percentage of the sum of the lengths of intact core pieces longer than 0.1 m relative to the chosen length of an evaluated drill core (RQD = ( the lengths of intact core pieces ≥ 10 cm per footage the length of an evaluated drill core (%)); however, the definition of RQD requires that the quality of core pieces must be "hard and sound", meaning it can only calculate the percentage of the core with a length beyond 10 cm in the total cored length, which is not suitable for integrity evaluations when characterizing the multi-joint rock mass.
This paper introduces the rock mass block index (RBI) and the absolute weathering index (AWI) to quantify the GSI. The RBI, proposed by Hu et al. [36], can thoroughly characterize the block dimension, texture type, and structural packing of blocks for the rock mass. The RBI defines the percentages of drilling cores with the measured lengths of 3-10 cm, 10-30 cm, 30-50 cm, 50-100 cm, and >100 cm as the weights, and the value of the RBI is the sum of the weights multiplied by the corresponding coefficients.
Although the rock mass has the same RQD value, the RBI values can be different. For example, the rock mass of RQD = 90% can be 10~30 cm with a mosaic structure, or 30~50 cm with a block structure, or 50~100 cm with a block structure, or even more than 100 cm with intact structure. Under the same RQD, the larger the RBI is (from a mosaic structure to an intact structure), the more intact the rock mass is, and the RBI can be regarded as an extension of RQD. The calculation formula of RBI is shown below: where C r3 , C r10 , C r30 , C r50 , and C r100 are the acquisition rates (weights in percentages) of the cores with the lengths of 3-10 cm, 10-30 cm, 30-50 cm, 50-100 cm, and >100 cm, respectively. The details of the rock mass texture representation by RBI, associated with the field description, are shown in Table 2. The mechanical properties of rock are also related to the weathering level. Chemical weathering can alter the mineral composition of rock and thus alter its physical and mechanical performance.
Parker [37] proposed an index of weathering for silicate rocks, based on the proportions of the alkali and alkaline earth metals present. The bond strengths of the most mobile of the major elements with oxygen are used as weighting factors in the index; the proposed index is defined by the following expression: where a(X) (X = Na, Mg, K, Ca) indicates the atomic proportion of element X, defined as the atomic percentage divided by the atomic weight, and the denominator is the bond strength of element X with oxygen, which represents the stability of the element in the weathering process. From fresh rock to weathered rock, the index WI is gradually reduced. The weathering index of rock mass AW I 0 [38] can be expressed as below: where WI stands for the weathering index of the weathered rock, whereas WI represents that of the fresh rock. Thus, a higher value of AWI indicates fresher rock with a lower level of chemical weathering and correspondingly better mechanical properties. On the contrary, a higher AWI value denotes a high level of chemical weathering. Obviously, AWI quantitatively characterizes the weathering conditions of the engineering rock mass. Based on the research on the relationship between the weathering index and the weathering level of the rock [39], the weathering level of rock based on AWI is shown in Table 3/ Figure 1. In order to conveniently determine the absolute weathering index AW I 0 of the rock mass, can use the ratio of uniaxial compressive strength of fresh rock to weathered rock. A correction factor λ is introduced to describe the effects of underground water on the weathering conditions of the jointed rock mass (Table 4). Accordingly, the rock mass weathering index, considering the effects of underground water, can be expressed as = . A correction factor λ is introduced to describe the effects of underground water on the weathering conditions of the jointed rock mass (Table 4). Accordingly, the rock mass weathering index, considering the effects of underground water, can be expressed as AW I = λAW I 0 . Table 4. Recommended values for the correction factor λ to capture the effects of underground water on the weathering conditions of rock mass.  The widely applied Nishihara model can characterize the decelerating, stead and accelerating creep phases of rock. Extensive experiments have shown that damage of the rock occurs only in the accelerating creep phase. The rock stren presented upon the failure of the rock, is the long-term strength for rock creep. Th elements of the rock mass model are shown in Figure 3. The creep equation can be expressed as below:

Production Status of Underground Water
where is the elastic modulus; is the viscoelastic modulus; and viscosity coefficients; is the load stress; is the total strain; is the lo strength calculated using Equation (1) (the Hoek-Brown failure criterion); and creep time.
The strength parameters of the laminated jointed rock mass at each position 4) can be obtained according to the Hoek-Brown failure criterion. Then, they ar The creep equation can be expressed as below: where E 0 is the elastic modulus; E 1 is the viscoelastic modulus; η 1 and η 2 are the viscosity coefficients; σ is the load stress; ε is the total strain; σ r is the long-term strength calculated using Equation (1)  into the numerical simulator via the creep constitutive model (Equation (12)) to calculate the deformation of the surrounding rock mass ( Figure 5).

Field Application
The stress state of the surrounding rock mass changes after starting the excavation of the tunnel. Hence, the support design based on the surrounding rock mass parameters obtained in the previous investigation is inappropriate, in particular for a tunnel engineering design with a high buried depth, large span, and highly-developed structural  (12)) to calcula the deformation of the surrounding rock mass ( Figure 5).

Field Application
The stress state of the surrounding rock mass changes after starting the excavation the tunnel. Hence, the support design based on the surrounding rock mass paramete

Field Application
The stress state of the surrounding rock mass changes after starting the excavation of the tunnel. Hence, the support design based on the surrounding rock mass parameters obtained in the previous investigation is inappropriate, in particular for a tunnel engineering design with a high buried depth, large span, and highly-developed structural planes. To provide a basis for theoretical reference and design, this paper inverts and validates the surrounding rock mass parameters by integrating the on-site monitoring with a BP neural network.
The BP neural network is a typical nonlinear algorithm, composed of the input, output, and several (one or more) hidden layers (Figure 6), and each layer has several nodes. The connection of nodes between layers is represented by the weight. The BP neural network with one hidden layer is the traditional shallow neural network, whereas that containing multiple hidden layers is the deep learning neural network (Figure 7).  The core steps for training the BP neural network are illustrated below which the solid line represents forward propagation and the dashed line re propagation. The forward propagation means that the data (information first imported into the input end, delivered along the network direction, a by the corresponding weight. The products of the input data and the c weights are then summed, and the results are input into the activation calculation. Then calculation results are delivered to the next node as t calculation is performed successively until the final output is obtained.  The core steps for training the BP neural network are illustrated below which the solid line represents forward propagation and the dashed line r propagation. The forward propagation means that the data (information first imported into the input end, delivered along the network direction, a by the corresponding weight. The products of the input data and the weights are then summed, and the results are input into the activatio calculation. Then calculation results are delivered to the next node as calculation is performed successively until the final output is obtained. The core steps for training the BP neural network are illustrated below (Figure 8), in which the solid line represents forward propagation and the dashed line represents back propagation. The forward propagation means that the data (information or signal) are first imported into the input end, delivered along the network direction, and multiplied by the corresponding weight. The products of the input data and the corresponding weights are then summed, and the results are input into the activation function for calculation. Then calculation results are delivered to the next node as the input. The calculation is performed successively until the final output is obtained.
propagation. The forward propagation means that the data (information or signal) first imported into the input end, delivered along the network direction, and multipl by the corresponding weight. The products of the input data and the correspond weights are then summed, and the results are input into the activation function calculation. Then calculation results are delivered to the next node as the input. T calculation is performed successively until the final output is obtained. Here we summarize the inversion workflow of the creep parameters of surrounding rock mass, as shown below (Figure 9). Here we summarize the inversion workflow of the creep parameters of the surrounding rock mass, as shown below (Figure 9).  The Ganbao tunnel project of the Wenchuan-Maerkang expressway was taken as an engineering case study. It is an ultra-long left-right separated tunnel project (the left tunnel is 4777 m long and the right one is 4796 m long). The longitudinal profile of the tunnel is illustrated in Figure 10, in which the red surrounding rock mass mainly represents the sericite phyllite; the yellow part represents the sericite phyllite interbedded by carbonaceous phyllite and metasandstone; the blue part shows the carbonaceous phyllite and sericite phyllite interbedded by metasandstone; the green part represents the sericite phyllite and carbonaceous phyllite; and the black part indicates the carbonaceous phyllite. The tunnel was excavated using the drilling and blasting method. A large deformation occurs at the mile sections K141 + 540~K141 + 860 and K143 + 485~K143 + 650, which results in the circumferential cracking of multiple points along the initially supported arch and local cracking and spalling of the sidewall and tunnel crown ( Figure  11). The monitoring section is located at K143 + 621 and starts monitoring after 12 h of The Ganbao tunnel project of the Wenchuan-Maerkang expressway was taken as an engineering case study. It is an ultra-long left-right separated tunnel project (the left tunnel is 4777 m long and the right one is 4796 m long). The longitudinal profile of the tunnel is illustrated in Figure 10, in which the red surrounding rock mass mainly represents the sericite phyllite; the yellow part represents the sericite phyllite interbedded by carbonaceous phyllite and metasandstone; the blue part shows the carbonaceous phyllite and sericite phyllite interbedded by metasandstone; the green part represents the sericite phyllite and carbonaceous phyllite; and the black part indicates the carbonaceous phyllite. The tunnel was excavated using the drilling and blasting method. A large deformation occurs at the mile sections K141 + 540~K141 + 860 and K143 + 485~K143 + 650, which results in the circumferential cracking of multiple points along the initially supported arch and local cracking and spalling of the sidewall and tunnel crown ( Figure 11). The monitoring section is located at K143 + 621 and starts monitoring after 12 h of excavation. The overburden near the monitoring section consists of collapsed debris deposits, sericite phyllite and carbonaceous phyllite (from top to bottom, Figure 12),the lithology presented on the face of the surrounding rock mass of the largely-deformed tunnel section ( Figure 13) is dominated by the carbonaceous phyllite and sericite phyllite, mainly the Grade-V surrounding rock mass. The rock mass is found with a laminated texture and has a buried depth of 646 m, with a vertical stress of 21.7 MPa and a horizontal stress of 8.3 MPa. Field observations indicate minimal effects of underground water, and thus the correction factor λ for the effects of underground water on the weathering conditions of the rock mass is set as one.     The pre-excavation boreholes are arranged on the tunnel face (Figure 14), the rock cores are collected via pre-excavation boreholes for investigation, and the RBI value of rock mass (Table 5) for each borehole is calculated using Equation (9) and the core length measured on site. For the monitoring of surrounding rock mass deformation after excavation, the on-site monitoring point placement is illustrated in Figure 14. The pre-excavation boreholes are arranged on the tunnel face (Figure 14), the rock cores are collected via pre-excavation boreholes for investigation, and the RBI value of rock mass (Table 5) for each borehole is calculated using Equation (9) and the core length measured on site. For the monitoring of surrounding rock mass deformation after excavation, the on-site monitoring point placement is illustrated in Figure 14. The pre-excavation boreholes are arranged on the tunnel face (Figure 14), the rock cores are collected via pre-excavation boreholes for investigation, and the RBI value of rock mass (Table 5) for each borehole is calculated using Equation (9) and the core length measured on site. For the monitoring of surrounding rock mass deformation after excavation, the on-site monitoring point placement is illustrated in Figure 14.

Cores
Taking Point 4 as an example, the field investigation shows that the R surrounding rock mass cores at Point 4 is about 2.8; the AWI is about 0.65. Acc Table 4, the GSI value is 33. For the intact phyllite rock, the Hoek-Brown cons 10 and the uniaxial compressive strength = 35 MPa. From Equation (5),     Table 6.  Taking Point 4 as an example, the field investigation shows that the RBI of the surrounding rock mass cores at Point 4 is about 2.8; the AWI is about 0.65. According to Table 4 Table 6.   Table 6. Table 6. Basic parameters of the surrounding rock mass at each monitoring position.  Taking Point 4 as an example, the field investigation shows that the RBI of the surrounding rock mass cores at Point 4 is about 2.8; the AWI is about 0.65. According to Table 4 Table 6. Table 6. Basic parameters of the surrounding rock mass at each monitoring position.   Table 6. Table 6. Basic parameters of the surrounding rock mass at each monitoring position.

4.331
Taking Point 4 as an example, the field investigation shows that the RBI of the surrounding rock mass cores at Point 4 is about 2.8; the AWI is about 0.65. According to Table 4 Table 6. The fitting of the creep constitutive model (Equation (11)) to the surrounding rock mass monitoring data ( Figure 15) was performed using the least square method to calculate the creep parameters of the rock mass at each monitoring point. The calculation results are shown in Table 7. The fitting of the creep constitutive model (Equation (11)) mass monitoring data ( Figure 15) was performed using the calculate the creep parameters of the rock mass at each monitori results are shown in Table 7.    Since the factors affecting the deformation of the surrounding rock mass are varied in the different positions of the tunnel, a neural network was constructed and trained for each point to improve the inversion accuracy. The single-variable method was used for training, which means when the neural network of one point is being trained, the parameters of the other points are constant. Here, Point 4 was taken as an example to demonstrate the training of the neural network. As stated above, the parameters of the other points were fixed, and the learning sample dataset of the surrounding rock mass mechanical properties was built based on the orthogonal test design (Figure 16). The sample parameters in the training dataset were substituted into the numerical simulation model developed using Midas (Figure 17), and the forward calculation was performed to obtain the surrounding rock mass deformation magnitude at each point (Table 8). We used 80% of the data in Table 8 as the training dataset, whereas we used the other 20% as the testing dataset.
ppl. Sci. 2021, 11,10033 by sericite and carbonaceous phyllite, were: average density ρ = 3.42 g/cm modulus E = 5.52 GPa, Poisson ratio ν = 0.38, internal cohesion c = 175 kPa angle φ = 27°.         Section K143 + 621 was chosen for numerical simulation. Its buried depth was 646 m, we set the size of the model as long × wide = 360 m × 150 m, the fixed end restraint was set at the bottom, the overlying pressure (γH = 33.6 KN/m 3 × (646 − 150) m = 16.67 MPa) was added above the model, the horizontal stress (8.3 MPa) was added on both sides of the model, and the surrounding rock masses all adopted the M-C failure criterion, as shown in Figure 17. The physical parameters of the surrounding rock mass, dominated by sericite and carbonaceous phyllite, were: average density ρ = 3.42 g/cm 3 , average elastic modulus E = 5.52 GPa, Poisson ratio ν = 0.38, internal cohesion c = 175 kPa, internal friction angle ϕ = 27 • .

Displacement Displacement Displacement Displacemen
Using the BP neural network toolkit in MATLAB, the forward training of the rock mass parameters and deformation at Point 4 (Table 7) was performed, which produced the neural network model of Point 4, referred to as net4. Similarly, the training samples for the other point were generated and so were their neural network models, namely net1, net2, net3, and net5, respectively. Finally, the deformation measured at each point was input into the corresponding neural network model and the surrounding rock mass parameters were inverted. For the purpose of validation, the inverted mechanical parameters of the surrounding rock mass were substituted into the Midas numerical model to determine the deformation at each point of the surrounding rock mass (Table 9). The relative error between the measured displacement and the numerically simulated displacement obtained through inversion was within 15%, indicating the good consistency and high applicability of the proposed method. Therefore, one can build a monitoring system for the largely-deformed surrounding rock mass section, according to the onsite information (such as the lithology combination, underground water status, bedding thickness, rock mass texture, surrounding rock mass integrity, maximum principal stress, and on-site measured rock deformation) for the purposes of gaining real-time information of the tunnel's surrounding rock mass parameter, adjusting the support plan in a timely manner, and providing references for design and construction.

Conclusions
Due to the concealment of tunnel engineering, it is difficult to obtain the parameters of the surrounding rock mass after excavation. In this study, we combined the traditional Nashihara model with the Hoek-Brown failure criterion and developed the creep constitutive model for the jointed rock mass. Then, the creep parameters of the surrounding rock mass of the Ganbao tunnel were inverted and validated by integrating the on-site monitoring and the BP neural network. The following conclusions were drawn: (1) By reviewing the typical cases of large deformation in soft rock tunnels, the main influential factors can be summarized as the lithology combination, weathering effect, and underground water status. With the classical rock mass failure criterion, it is hard to thoroughly incorporate the geological characteristics of the actual rock mass and therefore the semi-empirical semi-theoretical Hoek-Brown approach is more fit-for-purpose. (2) The geological characteristics of the engineering rock mass were quantitatively characterized using two indexes, namely, the rock mass block index (RBI) and the absolute weathering index (AWI). Following the Hoek-Brown criterion, the long-term strength σ r and elastic modulus E 0 of the rock were obtained and then substituted into the rock mass creep constitutive model based on the Nashihara model and the Hoek-Brown failure criterion. By doing so, the original five creep parameters that needed to be determined in the creep equation were reduced to three, which simplifies the calculation of the constitutive model, while well reflecting the engineering, geological, and creep characteristics of the rock mass on site. (3) Considering the fact that the actual engineering geology varies at different positions of the tunnel's surrounding rock mass, a specific BP neural network model was built for each monitoring point. Then, the rock mass parameters at each point were inverted from the on-site measured deformation. At last, the inverted parameters were input into the numerical model to calculate the deformation at each point, which was compared with the corresponding measured deformation. The resultant errors were all within 15%, satisfying the engineering requirements and demonstrating the reliability of the proposed method. (4) The values of the rock mass GSI (Figure 2) were all determined according to the engineering geology handbook, relevant standards, and existing literature. Therefore, these values can be adjusted as per the field condition. The inversion of the surrounding rock mass parameters is highly affected by the basic parameters of the engineering geological characteristics and the in-situ stress field of the tunnel's surrounding rock mass. Hence, during applications of the proposed method, the basic parameters need to be accurately measured to improve the accuracy of the inversion. (5) This method is suitable for tunneling in unsupported rock mass or plainly supported tunnels after excavation. Institutional Review Board Statement: Not involved.

Informed Consent Statement: Not involved.
Data Availability Statement: The study did not report any data.