Regular Article
Multiphase gasflow model of an electrical submersible pump
^{1}
Department of Computational Mechanics, UNICAMPUniversity of Campinas, São Paulo, Brazil
^{2}
StudentsAgreement Graduate  PECPG, CAPES/CNPq, Brazil
^{*} Corresponding author: dmartinez628@gmail.com
Received:
7
December
2016
Accepted:
2
July
2018
Various artificial lifting systems are used in the oil and gas industry. An example is the Electrical Submersible Pump (ESP). When the gas flow is high, ESPs usually fail prematurely because of a lack of information about the twophase flow during pumping operations. Here, we develop models to estimate the gas flow in a twophase mixture being pumped through an ESP. Using these models and experimental system response data, the pump operating point can be controlled. The models are based on nonparametric identification using a support vector machine learning algorithm. The learning machine’s hidden parameters are determined with a genetic algorithm. The results obtained with each model are validated and compared in terms of estimation error. The models are able to successfully identify the gas flow in the liquidgas mixture transported by an ESP.
© D. Martinez Ricardo et al., published by IFP Energies nouvelles, 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Different artificial lifting methods that ensure the flow of oil through a production line are used in oil and gas production systems to reduce production losses. One such method involves the use of Electrical Submersible Pumps (ESPs), which were developed in 1910 by Annals Arutunoff [1]. ESPs are centrifugal pumps that operate under multiphase flow conditions and are positioned at the end of the production piping. They are driven by a threephase electric motor and are suitable for use in wells with highviscosity fluids, high flow rates, high water content and a low gasoil ratio. To ensure smooth operation of the pump, the gasoil ratio should not exceed 10%. Exceeding this ratio can cause pumping to stop because of a phenomenon known as gaslock, a major limitation of ESPs [2]. To avoid this phenomenon and consequent pump failure, the percentage of gas in the ESP when it is operating under multiphase flow conditions (specifically a gasliquid mixture) must be known. One way to avoid the problem is by stopping the pump or changing the pump operating conditions. Various strategies for monitoring and analyzing the operation of ESPs have also been developed in an attempt to circumvent the gaslock problem. A traditional approach is to use ammeter charts, which can be used to identify operating problems. By comparing an ammeter chart showing the equipment operating normally and another recorded when an inspection is being carried out, the state of the equipment during the inspection can be determined, as this is reflected in changes in the current during the measurement period [1]. Another method for identifying a gaslock is to monitor the output flow from the pump so that the pump can be turned off when there is no fluid in it [3]. Monitoring of this kind, however, requires premature shutdown of the pump, potentially reducing well productivity. Other strategies that have been developed involve characterizing the flow by analyzing signals from measuring instruments installed in the pipeline [4] and studying the fluid frequency response by measuring the vibration of the pump [5, 6]. However, these methods can suffer from drawbacks depending on where and how the instrumentation is installed in the pipeline.
As the gaslock phenomenon in this context is not well understood and theoretical models are difficult to develop because of the complexity of the process, experimental studies are required to enable a model to be developed that allows the gas flow in these systems to be estimated. Such a model would also provide a better understanding of the real dynamics of the system of interest and be an important tool for estimating and analyzing the behavior of the system parameters and so preventing possible failures [7].
System identification is used to develop mathematical models based on a data set obtained from the system of interest. The methodology has been evolving for several decades [8]. In the 90s new areas of interest were considered, including system identification in the frequency domain, leading to nonparametric identification, which is currently widely used. The number of different methodologies being used is increasing steadily.
One of these new methodologies is based on machine learning, which involves various computational techniques that aim to increase the automation and efficiency of knowledge acquisition processes through the use of data processing. The main objective is to find the relationship between the system variables (input/output) using samples acquired from the system [9]. One of the many applications of machine learning is the identification of linear and nonlinear systems, a task at which Support Vector Machines (SVMs) are known to be particularly efficient compared with other learning methods [10]. Support Vector Machines are based on the structural risk minimization principle, which originated from statistical learning theory and was proposed by Vapnik. This principle has been shown to yield better results than the principle of empirical risk minimization used in neural networks [11].
In light of the increasing use of SVMs in different areas of research, the new methodologies that have been developed in an attempt to improve the operation and performance of SVMs and the good results obtained with these algorithms, an SVM was used here as the basis of a methodology for identifying the gas flow in an ESP.
2 Optimization of SVMR parameters with genetic algorithms
To develop a model that identifies the gas flow in an ESP operating with a twophase (liquidgas) fluid from data acquired directly from the system, a nonparametric identification method based on machine learning was developed. Specifically, a Support Vector Machine for Regression (SVMR) was used and the SVMR parameters were optimized with Genetic Algorithms (GAs).
Support Vector Machines are learning systems based on optimization tools that seek to minimize structural risk. They use a hypothesis space of linear functions in a highdimensional space created by a kernel function. These transformations can be used in various learning problems, usually either classification or regression problems [12].
Support Vector Machine for Regression aims to estimate a function f(x) where the output y _{i} will depend on an input x _{i}. Given a data approximation problem involving a set of training data x _{i}, y _{i} where x ∈ ℝ^{a} and x ∈ ℝ, the objective is to find a linear function f(x) that is an approximation to the system using a vector of minimum weights w,(1)where the training data sets are x = x _{1},x _{2},⋯, x _{ n } and y = y _{1},y _{2},⋯,y _{n}, the weight vector is w = w _{1},w _{2},⋯,w _{ n } and b is the bias.
To achieve this objective, the vector w must be minimized using the Euclidean norm  w ^{2}. The problem can then be defined as the optimization problem(2)subject to(3)
The system loss function is estimated by the parameter (ε), which introduces a degree of tolerance when samples are penalized. As the linear function f(x) may not be able to fit all the training data, slack variables (ξ _{i}) and are added to allow some errors subject to the following approximation conditions:(5)
Finally, the SVMR problem can be described by equation (6), where the goal is to find w and b that minimize w [13]:(6)where C > 0 is the parameter that penalizes permissible errors.
If the system is not linear, the training data belonging to space must be mapped to a higherdimension space using a function ϕ called a kernel:(7)
The problem discussed in this paper can be represented by a nonlinear system, and the kernel used to perform higherdimensional mapping is a Gaussian function. The parameter to be optimized with a GA is σ ^{2}, which determines the width of the function. The kernel function is given by:(8)
Figure 1 shows an example of transformation of the mapping space by the kernel function to a larger space where the nonlinear system behaves as if it were linear.
Fig. 1. Mapping the vectors into higherdimensional space. 
Equation (6) can be solved using dual programming [14]. The goal is to build an objective function by adding a set of variables α _{i} called Lagrange multipliers, where the dual formulation is given by(9)
After solving the dual problem the optimal decision is given by(10)
Correct selection of the hidden parameters of an SVM ensures that the models developed using the SVM perform satisfactorily. The parameters are C, which determines the curvature of the margin penalizing permissible errors, ϵ, which defines the insensitivity of the margin, and σ ^{2}, which determines the width of the kernel function. In this article, the use of a GA is proposed to search for the best hidden parameters of the SVMR.
Genetic Algorithms provide a robust method for performing searches and optimization. They mimic the evolutionary processes of living things and can provide solutions for a wide range of problems. Because of their adaptability, GAs have been used by many researchers to look for the hidden parameters of SVMRs for various problems [15, 16]. The flowchart of the GA used in this paper is shown in Figure 2.
Fig. 2. Flowchart of the genetic algorithm. 
To find the best parameters (C, ϵ,σ ^{2}), the algorithm must perform the following four steps. The first step consists of creating a population of individual elements or that are candidates for a possible solution. In the second step, the adaptation of each of the individual elements to the problem being studied is measured using a fitness function. In this particular case, the validation error of each of the models built with the SVMR is used as the fitness criterion, and the objective function is given by(11)where n _{val} is the number of validation samples.
In the third step, the individuals that are better adapted to the problem are selected. This is done by the roulette method, which not only ensures a higher probability of betteradapted individuals being selected, but also also allows individuals that are less well adapted to the problem to be selected, helping to maintain the diversity of the population of solutions and so avoid convergence to local minima.
In the fourth step, a crossover is performed between two or more individuals selected in the previous stage (parents) to generate descendants with part of the genetic material from each parent. For the crossover to be performed, all individuals must be coded in binary representation using the same number of bits for each parameter [e.g., P = (C, ∊, σ^{2}) = (2,4,3) → P_{bit} = (0010 0100 0011)]. Hence, if a parameter is represented by b bits, the genome size is l = 3b. In the crossover process, a random point k is selected and the bits to the right of this are swapped between the parents to create descendants [e.g., if the parents are P _{1} = (0010 0100 0011) and P _{2} = (0101 0001 1111), then for k = 5 we have S _{1} = (0010 0001 1111) and S _{2} = (0101 0100 0011) as descendants]. This crossover process is repeated until the population of descendants is the same size as the initial population. The algorithm includes a final step called the mutation operator that is intended to maintain population diversity. This operator randomly selects some of the individuals generated and performs a mutation in their genomes by changing one bit at random [e.g., a selected individual S = (0010 0001 1111) becomes S _{mutate} = (0010 0001 1011) after the mutation]. The above steps are repeated until the stop condition is satisfied.
3 Results
An SVMR was used to develop a nonparametric blackbox model for identifying the gas flow in a J200N ESP based on different characteristic system parameters. Casing vibration, total flow, torque and elevation were used as inputs, and a graphical representation of estimated gas flow in the ESP was used as output, since these are the data that are available in the ESP circuit tests. Several experiments were performed in the ESP test circuit in the LABPETRO laboratory at Unicamp using the initial operating conditions in Table 1. The ESP system was tested at 1800 rpm and 3000 rpm with manometric inlet pressures of 100 bar and 200 bar, respectively, and gas flows of 0–3 kg h^{−1}.
ESP testcircuit configuration for data acquisition.
After the data used to estimate the model parameters had been acquired, the signals were postprocessed. This is necessary because the acquired ESP casing vibration signals contain noise and are difficult to predict in the time domain even though the initial testcircuit conditions are known. A spectral representation of the signals was therefore used so that they could be characterized by their power spectral density. This was done with a finite Fourier transform and the autocorrelation function, a technique typically used to analyze measurement signals when working with digital equipment and discrete algorithms, as in this case [17].
After the spectral representation had been generated, the autocorrelation of the power spectral density was calculated for specific positive frequencies corresponding to the pump speeds and their respective multiples. For example, for a pump speed of 1800 rpm, the frequencies evaluated were 30, 60, 90 and 120 Hz, as these are the frequencies at which there are significant changes in the spectral representation when the gas flow in the ESP changes. In this way, a new set of data is acquired that can be used as inputs in the system identification. The patterns of the spectral representations of the signals associated with the casing vibration variable can be seen in Figure 3.
Fig. 3. Patterns of the spectral representation of casing vibration. 
Ten different models were built using a variety of signals as input, and once the parameters had been estimated the models were tested. The different sets of training samples used for each model are shown in Table 2. The objective of these 10 sets is to be able to select as input the best combination of the data acquired from the test system to be able to identify the gas flow in the EPS and thus obtain models that represent the variations of the percentage of gas in the pump.
Input and output data used with each model.
Having defined and acquired the model inputs and outputs, we used a genetic algorithm to estimate the best parameters (C, ϵ, σ ^{2}) and so identify the most representative model for each set of training samples. Various tests were performed for each of the ten models in Table 2. The configuration of the GA for each of the tests is shown in Table 3.
Configuration of the genetic algorithm in each test.
The results of the tests with each of the models using the GA configurations in Table 3 are shown in Table 4. Selection of the parameters to be used with the SVMR was based on analysis of the Mean Square Error (MSE) and autocorrelation coefficient. The parameters that yielded the lowest MSE and highest autocorrelation coefficient were selected and then used in the model to predict the ESP gas flow. Table 4 shows the configurations of the GA that yielded the best results and the corresponding SVMR parameters, MSE and autocorrelation coefficient for each model.
Selected SVMR parameters.
Once the input and output data for each model have been acquired and the SVMR parameters have been estimated by the GA, the gas flow in the ESP can be estimated using different characteristic system parameters. The power of generalization of each model was determined by measuring the performance of the model with a data set not included in the training data used to build the model. This test data set corresponded to 20% of the total data used and the rest of the data is used for training [18].
To evaluate the ten SVMR models, the actual gas flow was compared with the gas flow predicted by each model. For each of the four flows, a mean was estimated for all the test samples and the MSE was calculated for each model so that their performance could be compared. Table 5 shows the gas flows estimated by the SVMR models and the corresponding MSE.
Gas flow estimated by each SVMR model.
To evaluate the performance of the models, the gas flows estimated by each model were compared. Figure 4 shows the flows estimated by each of the models for actual flows of 0–3 kg h^{−1}. These results show that all the models have a good power of generalization but that the model with the best power of generalization is model 5, which is based on total flow and elevation. The gas flow estimated with this model has the lowest error of all the models.
Fig. 4. Gas flow estimates for each model. 
The model chosen to predict the gas flow in the ESP was the model with the lowest MSE and highest autocorrelation coefficient.
The results in Table 6 show that model 5, which uses total flow and elevation as input parameters, yielded the best results (MSE = 0.003 and autocorrelation coefficient = 0.997) and that the experimental data provided good information about changes in the gas flow in the pump, as the errors in the predicted values are small for most of the models. The findings also show that when the training data is of poor quality, the models do not produce good estimates.
MSE and autocorrelation coefficient of the SVMR models.
4 Conclusion
We have proposed different nonparametric VSMR models for estimating the gas flow in an ESP. The models use experimental vibration, totalflow, torque and elevation data collected from the pump and pipe system.
The hidden parameters of the SVMR were estimated using a GA, thus ensuring that the models have good power of generalization and, in turn, perform well. This was confirmed by testing the models with data that were not used during training. The models can be used to solve the problem of premature failure due to high gas flow in ESP systems. This problem arises because of insufficient information about the multiphase nature of the flow in the pump and can be easily avoided just by changing the pump operating conditions.
One of the most significant findings of this study is that total flow and elevation yielded the most accurate estimates of gas flow, as the models using these parameters had the lowest MSE (0.003) and highest autocorrelation coefficient (0.997).
References
 Tákacs G. (2009) Electrical submersible pumps manual: design, operations, and maintenance. Gulf Professional Publishing. [Google Scholar]
 Varon M.P. (2013) Estudo de uma bomba centrífuga submersa (BCS) como medidor de vazão, Msc Thesis, FEMUNICAMP, Campinas, Brazil. [Google Scholar]
 Snyder J.R., Dale R., Joe Haws H. (1989) Pump off/gas lock motor controller for electrical submersible pumps, EP Patent App. EP19,880,202,383, 03 May. [Google Scholar]
 Klein F.L., Seleghim, Jr. P., Eric H. (2004) Timefrequency analysis of intermittent twophase flows in horizontal piping, J. Braz. Soc. Mech. Sci. Eng. 26, 2, 174–179. [CrossRef] [Google Scholar]
 Kolpak M.K., Rock T.J. (1996) Measuring vibration of a fluid stream to determine gas fraction, U.S. Patent No. 5,524,475, 11 June. [Google Scholar]
 Henry M.P., Richard R.P. (2014) Multiphase flow metering system, U.S. Patent Application No 14/135,085, 19 Dec. [Google Scholar]
 Nelles O. (2001) Nonlinear system identification, SpringerVerlag, Berlin Heidelberg. [CrossRef] [Google Scholar]
 Ljung L. (1999) System identification: Theory for the user, Prentice Hall PTR, Upper Saddle River, NJ. [Google Scholar]
 Dietterich T.G. (1986) Learning at the knowledge level, Mach. Lear. 1, 287–315. [Google Scholar]
 Rychetsky M. (2001) Algorithms and architectures for machine learning based on regularized neural networks and support vector approaches, Shaker Verlag GmbH, Germany. [Google Scholar]
 Vapnik V. (1995) The nature of statistical learning theory, SpringerVerlag, New York. [CrossRef] [Google Scholar]
 Vapnik V., Golowich S., Steven E., Smola A. (1996) Support vector method for function approximation, regression estimation, and signal processing, Advances in Neural Information Processing Systems 9, pp. 281–287. [Google Scholar]
 Stoean R., Dumitrescu D., Preuss M., Stoean C. (2006) Evolutionary support vector regression machines, in: SYNASC '06 Proceedings of the Eighth International Symposium on Symbolic and Numeric Algorithms for Scientific Computing, IEEE Computer Society, Washington, DC, pp. 330–335. [Google Scholar]
 Zhang H., Yongmei L. (2013) BSPBased Support Vector Regression Machine Parallel Framework, Computer and Information Science (ICIS), 2013 IEEE/ACIS 12th International Conference on. IEEE, Niigata, Japan, pp. 329–334. [Google Scholar]
 Yu Q., Liu Y., Rao F. (2009) Parameter selection of support vector regression machine based on differential evolution algorithm. Fuzzy systems and knowledge discovery, de FSKD ’09. Sixth International Conference, pp. 596–598. [Google Scholar]
 Zhenyue H., Mei C. (2009) Soft sensor modeling using SVR based on genetic algorithm and akaike information criterion, intelligent humanmachine systems and cybernetics, Intelligent HumanMachine Systems and Cybernetics, 2009. IHMSC’09. International Conference on, Vol. 2, IEEE. [Google Scholar]
 Blackman R.B., Tukey J.W. (1959) The measurement of power spectra, Dover publications, New York. [Google Scholar]
 Michael K. (1996) A bound on the error of cross validation using the approximation and estimation rates, with consequences for the trainingtest split, Advances in Neural Information Processing Systems, pp. 183–189. [Google Scholar]
All Tables
All Figures
Fig. 1. Mapping the vectors into higherdimensional space. 

In the text 
Fig. 2. Flowchart of the genetic algorithm. 

In the text 
Fig. 3. Patterns of the spectral representation of casing vibration. 

In the text 
Fig. 4. Gas flow estimates for each model. 

In the text 