CN112379435A - 相控岩溶型缝洞集合体刻画方法及装置 - Google Patents
相控岩溶型缝洞集合体刻画方法及装置 Download PDFInfo
- Publication number
- CN112379435A CN112379435A CN202011197527.XA CN202011197527A CN112379435A CN 112379435 A CN112379435 A CN 112379435A CN 202011197527 A CN202011197527 A CN 202011197527A CN 112379435 A CN112379435 A CN 112379435A
- Authority
- CN
- China
- Prior art keywords
- attribute
- karst
- sedimentary facies
- target reservoir
- 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 62
- 208000035126 Facies Diseases 0.000 claims abstract description 86
- 238000010586 diagram Methods 0.000 claims abstract description 28
- 230000005012 migration Effects 0.000 claims abstract description 27
- 238000013508 migration Methods 0.000 claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 15
- 238000004590 computer program Methods 0.000 claims description 16
- 230000000694 effects Effects 0.000 claims description 15
- 239000011148 porous material Substances 0.000 claims description 14
- 238000009826 distribution Methods 0.000 claims description 13
- 239000010459 dolomite Substances 0.000 claims description 13
- 229910000514 dolomite Inorganic materials 0.000 claims description 13
- 230000009471 action Effects 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 8
- 235000019738 Limestone Nutrition 0.000 claims description 7
- 239000006028 limestone Substances 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000007621 cluster analysis Methods 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 2
- 238000013075 data extraction Methods 0.000 claims 1
- 230000004927 fusion Effects 0.000 abstract description 7
- 230000008569 process Effects 0.000 description 11
- 238000011161 development Methods 0.000 description 10
- 239000011435 rock Substances 0.000 description 10
- 239000007789 gas Substances 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 8
- 238000005457 optimization Methods 0.000 description 8
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 6
- 238000005553 drilling Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000012512 characterization method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000008021 deposition Effects 0.000 description 3
- 235000019994 cava Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 241001385887 Tachys Species 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000001973 epigenetic effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种相控岩溶型缝洞集合体刻画方法及装置，该方法包括：根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。通过相控、分频、重构、融合等地震手段，逐步间接实现了小尺度的相控岩溶型孔洞集合体的精细刻画。
Description
技术领域
本发明涉及石油及天然气勘探开发技术领域，尤其涉及一种相控岩溶型缝洞集合体刻画方法及装置。
背景技术
近些年来，随着计算机技术与地球物理技术的不断进步与完善，碳酸盐岩的油气勘探开发取得了重大突破，已成为油气增储上产最重要的领域之一。与碎屑岩油相比，碳酸盐岩通常具有地层年代古老、现今埋藏深、基质孔隙度低，岩溶改造期次多等特点。这些特点的存在为碳酸盐岩储层地震描述，特别是非均质性极强的岩溶孔洞型储层的刻画来了巨大的挑战。
目前，关于碳酸盐岩岩溶孔洞体的研究主要集中于灰岩断溶型储层，比较典型的如塔里木盆地塔河油田奥陶系。这类储层在形成过程中风化暴露时间长，岩溶孔洞单元的规模主要为米级。由数个至数十个孔洞单元所组成的孔洞集合体在地震剖面上表现出十分明显的“串珠状”或“羊肉串状”特征。地震预测过程中，通常首先开展OVT域深度偏移处理，提高孔洞单元成像与归位精度；然后基于孔洞体发育时地震振幅减弱的特点，采用常规振幅类属性体实现孔洞体预测。这套技术方法可以很好的指导大尺度灰岩岩溶气藏的勘探开发，但对于受原始沉积相影响更大的相控岩溶型小尺度孔洞体则往往预测结果较差。首先，相控岩溶型孔洞集合体的形成具有明显的沉积选择性，其分布除了与风化暴露期古地貌形态、古断裂有关外，还明显受到原始沉积相带展布的控制。通常水动力强、沉积坡度较缓的地区，所形成的沉积物规模更大、原始物性更好，后期岩溶改造最有利，对于这一因素现有的孔洞体预测方法并未考虑，导致相控岩溶型孔洞体的预测能力差。其次，相控岩溶型储层的风化暴露时间通常较短或者仅受到表生岩溶作用改造，这使得所形成的单个孔洞单元规模很小，主要表现为厘米级。在叠合地震数据上利用各类常规属性与反演手段预测此类缝洞单元是无法实现的，因而无法对相控岩溶型缝洞，特别是小尺度的缝洞进行刻画。
发明内容
本发明实施例提供一种相控岩溶型缝洞集合体刻画方法，用以对相控岩溶型小尺度缝洞进行刻画，有效提升对相控岩溶孔洞体的预测能力，该方法包括：
根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；
利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；
以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；
根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；
在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；
利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。
本发明实施例还提供一种相控岩溶型缝洞集合体刻画装置，用以对相控岩溶型小尺度缝洞进行刻画，有效提升对相控岩溶孔洞体的预测能力，该装置包括：
沉积相平面图预测模块，用于根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；
地震数据分频模块，用于利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；
地震数据选取模块，用于以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；
地震数据重构模块，用于根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；
地震属性提取模块，用于在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；
三维刻画模块，用于利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述相控岩溶型缝洞集合体刻画方法。
本发明实施例也提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述相控岩溶型缝洞集合体刻画方法的计算机程序。
本发明实施例中，通过根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。通过以地震波形聚类得到的沉积相平面图为约束，开展地震数据分频、优选与重构，在重构数据体上提取振幅差异属性与纹理方差属性，分别表征孔洞单元的垂向与侧向变化，采用二维色表方式，对振幅差异和纹理方差属性进行融合，实现相控岩溶型孔洞体精细刻画。考虑到地震资料的极限分辨率问题，本发明并未刻意预测每个厘米级的岩溶孔洞单元，而是通过相控、分频、重构、融合等地震手段，逐步间接实现了小尺度的相控岩溶型孔洞集合体的精细刻画。相比已有的地震预测方法，本发明以沉积相平面图为约束，强调了宏观沉积规律对于地震预测过程的控制作用，使得分频地震数据的优选和重构更加合理，提升了地震资料对于缝洞体的响应能力，从而有效提升对相控岩溶型孔洞体的预测能力。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中相控岩溶型缝洞集合体刻画方法的示意图。
图2为本发明具体实施例中步骤101的具体实现方法示意图。
图3为本发明具体实施例中步骤103的具体实现方法示意图。
图4为本发明具体实施例中步骤104的具体实现方法示意图。
图5为本发明具体实施例中步骤105的具体实现方法示意图。
图6为本发明一具体应用实施中的相控岩溶型小尺度缝洞集合体刻画示意图。
图7为本发明一具体应用实施中的基于波形聚类的沉积相带平面预测图。
图8为本发明一具体应用实施中的主频10Hz的分频体上所提取的岩溶储层顶面沿层均方根振幅属性图。
图9为本发明一具体应用实施中的相控岩溶型小尺度孔洞单元三维刻画结果。
图10为本发明一具体应用实施中的相控岩溶型小尺度孔洞单元地震预测过井剖面。
图11为本发明实施例中相控岩溶型缝洞集合体刻画装置的示意图。
图12为本发明具体实施例中沉积相平面图预测模块1101的结构示意图。
图13为本发明具体实施例中地震数据选取模块1103的结构示意图。
图14为本发明具体实施例中地震数据重构模块1104的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例提供了一种相控岩溶型缝洞集合体刻画方法，用以对相控岩溶型小尺度缝洞进行刻画，有效提升对相控岩溶孔洞体的预测能力，如图1所示，该方法包括：
步骤101：根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；
步骤102：利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；
步骤103：以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；
步骤104：根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；
步骤105：在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；
步骤106：利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。
由图1所示流程可以得知，本发明实施例中，通过根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。通过以地震波形聚类得到的沉积相平面图为约束，开展地震数据分频、优选与重构，在重构数据体上提取振幅差异属性与纹理方差属性，分别表征孔洞单元的垂向与侧向变化，采用二维色表方式，对振幅差异和纹理方差属性进行融合，实现相控岩溶型孔洞体精细刻画。考虑到地震资料的极限分辨率问题，本发明并未刻意预测每个厘米级的岩溶孔洞单元，而是通过相控、分频、重构、融合等地震手段，逐步间接实现了小尺度的相控岩溶型孔洞集合体的精细刻画。相比已有的地震预测方法，本发明以沉积相平面图为约束，强调了宏观沉积规律对于地震预测过程的控制作用，使得分频地震数据的优选和重构更加合理，提升了地震资料对于缝洞体的响应能力，从而有效提升对相控岩溶型孔洞体的预测能力。
具体实施时，首先根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图。具体实施时，如图2所示，包括：
步骤201：根据目标储层的测井数据，对目标储层进行岩电标定，确定敏感测井曲线；
步骤202：基于敏感测井曲线，采用Q型聚类分析，预测白云岩、灰岩以及硅质岩在单井上的分布，并通过白云岩厚度与地层厚度的比值划分沉积相类型；
步骤203：针对划分的不同沉积相类型开展井震标定，确定每种沉积相类型对应的地震波形分类；
步骤204：按照有监督分类的方法，对带有分类标签的每种沉积相类型对应的地震波形开展地震波形聚类，实现沉积相类型的平面预测，得到目标储层的沉积相平面图。
得到目标储层的沉积相平面图后，利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体，以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体，根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体。
其中，利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体，具体是指利用小波变换将目标储层的叠前偏移地震数据按照预定的间距分频为一系列的单频段数据体，预定的间距例如可以是10Hz的间距。以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体，具体实施过程如图3所示，包括：
步骤301：对多个分频体的振幅值进行归一化处理，得到归一化后的多个分频体；
步骤302：在归一化后的多个分频体上，分别提取岩溶储层顶部的沿层均方根振幅属性平面图；
步骤303：按照宏观分布规律相似原则，选出一张与目标储层的沉积相平面图最相似的沿层均方根振幅属性图，将所述沿层均方根振幅属性图对应的分频体确定为选取出的低频端数据体。
根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体，具体实施过程如图4所示，包括：
步骤401：根据岩溶孔洞体地层厚度和多个分频体的调谐厚度，在多个分频体中选取出高频端数据体；
步骤402：将低频端数据体与高频端数据体按照等比例原则进行合并，得到重构地震数据体。
具体实施例中，根据目标储层的实钻井岩溶孔洞钻遇结果，统计岩溶孔洞体平均厚度。选取调谐厚度小于孔洞体平均厚度且对应振幅属性图上随机噪声最少的分频体，作为优选出的高频端数据体。
得到重构地震数据体后，在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体。其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化。具体实施过程，如图5所示，包括：
步骤501：在重构数据体上，采用熵值分析法，计算振幅差异属性体集合与纹理方差属性体集合，得到振幅差异属性体集合中每个振幅差异属性体的重要性指标，与纹理方差属性体集合中每个纹理方差属性体的重要性指标；
步骤502：根据振幅差异属性体集合中每个振幅差异属性体的重要性指标，在振幅差异属性体集合中选取振幅差异属性体；
步骤503：根据纹理方差属性体集合中每个纹理方差属性体的重要性指标，在纹理方差属性体集合中选取纹理方差属性体。
具体实施例中，首先采用不同的参数计算振幅差异属性集合
和纹理方差属性
其中，i和j分别表示属性下标索引量，W和Q分别为振幅差异属性和纹理方差属性的参数选择个数，均需要事先由解释人员确定，为了后续讨论方便，我们以振幅差异属性
为例进行讨论，纹理方差属性
的选择方式与此一样。然后根据实钻井孔洞体标定结果，由解释人员将实钻目标孔洞体分成N类，表示为
每一类
是样本点的集合，Nc为每一个类别集合sc的样本个数，(c为属性的下标索引)。每一个xi是实钻井孔处的样本点，由上文提到的W个属性表示，即xi＝(xi,1,xi,2,...,xi,W+Q)。
计算属性之间的相似性度量矩阵，|*|表示取两者之间的L1范数：
任意两个属性Yi,Yj之间的相互熵Hij可以定义为：
考虑所有特征空间之间的特征对象，总体上相互熵EH定义为：
为了计算单个属性Yt的重要程度，将属性集合
中，剔除属性Yt后，重新计算总体上相互熵的EH(F-Yt)，那么属性Yt的重要性
可以表示为：
其中，(F-Yt)表示从属性集合中剔除掉属性Yt，(t为属性的下标索引)。
计算每一个特征的重要性，
并且对他们进行排序，数值大的放在前面。
对每一种类型孔洞体数据集合Sc，都进行计算相对应的EH和属性特征所对应的重要性指标集合
为了容易区分，将Sc对应的指标集合
记录为：
对于任何一个属性Yt，根据其重要性数值
可以得到其在类Sc中的对应重要程度编号
那么属性Yt在整个数据集合中的重要性指标
可以表示为：
(t为属性的下标索引)。
计算所有属性
所对应的
按照由大到小的顺序排序，从中选择出数值
最大的属性Ym(m为最大的属性的下标索引)，这个属性Ym将被选取为用于后续的融合与刻画操作的振幅差异属性体。
按照上述步骤，对纹理方差属性
进行如上相同操作，从中选择出
值最大的属性Tp(p为最大的属性的下标索引)，这个属性Tp将被选取为用于后续的融合与刻画操作的纹理方差属性体。
提取得到振幅差异属性体和纹理方差属性体后，利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。
下面给出一具体实例说明本发明实施例如何进行相控岩溶型缝洞集合体刻画。本例应用于四川盆地高石梯-磨溪地区震旦系灯影组碳酸盐岩气藏的岩溶孔洞预测。该气藏具有典型的构造-岩性气藏特点，储层分布受到白云岩丘滩体与短时间风化岩溶改造共同控制，储集单元为厘米级小尺度孔洞。
本具体实例以沉积相预测结果为约束，开展分频数据体优选与重构，进而将纹理方差属性体与振幅差异属性体融合，完成相控型小尺度孔洞预测体刻画。具体实施步骤如图6所示，包括：
首先基于研究区目的层岩相的古地理背景，开展岩心描述与岩电标定，优选出3条对于矿物成分最为敏感的测井曲线。基于敏感测井曲线，采用Q型聚类分析的方法，预测单井上白云岩、灰岩以及硅质岩分布。
在单井上，根据目的层段内部的白云岩厚度与地层厚度的比值划分单井沉积相类型。
针对不同沉积相类型开展井震标定，确定出一定时窗范围内各自所对应的地震波形分类。按照有监督分类的方法开展波形聚类，实现沉积相类型的平面预测，如图7所示为本具体实例的基于波形聚类的沉积相带平面预测图，其中3类代表的是白云岩厚度与地层厚度比值大于0.6的藻丘相；2类代表的是白云岩厚度与地层厚度比值介于0.3-0.6之间的颗粒滩相，1类代表的是白云岩厚度与地层厚度比值小于0.3的滩间海相。由3类到1类，沉积相原始物性逐渐变差，厚度逐渐减薄，后期被岩溶改造的能力逐渐降低。
采用连续小波变换，将原始地震数据按照10Hz-80Hz(间距为10Hz)8个中心频率进行分解，得到分频体。由于不同中心频率的分频体的能量差别较大，因此将所有分频体的振幅值进行归一化处理。
在归一化后的各个分频体上，分别提取岩溶储层顶部的沿层均方根振幅属性平面图。按照宏观分布规律相似原则，选出一张与沉积相平面图最为接近的振幅属性图，其所对应的分频体即为优选出的低频端数据体。如图8所示为主频10Hz的分频体上所提取的岩溶储层顶面沿层均方根振幅属性图。通过将所有低频端分频体上所提取的沿层均方根振幅图与沉积相图(图7)比对后认为，10Hz的分频体上所提取的图件与沉积相宏观分布规律最为相似，能够最好的放映出原始沉积相带对后期岩溶改造的控制作用。因此，10Hz分频体即为最终优选出的低频端数据体。
根据研究区实钻井岩溶孔洞钻遇结果，统计岩溶孔洞体平均厚度。选取调谐厚度小于孔洞体平均厚度并且对应振幅属性图上随机噪声最少的分频体作为优选出的高频端数据体。
将得到的低频端单频体与高频端单频体按照1:1的等比例关系进行合并，完成地震数据体重构。
在重构的地震数据体上，采用熵值分析法完成振幅差异属性体集合与纹理方差属性体集合的初选与排序。按照自动排序的原则，根据重要性指标从振幅差异属性体集合与纹理方差属性体集合，这两个数据集合中分别优选中最优化的振幅差异属性体与纹理方差属性体。
采用二维色表的方式将最终优选的振幅差异体与纹理方差体进行融合，其中振幅差异体取值域，纹理方差体取亮度，通过二者的值共同确定一个颜色。如图9所示，为本具体实例中相控岩溶型小尺度孔洞单元三维刻画结果。
如图10所示为相控岩溶型小尺度孔洞单元地震预测过井剖面，结合图9和图10能够看出：应用本方法在研究区共识别出3类孔洞体发育区，其中一类发育区孔洞单元最为发育，二类区其次，三类区最差。本次预测结果与该区整体地质认识吻合，与已钻井匹配程度高，大幅度提高了孔洞体预测精度。且利用该成果有效指导了研究区的开发建产，共提出建议井位目标28口，其中15口已实施钻探并全部获得日产50万方以上的高产气流，为该地区的天然气快速上产提供了地质依据。
通过围绕沉积相带的地震预测、地震数据的分频与重构、敏感地震属性的优选与融合等方面进行研究，可以有效提升现有地震资料对于小尺度岩溶孔洞体的预测能力，为此类油气田高效勘探开发提供可靠数据。相比已有的地震预测方法，本具体实例强调了宏观沉积规律对于地震预测过程的控制作用，使得分频地震数据的优选和重构更加合理，提升了地震资料对于缝洞体的响应能力。针对小尺度岩溶孔洞体，采用敏感地震属性融合的方式，实现了精细雕刻。
上述具体应用的实施仅为举例，其余实施方式不再一一赘述。
基于同一发明构思，本发明实施例还提供一种相控岩溶型缝洞集合体刻画装置，由于相控岩溶型缝洞集合体刻画装置所解决问题的原理与相控岩溶型缝洞集合体刻画方法相似，因此相控岩溶型缝洞集合体刻画装置的实施可以参见相控岩溶型缝洞集合体刻画方法的实施，重复之处不再赘述，具体结构如图11所示：
沉积相平面图预测模块1101，用于根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；
地震数据分频模块1102，用于利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；
地震数据选取模块1103，用于以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；
地震数据重构模块1104，用于根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；
地震属性提取模块1105，用于在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；
三维刻画模块1106，用于利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。
具体实施例中，沉积相平面图预测模块1101，如图12所示，包括：
敏感测井曲线确定单元1201，用于根据目标储层的测井数据，对目标储层进行岩电标定，确定敏感测井曲线；
测井解释单元1202，用于基于敏感测井曲线，采用Q型聚类分析，预测白云岩、灰岩以及硅质岩在单井上的分布，并通过白云岩厚度与地层厚度的比值划分沉积相类型；
地震波形分类分类单元1203，用于针对划分的不同沉积相类型开展井震标定，确定每种沉积相类型对应的地震波形分类；
沉积相平面图预测单元1204，用于按照有监督分类的方法，对带有分类标签的每种沉积相类型对应的地震波形开展地震波形聚类，实现沉积相类型的平面预测，得到目标储层的沉积相平面图。
具体实施例中，地震数据选取模块1103，如图13所示，包括：
归一处理单元1301，用于对多个分频体的振幅值进行归一化处理，得到归一化后的多个分频体；
振幅属性平面图提取单元1302，用于在归一化后的多个分频体上，分别提取岩溶储层顶部的沿层均方根振幅属性平面图；
低频端数据体选取单元1303，用于按照宏观分布规律相似原则，选出一张与目标储层的沉积相平面图最相似的沿层均方根振幅属性图，将所述沿层均方根振幅属性图对应的分频体确定为选取出的低频端数据体。
具体实施例中，地震数据重构模块1104，如图14所示，包括：
高频端数据体选取单元1401，用于根据岩溶孔洞体地层厚度和多个分频体的调谐厚度，在多个分频体中选取出高频端数据体；
地震数据重构单元1402，用于将低频端数据体与高频端数据体按照等比例原则进行合并，得到重构地震数据体。
具体实施例中，地震属性提取模块1105具体用于：
在重构数据体上，采用熵值分析法，计算振幅差异属性体集合与纹理方差属性体集合，得到振幅差异属性体集合中每个振幅差异属性体的重要性指标，与纹理方差属性体集合中每个纹理方差属性体的重要性指标；
根据振幅差异属性体集合中每个振幅差异属性体的重要性指标，在振幅差异属性体集合中选取振幅差异属性体；
根据纹理方差属性体集合中每个纹理方差属性体的重要性指标，在纹理方差属性体集合中选取纹理方差属性体。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述相控岩溶型缝洞集合体刻画方法。
本发明实施例还提供一种计算机可读存储介质，存储有执行上述相控岩溶型缝洞集合体刻画方法的计算机程序。
综上所述，本发明实施例提供的相控岩溶型缝洞集合体刻画方法及装置具有如下优点：
通过根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。通过以地震波形聚类得到的沉积相平面图为约束，开展地震数据分频、优选与重构，在重构数据体上提取振幅差异属性与纹理方差属性，分别表征孔洞单元的垂向与侧向变化，采用二维色表方式，对振幅差异和纹理方差属性进行融合，实现相控岩溶型孔洞体精细刻画。考虑到地震资料的极限分辨率问题，本发明并未刻意预测每个厘米级的岩溶孔洞单元，而是通过相控、分频、重构、融合等地震手段，逐步间接实现了小尺度的相控岩溶型孔洞集合体的精细刻画。相比已有的地震预测方法，本发明以沉积相平面图为约束，强调了宏观沉积规律对于地震预测过程的控制作用，使得分频地震数据的优选和重构更加合理，提升了地震资料对于缝洞体的响应能力，从而有效提升对相控岩溶型孔洞体的预测能力。
本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、装置和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (12)
1.一种相控岩溶型缝洞集合体刻画方法，其特征在于，包括：
根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；
利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；
以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；
根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；
在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；
利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。
2.如权利要求1所述的方法，其特征在于，根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图，包括：
根据目标储层的测井数据，对目标储层进行岩电标定，确定敏感测井曲线；
基于敏感测井曲线，采用Q型聚类分析，预测白云岩、灰岩以及硅质岩在单井上的分布，并通过白云岩厚度与地层厚度的比值划分沉积相类型；
针对划分的不同沉积相类型开展井震标定，确定每种沉积相类型对应的地震波形分类；
按照有监督分类的方法，对带有分类标签的每种沉积相类型对应的地震波形开展地震波形聚类，实现沉积相类型的平面预测，得到目标储层的沉积相平面图。
3.如权利要求1所述的方法，其特征在于，以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体，包括：
对多个分频体的振幅值进行归一化处理，得到归一化后的多个分频体；
在归一化后的多个分频体上，分别提取岩溶储层顶部的沿层均方根振幅属性平面图；
按照宏观分布规律相似原则，选出一张与目标储层的沉积相平面图最相似的沿层均方根振幅属性图，将所述沿层均方根振幅属性图对应的分频体确定为选取出的低频端数据体。
4.如权利要求1所述的方法，其特征在于，根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体，包括：
根据岩溶孔洞体地层厚度和多个分频体的调谐厚度，在多个分频体中选取出高频端数据体；
将低频端数据体与高频端数据体按照等比例原则进行合并，得到重构地震数据体。
5.如权利要求1所述的方法，其特征在于，在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体，包括：
在重构数据体上，采用熵值分析法，计算振幅差异属性体集合与纹理方差属性体集合，得到振幅差异属性体集合中每个振幅差异属性体的重要性指标，与纹理方差属性体集合中每个纹理方差属性体的重要性指标；
根据振幅差异属性体集合中每个振幅差异属性体的重要性指标，在振幅差异属性体集合中选取振幅差异属性体；
根据纹理方差属性体集合中每个纹理方差属性体的重要性指标，在纹理方差属性体集合中选取纹理方差属性体。
6.一种相控岩溶型缝洞集合体刻画装置，其特征在于，包括：
沉积相平面图预测模块，用于根据目标储层的测井数据，进行沉积相类型的平面预测，得到目标储层的沉积相平面图；
地震数据分频模块，用于利用小波变换将目标储层的叠前偏移地震数据进行分解，得到多个分频体；
地震数据选取模块，用于以目标储层的沉积相平面图为约束，在多个分频体中选取出低频端数据体；
地震数据重构模块，用于根据低频端数据体，重构叠前偏移地震数据，得到重构地震数据体；
地震属性提取模块，用于在重构地震数据体上提取得到振幅差异属性体和纹理方差属性体；其中，振幅差异属性体用于表征顺层岩溶作用下的孔洞体平面变化；纹理方差属性体用于表征垂向断裂岩溶作用下的孔洞体垂向变化；
三维刻画模块，用于利用二维色表融合振幅差异属性体和纹理方差属性体，得到相控岩溶型缝洞集合体的三维刻画。
7.如权利要求6所述的装置，其特征在于，所述沉积相平面图预测模块，包括：
敏感测井曲线确定单元，用于根据目标储层的测井数据，对目标储层进行岩电标定，确定敏感测井曲线；
测井解释单元，用于基于敏感测井曲线，采用Q型聚类分析，预测白云岩、灰岩以及硅质岩在单井上的分布，并通过白云岩厚度与地层厚度的比值划分沉积相类型；
地震波形分类分类单元，用于针对划分的不同沉积相类型开展井震标定，确定每种沉积相类型对应的地震波形分类；
沉积相平面图预测单元，用于按照有监督分类的方法，对带有分类标签的每种沉积相类型对应的地震波形开展地震波形聚类，实现沉积相类型的平面预测，得到目标储层的沉积相平面图。
8.如权利要求6所述的装置，其特征在于，所述地震数据选取模块，包括：
归一处理单元，用于对多个分频体的振幅值进行归一化处理，得到归一化后的多个分频体；
振幅属性平面图提取单元，用于在归一化后的多个分频体上，分别提取岩溶储层顶部的沿层均方根振幅属性平面图；
低频端数据体选取单元，用于按照宏观分布规律相似原则，选出一张与目标储层的沉积相平面图最相似的沿层均方根振幅属性图，将所述沿层均方根振幅属性图对应的分频体确定为选取出的低频端数据体。
9.如权利要求6所述的装置，其特征在于，所述地震数据重构模块，包括：
高频端数据体选取单元，用于根据岩溶孔洞体地层厚度和多个分频体的调谐厚度，在多个分频体中选取出高频端数据体；
地震数据重构单元，用于将低频端数据体与高频端数据体按照等比例原则进行合并，得到重构地震数据体。
10.如权利要求6所述的装置，其特征在于，所述地震属性提取模块具体用于：
在重构数据体上，采用熵值分析法，计算振幅差异属性体集合与纹理方差属性体集合，得到振幅差异属性体集合中每个振幅差异属性体的重要性指标，与纹理方差属性体集合中每个纹理方差属性体的重要性指标；
根据振幅差异属性体集合中每个振幅差异属性体的重要性指标，在振幅差异属性体集合中选取振幅差异属性体；
根据纹理方差属性体集合中每个纹理方差属性体的重要性指标，在纹理方差属性体集合中选取纹理方差属性体。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。
12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011197527.XA CN112379435B (zh) | 2020-10-30 | 2020-10-30 | 相控岩溶型缝洞集合体刻画方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011197527.XA CN112379435B (zh) | 2020-10-30 | 2020-10-30 | 相控岩溶型缝洞集合体刻画方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112379435A true CN112379435A (zh) | 2021-02-19 |
| CN112379435B CN112379435B (zh) | 2024-06-25 |
Family
ID=74576931
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011197527.XA Active CN112379435B (zh) | 2020-10-30 | 2020-10-30 | 相控岩溶型缝洞集合体刻画方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112379435B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115877460A (zh) * | 2023-02-28 | 2023-03-31 | 福瑞升(成都)科技有限公司 | 一种碳酸盐岩岩溶缝洞型储层的增强方法 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103527184A (zh) * | 2013-10-28 | 2014-01-22 | 北京大学 | 一种白云岩储层的预测方法和系统 |
| WO2014169499A1 (zh) * | 2013-04-19 | 2014-10-23 | 中国石油大学(华东) | 一种识别和解释碳酸盐岩古岩溶储层三维结构的方法 |
| CN104155689A (zh) * | 2014-08-27 | 2014-11-19 | 中国石油集团东方地球物理勘探有限责任公司 | 碳酸盐岩缝洞雕刻方法及装置 |
| CN105158795A (zh) * | 2015-08-27 | 2015-12-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用地层叠前纹理属性值来检测缝洞的方法 |
| US20160363691A1 (en) * | 2015-06-15 | 2016-12-15 | Petrochina Company Limited | Physical simulation method and experiment device of fracture-cavity carbonate reservoir hydrocarbon charge |
| CN107728230A (zh) * | 2017-09-11 | 2018-02-23 | 中国石油大港油田勘探开发研究院 | 一种基于地球物理三联技术的岩性油藏预测方法 |
| CN108196302A (zh) * | 2017-11-28 | 2018-06-22 | 中国石油天然气股份有限公司 | 白云岩缝洞发育储层甜点区预测方法及装置 |
| CN109387873A (zh) * | 2017-08-04 | 2019-02-26 | 中国石油化工股份有限公司 | 一种缝洞储集体反演方法及系统 |
| RU2017131249A (ru) * | 2016-09-07 | 2019-03-05 | Чайна Петролеум Энд Кемикал Корпорейшн | Способ и система автоматического распознавания центра залежи в карстовой пещере |
| CN109779624A (zh) * | 2019-03-11 | 2019-05-21 | 西安石油大学 | 基于岩溶参数综合评价裂缝—孔洞型碳酸盐岩储层的方法 |
| CN110308487A (zh) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | 一种断溶体型油藏定量表征方法 |
| CN110308497A (zh) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | 缝洞型油藏不同类型储集体模型融合方法 |
| CN110309518A (zh) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | 缝洞型碳酸盐岩油藏溶蚀孔洞分级分类建模方法 |
-
2020
- 2020-10-30 CN CN202011197527.XA patent/CN112379435B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014169499A1 (zh) * | 2013-04-19 | 2014-10-23 | 中国石油大学(华东) | 一种识别和解释碳酸盐岩古岩溶储层三维结构的方法 |
| CN103527184A (zh) * | 2013-10-28 | 2014-01-22 | 北京大学 | 一种白云岩储层的预测方法和系统 |
| CN104155689A (zh) * | 2014-08-27 | 2014-11-19 | 中国石油集团东方地球物理勘探有限责任公司 | 碳酸盐岩缝洞雕刻方法及装置 |
| US20160363691A1 (en) * | 2015-06-15 | 2016-12-15 | Petrochina Company Limited | Physical simulation method and experiment device of fracture-cavity carbonate reservoir hydrocarbon charge |
| CN105158795A (zh) * | 2015-08-27 | 2015-12-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用地层叠前纹理属性值来检测缝洞的方法 |
| RU2017131249A (ru) * | 2016-09-07 | 2019-03-05 | Чайна Петролеум Энд Кемикал Корпорейшн | Способ и система автоматического распознавания центра залежи в карстовой пещере |
| CN109387873A (zh) * | 2017-08-04 | 2019-02-26 | 中国石油化工股份有限公司 | 一种缝洞储集体反演方法及系统 |
| CN107728230A (zh) * | 2017-09-11 | 2018-02-23 | 中国石油大港油田勘探开发研究院 | 一种基于地球物理三联技术的岩性油藏预测方法 |
| CN108196302A (zh) * | 2017-11-28 | 2018-06-22 | 中国石油天然气股份有限公司 | 白云岩缝洞发育储层甜点区预测方法及装置 |
| CN110308487A (zh) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | 一种断溶体型油藏定量表征方法 |
| CN110308497A (zh) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | 缝洞型油藏不同类型储集体模型融合方法 |
| CN110309518A (zh) * | 2018-03-20 | 2019-10-08 | 中国石油化工股份有限公司 | 缝洞型碳酸盐岩油藏溶蚀孔洞分级分类建模方法 |
| CN109779624A (zh) * | 2019-03-11 | 2019-05-21 | 西安石油大学 | 基于岩溶参数综合评价裂缝—孔洞型碳酸盐岩储层的方法 |
Non-Patent Citations (10)
| Title |
|---|
| 吕心瑞 等: "《碳酸盐岩储层多尺度缝洞体分类表征――以塔河油田S80单元奥陶系油藏为例》", 《石油与天然气地质》, vol. 38, no. 4, 31 August 2017 (2017-08-31), pages 813 - 821 * |
| 商晓飞 等: "《断控岩溶主控的缝洞型碳酸盐岩内部溶蚀相带表征――以塔河油田10区奥陶系油藏为例》", 《石油学报》, vol. 41, no. 3, 31 March 2020 (2020-03-31), pages 329 - 341 * |
| 孙振涛: "《基于叠前分频振幅差异的溶洞识别技术及应用》", 《石油物探》, vol. 57, no. 3, 31 May 2018 (2018-05-31), pages 452 - 457 * |
| 常少英 等: "《海相碳酸盐岩储层地震预测技术进展及应用实效》", 《海相油气地质》, vol. 25, no. 1, 31 March 2020 (2020-03-31), pages 22 - 34 * |
| 徐红霞 等: "《多属性分析技术在碳酸盐岩断溶体 预测中的应用》", 石油地球物理勘探, vol. 52, no. 2, 31 December 2017 (2017-12-31), pages 158 - 163 * |
| 朱正平 等;: "《塔里木盆地轮古西地区奥陶系碳酸盐岩储层的地震识别》", 《海洋地质前沿》, vol. 29, no. 1, 31 January 2013 (2013-01-31), pages 36 - 42 * |
| 李海英 等: "《顺北地区走滑断裂与断溶体圈闭识别描述技术》", 《中国石油勘探》, vol. 25, no. 3, 31 May 2020 (2020-05-31), pages 107 - 120 * |
| 王红伟 等: "《鄂尔多斯盆地西北部地区奥陶系岩溶缝洞型储层发育特征及有利目标区分析》", 《现代地质》, vol. 25, no. 5, 31 October 2011 (2011-10-31), pages 917 - 924 * |
| 高利君 等: "《塔里木盆地深层岩溶缝洞型储层三维雕刻"五步法"定量描述技术研究与应用》", 《物探与化探》, vol. 44, no. 3, 30 June 2020 (2020-06-30), pages 691 - 697 * |
| 龙隆 等: "《基于OVT域数据各向异性处理的缝洞预测方法――以四川盆地云锦向斜茅口组为例》", 《天然气勘探与开发》, vol. 43, no. 2, 30 June 2020 (2020-06-30), pages 23 - 29 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115877460A (zh) * | 2023-02-28 | 2023-03-31 | 福瑞升(成都)科技有限公司 | 一种碳酸盐岩岩溶缝洞型储层的增强方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112379435B (zh) | 2024-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109061765B (zh) | 非均质薄砂岩互层油藏的圈闭评价方法 | |
| Roy et al. | Generative topographic mapping for seismic facies estimation of a carbonate wash, Veracruz Basin, southern Mexico | |
| CN105652323B (zh) | 一种储层预测方法 | |
| Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
| CN102918423A (zh) | 用于地震烃体系分析的方法 | |
| Gupta et al. | Rock typing in eagle ford, barnett, and woodford formations | |
| Zhu et al. | Rapid identification of high-quality marine shale gas reservoirs based on the oversampling method and random forest algorithm | |
| CN102736107A (zh) | 一种能量约束非均质储层厚度识别系统 | |
| CN103792573A (zh) | 一种基于频谱融合的地震波阻抗反演方法 | |
| Euzen et al. | Well log cluster analysis: an innovative tool for unconventional exploration | |
| CN105842733A (zh) | 页岩储层地震识别方法 | |
| CN105301647B (zh) | 区分灰质泥岩和砂岩的方法 | |
| CN105116449A (zh) | 一种弱反射储层的识别方法 | |
| CN110412662A (zh) | 基于地震多属性深度学习的薄互层储层预测方法 | |
| CN108226997A (zh) | 一种基于叠前地震数据的地震相划分方法 | |
| Iltaf et al. | Facies and petrophysical modeling of Triassic Chang 6 tight sandstone reservoir, Heshui oil field, Ordos basin, China | |
| Wethington et al. | Stratigraphic architecture of the Mississippian limestone through integrated electrofacies classification, Hardtner field area, Kansas and Oklahoma | |
| CN112711068B (zh) | 一种砂岩中油气有效储层预测方法及装置 | |
| CN112379435B (zh) | 相控岩溶型缝洞集合体刻画方法及装置 | |
| Hosseini et al. | A Meta attribute for reservoir permeability classification using well logs and 3D seismic data with probabilistic neural network. | |
| CN113393335A (zh) | 基于多地震属性优化的储层油气预测方法及装置 | |
| Al-Baldawi | Applying the cluster analysis technique in logfacies determination for Mishrif Formation, Amara oil field, South Eastern Iraq | |
| CN114482995B (zh) | 一种细粒沉积物泥质含量精细确定方法 | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| CN109581489B (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 |