Research article

Prediction of gas-lift performance using neural network analysis

  • Received: 24 November 2020 Accepted: 07 February 2021 Published: 09 March 2021
  • Integrated gas lift system optimization plays an indispensable role in production and economics by maximizing the revenue from a gas lift field. This requires: optimization of gas lift parameters, finding the best tuning of the completion and surface production systems to keep pace with the dynamic reservoir changes along with saving gas quantities and compression costs. Accordingly, a comprehensive study is being carried out to measure the capability of the Artificial Neural Network (ANN) and Machine Learning (ML) in the optimization of gas lift parameters. The results of this study show the power of two different mechanisms of neural network (NN) which are Radial Base Function (RBF) and Back Propagation Function (BPF) to predict the most three important factors of the process: optimal gas injection rate, bottom hole pressure and flow rate and compare the findings with conventional methods. In addition, this work provides 3 functional equations that can be utilized by applying the field data with no artificial intelligence (AI) expertise or software knowledge. This effort provides forth an industrial insight into the role of data-driven computational models for the production recognition scheme, not only to validate the well tests, but also to reduce the uncertainties in production optimization. The work was completed by generating an economic analysis to illustrate the understanding of potential benefits of implementing irregular gas lift mechanisms in the field to stand on both technical and economic aspects of the study.

    Citation: Ahmed A. Elgibaly, Mohamed Ghareeb, Said Kamel, Mohamed El-Sayed El-Bassiouny. Prediction of gas-lift performance using neural network analysis[J]. AIMS Energy, 2021, 9(2): 355-378. doi: 10.3934/energy.2021019

    Related Papers:

    [1] Khaled M. Bataineh, Assem N. AL-Karasneh . Direct solar steam generation inside evacuated tube absorber. AIMS Energy, 2016, 4(6): 921-935. doi: 10.3934/energy.2016.6.921
    [2] Sunbong Lee, Shaku Tei, Kunio Yoshikawa . Properties of chicken manure pyrolysis bio-oil blended with diesel and its combustion characteristics in RCEM, Rapid Compression and Expansion Machine. AIMS Energy, 2014, 2(3): 210-218. doi: 10.3934/energy.2014.3.210
    [3] Saad S. Alrwashdeh . Investigation of the effect of the injection pressure on the direct-ignition diesel engine performance. AIMS Energy, 2022, 10(2): 340-355. doi: 10.3934/energy.2022018
    [4] Xiaoyu Zheng, Dewang Chen, Yusheng Wang, Liping Zhuang . Remaining useful life indirect prediction of lithium-ion batteries using CNN-BiGRU fusion model and TPE optimization. AIMS Energy, 2023, 11(5): 896-917. doi: 10.3934/energy.2023043
    [5] Gitsuzo B.S. Tagawa, François Morency, Héloïse Beaugendre . CFD investigation on the maximum lift coefficient degradation of rough airfoils. AIMS Energy, 2021, 9(2): 305-325. doi: 10.3934/energy.2021016
    [6] Charalambos Chasos, George Karagiorgis, Chris Christodoulou . Technical and Feasibility Analysis of Gasoline and Natural Gas Fuelled Vehicles. AIMS Energy, 2014, 2(1): 71-88. doi: 10.3934/energy.2014.1.71
    [7] Abdelrahman Kandil, Samir Khaled, Taher Elfakharany . Prediction of the equivalent circulation density using machine learning algorithms based on real-time data. AIMS Energy, 2023, 11(3): 425-453. doi: 10.3934/energy.2023023
    [8] Tien Duy Nguyen, Trung Tran Anh, Vinh Tran Quang, Huy Bui Nhat, Vinh Nguyen Duy . An experimental evaluation of engine performance and emissions characteristics of a modified direct injection diesel engine operated in RCCI mode. AIMS Energy, 2020, 8(6): 1069-1087. doi: 10.3934/energy.2020.6.1069
    [9] Anbazhagan Geetha, S. Usha, J. Santhakumar, Surender Reddy Salkuti . Forecasting of energy consumption rate and battery stress under real-world traffic conditions using ANN model with different learning algorithms. AIMS Energy, 2025, 13(1): 125-146. doi: 10.3934/energy.2025005
    [10] Kunzhan Li, Fengyong Li, Baonan Wang, Meijing Shan . False data injection attack sample generation using an adversarial attention-diffusion model in smart grids. AIMS Energy, 2024, 12(6): 1271-1293. doi: 10.3934/energy.2024058
  • Integrated gas lift system optimization plays an indispensable role in production and economics by maximizing the revenue from a gas lift field. This requires: optimization of gas lift parameters, finding the best tuning of the completion and surface production systems to keep pace with the dynamic reservoir changes along with saving gas quantities and compression costs. Accordingly, a comprehensive study is being carried out to measure the capability of the Artificial Neural Network (ANN) and Machine Learning (ML) in the optimization of gas lift parameters. The results of this study show the power of two different mechanisms of neural network (NN) which are Radial Base Function (RBF) and Back Propagation Function (BPF) to predict the most three important factors of the process: optimal gas injection rate, bottom hole pressure and flow rate and compare the findings with conventional methods. In addition, this work provides 3 functional equations that can be utilized by applying the field data with no artificial intelligence (AI) expertise or software knowledge. This effort provides forth an industrial insight into the role of data-driven computational models for the production recognition scheme, not only to validate the well tests, but also to reduce the uncertainties in production optimization. The work was completed by generating an economic analysis to illustrate the understanding of potential benefits of implementing irregular gas lift mechanisms in the field to stand on both technical and economic aspects of the study.



    Abbreviations: AAPRE: Average Absolute Percent Relative Error; FLP: Flowline pressure; APRE: Average Percent Relative Error; (𝑋p): Importance of input variable (𝑋p); Max. PRE: Maximum Percent Relative Error; Min. PRE: Minimum Percent Relative Error; MSE: Mean Squared Error; nh: Number of neurons in the hidden layer; np: Number of input variables; Oj: Output layer weight for jth hidden layer; WHP, T: Wellhead pressure, temperature; Pnode: Node Pressure (Psi); Poutlet: Outlet Pressure (Psi); Pwh: Wellhead Pressure (Psi); qg: Gas Lift Injection Rate (MMSCFD); R: Correlation Coefficient; R2: Correlation Coefficient Squared; RMSE: Root Mean Squared Error; RP: Average Reservoir Pressure (Psi); SD: Standard Deviation

    ANN's inception started by Warren McCulloch and Walter Pitts in 1943 [1] constructed the first computational model for Neural Network (NN) based on mathematical algorithms called threshold logic which was used for binary cases. This model paved the way for research to be divided into two systems. First system relied on biological processes meanwhile the other relied on the application of NN to AI. Several technical researches have discussed the applications of ANN at oil and gas industry, Elgibaly et al. (1998) [2-3] employed ANN in calculating the Optimal Hydrate Inhibition Policies Ghahfarokhi et al. (2018) used ANN to predict gas production, Khan et al. (2018) [4] utilized the ANN to predict the optimum production rate, luo et al. (2018) utilized the ANN to optimize the production in several fields, Nande (2018) utilized the ANNs to reach lowest error in predicting closure pressure for hydraulic fracturing analysis, Tariq (2018) [5] used ANNs to predict flowing bottom hole pressure. Bahaa (2018) [6] Utilized ANN to forecast the Fluid Rate and Bottom Hole Flowing Pressure for Gas-lifted Oil Wells. Radhi (2020) [7] assigned the Numerical Simulation to optimize gas lift quantity using Genetic Algorithm for a Middle East Oil Field.

    In the field of mathematical models there are several functions for neural networks and this paper is focused on only two types which are backpropagation function network and radial basis function network. Continuous and timely lift optimization on a well is dependent on knowledge of the operational conditions of the well and its associated reservoir parameters. As the target is to maximize the economic production values from a well at a given point in time, the operational expenditure can be inversely associated with production maximization. generating set-point recommendations to optimize lift performance at timely intervals considering the above-mentioned setup can be a significant challenge. Popularly, the gas lift injection rate set-point is decided by operators based on historical experience estimates. The main target of gas lift optimization is expanding the current production from the well by allocating the relation between the rates of gas injection and oil production. The main factor to describe this relationship is the marginal increment in oil production rate per unit change in the gas lift injection rate. With a low quantity of injected G/L can result in production rate loss due to insufficient reduction in production fluids gravitational head. On the contrary, overvalues of injected gas rate led to high frictional head and wellhead pressure resulting in additional backpressure on the downhole formation, this results in a significant loss in production. (Radhi 2020) [8-9]. The major parameters controlling the previously mentioned balance includes but are not limited to: reservoir formation pressure or SBHP, Pi or inflow performance, WHP, depth of injection point, produced gas oil ratio (GOR), Water Cut %, API gravity, tubing size, and tubing roughness (such as friction factor). Some of these factors are not continuously measurable, and in many situations, estimation is provided based on old measurements or proximal well properties. The variation in the well behavior further complicates the equation. As well, the variation in production values, GOR, and WC % in a short span of evaluation may be due to the individual well's daily production rates not being physically measured but mathematically allocated. The bearing in production values and associated factors might be generated from the well natural decline or other interventions such as a frac hit, workover or re-stimulation. With respect to all of this, an injection rate which is optimum at a given point in time may not be so at a later stage. This challenge leads to attempt to understand the underlying state and represent the well using physics-based models (Khan, 2020) [10]. In the present study, the data of total 205 wells from more than 5 different fields that are artificially subjected to gas lift were selected for training and testing for the proposed model, wells are located in the Western, Eastern Desert and offshore field in Egypt. For each of these wells, the set of data were gathered for subsequent use in the study to predict three important variables output by using conventional commercial software then they were compared with two different neural network models "Radial Base and Back Propagation" built by MATLAB 2018A to predict the optimal values of gas injection quantities, bottom hole flowing pressure and well productivity. The ultimate target of this approach is to:

    ● Provide a new way to speed up the calculations;

    ● Dispense the use of tedious and time-consuming software training;

    ● Show an attempt to cut the cost values of the engineering software supply.

    Table 1.  Comparison between the proposed different neural network (NN) techniques.
    Description Step 1: Feedforward and Backpropagation Step 2: Radial Basis Function Network (RBFN)
    Signal transmission Feedforward Feedforward
    Process of building the model The lag is only visible in RBF in the output layer. Two different, independent stages at the first stage by means of radial basis functions the probability distribution is established; the network learns the relations between input x and output y at the second stage.
    Threshold Yes No
    Type of parameters Weights and thresholds Location and width of basis function and weights binding basis functions with output.
    Functioning time Faster Slower (bigger memory required)
    Learning Time Slower Faster
    Description in Details Neural network is an algorithm coming from the neurons of the brain. It is important to recognize the patterns from complex data. This network consists of layers that contain neurons which are connected togetherIn the hidden layer, each neuron is connected to all neurons of the previous layer by the weights and has an activations number according to this equation:w1a1 + w2a2 +... + wnan = new neuronThis layer has an activation function as sigmoid function that is used. This function tries to scale the number of each neuron between 0 and 1 according to its relation:Sigmoid = σ = 11 + e − x.until getting the output.In proposed network, it utilizes two hidden layers, each has 50 neurons and two different activation functions which are sigmoid and linear functions. The weights are initially set randomly and the input data are normalized to be between 0 and 1. RBFN is a type of an ANN that is commonly used with function approximation challenges due to its universal approximation and faster learning & processing rate. Radial basis function networks have many uses including but not limited to system control, classification, time series prediction and function approximation.They were first formulated by Broomhead and Lowe (1988).For the purpose of the paper, RBFN can be trained for every well based on the well test data. Later it can be used on a daily basis or in real time for well rate prediction based on other measurements such as WHP, FLP, WC%, …. etc.ϕ(r) = exp (−r2/2σ2) where σ > 0.

     | Show Table
    DownLoad: CSV

    Adopted ANN in this work depends primarily on wells' actual test data obtained from test separator and measuring devices installed on both mobile test package, flow lines and gas lift lines. In addition, as ANN can deal with limited or faulty data while acting with data with uncertainty which is considered a proven advantage of ANN over any analytical conventional methods, downhole data obtained from static and flowing surveys using downhole memory gauges, production logging tools (PLT) using electrical line, reservoir rock and fluid properties obtained from PVT lab analysis and core lab analysis were used in this work to run different models to reach the optimum results with satisfactory accuracy. The input data set was randomly divided into 70% for training, 15% for validation, and 15% for the primary test. Training data are used to improve the network according to their error. Validation data are used to evaluate network generalization, and to stop training when generalization stops improving. Test data do not affect training, so they provide an independent measure of network performance during and after training. To effectively model the unique behavior in gas-lift wells, data such as the gas injection rate and depth of gas injection would have been required. None of the records of the gas-lift wells, in the compilation of more than 180 wells, provided these data. Owing to the insufficiency of comprehensive gas-lift well records for training and testing, a neural network that would predict the temperature profile in gas-lift wells was not developed. All data for gas and gas-lift wells were discarded. Furthermore, because of the presence of outliers and anomalies the database was reduced to 50 wells from different fields. The data ranges including minimum and maximum values of the input elements parameters assigned in different generated ANN models are listed in Table 1. Based on the oil API gravity. Twenty-two wells were removed because of extremely low oil gravities ranging from 10 to 19°API. In addition, to make the oil viscosities reported for the databases compatible, the viscosities in Database were neglected. because most of wells have very close API values except 3 or 4 wells, in addition, after testing the importance of each inputs, it is noticed that API is of less importance in the contribution.

    A spread of values for that parameter in the databases must exist, this permits the NN to more easily approximate the function. The parameter must not be dependent on other input variables only. A parameter may be dependent on other input variables but must also be dependent on some parameter that is not an input variable. In this method, the variable will supply more data about the well that is not already provided and known by the other variables. All the variables met these criteria except for water specific gravity, water flow rate, and dead oil viscosity. The water specific gravity was removed because the values reported for most of the wells were the same. Although for many of the wells the water flow rate was varying, the water flow rate was retained because it is an important parameter in describing the hydrodynamics of the system. Dead oil viscosity was dependent on oil °API gravity and temperature. However, since the dead oil viscosity for each well was determined at 100 °F, then the dead oil viscosity was dependent on a single variable, the oil °API gravity. To observe the effect of each parameter on the predictive power of the network, two networks were developed: one is forward and backward propagation and one radial base function network.

    In this paper, more than 2900 production real test data of 16 elements were gathered & investigated for inconvenience and processed, 16 valuable elements were chosen as inputs and 3 elements were selected as outputs representing Well flow rate (BPD), Gas lift rate (MMSCFD) and bottom hole flowing pressure (PSIA). Before the data sets were fed to the networks, the values of the input variables were normalized by dividing each value of the data by the maximum absolute value of data. Using below equation, every field in the data sets was normalized, reducing the range of the values from 0 to −13, 000 to between 0 and +1. Note that the limits are the same as those of the output of the hidden layer tan-sigmoid transfer function: x=xmin(x)max(x)min(x).

    In this study, the selected 2 different algorithms were conducted to compare the performance, the first one is backpropagation learning algorithm and the reason behind choosing it for this training is because it was desirable to adjust the weights and biases based on the error signal arising from the output layer. The goal of this algorithm was to reduce the mean-squared error of the network outputs for all the input sets to a global minimum value. At the end of every epoch, the mean-squared error was calculated. Momentum was also added to the backpropagation algorithm to reduce the possibility of the learning process getting stuck in a local minimum. To further increase the speed of convergence of the mean-squared error to a global minimum, the Levenberg-Marquardt method was used to update the weights and biases for forward and backward propagation and Radbas function for radial base function NN. Traditionally training is halted when the network error reaches the performance goal of the mean-squared error-specified or when the network has converged. To mitigate network over-training, the early stopping method was applied. This required monitoring the mean-squared error of both the training data set and validation data set.

    The accuracy of the proposed neural network models (Backward propagation and Radbas) was tested by different data from other fields by comparing the predicted following three parameters (fluid rate, gas injection and bottom hole flowing pressure profile) to actual tested data from mobile test package and generated from the commercial software.

    For finding the optimum network design, trial-and-error attempts were undertaken, starting with one hidden layer and going to two hidden layers with trials many numbers of hidden units was set almost equal to the number of inputs multiplied by two. Hidden units were then gradually added. The maximum number of hidden units is rarely requiring to exceed more than 4 times the number of inputs. The architecture was retrained at least 3 times (up to 10 times is recommended) with different initial weight randomizations and only the best one was saved for comparison with other architecture. Nodes optimum number needed in the hidden layer is a contingent issue, depending on the intricacy of the input and output mapping, the amount of noise in the data and the amount of training data available. If nodes number in hidden layer is less than the optimum, the back-propagation algorithm will fail to reach to a minimum during training. Conversely, too many nodes will result in the network over fitting the training data, resulting in poor generalization performance. The developed neural network contained 16 input variables, one hidden layer with 50 neurons and 3 output variables. Log-Sigmoid function (logsig) was used as a transfer function in the hidden layer.

    1- In the applied field, each well has its own data represented in WHP, FLP, FLT, Pr, WC%, GOR, injected volume, and injected rate…etc. and by using a total of 16 inputs to build a model using commercial software "PIPESIM" to predict the target data (bottom hole flowing pressure and flow rate with optimum gas injection) which is essential in reallocation and optimization of the limited gas among the 52 wells in the field for testing then to be applied on different fields for accuracy testing.

    2- Building two different ANN models "REDBAS and Feed-forward, back propagation" and compared with actual data and commercial software output data to show the power of computational program over the conventional software in solving problems that are hard to be modeled analytically.

    3- Development of usable empirical equations for calculating the three desired outputs, can be applied without the necessity of AI software or expertise.

    4- The economic analysis model was generated to calculate the commercial benefits of implementing gas lift optimization technique in the field and to be applied in several fields.

    Using normalized data, one hidden layer with 50 neurons and Log-Sigmoid activation function and the results as shown below in Figures 1 and 2.

    Figure 1.  MSE and error histogram for 1 hidden layer.
    Figure 2.  Regression results for one hidden layer.

    Using normalized data, two hidden layers with 50 neurons and Log-Sigmoid activation function and the results as shown below in Figures 4 and 5.

    Figure 3.  Network architecture (2 hidden layers).
    Figure 4.  MSE and error histogram for 2 hidden layers network.
    Figure 5.  Regression results for 2 hidden layers network.

    Using normalized data with spread parameters equal one to reach zero error, one hidden layer with 50 neurons and Radbas activation function, and the results are as shown below in Figures 7 and 8.

    Figure 6.  Network architecture (Radial base NN).
    Figure 7.  Training and validation of RADBAS NN.
    Figure 8.  Testing results and MSE for radbas NN.

    Three Approaches' error histogram for the three targets: well fluid rate, bottom hole flowing pressure and optimum gas injection rate respectively, shows error histogram almost close to actual data as shown in Figures 1, 4 and 9 which resulted in a very good estimation in the running data. Moreover, the regression values appear very close to one in training, validation and testing with values up to 0.995 which have an impact on efficient trained neural network.

    Figure 9.  Error histogram for 3 outputs (Radbas NN model).

    Based on the trials carried out to reach the optimum developed ANN model, 60 points were assigned from different fields as a validation and approbation of the model efficiency in predicting optimum Well rate (BFPD), optimum gas injection rate (MMSCF) and bottom hole flowing pressure (Psia), Tables 3 and 4. Figures 14 and 15 proved a significant matching between the predicted results and actual results for the three targets.

    Table 2.  The data range for the parameters used in the ANN model.
    Parameter Min. Max. Avg.
    Gas Upstream Pressure (Psia) 800 1150 975
    Gas downstream Pressure/injection (Psia) 40 1000 520
    Well Head Pressure (Psia) 40 550 295
    Flow Line Pressure (Psia) 40 220 130
    Flow line Temperature © 15 65 40
    Separator Pressure (Psia) 40 90 65
    WC (%) 1 99.9 50.45
    Bubble point pressure (Psia) 900 1200 1050
    Downhole Temp © 35 80 57.5
    Well Depth (ft) 3000 9500 6250
    Depth of Injection (ft) 1500 7800 4650
    Current Reservoir Pressure (Psia) 500 1800 1150
    Initial Reservoir Pressure (Psia) 1500 1999 1749.5
    Productivity Index J 0.052 20 10.026
    Average Permeability (MD) 30 1200 615
    Average Porosity (Ø) 18 30 24

     | Show Table
    DownLoad: CSV
    Table 3.  Comparison summary of performance parameters for different neural network models.

     | Show Table
    DownLoad: CSV
    Figure 10.  Pwf comparisons with different models.
    Figure 11.  Well fluid rate comparison with different models.
    Figure 12.  Optimum gas injection comparison with different models.
    Figure 13.  Workflow for weights produced between artificial neurons of the NN model [15].
    Figure 14.  Input data importance to the target.
    Figure 15.  Optimum production and without gas lift (natural flowing well) from each of the 60 wells.

    The extensive study carried out between the commercial software "Pipesim" and different ANN models as shown on below Figures 10, 11 and 12 illustrates the results which showed a clear result that ANN models are matched properly with actual data, on the contrary, Pipesim results showed miss alignment with the best fit curve "actual data curve" while applying Pipesim software, a different correlations has been applied ended by applying Duns and Ros (modified) due to It has a flow regime map extended by the work of Gould et al. [11]. This includes a new transition region between bubble and slug flow, and an additional froth flow region at high flow rates. The holdup is considered as no-slip for froth flow, and is interpolated over the bubble-slug transition, the other holdup relationships are as for the standard Duns and Ros [12]. Friction is calculated by the method proposed by Kleyweg. This uses a monophasic friction factor rather than two-phase, but involves use of an average fluid velocity. This is claimed by Kleyweg to be a better method. Duns and Ros Modified gives the highest pressure drops in the slug flow regime for oil wells.

    Applying a statistical analysis for the three targets; Well rate, bottom hole flowing pressure and gas lift rate results respectively based on Average Percent Relative Error (APRE), Average Absolute Percent Relative Error (AAPRE), Minimum Percent Relative Error (Min PRE), Maximum Percent Relative Error (Max PRE), Root Mean Squared Error (RMSE), Mean Squared Error (MSE), Correlation Coefficient Squared (R2), Correlation Coefficient (R) and Error Standard Deviation (SD).

    As mentioned in Table 2, ANN Output Forward and Backward Propagation Results with two hidden layer models produce the most accurate results, as they depict the lowest APRE, AAPRE, MSE, RMSE and SD, and Highest R2 and R, which means the error of this method is the closest to zero error, on contrary Pipesim models produce the least accurate results as they depict the highest APRE, AAPRE, MSE, RMSE and SD, and lowest R2 and R, which means the error of this method is the furthest to zero error.

    The proposed study provides the contribution of all inputs by using input and output weights of the hidden layers which was presented previously by Garson (1991) and repeated by Goh (1995) [13-14] in different aspects. This process essentially requires extracting input and output weights of hidden layers. And the reason behind using Garson is based on applying 7 different methods to assess the variable relative importance and Garson law showed much accurate and better results than the other methods and matches with the results of MLR in terms of the partial regression, In this paper, to calculate the relative importance of the sixteen inputs on the three target outputs, assessment process based on the weight matrix of the proposed optimized network and Garson's modified equation have been used. The equation below as follows:

    Ij=m=Nhm=1((|wihjm|k=Nik=1|Wihkm|)×|Whomn|)k=Nik=1{m=Nhm=1(|Wihkm|/k=Nik=1|Wihkm|)×|Whomn|} (1)

    where, Ij represents the relative importance of the jth input on the output, Ni refers to the number of input neurons similarly, Nh refers to the number of hidden neurons in hidden layers, and the terms Wih and Who are connection weights between the input and the first hidden layers and Who is the connection weights between the second hidden layer and the output layer, subscripts k, n and m refer to input, hidden and output neurons, respectively. Table below demonstrates the weights produced between artificial neurons of the neural network model used in this study.

    Table 6 lists the relative importance of various input parameters on both oil rate output and gas lift rate output respectively. As can be generally seen, water cut has the greatest impact on oil rate and gas lift rate prediction followed by net pay thickness (or producing interval), most of the parameters have almost equal importance in the average range of (4%–8%), separator pressure and flow line pressure have the least impact on the output parameters in the developed ANN model.

    Table 4.  "Input" sample points used to validate the model.

     | Show Table
    DownLoad: CSV
    Table 5.  Actual vs predicted sample results from ANN models.

     | Show Table
    DownLoad: CSV
    Table 6.  Models relative importance of input parameters on gas lift injection, well fluid rate and bottom hole flowing pressure outputs.
    Parameters GL injectionRelative Importance % Gross Well Fluid Rate Relative Importance % Pwf Relative Importance % AverageRelative Importance %
    WC % 15% 6% 8% 9%
    Gas Upstream pressure 8% 7% 6% 7%
    DH Temp. 7% 7% 8% 8%
    Depth of injection 7% 7% 7% 7%
    Bubble Point press. 6% 6% 6% 6%
    Well Depth. 6% 9% 9% 8%
    Sep. Press. 6% 4% 5% 5%
    Current Reser. Press. 6% 6% 7% 6%
    Injection Pressure 6% 5% 4% 5%
    WHP 6% 4% 4% 4%
    Permeability 5% 9% 7% 7%
    FLP 5% 4% 4% 4%
    Initial Reser. Press. 5% 5% 5% 5%
    Porosity 4% 7% 6% 6%
    PI 4% 9% 8% 7%
    FLT 4% 7% 6% 5%

     | Show Table
    DownLoad: CSV

    The results of Garson calculation demonstrate the following points:

    ● Water cut has the greatest proven impact on the targets: well rate and optimum gas injection rate and bottom hole flowing pressure. Results in water cut increases in a well, the total pressure gradient in the well will increase due to the increase in liquid density as water is heavier than oil, thus causes a decrease in well rate and necessitates increasing gas lift rate to bring oil production rate to surface with previous value of gas injection.

    ● Well depth and downhole temperature along with well permeability come in second place in relative importance. As the well depth increases, the required gas injection values increase, to lift the liquid upward in parallel with downhole temperature which has a shown impact on the well fluid mobility, furthermore, the effect of permeability appears with same importance values of well depth and downhole temperature around 7% indeed that when permeability values are ranged between 200 up to 800 (md) has a great impact on well fluid rate (9% relative importance).

    ● The rest of the parameters related to wellbore and reservoir fluid properties has almost same impact on the desired target; well rate, optimum gas injection and bottom hole flowing pressure values.

    The least impact resulted from separator pressure and flow line pressure which have minor influence on the targets, so it could be neglected in ANN models.

    In addition to the previous work achieved for prediction of the well flow rate through ANN along with MATLAB, an actual and usable empirical model could be applied with the elimination of AI programs and technical expertise. The generation of this equation is based on utilizing a group of biases and weights which relates the layers from the output to the input. The weights and biases contributing to their neurons; those linked with input-hidden layer is characterized by w1, whereas, those linked with the hidden-output layer is called w2, moreover, the hidden and output layer biases are b1 and b2 respectively. The novel empirical correlation developed by this AI technique to estimate any of three outputs rates in gas lift oil wells is given by the equation below:

    O3=W31+e(W21+e(W1P+b1)+b2)+b3 (2)

    where O2 refers to the output from the second layer which is a vector containing three main outputs (gross fluid, Produced Gas and PWF). The term W2 represents weights of second layer as matrix has 3 rows corresponding to the number of outputs and 50 columns corresponding to the number of neurons. Similarly, W1 is a matrix with 50 rows and 16 columns corresponding to number of neurons and inputs respectively, this term refers to the weights of the first layer. The term P is the input vector, b1 and b2 are the biases of first and second layer.

    First, the dataset is organized in a column vector that has 16 values which is denoted in the equation by term "P". Since there are two hidden layers, there are subscripts 1 and 2 corresponding to them. In the first hidden layer, the number of neurons is 50 and each neuron is connected to all input data by the weights. So, the weights were organized in a matrix with 50 rows corresponding to number of neurons and 16 columns corresponding to number of data set. And the bias is organized in a column vector that has 50 values corresponding to number of neurons which is an integer number used to optimize the output value. The value of each neuron in this layer is the summation of the product weights with the dataset then the bias is added. The output from this calculation is a vector that has 50 values for 50 neurons. These values are scaled by sigmoid activation function according to this equation (x)=11+ex. Similarly, the second hidden layer has also 50 neurons connected to the neurons of the first hidden layer by weights. The weights of this layer are organized in matrix with 50 rows and 50 columns. The values of these neurons are calculated with the same manner and using the same activation function. Finally, the last or third layer has 3 neurons and its weights are organized by a matrix that has 3 rows and 50 columns. The values of these neurons also are calculated as before but it is used purelin activation function according to this equation σ(x)=x.

    The new neuron in the L layer comes from the summation the outputs of multiplying the weights connect this neuron with previous layer neurons by the activations of the previous layer, as shown in the following relation:

    wL1aL11+w2LaL12++WLnaL1n= new neuron in L layer  (3)

    The accuracy of well flow rate formula is within range of variables data which are used to derive it. For this reason, many researchers tried to derive a unique production estimation formula for a special reservoir or a layer. In this paper, a global well rate formula is generated and validated by a set of unfiltered measured data points. Besides of inherent errors in measurement devices (gauge pressures, flow rate meters) and calculation procedures for liquid calculations, human error is one of the important errors in field operation. By examining the above-mentioned formula, it is revealed that −20% error in reported input data would cause an amount of +/−14% error in estimated flow rate. Also, just 9% error in choke size diameter would result more than 7% error in rate.

    An economical study was carried out to evaluate the commercial results of implementing gas lift optimization technique in the field. This analysis is performed by the equation adopted from Huh et al. (2010), Nakashima and Camponogara (2006) [16] as shown below. There are some assumptions that this analysis is based on, whereas oil and gas price are $55/STB, and $5, 500/mmscf respectively. The water disposal cost, gas lift operational cost, other operation cost is assumed to be respectively $1/bbl of water, $3, 500/mmscf and $8/STB produced well (Huh et al. (2010) [17].

    NetProfit=ni[qoi(pocop)+(pgqg)(qwicw)(Qgi,injcginj)] (4)

    N Wells Number; qoi Well Oil production, STB/day; qg Gas production quantity, MMscf/day; Qgi, inj Gas lift injection quantity, MMscf/day; qwi Water production quantity, STB/day; po Oil price, $/STB; po Gas price, $/MMscf; cop Operational cost per each STB of oil, $/STB.

    The calculated daily production of oil, water, gas, gas injection rate, and net profit for each well are presented in Figures 15 and 16 for natural flowing wells and gas lift wells respectively. The daily net profit with and without gas lift and overall gain in daily oil production rate for 60 wells are compared in Figure 17. As can be seen in the Figures, the gas lift technique can substantially increase the daily oil production rate of every well in the field; and thus, the net profit as compared to naturally producing wells.

    Figure 16.  Daily production and net profit with gas lift from each of the 60 wells.
    Figure 17.  Comparison of daily net profit with and without gas lift and oil production gain for 60 wells.
    Figure 18.  Comparison of the developed ANN model with other published work regarding (AAPE).

    Several works were carried out and discussed to compare the power of artificial neural network along with empirical correlations on gas lift parameters, below table illustrates the results reach throughout the studies.

    Table 7.  Comparison of the developed ANN model with other published work.
    Author, (date) No. of input No. of output AAPE MSE Regression
    Gilbert, (1954) 7 1 18.55 385 0.904
    Baxendell, (1958) 8 1 31.05 1554 0.911
    Ros, (1960) 8 1 26.27 15564 0.899
    Achong I, (1961) 7 1 37.08 1588 0.858
    Osman and Dokla, (1990) 6 1 85.38 645222 0.889
    Alrumah and Bizanti, (2007) 10 1 114.08 358865 0.895
    E.Khamehchi1, (2009) 9 2 1.3 0.005 0.9192
    Ranjan et al (2015) 10 2 5.5 44764650.82 0.89523
    Mazn, (2018) 7 2 9.6 1210 0.89
    Khan, (2019) 14 1 2.56 330 0.96
    The paper results 16 3 0.18 .0012 0.995

     | Show Table
    DownLoad: CSV

    The paper results showed great results compared to the other papers observed in AAPE, MSE and regression as well due to enormous actual dataset used, different normalization technique applied to reach the optimum input data range and different model for neural functional networks.

    This study primarily focuses on exploring the feasibility of the implementation of ANN-based optimization technique in numerical modeling for optimizing the gas lift wells on daily basis in a large field with complex network system. Accordingly, ANN technique is utilized to optimize the allocation of the continuous gas lift injection rate for 60 wells. The principles of ANN, and mathematical model including the workflow for performing the simulation studies are comprehensively discussed in this paper. Sensitivity studies and sample economic analysis were also performed to get an insight into the benefit of implementing gas lift techniques for depletion drive reservoir, especially in the event of increasing the water cut and very low reservoir pressure.

    ANN technique appears to be an efficient technique with an ability to model a large number of wells produced concurrently in a network system for the prediction of optimally allocating the gas injection rate towards maximization of oil production rate while maintaining the optimum bottom hole flowing pressure.

    Gas lift technique is found to be more beneficial for wells with relatively higher water cut.

    In the event of reservoir pressure depletion, gas lift appears to be not only beneficial for improving the well production performance but also for increasing the field life cycle by allowing the well to continue production even at a very low flowing wellbore pressure (Pwf) at its given minimum well head pressure.

    In this work, various training algorithms for BP and FP networks and RBF networks were nominated to test the prediction of three important parameters in well life cycle. Based on the study carried out, the following conclusions were achieved:

    1-First approach performed was a feed forward and backward propagation network, that was tested by using different numbers of neurons, one or two hidden layers and different training process which are steepest descent, Levenberg–Marquardt and Bayesian to get the best results. The second approach was a radial basis network. Both systems can exhibit advantages and disadvantages when compared to one another and in both cases the results were quite satisfying.

    2-RBFs can be trained much faster than perceptron's and achieved less mean square error with high number of neurons.

    3-After applying test data from the different fields found the fastest training and testing error were achieved with BP and FP neural network which showed in statistical analysis compared with Pipesim models with accuracy up to 91%.

    4-Stability is achieved by the RBFs network and the network trained by Levenberg–Marquardt algorithm contrary to the networks trained with the Bayesian and steepest descent algorithms.

    5-The experimental results indicate that a strong matching between model predictions and observed values, since MSE is 0.0012. When performance results are compared, it was concluded that RBFNN-based model is a more reliable predictor, with MSE value of 0.003 and ARPE of 8.2. Therefore, the smallest MSE value indicates a creditable method for accuracy, while RBF finding illustrates best proposed model to analyze the output.

    6-It has been shown that there is a minimum number of inputs required to achieve an accurate model and any further decrease in input number would result in increasing MSE and reducing regression value which affects the whole model performance.

    7-Using Garson algorithm, sensitivity analysis processed, the rate and the way of input data distribution with the highest impact on the model output is determined. With this process, the trial-and-error steps in the design process can be reduced and the most important effective parameters can be identified.

    8-A comprehensive comparison with antecedent studies showed that this work presents more reliable results due to using large actual data with respect to pre-processing for the data resulted in normalization with different methods, then applying two different function neural networks along with generating of a new empirical equation.

    9-A new empirical model is proposed for estimating the three outputs, it can be applied to any dataset given the input parameters are within the model's range. This would not require the expertise of coding or using convoluted software.

    10-Sample economic analysis demonstrates that the gas lift technique can substantially increase the daily oil production rate of every well in the field; and thus, the net profit as compared to naturally producing wells.

    The authors would like to thank the Department of Petroleum Engineering, Faculty of Petroleum and Mining Engineering, Suez University for useful discussion and cooperation on this research.

    The authors declare no conflict of interest.



    [1] Farshad FF, Garber JD, Lorde JN (2000) Predicting temperature profiles in producing oil wells using artificial neural networks. Eng Comput 17: 735-754. doi: 10.1108/02644400010340651
    [2] Tariq Z, Mahmoud M, Abdulraheem A (2019) Core log integration: a hybrid intelligent data-driven solution to improve elastic parameter prediction. Neural Comput Appl 31: 8561-8581. doi: 10.1007/s00521-019-04101-3
    [3] Mahmoud M (2019) New correlation for the gas deviation factor for high-temperature and high-pressure gas reservoirs using neural networks. Energy Fuels 33: 2426-2436. doi: 10.1021/acs.energyfuels.9b00171
    [4] Olabisi OT, Atubokiki AJ, Babawale O (2019) Artificial neural network for prediction of hydrate formation temperature. Soc Pet Eng. Available from: https://onepetro.org/SPENAIC/proceedings-abstract/19NAIC/2-19NAIC/D023S026R009/219498.
    [5] Tariq Z, Abdulraheem A, Khan MR, et al. (2018a) New inflow performance relationship for a horizontal well in a naturally fractured solution gas drive reservoirs using artificial intelligence technique. In: Offshore technology conference Asia. Offshore Technology Conference.
    [6] Khan MR, Tariq Z, Abdulraheem A (2020) Application of artificial intelligence to estimate oil flow rate in gas-lift wells. Nat Resour Res 29: 4017-4029. doi: 10.1007/s11053-020-09675-7
    [7] Zubarev M, Zubarev D (2019) Use of radial basis function networks for efficient well production allocation. SPE Nigeria Annual Int Conf Exhib. Available from: https://onepetro.org/SPENAIC/proceedings-abstract/19NAIC/2-19NAIC/D023S007R002/219944.
    [8] Ekundayo J, Rezaee R (2019) Effect of equation of states on high pressure volumetric measurements of methane-coal sorption isotherms-part 1: Volumes of free space and methane adsorption isotherms. Energy Fuels 33: 1029-1033. doi: 10.1021/acs.energyfuels.8b04016
    [9] Tariq Z, Abdulraheem A, Khan MR, et al. (2018a) New inflow performance relationship for a horizontal well in a naturally fractured solution gas drive reservoirs using artificial intelligence technique. In: Offshore technology conference Asia.
    [10] El-Moniem MAA, El-Banbi AH (2018) Development of an expert system for selection of multiphase flow correlations. J Pet Explor Prod Technol 8: 1473-1485. doi: 10.1007/s13202-018-0442-7
    [11] Baziar S, Shahripour HB, Tadayoni M, et al. (2018) Prediction of water saturation in a tight gas sandstone reservoir by using four intelligent methods: a comparative study. Neural Comput Appl 30: 1171-1185. doi: 10.1007/s00521-016-2729-2
    [12] Tariq Z (2018) An automated flowing bottom-hole pressure prediction for a vertical well having multiphase flow using computational intelligence techniques. J Pet Explor Prod Technol 354: 463-472.
    [13] Chen W, Di Q, Ye F, et al. (2017) Flowing bottom hole pressure prediction for gas wells based on support vector machine and random samples selection. Int J Hydrogen Energy 42: 18333-18342. doi: 10.1016/j.ijhydene.2017.04.134
    [14] Bazargan H, Adibifard M (2017) A stochastic well-test analysis on transient pressure data using iterative ensemble Kalman filter. Neural Comput Appl 31: 3227-3243. doi: 10.1007/s00521-017-3264-5
    [15] Jahed AD, Shoib R, Faizi K, et al. (2017) Developing a hybrid PSO ANN model for estimating the ultimate bearing capacity of rock-socketed piles. Neural Comput Appl 28: 391-405. doi: 10.1007/s00521-015-2072-z
    [16] Al-Fatlawi O, Hossain MM, Osborne J (2017) Determination of best possible correlation for gas compressibility factor to accurately predict the initial gas reserves in gas-hydrocarbon reservoirs. Int J Hydrogen Energy 42: 25492-25508. doi: 10.1016/j.ijhydene.2017.08.030
    [17] Shokir EM, Hamed MMB, Ibrahim ESB, et al. (2017) Gas lift optimization using artificial neural network and integrated production modeling. Energy Fuels 31: 9302-9307. doi: 10.1021/acs.energyfuels.7b01690
  • This article has been cited by:

    1. Okorie Ekwe Agwu, Saad Alatefi, Reda Abdel Azim, Ahmad Alkouh, Applications of artificial intelligence algorithms in artificial lift systems: A critical review, 2024, 97, 09555986, 102613, 10.1016/j.flowmeasinst.2024.102613
    2. Okorie Ekwe Agwu, Ahmad Alkouh, Saad Alatefi, Reda Abdel Azim, Razaq Ferhadi, Utilization of machine learning for the estimation of production rates in wells operated by electrical submersible pumps, 2024, 14, 2190-0558, 1205, 10.1007/s13202-024-01761-3
    3. Onyebuchi Ivan Nwanwe, Nkemakolam Chinedu Izuwa, Nnaemeka Princewill Ohia, Anthony Kerunwa, Nnaemeka Uwaezuoke, Determining optimal controls placed on injection/production wells during waterflooding in heterogeneous oil reservoirs using artificial neural network models and multi-objective genetic algorithm, 2024, 1420-0597, 10.1007/s10596-024-10300-2
    4. Daniel Chuquin-Vasco, Dennise Chicaiza-Sagal, Cristina Calderón-Tapia, Nelson Chuquin-Vasco, Juan Chuquin-Vasco, Lidia Castro-Cepeda, Forecasting mixture composition in the extractive distillation of n-hexane and ethyl acetate with n-methyl-2-pyrrolidone through ANN for a preliminary energy assessment, 2024, 12, 2333-8334, 439, 10.3934/energy.2024020
    5. Mariea Marcu, Approximation methods used to build the gas-lift performance curve based on reduced datasets: a comprehensive statistical comparison, 2024, 1091-6466, 1, 10.1080/10916466.2024.2404625
    6. Sungil Kim, Tea-Woo Kim, Suryeom Jo, Artificial intelligence in geoenergy: bridging petroleum engineering and future-oriented applications, 2025, 15, 2190-0558, 10.1007/s13202-025-01939-3
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5158) PDF downloads(523) Cited by(6)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog