Next Article in Journal
Spatial Differences in Outdoor Thermal Comfort during the Transition Season in Cold Regions of China
Previous Article in Journal
Effect of Municipal Solid Waste Incineration Ash on Microstructure and Hydration Mechanism of Geopolymer Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Data-Driven Rock Strength Parameter Identification Using Artificial Bee Colony Algorithm

1
School of Civil and Architectural Engineering, Shandong University of Technology, Zibo 255000, China
2
State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
*
Author to whom correspondence should be addressed.
Buildings 2022, 12(6), 725; https://doi.org/10.3390/buildings12060725
Submission received: 29 April 2022 / Revised: 18 May 2022 / Accepted: 24 May 2022 / Published: 26 May 2022
(This article belongs to the Section Building Materials, and Repair & Renovation)

Abstract

:
Rock strength parameters are essential to understanding the rock failure mechanism and safely constructing rock excavation. It is a challenging problem for determining the rock failure criterion and its parameters due to the complexity of rock media. This study adopts an artificial bee colony (ABC) algorithm to determine the Hoek-Brown failure criterion, widely used in rock engineering practice, based on experimental data. The ABC-based approach is presented in detail and applied to a collection of experimental data collected from the literature. The ABC-based approach successfully determines the Hoek-Brown failure criterion, and the determined failure envelope is in excellent agreement with the measured curve. The maximum relative error obtained by ABC is only 2.15% and is far less than the 12.24% obtained by the traditional method. Then, the developed approach is applied to the Goupitan Hydropower Station, China, and determines the rheological parameters of soft rock based on the Burgers model. The deformation of an experiment located in the Goupitan Hydropower Station is evaluated based on obtained parameters by the developed approach. The predicted deformation matches the monitored displacement in the field. The obtained parameters of the failure criterion characterize the mechanical behavior of rock mass well. Thus, the method used provides a reliable and robust approach to determining the mechanical parameters of the failure criterion.

1. Introduction

Stability analysis is one of the most challenging issues in rock excavation activities, such as tunnelling, sloping, coal mining, etc. The rock failure criterion is critical to estimating the stability of rock structures. Various failure criteria have been developed for rock engineering in the last decades [1,2,3,4,5,6]. The Hoek-Brown failure criterion, derived for brittle and jointed rock masses by Hoek and Brown, has been widely used in practical rock engineering for the last several decades [7,8,9,10]. In order to estimate and predict the accuracy of the failure behavior for a rock mass, it is essential to determine the Hoek-Brown failure criterion, which is difficult to measure experimentally in laboratory and field tests. The rational determination of parameters for the model have drawn increased attention from researchers and engineers in rock engineering.
In order to understand and characterize the failure mechanism and deformation behavior of a rock mass, the strength failure criterion is necessary for stability analysis and construction in rock engineering. The Hoek-Brown failure criterion is a relatively simple function with a small number of parameters; it defines rock strength as a nonlinear function of confining pressure. The Hoek-Brown failure criterion provides a standardized method for determining rock mass strength parameters based on intact laboratory test data [11]. Mathematical and statistical methods are essential to approximate the failure criterion and its parameters [12,13]. Simple linear regression was used to determine the rock failure criterion, illustrated by three examples [14]. The simplex reflection method was applied to fit the rock failure criterion based on the laboratory test data [15]. A bonded block numerical model was developed to numerically quantify the strength of fractured rock mass [16]. A GSI-softening model was proposed to capture the strength behavior of thermally damaged rock based on the Hoek-Brown failure criterion [17]. With the development of computational technology, nonlinear regression was also used to determine the rock failure criterion [18,19]. The multiple regression method was developed to estimate the rock failure criterion [20]. Using a single empirical model is not enough to describe the failure envelope in tension and high compression zones. A bi-segmental Hoek-Brown failure criterion was developed to overcome the above problem by considering the failure mechanism in different zones [21]. A mathematical formula was developed to characterize the fatigue and strength of concrete using symbolic regression [22]. Recently, researchers applied Bayes theory to determine the rock failure criterion [23,24]. Bayesian inference was used to determine the rock strength anisotropy and evaluate the uncertainty of the Hoek–Brown failure criterion [25]. With the development of artificial intelligence, the soft computing method provides an optimal method for determining the parameters of the model [26,27]. The genetic algorithm was used to determine the Mohr-Coulomb constitutive model [28]. Symbolic regression-based data were developed to determine the fatigue equation [22]. Zhao et al. adopted particle swarm optimization (PSO) to determine the geomechanical parameters and perform reliability analysis for geotechnical engineering [29,30]. A fuzzy inference system, artificial neural network, and adaptive neuro-fuzzy inference system were utilized to predict the Brazilian tensile strength of rock samples [17]. The artificial Bee Colony (ABC) algorithm, which simulates the intelligent foraging behavior of honey bee swarms, was proposed by Karaboga for optimizing numerical problems [31]. Sonmez and Kang applied it to structural engineering [32,33]. Zhao et al. applied ABC to the back analysis in geotechnical engineering [34]. Villegas et al. applied ABC to optimize the values for the variables of a proportional integral controller for observing the behavior of the controller [35]. Caraveo et al. developed a modification of a bio-inspired algorithm based on ABC for optimizing fuzzy controllers [36]. An improved artificial bee colony algorithm was proposed to address the optimization of the injection scheme of steam flooding [37]. A hybrid model was proposed for predicting the cemented paste backfill strength based on an adaptive neuro-fuzzy inference system and artificial bee colony [38]. Yavuz et al. developed the artificial bee colony algorithm with distant savants for solving constrained optimization problems [39]. Uncertainty quantification was utilized to capture the strength of the reinforcement rock mass [40]. The maximum entropy probability density function was adopted to obtain the uncertainty of the shear strength of the intact rock in the Jinping II project [41]. In this study, ABC is adopted to determine rock strength parameters based on testing data for rock engineering.
It is not easy to determine the strength parameters due to the complexity and unclear failure mechanism of the rock mass. ABC is a straightforward, robust, and population-based stochastic optimization algorithm that does not depend on the specific problem and is suitable for the black-box problem. The significant advantage of ABC is it is easy to implement without any gradient information. Moreover, ABC offers a good balance between fast convergence and exploratory search, allowing for escaping from local minima and aiming for global instead. Meanwhile, ABC could provide a data-driven model capable of correctly determining the strength parameters of the rock mass based on laboratory tests. The remainder of this paper is organized as follows. Firstly, the Hoek-Brown failure criterion and the Burgers model of rock mass are introduced in detail. Secondly, the algorithm of ABC is described, and the procedure of determining the rock failure criterion and its parameters is presented briefly. Then, the proposed method is verified and illustrated using the test data of the uniaxial compression test, triaxial compression test, and tensile test, and is applied to the soft rock mass in a practical experimental tunnel. Finally, some conclusions are presented.

2. Rock Failure Criteria

2.1. Hoek-Brown Failure Criterion

Rock failure originates from micro-cracks or flaws in the rock mass. Griffith developed the tensile failure and deepened it to obtain a nonlinear compressive failure envelope for brittle materials [42]. On the basis of the above nonlinear Griffith failure criterion, Hoek and Brown proposed the generalized Hoek-Brown criterion for rock mass strength [7]. The generalized Hoek-Brown failure criterion could be expressed as follows:
σ 1 = σ 3 + σ c i ( m b σ 3 σ c i + s ) a
where σ1 and σ3 denote the major and minor principal stresses, respectively; σci denotes the unconfined compressive strength; mb, s, and a denote the rock mass material constants, determined as follows:
m b = m i e G S I 100 28 14 D
s = e G S I 100 9 3 D
a = 1 2 + 1 6 ( e G S I 15 e 20 3 )
where GSI denotes the geological strength index, which presents the failure criterion based on the engineering geology observations in the field; D denotes a factor that considers the degree of the rock mass disturbance by blast damage and stress relaxation. For the intact rock, the material constants are adopted by mi, s = 1, and a = 0.5. The Hoek-Brown failure criterion was employed to capture the results of a wide range of triaxial tests on intact rock samples.
σ 1 = σ 3 + σ c i m i σ 3 σ c i + 1
In order to characterize the Hoek-Brown failure criterion for intact rock mass, mi and σci are determined based on the experimental data. This study utilizes ABC to characterize the Hoek-Brown failure criterion and determine its material constants based on laboratory tests.

2.2. Burgers Model

The rheological properties of rock mass characterize the deformation behavior of the surrounding rock mass over time during rock excavation. The rock failure criterion of the combination model is commonly adopted in rheological numerical calculation and analysis. Its advantage is that the appropriate combination of idealized basic elements can be utilized to capture the complex rheological properties of the actual rock mass, and the mechanical concept and the physical meaning are also simple and clear. The Burgers model is generally adopted to capture the mechanical and deformation behavior of soft rock [43]. It consists of the Kelvin model and the Maxwell model (Figure 1). The Burgers model of rock failure can be expressed as follows [44,45,46]:
σ ¨ + ( E 1 η 2 + E 1 η 1 + E 2 η 1 ) σ ˙ + E 1 E 2 η 1 η 2 σ = E 1 ε ¨ + E 1 E 2 η 2 ε ˙
where σ and ε denote the stress and strain of the Burgers model, E1 and η1 denote the elastic modulus and viscosity for the Kelvin model, and E2 and η2 denote the elastic modulus and viscosity for the Maxwell model. Under constant load σ0, the Burgers equation can be changed into the following equation:
ε = σ 0 { 1 E 1 + t η 1 + 1 E 2 [ 1 exp ( E 2 t η 2 ) ] }
The Burgers model is utilized in a soft rock mass with instantaneous deformation, deceleration creep, constant velocity creep, relaxation, and elastic after-effect.

3. The Artificial Bee Colony Algorithm

The artificial bee colony (ABC) algorithm was first developed for unconstrained optimization problems, where it showed good performance [47]. ABC was inspired by the intelligent foraging behavior of bee swarms. There are various tasks done by specialized bee individuals in the actual bee colony. In the ABC algorithm, the colony of artificial bees includes three groups of bees, i.e., employed bees, onlooker bees, and scout bees. Employed bees are used to exploit the nectar sources previously explored and then transfer information to the other waiting bees in the hive about the quality of the food source they seek. Onlooker bees wait in the hive and build a food source to exploit based on the information gained from the employed bees. Scout bees seek an environment in order to obtain a new food source. Once the ABC algorithm is initialized, it needs a cycle of three phases: employed bee phase, onlooker bee phase, and scout bee phase.

3.1. Initialization Phase

Firstly, the ABC algorithm generates the initial population of SN solutions/bees randomly and calculates the fitness function of each solution/bee.
x ( i , j ) = x min j + r a n d ( 0 , 1 ) ( x max j x min j )
where SN is an integer; x(i,j) denotes the candidate solution of problem; i = 1, 2,..., SN/2 and SN/2 denotes the size of population; j = 1,2, …, Dim and Dim denotes the dimension number of each solution; rand(0,1) denotes a random number between [0, 1]; xjmin and xjmax denote the upper and lower bound of each solution.

3.2. Employed Bee Phase

While initialization is finished, employed bees generate the new solution based on the ABC algorithm (Figure 2). A candidate solution can be obtained using the memory of bees, as follows:
v ( i , j ) = x ( i , j ) + φ i j ( x ( i , j ) x ( k , j ) )
where k is different from i and denotes the indexes randomly chosen from {1, 2,..., SN/2}, j also denotes randomly chosen indexes from {1, 2,..., Dim}, φij denotes a random number in [−1, 1], characterizes the generation of neighbor food sources around x(i,j), and denotes the comparison of two food positions seen by a bee.

3.3. Onlooker Bee Phase

Onlooker bees determine a solution based on the fitness function gained from employed bees, generate the solution to be abandoned, and allocate its employed bees as scout bees. The probability of being selected for each fitness function pi can be computed as follows:
p i = f i t n e s s i n = 1 S N / n f i t n e s s n
where fitnessi denotes the fitness value of the solution.

3.4. Scout Bee Phase

In the ABC algorithm, scout bees randomly seek a new solution in the predetermined searching ranges. No further improved solution is assumed to be abandoned by onlookers in a predetermined number of cycles. The abandoned solution x(i,j) is updated by a new solution x’(i,j), generated by the scout bees using Equation (9). Each candidate solution v(i,j) generated by x(i,j) can be determined using the comparison between x(i,j) and its old solution. The old solution will be updated by the new solution when it equals or is better than the old solution. Otherwise, the old solution is unchanged in the memory.

3.5. Procedure of ABC

ABC mainly depends on three controlling parameters: the number of food sources equals the number of employed or onlooker bees (SN), the value of limit, and the maximum cycle number (MCN). The fitness function is essential to ABC algorithm, which depends on the specific problem. In this study, mean absolute error and mean relative error are utilized to determine the fitness function. The detailed formation is presented in Section 4.1. The procedure of ABC can be described in detail as follows:
Step1:
Set the value of control parameters SN, MCN, and “limit” of the ABC algorithm.
Step2:
Initialize the population x(i,j) using Equation (8) and calculate the fitness value of each solution.
Step3:
For each employed bee, generate a new solution v(i,j) using Equation (9) and calculate its fitness.
Step4:
Determine the probability pi for the solution x(i,j) using Equation (10).
Step5:
For each onlooker bee, determine a solution x(i,j) based on pi, generate a new solution v(i,j), and calculate the fitness.
Step6:
If there is an abandoned solution for the scout, it is in place of a new solution that will be randomly generated using Equation (8).
Step7:
Record the best solution.
Step8:
Repeat Step3 to Step7 until reaching the maximum cycle.

4. Determination of Hoek-Brown Failure Criterion Based on Artificial Bee Colony

Based on the laboratory test, the ABC-based approach is developed to characterize the Hoek-Brown failure criterion. The ABC is employed to seek the material constant mi and rock compressive strength σci based on the error function, which is the difference between the predicted and actual values.

4.1. Fitness Function

In order to use the ABC algorithm to characterize the Hoek-Brown failure criterion, as in any conventional optimal method, it is necessary to construct the error function/fitness function for ABC. In this study, the absolute error and relative error that Douglas recommended are defined as follows [20]:
f i t n e s s = {   σ 1 σ 1 n σ 1 > 3 σ 3   σ 3 σ 3 n σ 1 3 σ 3
f i t n e s s = {   σ 1 σ 1 σ 1 n σ 1 > 3 σ 3   σ 3 σ 3 σ 3 n σ 1 3 σ 3
where σ1′ and σ3′ denote the predicted stresses using the Hoek-Brown failure criterion, σ1 and σ3 denote the measured stress values from the experiment, and n denotes the number of experimental data. In this study, the root mean squared error is also utilized to illustrate the performance of ABC.
f i t n e s s = {   ( σ 1 σ 1 ) 2 n σ 1 > 3 σ 3   ( σ 3 σ 3 ) 2 n σ 1 3 σ 3
f i t n e s s = {   ( σ 1 σ 1 σ 1 ) 2 n σ 1 > 3 σ 3   ( σ 3 σ 3 σ 3 ) 2 n σ 1 3 σ 3

4.2. Procedure of Determination of Hoek-Brown Failure Criterion

Once the rock laboratory experiments are finished, the fitness value can be determined based on testing data and the error function in the above section. ABC is utilized to seek the material parameters of the Hoek-Brown failure criterion based on the error function. The analysis proceeds through these steps as follows (Figure 3):
  • Step 1: Collect the information on the rock and determine the parameters of ABC.
  • Step 2: Determine the test scheme based on rock mass property.
  • Step 3: Implement the rock test according to the test standard.
  • Step 4: Generate the test data.
  • Step 5: Calculate the fitness value based on the fitness function using the testing data above.
  • Step 6: Use the ABC algorithm to seek material constants of the Hoek-Brown failure criterion.
  • Step 7: Characterize the Hoek-Brown failure criterion.

4.3. Verification

The Hoek-Brown failure criterion is utilized to capture the rock properties based on experimental data. For a good estimation, the more experimental data, the more accurate the estimation of the rock properties will be. Sari collected 266 experimental data points, including the uniaxial compression test, Brazilian tensile test, and triaxial compression test [48]. The data points are very widely scattered (Figure 4). Sari adopted the linear regression method and non-linear Levenberg-Marquardt fitting algorithm to determine the Hoek-Brown failure criterion based on the above experimental data. This study utilizes ABC to characterize the Hoek-Brown failure criterion (Figure 5). Based on the Hoek-Brown criterion (Equation (1)) and fitness function (Equations (11)–(14)), material constant σci and mi were determined using ABC. The major principal stress σ1 and minor principal stress σ3 are the input parameters from 266 experimental data for the ABC model based on the fitness function. The results are compared with the different methods. The performance of ABC and estimation of rock properties are illustrated using the different parameters of ABC, searching range, and fitness/error function.
Hoek-Brown failure criterion σc and mi parameters are listed in Table 1 based on the 266 experimental data points using the different methods. The obtained strength parameters are closer to the measured parameters than the modified least square method (MLS). Based on the MLS, the maximum relative error is 12.24%, but the maximum relative error is only 2.15% based on ABC. Meanwhile, the obtained parameters are almost identical to the data, including and excluding tensile experiments using ABC. This shows that ABC can capture the strength property of rock mass better than MLS. Figure 6 shows the comparison of the obtained parameters using the different methods. The above investigation shows that the proposed method has an optimal global performance. The envelopes measured and obtained by ABC and MLS is shown in Figure 7. The envelope obtained using ABC is closer to the measured envelope than that obtained using MLS when using data included in the tensile experiment. The proposed method is obviously superior to MLS for the envelope excluding tensile experiment data. With the increase of the confining pressure, the difference between the measured parameters and those predicted using ABC will increase. However, the difference is very insignificant. This further proves that the proposed method reasonably determines the Hoek-Brown failure criterion using the experimental data points.
The convergence performance of the proposed method is shown in Figure 8. Figure 9 shows the variation process of parameters of the failure criterion. It shows that the proposed method has good convergence and can quickly determine the Hoek-Brown failure criterion based on the experimental data points. The fitness of the population is shown in Figure 10. The results are similar to the above results. This shows that ABC has a fast convergence performance. Based on the test data, the ABC-based approach can accurately determine the Hoek-Brown failure criterion.
The optimal global performance is an essential index to the optimal method. In order to verify the ABC algorithm, a different searching range is used to determine the parameters of the failure criterion. Table 2 lists the results in different searching ranges. The ABC-based approach can determine well the parameters of the failure criterion and find the minimum fitness. The proposed method has a good global optimal performance. The ABC-based approach can avoid the optimal local solution and determine the best solution for the Hoek-Brown failure criterion. Figure 11 shows the convergence in different ranges. The ABC-based approach can determine the solution at about Cycle 10 in a smaller range, and it can determine the solution at about Cycle 20 in a more extensive range. This shows that the proposed method is suitable and can quickly determine the solution in a smaller range.
The objective function drives the optimal method to search for the solution and is essential to the optimal method. Sari investigated the fitness function in multiple regression and found that the fitness function is essential for determining the Hoek-Brown failure criterion [48]. In this study, various fitness functions were used to verify the performance of the ABC-based approach. Figure 12 shows the Hoek-Brown failure envelope using a different fitness function. The fitness function has little influence on the Hoek-Brown failure criterion when using the ABC-based approach. It is important to determine the failure criterion because of the difficulty of determining the fitness function. This shows that the ABC-based approach is a reliable and robust method for determining the Hoek-Brown failure criterion.

5. Application

The Goupitan Hydropower Station is located in Wujiang, Guizhou Province, China. It is a landmark electricity transmission project ranging from the west to the east of the province [49]. In order to characterize the rheological property of the rock mass, rock samples of gray-green claystone and mauve claystone were analyzed using the uniaxial compression rheological test in the laboratory. The samples were taken as blocks on site and processed into ϕ50 mm × 100 mm. The test was implemented on an LYJ rheometer. The loading method was step-by-step loading, and the load direction was perpendicular to the rock layer. After loading the first stage load, the instantaneous displacement was measured immediately; then, the creep displacement was measured at 5 min, 10 min, 15 min, 30 min, 1 h, 2 h, 4 h, 8 h, 12 h, 16 h, and 24 h (cumulative time). After 24 h, the creep displacement was measured twice a day until the creep displacement was stable (24 h displacement was no more than 0.001 mm). Then, the next shear load was applied, and the above measurements were repeated until the specimen was rheologically damaged. The strain increment was determined based on the relationship between the displacement and time at the various stress levels to the step-by-step loading. The curve between time and the increment of strain is shown in Figure 13 for different stress levels.
Once the testing data were obtained, ABC could be used to determine the rheological parameters of the Burgers model. The size of the population and the maximum cycle of ABC were 50 and 100, respectively. The range of rheological parameters was determined based on the field measurements and laboratory tests (Table 3). The rheological parameters of grey-green claystone are shown in Figure 14 for the different stress levels. We can see from Figure 14 that the rock rheological parameters depend on the stress level (loading) and are not constant under different stress conditions. It is difficult to reveal the complex and uncertain relationship between rheological parameters and stress levels (loading). Generally, the average value of different stress levels is regarded as the identified rheological parameters in rock engineering. In this study, the rheological parameters were determined based on the stress conditions of the Goupitan experiment tunnel. The width, height, and buried depth of the tunnel were 2 m, 2 m, and 70 m, respectively. The experimental tunnel was excavated through rock layer S 2 h 1 1 ; the S 2 h 1 2 stratum lies about 3 m below the tunnel, and layer S 1 q 1   lies 30 m above it (Figure 15). According to the buried depth, the testing data at 1.48 MPa and 2.71 MPa was adopted to determine the rheological parameters of grey-green claystone and mauve claystone, respectively. The values of the rheological parameters are listed in Table 4. Figure 16 shows the comparison between testing data and predicted data based on the obtained rheological parameters using the Burgers model. The rheological parameters are determined based on the testing data of stress level 1.48 MPa. The predicted strain is in excellent agreement with the testing data (Figure 16b). This shows that ABC can better characterize the mechanical properties. For the other stress level, there is a higher error. The predicted curve is better closer to the 1.48 MPa stress level (Figure 16a,c,d) than further away (Figure 16e–g). There is a larger discrepancy between the test data and predicted values in Figure 16f because the stress level is 4.77 MPa, far removed from 1.48 MPa. This shows that the rock rheological parameters should be determined based on the practical in situ stress conditions. A comprehensive comparison of the predicted curve based on the full test data is shown in Figure 16h. It shows that the rheological parameters should be determined based on the field stress conditions. Figure 17 shows the results of mauve claystone at a 2.71 MPa stress level. This proves that the parameters obtained by ABC characterize well the deformation and mechanical properties of rock.
Figure 18 shows the changing process of population fitness with the increase of the cycle. The fitness of the population is converged to the constant with the cycle. Figure 19 shows the convergence of the ABC algorithm. This shows that the developed method has a good global optimal and convergence performance. This proves again that the developed method has a good performance and can quickly determine the solution in a more extensive range.
In the field, an experimental tunnel was excavated to observe the rheological behavior of the claystone. Some monitoring sites were set up to investigate the deformation of the surrounding rock during excavation. The experimental tunnel was analyzed based on the determined rheological parameters in this study. The width, height, and buried depth of the tunnel were 2 m, 2 m, and 70 m, respectively. Fast Lagrangian analysis of continua (FLAC) software was utilized to compute the deformations of the surrounding rock mass. The deformation of the surrounding rock mass was evaluated based on the rheological parameters of the Burgers model. The displacements of the monitored borehole were then evaluated by a numerical model that used the determined rheological parameters. These are compared in Figure 20, in which it is evident that the predicted displacement meets the laws of monitored displacement with increasing time. This proves that the testing data are helpful for characterizing the mechanical behavior of the rock mass and can be utilized for capturing the time-dependent properties of the surrounding rock mass. The geotechnical parameters obtained in this way can be helpful for design, stability analysis, construction, and safety management in rock excavation.

6. Conclusions

The rock failure criterion is critical to rock excavation. However, the determination of rock strength is a challenging issue because of the unclear mechanism of rock failure and the complexity of geology conditions. This study develops an ABC-based approach to determine the rock strength parameters by combing them with the experimental data.
(1)
In this study, ABC was utilized to determine the strength parameters, which characterize the rock failure mechanism and deformation behavior. Once the strength parameters were determined, the corresponding failure criterion could be used to evaluate the stability of the rock mass and determine the supporting pattern of the surrounding rock mass in practical engineering.
(2)
Laboratory testing is a common way to determine rock mass strength. This study developed an ABC-based approach to characterize the rock strength properties based on the Hoek-Brown and Burges models. Test data, hiding the rock failure mechanism, were utilized to capture the rock failure criterion based on the mathematical tool. The determined Hoek-Brown failure envelope was in excellent agreement with the experimental curve. The ABC-based approach provides a promising and scientific way to determine the rock failure criterion.
(3)
The fitness function is an essential component of the developed approach. The ABC-based approach can determine the optimal material constants using different fitness functions based on experimental data and avoids the limitations and disadvantages of the traditional methods. Meanwhile, the strength parameters obtained by the developed approach characterize well the deformation and strength properties of the rock mass.
(4)
It is challenging to understand and characterize the failure mechanism and deformation behavior of a rock mass. Thus, it is not easy to determine the rock strength criterion and its coefficient. The ABC-based approach has a good performance for global optimization; it can avoid the optimal local solution and provides an alternative tool to address it, which is illustrated by its performance. The developed approach focuses on the Hoek-Brown failure criterion and the Burger model for rock mass. It is worth noting that the proposed ABC approach can be employed to determine the strength parameters of other failure criteria in rock mechanics and engineering. Further study is necessary for various strength criteria in the future.
(5)
In this study, ABC was adopted to determine the strength parameters of rock mass based on laboratory data. ABC has been proven to have a strong global searching capability, which can significantly increase the efficiency of the strength parameter determination process. However, the efficiency depends on the number of laboratory data. With the increase in laboratory data, the computation time will increase. The developed method will be further studied by combining it with a new seeking strategy in the future.

Author Contributions

Conceptualization, H.Z.; methodology, H.Z. and B.C.; software, H.Z.; validation, H.Z. and M.W.; formal analysis, H.Z.; investigation, H.Z. and B.C.; data curation, H.Z. and M.W.; writing—original draft preparation, H.Z.; writing—review and editing, B.C. and M.W.; visualization, H.Z.; project administration, B.C.; funding acquisition, H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Open Research Fund of the State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, grant number Z020006.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Some or all data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors sincerely acknowledge the anonymous reviewers for their insights and comments to further improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hoek, E.; Brown, E.T. Underground Excavations in Rock; London Institution of Mining and Metallurgy: London, UK, 1980. [Google Scholar]
  2. Yang, Y.; Xu, D.; Liu, F.; Zheng, H. Modeling the entire progressive failure process of rock slopes using a strength-based criterion. Comput. Geotech. 2020, 126, 103726. [Google Scholar] [CrossRef]
  3. Sheorey, P.; Biswas, A.; Choubey, V. An empirical failure criterion for rocks and jointed rock masses. Eng. Geol. 1989, 26, 141–159. [Google Scholar] [CrossRef]
  4. Labuz, J.F.; Zeng, F.; Makhnenko, R.; Li, Y. Brittle failure of rock: A review and general linear criterion. J. Struct. Geol. 2018, 112, 7–28. [Google Scholar] [CrossRef]
  5. Ramamurthy, T. Shear strength response of some geological materials in triaxial compression. Int. J. Rock Mech. Min. Sci. 2001, 38, 683–697. [Google Scholar] [CrossRef]
  6. Zhao, L.; Yu, C.; Cheng, X.; Zuo, S.; Jiao, K. A method for seismic stability analysis of jointed rock slopes using Barton-Bandis failure criterion. Int. J. Rock Mech. Min. Sci. 2020, 136, 104487. [Google Scholar] [CrossRef]
  7. Hoek, E.; Brown, E.T. The Hoek-Brown failure criterion and GSI-2018 edition. J. Rock Mech. Geotech. Eng. 2019, 11, 445–463. [Google Scholar] [CrossRef]
  8. Park, D.; Michalowski, R.L. Three-dimensional stability assessment of slopes in intact rock governed by the Hoek-Brown failure criterion. Int. J. Rock Mech. Min. Sci. 2021, 137, 104522. [Google Scholar] [CrossRef]
  9. Lin, H.; Singh, S.; Oh, J.; Canbulat, I.; Kang, W.; Hebblewhite, B.; Stacey, T. A combined approach for estimating horizontal principal stress magnitudes from borehole breakout data via artificial neural network and rock failure criterion. Int. J. Rock Mech. Min. Sci. 2020, 136, 104539. [Google Scholar] [CrossRef]
  10. He, M.; Zhang, Z.; Zheng, J.; Chen, F.; Li, N. A New Perspective on the Constant mi of the Hoek–Brown Failure Criterion and a New Model for Determining the Residual Strength of Rock. Rock Mech. Rock Eng. 2020, 53, 3953–3967. [Google Scholar] [CrossRef]
  11. Hoek, E.; Brown, E.T. Practical estimates of rock mass strength. Int. J. Rock Mech. Min. Sci. 1997, 34, 1165–1186. [Google Scholar] [CrossRef]
  12. Li, Z.; Liu, J.; Xu, R.; Liu, H.; Shi, W. Study of grouting effectiveness based on shear strength evaluation with experimental and numerical approaches. Acta Geotech. 2021, 16, 3991–4005. [Google Scholar] [CrossRef]
  13. Su, G.-S.; Gan, W.; Zhai, S.-B.; Zhao, G.-F. Acoustic emission precursors of static and dynamic instability for coarse-grained hard rock. J. Cent. South Univ. 2020, 27, 2883–2898. [Google Scholar] [CrossRef]
  14. Hoek, E. Strength of jointed rock masses. Geotechniques 1983, 33, 187–223. [Google Scholar] [CrossRef] [Green Version]
  15. Shah, S.; Hoek, E. Simplex reflection analysis of laboratory strength data. Can. Geotech. J. 1992, 29, 278–287. [Google Scholar] [CrossRef]
  16. Zhang, Y.; Liu, X.; Guo, P.; Jia, P.; Yang, T.; Ren, F. Quantifying strength and permeability of fractured rock mass using 3D bonded block numerical model. J. Pet. Sci. Eng. 2021, 208, 109592. [Google Scholar] [CrossRef]
  17. Peng, J.; Tang, Z.C.; Hou, D. A GSI-softening model for characterizing strength behavior of thermally-damaged rock. Eng. Geol. 2021, 292, 106251. [Google Scholar] [CrossRef]
  18. Hoek, E.; Carranza-Torres, C.; Corkurn, B. Hoek-Brown failure criterion-2002 edition. In Proceeding of the 5th North American Rock Mechanics Symposium and the 17th Tunnelling Association of Canada Conference (NARMS-TAC 2002); University of Toronto Press: Toronto, ON, Canada, 2002; Volume 1. [Google Scholar]
  19. Langford, J.C.; Diederichs, M. Quantifying uncertainty in Hoek–Brown intact strength envelopes. Int. J. Rock Mech. Min. Sci. 2015, 74, 91–102. [Google Scholar] [CrossRef]
  20. Douglas, K.J. The Shear Strength of Rock Masses. Ph.D. Thesis, University of New South Wales, Sydney, Australia, 2002. [Google Scholar]
  21. Cai, M. Practical estimates of tensile strength and Hoek-Brown strength parameter mi of brittle rocks. Rock Mech. Rock Eng. 2010, 43, 167–184. [Google Scholar] [CrossRef]
  22. Ren, J.; Zhang, L.; Zhao, H.; Zhao, Z.; Wang, S. Determination of the fatigue equation for the cement-stabilized cold recycled mixtures with road construction waste materials based on data-driven. Int. J. Fatigue 2022, 158. [Google Scholar] [CrossRef]
  23. Contreras, L.F.; Brown, E.T.; Ruest, M. Bayesian data analysis to quantify the uncertainty of intact rock strength. J. Rock Mech. Geotech. Eng. 2018, 10, 11–31. [Google Scholar] [CrossRef]
  24. Bozorgzadeh, N.; Escobar, M.D.; Harrison, J.P. Comprehensive statistical analysis of intact rock strength for reliability-based design. Int. J. Rock Mech. Min. Sci. 2018, 106, 374–387. [Google Scholar] [CrossRef]
  25. Gomes, G.J.; Forero, J.H.; Vargas, E.A.; Vrugt, J.A. Bayesian inference of rock strength anisotropy: Uncertainty analysis of the Hoek–Brown failure criterion. Int. J. Rock Mech. Min. Sci. 2021, 148, 104952. [Google Scholar] [CrossRef]
  26. Pichler, B.; Lackner, R.; Mang, H.A. Back analysis of model parameters in geotechnical engineering by means of soft computing. Int. J. Numer. Methods Eng. 2003, 57, 1943–1978. [Google Scholar] [CrossRef]
  27. Ren, J.; Zhao, H.; Zhang, L.; Zhao, Z.; Xu, Y.; Cheng, Y.; Wang, M.; Chen, J.; Wang, J. Design optimization of cement grouting material based on adaptive boosting algorithm and simplicial homology global optimization. J. Build. Eng. 2022, 49, 104049. [Google Scholar] [CrossRef]
  28. Levasseur, S.; Malecot, Y.; Boulon, M.; Flavigny, E. Soil parameter identification using a genetic algorithm. Int. J. Numer. Anal. Methods Geomech. 2008, 32, 189–213. [Google Scholar] [CrossRef]
  29. Zhao, H.B.; Yin, S.D. Geomechanical parameters identification by particle swarm optimization and support vector ma-chine. Appl. Math. Model. 2009, 33, 3997–4012. [Google Scholar] [CrossRef]
  30. Zhao, H.; Ru, Z.; Chang, X.; Li, S. Reliability Analysis Using Chaotic Particle Swarm Optimization. Qual. Reliab. Eng. Int. 2014, 31, 1537–1552. [Google Scholar] [CrossRef]
  31. Karaboga, D. An Idea Based on Honey Bee Swarm for Numerical Optimization; Technical Report-TR06; Erciyes University, Engineering Faculty, Computer Engineering Department: Kayseri, Turkey, 2005. [Google Scholar]
  32. Sonmez, M. Artificial Bee Colony algorithm for optimization of truss structures. Appl. Soft Comput. 2011, 11, 2406–2418. [Google Scholar] [CrossRef]
  33. Kang, F.; Li, J.; Xu, Q. Structural inverse analysis by hybrid simplex artificial bee colony algorithms. Comput. Struct. 2009, 87, 861–870. [Google Scholar] [CrossRef]
  34. Zhao, H.B.; Ru, Z.; Li, S. Coupling relevance vector machine and response surface for geomechanics parameters identification. Geomech. Eng. 2018, 15, 1207–1217. [Google Scholar]
  35. Villegas, J.M.; Caraveo, C.; Mejía, D.A.; Rodríguez, J.L.; Vega, Y.; Cervantes, L.; Medina-Santiago, A. Intelligent Search of Values for a Controller Using the Artificial Bee Colony Algorithm to Control the Velocity of Displacement of a Robot. Algorithms 2021, 14, 273. [Google Scholar] [CrossRef]
  36. Caraveo, C.; Valdez, F.; Castillo, O. Optimization of fuzzy controller design using a new bee colony algorithm with fuzzy dynamic parameter adaptation. Appl. Soft Comput. 2016, 43, 131–142. [Google Scholar] [CrossRef]
  37. Ni, H.-M.; Liu, Y.-J.; Fan, Y.-C. Optimization of injection scheme to maximizing cumulative oil steam ratio based on improved artificial bee colony algorithm. J. Pet. Sci. Eng. 2018, 173, 371–380. [Google Scholar] [CrossRef]
  38. Qi, C.; Ly, H.-B.; Le, L.M.; Yang, X.; Guo, L.; Pham, B.T. Improved strength prediction of cemented paste backfill using a novel model based on adaptive neuro fuzzy inference system and artificial bee colony. Constr. Build. Mater. 2021, 284, 122857. [Google Scholar] [CrossRef]
  39. Yavuz, G.; Durmuş, B.; Aydın, D. Artificial Bee Colony Algorithm with Distant Savants for constrained optimization. Appl. Soft Comput. 2021, 116, 108343. [Google Scholar] [CrossRef]
  40. Zhang, L.; Wang, M.; Zhao, H.; Chang, X. Uncertainty quantification for the mechanical behavior of fully grouted rockbolts subjected to pull-out tests. Comput. Geotech. 2022, 145, 104665. [Google Scholar] [CrossRef]
  41. Deng, J.; Li, S.; Jiang, Q.; Chen, B. Probabilistic analysis of shear strength of intact rock in triaxial compression: A case study of Jinping II project. Tunn. Undergr. Space Technol. 2021, 111, 103833. [Google Scholar] [CrossRef]
  42. Griffith, A.A. The Theory of Rupture. In Proceedings of the 1st International Congress of Applied Mechanics, Delft, The Netherlands, 1 January 1924; pp. 55–63. [Google Scholar]
  43. Goodman, R.E. Introduction to Rock Mechanics, 2nd ed.; Wiley: New York, NY, USA, 1989. [Google Scholar]
  44. Fahimifar, A.; Tehrani, F.M.; Hedayat, A.; Vakilzadeh, A. Analytical solution for the excavation of circular tunnels in a vis-co-elastic Burger’s material under hydrostatic stress field. Tunn. Undergr. Space Technol. 2010, 25, 297–304. [Google Scholar] [CrossRef]
  45. Bai, Q.; Xia, Y.; Liu, X.; Yang, Z. Mining Roadway Support Technology Based on Rheological Analysis in Malin Coalmine. Energy Procedia 2012, 16, 1129–1134. [Google Scholar] [CrossRef] [Green Version]
  46. Du, Y.; Yan, N.; Kortschot, M.T. An experimental study of creep behavior of lightweight natural fiber-reinforced polymer composite/honeycomb core sandwich panels. Compos. Struct. 2013, 106, 160–166. [Google Scholar] [CrossRef]
  47. Karaboga, D.; Basturk, B. On the performance of artificial bee colony (ABC) algorithm. Appl. Soft Comput. 2008, 8, 687–697. [Google Scholar] [CrossRef]
  48. Sari, M. An improved method of fitting experimental data to the Hoek–Brown failure criterion. Eng. Geol. 2012, 127, 27–35. [Google Scholar] [CrossRef]
  49. Chen, B.R. Back Analysis of Rheological Parameters of Rock Mass Using Intelligent Method. Master’s Thesis, Northeastern University, Shenyang, China, December 2003. [Google Scholar]
Figure 1. Burgers model.
Figure 1. Burgers model.
Buildings 12 00725 g001
Figure 2. A simple position update scheme.
Figure 2. A simple position update scheme.
Buildings 12 00725 g002
Figure 3. The ABC-based approach for determination of Hoek-Brown failure criterion.
Figure 3. The ABC-based approach for determination of Hoek-Brown failure criterion.
Buildings 12 00725 g003
Figure 4. The experimental data.
Figure 4. The experimental data.
Buildings 12 00725 g004
Figure 5. The determination of the Hoek-Brown failure criterion using ABC.
Figure 5. The determination of the Hoek-Brown failure criterion using ABC.
Buildings 12 00725 g005
Figure 6. The comparison between the predicted and measured coefficients: (a) σci; (b) mi.
Figure 6. The comparison between the predicted and measured coefficients: (a) σci; (b) mi.
Buildings 12 00725 g006
Figure 7. The Hoek-Brown failure envelope determined using different methods: (a) tensile; (b) no tensile.
Figure 7. The Hoek-Brown failure envelope determined using different methods: (a) tensile; (b) no tensile.
Buildings 12 00725 g007
Figure 8. The convergence of the ABC-based approach.
Figure 8. The convergence of the ABC-based approach.
Buildings 12 00725 g008
Figure 9. The searching process of the parameters of the failure criterion: (a), σc, mi, and fitness; (b) σc, mi, and cycle; (c) σc, cycle, and fitness; (d) cycle, mi, and fitness.
Figure 9. The searching process of the parameters of the failure criterion: (a), σc, mi, and fitness; (b) σc, mi, and cycle; (c) σc, cycle, and fitness; (d) cycle, mi, and fitness.
Buildings 12 00725 g009
Figure 10. The fitness value of population with the cycle: (a) initial cycle, (b) cycle 20, (c) cycle 40, and (d) cycle 60.
Figure 10. The fitness value of population with the cycle: (a) initial cycle, (b) cycle 20, (c) cycle 40, and (d) cycle 60.
Buildings 12 00725 g010
Figure 11. The convergence in different searching ranges.
Figure 11. The convergence in different searching ranges.
Buildings 12 00725 g011
Figure 12. The Hoek-Brown failure envelope using different fitness functions.
Figure 12. The Hoek-Brown failure envelope using different fitness functions.
Buildings 12 00725 g012
Figure 13. The curve of the uniaxial compression rheological test: (a) grey-green claystone, (b) mauve claystone.
Figure 13. The curve of the uniaxial compression rheological test: (a) grey-green claystone, (b) mauve claystone.
Buildings 12 00725 g013
Figure 14. The rheological parameters based on different stress levels using ABC.
Figure 14. The rheological parameters based on different stress levels using ABC.
Buildings 12 00725 g014
Figure 15. The geology map of the Goupitan experimental tunnel.
Figure 15. The geology map of the Goupitan experimental tunnel.
Buildings 12 00725 g015
Figure 16. The curve comparison between test data and values predicted by the obtained rheological parameters using the Burgers model: (a) 0.9 MPa, (b) 1.48 MPa, (c) 2.08 MPa, (d) 2.98 MPa, (e) 3.57 MPa, (f) 4.77 MPa, (g) 5.39 MPa, and (h) predicted curve based on the full test data.
Figure 16. The curve comparison between test data and values predicted by the obtained rheological parameters using the Burgers model: (a) 0.9 MPa, (b) 1.48 MPa, (c) 2.08 MPa, (d) 2.98 MPa, (e) 3.57 MPa, (f) 4.77 MPa, (g) 5.39 MPa, and (h) predicted curve based on the full test data.
Buildings 12 00725 g016
Figure 17. The predicted and testing data at 2.71 MPa stress level for mauve claystone.
Figure 17. The predicted and testing data at 2.71 MPa stress level for mauve claystone.
Buildings 12 00725 g017
Figure 18. The varied process of fitness with the cycle.
Figure 18. The varied process of fitness with the cycle.
Buildings 12 00725 g018
Figure 19. The convergence of ABC algorithm.
Figure 19. The convergence of ABC algorithm.
Buildings 12 00725 g019
Figure 20. The displacement comparison between monitored data and values predicted by the numerical model based on obtained rheological parameters: (a) the monitored point at #4 borehole and (b) the monitored point at #6 borehole.
Figure 20. The displacement comparison between monitored data and values predicted by the numerical model based on obtained rheological parameters: (a) the monitored point at #4 borehole and (b) the monitored point at #6 borehole.
Buildings 12 00725 g020
Table 1. The obtained strengths using different methods.
Table 1. The obtained strengths using different methods.
Material ConstantsInclude TensileNo TensileMeasured
MLSABCMLSABC
σci74.0269.7772.3069.7668.30
Relative error (%)8.372.155.862.14-
mi9.2510.7210.3610.7210.54
Relative error (%)12.241.711.711.73-
Table 2. The parameters of the failure criterion and fitness value in different ranges.
Table 2. The parameters of the failure criterion and fitness value in different ranges.
Rangeσc (MPa)miFitness
[0, 100]67.284211.04399.99552
[0, 200]67.284811.04379.99552
[0, 400]67.282711.04429.99552
[0, 1000]67.283211.04499.99552
Table 3. The ranges of rheological parameters.
Table 3. The ranges of rheological parameters.
Rock TypeE1 (MPa)E2 (MPa)η1 (MPa·h)η2 (MPa·h)
Grey-green claystone5.0 × 102–4.5 × 1031.0 × 102–3.5 × 1033.6 × 104–8.4 × 1062.4 × 102–8.4 × 104
Mauve claystone1.0 × 103–1.5 × 1045.0 × 103–2.0 × 1043.6 × 105–1.08 × 1082.4 × 103–3.6 × 106
Table 4. The value of rheological parameters.
Table 4. The value of rheological parameters.
Rock TypeE1 (MPa)E2 (MPa)η1 (MPa·h)η2 (MPa·h)
Grey-green claystone1169.0985390.94481,041,173.48992500.9780
Mauve claystone2432.531710,800.772136,000,000.000024,000.1511
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, M.; Chen, B.; Zhao, H. Data-Driven Rock Strength Parameter Identification Using Artificial Bee Colony Algorithm. Buildings 2022, 12, 725. https://doi.org/10.3390/buildings12060725

AMA Style

Wang M, Chen B, Zhao H. Data-Driven Rock Strength Parameter Identification Using Artificial Bee Colony Algorithm. Buildings. 2022; 12(6):725. https://doi.org/10.3390/buildings12060725

Chicago/Turabian Style

Wang, Meng, Bingrui Chen, and Hongbo Zhao. 2022. "Data-Driven Rock Strength Parameter Identification Using Artificial Bee Colony Algorithm" Buildings 12, no. 6: 725. https://doi.org/10.3390/buildings12060725

APA Style

Wang, M., Chen, B., & Zhao, H. (2022). Data-Driven Rock Strength Parameter Identification Using Artificial Bee Colony Algorithm. Buildings, 12(6), 725. https://doi.org/10.3390/buildings12060725

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop