[发明专利]一种基于全局优化的叠前地震信号同相轴拉平方法

摘要

本发明公开了一种基于全局优化的叠前地震信号同相轴拉平方法，利用时窗中心点作为种子点，不断滑动时窗求得道集间相似度最大时的移动量作为种子点的移动量，对种子点移动量进行全局优化，最后通过三次样条函数以种子点移动量作为原始数据插值得到拉平后的叠前道集振幅。本发明通过近邻相似传播选取基准地震道，可以自动选取相似度高的道集作为基准道，不需要人工干预，相较于直接从相似度矩阵中相似度较高的道集中任选一道准确度更高。本发明在得到种子点移动量后设计新的目标函数和约束条件对种子点移动量进行全局优化，解决拉平过程中的“突跳”以及波形失真的问题。

主权项

一种基于全局优化的叠前地震信号同相轴拉平方法，其特征在于，包括以下步骤：S1、初始化道集拉平参数；S2、选取基准道；所述步骤S2包括以下分步骤：S21、计算任意两个道集的相似度矩阵C；S22、以相似度矩阵为基础初始化吸引度矩阵与归属度矩阵；S23、迭代更新吸引度矩阵与归属度矩阵；S24、计算使吸引度矩阵与归属度矩阵之和最大的道集k；S25、判断迭代次数是否达到指定次数，若是则进入步骤S3，否则进入步骤S26；S26、判断道集k是否与上次迭代时结果一致，若是则进入步骤S3，否则返回步骤S23；S3、计算种子点移动量；所述步骤S3包括以下分步骤：S31、求取两个道集的最大相似度矩阵Cmax并定义最优移动量矩阵S；S32、计算矩阵Cmax对应分位数阈值处的值cm；S33、统计大于cm值的道集个数，选择个数最多的行，该行相似度对应的移动量即为当前时窗种子点的移动量；S34、插值得到其余时窗种子点的移动量；S35、对各时窗种子点移动量进行全局优化；所述步骤S35包括以下分步骤：S351、计算道集间相似度最大的移动量矩阵；S352、计算水平方向和垂直方向的位移差分矩阵；S353、判断水平方向和垂直方向的位移差分矩阵是否满足约束条件，若是则进入步骤S4，否则选择相似度次优的移动量组成新的移动量矩阵并返回步骤S352；S4、同相轴拉平。