第五章 复杂形体重、磁异常正演问题综述.pdf

第32卷第1期 物探与化探 V01．32．No．1 2008年2月 GEOPHYSICAL＆GEOCHEMICALEXPLORA’n0N Feb．．2008 复杂形体重、磁异常正演问题综述 李焓，邱之云，王万银 (长安大学地质工程与测绘工程学院，陕西西安710054) 摘要：总结了复杂形体重、磁异常各种正演方法，并对各种方法所具备的优点和存在的缺点进行了分析说明。据此 讨论了复杂形体重、磁异常正演方法将来的发展方向。 关键词：复杂形体；重力异常；磁力异常；正演计算 中圈分类号：P631 文献标识码：A 文章编号：1000一8918(2008)01一0036—08 复杂形体重、磁异常正演是重、磁勘探中的一项 1．1．1．1点元法¨“o 重要工作。在勘探阶段，为了判断有无盲矿体的存 如图1所示，将任意形体用3组平行于直角坐 在，必须计算形体比较复杂的已知矿体所能产生的 标面的截平面进行分割，使物体分成许多具有规则 异常值，并与实测异常值进行比较而给予确定。正 形状的长方体元。用解析方法计算出所有这些长方 演问题是已知物性和形态求异常值，而反演问题是 体在计算点产生的异常，并累加求和，就可以得到整 已知异常值求物性或形态，所以正演问题是反演问 个形体在计算点引起的异常值。 题的基础。国内外许多学者都曾对不同形状的二度 体和三度体的重、磁异常正演问题进行了研究，并给 出了相应的计算公式，每种方法均有各自的应用条 件，也有各自的优越性。 复杂形体重、磁异常正演从形体角度出发可以 分为三度体和二度体；按求解域还可分为空间域和 (a)组合模型 (b)任一单元模型 频率域。笔者在查阅前人研究成果的基础上，对各 图l三度体点兀法示惹 种正演方法进行了归类和总结，并给出了进一步研 当长方体顶面的4个角点在计算平面上的垂直 究的方向。 投影点与网格点重合时，原重、磁异常表达式会出现 1三度体正演问题 分子、分母同时为零的情况，因而这些网格点成为无 法计算的解析“奇点”。郭志宏在文献[4]一文中对 1．1空间域正演方法 该问题作了详细的研究，并给出正确的表达式： 复杂形体重、磁异常正演在空间域中的计算又 △丁(戈，)，，z)= 可分为两大类：有限单元法¨娟1和边界单元 法一一引。空间域正演方法的基本优点是适用于任 意测网。 无31n[r+(f一彳)]+ 1．1．1有限单元法 蜘rctan瓦寻等等编+ 本类方法的基本思想是用不同的方式对复杂形 体进行分割，使之转化为一系列简单形体的组合，计 七sarctan 算该简单形体的重、磁异常再求和即可得整个复杂 iji：_歹i；‘：亏ji；{；：；编+ 形体的异常。这就包括了2个方面的近似：一是分 ‰arctan 割方式与实际形体的拟合程度；二是数值积分代替 g专掣《：：：!：!_～，c，， 解析积分的近似程度。 收稿日期：2006一12—14；修回日期：2007—01—10 万方数据 _L～————翌墼竺兰兰望望竺兰