矿产勘探地球物理技术：从原理到应用

山东地矿局总工研讨班—学术讲座之一

矿产勘探地球物理技术：从原理到应用

山东地矿局总工研讨班—学术讲座之一

汇报提纲

1引言—两个原理模型

2勘探技术—仪器设备

3勘探技术—新方法

4勘探趋势：从矿床到成矿系统

/2/

一、引言-两个模型

勘探地球物理：原理模型

LiYaoguo，2013

一、引言-两个模型

磁法

直流电法

电磁法

反射地震

天然场电法法

LiYaoguo，2013

一、引言-两个模型

原理模型的数学描述：

勘探地球物理：原理模型

数据处理：把一次场（源场）、地形、地表、路径仪器响

应、干扰等影响去掉，保留探测目标的响应（二次场）；

反演及解释：恢复目标体的空间形态和物质属性，给出合

理的地质解释

LiYaoguo，2013

一、引言-两个模型

勘探地球物理：工作模型

（工作流程）

几个值得注意的关键环节（贯

穿整个模型）：

u施工设计：采样间隔与网

格化

u数据处理（异常改正、分

离、去噪等）：经验与判

断

u反演解释：2D与3D反演与

地质解释

一、引言-两个模型

1、施工设计Nyquist采样定理：连接连续（时空）信号和离散信号的基本桥梁

采样定律告诉我们什么？

1、目标导向：明确或估计地

球物理调查目标的空间尺寸、

形态；

2、合理设计地球物理测量的

采样频率大于2倍的信号最高采样频率小于信号最高频率，不

点、线距，2D或3D频率，可以完全恢复信号能恢复信号，产生假异常

3、点、线距离都要小于调查

目标空间尺寸的2倍以上

Fc截至频率=1/2∆

Dentith，2013

一、引言-两个模型

两个实例：白线

采样

点、线距离设计对能否发现异常和线距

确定异常的真实形态致关重要

幅值、形态有很大影响，直接决定解不同点

线距对异常释结果（右图）

Dentith，2013

一、引言-两个模型

网格化间距大小的选择：

太小：引起串珠状或豆荚状异常

太大：引起低频假异常

原则：数据平均点距的1/2

测点太稀易产生“牛眼”异常

Dentith，2013

一、引言-两个模型

点、线距不匹配造成的串珠、点线距不匹配造

香肠状异常、台阶、梯状异常成假异常A和B

当线距与点距相

差较多时（1:10

或1:20），网格

化方法的选择极

为重要，一般以

1/5或1/3线距

为网格化间距

Dentith，2013

一、引言-两个模型

2、位场分离（区域-局部场分离）

u位场是一个三维体积效应，由上地幔和地壳所有物性差异产生场的总和

u异常通常淹没在由下列场源形成的“区域异常”中：

(i)更大的场源；(ii)更深的场源；(iii)或更远的场源.

LiYaoguo，2012

一、引言-两个模型

两个典型例子

Dentith，2013

一、引言-两个模型

位场分离方法：

u最小平方拟合:使用简单

多项式（二阶、三阶趋势）

拟合区域场

u滤波法：把长波场分量

当做区域场；在空间域和

频率域都有效

u“剥皮（Stripping）

法”：通过对区域场源的

正演模拟，去掉区域场

Dentith，2013

一、引言-两个模型

区域-局部异常分离

总场

理论模型

一、引言-两个模型

（1）最小平方拟合:

LS估计的区域异常理论区域异常

剩余异常

效果：

•基本近似

•剩余异常畸变

LiYaoguo，2012

一、引言-两个模型

矿体重力异常分离对认识矿体分布空间范围起到关键作用

u地下700米测量精度要求：磁法测量精度

泥河铁矿原始布格异常位场分离后的矿体异常

重力异常要求10-20nT;重力测量需要30-40微伽；

u地下1500米测量精度要求：磁法测量精度

要求5-10nT;重力测量要求10-20微伽

上图为安徽泥河铁矿重磁场及3D反演结果。精地面高精度磁异常

确定位矿体的深度和形态吕庆田等，2008

一、引言-两个模型

（2）频率域滤波法：

在频率域将能谱的高频部分与低频部分分开。基本假定：

长波长的能量代表区域异常，短波长能量代表局部异常，具体

计算流程如下：

低通滤波：提取区高通滤波：提取局

域场，通常与方向部场

导数一起组合

一、引言-两个模型

滤波估计的区域异常理论区域异常

效果：

•近似程度差剩余异常

•剩余异常畸变

一、引言-两个模型

滤波法的基本要求：

Ø数据必须是网格化数据

Ø平面数据观测假设，数据是平面数据

Ø在能谱上，区域异常和局部异常不能重叠过多

Ø数据范围要远大于拟提取的局部异常范围

滤波法的优缺点：

Ø使用FFT，快速易行

Ø快速处理大区域数据

Ø高度依赖于频谱特性

Ø经常使提取异常畸变

一、引言-两个模型

（3）基于反演的“剥皮”法（基于反演的位场分离）

B：扣除局部异常区

A：大尺度反演域，进行正演

D、反演解释

C=A-B：异常数据

AB

LiYaoguo，2012

一、引言-两个模型

局部场源

区域场源

LiYaoguo，2012

一、引言-两个模型

大范围反演真实模型

一、引言-两个模型

剥皮法结果理论值

一、引言-两个模型

3、地球物理反演：为什么？

这是一个测区的磁法数

为什么要开展反演解释？据，地下有几个磁性体？

u由于重磁、直流电场为静

态、弥散场，每个数据受地

下所有物质影响

u数据中的关于场源的信息

是混乱，不清楚的

u“模型域”的地质解释比

“数据域”容易很多

u了解深度信息

LiYaoguo，2012

一、引言-两个模型

1)数字模拟实验：单个良导体，观测数据加5%噪声

拟断面显示

LiYaoguo，2012

一、引言-两个模型

2)地表增加若干小的良导体干扰

一、引言-两个模型

实验数据:直流电阻率

拟断面大致反映地下良导体位置

-1000100120

2

N

-

s

p8

a

c

i

n

g45

干扰情况下，地下良导体位置在何处?

126

-1000100

2

N

-

s8

p

a4

c

i

n

g

一、引言-两个模型

2D反演结果:单一柱体

u单极-偶极;n=1,8;a=10m;

一、引言-两个模型

实例：地表加干扰的单一柱体反演结果

u单极-偶极;n=1,8;a=10m;

理论模型正演反演模型正演

一、引言-两个模型

反演的几点重要认识：

1、反演是一种重要地质解释手段；

2、反演不是数字游戏，随着地质约束的不断增加，多

数据联合反演解释技术的进步，反演逐渐演变为一种综

合地质解释方法；

3、反演的多解性将随着地质约束的增加、多数据的联

合反演，将逐渐消失。

一、引言-两个模型

地质约束反演:正确方向的跨越

江苏省地质调查院学术讲座之一

汇报提纲

1引言—两个原理模型

2勘探技术—仪器设备

3勘探技术—新方法

4勘探趋势：从矿床到成矿系统

/32/

二、勘探技术—仪器设备

Terraquest’sPocoHelicopter

航空地球物理技术新进展TDEMsysteminflight

u航空重力及重力梯度测量系统；GedexSystem公司宣称其高

分辨率重力梯度系统（HD-AGG）在空间分辨率为60m的情

况下，分辨率已经达到1EotvosHz-1/2。

u航空磁场总场测量、梯度及张量测量系统；

u直升机吊舱式和固定翼时间域航空电磁探测系统；

u半空间航空瞬变电磁测量系统；

u航空伽马能谱测量系统；

u无人机航磁、航放测量系统等

形成了测量参数多样、平台可选择、实时GPS定位的完整的航

空地球物理勘查技术体系。

二、勘探技术—仪器设备

地面及井中地球物理探测技术

u分布式三维激电探测系统。QuantecGeoscience公司ORION3D系统同

时采集直流DCIP和MT信息，像地震勘探一样实现快速、高效的3D覆盖；

u自动智能重磁测量系统。Scintrex公司的CG-6高精度重力仪分辨率达到1

微伽，同时实现自动测量、野外“蓝牙”传输和“遥控”测量；

u三维电磁探测系统。CroneGeophysics公司的Mark73DE-SCAN系统探

测深度达到4000m；

u高温超导传感器（HTS）和低温超导传感器（LTS）使TEM探测深度更

深、信噪比更高；

u超宽频带（BMT）MT感应式传感器使AMT和MT测量一次完成，大大

提高了探测效率；

u井中重力和全波形TEM技术已经成为常规勘探手段，极大提高了对“盲

矿体”的识别能力；

u全三维反演解释、多数据类型联合反演、岩性填图和综合建模技术极大

提高了深部资源探测的准确性，实现4000米甚至更深的探测能力。

二、勘探技术—仪器设备

1、航空重力及梯度梯度测量系统

u航空重力测量技术逐渐成熟，开始应

用于矿产、油气勘查，目前的精度更适合

区域地质调查

TAGS和地面重力比较

从重力仪或加速度计测量的比力

信息中扣除载体的运动加速度，

得到重力值。

nnnnnn

gfv(2ωieωen)v

二、勘探技术—仪器设备

Turnkey

Airborne

Gravity

System

(TAGS)

TAGSInstallationintoAires

NavajoTAGSSystem

搭载AIRGravSystem的

飞机AIRGravSystemGT-1A型航空重力仪

二、勘探技术—仪器设备

重力梯度测量系统:美国BellAerospace公司的Air-FTG（全张量梯度）系统，

澳大利亚BHP公司的Falcon航重梯度综合勘查系统，英国ARKeX公司的超导

重力梯度仪EGG等

2004-2008

BellGeospaceacquires

FTGdatainGOM1stARKeXcommencesbuild

commercial3D-FTGofnextgenerationsuper

systemavailableconductinggradiometer

the‘EGG’

Air-FTG3D重力梯度测

量系统

Geotech’sGT-2A

gravimeter.FUGRO

(FALCON)

二、勘探技术—仪器设备

航空重力梯度技术的总体现状：

Ø旋转加速度计重力梯度仪——商业应用，成熟技术

Ø超导重力梯度仪——研究热点，准实用技术

Ø原子干涉型重力梯度仪——探索阶段、试验样机，前沿研究

Gedex公司和马里力拓公司和西澳

ARKeX公司自行研制

兰大学联合研制大学联合研制

二、勘探技术—仪器设备

铅锌矿

油气田

煤矿

2mgal

铅锌矿

AdaptedfromFrankvanKann，2004

世界典型矿床（煤田和盐丘型油藏）重力异常与信号波长关系图；几个航空重力系统的噪声水平(i.e.,

Carson直升机系统1996,Carson1995,LCT1996和Sander2003固定翼系统)以及基于GPS的加速度

测量噪声水平

二、勘探技术—仪器设备

AdaptedfromFrankvanKann，2004

世界典型矿床的重力梯度与波长关系图.同图显示对应梯度测量系统的噪声水平

二、勘探技术—仪器设备

目前的主要系统对世界典型矿床的探测能力

AirborneGravimetry：

代表当前航空重力的探

测能力

FalconGravity：

代表当前5E精度旋

镍矿转加速度计航空重力

铜金氧化矿梯度系统的探测能力

铁矿

金刚石

煤矿VK-1Gravity：

基本金属代表下一代1E精度

超导航空重力梯度

系统的探测能力

AdaptedfromAirborneGravity2010

二、勘探技术—仪器设备

Kimberlite,Abnerregion,NorthernTerritory,Australia

找

矿从G数据中圈出了一个20E的异常，钻孔证

实DD

例实在地下10m发现了一个1.3公顷的主管状岩

二脉，旁侧发现了一个0.3公顷的次级岩脉。

：2005年打钻发现了一个金刚石矿。

寻

找

金

伯

利

岩

脉

二、勘探技术—仪器设备

实IOCG,ProminentHill,Australia

例

三

：

铁

氧

化

物

铜

金

型

矿上图是花了几年做的地面

床重力垂向一阶导数，下

勘

探图是FALCON™AGG系

统两周测量结果

2、航空电磁法技术：

Ø时间域航空电磁系统：向地下发射瞬变电磁场脉冲，接收地质体产生的二次场，研究地下

地质体的电性体分布特征，实现地质勘查目标。

Ø天然场航空电磁系统：利用天然场源，分别在空中接收垂向大地电磁信号、在地面接收水

平正交的大地电磁信号，计算倾子并反演解释，实现地质勘查目标。

陈斌等，2021

二、勘探技术—仪器设备

主动源电磁探测技术发展趋势：降低发射频率、增大发射磁矩、

采用直升机平台；高温SQUID磁传感器提高测量精度；全波场反演

提高反演可靠性。

AeroTEM

SkyTEMsystemVTEMsystemHeliGeotem

二、勘探技术—仪器设备

我国直升机航空电磁系统研制历程：

陈斌等，2021

二、勘探技术—仪器设备

与国外同类系统主要技术指标对比

硬支架系统直径发射电流发射磁矩动态噪声探测深度备注

CHTEM-I(863研发)12米490A250,000Am2±50nT/s280米“十一五”完成

CHTEM-II(863研发)12米550A280,000Am2±15nT/s320米“十二五”完成

CHTEM-III

12米624A352,860Am2±1.2nT/s400米“十三五”完成

(项目研发)

TM2Geotech公司软支架系统,因其磁矩、噪声

VTEM17.4米243A240,000Am±1nT/s350米和CHTEM-III相近因此一起比较

±0.243—Geotech公司软支架系统，指标为其在川

VTEMLITE17.6米329A320,000Am2400米

±0.96nT/s藏线实测数据指标。

软支架系统直径发射电流发射磁矩动态噪声探测深度备注

CHTEMs-I1,364,000

32米424A±0.4nT/s650米“十三五”完成

(项目研发)Am2

VTEMMax35米260A1,000,000Am2±0.8nT/s700米Geotech公司系统

Geotech公司

VTEMPlus26米295A625,000Am2±0.5nT/s650米

最常用系统

天然场系统直径动态噪声探测深度

CHMT-I(项目研发)8.5米优于1%1500米“十三五”完成

ZTEM7.4米优于1%2000米Geotech公司系统

陈斌等，2021

二、勘探技术—仪器设备

航空电磁法－－圈定咸水、半咸水、淡水、圈定地下卤水体

二、勘探技术—仪器设备

无人机航空物探（磁/放）综合测量系统

航磁：

CS-3铯光泵仪，动态噪声≤70pT，

转向差≤1.0nT，总精度优于

1.65nT

航放：

下视16升，256道，自动数字稳谱

仪，分辨率±0.22%，峰漂±0.67道

无

航测：人

航高≥120m地形起伏，偏航距机

≤15m，航程1600km/10h（单航

磁），测控范围200km（视距，地磁

面控制），800km（海事卫星）

二、勘探技术—仪器设备

无人机航空物探（磁/放）综合测量系统

航空放射性探测系统RS-601

无人直升机航磁系统多旋翼无人机航磁系统

二、勘探技术—仪器设备

1、数据质量优势

夜航测量时干扰小,飞控电脑精确

控制飞行高度和偏航距，不疲劳，

不受阳光等干扰，窗帘效应影响明（A）CH-3无人机(B)Y-12有人机

飞行高度统计

显由于传统航空物探

2、测量效率高

单架次8小时（约1600km），可夜

航，一飞行日可多架次执飞，可一

站多机同时测量，相同时间段内测

量效率是传统航空物探的1.5～2倍

（A）CH-3无人机(B)Y-12有人机

无人机和有人机航磁ΔT等值线阴影图对比（3000测线千米，1∶5万）

二、勘探技术—仪器设备

3、三维激发极化—3DDCIPIV

u传统同剖面二维DC/IP观测装置IV

–单极-单极（pole–pole）

VI

–单极-偶极（pole-dipole(PDR)

–单极-偶极（pole-dipole(PDL)IV

–偶极-偶极（dipole–dipole）I

中梯（温纳，）

–WennerV

缺点：

McDermott

u天然或人为噪声矿床IP反演

uAB与MN严格同剖面，空剖面

据LiYG等，

间分辨率差2012

u深度浅（n<6，受电缆

长度限制）

二、勘探技术—仪器设备

Newmont的NEWDAS3DIP采集系统

关键技术：

u无线采集与现代电子技术（现场

监控和分析、噪声压制）

uGPS时间同步

u3DDC/IP反演模拟技术（UBC-

GIF，DCIP3D）

NEWDAS系统特点：

u4-通道,无线或有线，高密度采样（实时数据监控和分析）

u分布式数据存储与回放，GPS控制所有时间，并提供观测位置信息

u使用远参考噪声压制，新型便携5kW高压发射

u使用电流监控和反褶积策略

u可以使用3D规则或随机网络观测和任意位置发射

二、勘探技术—仪器设备

智利SantaCecilia斑岩铜金矿床3D

DCIP-MT探测实例

§位于WesternCordillera,Maricunga

成矿带

§强烈热液蚀变

§地面磁测,CSAMT和ORION3D

DCIP/MTSanta

Cecilia

CerroDelMedio

----4600m

----3600m据N.Bournas等，2013

二、勘探技术—仪器设备

Qz-alunite-clay-PyPeripheral

propylitic

Silicification

Mz-intrusion

Intermediate

argillic

Advanced

argillic

Intensehydrothermalalteration

1km

据N.Bournas等，2013

二、勘探技术—仪器设备

地面磁测，线距：100米，显示环

状结构

垂向二阶倒数：高低梯度带大致

圈出蚀变范围

磁化率反演：蓝色区域指示因蚀

变而降低的磁化率

据N.Bournas等，2013

二、勘探技术—仪器设备

ORION3DSurvey(Fall2012)

•539Currentinjections(Transmits)

•300()receiverdipoles(150m)

•50DAUs(loggers)

Objective:3Dimagingofmineralization

andalteration

据N.Bournas等，2013

二、勘探技术—仪器设备

3D电阻率模型

据N.Bournas等N

1km据等，

Surveyresultsfromalarge-scale3-DchargeabilitysurveybyOuantec.N.Bournas2013

二、勘探技术—仪器设备

3D激发极化模型

据N.Bournas等N

1km据N.Bournas等，2013

二、勘探技术—仪器设备

据DNis.Bcoourvneasr等y，h2o01le3s(2011)Mineralizationover

1000m

Au0.2g/t

Cu0.25%

Mo80ppm

CSAMT

0

.

7

Discoveryholes

k

m(1400and1700m)

二、勘探技术—仪器设备

Astana3-DIP测量

据LiYG等，2013

二、勘探技术—仪器设备

AstanaDCIP3D结果

据LiYG等，2013ColorRange:Red>500,Purple<30ohmmeters

二、勘探技术—仪器设备

DCIP3D反演模拟结果

据LiYG等，2013

二、勘探技术—仪器设备

DCIP3D极化率模型

DCIP3D电阻率模型

据LiYG等，2013

二、勘探技术—仪器设备

3DDCIP+MT

加拿大Quantec公司Titan分布式

电磁系统：3DIP+MT

DCIP&MT的完

美组合

二、勘探技术—仪器设备

Titan24在萨德贝里Norman矿区试验结果，验证该技术有能力发现

1000－2000米深的矿层（Norman2000矿层）

山东地矿局总工研讨班—学术讲座之一

汇报提纲

1引言—两个原理模型

2勘探技术—仪器设备

3勘探技术—新方法

4勘探趋势：从矿床到成矿系统

/67/

三、勘探技术—新方法

1、金属矿地震技术：（1）物性前提

普通硫化物在200Mpa室温下的声学参数

密度波阻抗

矿物名称VpVs泊松比

（g/cm3）（km/s）（km/s）×105g/cm2s

Pyrite(黄铁矿)4.917.994.940.1939.2

Sphalerite（闪锌矿）4.075.442.800.3222.1

Chalcopyrite（黄铜矿）4.285.122.490.3521.9

Pyrrhotite（磁黄铁矿）4.714.602.730.2321.7

Pentlandite(镍黄铁矿)4.684.562.950.1421.3

Bornite（斑铜矿）5.074.402.070.3622.3

Galena（方铅矿）7.513.692.110.2627.7

Harvey,1997；Salisbury,etal.,1996

三、勘探技术—新方法

u黄铁矿与酸性脉石：

速度增加密度增加

u闪锌矿、黄铜矿、

磁黄铁矿与脉石：密

度增加，速度降低

u多数硫化物矿石，

相对于酸性岩石，可

成为反射体；一些甚

至对任意岩石，在几

何条件满足情况下，

硫化物矿石和硅酸盐岩石速度—密度关系曲线图，可称为反射体。

图上叠加波阻抗（Z）曲线

（Saliburyeral.,2003）

三、勘探技术—新方法

（2）分辨率问题。解决控矿构造，甚至直接探测矿体

波速（V）

波长（）=

频率（f）

垂向分辨率:

最小可分辨的构造，比如,薄

层的顶和底板

V

=50Hzwavelength

44f

LITHOPROBE反射地震的典型分辨能力:

震源频率：10–80Hz

速度范围：3–8km/sAfterClowes,2011

垂向分辨率范围：10–200m

三、勘探技术—新方法

横向分辨率：可分辨的水平两个反射体

Z取决于速度(v),

0+

频率(f)和深度(z)

4

Z0/

2zV

菲涅耳半径(r)=

f

FresnelZone

TypicalLITHOPROBEreflectionexperiment

•浅层地壳5km:<500m（水平分辨）

•Moho(35km):~2.5km（水平分辨）

三、勘探技术—新方法

频率:理论和实际探测分辨率

对4.3km/s的金伯利岩:

λ/4λ/20

频率波长分辨率理论极限探测极限

(Hz)(m)(m)(m)

5086224

10043112

2002251

300143.50.7

50092.250.5

问题：Approximate

Ø高频信号易衰减requirestuning

Ø震源与接收之间的耦合

Ø静校正(地形,横向结构和速度变化)

三、勘探技术—新方法

实例一：南非兰德盆地砾岩型金矿探测与开采设计

矿区名称：在数个矿区开展了几十到几百平方公里的3D反

射地震调查

uWUDL：WesternUltraDeepLevels:300km2

uTargetGoldProject:135km

3D反射地震：揭示结构框架、地层分布，追踪层状金矿层

（VentersdorpContactReef，VCR）；确定可能的矿

灾位置（断层及不稳定岩性分布）

属性分析：钻孔约束+地质观察+合成地震图：提供小断层

（10-20m）信息，砾岩层的岩性变化

三、勘探技术—新方法

CarletonvilleGF

WUDL

早太古代(37-31亿):古老的Kaapvaal克拉通(TTG)形成;绿岩型金矿

晚太古代(31-26亿):伸展作用形成Witwatersrand盆地;地层由老到新:

Dominion群；

Witwatersrand群,西兰德和中兰德组（砾岩型金矿主要层位）

Ventersdorp群（Klipviersberglava,Platberg组，Pneil沉积）

Transvaal群

三、勘探技术—新方法

WE

TargetGold

Project

从3D数据体中

抽出的东西向切片，

显示整体结构、构

造。

KRG是Klipriviersberglavas底部的反射，向西倾；下

伏的中兰德组在挤压变形中发生褶皱（向斜），发生逆冲和

背冲构造

三、勘探技术—新方法

Dreyeskuil砾岩层振幅从西到东的变化

Klipriviersberg熔岩底部Dreyerskuil

砾岩层的振幅水平切片(红色=高振幅,

蓝色=低振幅)。

注意：线状振幅低指示断裂延伸；东

南方向的低振幅表示顶部砾岩成分增加，

速度增大，造成与KRG的反射系数减

小；

三、勘探技术—新方法

实例2：反射地震技术寻找金伯利岩

—加拿大SnapLake地震项目

目的：发展探测低角度金伯利岩床喷出相

、岩席的新技术；降低钻探成本，对爆发相

岩席进行成像，确定：

Ø矿床的大致延伸

Ø为矿山开采设计提供结构依据复式岩席

管道相岩床

岩墙

AfterClowes,2011

三、勘探技术—新方法

granitoids

metavolcanics

矿山地质anorthosites,gabbros

metaturbidites

diabasedykes

faults

DeBeersCamp

D

ow

n-

dip

Seismicprofiles

Kirkleyetal.,2003

三、勘探技术—新方法

Line1:Down-dip(5.86km)

震源参数

•药量0.25kg

•480shots

•井深1-3m到基岩

•4-8-12-16-24mspacing

Explosives接收参数

•采样率0.5ms

•检波器主频28Hz

•816型单点检波器

•非对称变观观测系统

•2m道距:0-1.44km

•4m道距:1.44-5.86km

Vibroseis•最大偏移距:3.265km

AfterClowes,2011

三、勘探技术—新方法

Trackedrig

Towedor

Diamondheexli-supportedrigAfterClowes,2011

三、勘探技术—新方法

Kimberlite

Core

km

1.63m

从60米到1350米发现了

稳定的岩床反射MigratedStack

AfterClowes,2011

三、勘探技术—新方法

3D重磁反演

-500

泥河

小包庄

泥河铁矿纵1剖

面

重力异常小

-500包

庄

泥河铁

小包庄矿

2

线

剖

地磁异常面

长江中下游成矿带深部探测与找矿突

破

三、勘探技术—新方法

反射地震确定多级次火山岩的隆起部位，确定罗河铁

矿下面新的找矿目标

泥河

小包庄

上次火山岩

上次火山岩

下次火山岩

三、勘探技术—新方法

2、岩性填图技术：地球物理正反演

正问题：给定磁化率的分布kjj=1,M；计算位场值bii=1,

N

正演计算

反问题：给定观测值bii=1,N；寻求磁化率的分布kjj=1,

M

反演计算