CN105626007B - Water times number calculating method is crossed based on different parts in rock core scale oil reservoir - Google Patents
Water times number calculating method is crossed based on different parts in rock core scale oil reservoir Download PDFInfo
- Publication number
- CN105626007B CN105626007B CN201410641714.0A CN201410641714A CN105626007B CN 105626007 B CN105626007 B CN 105626007B CN 201410641714 A CN201410641714 A CN 201410641714A CN 105626007 B CN105626007 B CN 105626007B
- Authority
- CN
- China
- Prior art keywords
- streamline
- water
- different parts
- rock core
- grid
- 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
Abstract
The present invention provides one kind and crosses water times number calculating method based on different parts in rock core scale oil reservoir, should cross water times number calculating method based on different parts in rock core scale oil reservoir and include：Step 1, five-spot uniform soft soil base streamline equation is derived, obtains different parts streamline distribution；Step 2, flow allocation method between streamline is established between injection-production well, obtains different parts flow distribution；Step 3, it calculates and water multiple is crossed along streamline different parts based on the one-dimensional uniform soft soil base of rock core scale；And step 4, calculate five-spot uniform soft soil base crosses water multiple along streamline different parts.This based on different parts in rock core scale oil reservoir cross water times number calculating method make full use of laboratory test abundant materials accurate description oil reservoir different parts cross water multiple, the displacement degree of quantization signifying oil reservoir different parts establishes solid foundation to formulate mining site modification scenario technical measures.
Description
Technical field
The present invention relates to oil field development technical fields, especially relate to a numerical simulator or practical oil reservoir one
Determine to cross water times number calculating method based on rock core scale different parts under Injection volume.
Background technology
Amount of water injected is by total accumulative injection rate of oil reservoir and the ratio of reservoir pore volume, reflects oil reservoir macroscopic view water drive
Situation；It is total ratio for adding up injection rate and rock pore volume by conventional core (Φ 2.5cm × 5cm) to cross water multiple,
Reflect microcosmic water drive situation.The amount of water injected for crossing entire oil reservoir in water multiple and traditional sense under the conditions of indoor Physical Experiment is not
It is a concept, is not directly applicable to mining site Development Practice.
Under the conditions of certain amount of water injected, the water multiple of crossing of oil reservoir different parts is mainly influenced by both sides factor：When
The location of in flooding pattern, water multiple height is crossed at main stream line position, and it is low that water multiple is crossed at non-master streamline position；Second is that it calculates single
The size of first (hereinafter referred to as infinitesimal) scale, the water multiple excessively at the small same position of infinitesimal scale is big, infinitesimal scale position same greatly
Cross water multiple it is small.Since position of any infinitesimal in well pattern is fixed, and therefore, the selection of scale is to crossing water multiple meter
It plays a crucial role.
Rock core displacement test is conventional, is also important a kind of laboratory facilities.Many weights can be obtained by the experiment
The useful data wanted crosses water multiple and oil displacement efficiency correspondence such as different；But in these logging data applications to practical oil reservoir
When there are bigger difficulty.What reason was that conventional core displacement experiment obtains cross, and water multiple generally reaches hundreds of or even tens of thousands of
Times, oil displacement efficiency can reach more than 70%, be far longer than the amount of water injected and oil displacement efficiency of practical oil reservoir；And in practical oil reservoir not
Water amount with position also differs, if it is assumed that the rock sample of rock core size is placed in oil reservoir, when the Injection volume of oil reservoir reaches
During 1PV-2PV, the water multiple of crossing of rock core scale also can reach up to ten thousand times near well, and sealing core data is also indicated that near well
More than 70% can be reached by waiting the oil displacement efficiency of regional areas.It can be seen that the mistake of oil reservoir different parts is described based on rock core scale
Water multiple, oil displacement efficiency can preferably reflect the practical displacement situation of oil reservoir.
Therefore, the mathematical method that oil reservoir different parts cross water multiple is described based on rock core scale we have invented a kind of, filled
That divides the abundant materials accurate description oil reservoir different parts for utilizing laboratory test crosses water multiple, quantization signifying oil reservoir different parts
Displacement degree establishes solid foundation to formulate mining site modification scenario technical measures.
Invention content
The object of the present invention is to provide a kind of researchs to match with coreflooding experiment, quantization signifying oil reservoir different parts
Displacement degree number crosses water times number calculating method based on different parts in rock core scale oil reservoir.
The purpose of the present invention can be achieved by the following technical measures：Water times is crossed based on different parts in rock core scale oil reservoir
Number calculating method should cross water times number calculating method based on different parts in rock core scale oil reservoir and include：Step 1, five-spot is derived
Uniform soft soil base streamline equation obtains different parts streamline distribution；Step 2, flow allocation method between streamline is established between injection-production well, is obtained
To different parts flow distribution；Step 3, it calculates and water times is crossed along streamline different parts based on the one-dimensional uniform soft soil base of rock core scale
Number；And step 4, calculate five-spot uniform soft soil base crosses water multiple along streamline different parts.
The purpose of the present invention can be also achieved by the following technical measures：
In step 1, according to Darcy's law, the principle of stacking of gesture, the gesture of coordinates regional arbitrary point is obtained；Stream function and gesture
Meet Cauchy-Riemann condition between function, stream function expression formula is obtained to get to five-spot homogeneous mould according to potential function expression formula
Type streamline equation.
In step 1, the five-spot uniform soft soil base streamline equation obtained is：
In formula：
The flow of q-producing well, m3/d；
X-abscissa, m；
Y-ordinate, m；
The 1/2, m of a-oil well well spacing；
C-it is the integral constant determined by boundary condition.
In step 2, the flow on every streamline is acquired according to resistor coupled in parallel principle, i.e., the streamline two between same injection-production well
End has identical injection production pressure difference；Distance between streamline is different, and water saturation distribution is different, and filtrational resistance is different；
Assuming that from the N streamline that have of well, wherein the filtrational resistance of kth streamline is RkIf QtFor practical well
Total water injection rate, the then water injection rate that kth streamline distributes are：
In step 3, it is assumed that model meshes are turned to n grid, the length of each grid by One-Dimensional Piston formula waterflooding model
It spends for dx, model total length, that is, injector producer distance is L, model width B, thickness h, porosity Φ, VφFor fenestral porosity body
Product, VoVolume for moveable oil, it is assumed that constant speed is injected, injection rate Q, water injection time t, swcFor irreducible water saturation；
Leading edge locus is obtained according to material balance：
All injection water can all flow through first grid, therefore the water multiple of crossing of first grid is exactly with total water injection rate
Divided by the pore volume of grid to obtain the final product：
The water that second grid flows through is inflow-rate of water turbine of total injection water first the oily displacement in first grid after complete, because
This, the water multiple excessively of second grid：
Similarly, the water multiple excessively of n-th of grid：
In step 4, between injection-production well in different zones by stream number of lines it is different；Grid streamline is whole where well point
Pass through, and some folding corner regions pass through even without streamline；The streamline density of different zones is different；
With grid where well by stream number of lines nmaxFor radix, define different grids by streamline density be：
PviFor the streamline by i-th of grid among i-th of mesh point it is corresponding based on rock core scale infinitesimal cross water
Multiple.
Water times number calculating method is crossed based on different parts in rock core scale oil reservoir in the present invention, by laboratory test and mining site
Development Practice organically combines, and quantitatively characterizing displacement degree develops heterogeneous influence to mining site, is mining site modification scenario strategy
Formulation provide strong foundation.This realizes Efficient Development oil field for the ultra-high water-containing stage, has important theory and reality
Meaning.
Description of the drawings
Fig. 1 is the specific embodiment that water times number calculating method is crossed based on different parts in rock core scale oil reservoir of the present invention
Flow chart；
Fig. 2 is the five-spot uniform soft soil base streamline distribution figure in the specific embodiment of the present invention；
Fig. 3 is that the difference in the specific embodiment of the present invention crosses filtrational resistance and well spacing relational graph under water multiple；
Fig. 4 is that the one-dimensional uniform soft soil base in the specific embodiment of the present invention crosses the distribution of water multiple along streamline different parts
Figure；
Fig. 5 crosses water multiple distribution map for the five-spot uniform soft soil base difference portion in the specific embodiment of the present invention.
Specific embodiment
For enable the present invention above and other objects, features and advantages be clearer and more comprehensible, it is cited below particularly go out preferable implementation
Example, and coordinate institute's accompanying drawings, it is described in detail below.
As shown in FIG. 1, FIG. 1 is the streams that water times number calculating method is crossed based on different parts in rock core scale oil reservoir of the present invention
Cheng Tu.
In a step 101, five-spot uniform soft soil base streamline equation is derived, obtains different parts streamline distribution.According to darcy
The principle of stacking of law, gesture can obtain the gesture of coordinates regional arbitrary point；Meet Cauchy-Riemann condition between stream function and potential function,
Stream function expression formula can obtain to get to Five-point method pattern streamline equation according to potential function expression formula：
In formula：
The flow of q-producing well, m3/d；
X-abscissa, m；
Y-ordinate, m；
The 1/2, m of a-oil well well spacing；
C-it is the integral constant determined by boundary condition.
According to streamline equation, you can obtain streamline distribution.
In one embodiment, according to formula (1) Five-point method pattern streamline equation, a five-spot uniform soft soil base is calculated
Streamline distribution figure, as shown in Figure 2.
In step 103, flow allocation method between streamline is established between injection-production well, obtains different parts flow distribution.According to
Resistor coupled in parallel principle acquires the flow on every streamline.That is, the streamline both ends between same injection-production well have identical injection production pressure difference；
Distance between streamline is different, and water saturation distribution is different, and filtrational resistance is different.
Assuming that from the N streamline that have of well, wherein the filtrational resistance of kth streamline is RkIf QtFor practical well
Total water injection rate, the then water injection rate that kth streamline distributes are：
In one embodiment, the filtrational resistance on every streamline is acquired according to resistor coupled in parallel principle, as shown in Figure 3.It is found that
Filtrational resistance is gradually reduced with the increase of amount of water injected, and incipient stage decline is obvious, and the later stage is due to the injection rate of water
Greatly, the relative permeability of oil phase reduces to a certain extent, and after the relative permeability rising to a certain extent of water phase, variation is little,
Filtrational resistance variation range very little, and good logarithmic relationship is presented with amount of water injected in filtrational resistance.It is acquired often according to formula (2)
The water injection rate of one streamline distribution.
In step 105, it calculates and water multiple is crossed along streamline different parts based on the one-dimensional uniform soft soil base of rock core scale.Assuming that
Model meshes are turned to n grid by One-Dimensional Piston formula waterflooding model, and the length of each grid is dx, and model total length is noted
Well spacing is adopted as L, model width B, thickness h, porosity Φ, VφFor fenestral porosity volume, VoIt is false for the volume of moveable oil
If constant speed is injected, injection rate Q, water injection time t, swcFor irreducible water saturation.
Leading edge locus is obtained according to material balance：
All injection water can all flow through first grid, therefore the water multiple of crossing of first grid is exactly with total water injection rate
Divided by the pore volume of grid to obtain the final product：
The water that second grid flows through is inflow-rate of water turbine of total injection water first the oily displacement in first grid after complete, because
This, the water multiple excessively of second grid：
Similarly, the water multiple excessively of n-th of grid：
In one embodiment, the master data of one-dimensional uniform soft soil base calculated examples, as shown in table 1.According to formula (3)~
(6), the one-dimensional uniform soft soil base based on rock core scale is calculated and crosses water multiple along streamline different parts, as shown in Figure 4.It is found that
Amount of water injected is bigger, and the water multiple excessively of different parts is bigger.Under same amount of water injected, with the increase with well distance, water is crossed
Multiple linearly reduces.And it crosses water multiple and the amount of water injected of model is mutually far short of what is expected, hundreds of can be reached even by crossing water multiple
It is thousands of, and macroscopical amount of water injected only has several times.
1 calculated examples data of table are set
| Parameter | Numerical value | Parameter | Numerical value |
| Length | 200m | Injection rate | 10% |
| Thickness | 3m | Initial oil saturation | 0.7 |
| Width | 10m | Residual oil saturation | 0.25 |
| Porosity | 0.25 | Rock core length | 5cm |
In step 107, calculate five-spot uniform soft soil base crosses water multiple along streamline different parts.Not same district between injection-production well
In domain by stream number of lines it is different；Grid streamline where well point all passes through, and some folding corner regions are even without streamline
Pass through.The streamline density of different zones is different.
With grid where well by stream number of lines nmaxFor radix, define different grids by streamline density be：
PviFor the streamline by i-th of grid among i-th of mesh point it is corresponding based on rock core scale infinitesimal cross water
Multiple.
In one embodiment, it is assumed that five-spot uniform soft soil base length is 200m, wherein the rock core length intercepted is 5cm, rock
According to formula (7)~(8), the mistake of any streamline different parts in different zones between injection-production well is calculated in core diameter position 2.5cm
Water multiple, as shown in Figure 5.It is found that near oil well, well cross water multiple can reach more than 100000PV, from water injection well to
Producing well has a main stream line, is nearby the flow channel of advantage, and Injection volume is more than 40000PV.Other overwhelming majority areas
The water multiple crossed in domain is less than 10000PV, and for part corner areas Injection volume less than 500PV, different parts cross that water multiple is differential to be reached
To more than 200.
Claims (5)
1. cross water times number calculating method based on different parts in rock core scale oil reservoir, which is characterized in that should be based on rock core scale oil
Different parts are crossed water times number calculating method and are included in Tibetan：
Step 1, five-spot uniform soft soil base streamline equation is derived, obtains different parts streamline distribution；
Step 2, flow allocation method between streamline is established between injection-production well, obtains different parts flow distribution；
Step 3, it calculates and water multiple is crossed along streamline different parts based on the one-dimensional uniform soft soil base of rock core scale；And
Step 4, calculate five-spot uniform soft soil base crosses water multiple along streamline different parts；
In step 2, the flow on every streamline is acquired according to resistor coupled in parallel principle, i.e., the streamline both ends tool between same injection-production well
There is identical injection production pressure difference；Distance between streamline is different, and water saturation distribution is different, and filtrational resistance is different；
Assuming that from the N streamline that have of well, wherein the filtrational resistance of kth streamline is RkIf QtTotal note for practical well
Water, the then water injection rate that kth streamline distributes are：
2. according to claim 1 cross water times number calculating method based on different parts in rock core scale oil reservoir, feature exists
In in step 3, it is assumed that model meshes are turned to n grid by One-Dimensional Piston formula waterflooding model, and the length of each grid is
Dx, model total length, that is, injector producer distance be L, model width B, thickness h, porosity Φ, VφFor fenestral porosity volume, Vo
Volume for moveable oil, it is assumed that constant speed is injected, injection rate Q, water injection time t, swcFor irreducible water saturation；
Leading edge locus is obtained according to material balance：
All injection water can all flow through first grid, thus first grid cross water multiple be exactly with total water injection rate divided by
The pore volume of grid to obtain the final product：
The water that second grid flows through is inflow-rate of water turbine of total injection water first the oily displacement in first grid after complete, therefore,
The water multiple excessively of second grid：
Similarly, the water multiple excessively of n-th of grid：
3. according to claim 1 cross water times number calculating method based on different parts in rock core scale oil reservoir, feature exists
In, in step 4, between injection-production well in different zones by stream number of lines it is different；Grid streamline where well point all passes through,
And some folding corner regions pass through even without streamline；The streamline density of different zones is different；
With grid where well by stream number of lines nmaxFor radix, define different grids by streamline density be：
Different parts cross the calculation formula (6) of water multiple on known each streamline, then obtain well pattern any position rock core scale
Infinitesimal cross water multiple be：
PviFor the streamline by i-th of grid among i-th of mesh point it is corresponding based on rock core scale infinitesimal cross water multiple.
4. according to claim 1 cross water times number calculating method based on different parts in rock core scale oil reservoir, feature exists
In in step 1, according to Darcy's law, the principle of stacking of gesture, obtaining the gesture of coordinates regional arbitrary point；Stream function and potential function
Between meet Cauchy-Riemann condition, stream function expression formula is obtained to get to five-spot uniform soft soil base stream according to potential function expression formula
Line equation.
5. according to claim 4 cross water times number calculating method based on different parts in rock core scale oil reservoir, feature exists
In in step 1, obtained five-spot uniform soft soil base streamline equation is：
In formula：
The flow of q-producing well, m3/d；
X-abscissa, m；
Y-ordinate, m；
The 1/2, m of a-oil well well spacing；
C-it is the integral constant determined by boundary condition.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410641714.0A CN105626007B (en) | 2014-11-07 | 2014-11-07 | Water times number calculating method is crossed based on different parts in rock core scale oil reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410641714.0A CN105626007B (en) | 2014-11-07 | 2014-11-07 | Water times number calculating method is crossed based on different parts in rock core scale oil reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105626007A CN105626007A (en) | 2016-06-01 |
| CN105626007B true CN105626007B (en) | 2018-06-15 |
Family
ID=56041298
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410641714.0A Active CN105626007B (en) | 2014-11-07 | 2014-11-07 | Water times number calculating method is crossed based on different parts in rock core scale oil reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105626007B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109559247B (en) * | 2017-09-27 | 2020-11-20 | 中国石油化工股份有限公司 | Method for determining injection speed of water drive experiment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101899972A (en) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | Establishment method of physical model capable of forecasting waterflood development of sandstone reservoirs |
| CN102146788A (en) * | 2011-03-30 | 2011-08-10 | 中国石油化工股份有限公司 | Method for enhancing recovery ratio through three-dimensional streamline regulation of water drive reservoir |
| CN103498669A (en) * | 2013-09-04 | 2014-01-08 | 中国石油天然气股份有限公司 | Quantitative determination method of interbedded cross flows of heterogeneous rock core models |
| CN103939066A (en) * | 2014-04-04 | 2014-07-23 | 中国石油大学(华东) | Method for determining oil well liquid production capacities by determining water injection capacities of one-injection and multi-production well group |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2945879B1 (en) * | 2009-05-20 | 2011-06-24 | Inst Francais Du Petrole | METHOD FOR OPERATING A POROUS MEDIUM USING MODELING FLUID FLOWS |
| US8646525B2 (en) * | 2010-05-26 | 2014-02-11 | Chevron U.S.A. Inc. | System and method for enhancing oil recovery from a subterranean reservoir |
-
2014
- 2014-11-07 CN CN201410641714.0A patent/CN105626007B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101899972A (en) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | Establishment method of physical model capable of forecasting waterflood development of sandstone reservoirs |
| CN102146788A (en) * | 2011-03-30 | 2011-08-10 | 中国石油化工股份有限公司 | Method for enhancing recovery ratio through three-dimensional streamline regulation of water drive reservoir |
| CN103498669A (en) * | 2013-09-04 | 2014-01-08 | 中国石油天然气股份有限公司 | Quantitative determination method of interbedded cross flows of heterogeneous rock core models |
| CN103939066A (en) * | 2014-04-04 | 2014-07-23 | 中国石油大学(华东) | Method for determining oil well liquid production capacities by determining water injection capacities of one-injection and multi-production well group |
Non-Patent Citations (1)
| Title |
|---|
| 水驱油藏中微观注水倍数计算方法;科学技术与工程;《科学技术与工程》;20140831;第14卷(第23期);第184页-第186页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105626007A (en) | 2016-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103246820B (en) | A kind of numerical simulation for oil-gas reservoir computing method | |
| CN103939066B (en) | An a kind of note adopts the method that well group determines water injection rate determination oilwell produced fluid amount more | |
| CN102913233B (en) | Method for recognizing dominant flow channel based on zero dimension comparison plate | |
| AU2012375233B2 (en) | System and method for automatic local grid refinement in reservoir simulation systems | |
| CN103455667B (en) | Aeration administers the method for numerical simulation of artesian aquifer seawater invasion | |
| CN104179499A (en) | Numerical simulation method considering oil reservoir parameter time variant | |
| CN105547967B (en) | Fissuted medium system permeability tensor indoor measurement device | |
| CN109522634A (en) | A kind of tight gas multistage volume fracturing horizontal well numerical analysis method | |
| CN107038268A (en) | A kind of method for determining heterogeneous reservoir five-spot pattern water drive sweep efficiency | |
| CN106649963B (en) | Volume fracturing complexity seam net average crack length and equivalent fissure item number determine method | |
| CN104331537A (en) | Well placement optimization design method based on reservoir static factors | |
| CN105626007B (en) | Water times number calculating method is crossed based on different parts in rock core scale oil reservoir | |
| CN107169227B (en) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well | |
| CN106481315B (en) | Land sandstone oil reservoir individual well recoverable reserves quickly determines model and method for building up | |
| CN103498661B (en) | A kind of method determining oil reservoir physical data under high pressure | |
| CN109632604B (en) | Method for coarsening relative permeability of polymer flooding from pore size to core size | |
| Cao et al. | Study of single phase mass transfer between matrix and fracture in tight oil reservoirs | |
| Peng et al. | A numerical simulation of seepage structure surface and its feasibility | |
| Artus* et al. | Efficient Proxies for Numerical Simulation of Unconventional Resources | |
| CN110991084B (en) | Reservoir permeability calculation method based on streamline numerical value well test | |
| CN109505591A (en) | Determine the method and system of the unfilled solution cavity permeability boundary of fracture-pore reservoir | |
| CN106529199B (en) | A kind of determination method of Conglomerate Reservoir chemical flooding well spacing | |
| CN105631763A (en) | Underground water type drinking water source protection area division method | |
| CN109839342A (en) | The multiple dimensioned drill core permeability rate determination techniques of high density electrical parameter monitoring | |
| Bedrikovetsky et al. | Mathematical Theory for Two Phase Displacement with Hysteresis (With application to WAG injection) |
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 |