CN105370255B - Method and device for determining temperature distribution of gas injection and electric ignition shaft of in-situ combustion cage - Google Patents

Method and device for determining temperature distribution of gas injection and electric ignition shaft of in-situ combustion cage Download PDF

Info

Publication number
CN105370255B
CN105370255B CN201510864628.0A CN201510864628A CN105370255B CN 105370255 B CN105370255 B CN 105370255B CN 201510864628 A CN201510864628 A CN 201510864628A CN 105370255 B CN105370255 B CN 105370255B
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
represent
oil pipe
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
Application number
CN201510864628.0A
Other languages
Chinese (zh)
Other versions
CN105370255A (en
Inventor
马振
张洪君
刘利
李树全
曲绍刚
孙振宇
屈振哲
王智博
杨宝春
于晓聪
张成博
方梁锋
刘红芹
汪盈盈
关志忠
王河
钟满发
徐东
夏进军
吴超
崔冠麟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201510864628.0A priority Critical patent/CN105370255B/en
Publication of CN105370255A publication Critical patent/CN105370255A/en
Application granted granted Critical
Publication of CN105370255B publication Critical patent/CN105370255B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)

Abstract

the invention discloses a method and a device for determining temperature distribution of a gas injection and electric ignition shaft of an in-situ combustion cage system, wherein the method comprises the steps of 1, dividing a shaft into a plurality of shaft units dl in the axial direction, enabling L to be 0 and k to be 1, 2, calculating the temperature of oil pipe air heated by an electric igniter, 3, calculating formation thermal resistance, cement ring thermal resistance, sleeve wall thermal resistance, thermal resistance between oil sleeve annulus air and a sleeve, oil pipe wall thermal resistance and air thermal resistance in the oil pipe, sequentially enabling the shaft to comprise the oil pipe, the sleeve and the cement ring from inside to outside, enabling the position in the oil sleeve annulus corresponding to the top of an ignition gun to be connected with a packer, dividing the oil pipe into an upper section and a lower section, 4, calculating the total thermal resistance of the shaft in the radial direction, 5, calculating the heat loss of the shaft in the radial direction, 6, calculating the air temperature of the oil pipe, 7, enabling L to be L + dl and k to be k +1, repeatedly executing 3-6 to carry out iterative calculation until L is not less than the end, obtaining an oil pipe temperature curve.

Description

The determination method and device of the general gas injection electric ignition well bore temperature distribution of combustion in situ
Technical field
The present invention relates to combustion in situ field, more particularly to a kind of general gas injection electric ignition well bore temperature distribution of combustion in situ Determination method and device.
Background technology
Domestic combustion in situ electric ignition technique is largely burnt for individual layer at present, the general gas injection of individual layer igniting.Casing programme As shown in figure 1, include successively from the inside to the outside：Oil pipe, oil jacket annular space, sleeve pipe, cement sheath, stratum.There is electric igniter in oil pipe, by It is very big in the caloric value of burning torch, returned to not allow as much as possible on annular space heat, the oil pipe company of (in figure B at) at the top of burning torch Packer is connected to, by oil jacket annular isolation.In Fig. 1, AB sections are the cable length of igniter, and the caloric value of cable is comparatively It is smaller；BC sections are the length of burning torch, and burning torch about length is 50m.AB sections are plain tubing, and BC sections have three kinds of possibility： 1. using instlated tubular, thermal loss is avoided as much as, this mode can farthest improve the temperature of exit air, But this tubular column structure connect it is cumbersome because plain tubing connects with instlated tubular needs reducing；2. screen casing is used, This mode connects more convenient, but heat loss is maximum；3. using plain tubing, make upper-lower section oil pipe easy to connect, The heat loss of this mode is larger.General gas injection electric ignition can be realized by said structure.
It is the key for realizing above-mentioned general gas injection electric ignition to calculate well bore temperature distribution, still, not yet proposes to be directed at present The determination method of the gas injection well well bore temperature distribution of above-mentioned general gas injection electric ignition.
The content of the invention
The invention provides a kind of determination method and device of the general gas injection electric ignition well bore temperature distribution of combustion in situ, with At least solve the problems, such as that there has been no general gas injection ignition process well-sinking temperature field at present to determine method.
According to an aspect of the invention, there is provided a kind of general gas injection electric ignition well bore temperature distribution of combustion in situ is really Determine method, including：Step 1, pit shaft is divided into multiple pit shaft units in the axial direction, the length of each pit shaft unit is dl, makes l =0, k=1, wherein, l represents the length currently calculated, and k represents iterations；Step 2, the air warp for being injected into oil pipe is calculated The temperature T crossed after electric igniter heatings；Step 3, the thermal resistance R on stratum is calculated respectively1, cement sheath thermal resistance R2, sleeve pipe inside and outside wall Between thermal resistance R3, thermal resistance R between air and sleeve pipe in oil jacket annular space4, thermal resistance R between oil pipe inside and outside wall5And oil pipe The thermal convection current thermal resistance R of interior air6；Wherein, the pit shaft radially includes successively from the inside to the outside：Oil pipe, sleeve pipe and cement sheath, oil The position that the collar corresponds at the top of burning torch in the air is connected with packer, and using the packer as boundary, the oil pipe is divided into epimere Oil pipe and hypomere oil pipe, the hypomere oil pipe are instlated tubular, screen casing or plain tubing, are stratum outside the pit shaft；Step 4, According to R1To R6Calculate the entire thermal resistance of pit shaft diametrically；Step 5, according to the temperature Ts, the entire thermal resistance and formation temperature, Calculate the heat loss of pit shaft diametrically；Step 6, according to the temperature Ts, the heat loss and the electric igniter power, Calculate the air themperature of oil pipe；Step 7, l=l+dl, k=k+1 are made, according to the change of formation temperature, repeats above-mentioned steps 3 to step 6, is iterated calculating, and until l >=L, then iteration terminates, and obtains the air themperature distribution curve of the oil pipe, its In, L represents the total length of oil pipe.
According to another aspect of the present invention, there is provided a kind of general gas injection electric ignition well bore temperature distribution of combustion in situ Determining device, including：Division unit, for pit shaft to be divided into multiple pit shaft units, the length of each pit shaft unit in the axial direction Spend for dl, make l=0, k=1, wherein, l represents the length currently calculated, and k represents iterations；First computing unit, based on Calculate temperature T of the air for being injected into oil pipe after electric igniter heatss；Second computing unit, for calculating stratum respectively Thermal resistance R1, cement sheath thermal resistance R2, thermal resistance R between sleeve pipe inside and outside wall3, thermal resistance between air and sleeve pipe in oil jacket annular space R4, thermal resistance R between oil pipe inside and outside wall5And the thermal convection current thermal resistance R of oily inner air tube6；Wherein, the pit shaft is radially from interior To including successively outside：Oil pipe, sleeve pipe and cement sheath, the position corresponded at the top of burning torch in oil jacket annular space are connected with packer, Using the packer as boundary, the oil pipe is divided into epimere oil pipe and hypomere oil pipe, and the hypomere oil pipe is instlated tubular, screen casing or general Oil pipe, the pit shaft outside is stratum；3rd computing unit, for according to R1To R6Calculate the entire thermal resistance of pit shaft diametrically； 4th computing unit, for according to the temperature Ts, the entire thermal resistance and formation temperature, calculate the heat waste of pit shaft diametrically Lose；5th computing unit, for according to the temperature Ts, the heat loss and the electric igniter power, calculate oil pipe Air themperature；Unit is iterated to calculate, for making l=l+dl, k=k+1, according to the change of formation temperature, is calculated using second single Member is iterated calculatings to the 5th computing unit, and until l >=L, then iteration terminates, and obtains the air themperature distribution song of the oil pipe Line, wherein, L represents the total length of oil pipe.
By the determination method and device of the general gas injection electric ignition well bore temperature distribution of combustion in situ of the present invention, synthesis is examined The radially change of many factors along well depth such as heat transfer and stratum thermophysical property of well bore and oil pipe rod structure, pit shaft is considered, by pit shaft It is divided into some sections, obtains the physical parameter (thermal resistance, heat transfer coefficient) of correspondent section, part physical parameter is the function of temperature, is used Solution by iterative method, oil pipe Temperature Distribution is calculated.In the case of general gas injection ignition process can accurately be calculated, shape is arbitrarily flowed The Temperature Distribution of condition, any time along gas injection well shaft.Meanwhile calculating process is simple and convenient, there is higher precision, iteration time Number is low, and computational efficiency is high, has extraordinary stability and convergence.According to the Temperature Distribution of pit shaft, can effectively predict Up to the air themperature of oil reservoir, to adjust gas injection rate and igniter power, and then ensure the smooth implementation of ignition process.
Brief description of the drawings
Accompanying drawing described herein is used for providing a further understanding of the present invention, forms the part of the application, this hair Bright schematic description and description is used to explain the present invention, does not form limitation of the invention.In the accompanying drawings：
Fig. 1 is the structural representation of the general gas injection electric ignition of combustion in situ of the embodiment of the present invention；
Fig. 2 is the flow of the determination method of the general gas injection electric ignition well bore temperature distribution of combustion in situ of the embodiment of the present invention Figure；
Fig. 3 is the structure of the determining device of the general gas injection electric ignition well bore temperature distribution of combustion in situ of the embodiment of the present invention Block diagram；
Fig. 4 is that the BC sections of the embodiment of the present invention are the temperature distributing curve diagrams of oil pipe in the case of instlated tubular.
Embodiment
With reference to 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 Ground describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.Based on this The embodiment of invention, the every other implementation that those of ordinary skill in the art are obtained under the premise of creative work is not made Example, belongs to protection scope of the present invention.
The embodiments of the invention provide a kind of determination method of the general gas injection electric ignition well bore temperature distribution of combustion in situ, figure 2 be the flow chart of the determination method of the general gas injection electric ignition well bore temperature distribution of combustion in situ of the embodiment of the present invention, such as Fig. 2 institutes Show, this method includes steps S201 to step S207.
Step S201, pit shaft is divided into multiple pit shaft units in the axial direction, the length of each pit shaft unit is dl, makes l =0, k=1, wherein, l represents the length currently calculated, and k represents iterations.
Step S202, calculate temperature T of the air for being injected into oil pipe after electric igniter heatss。
Step S203, the thermal resistance R on stratum is calculated respectively1, cement sheath thermal resistance R2, thermal resistance R between sleeve pipe inside and outside wall3, oil Thermal resistance R between the collar aerial air and sleeve pipe4, thermal resistance R between oil pipe inside and outside wall5And the thermal convection current of oily inner air tube Thermal resistance R6；Wherein, pit shaft radially includes successively from the inside to the outside：Oil pipe, sleeve pipe and cement sheath, igniting is corresponded in oil jacket annular space Position at the top of rifle is connected with packer, and using packer as boundary, oil pipe is divided into epimere oil pipe and hypomere oil pipe, hypomere oil pipe be every Heat pipe, screen casing or plain tubing, pit shaft outside is stratum.
Step S204, according to R1To R6Calculate the entire thermal resistance of pit shaft diametrically.
Step S205, according to temperature Ts, entire thermal resistance and formation temperature, calculate the heat loss of pit shaft diametrically.
Step S206, according to temperature Ts, heat loss and electric igniter power, calculate the air themperature of oil pipe.
Step S207, makes l=l+dl, k=k+1, according to the change of formation temperature, repeats above-mentioned steps S203 to step Rapid S206, calculating is iterated, until l >=L, then iteration terminates, and obtains the air themperature distribution curve of oil pipe, wherein, L is represented The total length of oil pipe.
By the above method, well bore and oil pipe rod structure, pit shaft radially heat transfer and stratum thermophysical property etc. are considered Change of many factors along well depth, pit shaft is divided into some sections, obtains the physical parameter of correspondent section, part physical parameter is temperature Function, using solution by iterative method, the distribution of oil pipe air themperature is calculated.This method can accurately calculate general gas injection point firer In the case of skill, the air themperature distribution of any flow condition, any time along gas injection well shaft.Meanwhile this method is simple and convenient, With higher precision, iterations is low, and computational efficiency is high, has extraordinary stability and convergence.According to the temperature of pit shaft Degree distribution, can effectively predict the air themperature up to oil reservoir, to adjust gas injection rate and igniter power, and then ensure point firer The smooth implementation of skill.
Main assumption condition in the embodiment of the present invention is：
(1) for fluid flow state to stablize one-way flow, fluid is gas single-phase flow；
(2) heat transfer is steady heat transfer in pit shaft；
(3) stratum heat transfer is unsteady heat transfer, and obeys Ramey non dimensional time function；
(4) casing programme is as shown in Figure 1：Oil pipe-oil jacket annular space-sleeve pipe-cement sheath-stratum；
(5) heat loss in pit shaft and surrounding formation is radial direction, while is also contemplated for biography of the air flow along well depth direction Heat；
(6) formation temperature presses linear change, it is known that geothermal gradient and surface temperature；
(7) tubing and casing is concentric.
Casing programme as shown in Figure 1, it is the origin of coordinates to take well head, straight down for just.
In one embodiment, below equation can be used to calculate the thermal resistance R on stratum1：
Wherein, KeFormation thermal conductivity is represented, unit is W/ (mK)；A represents the average coefficient of heat transfer in stratum, unit m2/ d；T represents the oil well production time；rhRepresent wellbore radius (i.e. distance of the gas injection well axis to cement sheath outer wall), unit m.
In one embodiment, below equation can be used to calculate the thermal resistance R of cement sheath2：
Wherein, KcemCement sheath thermal conductivity factor is represented, unit is W/ (mK)；rhRepresent wellbore radius, unit m；rcoTable Show sleeve outer wall radius, unit m.
In one embodiment, below equation can be used to calculate the thermal resistance R between sleeve pipe inside and outside wall3：
Wherein, KcasSleeve pipe thermal conductivity factor is represented, unit is W/ (mK)；rciRepresent internal surface of sleeve pipe radius, unit m； rcoRepresent sleeve outer wall radius, unit m.
In one embodiment, below equation can be used to calculate the thermal resistance between air and sleeve pipe in oil jacket annular space R4：
Wherein, hc1The free convection heat transfer coefficient of air in oil jacket annular space is represented, unit is W/ (m2·K)；hr1Represent oil The heat radiation heat transfer coefficient of the aerial air of the collar, unit are W/ (m2·K)；rciRepresent internal surface of sleeve pipe radius.
Heat radiation heat transfer coefficient h is calculated using below equationr1：
Wherein, δ represents Stefan-Boltzmann (this special fence-Boltzmann) constant, and value is 2.189 × 10-8W/ (m2·K)；FtciRepresent oil-pipe external wall surface to internal surface of sleeve pipe surface emissivity coefficient of efficiency；TtoRepresent oil-pipe external wall temperature；TciTable Show internal surface of sleeve pipe temperature；εoRepresent oil-pipe external wall blackness；εciRepresent internal surface of sleeve pipe blackness；rtoRepresent oil-pipe external wall radius；
Free convection heat transfer coefficient h is calculated using below equationc1：
Wherein, GrRepresent Grashof numbers (grashof number)；PrRepresent Prandtl numbers (Prandtl number)；KhaRepresent oil jacket The thermal conductivity factor of the air of annular space, unit are W/ (mK)；G represents acceleration of gravity, unit m/s2；ρanRepresent oil jacket annular space Air in mean temperature TanUnder density, unit kg/m3；UanRepresent the air of oil jacket annular space in mean temperature TanUnder Viscosity, unit mPas；CanRepresent the air of oil jacket annular space in mean temperature TanUnder thermal capacitance, unit is J (m3·K)；β tables Show the thermal cubic expansion coefficient of air in oil jacket annular space, be a constant, value can be 1.78 × 10-3。
The thermal resistance of oil pipe can be divided to two sections to be calculated, as shown in figure 1, using B point packers as boundary, AB sections (i.e. epimere oil Pipe) it is plain tubing, BC sections (i.e. hypomere oil pipe) have three kinds of situations：Instlated tubular, screen casing or plain tubing.Illustrate its calculating below Process.
In one embodiment, the thermal resistance R between oil pipe inside and outside wall is calculated5As shown in formula (12) to (15).
For epimere oil pipe, its thermal resistance R is calculated using below equation5：
Wherein, KtubOil pipe thermal conductivity factor is represented, unit is W/ (mK)；rtoRepresent oil-pipe external wall radius, unit m； rtiRepresent tube inner wall radius, unit m.
If hypomere oil pipe is instlated tubular, its thermal resistance R is calculated using below equation5：
Wherein, KtubOil pipe thermal conductivity factor is represented, unit is W/ (mK)；roRepresent instlated tubular outer tube outer wall radius, unit For m；riRepresent instlated tubular outer tube wall radius, unit m；KinsInstlated tubular thermal conductivity factor is represented, unit is W/ (mK)；rto1 Represent instlated tubular outer wall of inner tube radius, unit m；rti1Represent instlated tubular inner tube wall radius, unit m.
If hypomere oil pipe is screen casing, its thermal resistance R is calculated using below equation5：
Wherein, KsieScreen casing thermal conductivity factor is represented, unit is W/ (mK)；rto2Represent screen casing exterior radius, unit m； rti2Represent screen casing inwall radius, unit m.
If hypomere oil pipe is plain tubing, its thermal resistance calculation is identical with epimere oil pipe, specifically, using below equation its Thermal resistance R5：
Wherein, KtubOil pipe thermal conductivity factor is represented, unit is W/ (mK)；rtoRepresent oil-pipe external wall radius, unit m； rtiRepresent tube inner wall radius, unit m.
In one embodiment, below equation can be used to calculate the thermal convection current thermal resistance R of oily inner air tube6：
Wherein, hfThe thermal conductivity factor coefficient of oily inner air tube is represented, value is 0.05W/ (mK)；rtiRepresent tube inner wall Radius；rtoRepresent oil-pipe external wall radius.
In one embodiment, according to R1To R6The entire thermal resistance of pit shaft diametrically is calculated using below equation：
R=R1+R2+R3+R4+R5+R6 (17)
In Fig. 1, pit shaft is divided into several pit shaft units, initial temperature (the i.e. oily inner air tube of oil pipe in the axial direction Temperature after electric igniter heats) it is the T being calculated by formula (22) and (23)s, main heat loss is radially On heat loss.
In one embodiment, according to temperature Ts, entire thermal resistance and formation temperature, calculate the heat loss of pit shaft diametrically, Including：According to law of conservation of energy, heat loss is calculated using below equation：
Wherein, Q represents pit shaft unit radial heat loss, unit W；TeFormation temperature is represented, unit is DEG C；R represents well Cylinder unit radial entire thermal resistance.
In one embodiment, according to temperature Ts, heat loss and electric igniter power, calculate the air themperature of oil pipe, Including：
Air themperature in epimere oil pipe (i.e. AB sections) is calculated using below equation：
CmTs- Q/1000+0.4P=CmTs' (19)
Air themperature in hypomere oil pipe (i.e. BC sections) is calculated using below equation：
CmTs- Q/1000+0.6P=CmTs' (20)
Wherein, Ts' the air themperature after changing in oil pipe is represented, unit is DEG C；C represents the specific heat capacity of air, and value is 1.0069kJ/ (kg ﹒ DEG C)；M represents the mass flow of air, unit kg/s；P represents the power of electric igniter；Q represents pit shaft Unit radial heat loss.
In one embodiment, can with step S207 use below equation calculate formation temperature change：
Te=Tins+αdl (21)
Wherein, TinsSurface temperature is represented, unit is DEG C；α represent geothermal gradient, unit for DEG C/m；TeRepresent stratum temperature Degree, unit are DEG C.
In one embodiment, it is contemplated that electric igniter cable probably has 40% along journey heat loss, and calculating is injected into oil pipe Temperature T of the air after electric igniter heatss, including：
Temperature T of the air of epimere oil pipe after electric igniter heats is calculated using below equations：
CmT+0.4P=CmTs (22)
The air of hypomere oil pipe is calculated by electric igniter temperature after heating change T using below equations：
CmT+0.6P=CmTs (23)
Wherein, T represents the initial temperature of air, and unit is DEG C；C represents the specific heat capacity of air；M represents the quality stream of air Amount；P represents the power of electric igniter.
Based on same inventive concept, the embodiment of the present invention additionally provides a kind of general gas injection electric ignition pit shaft temperature of combustion in situ The determining device of distribution is spent, can be used for realizing the method described by above-described embodiment, part is repeated and repeats no more.It is following to be made , term " unit " can realize the combination of the software and/or hardware of predetermined function.Although described by following examples is System is preferably realized with software, but hardware, or software and hardware combination realization and may and be contemplated.
Fig. 3 is the structure of the determining device of the general gas injection electric ignition well bore temperature distribution of combustion in situ of the embodiment of the present invention Block diagram, as shown in figure 3, the device includes：Division unit 31, the first computing unit 32, the second computing unit the 33, the 3rd calculate single First 34, the 4th computing unit 35, the 5th computing unit 36 and iterative calculation unit 37.The structure is specifically described below.
Division unit 31, for pit shaft to be divided into multiple pit shaft units in the axial direction, the length of each pit shaft unit is Dl, l=0, k=1 are made, wherein, l represents the length currently calculated, and k represents iterations；
First computing unit 32, temperature T of the air of oil pipe after electric igniter heats is injected into for calculatings；
Second computing unit 33, for calculating the thermal resistance R on stratum respectively1, cement sheath thermal resistance R2, between sleeve pipe inside and outside wall Thermal resistance R3, thermal resistance R between air and sleeve pipe in oil jacket annular space4, thermal resistance R between oil pipe inside and outside wall5It is and empty in oil pipe The thermal convection current thermal resistance R of gas6；Wherein, pit shaft radially includes successively from the inside to the outside：Oil pipe, sleeve pipe and cement sheath, in oil jacket annular space Packer is connected with corresponding to the position at the top of burning torch, using packer as boundary, oil pipe is divided into epimere oil pipe and hypomere oil pipe, under Section oil pipe is instlated tubular, screen casing or plain tubing, is stratum outside pit shaft；
3rd computing unit 34, for according to R1To R6Calculate the entire thermal resistance of pit shaft diametrically；
4th computing unit 35, for according to temperature Ts, entire thermal resistance and formation temperature, calculate the heat waste of pit shaft diametrically Lose；
5th computing unit 36, for according to temperature Ts, heat loss and electric igniter power, calculate the Air Temperature of oil pipe Degree；
Unit 37 is iterated to calculate, for making l=l+dl, k=k+1, according to the change of formation temperature, is calculated using second single First 33 to the 5th computing units 36 are iterated calculatings, and until l >=L, then iteration terminates, and obtains the air themperature distribution song of oil pipe Line, wherein, L represents the total length of oil pipe.
By said apparatus, well bore and oil pipe rod structure, pit shaft radially heat transfer and stratum thermophysical property etc. are considered Change of many factors along well depth, pit shaft is divided into some sections, obtains the physical parameter (thermal resistance, heat transfer coefficient) of correspondent section, portion Divide the function that physical parameter is temperature, using solution by iterative method, the distribution of oil pipe air themperature is calculated.The device can be counted accurately In the case of calculating general gas injection ignition process, the air themperature distribution of any flow condition, any time along gas injection well shaft.Together When, calculating process is simple and convenient, has higher precision, and iterations is low, and computational efficiency is high, have extraordinary stability and Convergence.According to the Temperature Distribution of pit shaft, the air themperature up to oil reservoir can be effectively predicted, to adjust gas injection rate and igniter Power, and then ensure the smooth implementation of ignition process.
In one embodiment, the second computing unit 33 is specifically used for the thermal resistance R that stratum is calculated using below equation1：
Wherein, KeRepresent formation thermal conductivity；A represents the average coefficient of heat transfer in stratum；T represents the oil well production time；rhRepresent Wellbore radius.
In one embodiment, the second computing unit 33 is specifically used for the thermal resistance R that cement sheath is calculated using below equation2：
Wherein, KcemRepresent cement sheath thermal conductivity factor；rhRepresent wellbore radius；rcoRepresent sleeve outer wall radius.
In one embodiment, the second computing unit 33 is specifically used for using between below equation calculating sleeve pipe inside and outside wall Thermal resistance R3：
Wherein, KcasRepresent sleeve pipe thermal conductivity factor；rciRepresent internal surface of sleeve pipe radius；rcoRepresent sleeve outer wall radius.
In one embodiment, the second computing unit 33 is specifically used for calculating the air in oil jacket annular space using below equation Thermal resistance R between sleeve pipe4：
Wherein, hc1Represent the free convection heat transfer coefficient of air in oil jacket annular space；hr1Represent the heat of air in oil jacket annular space Radiation heat transfer coefficient；rciRepresent internal surface of sleeve pipe radius；
Heat radiation heat transfer coefficient h is calculated using below equationr1：
Wherein, δ represents Stefan-Boltzmann constants；FtciRepresent oil-pipe external wall surface to internal surface of sleeve pipe surface emissivity Coefficient of efficiency；TtoRepresent oil-pipe external wall temperature；TciRepresent internal surface of sleeve pipe temperature；εoRepresent oil-pipe external wall blackness；εciRepresent sleeve pipe Inwall blackness；rtoRepresent oil-pipe external wall radius；
Free convection heat transfer coefficient h is calculated using below equationc1：
Wherein, GrRepresent Grashof numbers；PrRepresent Prandtl numbers；KhaRepresent the thermal conductivity factor of the air of oil jacket annular space；g Represent acceleration of gravity；ρanRepresent the air of oil jacket annular space in mean temperature TanUnder density；UanRepresent the air of oil jacket annular space In mean temperature TanUnder viscosity；CanRepresent the air of oil jacket annular space in mean temperature TanUnder thermal capacitance；β represents oil jacket annular space The thermal cubic expansion coefficient of middle air.
In one embodiment, the second computing unit 33 is specifically used for：
For epimere oil pipe, its thermal resistance R is calculated using below equation5：
Wherein, KtubRepresent oil pipe thermal conductivity factor；rtoRepresent oil-pipe external wall radius；rtiRepresent tube inner wall radius；
If hypomere oil pipe is instlated tubular, its thermal resistance R is calculated using below equation5：
Wherein, KtubRepresent oil pipe thermal conductivity factor；roRepresent instlated tubular outer tube outer wall radius；riRepresent instlated tubular outer tube wall Radius；KinsRepresent instlated tubular thermal conductivity factor；rto1Represent instlated tubular outer wall of inner tube radius；rti1Represent instlated tubular inner tube wall half Footpath；
If hypomere oil pipe is screen casing, its thermal resistance R is calculated using below equation5：
Wherein, KsieRepresent screen casing thermal conductivity factor；rto2Represent screen casing exterior radius；rti2Represent screen casing inwall radius；
If hypomere oil pipe is plain tubing, using its thermal resistance of below equation R5：
Wherein, KtubRepresent oil pipe thermal conductivity factor；rtoRepresent oil-pipe external wall radius；rtiRepresent tube inner wall radius.
In one embodiment, the second computing unit 33 is specifically used for calculating the hot right of oily inner air tube using below equation Flow thermal resistance R6：
Wherein, hfRepresent the thermal conductivity factor coefficient of oily inner air tube；rtiRepresent tube inner wall radius；rtoRepresent oil-pipe external wall Radius.
In one embodiment, the 3rd computing unit 34 is specifically used for calculating pit shaft diametrically total using below equation Thermal resistance：
R=R1+R2+R3+R4+R5+R6。
In one embodiment, the 4th computing unit 35 is specifically used for according to law of conservation of energy, using below equation meter Calculate heat loss：
Wherein, Q represents pit shaft unit radial heat loss；TeRepresent formation temperature；R represents pit shaft unit radial entire thermal resistance.
In one embodiment, the 5th computing unit 36 is specifically used for：
Air themperature in epimere oil pipe is calculated using below equation：
CmTs- Q/1000+0.4P=CmTs',
Air themperature in hypomere oil pipe is calculated using below equation：
CmTs- Q/1000+0.6P=CmTs',
Wherein, Ts' represent the air themperature after changing in oil pipe；C represents the specific heat capacity of air；M represents the quality of air Flow；P represents the power of electric igniter；Q represents pit shaft unit radial heat loss.
In one embodiment, the change that unit 37 is specifically used for calculating formation temperature using below equation is iterated to calculate：
Te=Tins+ α dl,
Wherein, TinsRepresent surface temperature；α represents geothermal gradient；TeRepresent formation temperature.
In one embodiment, the first computing unit 32 is specifically used for the air warp that epimere oil pipe is calculated using below equation The temperature T crossed after electric igniter heatings：
CmT+0.4P=CmTs,
The air of hypomere oil pipe is calculated by electric igniter temperature after heating change T using below equations：
CmT+0.6P=CmTs,
Wherein, T represents the initial temperature of air；C represents the specific heat capacity of air；M represents the mass flow of air；P is represented The power of electric igniter.
Certainly, above-mentioned Module Division is a kind of signal division, and the invention is not limited in this.As long as the present invention can be realized Purpose Module Division, protection scope of the present invention all should be belonged to.
In order to be carried out more to the determination method and device of the general gas injection electric ignition well bore temperature distribution of above-mentioned combustion in situ Clearly explain, illustrated with reference to specific embodiment, however, it should be noted that the embodiment is merely to more Illustrate the present invention well, do not form and the present invention is improperly limited.
(1) pit shaft is divided into several pit shaft units in the vertical, each pit shaft element length is dl, is counted since well head Calculate, make l=0, k=1, air injects in well head, and the initial temperature of air is T0, make Ts=T0。
(2) R is calculated1, R2, R3, R5, R6, make R4=0 (due to R4Relevant with heat transfer coefficient, the temperature of heat transfer coefficient and pipe has Close, and do not know channel temp value initially, therefore, R is first set4It is worth 0), entire thermal resistance R to be calculated by formula (17).
(3) heat loss is calculated by formula (18)
(4) oil-pipe external wall temperature T is calculatedto=Ts-(R5+R6)×Q/dl。
(5) internal surface of sleeve pipe temperature T is calculatedci=Te+(R1+R2+R3)×Q/dl。
(6) R is calculated by formula (5)~(11)4。
(7) entire thermal resistance R is calculated again by formula (17).
(8) heat loss is calculated again
(9) the temperature T for calculating oil pipe air is segmented by formula (19)~(20)s'。
(10) k=k+1, l=l+dl are made, calculating formation temperature by formula (21) changes Te=Tins+ adl, return to the (3) step continues to iterate to calculate；If l >=L (oil pipe total length), iteration terminates.The temperature distribution history of oil pipe is obtained, with BC sections It is that oil pipe temperature curve is as shown in Figure 4 exemplified by instlated tubular.
In summary, the embodiment of the present invention is for there has been no general gas injection ignition process well-sinking temperature field determination side at present The problem of method, it is proposed that a kind of determination method and device of the general gas injection electric ignition well bore temperature distribution of combustion in situ, for fire During oil firing layer electric ignition, the calculating of general gas injection gas injection well well bore temperature distribution.Consider well bore and oil pipe rod structure, pit shaft The radially change of many factors along well depth such as heat transfer and stratum thermophysical property, is divided into some sections by pit shaft, obtains correspondent section Physical parameter, part physical parameter are the functions of temperature, and using solution by iterative method, the distribution of oil pipe air themperature is calculated.Root It is distributed according to the air themperature of pit shaft, can effectively predicts the air themperature up to oil reservoir, to adjust gas injection rate and igniter power, And then ensure the smooth implementation of ignition process.
The present invention establishes corresponding mathematical modeling using thermal conduction study method, and has carried out computer programming to this method. When establishing temperature distribution model, it is assumed that the heat transfer in pit shaft is steady state heat transfer, and the heat transfer in pit shaft surrounding formation is unstable state Heat transfer, heat loss radially is not only allowed for when calculating well bore temperature distribution, it is also considered that air flow is along well depth direction The influence conducted heat to well bore temperature distribution, it is segmented according to each section of tubing string difference of the tubular column structure of oil pipe, different situations difference Calculated.Calculating process is simple and convenient, has higher precision, and iterations is low, and computational efficiency is high, has extraordinary steady Qualitative and convergence, is more suitable for computer programming.In the case of general gas injection ignition process can accurately be calculated, shape is arbitrarily flowed The Temperature Distribution of condition, any time along gas injection well shaft
Any process or method described otherwise above description in flow chart or herein is construed as, and represents to include Module, fragment or the portion of the code of the executable instruction of one or more the step of being used to realize specific logical function or process Point, and the scope of the preferred embodiment of the present invention includes other realization, wherein can not press shown or discuss suitable Sequence, including according to involved function by it is basic simultaneously in the way of or in the opposite order, carry out perform function, this should be of the invention Embodiment person of ordinary skill in the field understood.
It should be appreciated that each several part of the present invention can be realized with hardware, software, firmware or combinations thereof.Above-mentioned In embodiment, software that multiple steps or method can be performed in memory and by suitable instruction execution system with storage Or firmware is realized.If, and in another embodiment, can be with well known in the art for example, realized with hardware Any one of row technology or their combination are realized：With the logic gates for realizing logic function to data-signal Discrete logic, have suitable combinational logic gate circuit application specific integrated circuit, programmable gate array (PGA), scene Programmable gate array (FPGA) etc..
Particular embodiments described above, the purpose of the present invention, technical scheme and beneficial effect are carried out further in detail Describe in detail it is bright, should be understood that the foregoing is only the present invention specific embodiment, the guarantor being not intended to limit the present invention Scope is protected, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., should be included in this Within the protection domain of invention.

Claims (13)

  1. A kind of 1. determination method of the general gas injection electric ignition well bore temperature distribution of combustion in situ, it is characterised in that including：
    Step 1, pit shaft is divided into multiple pit shaft units in the axial direction, the length of each pit shaft unit is dl, makes l=0, k= 1, wherein, l represents the length currently calculated, and k represents iterations；
    Step 2, temperature T of the air for being injected into oil pipe after electric igniter heats is calculateds；
    Step 3, the thermal resistance R on stratum is calculated respectively1, cement sheath thermal resistance R2, thermal resistance R between sleeve pipe inside and outside wall3, oil jacket annular space In air and sleeve pipe between thermal resistance R4, thermal resistance R between oil pipe inside and outside wall5And the thermal convection current thermal resistance R of oily inner air tube6； Wherein, the pit shaft radially includes successively from the inside to the outside：Oil pipe, sleeve pipe and cement sheath, correspond to burning torch in oil jacket annular space The position at top is connected with packer, and using the packer as boundary, the oil pipe is divided into epimere oil pipe and hypomere oil pipe, under described Section oil pipe is instlated tubular, screen casing or plain tubing, is stratum outside the pit shaft；
    Step 4, according to R1To R6Calculate the entire thermal resistance of pit shaft diametrically；
    Step 5, according to the temperature Ts, the entire thermal resistance and formation temperature, calculate the heat loss of pit shaft diametrically；
    Step 6, according to the temperature Ts, the heat loss and the electric igniter power, calculate the air themperature of oil pipe；
    Step 7, l=l+dl, k=k+1 are made, according to the change of formation temperature, above-mentioned steps 3 is repeated to step 6, is changed In generation, calculates, and until l >=L, then iteration terminates, and obtains the air themperature distribution curve of the oil pipe, wherein, L represents the total of oil pipe Length.
  2. 2. according to the method for claim 1, it is characterised in that the thermal resistance R on stratum is calculated using below equation1：
    <mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mi>e</mi> </msub> </mrow> </mfrac> <mo>,</mo> </mrow>
    <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>l</mi> <mi>n</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <msqrt> <mrow> <mi>a</mi> <mi>t</mi> </mrow> </msqrt> </mrow> <msub> <mi>r</mi> <mi>h</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mn>0.29</mn> <mo>,</mo> </mrow>
    Wherein, KeRepresent formation thermal conductivity；A represents the average coefficient of heat transfer in stratum；T represents the oil well production time；rhRepresent pit shaft Radius.
  3. 3. according to the method for claim 1, it is characterised in that the thermal resistance R of cement sheath is calculated using below equation2：
    <mrow> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>m</mi> </mrow> </msub> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <msub> <mi>r</mi> <mi>h</mi> </msub> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>o</mi> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>
    Wherein, KcemRepresent cement sheath thermal conductivity factor；rhRepresent wellbore radius；rcoRepresent sleeve outer wall radius.
  4. 4. according to the method for claim 1, it is characterised in that the thermal resistance between sleeve pipe inside and outside wall is calculated using below equation R3：
    <mrow> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>c</mi> <mi>a</mi> <mi>s</mi> </mrow> </msub> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>o</mi> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>
    Wherein, KcasRepresent sleeve pipe thermal conductivity factor；rciRepresent internal surface of sleeve pipe radius；rcoRepresent sleeve outer wall radius.
  5. 5. according to the method for claim 1, it is characterised in that the air and set in oil jacket annular space are calculated using below equation Thermal resistance R between pipe4：
    <mrow> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msub> <mi>h</mi> <mrow> <mi>c</mi> <mn>1</mn> </mrow> </msub> <mo>+</mo> <msub> <mi>h</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> <mo>,</mo> </mrow>
    Wherein, hc1Represent the free convection heat transfer coefficient of air in oil jacket annular space；hr1Represent the heat radiation of air in oil jacket annular space Heat transfer coefficient；rciRepresent internal surface of sleeve pipe radius；
    Heat radiation heat transfer coefficient h is calculated using below equationr1：
    <mrow> <msub> <mi>h</mi> <mrow> <mi>r</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>&amp;delta;F</mi> <mrow> <mi>t</mi> <mi>c</mi> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <msubsup> <mi>T</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> <mrow> <mo>*</mo> <mn>2</mn> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>T</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> <mrow> <mo>*</mo> <mn>2</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msubsup> <mi>T</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> <mo>*</mo> </msubsup> <mo>+</mo> <msubsup> <mi>T</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> <mo>*</mo> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    <mrow> <msubsup> <mi>T</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> <mo>*</mo> </msubsup> <mo>=</mo> <msub> <mi>T</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <mo>+</mo> <mn>273.15</mn> <mo>,</mo> <msubsup> <mi>T</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> <mo>*</mo> </msubsup> <mo>=</mo> <msub> <mi>T</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> <mo>+</mo> <mn>273.15</mn> <mo>,</mo> </mrow>
    <mrow> <mfrac> <mn>1</mn> <msub> <mi>F</mi> <mrow> <mi>t</mi> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>&amp;epsiv;</mi> <mi>o</mi> </msub> </mfrac> <mo>+</mo> <mfrac> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
    Wherein, δ represents Stefan-Boltzmann constants；FtciRepresent that oil-pipe external wall surface is effective to internal surface of sleeve pipe surface emissivity Coefficient；TtoRepresent oil-pipe external wall temperature；TciRepresent internal surface of sleeve pipe temperature；εoRepresent oil-pipe external wall blackness；εciRepresent internal surface of sleeve pipe Blackness；rtoRepresent oil-pipe external wall radius；
    Free convection heat transfer coefficient h is calculated using below equationc1：
    <mrow> <msub> <mi>h</mi> <mrow> <mi>c</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mn>0.049</mn> <msup> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mi>r</mi> </msub> <msub> <mi>P</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <mn>0.33</mn> </msup> <msubsup> <mi>P</mi> <mi>r</mi> <mn>0.074</mn> </msubsup> <msub> <mi>K</mi> <mrow> <mi>h</mi> <mi>a</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <mi>ln</mi> <mfrac> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> </mfrac> </mrow> </mfrac> <mo>,</mo> </mrow>
    <mrow> <msub> <mi>G</mi> <mi>r</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>3</mn> </msup> <msubsup> <mi>g&amp;rho;</mi> <mrow> <mi>a</mi> <mi>n</mi> </mrow> <mn>2</mn> </msubsup> <mi>&amp;beta;</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <msubsup> <mi>U</mi> <mrow> <mi>a</mi> <mi>n</mi> </mrow> <mn>2</mn> </msubsup> </mfrac> <mo>,</mo> </mrow>
    <mrow> <msub> <mi>P</mi> <mi>r</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mrow> <mi>a</mi> <mi>n</mi> </mrow> </msub> <msub> <mi>U</mi> <mrow> <mi>a</mi> <mi>n</mi> </mrow> </msub> </mrow> <msub> <mi>K</mi> <mrow> <mi>h</mi> <mi>a</mi> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>
    Wherein, GrRepresent Grashof numbers；PrRepresent Prandtl numbers；KhaRepresent the thermal conductivity factor of the air of oil jacket annular space；G is represented Acceleration of gravity；ρanRepresent the air of oil jacket annular space in mean temperature TanUnder density；UanRepresent the air of oil jacket annular space flat Equal temperature TanUnder viscosity；CanRepresent the air of oil jacket annular space in mean temperature TanUnder thermal capacitance；β represents that oil jacket annular space is hollow The thermal cubic expansion coefficient of gas.
  6. 6. according to the method for claim 1, it is characterised in that calculate the thermal resistance R between oil pipe inside and outside wall5Including：
    For epimere oil pipe, its thermal resistance R is calculated using below equation5：
    <mrow> <msub> <mi>R</mi> <mn>5</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>t</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>
    Wherein, KtubRepresent oil pipe thermal conductivity factor；rtoRepresent oil-pipe external wall radius；rtiRepresent tube inner wall radius；
    If hypomere oil pipe is instlated tubular, its thermal resistance R is calculated using below equation5：
    <mrow> <msub> <mi>R</mi> <mn>5</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>t</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <msub> <mi>r</mi> <mi>o</mi> </msub> <msub> <mi>r</mi> <mi>i</mi> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>s</mi> </mrow> </msub> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <msub> <mi>r</mi> <mi>i</mi> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> <mn>1</mn> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>t</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>i</mi> <mn>1</mn> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>
    Wherein, KtubRepresent oil pipe thermal conductivity factor；roRepresent instlated tubular outer tube outer wall radius；riRepresent instlated tubular outer tube wall half Footpath；KinsRepresent instlated tubular thermal conductivity factor；rto1Represent instlated tubular outer wall of inner tube radius；rti1Represent instlated tubular inner tube wall radius；
    If hypomere oil pipe is screen casing, its thermal resistance R is calculated using below equation5：
    <mrow> <msub> <mi>R</mi> <mn>5</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>s</mi> <mi>i</mi> <mi>e</mi> </mrow> </msub> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> <mn>2</mn> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>i</mi> <mn>2</mn> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>
    Wherein, KsieRepresent screen casing thermal conductivity factor；rto2Represent screen casing exterior radius；rti2Represent screen casing inwall radius；
    If hypomere oil pipe is plain tubing, using its thermal resistance of below equation R5：
    <mrow> <msub> <mi>R</mi> <mn>5</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;K</mi> <mrow> <mi>t</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>
    Wherein, KtubRepresent oil pipe thermal conductivity factor；rtoRepresent oil-pipe external wall radius；rtiRepresent tube inner wall radius.
  7. 7. according to the method for claim 1, it is characterised in that the thermal convection current heat of oily inner air tube is calculated using below equation Hinder R6：
    <mrow> <msub> <mi>R</mi> <mn>6</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>o</mi> </mrow> </msub> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;h</mi> <mi>f</mi> </msub> <msub> <mi>r</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> <mo>,</mo> </mrow>
    Wherein, hfRepresent the thermal conductivity factor coefficient of oily inner air tube；rtiRepresent tube inner wall radius；rtoRepresent oil-pipe external wall half Footpath.
  8. 8. according to the method for claim 1, it is characterised in that according to R1To R6Pit shaft is calculated radially using below equation On entire thermal resistance：
    R=R1+R2+R3+R4+R5+R6。
  9. 9. according to the method for claim 1, it is characterised in that according to the temperature Ts, the entire thermal resistance and formation temperature, The heat loss of pit shaft diametrically is calculated, including：
    According to law of conservation of energy, the heat loss is calculated using below equation：
    <mrow> <mi>Q</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>e</mi> </msub> </mrow> <mi>R</mi> </mfrac> <mi>d</mi> <mi>l</mi> <mo>,</mo> </mrow>
    Wherein, Q represents pit shaft unit radial heat loss；TeRepresent formation temperature；R represents pit shaft unit radial entire thermal resistance.
  10. 10. according to the method for claim 1, it is characterised in that according to the temperature Ts, the heat loss and the electric point The power of firearm, the air themperature of oil pipe is calculated, including：
    Air themperature in epimere oil pipe is calculated using below equation：
    CmTs- Q/1000+0.4P=CmT 's,
    Air themperature in hypomere oil pipe is calculated using below equation：
    CmTs- Q/1000+0.6P=CmT 's,
    Wherein, T 'sRepresent the air themperature after changing in oil pipe；C represents the specific heat capacity of air；M represents the mass flow of air；P Represent the power of electric igniter；Q represents pit shaft unit radial heat loss.
  11. 11. according to the method for claim 1, it is characterised in that calculate the change of formation temperature in step 7 using below equation Change：
    Te=Tins+ α dl,
    Wherein, TinsRepresent surface temperature；α represents geothermal gradient；TeRepresent formation temperature.
  12. 12. according to the method for claim 1, it is characterised in that the air that calculating is injected into oil pipe adds by electric igniter Temperature T after heats, including：
    Temperature T of the air of epimere oil pipe after electric igniter heats is calculated using below equations：
    CmT+0.4P=CmTs,
    The air of hypomere oil pipe is calculated by electric igniter temperature after heating change T using below equations：
    CmT+0.6P=CmTs,
    Wherein, T represents the initial temperature of air；C represents the specific heat capacity of air；M represents the mass flow of air；P represents electric point The power of firearm.
  13. A kind of 13. determining device of the general gas injection electric ignition well bore temperature distribution of combustion in situ, it is characterised in that including：
    Division unit, for pit shaft to be divided into multiple pit shaft units in the axial direction, the length of each pit shaft unit is dl, makes l =0, k=1, wherein, l represents the length currently calculated, and k represents iterations；
    First computing unit, temperature T of the air of oil pipe after electric igniter heats is injected into for calculatings；
    Second computing unit, for calculating the thermal resistance R on stratum respectively1, cement sheath thermal resistance R2, thermal resistance between sleeve pipe inside and outside wall R3, thermal resistance R between air and sleeve pipe in oil jacket annular space4, thermal resistance R between oil pipe inside and outside wall5And the heat of oily inner air tube Thermal-convection resistance R6；Wherein, the pit shaft radially includes successively from the inside to the outside：Oil pipe, sleeve pipe and cement sheath, it is right in oil jacket annular space Packer should be connected with the position at the top of burning torch, using the packer as boundary, the oil pipe is divided into epimere oil pipe and hypomere Oil pipe, the hypomere oil pipe are instlated tubular, screen casing or plain tubing, are stratum outside the pit shaft；
    3rd computing unit, for according to R1To R6Calculate the entire thermal resistance of pit shaft diametrically；
    4th computing unit, for according to the temperature Ts, the entire thermal resistance and formation temperature, calculate the heat of pit shaft diametrically Loss；
    5th computing unit, for according to the temperature Ts, the heat loss and the electric igniter power, calculate oil pipe Air themperature；
    Unit is iterated to calculate, for making l=l+dl, k=k+1, according to the change of formation temperature, utilizes the second computing unit to the Five computing units are iterated calculating, and until l >=L, then iteration terminates, and obtains the air themperature distribution curve of the oil pipe, its In, L represents the total length of oil pipe.
CN201510864628.0A 2015-12-01 2015-12-01 Method and device for determining temperature distribution of gas injection and electric ignition shaft of in-situ combustion cage Active CN105370255B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510864628.0A CN105370255B (en) 2015-12-01 2015-12-01 Method and device for determining temperature distribution of gas injection and electric ignition shaft of in-situ combustion cage

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510864628.0A CN105370255B (en) 2015-12-01 2015-12-01 Method and device for determining temperature distribution of gas injection and electric ignition shaft of in-situ combustion cage

Publications (2)

Publication Number Publication Date
CN105370255A CN105370255A (en) 2016-03-02
CN105370255B true CN105370255B (en) 2018-01-05

Family

ID=55372813

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510864628.0A Active CN105370255B (en) 2015-12-01 2015-12-01 Method and device for determining temperature distribution of gas injection and electric ignition shaft of in-situ combustion cage

Country Status (1)

Country Link
CN (1) CN105370255B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105672968B (en) * 2016-03-24 2018-02-02 中国石油天然气股份有限公司 Oil reservoir exploitation method
CN105952443A (en) * 2016-05-31 2016-09-21 中国石油天然气股份有限公司 Method and device for determining temperature distribution of shaft of flue gas drive gas injection well
CN106948799B (en) * 2017-02-28 2019-05-07 中国石油天然气股份有限公司 Method and device for determining ignition time of oil reservoir
CN106894801A (en) * 2017-04-10 2017-06-27 中国石油天然气股份有限公司 Method and device for judging ignition of oil layer
CN108825218A (en) * 2018-04-27 2018-11-16 中国石油天然气股份有限公司 formation temperature testing method and device
CN112697313B (en) * 2020-12-08 2023-07-18 长江大学 Heat loss amount evaluation method and device in steam injection shaft heat transfer process

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4369842A (en) * 1981-02-09 1983-01-25 Occidental Oil Shale, Inc. Analyzing oil shale retort off-gas for carbon dioxide to determine the combustion zone temperature
CN103899295A (en) * 2014-03-10 2014-07-02 中国石油天然气股份有限公司 Method and device for detecting heat loss of vertical well shaft in thickened oil steam production and injection
CN104392092B (en) * 2014-10-10 2017-06-13 中国石油天然气股份有限公司 Temperature calculation and control method and device for mixed liquid of gravity fireflood production well
CN104806230B (en) * 2015-02-16 2018-01-05 中国石油天然气股份有限公司 Method for calculating wellbore temperature field and pressure field distribution of supercritical steam injection well

Also Published As

Publication number Publication date
CN105370255A (en) 2016-03-02

Similar Documents

Publication Publication Date Title
CN105370255B (en) Method and device for determining temperature distribution of gas injection and electric ignition shaft of in-situ combustion cage
CN103775058B (en) Method for determining heat loss of shaft
Song et al. Numerical analysis of the heat production performance of a closed loop geothermal system
Serageldin et al. Earth-Air Heat Exchanger thermal performance in Egyptian conditions: Experimental results, mathematical model, and Computational Fluid Dynamics simulation
Bouhacina et al. Numerical investigation of a novel tube design for the geothermal borehole heat exchanger
CN104806230B (en) Method for calculating wellbore temperature field and pressure field distribution of supercritical steam injection well
Ma et al. A coupled heat transfer model of medium-depth downhole coaxial heat exchanger based on the piecewise analytical solution
Cheng et al. An analysis of insulation of abandoned oil wells reused for geothermal power generation
CN104481482B (en) Horizontal well concentric double-pipe gas injection heat insulation analysis method and device
Luo et al. Proposing stratified segmented finite line source (SS-FLS) method for dynamic simulation of medium-deep coaxial borehole heat exchanger in multiple ground layers
CN104850681B (en) Method for calculating temperature, pressure and dryness of injected steam and nitrogen mixture
CN105422084A (en) Method and device for acquiring shaft temperature field of hot water circulation heating viscosity reduction process
Li et al. Heat transfer of a horizontal spiral heat exchanger under groundwater advection
CN104392092B (en) Temperature calculation and control method and device for mixed liquid of gravity fireflood production well
CN104389589A (en) Method and system for determining wellbore temperature field distribution based on hollow rod
Jiashu et al. The influence of thermal boundary conditions of wellbore on the heat extraction performance of deep borehole heat exchangers
CN105370266B (en) In-situ combustion layered electric ignition gas injection well shaft temperature distribution determination method and device
Cravero et al. Numerical simulation of regenerative chambers for glass production plants with a non-equilibrium heat transfer model
CN109033012B (en) Method and device for determining temperature field of hollow sucker rod hot water injection circulation shaft
CN102777159B (en) A kind of note CO 2wellbore of Gas Wells fluidised form is determined and parameter optimization method
CN105184061B (en) The method for numerical simulation of gas-producing well temperature, pressure distribution
CN106948799A (en) Method and device for determining ignition time of oil reservoir
Bravo et al. Design and experimental testing of a 150 kWh thermal battery using thermosiphons embedded in a concrete matrix for power plant flexible operation
Zhou et al. Simulation of thermal transport in aquifer: a GWHP system in Chengdu, China
CN113515863B (en) Method for calculating heat quantity of middle-deep sleeve type heat exchanger based on numerical inversion

Legal Events

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