CSAMT与AMT数据“拼接”处理——以湖南仁里铌钽矿床7号剖面为例

刘俊峰，程云涛，邓志强，周芳春，曹创华，刘翔，曾美强，李杰，黄志彪，陈虎

(1.湖南省核工业地质调查院，湖南 长沙 410011; 2.湖南省核工业地质局 311大队，湖南 长沙 410100; 3.湖南省地质调查院，湖南 长沙 410011; 4.湖南省生态环境事务中心，湖南 长沙 410014; 5.湖南省地质院，湖南 长沙 410014)

0 引言

AMT方法为音频段MT方法，通常使用的频率范围在0.01～10 000 Hz，勘探深度可达数千米(通常为2 km以浅)。AMT方法数据质量与高空大地电磁活动强弱有直接联系，北半球通常表现为夏强冬弱，夜强日弱，且一般在1～4 kHz 频带天然场信号强度通常很低，容易造成阻抗畸变，通常将这个频带称为“死频带”[1-2]。针对死频带干扰，A. Müller提出采用避开白天信号弱时段在夜间采集减小死频带干扰[3]；Gamble T等提出了用于改进采集点数据质量的远参考法[4]，可以在一定程度上减小畸变曲线，但通常参考点也受死频带影响，导致效果不佳；杨生和仇根根等提出相位校正法[5-6]，该方法要求相位质量相对较好，目前实际运用较少；Garcia X和 Jones A G提出采用连续小波变换的处理方法，但容易出现离散[7]；李红领等通过删除不满足相关度的时间序列去噪减小死频带干扰[8]，取得一定效果，但是时间序列的删选需要有一定的时间序列长度，否则处理效果将会受到影响。人工剔除飞点圆滑曲线也是“死频带”校正常用的方法[9]，但该方法通常人为影响很大，人工调整的畸变曲线不能真实反映实际情况。谭捍东等用Rhoplus校正方法处理大地电磁测深数据[10]；汤井田等和谭洁用Rhoplus校正方法处理校正“死频带”取得了较好的效果[11-13]；在改进观测方式方面，Garcia和Jones提出的T-MT法为一种基于大地电流法(TT)和电流—电磁法(T-MT)相结合的一种混合数据获取与处理方法[14]；乔宝强等展示了EH4仪器自带的场源增强信号的实例效果[15]。

CSAMT是利用人工场源的AMT方法，通过控制收发信号，能有效抑制噪声，减少“死频带”干扰。该法实际运用较多，但目前国内大部分仪器(GDP-32与V8)发射功率有限，收发距常设置偏小，频率进入近区过早，导致探测深度有限。针对近区校准，不少学者做了工作，但很少运用到实际数据处理当中，CSAMT的解释处理大多数情况下仍然是基于“远区”进行[16-17]。因此，如能结合AMT 与CSAMT各自优点，则有可能达到减少“死频带”干扰提高数据质量之目的。

已有数据表明，CSAMT的远区数据与AMT具有相似性，同测点的曲线形态基本一致[18]，CSAMT的远区数据与AMT数据拼接具可行性。本文通过均匀模型模拟计算证实上述两种方法响应数据拼接是可行的。以湖南仁里铌钽矿床7号剖面采集的CSAMT和AMT数据为例，展示了这种“拼接”的方法处理过程：将CSAMT近区和过渡带数据删除，同时将AMT中“死频带”畸变数据剔除，将两者拼接，经过静校正等数据处理后进行反演。与地质勘探剖面对比，表明这种数据处理方式对地电断面的浅层显示有一定的改善，且深部地电信息有更丰富的展示。

1 基本理论与数值模拟

1.1 基本理论

采用直角坐标系时，CSMAT远区各分量如下[19]：

(1)

Hy=Hrsinφ+Hφcosφ=

(2)

(3)

则视电阻率可以表示为：

(4)

(5)

式(5)也为AMT定义视电阻率公式。当收发距达到一定“足够”距离时，CSAMT远区数据与MT的偏差较小，可以近似等同[20-21]。

1.2 数值模拟

为了符合实际地层信息，曹创华等结合长株潭地区地层特点，进行了CSAMT收发距等相关实测参数的理论计算和实例分析[22]，得到了较好的结论。本文借鉴其思路，结合笔者团队在湘东北幕阜山仁里矿区实际观测的参数，设计CSAMT正演采用赤道偶极，收发距为5 000 m，电偶极源长为 1 000 m，接收点在中垂线上，正演频率为野外数据采集常用频率。

理论模拟计算采用均匀半空间电阻率模型，电阻率为2 000 Ω·m，采集频点共55个：分别为10 400、8 800、7 200、6 000、5 200、4 400、3 600、3 000、2 600、2 200、1 800、1 500、1 300、1 100、900、780、640、530、460、390、320、265、229、194、159、132、115、97、79、66、57、49、40、33、27.5、22.5、18.8、16.2、13.7、11.2、9.4、8.1、6.9、5.6、4.7、4.1、3.4、2.81、2.34、2.03、1.72、1.41、1.17、1.02、0.86 Hz。按照式(1～5)进行计算，得到正演模拟结果如图1所示。由图可见，CSAMT数据大于800 Hz为远区，AMT和CSAMT数据在远区电阻率与相位数值基本接近，误差较小，表明两种方法数据拼接具有可行性。

图1 理论模拟计算曲线Fig.1 Theoretical simulation calculation curve

2 实例验证

2.1 地质概况

实例数据在湖南仁里铌钽矿床7号剖面取得。该矿床为近年来华南地区发现的超大型、高品位花岗伟晶岩型铌钽矿，为我国华南地区新发现的超大型高品位花岗伟晶岩型铌钽矿，被评为2017年全国十大地质找矿成果之一。该矿床位于扬子陆块与华夏陆块交汇之江南隆起造山带中段北缘之湘东北断隆带，扬子陆块与华夏陆块的过渡部位，幕阜山岩体西南缘。区域上出露地层有青白口系、震旦系、寒武系、白垩系和第四系。矿区主要出露地层为冷家溪群，属于扬子陆块变质褶皱基底，为一套浅变质碎屑岩系，以片岩为主。区域构造表现为褶皱基底构造与盖层组成近东西向褶皱构造，断裂构造为燕山运动形成的NNE、NE向断裂构造[23-24]。

区内岩浆活动频繁，南部和东南部出露的雪峰期梅山岩体，梅仙岩体、三墩岩体、传梓源岩体和幕阜山岩体，其中以北面燕山期侵入为主的幕阜山大型复式岩体体量最大，产出面积为2 360 km2[25]。矿区地质简图见图2。

1—第四系;2—冷家溪群片岩;3—细粒花岗闪长岩;4—细粒二云母二长花岗岩;5—中粒似斑状黑云母二长花岗岩;6—粗中粒片麻状黑云母二长花岗岩;7—新元古代中细粒黑云母斜长花岗岩;8—伟晶岩脉及其编号; 9—主要断裂及编号;10—物探测线及编号1—Quaternary alluvium; 2—Lengjiaxi group schist; 3—fine-grained granodiorite; 4—fine-grained two-mica porphyritic biotite monzogranite; 5—medium-grained porphyritic biotite monzogranite；6—medium-grained gneissic biotite monzogranite; 7—Neoproterozoic two-mica plagioclase granite; 8—pegmatite and its serial number;9—main faultsand and its serial number; 10—measuring line and serial number图2 矿区地质与测线布置[24]Fig.2 Layout of geology and survey line in the mining area[24]

区内花岗质伟晶岩脉发育，主要呈脉状体产出于幕阜山复式岩体及其围岩冷家溪群片岩中，富含铌钽矿化的伟晶岩脉主要产于距离岩体接触带0.2～2 km 范围内的冷家溪群片岩中。典型的2、3、5、6 号脉走向总体呈NW向平行，倾向南西，规模大，其资源量占矿区Ta2O5资源量的97%，其中又以5号伟晶岩脉铌钽资源量最为显著，5号矿脉中的3个矿体(5-1,5-2，5-3)中Ta2O5资源量占全矿区Ta2O5资源量的67%[26-27]。

根据小比例尺航空磁测的研究推断[28-29]，幕阜山岩体西南缘(包含仁里段)存在隐伏磁性岩体，对中低温热液型多金属矿及与岩体相关的稀有金属矿产较为有利。相关研究[30-33]表明，含矿伟晶岩与花岗岩关系密切。

2.2 地球物理特征

矿区地层以冷家溪群浅变质碎屑岩系为主，主要为片岩，第四系在溪沟附近地表少量分布，厚度普遍较小。矿区北东为幕阜山岩体，南西为梅仙岩体。伟晶岩脉赋存于幕阜山复式岩体及其围岩冷家溪群片岩中。物性参数测试根据《湖南省区域地质志》相关物性参数章节[34]，并经地表对称小四极和岩心标本测试验证取得：第四系电阻率最低，为10～30 Ω·m；板溪群、冷家溪群变质片岩的电阻率为20～104Ω·m，大部分介于800～1 500 Ω·m之间。绝大部分伟晶岩和花岗岩电阻率接近，无统计学上明显差别，电阻率为中高阻—高阻，电阻率区间为102～106，大部分介于103～104之间，但值得注意的是少部分变质程度较深的变质岩与花岗岩(伟晶岩)电阻率接近。地表风化严重的伟晶岩和花岗岩(松散状)物性接近第四系电阻率。

矿区地质—地球物理特征与四川甲基卡锂矿田伟晶岩脉具有可类比性，即稀有金属矿床含矿体密度低而电阻率高，变质岩围岩电阻率低导电性高的特点。利用伟晶岩与围岩(片岩)之间电阻率的差异可指导找矿[35-36]。仁里矿区围岩为变质岩，含矿地质体为花岗质伟晶岩，具有开展电磁探测的前提条件，即可借助地电探测区分中低阻的片岩和中高阻的伟晶岩(或花岗岩)。

2.3 工作布置

针对典型剖面勘探7线进行了物探测线布置，物探测线与地质勘探剖面线重合，同时垂直于伟晶岩脉走向，测线方位角为40°。定义测线方向为y方向，垂直y方向为x方向, CSAMT观测方式为Ey/Hx模式(本测区为TM模式)，AMT的Ey/Hx模式与CSAMT的Ey/Hx模式电磁道数据方向一致。CSAMT采用GDP-32仪器进行采集，仪器采用电偶极源，发射偶极子长1 000 m，发射机功率为10 kW，发射频率为8 192～1 Hz，收发距约5 000 m，采集时间为2018年5月。AMT采用V8多功能电法工作站，采用矢量观测方式，观测时长为30～45 min，野外施工时间为2019年8～10月白天，根据已有资料,该段时间为北半球信号天然电磁场信号弱的时段，容易造成“死频带”干扰。AMT起点位置为0 m，终点位置为4 520 m，起点和终点位置均与岩体相接，起点为梅仙岩体，终点为幕阜山岩体；CSAMT测量位置与AMT测点在1 600～4 400 m段为重合位置，两种方法点距和极距均为40 m。测线布置见图2。

2.4 实测CASMT与AMT典型曲线特征

实测CSAMT由于收发距短，且由于该区地下整体电阻率较高，信噪比较高，曲线整体较圆滑，大部分曲线在1 200 Hz开始进入过渡带(图3)。取平均电阻率为2 000 Ω·m，可得勘探深度约为500 m。

图3 2 360 m测点处CSAMT电阻率与相位曲线Fig.3 Resistivity and phase curves of CSAMT at 2 360 m

AMT曲线整体质量较好，大部分测点在500～5 000 Hz频段，ρxy和ρyx均出现视电阻率曲线畸变的趋势[11]，主要表现为曲线不光滑，曲线鲁棒性变差。阻抗相位曲线在该频带亦出现不光滑现象。在功率谱上表现为该频段信号强度低，考虑为该段天然场信号弱，以2 360 m测点为例，实测曲线见图4。

图4 2 360 m测点处AMT电阻率与相位曲线Fig.4 Resistivity and phase curves of AMT at 2 360 m

2.5 拼接方法

数据拼接采用同测点拼接方式，步骤如下：

1) 取CSAMT远区数据，将过渡带和近场区数据剔除。

2) 将AMTEy/Hx模式死频带数据去除。

3) 将编辑后的CSAMTρyx与编辑后同方向的AMT数据ρyx数据拼接，相位拼接方式相同。

4) 静态校准：因多种因素影响，比如地形起伏、场源效应，数据采集时间段不同(地表岩土干湿状态不同)接地电阻不一致，CSAMT部分测点远区视电阻率与AMT相比呈现整体偏大或偏小情况，对此，可以整体平移到大体一致，然后进行5点滤波校准。

5) 数据插值：本次CSAMT 远区数据大部分为1 200～9 000 Hz频率段，而AMT“死频带”主要集中在500～5 000 Hz频率段，因此Ey/Hx模式方向数据拼接后可能仍然存在一定的数据空白段(如500～1 200 Hz)，可以适当用插值来解决。这时可以使用Ey/Hx模式数据进行反演(本次为TM模式)。

为了更加有效利用AMTEx/Hy模式数据，AMTEx/Hy模式“死频带”数据结合Ex/Hy模式曲线形态及同频段Ex/Hy模式数据进行人工调整编辑，使曲线较光滑，此段数据调整主要根据经验判断。拼接后数据示例见图5，可以看出，拼接后数据形态整体较好，较有效地减少了“死频带”的干扰。

图5 2 360 m测点处数据拼接曲线Fig.5 Data splicing curve at 2 360 m measuring point

2.6 反演成果及地质解释

图6 7线拼接数据反演结果Fig.6 Inversion diagram of line 7 mosaic data

地电成果剖面整体展现为中高阻—高阻，浅地表电阻率较低，深部电阻率大于地表。在剖面起始端0～300 m和4 000～4 520 m段电阻率呈现高阻，并往深部延伸，分别对应梅仙岩体和幕阜山岩体。300～4 000 m段电阻率总体可分为低—高—低—高结构：剖面标高0～500 m深度，浅地表为低阻，对应于片岩；中间深度有数十个向上延伸的“脉状”高阻侵入体，推测为侵入伟晶岩花岗岩；伟晶岩表现为中低阻，推测为Pt板溪群片岩并夹杂有伟晶岩(花岗岩)通道。

其中2 500～3 000 m段位置，“脉状”高阻与已知勘探剖面(图7)有一定的对应关系：ZK716、ZK708、ZK704的片岩，对应于剖面的中低阻，伟晶岩与花岗岩对应于地电剖面的中高阻。伟晶岩在反演图上分辨不明显，体现在高低阻的接触过渡带上。结合上面物性参数分析可知伟晶岩与花岗岩物性参数相近，无法区别，但结合地质解释，可以大致判别片岩与岩脉的界线；根据如图7b所示勘查成果[23]，AMT与CSAMT数据拼接处理后对应的地电剖面如图7a所示。可见拼接处理后地电剖面能较好地呈现伟晶岩(或花岗岩)形态，对该区勘查工作有较大的指导意义。

1—冷家溪群片岩;2—花岗岩;3—伟晶岩脉及其编号;4—矿体及编号;5—钻孔及编号;6—槽探及编号1—Lengjiaxi group schist; 2—granites; 3—pegmatite and its numbered; 4—ore body and its numbered; 5—drilling engineering and its numbered; 6—pit engineering and its numbered图7 7线已有钻孔段地电剖面(a)和地质勘探剖面(b)[23]Fig.7 Geoelectric profile of drilling verification section (a) and geological exploration section(b) of line 7 [23]

3 结论与讨论

本文从已有理论出发，分析了CSAMT与AMT数据拼接进行数据处理反演的正确性，表明CSAMT远区数据与AMT数据拼接是可行的。在实例上展示了数据拼接处理的过程和处理后的反演效果,利用在仁里7号勘探剖面采集的CSAMT与AMT数据，对AMT数据“死频带”数据去除和CSAMT近区数据剔除后拼接，通过静校正等处理手段后反演取得了较好的效果，在浅部地质体对照和深部岩体探测方面均取得效果，表明该组合方法既结合了CSAMT浅层分辨较高，AMT勘探深度大的优点，同时不受AMT“死频带”和CSAMT近区干扰的影响。

值得注意的是，受场源效应(CSAMT法)、收发距(CSAMT法)、地形(CSAMT法和AMT 法)等影响，实际数据与AMT仍然存在一定的偏差，因此实际数据处理过程中需要认真分析数据造成偏差的原因后进行处理，避免“硬拼接”。