论文部分内容阅读
音频大地电磁(AMT)技术在20世纪逐渐发展成一种重要的地球物理勘探手段。我国自“九五”开始从北美地区引进了数量较多的该类型测量设备,美国Geometric与Electromagnetic Instruments(EMI)所生产的STRATAGEM音频大地电磁(AMT)测量系统(俗称EH4)即是其中重要的一种设备,并且在国内AMT市场上有一个相当大的用户群。但是,所引进的STRATAGEM测量系统的阻抗张量估算技术仍停留在大地电磁资料处理技术的早期阶段(采用方法为最小二乘法),与AMT测量大范围应用及大地电磁资料处理技术的发展不相称。本文剖析了STRATAGEM音频大地电磁测量系统的数据组织与资料处理流程;从常相干度参数入手,剖析了系统的IMAGEM软件处理时序资料谱的筛选、叠加等关键步骤,在此基础上分析了该软件处理时序文件与阻抗估算的不足。大地电磁阻抗资料稳健估算是当前该领域主流的处理方法之一。最大似然估算(即M估算)是一种重要的稳健(Robust)处理方法。本文在介绍最大似然估算方法的基础上,借鉴前人的大地电磁阻抗稳健估算算法并融合STRATAGEM的时序文件数据结构,应用胡倍尔(Huber)提出的最大似然估算的损失函数、影响函数构建了具有稳健(Robust)估算大地电磁阻抗功能的TSP程序。在TSP的M估算实现过程中,其尺度参数采用以残差均方根估算与帽子矩阵融合的方法求得。最后,通过实例给出了由TSP与IMAGEM处理的阻抗曲线;对比显示,TSP处理结果抗干扰能力明显增强。