一种气云区空间展布范围刻画方法和系统技术方案

本发明专利技术属于海上油气田勘探技术领域，涉及一种气云区空间展布范围刻画方法和系统，包括以下步骤：S1对地震数据进行倾角导向中值滤波；S2通过经过滤波的地震数据计算出2阶方差体；S3通过经过滤波的地震数据计算出瞬时频宽属性体；S4将2阶方差体和瞬时频宽属性体融合得到方差

【技术实现步骤摘要】
一种气云区空间展布范围刻画方法和系统


[0001]本专利技术涉及一种气云区空间展布范围刻画方法和系统，属于海上油气田勘探


技术介绍

[0002]海上油田在世界原油产量中占有很大比重，但广泛分布的“气云区”已经成为制约海上油田勘探开发的重要难题。“气云区”是指地层中由于天然气的聚集和活动，造成纵波地震剖面上出现成像模糊带。我国渤海湾油田存在多个气云区，气云区内地震资料品质差，整体表现为同向轴畸变、破碎与地震有效带宽的明显衰减，这给构造解释与储层预测带来极大困难。
[0003]准确刻画“气云区”的空间展布范围是资料品质评估、井网部署、工程平台选址以及后续针对“气云区”成像品质改善研究的基础，具有重要研究价值。但目前本领域还没有准确刻画“气云区”的空间展布范围的方法。

技术实现思路

[0004]针对上述问题，本专利技术的目的是提供了一种气云区空间展布范围刻画方法和系统，其能够快速、准确得到“气云区”空间展布范围的刻画结果。
[0005]为实现上述目的，本专利技术采取以下技术方案：一种气云区空间展布范围刻画方法，包括以下步骤：S1对地震数据进行倾角导向中值滤波；S2通过滤波的地震数据计算出2阶方差体；S3通过滤波的地震数据计算出瞬时频宽属性体；S4将2阶方差体和瞬时频宽属性体融合得到方差-瞬时频宽融合体，获得气云区空间展布范围。
[0006]进一步，步骤S2中的2阶方差体计算过程为：S2.1以经过滤波的地震数据为输入得到原始数据方差体；S2.2以原始数据方差体为输入得到2阶方差体。
[0007]进一步，步骤S3中计算出瞬时频宽属性体后，需要对瞬时频宽属性体进行三维高斯平滑滤波处理。
[0008]进一步，步骤S4中方差-瞬时频宽融合体采用以下公式获得：
[0009][0010]其中，Var-Bw为融合属性体，Var为2阶方差属性体，Bw为采用三维高斯平滑滤波处理后的瞬时频宽属性体，C1、C2、C3权重常数。
[0011]进一步，步骤S4中通过对方差-瞬时频宽融合体进行三维高斯平滑滤波处理获得气云区空间展布范围。
[0012]进一步，三维高斯平滑滤波处理通过在主测线、联络测线和时间方向上各取若干个样点，在各样点设置三维高斯平滑滤波器实现。
[0013]本专利技术还公开了一种气云区空间展布范围刻画系统，包括：滤波模块，用于对地震数据进行倾角导向中值滤波；方差体获取模块，用于通过经过滤波的地震数据计算出2阶方差体；瞬时频宽属性体模块，用于通过经过滤波的地震数据计算出瞬时频宽属性体；融合模块，用于将2阶方差体和瞬时频宽属性体融合得到方差-瞬时频宽融合体，通过方差-瞬时频宽融合体获得气云区空间展布范围。
[0014]进一步，融合模块中方差-瞬时频宽融合体采用以下公式获得：
[0015][0016]其中，Var-Bw为融合属性体，Var为2阶方差属性体，Bw为采用三维高斯平滑滤波处理后的瞬时频宽属性体，C1、C2、C3权重常数。
[0017]进一步，融合模块中中通过对方差-瞬时频宽融合体进行三维高斯平滑滤波处理获得气云区空间展布范围。
[0018]进一步，三维高斯平滑滤波处理通过在主测线、联络测线和时间方向上各取若干个样点，在各样点设置三维高斯平滑滤波器实现。
[0019]本专利技术由于采取以上技术方案，其具有以下优点：本专利技术能够快速、准确得到“气云区”空间展布范围的刻画结果。通过2阶方差体能够过滤掉大部分断层引起的方差体异常，而保留“气云区”存在造成的方差体响应使结果更加准确。将方差与瞬时频宽融合得到的融合属性，利用瞬时频宽属性体对2阶方差体进行了改造，从而进一步提升对“气云区”分布范围的刻画精度。
附图说明
[0020]图1是本专利技术一实施例中气云区空间展布范围刻画方法的流程图；
[0021]图2是本专利技术一实施例中某海上油田地震资料倾角导向中值滤波处理前后剖面效果对比，图2(a)是滤波处理前的剖面效果图，图2(b)是滤波处理前的剖面效果图；
[0022]图3是本专利技术一实施例中2阶方差体的示意图，图3(a)是原始数据方差体的示意图，图3(b)是由原始数据方差体获得的2阶方差体的示意图；
[0023]图4是本专利技术一实施例中瞬时频宽属性体的示意图，图4(a)是瞬时频宽属性体的示意图，图4(b)是经过三维高斯平滑后的瞬时频宽属性体的示意图；
[0024]图5是本专利技术一实施例中方差-瞬时频宽融合体的示意图；
[0025]图6是本专利技术一实施例中经过三维高斯平滑后的方差-瞬时频宽融合体的示意图，图6(a)和图6(b)分别对应0.5秒、1秒时间切片的方差-瞬时频宽融合体的示意图。
具体实施方式
[0026]为了使本领域技术人员更好的理解本专利技术的技术方向，通过具体实施例对本专利技术进行详细的描绘。然而应当理解，具体实施方式的提供仅为了更好地理解本专利技术，它们不应该理解成对本专利技术的限制。在本专利技术的描述中，需要理解的是，所用到的术语仅仅是用于描述的目的，而不能理解为指示或暗示相对重要性。
[0027]实施例一
[0028]以渤海湾盆地某一受大面积“气云区”影响的油田地震资料为例，本实施例公开了一种气云区空间展布范围刻画方法，包括以下步骤：
[0029]S1对该油田用于雕刻“气云区”空间范围的地震数据进行倾角导向中值滤波，在保证结构不变的前提下，尽可能提升资料信噪比。某海上油田地震资料倾角导向中值滤波处理前后剖面效果对比，如图2所示，图2(a)是滤波处理前的剖面效果图，图2(b)是滤波处理前的剖面效果图。
[0030]S2通过经过滤波的地震数据计算出2阶方差体，如图3所示。
[0031]步骤S2中的2阶方差体计算过程为：
[0032]S2.1以经过滤波的地震数据为输入得到原始数据方差体，原始数据方差体的结构示意图如图3(a)所示；
[0033]S2.2以原始数据方差体为输入得到2阶方差体，2阶方差体的结构示意图如图3(a)所示。
[0034]步骤S2.1中计算的原始数据方差体刻画断层与“气云区”造成的方差体异常，而步骤S2.2中计算的方差体属性能够过滤掉大部分断层引起的方差体异常，而保留“气云区”存在造成的方差体响应。
[0035]S3通过经过滤波的地震数据计算出瞬时频宽属性体，如图4所示。“气云区”内的地震资料由于有效频带明显衰减从而在瞬时频宽属性体上具有较明显响应。在计算出瞬时频宽属性体后，需要对瞬时频宽属性体进行三维高斯平滑滤波处理。图4(a)是瞬时频宽属性体的示意图，图4(b)是经过三维高斯平滑后的瞬时频宽属性体的示意图。其中，三维高斯平滑滤波处理通过在主测线、联络测线和时间方向上各取10个样点，在各样点设置三维高斯平滑滤波器实现。
[0036]S4将2阶方差体和瞬时频宽属性体融合得到方差-瞬时频宽融合体，其具体结构图如图5所示。
[0037]步骤S4中方差-瞬时频宽融合体采用以下公式获得：
[0038][0本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种气云区空间展布范围刻画方法，其特征在于，包括以下步骤：S1对地震数据进行倾角导向中值滤波；S2通过滤波的地震数据计算出2阶方差体；S3通过滤波的地震数据计算出瞬时频宽属性体；S4将所述2阶方差体和瞬时频宽属性体融合得到方差-瞬时频宽融合体，获得气云区空间展布范围。2.如权利要求1所述的气云区空间展布范围刻画方法，其特征在于，所述步骤S2中的2阶方差体计算过程为：S2.1以所述经过滤波的地震数据为输入得到原始数据方差体；S2.2以所述原始数据方差体为输入得到2阶方差体。3.如权利要求1所述的气云区空间展布范围刻画方法，其特征在于，所述步骤S3中计算出所述瞬时频宽属性体后，需要对所述瞬时频宽属性体进行三维高斯平滑滤波处理。4.如权利要求1所述的气云区空间展布范围刻画方法，其特征在于，所述步骤S4中方差-瞬时频宽融合体采用以下公式获得：其中，Var-Bw为融合属性体，Var为2阶方差属性体，Bw为采用三维高斯平滑滤波处理后的瞬时频宽属性体，C1、C2、C3权重常数。5.如权利要求1所述的气云区空间展布范围刻画方法，其特征在于，所述步骤S4中通过对所述方差-瞬时频宽融合体进行三维高斯平滑滤波处理获得气云区空间展布范围。6.如权利要求3或5所述的气云区空间展布范围刻画方法，...

【专利技术属性】
技术研发人员：杜昕，范廷恩，张显文，张晶玉，董建华，高云峰，范洪军，聂妍，马良涛，
申请(专利权)人：中海石油中国有限公司北京研究中心，
类型：发明
国别省市：