第一作者简介: 张家嘉(1987—),男,高级工程师,主要从事地质矿产勘查工作。Email: 125313656@qq.com。
凤阳江山金矿是皖东地区规模最大的金矿床,近年来在其周边实施了多个勘查项目,但找矿效果不佳。通过开展重力扫面、重力剖面、音频大地电磁测深测量(audio magnetotelluric,AMT)等物探方法,总结区域重磁场特征,认为研究区存在隐伏构造岩浆岩带; 通过对比已知矿床的物探异常特征,认为AMT低阻异常带与含金层间破碎带具有良好的对应关系; 在总结分析成矿地质条件和成矿规律的基础上,认为鹿塘断裂是区内主要控矿构造; 通过重力反演,认为深部隐伏岩体顶界面埋深大部分约1 100 m; 基于“四步”式三维成矿定量预测方法体系,开展凤阳江山地区三维成矿预测研究,圈定了5个找矿靶区,为下步勘查工作提供依据。
Fengyang Jiangshan gold mine is the largest gold deposit in eastern Anhui Province. In recent years, multiple exploration projects have been implemented in its surrounding areas without any good prospecting results. Through the use of geophysical methods, such as gravity scanning, gravity magnetic profiling, and audio magnetotelluric depth sounding (AMT), the authors summarized the characteristics of regional gravity and magnetic fields and considered that there are hidden tectonic magmatic rocks in the study area. The AMT low resistance anomaly zone is believed to be corresponded well with the gold bearing interlayer fracture zone, after the analysis of the geophysical characteristics of known mineral deposits. Besides, Lutang fault is the main ore controlling structure in the area, on basis of the summary and analysis of the geological conditions and laws of mineralization. Through gravity inversion, it is believed that the top interface burial depth of hidden rock mass in the deep part is mostly around 1 100 m. The 3D mineralization prediction research in Jiangshan of Fengyang area was conducted by 3D mineralization quantitative prediction method system with 4 steps, and 5 prospecting target areas were identified, providing the basis for the exploration work in the future.
凤阳地区位于蚌埠隆起东段, 蚌埠隆起东段西起怀远县, 经蚌埠市、凤阳县, 止于五河县, 该区为由一系列近EW向和近SN向断裂切割形成的断块格局, 主体呈EW走向, 并具西窄东宽的特点[1, 2]。蚌埠隆起东段一直是安徽省金多金属矿勘查的重要远景区和找矿靶区, 目前已发现的矿种有金、铅、铁、钛、钼、铜、晶质石墨、重晶石等, 共有矿床(点)20余处, 其中金矿床有江山金矿、毛山金矿、中家山铅锌(金银)矿和上王庄金矿, 金总资源量近20 t, 但只有江山金矿达到中型规模, 其他均为小型矿床(点)。
近年来, 前人在凤阳地区开展了大比例尺物探工作, 如高精度磁法、激电中梯、激电测深和可控源音频大地电磁测深等方法, 但由于矿权分散, 使得完整的成矿区带被分割, 各矿业权人开展地质找矿时均受其矿业权边界限制, 受勘查范围不超过探矿权或采矿权范围的限制, 存在勘查范围小、勘查视野窄的弊端, 从而不能从更宏观的成矿地质背景、控矿条件找矿。这是造成近年来凤阳地区未能取得新的找矿成果的主要原因之一。
在新一轮找矿突破战略行动中, 为进一步明确凤阳地区找矿方向, 安徽省开展实施了“ 皖东五河— 凤阳矿集区三维综合探测与深部金铅锌多金属成矿预测” 项目, 实施了1∶ 5万高精度重力测量、1∶ 1万重力剖面测量、音频大地电磁测深测量(audio magneto-telluric, AMT)等物探方法。凤阳江山金矿最早为地表开采的小型铅锌矿, 随后在矿床东部勘查发现金矿化线索, 经过多年勘查, 截至2016年勘探累计探明金资源量约17 t, 是蚌埠隆起区东段目前已探明的最大金矿床[3], 作为典型的蚀变岩型金矿, 在矿床深部各矿体是否联通及深部找矿潜力等方面还需要进一步的预测和研究。
三维地质建模是实现深部矿产勘查突破的重要途径, 也是开展三维成矿预测工作的重要基础。通过对控制地质体和地质结构进行三维建模, 可以更加精细地刻画控矿要素之间的空间、成因和演化关系, 帮助理解成矿系统, 开展深部找矿预测。因此, 本文通过对比江山金矿矿床地质特征, 总结分析成矿地质条件和成矿规律, 圈定找矿靶区, 为凤阳地区下步找矿方向提供参考。
蚌埠隆起区位于华北陆块东南缘, 北以利辛— 五河断裂(F1)为界与淮北断褶带相邻; 南以刘府断裂(F2)为界, 与淮南褶断带相接, 东至郯庐断裂带, 向西经过夏邑— 固始断裂延伸至河南省(图1)。区域上还发育一系列NE— NNE向断裂, 以往研究认为这些断裂是郯庐断裂带的次级断裂, 也是本区的主要控矿构造[4]。蚌埠隆起区大部分被第四系覆盖, 基底岩系为太古代五河群变质岩, 包括表壳岩系、变质镁铁质岩系和变质变形花岗质侵入体, 《安徽省区域地质志》[1]将五河群由老到新分为下亚群和上亚群, 下亚群包括西堌堆组、庄子里组和峰山李组, 上亚群包括小张庄组和殷家涧组。
蚌埠隆起区岩浆活动强烈, 区内产出不同时代的花岗质侵入岩体。其中磨盘山和庄子里岩体形成于古元古代[5], 其岩性主要为变钾长花岗岩[6]; 其余侵入岩均为中生代岩浆活动的产物, 大致分为3期: 晚侏罗世岩体(约167~162 Ma)如荆山岩体和涂山岩体, 其岩性以混合花岗岩为主; 早白垩世早期岩体(131~127 Ma)如女山岩体、西芦山岩体等, 其岩性以花岗闪长岩为主; 早白垩世晚期岩体(117~110 Ma)如曹山岩体、锥子山岩体等, 其岩性以二长花岗岩为主[7, 8, 9, 10]。
凤阳江山金矿位于安徽省凤阳县城东南部, 矿区地表多被第四系覆盖, 根据钻孔揭露, 其下覆地层主要为庄子里组, 岩性主要为斜长角闪岩、黑云斜长片麻岩、浅粒岩和大理岩, 其次为峰山李组, 岩性主要为黑云斜长片麻岩、角闪斜长片麻岩。矿区位于近EW向官沟— 下纪家断裂和NNE向临淮关— 亮岗断裂的交汇部位。官沟— 下纪家断裂倾向南, 倾角约35° , 具有压性特征, 在矿区内长约3 000 m, 宽110~113 m; 临淮关— 亮岗断裂在矿区内长约4 500 m, 断裂早期为塑性流动变形, 表现为糜棱岩化、片理化, 后期为构造角砾岩化、碎裂岩化。矿区地表大面积出露正长岩, 钻孔揭露还存在石英闪长岩和闪长玢岩[11]。
金矿体主要沿着断层及次级断层分布, 严格受层间破碎带控制(图2)。矿体呈层状、似层状、透镜状产出, 走向NNE, 倾向南东, 产状平缓, 一般为0° ~25° 。赋矿围岩主要为碎裂状硅化大理岩、碎裂岩。主要的矿石矿物为黄铁矿, 脉石矿物主要为石英、绢云母和少量的碳酸盐矿物。矿石的构造主要为块状构造、细脉浸染状构造和浸染状构造等, 矿石的结构主要为自形— 半自形结构[12]。
![]() | 图2 江山矿区33勘探线地质剖面示意图(据参考文献[12]修改)Fig.2 Schematic diagram for geological profile of exploration line 33 in Jiangshan mining area (modified after reference [12]) |
研究收集了凤阳地区主要岩性的物性参数, 并在江山— 屈家湾地区钻孔中实测了部分岩性的物性参数。研究区主要地层为五河群庄子里组, 其岩性主要为黑云斜长片麻岩、斜长角闪片麻岩和斜长角闪岩。从表1和表2对比可发现, 研究区内3种主要岩性的密度参数明显高于区域上的数值。
![]() | 表1 研究区岩性物性参数统计 Tab.1 Statistics of regional lithological and physical property parameters in the study area |
![]() | 表2 研究区主要岩性物性参数统计 Tab.2 Statistics of lithological and physical properties parameters in the study area |
研究区总体位于NW向蚂蚁山— 总铺重力低值带与NNE向陈家户— 头铺重力低值带交汇区, 江山金矿处于NE向庞家户— 屈家湾重力高东凸部位, 呈鼻状。研究区及周边的第四系厚度约10~15 m, 埋深较为稳定; 重力场总体表现为重力低, 其北东方向后毛— 庙家、南西方向庞家户一带的重力高, 经钻探揭露为五河群变质岩系, 而区内主要岩性的密度参数又明显高于区域上的数值, 因此重力低显然不是由第四系引起, 而是深部酸性岩体的反映。推测NW向蚂蚁山— 总铺重力低值带是怀远— 蚌埠— 凤阳一线花岗质岩浆(位置见图1)向东南方向侵位所致, 而庞家户— 江山NE向的重力高则是五河岩群变质地层残留体的表征[13], 该NE向布格重力高有独立圈闭趋势, 并在江山金矿一带等值线向东凸出呈舒缓状分布(图3)。
为了进一步推断隐伏岩体的深度, 本次工作在1∶ 5 万布格重力异常图上选取了2条NNE剖面(图4), 利用反演软件对深部隐伏岩体开展了进一步反演。1号反演剖面显示隐伏岩体与围岩接触顶界面深度为900~1 900 m; 2号反演剖面隐伏岩体与围岩接触顶界面深度为 650~2 800 m。隐伏岩体顶界面在屈家湾— 江山一带埋深相对较深, 其他地方深度均在 1 100 m 以浅。
区域航磁异常特征总体表现为西北部和东北部为正磁场背景, 其上叠加了走向和形态不同的局部异常, 大多局部异常在其北侧伴生有明显的负异常(图5)。西北部蚂蚁山正磁异常呈近EW走向, 位于蚂蚁山重力低值区内; 该磁异常向东南方向延伸部分位于重力低的北侧梯级带上, 形成重梯磁高组合异常; 在蚂蚁山庄子里组大理岩接触带南北两侧, 各有一个沉积变质-热液改造型铁矿点, 重力低由深部隐伏中酸性岩体侵入引起。东北部临淮关— 杨家圩一带磁异常呈近EW向高低相间排列, 磁高是西堌堆组的反映, 而磁低则是低级次EW向断裂的表征; 后毛— 杨家圩一带发育十余个金矿点, 物探异常特征表现为重力高值区和磁异常边部。该区重力高相对完整反映了西堌堆组地层巨厚, 因此缺少足够的热源条件, 形成规模矿床可能性较小。
中南部航磁总体为平缓负磁背景区, 其上叠加了不同形态、不同走向和不同强度的局部磁高异常。如呈串珠状呈NE向展布的庞家户— 屈家湾磁异常, 磁化极异常极大值仅50 nT, 在其异常北东边部分布有大王府铅锌矿和江山金矿, 西侧边部分布有屈家湾金矿。
本次工作实施了AMT扫面工作, 线距4 km, 点距 500 m, 选取其中3条测线, 即102线、104线和106线进行初步解释, 3条测线方向均为SN向(图3)。
如前文所述, 江山金矿床主要受近水平的层间破碎带控制, 从图6(104线和106线)可以看出, AMT低阻异常带和含矿层间破碎带对应关系好, 且附近3条AMT剖面均有类似的特征, 表明浅部低阻异常带走向较为稳定, 其深度自西向东有逐渐变浅的趋势, 其平面的位置如图3所示, 整体走向为近EW— NNE向, 江山金矿和屈家湾金矿在300~350 m附近均发现有含金层间破碎带(图7)。因此, 利用AMT低阻异常带圈定找矿靶区具有较好的效果, 推测的低阻异常带上具有寻找江山式金矿的潜力。
本文总体基于“ 四步式” 三维成矿定量预测方法体系对凤阳地区开展三维成矿预测研究[14]。由三维地质数据库构建、三维地质建模、三维预测信息深度挖掘和三维数据融合及预测评价4个主要步骤组成。该方法能够有效融合多元、多维的地质数据并挖掘其中的特征和联系, 预测深部找矿靶区。
区内金矿体沿着断裂及次级断裂分布, 严格受层间破碎带控制。因此, 断裂是该区域内最重要的控矿条件。本文基于最新AMT资料成果, 通过融合江山金矿床最新矿产勘查数据资料, 利用三维隐式地质建模方法构建了区内断裂模型。然后结合三维空间分析和三维成矿过程数值模拟方法, 深入开展三维预测信息的挖掘。首先, 计算区内断裂三维距离场, 三维距离场分析能够反映控矿要素的影响范围, 计算出各个块体与控矿要素之间的距离; 其次, 由于有利倾角的断裂可以形成较为稳定的场所, 有助于矿体的形成和富集[15], 故对模型中的断裂倾角进行了计算。研究表明, 断裂在交汇部分往往是成矿最有利的部位[16], 提取断裂交点作为预测要素, 计算了其与研究区各块体的距离。此外, 断裂构造发育较密集区域也是矿体赋存的有利空间, 故提取了断裂有利密集区分析局部断裂发育位置, 并计算了三维距离场。根据地球物理数据推断的隐伏岩体, 作为控矿要素之一, 进行了三维距离场计算。基于三维成矿过程数值模拟结果, 提取体应变增量和剪切应变增量的高值区作为与成矿空间紧密相关的预测要素, 并计算它们的距离场。
基于上面构建完成的三维成矿预测信息集, 以江山金矿床主要矿体作为数据驱动端元, 采用隶属于机器学习范畴的 Logistic 回归方法对三维预测信息进行数据融合, 该方法能有效分析和度量不同三维控矿要素与已知矿体的相关性, 进而对整个模型的成矿概率进行计算[17]。逻辑回归方法的计算公式[18]为
式中: φ (d)为成矿概率; xi为第i种控矿要素; α 是一个常数; β i为由最大似然估计方法求得的回归因子。
本次研究选取了成矿概率0.01作为阈值, 提取了Logistic方法预测结果中的高成矿潜力单元, 并圈定了5处新的找矿靶区, 以供为进一步勘探提供方向和依据(图8(b))。
胶东金成矿区位于郯庐断裂带的东侧, 主要受5条断裂构造控制, 分别为三山岛断裂带、焦家断裂带、招平断裂带、栖霞断裂带和牟乳断裂带, 主要大、中型矿床均位于主干断裂附近, 研究认为这些主要的控矿断裂均为郯庐断裂带的次级断裂。主控矿构造是胶东金矿区最主要控矿因素, 其次大规模的岩浆活动特别是玲珑花岗岩为金矿的形成提供了热源及成矿物质来源[19]。蚌埠— 凤阳地区位于郯庐断裂带西侧, 与胶东地区具有相似的成矿地质背景, 金多金属矿床主要受郯庐断裂带及其次级断裂构造控制[20]。因此, 寻找主控矿构造及与成矿有关的岩体是本区实现找矿突破的重中之重。
研究区断裂构造发育, 主要有NE向、NNE向和NW向3组。其中NE向的鹿塘断裂(F2)规模较大, 布格重力异常在断裂中段表现为重力低, 其余部位则是重力梯级带(图3); 区域航磁异常错位明显, 在不同上延高度的重磁图中均有该重磁场特征出现[21]; 从AMT异常图上推测鹿塘断裂倾向NW(图6), 区内长约32 km。从图3可看出, 区内已发现的金多金属矿床点多位于鹿塘断裂北西侧, 尤其是北西侧的江山金矿构造破碎蚀变范围大, 金矿化明显, 矿体垂向上有向深部延伸的趋势, 今年新发现的屈家湾金矿体同样位于构造破碎带内。鹿塘断裂可能是凤阳地区主要的控矿构造。
蚂蚁山以东地区为蚌埠重力低, 向西一直延伸至怀远, 对应了怀远— 蚌埠一线的花岗质岩浆岩出露区(图1)。前文分析表明, 研究区存在隐伏岩体, 推测是蚌埠— 怀远花岗质岩体向南东方向延伸。
研究区由不同形态的局部重力高和局部重力低组成一个不完整的重力低值带, 主要为低背景负磁场, 相对磁力高形成的局部异常以NE向为主, 与区域性磁异常方向不一致。上述重磁场特征应是地下地质现象的客观反映, 表明该区岩浆和构造热液活动比较强烈, 各个小型局部重力低异常可能与侵入体有关, 重力低值带主要由中酸性岩浆侵入和断裂构造破坏结晶基底形成的。这就为形成矿产提供了足够的热源条件, 使地层的成矿元素活化转移、富集而形成矿体, 江山金矿以及多个金多金属矿床点便是例证。
(1)研究区重磁场异常复杂, 是地下复杂地质环境的反映, 区内存在隐伏的构造岩浆岩带。
(2)深部的岩浆活动为金多金属矿床提供了良好热源条件, 鹿塘断裂(F2)是区内主要的控矿断裂, 含矿热液沿着主控矿构造向上运移并在浅部低缓的破碎带处富集成矿。
(3)研究区还存在多个与江山金矿相似的物探异常, 浅部具有寻找江山式金矿的潜力。深部隐伏岩体与变质基底的接触部位可能是金矿重要的赋矿部位, 岩体内部是否有寻找斑岩型金矿的可能也值得探索。
(4)通过初步的三维成矿预测, 认为徐家湾— 二铺一带找矿潜力最大, 其次为总铺— 高陈家一带, 最后为潘家湾。
(责任编辑: 常艳)
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|