一种基于改进霍夫变换的自适应时频峰值滤波消噪方法技术

本发明专利技术涉及一种基于改进霍夫变换的自适应时频峰值滤波消噪方法，属于地球物理技术领域。针对直线型同相轴地震记录提出一种多层分数阶傅里叶霍夫，对地震数据进行同相轴包络识别，建立与同相轴走势相近的滤波轨线方程，实现时频峰值滤波最优轨线的自适应选取，从而充分降低滤波偏差。本发明专利技术能够在强噪声条件下快速得出检测结果，具有较高的准确性及较强的抗干扰能力。

An adaptive time frequency peak filtering denoising method based on improved Hof transform

The invention relates to an adaptive time frequency peak filtering noise elimination method based on an improved Hof transform, belonging to the geophysical technical field. According to the linear phase axis seismic records, this paper proposed a multi-layer fractional Fourier Hof for event recognition on the envelope of seismic data, a similar trend with the trajectory equation filtering phase axis, select adaptive frequency peak filtering to achieve optimal trajectory, so as to reduce the filter error. The invention can obtain the detection result quickly under the strong noise condition, and has higher accuracy and stronger anti-interference ability.

【技术实现步骤摘要】
一种基于改进霍夫变换的自适应时频峰值滤波消噪方法
本专利技术属于地球物理
，尤其是指针对包含直线型同相轴的地震记录提出了一种基于多层分数阶傅里叶霍夫变换的自适应时频峰值滤波消噪方法，该方法通过一种新型的霍夫变换获取地震数据中的同相轴包络斜率，进而建立具有该斜率的滤波轨线方程并沿着这些轨线对含噪地震数据进行样本重采样，其目的是使得采样后数据中有效信号线性度充分降低，最后对重采样后数据进行时频峰值滤波有效压制噪声，为地震勘探信号消噪处理及地震数据中有效信号识别提供了新方法。
技术介绍
地震勘探作为一种重要的物理勘探手段，在油田的开发中起着重要的作用。野外地震数据采集和地震资料的数字处理技术是地震勘探的两个重要环节。只有得到高质量的地震剖面图，才能为石油勘探提供可靠的依据。在利用检波器进行实际的野外数据采集过程中，有效地震信号经常会受到各种噪声的干扰(规则噪声和随机噪声)，这些都直接影响地震解释工作的准确性。相比之下，随机噪声是不规则的、没有规律的，且相邻道间彼此是互不相关的，它们没有固定的频率，几乎分布于整个频带，严重影响着地震记录的信噪比。因此随机噪声压制是地震数据处理中的重点和难点。近几年来，对于随机噪声压制，除了目前发展较为成熟的几种消噪算法，如自适应滤波、经验模态分解、多项式拟合、f-x预测滤波、小波变换、维纳滤波等，时频峰值滤波算法被成功地应用于陆地地震勘探随机噪声压制，并先后发展出多种改进算法，如局部径向轨线时频峰值滤波、二次轨线时频峰值滤波及分段轨线时频峰值滤波等。时频峰值滤波方法较上述成熟算法在强随机噪声、非平稳信号及先验信息未知等条件下存在明显的优势，但现有时频峰值滤波方法中轨线参数多针对某一地震记录人为设置，灵活性较差。
技术实现思路
本专利技术提供一种基于改进霍夫变换的自适应时频峰值滤波消噪方法，以解决现有轨线型时频峰值滤波方法中存在的轨线参数及形状固定，灵活性和适应性较差问题。目的在于针对包含直线型同相轴的地震记录，提出一种改进的霍夫变换，并将其应用与地震记录中的同相轴形态识别，从而实现时频峰值滤波方法中最优轨线的自适应选取。本专利技术采取的技术方案是，包括下列步骤：步骤一：一张含噪的、包含直线型同相轴的地震记录，记录中的反射同相轴就相当于图像中的边缘信息，为了便于对地震数据进行形态识别，首先对地震数据进行灰度处理，将数据值转换到0-255之间，处理后的地震数据为一张图像数据，接着，利用改进的霍夫变换实现地震记录中同相轴形态识别，具体过程为：先利用公式对灰度处理后的图像数据进行多层分数阶傅里叶变换，将地震数据变换到时域和频域之间的一个中间域，将其旋转到从时间轴指向频率轴的一个角度，得到多组不同阶数的傅里叶变换结果；再将变换结果转换成由模值ρ和相角θ构成的极坐标形式，在极坐标内，对数据进行共轭对称，并沿模值方向进行一维傅里叶反变换便可得到一张由多条正弦线构成的正弦谱图；实际上，从多层分数阶傅里叶变换到转换成极坐标，再到进行一维傅里叶反变换，整个过程相当于对地震数据进行了一次霍夫变换，得到的正弦谱图中多条正弦线的交点，即正弦谱矩阵中的数据最大值就对应了地震数据中的同相轴；对于直线型同相轴的地震记录，通过搜索正弦谱图中的峰值便可得到与其对应的地震数据中的某条同相轴在极坐标内的模值ρpeak和相角θpeak，再根据极坐标与原图像坐标之间的关系式ρ＝xcosθ+ysinθ便可以得到地震数据内同相轴包络线的斜率kl＝ctanθpeak；于是沿着多条斜率为kl的直线型滤波轨线对含噪地震记录进行样本重采样，有效提高采样后信号的线性度，降低滤波偏差；步骤二：根据检测到的同相轴包络斜率kl建立与同相轴走势接近的滤波轨线方程yl＝klxl+bl，其中xl为地震数据点的横坐标，代表炮检距方向，yl为地震数据点的纵坐标，代表时间方向，kl为轨线斜率，bl为截距，决定轨线与地震数据纵坐标的交点位置，改变截距bl的数值便会生成多条具有相同斜率的直线轨线，沿着这些轨线对地震数据进行样本重采样，便可实现对地震数据沿同相轴形态方向的数据提取；步骤三：对采样后的数据样本，利用时频峰值滤波方法进行噪声压制，由于对地震数据进行样本重采样时所采用的轨线斜率是通过检测同相轴形态得出的，因此轨线与同相轴趋势接近，沿着这些轨线重采样后的地震数据有效信号线性度充分提高，主频充分降低，因此对其进行时频峰值滤波时所采用的滤波窗长比传统时频峰值滤波方法中长很多，较长的窗长在完整保持有效成分的同时可以有效地压制随机噪声。本专利技术所述步骤一中对于直线型同相轴，地震数据中的某一点(x,y)与极坐标内的参数对(ρ，θ)满足如下关系：ρ＝xcosθ+ysinθ,θ∈[0,π]，其中θ为直线法线与x轴的夹角，ρ为图像空间中原点到直线的距离。本专利技术以提高时频峰值滤波方法的适用性和灵活性，有效压制地震记录中随机噪声为目的，针对包含直线型同相轴地震记录，基于多层分数阶傅里叶提出了一种改进的霍夫变换，并将其应用于同相轴形态检测，实现了时频峰值滤波方法中最优滤波轨线的自适应选取。改进的霍夫变换与现有霍夫变换相比具有较高的分辨率，没有梯度信息取向的损失，忽略了从笛卡尔到极坐标映射的零填充，在同相轴形态检测方面计算速度快，准确性高。多层分数阶傅里叶变换在频域内具有较高的分辨率，没有梯度信息取向的损失，忽略了从笛卡尔到极坐标映射的零填充，计算速度明显提高。而现有的随机霍夫变换有时不能检测到短线段；线段检测器对噪声比较敏感，在嘈杂的图像中检测性能明显下降。在重采样过程中由于位移-时间域内的样本格点与新作用域内的样本格点不一致，部分样本点不能落到具体的地震道上。如果丢弃这些“脱道”的点，采样得到的序列中就会出现一些断点和尖峰脉冲。本专利技术通过在采样前在相邻地震道间插入一些虚拟地震道的方式对地震数据进行扩充，使得原本“脱道”的点能够落在具体的地震道上，从而有效降低采样误差。本专利技术的有益效果是提出了一种基于多层分数阶傅里叶变换实现的改进霍夫变换方法，并将其应用于直线型地震记录中同相轴形态检测，获取同相轴斜率。改进的霍夫变换与现有霍夫变换相比具有较高的分辨率，没有梯度信息取向的损失，忽略了从笛卡尔到极坐标映射的零填充，在同相轴形态检测方面计算速度快，准确性高。同相轴形态的准确识别轨线型时频峰值滤波方法中最优滤波轨线的自适应选取为提供了先决条件。此外，沿滤波轨线进行的插值重采样方案通过在相邻地震道间插入虚拟地震道的方式，有效避免了采样点的脱道现象，降低了重采样过程造成的误差。以上方面均对充分降低滤波偏差，完整保持有效信号具有重要意义。附图说明图1(a)图像空间内直线各参数情况，其中点P1、P2和P3为直线上的三个点，直线与坐标原点距离为ρ，直线与y轴夹角为θ；图1(b)霍夫变换参数空间内正弦曲线谱图，其中三条正弦曲线与点P1、P2和P3三点对应，三条线相交于一点，该点的模值和相角即对应了图1a中直线与坐标原点的距离ρ及其与y轴夹角θ；图2(a)具有不同噪声强度的含噪地震记录，该记录中包含三条直线型同相轴，其中区域R1信噪比较低，有效信号被淹没其中难以辨识；区域R2包含较强的脉冲噪声；区域E1存在交叉轴；图2(b)变换后得到的参数空间的正弦谱图，其中存在三个明显的交点(黑色方框标本文档来自技高网...

【技术保护点】
一种基于改进霍夫变换的自适应时频峰值滤波消噪方法，其特征在于，包括下列步骤：步骤一：一张含噪的、包含直线型同相轴的地震记录，记录中的反射同相轴相当于图像中的边缘信息，为了便于对地震数据进行形态识别，首先对地震数据进行灰度处理，将数据值转换到0‑255之间，处理后的地震数据为一张图像数据，接着，利用改进的霍夫变换实现地震记录中同相轴形态识别，具体过程为：先利用公式

【技术特征摘要】
1.一种基于改进霍夫变换的自适应时频峰值滤波消噪方法，其特征在于，包括下列步骤：步骤一：一张含噪的、包含直线型同相轴的地震记录，记录中的反射同相轴相当于图像中的边缘信息，为了便于对地震数据进行形态识别，首先对地震数据进行灰度处理，将数据值转换到0-255之间，处理后的地震数据为一张图像数据，接着，利用改进的霍夫变换实现地震记录中同相轴形态识别，具体过程为：先利用公式(0≤α≤1)对灰度处理后的图像数据进行多层分数阶傅里叶变换，将地震数据变换到时域和频域之间的一个中间域，将其旋转到从时间轴指向频率轴的一个角度，得到多组不同阶数的傅里叶变换结果；再将变换结果转换成由模值ρ和相角θ构成的极坐标形式，在极坐标内，对数据进行共轭对称，并沿模值方向进行一维傅里叶反变换便可得到一张由多条正弦线构成的正弦谱图；实际上，从多层分数阶傅里叶变换到转换成极坐标，再到进行一维傅里叶反变换，整个过程相当于对地震数据进行了一次霍夫变换，得到的正弦谱图中多条正弦线的交点，即正弦谱矩阵中的数据最大值就对应了地震数据中的同相轴；对于直线型同相轴的地震记录，通过搜索正弦谱图中的峰值便可得到与其对应的地震数据中的某条同相轴在极坐标内的模值ρpeak和相角θpeak，再根据极坐标与原图像坐标之间的关系式ρ＝xcosθ+ysinθ便可以得到地震数据内同相轴包络线的斜率kl＝cta...

【专利技术属性】
技术研发人员：田雅男，李月，吴宁，于朋君，李默，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林,22