CN107066679A - One kind is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding - Google Patents
One kind is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding Download PDFInfo
- Publication number
- CN107066679A CN107066679A CN201710056182.8A CN201710056182A CN107066679A CN 107066679 A CN107066679 A CN 107066679A CN 201710056182 A CN201710056182 A CN 201710056182A CN 107066679 A CN107066679 A CN 107066679A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- deck
- double
- mfrac
- 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
- 229920000642 polymer Polymers 0.000 title claims abstract description 152
- 238000012360 testing method Methods 0.000 title claims abstract description 62
- 238000004458 analytical method Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000002347 injection Methods 0.000 claims abstract description 50
- 239000007924 injection Substances 0.000 claims abstract description 50
- 238000004088 simulation Methods 0.000 claims abstract description 19
- 238000011161 development Methods 0.000 claims abstract description 9
- 239000000243 solution Substances 0.000 claims description 24
- 230000015572 biosynthetic process Effects 0.000 claims description 17
- 230000035699 permeability Effects 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 7
- 230000001052 transient effect Effects 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 3
- 239000010410 layer Substances 0.000 description 51
- 210000004027 cell Anatomy 0.000 description 9
- 239000012530 fluid Substances 0.000 description 8
- 210000003491 skin Anatomy 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 150000001875 compounds Chemical class 0.000 description 3
- 210000002615 epidermis Anatomy 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000000518 rheometry Methods 0.000 description 3
- 238000010008 shearing Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000006116 polymerization reaction Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 102000003712 Complement factor B Human genes 0.000 description 1
- 108090000056 Complement factor B Proteins 0.000 description 1
- 206010068052 Mosaicism Diseases 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000002355 dual-layer Substances 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Geometry (AREA)
- Geophysics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
It is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding the present invention relates to one kind, specifically includes:1) the bottom pressure value of the double-deck channelling oil reservoir injection well of collection polymer flooding, and calculate the actual bottom pressure variation relation value for obtaining the double-deck channelling oil reservoir injection well of polymer flooding;2) according to the bottom pressure expression formula of the double-deck channelling oil reservoir injection well of polymer flooding and the analog parameter value of input, the simulation bottom pressure variation relation value for obtaining the double-deck channelling oil reservoir injection well of polymer flooding is calculated, and draws the theoretical plate of the double-deck channelling oil reservoir typical curve of polymer flooding;3) according to the well test analysis parameter value that the corresponding double-deck channelling oil reservoir of analog parameter value acquisition polymer flooding of bottom pressure variation relation value is simulated after amendment;4) reservoir situation is driven according to well test analysis parameter value analyzing polymers, takes corresponding measure to instruct oil field development.The present invention can be widely applied in the double-deck channelling oil reservoir well test analysis of polymer flooding.
Description
Technical field
The present invention relates to Well Test Technology field, and in particular to one kind is used for the double-deck channelling oil reservoir well test analysis system of polymer flooding
System and method.
Background technology
The water-drive pool Well Test Data Analysis Method of the well test analysis problem of Polymer Flooding Reservoirs and relative maturity has very
Big difference, its main cause is attributed to inject the rheology sex chromosome mosaicism of fluid, and water is Newtonian fluid, and viscosity does not change, but
It is polymer solution as a kind of non-newtonian fluid, when occurring seepage flow on stratum, not only there is shearing and viscoelastic effect, also simultaneously
There are other physics chemical actions, therefore the rheological characteristic of polymer solution is extremely complex, and then result in polymer flooding well testing solution
Release the difference with water-drive pool well test analysis.Well testing typical curve plots both domestic and external are nearly all at present all have ignored polymer
Interaction between formation rock, including shearing, diffusion, convection current, inaccessible pore volume and permeability decrease system
The problems such as number, existing flow model in porous media is caused can not really to reflect the true seepage flow mechanism of polymer solution in the earth formation, no
True pressure response characteristic of the polymer solution in formation flow can be described.
Actual oil reservoir is made up of the oil reservoir with different physical properties, and interlayer also has channelling, and many scholars are to examining
The stratified reservoir well test model for considering crossflow is studied, but the fluid that these researchs consider is Newtonian fluid, right
The multilayer channelling Study on Problems of non-newtonian fluid such as polymer is less.
The content of the invention
It is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding it is an object of the invention to provide one kind, with poly-
Based on compound rheology characteristic, it is considered to the interaction on polymer solution and stratum, the double-deck channelling oil reservoir of polymer flooding is established
Well test analysis mathematical modeling, the actual seepage flow mechanism of reflection polymer solution in the earth formation improves the double-deck channelling of polymer flooding
The degree of accuracy of oil reservoir well test analysis.
To achieve the above object, the present invention takes following technical scheme:One kind is used for the double-deck channelling oil reservoir examination of polymer flooding
Well analysis method, it is characterised in that comprise the following steps:1) bottom pressure of the double-deck channelling oil reservoir injection well of collection polymer flooding
Value, according to the corresponding acquisition time of bottom pressure value, calculates the real well bottom for obtaining the double-deck channelling oil reservoir injection well of polymer flooding
Pressure change relation value, and draw the actual bottom pressure variation relation curve of the double-deck channelling oil reservoir injection well of polymer flooding;2)
According to the bottom pressure expression formula of the double-deck channelling oil reservoir injection well of polymer flooding and the analog parameter value of input, calculating is polymerize
Thing drives the simulation bottom pressure variation relation value of double-deck channelling oil reservoir injection well, and it is poly- to draw the double-deck channelling oil reservoir injection of polymer flooding
The theoretical plate of the typical curve of well;
Wherein, bottom pressure expression formula is:
In formula, PwfFor bottom pressure, unit is MPa;K is time step, tnFor the n times at moment;
3) according to actual bottom pressure variation relation value, analog parameter value corresponding to simulation bottom pressure variation relation value
It is modified, and altered according to the corresponding analog parameter value acquisition polymer flooding bilayer of bottom pressure variation relation value is simulated after amendment
Flow the well test analysis parameter value of oil reservoir;4) it is poly- according to the well test analysis parameter value analysis of the double-deck channelling oil reservoir of obtained polymer flooding
Compound drives reservoir situation, takes corresponding measure to instruct oil field development.
The step 2) in, drawing the method for the theoretical plate of the double-deck channelling oil reservoir typical curve of polymer flooding includes following step
Suddenly:1. according to the double-deck channelling characteristics of reservoirs of polymer flooding, polymer solution viscosity model and the double-deck channelling oil of polymer flooding are set up
The physical model of Tibetan;2. according to the polymer solution viscosity model of foundation and the double-deck channelling physical models of reservoir of polymer flooding, build
The mathematical modeling of the vertical double-deck channelling oil reservoir of polymer flooding;3. the mathematical modeling of foundation is solved using finite difference method,
Obtain the bottom pressure changing value of the double-deck channelling oil reservoir injection well of polymer flooding;4. the characteristic of analog parameter value is defined to obtaining
The bottom pressure changing value of the double-deck channelling oil reservoir injection well of polymer flooding arrived carries out nondimensionalization, and it is double-deck to draw polymer flooding
Channelling oil reservoir typical curve theory plate.
The step 2. in, the method for setting up the mathematical modeling of the double-deck channelling oil reservoir of polymer flooding comprises the following steps:a、
Set up the transient seepage flow differential equation of the single-phase micro- compressible liquid of each layer;B, the double-deck channelling oil reservoir of given polymer flooding mathematics
The primary condition equation of model;C, the mathematical modeling of the double-deck channelling oil reservoir of given polymer flooding boundary condition equation.
The transient seepage flow differential equation is:
In formula, p1、p2Respectively one or two stressor layers, unit is MPa;K1、K2Respectively one or two layers permeability, unit is 10-3
μm2;h1、h2Respectively one or two thickness degree, unit is m;Ct1、Ct2Respectively one or two layers system compressibility, unit is MPa-1;
φ1、φ2Respectively one or two layers porosity, dimensionless;R is radial distance, and unit is m;T is the time, and unit is h;A is up and down two
Laminar flow amount exchange coefficient, μpFor polymer solution initial viscosity.
The primary condition equation is:
p1(r, 0)=p2(r, 0)=poj,
p1(rw, t)=p2(rw, t)=pw(t),
In formula, rwFor well radius, pwFor flowing bottomhole pressure (FBHP), pojFor original formation pressure;
The boundary condition equation is:
Internal boundary condition:
Outer Boundary Conditions:
p1(∞, t)=p2(∞, t)=po,
In formula, q is flow, and unit is m3/ d, B are volume factor, dimensionless;C is bottom-hole storage coefficient, and unit is m3/
MPa;pwfFor bottom pressure, unit is MPa;S1、S2Respectively one or two layers skin factor, dimensionless;poFor original formation pressure,
Unit is MPa;μ2For the viscosity of two layers of polymer thing solution.
The step 3. in, the mathematical modeling of foundation is solved using finite difference method, polymer flooding is obtained double
The specific method of bottom pressure changing value of layer channelling oil reservoir injection well is:A, the spatial domain to mathematical modeling and time-domain are carried out
Mesh generation;B, by the basis partial differential equation of foundation, boundary condition equation and primary condition equation carry out difference discretization, obtain
Each equation of discrete scheme;C, obtained each discrete scheme equation is subjected to simultaneous, obtains the double-deck channelling oil reservoir injection of polymer flooding
The pressure formula of poly- well:
The step 4. in, the characteristic of definition includes dimensionless pressure, nondimensional time, dimensionless bottom-hole storage system
Ratio and interporosity flow coefficient are held in number, formation capacity ratio, storage.
The expression of each characteristic is:
Dimensionless pressure:
Nondimensional time:
Dimensionless bottom-hole storage coefficient:
Formation capacity ratio:
Chu Rong ratios:
Interporosity flow coefficient:
It is a kind of realize methods described be used for the double-deck channelling oil reservoir well test analysis system of polymer flooding, it includes collection singly
Member, memory cell, the first computing unit, the second computing unit, input block, judging unit and amending unit, it is characterised in that:
The bottom-hole pressure model of the double-deck channelling oil reservoir injection well of the polymer flooding that is stored with second computing unit, its according to by
The analog parameter value of the input block input calculates the simulation bottom pressure for obtaining the double-deck channelling oil reservoir injection well of polymer flooding
Variation relation value, and it is sent to memory cell and judging unit;The memory cell is used to store the mould inputted by input block
Intend parameter value simulation bottom pressure variation relation value corresponding with the analog parameter value;
Wherein, bottom-hole pressure model is:
Wherein, PwfFor bottom pressure, K is time step.
The system also includes control unit, and described control unit is used to be obtained according to revised analog parameter value and polymerize
Thing drives the well test analysis parameter value of double-deck channelling oil reservoir, and is driven according to well test analysis parameter value analyzing polymers after effect, takes
Corresponding measure instructs oil reservoir development.
The present invention is due to taking above technical scheme, and it has advantages below:1st, the present invention is characterized as with polymeric rheology
Basis, establishes the physical model of the double-deck channelling oil reservoir of polymer flooding, it is considered to the interaction between polymer and formation rock,
The problems such as including shearing, diffusion, convection current, inaccessible pore volume and permeability decrease coefficient, it can really reflect polymerization
The true seepage flow mechanism of thing solution in the earth formation.2nd, the present invention is due to using actual bottom pressure changing value to simulation bottom pressure
Variation relation value is modified, and obtains poly- according to the corresponding analog parameter value of revised simulation bottom pressure variation relation value
Compound drives double-deck channelling oil reservoir well test analysis parameter, improves the accuracy of well test analysis.3rd, the present invention is according to obtained well testing
Explain that parameter is analyzed polymer flooding effect, and corresponding measure is formulated according to analysis result, be used to refer to lead oil reservoir development, have
It is of practical significance.Thus the present invention can be widely applied in the well test analysis of the double-deck channelling oil reservoir of polymer flooding, for improving
The degree of accuracy of well test analysis.
Brief description of the drawings
Fig. 1 is that the present invention is used for the double-deck channelling oil reservoir well test analysis method flow chart of polymer flooding;
Fig. 2 is the drafting flow chart of the double-deck channelling oil reservoir typical curve plots of polymer flooding in the present invention;
Fig. 3 is the simplification physical model of the double-deck channelling oil reservoir of polymer flooding in the present invention;
Fig. 4 is the double-deck channelling oil reservoir well test analysis typical curve plots of polymer flooding in the present invention;
Fig. 5 is the typical curve plots under the influence of different interporosity flow coefficients in the present invention;
Fig. 6 is the typical curve plots under the influence of Different Strata coefficient ratio in the present invention;
Fig. 7 is the different typical curve plots stored up under the influence of appearance ratio in the present invention;
Fig. 8 is the typical curve plots under the influence of different initial polymer concentrations in the present invention;
Fig. 9 is structured flowchart of the present invention for the double-deck channelling oil reservoir well test analysis system of polymer flooding;
Figure 10 is that the present invention is used for the double-deck channelling oil reservoir well test analysis systematic difference schematic diagram of polymer flooding;
Test data and theoretical plate matched curve in embodiment in Figure 11 present invention.
Embodiment
The present invention is described in detail with reference to the accompanying drawings and examples.
As shown in figure 1, the present invention is used for the double-deck channelling oil reservoir well test analysis method of polymer flooding, comprise the following steps:
1) the bottom pressure value of the double-deck channelling oil reservoir injection well of collection polymer flooding, according to the corresponding collection of bottom pressure value
Time, the actual bottom pressure variation relation value for obtaining the double-deck channelling oil reservoir injection well of polymer flooding is calculated, and draw polymer
Drive the actual bottom pressure variation relation curve of double-deck channelling oil reservoir injection well.
Wherein, the actual bottom pressure variation relation of the double-deck channelling oil reservoir injection well of polymer flooding for actual bottom pressure and
The relation or actual bottom pressure of correspondence acquisition time and the derivative of corresponding acquisition time relation.For convenience of calculation, using immeasurable
Guiding principle amount represents, namely actual bottom pressure variation relation using dimensionless bottom pressure and nondimensional time double-log relation or
The double-log relation of person's dimensionless bottom pressure derivative and nondimensional time is represented.The actual bottom pressure used in the present invention becomes
Change relation is the double-log relation of dimensionless bottom pressure and nondimensional time, the bottom pressure value that will also gather and the well
The corresponding acquisition time of base pressure force value, is plotted in log-log coordinate (as shown in Figure 4).
2) according to the bottom pressure expression formula of the double-deck channelling oil reservoir injection well of polymer flooding and the analog parameter value of input, meter
The simulation bottom pressure variation relation value for obtaining the double-deck channelling oil reservoir injection well of polymer flooding is calculated, and draws polymer flooding bilayer and is altered
Flow the theoretical plate of typical curve of oil reservoir injection well.
As shown in Fig. 2 drawing the method for the theoretical plate of typical curve of the double-deck channelling oil reservoir injection well of polymer flooding includes
Following steps:
1. as shown in figure 3, according to the double-deck channelling characteristics of reservoirs of polymer flooding, setting up polymer solution viscosity model with polymerizeing
Thing drives the physical model of double-deck channelling oil reservoir, and the physical model meets it is assumed hereinafter that condition:
A, oil reservoir contain two kinds of components of water and polymer, and water and polymer complete miscibility, each oil reservoir polymer solution
Matter is identical, and each reservoir physical property is different;
B, the flowing of oil reservoir follow darcy flow, ignore the influence of gravity;
C, fluid are micro- compressible, and flow process is isothermal seepage flow, and interlayer is quasi-stable state channelling;
D, the storage of each layer well, epidermis are constant.
2. according to the polymer solution viscosity model of foundation and the double-deck channelling physical models of reservoir of polymer flooding, polymerization is set up
Thing drives the mathematical modeling of double-deck channelling oil reservoir.
The method for setting up the mathematical modeling of the double-deck channelling oil reservoir of polymer flooding comprises the following steps:
A, the transient seepage flow differential equation for setting up the single-phase micro- compressible liquid of each layer:
In formula (1) and (2), p1、p2Respectively one or two stressor layers, unit is MPa;K1、K2Respectively one or two layers permeability,
Unit is 10-3μm2;h1、h2Respectively one or two thickness degree, unit is m;Ct1、Ct2Respectively one or two layers system compressibility, unit
For MPa-1;φ1、φ2Respectively one or two layers porosity, dimensionless;R is radial distance, and unit is m;T is the time, and unit is h;a
For two layers of flow exchange coefficient up and down;μpFor polymer solution initial viscosity.
B, the mathematical modeling of the double-deck channelling oil reservoir of given polymer flooding primary condition equation:
p1(r, 0)=p2(r, 0)=poj (3)
p1(rw, t)=p2(rw, t)=pw(t) (4)
In formula (3) and (4), rwFor well radius, pwFor flowing bottomhole pressure (FBHP), pojFor original formation pressure.
C, the mathematical modeling of the double-deck channelling oil reservoir of given polymer flooding boundary condition equation.
Internal boundary condition:
Outer Boundary Conditions:
p1(∞, t)=p2(∞, t)=po (7)
Formula (5) is arrived in formula (7), and q is flow, and unit is m3/ d, B are volume factor, dimensionless;C is bottom-hole storage coefficient,
Unit is m3/MPa;pwfFor bottom pressure, unit is MPa;S1、S2Respectively one or two layers skin factor, dimensionless;poTo be original
Strata pressure, unit is MPa;μ2For the viscosity of two layers of polymer thing solution.
3. the mathematical modeling of foundation is solved using finite difference method, obtains the double-deck channelling oil reservoir injection of polymer flooding
The bottom pressure changing value of poly- well.
When being solved using finite difference method to the double-deck channelling reservoir mathematical model of polymer flooding, including following step
Suddenly:
A, the spatial domain to mathematical modeling and time-domain carry out mesh generation.
B, to step, 2. the basis partial differential equation of middle foundation, boundary condition equation and primary condition equation carry out difference discrete
Change, the discrete scheme for obtaining basis partial differential equation, boundary condition equation and primary condition equation is as follows:
(i) discrete scheme of basis partial differential equation
First layer oil reservoir seepage flow diffusion equation:
In formula (8), T is to write a Chinese character in simplified form, and
Second layer oil reservoir seepage flow diffusion equation:
In formula (10), M is writes a Chinese character in simplified form, and M is:
Formula (8) is arrived in formula (11),Respectively n+1 moment when one layer in i+1, i and i-1 net
The pressure of lattice;Respectively n+1 moment when two layers in i+1, i and i-1 grid pressure;Respectively n+1 moment when one or two layers of i-th of grid permeability;Respectively n+1
The permeability of one or two layers of i+1 grid during the moment;Respectively n+1 moment when one layer in i-th and i-1
The viscosity of grid;Respectively n+1 moment when the second layer in i-th and i-1 grid viscosity;Point
Not Wei n moment when one or two layers in i-th of grid pressure;riFor the radial distance of i-th of grid;Δ x is space lattice size,
Unit is m;Δ t is time step, and unit is h.
(ii) discrete scheme of boundary condition equation and primary condition equation
First layer inner boundary:
Second layer inner boundary:
External boundary:
Formula (12) is arrived in formula (14), and Q is flow, and unit is m2/d;S1、S2Respectively one or two layers skin factor, dimensionless;One layer of the 1st grid pressure during for the n+1 moment,One layer of initial pressure during for the n+1 moment,During for n+1
Two layers of the 1st grid pressure during quarter,One layer of initial pressure during for the n moment,One layer of the 1st grid pressure during for the n moment,Two layers of initial pressure during for the n moment,Two layers of the 1st grid pressure, p during for the n momentiFor original formation pressure,For n
Moment external boundary pressure,One layer of original permeability during for the n+1 moment,Two layers of initial infiltration during for the n+1 moment
Rate,Polymer initial viscosity, Δ t during for the n+1 momentnRespectively n moment time steps.
C, the discrete scheme progress simultaneous by obtained basis partial differential equation, boundary condition equation and primary condition equation,
Obtain the pressure distribution of the double-deck channelling oil reservoir injection well of polymer flooding.
Simultaneous formula (8)~(14),
By internal boundary condition, deformation is arranged:
It is that can obtain the bottom pressure change of the double-deck channelling oil reservoir injection well of polymer flooding according to formula (16):
Formula (15) is arrived in formula (17), PwfFor bottom pressure, unit is MPa;For n moment initial pressure,Respectively n
The stressor layer of moment one;S is skin factor;Δ x is space lattice size, and unit is m;Δ t is time step, and unit is h;N is
Space lattice number;K is time step.
4. as shown in figure 4, the characteristic for defining analog parameter value is gathered to the obtained double-deck channelling oil reservoir injection of polymer flooding
The bottom pressure changing value of well carries out nondimensionalization, and draws the theoretical plate of the double-deck channelling oil reservoir typical curve of polymer flooding.
The analog parameter value being related in the present invention include permeability, skin factor, bottom-hole storage coefficient, interporosity flow coefficient,
Ratio etc. is held in layer coefficients ratio and storage.The characteristic thus introduced includes dimensionless pressure, nondimensional time, dimensionless bottom-hole storage
Ratio and interporosity flow coefficient are held in coefficient, formation capacity ratio, storage, and the expression of each characteristic is:
Dimensionless pressure:
Nondimensional time:
Dimensionless bottom-hole storage coefficient:
Formation capacity ratio:
Chu Rong ratios:
Interporosity flow coefficient:
Above-mentioned each characteristic is substituted into formula (17), the well of the double-deck channelling oil reservoir injection well of dimensionless polymer flooding is obtained
Bottom pressure varying type.
As shown in figure 4, according to solve to carry out the bottom pressure of the double-deck channelling oil reservoir injection well of dimensionless polymer flooding and become
Change formula, the double-deck channelling oil reservoir well test analysis typical curve plots of polymer flooding of drafting.As can be seen from Figure, there is channelling
Dual-layer reservoir polymer flooding typical curve can be divided into 5 flowing stages:Ith section is pure bottom-hole storage stage, pressure and pressure
Power derivative overlap, reflection be well store up the stage pressure-responsive feature;IIth section is changeover portion, and the description pure bottom-hole storage stage arrives
The pressure-responsive feature in inner region radial flow stage;IIIth section is mid-term radial-flow section, and the flowing of expression layer 1 and layer 2 reaches radially
Flowing;IVth section is channelling acting section of the less permeable layer to high permeability zone;Vth section is that total system reaches radial flow active phase,
Due to being influenceed by polymer solution Newtonian fluid property, curve shows to upwarp by a small margin.
3) according to actual bottom pressure variation relation value, analog parameter value corresponding to simulation bottom pressure variation relation value
It is modified, and altered according to the corresponding analog parameter value acquisition polymer flooding bilayer of bottom pressure variation relation value is simulated after amendment
Flow the well test analysis parameter value of oil reservoir.
Specific steps include:
1. judge whether the difference for simulating bottom pressure variation relation value and actual bottom pressure variation relation value is less than one
Individual predetermined comparative factor;
If 2. both differences are more than predetermined comparative factor, analog parameter corresponding to simulation bottom pressure variation relation value
Value is modified, until both differences are less than predetermined comparative factor;
3. it is the double-deck channelling of polymer flooding to obtain the corresponding analog parameter value of simulation bottom pressure variation relation value after amendment
The well test analysis parameter value of oil reservoir.
When Fig. 5~Fig. 8 show change different analog parameter values, obtained typical curve plots.Wherein, Fig. 5 is difference
Typical curve plots under the influence of interporosity flow coefficient, Fig. 6 is the typical curve plots under the influence of Different Strata coefficient ratio, and Fig. 7 is not
Hold the typical curve plots under the influence of ratio with storage, Fig. 8 is the typical curve plots under the influence of different initial polymer concentrations.Can be with
Find out, in typical curve plots, the different WELL TEST INTERPRETATION MODEL of different analog parameter value correspondences, namely contain some well testing solutions
The mathematical equation (17) of parameter value is released, so when adjusting a certain simulative relation curve, the corresponding mould of simulative relation curve
Intending parameter value namely well test analysis parameter values can also change therewith, that is, correct the process of well test analysis parameter, pass through and intend
Close, obtain and the preferable simulative relation curve of actual bottom pressure variation relation curves, then the simulative relation curve pair
The well test analysis parameter answered can serve as the well test analysis parameter of actual bottom pressure variation relation curve, produces actual note and gathers
The well test analysis parameter of oil reservoir.
4) according to the well test analysis Parameter analysis polymer flooding reservoir situation of the double-deck channelling oil reservoir of obtained polymer flooding, adopt
Corresponding measure is taken to instruct oil field development.
If specifically, skin factor is excessive and permeability decrease is obvious, the basic de-plugging such as acidifying, pressure break can be taken to arrange
Apply;If acid fracturing effect is bad, sidetracking measure can further be taken according to reserves situation, polymer flooding effect is improved with this
Really, oil field development is instructed.
As shown in figure 9, the present invention, which also provides one kind, is used for the double-deck channelling oil reservoir well test analysis control system of polymer flooding, its
Including:Collecting unit, memory cell, the first computing unit, the second computing unit, input block, judging unit, amending unit and
Control unit.
Wherein, collecting unit is used to gather injection well in closing well bottom pressure value interior for a period of time;First computing unit
For according to the corresponding acquisition time of bottom pressure value, obtaining the actual bottom pressure of the double-deck channelling oil reservoir injection well of polymer flooding
With the relation of time;Input block is used to input analog parameter value, and is sent to memory cell and the second computing unit;Second meter
Calculating unit is used for the bottom hole pressure bomb according to the analog parameter value of input and the double-deck channelling oil reservoir injection well of the polymer flooding of storage
Model is calculated, the simulation bottom pressure variation relation value for obtaining the double-deck channelling oil reservoir injection well of polymer flooding is calculated, and is sent to and deposits
Storage unit and judging unit;Memory cell is used for the analog parameter value and the corresponding simulation well of the analog parameter value for storing input
Bottom pressure variation relation value;Judging unit is used to calculate actual bottom pressure variation relation value and simulation bottom pressure variation relation
The difference of value, and judge that whether both differences are less than a predetermined comparative factor, and will determine that result is sent to correcting process list
Member;Amending unit is used for the analog parameter value stored in the judged result and memory cell according to judging unit, to simulation well
The corresponding analog parameter value of bottom pressure variation relation value is modified;Control unit is used to be obtained according to revised analog parameter value
The well test analysis parameter value of the double-deck channelling oil reservoir of polymer flooding is taken, and effect is driven according to well test analysis parameter value analyzing polymers
Afterwards, corresponding measure is taken to instruct oil reservoir development.
The computation model of the double-deck channelling oil reservoir injection well of polymer flooding stored in second computing unit, including bottom pressure
Expression formula is:
Wherein, PwfFor bottom pressure, K is time step.
In above-described embodiment, collecting unit is adopted using conventional reservoir pressure collecting unit or measurement apparatus to realize
Collection, collection result can be supplied to the first computing unit, such as portable instrument transmission, network transmission, website hair in several ways
Cloth etc..
As shown in Figure 10, during the embodiment of the present invention is implemented, the first computing unit, memory cell, second calculate
Unit, judging unit, amending unit realize its respective function by computer server, and memory cell can also independently go out
To realize its function.
The present invention is described further with reference to specific embodiment.
The injection well that the present embodiment is used starts note on March 1st, 2013 and gathered, the polymer concentration 1500mg/L of injection,
Two layers of implanted layer position, thickness 10m;8 days 15 June in 2013:00, which stops note, surveys static pressure, pressure measurement force data 3d, June 11 in 2013
Day 15:00 recovers stake.There is the basic parameter of closing well and oil reservoir, it is as shown in table 1 below.
As shown in figure 11, it is well test data and the theoretical plate matched curve, fitting parameter prime stratum of the present embodiment
Pressure Pi=16.17MPa, the permeability K of layer 11=0.3 μm2, the permeability K of layer 22=0.00055 μm2, the epidermis S of layer 11=9.54, layer
2 epidermis S2=16.26, bottom-hole storage coefficient C=0.54m3/MPa。
The well of table 1 and oil reservoir basic parameter
Basic parameter | Symbol | Value |
Injection rate, m3/d | Q | 100 |
1 core intersection, m | h1 | 10 |
2 core intersections, m | h2 | 8 |
Oil volume factor | B | 1.2 |
Porosity, decimal | φ | 0.2 |
Viscosity/mPas of water | μw | 0.5 |
Well radius, m | rw | 0.1 |
1 layer of total compression coefficient, MPa-1 | Ct1 | 0.004 |
2 layers of total compression coefficient, MPa-1 | Ct2 | 0.036 |
Particular embodiments described above, has been carried out further in detail to the purpose of the present invention, technical scheme and beneficial effect
Describe in detail it is bright, the foregoing is only the present invention specific embodiment, and without with limit protection scope of the present invention, it is every
The equivalents carried out on the basis of technical solution of the present invention and improvement, should not be excluded outside protection scope of the present invention.
Claims (10)
1. one kind is used for the double-deck channelling oil reservoir well test analysis method of polymer flooding, it is characterised in that comprise the following steps:
1) the bottom pressure value of the double-deck channelling oil reservoir injection well of collection polymer flooding, during collection corresponding according to bottom pressure value
Between, the actual bottom pressure variation relation value for obtaining the double-deck channelling oil reservoir injection well of polymer flooding is calculated, and draw polymer flooding
The actual bottom pressure variation relation curve of double-deck channelling oil reservoir injection well;
2) according to the bottom pressure expression formula of the double-deck channelling oil reservoir injection well of polymer flooding and the analog parameter value of input, calculate
To the simulation bottom pressure variation relation value of the double-deck channelling oil reservoir injection well of polymer flooding, and draw the double-deck channelling oil of polymer flooding
Hide the theoretical plate of typical curve of injection well;
Wherein, bottom pressure expression formula is:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mo>...</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
In formula, PwfFor bottom pressure, unit is MPa;K is time step, tnFor the n times at moment;
3) according to actual bottom pressure variation relation value, analog parameter value corresponding to simulation bottom pressure variation relation value is carried out
Amendment, and obtain the double-deck channelling oil of polymer flooding according to the corresponding analog parameter value of bottom pressure variation relation value is simulated after amendment
The well test analysis parameter value of Tibetan;
4) reservoir situation is driven according to the well test analysis parameter value analyzing polymers of the double-deck channelling oil reservoir of obtained polymer flooding, taken
Corresponding measure instructs oil field development.
2. it is as claimed in claim 1 a kind of for the double-deck channelling oil reservoir well test analysis method of polymer flooding, it is characterised in that:Institute
State step 2) in, the method for drawing the theoretical plate of the double-deck channelling oil reservoir typical curve of polymer flooding comprises the following steps:
1. according to the double-deck channelling characteristics of reservoirs of polymer flooding, polymer solution viscosity model and the double-deck channelling oil of polymer flooding are set up
The physical model of Tibetan;
2. according to the polymer solution viscosity model of foundation and the double-deck channelling physical models of reservoir of polymer flooding, polymer flooding is set up
The mathematical modeling of double-deck channelling oil reservoir;
3. the mathematical modeling of foundation is solved using finite difference method, obtains the double-deck channelling oil reservoir injection well of polymer flooding
Bottom pressure changing value;
4. the bottom pressure for defining the double-deck channelling oil reservoir injection well of polymer flooding of the characteristic of analog parameter value to obtaining becomes
Change value carries out nondimensionalization, and draws the theoretical plate of the double-deck channelling oil reservoir typical curve of polymer flooding.
3. it is as claimed in claim 2 a kind of for the double-deck channelling oil reservoir well test analysis method of polymer flooding, it is characterised in that:Institute
In stating step 2., the method for setting up the mathematical modeling of the double-deck channelling oil reservoir of polymer flooding comprises the following steps:
A, the transient seepage flow differential equation for setting up the single-phase micro- compressible liquid of each layer;
B, the mathematical modeling of the double-deck channelling oil reservoir of given polymer flooding primary condition equation;
C, the mathematical modeling of the double-deck channelling oil reservoir of given polymer flooding boundary condition equation.
4. it is as claimed in claim 3 a kind of for the double-deck channelling oil reservoir well test analysis method of polymer flooding, it is characterised in that:Institute
Stating the transient seepage flow differential equation is:
<mrow>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>p</mi>
</msub>
</mfrac>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>a</mi>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>p</mi>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&phi;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>t</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>p</mi>
</msub>
</mfrac>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>a</mi>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>p</mi>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&phi;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
In formula, p1、p2Respectively one or two stressor layers, unit is MPa;K1、K2Respectively one or two layers permeability, unit is 10-3μm2;
h1、h2Respectively one or two thickness degree, unit is m;Ct1、Ct2Respectively one or two layers system compressibility, unit is MPa-1;φ1、
φ2Respectively one or two layers porosity, dimensionless;R is radial distance, and unit is m;T is the time, and unit is h;A is two layers up and down
Flow exchange coefficient, μpFor polymer solution initial viscosity.
5. it is as claimed in claim 3 a kind of for the double-deck channelling oil reservoir well test analysis method of polymer flooding, it is characterised in that:Institute
Stating primary condition equation is:
p1(r, 0)=p2(r, 0)=poj,
p1(rw, t)=p2(rw, t)=pw(t),
In formula, rwFor well radius, pwFor flowing bottomhole pressure (FBHP), pojFor original formation pressure;
The boundary condition equation is:
Internal boundary condition:
<mrow>
<mi>q</mi>
<mi>B</mi>
<mo>=</mo>
<mi>C</mi>
<mfrac>
<mrow>
<msub>
<mi>dp</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mi>p</mi>
</msub>
</mfrac>
<mi>r</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mi>r</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mi>r</mi>
<mo>=</mo>
<msub>
<mi>r</mi>
<mi>w</mi>
</msub>
</mrow>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>p</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
<mi>r</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mi>r</mi>
<mo>=</mo>
<msub>
<mi>r</mi>
<mi>w</mi>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
<mi>r</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mi>r</mi>
<mo>=</mo>
<msub>
<mi>r</mi>
<mi>w</mi>
</msub>
</mrow>
</msub>
<mo>,</mo>
</mrow>
Outer Boundary Conditions:
p1(∞, t)=p2(∞, t)=po,
In formula, q is flow, and unit is m3/ d, B are volume factor, dimensionless;C is bottom-hole storage coefficient, and unit is m3/MPa;pwf
For bottom pressure, unit is MPa;S1、S2Respectively one or two layers skin factor, dimensionless;poFor original formation pressure, unit is
MPa;μ2For the viscosity of two layers of polymer thing solution.
6. it is as claimed in claim 2 a kind of for the double-deck channelling oil reservoir well test analysis method of polymer flooding, it is characterised in that:Institute
In stating step 3., the mathematical modeling of foundation is solved using finite difference method, the double-deck channelling oil reservoir of polymer flooding is obtained
The specific method of the bottom pressure changing value of injection well is:
A, the spatial domain to mathematical modeling and time-domain carry out mesh generation;
B, by the basis partial differential equation of foundation, boundary condition equation and primary condition equation carry out difference discretization, obtain discrete
Each equation of form;
C, obtained each discrete scheme equation is subjected to simultaneous, the pressure gauge for obtaining the double-deck channelling oil reservoir injection well of polymer flooding reaches
Formula:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mo>...</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
7. it is as claimed in claim 2 a kind of for the double-deck channelling oil reservoir well test analysis method of polymer flooding, it is characterised in that:Institute
In stating step 4., the characteristic of definition include dimensionless pressure, nondimensional time, dimensionless bottom-hole storage coefficient, series of strata
Ratio and interporosity flow coefficient are held in number ratio, storage.
8. it is as claimed in claim 7 a kind of for the double-deck channelling oil reservoir well test analysis method of polymer flooding, it is characterised in that:Respectively
The expression of the characteristic is:
Dimensionless pressure:
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>D</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mn>1.842</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<msub>
<mi>q&mu;</mi>
<mi>p</mi>
</msub>
<mi>B</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>o</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Nondimensional time:
<mrow>
<msub>
<mi>t</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3.6</mn>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&phi;C</mi>
<mi>t</mi>
</msub>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>j</mi>
</msub>
<msub>
<mi>&mu;</mi>
<mi>p</mi>
</msub>
<msubsup>
<mi>r</mi>
<mi>w</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mi>t</mi>
<mo>,</mo>
</mrow>
Dimensionless bottom-hole storage coefficient:
<mrow>
<msub>
<mi>C</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>C</mi>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&pi;r</mi>
<mi>w</mi>
<mn>2</mn>
</msubsup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&phi;C</mi>
<mi>t</mi>
</msub>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
Formation capacity ratio:
<mrow>
<mi>&chi;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
Chu Rong ratios:
<mrow>
<mi>&omega;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&phi;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>t</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>&phi;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>t</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&phi;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
Interporosity flow coefficient:
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>ar</mi>
<mi>w</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
9. a kind of realize if any one of claim 1-8 methods described is for the double-deck channelling oil reservoir well test analysis system of polymer flooding
System, it includes collecting unit, memory cell, the first computing unit, the second computing unit, input block, judging unit and amendment
Unit, it is characterised in that:The bottom pressure of the double-deck channelling oil reservoir injection well of the polymer flooding that is stored with second computing unit
Computation model, it obtains the double-deck channelling oil reservoir injection of polymer flooding according to the analog parameter value calculating inputted by the input block and gathered
The simulation bottom pressure variation relation value of well, and it is sent to memory cell and judging unit;The memory cell be used for store by
The analog parameter value simulation bottom pressure variation relation value corresponding with the analog parameter value of input block input;
Wherein, bottom-hole pressure model is:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>f</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mo>...</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein, PwfFor bottom pressure, K is time step.
10. it is as claimed in claim 9 a kind of for the double-deck channelling oil reservoir well test analysis system of polymer flooding, it is characterised in that:
The system also includes control unit, and described control unit is used to obtain polymer flooding bilayer according to revised analog parameter value
The well test analysis parameter value of channelling oil reservoir, and driven according to well test analysis parameter value analyzing polymers after effect, take corresponding measure
Instruct oil reservoir development.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710056182.8A CN107066679A (en) | 2017-03-09 | 2017-03-09 | One kind is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710056182.8A CN107066679A (en) | 2017-03-09 | 2017-03-09 | One kind is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107066679A true CN107066679A (en) | 2017-08-18 |
Family
ID=59598411
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710056182.8A Pending CN107066679A (en) | 2017-03-09 | 2017-03-09 | One kind is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107066679A (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526891A (en) * | 2017-08-24 | 2017-12-29 | 中海石油(中国)有限公司 | A kind of polymer flooding macropore oil reservoir well test analysis method |
CN109033519A (en) * | 2018-06-22 | 2018-12-18 | 中国石油天然气股份有限公司 | A kind of Well Test Data Analysis Method and device of abnormal high pressure carbonate rock volatile oil reservoir |
CN110175412A (en) * | 2019-05-28 | 2019-08-27 | 广东海洋大学 | Reservoir Data approximating method, Reservoir Data analysis method and device |
CN110532579A (en) * | 2018-05-25 | 2019-12-03 | 中国科学技术大学 | A kind of calculation method of parameters, device and equipment |
CN111764879A (en) * | 2020-07-08 | 2020-10-13 | 中国石油大学(华东) | Evaluation method for flow channeling of binary composite flooding agent |
CN111927413A (en) * | 2020-08-11 | 2020-11-13 | 中海石油(中国)有限公司 | Method for determining reasonable bottom hole pressure of constant-pressure injection of polymer injection well |
CN112343566A (en) * | 2019-08-08 | 2021-02-09 | 中国石油天然气股份有限公司 | Method and device for determining bottom hole flowing pressure of injection polymer well |
CN112347599A (en) * | 2019-08-07 | 2021-02-09 | 中国石油天然气股份有限公司 | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment |
CN112685971A (en) * | 2020-12-26 | 2021-04-20 | 合肥工业大学 | Parameter interpretation method and device for hydrophilic oil reservoir |
CN113704946A (en) * | 2020-05-09 | 2021-11-26 | 中国石油化工股份有限公司 | Method, device, electronic equipment and medium for identifying reservoir boundary of double-layer commingled production well |
CN116455946A (en) * | 2023-06-19 | 2023-07-18 | 安徽井上天华科技有限公司 | Cloud-based high-frequency wellhead pressure production data analysis method |
CN116562179A (en) * | 2023-03-16 | 2023-08-08 | 东北石油大学 | Semi-resolution method for dynamic characteristics of polymer flooding pressure by considering rheological characteristics |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090114009A1 (en) * | 2005-05-10 | 2009-05-07 | Schlumberger Technology Corporation | Method for analysis of pressure response in underground formations |
CN101446196A (en) * | 2008-04-14 | 2009-06-03 | 中国石油大学(北京) | Well test analysis method and device of treble medium oil pool branch horizontal well |
CN105335600A (en) * | 2014-08-08 | 2016-02-17 | 中国科学技术大学 | Method and system for obtaining shear-thinning characteristic of polymer solution in stratum |
CN105464644A (en) * | 2015-11-16 | 2016-04-06 | 中国海洋石油总公司 | Wellhead pressure monitoring method and control system |
-
2017
- 2017-03-09 CN CN201710056182.8A patent/CN107066679A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090114009A1 (en) * | 2005-05-10 | 2009-05-07 | Schlumberger Technology Corporation | Method for analysis of pressure response in underground formations |
CN101446196A (en) * | 2008-04-14 | 2009-06-03 | 中国石油大学(北京) | Well test analysis method and device of treble medium oil pool branch horizontal well |
CN105335600A (en) * | 2014-08-08 | 2016-02-17 | 中国科学技术大学 | Method and system for obtaining shear-thinning characteristic of polymer solution in stratum |
CN105464644A (en) * | 2015-11-16 | 2016-04-06 | 中国海洋石油总公司 | Wellhead pressure monitoring method and control system |
Non-Patent Citations (2)
Title |
---|
朱常玉等: "《聚合物驱双层窜流复合模型试井分析方法》", 《陕西科技大学学报》 * |
郭辉等: "《聚合物驱双层油藏试井分析方法》", 《科学技术与工程》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526891B (en) * | 2017-08-24 | 2020-11-10 | 中海石油(中国)有限公司 | Polymer flooding large-pore oil reservoir well testing analysis method |
CN107526891A (en) * | 2017-08-24 | 2017-12-29 | 中海石油(中国)有限公司 | A kind of polymer flooding macropore oil reservoir well test analysis method |
CN110532579A (en) * | 2018-05-25 | 2019-12-03 | 中国科学技术大学 | A kind of calculation method of parameters, device and equipment |
CN110532579B (en) * | 2018-05-25 | 2022-10-28 | 中国科学技术大学 | Parameter calculation method, device and equipment |
CN109033519A (en) * | 2018-06-22 | 2018-12-18 | 中国石油天然气股份有限公司 | A kind of Well Test Data Analysis Method and device of abnormal high pressure carbonate rock volatile oil reservoir |
CN110175412A (en) * | 2019-05-28 | 2019-08-27 | 广东海洋大学 | Reservoir Data approximating method, Reservoir Data analysis method and device |
CN110175412B (en) * | 2019-05-28 | 2023-03-28 | 广东海洋大学 | Oil reservoir data fitting method, oil reservoir data analysis method and device |
CN112347599A (en) * | 2019-08-07 | 2021-02-09 | 中国石油天然气股份有限公司 | Polymer flooding concentric double-pipe injection parameter calculation method and device and computer equipment |
CN112343566B (en) * | 2019-08-08 | 2022-12-02 | 中国石油天然气股份有限公司 | Method and device for determining bottom hole flowing pressure of injection polymer well |
CN112343566A (en) * | 2019-08-08 | 2021-02-09 | 中国石油天然气股份有限公司 | Method and device for determining bottom hole flowing pressure of injection polymer well |
CN113704946B (en) * | 2020-05-09 | 2022-09-30 | 中国石油化工股份有限公司 | Method, device, electronic equipment and medium for identifying reservoir boundary of double-layer commingled production well |
CN113704946A (en) * | 2020-05-09 | 2021-11-26 | 中国石油化工股份有限公司 | Method, device, electronic equipment and medium for identifying reservoir boundary of double-layer commingled production well |
CN111764879B (en) * | 2020-07-08 | 2021-06-15 | 中国石油大学(华东) | Evaluation method for flow channeling of binary composite flooding agent |
CN111764879A (en) * | 2020-07-08 | 2020-10-13 | 中国石油大学(华东) | Evaluation method for flow channeling of binary composite flooding agent |
CN111927413A (en) * | 2020-08-11 | 2020-11-13 | 中海石油(中国)有限公司 | Method for determining reasonable bottom hole pressure of constant-pressure injection of polymer injection well |
CN112685971A (en) * | 2020-12-26 | 2021-04-20 | 合肥工业大学 | Parameter interpretation method and device for hydrophilic oil reservoir |
CN112685971B (en) * | 2020-12-26 | 2022-12-02 | 合肥工业大学 | Parameter interpretation method and device for hydrophilic oil reservoir |
CN116562179A (en) * | 2023-03-16 | 2023-08-08 | 东北石油大学 | Semi-resolution method for dynamic characteristics of polymer flooding pressure by considering rheological characteristics |
CN116562179B (en) * | 2023-03-16 | 2024-01-09 | 东北石油大学 | Semi-resolution method for dynamic characteristics of polymer flooding pressure by considering rheological characteristics |
CN116455946A (en) * | 2023-06-19 | 2023-07-18 | 安徽井上天华科技有限公司 | Cloud-based high-frequency wellhead pressure production data analysis method |
CN116455946B (en) * | 2023-06-19 | 2023-09-19 | 安徽井上天华科技有限公司 | Cloud-based high-frequency wellhead pressure production data analysis method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107066679A (en) | One kind is used for the double-deck channelling oil reservoir well test analysis system and method for polymer flooding | |
CN107301306B (en) | Dynamic non-resistance flow prediction method for tight sandstone gas reservoir fractured horizontal well | |
CN104101564B (en) | A kind of method of unstable state High Temperature High Pressure test flow in low permeability core free-boundary problem | |
CN104533370B (en) | Pressure break horizontal well oil reservoir, crack, pit shaft coupled model method | |
Behmanesh et al. | Impact of distance-of-investigation calculations on rate-transient analysis of unconventional gas and light-oil reservoirs: new formulations for linear flow | |
CN106761621A (en) | A kind of three layers of acquisition methods of channelling oil reservoir well test analysis parameter of polymer flooding | |
CN109594968A (en) | Fracture parameters evaluation method and system after a kind of shale gas multistage pressure break horizontal well pressure | |
CN106814018A (en) | A kind of tight rock gas phase relative permeability measurement apparatus and method | |
CN105840187A (en) | Method for calculating staged fracturing productivity of compact reservoir horizontal well | |
Ren et al. | The transient pressure behavior model of multiple horizontal wells with complex fracture networks in tight oil reservoir | |
CN106503407A (en) | There is the well test analysis method and device of the linear water enchroachment (invasion) oil reservoir of part connection tomography | |
CN104712328B (en) | The method of single flow unit producing status in Fast Evaluation Complex Reservoir | |
CN114427445A (en) | Method and system for calculating infinite stratum dynamic capacity of fracture-cavity oil reservoir | |
CN112016212B (en) | Reservoir longitudinal heterogeneity evaluation method based on seepage control equation | |
CN105443093A (en) | Combined well mouth testing device for polymer injection well and method thereof | |
CN104316448B (en) | A kind of Forecasting Methodology of high-order coal petrography gas phase permeability dynamic change | |
CN111177946B (en) | Simulation calculation method for total skin and compacted zone skin coefficients of high-pressure dense gas reservoir perforation | |
CN114066666A (en) | Method for analyzing connectivity among wells through injection-production profile monitoring data | |
CN105929461A (en) | Dynamic and static rock mechanical parameter correction system | |
CN105844011A (en) | Capillary model-based permeability computing method | |
CN115452681B (en) | Method for determining relative permeability of reservoir oil gas | |
Miguel | Transport phenomena through porous screens and openings: from theory to greenhouse practice | |
CN112412411A (en) | Method and device for establishing and analyzing well testing model of fracturing well multi-well system | |
CN115017841B (en) | Method and system for determining fracture-cavity spatial structure of broken solution combined production reservoir | |
CN112377184B (en) | Physical property parameter analysis method and device for beaded fracture-cavity type carbonate reservoir |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170818 |
|
RJ01 | Rejection of invention patent application after publication |