CN104389589B - Method and system of determining shaft temperature field distribution based on hollow-rod - Google Patents
Method and system of determining shaft temperature field distribution based on hollow-rod Download PDFInfo
- Publication number
- CN104389589B CN104389589B CN201410534832.1A CN201410534832A CN104389589B CN 104389589 B CN104389589 B CN 104389589B CN 201410534832 A CN201410534832 A CN 201410534832A CN 104389589 B CN104389589 B CN 104389589B
- Authority
- CN
- China
- Prior art keywords
- thermal
- thermal resistance
- liquid
- radius
- oil
- 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.)
- Active
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000007788 liquid Substances 0.000 claims abstract description 212
- 239000003921 oil Substances 0.000 claims abstract description 187
- 239000012530 fluid Substances 0.000 claims abstract description 71
- 244000273618 Sphenoclea zeylanica Species 0.000 claims description 161
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 64
- 238000012546 transfer Methods 0.000 claims description 53
- 239000004568 cement Substances 0.000 claims description 51
- 239000010779 crude oil Substances 0.000 claims description 25
- 230000015572 biosynthetic process Effects 0.000 claims description 21
- 238000004519 manufacturing process Methods 0.000 claims description 18
- 239000003129 oil well Substances 0.000 claims description 17
- 239000000203 mixture Substances 0.000 claims description 14
- 230000005855 radiation Effects 0.000 claims description 11
- 239000000463 material Substances 0.000 claims description 3
- 239000000295 fuel oil Substances 0.000 abstract 2
- 238000005065 mining Methods 0.000 abstract 1
- 241000209094 Oryza Species 0.000 description 180
- 235000007164 Oryza sativa Nutrition 0.000 description 180
- 235000009566 rice Nutrition 0.000 description 180
- 239000013043 chemical agent Substances 0.000 description 13
- 238000010586 diagram Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000010438 heat treatment Methods 0.000 description 4
- 238000005086 pumping Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 230000008014 freezing Effects 0.000 description 2
- 238000007710 freezing Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000012188 paraffin wax Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000001993 wax Substances 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- 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/06—Measuring temperature or pressure
- E21B47/07—Temperature
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)
- Earth Drilling (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention provides a method and a system of determining shaft temperature field distribution by a hollow-rod, the method including the steps of: 1) obtaining all the data relevant to the hollow-rod and a shaft; 2) setting the step size according to the height of the shaft working fluid level and the shaft depth; 3) dividing the shaft and the hollow-rod into a plurality of shaft sections and hollow-rod sections according to the specified step size; 4) determining the temperature of fluid in the shaft sections and the hollow-rod sections according to the specified data; 5) forming the shaft temperature field distribution by combination of temperature of the fluid in the hollow-rod sections and the temperature of the liquid in the shaft sections. The method and the system can determine the shaft temperature field distribution and provide data basis for the reasonable follow-up determination of mixing discharge capacity at a shaft entrance and temperature so as to satisfy the requirement for conventional mining of thickened oil, extra heavy oil and super heavy oil.
Description
Technical field
The present invention with regard to exploration of oil and gas field technical field, especially with regard to the lifting skill of viscous crude, special thick oil and super-viscous oil
Art, is concretely a kind of method and system determining temperature in wellbore field distribution based on hollow stem.
Background technology
In the recovery process of the oil reservoir low in strata pressure, poor permeability, crude oil freezing point are high, paraffin content is high, viscosity is high,
Oil well paraffinication is serious, easily causes starting difficulty after pump leakage or termination of pumping, or even causes holddown accident, brings not to oil well production
Profit impact.
For solving the holddown phenomenon that causes because of viscous crude and pump barrel wax deposition, in prior art using than wide be hollow taking out
Beam hanger is mixed with hot water or the DP technology of mixing light oil or chemical agent.This technological process is：Mixed using existing pumping well ground
Hot water piping is connected with hollow stem, and the hot water of incorporation, through hollow stem endoporus vertical current to bottom water nozzle, is then ejected into open tubular column
Mix with crude oil in plug, mixed liquor enters the annular space of hollow stem and oil pipe, and will be given rise to ground.This technique has comprehensive
Close low cost, easy to maintenance, the advantages of efficiency of utilization is high.
During actually used, hollow rod is mixed with hot water or one of the key of oil recovery technique of mixing light oil or chemical agent
It is to determine that hot water or thin oil or the well head of chemical agent mix discharge capacity and temperature.This two parameter index is mainly subject to viscosity of crude with temperature
The impact of degree change.Therefore, to determine that rational well head mixes discharge capacity and temperature, be necessary for studying the temperature in the wellbore of crude oil
The distribution of field.
Therefore, how to determine the distribution in crude oil temperature field in the wellbore, and then select rational well head to mix row accordingly
Amount and temperature are this area technical barriers urgently to be resolved hurrily to meet the exploitation of existing viscous crude, special thick oil and super-viscous oil.
Content of the invention
In order to solve hollow stem electric-heating technology of the prior art due to Wellbore Temperature Field cannot be determined, and then it is difficult to
Select rational well head to mix discharge capacity and the exploitation that cannot meet existing viscous crude, special thick oil and super-viscous oil that temperature causes difficulty
Topic, the invention provides a kind of method and system determining temperature in wellbore field distribution based on hollow stem, is that one kind is accurately based on
Hollow stem determines the scheme of temperature in wellbore field distribution, by obtaining the data information related to hollow stem and pit shaft, according to setting
Pit shaft, hollow stem are divided into multiple wellbore sections, hollow bar segment by fixed step size, determine the temperature of liquid, sky in each wellbore section successively
The temperature of fluid in core bar section, has so then obtained temperature in wellbore field distribution, be subsequently selected rational well head mix discharge capacity and
Temperature provides data foundation with the exploitation meeting existing viscous crude, special thick oil and super-viscous oil.
It is an object of the invention to provide a kind of method that temperature in wellbore field distribution is determined based on hollow stem, including：Obtain
Take the data information related to hollow stem and pit shaft；Depth-set according to the height in hydrodynamic face and pit shaft in described pit shaft
Step-length；Described pit shaft and hollow stem are divided into by multiple wellbore sections, hollow bar segment according to described step-length；According to described number
Determine the temperature of fluid in the temperature of liquid in the plurality of wellbore section, multiple hollow bar segment according to data respectively；The plurality of sky
The temperature composition temperature in wellbore field distribution of liquid in the temperature of fluid and described multiple wellbore sections in core bar section.
It is an object of the invention to provide a kind of system determining temperature in wellbore field distribution based on hollow stem, including：
Data information acquisition device, for obtaining the data information related to hollow stem and pit shaft；Step size settings device, for basis
The depth-set step-length of the height in hydrodynamic face and pit shaft in described pit shaft；Sectioning, for according to described step-length by institute
The pit shaft stated and hollow stem are divided into multiple wellbore sections, hollow bar segment；Temperature determines device, for according to described data information
Determine the temperature of fluid in the temperature of liquid in the plurality of wellbore section, multiple hollow bar segment respectively；Thermo parameters method determines dress
Put, the temperature composition pit shaft temperature of liquid in the temperature for fluid in the plurality of hollow bar segment and described multiple wellbore sections
Degree field distribution.
The beneficial effects of the present invention is, there is provided a kind of the method for temperature in wellbore field distribution determined based on hollow stem and is
System, is a kind of accurate scheme determining temperature in wellbore field distribution based on hollow stem, by obtaining and hollow stem and pit shaft phase
Pit shaft, hollow stem are divided into multiple wellbore sections, hollow bar segment according to setting step-length, determine each well successively by the data information closing
In the temperature of liquid in cylinder section, hollow bar segment, the temperature of fluid, has so then obtained temperature in wellbore field distribution, has been subsequently selected conjunction
The well head of reason mixes discharge capacity and temperature provides data foundation with the exploitation meeting existing viscous crude, special thick oil and super-viscous oil, enters
And improve the efficiency of the exploitation of viscous crude, special thick oil and super-viscous oil.
It is that the above and other objects, features and advantages of the present invention can be become apparent, preferred embodiment cited below particularly,
And coordinate institute's accompanying drawings, it is described in detail below.
Brief description
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
Have technology description in required use accompanying drawing be briefly described it should be apparent that, drawings in the following description be only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, acceptable
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is a kind of flow process of method determining temperature in wellbore field distribution based on hollow stem provided in an embodiment of the present invention
Figure；
Fig. 2 is the particular flow sheet of step S104 in Fig. 1；
Fig. 3 is the particular flow sheet of step S201 in Fig. 2；
Fig. 4 is the particular flow sheet of step S304 in Fig. 3；
Fig. 5 is the particular flow sheet of step S305 in Fig. 3；
Fig. 6 is the particular flow sheet of the embodiment one of step S202 in Fig. 2；
Fig. 7 is the particular flow sheet of the embodiment two of step S202 in Fig. 2；
Fig. 8 is a kind of structural frames of system determining temperature in wellbore field distribution based on hollow stem provided in an embodiment of the present invention
Figure；
Fig. 9 is a kind of temperature being determined based on hollow stem in the system of temperature in wellbore field distribution provided in an embodiment of the present invention
Determine the concrete structure block diagram of device 104；
Figure 10 is a kind of heat being determined based on hollow stem in the system of temperature in wellbore field distribution provided in an embodiment of the present invention
The concrete structure block diagram of resistance determining module 201；
Figure 11 be provided in an embodiment of the present invention a kind of determined based on hollow stem in the system of temperature in wellbore field distribution
The concrete structure block diagram of one thermal resistance determining unit 304；
Figure 12 be provided in an embodiment of the present invention a kind of determined based on hollow stem in the system of temperature in wellbore field distribution
The concrete structure block diagram of two thermal resistance determining units 305；
Figure 13 is a kind of heat being determined based on hollow stem in the system of temperature in wellbore field distribution provided in an embodiment of the present invention
The concrete structure block diagram of the embodiment one of resistance coefficient determination module 202；
Figure 14 is a kind of heat being determined based on hollow stem in the system of temperature in wellbore field distribution provided in an embodiment of the present invention
The concrete structure block diagram of the embodiment two of resistance coefficient determination module 202；
Figure 15 is rod-pumped well hollow stem electrical heating process structural representation of the prior art.
Specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation description is it is clear that described embodiment is only a part of embodiment of the present invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of not making creative work
Embodiment, broadly falls into the scope of protection of the invention.
The present invention is directed to the exploitation of viscous crude, special thick oil and super-viscous oil it is proposed that a kind of hollow stem is mixed with hot water or mixing light oil
Or mix the computational methods that chemical agent calculates temperature in wellbore field distribution, and define a set of interpretative tool.By to after mixed with hot water
The decoupled method analysis of Wellbore Temperature Field, can calculate every section of viscosity of crude, further can calculate the friction of rod liquid
Power and pipe liquid frrction load, may finally calculate oil pumping machine polished rod load by roofbolt up stroke.The well head that hot water can be adjusted is mixed
Enter discharge capacity and temperature, reach the purpose optimizing hoisting system efficiency.
The basic assumption condition of the present invention includes：
(1), ignore the heat exchange in sucker rod, pit shaft and formation rock longitudinal direction；
(2), the pressure of well head production fluid, temperature keep constant；
(3), the annular space of oil pipe and sleeve pipe formation is full of low-pressure air；
(4), with sucker rod center line as symmetry axis, sucker rod, pit shaft and formation rock isotropism；
(5), the thermal physical property parameter in model system is temperature independent, that is, be considered permanent physical property；
(6), prime stratum temperature is linearly distributed；
(7), hot water is from the injection of hollow stem center.
Fig. 1 is a kind of particular flow sheet of method determining temperature in wellbore field distribution based on hollow stem proposed by the present invention,
As shown in Figure 1, described method includes：
S101：Obtain the data information related to hollow stem and pit shaft.
In the particular embodiment, the data information related to hollow stem and pit shaft includes formation thermal conductivity, stratum
Average coefficient of heat transfer, the oil well production time, wellbore radius, sleeve outer wall radius, cement sheath thermal conductivity factor, sleeve pipe thermal conductivity factor,
Internal surface of sleeve pipe radius, sleeve outer wall radius, annular space radiation heat transfer coefficient, annular space free convection heat transfer coefficient, oil-pipe external wall radius,
Oil pipe thermal conductivity factor, tube inner wall radius, crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem outer wall half
Footpath, relative oil density, hollow stem inwall radius, the thermal conductivity factor of hollow stem.
S102：Depth-set step-length according to the height in hydrodynamic face and pit shaft in described pit shaft.Figure 15 is prior art
In rod-pumped well hollow stem process structure schematic diagram.As shown in Figure 15,1 is stratum, and 2 is cement sheath, and 3 is sleeve pipe, and 4 is hydrodynamic
Face, 5 is heating cable, and 6 is oil pipe, and 7 is liquid in well, and 8 is oil reservoir.
S103：Described pit shaft and hollow stem are divided into by multiple wellbore sections, hollow bar segment according to described step-length.In tool
It is assumed that the total depth of pit shaft is 1000 meters in the embodiment of body, as shown in figure 15, the height in hydrodynamic face is 300 meters, setting
Step-length is 100 meters, then altogether pit shaft can be divided into 10 wellbore sections in this embodiment, be followed successively by 0- from shaft bottom to well head
100 meters, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice,
800-900 rice, 900-1000 rice.Altogether hollow stem can be divided into 10 hollow bar segment in this embodiment, from shaft bottom to well
Mouth be followed successively by 0-100 rice, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice,
700-800 rice, 800-900 rice, 900-1000 rice.
S104：The temperature of liquid in the plurality of wellbore section, multiple hollow stem are determined respectively according to described data information
The temperature of fluid in section.Fig. 2 is the particular flow sheet of step S104.
S105：The temperature composition of liquid in the temperature of fluid and described multiple wellbore sections in the plurality of hollow bar segment
Temperature in wellbore field distribution.
Fig. 2 is the particular flow sheet of step S104, and as shown in Figure 2, step S104 specifically includes：
S201：Determine each wellbore section, the thermal resistance of hollow bar segment successively.Fig. 3 is the particular flow sheet of step S201.
S202：Determine the thermal resistivity of each wellbore section, hollow bar segment according to described thermal resistance.
S203：Integral constant is determined according to wellhead temperature, bottom hole temperature (BHT).
S204：According to described integral constant and described thermal resistivity determine the temperature of liquid in each wellbore section,
The temperature of fluid in each hollow bar segment.
From the figure 3, it may be seen that step S201 specifically includes：
S301：Judge whether each described wellbore section, hollow bar segment are in described hydrodynamic face and described pit shaft pair successively
Between the well head answered；
S302：When being judged as YES, the multiple pit shafts between the well head corresponding with described pit shaft of described hydrodynamic face will be in
Section, hollow bar segment are set to first kind section；
S303：Otherwise, by the multiple wellbore sections being between the shaft bottom corresponding with described pit shaft of described hydrodynamic face, hollow stem
Section is set to Equations of The Second Kind section；
In a particular embodiment it is assumed that the total depth of pit shaft is 1000 meters, as shown in figure 15, the height in hydrodynamic face is
300 meters, pit shaft as 100 meters, then can be divided into 10 wellbore sections in this embodiment, from shaft bottom to well by step-length altogether that set
Mouth be followed successively by 0-100 rice, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice,
700-800 rice, 800-900 rice, 900-1000 rice.Altogether hollow stem can be divided into 10 hollow bar segment in this embodiment,
Be followed successively by from shaft bottom to well head 0-100 rice, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice,
600-700 rice, 700-800 rice, 800-900 rice, 900-1000 rice.
According to step S301 to step S303, hollow bar segment 0-100 rice, 100-200 rice, 200-300 rice are
Two class sections, hollow bar segment 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice, 800-900 rice,
900-1000 rice is first kind section.Wellbore section 0-100 rice, 100-200 rice, 200-300 rice are Equations of The Second Kind section, wellbore section
300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice, 800-900 rice, 900-1000 rice are
One class section.
S304：Determine the thermal resistance of described first kind section according to described data information；
S305：Determine the thermal resistance of described Equations of The Second Kind section according to described data information.
Fig. 4 is the particular flow sheet of step S304, and as shown in Figure 4, this step specifically includes：
S401：True according to described formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well production time, wellbore radius
Determine the thermal resistance on stratum, in a particular embodiment, the thermal resistance R on stratum1Represent, then：
Wherein, KeFor formation thermal conductivity, a is the average coefficient of heat transfer in stratum, and t is the oil well production time, rhFor pit shaft half
Footpath.
S402：Determine the thermal resistance of cement sheath according to described sleeve outer wall radius, cement sheath thermal conductivity factor, wellbore radius,
In a particular embodiment, the thermal resistance R of cement sheath2Represent, then：
Wherein, rcoFor sleeve outer wall radius, KcemFor cement sheath thermal conductivity factor, rhFor wellbore radius.
S403：Determine the heat of casing wall according to described sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall radius
Resistance, in a particular embodiment, stratum heat transfer thermal resistance R3Represent, then：
Wherein, KcasFor sleeve pipe thermal conductivity factor, rciFor internal surface of sleeve pipe radius, rcoFor sleeve outer wall radius.
S404：Determined according to described annular space radiation heat transfer coefficient, annular space free convection heat transfer coefficient, internal surface of sleeve pipe radius
Thermal convection current thermal resistance between air in oil jacket annular space and sleeve pipe, in a particular embodiment, between liquid and sleeve pipe, heat is right
Flowing liquid thermal resistance R4Represent, then：
Wherein, hrFor annular space radiation heat transfer coefficient, hcFor annular space free convection heat transfer coefficient, rciFor internal surface of sleeve pipe radius.
Annular space radiation heat transfer coefficient hrCalculated by equation below：
Wherein, δ is Stefan-Boltzmann constant, 2.189 × 10-8W/(m2·K)；FtciOutside for oil pipe or heat-insulated pipe
Wall surface is to internal surface of sleeve pipe surface emissivity coefficient of efficiency；εoFor adiabatic pipe outer wall blackness；εciFor internal surface of sleeve pipe blackness.For oil
The absolute temperature of pipe outer wall.TtoFor oil-pipe external wall temperature, TciFor internal surface of sleeve pipe temperature,Absolute temperature for internal surface of sleeve pipe.
Annular space free convection heat transfer coefficient hcCalculated by equation below：
Wherein, Gr is Grashof number, and Pr is Prandtl number, KhaFor the thermal conductivity factor of annular fluid, W/ (m K)；G is
Acceleration of gravity, m/s2；ρanFor annular fluid in mean temperature TanUnder density, kg/m3；UanFor annular fluid in mean temperature
TanUnder viscosity, mPa s；CanFor annular fluid in mean temperature TanUnder thermal capacitance, J (m3·K).TanFlat for oil jacket annular space
All temperature.
S405：According to described oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall radius determine oil pipe inside and outside wall it
Between thermal resistance, in a particular embodiment, heat transfer thermal resistance R between oil pipe inside and outside wall5Represent, then：
Wherein, KtubFor oil pipe thermal conductivity factor, rtiFor tube inner wall radius, rtoFor oil-pipe external wall radius.
S406：According to described crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem exterior radius, oil
Inside pipe wall radius determines the thermal convection current liquid thermal resistance between liquid and oil pipe, in a particular embodiment, in gas and oil pipe
Thermal convection current thermal resistance R between wall6Represent, then：
Wherein, λ0For crude oil thermal conductivity factor, fwFor liquid aqueous rate, λwFor the thermal conductivity factor of water, roFor hollow stem outer wall half
Footpath, rtiFor tube inner wall radius.
S407：According to described relative oil density, hollow stem inwall radius, hollow stem thermal conductivity factor KRodDetermine institute
State the thermal convection current liquid thermal resistance between the inside and outside wall of hollow stem.In a particular embodiment, between the inside and outside wall of hollow stem
Thermal convection current liquid thermal resistance R7Represent, then：
Wherein, riFor relative oil density, roFor hollow stem exterior radius, KRodThermal conductivity factor for hollow stem.
S408：Thermal convection current thermal resistance between being determined in hot fluid and hollow stem according to the thermal conductivity factor of described water.In tool
In the embodiment of body, in hot fluid and hollow stem between thermal convection current thermal resistance R8Represent, then：
Wherein, λwThermal conductivity factor for water.
Fig. 5 is the particular flow sheet of step S305, and as shown in Figure 5, this step specifically includes：
S501：True according to described formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well production time, wellbore radius
Determine the thermal resistance on stratum, in a particular embodiment, the thermal resistance R on stratum1Represent, then：
Wherein, KeFor formation thermal conductivity, a is the average coefficient of heat transfer in stratum, and t is the oil well production time, rhFor pit shaft half
Footpath.
S502：Determine the thermal resistance of cement sheath according to described sleeve outer wall radius, cement sheath thermal conductivity factor, wellbore radius,
In a particular embodiment, the thermal resistance R of cement sheath2Represent, then：
Wherein, rcoFor sleeve outer wall radius, KcemFor cement sheath thermal conductivity factor, rhFor wellbore radius.
S503：Determine the heat of casing wall according to described sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall radius
Resistance, in a particular embodiment, stratum heat transfer thermal resistance R3Represent, then：
Wherein, KcasFor sleeve pipe thermal conductivity factor, rciFor internal surface of sleeve pipe radius, rcoFor sleeve outer wall radius.
S504：According to described crude oil thermal conductivity factor, the thermal conductivity factor of water, liquid aqueous rate, oil-pipe external wall radius, sleeve pipe
Inwall radius determines the thermal convection current liquid thermal resistance between liquid and sleeve pipe, in a particular embodiment, between liquid and sleeve pipe
Thermal convection current liquid thermal resistance R '4Represent, determined by equation below：
Wherein, λo=0.01172 (1-0.00054T)/γo
λw=3.51153-0.04436 (T+273.15)+2.41233 × 10-4×(T+273.15)2-6.051×10-7×
(T+273.15)3+7.22766×10-10(T+273.15)4-3.3716×10-13(T+273.15)5
R′4For thermal convection current liquid thermal resistance, λ between liquid and sleeve pipeoFor crude oil thermal conductivity factor, γoFor relative oil density,
λwFor water thermal conductivity factor, fwFor liquid aqueous rate, rtoFor oil-pipe external wall radius, T is annular fluid temperature between liquid and sleeve pipe
Value.
S505：According to described oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall radius determine oil pipe inside and outside wall it
Between thermal resistance, in a particular embodiment, heat transfer thermal resistance R between oil pipe inside and outside wall5Represent, then：
Wherein, KtubFor oil pipe thermal conductivity factor, rtiFor tube inner wall radius, rtoFor oil-pipe external wall radius.
S506：According to described crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem exterior radius, oil
Inside pipe wall radius determines the thermal convection current liquid thermal resistance between liquid and oil pipe, in a particular embodiment, liquid and oil pipe it
Between thermal convection current liquid thermal resistance R6Represent, then：
Wherein, λ0For crude oil thermal conductivity factor, fwFor liquid aqueous rate, λwFor the thermal conductivity factor of water, roFor hollow stem outer wall half
Footpath, rtiFor tube inner wall radius.
S507：According to described relative oil density, hollow stem inwall radius, KRodDetermine the inside and outside wall of described hollow stem
Between thermal convection current liquid thermal resistance.In a particular embodiment, the thermal convection current liquid thermal resistance between the inside and outside wall of hollow stem is used
R7Represent, then：
Wherein, riFor relative oil density, roFor hollow stem exterior radius, KRodFor.
S508：Thermal convection current thermal resistance between being determined in hot fluid and hollow stem according to the thermal conductivity factor of described water.In tool
In the embodiment of body, in hot fluid and hollow stem between thermal convection current thermal resistance R8Represent, then：
Wherein, λwThermal conductivity factor for water.
Fig. 6 is the particular flow sheet of the embodiment one of step S202 in Fig. 2, in this embodiment, when described
When wellbore section, hollow bar segment are first kind section, step S202 includes：
S601：According to the heat between the thermal convection current liquid thermal resistance between described liquid and oil pipe, the inside and outside wall of hollow stem
Thermal convection current thermal resistance between in convective liquid thermal resistance, hot fluid and hollow stem determines the first thermal resistivity, in specific embodiment party
In formula, the first thermal resistivity K11Represent, then：
Wherein, R6Represent the thermal convection current liquid thermal resistance between liquid and oil pipe, R7Represent the heat between the inside and outside wall of hollow stem
Convective liquid thermal resistance, R8Thermal convection current thermal resistance between representing in hot fluid and hollow stem is used.
S602：According to the air in the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, oil jacket annular space
Thermal convection current liquid thermal resistance between thermal resistance between thermal convection current thermal resistance and sleeve pipe between, oil pipe inside and outside wall, liquid and oil pipe is true
Fixed second thermal resistivity, in a particular embodiment, the second thermal resistivity K12Represent, then：
Wherein, R1For the thermal resistance on stratum, R2For the thermal resistance of cement sheath, R3For the thermal resistance of casing wall, R4For in oil jacket annular space
Thermal convection current thermal resistance between air and sleeve pipe, R5For the thermal resistance between oil pipe inside and outside wall, R6For the thermal convection current between liquid and oil pipe
Liquid thermal resistance.
S603：True according to the thermal convection current liquid thermal resistance between the thermal resistance between described oil pipe inside and outside wall, liquid and oil pipe
Fixed 3rd thermal resistivity, in a particular embodiment, the 3rd thermal resistivity K13Represent, then：
Wherein, R5For the thermal resistance between oil pipe inside and outside wall, R6For the thermal convection current liquid thermal resistance between liquid and oil pipe.
S604：According to the air in the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, oil jacket annular space
Thermal convection current thermal resistance and sleeve pipe between determines the 4th thermal resistivity, in a particular embodiment, the 4th thermal resistivity K14Table
Show, then：
Wherein, R1For the thermal resistance on stratum, R2For the thermal resistance of cement sheath, R3For the thermal resistance of casing wall, R4For in oil jacket annular space
Thermal convection current thermal resistance between air and sleeve pipe.
Fig. 7 is the particular flow sheet of the embodiment two of step S202 in Fig. 2, in this embodiment, when described
When wellbore section, hollow bar segment are Equations of The Second Kind section, step S202 includes：
S701：According to the heat between the thermal convection current liquid thermal resistance between described liquid and oil pipe, the inside and outside wall of hollow stem
Thermal convection current thermal resistance between in convective liquid thermal resistance, hot fluid and hollow stem determines the first thermal resistivity, in specific embodiment party
In formula, the first thermal resistivity K11Represent, then：
Wherein, R6Represent the thermal convection current liquid thermal resistance between liquid and oil pipe, R7Represent the heat between the inside and outside wall of hollow stem
Convective liquid thermal resistance, R8Thermal convection current thermal resistance between representing in hot fluid and hollow stem is used.
S702：According between the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, liquid and sleeve pipe
Thermal convection current liquid thermal resistance between thermal resistance between thermal convection current liquid thermal resistance, oil pipe inside and outside wall, liquid and oil pipe determines the second heat
Resistance coefficient, in a particular embodiment, the second thermal resistivity K12Represent, then：
Wherein, R1For the thermal resistance on stratum, R2For the thermal resistance of cement sheath, R3For the thermal resistance of casing wall, R '4For liquid and sleeve pipe
Between thermal convection current liquid thermal resistance, R5For the thermal resistance between oil pipe inside and outside wall, R6For the thermal convection current liquid heat between liquid and oil pipe
Resistance.
S703：True according to the thermal convection current liquid thermal resistance between the thermal resistance between described oil pipe inside and outside wall, liquid and oil pipe
Fixed 3rd thermal resistivity, in a particular embodiment, the 3rd thermal resistivity K13Represent, then：
Wherein, R5For the thermal resistance between oil pipe inside and outside wall, R6For the thermal convection current liquid thermal resistance between liquid and oil pipe.
S704：According between the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, liquid and sleeve pipe
Thermal convection current liquid thermal resistance determines the 4th thermal resistivity, in a particular embodiment, the 4th thermal resistivity K14Represent, then：
Wherein, R1For the thermal resistance on stratum, R2For the thermal resistance of cement sheath, R3For the thermal resistance of casing wall, R '4For liquid and sleeve pipe
Between thermal convection current liquid thermal resistance.
As shown in Figure 2, step S104 also includes：
S203：Integral constant is determined according to wellhead temperature, bottom hole temperature (BHT).In a particular embodiment, integral constant is
C1、C2、C3、C4, can be calculated by following formula binding characteristic values (i.e. known wellhead temperature and bottom hole temperature (BHT)).
S204：According to described integral constant and described thermal resistivity determine the temperature of liquid in each wellbore section,
The temperature of fluid in each hollow bar segment.The present invention employs following energy-balance equation in calculating temperature field：Unit interval
Flow into the change of cell cube interior energy in the enthalpy+potential variation=unit interval of enthalpy-unit interval outlet unit body in unit body.
Thus derive that the temperature of liquid in wellbore section is calculated by equation below：
In hollow bar segment, the temperature of fluid is calculated by equation below：
Wherein, θ1For the temperature of liquid in wellbore section, θ2For the temperature of fluid in hollow bar segment, t is to mix hot water or thin oil
Or the temperature of chemical agent (temperature mixing the temperature of hot water or thin oil or chemical agent in well head, close to wellhead temperature, mixes hot water
Or the temperature in shaft bottom for the temperature of thin oil or chemical agent is close to bottom hole temperature (BHT)), K11It is outflow in hollow stem for the first thermal resistivity
Overall heat-transfer coefficient between body, W/m DEG C；K12It is between hollow stem and oil pipe annular fluid and surrounding formation for the second thermal resistivity
Overall heat-transfer coefficient, W/m DEG C；K13It is the overall heat-transfer coefficient between outer fluid in oil pipe for the 3rd thermal resistivity, W/m DEG C；K14For
4th thermal resistivity is the fluid in oil jacket annular space and the overall heat-transfer coefficient between bottom, and m is geothermal gradient, and l is along well depth direction
Length, g be acceleration of gravity, q be heat source strength, tsFor earth's surface thermostat layer temperature.W is the water equivalent of gas mixture, We
For mixing the water equivalent of hot water or thin oil or chemical agent, W/ DEG C；MfFor crude quality flow, CfFor crude oil specific heat, MgFor water quality
Flow, CgFor the specific heat of water, W can be calculated as below：W=MfCf+MgCg.
Fig. 8 is a kind of structural frames of system determining temperature in wellbore field distribution based on hollow stem provided in an embodiment of the present invention
Figure, as shown in Figure 8, described system includes：
Data information acquisition device 101, for obtaining the data information related to hollow stem and pit shaft.
In the particular embodiment, the data information related to hollow stem and pit shaft includes formation thermal conductivity, stratum
Average coefficient of heat transfer, the oil well production time, wellbore radius, sleeve outer wall radius, cement sheath thermal conductivity factor, sleeve pipe thermal conductivity factor,
Internal surface of sleeve pipe radius, sleeve outer wall radius, annular space radiation heat transfer coefficient, annular space free convection heat transfer coefficient, oil-pipe external wall radius,
Oil pipe thermal conductivity factor, tube inner wall radius, crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem outer wall half
Footpath, relative oil density, hollow stem inwall radius, the thermal conductivity factor of hollow stem.
Step size settings device 102, for the depth-set step-length according to the height in hydrodynamic face and pit shaft in described pit shaft.
Figure 15 is rod-pumped well hollow stem process structure schematic diagram of the prior art.As shown in Figure 15,1 is stratum, and 2 is cement sheath, 3
For sleeve pipe, 4 is hydrodynamic face, and 5 is heating cable, and 6 is oil pipe, and 7 is liquid in well, and 8 is oil reservoir.
Sectioning 103, for being divided into multiple wellbore sections, sky according to described step-length by described pit shaft and hollow stem
Core bar section.In a particular embodiment it is assumed that the total depth of pit shaft is 1000 meters, as shown in figure 15, the height in hydrodynamic face is
300 meters, pit shaft as 100 meters, then can be divided into 10 wellbore sections in this embodiment, from shaft bottom to well by step-length altogether that set
Mouth be followed successively by 0-100 rice, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice,
700-800 rice, 800-900 rice, 900-1000 rice.Altogether hollow stem can be divided into 10 hollow bar segment in this embodiment,
Be followed successively by from shaft bottom to well head 0-100 rice, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice,
600-700 rice, 700-800 rice, 800-900 rice, 900-1000 rice.
Temperature determines device 104, for determining liquid in the plurality of wellbore section respectively according to described data information
The temperature of fluid in temperature, multiple hollow bar segment.Fig. 2 is the particular flow sheet of step S104.
Thermo parameters method determines device 105, the temperature for fluid in the plurality of hollow bar segment and described multiple
The temperature composition temperature in wellbore field distribution of liquid in wellbore section.
Fig. 9 is a kind of temperature being determined based on hollow stem in the system of temperature in wellbore field distribution provided in an embodiment of the present invention
Determine the concrete structure block diagram of device 104, as shown in Figure 9, temperature determines that device 104 specifically includes：
Thermal resistance determining module 201, for determining the thermal resistance of each wellbore section, hollow bar segment successively.Figure 10 determines for thermal resistance
Module, for 201 concrete figure.
Thermal resistivity determining module 202, for determining the thermal resistance system of each wellbore section, hollow bar segment according to described thermal resistance
Number.
Integral constant determining module 203, for determining integral constant according to wellhead temperature, bottom hole temperature (BHT).
Temperature determination module 204, for determining each pit shaft according to described integral constant and described thermal resistivity
The temperature of fluid in the temperature of liquid in section, each hollow bar segment.
As shown in Figure 10, thermal resistance determining module 201 specifically includes：
Judging unit 301, for judge successively each described wellbore section, hollow bar segment whether be in described hydrodynamic face with
Between the corresponding well head of described pit shaft；
First kind section arranging unit 302, for when described judge module is judged as YES, will be in described hydrodynamic face with
Multiple wellbore sections between the corresponding well head of described pit shaft, hollow bar segment are set to first kind section；
Equations of The Second Kind section arranging unit 303, for when judge module is judged as NO, being in described hydrodynamic face and described well
Multiple wellbore sections between the corresponding shaft bottom of cylinder, hollow bar segment are set to Equations of The Second Kind section；
In a particular embodiment it is assumed that the total depth of pit shaft is 1000 meters, as shown in figure 15, the height in hydrodynamic face is
300 meters, pit shaft as 100 meters, then can be divided into 10 wellbore sections in this embodiment, from shaft bottom to well by step-length altogether that set
Mouth be followed successively by 0-100 rice, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice,
700-800 rice, 800-900 rice, 900-1000 rice.Altogether hollow stem can be divided into 10 hollow bar segment in this embodiment,
Be followed successively by from shaft bottom to well head 0-100 rice, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice,
600-700 rice, 700-800 rice, 800-900 rice, 900-1000 rice.
According to judging unit 301 to Equations of The Second Kind section arranging unit 03, hollow bar segment 0-100 rice, 100-200 rice,
200-300 rice is Equations of The Second Kind section, hollow bar segment 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800
Rice, 800-900 rice, 900-1000 rice are first kind section.Wellbore section 0-100 rice, 100-200 rice, 200-300 rice are second
Class section, wellbore section 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice, 800-900 rice, 900-
1000 meters are first kind section.
First thermal resistance determining unit 304, for determining the thermal resistance of described first kind section according to described data information；
Second thermal resistance determining unit 305, for determining the thermal resistance of described Equations of The Second Kind section according to described data information.
Figure 11 be provided in an embodiment of the present invention a kind of determined based on hollow stem in the system of temperature in wellbore field distribution
The concrete structure block diagram of one thermal resistance determining unit 304, as shown in Figure 11, the first thermal resistance determining unit specifically includes：
First thermal resistance determining unit 401, for according to described formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well
Production time, wellbore radius determine the thermal resistance on stratum, in a particular embodiment, the thermal resistance R on stratum1Represent.
Second thermal resistance determining unit 402, for according to described sleeve outer wall radius, cement sheath thermal conductivity factor, pit shaft half
Footpath determines the thermal resistance of cement sheath, in a particular embodiment, the thermal resistance R of cement sheath2Represent.
3rd thermal resistance determining unit 403, for according to described sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall
Radius determines the thermal resistance of casing wall, in a particular embodiment, stratum heat transfer thermal resistance R3Represent.
4th thermal resistance determining unit 404, for according to described annular space radiation heat transfer coefficient, annular space free convection heat transfer system
Number, internal surface of sleeve pipe radius determine the thermal convection current thermal resistance between air and sleeve pipe in oil jacket annular space, in a particular embodiment,
Thermal convection current liquid thermal resistance R between liquid and sleeve pipe4Represent.
5th thermal resistance determining unit 405, for according to described oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall
Radius determines the thermal resistance between oil pipe inside and outside wall, in a particular embodiment, heat transfer thermal resistance R between oil pipe inside and outside wall5
Represent.
6th thermal resistance determining unit 406, for according to described crude oil thermal conductivity factor, liquid aqueous rate, water heat conduction system
Number, hollow stem exterior radius, tube inner wall radius determine the thermal convection current liquid thermal resistance between liquid and oil pipe, implement specific
In mode, thermal convection current thermal resistance R between gas and tube inner wall6Represent.
7th thermal resistance determining unit 407, for according to described relative oil density, hollow stem inwall radius, hollow stem
Thermal conductivity factor determine the thermal convection current liquid thermal resistance between the inside and outside wall of described hollow stem.In a particular embodiment, hollow
Thermal convection current liquid thermal resistance R between the inside and outside wall of bar7Represent.
8th thermal resistance determining unit 408, for determining in hot fluid and hollow stem it according to the thermal conductivity factor of described water
Between thermal convection current thermal resistance.In a particular embodiment, the thermal convection current thermal resistance R between in hot fluid and hollow stem8Represent.
Figure 12 be provided in an embodiment of the present invention a kind of determined based on hollow stem in the system of temperature in wellbore field distribution
The concrete structure block diagram of two thermal resistance determining units 305, as shown in Figure 12, the second thermal resistance determining unit specifically includes：
First thermal resistance determining unit 501, for according to described formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well
Production time, wellbore radius determine the thermal resistance on stratum, in a particular embodiment, the thermal resistance R on stratum1Represent.
Second thermal resistance determining unit 502, for according to described sleeve outer wall radius, cement sheath thermal conductivity factor, pit shaft half
Footpath determines the thermal resistance of cement sheath, in a particular embodiment, the thermal resistance R of cement sheath2Represent.
Second thermal resistance determining unit 503, for according to described sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall
Radius determines the thermal resistance of casing wall, in a particular embodiment, stratum heat transfer thermal resistance R3Represent.
4th thermal resistance determining unit 504, for according to the described crude oil thermal conductivity factor, thermal conductivity factor of water, liquid aqueous
Rate, oil-pipe external wall radius, internal surface of sleeve pipe radius determine the thermal convection current liquid thermal resistance between liquid and sleeve pipe, in specific embodiment party
Thermal convection current liquid thermal resistance R ' in formula, between liquid and sleeve pipe4Represent.
5th thermal resistance determining unit 505, for according to described oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall
Radius determines the thermal resistance between oil pipe inside and outside wall, in a particular embodiment, heat transfer thermal resistance R between oil pipe inside and outside wall5
Represent.
6th thermal resistance determining unit 506, for according to described crude oil thermal conductivity factor, liquid aqueous rate, water heat conduction system
Number, hollow stem exterior radius, tube inner wall radius determine the thermal convection current liquid thermal resistance between liquid and oil pipe, implement specific
Thermal convection current liquid thermal resistance R in mode, between liquid and oil pipe6Represent.
7th thermal resistance determining unit 507, for according to described relative oil density, hollow stem inwall radius, hollow stem
Thermal conductivity factor determine the thermal convection current liquid thermal resistance between the inside and outside wall of described hollow stem.In a particular embodiment, hollow
Thermal convection current liquid thermal resistance R between the inside and outside wall of bar7Represent.
8th thermal resistance determining unit 508, for determining in hot fluid and hollow stem it according to the thermal conductivity factor of described water
Between thermal convection current thermal resistance.In a particular embodiment, the thermal convection current thermal resistance R between in hot fluid and hollow stem8Represent.
Figure 13 is a kind of heat being determined based on hollow stem in the system of temperature in wellbore field distribution provided in an embodiment of the present invention
The concrete structure block diagram of the embodiment one of resistance coefficient determination module 202, in this embodiment, when described wellbore section, sky
When core bar section is first kind section, thermal resistivity determining module 202 includes：
First thermal resistivity determining unit 601, for according to the thermal convection current liquid thermal resistance between described liquid and oil pipe,
Thermal convection current thermal resistance between in thermal convection current liquid thermal resistance between the inside and outside wall of hollow stem, hot fluid and hollow stem determines the first heat
Resistance coefficient, in a particular embodiment, the first thermal resistivity K11Represent.
Second thermal resistivity determining unit 602, for according to the thermal resistance on described stratum, the thermal resistance of cement sheath, casing wall
Thermal resistance, the air in oil jacket annular space and the thermal resistance between the thermal convection current thermal resistance between sleeve pipe, oil pipe inside and outside wall, liquid and oil pipe
Between thermal convection current liquid thermal resistance determine the second thermal resistivity, in a particular embodiment, the second thermal resistivity K12Represent.
3rd thermal resistivity determining unit 603, for according to the thermal resistance between described oil pipe inside and outside wall, liquid and oil pipe
Between thermal convection current liquid thermal resistance determine the 3rd thermal resistivity, in a particular embodiment, the 3rd thermal resistivity K13Represent.
4th thermal resistivity determining unit 604, for according to the thermal resistance on described stratum, the thermal resistance of cement sheath, casing wall
Thermal resistance, the air in oil jacket annular space and the thermal convection current thermal resistance between sleeve pipe determine the 4th thermal resistivity, in specific embodiment party
In formula, the 4th thermal resistivity K14Represent.
Figure 14 is a kind of heat being determined based on hollow stem in the system of temperature in wellbore field distribution provided in an embodiment of the present invention
The concrete structure block diagram of the embodiment two of resistance coefficient determination module 202, in this embodiment, when described wellbore section, sky
When core bar section is Equations of The Second Kind section, thermal resistivity determining module 202 includes：
Thermal resistivity the first determining unit 701, for according to the thermal convection current liquid thermal resistance between described liquid and oil pipe,
Thermal convection current thermal resistance between in thermal convection current liquid thermal resistance between the inside and outside wall of hollow stem, hot fluid and hollow stem determines the first heat
Resistance coefficient, in a particular embodiment, the first thermal resistivity K11Represent.
Thermal resistivity the second determining unit 702, for according to the thermal resistance on described stratum, the thermal resistance of cement sheath, casing wall
Thermal resistance, the thermal convection current liquid thermal resistance between liquid and sleeve pipe, the thermal resistance between oil pipe inside and outside wall, the heat between liquid and oil pipe
Convective liquid thermal resistance determines the second thermal resistivity, in a particular embodiment, the second thermal resistivity K12Represent.
Thermal resistivity the 3rd determining unit 703, for according to the thermal resistance between described oil pipe inside and outside wall, liquid and oil pipe
Between thermal convection current liquid thermal resistance determine the 3rd thermal resistivity, in a particular embodiment, the 3rd thermal resistivity K13Represent.
Thermal resistivity the 4th determining unit 704, for according to the thermal resistance on described stratum, the thermal resistance of cement sheath, casing wall
Thermal resistance, the thermal convection current liquid thermal resistance between liquid and sleeve pipe determine the 4th thermal resistivity, in a particular embodiment, the 4th
Thermal resistivity K14Represent.
As shown in Figure 9, temperature determines that device 104 also includes：
Integral constant determining module 203, for determining integral constant according to wellhead temperature, bottom hole temperature (BHT).Specifically real
Apply in mode, integral constant is C1、C2、C3、C4, formula binding characteristic value (i.e. known wellhead temperature and shaft bottom temperature can be passed through
Degree) calculated.
Temperature determination module 204, for determining each pit shaft according to described integral constant and described thermal resistivity
The temperature of fluid in the temperature of liquid in section, each hollow bar segment.The present invention employs following energy in calculating temperature field and puts down
Weighing apparatus equation：Single in the enthalpy+potential variation=unit interval of enthalpy-unit interval outlet unit body in unit interval inflow unit body
The change of the internal energy of unit.Thus derive that the temperature of liquid in wellbore section is calculated by equation below：
In hollow bar segment, the temperature of fluid is calculated by equation below：
Wherein, θ1For the temperature of liquid in wellbore section, θ2For the temperature of fluid in hollow bar segment, t is to mix hot water or thin oil
Or the temperature of chemical agent (temperature mixing the temperature of hot water or thin oil or chemical agent in well head, close to wellhead temperature, mixes hot water
Or the temperature in shaft bottom for the temperature of thin oil or chemical agent is close to bottom hole temperature (BHT)), K11It is outflow in hollow stem for the first thermal resistivity
Overall heat-transfer coefficient between body, W/m DEG C；K12It is between hollow stem and oil pipe annular fluid and surrounding formation for the second thermal resistivity
Overall heat-transfer coefficient, W/m DEG C；K13It is the overall heat-transfer coefficient between outer fluid in oil pipe for the 3rd thermal resistivity, W/m DEG C；K14For
4th thermal resistivity is the fluid in oil jacket annular space and the overall heat-transfer coefficient between bottom, and m is geothermal gradient, and l is along well depth direction
Length, g be acceleration of gravity, q be heat source strength, tsFor earth's surface thermostat layer temperature.W is the water equivalent of gas mixture, We
For mixing the water equivalent of hot water or thin oil or chemical agent, W/ DEG C；MfFor crude quality flow, CfFor crude oil specific heat, MgFor water quality
Flow, CgFor the specific heat of water, W can be calculated as below：W=MfCf+MgCg.
With reference to specific embodiment, technical scheme is discussed in detail.In a particular embodiment, false
If the length of the total depth of pit shaft, hollow stem is 1000 meters, as shown in figure 15, the height in hydrodynamic face is 300 meters, the step of setting
A length of 100 meters, then altogether pit shaft can be divided into 10 wellbore sections in this embodiment, be followed successively by 0-100 from shaft bottom to well head
Rice, 100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice, 800-
900 meters, 900-1000 rice.Altogether hollow stem can be divided into 10 hollow bar segment, be followed successively by from shaft bottom to well head 0-100 rice,
100-200 rice, 200-300 rice, 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice, 800-900
Rice, 900-1000 rice.According to step S301 to step S303, wellbore section 0-100 rice, 100-200 rice, 200-300 rice are equal
For Equations of The Second Kind section, wellbore section 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice, 800-900 rice,
900-1000 rice is first kind section.Hollow bar segment 0-100 rice, 100-200 rice, 200-300 rice are Equations of The Second Kind section, hollow stem
Section 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice, 800-900 rice, 900-1000 rice are
First kind section.Introduce separately below and how fluid in the fluid temperature of each wellbore section, hollow bar segment is determined according to data information
Temperature.
1st, Equations of The Second Kind section
Wellbore section 0-100 rice, 100-200 rice, 200-300 rice, hollow bar segment 0-100 rice, 100-200 rice, 200-300 rice
It is Equations of The Second Kind section, illustrate taking wellbore section 0-100 rice as a example.
(1) thermal resistance on stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, liquid and the set of wellbore section 0-100 rice, are determined
Thermal convection current liquid thermal resistance between thermal resistance between thermal convection current liquid thermal resistance between pipe, oil pipe inside and outside wall, liquid and oil pipe, sky
Thermal convection current thermal resistance between in thermal convection current liquid thermal resistance between the inside and outside wall of core bar, hot fluid and hollow stem.
(2), according to the thermal convection current liquid heat between thermal resistance, the thermal resistance of cement sheath, the thermal resistance of casing wall, liquid and sleeve pipe
Heat between thermal convection current liquid thermal resistance between thermal resistance between resistance, oil pipe inside and outside wall, liquid and oil pipe, the inside and outside wall of hollow stem
Thermal convection current thermal resistance between in convective liquid thermal resistance, hot fluid and hollow stem determine the first thermal resistivity, the second thermal resistivity, the
Three thermal resistivities, the 4th thermal resistivity.
(3), integral constant C is determined according to wellhead temperature, bottom hole temperature (BHT)1、C2、C3、C4.
(4), basisCalculate liquid in wellbore section
Temperature.
According toCalculate the temperature of fluid in hollow bar segment.
2nd, first kind wellbore section
Wellbore section 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-800 rice, 800-900 rice,
900-1000 rice is first kind section, hollow bar segment 300-400 rice, 400-500 rice, 500-600 rice, 600-700 rice, 700-
800 meters, 800-900 rice, 900-1000 rice be first kind section, illustrate taking lower wellbore section 300-400 rice as a example.
(1) thermal resistance on stratum of wellbore section 300-400 rice, the thermal resistance of cement sheath, the thermal resistance of casing wall, oil jacket ring, are determined
Thermal convection current between thermal resistance between thermal convection current thermal resistance between aerial air and sleeve pipe, oil pipe inside and outside wall, liquid and oil pipe
Thermal convection current thermal resistance between in thermal convection current liquid thermal resistance between liquid thermal resistance, the inside and outside wall of hollow stem, hot fluid and hollow stem.
(2), according to the heat between the air in thermal resistance, the thermal resistance of cement sheath, the thermal resistance of casing wall, oil jacket annular space and sleeve pipe
Thermal convection current liquid thermal resistance between thermal resistance between thermal-convection resistance, oil pipe inside and outside wall, liquid and oil pipe, hollow stem inside and outside wall it
Between thermal convection current liquid thermal resistance, in hot fluid and hollow stem between thermal convection current thermal resistance determine the first thermal resistivity, the second thermal resistance
Coefficient, the 3rd thermal resistivity, the 4th thermal resistivity.
(3), integral constant C is determined according to wellhead temperature, bottom hole temperature (BHT)1、C2、C3、C4.
(4), basisCalculate liquid in wellbore section
Temperature.
According toCalculate the temperature of fluid in hollow bar segment.
Without the iterative calculation carrying out total thermal conductivity factor, calculating speed faster, arbitrarily can arrange material calculation to this programme,
When step-length is less, computational accuracy is higher.The method has extraordinary stability and convergence, is more suitable for computer programming.
As follows：
1st, data prepares, and data mainly includes：Formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well production time, well
Cylinder radius, sleeve outer wall radius, cement sheath thermal conductivity factor, sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall radius, ring
Empty radiation heat transfer coefficient, annular space free convection heat transfer coefficient, oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall radius, former
Oily thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem exterior radius, relative oil density, hollow stem inwall half
Footpath, the thermal conductivity factor of hollow stem.
2. calculate thermal resistance.
3. calculate thermal resistivity K11, K12, K13, K14.
4. start to calculate from well head, l=0, k=1.
5. C1, C2, C3, C4 are calculated by known conditions.
6. calculate temperature θ of liquid in pit shaft1With fluid temperature (F.T.) θ in hollow stem2.
7.k=k+1, l=l+dl, return the 3rd step and continue iterative calculation.If l>Well depth, then iteration terminate.
In sum, the invention provides a kind of method and system determining temperature in wellbore field distribution based on hollow stem, it is
A kind of accurate scheme determining temperature in wellbore field distribution based on hollow stem, by obtaining the number related to hollow stem and pit shaft
According to data, according to setting step-length, pit shaft, hollow stem are divided into multiple wellbore sections, hollow bar segment, determine successively in each wellbore section
In the temperature of liquid, hollow bar segment, the temperature of fluid, has so then obtained temperature in wellbore field distribution, is subsequently selected rational well
Mouth mixes discharge capacity and temperature provides data foundation with the exploitation meeting existing viscous crude, special thick oil and super-viscous oil, and then improves
The efficiency of the exploitation of viscous crude, special thick oil and super-viscous oil.
The present invention be directed to the exploitation of viscous crude, special thick oil and super-viscous oil, the big feature of the main temperature influence of viscosity, carry
Go out a kind of method calculating hollow stem temperature in wellbore mixed with hot water field distribution, and define a set of interpretative tool.By to mixing heat
After water, the decoupled method analysis of Wellbore Temperature Field, can calculate every section of viscosity of crude, further can calculate rod liquid
Frictional force and pipe liquid frrction load, may finally calculate oil pumping machine polished rod load by roofbolt up stroke.The well of hot water can be adjusted
Mouth mixes discharge capacity and temperature, has reached the purpose optimizing hoisting system efficiency, cannot be only used for hollow stem temperature field mixed with hot water and divide
The calculating of cloth, it may also be used for mixing light oil, mix chemical agent, without the iterative calculation carrying out total thermal conductivity factor, calculating speed faster,
Material calculation can be arbitrarily set, and when step-length is less, computational accuracy is higher.There is extraordinary stability and convergence, more
Suitable computer programming.
The present invention is directed to the oil extraction that freezing point is high, paraffin content is high, viscosity is high, sets up description with thermodynamics general principle
The Mathematical Modeling of hollow stem temperature in wellbore mixed with hot water field distribution, and solved with numerical method, mix heat to understand and grasp
Bar, pit shaft and formation rock Temperature Distribution and variation tendency during water, are mixed with the rational well head of Instructing manufacture practical choice
Discharge capacity and temperature.
Computational methods involved in the present invention have good computational stability and higher computational accuracy, by this algorithm
And interpretative tool, can be very good hollow stem electric cable heating power and lifting technology parameter to be predicted and adjusts.
One of ordinary skill in the art will appreciate that realizing all or part of flow process in above-described embodiment method, Ke Yitong
Cross computer program to complete come the hardware to instruct correlation, described program can be stored in general computer read/write memory medium
In, this program is upon execution, it may include as the flow process of the embodiment of above-mentioned each method.Wherein, described storage medium can be magnetic
Dish, CD, read-only memory (Read-Only Memory, ROM) or random access memory (Random Access
Memory, RAM) etc..
Those skilled in the art are it will also be appreciated that various functions that the embodiment of the present invention is listed are by hardware or soft
Part is realizing the design requirement depending on specific application and whole system.Those skilled in the art can be for every kind of specific
Application, it is possible to use various methods realize described function, but this realization is understood not to protect beyond the embodiment of the present invention
The scope of shield.
Apply specific embodiment in the present invention principle of the present invention and embodiment are set forth, above example
Explanation be only intended to help and understand the method for the present invention and its core concept；Simultaneously for one of ordinary skill in the art,
According to the thought of the present invention, all will change in specific embodiments and applications, in sum, in this specification
Hold and should not be construed as limitation of the present invention.
Claims (10)
1. a kind of method determining temperature in wellbore field distribution based on hollow stem, is characterized in that, described method includes：
Obtain the data information related to hollow stem and pit shaft, described data information includes formation thermal conductivity, stratum is put down
All coefficient of heat transfer, oil well production time, wellbore radius, sleeve outer wall radius, cement sheath thermal conductivity factor, sleeve pipe thermal conductivity factor, sets
Inside pipe wall radius, annular space radiation heat transfer coefficient, annular space free convection heat transfer coefficient, oil-pipe external wall radius, oil pipe thermal conductivity factor, oil
Inside pipe wall radius, crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem exterior radius, relative oil density,
Hollow stem inwall radius, the thermal conductivity factor of hollow stem；
Depth-set step-length according to the height in hydrodynamic face and pit shaft in described pit shaft；
Described pit shaft and hollow stem are divided into by multiple wellbore sections, hollow bar segment according to described step-length；
Fluid in the temperature of liquid in the plurality of wellbore section, multiple hollow bar segment is determined respectively according to described data information
Temperature, this step includes determining each wellbore section, the thermal resistance of hollow bar segment successively；Each pit shaft is determined according to described thermal resistance
Section, the thermal resistivity of hollow bar segment；Integral constant is determined according to wellhead temperature, bottom hole temperature (BHT)；According to described integral constant with
And described thermal resistivity determines the temperature of fluid in the temperature of liquid in each wellbore section, each hollow bar segment；
The temperature composition temperature in wellbore of liquid in the temperature of fluid and described multiple wellbore sections in the plurality of hollow bar segment
Field distribution；
Wherein, determine that each wellbore section, the thermal resistance of hollow bar segment include successively：Judge each described wellbore section, hollow successively
Whether bar segment is between the well head corresponding with described pit shaft of described hydrodynamic face；When being judged as YES, described hydrodynamic face will be in
Multiple wellbore sections between well head corresponding with described pit shaft, hollow bar segment are set to first kind section；Otherwise, described hydrodynamic will be in
Multiple wellbore sections between the shaft bottom corresponding with described pit shaft of face, hollow bar segment are set to Equations of The Second Kind section；According to described data money
Material determines the thermal resistance of described first kind section；Determine the thermal resistance of described Equations of The Second Kind section according to described data information.
2. method according to claim 1, is characterized in that, determines the heat of described first kind section according to described data information
Resistance includes：
Determine the heat on stratum according to described formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well production time, wellbore radius
Resistance；
Determine the thermal resistance of cement sheath according to described sleeve outer wall radius, cement sheath thermal conductivity factor, wellbore radius；
Determine the thermal resistance of casing wall according to described sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall radius；
Oil jacket annular space is determined according to described annular space radiation heat transfer coefficient, annular space free convection heat transfer coefficient, internal surface of sleeve pipe radius
In air and sleeve pipe between thermal convection current thermal resistance；
Thermal resistance between oil pipe inside and outside wall is determined according to described oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall radius；
According to described crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem exterior radius, tube inner wall half
Footpath determines the thermal convection current liquid thermal resistance between liquid and oil pipe；
Determined in described hollow stem according to the thermal conductivity factor of described relative oil density, hollow stem inwall radius, hollow stem
Thermal convection current liquid thermal resistance between outer wall；
Thermal convection current thermal resistance between being determined in hot fluid and hollow stem according to the thermal conductivity factor of described water.
3. method according to claim 2, is characterized in that, when described wellbore section, hollow bar segment are first kind section, root
Determine each wellbore section according to described thermal resistance, the thermal resistivity of hollow bar segment includes：
According to the thermal convection current liquid heat between the thermal convection current liquid thermal resistance between described liquid and oil pipe, the inside and outside wall of hollow stem
Thermal convection current thermal resistance between in resistance, hot fluid and hollow stem determines the first thermal resistivity；
According between the air in the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, oil jacket annular space and sleeve pipe
Thermal convection current thermal resistance, the thermal resistance between oil pipe inside and outside wall, the thermal convection current liquid thermal resistance between liquid and oil pipe determine the second thermal resistance
Coefficient；
3rd thermal resistance is determined according to the thermal convection current liquid thermal resistance between the thermal resistance between described oil pipe inside and outside wall, liquid and oil pipe
Coefficient；
According between the air in the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, oil jacket annular space and sleeve pipe
Thermal convection current thermal resistance determine the 4th thermal resistivity.
4. method according to claim 1, is characterized in that, determines the heat of described Equations of The Second Kind section according to described data information
Resistance includes：
Determine the heat on stratum according to described formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well production time, wellbore radius
Resistance；
Determine the thermal resistance of cement sheath according to described sleeve outer wall radius, cement sheath thermal conductivity factor, wellbore radius；
Determine the thermal resistance of casing wall according to described sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall radius；
According to described crude oil thermal conductivity factor, the thermal conductivity factor of water, liquid aqueous rate, oil-pipe external wall radius, internal surface of sleeve pipe radius
Determine the thermal convection current liquid thermal resistance between liquid and sleeve pipe；
Thermal resistance between oil pipe inside and outside wall is determined according to described oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall radius；
According to described crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem exterior radius, tube inner wall half
Footpath determines the thermal convection current liquid thermal resistance between liquid and oil pipe；
Determined in described hollow stem according to the thermal conductivity factor of described relative oil density, hollow stem inwall radius, hollow stem
Thermal convection current liquid thermal resistance between outer wall；
Thermal convection current thermal resistance between being determined in hot fluid and hollow stem according to the thermal conductivity factor of described water.
5. method according to claim 4, is characterized in that, when described wellbore section, hollow bar segment are Equations of The Second Kind section, root
Determine each wellbore section according to described thermal resistance, the thermal resistivity of hollow bar segment includes：
According to the thermal convection current liquid heat between the thermal convection current liquid thermal resistance between described liquid and oil pipe, the inside and outside wall of hollow stem
Thermal convection current thermal resistance between in resistance, hot fluid and hollow stem determines the first thermal resistivity；
According to the thermal convection current liquid between the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, liquid and sleeve pipe
Thermal convection current liquid thermal resistance between thermal resistance between thermal resistance, oil pipe inside and outside wall, liquid and oil pipe determines the second thermal resistivity；
3rd thermal resistance is determined according to the thermal convection current liquid thermal resistance between the thermal resistance between described oil pipe inside and outside wall, liquid and oil pipe
Coefficient；
According to the thermal convection current liquid between the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall, liquid and sleeve pipe
Thermal resistance determines the 4th thermal resistivity.
6. a kind of system being determined temperature in wellbore field distribution based on hollow stem, be is characterized in that, described system includes：
Data information acquisition device, for obtaining the data information related to hollow stem and pit shaft, described data information bag
Include formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well production time, wellbore radius, sleeve outer wall radius, cement sheath heat conduction
Coefficient, sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, annular space radiation heat transfer coefficient, annular space free convection heat transfer coefficient, oil-pipe external wall
Outside radius, oil pipe thermal conductivity factor, tube inner wall radius, crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow stem
Wall radius, relative oil density, hollow stem inwall radius, the thermal conductivity factor of hollow stem；
Step size settings device, for the depth-set step-length according to the height in hydrodynamic face and pit shaft in described pit shaft；
Sectioning, for being divided into multiple wellbore sections, hollow bar segment according to described step-length by described pit shaft and hollow stem；
Temperature determines device, for determining the temperature of liquid in the plurality of wellbore section, many respectively according to described data information
The temperature of fluid in individual hollow bar segment, described temperature determines that device includes thermal resistance determining module, for determining each well successively
Cylinder section, the thermal resistance of hollow bar segment；Thermal resistivity determining module, for determining each wellbore section, hollow stem according to described thermal resistance
The thermal resistivity of section；Integral constant determining module, for determining integral constant according to wellhead temperature, bottom hole temperature (BHT)；Temperature determines
Module, for according to described integral constant and described thermal resistivity determine the temperature of liquid in each wellbore section, each
The temperature of fluid in hollow bar segment；
Thermo parameters method determines device, in the temperature for fluid in the plurality of hollow bar segment and described multiple wellbore sections
The temperature composition temperature in wellbore field distribution of liquid；
Wherein, described thermal resistance determining module includes judging unit, for judging each described wellbore section, hollow bar segment successively
Whether it is between the well head corresponding with described pit shaft of described hydrodynamic face；First kind section arranging unit, for when described judgement
When module is judged as YES, by the multiple wellbore sections being between the well head corresponding with described pit shaft of described hydrodynamic face, hollow bar segment
It is set to first kind section；Equations of The Second Kind section arranging unit, for when judge module is judged as NO, will be in described hydrodynamic face with described
Multiple wellbore sections between the corresponding shaft bottom of pit shaft, hollow bar segment are set to Equations of The Second Kind section；First thermal resistance determining module, for basis
Described data information determines the thermal resistance of described first kind section；Second thermal resistance determining module, for according to described data information
Determine the thermal resistance of described Equations of The Second Kind section.
7. system according to claim 6, is characterized in that, the first described thermal resistance determining module includes：
First thermal resistance determining unit, for during according to described formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well production
Between, wellbore radius determine the thermal resistance on stratum；
Second thermal resistance determining unit, for determining water according to described sleeve outer wall radius, cement sheath thermal conductivity factor, wellbore radius
The thermal resistance of mud ring；
3rd thermal resistance determining unit, for determining according to described sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall radius
The thermal resistance of casing wall；
4th thermal resistance determining unit, for according to described annular space radiation heat transfer coefficient, annular space free convection heat transfer coefficient, sleeve pipe
Inwall radius determines the thermal convection current thermal resistance between air and sleeve pipe in oil jacket annular space；
5th thermal resistance determining unit, for determining according to described oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall radius
Thermal resistance between oil pipe inside and outside wall；
6th thermal resistance determining unit, for according to described crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow
Bar exterior radius, tube inner wall radius determine the thermal convection current liquid thermal resistance between liquid and oil pipe；
7th thermal resistance determining unit, for according to described relative oil density, hollow stem inwall radius, hollow stem heat conduction system
Number determines the thermal convection current liquid thermal resistance between the inside and outside wall of described hollow stem；
8th thermal resistance determining unit is right for the heat between being determined in hot fluid and hollow stem according to the thermal conductivity factor of described water
Stream thermal resistance.
8. system according to claim 7, is characterized in that, when described wellbore section, hollow bar segment are first kind section, institute
The thermal resistivity determining module stated includes：
First thermal resistivity determining unit, for according to the thermal convection current liquid thermal resistance between described liquid and oil pipe, hollow stem
Inside and outside wall between thermal convection current liquid thermal resistance, in hot fluid and hollow stem between thermal convection current thermal resistance determine the first thermal resistance system
Number；
Second thermal resistivity determining unit, for according to the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall,
Between thermal resistance between thermal convection current thermal resistance between air in oil jacket annular space and sleeve pipe, oil pipe inside and outside wall, liquid and oil pipe
Thermal convection current liquid thermal resistance determines the second thermal resistivity；
3rd thermal resistivity determining unit, for according between the thermal resistance between described oil pipe inside and outside wall, liquid and oil pipe
Thermal convection current liquid thermal resistance determines the 3rd thermal resistivity；
4th thermal resistivity determining unit, for according to the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall,
Thermal convection current thermal resistance between air in oil jacket annular space and sleeve pipe determines the 4th thermal resistivity.
9. system according to claim 6, is characterized in that, the second described thermal resistance determining module includes：
First thermal resistance determining unit, for during according to described formation thermal conductivity, the average coefficient of heat transfer in stratum, oil well production
Between, wellbore radius determine the thermal resistance on stratum；
Second thermal resistance determining unit, for determining water according to described sleeve outer wall radius, cement sheath thermal conductivity factor, wellbore radius
The thermal resistance of mud ring；
3rd thermal resistance determining unit, for determining according to described sleeve pipe thermal conductivity factor, internal surface of sleeve pipe radius, sleeve outer wall radius
The thermal resistance of casing wall；
4th thermal resistance determining unit, for according to described crude oil thermal conductivity factor, the thermal conductivity factor of water, liquid aqueous rate, oil pipe
Exterior radius, internal surface of sleeve pipe radius determine the thermal convection current liquid thermal resistance between liquid and sleeve pipe；
5th thermal resistance determining unit, for determining according to described oil-pipe external wall radius, oil pipe thermal conductivity factor, tube inner wall radius
Thermal resistance between oil pipe inside and outside wall；
6th thermal resistance determining unit, for according to described crude oil thermal conductivity factor, liquid aqueous rate, the thermal conductivity factor of water, hollow
Bar exterior radius, tube inner wall radius determine the thermal convection current liquid thermal resistance between liquid and oil pipe；
7th thermal resistance determining unit, for according to described relative oil density, hollow stem inwall radius, hollow stem heat conduction system
Number determines the thermal convection current liquid thermal resistance between the inside and outside wall of described hollow stem；
8th thermal resistance determining unit is right for the heat between being determined in hot fluid and hollow stem according to the thermal conductivity factor of described water
Stream thermal resistance.
10. system according to claim 8, is characterized in that, when described wellbore section, hollow bar segment are Equations of The Second Kind section,
Described thermal resistivity determining module includes：
Thermal resistivity first determining unit, for according to the thermal convection current liquid thermal resistance between described liquid and oil pipe, hollow stem
Inside and outside wall between thermal convection current liquid thermal resistance, in hot fluid and hollow stem between thermal convection current thermal resistance determine the first thermal resistance system
Number；
Thermal resistivity second determining unit, for according to the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall,
Thermal convection current liquid between thermal resistance between thermal convection current liquid thermal resistance between liquid and sleeve pipe, oil pipe inside and outside wall, liquid and oil pipe
Body heat resistance determines the second thermal resistivity；
Thermal resistivity the 3rd determining unit, for according between the thermal resistance between described oil pipe inside and outside wall, liquid and oil pipe
Thermal convection current liquid thermal resistance determines the 3rd thermal resistivity；
Thermal resistivity the 4th determining unit, for according to the thermal resistance on described stratum, the thermal resistance of cement sheath, the thermal resistance of casing wall,
Thermal convection current liquid thermal resistance between liquid and sleeve pipe determines the 4th thermal resistivity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410534832.1A CN104389589B (en) | 2014-10-11 | 2014-10-11 | Method and system of determining shaft temperature field distribution based on hollow-rod |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410534832.1A CN104389589B (en) | 2014-10-11 | 2014-10-11 | Method and system of determining shaft temperature field distribution based on hollow-rod |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104389589A CN104389589A (en) | 2015-03-04 |
| CN104389589B true CN104389589B (en) | 2017-02-15 |
Family
ID=52607538
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410534832.1A Active CN104389589B (en) | 2014-10-11 | 2014-10-11 | Method and system of determining shaft temperature field distribution based on hollow-rod |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104389589B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105134179B (en) * | 2015-08-21 | 2018-02-02 | 中国石油天然气股份有限公司 | A kind of computational methods of natural gas well wellbore pressure and the distributed data of temperature |
| CN105422084B (en) * | 2015-11-26 | 2018-07-13 | 中国石油天然气股份有限公司 | A kind of the Wellbore Temperature Field acquisition methods and device of hot water circulation heated viscosity reduction technique |
| CN106021958B (en) * | 2016-06-14 | 2019-04-12 | 中国石油天然气股份有限公司 | A kind of determination method and device of gas injection well packer upper and lower ends temperature |
| CN106401568A (en) * | 2016-09-23 | 2017-02-15 | 中国石油天然气股份有限公司 | Method and device for determining temperature distribution of heat-insulation oil pipe wellhole |
| CN107842361B (en) * | 2017-10-30 | 2021-03-23 | 西南石油大学 | Method for measuring original formation temperature, empty wellbore static temperature, annulus static temperature and annulus dynamic temperature |
| CN108896807B (en) * | 2018-06-01 | 2021-08-03 | 清华大学深圳研究生院 | High-voltage bushing |
| CN109033012B (en) * | 2018-06-28 | 2023-01-06 | 中国石油天然气股份有限公司 | Method and device for determining temperature field of hollow sucker rod hot water injection circulation shaft |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2319834C1 (en) * | 2006-06-29 | 2008-03-20 | Общество с ограниченной ответственностью "ПермНИПИнефть" | Method for fluid flow control in well production strings |
| CN102682195A (en) * | 2012-03-09 | 2012-09-19 | 中国石油大学(华东) | Mine shaft temperature calculation method for semi-submersible type platform transient drilling well |
| CN103132990A (en) * | 2011-11-25 | 2013-06-05 | 普拉德研究及开发股份有限公司 | Dynamic prediction of downhole temperature distributions |
| CN203248144U (en) * | 2013-05-28 | 2013-10-23 | 中国石油天然气股份有限公司 | Digital temperature and pressure measuring apparatus for oil field well mouth fluid |
| CN103590818A (en) * | 2013-10-21 | 2014-02-19 | 中国石油天然气股份有限公司 | Shaft interior temperature distribution semi-analytical determination method and device |
| CN103775058A (en) * | 2013-12-31 | 2014-05-07 | 中国石油天然气股份有限公司 | Shaft heat loss determining method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8714239B2 (en) * | 2011-04-27 | 2014-05-06 | Luis Phillipe TOSI | Flow-induced electrostatic power generator for downhole use in oil and gas wells |
-
2014
- 2014-10-11 CN CN201410534832.1A patent/CN104389589B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2319834C1 (en) * | 2006-06-29 | 2008-03-20 | Общество с ограниченной ответственностью "ПермНИПИнефть" | Method for fluid flow control in well production strings |
| CN103132990A (en) * | 2011-11-25 | 2013-06-05 | 普拉德研究及开发股份有限公司 | Dynamic prediction of downhole temperature distributions |
| CN102682195A (en) * | 2012-03-09 | 2012-09-19 | 中国石油大学(华东) | Mine shaft temperature calculation method for semi-submersible type platform transient drilling well |
| CN203248144U (en) * | 2013-05-28 | 2013-10-23 | 中国石油天然气股份有限公司 | Digital temperature and pressure measuring apparatus for oil field well mouth fluid |
| CN103590818A (en) * | 2013-10-21 | 2014-02-19 | 中国石油天然气股份有限公司 | Shaft interior temperature distribution semi-analytical determination method and device |
| CN103775058A (en) * | 2013-12-31 | 2014-05-07 | 中国石油天然气股份有限公司 | Shaft heat loss determining method |
Non-Patent Citations (1)
| Title |
|---|
| 电加热井的井筒温度场数学模型;汪泓;《油气井测试》;20030630;第12卷(第03期);1-3 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104389589A (en) | 2015-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104389589B (en) | Method and system of determining shaft temperature field distribution based on hollow-rod | |
| CN113826099B (en) | Flow simulation and transient well analysis method based on generalized pipe flow seepage coupling | |
| Li et al. | Experimental and numerical studies on heat transfer characteristics of vertical deep-buried U-bend pipe to supply heat in buildings with geothermal energy | |
| CN105840187B (en) | Compactness reservoir horizontal well staged fracturing Productivity | |
| CN104453861B (en) | A kind of determination method of high-pressure gas well well bore temperature distribution and system | |
| CN102392634B (en) | Measuring device and measuring method for well-cementing annular weight loss of cement slurry | |
| CN105156102B (en) | Bottom water reservoir water energy three-dimensional physical simulation device and method | |
| CN107145696A (en) | A kind of analogy method of coal bed gas above and below ground couple solution | |
| CN105134179B (en) | A kind of computational methods of natural gas well wellbore pressure and the distributed data of temperature | |
| CN105422084B (en) | A kind of the Wellbore Temperature Field acquisition methods and device of hot water circulation heated viscosity reduction technique | |
| Furui et al. | A comprehensive model of high-rate matrix acid stimulation for long horizontal wells in carbonate reservoirs | |
| CN102682195B (en) | Semisubmersible platform transient state bored shaft temperature computation method | |
| CN112541287A (en) | Loose sandstone fracturing filling sand control production increase and profile control integrated design method | |
| CN109033012B (en) | Method and device for determining temperature field of hollow sucker rod hot water injection circulation shaft | |
| CN109667564A (en) | A kind of determination method of offshore viscous oil field steam soak exploitation orientation well capacity | |
| Du et al. | Experimental and numerical simulation research on heat transfer performance of coaxial casing heat exchanger in 3500m-deep geothermal well in Weihe Basin | |
| Sliwa et al. | Applicability of Borehole R-1 as BHE for Heating of a Gas Well | |
| CN108106687A (en) | A kind of bedrock underground water drift net containing soft interlayer probes into method and double-capsule water-stop | |
| Liao et al. | A novel 2-D equivalent numerical model of helix energy pile based on heat transfer characteristics of internal heat convection | |
| CN105184061B (en) | The method for numerical simulation of gas-producing well temperature, pressure distribution | |
| CN109930999A (en) | A kind of failure analysis and its manufacturing process | |
| CN206477824U (en) | The equipment of exploitation of geothermal energy | |
| CN104373117B (en) | Method and the system of temperature in wellbore field distribution is determined based on hollow stem electrical heating | |
| CN112347599B (en) | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment | |
| Meng et al. | Numerical simulation of fracture temperature field distribution during oil and gas reservoir hydraulic fracturing based on unsteady wellbore temperature field model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |