Next Article in Journal
Comparison of Emissions and Efficiency of Two Types of Burners When Burning Wood Pellets from Different Suppliers
Next Article in Special Issue
Fast Well Control Optimization with Two-Stage Proxy Modeling
Previous Article in Journal
A MIMO–ANFIS-Controlled Solar-Fuel-Cell-Based Switched Capacitor Z-Source Converter for an Off-Board EV Charger
Previous Article in Special Issue
Application of Digitalization in Real-Time Analysis of Drilling Dynamics Using Along-String Measurement (ASM) Data along Wired Pipes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Review of Different Methods for Identification of Transients in Pressure Measurements by Permanent Downhole Gauges Installed in Wells

1
Faculty of Science and Technology, Department of Electrical Engineering and Computer Science, University of Stavanger, 4021 Stavanger, Norway
2
NORCE, Norwegian Research Center, 4021 Stavanger, Norway
3
Institute of GeoEnergy Engineerin, School of Energy, Geoscience, Infrastructure and Society, Heriot-Watt University, Edinburgh EH14 4AS, UK
*
Author to whom correspondence should be addressed.
Submission received: 25 December 2022 / Revised: 10 January 2023 / Accepted: 18 January 2023 / Published: 8 February 2023

Abstract

:
Permanent downhole gauges (PDG) are massively installed in injection and production wells operated in different industries such as oil and gas, geological CO2 storage, and the geothermal industry. These gauges provide a vast amount of real-time pressure measurements. The pressure measurements may be divided into periods with a predominantly monotonic change of pressure in response to a sudden change of rate, called transients. These transients are caused by well operations, such as variation of injection or production rate and well shut-ins. Transient identification is one of the important steps in processing and interpreting the PDG data. Traditional transient identification is performed by processing and analyzing with human involvement, which is a step in post-operation well analysis. In modern well surveillance technology, permanent and reliable data transmission from the wellbore to the surface provide the possibility to analyze well performance in real time or proactively. So automated transient identification is a practical demand, but a challenge at the same time. This article starts with the definition of a transient, then reviews and compares seven methods for transient identification proposed by previous works available in the literature. A comparative analysis of these methods is carried out accounting for the detection algorithm and procedure, results of testing, and general positive and negative sides of performance and application of these methods. The results of this review facilitate further developments of field data interpretation techniques by the R&D community and academia and may help in the selection of a proper method for further application in well surveillance workflows developed in the industry.

1. Introduction

Permanent downhole gauges (PDG) are reliable well data sources in modern well surveillance systems, which are installed in wellbores near the reservoir sand-face, where the wells are connected to the reservoir via open-hole or perforation intervals, Figure 1 shows a possible position of the gauges in the wellbores. PDG provides a long-term continuous record of data, such as pressure and temperature, which can be recorded for a well’s life during production or injection, reflecting all changes in the well rate and well shut-ins. The PDG data may be used for many applications, such as detection of changes in reservoir properties [1,2], evaluation of production or injection capacity [3], effect of the well interventions [4], monitoring of well interference [5,6], and monitoring of changes in production capacity [7,8]. So, PDGs become a significant source of downhole data for pressure transient analysis in oil and gas production [9] and geological CO2 storage [10]. Athichanagorn suggested a multi-step workflow for processing and interpretation of the PDG data, which includes: (i) outlier removal, (ii) denoising, (iii) transient identification, (iv) data reduction, (v) flow-history reconstruction, (vi) behavioral filtering, and (vii) moving window analysis [11]. His workflow is a clear example to illustrate the procedure of performing reservoir analysis by using PDG.
PDG data are used to address some challenges in reservoir development. First, PDG provides nearly real-time data directly from the well site to reservoir and production engineer. Real-time data interpretation is a valuable resource for decision making on well operations and reservoir management. However, these vast volumes of data containing outliers and noise usually take much human effort to process, making extraction of the reservoir information difficult or provide the information with a delay. Second, to utilize the PDG dataset, it is essential to divide the huge volume of long-term PDG data into individual transients, such as flowing transients or shut-in transients; then, the data interpretation can be carried out for the identified transients. Automation of transient identification should help to facilitate analysis of big well surveillance data sets available in the industry, giving also a chance for on-the-fly analysis.
Reservoir studies, such as pressure transient analysis (PTA), were mainly carried out based on a shut-in data due to low quality and frequency of rate measurements during flowing periods. Installation of PDGs in combination with frequent and accurate flow-metering provided basis for inclusion of flowing periods in PTA applications [1]. It is a challenge to identify the desired transients from the original PDG data because (1) the volume of PDG data is huge, (2) the measurements usually contain outliers and noise, and (3) the rate data from the surface are not synchronized with downhole pressure data. The accuracy of transient identification may affect interpretation of reservoir properties. For example, inaccurate identification may lead to biased estimates of well and reservoir properties like well skin and reservoir permeability, due to inaccurate pressure drop and early time response of the transient in focus. In application of the deconvolution method, inaccurate transient identification may also mislead the interpretation of the reservoir behavior [12,13]. A few works have been published proposing solutions for transient identification issue.
In this article, we review seven methods for transients identification from the literature, discuss, and compare their performance. These methods use the pressure data from the PDG as the only data to identify transients because the flow rate is not always measured in line with pressure or may have a much lower sampling rate. Even though there are some data analysis and machine learning applications based on PDG data published, for example in [14,15,16,17], these studies were usually not concentrated on transient identification problem itself. For a clear understanding of the characteristics of each reviewed method, the terminology is first summarized using previous statements found in the literature and the terms of a ’transient’ and a ’break point’ are introduced. Then seven methods are categorized into four types, which are (1) wavelet decomposition methods, (2) pressure data trend methods, (3) pattern recognition methods, and (4) a segmentation method. A comparison of these methods is given to demonstrate the features of each method, considering the detection theory and procedure, results from the literature, and general positive and negative sides of the performance and application of these methods.
Three main objectives of reviewing these methods may be specified as:
  • Define the common terms used in processing pressure measurements from permanent downhole gauges (PDG);
  • Provide strategies to identify transients from PDG datasets for well surveillance or similar transient datasets;
  • Facilitate further developments of field data interpretation techniques by the R&D community and academia via helping in the selection of a proper method for further application in well surveillance workflows applied in industry.
This review is organized as below:
  • Introduce the terminology for transient identification;
  • Introduce the reviewed transient identification methods;
  • Compare the methods;
  • Discussions;
  • Conclusions.

Terminology

The term transient has a specific meaning when used in different fields. In terms of PDG data processing and interpretation, transients can be defined as pressure measurements’ responses to flow rate changes. An abrupt change of measured pressure usually occurs due to a sudden change of flow rate; then, distinct flow periods may be separated by detection of these abrupt changes. In this context, definitions of a ’transient’, a ’break point’, and transient identification as a general task proposed by different researchers are reviewed and summarized hereafter. Regarding the description of break point, Khong [18] regarded break points as sudden pressure changes in pressure signals caused by the flow rate change. Rai [19] described a break point as a high or low point in a detailed signal. Nomura [12] referred to break point as the starting point of each transient. Liu [20] described a break point as a point where a flow rate change event happened, and it usually indicates the end of the previous transient and the beginning of the next transient. Regarding the description of transients, Thomas [21] regarded transients as distinct flow periods when the pressure changes occur with flow rate changes and break points marked as the start and end of a single transient. Regarding the description of transient identification, Athichanagorn [22] referred to transient identification as a procedure to detect changes in pressure signal when the flow rate alters. Li [23] treated transient identification as a methodology aiming at automatically identifying the time at which the pressure trend changes because of a rate change. The critical issue in transient identification according to his study is to select the proper break point of every flow period. Liu [20] mentioned that transient identification requires break-point detection. Tian [24] referred to the process of detecting the points to separate multiple transients as break-point identification. We, therefore, tried to define these common terms by summarizing the concepts from the above statements in this literature review and putting the methods for transient identification into the same context. As a result, the following formulations for these terms are suggested:
  • A pressure transient is a predominantly monotonic change of pressure in response to a sudden change of rate;
  • A break point is a point in time separating two transients;
  • Transient identification is a process of dividing pressure time series into sequential transients based on the objectives of the data interpretation.
The first and second definitions above are in line with those suggested in the literature reviewed but slightly modified to account for different transients, usually observed in real data, namely flowing and shut-in periods. The definition of transient identification task includes an extension compared to common definitions found in the literature highlighting the fact that the same sequence of multiple rate transients may be considered as different transients or as a single transient following objective of the study. As an example, step rate testing generates sequences of transients, where each transient in the sequence may be separately interpreted [25], while the transient response to all the sequence of rate steps may also be considered as a single transient.

2. Methods for Transient Identification

The seven reviewed methods for transient identification are:
  • Wavelet modulus maxima method;
  • Data grouping method;
  • Pressure-data slope method;
  • Pressure derivative method;
  • Filter convolution method;
  • Image pattern method;
  • Segmentation method.
The detection of breakpoints is a crucial task in identifying transients. The number of breakpoints determines the number of transients present, and their location determines the duration of each transient. As flow rate data is not always available, the methods discussed in this article for identifying transients rely solely on pressure data from PDG. These methods are categorized into four types: wavelet decomposition, pressure data trend analysis, pattern recognition, and segmentation. Pre-processing and thresholding are common techniques used across all methods, but they require human intervention for optimal results. However, pre-processing PDG data under the assumption that noise and outliers are solely caused by noise and do not reflect the behavior of the reservoir may result in the loss of valuable information. Additionally, thresholding is often used to decide the number of transients, which can lead to over- or under-detection.

2.1. Wavelet Decomposition Methods

Wavelet transformation has been widely studied and applied in signal processing, with many researchers proposing and improving wavelet decomposition methods for PDG signal processing [18,22,26,27,28,29]. The wavelet transform is particularly useful for PDG data pre-processing and the decomposed detail signals can be used for subsequent detection of transients. Wavelet-based methods can be understood as the original data being processed through a low-pass filter and a high-pass filter. After filtering, the original data is decomposed into two signals—the approximated data from the low-pass filter and the detailed data from the high-pass filter. Breakpoints are detected by setting a criterion over a certain threshold in the detailed data, as shown in Figure 2. The number of identified transients depends on the threshold selection. Khong [18] has used Fourier analysis to filter out false breakpoints based on wavelet transform. Ouyang [28] further studied the mechanism behind transient detection by considering reservoir parameters.

2.1.1. Wavelet Modulus Maxima Method

a.
Detection theory
Athichanagorn processed the original data by using spline wavelet decomposition, which decomposes the data into approximation signal data by using the scaling function and detail signal data by using the spline wavelet function [22]. Depending on the resolution level, the scaling function provides an approximation of the original data, while the wavelet function provides the differences between the approximation data and the original data. Breakpoints are detected by finding the modulus maxima over a certain pressure slope threshold, denoted as τ , in the detailed signal at an intermediate decomposition level L.
b.
Detection procedure
The detection procedure can be defined as follows:
  • Pre-process data for outlier removal and denoising;
  • Determine the level of spline wavelet decomposition L, according to desired data resolution;
  • Try proper slope or pressure thresholds to detect the break points by wavelet modulus maxima higher than threshold τ ;
  • Confirm the position of break points by checking slope changes.
c.
Experiment result
The author conducted three experiments using this method, and Figure 3 shows the detection results based on the denoised data. The vertical lines indicate the positions of the detected breakpoints. The level of spline wavelet decomposition used was intermediate, and the pressure slope threshold, τ , was set at 10 psi/h.
d.
Advantages and Disadvantages
This method can effectively detect obvious breakpoints, but it may fail to detect breakpoints with small magnitudes. The final detection quantity of breakpoints is determined by the threshold selection. Pre-processing, including outlier removal and denoising, is necessary to avoid false breakpoint detection. The detection result is highly affected by the wavelet decomposition level and the threshold, which requires human intervention to achieve the desired result. Additionally, the final position of a breakpoint needs to be confirmed in the original data.

2.1.2. Data Grouping Method

a.
Detection theory
Li [23] presented that decomposition by Haar wavelet transforms, a discrete wavelet transform, can detect more accurate breakpoint positions compared to spline wavelet transform (used in the wavelet modulus maxima method). Based on Haar wavelet decomposition, Li proposed a method to detect breakpoints by grouping data in the detailed signal. The threshold τ is decided by statistical grouping with the assistance of the mean μ and standard deviation δ calculated from the detailed signal.
b.
Detection procedure
The detection procedure can be detailed as follows:
  • Pre-process includes removing overlap data, outliers, and denoising;
  • Apply wavelet decomposition and use the detail signal in Level 1 as default data resolution for subsequent processing;
  • Identify the break points by grouping the data by setting proper threshold τ and regarding the break points located in the data group in which the data values are bigger than the threshold;
  • Refine break points position.
c.
Experiment result
The author conducted five experiments using this method [23], and the significant transients were detected, but with over-detection on transients with small magnitudes. The periods with different color represent different flow periods, and the period with green color represents the final flow period (Figure 4).
d.
Advantages and Disadvantages
Li improved detection of the break points with a more precise position in the original data by using the Haar wavelet transform and detecting more break points with smaller magnitudes by determining the threshold in statistic consideration. However, pre-processing including outlier removal and denoising is necessary, and the threshold selection by a statistic method results in over-detection with redundant transients.

2.2. Pressure Data Trend Methods

Compared to wavelet decomposition methods that use detailed signal data to detect breakpoints, pressure data trend methods use the detailed data as the source to detect breakpoints by calculating the pressure slope [30,31] or pressure derivative [19] over a certain threshold, as shown in Figure 5. These methods are straightforward to implement and can detect significant breakpoints when the original PDG data are properly pre-processed. However, the desired threshold is selected through trial and error and the detected position still needs to be confirmed in the original data.

2.2.1. Pressure–Data Slope Method

a.
Detection theory
Viberti proposed a method that uses the slope of pressure data as the criterion to identify transients [30]. The original PDG data are first pre-processed by wavelet algorithms for outlier removal and denoising. Then, the detection process starts by setting the number of points k for calculating the slope in the pre-processed data. Finally, the calculated slope values greater or less than the thresholds τ 1 , τ 2 will be detected as breakpoints.
b.
Detection procedure
The detection procedure can be detailed as follows:
  • Pre-processing includes removing outliers and denoising;
  • Decide the number of points k and calculate the slope from the pre-processed data;
  • Set slope thresholds τ 1 , τ 2 for flowing and shut-in break points.
c.
Experiment result
The author made two experiments using this method, and Figure 6 shows the detection results based on the denoised data. The vertical line indicates the detected position of break points.
d.
Advantages and Disadvantages
The pressure-slope-based method is a simple method for detecting large magnitude break points in data. It relies on pre-processing the data and selecting an appropriate threshold. The effectiveness of the approach can be affected by the quality of pre-processing and threshold selection. Human involvement is necessary in these steps.

2.2.2. Pressure Derivative Method

a.
Detection theory
Rai proposed using the derivative of pressure as a criterion to detect breakpoints in a pressure curve [19]. The Savitzky-Golay filter is used to smooth and transform the discrete pressure data into a continuous curve. The derivative of this curve is then calculated and used as the source data for breakpoint detection. The derivative of pressure exhibits peaks when there are significant changes in the gradient of the data. These peaks in the first and higher derivatives indicate the vicinity of the breakpoints location.
b.
Detection procedure
The detection procedure can be detailed as follows:
  • Smoothen the original pressure data using the Savitzky–Golay filter by adjusting the size of the filtering window S and the order of the polynomial θ 1 ;
  • Calculate the pressure derivative and select the proper order of derivative θ 2 ;
  • Detect break points by derivative values bigger than the threshold τ 1 or smaller than the threshold τ 2 ;
  • Adjust break points position by checking slope changes.
c.
Experiment result
The author conducted one experiment using this method, and Figure 7 shows the detection results based on the denoised data. The red dot indicates the detected position of break points.
d.
Advantages and Disadvantages
The method can effectively detect significant breakpoints. However, it requires many parameters to be set and involves many human interactions. Additionally, the method may fail to detect a few breakpoints with small magnitudes, as seen in the author’s experiments.

2.3. Pattern Recognition Methods

The pattern recognition method is a widely used concept for detecting changes and identifying pressure transients in PDG data [32]. There are two specific methods that have been developed for this purpose. Both methods involve pre-defining two patterns for detection. Suzuki [33] uses the responses of the pressure derivative as patterns to detect breakpoints, while Olsen [34] represents the time window containing breakpoints as images. One of the benefits of using pattern recognition methods is that they require minimal human intervention once the patterns have been properly set. They can also be used for real-time detection. However, these methods can be difficult to generalize to different PDG datasets, and pre-processing must be done carefully to avoid over-smoothing, which can negatively impact detection performance.

2.3.1. Filter Convolution Method

a.
Detection theory
Suzuki’s method is a break-point detection method that uses filter convolution to detect specific patterns in the pressure derivative rather than just changes in frequency. This method can also be combined with noise-removal filters to improve its tolerance to noise, as described in the reference [33].
b.
Detection procedure
The detection procedure can be detailed as follows:
  • Pre-process data for outlier removal and denoising;
  • Calculate the pressure derivative;
  • Determine the criterion for patterns by deciding the size of time window S in the filter function according to the minimum buildup interval desired to detect;
  • Detect the break points by locating the maxima of the filter function bigger than threshold τ .
c.
Experiment result
The author conducted one experiment with oil well PDG data and two experiments with gas well PDG data by using this method. The author mentioned that the detection result was better with the oil well PDG data than with gas wells. Figure 8 shows the detection results based on the denoised data with oil well PDG data. The vertical line indicates the detected position of break points.
d.
Advantages and Disadvantages
The Suzuki method is designed to be tolerant to noise to some extent through the use of a filter function. However, the results of the experiment indicate that it is difficult to pre-define universal patterns for different datasets. The threshold criteria selection requires human interaction and is a trial-and-error process.

2.3.2. Image Pattern Method

a.
Detection theory
In the article [35], Olsen proposed a pattern recognition method for transient identification, which is further elaborated on in his Ph.D. thesis [34]. The method pre-defines two patterns for break points, one for build-up, as shown in Figure 9, and one for draw-down, and represents the break points as binary images (0 or 1). The break points will be detected when the shape made by the points belongs to the pre-defined transient pattern. To remove the smallest transient patterns, the author proposed a complementary measure that uses statistical validation criteria. This may cause some of the smallest transients to be missed, but it also decreases the number of faulty detected transients.
b.
Detection procedure
The detection procedure can be detailed as follows:
  • Pre-process data for outlier removal and denoising;
  • Pre-define patterns considering various reservoir models and the influence of noise;
  • Choose the number of measurements N that represent the image for pattern recognition; decide the pressure interval Δ P in each window;
  • Keep the detected break points by discarding the ones below a certain signal standard deviation value τ .
c.
Experiment result
The author conducted an experiment with synthetic data and two experiments with real data using this method. The detection rate was 86% for the synthetic data and 80% for the real data. The missed detections were mainly due to small pressure changes. The results of the experiment with real data are shown in Figure 10, where the red vertical line indicates the position of the detected break points. The results are based on denoised data.
d.
Advantages and Disadvantages
The Olsen’s method uses pre-defined patterns for break-point detection and can be used for real-time automatic detection once the patterns are properly defined, requiring less user interaction during the detection stage. The method can detect most break points, but may miss some with small magnitude, as seen in the author’s experiments. Data pre-processing is necessary before detection, but care must be taken to avoid oversmoothing the data.

2.4. Segmentation Method

a.
Detection theory
Rai proposed a segmentation method for break-point detection in [19] which requires only one user input parameter for identifying the break points. The method is based on time series segmentation from Wallace [36] and involves solving a sequence of maximum orthogonal (Euclidean) distance problems at strategic points, as illustrated in Figure 11. The author also attempted to improve the performance by incorporating flow rate data into the segmentation method as a variant.
b.
Detection procedure
The detection procedure can be detailed as follows:
  • Label the first and the last points in the dataset as two strategic points;
  • Find the third point with the maximum orthogonal distance from the line segment joining the two strategic points. Put the third point in the list of strategic points;
  • Perform this iterative calculation until the greatest orthogonal distance is smaller than a prescribed threshold τ .
c.
Experiment result
The author conducted the experiment using this method with the same data using the pressure derivative method mentioned in the pressure data trend type with the derivative method. Figure 12 shows the detection results. The red dots indicate the detected position of break points.
d.
Advantages and Disadvantages
The Rai’s method can detect significant break points without pre-processing or denoising, which is an advantage over other methods. Additionally, it requires less human intervention, as only one threshold is needed to decide the number of break points. However, it may fail to detect some break points with small magnitude, and there are no similar methods found to process PDG data without denoising, so more tests are needed to prove the method’s feasibility. Additionally, the author improved the method’s performance by incorporating flow rate data, but this may not always be available.

3. Comparison of the Methods

Seven methods for transient identification have been reviewed and grouped into four types. Figure 13 shows the workflow of wavelet decomposition methods and pressure trend methods. For these two types of methods, the detection procedure is similar, and the difference is mainly related to the step of break-point detection.
Figure 14 shows the workflow of pattern recognition methods. Compared with Workflow 1 (Figure 13), the workflow of pattern recognition has the same four steps, but the difference is that the pre-defined patterns’ function should be defined and work as a threshold. The pattern recognition workflow may be used for real-time break-point detection.
Figure 15 shows the workflow of the segmentation method. Compared with Workflow1 (Figure 13) and Workflow2 (Figure 14), the threshold is still the crucial influencer for the detection, and the difference is in avoiding pre-processing and position adjustment, which is an advantage for automation.
By comparison, the most common thing among the three workflows reviewed is the need for data pre-processing and threshold selection, which requires human interaction in all of the seven methods. The wavelet decomposition methods, pressure derivative methods, and segmentation method need to adjust the threshold in the detection stage, whereas the pattern recognition methods set the conditions in the pre-defined patterns acting similarly to the threshold. If the threshold is set too high, then only break points with high-pressure change amplitude are detected, which may ignore many lower-amplitude transients damaging subsequent analysis. If the threshold is set too low, more break points with small amplitude are detected, which may disturb the subsequent analysis with redundant transients identified. Pattern recognition methods pre-define the patterns, so it is more efficient and suited for real-time detection compared to the rest of the methods since human interaction is needed only before the analysis. All seven methods discussed above are compared with their main ideas, aspects of positive and negative sides, and key user inputs in Table 1.

4. Discussions

  • Pre-processing is considered a necessary step by most methods for transient identification. However, pre-processing should be cautiously used because the processed data may lose the signal of the original data. There is no clear answer to decide on the level of processing.
  • Selection of a threshold is crucial for transient identification in many methods reviewed. How to minimize its effect is still an unsolved challenge.
  • Transient identification includes considerations of shut-in and flowing transients. From the reviewed literature, both types of identification use the same methods. However, a combined identification methods, which mean identification shut-in transients and multi-rate transients by using different methods may help in improving the identification results.
  • The methods reviewed in this article need to be further tested on different datasets. The capabilities of the methods tested on gas or oil wells needs to be studied and further tested on different cases of wells producing or injecting different fluids (like liquid and gases).
  • Integrating the methods reviewed with machine learning (ML, including hybrid ML guided by physical models) and other relevant methods is worth to be tested in order to evaluate their suitability for automated transient identification.

5. Conclusions

The modern well surveillance systems now focus on stable and reliable transient identification methods capable of processing large volumes of PDG data and requiring minimum human interaction with low processing time for on-the-fly well monitoring and control. In this context, the literature review provided analysis and guidance on the most promising methods. The following results and conclusions may be reported based on the review and analysis of the transient identification methods available in the literature:
  • The common terms used for processing pressure measurements from permanent downhole gauges (PDG) are defined based on a critical review of the concepts and the terms available in the literature.
    • A pressure transient is a predominantly monotonic change of pressure in response to a sudden change of rate.
    • A break point is a point in time separating two transients.
    • Transient identification is a process of dividing pressure time series into sequential transients based on the objectives of the data interpretation.
  • Seven methods for transient identification have been reviewed and grouped into four types.
    • Wavelet decomposition methods can effectively detect significant transients, but it needs human interaction for better performance.
    • Pressure data trend methods can be easily implemented and understood but require high-quality data pre-processing.
    • Pattern recognition methods can be used for real-time transient identification without human interaction during the identification. Human interaction is however needed before the identification of pre-defined patterns.
    • The segmentation method can detect transients without pre-processing, but threshold selection is needed for better performance as for the wavelet decomposition above.
  • The segmentation and pattern recognition methods may be suggested as the best candidates in the context of automated workflow development. The segmentation method may be considered the best choice taking the simplicity of its application and reduced workload for pre-processing. The main advantage of the pattern recognition method is the ability to use it without human interaction, although based on human pre-defined patterns.

Author Contributions

Conceptualization, A.S. and B.C.; methodology, A.S., L.C. and B.C.; validation, V.D., N.Z. and A.S.; investigation, L.C., A.S. and B.C.; writing—original draft preparation, L.C., B.C., N.Z. and A.S.; writing—review and editing, B.C., N.Z., V.D., C.R. and A.S.; visualization, B.C.; supervision, A.S., V.D. and C.R.; project administration, C.R. and A.S.; funding acquisition, A.S. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the AutoWell project funded by the Research Council of Norway and the industry partners including ConocoPhillips Skandinavia, Lundin Energy Norway, Sumitomo Corporation Europe Norway Branch, Wintershall Dea Norge, and Aker BP (grant no. 326580) for support in preparation and publishing of this paper. The interpretations and conclusions presented in this paper are those of the authors and not necessarily those of the industry partners of the AutoWell project.

Data Availability Statement

This review article operates only with figures and data available in the original papers reviewed.

Acknowledgments

This paper is based on the study initiated within the MSc thesis of L. Chen [37] and continued in the Ph.D. research by B. Cui. The authors acknowledge K. Muradov and T. Wiktorski for valuable comments and feedback on the paper. The copyright holders are acknowledged for permissions to use separate figures from the original publications in this review article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shchipanov, A.A.; Berenblyum, R.A.; Kollbotn, L. Pressure Transient Analysis as an Element of Permanent Reservoir Monitoring. In Proceedings of the SPE Annual Technical Conference and Exhibition, Amsterdam, The Netherlands, 27–29 October 2014. [Google Scholar] [CrossRef]
  2. Soeriawinata, T.; Kelkar, M. Reservoir Management Using Production Data. In Proceedings of the SPE Oklahoma City Oil and Gas Symposium/Production and Operations Symposium, Oklahoma City, OK, USA, 28–31 March 1999. [Google Scholar] [CrossRef]
  3. Kamal, M.; Tian, C.; Suleen, S. Use of Transient Tests to Monitor Progress of Flooding in IOR/EOR Operations. In Proceedings of the SPE Annual Technical Conference and Exhibition, Dubai, United Arab Emirates, 26–28 September 2016. [Google Scholar] [CrossRef]
  4. Ribeiro, P.M.; Horne, R.N. Pressure and Temperature Transient Analysis: Hydraulic Fractured Well Application. In Proceedings of the SPE Annual Technical Conference and Exhibition, New Orleans, LA, USA, 2 October 2013. [Google Scholar] [CrossRef]
  5. Yousef, A.A.; Gentil, P.; Jensen, J.L.; Lake, L.W. A Capacitance Model To Infer Interwell Connectivity From Production- and Injection-Rate Fluctuations. SPE Reserv. Eval. Eng. 2006, 9, 630–646. [Google Scholar] [CrossRef]
  6. Sankaran, S.; Wright, D.; Gamblin, H.; Kumar, D. Creating Value by Implementing an Integrated Production Surveillance and Optimization System—An Operator’s Perspective. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 9–11 October 2017. [Google Scholar] [CrossRef]
  7. Ahmadi, M.A.; Ebadi, M.; Hosseini, S.M. Prediction breakthrough time of water coning in the fractured reservoirs by implementing low parameter support vector machine approach. Fuel 2014, 117, 579–589. [Google Scholar] [CrossRef]
  8. Chorneyko, D.M. Real-Time Reservoir Surveillance Utilizing Permanent Downhole Pressures—An Operator’s Experience. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 24–27 September 2006. [Google Scholar] [CrossRef]
  9. Horne, R.N. Listening to the Reservoir—Interpreting Data From Permanent Downhole Gauges. J. Pet. Technol. 2007, 59, 78–86. [Google Scholar] [CrossRef]
  10. Shchipanov, A.; Kollbotn, L.; Berenblyum, R. Characterization and monitoring of reservoir flow barriers from pressure transient analysis for CO2 injection in saline aquifers. Int. J. Greenh. Gas Control 2019, 91, 102842. [Google Scholar] [CrossRef]
  11. Athichanagorn, S.; Horne, R.N.; Kikani, J. Processing and Interpretation of Long-Term Data Acquired From Permanent Pressure Gauges. SPE Reserv. Eval. Eng. 2002, 5, 384–391. [Google Scholar] [CrossRef]
  12. Nomura, M. Processing and Interpretation of Pressure Transient Data from Permanent Downhole Gauges. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 2006. [Google Scholar]
  13. Cumming, J.A.; Wooff, D.A.; Whittle, T.; Gringarten, A.C. Multiwell Deconvolution. SPE Reserv. Eval. Eng. 2014, 17, 457–465. [Google Scholar] [CrossRef]
  14. Liu, Y.; Horne, R.N. Interpreting Pressure and Flow-Rate Data From Permanent Downhole Gauges by Use of Data-Mining Approaches. SPE J. 2012, 18, 69–82. [Google Scholar] [CrossRef]
  15. Liu, Y.; Horne, R.N. Interpreting Pressure and Flow Rate Data from Permanent Downhole Gauges Using Convolution-Kernel-Based Data Mining Approaches. In Proceedings of the SPE Western Regional Meeting, Monterey, CA, USA, 19–25 April 2013. [Google Scholar] [CrossRef]
  16. Tian, C.; Horne, R.N. Applying Machine Learning Techniques to Interpret Flow Rate, Pressure and Temperature Data From Permanent Downhole Gauges. In Proceedings of the SPE Western Regional Meeting, Garden Grove, CA, USA, 27–30 April 2015. [Google Scholar] [CrossRef]
  17. Tian, C.; Horne, R.N. Recurrent Neural Networks for Permanent Downhole Gauge Data Analysis. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 9–11 October 2017. [Google Scholar] [CrossRef]
  18. Kin, K. Permanent Downhole Gauge Data Interpretation. Master’s Thesis, Stanford University, Stanford, CA, USA, 2001. [Google Scholar]
  19. Rai, H.; Horne, R.N. Analyzing Simultaneous Rate and Pressure Data from Permanent Downhole Gauges. In Proceedings of the SPE Annual Technical Conference and Exhibition, Anaheim, CA, USA, 11–14 November 2007. [Google Scholar] [CrossRef]
  20. Liu, Y. Interpreting Pressure and Flow Rate Data from Permanent Downhole Gauges Using Data Mining Approaches. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 2013. [Google Scholar]
  21. Thomas, A.O. The Data as the Model: Interpreting Permanent Downhole Gauge Data without Knowing the Reservoir Model. Master’s Thesis, Stanford University, Stanford, CA, USA, 2002. [Google Scholar]
  22. Athichanagorn, S. Development of an Interpretation Methodology for Long-Term Pressure Data from Permanent Downhole Gauges. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 1999. [Google Scholar]
  23. Li, X. Processing and Analysis of Transient Data from Permanent Down-Hole Gauges (PDG). Ph.D. Thesis, Heriot-Watt University, Edinburgh, Scotland, 2009. [Google Scholar]
  24. Tian, C. Machine Learning Approaches for Permanent Downhole Gauge Data Interpretation. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 2018. [Google Scholar]
  25. Shchipanov, A.A.; Kollbotn, L.; Prosvirnov, M. Step Rate Test as a Way to Understand Well Performance in Fractured Carbonates. In Proceedings of the 79th EAGE Annual Conference & Exhibition, Paris, France, 12–15 June 2017. [Google Scholar] [CrossRef]
  26. Mallat, S.; Hwang, W. Singularity detection and processing with wavelets. IEEE Trans. Inf. Theory 1992, 38, 617–643. [Google Scholar] [CrossRef]
  27. Crouse, M.; Nowak, R.; Baraniuk, R. Wavelet-based statistical signal processing using hidden Markov models. IEEE Trans. Signal Process. 1998, 46, 886–902. [Google Scholar] [CrossRef] [Green Version]
  28. Ouyang, L.B.; Kikani, J. Improving Permanent Downhole Gauge (PDG) Data Processing via Wavelet Analysis. In Proceedings of the European Petroleum Conference, Aberdeen, UK, 29–31 October 2002. [Google Scholar] [CrossRef]
  29. Soliman, M.Y.; Ansah, J.; Stephenson, S.; Mandal, B. Application of Wavelet Transform to Analysis of Pressure Transient Data. In Proceedings of the SPE Annual Technical Conference and Exhibition, New Orleans, LA, USA, 30 September–3 October 2001. [Google Scholar] [CrossRef]
  30. Viberti, D.; Verga, F.; Delbosco, P. An Improved Treatment of Long-Term Pressure Data for Capturing Information. SPE Reserv. Eval. Eng. 2007, 10, 359–366. [Google Scholar] [CrossRef]
  31. Maragos, P. Slope transforms: Theory and application to nonlinear signal processing. IEEE Trans. Signal Process. 1995, 43, 864–877. [Google Scholar] [CrossRef]
  32. Culea-Florescu, A.; Aiordachioaie, D. Some results on change detection based on advanced signal processing paradigm. In Proceedings of the 2016 8th International Conference on Electronics, Computers and Artificial Intelligence (ECAI), Ploiesti, Romania, 30 June–2 July 2016. [Google Scholar] [CrossRef]
  33. Suzuki, S.; Chorneyko, D. Automatic Detection of Pressure-Buildup Intervals From Permanent Downhole Pressure Data Using Filter Convolution. In Proceedings of the SPE Annual Technical Conference and Exhibition, New Orleans, LA, USA, 4–7 October 2009. [Google Scholar] [CrossRef]
  34. Olsen, S. On Automatic Data Processing and Well-Test Analysis in Real-Time Reservoir Management Applications. Ph.D. Thesis, University of Bergen, Bergen, Norway, 2011. [Google Scholar]
  35. Olsen, S.; Nordtvedt, J.E. Automatic Filtering and Monitoring of Real-Time Reservoir and Production Data. In Proceedings of the SPE Annual Technical Conference and Exhibition, Dallas, TX, USA, 9–12 October 2005. [Google Scholar] [CrossRef]
  36. Wallace, W.E.; Kearsley, A.J.; Guttman, C.M. An Operator-Independent Approach to Mass Spectral Peak Identification and Integration. Anal. Chem. 2004, 76, 2446–2452. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, L. Transient Identification and Interpretation of Real-Time Pressure Data from Modern Well Surveillance Systems. Master’s Thesis, University of Stavanger, Stavanger, Norway, 2022. [Google Scholar]
Figure 1. A common location of Permanent downhole gauge(PDG) in the wellbore.
Figure 1. A common location of Permanent downhole gauge(PDG) in the wellbore.
Energies 16 01689 g001
Figure 2. Detection of break points in details data by wavelet decomposition.
Figure 2. Detection of break points in details data by wavelet decomposition.
Energies 16 01689 g002
Figure 3. Detection result using Wavelet Modulus Maxima from Athichanagorn [22].
Figure 3. Detection result using Wavelet Modulus Maxima from Athichanagorn [22].
Energies 16 01689 g003
Figure 4. Detection result using the data grouping method from Li [23].
Figure 4. Detection result using the data grouping method from Li [23].
Energies 16 01689 g004
Figure 5. Detection of break points in detail data by the pressure trend method.
Figure 5. Detection of break points in detail data by the pressure trend method.
Energies 16 01689 g005
Figure 6. Detection result using the pressure slope from Viberti [30].
Figure 6. Detection result using the pressure slope from Viberti [30].
Energies 16 01689 g006
Figure 7. Detection result using the pressure derivative from Rai [19].
Figure 7. Detection result using the pressure derivative from Rai [19].
Energies 16 01689 g007
Figure 8. Detection result using the filter convolution method from Suzuki [33].
Figure 8. Detection result using the filter convolution method from Suzuki [33].
Energies 16 01689 g008
Figure 9. Pre-define pattern for build-up break points.
Figure 9. Pre-define pattern for build-up break points.
Energies 16 01689 g009
Figure 10. Detection result using the image pattern method from Olsen [34].
Figure 10. Detection result using the image pattern method from Olsen [34].
Energies 16 01689 g010
Figure 11. Points with maximum orthogonal distance.
Figure 11. Points with maximum orthogonal distance.
Energies 16 01689 g011
Figure 12. Detection result using the segmentation method from Rai [19].
Figure 12. Detection result using the segmentation method from Rai [19].
Energies 16 01689 g012
Figure 13. Workflow1: Wavelet decomposition and pressure trend.
Figure 13. Workflow1: Wavelet decomposition and pressure trend.
Energies 16 01689 g013
Figure 14. Workflow2: Pattern recognition.
Figure 14. Workflow2: Pattern recognition.
Energies 16 01689 g014
Figure 15. Workflow3: Segmentation method.
Figure 15. Workflow3: Segmentation method.
Energies 16 01689 g015
Table 1. Comparison among different methods for identification.
Table 1. Comparison among different methods for identification.
Main IdeaAdvantagesDisadvantagesKey Input
Wavelet modulus maximaDetect break points by detailed signal higher than threshold after wavelet decomposition.Significant break points can be detected.Pre-processing is needed. Some false break points are detected. Human intervention to adjust threshold.Decomposition level L, Threshold τ
Data groupingDetect break points by statistical analysis in the detailed signal after wavelet decomposition.Both big and small break points can be detected. The detected position is more precise.Pre-processing is needed. Some false break points are detected. Human intervention to adjust threshold.Decomposition level L, Threshold τ
Pressure slopeDetect break points by slope value higher than threshold.Easy to implement. Significant break points can be detected.Pre-processing is needed. Some false break points are detected. Human intervention to adjust threshold. Some break points missed.Number of points N, Threshold τ
Pressure derivativeDetect break points by derivative value higher than threshold.Significant break points can be detected. Filtering is part of pre-processing.Pre-processing is needed. Many parameters to choose. Human interaction to adjust threshold. Some break points missedWindow size S, Polynomial order θ 1 , Derivative order θ 2 , Threshold τ
Filter convonlutionDetect break points by processed signal higher than threshold.Significant break points can be detected. Tolerant to low level noise. Less human interaction.Pre-processing is needed. Sensitive to high noise level. Human interaction to adjust threshold.Patterns criterion C Window size S, Threshold τ
Pattern recognitionDetect break points by the shape belongs to pre-defined patternsSignificant break points can be detected. Require less human interaction.Pre-processing is needed. Failed to detect break points with small magnitude. Caution to avoid oversmoothing.Patterns criterion C, Number of points N, Pressure difference Δ P
SegmentationDetect break points by calculating data distance smaller than threshold.Significant break points can be detected. Denoising may not be required. Less human interaction.Outliers affect result. Some break points missed. Human interaction to adjust threshold. Need to confirm the detected positionThreshold τ
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

Cui, B.; Chen, L.; Zhang, N.; Shchipanov, A.; Demyanov, V.; Rong, C. Review of Different Methods for Identification of Transients in Pressure Measurements by Permanent Downhole Gauges Installed in Wells. Energies 2023, 16, 1689. https://0-doi-org.brum.beds.ac.uk/10.3390/en16041689

AMA Style

Cui B, Chen L, Zhang N, Shchipanov A, Demyanov V, Rong C. Review of Different Methods for Identification of Transients in Pressure Measurements by Permanent Downhole Gauges Installed in Wells. Energies. 2023; 16(4):1689. https://0-doi-org.brum.beds.ac.uk/10.3390/en16041689

Chicago/Turabian Style

Cui, Boyu, Lejun Chen, Nan Zhang, Anton Shchipanov, Vasily Demyanov, and Chunming Rong. 2023. "Review of Different Methods for Identification of Transients in Pressure Measurements by Permanent Downhole Gauges Installed in Wells" Energies 16, no. 4: 1689. https://0-doi-org.brum.beds.ac.uk/10.3390/en16041689

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