CN107526891B  Polymer flooding largepore oil reservoir well testing analysis method  Google Patents
Polymer flooding largepore 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
 238000004458 analytical method Methods 0.000 title claims abstract description 11
 230000035699 permeability Effects 0.000 claims abstract description 46
 230000015572 biosynthetic process Effects 0.000 claims description 24
 238000005755 formation reaction Methods 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
 238000007906 compression Methods 0.000 claims description 12
 XLYOFNOQVPJJNPUHFFFAOYSAN water Substances data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwLjAnIGhlaWdodD0nMzAwLjAnIHg9JzAuMCcgeT0nMC4wJz4gPC9yZWN0Pgo8dGV4dCB4PScxMDAuNScgeT0nMTcwLjAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5IPC90ZXh0Pgo8dGV4dCB4PScxMjYuMScgeT0nMTg2LjAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MjZweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID4yPC90ZXh0Pgo8dGV4dCB4PScxMzguMCcgeT0nMTcwLjAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NS4wJyBoZWlnaHQ9Jzg1LjAnIHg9JzAuMCcgeT0nMC4wJz4gPC9yZWN0Pgo8dGV4dCB4PScxMy4zJyB5PSc1My42JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjIzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+SDwvdGV4dD4KPHRleHQgeD0nMjguMicgeT0nNjIuOScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToxNXB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPjI8L3RleHQ+Cjx0ZXh0IHg9JzM1LjAnIHk9JzUzLjYnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MjNweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8L3N2Zz4K O XLYOFNOQVPJJNPUHFFFAOYSAN 0.000 claims description 11
 238000000034 method Methods 0.000 claims description 10
 210000003491 Skin Anatomy 0.000 claims description 6
 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
 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
 206010008531 Chills Diseases 0.000 description 1
 210000003800 Pharynx Anatomy 0.000 description 1
 238000004220 aggregation Methods 0.000 description 1
 230000002776 aggregation Effects 0.000 description 1
 230000001276 controlling effect Effects 0.000 description 1
 230000000875 corresponding Effects 0.000 description 1
 230000001808 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
 239000000700 tracer Substances 0.000 description 1
Classifications

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06F—ELECTRIC DIGITAL DATA PROCESSING
 G06F30/00—Computeraided design [CAD]
 G06F30/20—Design optimisation, verification or simulation

 E—FIXED CONSTRUCTIONS
 E21—EARTH DRILLING; MINING
 E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
 E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
 E21B43/16—Enhanced recovery methods for obtaining hydrocarbons

 E—FIXED CONSTRUCTIONS
 E21—EARTH DRILLING; MINING
 E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
 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 CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE 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
Abstract
The invention discloses a polymer flooding largepore 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 largepore oil reservoir physical model according to the characteristics of the polymer flooding largepore oil reservoir; 2) determining a mathematical model of the polymer flooding largepore oil reservoir according to the physical model of the polymer flooding largepore oil reservoir obtained in the step 1); 3) solving the mathematical model of the polymer flooding largepore 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 largepore 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 largepore oil reservoir well testing analysis method, and belongs to the technical field of well testing.
Background
In the process of longterm 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 extrahigh 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 lowefficiency 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 largepore oil reservoir and the determination of the distribution condition of the polymer flooding largepore 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 loglog curve after water injection. The Kjeldahl provides an injection well radial model considering mutual coupling of a large pore passage and a nonlarge 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 waterflooding reservoir, and the waterflooding reservoir is not researched.
Disclosure of Invention
Aiming at the problems, the invention aims to provide a polymer flooding largepore 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 largepore oil reservoir well testing analysis method is characterized by comprising the following steps:
1) establishing a polymer flooding largepore oil reservoir physical model according to the characteristics of the polymer flooding largepore oil reservoir;
2) determining a mathematical model of the polymer flooding largepore oil reservoir according to the physical model of the polymer flooding largepore oil reservoir obtained in the step 1);
3) solving the mathematical model of the polymer flooding largepore 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 largepore 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 x_{f}(ii) a Large pore permeability of K_{f}The oil layer permeability is K, the ratio of the large pore permeability to the oil layer permeability is beta, K_{f}Is far greater than K; the thickness of the large pore passage is the oil layer thickness H, and the equivalent width is W_{f}(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 largepore oil reservoir in the step 2) is as follows:
determining an unstable seepage differential equation of the singlephase slightly compressible liquid:
wherein P is the formation pressure; k is the oil layer permeability; mu.s_{a}Is the initial viscosity of the polymer; c_{m}Is an oil layerSynthesizing the compression coefficient; phi is a_{m}Is 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 largechannel oil reservoir:
the initial condition equation of the mathematical model of the polymer flooding largepore oil reservoir is as follows:
P(x,y,t)_{t＝0}＝P_{0} (2)
in the formula, P_{0}Is the original formation pressure;
the inner boundary condition equation of the mathematical model of the polymer flooding largepore 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; p_{wf}Is bottom hole flowing pressure; p_{w}Is the pressure at the borehole wall; r is_{w}Is the wellbore radius; r is_{e}Is the outer boundary radius; x is the number of_{0}、y_{0}Respectively 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 largepore 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 largepore 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., φ_{f}Is the porosity of the macropores, when phi is subscript m, i.e. + . phi_{m}Porosity of the oil layer;
wherein in formula (6)Is determined according to harmonic mean, namely:
solving equation (6) by using Gaussseidel iterative method, wherein the iterative equation is as follows:
in the formula, a_{i,j}、b_{i,j}、c_{i,j}、d_{i,j}、e_{i,j}And g_{i,j}Are all intermediate variables:
c_{i,j}＝d_{i,j}b_{i,j}e_{i,j}a_{i,j}g_{i,j}；
g_{i,j}＝277.78H_{i,j}Δx_{i,j}Δy_{i,j}φC_{t}/Δt；
q_{i,j}is the well flow rate at grid coordinates (i, j), well point grid q_{i,j}Q, nonwell point grid q_{i,j}＝0；P_{i,j}Is 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:
P_{1,j}＝P_{m,j}＝P_{e}(j＝1,2…k) (13)
P_{i,1}＝P_{i,k}＝P_{e}(i＝1,2…m) (14)
in the formula, P_{1,j}Is the formation pressure at grid coordinates (1, j); p_{m,j}Is the formation pressure at grid coordinates (m, j); p_{i,1}Is the formation pressure at grid coordinates (i, 1); p_{i,k}Is the formation pressure at grid coordinates (i, k); p_{e}Is 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 largepore oil reservoir.
The grid is rectangular or square, when the grid is rectangular,when a square grid, r_{e}＝0.208Δx。
Step 4) is to perform nondimensionalization on the change relation of the pressure along with the time, and the nondimensional equation is as follows:
P_{wD}the polymer flooding largechannel reservoir bottom hole pressure dimensionless value is obtained; c_{D}Dimensionless values for wellbore reservoir coefficients; x is the number of_{D}Is a nondimensionalized numerical value of the abscissa; y is_{D}Dimensionless numerical values for the total coordinates; t is t_{D}Is a time dimensionless numerical value; c_{f}The large pore comprehensive compression coefficient; c_{m}The comprehensive compression coefficient of the oil layer; phi is a_{f}Porosity in macropores; beta is the permeability ratio.
The fitting process of the polymer flooding largepore 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 largepore 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 largepore oil reservoir, obtaining a numerical value solution of bottom hole pressure of the polymer flooding largepore oil reservoir by using a finite element difference method, drawing a typical curve theoretical plate of the polymer flooding largepore oil reservoir, fitting the typical curve theoretical plate of the polymer flooding largepore 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 largepore 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 largepore oil reservoir well testing analysis method, which comprises the following steps:
1) establishing a physical model of the polymer flooding largepore oil reservoir according to the characteristics of the polymer flooding largepore 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 x_{f}(ii) a Large pore permeability of K_{f}The oil layer permeability is K, the ratio of the large pore permeability to the oil layer permeability is beta, K_{f}Is far greater than K; the thickness of the large pore passage is the oil layer thickness H, and the equivalent width is W_{f}(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 largepore oil reservoir according to the physical model of the polymer flooding largepore oil reservoir obtained in the step 1), wherein the specific process is as follows:
determining an unstable seepage differential equation of the singlephase slightly compressible liquid:
wherein P is the formation pressure; k is the oil layer permeability; mu.s_{a}Is the initial viscosity of the polymer; c_{m}The comprehensive compression coefficient of the oil layer; phi is a_{m}Is 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 largechannel oil reservoir:
the initial condition equation of the mathematical model of the polymer flooding largepore oil reservoir is as follows:
P(x,y,t)_{t＝0}＝P_{0} (2)
in the formula, P_{0}Is the original formation pressure;
the inner boundary condition equation of the mathematical model of the polymer flooding largepore 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; p_{wf}Is bottom hole flowing pressure; p_{w}Is the pressure at the borehole wall; r is_{w}Is the wellbore radius; r is_{e}Is the outer boundary radius; x is the number of_{0}、y_{0}Respectively 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 largepore oil reservoir is as follows:
3) solving the mathematical model of the polymer flooding largepore 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 largepore 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. + _{f}Is the porosity of the macropores, when phi is subscript m, i.e. + . phi_{m}Porosity of the oil layer;
whereinIs determined according to harmonic mean, namely:
solving equation (6) by using Gaussseidel iterative method, wherein the iterative equation is as follows:
in the formula, a_{i,j}、b_{i,j}、c_{i,j}、d_{i,j}、e_{i,j}And g_{i,j}Are all intermediate variables:
c_{i,j}＝d_{i,j}b_{i,j}e_{i,j}a_{i,j}g_{i,j}；
g_{i,j}＝277.78H_{i,j}Δx_{i,j}Δy_{i,j}φC_{t}/Δt；
q_{i,j}is the well flow rate at grid coordinates (i, j), well point grid q_{i,j}Q, nonwell point grid q_{i,j}＝0；P_{i,j}Is 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 gridsFor a square grid r_{e}＝0.208Δx。
Discretizing the external boundary condition to obtain:
P_{1,j}＝P_{m,j}＝P_{e}(j＝1,2…k) (13)
P_{i,1}＝P_{i,k}＝P_{e}(i＝1,2…m) (14)
in the formula, P_{1,j}Is the formation pressure at grid coordinates (1, j); p_{m,j}Is the formation pressure at grid coordinates (m, j); p_{i,1}Is the formation pressure at grid coordinates (i, 1); p_{i,k}Is the formation pressure at grid coordinates (i, k); p_{e}Is 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 largepore oil reservoir.
4) Carrying out nondimensionalization on the change relation of the pressure along with the time, and drawing a typical curve theoretical plate of the polymer flooding largepore oil reservoir, wherein the nondimensionalization equation is as follows:
P_{wD}the polymer flooding largechannel reservoir bottom hole pressure dimensionless value is obtained; c_{D}Dimensionless values for wellbore reservoir coefficients; x is the number of_{D}Is a nondimensionalized numerical value of the abscissa; y is_{D}Dimensionless numerical values for the total coordinates; t is t_{D}Is a time dimensionless numerical value; c_{f}The large pore comprehensive compression coefficient; c_{m}The comprehensive compression coefficient of the oil layer; phi is a_{f}Porosity in macropores; beta is the permeability ratio;
FIG. 3 is a schematic diagram showing a typical curve chart comparison of a largepore 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 largepore 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 profilecontrolling 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, K_{f}＝8850mD，W_{f}3.5 m. According to the result of the interpretation,fully profilecontrolling 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 largepore oil reservoir well testing analysis method is characterized by comprising the following steps:
1) establishing a polymer flooding largepore oil reservoir physical model according to the characteristics of the polymer flooding largepore oil reservoir;
2) determining a mathematical model of the polymer flooding largepore oil reservoir according to the physical model of the polymer flooding largepore oil reservoir obtained in the step 1);
the specific process of determining the mathematical model of the polymer flooding largepore oil reservoir in the step 2) is as follows:
determining an unstable seepage differential equation of the singlephase slightly compressible liquid:
wherein P is the formation pressure; k is the oil layer permeability; mu.s_{a}Is the initial viscosity of the polymer; c_{m}The comprehensive compression coefficient of the oil layer; phi is a_{m}Is 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 largechannel oil reservoir:
the initial condition equation of the mathematical model of the polymer flooding largepore oil reservoir is as follows:
P(x,y,t)_{t＝0}＝P_{0} (2)
in the formula, P_{0}Is the original formation pressure;
the inner boundary condition equation of the mathematical model of the polymer flooding largepore 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; p_{wf}Is bottom hole flowing pressure; p_{w}Is the pressure at the borehole wall; r is_{w}Is the wellbore radius; r is_{e}Is the outer boundary radius; x is the number of_{0}、y_{0}Respectively 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 largepore oil reservoir is as follows:
3) solving the mathematical model of the polymer flooding largepore 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 largepore 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 x_{f}(ii) a Large pore permeability of K_{f}The oil layer permeability is K, the ratio of the large pore permeability to the oil layer permeability is beta, K_{f}Greater than K; the thickness of the large pore passage is the oil layer thickness H, and the equivalent width is W_{f}(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 largepore 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., φ_{f}Is the porosity of the macropores, when phi is subscript m, i.e. + . phi_{m}Porosity of the oil layer;
wherein in formula (6)Is determined according to harmonic mean, namely:
solving equation (6) by using Gaussseidel iterative method, wherein the iterative equation is as follows:
in the formula, a_{i,j}、b_{i,j}、c_{i,j}、d_{i,j}、e_{i,j}And g_{i,j}Are all intermediate variables:
c_{i,j}＝d_{i,j}b_{i,j}e_{i,j}a_{i,j}g_{i,j}；
g_{i,j}＝277.78H_{i,j}Δx_{i,j}Δy_{i,j}φC_{t}/Δt；
q_{i,j}is the well flow rate at grid coordinates (i, j), well point grid q_{i,j}Q, nonwell point grid q_{i,j}＝0；P_{i,j}Is 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:
P_{1,j}＝P_{m,j}＝P_{e}(j＝1,2…k) (13)
P_{i,1}＝P_{i,k}＝P_{e}(i＝1,2…m) (14)
in the formula, P_{1,j}Is the formation pressure at grid coordinates (1, j); p_{m,j}Is the formation pressure at grid coordinates (m, j); p_{i,1}Is the formation pressure at grid coordinates (i, 1); p_{i,k}Is the formation pressure at grid coordinates (i, k); p_{e}Is 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 largepore 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, r_{e}＝0.208Δx。
4. The method for analyzing the well test of the polymer flooding largepore 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:
P_{wD}the polymer flooding largechannel reservoir bottom hole pressure dimensionless value is obtained; c_{D}Dimensionless values for wellbore reservoir coefficients; x is the number of_{D}Is a nondimensionalized numerical value of the abscissa; y is_{D}Dimensionless numerical values for the total coordinates; t is t_{D}Is a time dimensionless numerical value; c_{f}The large pore comprehensive compression coefficient; c_{m}Is an oil layerSynthesizing the compression coefficient; phi is a_{f}Porosity 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 largepore 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)  20170824  20170824  Polymer flooding largepore oil reservoir well testing analysis method 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201710735830.2A CN107526891B (en)  20170824  20170824  Polymer flooding largepore oil reservoir well testing analysis method 
Publications (2)
Publication Number  Publication Date 

CN107526891A CN107526891A (en)  20171229 
CN107526891B true CN107526891B (en)  20201110 
Family
ID=60682023
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201710735830.2A Active CN107526891B (en)  20170824  20170824  Polymer flooding largepore oil reservoir well testing analysis method 
Country Status (1)
Country  Link 

CN (1)  CN107526891B (en) 
Families Citing this family (4)
Publication number  Priority date  Publication date  Assignee  Title 

CN110486004B (en) *  20180514  20220510  中国石油天然气股份有限公司  Method and device for identifying water flow dominant channel of sandstone reservoir 
CN109033519A (en) *  20180622  20181218  中国石油天然气股份有限公司  A kind of Well Test Data Analysis Method and device of abnormal high pressure carbonate rock volatile oil reservoir 
CN109441415B (en) *  20181219  20210330  中国石油天然气股份有限公司  Well testing interpretation method of polymer flooding oil reservoir test well based on adjacent well interference 
CN110175412A (en) *  20190528  20190827  广东海洋大学  Reservoir Data approximating method, Reservoir Data analysis method and device 
Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN106503407A (en) *  20161205  20170315  中海石油(中国)有限公司  There is the well test analysis method and device of the linear water enchroachment (invasion) oil reservoir of part connection tomography 
CN106761621A (en) *  20170208  20170531  中海石油(中国)有限公司  A kind of three layers of acquisition methods of channelling oil reservoir well test analysis parameter of polymer flooding 
CN107066679A (en) *  20170309  20170818  中海石油(中国)有限公司  One kind is used for the doubledeck channelling oil reservoir well test analysis system and method for polymer flooding 

2017
 20170824 CN CN201710735830.2A patent/CN107526891B/en active Active
Patent Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN106503407A (en) *  20161205  20170315  中海石油(中国)有限公司  There is the well test analysis method and device of the linear water enchroachment (invasion) oil reservoir of part connection tomography 
CN106761621A (en) *  20170208  20170531  中海石油(中国)有限公司  A kind of three layers of acquisition methods of channelling oil reservoir well test analysis parameter of polymer flooding 
CN107066679A (en) *  20170309  20170818  中海石油(中国)有限公司  One kind is used for the doubledeck channelling oil reservoir well test analysis system and method for polymer flooding 
Also Published As
Publication number  Publication date 

CN107526891A (en)  20171229 
Similar Documents
Publication  Publication Date  Title 

CN107526891B (en)  Polymer flooding largepore oil reservoir well testing analysis method  
CN108518212B (en)  Method for calculating unsteady state yield of shale gas reservoir complex fracture network  
CN107301306B (en)  Dynamic nonresistance flow prediction method for tight sandstone gas reservoir fractured horizontal well  
CN102953726B (en)  Method and device for water drive oilfield advantage channel recognition  
CN106599449B (en)  Well test interpretation method for karst cave volume calculation  
CN105893679B (en)  Well Test Data Analysis Method is corrected in low yield horizontal well afterflow  
CN108133080B (en)  Heterogeneous fractured shale gas reservoir numerical simulation method considering unsteady adsorption  
CN107462936B (en)  Utilize the method for pressure monitoring Data Inversion low permeability reservoir nonDarcy percolation law  
CN105715241B (en)  A kind of measuring method of Curve of Polymer Displacement Relative Permeability  
CN110598167A (en)  Processing method of oilwater relative permeability experimental data of lowpermeability reservoir  
CN111581786B (en)  Well test interpretation method of well test interpretation model for analyzing fracturecave seriesmode doublehole composite reservoir  
CN107045671A (en)  Waterproducing gas well hydrops Risk Forecast Method  
CN106522928B (en)  Well testing method for unstable pressure drop of well logging head by stopping pump after acidizing and fracturing  
CN106761621B (en)  A kind of acquisition methods of three layers of channelling oil reservoir well test analysis parameter of polymer flooding  
CN109441415B (en)  Well testing interpretation method of polymer flooding oil reservoir test well based on adjacent well interference  
Feng et al.  Experimental and numerical study of gel particles movement and deposition in porous media after polymer flooding  
Zeng et al.  Experimental and theoretical characterization of the natural gas migration and accumulation mechanism in lowpermeability (tight) sandstone cores  
CN112069690A (en)  Evaluation method for test productivity of deep water fault block oil reservoir long horizontal well multistage oil nozzle  
Wang et al.  PreDarcy flow behavior of CO2 HuffnPuff development in Fuyu tight formation: Experiment and numerical evaluation  
Li et al.  Skin effect on singlewell pushpull tests with the presence of regional groundwater flow  
CN111734394A (en)  Method for determining unsteady flow bottom pressure of tight reservoir fracturing well  
CN111353205A (en)  Method for calculating stratum pressure and dynamic capacity of waterproducing gas well of tight gas reservoir  
CN110593865B (en)  Well testing interpretation method for characteristic parameters of oil reservoir fracture hole  
AU2012382975B2 (en)  Method and apparatus for formation tester data interpretation with diverse flow models  
CN112036097A (en)  Capacity calculation method for waterlock gas well 
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 