CN104818982B - Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir - Google Patents
Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir Download PDFInfo
- Publication number
- CN104818982B CN104818982B CN201510187396.XA CN201510187396A CN104818982B CN 104818982 B CN104818982 B CN 104818982B CN 201510187396 A CN201510187396 A CN 201510187396A CN 104818982 B CN104818982 B CN 104818982B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- fracture
- crack
- attribute
- relation
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004458 analytical method Methods 0.000 claims abstract description 57
- 238000004445 quantitative analysis Methods 0.000 claims abstract description 41
- 238000004891 communication Methods 0.000 claims abstract description 24
- 238000011161 development Methods 0.000 claims abstract description 14
- 238000001615 p wave Methods 0.000 claims description 35
- 238000000605 extraction Methods 0.000 claims description 7
- 230000004927 fusion Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims 2
- 208000010392 Bone Fractures Diseases 0.000 description 125
- 206010017076 Fracture Diseases 0.000 description 125
- 238000010586 diagram Methods 0.000 description 11
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 239000011435 rock Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 230000003068 static effect Effects 0.000 description 3
- 208000006670 Multiple fractures Diseases 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 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
- 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
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 And Detection Of Objects (AREA)
Abstract
The invention provides a method and a device for quantitatively analyzing connectivity of a fracture-cavity reservoir stratum, wherein the method comprises the following steps: carrying out inversion by using seismic and logging data to obtain a longitudinal wave impedance body, and generating a porosity body according to the relation between logging porosity and longitudinal wave impedance; determining an effective reservoir by taking the lower limit value of the porosity of the effective reservoir as a critical value; judging the connectivity of the effective reservoirs according to the contact relation between the effective reservoirs, and determining the boundary of each fracture-cavity reservoir body; extracting an initial seismic fracture prediction attribute body by utilizing a coherence or curvature algorithm; determining the connectivity of the reservoirs among wells by a dynamic connectivity analysis method, and judging the connectivity relation among the fracture-cavity reservoirs according to the boundaries of the fracture-cavity reservoirs; determining the crack development degree according to the communication relation, acquiring a crack attribute threshold range, and determining an earthquake crack prediction attribute body; and generating a reservoir connectivity quantitative analysis result according to the contact relation between the fracture-cavity reservoir body and the seismic fracture prediction attribute body and the fracture attribute threshold range.
Description
Technical field
The present invention is on reservoir distribution Predicting Technique, specifically connective quantitative point on a kind of fractured-vuggy reservoir
Analyse method and device.
Background technology
Fracture-cavity units division is developed significant to Carbonate Reservoir, is how to judge to store up the problem of its most critical
Connectedness between layer.It can determine that it is connective using dynamic connectivity analysis method between well, specifically include oil reservoir
Production interference analytic approach, fluid properties discriminant analysis, inter-well interference well testing method and tracer method between pressure trend analytic approach, well
Deng.Reservoir beyond well control range, then need to predict that it is connective by data such as earthquakes.In view of influence fractured-vuggy reservoir
Two connective Dominated Factors are reservoir itself and crack, and conventional method is:Reservoir and crack are carried out using seismic properties
Prediction, by judging the contact relation between reservoir and crack, reach connective prediction purpose.Due to reservoir and FRACTURE PREDICTION master
If being carried out under static well data (rock core, imaging logging etc.) constraint, often lead to connective prediction result and bored with real later
Well inter-well dynamic connectivity analysis result is inconsistent.
For the above situation, existing Predicting Technique is that reservoir and FRACTURE PREDICTION are carried out under the constraint of connectedness.Specifically
It is to utilize seismic amplitude attribute forecast reservoir, a kind of AFE (main flow FRACTURE PREDICTION technology) predictions crack, both constantly regulates threshold
Value, makes reservoir and results of fracture prediction be matched with Production development, so that the reservoir without well controlling zone domain and FRACTURE PREDICTION knot
Fruit is more accurate.This method is specifically situation about being determined by which factor (reservoir or crack) not distinguishing reservoir communication
Under, it is general by dynamic connectivity analysis result, while constrain reservoir and FRACTURE PREDICTION, i.e., both adjustment threshold value (such as Fig. 1 simultaneously
It is shown).But the problem of above-mentioned analysis method is present mainly has:(1) multi-solution is too strong, is the equal of (dynamic with a known variables
State connectedness result) two known variables of constraint solving (amplitude and AFE attribute thresholds) are gone, it is probably specifically adjustment amplitude threshold
It is more, or adjustment AFE threshold values are more, or the two adjusts simultaneously, can be consistent with dynamic analysis result.(2) due to strong
It is most important connection factor to adjust reservoir itself, and next is only crack, does not analyse in depth the connectedness in crack, this also influences
The precision of final reservoir communication prediction result.
The content of the invention
For background technology multi-solution it is strong, fracture connection Journal of Sex Research not enough gos deep into the problems such as, the embodiment of the present invention
Main purpose is to provide a kind of fractured-vuggy reservoir connectedness quantitative analysis method and device, to reduce reservoir and fractue spacing point
The multi-solution of result is analysed, so as to improve the degree of accuracy of reservoir communication analysis.
To achieve these goals, the embodiment of the present invention provides a kind of fractured-vuggy reservoir connectedness quantitative analysis method, institute
The method of stating includes:Inverting is carried out using earthquake and log data, p-wave impedance body is obtained, according to well logging porosity and p-wave impedance
Between relation, the p-wave impedance body is converted into porosity cube;With in the porosity cube that is obtained by well log interpretation
The lower limit of effective reservoir porosity determines effective reservoir as critical value;According to the contact relation between the effective reservoir
Judge the connectedness of the effective reservoir, determine the border of each fracture and cave reservoir;It is initial using relevant or curvature algorithm, extraction
Seismic prediction of fracture attribute volume;The connectedness of reservoir between well is determined by dynamic connectivity analysis method, and according to described seam
Hole Reservoir Body border judges the connected relation of the fracture and cave reservoir between any two;Fracture development is determined according to the connected relation
Degree, obtains the scope of the crack attribute thresholds of the primary earthquake FRACTURE PREDICTION attribute volume, and determines Seismic prediction of fracture category
Property body;According to the contact relation and the crack attribute thresholds of the fracture and cave reservoir and the Seismic prediction of fracture attribute volume
Scope, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results.
In one embodiment, the relation between above-mentioned well logging porosity and p-wave impedance is:With the increase of porosity,
The linear relationship that p-wave impedance accordingly reduces.
In one embodiment, the connectedness of the effective reservoir is judged according to the contact relation between the effective reservoir,
Reservoir Body border is determined, including:Judge with the presence or absence of contact, the contact to include between the effective reservoir:Point contact, line connect
Touch face contact;If contact between the effective reservoir be present, judge between the effective reservoir for connection;Otherwise, sentence
It is not connect between the fixed effective reservoir.
In one embodiment, development degree of micro cracks in oil is determined according to the connected relation, it is pre- obtains the primary earthquake crack
The scope of the crack attribute thresholds of attribute volume is surveyed, including:When the connected relation for connection when, in units of a constant, from big
To the small numerical value for gradually reducing the crack attribute thresholds, until there is crack between two fracture and cave reservoirs, and connect
For a fracture-cavity units, the higher limit using the crack attribute thresholds now as the crack attribute thresholds;When the company
Clearance system is when not connecting, in units of the constant, gradually to increase the numerical value of the crack attribute thresholds from small to large, until
The fracture and cave reservoir separates, the lower limit using the crack attribute thresholds now as the crack attribute thresholds;According to
The higher limit and lower limit of the crack attribute thresholds obtain the scope of the crack attribute thresholds.
In one embodiment, according to the contact relation and institute of the fracture and cave reservoir and primary earthquake FRACTURE PREDICTION attribute volume
The scope of crack attribute thresholds is stated, judges the connectedness of fractured-vuggy reservoir, generates reservoir communication quantitative analysis results, including:
The effective reservoir is merged to one fused data body of generation with Seismic prediction of fracture attribute volume, inside the fused data body
Contact relation, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results.
The embodiment of the present invention also provides a kind of fractured-vuggy reservoir connectedness quantitative analysis device, and described device includes:Hole
Degree body generation unit, for carrying out inverting using earthquake and log data, p-wave impedance body is obtained, according to well logging porosity with indulging
Relation between wave impedance, the p-wave impedance body is converted into porosity cube;Effective reservoir determining unit, for being surveyed with passing through
Well explains that the lower limit of effective reservoir porosity in the porosity cube obtained as critical value, determines effective reservoir;Fracture hole
Reservoir Body border determining unit, for judging the connection of the effective reservoir according to the contact relation between the effective reservoir
Property, determine the fracture and cave reservoir border of each Reservoir Body;Primary earthquake FRACTURE PREDICTION attribute volume extraction unit, it is relevant for utilizing
Or curvature algorithm, extract primary earthquake FRACTURE PREDICTION attribute volume;Connected relation judging unit, for being analyzed by dynamic connectivity
Method determines the connectedness of reservoir between well, and judges the fracture and cave reservoir between any two according to described fracture and cave reservoir border
Connected relation;Threshold range acquiring unit, for determining development degree of micro cracks in oil according to the connected relation, obtain described initial
The scope of the crack attribute thresholds of Seismic prediction of fracture attribute volume, and determine Seismic prediction of fracture attribute volume;Analysis result generates
Unit, for the contact relation according to the fracture and cave reservoir and the Seismic prediction of fracture attribute volume and the crack attribute threshold
The scope of value, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results.
In one embodiment, the relation between above-mentioned well logging porosity and p-wave impedance is:With the increase of porosity,
The linear relationship that p-wave impedance accordingly reduces.
In one embodiment, above-mentioned fracture and cave reservoir border determining unit is specifically used for:Judge the effective reservoir it
Between with the presence or absence of contact, the contact includes:Point contact, linear contact lay and face contact;Connect if existed between the effective reservoir
Touch, then judge between the effective reservoir for connection;Otherwise, it is determined that it is not connect between the effective reservoir.
In one embodiment, above-mentioned threshold range acquiring unit is specifically used for:When the connected relation is connects, with
One constant is unit, gradually reduces the numerical value of the crack attribute thresholds from big to small, until two fracture and cave reservoirs it
Between there is crack, and it is a fracture-cavity units to connect, using the crack attribute thresholds now as the crack attribute thresholds
Higher limit;When the connected relation is not connect, in units of the constant, gradually increase the crack category from small to large
Property threshold value numerical value, until the fracture and cave reservoir separates, using the crack attribute thresholds now as the crack attribute
The lower limit of threshold value;The scope of the crack attribute thresholds is obtained according to the higher limit of the crack attribute thresholds and lower limit.
In one embodiment, above-mentioned analysis result generation unit is specifically used for:By the effective reservoir and the earthquake
FRACTURE PREDICTION attribute volume fusion one fused data body of generation, according to the contact relation inside the fused data body, judges fracture hole
The connectedness of type reservoir, generate reservoir communication quantitative analysis results.
The beneficial effects of the present invention are due to reference to many factors, and the qualitative threshold value model for calculating prediction attribute volume
Enclose, the multi-solution of reservoir and fractue spacing analysis result can be reduced, so as to improve the degree of accuracy of reservoir communication analysis.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, embodiment will be described below
In the required accompanying drawing used be briefly described, it should be apparent that, drawings in the following description be only the present invention some
Embodiment, for those of ordinary skill in the art, without having to pay creative labor, can also be according to these
Accompanying drawing obtains other accompanying drawings.
Fig. 1 is that the reservoir based on background technology adjusts process and analysis result figure with crack attribute thresholds;
Fig. 2 is the flow chart according to the fractured-vuggy reservoir connectedness quantitative analysis method of the embodiment of the present invention;
Fig. 3 is the prediction result diagrammatic cross-section of different method for predicting reservoir;
Fig. 4 is the schematic diagram of carbonate rock fractured cave body development models;
Fig. 5 is the effective reservoir quantitative analysis distribution schematic diagram according to the embodiment of the present invention;
Fig. 6 is the extraction result schematic diagram according to the primary earthquake FRACTURE PREDICTION attribute volume of the embodiment of the present invention;
Fig. 7 A and Fig. 7 B are the connection pattern schematic diagram between the fracture and cave reservoir according to the embodiment of the present invention;
Fig. 8 is the contact relation schematic diagram between the fracture and cave reservoir according to the embodiment of the present invention;
Fig. 9 is the fracture-cavity units division result schematic diagram according to the embodiment of the present invention;
Figure 10 A to Figure 10 D are the crack attribute thresholds scope division result schematic diagram according to the embodiment of the present invention;
Figure 11 is the analysis result schematic diagram of Seismic prediction of fracture attribute volume;
Figure 12 is the structural representation according to the fractured-vuggy reservoir connectedness quantitative analysis device of the embodiment of the present invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other under the premise of creative work is not made
Embodiment, belong to the scope of protection of the invention.
The embodiment of the present invention provides a kind of fractured-vuggy reservoir connectedness quantitative analysis method and device.Below in conjunction with accompanying drawing pair
The present invention is described in detail.
The embodiment of the present invention provides a kind of fractured-vuggy reservoir connectedness quantitative analysis method, as shown in Fig. 2 the analysis method
Including:
Step S101:Inverting is carried out using earthquake and log data, obtains p-wave impedance body, according to well logging porosity with indulging
Relation between wave impedance, p-wave impedance body is converted into porosity cube;
Step S102:Using the lower limit of effective reservoir porosity in the porosity cube obtained by well log interpretation as critical
Value, determines effective reservoir;
Step S103:The connectedness of effective reservoir is judged according to the contact relation between effective reservoir, determines each fracture hole
The border of Reservoir Body;
Step S104:Using relevant or curvature algorithm, primary earthquake FRACTURE PREDICTION attribute volume is extracted;
Step S105:The connectedness of reservoir between well is determined by dynamic connectivity analysis method, and according to fracture and cave reservoir
Border judges the connected relation of fracture and cave reservoir between any two;
Step S106:Development degree of micro cracks in oil is determined according to connected relation, obtains splitting for primary earthquake FRACTURE PREDICTION attribute volume
The scope of attribute thresholds is stitched, and determines Seismic prediction of fracture attribute volume;
Step S107:According to the contact relation and crack attribute thresholds of fracture and cave reservoir and Seismic prediction of fracture attribute volume
Scope, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results.
It is connective by effective reservoir itself on the basis of reservoir quantitative forecast by above-mentioned steps S101~S107
Analysis, determine fracture and cave reservoir border;Crack is divided into two parts again, a part is the crack inside fracture and cave reservoir,
Another part is the crack between fracture and cave reservoir, using dynamic connectivity analysis result, is connected for the crack between fracture hole body
The general character is specialized in, and determines crack attribute thresholds, so as to quantitative forecast crack, further according to contact relation therebetween,
Judge reservoir communication.By the fractured-vuggy reservoir connectedness quantitative analysis method of the embodiment of the present invention, can reduce reservoir with
The multi-solution of fractue spacing analysis result, so as to improve the degree of accuracy of reservoir communication analysis.
Carbonate rock fractured cave type reservoir prediction at present mainly has two methods to be:Reservoir prediction side based on amplitude attribute
Method, the method for predicting reservoir based on inverting p-wave impedance.Research shows that p-wave impedance can more accurately reflect compared with amplitude attribute
(as shown in figure 3, Fig. 3 (a) is seismic amplitude profile, Fig. 3 (b) is frequency dividing amplitude section for the locus of fracture hole body and form
Figure, Fig. 3 (c) is reflected intensity paving figure, and Fig. 3 (d) is inverting wave impedance section figure), and have definitely with physical properties of rock
Correlation, so effective reservoir quantitative forecast precision is higher.
For fracture and vug carbonate reservoir, its development models as shown in figure 4, because seldom well drills fracture hole body,
Conventional determinated back analysis, now can only establish initial model with a constant, it is difficult to an accurate low frequency trend is obtained,
Therefore inverting impedance precision is influenced whether.And geostatistical inversion can combine well more preferably with earthquake, make full use of random
Modeling technique integrates the data of different scale, and inversion accuracy is higher than determinated back analysis.In embodiments of the present invention, above-mentioned steps
S101, it is that wave impedance inversion is carried out by Inversion of geostatistics, so as to obtain p-wave impedance body.Further, root
According to the relation between well log interpretation porosity and wave impedance, p-wave impedance body is converted into porosity cube.
Above-mentioned steps S102, using the lower limit of effective reservoir porosity in the porosity cube obtained by well log interpretation as
Critical value, determine effective reservoir.Invalid reservoir is proposed by the critical value, so as to obtain effective reservoir, identified effectively storage
The result of layer is as shown in Figure 5.Wherein, the relation between above-mentioned well logging porosity and p-wave impedance is:With the increasing of porosity
Greatly, the linear relationship of the linear relationship that p-wave impedance accordingly reduces, e.g. y=ax+b, wherein y are p-wave impedance, and x is hole
Degree.But coefficient a, b numerical value may be different and otherwise varied according to the section for carrying out earthquake detection.
Above-mentioned steps S103:The connectedness of effective reservoir is judged according to the contact relation between effective reservoir, is determined each
The border of fracture and cave reservoir, is specifically included:Judge to whether there is contact relation between effective reservoir, contact relation includes:Point connects
Touch, linear contact lay and face contact;It is connection between judgement effective reservoir if contact relation between effective reservoir be present;If
Contact relation is not present between effective reservoir, then is not connect between judgement effective reservoir.
After the border of each fracture and cave reservoir is determined, above-mentioned steps S104 is performed, using relevant or curvature algorithm, is carried
Primary earthquake FRACTURE PREDICTION attribute volume is taken, as shown in Figure 6.But due to FRACTURE PREDICTION attribute thresholds now determination not according to
According to random stronger, institute's Prediction of fracture degree is not accurate enough.Therefore, the threshold value of fracture prediction attribute is also needed to enter to advance
The quantitative analysis of one step.
Specifically, the further quantitative analysis of threshold value of fracture prediction attribute is realized by above-mentioned steps S105, is led to
The connectedness that dynamic connectivity analysis method determines reservoir between well is crossed, and fracture and cave reservoir two is judged according to fracture and cave reservoir border
Connected relation between two.By dynamic connectivity analysis method, reservoir communication between the well of fixed well is determined, and combine step
The fracture and cave reservoir border determined in S103, it is not connect or only connected by crack to judge which fracture and cave reservoir, as adjustment
The foundation of crack attribute thresholds.Wherein, as shown in figs. 7 a and 7b, Fig. 7 A are illustrated that the connection pattern between fracture and cave reservoir
The connection pattern of wall scroll slit, Fig. 7 B are illustrated that the connection pattern of a plurality of short seam.
Above-mentioned steps S106:Development degree of micro cracks in oil is determined according to connected relation, obtains primary earthquake FRACTURE PREDICTION attribute volume
Crack attribute thresholds scope, and determine Seismic prediction of fracture attribute volume.Specifically, step S106 includes following steps:
When connected relation is connects, with a constant (being, for example, 1) for unit, gradually reduce the number of crack attribute thresholds from big to small
Value, until occur crack between two fracture and cave reservoirs, and it is a fracture-cavity units to connect, and crack attribute thresholds now are made
For the higher limit of crack attribute thresholds;When connected relation is not connect, still in units of the constant, gradually increase from small to large
The numerical value of crack attribute thresholds, until fracture and cave reservoir separates, using crack attribute thresholds now as crack attribute thresholds
Lower limit;The scope of crack attribute thresholds is obtained according to the higher limit of crack attribute thresholds and lower limit.Wherein, above-mentioned connection
Relation determines according to the contact relation between two fracture and cave reservoirs, contact relation such as Fig. 8 institutes between two fracture and cave reservoirs
Show, mainly include:Point contact, linear contact lay and face contact, it is connection to exist between two fracture and cave reservoirs of contact relation, connection
Fracture hole body can form an above-mentioned fracture-cavity units, and analysis result as shown in Figure 9 can be ultimately formed, wherein, it is same
Color represents multiple fracture and cave reservoirs of one fracture-cavity units of composition.
It is to adjust the mistake of crack attribute thresholds by dynamic connectivity analysis result to constrain as shown in Figure 10 A to Figure 10 D
Journey:Confirmed by dynamic connectivity analysis result:This four mouthfuls of well connections of R301, R3013, R3011C, R3014X, belong to same seam
Hole unit, R3017 wells do not connect with above-mentioned four mouthfuls of wells, belong to another fracture-cavity units.Wherein, Figure 10 A show crack attribute
Threshold value takes analysis result when 235, and Figure 10 B show analysis result when crack attribute thresholds take 236, and Figure 10 C, which are shown, to be split
Seam attribute thresholds take analysis result when 254, and Figure 10 D show analysis result when crack attribute thresholds take 253, thus may be used
See, when crack attribute thresholds scope is 236-253, static connectivity analysis result is upper consistent with dynamic.
And it is possible to the crack attribute thresholds scope determined according to primary earthquake FRACTURE PREDICTION attribute volume, foundation, further
Ground determines Seismic prediction of fracture attribute volume, is that identified Seismic Fracture is pre- when crack attribute thresholds are 236 as shown in figure 11
Survey the result schematic diagram of attribute volume.
Above-mentioned steps S107 belongs to according to the contact relation between effective reservoir and Seismic prediction of fracture attribute volume, and crack
Property threshold value scope, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results, including:To effectively it store up
Layer merges one fused data body of generation with Seismic prediction of fracture attribute volume, according to the contact relation inside fused data body, judges
The connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results.
Compared to the analysis method in background technology, the fractured-vuggy reservoir connectedness quantitative analysis method of the embodiment of the present invention
The advantages of be:
1) analyzed under Geological Mode guidance, greatly reduce multi-solution;
2) use inversion method, more accurately determine reservoir position, for it is accurate portray fracture and cave reservoir border lay it is solid
Basis;
3) fracture connectedness carries out being directed to Journal of Sex Research, and more rationally determination crack attribute thresholds, precision of prediction are higher.
The embodiment of the present invention also provides a kind of fractured-vuggy reservoir connectedness quantitative analysis device, as shown in figure 12, the analysis
Device includes:Porosity cube generation unit 1, effective reservoir determining unit 2, fracture and cave reservoir border determining unit 3, primary earthquake
FRACTURE PREDICTION attribute volume extraction unit 4, connected relation judging unit 5, threshold range acquiring unit 6 and analysis result generation unit
7。
Carbonate rock fractured cave type reservoir prediction at present mainly has two methods to be:Reservoir prediction side based on amplitude attribute
Method, the method for predicting reservoir based on inverting p-wave impedance.Research shows that p-wave impedance can more accurately reflect compared with amplitude attribute
(as shown in figure 3, Fig. 3 (a) is seismic amplitude profile, Fig. 3 (b) is frequency dividing amplitude section for the locus of fracture hole body and form
Figure, Fig. 3 (c) is reflected intensity paving figure, and Fig. 3 (d) is inverting wave impedance section figure), and have definitely with physical properties of rock
Correlation, so effective reservoir quantitative forecast precision is higher.
For fracture and vug carbonate reservoir, its development models as shown in figure 4, because seldom well drills fracture hole body,
Conventional determinated back analysis, now can only establish initial model with a constant, it is difficult to an accurate low frequency trend is obtained,
Therefore inverting impedance precision is influenced whether.And geostatistical inversion can combine well more preferably with earthquake, make full use of random
Modeling technique integrates the data of different scale, and inversion accuracy is higher than determinated back analysis.In embodiments of the present invention, above-mentioned hole
Porosity body generation unit 1, it is that wave impedance inversion is carried out by Inversion of geostatistics, so as to obtain p-wave impedance body.
Further, according to the relation between well log interpretation porosity and wave impedance, p-wave impedance body is converted into porosity cube.
Effective reservoir porosity in the porosity cube that above-mentioned effective reservoir determining unit 2 is used to obtain by well log interpretation
Lower limit as critical value, determine effective reservoir.Invalid reservoir, so as to obtain effective reservoir, institute are proposed by the critical value
The result schematic diagram of the effective reservoir of determination is as shown in Figure 5.Wherein, the relation between above-mentioned well logging porosity and p-wave impedance
For:With the increase of porosity, linear relationship that p-wave impedance accordingly reduces, e.g. y=ax+b linear relationship, wherein y
For p-wave impedance, x is porosity.But coefficient a, b numerical value may be different and otherwise varied according to the section for carrying out earthquake detection.
Fracture and cave reservoir border determining unit 3 is used for the company that effective reservoir is judged according to the contact relation between effective reservoir
The general character, the border of each fracture and cave reservoir is determined, is specifically included:Judge to whether there is contact relation between effective reservoir, contact
Relation includes:Point contact, linear contact lay and face contact;If contact relation between effective reservoir be present, judge effective reservoir it
Between for connection;It is not connect between judgement effective reservoir if contact relation is not present between effective reservoir.
After the border of each fracture and cave reservoir is determined, pass through the profit of primary earthquake FRACTURE PREDICTION attribute volume extraction unit 4
With relevant or curvature algorithm, primary earthquake FRACTURE PREDICTION attribute volume is extracted, as shown in Figure 6.But due to FRACTURE PREDICTION category now
Property threshold value determination have no basis, random stronger, institute's Prediction of fracture degree is not accurate enough.Therefore, also need fracture pre-
The threshold value for surveying attribute carries out further quantitative analysis.
Specifically, the further quantitative analysis of threshold value of fracture prediction attribute is by above-mentioned connected relation judging unit
5 realize, the connectedness of reservoir between well are determined by dynamic connectivity analysis method, and judge to stitch according to fracture and cave reservoir border
The connected relation of hole Reservoir Body between any two.By dynamic connectivity analysis method, reservoir communication between the well of fixed well is determined,
And the fracture and cave reservoir border determined in step S103 is combined, it is not connect or only connected by crack to judge which fracture and cave reservoir
It is logical, the foundation as adjustment crack attribute thresholds.Wherein, the connection pattern between fracture and cave reservoir as shown in figs. 7 a and 7b,
Fig. 7 A are illustrated that the connection pattern of wall scroll slit, and Fig. 7 B are illustrated that the connection pattern of a plurality of short seam.
Threshold range acquiring unit 6 is used to determine development degree of micro cracks in oil according to connected relation, and it is pre- to obtain primary earthquake crack
The scope of the crack attribute thresholds of attribute volume is surveyed, and determines Seismic prediction of fracture attribute volume.The threshold range acquiring unit 6 is specific
Perform following steps:When connected relation is connects, with a constant (being, for example, 1) for unit, gradually reduce split from big to small
The numerical value of attribute thresholds is stitched, until occur crack between two fracture and cave reservoirs, and it is a fracture-cavity units to connect, by now
Higher limit of the crack attribute thresholds as crack attribute thresholds;When connected relation is not connect, still in units of the constant, from
It is small to the big numerical value for gradually increasing crack attribute thresholds, until fracture and cave reservoir separates, using crack attribute thresholds now as
The lower limit of crack attribute thresholds;The scope of crack attribute thresholds is obtained according to the higher limit of crack attribute thresholds and lower limit.
Wherein, above-mentioned connected relation determines according to the contact relation between two fracture and cave reservoirs, between two fracture and cave reservoirs
Contact relation is as shown in figure 8, mainly include:Point contact, linear contact lay and face contact, exist contact relation two fracture and cave reservoirs it
Between be connection, the fracture hole body of connection can form an above-mentioned fracture-cavity units, and can ultimately form analysis knot as shown in Figure 9
Fruit, wherein, same color represents multiple fracture and cave reservoirs of one fracture-cavity units of composition.
It is to adjust the mistake of crack attribute thresholds by dynamic connectivity analysis result to constrain as shown in Figure 10 A to Figure 10 D
Journey:Confirmed by dynamic connectivity analysis result:This four mouthfuls of well connections of R301, R3013, R3011C, R3014X, belong to same seam
Hole unit, R3017 wells do not connect with above-mentioned four mouthfuls of wells, belong to another fracture-cavity units.Wherein, Figure 10 A show crack attribute
Threshold value takes analysis result when 235, and Figure 10 B show analysis result when crack attribute thresholds take 236, and Figure 10 C, which are shown, to be split
Seam attribute thresholds take analysis result when 254, and Figure 10 D show analysis result when crack attribute thresholds take 253, thus may be used
See, when crack attribute thresholds scope is 236-253, static state connectivity analysis result and dynamic are upper consistent.
Also, the threshold range acquiring unit 6 can be according to primary earthquake FRACTURE PREDICTION attribute volume, according to the crack determined
Attribute thresholds scope, Seismic prediction of fracture attribute volume is further determined, be when crack attribute thresholds are 236 as shown in figure 11
When, the result schematic diagram of identified Seismic prediction of fracture attribute volume.
Analysis result generation unit 7 is used for according to the contact relation between effective reservoir and Seismic prediction of fracture attribute volume,
And the scope of crack attribute thresholds, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results, bag
Include:Effective reservoir is merged to one fused data body of generation with Seismic prediction of fracture attribute volume, according to connecing inside fused data body
The relation of touching, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results.
Compared to the analysis method in background technology, the fractured-vuggy reservoir connectedness quantitative analysis device of the embodiment of the present invention
The advantages of analysis process realized, is:
1) analyzed under Geological Mode guidance, greatly reduce multi-solution;
2) use inversion method, more accurately determine reservoir position, for it is accurate portray fracture and cave reservoir border lay it is solid
Basis;
3) fracture connectedness carries out being directed to Journal of Sex Research, and more rationally determination crack attribute thresholds, precision of prediction are higher.
One of ordinary skill in the art will appreciate that realize that all or part of step in above-described embodiment method can lead to
Program is crossed to instruct the hardware of correlation to complete, the program can be stored in a computer read/write memory medium, such as
ROM/RAM, magnetic disc, CD etc..
Particular embodiments described above, the purpose of the present invention, technical scheme and beneficial effect are carried out further in detail
Describe in detail it is bright, should be understood that the foregoing is only the present invention specific embodiment, the guarantor being not intended to limit the present invention
Scope is protected, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., should be included in this
Within the protection domain of invention.
Claims (8)
1. a kind of fractured-vuggy reservoir connectedness quantitative analysis method, it is characterised in that methods described includes:
Inverting is carried out using earthquake and log data, p-wave impedance body is obtained, according between well logging porosity and p-wave impedance
Relation, the p-wave impedance body is converted into porosity cube;
Using the lower limit of effective reservoir porosity in the porosity cube obtained by well log interpretation as critical value, it is determined that having
Imitate reservoir;
The connectedness of the effective reservoir is judged according to the contact relation between the effective reservoir, determines each fracture and cave reservoir
Border;
Using relevant or curvature algorithm, primary earthquake FRACTURE PREDICTION attribute volume is extracted;
The connectedness of reservoir between well is determined by dynamic connectivity analysis method, and is judged according to described fracture and cave reservoir border
The connected relation of the fracture and cave reservoir between any two;
Development degree of micro cracks in oil is determined according to the connected relation, obtains the crack attribute of the primary earthquake FRACTURE PREDICTION attribute volume
The scope of threshold value, and determine Seismic prediction of fracture attribute volume;
According to the contact relation and the crack attribute thresholds of the fracture and cave reservoir and the Seismic prediction of fracture attribute volume
Scope, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results;
Development degree of micro cracks in oil is determined according to the connected relation, obtains the crack attribute of the primary earthquake FRACTURE PREDICTION attribute volume
The scope of threshold value, including:
When the connected relation is connects, in units of a constant, gradually reduce the crack attribute thresholds from big to small
Numerical value, until occur crack between two fracture and cave reservoirs, and it is a fracture-cavity units to connect, by the crack now
Higher limit of the attribute thresholds as the crack attribute thresholds;
When the connected relation is not connect, in units of the constant, gradually increase the crack attribute threshold from small to large
The numerical value of value, until the fracture and cave reservoir separates, using the crack attribute thresholds now as the crack attribute thresholds
Lower limit;
The scope of the crack attribute thresholds is obtained according to the higher limit of the crack attribute thresholds and lower limit.
2. fractured-vuggy reservoir connectedness quantitative analysis method according to claim 1, it is characterised in that the well logging hole
Spend p-wave impedance between relation be:With the increase of porosity, linear relationship that p-wave impedance accordingly reduces.
3. fractured-vuggy reservoir connectedness quantitative analysis method according to claim 1, it is characterised in that according to described effective
Contact relation between reservoir judges the connectedness of the effective reservoir, determines Reservoir Body border, including:
Judge with the presence or absence of contact, the contact to include between the effective reservoir:Point contact, linear contact lay and face contact;
If contact between the effective reservoir be present, judge between the effective reservoir for connection;Otherwise, it is determined that described have
It is not connect between effect reservoir.
4. fractured-vuggy reservoir connectedness quantitative analysis method according to claim 1, it is characterised in that according to the fracture hole
Reservoir Body and the contact relation of the Seismic prediction of fracture attribute volume and the scope of the crack attribute thresholds, judge that fracture hole type stores up
The connectedness of layer, reservoir communication quantitative analysis results are generated, including:
The effective reservoir is merged to one fused data body of generation with the Seismic prediction of fracture attribute volume, according to the fusion number
According to the contact relation in internal portion, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results.
5. a kind of fractured-vuggy reservoir connectedness quantitative analysis device, it is characterised in that described device includes:
Porosity cube generation unit, for carrying out inverting using earthquake and log data, p-wave impedance body is obtained, according to well logging hole
Relation between porosity and p-wave impedance, the p-wave impedance body is converted into porosity cube;
Effective reservoir determining unit, for under effective reservoir porosity in the porosity cube that is obtained by well log interpretation
Limit value determines effective reservoir as critical value;
Fracture and cave reservoir border determining unit, for judging the effective reservoir according to the contact relation between the effective reservoir
Connectedness, determine the fracture and cave reservoir border of each Reservoir Body;
Primary earthquake FRACTURE PREDICTION attribute volume extraction unit, for pre- using relevant or curvature algorithm, extraction primary earthquake crack
Survey attribute volume;
Connected relation judging unit, for determining the connectedness of reservoir between well by dynamic connectivity analysis method, and according to institute
The fracture and cave reservoir border stated judges the connected relation of the fracture and cave reservoir between any two;
Threshold range acquiring unit, for determining development degree of micro cracks in oil according to the connected relation, initially shattered described in acquisition
The scope of the crack attribute thresholds of seam prediction attribute volume, and determine Seismic prediction of fracture attribute volume;
Analysis result generation unit, for the contact relation according to the fracture and cave reservoir and the Seismic prediction of fracture attribute volume
And the scope of the crack attribute thresholds, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results;
The threshold range acquiring unit is specifically used for:
When the connected relation is connects, in units of a constant, gradually reduce the crack attribute thresholds from big to small
Numerical value, until occur crack between two fracture and cave reservoirs, and it is a fracture-cavity units to connect, by the crack now
Higher limit of the attribute thresholds as the crack attribute thresholds;
When the connected relation is not connect, in units of the constant, gradually increase the crack attribute threshold from small to large
The numerical value of value, until the fracture and cave reservoir separates, using the crack attribute thresholds now as the crack attribute thresholds
Lower limit;
The scope of the crack attribute thresholds is obtained according to the higher limit of the crack attribute thresholds and lower limit.
6. fractured-vuggy reservoir connectedness quantitative analysis device according to claim 5, it is characterised in that the well logging hole
Spend p-wave impedance between relation be:With the increase of porosity, linear relationship that p-wave impedance accordingly reduces.
7. fractured-vuggy reservoir connectedness quantitative analysis device according to claim 5, it is characterised in that the fracture hole preserves
Body border determining unit is specifically used for:
Judge with the presence or absence of contact, the contact to include between the effective reservoir:Point contact, linear contact lay and face contact;
If contact between the effective reservoir be present, judge between the effective reservoir for connection;Otherwise, it is determined that described have
It is not connect between effect reservoir.
8. fractured-vuggy reservoir connectedness quantitative analysis device according to claim 5, it is characterised in that the analysis result
Generation unit is specifically used for:
The effective reservoir is merged to one fused data body of generation with the Seismic prediction of fracture attribute volume, according to the fusion number
According to the contact relation in internal portion, judge the connectedness of fractured-vuggy reservoir, generate reservoir communication quantitative analysis results.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510187396.XA CN104818982B (en) | 2015-04-20 | 2015-04-20 | Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510187396.XA CN104818982B (en) | 2015-04-20 | 2015-04-20 | Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104818982A CN104818982A (en) | 2015-08-05 |
CN104818982B true CN104818982B (en) | 2017-12-05 |
Family
ID=53729434
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510187396.XA Active CN104818982B (en) | 2015-04-20 | 2015-04-20 | Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104818982B (en) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106703779A (en) * | 2016-12-02 | 2017-05-24 | 中国石油化工股份有限公司 | Injection-production well pattern construction method suitable for carbonate fractured-cave reservoirs |
CN106837297B (en) * | 2016-12-22 | 2020-04-10 | 中国石油天然气股份有限公司 | Method for identifying connectivity among wells and predicting oil-water dynamic state |
CN107045145B (en) * | 2016-12-30 | 2019-02-15 | 中国石油集团东方地球物理勘探有限责任公司 | Indication using prestack seismic amplitude under seismic sequence control changes detection fracture hole method with offset distance |
CN108572392B (en) * | 2017-03-09 | 2020-07-28 | 中国石油化工股份有限公司 | Method for depicting carbonate rock fracture-cave communicating body |
CN107143312B (en) * | 2017-06-20 | 2019-08-30 | 中国石油化工股份有限公司 | Fracture and vug carbonate reservoir energy self-balancing development approach |
CN107083939A (en) * | 2017-06-26 | 2017-08-22 | 中国石油化工股份有限公司 | Three-dimensional development method for the disconnected solution oil reservoir of carbonate rock |
CN107191175A (en) * | 2017-07-11 | 2017-09-22 | 中国石油化工股份有限公司 | Flooding pattern construction method for the disconnected solution oil reservoir of carbonate rock |
CN109116428B (en) * | 2018-07-02 | 2020-06-09 | 中国石油天然气股份有限公司 | Fracture-cavity carbonate reservoir uncertainty modeling method and device |
CN111173495B (en) * | 2018-10-23 | 2023-09-26 | 中国石油化工股份有限公司 | Fracture-cavity type oil reservoir fracture-cavity unit communication condition evaluation method |
CN109447532B (en) * | 2018-12-28 | 2021-10-22 | 中国石油大学(华东) | Oil reservoir inter-well connectivity determination method based on data driving |
CN109577959B (en) * | 2019-01-23 | 2022-03-11 | 四川富利斯达石油科技发展有限公司 | Method for measuring crack connectivity of adjacent fracturing sections by using tracer |
CN111830558B (en) * | 2019-04-17 | 2022-07-22 | 中国石油化工股份有限公司 | Fracture zone engraving method |
CN110259426B (en) * | 2019-07-02 | 2021-12-07 | 捷贝通石油技术集团股份有限公司 | Method for evaluating pressure channeling degree between unconventional platform wells |
CN113495295B (en) * | 2020-04-07 | 2024-01-30 | 中国石油天然气股份有限公司 | Weak-connectivity fracture-cavity reservoir identification method, device, equipment and storage medium |
CN111596351B (en) * | 2020-04-28 | 2023-04-25 | 中国石油天然气股份有限公司 | Carbonate rock conductor system quantitative evaluation method, system, device and storage medium |
CN112859171B (en) * | 2020-11-19 | 2024-01-30 | 中国石油天然气股份有限公司 | Method and apparatus for reservoir connectivity analysis using effective fractures |
CN113156499B (en) * | 2021-02-26 | 2023-05-26 | 中国石油天然气股份有限公司 | Seismic data post-stack quantitative prediction method for fractured reservoir in basin area |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102373923A (en) * | 2010-08-20 | 2012-03-14 | 中国石油天然气股份有限公司 | Reservoir stratum identification method |
CN102759757A (en) * | 2012-07-06 | 2012-10-31 | 中国石油天然气股份有限公司 | Diabase alteration zone reservoir prediction method based on well-seismic combination |
CN103413030A (en) * | 2013-07-24 | 2013-11-27 | 中国石油天然气股份有限公司 | Fracture-cavity type carbonate rock gas reservoir dynamic analysis method and system |
CN103869359A (en) * | 2014-02-25 | 2014-06-18 | 中国石油天然气股份有限公司 | Method and device for predicting cracks through seismic longitudinal wave multi-azimuth attribute elliptic fitting |
-
2015
- 2015-04-20 CN CN201510187396.XA patent/CN104818982B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102373923A (en) * | 2010-08-20 | 2012-03-14 | 中国石油天然气股份有限公司 | Reservoir stratum identification method |
CN102759757A (en) * | 2012-07-06 | 2012-10-31 | 中国石油天然气股份有限公司 | Diabase alteration zone reservoir prediction method based on well-seismic combination |
CN103413030A (en) * | 2013-07-24 | 2013-11-27 | 中国石油天然气股份有限公司 | Fracture-cavity type carbonate rock gas reservoir dynamic analysis method and system |
CN103869359A (en) * | 2014-02-25 | 2014-06-18 | 中国石油天然气股份有限公司 | Method and device for predicting cracks through seismic longitudinal wave multi-azimuth attribute elliptic fitting |
Non-Patent Citations (2)
Title |
---|
塔里木盆地奥陶系碳酸盐岩大型缝洞集合体定量描述;杨鹏飞等;《岩性油气藏》;20131231;第25卷(第6期);第89-94页 * |
碳酸盐岩油气藏缝洞单元多学科综合划分方法研究——以塔里木盆地轮古11 井区为例;邓兴梁等;《海相油气地质》;20130430;第18卷(第2期);第72-78页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104818982A (en) | 2015-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104818982B (en) | Method and device for quantitatively analyzing connectivity of fracture-cavity reservoir | |
US10310136B2 (en) | Lateral placement and completion design for improved well performance of unconventional reservoirs | |
US7756694B2 (en) | Method for interactive automation of fault modeling including a method for intelligently sensing fault-fault relationships | |
CN104977617B (en) | Reservoir Fracture recognition methods and imaging logging Reservoir Fracture recognition methods | |
CN105507890B (en) | A kind of Logging Evaluation Method suitable for shale gas GEOLOGICAL MODELS OF PETROLEUM RESERVOIR | |
CN106951660A (en) | Sea facies clastic rock horizontal well reservoir logging interpretation method and device | |
MX2015002141A (en) | System and method for performing stimulation operations. | |
CN102749647B (en) | Quantitative method for space of holes in reservoir | |
CN108397130B (en) | Boring method and device | |
CN108562938B (en) | Method, device and system for eliminating frequency dispersion effect | |
CN107292074A (en) | A kind of inter well connectivity determination methods | |
CN115166853A (en) | Shale gas reservoir natural fracture model establishing method and device, electronic equipment and medium | |
CN109425899A (en) | Method and device for predicting distribution of fractured zone of carbonate rock | |
KR20200027264A (en) | A method for selecting a sweet spot in shale gas | |
CN107229076B (en) | A method of temperature-responsive signature analysis is carried out based on well-log information | |
CN107861162B (en) | Microelectrode logging data-based natural crack identification method and system | |
US11073463B2 (en) | Methods for modeling permeability in layered rock formations | |
CN103282797B (en) | The method and system of computed tomography drop | |
CN115032062A (en) | Method and device for determining compressibility of shale oil reservoir and compressibility determining system | |
CN110208860A (en) | A kind of prediction technique and device of Igneous rock invasion range | |
Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
CN113325468B (en) | Reservoir fracture-cavity distribution range prediction method and device | |
CN110208861B (en) | Prediction method and device for constructing soft coal development area | |
CN109283574A (en) | Low frequency model construction method and computer readable storage medium | |
CN110308488B (en) | Method and system for determining cave filling degree |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |