Abstract
Uncertainty is an essential property of rock mechanics and engineering, which is of great significance to excavation, design, and control of rock engineering. In this study, an innovative framework of the reliability-based design was developed for the rock tunnel under uncertainty. The convergence-confinement method is used to characterize the interaction mechanism between the support structure and surrounding rock mass. Artificial bee colony (ABC) was adopted to solve the optimization problem in the reliability-based design. The probabilistic properties of rock strength and failure envelope were obtained based on the triaxial compression test data using the Bayesian method. The reliability of the tunnel and support structure was evaluated based on the abovementioned probabilistic properties of rock strength using the reliability analysis method. A circular tunnel was used to illustrate the developed framework, and the procedure was presented in detail. The time of rockbolt installed, the thickness of the shotcrete, length of rockbolt, circumferential space, and longitudinal space of rockbolt were determined and met the constraints of reliability index. Results show that the developed framework can consider the uncertainty for support design in the tunnel. It provides a good and promising way to support design considering the uncertainty of test data using the reliability-based design.
1. Introduction
The support structure (such as rockbolt, shotcrete, and steel bar) is necessary to control the deformation and instability of surrounding rock mass in rock engineering. Support type and parameters are affected by many factors, such as the mechanical properties of rock mass, the failure mechanism of rock, structure of rock mass, and interaction mechanism between rock and support, etc. Among them, the mechanical and strength property of surrounding rock mass are the essential factor. Rock tests such as uniaxial and triaxial compression tests are essential to obtain the mechanical and strength properties of surrounding rock mass. However, test data are random and scattered due to the complexity of rock mass and geological conditions. For the traditional deterministic approach of support design, the mean value of rock mechanical parameters obtained based on test data was used to determine the support type and parameters for the tunnel. The deterministic approach widely adopted in rock engineering is impractical and does not meet tunnel excavation requirements in practice. It is a challenging task to determine the support design of the tunnel under uncertainty.
Rock test data is one of the most basic information of tunnel support design. Various methods of determining mechanical and strength parameters of the rock mass were developed based on test data for the last few decades [1–9]. However, these methods do not consider the uncertainty and discreteness of natural rock mass and geological conditions. The probabilistic method was commonly used to handle the uncertainty in engineering systems and test data [10–16]. At the same time, the interaction analysis between rock and support has also made significant progress [17]. However, most of the interaction analysis between rock and support is based on a deterministic approach that cannot handle the uncertainty of rock properties and engineering loading. In order to reduce the instability risk of surrounding rock mass in tunnel excavation, the probabilistic approach has been developed for the design and evaluation of tunnels [18–20]. Lv et al. (2011) evaluated the reliability of rock-support interaction using the response surface method [19].
Recently, uncertainty optimization design has been widely used in various engineering fields to improve the performance of engineering systems [21–25]. The reliability-based design (RBD) is an essential tool to consider the uncertainty and is paid more attention in engineering systems. In the last decades, significant progress has been made in geotechnical and geological engineering using RBD [26–30]. Oreste estimated the probabilistic distribution of the safety factors of the support structure and determined the support structure and their parameters for the rock tunnel by considering the uncertainty of rock properties and support material [31]. Connor Langford and Diederichs (2013) adopted the reliability-based design to determine the composite lining by combining a modified reliability analysis method, Monte Carlo simulation, and finite element method [28]. The RBD has some evident advantages over the traditional deterministic design approach. So far, the uncertainty of tunnel excavation, including probabilistic analysis of rock mechanical parameters, reliability analysis, and RBD, has gained more attention. However, no framework integrated the uncertainty and scattered test data in the support design of tunnel excavation. In this study, a novel framework was developed by integrating the uncertainty analysis and evaluation of rock test data, stability analysis of tunnel, interaction mechanism between rock and support, and optimal design into the process of support design of the tunnel.
Uncertainty of test data was estimated to obtain the statistical feature of mechanical and strength parameters of rock mass based on the Bayesian method and the Hoek–Brown failure criterion. The stability of the tunnel was evaluated using FORM based on the above obtained statistical feature of the rock mass. And then, the support design of the tunnel was implemented based on the convergence-confinement method. The structure of this study is organized as follows. The Hoek–Brown failure criterion, rock-support interaction model, reliability analysis method, and optimal method are introduced briefly. The idea and procedure of the developed framework are presented in detail. The proposed method is illustrated and developed using a circular tunnel and its rock test data of the triaxial compression test. Some conclusions are made in the Conclusions section.
2. Rock-Support Interaction Designs for the Rock Tunnel
Support and reinforcement structure (rockbolt, shotcrete, steel set, etc.) is a significant measure that controls the deformation and prevents instability of surrounding rock mass induced by excavation. The support design of the tunnel is to minimize the risk of potential failure in the process of tunnel excavation [18]. Rock-support interaction analysis is an essential step in the support design of the tunnel. The convergence-confinement method is widely used in tunnel design. Moreover, it provides a fundamental approach for evaluating the interaction mechanism between rock and support and determining the support structure and parameters. The convergence-confinement method includes the three essential components in rock-support interaction analysis (i.e., the rock response curve (RSC), the support characteristic curve (SCC), and the longitudinal deformation profile (LDP)). Figure 1 illustrates and presents the three curves in a circular tunnel with hydrostatic stress.

2.1. Rock Response Curve
The rock response curve is used to present the relationship between the fictitious internal pressure pi and the radial displacement of the tunnel wall ur (seen in Figure 1). The analytical solution can be obtained in a circular tunnel with a hydrostatic stress field. The rock response curve can be determined from the elastoplastic analysis of stresses and strains for the surrounding rock mass base on the analytical solution. In this study, the rock response curve was determined using the analytical solution of a circular tunnel based on the Hoek–Brown failure criterion [32] and the convergence-confinement method. The corresponding code is implemented through VBA in Microsoft Excel.
2.2. Support Characteristic Curve
A support characteristic curve is used to characterize the feature and deformation mechanism of the support structure. The support stiffness and maximum support pressure are necessary to determine the support characteristic curve. For anchored rockbolt, the support stiffness and maximum support pressure can be obtained based on the pullout tests using the following equations:where L, db, Eb, and Q are the length, the diameter, the elastic modulus, and the load-deformation constant for the anchored rockbolt, respectively. Tbf is the ultimate failure load obtained from a pullout test. Sc and SL are the circumferential and longitudinal spacing of rockbolt, respectively. In this study, an Excel spreadsheet and VBA code are finished to determine the support stiffness and the maximum support pressure of the anchored rockbolt.
For shotcrete, the support stiffness and the maximum support pressure can be obtained by the following equations:where Ec, , tc, and σc are the elastic modulus, the Poisson ratio, the thickness of shotcrete, and the uniaxial compressive strength of the concrete or shotcrete, respectively.
For the combined support system, the following equations are used to determine the stiffness ktotal and the maximum support pressure ptotmax:where k1 and k2 are the stiffness of support system 1 and support system 2, respectively.
2.3. Longitudinal Deformation Profile
The longitudinal deformation profile (LDP) is another essential element of the convergence-confinement method. It is used to explain how quickly the support structure begins to prevent the deformation of the surrounding rock mass behind the face of the tunnel [33]. The LDP should be obtained from the measured and monitored data in the process of tunnel excavation. However, it is not easy to obtain this information for rock tunnels in practice. So, the numerical model provides an excellent way to obtain the LDP. In this study, an empirical LDP expression was adopted to determine the LDP curve.where ur is the radial deformation of the tunnel wall, is the maximum value of the radial deformation of the tunnel wall, x is the distance between position and the tunnel face, and R is the radius of the tunnel.
In this study, a circular tunnel under the hydrostatic pressure field was used to illustrate the design process of interaction between rock and support for rock tunnels based on the analytical solution. The Hoek–Brown failure criterion was adopted to obtain the analytical solution of deformation in surrounding rock mass. The analytical solution of the tunnel was used to generate the RSC. SCC was determined based on equations (1)–(9) in Section 2.2. LDP was determined using Hoek’s empirical equation (equation (10)). A brief procedure of the convergence-confinement method can be found in the literature [33].
3. Reliability-Based Design for Support of the Tunnel
3.1. First-Order Reliability Method (FORM)
The reliability analysis is an essential step in the reliability-based design for tunnel support design. In this study, the FORM was adopted to evaluate the reliability of tunnel deformation based on the optimization technology.
Low and Tang proposed a practical FORM approach by transferring the Hasofer-Lind index into a constrained optimization problem based on the perspective of an expanding ellipsoid in the original space of the basic random variables [34]. However, it needs to compute the equivalent normal means and equivalent normal standard deviations of random variables. A new efficient FORM algorithm was developed to avoid the above computation by introducing the dimensionless number ni.where R is the correlation matrix of random variables and n is a column vector of ni. While the value of ni is varied during constrained optimization, the corresponding value of random variables xi can be calculated using the following equation:
In this study, the reliability index is calculated using equation (11) based on Excel Solver that is regarded as an optimization tool. The computation procedure of the reliability index is coded in Excel VBA, which is called the Excel Solver.
3.2. Artificial Bee Colony Algorithm
Karaboga and Basturk (2008) were first inspired by the intelligent foraging behavior of bee swarms to develop the artificial bee colony (ABC) algorithm for solving the unconstraint optimization problem on where that it shows superior performance [35]. In the idea and algorithm of ABC, the colony of artificial bees is composed of employed bees, onlooker bees, and scout bees. Firstly, the ABC algorithm generates the initial population of SN solutions/bees randomly and determines the value of the fitness function for every solution/bee.where x (i, j) is the candidate solution of optimization problem; i = 1, 2, …, SN/2 and SN/2 is the number of populations; j = 1, 2, …, D, where D is the number of dimensions in each solution; rand (0, 1) denotes a random number generated between [0, 1]; and are the maximum and minimum values of each solution.
After the initialization population is generated, employed bees generate the updated solution using the ABC algorithm and calculate the value of the fitness function. The updated solution can be obtained by the employed bees, which is defined as follows:where k is an integer chosen randomly from {1, 2, …, SN/2} and is different from i, j is also integer generated randomly from {1, 2, …, D}. φij is generated randomly in [−1, 1], controls the generation of neighbor food sources around x (i, j), and represents the difference of two food positions found by a bee.
Onlooker bees generated a new solution based on the information obtained from the employed bees and selected an updated solution based on the value of the fitness function obtained from employed bees. And then, onlooker bees determine the abandoned solution and label the employed bees as scout bees. The probability pi of each selected bee can be calculated using the following equation:where fitnessi is the value of fitness function for the ith solution.
A new solution was generated randomly based on the searching of scout bees in the predetermined searching ranges to obtain the global solution. Onlookers will delete the solution in the predetermined number of cycles if it cannot be improved further. A new solution x′(i, j), determined by the scout bees using equation (14), replaces the abandoned solution x (i, j). The new candidate solution , which is generated based on x (i, j), can be evaluated based on the comparison of the fitness function between x (i, j) and . The new solution will be added to the population, and the old solution will be deleted when the new solution is equal or better than the old solution. Otherwise, the old solution is not changed in the population.
ABC algorithm will repeat the above three phases till the global solution was generated. The detailed information of ABC beyond the scope of this study, its algorithm, and procedure can reference some literature about ABC [35].
3.3. Reliability-Based Design Using ABC
Let d represent the design variables that need to determine and optimize for the engineering system. C (d) represents the objective function of an optimization problem that drives the design process for the engineering structure. The constraints of the reliability-based design are that the reliability index for the ith failure mode is larger than a target reliability index, . Let represent the ith element of d, the minimum, and maximum permissible values for di, respectively. The reliability-based design is used to determine the design variables that can minimize the objective function, C (d), and satisfy all reliability constraints. The reliability-based design can be regarded as an optimization problem with constraints in the following form:where C (d) is the objective function such as costs, weight, etc., βi (d) and are the reliability index and reliability constraint of the ith failure model such as the failure of the support structure and the radial deformation of the surrounding rock mass in the tunnel, di is the design variables such as length of rockbolt and the thickness of shotcrete, etc., are the minimum and maximum values of the ith design variables, and nm and nd are the number of reliability constraints and design variables, respectively. The reliability-based design can consider the uncertainty of the engineering system by involving the reliability constraints and is markedly different from determination optimization.
The procedure of reliability-based design includes inner and outer loops. The reliability index is calculated by reliability analysis in the inner loop, and the design variables were determined by solving an optimization problem in the outer loop. In this study, the reliability index is calculated using Excel Solver based on FORM in the inner loop. Then, the global solution is determined by searching the space of design variables using the ABC algorithm in the outer loop.
4. Framework of Tunnel Support Design considering Data Uncertainty
Rock laboratory and field tests are very critical to the support design for a tunnel. However, rock test data is scattered and uncertain. It is rational to consider the uncertainty of the support design for the tunnel. In this study, a framework of support design considering uncertainty was developed through combining test data, rock failure criterion, Bayesian method, reliability analysis, reliability-base design, and ABC optimization method (Figure 2). Bayesian method was used to analyze the probabilistic property of the Hoek–Brown failure criterion and its parameters. The reliability model of the tunnel was built based on the above probabilistic data and FORM. Then, the support design was determined by a reliability-based design and ABC. It provides a complete procedure of support design from test data to the final design of the tunnel by considering the test data uncertainty.

4.1. Hoek–Brown Failure Criterion
Hoek and Brown (2019) developed the generalized Hoek–Brown failure criterion to estimate rock mass strength based on the nonlinear Griffith failure criterion [8].where σ1 and σ3 are the first and third principal stresses, respectively. σci is the unconfined compressive strength and mb, s, and a are the material constants of rock mass obtained based on the following equations:where GSI is the geological strength index that represents the engineering geology observations of the field in the failure criterion and D is a coefficient that represents the degree of disturbance of the rock mass induced by excavation. For the intact rock, the material constants are denoted by mi, s = 1, and a = 0.5. Hoek–Brown failure criterion was used to characterize the failure mechanism of the intact rock mass in the following form:
To obtain the Hoek–Brown failure criterion for intact rock mass, mi and σci need to be determined based on the rock test data.
4.2. Determination of the Hoek–Brown Failure Criterion and Its Uncertainty Using the Bayesian Method
To determine the Hoek–Brown failure criterion and its uncertainty, equation (19) was considered the regression equation for Bayesian regression [36]. We are interested in the predicting outcome principal stress σ1 as normally distributed observations with an expected value σµ that is a nonlinear function of two predictor material constants σc and mi of Hoek–Brown failure criterion.
We will apply normal distribution priors with mean of Mσc and variance of Sσc to material constants σc and uniform [0, Vmi] to material constants mi, which corresponds to weak information regarding the actual parameter values.
Specifying the above model using the Bayesian method, the posterior estimates for the unknown material constant σc and mi can be obtained based on the Hoek–Brown failure criterion. In this study, MCMC sampling methods were used to determine the regression coefficient of the Hoek–Brown failure criterion.
4.3. Procedure of the Support Design of the Tunnel considering the Uncertainty
Once finished rock test and obtained the data of intact rock mass, the probabilistic property of the Hoek–Brown failure criterion and its parameters will be obtained based on the Bayesian method. Then, FORM was used to calculate the reliability of the tunnel. The reliability-based design of the tunnel support will be implemented by combining the convergence-confinement method, FORM, and ABC. The detailed procedures are as follows (Figure 3): Step 1. Implement the rock test and obtain the test data. Step 2. Compute the probabilistic properties of rock strength and material constant based on the Bayesian method and test data. Step 3. Build the analytical or numerical model of the tunnel. Step 4. Determine the RSC, SSC, and LDP for the tunnel. Step 5. Build the reliability model of the tunnel with support using FORM and Solver. Step 6. Determine the ABC algorithm and its parameters. Step 7. Determine the support design of the tunnel using an ABC-based reliability-based design based on the above data and model.

5. Application
A circular tunnel, where the surrounding rock mass is Carrara marble, was used to illustrate the developed support design framework considering uncertainty. The radius of the tunnel is 5 m. In situ stress is 30 MPa. The rock mass is assumed to isotropic elastoplastic. Hoek–Brown criterion and associated flow rules were adopted in the analytical model [37]. Fourteen test data of surrounding rock are listed in Table 1. Figure 4 shows that the test data is scattered. The GSI of the rock mass is 50. The Poisson ratio and density of rock mass are 0.25 and 0.027 kN/m3, respectively.

(a)

(b)
The strength parameters and material constants of rock mass and their uncertainty were estimated based on the Hoek–Brown failure criterion and the Bayesian method. The mean value of σi and mi are 92.7430 MPa and 7.9928, respectively. The standard derivation of σi and mi are 8.4160 MPa and 1.5523, respectively. The statistics property of rock strength parameters is shown in Figure 4. Figure 5 shows the failure envelope of rock mass based on the Hoek–Brown failure criterion.

Once obtaining the statistical property of the surrounding rock mass, a reliability analysis of the tunnel can be implemented using FORM based on the proposed procedure. In this study, only the rock strength parameters are uncertain. The objective function is the total weight of rockbolt in this study. The constraints of reliability are the reliability index of tunnel wall displacement, bearing capacity of rockbolt, and bearing capacity of shotcrete, respectively. The combined support of rockbolt and shotcrete was adopted to enforce the surrounding rock mass for the tunnel. The diameter of rockbolt db is 25 mm. The Young modulus Eb is 207 GPa. The ultimate failure load Tbf and load-deformation constant Q of rockbolt are 0.285 MN and 0.143 m/MN, respectively. The Young modulus Ec, uniaxial compressive strength σc, and Poisson’s ratio ν of shotcrete is 20.7 GPa, 34.5 MPa, and 0.25, respectively. The initial deformation before the rockbolt is installed, the thickness of the shotcrete, length of rockbolt, circumferential space, and longitudinal space of rockbolt need to be determined using the proposed framework. The ranges of searching and the spreadsheet of reliability-based design are shown in Figure 6. ABC was regarded as an optimal method in reliability-based design for the tunnel. Figure 6 lists the parameters of ABC. The initial deformation before rockbolt was installed, length of rockbolt, circumferential space, longitudinal space of rockbolt, and lining thickness are listed in Figure 6, respectively. Figure 7 shows the curve of rock-support interaction. The deformation of rock mass, rockbolt, and lining are safe and meet the request of the tunnel.


6. Conclusions
Support and enforcement structures are essential to the tunnel and are commonly used to control surrounding rock mass deformation and failure. Support design depends on rock strength parameters, in situ stress, geography conditions, support structure parameters, etc. However, there are amounts of uncertainty in tunnel engineering. The deterministic design method cannot meet the requirements of the tunnel with the development of the tunnel structure. In this study, a framework of support design was developed by combining the convergence-confinement method and ABC to determine the support parameters based on the above reliability model. The time of rockbolt installed, the thickness of the shotcrete, length of rockbolt, circumferential space, and longitudinal space of rockbolt were determined by the developed method. The results show that the developed method provided an excellent way to determine the support design with uncertainty:(1)The uncertainty is a vital risk resource of geotechnical and geological engineering and affects the stability and support design. Rock test is the primary uncertainty resource, and test data are scattered and uncertain. The proposed framework provides a tool for considering the uncertainty from the rock test for tunnel design. The statistical property of rock strength parameters can be determined rationally using the Bayesian method.(2)Rock-support interaction is an essential component of support design for the tunnel. The convergence-confinement method was adopted to characterize the rock-support interactive mechanism. Reliability-based support design can be obtained through combining reliability-based design and the convergence-confinement method.(3)The developed method provides a practical and rational framework for considering the uncertainty for the support design of the tunnel. It only considered the uncertainty of rock strength parameters. Other uncertainties, such as support structure, in situ stress, GSI, etc., can be considered by including these uncertainties in the reliability model.(4)The developed method provided a general framework to deal with the uncertainty. Hoek–Brown failure criterion, FORM, and ABC can be replaced by other failure criteria, reliability method, and optimal method, respectively.
Data Availability
The original data involved in the article already exists in the figures in the article. All the figures are produced by the original data and can be edited.
Conflicts of Interest
The authors declare that they have no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was funded by the National Natural Science Foundation of China (Grant no. 51874119). Grateful appreciation is expressed for the support.