CN112347599B - Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment - Google Patents
Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment Download PDFInfo
- Publication number
- CN112347599B CN112347599B CN201910724223.5A CN201910724223A CN112347599B CN 112347599 B CN112347599 B CN 112347599B CN 201910724223 A CN201910724223 A CN 201910724223A CN 112347599 B CN112347599 B CN 112347599B
- Authority
- CN
- China
- Prior art keywords
- pipe
- column
- measured
- tubular column
- tested
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
Abstract
A polymer flooding concentric double-pipe injection parameter calculation method and a device thereof specifically comprise the steps of respectively calculating and obtaining the relationship between the rheological index and the consistency coefficient of a polymer solution and the temperature of a pipe column to be measured according to the size of the pipe column to be measured, the injection amount of the pipe column to be measured and the related parameters of the polymer solution; calculating to obtain the shearing rate of the pipe column to be measured according to the relation and the injection amount of the pipe column to be measured; calculating according to the injection amount of the pipe column to be detected and the relation to obtain a discrimination value, and obtaining the flow state condition of the pipe column to be detected according to the discrimination value; calculating the on-way friction resistance of the pipe column to be detected through a corresponding flow equation according to the flow state condition of the pipe column to be detected; and calculating to obtain the injection pressure of the wellhead of the tubular column to be tested through the hole friction of the corresponding interval of the tubular column to be tested, the friction pressure in the vertical pipe flow process of the tubular column to be tested, the formation pressure of the corresponding interval of the tubular column to be tested, the injection additional ascending pressure value of the corresponding interval of the tubular column to be tested, the seepage pressure drop of the corresponding interval of the tubular column to be tested and the vertical liquid column pressure of the tubular column to be tested.
Description
Technical Field
The invention relates to the field of chemical flooding oil extraction, in particular to a method and a device for calculating polymer flooding concentric double-pipe injection parameters and computer equipment.
Background
The polymer flooding is characterized in that high molecular weight (large molecular weight) water-soluble polymer is added into water, the viscosity of the water is increased, the rock permeability of a water invasion zone is reduced in the injection process, the oil-water fluidity ratio is improved, and the liquid volume of the injected water in an oil layer is enlarged, so that the recovery rate is improved, and the polymer flooding belongs to the field of chemical flooding oil recovery. In the polymer injection process, injection parameters and an injection column structure are main factors for determining the implementation effect of the polymer flooding. The injection parameters comprise injection quantity, injection viscosity, pressure and other parameters, the injection quantity is generally given by numerical simulation of oil reservoir engineering, and the injection pressure and the injection viscosity need to meet the requirements of ground process, bottom hole injection pressure and bottom hole viscosity loss. The existing chemical flooding injection string comprises a general injection and a layered injection, wherein the layered injection comprises a bridge type eccentric injection and a concentric double pipe. The optimization design method of the bridge type eccentric injection parameters is researched by more technicians, while the optimization design method of the concentric double-tube injection parameters does not exist, the viscosity loss of the concentric double-tube separate injection implemented on site at present can only be completed through an actual measurement mode, if the viscosity loss is large, the viscosity of an injection agent needs to be adjusted, so that great troubles are brought to production management, and the implementation effect of polymer flooding is also influenced. In summary, in order to meet the facility requirements of the polymer flooding concentric double-pipe ground process pressure regulator on pressure difference adjustment, meet the requirements of bottom hole viscosity loss and ensure the implementation effect of polymer flooding, a polymer flooding concentric double-pipe injection parameter calculation method is needed, and the wellhead injection pressure and the bottom hole viscosity loss can be predicted in advance.
Disclosure of Invention
The invention aims to realize the prediction of wellhead injection pressure and bottom viscosity loss in advance and guide the polymer flooding field construction by providing a polymer flooding concentric double-pipe injection parameter calculation method, a polymer flooding concentric double-pipe injection parameter calculation device and polymer flooding concentric double-pipe injection parameter calculation computer equipment.
To achieve the above object, the present invention provides a method for calculating injection parameters of polymer flooding concentric dual pipes, the method comprising: respectively calculating and obtaining the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the size of the to-be-measured pipe column, the injection amount of the to-be-measured pipe column and the related parameters of the polymer solution; calculating according to the relationship between the rheological index and the temperature of the tubular column to be measured and the injection amount of the tubular column to be measured to obtain the shear rate of the tubular column to be measured; calculating according to the injection amount of the pipe column to be detected and the rheological index and the consistency coefficient of the polymer solution to obtain a discrimination value, and obtaining the flow state condition of the pipe column to be detected according to the discrimination value; calculating the on-way friction resistance of the pipe column to be tested through a corresponding flow state equation according to the flow state condition of the pipe column to be tested and the injection amount of the pipe column to be tested; calculating according to the on-way friction resistance of the tubular column to be tested to obtain the hole friction resistance of the corresponding layer section of the tubular column to be tested and the friction resistance pressure of the tubular column to be tested in the vertical pipe flow process; and calculating to obtain the injection pressure of the wellhead of the tubular column to be tested through the friction resistance of the hole of the corresponding section of the tubular column to be tested, the friction resistance pressure in the vertical pipe flow process of the tubular column to be tested, the formation pressure of the corresponding section of the tubular column to be tested, the injection additional ascending pressure value of the corresponding section of the tubular column to be tested, the seepage pressure drop of the corresponding section of the tubular column to be tested and the vertical liquid column pressure of the tubular column to be tested.
In the above method for calculating the injection parameters of the polymer flooding concentric double pipes, preferably, the tubular column to be tested comprises an inner pipe or an annulus between the inner pipe and the outer pipe.
In the above method for calculating the injection parameter of the polymer flooding concentric double-pipe, preferably, the step of calculating the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the pipe column to be measured according to the size of the pipe column to be measured, the injection amount of the pipe column to be measured and the related parameter of the polymer solution comprises: calculating to obtain temperature field data of the tubular column to be measured according to the size of the tubular column to be measured and the injection amount of the tubular column to be measured; and respectively calculating the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the corresponding relationship between the temperature field data of the to-be-measured pipe column and the related parameters of the polymer solution.
In the above method for calculating injection parameters of polymer flooding concentric double pipes, preferably, obtaining the flow state condition of the pipe column to be measured according to the discrimination value includes: calculating to obtain a discrimination value when the power law fluid in the pipe column to be measured flows according to the polymer solution related parameters, the size of the pipe column to be measured, and the rheological index and the consistency coefficient of the polymer solution, and comparing the discrimination value with a preset threshold value; when the flow state of the tubular column to be detected is smaller than the preset threshold value, the flow state of the tubular column to be detected is a laminar flow state; and when the flow state of the tubular column to be detected is larger than the preset threshold value, the flow state of the tubular column to be detected is a turbulent flow state.
In the above method for calculating the injection parameters of the polymer flooding concentric double-pipe, preferably, when the pipe column to be measured is an inner oil pipe, the corresponding flow state equation comprises a laminar flow state friction calculation equation and a turbulent flow state friction calculation equation;
the calculation equation of the laminar flow state friction resistance comprises the following steps:
the calculation equation of the friction resistance of the turbulent flow state comprises the following steps:
wherein:
in the above formula, P Friction resistance inner pipe The friction resistance pressure in the process of vertical pipe flow of the inner oil pipe is MPa; p Friction inner tube i The friction resistance pressure corresponding to unit step length of the inner oil pipe is MPa; a is Inner tube i 、b Inner tube i The Booth-type empirical coefficient corresponding to the ith step length of the inner oil pipe; k Inner tube i For the temperature T in the inner oil pipe Inner tube i The corresponding consistency coefficient, mPa.sn; d Inner pipe Designing the inner diameter of an inner oil pipe to be mm; n is Inner tube i For temperature T in inner oil pipe Inner tube i The corresponding rheological coefficient is treated, and the dimension is not required; v. of Inner tube The average flow velocity of the inner oil pipe is m/s; l is the calculation step length, 100m; rho is the fluid density, kg/m 3 。
In the above method for calculating the injection parameters of the polymer flooding concentric double-pipe, preferably, when the pipe column to be measured is an inner pipe and an outer pipe which are annularly empty, the corresponding flow state equation comprises a laminar flow state friction calculation equation and a turbulent flow state friction calculation equation;
the laminar flow state friction calculation equation comprises:
the calculation equation of the friction resistance in the turbulent flow state comprises the following steps:
in the above formula, P Friction resistance annular space The friction resistance pressure in the annular vertical pipe flow process is MPa; p Frictional annular space i The friction resistance pressure corresponding to the unit step length of the annulus is MPa; a is a Annulus i 、b Annulus i The experimental parameters are Bulazus type empirical coefficients corresponding to the ith step length of the annular space; k is Annulus i For temperature T in annulus Annulus i The corresponding consistency coefficient, mPa.sn; d Annular space Designing the inner diameter of the annulus in mm; n is a radical of an alkyl radical Annulus i For temperature T in annulus Annulus i The corresponding rheological coefficient is treated, and the dimension is not required; v. of Annular space The average flow velocity of the annulus is m/s; l is the calculation step length, 100m; rho is the fluid density, kg/m 3 。
The present invention also provides a polymer flooding concentric double pipe injection parameter calculation apparatus, comprising: the relation module is used for respectively calculating and obtaining the relation between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the size of the to-be-measured pipe column, the injection amount of the to-be-measured pipe column and the related parameters of the polymer solution; the shear rate calculation module is used for calculating and obtaining the shear rate of the tubular column to be measured according to the relation between the rheological index and the temperature of the tubular column to be measured and the injection amount of the tubular column to be measured; the flow state judgment module is used for calculating to obtain a discrimination value according to the injection amount of the pipe column to be detected and the rheological index and the consistency coefficient of the polymer solution, and obtaining the flow state condition of the pipe column to be detected according to the discrimination value; the friction calculation module is used for calculating the on-way friction of the pipe column to be measured through a corresponding flow state equation according to the flow state condition of the pipe column to be measured and the injection amount of the pipe column to be measured; the pressure prediction module is used for calculating and obtaining the hole friction resistance of the corresponding layer section of the tubular column to be tested and the friction resistance pressure of the tubular column to be tested in the vertical pipe flow process according to the on-way friction resistance of the tubular column to be tested; and calculating to obtain the injection pressure of the wellhead of the tubular column to be tested through the hole friction of the corresponding interval of the tubular column to be tested, the friction pressure in the vertical pipe flow process of the tubular column to be tested, the formation pressure of the corresponding interval of the tubular column to be tested, the injection additional ascending pressure value of the corresponding interval of the tubular column to be tested, the seepage pressure drop of the corresponding interval of the tubular column to be tested and the vertical liquid column pressure of the tubular column to be tested.
In the above polymer flooding concentric double-pipe injection parameter calculation apparatus, preferably, the pipe string to be tested includes an inner oil pipe or an inner and outer pipe annulus.
In the above apparatus for calculating polymer flooding concentric double pipe injection parameters, preferably, the relationship module further comprises: calculating according to the size of the pipe column to be measured and the injection amount of the pipe column to be measured to obtain temperature field data of the pipe column to be measured; and respectively calculating the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-detected pipe column according to the corresponding relationship between the temperature field data of the to-be-detected pipe column and the related parameters of the polymer solution.
In the above apparatus for calculating the polymer flooding concentric dual pipe injection parameter, preferably, the flow state determination module further comprises: calculating to obtain a discrimination value when the power law fluid in the pipe column to be measured flows according to the polymer solution related parameters, the size of the pipe column to be measured, and the rheological index and the consistency coefficient of the polymer solution, and comparing the discrimination value with a preset threshold value; when the flow state of the tubular column to be detected is smaller than the preset threshold value, the flow state of the tubular column to be detected is a laminar flow state; and when the flow state of the tubular column to be detected is larger than the preset threshold value, the flow state of the tubular column to be detected is a turbulent flow state.
The invention also provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method when executing the computer program.
The present invention also provides a computer-readable storage medium storing a computer program for executing the above method.
The beneficial technical effects of the invention are as follows: according to the parameters of the given injection amount, the polymer concentration, the polymer molecular weight, the sizes of the inner pipe and the outer pipe and the like, the injection pressure of a wellhead and the viscosity loss in a shaft are predicted in advance, the guidance on site construction is realized, and the oil displacement effect of the polymer is ensured.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
fig. 1 is a schematic flow chart illustrating a method for calculating polymer flooding concentric dual pipe injection parameters according to an embodiment of the present invention;
fig. 2 is a schematic flow chart of an inner oil pipe in the polymer flooding concentric dual-pipe injection parameter calculation method according to an embodiment of the present invention;
FIG. 3 is a schematic flow chart of an inner and outer pipe annulus in a method for calculating injection parameters of a polymer flooding concentric double pipe according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of a polymer flooding concentric dual-tube injection parameter calculation apparatus according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the present invention is further described in detail below with reference to the embodiments and the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
Referring to fig. 1, the method for calculating injection parameters of polymer flooding concentric dual-pipe provided by the present invention comprises: s101, respectively calculating according to the size of a pipe column to be measured, the injection amount of the pipe column to be measured and related parameters of the polymer solution to obtain the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the pipe column to be measured; s102, calculating according to the relation between the rheological index and the temperature of the tubular column to be tested and the injection amount of the tubular column to be tested to obtain the shear rate of the tubular column to be tested; s103, calculating according to the injection amount of the column to be tested and the rheological index and the consistency coefficient of the polymer solution to obtain a discrimination value, and obtaining the flow state condition of the column to be tested according to the discrimination value; s104, calculating according to the flow state condition of the pipe column to be measured and the injection allocation amount of the pipe column to be measured through a corresponding flow state equation to obtain the on-way friction resistance of the pipe column to be measured; s105, calculating according to the on-way friction resistance of the tubular column to be tested to obtain the hole friction resistance of the corresponding layer section of the tubular column to be tested and the friction resistance pressure in the vertical pipe flow process of the tubular column to be tested; s106, calculating to obtain the injection pressure of the wellhead of the tubular column to be tested through the hole friction of the corresponding section of the tubular column to be tested, the friction pressure in the vertical tubular flow process of the tubular column to be tested, the formation pressure of the corresponding section of the tubular column to be tested, the injection additional ascending pressure value of the corresponding section of the tubular column to be tested, the seepage pressure drop of the corresponding section of the tubular column to be tested and the vertical liquid column pressure of the tubular column to be tested; the tubular column to be tested comprises an inner oil pipe or an inner and outer pipe annulus.
In the above embodiment, the step S101 of respectively calculating and obtaining the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the column to be measured according to the size of the column to be measured, the dispensing amount of the column to be measured, and the related parameters of the polymer solution includes: calculating according to the size of the pipe column to be measured and the injection amount of the pipe column to be measured to obtain temperature field data of the pipe column to be measured; and respectively calculating the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the corresponding relationship between the temperature field data of the to-be-measured pipe column and the related parameters of the polymer solution. The shear rate in the step S102 is used for determining the injection pressure of the annular wellhead, and secondly, the shear rate of the tubular column to be tested is determined; since the method of use of the shear rate is a common technique in the art after it is achieved, the present invention will not be described in detail herein.
In step S103, obtaining the flow state condition of the string to be measured according to the discrimination value includes: calculating to obtain a discrimination value when the power law fluid in the pipe column to be measured flows according to the polymer solution related parameters, the size of the pipe column to be measured, and the rheological index and the consistency coefficient of the polymer solution, and comparing the discrimination value with a preset threshold value; when the flow state of the tubular column to be detected is smaller than the preset threshold value, the flow state of the tubular column to be detected is a laminar flow state; and when the flow state of the tubular column to be detected is larger than the preset threshold value, the flow state of the tubular column to be detected is a turbulent flow state.
Specifically, when the pipe column to be tested is an inner oil pipe, the corresponding flow state equation comprises a laminar flow state friction calculation equation and a turbulent flow state friction calculation equation;
the calculation equation of the laminar flow state friction resistance comprises the following steps:
the calculation equation of the friction resistance in the turbulent flow state comprises the following steps:
wherein:
in the above formula, P Friction resistance inner pipe The friction resistance pressure in the vertical pipe flow process of the inner oil pipe is MPa; p Friction inner tube i The friction resistance pressure corresponding to the unit step length of the inner oil pipe is MPa; a is Inner tube i 、b Inner tube i The Booth-type empirical coefficient corresponding to the ith step length of the inner oil pipe; k Inner tube i For temperature T in inner oil pipe Inner tube i The corresponding consistency coefficient, mPa.sn; d Inner pipe Designing the inner diameter of an inner oil pipe to be mm; n is Inner tube i For temperature T in inner oil pipe Inner tube i Corresponding rheological coefficient, dimensionless；v Inner pipe The average flow velocity of the inner oil pipe is m/s; l is the calculation step length, 100m; rho is fluid density, kg/m 3 。
When the pipe column to be tested is empty in an inner pipe and an outer pipe, the corresponding flow state equation comprises a laminar flow state friction resistance calculation equation and a turbulent flow state friction resistance calculation equation;
the calculation equation of the laminar flow state friction resistance comprises the following steps:
the calculation equation of the friction resistance in the turbulent flow state comprises the following steps:
in the above formula, P Friction resistance annular space The friction resistance pressure in the annular vertical pipe flow process is MPa; p Frictional annular space i The friction resistance pressure corresponding to the unit step length of the annulus is MPa; a is Annulus i 、b Annulus i The experimental parameters are Brachiis type empirical coefficients corresponding to the ith step length of the annulus; k is Annulus i For temperature T in annulus Annulus i The corresponding consistency coefficient, mPa.sn; d Annular space Designing the inner diameter of the annulus in mm; n is Annulus i For temperature T in annulus Annulus i The corresponding rheological coefficient is treated, and the dimension is not required; v. of Annular space The average flow velocity of the annulus is m/s; l is the calculation step length, 100m; rho is the fluid density, kg/m 3 。
In order to facilitate understanding of the method for calculating the injection parameters of the polymer flooding concentric dual-tube, the following takes a specific actual calculation flow as an example to integrally describe the above embodiments, specifically as follows:
the calculation method of the rated polymer flooding concentric double-pipe injection parameter provided by the invention actually comprises an inner pipe injection parameter calculation method and an inner and outer pipe annular injection parameter calculation method.
The method for calculating the inner tube injection parameters can be specifically shown in fig. 2, the whole method can be divided into six steps, and the detailed flow is as follows:
1. the inner tube injection parameter calculation method comprises the following steps:
(1) According to the given injection amount and the size of the inner oil pipe, calculating the temperature distribution of different positions in the inner pipe to obtain the temperature field of the inner oil pipe; the related formula comprises:
the temperature model of the oil pipe from the well mouth to the section I of the injection layer is as follows:
C＝2πr xo U xo /C p1 M 1 ；
in the formula:
-inner and outer pipe annulus wellhead fluid temperature, K;
-the well head fluid temperature in the inner tubing, K; v. of 1 -the fluid injection speed in the inner tubing, m/s; r is xo -the outer radius of the inner tube, m; u shape xo -heat transfer coefficient from inner and outer pipe annuli to inner oil pipe, J/(s m) 2 ·K)；C p1 The constant pressure specific heat capacity of the fluid injected into the inner oil pipe is J/(kg. K); m 1 -the mass flow of the fluid injected by the inner oil pipe is kg/s.
An inner oil pipe temperature model from an injection layer I section to an injection layer II section: after the injection zone I section is calculated, only the inner pipe is present when the injection zone II section is to be calculated, and the temperature of the rest well sections is calculated according to the following formula. According to the energy conservation equation, the temperature distribution of the inner oil pipe from the injection layer I to the injection layer II is as follows:
in the formula:
-formation temperature, K, at injection layer i;
the temperature of the fluid in the inner oil pipe at the injection layer I is K; u shape ho -heat transfer coefficient from inside the inner tubing to outside the cement sheath, J/(s · m) 2 ·K)；r xo -the outer radius of the inner tube, m; u shape xo -heat transfer coefficient from inner and outer pipe annuli to inner oil pipe, J/(s m) 2 ·K)；C p1 The constant pressure specific heat capacity of the fluid injected into the inner oil pipe is J/(kg. K); m 1 The mass flow of the injected fluid of the inner oil pipe is kg/s; k is a radical of formula e -formation thermal conductivity, J/(s.m.K).
(2) And obtaining the relation between the consistency coefficient of the polymer solution to be injected into the target block and the temperature of the inner tube and the relation between the rheological index and the temperature of the inner tube through an indoor test.
The apparent viscosity of the polymer fluid is:
k (consistency factor) and n (rheology index) were determined by laboratory experiments with the unit temperature field variation as an iterative variable.
k Inner pipe ＝aT Inner pipe +b；
n Inner tube ＝cT Inner pipe +d；
(3) And calculating the shear rate of the inner oil pipe according to the given inner pipe injection amount and the rheological index.
In the formula:
-the wall of the pipe corresponding to the ith step length in the inner pipeShear rate of (d), s -1 。
(4) And calculating a discrimination value according to the given inner pipe distribution injection amount, the rheological index and the consistency coefficient, and judging the flow state of the inner pipe according to the discrimination value.
K Inner tube i ＝a 1 T Inner tube i +b 1 ；
n Inner tube i ＝a 2 T Inner tube i +b 2 ；
In the formula: z Inner tube i Is a discrimination value and has no dimension; re Inner tube i Is Reynolds number, dimensionless; d Inner tube Designing the inner diameter of the inner pipe to be mm; rho is the fluid density, kg/m 3 ；Q Inner pipe Amount of inner tube dosed, m 3 /d；v Inner pipe The average flow velocity of the inner tube is m/s; k Inner tube i Is the temperature T in the inner tube Inner tube i The corresponding consistency coefficient, mPa.sn; n is Inner tube i Is the temperature T in the inner tube Inner tube i The corresponding rheological coefficient is treated, and the dimension is not required; t is Inner tube i The temperature corresponding to the ith l depth of the inner tube is DEG C; l is the calculation step length, 100m; a is 1 、b 1 、a 2 、b 2 The coefficients of K-T and n-T corresponding to the molecular weight and concentration of the injected polymer are obtained; ryan and Johnson indicate that the critical value of Z for power law fluid flow in the circular tube and annulus is 808, laminar flow for values of Z less than 808, and turbulent flow for values of Z greater than 808.
(5) And calculating the in-process friction resistance of the inner oil pipe according to the given inner pipe injection amount and the flow state of the fluid with the inner pipe diameter.
Under laminar flow conditions:
in a turbulent state:
in the formula: p is Friction resistance inner tube The friction resistance pressure in the process of the inner pipe vertical pipe flow is MPa; p Friction inner tube i The friction resistance pressure corresponding to the unit step length of the inner pipe is MPa; a is Inner tube i 、b Inner tube i And the coefficient is a Brazius empirical coefficient corresponding to the ith step length of the inner tube.
(6) Inner tube wellhead pressure calculation
Pressure drop across the orifice:
in the formula: p Eyelet inner tube The abrasion resistance pressure drop of the inner pipe corresponding to the perforation well section hole is MPa; n is a radical of hydrogen Inner pipe The number of effective holes of the inner pipe corresponding to the perforation well section is dimensionless; d 0 Is the diameter of the perforation hole m; c d inner tube The flow coefficient of the inner pipe corresponding to the perforation well section can be in a range of 0.56-0.6; mu.s s inner tube The viscosity of the inner pipe at the blast hole position corresponding to the perforation well section is mPa & s; k t inner pipe The consistency coefficient of the polymer solution at the blast hole position of the inner pipe corresponding to the perforation well section is mPa.sn; n is t inner pipe The rheological index of the polymer solution at the blast hole position of the inner pipe corresponding to the perforation well section is dimensionless; v 0 inner pipe Is the borehole flow velocity, m/s.
Seepage pressure drop:
in the formula: p Seepage inner pipe The seepage pressure drop of the corresponding layer section of the inner pipe is MPa;
porosity,%, of the inner tube corresponding to the interval; s Inner tube The water absorption index of the corresponding interval of the inner pipe is m < 3 >/(d.m.MPa); l is the injection-production well spacing, m; v is the injection speed of the interval corresponding to the inner pipe, PV/a.
Injection pressure rise value:
λ inner pipe ＝k Inner pipe /μ s inner tube ；
In the formula: p is Ascending inner pipe Is the capillary force and friction resistance rise value caused by injecting non-Newtonian fluid, MPa; h is Inner pipe The effective thickness m of the corresponding interval of the inner pipe is; lambda [ alpha ] Inner pipe The fluidity of the polymer solution in the corresponding interval of the inner pipe is mD/mPa & s; d w Is the inside diameter of the casing, m; k is a radical of Inner pipe Permeability, mD, of the corresponding interval of the inner tube.
Formation pressure:
P inner pipe of stratum ＝c Inner pipe H Inner tube ；
In the formula: p is Inner pipe of stratum The current stratum pressure of the corresponding interval of the inner pipe is MPa; c. C Inner pipe The current formation pressure coefficient of the corresponding interval of the inner pipe is dimensionless; h Inner tube The depth of the corresponding interval of the inner pipe is m.
Wellbore fluid column pressure:
P liquid column inner tube ＝ρgH Inner pipe ；
In the formula: p Liquid column inner pipe The inner tube is perpendicular to the liquid column pressure, MPa.
Predicting the pressure of the inner pipe wellhead:
P wellhead inner pipe ＝P Inner pipe of stratum +P Seepage inner pipe +P Ascending inner pipe +P Perforated inner pipe +P Friction resistance inner pipe -P Liquid column inner tube ；
2. The method for calculating the annular injection parameters of the inner pipe and the outer pipe is specifically shown in fig. 3.
(1) According to the given inner and outer tube annular distribution injection quantity, the size of the inner oil pipe and the size of the outer oil pipe, calculating the temperature distribution of different positions of the inner and outer tube annular to obtain an inner and outer tube annular temperature field:
let the initial liquid temperature:
the formation temperature is generally unaffected by heat transfer and is linearly distributed, i.e.
In the formula:
the temperature of the fluid at the annular wellhead of the inner pipe and the outer pipe is K;
is the wellhead fluid temperature in the inner tubing, K; c P2 Injecting fluid into the annular space of the inner pipe and the outer pipe with constant pressure and specific heat capacity, J/(kg.K); k is a radical of e Is the formation thermal conductivity, J/(s.m.K);
is the formation temperature, K; m 2 Injecting fluid mass flow in kg/s into the inner and outer pipe annular spaces; f (t) is the transient heat transfer function; v. of 2 The annular fluid injection speed of the inner pipe and the outer pipe is m/s; r is to Is the outer radius of the outer oil pipe, m; u shape to The heat transfer coefficient from the inner and outer tube ring spaces to the outside of the cement sheath is J/(s.m) 2 ·K)；M 1 Injecting fluid mass flow, kg/s, into the inner oil pipe; r is xo Is the outer radius of the inner oil pipe, m; u shape xo The heat transfer coefficient from the inner pipe annulus to the outer pipe annulus to the inner oil pipe is J/(s.m) 2 ·K)；T D Is the ground temperature, K; g T Is the geothermal gradient, K/m.
(2) And obtaining the relation between the consistency coefficient of the polymer solution to be injected into the target block and the annular temperature and the relation between the rheological index and the annular temperature through an indoor test.
Apparent viscosity of the polymer fluid is
Determining k by using the unit temperature field change as an iteration variable through an indoor experiment Annular space (coefficient of consistency) and n Annular space (rheology index).
k Annular space ＝aT Annular space +b；
n Annular space ＝cT Annular space +d；
(3) And obtaining the shearing rate of the inner and outer oil pipe annuluses according to the given inner and outer pipe annulus dosage and rheological index.
In the formula:
shear rate at the pipe wall corresponding to the ith step length in the annulus, s -1 。
(4) And judging the flow state of the inner and outer pipe annuluses according to the given calculation results of the inner and outer pipe annulus distribution injection amount, the rheological index and the consistency coefficient.
K Outer tube i ＝a 1 T Outer tube i +b 1 ；
n Outer tube i ＝a 2 T Outer tube i +b 2 ；
In the formula: z Annulus i Is a discrimination value and has no dimension; re Annulus i Is Reynolds number, dimensionless; d is a radical of Outer tube Designing the inner diameter of the outer pipe to be mm; d Inside and outside Designing the outer diameter of the inner pipe to be mm; q Annular space For annular space dosage, m 3 /d；v Annular space The average flow velocity of the annulus is m/s; k Annulus i For temperature T in annulus Annulus i The corresponding consistency coefficient, mPa.sn; n is Annulus i For temperature T in annulus Annulus i The corresponding rheological index is zero dimension; t is a unit of Annulus i The temperature corresponding to the ith depth of the annulus is in DEG C; ryan and Johnson indicate that the critical value of Z for power law fluid flow in the circular tube and annulus is 808, laminar flow for values of Z less than 808, and turbulent flow for values of Z greater than 808.
(5) And calculating to obtain the in-process friction resistance according to the given annular injection amount of the inner pipe and the outer pipe and the flow state of the fluid under the annular size.
Under laminar flow conditions:
in a turbulent state:
in the formula: p Friction resistance annular space The friction resistance pressure in the annular vertical pipe flow process is MPa; p Frictional annular space i The friction resistance pressure corresponding to the unit step length of the annulus is MPa; a is Annulus i 、b Annulus i And obtaining a Brazius empirical coefficient corresponding to the ith step length of the annulus.
(6) Annular wellhead pressure calculation
Pressure drop across the orifice:
in the formula: p Perforated annular space Friction resistance pressure drop of the annular corresponding perforation well section hole is MPa; n is a radical of Annular space The number of effective holes of the annular corresponding perforation well section is dimensionless; d 0 M is the diameter of the perforation hole; c d ring space The flow coefficient of the corresponding perforation well section of the annular space can be in a range of 0.56-0.6; mu.s s-ring space The viscosity of the annular space corresponding to the position of a blast hole of the perforation well section is mPa & s; k is t-ring space The consistency coefficient of the polymer solution at the blast hole position of the annular corresponding perforation well section is mPa.sn; n is t-ring space The rheological index of the polymer solution at the blast hole position of the annular corresponding perforation well section is dimensionless; v 0 annular space The annular borehole flow rate is m/s.
Seepage pressure drop:
in the formula: p Seepage annulus The seepage pressure drop of the corresponding layer section of the annulus is MPa; phi is a Inner pipe Porosity of the corresponding interval of the annulus,%; s Annular space The water absorption index of the corresponding interval of the annulus is m < 3 >/(d.m.MPa); v Annular space Injecting speed, PV/a, for the interval corresponding to the annulus;
injection pressure rise value:
λ annular space ＝k Annular space /μ s-ring space ；
In the formula: p Ascending annulus Is the capillary force and friction resistance rise value caused by injecting non-Newtonian fluid, MPa; h is Annular space The effective thickness of the corresponding layer section of the annulus, m; lambda [ alpha ] Annular space The polymer solution fluidity, mD/mPa & s, of the corresponding interval of the annulus; d is a radical of w Is the inside diameter of the casing, m; k is a radical of Annular space And permeability of the corresponding interval of the annulus, mD.
A formation pressure model:
P formation annulus ＝c Annular space H Annular space ；
In the formula: p Formation annulus The current stratum pressure of the interval corresponding to the inner pipe is MPa; c. C Annular space The current formation pressure coefficient of the corresponding interval of the inner pipe is dimensionless; h Annular space The depth of the corresponding interval of the inner pipe is m.
Wellbore fluid column pressure model:
P liquid column annulus ＝ρgH Annular space ；
In the formula: p Liquid column annulus The inner tube is under vertical liquid column pressure, MPa.
Predicting the wellhead pressure of the inner and outer pipe annular space:
P well head annular space ＝P Formation annulus +P Seepage annulus +P Ascending annulus +P Perforated annular space +P Friction resistance annular space -P Liquid column annulus ；
Referring to fig. 4, the present invention further provides a polymer flooding concentric dual pipe injection parameter calculation apparatus, comprising: the relation module is used for respectively calculating and obtaining the relation between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the size of the to-be-measured pipe column, the injection amount of the to-be-measured pipe column and the related parameters of the polymer solution; the shear rate calculation module is used for calculating and obtaining the shear rate of the tubular column to be measured according to the relation between the rheological index and the temperature of the tubular column to be measured and the injection amount of the tubular column to be measured; the flow state judgment module is used for calculating to obtain a discrimination value according to the injection amount of the pipe column to be detected and the rheological index and the consistency coefficient of the polymer solution, and obtaining the flow state condition of the pipe column to be detected according to the discrimination value; the friction calculation module is used for calculating the on-way friction of the pipe column to be measured through a corresponding flow state equation according to the flow state condition of the pipe column to be measured and the injection amount of the pipe column to be measured; the pressure prediction module is used for calculating and obtaining the hole friction resistance of the corresponding layer section of the tubular column to be tested and the friction resistance pressure of the tubular column to be tested in the vertical pipe flow process according to the on-way friction resistance of the tubular column to be tested; and calculating to obtain the injection pressure of the wellhead of the tubular column to be tested through the friction resistance of the hole of the corresponding section of the tubular column to be tested, the friction resistance pressure in the vertical pipe flow process of the tubular column to be tested, the formation pressure of the corresponding section of the tubular column to be tested, the injection additional ascending pressure value of the corresponding section of the tubular column to be tested, the seepage pressure drop of the corresponding section of the tubular column to be tested and the vertical liquid column pressure of the tubular column to be tested. Wherein, the tubular column that awaits measuring contains interior oil pipe or interior outer tube annular space.
In an embodiment of the present invention, the relationship module further includes: calculating to obtain temperature field data of the tubular column to be measured according to the size of the tubular column to be measured and the injection amount of the tubular column to be measured; and respectively calculating the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the corresponding relationship between the temperature field data of the to-be-measured pipe column and the related parameters of the polymer solution.
In an embodiment of the present invention, the flow state determining module further includes: calculating to obtain a discrimination value when the power law fluid in the pipe column to be measured flows according to the polymer solution related parameters, the size of the pipe column to be measured, and the rheological index and the consistency coefficient of the polymer solution, and comparing the discrimination value with a preset threshold value; when the flow state of the tubular column to be detected is smaller than the preset threshold value, the flow state of the tubular column to be detected is a laminar flow state; and when the flow state of the tubular column to be tested is greater than the preset threshold value, the flow state of the tubular column to be tested is in a turbulent flow state.
The invention also provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method when executing the computer program.
The present invention also provides a computer-readable storage medium storing a computer program for executing the above method.
The beneficial technical effects of the invention are as follows: according to the parameters of the given injection amount, the polymer concentration, the polymer molecular weight, the sizes of the inner pipe and the outer pipe and the like, the injection pressure of a wellhead and the viscosity loss in a shaft are predicted in advance, the guidance on site construction is realized, and the polymer oil displacement effect is ensured.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (12)
1. A method for calculating injection parameters of a polymer flooding concentric dual pipe, the method comprising:
respectively calculating and obtaining the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the size of the to-be-measured pipe column, the injection amount of the to-be-measured pipe column and the related parameters of the polymer solution;
calculating according to the relationship between the rheological index and the temperature of the tubular column to be measured and the injection amount of the tubular column to be measured to obtain the shear rate of the tubular column to be measured;
calculating according to the injection amount of the pipe column to be detected and the rheological index and the consistency coefficient of the polymer solution to obtain a discrimination value, and obtaining the flow state condition of the pipe column to be detected according to the discrimination value;
calculating the on-way friction resistance of the pipe column to be detected through a corresponding flow state equation according to the flow state condition of the pipe column to be detected and the injection amount of the pipe column to be detected;
calculating according to the on-way friction resistance of the tubular column to be tested to obtain the eyelet friction resistance of the corresponding section of the tubular column to be tested and the friction resistance pressure of the tubular column to be tested in the vertical pipe flow process;
and calculating to obtain the injection pressure of the wellhead of the tubular column to be tested through the friction resistance of the hole of the corresponding section of the tubular column to be tested, the friction resistance pressure in the vertical pipe flow process of the tubular column to be tested, the formation pressure of the corresponding section of the tubular column to be tested, the injection additional ascending pressure value of the corresponding section of the tubular column to be tested, the seepage pressure drop of the corresponding section of the tubular column to be tested and the vertical liquid column pressure of the tubular column to be tested.
2. The method of calculating polymer flooding concentric dual pipe injection parameters of claim 1, wherein the pipe string to be tested comprises an inner pipe or an inner and outer pipe annulus.
3. The method for calculating the polymer flooding concentric double-pipe injection parameter according to claim 2, wherein the step of respectively calculating and obtaining the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the pipe column to be measured according to the size of the pipe column to be measured, the injection amount of the pipe column to be measured and the related parameter of the polymer solution comprises the following steps:
calculating to obtain temperature field data of the tubular column to be measured according to the size of the tubular column to be measured and the injection amount of the tubular column to be measured;
and respectively calculating the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the corresponding relationship between the temperature field data of the to-be-measured pipe column and the related parameters of the polymer solution.
4. The method for calculating the polymer flooding concentric double-pipe injection parameter according to claim 2, wherein the step of obtaining the flow state condition of the pipe column to be measured according to the discrimination value comprises the following steps:
calculating to obtain a discrimination value when the power law fluid in the pipe column to be measured flows according to the polymer solution related parameters, the size of the pipe column to be measured, and the rheological index and the consistency coefficient of the polymer solution, and comparing the discrimination value with a preset threshold value; when the flow state of the tubular column to be detected is smaller than the preset threshold value, the flow state of the tubular column to be detected is a laminar flow state; and when the flow state of the tubular column to be detected is larger than the preset threshold value, the flow state of the tubular column to be detected is a turbulent flow state.
5. The polymer flooding concentric double-pipe injection parameter calculation method according to claim 4, wherein when the pipe column to be measured is an inner oil pipe, the corresponding flow state equations comprise a laminar flow state friction calculation equation and a turbulent flow state friction calculation equation;
the calculation equation of the laminar flow state friction resistance comprises the following steps:
the calculation equation of the friction resistance of the turbulent flow state comprises the following steps:
wherein:
in the above formula, P Friction resistance inner tube The friction resistance pressure in the vertical pipe flow process of the inner oil pipe is MPa; p Friction inner tube i The friction resistance pressure corresponding to unit step length of the inner oil pipe is MPa; a is Inner tube i 、b Inner tube i The Booth-type empirical coefficient corresponding to the ith step length of the inner oil pipe; k Inner tube i For the temperature T in the inner oil pipe Inner tube i Treating the corresponding consistency coefficient, mPa.sn; d Inner pipe Designing the inner diameter of an inner oil pipe to be mm; n is Inner tube i For temperature T in inner oil pipe Inner tube i The corresponding rheological coefficient is adopted, and the dimension is not required; v. of Inner pipe The average flow velocity of the inner oil pipe is m/s; l is the calculation step length, 100m; rho is fluid density, kg/m 3 。
6. The polymer flooding concentric double-pipe injection parameter calculation method according to claim 4, wherein when the pipe column to be measured is an inner pipe and an outer pipe which are annularly empty, the corresponding flow state equations comprise a laminar flow state friction calculation equation and a turbulent flow state friction calculation equation;
the laminar flow state friction calculation equation comprises:
the calculation equation of the friction resistance in the turbulent flow state comprises the following steps:
in the above formula, P Friction resistance annular space The friction resistance pressure in the annular vertical pipe flow process is MPa; p Frictional annular space i The friction resistance pressure corresponding to the unit step length of the annulus is MPa; a is a Annulus i 、b Annulus i The experimental parameters are Brachiis type empirical coefficients corresponding to the ith step length of the annulus; k Annulus i For temperature T in annulus Annulus i Treating the corresponding consistency coefficient, mPa.sn; d Annular space Designing the inner diameter of the annulus in mm; n is Annulus i For temperature T in annulus Annulus i The corresponding rheological coefficient is treated, and the dimension is not required; v. of Annular space The average flow velocity of the annulus is m/s; l is the calculation step length, 100m; rho is the fluid density, kg/m 3 。
7. A polymer flooding concentric dual tube injection parameter calculation apparatus, comprising:
the relation module is used for respectively calculating and obtaining the relation between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the size of the to-be-measured pipe column, the injection amount of the to-be-measured pipe column and the related parameters of the polymer solution;
the shear rate calculation module is used for calculating and obtaining the shear rate of the tubular column to be measured according to the relation between the rheological index and the temperature of the tubular column to be measured and the injection amount of the tubular column to be measured;
the flow state judgment module is used for calculating to obtain a discrimination value according to the injection amount of the pipe column to be detected and the rheological index and the consistency coefficient of the polymer solution, and obtaining the flow state condition of the pipe column to be detected according to the discrimination value;
the friction calculation module is used for calculating the on-way friction of the pipe column to be measured through a corresponding flow state equation according to the flow state condition of the pipe column to be measured and the injection amount of the pipe column to be measured;
the pressure prediction module is used for calculating and obtaining the eyelet friction resistance of the corresponding section of the tubular column to be tested and the friction resistance pressure of the tubular column to be tested in the vertical pipe flow process according to the on-way friction resistance of the tubular column to be tested; and calculating to obtain the injection pressure of the wellhead of the tubular column to be tested through the hole friction of the corresponding interval of the tubular column to be tested, the friction pressure in the vertical pipe flow process of the tubular column to be tested, the formation pressure of the corresponding interval of the tubular column to be tested, the injection additional ascending pressure value of the corresponding interval of the tubular column to be tested, the seepage pressure drop of the corresponding interval of the tubular column to be tested and the vertical liquid column pressure of the tubular column to be tested.
8. The polymer flooding concentric dual pipe injection parameter calculation device of claim 7, wherein the tubular string under test comprises an inner tubing or an inner and outer tubing annulus.
9. The polymer flooding concentric dual tube injection parameter calculation apparatus of claim 8, wherein the relationship module further comprises:
calculating to obtain temperature field data of the tubular column to be measured according to the size of the tubular column to be measured and the injection amount of the tubular column to be measured;
and respectively calculating the relationship between the rheological index and the consistency coefficient of the polymer solution and the temperature of the to-be-measured pipe column according to the corresponding relationship between the temperature field data of the to-be-measured pipe column and the related parameters of the polymer solution.
10. The polymer flooding concentric dual pipe injection parameter calculation apparatus according to claim 8, wherein the flow regime determination module further comprises:
calculating to obtain a discrimination value when the power law fluid in the pipe column to be measured flows according to the polymer solution related parameters, the size of the pipe column to be measured, and the rheological index and the consistency coefficient of the polymer solution, and comparing the discrimination value with a preset threshold value; when the flow state of the tubular column to be detected is smaller than the preset threshold value, the flow state of the tubular column to be detected is a laminar flow state; and when the flow state of the tubular column to be detected is larger than the preset threshold value, the flow state of the tubular column to be detected is a turbulent flow state.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 6 when executing the computer program.
12. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910724223.5A CN112347599B (en) | 2019-08-07 | 2019-08-07 | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910724223.5A CN112347599B (en) | 2019-08-07 | 2019-08-07 | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112347599A CN112347599A (en) | 2021-02-09 |
| CN112347599B true CN112347599B (en) | 2022-10-04 |
Family
ID=74366525
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910724223.5A Active CN112347599B (en) | 2019-08-07 | 2019-08-07 | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112347599B (en) |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107066679A (en) * | 2017-03-09 | 2017-08-18 | 中海石油(中国)有限公司 | One kind is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding |
| CN107145696B (en) * | 2017-06-29 | 2020-07-28 | 中国石油大学(北京) | Simulation method for coal bed gas ground and underground coupling solution |
| CN108252696B (en) * | 2017-12-27 | 2021-01-01 | 中国石油天然气股份有限公司 | Screening method of chemical flooding injection oil pipe |
| CN108829945A (en) * | 2018-05-29 | 2018-11-16 | 西南石油大学 | One kind assessing storey increase design volume method based on Fracturing Pressure in real time |
| CN108590634A (en) * | 2018-06-02 | 2018-09-28 | 青岛科技大学 | A kind of shale gas multistage fracturing horizontal well wellbore pressure and inbound traffics prediction technique |
-
2019
- 2019-08-07 CN CN201910724223.5A patent/CN112347599B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112347599A (en) | 2021-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Daccord et al. | Carbonate acidizing: toward a quantitative model of the wormholing phenomenon | |
| CN108131130A (en) | To the analysis method and its device of gas well mouth casing annulus pressure monitoring data | |
| US7536905B2 (en) | System and method for determining a flow profile in a deviated injection well | |
| CN105089639B (en) | A kind of coal bed gas well Wellbore Flow dynamic prediction method | |
| CN104389589B (en) | Method and system of determining shaft temperature field distribution based on hollow-rod | |
| CN102587887B (en) | Forecasting method for gas well condition | |
| Furui et al. | A comprehensive model of high-rate matrix acid stimulation for long horizontal wells in carbonate reservoirs | |
| CN115587674B (en) | Dynamic capacity prediction method for gas well in oil reservoir reconstruction gas storage capacity expansion and production process | |
| RU2460878C2 (en) | Method for determining profile of fluid influx and parameters of borehole environment | |
| Deng et al. | Influence of sand production in an unconsolidated sandstone reservoir in a deepwater gas field | |
| Liu et al. | Design and calculation of process parameters in bullheading and pressurized mud cap drilling | |
| CN104948149A (en) | Multi-media separated injection system suitable for mine field complex oil reservoir | |
| Czarnota et al. | Semianalytical horizontal well length optimization under pseudosteady-state conditions | |
| Krishna et al. | Simplified predictive model for downhole pressure surges during tripping operations using power law drilling fluids | |
| CN113338915B (en) | Method for judging whether gas well accumulates liquid and predicting liquid accumulation height | |
| CN110929447B (en) | Shaft temperature field numerical calculation method in thickening acid acidification process | |
| CN112347599B (en) | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment | |
| Lu et al. | Experimental study of reasonable drawdown pressure of horizontal wells in oil reservoir with bottom water | |
| Chen et al. | Modeling transient circulating mud temperature in the event of lost circulation and its application in locating loss zones | |
| CN112434394B (en) | Method and device for selecting size of polymer flooding concentric double-pipe injection oil pipe | |
| Tong et al. | A Novel Intermittent Gas Lifting and Monitoring System Toward Liquid Unloading for Deviated Wells in Mature Gas Field | |
| CN105184061B (en) | The method for numerical simulation of gas-producing well temperature, pressure distribution | |
| CN107130949A (en) | Deep-layer heavy crude CO_2 stimulation injection process phase-state change method of testing | |
| Abdelhafiz | Modelling and simulating the thermal interaction between the wellbore fluids and rock formation during drilling operations | |
| CN104948175B (en) | A kind of Experimental Method in Laboratory for monitoring interlayer liquid-producing capacity difference |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |