CN109403962A - Oil reservoir block Monitoring Indexes association analysis method - Google Patents
Oil reservoir block Monitoring Indexes association analysis method Download PDFInfo
- Publication number
- CN109403962A CN109403962A CN201811197578.5A CN201811197578A CN109403962A CN 109403962 A CN109403962 A CN 109403962A CN 201811197578 A CN201811197578 A CN 201811197578A CN 109403962 A CN109403962 A CN 109403962A
- Authority
- CN
- China
- Prior art keywords
- block
- association
- degree
- monitoring index
- oil
- 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.)
- Pending
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
Abstract
The invention discloses oil reservoir block Monitoring Indexes association analysis methods, the following steps are included: (a) extracts the dynamic test data of each individual well in block, pass through the dynamic test data of each individual well, calculation block monitoring index, the block monitoring index include: block water absorbing capacity, block oil productive capacity, block liquid-producing capacity, block percolation ability, block pressure holding level；(b) degree of association average value of each block monitoring index with respect to remaining block monitoring index is calculated, the degree of association between all block monitoring indexes is obtained；(c) according to calculation of relationship degree as a result, carrying out correlation analysis.The present invention is to solve the problems, such as that dynamic data can not make that overall contribution, each development evaluation index be mutually isolated, evaluation method has lag, passive in the prior art for oil reservoir block, it realizes and uses dynamic data into block global analysis, the connection between each evaluation index is established, initiative, the purpose of foresight of oil field development evaluation are improved.
Description
Technical field
The present invention relates to oil and gas development fields, and in particular to oil reservoir block Monitoring Indexes association analysis method.
Background technique
China's Dynamic Monitoring have passed through development in more than 40 years in oil and gas development field, and Dynamic Monitoring is on basis
Fruitful work has been carried out in experiment, theoretical research, log data acquisition, data processing explanation and application etc., and just
Step forms the distinctive mature Dynamic Monitoring of a set of tool, mainly tests including injection profile measuring technology, production profile
Technology, transient well test technology, pressure monitoring technology (stream pressure, static pressure test, dynamic and static liquid level test, oil pressure, casing pressure test),
Dynamic logging technology etc. between remaining oil saturation logging technique, engineering log technology and well.
The Eleventh Five-Year Plan period, with the complication of oil field at home development object, such as " three-lows " oil field, the big rule of deep natural gas
Mould exploitation, the diversification of development scheme, such as the scale of tertiary oil recovery technology, horizontal wells technology and heavy crude heat extraction technology
Using oil field development new situations propose stern challenge to Dynamic Monitoring.Currently, most domestic old filed has entered
The High water cut depletion stage of waterflooding extraction, the unfavorable factors such as oil reservoir is flat, section is particularly thorny, rate of water cut increase is fast restrict oil
Field production.
However in the prior art, oil field researcher is primarily focused on individual well data acquisition, the processing of information；Dynamically
Manager, and will be at monitoring by the explanation of treated monitoring index is widely used in individual well and single well measure adjustment etc.
Result is managed using the less of entire oil reservoir, yet there are no report as the research that entire block monitoring index is characterized is applied to
Road.
Analyzing an Exploitation Level of Oilfield and evaluating its development effectiveness needs many index synthesis to be measured, for example contains
Water climbing, strata pressure keep situation, producing degree of reservoir etc..But at present in oil reservoir development assessment of levels, analysis
When these indexs, often isolated each self-evaluation of carry out, does not set up a kind of connection between them, foothold is also
Overall merit is carried out to the development level that oil field is currently in；Dynamic analysis personnel also usually exist when carrying out dynamic analysis
After going wrong, starts to look for problem, proposes measure, adjustment improvement, even if some wells take measure, it is also difficult to be restored to pervious
Yield.The method of the above managemet of the developmet of oilfeild can not all accomplish that an oil field is allowed continually and steadily to keep developing, therefore, at present
Managemet of the developmet of oilfeild personnel are badly in need of a kind of a kind of system that is more proactive, can prejudging development situation, accomplish to oil field development
Manage more initiative, more proactive method.
Summary of the invention
The purpose of the present invention is to provide oil reservoir block Monitoring Indexes association analysis methods, to solve in the prior art
Dynamic data can not make that overall contribution, each development evaluation index be mutually isolated, evaluation method has lag, quilt for oil reservoir block
Dynamic problem is realized and uses dynamic data into block global analysis, and the connection between each evaluation index is established, and improves oil field
Initiative, the purpose of foresight of development evaluation.
The present invention is achieved through the following technical solutions:
Oil reservoir block Monitoring Indexes association analysis method, comprising the following steps:
(a) the dynamic test data for extracting each individual well in block passes through the dynamic test data of each individual well, calculation block prison
Index is surveyed, the block monitoring index includes: block water absorbing capacity, block oil productive capacity, block liquid-producing capacity, block seepage flow energy
Power, block pressure keep horizontal；
(b) degree of association average value of each block monitoring index with respect to remaining block monitoring index is calculated, is owned
The degree of association between block monitoring index；
(c) according to calculation of relationship degree as a result, carrying out correlation analysis.
It is mutually lonely overall contribution, each development evaluation index can not to be made for oil reservoir block for dynamic data in the prior art
Vertical, evaluation method has lag, passive problem, and the present invention proposes oil reservoir block Monitoring Indexes association analysis method,
Middle step (a) is the size for influencing contribution on block by individual well monitoring index, obtains the quantificational expression of block monitoring index.Step
Suddenly (b) successively determine a block monitoring index as dependent variable factor, remaining block monitoring index as independent variable factor,
Degree of being associated calculates and seeks its average value, after completing the calculating using each block monitoring index as dependent variable factor,
Carry out the calculation of relationship degree between all block monitoring indexes.The Computing Principle of the degree of association is as follows: setting each index factor as factor
The correlativity of each factor is analyzed and studied to curve by extracting a sample value of the moment curve on time series point, this
It is the multifactor Statistical Analysis Problems belonged in INFORMATION OF INCOMPLETE system.To seek the reality rule that this few data series are characterized
Rule, inventors herein proposes with the degree of association power of relationship between describing these index factors, size and the thought for being associated with order.Such as
Fruit sample data sequence reflects that the situation (including direction, size, speed etc.) changed between two factors is almost the same, then they it
Between correlation degree with regard to relatively high；Conversely, correlation degree is just smaller.That is by comparing because of the similitude one between prime sequences
Cause degree studies the correlation between factor.If two because the time graph shape of prime sequences is close, the pass between two indexes
Connection degree is with regard to larger；Conversely, the degree of association between two indexes is just smaller.(c) is entered step later, is monitored according to obtained all blocks
The degree of association between index carries out correlation analysis, to each block monitoring index be associated, so that each index no longer phase
It is mutually isolated, adequately the dynamic monitoring data of each individual well can be applied in the analysis of oil reservoir or even entire block, finally
Realize the initiative for improving oil field development evaluation, the purpose of foresight.
The calculation method of block water absorbing capacity are as follows:
Wherein, n1 --- the number of the surveyed water injection well of certain month block；M1 --- certain moon test water filling individual well explains the number of plies；
Certain month block water absorbing capacity is calculated using above formula.The block water absorbing capacity of calculating indicates entire block water absorbing capacity
The size of power, value represent the quality of block stratum water swelling elastomer.Variation by analyzing block some months water absorbing capacity becomes
Gesture may determine that whether block stratum water swelling elastomer changes.If trend stability, illustrate that block water swelling elastomer is good；If becoming
Gesture rises or falls, then illustrates that block water swelling elastomer has large change, it may be possible to have hypertonic water absorption layer, water absorption layer water flow to advance by leaps and bounds
Reason, it is also possible to which the water swelling of formation rock minerals, stratum impurity plug the reasons such as water flow duct and results in permeability
Decline.
The calculation method of block oil productive capacity are as follows:
Wherein, n2 --- the number of the surveyed producing well of certain month block；M2 --- certain moon tests individual well and explains the oil-producing number of plies；
Be firstly, it is all to block test individual well test data analyze, account for individual well by test layer oil-producing thickness
The ratio of total oil-producing thickness, to the contribution rate of individual well oil-producing, counterpart contribution rate will be multiplied by with respect to oil production again as the opposite oil production of layer
Superposition characterization individual well oil-producing coefficient；Then, using the ratio of the total oil production of oil output per well occupied area block testing well as individual well oil-producing pair
Individual well oil-producing coefficient is multiplied by corresponding contribution rate and is superimposed characterization block oil-producing coefficient again by the contribution rate of block oil-producing；Finally, will
Block oil-producing coefficient obtains block oil productive capacity multiplied by total oil-producing intensity of block.Certain month block oil-producing is calculated using above method
Ability.The oil productive capacity acquired indicates the size of entire block oil productive capacity, represents block oil-producing situation with the size of its value
Quality.Variation tendency by analyzing block some months oil productive capacity may determine that whether block stratum oil-producing situation becomes
Change.If trend stability, illustrate block oil-producing in order；If trend declines, illustrate that block oil-producing situation has large change,
May be that strata pressure supply is insufficient or formation rock mineral are in long-term water enchroachment (invasion) environment water swelling etc. leads to formation porosity
Permeability；Since a possibility that oil field is in the later development stage, and oil-producing condition trend obviously rises, is little.
The calculation method of block liquid-producing capacity are as follows:
Wherein, n3 --- the number of the surveyed producing well of certain month block；M3 --- certain moon tests producing well and explains fluid producing layer number.
Certain month block liquid-producing capacity is calculated using formula above.The liquid-producing capacity acquired indicates entire block liquid-producing capacity
Size, with the size of its value represent block producing well produce liquid condition quality, analyze block some months liquid-producing capacity variation
Trend may determine that block stratum produces whether liquid condition changes.If trend stability, illustrate that block produces liquid in order；If
Trend decline, then illustrate that block stratum condition has large change, it may be possible to which strata pressure keeps horizontal poor or formation rock stone ore
Substance leads to the reduction of in-place permeability in reasons such as long-term water enchroachment (invasion) environment water swellings；If in liquid-producing capacity variation tendency
It rises, it may be possible to cause due to production equipment increases lifting rate or so that strata pressure is maintained higher level etc. water filling.
The calculation method of block percolation ability are as follows:
Wherein, n4, m4 --- the number of transient well test in certain month；K --- in-place permeability (μm2)；
H --- core intersection (m)；μ --- fluid viscosity (mPas)；T --- the closed-in time (h)；
φ --- formation porosity (decimal)；Ct--- stratum system compressibility (MPa-1)；
--- the investigation area or drainage area (m of testing well2)；
The maximum magnitude that Well Test Technology is tested is namely by several mouthfuls of wells or the well group of more than ten mouthfuls of well constructions this small range
Stratum filtration situation, for that cannot reflect if understanding the variation of bigger region (such as certain fault block) stratum filtration situation.On oil field
The percolation ability for wouling have to consider entire block when in order to the problems such as analyzing fault block, boundary, with block institute testing well data come
Characterization block percolation ability is very important.Flow coefficient, effective permeability can be obtained by transient well test explanation.Other ginsengs
The number such as testing time is the closed-in time recorded when well testing, and porosity, viscosity, system compressibility etc. can pass through existing skill
Art is acquired as known conditions.Firstly, using the effective permeability of well test analysis, along with test effective thickness and formation fluid
Viscosity can indicate fluid fluid ability in the earth formation in institute's testing well influence area in the hope of formation flow coefficient；Then, with
The ratio between total testing well investigation area of the investigation area Zhan of testing well indicates the fluid ability of testing well to total testing well energy of flow
The contribution rate of power；Finally, with indicating that the flow coefficient of formation flow capacity of water is multiplied by corresponding contribution rate, then divided by surveyed tribute
The summation of rate is offered to characterize block percolation ability.Certain month block percolation ability is calculated using formula above.The block seepage flow acquired
Ability is to represent the size of entire block stratum filtration ability, represents the good of block producing well seepage flow situation with the size of its value
It is bad, it may determine that whether block formation rock and fluid flow level change with this.
Block pressure keeps horizontal calculation method are as follows:
Wherein, n5 --- certain moon test water filling well stream kill-job number；M5 --- certain moon tests water injection well stratum static pressure well number；
T1 --- certain moon tests producing well stratum static pressure well number；S --- certain moon test production well stream kill-job number.
Block pressure holding level can indicate block water filling supply to the supplement degree of block production pressure loss.If
Block pressure keeps horizontal higher, then illustrates the proper of pressure supply, and this production method is relatively reasonable；If block
Pressure keeps horizontal lower, then illustrates that pressure supply is insufficient, it is necessary to increase the intensity that various measures improve block supply.It utilizes
It is horizontal to calculate block pressure holding in certain month for formula above.The block pressure holding level acquired represents entire block strata pressure
The size of recharge degree illustrates that entire block note adopts pressure if several month block pressure keep horizontal size to tend towards stability
Level reaches best equilibrium state；Illustrate the stratum of water filling supplement if several month block pressure holding levels tend to rise
Energy has been more than the landing degree for producing well pressure；If some months, block pressure holding level tended to decline, illustrate water filling not
The landing that stratum energy can reasonably be supplemented has to find out that the water why injected cannot sufficiently supplement stratum energy in this case
The landing of amount, to direct the implementation of measure.
Further, the block monitoring index further include: employ coefficient, block potential advantages in block potential advantages layer position
Abundance；
Employ the calculation formula of coefficient in the block potential advantages layer position are as follows:
Wherein, d is that coefficient is employed in potential advantages layer position；M is toward annual data: the data of recent five years is taken when data are abundant, number
According to inadequate 5 years, add up data in recent years；N is the small number of plies of block；
Oil reservoir long-term injecting water, Oil Layer Reservoir physical property, pressure system etc. have all changed, and layer position is employed will also become therewith
Change, and in change procedure, substratum development degree is that possible reinforce, it is also possible to weaken, this be reservoir waterflooding process whether
A reasonable embodiment.The present invention refers to block potential advantages layer position and employs coefficient d, which can be to the conjunction of oil reservoir development
Rationality makees the characterization of an amount, can also make judgement appropriate to the potentiality of oil reservoir unreasonable in exploitation history, make oil reservoir pipe
Reason person is so quick that understand Reservoir behavior.Specifically, according to flow profile test data, since oil reservoir investment exploitation,
With 5 Nian Weiyi periods, statistics this period each substratum producing status.It is employed with the small position layer by layer most with the layer position of recent five years
Big use rate comparison (inadequate 5 years, accumulate data in recent years), then the value of each substratum is summed, as block
Employ coefficient in potential advantages layer position.D is bigger, and it is stronger to indicate that ability is employed in the block potential advantages layer position, illustrates block High water cut
Later period Exploitation Potential is larger, the characterization parameter can also be used for substratum, implements excavating resistance to substratum.
The calculation formula of the block potential advantages abundance are as follows:
Wherein, ΩoFor potential advantages abundance, 104t/km2；D is that coefficient is employed in potential advantages layer position；H is that reservoir is effectively thick
Degree, m；For porosity；ρoThe oil density for ground, t/m3；BoFor oil volume factor, m3/m3。
In high water-cut stage, the method for conventional characterization residue oilreserves is to calculate remaining oil reserves abundance or remaining oil to adopt storage
Measure abundance.This method reflects the remaining oil enrichment amount in block plane to a certain extent, but has ignored remaining oil
Fluid ability.The present invention uses potential advantages abundance Ωo, this characterization parameter gets on to reflect the surplus of underground from another meaning
Excess oil reserves, while the separation capacity of remaining oil is also reflected, so that it is determined that potentiality of remaining oil advantage field.With production development
It carries out, underground crude oil is unevenly distributed, and application advantage potentiality abundance is just that taping the latent power for high water cut viscous heavy crude oil remaining oil provides
Direction definitely.
Further, in step (b) degree of association calculation method are as follows:
(A) determine a block monitoring index as dependent variable factor, remaining block monitoring index as independent variable because
Element, if dependent variable data time sequence sequence is reference sequences x '0, each argument data time series, which is constituted, compares sequence x'0(i=1,
2,3 ..., n), this n+1 Variables Sequence is formed into matrix；
(B) nondimensionalization is carried out to Variables Sequence, obtained after nondimensionalization, each block monitoring index factor timing sequence
Arrange the matrix formed；
(C) the reference sequences absolute difference of sequence on the timing corresponding phase compared with each is calculated, absolute difference matrix is formed,
Extract maximum number Δ (max), the minimum number Δ (min) in absolute difference matrix；
(D) such as down conversion is made to the data in absolute difference matrix:
Obtain incidence coefficient matrix:
Wherein, ρ is resolution ratio, 0 < ρ < 1；N is the timing point length of dependent variable sequence
(E) calculating correlation:
(F) it sorts according to the size of the degree of association, determines the block monitoring index as dependent variable factor, become with as oneself
Between remaining block monitoring index of amount factor relationship power be associated with order.
The principle of the calculation method of the above-mentioned degree of association are as follows: set each index factor as moment curve, by extracting a factor song
The correlativity of each factor is analyzed and studied to sample value of the line on time series point, this is belonged in INFORMATION OF INCOMPLETE system
Multifactor Statistical Analysis Problems.To seek the reality rule that this few data series are characterized, the invention proposes with association
It spends to describe the power of relationship, size and the thought for being associated with order between these index factors.If sample data sequence reflects
The situation (including direction, size, speed etc.) changed between two factors is almost the same, then the correlation degree between them just compares
It is high；Conversely, correlation degree is just smaller.That is come between research factor by comparing because of the similitude consistent degree between prime sequences
Correlation.If two because the time graph shape of prime sequences is close, the degree of association between two indexes is with regard to larger；Conversely, two
The degree of association between index is just smaller.This programme can finally obtain the block monitoring index as dependent variable factor, with as from
Between remaining block monitoring index of Variable Factors relationship power and be associated with order, thus between analysis indexes curve correlation pass
System provides quantitative scientific basis.
Preferably, 0.1≤ρ≤0.5.Resolution ratio ρ takes lesser value, and the significant of difference between incidence coefficient can be improved
Property.
Preferably, the nondimensionalization that Variables Sequence carries out is carried out by equalization method or first value method.Equalization method and
First value method is the nondimensionalization calculation method received and had.
Wherein, the calculation formula of equalization method are as follows:
The calculation formula of first value method are as follows:
The matrix that each block monitoring index factor time series is formed after nondimensionalization are as follows:
The reference sequences absolute difference of sequence on the timing corresponding phase compared with each is calculated, following absolute difference square is formed
Battle array:
Wherein Δ0i(k)=| x0(k)-xi(k)
I=1,2 ..., n；K=1,2 ..., N
It is maximum difference and lowest difference to obtain maximum number and minimum number in absolute difference battle array:
Preferably, correlation analysis described in step (c) is correlation classification and/or correlation quantitative analysis.
Preferably, correlation classification the following steps are included:
(I) degree of association between all variables is constituted into matrix R=(rij)N×N, wherein 0 < rij<1；The variable includes
Dependent variable and independent variable；
(II) transitive closure matrix of matrix R is acquired using quadratic method；
(III) it is based on transitive closure matrix, makes dynamic clustering figure according to common equivalence relation.According to dendrogram
The classification of correlation required for excavating.
Preferably, the correlation quantitative analysis the following steps are included:
(1) degree of association size between all block monitoring indexes is calculated, is constructed according to calculated all degrees of association
Fuzzy Correlation matrix；
(2) fuzzy equivalent matrix is calculated according to fuzzy clustering principle；Each equivalence class is found out by correlation classification method；
(3) dynamic screening is carried out to the equivalence class that step (2) is found out, it is obviously partially strong removes the block monitoring index degree of association
Equivalence class obtains a new Fuzzy Correlation matrix；
(4) step (2)~(3) are repeated, the block monitoring index obtained in step (3) individually at class when, terminate dynamic
State screening；
(5) degree of association that block monitoring index class obtained in step (4) is mutual constitutes fuzzy similarity matrix, simultaneously
The fuzzy equivalent matrix between block monitoring index is calculated by fuzzy clustering principle；Block is found out by correlation classified calculating method
The equivalence class of monitoring index；
(6) degree of association in each block monitoring index equivalence class between each block monitoring index is calculated；It is supervised by each block
The size for surveying the degree of association in index equivalence class determines the representative block monitoring index of each block monitoring index.So that it is determined that going out
The prediction index system of oil field block monitoring index.
Compared with prior art, the present invention having the following advantages and benefits:
Oil reservoir block Monitoring Indexes association analysis method of the present invention, according between obtained all block monitoring indexes
The degree of association, carry out correlation analysis, so that each block monitoring index be associated so that each index is no longer mutually isolated,
Adequately the dynamic monitoring data of each individual well can be applied in the analysis of oil reservoir or even entire block, final realize is improved
Initiative, the purpose of foresight of oil field development evaluation.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below with reference to embodiment, the present invention is made
Further to be described in detail, exemplary embodiment of the invention and its explanation for explaining only the invention, are not intended as to this
The restriction of invention.
Embodiment 1:
Oil reservoir block Monitoring Indexes association analysis method, comprising the following steps: each individual well is dynamic in (a) extraction block
State test data, by the dynamic test data of each individual well, calculation block monitoring index, the block monitoring index includes: area
Block water absorbing capacity, block oil productive capacity, block liquid-producing capacity, block percolation ability, block pressure keep horizontal；(b) it calculates every
One block monitoring index obtains between all block monitoring indexes with respect to the degree of association average value of remaining block monitoring index
The degree of association；(c) according to calculation of relationship degree as a result, carrying out correlation analysis.
The wherein calculation method of the block water absorbing capacity are as follows:
Wherein, n1 --- the number of the surveyed water injection well of certain month block；M1 --- certain moon test water filling individual well explains the number of plies；
The calculation method of the block oil productive capacity are as follows:
Wherein, n2 --- the number of the surveyed producing well of certain month block；M2 --- certain moon tests individual well and explains the oil-producing number of plies；
The calculation method of the block liquid-producing capacity are as follows:
Wherein, n3 --- the number of the surveyed producing well of certain month block；M3 --- certain moon tests producing well and explains fluid producing layer number.
The calculation method of the block percolation ability are as follows:
Wherein, n4, m4 --- the number of transient well test in certain month；K --- in-place permeability (μm2)；
H --- core intersection (m)；μ --- fluid viscosity (mPas)；T --- the closed-in time (h)；
φ --- formation porosity (decimal)；Ct--- stratum system compressibility (MPa-1)；
--- the investigation area or drainage area (m of testing well2)；
The block pressure keeps horizontal calculation method are as follows:
Wherein, n5 --- certain moon test water filling well stream kill-job number；M5 --- certain moon tests water injection well stratum static pressure well number；
T1 --- certain moon tests producing well stratum static pressure well number；S --- certain moon test production well stream kill-job number.
By taking the YS fault block of certain domestic oil field as an example, the dynamic test data of each individual well is as shown in the table in the present embodiment:
The block water absorbing capacity calculated result of the YS fault block month in the oil field is as follows:
The block oil productive capacity calculated result of oil field YS fault block month is as follows:
The block liquid-producing capacity calculated result of oil field YS fault block month is as follows:
The block percolation ability calculated result of oil field YS fault block month is as follows:
The calculated result that the block pressure of oil field YS fault block month keeps horizontal is as follows:
To sum up the summary table for obtaining block monitoring index is as follows:
The present embodiment is calculating block water absorbing capacity as shown above, block oil productive capacity, block liquid-producing capacity, area
After block percolation ability and block pressure keep horizontal, in order to which prediction block development situation must also screen the biggish block of the degree of association
Monitoring index is used for the prediction to block monitoring index.
Most concerned oil field development production is exactly block oil productive capacity and block percolation ability, and the present embodiment is with the two areas
Block monitoring index makees correlation analysis with other block monitoring index features as Main Analysis object, by it, verifies itself and its
The size of his the block monitoring index degree of association.
By taking block oil productive capacity as an example, 4 block monitoring index features of itself and other are done into correlation analysis, can similarly be divided
Block water absorbing capacity, block liquid-producing capacity, block percolation ability and block pressure holding level and remaining four block are not made
The size of monitoring index feature association degree.
First using block oil productive capacity as dependent variable factor, remaining block monitoring index is as independent variable factor.Meter
It calculates block oil productive capacity and block water absorbing capacity, block liquid-producing capacity, block percolation ability and block pressure keeps horizontal pass
Connection degree γ0i, block oil productive capacity and the calculated result for being associated with angle value of remaining four block monitoring index is as follows:
The degree of association of block oil productive capacity and block water absorbing capacity: γBlock water absorbing capacity=0.5657；
The degree of association of block oil productive capacity and block liquid-producing capacity: γBlock liquid-producing capacity=0.5952；
The degree of association of block oil productive capacity and block percolation ability: γBlock percolation ability=0.3534；
Block oil productive capacity and block pressure keep the horizontal degree of association: γBlock pressure keeps horizontal=0.5792.
According to correlation analysis theory, degree of association very little between index factor or and absolute difference it is very big, then illustrate index
Correlation degree between factor is weaker, regards as and incidence relation is not present between two indices, illustrate that two indices are uncorrelated.By above
Block oil productive capacity and block water absorbing capacity, block liquid-producing capacity, block percolation ability and the holding of block pressure acquired is horizontal
The degree of association can be seen that block oil productive capacity and block water absorbing capacity, block liquid-producing capacity, block percolation ability, block pressure
The degree of association between water holding of trying hard to keep is flat is smaller, therefore this four indexs and block oil productive capacity correlation be not strong, it should take this
Five index prediction block exploitations.
Above-described specific embodiment has carried out further the purpose of the present invention, technical scheme and beneficial effects
It is described in detail, it should be understood that being not intended to limit the present invention the foregoing is merely a specific embodiment of the invention
Protection scope, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should all include
Within protection scope of the present invention.
Claims (9)
1. oil reservoir block Monitoring Indexes association analysis method, which comprises the following steps:
(a) the dynamic test data for extracting each individual well in block, by the dynamic test data of each individual well, calculation block monitoring refers to
Mark, the block monitoring index include: block water absorbing capacity, block oil productive capacity, block liquid-producing capacity, block percolation ability,
Block pressure keeps horizontal；
(b) degree of association average value of each block monitoring index with respect to remaining block monitoring index is calculated, all blocks are obtained
The degree of association between monitoring index；
(c) according to calculation of relationship degree as a result, carrying out correlation analysis.
2. oil reservoir block Monitoring Indexes association analysis method according to claim 1, which is characterized in that
The calculation method of the block water absorbing capacity are as follows:
Wherein, n1 --- the number of the surveyed water injection well of certain month block；M1 --- certain moon test water filling individual well explains the number of plies；
The calculation method of the block oil productive capacity are as follows:
Wherein, n2 --- the number of the surveyed producing well of certain month block；M2 --- certain moon tests individual well and explains the oil-producing number of plies；
The calculation method of the block liquid-producing capacity are as follows:
Wherein, n3--- the number of the surveyed producing well of certain month block；m3--- certain moon tests producing well and explains fluid producing layer number.
The calculation method of the block percolation ability are as follows:
Wherein, n4, m4 --- the number of transient well test in certain month；K --- in-place permeability (μm2)；
H --- core intersection (m)；μ --- fluid viscosity (mPas)；T --- the closed-in time (h)；
φ --- formation porosity (decimal)；Ct--- stratum system compressibility (MPa-1)；
--- the investigation area or drainage area (m of testing well2)；
The block pressure keeps horizontal calculation method are as follows:
Wherein, n5 --- certain moon test water filling well stream kill-job number；M5 --- certain moon tests water injection well stratum static pressure well number；
T1 --- certain moon tests producing well stratum static pressure well number；S --- certain moon test production well stream kill-job number.
3. oil reservoir block Monitoring Indexes association analysis method according to claim 1, which is characterized in that the block
Monitoring index further include: employ coefficient, block potential advantages abundance in block potential advantages layer position；
Employ the calculation formula of coefficient in the block potential advantages layer position are as follows:
Wherein, d is that coefficient is employed in potential advantages layer position；M is toward annual data: the data of recent five years is taken when data are abundant, data are not
Enough 5 years, add up data in recent years；N is the small number of plies of block；
The calculation formula of the block potential advantages abundance are as follows:
Wherein, ΩoFor potential advantages abundance, 104t/km2；D is that coefficient is employed in potential advantages layer position；H is reservoir effective thickness, m；For porosity；ρoThe oil density for ground, t/m3；BoFor oil volume factor, m3/m3。
4. oil reservoir block Monitoring Indexes association analysis method according to claim 1, which is characterized in that step (b)
The calculation method of the middle degree of association are as follows:
(A) determine a block monitoring index as dependent variable factor, remaining block monitoring index as independent variable factor, if
Dependent variable data time sequence sequence is reference sequences x '0, each argument data time series, which is constituted, compares sequence x'0(i=1,2,
3 ..., n), this n+1 Variables Sequence is formed into matrix；
(B) nondimensionalization is carried out to Variables Sequence, obtained after nondimensionalization, each block monitoring index factor time series shape
At matrix；
(C) the reference sequences absolute difference of sequence on the timing corresponding phase compared with each is calculated, absolute difference matrix is formed, is extracted
Maximum number Δ (max), minimum number Δ (min) in absolute difference matrix；
(D) such as down conversion is made to the data in absolute difference matrix:
Obtain incidence coefficient matrix:
Wherein, ρ is resolution ratio, 0 < ρ < 1；N is the timing point length of dependent variable sequence
(E) calculating correlation:
(F) according to the degree of association size sort, determine the block monitoring index as dependent variable factor, with as independent variable because
Element remaining block monitoring index between relationship power be associated with order.
5. oil reservoir block Monitoring Indexes association analysis method according to claim 4, which is characterized in that 0.1≤ρ≤
0.5。
6. oil reservoir block Monitoring Indexes association analysis method according to claim 4, which is characterized in that variable sequence
The nondimensionalization that column carry out is carried out by equalization method or first value method.
7. oil reservoir block Monitoring Indexes association analysis method according to claim 1, which is characterized in that step (c)
Described in correlation analysis be correlation classify and/or correlation quantitative analysis.
8. oil reservoir block Monitoring Indexes association analysis method according to claim 7, which is characterized in that the correlation
Property classification the following steps are included:
(I) degree of association between all variables is constituted into matrix R=(rij)N×N, wherein 0 < rij<1；
(II) transitive closure matrix of matrix R is acquired using quadratic method；
(III) it is based on transitive closure matrix, makes dynamic clustering figure according to common equivalence relation.
9. oil reservoir block Monitoring Indexes association analysis method according to claim 7, which is characterized in that the correlation
Property quantitative analysis the following steps are included:
(1) degree of association size between all block monitoring indexes is calculated, is constructed according to calculated all degrees of association fuzzy
Incidence matrix；
(2) fuzzy equivalent matrix is calculated according to fuzzy clustering principle；Each equivalence class is found out by correlation classification method；
(3) equivalence class found out to step (2) carries out dynamic screening, removes the block monitoring index degree of association obviously partially strong equivalence
Class obtains a new Fuzzy Correlation matrix；
(4) step (2)~(3) are repeated, the block monitoring index obtained in step (3) individually at class when, terminate dynamic and sieve
Choosing；
(5) degree of association that block monitoring index class obtained in step (4) is mutual constitutes fuzzy similarity matrix and presses mould
Paste cluster principle calculates the fuzzy equivalent matrix between block monitoring index；Block monitoring is found out by correlation classified calculating method
The equivalence class of index；
(6) degree of association in each block monitoring index equivalence class between each block monitoring index is calculated；Refer to by the monitoring of each block
The size of the degree of association determines the representative block monitoring index of each block monitoring index in mark equivalence class.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811197578.5A CN109403962A (en) | 2018-10-15 | 2018-10-15 | Oil reservoir block Monitoring Indexes association analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811197578.5A CN109403962A (en) | 2018-10-15 | 2018-10-15 | Oil reservoir block Monitoring Indexes association analysis method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109403962A true CN109403962A (en) | 2019-03-01 |
Family
ID=65467973
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811197578.5A Pending CN109403962A (en) | 2018-10-15 | 2018-10-15 | Oil reservoir block Monitoring Indexes association analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109403962A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111827967A (en) * | 2020-06-09 | 2020-10-27 | 北京永源思科技发展有限公司 | Oil and gas reservoir analysis method |
| CN111986041A (en) * | 2019-05-22 | 2020-11-24 | 中国石油天然气股份有限公司 | Oil reservoir development method and device |
| CN113494290A (en) * | 2021-08-02 | 2021-10-12 | 中国石油化工股份有限公司 | Oil layer thickness measurement system based on big data |
| CN114086928A (en) * | 2020-08-07 | 2022-02-25 | 中国石油化工股份有限公司 | Multi-thin-layer thick oil effective reservoir classification evaluation method |
| CN116066067A (en) * | 2021-12-30 | 2023-05-05 | 中国石油天然气集团有限公司 | Evaluation method for potential of residual oil in oil field and application thereof |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101042048A (en) * | 2006-03-24 | 2007-09-26 | 中国石油天然气股份有限公司 | Complicated fault block fluvial facies reservoir oil water well using situation split system |
| CN101295374A (en) * | 2007-04-28 | 2008-10-29 | 中国石油天然气股份有限公司 | Dynamic quantitative method for oil-field development influenced by multifactor |
| WO2009032924A2 (en) * | 2007-09-04 | 2009-03-12 | Schlumberger Canada Limited | Method and system for increasing production of a reservoir using lateral wells |
| CN101719231A (en) * | 2008-10-09 | 2010-06-02 | 西南石油大学 | Dynamic differential simulation forecasting method and data processing handheld device applying same |
| CN105158804A (en) * | 2015-09-17 | 2015-12-16 | 刘立峰 | Solution-pore-type carbonate reservoir oil and gas enrichment degree obtaining method and device |
| CN105626036A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Oil reservoir engineering calculation method for determining reasonable oil reservoir liquid yield |
| CN106803010A (en) * | 2016-08-10 | 2017-06-06 | 中国石油大学(华东) | For the Fuzzy Grey comprehensive evaluation method and device of low permeability reservoir quantitative assessment |
| CN107038516A (en) * | 2016-10-20 | 2017-08-11 | 中国石油化工股份有限公司 | Complex Fault Block Oil Reservoir water flooding effectiveness quantitative evaluation method is oozed in one kind |
| CN108240208A (en) * | 2018-02-05 | 2018-07-03 | 东北石油大学 | A kind of oilfield water flooding classification well group development effectiveness is to marking method |
-
2018
- 2018-10-15 CN CN201811197578.5A patent/CN109403962A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101042048A (en) * | 2006-03-24 | 2007-09-26 | 中国石油天然气股份有限公司 | Complicated fault block fluvial facies reservoir oil water well using situation split system |
| CN101295374A (en) * | 2007-04-28 | 2008-10-29 | 中国石油天然气股份有限公司 | Dynamic quantitative method for oil-field development influenced by multifactor |
| WO2009032924A2 (en) * | 2007-09-04 | 2009-03-12 | Schlumberger Canada Limited | Method and system for increasing production of a reservoir using lateral wells |
| CN101719231A (en) * | 2008-10-09 | 2010-06-02 | 西南石油大学 | Dynamic differential simulation forecasting method and data processing handheld device applying same |
| CN105626036A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Oil reservoir engineering calculation method for determining reasonable oil reservoir liquid yield |
| CN105158804A (en) * | 2015-09-17 | 2015-12-16 | 刘立峰 | Solution-pore-type carbonate reservoir oil and gas enrichment degree obtaining method and device |
| CN106803010A (en) * | 2016-08-10 | 2017-06-06 | 中国石油大学(华东) | For the Fuzzy Grey comprehensive evaluation method and device of low permeability reservoir quantitative assessment |
| CN107038516A (en) * | 2016-10-20 | 2017-08-11 | 中国石油化工股份有限公司 | Complex Fault Block Oil Reservoir water flooding effectiveness quantitative evaluation method is oozed in one kind |
| CN108240208A (en) * | 2018-02-05 | 2018-07-03 | 东北石油大学 | A kind of oilfield water flooding classification well group development effectiveness is to marking method |
Non-Patent Citations (10)
| Title |
|---|
| WEIXIAO ET AL.: "Multiobjective linear programming model on injection oilfield recovery system", 《COMPUTERS & MATHEMATICS WITH APPLICATIONS》 * |
| 丁亚军等: "建立油田开发动态预测指标体系的新方法", 《西部探矿工程》 * |
| 何勇明等: "《低渗透油藏渗流机理及应用》", 30 September 2014, 地质出版社 * |
| 唐然等: "六中区克下组高含水油藏剩余油分布及挖潜研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 * |
| 姚军等: "《油藏工程原理与方法》", 31 May 2016, 中国石油大学出版社 * |
| 宋万超等: "《高含水期油田开发技术和方法》", 15 March 2003, 地质出版社 * |
| 张平等: "不同类型油藏动态测试指标表征方法应用研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
| 方凌云等: "《砂岩油藏注水开发动态分析》", 31 July 1998, 石油工业出版社 * |
| 赵智勇等: "动态监测指标在预判区块开发中的应用研究", 《西南石油大学学报(自然科学版)》 * |
| 钟仪华等: "特高含水油田开发规划动态预测方法研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111986041A (en) * | 2019-05-22 | 2020-11-24 | 中国石油天然气股份有限公司 | Oil reservoir development method and device |
| CN111827967A (en) * | 2020-06-09 | 2020-10-27 | 北京永源思科技发展有限公司 | Oil and gas reservoir analysis method |
| CN114086928A (en) * | 2020-08-07 | 2022-02-25 | 中国石油化工股份有限公司 | Multi-thin-layer thick oil effective reservoir classification evaluation method |
| CN114086928B (en) * | 2020-08-07 | 2023-11-17 | 中国石油化工股份有限公司 | Multi-thin-layer thick oil effective reservoir classification evaluation method |
| CN113494290A (en) * | 2021-08-02 | 2021-10-12 | 中国石油化工股份有限公司 | Oil layer thickness measurement system based on big data |
| CN113494290B (en) * | 2021-08-02 | 2023-11-14 | 中国石油化工股份有限公司 | Oil layer thickness measurement system based on big data |
| CN116066067A (en) * | 2021-12-30 | 2023-05-05 | 中国石油天然气集团有限公司 | Evaluation method for potential of residual oil in oil field and application thereof |
| CN116066067B (en) * | 2021-12-30 | 2024-02-02 | 中国石油天然气集团有限公司 | Evaluation method for potential of residual oil in oil field and application thereof |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109403962A (en) | Oil reservoir block Monitoring Indexes association analysis method | |
| Fan et al. | Understanding gas-condensate reservoirs | |
| CN110321595B (en) | Fault sealing evaluation method for extracting static quality coefficient by logging | |
| CN110288258A (en) | A kind of high water-cut reservoir Tapping Residual Oil method | |
| CN109594968A (en) | Fracture parameters evaluation method and system after a kind of shale gas multistage pressure break horizontal well pressure | |
| CN104747185A (en) | Heterogeneous reservoir stratum synthetic classifying evaluation method | |
| Bhattacharya et al. | Analysis of production history for unconventional gas reservoirs with statistical methods | |
| Chiotoroiu et al. | Forecasting incremental oil production of a polymer-pilot extension in the Matzen field including quantitative uncertainty assessment | |
| CN108343430A (en) | Fault sealing property fuzzy mathematics comprehensive quantitative evaluation method | |
| CN111878073B (en) | Method and device for evaluating fracturing effect of tight reservoir | |
| Masalmeh et al. | Developing high resolution static and dynamic models for waterflood history matching and EOR evaluation of a Middle Eastern carbonate reservoir | |
| CN105138749A (en) | Method for recovering original organic carbon of highly-mature hydrocarbon source rocks under geological constraint | |
| Craig et al. | Permeability, pore pressure, and leakoff-type distributions in Rocky Mountain basins | |
| Kharrat et al. | Rock Typeand Permeability Prediction of a Heterogeneous Carbonate Reservoir Using Artificial Neural Networks Based on Flow Zone Index Approach | |
| Konwar et al. | Understanding mauddud waterflood performance in a heterogeneous carbonate reservoir with surveillance data and ensemble of analytical tools | |
| Izgec | Understanding waterflood performance with modern analytical techniques | |
| Wang et al. | Drainage type classification and key controlling factors of productivity for CBM wells in the Zheng Zhuang Area, Southern Qinshui Basin, North China | |
| Fenik et al. | Criteria for ranking realizations in the investigation of SAGD reservoir performance | |
| Parra et al. | Application of the Producer-Based Capacitance Resistance Model to Undersaturated Oil Reservoirs in Primary Recovery | |
| Zhang et al. | Primary Controlling Factors of Coalbed Methane Well Productivity and High Productive Well Patterns in Eastern Yunnan and Western Guizhou, China | |
| Bhatti et al. | Generation of Saturation Functions from Historic SCAL Data for Simulation of a Super-Giant Carbonate Reservoir | |
| CN112241835A (en) | Deep shaft project water inrush disaster multi-source information evaluation method | |
| Varayesi et al. | The Effect of the Drainage Area on Well and Reservoir Performance | |
| CN106401575B (en) | Ultra-high water cut stage core hole Water Flooding Layer initial oil saturation calculation method | |
| Cramer | Analyzing Well Performance in Hydraulically Fractured Gas Wells: Non-Ideal Cases |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190301 |
|
| RJ01 | Rejection of invention patent application after publication |