王品丰+李强强

摘 要：线性算法的二维大地电磁反演对初始模型的依赖性很强，反演结果的优劣很大程度取决于初始模型，而利用先验信息对初始模型的约束可以有效减少反演的多解性，提高反演精度与可靠性。文章利用MTsoft2D2.4软件，以经典二维模型为研究对象，开展正反演模拟计算，并以先验信息约束初始模型进行反演，得到较理想的效果。

关键词：大地电磁；初始模型；约束反演；二维反演

1 概述

地球物理反演是一个寻找地球物模型，使其响应与实测结果相匹配的过程。其本质是一个最优化的问题，即构建合适的误差目标函数，求取满足精度要求的最小值。地球物理反演本身不可避免的存在着多解性，在解决实际地质问题中，工区的地质、钻孔以及其他物化探先验信息可以帮助构建更合理的初始模型，有效消除假二维反演中的假异常，减少反演的多解性。

2 模型正演

模型参考“全国电法及电磁法勘探正反演软件推荐会（杭州）”组委会提供的二维理论模型，电阻率略有改动。模型（图1）大体为“低-中高-高”三层结构，第一层（0～-4km）模拟表层破碎带，断裂带较多，电性结构较为复杂，高低阻地层互层。中间层（-4～-7km）为中高阻，下层（-7～-12km）是高阻地层，结构较为单一。

模型正演结果（图2）可以看出TE模式的视电阻率结果较TM模式好，出露地表的两个高阻体有所反映，但深部延伸形态反映不够清晰，中高阻层分层情况也不明显；TM模式的左右两边都出现了条带状的“挂面条”现象，这是由于在地表电性突变导致的静位移，是TM模式的明显特征。在相位剖面中，TE和TM模式均可反映出明显的三层结构，但对表层薄层结构分辨不够。因此本文采取TE模式数据进行反演研究。

3 一维反演结果

尽管如今大地电磁的反演算法已经发展到三维反演，但在现今的实测质料处理解释中运用较为成熟的仍是大地电磁的一、二维反演，其中一维反演的结果可以为解释提供依据，更重要的是为二维反演提供初始模型，线性算法的二维反演的结果对初始模型的依赖性较强，这意味着一维反演的结果尤其重要。本文对模型TE数据进行了Bostick、Occam一维反演计算。

图3为模型的一维反演结果，总体来看，Occam反演结果较圆滑，层状反映不明显，浅层构造都基本没有反演出来；Bostick反演结果对左右两个小的斜向的高阻条带反映较为明显，且深部形态也有体现，但在浅部中间位置低高阻相间的层状结构没有反映出来，整体表现为一个低阻体，深部电性界面与实际模型有所差异。

4 约束初始模型及二维反演结果

本文采用稳定高速的NLCG反演，采用不同初始模型进行二维NLCG反演，对比约束初始模型前后的效果。约束反演是在已知地层结构的先验信息的基础上进行的，在反演过程中用先验信息对模型加以约束，以期望反演出符合地质构造的结果。在实际资料的处理中，很多工区都是多方法勘探，包括地震、钻井、重磁勘探等等，可以利用已知资料对大地电磁反演过程中的初始模型加以约束。本文理论模型的反演试算过程中，地层构造都是清楚的，在反演过程中可以根据理论模型对反演的初始模型加以约束。文本通过MTsoft2D2.4软件约束初始模型，先进行模型正演数据的一维反演，查看一维反演的结果与实际模型的差别，然后在一维反演结果的基础上做相应编辑，将修改后的模型作为二维NLCG反演的初始模型，经过迭代计算后得到图4-4中的结果。

图4为不同初始模型的二维NLCG反演结果，其中初始模型為均匀半空间结果比较粗糙；Bostick的二维反演与一维反演结果很相似，只反映出了很粗的结构，甚至掩埋了一维结果中的一些信息，浅部的低阻中层间结构没有反映出来，深部的电性变化的界面也与实际模型有所差异；以一维Occam为初始模型的反演结果中，浅部效果相对较好，界面明显，模型浅部中间段的高低阻相间的结果反映较好，但是深部却与模型差异很大；在一维Occam的反演基础上，利用已经模型约束后的NLCG结果，反演结果与实际模型对应较好，且电性界面清楚，反演效果很好。

5 结束语

二维反演中选择不同的初始模型，二维反演结果差别大。以均匀半空间为初始模型的二维反演结果能反映出模型宏观的特征，仅具有一定参考价值。尽管解决的是二维地质问题，但二维反演依赖初始模型，其结果可能与一维初始模型相似，甚至效果更差。

通过先验信息约束初始模型，可以有效消除反演中的假异常，降低反演多解性。在实际的资料处理解释中，如果有工区的地质地球物理信息，利用已知的信息对初始模型加以约束，然后进行反演计算，可以达到事半功倍的效果。

参考文献

[1]陈乐寿，王光愕.大地电磁测深法[M].地质出版社，1990.

[2]Railway[J].Coal Geology & Exploration，2012.

[3]汤井田，任政勇，化希瑞.地球物理学中的电磁场正演与反演[J].地球物理学进展，2007，22（4）：1181-1194.

[4]柳建新，董孝忠，郭荣文.大地电磁测深法勘探-资料处理、反演与解释[M].北京：科学出版社，2012.

[5]李强强，闵刚.大地电磁中覆盖层对断层正演响应的影响[J].科技创新与应用，2015（33）：51-52.

[6]万汉平.大地电磁测深的TE和TM极化模式对比研究[D].成都理工大学，2010.

[7]王家映.地球物理反演理论[M].北京：高等教育出版社，2012.

作者简介：王品丰（1990-），男，汉族，湖北荆州人，在读硕士研究生，从事电磁法勘探研究。