航空大地电磁法在辽宁省丹东地区的应用
The application of aeromagnetotelluric survey technology to Dandong area, Liaoning Province
责任编辑: 沈效群
收稿日期: 2019-10-21 修回日期: 2020-01-23 网络出版日期: 2020-08-20
基金资助: |
|
Received: 2019-10-21 Revised: 2020-01-23 Online: 2020-08-20
作者简介 About authors
王志宏(1973-),男,汉族,正高级工程师,毕业于中国地质大学(北京)地球物理工程专业,现主要从事航空电磁法生产研究工作。Email:
本文简要介绍了航空大地电磁法原理、测量系统、倾子总散度、总相位旋转参数及二维、三维反演等;利用正演程序,计算了横向电性异常体倾子响应,结果表明倾子资料对横向电性分界面反映明显;最后结合实际地质情况,利用实测倾子及反演电阻率资料,推断解释了区内岩体、控矿构造的分布特征。实际应用表明,航空大地电磁法在探测断裂、岩体等方面具有较好的效果。
关键词:
This paper briefly introduces the principle of aeromagnetotelluric method, the measurement system, the calculation of tiltepper parameters and the two-imensional and three-dimensional inversion. By using the forward program, the tipplter parameters of the geoelectric model of the transverse electrical interface were calculated. The results show that the tilter data obviously reflect the transverse electrical interface. Finally, combined with the actual geological situation, the distribution characteristics of rock mass and ore-control structure in the area were found by using the measured tilter and inversion resistivity data, and some results were obtained, which verify the feasibility of the method.
Keywords:
本文引用格式
王志宏, 江民忠, 彭莉红, 程莎莎.
WANG Zhi-Hong, JIANG Min-Zhong, PENG Li-Hong, CHENG Sha-Sha.
0 引言
航空大地电磁法是利用全球范围内雷电活动引起的电磁场信号,依据地下介质的导磁性和导电性的差异,通过飞行器来测量不同频率的电磁场的频率特性和空间分布,由此来获取地下地电结构信息的一种物探方法。该方法垂直磁场分量在空中测量,水平磁场分量则在地面某一基站观测,通过远参考技术来降低垂直与水平磁场之间的噪声相关性,通过倾子将磁场的垂直分量与水平分量联系起来[1]。该方法具有工作效率高、勘探深度大的特点,可以快速获得地下三维电性结构特征信息。
目前国际上一些学者利用航空大地电磁法探测深部的地质构造和矿体,并且已取得了较为丰富的成果[24]。国内学者关于倾子的研究大都集中在大地电磁测量方法,航空大地电磁倾子方面相对较少[59],且多侧重于理论、模拟数据研究等,较少涉及应用。许智博等研究了倾子资料起伏地形的二维正反演研究[10],李志强等研究了倾子数据的三维正反演[11],张铭等通过对航空大地电磁法倾子响应特征影响因素分析,提出了一种利用倾子及其梯度的异常位置和异常幅度对地下二维异常体进行快速、定量识别的方法[12]。核工业航测遥感中心近年来利用航空大地电磁法开展了大量的生产、研究工作[13],本文是在大量生产研究的基础上,对方法原理、测量系统、数据处理及应用等方面进行了探索与研究,为国内航空电磁法研究提供技术支撑,进而填充国内航空电磁法技术的空白。
1 航空大地电磁法原理
通常在三维条件下,电场的水平分量与磁场的水平分量之间存在着复系数线性关系,它们直接通过阻抗张量联系起来;而磁场的垂直分量和水平分量之间也存在复系数线性关系,二者通过倾子矢量联系起来[14]。
一般而言,大地电磁倾子表示地表同一测点处垂直磁场分量与两个水平磁场分量之间的关系,而航空大地电磁法倾子则表示地面以上空中同一平面不同位置的垂直磁场与地表某一基站处的两个水平磁场的关系,其关系表达式如下:
解方程组可得到倾子的计算表达式;
2 倾子响应模拟
倾子作为航空大地电磁法的主要参数,其能够作为复杂地质构造的表征,不依赖于绝对电导,对电阻率横向不均匀性反应非常灵敏,特别是在探测垂直或者倾斜异常体(蚀变带、断层带等)的空间分布情况时作用非常明显,在地质解释中起到十分重要的作用。在均匀大地中建立如图1所示的地电模型,其背景电阻率为50 Ω·m,模型长度2 000 m,测点41个(点距50 m),平距500 m处存在一250 Ω·m的高阻异常,埋深为400 m,厚度为150 m;平距 1 350 m处存在一25 Ω·m的低异常,异常体宽为100 m,下延深度大于1 500 m,倾向小号点。
图1
图2
3 应用实例
3.1 地质概况
图3
图3
工作区地质简图
1—第四系;2—侏罗-白垩系未分小岭组;3—盖县组;4—大石桥组;5—高家峪组;6—里尔峪组;7—白垩纪二长花岗岩;8—侏罗纪二长花岗岩;9—三叠纪二长花岗岩;10—古元古代花岗岩;11—花岗斑岩;12—实测正断层;13—实测逆断层;14—实测性质不明断裂;15—铅锌矿;16—金矿;17—航空大地电磁测线
Fig.3
Simple geology over survey of working area
1—Quaternary;2—Jurassic-Certaceous Xiaoling Group;3—Daixian Group;4—Dashiqiao Group;5—Gaojiayu Group;6—Lieryu Group;7—Cretaceous monzogranite;8—Jurassic monzogranite;9—Triassic monzogranite;10—Paleoproterozoic granite;11—granite porphyry;12—measured normal fault;13—measured reverse fault;14—unknown fault measured;15—lead zinc ore;16—gold mine;17—aeromagnetotelluric line
3.2 数据采集
图4
图5
航空大地电磁法测线布置为SN向,共布置60条,测线间距400 m。数据采集时直升机沿地形飞行,空中接收线圈平均离地高度97.2 m,飞行速度平均为95.04 km。对采集的航空大地电磁数据进行去噪、补偿及傅里叶变换后获得25、37、75、150、300、600 Hz六个频率沿测线方向、垂直测线方向倾子的实部、虚部数据,航空大地电磁数据采集同时获得高精度航磁数据。
3.3 定性解释
图6为航磁ΔT及航空大地电磁150 Hz总散度平面,箭头所指为断裂或岩体位置。可以看出,两种资料中断裂构造具有明显的线性磁异常、总散度线性高值异常等,如NE向二道沟断裂,NW向青城子、尖山子断裂,白云金矿北侧断裂等。中生代双顶沟、洼岭岩体在两种资料中也有明显的特征,古元古代方家隈子、大顶子岩体无明显的航磁异常,然而在航空大地电磁测量资料中则具有明显电磁异常,一般表现为岩体内部为低值特征,岩体与围岩接触带存在明显的高值电磁异常环。测量资料证明航空大地电磁、航磁两种资料具有互补性,也进一步反映了航空大地电磁测量对构造、岩体等的探测能力。
图6
图6
航磁ΔT、150 Hz TD等值线平面
a—航磁ΔT;b—150 Hz TD;1—断裂;2—岩体分布范围;3—航空大地电磁测线
Fig.6
Aeromagnetic ΔT and 150 Hz in-phase TD contour
a—ΔT;b—150 Hz TD;1—fracture; 2—rock mass distribution range; 3—aeronautical magnetotelluric survey line
由前文可知,对导体、电阻等电阻率异常体,倾子在电性不均匀体附近均会出现由正到负或由负到正的交叉异常,为了使倾子异常更加明显,对倾子分别进行了总散度(TD)及总相位旋转(TPR)计算。总散度(TD)计算是沿测线和垂直测线方向的倾子(实部或虚部分量)导数之和,计算公式为
对倾子沿测线和垂直测线方向(实部或虚部分量)进行90°相位旋转(PR),获得总相位旋转(TPR),TPR可以弥补总散度资料的不足,保留长波信息,在导体上方显示最大值,计算公式为
图7
图7
倾子实部总散度等值线平面
a—300 Hz;b—75 Hz;图例同
Fig.7
Tipper in-phase TD contour
a—300 Hz;b—75 Hz;The legend is the same as that in
图8
图8
倾子实部总相位旋转等值线平面
a—300 Hz;b—75 Hz;图例同
Fig.8
Tipper in-phase TPR contour
a—300 Hz;b—75 Hz;the legend is the same as that in
方家隈子、大顶子、双顶沟、石家岭、洼岭岩体在TD、TPR资料反映岩体内部为低值异常(高阻),岩体与围岩接触带则为明显的似环状高值异常带(低阻),反映了岩体与围岩接触带附近岩石破碎或蚀变发育等。方家隈子、大顶子岩体随着测量频率的降低,异常特征逐渐减弱,推测两岩体主要分布在浅部;双顶沟岩体随着测量频率的降低异常特征基本不变,推测该岩体由浅至深均有分布;洼岭岩体随着测量频率的降低异常特征有逐渐增强的特征,推测该岩体主要分布于深部。
3.4 反演处理及推断解释
3.4.1 反演处理
二维反演使用AV2DTOPO软件,使用正则化牛顿-高斯非线性算法进行反演,反演数据为沿测线方向(Tzx)的6个频率的倾子实部、虚部,反演时空中接收线圈与地面设置为空气层,数据误差设置为 6.6%,初始模型为半空间均匀模型,背景电阻率设置为750 Ω·m,经过5次迭代后模型收敛,得到反演结果。三维反演使用ZTEM_MT3Dinv软件,反演时把地形和电磁接收线圈离地高度作为反演约束参数,中心区网格尺寸为200 m×200 m×10 m,垂直网格大小随着深度增加几何增长,反演数据为25~600 Hz沿测线及垂直测量线分量的倾子实部、虚部数据,初始电阻率根据全区二维反演结果设为 1 000 Ω·m,电阻率范围为1~1 000 000 Ω·m,相对误差设置为20%,倾子噪声设置为1%。经过20次叠加得到反演结果。
3.4.2 推断解释
图9为区内不同深度的反演电阻率切片,可以看出电阻率切片反映的断裂构造、岩体分布等与航磁、航空大地电磁TD、TPR反映基本一致。青城子断裂、白云金矿北侧断裂随着深度的增加,电阻率异常特征逐渐增加,二道沟断裂、尖山子断裂由浅至深电阻率异常特征基本一致。古元古代方家隈子、大顶子岩体仅分布于浅部,随着深度增加至 1 000 m,其电性特征逐渐消失,而三叠纪双顶沟、新岭、姚家沟岩体由浅至深均反映为高阻,进一步证明了区内的铅锌矿体、金矿体与中生代岩体有密切关系。
图9
图9
反演电阻率深度切片等值线
a—深度200 m;b—深度1 000 m;图例同
Fig.9
Resistivity depth slice contour
a—200 m depth;b—1 000 m depth;the legend is the same as that in
图10
4 结论
正演模拟计算表明倾子资料的实部、虚部、振幅和相位响应能较好地分辨横向不均匀异常体,倾子响应的最大值与横向分界面对应;计算了实测倾子的总散度、总旋转相位等,结合航磁资料及实测倾子数据的2D、3D反演资料,推断解释了区内岩体、控矿构造等的分布特征。实测结果表明,航空大地电磁法在探测岩体、断裂等横向不均体方面具有明显效果。
参考文献
Geophysical exploration with audio frequency natural magnetic fields
[J]. ,
A brief analysis of ZTEM data from the Forrestania test site
[R]. ,
ZTEM airborne AFMAG survey results over low sulphidation epithermal gold-silver vein systems at Gold Springs,south eastern Nevada
[A]// ,
Z-TEM (Airborne AFMAG) tests over unconformity uranium deposits
[A]// ,
埋藏球体的倾子响应特征分析
[J]. ,
Analysis of the tipper response characteristics of buried sphere
[J].
磁倾子矢量的图示分析及其应用研究
[J]. ,
Analysis of tipper visual vectors and its application
[J].
利用有限单元法模拟二维MT倾子响应
[J]. ,
Calculating tipper response in two-dimensional magnetelluric using finite element method
[J].
大地电磁二维倾子和视倾子模拟及其应用研究
[J]. ,
Two-dimensional magnetotelluric tipper and apparent tipper simulation and application
[J].
倾子资料的三维共轭梯度反演研究
[J]. ,
Three-dimensional conjugate inversion of tipper data
[J].
ZTEM起伏地形二维正反演研究
[D]. ,
Research of 2D ZTEM forward modeling and inversion with uneven topography
[D].
ZTEM三维正反演研究
[D]. ,
Research on ZTEM three-dimensional forward modeling and inversion
[D].
基于航空大地电磁倾子特征的二维大地结构探测方法
[J]. ,
Two-dimensional ground structure detection method based on the tipper characteristics of airborne magnetotelluric
[J].
航空和地面天然场电磁法联合开展深部矿产资源勘探
[J]. ,
Deep mineral exploration by airborne and ground natural field electromagnetic method
[J].
二维及三维倾子响应和异常体识别
[D]. ,
2D&3D tipper response and distinguishing to anomalous bodies
[D].
Three-dimensional inversion of ZTEM data
[J]. ,
大地电磁法张量阻抗通用计算公式
[J]. ,
General calculation formula of tensor impedance in magnetotelluric method
[J].
Numerical modeling of ZTEM (airborne AFMAG) responses to guide exploration strategies
[A]// ,
3D finite-difference modeling algorithm and anomaly features of ZTEM
[J]. ,
青城子铅锌金银多金属矿田矿床地质特征及成矿系统分析
[D]. ,
Geological characteristics of the deposits and analysis of the mineralization system of Qingchenzi Pb-Zn-Au-Ag polymetallic ore field
[D].
青城子矿田构造变形结构及其控矿特征
[J]. ,
Tectonic deformation texture and ore-controlling features of Qingchengzi orefield
[J].
/
〈 | 〉 |