2. 中国科学院超导电子学卓越创新中心, 上海 200050;
3. 中国科学院大学, 北京 100049;
4. 吉林大学仪器科学与电气工程学院, 长春 130026
）和带宽大（>200 kHz）等特点, 勘探浅层分辨率高、探测深度大, 优势明显。本文将SQUID系统与商用EM67系统进行了异常环对比实验, 验证了以SQUID代替感应线圈作为接收装置的优势与实用性；将该方法应用于内蒙古四子王旗大井坡北部区域的电磁勘探工作, 获取了高质量探测数据, 反演结果显示浅层和中深层（> 2 000 m）为连续性低阻层；在已有地质资料和物探资料的基础上, 结合反演结果可知, 测线附近浅层500 m存在北东向断裂带, 并推断深部2 000 m处分布着同向断裂带。2. Center of Excellence in Superconducting Electronics, Chinses Academy of Sciences, Shanghai 200050, China;
3. University of Chinese Academy of Sciences, Beijing 100049, China;
4. College of Instrumentation & Electrical Engineering, Jilin University, Changchun 130026, China
) and large bandwidth (>200 kHz). As an important mineral reserve, Siziwangqi in Inner Mongolia is rich in mineral resources. The reliability of SQUID instead of induction coil was verified by comparing the experimental results of the anomalous rings obtained by superconducting receiving system and EM67 TEM instrument. The superconducting TEM method was applied to the northern area of Dajingpo, Siziwangqi, Inner Mongolia. The underground anomalous response information in shallow and deep layers (>2 000 m) was obtained to deduce a continuous low resistivity layer. Combined with the existing geological and geophysical data, the inversion results effectively reveal a north-east trending fault zone at 500 m in the shallow layer near the survey line, and then infer that there is a similar trending fault zone at 2 000 m in the deep layer.0 引言
瞬变电磁法(transient electromagnetic method, TEM)是一种时间域人工源电磁方法, 以探测介质间的导电率和导磁率为物性前提, 主要用于探测低阻目标, 研究地下地电结构。通常情况下, 它通过感应线圈或接地电极接收一次脉冲电磁场激发的二次电磁场信号, 并利用电磁理论建立的相关模型反演出地下介质的电阻率分布特征。
目前国内外市场上性能稳定、可靠的瞬变电磁仪有十几种, 如加拿大GEONIC公司的EM(electromagnetic)系列和吉林大学的ATEM(airborne transient electromagnetic)型瞬变电磁仪等[1-2]。这些TEM仪器均使用感应线圈作为接收装置, 测量磁场B随时间t变化产生的感应电动势A
超导量子干涉器件(superconducting quantum interference device, SQUID)是一种灵敏度极高的磁传感器, 能够测量10-15T量级的磁场波动[7]。同时, SQUID在频带范围内的灵敏度与频率无关, 是理想的瞬变电磁接收传感器[8]。根据工作温度的不同, SQUID分为高温SQUID(工作于77.0 K的液氮环境)和低温SQUID(工作于4.2 K的液氦环境)。由于液氦制备及使用条件复杂, 低温SQUID在TEM中的应用较少；但是低温SQUID噪声(5~7 fT/
本文对比了异常环实验中感应线圈与低温SQUID的衰减曲线, 并在内蒙古四子王旗大井坡北部测区内应用低温超导瞬变电磁系统进行探测实验, 完成2条测线的勘探, 反演深度大于2 000 m。本文揭示了测线附近的地质分布情况, 验证了浅层断裂带的位置和走向, 并推断出断裂带在深部的分布, 以期为评估该区域矿产资源潜力提供参考。
1 四子王旗大井坡地球物理特征本文选取我国首个航空地球物理综合实验区(大井坡实验区)的局部地区进行精细探测。从整体上看, 该实验区海拔从北至南逐渐升高, 依据地质资料可将全区由北西至南东划分为3个地质区(图 1)[9]。实验区内地层主要有中新元古界、中生界和新生界出露。主要地层有:中、新元古界为一套浅海碎屑岩-碳酸盐岩建造, 在Ⅱ区中部大面积出露；中生界除少量玄武岩、凝灰岩等火山碎屑岩外, 以陆相碎屑岩沉积为主, 主要分布于Ⅲ区东南部的大青山组；新生界在Ⅰ区内分布最广, 为一套河流相碎屑岩建造, 岩性以砂岩、泥岩及现代冲洪积物为主。Ⅱ区与Ⅲ区之间有一条天然沟带, 带内变质闪长岩及含阿牙登岩组的结晶灰岩夹板岩含量丰富[10], 这些岩石表现为低电阻率特性。同时, Ⅱ、Ⅲ区之间的天然沟带从测区中部穿过, 实地考察发现区域内铁矿石出露明显。
|
| 1.全新统湖积洪积砂砾层；2.上新统宝格达乌拉组泥岩、砂质泥岩；3.上侏罗统大青山组砾岩、砂岩、粉砂岩；4.新元古界呼吉尔图岩组石英岩、变质石英砂岩、变质粉砂岩、千枚岩、变质砂岩、微晶绿帘角闪岩及结晶灰岩；5.新元古界白音宝拉格岩组变质石英砂岩、石英岩；6.中元古界比鲁特岩组板岩、斑点板岩；7.中元古界哈拉霍圪特岩组灰岩, 夹板岩、变质含砾石英砂岩；8.中元古界尖山组灰岩-锈红色变质石英砂岩、板岩；9.中二叠世花岗岩；10.早泥盆世灰黑色辉长岩；11.时代不明超基性岩；12.花岗斑岩脉；13.花岗伟晶岩脉；14.闪长岩脉；15.石英脉；16.闪长玢岩脉；17.断裂；18.动力变质带；19.典型剖面名称及位置；20.测线1；21.测线2；22.根据地质图所划分的区域；23.测区位置；24.呼勒斯太以北磁铁矿；25.小南山含铂铜镍矿；26.疑似断裂。据文献[9]修改。 图 1 四子王旗大井坡地质分布及测区位置 Fig. 1 Geological distribution map of Dajingpo, Siziwangqi and the position of measuring areas |
|
|
北京勘察技术工程有限公司应用瞬变电磁方法在大井坡实验区进行了全范围的地质勘探[9], 图 2所示为t=1.417 1 ms的视电阻率平面图。从图 2可知, Ⅰ区到Ⅲ区的视电阻率依次为低阻、高阻、次高阻, 且Ⅱ区和Ⅲ区之间有一条连续的低阻带。根据工区地质情况, 此条低阻带与Ⅱ、Ⅲ区的天然沟带在同一地理位置, 且沟带内的地质特性以低电阻率的岩石为主。对照分析可知, 此条近东西向展布的低阻带处, 应为断裂破碎带的反映。
已有地质资料显示[11-12]:大井坡高阻区广泛分布中二叠世花岗岩和花岗伟晶岩脉, 电阻率大于3 000 Ω·m；低阻区分布全新统湖积洪积砂砾层和上新统宝格达乌拉组, 电阻率小于200 Ω·m。由上述资料可知, 该地区电阻率差异明显, 为电阻率法勘探提供了良好的物性前提。
2 超导瞬变电磁系统超导瞬变电磁接收系统与发射机通过GPS实时同步, 发射机工作期间, SQUID响应磁场变化, 并利用读出电路传输至接收机, 如图 3所示。与感应线圈相比, 超导接收系统具有更低的噪声水平和更大的带宽, 具体参数如表 1所示。
|
| 图 3 超导瞬变电磁接收系统 Fig. 3 Superconducting transient electromagnetic receiving system |
|
|
| 噪声水平/(fT/) |
转折频率/Hz | 系统带宽/kHz | 传感器带宽/MHz | 摆率/(mT/s) | 质量/kg | 工作时长/d |
| 7 | 10 | 200 | 2 | 3 | ＜15 | 3 |
在数据处理中, 通常采用正负叠加技术抑制工频噪声。由于SQUID测量磁场信号B, 系统还易受低频噪声(尤其是地球磁场波动)干扰。考虑到发射信号发射时间为80 ms, 关断时间为920 ms, 地磁场在发射信号的一个周期内可近似于线性变化, 可以利用线性插值分段拟合的方式去除低频噪声[13]:首先对发射时间晚期道噪声段进行线性拟合, 获得短时间内地球磁场的变化规律, 然后把该规律应用于整个关断阶段以去除地磁场的干扰。
地磁场校正前后的叠加结果如图 4所示。从图 4可以看出:由于地磁场的干扰, 校正前的原始数据在400 ms处偏离了原先的衰减趋势, 急剧衰减到零值以下；经过校正后叠加信号的晚期信号衰减时间延长至700 ms, 此前仍保持原先的衰减趋势。叠加完成后, 对数据进行抽道, 再通过全区视电阻率的定义得到该测点视电阻率与深度的关系, 最后将结果作为初始模型进行反演, 数据处理流程如图 5所示。
|
| 图 4 二次场信号校正前后对比 Fig. 4 Comparisons of secondary field signal before and after correction |
|
|
|
| 图 5 模型建立及数据处理流程图 Fig. 5 Flow chart of model building and data processing |
|
|
本文采用奇异值分解(singular value decomposition, SVD)方法直接对磁场进行反演。首先建立目标函数:
(1) 式中:E为拟合误差；N为数据总量；di为实测数据；mk为k次迭代后的模型参数；F(mk)为瞬变电磁法一维正演方程。
根据奇异值分解法, 可将式(1)中di-F(mk)转换为
(2) 式中: 
为验证SQUID系统的可靠性, 本文选择传统的商业化仪器EM67(接收装置为感应线圈)作对比。发射线框为100 m×100 m, 发射电流16 A, 基频25 Hz, 分别用SQUID与EM67接收机进行异常环实验。将SQUID采集的B信号做数值微分后与EM67所测
|
| 图 6 SQUID与感应线圈异常环实验数据对比 Fig. 6 Comparison of data acquired by the SQUID and induction coil |
|
|
结合图 1和图 2可知, Ⅱ区和Ⅲ区之间存在浅层低阻断裂带, 该断裂带具有一定的成矿条件；根据勘探需要, 选择对该断裂带应用超导瞬变电磁系统进行精细探测。
实验选择垂直于该低阻带的方式布置2条测线(图 1)。测线1和测线2各有6个测点, 2条测线间距200 m, 且每个测点间隔均为200 m。布置200 m×200 m线框, 发射电流80 A, 发射信号基频0.5 Hz, 其中发射时间为80 ms, 关断时间为920 ms, 每组采集时间约为9 min, 叠加次数500次。图 7为测线1、2的采集信号抽道曲线剖面图。从图 7中可以看到, 测线1、2均整条测线上数据一致性良好。
|
| a—d分别表示从早期到晚期(0.07 ~ 6.27 ms)不同时刻的抽道数据。 图 7 测线1(左)、2(右)感应磁场剖面曲线 Fig. 7 Induction magnetic field profile curve of Line 1 (left) and Line 2 (right) |
|
|
测线1和测线2中各个测点反演结果与实测数据拟合误差如表 2所示。各点拟合误差均小于10%, 拟合效果较好。
| % | |||||||||||||||||||||||||||||
| 测点1 | 测点2 | 测点3 | 测点4 | 测点5 | 测点6 | ||||||||||||||||||||||||
| 测线1 | 9.21 | 5.58 | 6.26 | 3.37 | 3.28 | 7.94 | |||||||||||||||||||||||
| 测线2 | 9.37 | 2.32 | 6.42 | 2.98 | 4.21 | 8.91 | |||||||||||||||||||||||
结合四子王旗大井坡地质、视电阻率资料及地面初步勘察情况, 可对2条测线的反演结果(图 8)进行综合解释。从整体上分析, 2条测线位于Ⅱ、Ⅲ区交界处, 靠近呼勒斯太磁铁矿床和小南山含铂铜镍矿床(图 1), 地质上为上新统宝格达乌拉组, 周边分布泥岩和砂质泥岩；同时测区内局部出露哈拉霍圪特岩组灰岩、夹板岩, 中二叠世花岗岩, 白音宝拉格岩组石英岩及辉长岩体, 地质情况特殊, 岩石种类丰富；该测区的电阻率变化层次分明, 符合地质资料中低阻断裂带的反映。根据反演结果, 我们对该区域的地质分布情况进行如下推断:
|
| 图 8 测线1和测线2反演结果 Fig. 8 Inversion results of Line 1 and Line 2 |
|
|
1) 测线1与测线2电阻率分布一致, 浅层400 m以内(图 8中1处)均表现为低电阻率特性(电阻率小于200 Ω·m), 与北京勘察技术工程有限公司探测结果(图 2)显示的低阻带分布吻合。推测该低阻特性是由断裂带引起, 且地表发现铁矿石出露, 该低阻断裂带可能属于呼勒斯太磁铁矿床。
2) 2条测线南北两侧电阻率较高, 与图 2中Ⅱ区、Ⅲ区的高阻特性一致。已有地质资料[9]显示, 该地区主要地层为中元古界哈拉霍圪特岩组变质石英砂岩和变泥灰岩互层, 其表现高阻特性, 与测量结果吻合。
3) 2条测线700 m处均显示地下约1 500 m(图 8中虚线2处)呈现较大范围的低阻异常, 水平宽度约200 m, 可能分布深大断裂带。其中测线1小于500 Ω·m的低阻区域延伸到深部2 000 m附近, 且测点之间一致性非常好。2条测线均显示在深度800 ~1 000 m附近, 深部低阻与浅部低阻带存在一定的隔离。根据当地地质资料[10], 已探明测区东部约3 km处为呼勒斯太磁铁矿床, 南部约5 km处为小南山矿床。该测区矿床地质构造复杂, 受近东西向深大断裂和北东东向、北西向断裂控制, 断裂带中古生代晚期辉长岩发育良好。辉长岩是铜镍矿、黄铁矿等金属矿物成矿母岩, 当含有金属元素时, 其表现为低电阻率特性。推测深部低阻异常为东西向深大断裂带中的基底层辉长岩引起, 可能具备一定成矿条件。
5 结论1) 本文对比了EM67瞬变电磁仪与超导瞬变电磁系统的异常环探测结果, 验证了SQUID应用于瞬变电磁法在探测精度和深度上的优势。根据断裂带与周围岩体的电阻率差异, 在内蒙古四子王旗大井坡采用低温SQUID接收和200 m发射线框的瞬变电磁勘探方案, 对该区域电阻率异常进行精准探测, 获得了2 000 m以上的有效反演深度。
2) 本文反演结果显示的低阻带与已有地质资料的断裂带相吻合, 并且进一步推断出深大断裂带分布在地下1 500~2 000 m深度。本文利用超导瞬变电磁法对该区域进行分析, 存在浅部明显断裂和深部疑似断裂带, 为四子王旗大井坡区域的地质构造研究以及矿产勘探提供了参考。
| [1] |
嵇艳鞠, 林君, 程德福, 等. ATEM-Ⅱ瞬变电磁仪数据处理软件的研制与应用[J]. 吉林大学学报(地球科学版), 2003, 33(2): 242-245. Ji Yanju, Lin Jun, Cheng Defu, et al. Development and Application of Data Processing Software of ATEM-ⅡTransient Electromagnetic Instrument[J]. Journal of Jilin University (Earth Science Edition), 2003, 33(2): 242-245. |
| [2] |
嵇艳鞠, 林君, 于生宝, 等. ATEM瞬变电磁系统在长春市活断层勘探中的应用[J]. 吉林大学学报(地球科学版), 2005, 35(增刊1): 103-107. Ji Yanju, Lin Jun, Yu Shengbao, et al. Sounding Application of ATEM Instrument in Changchun Active Fault[J]. Journal of Jilin University (Earth Science Edition), 2005, 35(Sup. 1): 103-107. |
| [3] |
Spies B R, Eggers D E. The Use and Misuse of Apparent Resistivity in Electromagnetic Methods[J]. Geophysics, 1986, 51(7): 1462-1471. DOI:10.1190/1.1442194 |
| [4] |
Clarke J. Geophysical Applications of SQUIDS[J]. IEEE Transactions on Magnetics, 1983, 19(3): 288-294. DOI:10.1109/TMAG.1983.1062524 |
| [5] |
Wang S G, Zhang L H, Wang C J, et al. Application of High Temperature Superconducting Quantum Interference Devices in Transient Electromagnetic Method for Magnetotelluric Soundings[J]. Physica C, 1997, 282/283/284/285/286/287: 411-414. |
| [6] |
陈晓东, 赵毅, 张杰, 等. 高温超导磁强计在瞬变电磁法中的应用研究[J]. 地球物理学报, 2012, 55(2): 702-708. Chen Xiaodong, Zhao Yi, Zhang Jie, et al. The Applications of HTc SQUID Magnetometer to TEM[J]. Chinese Journal of Geophysics, 2012, 55(2): 702-708. |
| [7] |
Clarke J, Braginski A I. The SQUID Handbook[M]. Weinheim: Wiley-Vch, 2004.
|
| [8] |
荣亮亮, 蒋坤, 裴易峰, 等. 适用于瞬变电磁勘探的低温超导磁传感器[J]. 仪器仪表学报, 2016, 37(12): 2671-2677. Rong Liangliang, Jiang Kun, Pei Yifeng, et al. Low-Tc SQUID Sensor for Time Domain Electromagnetic Prospecting[J]. Chinese Journal of Scientific Instrument, 2016, 37(12): 2671-2677. |
| [9] |
张婉, 刘英会, 朱卫平, 等. 大井坡航空地球物理试验场典型剖面综合地球物理解释[J]. 吉林大学学报(地球科学版), 2016, 46(2): 569-580. Zhang Wan, Liu Yinghui, Zhu Weiping, et al. Comprehensive Geophysical Interpretation of Typical Profiles of Aviation Geophysical Test Site at Dajingpo[J]. Journal of Jilin University (Geoscience Edition), 2016, 46(2): 569-580. |
| [10] |
彭芳苹, 董良, 葛志广. 航空物探试验场地面地球物理场特征分析[J]. 工程地球物理学报, 2014, 11(1): 118-122. Peng Fangping, Dong Liang, Ge Zhiguang. The Analysis of Geophysical Characteristics of Airborne Geophysical Testing Ground[J]. Chinese Journal of Engineering Geophysics, 2014, 11(1): 118-122. DOI:10.3969/j.issn.1672-7940.2014.01.022 |
| [11] |
李凤仁, 张赋. 内蒙古四子王旗夏尔楚鲁金矿床地质特征及成因[J]. 地质与资源, 2012, 21(4): 350-353. Li Fengren, Zhang Fu. Geology and Genesis of the Xiaerchulu Gold Deposit in Inner Mongolia[J]. Geology and Resources, 2012, 21(4): 350-353. DOI:10.3969/j.issn.1671-1947.2012.04.002 |
| [12] |
朱卫平, 熊盛青, 刘英会, 等. 大井坡航空地球物理试验场选场原则与方法[J]. 物探与化探, 2013, 37(1): 98-103. Zhu Weiping, Xiong Shengqing, Liu Yinghui, et al. Selection Principles and Methods of Dajingpo Airborme Geophysical Test Site[J]. Geophysical and Geochemical Exploration, 2013, 37(1): 98-103. |
| [13] |
裴易峰, 荣亮亮, 张懿, 等. 低温SQUID瞬变电磁法中去除地磁场影响的方法研究[J]. 地球物理学进展, 2019, 34(2): 212-217. Pei Yifeng, Rong Liangliang, Zhang Yi, et al. Study on the Method of Removing the Effect of Geomagnetic Field in Low Temperature SQUID Transient Electromagnetic Method[J]. Progress in Geophysics, 2019, 34(2): 212-217. |
| [14] |
Ji Y, Du S, Xie L, et al. TEM Measurement in a Low Resistivity Overburden Performed by Using Low Temperature SQUID[J]. Journal of Applied Geophysics, 2016, 135: 243-248. DOI:10.1016/j.jappgeo.2016.09.027 |
| [15] |
Du S, Zhang Y, Pei Y, et al. Study of Transient Electromagnetic Method Measurements Using a Superconducting Quantum Interference Device as B Sensor Receiver in Polarizable Survey Area[J]. Geophysics, 2018, 83(2): E111-E116. |