CN107526891B - Polymer flooding large-pore oil reservoir well testing analysis method - Google Patents
Polymer flooding large-pore oil reservoir well testing analysis method Download PDFInfo
- Publication number
- CN107526891B CN107526891B CN201710735830.2A CN201710735830A CN107526891B CN 107526891 B CN107526891 B CN 107526891B CN 201710735830 A CN201710735830 A CN 201710735830A CN 107526891 B CN107526891 B CN 107526891B
- Authority
- CN
- China
- Prior art keywords
- pore
- polymer flooding
- pressure
- oil reservoir
- reservoir
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 135
- 229920000642 polymer Polymers 0.000 title claims abstract description 110
- 238000012360 testing method Methods 0.000 title claims abstract description 31
- 238000004458 analytical method Methods 0.000 title claims abstract description 11
- 230000035699 permeability Effects 0.000 claims abstract description 46
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000013178 mathematical model Methods 0.000 claims abstract description 26
- 230000008859 change Effects 0.000 claims abstract description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 24
- 238000002347 injection Methods 0.000 claims description 18
- 239000007924 injection Substances 0.000 claims description 18
- 239000000243 solution Substances 0.000 claims description 13
- 230000006835 compression Effects 0.000 claims description 12
- 238000007906 compression Methods 0.000 claims description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 239000006185 dispersion Substances 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 239000003129 oil well Substances 0.000 claims description 3
- 239000010410 layer Substances 0.000 description 26
- 238000010586 diagram Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 6
- 239000003795 chemical substances by application Substances 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000011084 recovery Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 238000005259 measurement Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
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
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- 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)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention discloses a polymer flooding large-pore oil reservoir well testing analysis method, which is characterized by comprising the following steps of: it comprises the following steps: 1) establishing a polymer flooding large-pore oil reservoir physical model according to the characteristics of the polymer flooding large-pore oil reservoir; 2) determining a mathematical model of the polymer flooding large-pore oil reservoir according to the physical model of the polymer flooding large-pore oil reservoir obtained in the step 1); 3) solving the mathematical model of the polymer flooding large-pore oil reservoir obtained in the step 2) by adopting a finite difference method to obtain a numerical solution of bottom hole pressure; 4) carrying out dimensionless on the change relation of the pressure along with the time, and drawing a typical curve theoretical plate of the polymer flooding large-pore oil reservoir; 5) fitting the typical curve theoretical plate of the polymer flooding large pore oil reservoir obtained in the step 4) with an oil field measured data curve to obtain the equivalent width and permeability parameters of the large pore.
Description
Technical Field
The invention relates to a polymer flooding large-pore oil reservoir well testing analysis method, and belongs to the technical field of well testing.
Background
In the process of long-term water injection development, the offshore sandstone oil reservoir has large change in pore structure, and because the permeability of the reservoir is increased and the pore throat radius is increased, a high permeability strip and an extra-high permeability strip, namely a large pore channel, are easily formed in the reservoir. The existence of the large pore channel aggravates the interlayer contradiction, and the polymer injected by means of improving the recovery ratio is adopted in the later period, so that other parts in the reservoir layer are difficult to be affected along the low-efficiency or ineffective circulation in the large pore channel, the oil displacement efficiency is seriously influenced, the saturation difference of residual oil on the plane is obvious, and the polymer flooding effect is obviously poor. Therefore, the identification of the polymer flooding large-pore oil reservoir and the determination of the distribution condition of the polymer flooding large-pore oil reservoir have important significance for plugging the large pore, improving the polymer flooding effect and increasing the recovery ratio by adopting corresponding measures.
At present, the main means for identifying the large pore path are well testing data identification, production dynamic monitoring, tracer monitoring and well logging data identification. The most sensitive development parameters in the process of identifying the large pore passage are yield and pressure, and related researches are carried out on judging the existence of the large pore passage by using pressure drop test data through a well testing means. Steady history and Shi Ying adopt a double logarithmic curve of wellhead pressure drop of a water injection well to diagnose the existence of a large pore passage, and establish a theoretical model of water drive reservoir well testing with the large pore passage. Yangsheron utilizes a pressure drop well testing method to reflect the change condition of the permeability of the reservoir through the characteristics of a pressure drop log-log curve after water injection. The Kjeldahl provides an injection well radial model considering mutual coupling of a large pore passage and a non-large pore passage region, and the development area and the development multiple of the large pore passage are judged according to the bottom hole pressure change characteristic of the injection well. And the Liuhua flood establishes a dominant seepage channel well testing model by analyzing the basic form of the dominant seepage channel. The Li Chengyonguan establishes an asymmetric dominant seepage channel well testing interpretation mathematical model on the basis of history and rigor. However, the method established by the method is only suitable for identifying the large pore path of the water-flooding reservoir, and the water-flooding reservoir is not researched.
Disclosure of Invention
Aiming at the problems, the invention aims to provide a polymer flooding large-pore oil reservoir well testing analysis method, which can accurately explain measured data of an oil field to obtain related parameters of a large pore, provides data support for plugging the large pore and has high application value in a mine field.
In order to achieve the purpose, the invention adopts the following technical scheme: a polymer flooding large-pore oil reservoir well testing analysis method is characterized by comprising the following steps:
1) establishing a polymer flooding large-pore oil reservoir physical model according to the characteristics of the polymer flooding large-pore oil reservoir;
2) determining a mathematical model of the polymer flooding large-pore oil reservoir according to the physical model of the polymer flooding large-pore oil reservoir obtained in the step 1);
3) solving the mathematical model of the polymer flooding large-pore oil reservoir obtained in the step 2) by adopting a finite difference method to obtain a numerical solution of bottom hole pressure;
4) carrying out dimensionless on the change relation of the pressure along with the time, and drawing a typical curve theoretical plate of the polymer flooding large-pore oil reservoir;
5) fitting the typical curve theoretical plate of the polymer flooding large pore oil reservoir obtained in the step 4) with an oil field measured data curve to obtain the equivalent width and permeability parameters of the large pore.
The physical model established in the step 1) takes a shaft as a symmetric center, and two large channels are symmetrically distributed on two sides of the shaft; the polymer flooding large pore reservoir characteristics in the step 1) comprise: the reservoir is horizontal, equal thickness, homogeneous and isotropic; the equivalent extended length of the large pore channel is xf(ii) a Large pore permeability of KfThe oil layer permeability is K, the ratio of the large pore permeability to the oil layer permeability is beta, KfIs far greater than K; the thickness of the large pore passage is the oil layer thickness H, and the equivalent width is Wf(ii) a There is fluid exchange along the large channels and there is a pressure drop.
The specific process of determining the mathematical model of the polymer flooding large-pore oil reservoir in the step 2) is as follows:
determining an unstable seepage differential equation of the single-phase slightly compressible liquid:
wherein P is the formation pressure; k is the oil layer permeability; mu.saIs the initial viscosity of the polymer; cmIs an oil layerSynthesizing the compression coefficient; phi is amIs the oil layer porosity; t is the fluid flow time; x is the distance from the horizontal coordinate direction to the center of the well; y is the distance from the longitudinal coordinate direction to the center of the well;
determining an initial condition equation, an inner boundary condition equation and an outer boundary condition equation of the mathematical model of the polymer flooding large-channel oil reservoir:
the initial condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
P(x,y,t)|t＝0＝P0 (2)
in the formula, P0Is the original formation pressure;
the inner boundary condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
wherein H is the thickness of the oil layer; q is the well flow rate; c is a wellbore storage coefficient; b is a volume coefficient; s is the epidermal coefficient of the reservoir; pwfIs bottom hole flowing pressure; pwIs the pressure at the borehole wall; r iswIs the wellbore radius; r iseIs the outer boundary radius; x is the number of0、y0Respectively are the horizontal and vertical coordinates of the center of the bottom of the oil well; delta x, delta y differential infinitesimal; e is a constant, e ═ 2.7182818;
the outer boundary condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
in the step 3), the specific process of solving the mathematical model of the polymer flooding large oil reservoir obtained in the step 2) by adopting a finite difference method is as follows:
carrying out grid division on the polymer flooding large-pore oil reservoir physical model in the step 1) in space and time;
performing differential discretization on an unstable seepage differential equation, an initial condition equation, an inner boundary condition equation and an outer boundary condition equation, namely performing differential discretization on the formulas (1) to (5);
wherein, the dispersed seepage diffusion equation is as follows:
in the formula, i and j are the dispersion of space; n is the dispersion over time; Δ t is the time step, φ is the porosity, when φ is subscript f, i.e., φfIs the porosity of the macropores, when phi is subscript m, i.e. + -. phimPorosity of the oil layer;
wherein in formula (6)
Is determined according to harmonic mean, namely:
solving equation (6) by using Gauss-seidel iterative method, wherein the iterative equation is as follows:
in the formula, ai,j、bi,j、ci,j、di,j、ei,jAnd gi,jAre all intermediate variables:
ci,j＝-di,j-bi,j-ei,j-ai,j-gi,j；
gi,j＝277.78Hi,jΔxi,jΔyi,jφCt/Δt；
qi,jis the well flow rate at grid coordinates (i, j), well point grid qi,jQ, non-well point grid qi,j＝0；Pi,jIs the formation pressure at grid coordinate (i, j);
for the inner boundary condition, the well is taken as the convergent source item of the grid where the well is located for processing, and the differential grid of the inner boundary condition is subjected to linearization processing due to the fact that the pressure gradient near the bottom of the well is large, so that the following steps are obtained:
in the formula: μ is the viscosity of the polymer at any one time;
discretizing the external boundary condition to obtain:
P1,j＝Pm,j＝Pe(j＝1,2…k) (13)
Pi,1＝Pi,k＝Pe(i＝1,2…m) (14)
in the formula, P1,jIs the formation pressure at grid coordinates (1, j); pm,jIs the formation pressure at grid coordinates (m, j); pi,1Is the formation pressure at grid coordinates (i, 1); pi,kIs the formation pressure at grid coordinates (i, k); peIs the outer boundary pressure; m represents the number of meshes in the i direction, and k represents the number of meshes in the j direction;
to facilitate the differential calculation solution, the processing equation for the pressure gradient is as follows:
and thirdly, carrying out numerical iteration solution on the differential equation set, wherein the differential equation set is an equation set consisting of an unstable seepage differential equation, initial conditions and boundary conditions, and solving the numerical solution of the bottom hole pressure of the polymer flooding large-pore oil reservoir.
The grid is rectangular or square, when the grid is rectangular,
when a square grid, re＝0.208Δx。
Step 4) is to perform non-dimensionalization on the change relation of the pressure along with the time, and the non-dimensional equation is as follows:
PwDthe polymer flooding large-channel reservoir bottom hole pressure dimensionless value is obtained; cDDimensionless values for wellbore reservoir coefficients; x is the number ofDIs a non-dimensionalized numerical value of the abscissa; y isDDimensionless numerical values for the total coordinates; t is tDIs a time dimensionless numerical value; cfThe large pore comprehensive compression coefficient; cmThe comprehensive compression coefficient of the oil layer; phi is afPorosity in macropores; beta is the permeability ratio.
The fitting process of the polymer flooding large-pore oil reservoir typical curve theoretical plate and the oil field measured data curve in the step 5) is as follows: firstly, inputting basic data into a theoretical plate program of a polymer flooding large-pore oil reservoir, wherein the basic data comprises the thickness of an oil layer, the porosity of the oil layer, the permeability of the oil layer, a comprehensive compression coefficient, a skin coefficient, a shaft reservoir coefficient, the injection amount of a polymer injection well, a volume coefficient, water phase viscosity, initial concentration of a polymer, original formation pressure, well diameter, the width of a large pore and the permeability of the large pore; then, calculating to obtain theoretical pressure and a theoretical pressure derivative curve by adjusting the permeability of an oil layer, a skin coefficient, a shaft reservoir coefficient, the width of a large pore passage and the permeability of the large pore passage; then, fitting the theoretical pressure curve and the real well pressure curve, and the theoretical pressure derivative curve and the real well pressure derivative curve by using the theoretical pressure curve, the theoretical pressure derivative curve and the measured data of the oil field pressure; and finally obtaining the width and permeability parameters of the large pore channel according to the fitting result.
Due to the adoption of the technical scheme, the invention has the following advantages: the method comprises the steps of establishing a physical model and a mathematical model of the polymer flooding large-pore oil reservoir, obtaining a numerical value solution of bottom hole pressure of the polymer flooding large-pore oil reservoir by using a finite element difference method, drawing a typical curve theoretical plate of the polymer flooding large-pore oil reservoir, fitting the typical curve theoretical plate of the polymer flooding large-pore oil reservoir with an oil field actual measurement data curve to obtain parameters such as large pore width, permeability and the like, and providing data support for plugging large pores by selecting a proper plugging agent, improving the polymer flooding effect and improving the recovery ratio.
Drawings
FIG. 1 is a schematic overall flow diagram of the present invention;
FIG. 2 is a schematic diagram of a simplified physical model of a reservoir in accordance with the present invention;
FIG. 3 is a schematic diagram of a typical curve chart comparison of a large-pore reservoir in water flooding and polymer flooding according to the present invention;
FIG. 4 is a schematic illustration of a graphical representation of an exemplary well test interpretation of a reservoir under the influence of different initial polymer concentrations in accordance with the present invention;
FIG. 5 is a schematic diagram of a typical plot plate for reservoir well testing interpretation under the influence of different permeability ratios in accordance with the present invention;
FIG. 6 is a schematic diagram of a typical plot plate for reservoir well testing interpretation under the influence of different large pore widths according to the present invention;
FIG. 7 is a schematic diagram of a theoretical plate for reservoir well testing interpretation and a fitted curve of real well test data in the invention.
Detailed Description
The invention is described in detail below with reference to the figures and examples.
As shown in FIG. 1, the invention provides a polymer flooding large-pore oil reservoir well testing analysis method, which comprises the following steps:
1) establishing a physical model of the polymer flooding large-pore oil reservoir according to the characteristics of the polymer flooding large-pore oil reservoir, wherein the model takes a shaft as a symmetric center, and two large pores are symmetrically distributed on two sides of the shaft (as shown in figure 2);
wherein, the polymer flooding large pore oil reservoir characteristics comprise: the reservoir is horizontal, equal thickness, homogeneous and isotropic; the 2 large channels are symmetrical to the shaft, and the equivalent extension length of the large channels is xf(ii) a Large pore permeability of KfThe oil layer permeability is K, the ratio of the large pore permeability to the oil layer permeability is beta, KfIs far greater than K; the thickness of the large pore passage is the oil layer thickness H, and the equivalent width is Wf(ii) a There is fluid exchange along the large channels and there is a pressure drop.
2) Determining a mathematical model of the polymer flooding large-pore oil reservoir according to the physical model of the polymer flooding large-pore oil reservoir obtained in the step 1), wherein the specific process is as follows:
determining an unstable seepage differential equation of the single-phase slightly compressible liquid:
wherein P is the formation pressure; k is the oil layer permeability; mu.saIs the initial viscosity of the polymer; cmThe comprehensive compression coefficient of the oil layer; phi is amIs the oil layer porosity; t is the fluid flow time; x is the distance from the horizontal coordinate direction to the center of the well; y is the distance from the center of the well in the ordinate direction.
Determining an initial condition equation, an inner boundary condition equation and an outer boundary condition equation of the mathematical model of the polymer flooding large-channel oil reservoir:
the initial condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
P(x,y,t)|t＝0＝P0 (2)
in the formula, P0Is the original formation pressure;
the inner boundary condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
wherein H is the thickness of the oil layer; q is the well flow rate; c is a wellbore storage coefficient; b is a volume coefficient; s is the epidermal coefficient of the reservoir; pwfIs bottom hole flowing pressure; pwIs the pressure at the borehole wall; r iswIs the wellbore radius; r iseIs the outer boundary radius; x is the number of0、y0Respectively are the horizontal and vertical coordinates of the center of the bottom of the oil well; Δ x and Δ y are differential infinitesimals; e is a constant, e ═ 2.7182818;
the outer boundary condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
3) solving the mathematical model of the polymer flooding large-pore oil reservoir obtained in the step 2) by adopting a finite difference method to obtain a numerical solution of bottom hole pressure, wherein the finite difference method comprises the following specific processes:
carrying out grid division on the polymer flooding large-pore oil reservoir physical model in the step 1) in space and time;
secondly, carrying out differential discretization on an unstable seepage differential equation, an initial condition equation, an inner boundary condition equation and an outer boundary condition equation, namely carrying out differential discretization on equations (1) to (5):
wherein, the dispersed seepage diffusion equation is as follows:
in the formula, i and j are the dispersion of space; n is the dispersion over time; Δ t is the time step; phi is the porosity, when phi has the subscript f, i.e. + -fIs the porosity of the macropores, when phi is subscript m, i.e. + -. phimPorosity of the oil layer;
wherein
Is determined according to harmonic mean, namely:
solving equation (6) by using Gauss-seidel iterative method, wherein the iterative equation is as follows:
in the formula, ai,j、bi,j、ci,j、di,j、ei,jAnd gi,jAre all intermediate variables:
ci,j＝-di,j-bi,j-ei,j-ai,j-gi,j；
gi,j＝277.78Hi,jΔxi,jΔyi,jφCt/Δt；
qi,jis the well flow rate at grid coordinates (i, j), well point grid qi,jQ, non-well point grid qi,j＝0；Pi,jIs the formation pressure at grid coordinate (i, j).
For the inner boundary condition, the well is taken as the convergent source item of the grid where the well is located for processing, and the differential grid of the inner boundary condition is subjected to linearization processing due to the fact that the pressure gradient near the bottom of the well is large, so that the following steps are obtained:
in the formula: μ is the viscosity of the polymer at any one time; for rectangular grids
For a square grid re＝0.208Δx。
Discretizing the external boundary condition to obtain:
P1,j＝Pm,j＝Pe(j＝1,2…k) (13)
Pi,1＝Pi,k＝Pe(i＝1,2…m) (14)
in the formula, P1,jIs the formation pressure at grid coordinates (1, j); pm,jIs the formation pressure at grid coordinates (m, j); pi,1Is the formation pressure at grid coordinates (i, 1); pi,kIs the formation pressure at grid coordinates (i, k); peIs the outer boundary pressure; m represents the number of meshes in the i direction; k represents the number of meshes in the j direction;
to facilitate the differential calculation solution, the processing equation for the pressure gradient is as follows:
and thirdly, carrying out numerical iteration solution on the differential equation set, wherein the differential equation set is an equation set consisting of an unstable seepage differential equation, initial conditions and boundary conditions, and solving the numerical solution of the bottom hole pressure of the polymer flooding large-pore oil reservoir.
4) Carrying out non-dimensionalization on the change relation of the pressure along with the time, and drawing a typical curve theoretical plate of the polymer flooding large-pore oil reservoir, wherein the non-dimensionalization equation is as follows:
PwDthe polymer flooding large-channel reservoir bottom hole pressure dimensionless value is obtained; cDDimensionless values for wellbore reservoir coefficients; x is the number ofDIs a non-dimensionalized numerical value of the abscissa; y isDDimensionless numerical values for the total coordinates; t is tDIs a time dimensionless numerical value; cfThe large pore comprehensive compression coefficient; cmThe comprehensive compression coefficient of the oil layer; phi is afPorosity in macropores; beta is the permeability ratio;
FIG. 3 is a schematic diagram showing a typical curve chart comparison of a large-pore reservoir in water flooding and polymer flooding; FIG. 4 is a schematic diagram illustrating a typical plot plate for reservoir well testing under the influence of different initial polymer concentrations; as shown in fig. 5, a schematic diagram of a typical curve chart is explained for reservoir well testing under the influence of different permeability ratios; as shown in fig. 6, a schematic diagram of a typical plot plate is explained for reservoir well testing under the influence of different large pore widths.
5) Fitting the typical curve theoretical plate of the polymer flooding large-pore oil reservoir obtained in the step 4) with an oil field actual measurement data curve to obtain the equivalent width and permeability parameters (namely the permeability of the large pore) of the large pore, fully profile-controlling the stratum, and selecting a plugging agent with a proper particle size to plug the large pore, thereby improving the polymer flooding effect;
in step 5), basic data are input into a theoretical plate program of the polymer flooding large pore channel oil reservoir, wherein the basic data comprise the thickness of an oil layer, the porosity of the oil layer, the permeability of the oil layer, a comprehensive compression coefficient, a skin coefficient, a wellbore storage coefficient, the injection amount of a polymer injection well, a volume coefficient, water phase viscosity, initial concentration of a polymer, original formation pressure, a well diameter, the width of a large pore channel and the permeability of the large pore channel, and then theoretical pressure and a theoretical pressure derivative curve are calculated by adjusting the permeability of the oil layer, the skin coefficient, the wellbore storage coefficient, the width of the large pore channel and the permeability of the large pore channel; then, fitting the theoretical pressure curve and the real well pressure curve, and the theoretical pressure derivative curve and the real well pressure derivative curve by using the theoretical pressure curve, the theoretical pressure derivative curve and the measured data of the oil field pressure; and finally, obtaining the width and permeability parameters of the large pore passage according to the fitting result, fully adjusting and profiling the stratum by utilizing the explanation result, and selecting a plugging agent with a proper particle size to plug the large pore passage, thereby improving the polymer flooding effect and improving the recovery ratio.
Specific examples are listed below:
examples
Example data is obtained from pressure drop data of a certain polymer injection well in a Bohai sea B oil field, the oil reservoir stratum is strong in heterogeneity, the porosity is 0.31, the average permeability is 2000mD, a test well group carries out transfer polymer injection from 2013 to 8 months, the concentration of injected polymers is 1500mg/L, since polymer injection, a production well has high polymer injection speed, the polymer production concentration is high, the water content rises quickly, the polymer injection effect is reduced, and the fact that a large pore channel exists in the well group is preliminarily judged. The injection well was tested for pressure drop in 2016, 3, and 2 days.
The polymer flooding large pore channel oil reservoir well testing analysis method is adopted to perform well testing explanation on a certain polymer injection well in the Bohai sea B oil field, the curve (shown in figure 7) is fitted by the actually measured data of the oil field and a theoretical plate in the embodiment, and the model explanation obtains the large pore channel parameters as follows: beta 4.4, Kf＝8850mD，Wf3.5 m. According to the result of the interpretation,fully profile-controlling the stratum, selecting a plugging agent with proper particle size to plug a large pore, effectively controlling the water content of the well group at present, and reducing the aggregation concentration.
The above embodiments are only used for illustrating the present invention, and the structure, connection mode and the like of each component can be changed, and all equivalent changes and improvements made on the basis of the technical scheme of the present invention should not be excluded from the protection scope of the present invention.
Claims (5)
1. A polymer flooding large-pore oil reservoir well testing analysis method is characterized by comprising the following steps:
1) establishing a polymer flooding large-pore oil reservoir physical model according to the characteristics of the polymer flooding large-pore oil reservoir;
2) determining a mathematical model of the polymer flooding large-pore oil reservoir according to the physical model of the polymer flooding large-pore oil reservoir obtained in the step 1);
the specific process of determining the mathematical model of the polymer flooding large-pore oil reservoir in the step 2) is as follows:
determining an unstable seepage differential equation of the single-phase slightly compressible liquid:
wherein P is the formation pressure; k is the oil layer permeability; mu.saIs the initial viscosity of the polymer; cmThe comprehensive compression coefficient of the oil layer; phi is amIs the oil layer porosity; t is the fluid flow time; x is the distance from the horizontal coordinate direction to the center of the well; y is the distance from the longitudinal coordinate direction to the center of the well;
determining an initial condition equation, an inner boundary condition equation and an outer boundary condition equation of the mathematical model of the polymer flooding large-channel oil reservoir:
the initial condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
P(x,y,t)|t＝0＝P0 (2)
in the formula, P0Is the original formation pressure;
the inner boundary condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
wherein H is the thickness of the oil layer; q is the well flow rate; c is a wellbore storage coefficient; b is a volume coefficient; s is the epidermal coefficient of the reservoir; pwfIs bottom hole flowing pressure; pwIs the pressure at the borehole wall; r iswIs the wellbore radius; r iseIs the outer boundary radius; x is the number of0、y0Respectively are the horizontal and vertical coordinates of the center of the bottom of the oil well; delta x, delta y differential infinitesimal; e is a constant, e ═ 2.7182818;
the outer boundary condition equation of the mathematical model of the polymer flooding large-pore oil reservoir is as follows:
3) solving the mathematical model of the polymer flooding large-pore oil reservoir obtained in the step 2) by adopting a finite difference method to obtain a numerical solution of bottom hole pressure;
4) carrying out dimensionless on the change relation of the pressure along with the time, and drawing a typical curve theoretical plate of the polymer flooding large-pore oil reservoir;
5) fitting the typical curve theoretical plate of the polymer flooding large pore oil reservoir obtained in the step 4) with an oil field measured data curve to obtain the equivalent width and permeability parameters of the large pore;
the physical model established in the step 1) takes a shaft as a symmetric center, and two large channels are symmetrically distributed on two sides of the shaft; the polymer flooding large pore reservoir characteristics in the step 1) comprise: the reservoir is horizontal, equal thickness, homogeneous and isotropic; of the large pore channelEffective extension length of xf(ii) a Large pore permeability of KfThe oil layer permeability is K, the ratio of the large pore permeability to the oil layer permeability is beta, KfGreater than K; the thickness of the large pore passage is the oil layer thickness H, and the equivalent width is Wf(ii) a There is fluid exchange along the large channels and there is a pressure drop.
2. The polymer flooding large pore reservoir well testing analysis method of claim 1, characterized by: in the step 3), the specific process of solving the mathematical model of the polymer flooding large oil reservoir obtained in the step 2) by adopting a finite difference method is as follows:
carrying out grid division on the polymer flooding large-pore oil reservoir physical model in the step 1) in space and time;
performing differential discretization on an unstable seepage differential equation, an initial condition equation, an inner boundary condition equation and an outer boundary condition equation, namely performing differential discretization on the formulas (1) to (5);
wherein, the dispersed seepage diffusion equation is as follows:
in the formula, i and j are the dispersion of space; n is the dispersion over time; Δ t is the time step, φ is the porosity, when φ is subscript f, i.e., φfIs the porosity of the macropores, when phi is subscript m, i.e. + -. phimPorosity of the oil layer;
wherein in formula (6)
Is determined according to harmonic mean, namely:
solving equation (6) by using Gauss-seidel iterative method, wherein the iterative equation is as follows:
in the formula, ai,j、bi,j、ci,j、di,j、ei,jAnd gi,jAre all intermediate variables:
ci,j＝-di,j-bi,j-ei,j-ai,j-gi,j；
gi,j＝277.78Hi,jΔxi,jΔyi,jφCt/Δt；
qi,jis the well flow rate at grid coordinates (i, j), well point grid qi,jQ, non-well point grid qi,j＝0；Pi,jIs the formation pressure at grid coordinate (i, j);
for the inner boundary condition, the well is taken as the convergent source item of the grid where the well is located for processing, and because the pressure gradient near the bottom of the well is large, the difference grid of the inner boundary condition is subjected to linearization processing, so that the following steps are obtained:
in the formula: μ is the viscosity of the polymer at any one time;
discretizing the external boundary condition to obtain:
P1,j＝Pm,j＝Pe(j＝1,2…k) (13)
Pi,1＝Pi,k＝Pe(i＝1,2…m) (14)
in the formula, P1,jIs the formation pressure at grid coordinates (1, j); pm,jIs the formation pressure at grid coordinates (m, j); pi,1Is the formation pressure at grid coordinates (i, 1); pi,kIs the formation pressure at grid coordinates (i, k); peIs the outer boundary pressure; m represents the number of meshes in the i direction, and k represents the number of meshes in the j direction;
to facilitate the differential calculation solution, the processing equation for the pressure gradient is as follows:
and thirdly, carrying out numerical iteration solution on the differential equation set, wherein the differential equation set is an equation set consisting of an unstable seepage differential equation, initial conditions and boundary conditions, and solving the numerical solution of the bottom hole pressure of the polymer flooding large-pore oil reservoir.
3. The polymer flooding large pore reservoir well testing analysis method of claim 2, characterized in that: the grid is rectangular or square, when the grid is rectangular,
when a square grid, re＝0.208Δx。
4. The method for analyzing the well test of the polymer flooding large-pore oil reservoir according to claim 1, wherein the step 4) is a dimensionless method for the pressure change with time, and the dimensionless equation is as follows:
PwDthe polymer flooding large-channel reservoir bottom hole pressure dimensionless value is obtained; cDDimensionless values for wellbore reservoir coefficients; x is the number ofDIs a non-dimensionalized numerical value of the abscissa; y isDDimensionless numerical values for the total coordinates; t is tDIs a time dimensionless numerical value; cfThe large pore comprehensive compression coefficient; cmIs an oil layerSynthesizing the compression coefficient; phi is afPorosity of macropores; beta is the permeability ratio.
5. The method for analyzing the well testing of the polymer flooding large pore reservoir as claimed in claim 1, wherein the fitting process of the typical curve theoretical plate of the polymer flooding large pore reservoir and the measured data curve of the oil field in the step 5) is as follows: firstly, inputting basic data into a theoretical plate program of a polymer flooding large-pore oil reservoir, wherein the basic data comprises the thickness of an oil layer, the porosity of the oil layer, the permeability of the oil layer, a comprehensive compression coefficient, a skin coefficient, a shaft reservoir coefficient, the injection amount of a polymer injection well, a volume coefficient, water phase viscosity, initial concentration of a polymer, original formation pressure, well diameter, the width of a large pore and the permeability of the large pore; then, calculating to obtain theoretical pressure and a theoretical pressure derivative curve by adjusting the permeability of an oil layer, a skin coefficient, a shaft reservoir coefficient, the width of a large pore passage and the permeability of the large pore passage; then, fitting the theoretical pressure curve and the real well pressure curve, and the theoretical pressure derivative curve and the real well pressure derivative curve by using the theoretical pressure curve, the theoretical pressure derivative curve and the measured data of the oil field pressure; and finally obtaining the width and permeability parameters of the large pore channel according to the fitting result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710735830.2A CN107526891B (en) | 2017-08-24 | 2017-08-24 | Polymer flooding large-pore oil reservoir well testing analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710735830.2A CN107526891B (en) | 2017-08-24 | 2017-08-24 | Polymer flooding large-pore oil reservoir well testing analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107526891A CN107526891A (en) | 2017-12-29 |
| CN107526891B true CN107526891B (en) | 2020-11-10 |
Family
ID=60682023
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710735830.2A Active CN107526891B (en) | 2017-08-24 | 2017-08-24 | Polymer flooding large-pore oil reservoir well testing analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107526891B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110486004B (en) * | 2018-05-14 | 2022-05-10 | 中国石油天然气股份有限公司 | Method and device for identifying water flow dominant channel of sandstone reservoir |
| CN109033519A (en) * | 2018-06-22 | 2018-12-18 | 中国石油天然气股份有限公司 | Well testing interpretation method and device for abnormally high pressure carbonate rock volatile oil reservoir |
| CN109441415B (en) * | 2018-12-19 | 2021-03-30 | 中国石油天然气股份有限公司 | Well testing interpretation method of polymer flooding oil reservoir test well based on adjacent well interference |
| CN110175412B (en) * | 2019-05-28 | 2023-03-28 | 广东海洋大学 | Oil reservoir data fitting method, oil reservoir data analysis method and device |
| CN110630245B (en) * | 2019-09-20 | 2022-09-16 | 中国石油集团川庆钻探工程有限公司 | Single-well oil drainage area quantitative calculation method based on permeability advantage channel |
| CN112576243B (en) * | 2019-09-30 | 2024-05-31 | 中国石油化工股份有限公司 | Method for predicting pore blockage of oil and gas reservoir caused by injection of suspended matters in water |
| CN114060018B (en) * | 2020-08-04 | 2024-05-28 | 中国石油天然气股份有限公司 | Reservoir dynamic reserve determination method, system, equipment and readable storage medium |
| CN112395816A (en) * | 2020-11-16 | 2021-02-23 | 中国海洋石油集团有限公司 | Polymer flooding reservoir interference well testing analysis method |
| CN116562179B (en) * | 2023-03-16 | 2024-01-09 | 东北石油大学 | Semi-resolution method for dynamic characteristics of polymer flooding pressure by considering rheological characteristics |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106503407A (en) * | 2016-12-05 | 2017-03-15 | 中海石油(中国)有限公司 | There is the well test analysis method and device of the linear water enchroachment (invasion) oil reservoir of part connection tomography |
| CN106761621A (en) * | 2017-02-08 | 2017-05-31 | 中海石油(中国)有限公司 | A kind of three layers of acquisition methods of channelling oil reservoir well test analysis parameter of polymer flooding |
| 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 |
-
2017
- 2017-08-24 CN CN201710735830.2A patent/CN107526891B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106503407A (en) * | 2016-12-05 | 2017-03-15 | 中海石油(中国)有限公司 | There is the well test analysis method and device of the linear water enchroachment (invasion) oil reservoir of part connection tomography |
| CN106761621A (en) * | 2017-02-08 | 2017-05-31 | 中海石油(中国)有限公司 | A kind of three layers of acquisition methods of channelling oil reservoir well test analysis parameter of polymer flooding |
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107526891A (en) | 2017-12-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107526891B (en) | Polymer flooding large-pore oil reservoir well testing analysis method | |
| CN107301306B (en) | Dynamic non-resistance flow prediction method for tight sandstone gas reservoir fractured horizontal well | |
| CN106599449B (en) | Well test interpretation method for karst cave volume calculation | |
| CN111581786B (en) | Well test interpretation method of well test interpretation model for analyzing fracture-cave series-mode double-hole composite reservoir | |
| CN102953726B (en) | Method and device for water drive oilfield advantage channel recognition | |
| CN108518212A (en) | A method of calculating shale gas reservoir complex fracture network unstable state yield | |
| CN105893679B (en) | Well Test Data Analysis Method is corrected in low yield horizontal well afterflow | |
| CN110593865B (en) | Well testing interpretation method for characteristic parameters of oil reservoir fracture hole | |
| CN107462936B (en) | Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law | |
| CN110598167A (en) | Processing method of oil-water relative permeability experimental data of low-permeability reservoir | |
| CN107045671A (en) | Water-producing gas well hydrops Risk Forecast Method | |
| CN108133080B (en) | Heterogeneous fractured shale gas reservoir numerical simulation method considering unsteady adsorption | |
| CN111734394B (en) | Method for determining unsteady bottom-hole pressure of tight oil reservoir fracturing well | |
| CN107832574B (en) | Horizontal well water flooded layer interpretation method based on logging while drilling | |
| Feng et al. | Experimental and numerical study of gel particles movement and deposition in porous media after polymer flooding | |
| CN109441415B (en) | Well testing interpretation method of polymer flooding oil reservoir test well based on adjacent well interference | |
| CN111353205A (en) | Method for calculating stratum pressure and dynamic capacity of water-producing gas well of tight gas reservoir | |
| Zeng et al. | Experimental and theoretical characterization of the natural gas migration and accumulation mechanism in low-permeability (tight) sandstone cores | |
| CN110735633A (en) | early-stage identification method for low-permeability carbonate gas reservoir gas well shaft effusion | |
| Liu et al. | An inversion method of relative permeability curves in polymer flooding considering physical properties of polymer | |
| CN114372398B (en) | Polymer flooding well testing interpretation method considering reservoir fracture closure effect | |
| CN113294147B (en) | Single-hole type broken solution reservoir well testing interpretation method considering gravity factor influence | |
| CN112257349A (en) | Method for judging whether compact sandstone movable water-gas reservoir gas well has development value | |
| CN112069690A (en) | Evaluation method for test productivity of deep water fault block oil reservoir long horizontal well multistage oil nozzle | |
| Yu et al. | Insight into gas threshold pressure gradient and permeability of coal seam: Principle and method for field test |
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 |