Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Influence of Front-End Electronics on Metrological Performance of QCM Systems
Previous Article in Journal
Cam-Unet: Print-Cam Image Correction for Zero-Bit Fourier Image Watermarking
Previous Article in Special Issue
Robust Vector BOTDA Signal Processing with Probabilistic Machine Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Gas–Liquid Two-Phase Slug Flow Using Distributed Acoustic Sensing in Horizontal Pipes

by
Sharifah Ali
1,
Ge Jin
2 and
Yilin Fan
1,*
1
Petroleum Engineering Department, Colorado School of Mines, Golden, CO 80401, USA
2
Geophysics Department, Colorado School of Mines, Golden, CO 80401, USA
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(11), 3402; https://doi.org/10.3390/s24113402
Submission received: 24 April 2024 / Revised: 19 May 2024 / Accepted: 22 May 2024 / Published: 25 May 2024
(This article belongs to the Special Issue Advances in Fiber Optic Sensors for Energy Applications)

Abstract

:
This article discusses the use of distributed acoustic sensing (DAS) for monitoring gas–liquid two-phase slug flow in horizontal pipes, using standard telecommunication fiber optics connected to a DAS integrator for data acquisition. The experiments were performed in a 14 m long, 5 cm diameter transparent PVC pipe with a fiber cable helically wrapped around the pipe. Using mineral oil and compressed air, the system captured various flow rates and gas–oil ratios. New algorithms were developed to characterize slug flow using DAS data, including slug frequency, translational velocity, and the lengths of slug body, slug unit, and the liquid film region that had never been discussed previously. This study employed a high-speed camera next to the fiber cable sensing section for validation purposes and achieved a good correlation among the measurements under all conditions tested. Compared to traditional multiphase flow sensors, this technology is non-intrusive and offers continuous, real-time measurement across long distances and in harsh environments, such as subsurface or downhole conditions. It is cost-effective, particularly where multiple measurement points are required. Characterizing slug flow in real time is crucial to many industries that suffer slug-flow-related issues. This research demonstrated the DAS’s potential to characterize slug flow quantitively. It will offer the industry a more optimal solution for facility design and operation and ensure safer operational practices.

1. Introduction

Multiphase flow in pipes is a common but also complex phenomenon occurring in various fields, such as petroleum engineering, geothermal applications, nuclear engineering, etc. For example, in the petroleum industry, oil and gas flow is inevitable due to the pressure reduction and temperature variations along the wellbore and pipeline systems. As energy demand increases and the world explores deeper harsh zones to increase oil and gas recovery, multiphase-flow-related safety issues are becoming more complex and concerning. Critical parameter monitoring in real time is essential for better production management and optimization. Detecting and forecasting the circumstances of the well at earlier periods have a substantial effect on well control strategies; health, safety, and environment (HSE); and risk management [1]. This will give the operators the ability to take the appropriate action at the right time [2].
There are several types of point sensors that can characterize slug flow, such as γ-ray sensors, wire-mesh sensors, electrical impedance (capacitance and conductance) sensors, microwave sensors, optical sensors, etc. They are either not safe (γ-ray sensors have a radioactive source), intrusive (such as the wire-mesh sensors), or have special requirements on the type of fluids or flow pattern (see details in [3]). One significant drawback is that they all require installation directly into the piping system, necessitating the cutting of pipes and their subsequent connections which adds another failure factor. Moreover, these technologies can only provide point measurements at a specific location where the sensor is installed. Considering the cost of each sensor, it is not economically feasible to install one for each individual pipeline/well that requires detailed monitoring.
The fiber-optic sensing technology is developing rapidly and signifies the future of multiphase flow supervising and monitoring. Fiber optics can prevail over the limitations of conventional electrical sensor arrays due to their smaller size, non-intrusive features, corrosion resistance, isolation, and ability to function effectively in intense environmental situations [4]. The small size of these sensors enables them to be securely utilized over extended distances with minimized future maintenance operations [5]. One of the crucial benefits of fiber optics is that the same sensor can carry dual purposes. It can act in the same way as the sensing element for measuring the critical physical parameters and as a transmission medium for the detected signal [6].
Generally, distributed fiber-optic sensing (DFOS) can be classified into three types: (a) distributed temperature sensing (DTS), (b) distributed acoustic sensing (DAS), and (c) distributed strain sensing (DSS) [7]. DAS employs an optical time-domain reflectometer (OTDR) system to record high-frequency strain rates over long distances. The sensing range along the fiber length of a DAS system can be more than 10 km, with a spatial resolution less than 1 m and a temporal sampling frequency up to 10 kHz [8]. Since a well-protected sensing fiber can endure harsh environments, DAS is desirable in wellbore diagnostics to examine several aspects. It can serve as a noise log to detect leakage or estimate injection rate—e.g., [9,10], supporting the differentiation between the channeling, borehole flow, matrix flow, and fracture flow [11]. It can also be used for hydraulic fracture monitoring [12] and seismic acquisition [13]. DFOS is commonly used and has many applications in engineering and geophysics disciplines.
van der Horst et al. (2013) [14] reported using DAS for flow monitoring for tight gas well producers. The acoustic noise was monitored, and the production profiles were generated by the gas phase flowing from the perforations, which is then converted to flowrates using the DAS signal in a specified frequency band. This method has provided some quantitative production measurement but does not apply to oil well producers. A study shows the flow velocity estimation using low-frequency DAS signals for tracking the wellbore thermal slugging. Jin et al. (2019) [15,16] demonstrate the use of a combination of DTS and DAS for production monitoring. The method is used for oil producers with low producing rates where the DTS was utilized to measure steady-state borehole temperatures, and the DAS is used to measure transient borehole flow velocities by tracking the signals generated by thermal slugging. This research was further developed for liquid and gas two-phase flow, specifically for gas–oil two-phase slug flow [17,18,19]. The experiments’ findings using a vertical flow loop indicated that the acoustic and thermal signals measured by DAS are sensitive to low air injection. Also, the thermal signatures obtained from the DAS data mainly respond to the water phase in the flow, while the air phase provides unique characteristics in the acoustic domain.
Another study shows that DAS measurements for two-phase gas–liquid slug flow are used to estimate the velocity and extract velocity patterns. The research includes two methods, k-f transform, and distributed cross-correlation. The results are validated by a commercial conductance-based phasemeter [20]. Many DFOS data processing are based on qualitative approaches only. In this article, experiments are conducted to develop a quantitative analysis for multiphase flow in pipes using DAS, which can expand the applications of DFOS data and contribute to better multiphase flow characterization.
Flow patterns in horizontal two-phase flows can range from dispersed, stratified, intermittent, and annular flow patterns depending on the rates of the phases, fluid properties, and pipe diameter [21]. Slug flow is a type of intermittent flow pattern, in which the gas and liquid phases flow alternatively. It is also one of the most common flow patterns encountered in the oil and gas transportation system [22]. This type of flow pattern have been extensively studied in the laboratory [23], such as [24,25,26,27], to name a few. The authors of [21,28] provide a detailed review of this flow pattern and the relevant modeling studies. A typical slug unit consists of a slug body region and a film region (also called gas pocket region), as depicted in Figure 1. Some previous studies claimed that the slug region could be considered as a dispersed bubble flow, and the film region as segregated flow [29]. The entire slug unit is moving forward at the translational velocity (vT), which is also the slug front velocity. Figure 1 illustrates the typical characteristic parameters for slug flow, including the translation velocity. From Figure 1, LU is the slug unit length, LS is the liquid slug body length, Lf is the length of the liquid film region, vs. is the average mixture velocity inside the slug body, vF is the average liquid velocity of the liquid film region, and vC is the velocity of the gas pocket in the liquid film region.
It is essential to characterize the slug flow, not only because of its common occurrence in the oil and gas transportation system, but also due to its alternating behavior of liquid slugs and gas pockets which could pose a risk to the production system and compromise the safety of the pipelines, especially to the joints and elbows [30,31,32]. By monitoring the slug flow characteristics, such as frequency, translational velocity, and the lengths of the slug body and the liquid film regions, we will be able to assess their risks to the pipeline systems and the piping components, react immediately if the slug flow behavior is determined to be damageable, and assess the multiphase flow rates for individual pipelines.
In this research, low-frequency strain-rate analysis of the DAS is used to detect gas–oil two-phase slug flow and quantify its characteristics. Novel algorithms were developed to characterize the slug translational velocity, and the lengths of slug body and the liquid film regions that had never been discussed in previous studies. These parameters are essential for slug-induced risk assessment and flow rate estimation. This study also employed a high-speed camera next to the fiber cables for validation purposes, achieving a good correlation among the measurements under all conditions tested.

2. Materials and Methods

The two-phase gas–oil experiments were conducted in a three-phase flow loop at the Colorado School of Mines, consisting of a transparent polyvinyl chloride (PVC) horizontal test section measuring 14 m long, a 5.25 cm inner diameter, and a 6.03 cm outer diameter. Figure 2 shows the schematic of the facility. Figure 3 is a picture of the test section, high-speed camera, and the pipe helically wrapped with yellow fiber cable. A fiber cable was wrapped helically on a 35 cm long section and was connected to a Terra15 Treble DAS interrogator. The winding pitch is 0.112 cm. The data were acquired with a sampling rate of 10 kHz and a special sampling of 0.816 m along the fiber cable. The gauge length that was applied is 2.4 m. The flow loop had a differential pressure transducer to measure the pressure drop in the horizontal test section. Pressure and temperature transducers were installed at the test section to monitor the pressure and temperature during the tests. The sensors were all connected to a data acquisition system and recorded at a frequency of 10 Hz.
Mineral oil and compressed air were used as testing fluids in this study. The mineral oil had a density of 810 kg/m3 and a viscosity of 0.012 Pa∙s at 19 °C and atmospheric pressure. While the air had a density of 1.2 kg/m3 and a viscosity of 1.83 × 10−5 Pa∙s at 19 °C and atmospheric pressure. The oil was introduced to the system through a wye inlet connector, and then the building air from a compressor was introduced afterward. Various rates for the gas and liquid phases were examined by varying the control valve for each phase or the pump motor through the variable frequency drive (VFD). The superficial oil velocity (defined as the volumetric flow rate divided by the pipe cross-sectional area) varied from 0.2 to 0.8 m/s, while the gas phase varied from 0.16 to 0.78 m/s. The test matrix presented in this work is shown in Table 1.
The facility was equipped with a high-speed camera (Phantom VEO640, Wayne, NJ, USA) to acquire a side view of the fluid flow behavior in the acrylic test section (Figure 3a). The video was recorded at 200 Hz with a resolution of 1024 × 700. A light source was fixed behind the pipe to provide illumination. The camera was used to record the flow pattern of testing fluids that were circulated through the pipe. Camera videos were also processed and analyzed to obtain the slug flow characteristic parameters and used to validate the ones obtained from DAS.

3. Data Processing Workflow and Results

In this section, the data processing is discussed in two sections based on measurement technology. The first part, which is also the main one, describes the automated workflow developed in this study to detect and track the slugging signals in the gas–oil two-phase slug flow and obtain the slug flow characteristic parameters from the DAS data. It also includes a discussion of the workflow process and its accuracy. The second part shows the high-speed camera data processing of the tested cases and the slug parameters obtained. The comparison of the results from DAS and the high-speed camera is discussed afterward, followed by a subsection showing the slug flow variations with the flowing conditions.

3.1. DAS Data Processing Workflow

This workflow uses low-frequency strain-rate analysis of the DAS signal to demonstrate its capabilities of quantifying two-phase slug flow characteristic parameters, including the slug frequency, slug body length, slug translational velocity, slug unit, and the liquid film region lengths. The DAS can capture these signals due to the sensitivity of helically wrapped fiber to the slight pipe diameter change that is caused by the pressure differential exerted by a passing slug.
The raw DAS data are read into the workflow for the selected beginning and ending times. To reduce noises and enhance the analysis, the median value is computed and removed from the channel direction, and then a low-pass filter is used to filter the data, displaying signals that are lower than 20 Hz. Figure 4a shows the low-frequency data for gas–oil two-phase slug flow. The horizontal axis is time, while the vertical axis shows the channel number, which corresponds to the fiber length from the interrogator. The wrapped test section is from around 45 m to 110 m. A zoomed portion of the data over a fifteen-second period in Figure 4b, presented four slugs passing which are indicated by the red and blue inclined lines. tSU denotes the time interval between the two adjacent slugs at a position in the test section around 80 m in fiber length from the interrogator (near the middle of the fiber cable wrapped section).
Figure 5 is a visualization of the DAS signal of a single slug passing through and the identification of the slug front and tail. One way to obtain the slug translational velocity, vTB, is to calculate the slope of the blue line, which indicates the slug front, given in Equation (1):
v T B = L 2 L 1 t 2 t 1   ,  
where (L2L1) represents the distance between two points shown in the DAS data;   is a correction factor that converts the fiber cable distance to the actual distance in the test section; (t2t1) represents the time required for this slug traveling from L1 to L2.
tSS in Figure 5 is the time interval between slug front and tail, t S l u g T a i l t S l u g F r o n t , representing the time required for the slug body to pass a particular point on the test section. It can be used to determine the slug body length, given in Equation (2):
L S = v T B t S S = v T B t S l u g T a i l t S l u g F r o n t   ,
The slug unit length, LU, can be determined using Equation (3), where tSU is the time interval between two adjacent slugs (Figure 4b):
L U = v T B t S U   ,
The length of the film region, LF, are then calculated using Equations (4):
L F = L U L S   ,
In this paper, we propose a more automatic workflow for the determination of the translational velocity and the lengths of the slug unit, slug body, and the film region. The following paragraphs describe this workflow.
We selected the data between channels 45 to 110 in the data processing to eliminate the excessive noise beyond this range. These channels are associated with the wrapped fiber section near the camera section. A sum is then taken in the channel direction after taking the absolute value of the signal and plotted against the time to identify the time of each slug, which are presented as the peaks in Figure 6. Each of the slugs was given a unique slug index number (ID#) to be used in comparison with the camera data at later stages.
The peaks in Figure 6 were counted to determine the slug frequency, fS, which is the number of slugs divided by the recording time. Note that the peaks that are very close are counted as one peak, which should correspond to a single slug but probably with a more complex structure. This prevents the overcounting of slugs due to adjacent peaks when calculating the slug frequency.
To determine the slug translational velocity, we used semblance, a quantitative measure of the waveform similarity that estimates the consistency of the waveforms from different channels. The process is explained in Figure 7, which includes four plots: (a) slug waterfall plot, (b) semblance vs. velocity, (c) waterfall plot after applying linear moveout using velocity that gives the highest semblance value from the previous plot, (d) vertical stacking in the distance direction to determine the width of the negative peak for slug body length calculation.
The plot generation processes are described as follows:
(a)
The waterfall, a plot of distance versus time, is generated for each slug.
(b)
A linear moveout correction (LMO) is applied, which is a velocity correction to shift traces in time based on an assumed velocity. The best velocity reflects the most consistency of the waveform among traces in time after the correction. The following equation represents the linear moveout correction:
t L M O n , v = d n v   ,
where n represents the trace number in the channel direction, d(n) is the fiber distance of the channel n , v represents the testing velocity, and t L M O represents the time shift applied to the channel.
(c)
The semblance is a quantitative measure of the waveform similarity from different channels, which is a metric commonly used in seismic processing. It is calculated using the following equations:
S = t = t 0 t t = t 0 + t n = 1 N f n t t L M O n , v 2 ,
where S is the summation of all waveforms in the channel direction, which is then squared and summed again in the time direction. fn(t) represents the data value of channel n at time t . N is the total number of channels, Δ t is the half window length in time for the semblance calculation, and t 0 is the center of the time window.
E = t = t 0 t t = t 0 + t n = 1 N f n 2 ( t t L M O n , v )   ,
where E(t0) is the summation of the square of all waveforms in both channel and time directions.
R v = S E N 1   E   ,
Finally, R   is the normalized semblance value of all the channels for the testing velocity v . An array of the different testing velocities versus semblance can be obtained, and the best velocity is chosen at the maximum semblance value which represents the highest consistency of waveforms after linear moveout correction. This can be used as the estimation of the slug translational velocity after applying a fiber-length-to-pipe-length ratio of 169:1.
(d)
All the signals are stacked vertically in the channel direction after the linear moveout correction using the best velocity, and the negative peak is identified with its start and end times, which represent the slug’s front and tail.
The slug characteristic parameters, including the slug body length, slug unit length, and the length of the film region, are calculated using Equations (2)–(4), respectively. An example of one of the slugs with its calculated parameters from the process in this section is shown in Table 2.
To increase the accuracy of the quantitative analysis on the slug characteristic parameters, a threshold for the calculated semblance after correction is selected at 0.1 to filter out some of the bad signals that may not represent an actual or complete slug signal and do not give a good quantitative analysis. An example is shown in Figure 8, where the automated workflow does not predict a reliable measurement for slug#12. This is also observed in Table 3 where the semblance after correction is very low (highlighted in red). In this example, the DAS signal in the selected window only contains a little part of a slug body shown at the bottom left in the waterfall plot (Figure 8a), resulting in unreliable quantification of the slug characteristics as anticipated. Nonetheless, these “bad” slugs represent a low percentage of the total population of the detected slugs that have reliable quantitative analysis and velocity measurement. Table 4 summarizes the percentage of the slugs with reliable quantitative analysis. The characteristic parameters of the slug were analyzed for these slugs and evaluated with data from the high-speed camera. The comparison is discussed in the later sub-section.
The summary of the workflow steps in the DAS data processing and the slug characteristic parameters that can be obtained is presented in Figure 9.

3.2. High-Speed Camera Data Processing

The camera data processing involves the recorded videos from the high-speed camera of the flow at each of the gas–oil two-phase flow cases presented in Table 1 previously. Each video is recorded after reaching the stabilization of the flow for about three minutes at a frequency of 200 Hz. Figure 10 presents a capture of a slug passing by for a liquid superficial velocity at 0.2 m/s and gas superficial velocity at 0.16 m/s. A scale was attached to the pipe which was used as a point of entry for the slug’s front and tail. The videos are then analyzed entirely and each slug that is passing would have a unique slug index and the following data processing for it:
(a)
The time when the slug front starts is determined once reaching the scale, tSF.
(b)
The time when the slug tail reaches the scale is documented, tST.
(c)
The time when the slug front reaches the beginning of the camera exposure is recorded, tSF_In.
(d)
The time when the slug front reaches the end of the camera exposure is recorded, tSF_Out.
(e)
Translational velocity, vT, is obtained by finding the length of the horizontal section that is exposed to the camera divided by the duration of exposure of each slug, i.e., vT = L/(tSF_Out − tSF_In).
(f)
Slug unit length is determined by LU = vT (tSF1 − tSF2), where (tSF1 − tSF2) is the time interval between two adjacent slugs.
(g)
The length of the slug body is determined by LS = vT (tSF − tST), and the length of the film region is determined using Equation 4 as presented previously.
(h)
Slug frequency, fS, is determined by counting the number of slugs divided by the corresponding recording time.
(i)
For each of the slug characteristic parameters obtained, the average and the median values were calculated over the full three-minute duration of the recorded video.

3.3. Data Validation

The accuracy of the DAS workflow process is validated using the high-speed camera videos which serve as the ground truth. An example of the comparison is presented in Figure 11, which shows the comparison of vT, LU, LS, and LF, between the data obtained from DAS (red points) and the high-speed camera videos (black points). The oil superficial velocity is 0.2 m/s, and the gas superficial velocity is 0.16 m/s. A good match between the DAS data and the high-speed camera data is observed.
One might have noticed that the slug characteristics are not the same at the same flowing condition. This fluctuation is the nature of slug flow. Dozens of previous studies have tried to develop statistical models using laboratory data to predict the slug characteristics, especially the slug length which is a critical parameter for facility design [33,34,35]. With the development of DAS, we can now monitor the slug characteristics directly in the field in real time.
To further evaluate the performance of DAS and the data processing workflow, we generated cross plots for all the gas–oil two-phase experimental data points, and the comparison of the median values is shown in Figure 12. Considering the median value of each slug parameter in the 13 experiments over the three-minute period of recorded data, a very good match can be observed with the camera data with an error bar of less than 20% in the majority of the cases. One astonishing observation is that the DAS can detect the number of slugs or quantify the slug frequency, very successfully for the conditions investigated in the current study. For the other parameters, the error seems to increase with the increase in gas flow rate or decrease in liquid flow rate. One of the reasons could be the reduced number of slugs investigated as the gas flow rate increases or the liquid flow rate decreases, so a smaller number of measurements are averaged. Another reason could be the change in the physical properties of the slug body, i.e., more gas is present in the slug body as the gas flow rate increases or liquid flow rate decreases, and the liquid slugs become more chaotic, frothy, and shorter, making the signal less clear and more difficult to be processed [31,32,36]. This phenomenon is illustrated in the images from the high-speed camera in Figure 13. On the other hand, we also assessed the measurement uncertainty from the high-speed camera. The errors for translational velocity, slug unit, and body lengths range from 0.87% to 1.1%. These errors are not visible in Figure 12.

3.4. Slugs at Different Flowing Conditions

In this subsection, we will show the slug flow behavior at different flowing conditions. Figure 14 shows the DAS signal after applying the low-pass filter and the median value removal in the channel direction, for three different liquid flow rates and four gas flow rates. The time duration is 3 s for all the figures. It can be clearly noticed that the slug frequency increases as the liquid flow rate increases. Moreover, the lines become steeper as the gas or liquid velocity increases, indicating an increase in the slug translational velocity. This observation is consistent with other previous studies for gas–liquid flow in horizontal pipes [37]. The slug body length cannot be directly read from the plots, as it is a function of time duration and also the slug translational velocity.
The pictures of some slugs in Figure 14 are shown in Figure 15 and Figure 16 for 0.5 m/s and 0.8 m/s liquid superficial velocities, respectively. As the gas velocity increases, more gas bubbles are entrained inside the slug body due to the high turbulence in the mixing zone at the slug front [38]. This increases the noise in the DAS signal, as depicted in Figure 14. We suspect that the intensity of the DAS signal within the slug body may correlate with the turbulence and entrained gas bubbles within the slug body. Further investigations are necessary to better understand this relationship quantitatively.

4. Discussion

Slug flow characterization is of great importance to many field applications. In the oil and gas industry, knowing the slug flow characteristics, such as the slug body length, translational velocity, and frequency, is crucial for facility design and operation. For example, the size and internal design of the separator should be able to handle the slugs to avoid flooding or ineffective separation. Slug catcher design also requires information on the slug flow characteristic parameters to achieve optimal gas–liquid separation [39]. Topside slugging control for a riser system requires real-time slug flow monitoring that is automated with the control valve systems [40]. Moreover, slug flow can accelerate corrosion and erosion due to their high translational velocities and the highly turbulent zoom at the slug front that could possibly lead to cavitation [41,42]. Knowing the slug flow characteristics, particularly in real time, empowers us to more confidently assess the risk of slug-flow-induced or -accelerated corrosion and erosion and adjust the operating parameters accordingly to minimize the risks. The current most widely used method in the field for multiphase surface facilities design and flow assurance is based on predictive models that are developed based on laboratory data or very limited field data. Their accuracy becomes questionable as it is scaled up to the field application. Because the slug flow behavior is transient and dramatically diverse in both the spatial and time domains, even at the same flowing conditions [22,43], real-time monitoring in the field will empower us to more effectively track their behaviors, allowing us to better optimize facility design and operation.
DAS is advanced in many aspects compared with other point sensors, such as γ-ray sensors, wire-mesh sensors, electrical impedance (capacitance and conductance) sensors, microwave sensors, optical sensors, etc., as previously mentioned in the introduction. It is completely non-intrusive and does not disturb the existing facility. It can provide measurements at multiple different locations and over a long distance in real time, instead of a point measurement like the other types of sensors, making it more cost-effective. Another big advantage is that DAS can work in harsh environments like downhole, where the other type of sensor can hardly survive.
In this paper, we introduce a new method for characterizing the slug translational velocity using semblance, showing great alignment with the observations from the high-speed camera that serve as the ground truth. Moreover, we demonstrated a new algorithm to characterize lengths of slug body, slug unit, and the liquid film region, which are the most critical parameters required in field design and operation as discussed above. The algorithm requires the fiber to be densely wrapped around the pipe, which increases the sensitivity of the strain rate measurement due to the pipe diameter variations caused by the slug-induced pressure fluctuations (see Figure 5).
The slug flow behavior should be related to different types of fluids and compositions, and the DAS signal should respond differently depending on the fluid properties and the slug flow behaviors. In the next phase, we are going analyze the DAS signals for other different fluid compositions, including gas–water, oil–water, and gas–oil–water flows.
Furthermore, we will investigate the causes of the close peaks in Figure 6. There is a possibility that these close peaks are induced by large gas pockets within long slug bodies, but we will further investigate this phenomenon from both fluid flow dynamics and data processing perspectives. Frequency domain dynamic averaging (FDDA) and/or activation function dynamic averaging (AFDA) methods will be tested.

5. Conclusions

A series of experiments to investigate the capability of DAS to quantify the slug characteristics were conducted in a horizontal pipe under different flowing conditions. A low pass filter and median value removal were applied to better extract the slug signals. We introduce a new method for characterizing the slug translational velocity using semblance, showing great alignment with the observations from the high-speed camera that serve as the ground truth. Moreover, we demonstrated a new algorithm to characterize the lengths of the slug body, slug unit, and the liquid film region. The algorithm requires the fiber to be densely wrapped around the pipe, which increases the sensitivity of the strain rate measurement due to the pipe diameter variations caused by the slug-induced pressure fluctuations.
Real-time slug flow monitoring and characterization are crucial for many industries that involve slug slow, such as petroleum engineering, geothermal wells that have two-phase flow production, nuclear engineering, etc. This research demonstrates the capabilities of DAS to quantitatively measure the slug flow characteristic parameters, especially the lengths of the slug body, slug unit, and the liquid film region, which are the most critical parameters required in field facility design and operation. This technology will offer the industry a more optimal solution for facility design and ensure safer operational practices.

Author Contributions

Conceptualization, S.A., G.J. and Y.F.; methodology, S.A., G.J. and Y.F.; software, S.A., G.J. and Y.F.; validation, S.A., G.J. and Y.F.; formal analysis, S.A., G.J. and Y.F.; investigation, S.A., G.J. and Y.F.; resources, G.J. and Y.F.; data curation, S.A.; writing—original draft preparation, S.A.; writing—review and editing, G.J. and Y.F.; visualization, G.J. and Y.F.; supervision, G.J. and 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.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data may be provided upon request.

Acknowledgments

The authors acknowledge Ge Jin’s and Yilin Fan’s research groups’ support throughout this work. We also acknowledge Denghong Zhou, Kanat Karatayev, and Mouna-Keltoum Benabid for their significant help with running the experiments. We appreciate the scholarship provided by the Kuwait Petroleum Corporation (KPC), which supported Sharifah Ali’s Ph.D. study at the Colorado School of Mines.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Stalford, H.; Ahmed, R. Intelligent Casing-Intelligent Formation (ICIF) Design. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 5–8 May 2014; p. D011S012R008. [Google Scholar]
  2. Zhong, Z.Y.; Zhi, X.L.; Yi, W.J. Oil Well Real-time Monitoring with Downhole Permanent FBG Sensor Network. In Proceedings of the 2007 IEEE International Conference on Control and Automation, Guangzhou, China, 30 May–1 June 2007; pp. 2591–2594. [Google Scholar]
  3. Falcone, G.; Hewitt, G.; Alimonti, C. Multiphase Flow Metering: Principles and Applications; Elsevier: Amsterdam, The Netherlands, 2009. [Google Scholar]
  4. Mukhopadhyay, S.C. New Developments in Sensing Technology for Structural Health Monitoring; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  5. Hill, D.; Fellow, S. Managing Oil and Gas with Fibre-Optic Sensing; EIC Energy Focus: London, UK, 2014. [Google Scholar]
  6. Kinet, D.; Mégret, P.; Goossen, K.W.; Qiu, L.; Heider, D.; Caucheteur, C. Fiber Bragg grating sensors toward structural health monitoring in composite materials: Challenges and solutions. Sensors 2014, 14, 7394–7419. [Google Scholar] [CrossRef]
  7. Hartog, A.H. An Introduction to Distributed Optical Fibre Sensors; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  8. Boering, M.; Braal, R.; Cheng, L.K. Toward the next fiber optic revolution and decision making in the oil and gas industry. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 6–9 May 2013. [Google Scholar] [CrossRef]
  9. In’t Panhuis, P.; den Boer, H.; van der Horst, J.; Paleja, R.; Randell, D.; Joinson, D.; McIvor, P.; Green, K.; Bartlett, R. Flow monitoring and production profiling using DAS. In SPE Annual Technical Conference and Exhibition? SPE: Kuala Lumpur, Malaysia, 2014. [Google Scholar] [CrossRef]
  10. Raab, T.; Reinsch, T.; Aldaz Cifuentes, S.R.; Henninges, J. Real-time well-integrity monitoring using fiber-optic distributed acoustic sensing. SPE J. 2019, 24, 1997–2009. [Google Scholar] [CrossRef]
  11. Boone, K.; Ridge, A.; Crickmore, R.; Onen, D. Detecting leaks in abandoned gas wells with fibre-optic distributed acoustic sensing. In Proceedings of the IPTC 2014: International Petroleum Technology Conference, Doha, Qatar, 19–22 January 2014; European Association of Geoscientists & Engineers: Rotterdam, The Netherlands, 2014. [Google Scholar] [CrossRef]
  12. Jin, G.; Roy, B. Hydraulic-fracture geometry characterization using low-frequency DAS signal. Lead. Edge 2017, 36, 975–980. [Google Scholar] [CrossRef]
  13. Mateeva, A.; Lopez, J.; Mestayer, J.; Wills, P.; Cox, B.; Kiyashchenko, D.; Yang, Z.; Berlang, W.; Detomo, R.; Grandi, S. Distributed acoustic sensing for reservoir monitoring with, V.S.P. Lead. Edge 2013, 32, 1278–1283. [Google Scholar] [CrossRef]
  14. van der Horst, J.; den Boer, H.; in’t Panhuis, P.; Kusters, R.; Roy, D.; Ridge, A.; Godfrey, A. Fibre Optic Sensing for Improved Wellbore Surveillance. In Proceedings of the International Petroleum Technology Conference, Beijing, China, 26–28 March 2013. [Google Scholar] [CrossRef]
  15. Jin, G.; Friehauf, K.; Roy, B. The calibration of double-ended distributed temperature sensing for production logging purposes. In Proceedings of the SPE/AAPG/SEG Unconventional Resources Technology Conference, Denver, CO, USA, 22–24 July 2019. [Google Scholar] [CrossRef]
  16. Jin, G.; Friehauf, K.; Roy, B.; Constantine, J.J.; Swan, H.W.; Krueger, K.R.; Raterman, K.T. Fiber optic sensing-based production logging methods for low-rate oil producers. In Proceedings of the SPE/AAPG/SEG Unconventional Resources Technology Conference, Denver, CO, USA, 22–24 July 2019. [Google Scholar] [CrossRef]
  17. Titov, A.; Fan, Y.; Kutun, K.; Jin, G. Distributed Acoustic Sensing (DAS) Response of Rising Taylor Bubbles in Slug Flow. Sensors 2022, 22, 1266. [Google Scholar] [CrossRef]
  18. Titov, A.; Jin, G.; Fan, Y.; Tura, A.; Kutun, K.; Miskimins, J. Distributed fiber-optic sensing based production logging investigation: Flowloop experiments. First EAGE Workshop on Fibre Optic Sensing. Eur. Assoc. Geosci. Eng. 2020, 2020, 1–5. [Google Scholar] [CrossRef]
  19. Titov, A.; Fan, Y.; Jin, G.; Tura, A.; Kutun, K.; Miskimins, J. Experimental investigation of distributed acoustic fiber-Optic Sensing in production logging: Thermal slug tracking and multiphase flow characterization. In SPE Annual Technical Conference and Exhibition? Day 4 Thu, October 29, 2020; SPE: Kuala Lumpur, Malaysia, 2020. [Google Scholar] [CrossRef]
  20. Weber, G.H.; dos Santos, E.N.; Gomes, D.F.; Santana, A.L.B.; da Silva, J.C.C.; Martelli, C.; Pipa, D.R.; Morales, R.E.; de Camargo Júnior, S.T.; da Silva Junior, M.F.; et al. Measurement of Gas-Phase Velocities in Two-Phase Flow Using Distributed Acoustic Sensing. IEEE Sens. J. 2023, 23, 3597–3608. [Google Scholar] [CrossRef]
  21. Al-Safran, E.M.; Brill, J.P. Applied Multiphase Flow in Pipes and Flow Assurance: Oil and Gas Production; Society of Petroleum Engineers: Kuala Lumpur, Malaysia, 2017. [Google Scholar]
  22. Al-Kayiem, H.H.; Mohmmed, A.O.; Al-Hashimy, Z.I.; Time, R.W. Statistical assessment of experimental observation on the slug body length and slug translational velocity in a horizontal pipe. Int. J. Heat Mass Transf. 2017, 105, 252–260. [Google Scholar] [CrossRef]
  23. Raimondi, L. Gas/liquid two-phase flow in pipes: Slugs, classical flow-map, and 1D compositional simulation. SPE J. 2022, 27, 532–551. [Google Scholar] [CrossRef]
  24. Ekinci, S. Pipe Inclination Effects on Slug Flow Characteristics of High Viscosity Oil-Gas Two-Phase Flow. Master’s Thesis, The University of Tulsa, Tulsa, OK, USA, 2015. [Google Scholar]
  25. 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]
  26. Liu, Y.; Upchurch, E.R.; Ozbayoglu, E.M. Experimental Study of Single Taylor Bubble Rising in Stagnant and Downward Flowing Non-Newtonian Fluids in Inclined Pipes. Energies 2021, 14, 578. [Google Scholar] [CrossRef]
  27. 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]
  28. Shoham, O. Mechanistic Modeling of Gas-Liquid Two-Phase Flow in Pipes; Society of Petroleum Engineers: Kuala Lumpur, Malaysia, 2006. [Google Scholar]
  29. 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 Res. Technol. 2003, 125, 266–273. [Google Scholar] [CrossRef]
  30. Ludwig, E.E. Applied Process Design for Chemical and Petrochemical Plants: Volume 1; Elsevier Science: Amsterdam, The Netherlands, 1995. [Google Scholar]
  31. Fan, Y.; Soedarmo, A.; Pereyra, E.; Sarica, C. A comprehensive review of pseudo-slug flow. J. Pet. Sci. Eng. 2022, 217, 110879. [Google Scholar] [CrossRef]
  32. Fan, Y.; Pereyra, E.; Sarica, C. Experimental study of pseudo-slug flow in upward inclined pipes. J. Nat. Gas Sci. Eng. 2020, 75, 103147. [Google Scholar] [CrossRef]
  33. Bernicot, M.F.; Drouffe, J.-M. A Slug-Length Distribution Law for Multiphase Transportation Systems. SPE Prod. Eng. 1991, 6, 166–170. [Google Scholar] [CrossRef]
  34. Barnea, D.; Taitel, Y. A model for slug length distribution in gas-liquid slug flow. Int. J. Multiph. Flow 1993, 19, 829–838. [Google Scholar] [CrossRef]
  35. Al-Safran, E.M.; Sarica, C.; Zhang, H.Q.; Brill, J.P. Probabilistic/mechanistic modeling of slug-length distribution in a horizontal pipeline. SPE Prod. Facil. 2005, 20, 160–172. [Google Scholar] [CrossRef]
  36. Fan, Y.; Aljasser, M. Unified modeling of pseudo-slug and churn flows for liquid holdup and pressure gradient predictions in pipelines and wellbores. SPE Prod. Oper. 2022, 38, 350–366. [Google Scholar] [CrossRef]
  37. Abdulkadir, M.; Hernandez-Perez, V.; Lowndes, I.S.; Azzopardi, B.J.; Sam-Mbomah, E. Experimental study of the hydrodynamic behaviour of slug flow in a horizontal pipe. Chem. Eng. Sci. 2016, 156, 147–161. [Google Scholar] [CrossRef]
  38. Al-Safran, E.; Kora, C.; Sarica, C. Prediction of slug liquid holdup in high viscosity liquid and gas two-phase flow in horizontal pipes. J. Pet. Sci. Eng. 2015, 133, 566–575. [Google Scholar] [CrossRef]
  39. Suleimanov, R.I.; Ya Khabibullin, M. Optimization of the design of the scrubber separator slug catcher. IOP Conf. Ser. Mater. Sci. Eng. 2020, 905, 012090. [Google Scholar] [CrossRef]
  40. Nnabuife, S.G.; Tandoh, H.; Whidborne, J.F. Slug flow control using topside measurements: A review. Chem. Eng. J. Adv. 2022, 9, 100204. [Google Scholar] [CrossRef]
  41. Sun, J.; Jepson, W.P. Slug Flow Characteristics and Their Effect on Corrosion Rates in Horizontal Oil and Gas Pipelines; SPE: Kuala Lumpur, Malaysia, 1992. [Google Scholar] [CrossRef]
  42. Sarac, B.E.; Stephens, D.S.; Eisener, J.; Rosselló, J.M.; Mettin, R. Cavitation bubble dynamics and sonochemiluminescence activity inside sonicated submerged flow tubes. Chem. Eng. Process.—Process Intensif. 2020, 150, 107872. [Google Scholar] [CrossRef]
  43. Olbrich, M.; Bär, M.; Oberleithner, K.; Schmelter, S. Statistical characterization of horizontal slug flow using snapshot proper orthogonal decomposition. Int. J. Multiph. Flow 2021, 134, 103453. [Google Scholar] [CrossRef]
Figure 1. Physical structure of slug flow with its characteristic parameters.
Figure 1. Physical structure of slug flow with its characteristic parameters.
Sensors 24 03402 g001
Figure 2. The schematic of the three-phase multiphase flow loop.
Figure 2. The schematic of the three-phase multiphase flow loop.
Sensors 24 03402 g002
Figure 3. (a) Pictures of the testing section with the helically wrapped yellow fiber cables, and the high-speed camera taking videos. (b) A closer look at the wrapped fiber cable.
Figure 3. (a) Pictures of the testing section with the helically wrapped yellow fiber cables, and the high-speed camera taking videos. (b) A closer look at the wrapped fiber cable.
Sensors 24 03402 g003
Figure 4. (a) Low-frequency DAS signal for slug flow; (b) fifteen-second time selected to illustrate the passing slugs in the signal.
Figure 4. (a) Low-frequency DAS signal for slug flow; (b) fifteen-second time selected to illustrate the passing slugs in the signal.
Sensors 24 03402 g004
Figure 5. DAS signal of a single slug in the gas–oil two-phase slug flow.
Figure 5. DAS signal of a single slug in the gas–oil two-phase slug flow.
Sensors 24 03402 g005
Figure 6. The peak value of the slugs over a three-minute recorded time for the case of 0.2 m/s oil velocity and 0.16 m/s gas velocity.
Figure 6. The peak value of the slugs over a three-minute recorded time for the case of 0.2 m/s oil velocity and 0.16 m/s gas velocity.
Sensors 24 03402 g006
Figure 7. An example of quality check plots for workflow evaluation at slug#45. (a) Slug waterfall plot; (b) semblance as a function of testing velocity; (c) waterfall plot after applying linear moveout using velocity that gives the highest semblance value in (b); (d) vertical stacking in the distance direction.
Figure 7. An example of quality check plots for workflow evaluation at slug#45. (a) Slug waterfall plot; (b) semblance as a function of testing velocity; (c) waterfall plot after applying linear moveout using velocity that gives the highest semblance value in (b); (d) vertical stacking in the distance direction.
Sensors 24 03402 g007
Figure 8. An example of a bad slug (ID #12) that does not give good quantitative analysis. (a) Slug waterfall plot; (b) semblance as a function of testing velocity; (c) waterfall plot after applying linear moveout using velocity that gives the highest semblance value in (b); (d) vertical stacking in the distance direction.
Figure 8. An example of a bad slug (ID #12) that does not give good quantitative analysis. (a) Slug waterfall plot; (b) semblance as a function of testing velocity; (c) waterfall plot after applying linear moveout using velocity that gives the highest semblance value in (b); (d) vertical stacking in the distance direction.
Sensors 24 03402 g008
Figure 9. Workflow of the DAS data processing for slug flow characterization.
Figure 9. Workflow of the DAS data processing for slug flow characterization.
Sensors 24 03402 g009
Figure 10. Slug flow captured from a recorded video by a high-speed camera at 0.16 m/s gas velocity and 0.2 m/s oil velocity.
Figure 10. Slug flow captured from a recorded video by a high-speed camera at 0.16 m/s gas velocity and 0.2 m/s oil velocity.
Sensors 24 03402 g010
Figure 11. Comparison between the slug characteristic parameters from the DAS data and the camera for the case of 0.2 m/s oil superficial velocity and 0.16 m/s gas superficial velocity.
Figure 11. Comparison between the slug characteristic parameters from the DAS data and the camera for the case of 0.2 m/s oil superficial velocity and 0.16 m/s gas superficial velocity.
Sensors 24 03402 g011
Figure 12. Comparison between the median value from DAS and high-speed camera.
Figure 12. Comparison between the median value from DAS and high-speed camera.
Sensors 24 03402 g012
Figure 13. Flow pattern captured from the high-speed camera at 0.2 m/s oil superficial velocity and different gas flow rates.
Figure 13. Flow pattern captured from the high-speed camera at 0.2 m/s oil superficial velocity and different gas flow rates.
Sensors 24 03402 g013
Figure 14. DAS signals at different flowing conditions (vSL: liquid superficial velocity, vSG: gas superficial velocity. Superficial velocity is defined as the volumetric flow rate divided by the pipe cross-sectional area).
Figure 14. DAS signals at different flowing conditions (vSL: liquid superficial velocity, vSG: gas superficial velocity. Superficial velocity is defined as the volumetric flow rate divided by the pipe cross-sectional area).
Sensors 24 03402 g014
Figure 15. Pictures from high-speed camera videos for the slugs shown in Figure 14 for vSL = 0.5 m/s.
Figure 15. Pictures from high-speed camera videos for the slugs shown in Figure 14 for vSL = 0.5 m/s.
Sensors 24 03402 g015
Figure 16. Pictures from high-speed camera videos for the slugs shown in Figure 14 for vSL = 0.8 m/s.
Figure 16. Pictures from high-speed camera videos for the slugs shown in Figure 14 for vSL = 0.8 m/s.
Sensors 24 03402 g016
Table 1. Test matrix in the current study.
Table 1. Test matrix in the current study.
Case# 1Superficial Oil Velocity [m/s]Superficial Gas Velocity [m/s]
1–50.20.16, 0.31, 0.47, 0.62, 0.78
6–90.50.16, 0.31, 0.47, 0.62
10–130.80.16, 0.31, 0.47, 0.62
1 The case number # in each row corresponds to different superficial gas velocity. For example, Case#1 corresponds to an 0.2 m/s superficial oil velocity and 0.16 m/s superficial gas velocity; Case#3 corresponds to an 0.2 m/s superficial oil velocity and 0.47 m/s superficial gas velocity.
Table 2. Characteristic parameters for slug#45 determined from the DAS signals using the proposed approach.
Table 2. Characteristic parameters for slug#45 determined from the DAS signals using the proposed approach.
Parameters ValueParameters Value
Slug#ID45tSlugFront (s)0.406
Slug occurrence time (s)170.26tTail (s)0.668
Semblance before correction−0.008588tSS (s)0.262
Best velocity (m/s)265.52Slug Frequency (s−1)0.239
Semblance after correction0.602624Slug Body Length, LS (m) 0.410
Slug translational velocity, vT (m/s)1.565Liquid Film Region Length, LF (m)6.125
tSU (s)4.176Negative Peak Time (s)0.500
Slug Unit Length, LU (m)6.536Negative Peak Value−3.487 × 10−6
Table 3. Characteristic parameters for slug#12 determined from the DAS signals using the proposed approach that does not give good quantitative analysis.
Table 3. Characteristic parameters for slug#12 determined from the DAS signals using the proposed approach that does not give good quantitative analysis.
Parameters ValueParametersValue
Slug#ID412tSlugFront (s)0.000
Slug occurrence time (s)43.074tTail (s)0.673
Semblance before correction−0.004933tSS (s)0.673
Best velocity (m/s)279.31Slug Frequency (s−1)1.429
Semblance after correction0.024887Slug Body Length, LS (m) 1.108
Slug translational velocity, vT (m/s)1.6464Liquid Film Region Length, LF (m)0.046
tSU (s)0.7Negative Peak Time (s)0.045
Slug Unit Length, LU (m)1.154Negative Peak Value0.363 × 10−6
Table 4. Percentage of the slugs with reliable quantitative analysis from the DAS automated algorithm.
Table 4. Percentage of the slugs with reliable quantitative analysis from the DAS automated algorithm.
Case #vSO [m/s]vSG [m/s]P * [%]
10.20.1678
20.3184
30.4788
40.6273
50.3563
60.50.1681
70.3190
80.4788
90.6294
100.80.1694
110.3194
120.4795
130.6298
* P: Percentage of slugs with reliable quantitative analysis and used for slug characteristic parameter determination.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ali, S.; Jin, G.; Fan, Y. Characterization of Gas–Liquid Two-Phase Slug Flow Using Distributed Acoustic Sensing in Horizontal Pipes. Sensors 2024, 24, 3402. https://doi.org/10.3390/s24113402

AMA Style

Ali S, Jin G, Fan Y. Characterization of Gas–Liquid Two-Phase Slug Flow Using Distributed Acoustic Sensing in Horizontal Pipes. Sensors. 2024; 24(11):3402. https://doi.org/10.3390/s24113402

Chicago/Turabian Style

Ali, Sharifah, Ge Jin, and Yilin Fan. 2024. "Characterization of Gas–Liquid Two-Phase Slug Flow Using Distributed Acoustic Sensing in Horizontal Pipes" Sensors 24, no. 11: 3402. https://doi.org/10.3390/s24113402

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop