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 PDF

Info

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
Application number
CN201610446421.6A
Other languages
Chinese (zh)
Other versions
CN106126899A (en
Inventor
冯其红
陈红伟
张先敏
王相
邢祥东
胡鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201610446421.6A priority Critical patent/CN106126899B/en
Publication of CN106126899A publication Critical patent/CN106126899A/en
Application granted granted Critical
Publication of CN106126899B publication Critical patent/CN106126899B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject 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

Water-drive pool Encryption Well well location determines method
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:
d i = ( x i - x 0 , y i - y 0 ) ( x i - x 0 ) 2 + ( y i - y 0 ) 2 - - - ( 1 )
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:
D i = ( S w ‾ - S w i ) / | g i | - - - ( 2 )
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:
K i = 1 L i ∫ l k i ( x , y ) dx i - - - ( 4 )
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:
φ i = 1 L i ∫ l φ i ( x , y ) dx i - - - ( 5 )
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:
S w i = 1 L i ∫ l S w i ( x , y ) dx i - - - ( 6 )
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:
S w = 1 - [ ( 1 - S w i ) 1 - b + a ( b - 1 ) k ( p e - p i ) ( t - t i ) φμ o L 2 ] 1 1 - b - - - ( 8 )
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.
CN201610446421.6A 2016-06-20 2016-06-20 Method for determining well positions of infilled wells of water-drive reservoir Expired - Fee Related CN106126899B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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