重力学与固体潮(第3篇).ppt

四. 多解性问题 重力反演结果的可靠性是反演方法是否有效的最关键问题。然而，重力解释中固有的多解性，即无数个不同的场源体可以引起在测量精度范围内相同的异常，严重影响到计算结果的可靠性，使解释者难以从计算得到的许多模型体中选择出适当的模型体表示引起观测异常的地质体。因此，要研究反演，就必须涉及多解性（非唯一性）的问题。 自从Skeels（1947）的著名论文“重力解释中的多解性”发表以来，国外对这个问题有不少论述。本节将就多解性的表现，引起多解性的原因，以及限制多解性影响的方法和途径进行讨论。 【多解性的表现】 Skeels（1947）用一个理想化的理论例子表明了多解性的存在。假定在重力剖面上有17个重力观测值，观测误差为0.1×10-5m/s-2。同时假定，重力异常是真正的二维异常，即重力等值线是一些平行的直线；对观测值已经做过纬度及区域校正，因而这个异常完全是由基底表面的起伏或构造引起。设基底物质密度为2.6g/cm3，在基底上方是密度为2.4g/cm3的均匀的沉积物。根据重力资料，我们能否求出基底顶面的深度及起伏？ 给定在某个深度上的界面形状，用网格法计算这个界面在不同点处的重力效应；根据计算值和观测值的比较，再做一个新的界面；再计算，比较，直到得到一个其计算的重力效应同观测的重力效应的差值小于观测误差的构造为止。 根据上述的反演（正演模拟）方法得到的结果发现，处在不同深度的、具有不同起伏的界面（图2-7-2中的1，2，3…），还包括在这些界面之间的许多界面，都能够引起在测量精度范围内的相同异常。 Skeels认为，计算的重力异常剖面与观测的剖面一致，并不能确保于计算异常对应的模型体就是引起观测异常的地质体。幻想用直接的数学解释方法处理重力资料，以减少解的数目，也是不可能的。因为地球物理场的多解性是固有的，即不同的质量分布能引起相同的场，没有一种数据处理方法能够改变这种情况。 满足所给重力异常剖面的基底起伏的各种解释（引自Skeel，1947） 【实例】 内特尔顿（1987）指出，与一定的异常宽度对应的场源的可能最大深度，就是引起同样宽度异常的点源（球体）的深度。在这个最大深度和地面之间，存在一个可能源的锥形区。同一个异常可以有埋藏很浅的薄透镜体，或不太狭窄的较厚的物体，或球体引起；而且在球体上方较浅的深度上，却有无限多个可能的场源引起这个相同的异常。引起同一异常的不同的场源，有一个唯一的共同的特性，就是它们的剩余质量必须是相同的。 引起相同异常的可能源的锥形区（引自内特尔顿，1987） 【练习】怎样解释中心埋深相同，剩余质量相同的大球与小球产生相同的重力异常？ 由球体正演公式可以看出，重力异常与球体的剩余质量有关，不同半径，不同剩余密度的组合都可以得到相同的剩余质量，因此说，埋深相同的大球与小球（剩余质量一样）产生相同的异常，无法唯一确定，这也是重力反演的多解性。 第三节 密度分界面重力异常的正反问题 一、密度界面的确定 剩余密度σ＝σ1－σ2 剩余密度σ‘＝σ2－σ1 σ1 σ2 σ’＝σ2－σ1 σ＝σ1－σ2 以地面为起算点，逐层计算各密度界面的影响： 1、T3界面： σ＝σ1－σ2；abcd 2、T4界面： σ＝σ2－σ3；cdef 3、T5界面： σ＝σ3－σ4；efgh c d e f g h 二、多个密度分界面重力异常的计算 1、计算这两个密度界面引起的重力异常(图7－31(a))相当于计算两个界面 2、计算这两个密度界面重力效应的方法是把 间的密度层分为两层(图7－32)。 在实际的地质条件下，常出现在某一地区缺失某个地层的情况，即出现“尖灭”现象，对于这种情况的重力异常，也不难计算。 三、一个密度界面的正反演 σ1 σ2 线性公式反演方法 求S界面的深度 对于A点 在已知B点处深度HB时，就可计算出A点处的深度HA 2、方域计算法 在计算机上进行数值计算，采用方域网计算显得十分方便。 将测区以间隔为△x、△y的两组直线分成一系列小方块，每个结点就代表一个横截面为△x●△y，深度为Hij的直立矩形柱体。 将计算点作为坐标原点，计算范围取： X: -m △x~m △x; Y: -n △y~n △y； z: 0~h, h=h0+ △h 则在直角坐标系中有： 在当△h=h-h0＜＜h0时，对被积函数用泰勒级数展开，忽略△h二次以上的各项，可以得到： （四）三维密度界面的快速正演方法原理 单一密度界面 双密度界面 （五）变密度多界面快速反演方法 六、矩阵法 1、模型的正演计算 四 矩阵法 重力异常反演 一、计算地质体几何参数和物性参数的特征点法 特征点法（或任意点法）是根