Mechanism and Control of Geological Disasters in Deep Engineering Under High Temperature, Ground Stress and Water Pressure
View this Special IssueResearch Article  Open Access
Feisheng Feng, Pan Wang, Zhen Wei, Guanghui Jiang, Dongjing Xu, Jiqiang Zhang, Jing Zhang, "A New Method for Predicting the Permeability of Sandstone in Deep Reservoirs", Geofluids, vol. 2020, Article ID 8844464, 16 pages, 2020. https://doi.org/10.1155/2020/8844464
A New Method for Predicting the Permeability of Sandstone in Deep Reservoirs
Abstract
Capillary pressure curve data measured through the mercury injection method can accurately reflect the pore throat characteristics of reservoir rock; in this study, a new methodology is proposed to solve the aforementioned problem by virtue of the support vector regression tool and two improved models according to Swanson and capillary parachor parameters. Based on previous research data on the mercury injection capillary pressure (MICP) for two groups of core plugs excised, several permeability prediction models, including Swanson, improved Swanson, capillary parachor, improved capillary parachor, and support vector regression (SVR) models, are established to estimate the permeability. The results show that the SVR models are applicable in both high and relatively low porositypermeability sandstone reservoirs; it can provide a higher degree of precision, and it is recognized as a helpful tool aimed at estimating the permeability in sandstone formations, particularly in situations where it is crucial to obtain a precise estimation value.
1. Introduction
In the process of exploring and developing oil and gas fields, permeability has been recognized as one of the key parameters to understand the characteristics of reservoir percolation and to evaluate the productivity of oil and gas wells [1–3]. However, the pore structure of a reservoir is complex, particularly in a reservoir with a high degree of heterogeneity. The accurate prediction of permeability has always been a difficult problem in reservoir evaluation and oil and gas exploration. The research results of Leverett [4] show that mercury injection can be used to generate capillary pressure curves of core samples, which is significant in terms of studying the reservoir pore structure. Thus, a large amount of research has been carried out by scholars for prediction of the reservoir permeability based on mercury injection capillary pressure (MICP) data. Purcell [5] initially proposed that MICP data are useful for estimation of the permeability. If it is assumed that numerous parallel bundles of capillary tubes constitute the porous media, Poiseuille’s equation and Darcy’s law can estimate the permeability. In addition, there are a large number of methods and corresponding models for the sake of estimating the permeability as accurately as possible, including correlating the permeability with MICP parameters. For predicting the permeability from MICP data, there exist certain common models, such as the Swanson parameterbased model [6], the capillary parachor parameterbased model [7], R35 (the pore throat radius corresponding to 35.0% of mercury injection saturation) based model, R50 (the pore throat radius corresponding to 50.0% of mercury injection saturation) based model, and R10 (the pore throat radius corresponding to 10.0% of mercury injection saturation) based model [8–10]. By analyzing the mercury pressure data of cores with different levels of permeability, the parameters reflecting the characteristics of reservoir pore throats are proposed, and a statistical model for the permeability is established [11]. Generally, as is true for researchers who study the permeability estimation by virtue of mercury injection capillary pressure information, the models have only one goal: improvement of the prediction accuracy to meet the needs of reservoir evaluation and oil and gas exploration.
Recently, according to an increasing number of studies, intelligent systems are superior to practical and statistical methods with regard to the relative problems of geosciences and petroleum [12–14]. Therefore, it has been observed that most scholars tend to exploit artificial intelligence techniques to solve their problems in different fields of research and application [15–18]. This paper will attempt to apply artificial intelligence technology to permeability estimation research. Indeed, the SVR method exploits structural risk minimization (SRM) combined with empirical risk minimization (ERM). Under these circumstances, much more reliable results are produced by the SVR method compared to those of neural networks that focus on using the ERM principle. Therefore, the SVR method is used to estimate the permeability from MICP data. In addition, the classical methods proposed by previous researchers have also been improved in this study. Finally, the SVR method is superior to other methods based on comparison results.
2. Theory and Methodology
2.1. The Analysis of MercuryInjection Capillary Pressure (MICP) Curves
Mercury injection capillary pressure (MICP) curves are primarily utilized for several purposes, including the classification of rock types, calculation of the oil saturation according to the height above free water level (HAFWL) method, evaluation of the rock quality, and estimation of the relative permeability [10]. Purcell [5] initially suggested that the permeability estimation could be achieved through MICP data. Once it is assumed that numerous parallel bundles of capillary tubes constitute the porous medium, Poiseuille’s equation and Darcy’s law can be used to estimate the permeability, which is expressed as follows: where is a lithologic parameter; represents the interfacial tension of the twophase fluid, mN/m; denotes the surface contact angle, °; refers to the capillary pressure, MPa; and is the nonwetting phase fluid saturation, %.
Throughout field applications, it can be challenging to estimate the permeability with equation (1) from MICP data in that there are numerous input parameters that need to be obtained first. There are numerous methods and corresponding models for estimating the permeability as accurately as possible, including correlating the permeability with MICP parameters. For predicting the permeability from MICP data, there exist certain common models, such as the capillary parachor parameterbased model and the Swanson parameterbased model.
In most cases, semilog coordinates are employed to plot MICP curves; the linear coordinates of the axis represent the mercury injection saturation (), and the logarithmic coordinates of the axis represent the mercury injection pressure (). Nevertheless, according to Rahul et al. [19], it is said that while saturation data and capillary pressures become marked on a loglog scale, a smooth curve that relatively fits all points and resembles a hyperbola can be described mathematically as follows: where refers to the pressure of mercury intrusion, MPa; denotes the displacement pressure, MPa; is the mercury saturation, %; refers to the intrusion mercury saturation at infinite capillary pressure, %; and is the pore geometry factor.
As shown in Figure 1(a), the inflection point A in the capillary pressure curve necessarily refers to the vertex of the hyperbola on a loglog scale, which is closely related to the situation where part of the effective pore space, which controls the fluid flow, in effect becomes greatly dominated by the nonwetting phase.
(a)
(b)
The value of point A is called the Swanson parameter [6]. According to Guo et al. [7], under steady state conditions, equation (2) is decomposed, and it has been observed that once a chart is designed that uses as the abscissa and as the ordinate, the capillary parachor parameter can be regarded as the highest point, , at which the pore structure of the reservoir can also be reflected. Figure 1(b) shows typical curves of the relations of and versus the intrusion mercury saturation ; it is observed that the curves share similarities with the parabolic curves with a maximum.
2.2. Methodology of Support Vector Regression (SVR)
Vapnik [20] proposed that support vector regression could be regarded as a machine learning technique. Attributed to its superior ability with regard to successfully dealing with many nonlinear regression problems, the SVR method is regarded as a kind of arresting tool featuring promising applications. In addition to the supplementary empirical risk minimization principle that traditionally uses neural networks for the development of accurate models, the structural risk minimization principle is exploited by the SVR method [20–23]. The underlying structure of the SVR method is elaborated on in the following section. Actually, as for SVR regression, the final aim should be that the linear relation is determined between dimensional input vectors x ϵ R^{n} and variables y ϵ R, which tput is expressed in the following way: where and are the slope and offset of the regression line, respectively. To determine the values of and , the minimization of the following equation is necessary:
Introduced by Vapnik (1995), the loss function employed in the SVR method can be insensitive and is expressed as follows:
A dual space can be used to reformulate this problem, which is expressed as follows:
to maximize and subject to where refers to the positive Lagrange multipliers. denotes a regulated positive parameter to determine the tradeoff between the approximation error and the weight vector norm . After calculating the Lagrange multipliers and , training data points that satisfy can be employed to define the function of decision. Therefore, the most appropriate linear hyper surface regression can be provided as follows:
The desired weight vector of the regression hyper plane can be provided by:
In nonlinear regression, the input data are mapped onto a higher dimensional feature space through a kernel function, such that a linear regression hyper plane is generated. In the SVR method, some of the common kernel functions include the radial basis function (RBF), polynomial function, and sigmoid function. Under the condition of nonlinear regression, the identical approach is employed in the linear case to formulate the learning problem again; that is, the nonlinear hyperplane regression function is obtained as follows: where refers to the kernel function, which can be expressed as follows: where and refer to the projections of and in the feature space, respectively. Through research, the Gaussian radial basis function was selected as the kernel function for the SVR model construction. In the above equation, refers to the width of the kernel function. For the reader, a brief description of the SVR method has been provided. In addition, several papers and reviews offer additional information about the SVR method through detailed studies, and information on the SVR theory is contained in the references [24–27].
3. Data Processing and Analysis
3.1. Data Preparation
In mercury injection capillary pressure (MICP) data, several parameters exist that can represent the pore structure, such as the quality coefficient of the reservoir, displacement pressure , capillary pressure midvalue , and mean value of the pore throat radius . The seepage capacity of the pore can be characterized. The displacement pressure can be seen as the capillary pressure, which is related to the maximum interrelated pore throat in the pore system. Moreover, the displacement pressure may be closely connected with the rock permeability [28]. To a certain extent, if the rock specimen has a higher permeability, the value of the displacement pressure will decrease and vice versa. In addition, the displacement pressure belongs to a group of the key parameters by which the reservoir property of reservoir rocks is classified [29–31]. When the mercury saturation becomes 50%, the , regarded as the main calculation parameter of the capillary pressure distribution trend, refers to the median capillary pressure. The overall average pore throat size can be calculated according to the mean value of the pore throat radius [32–34].
Liu et al.[35] analyzed the MICP data, and the porosity and permeability of thirty core samples, as well as the Swanson parameter and capillary parachor parameter of every core sample, were calculated based on the method mentioned, as shown in Figure 1. Table 1 summarizes the results. It can be clearly observed that the 30 core samples are high porosity and permeability rock samples.

3.2. Data Analysis
To evaluate the simple relationships between the quality coefficient of the reservoir and MICP parameters, this paper exploited an analysis technique to highlight the sensitivity of the quality coefficient of the reservoir in the MICP data. This process could also be useful to select input variables in that many unrelated MICP parameters could be eliminated. In addition, the number of input parameters of the model is rather large, which may overwhelm the model and potentially generate a certain amount of noise rather than the signal. Furthermore, the simple linear regression test (cross plotting) was used, of which the correlation coefficient () is regarded as a significant parameter of investigation concerning the influence of various MICP parameters on the laboratorycalculated permeability results. The ratio of overall variance is represented by as calculated for the model, which is obtained as follows:
Cross plots of the MICP parameters of the 30 core samples are presented in Figure 2. It is observed from the linear regression that there was a positive relation between the capillary parachor parameter and permeability, which possessed the highest correlation coefficient of 0.9685. In addition, it seemed that the Swanson parameter, , and porosity had a positive relation with the permeability, and the relations had different correlation coefficients of 0.9169, 0.9092, and 0.3314, respectively. However, it seemed that and exhibited a negative relationship with the permeability. In addition, both accordingly had poor correlation coefficients of 0.2587 and 0.1088, respectively (see Figure 2).
(a)
(b)
(c)
(d)
(e)
(f)
4. Input Selection by Sensitivity Analysis
By virtue of a backpropagation neural network, a stable method was proposed by Dutta and Gupta [36] based on the partial derivative of the output in terms of the th input, which was aimed at capturing the related effort of every input in estimating the output. The equation below is exploited for evaluation of the partial derivative of the output in terms of the ith input: where is a partial derivative of the permeability in terms of the th input, refers to the weight between the output neuron and the th hidden neuron, and refers to the th neuron of the hidden layer. The sum of squares of the partial derivatives () is used for calculating the related effort of backpropagation neural network inputs, which can be expressed as follows:
The contribution of each input variable is given by: where is the relative contribution of the th input variable. The variable with the highest RC affects the output the most.
An improved approach was proposed, and the optimal number of inputs was calculated. First, as shown in Table 2, all of the available MICP parameters were used for the construction of a feed forward backpropagation neural network. Meanwhile, the RC value for every input could be computed by virtue of the sensitivity analysis. Despite the correlation coefficient that had been regarded as the qualitative standard with which to illustrate the relation between inputs and outputs, the sensitivity results represented quantitative standards that tended to be more dependable. The inputs were then ranked according to their RC values.

The optimal number of introduced inputs was an important parameter influencing the design of the SVR model. Thus, on the basis of their RC values, MICP parameters were introduced into the SVR model one by one, and the performance of the SVR model was assessed for every group of inputs. As shown in Figure 3, the results indicated that the optimal SVR model was realized when the top 4 inputs with the highest relative contribution values were used including the capillary parachor parameter, Swanson parameter, , and porosity.
As shown in Figure 4, the flowchart summarizes the aforementioned procedure. In this work, we omitted the detailed introduction of the theory of backpropagation neural networks (BPNNs). Readers who are unfamiliar with BPNNs can learn about the networks in detail with the help of the BPNN results achieved by Mohaghegh.
5. Results and Discussion
In terms of the evaluation and comparison of the performance of the suggested SVR model, certain earlier methods, such as the capillary parachor parameterbased model and Swanson parameterbased model, were employed to estimate permeability values by using the same dataset.
According to Swanson [6], the Swanson parameter has a suitable relationship with the absolute permeability of a core sample according to analysis of 58 excised core specimens that had common combinations of high porosity and high permeability. A relationship between the absolute permeability and the Swanson parameter was established by Swanson. Meanwhile, the following model satisfies the permeability equation as suggested by Swanson: where refers to the permeability, mD; denotes the mercury intrusion saturation, %; represents the mercury intrusion pressure, MPa; is the Swanson parameter, MPa^{1}; and and are statistically undetermined constants, which can be standardized through the use of data sets from mercury injection experiments with core specimens.
As mentioned before, the capillary parachor parameter refers to the maximum of the cross plot of the mercury injection saturation versus the capillary pressure squared ratio of . Figure 1(b) shows the rule of deciding the capillary parachor parameter from the MICP data. By virtue of mercury injection experimental data from eleven core specimens once presented by Purcell [5], a model for estimating the permeability was developed by Guo et al. [7] based on the capillary parachor parameter, which is expressed as follows: where refers to the capillary parachor parameter, MPa^{2}.
Furthermore, resorting to the method suggested by Guo et al. [7], Xiao et al. [37] analyzed the MICP data from 157 core samples covering three different oil fields, and absolute permeability ranges from 0.002 to 1150.0 mD. Moreover, a common formula for estimating the permeability from MICP data was proposed by Xiao et al. [37] that was in line with the capillary parachor parameter. As such, the proposed equation can be expressed as follows: where and refer to statistical model parameters, which are standardized through the usage of data from mercury injection experiments with core specimens.
Recently, Liu et al. [35] improved the permeability estimation model based on the capillary parachor parameter by adding a porosity factor, and the reconstructed model is as follows: where refers to the porosity, %; is the permeability, mD; represents the mercury intrusion saturation, %; is the mercury intrusion pressure, MPa; is the capillary parachor parameter, MPa^{2}; and , , and are undetermined parameters.
To construct a model aimed at the estimation of the permeability from the MICP data, an epsilon support vector regression (εSVR) algorithm was exploited, which included the capillary parachor parameter, Swanson parameter, , and porosity. To optimize the performance of kernel functions and enhance the ultimate precision of the estimation, all data, including inputs and outputs, were normalized within the range of [1,1] before the SVR model construction.
Model construction primarily relied on data from the group of 30 core samples (Table 1) mentioned before. Previous studies had demonstrated that the kernel function could be approximated by radial basis functions (RBFs) because of fewer parameters that required tuning and lower computational costs (Keerthi and Lin, 2003). Thus, an RBF served as a kernel function to construct the SVR model. The relevant parameters within the SVR model and kernel function (, γ, and ε) greatly influenced the performance of the SVR model. Thus, it was important to thoroughly examine these parameters. It had been proposed by You et al. [35] that an appropriate way to conduct this survey was to combine the grid search method with pattern search techniques on the grounds that the area within the optimal points would be determined by the grid search and the global optimal point would be found by a pattern search of the area discovered through grid searching. The particular search range for , γ, and ε was [0.05,500000], [0.000005,20], and [0.0005,50], respectively, whereas the extracted optimal scope parameters were 2.8284, 0.2537, and 0.0051173, respectively.
The predictive performance of the three models was evaluated through the correlation coefficient () and RMSE, which is shown in Figure 5. Figures 5(b) and 5(c) present the cross plots of the modelderived permeability and core assessment results. Among the three different methods employed for permeability estimation, the lowest error and highest correlation coefficient were provided by the SVR model. Furthermore, the Swanson parameterbased model had the largest error and lowest correlation coefficient. Finally, Figure 5(d) shows the relative errors, from which could be observed that the SVR model proposed in this paper outperformed all the other models, while the capillary parachor parameterbased model with the influencing porosity parameter was superior to the Swanson parameterbased model with regard to predicting the permeability. In addition, as shown in Figure 5(d), there is one sample which has maximum relative error for three methods which may be caused by the poor prediction correlation of this point.
(a)
(b)
(c)
(d)
As shown in Figure 5, the MICP data in a normal sandstone reservoir with a high porositypermeability could be used to accurately estimate the permeability. To investigate the performance of the suggested models in predicting the permeability in detail, an additional group of 22 core samples excised from well X4 in the Xujiahe gas sandstone formation with low porositypermeability values in the central Sichuan Basin, southwest China, were employed for assessment of the models proposed in this study.
Analysis was conducted with mercury injection capillary pressure (MICP), porosity, and permeability data of 22 core plugs. In addition, the method mentioned in Figure 6(b) was employed to calculate the Swanson and capillary parachor parameters of every core plug. Finally, the obtained data are shown in Table 3.
(a)
(b)

In the latter stage of this study, the wellknown Swanson and capillary parachor models were improved by adding porosity information since preceding studies had indicated that the porosity had a significant influence on predicting the permeability [19, 36, 38]. By using the data sets obtained from the 22 core samples, the Swanson and capillary parachor parameterbased models and their improved models were established. The final results of the four models constructed are visualized in Figure 7. At the same time, Figures 7(a) and 7(b) depict the Swanson and improved Swanson model, respectively. Figures 7(c) and 7(d) show the capillary parachor and improved capillary parachor model, respectively.
(a)
(b)
(c)
(d)
The following equations show the established model functions.
Swanson model:
Improved Swanson model:
Capillary parachor model:
Improved capillary parachor model: where denotes the permeability, mD; φ is the porosity, %; represents the Swanson parameter, MPa^{1}; and is the capillary parachor parameter, MPa^{2}.
After model construction, the obtained predicted permeability results for each model are presented in Table 3. Figure 8 shows the relation between the predicted permeability and measured permeability of the core specimens; Table 4 presents a set of relative errors of all models.
(a)
(b)
(c)
(d)
(e)
(f)

To assess the models’ performance, two significant concepts were used, including the correlation coefficient () and RMSE (root mean square error). Figures 8(a) and 8(b) show the results of the cross plots of the modelderived permeability and core assessment results; it was observed that the correlation coefficient between the predicted and measured permeability values was 0.9418 for the Swanson model and 0.9475 for the improved Swanson model, which demonstrated that the performance of the Swanson model had barely been improved by adding porosity information. Figures 8(c) and 8(d) show the results of the cross plots of the modelderived permeability and core assessment results, where the correlation coefficient between the predicted and measured permeability values was 0.9401 for the capillary parachor model and 0.9593 for the improved capillary parachor model; these results indicated that the performance of the capillary parachor model had improved by adding porosity information.
Once the value was greater than 0.9, the model performance was considered premium. Generally, the Swanson and capillary parachor models and their improved models had correlation coefficients that were greater than 0.9, demonstrating that permeability prediction had been successfully achieved. Moreover, according to the comprehensive analysis of Figures 5(b), 8(c), and 8(d), it was observed that the improved capillary parachor model could provide improved prediction accuracy compared to those of the previous Swanson and capillary parachor models, particularly in a high porositypermeability sand formation. From Figure 8, it was observed that the prediction accuracy varied among all proposed models, yet the SVR model could provide the lowest error (RMSE =4.3151) and highest correlation coefficient () among the five different models aimed at estimating the permeability. In addition, according to the results of the relative error of all models shown in Figure 8(f), the SVR model predicted permeability values had very low relative errors with regard to the measured values, which verified the robustness of the SVR model.
Finally, Figure 9 shows the statistical analysis results of the proposed models’ performance through error distribution information. From Figure 9(e), it could be observed that the lowest mean absolute error (0.3659) and standard deviation (4.4007) belonged to the SVR model, indicating that the SVR model resulted in an improved prediction performance compared to the other models. Figures 8(a) and 8(b) show that the improved Swanson model was not suitable for the permeability estimation of low porositypermeability sand formations due to the reduced precision in predicting the permeability; the latter is confirmed by its mean absolute error value of 0.9373 and standard deviation value of 8.3921, compared to the mean absolute error value of 0.7255 and standard deviation value of 8.5486 of the Swanson model. The mean absolute error and standard deviation values were 0.625 and 8.9007, respectively, for the capillary parachor parameterbased model, while the mean absolute error and standard deviation values of the improved capillary parachor parameterbased model were 0.4473 and 7.3365, respectively, which were lower than those of the capillary parachor parameterbased model. Thus, it was concluded that the improved capillary parachor parameterbased model had a higher accuracy than the capillary parachor parameterbased model; the latter illustrated that the improved capillary parachor parameterbased model was suitable for the permeability estimation of a low porositypermeability sand formation.
(a)
(b)
(c)
(d)
(e)
As illustrated in Figures 5(c) and 8(e), by comparing the measured and SVR model predicted permeability for both high porositypermeability and relatively low porositypermeability core samples, it was demonstrated that the SVR model was capable of estimating permeability values with a considerable degree of accuracy. Practical results have verified that the SVR model achieved better results than all other models and can be considered a method aimed at estimating the permeability in sandstone formations, particularly in situations where it is crucial to estimate precisely.
6. Conclusions
The permeability, as one of the most significant quality parameters of reservoirs, is capable of providing meaningful data for characterizing reservoirs and petro physical studies when used in combination with the porosity. In fact, certain researchers have attempted to determine the formation permeability by virtue of empirical correlations based on related experiments. Numerous studies have estimated the permeability based on mercury injection capillary pressure (MICP) data due to the significance of the calledfor permeability knowledge. However, the estimation requires methods with great precision. Throughout the paper, the support vector regression method and two improved models based on the Swanson model and capillary parachor parameterbased model were utilized in response to this requirement. The MICP data and porosity information were utilized in the SVR model, including the Swanson parameter, capillary parachor parameter, mean pore throat radius (R_{m}), and porosity, for permeability estimation. All of the models proposed in this study were established based on experimental analysis data from two sets of rock samples, one consisting of 30 core samples with high porositypermeability values and the other consisting of 22 core samples with relatively low porositypermeability values. The results implied that the performance of the SVR model was satisfactory and that the SVR model could extract the implied permeability knowledge from the MICP data and porosity information. A comparison between the SVR model and the other four models, including the Swanson model, improved Swanson model, capillary parachor model and improved capillary parachor model, verified the superiority of the SVR model in terms of the permeability estimation for high porositypermeability and relatively low porositypermeability formations. Through the comparison, it was observed that the SVR model had a higher correlation coefficient and lower root mean square error, mean absolute error, and standard deviation. In addition, for relatively low porositypermeability core samples, the improved capillary parachor model outperformed the capillary parachor model, while the improved Swanson model did not perform better than the Swanson model. The improved capillary parachor model was better than the Swanson and capillary parachor model. To summarize the comparative analysis, although the improved capillary parachor model had been greatly improved in accuracy, the SVR model achieved better results than all the other models and was recognized as a significant tool for estimating the permeability of sandstone formations, particularly in situations where it was crucial to estimate with high degree of precision.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors are very grateful to the key consulting project of the Chinese Academy of Engineering (2020XZ13), the National Key Research and Development Program (2019YFC1904304), and the research on the rheological properties of rock failure in deep mines and its impact on roadway stability (SKLMRDPC19ZZ08) for financial support.
References
 F. Q. Gong, J. Y. Yan, X. B. Li, and S. Luo, “A peakstrength strain energy storage index for rock burst proneness of rock materials,” International Journal of Rock Mechanics and Mining Science, vol. 117, pp. 76–89, 2019. View at: Google Scholar
 Z. Yin, Z. Chen, J. Chang, Z. Hu, H. Ma, and R. Feng, “Crack initiation characteristics of gascontaining coal under gas pressures,” Geofluids, 12 pages, 2019. View at: Google Scholar
 Z. J. Wen, S. L. Jing, and Y. J. Jiang, “Study of the fracture law of overlying strata under water based on the flowstressdamage model,” Geofluids, 12 pages, 2019. View at: Google Scholar
 M. Leverett, “Capillary behavior in porous solid,” Transactions of the AIME, vol. 142, pp. 152–169, 2013. View at: Google Scholar
 W. R. Purcell, “Capillary pressurestheir measurement using mercury and the calculation of permeability therefrom,” Journal of Petroleum Technology, vol. 186, pp. 39–48, 1949. View at: Google Scholar
 B. F. Swanson, “A simple correlation between permeabilities and mercury capillary pressures,” Journal of Petroleum Technology, vol. 33, pp. 2498–2504, 2013. View at: Google Scholar
 B. Guo, A. Ghalambor, and S. Duan, “Correlation between sandstone permeability and capillary pressure curves,” Journal of Petroleum Science & Engineering, vol. 43, no. 34, pp. 239–246, 2004. View at: Publisher Site  Google Scholar
 H. L. Zhang, M. Tu, and H. Cheng, “Breaking mechanism and control technology of sandstone straight roof in thin bedrock stope,” International journal of mining science and technology, vol. 30, no. 2, pp. 259–263, 2020. View at: Publisher Site  Google Scholar
 S. I. Lafage, An alternative to theWinland R35 method for determining carbonate reservoir quality, M.Sc. Thesis, Texas A&M University, USA, 2008.
 V. Nourani, M. T. Alami, and F. D. Vousoughi, “Selforganizing map clustering technique for ANNbased spatiotemporal modeling of groundwater quality parameters,” Journal of Hydroinformatics, vol. 18, pp. 288–309, 2016. View at: Google Scholar
 Z. G. Cheng, S. C. Luo, Z. W. Du, S. Chang, G. L. Li, and H. Li, “The method to calculate tight sandstone reservoir permeability using pore throat characteristic parameters,” Well Logging Technology, vol. 38, pp. 185–189, 2014. View at: Google Scholar
 M. Rajabi, B. Bohloli, and E. Ahangar, “Intelligent approach for prediction of compressional, shear and stoneley wave velocities from conventional well log data: a case study from the Sarvak carbonate reservoir in the Abadan plain (Southern Iran),” Comput. Geosci, vol. 36, pp. 647–664, 2009. View at: Google Scholar
 M. Zoveidavianpoor, A. Samsuri, and S. R. Shadizadeh, “Prediction of compressional wave velocity by an artificial neural network using some conventional well logs in a carbonate reservoir,” Journal of Geophysics and Engineering, vol. 10, pp. 1–13, 2013. View at: Google Scholar
 H. A. Nooruddin, M. E. Hossain, H. AlYousef, and T. Okasha, “Comparison of permeability models using mercury injection capillary pressure data on carbonate rock samples,” Journal of Petroleum Science & Engineering, vol. 121, pp. 9–22, 2014. View at: Publisher Site  Google Scholar
 E. Ebrahimi, M. Monjezi, M. R. Khalesi, and D. J. Armaghani, “Prediction and optimization of backbreak and rock fragmentation using an artificial neural network and a bee colony algorithm,” Bulletin of Engineering Geology & the Environment, vol. 75, no. 1, pp. 27–36, 2016. View at: Publisher Site  Google Scholar
 A. H. Fath, “Application of radial basis function neural networks in bubble point oil formation volume factor prediction for petroleum systems,” Fluid Phase Equilibria, vol. 437, pp. 14–22, 2017. View at: Publisher Site  Google Scholar
 P. Wang, S. Peng, and T. He, “A novel approach to total organic carbon content prediction in shale gas reservoirs with well logs data, Tonghua Basin, China,” Journal of Natural Gas Science & Engineering, vol. 55, pp. 1–15, 2018. View at: Publisher Site  Google Scholar
 R. Rezaee, A. Saeedi, and B. Clennell, “Tight gas sands permeability estimation from mercury injection capillary pressure and nuclear magnetic resonance data,” Journal of Petroleum Science & Engineering, vol. 8889, pp. 92–99, 2012. View at: Publisher Site  Google Scholar
 R. Dastidar, C. H. Sondergeld, and C. S. Rai, “An improved empirical permeability estimator from mercury injection for tight clastic rocks,” Petrophysics, vol. 48, pp. 186–190, 2007. View at: Google Scholar
 V. Vapnik, The Nature of Statistical Learning Theory, Springer, New York, 1995. View at: Publisher Site
 Z. You, Z. Yin, K. Han, D. Huang, and X. Zhou, “A semisupervised learning approach to predict synthetic genetic interactions by combining functional and topological properties of functional gene network,” BMC Bioinformatics, vol. 11, no. 1, pp. 343–355, 2010. View at: Publisher Site  Google Scholar
 B. T. Jiang and F. Y. Zhao, “Combination of support vector regression and artificial neural networks for prediction of critical heat flux,” International Journal of Heat and Mass Transfer, vol. 62, pp. 481–494, 2013. View at: Publisher Site  Google Scholar
 J. Xiang, M. Liang, and Y. He, “Experimental investigation of frequencybased multidamage detection for beams using support vector regression,” Engineering Fracture Mechanics, vol. 131, pp. 257–268, 2014. View at: Publisher Site  Google Scholar
 R. A. Mozumder, B. Roy, and A. I. Laskar, “Support vector regression approach to predict the strength of FRP confined concrete,” Arabian Journal for Science & Engineering, vol. 42, no. 3, pp. 1129–1146, 2017. View at: Publisher Site  Google Scholar
 P. Kang, L. Zhaopeng, Z. Quanle, Z. Zhenyu, and Z. Jiaqi, “Static and dynamic mechanical properties of granite from various burial depths,” Rock Mechanics and Rock Engineering, vol. 52, no. 10, pp. 3545–3566, 2019. View at: Publisher Site  Google Scholar
 K. Peng, Z. Liu, Q. Zou, Q. Wu, and J. Zhou, “Mechanical property of granite from different buried depths under uniaxial compression and dynamic impact: an energybased investigation,” Powder Technology, vol. 362, pp. 729–744, 2020. View at: Publisher Site  Google Scholar
 C. Chen, C. Yan, N. Zhao, B. Guo, and G. Liu, “A robust algorithm of support vector regression with a trimmed Huber loss function in the primal,” Soft Computing, vol. 21, no. 18, pp. 5235–5243, 2017. View at: Publisher Site  Google Scholar
 L. Bo, L. Wang, and L. Jiao, “Recursive finite Newton algorithm for support vector regression in the primal,” Neural Computation, vol. 19, no. 4, pp. 1082–1096, 2007. View at: Publisher Site  Google Scholar
 Y. Lee, S. Oh, and D. H. Choi, “Design optimization using support vector regression,” Journal of Mechanical Science & Technology, vol. 22, no. 2, pp. 213–220, 2008. View at: Publisher Site  Google Scholar
 D. S. Kim, S. W. Lee, and M. G. Na, “Prediction of axial DNBR distribution in a hot fuel rod using support vector regression models,” IEEE Transactions on Nuclear Science, vol. 58, no. 4, pp. 2084–2090, 2011. View at: Publisher Site  Google Scholar
 S. M. Mousavi, R. TavakkoliMoghaddam, B. Vahdani, H. Hashemi, and M. J. Sanjari, “A new support vector modelbased imperialist competitive algorithm for time estimation in new product development projects,” Robotics Computer Integrated Manufacturing, vol. 29, no. 1, pp. 157–168, 2013. View at: Publisher Site  Google Scholar
 H. L. Ritter, “Poresize distribution in porous materials. Pressure porosimeter and determination of complete macropore size distributions,” Industrial and Engineering Chemistry, vol. 17, pp. 782–791, 1945. View at: Google Scholar
 E. C. Donaldson, N. Ewall, and B. Singh, “Characteristics of capillary pressure curves,” Journal of Petroleum Science & Engineering, vol. 6, no. 3, pp. 249–261, 1991. View at: Publisher Site  Google Scholar
 D. Tiab and E. C. Donaldson, Petrophysics, Elsevier, USA, 1996, 270.
 L. JingQiang, Z. ChaoMo, and Z. Zhansong, “Combine the capillary pressure curve data with the porosity to improve the prediction precision of permeability of sandstone reservoir,” Journal of Petroleum Science and Engineering, vol. 139, pp. 43–48, 2016. View at: Publisher Site  Google Scholar
 S. Dutta and J. P. Gupta, “PVT correlations for Indian crude using artificial neural networks,” Journal of Petroleum Science & Engineering, vol. 72, no. 12, pp. 93–109, 2010. View at: Publisher Site  Google Scholar
 Z. X. Xiao, L. Xiao, and W. Zhang, “A new method for calculating sandstone permeability by using capillary pressure curves,” Geophysical Prospecting for Petroleum, vol. 47, pp. 204–207, 2014. View at: Google Scholar
 A. Timur, “An investigation of permeability, porosity, and residual water saturation relationship for sandstone reservoirs,” Log Analyst, vol. 9, pp. 8–17, 1968. View at: Google Scholar
Copyright
Copyright © 2020 Feisheng Feng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.