Next Article in Journal
General Fractional Calculus in Multi-Dimensional Space: Riesz Form
Next Article in Special Issue
Data-Driven Optimized Artificial Neural Network Technique for Prediction of Flyrock Induced by Boulder Blasting
Previous Article in Journal
Fourier Graph Convolution Network for Time Series Prediction
Previous Article in Special Issue
Refined Design and Optimization of Underground Medium and Long Hole Blasting Parameters—A Case Study of the Gaofeng Mine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predictive Modeling of the Uniaxial Compressive Strength of Rocks Using an Artificial Neural Network Approach

1
School of Mines, China University of Mining and Technology, Xuzhou 221116, China
2
The State Key Laboratory for Geo Mechanics and Deep Underground Engineering, China University of Mining & Technology, Xuzhou 221116, China
3
School of Mines and Civil Engineering, Liupanshui Normal University, Liupanshui 553001, China
4
Guizhou Guineng Investment Co., Ltd., Liupanshui 553001, China
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(7), 1650; https://doi.org/10.3390/math11071650
Submission received: 13 February 2023 / Revised: 15 March 2023 / Accepted: 27 March 2023 / Published: 29 March 2023

Abstract

:
Sedimentary rocks provide information on previous environments on the surface of the Earth. As a result, they are the principal narrators of the former climate, life, and important events on the surface of the Earth. The complexity and cost of direct destructive laboratory tests adversely affect the data scarcity problem, making the development of intelligent indirect methods an integral step in attempts to address the problem faced by rock engineering projects. This study established an artificial neural network (ANN) approach to predict the uniaxial compressive strength (UCS) in MPa of sedimentary rocks using different input parameters; i.e., dry density (ρd) in g/cm3, Brazilian tensile strength (BTS) in MPa, and wet density (ρwet) in g/cm3. The developed ANN models, M1, M2, and M3, were divided as follows: the overall dataset, 70% training dataset and 30% testing dataset, and 60% training dataset and 40% testing dataset, respectively. In addition, multiple linear regression (MLR) was performed for comparison to the proposed ANN models to verify the accuracy of the predicted values. The performance indices were also calculated by estimating the established models. The predictive performance of the M2 ANN model in terms of the coefficient of determination (R2), root mean squared error (RMSE), variance accounts for (VAF), and a20-index was 0.831, 0.27672, 0.92, and 0.80, respectively, in the testing dataset, revealing ideal results, thus it was proposed as the best-fit prediction model for UCS of sedimentary rocks at the Thar coalfield, Pakistan, among the models developed in this study. Moreover, by performing a sensitivity analysis, it was determined that BTS was the most influential parameter in predicting UCS.

1. Introduction

Sedimentary rocks provide information about the previous environment of the Earth’s surface. As such, they are the primary narrators of climate, life, and important events that occurred prior to the Earth’s surface being formed. Uniaxial compressive strength (UCS) is an essential rock strength parameter widely used in the design of rock structures [1,2]. UCS is an integral parameter in rock characterization, tunnel construction, slope stability analysis, construction, bridges, and other rock-related complications [3,4,5,6,7,8]. Direct estimation of UCS based on the principles of ISRM (International Society of Rock Mechanics) and ASTM (American Society for Testing and Materials) is a complex, time-consuming, and expensive procedure. It makes testing infeasible for engineering projects where large amounts of data are needed.
To overcome these shortcomings, this study establishes artificial neural network (ANN) predictive models for the estimation of UCS. Many research scholars have established predictive methods to deal with such complex problems using various statistical methods such as ANN and adaptive neuro-fuzzy interference system (ANFIS) [9,10,11,12,13,14,15,16,17]. Currently, intelligent methods such as ANN, ANFIS, PSO (particle swarm optimization), and GA (genetic algorithm) are frequently applied to solve problems related to rock structure design [2], and these methods are considered to be fast and economical, as well as to have achieved good agreement between the measured and predicted values of rock mechanical properties, i.e., UCS and E (modulus of elasticity in MPa), among others [13]. Torabi-Kaveh employed ANN and multiple regression methods to estimate UCS, and their findings indicated that the ANN method performed better [18]. Yagiz analyzed ANN and multiple regression for predicting UCS of carbonate rocks and found that the ANN method is in good agreement with traditional multiple regression [19]. Ceryan also employed the ANN and regression methods to predict UCS of carbonate rocks and proposed that the ANN results were significantly accurate [20]. Mohamad used a PSO-based ANN method to estimate UCS of soft rocks with input parameters of Brazilian tensile strength (BTS) in MPa, point load index (Is(50)) in MPa, and ultrasonic (Vp) in m/s, and demonstrated the high performance of the proposed model [21]. The ANN method has proved to be a key method among all intelligent methods and is thus mostly used to solve challenging problems that are reliant on laboratory experimental data because of their high efficiency and ability to learn from inputs [22]. Based on the reliable predictions of ANN methods, some researchers have estimated various mechanical properties of rocks by analyzing the correlation among various physical parameters [23,24]. Yin employed an ANN back-propagation algorithm, which has been considered as the best prediction method based on previous studies [25]. Skentou used hybrid ANN models for predicting UCS of granite rocks with optimal results. Similarly [26], Kaloop developed six hybrid ANN models to predict UCS of different rock types. Based on the performance indicators, such as R2 and RMSE [27], the multivariate adaptive regression splines (MARS) revealed ideal results compared with other models developed in the study. Xiang estimated the in situ rock strength from borehole geophysical logs using ANN models [28]. KÖKEN used different soft computing models including ANN for estimating the fracture toughness of rocks [29]. Table 1 shows previous studies using intelligent methods to predict UCS.
This study applied the ANN approach to estimate UCS with different input parameters such as dry density (ρd) in g/cm3, Brazilian tensile strength (BTS) in MPa, and wet density (ρwet) in g/cm3. A total of 78 sedimentary rock samples, i.e., claystone, sandstone, and siltstone, of each type of core rock were selected from Block IX of the Thar coalfield. For the developed ANN models, the dataset is distributed as follows: model 1 (M1) is the overall dataset, model 2 (M2) consists of 70% as the training dataset and 30% as the testing dataset, and model 3 (M3) consists of 60% as the training dataset and 40% as the testing dataset. Similarly, multiple linear regression (MLR) analyses are performed for comparison to the proposed ANN model to check the accuracy of the predicted values. The performance indices are also calculated by estimating the established models. Furthermore, to determine the effect of each variable on the estimated values of UCS, a sensitivity analysis was performed. The complexity and cost of direct destructive laboratory tests adversely affect the data scarcity problem, making the development of intelligent indirect methods an integral step in attempts to address the problem faced by rock engineering projects. In this study, we apply, for the first time, an intelligent prediction method to predict UCS of sedimentary rocks from Block IX of the Thar coalfield. To the best of the authors’ knowledge, there is no such application of intelligent prediction techniques.

2. Materials and Methods

2.1. Building Dataset

In this study, sedimentary rock samples, i.e., claystone, sandstone, and siltstone, were collected from Block IX of the Thar coalfield, Pakistan. Figure 1 represents the geological site of the collected rock samples [38]. Initially, a total of 78 core rock samples of each type were prepared and subdivided into standardized samples according to ISRM and ASTM standards to maintain the same rock core dimensions as well as geological and geotechnical features [39,40]. Next, these rock samples were tested in the laboratory at the Department of Mining Engineering, Mehran University of Engineering and Technology, to determine the physical and mechanical parameters, including ρd in g/cm3, BTS in MPa, ρwet in g/cm3, and UCS in MPa, using a universal testing machine (UTM), as shown in Figure 2a,b. Figure 2a,b represent the deformed rock core specimen for UCS and BTS tests, respectively. Table 2 presents the five heads and five tails of the dataset of physical and mechanical parameters. Table 3 shows the minimum, maximum, average, and standard deviation of parameters of rock samples determined in the laboratory.
Figure 3 represents histogram plots of the original dataset in this study: (a) dry density (g/cm3), (b) BTS (MPa), (c) wet density (g/cm3), and (d) UCS (MPa). Figure 4 presents the pairwise plot of the original dataset of different parameters and UCS under this study. Notably, none of the parameters are well-correlated to the UCS, thus all of the parameters are analyzed for UCS prediction. In addition, Figure 4 represents a moderate positive correlation of BTS with UCS; however, the dry density and wet density show a negative correlation with UCS.

2.2. Methods

The artificial neural network (ANN) approach was employed to predict UCS with three corresponding inputs: ρd (g/cm3), BTS (MPa), and ρwet (g/cm3). Figure 5 demonstrates the flow chart of the predictive modeling process for UCS. Owing to the small number of resources available for collecting samples, the current study used a limited dataset, that is, 78 samples divided for the established models, including M1, M2, and M3, as presented in Table 4. M1 means the model was trained on the overall dataset, M2 means the model was trained on 70% (55 datasets) of the dataset and tested on 30% (23 datasets) of the dataset, and M3 means the model was trained on 60% (47 datasets) of the dataset and tested on 40% (31 datasets) of the dataset. In addition, Taylor diagram representation was used, which explains a brief qualitative depiction of the best fit of the model to standard deviations and correlations. Moreover, cosine amplitude method (CAM)-based sensitivity analysis was carried out in order to estimate the influence of each input variable on output UCS.

2.2.1. Artificial Neural Network

The concept of ANN was originally introduced by Frank Rosenblatt in 1958 [41]. ANN is considered to be the most common and effective soft computing technique based on the function of the human brain’s nervous system [42,43,44,45,46,47]. This technique is mainly used to solve complex rock structure design problems, i.e., mining, civil, geotechnical, geological engineering, and so on. The ANN structure is an essential factor in designing the ultimate prediction model, as the structure affects the learning capability and performance when estimating the network data. The ANN is structured with three layers (i.e., input layer, hidden layer, and output layer) with a number of interrelated units, called neurons, and the method is used to classify the appropriate correlation between the specified input and output parameters [48]. Figure 6 shows the structure of the ANN to estimate UCS in this research. Because of the complexity of the problem, each neuron has sufficient neuron capacity, and each neuron is related to the weight of the next layer [49,50,51]. Equation (1) is used to evaluate the approximate number of neurons in the hidden layer, as the improper selection of the number of neurons in the hidden layer often leads to “under-fitting” and “over-fitting” and must be avoided.
N H 2 N 1 + 1
ANN toolbox in MATLAB package 2018a was used in this study to develop the feed-forward back propagation (FFBP) ANN model with 3-7-1. BP is the most commonly applied powerful learning algorithm in multilayer networks [52,53]. The predictive input parameters, ρd, BTS, and ρwet, were allocated to an input layer composed of three neurons to predict UCS of the output layer. The ANN models, M1, M2, and M3, were trained, tested, and validated. One hundred epochs were used to train the models and the minimum validation error was considered as a stopping point to prevent overfitting. Figure 7 represents the validation curves for the training performance of the ANN models of UCS. Therefore, model M2 demonstrates the best performance curve of UCS, with validation error equal to 0.14562, which is reached at 0 epochs. Figure 8 illustrates the training scatter plots of predicted UCS against measured UCS, as M1 for overall dataset and as M2 and M3 for the training and testing dataset, respectively.

2.2.2. Multiple Linear Regression

SPSS (version 23) was used to conduct a multiple linear regression (MLR) analysis to determine the existence of a linear relationship between the dependent variable and the independent variables. Regression analysis is used to determine the independent variables’ significance in determining the dependent variable’s values [54]. More precisely, the purpose of regression analysis in this study was to compare the performance of the ANN analysis to that of conventional linear regression. This approach has also been used in several recent studies on the application of ANNs and linear regression analysis [55]. The basic linear regression equation (Equation (2)), modified to include our dependent and independent variables, is as follows:
D = α + B 1 T 1 + B 2 T 2 + B 3 T 3 + B n T n + e
where D represents the dependent variable, α represents the regression constant, B represents the regression coefficient, and T represents the value of the independent variable.

2.2.3. Model Evaluation

This study used ANN and MLR methods. To verify the prediction results of models M1, M2, and M3, the performance indices were calculated. The outcomes of all established models are illustrated as measured and predicted values. Equations (3)–(6) were used to find the coefficient of determination (R2), root mean squared error (RMSE), variance accounts for (VAF), and a20-index of each model, respectively. Table 5 represents the performance indices of the ANN and MLR models for predicting UCS on the overall dataset, training dataset, and testing dataset.
R 2 = i = 1 n ( U C S o U C S o ) ( U C S p U C S p ) i = 1 n U C S o U C S o 2 U C S p U C S p 2
R M S E = i = 1 n ( U C S o U C S p ) 2 n
V A F = 1 v a r ( U C S o U C S p ) v a r ( U C S o ) × 100
In addition, to further assess the reliability of the model, a new engineering index, a20-index, was applied to the studied models.
a 20 i n d e x = m 20 N
where U C S o is the measured value; U C S p is the predicted value; U C S o and U C S p are the mean of the measured and predicted value, respectively; and n shows the number of the dataset. m 20 denotes the dataset with a value rate of measured UCS/predicted UCS between 0.80 and 1.20 and N represents the dataset number.

3. Prediction and Discussion of Uniaxial Compressive Strength

The main objective of this study is to investigate the capability of an intelligent model, i.e., ANN, for predicting UCS of sedimentary rocks. The actual and predicted output values were later collated and plotted to ease the performance analysis and correlation studies of these developed models. Various analytical metrics including R2, RMSE, VAF, and a20 index were used as performance criteria to examine the final output, to analyze and compare the expected models, and to evaluate the optimal model for data prediction. Model 1 (M1) is the overall dataset, model 2 (M2) consists of 70% as the training dataset and 30% as the testing dataset, and model 3 (M3) consists of 60% as the training dataset and 40% as the testing dataset.
Figure 9 indicates the predicted values of the ANN model M1 for UCS against the measured UCS for the overall dataset. The predicted correlation coefficient of M1 is R2 = 0.793. Based on the M1 predicted outputs, Figure 10a shows the aggregated comparison of predicted versus measured values for UCS. Figure 10b specifies the change in relative error between the measured and predicted values. The MSE value of model M1 achieved is 0.00599. Figure 10c illustrates the error histogram of the established model M1. Here, it can be considered that the distribution of the errors is approximately zero, which is in good agreement with the performance of model M1.
Figure 11 shows the predicted outputs of the ANN model M2 for UCS versus measured data for the training and testing data. For the training and testing data, the predicted R2 values of model M2 are 0.834 and 0.831, respectively. According to the M2 estimated results for the training data, Figure 12a displays the aggregated comparison of the predicted against measured values for UCS. Figure 12b shows the change in relative error between the measured and predicted values. The MSE value of model M2 is 0.00002. Figure 12c denotes the error histogram of model M2. It can be seen that the distribution of the errors is almost zero, which indicates that the performance of the proposed model M2 is satisfactory and reliable. Similarly, Figure 12d exhibits the aggregated comparison of the predicted against measured values for UCS of estimated outputs of M3 for the testing data. Figure 12e denotes the change in relative error between the measured and predicted values. The MSE value is achieved as 0.07657. Figure 12f represents the error histogram of model M3. Consequently, it can be seen that the distribution of the errors is nearly zero, which indicates that the performance of the proposed model M2 is acceptable.
In Figure 13, the predicted outputs of the ANN model M3 for UCS versus measured data for the training and testing data are presented. Thus, the predicted R2 values of model M3 are 0.807 and 0.775 for the training and testing data, respectively. Regarding the estimated results of M3 for the training data, Figure 14a shows the aggregated comparison of the predicted against measured values of UCS. Figure 14b shows the change in relative error between the measured and predicted values. The MSE value of M3 is 0.00015. Figure 14c signifies the error histogram of the developed model M3. Hence, it can be noted that the error distribution approaches zero, which shows that the performance of model M3 is adequate. Likewise, for predictive outputs of M3 for the testing data, Figure 14d reveals the aggregated comparison of the predicted against measured values for UCS. Figure 14e indicates the change in relative error between the measured and predicted values. The MSE value of M3 is 0.04541. Figure 14f presents the error histogram of model M3. Thus, the distribution of the errors is nearly zero, which indicates that the performance of the established model M3 is satisfactory.
The first step is to determine whether the data under consideration are appropriate for linear regression analysis. Numerous tests are suggested in the literature for this purpose. Apart from R2, another very commonly used test is the ANOVA test. In the first case, linear regression was used to determine the relationship between the dependent variable measured UCS and the three independent variables: ρ d , BTS, and ρ w e t . In Table 6, the R2 values of UCS are estimated using different equations of the MLR models, including M1, M2, and M3, for the overall dataset and training and testing data, i.e., 0.187 for M1, 0.292 and 0.066 for M2, and 0.425 and 0.062 for M3, respectively. Therefore, the R2 values of UCS are quite satisfactory in models M1, M2, and M2. Furthermore, the ANOVA test also rejected the null hypothesis at a significance value of p < 0.001.
Table 5. Performance indices of the ANN and MLR models for predicting UCS for the overall dataset, training dataset, and testing dataset.
Table 5. Performance indices of the ANN and MLR models for predicting UCS for the overall dataset, training dataset, and testing dataset.
ModelUCS
R2RMSEVAF (%)a20-index
ANNM1Overall dataset0.7930.077390.960.95
M2Train0.8340.004840.990.99
Test0.8310.276720.920.80
M3Train0.8070.012110.990.99
Test0.7750.213110.900.80
MLRM1Overall dataset0.1876.704040.981.07
M2Train0.2923.330670.770.80
Test0.0661.409500.810.99
M3Train0.4251.325180.821.05
Test0.0627.126920.990.99
Table 6. Multiple linear regression analysis for UCS in MPa; ρ d (g/cm3), BTS (MPa), and ρ w e t (g/cm3) are the dry density, Brazilian tensile strength, and wet density, respectively.
Table 6. Multiple linear regression analysis for UCS in MPa; ρ d (g/cm3), BTS (MPa), and ρ w e t (g/cm3) are the dry density, Brazilian tensile strength, and wet density, respectively.
Model CodeDatasetEquationR2
M1Overall U C S = 1.49 0.93 ρ d + 3.12 B T S + 0.26 ρ w e t 0.187
M2Train U C S = 1.04 1.11 ρ d + 4.35 B T S + 0.41 ρ w e t 0.292
Test U C S = 7.83 7.24 ρ w e t + 4.61 ρ d + 0.80 B T S 0.066
M3Train U C S = 0.72 1.80 ρ d + 0.80 ρ w e t + 6.17 B T S 0.425
Test U C S = 0.59 4.05 ρ w e t + 4.90 ρ d + 0.24 B T S 0.062

Taylor Diagram

The Taylor diagram provides a short numerical explanation of how the fit patterns match their connection and standard deviation. The Taylor diagram can be expressed as follows:
R = 1 Z z Z ( l n l ) ( m n m ) σ l σ m
where R denotes the correlation; Z denotes the discrete points; l n and m n represent two variables; σ l and σ m show the standard deviation of l and m, respectively; and l and m denote the average of σ l and σ m , respectively.
Figure 15 indicates the Taylor diagrammatic correlation between the R2, RMSE, and standard deviation of the original and predicted UCS for the M2 and M3 ANN and MLR models for the testing stage. The prediction of ANN model M3 is highly correlated with the original values and, compared with the other developed models, the standard deviation is similar to the original value. Thus, ANN model M2 with R2 = 0.831 is the most suitable for predicting UCS of sedimentary rocks in the Thar coalfield, Pakistan, among the developed models.
In an ideal scenario, the best-fit prediction model is considered as the one in which the R2 value is highest, the RMSE is lowest, the VAF is at a maximum, and the a20-index is reliable. Therefore, according to Figure 15, ANN model M2 for the testing dataset revealed the optimal results and is proposed as the best-fit prediction model for UCS in this study.

4. Sensitivity Analysis

It is crucial to accurately analyze the most important parameters that have a great influence on UCS of rock, which can certainly be problematic in the design of structures. Therefore, in this study, the cosine amplitude method was used to investigate the relative influence of the input parameters on the output [56,57]. The general formula of the adopted method can be expressed as follows:
r i j = k = 1 n ( U C S i k U C S j k ) k = 1 n U C S i k 2 k = 1 n U C S j k 2
where U C S i and U C S j are input and output values, respectively, and n denotes the dataset number during the testing stage. Finally, r i j ranges between 0 and 1, specifying additional evidence of the accuracy between each variable and the target. According to Equation (6), if the r i j of any parameter is 0, this indicates that there is no significant relationship between this parameter and the target. On the contrary, when r i j is equal to 1 or approximately 1, a significant relationship can be considered that can greatly influence UCS of the rocks.
Figure 16 shows the relationship between each input parameter ( ρ d , BTS, and ρ w e t ) of the developed model and the output (UCS). Therefore, it can be seen from the figure that BTS is the most influential parameter in predicting UCS. The corresponding coefficient values are ρ d = 0.0437, BTS = 0.485, and ρ w e t = 0.0435.

5. Conclusions

In this study, an intelligent method was used to predict the output, UCS, of sedimentary rocks collected from Block IX of the Thar coalfield, using ρd, BTS, and ρwet as input parameters. The physical and mechanical properties of rock samples were determined in a laboratory in accordance with ISRM and ASTM standards. This study determined the predictive performance of ANN and MLR models by determining the highest R2, the smallest RMSE, the highest VAF, and a reliable a20-index as follows:
For ANN models, R2, RMSE, VAF, and a20-index were 0.793, 0.07739, 0.96, and 0.95, respectively, for M1; 0.834 and 0.831, 0.00484 and 0.27672, 0.99 and 0.92, and 0.99 and 0.80, respectively, for the training and testing dataset of M2; and 0.807 and 0.775, 0.01211 and 0.21311, 0.99 and 0.90, and 0.99 and 0.80, respectively, for the training and testing dataset of M3.
In comparison, for the MLR models, R2, RMSE, VAF, and a20-index were 0.187, 6.70404, 0.98, and 1.07, respectively, for M1; 0.292 and 0.066, 3.33067 and 1.40950, 0.77 and 0.81, and 0.80 and 0.99, respectively, for the training and testing dataset of M2; and 0.425 and 0.062, 1.32518 and 7.12692, 0.82 and 0.99, and 1.05 and 0.99, respectively, for the training and testing dataset of M3.
Thus, the proposed ANN model M2 for the testing dataset yielded the optimal results and is proposed as the best-fit prediction model for UCS in this study.
Finally, by performing a sensitivity analysis, it was concluded that BTS was the most influential parameter in predicting UCS.
The current study used only ANN to predict UCS, in comparison with MLR, which could have produced more suitable results. However, future work will focus on predicting UCS using metaheuristic techniques and enhancing the accuracy of the model prediction and model performance in heterogeneous and big datasets. Moreover, the author plans to investigate UCS using optimized machine learning algorithms as well as hybrid and ensemble learning. Furthermore, some other influential attributes will be added to the UCS database to further understand the nature of this study area.

Author Contributions

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

Funding

This research was supported by the Science and Technology Innovation Project of Guizhou Province (Qiankehe Platform Talent [2019] 5620 to X.Z.). No additional external funding was received for this study.

Data Availability Statement

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

Conflicts of Interest

The author declares no competing interests.

Nomenclature

Abbreviation/SymbolParameter NameAbbreviation/SymbolParameter Name
UCSUniaxial compressive stregthR2Coefficient of determination
ISRMInternational Society of Rock MechanicsRMSERoot mean squared error
ASTMAmerican Society for Testing and MaterialsVAFVariance accounts for
ANNArtificial neural networkµPoisson’s ratio
ANFISAdaptive neuro-fuzzy interference systemρ and rDensity
PSOParticle swarm optimizationBTSBrazilian tensile strength
GAGenetic algorithmSHNSchmidt hardness
MARSmultivariate adaptive regression splinesρwetWet density
ICAImperialist competitive algorithmNPorosity
IsPoint load strength indexIs(50)Point load index
RnSchmidt hammer rebound numberVpP-wave velocity
BPIBlock punch indexAb and WabsWater absorption
DD, ρdDry density

References

  1. Madhubabu, N.; Singh, P.K.; Kainthola, A.; Mahanta, B.; Tripathy, A.; Singh, T.N. Prediction of compressive strength and elastic modulus of carbonate rocks. Measurement 2016, 88, 202–213. [Google Scholar] [CrossRef]
  2. Asheghi, R.; Shahri, A.A.; Zak, M.K. Prediction of uniaxial compressive strength of different quarried rocks using metaheuristic algorithm. Arab. J. Sci. Eng. 2019, 44, 8645–8659. [Google Scholar] [CrossRef]
  3. Abdi, Y.; Taheri-Garavand, A. Application of the ANFIS Approach for Estimating the Mechanical Properties of Sandstones. Emir. J. Eng. Res. 2020, 25, 1. [Google Scholar]
  4. Abdi, Y.; Garavand, A.T.; Sahamieh, R.Z. Prediction of strength parameters of sedimentary rocks using artificial neural networks and regression analysis. Arab. J. Geosci. 2018, 11, 587. [Google Scholar] [CrossRef]
  5. Shahri, A.A.; Asheghi, R.; Zak, M.K. A hybridized intelligence model to improve the predictability level of strength index parameters of rocks. Neural Comput. Appl. 2020, 33, 3841–3854. [Google Scholar] [CrossRef]
  6. Barzegar, R.; Sattarpour, M.; Deo, R.; Fijani, E.; Adamowski, J. An ensemble tree-based machine learning model for predicting the uniaxial compressive strength of travertine rocks. Neural Comput. Appl. 2019, 32, 9065–9080. [Google Scholar] [CrossRef]
  7. Gockceoglu, C.; Zorlu, K. A fuzzy model to predict the uniaxial compressive strength and the modulus of elasticity of a problematic rock. Eng. Appl. Artif. Intell. 2004, 17, 61–72. [Google Scholar] [CrossRef]
  8. Baykasoğlu, A.; Gϋllϋ, H.; Ҫanakҫi, A.; Ӧzbakir, A. Predicting of compressive and tensile strength of limestone via genetic programming. Expert Syst. Appl. 2008, 35, 111–112. [Google Scholar] [CrossRef]
  9. Tiryaki, B. Predicting intact rock strength for mechanical excavation using multivariate statistics, artificial neural networks and regression trees. Eng. Geol. 2008, 99, 51–60. [Google Scholar] [CrossRef]
  10. Ozcelik, Y.; Bayram, F.; Yasitli, N.E. Prediction of engineering properties of rocks from microscopic data. Arab. J. Geosci. 2013, 6, 3651–3668. [Google Scholar] [CrossRef]
  11. Rajesh-Kumar, B.; Vardhan, H.; Govindaraj, M.; Vijay, G.S. Regression analysis and ANN models to predict rock properties from sound levels produced during drilling. Int. J. Rock Mech. Min. Sci. 2013, 58, 61–72. [Google Scholar] [CrossRef]
  12. Kong, F.; Shang, J. A validation study for the estimation of uniaxial compressive strength based on index tests. Rock Mech. Rock Eng. 2018, 51, 2289–2297. [Google Scholar] [CrossRef]
  13. Teymen, A.; Mengüç, E.C. Comparative evaluation of different statistical tools for the prediction of uniaxial compressive strength of rocks. Int. J. Min. Sci. Technol. 2020, 30, 785–797. [Google Scholar] [CrossRef]
  14. Ajalloeian, R.; Kamani, M. An investigation of the relationship between Los Angeles abrasion loss and rock texture for carbonate aggregates. Bull. Eng. Geol. Env. 2019, 78, 1555–1563. [Google Scholar] [CrossRef]
  15. Cabalar, A.F.; Cevik, A.; Gokceoglu, C. Some applications of Adaptive Neuro-Fuzzy Inference System (ANFIS) in geotechnical engineering. Comput. Geotech. 2012, 40, 14–33. [Google Scholar] [CrossRef]
  16. Bashari, A.; Beiki, M.; Talebinejad, A. Estimation of deformation modulus of rock masses by sing fuzzy clustering-based modeling. Int. J. Rock Mech. Min. Sci. 2011, 48, 1224–1234. [Google Scholar] [CrossRef]
  17. Umrao, R.K.; Sharma, L.K.; Singh, R.; Singh, T.N. Determination of strength and modulus of elasticity of heterogenous sedimentary rocks: An ANFIS predictive technique. Measurement 2018, 126, 194–201. [Google Scholar] [CrossRef]
  18. Torabi-Kaveh, M.; Naseri, F.; Saneie, S.; Sarshari, B. Application of artificial neural networks and multivariate statistics to predict UCS and E using physical properties of Asmari limestones. Arab. J. Geosci. 2015, 8, 2889–2897. [Google Scholar] [CrossRef]
  19. Yagiz, S.; Sezer, E.A.; Gokceoglu, C. Artificial neural networks and nonlinear regression techniques to assess the influence of slake durability cycles on the prediction of uniaxial compressive strength and modulus of elasticity for carbonate rocks. Int. J. Numer. Anal. Methods Geomech. 2012, 36, 1636–1650. [Google Scholar] [CrossRef]
  20. Ceryan, N.; Okkan, U.; Kesimal, A. Prediction of unconfined compressive strength of carbonate rocks using artificial neural networks. Environ. Earth Sci. 2013, 68, 807–819. [Google Scholar] [CrossRef]
  21. Mohamad, E.T.; Armaghani, D.J.; Momeni, E.; Abad, S.V.A.N.K. Prediction of the unconfined compressive strength of soft rocks: A PSO-based ANN approach. Bull. Eng. Geol. Environ. 2015, 74, 745–757. [Google Scholar] [CrossRef]
  22. Aboutaleb, S.; Behnia, M.; Bagherpour, R.; Bluekian, B. Using non-destructive tests for estimating uniaxial compressive strength and static Young’s modulus of carbonate rocks via some modeling techniques. Bull. Eng. Geol. Environ. 2018, 77, 1717–1728. [Google Scholar] [CrossRef]
  23. Bejarbaneh, B.Y.; Bejarbaneh, E.Y.; Fahimifar, A.; Armaghani, D.J.; Abd Majid, M.Z. Intelligent modelling of sandstone deformation behaviour using fuzzy logic and neural network systems. Bull. Eng. Geol. Environ. 2018, 77, 345–361. [Google Scholar] [CrossRef]
  24. Fakir, M.; Ferentinou, M.; Misra, S. An investigation into the rock properties influencing the strength in some granitoid rocks of KwaZulu-Natal, South Africa. Geotech. Geol. Eng. 2017, 35, 1119–1140. [Google Scholar] [CrossRef]
  25. Han, H.; Yin, S. In-situ stress inversion in Liard Basin, Canada, from caliper logs. Petroleum 2020, 6, 392–403. [Google Scholar] [CrossRef]
  26. Skentou, A.D.; Bardhan, A.; Mamou, A.; Lemonis, M.E.; Kumar, G.; Samui, P.; Armaghani, D.J.; Asteris, P.G. Closed-Form Equation for Estimating Unconfined Compressive Strength of Granite from Three Non-destructive Tests Using Soft Computing Models. Rock Mech. Rock Eng. 2022, 11, 487–514. [Google Scholar] [CrossRef]
  27. Kaloop, M.R.; Bardhan, A.; Samui, P.; Hu, J.W.; Zarzoura, F. Computational intelligence approaches for estimating the unconfined compressive strength of rocks. Arab. J. Geosci. 2023, 16, 37. [Google Scholar] [CrossRef]
  28. Xiang, Z.; Yu, Z.; Kang, W.H.; Si, G.; Oh, J.; Canbulat, I. Estimation of in-situ rock strength from borehole geophysical logs in Australian coal mine sites. Int. J. Coal Geol. 2023, 269, 104210. [Google Scholar] [CrossRef]
  29. Köken, E.; Koca, T.K. A comparative study to estimate the mode I fracture toughness of rocks using several soft computing techniques. Turk. J. Eng. 2023, 7, 296–305. [Google Scholar] [CrossRef]
  30. Li, D.; Armaghani, D.J.; Zhou, J.; Lai, S.H.; Hasanipanah, M. A GMDH predictive model to predict rock material strength using three non-destructive tests. J. Nondestruct. Eval. 2020, 39, 81. [Google Scholar] [CrossRef]
  31. Dehghan, S.; Sattari, G.H.; Chelgani, S.C.; Aliabadi, M.A. Prediction of uniaxial compressive strength and modulus of elasticity for Travertine samples using regression and artificial neural networks. Min. Sci. Technol. 2010, 20, 41–46. [Google Scholar] [CrossRef]
  32. Yesiloglu-Gultekin, N.; Gokceoglu, C.; Sezer, E.A. Prediction of uniaxial compressive strength of granitic rocks by various nonlinear tools and comparison of their performances. Int. J. Rock Mech. Min. Sci. 2013, 62, 113–122. [Google Scholar] [CrossRef]
  33. Mohamad, E.T.; Armaghani, D.J.; Momeni, E.; Yazdavar, A.H.; Ebrahimi, M. Rock strength estimation: A PSO-based BP approach. Neural Comput. Appl. 2018, 30, 1635–1646. [Google Scholar] [CrossRef]
  34. Armaghani, D.J.; Amin, M.F.M.; Yagiz, S.; Faradonbeh, R.S.; Abdullah, R.A. Prediction of the uniaxial compressive strength of sandstone using various modeling techniques. Int. J. Rock Mech. Min. Sci. 2016, 85, 174–186. [Google Scholar] [CrossRef]
  35. Armaghani, D.J.; Mohamad, E.T.; Momeni, E.; Monjezi, M.; Narayanasamy, M.S. Prediction of the strength and elasticity modulus of granite through an expert artificial neural network. Arab. J. Geosci. 2016, 9, 48. [Google Scholar] [CrossRef]
  36. Heidari, M.; Mohseni, H.; Jalali, S.H. Prediction of uniaxial compressive strength of some sedimentary rocks by fuzzy and regression models. Geotech. Geol. Eng. 2018, 36, 401–412. [Google Scholar] [CrossRef]
  37. Yılmaz, I.; Yuksek, A.G. An example of artificial neural network (ANN) application for indirect estimation of rock parameters. Rock Mech. Rock Eng. 2008, 41, 781–795. [Google Scholar] [CrossRef]
  38. Shahani, N.M.; Zheng, X.; Guo, X.; Wei, X. MachineLearning-Based Intelligent Predictionof Elastic Modulus of Rocks at Thar Coalfield. Sustainability 2022, 14, 3689. [Google Scholar] [CrossRef]
  39. Brown, E.T. Rock Characterization Testing & Monitoring—ISRM Suggested Methods; ISRM—International Society for Rock Mechanics/Pergamon Press: London, UK, 2007; p. 211. [Google Scholar]
  40. ASTM—American Society for Tenting and Materials. Standard Practices for Preparing Rock Core as Cylindrical Test Specimens and Verifying Conformance to Dimensionaland Shape Tolerances; ASTM: West Conshohocken, PA, USA, 2013. [Google Scholar]
  41. Alexx, K. Artificial Neural Networks. Available online: https://www.computerworld.com/article/2591759/artificial-neural-networks.html#:~:text=One%20answer%20is%20to%20use,and%20learned%20to%20recognize%20objects (accessed on 3 July 2021).
  42. Alizadeh, M.; Alizadeh, E.; Asadollahpour Kotenaee, S.; Shahabi, H.; Beiranvand Pour, A.; Panahi, M.; Bin Ahmad, B.; Saro, L. Social vulnerability assessment using artificial neural network (ANN) model for earthquake hazard in Tabriz city, Iran. Sustainability 2018, 10, 3376. [Google Scholar] [CrossRef] [Green Version]
  43. Asteris, P.G.; Mokos, V.G. Concrete compressive strength using artificial neural networks. Neural Comput. Appl. 2019, 32, 11807–11826. [Google Scholar] [CrossRef]
  44. Pham, B.T.; Singh, S.K.; Ly, H.B. Using Artificial Neural Network (ANN) for prediction of soil coefficient of consolidation. Vietnam J. Earth Sci. 2020, 42, 311–319. [Google Scholar] [CrossRef]
  45. Pham, T.A.; Ly, H.B.; Tran, V.Q.; Giap, L.V.; Vu, H.L.T.; Duong, H.A.T. Prediction of pile axial bearing capacity using artificial neural network and random forest. Appl. Sci. 2020, 10, 1871. [Google Scholar] [CrossRef] [Green Version]
  46. Le, V.M.; Pham, B.T.; Le, T.T.; Ly, H.B.; Le, L.M. Daily rainfall prediction using nonlinear autoregressive neural network. Micro-Electron. Telecommun. Eng. 2020, 106, 213–221. [Google Scholar]
  47. Le, T.T.; Pham, B.T.; Le, V.M.; Ly, H.B.; Le, L.M. A robustness analysis of different nonlinear autoregressive networks using Monte Carlo simulations for predicting high fluctuation rainfall. Micro-Electron. Telecommun. Eng. 2020, 106, 205–212. [Google Scholar]
  48. Ly, H.B.; Le, T.T.; Vu, H.L.T.; Tran, V.Q.; Le, L.M.; Pham, B.T. Computational hybrid machine learning based prediction of shear capacity for steel fiber reinforced concrete beams. Sustainability 2020, 12, 2709. [Google Scholar] [CrossRef] [Green Version]
  49. Rashidian, V.; Hassanlourad, M. Application of an artificial neural network for modeling the mechanical behavior of carbonate soils. Int. J. Geomech. 2014, 14, 142–150. [Google Scholar] [CrossRef]
  50. Fidan, S.; Oktay, H.; Polat, S.; Ozturk, S. An artificial neural network model to predict the thermal properties of concrete using different neurons and activation functions. Adv. Mater. Sci. Eng. 2019, 2019, 3831813. [Google Scholar] [CrossRef] [Green Version]
  51. Gowida, A.; Moussa, T.; Elkatatny, S.; Ali, A. A hybrid artificial intelligence model to predict the elastic behavior of sandstone rocks. Sustainability 2019, 11, 5283. [Google Scholar] [CrossRef] [Green Version]
  52. Hajihassani, M.; Armaghani, D.J.; Sohaei, H.; Mohamad, E.T.; Marto, A. Prediction of airblast-overpressure induced by blasting using a hybrid artificial neural network and particle swarm optimization. Appl. Acoust. 2014, 80, 57–67. [Google Scholar] [CrossRef]
  53. Ekemen Keskin, T.; Özler, E.; Şander, E.; Düğenci, M.; Ahmed, M.Y. Prediction of electrical conductivity using ANN and MLR: A case study from Turkey. Acta Geophys. 2020, 68, 811–820. [Google Scholar] [CrossRef]
  54. Sajid, M.J. Modelling best fit-curve between China’s production and consumption-based temporal carbon emissions and selective socio-economic driving factors. IOP Conf. Series Earth Environ. Sci. 2020, 431, 012061. [Google Scholar] [CrossRef]
  55. Sajid, M.J. Machine Learned Artificial Neural Networks Vs Linear Regression: A Case of Chinese Carbon Emissions. IOP Conf. Series Earth Environ. Sci. 2020, 495, 012044. [Google Scholar] [CrossRef]
  56. Momeni, E.; Nazir, R.; Armaghani, D.J.; Maizir, H. Prediction of pile bearing capacity using a hybrid genetic algorithm-based ANN. Measurement 2014, 57, 122–131. [Google Scholar] [CrossRef]
  57. Ji, X.; Liang, S.Y. Model-based sensitivity analysis of machining-induced residual stress under minimum quantity lubrication. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2017, 231, 1528–1541. [Google Scholar] [CrossRef]
Figure 1. Geological site of the collected rock samples [38].
Figure 1. Geological site of the collected rock samples [38].
Mathematics 11 01650 g001
Figure 2. (a) Deformed rock core specimen for Brazilian tensile strength test and (b) deformed rock core specimen for UCS test.
Figure 2. (a) Deformed rock core specimen for Brazilian tensile strength test and (b) deformed rock core specimen for UCS test.
Mathematics 11 01650 g002
Figure 3. Histogram plots of the original dataset in this study: (a) dry density (g/cm3), (b) BTS (MPa), (c) wet density (g/cm3), and (d) UCS (MPa).
Figure 3. Histogram plots of the original dataset in this study: (a) dry density (g/cm3), (b) BTS (MPa), (c) wet density (g/cm3), and (d) UCS (MPa).
Mathematics 11 01650 g003
Figure 4. Correlation plot of inputs (dry density (g/cm3), BTS (MPa), and wet density (g/cm3)) and output (UCS (MPa)) of the original dataset in this study.
Figure 4. Correlation plot of inputs (dry density (g/cm3), BTS (MPa), and wet density (g/cm3)) and output (UCS (MPa)) of the original dataset in this study.
Mathematics 11 01650 g004
Figure 5. Flow chart of the predictive modeling process for UCS.
Figure 5. Flow chart of the predictive modeling process for UCS.
Mathematics 11 01650 g005
Figure 6. Structure of the artificial neural network.
Figure 6. Structure of the artificial neural network.
Mathematics 11 01650 g006
Figure 7. Validation performance curves of UCS at (a) M1, (b) M2, and (c) M3.
Figure 7. Validation performance curves of UCS at (a) M1, (b) M2, and (c) M3.
Mathematics 11 01650 g007
Figure 8. ANN training scatter plots of predicted UCS against measured UCS for (a) M1, (b) M2, and (c) M3.
Figure 8. ANN training scatter plots of predicted UCS against measured UCS for (a) M1, (b) M2, and (c) M3.
Mathematics 11 01650 g008
Figure 9. ANN model M1 results for UCS plotted against the measured data.
Figure 9. ANN model M1 results for UCS plotted against the measured data.
Mathematics 11 01650 g009
Figure 10. The demonstration of ANN model M1 for UCS. (a) Model M1 results aggregated with measured UCS. (b) The variation in error between the measured and predicted values. (c) Error histogram.
Figure 10. The demonstration of ANN model M1 for UCS. (a) Model M1 results aggregated with measured UCS. (b) The variation in error between the measured and predicted values. (c) Error histogram.
Mathematics 11 01650 g010
Figure 11. ANN model M2 results for UCS plotted against the measured data for the (a) training and (b) testing data.
Figure 11. ANN model M2 results for UCS plotted against the measured data for the (a) training and (b) testing data.
Mathematics 11 01650 g011
Figure 12. The demonstration of ANN model M2 for UCS. (a) The model M2 results aggregated with the measured UCS. (b) The variation in error between the measured and predicted values. (c) Error histogram for the training data and (d) model M2 results aggregated with the measured data. (e) The variation in error between the measured and predicted values. (f) Error histogram for the testing data.
Figure 12. The demonstration of ANN model M2 for UCS. (a) The model M2 results aggregated with the measured UCS. (b) The variation in error between the measured and predicted values. (c) Error histogram for the training data and (d) model M2 results aggregated with the measured data. (e) The variation in error between the measured and predicted values. (f) Error histogram for the testing data.
Mathematics 11 01650 g012
Figure 13. ANN model M3 results for UCS plotted against the measured data for the (a) training and (b) testing data.
Figure 13. ANN model M3 results for UCS plotted against the measured data for the (a) training and (b) testing data.
Mathematics 11 01650 g013
Figure 14. The demonstration of ANN model M3 for UCS. (a) Model M3 results aggregated with the measured UCS. (b) The variation in error between the measured and predicted values. (c) Error histogram for the training data and (d) model M3 results aggregated with the measured data. (e) The variation in error between the measured and predicted values. (f) Error histogram for the testing data.
Figure 14. The demonstration of ANN model M3 for UCS. (a) Model M3 results aggregated with the measured UCS. (b) The variation in error between the measured and predicted values. (c) Error histogram for the training data and (d) model M3 results aggregated with the measured data. (e) The variation in error between the measured and predicted values. (f) Error histogram for the testing data.
Mathematics 11 01650 g014
Figure 15. Demonstration of the Taylor diagram for the testing data based on ANN and MLR.
Figure 15. Demonstration of the Taylor diagram for the testing data based on ANN and MLR.
Mathematics 11 01650 g015
Figure 16. The effect of input variables on the result of the established model.
Figure 16. The effect of input variables on the result of the established model.
Mathematics 11 01650 g016
Table 1. Previous studies using intelligent methods to predict UCS.
Table 1. Previous studies using intelligent methods to predict UCS.
MethodInputOutputR2References
ANNn, Is, µ, ρ, VpUCS0.97(Madhubabu et al., 2016) [1]
ANNρ, n, Vp, AbUCS0.93(Abdi et al., 2018) [4]
ANNn, r, WabsUCS0.92(Kamani et al., 2020) [14]
ANNVp, Is(50), BTSUCS0.97(Mohamad et al., 2015) [21]
ANNRn, Vp, DDUCS0.82(Li et al., 2020) [30]
ANNIs, Vp, Rn, nUCS0.93(Dehghan et al., 2010) [31]
ANFISBTS, VpUCS0.60(Yesiloglu-Gultekin et al., 2013) [32]
PSO-BPDD, MC, Vp, Is(50), Id2UCS0.999(Mohamad et al., 2018) [33]
ICA-ANNRn, Vp, Is(50)UCS0.949(Armaghani et al., 2016a) [34]
ICA-ANNn, Rn, Vp, Is(50)UCS0.915(Armaghani et al., 2016b) [35]
MLRn, Is, µ, ρ, VpUCS0.91(Madhubabu et al., 2016) [1]
MLRρ, n, Vp, AbUCS0.88(Abdi et al., 2018) [4]
MLRVp, IS(50), SHN, BPIUCS0.91(Heidari et al., 2018) [36]
MLRId2, Is(50), N, éUCS0.58(Yılmaz et al., 2008) [37]
Table 2. Physical and mechanical parameters of the dataset.
Table 2. Physical and mechanical parameters of the dataset.
Datasetρd (g/cm3)BTS (MPa)ρwet (g/cm3)UCS (MPa)
11.910.3052.130.404
21.750.2172.010.491
31.770.3182.040.531
41.780.27120.579
51.760.2922.040.557
741.810.1782.10.541
751.840.1892.110.476
761.960.22.180.508
771.780.1082.090.511
781.840.1382.091.415
Table 3. The minimum, maximum, average, and standard deviation of the dataset.
Table 3. The minimum, maximum, average, and standard deviation of the dataset.
Parametersρd (g/cm3)BTS (MPa)ρwet (g/cm3)UCS (MPa)
Minimum1.220.0231.630.304
Maximum2.120.6272.33.55
Average1.760.322.041.38
Standard deviation0.220.130.150.98
Table 4. The dataset distribution for the ANN and MLR models.
Table 4. The dataset distribution for the ANN and MLR models.
Model CodeDatasetDataset Distribution (%)Total Dataset
Model 1 (M1)Overall10078
Model 2 (M2)Train7055
Test3023
Model 3 (M3)Train6047
Test4031
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wei, X.; Shahani, N.M.; Zheng, X. Predictive Modeling of the Uniaxial Compressive Strength of Rocks Using an Artificial Neural Network Approach. Mathematics 2023, 11, 1650. https://doi.org/10.3390/math11071650

AMA Style

Wei X, Shahani NM, Zheng X. Predictive Modeling of the Uniaxial Compressive Strength of Rocks Using an Artificial Neural Network Approach. Mathematics. 2023; 11(7):1650. https://doi.org/10.3390/math11071650

Chicago/Turabian Style

Wei, Xin, Niaz Muhammad Shahani, and Xigui Zheng. 2023. "Predictive Modeling of the Uniaxial Compressive Strength of Rocks Using an Artificial Neural Network Approach" Mathematics 11, no. 7: 1650. https://doi.org/10.3390/math11071650

APA Style

Wei, X., Shahani, N. M., & Zheng, X. (2023). Predictive Modeling of the Uniaxial Compressive Strength of Rocks Using an Artificial Neural Network Approach. Mathematics, 11(7), 1650. https://doi.org/10.3390/math11071650

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