CN111596351A - 碳酸盐岩输导体系定量评价方法、系统、装置及存储介质 - Google Patents
碳酸盐岩输导体系定量评价方法、系统、装置及存储介质 Download PDFInfo
- Publication number
- CN111596351A CN111596351A CN202010348654.9A CN202010348654A CN111596351A CN 111596351 A CN111596351 A CN 111596351A CN 202010348654 A CN202010348654 A CN 202010348654A CN 111596351 A CN111596351 A CN 111596351A
- Authority
- CN
- China
- Prior art keywords
- transmission
- dimensional
- work area
- data volume
- dimensional seismic
- 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 88
- 239000011435 rock Substances 0.000 title claims abstract description 71
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 55
- 238000011158 quantitative evaluation Methods 0.000 title claims abstract description 37
- 238000003860 storage Methods 0.000 title claims abstract description 10
- 230000005540 biological transmission Effects 0.000 claims abstract description 349
- 239000004020 conductor Substances 0.000 claims abstract description 248
- 238000009826 distribution Methods 0.000 claims abstract description 116
- 230000000704 physical effect Effects 0.000 claims abstract description 79
- 238000004364 calculation method Methods 0.000 claims abstract description 50
- 230000035699 permeability Effects 0.000 claims abstract description 24
- 238000012512 characterization method Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims description 64
- 238000004458 analytical method Methods 0.000 claims description 32
- 238000011161 development Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 230000004927 fusion Effects 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 8
- 230000026683 transduction Effects 0.000 claims description 8
- 238000010361 transduction Methods 0.000 claims description 8
- 239000011148 porous material Substances 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 239000011324 bead Substances 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 3
- 230000000452 restraining effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000005553 drilling Methods 0.000 claims 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 2
- 230000014509 gene expression Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 17
- 238000005516 engineering process Methods 0.000 description 12
- 230000008901 benefit Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 235000019994 cava Nutrition 0.000 description 3
- 239000003086 colorant Substances 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000010459 dolomite Substances 0.000 description 2
- 229910000514 dolomite Inorganic materials 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 239000006028 limestone Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6226—Impedance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6244—Porosity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6246—Permeability
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种碳酸盐岩输导体系定量评价方法、系统、装置及存储介质。该方法包括:获取工区的高精度三维地震振幅数据体和三维地震反演的波阻抗数据体;基于振幅数据体和/或波阻抗数据体确定输导体空间展布;获取工区不同类型的输导体对应的多种地震属性门槛值,基于振幅数据体,利用多种地震属性门槛值,在输导体空间展布范围内进行输导体类型划分,得到工区输导体系类型的三维空间分布;获取不同类型的输导体对应的物性计算公式,物性计算公式为物性与波阻抗的关系式,物性包括孔隙度和渗透率中的至少一种;基于波阻抗数据体,利用工区输导体系类型的三维空间分布计算工区输导体的物性,从而实现工区输导体的宏观输导性能表征。
Description
技术领域
本发明涉及石油勘探技术领域,特别涉及一种利用三维地震资料进行的碳酸盐岩输导体系定量评价方法、系统、装置及存储介质。
背景技术
输导体系是指含油气系统中所有的运移通道及相关围岩的总和。碳酸盐岩输导体系可分为溶洞、溶孔、裂缝和断裂共4种输导体类型。在某一系统中,输导体系并非为单一类型,而是多种类型的组合。输导体系与油气聚集相辅相成,含油气系统的输导体系不同则油气运聚的方式各异,反过来油气的不同运聚方式又可以改造甚至形成新的输导体系。由于输导体系具有相对独立性、时空性和复杂性,所以在研究输导体系时,应综合研究其三维空间分布及评价其输导能力。输导体系识别预评价是根据勘探开发的需要,应用现有的资料和认识,采用多种技术手段,对认为是有效输导体系的某些(或全部)特征进行预测,主要是以地震信息为主要依据,综合利用其他资料(地质、测井、岩石物理等)作为约束,对输导体系的空间分布、输导性能等进行预测的一项专门技术。
目前三维地震预测输导体系的主要工作内容大体分为四个方面:一是输导体系岩性预测;二是输导体系形态预测;三是输导体系物性预测;四是输导体系含油气性综合分析。目前输导体系地震预测技术主要有地震正反演、地震相技术、属性分析、频谱分析、时频分析、多波和烃类检测技术等。而随着技术的成熟,地震反演、地震属性分析、AVO分析等技术成为了主要的技术手段,是预测地下输导体系岩性、物性的技术关键。
尽管地震输导体系预测技术取得了骄人的成绩,但就目前而言,对输导体系研究还不能满足目前的需求。例如对于至关重要的孔隙度和渗透率这两个物性参数,仅仅能够对孔隙度做出比较准确的预测,而不能直接根据地震资料得出渗透率。目前行业内仍然缺乏真正把地震信息与地质、测井信息融合在一起,运用多种数学、软件工具进行可视化描述的高效方法。
发明内容
针对现有技术的不足,本发明的目的在于提供一种定量评价碳酸盐岩输导体系的方法,实现精确定量评价碳酸盐岩输导体系的物性。
为了实现上述目的,本发明提供了一种碳酸盐岩输导体系定量评价方法,其中,该方法包括:
获取工区的高精度三维地震振幅数据体;其中,所述高精度三维地震振幅数据体为主频不低于预设主频的三维地震振幅数据体,所述预设主频不低于30Hz;
获取工区的三维地震反演的波阻抗数据体;
基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布;
确定工区输导体系类型的三维空间分布:获取工区不同类型的输导体对应的多种地震属性门槛值;基于所述高精度三维地震振幅数据体,利用所述不同类型的输导体对应的多种地震属性门槛值,在确定的输导体的空间展布范围内进行输导体类型划分,得到工区输导体系类型的三维空间分布;其中,所述多种地震属性至少包括振幅变化率、主频和频率衰减因子;
确定工区输导体的宏观输导性能:获取不同类型的输导体对应的物性计算公式,其中,物性计算公式为物性与波阻抗的关系式,物性包括孔隙度和渗透率中的至少一种;基于所述三维地震反演的波阻抗数据体,利用不同类型的输导体对应的物性计算公式计算所述工区输导体系类型的三维空间分布中输导体的物性,从而实现工区输导体的宏观输导性能表征。
在上述碳酸盐岩输导体系定量评价方法中,优选地,所述预设主频不低于50Hz。更优选地,所述预设主频为50Hz-60Hz。
在上述碳酸盐岩输导体系定量评价方法中,优选地,所述获取工区的高精度三维地震振幅数据体包括:
获取工区的三维地震振幅数据体;
判断所述工区的三维地震振幅数据体的主频是否低于所述预设主频;
若所述工区的三维地震振幅数据体的主频不低于预设主频,则所述工区的三维地震振幅数据体即为所述工区的高精度三维地震振幅数据体;
若所述工区的三维地震振幅数据体的主频低于预设主频,则对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理产生主频不低于所述预设主频的所述工区的高精度三维地震振幅数据体。
在一具体实施方式中,采用采用超高分辨率处理方法(HFI,High FrequencyImaging)对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理。
在上述碳酸盐岩输导体系定量评价方法中,优选地,所述获取工区的三维地震反演的波阻抗数据体包括:基于获取的高精度三维地震振幅数据体进行反演得到所述工区的三维地震反演的波阻抗数据体。
在上述碳酸盐岩输导体系定量评价方法中,优选地,所述基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布包括:
基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布;
基于高精度三维地震振幅数据体识别输导体系横向展布,对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定。
该优选方法利用三维地震资料进行,将属性分析和相控反演相结合,更有助于提高三维输导体系的刻画精度。具体地,该优选方案中,在确定输导体的空间展布过程中,通过平面与剖面结合的双面工作模式,利用三维地震反演的波阻抗数据体以及三维地震振幅数据体的结合,将属性和反演一体化实现了输导体空间展布的精确刻画。将属性分析和相控反演有机结合进行输导体空间展布,可以兼顾属性和反演对三维输导体系刻画的优势。
在一具体实施方式中,基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布可以采用下述方式进行:利用三维地震反演的波阻抗数据体的反演数据值信息,通过门槛值划分,初步定量预测三维输导体系包括厚度在内的空间展布。
在一具体实施方式中,基于高精度三维地震振幅数据体识别输导体系横向展布可以采用下述方式进行:基于高精度三维地震振幅数据体,利用地震属性-泊松分频RGB融合属性识别输导体系横向展布。该优选方案将地震多属性分析技术、地震时频分析技术以及小波分析技术相结合,更有利于为地震资料的全三维解释和定量解释提供丰富信息。
在上述碳酸盐岩输导体系定量评价方法中,优选地,输导体的类型包括溶洞、溶孔、裂缝和断裂。
在上述碳酸盐岩输导体系定量评价方法中,优选地,所述获取工区不同类型的输导体对应的多种地震属性门槛值包括:
采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,在所述已钻井钻遇的输导体中确定出厚度超过所述高精度三维地震振幅数据体调谐厚度的输导体作为第一输导体,并划分第一输导体的输导体类型;
利用时深转换的技术方法,在所述高精度三维地震振幅数据体中标定所述第一输导体的位置;基于标定的所述第一输导体的位置,确定出第一输导体中不同类型的输导体对应的所述多种地震属性门槛值。
在上述碳酸盐岩输导体系定量评价方法中,优选地,所述获取不同类型的输导体对应的物性计算公式采用下述方式实现:
采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,划分已钻井钻遇的输导体的输导体类型;
获取输导体的空间展布范围内已钻井的岩心物性测试结果;
按照输导体类型,分别进行物性与波阻抗关系拟合,得到不同类型的输导体对应的物性计算公式。
本发明还提供了一种碳酸盐岩输导体系定量评价系统,其中,该系统包括:
高精度三维地震振幅数据体获取模块:用于获取工区的高精度三维地震振幅数据体;其中,所述高精度三维地震振幅数据体为主频不低于预设主频的三维地震振幅数据体,所述预设主频不低于30Hz;
三维地震反演的波阻抗数据体获取模块:用于获取工区的三维地震反演的波阻抗数据体;
输导体的空间展布分析模块:基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布;
输导体系类型的三维空间分布分析模块;其中,所述输导体系类型的三维空间分布分析模块包括多种地震属性门槛值获取子模块和输导体类型划分子模块;所述多种地震属性门槛值获取子模块用于获取工区不同类型的输导体对应的多种地震属性门槛值;所述输导体类型划分子模块用于基于所述高精度三维地震振幅数据体,利用所述不同类型的输导体对应的多种地震属性门槛值,在确定的输导体的空间展布范围内进行输导体类型划分,得到工区输导体系类型的三维空间分布;其中,所述多种地震属性至少包括振幅变化率、主频和频率衰减因子;
输导体的宏观输导性能分析模块;其中,输导体的宏观输导性能分析模块包括物性计算公式获取子模块和输导体物性计算子模块;所述物性计算公式获取子模块用于获取不同类型的输导体对应的物性计算公式,其中,物性计算公式为物性与波阻抗的关系式,物性包括孔隙度和渗透率中的至少一种;所述输导体物性计算子模块用于基于所述三维地震反演的波阻抗数据体,利用所述工区输导体系类型的三维空间分布计算工区输导体的物性,从而实现工区输导体的宏观输导性能表征。
在上述碳酸盐岩输导体系定量评价系统中,优选地,所述预设主频不低于50Hz。更优选地,所述预设主频为50Hz-60Hz。
在上述碳酸盐岩输导体系定量评价系统中,优选地,所述高精度三维地震振幅数据体获取模块包括:
第一获取单元:用于获取工区的三维地震振幅数据体;
第一判别单元:判断第一获取单元获取的所述工区的三维地震振幅数据体的主频是否低于所述预设主频;若第一判别单元判断得出所述工区的三维地震振幅数据体的主频不低于预设主频,则采用第一处理单元进行处理;若所述工区的三维地震振幅数据体的主频低于预设主频,则采用第二处理单元进行处理;
第一处理单元:用于实施将所述工区的三维地震振幅数据体作为所述工区的高精度三维地震振幅数据体的操作。
第二处理单元:用于实施对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理产生主频不低于所述预设主频的所述工区的高精度三维地震振幅数据体的操作。
在一具体实施方式中,所述第二处理单元采用超高分辨率处理方法(HFI,HighFrequency Imaging)对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理。
在上述碳酸盐岩输导体系定量评价系统中,优选地,所述三维地震反演的波阻抗数据体获取模块用于基于获取的高精度三维地震振幅数据体进行反演得到所述工区的三维地震反演的波阻抗数据体。
在上述碳酸盐岩输导体系定量评价系统中,优选地,所述输导体的空间展布分析模块包括:
第三处理单元:用于基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布;
第四处理单元:用于基于高精度三维地震振幅数据体识别输导体系横向展布,对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定。
在一具体实施方式中,所述第三处理单元用于利用三维地震反演的波阻抗数据体的反演数据值信息,通过门槛值划分,初步定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布。
在一具体实施方式中,所述第四处理单元用于基于高精度三维地震振幅数据体,利用地震属性-泊松分频RGB融合属性识别输导体系横向展布,对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定。
在上述碳酸盐岩输导体系定量评价系统中,优选地,输导体的类型包括溶洞、溶孔、裂缝和断裂。
在上述碳酸盐岩输导体系定量评价系统中,优选地,所述多种地震属性门槛值获取子模块包括:
第二获取单元:用于获取输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料;
第五处理单元:用于采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,在所述已钻井钻遇的输导体中确定出厚度超过所述高精度三维地震振幅数据体调谐厚度的输导体作为第一输导体,并划分第一输导体的输导体类型;
第六处理单元:用于利用时深转换的技术方法,在所述高精度三维地震振幅数据体中标定所述第一输导体的位置;基于标定的所述第一输导体的位置,确定出第一输导体中不同类型的输导体的所述多种地震属性门槛值。
在上述碳酸盐岩输导体系定量评价系统中,优选地,所述物性计算公式获取子模块包括:
第三获取单元:用于获取输导体的空间展布范围内已钻井的岩心物性测试结果
第七处理单元:用于采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,划分已钻井钻遇的输导体的输导体类型;
第八处理单元:用于按照输导体类型,分别进行物性与波阻抗关系拟合,得到不同类型的输导体对应的物性计算公式。
本发明还提供了一种碳酸盐岩输导体系定量评价装置,包括处理器及存储器;其中,
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现上述碳酸盐岩输导体系定量评价方法的步骤。
本发明还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现上述碳酸盐岩输导体系定量评价方法的步骤。
本发明提供的上述技术方案,将三维输导体系进行分类,建立地震信息属性体与三维输导体系类型参数间的相关性,实现地震资料的全三维解释得到工区输导体系类型的三维空间分布;基于物性与波阻抗的关系式,按照类别计算输导体的物性,实现工区输导体的宏观输导性能表征,将属性分析与相控反演有机结合,充分发挥两种技术的优势,更准确的表征输导体的宏观输导性能。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图:
图1为本发明一实施例提供的碳酸盐岩输导体系定量评价方法的流程示意图;
图2为本发明一实施例提供的碳酸盐岩输导体系定量评价方法中获取工区的高精度三维地震振幅数据步骤的优化示意图;
图3为本发明一实施例提供的碳酸盐岩输导体系定量评价方法的流程示意图;
图4为本发明一实施例提供的碳酸盐岩输导体系定量评价系统的结构示意图;
图5为本发明一实施例提供的碳酸盐岩输导体系定量评价系统中高精度三维地震振幅数据体获取模块的结构示意图;
图6为本发明一实施例提供的碳酸盐岩输导体系定量评价系统中输导体的空间展布分析模块的结构示意图;
图7为本发明一实施例提供的碳酸盐岩输导体系定量评价系统中多种地震属性门槛值获取子模块的结构示意图;
图8为本发明一实施例提供的碳酸盐岩输导体系定量评价系统中物性计算公式获取子模块的结构示意图;
图9为本发明一实施例提供的碳酸盐岩输导体系定量评价装置的结构示意图;
图10为本发明一实施例中输导体类型的三维空间分布图;
图11为本发明一实施例中孔隙度与波阻抗关系图;
图12为本发明一实施例中孔隙度与渗透率关系图;
图13为本发明一实施例中工区输导体的渗透率分布图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明的保护范围。
下面参考本发明的若干代表性实施方式,详细阐述本发明的原理和精神。
参见图1,本发明的一实施例提供了一种碳酸盐岩输导体系定量评价方法,该方法包括:
步骤S1:获取工区的高精度三维地震振幅数据体;其中,所述高精度三维地震振幅数据体为主频不低于预设主频的三维地震振幅数据体,所述预设主频不低于30Hz;
步骤S2:获取工区的三维地震反演的波阻抗数据体;
步骤S3:确定输导体的空间展布:基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布;
步骤S4:确定工区输导体系类型的三维空间分布:获取工区不同类型的输导体对应的多种地震属性门槛值;基于所述高精度三维地震振幅数据体,利用所述不同类型的输导体对应的多种地震属性门槛值,在确定的输导体的空间展布范围内进行输导体类型划分,得到工区输导体系类型的三维空间分布;其中,所述多种地震属性至少包括振幅变化率、主频和频率衰减因子;
步骤S5:确定工区输导体的宏观输导性能:获取不同类型的输导体对应的物性计算公式,其中,物性计算公式为物性与波阻抗的关系式,物性包括孔隙度和渗透率中的至少一种;基于所述三维地震反演的波阻抗数据体,利用不同类型的输导体对应的物性计算公式计算所述工区输导体系类型的三维空间分布中输导体的物性,从而实现工区输导体的宏观输导性能表征。
目前石油勘探领域的地震数据体采用25m×25m道间距,采集了0-6s的反射波振幅数据,采用2ms的垂向采样间隔。在本发明实施例的串珠体空间发育规律分析方法中,根据识别三维输导体系的技术需求,地震数据振幅数据体在油气藏发育层段的主频应超过预设主频,在本技术中最低的预设主频应30Hz以上,通常可以设定预设主频在50Hz以上,优选的预设主频可以设定为50Hz-60Hz。
参见图2,进一步,步骤S1中获取工区的高精度三维地震振幅数据体可以进一步包括如下步骤:
S11:获取工区的三维地震振幅数据体;
S12:判断所述工区的三维地震振幅数据体的主频是否低于所述预设主频;
S13:若所述工区的三维地震振幅数据体的主频不低于预设主频,则所述工区的三维地震振幅数据体即为所述工区的高精度三维地震振幅数据体;
S14:若所述工区的三维地震振幅数据体的主频低于预设主频,则对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理产生主频不低于所述预设主频的所述工区的高精度三维地震振幅数据体;
其中,对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理可以采用诸如超高分辨率处理方法(HFI,High Frequency Imaging)等方式进行。
进一步,在步骤S2中,所述获取工区的三维地震反演的波阻抗数据体可以通过下述方式实现:基于获取的高精度三维地震振幅数据体进行反演得到所述工区的三维地震反演的波阻抗数据体。此处,可以采用常规反演方法,利用高精度三维地震振幅数据体获取三维地震反演的波阻抗数据体。
进一步,在步骤S3中,所述基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布可以包括:
基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布;
基于高精度三维地震振幅数据体识别输导体系横向展布,对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定;
其中,所述基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布可以采用下述方式进行:利用三维地震反演的波阻抗数据体的反演数据值信息,通过门槛值划分,初步定量预测三维输导体系的空间分布(包括厚度)。三维地震反演的波阻抗数据体有较好的垂向分辨率,对三维输导体的厚度有很高的刻画能力。
利用三维地震反演的波阻抗数据体的反演数据值信息,通过门槛值划分,初步定量预测三维输导体系的空间分布可以采用本领域常用方式进行。例如,利用单井的测井数据识别井筒周围缝洞体,然后通过井震的时深转换,将深度域井筒的缝洞体标定到时间域的反演数据体上,根据缝洞体发育范围内的波阻抗数据,确定门槛值,然后通过门槛值初步定量预测三维输导体系的空间分布。
其中,所述基于三维地震振幅数据体识别输导体系横向展布可以采用下述方式进行:基于三维地震振幅数据体,利用地震属性-泊松分频RGB融合属性识别输导体系横向展布。高精度的三维地震振幅数据体的地震属性有好的横向分辨率,所以地震属性能对三维输导体系的外形和结构有较好的刻画,解释人员也容易解释输导体系展布方向和范围。
基于三维地震振幅数据体,利用地震属性-泊松分频RGB融合属性识别输导体系横向展布可以采用本领域常用方式进行。例如,在15、35、50Hz共3套频率调谐数据体基础上进行RGB混色,自然界中的可见光由红(R)、绿(G)、蓝(B)3种基色按照不同比例合成产生,每一种颜色按其亮度的不同分为256个等级,各参数的取值范围为0-255;在RGB色彩空间中,将3种基色按不同比例混合,形成各种不同的色彩图像,定义为映射函数S,S可对输入的3种基色比例进行调节,最终形成的融合图中每个值对应某一特定的颜色,以此识别输导体系横向展布。
进一步,输导体的类型可以包括溶洞、溶孔、裂缝和断裂。
进一步,在步骤S4中,所述获取工区不同类型的输导体对应的多种地震属性门槛值包括:
采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,在所述已钻井钻遇的输导体中确定出厚度超过所述高精度三维地震振幅数据体调谐厚度(例如5-10m)的输导体作为第一输导体,并划分第一输导体的输导体类型;
利用时深转换的技术方法,在所述高精度三维地震振幅数据体中标定所述第一输导体的位置;基于标定的所述第一输导体的位置,确定出第一输导体中不同类型的输导体的所述多种地震属性门槛值;
所述确定出第一输导体中不同类型的输导体的所述多种地震属性门槛值可以采用数据统计方式进行;例如支持向量机算法,但不限于次。
进一步,在步骤S5中,所述获取不同类型的输导体对应的物性计算公式可以采用下述方式实现:
采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,划分已钻井钻遇的输导体的输导体类型;
获取输导体的空间展布范围内已钻井的岩心物性测试结果;
按照输导体类型,分别进行物性与波阻抗关系拟合,得到不同类型的输导体对应的物性计算公式。
在一实施方式中,步骤S5包括:
采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,划分已钻井钻遇的输导体的输导体类型;
获取输导体的空间展布范围内已钻井的岩心渗透率和孔隙度测试结果;
按照输导体类型,分别进行孔隙度与波阻抗关系拟合,得到不同类型的输导体对应的孔隙度计算公式;
按照输导体类型分别进行渗透率与孔隙度关系拟合,得到不同类型的输导体对应的渗透率计算公式;
基于所述三维地震反演的波阻抗数据体,利用不同类型的输导体对应的孔隙度计算公式计算所述工区输导体系类型的三维空间分布中输导体的孔隙度;然后再利用不同类型的输导体对应的渗透率计算公式计算所述工区输导体系类型的三维空间分布中输导体的渗透率;从而实现工区输导体的宏观输导性能表征。
进一步,所述碳酸盐岩输导体系定量评价方法还可以包括:利用确定的输导体的空间展布计算输导体的体积。
输导体的体积也可以用来表征输导体的宏观输导性能;通常可以认为超过300-500m3的输导体系为有效输导体系。
中国深层碳酸盐岩储集层主要由古生代海相地层组成,具有很强的非均质性和较差的物性。塔中油田是中国深埋碳酸盐岩油藏的典型代表油田。塔中油田的主要油气发现集中在奥陶系缝洞型油藏中,深度超过5500米。本领域技术人员在先前的研究中心确定了奥陶纪地层层序地层格架,并分析了该油藏的沉积相;分析了层间岩溶储集层的发育特征,并将上鹰山组垂直划分为四个区域;利用碳酸盐岩储层的主要控制因素将其划分为礁滩型储层,风化壳和热液蚀变型储层。以前的文献表明,缝洞型储层在埋藏历史上经历了许多地质转变,包括地表岩溶,断层和深部热液岩溶作用,严重改善了古碳酸盐岩地层。缝洞型储层的储集空间主要由溶解的小孔隙和裂缝组成,总体上物理性能较差。因此,塔中油田奥陶纪缝洞型油藏的识别和物性表征是油气成藏研究的问题。应用本发明提供的碳酸盐岩输导体系定量评价方法进行塔中油田碳酸盐岩输导体系进行定量评价,能够有效实现输导体系的定量评价。
参见图3,本发明的一实施例提供了碳酸盐岩输导体系定量评价方法,该方法用以对塔中油田碳酸盐岩输导体系进行定量评价;该方法包括:
步骤S1,获取工区的高精度三维地震振幅数据体:
步骤S11,获取工区的三维地震振幅数据体;
步骤S12,判断所述工区的三维地震振幅数据体的主频是否低于所述预设主频;所述预设主频为50Hz-60Hz;
步骤S13,若所述工区的三维地震振幅数据体的主频不低于预设主频,则所述工区的三维地震振幅数据体即为工区的高精度三维地震振幅数据体;
步骤S14,若所述工区的三维地震振幅数据体的主频低于预设主频,则对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理产生主频不低于所述预设主频的工区的高精度三维地震振幅数据体;
步骤S2,获取工区的三维地震反演的波阻抗数据体:基于获取的高精度三维地震振幅数据体进行反演得到工区的三维地震反演的波阻抗数据体;
步骤S3,确定输导体的空间展布:
步骤S31,利用三维地震反演的波阻抗数据体的反演数据值信息,通过门槛值划分,初步定量预测三维输导体系的空间分布(包括厚度);
步骤S32,基于三维地震振幅数据体,利用地震属性-泊松分频RGB融合属性识别输导体系横向展布对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定;
步骤S4,确定工区输导体系类型的三维空间分布:
步骤S41,采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,在所述已钻井钻遇的输导体中确定出厚度超过所述高精度三维地震振幅数据体调谐厚度(例如5-10m)的输导体作为第一输导体,并划分第一输导体的输导体类型具体包括溶洞、灰岩溶孔、白云岩溶孔、裂缝、断裂和原岩共6种输导体类型(参见表1);
步骤S42,利用时深转换的技术方法,在所述高精度三维地震振幅数据体中标定所述第一输导体的位置;基于标定的所述第一输导体的位置,确定出第一输导体中溶洞、灰岩溶孔、白云岩溶孔、裂缝、断裂和原岩共6种输导体类型对应的振幅变化率、主频和频率衰减因子的门槛值作为工区不同类型的输导体对应的多种地震属性门槛值;
步骤S43,基于所述高精度三维地震振幅数据体,利用所述不同类型的输导体对应的多种地震属性门槛值,在确定的输导体的空间展布范围内进行输导体类型划分,得到工区输导体系类型的三维空间分布(如图10所示);
步骤S5,确定工区输导体的宏观输导性能:
步骤S51,采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,划分已钻井钻遇的输导体的输导体类型;
步骤S52,获取输导体的空间展布范围内已钻井的岩心渗透率和孔隙度测试结果;
步骤S53,按照输导体类型,分别进行孔隙度与波阻抗关系拟合,得到不同类型的输导体对应的孔隙度计算公式;按照输导体类型分别进行渗透率与孔隙度关系拟合,得到不同类型的输导体对应的渗透率计算公式(参见图11、图12以及表2);
步骤S54,基于所述三维地震反演的波阻抗数据体,利用不同类型的输导体对应的渗透率计算公式计算所述工区输导体系类型的三维空间分布中输导体的渗透率;然后再利用不同类型的输导体对应的孔隙度计算公式计算所述工区输导体系类型的三维空间分布中输导体的孔隙度;从而实现工区输导体的宏观输导性能表征(渗透率分布参见图13)。
表1
表2
序号 | 表示符号 | 孔隙度计算公式 | 渗透率计算公式 |
1 | SF1 | Por=-11.31·Ln(AI)+111.24 | K=0.0106·e<sup>(1.0178·Por)</sup> |
2 | SF2 | Por=-25.82·Ln(AI)+251.38 | K=0.002·e<sup>(1.1601·Por)</sup> |
3 | SF3 | Por=-43.97·Ln(AI)+429.96 | K=0.0079·e<sup>(0.4239·Por)</sup> |
4 | SF4 | Por=-14.51·Ln(AI)+140.57 | K=0.0808·e<sup>(1.0878·Por)</sup> |
5 | SF5 | Por=-15.35·Ln(AI)+149.48 | K=0.3651·e<sup>(1.5724·Por)</sup> |
6 | SF6 | Por=-38.33·Ln(AI)+369.32 | K=0.0042·e<sup>(0.6631·Por)</sup> |
其中,Por为孔隙度,单位%;AI为波阻抗,单位104g/cm3·m/s;K为渗透率,单位md。
本发明实施例还提供了一种碳酸盐岩输导体系定量评价系统,优选地,该系统用于实现上述的方法实施例。
图4是根据本发明实施例的串珠体空间发育规律分析系统的结构框图,如图4所示,该装置包括:高精度三维地震振幅数据体获取模块41、三维地震反演的波阻抗数据体获取模块42、输导体的空间展布分析模块43、输导体系类型的三维空间分布分析模块44、输导体的宏观输导性能分析模块45;
高精度三维地震振幅数据体获取模块41:用于获取工区的高精度三维地震振幅数据体;其中,所述高精度三维地震振幅数据体为主频不低于预设主频的三维地震振幅数据体,所述预设主频不低于30Hz;
三维地震反演的波阻抗数据体获取模块42:用于获取工区的三维地震反演的波阻抗数据体;
输导体的空间展布分析模块43:基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布;
所述输导体系类型的三维空间分布分析模块44包括多种地震属性门槛值获取子模块441和输导体类型划分子模块442;所述多种地震属性门槛值获取子模块441用于获取工区不同类型的输导体对应的多种地震属性门槛值;所述输导体类型划分子模块442用于基于所述高精度三维地震振幅数据体,利用所述不同类型的输导体对应的多种地震属性门槛值,在确定的输导体的空间展布范围内进行输导体类型划分,得到工区输导体系类型的三维空间分布;其中,所述多种地震属性至少包括振幅变化率、主频和频率衰减因子;
输导体的宏观输导性能分析模块45包括物性计算公式获取子模块451和输导体物性计算子模块452;所述物性计算公式获取子模块451用于获取不同类型的输导体对应的物性计算公式,其中,物性计算公式为物性与波阻抗的关系式,物性包括孔隙度和渗透率中的至少一种;所述输导体物性计算子模块452用于基于所述三维地震反演的波阻抗数据体,利用所述工区输导体系类型的三维空间分布计算工区输导体的物性,从而实现工区输导体的宏观输导性能表征。
上述高精度三维地震振幅数据体获取模块51中,预设主频可以设为不低于30Hz,例如50Hz-60Hz。
进一步,继续参见图5,上述高精度三维地震振幅数据体获取模块41可以包括:
第一获取单元411:用于获取工区的三维地震振幅数据体;
第一判别单元412:判断第一获取单元411获取的所述工区的三维地震振幅数据体的主频是否低于所述预设主频;若第一判别单元412判断得出所述工区的三维地震振幅数据体的主频不低于预设主频,则采用第一处理单元413进行处理;若所述工区的三维地震振幅数据体的主频低于预设主频,则采用第二处理单元414进行处理;
第一处理单元413:用于实施将所述工区的三维地震振幅数据体作为所述工区的高精度三维地震振幅数据体的操作。
第二处理单元414:用于实施对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理产生主频不低于所述预设主频的所述工区的高精度三维地震振幅数据体的操作。
其中,第二处理单元414可以采用超高分辨率处理方法(HFI,High FrequencyImaging)对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理。
进一步,三维地震反演的波阻抗数据体获取模块42用于基于获取的高精度三维地震振幅数据体进行反演得到所述工区的三维地震反演的波阻抗数据体。
进一步,继续参见图6,所述输导体的空间展布分析模块43还可以包括:
第三处理单元431:用于基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布;
第四处理单元432:用于基于高精度三维地震振幅数据体识别输导体系横向展布,对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定。
其中,所述第三处理单元431还可以用于利用三维地震反演的波阻抗数据体的反演数据值信息,通过门槛值划分,初步定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布。
其中,所述第四处理单元432还可以用于基于高精度三维地震振幅数据体,利用地震属性-泊松分频RGB融合属性识别输导体系横向展布,对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定。
进一步,输导体的类型还可以包括溶洞、溶孔、裂缝和断裂。
进一步,继续参见图7,多种地震属性门槛值获取子模块441可以包括:
第二获取单元4411:用于获取输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料;
第五处理单元4412:用于采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,在所述已钻井钻遇的输导体中确定出厚度超过所述高精度三维地震振幅数据体调谐厚度的输导体作为第一输导体,并划分第一输导体的输导体类型;
第六处理单元4413:用于利用时深转换的技术方法,在所述高精度三维地震振幅数据体中标定所述第一输导体的位置;基于标定的所述第一输导体的位置,确定出第一输导体中不同类型的输导体的所述多种地震属性门槛值。
进一步,继续参见图8,所述物性计算公式获取子模块451包括:
第三获取单元4511:用于获取输导体的空间展布范围内已钻井的岩心物性测试结果
第七处理单元4512:用于采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,划分已钻井钻遇的输导体的输导体类型;
第八处理单元4513:用于按照输导体类型,分别进行物性与波阻抗关系拟合,得到不同类型的输导体对应的物性计算公式。
进一步,该装置可以进一步包括:输导体体积计算模块46,该模块用于基于工区输导体系类型的三维空间分布计算工区输导体体积。
图9是根据本发明实施例的碳酸盐岩输导体系定量评价装置的示意图。图9所示的碳酸盐岩输导体系定量评价装置为通用数据处理装置,其包含通用的计算机硬件结构,其至少包含处理器1000、存储器1111;所述处理器1000用于执行所述存储器中存储的分子结构生成程序,以实现各方法实施例所述的碳酸盐岩输导体系定量评价方法(具体方法参见上述方法实施例的描述,在此不再赘述)。
本发明实施例还提供了一种计算机可读存储介质,所述存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现各方法实施例所述的碳酸盐岩输导体系定量评价方法(具体方法参见上述方法实施例的描述,在此不再赘述)。
以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的,因此权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外,由于本领域的技术人员容易想到很多修改和改变,因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作,而是可以涵盖落入其范围内的所有合适修改和等同物。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (20)
1.一种碳酸盐岩输导体系定量评价方法,其中,该方法包括:
获取工区的高精度三维地震振幅数据体;其中,所述高精度三维地震振幅数据体为主频不低于预设主频的三维地震振幅数据体,所述预设主频不低于30Hz(优选为不低于50Hz);
获取工区的三维地震反演的波阻抗数据体;
基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布;
确定工区输导体系类型的三维空间分布:获取工区不同类型的输导体对应的多种地震属性门槛值;基于所述高精度三维地震振幅数据体,利用所述不同类型的输导体对应的多种地震属性门槛值,在确定的输导体的空间展布范围内进行输导体类型划分,得到工区输导体系类型的三维空间分布;其中,所述多种地震属性至少包括振幅变化率、主频和频率衰减因子;
确定工区输导体的宏观输导性能:获取不同类型的输导体对应的物性计算公式,其中,物性计算公式为物性与波阻抗的关系式,物性包括孔隙度和渗透率中的至少一种;基于所述三维地震反演的波阻抗数据体,利用所述工区输导体系类型的三维空间分布计算工区输导体的物性,从而实现工区输导体的宏观输导性能表征。
2.根据权利要求1所述的方法,其中,所述获取工区的高精度三维地震振幅数据体包括:
获取工区的三维地震振幅数据体;
判断所述工区的三维地震振幅数据体的主频是否低于所述预设主频;
若所述工区的三维地震振幅数据体的主频不低于预设主频,则所述工区的三维地震振幅数据体即为所述工区的高精度三维地震振幅数据体;
若所述工区的三维地震振幅数据体的主频低于预设主频,则对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理产生主频不低于所述预设主频的所述工区的高精度三维地震振幅数据体。
3.根据权利要求1所述的方法,其中,所述获取工区的三维地震反演的波阻抗数据体包括:基于获取的高精度三维地震振幅数据体进行反演得到所述工区的三维地震反演的波阻抗数据体。
4.根据权利要求1所述的方法,其中,所述基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布包括:
基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布;
基于高精度三维地震振幅数据体识别输导体系横向展布,对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定。
5.根据权利要求4所述的方法,其中,所述基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布可以采用下述方式进行:利用三维地震反演的波阻抗数据体的反演数据值信息,通过门槛值划分,初步定量预测三维输导体系包括厚度在内的空间展布。
6.根据权利要求4所述的方法,其中,所述基于高精度三维地震振幅数据体识别输导体系横向展布可以采用下述方式进行:基于高精度三维地震振幅数据体,利用地震属性-泊松分频RGB融合属性识别输导体系横向展布。
7.根据权利要求1所述的方法,其中,所述输导体的类型包括溶洞、溶孔、裂缝和断裂。
8.根据权利要求1所述的方法,其中,所述获取工区不同类型的输导体对应的多种地震属性门槛值包括:
采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,在所述已钻井钻遇的输导体中确定出厚度超过所述高精度三维地震振幅数据体调谐厚度的输导体作为第一输导体,并划分第一输导体的输导体类型;
利用时深转换的技术方法,在所述高精度三维地震振幅数据体中标定所述第一输导体的位置;基于标定的所述第一输导体的位置,确定出第一输导体中不同类型的输导体的所述多种地震属性门槛值。
9.根据权利要求1所述的方法,其中,所述获取不同类型的输导体对应的物性计算公式采用下述方式实现:
采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,划分已钻井钻遇的输导体的输导体类型;
获取输导体的空间展布范围内已钻井的岩心物性测试结果;
按照输导体类型,分别进行物性与波阻抗关系拟合,得到不同类型的输导体对应的物性计算公式。
10.一种碳酸盐岩输导体系定量评价系统,其中,该系统包括:
高精度三维地震振幅数据体获取模块:用于获取工区的高精度三维地震振幅数据体;其中,所述高精度三维地震振幅数据体为主频不低于预设主频的三维地震振幅数据体,所述预设主频不低于30Hz(优选为不低于50Hz);
三维地震反演的波阻抗数据体获取模块:用于获取工区的三维地震反演的波阻抗数据体;
输导体的空间展布分析模块:基于所述高精度三维地震振幅数据体和/或所述工区的三维地震反演的波阻抗数据体确定输导体的空间展布;
输导体系类型的三维空间分布分析模块;其中,所述输导体系类型的三维空间分布分析模块包括多种地震属性门槛值获取子模块和输导体类型划分子模块;所述多种地震属性门槛值获取子模块用于获取工区不同类型的输导体对应的多种地震属性门槛值;所述输导体类型划分子模块用于基于所述高精度三维地震振幅数据体,利用所述不同类型的输导体对应的多种地震属性门槛值,在确定的输导体的空间展布范围内进行输导体类型划分,得到工区输导体系类型的三维空间分布;其中,所述多种地震属性至少包括振幅变化率、主频和频率衰减因子;
输导体的宏观输导性能分析模块;其中,输导体的宏观输导性能分析模块包括物性计算公式获取子模块和输导体物性计算子模块;所述物性计算公式获取子模块用于获取不同类型的输导体对应的物性计算公式,其中,物性计算公式为物性与波阻抗的关系式,物性包括孔隙度和渗透率中的至少一种;所述输导体物性计算子模块用于基于所述三维地震反演的波阻抗数据体,利用所述工区输导体系类型的三维空间分布计算工区输导体的物性,从而实现工区输导体的宏观输导性能表征。
11.根据权利要求10所述的系统,其中,所述高精度三维地震振幅数据体获取模块包括:
第一获取单元:用于获取工区的三维地震振幅数据体;
第一判别单元:判断第一获取单元获取的所述工区的三维地震振幅数据体的主频是否低于所述预设主频;若第一判别单元判断得出所述工区的三维地震振幅数据体的主频不低于预设主频,则采用第一处理单元进行处理;若所述工区的三维地震振幅数据体的主频低于预设主频,则采用第二处理单元进行处理;
第一处理单元:用于实施将所述工区的三维地震振幅数据体作为所述工区的高精度三维地震振幅数据体的操作;
第二处理单元:用于实施对所述工区的三维地震振幅数据体进行三维地震数据的拓频处理产生主频不低于所述预设主频的所述工区的高精度三维地震振幅数据体的操作。
12.根据权利要求10所述的系统,其中,所述三维地震反演的波阻抗数据体获取模块用于基于获取的高精度三维地震振幅数据体进行反演得到所述工区的三维地震反演的波阻抗数据体。
13.根据权利要求10所述的系统,其中,所述输导体的空间展布分析模块包括:
第三处理单元:用于基于所述三维地震反演的波阻抗数据体定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布;
第四处理单元:用于基于高精度三维地震振幅数据体识别输导体系横向展布,对初步得到的三维输导体系的空间展布中的横向展布进行约束,实现输导体的空间展布的最终确定。
14.根据权利要求13所述的系统,其中,所述第三处理单元用于利用三维地震反演的波阻抗数据体的反演数据值信息,通过门槛值划分,初步定量预测三维输导体系包括厚度在内的空间展布,初步得到三维输导体系的空间展布。
15.根据权利要求13所述的系统,其中,所述串珠体空间发育规律分析模块进一步包括:
第七处理单元:用于基于第六处理单元确定的串珠体的优势发育层位,进一步确定优势发育层中的串珠体的顶部发育深度、底部发育深度、东西向发育宽度、南北向发育宽度、最大振幅变化率数值、钻遇的串珠体产油量、钻遇的串珠体产气量、钻遇的串珠体产水量和钻遇的串珠体的油气水地球化学参数中的至少一种。
16.根据权利要求10所述的系统,其中,输导体的类型包括溶洞、溶孔、裂缝和断裂。
17.根据权利要求10所述的系统,其中,所述多种地震属性门槛值获取子模块包括:
第二获取单元:用于获取输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料;
第五处理单元:用于采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,在所述已钻井钻遇的输导体中确定出厚度超过所述高精度三维地震振幅数据体调谐厚度的输导体作为第一输导体,并划分第一输导体的输导体类型;
第六处理单元:用于利用时深转换的技术方法,在所述高精度三维地震振幅数据体中标定所述第一输导体的位置;基于标定的所述第一输导体的位置,确定出第一输导体中不同类型的输导体的所述多种地震属性门槛值。
18.根据权利要求10所述的系统,其中,所述物性计算公式获取子模块包括:
第三获取单元:用于获取输导体的空间展布范围内已钻井的岩心物性测试结果
第七处理单元:用于采用输导体的空间展布范围内已钻井的岩心和/或岩屑地质资料,划分已钻井钻遇的输导体的输导体类型;
第八处理单元:用于按照输导体类型,分别进行物性与波阻抗关系拟合,得到不同类型的输导体对应的物性计算公式。
19.一种碳酸盐岩输导体系定量评价装置,包括处理器及存储器;其中,
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现权利要求1-9中任一项所述的碳酸盐岩输导体系定量评价方法的步骤。
20.一种计算机可读存储介质,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现权利要求1-9中任一项所述的碳酸盐岩输导体系定量评价方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010348654.9A CN111596351B (zh) | 2020-04-28 | 2020-04-28 | 碳酸盐岩输导体系定量评价方法、系统、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010348654.9A CN111596351B (zh) | 2020-04-28 | 2020-04-28 | 碳酸盐岩输导体系定量评价方法、系统、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111596351A true CN111596351A (zh) | 2020-08-28 |
CN111596351B CN111596351B (zh) | 2023-04-25 |
Family
ID=72180441
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010348654.9A Active CN111596351B (zh) | 2020-04-28 | 2020-04-28 | 碳酸盐岩输导体系定量评价方法、系统、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111596351B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112698397A (zh) * | 2020-11-20 | 2021-04-23 | 中国石油天然气股份有限公司 | 一种台盆区走滑断裂空腔储层刻画方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090132169A1 (en) * | 2007-11-19 | 2009-05-21 | Schlumberger Technology Corporation | Methods and systems for evaluating fluid movement related reservoir properties via correlation of low-frequency part of seismic data with borehole measurements |
US20120195165A1 (en) * | 2011-01-31 | 2012-08-02 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
CN103527184A (zh) * | 2013-10-28 | 2014-01-22 | 北京大学 | 一种白云岩储层的预测方法和系统 |
CN103529475A (zh) * | 2013-04-19 | 2014-01-22 | 中国石油大学(华东) | 一种识别和解释碳酸盐岩古岩溶储层三维结构的方法 |
CN104818982A (zh) * | 2015-04-20 | 2015-08-05 | 中国石油天然气股份有限公司 | 一种缝洞型储层连通性定量分析方法及装置 |
CN104834007A (zh) * | 2015-05-04 | 2015-08-12 | 中国石油天然气股份有限公司 | 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 |
CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
CN105510991A (zh) * | 2015-12-01 | 2016-04-20 | 中国石油天然气股份有限公司 | 碳酸盐岩油气勘探方法和装置 |
CN109959965A (zh) * | 2017-12-25 | 2019-07-02 | 中国石油天然气股份有限公司 | 碳酸盐岩缝洞体内部结构的获取方法及装置 |
CN110308487A (zh) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | 一种断溶体型油藏定量表征方法 |
-
2020
- 2020-04-28 CN CN202010348654.9A patent/CN111596351B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090132169A1 (en) * | 2007-11-19 | 2009-05-21 | Schlumberger Technology Corporation | Methods and systems for evaluating fluid movement related reservoir properties via correlation of low-frequency part of seismic data with borehole measurements |
US20120195165A1 (en) * | 2011-01-31 | 2012-08-02 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
CN103529475A (zh) * | 2013-04-19 | 2014-01-22 | 中国石油大学(华东) | 一种识别和解释碳酸盐岩古岩溶储层三维结构的方法 |
CN103527184A (zh) * | 2013-10-28 | 2014-01-22 | 北京大学 | 一种白云岩储层的预测方法和系统 |
CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
CN104818982A (zh) * | 2015-04-20 | 2015-08-05 | 中国石油天然气股份有限公司 | 一种缝洞型储层连通性定量分析方法及装置 |
CN104834007A (zh) * | 2015-05-04 | 2015-08-12 | 中国石油天然气股份有限公司 | 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 |
CN105510991A (zh) * | 2015-12-01 | 2016-04-20 | 中国石油天然气股份有限公司 | 碳酸盐岩油气勘探方法和装置 |
CN109959965A (zh) * | 2017-12-25 | 2019-07-02 | 中国石油天然气股份有限公司 | 碳酸盐岩缝洞体内部结构的获取方法及装置 |
CN110308487A (zh) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | 一种断溶体型油藏定量表征方法 |
Non-Patent Citations (4)
Title |
---|
ZHIYAO ZHANG等: "Variations of diamondoids distributions in petroleum fluids during migration induced phase fractionation: A case study from the Tazhong area,NW China", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
余智超等: "塔河油田缝洞型油藏不同成因岩溶储集体表征", 《油气地质与采收率》 * |
朱光有等: "塔里木盆地8000m以深超深层海相油气勘探潜力与方向", 《天然气地球科学》 * |
李宗杰等: "塔北超深层碳酸盐岩储层预测方法和技术", 《石油与天然气地质》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112698397A (zh) * | 2020-11-20 | 2021-04-23 | 中国石油天然气股份有限公司 | 一种台盆区走滑断裂空腔储层刻画方法 |
CN112698397B (zh) * | 2020-11-20 | 2023-08-22 | 中国石油天然气股份有限公司 | 一种台盆区走滑断裂空腔储层刻画方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111596351B (zh) | 2023-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Qi et al. | Attribute expression of fault-controlled karst—Fort Worth Basin, Texas: A tutorial | |
CN103527184B (zh) | 一种白云岩储层的预测方法和系统 | |
CN105528522A (zh) | 一种基于空间网格计算连续性油气藏资源量的方法及装置 | |
Milad et al. | Impact of lithofacies variations and structural changes on natural fracture distributions | |
CN110632652B (zh) | 断块油气藏的油气运移方式确定方法及装置 | |
CN105388525A (zh) | 储层预测方法及装置 | |
Infante-Paez et al. | Seismic analysis of paleotopography and stratigraphic controls on total organic carbon: Rich sweet spot distribution in the Woodford Shale, Oklahoma, USA | |
CN104834006B (zh) | 结构成因沉积参数图的处理方法和装置 | |
Schneider et al. | Interpretation of fractured zones using seismic attributes—Case study from Teapot Dome, Wyoming, USA | |
CN111596351B (zh) | 碳酸盐岩输导体系定量评价方法、系统、装置及存储介质 | |
CN110910267B (zh) | 海相碳酸盐岩隔夹层的确定方法及装置 | |
CN106568918B (zh) | 页岩有机碳含量toc预测方法 | |
da Fonseca Aguiar et al. | Analysis of seismic attributes to recognize bottom simulating reflectors in the Foz do Amazonas basin, Northern Brazil | |
Patel et al. | Seismic illumination of small-throw seismogenic faults, Anadarko Basin, Oklahoma | |
Naseer | Application of instantaneous spectral decomposition-based porosity simulations for imaging shallow-marine stratigraphic traps of Lower-Eocene carbonates sequences of Indus Basin, Onshore Pakistan | |
CN113484907B (zh) | 一种预测不同类型储层平面上分布的方法 | |
Deng et al. | Machine learning based stereoscopic triple sweet spot evaluation method for shale reservoirs | |
CN112505761A (zh) | 一种储层含气性检测方法和装置 | |
Fournier et al. | Quantitative lithostratigraphic interpretation of seismic data for characterization of the Unayzah Formation in central Saudi Arabia | |
RU22830U1 (ru) | Автоматизированная система поиска, разведки и оценки эксплуатационных свойств залежей и месторождений полезных ископаемых и прогноза тектонических и литофизических свойств геологических сред | |
CN114755740A (zh) | 岩石分布确定方法、装置、设备及介质 | |
CN107435537B (zh) | 一种录井目标层自动划分的方法及装置 | |
RU2206911C2 (ru) | Способ поиска, разведки, исследования и создания модели месторождения полезных ископаемых | |
CN117345208B (zh) | 一种压裂优势区定量表征方法、装置、电子设备及介质 | |
CN113031056B (zh) | 构造约束下断层封闭性分析方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |