三维大地电磁测：数据空间法

三维大地电磁反演：数据空间法

摘要目前，一种三维大地电磁最小模型反演算法已经提出，这种算法是奥卡姆反演方法的变种，其主要是思想是基于数据空间的反演算法。由于模型空间矩阵的计算时间相对较长，使得基于模型空间的奥卡姆法三维大地电磁反演并不实用。这些困难能够用基于数据空间的奥卡姆反演算法来解决，在这种方法中，矩阵维数依赖于数据空间的大小，而不是模型空间的取值。通过将模型空间转换到数据空间，从而使得奥卡姆方法能够在PC机上进行反演计算。为了减小计算时间，一种宽松的集中规则被用于计算灵敏度矩阵的迭代正演模型标准。这种规则使得计算时间压缩了70%，且不影响反演结果。通过模拟数据的数值计算试验表明通过少量的迭代次数中就能够得到满意的结果，在之后的迭代中需要取消不必要的结构并找到最小标准的地电模型。

关键词大地电磁；数据空间法；三维反演；奥卡姆反演

1引言

进行三维大地电磁反演的常规运算是对大地电磁方法的未来发展的需求，由于二维解释常常并不能解释复杂的地质区域中场数据集呈现的重要特征。近年来一些人在发展三维大地电磁反演算法做出努力，使用了合理的大范围趋近方法（例如：Mackie 和Madden在1993年,Newman和Alumbaugh在2000年，Farquharon等在2002年）。这些方法已经能够去合理的恢复电导率变化，至少在理论数据的案例测试中得到验证。然而，三维大地电磁反演问题还远没有解决。高配置终端工作站或者并行计算机的需求依然阻碍三维程序运行的应用，计算机问题与实际数据的真实性影响着所有的设想方法。所以，提高实施三维反演算法的效率受到了高度的关注。

基于快速相似模型的规则有望提高效率，比如你拟线性或拟解析相似模型（Torress-Verdin和Habashy在1994年，Zhdanov和方在1996上半年，Tseng等在2003年）。因为这些相似性模型响应是正比于修改的电导率张量拟线性函数，是可能

1 序言

第 - 2 - 页 （共 16 页）

的简化反演方法（有如：Zhdanov 和方在1996下半年，Zhdanov 等在2000下半年）。依据被使用的简化程度，像这样的方法能够快速对地球结构成像。但是，这些相似性方法也有各自的局限性，比如，当电导率相对较小时效果最好，而且一般反演方法的真实性和准确性对所有问题都比较通用。尽管这些快速反演规则没有问题的求解，但这些基于电磁感应方程的所有解答方法依然有诸多用途。

在大部分的三维大地电磁反演方法采用了更多经典的反演接近方式中；有的是将函数和数据的误差与粗糙模型链接起来（像Parker 在1994年），有的是将被用于计算的惩罚函数梯度最小化，适应模型规则，并用满足模型的趋近方式数字法解解决相关的正演问题。例如，Mackie(2002年在个性化交流中)已经将非线性共轭梯度方法向三维扩展（NLCG,Rodi 和Mackie 在2001年），在PC 机组的信息通行界面（MPI ）的辅助下，Newman 和Alumbaugh(2000年)就向三维反演使用了类似的技术，不过在大量平行系统的帮助下，Sasaki(2001年)和Farquharson 等人（2002年）就根据高斯-牛顿（GN ）发展了三维反演方法。所有的这些规则能够成成模型空间反演，并于此研究优化M 维模型空间电导率的形成。

用模型空间方法，方程的数量和可使用内存（RAM ）都强烈的依赖于模型组M 的大小。特别地，在模型空间中最直接最小接近惩罚函数（比如高斯-牛顿法）时要生成并解决M M ⨯组的常规线性方程。对于三维反演，M 的变大在关系到计算时间和特别的是可写入内存使得计算机不可能运行。关于解决这个困难的一种近似方法（即，Sssaki 在2001年，Newman 和Alumbaugh 在2000年）是使用一种电导率变化的粗糙或原始方案，使得M 依然较小。但是正演的结果却强烈的依赖于模型方案的选择，除非地球结构有非常强的优先限制，否则这样的计算结果可能会误导我们。对于经过迭代的反演方法，比如共轭梯度（CG,Mackie 和Madden 在1993年）或者非线性共轭梯度（Newman 和Alumbaugh 在2000年）避免了过于严格的形式并为常规方程储存

M M ⨯阶系数矩阵，

而且这种方法能够容许更多一般的地质真实模型规则。这种一般的迭代趋近方法在只有为已有规则的最小三维反演结构的基础上实际计算规则而且M 较大时逐渐引起重视。

这篇文章里，我们引出一种新兴的数据空间三维大地电磁反演算法，并将惩罚函数的最小值在N 维数据空间中形成。并用数据空间方法趋近原始空间（即：部分范围空间在数据中没有影响）在输出时被取消，同时M M ⨯组常规方程被N N ⨯组方程

三维大地电磁反演：数据空间法

第 - 3 - 页 （共 16 页）

取代。这样一来，所有计算的数量和需要的数组首要的依赖于独立数据N 的大小，这对三维地质仿真模型将远远小于M 。实际上数据空间方法被广泛的应用于反演地质问题（即：Parker 在1994年）和其他的物理场（Egbert 等在1994，Chua 和bennett 在2001年）。如果没有特殊的限制，数据空间趋近法不考虑共轭梯度法而考虑反演算法。我们准确认为这类趋近方法是在奥卡姆准则上的数据空间转化。原始的奥卡姆反演（Constable 等于1987年，deGroot-Hedlin 和Constable 于1990年）方法是在模型空间列出方程。现在，我们重新列出在数据空间的奥卡姆方程等式就像在Siripunraraporn 和Egbert(2000)中解决的二维大地电磁反演问题。如同在Siripunvaraporn 和Egbert(2000),数据空间算法利用减小的基本接近更加有效。在这里，我们只考虑大部分直接在奥卡姆准则下的数据空间转化执行情况，并作为数据空间趋近中实际意义的一个实例。

2 奥卡姆反演

奥卡姆反演寻找最平滑或者最小的标准模型满足数据（Constable 于1987）。在数学上，只要能找到常规泛函数(,)U m λ:

{}1112*00(,)()()([])([])T T m d U m m m m m d F m d F m C C X λλ---=--+--- (1)

这里m 是电阻率模型, 0m 是初始模型，m C 是模型的协方差矩阵定义的模型标准，d 是已经观测的数据，[]F m 是模型响应，d C 是数据响应矩阵，X *是期望误差，1λ-是拉格朗日系数。假设数据协方差特别的准确，X *应该是理论上的1（Constable 等于1987），同时我们在模拟数据实验上使用这个值。对于真实数据*X 在一般情况下必须设置一个高一点的值。

为了找到方程（1）中的驻点，为了不直接求解（1），我们对求惩罚函数()W m λ。当λ固定时，U 和W λ都有同样的驻点。

{}11100()()()([])([])T T m d W m m m m m d F m d F m C C λλ---=--+-- (2)

用一系列λ使W λ最小，则能够得到U 的驻点值。（例如，找到λ使得数据误差是2*X ）。

2.1 模型空间方法

由于是大地电磁反演非线性问题，需要采取迭代法逼近，考虑到线性函数[]F m 是

11[][][]()k k k k k k F m F m m F m j m m ++=+∆=+- (3)