Prediction of Uniaxial Compression Strength of Limestone Based on the Point Load Strength and SVM Model

: Uniaxial compression strength (UCS) is a fundamental parameter to carry out geotechnical engineering design and construction. It is simple and efﬁcient to predict UCS using point load strength (PLS) at engineering sites. However, the high dispersion of rock strength limits the accuracy of traditional ﬁtting prediction methods. In order to improve the UCS prediction accuracy, 30 sets of regular cylindrical specimen tests between PLS and UCS are conducted on limestone mines. The correlation relationship between PLS and UCS is found by using four basic ﬁtting functions. Then, a prediction model is established by using SVM algorithm. Multiple training test data are used to achieve high-precision prediction of UCS and the results show it is less different from the actual values. Especially, the R 2 coefﬁcient reached 0.98. The SVM model prediction performance is signiﬁcantly better than the traditional ﬁtting function. The constructed SVM model in this study can accurately predict the UCS using the PLS obtained in the ﬁeld, which has a great signiﬁcance to the rock stability judgment in the actual construction environment.


Introduction
Construction projects such as water and hydropower, transportation, mining, industrial and civil construction engineering all involve rock engineering. In order to accurately and promptly assess the stability of the mountains involved in the construction project and ensure the safety of the structure, rock classification methods are often used in the rock engineering project and design. To evaluate the stability of engineering rock mass, the rock classification divides engineering rock mass into several levels with different basic quality and strength. In this way, it provides an important basis and safety protection selection for the technical planning and various construction engineering [1].
Quantitative indicators of strength in rock classification are generally expressed in terms of uniaxial compression strength of rocks. However, the determination of the uniaxial compression strength usually needs to be measured in the laboratory, which cannot meet the requirements of obtaining data quickly at the engineering site [2]. In 1972, Broch and Franklin [3] gave the first point load strength formula by studying the damage model analysis of cylindrical specimens analyzed along the diameter direction and found that the uniaxial compression strength values of certain rocks were approximately 24 times higher than those of rock specimens tested with a 50 mm diameter point load. This discovery provided an alternative method of simple operation and high accuracy for testing the uniaxial compression strength on the engineering site. Compared with the uniaxial compression strength, testing the point load strength of the rock has the advantages of convenience and efficiency, and the strength index of the rock can be obtained in a shorter time period. This makes the point load strength very important in engineering practice. It is widely used in practical scenarios such as tunnel engineering, slope stability analysis, mining roadway support, etc [4,5]. Subsequently, many scholars and scientific research institutions have launched a large number of related research.
Bieniawski [6] stated that the uniaxial compression strength of rocks is about 23 times the point load strength. The International Society of Rock Mechanics [7] and the American Society for Testing Materials [8] proposed a general test method and procedure for point load strength and gave a general formula R = (20~25) I s (50) between them in 1985. Tsiambaos and Sabatakakis [9] studied the ratio of uniaxial compression strength to point load strength for intact sedimentary rocks. The test results showed that this ratio varied with different point load strength distributions. For example, different ratios existed from soft sedimentary rocks with point load strength less than 2 MPa to hard rocks with point load strength more than 5 MPa. A. Basu and M. Kamran [10] have explored the feasibility of point load tests to evaluate anisotropic rocks uniaxial compression strength such as schists based on successful estimation of uniaxial compression strength by point load tests on isotropic rocks. Through evaluating uniaxial compression strength of these rocks, relevant fitting formulas have been obtained with good prediction results. Robina H. C. Wong [11] investigated the correlation between point load strength and uniaxial compression strength of irregular volcanic rock masses from 16 different locations and obtained a good linear fit based on the analysis of all test data. Accordingly, they found a suitable size correction factor by least squares method.
Based on the summary of a large number of literatures, these research results are arranged in the table by year, as shown in Table 1. It can be identified that many scholars [12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27] have found the correlation between point load strength and uniaxial compression strength by experimental methods. However, the data processing of these research results is carried out by function fitting, and there are large deviations in predicting the test points that deviate from the data center. Moreover, it can be found from Table 1 that most of the derived correlation formulas are linear, but some of them are nonlinear functions such as exponential, logarithmic, and power functions. These results also imply that there exists a nonlinear relationship of limestone between the point load strength and the uniaxial compression strength [14,26]. This shows that large deviation will be produced if the general linear formulas are simply used to perform the uniaxial compression strength estimation. So linear relationship cannot be blindly extended to all engineering scenarios for undeveloped projects [28,29]. The shortcomings shown above indicate that the currently derived conversion formulas still have large errors for engineering applications, and some of them even fail to meet the engineering accuracy requirements for direct application in actual engineering. In view of this, it is necessary to further study and establish the correlation model between rock uniaxial compression strength and point load strength to improve the prediction accuracy and calculation.
With the continuous progress in computational science, techniques such as machine learning and artificial intelligence, many scholars have applied these to the field of civil and mining engineering to solve many engineering challenges encountered. For example, Yilmaz and Yuksek [16] developed a high-performance adaptive neural network model for predicting gypsum strength and elastic modulus to reduce rock engineering project uncertainty. Ruhul Amin Mozumder [30] used an artificial neural network model to analyze the effect of different geopolymer ratios on uniaxial compression strength and predicted uniaxial compression strength by this model. The usage of these methods improved the prediction model reliability and accuracy to a greater extent and were more accurate in predicting discrete points that deviate from the data center region compared to traditional fitting methods.
Based on these research experiences, a prediction model based on machine learning algorithm can be developed in conjunction with the correlation relationship existing between point load strength and uniaxial compression strength, which can further enhance the prediction accuracy. In this study, 30 groups of corresponding point load tests and uniaxial compression strength tests were conducted, and four basic function types were used to fit and regress the data in order to investigate the correlation relationship between uniaxial compression strength and point load strength. Then a computational model with high accuracy was established for predicting uniaxial compression strength based on support vector machine algorithm. The model can be used to quickly determine the rock strength level at engineering sites, which provides a reference for construction design and safe production.

Rock Sampling and Fabrication
Point load strength tests can be performed on all types of rocks in nature, and test specimens are mainly divided into regular and irregular blocks according to the specimen's integrity. In general, the use of irregular specimens makes the field operation procedure more convenient, but the accuracy of this method is relatively low, so regular blocks for point load strength parameter determination is more in line with the practical needs [31][32][33]. In order to improve accuracy of the point load strength index, this study uses regular cylindrical specimens for testing and analysis.
As we all know, rock is a kind of non-uniform material, and there are pores and micro cracks in it. These factors will affect the macro mechanical strength of rock such as UCS and PLS [34], so complete specifications and no obvious fissures rocks were selected when collecting. The rock samples in this study were obtained from the limestone mine of Sichuan Red Lion Cement in China. The main component of the rock is calcium carbonate, which locally contains 10%~25% argillaceous limestone sand debris. The main density distribution of limestone in the mine is 2.65~2.83 g/cm 3 . The nature of the rock is limestone with uniform block size, good integrity and regular shape, which can meet the test requirements. The rock specimens were prepared on this basis.
In order to prepare test specimens, holes were drilled in the collected rocks to obtain initial specimens with a 50 mm diameter, which were marked. About 4 to 6 initial specimens were obtained on a single rock. Subsequently, these initial specimens were cut and polished to produce UCS and PLS rock specimens with smooth and flat surfaces and no obvious fissures. Ultimately, the uniaxial compression strength test specimens were prepared as standard cylindrical specimens with 100 mm height and 50 mm diameter. The point load test specimens were cylindrical in overall shape with 50 mm diameter, the diametral loading specimens' height was 70 mm, and the axial loading specimen's height was 30 mm.
The rock specimens and equipment used in the test are shown in Figure 1. In this study, 30 specimens were separately prepared for axially loaded point load test, diametral loaded point load test and uniaxial compression test, totaling 90 rock specimens. The water content of the specimens was natural and the tests were run at room temperature. which locally contains 10~25% argillaceous limestone sand debris. The main density distribution of limestone in the mine is 2.65~2.83 g/cm 3 . The nature of the rock is limestone with uniform block size, good integrity and regular shape, which can meet the test requirements. The rock specimens were prepared on this basis. In order to prepare test specimens, holes were drilled in the collected rocks to obtain initial specimens with a 50 mm diameter, which were marked. About 4 to 6 initial specimens were obtained on a single rock. Subsequently, these initial specimens were cut and polished to produce UCS and PLS rock specimens with smooth and flat surfaces and no obvious fissures. Ultimately, the uniaxial compression strength test specimens were prepared as standard cylindrical specimens with 100 mm height and 50 mm diameter. The point load test specimens were cylindrical in overall shape with 50 mm diameter, the diametral loading specimens' height was 70 mm, and the axial loading specimen's height was 30 mm.
The rock specimens and equipment used in the test are shown in Figure 1. In this study, 30 specimens were separately prepared for axially loaded point load test, diametral loaded point load test and uniaxial compression test, totaling 90 rock specimens. The water content of the specimens was natural and the tests were run at room temperature.

Point Load Strength Test
Point load test applies a concentrated load at the tip of the rock specimen until the specimen is destroyed, and the point load strength of the specimen can be obtained by calculation. The larger the concentrated load on the specimen, the greater the point load strength obtained. The concentrated loading methods mainly include diametral and axial, from which the diametral point load strength (DPLS) and axial point load strengths (APLS) are obtained, both of which have comparable reliability and same applied [35,36].
The test method procedure in common use today is mainly the recommended method for point load strength testing proposed by the International Society of Rock Mechanics in 1985. The ISRM method uses the concept of equivalent diameter De to calculate the point load strength. Moreover, to make the resultant data more comparable, the calculated raw Is is generally converted to the point load strength Is (50) at De = 50 mm as follows:

Point Load Strength Test
Point load test applies a concentrated load at the tip of the rock specimen until the specimen is destroyed, and the point load strength of the specimen can be obtained by calculation. The larger the concentrated load on the specimen, the greater the point load strength obtained. The concentrated loading methods mainly include diametral and axial, from which the diametral point load strength (DPLS) and axial point load strengths (APLS) are obtained, both of which have comparable reliability and same applied [35,36].
The test method procedure in common use today is mainly the recommended method for point load strength testing proposed by the International Society of Rock Mechanics in 1985. The ISRM method uses the concept of equivalent diameter De to calculate the point load strength. Moreover, to make the resultant data more comparable, the calculated raw I s is generally converted to the point load strength I s (50) at De = 50 mm as follows: In Formula (1), P represents the point load, De is the equivalent core diameter (mm).
In Formula (2), F represents the correction factor.
In Formula (3), m is the correction index, generally desirable 0.40~0.45, or according to the value of similar rocks to determine the experience.
In this study, point load strength tests were respectively conducted using both diametral loading and axial loading. When the diametral test was used, the equivalent core diameter was calculated according to the following equation: In Formula (4), D is the loading point spacing (mm). When the axial test is used, the equivalent core diameter is calculated according to Formula (5): In Formula (5), W is the width or the smallest section average width through the two loading points (mm).
The test procedure was carried out according to ISRM. During the test, specimens were placed between the upper and lower loading cone heads of the point load strength tester, and the lower cone head was raised by cranking the hand crank oil pump to maintain the loading speed between 0.05~0.1 MPa. The damage occurred within 10~60 s. The rock specimen is damaged between the two loading heads, forming a damage surface through the upper and lower loading points. A typical damage diagram is shown in Figure 1.

Uniaxial Compression Strength Test
Uniaxial compression strength is the main index to evaluate rock stability and rock classification. Put the rock sample between the upper and lower pressing plates of the press. When the sample is crushed, its pressure strength value can be measured, which is the uniaxial compression strength. The greater pressure the rock can bear, the higher the uniaxial compression strength is [37].
Uniaxial compression strength test also refers to the ISRM standard. A total of 30 Rock standard specimens are used in the test. The specimens are tested on the HCT-206A rock true triaxial compression testing machine ( Figure 1) and loaded to failure at the speed of 0.6 MPa/s. The ultimate stress value on the compression surface at the moment of failure is the uniaxial compression strength of the rock, and the specimens will undergo monotonic failure under axial load. The uniaxial compression strength is calculated according to Formula (6): In Formula (6), R is the uniaxial compression strength of the rock, MPa; P is the peak breaking load, N; A is the initial cross-sectional area of the specimen, mm 2 .

Data Analysis and Inspection
In the point load strength test and uniaxial compression strength test, data distortions caused by improper test operations often occur. In this study, such apparently abnormal data were eliminated, and Grubbs' criterion was used to test whether there were suspicious data, and when the suspicious data satisfied Formula (7), they should be eliminated from the group.
In Formula (7), X p is the data being tested, X is the average value, S is the standard deviation, Gp(n) is the test threshold, where p is the confidence probability, taken as 0.95, and n is the number of data, taken as 30.
After rigorous calculations, the Grubbs test value of diametral point load strength data is 1.97, the Grubbs test value of axial loading point load strength data is 2.77, and the Grubbs test value of uniaxial compression strength test data is 2.15, which has no abnormal value at 5% significance level and satisfies the condition of no abnormal value of test data. The above analysis indicates that test data obtained from the experiment has high reliability.

Experimental Data Results
After the test according to the scheme, all the test data were measured in this study, which can be obtained from supplementary materials (Tables S1-S3). After the test according to the scheme, the test data were measured in this study. Then, we drew a Q-Q chart to examine the normality of the data, as shown in Figure 2. The data distribution, mean, and standard deviation of the three experiments also can be obtained in Figure 2 data were eliminated, and Grubbs' criterion was used to test whether there were suspicious data, and when the suspicious data satisfied Formula (7), they should be eliminated from the group.
In Formula (7), Xp is the data being tested, X is the average value, S is the standard deviation, Gp(n) is the test threshold, where p is the confidence probability, taken as 0.95, and n is the number of data, taken as 30.
After rigorous calculations, the Grubbs test value of diametral point load strength data is 1.97, the Grubbs test value of axial loading point load strength data is 2.77, and the Grubbs test value of uniaxial compression strength test data is 2.15, which has no abnormal value at 5% significance level and satisfies the condition of no abnormal value of test data. The above analysis indicates that test data obtained from the experiment has high reliability.

Experimental Data Results
After the test according to the scheme, the test data were measured in this study. Then, we drew a Q-Q chart to examine the normality of the data, as shown in Figure 2. The data distribution, mean, and standard deviation of the three experiments also can be obtained in Figure  It can be seen from the figure that the three experiments data are in accordance with the normal distribution law within the confidence interval of 95% level. The sample average can well represent the overall data level and data concentration degree [38].

Discussion and Analysis
Firstly, the relationship between APLS and UCS is found. As shown in Figure 3a, the correlation coefficient by linear fitting is 0.8883, logarithmic fitting is 0.7768, exponential fitting is 0.9485, and power function fitting is 0.9485. Then four relations between the APLS and UCS are obtained, as shown in Figure 3b. The correlation coefficient of linear fitting is 0.6873, logarithmic fitting is 0.5313, exponential fitting is 0.8846, and power function fitting is 0.7680. It can be seen from the figure that the three experiments data are in accordance with the normal distribution law within the confidence interval of 95% level. The sample average can well represent the overall data level and data concentration degree [38].

Discussion and Analysis
Firstly, the relationship between APLS and UCS is found. As shown in Figure 3a, the correlation coefficient by linear fitting is 0.8883, logarithmic fitting is 0.7768, exponential fitting is 0.9485, and power function fitting is 0.9485. Then four relations between the APLS and UCS are obtained, as shown in Figure 3b. The correlation coefficient of linear fitting is 0.6873, logarithmic fitting is 0.5313, exponential fitting is 0.8846, and power function fitting is 0.7680. From these fitting results, the correlation coefficients using exponential function fitting are greater than other function fitting, whether axial loading or diametral loading, which has better fitting effect, indicating that the relationship between PLS and UCS in this mine is more aptly described in the form of exponential function.
According to the exponential function relationship, the predicted value by using this formula is compared with the experimental value, and the error value is presented in Figure 4. It can be seen from the prediction effect of the final fitting function that even if the exponential function with the largest correlation coefficient is used for strength prediction, the calculation accuracy still deviates too much from the actual requirements, especially in the high-strength rocks. The maximum deviation reached 47.53 MPa, and concentrated in the axial point load strength area. This shows that the traditional fitting method has poor accuracy and cannot adapt to the high dispersion disadvantages of rock strength [39].
Through the above tests and data analysis, it can be known that there is an obvious non-linear relationship between the point load strength of the mine's limestone and the uniaxial compression strength. On this basis, the artificial intelligence algorithm can be used to construct the relationship model between both. Data training can improve the accuracy of the model to meet engineering needs.    From these fitting results, the correlation coefficients using exponential function fitting are greater than other function fitting, whether axial loading or diametral loading, which has better fitting effect, indicating that the relationship between PLS and UCS in this mine is more aptly described in the form of exponential function.
According to the exponential function relationship, the predicted value by using this formula is compared with the experimental value, and the error value is presented in Figure 4. It can be seen from the prediction effect of the final fitting function that even if the exponential function with the largest correlation coefficient is used for strength prediction, the calculation accuracy still deviates too much from the actual requirements, especially in the high-strength rocks. The maximum deviation reached 47.53 MPa, and concentrated in the axial point load strength area. This shows that the traditional fitting method has poor accuracy and cannot adapt to the high dispersion disadvantages of rock strength [39].  From these fitting results, the correlation coefficients using exponential function fitting are greater than other function fitting, whether axial loading or diametral loading, which has better fitting effect, indicating that the relationship between PLS and UCS in this mine is more aptly described in the form of exponential function.
According to the exponential function relationship, the predicted value by using this formula is compared with the experimental value, and the error value is presented in Figure 4. It can be seen from the prediction effect of the final fitting function that even if the exponential function with the largest correlation coefficient is used for strength prediction, the calculation accuracy still deviates too much from the actual requirements, especially in the high-strength rocks. The maximum deviation reached 47.53 MPa, and concentrated in the axial point load strength area. This shows that the traditional fitting method has poor accuracy and cannot adapt to the high dispersion disadvantages of rock strength [39].
Through the above tests and data analysis, it can be known that there is an obvious non-linear relationship between the point load strength of the mine's limestone and the uniaxial compression strength. On this basis, the artificial intelligence algorithm can be used to construct the relationship model between both. Data training can improve the accuracy of the model to meet engineering needs.     Through the above tests and data analysis, it can be known that there is an obvious non-linear relationship between the point load strength of the mine's limestone and the uniaxial compression strength. On this basis, the artificial intelligence algorithm can be used to construct the relationship model between both. Data training can improve the accuracy of the model to meet engineering needs.

Theoretical Basis
Support vector machine (SVM) is a classical artificial intelligence algorithm. SVM model can theoretically approximate any nonlinear function in a global sense. It improves the generalization ability of the model according to the principle of structural risk minimization and can obtain better statistical laws with a small statistical sample size [40]. SVM has the significant advantage of transforming the problem to a high-dimensional feature through nonlinear transformation space, thus avoiding catastrophic data phenomena, local miniaturization, and other problems that often occur in the prediction process [41][42][43].
SVM is a machine learning mechanism. The essence of machine learning is to obtain the relationship estimation between input and output of the system from a given finite sample [44]. Thus, it can make a more accurate prediction of the unknown output. In terms of development, support vector machine theory is developed based on statistical learning theory, which overcomes many shortcomings of classical statistical methods and empirical nonlinear methods [45]. The rock test has the nature of small samples. Each strength index has an obvious non-linear relationship, so it is very suitable to use the SVM model to solve the data prediction problem related to rock specimens. The principle of SVM nonlinear regression is briefly described below [46,47].
For a set of real numbers H = {(x i , y i )}, where i = 1, 2,. . . , n, it is assumed that the functional relation y = g(x) is fitted between x i and y i , where x i is the input quantity in this relation, y i is the output quantity, and n is the sample capacity. For this functional relationship, the fitted relationship equation used in SVM is: In Formula (8), f (x) is the fit function, ϕ(x) is the nonlinear mapping function, ω is the weight vector, and b is the partiality constant. Support vector machines use a structural risk minimization strategy by which ω and b can be obtained: In Formula (9), L ε (y i , f (x i )) is the empirical risk function, 1/2 ω 2 is the regularization term, C is the penalty coefficient, and L ε is the insensitive loss function of ε.
S m can prevent the decision function from over-fitting, and then continue to introduce slack variables ξ i and ξ * i , so that the minimization of S m is abbreviated as: Subject to To determine the range of ω, the solution is performed by Lagrange optimization method, which leads to: In Formula (12), a i and a * i are Lagrange multipliers corresponding to sample of i. The final decision function is obtained by substituting Formula (12) into Formula (8): Minerals 2021, 11, 1387 9 of 14 After establishing the model, R 2 and mean square error (MSE) are used to evaluate the performance of established model. R 2 can be used to measure how well the predicted values fit the true values. In the ideal case, the value of R 2 is 1. The more R 2 converges to 1, the better the prediction performance of the model, which is calculated as in Formula (15): In Formula (15), U p is the predicted value of uniaxial compression strength, U i is the true value of uniaxial compression strength, U represents the average value of the true uniaxial compression strength of the group of tests, and n is the sample capacity.
The MSE refers to the average value of the difference between the predicted value and the actual value after the sum of the squares. It is often used to measure the prediction model accuracy and evaluate the data change degree. The smaller the value of the MSE, the better the model's ability to fit experimental data. The calculation formula is as follows:

Model and Results
Based on the above theoretical basis, the SVM method is used to establish the relationship model between PLS and UCS. As shown in Figure 5, the model uses APLS and DPLS as input vectors. Kernel function and partiality constant b is used as the hidden layer, and y represents the measured UCS value.
After establishing the model, R 2 and mean square error (MSE) are used to evaluate the performance of established model. R 2 can be used to measure how well the predicted values fit the true values. In the ideal case, the value of R 2 is 1. The more R 2 converges to 1, the better the prediction performance of the model, which is calculated as in Formula (15): In Formula (15), Up is the predicted value of uniaxial compression strength, Ui is the true value of uniaxial compression strength, U represents the average value of the true uniaxial compression strength of the group of tests, and n is the sample capacity.
The MSE refers to the average value of the difference between the predicted value and the actual value after the sum of the squares. It is often used to measure the prediction model accuracy and evaluate the data change degree. The smaller the value of the MSE, the better the model's ability to fit experimental data. The calculation formula is as follows:

Model and Results
Based on the above theoretical basis, the SVM method is used to establish the relationship model between PLS and UCS. As shown in Figure 5, the model uses APLS and DPLS as input vectors. Kernel function and partiality constant b is used as the hidden layer, and y represents the measured UCS value. The introduction of the kernel function makes the SVM perfectly solve the problem of sample nonlinearity, so that the inner product operation of the vector in the high-dimensional space is transformed into the function of the vector in the original space, which effectively handles the input of the high-dimensional vector and can directly affect the The introduction of the kernel function makes the SVM perfectly solve the problem of sample nonlinearity, so that the inner product operation of the vector in the highdimensional space is transformed into the function of the vector in the original space, which effectively handles the input of the high-dimensional vector and can directly affect the SVM model performance. At present, the kernel functions commonly used in SVM models include linear kernel function, polynomial kernel function, Gaussian kernel function, Laplacian kernel function, and Sigmoid kernel function. Among them, the Gaussian kernel function has strong learning ability and can be used in the absence of research samples. A good predictive effect is obtained when the prior knowledge is obtained [48]. The SVM constructed in this study uses a Gaussian kernel function, and the specific equation is shown in Formula (17), where γ is the kernel parameter of the kernel function [49]: The SVM is used to establish the uniaxial compression strength prediction model M1, as shown in Figure 6. The input vector currently is the APLS and the corresponding DPLS, and the output vector is the UCS. After sorting out the experimental data, randomly divide 25 groups of data into the training set, and the remaining 5 groups of data into the prediction set. Then run the model. It can be seen from the two model evaluation indexes that the accuracy of model M2 is higher than that of model M3, and the evaluation indexes of the M2 are better. Therefore, for this mining project, it is more suitable to use diametral loading for point load strength test in the field, which can quickly predict the uniaxial compression strength and enhance the efficiency and geotechnical design capability.
Combining the performance results of the three models, the overall performance of model M1 is better and the accuracy is the highest. This shows that the prediction method of using APLS and DPLS to predict UCS is better than using only a single type of point load strength. When the engineering economic conditions permit, the axial point load test and the diametral point load test should be carried out at the same time to obtain the best prediction effect. After a large amount of training with data in the training set, the model accuracy and generalization ability continue to improve. According to the model M1 prediction results, the MSE in the training set is 62.66, R 2 reaches 0.96937, the MSE of the prediction set is 50.05, and R 2 is 0.98, which shows the SVM model has a high fit degree and a small deviation from the average. It is found that the sample area with UCS lower than 100 MPa has higher prediction accuracy than the sample area with UCS higher than 100 MPa. Compared with some existing studies [9,15,26], the R 2 coefficient of model M1 is higher and the prediction performance is better.
According to the practical application experience of SVM model, the main parameters affecting the prediction results are penalty factor C and sensitive parameter γ. By finding the optimal C and γ, the applicability of the model can be optimized. The parameters used in this study are shown in Table 2. Subsequently, we establish the relationship model M2 for predicting UCS by using DPLS. As shown in Figure 7, the results show that the MSE in the training set is 48.24, R 2 is 0.97, the MSE in the prediction set is 76.44, and R 2 is 0.99. There is not much difference between the predicted value and the true value in the prediction set. enhance the efficiency and geotechnical design capability.
Combining the performance results of the three models, the overall performance of model M1 is better and the accuracy is the highest. This shows that the prediction method of using APLS and DPLS to predict UCS is better than using only a single type of point load strength. When the engineering economic conditions permit, the axial point load test and the diametral point load test should be carried out at the same time to obtain the best prediction effect.

Limitations
The main content of this research is the correlation between the PLS and the UCS, and exploring the feasibility and predictive performance of the SVM model in predicting the uniaxial compressive strength by using the point load strength. Although this model has significant performance and advantages in predicting the uniaxial compressive strength of rocks, it also has some limitations that require deeper research.
1. The mechanical strength of rock has a high degree of dispersion, and the amount of data used in this model is small, which may easily lead to errors in the prediction of the strength of individual specimens [50].  Under the same grouping conditions, we establish the relationship model M3 to predict UCS by APLS. As shown in Figure 7, the MSE also obtained in training set is 19.81, R 2 is 0.98, and the MSE in the prediction set is 176.48, R 2 reached 0.96. It can be found that the predicted value of sample No. 4 and sample No. 5 has a large error with the true value. The main reason is that the training data amount is small, and the strength discreteness interferes with the prediction accuracy.
It can be seen from the two model evaluation indexes that the accuracy of model M2 is higher than that of model M3, and the evaluation indexes of the M2 are better. Therefore, for this mining project, it is more suitable to use diametral loading for point load strength test in the field, which can quickly predict the uniaxial compression strength and enhance the efficiency and geotechnical design capability.
Combining the performance results of the three models, the overall performance of model M1 is better and the accuracy is the highest. This shows that the prediction method of using APLS and DPLS to predict UCS is better than using only a single type of point load strength. When the engineering economic conditions permit, the axial point load test and the diametral point load test should be carried out at the same time to obtain the best prediction effect.

Limitations
The main content of this research is the correlation between the PLS and the UCS, and exploring the feasibility and predictive performance of the SVM model in predicting the uniaxial compressive strength by using the point load strength. Although this model has significant performance and advantages in predicting the uniaxial compressive strength of rocks, it also has some limitations that require deeper research.

1.
The mechanical strength of rock has a high degree of dispersion, and the amount of data used in this model is small, which may easily lead to errors in the prediction of the strength of individual specimens [50]. In order to indicate the accuracy and stability of the model, further research is required of more relevant test data sets for verification.

2.
This study only focuses on the model relationship between the point load strength and the uniaxial compressive strength. Other indicators involved in the test, such as specimen loading point spacing D and the minimum cross-sectional width W, have not been involved. The specific relationship between other variables and the uniaxial compressive strength of the rock is not clear.

3.
In order to better improve the prediction performance and accuracy of the model, it is necessary to better optimize and search the built-in parameters of the established SVM model, such as the penalty factor C and kernel parameter γ. The calculation work in this area can be enhanced in the following research.

Conclusions
Through the test results analysis, it was found that the PLS of the limestone in the mine area had an exponential relation with the UCS. Then based on the exponential relationship, the SVM models for predicting UCS by DPLS and APLS are established respectively, and the following conclusions are drawn:

1.
Four fitting methods are both used to obtain the formula with the highest correlation coefficient. The most accurate fitting formula and function for APLS and UCS is R = 6.46e 0.56I s (50) , where R 2 is 0.8846, and the most accurate fitting formula for DPLS and UCS is R = 9.80e 0.47I s (50) , where R 2 is 0.9485. The above formulas are all expressed as exponential functions.

2.
The proposed SVR model can accurately predict the UCS through the PLS. The predicted value is basically the same as the true value, and the maximum error rate is below 15%. The SVM model prediction result is more accurate compared with the traditional fitting methods and it also proves that the SVM algorithm can be better applied to the rock strength prediction scenarios under the condition of small samples.

3.
For the first time, the APLS and DPLS are both used as an entirety to predict the UCS, and the prediction model M1 was established. By comparing the three models, it is found that in the R 2 performance: the correlation coefficient of model M2 is the largest and the correlation coefficient of M3 is the smallest. In the MSE performance, the MSE of model M1 is the smallest and the M3 model is the largest. Combining the R 2 and MSE judgment indexes, model M1 has the best overall performance, which indicates that the accuracy of predicting UCS by using APLS together with DPLS is the highest.

4.
Through comparison, it is found that model M2 has better prediction performance than model M3, which indicates from the side that the prediction of UCS by using DPLS is more accurate than using APLS. If DPLS and APLS cannot be obtained simultaneously at the engineering site, the point load test method of diametral loading can be preferred. Author Contributions: S.L. and Y.W. provided experimental ideas and carried out the experiment; X.X. directed and assisted the experiment; S.L. wrote the paper. All authors have read and agreed to the published version of the manuscript.