一种基于动态时间规整的地震数据加权叠加方法技术

本发明专利技术提供一种基于动态时间规整的地震数据加权叠加方法，包手以下步骤：S1、选择一个参考地震道，该参考地震道设置为叠前地震数据矩阵在空间方向上的算术均值；S2、通过下式计算每个叠前地震道与参考地震道之间的DTW匹配点，在CMP道集中寻找与参考地震道的各个采样点相关联的采样点；S3、根据各相关采样点的噪声设置各样点的权值；S4、将参考地震道的每个采样点更新为与之相关的采样点的加权平均值；S5、重复步骤S2～S4直到收敛，将参考地震道输出为最终的叠加地震道。本申请结合了DTW方法与加权叠加方法的优点，克服了地震道时间偏移误差和地震数据噪声对数据叠加的影响，通过在合成数据和真实地震数据上的测试，获得了良好的应用效果。

【技术实现步骤摘要】
一种基于动态时间规整的地震数据加权叠加方法
本专利技术涉及地震资料
，具体涉及一种基于动态时间规整的地震数据加权叠加方法。
技术介绍
共中心点道集(CMP)的叠加可以提高地震资料的信噪比，在地震资料解释中起着重要的作用。最简单的叠加方法是沿偏移距方向计算叠前数据的算术平均值，若叠前地震道集时间校准准确，并且仅包含高斯白噪声，简单的平均可获得最优的叠加结果。然而，由于静校准和速度拾取的不准确性，在实际的CMP道集中经常出现异常的偏移、不规则和非高斯随机噪声，平均叠加方法不能有效抑制存在的噪声，从而不能满足实际需要。为了克服上述问题，加权平均的叠加方法应运而生。叠前地震道在叠加过程中的重要性不同，其大小由其噪声水平决定。例如，Neelamani等人提出了一种迭代算法(Leave-Me-Out)来估计数据的信噪比，从而实现加权的地震道叠加；Rashed提出了智能叠加算法，该算法通过排除叠加中的有害样本，并对叠前数据偏移距较小的部分施加更大的权重，从而优化叠加信号的幅值；Liu等人利用基于局部相关算法的时变平滑权值进行叠加；Deng等人将AB相干和局部相似技术结合起来，处理存在极化反转异常时的CMP道集加权叠加问题；Xie等人利用主成分分析(PCA)获得了一种理想的NMO校正的数据矩阵，以减轻加权叠加中异常地震道带来的负面影响；Rashed引入了去除离群值的叠加技术(OIO)，从每个时间采样点中设置一个独特的离群值范围，这意味着相应的权重为零。然而，当叠前数据中存在较多的时间偏移异常地震道时，目前的加权叠加方法通常采用一些权重值很小甚至为零的权重来减小这些异常地震道对叠加结果的影响，但本专利技术的专利技术人经过研究发现，这些方法不能充分利用这些未对齐地震道的反射信息，即丢弃了过多的地层反射信息，因此亟需提出一种新的方法来解决这一问题。
技术实现思路
针对现有加权叠加方法通常采用一些权重值很小甚至为零的权重来减小这些异常地震道对叠加结果的影响，但这些方法不能充分利用这些未对齐地震道反射信息的技术问题，本专利技术提供一种基于动态时间规整的地震数据加权叠加方法。为了解决上述技术问题，本专利技术采用了如下的技术方案：一种基于动态时间规整的地震数据加权叠加方法，包括以下步骤：S1、选择一个参考地震道，该参考地震道设置为叠前地震数据矩阵在空间方向上的算术均值，即：其中，Ai(t)为第i个叠前地震道Ai在时间t时的采样值，C为叠前地震道数目，R(t)为参考地震道R在时间t时的采样值；S2、通过下式计算每个叠前地震道与参考地震道之间的动态时间规整即DTW匹配点，在CMP道集中寻找与参考地震道的各个采样点相关联的采样点：其中，γk是第k个点对pk＝(rn,aim)的差异，该差异表示连接参考地震道R＝[r1,r2,…,rn,…rN]与叠前地震道Ai＝[ai1,ai2,…,aim,…,aiM]两个道中采样点的连接线，具体用(rn-aim)2来计算，索引n并不一定等于m，且每个点对代表了在以下约束下的一条对齐路径：单调性：对齐路径pk＝(rn,aim)和pk-1＝(rn',aim')的任意两个相邻元素，满足不等式n-n'≥0和m-m'≥0；连续性：对齐路径pk＝(rn,aim)和pk-1＝(rn',aim')的任意两个相邻元素，满足不等式n-n'≤1和m-m'≤1；边界：第一个对齐路径是p1＝(r1,ai1)，最后一个对齐路径是pk＝(rN,aiM)；S3、根据各相关采样点的噪声设置各样点的权值：wp＝1/(dp)sdp＝|R(t)-Xp|其中，wp为Xp的权重，dp为Xp与参考地震道R(t)之间的距离，t为时间，p＝1，2，...，P，P为关联点数目，Xp为叠前地震道中与参考地震道中采样点相关联的点，s为用于控制权重函数的用户定义常量；S4、将参考地震道的每个采样点更新为与之相关的采样点的加权平均值：其中，R(t)即为由加权平均值更新后的参考地震道；S5、重复步骤S2～S4直到收敛，将参考地震道输出为最终的叠加地震道。与现有技术相比，本专利技术提供的基于动态时间规整的地震数据加权叠加方法，使用动态时间规整即DTW在全局时间约束下对齐地震反射事件，对地震反射同相轴进行时间偏差校正，迭代计算这些校正后的采样点的加权平均值，整个过程迭代执行以获得更好的叠加效果，因此本方法给合了DTW方法与加权叠加方法的优点，克服了地震道时间偏移误差和地震数据噪声对数据叠加的影响，能充分利用地震道的反射信息，通过在合成数据和真实地震数据上的测试，获得了良好的应用效果，有效证实了本申请方法的实用性。进一步，所述步骤S3中，用于控制权重函数的用户定义常量s取值为0～2。附图说明图1是本专利技术提供的基于动态时间规整的地震数据加权叠加方法流程示意图。图2a为用DTW对齐的两条两条地震道，图2b为使用DTW叠加的示意图。图3a为叠前地震道集，图3b为叠加真值，图3c为使用传统等权叠加法叠加效果示意图，图3d为使用局部相关加权叠加法叠加效果示意图，图3e为使用基于PCA的加权叠加法叠加效果示意图，图3f为使用本申请方法叠加效果示意图。图4a为CMP道集，图4b为使用传统等权叠加法叠加地震道结果示意图，图4c为使用局部相关加权叠加法叠加地震道结果示意图，图4d为使用基于PCA的加权叠加法叠加地震道结果示意图，图4e为使用本申请方法叠加地震道结果示意图。具体实施方式为了使本专利技术实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体图示，进一步阐述本专利技术。请参考图1、图2a和图2b所示，本专利技术提供一种基于动态时间规整的地震数据加权叠加方法，包括以下步骤：S1、选择一个参考地震道，该参考地震道设置为叠前地震数据矩阵在空间方向上的算术均值，即：其中，Ai(t)为第i个叠前地震道Ai在时间t时的采样值，C为叠前地震道数目，R(t)为参考地震道R在时间t时的采样值；如在图2b中，叠前地震道数目C＝4，采用本步骤中的公式计算得到的参考地震道如图2b中的Stackedtrace。S2、通过下式计算每个叠前地震道与参考地震道之间的动态时间规整即DTW匹配点，在CMP道集中寻找与参考地震道的各个采样点相关联的采样点：其中，γk是第k个点对pk＝(rn,aim)的差异，该差异表示连接参考地震道R＝[r1,r2,…,rn,…rN]与叠前地震道Ai＝[ai1,ai2,…,aim,…,aiM]两个道中采样点的连接线，即图2a中连接参考地震道R与叠前地震道Ai两个道中采样点的粗实线与细实线，具体用(rn-aim)2来计算，索引n并不一定等于m，且每个点对代表了在以下约束下的一条对齐路径：单调性：对齐路径pk＝(rn,aim)和pk-1＝(rn',aim')的任意两个相邻元素，满足不等式n-n'≥0和m-m'≥0，该约束本文档来自技高网...

【技术保护点】
1.一种基于动态时间规整的地震数据加权叠加方法，其特征在于，包括以下步骤：/nS1、选择一个参考地震道，该参考地震道设置为叠前地震数据矩阵在空间方向上的算术均值，即：/n

【技术特征摘要】
1.一种基于动态时间规整的地震数据加权叠加方法，其特征在于，包括以下步骤：
S1、选择一个参考地震道，该参考地震道设置为叠前地震数据矩阵在空间方向上的算术均值，即：



其中，Ai(t)为第i个叠前地震道Ai在时间t时的采样值，C为叠前地震道数目，R(t)为参考地震道R在时间t时的采样值；
S2、通过下式计算每个叠前地震道与参考地震道之间的动态时间规整即DTW匹配点，在CMP道集中寻找与参考地震道的各个采样点相关联的采样点：



其中，γk是第k个点对pk＝(rn,aim)的差异，该差异表示连接参考地震道R＝[r1,r2,…,rn,…rN]与叠前地震道Ai＝[ai1,ai2,…,aim,…,aiM]两个道中采样点的连接线，具体用(rn-aim)2来计算，索引n并不一定等于m，且每个点对代表了在以下约束下的一条对齐路径：
单调性：对齐路径pk＝(rn,aim)和pk-1＝(rn',aim')的任意两个相邻元素，满足不等式n-n'≥0和m-m'≥0；
连续...

【专利技术属性】
技术研发人员：宋承云，李凌轩，陈韵清，武子腾，
申请(专利权)人：重庆理工大学，
类型：发明
国别省市：重庆;50