CN111090117B - 一种相控正演约束下的有效储层预测方法及系统 - Google Patents
一种相控正演约束下的有效储层预测方法及系统 Download PDFInfo
- Publication number
- CN111090117B CN111090117B CN201811244512.7A CN201811244512A CN111090117B CN 111090117 B CN111090117 B CN 111090117B CN 201811244512 A CN201811244512 A CN 201811244512A CN 111090117 B CN111090117 B CN 111090117B
- Authority
- CN
- China
- Prior art keywords
- frequency
- seismic
- minimum
- area
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000002349 favourable effect Effects 0.000 claims abstract description 47
- 239000012530 fluid Substances 0.000 claims abstract description 18
- 230000008021 deposition Effects 0.000 claims abstract description 13
- 230000008859 change Effects 0.000 claims abstract description 12
- 239000004576 sand Substances 0.000 claims description 36
- 239000011229 interlayer Substances 0.000 claims description 26
- 238000005553 drilling Methods 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 15
- 230000009286 beneficial effect Effects 0.000 claims description 12
- 230000002159 abnormal effect Effects 0.000 claims description 10
- 238000013508 migration Methods 0.000 claims description 10
- 230000005012 migration Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 238000002360 preparation method Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000011946 reduction process Methods 0.000 claims 1
- 125000006850 spacer group Chemical group 0.000 claims 1
- 238000001228 spectrum Methods 0.000 abstract description 10
- 208000035126 Facies Diseases 0.000 abstract description 6
- 230000004044 response Effects 0.000 abstract description 3
- 238000013459 approach Methods 0.000 abstract description 2
- 239000012071 phase Substances 0.000 description 16
- 229930195733 hydrocarbon Natural products 0.000 description 11
- 150000002430 hydrocarbons Chemical class 0.000 description 11
- 239000004215 Carbon black (E152) Substances 0.000 description 8
- 238000001514 detection method Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 5
- 239000010410 layer Substances 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000033764 rhythmic process Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000002637 fluid replacement therapy Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000007790 solid phase Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000012800 visualization 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/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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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
技术领域
本发明涉及一种相控正演约束下的有效储层预测方法及系统,属于地震勘探储层预测技术领域。
背景技术
在储层预测地球物理技术方面,近年来物探技术的不断进步,以烃类检测为核心的有效储层预测技术已被广泛应用于石油勘探和开发的各个阶段。在储集体空间刻画方面,地震沉积学技术方法为沉积体系演化及岩性圈闭赋存规律等的研究提供了重要技术手段。地震地貌学和地震岩性学是地震沉积学的两个重要分支学科,其中地震地貌学通过三维地震资料平面及空间特征恢复古地貌为沉积演化、指导勘探发挥了重要作用,地震岩性学为如何将地震数据信息转化为岩性信息提供了理论技术依据。在岩性圈闭的刻画识别上,多种地震属性的应用、地质统计学反演等多种反演手段的应用以及基于虚拟现实系统三维可视化解释系统的出现为储集体三维空间刻画提供了多种技术手段。在有效储层预测方面,结合烃类检测技术进行有效储层预测是提高储层预测精度的重要手段之一。由于高频衰减、低频共振烃类检测理论的不断创新和一系列软件的推出和应用,烃类检测实现了由单纯的叠后地震资料预测向叠前叠后联合预测的发展,含油气范围圈定的精度不断提高,提升了勘探开发井的钻探成功率。
双相介质Biot理论认为,地震波穿过双相介质时,固相和流相之间产生相对位移并发生相互作用,产生第二纵波。第二纵波速度很低,且极性与第一纵波相反,由于流体的存在引起地震波的“低频增加、高频衰减”。但双相介质包括了若干必需的参数,且孔隙、岩性结构的不同也会形成低速层,若形成相同的波阻抗差就会形成假亮点,造成烃类检测技术存在多解性。钟国城在《塔里木盆地新和地区中生界地震烃类检测及有利圈闭预测》一文中,提出了通过离散傅里叶变换将时间域的三维地震数据体分为频率域间隔5Hz的离散能量数据体,以开展针对中生界目的层位的烃类检测。研究成果在研究区具有适用性验证及其预测结果对该地区寻找隐蔽圈闭油气藏具有重要的指导意义,但傅里叶变换存在时间和频率分辨率相对较低的局限性,造成烃类检测的误差增大。
公开号为CN107193041A的中国专利申请文件公开了一种基于双相介质的有效储层正演模拟的方法及系统,该方法包括如下步骤:确定储层的成藏物性下限,高于成藏物性下限的储层划分为有效储层;计算有效储层内岩石骨架的纵波速度、流体的纵波速度;根据岩石骨架的纵波速度和流体的纵波速度,分别计算岩石骨架的反射系数和流体的反射系数;将岩石骨架的反射系数和流体的反射系数分别与地震的实际子波褶积,然后将褶积的各正演结果进行叠加,得到有效储层的地震正演模型。该方法虽然提高了储层描述中正演精度和吻合度,但这只是储层预测中的过程之一,并没有给出模型在应用于勘探实践中如何获得参数、以及如何根据地震正演模型预测有效储层。
发明内容
本发明的目的是提供一种相控正演约束下的有效储层预测方法及系统,以解决目前有效储层预测存在的预测精度低的问题。
本发明为解决上述技术问题而提供一种相控正演约束下的有效储层预测方法,该预测方法包括以下步骤:
1)搜集目标区域三维叠后地震数据、测井数据和岩心数据,进行数据准备;
2)利用获取的地震数据确定地震可识别最小等时地层单元及与之匹配的最小等时沉积单元,并利用最小等时沉积单元的测井数据圈定最小地震层序格架内有利砂体分布区;
3)对地震可识别最小等时地层单元的古地貌进行定量恢复,得到地震可识别最小等时地层单元古地貌,以古地貌中坡折线和坡脚线之间的区域作为砂体卸载区;
4)将最小地震层序格架内有利砂体分布区和砂体卸载区的重合部分作为进行砂体预测的有利区;
5)利用目标区域已钻井进行实钻井正演,确定合并砂泥岩互层的最佳隔层厚度参数及地震波主频、有效频带范围和低频共振点频率;
6)对有利区的叠后偏移地震数据进行分频处理,得到有效频带范围内高频地震分频数据和低频共振点频率的地震分频数据;
7)对所述的地震可识别最小等时地层单元进行分频属性反演,确定所述的最小等时沉积单元内有利储层分布,得到有利储层分布区;
8)以低频共振点频率为轴线在设定的低频范围内计算低频共振点频率地震分频数据的极值与平均值之差,确定低频异常区;
9)将确定的有利储层分布区、有利区和低频异常区进行叠加,将三者重合范围作为有效储层分布范围。
本发明首先通过沉积相带及古地貌结合锁定有利区,使用分频地震数据确定地震可识别最小等时地层单元及与之匹配的最小等时沉积单元,确保沉积界面的等时性,使地震属性提取及反演的可靠性得到了保证;其次通过模型正演得到地震频谱变化参数,作为有利区目的层段的井旁地震频谱信息对比依据,圈定未知区有利储层;再根据储层含流体后地震低频谐振现象,突出油气信息形成的地震响应,叠加后得到最终的有效储层分布区。本发明以“沉积-古地貌-地震可识别最小等时地层单元-有利储层-有效储层”步步逼近,同时着眼于储层含流体后地震低频谐振现象,大大提高了有效储层预测精度。
进一步地,本发明还给出了如何确定可识别最小等时地层单元,所述的地震可识别最小等时地层单元的确定过程为:根据分频地震数据体确定随着频率变化而地震相位不变的相邻两个波峰或波谷,所述相邻两个波峰或波谷之间的沉积地层为地震可识别最小等时地层单元。
进一步地,本发明还给出了如何确定最佳隔层厚度参数,所述的最佳隔层厚度参数的确定过程为:根据目标区内实钻井信息,合并砂泥岩薄互层同类项为砂泥岩互层后进行实钻井正演,将正演模型与地震数据进行匹配,调整合并的隔层厚度参数,使正演模型与实际地震数据的频率和波形达到最佳匹配,此时的隔层厚度参数为最佳隔层厚度参数。
进一步地,本发明还给出了如何确定低频共振点频率,所述的低频共振点频率的确定过程为:根据正演模型与实际地震数据的频率和波形达到最佳匹配的砂泥岩互层正演的最佳匹配模型,进行流体弹性参数正演,以调谐能量最强的频率作为低频共振点频率。
进一步地,本发明还给出了如何进行古地貌定量恢复,所述步骤3)中古地貌定量恢复采用沉积速率、残余厚度法实现。
进一步地,为了保证所确定的有效储层精度更高,该方法还包括根据最佳隔层厚度参数对目标区域测井曲线进行降频处理,对有利区的叠后偏移地震数据进行高频成分补偿,使测井数据和地震数据相互匹配,以完成井震储层精细标定。
本发明还提供了一种相控正演约束下的有效储层预测系统,该预测系统包括处理器和存储器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
1)搜集目标区域三维叠后地震数据、测井数据和岩心数据,进行数据准备;
2)利用获取的地震数据确定地震可识别最小等时地层单元及与之匹配的最小等时沉积单元,并利用最小等时沉积单元的测井数据圈定最小地震层序格架内有利砂体分布区;
3)对地震可识别最小等时地层单元的古地貌进行定量恢复,得到地震可识别最小等时地层单元古地貌,以古地貌中坡折线和坡脚线之间的区域作为砂体卸载区;
4)将最小地震层序格架内有利砂体分布区和砂体卸载区的重合部分作为进行砂体预测的有利区;
5)利用目标区域已钻井进行实钻井正演,确定合并砂泥岩互层的最佳隔层厚度参数及地震波主频、有效频带范围和低频共振点频率;
6)对有利区的叠后偏移地震数据进行分频处理,得到有效频带范围内高频地震分频数据和低频共振点频率的地震分频数据;
7)对所述的地震可识别最小等时地层单元进行分频属性反演,确定所述的最小等时沉积单元内有利储层分布,得到有利储层分布区;
8)以低频共振点频率为轴线在设定的低频范围内计算低频共振点频率地震分频数据的极值与平均值之差,确定低频异常区;
9)将确定的有利储层分布区、有利区和低频异常区进行叠加,将三者重合范围作为有效储层分布范围。
进一步地,所述的地震可识别最小等时地层单元的确定过程为:根据分频地震数据体确定随着频率变化而地震相位不变的相邻两个波峰或波谷,所述相邻两个波峰或波谷之间的沉积地层为地震可识别最小等时地层单元。
进一步地,所述的最佳隔层厚度参数的确定过程为:根据目标区内实钻井信息,合并砂泥岩薄互层同类项为砂泥岩互层后进行实钻井正演,将正演模型与地震数据进行匹配,调整合并的隔层厚度参数,使正演模型与实际地震数据的频率和波形达到最佳匹配,此时的隔层厚度参数为最佳隔层厚度参数。
进一步地,所述的低频共振点频率的确定过程为:根据正演模型与实际地震数据的频率和波形达到最佳匹配的砂泥岩互层正演的最佳匹配模型,进行流体弹性参数正演,以调谐能量最强的频率作为低频共振点频率。
附图说明
图1是本发明相控正演约束下的有效储层预测方法的流程图;
图2是本发明实施例中A地区地震可识别最小等时地层单元古地貌图;
图3是本发明实施例中A地区确定的有利区示意图;
图4是本发明实施例中A地区实钻井正演模型与实际地震资料匹配图;
图5是本发明实施例中A地区低频共振点Y频率的地震分频剖面示意图;
图6是本发明实施例中A地区有利储层分布区示意图;
图7是本发明实施例中A地区有效储层分布区示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明。
本发明相控正演约束下的有效储层预测方法的实施例
本发明首先使用分频地震数据确定地震可识别最小等时地层单元及与之匹配的最小等时沉积单元,确保沉积界面的等时性,使地震属性提取及反演的可靠性得到了保证;然后通过沉积相带及古地貌结合剔除非有利砂体沉积区,锁定有利区;其次通过模型正演得到与地震道精细标定的测井降频参数及砂体含烃后地震频谱变化参数,作为有利区目的层段的井旁地震频谱信息对比依据,圈定未知区有利储层;再根据储层含流体后地震低频谐振现象,突出油气信息形成的地震响应,叠加后得到最终的有效储层分布区。下面以A地区为例对本发明有效储层预测方法的过程进行说明,具体实施流程如图1所示,具体过程如下。
1、数据准备
采集A地区三维叠后偏移地震数据、测井数据和岩心数据,并进行整理与分类。
2、利用A地区的地震数据、测井数据和岩心数据,圈定范围,建立底图。
本发明通过统计最小等时沉积单元测井曲线形态、幅度、韵律变化,分析地层岩性、砂体厚度、沉积物粒度变化、流体性质等,结合岩心颜色、沉积构造、古生物标志等,指示出有利的沉积亚相(及微相)发育区,即为最小地震层序格架内有利砂体分布区,过程如下:
2.1确定研究范围和目的层段,根据测井数据和岩心数据对目的层段测井曲线进行频谱分析,分解出不同级别的沉积旋回,建立A地区1#-12#井单井、连井的不同级别的沉积旋回对比剖面;
2.2根据采集的叠后偏移地震数据的分频地震数据体确定随着不同频率的变化、而地震相位始终不变的相邻两个稳定地震波形(波峰或者波谷),作为顶底界面,其间沉积地层即为地震可识别最小等时地层单元;
2.3根据步骤2.2中确定的地震可识别最小等时地层单元卡定测井与之匹配的短期沉积旋回,本发明称之为最小等时沉积单元;
2.4根据步骤1中得到的测井数据、岩心数据分析步骤2.3中最小等时沉积单元内的测井相、岩心相、岩相、沉积韵律等沉积特征,建立最小等时沉积单元内的单井、连井的沉积剖面;
2.5结合区域物源、沉积体系,根据步骤2.4得到的最小等时沉积单元内的单井、连井的沉积剖面,编制最小等时沉积单元平面沉积相图、圈定最小地震层序格架内有利砂体分布区,作为底图。
3、进行古地貌恢复,圈定砂体卸载区。
3.1根据沉积速率、残余厚度法等进行步骤2.2中所述的地震可识别最小等时地层单元的古地貌定量恢复,获得地震可识别最小等时地层单元古地貌图(如图2所示);
3.2在图2所示的古地貌图上圈定坡降大的地貌陡坎和坡折,勾画坡折线和坡脚线,两线之间是侵蚀和砂体沉积作用活跃区域,易发生河道下切等作用形成的厚层砂体及坡移扇砂体沉积,以两线之间圈出砂体卸载区。
4、确定有利区。
将步骤2得到的底图与步骤3.2得到的砂体卸载区进行叠合,得到可能厚层砂体分布区,本发明称之为有利区如图3(此步骤旨在剔除非有利砂体沉积区,减少因不同物源、不同相带、不同岩性配比等形成的其他因素的干扰)。
5、进行A地区1#-12#实钻井正演,记录合并砂泥岩互层的最佳隔层厚度参数及砂体含烃类后低频能量变化参数。
5.1根据1#-12#实钻井,合并砂泥岩薄互层同类项(合并砂岩或者泥岩集中发育段)变为砂泥岩互层后进行实钻井正演,把正演模型与实际地震资料进行匹配如图4所示,此过程中不断的调整合并的隔层厚度参数,直至得到的正演模型与实际地震资料的频率和波形达到最佳匹配,即正演模型的频率和波形与实际地震资料的频率和波形最接近,最吻合。此时记录达到最佳匹配的隔层厚度参数,记录为X,此时的正演模型为砂泥岩互层正演的最佳匹配模型。
5.2根据步骤5.1中的砂泥岩互层正演的最佳匹配模型,得到地震波的主频、有效频带范围;再利用流体替换进行流体弹性参数正演,观察地震记录频谱变化(含油气后频率域能量向低频移动),记录调谐能量最强的频率,以该频率作为低频增强的频率点Y。本实施例在进行流体弹性参数正演时,在18Hz低频段是调谐能量最强的优势频率,这里的“低频”范围是15-20Hz,优势频率是18Hz,即Y为18Hz。
为了验证低频增强的频率点Y的正确性,本发明对地震数据体进行Morlet小波频谱成像,根据成像结果来进行验证。
6、对步骤4中有利区的叠后偏移地震数据进行分频处理,得到不同频率的分频地震数据,得到步骤5.2中所述的有效频带范围内高频地震分频数据及低频共振点频率Y的地震分频数据如图5所示,高频地震分频数据中的高频的范围与原始地震资料有关,受控于原始地震资料主频与有效频宽,是在有效频宽范围内的相对高频值。
7、圈定最小地震层序格架内有效储层分布区。
7.1根据步骤5.1中所得到的合并同类项中最佳隔层厚度参数X,对目的层段测井曲线进行降频处理,对步骤4中有利区的叠后偏移地震数据适度补偿高频成分,使两种数据(测井数据和地震数据)相互匹配,完成井震储层精细标定;该过程是利用测井资料计算出的反射系数序列,与估算子波不断褶积,得到一个与井旁地震道相匹配的反射系数序列,其同时具有测井资料的高频成分信息和地震资料的低频成分信息;井震储层精细标定是明确有效储层的地震反射波形特征,是提高储层预测精度、提高反演可靠性的前提;
7.2对步骤2.2中所述的地震可识别最小等时地层单元进行分频属性反演,得到步骤2.3中所述的最小等时沉积单元内有利储层分布预测图,圈出有利储层分布区如图6所示;
7.3对经过步骤6中得到的低频共振点频率Y的地震分频数据体,通过时频分析技术得到单道二维谱,在以步骤5.2中所述的低频共振点频率Y为轴线的一定低频范围内,计算极值与平均值之差,得到平面ABV吸收属性图,圈出低频异常区;
7.4将步骤7.2中得到的预测有利储层分布区、步骤7.3中得到的低频异常区与步骤4中得到的有利区叠加,其三者重合范围即为有效储层分布范围(如图7所示)。
本发明相控正演约束下的有效储层预测系统的实施例
本发明的相控正演约束下的有效储层预测系统为计算机等具备数据处理能力的设备,该预测系统包括:存储器、处理器及网络模块。存储器、处理器以及网络模块相互之间直接或间接地电性连接,以实现数据的传输或交互。例如,这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。存储器包括至少一个可以软件或固件(firmware)的形式存储于存储器中的软件功能模块,处理器通过运行存储在存储器内的软件程序以及模块,从而执行各种功能应用以及数据处理,即实现本发明方法实施例中的相控正演约束下的有效储层预测方法。
其中,存储器可以是,但不限于,随机存取存储器(Random Access Memory,RAM),只读存储器(Read Only Memory,ROM),可编程只读存储器(Programmable Read-OnlyMemory,PROM),可擦除只读存储器(Erasable Programmable Read-Only Memory,EPROM),电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory,EEPROM)等。其中,存储器用于存储程序,处理器在接收到执行指令后,执行所述程序。
处理器可能是一种集成电路芯片,具有信号的处理能力。上述的处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(NetworkProcessor,NP)等。还可以是数字信号处理器(DSP))、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规处理器等。
网络模块用于通过网络建立电子设备与外部通信终端之间的通信连接,实现网络信号及数据的收发操作。上述网络信号可包括无线信号或者有线信号。
本发明首先使用分频地震数据确定地震可识别最小等时地层单元及与之匹配的最小等时沉积单元,确保沉积界面的等时性,使地震属性提取及反演的可靠性得到了保证;其次针对测井和地震子波非线性复杂的影射关系中采用分频属性反演砂体及物性,对于砂泥岩波阻抗差异小、复杂地质背景地区尤其适用;第三对地震数据体采用成熟的Morlet小波频谱成像技术,对分析局部信号、小于λ/4的地质体较傅立叶变换等更具有优势。本发明以“沉积-古地貌-地震可识别最小等时地层单元-有利储层-有效储层”步步逼近的技术思路预测甜点分布,同时着眼于储层含流体后地震低频谐振现象,大大提高了有效储层预测精度。
本发明在某油田探区已知油藏区及预探区经过实践检验,预测精度达到85%以上,较以往的同类技术指导钻探成功率提高了10%,高效指导了油气勘探,满足小型断陷湖盆砂泥岩薄互储层预测精度要求。
Claims (6)
1.一种相控正演约束下的有效储层预测方法,其特征在于,该预测方法包括以下步骤:
步骤1)搜集目标区域三维叠后地震数据、测井数据和岩心数据,进行数据准备;
步骤2)利用获取的地震数据确定地震可识别最小等时地层单元及与之匹配的最小等时沉积单元,并利用最小等时沉积单元的测井数据圈定最小地震层序格架内有利砂体分布区;
步骤3)对地震可识别最小等时地层单元的古地貌进行定量恢复,得到地震可识别最小等时地层单元古地貌,以古地貌中坡折线和坡脚线之间的区域作为砂体卸载区;
步骤4)将最小地震层序格架内有利砂体分布区和砂体卸载区的重合部分作为进行砂体预测的有利区;
步骤5)利用目标区域已钻井进行实钻井正演,确定合并砂泥岩互层的最佳隔层厚度参数及地震波主频、有效频带范围和低频共振点频率;
步骤6)对有利区的叠后偏移地震数据进行分频处理,得到有效频带范围内高频地震分频数据和低频共振点频率的地震分频数据;
步骤7)对所述的地震可识别最小等时地层单元进行分频属性反演,确定所述的最小等时沉积单元内有利储层分布,得到有利储层分布区;
步骤8)以低频共振点频率为轴线在设定的低频范围内计算低频共振点频率地震分频数据的极值与平均值之差,确定低频异常区;
步骤9)将确定的有利储层分布区、有利区和低频异常区进行叠加,将三者重合范围作为有效储层分布范围;
所述的最佳隔层厚度参数的确定过程为:根据目标区内实钻井信息,合并砂泥岩薄互层同类项为砂泥岩互层后进行实钻井正演,将正演模型与地震数据进行匹配,调整合并的隔层厚度参数,使正演模型与实际地震数据的频率和波形达到最佳匹配,此时的隔层厚度参数为最佳隔层厚度参数;
所述的低频共振点频率的确定过程为:根据正演模型与实际地震数据的频率和波形达到最佳匹配的砂泥岩互层正演的最佳匹配模型,进行流体弹性参数正演,以调谐能量最强的频率作为低频共振点频率。
2.根据权利要求1所述的相控正演约束下的有效储层预测方法,其特征在于,所述的地震可识别最小等时地层单元的确定过程为:根据分频地震数据体确定随着频率变化而地震相位不变的相邻两个波峰或波谷,所述相邻两个波峰或波谷之间的沉积地层为地震可识别最小等时地层单元。
3.根据权利要求1所述的相控正演约束下的有效储层预测方法,其特征在于,所述步骤3)中古地貌定量恢复采用沉积速率、残余厚度法实现。
4.根据权利要求1所述的相控正演约束下的有效储层预测方法,其特征在于,该方法还包括根据最佳隔层厚度参数对目标区域测井曲线进行降频处理,对有利区的叠后偏移地震数据进行高频成分补偿,使测井数据和地震数据相互匹配,以完成井震储层精细标定。
5.一种相控正演约束下的有效储层预测系统,其特征在于,该预测系统包括处理器和存储器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
1)搜集目标区域三维叠后地震数据、测井数据和岩心数据,进行数据准备;
2)利用获取的地震数据确定地震可识别最小等时地层单元及与之匹配的最小等时沉积单元,并利用最小等时沉积单元的测井数据圈定最小地震层序格架内有利砂体分布区;
3)对地震可识别最小等时地层单元的古地貌进行定量恢复,得到地震可识别最小等时地层单元古地貌,以古地貌中坡折线和坡脚线之间的区域作为砂体卸载区;
4)将最小地震层序格架内有利砂体分布区和砂体卸载区的重合部分作为进行砂体预测的有利区;
5)利用目标区域已钻井进行实钻井正演,确定合并砂泥岩互层的最佳隔层厚度参数及地震波主频、有效频带范围和低频共振点频率;
6)对有利区的叠后偏移地震数据进行分频处理,得到有效频带范围内高频地震分频数据和低频共振点频率的地震分频数据;
7)对所述的地震可识别最小等时地层单元进行分频属性反演,确定所述的最小等时沉积单元内有利储层分布,得到有利储层分布区;
8)以低频共振点频率为轴线在设定的低频范围内计算低频共振点频率地震分频数据的极值与平均值之差,确定低频异常区;
9)将确定的有利储层分布区、有利区和低频异常区进行叠加,将三者重合范围作为有效储层分布范围;
所述的最佳隔层厚度参数的确定过程为:根据目标区内实钻井信息,合并砂泥岩薄互层同类项为砂泥岩互层后进行实钻井正演,将正演模型与地震数据进行匹配,调整合并的隔层厚度参数,使正演模型与实际地震数据的频率和波形达到最佳匹配,此时的隔层厚度参数为最佳隔层厚度参数;
所述的低频共振点频率的确定过程为:根据正演模型与实际地震数据的频率和波形达到最佳匹配的砂泥岩互层正演的最佳匹配模型,进行流体弹性参数正演,以调谐能量最强的频率作为低频共振点频率。
6.根据权利要求5所述的相控正演约束下的有效储层预测系统,其特征在于,所述的地震可识别最小等时地层单元的确定过程为:根据分频地震数据体确定随着频率变化而地震相位不变的相邻两个波峰或波谷,所述相邻两个波峰或波谷之间的沉积地层为地震可识别最小等时地层单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811244512.7A CN111090117B (zh) | 2018-10-24 | 2018-10-24 | 一种相控正演约束下的有效储层预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811244512.7A CN111090117B (zh) | 2018-10-24 | 2018-10-24 | 一种相控正演约束下的有效储层预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111090117A CN111090117A (zh) | 2020-05-01 |
CN111090117B true CN111090117B (zh) | 2021-11-19 |
Family
ID=70391555
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811244512.7A Active CN111090117B (zh) | 2018-10-24 | 2018-10-24 | 一种相控正演约束下的有效储层预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111090117B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113945970B (zh) * | 2020-07-15 | 2024-06-07 | 中国石油化工股份有限公司 | 一种致密砂岩储层预测方法 |
CN114252935B (zh) * | 2020-09-25 | 2024-05-28 | 中国石油天然气股份有限公司 | 火山岩风化壳风化淋滤的波及深度预测方法及装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6829537B2 (en) * | 2003-03-10 | 2004-12-07 | Conocophillips Company | Contour-based attribute seismic mapping |
CN103439741B (zh) * | 2013-09-16 | 2015-10-14 | 中国石油大港油田勘探开发研究院 | 一种零值法单砂体逐层剥离预测方法 |
CN105467447B (zh) * | 2014-09-04 | 2018-11-20 | 中国石油化工股份有限公司 | 相控趋势能量匹配的地震保幅评价方法 |
CN106707337B (zh) * | 2016-11-15 | 2018-08-21 | 西安石油大学 | 基于分频地震反射能量地层切片进行沉积体系解释的方法 |
CN107121699B (zh) * | 2017-02-28 | 2019-06-21 | 中国石油化工股份有限公司 | 一种地震相控制下的沉积微相识别方法 |
CN107728230B (zh) * | 2017-09-11 | 2019-06-25 | 中国石油大港油田勘探开发研究院 | 一种基于地球物理三联技术的岩性油藏预测方法 |
CN108680951A (zh) * | 2018-03-22 | 2018-10-19 | 中国地质大学(北京) | 一种基于地震信息判断煤层气富集沉积控制作用的方法 |
-
2018
- 2018-10-24 CN CN201811244512.7A patent/CN111090117B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111090117A (zh) | 2020-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2018340369B2 (en) | Method and device for determining thin interlayer | |
CN102707317B (zh) | 一种利用地震波吸收衰减特征进行储层分析的方法 | |
CN109425896A (zh) | 白云岩油气储层分布预测方法及装置 | |
EP3488267B1 (en) | Seismic spectral balancing | |
CN104237945B (zh) | 一种地震资料自适应高分辨处理方法 | |
CN102109611B (zh) | 一种利用低频地震属性预测油藏优质储层的方法 | |
CN103527184A (zh) | 一种白云岩储层的预测方法和系统 | |
CN101349764B (zh) | 一种地震旋回分析方法 | |
Xue et al. | Recent developments in local wave decomposition methods for understanding seismic data: application to seismic interpretation | |
CN105116449B (zh) | 一种弱反射储层的识别方法 | |
CN109425900A (zh) | 一种地震储层预测方法 | |
CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
CN108508489B (zh) | 一种基于波形微变化匹配的地震反演方法 | |
CN111090117B (zh) | 一种相控正演约束下的有效储层预测方法及系统 | |
CN107339099B (zh) | 一种确定储层岩性的方法和装置 | |
CN105092343A (zh) | 一种基于叠前道集去除薄层调谐效应的方法 | |
CN113109875A (zh) | 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法 | |
CN106501851B (zh) | 一种地震属性优化方法及装置 | |
CN115629417B (zh) | 基于地震沉积学的多尺度融合及相控颗粒滩刻画方法 | |
Ningkai et al. | Stepped and detailed seismic prediction of shallow-thin reservoirs in Chunfeng oilfield of Junggar Basin, NW China | |
CN116068663A (zh) | 基于磁震联合低频建模的火成岩波阻抗反演方法 | |
CN104345347B (zh) | 一种用于致密含气砂岩储层预测的测井曲线恢复方法 | |
CN113009567B (zh) | 泥包砂型薄储层地震特征的识别方法 | |
Wang et al. | Application of multiple attributes fusion technology in the Su-14 Well Block | |
CN117111157A (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 |