Next Article in Journal
Experimental Study on Hot Spot Stresses of Curved Composite Twin-Girder Bridges
Previous Article in Journal
Hyperbranched Polymer Network Based on Electrostatic Interaction for Anodes in Lithium-Ion Batteries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting the Engineering Properties of Rocks from Textural Characteristics Using Some Soft Computing Approaches

1
School of Earth Sciences, Damghan University, Damghan 36716-41167, Iran
2
Department of Geology and Pole of CGeo—Geoscience Center, University of Trás-os-Montes e Alto Douro, 5000-801 Vila Real, Portugal
*
Author to whom correspondence should be addressed.
Materials 2022, 15(22), 7922; https://doi.org/10.3390/ma15227922
Submission received: 26 September 2022 / Revised: 17 October 2022 / Accepted: 3 November 2022 / Published: 9 November 2022
(This article belongs to the Section Construction and Building Materials)

Abstract

:
Rock is used as a foundation and building material in many engineering projects and it is important to determine/predict its engineering properties before project construction. Petrographic and textural characteristics are useful parameters for predicting engineering properties of rocks in such applications. In this research, fifteen rock samples were taken and their engineering characteristics, namely dry and saturated unit weights, porosity, water absorption, slake durability index (SDI), Schmidt rebound hardness (SRH), ultrasonic P-wave velocity (UPV), and uniaxial compressive strength (UCS), were measured in the laboratory. Petrographic and textural characteristics of the rocks, determined from thin section and X-ray diffraction investigations, led to the evaluation of the texture coefficient (TC). Based on simple regression analysis (SRA), the TC values have direct relationships with density, SDI, SRH, UPV, and UCS, and inverse relationships with porosity and water absorption. Experimental models were developed using multiple regression analysis (MRA) and artificial neural network (ANN) to predict Id2, SRH, UPV, and UCS of the tested rocks from the values of TC. Some statistical parameters including Pearson regression coefficient (R), coefficient values account for (VAF), root mean square error (RMSE), mean absolute percentage error (MAPE), and performance index (PI) were calculated to assess the performances of the MRA and ANN models. The correlations between experimental and calculated values of Id2, SRH, UPV, and UCS indicated that predicted values of the ANN models are more valid than the MRA. Additionally, the residual error of the ANN models varies less than the MRA. Finally, it has been concluded that the SRA, MRA, and ANN methods can successfully predict the rock engineering properties from the TC.

1. Introduction

Natural stones are used worldwide at an increasing rate because of their availability, excellent physical-mechanical properties, and diversity of textural characteristics [1,2]. Several factors are considered during an exploration campaign of building stones, such as the fracturing in the quarry, which defines if a rock mass can be exploited, or the durability, which delimits the possible applications [3,4,5,6,7,8]. The evaluation of the geotechnical properties of building stones, such as color and brightness, is a key factor to define their utilization [9,10,11,12,13,14]. Furthermore, engineering characteristics of different rocks are essential parameters for design and construction of various engineering projects such as dams, tunnels, caverns, and foundations either on or inside rocks. These characteristics could be determined practically through direct methods as described in ISRM [15]. Applying direct tests to assess these properties is usually expensive, time consuming, and requires regularly shaped samples. Therefore, in such cases, experimental approaches and empirical equations can be used to predict the engineering properties of rocks from simple and index test results.
Several researchers have attempted to develop various soft computing models for predicting different parameters from others in material sciences [16,17,18,19,20,21,22,23,24,25] and engineering properties of different rock types from their petrographic characteristics in engineering geology and rock mechanics [26,27,28,29,30,31,32,33,34,35]. In the recent years, some research works have been performed to assess correlations between mineralogical and textural characteristics and mechanical properties of different rocks by using statistical analyses and different soft computing approaches such as genetic programing (GP), artificial neural network (ANN), adaptive neuro-fuzzy inference system (ANFIS), and support vector machine (SVM) [36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52]. The main advantages of these approaches are that (i) they have made it possible to solve nonlinear problems, in which mathematical models are not available, and (ii) they have introduced human knowledge such as cognition, recognition, understanding, learning, and others in the fields of computing [53]. In this regard, Brace [54] recognized that fine-grained rocks have higher strength than coarse-grained ones. Ulusay et al. [55] tested sandstones in a laboratory and found that textural properties are important parameters for predicting engineering properties of these rocks. Tugrul and Zarif [56] stated that textural characteristics are more important than mineral content for estimating the engineering properties of rocks. Eberli et al. [57] observed that the size and shape of grains, grain size distribution, and the ratio of grain to matrix have an influence on the acoustic wave velocity through the unconsolidated carbonate sediments. Meng and Pan [58] and Khanlari et al. [59] indicated that physical properties (i.e., density, porosity, permeability) and mechanical characteristics (such as strength, deformability, durability, and hardness) of rocks are affected by mineral composition, grain size, grain contact, and rock cement. Jensen et al. [38] found that porosity, grain size, weakness planes, and microfractures presented in limestones’ texture have important effects on the strength of the rocks.
The texture coefficient (TC) was developed by Howarth and Rowlands [60] to understand the effect of textural characteristics on the physical and mechanical properties of rocks. They investigated the relations between the mechanical properties and TC for igneous rocks, marbles, and sandstones and found close relationships between mechanical rock properties and the TC. In many research studies, the TC has been applied to predict the engineering properties of rocks [61,62,63,64]. Singh et al. [36] used artificial neural networks to estimate the strength of schistose rocks from their petrographic characteristics. In their research, the network was trained to predict axial point load, uniaxial compression, and tensile strengths from mineralogical content, grain size, aspect ratio, surface weight, form factor, and foliation orientation. Alber and Kahraman [65] used regression analysis to predict the elasticity modulus and uniaxial compressive strength of a fault breccia from the TC. The results of their research remarked that the uniaxial compressive strength of the tested breccia could be estimated from the texture coefficient. Manouchehrian et al. [42] applied artificial neural networks and multivariate regression for assessing the uniaxial compressive strength based on textural characteristics of rocks. Ersoy and Acar [66] studied the petrographic and textural properties’ influences on the strength of granites and found that the size of mineral has a higher effect than mineral type on rock strength. Aligholi et al. [67] evaluated relationships between drillability indices on one hand and index properties and petrographic data on the other hand of hard igneous rocks. They revealed that the multiple regression models prepared by using petrographic features provide a better prediction of drillability in comparison to those prepared by using index properties. Akram et al. [68], by assessing petrographic and mechanical characteristics of Sakesar limestone, attempted to establish the relationships between the parameters and the strength parameters such as uniaxial compressive strength, point load index, and Schmidt rebound hardness values. In their research, regression analyses were used to establish correlations between limestone constituents and the parameters of strength. Kolay and Baser [69] focused on the relationships between internal structure and engineering parameters of basalt rocks by digitizing the textural properties. In this research, simple regression analysis was performed using the laboratory results incorporating texture coefficient and engineering parameters including dry unit weight, P-wave velocity, Schmidt rebound hardness, point load strength index, and uniaxial compressive strength of the rocks. Their results showed that the strength of the basalts has direct relation by their texture coefficient.
All the above studies show that the petrographic characteristics and rock texture affect the engineering properties of different rocks. These properties can be easily assessed by routine thin section studies in a laboratory. As sedimentary rocks have a very diverse nature and suitable texture in the presence of grains with different sizes and shapes in microscopic studies, the effect of petrographic and textural characteristics on engineering properties of sedimentary rocks is sensible. Therefore, this research focuses on the assessment of engineering properties and the textural coefficient of the rocks to predict the engineering properties from petrographic and textural characteristics using statistical methods and artificial neural networks. So, the main objective of the research is to establish relationships between the textural characteristics of rocks as a primary and essential parameter from one side and engineering properties as advanced and perfect parameters from the other by using careful approaches such as soft computing methods. Additionally, another benefit of this research is that the use of soft computing methods for predicting rock engineering characteristics is a cheap and quick method, especially for soft rocks such as sedimentary rocks, for which laboratory tests are difficult. On the other hand, this research can be considered as a model for future research and its results can be used by other researchers. These issues motivated this research. Therefore, fifteen sedimentary rock samples including limestones and sandstones have been studied to predict their engineering properties from petrographic and textural characteristics using statistical analysis and soft computing methods.

2. Geological Setting and Methodology

The study area and sampling locations are placed in the north and northwest of the Damghan region in the north of Iran. This region has a continental climate and an irregular topography related to its geological history. The area is situated in the Alborz-Azarbayejan structural zone and contains Paleozoic to Neogene outcropped rocks. Suitable sampling locations were found from the 1:100,000 geological map presented by the Geological Society of Iran [70]. Fifteen rock samples including limestones and sandstones were collected for performing the study. Sampling was carried out from quarries, trenches, and road cuttings during field investigations. The limestone samples include Sibzar (SBZ), Elika (ELK), Dalichai (DLC), and Lar (LAR) and the sandstone samples contain Padeha (PDH), Barut (BRT), and Shemshak (SMK). SBZ samples from the Middle Devonian age were collected from the Bahram and Sibzar formations from the northern part of a ballast mine in the region. ELK samples from the Lower Triassic age were collected from the Elika formation from a new trench near the Astaneh village. DLC and LAR samples, which depend on the Dalichai and Lar formations of the Upper Jurassic age were collected from road cuttings. The samples of PDH were obtained from a ballast mine 10 km from the Chesmeh-Ali road. The age of the Padeha formation is Devonian. The samples of BRT were obtained from the Barut formation with a Lower Cambrian age. The samples of SMK were collected from the Shemshak formation outcropped near a travertine quarry. These samples are from the Middle Jurassic age. Figure 1 illustrates the geological properties and the sampling locations.
After the field investigations, a regular laboratory test program was carried out to assess engineering properties of the collected rocks. XRD analyses and microscopic investigations were performed based on ASTM [71] and ISRM [15] to determine the rocks’ constituent minerals, petrography properties, textural characteristics, and finally texture coefficient (TC). X-ray diffraction (XRD) analysis was carried out by X’Pert HighScore plus (v.3.0.5) software. Composition materials of the studied rocks were identified in the 2θ ranges from 4° to 70°. The dry unit weight, saturated unit weight, effective porosity, water absorption, P-wave velocity, Schmidt hardness, slake durability index, and uniaxial compressive strength of the rocks were measured. The test plan for the present study is summarized in Figure 2. These tests were conducted on standard NX core specimens (prepared based on ISRM [15]) presented in Figure 3. Additionally, the used apparatuses and testing procedures in this study are shown in Figure 4. Physical properties including dry and saturated unit weights (γdry, γsat), effective porosity (n), and water absorption (Wa) were determined in accordance with ISRM [15]. The slake durability index (SDI) test was carried out in accordance with ISRM [15] and ASTM [72] up to 3 cycles. This test was performed on 150 irregular or lump prepared specimens. The Schmidt rebound hardness (SRH) test was carried out according to ISRM [73] and ASTM [74] on blocks of the rock samples. The ultrasonic P-wave velocity (UPV) test was conducted based on ASTM [75] in the laboratory. The uniaxial compressive strength test was conducted in accordance with ISRM [15] and ASTM [76]. Five experiments for each rock sample were performed on cylindrical specimens with a length to diameter ratio of 2:2.5 to determine the average values of UCS. Finally, the estimation of engineering properties of the rocks was performed from petrographic and textural characteristics using statistical methods and artificial neural networks.

3. Results of the Laboratory Investigations

3.1. Mineralogical and Petrographic Studies

Based on the obtained results, the studied limestones are generally composed of calcite, and the studied sandstones are quartz-rich (see Table 1 and Figure 5). The SBZ samples are sparry limestone (type I) or dolomitic limestone. These samples mainly contain calcite, dolomite, and quartz. The original calcite matrix has been replaced by dolomite, but the micritic allochems are partly replaced. Where replacement is incomplete, euhedral rhombe-shaped crystals are visible. ELK samples are mainly composed of ooids, peloids, and many fossils such as alga and echinoderms. The allochems are a mixture of ooids and bioclasts. Therefore, these samples are oosparite (type I) according to the Folk or ooid grainstone according to the Dunham classifications. Since the allochems are rounded, it would be a rounded oosparite using the Folk textural spectrum. DLC samples are sparry allochemical limestone (type I) or sandy grainstone. LAR samples are microcrystalline limestone (type III) without allochems. These samples are mainly composed of ooids, peloids, alga, and echinoderms. According to the Folk classification, PDH samples are typical greywacke that are dominantly composed of calcite, quartz, and dolomite. Poor sorting, abundant quartz monocrystalline, sedimentary origin, and the existence of rhombohedra dolomite are the main textural characteristics of these samples. Additionally, the rock matrix is composed of a hematite and dolomite fill between the fine-grain quartzes in these samples. BRT samples are litharenite and are composed of quartz, calcite, orthoclase, muscovite and others rock fragments such as chert. Quartz is the main mineral in these samples and calcite and feldspar are accessory minerals. SMK samples are classified as sub-litharenite and are composed of quartz (more than 70%) as the main mineral, and calcite and feldspar as accessory minerals. In these samples, the matrix fill is calcite and hematite between quartzes. The mineralogical and petrographic properties of the studied rocks are summarized in Table 1.

3.2. XRD Analysis

X-ray diffraction (XRD) analysis was performed and the provided diffractograms are shown in Figure 6. The results indicated that calcite and quartz are the two main minerals in all studied rocks. Other minerals such as hematite are present as secondary minerals in the rock samples. Similar to the thin section results, the limestones and sandstones are commonly composed of calcite and quartz, respectively. So, the results of both thin section studies and XRD analysis are completely coordinated with each other.

3.3. Engineering Properties Evaluation

The values of all physical properties of the studied rock samples including dry and saturated unit weights (γdry, γsat), effective porosity (n), and water absorption (Wa) are presented in Table 2. The values of γdry are between 2.35 and 2.60 for the samples of PDH3 and SMK2, respectively. Additionally, the values of porosity are between 2 and 10.32 for the samples of SMK2 and PDH3, respectively. Based on the IAEG [77] rock classification, most of the samples were placed in the rock groups with moderate to high unit weight and low to medium porosity. The samples of PDH3 and SMK2 have the lowest and highest values of the slake durability index equal to 95.54% and 99.21% at the third cycle (Id3), respectively (Table 2). According to the Gamble [78] classification, the tested samples were placed in durable and very durable rock groups. The lowest and highest values of the Schmidt rebound hardness (SRH) are between 24 and 49, which belong to the samples of PDH3 and SMK2, respectively. The value of ultrasonic P-wave velocity (UPV) ranges from 3.09 to 5.75 km/s. The samples PDH3 and LAR2 have the highest and lowest values of wave velocities, respectively. The minimum and maximum values of the uniaxial compressive strength (UCS) are between 16.09 and 127.43 MPa, which belong to the samples of PDH3 and SMK2, respectively. According to Broch and Franklin [79], the tested samples were placed in weak to extremely strong rock groups. The values of engineering properties for the studied rock samples are presented in Table 2.
Some comprehensive linear correlations were obtained between the physical and mechanical properties, which are presented in Table 3. Regression coefficients (R) as Pearson correlations were found between 0.77 and 0.98. Significant level values (Sig.) were between 0.000 and 0.001, and the high values of R show that the determined engineering parameters are significantly related to each other. Based on the results, Gs, γdry, γsat, SDI, SRH, UPV, and UCS have direct linear relations, whereas n and Wa have inverse linear relations to the other parameters. These results are comparable to the results presented by other authors [34,67,69,80].

3.4. Texture Coefficient (TC) Calculation

The texture of rocks was defined as “the degree of crystallinity, grain size or granularity, and the fabric or geometrical relationships between the constituents of a rock” [81]. The texture coefficient (TC) proposed by Howard and Rowlands [60] is a very good factor to quantify rock texture. The values of TC, which is the so-called rock fabric coefficient by Prikryl [82], were evaluated by performing four processes, namely (1) measuring grain circularity, (2) measuring grain elongation, (3) measuring and quantifying grain orientation, and (4) weighting of results based on the degree of grain packing [61]. After passing the processes, the TC is calculated from the following equation [61]:
TC = AW [ ( N 0 N 0 + N 1 × 1 FF 0 ) + ( N 1 N 0 + N 1 × AR 1 × AF 1 ) ]
Here, TC is the texture coefficient, AW is the grain packing weighting, N0 is the number of grains whose aspect ratio is below a pre-set discrimination level, N1 is the number of grains whose aspect ratio is above a pre-set discrimination level, FF0 is the arithmetic mean of discriminated form factors, AR1 is the arithmetic mean of discriminated aspect ratios, and AF1 is the angle factor, quantifying grain orientation.
The AW, AR, and FF are calculated by the following equations [61]:
AW =   ( grain   areas   within   the   reference   area   boundary ) ( area   boundary   by   the   reference   area   boundary )
AR = L W
FF = 4 π A P 2
In these equations, L is length, W is width, P is the perimeter, and A is the area. The angle factor (AF) is calculated by the equations below [61]:
AF = i = 1 9 ( x i N ( N 1 ) 2 ) × i
AF 1 = AF 5
where xi is the number of angular differences in each class, N is the total number of elongated grains, and i is the weighting factor and class number.
Calculation of TC is possible from TIFF format of microscopic images of rocks using the JMicroVision (v1. 27) program. This is a valuable software program for determining the size, shape, orientation, and texture properties of different rocks. The JMicroVision (v1. 27) environment for calculating TC is presented in Figure 7. The values of TC calculated from Equation (1) for the studied rock samples are presented in Table 4. Additionally, requiremed parameters such as AW, AR, FF, and AF were determined and are summarized in this table. Based on the results, TC values were obtained from “0.42 to 1.06”. The samples of PDH3 and LAR2 have the lowest and highest values of TC, respectively. A comparative diagram of the values of TC for the studied rocks is shown in Figure 8a. Based on the histogram, the TC has an average value of 0.79, which is presented in Figure 8b. The standard deviation (Std. Dev.) values are very low, which shows a close dispersion of the TC values (Figure 8b). The parameters of AW, AR1, and TC have direct relations. Therefore, the value of TC will increase by increasing AW and AR1. The results are confirmed by the results presented by Howarth and Rowlands [60], Howarth and Rowlands [61], Ersoy and Waller [62], Gupta and Sharma [63], Tandon and Gupta [64], and Kolay and Baser [69].

4. Data Analysis

4.1. Simple Regression Analysis (SRA)

The statistical and regression analyses were conducted by using the IBM SPSS Statistics v. 24 [83] computer program. The index, physical, and mechanical properties including dry and saturated unit weight, effective porosity, water absorption, slake durability index, Schmidt hardness, P-wave velocity, and uniaxial compressive strength were considered as independent variables and the TC was used as the dependent variable. Graphs of these correlations are presented in Figure 9. Additionally, the obtained predictive models and coefficient of determination (R2) are presented in Table 5. The R2 shows how the values will fit the data. It is calculated by the following equation:
R 2 = 1 i = 1 N ( y   y ) 2 i = 1 N ( y   y ˜ ) 2
The results of this evaluation show that R2 is between 0.77 and 0.94 (similar to Gupta and Sharma [63] and Khanlari et al. [59]). Figure 8a shows direct exponential relationships between dry and saturated unit weights (γd, γs) and TC. High coefficients of determination (R2) of 0.80 and 0.83 were obtained from these evaluations. This means that the density of the studied rock samples is greatly affected by the values of TC. Relationships between porosity and water absorption and TC are inverse linear relations (Figure 9b) with high coefficients of determination (R2) of 0.86 and 0.88, respectively. Figure 9c shows relations between slake durability indexes and TC. For this purpose, Id1, Id2, and Id3 were separately correlated to TC. There were direct power relations with good coefficients of determination (R2) of 0.77, 0.78, and 0.82, respectively. As can be seen, when the cycles of the slake durability test increase, the coefficient of determination increases. Therefore, correlations between SDI and TC present better results in upper cycles of the test. As can be seen in Figure 9d, a good direct logarithmic relationship is recognizable between SRH and TC with a high coefficient of determination (R2) of 0.82. A direct power relation was found between TC and UPV with a coefficient of determination (R2) of 0.85, which is shown in Figure 9e. Figure 9 shows a direct linear relation between UCS and TC with a very good coefficient of determination (R2) of 0.94. Direct relationships of SDI, SRH, UPV, and UCS with TC indicated that the engineering properties of the rocks are well correlated to texture coefficients (TC). A three-dimensional view of the correlations among texture coefficient (TC), water absorption (Wa), porosity (n), slake durability index (Id2), Schmidt rebound hardness (SRH), and uniaxial compressive strength (UCS) is shown in Figure 10. This figure can simultaneously show the relationship between the TC and two other parameters and expose the simultaneous relationship of these parameters with each other.
To assess the performances of the models, the statistical parameters such as the regression coefficient (R), the root mean square error (RMSE), the coefficient values accounted for (VAF), the mean absolute percentage error (MAPE), and the performance index (PI) were evaluated by the following equations [84]:
R = n (   xy ) (   x ) (   y ) [ n (   x 2 ) (   x ) 2 ] [ n (   y 2 ) (   y ) 2 ]
RMSE = ( 1 N ) × i = 1 N ( y   y ) 2
VAF = [ 1 var y   y var y ] × 100
MAPE = [ 1 N i = 1 N | y   y y | ] × 100
PI = [ R 2 + ( VAF 100 ) RMSE ]
In these equations, y and y are the experimental and calculated values, respectively, and N is the total number of points. The models are excellent (the predicted values from the regression equation are equal to the experimental values) if R = 1, RMSE = 0, VAF = 100, MAPE = 0, and PI = 100. Ten predictive models were developed by simple regression analysis (Table 5). The values of parameters related to the model and significance level (Sig.) are presented in Table 6. R values obtained between 0.88 and 0.97 show the high validity of the experimental equations. RMSE values are very low in models 1 and 2, equal to 0.03, that correlate γd and γs to the TC. High VAF values were obtained from the models of SRH and UCS. However, mean absolute percentage error values are very low in the SDI models, especially in the Id2 model. Performance assessment by PI indicates model 8 has best performance among all models. This model correlates SRH to TC. The significance level for all developed models was calculated equal to 0.000. This indicates that the obtained predictive models by SRA with a 99% confidence limit predict the values of the independent variable when the TC is considered as an input. The results the simple regression analysis revealed that texture coefficients are good parameters for evaluating and predicting the physical, index, and mechanical properties of the studied rock samples.

4.2. Multiple Regression Analysis (MRA)

The multiple regression analysis was performed by IBM SPSS Statistics v. 24 [83]. The multiple regression analysis creates a correlation between more than two parameters, and it results in a more accurate correlation with higher values of correlation coefficients (R) and coefficients of determination (R2) than simple regression analysis. In the multiple regression analysis, physical properties, namely, dry unit weight and porosity that are determined by a simple and cheap test, and also texture coefficient values were developed in this research and were considered as inputs or independent variables. Other engineering properties including Id2, SRH, UPV, and UCS were considered as dependent variables or outputs. Therefore, by three inputs of γd, n, and TC in four models, the values of Id2, SRH, UPV, and UCS were predicted.
The summary of coefficients related to multiple regression models, namely, Beta coefficient, standard error of Beta, significance level, lower and upper confidence interval, tolerance, and VIF, is presented in Table 7. The variance inflation factor (VIF) and tolerance are the statistical parameters that show multi-collinearity of models. Based on the results, tolerance values are obtained between 0 and 1. Therefore, the issue of collinearity does not exist in MR models.
Table 8 shows the Rs values such as R, R square, adjusted R, standard error, R Square Change, and Dubrin–Watson for the models. The results in this table show that model 4 has greater R, which estimates UCS with γdry, n, and TC. The Dubrin–Watson (DW) statistic is a number that tests for autocorrelation in the residual form of a statistic regression analysis. The DW statistic is always between 0 and 4. When the value is 0, it indicates negative autocorrelation.
Predictive models that were obtained from the multiple regression analysis and statistical parameters for model performance are presented in Table 9. The predictive model for Id2 has lower values of RMSE and MAPE. Model 4, which predicted UCS, has better values of VAF and PI, whereas this model has higher values of RMSE and MAPE. Significance level values indicated that all models predicted output values with a 99% confidence limit. γdry, n, and TC are good parameters to predict the engineering characteristics of the studied rocks based on the regression coefficients of multiple regression analysis. This result is similar to the results of Singh et al. [85], Alber and Kahraman [65], Manouchehrian et al. [42], Ersoy and Acar [66], and Akram et al. [68].

4.3. Artificial Neural Network (ANN)

The artificial neural network is a mathematical and computational approach that is taken from the structural and functional aspects of biological neural networks [86]. It has been widely adopted by scientists due to its ability to develop complex nonlinear models and is applied for solving a wide variety of applications [87]. In the present research, the feed-forward neural network was used to estimate the engineering characteristics of the studied rocks. Similar to the multiple regression models, γdry, n, and TC are constant input parameters in all models. Id2, SRH, UPV, and UCS were considered as outputs that were predicted by models 1 to 4, respectively. All data were divided into three data sets for training, validation, and testing. Architectures of ANNs that were used in this research are presented in Figure 11. An ANN toolbox of MATLAB Version 9.1 (R2016b) [88] was used for neural network analyses.
Figure 12 shows the training state for the ANN models for predicting the outputs containing Id2, SRH, UPV, and UCS established in the MATLAB environment. As can be seen from this figure, in models 1 to 4 after the epoch numbers of 4, 2, 3, and 2, the errors occur three times, and the test is stopped at different epoch numbers. Since the errors happen three times before stopping the test, the validation checks would be equal to three. Figure 13 shows the validation curves in ANN models for predicting the outputs containing Id2, SRH, UPV, and UCS, which are established in the MATLAB (v 9.1) environment. Based on Figure 13a, which shows the predicted Id2, the best validation performance is happening at epoch 4, and the test is stopped at epoch 7 after three error repetitions. In models 2 and 4, for predicting SRH and UCS (Figure 13b,d), the test is stopped at epoch 2 after three error repetitions. Figure 13c shows the best validation performance is happening at epoch 3 (in model 3), and the test is stopped at epoch 6 after three error repetitions. The values of mean squared error (MSE), which is found on the y-axis in Figure 13, are approximately similar in four models, and in model 1 is less than other models. The plots of regression for training, testing, and validation of the proposed ANN model were established in the MATLAB (v. 9.1) environment. Very accurate correlation coefficients were obtained, as illustrated in Figure 14. In all models, R higher than 0.99 shows the high-performance validity of established models by ANN.
To assess the performance of the results of the ANN models, the statistical parameter was calculated and is presented in Table 10. As can be seen, R values in all models are equal to 0.99. RMSE values, except in model 4 which predicted UCS, are very low, and in the predictive model of UPV is equal to 0.06. The performance index of the UPV model is equal to 99.94, which is high. Despite high values of RMSE and MAPE in model 4, the VAF value is higher than other models. The significance level in all models is equal to 0.000, indicating that predictive models by ANN estimated the output parameters with a 99% confidence limit.

5. Discussion

Statistical parameters, namely R, RMSE, VAF, MAPE, and PI, were evaluated for the developed models by SRA, MRA, and ANN to assess their accuracy. High values of the determination coefficient obtained from the SRA method indicate TC is an able coefficient for engineering properties estimation of the studied rocks. In MRA and ANN models, γdry, n, and TC were considered as inputs and four models were created for estimating Id2, SRH, UPV, and UCS of the studied samples. Then, the statistical parameter assesses the results of the MR and ANN models. A comparative diagram of R, RMSE, VAF, MAPE, and PI for the results of MRA and ANN models is presented in Figure 15. It is clear that the ANN models have better performance than the MR models. Except in model 4, for predicting UCS that shows approximately similar results for MRA and ANN models, the other models have valuable differences in statistical parameters. A higher R value for the ANN models means that there is a better relationship between experimental and calculated values in comparison to the MRA model. This result is confirmed by the research of Alber and Kahraman [65], Manouchehrian et al. [42], and Ersoy and Acar [66].
Table 11 presents the experimental and calculated values of Id2, SRH, UPV, and UCS provided by the laboratory tests and models of MR and ANN together with the obtained standard deviation (S.D.) for all models. The results show that the values of the mean and standard deviation in ANN models are approximately equal to experimental values, which indicates less variety than MR models.
The experimental and calculated values of Id2, SRH, UPV, and UCS have been plotted in diagrams of Figure 15 to compare the results of MRA and ANN for predicting the engineering characteristics of the studied rocks. The adaption rate of trend lines of MRA and ANN with the line y = x show the validity of the models. The correlations between experimental and calculated values of Id2 that were obtained from γdry, n, and TC are shown in Figure 16a. Based on these graphs, the trend line of ANN completely fits the line y = x, whereas the MRA trend line approximately intercepts the line y = x. This shows the moderate validity of the equations of the MRA for estimating Id2. The correlations between experimental and calculated values of SRH are shown in Figure 16b. As can be seen, the MRA and ANN trend lines approximately overlap together and show that both methods in this model obtained the same results for predicting Schmidt rebound hardness. Figure 16c shows the relationships between experimental results and calculated values of UPV. It is obvious that ANN estimated values of UPV are closer to the corresponding calculated values than MRA models. The trend line of the MRA model does not completely fit the 45° line and the results of the ANN are better than those. The correlations between experimental and calculated values of UCS are shown in Figure 16d. It shows that MRA and ANN trend lines of these diagrams entirely overlap together. So, both methods in this model have the same results for estimating the uniaxial compressive strength. Therefore, based on the results of different models, the predicted values from ANN models highly fit 45° lines and show the high validity of this method.
For the final assessment, residual error values were calculated for the results of MR and ANN models (Table 12). These values were calculated by subtraction of experimental and calculated values based on MRA and ANN created models. To compare the obtained values of residual error, their graphs were depicted for Id2, SRH, UPV, and UCS models, respectively (Figure 17). Very low values of S.D. of ANN models show the validity of this method for predicting engineering characteristics of the studied rocks. Additionally, it was found that the average values of residual error in ANN models had a tendency to be 0.00, for instance, in model 3. Based on Figure 17, the ANN model has a lower residual error than the MR model. It can be seen that the error trend lines of ANN are more fit to the horizontal line (x-axis).

6. Conclusions

In this research, various limestone and sandstone rock samples were investigated to develop predictive models between their engineering properties including dry unit weight (γdry), saturated unit weight (γsat), water absorption (Wa), porosity (n), Schmidt rebound hardness (SRH), ultrasonic P-wave velocity (UPV), and uniaxial compressive strength (UCS). The texture investigations carried out by using JMicroVision (v1. 27) software led to texture coefficients (TC). Predictive models were developed by statistical and soft computing methods including simple regression analysis (SRA), multiple regression analysis (MRA), and artificial neural network (ANN). Ten experimental equations were developed by the SRA for predicting γd, γs, n, Wa, SDI, SRH, UPV, and UCS from the TC values with very high values of Pearson regression coefficient (R). It was found that the TC has direct relations with γd, γs, SDI, SRH, UPV, and UCS, and inverse relations with n and Wa. Additionally, four predictive models were developed by the MRA and ANN methods where γd, n, and TC were considered as inputs and Id2, SRH, UPV, and UCS were considered as outputs. The values of the obtained correlation coefficients were very high from the MRA. In the ANN, the feed-forward neural network was used to establish predictive models, and the R value was obtained equal to 0.99 in all models. Statistical parameters including the RMSE, VAF, MAPE, and PI were evaluated and compared to each other for assessing the performance of the MRA and ANN models. The calculated statistical parameters indicated that the ANN models have better performance than the MR models. The experimental and calculated values of Id2, SRH, UPV, and UCS obtained from the laboratory tests and predicted by the MRA and ANN models were approximately equal to each other. Additionally, the results of residual error analysis show lower standard deviations for predicted data by the ANN models than the MRA ones. The results of this research indicated texture coefficient analysis is a reliable method (as mentioned in Ersoy and Waller [62], Singh et al. [36], Alber and Kahraman [65], Gupta and Sharma [63], Manouchehrian et al. [42], Tandon and Gupta [64], Akram et al. [68], Kolay and Baser [69]) to predict engineering properties of the studied rocks. Therefore, by using the obtained results from the research and statistical methods as well as artificial neural networks, it is successfully passible to establish credible models for predicting engineering properties of other rocks. This can be an advantage of the present research. On the other hand, indirect estimation of engineering properties of rocks is not always as accurate as direct laboratory methods and this can be a limitation for this research and similar research.

Author Contributions

Investigation, D.F.; supervision, L.S.; writing—original draft, D.F.; writing—review & editing, L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Fundação para a Ciência e a Tecnologia in the frame of the UIDB/00073/2020 and UIDP/00073/2020 projects of the I&D unit Geosciences Center (CGEO).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the official support of the Engineering Geology and Rock Mechanics Laboratory of Damghan University for performing all laboratory tests of the research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carvalho, J.M.F.; Lisboa, J.V.; Moura, A.C.; Carvalho, C.; Sousa, L.M.O.; Leite, M.M. Evaluation of the Portuguese ornamental stone resources. Key Eng. Mater. 2013, 548, 3–9. [Google Scholar] [CrossRef] [Green Version]
  2. Siegesmund, S.; Sousa, L.; López-Doncel, R.A. Editorial to the topical collection in Environmental Earth Sciences “Stone in the architectural heritage: From quarry to monuments-environment, exploitation, properties and durability”. Environ. Earth Sci. 2018, 77, 730. [Google Scholar] [CrossRef] [Green Version]
  3. Siegesmund, S.; Dürrast, H. Physical and mechanical properties of the rocks. In Stone in Architecture. Properties, Durability, 5th ed.; Siegesmund., S., Snethlage, R., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 97–224. [Google Scholar]
  4. Mustafa, S.; Khan, M.A.; Kha, M.R.; Sousa, L.M.O.; Hameed, F.; Mughal, M.S.; Niaz, A. Building stone evaluation—A case study of the sub-Himalayas, Muzafarabad region, Azad Kashmir, Pakistan. Eng. Geol. 2016, 209, 56–69. [Google Scholar] [CrossRef]
  5. Santos, I.; Sousa, L.; Lourenço, J. Granite resources evaluation-example of an extraction area in North of Portugal. Environ. Earth Sci. 2018, 77, 608. [Google Scholar] [CrossRef]
  6. Yarahmadi, R.; Bagherpour, R.; Taherian, S.-G.; Sousa, L.M.O. A new quality factor for the building stone industry: A case study of stone blocks, slabs, and tiles. Bull. Eng. Geol. Environ. 2019, 78, 533–542. [Google Scholar] [CrossRef]
  7. Ahmed. I.; Basharat, M.; Sousa, L.; Mughal, M.S. Evaluation of building and dimension stone using physico-mechanical and petrographic properties: A case study from the Kohistan and Ladakh batholith, Northern Pakistan. Environ. Earth Sci. 2021, 80, 759. [Google Scholar] [CrossRef]
  8. Bogdanowitsch, M.; Sousa, L.; Siegesmund, S. Building stone quarries: Resource evaluation by block modelling and unmanned aerial photogrammetric surveys. Environ. Earth Sci. 2022, 81, 16. [Google Scholar] [CrossRef]
  9. Sousa, L.M.O.; Gonçalves, B.M.M. Differences in the quality of polishing between sound and weathered granites. Environ. Earth Sci. 2013, 69, 1347–1359. [Google Scholar] [CrossRef]
  10. Sousa, L.; Siegesmund, S.; Wedekind, W. Salt weathering in granitoids: An overview on the controlling factors. Environ. Earth Sci. 2018, 77, 502. [Google Scholar] [CrossRef]
  11. Menningen, J.; Siegesmund, S.; Lopes, L.; Martins, R.; Sousa, L. The Estremoz marbles: An updated summary on the geological, mineralogical and rock physical characteristics. Environ. Earth Sci. 2018, 77, 191. [Google Scholar] [CrossRef]
  12. Vazquez, P.; Sánchez-Delgado, N.; Carrizo, L.; Thomachot-Schneider, C.; Alonso, F.J. Statistical approach of the influence of petrography in mechanical properties and durability of granitic stone. Environ. Earth Sci. 2018, 77, 287. [Google Scholar] [CrossRef]
  13. Sousa, L. Behaviour of hard stones submitted to different foot traffic. Environ. Earth Sci. 2019, 78, 680. [Google Scholar] [CrossRef]
  14. Freire-Lista, D.M.; Sousa, L.; Carter, R.; Al-Na‘īmī, F. Petrography and petrophysical characterisation of the heritage stones of Fuwairit Archaeological Site (NW Qatar) and their historical quarries: Implications for heritage conservation. Episodes 2021, 44, 43–58. [Google Scholar] [CrossRef]
  15. ISRM. The Blue Book: The Complete ISRM Suggested Methods for Rock Characterization, Testing and Monitoring, 1974–2006; Compilation Arranged by the ISRM Turkish National Group, Ankara, Turkey; Ulusay, R., Hudson, J.A., Eds.; Kazan Offset Press: Ankara, Turkey, 2007. [Google Scholar]
  16. AbuShanab, W.S.; Abd Elaziz, M.; Ghandourah, E.I.; Moustafa, E.B.; Elsheikh, A.H. A new fine-tuned random vector functional link model using Hunger games search optimizer for modeling friction stir welding process of polymeric materials. J. Mater. Res. Technol. 2021, 14, 1482–1493. [Google Scholar] [CrossRef]
  17. Elsheikh, A.H.; Abd Elaziz, M.; Vendan, A. Modeling ultrasonic welding of polymers using an optimized artificial intelligence model using a gradient-based optimizer. Weld. World 2022, 66, 27–44. [Google Scholar] [CrossRef]
  18. Elsheikh, A.H.; Muthuramalingam, T.; Shanmugan, S.; Ibrahim, A.M.; Ramesh, B.; Khoshaim, A.B.; Moustafa, E.B.; Bedairi, B.; Panchal, H.; Sathyamurthy, R. Fine-tuned artificial intelligence model using pigeon optimizer for prediction of residual stresses during turning of Inconel 718. J. Mater. Res. Technol. 2021, 15, 3622–3634. [Google Scholar] [CrossRef]
  19. Elsheikh, A.H.; Panchal, H.; Ahmadein, M.; Mosleh, A.O.; Sadasivuni, K.K.; Alsaleh, N.A. Productivity forecasting of solar distiller integrated with evacuated tubes and external condenser using artificial intelligence model and moth-flame optimizer. Case Stud. Therm. Eng. 2021, 28, 101671. [Google Scholar] [CrossRef]
  20. Khoshaim, A.B.; Moustafa, E.B.; Bafakeeh, O.T.; Elsheikh, A.H. An optimized multilayer perceptrons model using grey wolf optimizer to predict mechanical and microstructural properties of friction stir processed aluminum alloy reinforced by nanoparticles. Coatings 2021, 11, 1476. [Google Scholar] [CrossRef]
  21. Thangaraj, M.; Ahmadein, M.; Alsaleh, N.A.; Elsheikh, A.H. Optimization of abrasive water jet machining of SiC reinforced aluminum alloy based metal matrix composites using Taguchi–DEAR technique. Materials 2021, 14, 6250. [Google Scholar] [CrossRef]
  22. Elsheikh, A.H.; Abd Elaziz, M.; Ramesh, B.; Egiza, M.; Al-qaness, M.A. Modeling of drilling process of GFRP composite using a hybrid random vector functional link network/parasitism-predation algorithm. J. Mater. Res. Technol. 2021, 14, 298–311. [Google Scholar] [CrossRef]
  23. Elsheikh, A.H.; Abd Elaziz, M.; Das, S.R.; Muthuramalingam, T.; Lu, S. A new optimized predictive model based on political optimizer for eco-friendly MQL-turning of AISI 4340 alloy with nano-lubricants. J. Manuf. Process. 2021, 67, 562–578. [Google Scholar] [CrossRef]
  24. Moustafa, E.B.; Hammad, A.H.; Elsheikh, A.H. A new optimized artificial neural network model to predict thermal efficiency and water yield of tubular solar still. Case Stud. Therm. Eng. 2022, 30, 101750. [Google Scholar] [CrossRef]
  25. Abdolrasol, M.G.; Hussain, S.S.; Ustun, T.S.; Sarker, M.R.; Hannan, M.A.; Mohamed, R.; Ali, J.A.; Mekhilef, S.; Milad, A. Artificial neural networks based optimization techniques: A review. Electronics 2021, 10, 2689. [Google Scholar] [CrossRef]
  26. Deere, D.U.; Miller, R.P. Engineering Classification and Index Properties for Intact Rock; Air Force Weapons Lab Tech Report AFWL-TR 65–116, Kirtland Base; US Air Force Weapons Laboratory, Kirtland AFB: Albuquerque, NM, USA, 1966. [Google Scholar]
  27. Bell, F.G. Physical and mechanical properties of the Fell sandstones, Northumberland, England. Eng. Geol. 1978, 12, 1–29. [Google Scholar] [CrossRef]
  28. Coggan, J.S.; Stead, D.; Howe, J.H.; Faulks, C.I. Mineralogical controls on the engineering behavior of hydrothermally altered granites under uniaxial compression. Eng. Geol. 2013, 160, 89–102. [Google Scholar] [CrossRef]
  29. Alikarami, R.; Torabi, A.; Kolyukhin, D.; Skurtveit, E. Geostatistical relationships between mechanical and petrophysical properties of deformed sandstone. Int. J. Rock Mech. Min. Sci. 2013, 63, 27–38. [Google Scholar] [CrossRef]
  30. Cantisani, E.; Garzonio, C.A.; Ricci, M.; Vettori, S. Relationships between the petrographical, physical and mechanical properties of some Italian sandstones. Int. J. Rock Mech. Min. Sci. 2013, 60, 321–332. [Google Scholar] [CrossRef]
  31. Abdlmutalib, A.; Abdullatif, O.; Korvin, G.; Abdulraheem, A. The relationship between lithological and geomechanical properties of tight carbonate rocks from Upper Jubaila and Arab-D Member outcrop analog, Central Saudi Arabia. Arab. J. Geosci. 2015, 8, 11031–11048. [Google Scholar] [CrossRef]
  32. Koralegedara., N.H.; Maynard, J.B. Chemical, mineralogical and textural properties of the Kope Formation mudstones: How they affect its durability. Eng. Geol. 2017, 228, 312–322. [Google Scholar] [CrossRef] [Green Version]
  33. Cowie, S.; Walton, G. The effect of mineralogical parameters on the mechanical properties of granitic rocks. Eng. Geol. 2018, 240, 204–225. [Google Scholar] [CrossRef]
  34. Fereidooni, D.; Khajevand, R. Determining the geotechnical characteristics of some sedimentary rocks from Iran with an emphasis on the correlations between physical index and mechanical properties. Geotech. Test. J. 2018, 41, 555–573. [Google Scholar] [CrossRef]
  35. Lawal, A.I.; Oniyide, G.O.; Kwon, S.; Onifade, M.; Köken, E.; Ogunsola, N.O. Prediction of Mechanical Properties of Coal from Non-destructive Properties: A Comparative Application of MARS, ANN, and GA. Nat. Resour. Res. 2021, 30, 4547–4563. [Google Scholar] [CrossRef]
  36. Singh, V.K.; Singh, D.; Singh, T.N. Prediction of strength properties of some schistose rocks from petrographic properties using artificial neural networks. Int. J. Rock Mech. Min. Sci. 2001, 38, 269–284. [Google Scholar] [CrossRef]
  37. Kilic, A.; Teymen, A. Determination of mechanical properties of rocks using simple methods. Bull. Eng. Geol. Environ. 2008, 67, 237–244. [Google Scholar] [CrossRef]
  38. Jensen, L.R.; Friis, H.; Fundal, E.; Moller, P.; Jespersen, M. Analysis of limestone micromechanical properties by optical microscopy. Eng. Geol. 2010, 110, 43–50. [Google Scholar] [CrossRef]
  39. Asadi, M.; Eftekhari, M.; Bagheripour, M.H. Evaluating the strength of intact rocks through genetic programming. Appl. Soft Comput. 2011, 11, 1932–1937. [Google Scholar] [CrossRef]
  40. Cevik., A.; Sezer, E.A.; Cabalar, A.F.; Gokceoglu, C. Modeling of the uniaxial compressive strength of some clay-bearing rocks using neural network. Appl. Soft Comput. 2011, 11, 2587–2594. [Google Scholar] [CrossRef]
  41. Singh, T.N.; Verma, A.K. Comparative analysis of intelligent algorithms to correlate strength and petrographic properties of some schistose rocks. Eng. Comput. 2012, 28, 1–12. [Google Scholar] [CrossRef]
  42. Manouchehrian, A.; Sharifzadeh, M.; Moghadam, R.H. Application of artificial neural networks and multivariate statistics to estimate UCS using textural characteristics. Int. J. Min. Sci. Technol. 2012, 22, 229–236. [Google Scholar] [CrossRef]
  43. 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]
  44. Esamaldeen, A.; Guang, W. Selection of influential microfabric properties of anisotropic amphibolite rocks on its uniaxial compressive strength (UCS): A comprehensive statistical study. J. Appl. Math. Phys. 2014, 2, 1130–1138. [Google Scholar]
  45. Liu, Z.; Shao, J.; Xu, W.; Zhang, Y.; Chen, H. Prediction of elastic compressibility of rock material with soft computing techniques. Appl. Soft Comput. 2014, 22, 118–125. [Google Scholar] [CrossRef]
  46. Chen, W.; Konietzky, H.; Abbas, S.M. Numerical simulation of time-independent and -dependent fracturing in sandstone. Eng. Geol. 2015, 193, 118–131. [Google Scholar] [CrossRef]
  47. Ajalloeian, R.; Mansouri, H.; Baradaran, E. Some carbonate rock texture effects on mechanical behavior, based on Koohrang tunnel data, Iran. Bull. Eng. Geol. Environ. 2016, 76, 295–307. [Google Scholar] [CrossRef]
  48. Fereidooni, D. Determination of the geotechnical characteristics of hornfelsic rocks with a particular emphasis on the correlation between physical and mechanical properties. Rock Mech. Rock Eng. 2016, 49, 2595–2608. [Google Scholar] [CrossRef]
  49. Ahmad, M.; Ansari, M.K.; Singh, R.; Sharma, L.K.; Singh, T.N. Assessment of durability and weathering state of some igneous and metamorphic rocks using micropetrographic index and rock durability indicators: A case study. Geotech. Geol. Eng. 2017, 35, 827–842. [Google Scholar] [CrossRef]
  50. Germinario, L.; Siegesmund, S.; Maritan, L.; Mazzoli, C. Petrophysical and mechanical properties of Euganean trachyte and implications for dimension stone decay and durability performance. Environ. Earth Sci. 2017, 76, 739. [Google Scholar] [CrossRef]
  51. Yalcinalp, B.; Aydin, Z.O.; Ersoy, H.; Seren, A. Investigation of geological, geotechnical and geophysical properties of Kiratli (Bayburt, NE Turkey) travertine. Carbonates Evaporites 2018, 33, 421–429. [Google Scholar] [CrossRef]
  52. Matin, S.S.; Farahzadi, L.; Makaremi, S.; Chelgani, S.C.; Sattari, G.H. Variable selection and prediction of uniaxial compressive strength and modulus of elasticity by random forest. Appl. Soft Comput. 2018, 70, 980–987. [Google Scholar] [CrossRef]
  53. Brace, W.F. Dependence of fracture strength of rocks on grain size. In 4th US Symposium on Rock Mechanics (USRMS); OnePetro: University Park, PA, USA, 1961; pp. 99–103. [Google Scholar]
  54. Das, S.K.; Kumar, A.; Das, B.; Burnwal, A.P. On soft computing techniques in various areas. Comput. Sci. Inf. Technol 2013, 3, 166. [Google Scholar]
  55. Ulusay, R.; Tureli, K.; Ider, M.H. Prediction of engineering properties of selected litharenite sandstone from its petrographic characteristics using correlation and multivariate statistical techniques. Eng. Geol. 1994, 37, 135–157. [Google Scholar] [CrossRef]
  56. Tugrul, A.; Zarif, I.H. Correlation of mineralogical and textural characteristics with engineering properties of selected granitic rocks from Turkey. Eng. Geol. 1999, 51, 303–317. [Google Scholar] [CrossRef]
  57. Eberli, G.P.; Baechle, G.T.; Anselmetti, F.S.; Incze, M.L. Factors controlling elastic properties in carbonate sediments and rocks. Lead. Edge 2003, 22, 654–660. [Google Scholar] [CrossRef]
  58. Meng, J.; Pan, J. Correlation between petrographic characteristics and failure duration in clastic rocks. Eng. Geol. 2007, 89, 258–265. [Google Scholar] [CrossRef]
  59. Khanlari, G.R.; Heidari, M.; Noori, M.; Momeni, A. The effect of petrographic characteristics on engineering properties of conglomerates from famenin region, northeast of Hamedan, Iran. Rock Mech. Rock Eng. 2016, 49, 2609–2621. [Google Scholar] [CrossRef]
  60. Howarth, D.F.; Rowlands, J.C. Development of an index to quantify rock texture for qualitative assessment of intact rock properties. Geotech. Test. J. 1986, 9, 169–179. [Google Scholar]
  61. Howarth, D.F.; Rowlands, J.C. Quantitative assessment of rock texture and correlation with drillability and strength properties. Rock Mech. Rock Eng. 1987, 20, 57–85. [Google Scholar] [CrossRef]
  62. Ersoy, A.; Waller, M.D. Textural characterization of rocks. Eng. Geol. 1995, 39, 123–136. [Google Scholar] [CrossRef]
  63. Gupta, V.; Sharma, R. Relationship between textural, petrophysical and mechanical properties of quartzites: A case study from northwestern Himalaya. Eng. Geol. 2012, 135, 1–9. [Google Scholar] [CrossRef]
  64. Tandon, R.S.; Gupta, V. The control of mineral constituents and textural characteristics on the petrophysical and mechanical (PM) properties of different rocks of the Himalaya. Eng. Geol. 2013, 153, 125–143. [Google Scholar] [CrossRef]
  65. Alber, M.; Kahraman, S. Predicting the uniaxial compressive strength and elastic modulus of a fault breccia from texture coefficient. Rock Mech. Rock Eng. 2009, 42, 117–127. [Google Scholar] [CrossRef]
  66. Ersoy, H.; Acar, S. Influences of petrographic and textural properties on the strength of very strong granitic rocks. Environ. Earth Sci. 2016, 75, 1461–1476. [Google Scholar] [CrossRef]
  67. Aligholi, S.; Lashkaripour, G.R.; Ghafoori, M.; Azali, S.T. Evaluating the relationships between NTNU/SINTEF drillability indices with index properties and petrographic data of hard igneous rocks. Rock Mech. Rock Eng. 2017, 50, 2929–2953. [Google Scholar] [CrossRef]
  68. Akram, M.S.; Farooq, S.; Naeem, M.; Ghazi, G. Prediction of mechanical behaviour from mineralogical composition of Sakesar limestone, Central Salt Range, Pakistan. Bull. Eng. Geol. Environ. 2017, 76, 601–615. [Google Scholar] [CrossRef]
  69. Kolay, E.; Baser, T. The effect of the textural characteristics on the engineering properties of the basalts from Yozgat region, Turkey. J. Geol. Soc. India 2017, 90, 102–110. [Google Scholar] [CrossRef]
  70. Geological Society of Iran (GSI). Geological Quadrangle Map of Iran. No. D6, Scale 1:100,000; Offset Press: Tehran, Iran, 1977. [Google Scholar]
  71. ASTM. Standard Guide for Petrographic Examination of Dimension Stone (C1721); Book Standards vol 04.07; ASTM International: West Conshohocken, PA, USA, 2009. [Google Scholar] [CrossRef]
  72. ASTM. Standard test method for slake-durability of shales and similar weak rocks (D-4644). In Annual Book of ASTM Standards; ASTM International: West Conshohocken, PA, USA, 1990; Volume 4.08, pp. 863–865. [Google Scholar]
  73. ISRM. Suggested Methods for Determining Hardness and Abrasiveness of Rocks; Commission on Standardization of Laboratory and Field Test; International Society for Rock Mechanics: Salzburg, Austria, 1978; Volume 15, pp. 89–97. [Google Scholar]
  74. ASTM. Standard Test Method for Determination of Rock Hardness by Rebound Hammer Method; ASTM standards on disc 04.09; ASTM International: West Conshohocken, PA, USA, 2001; pp. D5873–D5880. [Google Scholar]
  75. ASTM. Standard Test Method for Laboratory Determination of Pulse Velocities and Ultrasonic Elastic Constants of Rock; ASTM International: West Conshohocken, PA, USA, 1996; pp. D2845–D2895. [Google Scholar]
  76. ASTM. Standard Test Method for Unconfined Compressive Strength of Intact Rock Core Specimens; ASTM standards on disc 04.08; ASTM International: West Conshohocken, PA, USA, 1995; p. D2938. [Google Scholar]
  77. IAEG. Classification of rocks and soils for engineering geological mapping, Part 1: Rock and soil materials. Bull. Int. Assoc. Eng. Geol. 1979, 19, 364–371. [Google Scholar] [CrossRef]
  78. Gamble, J.C. Durability-Plasticity Classification of Shales and Other Argillaceous Rocks. Ph.D. Thesis, University of Illinois, Urbana-Champaign, IL, USA, 1971; p. 161. [Google Scholar]
  79. Broch, E.; Franklin, J.A. The point load strength test. Int. J. Rock Mech. Min. Sci. 1972, 9, 669–697. [Google Scholar] [CrossRef]
  80. Sousa, L.; Menningen, J.; López-Doncel, R.; Siegesmund, S. Petrophysical properties of limestones: Influence on behaviour under different environmental conditions and applications. Environ. Earth Sci. 2021, 80, 814. [Google Scholar] [CrossRef]
  81. Williams, H.; Turner, F.J.; Gilber, C.M. Petrography: An Introduction to the Study of Rocks in Thin Section; W.H. Freeman Company: San Francisco, CA, USA, 1954; p. 406. [Google Scholar]
  82. Prikiryl, R. Assessment of rock geomechanical quality by quantitative rock fabric coefficients: Limitations and possible source of misinterpretations. Eng. Geol. 2006, 87, 149–162. [Google Scholar] [CrossRef]
  83. IBM Corp. Released. IBP SPSS Statistics for Windows; Version 24.0; IBM Corp: Armonk, NY, USA, 2016. [Google Scholar]
  84. Fereidooni, D. Importance of the mineralogical and textural characteristics in the mechanical properties of rocks. Arab. J. Geosci. 2022, 15, 637. [Google Scholar] [CrossRef]
  85. Singh, R.; Kainthola, A.; Singh, T.N. Estimation of elastic constant of rocks using an ANFIS approach. Appl. Soft Comput. 2012, 12, 40–45. [Google Scholar] [CrossRef]
  86. Sirdesai, N.N.; Singh, A.; Sharma, L.K.; Singh, R.; Singh, T.N. Development of novel methods to predict the strength properties of thermally treated sandstone using statistical and soft-computing approach. Neural Comput. Appl. 2017, 31, 2841–2867. [Google Scholar] [CrossRef]
  87. Suparta, W.; Alhasa, K.M. Modeling of Tropospheric Delays Using ANFIS. Springer Briefs in Meteorology; Springer: Cham, Switzerland, 2016. [Google Scholar]
  88. MATLAB and Statistical Toolbox; The Mathworks, Inc.: Natick, MA, USA, 2016.
Figure 1. Sampling locations and geological map of the study area [70].
Figure 1. Sampling locations and geological map of the study area [70].
Materials 15 07922 g001
Figure 2. Test plan in this study for determining different properties of the rocks.
Figure 2. Test plan in this study for determining different properties of the rocks.
Materials 15 07922 g002
Figure 3. Block samples in the field and core specimens of the rocks.
Figure 3. Block samples in the field and core specimens of the rocks.
Materials 15 07922 g003
Figure 4. Testing procedure and the used apparatuses for the study. (a) P-wave velocity test, (b) uniaxial compressive test, and (c) slake durability test.
Figure 4. Testing procedure and the used apparatuses for the study. (a) P-wave velocity test, (b) uniaxial compressive test, and (c) slake durability test.
Materials 15 07922 g004
Figure 5. Microscopic images of the studied rock samples. (Qtz.: Quartz (SiO2); Cal.: Calcite (CaCO3); Dol.: Dolomite (Ca.Mg(CO3)2); Fld.: Feldspar; Fos.: Fossil).
Figure 5. Microscopic images of the studied rock samples. (Qtz.: Quartz (SiO2); Cal.: Calcite (CaCO3); Dol.: Dolomite (Ca.Mg(CO3)2); Fld.: Feldspar; Fos.: Fossil).
Materials 15 07922 g005
Figure 6. XRD graphs of the studied rock samples.
Figure 6. XRD graphs of the studied rock samples.
Materials 15 07922 g006
Figure 7. JMicroVision (v1. 27) software environment for measuring TC of the rocks.
Figure 7. JMicroVision (v1. 27) software environment for measuring TC of the rocks.
Materials 15 07922 g007
Figure 8. (a) Diagram of the texture coefficient (TC) values and (b) histogram of the TC for the studied rock samples.
Figure 8. (a) Diagram of the texture coefficient (TC) values and (b) histogram of the TC for the studied rock samples.
Materials 15 07922 g008
Figure 9. Relations between texture coefficient (TC) and (a) dry and saturated unit weights (γd, γs), (b) effective porosity (n) and water absorption (Wa), (c) slake durability index (SDI), (d) Schmidt rebound hardness (SRH), (e) P-wave velocity (UPV), (f) uniaxial compressive strength (UCS) for the rocks.
Figure 9. Relations between texture coefficient (TC) and (a) dry and saturated unit weights (γd, γs), (b) effective porosity (n) and water absorption (Wa), (c) slake durability index (SDI), (d) Schmidt rebound hardness (SRH), (e) P-wave velocity (UPV), (f) uniaxial compressive strength (UCS) for the rocks.
Materials 15 07922 g009
Figure 10. Three-dimensional view of the correlations among (a) texture coefficient (TC), slake durability index (Id2) and water absorption (Wa); (b) texture coefficient (TC), Schmidt rebound hardness (SRH) and porosity; (c) texture coefficient (TC), uniaxial compressive strength (UCS), and Schmidt rebound hardness (SRH).
Figure 10. Three-dimensional view of the correlations among (a) texture coefficient (TC), slake durability index (Id2) and water absorption (Wa); (b) texture coefficient (TC), Schmidt rebound hardness (SRH) and porosity; (c) texture coefficient (TC), uniaxial compressive strength (UCS), and Schmidt rebound hardness (SRH).
Materials 15 07922 g010
Figure 11. Schematic representation of the constructed neural network.
Figure 11. Schematic representation of the constructed neural network.
Materials 15 07922 g011
Figure 12. Training state in ANN models for prediction of (a) Id2, (b) SRH, (c) UPV, and (d) UCS.
Figure 12. Training state in ANN models for prediction of (a) Id2, (b) SRH, (c) UPV, and (d) UCS.
Materials 15 07922 g012
Figure 13. Validation curves in ANN models for prediction of (a) Id2, (b) SRH, (c) UPV, and (d) UCS.
Figure 13. Validation curves in ANN models for prediction of (a) Id2, (b) SRH, (c) UPV, and (d) UCS.
Materials 15 07922 g013
Figure 14. Regression plots of ANN models for predicting (a) Id2, (b) SRH, (c) UPV, and (d) UCS.
Figure 14. Regression plots of ANN models for predicting (a) Id2, (b) SRH, (c) UPV, and (d) UCS.
Materials 15 07922 g014
Figure 15. Comparative diagrams of statistical coefficients for results of MRA and ANN models (a) R, (b) RMSE, (c) VAF, (d) MAPE, and (e) PI.
Figure 15. Comparative diagrams of statistical coefficients for results of MRA and ANN models (a) R, (b) RMSE, (c) VAF, (d) MAPE, and (e) PI.
Materials 15 07922 g015
Figure 16. Plots of experimental and calculated values of (a) Id2, (b) SRH, (c) UPV, and (d) UCS from the MRA and ANN.
Figure 16. Plots of experimental and calculated values of (a) Id2, (b) SRH, (c) UPV, and (d) UCS from the MRA and ANN.
Materials 15 07922 g016
Figure 17. Residual error plots of MRA and ANN models to predict (a) Id2, (b) SRH, (c) UPV, and (d) UCS.
Figure 17. Residual error plots of MRA and ANN models to predict (a) Id2, (b) SRH, (c) UPV, and (d) UCS.
Materials 15 07922 g017
Table 1. Mineralogical and petrographic properties of the studied rock samples.
Table 1. Mineralogical and petrographic properties of the studied rock samples.
SampleLithologyMinerals (%)
Qtz.Cal.Dol.Fld.Fos.Other Minerals
SBZ1Dolomitic limestone106025--5
SBZ2Dolomitic limestone204035--5
ELK1Ooid grainstone (limestone)510--805
ELK2Ooid grainstone (limestone)520--705
DLC1Sandy grainstone (limestone)20655--10
DLC2Sandy grainstone (limestone)25605--10
LAR1Microcrystalline limestone1570--123
LAR2Microcrystalline limestone2063--152
PDH1Greywacke (sandstone)701510--5
PDH2Greywacke (sandstone)553010--5
PDH3Greywacke (sandstone)553010--5
BRT1Litharenite (sandstone)7010-15-5
BRT2Litharenite (sandstone)6015-20-5
SMK1Sub litharenite (sandstone)707-10-13
SMK2Sub litharenite (sandstone)807-5-8
Note: Qtz.: Quartz (SiO2); Cal.: Calcite (CaCO3); Dol.: Dolomite (Ca.Mg(CO3)2); Fld.: Feldspar; Fos.: Fossil.
Table 2. Values of engineering properties of the tested rock samples.
Table 2. Values of engineering properties of the tested rock samples.
Sampleγdry (g/cm3)γsat (g/cm3)n (%)Wa (%)Id1
(%)
Id2
(%)
Id3
(%)
SRHUPV (Km/s)UCS
(MPa)
SBZ12.452.533.741.1299.0298.6998.46434.8981.99
SBZ22.452.534.981.7899.2598.8798.52394.8067.40
ELK12.482.558.153.4998.9198.4397.75314.5945.54
ELK22.452.527.012.8799.1098.1597.89354.6156.58
DLC12.502.616.502.4099.4198.6698.21434.9185.59
DLC22.502.627.212.9899.1598.5798.02404.3773.67
LAR12.512.733.961.2499.3798.9998.94455.4495.45
LAR22.572.702.850.6499.5699.3499.18475.75108.55
PDH12.442.537.583.6998.1097.2596.56303.8546.29
PDH22.452.528.413.4298.1997.2696.42263.5232.46
PDH32.352.4610.324.8298.0497.0495.54243.0916.09
BRT12.592.672.150.2699.5299.0598.76485.53119.82
BRT22.562.623.330.9099.0398.7197.87435.1090.12
SMK12.562.673.100.7799.6199.2999.15464.91105.02
SMK22.602.682.000.1899.6599.3899.21495.16127.43
Table 3. Descriptive results of linear regression analysis for correlating the engineering properties of the rocks.
Table 3. Descriptive results of linear regression analysis for correlating the engineering properties of the rocks.
Correlations
γdryγsatnWaId1Id2Id3SRHUPVUCS
GsR
Sig.
γdryR1
Sig.
γsatR0.8751
Sig.0.000
nR−0.845−0.7681
Sig.0.0000.001
WaR−0.855−0.7740.9941
Sig.0.0000.0010.000
Id1R0.7960.784−0.803−0.8331
Sig.0.0000.0010.0000.000
Id2R0.8200.801−0.869−0.8890.9711
Sig.0.0000.0000.0000.0000.000
Id3R0.7920.794−0.868−0.8840.9640.9781
Sig.0.0000.0000.0000.0000.0000.000
SRHR0.8490.843−0.930−0.9370.9230.9400.9321
Sig.0.0000.0000.0000.0000.0000.0000.000
UPVR0.8130.810−0.878−0.8870.9000.9240.9290.9171
Sig.0.0000.0000.0000.0000.0000.0000.0000.000
UCSR0.9140.873−0.945−0.9480.8910.9090.9020.9800.8911
Sig.0.0000.0000.0000.0000.0000.0000.0000.0000.000
Note: R: Pearson correlation coefficient and Sig.: significant level, Sig. (two-tailed).
Table 4. Texture coefficient derivations for the rocks.
Table 4. Texture coefficient derivations for the rocks.
SampleAW N 0 N 0 + N 1 1 FF 0 N 0 N 0 + N 1 AR1AF1TC
SBZ10.660.811.290.191.730.50.8
SBZ20.660.821.130.181.650.50.7
ELK10.470.921.190.081.560.50.55
ELK20.60.841.210.161.740.50.7
DLC10.80.921.060.081.550.50.83
DLC20.710.91.10.11.530.50.75
LAR10.860.811.180.191.750.50.97
LAR20.980.831.140.171.690.51.06
PDH10.570.911.030.091.460.50.57
PDH20.580.971.040.031.470.50.6
PDH30.380.961.10.041.520.50.42
BRT10.870.831.250.171.80.51.04
BRT20.770.831.160.171.850.50.86
SMK10.850.851.150.151.680.50.94
SMK20.90.821.170.181.790.51.01
Note: The meaning of the acronyms and letters in the table were previously presented.
Table 5. Simple regression models for texture coefficient (TC) of the studied rocks.
Table 5. Simple regression models for texture coefficient (TC) of the studied rocks.
ModelParameterPredictive ModelEquation TypeR2
1TC—γdγd = 2.261 × e0.1259TCExponentialR2 = 0.80
2TC—γsγs = 2.3149 × e0.1451TCExponentialR² = 0.83
3TC—nn = −12.328TC + 15.118LinearR² = 0.86
4TC—WaWa = −6.8834TC + 7.4529LinearR² = 0.88
5TC—Id1Id1 = 99.541TC0.0178PowerR² = 0.77
6TC—Id2Id2 = 99.196TC0.0255PowerR² = 0.78
7TC—Id3Id3 = 99.033TC0.0375PowerR² = 0.82
8TC—HSSRH = 28.948ln(TC) + 47.149LogarithmicR² = 0.92
9TC—VPUPV = 5.423TC0.5718PowerR² = 0.82
10TC—UCSUCS = 161.08TC − 49.913LinearR² = 0.94
Table 6. Statistical parameters for performance of SRA models.
Table 6. Statistical parameters for performance of SRA models.
ModelRRMSEVAF (%)MAPE (%)PISig. (Two-Tailed)
10.890.0379.732.0679.770.000
20.910.0382.472.7683.790.000
30.930.9585.9240.5285.910.000
40.940.4887.9874.0888.400.000
50.880.2575.850.1377.510.000
60.880.3578.270.0678.440.000
70.910.4581.700.2683.360.000
80.962.1792.425.3790.750.000
90.910.3082.202.3083.520.000
100.977.5694.263.7187.380.000
Table 7. Multiple regression coefficients for the developed models.
Table 7. Multiple regression coefficients for the developed models.
Variable
ModelUnstandardized CoefficientsStandardized CoefficientstSig.95.0% Confidence Interval for BCollinearity Statistics
BStd. ErrorBetaLower BoundUpper BoundToleranceVIF
1Id293.6357.982 11.7310.00076.068111.203
γd1.7433.3920.1570.5140.618−5.7239.2090.2004.991
n−0.1150.107−0.391−1.0680.308−0.3510.1220.1397.183
TC1.4571.6940.3730.8600.408−2.2725.1850.09910.114
2SRH35.30750.731 0.6960.501−76.351146.965
γd−5.28321.559−0.045−0.2450.811−52.73542.1690.2004.991
n−1.0130.682−0.325−1.4850.166−2.5150.4880.1397.183
TC28.78510.7670.6942.6730.0225.08752.4830.09910.114
3UPV2.9697.110 0.4180.684−12.68118.619
γd0.1313.0220.0120.0430.966−6.5206.7820.2004.991
n−0.0820.096−0.288−0.8530.412−0.2920.1290.1397.183
TC2.3491.5090.6241.5560.148−0.9735.6700.09910.114
4UCS−215.929138.709 −1.5570.148−521.22589.367
γd98.47458.9480.2081.6710.123−31.269228.2170.2004.991
n−3.6801.865−0.295−1.9730.074−7.7860.4250.1397.183
TC84.86029.4400.5112.8830.01520.064149.6560.09910.114
Note: γd, n, and TC are independent variables (predictors); Id2, SRH, UPV, and UCS are dependent variables.
Table 8. Various Rs values for the selected models.
Table 8. Various Rs values for the selected models.
Model Summary a
ModelRR SquareAdjusted R SquareStd. Error of the EstimateChange StatisticsDurbin–Watson
R Square ChangeF Changedf1df2Sig. F Change
10.892 b0.7950.7390.392130.79514.2203110.0000.919
20.963 c0.9270.9072.492320.92746.3973110.0001.206
30.908 d0.8250.7780.349320.82517.3153110.0000.814
40.983 e0.9660.9566.814530.966103.4213110.0001.931
a. Predictors (constants): TC, γd, n; b. Dependent variable: Id2; c. Dependent variable: SRH; d. Dependent variable: UPV; e. Dependent variable: UCS.
Table 9. Predictive models from multiple regression analysis and statistical parameters for performance of the models.
Table 9. Predictive models from multiple regression analysis and statistical parameters for performance of the models.
ModelPredictive ModelRMSEVAF (%)MAPE (%)PISig.
(Two-Tailed)
1Id2 = 93.64 + 1.74γd − 0.12n + 1.46TC0.3479.470.0780.460.000
2SRH = 35.31 − 5.28γd − 1.01n + 28.79TC2.1392.683.1991.790.000
3UPV = 2.97 + 0.13γd − 0.08n + 2.35TC0.3082.520.3483.530.000
4UCS = −215.93 + 98.47γd − 3.68n + 84.86TC5.8496.583.1092.130.000
Table 10. Statistical parameters for performance of ANN models.
Table 10. Statistical parameters for performance of ANN models.
ModelRRMSEVAF (%)MAPE (%)PISig. (Two-Tailed)
10.990.0998.600.3699.890.000
20.990.4899.630.4299.220.000
30.990.0699.390.8599.940.000
40.992.0599.672.4799.670.000
Table 11. Experimental and calculated values of Id2, SRH, UPV, and UCS obtained from the MRA and ANN models.
Table 11. Experimental and calculated values of Id2, SRH, UPV, and UCS obtained from the MRA and ANN models.
SampleModel 1 (Id2)Model 2 (SRH)Model 3 (UPV)Model 4 (UCS)
Ex.MRAANNEx.MRAANNEx.MRAANNEx.MRAANN
SBZ198.6998.6299.044341.6343.184.894.874.9381.9979.4584.02
SBZ298.8798.3398.873937.5039.144.804.544.8067.4066.467.88
ELK198.4397.7898.433129.8230.744.593.934.5745.5444.9646.14
ELK298.1598.0898.153535.4534.894.614.374.6256.5858.9356.67
DLC198.6698.4298.664339.4442.914.914.734.9385.5976.7685.86
DLC298.5798.2298.564036.4239.954.374.484.3773.6767.3674.01
LAR198.9998.9598.984545.9845.155.445.265.4895.4598.9799.09
LAR299.3499.3299.374749.3847.905.755.575.83108.55116.6111.65
PDH197.2597.8197.253031.1829.883.854.023.8446.2944.8141.78
PDH297.2697.7797.262631.1525.823.524.033.5132.4645.2933.38
PDH397.0497.1097.042424.5723.253.093.443.0716.0913.1418.22
BRT199.0599.4199.064849.4048.305.535.585.35119.82119.45120.68
BRT298.7198.9598.744343.1943.845.105.065.1390.1296.8891.46
SMK199.2999.0999.324645.7245.034.915.264.96105.09104.51107.73
SMK299.3899.4099.384948.6448.675.165.525.17127.43118.44127.97
Mean98.5198.4898.5439.2739.3039.244.704.714.7076.8076.8077.77
S.D.0.770.700.778.167.858.360.740.670.7432.6432.0833.30
Table 12. Residual error values of MR and ANN models.
Table 12. Residual error values of MR and ANN models.
SampleModel 1 (Id2)Model 2 (SRH)Model 3 (UPV)Model 4 (UCS)
MRAANNMRAANNMRAANNMRAANN
SBZ10.068−0.3531.371−0.1830.017−0.0422.544−2.027
SBZ20.5430.0031.503−0.1410.261−0.0051.003−0.479
ELK10.6500.0001.1810.2640.6590.0210.583−0.599
ELK20.0660.000−0.4470.1070.240−0.005−2.347−0.087
DLC10.238−0.0013.5590.0890.188−0.0218.831−0.269
DLC20.3500.0143.5800.049−0.115−0.0016.313−0.339
LAR10.0420.007−0.984−0.1550.182−0.042−3.521−3.639
LAR20.023−0.029−2.379−0.9000.186−0.080−8.051−3.104
PDH1−0.5580.000−1.1810.123−0.1660.0121.4774.508
PDH2−0.5100.000−5.1540.176−0.5090.010−12.829−0.919
PDH3−0.0640.000−0.5710.752−0.3440.0272.952−2.133
BRT1−0.357−0.011−1.405−0.295−0.0490.1770.370−0.857
BRT2−0.240−0.030−0.189−0.8400.045−0.032−6.758−1.339
SMK10.195−0.0290.2750.970−0.358−0.0540.576−2.645
SMK2−0.0190.0020.3600.335−0.361−0.0088.989−0.536
Mean0.03−0.03−0.030.02−0.010.000.01−0.96
S.D.0.340.092.130.480.300.065.841.81
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fereidooni, D.; Sousa, L. Predicting the Engineering Properties of Rocks from Textural Characteristics Using Some Soft Computing Approaches. Materials 2022, 15, 7922. https://doi.org/10.3390/ma15227922

AMA Style

Fereidooni D, Sousa L. Predicting the Engineering Properties of Rocks from Textural Characteristics Using Some Soft Computing Approaches. Materials. 2022; 15(22):7922. https://doi.org/10.3390/ma15227922

Chicago/Turabian Style

Fereidooni, Davood, and Luís Sousa. 2022. "Predicting the Engineering Properties of Rocks from Textural Characteristics Using Some Soft Computing Approaches" Materials 15, no. 22: 7922. https://doi.org/10.3390/ma15227922

APA Style

Fereidooni, D., & Sousa, L. (2022). Predicting the Engineering Properties of Rocks from Textural Characteristics Using Some Soft Computing Approaches. Materials, 15(22), 7922. https://doi.org/10.3390/ma15227922

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