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 PDF

Info

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
Application number
CN201510219252.8A
Other languages
Chinese (zh)
Other versions
CN104834007B (en
Inventor
邢翔
张义杰
王召明
段书府
杨海军
贾进华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201510219252.8A priority Critical patent/CN104834007B/en
Publication of CN104834007A publication Critical patent/CN104834007A/en
Application granted granted Critical
Publication of CN104834007B publication Critical patent/CN104834007B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

The method of carbonate rock fractured cave type reservoir filling operation is calculated in seismic inversion process
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.
CN201510219252.8A 2015-05-04 2015-05-04 Method for calculating filling degree of carbonate fracture-cave reservoir in seismic inversion process Active CN104834007B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
苏劲 等: "地震资料解释与盆地模拟在油气成藏研究中的联合运用", 《石油地球物理勘探》 *
赵军 等: "碳酸盐岩洞穴充填物及其充填程度的测井判别方法", 《中国岩溶》 *

Cited By (8)

* Cited by examiner, † Cited by third party
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