Next Article in Journal
Laboratory Investigation on Hydrodynamic Performance of an Innovative Aeration Device with a Wave-Driven Heaving Buoy
Previous Article in Journal
A Study on the Multi-Agent Based Comprehensive Benefits Simulation Analysis and Synergistic Optimization Strategy of Distributed Energy in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating the Predictive Performance of Gaussian Process Regression in Evaluating Reservoir Porosity and Permeability

by
Solomon Asante-Okyere
1,2,
Chuanbo Shen
1,2,*,
Yao Yevenyo Ziggah
3,
Mercy Moses Rulegeya
2 and
Xiangfeng Zhu
1,2
1
Key Laboratory of Tectonics and Petroleum Resources, Ministry of Education, China University of Geosciences, Wuhan 430074, China
2
Department of Petroleum Geology, Faculty of Earth Resources, China University of Geosciences, Wuhan 430074, China
3
Department of Geomatic Engineering, Faculty of Mineral Resource Technology, University of Mines and Technology, Tarkwa 00233, Ghana
*
Author to whom correspondence should be addressed.
Submission received: 2 November 2018 / Revised: 15 November 2018 / Accepted: 20 November 2018 / Published: 23 November 2018

Abstract

:
In this paper, a new predictive model based on Gaussian process regression (GPR) that does not require iterative tuning of user-defined model parameters has been proposed to determine reservoir porosity and permeability. For this purpose, the capability of GPR was appraised statistically for predicting porosity and permeability of the southern basin of the South Yellow Sea using petrophysical well log data. Generally, the performance of GPR is deeply reliant on the type covariance function utilized. Therefore, to obtain the optimal GPR model, five different kernel functions were tested. The resulting optimal GPR model consisted of the exponential covariance function, which produced the highest correlation coefficient (R) of 0.85 and the least root mean square error (RMSE) of 0.037 and 6.47 for porosity and permeability, respectively. Comparison was further made with benchmark methods involving a back propagation neural network (BPNN), generalized regression neural network (GRNN), and radial basis function neural network (RBFNN). The statistical findings revealed that the proposed GPR is a powerful technique and can be used as a supplement to the widely used artificial neural network methods. In terms of computational speed, the GPR technique was computationally faster than the BPNN, GRNN, and RBFNN methods in estimating reservoir porosity and permeability.

1. Introduction

The focus of oil and gas exploration in complex and challenging fields has intensively increased in recent years [1]. These complex hydrocarbon fields exhibit a high degree of heterogeneity and also non-uniform spatial distribution of reservoir properties [2,3]. The characterization of such reservoir properties like porosity and permeability cannot be accurately accomplished by developing empirical models from well log parameters since they are inadequate to fully account for the heterogeneous nature of the complex well conditions. An ideal and a more precise approach to directly measure these parameters from core and well testing analysis. However, it is known to be very expensive and time-consuming when conducted on all wells.
This is the major reason why the petroleum industry has embraced the application of computational intelligence techniques, particularly in reservoir characterization, as a result of its inexpensive approach in achieving reliable values [4]. Most studies [5,6,7,8,9,10,11,12,13,14] have reported the capability of various forms and combination of computational intelligence methods, including artificial neural network (ANN), in generating porosity and permeability values. However, most of these computational intelligence methods (including ANN) are slow in generating results mainly due to the iterative tuning of the models’ user-defined parameters and the steepest-gradient training algorithms utilized. Also, when performing prediction, supervised techniques like back propagation neural network (BPNN) and polynomial regression methods mainly engage in error reduction with less regard for the generalization or predictive ability of the model. Generally, this causes overfitting where the computational intelligence model performs well during training but cannot replicate the same results when tested on completely unseen data [15].
To overcome these issues, a non-parametric approach to regression known as Gaussian process regression (GPR) avoids overfitting by defining a function distribution and setting a prior distribution of unlimited possibilities over the function directly [15,16,17,18]. GPR is also known to generalize well due to its preference to a smooth function that accurately explains the training data without manual parameter tuning as has been the case of ANN [19]. This has led to its application in the oil and gas exploration studies. Rawlinson and Vasudeven [20] used GPR to successfully model a deep induction resistivity well log. An observation was made that GPR could predict better than linear interpolation fitting and inverse distance weighting when evaluating and quantifying uncertainties at locations with no data. A Gaussian process as a machine learning tool was also adopted by Yu et al. [15] to predict total organic carbon (TOC) data of the Yangchang formation in the south-eastern Ordos basin. Yu et al. [15] realized that the GPR model selected well logs that had improved correlation with the TOC data to generate predictions with the least error as compared to conventional methods like the Schmoker and Parsey approach.
Based on the computational benefits realized in the aforementioned discussion, the main aim of this study is to extend the application of Gaussian process regression (GPR) and assess its viability as a predictive intelligent model for the evaluation of reservoir porosity and permeability from petrophysical well logs. Furthermore, comparative analysis with benchmark ANN methods of back propagation, generalized regression, and radial basis function was conducted. This was to ascertain the predictive performance of GPR in acquiring porosity and permeability data of the southern basin of the South Yellow Sea.
The South Yellow Sea is the southern major sedimentary basin of the Yellow Sea is located between the Chinese continental landmass and the Korean peninsula (Figure 1a). A central uplift divides the South Yellow Sea basin into the northern and southern basin. This research will be focused on the southern basin due to the recent oil discovery. This article is structured as follows: a detail description of the data used and the GPR will be presented in Section 2. Results and discussion of the predictions of GPR and ANN models will be given in Section 3. The final conclusion is made in Section 4.

2. Data and Methods

2.1. Data Description

The South Yellow Sea basin is a south-west-west oriented rift depression basin located between the Subei basin and Korea peninsular as shown in Figure 1b. It was formed as a result of multi-cycle tectonic action, erosion, and uplift events from the Neoproterozoic through to the Cenozoic era [21,22,23]. It is worth noting that significant oil has been discovered in the Funing, Dainan, and Sanduo Formations of the southern basin. Funing, Dainan, and Sanduo Formations are composed mainly of an interbedding of sandstone and mudstone as indicated in Figure 1c. Muddy limestone can be found in the middle and upper part of the Funing Formation. Also, Figure 1c shows a layer of coal present in the lower portion of the Dainan Formation and muddy limestone in its upper portion.
Two wells (i.e., Well A and Well B) that target the oil-bearing Formations in the southern basin of the South Yellow Sea were considered in the present study (Figure 1b). Well logs, core porosity, and permeability data of Well A were used to train the computational intelligence models while the models’ predictive capabilities were examined on the data of Well B. The input well log parameters were gamma ray (GR), sonic travel time (DT), resistivity (RT), and spontaneous potential (SP). The output parameters were core porosity and permeability. A total of 727 sample data of Well A was used for training and 311 sample data of Well B was used as testing data. Well logs adopted in the study are seen in Figure 2. The statistical description of the well log parameters, porosity and permeability considered for this study are listed in Table 1 and Table 2.

2.2. Gaussian Process Regression

A Gaussian process (GP) is an infinite group of random variables of which any of the finite subsets has a constant joint Gaussian distribution [24,25,26]. A GP is represented by a mean function and a covariance function. Since the GP is a linear combination of random variables having a normal distribution, by simplicity, the mean function is usually assumed to be zero. Assuming a training set y of n number of parameters and having an input matrix x R n and output variable y R , which is expressed to be porosity or permeability. The Gaussian process is therefore represented in Equation (1) as:
y * GP ( m ( x ) , k ( x , x ) )  
where GP is the Gaussian process, m(x) is the mean function and k ( x , x ) is the covariance function. The m(x) in the GP represents the expected value of the function y * at the input matrix point x as expressed in Equation (2):
m ( x ) = E [ f ( x ) ]  
The k ( x , x ) is a measure of the confidence level for m(x) as represented in Equation (3). The covariance function takes any two arguments such that it generates a non-negative covariance matrix K.
k ( x , x ) = E [ ( f ( x ) m ( x ) ) ( f ( x ) m ( x ) ) ]  
The covariance function helps to implicitly specify certain aspects of the model such as smoothness, periodicity, stationarity, etc. The basic and widely used Gaussian process regression (GPR) is composed of a simple zero mean and squared exponential covariance function as expressed in Equation (4):
k ( x , x ) = σ f 2 exp [ r 2 ]  
where
r = | x x | 2 l 2  
σ f and l are the hyperparameters and they influence the performance of GP. σ f is the model noise and l is the length scale. The covariance is close for any set of inputs that are in close proximity; however, the covariance is exponentially decreased as the distance between input parameters increases.
There are various covariance functions (kernel functions) that can be employed in a GPR as denoted in Equations (6)–(9) [16,20]:
  • Matern 3/2:
    k ( x , x ) = σ f 2 ( 1 + 3 r ) exp [ 3 r ]  
  • Matern 5/2:
    k ( x , x ) = σ f 2 ( 1 + 5 r + 5 r 3 ) exp [ 5 r ]  
  • Exponential:
    k ( x , x ) = σ f 2 exp [ r ]  
  • Rational quadratic:
    k ( x , x ) = σ f 2 ( 1 + 1 2 α l 2 ( x x ) Τ ( x x ) ) α  
To estimate the expected function value ( f * ), which is the joint Gaussian prior distribution given the test input ( x * ), Equation (10) can be used:
y * | y N ( f ¯ * , cov ( f * ) )  
f ¯ * is the mean value of the prediction and it gives the best estimate for y * . The variance cov ( f * ) is an indication of the uncertainty of the prediction. The mean prediction, f ¯ * in Equation (11) is a linear combination of the target y while the variance cov ( f * ) is not dependent on the target but only inputs:
( y f * ) N ( m ( X ) m ( X * ) , [ K ( X , X ) + σ n 2 I K ( X , X ) K ( X , X * ) K ( X * , X * ) ] )  
where K ( X , X ) is the covariance matrix of the training dataset. K ( X * , X * ) represents the covariance matrix of the testing data, which represents the N × N * covariance matrix obtained from training and testing data, K ( X * , X ) = K ( X , X * ) Τ . The marginal likelihood over f * is expressed in Equation (12) as:
P ( f * | X , y , x * ) N ( f ¯ * , cov ( f * ) )  
where
f ¯ * = E ( f * | X , y , x * ) = m ( X * ) + K ( X * , X ) [ K ( X , X ) + σ n 2 I ] 1 [ y m ( x ) ]  
cov ( y * ) = K ( X * , X * ) K ( X * , X ) [ K ( X , X ) + σ n 2 I ] 1 K ( X , X * )  
Using the logarithmic identifier to simplify the integral expression of Equation (12) generates the log marginal likelihood given in Equation (15):
log P ( y | X , θ ) = 1 2 y Τ ( K ( X , X ) 1 y ) 1 2 log | K ( X , X ) | n 2 log 2 π  
where θ is the set of hyperparameters needed for a given covariance function. The minimum posterior hyperparameter in the covariance function can be obtained by maximizing the marginal likelihood P ( y | X , θ ) or minimizing the negative log marginal likelihood. The output of the GPR model is presented in terms of its mean and variance.

3. Results

3.1. Developing GPR Models

Different base kernel (covariance) functions, namely squared exponential (covSE), rational quadratic (covRQ), materniso 5/2 (covMatern5), materniso 3/2 (covMatern3), and exponential (CovExp), were assessed when developing the GPR models for porosity and permeability prediction. The idea here was to determine the optimal covariance function that could produce reliable predictions of porosity and permeability that were closely related to the actual values of Well B. For a simple GPR, the constant mean function is known to generate better results than both zero and linear mean functions [16]. The study, therefore, adopted the constant mean function for all the GPR models created.
Correlation coefficient (R) and root mean square error (RMSE) were the statistical indices employed to quantify the performance of the developed models. Here, R examined whether there existed a linear relationship between the predicted and the measured porosity and permeability values. RMSE showed the difference between the predicted outcome and the measured porosity and permeability values. The mathematical expression for R and RMSE are given in Equations (14) and (15), respectively. The GPR models were coded and implemented in MATLAB software (R2016a version, MathWorks Inc., Natick, MA, USA).
R = i = 1 n ( t i t ¯ i ) ( p i p ¯ i ) i = 1 n ( t i t ¯ i ) 2 × i = 1 n ( p i p ¯ i ) 2  
Root   mean   square   error   ( RMSE ) = 1 n i = 1 n ( t i p i ) 2  
where t is the measured parameter value, p is the predicted parameter value, t ¯ is the mean measured parameter value, p ¯ is the mean predicted parameter value, and n is the total number of data points [27].

3.2. Performance of Porosity Models

The prediction results of the developed GPR porosity models based on different covariance functions are presented in Table 3. To do such a performance evaluation, the completely unseen porosity data of Well B was introduced into the GPR porosity models. We can observe that covSE performed worse compared to the other covariance functions. This assertion can also be confirmed by Figure 3A as the correlation between measured porosity and prediction from covSE-GPR was 0.8334. Improved porosity results compared to that of covSE was realized in the cases of covRQ, covMatern5, and covMatern3 as they produced R-values of 0.8415, 0.8412, and 0.8379, respectively (see Figure 3B–D). However, the highest R-value (0.85) was from the covExp porosity model. From Figure 3E, covExp was identified to be the best-ranked GPR porosity model with the least prediction RMSE of 0.0374 as compared to 0.0387 provided by covSE, covRQ, covMatern5, and covMatern3 (Table 3).

3.3. Performance of Permeability Models

The various GPR permeability models based on the different covariance functions were assessed based on the RMSE and the R criteria. In Figure 4A, it can be observed that covSE generalized poorly among the candidate models with R-value of 0.8402. Improved R-values of 0.8428, 0.8429, and 0.8439 were however produced by covRQ, covMatern5, and covMatern3, respectively, as compared to covSE. The correlation plot of prediction from covRQ, covMatern5, and covMatern3 as a function of measured permeability data are shown in Figure 4B–D. Figure 4E reveals covExp scoring the highest R-value of 0.85.
We can also see from Table 4 that covExp performed slightly better than the rest of the kernel functions in predicting the permeability data of Well B. The prediction from the covExp kernel function generated an error of 6.4717, which was marginally lower than the error obtained from covMatern3 (6.5101), covMatern5 (6.5295), and covRQ (6.5309) as listed in Table 4. However, covSE attained a significantly high prediction RMSE of 6.5780 (Table 2). Therefore, based on the covariance function assessed for predicting permeability, covExp was identified to be the best-performing kernel function.

3.4. Comparative Analysis with Artificial Neural Network (ANN)

The covExp-GPR model, which had the best performance in predicting porosity and permeability, was further compared with ANN methods of back propagation, generalized regression, and radial basis functions, which have been widely used for such a purpose. To do a proper competitive comparison, the same petrophysical well log data for Well A was used to train the ANN methods and Well B was used to analyze the optimal ANN models.
A single layer back propagation neural network (BPNN) model was trained using the Levenberg–Marquardt learning algorithm. The optimal model structure, which is highly dependent on the number of hidden neurons, was achieved using a sequential trial and error method based on the lowest RMSE and highest R criteria. It is important to note that the training process was done for 1000 epochs with a learning rate of 0.03 and momentum coefficient of 0.7. In the case of a generalized regression neural network (GRNN) and radial basis function neural network (RBFNN), their output is highly based on the value of the width parameter. Therefore, the optimal width parameter value for GRNN and RBFNN were also achieved through a sequential trial and error approach. Also, a gradient descent rule was adopted to train models of GRNN and RBFNN. The ANN models were coded and implemented in MATLAB R2016a software.

3.5. Comparing CovExp-GPR Porosity Results with ANN

After numerous trials, the BPNN model structure that gave the highest correlation (R) and the least porosity prediction RMSE had 13 hidden neurons. Hence the optimal BPNN porosity model was 4 inputs, 7 hidden neurons, and 1 output variable. On the other hand, the optimal GRNN and RBFNN porosity models with the least RMSE and highest R values had a spread parameter of 0.03. This means that the best performing GRNN porosity model was 4 inputs, a spread parameter of 0.03, and 1 output, while the optimal RBFNN model for predicting porosity composed of 4 inputs, a hidden layer with a maximum of 50 hidden neurons, and a width parameter of 0.3.
The best performing BPNN model took 265.79 s on a Windows AMD Ryzen 5 @ 2 GHz machine (pavilion 15, Hewlett-Packard, Palo Alto, CA, USA) and had a correlation of 0.84 and prediction error of 0.038, as listed in Table 5. Table 5 also shows that the CovExp-GPR model produced comparable porosity prediction results to the widely used BPNN, GRNN, and RBFNN. This can be seen from the R and RMSE results produce by each approach. However, in terms of computational speed, the proposed CovExp-GPR model produced the fastest results. Here, the CovExp-GPR was produced in approximately 22.01 s, while BPNN, GRNN, and RBFNN took 265.79, 29.66, and 96.58 s on a machine having Windows AMD Ryzen 5 @ 2 GHz (see Table 5).

3.6. Comparing Permeability Results with ANN

The best performing BPNN permeability model with improved statistical metrics was observed at a structure of 4 inputs, a hidden layer with 8 neurons, and 1 output. While the optimal GRNN permeability model structure was identified at 4 inputs, a width parameter of 0.03, and 1 output. Similarly, the optimal RBFNN architecture for permeability was achieved with 4 inputs, a hidden layer with a maximum of 50 hidden neurons, a width parameter of 0.3, and 1 output.
From Table 6, it is seen that the proposed CovExp-GPR model produced similar results to the other methods. The reason is related to the reported statistical measures (R and RMSE) in Table 6. Analysis of Table 6 indicates that the proposed CovExp-GPR model predictions were closely related to the observed permeability with a prediction accuracy of 85%. The same was observed for the RBFNN while 86% prediction accuracy was noticed for BPNN and GRNN. The interpretation here shows that a marginal prediction limitation of 15% (GPR) and 14% (BPNN, GRNN, and RBFNN) could not be explained using the respective models. On the basis of the computational speed, the GPR produced the results using a computational time of 20.72 s, while BPNN, GRNN, and RBFNN produced results in a time of 190.04, 28.21, and 100.14 s (Table 6) on a Windows machine with AMD Ryzen 5 @ 2 GHz.

4. Conclusions

In the present research, the suitability of GPR technique was evaluated for the prediction of reservoir porosity and permeability. For this aim, the southern basin of the South Yellow Sea situated between China and the Korean peninsula was selected as a case study application. GPR was developed using five kernel (covariance) functions: covSE, covRQ, covMatern5, covMatern3, and covExp. The motivation was to identify the best-performing (optimal) kernel function suitable for accurately predicting porosity and permeability. It was confirmed from the statistical analysis that the exponential covariance function (covExp) achieved the highest correlation coefficient and least prediction root mean square error when estimating both porosity and permeability.
The performance of the optimal GPR was verified against the benchmark ANN methods of back propagation, generalized regression, and radial basis function. It was further revealed that GPR porosity and permeability prediction models generated faster and comparable results to the widely used ANN methods. The results of the comparative study strongly suggest that the proposed GPR is an efficient approach in obtaining reliable porosity and permeability values and can therefore be used for other wells of the southern basin of the South Yellow Sea for which there are no available porosity and permeability data.

Author Contributions

S.A.-O. conceived and designed the research and prepared the first draft. C.S. supervised the research and analyzed data and results. Y.Y.Z. prepared the data and implemented the algorithms. M.M.R. analyzed the results and performed editing and revision of the first draft. X.Z. contributed to the design of the experiment, analysis and interpretation of data.

Funding

This work was supported by the Major National Science and Technology Programs in the “Thirteenth Five-Year” Plan period (No. 2016ZX05024-002-005, 2017ZX05032-002-004), the Outstanding Youth Funding of Natural Science Foundation of Hubei Province (No. 2016CFA055), the Program of Introducing Talents of Discipline to Universities (No. B14031) and the Wuhan Science and Technology Project (No. 2016070204010145).

Conflicts of Interest

The authors declare no competing interest.

References

  1. Ehsan, M.; Gu, H.; Akhtar, M.M.; Abbasi, S.S.; Ullah, Z. Identification of Hydrocarbon Potential of Talhar Shale: Member of Lower Goru Formation by Using Well Logs Derived Parameters, Southern Lower Indus Basin, Pakistan. J. Earth Sci. 2018, 29, 587–593. [Google Scholar] [CrossRef]
  2. Gao, D.; Cheng, R.; Shen, Y.; Wang, L.; Hu, X. Weathered and Volcanic Provenance-Sedimentary System and Its Influence on Reservoir Quality in the East of the Eastern Depression, the North Yellow Sea Basin. J. Earth Sci. 2018, 29, 353–368. [Google Scholar] [CrossRef]
  3. Mohaghegh, S.; Arefi, R.; Ameri, S.; Aminiand, K.; Nutter, R. Petroleum reservoir characterization with the aid of artificial neural networks. J. Pet. Sci. Eng. 1996, 16, 263–274. [Google Scholar] [CrossRef]
  4. Anifowose, F.A.; Labadin, J.; Abdulraheem, A. Hybrid intelligent systems in petroleum reservoir characterization and modeling: The journey so far and the challenges ahead. J. Pet. Explor. Prod. Technol. 2017, 7, 251–263. [Google Scholar] [CrossRef]
  5. Karimpouli, S.; Fathianpour, N.; Roohi, J. A new approach to improve neural networks’ algorithm in permeability prediction of petroleum reservoirs using supervised committee machine neural network (SCMNN). J. Pet. Sci. Eng. 2010, 73, 227–232. [Google Scholar] [CrossRef]
  6. Iturrarán-Viveros, U.; Parra, J.O. Artificial Neural Networks applied to estimate permeability, porosity and intrinsic attenuation using seismic attributes and well-log data. J. Appl. Geophys. 2014, 107, 45–54. [Google Scholar] [CrossRef]
  7. Bagheripour, P. Committee Neural Network Model for Rock Permeability Prediction. J. Appl. Geophys. 2014, 104, 142–148. [Google Scholar] [CrossRef]
  8. Na’imi, S.R.; Shadizadeh, S.R.; Riahi, M.A.; Mirzakhanian, M. Estimation of Reservoir Porosity and Water Saturation Based on Seismic Attributes Using Support Vector Regression Approach. J. Appl. Geophys. 2014, 107, 93–101. [Google Scholar] [CrossRef]
  9. Plastino, A.; Gonçalves, E.C.; da Silva, P.N.; Carneiro, G.; Azeredo, R.B. Combining Classification and Regression for Improving Permeability Estimations from 1H-NMR Relaxation Data. J. Appl. Geophys. 2017, 146, 95–102. [Google Scholar] [CrossRef]
  10. Anifowose, F.; Abdulraheem, A. Fuzzy logic-driven and SVM-driven hybrid computational intelligence models applied to oil and gas reservoir characterization. J. Nat. Gas Sci. Eng. 2011, 3, 505–517. [Google Scholar] [CrossRef]
  11. Helmy, T.; Fatai, A.; Faisal, K. Hybrid computational models for the characterization of oil and gas reservoirs. Expert Syst. Appl. 2010, 37, 5353–5363. [Google Scholar] [CrossRef]
  12. Anifowose, F.; Labadin, J.; Abdulraheem, A. Improving the prediction of petroleum reservoir characterization with a stacked generalization ensemble model of support vector machines. Appl. Soft Comput. 2015, 26, 483–496. [Google Scholar] [CrossRef] [Green Version]
  13. Ali, D.; Ebrahim, S. Physical properties modeling of reservoirs in Mansuri oil field, Zagros region, Iran. Pet. Explor. Dev. 2016, 43, 611–615. [Google Scholar] [CrossRef]
  14. Baziar, S.; Tadayoni, M.; Nabi-Bidhendi, M.; Khalili, M. Prediction of permeability in a tight gas reservoir by using three soft computing approaches: A comparative study. J. Nat. Gas Sci. Eng. 2014, 21, 718–724. [Google Scholar] [CrossRef]
  15. Yu, H.; Wang, Z.; Rezaee, R.; Zhang, Y.; Xiao, L.; Luo, X.; Wang, X.; Zhang, L. The Gaussian Process Regression for TOC Estimation Using Wireline Logs in Shale Gas Reservoirs. In Proceedings of the International Petroleum Technology Conference, Bangkok, Thailand, 14–16 November 2016. [Google Scholar] [CrossRef]
  16. Aye, S.A.; Heyns, P.S. An integrated Gaussian process regression for prediction of remaining useful life of slow speed bearings based on acoustic emission. Mech. Syst. Signal Process. 2017, 84, 485–498. [Google Scholar] [CrossRef]
  17. Bishop, C.M. Pattern Recognition and Machine Learning; Springer: New York, NY, USA, 2006. [Google Scholar]
  18. Richardson, R.R.; Osborne, M.A.; Howey, D.A. Gaussian process regression for forecasting battery state of health. J. Power Sources 2017, 357, 209–219. [Google Scholar] [CrossRef]
  19. Heyns, T.; de Villiers, J.P.; Heyns, P.S. Consistent haul road condition monitoring by means of vehicle response normalisation with Gaussian processes. Eng. Appl. Artif. Intell. 2012, 25, 1752–1760. [Google Scholar] [CrossRef] [Green Version]
  20. Rawlinson, A.A.; Vasudevan, S. Gaussian Process Modeling of Well Logs. In Proceedings of the IEEE European Modelling Symposium, Madrid, Spain, 6–8 October 2015. [Google Scholar] [CrossRef]
  21. Yi, S.; Yi, S.; Batten, D.J.; Yun, H.; Park, S.J. Cretaceous and Cenozoic non-marine deposits of the Northern South Yellow Sea Basin, offshore western Korea: Palynostratigraphy and palaeoenvironments. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2003, 191, 15–44. [Google Scholar] [CrossRef]
  22. Wu, S.; Ni, X.; Cai, F. Petroleum geological framework and hydrocarbon potential in the Yellow Sea. Chin. J. Oceanol. Limnol. 2008, 26, 23–34. [Google Scholar] [CrossRef]
  23. Pang, Y.; Zhang, X.; Xiao, G.; Wen, Z.; Guo, X.; Hou, F.; Zhu, X. Structural and geological characteristics of the South Yellow Sea Basin in Lower Yangtze Block. Geol. Rev. 2016, 62, 604–616. (In Chinese) [Google Scholar]
  24. Fang, D.; Zhang, X.; Yu, Q.; Jin, T.C.; Tian, L. A novel method for carbon dioxide emission forecasting based on improved Gaussian processes regression. J. Clean. Prod. 2018, 173, 143–150. [Google Scholar] [CrossRef]
  25. Kong, D.; Chen, Y.; Li, N. Gaussian process regression for tool wear prediction. Mech. Syst. Signal Process. 2018, 104, 556–574. [Google Scholar] [CrossRef]
  26. Rasmussen, C.E.; Williams, C.K.I. Gaussian Processes for Machine Learning; The MIT Press: Cambridge, MA, USA; London, UK, 2006. [Google Scholar]
  27. Brantson, E.T.; Ju, B.; Ziggah, Y.Y.; Akwensi, P.H.; Sun, Y.; Wu, D.; Addo, B.J. Forecasting of Horizontal Gas Well Production Decline in Unconventional Reservoirs using Productivity, Soft Computing and Swarm Intelligence Models. Nat. Resour. Res. 2018. [Google Scholar] [CrossRef]
Figure 1. (a) Map showing the location of the South Yellow Sea. (b) Location of the wells considered in the southern basin of South Yellow Sea. (c) Lithology of the formations in the southern basin of the South Yellow Sea.
Figure 1. (a) Map showing the location of the South Yellow Sea. (b) Location of the wells considered in the southern basin of South Yellow Sea. (c) Lithology of the formations in the southern basin of the South Yellow Sea.
Energies 11 03261 g001
Figure 2. Petrophysical well logs of the southern basin of the South Yellow Sea. RT: Resistivity; GR: Gamma ray; DT: Sonic travel time; SP: Spontaneous potential.
Figure 2. Petrophysical well logs of the southern basin of the South Yellow Sea. RT: Resistivity; GR: Gamma ray; DT: Sonic travel time; SP: Spontaneous potential.
Energies 11 03261 g002
Figure 3. Plot of (A) covSE, (B) covRQ, (C) covMatern5, (D) covMatern3, and (E) covExp GPR prediction against measured porosity.
Figure 3. Plot of (A) covSE, (B) covRQ, (C) covMatern5, (D) covMatern3, and (E) covExp GPR prediction against measured porosity.
Energies 11 03261 g003
Figure 4. Plot of (A) covSE, (B) covRQ, (C) covMatern5, (D) covMatern3, and (E) covExp GPR prediction against measured permeability.
Figure 4. Plot of (A) covSE, (B) covRQ, (C) covMatern5, (D) covMatern3, and (E) covExp GPR prediction against measured permeability.
Energies 11 03261 g004
Table 1. Statistical description of well log, porosity, and permeability data of Well A.
Table 1. Statistical description of well log, porosity, and permeability data of Well A.
Well AGR (api)DT (μs/m)RT (Ωm)SP (mv)PorosityPermeability (mD)
Min33.84559.9790.73349.4960.024110.001
Max89.872147.33213.721387.1450.2947859.874
Mean65.364103.0162.166365.8270.174989.2888
SD15.43111.6480.92410.8850.0715211.1292
Table 2. Statistical description of well log, porosity and permeability data of Well B.
Table 2. Statistical description of well log, porosity and permeability data of Well B.
Well BGR (api)DT (μs/m)RT (Ωm)SP (mv)PorosityPermeability (mD)
Min36.58156.2190.732348.9720.023780.001
Max89.848142.1365.575387.0430.2911754.216
Mean64.046102.6182.187367.3560.1856910.9582
SD15.65010.11520.649810.7310.0708112.1387
Table 3. Performance indicators of GPR covariance functions in predicting porosity.
Table 3. Performance indicators of GPR covariance functions in predicting porosity.
Statistical IndexcovSEcovRQcovMarten 5covMatern 3covExp
R0.83340.84150.83790.84120.85
RMSE0.03870.03870.03870.03870.0374
Table 4. Performance indicators of GPR covariance functions in predicting permeability.
Table 4. Performance indicators of GPR covariance functions in predicting permeability.
Statistical IndexcovSEcovRQcovMarten 5covMatern 3covExp
R0.84020.84280.84290.84390.85
RMSE6.57806.53096.52956.51016.4717
Table 5. Performance of porosity models when tested on well B.
Table 5. Performance of porosity models when tested on well B.
Performance IndexCovExp-GPRBPNNGRNNRBFNN
R0.850.840.860.86
RMSE0.0360.0380.0370.036
Computational time (s)22.01265.7929.6696.58
Table 6. Performance of permeability models when tested on well B.
Table 6. Performance of permeability models when tested on well B.
Performance IndexCovExp-GPRBPNNGRNNRBFNN
R0.850.860.860.85
RMSE6.476.276.146.47
Computational time (s)20.72190.0428.21100.14

Share and Cite

MDPI and ACS Style

Asante-Okyere, S.; Shen, C.; Yevenyo Ziggah, Y.; Moses Rulegeya, M.; Zhu, X. Investigating the Predictive Performance of Gaussian Process Regression in Evaluating Reservoir Porosity and Permeability. Energies 2018, 11, 3261. https://0-doi-org.brum.beds.ac.uk/10.3390/en11123261

AMA Style

Asante-Okyere S, Shen C, Yevenyo Ziggah Y, Moses Rulegeya M, Zhu X. Investigating the Predictive Performance of Gaussian Process Regression in Evaluating Reservoir Porosity and Permeability. Energies. 2018; 11(12):3261. https://0-doi-org.brum.beds.ac.uk/10.3390/en11123261

Chicago/Turabian Style

Asante-Okyere, Solomon, Chuanbo Shen, Yao Yevenyo Ziggah, Mercy Moses Rulegeya, and Xiangfeng Zhu. 2018. "Investigating the Predictive Performance of Gaussian Process Regression in Evaluating Reservoir Porosity and Permeability" Energies 11, no. 12: 3261. https://0-doi-org.brum.beds.ac.uk/10.3390/en11123261

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