CN107829718B - Oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory - Google Patents
Oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory Download PDFInfo
- Publication number
- CN107829718B CN107829718B CN201710074025.XA CN201710074025A CN107829718B CN 107829718 B CN107829718 B CN 107829718B CN 201710074025 A CN201710074025 A CN 201710074025A CN 107829718 B CN107829718 B CN 107829718B
- Authority
- CN
- China
- Prior art keywords
- well
- note
- injection
- optimization
- adopts
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000013461 design Methods 0.000 title claims abstract description 21
- 238000005457 optimization Methods 0.000 claims abstract description 72
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 31
- 238000006073 displacement reaction Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000012407 engineering method Methods 0.000 claims abstract description 9
- 238000002360 preparation method Methods 0.000 claims abstract description 4
- 238000002347 injection Methods 0.000 claims description 30
- 239000007924 injection Substances 0.000 claims description 30
- 238000009826 distribution Methods 0.000 claims description 25
- 230000035699 permeability Effects 0.000 claims description 20
- 238000004519 manufacturing process Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- 239000004519 grease Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000004576 sand Substances 0.000 claims description 7
- 238000009795 derivation Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 238000010845 search algorithm Methods 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims 1
- 229920006395 saturated elastomer Polymers 0.000 claims 1
- 125000006850 spacer group Chemical group 0.000 claims 1
- 238000003860 storage Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 10
- 238000011084 recovery Methods 0.000 abstract description 2
- 239000009671 shengli Substances 0.000 description 21
- 238000010586 diagram Methods 0.000 description 14
- 239000007788 liquid Substances 0.000 description 9
- 238000004088 simulation Methods 0.000 description 6
- 239000003129 oil well Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000012913 prioritisation Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 206010016256 fatigue Diseases 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 239000011550 stock solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/30—Specific pattern of wells, e.g. optimising the spacing of wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/20—Displacing by water
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory, comprising: step 1, collects and arrange block geology and exploitation related data；Step 2, setting optimization relevant parameter completes well pattern note and adopts optimization preparation；Step 3, using reservoir engineering method, predict that current well location and note adopt displacement situation and quantitative assessment in all directions under parameter；Step 4, using global random searching algorithm, optimization generates new well location/note and adopts parameter, predicts and evaluate new well location note to adopt the displacement situation respectively infused and adopted on direction under parameter；Step 5, optimum results are arranged, well pattern note is formed and adopts design scheme, put into field conduct.This method is combined by the waterflooding development index calculating method established based on reservoir engineering theories with Optimum Theory, use a kind of global random searching algorithm automatic calculation, it ensure that and obtain the well pattern to match with practical oil reservoir and injection-production program while computational efficiency, improve oil recovery.
Description
Technical field
The present invention relates to oil-gas field development field, especially relate to a kind of oil reservoir well pattern based on balanced water drive theory and
Injection-production program optimum design method.
Background technique
Well pattern and injection-production program are the important contents in oil reservoir development scheme, directly influence effect of reservoir development.How
Making optimal well pattern and injection-production program is one of the problem in current oilfield development program design.Currently used well pattern and
Injection-production program optimum design method mainly has: carrying out well pattern and injection-production program design by artificial experience and is based on Optimum Theory
Well pattern and injection-production program optimization design.Well pattern and injection-production program design are carried out by artificial experience, needs Consideration numerous,
Subjective factor influences greatly, while being difficult to obtain real optimal well pattern and injection-production program；Well pattern and note based on Optimum Theory
Scheme optimization is adopted, using the superiority and inferiority of each scheme in numerical simulation evaluation optimization process, for practical extensive oil reservoir, calculation amount
Greatly, time-consuming, and conventional computer is difficult to realize.And method is broadly divided into well net optimization method at present and note adopts optimization method two
Mutually independent part, lacks combination.When being solved using gradient algorithm, gradient seeks complexity, and Restriction condition treat is complicated, together
When be easy to fall into local optimum and real optimal case can not be found.We have invented a kind of thus based on balanced water drive theory
Oil reservoir well pattern and injection-production program optimum design method, solve the above technical problem.
Summary of the invention
The object of the present invention is to provide obtain the well pattern and note to match with practical oil reservoir while a kind of computational efficiency to adopt
Scheme improves the oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory of oil recovery.
The purpose of the present invention can be achieved by the following technical measures: oil reservoir well pattern and note based on balanced water drive theory are adopted
Scheme optimization design method, which includes: step
1, collect and arrange block geology and exploitation related data；Step 2, setting optimization relevant parameter, it is quasi- that completion well pattern note adopts optimization
Standby work；Step 3, using reservoir engineering method, predict that current well location and note adopt under parameter in all directions displacement situation and quantitative
Evaluation；Step 4, using global random searching algorithm, optimization generates new well location note and adopts parameter, predicts and evaluates new well location/note
Adopt the displacement situation respectively infused and adopted on direction under parameter；Step 5, optimum results are arranged, well pattern note is formed and adopts design scheme, investment is existing
Implement field.
The purpose of the present invention can be also achieved by the following technical measures:
In step 1, the data of collection includes: static data: Permeability Distribution, porosity distribution, sand thickness, net hair
Than distribution；Dynamic data: saturation distribution, pressure distribution, grease density, grease viscosity, permeability saturation curve, existing oil
Well well location.
In step 2, setting optimization relevant parameter includes: setting optimization type, and optimization type includes well net optimization, infuses and adopt
Optimization, well pattern note adopt combined optimization；Well to be optimized is specified, optimized variable is generated；The constraint condition for specifying each optimized variable, for
Well net optimization, constraint condition include section locating for section, Y-coordinate locating for the X-coordinate of well, reservoir boundary constraint；Note is adopted excellent
Change, constraint condition includes the overall note amount of adopting, individual well infuses the amount of the adopting upper limit, individual well infuses the amount of adopting lower limit；Set initial well pattern/note side of adopting
Case；Set derivation algorithm relevant parameter, including initial ranging step-length, initial sample number, termination condition.
The oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory further include, after step 2,
Initial well location/the note for generating well to be optimized adopts parameter, comprising: judges whether is the initial well pattern/injection-production program of setting in step 2
It executes；If executing, using its input scheme as initial scheme；If being not carried out, set according to the constraint condition in step 2,
Each optimized variable is generated at random, as initial scheme.
In step 3, according to oil-water well well location, injection-production relation matrix is generated；According to the angle of two groups of injection-production well lines
Oil reservoir is splitted into multiple notes and adopts control unit by bisector；Each note adopts control unit intrinsic parameter equivalent process, and parameter includes infiltration
Rate, porosity, net-gross ratio, sand thickness, saturation degree, equivalent way are weighted average；Use theoretical calculation of reservoir engineering
Current time walks each equivalent filtrational resistance of control unit；Each control unit is walked using theoretical calculation of reservoir engineering current time to inject
Speed；Calculate each control unit average staturation after current time walks；Renewal time, t=t+ Δ t, wherein t is current time,
Δ t is time step, and it is average to compute repeatedly the equivalent filtrational resistance of each control unit, each each control unit of control unit injection rate
These steps of saturation degree, until reaching the prediction end time；Calculate the standard deviation of each flooding unit average staturation.
In step 3, current time is calculated using following formula walk each equivalent filtrational resistance of control unit:
In formula: RiFor the filtrational resistance between water injection well and i-th mouthful of producing well, mPas/ (μm2·cm)；I is producing well volume
Number；For the mean permeability between water injection well and i-th mouthful of producing well, 10-3μm2；It is flat between water injection well and i-th mouthful of producing well
Equal reservoir thickness, m；For the angle between two adjacent groups injection-production well line and the angular bisector of the injection-production well line；rfTo drive
For leading edge distance, m；rwFor the radius of water injection well, m；KroFor oil relative permeability；μoFor oil viscosity, mPas；Krw
For water phase relative permeability；μwFor water flooding viscosity, mPas；SwcFor irreducible water saturation；diIt is produced for water injection well and i-th mouthful
The well spacing of well, m；SweFor exit-end water saturation.
In step 3, current time is calculated using following formula walk each control unit injection rate:
qi=Δ pi/Ri
In formula: Δ piFor the pressure difference between water injection well and i-th mouthful of producing well, MPa；qiFor flooding unit where i-th mouthful of producing well
Injection rate, m3/d。
In step 3, each control unit average staturation after current time walks is calculated:
In formula: r is displacement distance, m；SwFor water saturation；T is displacement time, d；fwFor moisture content.
In step 4, the standard deviation of saturation degree under each scheme is calculated；It sorts from small to large by standard deviation；According to saturation degree
The smallest scheme of standard deviation generates well pattern/note side of adopting of new generation by the filial generation generation strategy in global random searching algorithm
Case, global random searching algorithm include genetic algorithm, particle swarm algorithm, covariance matrix evolution algorithm；Judge newly-generated well
Whether net/injection-production program meets constraint condition, recalculates if not satisfied, going to previous step；Judge newly-generated well pattern/note
Adopt whether scheme reaches optimization termination condition, when meeting optimization termination condition, optimization terminates, and otherwise returns and calculates current iteration
Walk the standard deviation of average staturation under each scheme.
The oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory in the present invention, based on optimization
Theoretical and automatic technology, the influence of the subjective factor of people, mitigates the working strength of people when avoiding solution formulation；This method is based on
Balanced waterflooding development theory carries out scheme evaluation using reservoir engineering method substitution numerical simulation, avoids in numerical simulation
Large-scale matrix iterative calculation, shortens scheme evaluation time, such that extensive oil reservoir optimizes；This method uses global
Random search algorithm seeks optimal case, can easily handle all kinds of constraint conditions, and the related mechanism in searching algorithm can be protected
Card obtains globally optimal solution；This method optimizes note on the basis of optimized well pattern and adopts parameter, and then obtains the optimal well pattern note side of adopting
It is mutually coordinated between parameter to ensure that well pattern and note are adopted for case.Each factor of oil reservoir can be comprehensively considered by this method, automatically
Optimal well pattern and injection-production program are obtained, technical support is provided for oil field development, improves waterflooding development effect.
Detailed description of the invention
Fig. 1 is that the present invention is based on the oil reservoir well patterns of balanced water drive theory and note to adopt optimum design method flow diagram；
Fig. 2 is Optimal Parameters setting procedure flow diagram；
Fig. 3 is that reservoir engineering calculates development index and evaluates displacement balance degree flow diagram；
Fig. 4 adopts control unit for note and splits a point schematic diagram；
Fig. 5 is that Optimized model solves flow diagram；
Fig. 6 is Shengli Oil Field A block part data graph；
Fig. 7 is Shengli Oil Field A block well well location restriction range schematic diagram to be optimized；
Fig. 8 is the initial well pattern schematic diagram of Shengli Oil Field A Block Set；
Fig. 9 is that method for numerical simulation and reservoir engineering method calculate time comparison diagram under different scales oil reservoir；
Interative computation is respectively for minimum saturation standard deviation figure when Figure 10 is Shengli Oil Field A block well net optimization；
Interative computation is respectively for minimum saturation standard deviation figure when Figure 11 adopts optimization for Shengli Oil Field A block note；
Figure 12 is well pattern arrangement schematic diagram after the optimization of Shengli Oil Field A block well location；
Figure 13 is that Shengli Oil Field A block oil well note adopts optimization front and back liquid measure comparison diagram；
Figure 14 is that Shengli Oil Field A block well note adopts optimization front and back liquid measure comparison diagram；
Figure 15 is Shengli Oil Field A block well pattern and note adopts prioritization scheme and oil-producing comparison diagram is tired out in artificial program prediction.
Specific embodiment
To enable above and other objects, features and advantages of the invention to be clearer and more comprehensible, embodiment is cited below particularly out, and
Cooperate institute's accompanying drawings, is described in detail below.
As shown in FIG. 1, FIG. 1 is the oil reservoir well pattern of the invention based on balanced water drive theory and injection-production program optimization design sides
The flow chart of one specific embodiment of method.
In step 101, block geology and exploitation related data are collected and arranged.The data of collection includes: static data: being seeped
Saturating rate distribution, porosity distribution, sand thickness, net-gross ratio distribution；Dynamic data: saturation distribution, pressure distribution, grease are close
Degree, grease viscosity, permeability saturation curve, existing oil-water well well location.Process enters step 102.
In step 102, setting optimization relevant parameter completes well pattern/note and adopts optimization preparation.As shown in Fig. 2, setting is excellent
Changing relevant parameter includes: setting optimization type, and optimization type includes well net optimization, note adopts optimization, well pattern note adopts combined optimization；Refer to
Fixed well to be optimized, generates optimized variable；The constraint condition for specifying each optimized variable, for well net optimization, constraint condition includes well
X-coordinate locating for section, section locating for Y-coordinate, reservoir boundary constraint；Optimization is adopted for note, constraint condition includes that overall note is adopted
Amount, the individual well note amount of the adopting upper limit, individual well infuse the amount of adopting lower limit；Set initial well pattern/injection-production program (optional)；It is related to set derivation algorithm
Parameter, including initial ranging step-length, initial sample number, termination condition.Process enters step 103.
In step 103, the initial well location/note for generating well to be optimized adopts parameter.The setting initial well of judgement in a step 102
Whether net/injection-production program (optional) executes；If executing, using its input scheme as initial scheme；If being not carried out, according to step
Constraint condition setting in rapid 102, generates each optimized variable, as initial scheme at random.Process enters step 104.
In step 104, using reservoir engineering method, predict that current well location and note adopt under parameter in all directions displacement situation simultaneously
Quantitative assessment.As shown in figure 3, generating injection-production relation matrix according to oil-water well well location；As shown in figure 4, being adopted according to two groups of notes
Oil reservoir is splitted into multiple notes and adopts control unit by the angular bisector of well line；Each note adopts control unit intrinsic parameter equivalent process, joins
Number includes permeability, porosity, net-gross ratio, sand thickness, saturation degree, and equivalent way is weighted average；Use oil reservoir work
Journey theoretical calculation current time walks each equivalent filtrational resistance of control unit；Each control is walked using theoretical calculation of reservoir engineering current time
Unit injection rate processed；Calculate each control unit average staturation after current time walks；Renewal time, t=t+ Δ t, wherein t be
Current time, Δ t are time step, compute repeatedly the equivalent filtrational resistance of each control unit, each control unit injection rate is respectively controlled
These steps of cell-average saturation degree processed, until reaching the prediction end time；Calculate the mark of each flooding unit average staturation
It is quasi- poor.
Current time, which is calculated, using following formula walks each equivalent filtrational resistance of control unit:
In formula: RiFor the filtrational resistance between water injection well and i-th mouthful of producing well, mPas/ (μm2·cm)；I is producing well volume
Number；For the mean permeability between water injection well and i-th mouthful of producing well, 10-3μm2；Between water injection well and i-th mouthful of producing well
Average reservoir thickness, m；For the angle between two adjacent groups injection-production well line and the angular bisector of the injection-production well line；rfFor
Displacing front distance, m；rwFor the radius of water injection well, m；KroFor oil relative permeability；μoFor oil viscosity, mPas；
KrwFor water phase relative permeability；μwFor water flooding viscosity, mPas；SwcFor irreducible water saturation；diFor water injection well and i-th mouthful of life
Produce the well spacing of well, m；SweFor exit-end water saturation.
Current time, which is calculated, using following formula walks each control unit injection rate:
qi=Δ pi/Ri
In formula: Δ piFor the pressure difference between water injection well and i-th mouthful of producing well, MPa；qiFor flooding unit where i-th mouthful of producing well
Injection rate, m3/d。
Calculate each control unit average staturation after current time walks:
In formula: r is displacement distance, m；SwFor water saturation；T is displacement time, d；fwFor moisture content.
Process enters step 105.
In step 105, judge whether to reach optimization termination condition, if not satisfied, process enters step 106；If satisfied, stream
Journey enters step 107；
In step 106, as shown in figure 5, optimization generates new well location/note and adopts parameter, in advance using global random searching algorithm
It surveys and evaluates new well location/note and adopt the displacement situation respectively infused and adopted on direction under parameter.Calculate the standard deviation of saturation degree under each scheme；It presses
Standard deviation sorts from small to large；It is raw by the filial generation in global random searching algorithm according to the smallest scheme of saturation degree standard deviation
At strategy, well pattern/injection-production program of new generation is generated, global random searching algorithm includes but is not limited to genetic algorithm, population calculation
Method, covariance matrix evolution algorithm；Judge whether newly-generated well pattern/injection-production program meets constraint condition, judges newly-generated
Whether well pattern/injection-production program reaches optimization termination condition, and when meeting optimization termination condition, optimization terminates, and otherwise returns and calculates
Current iteration walks the standard deviation of average staturation under each scheme.Process enters step 104.
In step 107, optimum results are arranged, well pattern/note is formed and adopts design scheme, put into field conduct.Process terminates.
To enable above content of the invention to be clearer and more comprehensible, below by taking Shengli Oil Field A block as an example, using based on equilibrium
The oil reservoir well pattern and injection-production program optimum design method of waterflooding development theory, are described in detail below:
1, collect and arrange the block geology and exploitation related data:
Data source can provide in derived digitlization there are many form from existing geological model or numerical model
Material can also carry out digitized processing according to related geologic development graph and obtain.
The data of collection includes:
(1) static data: Permeability Distribution, porosity distribution, sand thickness, net-gross ratio distribution；
(2) dynamic data: saturation distribution, grease density, grease viscosity, permeability saturation curve, has been deposited at pressure distribution
In oil-water well well location.
Shengli Oil Field A block Permeability Distribution, porosity distribution, saturation distribution and existing well well location as shown in fig. 6,
Wherein, Fig. 6 a is Permeability Distribution figure, and Fig. 6 b is porosity distribution map, and Fig. 6 c is saturation distribution figure and existing well well location.
2, setting optimization relevant parameter completes well pattern/note and adopts optimization preparation.Its specific implementation step is as follows:
(1) setting optimization type:
Shengli Oil Field A block should carry out well net optimization, carry out note again and adopt optimization, therefore optimizing type set is well pattern
Note adopts combined optimization；
(2) well to be optimized is specified, optimized variable is generated:
After preliminary analysis, determination needs 4 mouthfuls of drilling new well, including 2 mouthfuls of new well, 2 mouthfuls of grease hole.Well location optimizes to excellent
Change well and be appointed as newly bore 4 mouthfuls of wells, liquid measure optimizes well to be optimized and is appointed as block whole oil-water well.
(3) each optimized variable constraint condition is specified:
In conjunction with oil reservoir reality, specify the well spacing range of 4 mouthfuls of wells in well net optimization as shown in Figure 7.Note adopts each well in optimization
Injection rate lower limit be 0, the upper limit is 2 times of current maximum individual-well injection rate；The produced quantity lower limit of each oil well is 0, and the upper limit is to work as
2 times of preceding maximum individual well produced quantity.
(4) initial well pattern/injection-production program is set:
Initial scheme can not be set, and for Shengli Oil Field A block, set initial well location such as Fig. 8 based on previous work.If
The fixed old initial liquid measure of well is that the amount of the stock solution is constant, sets the initial liquid measure of new well and is averaged liquid measure as block.
(5) derivation algorithm relevant parameter is set:
For Shengli Oil Field A block, the particle swarm algorithm in global random searching algorithm is selected, sets initial population number
50 generations, other parameters take algorithm default value itself.
3, the initial well location/note for generating well to be optimized adopts parameter:
Shengli Oil Field A block waits for that well location optimizes and amounts to 4 mouthfuls of wells, and each well includes two X-coordinate, Y-coordinate variables, amounts to 8
A optimized variable；The total 13 mouthfuls of wells of optimization are adopted wait infuse, optimized variable is liquid measure, amounts to 13 variables.It is set according in step 2- (4)
Fixed initial scheme carries out initialization assignment to each variable.
4, using reservoir engineering method, predict that current well location and note are adopted under parameter displacement situation in all directions and quantitatively commented
Valence:
Shengli Oil Field A block includes 13 mouthfuls of wells, and injection-production relation matrix is the 0-1 matrix of 13*13, wherein 0 value represents
Without corresponding relationship between two mouthfuls of wells, 1 value, which represents between two mouthfuls of wells, corresponding relationship.Predicted time is set as 15 years, iteration time
Step-length is 30 days, solves each note and adopts average staturation on line, and takes standard deviation.Fig. 9 is reservoir engineering method and numerical simulation
Method predicts operation time under different scales oil reservoir, it can be seen that reservoir engineering method efficiency is significantly higher than numerical simulation side
Method.The saturation degree variance that initial scheme is calculated is shown in Table 1.
The saturation degree variance table of 1 Shengli Oil Field A block initial scheme of table
5, using global random searching algorithm, optimization generates new well location/note and adopts parameter, predicts and evaluate new well location/note
Adopt the displacement situation respectively infused and adopted on direction under parameter:
(1) raw by the filial generation generation strategy in global random searching algorithm according to the smallest scheme of saturation degree standard deviation
At well pattern/injection-production program of new generation:
The global random searching algorithm that Shengli Oil Field A block is chosen is particle swarm algorithm, and population invariable number 50 generates filial generation
When population, the smallest individual of saturation degree variance can be retained, according to movement speed, Xiang Bendai optimum individual leans on other 49 individuals
Hold together, generate certain displacement, to generate offspring individual, and forms progeny population.
(2) judge whether newly-generated well pattern/injection-production program meets constraint condition, if not satisfied, going to step 5- (1) weight
It is new to calculate, until meeting constraint condition；
(3) it calculates current iteration and walks average staturation standard deviation under each scheme.
6, interative computation, until seeking obtaining optimal case, the specific steps are as follows:
(1) judge whether to reach optimization termination condition, if not satisfied, repeating step 5-6, terminate item until reaching optimization
Part；
(2) optimum results are arranged, well pattern/note is formed and adopts design scheme, put into field conduct.
Interative computation is respectively for minimum saturation standard deviation figure such as Figure 10 when Shengli Oil Field A block well location optimizes；Shengli Oil Field A
Interative computation is respectively for minimum saturation standard deviation figure such as Figure 11 when block note adopts optimization；After optimization, the optimal well spacing of formation
Scheme such as Figure 12；The specific data of well location are shown in Table 2 after optimization；Optimization front and back each oil well liquid measure comparison diagram such as Figure 13；Optimization front and back is each
Well liquid measure comparison diagram such as Figure 14；Compared to the scheme by engineer, prioritization scheme prediction can increase oily 1.87 ten thousand steres more,
Amplification 37.54%, significant effect, as shown in figure 15.
Well location table after the optimization of 2 Shengli Oil Field A block of table
Claims (1)
1. oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory, which is characterized in that should be based on equilibrium
The oil reservoir well pattern and injection-production program optimum design method of water drive theory include:
Step 1, collect and arrange block geology and exploitation related data；
Step 2, setting optimization relevant parameter completes well pattern note and adopts optimization preparation；
Step 3, using reservoir engineering method, predict that current well location and note are adopted under parameter displacement situation in all directions and quantitatively commented
Valence；
Step 4, using global random searching algorithm, optimization generates new well location note and adopts parameter, predicts and evaluate new well location/note to adopt
Respectively note adopts the displacement situation on direction under parameter；
Step 5, optimum results are arranged, well pattern note is formed and adopts design scheme, put into field conduct；
In step 1, the data of collection includes: static data: Permeability Distribution, porosity distribution, sand thickness, net-gross ratio point
Cloth；Dynamic data: saturation distribution, pressure distribution, grease density, grease viscosity, permeability saturation curve, existing oil-water well
Well location；
In step 2, setting optimization relevant parameter include: setting optimization type, optimization type include well net optimization, note adopt it is excellent
Change, well pattern note adopts combined optimization；Well to be optimized is specified, optimized variable is generated；The constraint condition for specifying each optimized variable, for well
Network optimization, constraint condition include section locating for section, Y-coordinate locating for the X-coordinate of well, reservoir boundary constraint；Optimization is adopted for note,
Constraint condition includes the overall note amount of adopting, the individual well note amount of the adopting upper limit, the individual well note amount of adopting lower limit；Set initial well pattern/injection-production program；If
Determine derivation algorithm relevant parameter, including initial ranging step-length, initial sample number, termination condition；After step 2, it generates to excellent
Initial well location/the note for changing well adopts parameter, comprising: judges whether the initial well pattern/injection-production program of setting in step 2 executes；If
It executes, then using its input scheme as initial scheme；If being not carried out, set according to the constraint condition in step 2, it is random to generate
Each optimized variable, as initial scheme；
In step 3, according to oil-water well well location, injection-production relation matrix is generated；According to the angle bisection of two groups of injection-production well lines
Oil reservoir is splitted into multiple notes and adopts control unit by line；Each note adopts control unit intrinsic parameter equivalent process, and parameter includes permeability, hole
Porosity, net-gross ratio, sand thickness, saturation degree, equivalent way are weighted average；When current using theoretical calculation of reservoir engineering
The equivalent filtrational resistance of each control unit of spacer step；Each control unit injection rate is walked using theoretical calculation of reservoir engineering current time；
Calculate each control unit average staturation after current time walks；Renewal time, t=t+ Δ t, wherein t is current time, and Δ t is
Time step computes repeatedly the equivalent filtrational resistance of each control unit, each each control unit of control unit injection rate is averagely saturated
These steps are spent, until reaching the prediction end time；Calculate the standard deviation of each flooding unit average staturation；
Current time, which is calculated, using following formula walks each equivalent filtrational resistance of control unit:
In formula: RiFor the filtrational resistance between water injection well and i-th mouthful of producing well, mPas/ (μm2·cm)；I is producing well number；
For the mean permeability between water injection well and i-th mouthful of producing well, 10-3μm2；For the average storage between water injection well and i-th mouthful of producing well
Thickness degree, m；For the angle between two adjacent groups injection-production well line and the angular bisector of the injection-production well line；rfBefore displacement
Edge distance, m；rwFor the radius of water injection well, m；KroFor oil relative permeability；μoFor oil viscosity, mPas；KrwFor water
Phase relative permeability；μwFor water flooding viscosity, mPas；SwcFor irreducible water saturation；diFor water injection well and i-th mouthful of producing well
Well spacing, m；SweFor exit-end water saturation；
Current time, which is calculated, using following formula walks each control unit injection rate:
qi=Δ pi/Ri
In formula: Δ piFor the pressure difference between water injection well and i-th mouthful of producing well, MPa；qiFor the note of flooding unit where i-th mouthful of producing well
Enter speed, m3/d；
Current time, which is calculated, using following formula walks each control unit average staturation:
In formula: r is displacement distance, m；SwFor water saturation；T is displacement time, d；fwFor moisture content；
In step 4, the standard deviation of saturation degree under each scheme is calculated；It sorts from small to large by standard deviation；According to saturation degree standard
The smallest scheme of difference generates well pattern/injection-production program of new generation, entirely by the filial generation generation strategy in global random searching algorithm
Office's random search algorithm includes genetic algorithm, particle swarm algorithm, covariance matrix evolution algorithm；Judge newly-generated well pattern/note
It adopts whether scheme meets constraint condition, is recalculated if not satisfied, going to previous step；Judge newly-generated well pattern/note side of adopting
Whether case reaches optimization termination condition, and when meeting optimization termination condition, optimization terminates, and it is each otherwise to return to calculating current iteration step
The standard deviation of average staturation under scheme.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710074025.XA CN107829718B (en) | 2017-02-10 | 2017-02-10 | Oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710074025.XA CN107829718B (en) | 2017-02-10 | 2017-02-10 | Oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107829718A CN107829718A (en) | 2018-03-23 |
| CN107829718B true CN107829718B (en) | 2019-11-08 |
Family
ID=61643063
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710074025.XA Active CN107829718B (en) | 2017-02-10 | 2017-02-10 | Oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107829718B (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109165783A (en) * | 2018-08-15 | 2019-01-08 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | Oilfield development program optimization method and device |
| CN109146185B (en) * | 2018-08-29 | 2022-02-22 | 中国石油化工股份有限公司 | Injection-production well pattern differential water flooding optimization adjustment method |
| CN109816182B (en) * | 2019-03-27 | 2023-06-20 | 常州大学 | Progressive splitting optimization method for injection and production dynamic regulation and control of water-flooding oil reservoir |
| CN110941890B (en) * | 2019-09-27 | 2022-11-01 | 中国海洋石油集团有限公司 | Offshore oil reservoir dynamic real-time production optimization method based on optimal control theory |
| CN110924908B (en) * | 2019-11-08 | 2021-10-08 | 中国石油大学(华东) | Water drive reservoir injection-production parameter determination method and computer readable storage medium |
| CN111119815B (en) * | 2020-01-03 | 2022-07-26 | 中国石油化工股份有限公司 | Method for determining production well production allocation ratio through balanced displacement |
| CN113137212A (en) * | 2020-01-17 | 2021-07-20 | 中国石油天然气股份有限公司 | Multi-layer heterogeneous oil reservoir water injection method |
| CN111222252B (en) * | 2020-01-20 | 2023-05-12 | 东北石油大学 | Method and system for predicting oil-water two-phase post-pressure capacity of low-saturation oil reservoir |
| CN111291981B (en) * | 2020-01-21 | 2023-08-25 | 常州大学 | Oil reservoir well pattern injection and production well group perfection evaluation method |
| CN113898331A (en) * | 2020-06-22 | 2022-01-07 | 中国石油化工股份有限公司 | Well spacing optimization method for oil reservoir balanced displacement differentiation in high water cut period |
| CN112070269B (en) * | 2020-07-31 | 2024-05-14 | 中海油能源发展股份有限公司 | Injection and production balance adjustment method based on depth traversal |
| CN114165227B (en) * | 2020-09-11 | 2024-03-01 | 中国石油化工股份有限公司 | Quantitative evaluation method for flow field adjustment feasibility based on oil reservoir balanced displacement |
| CN112360411B (en) * | 2020-11-30 | 2022-02-01 | 河海大学 | Local well pattern water injection development optimization method based on graph neural network |
| CN115906675B (en) * | 2023-03-09 | 2023-05-26 | 中国石油大学(华东) | Well position and injection and production parameter joint optimization method based on time sequence multi-target prediction model |
| CN117291125B (en) * | 2023-11-24 | 2024-02-02 | 中国石油大学(华东) | Comprehensive adjustment optimization method and system for water-flooding well pattern and perforation interval |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606120A (en) * | 2012-01-12 | 2012-07-25 | 倪红梅 | Steam flooding injection and production scheme determining method based on stochastic disturbance particle swarm optimization |
| CN103924966A (en) * | 2013-01-10 | 2014-07-16 | 中国石油化工股份有限公司 | Bed set well pattern reconstruction method based on reservoir physical property time varying model |
| CN104632157A (en) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | Low permeability reservoir equilibrium displacement method |
| CN105756654A (en) * | 2015-12-30 | 2016-07-13 | 中国石油大学(华东) | Water-drive reservoir well network optimization design method |
| CN106126899A (en) * | 2016-06-20 | 2016-11-16 | 中国石油大学(华东) | Water-drive pool Encryption Well well location determines method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AR102540A1 (en) * | 2014-11-07 | 2017-03-08 | Bp Corp North America Inc | METHODS TO HANDLE THE REPLACEMENT OF VACUUM TRAINING IN WATER INJECTION PRODUCTION OPERATIONS TO INCREASE OIL RECOVERY |
-
2017
- 2017-02-10 CN CN201710074025.XA patent/CN107829718B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606120A (en) * | 2012-01-12 | 2012-07-25 | 倪红梅 | Steam flooding injection and production scheme determining method based on stochastic disturbance particle swarm optimization |
| CN103924966A (en) * | 2013-01-10 | 2014-07-16 | 中国石油化工股份有限公司 | Bed set well pattern reconstruction method based on reservoir physical property time varying model |
| CN104632157A (en) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | Low permeability reservoir equilibrium displacement method |
| CN105756654A (en) * | 2015-12-30 | 2016-07-13 | 中国石油大学(华东) | Water-drive reservoir well network optimization design method |
| CN106126899A (en) * | 2016-06-20 | 2016-11-16 | 中国石油大学(华东) | Water-drive pool Encryption Well well location determines method |
Non-Patent Citations (1)
| Title |
|---|
| 窄条带砂体随机模拟井网优化部署方法;李莉等;《大庆石油地质与开发》;20001031;第19卷(第5期);第15-19页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107829718A (en) | 2018-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107829718B (en) | Oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory | |
| CN105095986B (en) | The method of stratified reservoir overall yield prediction | |
| CN105808311B (en) | A kind of reservoir modeling Fast Fitting method based on Dimension-reduced Strategy | |
| CN102041995B (en) | System for monitoring complicated oil deposit flooding conditions | |
| Wang et al. | A multilevel coordinate search algorithm for well placement, control and joint optimization | |
| CN104750896B (en) | A kind of fractured-cavernous carbonate reservoir method for numerical simulation | |
| RU2614338C1 (en) | Method of real-time control of reservoir flooding | |
| CN106294282B (en) | Black oil reservoir simulation method and device | |
| CN108603402A (en) | The variation of capillary pressure and relative permeability in porous media caused by mineral deposits and dissolving is modeled and predicted | |
| WO2011149779A1 (en) | System and method for enhancing oil recovery from a subterranean reservoir | |
| CN104453834A (en) | Injection-production relation optimizing and adjusting method for well group | |
| CN106529667A (en) | Logging facies identification and analysis method based on fuzzy depth learning in big data environment | |
| RU2715593C1 (en) | Method of operative control of water flooding of formations | |
| CN109386272A (en) | Ultra deep reef flat facies gas reservoir rational spacing between wells Multipurpose Optimal Method | |
| Zhang et al. | Surrogate-assisted multiobjective optimization of a hydraulically fractured well in a naturally fractured shale reservoir with geological uncertainty | |
| Khalid et al. | Optimal extraction of groundwater in Gaza coastal aquifer | |
| CN108729912A (en) | Production split method suitable for reservoir numerical simulation | |
| CN107169227B (en) | A kind of the coarse grid analogy method and system of staged fracturing horizontal well | |
| CN110894788B (en) | Oil reservoir production prediction data-based oil reservoir development mode determination method and device | |
| CN108133286B (en) | Underground water multi-target calculation method based on ground settlement substitution model | |
| CN106846470A (en) | A kind of high accuracy oil-gas migration analogy method based on Corner-point Grids | |
| RU2669322C1 (en) | Method for determining perspective sites and optimum technological parameters for use of cyclic water flooding | |
| Feng et al. | A new method for evaluating the effectiveness of hydraulic fracturing in tight reservoirs | |
| Leemput et al. | Full-field reservoir modeling of Central Oman gas-condensate fields | |
| RU2558093C1 (en) | Control method of oil pool development |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Huang Yingsong Inventor after: Liu Weiwei Inventor after: Wang Bo Inventor after: Guo Longfei Inventor after: Wang Duanping Inventor after: Lai Shumin Inventor after: Feng Qihong Inventor after: Zhang Yigen Inventor after: Liu Zhihong Inventor after: Wei Ming Inventor after: Wang Xiang Inventor after: Yang Shengbo Inventor before: Wang Duanping Inventor before: Wang Bo Inventor before: Lai Shumin Inventor before: Yang Shengbo Inventor before: Feng Qihong Inventor before: Zhang Yigen Inventor before: Wang Xiang Inventor before: Liu Zhihong Inventor before: Wei Ming Inventor before: Huang Yingsong Inventor before: Guo Longfei Inventor before: Liu Weiwei |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |