基于MRAS证据权重模型的河南老湾地区金矿成矿预测
Metallogenic prediction of gold deposits in Laowan area, Henan Province using the weight of evidence model and MRAS
通讯作者: 陈建立(1966-),男,河南长葛人,教授级高级工程师。Email:842186572@qq.com
责任编辑: 蒋实
收稿日期: 2021-06-11 修回日期: 2021-08-26
基金资助: |
|
Received: 2021-06-11 Revised: 2021-08-26
作者简介 About authors
魏从玲(1979- ),女,河南南阳人,硕士,高级工程师,主要从事地质勘查及研究工作。Email:
河南省老湾地区是桐柏—大别有色贵金属成矿带上重要的金成矿区。在综合分析该区成矿因素的基础上,提取构置了地层、韧性剪切域、脆性断束、岩浆岩、与Au有关的蚀变带和Au、As、Sb、Ag水系沉积物单元素异常等9个预测变量,基于MRAS软件,采用证据权重模型对该区金矿进行了成矿远景预测工作。本次采用0.10 km×0.10 km的网格单元进行预测,经过各预测变量的先验概率计算、权重值统计、独立性检验和后验概率的计算等过程,划分出A、B、C三个级别的成矿远景区块,根据预测网格分布特点及该区金矿地质特征,圈出了Ⅰ级找矿远景区4处,Ⅱ级找矿远景区4处。综合分析认为:已知金矿床或矿点多数落入A级或B级成矿远景区块,C级成矿远景区块主要是在A级或B级成矿远景区块的周边,与成矿区带从中心向边缘渐变的特征相符合;圈出的Ⅰ级找矿远景区内均已发现金矿床(点),并且后验概率平均值较高,具有较好的金矿找矿潜力。
关键词:
The Laowan area in Henan Province is an important gold minerogenetic area in the Tongbai-Dabie noble nonferrous metal metallogenic belt. Based on the comprehensive analysis of the metallogenic factors of the area, the authors extracted and established nine predictive variables, i.e., strata, ductile shear zones, brittle fault zones, magmatic rocks, gold-related alteration zones, and single-element anomalies of gold, arsenic, antimony, and sliver in stream sediments. Then, the authors predicted metallogenic prospect areas of gold in this area using the weight of evidence model and the MRAS software. Grid cells with a size of 0.10 km×0.10 km were used for prediction, and they were divided into classes A, B, and C in terms of prospecting prediction through the prior probability calculation, weight statistics, independence tests, and posterior probability calculation of each prediction variable. Based on the distribution characteristics of various classes of prediction grid cells and the geological characteristics of gold deposits in the study area, this study delineated four first- and four second-order prospecting target areas. According to comprehensive analysis, most of the known gold deposits or gold ore occurrences fall into the prediction blocks of classes A and B and the prediction blocks of class C lie around the prediction blocks of classes A or B. This analytical result is consistent with the gradual changes of the metallogenic belt from the center to the edges. Gold deposits or ore occurrences have been discovered in all the delineated first-order prospecting target areas, and the posterior probabilities have a high average. Therefore, first-order prospecting target areas have great gold prospecting potential.
Keywords:
本文引用格式
魏从玲, 陈建立, 郭鹏.
WEI Cong-Ling, CHEN Jian-Li, GUO Peng.
0 引言
河南省桐柏县老湾地区是河南省一处重要的黄金产地,处于秦岭褶皱系东段,桐柏—大别有色贵金属成矿带上。多年来,众多勘查、科研单位先后在此做了大量工作,对区内金矿床地质特征、金矿物的赋存状态和成矿物质来源、成岩成矿年代、成矿因素和控矿条件、构造体系演化等做了一系列的研究,并取得了诸多成果[11⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓-25]。目前,老湾地区已发现了10多处金矿床(点),并已提交一处超大型金矿床;对区内金矿床的地质特征以及控矿因素等也有了较为全面的认识,构建了区内金矿找矿模型,并基于找矿模型对区内的找矿前景进行了评价和预测。本次研究通过分析老湾地区的控矿因素,收集老湾地区地质、化探等资料,首次利用MRAS软件中的证据权重模型,采用定量化的方法,对老湾地区的金矿资源潜力进行远景预测和评价。
1 区域地质背景
老湾地区金成矿带位于桐柏北部地区,处于北为松扒断裂、南为老湾断裂、西为泌阳凹陷、东为吴城盆地,呈NWW向展布的一条狭长地带中(图1)。
图1
区内主要赋矿地层是中元古界龟山岩组, 主要岩性为二云石英片岩、斜长角闪片岩、斜长角闪岩和变质石英岩等。它是一套具有较强构造混杂作用的带状无序变质地层,呈层状、似层状或透镜状,整体呈NWW—SEE向平行于区域构造线展布,北侧以松扒断裂与秦岭岩群呈断层接触,南侧以老湾断裂与肖家庙岩组呈断层接触。
区内的金矿主要赋存在NWW向展布的韧性剪切带中,区内地层经长期的多期次变质、变形等作用,形成了许许多多的构造岩片及混杂岩带,除控制矿带南北两侧边界的老湾、松扒两个区域性断裂构造外,带内低序次的剪切、走滑构造也极为发育,岩层中的片理、劈理、糜棱构造十分发育,加上多期次的成矿作用,使区内地质构造显得十分复杂。
2 控矿因素分析
控矿因素和成矿条件研究是矿产预测的重要基础。老湾地区金成矿带经历了华北板块与扬子板块分离、汇聚、俯冲和碰撞的演变过程,形成了以韧性构造变形为主,脆韧性构造变形叠加,并伴有岩浆侵入的复杂构造变形带。成矿作用伴随着构造的演化发生和发展,形成了具有韧性剪切变质热液和岩浆热液叠加的复合型矿化类型;矿化受地层、构造、岩浆岩等多因素控制[19]。
1) 地层:区内所有的金矿床和金矿点均产于老湾韧性剪切带所夹持的龟山岩组地层中。而龟山岩组地层金平均含量11.17×10-9,是地壳丰度值的3倍多,龟山岩组既是金矿的主要赋矿围岩,也是成矿物质的主要来源,具有矿源层特点。
2) 构造:区内构造活动贯穿于整个成矿过程,是成矿和控矿的重要因素。不同层次、不同类型的韧性和脆性断裂构造控制了区内金矿脉的分布和产出特征。其中南部的老湾断裂和北侧的松扒断裂规模大、切割深,具有多期活动特点,是重要的控岩控矿构造,两条断裂构造的形态和产状控制了区内矿脉的分布范围、展布方向和形态特征。
区内断裂构造的多期次活动和走滑运动造成赋矿地层龟山岩组早期发生深层次的韧性剪切变形,伴随着变质作用促使龟山岩组地层中成矿物质活化,金成矿物质初步富集;晚期发生浅层次的脆—韧性剪切变形为含金成矿流体提供了运移通道和容矿构造空间;主断裂之间的一系列NWW向次级断裂为岩脉侵位和热液矿脉的就位提供了充足的空间,对金矿成矿作用控制效应更为直接和明显,控制了矿床和矿点的分布范围及产出位置。
4) 围岩蚀变:区内围岩蚀变主要有硅化、绢云母化、绿泥石化、碳酸盐化,其中由硅化生成的含金属硫化物石英脉与金矿化的关系最为密切,往往形成金矿体。绢云母化由白云母蚀变而来,代表着热液改造后的结果,往往与金矿化关系密切。绿帘石化常常伴随有硅化以及绢云母化。碳酸盐化多发生在成矿晚期,标志着热液成矿作用的结束,亦是找矿的良好标志。
5) 化探异常:区内矿致异常元素组合复杂,Au-As-Sb-Ag是区内的特征元素组合,以Au、As、Sb、Ag为主的多金属元素组合的化探异常往往是工业矿体存在的有利地段,有时伴生Cu、Pb、Zn组成的综合异常。
3 MRAS证据权重模型
MRAS软件的空间数据库是以矢量数据结构为主体形成的,可根据空间数据库对成矿信息进行提取分析,并对多元空间成矿信息进行资源潜力的综合评价[2]。
证据权重模型进行成矿预测有5个基本步骤:
1) 通过控矿因素和成矿条件的研究,确立地质因素与矿化的关系,提取构置预测变量。
2) 计算各预测变量的先验概率,研究预测变量存在时,矿床(点)出现和不出现的条件概率,以及预测变量不存在时,矿床(点)出现和不出现的条件概率。
3) 对各预测变量和矿床(点)做空间分析处理,计算预测变量的权重值W+、W-和C值,对预测变量进行优化选择。
4) 对优选出的预测变量进行条件独立性检验,进一步对预测变量进行优化,剔除相关性过高的预测变量。
5) 计算各预测单元的后验概率值,生成成矿预测远景潜力图。
4 老湾地区金矿成矿预测
4.1 预测变量的分析与提取
根据老湾金矿带的控矿因素和成矿条件的分析研究,提取以下控矿因素作为预测变量(图2)。
图2
图2
老湾地区金矿成矿预测变量示意
Fig.2
Variable map of metallogenic prediction in Laowan Gold Belt
1) 地层:区内的金矿全部分布于龟山岩组地层中,本次成矿预测将龟山岩组(1:5万)作为地层变量(图2a)。
5) 与Au有关的蚀变带:区内的围岩蚀变伴随整个成矿过程,蚀变带是找矿的良好标志。因此选择与Au有关的蚀变带(1:1万)作为金矿预测的变量(图2e)。
4.2 先验概率的计算统计
表1 研究区各预测变量先验概率统计结果
Table 1
预测变量 | P1 | P2 | P3 | P4 |
---|---|---|---|---|
地层 | 1.000 000 | 0.042 279 | 0.000 000 | 0.957 720 |
韧性剪切域 | 0.857 142 | 0.022 780 | 0.142 857 | 0.977 219 |
脆性断束50 m缓冲区 | 0.571 428 | 0.010 408 | 0.428 571 | 0.989 591 |
岩浆岩150 m缓冲区 | 0.500 000 | 0.093 743 | 0.500 000 | 0.906 256 |
与Au有关的蚀变带 | 0.285 714 | 0.009 333 | 0.714 285 | 0.990 666 |
Au异常 | 0.928 571 | 0.143 670 | 0.071 428 | 0.856 329 |
As异常 | 0.785 714 | 0.214 160 | 0.214 285 | 0.785 839 |
Sb异常 | 0.642 857 | 0.197 109 | 0.357 142 | 0.802 890 |
Ag异常 | 0.785 714 | 0.244 137 | 0.214 285 | 0.755 862 |
注:P1表示金矿床(点)存在时预测变量出现的概率;P2表示金矿床(点)不存在时预测变量出现的概率;P3表示金矿床(点)存在时预测变量没有出现的概率;P4表示金矿床(点)不存在时预测变量没有出现的概率。
4.3 各预测变量的权重值
预测变量的权重计算是为后验概率做数据上的准备,计算公式如下:
其中:
式中:W+表示各预测变量出现时的权重值;W-表示各预测变量不存在时的权重值;C表示正负权重值的差值,C值越大该预测变量越能很好地指示找矿,C值越小表示该预测变量对找矿不能有很好的指示。C>0表明该预测变量对成矿有利,C<0表明该预测变量对成矿不利,C=0表明该预测变量对找矿缺乏指示意义。证据权重值统计结果见表2。
表2 研究区各预测变量证据权重值统计结果
Table 2
预测变量 | W+ | W- | C | W+方差 | W-方差 | C值排序 |
---|---|---|---|---|---|---|
地层 | 3.163 443 | 0.000 000 | 3.163 443 | 0.071 701 | 0.000 012 | 5 |
韧性剪切域 | 3.627 708 | -1.922 866 | 5.550 574 | 0.083 840 | 0.500 011 | 1 |
脆性断束50 m缓冲区 | 4.005 541 | -0.836 835 | 4.842 376 | 0.126 109 | 0.166 678 | 2 |
岩浆岩150 m缓冲区 | 1.674 045 | -0.594 714 | 2.268 759 | 0.142 980 | 0.142 869 | 7 |
与Au有关的蚀变带 | 3.421 337 | -0.327 094 | 3.748 431 | 0.251 237 | 0.100 011 | 4 |
Au异常 | 1.866 123 | -2.483 957 | 4.350 080 | 0.077 003 | 1.000 013 | 3 |
As异常 | 1.299 868 | -1.299 442 | 2.599 310 | 0.090 963 | 0.333 348 | 6 |
Sb异常 | 1.182 161 | -0.810 082 | 1.992 243 | 0.111 169 | 0.200 014 | 9 |
Ag异常 | 1.168 861 | -1.260 549 | 2.429 410 | 0.090 956 | 0.333 348 | 8 |
4.4 预测变量独立性检验
证据权重模型进行成矿预测要求参与预测的各变量都是条件独立或弱独立的,如果两个或几个预测变量具有相关性,会导致后验概率失真,因此需对参与预测的变量进行独立性检验。
本次对各预测变量做显著水平为0.05下的条件独立性检验,结果显示:参与预测的9个预测变量均相互独立,全部可以参与后验概率计算。预测变量的独立性检验结果见表3。
表3 研究区各预测变量独立性检验结果
Table 3
预测变量 | 地层 | 韧性剪 切域 | 脆性断束 50 m缓冲区 | 岩浆岩150 m 缓冲区 | 与Au有关 的蚀变带 | Au异常 | As异常 | Sb异常 | Ag异常 |
---|---|---|---|---|---|---|---|---|---|
地层 | — | ||||||||
韧性剪切域 | -40.25独立 | — | |||||||
脆性断束50 m缓冲区 | -34.55独立 | -30.20独立 | — | ||||||
岩浆岩150 m缓冲区 | -15.58独立 | -20.94独立 | -12.81独立 | — | |||||
与Au有关的蚀变带 | -37.23独立 | -37.30独立 | -32.38独立 | -12.52独立 | — | ||||
Au异常 | -37.11独立 | -43.73独立 | -56.12独立 | -40.98独立 | -18.59独立 | — | |||
As异常 | -36.82独立 | -35.25独立 | -29.49独立 | -37.91独立 | -27.72独立 | -9.16独立 | — | ||
Sb异常 | -43.77独立 | -53.54独立 | -49.74独立 | -34.05独立 | -29.33独立 | -18.82独立 | -19.80独立 | — | |
Ag异常 | -38.82独立 | -35.34独立 | -21.08独立 | -16.02独立 | -14.72独立 | -19.54独立 | -10.75独立 | -15.30独立 | — |
4.5 后验概率的计算及预测成果的表达
在完成前述工作的基础上,利用MRAS软件计算出每个网格单元成矿的后验概率。预测单元的后验几率Q计算公式如下:
其中:
后验概率:P=Q/(1+Q)。
本次预测计算出的后验概率最大值为0.999 997,最小值为0。根据后验概率累积频率(图3),区内的后验概率可分为5个区间,但是后验概率处于0.457~0.573之间的网格单元数量较少(图4),且较为分散;经过对比,把该部分网格单元并入后验概率为0.573~0.871的预测网格单元,更有利于圈出合理的成矿远景区,因此,选择处于0.185~0.457、0.457~0.871和>0.871的3个后验概率值较大区间的预测网格单元作为预测区块。其中网格后验概率P>0.871的为A级成矿远景区块,0.457<P<0.871的为B级成矿远景区块,0.185<P<0.457的为C级在矿远景区块。图5为用色块图表示的预测成果。
图3
图3
后验概率累积频率
Fig.3
Distribution of posterior probability’s cumulative frequency
图4
图4
预测单元后验概率分布
Fig.4
Histogram of posterior probability distribution of prediction unit
图5
从老湾金矿带成矿预测后验概率色块(图5)可以看出,已知金矿床或矿点多数落入A级或B级成矿远景区块,C级成矿远景区块主要是在A级或B级成矿远景区块的周边,这与成矿区带从中心向边缘渐变的特征相吻合。
4.6 找矿远景区的圈定
根据预测网格分布区块特点及老湾地区金矿地质特征,本次预测圈出了8个找矿远景区(见表4)。根据成矿远景区块后验概率平均值把找矿远景区划分成Ⅰ级找矿远景区4处,Ⅱ级找矿远景区4处。
表4 老湾地区金矿找矿远景区
Table 4
编号 | 后验概率最大值 | 后验概率平均值 | 面积/km2 | 已发现矿床点数 |
---|---|---|---|---|
Ⅰ1 | 0.999997 | 0.860023 | 7.17 | 4 |
Ⅰ2 | 0.999997 | 0.852098 | 2.79 | 2 |
Ⅰ3 | 0.999997 | 0.876286 | 3.47 | 2 |
Ⅰ4 | 0.999868 | 0.874190 | 1.47 | 1 |
Ⅱ1 | 0.996198 | 0.773105 | 0.80 | 1 |
Ⅱ2 | 0.99997 | 0.624653 | 0.96 | 2 |
Ⅱ3 | 0.999965 | 0.781311 | 1.24 | 0 |
Ⅱ4 | 0.997268 | 0.598727 | 1.22 | 1 |
5 结论
1) 研究分析了老湾地区金矿的控矿因素和成矿条件,确定了该区金矿化受地层、构造、岩浆岩、围岩蚀变和地球化学异常等多因素控制。
2) 提取构置了9个预测变量,分别为地层、韧性剪切域、脆性断束50 m缓冲区、岩浆岩150 m缓冲区、与Au有关的蚀变带以及Au、As、Sb、Ag 1:5万水系沉积物单元素异常;基于证据权重模型,采用MRAS软件对老湾地区金矿进行了成矿预测,划分出A、B、C三级成矿远景区块,圈出Ⅰ级找矿远景区4处,Ⅱ级找矿远景区4处。
3) 结合老湾地区金矿地质特征,综合分析后认为,圈出的找矿远景区中Ⅰ级找矿远景区具有较好的金矿找矿潜力,特别是彭家老庄至北杨庄的Ⅰ1找矿远景区,后验概率平均值高,成矿地质条件好,面积大,已发现矿床点数多,是老湾地区今后金矿勘查的重点区域。
参考文献
基于MRAS证据权重法的辽西锦州——阜新金矿化带潜力评价
[J].
Potential gold evaluation with weighting of evidence based on MRAS in Jinzhou-Fuxin gold metallogenic belt, western Liaoning Province, China
[J].
矿产资源评价系统(MRAS)在全国矿产资源潜力评价中的应用
[J].
Application of MRAS in national mineral resource potential assessment
[J].
基于证据权重法的海南岛金多金属矿成矿预测研究
[J].
Study on metallogenic prognosis for gold polymetallic deposits in Hainan island based on evidence weighting method
[J].
基于MRAS证据权重法的湖南怀化地区金矿成矿预测
[J].
Metallogenic prediction of gold deposits with weighting of evidence based on MRAS in Huaihua area, Hunan Province
[J].
基于证据权模型的赤峰——呼伦贝尔成矿带找矿预测
[J].
Prospecting prediction of Chifeng-Hulunbeier metallogenic belt based on weight of evidence model
[J].
基于MRAS的鞍山——本溪地区铁矿资源潜力分析
[J].
Analysis on the iron resource potential in Anshan-Benxi area based on MRAS
[J].
基于证据权重法的雪峰山中段金矿矿集区找矿远景预测
[J].
The gold metallogenic prognosis based on the weighting method of evidence in the middle Xuefeng mountain area
[J].
巴音戈壁盆地塔木素地区砂岩型铀矿预测及找矿方向
[J].
Prognosis and prospecting direction of sandstone-type uranium deposit in Tamusu area, Bayingebi basin
[J].
基于GIS和证据权重法的那仁宝力格地区铜多金属矿成矿预测
[J].
Metallogenic prognosis of copper polymetallic mineral resources in Narenbaolige area on basis of weights of evidence method and GIS
[J].
河南桐柏老湾金矿床氢氧氦同位素地球化学及成矿流体来源
[J].
Isotopic geochemistry of hydrogen, oxygen and helium, and ore-forming fluid sources of Laowan gold deposit in Tongbai, Henan Province
[J].
河南桐柏老湾金矿床地质特征及成因探讨
[J].
Geological features and genesis of the Laowan gold deposit in the Tongbai region,Henan
[J].
北淮阳构造带老湾金矿区构造与成矿作用的关系
[J].
The relationship between structures and mineralization of Laowan gold minning area in Northern Huaiyang tectonic belt
[J].
河南省老湾金矿床地球化学特征及矿床成因
[J].
Geochemical characteristics and genesis of the Laowan gold deposit in Henan Province
[J].
桐柏老湾金矿带成矿构造体系演化研究及其地质意义
[J].
The study on the evolution of metallogenic tectonic system and its significance of Laowan gold metallogenic belt in Tongbai county
[J].
老湾花岗岩体与金成矿关系新认识及其找矿意义
[J].
New understanding of the Laowan granite body-Au mineralization relation and the significance to ore-search breakthrough
[J].
桐柏山老湾金矿带右行走滑断裂控矿体系的构建及其意义
[J].
Establishment and significance of dextral strike-slip fault ore-controlling system of the Laowan gold belt, Tongbai mountains
[J].
河南桐柏老湾花岗岩体锆石SHRIMP U-Pb年龄及其地质意义
[J].
Zircon SHRIMP U-Pb dating of Laowan granite in Tongbai area, Henan Province, and its geological implications
[J].
桐柏地区金银多金属矿集区找矿前景预测研究报告
[R].
The prediction and research report of prospecting potential in gold silver polymetallic ore concentration area, Tongbai
[R].
桐柏地区变质杂岩和侵入岩类Pb同位素组成特征及其地质意义
[J].
Pb isotopic compositions of metamorphic rocks and intrusive rocks in Tongbai region and their geological implication
[J].
桐柏山老湾金矿带与燕山期岩浆作用有关的岩浆热液金多金属矿床成矿作用——来自地球化学、年代学证据及控矿构造地质约束
[J].
Magmatic hydrothermal gold and polymetallic metallogenesis related to Yanshanian magmatism of Laowan gold belt, Tongbai Mountain:Evidence from geochemistry, geochronology and ore-controlling structural geological constraints
[J].
老湾花岗岩及成矿物质与岩浆岩的同源性研究
[J].
The research of homology of Laowan granite, metallogenic material and magmatic rock
[J].
河南老湾金矿床流体包裹体及稀土元素地球化学研究
[J].
Geochemistry of fluid inclusions and rare earth elements from Lanwan gold deposit in Henan Province
[J].
北秦岭老湾金矿带变质岩原岩恢复及其形成过程
[J].
Restoration and formation process of metamorphic rocks in the Laowan gold ore belt, North Qinling
[J].
豫西南老湾金矿原生晕地球化学特征及深部成矿预测
[J].
Geochemical characteristics of primary halo and deep metallogenic prediction of Laowan gold deposit in Southwestern Henan Province
[J].
基于MRAS的证据权重法在香格里拉地区的综合信息成矿预测
[J].
Multi-source information metallogenic prognosis with weighting of evidence based on MRAS in Shangri-La
[J].
河南省桐柏县老湾金矿深部及外围普查报告
[R].
Deep and external prospecting report on the Laowan gold deposit in Tongbai County, Henan Province
[R].
构造叠加晕法是危机金矿山寻找接替资源的有效新方法
[J].
Method of structural overprinting geochemical halo, a new and effective approach to prospecting for succeeding resources in gold deposits with resources crisis
[J].
基于证据权重法的成矿预测——以广东省钨锡矿的成矿预测为例
[J].
Minerogenetic prediction based on the weight-of-evidence approach: A case study of the prediction of tungsten and tin deposits in Guangdong, China
[J].
/
〈 |
|
〉 |
