CN106021958A - Method and device for determining temperatures of upper end and lower end of gas injection well packer - Google Patents

Method and device for determining temperatures of upper end and lower end of gas injection well packer Download PDF

Info

Publication number
CN106021958A
CN106021958A CN201610416823.1A CN201610416823A CN106021958A CN 106021958 A CN106021958 A CN 106021958A CN 201610416823 A CN201610416823 A CN 201610416823A CN 106021958 A CN106021958 A CN 106021958A
Authority
CN
China
Prior art keywords
temperature
packer
instlated tubular
steam
wall
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.)
Granted
Application number
CN201610416823.1A
Other languages
Chinese (zh)
Other versions
CN106021958B (en
Inventor
马振
张洪君
孙振宇
曲绍刚
韦志华
王智博
杨宝春
张成博
李君�
刘锦
王磊
刘洪芹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201610416823.1A priority Critical patent/CN106021958B/en
Publication of CN106021958A publication Critical patent/CN106021958A/en
Application granted granted Critical
Publication of CN106021958B publication Critical patent/CN106021958B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/06Measuring temperature or pressure
    • E21B47/07Temperature
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The invention relates to the field of steam injection thermal recovery of thickened oil, in particular to a method and device for determining the temperatures of the upper end and the lower end of a gas injection well packer. The method includes the steps that the temperature field distribution of steam in a gas injection well heat insulation pipe is acquired; the temperature field distribution of the outer wall of the heat insulation pipe and the temperature field distribution of the inner wall of a casing pipe are acquired according to the temperature field distribution of the steam in the heat insulation pipe; the temperature of the lower end of the packer and the temperature of the upper end of the packer are acquired according to the temperature field distribution of the outer wall of the heat insulation pipe and the temperature field distribution of the inner wall of the casing pipe. According to the method, the temperature of the upper end of the packer and the temperature of the lower end of the packer are calculated, and therefore the position of the packer can be optimized; the method is of great guiding significance on improvement of the heat insulation effect of a gas injection well, prolonging of the thermal recovery life of an oil well and improvement of the quality of steam and the oil well thermal recovery effect.

Description

A kind of determination method and device of gas injection well packer upper and lower two ends temperature
Technical field
The present invention relates to thick oil filling steam oil recovery by heating field, particularly relate to gas injection well packer upper and lower two ends temperature really Determine method and device.
Background technology
Steam injection oil recovery by heating is based on thermodynamics and the theory of thermal conduction study and method, coagulates oil based on viscous crude and height Object to be exploited, reduces Crude viscosity by mode of heating (injection steam), releases oil layer blocking, improves stratum filtration special Property, thus improve the production technique of oil recovery factor.At present, the heat insulation mode of land steam thermal recovery of heavy oil field power DP technology makes mostly By the heat insulation technique of " instlated tubular+packer ", under this technique, owing to the heat enthalpy value of steam is much larger than the air of annular space, because of This some steam enters annular space and causes annular space temperature to raise, thus sleeve pipe can be caused fatigue damage;Packer be for Sealing oil jacket annular space, such that it is able to avoid thermal stress to act directly on sleeve pipe, and the position of packer can direct shadow Ring protected effect.But, current packer installation site is not fixed, but the most arbitrarily lays, thus have impact on envelope Protected effect every device.
Summary of the invention
The embodiment of the present application provides the determination method and device of a kind of gas injection well packer upper and lower two ends temperature, to optimize Packer location, improves well thermal recovery effect.
For reaching above-mentioned purpose, on the one hand, the embodiment of the present application provides gas injection well packer upper and lower two ends temperature really Determining method, described method includes:
The thermo parameters method of steam in acquisition gas injection well instlated tubular;
According to the thermo parameters method of steam in described instlated tubular, obtain thermo parameters method and the internal surface of sleeve pipe of instlated tubular outer wall Thermo parameters method;
Thermo parameters method according to described instlated tubular outer wall and the thermo parameters method of described internal surface of sleeve pipe, obtain under packer The temperature of end and the temperature of packer upper end.
Further, the described temperature of acquisition packer lower end and the temperature of packer upper end include:
Calculate the mean heat loss of packer lower end;
Under thermo parameters method, the thermo parameters method of described internal surface of sleeve pipe and described packer according to described instlated tubular outer wall The mean heat loss of end, calculates the temperature of packer lower end;
Calculate the mean heat loss of packer;
Temperature according to described packer lower end and the mean heat loss of described packer, calculate the temperature of packer upper end Degree.
Further, the temperature of described packer lower end calculates according to the following equation:
TF=TE-Q1*R′/dh
Wherein,
R ′ = 1 2 πh f r c i
In formula: TFTemperature for packer lower end;TEFor the temperature at instlated tubular end outlet;Q1For instlated tubular lower end Mean heat loss;Dh is steam (vapor) outlet to the distance of packer lower end;hfFor moisture film heat transfer coefficient;rciFor internal surface of sleeve pipe radius.
Further, the temperature of described packer upper end calculates according to the following equation:
TG=TF-Q2*R″/dh
Wherein,
R ′ ′ = 1 2 πK f l n d h + L f d h
In formula: TGTemperature for packer upper end;TFTemperature for packer lower end;Q2Mean heat loss for instlated tubular; Dh is steam (vapor) outlet to the distance of Thermal packer lower end;KfFor packer axial thermal conductivity coefficient;LfAxial effective for packer Length.
Further, in described instlated tubular, the thermo parameters method of steam calculates according to the following equation:
Ts=195.94P0.225-17.8
Wherein, TsFor the temperature of steam in instlated tubular, P is the pressure of steam in instlated tubular.
Further, the thermo parameters method of described instlated tubular outer wall calculates according to the following equation:
To=Ts-(R1+R2+R3+R4)dq/dl
Wherein,
R 1 = 1 2 πh f r t i
R 2 = 1 2 πK t u b l n r t o r t i
R 3 = 1 2 πK i n s l n r i r t o
R 4 = 1 2 πK t u b l n r o r i
In formula: T0For instlated tubular outside wall temperature;TsFor the temperature of steam in instlated tubular;Dl is well segment length, and dq is well segment length Heat loss on degree dl;hfFor moisture film heat transfer coefficient;rtiFor insulated tubing inner tube wall radius;KtubFor oil pipe heat conductivity;rto For insulated tubing outer wall of inner tube radius;KinsFor heat insulation layer heat conductivity;riFor insulated tubing outer tube wall radius;roFor insulation oil Pipe outer tube outer wall radius.
Further, the thermo parameters method of described internal surface of sleeve pipe calculates according to the following equation:
Tci=Te+(R6+R7+R8)dq/dl
Wherein,
R 6 = 1 2 πK c a s l n r c o r c i
R 7 = 1 2 πK c e m l n r h r c o
R 8 = f ( t ) 2 πK e
In formula, TciFor instlated tubular outside wall temperature;TeFor formation temperature;Dl is well segment length, and dq is in well segment length dl Heat loss;rciFor internal surface of sleeve pipe radius;rcoFor sleeve outer wall radius;rhFor cement sheath radius;KcasFor sleeve pipe heat conductivity;Kcem For cement sheath heat conductivity;KeFor formation thermal conductivity;F (t) is Ramey time function.
On the other hand, the embodiment of the present application additionally provides the determination device of a kind of gas injection well packer upper and lower two ends temperature, Described device includes:
First computing unit, the thermo parameters method of steam in obtaining gas injection well instlated tubular;
Second computing unit, for according to the thermo parameters method of steam in described instlated tubular, obtaining the temperature of instlated tubular outer wall Degree field distribution and the thermo parameters method of internal surface of sleeve pipe;
3rd computing unit, for the thermo parameters method according to described instlated tubular outer wall and the temperature field of described internal surface of sleeve pipe Distribution, obtains temperature and the temperature of packer upper end of packer lower end.
Further, described 3rd computing unit includes:
Packer lower end heat loss computation subunit, for calculating the mean heat loss of packer lower end;
Packer lower temperature computation subunit, for according to the thermo parameters method of described instlated tubular outer wall, described sleeve pipe The thermo parameters method of inwall and the mean heat loss of described packer lower end, calculate the temperature of packer lower end;
Packer heat loss computation subunit, for calculating the mean heat loss of packer;
Packer upper end temperature computation subelement, flat for according to the temperature of described packer lower end and described packer All heat loss, calculate the temperature of packer upper end.
Further, the temperature of described packer lower end calculates according to the following equation:
TF=TE-Q1*R′/dh
Wherein,
R ′ = 1 2 πh f r c i
In formula: TFTemperature for packer lower end;TEFor the temperature at instlated tubular end outlet;Q1For instlated tubular lower end Mean heat loss;Dh is steam (vapor) outlet to the distance of packer lower end;hfFor moisture film heat transfer coefficient;rciFor internal surface of sleeve pipe radius.
Further, the temperature of described packer upper end calculates according to the following equation:
TG=TF-Q2*R″/dh
Wherein,
R ′ ′ = 1 2 πK f l n d h + L f d h
In formula: TGTemperature for packer upper end;TFTemperature for packer lower end;Q2Mean heat loss for instlated tubular; Dh is steam (vapor) outlet to the distance of Thermal packer lower end;KfFor packer axial thermal conductivity coefficient;LfAxial effective for packer Length.
Further, in described instlated tubular, the thermo parameters method of steam calculates according to the following equation:
Ts=195.94P0.225-17.8
Wherein, TsFor the temperature of steam in instlated tubular, P is the pressure of steam in instlated tubular.
Further, the thermo parameters method of described instlated tubular outer wall calculates according to the following equation:
To=Ts-(R1+R2+R3+R4)dq/dl
Wherein,
R 1 = 1 2 πh f r t i
R 2 = 1 2 πK t u b l n r t o r t i
R 3 = 1 2 πK i n s l n r i r t o
R 4 = 1 2 πK t u b l n r o r i
In formula: T0For instlated tubular outside wall temperature;TsFor the temperature of steam in instlated tubular;Dl is well segment length, and dq is well segment length Heat loss on degree dl;hfFor moisture film heat transfer coefficient;rtiFor insulated tubing inner tube wall radius;KtubFor oil pipe heat conductivity;rto For insulated tubing outer wall of inner tube radius;KinsFor heat insulation layer heat conductivity;riFor insulated tubing outer tube wall radius;roFor insulation oil Pipe outer tube outer wall radius.
Further, the thermo parameters method of described internal surface of sleeve pipe calculates according to the following equation:
Tci=Te+(R6+R7+R8)dq/dl
Wherein,
R 6 = 1 2 πK c a s l n r c o r c i
R 7 = 1 2 πK c e m l n r h r c o
R 8 = f ( t ) 2 πK e
In formula, TciFor instlated tubular outside wall temperature;TeFor formation temperature;Dl is well segment length, and dq is in well segment length dl Heat loss;rciFor internal surface of sleeve pipe radius;rcoFor sleeve outer wall radius;rhFor cement sheath radius;KcasFor sleeve pipe heat conductivity;Kcem For cement sheath heat conductivity;KeFor formation thermal conductivity;F (t) is Ramey time function.
The determination method and device of a kind of gas injection well packer upper and lower two ends temperature that the embodiment of the present application provides, first obtains Take the thermo parameters method of steam in gas injection well instlated tubular, then determine that the thermo parameters method of instlated tubular outer wall and internal surface of sleeve pipe Thermo parameters method, calculates under packer finally according to the thermo parameters method of instlated tubular outer wall and the thermo parameters method of internal surface of sleeve pipe The temperature of end and the temperature of packer upper end.The application calculates temperature and the temperature of packer lower end of packer upper end, from And can according to this guide field operation, optimize packer location, thus to the effect of heat insulation of gas injection well, the thermal recovery life-span of oil well with And the raising of quality of steam and well thermal recovery effect has important directive significance.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this Some embodiments described in application, for those of ordinary skill in the art, in the premise not paying creative work Under, it is also possible to other accompanying drawing is obtained according to these accompanying drawings.
Fig. 1 is the structural representation of the gas injection well pit shaft of the embodiment of the present application;
Fig. 2 is the flow chart of the determination method of the gas injection well packer upper and lower two ends temperature of the embodiment of the present application;
Fig. 3 is the structure chart of the determination device of the gas injection well packer upper and lower two ends temperature of the embodiment of the present application.
Detailed description of the invention
For the technical scheme making those skilled in the art be more fully understood that in the application, real below in conjunction with the application Execute the accompanying drawing in example, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described enforcement Example is only some embodiments of the present application rather than whole embodiments.Based on the embodiment in the application, this area is common The every other embodiment that technical staff is obtained under not making creative work premise, all should belong to the application protection Scope.
Below in conjunction with the accompanying drawings, the detailed description of the invention of the embodiment of the present application is described in further detail.
With reference to Fig. 2, the embodiment of the present application provides the determination method of a kind of gas injection well packer upper and lower two ends temperature, the party Method includes that steps S1 is to step S3.
The thermo parameters method of steam in S1, acquisition gas injection well instlated tubular.
S2, according to the thermo parameters method of steam in described instlated tubular, obtain thermo parameters method and the sleeve pipe of instlated tubular outer wall The thermo parameters method of inwall.
S3, according to the thermo parameters method of described instlated tubular outer wall and the thermo parameters method of described internal surface of sleeve pipe, obtain packing The temperature of device lower end and the temperature of packer upper end.
By said method, the temperature of steam, the temperature of instlated tubular outer wall and the temperature of internal surface of sleeve pipe are considered Etc. factor along the change of well depth, such that it is able to calculate the temperature at the upper and lower two ends of packer accurately, and then packer can be passed through The up and down temperature optimization packer location at two ends, to the effect of heat insulation of steam injection well, the thermal recovery life-span of oil well and quality of steam and The raising of well thermal recovery effect has important directive significance.
Main assumption condition in the embodiment of the present invention is:
(1) flowing for one-dimensional stable and heat transfer in insulated tubing, flow velocity, pressure, temperature are the most axially varying, diametrically Identical;
(2) stratum axial symmetry distribution centered by pit shaft axis, the heat transfer of cement sheath outer rim to stratum is one-dimensional and unsteady state Heat transfer, obeys the non dimensional time function of Ramey, and does not consider the heat transfer along well depth direction;
(3) oil jacket annular space is full of air and part water (will be partially converted into water vapour during steam injection);
(4) ignore the change along well depth direction of formation thermal conductivity and geothermal gradient, and be considered a constant;
(5) packer sealing condition is good, No leakage phenomenon;
(6) packer is compression, main heat-barrier material be politef and graphite meet material, effective length is big It is about 150~200mm.
It is gas injection well structural representation as described in Figure 1, in one embodiment, before step S1, also includes:
Obtain the associated parameter data of gas injection well pit shaft.
Concrete, described associated parameter data, including:
Casing programme and relevant thermophysical property: cement heat conductivity Kcem;Well radius rh;Instlated tubular degree of depth h;Heat insulation Pipe heat conductivity Kins;Instlated tubular inner tube diameter rti;Instlated tubular inner tube external diameter rto;Instlated tubular outer tube external diameter ro;In instlated tubular outer tube Footpath ri;Casing inner diameter rci;Sleeve outer rco;Packer axial thermal conductivity COEFFICIENT Kf;Axial effective length L of packerf;Steam Vapor outlet is to distance dh of packer lower end;Steam injection pressure;Steam injection mass dryness fraction;Steam injection speed;The steam injection time;Geothermal gradient;Earth's surface Temperature;Formation thermal conductivity.
In one embodiment, the thermo parameters method of steam in described acquisition gas injection well instlated tubular, including:
The pressure field distribution of steam in calculating instlated tubular;
According to the pressure field distribution of steam in described instlated tubular, calculate the thermo parameters method of steam in described instlated tubular.
Concrete, the embodiment of the present application, when calculating the field distribution of steam in instlated tubular, can use the side that segmentation calculates Formula, will be divided into the infinitesimal section of several a length of dl by gas injection well pit shaft, may then pass through the side of integration from well head to shaft bottom Method tries to achieve the field distribution of steam in instlated tubular.
Steam in instlated tubular can be considered as steam water two phase flow, and in instlated tubular, the pressure drop of steam water two phase flow is that friction is damaged Mistake, potential variation and the synthesis result of kinetic energy change.The embodiment of the present application can obtain following formula according to momentum balance principle:
Dp=ρmgdl+ρmvmdl-τf
Wherein, p is the pressure of steam, τ in described instlated tubularfFor friction loss gradient, ρmFor the density of multiphase mixture, g For acceleration of gravity, vmFor the flow velocity of multiphase mixture, dl is the infinitesimal section of described instlated tubular.
Three pressure drops of attaching most importance to respectively on the right of equal sign in above formula, accelerate pressure drop and friction pressure drop, permissible when calculating It is positive direction in orientation, may thereby determine that the symbol before described three.
Further, the embodiment of the present application can use Beggs-Brill method, derives soda pop two on the basis of above formula The pressure drop computing formula flowed mutually:
d p d l = ρ m g sin θ + ρ m v m dv m d l - f m ρ m d v m 2 2 - - - ( 1 )
Wherein, fmCoefficient of frictional resistance when flowing for multiphase mixture, d is the caliber of described instlated tubular, and θ is hole angle Complementary angle.Further, above formula refinement can be arranged and obtain calculating the pressure of steam in described instlated tubular by the embodiment of the present application The formula of field distribution:
d p d l = - [ ρ l H l + ρ g ( 1 - H l ) ] g sin θ + f m Gv m 2 d A 1 - [ ρ l H l + ρ g ( 1 - H l ) ] v m v s g p - - - ( 2 )
Wherein, p is the pressure of steam in described instlated tubular, and dl is the infinitesimal section of described instlated tubular, ρlFor density of liquid phase, ρg For density of gas phase, HlFor liquid holdup, A is the sectional area of described instlated tubular, and g is acceleration of gravity, vmStream for multiphase mixture Speed, fmCoefficient of frictional resistance when flowing for multiphase mixture, d is the caliber of described instlated tubular, and θ is the complementary angle of hole angle, and G is Inject the mass flow of steam, vsgFor gas superficial flow velocity.
Calculate in described instlated tubular after the pressure field distribution of steam, can calculate according to saturated vapor temperature and pressure computing formula Go out the thermo parameters method of steam in instlated tubular.Specifically, in described instlated tubular, the thermo parameters method of steam can be according to following public affairs Formula calculates:
TS=195.94P0.225-17.8 (3)
Wherein, TsFor the temperature of steam in instlated tubular, p is the pressure of steam in described instlated tubular.
In one embodiment, the thermo parameters method of described instlated tubular outer wall calculates according to the following equation:
To=Ts-(R1+R2+R3+R4)dq/dl (4)
Wherein,
R 1 = 1 2 πh f r t i - - - ( 5 )
R 2 = 1 2 πK t u b l n r t o r t i - - - ( 6 )
R 3 = 1 2 πK i n s l n r i r t o - - - ( 7 )
R 4 = 1 2 πK t u b l n r o r i - - - ( 8 )
In formula: T0For instlated tubular outside wall temperature;TsFor the temperature of steam in instlated tubular;Dl is well segment length, and dq is well segment length Heat loss on degree dl;hfFor moisture film heat transfer coefficient;rtiFor insulated tubing inner tube wall radius;KtubFor oil pipe heat conductivity;rto For insulated tubing outer wall of inner tube radius;KinsFor heat insulation layer heat conductivity;riFor insulated tubing outer tube wall radius;roFor insulation oil Pipe outer tube outer wall radius.
In one embodiment, the thermo parameters method of described internal surface of sleeve pipe calculates according to the following equation:
Tci=Te+(R6+R7+R8)dq/dl (9)
Wherein,
R 6 = 1 2 πK c a s l n r c o r c i - - - ( 10 )
R 7 = 1 2 πK c e m l n r h r c o - - - ( 11 )
R 8 = f ( t ) 2 πK e - - - ( 12 )
In formula, TciFor instlated tubular outside wall temperature;TeFor formation temperature;Dl is well segment length, and dq is in well segment length dl Heat loss;rciFor internal surface of sleeve pipe radius;rcoFor sleeve outer wall radius;rhFor cement sheath radius;KcasFor sleeve pipe heat conductivity;Kcem For cement sheath heat conductivity;KeFor formation thermal conductivity;F (t) is Ramey time function.At present, time dependent heat conduction passes Heat content f (t) formula has a lot, the most frequently used WHAP model having Ramey model and K.Chiu et al..T worked as by Ramey model > During 11d very accurately, also can use during steam injection only 1d, error only 11%.WHAP model can be used for any steam injection time, but calculates error Up to 15%.Use Ramey model herein.
Concrete, heat loss dq in well segment length dl can be by calculating according to the following equation:
d q d l = 2 πr t o U t o K e K e + r t o U t o f ( t ) ( T s - T e ) - - - ( 13 )
In formula, UtoFor total heat conductivity;KeFor formation thermal conductivity;rtoFor insulated tubing outer wall of inner tube radius;TsFor every Vapor (steam) temperature in heat pipe;TeFor formation temperature;F (t) is Ramey time function.
Total heat conductivity:
U t o = 1.428 × 1 R - - - ( 14 )
Wherein entire thermal resistance:
R=R1+R2+R3+R4+R5+R6+R7+R8 (15)
R 5 = 1 2 π ( h c + h r ) r o - - - ( 16 )
In formula: R is entire thermal resistance, roFor insulated tubing outer tube outer wall radius;rciFor internal surface of sleeve pipe radius;rcoFor sleeve outer wall Radius;rhFor cement sheath radius;hrFor oil jacket annular space radiation heat transfer coefficient;hcFor oil jacket annular space free convection heat transfer coefficient.
In one embodiment, the described thermo parameters method according to described instlated tubular outer wall and the temperature of described internal surface of sleeve pipe Field distribution, obtains temperature and the temperature of packer upper end of packer lower end, including:
Calculate the mean heat loss of packer lower end;
Under thermo parameters method, the thermo parameters method of described internal surface of sleeve pipe and described packer according to described instlated tubular outer wall The mean heat loss of end, calculates the temperature of packer lower end;
Calculate the mean heat loss of packer;
Temperature according to described packer lower end and the mean heat loss of described packer, calculate the temperature of packer upper end Degree.
Concrete, can be according to below equation calculating pit shaft radially heat loss:
Q = T s - T e R d l - - - ( 17 )
In formula, Q is pit shaft radially heat loss;TsFor vapor (steam) temperature in instlated tubular;TsFor formation temperature;R is entire thermal resistance.
As it is shown in figure 1, calculate radial direction heat loss Q at interface A C by formula (17)ACWith the radial direction heat waste at the BD of interface QBD, then the mean heat loss Q of packer lower end1Can calculate according to below equation:
Q 1 = Q A C + Q B D 2 - - - ( 18 )
In one embodiment, as it is shown in figure 1, A point internal surface of sleeve pipe temperature can be calculated by formula (4), (9) TA, B point internal surface of sleeve pipe temperature TB, M point internal surface of sleeve pipe temperature TM, C point instlated tubular outside wall temperature TC, D point instlated tubular outside wall temperature TD, N point instlated tubular outside wall temperature TN
In one embodiment, the temperature of described packer lower end can calculate according to the following equation:
TF=TE-Q1*R′/dh (19)
Wherein,
R ′ = 1 2 πh f r c i - - - ( 20 )
T E = T A + T C 2 - - - ( 21 )
In formula: TFTemperature for packer lower end;TEFor the temperature at instlated tubular end outlet;Q1For instlated tubular lower end Mean heat loss;Dh is steam (vapor) outlet to the distance of packer lower end;hfFor moisture film heat transfer coefficient;rciFor internal surface of sleeve pipe radius.
In one embodiment, radial direction heat loss Q at the MN of interface can be calculated by formula (17)MNAt the BD of interface Radially heat loss QBD, then the mean heat loss Q of packer2Can calculate according to below equation:
Q 2 = Q M N + Q B D 2 - - - ( 22 )
The temperature of described packer upper end calculates according to the following equation:
TG=TF-Q2*R″/dh (23)
Wherein,
R ′ ′ = 1 2 πK f l n d h + L f d h - - - ( 24 )
In formula: TGTemperature for packer upper end;TFTemperature for packer lower end;Q2Mean heat loss for instlated tubular; Dh is steam (vapor) outlet to the distance of packer lower end, i.e. dh=AB=CD=EF;KfFor packer axial thermal conductivity coefficient;LfFor envelope Axial effective length every device.
In one embodiment, after step s 3, also include:
The temperature of described packer lower end and the temperature of packer upper end are modified, it is thus achieved that packer lower end after correction Temperature and the temperature of packer upper end.
Concrete, according to below equation, the temperature of described packer lower end can be modified:
T F ′ = T F + T B + T D 3 - - - ( 25 )
In formula, TF' be revise after the temperature of packer lower end.
Specifically, according to below equation, the temperature of described packer upper end can be modified:
T G ′ = T G + T M + T N 3 - - - ( 26 )
The embodiment of the present application is by calculating temperature and the temperature of packer lower end of packer upper end, when temperature does not meets pre- If during condition, packer location can be adjusted, calculate the temperature of packer upper and lower side the most again, through meet pre-conditioned be Only, such that it is able to guide field operation according to this, it is to avoid seal because of too high and sleeve pipe is caused damage every device upper and lower side temperature;Lead to simultaneously Cross packer upper and lower side temperature optimization packer location, to the effect of heat insulation of gas injection well, the thermal recovery life-span of oil well and steam matter The raising of amount and well thermal recovery effect has important directive significance.
As it is shown on figure 3, the embodiment of the present application additionally provides the determination device of a kind of gas injection well packer upper and lower two ends temperature, Described device includes:
First computing unit 21, the thermo parameters method of steam in obtaining gas injection well instlated tubular;
Second computing unit 22, for according to the thermo parameters method of steam in described instlated tubular, obtaining instlated tubular outer wall Thermo parameters method and the thermo parameters method of internal surface of sleeve pipe;
3rd computing unit 23, for the thermo parameters method according to described instlated tubular outer wall and the temperature of described internal surface of sleeve pipe Field distribution, obtains temperature and the temperature of packer upper end of packer lower end.
In one embodiment, described 3rd computing unit includes:
Packer lower end heat loss computation subunit, for calculating the mean heat loss of packer lower end;
Packer lower temperature computation subunit, for according to the thermo parameters method of described instlated tubular outer wall, described sleeve pipe The thermo parameters method of inwall and the mean heat loss of described packer lower end, calculate the temperature of packer lower end;
Packer heat loss computation subunit, for calculating the mean heat loss of packer;
Packer upper end temperature computation subelement, flat for according to the temperature of described packer lower end and described packer All heat loss, calculate the temperature of packer upper end.
Each ingredient of the device of the present embodiment is respectively used to realize each step of the method for previous embodiment, due to In embodiment of the method, each step is described in detail, has not repeated them here.
The determination method and device of a kind of gas injection well packer upper and lower two ends temperature that the embodiment of the present application provides, first counts Calculate the thermo parameters method of steam in gas injection well instlated tubular, then determine that the thermo parameters method of instlated tubular outer wall and internal surface of sleeve pipe Thermo parameters method, calculates under packer finally according to the thermo parameters method of instlated tubular outer wall and the thermo parameters method of internal surface of sleeve pipe The temperature of end and the temperature of packer upper end.The application calculates temperature and the temperature of packer lower end of packer upper end, from And can according to this guide field operation, optimize packer location, thus to the effect of heat insulation of gas injection well, the thermal recovery life-span of oil well with And the raising of quality of steam and well thermal recovery effect has important directive significance.
In one or more exemplary designs, the above-mentioned functions described by the embodiment of the present application can be at hardware, soft The combination in any of part, firmware or this three realizes.If realized in software, these functions can store and computer-readable On medium, or it is transmitted on the medium of computer-readable with one or more instructions or code form.Computer readable medium includes electricity Brain stores medium and is easy to so that allowing computer program transfer to the telecommunication media in other place from a place.Storage medium is permissible It is that any general or special computer can be with the useable medium of access.Such as, such computer readable media can include but It is not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage device, or other What may be used for carrying or storage can be by general or special computer or general or special handling with other with instruction or data structure Device reads the medium of the program code of form.
Particular embodiments described above, has been carried out the most in detail purpose, technical scheme and the beneficial effect of the application Describe in detail bright, be it should be understood that the specific embodiment that the foregoing is only the embodiment of the present application, be not used to limit this Shen Protection domain please, all within spirit herein and principle, any modification, equivalent substitution and improvement etc. done, all should wrap Within being contained in the protection domain of the application.

Claims (14)

1. the determination method of a gas injection well packer upper and lower two ends temperature, it is characterised in that described method includes:
The thermo parameters method of steam in acquisition gas injection well instlated tubular;
According to the thermo parameters method of steam in described instlated tubular, obtain thermo parameters method and the temperature of internal surface of sleeve pipe of instlated tubular outer wall Degree field distribution;
Thermo parameters method according to described instlated tubular outer wall and the thermo parameters method of described internal surface of sleeve pipe, obtain packer lower end Temperature and the temperature of packer upper end.
2. the method for claim 1, it is characterised in that the temperature of described acquisition packer lower end and packer upper end Temperature includes:
Calculate the mean heat loss of packer lower end;
Thermo parameters method, the thermo parameters method of described internal surface of sleeve pipe and described packer lower end according to described instlated tubular outer wall Mean heat loss, calculates the temperature of packer lower end;
Calculate the mean heat loss of packer;
Temperature according to described packer lower end and the mean heat loss of described packer, calculate the temperature of packer upper end.
3. method as claimed in claim 2, it is characterised in that the temperature of described packer lower end calculates according to the following equation:
TF=TE-Q1*R′/dh
Wherein,
R ′ = 1 2 πh f r c i
In formula: TFTemperature for packer lower end;TEFor the temperature at instlated tubular end outlet;Q1Average for instlated tubular lower end Heat loss;Dh is steam (vapor) outlet to the distance of packer lower end;hfFor moisture film heat transfer coefficient;rciFor internal surface of sleeve pipe radius.
4. method as claimed in claim 3, it is characterised in that the temperature of described packer upper end calculates according to the following equation:
TG=TF-Q2*R″/dh
Wherein,
R ′ ′ = 1 2 πK f l n d h + L f d h
In formula: TGTemperature for packer upper end;TFTemperature for packer lower end;Q2Mean heat loss for instlated tubular;Dh is Steam (vapor) outlet is to the distance of Thermal packer lower end;KfFor packer axial thermal conductivity coefficient;LfAxial effective for packer is long Degree.
5. the method for claim 1, it is characterised in that in described instlated tubular, the thermo parameters method of steam is according to following public affairs Formula calculates:
Ts=195.94P0.225-17.8
Wherein, TsFor the temperature of steam in instlated tubular, P is the pressure of steam in instlated tubular.
6. the method for claim 1, it is characterised in that the thermo parameters method of described instlated tubular outer wall is according to the following equation Calculate:
To=Ts-(R1+R2+R3+R4)dq/dl
Wherein,
R 1 = 1 2 πh f r t i
R 2 = 1 2 πK t u b ln r t o r t i
R 3 = 1 2 πK i n s l n r i r t o
R 4 = 1 2 πK t u b l n r o r i
In formula: T0For instlated tubular outside wall temperature;TsFor the temperature of steam in instlated tubular;Dl is well segment length, and dq is well segment length dl On heat loss;hfFor moisture film heat transfer coefficient;rtiFor insulated tubing inner tube wall radius;KtubFor oil pipe heat conductivity;rtoFor every Deep fat pipe outer wall of inner tube radius;KinsFor heat insulation layer heat conductivity;riFor insulated tubing outer tube wall radius;roOutside for insulated tubing Pipe exterior radius.
7. the method for claim 1, it is characterised in that the thermo parameters method of described internal surface of sleeve pipe is counted according to the following equation Calculate:
Tci=Te+(R6+R7+R8)dq/dl
Wherein,
R 6 = 1 2 πK c a s l n r c o r c i
R 7 = 1 2 πK c e m l n r h r c o
R 8 = f ( t ) 2 πK e
In formula, TciFor instlated tubular outside wall temperature;TeFor formation temperature;Dl is well segment length, and dq is the heat waste in well segment length dl Lose;rciFor internal surface of sleeve pipe radius;rcoFor sleeve outer wall radius;rhFor cement sheath radius;KcasFor sleeve pipe heat conductivity;KcemFor water Mud ring heat conductivity;KeFor formation thermal conductivity;F (t) is Ramey time function.
8. the determination device of a gas injection well packer upper and lower two ends temperature, it is characterised in that described device includes:
First computing unit, the thermo parameters method of steam in obtaining gas injection well instlated tubular;
Second computing unit, for according to the thermo parameters method of steam in described instlated tubular, obtaining the temperature field of instlated tubular outer wall Distribution and the thermo parameters method of internal surface of sleeve pipe;
3rd computing unit, for dividing according to the thermo parameters method of described instlated tubular outer wall and the temperature field of described internal surface of sleeve pipe Cloth, obtains temperature and the temperature of packer upper end of packer lower end.
9. device as claimed in claim 8, it is characterised in that described 3rd computing unit includes:
Packer lower end heat loss computation subunit, for calculating the mean heat loss of packer lower end;
Packer lower temperature computation subunit, for according to the thermo parameters method of described instlated tubular outer wall, described internal surface of sleeve pipe Thermo parameters method and the mean heat loss of described packer lower end, calculate packer lower end temperature;
Packer heat loss computation subunit, for calculating the mean heat loss of packer;
Packer upper end temperature computation subelement, for the temperature according to described packer lower end and the evenly heat of described packer Loss, calculates the temperature of packer upper end.
10. device as claimed in claim 9, it is characterised in that the temperature of described packer lower end calculates according to the following equation:
TF=TE-Q1*R′/dh
Wherein,
R ′ = 1 2 π h f γ ci
In formula: TFTemperature for packer lower end;TEFor the temperature at instlated tubular end outlet;Q1Average for instlated tubular lower end Heat loss;Dh is steam (vapor) outlet to the distance of packer lower end;hfFor moisture film heat transfer coefficient;rciFor internal surface of sleeve pipe radius.
11. devices as claimed in claim 10, it is characterised in that the temperature of described packer upper end is counted according to the following equation Calculate:
TG=TF-Q2*R″/dh
Wherein,
R ′ ′ = 1 2 πK f l n d h + L f d h
In formula: TGTemperature for packer upper end;TFTemperature for packer lower end;Q2Mean heat loss for instlated tubular;Dh is Steam (vapor) outlet is to the distance of Thermal packer lower end;KfFor packer axial thermal conductivity coefficient;LfAxial effective for packer is long Degree.
12. devices as claimed in claim 8, it is characterised in that in described instlated tubular, the thermo parameters method of steam is according to following Formula calculates:
Ts=195.94P0.225-17.8
Wherein, TsFor the temperature of steam in instlated tubular, P is the pressure of steam in instlated tubular.
13. devices as claimed in claim 8, it is characterised in that the thermo parameters method of described instlated tubular outer wall is according to following public affairs Formula calculates:
To=Ts-(R1+R2+R3+R4)dq/dl
Wherein,
R 1 = 1 2 πh f r t i
R 2 = 1 2 πK t u b l n r t o r t i
R 3 = 1 2 πK i n s l n r i r t o
R 4 = 1 2 πK t u b l n r o r i
In formula: T0For instlated tubular outside wall temperature;TsFor the temperature of steam in instlated tubular;Dl is well segment length, and dq is well segment length dl On heat loss;hfFor moisture film heat transfer coefficient;rtiFor insulated tubing inner tube wall radius;KtubFor oil pipe heat conductivity;rtoFor every Deep fat pipe outer wall of inner tube radius;KinsFor heat insulation layer heat conductivity;riFor insulated tubing outer tube wall radius;roOutside for insulated tubing Pipe exterior radius.
14. devices as claimed in claim 8, it is characterised in that the thermo parameters method of described internal surface of sleeve pipe is according to the following equation Calculate:
Tci=Te+(R6+R7+R8)dq/dl
Wherein,
R 6 = 1 2 πK c a s l n r c o r c i
R 7 = 1 2 πK c e m l n r h r c o
R 8 = f ( t ) 2 πK e
In formula, TciFor instlated tubular outside wall temperature;TeFor formation temperature;Dl is well segment length, and dq is the heat waste in well segment length dl Lose;rciFor internal surface of sleeve pipe radius;rcoFor sleeve outer wall radius;rhFor cement sheath radius;KcasFor sleeve pipe heat conductivity;KcemFor water Mud ring heat conductivity;KeFor formation thermal conductivity;F (t) is Ramey time function.
CN201610416823.1A 2016-06-14 2016-06-14 A kind of determination method and device of gas injection well packer upper and lower ends temperature Active CN106021958B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610416823.1A CN106021958B (en) 2016-06-14 2016-06-14 A kind of determination method and device of gas injection well packer upper and lower ends temperature

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610416823.1A CN106021958B (en) 2016-06-14 2016-06-14 A kind of determination method and device of gas injection well packer upper and lower ends temperature

Publications (2)

Publication Number Publication Date
CN106021958A true CN106021958A (en) 2016-10-12
CN106021958B CN106021958B (en) 2019-04-12

Family

ID=57087922

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610416823.1A Active CN106021958B (en) 2016-06-14 2016-06-14 A kind of determination method and device of gas injection well packer upper and lower ends temperature

Country Status (1)

Country Link
CN (1) CN106021958B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112796704A (en) * 2019-10-28 2021-05-14 中国石油化工股份有限公司 Optimization and parameter optimization method for oil well hot washing paraffin removal mode
CN114059957A (en) * 2020-08-05 2022-02-18 中石化石油工程技术服务有限公司 Method for improving sealing performance of cement sheath at overlapped section of casing

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103310123A (en) * 2013-07-10 2013-09-18 华东理工大学 Coupling modeling method for verifying and optimizing design of industrial ethylene steam cracking furnace
CN102901748B (en) * 2012-09-18 2014-10-22 成都思驰科技有限公司 Nondestructive testing device and method based on pipeline temperature field distribution
CN104392092A (en) * 2014-10-10 2015-03-04 中国石油天然气股份有限公司 Temperature calculation and control method and device for mixed liquor of gravity fire flooding production well
CN104389589A (en) * 2014-10-11 2015-03-04 中国石油天然气股份有限公司 Method and system of determining shaft temperature field distribution based on hollow-rod
CN104806230A (en) * 2015-02-16 2015-07-29 中国石油天然气股份有限公司 Supercritical steam injection well bore temperature field and pressure field distribution calculation method
CN104820748A (en) * 2015-05-07 2015-08-05 北京宇航系统工程研究所 Method for determining temperature field distribution of cabin of carrier rocket in flight phase in atmospheric layer
CN105422084A (en) * 2015-11-26 2016-03-23 中国石油天然气股份有限公司 Wellbore temperature field obtaining method and device for hot water circulating heating viscosity reduction process

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102901748B (en) * 2012-09-18 2014-10-22 成都思驰科技有限公司 Nondestructive testing device and method based on pipeline temperature field distribution
CN103310123A (en) * 2013-07-10 2013-09-18 华东理工大学 Coupling modeling method for verifying and optimizing design of industrial ethylene steam cracking furnace
CN104392092A (en) * 2014-10-10 2015-03-04 中国石油天然气股份有限公司 Temperature calculation and control method and device for mixed liquor of gravity fire flooding production well
CN104389589A (en) * 2014-10-11 2015-03-04 中国石油天然气股份有限公司 Method and system of determining shaft temperature field distribution based on hollow-rod
CN104806230A (en) * 2015-02-16 2015-07-29 中国石油天然气股份有限公司 Supercritical steam injection well bore temperature field and pressure field distribution calculation method
CN104820748A (en) * 2015-05-07 2015-08-05 北京宇航系统工程研究所 Method for determining temperature field distribution of cabin of carrier rocket in flight phase in atmospheric layer
CN105422084A (en) * 2015-11-26 2016-03-23 中国石油天然气股份有限公司 Wellbore temperature field obtaining method and device for hot water circulating heating viscosity reduction process

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
万仁溥主编,: "《采油工程手册(精要本) 2009年8月第1版》", 31 August 2009, 石油工业出版社 *
刘喜林主编: "《难动用储量开发稠油开采技术2005年11月第1版》", 30 November 2005, 石油工业出版社, *
张智,等;: "多封隔器密闭环空热膨胀力学计算方法及应用", 《钻井工程》 *
张毅,等;: "注汽工艺管柱对热采井套损的影响", 《石油机械》 *
杨雄文,等;: "预测井筒流动温度分布的新方法", 《大庆石油地质与开发》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112796704A (en) * 2019-10-28 2021-05-14 中国石油化工股份有限公司 Optimization and parameter optimization method for oil well hot washing paraffin removal mode
CN112796704B (en) * 2019-10-28 2023-07-28 中国石油化工股份有限公司 Optimization method for oil well hot washing paraffin removal mode and parameter optimization method
CN114059957A (en) * 2020-08-05 2022-02-18 中石化石油工程技术服务有限公司 Method for improving sealing performance of cement sheath at overlapped section of casing
CN114059957B (en) * 2020-08-05 2024-04-26 中石化石油工程技术服务有限公司 Method for improving sealing performance of cement sheath at overlapping section of sleeve

Also Published As

Publication number Publication date
CN106021958B (en) 2019-04-12

Similar Documents

Publication Publication Date Title
CN104806230B (en) The Wellbore Temperature Field of overcritical steam injection well and the computational methods of pressure field distribution
CN103775058B (en) A kind of determination method of wellbore heat loss
CN103590812B (en) The method of calculation of a kind of gas well plug-ging amount, calculating device and defining method
CN102011573B (en) Method for uniformly injecting multi-component thermal fluid in horizontal well
CN103899295A (en) Method and device for detecting heat losses of plumb shaft barrel in thickened oil steam production and injection process
CN101684727A (en) Optimization method for determining proportion of mixing light oil into heavy oil of ultra-deep well and light hydrocarbon mixer thereof
CN107145705A (en) A kind of method and device for obtaining circulating temperature
CN105134179B (en) A kind of computational methods of natural gas well wellbore pressure and the distributed data of temperature
Boban et al. Application and design aspects of ground heat exchangers
CN106021958A (en) Method and device for determining temperatures of upper end and lower end of gas injection well packer
CN104481482A (en) Concentric double-pipe gas injection and heat insulation analyzing method and device for horizontal well
CN103900737A (en) Method and device for detecting heat losses of surface pipeline in thickened oil steam production and injection
CN104462750B (en) Steam quality determines method and device
CN102682195A (en) Mine shaft temperature calculation method for semi-submersible type platform transient drilling well
Kang et al. Heat transfer modeling of wellbore with vacuum-insulated-tubing (VIT)
Least et al. Steam flow tests for comparing performance of nozzle, tube, and fluidic diode autonomous ICDs in SAGD wells
CN204536240U (en) A kind of experimental provision measuring Oil/gas Well wellbore fluids, oil well pipe thermodynamic parameter
Pei et al. Heat transfer analysis of underground U-type heat exchanger of ground source heat pump system
CN102108851A (en) Horizontal well steam injection method and system
CN105184061B (en) The method for numerical simulation of gas-producing well temperature, pressure distribution
Jia et al. Experimental and numerical study on the self-balancing heating performance of a thermosyphon during the process of oil production
CN106680177B (en) Scene determines the inner tube drawdown double-tube method of inner tube aquiclude hydrogeological parameter
Koenig et al. Development of a thermal resistance model to evaluate wellbore heat exchange efficiency
CN105952443A (en) Determination method and device for wellbore temperature distribution of flue gas flooding gas-injection well
CN105114060A (en) Method and device for calculating wellbore temperature distribution of gas injection well

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant