一种叠前地震道集AVO特征保持的评价方法技术

本发明专利技术涉及地球物理探测技术领域，特别涉及一种叠前地震数据AVO特征保持的评价方法，包括如下步骤，S1，获取叠前地震道集，计算叠前地震道集的截距P(t)和梯度G(t)，其中t为获取叠前地震道集的时间；S2，根据截距P(t)和梯度G(t)计算人工合成地震道集；S3，计算叠前地震道集和所述人工合成地震道集的一致性指示因子；S4，所述一致性指示因子与预设阈值的差值越小则AVO特征的保持性越好，所述一致性指示因子与预设阈值的差值越大则AVO特征的保持性越差。量化的评估结果能够直接反映出道集AVO特征的保持情况，提高对资料AVO特征保持评价的准确性及工作效率，通过对比处理前后指示因子的变化趋势，判断该步处理对道集AVO特征的影响情况，帮助优化处理流程与参数。帮助优化处理流程与参数。帮助优化处理流程与参数。

【技术实现步骤摘要】
一种叠前地震道集AVO特征保持的评价方法


[0001]本专利技术涉及地球物理探测
，特别涉及一种叠前地震道集AVO特征保持的评价方法。

技术介绍

[0002]叠前地震数据的AVO属性包含更多的地下介质信息，可用来反演地下介质的物性属性，是地震数据的重要属性之一。在岩性气藏中，通过反演技术从AVO特征中提取地下介质的物性属性是寻找、优选油气藏的重要依据。
[0003]目前的AVO特征保持评价方法多是定性评价，如残差法：将处理前和处理后的数据进行相减，观察相减得到的残差中是否包含有效信号；沿震幅分析法：观察标志层位地震波能量随偏移距的变化趋势是否在处理前后保持一致；井旁道正演对比法：利用测井数据建立井旁合成道集，对比合成道集的AVO特征与地震数据的AVO特征来判断AVO特征的保持情况。上述方法在实际使用中都存在一定的局限性：残差法和沿震幅分析法不仅只能定性判断，而且不适用于反褶积等提频处理前后道集AVO特征保持的评价；井旁道正演对比法需要工区内有测井数据，且该方法只适用于偏移处理后的地震道集的评价。

技术实现思路

[0004]本专利技术的目的在于克服现有技术中对叠前地震道集AVO特征在地震资料处理过程中是否得到保持只能定性进行评价的不足，提出了一种基于地震数据本身的量化评价道集AVO特征保持情况的评价方法。
[0005]一种叠前地震道集AVO特征保持的评价方法，包括如下步骤：
[0006]S1，获取叠前地震道集，计算所述叠前地震道集的截距P(t)和梯度G(t)，其中t为所述获取叠前地震道集的时间；
[0007]S2，根据所述截距P(t)和梯度G(t)计算人工合成地震道集；
[0008]S3，计算所述叠前地震道集和所述人工合成地震道集的一致性指示因子；
[0009]S4，根据所述一致性指示因子与预设阈值的差值判断AVO特征的保持性。
[0010]进一步的，步骤S3所述一致性指示因子由以下公式计算得到：
[0011][0012]其中，a表示所述叠前地震道集，b表示所述人工合成地震道集，CC是相关系数，DIS是欧几里得距离。
[0013]进一步的，所述欧几里得距离的计算公式为：
[0014][0015]其中，m为采样计算的时窗起始时刻，n为采样计算的时窗结束时刻，a
k
为第k个采样时刻的叠前地震道集，b
k
为第k个采样时刻的人工合成地震道集。
[0016]进一步的，所述相关系数的计算公式为：
[0017][0018]其中，m为采样计算的时窗起始时刻，n为采样计算的时窗结束时刻，a
k
为第k个采样时刻的叠前地震道集，为采样计算的时窗内叠前地震道集的平均值，b
k
为第k个采样时刻的人工合成地震道集，为采样计算的时窗内人工合成地震道集的平均值。
[0019]进一步的，步骤S1中所述叠前地震道集的截距P(t)和梯度G(t)的计算公式为二阶Aki—Richards反射系数公式。
[0020]进一步的，步骤S1还包括，对计算得到叠前地震道集的截距P(t)和梯度G(t)进行平滑处理。
[0021]进一步的，步骤S2中根据所述截距P(t)和梯度G(t)计算人工合成地震道集的计算公式为二阶Aki—Richards反射系数公式。
[0022]优选地，所述获取叠前地震道集的步骤包括，通过检波器获得纵波地震数据，采用地震波纵波处理方法对纵波地震数据进行处理，得到动校正后切除动校正拉伸畸变的角度域CDP道集。
[0023]优选地，所述角度域CDP道集包含的角度个数应大于3个，并且相邻的角度之间的差值相等。
[0024]基于相同的专利技术构思，提出一种叠前地震道集AVO特征保持的评价装置，其特征在于，包括至少一个处理器，以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行上述任一项所述的方法。
[0025]与现有技术相比，本专利技术的有益效果：
[0026]1、本专利技术求取叠前地震道集的截距和梯度属性，并根据截距和梯度属性计算出人工合成地震道集数据，通过组合实际地震道集和人工合成道集之间的相似性的相关系数和距离计算出一致性指示因子，一致性指示因子为一量化值，其取值分为在
‑
2至1之间，越接近1表示当两者越相似且震幅值大小接近，越小于1表示两者至少具有相似性差和震幅值差异大中的一项。因此指示因子越接近1，说明地震数据AVO特征保持性好，越远离1，说明AVO特征保持性差，量化的评估结果能够直接、清晰的反映出道集AVO特征的保持情况，极大的提高地震资料处理过程中对资料AVO特征保持的评价的准确性以及工作效率；
[0027]2、在地震资料处理过程中也能发挥积极的作用，可以在资料处理结束后对最终成果道集进行评估，也可在资料处理过程中对某一步骤或某几步处理步骤前后的道集数据进行动态的评估，指示因子变大，说明该步处理不仅保持之前的AVO特征，且道集质量得到提升。相反，指示因子变小，说明该步处理破坏了原有的AVO特征，通过对比处理前后指示因子的变化趋势，判断该步处理对道集AVO特征的影响情况，帮助优化处理流程与处理参数；
[0028]3、指示因子的变化率用于优化处理流程，减少低效环节，例如对偏移后的道集通
常还需要做一些处理，比如剩余时差校正、能量均衡、去噪等等，这些处理对道集质量的影响程度通过指示因子的变化率反应出来，保留变化率大的步骤，剔除变化率小的步骤，这样既保证了数据的质量，也进一步的减少处理步骤，提高了处理效率。
附图说明：
[0029]图1为一种基于地震道集AVO特征保持的评价方法的流程图；
[0030]图2为实施例1步骤S1得到的角度域CDP道集图；
[0031]图3为实施例1步骤S1通过二阶Aki—Richards反射系数公式求取的截距剖面图；
[0032]图4为实施例1步骤S1通过二阶Aki—Richards反射系数公式求取的梯度剖面图；
[0033]图5为实施例1步骤S2得到的人工合成地震道集图；
[0034]图6为实施例1步骤S3求得一致性指示因子的IF平面图。
具体实施方式
[0035]下面结合试验例及具体实施方式对本专利技术作进一步的详细描述。但不应将此理解为本专利技术上述主题的范围仅限于以下的实施例，凡基于本
技术实现思路
所实现的技术均属于本专利技术的范围。
[0036]实施例1
[0037]本实施例提供了一种基于地震数据本身的量化评价道集AVO特征保持情况的评价方法，包括如下步骤，
[0038]S1，获取叠前地震道集，计算所述叠前地震道集的截距P(t)和梯度G(t)，其中t为所述获取叠前地震道集的时间；
[0039]其中，所述获取叠前地震道集的步骤包括，通过检波器获得纵波地震数据，采用地震波纵波处理方法对纵本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种叠前地震道集AVO特征保持的评价方法，其特征在于，包括如下步骤：S1，获取叠前地震道集，计算所述叠前地震道集的截距P(t)和梯度G(t)，其中t为所述获取叠前地震道集的时间；S2，根据所述截距P(t)和梯度G(t)计算人工合成地震道集；S3，计算所述叠前地震道集和所述人工合成地震道集的一致性指示因子；S4，根据所述一致性指示因子与预设阈值的差值判断AVO特征的保持性。2.如权利要求1所述的一种叠前地震道集AVO特征保持的评价方法，其特征在于，步骤S3所述一致性指示因子由以下公式计算得到：其中，a表示所述叠前地震道集，b表示所述人工合成地震道集，CC是相关系数，DIS是欧几里得距离。3.如权利要求2所述的一种叠前地震道集AVO特征保持的评价方法，其特征在于，所述欧几里得距离的计算公式为：其中，m为采样计算的时窗起始时刻，n为采样计算的时窗结束时刻，a
k
为第k个采样时刻的叠前地震道集，b
k
为第k个采样时刻的人工合成地震道集。4.如权利要求2所述的一种叠前地震道集AVO特征保持的评价方法，其特征在于，所述相关系数的计算公式为：其中，m为采样计算的时窗起始时刻，n为采样计算的时窗结束时刻，a
k
为第k个采样时刻的叠前地震道集，为采样计算的时窗内叠前地震道集的平均值，b
k
为第k个采样时刻的人工合成地震道集，为采样计算的时窗内人工合成地震道集的...

【专利技术属性】
技术研发人员：张剑飞，马昭军，凌航，王延坤，王金龙，曹烈，沈杰，
申请(专利权)人：中国石油化工股份有限公司西南油气分公司，
类型：发明
国别省市：