Parker-Oldenburg-方法.pdfVIP

Parker-Oldenburg 方法 本文旨在简单对界面反演方法中的经典方法：Parker-Oldenburg 方法做一个总结, 以方便自己参考使用。更多的拓展和细节还是老老 实实以文献为主吧。 Parker-Oldenburg 方法起源于Parker （1973）年提出的界面势场 异常频率域快速正演公式，之后由他的学生Oldenburg （1974）在它 的基础上提出了反演的方法。与之前的方法相比（如压缩质面方法等）， 这种方法由于采用FFT （快速傅里叶变换），计算速度非常快，并且在 一定的条件下，计算的收敛性也可满足，因此成为了势场界面反演的 经典方法。 1. Introduction 首先简单介绍一下重力界面反演问题：假设地下有一个密度分界 面，分界面上下介质有密度差，我们在地面上获取到重力异常Δg, 并 且已知密度界面上下的密度分布（通过其他地球物理场来获取）。我 们希望通过这些信息来获得界面的形状。 下面的讨论中我们假设界面的平均深度为z0 ，界面相对于平均深 度的起伏为h(x, y) ，并假设上下界面的密度差为ρ(x, y) 。并假设相对 界面的起伏总在界面的一侧，且起伏不大。 在直角坐标系中，坐标轴的选取如图1 所示，以地面上某一点O 作为坐标原点，x-y 平面置于水平,z 轴垂直向下, 即沿重力方向。场源 坐标系和密度界面图示分别见图1 和图2 。 图1.场源坐标系 图2. 密度界面示意图 2. Parker 正演公式 对于地下的某一密度分布ρ(x, y),在地面(水平面，z = 0)上观测到 的重力异常Δg(x, y, 0)可以表达为： ( ) Δg = ∭ Gρ x′, y′ z′ 3 dx′dy′dz′ [(x′ − x)2 + (y′ − y)2 + z′2]2 这里我们认为h 0 , 即起伏界面在z0 下方，因此积分区域为： ( ) ( ) x′: −∞, +∞ , y′ : −∞, +∞ , z′: (z , z + h) 。 0 0 对Δg 中的x, y 变量做二维傅里叶变换： (+∞)2 FT(Δg) = ∬ Δg e−iux−ivy (−∞)2 z0 +ℎ 1 ( ′ ′) ′ ′ ′ ′ = G ∬ ρ x , y dx dy ∫ FT[(( − ′)2 + ( − ′)2 + ′2)3/2] z0 上式中右端的FT[.]一项由傅里叶变换的延迟定理可知，它等于： 1 ′ ′ FT [ 2 3 ] − − ( 2 + 2 + ′ )2 所以： ′ ′ z0 +ℎ 1 ( ) ( ′ ′) − − ′ ′ ′ ′ FT Δg = G ∬ ρ x , y dx