CN107291667A - A kind of interwell communication degree determines method and system - Google Patents
A kind of interwell communication degree determines method and system Download PDFInfo
- Publication number
- CN107291667A CN107291667A CN201610203017.6A CN201610203017A CN107291667A CN 107291667 A CN107291667 A CN 107291667A CN 201610203017 A CN201610203017 A CN 201610203017A CN 107291667 A CN107291667 A CN 107291667A
- Authority
- CN
- China
- Prior art keywords
- well
- mrow
- percent continuity
- water injection
- msub
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
Method is determined the invention discloses a kind of interwell communication degree, including:Inter well connectivity determines step;Percent continuity determines step, by founding mathematical models and specifying constraint, determines the percent continuity between connected well;Evaluation criterion establishment step, according to the percent continuity and statistical sample distribution, obtains the evaluation criterion of connecting degree between connected well;Interwell communication degree determines step, and interwell communication degree is determined according to the percent continuity and the evaluation criterion.The present invention can quick and precisely ask for interwell communication coefficient so that it is more accurate that interwell communication degree is calculated.
Description
Technical field
The invention belongs to oil reservoir development technical field, specifically, it is related to a kind of interwell communication degree and determines method
And system.
Background technology
Fracture-pore reservoir is the main Oil Reservoir Types of northwest China, with continuing for exploitation, is initially entered comprehensively
Waterflood stage, its subject matter faced be flooding pattern how to set up and water logging after flooding pattern how
Adjustment.For problem above, the research to fracture-pore reservoir inter well connectivity is broadly divided into two aspects:One is
Whether the connection situation understood fully between well, i.e., connect between injection-production well or even between producing well, mainly solve flooding pattern
The problem of how setting up;Another is clearly to know that connection between injection-production well is strong and weak, is to solve water logging
The problem of water drive is adjusted afterwards.At present, the method that inter well connectivity quantitative assessment technology is mainly used includes following several
Kind.
Grey relational grade evaluate inter well connectivity general principle be:The change of a bite oil well output, meeting in well group
There is close association type or interference with the water injection rate of a certain water injection well around, this property has uncertainty
(namely ash), and grey relational grade analysis can just solve Similar Problems.But, this method is deposited
Problem is determined in the degree of association, and many scholars have found that the computational methods of the degrees of association can not reflect two variables exactly
Correlation, while the value of resolution ratio is a constant interval, value difference can cause analysis result difference very
Greatly.
When reservoir engineering method evaluates interwell communication degree, often using multiple Production development indexs, work as difference
When contradiction occurs in metrics evaluation result, Fuzzy Evaluation method just shows its advantage.Fuzzy mathemacical evaluation
The method for evaluating inter well connectivity summarizes the numerical value of each different evaluation index by constructing membership function, draws
One comprehensive evaluation result.But, membership function has in many types, current oil field using being all most simple
Single trapezoid formula.Whether the trapezoid formula is applied to evaluate connective, need to be proved.Another aspect weight
Calculating there is also problem, analytic hierarchy process (AHP) equal weight determines that method is all the subjective determination by engineer experience
Method, accuracy is low.
Method for numerical simulation inverting inter well connectivity technology is to be based on static data, the porosity of such as oil reservoir, infiltration
The parameter fields such as rate, saturation degree, reservoir thickness, set up geological model of oil accumulation;Then the different individual wells of oil reservoir are fitted
The production histories such as bottom pressure, aquifer yield, oil production, verify geological model;Finally on the basis of this model
The percent continuity for obtaining oil reservoir is calculated, connecting degree is evaluated.But, Method for Numerical needs numerous sound state moneys
Material, will collect the authentic data for just meeting oil reservoir extremely difficult.In addition, fracture-pore reservoir modeling method is deposited in itself
In problem, the uncertainty of model is very big, and dynamic history match is also a cumbersome, time-consuming job, itself
There is also multi-solution.Therefore, method for numerical simulation inverting determines that the method reliability of inter well connectivity is not high.
Multiple linear regression derives from the thought of balanced flood, it is believed that the water injection rate of all wells is equal to all in well group
The production liquid of well adds the intrusion volume of bottom water (or side water).The estimation of weight coefficient is asked for using least square method
Value, estimate now is just the connecting degree between well.Scholars before use multiple linear regression analysis method meter
When calculating inter well connectivity, the feature mainly with two aspects:One is to apply balanced flood, interwell communication system
Showing of digital display is expressed as the weight of different producing well Liquid outputs;Two (curve matching and are forced using least square method
Nearly method), minimum by the quadratic sum of error, partial derivative is zero, solves system of linear equations and obtains percent continuity.
But, for substantial amounts of Production development data, calculated using least square method, huge data matrix can be caused,
Add the difficulty of solution.In addition, the percent continuity solved is easily negative value, with actual oil reservoir physical significance not
Clothes, accuracy needs to verify.
The content of the invention
To solve problem above, method and system are determined the invention provides a kind of interwell communication degree, for fast
Speed accurately asks for interwell communication coefficient so that it is more accurate that interwell communication degree is calculated.
Method is determined there is provided a kind of interwell communication degree according to one embodiment of present invention, including:
Inter well connectivity determines step;
Percent continuity determines step, by founding mathematical models and specifying constraint, determines the company between connected well
Logical coefficient;
Evaluation criterion establishment step, according to the percent continuity and statistical sample distribution, is obtained between connected well
The evaluation criterion of connecting degree;
Interwell communication degree determines step, and interwell communication journey is determined according to the percent continuity and the evaluation criterion
Degree.
According to one embodiment of present invention, percent continuity determines that step further comprises:
Based on situation is adopted a note in injection-production well group, setting up the balanced flood mathematical modeling with constraints, institute more
Stating mathematical modeling includes the percent continuity;
Percent continuity is primarily determined that using monte carlo method;
The percent continuity primarily determined that is normalized to obtain final percent continuity.
According to one embodiment of present invention, the step of primarily determining that percent continuity using monte carlo method enters one
Step includes:
Reality in acquisition well group connected relation, well group in the production liquid maximum fluctuation value of each producing well and water injection time section
Border water injection rate data;
Using gauss of distribution function, to connect the production liquid maximum fluctuation value of each producing well in well group as constraints,
One group of random number is produced for each percent continuity;
The corresponding random number of each percent continuity is substituted into the mathematical modeling, water injection rate valuation is calculated;
According to water injection rate valuation and actual water injection rate data, obtain and the actual immediate water injection rate of water injection rate data
And corresponding random number;
Using this group of random number as the axle value of the Gaussian Profile of corresponding well, the new of corresponding well percent continuity is produced
One group of random number, return to water injection rate valuation step, until water injection rate valuation is equal with actual water injection rate, it is determined that
Percent continuity.
According to one embodiment of present invention, the mathematical modeling is:
Wherein,Represent the water injection rate of i-th mouthful of water injection well;qj(n) producing well around water injection well is represented
Liquid output;β0jNon-equilibrium constant is represented, physical significance is the water influx of bottom water;βijRepresent jth mouthful production
The percent continuity of well and i-th mouthful of water injection well;J represents the producing well quantity in i-th mouthful of water filling well group,
Constraints is:
βij∈(wj-d,wj+d)
Wherein, wjThe zero dimension crest value of the creation data of jth mouthful producing well is represented, d represents that is permitted a change
The threshold value of change.
According to one embodiment of present invention, evaluation criterion establishment step further comprises:
Count different injection-production well groups, different notes and adopt the frequency that the percent continuity in period occurs, and draw cumulative frequency
Distribution map;
Read the percent continuity of the different accumulative perceptions of correspondence respectively in the cumulative frequency distributing graph;
The evaluation criterion of interwell communication degree is obtained according to the corresponding percent continuity of different accumulative perceptions.
According to one embodiment of present invention, inter well connectivity determines that step includes sentencing by tracer test result
Disconnected inter well connectivity.
According to one embodiment of present invention, inter well connectivity determines that step includes passing through oil well production dynamic data
Inter well connectivity is determined, including:
The Production development data of each producing well are obtained from creation data table;
According to the production curve situation of change of producing well during water filling, judge connective between injection-production well, wherein,
The interior pressure of producing well, production liquid occur in that obvious rising, or the moisture content of producing well go out during such as water filling
Show significant change, then illustrate connection between injection-production well.
According to one embodiment of present invention, if during water filling the production curve of producing well change it is same
When, working system changes, then it is assumed that do not connected between injection-production well.
According to one embodiment of present invention, place is normalized to originally determined percent continuity based on following formula
Reason:
Wherein, n is well quantity, and m is the quantity of oil well in i-th mouthful of water injection well well group, βijRepresent jth mouthful
The percent continuity of producing well and i-th mouthful of water injection well.
According to another aspect of the present invention, additionally provide a kind of interwell communication degree and determine system, including,
Inter well connectivity determining module, for determining the connectedness between well;
Percent continuity determining module, by founding mathematical models and specifying constraint, determines the company between connected well
Logical coefficient;
Evaluation criterion sets up module, according to the percent continuity and statistical sample distribution, obtains between connected well
The evaluation criterion of connecting degree;
Interwell communication degree determining module, interwell communication journey is determined according to the percent continuity and the evaluation criterion
Degree.
Beneficial effects of the present invention:
The present invention is using taking a little at random near Gaussian Profile axle value so that axle value Step wise approximation actual value is tried to achieve
Percent continuity, compared with the general uniform thought taken a little, accelerates the velocity of approch to actual value;The opposing party
Face, according to the characteristics of actual production dynamic fluctuation, row constraint is entered to percent continuity, compared with least square method,
Be not in the situation that percent continuity is negative value so that result of calculation is more accurate.
Other features and advantages of the present invention will be illustrated in the following description, also, partly from specification
In become apparent, or by implement the present invention and understand.The purpose of the present invention and other advantages can pass through
Specifically noted structure is realized and obtained in specification, claims and accompanying drawing.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment
Or the accompanying drawing required in description of the prior art does simple introduction:
Fig. 1 is method flow diagram according to an embodiment of the invention;
Fig. 2 is Monte carlo algorithm schematic flow sheet according to an embodiment of the invention;
Fig. 3 is cumulative probability distribution schematic diagram according to an embodiment of the invention;
Fig. 4 is that INJ-WELL6 notes according to an embodiment of the invention adopt response curve schematic diagram;
Fig. 5 is percent continuity cumulative probability distribution schematic diagram according to an embodiment of the invention;And
Fig. 6 is that a point index contrast analysis chart is split in water filling according to an embodiment of the invention.
Embodiment
Describe embodiments of the present invention in detail below with reference to drawings and Examples, whereby to the present invention how
Application technology means solve technical problem, and reach the implementation process of technique effect and can fully understand and real according to this
Apply.As long as it should be noted that conflict is not constituted, in each embodiment and each embodiment in the present invention
Each feature can be combined with each other, and the technical scheme formed is within protection scope of the present invention.
The present invention is based on balanced flood principle, regards the duct in oil reservoir between water injection well, producing well and well as one
Individual complete system, wherein, the injection rate of water injection well can regard the excitation of system as, and the Liquid output of producing well is should
The response output of system, the total yield output of producing well is equal to the water influx that total water injection rate in oil reservoir adds bottom water.This
Invention is special using covering using multiple linear regression analysis method founding mathematical models from analysis interwell communication relation
Carlow random distribution asks for interwell communication coefficient, finally evaluates the power of connecting degree between injection-production well.
It is the method flow diagram according to one embodiment of the present of invention as shown in Figure 1, comes below with reference to Fig. 1 to this
Invention is described in detail.
It is that step S110 inter well connectivities determine step first.Specifically, in this step, generally use with
Lower two methods determine the connectedness between well.One of which method is the injection-production well for doing tracer test
Group, can directly judge the connectedness between well using tracer test result.
Another method is the injection-production well group for not doing tracer test, can be by contrasting water filling
Afterwards, whether the production curve of surrounding producing well changes to judge whether connect between well.The finger that this method judges
The data such as the miscarriage pressure of producing well in injection-production well group, liquid, oil-producing, aqueous are designated as, these data sources are in oil field
Production development tables of data.Determine that the detailed process of inter well connectivity is as follows using this method.
The Production development data of each producing well are obtained from creation data table first.The Production development packet of acquisition
Include water injection well water injection time, water injection rate, the pressure of producing well, production liquid, aqueous etc..
Then analysis water filling during producing well production curve situation of change, according to the situation of change of production curve
Judge the connectedness between injection-production well.Specifically, the interior pressure of producing well, production liquid are occurred in that substantially during such as water filling
Rise;Or, the moisture content of producing well occurs in that significant change, such as substantially rising, illustrates that water filling causes water logging,
Such as it is decreased obviously, illustrates that water filling causes oil displacement efficiency.Appearance case above is connection between all illustrating injection-production well.
In addition, it should be noted that the influence due to working system change to producing well production curve is very big, therefore,
If the production curve of producing well changes during water filling, working system (oil nozzle, stroke and jig frequency
Deng) also changed, then it is assumed that the change of production curve is that caused by the influence of itself working system, note is adopted
Do not connected between well, should give and exclude.
Step is determined followed by step S120 percent continuities, by founding mathematical models and specifying constraint,
Determine the percent continuity between connected well.Specifically, percent continuity determines that step further comprises following steps.
First, based on situation is adopted a note in injection-production well group, set up the balanced flood mathematical modeling with constraints more,
The mathematical modeling includes percent continuity.In view of actual field use, the mathematical modeling is expressed as follows:
Wherein,Represent the water injection rate of i-th mouthful of water injection well;qj(n) producing well around water injection well is represented
Liquid output;β0jNon-equilibrium constant is represented, physical significance is the water influx of bottom water;βijRepresent jth mouthful production
The percent continuity of well and i-th mouthful of water injection well;J represents the producing well quantity in i-th mouthful of water filling well group.
According to actual conditions, multivariate model needs to change in the range of with actual physical meaning, beyond the model
Practical significance will not had by enclosing, and setting constraints is as follows:
βij∈(wj-d,wj+d) (2)
Wherein, wjThe zero dimension crest value of the creation data of jth mouthful producing well is represented, i.e. (jth during water filling
The maximum crest value of the Liquid output of mouth producing well)/(sum of the maximum crest value of all producing wells in well group),
It is reflected as the degree of jth mouthful producing well fluctuation, and d represents a threshold value for permitting change, and default value is 0.1.
The physics meaning is exactly constraint betaijChange in the range of creation data fluctuation, it is impossible to too general.
Next percent continuity is primarily determined that using monte carlo method.Monte carlo method is come using random number
The method for solving computational problem, main thought is a series of random points of production Gaussian Profile, and random point is distinguished
Make comparisons, can obtain and the immediate random point of desired value, then again using new random point as Gauss with desired value
The axle point of distribution, produces one group of new random point and is compared again with desired value, so go down repeatedly, instruct random point
Untill equal with target point.
It should be noted at following 2 points using above method:One is the probability of happening maximum near Gaussian Profile axle point,
Produce one group of new random point as the average value of Gaussian Profile using the random value close to desired value every time, understand with
Machine point is concentrated near advantageous point so that the speed for approaching desired value is accelerated;Two be that percent continuity is bent using production liquid
Line fluctuation characteristic prevents random number from exceeding the scope of actual physical meaning as constraint.Using computer iterations
Percent continuity is solved, specific steps are as shown in Figure 2.
Well group connected relation is obtained first, in well group in the production liquid maximum fluctuation value and water injection time section of each producing well
Actual water injection rate data.(or tracer result) obtains the connection of well group during response results can be adopted from note
Relation, obtains the water injection rate data in Production development data and water injection time section from creation data table.
Then, using gauss of distribution function, to connect the production liquid maximum fluctuation value of each producing well in well group as constraint
Condition, one group of random number is produced for each percent continuity.Specifically, using gauss of distribution function to each
Individual percent continuity βijAll randomly generate one group of random number, when producing random number with the production liquid of each producing well most
Great fluctuation process value is used as constraint so that value will not be too general.
Then the valuation that formula (1) calculates water injection rate is substituted into respectively.Specifically, being asked for using multiple linear regression
The difference of water filling estimate and actual value.If for example, percent continuity is n, each producing m random value,
M can so be calculatednIndividual water injection rate valuation.In this process, it is necessary to carry out MxN wheel calculating.
M is contrasted respectivelynIndividual water injection rate valuation and the size of actual water injection rate, are obtained closest with actual water injection rate
Water injection rate, and its corresponding one group of random number.If the valuation of water injection rate and the difference convergence of actual value are (small
In certain precision), this group of random number is final percent continuity.Otherwise immediate water injection rate valuation is selected,
Using its corresponding one group of random number as the axle value of the Gaussian Profile of corresponding well, corresponding well percent continuity is produced
One group of new random number, water injection rate valuation is recalculated, until water injection rate valuation is equal with actual value.
Finally, due in well group between different injection-production wells percent continuity and should be 1, it should by connection derived above
Coefficient is normalized.Specifically, obtaining final percent continuity after being handled based on following formula
Wherein, n is the quantity of well, and m is the quantity of oil well in i-th mouthful of water injection well well group.
Followed by step S130 evaluation criterion establishment steps, it is distributed according to the percent continuity and statistical sample
Scope, obtains the evaluation criterion of connecting degree between connected well.Specifically, counting different injection-production well groups, different notes
The frequency that the percent continuity in period occurs is adopted, and draws cumulative frequency distributing graph, as shown in Figure 3.Specifically,
The different injection-production well groups that are calculated by step S120, different notes adopt the percent continuity in period and are used as sample, count
The frequency that different percent continuities occur, draws cumulative frequency distribution curve.Then in the cumulative frequency distributing graph
The percent continuity of the different accumulative perceptions of correspondence is read respectively, according to the corresponding percent continuity of different accumulative perceptions
Obtain the evaluation criterion of interwell communication degree.For example, being counted by substantial amounts of sample point, the distribution of sample point
Scope quadrisection, it is 25%, 50%, 75% to read accumulative perception respectively in cumulative frequency distribution curve
Percent continuity size, numerical value is respectively a, b, c, and now weak connectedness is percent continuity<A, medium company on the weak side
Lead to for percent continuity (a~b), medium inclined strong continune is percent continuity (b~c), and strong continune is percent continuity
>C, that is, obtain the evaluation criterion of connecting degree.
It is finally that step S140 interwell communication degree determines step, well is determined according to percent continuity and evaluation criterion
Between connecting degree.Specifically, interwell communication degree can be judged according to the value of percent continuity and evaluation criterion.
Come to carry out checking explanation to the present invention below by way of a specific example.Experiment well group is certain oil field S1
The injection-production well group G1 of fracture-cavity units one.Oil reservoir is the special oil reservoir based on solution cavity, fracture hole at this.Hole, hole,
Seam is in different ways and scale will constitute a variety of Storage categories, and such oilfield reservoir has seriously non-
Matter, static state is connected and the disconnected situation of dynamic between well often occurs, accordingly, it would be desirable to using this patent
Method calculates that interwell communication degree is strong and weak, understands fully that water injection well splits dispensing water to every mouthful of well producing well, after guidance
Phase waterflooding extraction is adjusted.
That tests comprises the following steps that:
(1) read in water injection well INJ water injection time, water injection rate, well group and produce from Production development tables of data
Dynamic creation data (production liquid, fluctuation severe degree etc. of well.
(2) the connection situation between G1 well group injection-production wells is obtained from tracer test report, it is known that water filling
Well INJ is connected with WELL1, WELL2, WELL3, WELL4, WELL5.Due to tracer test
Limitation, the connectedness of low water cut well cannot be obtained, response analysis producing well WELL6 is adopted using note for this
Pressure, production liquid, oil-producing there occurs obvious change during INJ water fillings, as shown in Figure 4.Moreover, WELL6
Working system does not change, therefore is also what is connected between INJ and WELL6 wells, so as to understand fully the company of whole well group
Clearance system.
(3) connecting degree between injection-production well is determined.Based on Monte Carlo method circulation M wheels, n times are often taken turns, are often taken turns
It is secondary to take one group of number at random, and respectively with the production liquid peak value of each producing well as constraint, using this group of nonce count
The estimate of water injection rate and the difference of actual value are calculated, after n times terminate, compares these differences and preserves wherein minimum
Value, carries out the calculating of next one, whether the minimum value for comparing two rounds restrains, if not restraining, by this group of number
Respectively as the average value of the Gaussian Profile of corresponding well, Gaussian Profile exploitation again is substituted into.M takes turns can
With the corresponding one group of random number of the convergency value and convergency value that obtain the minimum value.
(4) this group of random number is utilized into formula (3) zero dimension, obtains percent continuity.WELL1,
Between WELL2 ... ..., WELL6 and INJ percent continuity be respectively 0.07,0.13,0.22,0.25,0.15,
0.18。
Experiment show:
(1) TK663 well groups connecting degree is contrasted
According to the cumulative frequency distribution curve of all well groups of oil reservoir, as shown in figure 5, percent continuity is (right less than 0.1
Answer cumulative frequency 25%) be weak connectedness, 0.1~0.15 (correspondence cumulative frequency 50%) be it is medium on the weak side, 0.15~
0.23 (correspondence cumulative frequency 75%) is medium partially strong, is strong continune more than 0.23.Table 1 is based on Meng Teka
Lip river image factoring connects degree of strength contrast table with the G1 well groups that tracer test is obtained.Can from experimental result
To find out that result of calculation of the present invention and tracer result are basically identical.
Table 1TK663 well group result of calculation contrast tables
(2) water injection rate splits a point situation contrast
Percent continuity is applied in waterflooding extraction be exactly that a point coefficient is split in water filling, Fig. 6 for result of calculation b of the present invention and
Tracer test result a is contrasted, it can be seen that the present invention calculates obtained water filling and splits point coefficient and tracer test
As a result more it coincide, it is also basically identical that water injection rate splits point scale order.Simultaneously because tracer test can not
The connectedness of low wet producer is tested, unappropriated water injection rate is there is, and the present invention can be finally inversed by well6
Water filling split component.Therefore show that this method can comprehensively reflect the truth of fracture-pore reservoir, can be in reality
Middle application.
According to another aspect of the present invention, a kind of interwell communication degree is additionally provided to determine between system, including well
Connective determining module, percent continuity determining module, evaluation criterion set up module and interwell communication degree determines mould
Block.
Wherein, inter well connectivity determining module is used to determine the connectedness between well;Percent continuity determining module passes through
Founding mathematical models and specifying constraint, determine the percent continuity between connected well;Evaluation criterion sets up module,
According to percent continuity and statistical sample distribution, the evaluation criterion of connecting degree between connected well is obtained;Connect between well
Logical degree determining module, interwell communication degree is determined according to percent continuity and evaluation criterion.
The present invention is using taking a little at random near Gaussian Profile axle value so that axle value Step wise approximation actual value is tried to achieve
Percent continuity, compared with the general uniform thought taken a little, accelerates the velocity of approch to actual value;The opposing party
Face, according to the characteristics of actual production dynamic fluctuation, row constraint is entered to percent continuity, compared with least square method,
Be not in the situation that percent continuity is negative value so that result of calculation is more accurate.
While it is disclosed that embodiment as above, but described content is only to facilitate understand the present invention
And the embodiment used, it is not limited to the present invention.Technology people in any the technical field of the invention
Member, on the premise of spirit and scope disclosed in this invention are not departed from, can implementation formal and details
On make any modification and change, but the scope of patent protection of the present invention still must be with appended claims institute
The scope defined is defined.
Claims (10)
1. a kind of interwell communication degree determines method, including:
Inter well connectivity determines step;
Percent continuity determines step, by founding mathematical models and specifying constraint, determines the company between connected well
Logical coefficient;
Evaluation criterion establishment step, according to the percent continuity and statistical sample distribution, is obtained between connected well
The evaluation criterion of connecting degree;
Interwell communication degree determines step, and interwell communication journey is determined according to the percent continuity and the evaluation criterion
Degree.
2. determination method according to claim 1, it is characterised in that percent continuity determines that step enters one
Step includes:
Based on situation is adopted a note in injection-production well group, setting up the balanced flood mathematical modeling with constraints, institute more
Stating mathematical modeling includes the percent continuity;
Percent continuity is primarily determined that using monte carlo method;
The percent continuity primarily determined that is normalized to obtain final percent continuity.
3. determination method according to claim 2, it is characterised in that preliminary using monte carlo method
The step of determining percent continuity further comprises:
Reality in acquisition well group connected relation, well group in the production liquid maximum fluctuation value of each producing well and water injection time section
Border water injection rate data;
Using gauss of distribution function, to connect the production liquid maximum fluctuation value of each producing well in well group as constraints,
One group of random number is produced for each percent continuity;
The corresponding random number of each percent continuity is substituted into the mathematical modeling, water injection rate valuation is calculated;
According to water injection rate valuation and actual water injection rate data, obtain and the actual immediate water injection rate of water injection rate data
And corresponding random number;
Using this group of random number as the axle value of the Gaussian Profile of corresponding well, the new of corresponding well percent continuity is produced
One group of random number, return to water injection rate valuation step, until water injection rate valuation is equal with actual water injection rate, it is determined that
Percent continuity.
4. the determination method according to any one of claim 1-3, it is characterised in that the mathematical modeling
For:
<mrow>
<msub>
<mover>
<mi>i</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mn>0</mn>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>q</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
Wherein,Represent the water injection rate of i-th mouthful of water injection well;qj(n) producing well around water injection well is represented
Liquid output;β0jNon-equilibrium constant is represented, physical significance is the water influx of bottom water;βijRepresent jth mouthful production
The percent continuity of well and i-th mouthful of water injection well;J represents the producing well quantity in i-th mouthful of water filling well group,
Constraints is:
βij∈(wj-d,wj+d)
Wherein, wjThe zero dimension crest value of the creation data of jth mouthful producing well is represented, d represents that is permitted a change
The threshold value of change.
5. the determination method according to any one of claim 1-4, it is characterised in that evaluation criterion is set up
Step further comprises:
Count different injection-production well groups, different notes and adopt the frequency that the percent continuity in period occurs, and draw cumulative frequency
Distribution map;
Read the percent continuity of the different accumulative perceptions of correspondence respectively in the cumulative frequency distributing graph;
The evaluation criterion of interwell communication degree is obtained according to the corresponding percent continuity of different accumulative perceptions.
6. determination method according to claim 1, it is characterised in that inter well connectivity determines step bag
Include and inter well connectivity is judged by tracer test result.
7. determination method according to claim 1, it is characterised in that inter well connectivity determines step bag
Include and inter well connectivity is determined by oil well production dynamic data, including:
The Production development data of each producing well are obtained from creation data table;
According to the production curve situation of change of producing well during water filling, judge connective between injection-production well, wherein,
The interior pressure of producing well, production liquid occur in that obvious rising, or the moisture content of producing well go out during such as water filling
Show significant change, then illustrate connection between injection-production well.
8. determination method according to claim 7, it is characterised in that if producing well during water filling
While production curve changes, working system changes, then it is assumed that do not connected between injection-production well.
9. determination method according to claim 2, it is characterised in that based on following formula to originally determined
Percent continuity is normalized:
<mrow>
<msub>
<mover>
<mi>&beta;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<mo>...</mo>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>.....</mn>
<mo>,</mo>
<mi>n</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>......</mn>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
Wherein, n is well quantity, and m is the quantity of oil well in i-th mouthful of water injection well well group, βijRepresent jth mouthful
The percent continuity of producing well and i-th mouthful of water injection well.
10. a kind of interwell communication degree determines system, including,
Inter well connectivity determining module, for determining the connectedness between well;
Percent continuity determining module, by founding mathematical models and specifying constraint, determines the company between connected well
Logical coefficient;
Evaluation criterion sets up module, according to the percent continuity and statistical sample distribution, obtains between connected well
The evaluation criterion of connecting degree;
Interwell communication degree determining module, interwell communication journey is determined according to the percent continuity and the evaluation criterion
Degree.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610203017.6A CN107291667B (en) | 2016-04-01 | 2016-04-01 | Method and system for determining communication degree between wells |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610203017.6A CN107291667B (en) | 2016-04-01 | 2016-04-01 | Method and system for determining communication degree between wells |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107291667A true CN107291667A (en) | 2017-10-24 |
CN107291667B CN107291667B (en) | 2020-11-13 |
Family
ID=60087487
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610203017.6A Active CN107291667B (en) | 2016-04-01 | 2016-04-01 | Method and system for determining communication degree between wells |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107291667B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110259426A (en) * | 2019-07-02 | 2019-09-20 | 北京捷贝通石油技术股份有限公司 | Pressure alters the evaluation method of degree between a kind of unconventional platform well well |
CN110410065A (en) * | 2019-06-24 | 2019-11-05 | 中国石油化工股份有限公司 | The input parameters analysis method and device of multilayer inter-well dynamic connectivity model |
CN110965970A (en) * | 2018-09-29 | 2020-04-07 | 北京国双科技有限公司 | Method and device for determining correlation between water injection well and oil production well |
CN111027249A (en) * | 2019-12-10 | 2020-04-17 | 北京科技大学 | Machine learning-based inter-well connectivity evaluation method |
CN111749688A (en) * | 2020-08-10 | 2020-10-09 | 西南石油大学 | Method for predicting development position and direction of dominant seepage channel |
CN112360411A (en) * | 2020-11-30 | 2021-02-12 | 河海大学 | Local well pattern water injection development optimization method based on graph neural network |
CN113495295A (en) * | 2020-04-07 | 2021-10-12 | 中国石油天然气股份有限公司 | Method, device, equipment and storage medium for identifying weakly connected fracture-cavity reservoir |
CN114060003A (en) * | 2021-11-18 | 2022-02-18 | 中海石油(中国)有限公司海南分公司 | Characterization method for interwell connectivity of offshore complex fault block oil reservoir |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101725346A (en) * | 2009-12-15 | 2010-06-09 | 中国石油大学(华东) | Oil deposit inter-well dynamic connectivity inverting method |
CA2774181A1 (en) * | 2009-10-20 | 2011-04-28 | Exxonmobil Upstream Research Company | Method for quantitatively assessing connectivity for well pairs at varying frequencies |
CN103670369A (en) * | 2013-12-12 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and device for judging communication state between injection wells and production wells |
CN105089649A (en) * | 2014-04-17 | 2015-11-25 | 中国石油化工股份有限公司 | Method for establishing inter-well connectivity model |
CN105389467A (en) * | 2015-11-02 | 2016-03-09 | 中国地质大学(武汉) | Method and apparatus of acquiring inter-well communication relationship |
CN105404726A (en) * | 2015-11-02 | 2016-03-16 | 中国地质大学(武汉) | Method and device for reverting inter-well connectivity through Gaussian distribution-based capacitance model |
-
2016
- 2016-04-01 CN CN201610203017.6A patent/CN107291667B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2774181A1 (en) * | 2009-10-20 | 2011-04-28 | Exxonmobil Upstream Research Company | Method for quantitatively assessing connectivity for well pairs at varying frequencies |
CN101725346A (en) * | 2009-12-15 | 2010-06-09 | 中国石油大学(华东) | Oil deposit inter-well dynamic connectivity inverting method |
CN103670369A (en) * | 2013-12-12 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and device for judging communication state between injection wells and production wells |
CN105089649A (en) * | 2014-04-17 | 2015-11-25 | 中国石油化工股份有限公司 | Method for establishing inter-well connectivity model |
CN105389467A (en) * | 2015-11-02 | 2016-03-09 | 中国地质大学(武汉) | Method and apparatus of acquiring inter-well communication relationship |
CN105404726A (en) * | 2015-11-02 | 2016-03-16 | 中国地质大学(武汉) | Method and device for reverting inter-well connectivity through Gaussian distribution-based capacitance model |
Non-Patent Citations (3)
Title |
---|
关振良 等: "海上油井井间流动单元预测方法", 《海洋石油》 * |
苏群英: "井间连通性分析模型及其应用研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
邓英尔 等: "井间连通性的综合分析法", 《断块油气田》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110965970A (en) * | 2018-09-29 | 2020-04-07 | 北京国双科技有限公司 | Method and device for determining correlation between water injection well and oil production well |
CN110965970B (en) * | 2018-09-29 | 2022-02-11 | 北京国双科技有限公司 | Method and device for determining correlation between water injection well and oil production well |
CN110410065A (en) * | 2019-06-24 | 2019-11-05 | 中国石油化工股份有限公司 | The input parameters analysis method and device of multilayer inter-well dynamic connectivity model |
CN110410065B (en) * | 2019-06-24 | 2022-05-03 | 中国石油化工股份有限公司 | Input parameter analysis method and device for dynamic connectivity model among multiple wells |
CN110259426B (en) * | 2019-07-02 | 2021-12-07 | 捷贝通石油技术集团股份有限公司 | Method for evaluating pressure channeling degree between unconventional platform wells |
CN110259426A (en) * | 2019-07-02 | 2019-09-20 | 北京捷贝通石油技术股份有限公司 | Pressure alters the evaluation method of degree between a kind of unconventional platform well well |
CN111027249A (en) * | 2019-12-10 | 2020-04-17 | 北京科技大学 | Machine learning-based inter-well connectivity evaluation method |
CN111027249B (en) * | 2019-12-10 | 2021-02-26 | 北京科技大学 | Machine learning-based inter-well connectivity evaluation method |
CN113495295A (en) * | 2020-04-07 | 2021-10-12 | 中国石油天然气股份有限公司 | Method, device, equipment and storage medium for identifying weakly connected fracture-cavity reservoir |
CN113495295B (en) * | 2020-04-07 | 2024-01-30 | 中国石油天然气股份有限公司 | Weak-connectivity fracture-cavity reservoir identification method, device, equipment and storage medium |
CN111749688B (en) * | 2020-08-10 | 2022-03-18 | 西南石油大学 | Method for predicting development position and direction of dominant seepage channel |
CN111749688A (en) * | 2020-08-10 | 2020-10-09 | 西南石油大学 | Method for predicting development position and direction of dominant seepage channel |
CN112360411A (en) * | 2020-11-30 | 2021-02-12 | 河海大学 | Local well pattern water injection development optimization method based on graph neural network |
CN114060003A (en) * | 2021-11-18 | 2022-02-18 | 中海石油(中国)有限公司海南分公司 | Characterization method for interwell connectivity of offshore complex fault block oil reservoir |
CN114060003B (en) * | 2021-11-18 | 2024-02-23 | 中海石油(中国)有限公司海南分公司 | Characterization method for inter-well connectivity of offshore complex fault block oil reservoir |
Also Published As
Publication number | Publication date |
---|---|
CN107291667B (en) | 2020-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107291667A (en) | A kind of interwell communication degree determines method and system | |
CN105117511B (en) | A kind of characterizing method of fracture hole oil reservoir interwell communication passage and flow parameter | |
CN110952978B (en) | Drilling leakage fracture width prediction method based on neural network data mining | |
CN104879103B (en) | A kind of seperated layer water injection effect analysis method | |
CN104895550A (en) | Tight gas fracturing horizontal well numerical value well testing model building and solving method | |
CN104278989B (en) | Method for obtaining saturability index of low porosity and low permeability reservoir | |
CN104899411B (en) | A kind of reservoir productivity prediction model method for building up and system | |
CN104060985B (en) | A kind of stratified deposit profile control and water plugging blocking agent penetration depth method of testing and system | |
CN103912248A (en) | Method for predicting water contents of water-drive oilfields | |
CN104091069A (en) | Method for determining oil driving efficiency and sweep coefficient of all layers and positions of heterogeneous reservoir stratum | |
CN106150477A (en) | A kind of method determining single well controlled reserves | |
CN104915530B (en) | The method for building up of oil reservoir interwell communication relation | |
CN110006568A (en) | It is a kind of to obtain the method for three-dimensional ground stress using core and its obtain system | |
CN104612658A (en) | Well test analysis control system and method used for horizontal well subsection liquid producing identification | |
CN106407503A (en) | Reservoir fracture prediction method and apparatus | |
CN107895092A (en) | A kind of interwell communication quantitative evaluation method that modeling is adopted based on complex nonlinear note | |
CN104632150B (en) | The method that the offshore oilfield difference reasonable Liquid output of well group determines | |
CN106481315A (en) | Land sandstone oil reservoir individual well recoverable reserves quickly determines model and method for building up | |
CN110390154A (en) | A method of improving Complex reservoir reservoir numerical simulation efficiency | |
CN116992694B (en) | Tunnel blasting quality evaluation and optimization method | |
CN112160734B (en) | Injection and production well correlation analysis method and device, storage medium and computer equipment | |
Du et al. | A deep learning framework using graph convolutional networks for adaptive correction of interwell connectivity and gated recurrent unit for performance prediction | |
CN105260493A (en) | Semantics based oil well indicator diagram measuring method | |
CN103334740A (en) | Method for confirming oil drainage front edge with consideration of starting pressure gradient | |
CN110159260A (en) | For the main inflow direction method of discrimination of rift portion closure pressure break straight well and device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |