This application claims priority to Russian Application No. 2014135163 filed Aug. 28, 2014, which is incorporated herein by reference in its entirety.
The disclosure relates to the field of studying oil and gas wells and is designed to correct results of pressure measurements in high-rate wells, performed during well tests.
At present, significant improvement in the accuracy and reliability of pressure sensors takes place. The pressure sensors have a resolution better than 0.01 psi and an absolute measurement accuracy of the order of 1 psi. This is especially significant in study of high-productive formations when rapid increase in a pressure is followed by very slow increase in the pressure for a long time. In this case, temperature variations which occur during a well test can change even the nature of pressure variation during the second—slow—stage: the pressure can decrease and not increase in a shut-in well.
Investigation of oil/gas formation properties by analyzing dynamics of pressure variation in a well at a change in the productivity of the well (including said change when the well is shut in) is called as Pressure Transient Analysis (PTA). The traditional techniques for processing data obtained in the PTA do not take into account a possible influence of non-isothermal effects and assume that a pressure sensor is at a fixed distance (50 to 100 m) from an upper boundary of a formation to be studied. Modern techniques for processing the PTA data (cf., Kuchuk F. J., Onur M., Hollaender F. Monograph series: Vol. 57. Pressure Transient Formation and Well Testing (1st ed.). ELSEVIER, 2010, pp. xv-xx, 23-26) comprise stimulating an inflow of a fluid from a formation, taking a sample and registering a flow rate and a pressure, preferably a bottom hole pressure. The PTA data is currently interpreted using analytical solutions of the piezoconductivity equation at various well completion schemes and various boundary conditions. The developed analytical techniques are applicable to fractured, layered, horizontally and radially composite formations. Numerical simulation techniques are used for even more complicated heterogeneous systems and multiphase flows. In the majority of cases, however, all analytical and numerical methods assume the flow of a fluid under isothermal conditions. With this, reduction in a pressure, as logged in some high-rate wells, shows that the traditional isothermal techniques for interpretation of measurement results are not applicable in such wells.
The technical result provided by the invention is an enhanced accuracy of downhole pressure measurement due to taking non-isothermal effects into account.
In accordance with the method, a pressure in a wellbore and a temperature at least at one point of the wellbore are measured during wellbore testing. Then, transient profiles of temperature along the wellbore are determined during the wellbore testing and changes in a density of a downhole fluid and in a length of a tubing when the wellbore is shut in are calculated. The pressure measurement results are corrected on the basis of the calculated changes in the density of the downhole fluid and in the length of the tubing.
The transient profiles of the temperature along the wellbore are determined by measuring the temperature along the wellbore using a system of sensors distributed on the wellbore at different depths or by numerical or analytical simulating the temperature profile. If necessary, the continuous transient profiles of the temperature can be obtained by interpolating the measured temperatures.
The pressure is measured using at least one sensor disposed at a fixed depth within the wellbore.
The disclosure is illustrated by the drawings wherein:
a shows change in an average temperature of a downhole fluid between a tubing anchorage point and a pressure sensor;
b shows a density of the downhole fluid between the sensor and an upper boundary of a reservoir;
c shows variation in a tubing length (the right graph) in time when the wellbore is shut in;
A wellbore testing is divided into two periods: an inflow period and a pressure recovery period. In the first period, the wellbore is opened to operate at a constant or variable flow rate, wherein a downhole pressure drops. For a time depending upon targets of the wellbore testing, formation properties, a bottom-hole zone state, reservoir fluid properties, and a pressure draw-down, the wellbore is shut in, and the pressure is recovered up to its original value.
Pressure and temperature sensors are usually placed on a tubing approximately at a distance L0=100 M above an upper boundary of a reservoir. Since a temperature of a fluid to be produced is essentially higher than an average temperature of overlying rocks, heating of rocks near the wellbore takes place. After shutting the wellbore in, a temperature in the vicinity of the borehole decreases, a fluid filling the wellbore is cooled down to a geothermal temperature, and an average temperature of said fluid decreases. Calculations show that the average temperature of the downhole fluid can decrease down to 30-40° C. in case of deep wells (4,000-6,000 m).
This circumstance can result in effects as follows:
1) a density of the fluid filling the wellbore in an interval between the pressure sensor and the reservoir increases;
2) a length of the tubing (anchored at the surface) changes and a position of the sensor relative to the reservoir changes.
Both said effects leads to reduction in measured pressure.
A pressure change related to a downhole temperature change can be calculated by the formula (1):
ΔP(t)=g·[
re ΔP(t) is a pressure change, Pa;
An increase of the fluid density below the pressure sensor is determined by a volumetric thermal expansion coefficient of the fluid, αf, K−1, and by a change in the average temperature of the fluid below the sensor in the shut-in well, ΔT(t), K:
A decrease in the tubing length, ΔL(t), and respective change in the sensor position is determined by a linear thermal expansion coefficient of the tubing, αt, K−1, and by a variation in a temperature profile through a depth (from an surface anchorage point to the sensor) in the shut-in well, m; g is a gravity acceleration, m/s2. To determine the value ΔL(t), it is proposed to divide a distance between the tubing anchorage point and the pressure sensor into n equal sections and calculate a reduction in the tubing length by the formula (3):
re Ti(t) is a temperature of ith section at a time
t (i=1, . . . n) in the shut-in well,
; Toi a temperature of ith section immediately after shutting the wellbore in,
.
For quantitative estimations of both effects, it is necessary to know a dependency of the downhole fluid temperature upon the depth at different times after shutting the wellbore in. Transient (non-stationary) profiles of the downhole temperature should be measured during the wellbore testing or obtained as a result of the numerical or analytical simulation.
We propose to calculate a corrected downhole pressure value (from which the influence of the fluid density and tubing length variation is excluded) by a formula (4):
P
g
r(t)=Pg(t)+ΔP(t), (4)
where Pg (t) is a measured pressure, Pa, ΔP(t) is a correction taking an influence of temperature effects into account and being calculated by the formulae (1)-(3).
In case of high-productive formations, a pressure in a shut-in well rapidly increases up to short of an initial pressure. Hereupon, there is a long and slow increase in the pressure, and the temperature effects considered above can have an essential influence upon the dynamics of said increase. In some cases, the measured pressure can decrease in time.
A solid line (Non-isothermal) shows calculation results for the same wellbore testing but with taking the non-isothermal effects into account. It was supposed that the pressure sensor is on the tubing at a depth of 100 m above the upper boundary of the reservoir. According to the calculation, an average temperature of a fluid filling the tubing has been decreased by 37° C. during the pressure recovery, while a density of a fluid in the wellbore has been increased by 10 kg/m3.
The decrease in the average temperature of the downhole fluid in the shut-in wellbore in the present case is 37° C. The processing of the pressure curve (Non-isothermal) resulted from measurement by standard isothermal techniques in the shut-in wellbore at such a pressure change gives an incorrect formation model, its permeability and a well skin factor.
Thus, to exclude the influence of the thermal effects, it is necessary to correct the results of measuring the downhole pressure.
The possibility to correct the measured pressure using the method of the disclosure is shown by synthetic examples prepared using the numerical simulator T-Mix (Ramazanov, A. Sh., et al. Termogidravlicheskie issledovania v skvazhine dlya opredelenia parametrov priskvazhinnoi zony plasta i debitov mnogoplastovoi sistemy” (Thermohydrodynamic Studies In Well For Determining Parameters Of Formation Nearfield and Multilayered System Flow Rates), 2010, SPE 136256). It is a code allowing simulation of transient pressure and temperature distributions when a single-phase fluid flows in a formation and in a wellbore, said code being able to reproduce an arbitrary sequence of production operations in the wellbore: the production beginning, the flow rate variation, the wellbore shutting in, etc.
The transient pressure distribution in the formation is simulated using the Darcy filtration law for a cylindrical gas flow or a poorly-compressible liquid in a layered medium. The downhole pressure is calculated using the quasi-stationary law of conservation of momentum with taking into account a friction force, a gravity force, an acceleration, and an effect of a wellbore filled with a compressible fluid.
The transient formation temperature field is calculated taking into account the conductive and convective heat transfer, the adiabatic effect, and the Joule-Thomson effect. The transient thermal model of the wellbore takes into account the fluid mixing effect, the heat exchange between the well and surrounding rocks as well as the adiabatic effect and the fluid heating due to viscous friction forces.
According to the method described above, a transient temperature distribution along the wellbore is measured or simulated; adequacy of the simulation is controlled by matching of simulated transient pressure, temperature and flow rate with available measurement results of said quantities for the whole testing duration.
The temperature, the pressure, and the flow rate were calculated for values of parameters of the wellbore, the reservoir and the sequence of operations described above, using the numerical simulator T-Mix.
The simulation had the following parameters:
properties of the reservoir: homogeneous, the reservoir thickness—100 m; the initial reservoir pressure at the outer boundary of the radius of 1500 m was fixed to be 7,251.89 psi; the reservoir permeability—2 d; the temperature—120° C.; the well skin—5; the wellbore depth—4,000 m;
properties of the fluid: oil having the density of 800 kg/m3 under formation conditions, thermal conductivity of 0.14 W/m/K, specific heat capacity of 2000 J/kg/K, viscosity of 1 cP, compressibility of 6.9×10−6 psi−1.
The sequence of production operations in the wellbore is as follows: circulation for 70 hrs; shut-in for 70 hrs; production at the flow rate of 2,000 M3 per day for 20 hrs (70 to 90 hrs,
The second step is to obtain the transient downhole temperature profiles when the wellbore is shut in, said profiled being derived from numerical calculations using a model having input parameters that give the best match with available measurements.
The resulted transient (non-stationary) temperature profiles are used to calculate change in the downhole fluid density and the tubing length when the wellbore is shut in. The distance between the tubing anchorage point and the pressure sensor was divided into 78 equal sections. The average downhole fluid temperature change in time,
between the tubing anchorage point and the pressure sensor, the downhole fluid density variation between the sensor and the upper boundary of the reservoir, and the tubing length reduction in the shut-in wellbore were calculated by the formulae (2) (3). The results of simulating the wellbore shut-in mode for 90 to 120 hr are shown in
The final step is correction of the downhole pressure measurement results by the formula (4) taking into account the obtained results with respect to the downhole fluid density and the tubing length in order to exclude the influence of non-isothermal effects.
Number | Date | Country | Kind |
---|---|---|---|
2014135163 | Aug 2014 | RU | national |