A Novel Neural Network Framework For The Prediction of Drilling Rate of Penetration
A Novel Neural Network Framework For The Prediction of Drilling Rate of Penetration
A Novel Neural Network Framework For The Prediction of Drilling Rate of Penetration
net/publication/316963557
CITATIONS READS
0 71
2 authors:
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Saurabh Tewari on 16 May 2017.
MODEL DEVELOPMENT
Bourgoyne and Young ROP model
Bourgoyne and Young have proposed a simple ROP prediction model for roller cone bit. The
drilling ROP can be given as:
8
ROP Exp a1 ai xi (1)
i 2
where xi are various functional drilling parameters which effect the ROP and a1 to a8 are
constants (Table 1 & 2). Drilling behavior of a given formation type is modelled through the
selection of constants from a1 to a8 in equation 1. BY determined the value of constants a1 to a8
through the MRA of field data.
Table 1. Functional relationships associated with ROP model [3].
Serial Equations Effect description
No
1 nX1 n is number of data point
2 X 10, 000.0 De Depth
2
3 1 Formation Compaction
X De 0.69 G 9.0
3 p 2
4 Differential pressure
X D e G P
4 p c
5 W W W Bit weight and bit diameter
X ln 4.0
5
d d t d t
6 X ln N 100 Rotary speed
6
7 X H Tooth wear
7
8 X pq 350 d Fluid properties
8 n
The first constant a1 symbolizes the effects of formation strength, bit types, mud types and solid
content which are not included in ROP model. The terms a2x2 and a3x3 represent the effect of
compaction on the drilling rate. The effect of overbalance on penetration rate is represented by
the a4x4 term. The effect of bit weight and rotary speed is included in the terms a5x5 and a6x6.
Remaining terms a7x7 and represents the effect of tooth wear and characterizes the effects of bit
hydraulics on ROP [3]. The functional relationships associated with equation 1 are listed in
Table 1.The coefficients a1 to a8 must be determined for each formation because of their
APCEC17 INTERNATIONAL CONFERENCE ON ADVANCES IN PETROLEUM CHEMICAL
& ENERGY CHALLENGES
dependency toward local drilling conditions [3]. BY has also suggested upper and lower for each
of eight coefficients after testing the different formations as shown in Table 2 [3].
Table 2. Recommended upper and lower bound for Bourgoyne and Young ROP model [3].
Coefficients a1 a2 a3 a4 a5 a6 a7 a8
Lower 0.5 0.0000001 0.0000001 0.000001 0.5 0.4 0.3 0.3
Bound
Upper 1.9 0.0005 0.0009 0.001 0.2 0.1 1.5 0.6
Bound
BY has determined the values of eight coefficients by utilizing the concept of MRA and
recommended that minimum thirty data points are required [3]. MRA gives results out of
recommended bounds in certain conditions results in zero or negative values also. Zero or
negative values of coefficients are practically worthless in nature. This method is also limited to
the recommended ranges for drilling parameters and results are not highly reliable in nature.
Therefore, new technique is required which can predict ROP in highly accurate way.
All the xi which are misclassified by previous iteration is the learning rate and is the error
associated with input. Many nodes combined together to form layers and these layers combine to
form network. The layers are divided into three types i.e. input layer, hidden layer and output
layer. Number of neurons and layers depends upon the complexity of problems. For training
purpose back propagation algorithm is most suitable for multilayer perceptron neural network.
The data is processed through input layer, then to hidden layer and finally to output layer. The
APCEC17 INTERNATIONAL CONFERENCE ON ADVANCES IN PETROLEUM CHEMICAL
& ENERGY CHALLENGES
error between the output and target value are measures and passed backward by updating
individual weights and bias neurons layers.
where XMin and XMax are minimum and maximum values of the corresponding datasets. Outlier
analysis is also done to detect and remove unwanted noisy data points. In this research work, for
the estimation of ROP, ANN models are developed with eight inputs (BY parameters) and ROP
as a single output. Optimum number of neurons and hidden layers essential for the prediction of
ROP were determined by minimum RMSE value. The different ANN network architectures were
developed and tested on available data and recorded in Table 3. Early stopping was implemented
to avoid the problem of over fitting and improvise the generalization performance of ANN. The
learning and generalization ability of ANN model is evaluated on the basis of three statistical
parameters i.e.
APCEC17 INTERNATIONAL CONFERENCE ON ADVANCES IN PETROLEUM CHEMICAL
& ENERGY CHALLENGES
N
1 Tm OP
MAPE
N
I 0 Tm
*100 (5)
N
1 2
RMSE
N
T
i 1
m Op (6)
N
2
T m Op
R2 1 i 1
2 (7)
N N
1
Tm N Tm i
i 1 i 1 i
where Tm and Op is the measured and predicted value of ROP. The 8-20-1 feed forward ANN
architecture is found to be optimum for prediction of ROP with least RMSE (0.39964) and
highest coefficient of correlation R2 (0.86174). The comparison between BY and ANN model is
done in table 5 using three statistical parameters. The test results clearly indicate that ANN has
better prediction performance as compared with conventional Bourgoyne and Young model.
Table 3. Comparison between various neural network models.
Table 5: Sensitivity analysis of conventional ROP models and ANN based ROP estimator
models.
Figure 3. Optimum Feedforward Back propagation ANN developed for ROP prediction.
Figure 4. Plot of predicted versus measured Bourgoyne and Young/ANN ROP model value.
APCEC17 INTERNATIONAL CONFERENCE ON ADVANCES IN PETROLEUM CHEMICAL
& ENERGY CHALLENGES
CONCLUSIONS
A real time ANN estimator model is successfully developed and tested on Kinabalu field data for
the prediction of drilling rate of penetration. The results indicate high degree of estimation
accuracy as compared with conventional Bourgoyne and Young ROP model. The proposed ANN
model can be adjusted according to the changing wellbore conditions and availability of data in
real time operation. The drilling parameters such as weight on bit and bit rotation speed are
adjusted to achieve most effective rate of penetration for drilling operations. Real time analysis
of ROP helps the drilling engineer for proper selection of various drilling and hydraulic
parameters easily. The results clearly show that ANN based ROP model is more efficient for
drilling oil and gas fields and can be used for proper planning and successful execution of real
time drilling operations.
APCEC17 INTERNATIONAL CONFERENCE ON ADVANCES IN PETROLEUM CHEMICAL
& ENERGY CHALLENGES
ACKNOWLEDGEMENT
We thank Rajiv Gandhi Institute of Petroleum Technology for providing us financial
support and research infrastructure.
NOMENCLETURE
a8 = hydraulic exponent
b = bearing constant
B = fractional bearing wear
d = bit diameters, in
D= well depth, ft
[5]. J. W. Graham and N.L. Muench, “ Analytical determination of optimum bit weight and rotary
speed combinations” in SPE-AIME 34th Annual Fall Meeting, SPE 1349-G, Dallas Oct. 4-7,
1959
[6]. S. Irawan, A.M.A. Rahman, and S.Q. Tunio, “Optimization of weight on bit during drilling
operation based on rate of penetration model” in Res. J. of App.Sci.Vol4 (12) pp 1690-1695,
2012.
[7]. H. Rahimzadeh, M. Mostofi, A. Hashemi, and Salahshoor, “Comparison of the penetration rate
models using field data for one of the gas field in Persian Gulf area” in SPE International Oil and
Gas Conf. and Exhi., SPE 131253, Beijing, China, 2012.
[8]. J.W. Speer, “A method for determining Optimum drilling Technique” in Drill. and Prod. Prac. ,
API, 1958.