CN108952652A - Superheated steam is handled up efficiency of heating- utilization evaluation method and system - Google Patents
Superheated steam is handled up efficiency of heating- utilization evaluation method and system Download PDFInfo
- Publication number
- CN108952652A CN108952652A CN201810722673.6A CN201810722673A CN108952652A CN 108952652 A CN108952652 A CN 108952652A CN 201810722673 A CN201810722673 A CN 201810722673A CN 108952652 A CN108952652 A CN 108952652A
- Authority
- CN
- China
- Prior art keywords
- steam
- heat dissipation
- heating
- efficiency
- utilization
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 49
- 230000017525 heat dissipation Effects 0.000 claims abstract description 101
- 238000010793 Steam injection (oil industry) Methods 0.000 claims abstract description 47
- 238000000034 method Methods 0.000 claims abstract description 18
- 230000003068 static effect Effects 0.000 claims abstract description 14
- 239000011435 rock Substances 0.000 claims description 20
- 238000009792 diffusion process Methods 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 229920006395 saturated elastomer Polymers 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 15
- 230000000694 effects Effects 0.000 claims description 9
- 230000008016 vaporization Effects 0.000 claims description 6
- 238000009834 vaporization Methods 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 2
- 239000003921 oil Substances 0.000 description 64
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 3
- 238000013210 evaluation model Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 238000010025 steaming Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 239000010779 crude oil Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/24—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
It handles up efficiency of heating- utilization evaluation method and system the present invention provides a kind of superheated steam, the method includes: obtaining oil reservoir static parameter and steam injection parameter;According to conservation of energy principle, overheated zone heat dissipation range computation model, steam dome heat dissipation range computation model and hydrothermal area's heat dissipation range computation model are established;It is handled up efficiency of heating- utilization model according to overheated zone range computation model, steam dome heat dissipation range computation model and the hydrothermal area range computation model foundation superheated steam that radiates that radiates;It calculates according to the oil reservoir static parameter and steam injection parameter and the superheated steam efficiency of heating- utilization model of handling up and obtains superheated steam and handle up efficiency of heating- utilization.
Description
Technical field
The present invention relates to Oil-Gas Field Development Engineering field, espespecially a kind of superheated steam for considering steam overlap is handled up heat
Utilization efficiency evaluation method and system.
Background technique
Compared with common moist steam, superheated steam institute heat content under uniform pressure is more, specific volume is big, has more in oil reservoir
Big heated perimeter, and more strong physical-chemical reaction occurs for superheated steam and Reservoir Minerals and fluid, so that it is to storage
The improvement result of layer seepage environment is more significant.Therefore, superheated steam injection, which has been increasingly becoming, improves heavy crude reservoir Simulation on whole pay zones effect
One of effective measures of fruit.During superheated steam injection, heat a part of well head injection is lost to pit shaft week through tubing string
Exclosure layer, only partial heat enter oil reservoir, and the heat entered in oil reservoir still some be lost to by pushing up bottom
In country rock, only partial heat is used for heated oil reservoir.Heat effect of the superheated steam injection to oil reservoir in order to make full use of, it is necessary to
Increase the efficiency of heating- utilization during superheated steam injection to the maximum extent.Due to the difference of vapour-liquid density, in superheated steam injection mistake
Steam is often assembled to top of oil horizon in journey, forms steam overlap phenomenon.The oil reservoir heat utilization efficiency of superheated steam injection not only with note
Vapour relating to parameters, is also influenced by steam overlap.Marx-Langenheim method is most common efficiency of heating- utilization
Calculation method, but this method is to be obtained based on common moist steam derivation is infused, and do not account for the influence of steam overlap,
To cause the efficiency of heating- utilization for calculating superheated steam injection using this method not accurate enough.
Summary of the invention
It is theoretical according to steam overlap present invention aims at the reservoir temperature distribution situation during combination superheated steam injection
And conservation of energy principle, building considers that the superheated steam of steam overlap is handled up efficiency of heating- utilization evaluation model, to improve
Hot steam is handled up efficiency of heating- utilization precision of prediction, provides theoretical foundation for superheated steam relevant parameter optimization of handling up.
It handles up efficiency of heating- utilization evaluation side in order to achieve the above object, one embodiment of the invention specifically provides a kind of superheated steam
Method, the method include: obtaining oil reservoir static parameter and steam injection parameter;According to conservation of energy principle, overheated zone heat dissipation model is established
Enclose computation model, steam dome heat dissipation range computation model and hydrothermal area's heat dissipation range computation model;It is radiated according to the overheated zone
Range computation model, steam dome heat dissipation range computation model and hydrothermal area heat dissipation range computation model foundation overheat are steamed
Vapour is handled up efficiency of heating- utilization model;According to the oil reservoir static parameter and steam injection parameter and the superheated steam handle up heat utilization effect
Rate model, which calculates, to be obtained superheated steam and handles up efficiency of heating- utilization.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that radiated range according to the overheated zone
Computation model, steam dome heat dissipation range computation model and hydrothermal area heat dissipation range computation model foundation superheated steam gulp down
Spitting efficiency of heating- utilization model includes: being defined based on efficiency of heating- utilization, according to the overheated zone heat dissipation range computation model, the steaming
Vapour area heat dissipation range computation model and the hydrothermal area range computation model foundation superheated steam that radiates are handled up efficiency of heating- utilization mould
Type.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that according to conservation of energy principle, established
Hot-zone heat dissipation range computation model includes: establishing the thermal energy charge velocity of overheated zone according to instantaneous heat balance principle and to push up bottom total
Heat loss rate and oil reservoir thermal energy are advanced the speed the first relation equation between three;Pass through drawing according to first relation equation
Laplace transform and inverse transformation obtain overheated zone heat dissipation range computation model.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that first relation equation includes:
In above formula, isFor steam charge velocity;hSh, z=DFor the heat enthalpy value of shaft bottom superheated steam;hsFor saturated vapor
Heat enthalpy value;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MRFor oil reservoir thermal capacitance;H is core intersection;
T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original reservoir temperature;D is well depth;
Wherein rbshFor overheated zone bottom heat dissipation radius;TshFor shaft bottom superheat steam temperature;TsFor saturated-steam temperature;k1For overheated zone
Temperature gradient.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that radiate range computation for the overheated zone
Model includes:
In above formula, tDFor non dimensional time, andisFor steam charge velocity;hSh, z=DFor shaft bottom mistake
The heat enthalpy value of hot steam;hsFor the heat enthalpy value of saturated vapor;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion for pushing up bottom
Coefficient;MRFor oil reservoir thermal capacitance;H is core intersection;T is the steam injection time;isFor steam charge velocity;A1Position overheated zone heat dissipation area.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that according to conservation of energy principle, establish and steam
Vapour area heat dissipation range computation model includes: obtaining the second relation equation according to steam overlap theory and instantaneous heat balance principle;Root
Heat dissipation range computation model in steam dome is obtained by Laplace transform and inverse transformation according to second relation equation.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that second relation equation includes:
In above formula: As1For steam dome top layer heat dissipation area;As2For steam dome bottom heat dissipation area;LvFor saturated vapor
The latent heat of vaporization;VsFor steam dome volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the heat for pushing up bottom
Diffusion coefficient;MRFor oil reservoir thermal capacitance;H is core intersection;T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;Ti
For original reservoir temperature;TsFor saturated-steam temperature.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that radiate range computation for the steam dome
Model includes:
In above formula, β is non dimensional time, andisFor steam charge velocity;LvFor saturated vapor
The latent heat of vaporization;VsFor steam dome volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor top bottom
Thermal diffusion coefficient;MRFor oil reservoir thermal capacitance;H is core intersection;tDFor non dimensional time;M is steam overlap coefficient;TiFor original oil
Layer temperature;TsFor saturated-steam temperature;As2For steam dome bottom heat dissipation area.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that according to conservation of energy principle, establish heat
Pool heat dissipation range computation model includes: establishing third relation equation according to instantaneous heat balance principle;According to the third relationship
Equation obtains overheated zone heat dissipation range computation model by Laplace transform and inverse transformation.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that the third relation equation includes:
In above formula: hwsFor the heat enthalpy value of hot water under saturated-steam temperature;hwrFor the heat enthalpy value of hot water under reservoir temperature;Vh
For hydrothermal area's volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MR
For oil reservoir thermal capacitance;T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original reservoir temperature;TsFor saturation
Vapor (steam) temperature;M is steam overlap coefficient;rbsFor steam dome bottom heat dissipation radius;R is Oil reservoir radius;k2For hydrothermal area's temperature ladder
Degree.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that radiate range computation for the overheated zone
Model includes:
In above formula: hwsFor the heat enthalpy value of hot water under saturated-steam temperature;hwrFor the heat enthalpy value of hot water under reservoir temperature;is
For steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MRFor oil reservoir thermal capacitance;tDFor
Non dimensional time;H is core intersection;A2For hydrothermal area's heat dissipation area.
It handles up in efficiency of heating- utilization evaluation method in above-mentioned superheated steam, it is preferred that the hot steam is handled up heat utilization effect
Rate model includes:
In above formula, ηtIt handles up efficiency of heating- utilization for hot steam;QihFor well head thermal energy charge velocity;QibFor shaft bottom oil reservoir heat
It can charge velocity;QoshFor overheated zone oil reservoir thermal energy incrementss;QosFor steam dome oil reservoir thermal energy incrementss;QohFor hydrothermal area's oil reservoir
Thermal energy incrementss;T is gas injection time.
It handles up efficiency of heating- utilization evaluation system the present invention also provides a kind of superheated steam, the evaluation system is obtained comprising data
Modulus block, model generation unit and evaluation unit;The data acquisition module is for obtaining oil reservoir static parameter and steam injection parameter;
The model generation unit is used to establish overheated zone heat dissipation range computation model, steam dome heat dissipation model according to conservation of energy principle
Enclose computation model and hydrothermal area's heat dissipation range computation model;And according to the overheated zone heat dissipation range computation model, the steaming
Vapour area heat dissipation range computation model and the hydrothermal area range computation model foundation superheated steam that radiates are handled up efficiency of heating- utilization mould
Type;The evaluation unit is used for efficiency of heating- utilization of handling up according to the oil reservoir static parameter and steam injection parameter and the superheated steam
Model, which calculates, to be obtained superheated steam and handles up efficiency of heating- utilization.
The beneficial technical effect of the present invention lies in: it is handled up efficiency of heating- utilization evaluation model according to the superheated steam of derivation, energy
The efficiency of heating- utilization enough handled up using analytic method accurate evaluation heavy crude reservoir superheated steam, it is simple and efficient;It is also convenient for analysis viscous crude
Oil reservoir superheated steam is handled up the influence factor of efficiency of heating- utilization, and optimizes related steam injection parameter.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present invention, constitutes part of this application, not
Constitute limitation of the invention.In the accompanying drawings:
Fig. 1 be one embodiment of the invention provided by superheated steam handle up efficiency of heating- utilization evaluation method process signal
Figure;
Fig. 2 is that superheated steam provided by one embodiment of the invention is handled up heating model schematic diagram;
Fig. 3 is efficiency of heating- utilization provided by one embodiment of the invention and steam injection time history;
Fig. 4 is that superheated steam provided by one embodiment of the invention is handled up the structural representation of efficiency of heating- utilization evaluation system
Figure.
Specific embodiment
Understand in order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below with reference to embodiment and attached
Figure, is described in further details the present invention.Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention,
But it is not as a limitation of the invention.
It please refers to shown in Fig. 1, one embodiment of the invention specifically provides a kind of superheated steam and handles up efficiency of heating- utilization evaluation side
Method, the method include: S101 obtains oil reservoir static parameter and steam injection parameter;S102 establishes overheat according to conservation of energy principle
Radiate range computation model, steam dome heat dissipation range computation model and hydrothermal area's heat dissipation range computation model in area;S103 is according to institute
State overheated zone heat dissipation range computation model, steam dome heat dissipation range computation model and hydrothermal area heat dissipation range computation mould
Type establishes superheated steam and handles up efficiency of heating- utilization model;S104 is according to the oil reservoir static parameter and steam injection parameter and the overheat
Steam soak efficiency of heating- utilization model, which calculates, to be obtained superheated steam and handles up efficiency of heating- utilization.In the above-described embodiments, step S103
Described dissipates according to overheated zone heat dissipation range computation model, steam dome heat dissipation range computation model and the hydrothermal area
Heat rating computation model establishes superheated steam efficiency of heating- utilization model of handling up: being defined based on efficiency of heating- utilization, according to institute
State overheated zone heat dissipation range computation model, steam dome heat dissipation range computation model and hydrothermal area heat dissipation range computation mould
Type establishes superheated steam and handles up efficiency of heating- utilization model.
In the above-described embodiments, according to conservation of energy principle, establishing overheated zone heat dissipation range computation model includes: according to wink
When heat balance principle establish the thermal energy charge velocity of overheated zone and top bottom Total heat loss rate and oil reservoir thermal energy advance the speed three
The first relation equation between person;Overheated zone is obtained by Laplace transform and inverse transformation according to first relation equation to dissipate
Heat rating computation model.
First relation equation includes:
In above formula, isFor steam charge velocity;hSh, z=DFor the heat enthalpy value of shaft bottom superheated steam;hsFor saturated vapor
Heat enthalpy value;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MRFor oil reservoir thermal capacitance;H is core intersection;
T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original reservoir temperature;D is well depth;
Wherein rbshFor overheated zone bottom heat dissipation radius;TshFor shaft bottom superheat steam temperature;TsFor saturated-steam temperature;k1For overheated zone
Temperature gradient.
The overheated zone heat dissipation range computation model includes:
In above formula, tDFor non dimensional time, andisFor steam charge velocity;hSh, z=DFor shaft bottom mistake
The heat enthalpy value of hot steam;hsFor the heat enthalpy value of saturated vapor;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion for pushing up bottom
Coefficient;MRFor oil reservoir thermal capacitance;H is core intersection;T is the steam injection time;isFor steam charge velocity;A1Position overheated zone heat dissipation area.
In the above-described embodiments, according to conservation of energy principle, establishing steam dome heat dissipation range computation model includes: according to steaming
Vapour onlap theory and instantaneous heat balance principle obtain the second relation equation;Become according to second relation equation by Laplce
It changes and inverse transformation obtains steam dome heat dissipation range computation model.
Second relation equation includes:
In above formula: As1For steam dome top layer heat dissipation area;As2For steam dome bottom heat dissipation area;LvFor saturated vapor
The latent heat of vaporization;VsFor steam dome volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the heat for pushing up bottom
Diffusion coefficient;MRFor oil reservoir thermal capacitance;H is core intersection;T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;Ti
For original reservoir temperature;TsFor saturated-steam temperature.
The steam dome heat dissipation range computation model includes:
In above formula, β is non dimensional time, andAs1For steam dome top layer heat dissipation area;As2For
Steam dome bottom heat dissipation area;LvFor the latent heat of vaporization of saturated vapor;VsFor steam dome volume;isFor steam charge velocity;λsFor
Push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MRFor oil reservoir thermal capacitance;H is core intersection;When t is steam injection
Between;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original reservoir temperature;TsFor saturated-steam temperature;TsIt is steamed for saturation
Stripping temperature;As2For steam dome bottom heat dissipation area.
In the above-described embodiments, according to conservation of energy principle, establishing hydrothermal area's heat dissipation range computation model includes: according to wink
When heat balance principle establish third relation equation;It is obtained according to the third relation equation by Laplace transform and inverse transformation
Radiate range computation model for overheated zone.
The third relation equation includes:
In above formula: hwsFor the heat enthalpy value of hot water under saturated-steam temperature;hwrFor the heat enthalpy value of hot water under reservoir temperature;Vh
For hydrothermal area's volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MR
For oil reservoir thermal capacitance;T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original reservoir temperature;TsFor saturation
Vapor (steam) temperature;M is steam overlap coefficient;rbsFor steam dome bottom heat dissipation radius;R is Oil reservoir radius;k2For hydrothermal area's temperature ladder
Degree.
The overheated zone heat dissipation range computation model includes:
In above formula: hwsFor the heat enthalpy value of hot water under saturated-steam temperature;hwrFor the heat enthalpy value of hot water under reservoir temperature;Vh
For hydrothermal area's volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MR
For oil reservoir thermal capacitance;T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original reservoir temperature;TsFor saturation
Vapor (steam) temperature;tDFor non dimensional time;H is core intersection;A2For hydrothermal area's heat dissipation area.
Above-mentioned three's computing module is combined with this, available hot steam is handled up efficiency of heating- utilization model, wherein the heat is steamed
Vapour efficiency of heating- utilization model of handling up includes:
In above formula, ηtIt handles up efficiency of heating- utilization for hot steam;QihFor well head thermal energy charge velocity;QibFor shaft bottom oil reservoir heat
It can charge velocity;QoshFor overheated zone oil reservoir thermal energy incrementss;QosFor steam dome oil reservoir thermal energy incrementss;QohFor hydrothermal area's oil reservoir
Thermal energy incrementss;T is gas injection time.
It handles up efficiency of heating- utilization evaluation method for clearer explanation superheated steam provided by the present invention, below with entirety
Process is further explained explanation to above-mentioned each link, and relevant technical staff in the field is when it is found that following explanation is only to help to manage
It solves superheated steam provided by the present invention to handle up efficiency of heating- utilization evaluation method, any restriction is not constituted to it.
Please refer to shown in Fig. 2, basic assumption condition includes: (1) oil reservoir homogeneous, uniform thickness, and oil reservoir thermal physical property parameter not with
Temperature change.(2) Wellhead steam injection speed and steam superheat are invariable.(3) overheated zone is the cylinder centered on steam injection well,
The tilted interface of skew lines is formed between steam dome and hydrothermal area, steam dome is the rotary table centered on steam injection well.(4) overheated zone
For temperature by shaft bottom superheat steam temperature linear decrease to saturated-steam temperature, constant steam dome temperature is saturated-steam temperature, heat
Pool temperature is by saturated-steam temperature linear decrease to original reservoir temperature.(5) on hydrothermal area's difference vertical position radially
Temperature gradient it is equal, i.e. steam dome leading edge and hydrothermal area's leading edge is parallel to each other.
Under above-mentioned assumed condition, overheated zone heat dissipation area computation model can according to instantaneous heat balance principle, overheated zone
Thermal energy charge velocity should be equal to area top bottom Total heat loss rate and oil reservoir thermal energy the sum of be advanced the speed, it may be assumed that
In formula: isFor steam charge velocity, kg/s;hSh, z=DFor the heat enthalpy value of shaft bottom superheated steam, J/kg;hsFor saturation
The heat enthalpy value of steam, J/kg;λsTo push up seat rock thermal coefficient, W/ (m DEG C);αsFor the thermal diffusion coefficient for pushing up bottom, m2/
h;MRFor oil reservoir thermal capacitance, J/ (m3·℃);H is core intersection, m;T is steam injection time, h;δ reaches certain point institute in oil reservoir for thermal energy
It takes time, h;TiFor original reservoir temperature, DEG C;D is well depth, m;Wherein rbshFor overheated zone bottom radiate radius,
m;TshFor shaft bottom superheat steam temperature, DEG C;TsFor saturated-steam temperature, DEG C.
Enable dA1=2 π r (k1r+Tsh-Ti) dr, Laplace transformation and inverse transformation are used to above formula, can be obtained:
In formula: tDFor non dimensional time, and
The steam dome heat dissipation area computation model is as follows: in steam dome, due to saturated vapor and condensed liquid density variation
It is very big, cause steam to be assembled to top of oil horizon, forms skewed displacing front.It is theoretical according to Van Lookeen steam overlap,
Steam dome leading edge equation are as follows:
In formula: h and hstRespectively core intersection and steam thickness, m;resAnd rbsRespectively steam dome top layer heat dissipation radius
With steam dome bottom heat dissipation radius, m;M* is pseudo mobility ratio, zero dimension;ARDFor form factor, zero dimension.
Defining steam overlap Coefficient m is that the ratio between bottom heat dissipation radius, i.e. m=r are pushed up in steam domees/rbs(m≥1).Work as r=rbs
When, hst=h, then have:
Instantaneous heat balance principle is utilized to steam dome, can be obtained:
In formula: As1For steam dome top layer heat dissipation area, m2;As2For steam dome bottom heat dissipation area, m2;LvFor saturated vapor
The latent heat of vaporization, J/kg;VsFor steam dome volume, m3。
Because of dAs1=m2dAs2, and as it is assumed that steam dome is the rotary table centered on steam injection well, so havingIt is substituted into formula (12), and with Laplace transformation and inverse transformation, can be obtained:
In formula: β is non dimensional time, and
Hydrothermal area's heat dissipation area computation model is as follows: instantaneous heat balance principle utilized to hydrothermal area, can be obtained,
In formula: hwsFor the heat enthalpy value of hot water under saturated-steam temperature, J/kg;hwrFor the heat enthalpy value of hot water under reservoir temperature,
J/kg;VhFor hydrothermal area's volume, m3。
Because of dVs=π h (2r+mrbs-rbs) dr, and enable dA2=π (2r+mrbs-rbs)[k2(r-rbs)+Ts-Ti] dr, by its generation
Enter formula (22), and with Laplace transformation and inverse transformation, can obtain:
The efficiency of heating- utilization computation model is as follows:
Define pit shaft heat utilization efficiency η1Heat, which is carried, for shaft bottom steam accounts for the percentage that well head steam carries heat, it may be assumed that
In formula: QihFor well head thermal energy charge velocity, J/s;QihFor shaft bottom oil reservoir thermal energy charge velocity, J/s;hSh, z=jFor well
Superheated steam heat enthalpy value (j=0) at a temperature of mouthful or the superheated steam heat enthalpy value (j=D) under bottom hole temperature (BHT), J/kg;D is well
It is deep, m.
Define oil reservoir heat utilization efficiency η2The percentage that oil reservoir injects heat is accounted for for heating zone oil reservoir thermal energy incrementss, it may be assumed that
In formula: QoshFor overheated zone oil reservoir thermal energy incrementss, J/s;QosFor steam dome oil reservoir thermal energy incrementss;QohFor hot water
Area's oil reservoir thermal energy incrementss.
Define total heat utilization efficiency ηtThe percentage that well head steam carries heat is accounted for for heating zone oil reservoir thermal energy incrementss, it may be assumed that
With this, by above-mentioned each model can the efficiency of heating- utilization handled up of accurate evaluation heavy crude reservoir superheated steam, it is easy
Fast, heavy crude reservoir superheated steam can also be analyzed to handle up the influence factor of efficiency of heating- utilization, and optimize related steam injection parameter.
It please refers to shown in Fig. 4, handles up efficiency of heating- utilization evaluation system the present invention also provides a kind of superheated steam, the evaluation
System includes data acquisition module, model generation unit and evaluation unit;The data acquisition module is for obtaining oil reservoir static state
Parameter and steam injection parameter;The model generation unit is used to establish overheated zone heat dissipation range computation mould according to conservation of energy principle
Type, steam dome heat dissipation range computation model and hydrothermal area's heat dissipation range computation model;And according to overheated zone heat dissipation range
Computation model, steam dome heat dissipation range computation model and hydrothermal area heat dissipation range computation model foundation superheated steam gulp down
Spit efficiency of heating- utilization model;The evaluation unit is used for according to the oil reservoir static parameter and steam injection parameter and the superheated steam
Efficiency of heating- utilization of handling up model, which calculates, to be obtained superheated steam and handles up efficiency of heating- utilization.
Handle up the effect of efficiency of heating- utilization evaluation method and system to verify superheated steam provided by the present invention, below with
Specific example carries out validation verification;
The oil reservoir heat utilization efficiency of its certain a bite superheated steam injection well is carried out based on Kazakhstan's oil deposit parameter
Example calculation.Casing programme parameter and reservoir physical property parameter are as shown in 1 casing programme parameter of table and reservoir physical property parameter, steam injection
Speed is 6t/h, and shaft bottom vapor (steam) temperature is 275 DEG C, and shaft bottom steam superheat is 30 DEG C.
Table 1
It is as shown in Figure 3 that efficiency of heating- utilization changes over time relationship.As can be seen from Figure 3: new model calculated result ratio Marx-
It is less than normal that Langenheim method calculates checkout result, but has preferable consistency with CMG STARS reservoir numerical simulation result,
Maximum relative error is 4.8%, shows that new model has higher computational accuracy compared with Marx-Langenheim method.
The beneficial technical effect of the present invention lies in: it is handled up efficiency of heating- utilization evaluation model according to the superheated steam of derivation, energy
The efficiency of heating- utilization enough handled up using analytic method accurate evaluation heavy crude reservoir superheated steam, it is simple and efficient;It is also convenient for analysis viscous crude
Oil reservoir superheated steam is handled up the influence factor of efficiency of heating- utilization, and optimizes related steam injection parameter.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the present invention, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection scope of invention.
Claims (13)
- The efficiency of heating- utilization evaluation method 1. a kind of superheated steam is handled up, which is characterized in that the method includes:Obtain oil reservoir static parameter and steam injection parameter;According to conservation of energy principle, overheated zone heat dissipation range computation model, steam dome heat dissipation range computation model and hot water are established Radiate range computation model in area;According to overheated zone heat dissipation range computation model, steam dome heat dissipation range computation model and hydrothermal area heat dissipation Range computation model foundation superheated steam is handled up efficiency of heating- utilization model;According to the oil reservoir static parameter and steam injection parameter and the superheated steam handle up efficiency of heating- utilization model calculating obtained Hot steam is handled up efficiency of heating- utilization.
- The efficiency of heating- utilization evaluation method 2. superheated steam according to claim 1 is handled up, which is characterized in that according to the mistake Hot-zone heat dissipation range computation model, steam dome heat dissipation range computation model and hydrothermal area heat dissipation range computation model are built Vertical superheated steam efficiency of heating- utilization model of handling up includes: being defined based on efficiency of heating- utilization, is radiated range meter according to the overheated zone Model, steam dome heat dissipation range computation model and hydrothermal area heat dissipation range computation model foundation superheated steam is calculated to handle up Efficiency of heating- utilization model.
- The efficiency of heating- utilization evaluation method 3. superheated steam according to claim 1 is handled up, which is characterized in that kept according to energy Permanent principle, establishing overheated zone heat dissipation range computation model includes: being injected according to the thermal energy that instantaneous heat balance principle establishes overheated zone Rate and top bottom Total heat loss rate and oil reservoir thermal energy are advanced the speed the first relation equation between three;According to described first Relation equation obtains overheated zone heat dissipation range computation model by Laplace transform and inverse transformation.
- The efficiency of heating- utilization evaluation method 4. superheated steam according to claim 3 is handled up, which is characterized in that described first closes It is that equation includes:In above formula, isFor steam charge velocity;hSh, z=DFor the heat enthalpy value of shaft bottom superheated steam;hsFor the heat content of saturated vapor Value;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MRFor oil reservoir thermal capacitance;H is core intersection;T is The steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original reservoir temperature;D is well depth;Its Middle rbshFor overheated zone bottom heat dissipation radius;TshFor shaft bottom superheat steam temperature;TsFor saturated-steam temperature;k1For overheated zone temperature Spend gradient.
- The efficiency of heating- utilization evaluation method 5. superheated steam according to claim 4 is handled up, which is characterized in that the overheated zone Heat dissipation range computation model includes:In above formula, tDFor non dimensional time, andisFor steam charge velocity;hSh, z=DIt overheats and steams for shaft bottom The heat enthalpy value of vapour;hsFor the heat enthalpy value of saturated vapor;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom; MRFor oil reservoir thermal capacitance;H is core intersection;T is the steam injection time;isFor steam charge velocity;A1Position overheated zone heat dissipation area.
- The efficiency of heating- utilization evaluation method 6. superheated steam according to claim 1 is handled up, which is characterized in that kept according to energy Permanent principle, establishing steam dome heat dissipation range computation model includes: obtaining the according to steam overlap theory and instantaneous heat balance principle Two relation equations;Heat dissipation range computation in steam dome is obtained by Laplace transform and inverse transformation according to second relation equation Model.
- The efficiency of heating- utilization evaluation method 7. superheated steam according to claim 6 is handled up, which is characterized in that described second closes It is that equation includes:In above formula: As1For steam dome top layer heat dissipation area;As2For steam dome bottom heat dissipation area;LvFor the vaporization of saturated vapor Latent heat;VsFor steam dome volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion for pushing up bottom Coefficient;MRFor oil reservoir thermal capacitance;H is core intersection;T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original Beginning reservoir temperature;TsFor saturated-steam temperature.
- The efficiency of heating- utilization evaluation method 8. superheated steam according to claim 7 is handled up, which is characterized in that the steam dome Heat dissipation range computation model includes:In above formula, β is non dimensional time, andisFor steam charge velocity;LvFor the vapour of saturated vapor Change latent heat;VsFor steam dome volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the thermal expansion for pushing up bottom Dissipate coefficient;MRFor oil reservoir thermal capacitance;H is core intersection;tDFor non dimensional time;M is steam overlap coefficient;TiFor original reservoir temperature Degree;TsFor saturated-steam temperature;As2For steam dome bottom heat dissipation area.
- The efficiency of heating- utilization evaluation method 9. superheated steam according to claim 1 is handled up, which is characterized in that kept according to energy Permanent principle, establishing hydrothermal area's heat dissipation range computation model includes: establishing third relation equation according to instantaneous heat balance principle;According to The third relation equation obtains overheated zone heat dissipation range computation model by Laplace transform and inverse transformation.
- The efficiency of heating- utilization evaluation method 10. superheated steam according to claim 9 is handled up, which is characterized in that the third Relation equation includes:In above formula: hwsFor the heat enthalpy value of hot water under saturated-steam temperature;hwrFor the heat enthalpy value of hot water under reservoir temperature;VhFor heat Pool volume;isFor steam charge velocity;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MRFor oil Layer thermal capacitance;T is the steam injection time;δ is that thermal energy reaches in oil reservoir the time required to certain point;TiFor original reservoir temperature;TsFor saturated vapor Temperature;M is steam overlap coefficient;rbsFor steam dome bottom heat dissipation radius;R is Oil reservoir radius;k2For hydrothermal area's temperature gradient.
- The efficiency of heating- utilization evaluation method 11. superheated steam according to claim 10 is handled up, which is characterized in that the overheat Area's heat dissipation range computation model includes:In above formula: hwsFor the heat enthalpy value of hot water under saturated-steam temperature;hwrFor the heat enthalpy value of hot water under reservoir temperature;isTo steam Vapour charge velocity;λsTo push up seat rock thermal coefficient;αsFor the thermal diffusion coefficient for pushing up bottom;MRFor oil reservoir thermal capacitance;tDFor it is no because The secondary time;H is core intersection;A2For hydrothermal area's heat dissipation area.
- The efficiency of heating- utilization evaluation method 12. superheated steam according to claim 1 is handled up, which is characterized in that the heat is steamed Vapour efficiency of heating- utilization model of handling up includes:In above formula, ηtIt handles up efficiency of heating- utilization for hot steam;QihFor well head thermal energy charge velocity;QibFor shaft bottom oil reservoir thermal energy note Enter rate;QoshFor overheated zone oil reservoir thermal energy incrementss;QosFor steam dome oil reservoir thermal energy incrementss;QohFor hydrothermal area's oil reservoir thermal energy Incrementss;T is gas injection time.
- The efficiency of heating- utilization evaluation system 13. a kind of superheated steam is handled up, which is characterized in that the evaluation system includes data acquisition Module, model generation unit and evaluation unit;The data acquisition module is for obtaining oil reservoir static parameter and steam injection parameter;The model generation unit is used to establish overheated zone heat dissipation range computation model, steam dome dissipates according to conservation of energy principle Heat rating computation model and hydrothermal area's heat dissipation range computation model;And according to the overheated zone heat dissipation range computation model, institute State steam dome heat dissipation range computation model and the hydrothermal area radiate range computation model foundation superheated steam handle up heat utilization effect Rate model;The evaluation unit is used for heat utilization effect of handling up according to the oil reservoir static parameter and steam injection parameter and the superheated steam Rate model, which calculates, to be obtained superheated steam and handles up efficiency of heating- utilization.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810722673.6A CN108952652A (en) | 2018-07-04 | 2018-07-04 | Superheated steam is handled up efficiency of heating- utilization evaluation method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810722673.6A CN108952652A (en) | 2018-07-04 | 2018-07-04 | Superheated steam is handled up efficiency of heating- utilization evaluation method and system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108952652A true CN108952652A (en) | 2018-12-07 |
Family
ID=64485524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810722673.6A Pending CN108952652A (en) | 2018-07-04 | 2018-07-04 | Superheated steam is handled up efficiency of heating- utilization evaluation method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108952652A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115341885A (en) * | 2022-08-29 | 2022-11-15 | 中国石油大学(北京) | Method and device for determining volume of steam-driven channeling channel |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130186627A1 (en) * | 2011-06-13 | 2013-07-25 | Kimberly Jantunen Cross | Additives for improving hydrocarbon recovery |
CN103775058A (en) * | 2013-12-31 | 2014-05-07 | 中国石油天然气股份有限公司 | Shaft heat loss determining method |
CN103899295A (en) * | 2014-03-10 | 2014-07-02 | 中国石油天然气股份有限公司 | Method and device for detecting heat losses of plumb shaft barrel in thickened oil steam production and injection process |
CN104850681A (en) * | 2015-04-16 | 2015-08-19 | 中国石油天然气股份有限公司 | Method for calculating temperature, pressure and dryness of injected steam and nitrogen mixture |
CN106761733A (en) * | 2017-01-06 | 2017-05-31 | 中国海洋石油总公司 | A kind of horizontal wells in heavy oil reservoir steam soak initial productivity Forecasting Methodology |
-
2018
- 2018-07-04 CN CN201810722673.6A patent/CN108952652A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130186627A1 (en) * | 2011-06-13 | 2013-07-25 | Kimberly Jantunen Cross | Additives for improving hydrocarbon recovery |
CN103775058A (en) * | 2013-12-31 | 2014-05-07 | 中国石油天然气股份有限公司 | Shaft heat loss determining method |
CN103899295A (en) * | 2014-03-10 | 2014-07-02 | 中国石油天然气股份有限公司 | Method and device for detecting heat losses of plumb shaft barrel in thickened oil steam production and injection process |
CN104850681A (en) * | 2015-04-16 | 2015-08-19 | 中国石油天然气股份有限公司 | Method for calculating temperature, pressure and dryness of injected steam and nitrogen mixture |
CN106761733A (en) * | 2017-01-06 | 2017-05-31 | 中国海洋石油总公司 | A kind of horizontal wells in heavy oil reservoir steam soak initial productivity Forecasting Methodology |
Non-Patent Citations (2)
Title |
---|
何聪鸽 等: ""稠油油藏注过热蒸汽加热半径计算新模型"", 《大庆石油地质与开发》 * |
赖令彬 等: ""基于蒸汽超覆的稠油蒸汽驱地层热效率计算模型"", 《重庆大学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115341885A (en) * | 2022-08-29 | 2022-11-15 | 中国石油大学(北京) | Method and device for determining volume of steam-driven channeling channel |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104806230B (en) | The Wellbore Temperature Field of overcritical steam injection well and the computational methods of pressure field distribution | |
Han et al. | Fast thermal simulation of a heated crude oil pipeline with a BFC-Based POD reduced-order model | |
Croucher et al. | Waiwera: A parallel open-source geothermal flow simulator | |
Gui et al. | Wave impact simulations by an improved ISPH model | |
BR102013025346A2 (en) | METHOD, SYSTEM, AND ONE OR MORE COMPUTER-LEGIBLE MEANS | |
González‐Albuixech et al. | Convergence of domain integrals for stress intensity factor extraction in 2‐D curved cracks problems with the extended finite element method | |
CN104091036A (en) | Heat conduction modeling and calculating method for natural circulation vapor generator | |
Fambri et al. | An efficient semi‐implicit method for three‐dimensional non‐hydrostatic flows in compliant arterial vessels | |
Yu et al. | Efficient finite-volume model for shallow-water flows using an implicit dual time-stepping method | |
CN108952652A (en) | Superheated steam is handled up efficiency of heating- utilization evaluation method and system | |
Lai et al. | Discontinuous Galerkin method for 1D shallow water flow in nonrectangular and nonprismatic channels | |
Zhang et al. | A high-order well-balanced positivity-preserving moving mesh DG method for the shallow water equations with non-flat bottom topography | |
Panou et al. | Direct and indirect estimation of the variance–covariance matrix of the parameters of a fitted ellipse and a triaxial ellipsoid | |
Zhu et al. | Thermal streamline simulation: steam floods | |
Xu et al. | Numerical analysis on flow instability of parallel channels in steam generator for sodium‐cooled fast reactor | |
Zhang et al. | One-dimensional coupled model of surface water flow and solute transport for basin fertigation | |
Behzadi et al. | Modeling of boiling flow in a channel under static along rolling and heaving motions | |
CN106321039B (en) | Method and apparatus for obtaining shaft bottom steam parameter in steam injection recovery | |
Wang et al. | A component-based Eulerian--Lagrangian formulation for multicomponent multiphase compositional flow and transport in porous media | |
Choi et al. | A horizontally curvilinear non‐hydrostatic model for simulating nonlinear wave motion in curved boundaries | |
Fan et al. | A fully‐implicit numerical algorithm of two‐fluid two‐phase flow model using Jacobian‐free Newton–Krylov method | |
Li et al. | A Complex Variable Solution for Rectangle Pipe Jacking in Elastic Half‐Plane | |
Wang et al. | A highly accurate regular domain collocation method for solving potential problems in the irregular doubly connected domains | |
Xu et al. | A 3D staggered Lagrangian scheme for ideal magnetohydrodynamics on unstructured meshes | |
Demazière et al. | Development of three-dimensional capabilities for modelling stationary fluctuations in nuclear reactor cores |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181207 |
|
RJ01 | Rejection of invention patent application after publication |