CN104989351A - Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process - Google Patents
Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process Download PDFInfo
- Publication number
- CN104989351A CN104989351A CN201510293979.0A CN201510293979A CN104989351A CN 104989351 A CN104989351 A CN 104989351A CN 201510293979 A CN201510293979 A CN 201510293979A CN 104989351 A CN104989351 A CN 104989351A
- Authority
- CN
- China
- Prior art keywords
- partiald
- rho
- temperature
- lambda
- equation
- 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
Landscapes
- Complex Calculations (AREA)
Abstract
The invention belongs to the technical field of oil-gas reservoir development engineering management, and relates to a dryness, temperature and pressure coupling predicting method in an oil-gas well gas injection process, in particular to the coupling analysis of dryness and pressure of oil-gas two-phase instantaneous flow and shaft-stratum unsteady heat conduction, the mathematical modeling and the numerical simulation so as to predict dryness, temperature and pressure distribution situations of an ultra-deep well. The coupling predicting method is used for solving the problem that in a traditional scheme, dryness, temperature and pressure are not accurately predicted. The method includes the implementation steps of A, establishing a stratum heat conduction model; B, establishing a differential equation coupling model; C, solving the differential equation coupling model. The method is suitable for oil-gas reservoir development.
Description
Technical field
The invention belongs to development of oil and gas reservoir administrative skill field, it relates to mass dryness fraction in Oil/gas Well gas injection process, temperature and coupling pressure Forecasting Methodology, be specifically related in HTHP Oil/gas Well gas injection process, the mass dryness fraction of oil gas two-phase transient state stream, pressure and pit shaft-stratum Geometry symmetry mode coupling analysis, mathematical modeling and numerical simulation, prediction ultradeep well mass dryness fraction, temperature, pressure distribution situation.
Background technology
In the gas injection operating mode of Oil-Gas Well Engineering, conventional method considers micro unit body conservation of mechanical energy and change of internal energy usually, but often ignore the interchange of heat with external environment condition, think that the change of internal steam work equals the change of outside, and also have ignored the friction loss of the steam of gas injection process, then calculate wellbore pressure and temperature, thus cause predicting extremely inaccurate to mass dryness fraction, temperature, pressure distribution in pit shaft.
Usually have ignored when building energy-balance equation the energy loss that frictional resistance causes and think that formation temperature is constant, just be subject to the impact of the degree of depth and geothermal gradient, the formation temperature that the heat trnasfer that also have ignored internal steam work change and the surrounding environment produced by frictional resistance causes can because of the change of conducting heat between pit shaft and stratum.The unsteady-state heat transfer on the friction loss in gas injection process and pit shaft-stratum is the key factor building differential equation coupling system model, gas-bearing formation infringement or improve the important parameter of degree, significant for mass dryness fraction in accurately predicting pit shaft, temperature, pressure distribution.
Oil/gas Well thermodynamic analysis: before Modling model, the hypothesis about steam flow and heat transfer is as follows:
(1) fluid physical property and formed by the impact of the degree of depth and temperature.
(2) all parameters (speed, pressure, temperature and aridity) of steam are constants.
(3) be stable from pipeline to the heat transfer of the second contact surface, from the second contact surface to the heat transfer on stratum be non-stable.
(4) vapor stream is regarded as one dimension two-phase homogeneous fluid.
(5) physical property (steam and some heat-barrier materials except) of other materials is not by the impact of time and temperature.
(6) there is a kind of linear relationship between the thermal conductivity of heat-barrier material and temperature.
(7) geothermal gradient is constant.
The hypothesis of formation temperature field is as follows:
(1) temperature in stratum around pit shaft distributes axisymmetricly.
(2) when injecting steam toward pit shaft, the formation temperature away from pit shaft axle center is identical.
(3) formation temperature field does not have inner-growth model and is transient state conduction.
Barometric gradient: because steam injection is constant concrete mass flow, that is: flow into micro unit quality and equal to flow out micro unit quality (see Fig. 1), it follows mass-conservation equation:
M=ρ
1υ
1A=ρ
2υ
2A=ρ
mυ
mA (1)
It is followed
in a period of time dt, the momentum that micro-body is subject to is Fdt=P
1adt+ ρ
mag cos θ dzdt-P
2adt-τ
fdt, in dt during this period of time, the change of momentum is ρ
2υ
2adt υ
2-ρ
1υ
1adt υ
1=Δ (m υ).It follows momentum theorem:
Gas mass dryness fraction: in view of the energy loss caused by vapor stream frictional resistance, we draw following energy conservation equation:
In this equation, dW represents that, because of the energy loss that frictional resistance causes between steam and tube wall, Hm represents the enthalpy of mixed flow, and it can be defined as:
H
m=H
sX+H
w(1-X) (4)
In this equation, Hs represents the enthalpy of saturated vapour, and Hw represents saturation water, can be derived by equation (4):
Because
And
So can derive:
In addition, the rate equation of mixed flow is:
So then have,
In this equation,
Because the direction of steam flow and frictional force is contrary, so the energy loss dW caused because of frictional resistance between steam and tube wall is negative value.In unit interval, the frictional force of dz is:
Equation (1) gives the definition of quality, thus we obtain the aridity equation of following gas model:
In this equation,
Wellbore heat: note, assuming that from pipeline to the heat transfer of the second contact surface be stable (see Fig. 2).It is followed
Second contact surface to the radial direction heat transfer of surrounding formation is:
In conjunction with equation (9) and (10), thus draw the heat transfer model amount between vapor stream and surrounding formation, as follows:
If α=π is D
tou
tok
e/ (0.5D
tou
tof (t
d)+K
e), then
Summary of the invention
Technical problem to be solved by this invention is: propose mass dryness fraction, temperature and coupling pressure Forecasting Methodology in a kind of Oil/gas Well gas injection process, solves in traditional scheme the inaccurate problem of mass dryness fraction, temperature and pressure prediction.
The technical solution adopted for the present invention to solve the technical problems is: mass dryness fraction, temperature and coupling pressure Forecasting Methodology in Oil/gas Well gas injection process, comprising:
A. stratum heat transfer model is built;
B. differential equation coupling model is built;
C. differential equation coupling model is solved.
Further, based on from the second contact surface to the hypothesis of stratum unsteady heat transfer in steps A, build stratum heat transfer model, specifically comprise:
Primary condition:
T
e=T
α+ γ
zif, t=0.
Fringe conditions:
Nondimensional number
with
substitute into equation (12),
Fringe conditions is converted into:
Relation between saturated saturated-steam temperature and pressure:
Further, in step B, described differential equation coupling model considers friction loss and pit shaft-stratum unsteady-state heat transfer, and the differential equation coupling model of structure is:
Further, in step C, before solving model, the parameter in model is handled as follows:
1) vapor (steam) velocity V
m:
2) vapour density Pm: because water vapour air-flow is biphase gas and liquid flow, therefore applies the averag density that Beggs-Brill method calculates mixture;
3) frictional force T
f:
4) the friction factor f of gas-liquid mixture
m, f
mthe function about Reynolds number and absolute roughness e,
5) thermal transmittance U
to:
Wherein, λ
insand λ
cembe respectively the thermal conductivity of thermal insulation material and cement sheath.Hc and hr is respectively convective heat-transfer coefficient and radiation heat transfer coefficient;
6) nondimensional time function:
Further, in step C, adopt interactive quadravalence Runge-Kutta finite difference method to solve differential equation coupling model, comprising:
Step 1. provides θ respectively
0, T
0, P
0, X
0, Q
0and T
e0initial value;
Step 2. calculates all coefficients in coupling model;
Step 3. makes the differential method be converted into function f
i(i=1,2,3,4), can obtain following coupling function system:
The equation group that step 4. solves above obtains P
k, X
k, Q
k, T
k, and T
ek, can coefficient be drawn:
a
i=f
i(P
k,X
k,Q
k,T
k),j=1,2,3,4;
Step 5. makes
te is obtained by separating following equation,
As time j, radial i when degree of depth z, T
j irepresent temperature, i=1; 2 ... M, j=1,2 ... N; Wherein M and N represents time and last radial node respectively; Finite difference method is adopted to carry out discrete obtaining to equation (23):
In this equation, τ represents the time interval, and ζ represents spaced radial; It can convert canonical form to, as follows:
Then calculus of finite differences discrete boundary condition is used; If rD=1, then:
It meets
If rD=N, then:
In conjunction with equation (25), (27) and (28), the numerical solution of formation temperature Te can be calculated;
Step 6., as rD=1, substitutes into equation system (23) Te, can obtain:
Step 7. calculates
With
d
j=f
j(P
k+hc
1,X
k+hc
2,Q
k+hc
3,T
k+hc
4)(j=1,2,3,4)
Step 8. calculates gas-liquid mixture in (K+1) this point, and mass dryness fraction, the pressure and temperature of gas can obtain:
Step 9. is T=T
k+1substitute into the fringe conditions of equation (23), can T be obtained by finite difference method
e,
k+1;
Step 10. repeats step 2-9, until calculate P
nx
nq
n, T
nand T
en.
The invention has the beneficial effects as follows: when setting up coupling model, considering formation temperature that the internal steam work change that produced by frictional resistance and the heat trnasfer of surrounding environment cause can because of the change of conducting heat between pit shaft and stratum, thus accurately predicting ultradeep well mass dryness fraction, temperature, pressure distribution situation.
Accompanying drawing explanation
Fig. 1 is the force analysis figure of micro unit;
Fig. 2 is shaft structure schematic diagram;
Fig. 3 is Forecasting Methodology flow chart of the present invention;
Fig. 4 is the change curve of pressure with injection length of steam;
Fig. 5 is the change curve of mass dryness fraction with injection length of steam;
Fig. 6 is the change curve of the temperature in well with injection length;
Fig. 7 is the change curve of near-bottom temperature with injection length.
Detailed description of the invention
The present invention is intended to propose mass dryness fraction, temperature and coupling pressure Forecasting Methodology in a kind of Oil/gas Well gas injection process, solves in traditional scheme the inaccurate problem of mass dryness fraction, temperature and pressure prediction.
As shown in Figure 3, this Forecasting Methodology comprises:
A. stratum heat transfer model is built;
B. differential equation coupling model is built;
C. differential equation coupling model is solved.
Be described in detail for each step from concrete enforcement below:
Conduct heat in stratum: according to from the second contact surface to the hypothesis of stratum unsteady heat transfer, and we draw with sub-surface heat transfer model:
Primary condition:
T
e=T
α+ γ
zif, t=0.
Fringe conditions:
Nondimensional number
with
substitute into equation (12),
Fringe conditions is converted into:
In addition, Ni etc. (2005) propose the relation between saturated-steam temperature and pressure:
Derivative scalar coupling builds: the present invention is that in prediction HTHP Oil/gas Well gas injection process, mass dryness fraction, temperature, pressure distribution propose a kind of differential equation coupling system model considering friction loss and pit shaft-stratum unsteady-state heat transfer:
Interactive quadravalence Runge-Kutta finite difference method: before proposition solving model, partial parameters is handled as follows:
(1) vapor (steam) velocity V
m.
(2) vapour density Pm. is because water vapour air-flow is biphase gas and liquid flow, and we apply the averag density that Beggs-Brill (Beggs and Brill, 1973) method calculates mixture.
(3) frictional force T
f.It can be calculated by following equations:
(4) the friction factor f of gas-liquid mixture
m.F
mthe function about Reynolds number and absolute roughness e,
(5) thermal transmittance U
to(from the diverse location of pit shaft axle to the second contact surface).
The right half part of equation (20) is the thermal resistance of insulation tube, is respectively target chamber and cement sheath.λ
insand λ
cembe respectively the thermal conductivity of thermal insulation material and cement sheath.Hc and hr is respectively convective heat-transfer coefficient and radiation heat transfer coefficient.
(6) nondimensional time function.It can be calculated by equation (21).
Because differential equation coupling system model not only comprises ODE, and comprise partial differential equation, so we propose interactive quadravalence Runge-Kutta finite difference method solve problem.Detailed algorithm can be summarized as with under type:
Step 1. provides θ respectively
0, T
0, P
0, X
0, Q
0and T
e0initial value.
Step 2. calculates all coefficients in coupled system.
Step 3. makes the differential method be converted into function f
i(i=1,2,3,4).Then we obtain following coupling function system:
The equation group that step 4. solves above obtains P
k, X
k, Q
k, T
k, and T
ek, then we draw coefficient
a
i=f
i(P
k,X
k,Q
k,T
k),j=1,2,3,4。
Step 5. makes
then by separating following equation, we obtain Te,
As time j, radial i when degree of depth z, T
j irepresent temperature, i=1; 2 ... M, j=1,2 ... N.Wherein M and N represents time and last radial node respectively.We apply finite difference method and carry out discrete obtaining to equation (23):
In this equation, τ represents the time interval, and ζ represents spaced radial.It can convert canonical form to, as follows:
Then calculus of finite differences discrete boundary condition is used.If rD=1, then:
It meets
If rD=N, then:
In conjunction with equation (25), (27) and (28), we can calculate the numerical solution of formation temperature Te.
Step 6. is as rD=1, and Te is substituted into equation system (23), and we obtain:
The similar method of step 7. calculates
With
d
j=f
j(P
k+hc
1,X
k+hc
2,Q
k+hc
3,T
k+hc
4)(j=1,2,3,4)
Step 8. calculates gas-liquid mixture in (K+1) this point, and mass dryness fraction, the pressure and temperature of gas can obtain:
Step 9. is T=T
k+1substitute into the fringe conditions of equation (23), by finite difference method, we can obtain T
e, k+1.
Step 10. repeats step 2-9, until calculate P
nx
nq
n, T
nand T
en.
Embodiment: for Y gas well at HTHP, utilizes the above model set up, calculates parameters such as the mass dryness fraction in gas injection process, pressure, temperature.As described in above-mentioned model analysis and solution procedure, well section is split into some unit from bottom, the time also will do similar division, then calculates according to above-mentioned calculation procedure.
Relevant parameter and survey data: internal flow density is 1000kg/m3; External fluid density is 1000kg/m3; The degree of depth is 7100 meters; Surface temperature is 20c; Ground level heat conductance parameter is 2.06; Surface temperature gradient is 0.0218
1c/m; The length of every section be 1m. in addition, pipe parameter, inclined shaft parameter, inclination, orientation and vertical depth parameter, respectively at table 1, are listed in table 2 and table 3.
Table 1
Table 2
Table 3
Result of calculation trend analysis: based on Runge-Kutta method and finite difference method, we obtain a series of result about this well.Pressure and temperature and the formation temperature of steam quality, pit shaft are as shown in table 4.
Table 4
Sensitivity analysis: main result as shown in FIG. 4,5,6, 7.The pressure of steam, temperature and mass dryness fraction increase along with the increase of injection length.In fact, it becomes close to actual conditions.Along with the increase of injecting quantity of steam, the water in pit shaft reduces gradually, so the aridity of steam increases.Meanwhile, along with the increase of heat transfer, the temperature on stratum also can increase, but the speed increased becomes slow.In order to study gas-liquid mixture injection rate and surface temperature gradient to the impact of gas aridity, pressure and temperature, we use different gas-liquid mixture injection rates and surface temperature gradient to reform former calculating.
Claims (5)
1. mass dryness fraction, temperature and coupling pressure Forecasting Methodology in Oil/gas Well gas injection process, is characterized in that, comprising:
A. stratum heat transfer model is built;
B. differential equation coupling model is built;
C. differential equation coupling model is solved.
2. mass dryness fraction, temperature and coupling pressure Forecasting Methodology in Oil/gas Well gas injection process as claimed in claim 1, is characterized in that,
Based on from the second contact surface to the hypothesis of stratum unsteady heat transfer in steps A, build stratum heat transfer model, specifically comprise:
Primary condition:
T
e=T
α+ γ z, if t=0.
Fringe conditions:
Nondimensional number
With
Substitute into equation (12),
Fringe conditions is converted into:
Relation between saturated saturated-steam temperature and pressure:
3. mass dryness fraction, temperature and coupling pressure Forecasting Methodology in Oil/gas Well gas injection process as claimed in claim 2, is characterized in that,
In step B, described differential equation coupling model considers friction loss and pit shaft-stratum unsteady-state heat transfer, and the differential equation coupling model of structure is:
4. mass dryness fraction, temperature and coupling pressure Forecasting Methodology in Oil/gas Well gas injection process as claimed in claim 3, is characterized in that,
In step C, before solving model, the parameter in model is handled as follows:
1) vapor (steam) velocity V
m:
2) vapour density Pm: because water vapour air-flow is biphase gas and liquid flow, therefore applies the averag density that Beggs-Brill method calculates mixture;
3) frictional force T
f:
4) the friction factor f of gas-liquid mixture
m, f
mthe function about Reynolds number and absolute roughness e,
5) thermal transmittance U
to:
Wherein, λ
insand λ
cembe respectively the thermal conductivity of thermal insulation material and cement sheath; Hc and hr is respectively convective heat-transfer coefficient and radiation heat transfer coefficient;
6) nondimensional time function:
5. mass dryness fraction, temperature and coupling pressure Forecasting Methodology in Oil/gas Well gas injection process as claimed in claim 4, is characterized in that,
Adopt interactive quadravalence Runge-Kutta finite difference method to solve differential equation coupling model, comprising:
Step 1. provides θ respectively
0, T
0, P
0, X
0, Q
0and T
e0initial value;
Step 2. calculates all coefficients in coupling model;
Step 3. makes the differential method be converted into function f
i(i=1,2,3,4), can obtain following coupling function system:
The equation group that step 4. solves above obtains P
k, X
k, Q
k, T
k, and T
ek, can coefficient be drawn:
a
j=f
i(P
k,X
k,Q
k,T
k),j=1,2,3,4;
Step 5. makes
te is obtained by separating following equation,
As time j, radial i when degree of depth z, T
j irepresent temperature, i=1; 2 ... M, j=1,2 ... N; Wherein M and N represents time and last radial node respectively; Finite difference method is adopted to carry out discrete obtaining to equation (23):
In this equation, τ represents the time interval, and ζ represents spaced radial; It can convert canonical form to, as follows:
Then calculus of finite differences discrete boundary condition is used; If rD=1, then:
It meets
If rD=N, then:
In conjunction with equation (25), (27) and (28), the numerical solution of formation temperature Te can be calculated;
Step 6., as rD=1, substitutes into equation system (23) Te, can obtain:
Step 7. calculates
With
Step 9. is T=T
k+1substitute into the fringe conditions of equation (23), can T be obtained by finite difference method
e, k+1;
Step 10. repeats step 2-9, until calculate P
nx
nq
n, T
nand T
en.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510293979.0A CN104989351A (en) | 2015-06-01 | 2015-06-01 | Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510293979.0A CN104989351A (en) | 2015-06-01 | 2015-06-01 | Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104989351A true CN104989351A (en) | 2015-10-21 |
Family
ID=54301220
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510293979.0A Pending CN104989351A (en) | 2015-06-01 | 2015-06-01 | Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104989351A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105389439A (en) * | 2015-12-03 | 2016-03-09 | 四川大学 | Natural gas well shut-in wellbore modeling simulation implementation method |
CN105484687A (en) * | 2015-11-27 | 2016-04-13 | 中国海洋石油总公司 | Perforation transient well fluid pressure prediction method |
CN105512433A (en) * | 2016-01-12 | 2016-04-20 | 北京航空航天大学 | Fluid-solid node two-phase flow modeling method |
CN107133373A (en) * | 2016-02-29 | 2017-09-05 | 中国石油化工股份有限公司 | A kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network |
CN107463761A (en) * | 2016-06-02 | 2017-12-12 | 中国石油天然气股份有限公司 | gas storage gas injection well control dynamic evaluation method |
CN114638089A (en) * | 2022-02-25 | 2022-06-17 | 海仿(上海)科技有限公司 | Heat transfer calculation method applied to deep sea oil and gas well |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6028992A (en) * | 1996-11-18 | 2000-02-22 | Institut Francais Du Petrole | Method for constituting a model representative of multiphase flows in oil production pipes |
CN101769147A (en) * | 2010-01-12 | 2010-07-07 | 蒋海岩 | Method for evaluating oilfield fireflood scheme |
-
2015
- 2015-06-01 CN CN201510293979.0A patent/CN104989351A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6028992A (en) * | 1996-11-18 | 2000-02-22 | Institut Francais Du Petrole | Method for constituting a model representative of multiphase flows in oil production pipes |
CN101769147A (en) * | 2010-01-12 | 2010-07-07 | 蒋海岩 | Method for evaluating oilfield fireflood scheme |
Non-Patent Citations (2)
Title |
---|
JIUPING XU,ET AL: "A Direct Eulerian GRP Scheme for the Prediction of Gas-Liquid", 《ABSTRACT AND APPLIED ANALYSIS》 * |
JIUPING XU,ET AL: "Numerical modelling of steam quality in deviated wells with variable(T, P) fields", 《CHEMICAL ENGINEERING SCIENCE》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105484687A (en) * | 2015-11-27 | 2016-04-13 | 中国海洋石油总公司 | Perforation transient well fluid pressure prediction method |
CN105484687B (en) * | 2015-11-27 | 2018-07-06 | 中国海洋石油总公司 | A kind of Forecasting Methodology of the instantaneous well liquid pressure of perforation |
CN105389439A (en) * | 2015-12-03 | 2016-03-09 | 四川大学 | Natural gas well shut-in wellbore modeling simulation implementation method |
CN105512433A (en) * | 2016-01-12 | 2016-04-20 | 北京航空航天大学 | Fluid-solid node two-phase flow modeling method |
CN105512433B (en) * | 2016-01-12 | 2019-04-09 | 北京航空航天大学 | The node two phase flow modeling method of fluid-solid |
CN107133373A (en) * | 2016-02-29 | 2017-09-05 | 中国石油化工股份有限公司 | A kind of coupled simulation method of shale gas reservoir, pit shaft and surface pipeline network |
CN107463761A (en) * | 2016-06-02 | 2017-12-12 | 中国石油天然气股份有限公司 | gas storage gas injection well control dynamic evaluation method |
CN114638089A (en) * | 2022-02-25 | 2022-06-17 | 海仿(上海)科技有限公司 | Heat transfer calculation method applied to deep sea oil and gas well |
CN114638089B (en) * | 2022-02-25 | 2024-06-07 | 海仿(上海)科技有限公司 | Heat transfer calculation method applied to deep sea oil-gas well |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104989351A (en) | Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process | |
Sun et al. | A numerical study on the non-isothermal flow characteristics of superheated steam in ground pipelines and vertical wellbores | |
Huang et al. | Experimental investigation on friction factor in pipes with large roughness | |
CN104895560B (en) | A kind of deep water test wellbore pressure, temperature field simulation and Hydrate Prediction method | |
CN102272418B (en) | Method for estimation of sagd process characteristics | |
Gao et al. | Two phase flow heat transfer analysis at different flow patterns in the wellbore | |
Lu et al. | The transient behaviour of CO2 flow with phase transition in injection wells during geological storage–Application to a case study | |
Ding et al. | Simulation and feedback analysis of the temperature field in massive concrete structures containing cooling pipes | |
CN106321065B (en) | Method for quantitatively explaining output profile of horizontal gas well | |
CN106802992A (en) | A kind of natural gas well head flow based on oil nozzle model determines method | |
Salilih et al. | Analysis of double U-tube ground heat exchanger for renewable energy applications with two-region simulation model by combining analytical and numerical techniques | |
CN107575214A (en) | Prediction method of temperature and pressure in well bore for injection-production process | |
CN110377879A (en) | A kind of calculation method of the defeated radius of mixed transportation of oil and gas pipeline room temperature collection | |
Uilhoorn | Evaluating the risk of hydrate formation in CO2 pipelines under transient operation | |
Wei et al. | Drilling parameter optimizing strategies to prevent hydrate decomposition risks | |
Spesivtsev et al. | Study of the Transient Terrain–Induced and Severe Slugging Problems by Use of the Drift–Flux Model | |
CN102108851B (en) | Horizontal well steam injection method and system | |
Xie et al. | Mathematic modeling of the salt deposits growing on the tubing during debrining for gas storage salt cavern | |
Li et al. | Effects of groundwater flow on thermal response test of deep borehole heat exchanger | |
CN104866681A (en) | Temperature pressure numerical simulation method in closing process of high-temperature high-pressure oil gas inclined shaft | |
CN105184061B (en) | The method for numerical simulation of gas-producing well temperature, pressure distribution | |
TWI734029B (en) | Production capacity estimation method for geothermal well | |
CN116738113A (en) | Construction method of mathematical model for predicting scaling position in geothermal well bore | |
Zhang et al. | Study on temperature distribution of perforated horizontal wellbore | |
Hou et al. | Discussion of “rigid water column model for simulating the emptying process in a pipeline using pressurized air” by Oscar E. Coronado-Hernández, Vicente S. Fuertes-Miquel, Pedro L. Iglesias-Rey, and Francisco J. Martínez-Solano |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151021 |