CN108868712A

CN108868712A - A kind of oil reservoir development production optimization method and system based on connectivity method

Info

Publication number: CN108868712A
Application number: CN201711286901.1A
Authority: CN
Inventors: 赵辉; 曹琳; 李国浩; 史永波; 刘伟
Original assignee: Yangtze University
Current assignee: Yangtze University
Priority date: 2017-12-07
Filing date: 2017-12-07
Publication date: 2018-11-23
Anticipated expiration: 2037-12-07
Also published as: CN108868712B

Abstract

The present invention relates to oil reservoir production method, in particular to a kind of oil reservoir development production optimization method and system based on connectivity method.The present invention establishes accurate leading edge method for tracing as object using interwell communication unit and carries out saturation computation, obtains each layer water-oil phase output dynamic at well point；By being connected to model parameter to the dynamic auto fitting inverting of oil reservoir history, obtain conductivity between injection-production well, flow is split point and the information such as water-injection efficiency, and layering dynamic production and injection proration Automatic Optimal Design is carried out to oil reservoir by iterative calculation on this basis, inefficient water drive direction flow is reduced, improves note and adopts contradiction.It is shown using example, the calculating of institute's construction method is quick, can accurately carry out saturation degree tracking and calculates output dynamic, gained seperated layer water injection scheme preferably realizes decreasing water cut and increasing oil effect, increases oil nearly 50,000 side prediction year after practical oil reservoir optimization, provides new approaches for oil field injection and extraction scheme real-time optimization.

Description

A kind of oil reservoir development production optimization method and system based on connectivity method

Technical field

The present invention relates to oil reservoir production method, in particular to a kind of oil reservoir development Optimizing manufacture side based on connectivity method Method and system.

Background technique

The all preservations of crude oil in China reserves 90% are in clastic reservoir rock, since longitudinal span is big, heterogeneity is strong, interlayer It is particularly thorny.Seperated layer water injection is to mitigate interlayer contradiction, increases production the important measures of stable yields, and Oil Field tentatively realizes layering at present Real-time commissioning process technology.But oil reservoir seperated layer water injection Design Method is urgently studied, and realizes that waterflood project and process matched therewith are neat Head is gone forward side by side, and actual field application could be more preferably met.

The design of artificial recharge scheme optimization depends on reservoir engineering method or numerical simulation technology, and design scheme is random Strong, time-consuming and laborious and easy " missing " optimal case of property.The real-time Optimizing manufacture theory of Reservoir behavior is that current oil field is infused automatically and adopted The important research hot spot of conceptual design adopts oil-water well note mainly in conjunction with Research Numerical Simulation Techique and Optimum Theory The design of parameter is converted into Optimal Control Model solution, and use is such as with gradient, Stochastic Perturbation Gradient algorithm, heuritic approach Deng maximizing pattern function, and then the optimal working system of automatic calculation, but since to solve difficulty higher for gradient, digital-to-analogue operand compared with Greatly, the dimension of actual optimization problem is higher, brings the inefficiency problem of optimization algorithm, and this method is premature apart from practical application.

Inter well connectivity is the important evidence for filling the water optimization design, has adopted the inter well connectivity model of dynamic data based on note Gradually from it is single-phase develop to water-oil phase, single layer develop to multilayer predict new stage, have calculate quickly, can quantitatively characterizing well Between connected relation the advantages that, oil reservoir scheme evaluation design in be gradually applied.But current connectivity model carries out complicated oil When hiding grease dynamic prediction, saturation degree method for tracing is not yet applied in oil reservoir demixing injection-production prioritization scheme, therefore cannot be practical Ground optimizes current oil reservoir recovery scheme, to improve oil reservoir production efficiency.

Summary of the invention

The present invention provides a kind of oil reservoir development production optimization methods and system based on connectivity method, solve above The technical problem.

The technical solution that the present invention solves above-mentioned technical problem is as follows：A kind of oil reservoir development production based on connectivity method Optimization method includes the following steps：

Step 1, it establishes for simulating between well grease dynamically initial connectivity modeling, the initial connectivity modeling will be oily Hiding injection and extraction system is separated into the connected unit between well and well, and each connected unit is joined by conductivity and control two features of volume Number characterization；

Step 2, based on connected unit and by matter balance equation and waterflood front promote equation carry out well point pressure and Flow rate calculation between well, and saturation degree tracking is carried out, it generates the current oil aquatic products at well point and goes out dynamic indicator；

Step 3, go out dynamic indicator to the current oil aquatic products using the history grease output dynamic indicator of oil reservoir to carry out certainly Dynamic fitting inverting, adjusts the characteristic parameter of the initial connectivity modeling, and generate the connectivity modeling after optimization；

Step 4, last moment all water injection wells are calculated in the water filling of each interval by the connectivity modeling after the optimization Efficiency, and the block for obtaining the oil reservoir development history last moment is averaged water-injection efficiency, according to the water-injection efficiency of each interval and described The be averaged comparison result of water-injection efficiency of block carries out the first suboptimization to the current development and production method of water injection well；After optimization Water injection well development and production method calculate water filling of each well to surrounding oil well and split a point coefficient, and a point coefficient is split according to the water filling Second suboptimization is carried out to the water injection well development and production method after optimization.

The beneficial effects of the invention are as follows：The comprehensive interwell communication relationship of technical solution of the present invention and water-injection efficiency inversion result A kind of new fine water filling optimization method of stratified reservoir layering is established, is to reduce inefficient water breakthrough flow and improve water-injection efficiency Target solves automatic formulation oil-water well demixing injection-production scheme by iteration optimization, and the injection allocation that realization stratified reservoir dynamic matches production is set Meter.Design method calculates quickly simultaneously, can accurately carry out saturation degree tracking and calculate output dynamic, gained seperated layer water injection side Case preferably realizes decreasing water cut and increasing oil effect, increases oil nearly 50,000 side prediction year after practical oil reservoir optimization, is that oil field injection and extraction scheme is real-time Optimization provides new approaches.

Based on the above technical solution, the present invention can also be improved as follows.

Further, the step 2 specifically includes following steps：

S201 establishes matter balance equation using the connected unit as object：

Wherein, t is the time, and i, j are well serial number, and k is oil reservoir serial number, N_lFor oil reservoir number, N_wFor injection-production well number, p_i, p_jIt is Average pressure in i well and jJing draining area, q_iFor the i-th well yield, injection is positive, output is negative, C_tkFor the synthesis of kth layer The compressed coefficient, V_pikFor the drainage volume of the i-th well of kth layer；

S202 carries out implicit difference solution to the matter balance equation using liquid and level pressure two ways is determined, generates well Between flow between connected unit

Wherein, n is current time, q_ijkFor kth layer i well and j interwell communication unit inner stream flow, p_i, p_jFor the i-th well and Average pressure in j well draining area, T_ijkFor the average conduction rate between kth layer i well and j well；

S203, when upstream well point j saturation degree is between leading edge saturation degree and irreducible water saturation and downstream well point i is beam When tiing up water saturation state, using the method for solving of shock wave, fltting speed of the saturation degree face in connected unit is denoted as ν_ijk：

If interwell communication element length is L_ijk, saturation degree face is in connected unit from the time of upstream node j movement For Δ n, well point i saturation computation result is at this time：

Wherein, n is current time, and n ' is last moment, f_wFor moisture content, f_w′(s_wijk) it is to be tracked in kth layer from jth well To the derivate of water content ratio of the i-th Jing Chu, f_w′(S_wjk) it is derivate of water content ratio of the jth well in kth layer, F_vijkFor from jth well stream to i-th Zero dimension integrated flow of the well in kth layer, L_ijkIt is i well and j well in kth layer distance, S_wijkFor kth layer i well and j interwell communication The water saturation of unit, S_wikWater saturation for the i-th well in kth layer, V_ijkFor kth layer i well and j interwell communication volume, ν_ijkFor saturation degree face fltting speed inside kth layer i well and j well connected unit.

Beneficial effect using above-mentioned further scheme is：Saturation degree method for tracing is carried out in this further technical solution Improve, solve in displacement process when upstream well point saturation degree is between irreducible water saturation and displacing front saturation degree When, the saturation degree method for tracing of the prior art is difficult to the problem of accurately solving saturation distribution state, and saturation computation method is more Add accurately, each layer water-oil phase output dynamic indicator is also more accurate at the well point of generation, further improves oil reservoir development Effect.

Further, the step 4 is specially：

S401, the connectivity modeling after being optimized by step 3 calculate last moment all water injection wells in the water filling of each interval Efficiency, and the block for acquiring the oil reservoir development history last moment is averaged water-injection efficiency；

S402 evaluates each interval water-injection efficiency of water injection well, if the water-injection efficiency of water injection well is flat less than the block Equal water-injection efficiency then carries out drop note, otherwise carries out augmented injection, and generate the oil reservoir development production method after the first suboptimization, specifically Each layer injection rate of water injection well after optimizationFor：

Wherein, η_ikWeight is adjusted in kth layer for the well, specific formula for calculation is as follows：

Wherein Wef is the water-injection efficiency of block, Wef_ikIt is i-th mouthful of water injection well in kth layer water-injection efficiency, Wef_max, Wef_min For minimum and peak efficiency in current water injection well, α is to produce liquid adjustment index, η_ikLiquid measure for i-th mouthful of water injection well in kth layer adjusts Coefficient, η_max, η_minFor the bound of liquid measure regulation coefficient, m is last moment optimizing phase；

S403 calculates last moment each well to surrounding oil well using the oil reservoir development production method after the first suboptimization A point coefficient is split in water filling, and is split point coefficient according to the water filling and carried out the second suboptimization to oil reservoir development production method, is specifically matched and is produced Optimizing calculation formula is：

Wherein A_ijkA point coefficient is split in the water filling of kth layer to jth mouth well for i-th mouthful of water injection well, m is last moment optimizing phase, η_ikFor i-th mouthful of water injection well kth layer liquid measure regulation coefficient,For each layer injection rate of water injection well after optimization, N_IFor water injection well number, q_iFor the i-th well yield, injection is positive, output is negative；

S404 after simulating preset time using the oil reservoir development production method after the second suboptimization, repeats step S401- S403 again optimizes current oil reservoir development and production method.

Beneficial effect using above-mentioned further scheme is：Further technical solution can reduce inefficient water drive direction stream for this Amount improves block water-injection efficiency, while realizing that oil reservoir injection adopts policy dynamic optimization.

Further, η_min=-0.5, η_max=0.5, α value are 2.

In order to solve technical problem of the invention, a kind of oil reservoir development Optimizing manufacture based on connectivity method is additionally provided System, including：

Model building module, it is described initially to connect for establishing for simulating between well grease dynamically initial connectivity modeling Oil reservoir injection and extraction system is separated into the connected unit between well and well by general character model, and each connected unit is by conductivity and control volume Two characteristic parameter characterizations of product；

Computing module, for promoting equation to carry out well point based on connected unit and by matter balance equation and waterflood front Flow rate calculation between pressure and well, and saturation degree tracking is carried out, it generates the current oil aquatic products at well point and goes out dynamic indicator；

Inverting module, for going out dynamic indicator to the current oil aquatic products using the history grease output dynamic indicator of oil reservoir Automatic Fitting inverting is carried out, the characteristic parameter of the initial connectivity modeling is adjusted, and generates the connectivity modeling after optimization；

Optimization module, for calculating last moment all water injection wells in each interval by the connectivity modeling after the optimization Water-injection efficiency, and the block for obtaining the oil reservoir development history last moment is averaged water-injection efficiency, according to the water-injection efficiency of each interval The first suboptimization is carried out to the current development and production method of water injection well with the be averaged comparison result of water-injection efficiency of the block；Using Water injection well development and production method after optimization calculates water filling of each well to surrounding oil well and splits a point coefficient, and is split according to the water filling Coefficient is divided to carry out the second suboptimization to the water injection well development and production method after optimization.

Further, the computing module specifically includes：

Model foundation unit, for establishing following matter balance equation by object of the connected unit：

First computing unit determines liquid and level pressure two ways to matter balance equation progress implicit difference for using It solves, generates the flow between interwell communication unitFor：

Second computing unit, for when upstream well point j saturation degree between leading edge saturation degree and irreducible water saturation and When downstream well point i is irreducible water saturation state, using the method for solving of shock wave, by propulsion of the saturation degree face in connected unit Speed is denoted as ν_ijk：

And interwell communication element length is set as L_ijk, saturation degree face in connected unit from upstream node j movement when Between be Δ n, generate well point i saturation computation result be：

Further, the optimization module specifically includes：

Third computing unit, for calculating last moment all water injection wells in each interval by the connectivity modeling after optimization Water-injection efficiency, and the block for acquiring the oil reservoir development history last moment is averaged water-injection efficiency；

First optimization unit, for evaluating each interval water-injection efficiency of water injection well, if the water-injection efficiency of water injection well is small It is averaged water-injection efficiency in the block, then carries out drop note, otherwise carry out augmented injection, and it is raw to generate the oil reservoir development after the first suboptimization Production method, each layer injection rate of water injection well after specific optimizationFor：

Second optimization unit, for calculating last moment each well using the oil reservoir development production method after the first suboptimization A point coefficient is split in water filling to surrounding oil well, and is split point coefficient according to the water filling and carried out the second suboptimum to oil reservoir development production method Change, it is specific to be with production optimization calculation formula：

Wherein A_ijkA point coefficient is split in the water filling of kth layer to jth mouth well for i-th mouthful of water injection well, m is last moment optimizing phase, η_ikFor i-th mouthful of water injection well kth layer liquid measure regulation coefficient,For each layer injection rate of water injection well after optimization, N_IFor water injection well number, q_iFor the i-th well yield, injection is positive, output is negative.

Further, η_min=-0.5, η_max=0.5, α value are 2.

Detailed description of the invention

Fig. 1 is that the process for the oil reservoir development production optimization method based on connectivity method that the embodiment of the present invention 1 provides is shown It is intended to；

Fig. 2 is that saturation degree tracks schematic diagram；

Fig. 3 is concept reservoir model well location schematic diagram；

Fig. 4 is that W2 well produces water curve comparison figure；

Fig. 5 a is the first layer permeability of concept Reservoir Permeability field in embodiment 2；

Fig. 5 b is the second layer permeability of concept Reservoir Permeability field in embodiment 2；

Fig. 6 a is block oil-producing fitting result in embodiment 2；

Fig. 6 b is individual well production rate fitting result in embodiment 2；

Fig. 7 a is first layer reservoir model parametric inversion result in embodiment 2；

Fig. 7 b is second layer reservoir model parametric inversion result in embodiment 2；

Fig. 8 a is that a point situation is split in first layer water filling in embodiment 2；

Fig. 8 b is that a point situation is split in second layer water filling in embodiment 2；

Fig. 9 is seperated layer water injection efficiency rating in embodiment 2；

Figure 10 a is water injection well prioritization scheme in embodiment 2；

Figure 10 b is producing well prioritization scheme in embodiment 2；

Figure 11 is optimization front and back block production target variation in embodiment 2；

Figure 12 a is each layer remaining oil comparison before first layer optimization in embodiment 2；

Figure 12 b is each layer remaining oil comparison after first layer optimization in embodiment 2；

Figure 12 c is each layer remaining oil comparison before second layer optimization in embodiment 2；

Figure 12 d is each layer remaining oil comparison after second layer optimization in embodiment 2；

Figure 13 a is that block tires out oil-producing fitting result in embodiment 3；

Figure 13 b is block moisture content fitting result in embodiment 3；

Figure 14 a is that component curve figure is split in first layer water filling in embodiment 3；

Figure 14 b is that component curve figure is split in second layer water filling in embodiment 3；

Figure 14 c is that component curve figure is split in third layer water filling in embodiment 3；

Figure 15 is seperated layer water injection efficiency chart in embodiment 3；

Figure 16 a is water injection well prioritization scheme in embodiment 3；

Figure 16 b is producing well prioritization scheme in embodiment 3；

Figure 17 is block optimization front and back index comparison diagram in embodiment 3；

Figure 18 is that the structure for the oil reservoir development production optimization system based on connectivity method that the embodiment of the present invention 4 provides is shown It is intended to.

Specific embodiment

The principle and features of the present invention will be described below with reference to the accompanying drawings, and the given examples are served only to explain the present invention, and It is non-to be used to limit the scope of the invention.

Fig. 1 is that the process for the oil reservoir development production optimization method based on connectivity method that the embodiment of the present invention 1 provides is shown It is intended to, as shown in Figure 1, including the following steps：

The overall evaluation of a technical project interwell communication relationship and water-injection efficiency inversion result of the embodiment of the present invention establish a kind of new The fine water filling optimization method of stratified reservoir layering, to reduce inefficient water breakthrough flow and improve water-injection efficiency as target, by repeatedly It formulates oil-water well demixing injection-production scheme automatically for Optimization Solution, realizes that stratified reservoir dynamic matches the injection allocation design of production.It designs simultaneously Method calculates quickly, can accurately carry out saturation degree tracking and calculate output dynamic, gained seperated layer water injection scheme is preferably real Existing decreasing water cut and increasing oil effect increases oil nearly 50,000 side prediction year after practical oil reservoir optimization, provides for oil field injection and extraction scheme real-time optimization New approaches.

Above scheme is described in detail below.

In the step 1 of above-described embodiment, data foundation is adopted according to note first and is dynamically initially connected for grease between simulating well General character model, basic principle be non-homogeneous multi-layer oil reservoir is separated into it is a series of by conductivity (T between well_ijk) and connected volume (V_pijk) parameter characterization interwell communication unit, establish matter balance equation by object of connected unit：

Two kinds of production methods of liquid level pressure are determined in consideration, after solving to formula (1) implicit difference, can be obtained between interwell communication unit Flow

, can be theoretical using gram Lay displacement of reservoir oil after acquiring flow in another embodiment, by the grease inside connected unit Flowing is regarded as one-dimensional water drive oil problem, is clearing object with well point, is weighed shown in method such as formula (3) using upstream in connected unit Saturation degree is calculated, and then calculates grease output dynamic：

In formula：For the derivate of water content ratio for tracking the i-th Jing Chu from jth well in kth layer；Exist for jth well The derivate of water content ratio of kth layer；For from jth well stream to the i-th well from 0 to n moment zero dimension integrated flux.

In other embodiments, it is handling in actual oil reservoir calculating process, due to shutting down well, note adopts the oil fields such as conversion Measure leads to shaft bottom flow pattern drastic change, and this method can not be handled, can be using after fluid diversion in order to guarantee saturation degree stability Saturation degree backward tracing calculation method, it is primary all in accordance with the tracking of positive and direction during calculating saturation degree each time, take Its minimum value guarantees all to realize accurately calculating for saturation degree in the case where frequent fluid diversion, derivate of water content ratio calculation formula It is as follows：

In formula,In its k of ijWithIt represents from 0 moment until the i well stream at n ' moment and n moment is to j The zero dimension integrated flow of well；It represents from the n ' moment to the n moment from j well stream to the reversed zero dimension integrated flow of i well.

Isosaturation surface movement law inside connected unit is similar with the propagation law of wave, and grease saturation distribution calculates It can be converted into Riemannian problem, solved using shock theory.According to the inferior saturation discontinuity face of initial condition or so saturation degree situation, The mobile equation of grease isosaturation surface can be decomposed into shock wave, three kinds of forms of wave and complex wave are evacuated and solve.

By Berkeley displacement of reservoir oil theory it is found that any isosaturation surface fltting speed is equal to the point moisture content to its saturation degree Derivative.Such as some connected unit upstream well point j and downstream well point i, initial saturation degree is all larger than leading edge saturation degree, such as Fig. 2 institute Show.On the connected unit saturation distribution, the corresponding fltting speed in two saturation degree faces being arbitrarily designated distinguishes f_w'(S_w1) and f_w'(S_w2), the f known to its saturation degree_w'(S_w1) < f_w'(S_w2), so that two isosaturation surface distances are as time goes by, by Step is become estranged, therefore the saturation degree tracing computation belongs to evacuation wave.Before some connected unit upstream and downstream well point saturation degree is respectively less than Edge saturation degree is as shown in red lines, then the saturation degree face fltting speed above and below discontinuity surface is f_w'(S_w3) and f_w'(S_wc), by it F known to saturation degree_w'(S_w3) > f_w'(S_wc), two saturation degree faces can meet as time goes by, and then with f_w'(S_w3) speed Degree downstream promotes, and the region involved, saturation degree becomes S_w3, this seed type is shock wave.And complex wave be then by evacuation wave and Shock wave is composed, as shown in Fig. 2, upstream well point saturation degree is greater than leading edge saturation degree, downstream is less than leading edge saturation degree.From upstream At well point to oil water front, to evacuate wave, and oil water front to downstream well point is then shock wave.

Saturation degree point in connected unit when formula (4) can solve upstream well point water saturation greater than leading edge saturation degree very well Cloth problem, two species saturations distribution situation as shown in Figure 2.However in displacement process, if upstream well point saturation degree between When between irreducible water saturation and displacing front saturation degree, in the affiliated communicating downstream unit in the upstream well point, actually still with low It is promoted in leading edge saturation degree shock wave downstream, as shown in Figure 2, formula (4) can not describe this phenomenon, before aqueous Edge, which is broken through, can have hysteresis effect compared to actual conditions, aqueous after breakthrough to rise compared to reality faster.

For this problem, the technical solution of the present embodiment improves saturation degree method for tracing, for there are upstreams Well point j saturation degree is between leading edge saturation degree and irreducible water saturation, when downstream well point i is irreducible water saturation state, benefit With the method for solving of shock wave, fltting speed of the saturation degree face in connected unit is denoted as ν_ijk：

If interwell communication element length is L_ijk, upstream well point j is that coordinate remembers that saturation degree face moves Δ from upstream node The n time, so that well point i saturation computation results expression is as follows at this time：

Work as v_ijkΔ n > L_ijkWhen, then at the saturation degree face displacement to downstream well point i, at this point, moisture content can be pressed at the i of well point Flow rate equations are calculated, and otherwise well point output is still pure oil flow flowing.

Saturation computation process in summary, last improved well point saturation computation formula are as follows：

In order to verify to improved saturation degree method for tracing, using ECLIPSE simulator and the full of front and back is improved It is compared with degree method for tracing, devises " fourth " shape concept oil reservoir example, well location distribution is as shown in figure 3, totally 4 mouthfuls of wells, I1, I2 For two mouthfuls of well saliva wells, W1 closing well, the production of W2 well, wherein I1, I2 injection rate are 15 sides/day, 30 sides of W2 speed of production/day, grease Viscosity is respectively 20mPas, 1mPas, in-place permeability 1000mD, is produced 100 days altogether.As can be seen from Figure 4 after improving Saturation front method can be almost the same with the water breakthrough rule of business software oil reservoir grid model, it was demonstrated that improved side Method calculates more accurate.It is distributed by well location it is found that since I1 well waterflood front can preferentially arrive at W2 well, the waterflood front of subsequent I2 Arrive at W2 well.Therefore W2 well moisture content will appear two ladder-like transition, according to saturation computation method before improving, when W1 well meter When calculating well point saturation degree lower than aqueous leading edge saturation degree, then defaulting inside W1 to W2 connected unit is always pure oil flow, to lead W2 well water breakthrough is caused to break through lag, and improved method then can accurately simulate this problem.

To sum up improved saturation degree method for tracing can take into account low Water-cut Period saturation degree tracing computation.Saturation degree after improvement Method for tracing is still semi analytic, and pressure equation dimension is low, and whole process operation cost is far smaller than traditional reservoir modeling, and counts It calculates and stablizes, be able to satisfy live dynamic playback and program decisions in real time.

It, can be using the methods of random perturbation approximation method and Projected method to being initially connected in specific step 3 Property model carry out inverting, adjust the characteristic parameter of the initial connectivity modeling, specific method is on the books in disclosed document, Herein without being described in detail.

Using connectivity modeling after step 3 of the present invention optimization, each moment oil well delamination can be exported and produce liquid, oil-producing and well The information such as point coefficient are split in seperated layer water injection, interaction relationship between injection-production well in announcement oil reservoir plane, longitudinal direction.It is specific as follows,

In formula,The day oil-producing that the daily output liquid of respectively the i-th well kth layer splits point coefficient and kth layer is split point A point coefficient is split in water filling between coefficient, the i-th well and jth well in kth layer.

Every layer of water-injection efficiency of well can be further accurately calculated using above- mentioned information, it refers to water injection well in each layer to week Side oil well water supply displacement go out crude oil total amount and its in this layer of water injection rate ratio, i.e.,：

By formula (10) it is found that water injection well water-injection efficiency can reflect connected unit displacement shape in wet well and periphery oil well Condition, water-injection efficiency is smaller, and displacement goes out that crude oil is less, and water consumption rate is larger in the case of illustrating unit injection water, and around there may be excellent Potential barrier road causes certain non-available water to recycle.

Technical solution of the present invention proposes the basic think of that seperated layer water injection dynamic optimization is proposed based on connectivity modeling as a result, Think：It is compared in each interval water-injection efficiency with block water-injection efficiency by calculating last moment all water injection wells, water-injection efficiency Augmented injection is then carried out higher than block water-injection efficiency, on the contrary then reduced, the then base of the working system of water injection well after the adjustment On plinth, a point coefficient is split to the water filling of each oil well according to each well, the working system of oil well is adjusted, to reduce inefficient Water drive direction flow improves block water-injection efficiency.After producing a period of time, water-injection efficiency is reappraised, adjustment note is adopted again Working system realizes that oil reservoir injection adopts policy dynamic optimization.

Seperated layer water injection dynamic optimization method specific workflow is as follows：

(1) it obtains the block at oil reservoir development history last m moment to be averaged water-injection efficiency, i.e. block oil production is infused divided by block Water, calculation formula are as follows:

(2) each layer water-injection efficiency of well is evaluated, if water injection well water-injection efficiency Wef_ik< Wef so then needs to drop Note, on the contrary then augmented injection, each layer injection rate of well after specific optimizationFor：

Usually consider the constraint of the actual block condition of production, the adjustment of injection-production well liquid measure is limited, generally recommendation η_min=-0.5, η_max =0.5, i.e. liquid measure amplitude fluctuation is no more than 50%, Wef_minWith Wef_maxFor water-injection efficiency in individual well is minimum and peak, α is to refer to It is several, it is general to recommend value 2, to constrain η_ikSize.

(3) after the completion of seperated layer water injection scheme optimization, according to knots modification and last moment after the optimization of well water injection rate to surrounding Point coefficient of splitting of oil well optimizes the Liquid output of oil-producing well.It is as follows with calculation formula is produced：

So far, oil-water well next stage note adopts prioritization scheme and has formulated completion, brings stratified reservoir inter well connectivity model into In carry out simulation a period of time, recalculate and split point coefficient and water-injection efficiency, repeat above step (1)~(3), formulate again Note adopts prioritization scheme, realizes dynamic production and injection proration optimization.

The paraphrase respectively met in above-mentioned formula is as follows：

A_ik--- i-th mouthful produces liquid in kth layer and splits component, f；A_ijk--- i-th mouthful of water injection well is filled the water to jth mouth well in kth layer Split a point coefficient, f；A_oik--- i-th mouthful of producing well splits component, f in kth layer oil-producing；C_tk--- the system compressibility of kth layer, MPa^-1；f_w--- moisture content, f；f_w′(s_wijk) --- the derivate of water content ratio of the i-th Jing Chu is tracked from jth well in kth layer；f_w′ (S_wjk) --- derivate of water content ratio of the jth well in kth layer；F_vijk--- it is accumulative in the zero dimension of kth layer from jth well stream to the i-th well Flow；L_ijk--- i well and j well are in kth layer distance；N_l--- oil reservoir number；N_w--- injection-production well number；N_I--- water injection well number；p_i, p_j--- the average pressure in the i-th well and jJing draining area, MPa；q_i--- the i-th well yield, injection is positive, output is negative, m³/ s；q_ijk--- kth layer i well and j interwell communication unit inner stream flow, m³/s；S_wc--- irreducible water saturation, f；S_wf--- water drive Leading edge saturation degree, f；S_wijk--- the water saturation of kth layer i well and j interwell communication unit, f；S_wik--- the i-th well is in kth The water saturation of layer, f；T --- time, s；T_ijk--- the average conduction rate between kth layer i well and j well, m³/(s·MPa)； V_ijk--- kth layer i well and j interwell communication volume, m³；V_pik--- the drainage volume of the i-th well of kth layer, m³；ν_ijk--- kth layer i Saturation degree face fltting speed, m/s inside well and j well connected unit；The water-injection efficiency of Wef --- block, f；Wef_ik--- i-th mouthful Water injection well is in kth layer water-injection efficiency；Wef_max, Wef_min--- minimum and peak efficiency in current water injection well, f；α --- produce liquid tune Whole index；η_ik--- liquid measure regulation coefficient of i-th mouthful of water injection well in kth layer；η_max, η_min--- above and below liquid measure regulation coefficient Limit.Subscript：M --- last moment optimizing phase；N --- current time, s；N ' --- last moment, s.Subscript：I, j --- well Serial number；K --- oil reservoir serial number.

It is illustrated below by way of specific embodiment 2.

Embodiment 2 by Research Numerical Simulation Techique construct a grid dividing be 21 × 21 × 2 there are hypertonic bands Positive rhythm reservoir, sizing grid DX=DY=20m, DZ=10m, 9 mouthfuls of wells of oil reservoir, 5 note 4 are adopted, separated injection and combined recovery, permeability fields point Cloth is as shown in Fig. 5 a, Fig. 5 b, and initial reservoir saturation degree is 0.2, and grease viscosity is respectively 1.0 and 20.0mPas.Using maturation Digital-to-analogue software Production development simulation trial carried out to the oil reservoir, the simulated production time is 900 days, and oil reservoir integrated injection is adopted flat Weighing apparatus, final block moisture content reach 82.0%.

Preceding 900 days Production developments to be simulated using stratified reservoir inter well connectivity, the primary time-consuming of operation is 0.22 second, Primary time-consuming is simulated 20 seconds compared to ECLIPSE, and speed-raising is nearly a hundred times.Then, by connectivity modeling simulation dynamic and ECLIPSE model Production development carry out history matching, calculating is optimized using random perturbation approximate gradient algorithm in fit procedure, by 50 The convergence of multi-Step Iterations Optimal Fitting.Fig. 6 a, Fig. 6 b are that block and individual well Production development fitting result, fitting effect are preferable.Fig. 7 a, Fig. 7 b conductivity inversion result between oil reservoir well, a numerical value in bracket are conductivity, m³/ (MPas-1), second value To control volume, 104m³, wherein one be high infiltration strip, one be moderate connection, another for connecting degree relatively It is low, be consistent substantially with the permeability fields of concept oil reservoir, demonstrate model parameter after inverting can hold oil reservoir main geologic it is special Sign.Each layer well splits a point coefficient results after inverting shown in Fig. 8 a, Fig. 8 b, and W5 well splits point size point to P1~P4 well in first layer Not Wei 0.48,0.15,0.27,0.10, W5 well the second layer to point size of splitting of P1~P4 well be 0.27,0.29,0.30, 0.14, wherein the triangular arrowheads the big, shows to split shunt volume between the oil-water well bigger, and water flow is split point mainly along high infiltration strip stream The dynamic Analysis on Mechanism with practical geologic feature is consistent, it can be seen that connectivity modeling can characterize practical oil reservoir stream after inverting The influence that dynamic rule and high infiltration strip flow grease.

According to formula (11), (12), the water-injection efficiency of each layer of oil reservoir individual well and block to the last production moment is counted It calculates, as shown in figure 9, the water-injection efficiency of block is 0.17, i.e., block injects side's water about 0.17 side of output oil.Due to Geological Model Second layer permeability grade is compared higher with first layer in type, and for water flow more easily along high infiltration strip channelling, non-available water circulation is existing As serious, second layer individual well water-injection efficiency is caused to be decreased obviously, information obtained by inverting is consistent with Analysis on Mechanism.

According to seperated layer water injection dynamic optimization workflow, working system is adopted to the note of oil reservoir according to formula (13)~(15) and is protected It holds injection-production ratio to optimize, optimize 540 days altogether, point 3 regulation steps, every step 180 days, working system such as Figure 10 a, figure after optimization Shown in 10b, since first layer water injection well efficiency is all larger than block water-injection efficiency, first layer for the first time infuse by regulation in prioritization scheme Water is significantly increased, and water injection well is below block water-injection efficiency in second layer water-injection efficiency, and then injection rate substantially reduces the second layer, To reduce the flow of second layer high permeability formation non-available water circulation, it is suppressed that highly water intensive drive improves interlayer and note adopts contradiction, real The uniform Flooding Efficiency of stratified reservoir is showed.For oil well, it can be seen that due to the remaining main enrichment of the presence of block high infiltration strip Between producing well P2 and P4, this corresponding Liquid output of two mouthfuls of wells is substantially improved, and P1 keeps stablizing, the P3 then steady water of descending liquid.

Block index after optimization is as shown in figure 11, can be seen that from block production rate curve by being layered fine note Water optimizes and revises working system, plays the effect of decreasing water cut and increasing oil.It is excellent by maintaining within 540 days the working system adjustment of injection-production ratio Block after change tires out oil-producing and increases 3660 sides, and moisture content then has dropped 0.5% compared to before optimization.By the work after optimization System brings ECLIPSE reservoir simulator into, and compared the reservoir model remaining oil feelings produced always according to original injection-production program Condition, as shown in Figure 12 a, Figure 12 b, Figure 12 c, Figure 12 d.Two layers of remaining oil is decreased obviously after optimization, is illustrated by utilizing company The water filling of general character method simulation is split point and water-injection efficiency information, realizes only to pass through to infuse under conditions of maintaining flooding pattern and adopts liquid Stratified reservoir Tapping Residual Oil is realized in amount adjustment.

Embodiment 3 is related to one piece of overseas multilayer water controlled field, and 12 mouthfuls of certain test block well, 6 note 6 is adopted, and carries out simulation calculating, The oil reservoir is always divided into three layers, and separated injection and combined recovery, by 2011 days before automatic history matching, inverse model parameter exported the last moment Water filling is split point and the information such as oil displacement efficiency, and then is adjusted to oil reservoir development scheme.

Block tires out oil-producing and the history matching of moisture content Production development as shown in Figure 13 a, Figure 13 b, which shows oil reservoir grease Dynamic has obtained preferable fitting, and fitting correlation coefficient reaches 95% or more.Meanwhile according to the parameter field after inverting as a result, mould Point information and water-injection efficiency are split in the quasi- each layer water filling for outputing the last moment, and wherein block water-injection efficiency is 0.34, is such as schemed 14a, Figure 14 b, shown in Figure 14 c and Figure 15.Seperated layer water injection dynamic optimization work in 540 days has been carried out to the block according to the information, It is divided into note three times and adopts regulation, every step 180 days, the note amount of the adopting working system after optimization is as shown in Figure 16 a, Figure 16 b, it can be seen that And the block water filling effect biggish water injection well injection rate of deviation has obtained significantly adjusting, similar well, and injection rate fluctuation is little. It is configured by seperated layer water injection dynamic optimization three times, the production target after block optimization is as shown in figure 17, and prioritization scheme plays bright Aobvious decreasing water cut and increasing oil effect, moisture content have dropped 1.6%, and tired oil-producing increases 7.05 ten thousand sides.

To sum up, oil reservoir development production optimization method provided by the invention has the advantages that：

(1) improved saturation degree method for tracing preferably features each species saturation face progradation in connected unit, Grease flowing law complicated in stratified reservoir can be more preferably simulated, the low water-cut stage especially in；

(2) model built achieves preferable automatic history matching effect, gained inter well connectivity parameter and conceptual model Geologic feature matches, and effectively discloses inefficient water breakthrough direction, and seperated layer water injection optimization method improves note and adopts contradiction, improves water filling Utilization efficiency, moisture content decline 1.6% after practical oil reservoir optimization, tires out oil-producing and increases by 7.05 ten thousand sides；

(3) essence of the seperated layer water injection optimization method proposed is that situation and flowing pass are connected between quantitative recognition injection-production well On the basis of system, by reducing the flow distribution of inefficient water drive direction, inhibits highly water intensive expelling water to alter, realize the effect of decreasing water cut and increasing oil Fruit.This method calculates quickly, does not depend on complex geologic model-building, can blend with the existing real-time water injection technology technology of layering, into one Step meets the oil field development production decision demand of big data.

Figure 18 is that the structure for the oil reservoir development production optimization system based on connectivity method that the embodiment of the present invention 4 provides is shown It is intended to, as shown in figure 18, including：

Comprehensive interwell communication relationship and water-injection efficiency inversion result establish a kind of new stratified reservoir in above-described embodiment The fine water filling optimization method of layering is solved using reducing inefficient water breakthrough flow and improving water-injection efficiency as target by iteration optimization It is automatic to formulate oil-water well demixing injection-production scheme, realize that stratified reservoir dynamic matches the injection allocation design of production.Design method calculates fast simultaneously Speed can accurately carry out saturation degree tracking and calculate output dynamic, and gained seperated layer water injection scheme preferably realizes decreasing water cut and increasing oil Effect increases oil nearly 50,000 side prediction year after practical oil reservoir optimization, provides new approaches for oil field injection and extraction scheme real-time optimization.

In preferred embodiment, the computing module is specifically included：

Saturation degree method for tracing is improved in above preferred embodiment, is solved in displacement process when upstream well When point saturation degree is between irreducible water saturation and displacing front saturation degree, the saturation degree method for tracing of the prior art is difficult to standard The problem of really solving saturation distribution state, saturation computation method is more accurate, and each layer water-oil phase produces at the well point of generation Dynamic indicator is also more accurate out, further improves the effect of oil reservoir development.

In another preferred embodiment, the optimization module is specifically included：

Wherein Wef is the water-injection efficiency of block, Wef_ikIt is i-th mouthful of water injection well in kth layer water-injection efficiency, Wef_max, Wef_min For minimum and peak efficiency in current water injection well, α is to produce liquid adjustment index, η_ikLiquid measure for i-th mouthful of water injection well in kth layer adjusts Coefficient, η_max, η_minFor the bound of liquid measure regulation coefficient, specific η_min=-0.5, η_max=0.5, α value are that 2, m is optimization Last moment in stage；

Above preferred embodiment can reduce inefficient water drive direction flow, improve block water-injection efficiency, while realizing oil reservoir Note adopts policy dynamic optimization.

Although the embodiments of the present invention has been shown and described above, it is to be understood that above-described embodiment is example Property, it is not considered as limiting the invention, those skilled in the art within the scope of the invention can be to above-mentioned Embodiment is changed, modifies, replacement and variant.

Claims

1. a kind of oil reservoir development production optimization method based on connectivity method, which is characterized in that include the following steps：

Step 1, it establishes for simulating between well grease dynamically initial connectivity modeling, the initial connectivity modeling is by oil reservoir injection Extraction system is separated into the connected unit between well and well, and each connected unit is by conductivity and control two characteristic parameter tables of volume Sign；

Step 2, equation is promoted to carry out between well point pressure and well based on connected unit and by matter balance equation and waterflood front Flow rate calculation, and saturation degree tracking is carried out, it generates the current oil aquatic products at well point and goes out dynamic indicator；

Step 3, go out dynamic indicator to the current oil aquatic products using the history grease output dynamic indicator of oil reservoir to be intended automatically Inverting is closed, the characteristic parameter of the initial connectivity modeling is adjusted, and generates the connectivity modeling after optimization；

Step 4, last moment all water injection wells are calculated in the water-injection efficiency of each interval by connectivity modeling after the optimization, And the block for obtaining the oil reservoir development history last moment is averaged water-injection efficiency, it is flat according to the water-injection efficiency of each interval and the block The comparison result of equal water-injection efficiency carries out the first suboptimization to the current development and production method of water injection well；Using the water filling after optimization Well development production method calculates water filling of each well to surrounding oil well and splits a point coefficient, and splits point coefficient to optimization according to the water filling Water injection well development and production method afterwards carries out the second suboptimization.

2. the oil reservoir development production optimization method according to claim 1 based on connectivity method, which is characterized in that described Step 2 specifically includes following steps：

S201 establishes matter balance equation using the connected unit as object：

Wherein, t is the time, and i, j are well serial number, and k is oil reservoir serial number, N_lFor oil reservoir number, N_wFor injection-production well number, p_i, p_jFor the i-th well and Average pressure in jJing draining area, q_iFor the i-th well yield, injection is positive, output is negative, C_tkFor the comprehensive compression system of kth layer Number, V_pikFor the drainage volume of the i-th well of kth layer；

S202 carries out implicit difference solution to the matter balance equation using liquid and level pressure two ways is determined, connects between generation well Flow between logical unit

Wherein, n is current time, q_ijkFor kth layer i well and j interwell communication unit inner stream flow, p_i, p_jFor the i-th well and jth well Average pressure in draining area, T_ijkFor the average conduction rate between kth layer i well and j well；

S203, when upstream well point j saturation degree is between leading edge saturation degree and irreducible water saturation and downstream well point i is irreducible water When saturation degree state, using the method for solving of shock wave, fltting speed of the saturation degree face in connected unit is denoted as ν_ijk：

If interwell communication element length is L_ijk, saturation degree face is Δ from the time of upstream node j movement in connected unit N, at this time well point i saturation computation result be：

Wherein, n is current time, and n ' is last moment, f_wFor moisture content, f '_w(s_wijk) it is to track the from jth well in kth layer The derivate of water content ratio of i Jing Chu, f '_w(S_wjk) it is derivate of water content ratio of the jth well in kth layer, F_vijkTo exist from jth well stream to the i-th well The zero dimension integrated flow of kth layer, L_ijkIt is i well and j well in kth layer distance, S_wijkFor kth layer i well and j interwell communication unit Water saturation, S_wikWater saturation for the i-th well in kth layer, V_ijkFor kth layer i well and j interwell communication volume, ν_ijkFor Saturation degree face fltting speed inside kth layer i well and j well connected unit.

3. the oil reservoir development production optimization method according to claim 1 or 2 based on connectivity method, which is characterized in that The step 4 is specially：

S401, by step 3 optimize after connectivity modeling calculate last moment all water injection wells each interval water-injection efficiency, And the block for acquiring the oil reservoir development history last moment is averaged water-injection efficiency；

S402 evaluates each interval water-injection efficiency of water injection well, if the water-injection efficiency of water injection well is less than the block and averagely infuses Water efficiency then carries out drop note, otherwise carries out augmented injection, and generate the oil reservoir development production method after the first suboptimization, specific to optimize Each layer injection rate of water injection well afterwardsFor：

Wherein Wef is the water-injection efficiency of block, Wef_ikIt is i-th mouthful of water injection well in kth layer water-injection efficiency, Wef_max, Wef_minTo work as Minimum and peak efficiency in preceding water injection well, α are to produce liquid adjustment index, η_ikLiquid measure for i-th mouthful of water injection well in kth layer adjusts system Number, η_max, η_minFor the bound of liquid measure regulation coefficient, m is last moment optimizing phase；

S403 calculates water filling of the last moment each well to surrounding oil well using the oil reservoir development production method after the first suboptimization A point coefficient is split, and point coefficient is split according to the water filling, the second suboptimization is carried out to oil reservoir development production method, it is specific with production optimization Calculation formula is：

Wherein A_ijkA point coefficient is split in the water filling of kth layer to jth mouth well for i-th mouthful of water injection well, m is last moment optimizing phase, η_ikFor I-th mouthful of water injection well kth layer liquid measure regulation coefficient,For each layer injection rate of water injection well after optimization, N_IFor water injection well number, q_iFor I-th well yield, injection is positive, output is negative；

S404 after simulating preset time using the oil reservoir development production method after the second suboptimization, repeats step S401-S403, Current oil reservoir development and production method is optimized again.

4. the oil reservoir development production optimization method according to claim 3 based on connectivity method, which is characterized in that η_min =-0.5, η_max=0.5, α value are 2.

5. a kind of oil reservoir development production optimization system based on connectivity method, which is characterized in that including：

Model building module, for establishing for simulating between well grease dynamically initial connectivity modeling, the initial connectivity Oil reservoir injection and extraction system is separated into the connected unit between well and well by model, and each connected unit is by conductivity and control volume two A characteristic parameter characterization；

Computing module, for promoting equation to carry out well point pressure based on connected unit and by matter balance equation and waterflood front The flow rate calculation between well, and saturation degree tracking is carried out, it generates the current oil aquatic products at well point and goes out dynamic indicator；

Inverting module goes out dynamic indicator to the current oil aquatic products for the history grease output dynamic indicator using oil reservoir and carries out Automatic Fitting inverting, adjusts the characteristic parameter of the initial connectivity modeling, and generates the connectivity modeling after optimization；

Optimization module, for calculating last moment all water injection wells in the note of each interval by the connectivity modeling after the optimization Water efficiency, and the block for obtaining the oil reservoir development history last moment is averaged water-injection efficiency, according to the water-injection efficiency of each interval and institute It states the be averaged comparison result of water-injection efficiency of block and the first suboptimization is carried out to the current development and production method of water injection well；Using optimization Water injection well development and production method afterwards calculates water filling of each well to surrounding oil well and splits a point coefficient, and is split according to the water filling and point be Water injection well development and production method after several pairs of optimizations carries out the second suboptimization.

6. the oil reservoir development production optimization system according to claim 5 based on connectivity method, which is characterized in that described Computing module specifically includes：

First computing unit is determined liquid and level pressure two ways for use and is asked matter balance equation progress implicit difference Solution generates the flow between interwell communication unitFor：

Second computing unit, for when upstream well point j saturation degree is between leading edge saturation degree and irreducible water saturation and downstream When well point i is irreducible water saturation state, using the method for solving of shock wave, by fltting speed of the saturation degree face in connected unit It is denoted as ν_ijk：

And interwell communication element length is set as L_ijk, saturation degree face is from the time of upstream node j movement in connected unit Δ n, generating well point i saturation computation result is：

7. the oil reservoir development production optimization system according to claim 5 or 6 based on connectivity method, which is characterized in that The optimization module specifically includes：

Third computing unit, for calculating last moment all water injection wells in the note of each interval by the connectivity modeling after optimization Water efficiency, and the block for acquiring the oil reservoir development history last moment is averaged water-injection efficiency；

First optimization unit, for evaluating each interval water-injection efficiency of water injection well, if the water-injection efficiency of water injection well is less than institute It states block to be averaged water-injection efficiency, then carries out drop note, otherwise carry out augmented injection, and generate the oil reservoir development producer after the first suboptimization Method, each layer injection rate of water injection well after specific optimizationFor：

Second optimization unit, for calculating last moment each well to week using the oil reservoir development production method after the first suboptimization A point coefficient is split in the water filling for enclosing oil well, and is split point coefficient according to the water filling and carried out the second suboptimization to oil reservoir development production method, It is specific to be with production optimization calculation formula：

Wherein A_ijkA point coefficient is split in the water filling of kth layer to jth mouth well for i-th mouthful of water injection well, m is last moment optimizing phase, η_ikFor I-th mouthful of water injection well kth layer liquid measure regulation coefficient,For each layer injection rate of water injection well after optimization, N_IFor water injection well number, q_iFor I-th well yield, injection is positive, output is negative.

8. the oil reservoir development production optimization system according to claim 7 based on connectivity method, which is characterized in that η_min =-0.5, η_max=0.5, α value are 2.