强干扰矿区中广域电磁法与CSAMT探测效果对比
李帝铨, 胡艳芳
中南大学 地球科学与信息物理学院,湖南 长沙 410083

作者简介: 李帝铨(1982-),男,副教授,博士,主要研究方向为电磁法正反演。E-mail:lidiquan@csu.edu.cn

摘要

安徽某铜矿区由于电磁、人文干扰严重,传统的电磁法无法取得有效的数据,抗干扰能力较强的CSAMT法在该区检查点数据均方相对误差达到了53.6%,勘探效果较差。将广域电磁法应用到该矿区,取得了质量较高的数据,平均相对均方误差为5.6%。与相同剖面上的CSAMT法反演结果相比,广域电磁法的反演结果与地质资料吻合,深度误差小于10%,矿体的位置与钻孔信息相符。试验结果表明,广域电磁法具有较强的抗干扰能力,可用于强干扰矿区探测。

关键词: 广域电磁法; 可控源音频大地电磁法; 抗干扰能力; 铜矿
中图分类号:P631 文献标志码:A 文章编号:1000-8918(2015)05-0967-06 doi: 10.11720/wtyht.2015.5.15
A comparison of wide field electromagnetic method with CSAMT method in strong interferential mining area
LI Di-Quan, HU Yan-Fang
School of Geosciences and Info-physics,Central South University,Changsha 410083, China
Abstract

Because of the severe electromagnetic interference and the human interference in a copper mine in Anhui Province, the traditional electromagnetic method is unable to obtain valid data. CSAMT method has stronger anti-interference capability, but the mean square relative error of CSAMT checkpoints data reaches 53.6% and the exploration effect is poorer. In this paper, wide field electromagnetic method was applied to the mining area and led to the obtaining of the high quality data, with the average relative mean square error being 5.6%. Compared with the inversion result of CSAMT method, wide-area electromagnetic inversion result is in accordance with the geological data, the depth error is less than 10%, and the position of the orebody is consistent with the drilling information. The test results show that wide field electromagnetic method has stronger anti-interference capability and can be used in strong interference field detection.

Keyword: wide field electromagnetic method; CSAMT; anti-interference capability; copper deposit

安徽省某铜矿埋深大、规模大, 似层状, 其规模和矿体分布范围基本上由地表钻孔控制。矿区内电力线、民房、管线横跨, 部分矿体所在位置的地面被民用与矿山建筑所覆盖, 矿床南段正在施工开采, 电动车的运行产生很大的电磁场噪声, 矿区的人文干扰强而复杂, 而现有的商用仪器(如GDP32、V8等)抗干扰能力有限, 在采矿区上方, 均不能获得符合规范要求的数据, 所有商用软件的抗干扰方法均为圆滑、滤波, 对矿区通过软件也不能消除干扰, 勘探效果差。

广域电磁法是由何继善教授提出的一种新的人工源频率域电磁测深法[1, 2], 该方法采用伪随机 1214信号发射, 在保证观测信息量的前提下, 只需要观测电磁场的一个分量, 因此具有探测速度快、成本低、抗干扰能力强、勘探深度大的优点。为了比对广域电磁法的抗干扰能力、强干扰地区数据采集的可靠性, 以及反演成像对地下电性结构的响应能力, 中南大学在该矿区开展了广域电磁法试验工作。

1 试验区概况
1.1 地质概况

矿区内地表除部分为第四系覆盖层外, 还出露有三叠系中、下统, 局部有上统零星分布。深部工程揭露的地层有二叠系、石炭系、泥盆系上统等地层[15]

矿区内岩浆岩比较发育, 以中酸性小型岩体及浅成岩脉为主。矿床的生成与岩浆岩活动紧密相关, 在岩浆岩与钙质岩类接触部位往往出现矽卡岩化, 矿体赋存在矽卡岩内[16]

矿床主矿体严格受容矿层位岩性控制, 属层控矽卡岩型铜矿床。在矿床主矿体中, 矿石类型复杂, 其中脉石矿物主要有蛇纹石、滑石等蚀变矿物和石榴石、透辉石等矽卡岩型矿物; 矿石矿物主要有磁黄铁矿、黄铁矿、黄铜矿和磁铁矿等。主矿体赋存于某背斜的轴部及两翼, 严格受石炭系中、上统黄龙— — 船山组层位控制, 形态简单, 呈似层状, 产状与控矿岩层近于一致。矿体在空间上以背斜隆起部位的赋存标高最高, 呈一个不完整的穹隆状, 沿走向及倾角均显舒缓波状起伏。

1.2 地球物理条件

表1所示, 矿区除二叠系地层以外, 其余地层具有高阻、低极化的电性特征, 电阻率平均值在 1 000 Ω · m, 极化率为2%~3%。而二叠系大隆组和孤峰组电阻率为几十到几百Ω · m, 极化率大于40%。这两组地层可视为电磁法勘探的最大地质干扰层。侵入岩(闪长岩类)为中等电阻率, 多为几百Ω · m, 极化率为2%~4%。矿体电性随矿化性质及矿化程度而异, 主要表现为低阻高极化特征, 电阻率一般小于100 Ω · m, 极化率大于60%。

表1 测区岩(矿)石物性参数
1.3 电磁干扰分析

由于历史原因, 该矿区的生活区和矿山作业区已基本相联在一起。随着人类活动的频繁, 其激发的电磁场在空间分布上日益增加, 在时间域上覆盖加大, 在频率域上频点增多, 而且不同频谱的信号能量也越来越强。主要表现为50 Hz的工频干扰, 其中矿区的输电高压线、井口附近的大型用电设备、巷道运矿石车的供电线是主要的磁场干扰源, 用电设备的接地线是主要的电场干扰源。在该矿区进行广域电磁法试验前, 采用V5-2000记录了该区的天然场时间序列信号, 经分析, 该区的强干扰类型主要包括:尖脉冲干扰、类阶跃干扰和似放电干扰等[17]。在这种强干扰情况下, 无论是电场信号还是磁场信号, 几乎都淹没在50 Hz工频干扰信号中(图1)。

图1 V5-2000记录的天然电磁场时间序列信号

2 试验剖面与试验参数选择

为了对比广域电磁法与CSAMT法在该铜矿区的勘探效果, 选择了曾经开展过CSAMT勘探的BD-1号剖面进行广域电磁法试验。该剖面由6个钻孔揭露, 为已知地质剖面, 具备开展广域电磁法与CSAMT法对比的试验前提。

图2为BD-1号剖面101号测点所测得的广域电磁测深降流数据对比图, 本次降流试验采用7频波发射, 最大电流分别为100 A、50 A、20 A和10 A。如图2所示, 广域电磁法在发射电流为100 A、50 A以及20 A的情况下, 观测到的电阻率曲线形态一致。相比之下, 当发射电流为10 A时, 所观测到的数据不太稳定, 在低频段, 曲线形态变化较大, 电阻率曲线整体有上移趋势, 主要是由于电流较小, 难以压制干扰引起, 整体曲线形态变化不大。

图2 101点广域电磁法降流数据对比

通过本次降流试验可知, 广域电磁法采用伪随机信号发射, 保留了频率域测总场要求的供电电流不大的优点, 即使采用7频波, 电流减小至20 A时, 所测得的广域视电阻率数据仍然可靠。

3 试验结果对比
3.1 原始曲线对比

为试验广域电磁法在该强干扰区的勘探效果, 在该剖面上布设了广域电磁法剖面, 测点距为 50 m, 测线完全重合前人开展CSAMT的勘探线。

前人在该剖面分别采用GDP-32和V8两种仪器进行CSAMT勘探。GDP-32的工作频率范围为4~4 096 Hz, 收发距为4 800 m, 发射电流为17 A, MN极距为50 m; V8的工作频率范围为1~7 680 Hz, 收发距为8 246 m, 发射电流为18 A, MN极距为100 m。为了和前人在该矿区开展的CSAMT进行有效的对比, 本次试验根据前人开展CSAMT勘探时的测点位置坐标, 进行广域电磁法野外数据采集。试验装备为中南大学自主研发的广域电磁法观测系统, 施工采用伪随机信号的7频波的0、1、2、3频组供电, 工作频率为0.75~8 192 Hz, 共28个频点, 发射电流100 A, 发收距为7 850 m。

图3所示为BD-1剖面GDP32测得的检查点卡尼亚电阻率[20]与广域电阻率对比。GDP32在该测点获得的检查点与初测曲线的形态基本相同, 但呈锯齿状, 相邻频点的电阻率变化剧烈且无规律, 数据可靠性低; 数据误差较大, 视电阻率均方相对误差达53.6%, 曲线整体产生位移。这主要是由该矿区的强电磁干扰以及复杂的人文干扰所引起。相比之下, 广域电磁法采用伪随机信号发射, 具有抗干扰能力强的优点, 获得的广域视电阻率曲线圆滑, 电阻率数值与工区的电性特征对应。

图3 88点广域电阻率与卡尼亚电阻率对比曲线

从图4可以看出, V8仪器在BD-1号剖面102点所观测到的检查数据与初始数据除部分畸变点外, 吻合度较高, 曲线形态一致, 但电阻率范围为106~108 Ω · m, 与工区的实际电性特征不符。广域电磁法所测的视电阻率相对于CSAMT法所测的电阻率曲线要相对平滑, 基本上没有畸变点, 电阻率数值与工区的电性特征相对应。

图4 102点广域电阻率与卡尼亚电阻率对比曲线

3.2 反演效果对比

图5为BD-1线地质剖面和CSAMT勘探与广域电磁法勘探获得的电阻率二维反演结果[20]。地质剖面为钻孔控制, 已知信息较多, 矿体的位置与钻孔附近地层信息准确。前人推断在86~96点间900 m以浅有一较宽的低阻区, 900 m以下此低阻区主体向左下方倾斜延伸, 等值线陡立密集, 推断此处有一断裂破碎带存在, 但对矿体位置的推断过于牵强。反演结果整体而言, 对矿体几乎没有反映, 地层的似层状特性也没有得到体现(图5b)。

图5 BD-1测线地质剖面及2种方法反演结果

图5c为BD-1线的广域电磁法电阻率二维反演剖面。地层的似层状体现非常明显, 矿体横向和深度范围内出现了低阻异常, 识别出了两条断层, 分别位于86~96点和116~119点之间。在82~101点1 000 m以下电阻率较低, 推断为石炭系地层。该矿区铜矿体受控于泥盆系和石炭系之间的层间滑脱构造带, 主要赋存于石炭系地层内, 与已知信息一致。推断测线小号点还存在成矿的空间, 该处低阻异常可能为矿体或矿化引起。

通过对比可以看出, 在已知矿体位置上广域电阻率反演结果呈现明显的低阻异常, 相比之下, CSAMT法沿用了CAGNIARD视电阻率的计算公式, 对低阻的响应相对较弱。图5c中异常区域的位置与实际地质情况吻合较好, 这说明广域电磁法能较好地反映该剖面的主要矿体位置, 也说明该方法在该区勘探是可行的。对比2种方法的勘探结果可发现, 相对于CSAMT法而言, 广域电磁法具有抗干扰能力强的优点, 同时也说明在强干扰矿区, 广域电磁法能获得比CSAMT法更可靠的数据和勘探效果。

4 结论

(1) 可控源音频大地电磁法采用波区视电阻率定义, 计算简便, 但在非波区会产生畸变, 影响测深曲线的解释。广域电磁法从场的统一性出发, 采用全区视电阻率定义, 改善了非远区的畸变效应, 使得测深能在广大的、不局限远区的全域进行, 在同等收发距上增大了勘探深度。

(2) 广域电磁法采用伪随机信号发射, 保留了在频率域观测总场供电电流不大的要求, 即使采用7频波, 电流减小至20 A时, 所测得的数据还是可靠的, 而且由于多频波同时供、测, 提高了观测精度, 提高了抗干扰的能力。

(3) 广域电磁法只需要测量电场数据进行视电阻率的计算, 从而避免了磁场信号的干扰, 拓展了人工源频率域电磁法的观测范围, 提高了观测速度、精度, 加大了勘探深度。

(4) 可控源音频大地电磁法电阻率测深曲线在高频段和低频段出现了畸变点, 主要是由电磁干扰引起的, 相比之下, 广域电磁法电阻率曲线平滑。说明在强干扰矿区, 广域电磁法采集的数据质量较高。

(5) 相对于CSAMT, 广域电磁法反演结果清楚地反应了矿体的位置和埋深。广域电磁法是探测大深度复杂的、隐伏矿体的有效工具。

The authors have declared that no competing interests exist.

参考文献
[1] 何继善. 广域电磁法和伪随机信号电法[M]. 北京: 高等教育出版社, 2010: 1-7. [本文引用:1]
[2] 何继善. 广域电磁测深法研究[J]. 中南大学学报: 自然科学版, 2010, 41(3): 1065-1072. [本文引用:1]
[3] 汤井田, 何继善. 可控源音频大地电磁法及应用[M]. 长沙: 中南大学出版社, 2005: 250-285. [本文引用:1]
[4] 刘春明, 佟铁钢, 何继善. 多种电磁法在某金矿的野外勘探应用[J]. 中国有色金属学报, 2013, 23(9): 2422-2429. [本文引用:1]
[5] 何继善. 可控源音频大地电磁法[M]. 长沙: 中南工业大学出版社, 1990. [本文引用:1]
[6] 李帝铨, 谢维, 程党性. E-Ex广域电磁法三维数值模拟[J]. 中国有色金属学报, 2013, 23(9): 2459-2470. [本文引用:1]
[7] 于昌明. CSAMT方法在寻找隐伏金矿中的应用[J]. 地球物理学报, 1998, 41(1): 133-138. [本文引用:1]
[8] 张国鸿, 李仁和. 可控源音频大地电磁法深部找矿实验效果[J]. 物探与化探, 2010, 34(1): 66-70. [本文引用:1]
[9] 邓锋华, 杨洋, 李帝铨. 广域电磁法在隐伏金矿中的应用[J]. 工程地球物理学报, 2013, 10(3): 357-362. [本文引用:1]
[10] 王大勇, 李桐林, 高远, . CSAMT法和TEM法在铜陵龙虎山地区隐伏矿勘探中的应用[J]. 吉林大学学报: 地球科学版, 2009, 39(6): 1134-1140. [本文引用:1]
[11] 汤井田, 葛伟男. 三维CSAMT 中的阴影和场源附加效应[J]. 物探化探计算技术, 2012, 34(1): 19-27. [本文引用:1]
[12] 何继善, 柳建新. 伪随机多频相位法及其应用简介[J]. 中国有色金属学报, 2002, 12(2): 374-376. [本文引用:1]
[13] 何继善. 三元素集合中的自封闭加法与2n系列伪随机信号编码[J]. 中南大学学报: 自然科学版, 2010, 41(2): 631-637. [本文引用:1]
[14] 何继善, 佟铁钢, 柳建新. an序列伪随机多频信号数学分析及实现[J]. 中南大学学报: 自然科学版, 2009, 40(6): 1666-1671. [本文引用:1]
[15] 杨爽. 安徽铜陵冬瓜山层控矽卡岩型铜矿床特征和成因[D]. 北京: 中国地质大学, 2012. [本文引用:1]
[16] 蔡运胜, 张进国, 杨学明. EH4电磁测量在安徽冬瓜山铜矿区找矿对比中的应用效果[J]. 矿产勘查, 2011, 2(1): 75-80. [本文引用:1]
[17] 汤井田, 徐志敏, 肖晓, . 庐枞矿集区大地电磁测深强噪声的影响规律[J]. 地球物理学报, 2012, 55(12): 4147-4159. [本文引用:1]
[18] 闫述, 陈明生. CSAMT测深中的阴影和场源复印效应问题[J]. 石油地球物理勘探, 2004, 39(S): 8-10. [本文引用:1]
[19] 陈明生, 闫述. CSAMT勘探中场区、记录规则、阴影及场源复印效应的解析研究[J]. 地球物理学报, 2005, 48(4): 951-958. [本文引用:1]
[20] 中国地质科学院地球物理地球化学勘查研究所. 危机矿山电(磁)方法有效性比对成果报告[R]. 2012: 3-48. [本文引用:2]