Characteristics of Magnetic Anomalies and Geologic Significance in the Southern Hongshigang Area of the Eastern Tianshan Mountains, China
-
摘要:
近些年,金属矿产勘查重点逐渐向岩体未出露的覆盖区隐伏矿床转移,但是在前期的调查工作中,由于覆盖层较厚很难得到地下异常体信息。位于东天山铜镍成矿带的红石岗一带出露的镁铁-超镁铁质岩体显示了良好的找矿潜力,1∶5万航磁和重力资料显示其南部第四系覆盖区分布一高磁高重力套合异常,目前对该异常尚无开展调查研究。因此本研究在红石岗南调查区进行了高精度地面磁测工作,采用磁异常化极、垂向一阶导数、归一化总梯度法、二维剖面模拟和三维磁异常反演方法对该区域内隐伏异常体进行初步勘探研究,发现调查区内存在的5处显著磁异常,隐伏异常体的平均磁化率约
6000 × 10−5 SI,埋深约200~300 m,存在5个主要的小异常体,走向为近SW—NE向,岩层向北侧倾斜。结合地质资料和岩石物性判断,该处异常体可能是镁铁-超镁铁质岩体,且磁化率较高。同时利用磁异常指出与前人推断位置一致的近SW—NE向黄山-镜儿泉深大断层和调查区内次级断裂,可为铜镍矿床提供有利的成矿环境。因此,红石岗南调查区的隐伏异常体存在较大的铜镍矿找矿潜力,同时本研究为后续在该区域所开展的勘探工作,提供了可信的地球物理资料。Abstract:In recent years, the emphasis of metal mineral exploration has gradually shifted to the hidden deposits in the covered areas where the rock mass is not exposed. However, it is difficult to obtain the information of underground abnormal bodies in the early investigation work because of the thick overlying layer. The Mafic and ultramafic rock mass outlying in the Hongshigang area of the Eastern Tianshan Cu-Ni metallogenic belt shows good prospecting potential. The 1:50 000 aeromagnetic data and gravity data show that there is a high magnetic and high gravity coupling anomaly in the southern Quaternary covered area, which has not been investigated yet. Therefore, this study conducted high-precision ground magnetic survey work in the southern Hongshigang area and used magnetic anomaly reduced to the pole, vertical first derivative, normalized total gradient method, two-dimensional profile simulation, and three-dimensional magnetic anomaly inversion methods to conduct preliminary exploration and research on the underground anomalous bodies in the area. It was found that there are five significant magnetic anomalies in the investigation area, and the average magnetic susceptibility of the hidden anomalous bodies is about
6000 × 10−5 SI, with a depth of about 200-300 meters. There are five main small anomalous bodies, which strike nearly southwest-northeast and dip to the north. Combined with geological data and rock physical properties, the anomalous bodies in this area may be mafic and ultramafic rock bodies with high magnetic susceptibility. Magnetic anomalies are used to indicate the Huangshan-Jingerquan faults in the southwest-northeast direction which are consistent with the inferred position of predecessors and the secondary fault in the investigated area, which could provide a favorable metallogenic environment for copper-nickel deposits. Therefore, the underground anomalous bodies in the southern Hongshigang area have great prospecting potential for copper-nickel deposits, and this study provides reliable geophysical data for subsequent exploration work in the area. -
新疆东天山铜镍成矿带是中国重要的铜镍矿集区(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)。二维剖面模拟的结果与归一化总梯度法确定的场源深度基本是一致的。三维反演模型和可视化处理更好地揭示了隐伏异常体的形态特征和分布规律(图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处断层具有明显与区域大构造走向一致的特征,与地质背景相符合,为解译黄山东至图拉尔根成矿带之间的浅覆盖区地质构造提供了重要线索。磁异常特征和反演结果显示,红石岗南区域可能存在潜在的镁铁-超镁铁质岩体,这些岩体具有较高的磁化率。此外,推测的断层也为铜镍矿床形成提供了有利条件。因此,红石岗南浅覆盖区内具有进一步寻找基性-超基性岩型铜镍矿床的潜力。
致谢:感谢各位审稿专家提出的宝贵意见!
-
图 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
图 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
图 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 (Altitude
1180 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表 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 图拉尔根岩心样品 -
宫辰 . 哈密红石岗铜镍矿矿床地质特征及找矿前景[J]. 中国金属通报,2020 (11 ):49 −50 . doi: 10.3969/j.issn.1672-1667.2020.21.023韩宝福, 季建清, 宋彪, 等 . 新疆喀拉通克和黄山东含铜镍矿镁铁-超镁铁杂岩体的SHRIMP锆石U-Pb年龄及其地质意义[J]. 科学通报,2004 (22 ):2324 −2328 . doi: 10.3321/j.issn:0023-074X.2004.22.012侯朝勇, 蔡厚安, 裴森龙 . 综合物化探方法在新疆哈密月牙湾铜镍矿勘查中的应用[J]. 矿产与地质,2021 ,35 (6 ):1116 −1123 .惠卫东, 赵鹏大, 秦克章, 等 . 东天山图拉尔根铜镍硫化物矿床综合信息找矿模型的应用[J]. 地质与勘探,2011 ,47 (3 ):388 −399 .李彤泰 . 新疆哈密市黄山基性-超基性岩带铜镍矿床地质特征及矿床成因[J]. 西北地质,2011 ,44 (1 ):54 −60 . doi: 10.3969/j.issn.1009-6248.2011.01.007刘璎, 孟贵祥, 严加永, 等 . 重磁3D物性反演技术在金属矿勘探中的应用[J]. 地质与勘探,2011 ,47 (3 ):448 −455 .刘隆, 周建平, 吴涛, 等 . 大洋中脊玄武岩磁性特征[J]. 地球物理学进展,2021 ,36 (5 ):1880 −1890 . doi: 10.6038/pg2021EE0403乔天成 . 高精度磁法在铜镍多金属矿普查工作中的应用[J]. 新疆有色金属,2016 ,39 (1 ):36 −39 .邵行来. 东天山黄山—镜儿泉超镁铁岩带地球物理特征研究及找矿应用[D]. 中国地质大学(北京), 2012. 邵行来, 薛春纪, 戴德文, 等 . 新疆哈密葫芦岩浆Cu-Ni矿勘查地球物理异常特征[J]. 现代地质,2010 ,24 (2 ):383 −391 . doi: 10.3969/j.issn.1000-8527.2010.02.025邵行来, 薛春纪, 周耀明 . 哈密图拉尔根镁铁-超镁铁岩磁法异常解释[J]. 新疆地质,2012 ,30 (4 ):425 −429 . doi: 10.3969/j.issn.1000-8845.2012.04.009师震, 陈宏骏, 钱壮志, 等 . 东天山红石岗镁铁—超镁铁质岩体成因及铜镍成矿潜力[J]. 地球科学与环境学报,2019 ,41 (2 ):156 −169 . doi: 10.3969/j.issn.1672-6561.2019.02.003宋谢炎, 邓宇峰, 颉炜, 等 . 新疆黄山-镜儿泉铜镍硫化物成矿带岩浆通道成矿特征及其找矿意义[J]. 矿床地质,2022 ,41 (6 ):1108 −1123 .王成 . 重、磁、电综合勘探方法在寻找铜镍矿中的应用[J]. 新疆有色金属,2018 ,41 (4 ):13 −17 .王庆功 . 综合物探方法在新疆某铜镍矿勘探中的应用[J]. 甘肃冶金,2021 ,43 (2 ):96 −98+104 . doi: 10.3969/j.issn.1672-4461.2021.02.028王志福, 吴飞, 谭克彬, 等 . 哈密红石岗铜镍矿矿床地质特征及找矿前景[J]. 新疆地质,2012 ,30 (3 ):307 −311 . doi: 10.3969/j.issn.1000-8845.2012.03.013王亚磊, 张照伟, 陈寿波, 等 . 新疆东天山红石岗北铜镍矿化镁铁质岩体岩石成因及成矿潜力分析[J]. 地质学报,2017 ,91 (4 ):776 −791 . doi: 10.3969/j.issn.0001-5717.2017.04.006吴功成. 新疆白石泉铜镍矿矿床地质与找矿预测[D]. 中国地质大学(北京), 2018. . 肖丹, 宋泽友, 宋维国 . 我国岩浆硫化物型镍矿床伴生矿产综合勘查评价指标探讨[J]. 国土资源导刊,2022 ,19 (1 ):48 −53 . doi: 10.3969/j.issn.1672-5603.2022.01.011杨大欢, 古志宏 . 广东与燕山期岩浆作用有关矿产资源的区域成矿分带特征及成因[J]. 矿产与地质,2021 ,35 (4 ):603 −609 .尹希文 . 新疆香山铜镍硫化物矿床岩浆深部过程与找矿方向探讨[J]. 西北地质,2015 ,48 (3 ):22 −30 . doi: 10.3969/j.issn.1009-6248.2015.03.003Cockett R, Kang S, Heagy L J, et al . SimPEG: An open source framework for simulation and gradient based parameter estimation in geophysical applications[J]. Computers & Geosciences,2015 ,85 :142 −154 .Ekwok S E, Achadu O-I M, Akpan A E, et al . Depth Estimation of Sedimentary Sections and Basement Rocks in the Bornu Basin, Northeast Nigeria Using High-Resolution Airborne Magnetic Data[J]. Minerals,2022 ,12 (3 ):285 .Ekwok S E, Akpan A E, Ebong E D . Enhancement and modelling of aeromagnetic data of some inland basins, southeastern Nigeria[J]. Journal of African Earth Sciences,2019 ,155 :43 −53 . doi: 10.1016/j.jafrearsci.2019.02.030Elysseieva I s., Pašteka R . Review Paper: Historical development of the total normalized gradient method in profile gravity field interpretation[J]. Geophysical Prospecting,2019 ,67 (1 ):188 −209 . doi: 10.1111/1365-2478.12704Feng Y, Qian Z, Duan J, et al . Geochronological and geochemical study of the Baixintan magmatic Ni-Cu sulphide deposit: New implications for the exploration potential in the western part of the East Tianshan nickel belt (NW China)[J]. Ore Geology Reviews,2018 ,95 :366 −381 . doi: 10.1016/j.oregeorev.2018.02.023He J, Fan Z, Xiong S, et al . Geophysical prospecting of copper-nickel deposits in Beishan rift zone, Xinjiang[J]. China Geology,2021 ,4 (1 ):126 −146 .Jahren C E . Magnetization of keweenawan rocks near duluth, minnesota[J]. Geophysics,1965 ,30 (5 ):858 . doi: 10.1190/1.1439660Koenigsberger J G . Natural residual magnetism of eruptive rocks[J]. Terrestrial Magnetism and Atmospheric Electricity,1938 ,43 (3 ):299 −320 . doi: 10.1029/TE043i003p00299Mao Y-J, Qin K-Z, Li C, et al. Petrogenesis and ore genesis of the Permian Huangshanxi sulfide ore-bearing mafic-ultramafic intrusion in the Central Asian Orogenic Belt, western China[J]. Lithos, 2014, 200–201: 111–125.
Mao Y-J, Qin K-Z, Tang D-M, et al . Crustal contamination and sulfide immiscibility history of the Permian Huangshannan magmatic Ni-Cu sulfide deposit, East Tianshan, NW China[J]. Journal of Asian Earth Sciences,2016 ,129 :22 −37 . doi: 10.1016/j.jseaes.2016.07.028Meng Z . Three-dimensional potential field data inversion with L0 quasinorm sparse constraints[J]. Geophysical Prospecting,2018 ,66 (3 ):626 −646 .Parkinson W D, Barnes C D . In situ determination of Koenigsberger ratio[J]. Australian Journal of Earth Sciences,1985 ,32 (1 ):1 −5 . doi: 10.1080/08120098508729308Sun J, Li Y . Adaptive Lp inversion for simultaneous recovery of both blocky and smooth features in a geophysical model[J]. Geophysical Journal International,2014 ,197 (2 ):882 −899 . doi: 10.1093/gji/ggu067Sun J, Li Y . Multidomain petrophysically constrained inversion and geology differentiation using guided fuzzy c-means clustering: 4[J]. Geophysics,2015 ,80 (4 ):ID1 −ID18 .Utsugi M. 3-D inversion of magnetic data based on the L1–L2 norm regularization[J]. Earth, Planets and Space, 2019, 71(1).
Vanzon T, Roy-Chowdhury K. Structural inversion of gravity data using linear programming: 3[J]. Geophysics, 2006, 71(3).
Wardinski I, Saturnino D, Amit H, et al . Geomagnetic core field models and secular variation forecasts for the 13th International Geomagnetic Reference Field (IGRF-13)[J]. Earth Planets and Space,2020 ,72 (1 ):155 . doi: 10.1186/s40623-020-01254-7Xiangjin R, Linfu X, Yanyan Z, et al. The 3D Visualization of 2D GM-SYS Gravity-Magnetic Inversion Sections Based on GoCAD[C]//2017 International Conference on Robots & Intelligent System (ICRIS). Huai An City, China: IEEE, 2017. 325–328.
Xiao F, Wang Z . Geological interpretation of Bouguer gravity and aeromagnetic data from the Gobi-desert covered area, Eastern Tianshan, China: Implications for porphyry Cu-Mo polymetallic deposits exploration[J]. Ore Geology Reviews,2017 ,80 :1042 −1055 . doi: 10.1016/j.oregeorev.2016.08.034Yurichev A N, Chernyshov A I . New Ore Minerals from the Kingash Ultramafic Massif, Northwestern Eastern Sayan[J]. Geology of Ore Deposits,2017 ,59 (7 ):626 −631 . doi: 10.1134/S107570151707011X石煜, 王玉往, 王京彬, 等 . 东天山黄山东和黄山西铜镍硫化物矿床含矿超镁铁岩的成岩-成矿作用机制: 来自斜长石成分的约束[J]. 地球科学,2022 ,47 (9 ):3244 −3257 . doi: 10.3321/j.issn.1000-2383.2022.9.dqkx202209009