统一计算设备架构下的F-X域预测滤波并行算法

杨先凤，贵红军*，傅春常

（1.西南石油大学计算机科学学院，成都 610500；2.西南民族大学计算机科学与技术学院，成都 610041）

（*通信作者电子邮箱gui18328655853@163.com）

0 引言

提高地震资料信噪比是进行高分辨率地震资料处理的关键，而恰当地去除随机噪声是提高地震资料信噪比的重要环节。目前在地震资料处理流程中，F-X 域预测滤波是去除随机噪声的有效方法。该方法由Canals Luis 于1984 年提出，经Gulunay［1］、HornBostel［2］等改进后逐渐投入生产使用。国内学者国九英等在文献［3］与文献［4］中分别提出了F-X 域预测与二维滤波并用以及F-X 域预测在三维场景下的使用，使F-X域预测技术进一步趋向成熟。

近年来，为了实现对地下复杂结构的精准勘探，地震资料随机噪声去除技术进一步向着高分辨率方向发展。主要有Oropeza 等［5］提出的多道奇异谱分析，以及Bekara 等［6］把F-X域经验模态分解应用于地震资料随机噪声的去除。以上两者处理非稳态信号的能力都强于F-X 域预测滤波技术，但计算复杂度更高。

随着现代探测设备与新兴处理技术的发展，地震资料处理精度已经能满足对地下复杂结构进行处理的需要，但是高分辨率技术往往意味着大量的计算时间。近十年来，统一计算设备架构（Compute Unified Device Architecture，CUDA）在高性能计算领域发展迅速，并且在地震资料处理领域中应用广泛［7-9］。在F-X域去噪领域，目前关于算法提速的研究较少，其中邴萍萍等［10］提出了改进的F-X 域经验模态分解技术及其分布式并行实现，但是并行的实现基于Matlab平台，在实际应用中具有的意义有限。……