[发明专利]置信叠前扩散滤波方法

摘要

置信叠前扩散滤波方法，包括循环进行的如下S1-S3步骤：S1.对读入数据利用①式进行高斯滤波，高斯滤波完成后求取叠前数据；在X方向和Y方向上的一阶导数。S2.构建结构张量，利用JACOBI方法求解出特征值及对应的特征向量；S3.利用S1中计算出的一阶导数及S2中计算出的扩散系数a、b、c反推得到u_O=u_O（x,y,t+1），以u_O作为步骤S1中新的读入数据，重复步骤S1-S3，直至达到预先设定的循环终止条件；输出最后一次循环得到的数据u_O，作为滤波后的地震剖面叠前道集。该置信叠前扩散滤波方法不仅去噪效果较好，而且有一定保幅性，在同相轴的保持方面优于前人的扩散滤波效果。

主权项

置信叠前扩散滤波方法，其特征在于，包括循环进行的如下S1‑S3步骤： S1.对读入数据利用①式进行高斯滤波，高斯滤波完成后求取叠前数据在X方向和Y方向上的一阶导数和；,‑‑‑‑‑‑①其中高斯滤波函数；u为地震振幅，为地震振幅的梯度，为尺度函数,t表示进程中的编号；S2.利用②式构建结构张量，‑‑‑‑②；当定义结构张量D的特征向量、，对应特征值时，扩散系数对实对称矩阵D，利用JACOBI方法并结合③式求解出D的特征值u=(u1，u2)及对应的特征向量v1、v2；‑‑‑‑‑③其中k为相干系数k=(u1‑u2)³,为最大向量方向的波动扩散率，线状置信度；根据求出的特征向量和系数来计算出扩散系数a、b、c，求解①式；S3.将S1中计算出的一阶导数和及S2中计算出的扩散系数a、b、c；带入④式， 解出等号左边的一阶导数；‑‑‑‑‑‑‑‑‑‑‑‑④反推原函数，得到u_O=u_O（x,y,t+1），以u_O作为步骤S1中新的读入数据，重复步骤S1‑S3 ，直至达到预先设定的循环终止条件；输出最后一次循环得到的数据u_O，作为滤波后的地震剖面叠前道集。