CN105373643B - A kind of leading edge simulation monitoring method that injection of polymer promotes - Google Patents
A kind of leading edge simulation monitoring method that injection of polymer promotes Download PDFInfo
- Publication number
- CN105373643B CN105373643B CN201510537423.1A CN201510537423A CN105373643B CN 105373643 B CN105373643 B CN 105373643B CN 201510537423 A CN201510537423 A CN 201510537423A CN 105373643 B CN105373643 B CN 105373643B
- Authority
- CN
- China
- Prior art keywords
- polymer
- concentration value
- injection
- value
- sample point
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to the leading edges that a kind of injection of polymer promotes to simulate monitoring method, belongs to oil recovering technical field.The present invention on experimental model sample point first by acquiring polymer solution concentration value；Then it is directed to the polymer concentration value that polymer solution concentration value on each sample point uses other positions of Kriging regression algorithm prediction model；The monitoring simulated to polymer advancing front is realized according to the polymer concentration value at predicted other positions of model.The present invention is in the way of mathematical geology, with Kriging regression method, according to the concentration value of several information samples inside and outside a region, a kind of estimation of linear unbias and mini-mum estimation variance is made to the region, to realize the monitoring simulated to whole region injection of polymer advancing front, and it can be adapted for the monitoring of cross-linked polymer advancing front simulation.
Description
Technical field
The present invention relates to the leading edges that a kind of injection of polymer promotes to simulate monitoring method, belongs to oil recovering technical field.
Background technology
Continental deposit oil reservoir by accumulating condition because influenced, and there is one for physical property distribution of the reservoir in interlayer and layer
Determine the difference of degree.The oil reservoir of waterflooding extraction, in development late stage reservoir for a long time by the water eroding and flooding of outside fluid, degree of consolidation phase
It migrates therewith to poor clay mineral and part skeleton particle and exacerbates the anisotropism of reservoir, in addition unfavorable grease stream
Ratio is spent, causes injection liquid to advance by leaps and bounds along high permeability zone, macropore channelling, leads to the too early water logging of oil reservoir, remaining oil saturation is relatively
High less permeable layer cannot be utilized reasonably, reduce water flooding effectiveness.Such as Henan Oil Field priority since " 95 "
In Shuan He, Xia Ermen, ancient city etc., totally 16 series of strata have carried out polymer flooding and secondary poly- drive industrial applications, at present polymer flooding
And one of the agent technology that crosslinked polymer flooding technology has become oil field development stable yields and delays production decline.But to polymer
Or the numerical simulation study of crosslinked polymer flooding, it establishes on the basis of two-dimentional two-phase and D beam element mostly, function is simpler
It is single, it is difficult to meet requirement of the Oil Field with the extensive fine reservoir polymer leading edge monitoring that block is entirety.
Invention content
The object of the present invention is to provide the leading edges that a kind of injection of polymer promotes to simulate monitoring method, to meet Oil Field
With the requirement that block is whole extensive fine reservoir polymer leading edge monitoring.
The present invention provides a kind of leading edge simulation monitoring method that injection of polymer promotes to solve above-mentioned technical problem, should
Monitoring method includes the following steps：
1) polymer solution concentration value is acquired on three dimensional physical cast experimental model sample point；
2) it is directed to polymer solution concentration value on each sample point and uses other positions of Kriging regression algorithm prediction model
Polymer concentration value；
3) prison simulated to polymer advancing front is realized according to the polymer concentration value at predicted other positions of model
It surveys.
The calculating process of Kriging method is as follows in the step 2)：
A. using the polymer solution concentration value acquired on sample point as initial data, experiment half is solved in a single direction and is become
Different function gamma*(h), by calculating γ*(h) value is mapped, and chooses theoretical model, is found out base station and is become journey；
B. the division that super block and time block are carried out according to the information state in region to be estimated, carries out returning block to Given information point；
C. centered on super block to be estimated, increase several super blocks to the periphery as a neighborhood, carry out sliding estimation；
D. according to parameter base station, become the model of journey and selection, golden equation in foundation gram；
E. it is that constraint is adjusted with unbiasedness and mini-mum estimation variance, golden equation, makes estimated value in solution gramMost preferably, it solves
It is region to be estimated to obtain result
The semivariable function γ*(h) formula is：
Wherein Z (xi) it is xiThe observation at place, i.e., in sample point xiLocate the polymer solution concentration value of acquisition；N (h) is meter
Calculate the number of numerical point, the i.e. number of sample point.
Equation employed in the step D is：
WhereinTo be estimated region true value ZVEstimated value；λiFor gram in golden weight coefficient；ZiFor in estimated value neighborhood
Observation.
The step 1) is to simulate reservoir plane heterogeneity using plain heterogeneity rock core, in 1 mouthful of extraction well direction
High infiltration strip is set, and injection of polymer (or cross-linked polymer) solution carries out the displacement of reservoir oil, dense in each sample point acquisition polymer along journey
Angle value simultaneously monitors pressure.
The beneficial effects of the invention are as follows:The present invention on experimental model sample point first by acquiring polymer solution concentration
Value；Then it is dense using the polymer at other positions of Kriging regression algorithm prediction model to be directed to polymer solution value on each sample point
Angle value；The monitoring simulated to polymer advancing front is realized according to the polymer concentration value at predicted other positions of model.This
Invention is in the way of mathematical geology, with Kriging regression method, according to the concentration of several information samples inside and outside a region
Value, a kind of estimation of linear unbias and mini-mum estimation variance is made to the region, to realize to whole region injection of polymer
The monitoring of (or cross-linked polymer) advancing front simulation.
Description of the drawings
Fig. 1 is the schematic diagram that selected three dimensional physical pours model in the embodiment of the present invention；
Fig. 2 is the note cross-linked polymer pressure-plotting at ten o'clock in the embodiment of the present invention；
Fig. 3 is the note cross-linked polymer pressure-plotting on ten fours in the embodiment of the present invention；
Fig. 4 is the note cross-linked polymer pressure-plotting of ten eights in the embodiment of the present invention；
Fig. 5 is the note cross-linked polymer pressure-plotting of two twelve-hours in the embodiment of the present invention；
Fig. 6 is the note cross-linked polymer pressure-plotting on two in the embodiment of the present invention；
Note cross-linked polymer concentration profile when Fig. 7 is the first sub-sampling in the embodiment of the present invention；
Note cross-linked polymer concentration profile when Fig. 8 is the second sub-sampling in the embodiment of the present invention；
Note cross-linked polymer concentration profile when Fig. 9 is third sub-sampling in the embodiment of the present invention.
Specific implementation mode
The specific implementation mode of the present invention is described further below in conjunction with the accompanying drawings.
The features such as it is high that the present invention be directed to oil field permeabilities, and anisotropism is strong, in order to understand fully in polymer flooding injection process
Displacing fluid involves range problem, the leading edge simulation monitoring method that a kind of injection of polymer proposed promotes.This method is used
Kriging regression method makes one kind according to certain characteristic data values of several information samples inside and outside a region to the region
The estimation of linear unbias and mini-mum estimation variance, to realize the monitoring to whole region.
Below with plain heterogeneity rock core (60*60*4.5cm) simulation reservoir plane heterogeneity (note four is adopted), this reality
It is that the three dimensional physical that a note four is adopted pours model to apply the model chosen in example, and model is as shown in Figure 1, base station and change journey are ground
Basic conception in matter statistics, base station refer to the particular value of experimental data, can voluntarily refer to according to the requirement of experiment of oneself
The fixed numerical value (base station in the present embodiment refers to the concentration value that sampled point measures), it is other regions (in addition to sampled point) warp to become journey
Cross the difference of centrifugal pump and base station value that Kriging method obtains.Permeability grade investigates 2,4 times two kinds, thereto 1 mouthful of extraction
High infiltration strip is arranged in well direction, and model high infiltration strip and basic permeability choose 1200/600 × 10-3μm2, injection intensity meets
It is required that cross-linked polymer carry out the displacement of reservoir oil, by along journey close sampling analyze and monitor pressure, determine that cross-linked polymer slug pushes away
Into leading edge and its regularity of distribution.
1. the polymer solution concentration value (mg/l) acquired on experimental model sample point
(1) prepare rock core and installation equipment
1. with sand paper polishing rock core, ensure that core surface is smooth, cleans surface with alcohol, and dry rock core, measure rock core
Size calculates apparent volume；
2. in core surface brush resinous coat, for preventing resin in casting process from largely penetrating into rock core；
3. in 5 mouthfuls of wells of upper surface drilling well of rock core, 8 pressure taps, as shown in Figure 1；
4. then rinsing rock core wellhole with clear water, and by rock core with hydrofluoric acid treatment pit shaft to prevent powder from blocking the borehole wall
It is placed on drying in oven；
5. rock core well location is consolidated, rock core is put into clamper, thereto casting resin, is poured and covers clamper
Model is pushed into baking oven, waits for hardening of resin by lid wherein 8 pressure taps connect sensor.
(2) initial oil, water saturation are established
1. vacuum pump is shifted to another a well, until all well locations by arbitrary a well to rock core vacuumizing to -0.1MPa
It is evacuated to -0.1MPa；
2. hand pump is connected to injection well, other four mouthfuls of wells are closed, are saturated water flooding, presses and stops until pressure gauge rises
After injecting a period of time, pressure does not decline, and stops injection, calculates the volume of injection water, i.e. rock pore volume；
3. being outlet water filling with constant flow pump by entrance and arbitrary a bite producing well of injection well, simultaneously closes off other three mouthfuls and adopt
Oil well, the water filling into rock core record smooth pressure until pressure is steady, calculate rock core along this extraction well direction matrix permeability, separately
Outer three mouthfuls of wells also repeat top-operation and measure permeability；
4. carrying out saturated oils operation, simulation oil is injected in injection well constant flow pump, waits for that arbitrary a bite extraction well is no longer discharged
And it is constant it is fuel-displaced after close the well, until whole well locations saturation is fully completed, calculate saturation oil volume, and under the conditions of 81 DEG C
Aging is more than 72 hours.
(3) experimental section
1. after abundant aging, a water drive is carried out with the flow of 1mL/min, each well connects Produced Liquid with test tube, and calculating contains
Water rate changes, while all extraction well outlet ends are placed on same test tube every half an hour during water drive and measure rock core synthesis
Moisture content, until composite water cut reaches 98%, one time water drive terminates；
2. injecting 0.1PV cross-linked polymer solutions with the flow of 0.5mL/min, solidifying, each well use in experimentation is fully waited
Test tube connects Produced Liquid, calculates water-cut variation, while pressure point pressure change, and close sampling are respectively adopted in monitoring, is surveyed per half an hour
Determine rock core composite water cut；
3. carrying out subsequent waterflooding with the flow of 1mL/min, monitored during water drive and respectively adopt pressure point pressure change, and intensive taken
Sample, while comprehensive rock core composite water cut is measured every half an hour, until reaching 98%, experiment terminates；
4. measuring polymer concentration in each pressure tap and the sampling of each well, matched curve, as Figure 7-9, concentration value is really
Surely be CdI_2-starch method measure sample point Produced Liquid polymer concentration value.
Respectively (10：00-10：15,14：00-14：15,18：00-18：15,22：00-22：15,2：00-2：15) it opens
No. 1, No. 3, No. 5, No. 7 inner ring pressure taps, take liquid 15 minutes；(10：15-10：30,14：15-14：30,18：15-18：30,22：
15-22：30,2：15-2：30) No. 2, No. 4, No. 6,8 extra ring pressure taps are opened, take liquid 15 minutes；Summarize drafting injection pressure
Distribution map, as shown in figures 2-6.
2. being finally inversed by the poly- of other positions of model using Kriging regression algorithm according to polymer solution value on each sample point
Close object concentration value.
1) to obtain the polymer concentration value of each sample point as initial data, which is carried out at extra-high value
Reason, checks whether Normal Distribution；
2) experiment with computing semivariable function, n direction；
Variation function is basic tool possessed by geostatistics, it both can be used for describing the space of regionalized variable
Structure change, and its random variation can be described.Assuming that spatial point x only changes in one-dimensional x-axis, regionalized variable Z (x)
The half of the difference of value at 2 points of x, x+h is defined as the semivariable functions of Z (x) in the direction of the x axis, i.e.,：
γ (x, h)=E (Z (x)-Z (x+h))/2
It is realistic on single direction to test semivariable function γ when meeting stationary hypothesis*(h) calculation formula is：
Z (x in formulai)—xiThe observation at place；
The number of N (h)-evaluation points
By calculating γ*(h) value is mapped, and chooses theoretical model (spherical, exponential model etc.), and find out base station, becomes journey.
3) division of super block and time block, and the just starting per places block and end are carried out according to the information state in region to be estimated
Only number file is established, and to carrying out returning block to Given information point.
4) centered on super block to be estimated, increase several super blocks to the periphery as a neighborhood, be 8 in the present embodiment
It is a, carry out sliding estimation；
5) according to parameter base station, become the model of journey and selection, seek γ (xi,xj) or r (xi, V), establish equation group.
Linear estimator is：
WhereinTo be estimated region true value ZVEstimated value；λiFor gram in golden weight coefficient；ZiFor in estimated value neighborhood
Observation；Make estimated valueMost preferably, it is necessary to meet following two condition：
Unbiasedness
Estimate varianceIt is minimum
Pass through equation group：
It acquires：
Wherein λiFor gram in golden weight coefficient；γ(xi,xj) it is observation xiWith xjBetween semivariable function value；r(xi,V)
For observation and the average variation function being estimated in area V between all the points；R (V, V) is to be estimated being averaged for all the points in area V
Semivariable function；μ is Lagrangian constant.
Claims (4)
1. the leading edge that a kind of injection of polymer promotes simulates monitoring method, which is characterized in that the monitoring method includes the following steps：
1) polymer solution concentration value is acquired on three dimensional physical cast experimental model sample point；
2) it is directed to the polymerization that polymer solution concentration value on each sample point uses other positions of Kriging regression algorithm prediction model
Object concentration value；
3) monitoring simulated to polymer advancing front is realized according to the polymer concentration value at predicted other positions of model；
The calculating process of Kriging method is as follows in the step 2)：
A. using the polymer solution concentration value acquired on sample point as initial data, experiment semivariation letter is solved in a single direction
Number γ*(h), by calculating γ*(h) value is mapped, and chooses theoretical model, is found out base station and is become journey；
B. the division that super block and time block are carried out according to the information state in region to be estimated, carries out returning block to Given information point；
C. centered on super block to be estimated, increase several super blocks to the periphery as a neighborhood, carry out sliding estimation；
D. according to parameter base station, become the model of journey and selection, golden equation in foundation gram；
E. it is that constraint is adjusted with unbiasedness and mini-mum estimation variance, golden equation, makes estimated value in solution gramMost preferably, solution obtains
As a result it is region to be estimated.
2. the leading edge that injection of polymer according to claim 1 promotes simulates monitoring method, which is characterized in that described half becomes
Different function gamma*(h) formula is：
Wherein Z (xi) it is xiThe observation at place, i.e., in sample point xiLocate the polymer solution concentration value of acquisition；N (h) is to calculate number
It is worth the number of point, the i.e. number of sample point.
3. the leading edge that injection of polymer according to claim 2 promotes simulates monitoring method, which is characterized in that the step
Equation employed in D is：
WhereinTo be estimated region true value ZVEstimated value；λiFor gram in golden weight coefficient；ZiFor the observation in estimated value neighborhood
Value.
4. the leading edge that injection of polymer according to claim 1 promotes simulates monitoring method, which is characterized in that the step
1) it is that reservoir plane heterogeneity is simulated using plain heterogeneity rock core, high infiltration strip is set in extraction well direction, injection is strong
It spends satisfactory transfer drive polymerization-filling object and carries out the displacement of reservoir oil, polymer concentration value is acquired in each sample point along journey and monitor pressure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510537423.1A CN105373643B (en) | 2015-08-27 | 2015-08-27 | A kind of leading edge simulation monitoring method that injection of polymer promotes |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510537423.1A CN105373643B (en) | 2015-08-27 | 2015-08-27 | A kind of leading edge simulation monitoring method that injection of polymer promotes |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105373643A CN105373643A (en) | 2016-03-02 |
| CN105373643B true CN105373643B (en) | 2018-09-04 |
Family
ID=55375839
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510537423.1A Active CN105373643B (en) | 2015-08-27 | 2015-08-27 | A kind of leading edge simulation monitoring method that injection of polymer promotes |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105373643B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112983401B (en) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | Boundary determination method for water invasion of boundary water gas reservoir |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104615862A (en) * | 2015-01-14 | 2015-05-13 | 中国石油天然气股份有限公司 | High water-cut oilfield well position determining method based on evolutionary algorithm |
| CN104809363A (en) * | 2015-05-25 | 2015-07-29 | 中国石油大学(华东) | Three-component combination drive optimization method based on BK approximate dynamic programming |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2001280934A1 (en) * | 2000-07-28 | 2002-02-13 | Alliance Pharmaceutical Corp. | Methods and compositions to upregulate, redirect or limit immune responses to bioactive compounds |
-
2015
- 2015-08-27 CN CN201510537423.1A patent/CN105373643B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104615862A (en) * | 2015-01-14 | 2015-05-13 | 中国石油天然气股份有限公司 | High water-cut oilfield well position determining method based on evolutionary algorithm |
| CN104809363A (en) * | 2015-05-25 | 2015-07-29 | 中国石油大学(华东) | Three-component combination drive optimization method based on BK approximate dynamic programming |
Non-Patent Citations (2)
| Title |
|---|
| Mapping soil organic matter using the topographic wetness index: A comparative study based on different flow-direction algorithms and kriging methods;Pei T 等;《Ecological Indicators》;20101231;第10卷(第3期);第610-619页 * |
| 油价变化下聚合物驱提高原油采收率动态规划方法研究;彭绍文;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20150615(第06期);B019-173 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105373643A (en) | 2016-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lu et al. | Horizontal hydraulic conductivity of shallow streambed sediments and comparison with the grain‐size analysis results | |
| CN107462936B (en) | Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law | |
| CN103352680B (en) | Foam based on the integration of pit shaft oil reservoir is handled up sediment outflow experimental facilities and method | |
| CN106522928B (en) | Well testing method for unstable pressure drop of well logging head by stopping pump after acidizing and fracturing | |
| CN110096718B (en) | Method for obtaining volume of karst cave in carbonate reservoir | |
| CN110263453B (en) | Modeling method of self-adaptive barrier lake gradual collapse coupling prediction model | |
| CN102809642A (en) | Method for determining hydrogeological parameters of aquitard | |
| CN110107277B (en) | Method for obtaining volume of karst cave encountered by well drilling in carbonate reservoir | |
| CN105369812B (en) | The pressure water head of foundation ditch pumping with constant rate determines method under suspension type water-stop curtain | |
| CN111189507B (en) | Karst water stratification discrimination and underground water level observation method | |
| CN109916447A (en) | A kind of set composite and its method for monitoring side slope saturation and sliding deformation | |
| CN106501156B (en) | Scene determines the outer tube drawdown double-tube method of outer tube aquiclude hydrogeological parameter | |
| CN108241772A (en) | Consider Tunnel Gushing method for determination of amount in the crack artesian aquifer of multiple factors | |
| CN209432184U (en) | A kind of set composite for monitoring side slope saturation and sliding deformation | |
| CN105373643B (en) | A kind of leading edge simulation monitoring method that injection of polymer promotes | |
| Rushton et al. | Significance of hydraulic head gradients within horizontal wells in unconfined aquifers of limited saturated thickness | |
| Lu et al. | Determination of the anisotropy of an upper streambed layer in east-central Nebraska, USA | |
| CN105386430B (en) | A kind of lower water-stop curtain both sides water level method for determining difference of water-stop curtain effect | |
| CN106761674B (en) | It is a kind of integrate draw water and the ground of unrestrained water test lead directly to down-hole drilling device | |
| CN105297739B (en) | The determination method of pressure water head when the default water level of foundation ditch draws water under water-stop curtain | |
| CN103195046A (en) | Method for embedding multiple pore water pressure gauges through one hole | |
| CN106680175B (en) | Scene determines the outer tube drawdown double-tube method of inner tube aquiclude hydrogeological parameter | |
| JP5882159B2 (en) | Simple measurement method for soil saturation | |
| CN113945497B (en) | Method for evaluating mobility of reservoir fluid of oil and gas reservoir | |
| CN106680177A (en) | Inner tube drawdown double tube method of determining hydrogeological parameters of aquitard of inner tube on site |
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 |