舒幸宁

（长江大学（武汉）地球物理与石油资源学院，湖北 武汉 430100）

大地电磁测深法自20世纪50年代提出以来，以其原理简单、勘探成本低、不受高阻屏蔽、勘探深度大等优点，广泛应用于油气勘探、深部构造、地震监测等领域。大地电磁的源为天然源，信号微弱，极易受到干扰。随着人类文明的发展，大地电磁野外实测的资料中含有的噪声愈发严重，噪声压制研究成为无法回避的问题。从含噪信号中去除噪声并提取有用信号是电磁测深法的一大挑战。采集的资料中往往含有许多脉冲噪声，这些脉冲噪声中含有丰富的频率信息，同时还会向主频两侧泄露，影响频带甚宽，甚至造成信号失真，以致无法得到真实的视电阻率和相位，从而导致无法得到真实的地电模型。张春虹等[1]将AR（p）模型用于解决大地电磁原始资料的不连续性，提高了大地电磁资料的利用率。唐恒专[2]将AR（p）模型应用到地震信号类型的分类与识别中，通过预测误差比作为分类标准，对核爆地震和天然地震进行了有效识别及分类。AR（p）模型是自回归-滑动平均（ARMA）模型的特例，无源激励的AR（p）模型只需要已知前N个数据，就可以得到AR（p）的预测模型。本文将基于均方误差最小的FPE准则用于确定AR（p）模型阶数，采用最小二乘法确定模型参数。本文将AR（p）预测模型引入大地电磁资料处理中，用以压制大地电磁资料中的脉冲噪声，以期提高大地电磁资料的质量。

1 AR模型

对于受脉冲噪声污染的大地电磁信号，我们可以将受突变噪声污染的时间段当做数据缺失处理，对于缺失数据的预测，建立合理的预测模型十分有必要。……