[发明专利]一种频域电磁测深资料二维结构化反演方法

说明书

本发明公开一种频域电磁测深资料二维结构化反演方法，包括：获取沿测线方向的频域水平电场数据；采用以静电场公式进行空间分布处理以模拟异常电场分布；对异常电场中电场异常数据进行分频信息提取；对提取处理后的单频数据进行电场异常信息分类预判；对各预判电场异常做界定范围内的局部最优化反演；外围节点和格架节点进行联合概率成像；将高频、中频和低频频点的反演成像结果进行融合处理；根据多频融合成像结果设计结构化电阻率初始模型；将结构化初始模型提供至二维反演程序，以得到最终电性分布断面成果。本发明的技术方案，以解决电磁场成像对感应源定位不够准确以及感应源对观测数据所造成的畸变影响问题。

技术领域

本发明属于地球物理勘探领域，尤其涉及基于异常体识别技术的频域电磁测深资料二维结构化反演方法。

背景技术

电磁场反演与成像解释方法在近几年都有不同程度的发展，增加了人们对复杂地质构造电磁场的认识。但是，在实际资料的反演与成像过程中仍然面临着诸多问题。比如，近似反演中三维电磁场模拟仍局限在规则且较简单的模型，或者采用一维或二维形式化解释替代方案,否则模拟难度、计算量和成本仍是主要问题。在成像反演中，尽管可以做到空间较连续的成像，但是对数据采集精度和采集密度提出了更高的要求，也使勘探成本大为提高。鉴于这些原因，人们寄希望于开拓新的反演成像方法。

在频域电磁测深实测资料的处理中，对于由浅部电性不均匀体所引起的静态效应，一般都是采取了先校正然后再进行电阻率反演的策略。在技术上，张量阻抗分解技术的应用可以有效地进行校正，并能压制其它干扰(Groom,R.W.,Bailay,R,C，1989，赵国泽、汤吉等1996，魏胜，王家映等1996，F.E.M.Lilley,1998,G.W.McNeice andA.G.Jones 2001)。基于静态位移在剖面上为“高频噪声”的认识，空间滤波和电磁阵列法(EMAP)可以较有效地克服静态位移影响(王家映等1990，1992)。

然而对静态效应的校正，实质上就是人为将客观存在的电性不均匀体给抹掉了，这样处理虽然保证了对断面中、深部电性分布反演结果的相对正确性，但却破坏了方法本身对浅层大地介质分布的高分辨特性。并且，能够引起静态效应的也不仅仅只是浅部不均匀体，模拟实验表明位于断面中、深部的电性不均匀体，同样可以对观测数据造成畸变影响，只不过影响的频段发生了后移，不再表现为简单的全频段平移影响特征！因此，对于这些畸变影响的处理，以往所研发提出的各种静校正方法已不再适用。

发明内容

本发明针对现有技术中电磁场成像对感应源定位不够准确以及感应源对观测数据所造成的畸变影响问题，提出了一种基于异常体识别技术的频域电磁测深资料二维结构化反演方法。

为实现对感应源的精确定位以及有效消除由其引起的对反演结果所带来的不良影响，本发明采用以下技术方案：

一种频域电磁测深资料二维结构化反演方法，包括如下S：

S1、获取沿测线方向的频域水平电场数据，然后将其转换为仅对感应源存在的空间位置敏感的数据形式；

S2、采用以静电场公式进行空间分布处理以模拟异常电场分布；

S3、对异常电场中电场异常数据进行分频信息提取；

S4、对提取处理后的单频数据进行电场异常信息分类预判；

S5、对各预判电场异常做界定范围内的局部最优化反演，获得感应源分布的空间位置；

S6、在确定感应源分布的主体控制格架后，对外围节点和格架节点进行联合概率成像；

S7、将高频、中频和低频频点的上述反演成像结果进行融合处理；

S8、根据多频融合成像结果设计结构化电阻率初始模型；

S9、将结构化初始模型提供至二维反演程序，进行迭代修正，以得到最终的电性分布断面成果。