An Investigation of Recent Deconvolution Methods for Well-Test Data Analysis
- Mustafa Onur (Istanbul Technical U.) | Murat Cinar (Istanbul Technical U.) | Dilhan Ilk (Texas A&M University) | Peter P. Valko (Texas A&M University) | Thomas A. Blasingame (Texas A&M University) | Peter S. Hegeman (Schlumberger)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- June 2008
- Document Type
- Journal Paper
- 226 - 247
- 2008. Society of Petroleum Engineers
- 5.6.4 Drillstem/Well Testing, 5.3.2 Multiphase Flow, 5.5.8 History Matching, 5.1.2 Faults and Fracture Characterisation, 4.3.4 Scale, 3.3 Well & Reservoir Surveillance and Monitoring, 4.1.2 Separation and Treating, 5.9.2 Geothermal Resources, 3.2.3 Hydraulic Fracturing Design, Implementation and Optimisation, 4.1.5 Processing Equipment, 4.6 Natural Gas, 5.1.5 Geologic Modeling
- 2 in the last 30 days
- 1,422 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 5.00|
|SPE Non-Member Price:||USD 35.00|
In this work, we present an investigation of recent deconvolution methods proposed by von Schroeter et al. (2002, 2004), Levitan (2005) and Levitan et al. (2006), and Ilk et al. (2006a, b). These works offer new solution methods to the long-standing deconvolution problem and make deconvolution a viable tool for well-test and production-data analysis. However, there exists no study presenting an independent assessment of all these methods, revealing and discussing specific features associated with the use of each method in a unified manner. The algorithms used in this study for evaluating the von Schroeter et al. and Levitan methods represent our independent implementations of their methods based on the material presented in their papers, not the original algorithms implemented by von Schroeter et al. and Levitan. Three synthetic cases and one field case are considered for the investigation.
Our results identify the key issues regarding the successful and practical application of each method. In addition, we show that with proper care and attention in applying these methods, deconvolution can be an important tool for the analysis and interpretation of variable rate/pressure reservoir performance data.
Applying deconvolution for well-test and production data analysis is important because it provides the equivalent constant rate/pressure response of the well/reservoir system affected by variable rates/pressures (von Schroeter et al. 2002, 2004; Levitan 2005; Levitan et al. 2006; Ilk et al. 2006a, b; Kuchuk et al. 2005). With the implementation of permanent pressure and flow-rate measurement systems, the importance of deconvolution has increased because it is now possible to process the well test/production data simultaneously and obtain the underlying well/reservoir model (in the form of a constant rate pressure response). New methods of analyzing well-test data in the form of a constant-rate drawdown system response and production data in the form of a constant-pressure rate system response have emerged with development of robust pressure/rate (von Schroeter et al. 2002, 2004; Levitan 2005; Levitan et al. 2006; Ilk et al. 2006a, b) and rate/pressure (Kuchuk et al. 2005) deconvolution algorithms. In this work, we focus on the pressure/rate deconvolution for analyzing well-test data.
For over a half century, pressure/rate deconvolution techniques have been applied to well-test pressure and rate data as a means to obtain the constant-rate behavior of the system (Hutchinson and Sikora 1959; Coats et al. 1964; Jargon and van Poollen 1965; Kuchuk et al. 1990; Thompson and Reynolds 1986; Baygun et al. 1997). A thorough review and list of the previous deconvolution algorithms can be found in von Schroeter et al. (2004). The primary objective of applying pressure/rate deconvolution is to convert the pressure data response from a variable-rate test or production sequence into an equivalent pressure profile that would have been obtained if the well were produced at a constant rate for the entire duration of the production history.
If such an objective could be achieved with some success, then, as stated by Levitan, the deconvolved response would remove the constraints of conventional analysis techniques (Earlougher 1977; Bourdet 2002) that have been built around the idea of applying a special time transformation [e.g., the logarithmic multirate superposition time (Agarwal 1980)] to the test pressure data so that the pressure behavior observed during individual flow periods would be similar in some way to the constant-rate system response. As also stated by Levitan, the superposition-time transform does not completely remove all effects of previous rate variations and often complicates test analysis because of residual superposition effects.
Unfortunately, deconvolution is an ill-posed inverse problem and will usually not have a unique solution even in the absence of noise in the data. Even if the solution is unique, it is quite sensitive to noise in the data, meaning that small changes in input (measured pressure and rate data) can lead to large changes in the output (deconvolved) result. Therefore, this ill-posed nature of the deconvolution problem combined with errors that are inherent in pressure and rate data makes the application of deconvolution a challenge, particularly so in terms of developing robust deconvolution algorithms which are error-tolerant. Although there exists a variety of different deconvolution algorithms proposed in the past, only those developed by von Schroeter et al., Levitan, and Ilk et al. appear to offer the necessary robustness to make deconvolution a viable tool for well-test and production data analysis. In this paper, our objectives are to conduct an investigation of these three deconvolution methods and to establish the advantages and limitations of each method.
As stated in the abstract, the algorithms used in this study for evaluating the von Schroeter et al. and Levitan methods represent our independent implementations based on the material presented in their papers; therefore, our implementations may not be identical to their versions. However, as is shown later, validation conducted on the simulated (test) data sets (von Schroeter et al. 2004; Levitan 2005) sent to us directly by von Schroeter and Levitan shows that our implementations reproduce almost exactly the same results generated by their original algorithms for these simulated data sets.
The paper is organized as follows: First, we describe the pressure/rate deconvolution model and error model considered in this work. Then, we provide the mathematical background of the von Schroeter et al., Levitan, and Ilk et al. methods together with their specific features. We compare the performance of each method by considering three synthetic and one field well-test data sets. Finally, we provide a discussion of our results obtained from this investigation.
|File Size||4 MB||Number of Pages||22|
Abate, J. and Valkó, P.P. 2004. Multi-Precision Laplace TransformInversion. International Journal for Numerical Methods inEngineering 60 (5): 979-993. doi: 10.1002/nme.995.
Abbaszadeh, M. and Cinco-Ley, H. 1995. Pressure-Transient Behavior in aReservoir With a Finite-Conductivity Fault. SPEFE 10 (1):26-32; Trans., SPE, 299. SPE-24704-PA doi: 10.2118/24704-PA
Agarwal, R.G. 1980. A New MethodTo Account for Producing Time Effects When Drawdown Type Curves Are Used ToAnalyze Buildup and Other Test Data. Paper SPE 9289 presented at the SPEAnnual Technical Conference and Exhibition, Dallas, 21-24 September. doi:10.2118/9289-MS
Bard, Y. 1974. Nonlinear Parameter Estimation. New York City:Academic Press, Inc.
Baygun, B., Kuchuk, F.J., and Arikan, O. 1997. Deconvolution Under NormalizedAutocorrelation Constraints. SPEJ 2 (3): 246-253.SPE-28405-PA doi: 10.2118/28405-PA
Björck, Å. 1996. Numerical Methods for Least Square Problems.Philadelphia, Pennsylvania: Society for Industrial and Applied Mathematics.
Bourdet, D. 2002. Well Test Analysis: The Use of Advanced InterpretationModels. Amsterdam: Elsevier Science B.V.
Bourdet, D., Ayoub, J.A., and Pirard, Y.M. 1989. Use of Pressure Derivative inWell-Test Interpretation. SPEFE 4 (2): 293-302;Trans., SPE, 287. SPE-12777-PA doi: 10.2118/12777-PA
Carvalho, R., Thompson, L.G., Redner, R., and Reynolds, A.C. 1996. Simple Procedures for ImposingConstraints for Nonlinear Least Squares Optimization. SPEJ 1(4): 395-402. SPE-29582-PA doi: 10.2118/29582-PA
Cheney, E.W. and Kincaid, D.R. 2003. Numerical Mathematics andComputing. Pacific Grove, California: Brooks Cole Publishing Company.
Çinar, M., Ilk, D., Onur, M., Valkó, P.P., and Blasingame, T.A. 2006. A Comparative Study of Recent RobustDeconvolution Algorithms for Well-Test and Production-Data Analysis. PaperSPE 102575 presented at the SPE Annual Technical Conference and Exhibition, SanAntonio, Texas, 24-27 September. doi: 10.2118/102575-MS
Coats, K.H., Rapoport, L.A., McCord, J.R., and Drews, W.P. 1964. Determination of Aquifer InfluenceFunctions From Field Data. JPT 16 (12): 1417-1424;Trans., AIME, 231. SPE-897-PA doi: 10.2118/897-PA
Dennis, J.E. and Schnabel, R.B. 1983. Numerical Methods for UnconstrainedOptimization and Nonlinear Equations. Englewood Cliffs, New Jersey:Prentice-Hall Inc.
Earlougher, C.R. Jr. 1977. Advances in Well Test Analysis. MonographSeries, SPE, Richardson, Texas, 5.
Golub, G.H., Hansen, P.C., and O'Leary, D.P. 1999. Tikhonov Regularization andTotal Least Squares. SIAM J. Matrix Anal. Appl. 21 (1):185-194. doi: 10.1137/S0895479897326432.
Hegeman, P.S., Hallford, D.L, and Joseph, J.A. 1993. Well Test Analysis With ChangingWellbore Storage. SPEFE 8 (3): 201-207; Trans., SPE,295. SPE-21829-PA doi: 10.2118/21829-PA
Hutchinson, T.S. and Sikora, V.J. 1959. A GeneralizedWater-Drive Analysis. Trans., AIME 216: 169-178.SPE-1123-G.
Ilk, D. 2005. Deconvolution of Variable Rate Reservoir Performance DataUsing B-Splines. MS thesis, College Station, Texas: Texas A&MUniversity.
Ilk, D., Valko, P.P., and Blasingame, T.A. 2006a. Deconvolution of Variable-RateReservoir-Performance Data Using B-Splines. SPEREE 9 (4):582-595. SPE-95571-PA doi: 10.2118/95571-PA
Ilk, D., Anderson, D.M., Valko, P.P., and Blasingame, T.A. 2006b. Analysis of Gas-Well ReservoirPerformance Data Using B-Spline Deconvolution. Paper SPE 100573 presentedat the SPE Gas Technology Symposium, Calgary, 15-17 May. doi:10.2118/100573-MS
Jargon J.R. and Van Poollen, H.K. 1965. Unit Response Function FromVarying-Rate Data. JPT 17 (8): 965-969; Trans., AIME,234. SPE-981-PA doi: 10.2118/981-PA
Kuchuk, F.J. 1990. Applicationsof Convolution and Deconvolution to Transient Well Tests. SPEFE5 (4): 375-384; Trans., SPE, 289. SPE-16394-PA doi:10.2118/16394-PA
Kuchuk, F.J., Carter, R.G., and Ayestaran, L. 1990. Deconvolution of Wellbore Pressureand Flow Rate. SPEFE 5 (1): 53-59. SPE-13960-PA doi:10.2118/13960-PA
Kuchuk, F.J., Hollaender, F., Gok, I.M., and Onur, M. 2005. Decline Curves from Deconvolution ofPressure and Flow-Rate Measurements for Production Optimization andPrediction. Paper SPE 96002 presented at the SPE Annual TechnicalConference and Exhibition, Dallas, 9-12 October. doi: 10.2118/96002-MS
Levitan, M.M. 2005. PracticalApplication of Pressure/Rate Deconvolution to Analysis of Real Well Tests.SPEREE 8 (2): 113-121. SPE-84290-PA doi: 10.2118/84290-PA
Levitan, M.M., Crawford, G.E., and Hardwick, A. 2006. Practical Considerations forPressure-Rate Deconvolution of Well-Test Data. SPEJ 11 (1):35-47. SPE-90680-PA doi: 10.2118/90680-PA
Onur, M. et al. 2007. Analysis of Well Tests in Afyon Ömer-Gecek GeothermalField, Turkey. Proc., 32nd Workshop on Geothermal Reservoir Engineering,Stanford University, Stanford, California, 22-24 January.
Press, H. Flannery, B.P., Teukolsky, S.A., and Vetterling, W.T. 1986.Numerical Recipes: The Art of Scientific Computing. Cambridge, UK:Cambridge University Press.
Saphir Ver. 3.20.10. software. 2005. Sophia Antipolis, France: KappaEngineering.
Thompson, L.G. and Reynolds, A.C. 1986. Analysis of Variable-Rate Well-TestPressure Data Using Duhamel's Principle. SPEFE 1 (5):453-469. SPE-13080-PA doi: 10.2118/13080-PA
van Everdingen, A.F. and Hurst, W. 1949. TheApplication of the Laplace Transformation to Flow Problems in Reservoirs.Trans., AIME 186: 305-324. SPE-949305-G.
von Schroeter, T., Hollaender, F., and Gringarten, A.C. 2004. Deconvolution of Well Test Data as aNonlinear Total Least-Squares Problem. SPEJ 9 (4): 375-390.SPE-77688-PA doi: 10.2118/77688-PA