重磁数据处理与解释

重、磁异常数据处理与解释部分
重、磁异常解释： 定性解释 定量解释— 重磁异常反演
问题： 1. 观测重、磁异常为叠加异常 2. 重、磁异常反演的不唯一性
重、磁异常反问题的多解性
例子：重力异常
解决途径：
采用综合分析方法，利用各种地球物理 资料与地质、钻井资料，以增加约束条件， 减少解释的多解性。
将复杂异常进行简化 — 重、磁异常的划分
2. 重、磁异常划分的任务： 根据不同的研究目的，需从叠加的异常中
提取出我们的研究对象，作为有用信息。
3. 重、磁异常划分方法：
分解法： 如 异常曲线平滑法、平均场法、 小波变换方法
场变换法：如 重磁异常导数变换、解析延拓 剥层法：
重、 磁异常资料的处理常规方法
重、磁异常数据
向上延拓：
由地表值，换算到 空中或地下某一深度的 重、磁场值
1）计算原理（以重力为例）
重、磁场满足拉普拉斯方程
2W x2

2W y 2

2W z 2
0
对 z 求偏导数，得
3W x2z

3W y 2z

3W z 3
0
2W 0
利用 g W z
2 g x2

2 g y 2

2 g z 2
（1）突出浅而小异常体的异常特征，压制区域 性深部物质引起的异常特征
一定程度上划分了不同深度和大小的异常源 产生的叠加异常
（2）导数阶次越高，对浅部异常体反映越敏锐
不同阶次的重力导数对不同埋深地质体的 反映不同
（3） 提高了对异常的分辨能力
2）水平导数的作用 突出线性异常带
2.垂向二阶导数计算
重、磁异常是叠加异常，来源于地下不同的 物质源，解释中希望将不同场源的异常分开
2. 重、磁异常数据处理的目的
将各种场源引起的异常分开，用于定量反 演计算与定性解释
3. 数据处理的思路
根据重磁异常特点
异常体埋深、规模大，异常宽缓，异常 值幅度大，在频率域中表现为低频成分多
一般异常体规模、埋深小，异常宽度窄， 幅值变化大，在频率域中表现为高频成分多
18 17 16 15 14 13 12 11 10 9 8 7 6 5 4
东海长江坳陷及邻区 布格重力异常图
东海长江坳陷及邻区重力异常 向上延拓 25km 图
2.4 重、磁异常数字滤波原理及其特性
2.4.1 基本原理
设空间域重（磁）异常值 Δg(x)或 Δg(x,y) 的傅氏变换为

g(u) g(x)ei2uxdx
第1章 重、磁异常解释的标志
1.1 断裂带的重磁异常标志：
1.1.1 重力异常主要标志： 1. 梯度带： 表现为某一方向递增（减）的等值线密集 带，反映了密度界面深度上的突然变化。
黄骅坳陷及邻区布格重力异常与沧东断裂
ES
ES
W92-96 ES
K85-120 ES
S87-229
K85-70 ES
基底断裂
串珠状磁异常：
塔里木盆地东部 Za（化极）垂向二阶导数
库车
轮台
库尔勒
nT/k
草湖
孔 尉犁 群克尔
辛格尔 兴地 库 鲁 克 塔 格
陷车泉
20
学堂
雀
团结村向阳村
铁干里克
河
满参1
阿拉干
罗
布
泊
10
红柳井
河 罗布庄
米兰
0
臣
若羌
塔中1
尔
阿尔金山
车
-1
且末
休堂 安迪尔
-2
0
100km
ddz-2
Z塔a（里化木盆极地）东垂南向边一缘阶导数
按照异常特征划分：局部异常、区域异常
重、磁异常分布一般规律
地质异常体
重力异常
有明显走向 平面等值线沿一定方向延伸，长 （二度体） 轴方向与地质走向一致。
磁异常 有相似的特征。
无明显走向 平面等值线也无明显走向， （三度体） 称等轴状异常。
物性、规模 剩余密度大，异常值大。 异常体规模大，异常规模大。
h
2


g ( ( 2
,,0)
2 h2
dd
)
z =-h
取柱坐标，作坐标变换：
(x,y,-h)
2 2 2 dd dd
o
ρ
x
换元得：
(ξ,η,0)
y
z
g(0,0,h)
1
2
2 0

g
0
(
,
,0)
(
2
h
h2
空间域划分
频率域划分
区域异常 局部异常 低频异常 高频异常
解释深部的、 大范围的地质体
解释浅部的、 小范围的地质体
2.2 重、磁异常划分的一般概念 1. 局部异常与区域异常：
区域异常： 指分布范围大，幅度大， 变化幅度比较平缓的异常
局部异常：指分布范围小，幅度小， 变化梯度较大的异常
例如：
g g区＋g局

g
(r
)

g
(o)

1 n
n

i1
gi
(3) 最佳半径选择
r R(最佳半径) g(r) g区
3. 面积平均法
实例：
实例：
2.3.2 重、磁异常的导数变换
1. 导数变换的作用
导数变换可包括： 1）垂向导数（Vzz、Vzzz） 2）水平导数（水平梯度）（ Vxz）
1）垂向导数的作用
K(ρi,2) 0.04034 0.12988 0.07588 0.14559 0.07651 0.09902
K(ρi,3) 0.01961 0.06592 0.05261 0.10563 0.07146 0.10226
K(ρi,4) 0.01141 0.03908 0.03566 0.07449 0.05841 0.09173

g ( x, x
y)
i

g ( x, y
y)
j
n

(i ,
j)
水平梯度模数
g n


g
(
x

x, y) x
-
g
(
x,
y
)

2


g
(
x,
y

y) y

g
(
x,
y)

2
1

2
或
g g(x x, y) - g(x, y)

0
2 g z 2


2 g x 2

2 g y 2

2）垂向二阶导数计算举例
近似公式：
2 g z2

1 r2
[ A1
g(0)
A2
g(r) A3
g(
2r) A4 g(
5r)]
其中： g(nr) 为nr半径圆周上重力值之和
gzz (0)
实例：
鄂尔多斯及邻区 航磁ΔT 与上延30km图
实例：
闽 浙
隆
美人峰1井
虎皮礁隆起
起 长江坳陷
带
海礁隆起
西湖凹陷
10 g.u.
28 26 24 22 20 18 16 14 12 10 8 6 4 2 0 -2 -4 -6 -8 -10
闽 浙
隆
美人峰1井
虎皮礁隆起
起 长江坳陷
带
海礁隆起
西湖凹陷
10 g.u.
1.2 岩体的重、磁场特征
1. 火成岩侵入体的重磁异常特征： 异常形态规则、幅度大，明显存在磁源重力
异常
2. 火山岩的重磁场特征： 异常零乱，幅度变化大，垂向导数异常明显，磁
源重力异常不大
利用航磁异常 圈定了花岗闪长岩 侵入体，并区别了 周围的熔岩。
3. 岩体的重力异常特征 — 异常取决于岩体与围岩的密度
g(x, y,0) 数学变换g（x, y, h)
Z=-h x
Z=h
2. 解析延拓的作用
(1) 向上延拓突出埋藏较深、体积较大的场源 引起的异常（突出区域异常）。
(2) 向下延拓使浅部异常的影响相对加强 （突出局部异常）。
向上延拓抑制浅部磁性体的异常
3. 向上延拓计算举例：
g(0,0,h)
沿休堂、且末至罗布庄也存在明显正磁异常带， 与重力对应很好，孔有雁列分布特征（车尔臣断裂）
库鲁克塔格
罗布庄以东有转为NEE向雀 正磁异常带
河
盆地边缘还存在几乎平行的另一条高磁异罗布 常带， 向南有两条NEE向高磁异常带在且末的南部与泊之相交
盆地东北边缘Za（化极）垂向二阶导数
17
有四条高磁异常带
从尉犁向SEE延 伸经兴地南转向近EW －NEE向；
由于断裂切割深， 岩浆多次侵入或喷 出，形成主要由火 成岩引起的剧烈变 化的强磁异常带。
基底断裂
3. 磁异常带轴向的水平错动。
实例： 由于断裂使磁
异常的轴线发生 水平错动
4. 挤压破碎带可表现为狭长的负异常带
小结： 断裂带的重、磁异常主要标志
（1） 梯度带 （2） 狭长异常带，或线性异常带 （3） 异常轴的平面位移 （4） 等值线有规律的扭曲 （5） 两侧异常明显不同的分界线 （6） 串珠状磁异常
有相似的特征。
磁化强度大，磁异常值大 有相似的特征。
埋藏深度大 异常值变小，异常曲线宽度变大， 有相似的特征。 梯度变小。
球体：
Δσ>0，Δg正异常（无负值） Δσ<0， Δg负异常（无正值）