重力异常的正则化等维反演的研究.pdfVIP

重力异常的正则化等维反演研究 张贵宾晋风明孟令媛 中国地质大学(北京)地球物理与信息技术学院100083 1引言 重力反演是根据观测的重力异常数据求解地下矿体模型。实际上，地下矿体产生的重力异常是全 空间连续分布的。然而，地面重力勘探只能在地表沿有限长剖面或在有限大小的面积上进行离散的 数据观测。所谓的重力反演只能是利用这有限的、离散的测量数据求解地下矿体的密度分布。这就 造成了重力反演结果具有多解性和在深度方向上分辨率降低。 相对于重磁反演而言，某些其它地球物理方法具有较好的深度分辨能力。原因之一是它们属于等 维反演，即，异常数据维数与待求解的地下矿体模型维数相同。例如电测深法除了在地表观测以外， 它还在同一测点上改变极距(直流)或频率(交流)地进行测量以便获得与矿体深度电性分布相关 的一维数据，形成了利用lD／2D／3D异常数据反演1D／2D／3D矿体模型的反演方式。而通常的重磁反演 它具有更严重的低深度分辨率的特点。本文将介绍一种重力正则化等维反演方法，即，利用三维异 常数据反演三维矿体模型。 宾等(1993)将BG理论与吉洪诺夫正则化方法相结合，提出了位场正则化反演方法，应用于二维重 磁反演研究表明正则化方法的使用能够稳定并改善非唯一性的构解。张贵宾等(1996，1997)将该 方法应用于二维自然电位异常反演解释中，分别对意大利Etna火山地区和墨西哥FuegodeC01ima 火山地区的自然电位异常反演结果表明它能有效的探查地下热液活动。Fedi(1999)提出一种有深度 分辨率的三维重磁反演方法，其核心是除了需要地面异常数据以外，还需要数个不同高度平面上的 的反演思想相结合提出了二维重力正则化等维反演方法，即i利用地面和几个高度上的重力数据反 演求解地下模型。研究结果表明等维反演能有效地改善重力反演深度分辨率低的问题，参于反演的 不同高度的异常数据越多，反演精度越高。 本文将张永明(2007)二维重力正则化等维反演方法推广至三维反演，即，三维重力正则化等维 反演方法，目的是在三维情况下研究正则化和等维的结合是否能像二维反演时一样有效。关于等维 反演所需要的三维异常数据可以通过两种方式获得。一是利用延拓方法将地面重力异常数据或航空 重力异常数据延拓到不同高度上，一形成三维异常数据。二是联合使用地面和航空重力异常数据。 2位场正则化线性反演理论 反演方法。在位场反演问题中，观测数据和地下模型之间存在线性泛函关系，并可以通过第一类弗 雷德霍姆方程表示。它是异常和模型之间，通过核函数建立起来的线性关系方程 Ef=IG，(x)聊(x)出 (1) 其中户1，2，…，Ⅳ(异常数据个数)。针对欠定问题，在线性关系方程约束下构制使解的欧几里德 长度取极小的模型，这就把反演问题转化为条件极值问题。将条件极值问题转化为无条件极值问题， 即，求取地下模型脚(力，使得如下目标函数取极小， (2) ①={r聊2(x)w2(x)出+芝哆[弓一eq(x)聊(x)出] _ J=I 其中矿(曲为权函数。通过变分原理，可以获得解模型如下， (3) W—IX)百。。 吣)2南荟哆q(x) 其中ai是如下方程的解向量反， r西=E (4) 其中矩阵r的元素为 11。=r兰兰裂出，(f，v『=1，2，…Ⅳ) (5) 考虑r的特征值分解，可以将方程(4)的解表达成 拈扣鲁砂．．+等式 ㈦