为了获得地下的结构和岩性，需要采集相关的地球物理数据。天然地震记录、人工地震和电磁勘探方法采集得到的地球物理数据含有大量的噪声，通常在获得地下有用的信息前，需要压制原始数据的噪声。目前已经开发了很多方法压制地球物理数据的噪声，如拉动变换、小波滤波、F-K滤波、F-X反折积和中值滤波等。但是，目前已有的这些方法在保持数据边界，有用信号能量恢复等方面均存在自身缺陷。

中科院地质与地球物理研究所“地震波传播与成像”学科组刘伊克研究员最近提出了地球物理资料的矢量中值滤波方法来压制噪声。该方法将地球物理数据看成多维矢量，通过矢量空间的一范数与二范数计算分量之间的相关性和矢量之间的依随性。研究结果表明，地球物理资料的矢量中值滤波方法，通过考虑分量的相关性可以有效地保持断层的边界。图1是一个地震偏移剖面，箭头是通过结构张量计算出的构造走向（event）。为了看清楚中值矢量滤波效果，将白色矩形区域放大（图2），图3是平均方法的滤波结果，可以看到将局部噪声去除，走向形成南北东西2个大方向，但白线划定的不整合面两侧模糊，没有保持断层的边界，图4是中值矢量方法，能很好地保持不整合面的边界。

该方法在噪声压制、裂缝方向和断层方向的拾取、速度分析、地震资料振幅与相位识别、多分量数据计算等方面均有广泛的应用前景。该研究成果近期发表在国际地球物理期刊Geophysics上（Yike Liu. Noise reduction by vector median filtering. Geophysics, 2013, 78(3): V79-V86）。

原文链接

图1 地震偏移剖面，在每一个网格点通过结构张量计算构造走向箭头表示构造走向，总的走向是东西与南北2个大方向，白色框是局部放大区域。

图2  表示图1的放大区域，可以看到构造走向的细节总的走向是南北与东西，但含有噪声。

图3  经过平均方法获得的图2滤波结果

图4  矢量中值方法的滤波结果