Deconvolution of Well-Test Data as a Nonlinear Total Least-Squares Problem
- T. von Schroeter (Imperial College) | F. Hollaender (Imperial College) | A.C. Gringarten (Imperial College)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- December 2004
- Document Type
- Journal Paper
- 375 - 390
- 2004. Society of Petroleum Engineers
- 5.6.4 Drillstem/Well Testing, 1.10 Drilling Equipment, 7.2.2 Risk Management Systems, 4.1.5 Processing Equipment, 4.1.2 Separation and Treating, 5.1.2 Faults and Fracture Characterisation, 5.3.2 Multiphase Flow, 5.8.8 Gas-condensate reservoirs
- 8 in the last 30 days
- 1,527 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
We present a new time-domain method for the deconvolution of well test data which is characterized by three novel features: (1) Instead of the rate-normalized pressure derivative itself, we estimate its logarithm, which makes explicit sign constraints necessary; (2) the formulation accounts for errors in both rate and pressure data, and thus amounts to a Total Least Squares (TLS) problem; and (3) regularization is based on a measure of the overall curvature of its graph. The resulting separable nonlinear TLS problem is solved using the Variable Projection algorithm. A comprehensive error analysis is given. The paper also includes tests with a simulated example and an application to a large field example.
With current trends towards permanent downhole instrumentation, continuous bottomhole well pressure monitoring is becoming the norm in new field developments. The resulting well-test data sets, recorded mainly during production, can consist of hundreds of flow periods and millions of pressure data points stretched over thousands of hours of recording time. Such data sets contain information about the reservoir at distances from the well which can be several orders of magnitude larger than the radius of investigation of a single flow period.
Conventional derivative analysis is thus ill equipped to access the full potential information content. What is required is an analysis method which can extract the response which the reservoir would exhibit when subjected to a single drawdown at constant rate over any period of time up to the entire production period. In mathematical terms, this is a deconvolution problem. Since its first formulation by Hutchinson and Sikora in 1959,1 it has received sporadic, but recurring attention.2-17
This paper presents a new approach which is based on a regularized, nonlinear TLS formulation. It is an update on earlier versions which were presented at the SPE Annual Meetings in 200118 and 200219 (henceforth referred to as "Paper I" and "Paper II"). More recently, our approach was taken up by Levitan,20 who subjected it to a critical evaluation and suggested some modifications.
In terms of the usual classification into time-domain and spectral approaches, ours is a time-domain approach. It differs from earlier approaches in this category in three important ways:
1. The solution is encoded in terms of the logarithm of the rate-normalized pressure derivative, which automatically ensures strict positivity of the derivative itself at the expense of rendering the problem nonlinear. However, we are thus able to avoid explicit constraints on the solution space which made previous constrained approaches so difficult, yet still cannot prevent zeros in the deconvolved derivative.
2. A new error measure accounts for uncertainties not only in the pressure, but also in the rate data, which are usually much less accurately known. Thus, provided sufficient data are available, our method can provide a joint estimate of initial pressure, rates, and response parameters; the time-consuming manual correction of rate errors is rendered obsolete. The mathematical formulation is an instance of what is known as a TLS problem in the numerical analysis literature and as an "Errors-In-Variables" problem in statistics. TLS has become a standard approach in parameter estimation problems, but its application to well-test analysis seems to be new.
3. Regularization is based on a measure of the total curvature of the deconvolved pressure derivative, instead of its average slope, as in an earlier approach15 and Paper I. Here, the motivation is that slopes provide important information about the flow regime and should therefore be preserved as much as possible.
The paper is organized as follows: The first two sections are introductory and give a summary of the deconvolution problem in well-test analysis and a concise survey of its treatment in the petroleum engineering and hydrology literature. Based on the mathematical framework developed in these sections, we then give a comprehensive account of our own approach. We also derive analytic expressions for bias and variance of the estimated parameter set based on simple Gaussian models for the measurement errors in pressure and rate signals. We illustrate our method with a small simulated data set, demonstrating the effect of varying levels of regularization on the confidence intervals. The final section presents an application to a large field example which allows a direct comparison of our method with conventional derivative analysis.
|File Size||635 KB||Number of Pages||16|