Next Article in Journal
A Framework of Runge–Kutta, Discontinuous Galerkin, Level Set and Direct Ghost Fluid Methods for the Multi-Dimensional Simulation of Underwater Explosions
Next Article in Special Issue
Air Entrainment in Drop Shafts: A Novel Approach Based on Machine Learning Algorithms and Hybrid Models
Previous Article in Journal
Development and Thermophysical Profile of Cetyl Alcohol-in-Water Nanoemulsions for Thermal Management
Previous Article in Special Issue
Development of a Methodology and Software Package for Predicting the Formation of Organic Deposits Based on the Results of Laboratory Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning Augmented Two-Fluid Model for Segregated Flow

1
BP America Inc., Denver, CO 80202, USA
2
Petroleum Engineering Department, Colorado School of Mines, Golden, CO 80401, USA
*
Author to whom correspondence should be addressed.
Fluids 2022, 7(1), 12; https://doi.org/10.3390/fluids7010012
Submission received: 7 December 2021 / Revised: 24 December 2021 / Accepted: 27 December 2021 / Published: 29 December 2021
(This article belongs to the Special Issue Multiphase Flow in Pipes with and without Porous Media)

Abstract

:
Segregated flow, including stratified and annular flows, is commonly encountered in several practical applications such as chemical, nuclear, refrigeration, and oil and gas industries. Accurate prediction of liquid holdup and the pressure gradient is of great importance in terms of system design and optimization. The current most widely accepted model for segregated flow is a physics-based two-fluid model that treats gas and liquid phases separately by incorporating mass and momentum conservation equations. It requires empirically derived closure relationships that have the limitation of being applicable only under a narrow range of input parameters under which they were developed. In this paper, we proposed a more generalized machine learning augmented two-fluid model, using a database that spans the range of various flowing conditions and fluid properties. Machine learning algorithms such as random forest, neural networks, and gradient boosting were tested for the best performing data-driven predictive model. The new model proposed in this work successfully captures the complex, dynamic, and non-linear relationships between the friction factor and flowing conditions. A comprehensive model evaluation against nineteen existing correlations shows the best results from the proposed model.

1. Introduction

Segregated flow, which includes stratified smooth, stratified wavy, and annular flows, is commonly encountered in several practical applications such as chemical, nuclear, refrigeration, and oil and gas industries. Accurate prediction of liquid holdup and the pressure gradient is of great importance in terms of system design and optimization. This section first reviews the two-fluid hydraulic model that is widely used in the industry for segregated flow pressure gradient and liquid holdup predictions, followed by a review of the common application of machine learning (ML) in multiphase flow modeling.

1.1. Review of Two-Fluid Hydraulic Models for Segregated Flow

The current most widely accepted two-fluid hydraulic models for steady-state gas-liquid two-phase segregated flow are primarily based on the mass and momentum conservation equations for gas and liquid phases separately [1,2], which is also referred to as a two-fluid model in the current study. Separate governing equations are derived by treating each phase or component as a separate entity. Figure 1 depicts a two-phase flow system, as well as the forces acting on the gas and liquid phase. Equations (1) and (2) provide the integral form of conservation of momentum equations for gas and liquid under steady-state conditions.
A L ( d p d L ) L τ W L S L + τ I S I A L ρ L g sin θ = 0 ,
A C ( d p d L ) C τ W C S C τ I S I A C ρ C g sin θ = 0 ,
where AL is the area occupied by the liquid phase, AC is the area occupied by the gas core, τWL is the liquid wall shear stress, τWC is the wall shear stress for the gas core, τI is the interfacial shear stress, SL is the liquid wetted perimeter, SC is the gas core wetted perimeter, and SI is the interfacial perimeter. dp/dL is the pressure gradient, and ρL and ρC are liquid and gas core densities respectively. The subscript L represents the liquid phase, C refers to the gas core with the entrained liquid droplet. θ is the inclination angle from the horizontal. vL and vC are the actual in-situ velocities for the liquid and gas core, respectively.
Equation (3) is the combined momentum equation of Equations (1) and (2) by canceling out the pressure gradient term, which can be used to solve for the liquid holdup, while Equations (4)–(6) provide the formulae to calculate shear stresses [1,3].
τ W L S L A L τ W C S C A C τ I ( S I A L + S I A C ) + ( ρ L ρ C ) g sin θ = 0 ,
τ W L = f L ρ L v L | v L | 2 ,
τ W C = f C ρ C v C | v C | 2 ,
τ I = f I ρ C ( v C v L ) | v C v L | 2 ,
To solve Equation (3), several closure relationships are needed, including liquid and gas wall friction factors (fL and fC), interfacial friction factor (fI), geometrical parameters such as SL, SC, SI, AC, and AL, and liquid entrainment (FE) defined by the entrained liquid mass flow rate in the gas core over the total liquid mass flow rate.
There are several different models for the geometrical parameters, such as the flat interface assumption that is widely used by various studies thereafter [3,4,5,6,7], apparent rough surface that assumes a liquid film with a constant thickness [8,9], and unified “double-circle” model and its simplified version applied to all inclination angles [10,11]. The geometrical parameters can be easily extended to annular flow by assuming a uniform film thickness [12,13].
Previous studies have shown that the interfacial friction factor plays the most important role in model performance. There are dozens of different interfacial friction factor correlations available in the literature, which were developed for specific experimental or field conditions. A previous study discussed the importance of the interfacial friction factor in gas–liquid stratified flow modeling and also emphasized the basic pitfalls in the commonly used single-phase closure relationships [1]. The primary reason behind the inaccuracy was attributed to the empirical nature through which the correlations were developed and the narrow range considered in the correlation development. Using the correlations outside the conditions originally used to develop the model can result in dramatic errors. Another problem is that the model prediction is discontinuous when switching from one correlation to another, which may potentially cause problems in production design. Due to the absence of a correlation built for a wide range of flow conditions, there is a desire for the industry to look for generalized correlations or approaches that can capture all the parameters that affect the pressure gradient and liquid holdup prediction, such as gas and liquid flow rates, inclination angle, fluid properties, pipe diameter, etc.
The wall friction factors in Equations (4) and (5) are commonly calculated using correlations for single-phase flow as functions of the Reynolds number, given in Equations (7) and (8), using hydraulic diameters. The formulae for the liquid and gas hydraulic diameters are given in Equations (9) and (10) by assuming open-channel flow for the liquid phase and closed-duct flow for the gas phase, respectively [3,14].
R e L = ρ L v L d L μ L ,
R e C = ρ C v C d C μ C ,
d L = 4 A L S L ,
d C = 4 A C S I + S C ,
Experimental data show that the use of the hydraulic diameter and single-phase correlation for the wall friction factor may also result in inaccuracy [15]. This conclusion can be obtained by comparing the experimental measured frictional pressure gradient and the one calculated from correlations, while assuming negligible errors from geometrical and entrainment correlations. The experimental measured frictional pressure gradient can be calculated from Equation (12), which is derived from Equations (1) and (2) by canceling out the interfacial shear stress terms, shown in Equation (11). The frictional pressure gradient from correlations can be obtained from Equation (13), in which the liquid and gas core wall shear stresses are calculated from Equations (4), (5) and (7)–(10). In this study, we used the liquid holdup obtained from experiments to calculate the geometrical parameters using a unified geometrical prediction model developed for all inclination angles [11], with a zero liquid entrainment assumption. The comparison of the two calculated frictional pressure gradients is plotted in Figure 2a for various datasets from [16,17,18,19] that cover all inclination angles from 0° to 90° and wide pressure and liquid flow rate ranges, showing unsatisfactory results. The sensitivity study on the single-phase wall friction factor correlation from [20,21,22,23,24] also shows that it does not have a significant impact on the conclusion (Figure 2b).
d p d L ) t o t a l = d p d L ) f + d p d L ) g = τ W L S L A p + τ W C S C A p f r i c t i o n a l   t e r m + [ ρ L H L + ρ G ( 1 H L ) ] g sin θ g r a v i t a t i o n a l   t e r m ,
d p d L ) f = d p d L ) t o t a l [ ρ L H L + ρ G ( 1 H L ) ] g sin θ g r a v i t a t i o n a l   t e r m ,
d p d L ) f = τ W L S L A p + τ W C S C A p ,

1.2. Application of Machine Learning in Multiphase Flow Modeling

As discussed in the previous section, closure relationships are used to estimate wall shear stresses and interfacial shear stress in a two-fluid model by empirically derived friction factor correlations [25]. Modeling of interfacial shear stress is one of the most critical issues in gas–liquid stratified flow models [1]. This is primarily due to the drag forces acting between the gas and liquid phases, which cause a considerable increase in the pressure gradient. Although mechanistic models are based on the fundamental laws of natural sciences, the closure relationships used in the process have a limitation of being applicable only under the conditions used for model development for acceptable accuracy [26,27].
In the past few years, multiphase flow solutions based on artificial intelligence were also developed and have proven to be a promising solution to model flow behavior. The advantage of using data-driven models lies in the reduced computational cost and a more simplistic approach without any user-induced bias in the form of model assumptions. A primary focus in the application of machine learning solutions in the past has been the prediction of flow patterns using pressure drop signal data as discussed by Mi et al. (2001), Shaban and Tavoularis (2014) [28,29]. Multiple other studies are based on fluid property experimental data using support vector machines, neural networks, or deep learning as a machine learning tool. Some of the related work was conducted by Osman (2004), Li et al. (2014), Ezzatabadipour et al. (2017), Rondon-Guillen et al. (2018), Mask et al. (2019), Hernandez et al. (2019), and Seong et al. (2020) [30,31,32,33,34,35,36]. The studies prove machine learning is an effective tool for predicting flow patterns based on pipe characteristics, superficial velocities, and other fluid properties. Mohammadi (2020) proposed an unsupervised learning model to determine clusters for closure relationships and used a genetic algorithm to find a near-optimal set of closure relationships and validate the results using statistics [37]. Lin (2020) incorporated physics by introducing a partial differential equation (PDE) aware deep learning model but was focused on solving the physical phenomenon of the saturation front [38].
Most of the previous studies are predominantly data-driven and did not capture the physics of the multiphase flow behavior, which limits their application range to the condition where the models were trained. In order to address the limitations of previous studies, this paper focuses on developing a machine learning model for friction factors and combining them with a physics-based two-fluid model to predict pressure drop and liquid holdup for segregated flow for a wide range of flow conditions and fluid properties.

2. Methods

In the following section, a hybrid physics-data-driven algorithm is proposed to model segregated flow using machine learning. Specifically, the new approach adopts the framework of a two-fluid model, incorporating the physics of gas and liquid segregated flow, and develops new generalized data-driven correlations for the liquid wall friction factor and interfacial friction factor using machine learning.

2.1. Model Development

Considering the inaccuracy using wall friction factor correlations developed for single-phase pipe flow as discussed previously (Figure 2), we proposed a new liquid wall friction factor (fL) data-driven model for two-phase segregated flow, which is a function of the liquid wall friction factor for single-phase pipe flow (fL-SP), given in Equation (14). The reason for choosing the liquid wall friction factor instead of the gas wall friction factor is the consideration of annular flow in which gas has zero contact with the pipe wall. The new data-driven model for fL is applicable to both stratified and annular flows.
f L = ϕ f L S P ,
With Equations (4)–(6) and (14), the gas and liquid conservation of momentum equations at steady-state becomes:
A L ( d p d L ) L ϕ f L S P ρ L v L | v L | 2 S L + f I ρ C ( v C v L ) | v C v L | 2 S I A L ρ L g sin θ = 0
A G ( d p d L ) G f G S P ρ C v C | v C | 2 S C f I ρ C ( v C v L ) | v C v L | 2 S I A C ρ C g sin θ = 0
where fL-SP is the liquid wall friction factor that can be calculated from correlations developed for single-phase pipe flow such as [20,21,22,23,24], using hydraulic diameter given in Equation (9); fG-SP is the gas wall friction calculated from correlations [20,21,22,23,24], using hydraulic diameter given in Equation (10).
In addition to ϕ, we also proposed a new data-driven machine learning model for fI, that is supposed to capture the complex behaviors under various flowing conditions and fluid properties. ϕ and fI can be directly calculated from experimental data for training purposes if both the pressure gradient and liquid holdup are known. The formulas to determine ϕ and fI from experimental data sets are given in Equations (17) and (18), which are derived from Equations (15) and (16).
ϕ = 2 A p ( d p d L ) L f G S P ρ C v C | v C | 2 S C [ ρ L H L + ρ G ( 1 H L ) ] A p g sin θ f L S P ρ L v L | v L | S L ,
f I = 2 A p H L ( d p d L ) L + ϕ f L S P ρ L v L | v L | 2 S L A p + ρ L H L g sin θ ρ C ( v C v L ) | v C v L | S I ,
In the current study, a data-driven methodology is proposed to model ϕ and fI to capture their non-linear relationships with flowing conditions and fluid properties for determining pressure drops and liquid holdups. The data source used in the model development covers a wide range of flowing conditions and fluid properties, discussed in the next section. The back-calculated interfacial friction factor, fI, and the coefficient, ϕ, from experimental measurement, become the target/response variable, and machine learning can be used to generate a correlation for accurately predicting this response variable while using fluid properties and flowing parameters as input variables. Figure 3 shows the stepwise workflow of the model development using machine-learning (Steps 1–5) and the implementation of the new model to predict pressure gradient and liquid holdup for segregated flow (Step 6).
The biggest advantage of the proposed method is that it preserves the physics of two-phase segregated flow while generating more robust and generalized interfacial and liquid wall friction factor correlations.

2.2. Dataset Description

The dataset was extracted from the open literature, including multiple studies where pressure gradient and liquid holdup were both presented, such as [16,18,19,39,40,41,42,43,44,45,46,47,48]. The two response variables, the coefficient for liquid wall friction factor, ϕ, and the interfacial friction factor, fI, were estimated from the experimental data first. These two parameters were correlated with the input variables using three machine learning algorithms—random forest, gradient boosting, and an artificial neural network. The input parameters include liquid and gas superficial velocities (vSL and vSg), liquid and gas densities and viscosities (ρL, ρG, μL, and μG), surface tension (σ), pipe diameter (d), and inclination angle (θ). The total number of filtered data points is approximately 1500.
Exploratory Data Analysis (EDA) is a critical step for understanding quantitative variables in an experimental dataset and is often used as a visual tool to understand the high-dimensional dataset. This step plays a critical role in learning from the data during the scientific process of model building and testing and discovering patterns in the dataset. In this section, the dataset is analyzed in a univariate and multivariate graphical sense by looking at statistical distributions in the form of histograms, scatter plots, and box plots. Data pre-processing includes the process of data wrangling by identifying the irregularities, cleaning the data, identifying, and removing anomalies and outliers as well as data transformation and feature selection.
There are four basic questions about the data that need to be asked before starting with data analysis. The first question relates to whether the data are discrete or continuous. The second step involves looking at the symmetry using distributions to identify if there is skewness present in the dataset. The third step is related to the upper and lower bounds of the data, and the final question determines the likelihood of observing extreme values in the distribution. The data used in the study are primarily of continuous form with values in a finite interval. Figure 4 shows the univariate distribution of all the variables that form the dataset. This figure also helps to determine the wide range of each variable, which is essential to create a model with high generalization capability. As an example, it can be observed from the inclination distribution plot that a majority of studies have been conducted on horizontal pipes.
A multivariable bivariate distribution leads to an identification of significant skewness in liquid density, gas density, liquid viscosity, and gas viscosity, but it is not necessarily an outlier, rather another subset of data points with lower frequency.
Box plots are often used as a graphical tool to display information for continuous univariate data and to identify highly influential points or possible outliers. The classical method from Tukey (1977) [49] was used to remove outliers by using a threshold of 1.5×interquartile range (IQR), which is defined as the difference between the first and third quartile [50]. Predictions for ϕ and fI using the two-fluid model were analyzed to identify and remove outliers. The results for the case of ϕ and fI before and after the outlier removal are shown in Figure 5.

2.3. Machine Learning Algorithms

Three different machine learning algorithms are implemented for the dataset—Random Forest [51], eXtreme Gradient Boosting—XGBoost [52], and Artificial Neural Networks [53]. The objective is to generate a predictive model for ϕ and fI with the highest accuracy. Although these three algorithms could be used for both regression and classification problems, in this work the focus will be supervised regression since the dataset has continuous variables, with no categorical information. A primary reason for the selection of these algorithms is their high performance to account for the non-linearity associated with the predictor variables or features [54]. A detailed explanation of each algorithm is readily available in the literature and is omitted from this paper.

2.4. Model Evaluation—Accuracy Metrics

The model performance is comprehensively evaluated with other existing methods available in the literature to show improvement while predicting pressure gradient and liquid holdup. To quantify the performance of different regression models, accuracy metrics such as the Root Mean Squared Error (RMSE), Absolute Average Error (AAE), and R2 are commonly used. The mathematical explanation of some of these error metrics is provided below. The error for each prediction is defined as the difference between the measured and predicted values.
Root Mean Squared Error (RMSE)—This metric is defined as the square root of the average square distance between the measured value and the prediction. It represents the sample standard deviation of the residuals, which is a measure of how concentrated the data are around the line of best fit. It can be calculated by Equation (19).
RMSE = 1 n j = 1 n ( y j y ^ j ) 2 ,
Mean Absolute Error (MAE)—The metric is the absolute of the difference between the predicted value and observed value and provides a linear score where all the individual differences are weighted equally in the average. The metric can be calculated by Equation (20).
MAE = 1 n j = 1 n | y j y ^ j | ,
R2—Defined as the coefficient of determination, R2 provides the goodness of fit of a model in the form of an intuitive scale varying between 0 and 1. Another way to define the metric is with respect to the proportion of the explainable variance in the dependent variable with the independent variable(s). A lower R2 indicates a low level of correlation while a score of 1 indicates a perfect correlation between the two variables with no variance. Equation (21) can be used to determine the metric.
R 2 = ( y ^ j y ¯ ) 2 ( y j y ¯ ) 2 ,
Average Absolute Percentage Relative Error (AAPRE): This metric can be derived from the absolute error and is an indication of how good a prediction is with respect to the observed value. Mathematically the metric can be calculated from Equation (22) below.
AAPRE = 1 n j = 1 n | y j y ^ j y j × 100 | ,
In the equations above, y j represents the experimentally measured value for the jth data point; y ^ j   represents the predicted value from the machine learning model; and y ¯ denotes the mean of the dataset, with n being the total number of data points.

3. Results

Predictive modeling results for different regression models along with different accuracy metrics are discussed in this section. Table 1 shows the results in a tabulated format, broken down into the base model and random search cross-validated results, for the three machine learning algorithms analyzed in this work.
Based on the results for different machine learning models, XGBoost performed better when average performance was evaluated between predictions for both ϕ and fI, in terms of RMSE for the test dataset, and was selected as the final model. The learning curve after hyper-parameter tuning for interfacial friction factor is shown in Figure 6. An ‘early stopping’ approach is used to prevent overfitting. The process enables the input of the number of iterations (epochs) after which the algorithm stops if the validation score does not increase. Cross plots between the known value and its corresponding prediction are also shown in Figure 7. Hyper-parameter tuning was carried out for each machine learning algorithm, and predictions for ϕ and fI were made.
Once the predictions are made for the dataset, the error can be calculated between the experimental and predicted values.
XGBoost also provides the capability to determine the feature importance, which plays a critical role to identify predictors with the highest impact on the model. Figure 8 shows the results of feature importance for the current dataset. The metric shown in the figure is the ‘F-score’—weight, which is determined simply by how many times each feature is used to split the data across all trees.
The next step is to evaluate the model in terms of the pressure gradient and liquid holdup predictions, using the physics-based two-fluid model, as shown in the workflow in Step 6 of Figure 3. Figure 9 and Figure 10 show the comparison between experimentally determined and predictions from the proposed model for pressure drop and liquid holdup, respectively, for the entire dataset on the left and each individual data source on the right. In general, the model gives good predictions for all data sets.
The results were also compared with other existing models in the literature. The primary equations of these models are based on the two-fluid model but incorporated with different correlations for geometrical parameters and interfacial friction factor (fI). Two wetted wall fraction correlations tested were Zhang and Sarica (2011) [11], and Taitel and Dukler (1976) [3]. Churchill (1977) was used to determine wall friction factors [21]. In total, 19 different interfacial friction factor correlations from the literature were tested, including [3,5,7,8,12,55,56,57,58,59,60,61,62,63,64], which led to a total of 38 different combinations of wetted wall friction and interfacial friction factor correlations. The list was divided into different cases, referred to in Figure 11 and Figure 12, and is included in the Appendix A.
Model performance is gauged through metrics such as the average absolute relative error (AAPRE) and root mean squared error (RMSE). Figure 11 shows the performance in the form of bar charts for both metrics during the prediction of the pressure gradient. The comparison shows that the new model has the lowest average absolute relative error and RMSE.
Predictions are also made for liquid holdup, and comparisons of results using the same metric are shown in Figure 12. Based on the error metrics, the new model proposed in this work provides the lowest error in comparison to two-fluid models that use empirically derived interfacial friction factor correlations.
In addition to liquid holdup and pressure gradient predictions, the new model developed from this study was incorporated into a state-of-art mechanistic model for onset of liquid accumulation prediction proposed in [65]. That model is based on liquid film reversal of segregated flow and requires pressure gradient and liquid holdup prediction for segregated flow. The critical gas velocity determined from the coupled mechanistic model and the proposed ML algorithm from this study was further evaluated against existing droplet and film reversal models from [25,66,67,68,69,70,71,72,73,74], as listed in Table 2. It is worth mentioning that the critical gas velocity refers to the minimum gas superficial velocity that maintains segregated flow in upward inclined pipes. The experimental data used in model evaluation were extracted from [16,17,19,48,64,65,74,75,76]. Details of the mechanistic model and experimental datasets are well described in [65] and were omitted from this paper.
As expected, the liquid film models provide a better prediction than the liquid droplet models, which is consistent with the findings from previous studies [65,74]. Overall, the new model gives better predictions for critical gas velocity compared to other models evaluated in this study. This is illustrated with the lowest AARE score of 12.78 units, shown in Figure 13 between the existing models and the proposed model in this paper.

4. Summary

Segregated flow is one of the most commonly encountered flow patterns in several practical applications such as chemical, nuclear, refrigeration, and oil and gas industries. Accurate prediction of its pressure gradient and the liquid holdup is an important consideration for facility design and operations. It is also crucial for predicting the onset of liquid accumulation in the petroleum industry, which needs to be prevented and carefully managed during oil and gas production. It has been found from previous studies that the interfacial friction factor plays a crucial role in two-fluid model performance. The previous approach is to use empirically derived correlations that are only applicable to the small subset of working conditions under which those solutions are derived.
Using predictions from a mechanistic model augmented with data-driven analysis, a novel and more generalized approach is proposed that combines a physics-based two-fluid model and machine learning algorithms to accurately predict pressure drop and liquid holdup. The new model is comprehensively evaluated with existing models with various closure relationships and shows a significant improvement in terms of prediction accuracy. A comparison of the average absolute percentage error was also made for critical gas velocity between the experimentally determined values and model predictions, which include both droplet and film-reversal based models. The hybrid model gives the best prediction for the critical gas velocity compared with other existing modeling approaches.
The proposed modeling workflow helps to reduce the dependence on empirically derived correlations and the need for interfacial friction factor correlation selection. Considering the gradual transition to the digital oil field, more field data will be available in the future, which can be used to augment the physics-based mechanistic multiphase flow modeling. The proposed approach in this study will add significant value in segregated flow modeling in terms of combining physics-based and data-driven models, and therefore the optimal design of production systems.

5. Discussion

A hybrid machine learning and mechanistic model to determine the pressure gradient and liquid holdup, as well as the critical gas velocity, presented in this work, provides better accuracy in comparison to studies derived from literature, which are primarily mechanistic or empirical-based. This has been made possible by leveraging the predictive capabilities of a machine learning algorithm along with the causal understanding of fluid flow mechanisms through a mechanistic model. However, this conclusion is accompanied by a caveat that significant, good-quality data are available before the implementation of any machine learning model. It is important to understand that machine learning and mechanistic modeling are two different paradigms and should never be considered to compete with or replace one other.
According to Baker (2018) [77], if large-scale datasets are available, machine learning algorithms can be an efficient and scalable modeling approach, as well as provide the ability to avoid the need to understand complex mechanisms. Advanced algorithms in the field of machine learning have been proven to identify hidden correlations between parameters that are not easy to identify through conventional approaches. One of the biggest drawbacks of a pure data-driven model is the applicability and generalization of using the model outside the boundary conditions in which the model was trained. This is where a mechanistic model has an advantage since the predictions from a mechanistic model, based on fundamental laws of nature, can be applied under conditions where experiments are either costly or difficult to perform. As discussed earlier, a mechanistic model relies on the generation of novel hypotheses and is usually built to mimic real-life events with certain assumptions. In some cases, these oversimplified assumptions and their extreme specific nature, commonly found in empirical models, prevent the model predictions to be applicable on a wider scale. This paper shows how a symbiotic relationship can be created between machine learning and mechanistic modeling approaches by harnessing the positive aspects of both approaches.

Author Contributions

Conceptualization, A.R. and Y.F.; methodology, A.R. and Y.F.; software, A.R. and Y.F.; validation, A.R. and Y.F.; formal analysis, A.R. and Y.F.; investigation, A.R. and Y.F.; resources, A.R. and Y.F.; data curation, A.R. and Y.F.; writing—original draft preparation, A.R.; writing—review and editing, Y.F.; supervision, Y.F.; project administration, Y.F.; funding acquisition, Y.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used in the model development and evaluation were obtained from open literature, and the references are provided in the Reference section. The software packages developed in this study and collected data are proprietary and cannot be released.

Acknowledgments

We thank Yilin Fan’s research group for the financial support of this project.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

ACGas core occupied cross-sectional area, m2
ALLiquid film occupied cross-sectional area, m2
ApPipe cross-sectional area, m2
dPipe inner diameter, m
dLLiquid hydraulic diameter, m
dCGas core hydraulic diameter, m
(dp/dL)totalTotal pressure gradient, Pa/m
(dp/dL)LLiquid phase pressure gradient, Pa/m
(dp/dL)CGas core pressure gradient, Pa/m
(dp/dL)fFrictional component of pressure gradient, Pa/m
(dp/dL)gGravitational component of pressure gradient, Pa/m
fCGas wall friction factor, -
fG-SPGas wall friction factor for single phase pipe flow, -
fIInterfacial friction factor, -
fLLiquid wall friction factor, -
fL-SPLiquid wall friction factor for single phase pipe flow, -
FEEntrainment fraction, -
gGravitational acceleration, m/s2
hLLiquid film thickness, m
HLTotal liquid holdup, -
ReReynolds number, -
SCGas wetted perimeter, m
SIInterfacial perimeter, m
SLLiquid wetted perimeter, m
vLAverage liquid velocity, m/s
vCAverage gas core velocity, m/s
vSgGas superficial velocity, m/s
vSLLiquid superficial velocity, m/s
εPipe roughness, m
θInclination angle from horizontal, ° (degree)
μGGas viscosity, Pa·s
μLLiquid viscosity, Pa·s
ρCGas core density, kg/m3
ρgGas density, kg/m3
ρLLiquid density, kg/m3
σSurface tension, N/m
τWCGas core wall shear stress, Pa
τIInterfacial shear stress, Pa
τWLAverage liquid wall shear stress, Pa
ϕLiquid wall friction factor coefficient, -
AAEAbsolute Average Error
AAPREAverage Absolute Percent Relative Error
ANNArtificial Neural Network
CVCross-Validation
EDAExploratory Data Analysis
IQRInter-Quartile Range
MAEMean Absolute Error
MLMachine Learning
PDEPartial Differential Equation
RFRandom Forest
RMSERoot Mean Squared Error
XGBoosteXtreme Gradient Boosting

Appendix A

Table indicating case number and the corresponding combination of geometrical parameters, wall friction factor, and interfacial friction factor used for model validation.
Table A1. Case number and the corresponding correlations used for model validation in Figure 11 and Figure 12.
Table A1. Case number and the corresponding correlations used for model validation in Figure 11 and Figure 12.
Case #Geometrical ParameterWall Friction Factor, fLInterfacial Friction Factor, fI
New ModelZhang and Sarica (2011)New ML Model from Current Study New ML Model from Current Study
Case 02Taitel and Dukler (1976)Churchill (1977)Cohen and Hanratty (1968)
Case 03Taitel and Dukler (1976)Churchill (1977)Hart et al. (1989)
Case 04Taitel and Dukler (1976)Churchill (1977)Kowalski (1987)
Case 05Taitel and Dukler (1976)Churchill (1977)Taitel and Dukler (1976)
Case 06Taitel and Dukler (1976)Churchill (1977)Vlachos et al. (1997)
Case 07Taitel and Dukler (1976)Churchill (1977)Wallis (1969)
Case 08Taitel and Dukler (1976)Churchill (1977)Wallis Modified (1969)
Case 09Taitel and Dukler (1976)Churchill (1977)Whalley and Hewitt (1978)
Case 10Taitel and Dukler (1976)Churchill (1977)Oliemans et al. (1986)
Case 11Taitel and Dukler (1976)Churchill (1977)Fore et al. (2000)
Case 12Taitel and Dukler (1976)Churchill (1977)Dallman et al. (1979)
Case 13Taitel and Dukler (1976)Churchill (1977)Ambrosini et al. (1991)
Case 14Taitel and Dukler (1976)Churchill (1977)Hamersma and Hart (1987)
Case 15Taitel and Dukler (1976)Churchill (1977)Chen et al. (1997)
Case 16Taitel and Dukler (1976)Churchill (1977)Andritsos and Hanratty (1987)
Case 17Taitel and Dukler (1976)Churchill (1977)Andritsos et al. (2008)
Case 18Taitel and Dukler (1976)Churchill (1977)Zhang et al. (2003)
Case 19Taitel and Dukler (1976)Churchill (1977)Brito (2015)
Case 20Taitel and Dukler (1976)Churchill (1977)Grolman and Fortuin (1997)
Case 21Zhang and Sarica (2011)Churchill (1977)Cohen and Hanratty (1968)
Case 22Zhang and Sarica (2011)Churchill (1977)Hart et al. (1989)
Case 23Zhang and Sarica (2011)Churchill (1977)Kowalski (1985)
Case 24Zhang and Sarica (2011)Churchill (1977)Taitel and Dukler (1976)
Case 25Zhang and Sarica (2011)Churchill (1977)Vlachos et al. (1997)
Case 26Zhang and Sarica (2011)Churchill (1977)Wallis (1997)
Case 27Zhang and Sarica (2011)Churchill (1977)Wallis Modified (1969)
Case 28Zhang and Sarica (2011)Churchill (1977)Whalley and Hewitt (1978)
Case 29Zhang and Sarica (2011)Churchill (1977)Oliemans et al. (1986)
Case 30Zhang and Sarica (2011)Churchill (1977)Fore et al. (2000)
Case 31Zhang and Sarica (2011)Churchill (1977)Dallman et al. (1979)
Case 32Zhang and Sarica (2011)Churchill (1977)Ambrosini et al. (1991)
Case 33Zhang and Sarica (2011)Churchill (1977)Hamersma and Hart (1987)
Case 34Zhang and Sarica (2011)Churchill (1977)Chen et al. (1997)
Case 35Zhang and Sarica (2011)Churchill (1977)Andritsos and Haratty (1987)
Case 36Zhang and Sarica (2011)Churchill (1977)Andritsos et al. (2008)
Case 37Zhang and Sarica (2011)Churchill (1977)Zhang et al. (2003)
Case 38Zhang and Sarica (2011)Churchill (1977)Brito (2015)
Case 39Zhang and Sarica (2011)Churchill (1977)Grolman and Fortuin (1997)

References

  1. Ullmann, A.; Brauner, N. Closure Relations for Two-fluid Models for Two-phase Stratified Smooth and Stratified Wavy Flows. Int. J. Multiph. Flow 2006, 32, 82–105. [Google Scholar] [CrossRef]
  2. Shoham, O. Mechanistic Modeling of Gas-Liquid Two-Phase Flow in Pipes; Society of Petroleum Engineers: Houston, TX, USA, 2006. [Google Scholar]
  3. Taitel, Y.; Dukler, A.E. A Model for Predicting Flow Regime Transitions in Horizontal and Near Horizontal Gas-liquid Flow. AIChE J. 1976, 22, 47–55. [Google Scholar] [CrossRef]
  4. Shoham, O.; Taitel, Y. Stratified Turbulent-turbulent Gas-liquid Flow in Horizontal and Inclined Pipes. AIChE J. 1984, 30, 377–385. [Google Scholar] [CrossRef]
  5. Kowalski, J.E. Wall and Interfacial Shear Stress in Stratified Flow in a Horizontal Pipe. AIChE J. 1987, 33, 274–281. [Google Scholar] [CrossRef]
  6. Andritsos, N.; Hanratty, T. Interfacial Instabilities for Horizontal Gas-liquid Flows in Pipelines. Int. J. Multiph. Flow 1987, 13, 583–603. [Google Scholar] [CrossRef]
  7. Andritsos, N.; Hanratty, T.J. Influence of Interfacial Waves in Stratified Gas-liquid Flows. AIChE J. 1987, 33, 444–454. [Google Scholar] [CrossRef]
  8. Hart, J.; Hamersma, P.; Fortuin, J. Correlations Predicting Frictional Pressure Drop and Liquid Holdup during Horizontal Gas-liquid Pipe Flow with a Small Liquid Holdup. Int. J. Multiph. Flow 1989, 15, 947–964. [Google Scholar] [CrossRef]
  9. Hamersma, P.; Hart, J. A Pressure Drop Correlation for Gas/liquid Pipe Flow with a Small Liquid Holdup. Chem. Eng. Sci. 1987, 42, 1187–1196. [Google Scholar] [CrossRef]
  10. Chen, X.T.; Cal, X.D.; Brill, J.P. Gas-Liquid Stratified-Wavy Flow in Horizontal Pipelines. J. Energy Resour. Technol. 1997, 119, 209–216. [Google Scholar] [CrossRef]
  11. Zhang, H.-Q.; Sarica, C. A Model for Wetted-Wall Fraction and Gravity Center of Liquid Film in Gas/Liquid Pipe Flow. SPE J. 2011, 16, 692–697. [Google Scholar] [CrossRef]
  12. Oliemans, R.; Pots, B.; Trompé, N. Modelling of Annular Dispersed Two-phase Flow in Vertical Pipes. Int. J. Multiph. Flow 1986, 12, 711–732. [Google Scholar] [CrossRef]
  13. Alves, I.; Caetano, E.; Minami, K.; Shoham, O. Modeling Annular Flow Behavior for Gas Wells. SPE Prod. Eng. 1991, 6, 435–440. [Google Scholar] [CrossRef]
  14. Gazley, C. Interfacial Shear and Stability in Two-Phase Flow. Ph.D. Thesis, University of Delaware, Newark, DE, USA, 1948. [Google Scholar]
  15. Rastogi, A. Onset of Liquid Loading in Large Diameter Inclined Pipes. Ph.D. Thesis, Colorado School of Mines, Golden, CO, USA, 2020. [Google Scholar]
  16. Alsaadi, Y.; Pereyra, E.; Torres, C.; Sarica, C. Liquid loading of highly deviated gas wells from 60° to 88°. In Proceedings of the SPE Annual Technical Conference and Exhibition, Houston, TX, USA, 28–30 September 2015. [Google Scholar] [CrossRef]
  17. Guner, M.; Pereyra, E.; Sarica, C.; Torres, C.F. An Experimental Study of Low Liquid Loading in Inclined Pipes from 90° to 45°. In Proceedings of the SPE Production and Operations Symposium, Oklahoma City, OK, USA, 1–5 March 2015. [Google Scholar] [CrossRef]
  18. Fan, Y. A Study of the Onset of Liquid Accumulation and Pseudo-Slug Flow Chracterization. Ph.D. Thesis, University of Tulsa, Tulsa, OK, USA, 2017. [Google Scholar]
  19. Rodrigues, H.T.; Pereyra, E.; Sarica, C. Pressure Effects on Low-Liquid-Loading Oil/Gas Flow in Slightly Upward Inclined Pipes: Flow Pattern, Pressure Gradient, and Liquid Holdup. SPE J. 2019, 24, 2221–2238. [Google Scholar] [CrossRef]
  20. Blasius, H. Grenzschichten in Flüssigkeiten mit kleiner Reibung. Zeits. Math. Phys. 1907, 56, 1–37. [Google Scholar]
  21. Churchill, S.W. Friction-Factor Equation Spans All Fluid-Flow Regimes. Chem. Eng. J. 1977, 84, 91–92. [Google Scholar]
  22. Hall, N.A. Thermodynamics of Fluid Flow; Longmans Green: New York, NY, USA, 1957. [Google Scholar]
  23. Massey, B.S. Mechanics of Fluids; Chapman and Hall: Boston, MA, USA, 1989. [Google Scholar]
  24. Colebrook, C.F. Turbulent Flow in Pipes, with Particular Reference to the Transition Region between the Smooth and Rough Pipe Laws. J. Inst. Civ. Eng. 1939, 11, 133–156. [Google Scholar] [CrossRef]
  25. Barnea, D. Transition from Annular Flow and from Dispersed Bubble Flow—Unified Models for the Whole Range of Pipe Inclinations. Int. J. Multiph. Flow 1986, 12, 733–744. [Google Scholar] [CrossRef]
  26. Kanin, E.; Osiptsov, A.; Vainshtein, A.; Burnaev, E. A Predictive Model for Steady-state Multiphase Pipe Flow: Machine Learning on Lab Data. J. Pet. Sci. Eng. 2019, 180, 727–746. [Google Scholar] [CrossRef] [Green Version]
  27. Buist, J.; Sanderse, B.; Van Halder, Y.; Koren, B.; Van Heijst, G. Machine Learning for Closure Models in Multiphase Flow Applications. In Proceedings of the 3rd International Conference on Uncertainty Quantification in Computational Sciences and Engineering (UNCECOMP 2019), Crete, Greece, 24–26 June 2019. [Google Scholar] [CrossRef] [Green Version]
  28. Mi, Y.; Ishii, M.; Tsoukalas, L. Flow Regime Ddentification Methodology with Neural Networks and Two-phase Flow Models. Nucl. Eng. Des. 2001, 204, 87–100. [Google Scholar] [CrossRef]
  29. Shaban, H.; Tavoularis, S. Measurement of Gas and Liquid Flow Rates in Two-phase Pipe Flows by the Application of Machine Learning Techniques to Differential Pressure Signals. Int. J. Multiph. Flow 2014, 67, 106–117. [Google Scholar] [CrossRef]
  30. Osman, E.-S.A. Artificial Neural Network Models for Identifying Flow Regimes and Predicting Liquid Holdup in Horizontal Multiphase Flow. SPE Prod. Facil. 2004, 19, 33–40. [Google Scholar] [CrossRef]
  31. Li, X.; Miskimins, J.L.; Sutton, R.P.; Hoffman, B.T. Multiphase Flow Pattern Recognition in Horizontal and Upward Gas-Liquid Flow Using Support Vector Machine Models. In Proceedings of the SPE Annual Technical Conference and Exhibition, Amsterdam, The Netherlands, 27–29 October 2014. [Google Scholar] [CrossRef]
  32. Ezzatabadipour, M.; Singh, P.; Robinson, M.D.; Guillen-Rondon, P.; Torres, C. Deep Learning as a Tool to Predict Flow Patterns in Two-Phase Flow. arXiv 2017, arXiv:1705.07117. Available online: Https://arxiv.org/pdf/1705.07117 (accessed on 1 August 2021). [Google Scholar]
  33. Guillen-Rondon, P.; Robinson, M.D.; Torres, C.; Pereya, E. Support Vector Machine Application for Multiphase Flow Pattern Prediction. arXiv 2018, arXiv:1806.05054. Available online: Http://arxiv.org/abs/1806.05054 (accessed on 1 August 2021). [Google Scholar]
  34. Mask, G.; Wu, X.; Ling, K. An Improved Model for Gas-liquid Flow Pattern Prediction Based on Machine Learning. J. Pet. Sci. Eng. 2019, 183, 106370. [Google Scholar] [CrossRef]
  35. Hernandez, J.S.; Valencia, C.; Ratkovich, N.; Torres, C.F.; Muñoz, F. Data Driven Methodology for Model Selection in Flow Pattern Prediction. Heliyon 2019, 5, e02718. [Google Scholar] [CrossRef] [Green Version]
  36. Seong, Y.; Park, C.; Choi, J.; Jang, I. Surrogate Model with a Deep Neural Network to Evaluate Gas–Liquid Flow in a Horizontal Pipe. Energies 2020, 13, 968. [Google Scholar] [CrossRef] [Green Version]
  37. Mohammadi, S. A Machine Learning Model to Find a Near-Optimal Set of Closure Relationships for the Multiphase Unified Model. Ph.D. Thesis, The University of Tulsa, Tulsa, OK, USA, 2020. [Google Scholar]
  38. Lin, Z.; Liu, X.; Lao, L.; Liu, H. Prediction of Two-phase Flow Patterns in Upward Inclined Pipes via Deep Learning. Energy 2020, 210, 118541. [Google Scholar] [CrossRef]
  39. Mukherjee, H. An Experimental Study of Inclined Two-Phase Flow. Master’s Thesis, The University of Tulsa, Tulsa, OK, USA, 1979. [Google Scholar]
  40. Magrini, K.L. Liquid Entrainment in Annular Gas-liquid Flow in Inclined Pipes. Master’s Thesis, The University of Tulsa, Tulsa, OK, USA, 2009. [Google Scholar]
  41. Yuan, G. Liquid Loading of Gas Wells. Master’s Thesis, The University of Tulsa, Tulsa, OK, USA, 2011. [Google Scholar]
  42. Fan, Y. An Investigation of Low Liquid Loading Gas-Liquid Stratified Flow in Near-Horizontal Pipes. Ph.D. Thesis, The University of Tulsa, Tulsa, OK, USA, 2005. [Google Scholar]
  43. Gokcal, B. An Experimental and Theoretical Investigation of Slug Flow for High Oil Viscosity in Horizontal Pipes. Ph.D. Thesis, The University of Tulsa, Tulsa, OK, USA, 2008. [Google Scholar]
  44. Brito, R. Effect of Medium Oil Viscosity on Two-Phase Oil-Gas Flow Behavior in Horizontal Pipes. Master’s Thesis, The University of Tulsa, Tulsa, OK, USA, 2012. [Google Scholar]
  45. Brill, J.P.; Chen, X.T.; Flores, J.G.; Marcano, R. Transportation of Liquids in Multiphase Pipelines under Low Liquid Loading Conditions; Final Report; The Pennsylvania State University: University Park, PA, USA, 1995. [Google Scholar]
  46. Guner, M. Liquid Loading of Gas Wells with Deviations from 0° to 45°. Master’s Thesis, The University of Tulsa, Tulsa, OK, USA, 2012. [Google Scholar]
  47. Meng, W. Low Liquid Loading Gas-Liquid Two-Phase Flow in Near-Horizontal Pipes. Ph.D. Thesis, The University of Tulsa, Tulsa, OK, USA, 1999. [Google Scholar]
  48. Langsholt, M.; Holm, H. Liquid Accumulation in Gas-Condensate Pipelines—An Experimental Study. In Proceedings of the 13th International Conference on Multiphase Production Technology, Edinburgh, UK, 13–15 June 2007. [Google Scholar]
  49. Exploratory Data Analysis. The Concise Encyclopedia of Statistics; Springer: New York, NY, USA, 2008; pp. 192–194. [Google Scholar]
  50. Seo, S. A Review and Comparison of Methods for Detecting Outliers in Univariate Datasets. Master’s Thesis, University of Pittsburgh, Pittsburgh, PA, USA, 2006. [Google Scholar]
  51. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  52. Chen, T.; Guestrin, C. XGBoost: A Scalable Tree Boosting System. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; ACM: New York, NY, USA, 2016; Volume 10, pp. 785–794. [Google Scholar]
  53. Géron, A. Hands-On Machine Learning with Scikit-Learn, Keras, and TensorFlow: Concepts, Tools, and Techniques to Build Intelligent Systems; O’Reilly Media: Newton, MA, USA, 2019. [Google Scholar]
  54. Pandey, Y.N.; Rastogi, A.; Kainkaryam, S.; Bhattacharya, S.; Saputelli, L. Machine Learning in the Oil and Gas Industry: Including Geosciences, Reservoir Engineering, and Production Engineering with Python; Apress: New York, NY, USA, 2020. [Google Scholar]
  55. Cohen, L.S.; Hanratty, T.J. Effect of Waves at a Gas—Liquid Interface on a Turbulent Air Flow. J. Fluid Mech. 1968, 31, 467–479. [Google Scholar] [CrossRef]
  56. Vlachos, N.; Paras, S.; Karabelas, A. Liquid-to-wall Shear Stress Distribution in Stratified/atomization Flow. Int. J. Multiph. Flow 1997, 23, 845–863. [Google Scholar] [CrossRef]
  57. Wallis, G.B. One-Dimensional Two-Phase Flow; McGraw-Hill: New York, NY, USA, 1969. [Google Scholar]
  58. Whalley, P.B.; Hewitt, G.F. The Correlation of Liquid Entrainment Fraction and Entrainment Rate in Annular Two-Phase Flow; UKAEA Atomic Energy Research Establishment: Harwell, UK, 1978. [Google Scholar]
  59. Fore, L.; Beus, S.; Bauer, R. Interfacial Friction in Gas–liquid Annular Flow: Analogies to Full and Transition Roughness. Int. J. Multiph. Flow 2000, 26, 1755–1769. [Google Scholar] [CrossRef] [Green Version]
  60. Dallman, J.C.; Jones, B.G.; Hanratty, T.J. Interpretation of Entrainment Measurements in Annular Gas-Liquid Flows. In Two-phase Momentum, Heat and Mass Transfer in Chemical, Process, and Energy Engineering Systems; Hemisphere Publishing Corporation: London, UK, 1979; Volume 2, pp. 681–693. [Google Scholar]
  61. Ambrosini, W.; Andreussi, P.; Azzopardi, B. A Physically Based Correlation for Drop Size in Annular Flow. Int. J. Multiph. Flow 1991, 17, 497–507. [Google Scholar] [CrossRef]
  62. Andritsos, N.; Tzotzi, C.; Hanratty, T.J. Interfacial Shear Stress in Wavy Stratified Gas-liquid Two-phase Flow. In Proceedings of the 5th European Thermal-Sciences Conference, Eindhoven, The Netherlands, 18–22 May 2008. [Google Scholar]
  63. Zhang, H.-Q.; Wang, Q.; Sarica, C.; Brill, J.P. Unified Model for Gas-Liquid Pipe Flow via Slug Dynamics—Part 1: Model Development. J. Energy Resour. Technol. 2003, 125, 266–273. [Google Scholar] [CrossRef]
  64. Brito, R.M. Effect of Horizontal Well Trajectory on Two-Phase Gas-Liquid Flow Behavior. Ph.D. Theis, The University of Tulsa, Tulsa, OK, USA, 2015. [Google Scholar]
  65. Rastogi, A.; Fan, Y. Experimental and Modeling Study of Onset of Liquid Accumulation. J. Nat. Gas Sci. Eng. 2020, 73, 103064. [Google Scholar] [CrossRef]
  66. Turner, R.G.; Hubbard, M.G.; Dukler, A.E. Analysis and Prediction of Minimum Flow Rate for the Continuous Removal of Liquids from Gas Wells. J. Pet. Technol. 1969, 21, 1475–1482. [Google Scholar] [CrossRef]
  67. Coleman, S.B.; Clay, H.B.; McCurdy, D.G.; Norris, L.H.I. A New Look at Predicting Gas-Well Load-Up. J. Pet. Technol. 1991, 43, 329–333. [Google Scholar] [CrossRef]
  68. Li, M.; Lei, S.; Li, S. New View on Continuous-Removal Liquids from Gas Wells. In Proceedings of the SPE Permian Basin Oil and Gas Recovery Conference, Midland, TX, USA, 15–17 May 2001. [Google Scholar]
  69. Wang, Y.Z.; Liu, Q.W. A New Method to Calculate the Minimum Critical Liquid Carrying Flow Rate for Gas Wells. Pet. Geol. Oilfield Dev. Daqing 2007, 6, 82–85. [Google Scholar]
  70. Belfroid, S.; Schiferli, W.; Alberts, G.; Veeken, C.A.M.; Biezen, E. Predicting Onset and Dynamic Behavior of Liquid Loading Gas Wells. In Proceedings of the SPE Annual Technical Conference and Exhibition, Denver, CO, USA, 21–24 September 2018. [Google Scholar] [CrossRef]
  71. Zhou, D.; Yuan, H. A New Model for Predicting Gas-Well Liquid Loading. SPE Prod. Oper. 2010, 25, 172–181. [Google Scholar] [CrossRef]
  72. Luo, S.; Kelkar, M.; Pereyra, E.; Sarica, C. A New Comprehensive Model for Predicting Liquid Loading in Gas Wells. SPE Prod. Oper. 2014, 29, 337–349. [Google Scholar] [CrossRef]
  73. Shekhar, S.; Kelkar, M.; Hearn, W.J.; Hain, L.L. Improved Prediction of Liquid Loading in Gas Wells. SPE Prod. Oper. 2017, 32, 539–550. [Google Scholar] [CrossRef]
  74. Fan, Y.; Pereyra, E.; Sarica, C. Onset of Liquid-Film Reversal in Upward-Inclined Pipes. SPE J. 2018, 23, 1630–1647. [Google Scholar] [CrossRef]
  75. Nair, J.V. Gas Lift Application and Severe Slugging in Toe-Down Horizontal Wells. Ph.D. Thesis, The University of Tulsa, Tulsa, OK, USA, 2017. [Google Scholar]
  76. Espedal, M. An Experimental Investigation of Stratified Two-Phase Pipe Flow at Small Inclinations. Ph.D. Thesis, Norwegian University of Science and Technology (NTNU), Trondheim, Norway, 1998. [Google Scholar]
  77. Baker, R.E.; Peña, J.-M.; Jayamohan, J.; Jérusalem, A. Mechanistic Models versus Machine Learning, a Fight worth Fighting for the Biological Community? Biol. Lett. 2018, 14, 20170660. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Parameters and forces in two-fluid model for segregated flow.
Figure 1. Parameters and forces in two-fluid model for segregated flow.
Fluids 07 00012 g001
Figure 2. Comparison of experimentally measured frictional pressure gradient with calculated values from correlations for (a) different datasets with various flowing conditions and (b) various wall friction correlations.
Figure 2. Comparison of experimentally measured frictional pressure gradient with calculated values from correlations for (a) different datasets with various flowing conditions and (b) various wall friction correlations.
Fluids 07 00012 g002
Figure 3. Workflow of model development (Steps 1–5) and model implementation (Step 5).
Figure 3. Workflow of model development (Steps 1–5) and model implementation (Step 5).
Fluids 07 00012 g003
Figure 4. Univariate distribution plot for the entire dataset.
Figure 4. Univariate distribution plot for the entire dataset.
Fluids 07 00012 g004
Figure 5. Comparison of the distribution for ϕ (top) and fI (bottom) before (left column) and after (right column) the outliers were removed.
Figure 5. Comparison of the distribution for ϕ (top) and fI (bottom) before (left column) and after (right column) the outliers were removed.
Fluids 07 00012 g005
Figure 6. Learning curve for the XGBoost model with respect to RMSE for 1000 epochs using Randomized Search CV from scikit-learn.
Figure 6. Learning curve for the XGBoost model with respect to RMSE for 1000 epochs using Randomized Search CV from scikit-learn.
Fluids 07 00012 g006
Figure 7. Prediction performance using XGBoost as an algorithm for ϕ and fI.
Figure 7. Prediction performance using XGBoost as an algorithm for ϕ and fI.
Fluids 07 00012 g007
Figure 8. Feature importance of predictors determined using XGBoost.
Figure 8. Feature importance of predictors determined using XGBoost.
Fluids 07 00012 g008
Figure 9. Comparison of experimental measurement and predictions for the pressure gradient.
Figure 9. Comparison of experimental measurement and predictions for the pressure gradient.
Fluids 07 00012 g009
Figure 10. Comparison of experimental measurement and predictions for liquid holdup.
Figure 10. Comparison of experimental measurement and predictions for liquid holdup.
Fluids 07 00012 g010
Figure 11. Model evaluation for pressure gradient predictions in terms of average absolute relative error and RMSE.
Figure 11. Model evaluation for pressure gradient predictions in terms of average absolute relative error and RMSE.
Fluids 07 00012 g011
Figure 12. Model evaluation for liquid holdup predictions in terms of average absolute relative error and RMSE.
Figure 12. Model evaluation for liquid holdup predictions in terms of average absolute relative error and RMSE.
Fluids 07 00012 g012
Figure 13. Average absolute relative error for different model predictions of critical gas velocity.
Figure 13. Average absolute relative error for different model predictions of critical gas velocity.
Fluids 07 00012 g013
Table 1. Comparison of the accuracy metric for different machine learning models.
Table 1. Comparison of the accuracy metric for different machine learning models.
Model and Response VariablesDataMAER2RMSEDataR2RMSE
RF (ϕ)Base ModelTrain0.0760.9490.129Entire Data0.8680.207
Test0.2150.6200.348
Random Search CVTrain0.0520.9720.094Entire Data0.8980.182
Test0.1920.6690.325
RF (fI)Base ModelTrain0.00150.9490.0028Entire Data0.8850.0043
Test0.00370.7040.0069
Random Search CVTrain0.00100.9720.0020Entire Data0.9010.0039
Test0.00370.7030.0070
XGB (ϕ)Base ModelTrain0.2070.7730.304Entire Data0.7410.3270
Test0.2580.6510.386
Random Search CVTrain0.0720.9740.101Entire Data0.9100.1920
Test0.2010.7290.340
XGB (fI)Base ModelTrain0.00390.8000.0057Entire Data0.7670.0063
Test0.00500.6860.0079
Random Search CVTrain0.00130.9790.0018Entire Data0.9190.0037
Test0.00380.7760.0067
ANN (ϕ)Base Model Train0.3020.4250.436Entire Data0.4160.436
Test0.2950.3800.439
Grid Search CVTrain0.2800.4480.427Entire Data0.4390.428
Test0.2840.4000.432
ANN (fI)Base ModelTrain0.00480.6650.0073Entire Data0.6480.0074
Test0.00520.5930.0076
Grid Search CVTrain0.00440.6980.0069Entire Data0.6830.007
Test0.00470.6310.0073
Table 2. Start-of-the-art droplet and film reversal-based modeling studies for critical gas velocity estimation.
Table 2. Start-of-the-art droplet and film reversal-based modeling studies for critical gas velocity estimation.
No.Models Based on Liquid Droplet RemovalNo.Models Based on Liquid Film
Dynamics/Reversal
1Turner et al. (1969)7Barnea (1986)
2Coleman et al. (1991)8Luo et al. (2014)
3Li et al. (2002)9Shekhar et al. (2017)
4Wang and Liu (2007)10Fan et al. (2018)
5Belfroid et al. (2008)11New Model (Rastogi and Fan (2019) with new ML algorithm from this study)
6Zhou and Yuan (2010)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rastogi, A.; Fan, Y. Machine Learning Augmented Two-Fluid Model for Segregated Flow. Fluids 2022, 7, 12. https://doi.org/10.3390/fluids7010012

AMA Style

Rastogi A, Fan Y. Machine Learning Augmented Two-Fluid Model for Segregated Flow. Fluids. 2022; 7(1):12. https://doi.org/10.3390/fluids7010012

Chicago/Turabian Style

Rastogi, Ayush, and Yilin Fan. 2022. "Machine Learning Augmented Two-Fluid Model for Segregated Flow" Fluids 7, no. 1: 12. https://doi.org/10.3390/fluids7010012

APA Style

Rastogi, A., & Fan, Y. (2022). Machine Learning Augmented Two-Fluid Model for Segregated Flow. Fluids, 7(1), 12. https://doi.org/10.3390/fluids7010012

Article Metrics

Back to TopTop