Next Article in Journal
Congestion Control and Traffic Differentiation for Heterogeneous 6TiSCH Networks in IIoT
Next Article in Special Issue
Prognosis of Bearing and Gear Wears Using Convolutional Neural Network with Hybrid Loss Function
Previous Article in Journal
Spectral Analysis of Macro-Fiber Composites Measured Vibration of Double-Panel Structure Coupled with Solenoids
Previous Article in Special Issue
NLP-Based Approach for Predicting HMI State Sequences Towards Monitoring Operator Situational Awareness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Time Prediction of Rate of Penetration in S-Shape Well Profile Using Artificial Intelligence Models

by
Salaheldin Elkatatny
Petroleum Engineering Department, College of Petroleum Engineering & Geosciences, King Fahd University of Petroleum & Minerals, Box: 5049, Dhahran 31261, Saudi Arabia
Sensors 2020, 20(12), 3506; https://doi.org/10.3390/s20123506
Submission received: 22 May 2020 / Revised: 14 June 2020 / Accepted: 19 June 2020 / Published: 21 June 2020
(This article belongs to the Special Issue Artificial Intelligence and Sensors)

Abstract

:
Rate of penetration (ROP) is defined as the amount of removed rock per unit area per unit time. It is affected by several factors which are inseparable. Current established models for determining the ROP include the basic mathematical and physics equations, as well as the use of empirical correlations. Given the complexity of the drilling process, the use of artificial intelligence (AI) has been a game changer because most of the unknown parameters can now be accounted for entirely at the modeling process. The objective of this paper is to evaluate the ability of the optimized adaptive neuro-fuzzy inference system (ANFIS), functional neural networks (FN), random forests (RF), and support vector machine (SVM) models to predict the ROP in real time from the drilling parameters in the S-shape well profile, for the first time, based on the drilling parameters of weight on bit (WOB), drillstring rotation (DSR), torque (T), pumping rate (GPM), and standpipe pressure (SPP). Data from two wells were used for training and testing (Well A and Well B with 4012 and 1717 data points, respectively), and one well for validation (Well C) with 2500 data points. Well A and Well B data were combined in the training-testing phase and were randomly divided into a 70:30 ratio for training/testing. The results showed that the ANFIS, FN, and RF models could effectively predict the ROP from the drilling parameters in the S-shape well profile, while the accuracy of the SVM model was very low. The ANFIS, FN, and RF models predicted the ROP for the training data with average absolute percentage errors (AAPEs) of 9.50%, 13.44%, and 3.25%, respectively. For the testing data, the ANFIS, FN, and RF models predicted the ROP with AAPEs of 9.57%, 11.20%, and 8.37%, respectively. The ANFIS, FN, and RF models overperformed the available empirical correlations for ROP prediction. The ANFIS model estimated the ROP for the validation data with an AAPE of 9.06%, whereas the FN model predicted the ROP with an AAPE of 10.48%, and the RF model predicted the ROP with an AAPE of 10.43%. The SVM model predicted the ROP for the validation data with a very high AAPE of 30.05% and all empirical correlations predicted the ROP with AAPEs greater than 25%.

1. Introduction

The rate of penetration refers to the speed under which the rock is broken under the bit. It is a measure of the progress or speed of drilling per unit area per unit time and reported in feet per hour (ft/hr) in the oil field unit [1]. In essence, a faster rate of penetration is desirable because it is cost-effective, thereby saving a lot of time and resources that could have been expended. However, other problems are associated with a faster rate of penetration. These include poor hole cleaning and excessive vibrations leading to more complications such as losing the bottom hole assembly (BHA) or wellbore instability [2].
The rate of penetration is affected by several factors which are inseparable. These factors have been categorized into five categories and they include formation characteristics, the efficiency of the rig, properties of mud, hydraulic factors, as well as mechanical factors. As asserted by Hossain and Al-Majed [3], these factors can be grouped into controllable and uncontrollable factors. On the one hand, the controllable factors are those that can be altered quickly without necessarily causing significant operational economies such as the design of the bit, the revolution of string per minute, and the weight of bit. On the other hand, the uncontrollable factors are impossible to alter due to geological and economic reasons which include weight and type of mud, underbalanced in formation in pressure, and the size of the bit. As it relates to the properties of fluids and the effects of the rate of penetration, it is hard to alter a single property without causing an impact on others. Therefore, it makes it difficult to assess the real impact of a specific parameter as it relates to the penetration rate [4]. While drilling a well, three components must be considered [5], namely the weight on bit (WOB), drillstring rotation (DSR), and the pumping rate (GPM) [6].
To help determine the effects of drilling parameters, as well as the properties of mud, on the penetration rate, there are established methods that have been proven to be effective in doing so. Among these established methods are included the basic mathematical and physics equations which help derive the relationships between these items [7]. In addition, the use of correlation and connecting them to obtain such relations has proven to be critical. However, there is no reliable or solid method that currently exists due to the level of complexity associated with the process of drilling. It is also not easy to capture each factor for predicting the penetration rate. Ricardo et al. [8], posited that there was no complete model that could be used to capture all parameters during drilling, as well as the properties of mud which affect the rate of penetration. Therefore, it is important to treat them independently and come up with individual correlations for reliable results.
The S-shape well profile undergoes a build, hold, and then drop section [9]. The well begins as vertical, then at the desired depth, directional BHA is picked to build an angle, then, hold it until a predefined departure is achieved, and then the angle is dropped back to zero [10]. In the upcoming sections, rotary BHA is used, however, when the horizontal displacement is very far, then, a performance motor is selected to allow downhole rotation. This reduces the surface DSR due to torque limitation. In addition, a casing protector needs to be used to reduce casing wear [11]. This profile is selected if there is a surface location such as having a mountain above the reservoir or to get away from subsurface trouble zones such as fractures or salt domes [3]. Different formations are drilled across the tangent which can cause wellbore stability issues across some formations. The sidewall forces keep changing and cycling depending on the BHA position. Torque and drag are present as well, causing some forces to be lost such as vertical forces (WOB) and rotational forces (DSR). Hole cleaning is a big challenge since the cuttings can easily avalanche downward [12,13].

Effect of Drilling Parameters on the Rate of Penetration

Maurer [14] adopted the concept of optimum hole cleaning and derived his rate of penetration (ROP) model for tricone bit type. He assumed that cuttings were being removed as the bit tooth impacted the formation rock. Bingham [15] performed multiple laboratory experiments in order to develop his ROP model. He assumed that the weight on bit threshold could be neglected. This resulted in having the ROP as a function of the rotary speed (DSR) and applied weight on bit. In his model, he developed a WOB exponent which could be calculated through extermination. Bourgoyne and Young [16] presented a mathematical regression model that used existing drilling data to calculate multiple exponents that were required for developing the full model. Each exponent captured a certain physical or mechanical meaning for the drilling process such as the effect of overbalance, overburden pressure, and bit tooth wear. Warren [17] developed his model considering an optimum cleaning scenario for tricone bit type, where the removal of cuttings rate under the bit equaled the rate of generating new cuttings. Al-AbdulJabbar [18] developed a new ROP model and took into consideration drilling mechanical and hydraulic parameters, as well as mud properties. Using nine inputs, two exponents were calculated which were bit exponent and formation compressive strength. Each formation type had its compressive strength coefficient.
Different previous studies have suggested the use of artificial intelligence (AI) to improve the predictability of different parameters related to the oil industry [19,20,21,22,23,24,25,26,27]. Bilgesu et al. [28] suggested the use of AI techniques for ROP prediction. They developed two artificial neural networks (ANN) models for predicting ROP, while drilling through various nine formations in different vertical wells. Amar and Ibrahim [29] developed two ANN models to evaluate the ROP based on the formation depth, ECD, WOB, DSR, pore pressure gradient, drill bit’s tooth wear, and Reynolds number function. A comparison of the prediction power of the developed ANN-based models with the available empirical equations showed that both ANN-based models were highly accurate for estimating the ROP as compared with the empirical equations.
Elkatatny [30] used the ANN feedforward network to predict ROP on three wells. Using two wells, the model was trained on 3333 data points with a correlation coefficient of 0.99 and an average absolute percentage error of 5%. Then, using 2700 unseen data from the third well, the model was able to predict the rate of penetration with a correlation coefficient of 0.99 and an average absolute percentage error of 4%. Al-AbdulJabbar et al. [31] used a feedforward ANN to predict ROP on three well. Using 1500 data points from only single well, the model was able to predict the rate of penetration with a correlation coefficient of 0.92. Later on, the model was used to predict the other two wells with unseen data with a correlation coefficient of 0.95 and 0.94, respectively. Only one well was used in building the model which showed the power of AI in modeling and prediction. Elkatatny et al. [32] demonstrated that once the ANN model was optimized and an empirical correlation was developed, the model could be converted from a black box to a white box making it flexible to deploy in real field applications and environments. Ahmed et al. [33] developed a rate of penetration support vector machine AI model. Using 10 inputs representing drilling mechanical parameters and fluid properties, he developed a resilient ROP model with an AAPE of 2.83%. Al-AbdulJabbar et al. [34] used ANN coupled with self-adaptive differential evolution (SaDE) to predict ROP in horizontal carbonate reservoirs. Using six inputs that coupled drilling mechanical parameters and formation petrophysical properties, such as gamma ray, resistivity, and bulk density, he achieved a strong correlation coefficient of 0.96 and an AAPE of 5.12% after building the model. Using another well with unseen data, he obtained R and AAPE values of 0.95 and 5.8% respectively. The ROP model was turned from a black box to a white box through extracting the weights and biases in a matrix form. A summary of these models including inputs and equations are presented in Table A1, Appendix B.
The main objective of this paper was to build new ROP models for the first time for the S-shape well profile based on the optimized fuzzy inference system (ANFIS), functional neural networks (FN), random forests (RF), and support vector machine (SVM). These models were built to enable a real-time ROP estimation based on the obtained data from the rig real-time sensors such as WOB, DSR, SPP, GPM, and T.

2. Artificial Intelligence Models Theory

ANFIS is the first model used in this study which is a fuzzy subtractive clustering-based fuzzy inference system. The fuzzy inference system consists of a multilayer feedforward adaptive network, and in this network, a specific function is applied to the incoming signal through the training nodes. The model training is conducted in the following two stages: First, the forward pass where the functional signals of the input training data going forward and the parameter in the output is identified through the least square formula, and secondly, the backward pass where the input parameters are updated using the gradient method while the error rates propagate in the opposite [35].
FN is the second model used in this work; this model is very similar to the usual ANN model, but it uses a generalized functional model while ANN uses the sigmoidal common model. In addition, the neuron’s functions of the FN are learned based on the existing training data, and this means the weights associated with these neurons are not needed [36]. The FN model is also characterized by the presence of different arguments in neural functions as compared with the ANN which has one argument [37].
The third model considered, in this study, is the RF which was developed to perform classification and regression tasks [38]. This model combines hundreds of decision trees; every decision tree is trained on different observations, every tree consists of several nodes, and every node considers several features. Then, the final prediction of the random forest model is defined as the average of every tree [39].
The last artificial intelligence model considered in this study is the SVM which was developed earlier in the framework of statistical learning theory as a classifying algorithm. This model uses a multidimensional hyperplane that helps to divide or classify the data into two or more divisions based on the kernel, margin, gamma, and the regularization parameter (C) [40].

3. Data Overview and Preparation

In this work, the field data were obtained from three S-shape wells that share the same hole size and intersect the same geological lithology. These wells were drilled using directional BHA and measurement while drilling (MWD). All the well profiles started as vertical, then, the section was kicked off until the tangent section was reached. The tangent was held at ±25–30°, and then the well was dropped back to zero. The surface data were acquired through real-time sensors, which were recorded on a footage base.
It is very important to recognize that the MWD tool provides information about inclination and azimuth based on a mud plus system which requires no circulation during the measurement. The MWD data are transferred to the surface approximately every 100 ft, while the connection time for the drill pipe stands. Therefore, it is not recommended to include the inclination data as an input parameter to be able to predict the ROP on a footage basis. The geological data, which are obtained from the logging while drilling (LWD) tool, was not available in this study, because it is not common to run the LWD tool outside the reservoir section. At the same time, the changes in pipe speed and torque data compensate for the effect of the geological data on the ROP prediction. The inclination data which is available every 100 ft (this is not a real-time record as the drilling mechanical parameters) did not affect the accuracy of the developed models for ROP prediction. This was confirmed, since the accuracies of the developed models without including the inclination were very high, which could be explained, as shown in Figure 1, by the fact that the change in the ROP took place even when the inclination was constant as compared with the torque which changed consistent with the change of the ROP.
The collected data initially included all types of drilling operations performed in the 12 inches section which included tripping, drilling, and deploying casing. The important and necessary drilling data that were used included weight on bit (WOB), drillstring rotation (DSR), torque (T), pumping rate (GPM), and standpipe pressure (SPP) [18].
Data cleaning was conducted to remove all unrealistic values and outliers. The first step was only to extract the portion where new drilling footage was made while discarding the remaining data, which required a human interface using data filtering and elimination. The second step was to clean the data based on the standard deviation, where all the data values without the range of ±3.0 standard deviation were removed from the data. Although the rotary steerable system (RSS) was used, somehow, the data were noisy. The RSS drive resulted in a much smoother hole as compared to using a mud motor. However, different formations were being drilled with ±25–30° inclination, which caused wellbore stability issues across some formations. In addition, sidewall forces kept changing and cycling depending on the BHA position. This introduced vibrations to the BHA, as well as excessive torque and drag. All of these affected the data quality, especially when the well inclination was approaching zero from its maximum value.
Figure 2 compares the correlation coefficients (R) among the ROP and all training parameters after the cleaning process. As indicated in this figure, the ROP is strongly affected by the T and WOB with Rs of 0.89 and 0.85, respectively, and the ROP has moderate functions with Q, DSR, and SPP with Rs of 0.43, 0.53, and 0.70, respectively.
At the end of a data cleaning process, the data consists of 4012, 1717, and 2500 data points from Well A, Well B, and Well C. Table 1 shows the statistical analysis for the training data obtained from Well A and Well B after the cleaning process. As indicated in this table, Q is between 559 and 993 GPM, DSR is ranging from 99 to 159 rpm, SPP is from 912 to 2490 psi, T is between 4.11 and 15.9 klbf-ft, WOB is from 5.14 to 57.2 klbf, and ROP is ranging from 6.34 to 64.5 ft/hr.

4. Artificial Intelligence Models Optimization

Since the S-Shape profile, with the schematic shown in Figure 3, had dual inclination phenomena (build then drop), two wells were used for training (Well A and Well B), and one for validation (Well C) with a ratio of 2.29:1. The two wells with the highest and lowest curvature were used to predict the well in between. The AI models were built using the ANFIS, FN, RF, and SVM models which were trained using the drilling parameters of the weight on bit (WOB), drillstring rotation (DSR), torque (T), pumping rate (GPM), and standpipe pressure (SPP). The data were loaded into the AI modeling software as inputs and output vectors separately. The data were randomly divided into a 70:30 ratio to be used in the training and testing stages. Well A and Well B were combined in the training-testing phase, where multiple AI model design parameters were varied on a trial-and-error basis to achieve the optimum prediction of the rate of penetration with the least amount of error. The design parameters of the ANFIS model optimized during this step were cluster radius and number of iterations, whereas, for the FN model, the design parameters of the training method and training function type were optimized. The maximum depth and maximum features of the RF model and the kernel, gamma, number of iterations, and C parameter of the SVM were optimized, in this study, to improve the ROP prediction. The governing factors for the selection of the optimum parameters for the AI models were the correlation coefficient (R) and the average absolute percentage error (AAPE), which are defined in Appendix A.
Table 2 summarizes the optimized artificial intelligence models. As indicated in this table, the optimized ANFIS model has a cluster with 0.1 radius and 300 iterations. The FN model has a training method and function type of nonlinear function without iteration terms and forward selection method, respectively. The optimized RF model is characterized by a maximum depth of 29 and maximum features of sqrt. The optimum design parameters of the SVM model are the radial basis function kernel, gamma scale, 1000 iterations, and C of 100.
Using the optimized parameters for the ANFIS-ROP model, sensitivity analysis was performed to evaluate the effect of including the inclination as an input parameter. Table 3 shows that for the training and testing data, there was little effect of including the inclination as an input parameter, where there was an increase of the AAPE from 9.57% to 13.81%, whereas for the validation data there was a significant negative effect of including the hole inclination as an input where the error was increased from 9% to 51%. On the basis of this sensitivity, it was decided to exclude the inclination from the input parameters. Please note that interpolation was performed to have a complete profile for the inclination as it was recorded every 100 ft.

5. Results and Discussion

5.1. Training the Artificial Intelligence Models

Drilling parameters were used as dependent parameters for the rate of penetration prediction using the ANFIS, FN, RF, and SVM models. The developed models were trained using the drilling parameters of the Q, DSR, SPP, T, and WOB.
Figure 4 shows the prediction of the rate of penetration by the developed AI models for the training data. The results showed correlation coefficients (Rs) of 0.97, 0.94, 0.996, and 0.83 between the actual and predicted ROP with average absolute percentage errors (AAPEs) of 9.50%, 13.44%, 3.25%, and 31.72% for the ANFIS, FN, RF, and SVM models, respectively. These results show that the ANFIS, FN, and RF models’ performances were very satisfactory, whereas the SVM. model was not able to accurately predict the ROP.

5.2. Testing Artificial Intelligence Models

The optimized ANFIS, FN, RF, and SVM models were, then, tested on 1717 data points. Figure 5 shows the prediction of the rate of penetration by the developed AI models for the testing data. These results confirmed the high accuracy of the optimized AI models in estimating the ROP for the testing data. The ANFIS model predicted the ROP with an R of 0.97 and an AAPE of 9.57%. The FN model predicted the ROP with an R and AAPE of 0.96 and 11.20%, respectively. The RF model was the best to estimate the ROP for the testing dataset with an R and AAPE of 0.98 and 8.37% (Figure 5c), respectively, whereas the SVM model predicted the ROP with the least accuracy among the other models. The SVM model estimated the ROP with an R of 0.83 only and a very high AAPE of 31.73% (Figure 5d).

5.3. Validating the Artificial Intelligence Models and Comparison with the Published ROP Correlations

To validate the developed AI models, and ensure that they did not memorize the data, Well C was used to predict the rate of penetration with more than 2500 unseen data points. Figure 6 compares the predicted versus an actual rate of penetration for the validation data. The ANFIS model overperformed all AI models developed in this study and all other correlations in estimating the ROP for the validation data collected from Well C. The R and AAPE for estimating the ROP using the ANFIS model were 0.97 and 9.06%, respectively. The FN model also overperformed the previously developed empirical correlations for estimating the ROP with an R of 0.96 and an AAPE of 10.48%. The RF model also overperformed all correlations for estimating the ROP for the validation data where it predicted the ROP with an R of 0.96 and an AAPE of 10.43%. Whereas the SVM model predicted the ROP for the validation data with very low accuracy in terms of an R and AAPE of 0.82 and 30.05% (Figure 6), respectively.
All available empirical correlations predicted the ROP with very high AAPEs and low Rs. Bingham, Maurer, and Bourgoyne–Young correlations estimated the ROP with Rs of 0.87, 0.85, and 0.84 and AAPEs of 24.83%, 27.43%, and 29.55%, respectively. These results confirmed the high accuracy of the developed AI models as compared with the available empirical correlations. The results in Figure 7 show the comparison of the AAPEs for the validation data, and confirms the high accuracy of the developed ANFIS, FN, and RF models as compared with the SVM model and other available empirical correlations. Table A2 (Appendix C) shows the comparison of the numerical values of the ROP estimated using the empirical correlations with the ROP values estimated with the AI models for 20 data points selected randomly from Well C. The values reported in this table confirm that the values predicted with the ANFIS, FN, and RF models are closer to the actual ROP as compared with the values predicted using the SVM model or the empirical correlations.

6. Conclusions

Actual field data from three S-shape wells were used for the rate of penetration prediction using AI models of the ANFIS, FN, RF, and SVM based on five inputs of the DSR, SPP, Q, T, and WOB. On the basis of the results, the following can be concluded:
  • The ANFIS, FN, and RF models could effectively predict the ROP from the drilling parameters in the S-shape well profile for training, testing, and validation data, whereas the SVM model showed very low accuracy in estimating the ROP.
  • The developed ANFIS, FN, and RF models predicted the ROP with AAPEs of 9.50%, 13.44%, and 3.25%, respectively, for the training data.
  • For the testing data, the optimized ANFIS, FN, RF models estimated the ROP with AAPEs of 9.57%, 11.20%, and 8.37%, respectively.
  • The developed ANFIS, FN, and RF models outperformed the SVM model and the three published empirical correlations for estimating the ROP for the validation data.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

ANNArtificial neural networks
AAPEAverage absolute percentage error
AIArtificial intelligence
BHABottom hole assembly
FNFunctional neural networks
GPMGallon per minute
GRGamma ray
LWDLogging while drilling
MDMeasured depth
MWMud weight
MWDMeasurement while drilling
MFVMarsh funnel viscosity
PVPlastic viscosity
QFlow rate
RCorrelation coefficient
RFRandom forests
RMSERoot mean square error
ROPRate of penetration
DSRDrillstring rotation
SaDESelf-adaptive differential evolution
SPPStandpipe pressure
SVMSupport vector machine
TFATotal flow area
TVDTrue vertical depth
UCSUnconfined compressive strength
WOBWeight on bit

Appendix A

The correlation coefficient (R) and the average absolute percentage error (AAPE) are defined by the following equations, respectively:
R = 1 i N ( y i a c t y i p r e d ) 2 i N ( y i a c t ) 2 i N ( y i p r e d ) 2 N
A A P E = 1 N i = 1 N | y i a c t y i p r e d | y i a c t
where yiact is the actual value, yipred is the predicted value, and N is the number of data points. Higher accuracy and better performance of prediction are indicated by higher values of R and lower values of the root mean square error (RMSE).

Appendix B

The table below summarizes multiple published ROP models including ones that were developed empirically or using artificial intelligence techniques.
Table A1. Summary of previous ROP models.
Table A1. Summary of previous ROP models.
ReferenceMethodInput ParametersFormulaConstants
Maurer [14]EmpiricalDSR, WOB, bit size, uniaxial compressive strength (UCS) R = k N W 2 D 2 S 2 K = 60,087 × 103
Bingham [15]EmpiricalDSR, WOB, bit size R = K ( W d b ) a N c K = 1.04; a = 0.081; c = 1.08
Bourgoyne and Young’s [16]EmpiricalTrue vertical depth (TVD), mud weight (MW), pore pressure, WOB, DSR, bit wear, Q, TFA d d t ( R ) = e ( a 1 + i = 2 8 a i x i ) a1 = 14.35; a2 = −1.5 × 10−3; a3 = −9 × 10−4; a4 = 3 × 10−4; a5 = 1.496; a6 = 1.269; a7 = 0.118
Warren [17]EmpiricalWOB, DSR, bit size, UCS R = [ a S 2 d b 2 N b W 2 + c N d b ] 1
Osgouei [4]EmpiricalTVD, MW, pore pressure, WOB, DSR, bit wear, Q, TFA, bit type, inclination d d t ( R ) = e ( a 1 + i = 2 11 a i x i )
Al-AbdulJabbar [18]EmpiricalDSR, WOB, torque, SPP, Q, bit size, UCS, MW, drilling fluid plastic viscosity (PV) R O P = W O B a × R P M × T × S P P × G P M d b 2 × ρ × P V × U C S b
Elkatatny et al. [41]ANNQ, SPP, DSR, WOB, MW, torque, PV-
Al-AbdulJabbar et al. [31]ANNWOB, DSR, SPP, Q, torque-
Ahmed et al. [33]SVMQ, SPP, DSR, WOB, MW, torque, drilling fluid yield point, marsh funnel viscosity, PV, solids volume-
Al-AbdulJabbar et al. [34]ANN SaDEDSR, torque, WOB, GR, formation deep resistivity, formation bulk density-

Appendix C

Table A2. Shows the main input parameters for ROP prediction using previous empirical correlations, and the actual and predicted ROP values.
Table A2. Shows the main input parameters for ROP prediction using previous empirical correlations, and the actual and predicted ROP values.
S-Shape Measured Length (ft)InputsActual ROP, ft/hr* Empirical Correlations-Based ROP, ft/hrAI-Based ROP, ft/hr
GPMRPMWOB, KlbfTorque, Klbf-ftSPP, psiDensity, PCFBit Size, inPVUCS, psiBingham [15]Maurer [14]B&Y [16]ANFISFNRFSVM
116843.9113132.4610.41715.7767123035,00036.0432.0623.0416.7736.7536.8436.6323.42
133848.33132319.881704.9967123035,00034.2430.8121.1715.8733.9132.8136.6222.52
165838.7913132.818.031667.6567123035,00026.0232.4223.5417.2126.4323.7628.2622.02
241843.7613626.397.321708.8167123025,00019.2526.9315.8113.3619.9518.5520.2221.74
244841.1813626.268.291706.7667123025,00024.2226.7915.6513.2724.2522.9925.1421.88
249844.8913725.797.821697.2967123025,00017.8826.5015.2113.0517.8820.4020.8921.28
254845.4913626.667.741681.8167123025,00020.6927.2116.1313.6020.8220.0220.6520.93
288844.3113627.627.311675.3267123025,00019.6128.2317.3214.4818.6418.4520.1120.96
306843.3313627.678.261690.0267123025,00022.8528.2817.3814.5723.8923.0323.1121.56
362940.8113727.318.22027.2767123030,00023.0728.1217.0514.8624.1324.5723.6328.72
365948.9813827.668.342050.4967123030,00023.528.7217.6215.3123.1225.5522.0229.39
388917.1413628.796.571961.9467123030,00018.0929.4718.8115.9717.5217.8318.8627.45
1076698.3812831.468.31371.9967123025,00027.6130.2721.1419.8126.6925.0826.2321.28
1326892.4612841.828.52077.7767123025,00035.8340.6637.3633.2837.2535.1736.0136.55
1488898.0312838.497.952065.3767123025,00033.8137.3131.6530.4833.7030.7034.6434.88
1720887.6912830.4711.382040.2767123025,00042.8229.2819.8422.4140.9246.1043.6632.95
1723890.2114034.1610.452060.7667123025,00044.9136.3127.2729.8344.0843.1248.1334.97
1729891.1914035.5310.892048.167123025,00049.9637.8229.5031.6748.3346.2054.5534.72
1827856.2812836.489.272058.2167123025,00037.9235.2928.4329.9536.9237.9837.8736.65
1831851.5714637.619.472075.9267123025,00048.2341.9734.4737.2046.0841.2847.4738.81
* The inputs for Bingham, RPM and WOB; for Maurer, RPM, WOB, UCS, and bit size; for Bourgoyne–Young, RPM, WOB, UCS, bit size, GPM, torque, SPP, density, and PV.

References

  1. Bourgoyne, A.T., Jr.; Millheim, K.K.; Chenevert, M.E.; Young, F.S. Applied Drilling Engineering; Society of Petroleum Engineers: Houston, TX, USA, 1986; Volume 1. [Google Scholar]
  2. Akgun, F. How to Estimate the Maximum Achievable Drilling Rate without Jeopardizing Safety. Soc. Pet. Eng. 2002. [Google Scholar] [CrossRef]
  3. Hossain, M.E.; Al-Majed, A.A. Fundamentals of Sustainable Drilling Engineering; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015. [Google Scholar] [CrossRef]
  4. Osgouei, R.E. Rate of Penetration Estimation Model for Directional and Horizontal Wells; The Graduate School, Middle East Technical University: Ankara, Turkey, 2007. [Google Scholar]
  5. Lyons, W.; Gary JPlisga, B.; Lorenz, M. Standard Handbook of Petroleum and Natural Gas Engineering; Gulf Professional Publishing; Elsevier: Amsterdam, The Netherlands, 2004. [Google Scholar]
  6. Mitchell, R.F.; Miska, S.Z. Fundamentals of Drilling Engineering; Society of Petroleum Engineers: Richardson, TX, USA, 2011. [Google Scholar]
  7. Hegde, C.; Daigle, H.; Millwater, H.; Gray, K. Analysis of rate of penetration (ROP) prediction in drilling using physics-based and data-driven models. J. Pet. Sci. Eng. 2017, 159, 295–306. [Google Scholar] [CrossRef]
  8. Ricardo, J.; Mendes, P.; Fonseca, T.C.; Serapaio, A.B.S. Applying a Neuro-model Reference Adaptive Controller in Drilling Optimization. World Oil Magazine, 1 October 2007; Volume 228, 29–38. [Google Scholar]
  9. Allouche, E.N.; Ariaratnam, S.T.; Lueke, J.S. Horizontal directional drilling: Profile of an emerging industry. J. Constr. Eng. Manag. 2000, 126, 68–76. [Google Scholar] [CrossRef]
  10. Short, J.A. Introduction to Directional and Horizontal Drilling; Pennwell Corporation: Tulsa, OK, USA, 1993. [Google Scholar]
  11. Inglis, T. Directional Drilling; Springer Science & Business Media: Berlin, Germany, 2013. [Google Scholar] [CrossRef]
  12. Rabia, H. Well Engineering and Construction; Entrac Consulting: London, UK, 2001. [Google Scholar]
  13. Mahmoud, A.A.; Elzenary, M.; Elkatatny, S. New Hybrid Hole Cleaning Model for Vertical and Deviated Wells. J. Energy Resour. Technol. 2020, 142, 034501. [Google Scholar] [CrossRef]
  14. Maurer, W. The “Perfect—Cleaning” Theory of Rotary Drilling. J. Pet. Technol. 1962, 14, 1270–1274. [Google Scholar] [CrossRef]
  15. Bingham, M.G. A New Approach to Interpreting Rock Drillability; Petroleum Publishing Co.: Houston, TX, USA, 1965. [Google Scholar]
  16. Bourgoyne, A.T.; Young, F.S. A Multiple Regression Approach to Optimal Drilling and Abnormal Pressure Detection. Soc. Pet. Eng. J. 1974, 14, 371–384. [Google Scholar] [CrossRef]
  17. Warren, T.M. Penetration Rate Performance of Roller Cone Bits. Soc. Pet. Eng. J. 1987. [Google Scholar] [CrossRef]
  18. Al-AbdulJabbar, A.M. Utilizing Field Data to Understand the Effect of Drilling Parameters and Mud Rheology on Rate of Penetration in Carbonate Formations. King Fahd University of Petroleum & Minerals (KFUPM). 2017. Available online: https://eprints.kfupm.edu.sa/id/eprint/140167 (accessed on 25 February 2020).
  19. Mahmoud, A.A.; Elkatatny, S.; Al-Shehri, D. Application of Machine Learning in Evaluation of the Static Young’s Modulus for Sandstone Formations. Sustainability 2020, 12, 1880. [Google Scholar] [CrossRef] [Green Version]
  20. Mahmoud, A.A.; Elkatatny, S.; Chen, W.; Abdulraheem, A. Estimation of Oil Recovery Factor for Water Drive Sandy Reservoirs through Applications of Artificial Intelligence. Energies 2019, 12, 3671. [Google Scholar] [CrossRef] [Green Version]
  21. Mahmoud, A.A.; Elkatatny, S.; Mahmoud, M.; Abouelresh, M.; Abdulraheem, A.; Ali, A. Determination of the total organic carbon (TOC) based on conventional well logs using artificial neural network. Int. J. Coal Geol. 2017, 179, 72–80. [Google Scholar] [CrossRef]
  22. Ahmed, S.A.; Mahmoud, A.A.; Elkatatny, S. Fracture Pressure Prediction Using Radial Basis Function. In Proceedings of the AADE National Technical Conference and Exhibition, AADE-19-NTCE-061, Denver, CO, USA, 9–10 April 2019. [Google Scholar]
  23. Mahmoud, A.A.; Elkatatny, S.; Abdulraheem, A.; Mahmoud, M.; Ibrahim, O.; Ali, A. New Technique to Determine the Total Organic Carbon Based on Well Logs Using Artificial Neural Network (White Box). In Proceedings of the 2017 SPE Kingdom of Saudi Arabia Annual Technical Symposium and Exhibition, Dammam, Saudi Arabia, 24–27 April 2017. Paper SPE-188016-MS. [Google Scholar] [CrossRef]
  24. Ahmed, S.A.; Mahmoud, A.A.; Elkatatny, S.; Mahmoud, M.; Abdulraheem, A. Prediction of Pore and Fracture Pressures Using Support Vector Machine. In Proceedings of the International Petroleum Technology Conference, Beijing, China, 26–28 March 2019. IPTC-19523-MS. [Google Scholar] [CrossRef]
  25. Elkatatny, S.; Al-AbdulJabbar, A.; Mahmoud, A.A. New Robust Model to Estimate the Formation Tops in Real-Time Using Artificial Neural Networks (ANN). Petrophysics 2019, 60, 825–837. [Google Scholar] [CrossRef]
  26. Mahmoud, A.A.; Elkatatny, S.; Ali, A.; Abouelresh, M.; Abdulraheem, A. Evaluation of the Total Organic Carbon (TOC) Using Different Artificial Intelligence Techniques. Sustainability 2019, 11, 5643. [Google Scholar] [CrossRef] [Green Version]
  27. Mahmoud, A.A.; Elkatatny, S.; Ali, A.; Moussa, T. Estimation of Static Young’s Modulus for Sandstone Formation Using Artificial Neural Networks. Energies 2019, 12, 2125. [Google Scholar] [CrossRef] [Green Version]
  28. Bilgesu, H.; Tetrick, L.; Altmis, U.; Mohaghegh, S.; Ameri, S. A new approach for the prediction of rate of penetration (ROP) values. In Proceedings of the SPE Eastern Regional Meeting, Lexington, Kentucky, 22–24 October 1997. SPE-39231-MS. [Google Scholar] [CrossRef]
  29. Amar, K.; Ibrahim, A. Rate of penetration prediction and optimization using advances in artificial neural networks, a comparative study. In Proceedings of the 4th International Joint Conference on Computational Intelligence, Barcelona, Spain, 5–7 October 2012; pp. 647–652. [Google Scholar] [CrossRef] [Green Version]
  30. Elkatatny, S. New Approach to Optimize the Rate of Penetration Using Artificial Neural Network. Arabian J. Sci. Eng. 2017, 43, 6297–6304. [Google Scholar] [CrossRef]
  31. Al-AbdulJabbar, A.; Elkatatny, S.; Mahmoud, M.; Abdulraheem, A. Predicting Rate of Penetration Using Artificial Intelligence Techniques. In Proceedings of the SPE Kingdom of Saudi Arabia Annual Technical Symposium and Exhibition, Dammam, Saudi Arabia, 23–26 April 2018. [Google Scholar] [CrossRef]
  32. Elkatatny, S.; Tariq, Z.; Mahmoud, M. Real time prediction of drilling fluid rheological properties using Artificial Neural Networks visible mathematical model (white box). J. Pet. Sci. Eng. 2016, 146, 1202–1210. [Google Scholar] [CrossRef]
  33. Ahmed, A.S.; Elkatatny, S.; Abdulraheem, A.; Mahmoud, M.; Ali, A.Z.; Mohamed, I.M. Prediction of Rate of Penetration of Deep and Tight Formation Using Support Vector Machine. In Proceedings of the SPE Kingdom of Saudi Arabia Annual Technical Symposium and Exhibition, Dammam, Saudi Arabia, 23–26 April 2018. [Google Scholar] [CrossRef]
  34. Al-AbdulJabbar, A.; Elkatatny, S.; Mahmoud, A.A.; Moussa, T.; Al-Shehri, D.; Abughaban, M.; Al-Yami, A. Prediction of the rate of penetration while drilling horizontal carbonate reservoirs using the self-adaptive artificial neural networks technique. Sustainability 2020, 12, 1376. [Google Scholar] [CrossRef] [Green Version]
  35. Jang, R. ANFIS Adaptive Network Based Fuzzy Inference System. IEEE Trans. Syst. Man Cybern. 1993, 23. [Google Scholar] [CrossRef]
  36. Bello, O.; Asafa, T. A Functional Networks Softsensor for Flowing Bottomhole Pressures and Temperatures in Multiphase Production Wells. In Proceeding of the SPE Intelligent Energy Conference & Exhibition, Utrecht, The Netherlands, 1–3 April 2014. [Google Scholar] [CrossRef]
  37. 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]
  38. Zhang, C.; Ma, Y. Ensemble Machine Learning; Springer: Boston, MA, USA, 2012. [Google Scholar]
  39. Hedge, C.; Wallace, S.; Gray, K. Using Trees, Bagging, and Random Forest to Predict Rate of Penetration during Drilling. In Proceeding of the SPE Middle East Intelligence Oil and Gas Conference and Exhibition, Dubai, UAE, 15–16 September 2015. SPE-176792-MS. [Google Scholar] [CrossRef]
  40. Rui, J.; Zhang, H.; Zhang, D.; Han, F.; Guo, Q. Total organic carbon content prediction based on support-vector-regression machine with particle swarm optimization. J. Pet. Sci. Eng. 2019, 180, 699–706. [Google Scholar] [CrossRef]
  41. Elkatatny, S.M.; Tariq, Z.; Mahmoud, M.A.; Al-AbdulJabbar, A. Optimization of Rate of Penetration using Artificial Intelligent Techniques. In Proceedings of the 51st U.S. Rock Mechanics/Geomechanics Symposium, San Francisco, CA, USA, 25–28 June 2017; ARMA-2017-0429. Available online: https://www.onepetro.org/download/conference-paper/ARMA-2017-0429?id=conference-paper%2FARMA-2017-0429 (accessed on 20 January 2020).
Figure 1. The correlation coefficient between the rate of penetration (ROP) and the training parameters.
Figure 1. The correlation coefficient between the rate of penetration (ROP) and the training parameters.
Sensors 20 03506 g001
Figure 2. The correlation coefficient between the ROP and the training parameters.
Figure 2. The correlation coefficient between the ROP and the training parameters.
Sensors 20 03506 g002
Figure 3. The schematic of the S-shape well trajectory.
Figure 3. The schematic of the S-shape well trajectory.
Sensors 20 03506 g003
Figure 4. Predicted and actual ROP along the S-shape section for the training data (4012 data points) using (a) ANFIS; (b) functional neural networks (FN); (c) random forests (RF); and (d) support vector machine (SVM) models.
Figure 4. Predicted and actual ROP along the S-shape section for the training data (4012 data points) using (a) ANFIS; (b) functional neural networks (FN); (c) random forests (RF); and (d) support vector machine (SVM) models.
Sensors 20 03506 g004
Figure 5. Predicted and actual ROP along the S-shape section for the testing data (1717 data points) using (a) ANFIS; (b) FN; (c) RF; and (d) SVM models.
Figure 5. Predicted and actual ROP along the S-shape section for the testing data (1717 data points) using (a) ANFIS; (b) FN; (c) RF; and (d) SVM models.
Sensors 20 03506 g005
Figure 6. Predicted and actual ROP along the S-shape section for the validation data (2500 unseen data points) using (a) ANFIS; (b) FN; (c) RF; (d) SVM; (e) Bingham correlation; (f) Maurer correlation; and (g) Bourgoyne–Young correlation.
Figure 6. Predicted and actual ROP along the S-shape section for the validation data (2500 unseen data points) using (a) ANFIS; (b) FN; (c) RF; (d) SVM; (e) Bingham correlation; (f) Maurer correlation; and (g) Bourgoyne–Young correlation.
Sensors 20 03506 g006
Figure 7. Comparison of the AAPE for ROP estimated using the ANFIS, FN, RF, and SVM models, as well as Bingham correlation, Maurer correlation, and Bourgoyne–Young correlation for the validation data.
Figure 7. Comparison of the AAPE for ROP estimated using the ANFIS, FN, RF, and SVM models, as well as Bingham correlation, Maurer correlation, and Bourgoyne–Young correlation for the validation data.
Sensors 20 03506 g007
Table 1. Well A and Well B training data statistical parameters (model training data).
Table 1. Well A and Well B training data statistical parameters (model training data).
ParameterQ (gpm)DSR (rpm)SPP (psi)T (klbf-ft)WOB (klbf)ROP (ft/h)
Minimum559999124.115.146.34
Maximum993159249015.957.264.5
Range43460157911.852.158.1
Standard deviation94.010.83472.2611.415.0
Skewness−0.596−0.976−0.3490.5240.0450.184
Kurtosis−0.6922.643−0.952−0.140−0.941−1.193
Table 2. The optimized artificial intelligence models.
Table 2. The optimized artificial intelligence models.
ANFIS
Cluster radius0.1
Ep-size300
FN
Training methodNonlinear function without iteration terms
Training function typeForward selection method
RF
Maximum depth29
Maximum featureSqrt
SVM
KernelRadial basis function
GammaScale
No. of iterations1000
C100
Table 3. Effect of including the hole inclination as an input parameter using the optimized fuzzy inference system (ANFIS).
Table 3. Effect of including the hole inclination as an input parameter using the optimized fuzzy inference system (ANFIS).
TrainingTestingValidation
RAAPE %RAAPE %RAAPE %
Without inclination0.979.050.979.570.979.06
With inclination0.9412.850.9213.810.1651.29

Share and Cite

MDPI and ACS Style

Elkatatny, S. Real-Time Prediction of Rate of Penetration in S-Shape Well Profile Using Artificial Intelligence Models. Sensors 2020, 20, 3506. https://doi.org/10.3390/s20123506

AMA Style

Elkatatny S. Real-Time Prediction of Rate of Penetration in S-Shape Well Profile Using Artificial Intelligence Models. Sensors. 2020; 20(12):3506. https://doi.org/10.3390/s20123506

Chicago/Turabian Style

Elkatatny, Salaheldin. 2020. "Real-Time Prediction of Rate of Penetration in S-Shape Well Profile Using Artificial Intelligence Models" Sensors 20, no. 12: 3506. https://doi.org/10.3390/s20123506

APA Style

Elkatatny, S. (2020). Real-Time Prediction of Rate of Penetration in S-Shape Well Profile Using Artificial Intelligence Models. Sensors, 20(12), 3506. https://doi.org/10.3390/s20123506

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