CN106126899B - Method for determining well positions of infilled wells of water-drive reservoir - Google Patents
Method for determining well positions of infilled wells of water-drive reservoir Download PDFInfo
- Publication number
- CN106126899B CN106126899B CN201610446421.6A CN201610446421A CN106126899B CN 106126899 B CN106126899 B CN 106126899B CN 201610446421 A CN201610446421 A CN 201610446421A CN 106126899 B CN106126899 B CN 106126899B
- Authority
- CN
- China
- Prior art keywords
- well
- note
- encryption
- location
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 63
- 238000006073 displacement reaction Methods 0.000 claims abstract description 41
- 230000035699 permeability Effects 0.000 claims abstract description 25
- 238000004519 manufacturing process Methods 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 6
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- 238000009826 distribution Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000012417 linear regression Methods 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims description 2
- 239000012530 fluid Substances 0.000 claims description 2
- 238000011835 investigation Methods 0.000 claims description 2
- 238000005325 percolation Methods 0.000 claims description 2
- 238000002347 injection Methods 0.000 abstract description 12
- 239000007924 injection Substances 0.000 abstract description 12
- 238000011161 development Methods 0.000 abstract description 11
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 2
- 230000002194 synthesizing effect Effects 0.000 abstract 1
- 230000008719 thickening Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000013475 authorization Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- 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 invention relates to a method for determining well positions of infilled wells of a water-drive reservoir. The method comprises the following steps of collecting and arranging geological data of the reservoir, a current development status, a region required to be infilled, a well number and well types of infilled wells, and a flowing bottomhole pressure of an injection and production well; giving an initial infilled well position, and performing equivalent processing on the permeability, porosity and initial water saturation, in each injection and production connection line direction, of the infilled well position; establishing a relationship between the water saturation in the injection and production direction and a well spacing, and calculating the water saturation in each injection and production direction and a mean value of the water saturation when infilling effect evaluation time is up in the infilled well position; determining a well position displacement in each injection and production direction; synthesizing an infilled well position displacement, determining a moving distance and a moving direction of the infilled well position, and determining a well position of the next infilled well; and calculating the water saturation in each injection and production direction in the infilled well position, and performing iterative calculation until an optimal infilled well position is obtained. According to the method, the infilled well position with maximized equilibrium displacement of the water-drive reservoir can be found and the optimization efficiency is improved.
Description
Technical field
The invention belongs to oil-gas field development field, in particular it relates to a kind of water-drive pool Encryption Well well location determines method,
When carrying out well pattern thickening in water-drive pool development process, under conditions of known encryption region and Encryption Well number, to realize oil reservoir
Equilibrium displacement maximum turns to target, determines Encryption Well well location, improves the recovery ratio of oil reservoir, improves the development effectiveness of oil reservoir.
Background technology
Well pattern thickening is the most direct effective means employed oil reservoir remaining oil and improve recovery ratio, and encryption is bored in position
Well, can improve injection-production relation, increase water drive sweep efficiency, improve the development effectiveness of oil reservoir.Need when determining oil reservoir Encryption Well well location
Consider the development features of oil reservoir and the anisotropism of reservoir, thus determine that optimum Encryption Well well location is very difficult.
By retrieval, the patent of invention of application number 201210076699.0 (Authorization Notice No. is CN102628353A) is disclosed
A kind of Well pattern edjustment and well pattern thickening method of the marine development of low-permeability oil reservoir of application, is set up using horizontal well and man-made fracture
Injection and extraction system, the water injection oil extraction method of Well pattern edjustment and well pattern thickening is realized by single-well injection-production;Application number
The patent of invention of 201510018306.4 (Authorization Notice No. is CN104615862A) discloses a kind of height based on evolution algorithm
The method that aqueous oil field determines well location, the method determines the well location dominant area of oil reservoir according to the remaining oil distribution of oil reservoir, utilizes
Kriging method and evolution algorithm determine optimum well location.
At present both at home and abroad method of the research worker using optimized algorithm in combination with reservoir numerical simulation determines Encryption Well well
Position, but the method number realization it is many, calculate the time it is long, and determine Encryption Well well location be numerical reservoir model grid position,
Actual encrypted position can not be accurately determined.Therefore, need further to study the determination method of Encryption Well well location, with compared with
The short calculating time obtains accurate Encryption Well well location.
The content of the invention
For the defect for overcoming current technology to exist, the present invention provides a kind of water-drive pool infilliing well well location determination side
Method.The method proposes the concept of Encryption Well Bit Shift, it is considered to oil reservoir development feature and reservoir heterogeneity, by iteration come fast
Speed, the Encryption Well well location for accurately finding optimum, at utmost to realize the equilibrium displacement of oil reservoir, improve the exploitation effect of oil reservoir
Really.
To achieve these goals, the present invention is employed the following technical solutions:
Water-drive pool Encryption Well well location determines method, it is characterised in that comprise the steps:
Step 1, collect with arrange the geologic information of oil reservoir, current Developing status and the region that need to be encrypted, plus
The well number of Mi Jing and well type, the flowing bottomhole pressure (FBHP) of injection-production well;
Step 2, at random given initial encryption well location, to each note of the encryption well location permeability in line direction, hole are adopted
Degree, initial water saturation carry out equivalent process;
Step 3, foundation note adopt the relation between direction water saturation and well spacing, and calculating is issued in the encryption well location and adds
Each note adopts the water saturation and its meansigma methodss in direction during the close effect assessment time;
Step 4, determine that each note adopts the well location displacement in direction;
Step 5, synthesis Encryption Well Bit Shift, it is determined that the displacement and moving direction of encryption well location, it is determined that next add
The well location of Mi Jing;
Each note adopts the water saturation in direction under the new encryption well location of step 6, calculating, then iterates to calculate, until obtaining most
Excellent Encryption Well well location.
Compared with prior art, the present invention has the advantages that:
1st, the displacement and moving direction that well location is encrypted in iterative process is determined according to Encryption Well Bit Shift, can be found
Realize the maximized Encryption Well well location of water-drive pool equilibrium displacement.
2nd, large-scale numerical simulation calculation need not be carried out during determining Encryption Well well location displacement, it is not necessary to set up
Numerical reservoir model, can quickly and accurately determine Encryption Well well location, improve optimization efficiency.
Description of the drawings
Fig. 1 determines flow chart for Encryption Well well location;
Fig. 2 is the permeability isopleth map of oil reservoir;
Fig. 3 is the oil-water relative permeability curve of oil reservoir;
Fig. 4 is the well location map of current oil reservoir;
Fig. 5 is the current water saturation isopleth map of oil reservoir;
Fig. 6 is the initial Encryption Well well location schematic diagram of oil reservoir;
Fig. 7 is the final Encryption Well well location schematic diagram of oil reservoir;
Fig. 8 is the oil saturation isopleth map of initial encryption scheme;
Fig. 9 is the water saturation isopleth map of final encipherment scheme.
Specific embodiment
As shown in figure 1, water-drive pool Encryption Well well location determines method, comprise the steps:
Step 1, collect with arrange the geologic information of oil reservoir, current Developing status and the region that need to be encrypted, plus
The well number of Mi Jing and well type, the flowing bottomhole pressure (FBHP) of injection-production well, comprise the following steps that:
Step 101, from region, carry out the overall investigation of oil reservoir, understand the geological condition of oil reservoir, collect and arrange oil
The geologic information of Tibetan, including:The permeability k distribution of oil reservoir, porosity φ distributions;Collect and arrange reservoir fluid physical characterization data,
Including oil phase viscosity μo, phase percolation curve;It is determined that the remaining oil distribution situation of current oil reservoir;
Step 102, collection and arrangement oil reservoir have at present the mouth coordinate and well type of well, delimit what oil reservoir need to be encrypted
Region, determines well number, the well type of Encryption Well, the flowing bottomhole pressure (FBHP) p of injection-production welle、piAnd cipher round results evaluation time t and encryption effect
Fruit starts evaluation time ti;
Step 2, at random given initial encryption well location, to each note of the encryption well location permeability in line direction, hole are adopted
Degree, initial water saturation carry out equivalent process, comprise the following steps that:
Step 201, at random given initial encryption well well location (xo,yo), determine in search radius R and constituted with the Encryption Well
The well of injection-production relation, it is determined that note adopts line number N and each note adopts length L of line1, L2..., LN;
Step 202, each note is adopted by the reservoir parameter of line is carried out equivalent process, equivalent method is as follows:
Each note is adopted by the permeability of line is carried out equivalent, computing formula is as follows:
Wherein, i adopts direction for i-th note;KiDirection equivalent permeability is adopted for i-th note;ki(x, y) is i-th note side of adopting
To the permeability at (x, y) place;
Each note is adopted by the porosity in direction is carried out equivalent, computing formula is as follows:
Wherein, φ adopts direction equivalent porosity for i-th note;φi(x, y) is the hole that i-th note adopts direction (x, y) place
Degree;
Each note is adopted by the initial water saturation in direction is carried out equivalent, computing formula is as follows:
Wherein, SwiThe initial equivalent water saturation in direction is adopted for i-th note;Swi(x, y) adopts direction (x, y) for i-th note
The initial water saturation at place;
Step 3, foundation note adopt the relation between direction water saturation and well spacing, and calculating is issued in the encryption well location and adds
Each note adopts the water saturation and its meansigma methodss in direction during the close effect assessment time, comprises the following steps that:
Step 301, oil-water relative permeability curve is fitted, fitting formula is as follows:
kro(Sw)=a (1-Sw)b (4)
Wherein, kro(Sw) for oil phase relative permeability;SwFor water saturation;A, b are fitting coefficient;
Step 302, determine it is each note the relation adopted between direction water saturation and well spacing, computing formula is as follows:
Step 303, relevant parameter is substituted into formula (5), calculate each note after producing the t times under encryption well location and adopt direction
Water saturation Sw1, Sw2..., SwN, and its meansigma methodss
Step 4, determine that each note adopts the well location displacement in direction, comprise the following steps that:
Step 401, each note is adopted into directional well-space it is divided into M deciles in the form of arithmetic progression, note adopts the heterogeneous body ginseng in direction
Number keeps constant, each note is calculated respectively and adopts water saturation numerical value of the direction under each well spacing, i-th note side of adopting of linear regression
Water saturation S upwardswiWith well spacing LiRelational expression Swi=gLi+ h, determines that each note adopts direction equilibrium displacement parameter with well spacing
Average rate of change Δ Swi/ΔLi=gi;
Step 402, determine that each note adopts the direction vector of line, computing formula is as follows:
Wherein, diFor the direction vector that i-th note adopts direction;(xo,yo) it is Encryption Well well location;(xi,yi) be and Encryption Well
Constitute the well location that i-th note adopts direction;
Step 403, according to the difference of each water saturation and average water saturation noted and adopt direction, and step 401
It is determined that each note adopt rate of change of the direction water saturation with well spacing, determine it is each note adopt the well location distance that direction need to move, calculate
Formula is as follows:
Wherein, DiFor the distance that Encryption Well need to adopt direction movement along i-th note;
Step 404, the displacement that the note that step 402 determines is adopted the well location that line direction vector and step 403 determine,
The Encryption Well well location displacement that note adopts direction is constituted, computing formula is as follows:
Di=Di·di (8)
Wherein, DiFor the displacement that Encryption Well adopts direction along i-th note;
Step 5, synthesis Encryption Well Bit Shift, it is determined that the displacement and moving direction of encryption well location, it is determined that next add
The well location of Mi Jing, comprises the following steps that:
Step 501, the moving displacement for adopting direction along each note according to Encryption Well, according to the mode of Vector modulation, synthesis encryption
The well location displacement of well, computing formula is as follows:
S=D1+D2+…+Dn (9)
Wherein, S is the well location displacement of Encryption Well;
Step 502, according to the well location displacement of current Encryption Well, it is determined that next Encryption Well well location, computing formula is as follows:
(x', y')=(x0,y0)+S (10)
Wherein, (x ', y ') it is next Encryption Well well location;
Step 503, Encryption Well well location is judged whether in the range of encryption, if in the range of encryption, carrying out step 6,
Otherwise, with well location (xo,yo) with the line of (x ', y ') and the intersection point of encryption range boundary as next Encryption Well well location.
Each note adopts the water saturation in direction under the new encryption well location of step 6, calculating, then iterates to calculate, until obtaining most
Excellent Encryption Well well location.
The water-drive pool Encryption Well well location for describing the present invention in detail so that the Encryption Well well location of certain oil reservoir determines as an example below is true
Determine method, but be not limited to the practical range of the present invention.
The embodiment oil reservoir size is 360m × 360m × 360m, and whole oil reservoir divides a set of straticulate exploration, reservoir pore
Spend for 0.3, permeability is inhomogeneous distribution, permeability isopleth map as shown in Figure 2 (unit, 10-3μm2), oil viscosity is
5mPs, residual oil saturation is 0.2, and oil-water relative permeability curve is as shown in Figure 3.The embodiment oil reservoir is at present using rule
Inverted nine-spot pattern well spacing, current well location map is as shown in 4 figures, and the remaining oil distribution isopleth map of oil reservoir is as shown in Figure 5.
Water injection well I turns producing well in the middle of embodiment, encrypts 4 mouthfuls of water injection wells I1, I2, I3, I4 so as to constitute 45 wells
Net, with the development effectiveness that oil reservoir after realizing 3 years reaches equilibrium displacement.All injection-production wells are produced using stable bottom hole pressure, producing well
Flowing bottomhole pressure (FBHP) is 18MPa, and water injection well flowing bottomhole pressure (FBHP) is 22MPa.
Given 4 initial encryption well locations for (90,90) (90,270) (270,90) (270,270), as shown in fig. 6, encryption
Well scope be respectively (1~180,1~180) (1~180,180~360) (180~360,1~180) (180~360,180~
360), 4 mouthfuls of producing wells of each initial encryption well well location and surrounding constitute injection-production relation, then each encryption well location has 4 note companies of adopting
Line, per bar, note adopts wire length and is 126m.
Adopting the permeability on line, initial water saturation to each note according to formula (1), formula (3) carries out equivalent process,
Be calculated it is each note adopt line equivalent parameterss numerical value it is as shown in table 1.
Respectively note adopts line permeability and initial water saturation equivalent numerical value to table 1
The oil relative permeability in Fig. 3 oil-water relative permeability curves is fitted according to formula (4), is fitted
Expression formula is kro(Sw)=1.99 × (1-Sw)3.01;
The relevant parameter that I1-P1 notes adopt direction is substituted into into formula (5), I1-P1 notes adopt the aqueous of direction after being calculated 3 years
Saturation is:
In the same manner, the I1-P2 notes for being calculated I1 well groups using formula (5) adopt direction, I1-P4 notes and adopt direction, I 1-I and note and adopt
The water saturation in direction is respectively 0.8210,0.8822,0.6162.
Calculate I1 well groups 4 note the average water saturation for adopting direction:
In the same manner, each note for being calculated I2, I3, I4 injection-production well group using formula (5) adopts direction water saturation, and calculating contains
The meansigma methodss of water saturation, result of calculation is as shown in table 2.
Respectively note adopts direction water saturation to table 2
Each note is adopted into directional well-space and is divided into 20 deciles in the form of arithmetic progression, note is adopted the heterogeneity parameter in direction and kept not
Become, each note calculated respectively and adopts water saturation of the direction under each well spacing, linear regression respectively note adopt on direction containing water saturation
Relational expression S of degree and well spacingwi=gLi+ h, determines that each note adopts average rate of change Δ S of the direction water saturation with well spacingwi/ΔLi
=gi, the rate of change for returning the relational expression that obtains and water saturation is as shown in table 3.
Respectively note adopts direction water saturation linear regression formula and rate of change to table 3
Determine that I1-P1 notes adopt the direction vector of line according to formula (6):
Determine that I1-P1 notes adopt the distance that direction need to move according to formula (7):
DI1-P1=(0.7305-0.6210)/| -2.0804 × 10-5|=5263cm=52.63m
I1-P1 notes are adopted by line direction vector according to formula (8) and adopts line direction with the displacement composition note of well location
Encryption Well well location displacement:
DI1-P1=52.63 × (0.7071,0.7071)=(37.21,37.21)
In the same manner, the Encryption Well well location displacement that other 15 notes adopt direction is calculated according to formula (6), formula (7), formula (8),
As a result it is as shown in table 4.
Respectively note adopts direction Encryption Well well location displacement to table 4
According to the well location displacement of formula (9) Vector modulation I1 wells:
SI1=(37.21,37.21)+(48.04, -44.83)+(- 110.05,117.91)+(- 46.25, -46.25)=(-
71.05,64.04)
In the same manner, using formula (9) Vector modulation I2, I3, I4 well location displacement be respectively (12.89,106.49),
(141.88,-41.71)、(26.31,0.94)。
The next iterative cryptographic well well location of Encryption Well I1 is calculated according to formula (10):
(xI1',yI1')=(90,90)+(- 71.05,64.04)=(18.95,154.04)
Due to the well location encryption scope (1~180,1~180) in, so need not be processed.
In the same manner, the next iterative cryptographic well well location for being calculated Encryption Well I2, I3, I4 using formula (10) is respectively
(102.88,376.49)、(411.88,48.28)、(296.31,269.91)。
Due to encrypting well location I2 (102.88,376.49) and encryption well location I3 (411.88,48.28) beyond encryption model
Enclose, BORDER PROCESSING need to be carried out to the two encryption well locations.
Calculate the line of the initial well locations (90,270) of I2 and encryption well location (102.88,376.49) and cryptographic boundary (1~
180,180~intersection point 360), determine I2 wells next one iterative cryptographic well well location for (100.89,360).
In the same manner, calculate the initial well locations of I3 (270,90) with encryption well location (411.88, line 48.28) and cryptographic boundary (1
~180,180~intersection point 360), determine I3 wells next one iterative cryptographic well well location for (360,63.54).
Step 6, with I1, I2, I3, I4 it is new encryption well location repeat step 3-6 iterative calculation, until obtain optimum encryption
Well well location.Finally give realize oil reservoir equilibrium displacement maximized Encryption Well I1, I2, I3, I4 optimal well location be I1 for (78,
90), I2 for (84,114), I3 for (210,90), I4 for (222,222), as shown in Figure 7.
The remaining oil saturation point under remaining oil saturation distribution and final encryption well location under contrast initial rules encryption
Cloth, as shown in Figure 8, Figure 9, it can be seen that when the Encryption Well well location obtained with iteration optimization is produced, the displacement of oil reservoir can be made
It is more balanced, the development effectiveness of oil reservoir can be improved.Therefore data above can be used as the embodiment oil reservoir development design
Theoretical foundation.
The embodiment of the present invention selects water saturation as the parameter for characterizing oil reservoir equilibrium displacement, without departing substantially from essence of the invention
In the case of god and essence, those of ordinary skill in the art also can select oil saturation, recovery percent of reserves, accumulation injection hole
Volume, moisture content can make various corresponding deformations and change as equilibrium displacement parameter is characterized, but all should cover in the present invention
Protection domain within.Therefore, within the protection domain that protection scope of the present invention should be defined with claim.
Claims (5)
1. a kind of water-drive pool Encryption Well well location determines method, it is characterised in that comprise the steps:
Step 1, the geologic information for collecting and arranging oil reservoir, current Developing status and the region that need to be encrypted, Encryption Well
Well number and well type, the flowing bottomhole pressure (FBHP) of injection-production well;
Step 2, at random given initial encryption well location, the permeability in line direction, porosity, just are adopted to each note of the encryption well location
Beginning water saturation carries out equivalent process;
Step 3, foundation note adopt the relation between direction water saturation and well spacing, calculate and are issued to encryption effect in the encryption well location
Each note adopts the water saturation and its meansigma methodss in direction during fruit evaluation time;
Step 4, determine that each note adopts the well location displacement in direction, comprise the following steps that:
Step 401, each note is adopted into directional well-space it is divided into M deciles in the form of arithmetic progression, note is adopted the heterogeneity parameter in direction and protected
Hold constant, each note is calculated respectively and adopts water saturation numerical value of the direction under each well spacing, i-th note of linear regression is adopted on direction
Water saturation SwiWith well spacing LiRelational expression Swi=gLi+ h, determines that each note adopts direction equilibrium displacement parameter with the flat of well spacing
Rate of change Δ Swi/ΔLi=gi;Wherein, g is the slope of regression straight line, and h is the intercept of regression straight line;
Step 402, determine that each note adopts the direction vector of line, computing formula is as follows:
Wherein, diFor the direction vector that i-th note adopts direction;(xo,yo) it is Encryption Well well location;(xi,yi) it is to constitute with Encryption Well
I-th note adopts the well location in direction;
Step 403, according to the difference of each water saturation and average water saturation noted and adopt direction, and step 401 determines
Each note adopt rate of change of the direction water saturation with well spacing, determine it is each note adopt the well location distance that direction need to move, computing formula
It is as follows:
Wherein, DiFor the distance that Encryption Well need to adopt direction movement along i-th note;
Step 404, the displacement that the note that step 402 determines is adopted the well location that line direction vector and step 403 determine, are constituted
Note adopts the Encryption Well well location displacement in direction, and computing formula is as follows:
Di=Di·di (3)
Wherein, DiFor the displacement that Encryption Well adopts direction along i-th note;
Step 5, synthesis Encryption Well Bit Shift, it is determined that the displacement and moving direction of encryption well location, it is determined that next Encryption Well
Well location;
Each note adopts the water saturation in direction under the new encryption well location of step 6, calculating, then iterates to calculate, until obtaining optimum
Encryption Well well location.
2. water-drive pool Encryption Well well location according to claim 1 determines method, it is characterised in that:The concrete step of step 1
It is rapid as follows:
Step 101, from region, carry out the overall investigation of oil reservoir, understand the geological condition of oil reservoir, collect and arrange oil reservoir
Geologic information, including:The permeability k distribution of oil reservoir, porosity φ distributions;Collect and arrange reservoir fluid physical characterization data, including
Oil phase viscosity μo, phase percolation curve;It is determined that the remaining oil distribution situation of current oil reservoir;
Step 102, collection and arrangement oil reservoir have at present the mouth coordinate and well type of well, delimit the area that oil reservoir need to be encrypted
Domain, determines well number, the well type of Encryption Well, the flowing bottomhole pressure (FBHP) p of injection-production welle、piAnd cipher round results evaluation time t and cipher round results
Start evaluation time ti。
3. water-drive pool Encryption Well well location according to claim 2 determines method, it is characterised in that:The concrete step of step 2
It is rapid as follows:
Step 201, at random given initial encryption well well location (xo,yo), determine that constituting note with the Encryption Well adopts in search radius R
The well of relation, it is determined that note adopts line number N and each note adopts length L of line1, L2..., LN;
Step 202, each note is adopted by the reservoir parameter of line is carried out equivalent process, equivalent method is as follows:
Each note is adopted by the permeability of line is carried out equivalent, computing formula is as follows:
Wherein, i adopts direction for i-th note;KiDirection equivalent permeability is adopted for i-th note;ki(x, y) adopts direction for i-th note
The permeability at (x, y) place;
Each note is adopted by the porosity in direction is carried out equivalent, computing formula is as follows:
Wherein, φ adopts direction equivalent porosity for i-th note;φi(x, y) is the porosity that i-th note adopts direction (x, y) place;
Each note is adopted by the initial water saturation in direction is carried out equivalent, computing formula is as follows:
Wherein, SwiThe initial equivalent water saturation in direction is adopted for i-th note;Swi(x, y) adopts direction (x, y) place for i-th note
Initial water saturation.
4. water-drive pool Encryption Well well location according to claim 3 determines method, it is characterised in that:The concrete step of step 3
It is rapid as follows:
Step 301, oil-water relative permeability curve is fitted, fitting formula is as follows:
kro(Sw)=a (1-Sw)b (7)
Wherein, kro(Sw) for oil phase relative permeability;SwFor water saturation;A, b are fitting coefficient;
Step 302, determine it is each note the relation adopted between direction water saturation and well spacing, computing formula is as follows:
Wherein, L is well spacing;
Step 303, relevant parameter is substituted into formula (5), calculate each note after producing the t times under encryption well location and adopt the aqueous of direction
Saturation Sw1, Sw2..., SwN, and its meansigma methodss
5. water-drive pool Encryption Well well location according to claim 4 determines method, it is characterised in that:The concrete step of step 5
It is rapid as follows:
Step 501, the moving displacement for adopting direction along each note according to Encryption Well, according to the mode of Vector modulation, synthesize Encryption Well
Well location displacement, computing formula is as follows:
S=D1+D2+…+Dn (9)
Wherein, S is the well location displacement of Encryption Well;
Step 502, according to the well location displacement of current Encryption Well, it is determined that next Encryption Well well location, computing formula is as follows:
(x', y')=(x0,y0)+S (10)
Wherein, (x ', y ') it is next Encryption Well well location;
Step 503, Encryption Well well location is judged whether in the range of encryption, if in the range of encryption, carrying out step 6, otherwise,
With well location (xo,yo) with the line of (x ', y ') and the intersection point of encryption range boundary as next Encryption Well well location.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610446421.6A CN106126899B (en) | 2016-06-20 | 2016-06-20 | Method for determining well positions of infilled wells of water-drive reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610446421.6A CN106126899B (en) | 2016-06-20 | 2016-06-20 | Method for determining well positions of infilled wells of water-drive reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106126899A CN106126899A (en) | 2016-11-16 |
CN106126899B true CN106126899B (en) | 2017-05-17 |
Family
ID=57470277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610446421.6A Expired - Fee Related CN106126899B (en) | 2016-06-20 | 2016-06-20 | Method for determining well positions of infilled wells of water-drive reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106126899B (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107829718B (en) * | 2017-02-10 | 2019-11-08 | 中国石油化工股份有限公司 | Oil reservoir well pattern and injection-production program optimum design method based on balanced water drive theory |
CN109753671B (en) * | 2017-11-07 | 2023-06-02 | 中国石油大学(华东) | Oil reservoir accurate well position optimization method based on fish swarm algorithm |
CN108561115B (en) * | 2018-01-10 | 2019-02-22 | 中国石油大学(华东) | A kind of irregular flood pattern unit design method of heterogeneous reservoir |
CN109403935B (en) * | 2018-09-30 | 2020-10-09 | 中海石油(中国)有限公司 | Oil yield calculation method for oil sand SAGD encrypted well |
CN111119840B (en) * | 2018-11-01 | 2023-05-26 | 中国石油天然气股份有限公司 | Well position optimization method for high water-cut period oil reservoir |
CN109632557B (en) * | 2019-01-22 | 2021-11-16 | 中国矿业大学 | Gas-liquid two-phase saturated coal rock sample experimental device and saturation testing method |
CN110485988A (en) * | 2019-08-19 | 2019-11-22 | 中国石油天然气股份有限公司大港油田分公司 | A kind of low permeability reservoir simulation grid refinement method and device |
CN112561111B (en) * | 2019-09-10 | 2024-03-05 | 中国石油天然气股份有限公司 | Oilfield water flooding wave and coefficient prediction method and device |
CN111396021B (en) * | 2020-03-04 | 2022-03-29 | 中国石油天然气股份有限公司 | Well arrangement optimization method for newly drilled and encrypted well in old oil field developed by water injection |
CN113898331A (en) * | 2020-06-22 | 2022-01-07 | 中国石油化工股份有限公司 | Well spacing optimization method for oil reservoir balanced displacement differentiation in high water cut period |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150051838A1 (en) * | 2013-08-15 | 2015-02-19 | Chevron U.S.A. Inc. | System and method of determining and optimizing waterflood performance |
CN104615862B (en) * | 2015-01-14 | 2017-09-08 | 中国石油天然气股份有限公司 | Method for determining well position of high-water-content oil field based on evolutionary algorithm |
-
2016
- 2016-06-20 CN CN201610446421.6A patent/CN106126899B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106126899A (en) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106126899B (en) | Method for determining well positions of infilled wells of water-drive reservoir | |
CN103939066B (en) | An a kind of note adopts the method that well group determines water injection rate determination oilwell produced fluid amount more | |
CN107829718A (en) | Oil reservoir well pattern and injection-production program Optimization Design based on balanced water drive theory | |
CN106437674B (en) | Imitative water injection of horizontal well well pattern adaptation method | |
CN108868748A (en) | A kind of calculation method of shale gas horizontal well refracturing crack cracking pressure | |
Ogbeiwi et al. | An approach to waterflood optimization: case study of the reservoir X | |
CN104453834A (en) | Injection-production relation optimizing and adjusting method for well group | |
CN104632202B (en) | Method and device for determining dry clay three-porosity logging parameter values | |
CN108343420A (en) | Multi-factor collaborative analysis industrial large well group well arrangement method | |
CN103867175A (en) | Steam flooding well pattern structure and steam flooding development method thereof | |
CN105652342A (en) | Sectional type ground temperature gradient fitting method based on stratigraphic unconformity surface | |
CN106547973A (en) | A kind of quantitative description of injection for heavy oil reservoir steam course channeling passage | |
CN107448197A (en) | A kind of method of marine facies sandstone oil reservoir fluid-channeling channel space development quantitative judge | |
CN106372810A (en) | TBM (tunnel boring machine) quantitative type selection method by use of efficacy coefficient method | |
CN104564069B (en) | The dynamic settlement prediction in a kind of ground based on square mesh method and reclamation method | |
CN107818518A (en) | A kind of method in quantization signifying oil shale in-situ exploitation constituency | |
CN103390108A (en) | Well type design method and device based on compressibility of fractured geologic body | |
Sennhauser et al. | A practical numerical model to optimize the productivity of multistage fractured horizontal wells in the cardium tight oil resource | |
CN105089659B (en) | Conglomerate oil reservoir seepage unit identification method | |
CN107219564A (en) | Method and device for processing vertical well connecting profile | |
CN110500089A (en) | Shale gas horizontal well stratification fracture evaluation method based on array induction logging data | |
Feizabadi et al. | An integrated approach to building history-matched geomodels to understand complex long lake oil sands reservoirs, part 2: Simulation | |
CN107766639A (en) | The computational methods of the natural gas lateral migration ultimate range of coefficient are reduced based on pressure | |
Leemput et al. | Full-field reservoir modeling of Central Oman gas-condensate fields | |
CN104459820B (en) | The method for building up of multiple dimensioned pressure field |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170517 |