基于频率域滤波的面波信号压制.docxVIP

基于频率域滤波的面波信号压制 面波是地震勘探中遇到的主要干扰波之一。面波具有能量强、频率低但频带变化范围相对较大的特点。在野外采集过程中, 通常采用检波器组合的方式来压制面波能量。煤田地震勘探中常用主频60 Hz的高频检波器, 可对面波能量起到一定程度的压制作用。由于面波的频带往往与有效波频带形成交集, 采用常规的一维频率域滤波 (如带通滤波) 无法保证在不伤害有效波的情况下压制面波能量。实际工作中压制面波的手段主要有内切带通滤波、f-k滤波以及自适应滤波。内切滤波原理主要是通过定义面波的发育范围, 仅在该范围内进行带通滤波来压制面波能量, 从而尽可能地影响到该范围之外的有效波能量。自适应滤波是Grisys处理系统中比较有特色的压制面波的模块, 它通过定义面波与有效波的主频, 以及面波的最大视速度范围, 对记录内的面波分布范围进行自适应滤波来压制面波能量。f-k滤波是利用有效波与干扰波的视速度差异, 通过二维傅里叶变换将地震数据从t-x域 (时间—空间域) 变换至f-k域 (频率—波数域) , 从而使有效波与干扰波在频率波数域分开, 进而可以保证在不伤害有效波的情况下将干扰波去除。面波通常具有线性特征, 视速度很低, 一般每秒在几百米, 而有效波视速度很高, 一般每秒都在几千米以上, 利用f-k滤波可以利用两者视速度的差异较好地压制面波能量, 提高叠前资料的信噪比。本文着重介绍f-k滤波压制面波的方法及经验。 1 频率波数谱设计 f-k滤波方法是基于二维傅里叶变换原理 F (f, k) 通常称为频率波数谱 (图1) 。 在频率波数谱上, 每一条通过原点直线上的f/k的值都是一个常数, 它表示在t-x域上的某一个视速度值。利用有效信号与规则干扰视速度差异这一性质, 可以在频率波数谱上设计滤波器: D是频率波数谱上的某一区域 (图1中有效信号区、规则干扰区) 。通过设计这样一个频率波数谱应用在地震数据的频率波数谱F (f, k) , 可以得到输出频率波数谱Fapos; (f, k) , 输出频率波数谱值包含了输入频率波数谱F (f, k) 上区域D上所包含的相关信息, 即: 在实际应用中, 根据有效信号的频率波数谱和干扰噪声的频率波数谱在f-k平面上的分布情况, 来确定区域D的形状, 确定其边界, 从而得到各种不同类型的f-k域滤波器。 2 面波压制范围 图2是豫东某地区野外原始单炮记录, 从单炮记录上可以看出, 主要干扰波面波比较发育, 其视速度在100~400 m/s范围之间, 主频在12 Hz左右, 而有效波视速度较面波速度高得多, 均在3 000 m/s以上, 主频在30~60 Hz之间。 原始单炮记录的f-k谱如图3所示, 从f-k谱上分析可知, 面波干扰的线性特征保持较好, 视速度约在400 m/s。实际处理中使用CGG软件的f-k滤波模块FKFIL对面波噪声进行压制。选择压制速度带:50—100—1 000—1 500, 速度压制范围较分析得到的面波速度范围稍大些, 可以保证压制效果。作业流程中首先在共炮点道集上应用了FKFIL模块, 随后抽道集到共接收点道集又应用了一次FKFIL模块, 通过在共炮点道集与共接收点道集上分别应用f-k滤波, 可以提高对线形干扰的压制效果 在资料处理去噪环节中, 经常需要根据噪声的发育特征变换到不同域中处理, 有时共炮点道集上噪声特征更明显, 有时在共接收点道集或者CMP道集上特征更明显, 实际工作中可以根据噪声在不同域的特征发育情况, 在不同域的道集上分别应用f-k滤波, 这里分别在共炮点道集以及共接收点道集分别应用FKFIL模块, 这样做可以提高f-k滤波的效果, 但也可能会加重f-k滤波的混波效应, 故在实际中应根据资料的具体情况灵活把握尺度。需要注意的是在应用f-k滤波时, 输入的道集 (共炮点道集或共接收点道集) 不能有过多的缺道, 否则将影响f-k滤波压制线性干扰波的效果。在FKFIL模块中可通过道内插来减弱此类影响。 3 压制其他域f-k滤波 综上所述, 可以看出f-k滤波在压制具有线性特征的面波噪声方面可以获得比较好的效果, 对于其他不同类型的线性相干噪声, f-k滤波方法也同样适用。在实际应用f-k滤波压制线性相干噪声时, 要注意选择适当的视速度压制带, 叠前共炮点道集上应用f-k滤波后如果效果不是太理想的情况下, 可以变换到其他域, 如共接收点道集或者共中心点道集多次应用f-k滤波, 通常可以取得更好的效果。另外, f-k滤波不仅在压制线性相干噪声中可以取得不错的效果, 在压制其他类型噪声 (如多次波、随机噪声) 中也可以起到不错的效果。实际上, 对于其他各种不同类型的噪声, 只要其与有效波在频率波数谱上可以有效区分, 均可以利用f-k滤波的方法进行压制。由此可见