SinoProbe大地电磁测深阵列观测资料的反演研究.pdfVIP

中国地质学会 2011 年学术年会 SinoProbe 大地电磁测深阵列观测资料 的反演研究 1，2，3）1 1，2，3） 1，2，3） 1，2，3） 1，2，3） 1，2，3） 董浩 ， 魏文博 ， 叶高峰 ， 金胜 ， 张乐天 ， 张帆 ， 谢成良 1，2，3） 1，2，3） ， 尹耀田 ， 以及 Sinoprobe01 项目组成员 1） 中国地质大学地下信息探测技术与仪器教育部重点实验室, 北京，100083； 2 ） 地质过程与矿产资源国家重点实验室， 北京， 100083 ； 3 ） 中国地质大学地球物理与信息技术学院, 北京， 100083 关键词： Sinoprobe ； 大地电磁测深； 反演方法 为了深入了解地壳和上地幔的导电性结构，SinoProbe 专项 01 项目在华北布设了 1°×1° 网度的大地电磁测深（MT ）阵列式示范观测，这包括约 1400 个大地电磁宽频和长周期测深点的 数据。显然，以往针对 MT 长剖面探测数据的二维反演方法已不能简单套用于面积性的 MT 阵列 数据反演。 大地电磁场二维反演是目前常用的 MT 反演方法。进行 MT 二维反演的首要问题是地下导电 结构的维性和电性主轴方向的判定，以及剖面方向的选择。但是，对于大区域的阵列观测，当截 取的 MT 长剖面跨越多个地质构造单元时，或是研究区地质构造复杂多变，电性主轴方向的分析 结果往往很难有比较一致的趋势，这给 MT 二维反演带来很大的困难。而对于面积性的 MT 阵列 观测，获得的是数据量巨大的 MT 三维数据体；如果考虑到地下介质的几何结构或导电性结构通 常具有很强的三维性，那么针对 MT 阵列观测数据体采用三维反演方法将是最合理的选择；但是， 就目前所具备的计算条件而言，MT 三维反演计算需要耗费大量计算机时间和内存量；由于 MT 阵列观测数据体过于庞大，即使运用我们现有的计算机服务器平台(HP BL460c G7)和并行算法， 所能达到的 MT 三维模型剖分网度和数据密度与二维反演相比仍极为有限，这可能导致三维反演 模型的分辨率降低，达不到预期效果。实际上，对于如此巨大数据量的 MT 阵列探测，在中国大 陆还是前所未有的，必需研究在现有计算条件下，如何实施MT 阵列观测数据体的三维反演？同 时研究MT 二维反演方法对于阵列观测数据反演的适用性。这无疑对于 SinoProbe-01 项目的成功 至关重要。 为此，我们首先开展理论研究，根据对实际地质构造条件的物性分析，设计一系列电性结构 模型进行三维正演计算，并把计算结果叠加 5%噪声数据用以模拟 MT 实测数据，构成阵列式观 测模拟数据体。利用这模拟数据体截取经度和纬度方向，以及垂直理论模型主构造方向的 MT 长 剖面数据，进行 MT 二维反演并分析反演模型之间的差异，以及与所设计的正演模型之间的近似 1 作者简介： 董浩，男，1984 年生。中国地质大学（北京）地球物理专业在读博士生。主要从事大地电磁测深 与反演技术方向的研究工作。电话：010； E-mail: donghao.cugb@ 。 38 中国地质学会 2011 年学术年会 程度。 此外，利用这模拟数据体研究了三维模型不同剖分网度和不同站点、频点密度的三维反演计 算时间、内存量和三维反演模型与理论模型之间的近似程度。在此基础上，对鄂尔多斯地块的 MT 阵列观测数据进行了三维反演。 研究结果表明，在选择剖面时，当电性主轴方向接近南北或东西方向的情况下，截取纬度或 经度方向的剖面应与主构造垂直方向相差不大，因此其二维反演结果都可以较好的还原正演模 型。在多数情况下，正演模型的主要结构都可以很好的被二维和三维反演解析出来。在现有的计 算条件下，由于三维反演使用的模型网格数与数据量的限制，二维的反演结果有时可能会比三维 的反演结果的边界更清晰，更精确。然而，即使经过维度分析和畸