CN106055843A - Quantitative analysis method and device of polymer production curve in polymer flooding process - Google Patents

Quantitative analysis method and device of polymer production curve in polymer flooding process Download PDF

Info

Publication number
CN106055843A
CN106055843A CN201610515517.3A CN201610515517A CN106055843A CN 106055843 A CN106055843 A CN 106055843A CN 201610515517 A CN201610515517 A CN 201610515517A CN 106055843 A CN106055843 A CN 106055843A
Authority
CN
China
Prior art keywords
poly
polymer
unit
polymer concentration
passage
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.)
Granted
Application number
CN201610515517.3A
Other languages
Chinese (zh)
Other versions
CN106055843B (en
Inventor
周建
刘金菊
王佳
江礼武
成杰
刘同敬
第五鹏祥
刘睿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Jinshi Source Technology Co Ltd
China University of Petroleum Beijing
Original Assignee
Beijing Jinshi Source Technology Co Ltd
China University of Petroleum Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Jinshi Source Technology Co Ltd, China University of Petroleum Beijing filed Critical Beijing Jinshi Source Technology Co Ltd
Priority to CN201610515517.3A priority Critical patent/CN106055843B/en
Publication of CN106055843A publication Critical patent/CN106055843A/en
Application granted granted Critical
Publication of CN106055843B publication Critical patent/CN106055843B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Abstract

The invention relates to a quantitative analysis method and device of a polymer production curve in a polymer flooding process. The method comprises the following steps: carrying out statistics on a relationship between the polymer concentrations of different oil well output liquid and monitoring time to obtain a polymer production curve, and classifying the polymer production curve according to curve characteristics; establishing a polymer concentration explanation model when a polymer gathering channel is not considered; utilizing the polymer concentration explanation model when the polymer gathering channel is not considered to establish a polymer concentration explanation model when the polymer gathering channel is considered; and detecting the polymer concentration of a producing well, fitting a polymer production curve for the detected polymer concentration according to the polymer concentration explanation model when the polymer gathering channel is considered, utilizing a classification result of the polymer production curve to classify the fit polymer production curve, and applying the polymer concentration explanation model when the polymer gathering channel is considered to analyze the fit polymer production curve.

Description

A kind of poly-drive during produce the quantitative analysis method and device of poly-curve
Technical field
The present invention relates to produce poly-curve interpretation technical field, particularly to a kind of poly-drive during produce the quantization of poly-curve and divide Analysis method and device.
Background technology
For waterflooding oil field, the primary heterogeneous body of oil reservoir and long-period water drive make anisotropism be further exacerbated by, in oil reservoir Gradually form high bleed-through road or macropore, form water breakthrough passage between oil-water well, have a strong impact on reservoir water drive development effectiveness.Therefore, Effective water and oil control is current main trend.
Polymer flooding is conventional tertiary oil recovery method, under proper condition, it is possible to significantly play decreasing water cut and increasing oil effect, But individual well takes effect in Polymer flooding block, seeing that poly-degree difference is relatively big, see that poly-rule is different, this is owing to reservoir property, local are oozed Rate, the oil many factors such as saturation and injection-production relation impact thoroughly.
Summary of the invention
For solving problem of the prior art, the present invention propose a kind of poly-drive during produce poly-curve quantitative analysis method and Device, to produce poly-curvilinear characteristic explanation difficulty during solving poly-driving, fluid-channeling channel and fluid parameter are difficult to only be produced by analysis Poly-feature just can be with the problem of quantitatively characterizing.
For achieving the above object, the invention provides a kind of poly-drive during produce the quantitative analysis method of poly-curve, including:
The different produced liquid in oil well polymer concentration of statistics and the relation of monitoring time, it is thus achieved that produce poly-curve, and according to curve Feature curve poly-to described product is classified;
Set up polymer concentration interpretation model when not considering to alter poly-passage;
The polymer utilizing polymer concentration interpretation model when not considering to alter poly-passage to set up when considering to alter poly-passage is dense Degree interpretation model;
The polymer concentration of detection producing well, according to the polymer concentration interpretation model considered when altering poly-passage to detecting Polymer concentration matching produce poly-curve, and utilize the classification results of the poly-curve of the described product poly-curve of product to simulating to carry out Classification, the poly-curve of product to simulating of polymer concentration interpretation model when applying described consideration to alter poly-passage is analyzed.
Preferably, described foundation considers that the step of polymer concentration interpretation model when altering poly-passage includes:
Determine the non-extension rate expression formula altering water relative polymer concentrations in poly-passage;
Polymer concentration interpretation model when not considering to alter poly-passage is multiplied with described extension rate, it is thus achieved that consideration is altered poly- Polymer concentration interpretation model during passage.
Preferably, described extension rate expression formula is:
α = T k × b μ p T k × b μ p + 1 × ( 1 - b ) μ w
In formula, α represents the extension rate of polymer concentration, and unit is f;Tk represents that altering poly-channel permeability rate advances by leaps and bounds coefficient, Unit is f;B represents that altering poly-passage takes up an area layer thickness ratio, and unit is f;μpRepresenting polymer solution viscosity, unit is mPa s;μwRepresenting aqueous viscosity, unit is mPa s.
Preferably, the expression formula of polymer concentration interpretation model when poly-passage is altered in described consideration is:
C C 0 = α 1 2 e r f c ( x - ut ′ 2 Dt ′ ) t ≤ Δ t α [ 1 2 e r f c ( x - ut ′ 2 Dt ′ ) - 1 2 e r f c ( x - ut ′ ′ 2 Dt ′ ′ ) ] t > Δ t
In formula,D is diffusion coefficient;C For polymer concentration, unit is mg/L;U is polymer percolation flow velocity, and unit is m/d;φ is porosity, and unit is f;φfFor Flowing porosity, unit is f;SoFor oil saturation, unit is f;SwcFor irreducible water saturation, unit is f;A is absorption system Number;ρrFor rock density, unit is g/cm3;T is the testing time, and unit is d;α represents output polymer concentration extension rate, single Position is f;Δ t is polymer slug injection length.
For achieving the above object, present invention also offers a kind of poly-drive during produce the quantitative analysis device of poly-curve, bag Include:
Taxon, for adding up the relation of different produced liquid in oil well polymer concentration and monitoring time, it is thus achieved that produce poly-song Line, and classify according to curvilinear characteristic curve poly-to described product;
First polymer concentration interpretation model sets up unit, for setting up polymer concentration solution when not considering to alter poly-passage Release model;
Second polymer concentration interpretation model sets up unit, explains for the polymer concentration when not considering to alter poly-passage On the basis of model, set up polymer concentration interpretation model when considering to alter poly-passage;
Interpretive analysis unit is for detecting the polymer concentration of producing well, dense according to polymer when considering to alter poly-passage Poly-curve is produced in the polymer concentration matching detected by degree interpretation model, and utilizes the classification results of the poly-curve of described product to plan The poly-curve of product closed out is classified, the product to simulating of polymer concentration interpretation model when applying described consideration to alter poly-passage Poly-curve is analyzed.
Preferably, described second polymer concentration interpretation model is set up unit and is included:
Extension rate determines module, for determining that the non-extension rate altering water relative polymer concentrations in poly-passage is expressed Formula;
Set up model module, polymer concentration interpretation model during for not considering to alter poly-passage and described extension rate It is multiplied, it is thus achieved that consider to alter polymer concentration interpretation model during poly-passage.
Preferably, described extension rate determines that the extension rate expression formula that module determines is:
α = T k × b μ p T k × b μ p + 1 × ( 1 - b ) μ w
In formula, α represents output polymer concentration extension rate, and unit is f;Tk represents to alter poly-channel permeability rate and advance by leaps and bounds Number, unit is f;B represents that altering poly-passage takes up an area layer thickness ratio, and unit is f;μpRepresenting polymer solution viscosity, unit is mPa·s;μwRepresenting aqueous viscosity, unit is mPa s.
Preferably, described second polymer concentration interpretation model sets up polymer when poly-passage is altered in the consideration that unit sets up The expression formula of concentrations account model is:
C C 0 = α 1 2 e r f c ( x - ut ′ 2 Dt ′ ) t ≤ Δ t α [ 1 2 e r f c ( x - ut ′ 2 Dt ′ ) - 1 2 e r f c ( x - ut ′ ′ 2 Dt ′ ′ ) ] t > Δ t
In formula,D is diffusion coefficient;C For polymer concentration, unit is mg/L;U is polymer percolation flow velocity, and unit is m/d;φ is porosity, and unit is f;φfFor Flowing porosity, unit is f;SoFor oil saturation, unit is f;SwcFor irreducible water saturation, unit is f;A is absorption system Number;ρrFor rock density, unit is g/cm3;T is the testing time, and unit is d;α represents the polymer concentration extension rate of output, Unit is f;Δ t is polymer slug injection length.
Technique scheme has the advantages that
The technical program, from producing poly-curvilinear characteristic, is gathered output polymer concentration and the relation of time during driving, is led to Crossing and set up polymer concentration interpretation model, poly-curve is produced in matching, obtains fluid-channeling channel and the crucial ginseng such as Polymer adsorption, diffusion Number, effectively can recognize oil reservoir by guide field so that be achieved by producing poly-log inversion fluid-channeling channel key parameter, formed Just can quickly identify the quantitative analysis method of fluid-channeling channel by altering poly-curve, for understanding oil reservoir fluid-channeling channel feature and Follow-up has stronger directive significance.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to Other accompanying drawing is obtained according to these accompanying drawings.
Fig. 1 be that the present invention proposes a kind of poly-drive during produce the quantitative analysis method flow diagram of poly-curve;
Fig. 2 is that flatbed produces poly-curve synoptic diagram;
Fig. 3 is that faint type produces poly-curve synoptic diagram;
Fig. 4 is that ramp type produces poly-curve synoptic diagram;
The geological model schematic diagram that Fig. 5 simplifies when being to consider to alter poly-passage;
Fig. 6 is one of poly-curve synoptic diagram of product that the present embodiment simulates;
Fig. 7 is the two of the poly-curve synoptic diagram of product that the present embodiment simulates;
Fig. 8 is the three of the poly-curve synoptic diagram of product that the present embodiment simulates;
Fig. 9 is that the advance by leaps and bounds dependency of coefficient of poly-channel thickness and permeability of altering that the poly-curve interpretation of dissimilar product obtains shows It is intended to;
Figure 10 be the poly-curve interpretation of dissimilar product obtain to alter poly-channel thickness relevant to alter poly-channel thickness ratio Property schematic diagram;
Figure 11 be that the present invention proposes a kind of poly-drive during produce the quantitative analysis device block diagram of poly-curve.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Describe, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments wholely.Based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under not making creative work premise Embodiment, broadly falls into the scope of protection of the invention.
The operation principle of the technical program is: the technical program, from producing poly-curvilinear characteristic, is produced according to during poly-driving Go out the relation of amount of polymers and time, classify to producing poly-curve, set up polymer concentration interpretation model, root on this basis Producing poly-curve according to the polymer concentration matching detected, according to classification results, curve poly-to the product simulated utilizes polymer Concentrations account model resolves, and has carried out quantitatively retouching to altering poly-passage, Polymer adsorption, polymer diffusion and seepage parameters State, for recognizing oil reservoir fluid-channeling channel feature and instructing follow-up corresponding regulating measures to have important theoretical and practical significance.
Based on above-mentioned operation principle, the present invention propose a kind of poly-drive during produce the quantitative analysis method of poly-curve, such as Fig. 1 Shown in.Including:
Step 101): the different produced liquid in oil well polymer concentration of statistics and the relation of monitoring time, it is thus achieved that produce poly-curve, and Classify according to curvilinear characteristic curve poly-to described product;
In a step 101, during adding up poly-driving, produce poly-concentration, according to curve shape by poly-for product curve classification, particularly as follows: The polymer concentration of the different produced liquid in oil well of statistics and the relation of monitoring time, abscissa is the production time, and vertical coordinate is polymerization Substrate concentration, is then divided into 3 classes by curvilinear characteristic: flatbed produces poly-curve, as shown in Figure 2 by poly-for product characteristic curve;Faint type produces Poly-curve, as shown in Figure 3;And the ramp type poly-curve of product, as shown in Figure 4.
Step 102): set up polymer concentration interpretation model when not considering to alter poly-passage;
Step 102 is specifically divided into following 4 processes:
Process 1: in the case of setting up One-Dimensional flows, it is considered to adsorbing polymer mass transfer diffusion background mathematics model, tool Body is:
D ∂ 2 C ∂ x 2 - u ∂ C ∂ x = [ φ ( 1 - S o ) + a ( 1 - φ ) ρ r ] φ f ( 1 - S o - S w c ) ∂ C ∂ t - - - ( 1 )
In formula (1), D is diffusion coefficient;C is polymer concentration, mg/L;U is polymer percolation flow velocity, m/d;φ is hole Degree, f;φfFor flowing porosity, f;SoFor oil saturation, f;SwcFor irreducible water saturation, f;A is adsorption coefficient;ρrFor rock Stone density, g/cm3;T is the testing time, d.
Process 2: during injection of polymer, it is considered to polymer starts continuous-stable from the t=0 moment and injects, one-dimensional feelings Under condition, boundary condition and initial condition are set to:
C ( x , 0 ) = C 0 x ≤ 0 0 x > 0 C ( 0 , t ) = C 0 t > 0 C ( ∞ , t ) = C 0 t > 0 - - - ( 2 )
Process 3: utilize Laplace conversion to solve polymer mass transfer diffusion background mathematics model in the case of One-Dimensional flows.
t ′ = φ f ( 1 - S o - S w c ) [ φ ( 1 - S o ) + a ( 1 - φ ) ρ r ] t - - - ( 3 )
One-dimensional Polymers mass transfer diffusion background mathematics model conversation is:
D ∂ 2 C ∂ x 2 - u ∂ C ∂ x = ∂ C ∂ t ′ - - - ( 4 )
The border initial condition of one-dimensional Polymers mass transfer diffusion basic data model is:
C ( x , 0 ) = C 0 x ≤ 0 0 x > 0 ; C ( 0 , t ) = C 0 t ′ > 0
Process 4: carry out mathematical model converting for the Laplace of time t', then solve and obtain Laplace space solution After, utilize Laplace inverse transformation to carry out inverting, obtaining analytic solutions is:
C C 0 = 1 2 e r f c ( x - ut ′ 2 Dt ′ ) + 1 2 e u x D e r f c ( x + ut ′ 2 Dt ′ ) - - - ( 5 )
It is typically small that formula (5) the right Section 2 compares Section 1, oil-producing well polymer when being then continuously injected into polymer Concentration is:
C C 0 = 1 2 e r f c ( x - ut ′ 2 Dt ′ ) - - - ( 6 )
For injection of polymer slug, from starting injecting timing, slug injection length is Δ t, then output polymer is dense Degree change expression formula is:
C C 0 = 1 2 e r f c ( x - ut ′ 2 Dt ′ ) - 1 2 e r f c ( x - ut ′ ′ 2 Dt ′ ′ ) - - - ( 7 )
In formula (7):
t ′ ′ = φ f ( 1 - S o - S w c ) [ φ ( 1 - S o ) + a ( 1 - φ ) ρ r ] ( t - Δ t ) - - - ( 8 )
Step 103): utilize polymer concentration interpretation model when not considering to alter poly-passage to set up when considering to alter poly-passage Polymer concentration interpretation model;
Step 103 specifically includes:
Consider that the geological model simplified when altering poly-passage is as shown in fig. 5, it is assumed that polymer is dense in Z-direction in altering poly-passage Spending identical, and think that polymer is migrated mainly along high bleed-through road, altering poly-channel thickness occupation of land layer thickness ratio is b, alters poly-logical Permeability coefficient of advancing by leaps and bounds in road is Tk.
After above simplification, it is considered to the polymer output concentrations account result altering poly-passage has only to consider non-altering The diluting effect of water relative polymer concentrations in poly-passage, extension rate expression formula is:
α = T k × b μ p T k × b μ p + 1 × ( 1 - b ) μ w - - - ( 9 )
In formula (9), α represents output polymer concentration extension rate, f;Tk represents that altering poly-channel permeability rate advances by leaps and bounds coefficient, f;B represents that altering poly-passage takes up an area layer thickness ratio, f;μpRepresent polymer solution viscosity mPa s;μwRepresent aqueous viscosity, mPa·s。
Utilize formula (9) that polymer output concentration expression formula can be obtained to be:
C C 0 = α 1 2 e r f c ( x - ut ′ 2 Dt ′ ) t ≤ Δ t α [ 1 2 e r f c ( x - ut ′ 2 Dt ′ ) - 1 2 e r f c ( x - ut ′ ′ 2 Dt ′ ′ ) ] t > Δ t
Step 104): the polymer concentration of detection producing well, explain mould according to polymer concentration when considering to alter poly-passage Poly-curve is produced in the polymer concentration matching detected by type, and utilizes the classification results of the poly-curve of the described product product to simulating Poly-curve is classified, and the poly-curve of product to simulating of polymer concentration interpretation model when applying described consideration to alter poly-passage enters Row is analyzed.
Accordingly, based on above-mentioned operation principle, the present invention also propose a kind of poly-drive during produce the quantitative analysis of poly-curve Device, as shown in figure 11.Including:
Taxon a, for adding up the relation of different produced liquid in oil well polymer concentration and monitoring time, it is thus achieved that produce poly-song Line, and classify according to curvilinear characteristic curve poly-to described product;
First polymer concentration interpretation model sets up unit b, for setting up polymer concentration when not considering to alter poly-passage Interpretation model;
Second polymer concentration interpretation model sets up unit c, for utilizing polymer concentration when not considering to alter poly-passage Interpretation model sets up polymer concentration interpretation model when considering to alter poly-passage;
Interpretive analysis unit d is for detecting the polymer concentration of producing well, dense according to polymer when considering to alter poly-passage Poly-curve is produced in the polymer concentration matching detected by degree interpretation model, and utilizes the classification results of the poly-curve of described product to plan The poly-curve of product closed out is classified, the product to simulating of polymer concentration interpretation model when applying described consideration to alter poly-passage Poly-curve is analyzed.
In order to the technical program application effect has understanding more intuitively, existing to use, certain Conglomerate Reservoir of said method is poly-to be driven The detailed description of the invention of the present invention is described as a example by process.
Table 1 gather for the Conglomerate Reservoir that embodiment provides drive during the underlying parameter of 3 mouthfuls of typical oil-producing wells, remaining parameter For: well spacing is all 125m, and formation water viscosity is 0.6mPa.s, and rock density is 2.4g/cm3, porosity 0.18, flowing porosity 0.144, irreducible water saturation 0.35, current oil saturation 0.5:
Table 1
Well-name Polymer viscosity mPa.s In the detection time, average note adopts speed t/d Individual well average moisture content f Perforation effective thickness m
A1 1294.64 13.6 0.882 4.25
A2 1294.64 8.6 0.779 3.5
A3 1311.43 6.6 0.624 1
The polymer concentration detecting producing well is fitted, it is thus achieved that the product poly-curve signal shown in Fig. 6, Fig. 7 and Fig. 8 Figure.
Table 2 is to resolve the key parameter about Polymer adsorption, diffusion and fluid-channeling channel obtained:
Table 2
Well-name Adsorption coefficient a Diffusion coefficient D Alter poly-channel thickness ratio b/f Alter poly-channel permeability rate to advance by leaps and bounds coefficient/f
A1 0.025 12 0.3 40
A2 0.015 13 0.2 20
A3 0.03 13 0.7 60
As can be seen from Table 2: flatbed alter poly-curve corresponding alter poly-channel thickness ratio and permeability is advanced by leaps and bounds coefficient Medium, and ramp type alter poly-curve corresponding alter poly-channel thickness ratio and permeability coefficient of advancing by leaps and bounds is relatively large.
In order to the application effect of the technical program has understanding more intuitively, existing to use certain conglomerate oil of technique scheme Hide poly-drive process as a example by the detailed description of the invention of the present invention is described.
Fig. 9 and Figure 10 is poly-to drive the poly-curve quantitative analysis of producing of 121 wells time in trial zone and explain the fluid-channeling channel ginseng obtained Number and the dependency altering poly-curve type.Wherein,
In fig .9, alter poly-channel thickness and the permeability that demonstrate that the poly-curve interpretation of dissimilar product obtains are advanced by leaps and bounds coefficient Dependency, as seen from the figure: produce poly-curve type with to alter poly-channel permeability rate and advance by leaps and bounds Relativity of Coefficients preferable: advance by leaps and bounds coefficient great Yi Forming ramp type and alter poly-rule, coefficient of advancing by leaps and bounds is little to be easily formed faint type and alters poly-curve.
In Fig. 10, demonstrate that what the poly-curve interpretation of dissimilar product obtained alters poly-channel thickness and alter poly-channel thickness ratio The dependency of example, as seen from the figure: produce poly-curve type with to alter poly-channel thickness ratio dependency preferable: thickness proportion great Yi shape Poly-rule altered by sloping type, and thickness proportion is little to be easily formed faint type and alter poly-curve.
It will be apparent from the above that: flatbed produces poly-curve and mainly reflects that altering poly-passage grows relatively;It is main that faint type produces poly-curve It is relatively weak that poly-developmental pathway is altered in reflection;Ramp type is produced poly-curve and is mainly reflected that alter poly-passage grows very much, sequent water flooding process In, polymer concentration the most quickly declines.
The technical program solve poly-during driving by produce poly-curve quantitative analysis alter poly-passage and Polymer adsorption and The key parameters such as diffusion so that be achieved by crucial parameter during analyzing the poly-curve quantitative analysis polymer flooding of product.
Above detailed description of the invention, has been carried out the most specifically the purpose of the present invention, technical scheme and beneficial effect Bright, be it should be understood that the detailed description of the invention that these are only the present invention, the protection model being not intended to limit the present invention Enclose, all within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. done, should be included in the present invention Protection domain within.

Claims (8)

1. one kind poly-drive during produce the quantitative analysis method of poly-curve, it is characterised in that including:
The different produced liquid in oil well polymer concentration of statistics and the relation of monitoring time, it is thus achieved that produce poly-curve, and according to curvilinear characteristic Curve poly-to described product is classified;
Set up polymer concentration interpretation model when not considering to alter poly-passage;
Polymer concentration interpretation model when not considering to alter poly-passage is utilized to set up polymer concentration solution when considering to alter poly-passage Release model;
The polymer concentration of detection producing well, poly-to detect according to the polymer concentration interpretation model considered when altering poly-passage Poly-curve is produced in compound concentration matching, and utilizes the classification results of the poly-curve of the described product poly-curve of product to simulating to carry out point Class, the poly-curve of product to simulating of polymer concentration interpretation model when applying described consideration to alter poly-passage is analyzed.
2. the method for claim 1, it is characterised in that polymer concentration when described foundation considers to alter poly-passage is explained The step of model includes:
Determine the non-extension rate expression formula altering water relative polymer concentrations in poly-passage;
Polymer concentration interpretation model when not considering to alter poly-passage is multiplied with described extension rate, it is thus achieved that consider to alter poly-passage Time polymer concentration interpretation model.
3. method as claimed in claim 2, it is characterised in that described extension rate expression formula is:
α = T k × b μ p T k × b μ p + 1 × ( 1 - b ) μ w
In formula, α represents the extension rate of polymer concentration, and unit is f;Tk represents that altering poly-channel permeability rate advances by leaps and bounds coefficient, unit For f;B represents that altering poly-passage takes up an area layer thickness ratio, and unit is f;μpRepresenting polymer solution viscosity, unit is mPa s;μwTable Showing aqueous viscosity, unit is mPa s.
4. method as claimed in claim 2 or claim 3, it is characterised in that polymer concentration when poly-passage is altered in described consideration is explained The expression formula of model is:
C C 0 = α 1 2 e r f c ( x - ut ′ 2 Dt ′ ) t ≤ Δ t α [ 1 2 e r f c ( x - ut ′ 2 Dt ′ ) - 1 2 e r f c ( x - ut ′ ′ 2 Dt ′ ′ ) ] t > Δ t
In formula,D is diffusion coefficient;C is poly- Compound concentration, unit is mg/L;U is polymer percolation flow velocity, and unit is m/d;φ is porosity, and unit is f;φfFor flowing Porosity, unit is f;SoFor oil saturation, unit is f;SwcFor irreducible water saturation, unit is f;A is adsorption coefficient;ρr For rock density, unit is g/cm3;T is the testing time, and unit is d;α represents the polymer concentration extension rate of output, unit For f;Δ t is polymer slug injection length.
5. one kind poly-drive during produce the quantitative analysis device of poly-curve, it is characterised in that including:
Taxon, for adding up the relation of different produced liquid in oil well polymer concentration and monitoring time, it is thus achieved that produce poly-curve, and Classify according to curvilinear characteristic curve poly-to described product;
First polymer concentration interpretation model sets up unit, explains mould for setting up polymer concentration when not considering to alter poly-passage Type;
Second polymer concentration interpretation model sets up unit, for the polymer concentration interpretation model when not considering to alter poly-passage On the basis of, set up polymer concentration interpretation model when considering to alter poly-passage;
Interpretive analysis unit, for detecting the polymer concentration of producing well, alters polymer concentration solution during poly-passage according to consideration Release model and the polymer concentration matching detected is produced poly-curve, and utilize the classification results of the poly-curve of described product to simulating The poly-curve of product classify, the poly-song of the product to simulating of polymer concentration interpretation model when applying described consideration to alter poly-passage Line is analyzed.
6. device as claimed in claim 5, it is characterised in that described second polymer concentration interpretation model sets up unit bag Include:
Extension rate determines module, for determining the non-extension rate expression formula altering water relative polymer concentrations in poly-passage;
Set up model module, polymer concentration interpretation model during for not considering to alter poly-passage and described extension rate phase Take advantage of, it is thus achieved that consider to alter polymer concentration interpretation model during poly-passage.
7. device as claimed in claim 6, it is characterised in that described extension rate determines that the extension rate that module determines is expressed Formula is:
α = T k × b μ p T k × b μ p + 1 × ( 1 - b ) μ w
In formula, α represents output polymer concentration extension rate, and unit is f;Tk represents that altering poly-channel permeability rate advances by leaps and bounds coefficient, single Position is f;B represents that altering poly-passage takes up an area layer thickness ratio, and unit is f;μpRepresenting polymer solution viscosity, unit is mPa s;μw Representing aqueous viscosity, unit is mPa s.
Device the most as claimed in claims 6 or 7, it is characterised in that described second polymer concentration interpretation model sets up unit The expression formula of polymer concentration interpretation model when poly-passage is altered in the consideration set up is:
C C 0 = α 1 2 e r f c ( x - ut ′ 2 Dt ′ ) t ≤ Δ t α [ 1 2 e r f c ( x - ut ′ 2 Dt ′ ) - 1 2 e r f c ( x - ut ′ ′ 2 Dt ′ ′ ) ] t > Δ t
In formula,D is diffusion coefficient;C is poly- Compound concentration, unit is mg/L;U is polymer percolation flow velocity, and unit is m/d;φ is porosity, and unit is f;φfFor flowing Porosity, unit is f;SoFor oil saturation, unit is f;SwcFor irreducible water saturation, unit is f;A is adsorption coefficient;ρr For rock density, unit is g/cm3;T is the testing time, and unit is d;α represents output polymer concentration extension rate, and unit is f;Δ t is polymer slug injection length.
CN201610515517.3A 2016-07-01 2016-07-01 A kind of poly- drive produces the quantitative analysis method and device of poly- curve in the process Expired - Fee Related CN106055843B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610515517.3A CN106055843B (en) 2016-07-01 2016-07-01 A kind of poly- drive produces the quantitative analysis method and device of poly- curve in the process

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610515517.3A CN106055843B (en) 2016-07-01 2016-07-01 A kind of poly- drive produces the quantitative analysis method and device of poly- curve in the process

Publications (2)

Publication Number Publication Date
CN106055843A true CN106055843A (en) 2016-10-26
CN106055843B CN106055843B (en) 2019-05-14

Family

ID=57201702

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610515517.3A Expired - Fee Related CN106055843B (en) 2016-07-01 2016-07-01 A kind of poly- drive produces the quantitative analysis method and device of poly- curve in the process

Country Status (1)

Country Link
CN (1) CN106055843B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110857627A (en) * 2018-08-22 2020-03-03 中国石油天然气股份有限公司 Unconsolidated sandstone reservoir seepage unit identification method and device

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090194276A1 (en) * 2008-01-31 2009-08-06 Total E&P Usa, Inc. Determination of an actual optimum salinity and an actual optimum type of microemulsion for surfactant/polymer flooding
CN104564049A (en) * 2015-01-22 2015-04-29 中国石油大学(华东) Detection method of polymer flooding and gathering channel
CN104747149A (en) * 2015-01-22 2015-07-01 中国石油大学(华东) Detecting method of binary combined flooding agent migration channel
CN105715241A (en) * 2016-01-15 2016-06-29 中国石油大学(华东) Method for measuring polymer flooding relative permeability curve

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090194276A1 (en) * 2008-01-31 2009-08-06 Total E&P Usa, Inc. Determination of an actual optimum salinity and an actual optimum type of microemulsion for surfactant/polymer flooding
CN104564049A (en) * 2015-01-22 2015-04-29 中国石油大学(华东) Detection method of polymer flooding and gathering channel
CN104747149A (en) * 2015-01-22 2015-07-01 中国石油大学(华东) Detecting method of binary combined flooding agent migration channel
CN105715241A (en) * 2016-01-15 2016-06-29 中国石油大学(华东) Method for measuring polymer flooding relative permeability curve

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘同敬等: "井间示踪剂测试半解析方法体系数学模型", 《石油学报》 *
程林松: "《高等渗流力学》", 30 November 2011 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110857627A (en) * 2018-08-22 2020-03-03 中国石油天然气股份有限公司 Unconsolidated sandstone reservoir seepage unit identification method and device
CN110857627B (en) * 2018-08-22 2022-10-04 中国石油天然气股份有限公司 Method and device for identifying seepage unit of unconsolidated sandstone reservoir

Also Published As

Publication number Publication date
CN106055843B (en) 2019-05-14

Similar Documents

Publication Publication Date Title
CN106837297B (en) Method for identifying connectivity among wells and predicting oil-water dynamic state
CN104100263B (en) Method for determining residual oil saturation of non-uniform water flooded layer
EP3800324B1 (en) Method and apparatus for determining integrated development approach for shale and adjacent oil layers
Parekh et al. A case study of improved understanding of reservoir connectivity in an evolving waterflood with surveillance data
CN106932324B (en) A kind of method of determining High water cut sandstone oil reservoir reservoir permeability changing rule
CN104632157A (en) Low permeability reservoir equilibrium displacement method
CN105631078B (en) The method for numerical simulation of Natural fractures system self-adapting medium
CN104060985A (en) Method and system for testing entering depth of stratified oil deposit profile control water plugging agent
CN106703768B (en) A kind of determination method of simulation high-concentration and low-concentration polymer flooding alternating injection timing
CN104912537A (en) Well network deployment method of giant thick fractured and buried hill reservoir
CN107448197A (en) A kind of method of marine facies sandstone oil reservoir fluid-channeling channel space development quantitative judge
Todd et al. An Evaluation of EOR Potential in the Elm Coulee Bakken Formation, Richland County, Montana
CN103334725A (en) Method and device for evaluating low permeability reservoir displacement effectiveness
CN106779229A (en) A kind of High water cut water plugging of horizontal well Potential Prediction method
CN106055843A (en) Quantitative analysis method and device of polymer production curve in polymer flooding process
CN111155980B (en) Water flow dominant channel identification method and device
CN104834934B (en) A kind of nucleome capture method for being used to identify reservoir fluid
CN105089659B (en) A kind of Conglomerate Reservoir permeable unit recognition methods
CN104570104B (en) A kind of vertical shear wave earthquake quality factor extracting method based on two-step method AVF
CN105092785B (en) A kind of injection well injectability evaluation method considering polymer viscoelastic
CN107832900A (en) Sandy conglomerate bodies Evaluation of Affusion Effect method
CN106593409A (en) Recognition method for futile cycle in oilfield injection-production relation
CN113971528A (en) Method for identifying dominant seepage channel
Langenbruch et al. Probability of brittle rock failure during hydraulic fracturing of conventional and unconventional reservoirs
CN113032996B (en) Water channeling channel identification method for hypotonic fractured reservoir horizontal well

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190514