一种基于变分原理的地震数据信号去噪方法技术

本发明专利技术公开了一种基于变分原理的地震数据信号去噪方法，包括以下步骤：S1、在地震信号处理模型中输入测试信号f和参考信号g并选取采样点；S2、将测试信号和参考信号进行归一化处理，构造成为差异矩阵；S3、将差异矩阵转换为慢度矩阵；S4、基于费马原理和惠更斯‑菲涅尔原理，采用最小路径射线追踪方法在慢度矩阵中提取全局最小旅行时及对应的最短路径；S5、将获得的最短路径中每个采样点的旅行时时移量施加在测试信号上，获得去噪后的地震数据信号。本方法从理论和实现两个方面进行了修改，不仅提高了计算效率，增强了最短路径求取的稳定性，而且为非线性信号测量打开了新的视角。

【技术实现步骤摘要】
一种基于变分原理的地震数据信号去噪方法
本专利技术涉及石油勘探领域，尤其涉及一种基于变分原理的地震数据信号去噪方法。
技术介绍
地震信号测量的方法不是很多且各有优缺点，目前检索到的估计地震信号的非线性差异测量技术有以下几种：1、互相关方法该方法是利用随机信号的平稳性，计算两道信号之间的互相关，通过提取互相关的最大值确定信号的时移量。该方法只适应与平稳信号，计算量大。2、多项式拟合方法该方法假设地震信号在局部邻域道内的时差符合一个待定的多项式，通过拟合该多项式提取不同的时间差。该方法需要较大的计算量来估计未知参数，并且假定了不同道之间的时差关系是低阶平滑的，无法处理强变化的时移情况。3、非线性beamforming方法该方法利用局部的抛物二阶时差近似去搜寻局部的相关事件，同样地，该方法需要估计未知的非线性系数，计算量巨大，并且基于二阶近似的方法无法适用非规则的信号时移差异。4、常规的DTW方法。传统的DTW方法需要通过累积和回溯两步寻找到最短路径，可以提取非线性的时移量，但是计算量大，并且易受噪音影响，陷入局部极值中，使得结果不稳定。这些方法各有优缺点，但其仍然无法达到人们对地震信号测量操作的效率及稳定性等方面的需求。
技术实现思路
本专利技术目的是针对上述问题，提供一种基于变分原理的地震数据信号去噪方法。为了实现上述目的，本专利技术的技术方案是：一种基于变分原理的地震数据信号去噪方法，包括以下步骤：S1、在地震信号处理模型中输入测试信号f和参考信号g并选取采样点；S2、将测试信号和参考信号进行归一化处理，构造成为差异矩阵；S3、将差异矩阵转换为慢度矩阵；S4、基于费马原理和惠更斯-菲涅尔原理，采用最小路径射线追踪方法在慢度矩阵中提取全局最小旅行时及对应的最短路径；S5、将获得的最短路径中每个采样点的旅行时时移量施加在测试信号上，获得去噪后的地震数据信号。进一步的，所述步骤S1中的测试信号f和参考信号g中有效信号的个数要求相同，其两侧端点也相同并且要求测试信号f和参考信号g均具有次序性。进一步的，所述步骤S2中差异矩阵的构造公式为：d(x,y)＝[f(x)-g(y)]2,x＝1,...,MT；y＝1,...,NT其中，MT是测试信号的采样点数，NT是参考信号的采样点数。进一步的，所述步骤S3中差异矩阵转换为慢度矩阵的转换公式为：其中，vmin和vmax分别为最小和最大速度值，dmax为差异矩阵的最大值。进一步的，所述步骤S4中在慢度矩阵中提取全局最小旅行时及对应的最短路径的方法包括以下步骤：S41、将慢度矩阵分为矩形网格，在矩形网格单元顶点处建立节点，并确定起始节点和终止节点；S42、将起始节点设为初始震源节点s0，围绕s0搜索所有可能的临近计算节点S43、计算从当前震源节点s0到计算节点的旅行时，从中选择最小的旅行时记为对应的计算节点记为g1，则：S44、把当前计算节点g1设为新的震源节点s1，重新搜索当前震源节点的计算节点并计算对应的旅行时，最小旅行时记为对应的计算节点记为g2；S45、从初始震源节点s0到当前选择的最新计算节点的总旅行时为对应的路径为S46、重复上述过程，直到最终获得的计算节点为终止节点，获得全局最小旅行时及对应的路径：与现有技术相比，本专利技术具有的优点和积极效果是：本专利技术中采用基于变分原理的DTW方法测量两段地震信号的差异，有着其他方法不具备的优势，其具体优势和特点表现在以下几个方面：第一、方法效果的可靠性。本方法通过形成的慢度场进行射线追踪，提取最小旅行时和最短路径，可以有效地避免局部极值的产生，并可有效地压制噪音的干扰，求取的最短路径更稳健。另外，同常规的DTW方法相比，该方法计算效率更高。第二、理论上进行了提升。将DTW问题转换为两个核心问题，一个是差异矩阵的构建，一个是寻找最优路径问题。差异矩阵的构建可以选择信息的属性进行提取，并选择其他抽象距离来定义；最优路径求取在变分原理的框架下求取，提供了多种方法来求取信号的扰动时差，为非线性时差测量提供新了求解思路。第三、具有非常广泛的应用价值。本方法处理两段信号的差异，在地震勘探中有诸多这样需求，比如叠前地震数据叠加去噪、插值；成像剖面处理；成像道集后处理；参数反演中模拟数据与实际数据的差异测量等等方面。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为计算节点搜索方式示意图；图2为本专利技术的方法流程图；图3为待处理地震信号的示意图；其中左侧A信号为参考信号，右侧B信号为待处理信号；图4a为构建的差异矩阵示意图；图4b为构建的慢度矩阵示意图；图5a为前两次获得的最短路径示意图；图5b为最后获得的全局最短路径示意图；图6为处理后的地震信号示意图；其中左侧A信号为参考信号，右侧B信号为处理后的信号；图7a为原始的共炮道集；图7b为处理后的炮道集。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，所作的任何修改、等同替换、改进等，均应包含在本专利技术的保护范围之内。在地震信号处理中，往往需要对两道信号进行波形匹配，或测量两段信号的差异。一种有效的方法是采用动态时间规划方法测量两段信号的差异，该方法是一个典型的优化问题，其问题在于用满足一定条件的差异矩阵描述测试信号和参考信号时间差异关系，求解这两个信号匹配时，累计距离最小时所对应的最短路径。常规的DTW方法需要通过累积和回溯两步寻找到最短路径，可以提取非线性的时移量，但是计算量大，并且容易陷入局部极值中，不稳定。本专利技术提出基于变分原理求解DTW，将DTW技术主要分为差异矩阵形成和求取最短路径两个方面。差异矩阵的形成采用信号之间的欧式距离来实现；最短路径的求取则基于变分原理，采用最小路径射线追踪方法提取全局最小旅行时和最短路径。从理论和实现两个方面进行了修改，不仅提高了计算效率，增强了最短路径求取的稳定性，而且为非线性信号测量打开了新的视角。本专利技术目的是提供一种基于变分原理的DTW技术测量两段信号的差异的方法，并将该方法应用到叠前地震数据去噪，优化叠加中。本专利技术具体包括以下步骤：1、差异矩阵形成：首先，差异矩阵的构造前提要求信号满足以下条件：a：一致性条件：测试信号f(t)和参考信号g(t)两段信号中包含相同的有效信息，即有效信号的个数是一致的；b：单调性条件：有效的地震信号具有次序性；c：边界条件：地震信号两侧端点一致。在满足上述要求的假设上，差异矩阵可有下式构造：d(x,y)＝[f(x)-g(y)]2,x＝1,...,MT；y＝1,...,NT其中，MT和NT分别是两段信号的采样点数。2、慢度矩阵的构建：将差异矩阵转换为慢度矩阵：其中，vmin和vmax为给定的最小和最大速度值，dmax为差异矩阵的最大值。3、构造变分问题：将求取两段信号的非线性时移量，转换为求取慢度场中最短旅行时本文档来自技高网...

【技术保护点】
1.一种基于变分原理的地震数据信号去噪方法，其特征在于：包括以下步骤：S1、在地震信号处理模型中输入测试信号f和参考信号g并选取采样点；S2、将测试信号和参考信号进行归一化处理，构造成为差异矩阵；S3、将差异矩阵转换为慢度矩阵；S4、基于费马原理和惠更斯‑菲涅尔原理，采用最小路径射线追踪方法在慢度矩阵中提取全局最小旅行时及对应的最短路径；S5、将获得的最短路径中每个采样点的旅行时时移量施加在测试信号上，获得去噪后的地震数据信号。

【技术特征摘要】
1.一种基于变分原理的地震数据信号去噪方法，其特征在于：包括以下步骤：S1、在地震信号处理模型中输入测试信号f和参考信号g并选取采样点；S2、将测试信号和参考信号进行归一化处理，构造成为差异矩阵；S3、将差异矩阵转换为慢度矩阵；S4、基于费马原理和惠更斯-菲涅尔原理，采用最小路径射线追踪方法在慢度矩阵中提取全局最小旅行时及对应的最短路径；S5、将获得的最短路径中每个采样点的旅行时时移量施加在测试信号上，获得去噪后的地震数据信号。2.如权利要求1所述的基于变分原理的地震数据信号去噪方法，其特征在于：所述步骤S1中的测试信号f和参考信号g中有效信号的个数要求相同，其两侧端点也相同并且要求测试信号f和参考信号g均具有次序性。3.如权利要求2所述的基于变分原理的地震数据信号去噪方法，其特征在于：所述步骤S2中差异矩阵的构造公式为：d(x,y)＝[f(x)-g(y)]2,x＝1,...,MT；y＝1,...,NT其中，MT是测试信号的采样点数，NT是参考信号的采样点数。4.如权利要求3所述的...

【专利技术属性】
技术研发人员：葛大明，梁鸿贤，王延光，王常波，高丽，孔青华，揭景荣，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：北京,11