US20150053398A1 - Method for determining an inflow profile of multilayer reservoir fluids in a wellbore - Google Patents
Method for determining an inflow profile of multilayer reservoir fluids in a wellbore Download PDFInfo
- Publication number
- US20150053398A1 US20150053398A1 US14/465,978 US201414465978A US2015053398A1 US 20150053398 A1 US20150053398 A1 US 20150053398A1 US 201414465978 A US201414465978 A US 201414465978A US 2015053398 A1 US2015053398 A1 US 2015053398A1
- Authority
- US
- United States
- Prior art keywords
- wellbore
- perforation
- temperature
- zone
- measured
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004519 manufacturing process Methods 0.000 claims abstract description 61
- 238000009529 body temperature measurement Methods 0.000 claims abstract description 16
- 230000000694 effects Effects 0.000 claims abstract description 15
- 238000001816 cooling Methods 0.000 claims abstract description 12
- 238000004088 simulation Methods 0.000 claims description 9
- 238000005553 drilling Methods 0.000 claims description 6
- 239000011435 rock Substances 0.000 description 14
- 239000000243 solution Substances 0.000 description 12
- 230000015572 biosynthetic process Effects 0.000 description 11
- 238000005755 formation reaction Methods 0.000 description 11
- 230000035699 permeability Effects 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 6
- 238000012546 transfer Methods 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 2
- 230000036962 time dependent Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical compound C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- HJUFTIJOISQSKQ-UHFFFAOYSA-N fenoxycarb Chemical compound C1=CC(OCCNC(=O)OCC)=CC=C1OC1=CC=CC=C1 HJUFTIJOISQSKQ-UHFFFAOYSA-N 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- E21B47/0003—
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B36/00—Heating, cooling or insulating arrangements for boreholes or wells, e.g. for use in permafrost zones
- E21B36/001—Cooling arrangements
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
- E21B47/103—Locating fluid leaks, intrusions or movements using thermal measurements
Definitions
- the disclosure relates to the field of geophysical exploration of oil and gas wells, namely to determining a profile of fluid inflow to a wellbore from producing layers of multilayer reservoirs.
- temperature is measured along a wellbore under quasi-steady production conditions and formation temperature near the wellbore is supposed to be close to a temperature of the undisturbed formation. It is assumed that rocks temperature near the wellbore is equal to the undisturbed rock mass temperature.
- a production rate of a producing layer is usually determined from flow temperature measured below and above the layer and from temperature of a fluid flowing from the formation into the wellbore (see, for example, Hill, A.D., 2002. Production Logging—Theoretical and Interpretive Elements. SPE Monograph Series).
- RU Patent No. 2474687 describes a method to determine the fluids inflow profile in multilayer wells based on results of measuring temperatures in the well immediately after perforation when temperature of rocks near the well is reduced due to previous technological operations (well drilling and circulation). According to this method, production rates of individual layers are calculated based on temperature rates of change measured in the wellbore.
- a theoretical basis for this method is a linear dependence established between the produced fluid temperature change rate and a formation specific production rate, which assumes cooling of a bottom-hole formation zone and a shut-in between the fluid circulation in the wellbore and perforation.
- a disadvantage of this method is a relatively short period of production, when this linear dependence between the temperature change rate and production rate is retained, which can restrict applicability of this method to determine the inflow profile and a need to measure the temperature change rate in the well before its perforation, which increases a probable error and imposes additional restrictions on the perforation schedule.
- the disclosure provides improved accuracy and reliability of determining an inflow profile in a multilayer wellbore at an initial production stage immediately after the wellbore perforation. In this case, there is no need for a shut-in between well circulation and perforation and it is not required to measure temperature change rate in the wellbore before its perforation.
- the method comprises cooling of a bottom-hole zone before perforation of a wellbore. Then the wellbore is perforated, and the flow temperature is measured in the wellbore above each perforation zone.
- a production rate of each pay zone is determined taking into account thicknesses of the perforation zones and using the results of temperature measurements acquired between the end of an initial production stage characterized by strong impact of a volume of the wellbore and quick flow temperature changes in the wellbore, and a time when the bottom-hole zone cooling impact on the temperature measurements becomes insignificant.
- the bottom-hole zone can be cooled by drilling or well circulation.
- the flow temperature in the wellbore is measured by means of sensors installed on a tubing string used for the perforation, above each perforation zone.
- T n,a (t) is a flow temperature measured in the wellbore above the n perforation zone
- t i is time moments equally distributed within an interval t 1 ⁇ t ⁇ t 2 , where t 1 is the end time of the initial production stage characterized by strong impact of a volume of the wellbore and quick flow temperature changes in the wellbore, t 2 is a time when the bottom-hole zone cooling impact on the temperature measurements becomes insignificant; a number of the time moments within this interval is equal to a number of the temperature measurements.
- Q 1 Q 1 + Y 2 1 - Y 2 + Y 3 ( 1 - Y 2 ) ⁇ ( 1 - Y 3 ) + ... + Y m ( 1 - Y 2 ) ⁇ ( 1 - Y 3 ) ⁇ ... ⁇ ( 1 - Y m )
- Q 2 Q 1 ⁇ Y 2 1 - Y 2
- Q 3 Q 1 ⁇ Y 3 ( 1 - Y 2 ) ⁇ ( 1 - Y 3 )
- Q n Q 1 ⁇ Y n ( 1 - Y 2 ) ⁇ ( 1 - Y 3 ) ⁇ ... ⁇ ( 1 - Y n )
- the wellbore flow temperature T n b (t) is measured under each perforation zone, with F(Y n , t) determined as follows:
- pressure is additionally measured in the wellbore, below all perforation zones; and to solve an inverse task, a numerical simulation is used with such numerical model parameters which ensure coincidence of measured and calculated temperature-time dependencies; the measured wellbore pressure is applied as a boundary condition for simulation of pressure and temperature fields in the productive layers.
- temperature measurement results are used, which are acquired preferably within an interval of 1 to 10 hours after start of production.
- FIG. 1 shows a scheme of a wellbore with two perforation zone and three temperature sensors
- FIG. 2 a shows calculated radial temperature distributions after the well circulation and shut-in period
- FIG. 2 b shows calculated temperatures of fluids produced from layers with different specific production rates (full lines) and a temperature of a layer with a high flow rate after scaling of the timescale (markers);
- FIG. 3 and FIG. 4 show calculated temperatures for various values of layer permeability
- FIG. 5 and FIG. 6 show dependencies of mis-ties on the formation relative production rate for various values of formation permeability
- FIG. 7 b shows dependence of mis-tie S on value Y 2 .
- Solution Y 2min 0.2;
- FIG. 8 b shows dependence of mis-tie S on value Y 2 .
- Solution Y 2min 0.69;
- FIG. 9 a shows values of selected permeability to solve an inverse task.
- FIG. 9 b shows specified temperatures and temperatures calculated when solving the inverse task.
- k 1 96 mD
- k 2 29 mD
- FIG. 10 a shows values of selected permeability to solve the inverse task.
- FIG. 10 b shows specified temperatures and temperatures calculated when solving the inverse task.
- k 1 29 mD
- k 2 96 mD
- the proposed method can be applied in conditions of typical perforation with a tubing string.
- a bottom-hole zone is cooled during a well drilling. If drilling was completed recently (several days before perforation), special cleanout of the wellbore is not required. If the wellbore was drilled a long time ago, special cleanout may be performed before the perforation, which ensures significant (by 5-10° C. and over) decrease in temperature of rocks near the wellbore as compared with the undisturbed rock temperature.
- the proposed method uses the fact that an initial radial distribution of temperature is approximately the same in all productive layers under consideration T ⁇ T 0 (r).
- V q 2 ⁇ ⁇ ⁇ r ( 2 )
- ⁇ f c f a fluid volumetric heat capacity
- ⁇ m c m a volumetric heat capacity of rock matrix
- ⁇ reservoir porosity
- r a distance to an axis of the wellbore.
- T ⁇ ( r , t ) T 0 ( r 2 + ⁇ ⁇ ⁇ q ⁇ t ) , ( 3 )
- T 0 (r) is a radial distribution of temperature in the reservoir before start of production, and q is a specific production rate of the reservoir.
- T in ⁇ ( t ) T 0 ( r w 2 + ⁇ ⁇ ⁇ q ⁇ t ) , ( 4 )
- r w is a radius of the wellbore
- a relative production rate of an upper productive layer Y 2 can be found by selecting such a scaling factor ‘a 2 ’ that will provide coincidence of time-dependent temperatures of fluids flowing from different productive layers.
- FIG. 2 This is illustrated ( FIG. 2 ) by results of numerical calculations made by means of software COMSOL MULTIPHYSICS 3.5TM. Calculated radial distributions of temperatures after the well circulation (24 hours of circulation at initial rock temperature of 100° C. and circulating fluid temperature of 50° C. (the dotted curve), and its shut-in for 24 hours (solid curve)) are shown on FIG. 2 a .
- a thick solid line matches a production rate 2 ⁇ q. Markers that actually coincide with the thin curve show a result of twofold time-based expansion for the curve matching the flow rate of 2 ⁇ q.
- Formulas (5) and (6) enable to find production rates of individual perforation zones, if all time dependent temperatures are measured for fluids flowing to the wellbore from different layers. In practice, only a wellbore flow temperature is measured. A temperature measured above a lower perforation zone is approximately equal to a temperature of fluid flowing from the lower zone, but for all other zones, temperature T in is normally unknown.
- the disclosure proposes the following method to solve the problem.
- relative production rates of productive layers can be approximately determined by means of the stationary energy conservation law that expresses a balance of thermal energy coming into the inflow zone and energy leaving this zone:
- ⁇ c is a fluid volumetric conductivity
- T in,2 ( t ) T in,1 ( a 2 ⁇ t ) ⁇ T 1,a ( a 2 ⁇ t )
- this equation does not include an unknown temperature of produced fluid T in2 (t) and an unknown temperature of undisturbed rocks T f .
- the unknown value Y 2 is determined solely by results of wellbore temperature measurements: T 1,a , T 2,b , T 2,a
- a relative production rate of the upper layer can be found by equation (9) for any fixed time moment ‘t’; however not all time moments are equivalent.
- T 1a , T 2b and T 2a are measured with a certain error. This error is connected with sensors and measurement method errors.
- calculation formula (10) includes bulk flow temperatures while the sensor measures a temperature in a certain flow point and there is some temperature distribution along the flow cross-section. Therefore, to reduce impacts of various factors on results of the inflow profile determination, it is expedient to use all data acquired within time interval t 1 ⁇ t ⁇ t 2 .
- the described above method to determine the inflow profile for two productive layers is easily generalized and can be used for arbitrary number m of layers.
- T n,a (t) is the fluid temperature measured in the wellbore above the n perforation zone
- t i is time moments equally distributed within interval t 1 ⁇ t ⁇ t 2
- t 1 is the end time of the initial production stage characterized by strong impact of the volume of the wellbore and quick fluid temperature changes in the wellbore
- t 2 is a time when the bottom-hole zone cooling effect on the temperature measurements becomes insignificant (of the same order that Joule-Thomson and geotherm effects); besides, a number of time moments in this interval is equal to a number of temperature measurements.
- Production rates of individual productive layers determined by means of formulas (12)-(15) can be considered as an approximate solution of the problem. In some instances, when the above conditions of the analytical model applicability are violated, an error in determination of individual production rates can reach 10-20% and more (see below).
- a possibility to determine the inflow profile using the proposed method was demonstrated based on synthetic test cases that were generated by means of the transient numerical model of associated processes of heat and mass transfer in the layers and the wellbore.
- the model functionality enables simulate arbitrary alternation of various processes: well circulation, shut-in, production, and fluid injection.
- the wellbore numerical model was tested for many years on analytical solutions, using commercial simulators (COMSOL® and ECLIPSE®) and was successfully applied for simulation and interpretation of complicated field cases.
- brackets are relative production rates for the upper zone in the cases under consideration.
- FIGS. 3 , 4 show calculated temperatures T 1,a , T 2,b , T 2,a , and T in,2 for these options.
- Temperature T 1,a is equal to a flow temperature in the point which is 1 m above the upper boundary of the perforation zone, T 2,b ⁇ 1 m below the lower boundary of the perforation zone. It is seen that in this case a difference between the temperatures T 1,a and T 2,b is very low since the well production rate is rather high, the layers are close to each other, and fluid energy losses to the surrounding rock are insignificant. This means that in this case, there is no need in a temperature sensor below the upper formation (T 2,b ), and formula (11) can be used for function F (Y 2 ,t).
- mis-tie S (Y 2 ) (11) can be calculated for a time interval of 3 to 7 hours.
- FIGS. 7 , 8 show that the proposed method to determine the inflow profile is resistant to random errors of wellbore temperature measurements.
- Solution of inverse tasks by using noisy temperature data (with a standard deviation of 0.1 K) produces practically the same result as with using the initial data.
- accuracy of the proposed method to determine the inflow profile based on temperature data can be insufficient. It can be connected with violation of conditions that were used to derive calculation formulas. For example, with long distances between layers, initial radial distributions of temperature in the layers will be different, geothermal gradient will have a significant impact; with high pressure differentials between a layer and the wellbore, Joule-Thomson effect may not be neglected; with low-value times of observations, wellbore volume impact may not be neglected, etc.
- FIGS. 9 , 10 illustrate application of this method to determine the inflow profile.
- permeability values of two productive layers were selected automatically, to ensure coincidence of measured and calculated temperatures.
- Permeability values acquired while applying the simplified model were used as initial approximations. Based on FIG. 9 , 10 , it is seen that practically full agreement of temperatures was obtained after about 150 solutions of the direct task. Obtained values of permeability and flow rates actually match the target ones:
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Measuring Volume Flow (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
Methods for determining an inflow profile of multilayer reservoir fluids in a wellbore are described herein. A bottom-hole zone is cooled before well perforation. Then, the wellbore is perforated, and a flow temperature is measured in the wellbore above each perforation zone. A production rate of each productive layer is determined taking into account thicknesses of the perforation zones and using results of temperature measurements acquired in a period between an end of an initial production stage characterized by strong impact of a volume of the wellbore and quick flow temperature changes in the wellbore, and a time when the bottom-hole zone cooling effect on the temperature measurements becomes insignificant.
Description
- This application claims priority to Russian Application No. 2013139149 filed Aug. 23, 2013, which is incorporated herein by reference in its entirety.
- The disclosure relates to the field of geophysical exploration of oil and gas wells, namely to determining a profile of fluid inflow to a wellbore from producing layers of multilayer reservoirs.
- Typically, when estimating production rates of individual producing layers based on temperature data, temperature is measured along a wellbore under quasi-steady production conditions and formation temperature near the wellbore is supposed to be close to a temperature of the undisturbed formation. It is assumed that rocks temperature near the wellbore is equal to the undisturbed rock mass temperature. A production rate of a producing layer is usually determined from flow temperature measured below and above the layer and from temperature of a fluid flowing from the formation into the wellbore (see, for example, Hill, A.D., 2002. Production Logging—Theoretical and Interpretive Elements. SPE Monograph Series).
- Traditional methods to determine an inflow profile based on temperature data apply the following simplifying assumptions: a quasi-steady-state fluid flow in the wellbore, a Joule-Thomson effect constant, and a temperature of rocks near the wellbore is determined by an undisturbed geotherm.
- These assumptions are not valid if temperature in the wellbore is measured at the initial production stage immediately after perforation. In the first turn, the assumption about the temperature of rocks near the wellbore is wrong. As a rule, during perforation, the temperature near the wellbore is significantly less than the temperature of undisturbed rocks. This is due to cooling effects of previous technological operations in the well: well drilling and circulation.
- RU Patent No. 2474687 describes a method to determine the fluids inflow profile in multilayer wells based on results of measuring temperatures in the well immediately after perforation when temperature of rocks near the well is reduced due to previous technological operations (well drilling and circulation). According to this method, production rates of individual layers are calculated based on temperature rates of change measured in the wellbore. A theoretical basis for this method is a linear dependence established between the produced fluid temperature change rate and a formation specific production rate, which assumes cooling of a bottom-hole formation zone and a shut-in between the fluid circulation in the wellbore and perforation.
- A disadvantage of this method is a relatively short period of production, when this linear dependence between the temperature change rate and production rate is retained, which can restrict applicability of this method to determine the inflow profile and a need to measure the temperature change rate in the well before its perforation, which increases a probable error and imposes additional restrictions on the perforation schedule.
- The disclosure provides improved accuracy and reliability of determining an inflow profile in a multilayer wellbore at an initial production stage immediately after the wellbore perforation. In this case, there is no need for a shut-in between well circulation and perforation and it is not required to measure temperature change rate in the wellbore before its perforation.
- The method comprises cooling of a bottom-hole zone before perforation of a wellbore. Then the wellbore is perforated, and the flow temperature is measured in the wellbore above each perforation zone. A production rate of each pay zone is determined taking into account thicknesses of the perforation zones and using the results of temperature measurements acquired between the end of an initial production stage characterized by strong impact of a volume of the wellbore and quick flow temperature changes in the wellbore, and a time when the bottom-hole zone cooling impact on the temperature measurements becomes insignificant.
- The bottom-hole zone can be cooled by drilling or well circulation.
- The flow temperature in the wellbore is measured by means of sensors installed on a tubing string used for the perforation, above each perforation zone.
- Starting from a second perforation zone relative production rates of the perforation zones Yn, n=2, 3, . . . m are calculated, successively upward, by means of minimization
-
- where
-
- Qn (n=2, 3, . . . m) are the production rates of the perforation zones
-
- hn (n=1, 2, . . . m) is a thickness of a n perforation zone, Tn,a(t) is a flow temperature measured in the wellbore above the n perforation zone, ti is time moments equally distributed within an interval t1<t<t2, where t1 is the end time of the initial production stage characterized by strong impact of a volume of the wellbore and quick flow temperature changes in the wellbore, t2 is a time when the bottom-hole zone cooling impact on the temperature measurements becomes insignificant; a number of the time moments within this interval is equal to a number of the temperature measurements.
- The production rates Qn (n=2, 3, . . . m) of the perforation zones are determined by formulas:
-
- where Q is a total volumetric well production rate.
- According to one of the embodiments of the disclosure, the wellbore flow temperature Tn b(t) is measured under each perforation zone, with F(Yn, t) determined as follows:
-
- According to another embodiments of the disclosure n, pressure is additionally measured in the wellbore, below all perforation zones; and to solve an inverse task, a numerical simulation is used with such numerical model parameters which ensure coincidence of measured and calculated temperature-time dependencies; the measured wellbore pressure is applied as a boundary condition for simulation of pressure and temperature fields in the productive layers.
- To determine the production rates of the productive layers, temperature measurement results are used, which are acquired preferably within an interval of 1 to 10 hours after start of production.
- The disclosure is illustrated by drawings where:
-
FIG. 1 shows a scheme of a wellbore with two perforation zone and three temperature sensors; -
FIG. 2 a shows calculated radial temperature distributions after the well circulation and shut-in period; -
FIG. 2 b shows calculated temperatures of fluids produced from layers with different specific production rates (full lines) and a temperature of a layer with a high flow rate after scaling of the timescale (markers); -
FIG. 3 andFIG. 4 show calculated temperatures for various values of layer permeability; -
FIG. 5 andFIG. 6 show dependencies of mis-ties on the formation relative production rate for various values of formation permeability; -
FIG. 7 a shows calculated temperatures with normal noise CKO=0.1 K. Option k1=100 mD, k2=30 mD; -
FIG. 7 b shows dependence of mis-tie S on value Y2. Solution Y2min=0.2; -
FIG. 8 a shows calculated temperatures with normal noise CKO=0.1 K. Option k1=30 mD, k2=100 mD; -
FIG. 8 b shows dependence of mis-tie S on value Y2. Solution Y2min=0.69; -
FIG. 9 a shows values of selected permeability to solve an inverse task. Option k1=100 mD, k2=30 mD; -
FIG. 9 b shows specified temperatures and temperatures calculated when solving the inverse task. k1=96 mD, k2=29 mD, Y2≈0.23; -
FIG. 10 a shows values of selected permeability to solve the inverse task. Option k1=30 mD, k2=100 mD; -
FIG. 10 b shows specified temperatures and temperatures calculated when solving the inverse task. k1=29 mD, k2=96 mD, Y2≈0.75. - The proposed method can be applied in conditions of typical perforation with a tubing string.
- Normally, a bottom-hole zone is cooled during a well drilling. If drilling was completed recently (several days before perforation), special cleanout of the wellbore is not required. If the wellbore was drilled a long time ago, special cleanout may be performed before the perforation, which ensures significant (by 5-10° C. and over) decrease in temperature of rocks near the wellbore as compared with the undisturbed rock temperature.
- The proposed method uses the fact that an initial radial distribution of temperature is approximately the same in all productive layers under consideration T≈T0(r).
- This assumption is valid if:
- (a) the productive layers are close to each other (at a distance of 30-50 m); and
(b) all productive layers have approximately the same thermal properties. - After start of production, the radial temperature profile in a reservoir and a temperature of fluid flowing to the wellbore from the formation are mainly determined by convection heat transfer (1):
-
- where
-
- is a fluid filtration rate, q [m3/m/s] is a production rate per 1 m of perforation zone, ρfcf is a fluid volumetric heat capacity, ρrcr=φ·ρfcf+(1−φ)·ρmcm is a volumetric heat capacity of the fluid saturated reservoir, ρmcm is a volumetric heat capacity of rock matrix, φ is reservoir porosity, and r is a distance to an axis of the wellbore.
- Equation (1) does not take into account impact of conductive heat transfer, Joule-Thomson effect, and adiabatic effect. Negligible impact of conductive heat transfer in the inflow zone for typical specific production rates was demonstrated by means of detailed numeric simulation. Relative impact of Joule-Thomson and adiabatic effects is determined by pressure differential between the reservoir and the wellbore and by typical decrease of temperature in the bottom-hole area. Since pressure differentials at the initial production stage are normally insignificant (10-30 bar), and the reservoir cooling reaches 10 K and over, these effects may be neglected in the first approximation. The inflow profile can be determined more accurately based on transient temperature logging data by means of numerical simulation (see below).
- Solution of equation (1) is as follows:
-
- where
-
- T0(r) is a radial distribution of temperature in the reservoir before start of production, and q is a specific production rate of the reservoir.
- From (3), a temperature-time dependence expression follows for the fluid flowing into the wellbore:
-
- where rw is a radius of the wellbore.
- Let us consider a wellbore with two productive layers (
FIG. 1 ). Temperatures of fluids flowing from different layers with specific production rates q1 and q2 are as follows (respectively): -
- It is obvious that temperature Tin,2(t) can be written as
-
- According to formula (6), a relative production rate of an upper productive layer Y2 can be found by selecting such a scaling factor ‘a2’ that will provide coincidence of time-dependent temperatures of fluids flowing from different productive layers.
- This is illustrated (
FIG. 2 ) by results of numerical calculations made by means of software COMSOL MULTIPHYSICS 3.5™. Calculated radial distributions of temperatures after the well circulation (24 hours of circulation at initial rock temperature of 100° C. and circulating fluid temperature of 50° C. (the dotted curve), and its shut-in for 24 hours (solid curve)) are shown onFIG. 2 a. A thin solid curve onFIG. 2 b shows calculated temperature of fluid flowing into the wellbore at a production rate q (q=10 m3/m/day). A thick solid line matches aproduction rate 2·q. Markers that actually coincide with the thin curve show a result of twofold time-based expansion for the curve matching the flow rate of 2·q. - Formulas (5) and (6) enable to find production rates of individual perforation zones, if all time dependent temperatures are measured for fluids flowing to the wellbore from different layers. In practice, only a wellbore flow temperature is measured. A temperature measured above a lower perforation zone is approximately equal to a temperature of fluid flowing from the lower zone, but for all other zones, temperature Tin is normally unknown.
- The disclosure proposes the following method to solve the problem.
- Let us consider a case when temperature sensors are installed on a tubing above and below each perforation zone (
FIG. 1 ). - In this case, relative production rates of productive layers can be approximately determined by means of the stationary energy conservation law that expresses a balance of thermal energy coming into the inflow zone and energy leaving this zone:
-
ρc·Q 1 ·T 2,b +ρc·Q 2 ·T in,2 =ρc·(Q 1 +Q 2)·T 2,a (7) - where ρc is a fluid volumetric conductivity.
- Since we consider transient processes and wellbore fluid temperature varies with time, using this equation is not entirely correct; however, as the below calculations demonstrate, it can be used for an approximate solution of the problem.
- From equation (7) the formula follows for a relative production rate of an upper layer Y2:
-
- Taking into account that
-
T 1,a(t)≈T in,1(t) -
and using formula (5) -
T in,2(t)=T in,1(a 2 ·t)≈T 1,a(a 2 ·t) - we find an equation for the unknown value Y2:
-
- It is essential that this equation does not include an unknown temperature of produced fluid Tin2(t) and an unknown temperature of undisturbed rocks Tf. The unknown value Y2 is determined solely by results of wellbore temperature measurements: T1,a, T2,b, T2,a
- If a distance between layers is short (about 10 m), then, as calculations demonstrate T2,b(t)≈T1,a(t), equation (10) will look as follows:
-
- A relative production rate of the upper layer can be found by equation (9) for any fixed time moment ‘t’; however not all time moments are equivalent.
- At short times (t<t1), the wellbore volume which is not considered by the simplified model described above can significantly impact the calculation results. Besides, at the initial stage, the wellbore fluid temperature quickly changes in the perforation zone and using stationary energy equation (7) can result in significant errors in determination of relative production rates of the layers.
- At long times (t>t2), when the bottom-hole zone cooling effect ceases and the fluid flowing to the wellbore heats the near wellbore rocks actually to the reservoir temperature, Joule-Thomson and geotherm effects become more significant. Both of these effects are not taken into account in the above simplified model, which may result in errors in determination of the inflow profile. Therefore, when the above simplified model is used, data acquired at relatively short times t<t2 is used.
- Values T1a, T2b and T2a are measured with a certain error. This error is connected with sensors and measurement method errors. In reality, calculation formula (10) includes bulk flow temperatures while the sensor measures a temperature in a certain flow point and there is some temperature distribution along the flow cross-section. Therefore, to reduce impacts of various factors on results of the inflow profile determination, it is expedient to use all data acquired within time interval t1<t<t2.
- It is proposed to find value Y2 by minimization of function S (Y2) (11), which is calculated as a sum of squared mis-ties for all time moments in the interval t1<t<t2, for which measured temperatures are available:
- The proposed method to determine the inflow profile was tested in several typical synthetic instances. It was shown, that values t1=1÷2 hours and t2=8÷10 hours ensure a satisfactory accuracy of the inverse task solution.
- The described above method to determine the inflow profile for two productive layers is easily generalized and can be used for arbitrary number m of layers.
- Starting from the second perforation zone relative production rates of perforation zones Yn, n=2, 3, . . . m are calculated, successively upward, by means of minimization
-
- where
-
- Qn (n=2, 3, . . . m) are production rates of perforation zones
-
- hn (n=1, 2, . . . m) is a thickness of a n perforation zone, Tn,a(t) is the fluid temperature measured in the wellbore above the n perforation zone, ti is time moments equally distributed within interval t1<t<t2, t1 is the end time of the initial production stage characterized by strong impact of the volume of the wellbore and quick fluid temperature changes in the wellbore, t2 is a time when the bottom-hole zone cooling effect on the temperature measurements becomes insignificant (of the same order that Joule-Thomson and geotherm effects); besides, a number of time moments in this interval is equal to a number of temperature measurements.
- Production rates (n=2, 3, . . . m) of perforation zones are determined by formulas:
-
- where Q is a total volumetric well production rate.
- Production rates of individual productive layers determined by means of formulas (12)-(15) can be considered as an approximate solution of the problem. In some instances, when the above conditions of the analytical model applicability are violated, an error in determination of individual production rates can reach 10-20% and more (see below).
- A possibility to determine the inflow profile using the proposed method was demonstrated based on synthetic test cases that were generated by means of the transient numerical model of associated processes of heat and mass transfer in the layers and the wellbore. The model functionality enables simulate arbitrary alternation of various processes: well circulation, shut-in, production, and fluid injection. The wellbore numerical model was tested for many years on analytical solutions, using commercial simulators (COMSOL® and ECLIPSE®) and was successfully applied for simulation and interpretation of complicated field cases.
- For verification of the proposed method to determine production rates of individual productive layers by means of the numerical model, the following sequence of technological operations in the wellbore was simulated:
- 1. Well circulation for 24 hours. The wellbore circulating fluid temperature at a depth is assumed 100° C., undisturbed rock temperature is 123° C.
2. Shut-in for 24 hours. This stage is not obligatory for using the proposed method but as a rule there is a certain time interval between the well circulation and perforation.
3. Wellbore perforation and oil production for 24 hours with a total production rate of 100 m3/day. - Thickness of perforation zones is h1=h2=10 m, distance between these zones is 10 m.
- Two cases of formation permeability values were considered (skin factor was taken equal to zero):
-
k 1=100 mD,k 2=30 mD(Y≈0.25) (1) -
k 1=30 mD,k 2=100 mD(Y≈0.75). (2) - In brackets are relative production rates for the upper zone in the cases under consideration.
-
FIGS. 3 , 4 show calculated temperatures T1,a, T2,b, T2,a, and Tin,2 for these options. Temperature T1,a is equal to a flow temperature in the point which is 1 m above the upper boundary of the perforation zone, T2,b−1 m below the lower boundary of the perforation zone. It is seen that in this case a difference between the temperatures T1,a and T2,b is very low since the well production rate is rather high, the layers are close to each other, and fluid energy losses to the surrounding rock are insignificant. This means that in this case, there is no need in a temperature sensor below the upper formation (T2,b), and formula (11) can be used for function F (Y2,t). - Calculations show that in all cases, mis-tie S (Y2) (11) can be calculated for a time interval of 3 to 7 hours.
FIG. 5 shows that in case No. 1, mis-tie S (Y2) has a single minimum at Y2≈0.233. This value is very close to the accurate solution Y2=0.25 (error ˜7%). In case No. 2 (FIG. 6 ) mis-tie S (Y2) also has a single minimum at Y2min=0.68 (error ˜9%). -
FIGS. 7 , 8 show that the proposed method to determine the inflow profile is resistant to random errors of wellbore temperature measurements. Solution of inverse tasks by using noisy temperature data (with a standard deviation of 0.1 K) produces practically the same result as with using the initial data. - In some cases, accuracy of the proposed method to determine the inflow profile based on temperature data can be insufficient. It can be connected with violation of conditions that were used to derive calculation formulas. For example, with long distances between layers, initial radial distributions of temperature in the layers will be different, geothermal gradient will have a significant impact; with high pressure differentials between a layer and the wellbore, Joule-Thomson effect may not be neglected; with low-value times of observations, wellbore volume impact may not be neglected, etc.
- In these cases, to improve accuracy of the inflow profile determination, it is possible, for solution of the inverse task, to use numerical simulation and to select such parameters of the numerical model, which would ensure agreement of measured and calculated time-temperature dependencies. However, the wellbore measured pressure not used above, can be used as a boundary condition for simulation of pressure and temperature fields in pay zones.
-
FIGS. 9 , 10 illustrate application of this method to determine the inflow profile. For the above two synthetic cases, permeability values of two productive layers were selected automatically, to ensure coincidence of measured and calculated temperatures. Permeability values acquired while applying the simplified model were used as initial approximations. Based onFIG. 9 , 10, it is seen that practically full agreement of temperatures was obtained after about 150 solutions of the direct task. Obtained values of permeability and flow rates actually match the target ones: -
k 1=96 mD,k 2=29 mD(Y 2≈0.23). (1) - Target values: k1=100 mD, k2=30 mD (Y2=0.25)
-
k 1=29 mD,k 2=96 mD(Y 2≈0.75). (2) - Target values: k1=30 mD, k2=100 mD (Y2=0.75)
- In this case, two parameters of permeability values were selected while solving the inverse task. Numerical experiments show that a single-value solution of the inverse task can be obtained with an increased number of selected parameters to 6 and more. In particular, along with permeability of individual formations, their skin factors can also be determined.
Claims (8)
1. A method for determining an inflow profile of multilayer reservoir fluids in a wellbore comprising:
cooling of a bottom-hole zone,
perforating the wellbore,
measuring flow temperature in the wellbore above each perforation zone, and
determining a production rate of each productive layer taking into account thicknesses of the perforation zones and using results of temperature measurements acquired between (i) an end of an initial production stage and (ii) a time starting from which the bottom-hole zone cooling effect on the temperature measurements becomes insignificant.
2. The method of claim 1 , wherein the bottom-hole zone is cooled by drilling.
3. The method of claim 1 , wherein the bottom-hole zone is cooled by well circulation.
4. The method of claim 1 , wherein the flow temperature in the wellbore is measured by sensors installed on a tubing string used for perforation above the each perforation zone.
5. The method of claim 1 , wherein for determining the production rates of the productive layers, temperature measurement results are used, which are acquired within an interval of 1 to 10 hours after the start of production.
6. The method of claim 1 , wherein starting from a second perforation zone relative production rates of perforation zones Yn, n=2, 3, . . . m are calculated successively upward by means of minimization according to:
where
(n=2, 3, . . . m) is the production rate of each perforation zone,
hn (n=1, 2, . . . m) is a thickness of a n perforation zone, Tn,a(t) is a flow temperature measured in the wellbore above the n perforation zone, ti is time moments equally distributed within the interval t1<t<t2, where t1 is the end time of the initial production stage, t2 is the time starting from which the bottom-hole zone cooling effect on the temperature measurements becomes insignificant, a number of the time moments in the interval t1<t<t2 is equal to a number of temperature measurements, and the production rates Qn (n=2, 3, . . . m) of the perforation zones are determined formulas according to:
where Q is a total volumetric wellbore production rate.
7. The method of claim 1 , wherein wellbore flow temperature is additionally measured under each perforation zone, and F(Yn, t) is determined according to:
where Tn,b(t) and Tn,a(t) are temperatures measured in the wellbore below and above the n perforation zone.
8. The method of claim 1 , wherein pressure is additionally measured in the wellbore below all perforation zones and an inverse task is solved using a numerical simulation with such numerical model parameters selected which ensure matching of measured and calculated temperature-time dependencies, and the measured wellbore pressure is applied as a boundary condition for simulation of pressure and temperature fields in the production layers.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013139149/03A RU2531499C1 (en) | 2013-08-23 | 2013-08-23 | Method of determining fluid movement profile of stacked pools in well |
RU2013139149 | 2013-08-23 |
Publications (1)
Publication Number | Publication Date |
---|---|
US20150053398A1 true US20150053398A1 (en) | 2015-02-26 |
Family
ID=52479315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/465,978 Abandoned US20150053398A1 (en) | 2013-08-23 | 2014-08-22 | Method for determining an inflow profile of multilayer reservoir fluids in a wellbore |
Country Status (2)
Country | Link |
---|---|
US (1) | US20150053398A1 (en) |
RU (1) | RU2531499C1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170016315A1 (en) * | 2014-05-02 | 2017-01-19 | Halliburton Energy Services, Inc. | Model for one-dimensional temperature distribution calculations for a fluid in a wellbore |
RU2645692C1 (en) * | 2016-12-21 | 2018-02-27 | Шлюмберже Текнолоджи Б.В. | Method for determining profile of fluid influx in multi-pay well |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030140711A1 (en) * | 2000-03-30 | 2003-07-31 | Brown George A | Method and apparatus for flow measurement |
US20070158064A1 (en) * | 2003-12-24 | 2007-07-12 | Pribnow Daniel F C | Method of determining a fluid inflow profile of wellbore |
US20090272526A1 (en) * | 2008-04-18 | 2009-11-05 | David Booth Burns | Electrical current flow between tunnels for use in heating subsurface hydrocarbon containing formations |
US20110154896A1 (en) * | 2009-12-31 | 2011-06-30 | Schlumberger Technology Corporation | Apparatus for a borehole influx fluid parameters measurement |
WO2011081552A1 (en) * | 2009-12-31 | 2011-07-07 | Шлюмберже Холдингс Лимитед | Method for determining the profile of an inflow and the parameters of a well-surrounding area in a multipay well |
WO2013062446A1 (en) * | 2011-10-26 | 2013-05-02 | Шлюмберже Холдингс Лимитед | Method for determining the inflow profile of fluids of multilayer deposits |
US20140326063A1 (en) * | 2010-02-22 | 2014-11-06 | John R. Lovell | Virtual flowmeter for a well |
US20140365130A1 (en) * | 2011-12-20 | 2014-12-11 | Bp Exploration Operating Company Limited | Estimating flow rates from multiple hydrocarbon reservoir layers into a production well |
US20160177712A1 (en) * | 2014-12-19 | 2016-06-23 | Schlumberger Technology Corporation | Method for determining a water intake profile in an injection well |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1421858A1 (en) * | 1986-11-19 | 1988-09-07 | Среднеазиатский научно-исследовательский институт природного газа | Method and apparatus for determining the profile of fluid inflow in operating gas well |
RU2121572C1 (en) * | 1997-08-13 | 1998-11-10 | Башкирский государственный университет | Method for investigation of injection wells |
RU2154161C1 (en) * | 1999-01-05 | 2000-08-10 | Башкирский государственный университет | Method of injection wells research (versions) |
GB9916022D0 (en) * | 1999-07-09 | 1999-09-08 | Sensor Highway Ltd | Method and apparatus for determining flow rates |
RU2455482C2 (en) * | 2010-09-30 | 2012-07-10 | Шлюмберже Текнолоджи Б.В. | Method of determination of fluid-movement profile and parameters of near-wellbore |
-
2013
- 2013-08-23 RU RU2013139149/03A patent/RU2531499C1/en active
-
2014
- 2014-08-22 US US14/465,978 patent/US20150053398A1/en not_active Abandoned
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030140711A1 (en) * | 2000-03-30 | 2003-07-31 | Brown George A | Method and apparatus for flow measurement |
US20070158064A1 (en) * | 2003-12-24 | 2007-07-12 | Pribnow Daniel F C | Method of determining a fluid inflow profile of wellbore |
US20090272526A1 (en) * | 2008-04-18 | 2009-11-05 | David Booth Burns | Electrical current flow between tunnels for use in heating subsurface hydrocarbon containing formations |
US20110154896A1 (en) * | 2009-12-31 | 2011-06-30 | Schlumberger Technology Corporation | Apparatus for a borehole influx fluid parameters measurement |
WO2011081552A1 (en) * | 2009-12-31 | 2011-07-07 | Шлюмберже Холдингс Лимитед | Method for determining the profile of an inflow and the parameters of a well-surrounding area in a multipay well |
US20130138348A1 (en) * | 2009-12-31 | 2013-05-30 | Schlumberger Technology Corporation | Method for determining the profile of an inflow and the parameters of a well-surrounding area in a multipay well |
US20140326063A1 (en) * | 2010-02-22 | 2014-11-06 | John R. Lovell | Virtual flowmeter for a well |
US20140288836A1 (en) * | 2011-10-23 | 2014-09-25 | Schlumberger Technology Corporation | Method for determining the inflow profile of fluids of multilayer deposits |
WO2013062446A1 (en) * | 2011-10-26 | 2013-05-02 | Шлюмберже Холдингс Лимитед | Method for determining the inflow profile of fluids of multilayer deposits |
US20140365130A1 (en) * | 2011-12-20 | 2014-12-11 | Bp Exploration Operating Company Limited | Estimating flow rates from multiple hydrocarbon reservoir layers into a production well |
US20160177712A1 (en) * | 2014-12-19 | 2016-06-23 | Schlumberger Technology Corporation | Method for determining a water intake profile in an injection well |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170016315A1 (en) * | 2014-05-02 | 2017-01-19 | Halliburton Energy Services, Inc. | Model for one-dimensional temperature distribution calculations for a fluid in a wellbore |
RU2645692C1 (en) * | 2016-12-21 | 2018-02-27 | Шлюмберже Текнолоджи Б.В. | Method for determining profile of fluid influx in multi-pay well |
Also Published As
Publication number | Publication date |
---|---|
RU2531499C1 (en) | 2014-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ramazanov et al. | Thermal modeling for characterization of near wellbore zone and zonal allocation | |
US20150083405A1 (en) | Method of conducting diagnostics on a subterranean formation | |
RU2455482C2 (en) | Method of determination of fluid-movement profile and parameters of near-wellbore | |
US10174612B2 (en) | Method for determining a water intake profile in an injection well | |
Tabatabaei et al. | Fracture-stimulation diagnostics in horizontal wells through use of distributed-temperature-sensing technology | |
US8606523B2 (en) | Method to determine current condensate saturation in a near-wellbore zone in a gas-condensate formation | |
Valiullin et al. | Interpretation of non-isothermal testing data based on the numerical simulation | |
Hoang et al. | Interpretation of wellbore temperatures measured using distributed temperature sensors during hydraulic fracturing | |
Kalia et al. | Wellbore Monitoring in Unconventional Reservoirs: Value of Accurate DTS Interpretation and Risks Involved | |
Sidorova et al. | Do Not Let Temperature Transients Hinder Your Build-up Pressure Interpretation–Proper Gauge Placement in Highly Productive Reservoirs in Well Testing Operations | |
US20150053398A1 (en) | Method for determining an inflow profile of multilayer reservoir fluids in a wellbore | |
EP2772610B1 (en) | Method for determining the inflow profile of fluids of multilayer deposits | |
US20160061025A1 (en) | Method for determining downhole pressure | |
Muradov et al. | Some case studies of temperature and pressure transient analysis in Horizontal, multi-zone, intelligent wells | |
US11236608B2 (en) | Method for injectivity profiling of injection wells | |
McCullagh et al. | Coupling distributed temperature sensing (DTS) based wellbore temperature models with microseismic data for enhanced characterization of hydraulic fracture stimulation | |
Meng et al. | Numerical simulation of fracture temperature field distribution during oil and gas reservoir hydraulic fracturing based on unsteady wellbore temperature field model | |
US20220010672A1 (en) | The method of determining a production well flow profile, including determination of hydrodynamic characteristics of reservoir pay zone | |
RU2645692C1 (en) | Method for determining profile of fluid influx in multi-pay well | |
Chevarunotai | Analytical Models for Flowing-Fluid Temperature Distribution in Single-Phase Oil Reservoirs Accounting for Joule-Thomson Effect | |
Zhang et al. | Flow Profile Determination from Inversion of Distributed Temperature Measurements | |
Maltsev et al. | Evaluating Efficiency of Multilateral Producing Wells in Bottom Water-Drive Reservoir with a Gas Cap by Distributed Fiber-Optic Sensors and Continuous Pressure Monitoring | |
Donovan et al. | Improved Production Profiling Using Thermal Balance and Statistical Modeling in the Pinedale Anticline of the US Rocky Mountains | |
Lavery et al. | Determining Produced Fluid Properties for Accurate Production Profiling During a Drill Stem Test Using Thermal Imaging Technology. | |
US20230194320A1 (en) | Virtual flow rate test |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SCHLUMBERGER TECHNOLOGY CORPORATION, TEXAS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:SHAKO, VALERY VASILIEVICH;PIMENOV, VYACHESLAV PAVLOVICH;THEUVENY, BERTRAND;AND OTHERS;REEL/FRAME:034276/0575 Effective date: 20141119 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO PAY ISSUE FEE |