【技术实现步骤摘要】
本专利技术属于油气地球物理勘探领域中的数据处理技术，涉及一种置信叠前扩散滤 波方法。
技术介绍
在地震资料处理中，噪声的压制通常采用滤波的方法，针对噪声的不同特征 及叠前、叠后处理，通常采用多种滤波方法，其中具有代表性方法包括：奇异值分解 (SVD) (Ursin 和 Zheng, 1985、Wenkai, 2006) [1] [2]、F_xy 域图像特征值滤波（Ulrych et al.，1999、Trickett, 2003) [3]、F_k 域滤波（March and Bailey, 1983) [4]、τ -p 域滤波 (Tatham et al.，1983) [5]、神经网络等一些模糊滤波去噪（Hashemi et al.，2008) [6]、扩 散滤波（孙夕平，2004;石岩峻，2010、杨威，2011; Kamel Baddari, Jalal Ferahtia, Tahar Aifa，etal, 2011) [7-10]等。 1984年，由Koenderink[ll]发现扩散滤波等价于一个热传导的物理过程。非线 性各向异性扩散方程是由Perona和Malik[12]在1990年把它作为图像处理中有选择性地 保留边界信息的工具而首先提出。P-Μ方程有可能是病态的，1996年You Y L[13]通过正 则化的方法消除了这种病态性。1998年，Weickert [14-15]对Perona和Malik提出的公式 进行了改进，引进了结构张量D，2002年提出最优旋转不变性的新型算法。 以上算法虽然可以用于对叠...
【技术保护点】
置信叠前扩散滤波方法，其特征在于，包括循环进行的如下S1‑S3步骤： S1.对读入数据利用①式进行高斯滤波，高斯滤波完成后求取叠前数据在X方向和Y方向上的一阶导数和；,‑‑‑‑‑‑①其中高斯滤波函数；u为地震振幅，为地震振幅的梯度，为尺度函数,t表示进程中的编号；S2.利用②式构建结构张量，‑‑‑‑②；当定义结构张量D的特征向量、，对应特征值时，扩散系数对实对称矩阵D，利用JACOBI方法并结合③式求解出D的特征值u=(u1，u2)及对应的特征向量v1、v2；‑‑‑‑‑③其中k为相干系数k=(u1‑u2)3,为最大向量方向的波动扩散率，线状置信度；根据求出的特征向量和系数来计算出扩散系数a、b、c，求解①式；S3.将S1中计算出的一阶导数和及S2中计算出的扩散系数a、b、c；带入④式， 解出等号左边的一阶导数；‑‑‑‑‑‑‑‑‑‑‑‑④反推原函数，得到uO=uO（x,y,t+1），以uO作为步骤S1中新的读入数据，重复步骤S1‑S3 ，直至达到预先设定的循环终止条件；输出最后一次循环得到的数据uO，作为滤波后的地震剖面叠前道集。
【技术特征摘要】
1. 置信叠前扩散滤波方法，其特征在于，包括循环进行的如下S1-S3步骤： S1.对读入数据利用①式进行高斯滤波，高斯滤波完成后求取叠前数据 在X方向和Y方向上的一阶导数和：；其中高斯滤波函数；U为地震振幅，Vi为地震振幅的梯度， ?为尺度函数，t表示进程中的编号；52. 利用②式构建结构张量，当定乂结构张量0的特征向量= h ，对应特征值u = (ul, u2)时， 扩散系数对实对称矩阵D，利用JACOBI方法并结合③式求解出D的特征值u=(ul，u2)及对应的 特征向量vl、v2其中k为相干系数k=(ul-u2)3，α为最大向量方向的波动扩散率，线状置信度根据求出的特征向量▼和系数i来计算出扩散系数a、b、c，求解①式；53. 将S1中计算出的一阶导数：和·及S2...
【专利技术属性】
技术研发人员：刘力辉，
申请(专利权)人：成都晶石石油科技有限公司，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。