不同类型低阻异常航电时间常数的特征分析
骆燕1,2, 江民忠1,2, 宁媛丽1,2, 彭莉红1,2, 朱琳1,2
1.核工业航测遥感中心,河北 石家庄 050002
2.中核集团 铀资源地球物理勘查技术中心重点实验室,河北 石家庄 050002

作者简介: 骆燕(1981-),女,工程师,2007年毕业于东华理工大学,获硕士学位,主要从事航空电磁法及航空磁法应用与研究工作。

摘要

时间域航空电磁系统工作原理是机载发射线圈向地下发射脉冲电磁波,接收线圈观测断电后来自地下低阻体的二次感应电磁场,通过对观测数据进行处理解释达到探测地下低阻体的目的。结合近年来从事时间域航空电磁法探矿数据和资料解释的实践,分别对电力干扰、表层强低阻体、表层弱低阻体和地下低阻体航电异常响应的时间常数 τ的特征进行了综合分析, 总结出以下规律:电力干扰在 τ图上一般显示为不连续的线性异常,电磁响应呈跳跃状的无规律衰减;低阻异常在 τ图上一般显示为高值区,电磁响应呈指数衰减,导电性越好规模越大的低阻体衰减越慢。

关键词: 低阻异常; 时间常数; 时间域航空电磁法; 电磁响应; 感应磁场
中图分类号:P631 文献标志码:A 文章编号:1000-8918(2016)05-0991-07 doi: 10.11720/wtyht.2016.5.24
The characteristics of time constant of different conductive anomalies from airborne TEM data
LUO Yan1,2, JIANG Min-Zhong1,2, NING Yuan-Li1,2, PENG Li-Hong1,2, ZHU Lin1,2
1. Airborne Survey and Remote Sensing Center of Nuclear Industry, Shijiazhuang 050002, China
2. Key Laboratory for Geophysical Exploration Technology Center of Uranium Resource CNNC, Shijiazhuang 050002, China
Abstract

The principle of airborne time domain electromagnetic system is that Tx coil transmits electromagnetic pulse into ground, and Rx coil receives the induced secondary EM fields from conductive bodies in the earth. The secondary fields generated by induced current of conductive bodies contain geological information. Data processing and interpretation help to achieve prospecting for multi-metal mineral resources. The authors analyzed features of time constant reflection to power line interference under the conditions of near surface conductive bodies, near surface weak conductive bodies, and deep conductive bodies respectively and summarized the characteristics of time constant of different geological bodies based on knowledge of known mineral formation and two years' experience in airborne EM surveys. Some conclusions have been reached: the power line interference generally shows a discontinuous linear anomaly in τ map and electromagnetic response is jumping like irregular attenuation; low resistivity anomaly generally shows high value area in τ map and the electromagnetic response decays exponentially. The low resistance body decays more slowly when the conduction is better and the body is bigger.

Keyword: conductive anomalies; time constant; airborne time domain electromagnetic; electromagnetic response; induced magnetic field

时间域航空电磁法(航空瞬变电磁法)是以飞行器为载体, 通过机载或在空中拖曳发射和接收装置进行探测的方法, 是航空电磁法的一个重要分支[1]。该方法测量速度快、效率高、受地形影响小, 是多金属矿产勘查的重要方法之一[2, 3], 在国外得到了广泛的应用。随着我国矿产资源勘查难度的加大, 新方法新技术的应用变得尤为重要。2013年核工业航测遥感中心与加拿大Geotech公司合作在已知VMS矿区开展航电测量试验, 试验效果显著[4, 5]。2014年, 在黑龙江、内蒙、新疆等地成功开展了航空电磁测量工作; 2015年初完成的青海航电测量, 是我国5 km以上高海拔地区第一个高精度航空时间域电磁测量项目, 在世界范围内同类地区也尚属首次。探测结果表明航空电磁异常与已知矿区具有较高吻合度, 新发现了数十处具有找矿意义的航电异常。

航电VTEM系统观测的数据是z分量和x分量的dB/dt[6](感应磁场随时间变化率, 简称磁场导数)和B场(感应磁场, 简称磁场)[6], 同时通过“ 移动窗口” 技术[7]计算z分量dB/dtB场的电磁响应时间常数[8, 9]。通过分析研究不同低阻体和干扰引起的航电异常响应和τ 特征, 指导航电资料的解释, 为圈定成矿靶区提供基础。

1 航空瞬变电磁法时间常数的理论和计算
1.1 时间常数的物理意义

根据地下电磁波的传播机理, 供电产生的一次磁场在传播过程中, 如遇到地下导电体, 将在其内部激发产生感应电流, 又称涡流或二次电流。二次电流在其周围又产生新的磁场, 称为二次磁场。由于良导电体内感应电流的热耗损, 二次场大致按指数规律随时间衰减。二次场主要来源于良导电体的感应电流, 因此它包含着与矿体有关的地质信息。根据二次场的衰减规律进行分析处理, 可以解释地下矿体及相关物理参数[10]

矿体内感应电流及与之对应的二次场随时间变化率, 取决于矿体的岩性和物性参数(如电导率、磁导率、矿体尺寸大小及形状等)。瞬变电磁曲线变化速率的大小, 取决于异常体的电性和几何形状, 其电性和几何形状反应在它的时间常数上, 因此, 异常体的时间常数反映瞬变场衰减速率的大小, 是推断异常体电性的重要参数。几何形状相同的矿体, 导电性越好, 其时间常数越大。一般导体的时间常数可表示为[3]:

τ=KAσμ,

式中:K为随异常体而变的常数, A为与异常体几何形状有关的参数, σ 为电导率, μ 为磁导率(常以真空磁导率μ 0近似代替, 单位ms)。

根据电磁感应原理, 导体磁通量的变化会产生与导电体一次场磁场变化率成比例的电动势, 此电动势会在导体中形成一种特殊的瞬变衰减涡流。其时间常数τ 是目标体导电率(或导电性)和几何形态(规模)的函数。衰减的电流会产生相应的二次场, 其变化率可以由接收线圈测量断电后的感应电压而得到。

接收线圈测到的感应电压与二次场时间变化率成比例:

V(t)(1/τ)e-t/τ

其中:τ =L/R, 为时间常数, R为电阻, L为电感。

根据加拿大多伦多大学A.V.Dyck等人的研究, 瞬变场晚期衰减规律其衰减曲线可近似表示为

A(t)=A0e-(t/τ)

由上式可以看出, 异常体瞬变电磁响应的异常幅度、衰减速度以及时间常数τ 都是表征瞬变电磁响应特征的参数。对于良导体, 二次场初值较小, 但衰减速度慢, 时间常数τ 值较大; 反之, 对于不良导体, 二次场初值较大, 但其衰减速度较快, 时间常数τ 值较小。如图1所示[7]

图1 良导体和不良导体的dB/dt衰减曲线

1.2 时间常数τ 的计算

航空瞬变电磁法时间常数τ 既可以用dB/dt衰减数据来计算, 也可以由B场衰减数据来计算, 但是dB/dt信噪比更好更稳定。一般情况下, 早期时间道响应的τ 反映的是近地表覆盖层及不良导体, 晚期道时间响应的τ 反映的是埋藏相对较深的目标体或良导体。因此, 根据dB/dt计算的晚期时间响应的τ 可以帮助推断航电地下低阻异常区。

时间常数τ 的计算方法有很多种, 对于VTEM系统观测的数据, 通常采用“ 移动窗口” 技术(Geotech公司自编软件)来计算dB/dtB场的电磁响应时间常数τ 。计算原理是沿电磁响应的衰减曲线滑动时间窗口, 根据噪声水平和预设的信号阈值水平, 确定晚期的4个时间道(图2); 然后对这4个时间道进行最小二乘拟合, 获得时间常数τ

图2 VTEM数据中典型的dB/dt衰减

2 航电异常响应的时间常数特征

地下不同深度, 相对复杂的地质构造, 其时间常数τ 是变化的。在瞬变电磁晚期分时窗进行时间常数求取, 得到晚期不同时间点, 即不同深度处的时间常数, 以得出在晚期对应深度上的异常体分布规律。由于良导电性异常体的响应有延迟效应, 晚期道时间常数不但有时间窗内异常体反映, 同时还包含了部分早期良导电性矿体响应信息, 因此时间常数是一个地下区域所有矿体的综合响应。

根据近两年的航电测量成果, 分别对电力干扰、表层强低阻体、表层弱低阻体和地下低阻体航电异常响应的时间常数τ 的特征进行了综合分析。

2.1 电力线干扰响应的时间常数特征

航电测量通常会在具有一定勘查程度、找矿前景较好的区域或已知矿区开展, 这就难免会遇到工作区内存在电力线的情况。图3为某工作区z分量dB/dt计算的时间常数τ 等值平面, 由图可以看出该区有一NNW向线性异常带。由电力线电压监测平面图(图4)可以看出, 该线性异常带为电力线干扰的异常带。

图3 某区时间常数τ 平面

图4 某区电力线电压监测平面

将有电力干扰异常的衰减曲线和正常响应的衰减曲线进行对比(图5), 可以看出干扰异常上的衰减曲线较为跳跃, 尤其是中晚期时间道, dB/dt响应在零值上下剧烈跳动; 而正常响应的曲线较为圆滑, 呈指数衰减, 该处无低阻异常, 因此衰减较快。

图5 dB/dt衰减曲线

由上可知, 电力线干扰的区域通常有如下特征:

在时间常数τ 等值图上呈线性异常带分布; 通过对比正常响应的dB/dt衰减曲线, 干扰异常衰减曲线, 响应较为跳跃特别是在中晚期时间道。在今后资料解释中如遇到该特征的航电异常需谨慎对待。

2.2 地下低阻体和表层低阻体的时间常数特征

内蒙古小营子铅锌矿出露地层为温都尔庙群哈尔哈达组、二叠系浅变质岩系。其中, 温都尔庙群哈尔哈达组条带状大理岩是主要的赋矿围岩。岩浆岩主要有中侏罗世中酸性中深成相— — 超浅成相侵入岩。中深成侵入体斜长花岗岩和石英闪长岩分布于矿区东西两侧, 呈北西撒开, 南东收拢的“ 蟹螯状” , 向深部两岩体倾斜变缓, 呈不对称状。沿走向收拢或沿倾斜变缓翘起部位, 均为成矿有利地段, 分别赋存富厚矿体或矿脉群。

该已知矿床为渗滤式矽卡岩型矿床, 主要蚀变岩石为绿帘石透辉石岩、绿泥石绿帘石岩及角岩, 主要金属矿物为方铅矿、闪锌矿, 其嵌布颗粒较粗、属易选硫化物矿石[11]

在该矿区上方完成3条测试线飞行, 测线间距为500 m。图6为根据z分量dB/dt计算的时间常数τ 等值平面。图中高亮区为τ 高值区(> 1.1 ms), 该区为连续的高导电区, 与已知矿区范围吻合良好, 推测为铅锌矿区范围。由于该矿正在开采, 因此推测东南测线上快速变弱的τ 值区为铅锌矿的尖灭区或采空区[12]

图6 小营子矿时间常数τ 等值线平面

图7为L1线的dB/dt曲线和视电阻率深度成像(RDI), 中间是地下低阻体(2号), 两侧有2个表层低阻体(1号、3号)。从RDI上可以看出每个目标体有不同的埋深。3个目标体都是良导体, 相对地下低阻体, 两侧表层低阻体的dB/dt信号早期道很强, 随时间衰减迅速, 晚期道相对较弱, 时间常数相对较小, 也只有地下低阻体的1/14。该实例突出了时间常数τ 分析在分辨不同目标体中的优势。

图7 小营子矿区L1线dB/dt曲线(上)和视电阻率深度成像示意(下)

图8为上述3个异常的dB/dt衰减曲线。由图可以看出:1号和3号异常早期道幅值相对较大而衰减较快, 2号异常早期道幅值偏小, 但衰减慢。对比不同异常衰减曲线发现, 地下低阻体上具有早期道幅值偏小、衰减较慢的特征; 而表层低阻体具有早期道幅值相对较大、但衰减较快的特征。

图8 L1线3个异常上的dB/dt衰减曲线

通过分析航电测量结果发现, 这种富含铅锌的的多金属矿具有高导电性的特点, 且矿体与围岩电阻率差异较大。这种特征可用于寻找周边类似特征的未知矿床。

2.3 地下强低阻体的时间常数特征

康古尔金矿位于塔里木板块北部被动大陆边缘阿奇山— 雅满苏裂谷带之北缘。矿区出露地层主要为下石炭统雅满苏组火山碎屑岩。矿区内有3条EW— NEE向韧— 脆性剪切带, 剪切带内安山岩、安山质凝灰角砾岩和英安质凝灰角砾岩均发生了糜棱岩化, 矿体就赋存在这些剪切带中。矿石中主要金属矿物为黄铁矿、磁铁矿、黄铜矿、方铅矿、闪锌矿、磁黄铁矿、赤铁矿。成矿作用分两个阶段, 即富硫化物含金石英脉阶段(主成矿阶段)和石英碳酸盐脉或碳酸盐脉阶段(晚成矿阶段)[13]。该矿床是一个受韧性挤压带控制的晚华力西期中低温岩浆变质热液型多金属矿床[14]

在该矿区上方完成测线间距为250 m的航电测量。图9为根据z分量dB/dt计算的时间常数τ 等值平面。矿区范围内的τ 高值区, τ 值一般大于0.72ms, 最高可达4.08 ms, 推测地下为强低阻体。τ 高值区呈团块状, 周围被一圈τ 低值区包围。该区为连续的高导电区, 位于一近EW向含矿蚀变带上, 推测四周快速变弱的τ 值区为矿的尖灭区或采空区。

图9 康古尔金矿时间常数等值线平面

图10为选取的L2线的综合剖面, 可以看出该剖面dB/dt曲线在矿区上方出现了不同异常响应的波形。同时, 在电阻率深度成像图上的相应位置出现一个视电阻率值< 20 Ω · m的低阻异常区。该异常区宽约200 m。图中左侧为矿区外的低阻异常, 预测与地质构造有关, 这里不作解释。

图10 康古尔金矿区L2线综合剖面示意

图11为该异常上的dB/dt衰减曲线。可以看出, 这种地下强低阻体的异常具有早期道幅值相对较小、衰减特别慢的特征。

图11 L2线强低阻异常上dB/dt衰减曲线

该区电力干扰不大, 对测量结果影响不大。

由上可知, 这种含金富硫化物具有高导电性的特点, 且矿体与围岩电阻率差异较大; 在规模上, 低阻区宽度不是很大。这种特征可用于寻找周边类似特征的未知矿床。

2.4 表层强低阻体的时间常数特征

新疆东部是中国主要的盐类矿床聚集区, 在该区域内开展航电测量, 发现石盐矿上方有明显的异常响应。图12为根据z分量dB/dt计算的时间常数τ 等值线平面。图中蓝色线范围内出露地层为第四系全新统, 主要有石盐、盐渍化亚砂土, 该区内有一高亮区为τ 高值区, τ 值一般大于0.77 ms, 最高可达1.68 ms, 推测为表层强低阻体, 由石盐矿引起。该异常区呈团块状, 由于处在工作区边界, 因此异常未闭合。

图12 新疆东部某石盐矿时间常数等值线平面

L3穿过该石盐矿区, 由图13可以看出该剖面dB/dt曲线在石盐矿上方出现幅度较大的单峰异常, 同时, 在电阻率深度成像图上的相应位置出现一个视电阻率值< 20 Ω · m的低阻异常区。该异常区宽约1 km。

图13 新疆东部某石盐L3线综合剖面

图14为L3线中部石盐矿响应衰减曲线和L4线南部可能矿致异常的衰减曲线, 可以看出表层强低阻的石盐矿衰减曲线几乎为一条直线且早期幅值特别高, 可能的矿致异常上早期幅值较小, 并呈指数衰减。该区无电力干扰。由航电成果资料分析可知这种石盐矿具有高导电性, 异常响应幅度较大, 相对一般的多金属矿产, 其规模较大。这种特征可用于寻找周边类似特征的石盐或卤水区。

图14 石盐矿及可能矿致异常上dB/dt衰减曲线

3 结论

时间常数是一个综合参数, 它覆盖时间域和所有导电区, 对许多由不同导电目标源组成的复杂导电区域具有优良的分辨率, 与良导目标体显示信号相应的深度和水平响应会被各自放大和增强; 因此, 时间常数常常会被用来圈定靶区。

地下低阻地质体的航电异常响应具有早期幅值不高、衰减缓慢的特征, 而时间常数表现为高值区。

本次对电力干扰、表层强低阻体、表层弱低阻体和地下低阻体的dB/dt航电异常响应及时间常数τ 的特征进行了综合分析, 可供从事航电资料解释的同行参考。

The authors have declared that no competing interests exist.

参考文献
[1] 梁盛军, 张力卡, 曹学峰, . 时间域航空电磁法研究进展综述[J]. 地质与勘探, 2014, 50(4): 735-740. [本文引用:1]
[2] Annan P. Airborne EM-Lessons learnt and future opportunities[C]//13rd SAGA Biennial Conference and Exhibition, and 6th International AEM Conference and Exhibition, Conference proceedings, 2013: 36-37. [本文引用:1]
[3] Vallée M A, Smith R S, Keating P. Metalliferous mining geophysics- State of the art after a decade in the new millennium[J]. Geophysics, 2011, 76(4): W31-W50. [本文引用:2]
[4] 骆燕, 曾阳, 石岩, . 航空电磁法在火山岩型块状硫化物矿区的试验[J]. 物探与化探, 2014, 38(4): 840-845. [本文引用:1]
[5] 李怀渊, 张景训, 江民忠, . 航空时间域电磁法测量系统VTEM应用效果分析[J]. 地质论评, 2015, 61(s): 665-666. [本文引用:1]
[6] 裴易峰, 殷长春, 刘云鹤, . 时间域航空电磁磁场计算与应用[J]. 地球物理学进展, 2014, 29(5): 2191-2196. [本文引用:2]
[7] Prikhodko A, Geo P. EM time constant(tau) analysis[M]. 2010. [本文引用:2]
[8] 刘冲, 王宇航, 皇健. 瞬变电磁视时间常数tau成像分析与应用研究[J]. 物探化探计算技术, 2014, 36(1): 28-34. [本文引用:1]
[9] 王亚冉. 时间域航空电磁数据的Tau域分析与研究[D]. 长春: 吉林大学, 2014. [本文引用:1]
[10] 静恩杰, 李志聃. 瞬变电磁法基本理论[J]. 中国煤田地质, 1995, 7(2): 83-87. [本文引用:1]
[11] 内蒙古地质工程有限责任公司, 核工业航测遥感中心. 内蒙古翁牛特旗南敖包山-大朝阳沟一带浅覆盖区铅锌多金属矿找矿方法实验及验证续作项目评估报告[R]. 内蒙古自治区地质勘查基金管理中心, 2013. [本文引用:1]
[12] Luo Y, Zhu L, Zhao Z, et al. Application prospects of VTEM helicopter time-domain electromagnetics in China[C]//International Workshop on Gravity, Electrical & Magnetic Methods and their Applications, Chengdu, China, 2015: 285-288. [本文引用:1]
[13] 王志良, 毛景文, 吴淦国, . 东天山康古尔金矿成矿晚阶段地幔流体参与成矿作用的碳氢氧同位素证据[J]. 地质学报, 2004, 78(2): 195-202. [本文引用:1]
[14] 马天林, 孙立, 徐兴旺. 新疆东天山康古尔金矿控矿构造特征[J]. 地质力学学报, 1998, 4(2): 45-52. [本文引用:1]