利用抗假频凸集投影算法的规则缺失地震数据重建

发布时间：2021-11-17 13:21

　　在地震数据现场采集过程中,因受工区复杂地形等条件限制,采集到的地震数据大多不完整、不规则,且会呈现不同程度的空间假频现象。在均匀空间网格下,规则地震道缺失引起的假频与真实频谱相似,而常规基于傅里叶变换和凸集投影（POCS）的重建方法不能实现反假频重建。为此,针对线性同相轴规则缺失地震数据,通过在f-k域中采用倾角扫描方法拾取有效波能量,建立相应的选择函数,将其引入到POCS算法对数据进行反演计算,重建得到无假频数据。数值实验结果表明,文中所提方法在重建规则缺失地震数据的同时有效地消除了假频信息。 

【文章来源】：石油地球物理勘探. 2020,55(04)北大核心EICSCD

【文章页数】：10 页

【部分图文】：

25%规则采样重建结果及其f-k谱

设d(t, x)为原始地震数据,D(f, k)为其相应的f-k域振幅谱。为了简化,将频率和波数做归一化处理,令时间轴和距离轴上采样间隔都为1, 且f-k频谱关于频率轴对称,只需考虑正频率的频谱,得到归一化频率范围f∈(0,0.5)、归一化波数范围k∈(-0.5,0.5)(图1)。设定每条角度射线起始位置为频域原点(f,k)∈(0,0),从图1看出,可沿着角度射线方向求取射线上振幅能量之和。则每条射线所扫描的能量值为

为体现本文阈值参数的优势,将其与文献[28]中的线性阈值与指数阈值做重建对比,其重建后的信噪比公式采用 R S/Ν =20? lg ∥x 0 ∥ 2 ∥x-x 0 ∥ 2 ,其中x0表示原始模型数据,x表示重建结果。显然信噪比越高, 表明重建效果越好。 计算结果如图2所示。 可见在取得相同信噪比情况下,本文阈值参数所用迭代次数少,计算效率高;在相同迭代次数下,本文阈值参数重建信噪比最高。因此,本文阈值明显优于线性阈值和指数阈值。将选择函数W代入到传统的POCS迭代重建过程中,即把式(3)代入式(4)中,最终的抗假频POCS重建算法为

【参考文献】：
期刊论文
[1]利用压缩感知技术的离散正交S变换地震数据重建[J]. 赵子越,李振春,张敏.  石油地球物理勘探. 2020(01)
[2]应用二维非均匀曲波变换压制地震随机噪声[J]. 张华,刁塑,温建亮,黄光南,朱雯婷,白敏.  石油地球物理勘探. 2019(01)
[3]双曲Radon-ASVD方法压制叠前地震数据随机噪声[J]. 徐彦凯,曹思远,何元.  石油地球物理勘探. 2017(03)
[4]基于高维小波变换的高抗噪性边缘检测技术[J]. 王清振,张金淼,姜秀娣,翁斌,朱振宇,丁继才.  石油地球物理勘探. 2016(05)
[5]地面微地震资料τ-p变换噪声压制[J]. 宋维琪,刘太伟.  石油地球物理勘探. 2015(01)
[6]半径-斜率域加权反假频地震数据重建[J]. 黄小刚,王一博,刘伊克,常旭.  地球物理学报. 2014(07)
[7]基于jitter采样和曲波变换的三维地震数据重建[J]. 张华,陈小宏.  地球物理学报. 2013(05)
[8]基于seislet变换的反假频迭代数据插值方法[J]. 刘财,李鹏,刘洋,王典,冯晅,刘殿秘.  地球物理学报. 2013(05)
[9]地震数据反假频规则化方法研究[J]. 石颖,张振,王建民,郭雪豹.  地球物理学进展. 2013(01)
[10]二维小波变换与中值滤波联合去噪方法研究[J]. 张胤彬,张华,刘松.  中国煤炭地质. 2011(05)



本文编号：3500999