Molybdenite Re-Os isochron age of the Chiwawu Mo-Cu deposit in Fuping County, Hebei Province in northern Taihang Mountains and its prospecting significance
-
摘要:
赤瓦屋斑岩型钼铜矿床是太行山北段斑岩-矽卡岩成矿带内探明的一个中型钼铜矿床，钼铜矿体主要赋存于赤瓦屋岩体内部相的斑状花岗闪长岩中。为厘定赤瓦屋钼铜矿床的成矿时代和成矿物质来源，利用辉钼矿Re-Os同位素定年方法对赤瓦屋钼铜矿床中7件辉钼矿样品进行了成矿年代学测定，获得的模式年龄为129.1~130.6 Ma，年龄加权平均值为129.7±0.7 Ma，对应的等时线年龄为128.7±4.4 Ma，模式年龄和等时线年龄在误差范围内基本一致，指示赤瓦屋钼铜矿床形成于早白垩世。赤瓦屋钼铜矿床辉钼矿样品的Re含量为25.9×10-6~37.1×10-6，表明其成矿物质来源于壳幔混源。结合区域年代学资料，认为太行山北段成矿带在134~124 Ma存在一期岩浆热液成矿事件，斑岩铜矿型钼铜矿床的找矿潜力较大。
Abstract:The Chiwawu porphyry Mo-Cu deposit in the northern Taihang Mountains porphyry-skarn metallogenic belt is a medium-sized one discovered in recent years.The Mo-Cu orebodies are mainly hosted within porphyritic granodiorite of internal facies in Chiwawu pluton.In order to determine the age of mineralization and the source of molybdenum, seven molybdenite samples were collected from the Chiwawu Mo-Cu deposit for Re-Os dating.The Re-Os dating yields model ages ranging from129.1 Ma to130.6 Ma, with a weighted mean age of 129.7±0.7 Ma, and obtains an isochron age of 128.7±4.4 Ma.Evidently, model ages are consistent with isochron age within the error ranges, which constrains the ore-forming age of Chiwawu Mo-Cu deposit at the Early Cretaceous.The Re values of molybdenum vary in the range of 25.9×10-6~37.1×10-6, indicating that they were derived from mixing of crust and mantle sources.Combined with available chronologic data, it is recognized that a regional magmatic-metallogenic event might exist during the 134~124 Ma in northern Taihang Mountains.The prospecting potentiality of porphyry Cu-Mo deposits is great in this area.
-
IOCG型矿床，即铁氧化物-铜-金矿床，指铁氧化物含量大于20%的铜-金矿床，主要分布在澳大利亚、巴西、智利、加拿大、美国、南非等国及中国内蒙古、海南、四川、长江中下游等地。其一般具有规模大、品位高、元素多、埋藏浅、易采选等特点，是近年来颇受国内外矿业界、勘查界、矿床界等地质领域重视的一种矿床组合类型[1-7]。由于IOCG型矿床富含铁氧化物，常缺失硫化物，且蚀变范围广阔，地球物理特征明显(成矿区的磁场和重力效应明显，具有重力高、中-高幅磁异常为标志)。因此，高精度磁测往往是寻找IOCG矿床重要的有效手段之一。
对于磁测数据的处理和解译是找矿预测的重要环节。其中RGIS是中国地质调查局发展研究中心牵头研发的物探重磁电数据处理与解释软件系统。与国内外重磁专业数据处理软件相比，RGIS系统具有重磁数据整理齐备、规范，数据预处理功能丰富，重磁数据处理功能全面、实用，数据范围和格式广泛，输入输出方便、快捷、规范等特点。目前已成为中国用户最多、覆盖面最广、产学研多方应用率最高的物探专业软件，在各类地质矿产和油气勘查与评价工作中发挥着重要作用[8-22]。
为探讨高精度磁测在IOCG型铁矿勘查中的应用效果，本文以智利英格瓦塞铁矿为例，在矿区1:5000地面高精度磁测基础上，利用RGIS软件重点针对中部矿区的磁测数据，进行面积上延、面积下延、精测剖面下延及2.5D反演拟合处理，初步确定了矿体的空间分布，为后期地质找矿勘探工作奠定了很好的基础。同时，通过初步对比海南石碌铁矿，发现IOCG型铁矿皆具有富含铁氧化物的特征，磁场特征明显，因而该方法也可进一步推广应用于中国IOCG型铁矿勘查。
1. 区域地质背景
智利英格瓦塞铁矿位于美洲大陆西缘安第斯褶皱山系南美褶皱带构造区内中生代格拉岛弧弧前盆地。区域出露地层主要为中生界上三叠统—下侏罗统水之歌组、下白垩统大十二弦琴组和查尼拉尔组、上白垩统塞里约斯组、新生界古近系—新近系和第四系。区域构造主要由褶皱和断层构成。其中褶皱主要分布在区域西北部；区域内断裂构造相对发育，主要为SN向断层(区域阿塔卡玛断裂近SN向贯穿本区)，其次为NE向、NW向断层。区域内岩浆活动强烈，主要为燕山期中酸性侵入体，岩体贯穿区域南北，规模巨大，呈岩基产出，出露面积占区域面积30%以上。区域内矿产资源丰富，主要矿种为铜、金、铁。铁氧化物铜金型(IOCG)矿床分布广泛(图 1)。
![]()
图 1 智利英格瓦塞铁矿中部矿区地质图(区域位置图据参考文献[1]修改)Figure 1. Geological map of the central mining area of the Incaguasi iron deposit in Chile2. 矿区地质特征
智利英格瓦塞铁矿属于IOCG型矿床，由中部、北部、西南部3个矿区构成，中部矿区是其主要部分，位于智利阿塔卡玛大区巴耶纳尔市英格瓦塞镇北西约1.5 km处，隶属英格瓦塞镇管辖。
矿区地层出露主要为中生界下白垩统大十二弦琴组(其岩性为安山岩)和第四系(冲积物)。矿区地层总体为一单斜构造，走向NE，倾向SE，倾角约为16°，局部有扭曲现象，产状略有变化。由于本区工作程度尚低，未见明显的断裂等构造。区内岩浆岩极发育，主要发育花岗闪长岩岩体、花岗斑岩岩脉及辉石角闪石岩岩体。其中，区内广泛发育花岗闪长岩，其侵入大十二弦琴组中，使地层出露分布较凌乱，同时亦侵入辉石角闪石岩岩体中，野外可见辉石角闪石岩岩体以捕虏体形式赋存于花岗闪长岩岩体中，两者之间界线清晰。辉石角闪石岩主要分布在矿区中部，局部根据含磁铁矿异同，可分为含磁铁辉石角闪石岩、磁铁角闪石岩、角闪磁铁岩。花岗斑岩岩脉主要分布在中部矿区的东北部，总体呈NW向产出，侵入于花岗闪长岩及含磁铁角闪石岩中。此外还出露少量正长岩，呈岩株产出(图 1)。
矿区目前发现4条规模较大的矿体，主要产于含磁铁辉石角闪石岩中，围岩基本为花岗闪长岩。矿体总体呈带状，走向NE，倾向NW，倾角较陡，多70°~90°，单个矿体长约几百米至上千米，宽约几十米至上百米。矿石类型可分为含磁铁辉石角闪石岩、磁铁角闪石岩，含角闪磁铁岩。矿石主要为细粒结构，块状构造。矿石矿物成分较简单，主要为磁铁矿，其次为赤铁矿，局部见孔雀石、蓝铜矿。脉石矿物以角闪、辉石为主。地表浅部围岩矿化与蚀变主要有褐铁矿化、磁铁矿化、赤铁矿化、镜铁矿化、孔雀石化、蓝铜矿化(多在岩石裂隙面发育)、高岭土化、绿泥石化等。
3. 中部矿区磁异常特征
从区域物性参数K(磁化率，单位为4π×10-6SI)和Jr(剩磁强度，单位为10-3A/m)可知，本区花岗闪长岩磁性一般为K=37.5~3046，Jr=14.74~350；花岗斑岩磁性一般为K=23.46~3792，Jr=12.24~228；磁铁角闪石岩为K=2000~120000，Jr=2000~160000；正长斑岩K=0.41~2000，Jr=0.01~30；安山岩K=5~3000，Jr=5~400；含磁铁辉石角闪石岩一般为K=10~10000，Jr=10~3000；辉石角闪石岩一般为K=0.1~1000，Jr=0.1~400。磁性最强的为磁铁角闪石岩，其次为含磁铁辉石角闪石岩，与围岩有明显的磁性差异，因此，本区采用磁法手段勘查非常有效。
根据1:5000地面高精度磁测ΔT等值线图，结合本地区地质及岩石磁性特征，有意义的异常区圈定为大于或等于1600 nT。该工区内的ΔT最大值为5173.3 nT，等值线总体呈NE向延伸。该磁异常区内含7个大于1600 nT的规模较大的磁异常圈，编号分别为Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ、Ⅶ(图 2)。
![]()
图 2 智利英格瓦塞中部矿区地面磁测ΔT等值线图Figure 2. ΔT contour map of ground magnetic survey for the central mining area of the Incaguasi iron deposit in ChileⅠ号磁异常位于中部矿区的北部，异常中心ΔT最大值为3273.1 nT。该磁异常呈不规则椭圆状，长轴呈NE向，长321 m，短轴最大宽度123 m，规模为40936 m2。该磁异常东南翼等值线较密，梯度变化较大，北西翼梯度变化较小，定性判断引起该磁异常的磁性地质体呈NE走向、向NW向倾斜。
Ⅱ号磁异常位于Ⅰ号磁异常的西南部，异常中心ΔT最大值为2229.3 nT。该磁异常呈不规则椭圆状，长轴呈EW向，长轴长340 m，短轴最大宽度为97 m，规模为39856 m2。该磁异常两翼较对称，初步认定引起该磁异常的磁性地质体呈EW走向、倾角较直立。
Ⅲ号磁异常位于矿区的中西部，异常中心ΔT最大值为2423.8 nT。该异常未完全封闭，长轴呈NE向，长轴长500 m，短轴最大宽度180 m，规模为98163 m2。该异常两翼较对称，初步认定引起该磁异常的磁性地质体呈NE走向、倾角较直立。
Ⅳ号磁异常位于Ⅲ号磁异常的南部，异常中心ΔT最大值为2814.6 nT。该磁异常呈不规则长圆状，长轴呈NE向，长370 m，短轴最大宽度120 m，规模为53993 m2。该磁异常南部为负异常且梯度变化大，初步判定引起该磁异常的磁性地质体呈NE走向、向NW倾斜。
Ⅴ号磁异常位于矿区的南部，异常中心ΔT最大值为2511.3 nT。该异常呈不规则长圆状，长轴呈NE向，长190 m，短轴最大宽度64 m，规模为14396 m2。该磁异常两翼较对称，初步认定引起该磁异常的磁性地质体呈NE走向、倾角较直立。
Ⅵ号磁异常位于矿区的东南部，异常中心ΔT最大值为4339.6 nT。该磁异常呈不规则椭圆状，长轴呈NE向，长660 m，短轴最大宽度170 m，规模为75645 m2。该磁异常南部为负异常且梯度变化大，初步判定引起该磁异常的磁性地质体呈NE走向、向NW倾斜。
Ⅶ号磁异常位于矿区的东北部，异常中心ΔT最大值为5173.3 nT。该磁异常呈椭圆状，长轴呈NE向，长300 m，短轴最大宽度130 m，规模为34740 m2。该磁异常西南翼等值线较密，梯度变化较大，东北部梯度变化较小，初步判定引起该磁异常的磁性地质体呈NE走向、向NW倾斜。
4. 中部矿区磁异常解释
磁异常是各种不同性质、不同深度、不同形态、不同规模的磁性体磁场叠加在一起反应的结果。但是在原始图件上很难识别，给地质解释及验证工作带来了难度。为了进一步提高对磁异常的分辨能力，突出更多有用信息，根据测区磁异常特征和地质解释需要，利用RGIS2010对智利英格瓦塞中部矿区磁测数据进行了面积上延、面积下延、剖面下延及2.5D反演拟合处理。
4.1 面积上延解释
面积上延计算主要目的是压制浅部干扰异常或规模较小的局部异常，突出埋藏较深和规模较大的磁异常。对英格瓦塞中部矿区面积磁测分别进行了向上延拓25 m、50 m、75 m、100 m的数据处理，经统计各矿体ΔT等值线分布情况，如表 1所示。
表 1 各地磁异常上延中心变化统计结果Table 1. Statistical table of the center variation of magnetic anomalies'upward continuation磁异常编号 延拓高度/m 中心异常/nT 与原始曲线差值/nT 每25 m变化值/nT 上延解译 Ⅰ 0 3273.1 0 0 Ⅰ号磁异常区峰值衰减速度快，说明引起该异常的磁性地质体埋藏较浅，规模较小。初步认为该磁性地质体为有限延深 +25 2024.9 -1248.2 -1248.2 +50 1548.2 -1724.9 -476.7 +75 1315.1 -1958 -233.1 +100 1175.3 -2097.8 -139.8 Ⅱ 0 2229.3 0 0 Ⅱ号磁异常区峰值衰减速度较慢，说明引起该异常的磁性地质体埋藏有一定的深度，磁性体有一定的规模。初步认为该磁性地质体为有限延深 +25 1725 -504.3 -504.3 +50 1505.9 -723.4 -219.1 +75 1337.3 -892 -168.6 +100 1277 -952.3 -60.3 Ⅲ 0 2423.8 0 0 Ⅲ号磁异常区峰值衰减速度慢，说明引起该异常的磁性地质体埋藏有深度，磁性体规模较大。初步认为该磁性地质体为有限延深 +25 1934.7 -489.1 -489.1 +50 1710.4 -713.4 -224.3 +75 1557 -866.8 -153.4 +100 1443.7 -980.1 -113.3 Ⅳ 0 2814.6 0 0 Ⅳ号磁异常区峰值衰减速度较慢，说明引起该异常的磁性地质体埋藏有一定深度，磁性体有一定的规模。初步认为该磁性地质体为有限延深 +25 2000.2 -814.4 -814.4 +50 1619.5 -1195.1 -380.7 +75 1410.9 -1403.7 -208.6 +100 1283.8 -1530.8 -127.1 Ⅴ 0 2511.3 0 0 Ⅴ号磁异常区峰值衰减速度非常快，说明引起该异常的磁性地质体埋藏较浅，规模较小。初步认为该磁性地质体为有限延深 +25 1751.2 -760.1 -760.1 +50 异常消失 +75 异常消失 +100 异常消失 Ⅵ 0 4339.6 0 0 Ⅵ号磁异常区峰值衰减速度慢，说明引起该异常的磁性地质体埋藏有深度，磁性体规模大 +25 2667 -1672.6 -1672.6 +50 1933.2 -2406.4 -733.8 +75 1641.1 -2698.5 -292.1 +100 1425.8 -2913.8 -215.3 Ⅶ 0 5173.3 0 0 Ⅶ号磁异常区峰值衰减速度快，说明引起该异常的磁性地质体埋藏浅，磁性体规模大 +25 2946.1 -2227.2 -2227.2 +50 1916.3 -3257 -1029.8 +75 1352.1 -3821.2 -564.2 +100 1013 -4160.3 -338.8 从地面磁测上延25 m、50 m、75 m、100 m的ΔT等值线统计结果(表 1)明显看出，随着向上高度的增加，Ⅰ、Ⅴ磁异常中心值迅速减小，且磁异常消失，说明引起Ⅰ、Ⅴ磁异常的磁性地质体规模小，埋藏浅且为有限延深；Ⅱ、Ⅳ、Ⅶ磁异常中心值衰减较慢，且磁异常有一定的规模，说明引起Ⅱ、Ⅳ、Ⅶ磁异常的磁性地质体有一定的规模，埋藏浅且为有限延深；Ⅲ、Ⅵ磁异常中心值衰减慢，且磁异常规模较大，说明引起Ⅲ、Ⅵ磁异常的磁性地质体有一定规模，埋藏有深度且有一定的延深。
4.2 面积下延解释
向下延拓是由实测磁场向磁源方向的延拓，利用向下延拓可以处理旁侧叠加异常，将多个单独异常从叠加异常中区分出来，同时延拓也可以用于分离上下叠加异常。对英格瓦塞中部矿区面积磁测分别进行了向下延拓25 m、50 m、75 m的数据处理，其磁异常分布情况如下(图 3、图 4、图 5)。
![]()
图 3 智利英格瓦塞中部矿区地面磁测ΔT下延25 m等值线图Figure 3. ΔT contour map with downward continuation 25 m of ground magnetic survey for the central mining area of the Incaguasi iron deposit in Chile![]()
图 4 智利英格瓦塞中部矿区地面磁测ΔT下延50 m等值线图Figure 4. ΔT contour map with downward continuation 50 m of ground magnetic survey for the central mining area of the Incaguasi iron deposit in Chile![]()
图 5 智利英格瓦塞中部矿区地面磁测ΔT下延75 m等值线图Figure 5. ΔT contour map with downward continuation 75 m of ground magnetic survey for the central mining area of the Incaguasi iron deposit in Chile从图 3可以看出，Ⅰ号、Ⅱ号、Ⅲ号、Ⅳ号、Ⅴ号、Ⅵ号磁异常中心在下延25 m时，其磁异常中心值增加速度快，说明延拓深度接近磁性地质体深度；而Ⅶ号磁异常中心在下延25 m时，其中心值增加速度快，且磁异常区及周边场源明显出现凌乱，定性分析已经接近磁性地质体场源上表。由此推断，引起Ⅶ号磁异常的磁性地质体顶板埋深在25 m左右。
从图 4可以看出，Ⅰ号、Ⅳ号磁异常中心在下延50 m时，磁异常区及周边场源明显出现凌乱，说明其已经接近磁性地质体场源上表，结合图 3推断，它们的磁性地质体顶板埋深在25~50 m之间；Ⅱ号、Ⅴ号、Ⅵ号磁异常中心在下延50 m时，磁异常区及周边场源开始出现凌乱，说明延拓深度接近磁性地质体深度；Ⅲ号磁异常中心在下延50 m时，其磁异常中心值增加速度快，说明延拓深度向磁性地质体深度靠近。
从图 5可以看出，Ⅰ号、Ⅳ号磁异常中心在下延75 m时，磁异常周围非常凌乱，证实延拓深度超过磁性地质体顶板；Ⅱ号、Ⅴ号、Ⅵ号磁异常中心在下延75 m时，磁异常区及周边场源出现明显凌乱，定性分析已经接近至磁性地质体场源上表。因此推断引起Ⅱ号、Ⅴ号、Ⅵ号磁异常的磁性地质体顶板埋深在75 m左右；Ⅲ号磁异常中心在下延75 m时，磁异常区及周边场源开始出现凌乱，说明延拓深度接近磁性地质体深度，推断引起Ⅲ号磁异常的磁性地质体顶板埋深在75 m以下。
磁异常面积下延综合表明，Ⅰ号、Ⅳ号磁性地质体顶板埋深在25~50 m之间；Ⅱ号、Ⅴ号、Ⅵ号磁性地质体顶板埋深在75 m左右；Ⅲ号磁性地质体顶板埋深在75 m以下；Ⅶ号磁性地质体顶板埋深在25 m左右。此外，地面磁测ΔT下延25 m、50 m、75 m等值线图显示，有一条NE向呈连续的串珠状的等值线，推测为一条断裂构造(TF1)。
4.3 精测剖面下延与2.5D反演拟合
通过对本区的7个异常8条(242线、246线、254线、274线、290线、292线、296线、336线)精测剖面进行曲线圆滑，高频滤波等数据处理，去掉了地表干扰和一些突跳点。利用RGIS2010对精测剖面进行向下延拓拟推断面处理和2.5D反演拟合，从下延拟推断面图和2.5D拟合结果看，8条精测剖面下延拟推与2.5D拟合结果基本一致(表 2)。
表 2 精测剖面下延拟推和2.5D反演拟合各地磁体空间分布统计结果Table 2. Statistical table of each magnetic bodies' spatial distribution by inferred profile with downward continuation and 2.5D inversion fitting线号 对应磁异常 推测磁性体顶板埋深/m 推测磁性体倾角/° 精测剖面下延拟推 2.5D反演拟合 精测剖面下延拟推 2.5D反演拟合 242 Ⅲ 118 62 75 78 246 Ⅳ 56 60 78 68 254 Ⅲ 60 55 75 75 274 Ⅱ 80 77 70 72 290 Ⅰ 37 42 75 69 292 Ⅵ 70 70 75 73 296 Ⅵ 60 80 75 70 336 Ⅶ 20 18 80 76 例如246线剖面，在对精测剖面进行向下延拓时，首先考虑到随着延拓深度的加大，在有意义的异常得到“放大”的同时，也把某些浅部的干扰或误差随着延拓深度的加大而放大，使延拓曲线发生剧烈的跳动，出现震荡现象，为克服延拓引起的震荡影响，在延拓前先对剖面异常进行“光滑”(即滤波处理)。向下延拓深度与精测剖面的取样点距(即采用数据的剖面长度)有直接关系，一般要求延拓的点距大于延拓的深度，本次取样点距为450 m。采用泊松积分公式进行向下延拓，结合向下延拓实践中出现的规律性，在向下延拓至100 m时，异常出现起伏跳动的震荡现象，就停止向下继续延拓，即到达延拓的深度。同时结合本次地面高精度磁测以大于等于1600 nT的异常确定为有意义的磁异常，推测磁性体顶部埋深56 m，倾角78°。在本剖面野外地质特征和磁异常的基础上，进行了2.5D反演拟合，剖面方位为150°，磁化强度为11271.3×10-2，模型密度为3.35 g/cm3，地磁倾角为为-26.2231°、磁偏角为-1.061°，反演结果磁铁矿体顶部埋深60 m，倾角68°。由此可见，下延拟推断面图和2.5D拟合结果基本一致(图 6、图 7)。
![]()
图 6 智利英格瓦塞中部矿区246线精测剖面曲线及下延拟推断面图Figure 6. Fine profile curve and inferred profile with downward continuation of the 246 surveying line of the central mining area of the Incaguasi iron deposit in Chile![]()
图 7 智利英格瓦塞中部矿区246线精测剖面2.5D反演拟合剖面图Figure 7. Fine profile from 2.5D inversion fitting of the 246 surveying line of central mining area of the Incaguasi iron deposit in Chile需特别说明的是，242、246、254、274、292、336号线的2.5D反演拟合磁异常剖面曲线推测为由2个到几个不等的磁性体引起，若要开展下一步地质工作，应首先对引起磁异常峰值的磁性地质体进行验证。
5. 中国IOCG型矿床勘查简论
自20世纪70年代，在探明的澳大利亚奥林匹克坝铜-铁-金-铀矿床中发现了富铁氧化物这一新特征，受到大家的关注。Hitzman等[23]通过研究澳大利亚奥林匹克坝矿床、加拿大韦尔内克山矿床、瑞典基鲁钠铁矿、美国密苏里东南铁矿及中国白云鄂博矿床，发现它们具有富含铁氧化物矿物组合的共同特征，于1992年将其统称为元古宙铁氧化物(铜-铀-金-稀土)矿床。因该类矿床巨大的潜在经济价值和理论研究意义，也逐渐受到了中国专家学者的关注，并对IOCG矿床进行了研究现状评述[2-3, 24-27]。聂凤军等[3]认为，IOCG型矿床主要产于大陆边缘俯冲带或克拉通盆地(裂谷、张裂带)，容矿围岩多为元古宙、中生代和新生代火山-沉积岩，成矿与A型、I型或磁铁矿系列侵入岩有关，矿体主要为角砾状、脉型、浸染状或块状矿石构成，矿石矿物富铁氧化物(磁铁矿和赤铁矿)而相对贫铁硫化物，围岩蚀变从下往上为钠化带、钾化带和硅化-绢云母化-赤铁矿化带，其矿床成因主要有岩浆流体说、地表/盆地流体说和变质流体说3类[3]。尽管中国对IOCG型矿床的研究刚刚起步，但是一些专家根据其地质特征，建议将长江中下游宁芜庐枞地区的玢岩铁矿、海南石碌铁矿、东准格尔地区的老山口和乔夏哈拉铁-铜-金矿床、内蒙古白云鄂博铁矿、四川拉拉矿床等矿床归为IOCG型矿床[2-3, 27]，展示了很好的资源前景。
随着传统硫化物铁型铜多金属矿床浅部资源的耗竭和深部勘查难度的增大，寻找以氧化铁为主体的铜金矿床(IOCG型矿床)则是一新的勘查方向。富铁氧化物和贫铁硫化物是IOCG型矿床的重要特征，对于勘查IOCG型铁矿而言，利用高精度磁测可以快速地圈定因富含磁铁矿引起的局部磁异常，这也是国内外勘查磁铁矿床最常见、有效的勘查方法。再利用RGIS软件对磁异常进行反演处理分析，推测和圈定矿体埋深和产出状态，能有效地减少中国此类矿床的后期勘查风险，具有很好的应用前景。
比如，中国海南石碌铁矿，其大地构造位于华南褶皱系五指山褶皱带的西段。区内地层主要为青白口系石碌群、震旦系石灰顶组，东南部分布有石炭系南好组、二叠系南龙组、峨查组-鹅顶组等。其中，石碌群中段岩性为条带状透辉透闪石化白云岩、条带状透辉透闪岩、铁质千枚岩等，是主要的赋铁矿围岩。矿区北部、西部、南部三面均被花岗岩侵入。区内褶皱、断裂发育，对赋矿地层和矿体起着重要的控制或改造作用。主矿床主要呈层状、似层状产出。铁矿石中，赤铁矿和石英分别是富铁矿石主要的矿石和脉石矿物，而贫铁矿石矿物主要为赤铁矿和磁铁矿，脉石矿物由石英、透辉石、透闪石、白云石等组成[26, 28-29]。智利英格瓦塞铁矿大地构造位于美洲大陆西缘安第斯褶皱山系南美褶皱带构造区内中生代格拉岛弧弧前盆地，矿区地层主要为中生界下白垩统大十二弦琴组(其岩性为安山岩)，铁矿体主要产于含磁铁辉石角闪石岩中，赋矿围岩基本为花岗闪长岩。矿体呈带状、层状，矿石矿物主要为磁铁矿，次为赤铁矿。虽然中国海南石碌铁矿与智利英格瓦塞铁矿在产状背景、赋矿地层、赋矿围岩、成因有很大差别，但是两者的矿石皆具有富铁氧化物的特征，与围岩的磁性参数差异较大，为高精度磁测奠定了很好的基础。例如，蒋立婷[30]通过高精度磁测在中国海南石碌铁矿矿区东南部如翁—牙加地区发现了7个磁异常，确定了2个铁多金属重点找矿靶区。
6. 结论
(1) 通过对智利英格瓦塞铁矿进行矿区地质、地球物理特征综合研究，利用RGIS软件对高精度磁测数据进行了面积上延、面积下延、精测剖面下延及2.5D反演拟合处理。从结果看，英格瓦塞中部矿区磁异常总规模为371585 m2，引起磁异常的磁性地质体整体呈NE向延伸、倾向NW、倾角较陡。除Ⅰ号和Ⅴ号磁异常规模较小外，其他几个磁异常中心找矿前景较大，有规模的磁异常区大多埋深在50~120 m之间。为本区下一步地质找矿方向提供了很好的依据。
(2) 虽然对IOCG型铁矿床的产出背景、矿床特征、成因机制还有待进一步认识，但是该类矿床富含铁氧化物和明显的地球物理异常却是其显著特征，这为高精度磁测奠定了很好的基础。利用RGIS软件对高精度磁测数据进行处理和解译，可获得磁性体相对可靠的空间产状和分布特征。因而，结合RGIS软件，利用高精度磁测对IOCG型铁矿床的找矿勘查是一种有效的手段，建议在中国大力推广应用。
致谢: 室内工作得到国家地质实验测试中心Re-Os同位素实验室李超博士和李欣尉老师的帮助与指导，两位审稿专家提出了宝贵的修改意见，让本文质量有了很大提高，在此一并表示衷心的感谢。 -
![]()
图 1 太行山北段地质简图(显示主要岩体和重要矿床分布，据参考文献[11]修改)
Figure 1. Simplified geological map of northern Taihang Mountains showing the distribution of major plutons and deposits
![]()
图 2 赤瓦屋钼铜矿地质简图(a)和0号勘探线地质剖面(b)
Figure 2. Sketch geological map(a)and section along No.0 exploration line(b)of Chiwawu Mo-Cu deposit
![]()
图版1
a、b.黄铜矿和辉钼矿赋存于石英细脉中; c、d.石英、黄铜矿、黄铁矿脉切穿斑状花岗闪长岩; e、f.黄铜矿和辉钼矿在石英脉中呈浸染状产出; g.浸染状辉钼矿，黄铜矿和磁黄铁矿共生; h.辉钼矿和黄铜矿共生; i.黄铜矿和黄铁矿、闪锌矿共生；Kfs—钾长石; Qz—石英; Ccp—黄铜矿; Py—黄铁矿; Mo—辉钼矿; Po—磁黄铁矿; Sp—闪锌矿
图版1.
![]()
图 3 赤瓦屋钼铜矿床Re-Os同位素模式年龄(a)及等时线年龄(b)
Figure 3. Molybdenite model ages(a)and Re-Os isochron age(b)of the Chiwawu Mo-Cu deposit
表 1 赤瓦屋钼铜矿床Re-Os同位素测试结果
Table 1 Molybdenite Re-Os isotope data of the Chiwawu Mo-Cu deposit
样品号 样重/g Re /10-6 普Os /10-9 187Re /10-9 187Os /10-9 模式年龄/Ma 测定值 2σ 测定值 2σ 测定值 2σ 测定值 2σ 测定值 2σ CWb31 0.03097 31.6 0.35 0.0062 0.0162 19831 222 43.20 0.26 130.6 2.1 CWb3 0.02012 37.1 0.29 0.0011 0.0602 23344 184 50.27 0.41 129.1 1.9 CWb9 0.02025 30.3 0.29 0.0011 0.0481 19067 182 41.07 0.25 129.2 1.9 CWb18 0.02018 26.2 0.21 0.0011 0.0476 16445 133 35.63 0.21 129.9 1.8 CWb24-1 0.02047 32.5 0.32 0.0011 0.0474 20409 200 44.16 0.30 129.7 2.0 CWk3 0.02004 27.2 0.21 0.0011 0.0141 17115 135 37.02 0.25 129.7 1.9 CWk6 0.02018 25.9 0.19 0.0011 0.0137 16278 119 35.19 0.24 129.6 1.8 
下载: 导出CSV
表 2 太行山北段赤瓦屋和麻棚岩体及相关矿床测年数据
Table 2 Chronology data of Chiwawu and Mapeng intrusions and their related ore deposits in the northern Taihang Mountains
矿床及相关岩体名称 年龄/Ma 测试方法 参考文献 赤瓦屋钼铜矿床 128.7±4.4 辉钼矿Re-Os等时线法 本文数据 秋树林钼矿 129.06±1.13 辉钼矿Re-Os模式年龄 [10] 阎家沟钼矿 126.7±1.1 辉钼矿Re-Os等时线法 [10] 银洞钼矿 127.6 辉钼矿Re-Os模式年龄 [10] 赤瓦屋岩体 134±1 锆石LA-ICP-MS U-Pb法 [11] 133±1 锆石LA-ICP-MS U-Pb法 [11] 131±1 锆石LA-ICP-MS U-Pb法 11] 128±1 锆石LA-ICP-MS U-Pb法 [11] 134.0±5.3 锆石SHRIMP U-Pb法 [12] 139.8±3.1 锆石SHRIMP U-Pb法 [12] 126.4±2.4 锆石LA-ICP-MS U-Pb法 [13] 130±1.0 锆石LA-ICP-MS U-Pb法 [14] 128.2±1.5 锆石LA-ICP-MS U-Pb法 [14] 麻棚岩体 131±2 锆石LA-ICP-MS U-Pb法 [9] 131±2 锆石LA-ICP-MS U-Pb法 [9] 125±3 锆石SHRIMP U-Pb法 [36] 125.4±2 锆石LA-ICP-MS U-Pb法 [13] 126.2±2 锆石LA-ICP-MS U-Pb法 [13] 128.3±1.6 锆石LA-ICP-MS U-Pb法 [37] 125.1±1.1 锆石LA-ICP-MS U-Pb法 [37] 124.3±1.3 锆石LA-ICP-MS U-Pb法 [37]
下载: 导出CSV
-
徐义刚. 太行山重力梯度带的形成与华北岩石圈减薄的时空差异性有关[J]. 地球科学, 2006, 31: 14-22. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200601003.htm 陈斌, 翟明国, 邵济安. 太行山北段中生代岩基的成因和意义: 主要和微量元素地球化学证据[J]. 中国科学(D辑), 2002, 32(11): 896-907. https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200211002.htm Chen B, Chen Z C, Jahn B M. Origin of mafic enclaves from the Taihang Mesozoic orogen, north China Craton[J]. Lithos, 2009, 110: 343-358. doi: 10.1016/j.lithos.2009.01.015
朱日祥, 陈凌, 吴福元, 等. 华北克拉通破坏的时间、范围与机制[J]. 中国科学(D辑), 2011, 41(5): 583-592. https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201105001.htm Mao J W, Pirajno F, Lehmann B, et al. Distribution of porphyry deposits in the Eurasian continent and their corresponding tectonic settings[J]. Journal of Asian Earth Sciences, 2014, 79: 576-584. doi: 10.1016/j.jseaes.2013.09.002
高永丰, 魏瑞华, 侯增谦, 等. 木吉村斑岩铜矿成矿作用: 华北克拉通中生代岩石圈减薄的响应[J]. 矿床地质, 2011, 30: 890-902. doi: 10.3969/j.issn.0258-7106.2011.05.010 者萌, 胡建中, 周伟, 等. 河北省安妥岭钼矿床地质特征及辉钼矿Re-Os同位素年龄[J]. 现代地质, 2014, 28: 339-347. doi: 10.3969/j.issn.1000-8527.2014.02.012 黄典豪, 杜安道, 吴澄宇, 等. 华北地台钼(铜)矿床成矿年代学研究——辉钼矿铼-锇年龄及其地质意义[J]. 矿床地质, 1996, 15(4): 365-373. https://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ604.008.htm Li S R, Santosh M, Zhang H F, et al. Inhomogeneous lithospheric thinning in the central North China Craton: Zircon U-Pb and S-He-Ar isotopic record from magmatism and metallogeny in the Taihang Mountains[J]. Gondwana Research, 2013, 23: 141-160. doi: 10.1016/j.gr.2012.02.006
孙文燕. 太行山中段银钼成矿系统结构及其深部找矿远景[D]. 中国地质大学(北京)博士学位论文, 2014. 李瑞玲, 段超, 陈志宽, 等. 太行山北段赤瓦屋铜钨矿化区花岗质岩石的锆石U-Pb年龄及其地质意义[J]. 中国地质, 2016, 43(5): 1761-1770. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201605024.htm 刘阳, 李程明, 穆一青, 等. 太行山北段赤瓦屋岩体锆石SHRIMP U-Pb年龄及其意义[J]. 地质与勘探, 2010, 46(3): 442-447. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201003010.htm 李林林, 韩宝福, 苗国均, 等. 太行山阜平杂岩中麻棚赤瓦屋岩体的时代、侵位深度及构造意义[J]. 岩石矿物学杂志, 2012, 31(3): 289-306. doi: 10.3969/j.issn.1000-6524.2012.03.001 He X F, Santosh M. Crustal recycling through intraplate magmatism: Evidence from the Trans-North China Orogen[J]. Journal of Asian Earth Sciences, 2014, 95: 147-163. doi: 10.1016/j.jseaes.2014.02.011
Zhao G C, Cawood P A, Li S Z, et al. Amalgamation of the North China Craton: Key issues and discussion[J]. Precambrian Research, 2012, 222/223: 55-76. doi: 10.1016/j.precamres.2012.09.016
Zhao G C, Wilde S A, Cawood P A, et al. Archean blocks and their boundaries in the North China Craton: lithological, geochemical, structural and P-T path constraints and tectonic evolution[J]. Precambrian Research, 2001, 107: 45-73. doi: 10.1016/S0301-9268(00)00154-6
段超, 毛景文, 谢桂青, 等. 太行山北段木吉村髫髻山组安山岩锆石U-Pb年龄和Hf同位素特征及其对区域成岩成矿规律的指示[J]. 地质学报, 2016, 90(2): 250-266. doi: 10.3969/j.issn.0001-5717.2016.02.005 邓晋福, 莫宣学, 赵海玲, 等. 中国东部岩石圈根/去根作用与大陆"活化"——东亚型大陆动力学模式研究计划[J]. 现代地质, 1994, 8(3): 349-356. https://www.cnki.com.cn/Article/CJFDTOTAL-XDDZ403.017.htm 河北省地质矿产局. 河北省北京市天津市区域地质志[M]. 北京: 地质出版社, 1989: 125-135. 喻学惠, 任建业, 张俊霞. 太行山中段铜-金成矿条件及找矿方向[M]. 北京: 地质出版社, 1996: 11-30. 王雪飞. 综合物化探方法在河北保定赤瓦屋铜多金属矿区勘查中的应用研究[D]. 吉林大学硕士学位论文, 2020. 杜安道, 何红寥, 殷宁万, 等. 辉相矿的铼-锇同位素地质年龄测定方法研究[J]. 地质学报, 1994, 68(4): 339-347. doi: 10.3321/j.issn:0001-5717.1994.04.005 杜安道, 屈文俊, 李超, 等. 铼-锇同位素定年方法及分析测试技术的进展[J]. 岩矿测试, 2009, 28(3): 288-304. doi: 10.3969/j.issn.0254-5357.2009.03.019 杜安道, 屈文俊, 王登红, 等. 铼-锇法及其在矿床学中的应用[M]. 北京: 地质出版社, 2012: 11-80. 屈文俊, 杜安道. 高温密闭溶样电感耦合等离子体质谱准确测定辉钼矿铼-锇地质年龄[J]. 岩矿测试, 2003, 22(4): 254-262. doi: 10.3969/j.issn.0254-5357.2003.04.003 屈文俊, 杜安道, 任静. 过氧化氢在黄铁矿的溶解过程中对铼-锇信号强度及年龄的影响研究[J]. 分析化学, 2008, 8(2): 223-226. https://www.cnki.com.cn/Article/CJFDTOTAL-FXHX200802018.htm 屈文俊, 杜安道, 李超. 金川铜镍硫化物样品中锇同位素比值的高精度分析[J]. 岩矿测试, 2009, 28(3): 219-222. doi: 10.3969/j.issn.0254-5357.2009.03.005 李超, 杜安道, 屈文俊. 大颗粒辉钼矿Re-Os同位素失耦现象及187Os迁移模式研究[J]. 矿床地质, 2009, 28(5): 707-712. doi: 10.3969/j.issn.0258-7106.2009.05.016 李超, 屈文俊, 杜安道. 铼-锇同位素定年法中丙酮萃取铼的系统研究[J]. 岩矿测试, 2009, 28(3): 233-238. doi: 10.3969/j.issn.0254-5357.2009.03.008 李超, 屈文俊, 周利敏, 等. Carius管直接蒸馏快速分离锇方法研究[J]. 岩矿测试, 2010, 29(1): 14-16. https://www.cnki.com.cn/Article/CJFDTOTAL-YKCS201203007.htm Ludwig K. Isoplot/Ex, version 2.0: A geochronogical toolkit for Micros of Excel[C]//Berkeley Geochronology Center, Berkeley, CA, 2009.
Smoliar M I, Walker R J, Morgan J W. Re-Os ages of groupⅡA, ⅢA, ⅣA and ⅥB iron meteorites[J]. Science, 1996, 271: 1099-1109. doi: 10.1126/science.271.5252.1099
Stein H J, Markey R J, Morgan J W, et al. The remarkable Re-Os chronometer in molybdenite: How and why it works[J]. Terra Nova, 2001, 13(6): 479-486. doi: 10.1046/j.1365-3121.2001.00395.x
Foster J G, Lambert D D, Frick L R, et al. Re-Os isotopic evidence for genesis of Archaean nickel ores from uncontaminated komatiites[J]. Nature, 1996, 382(6593): 703-706.
李文昌, 余海军, 尹光侯, 等. 滇西北铜厂沟钼多金属矿床辉钼矿Re-Os同位素年龄及其成矿环境[J]. 矿床地质, 2012, 31(2): 282-292. https://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ201202010.htm 刘阳, 李程明, 郑杰, 等. 太行山北段麻棚岩体锆石SHRIMP U-Pb年龄及其意义[J]. 地质与勘探, 2010, 46: 622-627. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201004007.htm Li Q, Santosh M, Li S R, et al. The formation and rejuvenation of continental crust in the central North China Craton: Evidence from zircon U-Pb geochronology and Hf isotope[J]. Journal of Asian Earth Sciences, 2014, 95: 17-32. http://www.onacademic.com/detail/journal_1000036245042610_7ca5.html
Li S R, Santosh M. Metallogeny and craton destruction: Records from the North China Craton[J]. Ore Geology Reviews, 2014, 56: 376-414. http://www.sciencedirect.com/science/article/pii/s0169136813000863
毛景文, 谢桂青, 张作衡, 等. 中国北方中生代大规模成矿作用的期次及其地球动力学背景[J]. 岩石学报, 2005, 31(1): 169-188. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200501018.htm Mao J W, Xie G Q, Duan C, et al. A tectono-genetic model for porphyry-skarn-stratabound Cu-Au-Mo-Fe and magnetite-apatite deposits along the Middle-Lower Yangtze River Valley, Eastern China[J]. Ore Geology Reviews, 2011, 43: 294-314. http://www.ingentaconnect.com/content/el/01691368/2011/00000043/00000001/art00024
毛景文, 段超, 刘佳林, 等. 陆相火山-侵入岩有关的铁多金属矿成矿作用及矿床模型——以长江中下游为例[J]. 岩石学报, 2012, 28(1): 1-14. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201201003.htm 周涛发, 范裕, 袁峰, 等. 长江中下游成矿带火山岩盆地的成岩成矿作用[J]. 地质学报, 2011, 85(5): 712-730. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE201105010.htm Mao J W, Xie G Q, Bierlein F, et al. Tectonic implications from Re-Os dating of Mesozoic molybdenum deposits in the East Qinling-Dabie orogenic belt[J]. Geochimica Cosmochimica Acta, 2008, 72: 4607-4626. http://www.sciencedirect.com/science/article/pii/S0016703708004043
毛景文, 罗茂澄, 谢桂青, 等. 斑岩铜矿床的基本特征和研究勘查新进展[J]. 地质学报, 2014, 88(12): 2153-2175. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE201412002.htm -
期刊类型引用(2)
1. 安国强，鲁宝亮，高新宇，朱武，李柏森. 边缘特征和深度加权约束的重力三维相关成像反演. 物探与化探. 2024(01): 113-124 . 
百度学术
2. 闫顺尚，王玲，董方营，刘晓，吕晓磊. 基于电-磁法联合勘探的火烧区范围及富水性的预测. 煤田地质与勘探. 2022(02): 132-139 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 2283
- HTML全文浏览量: 520
- PDF下载量: 1759
- 被引次数: 5
