稳健数据平滑在电磁测深数据处理中应用.pdfVIP

第33卷增刊 物探化探计算技术 2011年7月 文章编号：lool—1749(2011)增刊(1)003704 稳健数据平滑在电磁测深数据处理中的应用 化希瑞1’2 (1．中铁第四勘察设计院集团有限公司，武汉430063； 2．中南大学地球科学与信息物理学院，长沙410083) 摘要：电磁类地球物理方法由于极易受到各类噪声的干扰，使得估计的视电阻率曲线或相位 曲线发生畸变，严重地影响了反演解译精度，如何对这类观测曲线进行合理有效地平滑，是目 前数据处理中的重点，也是难点。作者在本文采用自适应robust平滑方法处理电磁测深曲线， 该方法基于惩罚最小二乘估计，利用广义交叉验证估算平滑因子，通过robust估计有效地压制了 噪声与“飞点”的干扰，恢复了曲线的基本形态，有效地保持了曲线基本特征。 关键词：数据平滑；惩罚最小二乘；稳健；数据处理 631．3+25 中图分类号：P 文献标识码：A 自动robust平滑方法对电磁测深曲线进行平滑处 0前言 理，有效地压制了噪声及“飞点”的干扰，恢复了 曲线的基本形态，实验结果显示了方法的有效性。 频率域电磁测深法由于工作效率高、分辨能力 强、探测深度范围大等特点，在固体地球物理、深部 1 robust平滑理论 地球物理、工程地球物理及石油、煤田资源探测等 应用中勘探效果良好。然而由于众所周知的原因， 1．1惩罚最小二乘回归 电磁类方法极易受到各类噪声的干扰，使得估计 假设数据初始模型如下： 的视电阻率曲线或相位曲线发生畸变，严重地影响 Y(t)=儿(t)+占(t) (1) 了反演解译精度。其中最常见的现象之一就是观 其中 Y(t)为观测数据；Y0(t)为一个光滑曲线； 测曲线的不平滑、不连续，甚至有单个“飞点” 占(￡)为白噪声过程，根据混杂有白噪声的观察资 料Y；，i=1，…，n去重构光滑函数g，这在统计学 (outliers)或连续“飞点”存在，如何对这类曲线进 行合理有效地平滑，是目前数据预处理的重点旧j。 中叫非参数回归，在数字信号处理中则叫信噪分 常用的各类插值平滑方法中由于平滑效果对于窗 离，非参数回归方法有局部线性估计、移动平均滤 口大小选择极为敏感，不仅不同曲线需选用不同的 波、近邻估计、多项式样条估计和惩罚最小二乘估 窗口，即使在同一曲线中，由于噪声类型或强度不 计等，其中惩罚最／J、--乘回归，是在最／b乘法的 同，在不同的频段也需选择不同的窗口，因此这类 目标函数上增加一个非参数部分的惩罚项(P)，即 方法平滑效果及效率极为有限。电磁探测作为一 求解如下条件极值问题【l。4】： 种体积探测方法，不同频率的探测结果应是一条光 RSS+s尸(儿)=0儿一Y0 2+ 滑的曲线，曲线的极值点、拐点或者说曲线的特征 sP(y,)=rain (2) 点反映了观测对象性质的变化。噪声的存在往往 其中 ”0是欧氏范数；3为平滑因子，Weinert【2】 会影响这些特征点的确定，平滑的目的就是为了压 提出了惩罚函数(P)简单表达的形式： 制噪声干扰并增强特征对比度。作者在本文利用 P(儿)=0or,112 (3) 收稿日期：2011—06—07 38 物探化探计算技术