一种深度域储层能量密度的直接计算方法技术

本发明专利技术属于油气地震勘探领域，是提供一种直接利用深度域地震数据计算储层能量密度实现高精度储层识别的解释技术。本发明专利技术方法首先对叠后地震道数据进行深

【技术实现步骤摘要】
一种深度域储层能量密度的直接计算方法


[0001]本专利技术属于油气地震勘探解释领域，是一种直接通过深度域地震数据计算储层能量密度实现高精度储层预测的解释技术。

技术介绍

[0002]叠前深度偏移地震成像技术在地下复杂地质构造的高精度成像方面具有突出优势，相较于时间域地震资料，其不仅能够提供更为准确的构造信息，而且还能直接与测井资料进行对比，能够有效降低勘探风险，业界对深度域地震数据的反演和解释方法技术提出了日益迫切的需求。因此，探索深度域的储层预测方法具有重要意义。目前，常规的深度域储层预测方法，通常需要将深度域地震数据通过速度场转换至时间域，在时间域完成处理和解释，再将结果转回至深度域。这样的流程主要存在以下问题：
①
将数据从深度域转换至时间域的过程中，会损失其中的有效特征信息；
②
数据在不同域之间的转换过程中，也会在有效数据中引入误差；
③
域转换过程增加了计算成本和时间成本，不利于提高勘探开发工作的效率。因此，常规的深度域储层预测方法不能充分利用深度域地震数据中反映储层结构、岩性和流体情况等的有效信息，不利于挖掘深度域地震数据的应用价值和潜力，也不利于实现储层的精细描述。因此，直接利用深度域地震数据进行储层预测的方法亟待发展。

技术实现思路

[0003]本专利技术是要提供一种深度域储层能量密度的直接计算方法，通过在深度域直接计算地震数据的深度
‑
波数谱，利用低波数段的振幅能量数据实现储层能量密度的计算，并进行储层预测。/>[0004]本专利技术的目的在于提供一种深度域储层能量密度的直接计算方法，所述方法主要包括以下步骤：
[0005](1)对叠后地震道数据进行深
‑
波(深度
‑
波数)变换，求取对应的深度域波数谱：
[0006]X(d,k)＝DT[x(d)][0007]式中，DT[
·
]表示深
‑
波变换，x(d)表示深度域的地震道数据，X(d,k)表示进行深波变换后的深度
‑
波数谱，d表示深度，k表示波数；
[0008]⑵
选取最优参考波数k0，并利用该波数下的数据利用下式直接计算深度域的储层能量密度E：
[0009][0010]式中，k
opt
为最优参考波数，A
opt
为X(d,k)数据中最优参考波数处的瞬深振幅谱能量，其中最优参考波数按如下步骤确定：
[0011]a)选取叠后井旁道数据，并在地震低波数段范围内(约1～30km
‑1)，分别利用所有的波数作为参考波数进行能量密度的计算；
[0012]b)将所有波数下得到的能量映射成(深度d
‑
波数k)的二维分布图E(d,k)：
[0013]E(d,k)＝E[(d,k1)；(d,k2)；(d,k3)；...；(d,k
m
)][0014]式中，[k1,k2,k3,
…
,k
m
]为低波数范围内的参考波数，m为波数的个数；
[0015]c)利用测井解释结果进行约束，选取含流体储层位置处的能量曲线进行优势波数的确定，当能量值达到最大值时，即可获得最优的参考波数，即：
[0016][0017]式中，d
r
表示含流体储层位置处的深度点，表示求取k中的最大值；
[0018]⑶
对叠后地震剖面的多有地震道数据利用步骤(1)至
⑵
进行计算，获取叠后地震剖面对应的深度域储层能量密度，进而利用其结果识别储层；
[0019]⑷
重复上述步骤(1)至
⑶
，直到处理完整个工区的所有叠后地震道集，即可获得整个工区的三维储层预测结果，该预测结果可为储层的高精度预测提供支撑。
[0020]本专利技术的一种深度域储层能量密度的直接计算方法，具有如下优越性：
[0021](1)计算时直接利用了深度域的地震数据，避免了时深转换过程中速度误差引起的不确定性，求取的储层能量密度结果更加准确，可直接进行深度域的储层预测；
[0022](2)采用了最优波数高效交会确定方法，且利用了测井解释结果进行约束，基于此波数计算的储层能量密度可最大限度地突出储层的能量异常；
[0023](3)在无井或少井信息的情况下，该方法也可以获得储层的相对能量密度属性，为不依赖井而直接进行储层预测提供了一种可行的方法。
附图说明
[0024]图1是本专利技术技术实施例的深度域原始叠加地震剖面图。
[0025]图2是本专利技术技术确定的最佳参考波数示意图。
[0026]图3是本专利技术实施例深度域地震剖面计算的储层能量密度结果。
[0027]图4是利用本专利技术技术利用过井地震道计算的储层能量密度结果与测井曲线的对比图。
具体实施方式
[0028]本专利技术的具体实施方式如下：
[0029](1)输入某工区的三维叠后深度域地震数据；
[0030]⑵
利用控制井的测井信息及其对应的井旁地震道数据，给定参波数k范围内的一组向量k＝[k1,k2,k3,
…
,k
m
]，k
m
为波数的上限，按如下步骤首先确定最佳参考波数k
opt
：
[0031]a)在k＝k
i
时，按如下公式计算波数为k
i
时的储层能量密度E
i
：
[0032][0033]式中，A
i
为参考波数k
i
处的瞬深振幅谱能量；
[0034]b)将所有波数下得到的能量映射成(深度d
‑
波数k)的二维分布图E(d,k)：
[0035]E(d,k)＝E[(d,k1)；(d,k2)；(d,k3)；...；(d,k
m
)][0036]式中，[k1,k2,k3,
…
,k
m
]为低波数范围内的参考波数，m为波数的个数；
[0037]c)利用测井解释结果进行约束，选取含流体储层位置处的能量曲线进行优势波数的确定，当能量值达到最大值时，即可获得最优的参考波数，即：
[0038][0039]式中，d
r
表示含流体储层位置处的深度点，表示求取k中的最大值；
[0040]⑶
利用步骤
⑵
中确定的最优参考波数k
opt
直接计算深度域的储层能量密度E：
[0041][0042]式中，A
opt
为最优参考波数处的瞬深振幅谱能量。
[0043]⑷
重复步骤(1)至
⑶
即可获得所有地震道数据的深度域储层能量密度；
[0044]⑸
利用步骤
⑷
得到的储层能量密度结果进行储层的识别；<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种深度域储层能量密度的直接计算方法，其特征在于采用以下具体步骤：(1)对叠后地震道数据进行深
‑
波(深度
‑
波数)变换，求取对应的深度域波数谱：X(d,k)＝DT[x(d)]式中，DT[
·
]表示深
‑
波变换，x(d)表示深度域的地震道数据，X(d,k)表示进行深波变换后的深度
‑
波数谱，d表示深度，k表示波数；
⑵
选取最优参考波数k0，并利用该波数下的数据利用下式直接计算深度域的储层能量密度E：式中，k
opt
为最优参考波数，A
opt
为X(d,k)数据中最优参考波数处的瞬深振幅谱能量，其中最优参考波数按如下步骤确定：a)选取叠后井旁道数据，并在地震低波数段范围内(约1～30km
‑1)，分别利用所有的波数作为参考波数进行能量密度的计算；b)将所有波数下得到的能量映射成(深度d
‑
波数k)的二维分布图E(d,k)：E(d,k)＝E[(d,k1)；(d,k2)；(d,k3)；...；(d,k
m
)]式中，[k1,k2,k3,

【专利技术属性】
技术研发人员：罗鑫，陈学华，吕丙南，蒋伟，王胜艺，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：