重力异常正反演问题
球体重力异常正演程序报告
球体重力异常正演程序报告球体重力异常正演是地球物理学中的一种重要方法，用于研究地下物质分布和地球内部结构。
本报告将重点介绍球体重力异常正演程序的原理、步骤和应用。
一、原理球体重力异常正演是基于牛顿引力定律和球体模型的数学计算方法。
根据牛顿引力定律，在球体表面上的任意一点，重力加速度可以表示为：g = G * (M / r^2)其中，g为重力加速度，G为引力常数，M为球体的质量，r为球心到该点的距离。
根据球体模型，球体的质量可以表示为：M = (4/3) * π * ρ * R^3其中，ρ为球体的密度，R为球体的半径。
将质量公式代入重力加速度公式，可得到球体表面上的重力加速度公式：g = (4/3) * G * π * ρ * R / r^2二、步骤球体重力异常正演程序的步骤如下：1. 确定观测点的位置和高度，以及球体模型的半径和密度。
2. 计算球体表面上的重力加速度，根据上述公式进行计算。
3. 根据观测点与球心的距离，计算球体表面上的重力加速度的投影值。
4. 重复步骤3，直到计算出所有观测点的重力加速度投影值。
5. 计算观测点的球体重力异常值，即观测点的重力加速度减去球体表面上的重力加速度投影值。
三、应用球体重力异常正演程序在地球物理勘探中有广泛的应用，主要包括以下几个方面：1. 地质勘探：通过球体重力异常正演，可以对地下的岩石密度分布进行推测，从而帮助地质勘探人员确定地质构造和找到潜在的矿产资源。
2. 油气勘探：油气藏通常与地下的密度异常有关，通过球体重力异常正演，可以对潜在的油气藏进行初步判断，指导油气勘探的方向和深度。
3. 地壳构造研究：地球内部的构造和演化与地下岩石的密度分布密切相关，通过球体重力异常正演，可以揭示地壳的变形和演化过程，为地壳构造研究提供重要的参考依据。
4. 火山和地震研究：火山和地震活动通常与地下的岩浆和断层有关，球体重力异常正演可以帮助科学家们理解火山和地震的发生机制，预测可能的灾害风险。
重力异常正演资料
单位长度， dm d
• 若水平圆柱体有限长， 则
密度均匀的水平圆柱体
l
g G
d
l [( x)2 ( y)2]3/2
(x
2Gh0l
h0 )(x2 l2
h02 )3/2
密度均匀的水平圆柱体
• 当 l 时，
g 2Gh0
密度均匀的球体Vg VFra bibliotekzG
v
( z)d dd [( x)2 ( y)2 ( z)2 ]3/2
密度均匀的球体
密度均匀的球体
Vg
GM
[x2
h0 y2
h02 ]3/2
密度均匀的球体
Vg
GMh0 ( x2 h02 )3/2
球体重力异常图
球体重力异常图
利用已知异常计算球体参数
重力异常正演
正问题与反问题
正问题也称为正演计算(Forward Calculation) 已知地质体的形状、产状和剩余密度等，通过理 论计算来求得异常的分布和规律。
正问题与反问题
• 反问题也称为反演(Inversion) • 已知异常的分布特征和变化规律，求场源的赋存
状态（如产状、形状和剩余密度等）
正问题与反问题
正演计算是解反问题的基础，解反 问题（反演）是目的
正问题与反问题
简单规则几何形体的异常
• 为了简化，假设地质形体孤立存在，密度均匀， 地面水平，所取剖面为中心剖面。
• 规则形体：球体、水平圆柱体、垂直台阶、脉状 体……
密度均匀的球体(点质量)
• 自然界中，一些近于等轴状的地质体， 如矿巢、矿囊、岩株、穹窿构造等， 都可以近似当作球体来计算它们的重 力异常，特别当地质体的水平尺寸小 于它的埋藏深度时，效果更好。
第一章重力勘探：5重力异常的推断解释
导数等
（1）重力异常的解析延拓
向上延拓： 将观测平面上的实测异常值，换算到观测平
面以上某一高度上的异常——称为向上延拓。 目的：消弱局部异常，突出深部异常
向下延拓： 将观测平面上的实测异常值，换算到观测平
面以下场源以外的某个深度上——称为向下延拓。 目的：压制深部的区域异常，突出浅部物质
产生的局部异常
地质解释时包括数量 方面的内容。如对异 常源的大小、产状、 埋藏深度给出具体解 答
正演问题： 已知地质体的形状、产状、物性参数，求场
（异常）的分布；
反演问题： 已知场（异常）的分布特征及变化规律，求
场源的赋存状态（如产状、物性参数、埋深等）
须指出的是：
正演问题的解是唯一的，而反演问题的解， 则具有“多解性”
A•
（二）重力异常的划分
所谓异常的划分，就是从迭加场中将区域场与 局部异常分离开来。
1、图解法（异常曲线平滑法）
布格异常平面等值线图
局部异常平异常平面等值线图
2、圆周平均法（多边形法）
r
3、位场转换
什么是位场转换？ 根据观测平面上的实测异常值，利用数学
的方法求场的导数或任意高度场的分布。 通常有解析延拓（上延和下延）、求高次
一、重力异常的识别和划分 （一）重力异常的识别
1、区域异常与局部异常的概念 区域异常 场源：大而深的岩体或地质构造
异常特征：幅值大、异常范围大、变化平缓 场源：小而浅的岩体、矿体或地质构造 局部异常 异常特征：幅值小、异常范围小、变化大
A B C
2、局部重力异常的识别 局部重力高—异常等值线在局部范围内向区域背景场 降低的方向凸出。
向下延拓
向上延拓
分
析
重力异常正反演问题
整个水平圆柱体在P点产生的重力异常为无 穷多个柱体在该点产生的重力异常之和，即：
g G h0
2 h0 dy G 2 2 2 2 3/ 2 2 ( x y h0 ) x h0
2、水平圆柱体
规 则 形 体 的 正 、 反 演 问 题
讨论：
3、垂直台阶
2πG△σ△h
由图可见：无论台阶产 状如何，异常的形态相 似，仅原点处的异常值 不同。 当台阶直立时：
P(x,0)
πG△σ△h
o
△σ
●
x
△g(0)= πG△σ△h
△h
当台阶面向台阶外侧倾 斜时：
△g(0) > πG△σ△h
当台阶面向台阶内侧倾 斜时：
△g(0) < πG△σ△h
二、不规则三度体的正演问题
（4-3）
如果还有其他的物性层界面存在，则可仿照以上公式（4-3）进行迭加，以 求多重界面的Δ g。将Δ g 进行傅立叶变换，便得空间域的Δ g(x,y,0)，即完成 正演计算。
双密度界面
（五）变密度多界面快速反演方法
1．反演问题的基本原理
现假设已知重力场为△g（r0,z0） ，其频谱记为 F[△g]，又假设已知密度函数ρ (r)的一个值为 ρ
重力异常的正反演
• 1. 重力异常的正演问题、反演问题； • 2. 均匀密度球体、水平圆柱体、台阶的重 力异常正演方法，异常特征，反演方法； • 3. 密度界面的剩余密度的确定方法； • 4. 单一密度界面异常的特征及反演解释方 法（近似解法、矩阵法）； • 5. 解复杂密度体正演问题的基本思想； • 6. 最优化选择法的基本思想；
用解析公式计算出每个小长方
反演问题正断层逆断层3
V XZ 3GM
Dx (x 2 D 2 )5/2
VZZ
2D2 x 2 GM 2 (x D 2 )5/2
2D2 3x2 3GMD 2 (x D 2 )7/2
VZZZ
2、水平圆柱体
规 则 形 体 的 正 、 反 演 问 题
D
讨论：
g 2GD x2 D2
重力作业（5月9日交）
• 1.请解释重力异常的实质. • 2.岩矿石的密度有哪些特征? • 3.画出球体重力异常的剖面特征与平面特征，它 与水平圆柱体重力异常有何不同？ • 4.什么是相对布格重力异常？写出其表达式. • 5.什么是重力异常的解析延拓?向上与向下延拓各 有什么作用? • 6.什么是重力的导数法?重力高次导数有什么作用? • 7.举例说明重力勘探的应用.
3、垂直台阶
规 则 形 体 的 正 、 反 演 问 题
不同埋深的台阶剖面（a）和铅垂台阶的Vxz、Vzz、Vzzz
断层或不同岩层的接触带都可作为台阶处理
x2 h12 1 x 1 x g G (h x ln 2 2h1tg 2h2tg ) 2 x h2 h1 h2
由公式可见： 当 x 时，
g 2G h
2G h
G h
球体参数：半径30m，中心埋深50m，密度 2.0g/cm3，重力异常单位g.u.
1、球体
规 则 形 体 的 正 、 反 演 问 题
4 3.5 3 2.5 2 1.5 1 0.5
x 10
-3
P(x,0)
-40 -30 -20 -10 0
0 -50
D
10
20
30
40
50
讨论：
简述重力场的正反演问题
简述重力场的正反演问题
重力场的正反演问题涉及重力异常的正演和反演。
正演问题是给定地下某种地质体的形状、产状和剩余密度等，通过理论计算来求得它在地面上产生的异常大小、特征和变化规律，这是正向思维的问题。
反演问题则是依据已获得的异常特征、数值大小、分布情形等并结合物性资料来求解地下地质体的形状和空间位置等，这是逆向思维的问题。
重力正演是指根据地下地质体的形状、大小、密度等物理参数，利用重力场理论计算其在地球表面产生的重力异常。
重力反演则是根据实测的重力异常数据，结合物性资料，推断地下地质体的形状、大小、空间位置等信息。
重力正演是解决正问题的过程，它从地下地质体的物理参数出发，预测其在地球表面产生的重力异常。
重力反演则是解决反问题的过程，它从实测的重力异常数据出发，推断地下地质体的形状、大小、空间位置等信息。
重力场的正反演问题在地球物理学中具有重要的应用价值，例如在矿产资源勘探、地质构造研究、地下水资源调查等领域都有广泛的应用。
通过正反演问题的解决，可以更好地理解地球内部结构和动力学过程，为资源开发和环境保护提供科学依据。
重力正演、反演
2)当σ>o时，极大值一侧对应着上升盘，极小 值一侧对应着下降盘，在极小值十分清晰且大 干极大值的绝对值时，属正断层类型，反之则 属逆断层类型。
二度铅垂柱体 对于沿水平方向延伸较长而横截面近于矩形的 矿脉，可以当成二度铅垂柱体来研究。在正演 它的异常时，坐标系及有关参数的选取见图,用 (x+α)与(x一α)分别代替铅垂台阶各公式中的 x，并将结果相减，即获得这一形体的重力异 常及各阶导数异常的公式:
当柱体的下底 H→+∞ 时，便可获得底部无限延 伸的铅垂脉的相应公式Δg→∞
( x − a) 2 + h 2 V xz = Gσ ln ( x + a) 2 + h 2 h h 2ah V zz = 2Gσ (tg −1 − tg −1 ) = 2Gσtg −1 2 x−a x+a x + h2 − a2 ⎡ ⎤ x+a x−a 2a ( a 2 + h 2 − x 2 ) V zzz = 2Gσ ⎢ = 2Gσ 2 − 2 2 2 2 ⎥ ( x + a) + h ⎦ ( x + a 2 + h 2 ) 2 − 4a 2 x 2 ⎣ ( x + a) + h
GM GMD = 2 2 nD ( x1 / n + D 2 ) 3 / 2
x 1/n = ± D n 2 / 3 − 1
取n=2，得x1/2=0.766D(X正半轴)和x’1/2=－0.766 D (X负半轴)，说明异常半极值点的横坐标为球心 深的0.766倍
4、当D不变，使M加大m倍时，异常也同样加大
[( x + a ) 2 + H 2 ][( x − a ) 2 + h 2 ] V xz = Gσ ln [( x + a ) 2 + h 2 ][( x − a ) 2 + H 2 ] H h H h ) − tg −1 − tg −1 + tg −1 V zz = 2Gσ (tg −1 x+a x+a x−a x−a ⎡ ⎤ x+a x+a x−a x−a − + − V zzz = 2Gσ ⎢ ⎥ 2 2 ( x + a) 2 + H 2 ( x − a) 2 + h 2 ( x − a) 2 + H 2 ⎦ ⎣ ( x + a) + h
第四节 地质体参数的计算 重力勘探5-正反演
i ) ln
2 i 1 i2
2 i 1
2 i
(i1
i
)
tg
1
i i
tg1
i1 i1
（二）任意形状三度体
1、线元法
➢用一组垂直于y轴的平面
和一组垂直于X轴的平面分 别切割地质体，则任意两 个平面的交线包合在地质 体之内的部分形成一个线 元。
x 时， g Gf h
1
P(x,0)
●x
h2
h 1 △σ △h
2
△σ △h
主剖面异常曲线单调变化，断层正上方梯度最大；平面异常等值 线呈条带状分布，与断层线平行。
在前述三个特征点上，异常值与埋深无关； 异常形态与埋深有关，埋藏越浅，水平梯度越大。
等值线为一系列平行台阶走向的直线，在断面附近等值线最密， 称为“重力梯级带”，且异常向台阶延伸方向单调增大。
第四节 地质体参数的计算
正演与反演
正问题也称正演，是指给定地质异常体的形状、产状 和剩余密度分布，通过计算得出重力异常的大小、特 征和变化规律等。
反问题也称反演，是指根据重力异常的数值大小、变 化规律等场的特征，结合已知的地质资料和地质体的 物性参数，求解地质体的形状和空间位置等。
正问题从给定地球物理模型，通过数值计算或物理模拟，得 出相应地球物理场的过程，目的是认识和掌握地球物理场的 特征与场源之间的对应关系；
当α=90°（垂直断层）时，重力异常极大值 与极小值绝对值相等，曲线以原点O为中心对 称
当α<90°（正断层）时，下降盘一侧异常极 小值明显
当α>90°（逆断层）时，上升盘一侧异常极 大值明显
16重力勘探-重力异常正反演解析
△gz
△g
FHale Waihona Puke rh1R0
测量垂直梯度原理 gz h2 h1
g h2 g h1
g ( z z ) g ( z )) g h1 g h2 z h2 h1
△g △g(x+△x) △g(x-△x)
△gx
A（x,0,0）
△g
F r
h
1
R
0
2 x1
2 3
2 3
2 3
h2
n 2 3
（n 1 ）
x1 h
n
（n 1）
2 3
• （3）反演剩余质量
m g max G 2 h
• （4）反演半径
g max h 2 m G
g max h 2 m G
4 3 m R 3
3 m R( ) 4
• 重力异常的正反演（正反演问题的关系：异常场源（地
形状
质因素产生的剩余质量）和重力异常之间的对应关系（互相关系）包 括数量上关系。 ）
大小
异常场源 位置 产状 深度 物性
根据数学物理方法： 万有引力 重力异常的 推断：定性或者定量 △g
分布规律 形态特征 幅度大小
A F △g
1）正问题是反问题的基础； 2）反问题强烈依赖于正问题。
2 7 2 2
②均匀的水平圆柱体（二维水平柱状体）
• • • • • •
在实际的地质现象中，如长轴背斜、向斜等， 可以近似看成水平圆柱体来讨论。 水平柱状体：向两端无限延伸 半径：R 埋深：h 延伸方向：y 剩余密度
定义：线密度
S R 2 dd
hx g x 4G 2 ( x h2 )2 h2 x2 g z 2G 2 ( x h2 )2 h 2 3x 2 g zz 4Gh 2 ( x h 2 )3
三维变密度体重力异常的正演
三维变密度体重力异常的正演
本报告旨在详细分析三维变密度体重力异常正演的结果。
首先，研究者使用三维变密度体模型对研究区进行建模，并以文献中提及的标准精度和公式进行重力正演计算，得出体重力场异常数据。
其次，分析者对正演的体重力场异常结果进行加工处理，构建具有原型性的重力场异常模型，同时运用高斯拟合方法优化模型，得出研究变量之间的关系及微地重力场异常等物理成分。
再次，分析者结合总气压系统，通过精细观测获得三维体重力异常数据，并综合各不同方位观测面上的实测异常值，进行拟合分析得出：研究区域的重力场异常显示出极为复杂的分布特征。
最后，回归分析结果表明：研究性变密体的重力异常总体变化趋势表明，在研究区域重力异常值的变化整体上表现为先增后减的趋势，总体变化符合模型得出的变化趋势，说明三维变密度体重力异常正演是一个有效的研究方法。
综上所述，本研究对三维变密度体重力异常正演进行了详细的分析，并结合高斯拟合方法，进行了定量\定性分析，验证了
三维变密度体重力异常正演计算的有效性和可行性。
《重力正演计算》PPT课件
g G [( x ) ln ( x )2 H 2 (x )2 h2
(x
)
ln
(x )2 H 2 (x )2 h2
2 H (arctan x arctan x )
H
H
2h(arctan x arctan x )]
h
h
第二节 规则形状异常密度体的正演问题 第三章 重力异常的正演问题
P(x,0)
x ●
△g(0) > πG△σ△h
➢ 当台阶面向台阶内侧倾斜时：
△h
△g(0) < πG△σ△h
第二节 规则形状异常密度体的正演问题 第三章 重力异常的正演问题
3.2.4 倾斜台阶
倾斜台阶基本特征
（2）Vxz:在 90时，曲线才是以纵轴为对 称轴，极大值位置正对应断面位置;而在
2πG△σ△9时0h，曲线为不对称形状，极大值偏向 斜面所在一方。
V zz
G
{1 2
sin
2
ln
(H (h
cot cot
x)2 x)2
H2 h2
H x cos
h x cos
2 sin 2 [arctan sin
arctan sin
]}
x sin
x sin
V zzz
2G
sin 2
[ (h
x 2h cot cot x)2
h2
(H
x 2H cot
✓ 当σ>0时，极大值一侧对应着上升盘，极小值一侧对应着 下降盘，在极小值十分清晰且大于极大值的绝对值时，属正断 层类型，反之则属逆断层类型。
第二节 规则形状异常密度体的正演问题 第三章 重力异常的正演问题
3.2.5 二度长方体（铅垂柱体）
16重力勘探-重力异常正反演
] 3Gm
mh
(x h ) h
2
3 2 2
] G m 2h 2 x 2 (x h )
2 5 2 2
g zz
g x
g z
[Gm g z g z g zz z h
2h 2 x 2 ( x2 h ) h
5 2 2
] 3Gmh 2h 2 3 x 2 (x h )
x1
h1 h2 h 2 )
h x 1
4
1 1 4 2tg 2 h x1 1 1 4 2tg 2 h x1 1 1 4 tg 4 h x1 1 1 1 4 tg tg ( 1) 4 h h x1
n
2G n 1
（3）反演圆柱体半径
R
R
2
• 如果剩余密度已知，圆柱体半径：
（4）导数异常
g z max h2 x2 2G 2 (x h 2 ) 2
x 0
h 2 2G 2G 4 2 h h
△gz
h2 x2 g z 0 2G 2 h x0 2 2 (x h )
g ( x x) g ( x x)) 测量水平梯度原理 g x 2x
△g • 反演： • （1）等值线的球心（x0,y0,0）→异常体中心 在地表的投影 A（x,0,0） • （2 ）反演深度
1 g 1 g max n n
△g
F r
h
g 1 G
n n
m h
h h2 h1
h h h h1 h2
2
O P(x,0,0) 1 2 Nhomakorabeal 2 l1 2 1
重力异常正演
倾斜台阶
倾斜台阶
倾斜台阶
两种延伸方向台阶异常的对比图
垂直台阶
垂直台阶
垂直台阶
垂直台阶
垂直台阶
垂直台阶
垂直台阶X方向水平梯度
垂直台阶高阶导数
垂直台阶
垂直台阶水平和垂直方向一阶导数
垂直台阶
垂直台阶垂直方向重力二阶导数
断层重力异常
断层重力异常
铅垂柱体
铅垂柱体
铅垂柱体
铅垂柱体
密度均匀的球体
g V z G v[( x )2 (( z ) y d )2 d (d z )2 ]3 /2
密度均匀的球体
密度均匀的球体
gGM[x2y2 h0 h02]3/2
密度均匀的球体
g
GMh0 (x2 h02)3/2
球体重力异常图
球体重力异常图
利用已知异常计算球体参数
密度均匀的水平圆柱体
VXZ4G [( (x )2x )((hh00 zz))2]2
VZZ2G[((h0x)z2)2(h(0zx))22]2
VZZZ4G 3([( x )2x()h20 (zh)0 (zh)02 ]3z)3
密度均匀的水平圆柱体
• 若以水平圆柱体的轴 线作为Y轴，Z轴垂直 向下，在轴线上取一
利用已知异常计算球体参数
三度球体引力位高阶导数
球体引力位高阶导数
球体引力位高阶导数（主剖面）
一阶水平和垂直导数
二阶垂直导数
两个球体组合模型理论重力异常
两个球体组合模型引力高阶导数
密度均匀的水平圆柱体
对于某些横截面近于圆形、沿水平方向延 伸较长的地质体，如扁豆状矿体、两翼较 陡的长轴背斜及向斜构造等，研究它们的 异常时，在一定精度要求内，可以当成水 平圆柱体的异常来对待。 对于无限长水平圆柱体所引起的异常，完 全可以当作质量集中在轴线上的物质线看 待。
重力异常正演计算知识讲解
重力异常
• 1、重力异常的概念
• 地下物质密度分布不均匀引起重力随空间位置的 变化。
• 在重力勘探中，将由于岩石、矿物分布不均匀所 引起的重力变化，或地质体与围岩密度的差异引 起的重力变化，称为重力异常。
重力异常
重力异常
• 在观测的重力值中，包含了重力正常值及重力异 常值两部分。
• 用实测重力值减去该点的正常值，也能够得到重 力异常。
VXZ4G S [( (x )2x )(( zz))2]2dd
密度均匀的水平圆柱体
• 对于剩余密度均匀的无限长水平圆柱体，可视为 质量集中在轴线上的物质线
d d S （S 是水平圆柱体的横截面积 ）
S
g2G(x)h 20(h z0z)2
密度均匀的水平圆柱体
VXZ4G[( (x )2x)((hh00zz))2]2
VZZ2G[((h0x)z2)2(h(0zx))22]2
VZZZ4G 3([( x )2x()h20 (zh)0 (zh)02 ]3z)3
密度均匀的水平圆柱体
• 若以水平圆柱体的轴 线作为Y轴，Z轴垂直 向下，在轴线上取一
单位长度， dmd
• 若水平圆柱体有限长 ，则
重力场的等价性
• 重力场的等价性： • 地下不同深度、形状、密度的地质体
在地表面可引起同样的重力异常。 • 重力场的等价性给重力异常的解释带来一
定的困难。
简单规则几何形体的异常
• 为了简化，假设地质形体孤立存在，密度均匀， 地面水平，所取剖面为中心剖面。
• 规则形体：球体、水平圆柱体、垂直台阶、脉状 体……
V
(z)d d d
g z V zG v[(x)2 (y)2 (z)2]32
g V z V zG v[(x) 2( ( z)y d )2d (d z)2]32
第八讲 重力异常反演课件
应用重力学第八讲重力异常反演d?解正问题是解反问题的基础，解反问题是目的。
仅从地质角度，解重力反演问题的目标9矿体类问题:寻找、研究或推断金属或非金属矿体；9构造类问题:研究地质构造，包括控矿构造，如含石油、天然气、煤的构造以及区域性的深部构造等。
从地球物理角度，解重力反演问题的目标9矿体类问题:确定地质体的几何和物性参数；9构造类问题:确定物性分界面的深度及起伏；9密度分布问题:确定密度的分布。
一、计算地质模型体的几何及物性参数(一)直接法直接利用由反演目标引起的局部异常，通过某种积分运算和函数关系，求得与异常分布有关地质体的某些参量。
(二)特征点法根据异常曲线上的一些点或特征点（如极大值点、零值点、拐点）的异常值及相应的坐标求取场源体的几何或物性参数;仅适用于剩余密度为常数的几何形体。
异常曲线形态分类第一类是单峰异常，零值点在无穷远处如球体的Δg曲线、台阶的Vxz曲线等；第二类是具有极大值、极小值和一个零值点如球体的Vxz曲线、台阶的Vzz、Vzzz曲线；第三类是具有一个极大值、两个极小值和两个零值点如球体、水平圆柱体的Vzz和Vzzz曲线；第四类是台阶的Δg曲线，一边高一边低的形态应用条件对异常作平滑处理，尽量准确确定原点的位置； 对异常曲线作分离处理，获得单纯由研究对象引起的异常；对剩余（局部）异常进行分类，判明该异常的场源体接近于何种可能的几何形体，然后选用相应的反演公式。
2223/2212()(GMDGMDg x D x D Δ=++)(6524.2/12/1x x ′−)(4811.03/13/1x x ′−)(4056.04/14/1x x ′−{{{ D1/2)nGD πμ=(三)选择法根据异常分布和变化特征，结合地质和其他地球物理和物性等资料，给出初始地质体模型;进行正演计算，将理论异常与实测异常对比; 若两者偏差较大，对模型进行修改，重算其理论异常计算，再次进行对比……;如此反复进行，直至两种异常的偏差达到事前要求的误差范围为止，则这最后的理论模型就可作为所求的解答了。
重力异常正演总结
一、 正方体重力异常正演推导过程（蒋）直立长方体外任意一点的重力异常公式推导如下：p(x,y,z)： 长方体外任一点ρ： 长方体密度θ(ε,η,τ): 长方体内任一点坐标222111222321112221()(,,)[()()()]a b h a b h a b h a b h GM V G d d r V z g x y z G d d d Z x y z ρεηττρεητεητ==∂-∴∆==∂-+-+-⎰⎰⎰⎰⎰⎰设：r =；ln[()]u y r η=-+；则1.()xdu d y r r εεη-=-+要求原式：222111ln[()h b a a b hG y d ρηε--+⎰在此先求：ln[()y d ηε-+⎰22.()11.()()...()()().ln[()].[()].u d x u x x d y r rx x y r d y r r εεεεηεεηεη=-=----+-=--+--+⎰⎰⎰上式第二项：2().[()].x d y r r εεη--+⎰222222222().[()].[()()].()()().()()...()()[()()].x r y d x z rz y y z d x d d d x z r x z rεηεεττηητεεεεετετ---=-+-----=---+-+--+-⎰⎰⎰⎰⎰上式最后一项：222().().[()()].y z d x z r ητεετ---+-⎰2222222222222222222().().()2().()222().().[()()]..[()()].[()()]().().[()]..[().().()]().()..()().().y z x r x y ry z y z d r x z y z y z r x d r x y r z y z d r z x y ητεεηητητεετητητεεεητητετεη--------+-=-+--+-----=--+----=+---=-⎰⎰⎰212222().[.()]()..[().().()]().()rz x r z d x y r z y x τετεεητηε------+---⎰212().[.()]()..().()().()rz x r z r d d y x y x τετεηεηε⎡⎤----=⎢⎥----⎣⎦⎰2222.()().()1().().().()1().().arctan ().()r z x y z r z d y x z r z y x τεηττηεττηε---⎡⎤-=--⎢⎥--+⎣⎦⎡⎤-=--⎢⎥--⎣⎦⎰由上述推导过程，可以得出：ln[()y d ηε-+⎰().ln[()]().arctanx x y r z z εεηεττ-=--+-+--().().ln[()]().arctan ().()z r y x r z y x τηετηε⎡⎤-+--++-⎢⎥--⎣⎦ 222111(,,)ln[()]h b a a b hg x y z G y d ρηε∴∆=--+⎰222111().().ln[()]().arctan ().ln[()]().arctan ().()h b a ab hxz r G x y r z y x r z z y x ετρεηετηεττηε⎡⎤--=---+-+-+--++-⎢⎥---⎣⎦将上式完全展开后，().arctan xz z ετετ----可以消除。
第八讲 重力异常反演
重力异常反演d仅从地质角度，解重力反演问题的目标9矿体类问题:寻找、研究或推断金属或非金属矿体；9构造类问题:研究地质构造，包括控矿构造，如含石油、天然气、煤的构造以及区域性的深部构造等。
从地球物理角度，解重力反演问题的目标9矿体类问题:确定地质体的几何和物性参数；9构造类问题:确定物性分界面的深度及起伏；9密度分布问题:确定密度的分布。
一、计算地质模型体的几何及物性参数(一)直接法直接利用由反演目标引起的局部异常，通过某种积分运算和函数关系，求得与异常分布有关地质体的某些参量。
(二)特征点法根据异常曲线上的一些点或特征点（如极大值点、零值点、拐点）的异常值及相应的坐标求取场源体的几何或物性参数;仅适用于剩余密度为常数的几何形体。
异常曲线形态分类第一类是单峰异常，零值点在无穷远处如球体的Δg曲线、台阶的Vxz曲线等；第二类是具有极大值、极小值和一个零值点如球体的Vxz曲线、台阶的Vzz、Vzzz曲线；第三类是具有一个极大值、两个极小值和两个零值点如球体、水平圆柱体的Vzz和Vzzz曲线；第四类是台阶的Δg曲线，一边高一边低的形态应用条件对异常作平滑处理，尽量准确确定原点的位置； 对异常曲线作分离处理，获得单纯由研究对象引起的异常；对剩余（局部）异常进行分类，判明该异常的场源体接近于何种可能的几何形体，然后选用相应的反演公式。
2223/2212()(GMD GMDg x D x D Δ=++)(6524.2/12/1x x ′−)(4811.03/13/1x x ′−)(4056.04/14/1x x ′−{{{ D1/2)nGD πμ=(三)选择法根据异常分布和变化特征，结合地质和其他地球物理和物性等资料，给出初始地质体模型;进行正演计算，将理论异常与实测异常对比;若两者偏差较大，对模型进行修改，重算其理论异常计算，再次进行对比……;如此反复进行，直至两种异常的偏差达到事前要求的误差范围为止，则这最后的理论模型就可作为所求的解答了。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
设：
（D.1）
所以：
按式（D.1）的形式累加起来，最后只需要求一次反正切函 数，这样处理后，计算速度提高一倍以上。
1.2.2 直立“线元”法
某工区物探、勘探工作布置示意图
正演问题的定义： 根据巳知的、具有剩余质量的地质体的形状、产状和剩余密度 分布，通过理沦计算，研究它们所引起的异常及其各阶导数异 常的数值大小、空间分布和变化规律。 反演问题的定义： （1）由观测上重力异常的分布，在给定物体边界位置函数的条 件下，求解物体的密度分布函数；（物性反演） （2）由观测面上重力异常分布，在给定物体密度函数的条件下， 求解物体的边界位置的数值；（几何反演） （3）由观测面上重力异常分布。在给定特殊约束（如设物体密 度均匀、形态规则)条件下，求解物体密度参数和几何参数。 给定的函数和特殊约束称为反演问题的定解条件。
什么是正问题与反问题？
反问题：m=G-1d
观测数据d
地质模型 m
正问题：d=Gm
（一）规则形体的正、反演问题
为了简化，假设地质形体孤立存在，密度均匀，地 面水平，所取剖面为中心剖面
规则形体：球体、水平圆柱体、垂直台阶、脉状体……
1、球体
规 则 形 体 的 正 、 反 演 问 题
近似于等轴状地质体，如盐丘、矿巢、溶洞等
lim g 0; g max
g max h02 m G 2 ; x1/ 2 0.766h0 ; m h0 G
2、水平圆柱体 2、水平圆柱体（线质量）
规 则 形 体 的 正 、 反 演 问 题
小柱体元在P(x,0,0)点产生的重力异常为
g G
h0 dy
(x y h )
用一组垂直于y轴的平面
和一组垂直于X轴的平面分 别切割地质体，则任意两 个平面的交线包合在地质 体之内的部分形成一个线 元。 用解析式计算每一个线 元在计算点产生的重力异 常作用值。 对所有钱元的作用值依 次进行X方向和Y方向的数 值积分，便得到整个地质 体在计算点所产生的重力 异常值。
1、线元法
规 则 形 体 的 正 、 反 演 问 题
断层或不同岩层的接触带都可作为台阶处理
x 2 h12 1 x 1 x g G (h x ln 2 2h1tg 2h2tg ) 2 x h2 h1 h2
由公式可见：
当X→∞时， △g =2πG△σ△h 当X→ －∞时， △g =0
球体参数：半径30m，中心埋深20m，剩余密度 2.0g/cm3，重力异常单位mGal
1、球体
规 则 形 体 的 正 、 反 演 问 题
4 3.5 3 2.5 2 1.5 1 0.5
x 10
-3
P(x,0,0)
-40 -30 -20 -10
0 -50
h0
0
10
20
30
40
50
σ0
ΔE
Δg
σ
讨论：
x
重力异常的正反演
• 1. 重力异常的正演问题、反演问题； • 2. 均匀密度球体、水平圆柱体、台阶的重 力异常正演方法，异常特征，反演方法； • 3. 密度界面的剩余密度的确定方法； • 4. 单一密度界面异常的特征及反演解释方 法（近似解法、矩阵法）； • 5. 解复杂密度体正演问题的基本思想； • 6. 最优化选择法的基本思想；
用解析公式计算出每个小长方
最后，将所有长方体的重力异
常值累加，以求得整个地质体在 计算点的异常值。
体在计算点所产生的重力异常值。
点元法 “点元”法所取的各个点元的体积可以相同，也可不同。各 个点元的物性可以相同，也可不同。通常是将勘探剖面之间的 地质体用适当的长方体或立方体来近似，确定出各个点元的角 点坐标，即可计算出该点元的三重积分值。 对于一个点元而言，其计算公式如下：
2 2 2 3/ 2 0
整个水平圆柱体在P点产生的重力异常为无 穷多个柱体在该点产生的重力异常之和，即：
g G h0
2 h0 dy G 2 2 2 2 3/ 2 2 ( x y h0 ) x h0
2、水平圆柱体
规 则 形 体 的 正 、 反 演 问 题
讨论：
3、垂直台阶
2．面元法
用一组垂直于z轴的平面
或者垂直于X轴、y轴的平 面切割地质体，地质体与平 面相交形成一系列的裁面。
近似地用一个多边形代
替每一个截面，用解析表达 式计算出计算点的重力作用 值。
将所有面的作用值用数
值积分求得整个地质体产生 的重力异常值。
3．长方体元法
用一组垂直于X轴的平面、一
组垂直于y轴的平面和一组垂直于 Z轴的平面切割地质体，于是地 质体被划分成许多小长方体。
式中：
对所有点元求和：
Fi f ij
k 1
L
i 1,2, M
1、对计算ln项的简化 从式(1.2.1)式的V0，V2，V3．V5的表达式可以看出，需 要计算一系列带ln的项,然后求和，由于大量调用ln标准子 程序会花费很多计算时间并影响计算精度，所以有必要作 适当的简化。 利用： 对L个点元有：
当X=0 时， △g = πG△σ△h
2πG△σ△h πG△σ△h
P(x,0)
o 1 2
h1
●
x
h2
△σ
△h
△σ
△h
3、垂直台阶
规 则 形 体 的 正 、 反 演 问 题
平面异常特征：
等值线为一系列平行台阶走向的直线，在断面附 近等值线最密，称为“重力梯级带”，且异常向 台阶延伸方向单调增大。
当台阶倾斜时：
一、重力正、反问题的解法
任意形体重力异常计算公式：
g G
V
( z ) d d d ( x) ( y ) ( z )
2 2 2 3/ 2
式中：G为万有引力常数，(x, y, z)为观测点P的坐 标，V为地质体体积，Δσ为剩余密度， (ξ,η,ζ)为剩 余质量元 m d d d 的坐标
2πG△σ△h
由图可见：无论台阶产 状如何，异常的形态相 似，仅原点处的异常值 不同。 当台阶直立时：
P(x,0)
πG△σ△h
o
△σ
●
x
△g(0)= πG△σ△h
△h
当台阶面向台阶外侧倾 斜时：
△g(0) > πG△σ△h
当台阶面向台阶内侧倾 斜时：
△g(0) < πG△σ△h
二、不规则三度体的正演问题