[发明专利]基于主频波数域均值滤波的地震数据去噪方法及系统

权利要求书

1.一种基于主频波数域均值滤波的地震数据去噪方法，其特征在于，包括：

对获得的地震数据进行分段处理；

基于傅里叶频谱计算，统计每段数据中有效频带的范围；

截取每段数据中主频波数域成分，并将所述主频波数域成分变换到曲波域，获得曲波系数矩阵；所述截取每段数据中主频波数域成分，具体为利用二维傅里叶变换将时空域数据变换到频率波数域，根据频率波数域的共轭对称性和确定的有效频带范围，截取正主频波数域成分；所述将所述主频波数域成分经变换到曲波域，具体为：对主频波数域成分采用2维 FFT 变换到频率域，在频率域采用插值方法获得不同尺度及不同方向的成分，利用抛物窗对各个尺度、方向的成分进行局部化，通过周期化技术将局部化后的成分映射到原点的仿射区域，最后对处理后数据做2维反傅里叶变换得到曲波系数；

利用能量均值滤波算法对所述曲波系数矩阵进行随机噪声的压制，避免对弱同相轴的破坏；

对于去噪后的曲波系数，依次进行曲波逆变换和二维傅里叶逆变换得到每段时空域去噪地震数据。

2.如权利要求1所述的一种基于主频波数域均值滤波的地震数据去噪方法，其特征在于，所述基于傅里叶频谱计算，统计每段数据中有效频带的范围，具体为求取各道数据的傅里叶频谱，构建整个地震数据的平均傅里叶频谱，获得符合地震数据的有效频带范围，对主频成分截取提供依据。

3.如权利要求1所述的一种基于主频波数域均值滤波的地震数据去噪方法，其特征在于，所述利用能量均值滤波算法对所述曲波系数矩阵进行随机噪声的压制，具体包括：在长度为L个采样点、宽度为W个采样点的曲波系数矩阵c_mn中移动一个长度为a个采样点的方窗，然后比较窗口中的均值E_w和矩阵的均值E_cm，如果E_wE_cm，将窗口内的曲波系数置零；如果E_wE_cm，保留窗口内的曲波系数。

4.如权利要求1所述的一种基于主频波数域均值滤波的地震数据去噪方法，其特征在于，所述对获得的地震数据进行分段处理，具体为，根据预设道数和时间采样点数，确定横向和纵向的段数，进而实现地震数据的分段。

5.一种基于主频波数域均值滤波的地震数据去噪系统，其特征在于，包括：

分段单元，其用于对获得的地震数据进行分段处理；

有效频带统计单元，其用于基于傅里叶频谱计算，统计每段数据中有效频带的范围；

曲波系数计算单元，其用于截取每段数据中主频波数域成分，并将所述主频波数域成分变换到曲波域，获得曲波系数矩阵；其中，所述截取每段数据中主频波数域成分，具体为利用二维傅里叶变换将时空域数据变换到频率波数域，根据频率波数域的共轭对称性和确定的有效频带范围，截取正主频波数域成分；所述将所述主频波数域成分经变换到曲波域，具体为：对主频波数域成分采用2维FFT变换到频率域，在频率域采用插值方法获得不同尺度及不同方向的成分，利用抛物窗对各个尺度、方向的成分进行局部化，通过周期化技术将局部化后的成分映射到原点的仿射区域，最后对处理后数据做2维反傅里叶变换得到曲波系数；

噪声压制单元，其用于利用能量均值滤波算法对所述曲波系数矩阵进行随机噪声的压制，避免对弱同相轴的破坏；

去噪数据获取单元，其用于对于去噪后的曲波系数，依次进行曲波逆变换和二维傅里叶逆变换得到每段时空域去噪地震数据。