CN104834007A - Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process - Google Patents
Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process Download PDFInfo
- Publication number
- CN104834007A CN104834007A CN201510219252.8A CN201510219252A CN104834007A CN 104834007 A CN104834007 A CN 104834007A CN 201510219252 A CN201510219252 A CN 201510219252A CN 104834007 A CN104834007 A CN 104834007A
- Authority
- CN
- China
- Prior art keywords
- seismic
- disparity
- stuff
- filling operation
- reservoir
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 30
- 239000011435 rock Substances 0.000 claims abstract description 28
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 238000011160 research Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000000284 extract Substances 0.000 claims description 9
- 230000001186 cumulative effect Effects 0.000 claims description 7
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 abstract description 5
- 239000004576 sand Substances 0.000 abstract 1
- 239000011148 porous material Substances 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for calculating the filling degree of a carbonate fracture-cave reservoir in a seismic inversion process, which comprises the following steps: extracting a seismic inversion wave impedance plane graph along a seismic interpretation horizon from a research area; obtaining filling materials and filling degree information in the carbonate rock fracture-cave type reservoir; selecting sample points with similar filling materials but different filling degrees and positioning geodetic coordinates of the sample points; determining effective sample points and calculating seismic wave impedance difference values among the effective sample points; performing intersection analysis on the sand filling quantity difference value and the seismic wave impedance difference value between the effective sample points to obtain a functional relation; calculating the filling degree of the fracture-cavity reservoir at each effective sample point, and obtaining an operational relational expression of the seismic wave impedance difference value and the filling degree of the sandy filling material according to the functional relational expression and the filling degree; and selecting a seismic reflection area with beaded reflection on the seismic section, and calculating the filling degree of the fracture-cavity reservoir in the seismic reflection area by using an operational relation.
Description
Technical field
The present invention relates to technical field of physical geography, particularly relate to a kind of method calculating carbonate rock fractured cave type reservoir filling operation in seismic inversion process.
Background technology
Seismic inversion be make full use of drilling well, the abundant information of the aspect such as relevant configuration, lithology, layer position that well logging, geologic information provide, carry out by common seismic section the information such as density, wave impedance, speed, permeability, factor of porosity, pressure, Sandstone Percentage that subsurface formations is derived in inverting.Wave impedance inversion refers to the earthquake special processing interpretation technique utilizing seismic data inverting to derive formation wave impedance (or speed).Compared with seismic model determination predicting oil/gas, neural network prediction formation parameter, amplitude matching Reservoir Prediction thickness etc. statistically method, wave impedance inversion processing interpretation technique has very clear and definite physical significance, it is the Deterministic Methods of reservoir lithology prediction, reservoir characterization, achieve significant geological effect in actual applications, therefore seismic inversion refers in particular to wave impedance inversion usually.
Current carbonate rock fractured cave type reservoir reserves evaluation work mainly uses this parameter of average pore, there is larger error in average pore method, cause certain adverse effect to the exploratory development production work of reality in the process calculating seam hole type hydrocarbon-bearing pool.For the problems referred to above, at present effective solution is not yet proposed.
Summary of the invention
The invention provides a kind of method calculating carbonate rock fractured cave type reservoir filling operation in seismic inversion process, use average pore method at least to solve current carbonate rock fractured cave type reservoir reserves evaluation work, there is the problem compared with big error.
According to an aspect of the present invention, provide a kind of method calculating carbonate rock fractured cave type reservoir filling operation in seismic inversion process, comprise: step 1, the geological data in research on utilization district, log data and logging data complete seismic interpretation and seismic inversion to described study area, extract seismic inversion wave impedance planimetric map along seismic interpretation layer position on inversion result; Step 2, the detail correlation of reservoir bed is carried out on the basis of described study area being carried out to core observation, logging explanation report analysis, when waiting in screen work in conjunction with actual production data, obtain the information of the inner stuff of carbonate rock fractured cave type reservoir and filling operation; Observe and data statistics according to core hole, select stuff similar but the sample spot that filling operation is different; And the terrestrial coordinate of described sample spot is jointly located according to well location coordinate, hole deviation data and layer bit depth data; Step 3, described seismic inversion wave impedance planimetric map being pressed numeric format derives, and obtains the seismic impedance numerical value along seismic interpretation layer position of described study area; Determine the sample spot identical with the terrestrial coordinate in described seismic impedance numerical value, as effective sample spot; And the seismic impedance values of disparity calculated between described effective sample spot; Step 4, according to research object accuracy requirement, setting fill-up minimum calculation unit, chiltern stuff fill-up values of disparity between described effective sample spot and described seismic impedance values of disparity are carried out cross analysis, obtains the functional relation of seismic impedance values of disparity and chiltern stuff fill-up values of disparity; Step 5, according to research object accuracy requirement, setting filling operation minimum calculation unit, the cumulative volume of each fractured-vuggy reservoir at each effective sample spot place is asked for according to described seismic impedance numerical value and seismic section, calculate the ratio of chiltern stuff fill-up and cumulative volume in each fractured-vuggy reservoir, obtain the filling operation of described fractured-vuggy reservoir; According to described seismic impedance values of disparity and the described functional relation of chiltern stuff fill-up values of disparity and the filling operation of described fractured-vuggy reservoir, obtain the operation relation formula of described seismic impedance values of disparity and chiltern stuff filling operation; Step 6, choose the seismic reflection district described seismic section with beading reflection, use the operation relation formula of described seismic impedance values of disparity and described chiltern stuff filling operation to carry out computing to the described seismic impedance numerical value of deriving along described seismic interpretation floor position in the described seismic reflection district chosen, obtain the filling operation of the fractured-vuggy reservoir in described seismic reflection district.
In one embodiment, described step 1 comprises: according to horizon calibration in described log data and described logging data completed well; Seismic interpretation is completed in conjunction with Seismic reflection character according to aboveground horizon calibration; The secondary computing of seismic inversion is carried out according to described geological data; Along seismic interpretation layer position, the inverting Wave Impedance Data Volume after secondary computing is cut into slices, obtain the described seismic inversion wave impedance planimetric map along described seismic interpretation layer position.
In one embodiment, the functional relation of described seismic impedance values of disparity and chiltern stuff fill-up values of disparity is Y=a × X+b, and wherein, Y is chiltern stuff fill-up values of disparity, X is seismic impedance values of disparity, and a, b are cross analysis gained constant.
In one embodiment, the operation relation formula of described seismic impedance values of disparity and chiltern stuff filling operation is Z=(c × X+d)/V, wherein, Z is chiltern stuff filling operations, X is seismic impedance values of disparity, c, d are cross analysis gained computing constant, and V is the volume of fractured-vuggy reservoir.
In one embodiment, in described step 2, in the different tectonic position in tectonic element and different tectonic element, extract described sample spot.
In one embodiment, in described step 4 and described step 5, described fill-up minimum calculation unit and the described filling operation minimum calculation unit of setting are all not less than research object precision prescribed.
By the method for calculating carbonate rock fractured cave type reservoir filling operation of the present invention, in seismic inversion process, quantize the inner filling operation of carbonate rock fractured cave type reservoir in seismic inversion wave impedance sign numerically, the objective numerical value embodying the inner true hole porosity of reservoir, avoid the error using average pore method to cause oil and gas reserves estimation, thus clearly can react the true geologic reserve of carbonate rock fractured cave type reservoir objectively.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, and form a application's part, schematic description and description of the present invention, for explaining the present invention, does not form limitation of the invention.In the accompanying drawings:
Fig. 1 is the process flow diagram calculating the method for carbonate rock fractured cave type reservoir filling operation in the seismic inversion process of the embodiment of the present invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on embodiments of the invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to protection scope of the present invention.
Embodiments provide a kind of method calculating carbonate rock fractured cave type reservoir filling operation (filling operation is also referred to as pack completeness) in seismic inversion process, Fig. 1 is the process flow diagram calculating the method for carbonate rock fractured cave type reservoir filling operation in the seismic inversion process of the embodiment of the present invention.As shown in Figure 1, the method comprises:
Step S101, the geological data in research on utilization district, log data and logging data complete seismic interpretation and seismic inversion to study area, extract seismic inversion wave impedance planimetric map along seismic interpretation layer position on inversion result.Wherein, geological data is the data of segy form.Extract seismic inversion wave impedance planimetric map along seismic interpretation layer position, require seismic horizon and fault interpretation, then synthetic record and time and depth transfer formula.
Step S102, the detail correlation of reservoir bed is carried out on the basis of study area being carried out to core observation, logging explanation report analysis, when waiting in screen work in conjunction with actual production data, obtain the information of the inner stuff of carbonate rock fractured cave type reservoir and filling operation; Observe and data statistics according to core hole, select stuff similar but the sample spot that filling operation is different; And according to the terrestrial coordinate of the common localizing sample point of well location coordinate, hole deviation data and layer bit depth data.Meanwhile, also can add up sample spot data, sample spot data comprise main stuff and filling operation.
Step S103, seismic inversion wave impedance planimetric map being pressed numeric format derives, and obtains the seismic impedance numerical value along seismic interpretation layer position of study area; Determine the sample spot identical with the terrestrial coordinate in seismic impedance numerical value, as effective sample spot; And the seismic impedance values of disparity calculated between effective sample spot.In this step, terrestrial coordinate in seismic impedance numerical value is carried out corresponding with the terrestrial coordinate of sample spot, requiring two coordinates as far as possible close to even overlapping, requiring that in plane, distance is the shortest, need consider that in sample spot coordinate, position angle and offset distance can cause the various factors of error.
Step S104, according to research object accuracy requirement, setting fill-up minimum calculation unit, chiltern stuff fill-up values of disparity between effective sample spot and seismic impedance values of disparity are carried out cross analysis, obtains the functional relation of seismic impedance values of disparity and chiltern stuff fill-up values of disparity.Wherein, in this functional relation, seismic impedance values of disparity increases with chiltern stuff fill-up and increases.
Step S105, according to research object accuracy requirement, setting filling operation minimum calculation unit, base area seismic wave impedance value and seismic section ask for the cumulative volume of each fractured-vuggy reservoir at each effective sample spot place, calculate the ratio of chiltern stuff fill-up and cumulative volume in each fractured-vuggy reservoir, obtain the filling operation of fractured-vuggy reservoir; The functional relation of base area seismic wave resistance difference numerical value and chiltern stuff fill-up values of disparity and the filling operation of fractured-vuggy reservoir, obtain the operation relation formula of seismic impedance values of disparity and chiltern stuff filling operation.
Step S106, choose seismic reflection district seismic section with beading reflection, use the operation relation formula of seismic impedance values of disparity and chiltern stuff filling operation to carry out computing to the seismic impedance numerical value of deriving along seismic interpretation floor position in the seismic reflection district chosen, obtain the filling operation of the fractured-vuggy reservoir in seismic reflection district.Computing in this step refers to that the operation relation formula utilizing step 5 gained carries out computing again to seismic event impedance value, and seismic impedance numerical value is converted to stuff numerical value.
Pass through said method, quantize the inner filling operation of carbonate rock fractured cave type reservoir in seismic inversion wave impedance sign numerically, the objective numerical value embodying the inner true hole porosity of reservoir, avoid the error using average pore method to cause oil and gas reserves estimation, thus clearly can react the true geologic reserve of carbonate rock fractured cave type reservoir objectively.Different main stuff and the seismic inversion Acoustic Impedance Data of pack completeness are discrepant, therefore, the present invention analyzes the geophysical response in the seam hole of different filling operation and utilizes above-mentioned difference to launch in the plane to release the filling operation without wellblock fractured-vuggy reservoir.
In one embodiment, step S101 can comprise: according to horizon calibration in log data and logging data completed well; Seismic interpretation is completed in conjunction with Seismic reflection character according to aboveground horizon calibration; The secondary computing of seismic inversion is carried out according to geological data; Along seismic interpretation layer position, the inverting Wave Impedance Data Volume after secondary computing is cut into slices, obtain the seismic inversion wave impedance planimetric map along seismic interpretation layer position.
In one embodiment, the functional relation of seismic impedance values of disparity and chiltern stuff fill-up values of disparity is Y=a × X+b, and wherein, Y is chiltern stuff fill-up values of disparity, X is seismic impedance values of disparity, and a, b are cross analysis gained constant.
In one embodiment, the operation relation formula of seismic impedance values of disparity and chiltern stuff filling operation is Z=(c × X+d)/V, wherein, Z is chiltern stuff filling operations, X is seismic impedance values of disparity, c, d are cross analysis gained computing constant, and V is the volume of fractured-vuggy reservoir.
In one embodiment, in step 2, in the different tectonic position in tectonic element and different tectonic element, extract sample spot.
In one embodiment, in step 4 and step 5, fill-up minimum calculation unit and the filling operation minimum calculation unit of setting are all not less than research object precision prescribed.
From the above, the method calculating carbonate rock fractured cave type reservoir filling operation in the seismic inversion process of the embodiment of the present invention comprises: utilize existing seismic volume segy data, log data and logging data to complete seismic interpretation and seismic inversion to study area, extract seismic inversion wave impedance planimetric map along seismic interpretation layer on inversion result wave resistance body; Bore according to practical condition statistics drilling well institute and meet filling operation in fractured-vuggy reservoir and main stuff; By x, y terrestrial coordinate to well point production data with extract seismic impedance coordinate and seek corresponding relation; Obtain seismic impedance numerical value to increase and the relational expression of increase with the total fill-up of chiltern stuff; The relational expression obtained is used to carry out secondary computing along on seismic interpretation layer bit plane wave impedance numerical value basis, filling operation in other carbonate rock fractured cave type reservoirs of study area is estimated, obtain fill-up true and quantitative in carbonate rock fractured cave type reservoir, and obtain reservoir filling operation in conjunction with cumulative volume, thus comparatively objectively can embody the inner true hole porosity of its reservoir.
In sum, the present invention by quantifying the inner filling operation of carbonate rock fractured cave type reservoir in seismic inversion wave impedance sign numerically, the objective numerical value embodying the inner true hole porosity of reservoir, avoid the error using average pore method to cause oil and gas reserves estimation, thus clearly can react the true geologic reserve of carbonate rock fractured cave type reservoir objectively.Calculating chiltern stuff in carbonate rock fractured cave type reservoir after filling operation, think that the net porosity in carbonate rock fractured cave type reservoir should be between the factor of porosity of stuff pack portion and 100% factor of porosity of non-pack portion, according to the filling operation of stuff determine weight coefficient and the practical porosity that distributes, instead of the average pore that uses of oil field at present.Therefore, utilize the method calculating carbonate rock fractured cave type reservoir filling operation in seismic inversion process, the numerical value of net porosity in carbonate rock fractured cave type reservoir can be embodied the most truly, for oil field Resources calculation provides foundation.
Describe and can be understood in process flow diagram or in this any process otherwise described or method, represent and comprise one or more for realizing the module of the code of the executable instruction of the step of specific logical function or process, fragment or part, and the scope of the preferred embodiment of the present invention comprises other realization, wherein can not according to order that is shown or that discuss, comprise according to involved function by the mode while of basic or by contrary order, carry out n-back test, this should understand by embodiments of the invention person of ordinary skill in the field.
In the description of this instructions, specific features, structure, material or feature that the description of reference term " embodiment ", " some embodiments ", " example ", " concrete example " or " some examples " etc. means to describe in conjunction with this embodiment or example are contained at least one embodiment of the present invention or example.In this manual, identical embodiment or example are not necessarily referred to the schematic representation of above-mentioned term.And the specific features of description, structure, material or feature can combine in an appropriate manner in any one or more embodiment or example.
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (6)
1. calculate a method for carbonate rock fractured cave type reservoir filling operation in seismic inversion process, it is characterized in that, comprising:
Step 1, the geological data in research on utilization district, log data and logging data complete seismic interpretation and seismic inversion to described study area, extract seismic inversion wave impedance planimetric map along seismic interpretation layer position on inversion result;
Step 2, the detail correlation of reservoir bed is carried out on the basis of described study area being carried out to core observation, logging explanation report analysis, when waiting in screen work in conjunction with actual production data, obtain the information of the inner stuff of carbonate rock fractured cave type reservoir and filling operation; Observe and data statistics according to core hole, select stuff similar but the sample spot that filling operation is different; And the terrestrial coordinate of described sample spot is jointly located according to well location coordinate, hole deviation data and layer bit depth data;
Step 3, described seismic inversion wave impedance planimetric map being pressed numeric format derives, and obtains the seismic impedance numerical value along seismic interpretation layer position of described study area; Determine the sample spot identical with the terrestrial coordinate in described seismic impedance numerical value, as effective sample spot; And the seismic impedance values of disparity calculated between described effective sample spot;
Step 4, according to research object accuracy requirement, setting fill-up minimum calculation unit, chiltern stuff fill-up values of disparity between described effective sample spot and described seismic impedance values of disparity are carried out cross analysis, obtains the functional relation of seismic impedance values of disparity and chiltern stuff fill-up values of disparity;
Step 5, according to research object accuracy requirement, setting filling operation minimum calculation unit, the cumulative volume of each fractured-vuggy reservoir at each effective sample spot place is asked for according to described seismic impedance numerical value and seismic section, calculate the ratio of chiltern stuff fill-up and cumulative volume in each fractured-vuggy reservoir, obtain the filling operation of described fractured-vuggy reservoir; According to described seismic impedance values of disparity and the described functional relation of chiltern stuff fill-up values of disparity and the filling operation of described fractured-vuggy reservoir, obtain the operation relation formula of described seismic impedance values of disparity and chiltern stuff filling operation;
Step 6, choose the seismic reflection district described seismic section with beading reflection, use the operation relation formula of described seismic impedance values of disparity and described chiltern stuff filling operation to carry out computing to the described seismic impedance numerical value of deriving along described seismic interpretation floor position in the described seismic reflection district chosen, obtain the filling operation of the fractured-vuggy reservoir in described seismic reflection district.
2. method according to claim 1, is characterized in that, described step 1 comprises:
According to horizon calibration in described log data and described logging data completed well;
Seismic interpretation is completed in conjunction with Seismic reflection character according to aboveground horizon calibration;
The secondary computing of seismic inversion is carried out according to described geological data;
Along seismic interpretation layer position, the inverting Wave Impedance Data Volume after secondary computing is cut into slices, obtain the described seismic inversion wave impedance planimetric map along described seismic interpretation layer position.
3. method according to claim 1, it is characterized in that, the functional relation of described seismic impedance values of disparity and chiltern stuff fill-up values of disparity is Y=a × X+b, wherein, Y is chiltern stuff fill-up values of disparity, X is seismic impedance values of disparity, and a, b are cross analysis gained constant.
4. method according to claim 1, it is characterized in that, the operation relation formula of described seismic impedance values of disparity and chiltern stuff filling operation is Z=(c × X+d)/V, wherein, Z is chiltern stuff filling operations, X is seismic impedance values of disparity, and c, d are cross analysis gained computing constant, and V is the volume of fractured-vuggy reservoir.
5. method according to claim 1, is characterized in that, in described step 2, extracts described sample spot in the different tectonic position in tectonic element and different tectonic element.
6. method according to claim 1, is characterized in that, in described step 4 and described step 5, described fill-up minimum calculation unit and the described filling operation minimum calculation unit of setting are all not less than research object precision prescribed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219252.8A CN104834007B (en) | 2015-05-04 | 2015-05-04 | Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219252.8A CN104834007B (en) | 2015-05-04 | 2015-05-04 | Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104834007A true CN104834007A (en) | 2015-08-12 |
CN104834007B CN104834007B (en) | 2017-09-26 |
Family
ID=53812003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510219252.8A Active CN104834007B (en) | 2015-05-04 | 2015-05-04 | Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104834007B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572746A (en) * | 2015-12-11 | 2016-05-11 | 中国石油天然气股份有限公司 | Method and device for determining drilling position |
CN109143397A (en) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | Carbonate reservoir fracture hole charges recognition methods and system |
CN109441541A (en) * | 2018-11-06 | 2019-03-08 | 中国矿业大学 | A kind of coal mine gob filling body carrying compression ratio monitoring system and its monitoring method |
CN110244363A (en) * | 2018-03-09 | 2019-09-17 | 中国石油化工股份有限公司 | A method of for predicting fracture hole type reservoir stock number |
CN111596351A (en) * | 2020-04-28 | 2020-08-28 | 中国石油天然气股份有限公司 | Quantitative evaluation method, system and device for carbonate rock transportation and conduction system and storage medium |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011154762A1 (en) * | 2010-06-07 | 2011-12-15 | Total Sa | Method for analyzing seismic data |
CN103064118A (en) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | Method of acoustic logging and quantifying cavern filling degree |
CN103116192A (en) * | 2012-07-30 | 2013-05-22 | 中国石油化工股份有限公司 | Fissure-cave type carbonate-reservoir reservoir body modeling method |
CN103529475A (en) * | 2013-04-19 | 2014-01-22 | 中国石油大学(华东) | Method for identifying and interpreting carbonate rock ancient karst reservoir layer three-dimensional structure |
-
2015
- 2015-05-04 CN CN201510219252.8A patent/CN104834007B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011154762A1 (en) * | 2010-06-07 | 2011-12-15 | Total Sa | Method for analyzing seismic data |
CN103116192A (en) * | 2012-07-30 | 2013-05-22 | 中国石油化工股份有限公司 | Fissure-cave type carbonate-reservoir reservoir body modeling method |
CN103064118A (en) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | Method of acoustic logging and quantifying cavern filling degree |
CN103529475A (en) * | 2013-04-19 | 2014-01-22 | 中国石油大学(华东) | Method for identifying and interpreting carbonate rock ancient karst reservoir layer three-dimensional structure |
Non-Patent Citations (2)
Title |
---|
苏劲 等: "地震资料解释与盆地模拟在油气成藏研究中的联合运用", 《石油地球物理勘探》 * |
赵军 等: "碳酸盐岩洞穴充填物及其充填程度的测井判别方法", 《中国岩溶》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572746A (en) * | 2015-12-11 | 2016-05-11 | 中国石油天然气股份有限公司 | Method and device for determining drilling position |
CN109143397A (en) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | Carbonate reservoir fracture hole charges recognition methods and system |
CN109143397B (en) * | 2017-06-28 | 2020-05-19 | 中国石油化工股份有限公司 | Carbonate reservoir fracture-cave filling identification method and system |
CN110244363A (en) * | 2018-03-09 | 2019-09-17 | 中国石油化工股份有限公司 | A method of for predicting fracture hole type reservoir stock number |
CN110244363B (en) * | 2018-03-09 | 2020-10-23 | 中国石油化工股份有限公司 | Method for predicting fracture-cavity reservoir resource amount |
CN109441541A (en) * | 2018-11-06 | 2019-03-08 | 中国矿业大学 | A kind of coal mine gob filling body carrying compression ratio monitoring system and its monitoring method |
CN109441541B (en) * | 2018-11-06 | 2020-01-03 | 中国矿业大学 | Coal mine goaf filling body bearing compression rate monitoring system and monitoring method thereof |
CN111596351A (en) * | 2020-04-28 | 2020-08-28 | 中国石油天然气股份有限公司 | Quantitative evaluation method, system and device for carbonate rock transportation and conduction system and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN104834007B (en) | 2017-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2960680B1 (en) | Fracturing and reactivated fracture volumes | |
US7933750B2 (en) | Method for defining regions in reservoir simulation | |
Ramírez et al. | Screening CO2 storage options in the Netherlands | |
Jenkins et al. | Quantifying and predicting naturally fractured reservoir behavior with continuous fracture models | |
Radwan | Three-dimensional gas property geological modeling and simulation | |
CN105223612B (en) | A kind of coal mine flood prediction and evaluation method based on earthquake information | |
RU2661489C1 (en) | Method of integrating initial data to update filtration structure of non-uniform carbonate reservoirs | |
Dance et al. | Illuminating the geology: Post-injection reservoir characterisation of the CO2CRC Otway site | |
CN102681013B (en) | Method and device for establishing carbonate reservoir space model | |
CN104992468A (en) | Three-dimensional geological modeling method for fracture-cave carbonate reservoir | |
CN104834007A (en) | Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process | |
CN104834008B (en) | Method for calculating filling degree of fractured-vuggy carbonate reservoir | |
CN105182424A (en) | Method and device for quantitatively predicting reservoir porosity based on patch saturation model | |
CN104977617A (en) | Reservoir fracture identification method and imaging logging reservoir fracture identification method | |
KR20200027264A (en) | A method for selecting a sweet spot in shale gas | |
EP3338115A1 (en) | Reservoir simulations with fracture networks | |
RU2586821C1 (en) | Distribution of petrographic facies using analytical simulation | |
CN108665545A (en) | Logging parameters three-dimensional geological model method for building up | |
Wilson et al. | Fracture model of the Upper Freeport coal: Marshall County West Virginia pilot ECBMR and CO2 sequestration site | |
Lyu et al. | Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China | |
Shen et al. | Hierarchical approach to modeling karst and fractures in carbonate karst reservoirs in the Tarim Basin | |
CN112184033A (en) | Carbonate rock stratum fault sealing evaluation method | |
Masoud et al. | Reservoir characterization and geostatistical model of the Cretaceous and Cambrian-Ordovician reservoir intervals, Meghil field, Sirte basin, Libya | |
Neele et al. | Options for CO2 sequestration in Kuwait | |
Price et al. | Dynamic calibration of the Shaikan Jurassic full-field fractured reservoir model through single-well DST and multi-well interference discrete fracture network simulation |
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 |