Next Article in Journal
Aeroelastic Analysis of a Coplanar Twin-Rotor Wind Turbine
Next Article in Special Issue
Estimation of Static Young’s Modulus for Sandstone Formation Using Artificial Neural Networks
Previous Article in Journal
Does Energy-Growth and Environment Quality Matter for Agriculture Sector in Pakistan or not? An Application of Cointegration Approach
Previous Article in Special Issue
A Novel Nonintrusive Load Monitoring Approach based on Linear-Chain Conditional Random Fields
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Data-Driven Framework to Predict the Rheological Properties of CaCl2 Brine-Based Drill-in Fluid Using Artificial Neural Network

by
Ahmed Gowida
1,
Salaheldin Elkatatny
1,*,
Emad Ramadan
2 and
Abdulazeez Abdulraheem
1
1
Petroleum department, King Fahd University of Petroleum & Minerals, Dhahran 31261 Box 5049, Saudi Arabia
2
Information & Computer Science Department, King Fahd University of Petroleum & Minerals, Dhahran 31261 Box 5049, Saudi Arabia
*
Author to whom correspondence should be addressed.
Energies 2019, 12(10), 1880; https://doi.org/10.3390/en12101880
Submission received: 18 April 2019 / Revised: 12 May 2019 / Accepted: 13 May 2019 / Published: 17 May 2019

Abstract

:
Calcium chloride brine-based drill-in fluid is commonly used within the reservoir section, as it is specially formulated to maximize drilling experience, and to protect the reservoir from being damaged. Monitoring the drilling fluid rheology including plastic viscosity, P V , apparent viscosity, A V , yield point, Y p , flow behavior index, n , and flow consistency index, k , has great importance in evaluating hole cleaning and optimizing drilling hydraulics. Therefore, it is very crucial for the mud rheology to be checked periodically during drilling, in order to control its persistent change. Such properties are often measured in the field twice a day, and in practice, this takes a long time (2–3 h for taking measurements and cleaning the instruments). However, mud weight, M W , and Marsh funnel viscosity, M F , are periodically measured every 15–20 min. The objective of this study is to develop new models using artificial neural network, ANN, to predict the rheological properties of calcium chloride brine-based mud using M W and M F measurements then extract empirical correlations in a white-box mode to predict these properties based on M W and M F . Field measurements, 515 points, representing actual mud samples, were collected to build the proposed ANN models. The optimized parameters of these models resulted in highly accurate results indicated by a high correlation coefficient, R, between the predicted and measured values, which exceeded 0.97, with an average absolute percentage error, AAPE, that did not exceed 6.1%. Accordingly, the developed models are very useful for monitoring the mud rheology to optimize the drilling operation and avoid many problems such as hole cleaning issues, pipe sticking and loss of circulation.

1. Introduction

Drilling Fluids are considered a key element in the drilling operation. Conventional drilling fluids are water-based, oil-based or synthetic-based fluid systems, which are used in the drilling process to give the best performance under certain temperatures and pressures experienced downhole [1]. Drilling the section from the sea-bed/land to the top of the reservoir is different, regarding the economic value of the final project, compared to the reservoir section. As in the top sections, the concerns are to seal the permeable formations, and help sustain the wellbore stability.
Special measures will be taken into consideration while drilling the reservoir section to avoid damaging the reservoir and plugging the reservoir pores. For that target, special drilling fluids are used, called reservoir drill-in fluids (RDFs), which are specially formulated to maximize drilling experience and protect the reservoir from being damaged until the completion process is proceeded [2].
There are many types of RDFs with different chemical compositions, but the concern of this study is about the clear brine-based mud which is often used within completions, as the presence of solids is a major contributor to formation damage [3]. However, when used as drilling fluid, the solids-free nature of brine operationally improves the rate of penetration (ROP), stabilization of sensitive formations, density, and abrasion or friction [4]. Clear brine fluids properties are easier to maintain than conventional solids-laden fluid systems, so that when properly run, clear systems require very little maintenance, because many functional issues are inherently solved by the dissolved salts. Clear brine fluids also allow for drill site cost reductions because of the ability to reuse the fluid [5].
Brine fluids can be prepared with one salt or a combination of salts. All salts provide unique properties to the base fluid. Saturated brines fluids provide excellent inhibitive properties and lubricity, as compared to conventional aqueous fluids. With optimal heat transference characteristics, they can greatly improve bit life, and increase the rate of penetration in hard rock drilling. Among the different salts used for clear-brine systems, calcium chloride has been selected, as it is considered one of the most economic brine systems, with its broad range of densities (from 9.0 to 11.6 ppg), availability, low cost, and its ability to reduce the water activity of the fluid [6].

1.1. Drilling Fluid Rheology

Drilling fluid rheology plays a key role in optimizing drilling performance [7]. These properties significantly affect the efficiency of the hole cleaning and the drilling rate [8], which are critical factors controlling the performance of drilling operation [9]. These rheological properties include mud density to provide the control on the formation pressure, while P V ,   Y P ,   A V ,   n and k are used for controlling hole cleaning and optimizing the drilling performance [10]. Plastic viscosity of the drilling fluid is crucial for optimizing the drilling operation [11]. It is an indication of the solid content in the drilling fluid which may negatively affect the drilling performance when it exceeds critical limits, and can cause many problems like pipe sticking and decreasing the rate of penetration [12]. Yield point can be simply stated as the attractive forces among colloidal particles in drilling fluid [11]. The optimization of yield point is central to controlling the efficiency of hole cleaning [10,13]. Moreover, apparent viscosity is considered a key factor in the optimization of mud hydraulics while drilling [8]. In addition, the parameters k and n can be used for evaluating hole cleaning during the drilling operation [14].
The rheological properties can be measured in the laboratory using mud balance and viscometer. The mud balance is used to measure the mud weight while the rheometer is used to measure ( P V ,   Y P ,   A V ,   n and k ). However, this process takes a relatively long time (2–3 h for taking measurements and cleaning the instruments) which makes it difficult to be performed periodically and practically in the field. Therefore, it is taken as a common procedure that only density and Marsh funnel viscosity are measured periodically every 15–20 min, using mud balance and Marsh funnel devices. On the other hand, a complete mud test (including all the drilling fluid properties), using the mud balance and viscometer, is performed twice a day. Marsh funnel viscosity provides an indication of the changes in the rheology of the drilling fluid. This funnel was first introduced by Marsh [15]. This tool is cheap and takes a short time, so it can be utilized to give field measurements frequently and estimate some parameters like yield stress [16]. Based on the literature, there are two models developed to predict the drilling fluid viscosity from mud density and Marsh funnel measurements. These two measurements were used as inputs to calculate the effective viscosity of the drilling fluid as stated by Pitt [17] in Equation (1). Then a modification on the previous model was introduced by Almahdawi et. al. [18], who figured out that changing the value of the constant to 28 in Equation (2) instead of 25 presented by Pitt [17], is more effective give more accurate results, as compared to Equation (1).
A V = D ( T 25 )
A V = D ( T 28 )
where A V is the apparent viscosity in cP, D is the fluid density in g / c m 3 , T is the Marsh funnel time in sec.
Several mathematical models have been mentioned in the literature for estimating the fluid rheological properties using Marsh funnel devices. Some of them suggested using the temporal variations in the fluid height in the funnel to determine different rheological parameters such as P V ,   Y P   and   A V [19,20,21,22]. They introduced a methodology to determine the shear rate and the shear stress on the walls of the Marsh funnel from the measured discharged fluid volume of the Marsh funnel at different points. Then several rheological parameters have been related to the obtained shear rate and shear stress. Abdulrahman et al. [23] investigated different water-based drilling fluids using the Marsh funnel and showed that P V   and   A V can be estimated using consistency plots and the methodology described in [19]. However, these models showed considerable discrepancies between the results obtained from the Marsh funnel and the standard viscometers. Other studies tried to model the fluid volume flow in the Marsh funnel with higher order polynomial functions, rather than the simplified functions used in the previously mentioned studies [24,25]. This attempt was to simulate the fluid temporal height in the Marsh funnel more properly, and to get closer results of rheological parameters to those obtained from the standard viscometers.
The objective of this work is to develop new models using artificial neural networks, ANN, to predict the rheological properties of the CaCl2 brine-based drilling fluid depending on frequent measurements of M W and M F . The real-time measurements of these parameters are very helpful for identifying the efficiency of the hole cleaning, optimizing the drilling fluid hydraulics, equivalent circulating density calculations and swab and surge pressure determination.

1.2. Artificial Neural Network (ANN)

Artificial intelligence, (AI), can be simply defined as the computer science branch for creating intelligent machines [26] to exhibit human brains to make predictions and help take the right decisions for the future scenarios [27]. Recently, different AI methods such as fuzzy logic, FL, support vector machine, SVM, genetic algorithm and artificial neural network, ANN, have been applied in petroleum engineering, and specifically in the field of drilling fluid engineering. Some of these applications include fluid flow patterns prediction in wellbore annulus [28], stuck pipe prediction [29], drilling hydraulics optimization [30], frictional pressure loss estimation [31], hole cleaning and prediction of cutting concentration [32], estimation of the static Poisson’s ratio from log data [33].
ANN is one of the most common AI techniques which has the ability to deal with different engineering problems with high complexity that exceed the computational capability of classical mathematics and procedures [34]. It is based on analogy with biological neural networks to simulate the performance of the human biological neural system [35]. The elementary units for ANN are neurons [36]. The structure of the ANN consists of three main types of layers. The first one is for the input parameters. The second one is called hidden layers, which include the neurons assigned with the transfer functions between the inputs and the outputs. The third type is for the outputs. These layers, with the suitable training algorithm, describe the nature of the problem [37]. The performance of the network is controlled by key parameters including the number of neurons, weights and biases [38]. To optimize the weight and biases, the network is trained using different algorithms to achieve the lowest possible error. Among these algorithms is Levenberg-Marquardt (LM), which is an iterative, curve fitting algorithm. This algorithm proved its outstanding performance in solving non-linear least-squares problems [39].
There are many ANN applications of ANN in the field of drilling fluid in the last few years. Some of these researches are the prediction of filtration volume and mud cake permeability of water-based mud (WBM) [40], drill cutting settling velocity prediction [41], prediction of differential pipe sticking [42], lost circulation prediction [43], hole cleaning efficiency of foam fluid [44], rheological properties of invert emulsion mud [45], invert emulsion mud rheology [46] and spud mud rheology prediction [47], generating geomechanical well logs [48], prediction of oil PVT properties [49].
Based on the literature, more than 50 percent of the applications in the drilling fluid area used ANN for the predictions and got high accurate results. Accordingly, ANN has been selected for building the proposed models in this study [26].

2. Methodology

2.1. Data Description

A typical sample of the data for CaCL2 brine-based drill-in fluid (515 field data for actual mud samples) is listed in Table 1, including ( M W ,   M F ,   P V ,   and   Y p ) . The drilling fluid samples are collected after the mud was cleaned from the cuttings by using the shale shaker M W and M F are measured in the field using a mud balance and Marsh funnel, respectively. The rheometer is used for measuring the rheology of the mud, namely P V ,   and   Y p at atmospheric pressure and 120 o F . The collected data have a wide range as follows: M W ranges from 43 to 119 I b / f t 3 , M F ranges from 26 to 135 s / q u a r t , P V ranges from 10 to 54 cP, and Y P ranges from 8 to 41 I b / 100   f t 2 . Figure 1 shows that M W has R of 0.36 and 0.76 with Y P and P V respectively while M F has R of 0.86 with Y P and 0.36 with P V .
For better prediction using AI models, data should be analyzed and filtered [50]. Therefore, the selected data have been cleaned from any noise and false values for higher representation quality. The filtration process included eliminating all the values that cannot be representative, like negative values. Finally removing the outliers that show significant deviation from the other values of a variable, the outliers were removed using a box and whisker plot, in which top whisker represents the upper limit of the data, and the bottom whisker represents the lower limit of the data, then any value beyond these limits is considered an outlier and removed [51]. These limits are determined by dividing the data into four equal divisions (quartiles) along with using the minimum, maximum, mean and median parameters [52] obtained from the statistical analysis of the data listed in Table 2.

2.2. Development of ANN Models

The collected data were used to calculate R 600 and R 300 (rheometer readings at 600 and 300 rpm, respectively) using Equations (3) and (4). These two parameters are very crucial for identifying fluid properties and flow regimes. Then, the apparent viscosity, A V , flow behavior index, n , and flow consistency index, k , are calculated using Equations (5)–(7) respectively.
R 600 = P V + R 300
R 300 = P V + Y p
A V = R 600 2
n = 3.32 × log ( R 600 R 300 )
k = R 600 1022 n
For all the upcoming developed models, different scenarios have been performed to optimize the ANN variables to reach the highest accuracy with the lowest possible error for prediction using different combinations of the available options of the ANN variables. The optimized parameters obtained from the tuning process of these parameters are summarized in Table 3. The chosen architecture for the developed models includes three layers:
Input layer: It contains input features which are M W and M F .
(One) Hidden layer: It contains the optimized number of neurons which was found to be 20 neurons.
Output layer: It contains the output parameters, which are ( P V ,   Y P ,   A V ,   n and k individually).
The network was trained using the Levenberg-Marquardt (LM) algorithm to get the optimized weights and biases. The neurons are arranged to be trained using a learning rate of 0.12. Activation function of the tan-sigmoidal type (tansig) was assigned between the input and hidden layers while the pure-linear function was assigned between the hidden and output layers. Figure 2 shows a typical schematic of the architecture of the developed ANN model.

3. Results and Discussion

3.1. Yield Point ( Y P ) Model

An ANN-Based model was developed using M W and M F as inputs to predict the Y P values. The obtained data were divided into ratios 70:30 for training and testing the model, respectively. Figure 3 shows the high match between the measured and the predicted Y P values from the developed ANN model in terms of R of 0.97 and AAPE of 3.9%. Thereafter, a new correlation has been developed using the ANN model to predict Y P based on M W and M F . First, the inputs should be normalized using Equations (8) and (9) to substitute the values M W n and M F n in Equation (10); where M W n refers to the first normalized input, and M F n refers to the second normalized input.
M W n = 0.036 ( M F 64 ) + 1
M F n = 0.133 ( M F 26 ) + 1
Then, the normalized value Y P n is calculated using Equation (9) with its optimized coefficients listed in Table 4.
Y P n = [ i = 1 N w 2 , i ( 2 1 + exp ( 2 ( M W n × w 1 i , 1 + M F n × w 1 i , 2 + b 1 , i ) ) 1 ) ] + b 2
where i is the index of the neuron in the hidden layer, N is the optimized number of neurons for only one hidden layer, which is found to be 20 , w 1 is the weight vector linking the input and the hidden layer, w 2 is the weight vector linking the hidden and output layer, b 1 is the biases vector for the input layer, and b 2 for the output layer.
Finally, the required Y P value can be obtained by denormalizing Y P n using Equation (11).
Y p = 16.5 ( Y p n + 1 ) + 8

3.2. Apparent Viscosity ( A V ) Model

Similarly, A V was predicted using ANN, based on M W and M F . The model was trained using 70% of the available data, while 30% of the data were used for testing the model. Figure 4 shows the high R between the predicted and the measured A V values, which is 0.99 with AAPE of 3.2%. Afterward, a new correlation for predicting A V was extracted from the developed ANN model. To use this correlation, the inputs should be normalized at first using Equations (12) and (13) to substitute the values M W n and M F n in Equation (14); where M W n refers to the first normalized input, and M F n refers to the second normalized input.
M W n = 0.036 ( M F 64 ) + 1
M F n = 0.053 ( M F 35 ) + 1
Then, the normalized value A V n is calculated using Equation (14) with its optimized coefficients listed in Table 5.
A V n = [ i = 1 N w 2 , i ( 2 1 + exp ( 2 ( M W n × w 1 i , 1 + M F n × w 1 i , 2 + b 1 , i ) ) 1 ) ] + b 2
Finally, A V can be predicted by denormalizing A V n using Equation (15).
A V = 27 ( A V n + 1 ) + 19

3.3. Plastic Viscosity ( P V ) Model

For P V , another ANN model was developed based on M W and M F . For building the model, the ratio of the training to testing points is 70:30. The model gave high accurate results indicated by a high R of 0.98 between the predicted and the measured P V values and maximum AAPE of 6.1% as shown in Figure 5. A new correlation has been extracted from the model to predict P V without the need to run the ANN model.
First, the inputs should be normalized using Equations (16) and (17) to substitute the values M W n and M F n in Equation (18); where M W n refers to the first normalized input, and M F n refers to the second normalized input.
M W n = 0.037 ( M F 64 ) + 1
M F n = 0.105 ( M F 35 ) + 1
Then, the normalized value P V n is calculated using Equation (18) with its optimized coefficients listed in Table 6.
P V n = [ i = 1 N w 2 , i ( 2 1 + exp ( 2 ( M W n × w 1 i , 1 + M F n × w 1 i , 2 + b 1 , i ) ) 1 ) ] + b 2
Finally, P V can be predicted by denormalizing P V n using Equation (19).
P V = 22 ( P V n + 1 ) + 10

3.4. Prediction Power Law Model Parameters (n and k)

Following the same procedure, another two models have been developed using ANN to predict n and k based on M W and M F . For the prediction of n, the R between the measured and the predicted values was 0.98 with AAPE of 2.4% as shown in Figure 6. While for the prediction of k, the R was 0.99 with AAPE of 3.6%, as indicated in Figure 7. Then new correlations for estimating n and k were extracted from the developed ANN models.
In the beginning, the inputs should be normalized using Equations (20) and (21) for the correlation of n and Equations (22) and (23) for the correlation of k in order to substitute the values M W n and M F n in Equations (24) and (25); where M W n refers to the first normalized input and M F n refers to the second normalized input.
For the Correlation of Parameter n:
M W n = 0.026 ( M F 43 ) + 1
M F n = 0.018 ( M F 26 ) + 1
For the Correlation of Parameter k:
M W n = 0.036 ( M F 64 ) + 1
M F n = 0.022 ( M F 30 ) + 1
Subsequently, the normalized values n n and k n can be estimated using Equations (24) and (25), respectively, with their optimized coefficients listed in Table 7 and Table 8, respectively.
n n = [ i = 1 N w 2 , i ( 2 1 + exp ( 2 ( M W n × w 1 i , 1 + M F n × w 1 i , 2 + b 1 , i ) ) 1 ) ] + b 2
k n = [ i = 1 N w 2 , i ( 2 1 + exp ( 2 ( M W n × w 1 i , 1 + M F n × w 1 i , 2 + b 1 , i ) ) 1 ) ] + b 2
Eventually, the predicted values of n and k can be estimated using Equations (26) and (27), respectively.
n = 0.244 ( n n + 1 ) + 0.263
k = 2.78 ( k n + 1 ) + 0.731

3.5. Validation of the Apparent Viscosity ( A V ) Model vs the Models in the Literature

As mentioned in the introduction, Pitt [17] introduced a numerical model to calculate the apparent viscosity using Equation (1). After using the collected data, the results obtained using Equation (1) showed a coefficient of determination R 2 of 0.5 and AAPE of 32.2%. Also, Almahdawi et al. [18] concluded that Equation (2) using the constant 28 is more appropriate than 25 in Equation (1), and the results obtained by applying Equation (2) to estimate A V using M F readings give R 2 of 0.5 and AAPE of 23.6%, as shown in Figure 8. However, the developed correlation for the ANN model gives highly accurate results, as shown in Figure 9 with R 2 of 0.98, and AAPE does not exceed 3.2%.

4. The Value of Predicting the Drilling Fluid Rheology in Real-Time

For drilling optimization, it is very important to have periodic monitoring of the parameters affecting the drilling process. Mud system design and hole cleaning processes are affected by the pressure losses within the system which rely on the properties of the drilling fluid used, and the efficiency of the cuttings removal from the hole. Pressure losses can be obtained once the parameters of the Bingham model Y P and A V and power low model (n and k) are obtained. Annular pressure losses can be calculated by Equation (28) based on the real-time values of ( Y P and A V ), which can be obtained from the developed ANN model.
Also, equivalent circulating density, ECD, can be calculated from Equation (29) [53] using the obtained pressure loss value, so that surge and swab pressures can be determined to help predict critical drilling problems such as pipe sticking and well control issues [54].
Δ p = ( P V × v 1000 ( d 2 d 1 ) 2 + Y p 200 ( d 2 d 1 ) ) L
E C D = M W + Δ p 0.052 × h
where Δ p is the annular pressure loss (in psi), P V is the predicted plastic viscosity ( i n   c P ), Y P is the predicted yield point ( i n   I b / 100 f t 2 ), v is the average annular velocity (in ft/s), d 1 is the inside diameter of the hole or casing, (in inches), d 2 is the outside diameter of the drill pipe, (inches), L is the drill pipe, or drill collar length (in ft), M W is the mud density (in ppg), h is the hole depth (in ft), and ECD is the equivalent circulation density (in ppg).
Accordingly, the ability of the prediction of the rheological properties in real-time can help avoid many problems during drilling with early detection of these problems by identifying the anomaly in normal behavior trends. This will optimize the drilling operation and save money by minimizing the drilling time.

5. Conclusions

In this work, new models have been developed using ANN to predict the rheological properties of CaCL2 brine-based drill-in fluid in a real-time (15–20 min) including ( P V ,   Y P ,   A V ,   n and k ) using 515 field data measurements of M W and M F in ratios 70:30 for training and validating the ANN models respectively. Accordingly, the following conclusions can be drawn:
(1)
The new ANN models can predict the rheological parameters ( P V ,   Y p ,   A V ,   n ,   a n d   k ) in real time based on M W and M F with high accuracy (R was greater than 0.97 and AAPE was less than 6.1%).
(2)
The optimization process for the ANN models showed that the optimized parameters yielding the highest accuracy and the lowest error were 20 neurons for only one hidden layer, the Levenberg-Marquardt algorithm of learning rate 0.12. The activation function linking the input and hidden layers was the tan-sigmoidal function, while a linear function was used for linking the hidden and output layers.
(3)
The extracted correlations from the developed ANN models provide the ability to estimate the rheological properties of CaCL2 brine-based mud directly without the need to run the models.
(4)
These models are very helpful in the calculations of rig hydraulics, surge and swab pressures, and ECD.
(5)
The developed correlations can help in predicting several drilling problems by providing the ability for real-time monitoring of the hole cleaning performance, and detecting any abnormal changes in the normal trends to avoid interrupting problems like sticking. As a result, this will save on the drilling cost, and it optimizes the drilling operation.

Author Contributions

Conceptualization, S.E. and A.G.; methodology, E.R., A.A.; validation, S.E., A.G.; formal analysis, A.G.; investigation, A.A.; resources, S.E.; data curation, A.G.; writing—original draft preparation, S.E. and A.G.; writing—review and editing, A.G., A.A.; visualization, A.A., E.R.; supervision, S.E., A.A.

Funding

This research received no external funding.

Acknowledgments

The authors wish to acknowledge King Fahd University of Petroleum and Minerals (KFUPM) for utilizing the various facilities in carrying out this research. Many thanks are due to the anonymous referees for their detailed and helpful comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Caenn, R.; Darley, H.C.H.; Gray, G.R. Composition and Properties of Drilling and Completion Fluids, 6th ed.; Gulf Professional Publishing: Waltham, MA, USA, 2011. [Google Scholar]
  2. Knox, D.A.; Bradshaw, R.J.; Svoboda, C.F.; Hodge, R.M.; Wolf, N.O.; Hudson, C.E. Reservoir Drilling Fluids-Designing for Challenging Wells in a Remote Location. In Proceedings of the SPE Annual Technical Conference and Exhibition, Dallas, TX, USA, 9–12 October 2005. [Google Scholar]
  3. Conners, J.H. Use of Clear Brine Completion Fluids as Drill-In Fluids. In Proceedings of the SPE Annual Technical Conference and Exhibition, Las Vegas, NV, USA, 23–26 September 1979. [Google Scholar]
  4. Doty, P.A. Clear Brine Drilling Fluids: A Study of Penetration Rates, Formation Damage, and Wellbore Stability in Full-Scale Drilling Tests. SPE Drill. Eng. 1986, 1, 17–30. [Google Scholar] [CrossRef]
  5. Ballantine, W.T. Drill-site Cost Savings Through Waste Management. In Proceedings of the SPE 68th Annual Technical Conference and Exhibition of the Society of Petroleum Engineers, Houston, TX, USA, 3–6 October 1993. [Google Scholar]
  6. Redburn, M.; Heath, G. Improved Fluid Characteristics with Clear Calcium Chloride Brine Drilling Fluid. In Proceedings of the Offshore Mediterranean Conference and Exhibition, Ravenna, Italy, 29–31 March 2017. [Google Scholar]
  7. Power, D.; Zamora, M. Drilling fluid yield stress: Measurement techniques for improved understanding of critical drilling fluid parameters. In Proceedings of the AADE Technical Conference, Houston, TX, USA, 1–3 April 2003. [Google Scholar]
  8. Guo, B.; Liu, G. Mud Hydraulics Fundamentals. In Applied Drilling Circulation Systems; Gulf Professional Publishing: Houston, TX, USA, 2011; pp. 19–59. [Google Scholar]
  9. Mitchell, R.F.; Miska, S.Z. Fundamentals of Drilling Engineering; Society of Petroleum Engineers: Houston, TX, USA, 2011; pp. 87–90. [Google Scholar]
  10. Feng, Y.; Gray, K.E. Review of fundamental studies on lost circulation and wellbore strengthening. J. Pet. Sci. Eng. 2017, 152, 511–522. [Google Scholar] [CrossRef]
  11. Weikey, Y.; Sinha, S.L.; Dewangan, S.K. Role of additives and elevated temperature on rheology of water-based drilling fluid: A review paper. Int. J. Fluid Mech. Res. 2018, 45, 37–49. [Google Scholar] [CrossRef]
  12. Paiaman, A.M.; Al-Askari, M.K.G.; Salmani, B.; Al-Anazi, B.D.; Masihi, M. Effect of drilling fluid properties on the rate of penetration. NAFTA 2009, 60, 129–134. [Google Scholar]
  13. Saasen, A.; Løklingholm, G. The Effect of Drilling Fluid Rheological Properties on Hole Cleaning. In Proceedings of the IADC/SPE Drilling Conference, Dallas, TX, USA, 26–28 February 2002. [Google Scholar]
  14. Robinson, l.; Morgan, M. Effect of hole cleaning on drilling rate performance. In Proceedings of the AADE Drilling Fluid Conference Held at the Radisson, Houston, TX, USA, 6–7 April 2004. [Google Scholar]
  15. Marsh, H. Properties and treatment of rotary mud. Trans. AIME 1931, 92, 234–251. [Google Scholar] [CrossRef]
  16. Balhoff, M.T.; Lake, L.W.; Bommer, P.M.; Lewis, R.E.; Weber, M.J.; Calderin, J.M. Rheological and yield stress measurements of non-Newtonian fluids using a Marsh funnel. J. Pet. Sci. Eng. 2011, 77, 393–402. [Google Scholar] [CrossRef]
  17. Pitt, M.J. The Marsh Funnel and Drilling Fluid Viscosity: A New Equation for Field Use. SPE Drill. Completion 2000, 15, 3–6. [Google Scholar] [CrossRef]
  18. Almahdawi, F.H.; Al-Yaseri, A.Z.; Jasim, N. Apparent viscosity direct from Marsh funnel test. Iraqi J. Chem. Pet. Eng. 2014, 15, 51–57. [Google Scholar]
  19. Guria, C.; Kumar, R.; Mishra, P. Rheological analysis of drilling fluid using Marsh Funnel. J. Pet. Sci. Eng. 2013, 105, 62–69. [Google Scholar] [CrossRef]
  20. Liu, K.; Qiu, Z.; Luo, Y.; Zhou, G. Measure rheology of drilling fluids with Marsh funnel viscometer. Drill. Fluid Completion I Fluid 2014, 31, 60–62. (In Chinese) [Google Scholar]
  21. Schoesser, B.; Thewes, M. Marsh Funnel testing for rheology analysis of bentonite slurries for Slurry Shields. In Proceedings of the ITA WTC 2015 Congress and 41st General Assembl, Dubrovnik, Croatia, 22–28 May 2015. [Google Scholar]
  22. Schoesser, B.; Thewes, M. Do we tap the full potential from Marsh funnel tests? Rheology testing for bentonite suspensions. Tunn. J. 2016. [Google Scholar] [CrossRef]
  23. Abdulrahman, H.A.; Jouda, A.S.; Mohammed, M.M.; Mohammed, M.M.; Elfadil, M.O. Calculation Rheological Properties of Water Base Mud Using Marsh Funnel. Ph.D. Thesis, Sudan University of Science and Technology, Khartoum State, Sudan, 2015. [Google Scholar]
  24. Sedaghat, A.; Omar, M.A.A.; Damrah, S.; Gaith, M. Mathematical modelling of the flow rate in a marsh funnel. J. Energy Technol. Res. 2016, 1, 1–12. [Google Scholar] [CrossRef]
  25. Sedaghat, A. A novel and robust model for determining rheological properties of Newtonian and non-Newtonian fluids in a marsh funnel. J. Pet. Sci. Eng. 2017, 156, 896–916. [Google Scholar] [CrossRef]
  26. Agwu, O.E.; Akpabio, J.U.; Alabi, S.B.; Dosunmu, A. Artificial intelligence techniques and their applications in drilling fluid engineering: A review. J. Pet. Sci. Eng. 2018, 167, 300–315. [Google Scholar] [CrossRef]
  27. Ali, J.K. Neural networks: A new tool for the petroleum industry? In Proceedings of the European Petroleum Computer Conference, Aberdeen, UK, 15–17 March 1994. [Google Scholar]
  28. Oladunni, O.O.; Trafalis, T.B. Single-phase fluid ow classification via learning models. Int. J. Gen. Syst. 2011, 40, 561–576. [Google Scholar] [CrossRef]
  29. Rostami, H.; Khaksar Manshad, A. A New Support Vector Machine and Artificial Neural Networks for Prediction of Stuck Pipe in Drilling of Oil Fields. J. Energy Resour. Technol. 2014, 136, 024502. [Google Scholar] [CrossRef]
  30. Wang, Y.; Salehi, S. Application of Real-Time Field Data to Optimize Drilling Hydraulics Using Neural Network Approach. J. Energy Resour. Technol. 2015, 137, 062903. [Google Scholar] [CrossRef]
  31. Shahdi, A.; Arabloo, M. Application of SVM algorithm for frictional pressure loss calculation of three phase flow in inclined annuli. J. Pet. Environ. Biotechnol. 2016, 5, 1. [Google Scholar] [CrossRef]
  32. Al-Azani, K.; Elkatatny, S.; Abdulraheem, A.; Mahmoud, M.; Ali, A. Prediction of Cutting Concentration in Horizontal and Deviated Wells 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]
  33. Elkatatny, S. Application of Artificial Intelligence Techniques to Estimate the Static Poisson’s Ratio Based on Wireline Log Data. J. Energy Resour. Technol. 2018, 140, 072905. [Google Scholar] [CrossRef]
  34. Shahin, M.A.; Jaksa, M.B.; Maier, H.R. Recent advances and future challenges for artificial neural systems in geotechnical engineering applications. Adv. Artif. Neural Syst. 2009, 5. [Google Scholar] [CrossRef]
  35. Omosebi, A.; Osisanya, S.; Chukwu, G.; Egbon, F. Annular pressure prediction during well control. In Proceedings of the Nigeria Annual International Conference and Exhibition, Lagos, Nigeria, 6–8 August 2012. [Google Scholar]
  36. Nakamoto, P. Neural Networks and Deep Learning: Deep Learning Explained to Your Granny a Visual Introduction for Beginners Who Want to Make Their Own Deep Learning Neural Network (Machine Learning); CreateSpace Independent Publishing Platform: Scotts Valley, CA, USA, 2017. [Google Scholar]
  37. Lippmann, R. An introduction to computing with neural nets. IEEE ASSP Mag. 1987, 4, 4–22. [Google Scholar] [CrossRef]
  38. Hinton, G.; Osindero, S.; Teh, Y. A fast learning algorithm for deep belief nets. Neural Comput. 2006, 18, 1527–1554. [Google Scholar] [CrossRef] [PubMed]
  39. Lourakis, M.I. A brief description of the Levenberg-Marquardt algorithm implemented by levmar. Found. Res. Technol. 2005, 4, 1–6. [Google Scholar]
  40. Jeirani, Z.; Mohebbi, A. Artificial neural networks approach for estimating the filtration properties of drilling fluids. J. Jpn. Pet. Inst. 2005, 49, 65–70. [Google Scholar] [CrossRef]
  41. Agwu, O.E.; Akpabio, J.U.; Alabi, S.B.; Dosunmu, A. Settling velocity of drill cuttings in drilling fluids: A review of experimental, numerical simulations and artificial intelligence studies. Powder Technol. 2018, 339, 728–746. [Google Scholar] [CrossRef]
  42. Siruvuri, C.; Nagarakanti, S.; Samuel, R. Stuck pipe prediction and avoidance: A convolutional neural network approach. In Proceedings of the IADC/SPE Drilling Conference, Miami, FL, USA, 21–23 February 2006. [Google Scholar]
  43. Moazzeni, A.R.; Nabaei, M.; Jegarluei, S.G. Prediction of lost circulation using virtual intelligence in one of Iranian oil fields. In Proceedings of the Annual SPE International Conference and Exhibition, Tinapa-Calabar, Nigeria, 31 July–7 August 2010. [Google Scholar]
  44. Rooki, R. Estimation of pressure loss of Herschel Bulkley drilling fl uids during horizontal annulus using artificial neural network. J. Dispers. Sci. Technol. 2015, 36, 161–169. [Google Scholar] [CrossRef]
  45. 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]
  46. Elkatatny, S. Real-time prediction of rheological parameters of KCL water-based drilling fluid using artificial neural networks. Arab. J. Sci. Eng. 2017, 42, 1655–1665. [Google Scholar] [CrossRef]
  47. Abdelgawad, K.; Elkatatny, S.; Moussa, T.; Mahmoud, M.; Patil, S. Real-Time Determination of Rheological Properties of Spud Drilling Fluids Using a Hybrid Artificial Intelligence Technique. J. Energy Resour. Technol. 2019, 141, 032908. [Google Scholar] [CrossRef]
  48. Parapuram, G.; Mokhtari, M.; Ben Hmida, J. An Artificially Intelligent Technique to Generate Synthetic Geomechanical Well Logs for the Bakken Formation. Energies 2018, 11, 680. [Google Scholar] [CrossRef]
  49. Elkatatny, S.; Moussa, T.; Abdulraheem, A.; Mahmoud, M. A Self-Adaptive Artificial Intelligence Technique to Predict Oil Pressure Volume Temperature Properties. Energies 2018, 11, 3490. [Google Scholar] [CrossRef]
  50. Hemphill, T.; Bern, P.; Rojas, J.; Ravi, K. Field validation of drillpipe rotation effects on equivalent circulating density. In Proceedings of the SPE annual technical conference and exhibition, Anaheim, CA, USA, 11–14 November 2007. [Google Scholar]
  51. Dawson, R. How Significant Is A Boxplot Outlier? J. Stat. Educ. 2011, 19. [Google Scholar] [CrossRef]
  52. Chandrasegar, T.; Vignesh, M.; Balaji, R. Data Analysis Using Box and Whisker Plot for Lung Cancer. In Proceedings of the Innovations in Power and Advanced Computing Technologies (i-PACT) Conference, Vellore, India, 21–22 April 2017. [Google Scholar]
  53. Guo, B.; Liu, G. Applied Drilling Circulation Systems: Hydraulics, Calculations, and Models; Gulf Professional Publishing: Houston, TX, USA, 2011. [Google Scholar]
  54. Burkhardt, J.A. Wellbore pressure surges produced by pipe movement. J. Pet. Technol. 1961, 13, 595–605. [Google Scholar] [CrossRef]
Figure 1. The relative importance of M W and M F with the rheological properties ( Y P and P V ) of CaCL2 brine-based mud in terms of the correlation coefficient, R.
Figure 1. The relative importance of M W and M F with the rheological properties ( Y P and P V ) of CaCL2 brine-based mud in terms of the correlation coefficient, R.
Energies 12 01880 g001
Figure 2. Typical schematic for the architecture of the developed ANN models.
Figure 2. Typical schematic for the architecture of the developed ANN models.
Energies 12 01880 g002
Figure 3. Measured Y P vs. Predicted Y P from the ANN model.
Figure 3. Measured Y P vs. Predicted Y P from the ANN model.
Energies 12 01880 g003
Figure 4. Measured A V vs. Predicted A V from the ANN model.
Figure 4. Measured A V vs. Predicted A V from the ANN model.
Energies 12 01880 g004
Figure 5. Measured P V vs. Predicted P V from the ANN model.
Figure 5. Measured P V vs. Predicted P V from the ANN model.
Energies 12 01880 g005
Figure 6. Measured n vs. Predicted n from the ANN model.
Figure 6. Measured n vs. Predicted n from the ANN model.
Energies 12 01880 g006
Figure 7. Measured k vs. Predicted k from the ANN model.
Figure 7. Measured k vs. Predicted k from the ANN model.
Energies 12 01880 g007
Figure 8. Prediction of A V using Pitt’s [17] and Almahdawi’s et al. [18] correlations.
Figure 8. Prediction of A V using Pitt’s [17] and Almahdawi’s et al. [18] correlations.
Energies 12 01880 g008
Figure 9. Prediction of A V using the developed ANN model (AAPE of 3.2%).
Figure 9. Prediction of A V using the developed ANN model (AAPE of 3.2%).
Energies 12 01880 g009
Table 1. A typical sample of the CaCl2 brine-based drilling fluid collected data.
Table 1. A typical sample of the CaCl2 brine-based drilling fluid collected data.
M W ,   I b / f t 3 M F ,   s / q u a r t P V ,   c P Y P ,   I b / 100   f t 2
7862821
8045822
88621218
73501120
65481121
75441220
Table 2. Statistical Analysis of the CaCl2 brine-based mud collected data.
Table 2. Statistical Analysis of the CaCl2 brine-based mud collected data.
Parameter M W ,   I b / f t 3 M F ,   s / q u a r t P V ,   c P Y P ,   I b / 100   f t 2
Min.4326108
Max.1191355441
Mean85.4556.3322.5925.88
Mode761094433
Range72501924
Skewness0.501.431.290.83
Table 3. Summary of the optimized parameters for the developed ANN models.
Table 3. Summary of the optimized parameters for the developed ANN models.
Neural Network ParameterTypes and Range
Training AlgorithmLevenberg Marquardt
Number of neurons20
Number of hidden layer(s)1
Learning rate0.12
The hidden layer transfer functionTan-sigmoidal
The outer layer transfer functionPure-linear
Training ratio70%
Testing ratio30%
Table 4. The optimized coefficients for estimating the normalized Y P n in Equation (10).
Table 4. The optimized coefficients for estimating the normalized Y P n in Equation (10).
Neuron IndexInput Layer WeightsHidden Layer WeightsInput Layer BiasesOutput Layer Bias
i w 1 i , 1 w 1 i , 2 w 2 , i b 1 , i b 2
1−4.2515.585−0.9505.530−0.508
2−0.709−6.857−0.1554.936-
32.6315.5390.168−4.752-
4−0.7436.4111.0053.899-
5−4.9865.903−0.9323.661-
6−5.203−0.250−1.0223.135-
74.8594.645−0.410−2.792-
8−1.185−6.1920.7212.879-
94.188−3.646−1.697−3.015-
103.238−5.0800.297−0.585-
110.708−7.849−0.3800.213-
12−4.893−9.220−0.428−1.489-
132.227−6.971−1.1732.051-
143.1015.504−1.0461.840-
15−6.059−1.558−0.030−3.063-
16−5.0203.702−0.902−3.873-
172.8925.503−0.2604.287-
18−0.736−5.6680.190−5.818-
194.290−4.5920.6395.571-
20−4.2904.570−0.686−6.252-
Table 5. The optimized coefficients for estimating the normalized A V n in Equation (14).
Table 5. The optimized coefficients for estimating the normalized A V n in Equation (14).
Neuron IndexInput Layer WeightsHidden Layer WeightsInput Layer BiasesOutput Layer Bias
i w 1 i , 1 w 1 i , 2 w 2 , i b 1 , i b 2
12.9607.136−0.950−5.6671.535
2−5.586−7.703−0.1552.751-
34.0375.3550.168−1.886-
4−3.3623.7431.0052.292-
56.295−5.066−0.932−8.110-
60.4067.091−1.0226.492-
7−10.26−8.654−0.410−0.995-
8−0.572−9.0220.721−0.909-
9−7.5654.329−1.6974.884-
10−4.2563.8550.2972.094-
116.4581.765−0.3801.786-
124.537−4.152−0.4281.324-
13−5.4103.103−1.173−2.553-
14−4.859−2.202−1.046−1.924-
15−7.1901.704−0.030−4.783-
162.196−5.993−0.9022.408-
17−0.5766.113−0.260−4.569-
18−2.889−4.7820.190−4.645-
19−3.799−7.5880.639−2.337-
20−3.8774.861−0.686−6.301-
Table 6. The optimized coefficients for estimating the normalized P V n in Equation (18).
Table 6. The optimized coefficients for estimating the normalized P V n in Equation (18).
Neuron IndexInput Layer WeightsHidden Layer WeightsInput Layer BiasesOutput Layer Bias
i w 1 i , 1 w 1 i , 2 w 2 , i b 1 , i b 2
1−3.7404.0011.9839.616−1.831
2−3.304−5.296−1.243−5.482-
3−11.57−3.7883.3807.537-
4−6.403−2.376−4.8334.301-
51.3087.156−0.3072.069-
60.457−12.031.1821.413-
7−3.68410.384−0.1413.236-
81.511−3.8871.414−3.601-
9−7.490−1.8481.7886.460-
10−5.945−6.028−0.9855.412-
112.211−1.365−1.087−1.030-
126.136−3.4091.0933.100-
13−0.450−2.7590.5601.993-
1415.10415.3360.6121.763-
158.4232.7740.5017.032-
16−6.361−2.459−0.544−1.495-
17−5.2525.0031.075−1.282-
18−4.470−4.5471.533−5.938-
19−3.4576.2101.761−2.930-
203.769−5.5431.0145.828-
Table 7. The optimized coefficients for estimating the normalized n n in Equation (24).
Table 7. The optimized coefficients for estimating the normalized n n in Equation (24).
Neuron IndexInput Layer WeightsHidden Layer WeightsInput Layer BiasesOutput Layer Bias
i w 1 i , 1 w 1 i , 2 w 2 , i b 1 , i b 2
111.343−1.862−0.798−9.530−0.867
22.093−7.313−1.517−7.827-
3−0.1977.842−0.9375.778-
46.518−3.3192.163−5.132-
54.9671.743−1.977−3.641-
6−5.643−4.7880.7892.501-
77.452−2.4911.194−1.205-
8−8.8731.1681.4950.426-
9−3.8210.042−8.5162.781-
10−4.8961.7596.4603.734-
11−6.355−8.0890.160−0.850-
12−12.053.917−1.539−1.785-
139.180−0.935−2.1992.377-
142.500−6.815−2.2033.462-
15−3.1054.973−0.752−3.526-
16−4.0684.687−1.196−3.731-
178.0379.2910.5929.212-
186.7265.9790.9633.974-
194.042−5.0581.3765.431-
204.585−3.606−0.7776.701-
Table 8. The optimized coefficients for estimating the normalized k n in Equation (25).
Table 8. The optimized coefficients for estimating the normalized k n in Equation (25).
Neuron IndexInput Layer WeightsHidden Layer WeightsInput Layer BiasesOutput Layer Bias
i w 1 i , 1 w 1 i , 2 w 2 , i b 1 , i b 2
1−6.753−1.1031.041−9.530−0.106
28.4343.590−2.501−7.827-
3−5.5412.8700.1115.778-
4−2.502−4.938−0.160−5.132-
51.257−4.5510.671−3.641-
6−6.8860.1570.5232.501-
72.427−3.9041.129−1.205-
83.711−4.231−0.6800.426-
9−4.383−2.456−4.2282.781-
103.7812.6280.7813.734-
114.197−0.920−0.658−0.850-
12−5.9867.171−0.378−1.785-
135.4294.213−2.2852.377-
143.700−7.289−1.8253.462-
154.0373.7232.991−3.526-
165.4322.211−1.677−3.731-
177.672−3.6915.3469.212-
18−1.7576.1142.9713.974-
192.7192.9062.7675.431-
20−7.991−2.6371.7336.701-

Share and Cite

MDPI and ACS Style

Gowida, A.; Elkatatny, S.; Ramadan, E.; Abdulraheem, A. Data-Driven Framework to Predict the Rheological Properties of CaCl2 Brine-Based Drill-in Fluid Using Artificial Neural Network. Energies 2019, 12, 1880. https://doi.org/10.3390/en12101880

AMA Style

Gowida A, Elkatatny S, Ramadan E, Abdulraheem A. Data-Driven Framework to Predict the Rheological Properties of CaCl2 Brine-Based Drill-in Fluid Using Artificial Neural Network. Energies. 2019; 12(10):1880. https://doi.org/10.3390/en12101880

Chicago/Turabian Style

Gowida, Ahmed, Salaheldin Elkatatny, Emad Ramadan, and Abdulazeez Abdulraheem. 2019. "Data-Driven Framework to Predict the Rheological Properties of CaCl2 Brine-Based Drill-in Fluid Using Artificial Neural Network" Energies 12, no. 10: 1880. https://doi.org/10.3390/en12101880

APA Style

Gowida, A., Elkatatny, S., Ramadan, E., & Abdulraheem, A. (2019). Data-Driven Framework to Predict the Rheological Properties of CaCl2 Brine-Based Drill-in Fluid Using Artificial Neural Network. Energies, 12(10), 1880. https://doi.org/10.3390/en12101880

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