CN107965318A - A kind of method of Volcanic Reservoir effective reservoir quantitative classification - Google Patents
A kind of method of Volcanic Reservoir effective reservoir quantitative classification Download PDFInfo
- Publication number
- CN107965318A CN107965318A CN201711275674.2A CN201711275674A CN107965318A CN 107965318 A CN107965318 A CN 107965318A CN 201711275674 A CN201711275674 A CN 201711275674A CN 107965318 A CN107965318 A CN 107965318A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- volcanic
- effective
- classification
- discriminant function
- 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
- 238000000034 method Methods 0.000 title claims abstract description 64
- 239000011435 rock Substances 0.000 claims abstract description 60
- 238000007621 cluster analysis Methods 0.000 claims abstract description 7
- 238000012886 linear function Methods 0.000 claims description 12
- 230000035699 permeability Effects 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 16
- 238000011156 evaluation Methods 0.000 abstract description 6
- 230000002349 favourable effect Effects 0.000 abstract description 3
- 230000018109 developmental process Effects 0.000 description 15
- 238000004519 manufacturing process Methods 0.000 description 10
- 208000035126 Facies Diseases 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 7
- 239000012530 fluid Substances 0.000 description 7
- 238000011160 research Methods 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 229910052770 Uranium Inorganic materials 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000004880 explosion Methods 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 239000004575 stone Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 239000000428 dust Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 2
- WKBPZYKAUNRMKP-UHFFFAOYSA-N 1-[2-(2,4-dichlorophenyl)pentyl]1,2,4-triazole Chemical compound C=1C=C(Cl)C=C(Cl)C=1C(CCC)CN1C=NC=N1 WKBPZYKAUNRMKP-UHFFFAOYSA-N 0.000 description 1
- 244000144725 Amygdalus communis Species 0.000 description 1
- 235000011437 Amygdalus communis Nutrition 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 208000010392 Bone Fractures Diseases 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 235000020224 almond Nutrition 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 244000144992 flock Species 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000011031 topaz Substances 0.000 description 1
- 229910052853 topaz Inorganic materials 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK 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
- E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
- E21B49/087—Well testing, e.g. testing for reservoir productivity or formation parameters
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The present invention provides a kind of method of Volcanic Reservoir effective reservoir quantitative classification.This method comprises the following steps:Obtain the reservoir quality parameter with reference to volcanic rock reservoir;Split clustering procedure using K means and cluster analysis is carried out to the reservoir quality parameter, to establish the criteria for classification of Volcanic Reservoir effective reservoir;Based on the criteria for classification of the Volcanic Reservoir effective reservoir, the linear discriminant function of all kinds of effective reservoirs is established using linear discriminant function method;Based on the linear discriminant function, the type of target volcanic rock reservoir is determined.Technical solution provided by the invention further for the deployment of such oil reservoir Efficient Development can provide technical support and be the favourable reference of perfect offer of Volcanic Reservoir effective reservoir evaluation of classification standard and technical system.
Description
Technical field
The present invention relates to a kind of method of Volcanic Reservoir effective reservoir quantitative classification, belongs to shale oil and hides development technique neck
Domain.
Background technology
In alive boundary stone exploration activity and Development history, clastic rock is often reservoir exploration and exploitation with carbonate reservoir
Main object, and relatively deep volcanic rock reservoir is buried, because its complicated characteristics of lithology and lithofacies, RESERVOIR RECOGNITION are pre- with distribution
Survey is often ignored there are larger difficulty by us.
Since early 20th century, with the development and the continuous improvement of exploration and development technology of petroleum industry, domestic and international volcano
Oil gas pool exploration and development constantly succeeds, such as the U.S., Cuba, Japan, state of Ghana, domestic huge port, the Liaohe River, Xinjiang,
The oil fields such as grand celebration, volcanic reservoirs also progressively become oil-gas exploration and development key areas, its reservoir characteristic becomes particular lithologic
The emphasis of reservoir study.
At present, domestic and foreign scholars are in Volcanic uranium deposit Lithofacies Identification, reservoir space and volcanic rock reservoir log interpretation model
Etc., great progress has been achieved, has been achieved in Volcanic uranium deposit, physical property and reservoir space etc. extensive common
Know, but in terms of Reservoir Classification, the object of forefathers' research is mostly to be directed to clastic rock oil reservoir, accidental effectively to be stored up on volcanic gas reservoir
The effective reservoir of system has been carried out Song-liao basin Xu man defensive wall volcanic gas reservoir in the research of layer classification, such as topaz imperial (2010)
Research.
In general, effective reservoir refer to possess preserve with seepage flow movable fluid (based on hydrocarbon fluid) ability, existing
There is the reservoir of commercial mining value under stage process technology and economic condition;And invalid reservoir then refers to that reservoir properties are poor,
There is no liquid-producing capacity or the reservoir less than dried layer standard.
At present, for the evaluation of classification of volcanic rock reservoir effective reservoir, ripe and effective method and skill are still lacked
Art, thus constrain the Efficient Development of such oil reservoir High-quality Reservoir.Therefore, volcanic rock reservoir effective reservoir evaluation of classification side is established
Method becomes one of this area technical problem urgently to be resolved hurrily.
The content of the invention
In order to solve the above technical problems, it is an object of the invention to provide a kind of Volcanic Reservoir effective reservoir quantitative classification
Method.This method can quickly identify volcanic rock reservoir effective reservoir, and provide technological guidance for volcanic rock reservoir Efficient Development
With reference.
To reach above-mentioned purpose, the present invention provides a kind of method of Volcanic Reservoir effective reservoir quantitative classification, it is wrapped
Include following steps:
Obtain the reservoir quality parameter with reference to volcanic rock reservoir;
Split clustering procedure using K-means and cluster analysis is carried out to the reservoir quality parameter, to establish volcanic rock reservoir
The criteria for classification of effective reservoir;
Based on the criteria for classification of the Volcanic Reservoir effective reservoir, all kinds of effectively storages are established using linear discriminant function method
The linear discriminant function of layer;
Based on the linear discriminant function, the type of target volcanic rock reservoir is determined.
Volcanic rock reservoir is different from conventional sandstone reservoir, its reservoir quality is mainly by matrix pores and crack dual media
Influence.Technical solution provided by the invention is combined cluster is split with linear discriminant analysis method, preferably reservoir quality parameter,
Quantitative classification algorithm is established, ensures data order and stability, so as to carry out quick evaluation of classification to reservoir.
In technical solution provided by the invention, the effective reservoir, which refers to possess, to be preserved with seepage flow movable fluid (with hydrocarbon
Based on class fluid) ability, there is the reservoir of commercial mining value under technology at this stage and economic condition.
Reservoir quality is primarily referred to as the ability that reservoir preserves and is percolated fluid.In technical solution provided by the invention, institute
Stating reservoir quality parameter includes porosity, permeability, density, fluidized bed index specific productivity index.These parameters not only allow for
Electrical parameter otherness between volcanic rock reservoir difference Volcanic uranium deposit, while rock seepage flow is reflected using FZI fluidized beds index
Ability, and emphasis considers the output capacity of different type Volcanic uranium deposit, and the dynamic and static parameter of summary is as volcanic rock
The favourable parameters of reservoir effective reservoir classification.These parameters can be obtained by existing rock core information.
In the above-mentioned methods, it is preferable that the fluidized bed index FZI is counted after being deformed by Kozeny-carman equations
Obtain, specifically as shown in formula 1 to formula 3:
By formula 1 to formula 3, can be derived from:Lg RQI=lg FZI+lg Φ z.
For formula 1 into formula 3, RQI is oil reservoir qualitative index;Φ z are the ratio between pore volume and particle volume;FZI is fluidized bed
Index;K is permeability (× 10-3μm2);Ф is porosity (decimal).
In the above-mentioned methods, it is preferable that split clustering procedure using K-means and cluster point is carried out to the reservoir quality parameter
Analysis includes:
Determine classification number k;
Sample point is categorized into k cluster, so that distance ratio of central point of cluster arrives belonging to each sample point to its
The distance of the central point of other clusters is small.
Technical solution provided by the invention splits clustering procedure to the target volcanic rock reservoir mass parameter using K-means
Cluster analysis is carried out, wherein, the core algorithm thought of K-means segmentation clustering procedures is exactly that n observation Exemplary classes are poly- to k
In class so that each central point of the observed value where it than other cluster centre points apart from smaller, apart from smaller, get over phase
Seemingly, diversity factor is smaller.On the premise of classification number k is provided, K-means can quickly handle density data collection, its center is counted
Method improves the discrete and mixing of serial number attribute and clusters, preferable for numerous and jumbled geologic data clustering effect, and should
Method influenced by singular value, Similar measure and inappropriate clustering variable it is smaller, can be into for inappropriate preliminary classification
Row adjusts iteration repeatedly, and is finally reached convergence.
In the above-mentioned methods, classification number k can be set according to the needs of goal in research, technical side provided by the invention
K is set to 3 by case, and expression can be divided into 3 classes, can represent respectively, in, it is poor.
In the above-mentioned methods, it is preferable that sample point is categorized into k cluster, so that poly- belonging to each sample point to its
The distance of central point of the distance of the central point of class than being clustered to other is small can to include procedure below:
It is assumed that sample set is { x(1), x(2)... x(m), x(i)∈R(n), according to the Euclidean distance calculation formula meter shown in formula 4
Calculate distance
K cluster centre is chosen, the affiliated classification c of each sample i is calculated according to the formula shown in formula 5(i)
For every a kind of central point μj, as shown in Equation 6, final class central point can be calculated by renewal
Recurring formula 5 and formula 6, until the quadratic sum of the Euclidean distance of the class central point of cluster belonging to all samples to its
Minimum, that is, it is minimum to reach in formula 7 function J (c, μ), and finally restrains, so as to obtain the poly- of global optimum or local optimum
Class gathering
In the above-mentioned methods, it is preferable that the criteria for classification of the volcanic rock reservoir effective reservoir such as table 1-1 and table 1-2 institutes
Show
Table 1-1
Note:Data above table 1-1 average values are minimum value to maximum;Such as the porosity of I class effective reservoirs is most
Small value is 5.6, maximum 19.3, average value 11.3.
Table 1-2
In the above-mentioned methods, the effective reservoir of the volcanic rock reservoir is divided into 3 classes, is I classes effective reservoir, II classes respectively
Effective reservoir and Group III effective reservoir;Wherein, the lithology of the I classes effective reservoir is mainly by pacifying mountain matter ignimbrite, dust
One or more of compositions in lava and stomata andesite;The lithology of the II classes effective reservoir is mainly by almond basalt, peace
One or more of compositions in mountain rock, rhyolite, peace mountain matter vent breccia and dust tufa stone;The Group III effective reservoir
Lithology mainly by tufa stone, dacite, English peace matter andesite and fine grain rhyolite in one or more form.
Volcanic Facies can be typically divided between explosion facies, overflow facies, volcanic vent phase and volcanic deposit phase.Wherein, volcano
Rock overflow facies can be generally divided into, in and the parfacies of lower part three;Volcanic rock explosion facies are mainly that tephre is various forms of
Assembly, can be typically divided between sky according to volcaniclastic type and combination and falls parfacies, hot radical wave parfacies, hot chip stream
Parfacies and splash down parfacies.The lithofacies that can be seen that I class effective reservoirs from table 1-2 are mainly overflow facies-top parfacies;II classes have
The lithofacies for imitating reservoir are mainly overflow facies-lower part parfacies, alternatively, explosion facies;The lithofacies of Group III effective reservoir are predominantly broken out
Phase-hot chip stream parfacies, alternatively, overflow facies-middle part parfacies (densification).
In the criteria for classification of above-mentioned Volcanic Reservoir effective reservoir, the reservoir quality of the I classes effective reservoir, seepage flow energy
Power is fabulous, and formation testing specific productivity index is also higher.Compared with the I classes effective reservoir, the reservoir properties of the II classes effective reservoir
Take second place, FZI values are opposite to diminish, and specific productivity index is then not much different.Compared with the II classes effective reservoir, the Group III is effective
The reservoir properties of reservoir are deteriorated, but still have certain seepage flow and output capacity after pressure break, between potentiality effectively and effective reservoir
Between.What it is less than Group III criteria for classification is then invalid reservoir, and fracturing production does not go out substantially or pumping only water outlet, and liquid measure also compared with
It is low.
In general, effective reservoir just refer to possess preserve with seepage flow movable fluid (based on hydrocarbon fluid) ability,
There is the reservoir that commercial mining is worth under technology and economic condition at this stage.And invalid reservoir is then primarily referred to as reservoir thing
Poor, the reservoir without liquid-producing capacity or less than dried layer standard of property.
Technical solution provided by the invention based on cluster analysis as a result, using linear discriminant function method establish it is all kinds of effectively
The linear discriminant function of reservoir.Linear discriminant function is the criterion function for being used directly to classify to pattern, and linear discriminant divides
The purpose of analysis is that the low-dimensional feature for most having classification information is found out from high-dimensional feature space, these features can will be same category of
Sample flocks together, and makes different classes of sample separated as much as possible.The following institute of general type of N-dimensional linear discriminant function
Show:
D (X)=w1x1+w2x2+…+wnxn+w0=WTX+w0In formula, X=[x1, x2..., xn]T, W=[w1, w2...,w]TFor
Weight vector;
To M linear separability pattern class, can be usedDichotomy method, will belong to ω with linear discriminant functioniThe mould of class
Formula is not belonging to ω with remainingiThe pattern of class separates.
In the above-mentioned methods, when classification is identified using linear discriminant function, if only di(X)>0, other d (X) are equal<
0, then it is judged to ωiClass;Wherein,
In the above-mentioned methods, it is preferable that the linear discriminant function of all kinds of effective reservoirs is as shown in table 2.
Table 2
The type of effective reservoir | Linear discriminant function |
I | Y1=0.145 × FZI+229.9 × DEN+0.5 × Jos-273.02 |
II | Y2=-1.423 × FZI+238.9 × DEN+0.15 × Jos-282.53 |
III | Y3=-1.184 × FZI+241.6 × DEN+0.03 × Jos-300.66 |
In the above-mentioned methods, it is preferable that based on the linear discriminant function, determining the type of target volcanic rock reservoir includes
Procedure below:
The reservoir quality parameter of target volcanic rock reservoir is obtained, and is substituted into the linear discriminant function;
As the functional value > 0 of the linear function of acquisition, determine that the type of target volcanic rock reservoir belongs to the linear function institute
The type of representative;Otherwise, it determines it is not belonging to the type representated by linear function;
When acquisition linear function value≤0, determine to belong to invalid reservoir.
In a detailed embodiment, based on the linear discriminant function, determine that the type of target volcanic rock reservoir can
With including procedure below:
By the reservoir quality parameter of target volcanic rock reservoir, linear discriminant function Y is substituted into1In, such as functional value of acquisition>0,
Determine to belong to Y1Representative type;Otherwise, it determines it is not belonging to Y1Representative type, substitutes into next linear discriminant function at this time
Y2In, such as functional value of acquisition>0, determine to belong to Y2Representative type;Otherwise, it determines it is not belonging to Y2Representative type, this
When substitute into next linear discriminant function Y3In, such as functional value of acquisition>0, determine to belong to Y3Representative type;Otherwise, it determines not
Belong to Y2Representative type, but belong to invalid reservoir.
In the above-mentioned methods, it is preferable that the reservoir quality parameter of the target volcanic rock reservoir and the linear discriminant letter
Variable in number is corresponding.It is highly preferred that the reservoir quality parameter of the target volcanic rock reservoir includes fluidized bed index, density
And specific productivity index.
Beneficial effects of the present invention:
Technical solution provided by the invention has important theory and practical significance, can be further efficient for such oil reservoir
Development deployment provides technical support, and can be the perfect offer of volcanic rock reservoir effective reservoir evaluation of classification standard and technical system
Favourable reference.
Brief description of the drawings
Fig. 1 is the typical clustering and discriminant classification chart that embodiment of the present invention provides;
Fig. 2 is the classification chart of certain oil field individual well volcanic reservoirs;
Fig. 3 is the classification plan of certain oil field volcanic reservoirs;
Fig. 4 is 1 year cumulative production prediction plan of certain oil field volcanic rock reservoir;
Fig. 5 is certain oil field volcanic rock reservoir individual well Production development curve map.
Embodiment
In order to which technical characteristic, purpose and the beneficial effect of the present invention is more clearly understood, now to the skill of the present invention
Art scheme carry out it is described further below, but it is not intended that to the present invention can practical range restriction.
Present embodiments provide for a kind of method of Volcanic Reservoir effective reservoir quantitative classification.This method is with certain volcanic rock
The effective reservoir of the 20 mouthfuls of individual wells in oil field is research object, and 250 samples therein are studied.Detailed process is as follows:
Step S101:Obtain the reservoir quality parameter with reference to volcanic rock reservoir
It can include porosity, permeability, density, fluidized bed index and ratio with reference to the reservoir quality parameter of volcanic rock reservoir
Productivity index;These parameters can be obtained by existing rock core information.
It is preferred that collecting volcanic rock reservoir mass parameter sample, parameter sample set is formed, wherein, fluidized bed index FZI can
With by being calculated after the deformation of Kozeny-carman equations, calculating process is as shown in formula 1 to formula 3:
Formula 1 is into formula 3:RQI is oil reservoir qualitative index;Φ z are the ratio between pore volume and particle volume;FZI is fluidized bed
Index;K is permeability (× 10-3μm2);Ф is porosity (decimal).
Step S102:Split clustering procedure using K-means and cluster analysis is carried out to the reservoir quality parameter, to establish fire
The criteria for classification of mountain rock reservoir effective reservoir
First, priori classification number k is provided, k=3 in the present embodiment.
Then, by sample classification to be clustered into k cluster, so that the central point of cluster belonging to each sample to its
The distance of central point of the distance than being clustered to other is small, and process is as follows:
It is assumed that sample set to be clustered is { x(1), x(2)... x(m), x(i)∈R(n), according to the Euclidean distance meter shown in formula 4
Calculate formula and calculate distance
The affiliated classification c of each sample i is calculated according to the formula shown in formula 5(i)
For every a kind of central point μj, as shown in Equation 6, final class central point can be calculated by renewal
Recurring formula 5 and formula 6, until the quadratic sum of the Euclidean distance at the center of cluster belonging to all samples to its is most
It is small, that is, reach function J (c, μ) value minimum in formula 7, and finally restrain, so as to obtain the clustering cluster of global optimum or local optimum
Collection;Wherein, the calculation formula of the J (c, μ) is as shown in Equation 7
Table 1-1
Note:Data above table 1-1 average values are minimum value to maximum;Such as the porosity of I class effective reservoirs is most
Small value is 5.6, maximum 19.3, average value 11.3.
Table 1-2
Using K-means Fast Segmentation clustering method statistical analysis parameter samples, volcanic rock reservoir effective reservoir point is established
Class evaluation criterion, as shown in table 1-1 and table 1-2.
Step S103:Based on the criteria for classification of the Volcanic Reservoir effective reservoir, established using linear discriminant function method
The linear discriminant function of all kinds of effective reservoirs
Established with fluidized bed index FZI, density DEN, specific productivity index Jos as variable using linear discriminant function method
The linear discriminant function of all kinds of effective reservoirs, the results are shown in Table 2
Table 2
The type of effective reservoir | Linear discriminant function |
I | Y1=0.145 × FZI+229.9 × DEN+0.5 × Jos-273.02 |
II | Y2=-1.423 × FZI+238.9 × DEN+0.15 × Jos-282.53 |
III | Y3=-1.184 × FZI+241.6 × DEN+0.03 × Jos-300.66 |
Step S104:Based on the linear discriminant function, the type of target volcanic rock reservoir is determined.
It is corresponding with the variable in linear discriminant function, the reservoir quality parameter of target volcanic rock reservoir is obtained, these ginsengs
Number includes density, fluidized bed index and specific productivity index;
The mass parameter of target volcanic rock reservoir is substituted into the linear discriminant function shown in table 2;
As the functional value > 0 of the linear function of acquisition, determine to belong to the type representated by the linear function;Otherwise, do not belong to
In the type representated by the linear function;
When acquisition linear function functional value≤0, determine to belong to invalid reservoir.
Differentiate result as shown in Fig. 1 and table 3.
Table 3
Type | Cluster sample |
I | 6 |
II | 201 |
III | 43 |
Table 4
Application analysis
In the 20 mouthfuls of individual well effective reservoir quantitative assessments in certain volcanic rock oil field classification, the volcanic rock that will provide according to the present invention
It is that the method for oil reservoir effective reservoir quantitative classification obtains as a result, the classification results obtained with individual well formation testing pilot production carry out contrast tests
Card, the results are shown in Table 4.In 250 samples of discriminant analysis, Volcanic Reservoir effective reservoir quantitative classification provided by the invention
The right judging rate of method reach 92%, it is seen that the result reliability of technical solution provided by the invention is high, can be with to non-core hole
Phase reservoir quality parameter (with fluidized bed index FZI, density DEN, specific productivity index Jos) is asked for according to well log interpretation, is substituted into each
The linear discriminant function of class effective reservoir carries out the Classification and Identification of effective reservoir, can rapidly and accurately to individual well volcanic rock reservoir into
The quantitative division (attached drawing 2, attached drawing 3) of row.
The present invention is according to different well point volcanic rock interval effective reservoir quantitative classifications and its thickness statistical result, comprehensive production
Dynamic data, takes production forecast time term, final calculating and plotting draws research area's effective reservoir production forecast distribution map, excellent
Choosing determines development deployment dominant area (attached drawing 4).As can be seen that research area the north J214-J215 well area lists from prognostic chart
One term of well tire out production it is relatively low, the current stage be not suitable for guide exploitation, J208, JL2008, J213, J204 well area individual well one
It is higher that term tires out production prediction, can the of a relatively high J208 wells fault block of prioritized deployment well control degree and J204 well areas remaining region conducts
REGION OF WATER INJECTION OILFIELD is disposed in lower step rolling test.
So far, deployment 1 mouthful of water horizontal well of finishing drilling in the J208 well areas of prioritized deployment, oil-producing 33.5 at initial stage day of going into operation
Ton day, produces 230 days, adds up 7700 tons of oil-producing.1 mouthful of JL2011 well areas deployment finishing drilling straight well, oil-producing 13.5 at initial stage day of going into operation
Ton day, produces 584 days, adds up 6311 tons of oil-producing, average 11 ton day of individual well day oil-producing.Horizontal well daily output is that straight well is made in Japan
(attached drawing 5) 1.9-5.8 again.As it can be seen that on the basis of volcanic rock reservoir effective reservoir quantitative judge with classification understanding, can be effectively
Instruct the efficient well site deployment in scene.
Shown by application analysis, classification results of the present invention and live actual effect matching degree are high, are volcanic rock
The deployment of reservoir Efficient Development provides strong technical support, has very strong practicality.
It is comprehensive volcanic rock reservoir lithology lithofacies and its electrical response parameter, microcosmic in conclusion method provided by the invention
Feature, fracture development parameter and Production development parameter, are known using clustering method and typical linear techniques of discriminant analysis scheduling theory
Know and examined with oil reservoir expert, being association of activity and inertia realizes that reservoir quality parameter is preferred, and the volcanic rock reservoir effective reservoir of foundation quantifies
Distinguished number, has preliminarily formed effective reservoir quickly identification and classification, to instruct the Efficient Development of volcanic rock reservoir.
Claims (10)
1. a kind of method of Volcanic Reservoir effective reservoir quantitative classification, it comprises the following steps:
Obtain the reservoir quality parameter with reference to volcanic rock reservoir;
Split clustering procedure using K-means and cluster analysis is carried out to the reservoir quality parameter, it is effective to establish Volcanic Reservoir
The criteria for classification of reservoir;
Based on the criteria for classification of the Volcanic Reservoir effective reservoir, all kinds of effective reservoirs are established using linear discriminant function method
Linear discriminant function;
Based on the linear discriminant function, the type of target volcanic rock reservoir is determined.
2. according to the method described in claim 1, wherein, the reservoir quality parameter with reference to volcanic rock reservoir includes hole
Degree, permeability, density, fluidized bed index and specific productivity index.
3. according to the method described in claim 2, wherein, the calculation formula of the fluidized bed index is as shown in formula 1 to formula 3
In formula 1- formulas 3, RQI is oil reservoir qualitative index;Φ z are the ratio between pore volume and particle volume;FZI is fluidized bed index;K
For permeability;Ф is porosity.
4. according to the method described in claim 1, wherein, using K-means split clustering procedure to the reservoir quality parameter into
Row cluster analysis includes:
Determine classification number k;
Sample point is categorized into k cluster, so that distance ratio of central point of cluster arrives other belonging to each sample point to its
The distance of the central point of cluster is small.
5. according to the method described in claim 4, wherein, the classification number k=3.
6. according to the method described in claim 1, wherein, the criteria for classification such as table 1-1 of the Volcanic Reservoir effective reservoir and
Shown in table 1-2:
Table 1-1
Table 1-2
7. according to the method described in claim 6, wherein, the linear discriminant function of all kinds of effective reservoirs is respectively:
Y1=0.145 × FZI+229.9 × DEN+0.5 × Jos-273.02;
Y2=-1.423 × FZI+238.9 × DEN+0.15 × Jos-282.53;
Y3=-1.184 × FZI+241.6 × DEN+0.03 × Jos-300.66;Wherein,
Y1Represent the linear discriminant function of I class effective reservoirs, Y2Represent the linear discriminant function of II class effective reservoirs, Y3Represent III
The linear discriminant function of class effective reservoir;In Y1、Y2And Y3In, FZI represents fluidized bed index, and DEN represents density, and Jos represents ratio
Productivity index.
8. the method according to claim 1 or 7, wherein, based on the linear discriminant function, determine target volcanic rock reservoir
Type include procedure below:
The reservoir quality parameter of target volcanic rock reservoir is obtained, and is substituted into the linear discriminant function;
As the functional value > 0 of the linear function of acquisition, determine to belong to the type representated by the linear function;Otherwise, it is not belonging to this
Type representated by linear function;
When acquisition linear function functional value≤0, determine to belong to invalid reservoir.
9. according to the method described in claim 8, wherein, the reservoir quality parameter of the target volcanic rock reservoir with it is described linear
Variable in discriminant function is corresponding.
10. according to the method described in claim 9, wherein, the reservoir quality parameter of the target volcanic rock reservoir includes flowing
Layer index, density and specific productivity index.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711275674.2A CN107965318A (en) | 2017-12-06 | 2017-12-06 | A kind of method of Volcanic Reservoir effective reservoir quantitative classification |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711275674.2A CN107965318A (en) | 2017-12-06 | 2017-12-06 | A kind of method of Volcanic Reservoir effective reservoir quantitative classification |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107965318A true CN107965318A (en) | 2018-04-27 |
Family
ID=61998399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711275674.2A Pending CN107965318A (en) | 2017-12-06 | 2017-12-06 | A kind of method of Volcanic Reservoir effective reservoir quantitative classification |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107965318A (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109555518A (en) * | 2018-12-14 | 2019-04-02 | 中国石油大港油田勘探开发研究院 | A kind of alluvial fan individual well configuration recognition methods based on cluster and distinguished number |
CN110826936A (en) * | 2019-11-22 | 2020-02-21 | 中国地质大学(北京) | Shale oil and gas resource grading evaluation method |
CN110969273A (en) * | 2018-09-28 | 2020-04-07 | 北京国双科技有限公司 | Production well yield prediction method and device |
CN110987751A (en) * | 2019-11-15 | 2020-04-10 | 东北石油大学 | Quantitative grading evaluation method for pore throat of compact reservoir in three-dimensional space |
CN111206921A (en) * | 2018-11-22 | 2020-05-29 | 中石化石油工程技术服务有限公司 | Description method suitable for favorable reservoir stratum of volcanic overflow phase |
CN112085109A (en) * | 2020-09-14 | 2020-12-15 | 电子科技大学 | Phase-controlled porosity prediction method based on active learning |
CN112348330A (en) * | 2020-10-27 | 2021-02-09 | 中国海洋石油集团有限公司 | Method for determining quantitative classification of altered pyroclastic glutenite reservoir |
CN112541607A (en) * | 2019-09-23 | 2021-03-23 | 中国石油天然气股份有限公司 | Volcanic oil reservoir productivity prediction method and device |
CN112983406A (en) * | 2021-03-15 | 2021-06-18 | 西南石油大学 | Natural gas hydrate reservoir parameter index evaluation method |
CN113688922A (en) * | 2021-08-31 | 2021-11-23 | 西南石油大学 | Fluid type identification method of agglomeration clustering unsupervised algorithm |
-
2017
- 2017-12-06 CN CN201711275674.2A patent/CN107965318A/en active Pending
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110969273B (en) * | 2018-09-28 | 2020-11-13 | 北京国双科技有限公司 | Production well yield prediction method and device |
CN110969273A (en) * | 2018-09-28 | 2020-04-07 | 北京国双科技有限公司 | Production well yield prediction method and device |
CN111206921A (en) * | 2018-11-22 | 2020-05-29 | 中石化石油工程技术服务有限公司 | Description method suitable for favorable reservoir stratum of volcanic overflow phase |
CN109555518B (en) * | 2018-12-14 | 2022-07-05 | 中国石油天然气股份有限公司 | Alluvial fan single well configuration identification method based on clustering and discrimination algorithm |
CN109555518A (en) * | 2018-12-14 | 2019-04-02 | 中国石油大港油田勘探开发研究院 | A kind of alluvial fan individual well configuration recognition methods based on cluster and distinguished number |
CN112541607A (en) * | 2019-09-23 | 2021-03-23 | 中国石油天然气股份有限公司 | Volcanic oil reservoir productivity prediction method and device |
CN112541607B (en) * | 2019-09-23 | 2024-05-28 | 中国石油天然气股份有限公司 | Volcanic oil reservoir productivity prediction method and device |
CN110987751A (en) * | 2019-11-15 | 2020-04-10 | 东北石油大学 | Quantitative grading evaluation method for pore throat of compact reservoir in three-dimensional space |
CN110826936A (en) * | 2019-11-22 | 2020-02-21 | 中国地质大学(北京) | Shale oil and gas resource grading evaluation method |
CN110826936B (en) * | 2019-11-22 | 2022-09-20 | 中国地质大学(北京) | Shale oil and gas resource grading evaluation method |
CN112085109A (en) * | 2020-09-14 | 2020-12-15 | 电子科技大学 | Phase-controlled porosity prediction method based on active learning |
CN112348330A (en) * | 2020-10-27 | 2021-02-09 | 中国海洋石油集团有限公司 | Method for determining quantitative classification of altered pyroclastic glutenite reservoir |
CN112348330B (en) * | 2020-10-27 | 2022-12-13 | 中国海洋石油集团有限公司 | Method for determining quantitative classification of altered pyroclastic glutenite reservoir |
CN112983406B (en) * | 2021-03-15 | 2022-03-25 | 西南石油大学 | Natural gas hydrate reservoir parameter index evaluation method |
CN112983406A (en) * | 2021-03-15 | 2021-06-18 | 西南石油大学 | Natural gas hydrate reservoir parameter index evaluation method |
CN113688922A (en) * | 2021-08-31 | 2021-11-23 | 西南石油大学 | Fluid type identification method of agglomeration clustering unsupervised algorithm |
CN113688922B (en) * | 2021-08-31 | 2023-07-25 | 西南石油大学 | Fluid type identification method of condensation clustering unsupervised algorithm |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107965318A (en) | A kind of method of Volcanic Reservoir effective reservoir quantitative classification | |
CN109061765B (en) | Trap evaluation method for heterogeneous thin sandstone interbed reservoir | |
CN104965979B (en) | A kind of tight sand effective reservoir recognition methods | |
CN110276827B (en) | Effectiveness evaluation method based on shale reservoir | |
CN105242026B (en) | A kind of gas reservoir source title method | |
CN104020509B (en) | Dam, chiltern beach based on Bayes discriminant analysis sedimentary micro Logging Identification Method | |
CN105160414B (en) | Predict the method and device of full oil reservoir producing region type | |
CN103510947B (en) | The method for building up of dam, beach sandstone microfacies recognition mode | |
Zifei et al. | A study on remaining oil distribution in a carbonate oil reservoir based on reservoir flow units | |
CN105572747B (en) | A method of identification waterflooding reservoir clastic rock lithology in the areas Fu Jia with high salt | |
CN105447762B (en) | A kind of calculation method of the low-permeability oil deposit water logging information of fluid replacement | |
CN106501870B (en) | A kind of lacustrine facies densification shelly limestone comparative good-quality reservoir stratum identification method | |
CN112147301A (en) | Quantitative evaluation method for effectiveness of continental-phase fresh water lake basin compact oil hydrocarbon source rock | |
CN107505663B (en) | A kind of method for building up of carbonate reservoir classification plate and application | |
CN109031424A (en) | A method of based on well logging Multiparameter low permeability reservoir Diagenetic Facies | |
Wang et al. | Determine level of thief zone using fuzzy ISODATA clustering method | |
CN109403960B (en) | Method for judging reservoir fluid properties by logging gas peak-logging state | |
CN108805158A (en) | A kind of fine and close oily reservoir diagenetic phase division methods | |
Orodu et al. | Hydraulic (flow) unit determination and permeability prediction: a case study of block Shen-95, Liaohe Oilfield, North-East China | |
CN110656922B (en) | Shale isochronous stratum logging dividing method and system based on pencils and stone belt characteristics | |
CN111460725A (en) | Shale gas sweet spot prediction based on multi-level fuzzy recognition | |
CN105259576B (en) | A kind of oil-gas reservoir identification method using earthquake statistics feature | |
CN115387785A (en) | Sea-facies carbonate-cuttings limestone reservoir high-permeability strip identification method and device | |
CN114076992A (en) | Classification evaluation method for micro pore throats of meandering stream reservoir based on lithofacies | |
CN114139242A (en) | Water flooded layer well logging evaluation method based on lithofacies |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180427 |
|
WD01 | Invention patent application deemed withdrawn after publication |