Multi-Source SATEM 1D Inversion and its Application to the 3D Model Data
-
摘要:
为了避免长导线源地空瞬变电磁装置的体积效应影响，并发挥装置探测深度大工作效率高的优势，分别设计了单个激发场源和多个激发场源的地空瞬变电磁数值模拟，对比分析场源分布对瞬变场的影响以及对地下模型的分辨特征。将一维反演应用于三维地电模型数据的解释之中，讨论复杂的激发源、简单的解释技术实现三维复杂目标反演的可行性。首先，采用三维数值模拟实现复杂激励源地空瞬变电磁三维正演模拟，分析多激励源瞬变场特征，证明可以通过改变源的布设方式减少电性源体积效应影响。然后，利用常规一维反演方法对数据进行解释，从而证明多源发射简单的解释方法也可以提高解释的分辨率。最后，对甘肃省某煤田采空区的野外数据进行一维反演解释，结果表明相较于单辐射源瞬变电磁的反演结果，使用多辐射源瞬变电磁探测方法可以得到更为精准的采空区分布信息。数值模型与实测数据解释结果充分说明复杂激发源即使采用简单的反演方法也能够有效提高解释结果的分辨率，这为提高瞬变电磁解释精度提供了新的思路。
Abstract:In order to avoid the volume effect of long-wire sources in semi-airborne transient electromagnetic device and leverage the advantages of high detection depth and working efficiency, numerical simulations were conducted using single and multiple field sources. The effect of source distribution on the transient fields and resolution characteristics of underground models was analyzed. The feasibility of achieving 3D complex target inversion using simple explanation techniques was discussed by applying the 1D inversion to interpret 3D geoelectric model data. First, the 3D FEM is used to realize the 3D forward modeling of multi-source semi-airborne TEM, analyze the characteristics of the multi-sources transient field, and prove that the volume effect of electrical sources can be reduced by changing the source layout. Then, the 1D inversion method is used to interpret the 3D model data to prove that the simple interpretation method of multi-source device can also improve the resolution of the result. Finally, 1D inversion interpretation of survey data from a coal mine goaf in Gansu Province is carried out. The results show that compared with the results of single-radiation source survey data, the multi-source survey data can be more accurate on the distribution of water zone. The interpretation of the synthetic model and the survey data demonstrate that the resolution of the results can be effectively improved even if simple inversion methods are used for complex excitation sources, which provides new ideas and useful explorations for improving the accuracy of TEM interpretation.
-
新疆东天山铜镍成矿带是中国重要的铜镍矿集区（Feng et al., 2018; 师震等，2019），分布有香山、黄山、图拉尔根、葫芦等多个铜镍硫化物矿床，其空间分布可划分为两个带，分别受康古尔大断裂和黄山大断裂控制。位于黄山东至图拉尔根成矿带之间的浅覆盖区，是黄山断裂成矿带之间的空白区，成矿潜力巨大。新疆地矿局物化探大队、中国国土资源航空物探遥感中心、中国地质调查局西安地质调查中心等单位在该区域开展过1∶5万综合物化探调查和1∶5万、1∶20万重力、航磁调查以及1∶5万矿产地质调查等。新疆地质矿产勘查开发局第六地质大队根据区域航空磁测成果，在大黄山浅覆盖区发现多个航磁异常，利用综合物探方法推测出类似黄山东的基性-超基性岩体，并在钻孔中发现了含铜镍矿超基性岩体，初步估算矿体规模达到中型。尽管取得了一定的进展，但是该区域依然存在很大的找矿潜力，因此依据大黄山的找矿成果，对大黄山一带浅覆盖区内的物探异常进行筛选，发现红石岗南区域存在1∶5万航磁异常与重力异常套合较好，同时前人研究表明，该浅覆盖区域北侧的红石岗岩体存在形成铜镍硫化物矿床的潜力（王志福等，2012; 王亚磊等，2017; Feng et al., 2018），目前对红石岗南区域调查研究程度不高，该区存在大规模第四系覆盖，无岩体出露，1∶5万航磁数据表明在该区的负异常背景场上叠加少量局部正磁异常，异常形态多呈椭圆状，磁异常幅值较小，因此很难评估该区域存在的隐伏异常体，难以对其成矿潜力进行评价。因此本文在附近区域矿床及岩体研究的基础上，开展高精度地面磁测，对该调查区磁异常数据采用磁异常化极、垂向一阶导数、归一化总梯度法、二维剖面模拟和三维磁异常反演方法进行研究，实现对区内隐伏异常体和断裂构造的初步勘探研究，确定了异常体的数量、磁化强度、埋深、产状要素、规模及断裂构造的位置，为探测隐伏异常体和后续在该区域所开展的勘探工作，提供了可信的地球物理依据和勘查实例。
1. 区域地质和地球物理背景
东天山红石岗南调查区位于新疆维吾尔自治区哈密市东南约160 km，整体上大地构造位于准噶尔洋壳板块与塔里木陆壳板块的聚合处，中天山隆起和吐哈盆地之间（王亚磊等，2017），区内地层出露主要以泥盆系和石炭系为主，少量奥陶系、志留系和侏罗系（Xiao & Wang, 2017）。红石岗南调查区在康古尔－黄山韧性剪切带内（宫辰，2020），处于东天山铜镍成矿带东段，附近存在多处铜镍矿床，如图拉尔根、葫芦、黄山、黄山东、香山铜镍矿等，红石岗镁铁-超镁铁质岩体位于土墩-镜儿泉断裂和黄山-镜儿泉断裂之间，而红石岗南调查区位于黄山-镜儿泉断裂以南，整体处于镜儿泉矿集区（图1a）。在黄山-康古尔断裂的北部，主要岩层为梧桐窝子组深灰色-灰绿色的海相喷发岩，而南侧则出露干墩组深灰-灰黑色浅变质的硅质岩、泥质岩等。区内铜镍矿和镁铁-超镁铁质岩体主要分布于这两套地层中（王志福等，2012）。调查区内侵入岩以花岗岩类为主，镁铁-超镁铁质岩为辅，镁铁-超镁铁质岩主要为晚石炭世-早二叠世（王亚磊等，2017），铜镍矿床也主要赋存于这些岩体中，整体呈南西-北东向分布，年代学研究表明矿集区内铜镍硫化物矿床均形成于早二叠世（韩宝福等，2004; Mao et al., 2016; 师震等，2019）。
![]()
图 1 （a） 东天山地区主要铜镍矿床及镁铁-超镁铁岩体分布图（据Mao et al., 2014修改）；（b） 东天山地区航磁异常图；（c） 红石岗南地区航磁异常图；（d） 红石岗南地区布格重力异常图Figure 1. (a) The distribution map of important copper-nickel deposits and mafic-ultramafic intrusion in the Eastern Tianshan Mountains ; (b) Aeromagnetic anomaly map of the Eastern Tianshan Mountains; (c) Aeromagnetic anomaly map of the southern Hongshigang area; (d) Bouguer gravity anomaly map of the southern Hongshigang area区域航磁调查成果显示（图1b）黄山-镜儿泉断裂以北存在诸多正磁异常，异常幅值大于200 nT，其中位于黄山和香山附近，夹于土墩-镜儿泉断裂和黄山-镜儿泉断裂的区域存在一片高磁区域，异常最强处幅值高于800 nT。该区域存在的矿床普遍位于土墩-镜儿泉断裂和黄山-镜儿泉断裂附近，矿床附近区域磁异常整体小于200 nT。雅满苏断裂形成连续的串珠状磁异常（图1b），而雅满苏断裂以北，大范围新生代沉积物覆盖区，磁异常很低。南侧中北天山北缘断裂东北部存在条带状高磁异常带，异常幅值高于200 nT。红石岗南调查区附近磁异常相对香山和黄山区域整体较弱（图1b~图1c），磁异常呈扁椭圆状，长轴方向约40°，长约8 km，宽约4 km，异常值140~250 nT，极大值282 nT，北侧伴生似圆形正磁异常，长约3.5 km，异常极大值约252 nT（图1c）。整个调查区对应在高重力和低重力异常之间的梯度带上，等值线呈密集束状，梯度变化率4.4 mGal/km（图1d）。
整个红石岗南调查区表层覆盖新生代沉积物（图1a），沉积层较厚，地形平缓，受控于区域性主干断裂带所派生的次级断裂或裂隙系统，可能为镁铁-超镁铁质岩浆的就位与成矿提供了良好的导矿和储矿空间（王亚磊等，2017）。由于岩浆活动是矿产形成的关键性因素之一（杨大欢等，2022; 肖丹等，2022），因此红石岗南调查区可能会存在铜镍硫化物矿床赋存于镁铁-超镁铁质岩体中。
2. 岩矿石磁性特征
由于调查区内覆盖层较厚，且针对该区域还未开展任何钻探工作，因此需要使用调查区附近岩石样品的磁性数据参考，根据前人的研究，区内岩石出露的主要岩石类型为闪长岩、辉长岩和苏长岩等（王亚磊等，2017; 师震等，2019），因此统计了附近铜镍矿床岩石标本的体积磁化率、剩磁和Q比（图2，表1），其中Q比是柯尼希斯贝格比（Koenigsberger ratio, Q比）（Koenigsberger, 1938），是剩余磁化和感应磁化的比值，$ Q={M}_{R}/{M}_{i}=\mathrm{N}\mathrm{R}\mathrm{M}[A/m]/\kappa \left[SI\right]\times H[A/m] $，其中$ {M}_{R} $表示剩余磁化，$ {M}_{i} $表示感应磁化，NRM是自然剩余磁化强度（Natural remanent magnetization），$ \kappa $是岩石样品的磁化率，H是当地地磁场强度，在文中计算时使用的地磁场强度为45.68 A/m （
57406.2 nT）。Q比大于1表明磁异常场主要是由剩余磁化引起的，Q比小于1主要是由感应磁化引起的，Q比等于1表明感应磁化和剩余磁化起同样的作用。一般而言大陆喷出岩具有较高的Q比（Parkinson & Barnes, 1985），但具有商业价值的铁矿床几乎都具有低Q比（Jahren, 1965）。虽然这些统计数据并不能代表所有岩性单元的真实值，但根据区内标本的实测值可以为磁异常的正反演提供有效的约束。![]()
表 1 区域岩矿石磁性范围统计表Table 1. Statistical table of magnetic range of rock and ore岩石类型 磁化率（10−5 SI） 剩磁（10−3 A/m） Q比 备注 区域变质岩 14-72 13-25 0.4-3.91 包括变粒岩、片岩等 沉积岩 50-101 50-100 1.10-4.44 包括砂岩、粉砂岩 花岗岩 19-50 10-50 1.15-2.22 图拉尔根、黄山东样品 闪长岩 12-458 11-150 0.17-2.22 图拉尔根、图拉尔根钻孔、黄山东样品 凝灰岩 40-290 11-50 0.15-2.22 图拉尔根、图拉尔根钻孔、黄山东样品 辉长岩 91-165 14-150 0.29-2.02 图拉尔根钻孔、黄山东样品 辉石岩 37-806 11-226 0.29-0.65 图拉尔根钻孔样品 橄榄岩 2160 -3610 196-792 0.20-0.45 黄山东、图拉尔根岩心、土墩井中、黄山井中样品 矿化橄榄岩 7470 -17800 1290 -12200 0.13-1.91 黄山东、土墩井中、黄山井中样品 铜镍矿石 3560 -18800 417- 7910 0.23-1.88 图拉尔根岩心样品 统计结果表明（图2），该区域的岩矿石磁化率的常见值处于101× 10−5~105 × 10−5 SI之间，剩磁常见值则为101× 10−3~105 × 10−3 A/m。橄榄岩、矿化岩石和铜镍矿石普遍具有高磁化率和高剩磁，与围岩存在显著的磁性差异（图2），和矿化有关的这些岩石磁化率高于103 × 10−5 SI，绝大多数样品的Q比小于1，这表明在矿化区域，岩矿石的感应磁化可能主导着区域磁异常。而调查区附近的区域变质岩、沉积岩和花岗岩的磁性偏弱（图2和表1），闪长岩、凝灰岩、辉长岩和辉石岩稍强，但磁化率整体小于103 × 10−5 SI，剩磁强度小于300 × 10−3 A/m。该区域所有岩矿石的Q比处于0.1~10之间，大部分样品的Q比小于1（图2），这表明该区域的围岩磁性可能也是由感应磁化主导的，当然Q比也与岩矿石中的磁性矿物的浓度和粒度有关（刘隆等，2021），需要进一步研究。因此，根据图2和表1的结果，含矿超基性岩体的磁化率普遍大于
6000 × 10−5 SI，因此在后续二维剖面模拟和三维反演中，假设异常体平均磁化率为6000 × 10−5 SI进行物性约束。3. 数据采集与处理
2022年由中国地质调查局西安地质调查中心在红石岗南调查区（图1）进行了100 m × 40 m网格的地面高精度磁测，测区约80 km2。数据采集和处理流程见图3，在获得地磁场总强度T后，通过日变改正和基于国际地磁参考场IGRF13th模型（Wardinski et al., 2020）的正常场改正（IGRF改正）之后得到地磁场异常$ \mathrm{\Delta }T $，整个调查区的研究区总磁场强度、磁倾角和磁偏角的平均值分别为57 335.79 nT、63.56°和−0.35°，磁测总精度为2.44 nT。
对磁异常数据进行化极处理，以消除斜磁化影响（Baranov, 1957; 骆遥，2013），并进行垂向一阶导数计算，更好地确定磁性体边界。而为了确定异常体的场源信息，对化极后的磁异常数据剖面进行归一化总梯度法计算。归一化总梯度法通过在跨越磁场源的截面上构造一个特殊变换的场（总归一化梯度），使其可以用来探测源位置（Elysseieva, 2019）。另一方面，为了估算到异常体的深度并圈定异常体形态，使用了Geosoft公司的Oasis Montaj软件中的2D GM-SYS建模模块，创建一个假设的地质模型，并计算磁响应（Xiangjin et al., 2017; Ekwok et al., 2019, 2022）。用于建模的两条剖面A’-A和B’-B数据来源于原始磁异常剖面数据，剖面走向分别为315°和334°，模拟中使用真实地形数据，观测高度为实际仪器高度，约2 m。
最后进行三维反演和可视化，在本研究中为了解决正则化平滑反演难以揭示真实的地质情况，模糊的边界不能很好地反映地质结构（Vanzon, 2006; Sun , 2014; Utsugi, 2019）和平滑度反演通常会低估恢复的磁化率的值（Sun, 2015）的情况，所以采取稀疏范式的三维反演方法（Meng, 2018）。在正则化中使用混合$ {\mathit{L}}_{p} $范数解决反演问题。利用迭代重加权最小二乘法（Iterative Reweight Least Square, IRLS）进行离散和评估，改写最小的模型组分，通过有限差分算子代替梯度项等步骤，得到最终的正则化函数。$ {\mathit{L}}_{p} $范数正则化是高度非线性的，随着迭代过程的进行和阈值趋近于0，正则化函数的重点关注模型值的范围，直到迭代持续到算法达到预定义的收敛标准。开源框架SimPEG（Cockett et al., 2015）已实现混合$ {\mathit{L}}_{p} $范数反演，本文中的三维模型反演工作均在SimPEG中进行，反演中使用地形数据为真实地形数据。
通过对磁异常数据的处理、二维剖面模拟、三维反演及可视化，可以确定异常体的数量、磁化强度、埋深、产状要素等。进一步结合区域地质资料和物性资料，判断异常体的岩性、规模和形态等（图3）。
4. 结果
4.1 平面磁异常解释
原始磁异常等值线图表现磁异常带总体呈现近SW—NE走向，异常中心存在高值，可达220 nT，异常高值呈条带状（图4a），主要异常位于调查区南侧，呈现近SW—NE走向。化极后磁异常形态突出，圈定了5个主要的磁异常（图4b）。C-1异常位于红石岗南调查区西部，整体异常幅值较小，幅值约210 nT，呈现近SW—NE走向，为条带状正磁异常（图4b）。因此该处异常体为近SW—NE走向，北侧等值线下降较缓慢，因此异常体埋藏可能较深并北倾。C-2异常是整个调查区的主要异常，位于调查区中心，呈现近SW—NE走向，椭圆状正磁异常，异常幅值约为160 nT（图4b），因此异常体为近SW—NE走向。通过对等值线异常分析，其正值范围很大，西北侧等值线下降缓慢，而东南侧下降很快，表明异常体向西北侧倾斜且埋藏很深。C-3异常位于调查区北部，为圆状正磁异常，幅值约160 nT（图4b）。C-4异常则处于调查区东北部，异常整体走向为近SW—NE走向，异常幅值约110 nT，东北侧等值线下降缓慢，而西南侧下降很快并出现较低值，表明异常体向东北侧倾斜且埋藏很深（图4b）。C-5异常幅值与C-4异常接近，但异常等值线西侧等值线下降缓慢，而东北侧下降很快并出现较低值，所以异常体可能向西侧倾斜（图4b）。同时化极磁异常显示在调查区北侧和南侧存在串珠状磁异常，与存在的断裂构造相关，其中北侧异常走向为南西-北东（图4b-c），垂向一阶导数更清晰地展示了该处异常（图4c蓝色箭头所指推测断层），可能是黄山-镜儿泉断裂。南部异常（图4c紫色箭头所指推测断层）走向一致，可能为区内次级断裂。垂向一阶导数结果也显示出其他5个异常体的顶面局部隆起处的地质边界（图4c）。
![]()
图 4 （a）调查区原始磁异常平面等值线图；（b）调查区磁异常化极平面等值线图；（c）调查区化极磁异常垂向一阶导数；（d）剖面A’-A和B’-B的磁异常归一化总梯度计算结果Figure 4. (a) Contour map of the RTP magnetic anomaly in the survey area; (b) Contour map of the RTP magnetic anomaly in the survey area; (c) The first vertical derivative of the RTP magnetic anomaly in the survey area; (d) Normalized total gradient of magnetic anomalies for profiles A’-A and B’-B选取的红石岗南调查区磁异常数据两条剖面A’-A和B’-B（图4a显示剖面位置）应用归一化总梯度法进行计算，极值最大值对应的深度和埋藏异常体的中心深度相关。从图4d的结果来看穿越C-1异常体的剖面A’-A下方异常中心深度约200 m，C-1异常体中心对应的埋藏深度可达200 m，穿越C-2异常体的剖面B’-B下方异常深度与剖面A’-A类似，因此C-2异常体埋藏深度可达200-300 m，但规模大于C-1异常体。
4.2 二维剖面模拟和三维反演模型的地球物理解译
进行二维剖面模拟的2条剖面A’-A和B’-B沿着东南-西北方向（图4a），贯穿调查区西侧C-1异常体和中心地区C-2异常体，剖面走向分别为315°和334°（图4f），模拟误差均小于2.62 nT（图5）。在跨越C-1异常体的剖面A’-A观察到一个规模相对较小的磁性异常体，由于上层新生代沉积物覆盖层的磁性接近0，背景磁化率设置为1× 10−5 SI，剖面A’-A的磁异常主要来源于C-1异常体。C-1异常幅值较小，假设异常体磁化率为
6000 × 10−5 SI，得到的二维总磁强度模型显示在约240 m的深度存在约120 m厚度的异常体（图5a）。而跨越C-2异常体的剖面B’-B观察到一个规模较大的磁性异常体，假设异常体磁化率为6000 × 10−5 SI时，获得的二维总磁强度模型显示在约280 m的深度存在约160 m厚度的异常体（图5b）。二维剖面模拟的结果与归一化总梯度法确定的场源深度基本是一致的。![]()
图 5 剖面A’-A和B’-B的二维总磁强度模型Figure 5. The two-dimensional total magnetic intensity model of profiles A-A' and B'-B.三维反演模型和可视化处理更好地揭示了隐伏异常体的形态特征和分布规律（图6），恢复的观测数据（图6a）和实际数据拟合的很好（图4a）。在红石岗南调查区，三维反演异常体的最大磁化率为
10000 × 10−5 SI，假设剩余磁化与当前地磁场方向一致时，将剩余磁化强度转换为磁化率，根据图2中磁性强度最强的矿化样品计算，该区域样品磁化率最大值为17000 × 10−5 SI，反演恢复的磁化率结果在这个范围内，因此是可信的。通过反距离插值法对反演结果进行处理，生~成三维可视化模型（图6b~图d）。可观察到异常体走向为近SW—NE向（图6b~图d），存在5个主要的异常体，异常体磁化率接近6000 × 10−5 SI，主要分布在红石岗南调查区中部和北部，形态接近长轴圆柱体，中部为调查区最大异常体（即C-2异常体），北部存在多个较小的异常体，异常体整体向北倾，总体上异常体深度可能位于新生代沉积物覆盖层之下，约200~300 m之间。![]()
图 6 （a）三维反演恢复的磁异常；（b）三维反演磁化率模型200 m深度（海拔高度1180 m）切片；（c）三维反演磁化率模型264 m深度（海拔高度1116 m）切片；（d）红石岗南调查区隐伏异常体三维可视化图像，显示阈值0.06 SIFigure 6. (a) Prediction of magnetic anomalies; (b) 200 m depth (Altitude1180 m) slice of 3D inversion susceptibility model; (c) 264 m depth (Altitude1116 m) slice of 3D inversion susceptibility model; (d) Three-dimensional visualization image showed a threshold value of 0.06 SI of the underground anomalous body in the southern Hongshigang area需要注意的是，反演时模型区域进行适当的扩大以减小边界效应，对反演数据进行补空，边界处未出现畸变（图6b~图6d），但反演结果垂向分辨率有待提高，存在一定的细节缺失。整体上，多种方法分析磁异常的结果对调查区隐伏异常体的水平和垂直位置限定相一致。
5. 找矿潜力分析
通过在红石岗南调查区开展的地面高精度磁测，并利用磁异常化极、归一化总梯度法、二维剖面模拟、三维磁异常反演确定了红石岗南调查区异常体的数量、磁化强度、埋深、产状要素及规模。化极数据显示存在5个正磁异常，走向普遍为近SW—NE向，异常体向东北侧倾斜且埋藏很深（图4a、图4b），北部和南部的串珠状异常指示调查区内存在2条近SW—NE向的断裂（图4b-c），其中北部为前人推断的黄山-镜儿泉断裂，南部为调查区内次级断裂（图4c）。利用归一化总梯度法确定了异常体的埋深可能达到200~300 m（图4d），二维模拟和三维可视化模型进一步对异常体的形态进行了精细刻画，异常体平均磁化率值约
6000 × 10−5SI，走向主要为近南西-北东向并向北侧倾斜。在对红石岗南异常体的岩性进行判断时，需要考虑该区域存在大规模新生代沉积物的覆盖，根据对北山裂谷带梧桐窝子组、甘泉组沉积物（He et al., 2021）和黄山~镜儿泉侏罗系煤窑沟群底砾岩（邵行来，2012）的磁性测试，其磁化率平均值小于100 × 10−5 SI，几乎没有磁性，因此红石岗南调查区的异常体与新生代沉积物和石炭统沉积岩无关。尽管该区域存在泥盆统火山岩，但除泥盆系下统大南湖组第四亚组凝灰岩、中基性火山凝灰岩显著高值（
13000 × 10−5 SI）外，泥盆系下统大南湖组的凝灰岩、火山角砾、粉砂岩大部分磁化率为100 × 10−5 SI和1000 × 10−5 SI量级，且此区域范围内泥盆统火山岩主要出露在土墩~镜儿泉断裂以北（图1b），所以红石岗北和红石岗的磁性异常体不是泥盆统火山岩。因此在红石岗区域异常体可能与闪长岩、花岗岩类及镁铁~超镁铁质岩体相关（图1b）。根据表1的统计结果，图拉尔根和黄山东的闪长岩的磁化率为12× 10−5~458 × 10−5 SI，剩磁量级为11× 10−3~150 × 10−3 A/m，远小于二维模拟和三维反演得到的异常体平均磁化率6000 × 10−5 SI。而花岗岩类的磁化率为19× 10−5~50 × 10−5 SI，剩磁量级为10× 10−3~50 × 10−3 A/m（表1），也远小于异常体平均磁化率6000 × 10−5 SI，东天山成矿带其他区域的闪长岩和花岗岩类也存在类似的磁性（邵行来等，2010; 刘璎等，2011; 惠卫东等，2011; 2012;Xiao, 2017; 吴功成，2018; 侯朝勇等，2021）。因此，红石岗南调查区的异常体不是闪长岩和花岗岩类。根据表1的结果，橄榄岩和矿化样品具有显著的高磁化率和剩磁值，结合区域地质情况，推测红石岗南调查区的隐伏异常体可能是镁铁-超镁铁质岩体，其平均磁化率约为
6000 × 10−5 SI，该区域新生代沉积物覆盖层厚度约50-140 m，下石炭统沉积岩厚度较大（李彤泰，2011），岩体埋深约200~300 m，走向主要为近SW—NE向，岩层向北侧倾斜。同时前人推断黄山-镜儿泉断裂通过红石岗南调查区（图1b），化极磁异常和垂向一阶导数的结果（图4b~图4c）指示了黄山-镜儿泉断裂断层的存在，同时发现了调查区内与黄山-镜儿泉断裂走向一致的次级断裂（图4c）。区域断裂构造意味着在基底内的薄弱地带，其为侵入体创造了侵入条件，超镁铁质岩等岩石的侵入对于矿床形成具有重要的控制作用，当超镁铁质岩侵入到地下时，其高温和高压作用会促进地下矿物的熔融和矿床的形成，同时改变地下流体的成分和运移路径，影响矿床的分布（Yurichev , 2017; 石煜等，2022），已在黄山-镜儿泉断裂带上发现多个铜镍矿床（图1a）。因此，在红石岗南地区可能存在镁铁-超镁铁岩和岩浆对中泥盆统火山岩和石炭统沉积岩（图1b）地广泛入侵，源区由岩石圈地幔和软流圈地幔组成，岩浆入侵期间遭受地壳物质的混染（王亚磊等，2017; 师震等，2019）。而侵位于干墩组的镁铁-超镁铁质岩体普遍规模较大，且含矿性较好（尹希文，2015; 宋谢炎等，2022），因此红石岗南调查区存在的超镁铁质岩体具有较好的找矿潜力。6. 结论
通过对红石岗南的高精度地面磁测数据的处理，利用磁异常化极、垂向一阶导数、归一化总梯度法、二维剖面模拟、三维磁异常反演等方法，发现红石岗南调查区呈现出5个磁性隐伏异常体，其平均磁化率约为
6000 × 10−5 SI，埋深大约在200~300 m范围内。异常体主要走向为近SW—NE向，并且显示出岩体向北侧倾斜的特征。磁异常反演结果揭示了该区域存在规模较大的异常体和延伸较长的断层，其中调查区中心存在C-2主要磁异常体，处于浅覆盖层之下，厚约160 m，与南部推测断层接壤。异常体可能由磁化率较高的镁铁-超镁铁质岩体构成。而推断出的2处断层具有明显与区域大构造走向一致的特征，与地质背景相符合，为解译黄山东至图拉尔根成矿带之间的浅覆盖区地质构造提供了重要线索。磁异常特征和反演结果显示，红石岗南区域可能存在潜在的镁铁-超镁铁质岩体，这些岩体具有较高的磁化率。此外，推测的断层也为铜镍矿床形成提供了有利条件。因此，红石岗南浅覆盖区内具有进一步寻找基性-超基性岩型铜镍矿床的潜力。
致谢：感谢各位审稿专家提出的宝贵意见！
-
![]()
图 10 H型模型电性源地空瞬变电磁法反演结果
a.反演结果图；b.反演拟合曲线
Figure 10. Results of H model of ground-airborne transient electromagnetic method with a long wire source
![]()
图 11 HA型模型电性源地空瞬变电磁法反演结果
a.反演结果图；b.反演拟合曲线
Figure 11. Results of HA model of ground-airborne transient electromagnetic method with a long wire source
![]()
图 16 测区地形图（来源google earth 2019）
Figure 16. Topographic map of survey area (from google earth 2019)
![]()
图 18 多源实测数据反演电阻率断面图
蓝色等值线起始值15 Ω·m；红色等值线终值165 Ω·m；间距15 Ω·m
Figure 18. Resistivity profile of multi-source data
-
范涛, 王秀臣, 李貅, 等. 瞬变电磁方法在探测煤矿浅层高阻采空区中的应用[J]. 西北地质, 2010, 43（2）: 156−162. FAN Tao, WANG Xiuchen, LI Xiu, et al. Application of TEM in Detecting Goaf of Coal Mine with High-resistivity and Shallow-layer[J]. Northwestern Geology,2010,43(2):156−162.
方涛, 张建军, 付成群, 等. 无人机地空瞬变电磁系统在冶山地下巷道探测中的应用[J]. 地球物理学进展, 2015, 30（5）: 2366−2372. FANG Tao, ZHANG Jianjun, FU Chengqun, et al. Using ground-airborne transient electromagnetic system on unmanned aerial vehicle detecting Yeshan underground tunnels[J]. Progress in Geophysics,2015,30(5):2366−2372.
嵇艳鞠, 王远, 徐江, 等. 无人飞艇长导线源时域地空电磁勘探系统及其应用[J]. 地球物理学报, 2013, 56（11）: 3640−3650. JI Yanju, WANG Yuan, XU Jiang, et al. Development and application of the grounded long wire source airborne electromagnetic exploration system based on an unmanned airship[J]. Chinese Journal of Geophysics,2013,56(11):3640−3650.
李肃义, 林君, 阳贵红, 等. 电性源时域地空电磁数据小波去噪方法研究[J]. 地球物理学报, 2013, 56（9）: 3145−3152. LI Suyi, LIN Jun, YANG Guihong, et al. Ground-Airborne electromagnetic signals de-noising using a combined wavelet transform algorithm[J]. Chinese Journal of Geophysics,2013,56(9):3145−3152.
李貅, 胡伟明, 薛国强. 多辐射源地空瞬变电磁响应三维数值模拟研究[J]. 地球物理学报, 2021, 64（2）: 716−723. LI Xiu, HU Weiming, XUE Guoqiang. 3D modeling of multi-radiation source semi-airborne transient electromagnetic response[J]. Chinese Journal of Geophysics,2021,64(2):716−723.
李貅, 张莹莹, 卢绪山, 等. 电性源瞬变电磁地空逆合成孔径成像[J]. 地球物理学报, 2015, 58（1）: 277−288. LI Xiu, ZHANG Yingying, LU Xushan, et al. Inverse Synthetic Aperture Imaging of Ground-Airborne transient electromagnetic method with a galvanic source[J]. Chinese Journal of Geophysics,2015,58(1):277−288.
李貅. 瞬变电磁测深的理论与应用[M]. 西安: 陕西科学技术出版社, 2002. LI Xiu. Theory and application of transient electromagnetic sounding [M]. Xi'an: Shaanxi Science Technology Press, 2002.
王凯, 刘宽厚, 张旭. 瞬变电磁法寻找铝土矿效果与分析[J]. 西北地质, 2010, 43（2）: 80−84. WANG Kai, LIU Kuanhou, ZHANG Xu. Effect and Analysis of Searching for Bauxite Deposit Using Transient Electromagnetic Method[J]. Northwestern Geology,2010,43(2):80−84.
薛国强, 李貅, 底青云. 瞬变电磁法正反演问题研究进展[J]. 地球物理学进展, 2008, 23（4）: 1165−1172. XUE Guoqiang, LI Xiu, DI Qingyun. Research progress in TEM forward modeling and inversion calculation[J]. Progress in Geophysics,2008,23(4):1165−1172.
张旭, 刘宽厚, 李貅. CSAMT的静校正应用——联合反演法[J]. 西北地质, 2010, 43（2）: 38−43. ZHANG Xu, LIU Kuanhou, LI Xiu. Application of Static Correction CSAMT——Joint inversion method[J]. Northwestern Geology,2010,43(2):38−43.
张莹莹, 李貅, 李佳, 等. 多辐射场源地空瞬变电磁法快速成像方法研究[J]. 地球物理学进展, 2016, 31（2）: 869−876. ZHANG Yingying, LI Xiu, LI Jia, et al. Fast imaging technique of multi-source ground-airborne transient electromagnetic method[J]. Progress in Geophysics,2016,31(2):869−876.
周道卿, 谭林, 谭捍东, 等. 频率域吊舱式直升机航空电磁资料的马奎特反演[J]. 地球物理学报, 2010, 53（2）: 421−427. ZHOU Daoqing, TAN Lin, TAN Handong, et al. Inversion of frequency domain helicopter-borne electromagnetic data with Marquardt’s method[J]. Chinese Journal of Geophysics,2010,53(2):421−427.
Allah S A, Mogi T, Ito H, et al. Three-dimensional resistivity modeling of GREATEM survey data from Kujukuri beach, Japan[J]. Proceedings of the 10th SEGJ International Symposium, 2011, 314-317.
Allah S A, Mogi T, Ito H, et al. Three-dimensional resistivity characterization of a coastal area: application of grounded electrical source airborne transient electromagnetic(GREATEM) survey data from Kujukuri beach, Japan[J]. Journal of Applied Geophysics,2013,99:1−11. doi: 10.1016/j.jappgeo.2013.09.011
Bryan W, Roger E. Trenton. Resistivity Arrays as an Early Warning System for Monitoring Runoff Holding Ponds[J]. Journal of Environmental and Engineering Geophysics,2015,20:319−335. doi: 10.2113/JEEG20.4.319
Cedrina L, Bonomo N, Osella A. An application of the synthetic emitter-array method to improve GPR signals[J]. Journal of Applied Geophysics,2010,70:237−244. doi: 10.1016/j.jappgeo.2010.01.004
Ito H, Mogi T, Jomori A, et al. Further invertigations of underground resistivity structures in coastal areas using grounded source airborne edectromagnetics[J]. Earth, planets and space,2011,63(8):9−12. doi: 10.5047/eps.2011.08.003
Ito H, Kaieda H, Mogi T, et al. Grounded electrical-source airborne transient electromagnetic (GREATEM) survey of Aso Volcano, Japan[J]. Exploration Geophysics (CSIRO PUBLISHING),2014,45(1):43−48. doi: 10.1071/EG12074
Li H, QI Z, et al. Numerical modelling analysis of multi-source semi-airborne TEM systems using a TFEM[J]. Journal of Geophysics and Engineering,2020,17:399−410. doi: 10.1093/jge/gxz119
Li J, Hu X, Zeng S, et al. Three-dimensional forward calculation for loop source transient electromagnetic method based on electric field Helmholtz equation, Chinese journal of geophysics, 2013, 56: 4256–4267.
Mogi T, Tanaka Y, Kusunoki K, et al. Development of grounded electrical-source airborne transient EM (GREATEM)[J]. Exploration Geophysics,1998,29:61−64. doi: 10.1071/EG998061
Mogi T, Kusunoki K, Kaieda H, et al. Grounded electrical-source airborne transient electromagnetic (GREATEM) survey of Mount Bandai, north-eastern Japan[J]. Exploration Geophysics,2009,40:1−7. doi: 10.1071/EG08115
Smith R S, Annan A P, McGowan P D. A comparison of data from airborn, semi-airborne, and ground electromagnetic sounding method[J]. Geophysics,2001,66(5):1379−1385. doi: 10.1190/1.1487084
Wright D A, Ziolkowski. Hydrocarbon detection and monitoring with a multicomponent transient electromagnetic (MTEM) survey[J]. The Leading Edge,2002,21:852−864. doi: 10.1190/1.1508954
ZHANG Xu, HAN Yuan, ZhANG Yongxing. Image Identification and Land Use Management of Anomalous Minerals from Geophysical Information and Deep Convolutional Networks[J]. The Journal of Geology,2023,131(4):325−337.
下载:
计量
- 文章访问数: 8
- HTML全文浏览量: 0
- PDF下载量: 4
