基于浓度-面积分形模型和模糊证据权的铜矿资源潜力评价：以江西九瑞地区为例.pdf

引文格式：林俞亨，王立立，欧阳永棚，等.基于浓度面积分形模型和模糊证据权的铜矿资源潜力评价：以江西九瑞地区为例J.西北地质，2024，57（1）：165178.DOI：10.12401/j.nwg.2023199Citation：LIN Yuheng，WANG Lili，OUYANG Yongpeng，et al.Evaluation of Copper Mineral Resource Potential Us-ing ConcentrationArea Fractal Model and Fuzzy Evidence Weighting:A Case Study of the Jiurui Region inJiangxiJ.Northwestern Geology，2024，57（1）：165178.DOI：10.12401/j.nwg.2023199基于浓度面积分形模型和模糊证据权的铜矿资源潜力评价：以江西九瑞地区为例林俞亨1，王立立1，欧阳永棚1,2,3，*，李增华1，曾闰灵2,3，陈祺2,3，邓友国2,3（1.东华理工大学地球科学学院，江西 南昌330013；2.江西省地质局第十地质大队，江西 鹰潭335001；3.鹰潭市紧缺及优势矿产勘查与研究重点实验室，江西 鹰潭335001）摘要：中国江西省的九瑞地区是长江中下游成矿带中最重要的铜矿产地之一，其中花岗闪长斑岩与铜成矿关系密切。基于水系沉积物与矿化相关的信息，采用因子分析（FA）、浓度面积分形法（CA）和模糊证据权方法（FWofE）相结合建立成矿潜力预测模型。使用因子分析处理包含32 个元素的 255 份水系沉积物样本数据，找到能够指示铜矿化的组合元素（即主因子）。采用多重分形反距离加权插值法（MIDW）创建主因子得分栅格图并用 CA 分形模型提取与铜矿化相关的地化异常。将得到和铜矿化相关的地球化学异常图与地质、遥感解译数据相结合，应用模糊证据权方法建立预测模型。结果表明：已知铜矿床位于圈定预测概率高值区，且受花岗闪长斑岩和断裂的分布共同控制；除已知铜矿床区域外，圈定的 3 个一级远景区域内也具有较高的概率，值得进一步铜勘查找矿工作的进行。关键词：矿产资源评价；铜矿；模糊证据权；CA 分形模型；九瑞地区中图分类号：P618.41文献标志码：A文章编号：1009-6248（2024）01-0165-14Evaluation of Copper Mineral Resource Potential Using ConcentrationArea Fractal Model andFuzzy Evidence Weighting:A Case Study of the Jiurui Region in JiangxiLIN Yuheng1，WANG Lili1，OUYANG Yongpeng1,2,3，*，LI Zenghua1，ZENG Runling2,3，CHEN Qi2,3，DENG Youguo2,3（1.School of Earth Sciences,East China University of Technology,Nanchang 330013,Jiangxi,China；2.10th Geological Party ofJiangxi Bureau of Geology,Yingtan 335001,Jiangxi,China；3.Yingtan Key Laboratory of Exploration and Research ofScarce and Advantage Minerals,Yingtan 335001,Jiangxi,China）Abstract：The Jiurui region in Jiangxi Province,China,is one of the most significant copper mining areas inthe middle and lower reaches of the Yangtze River mineralization belt,with a close relationship between gran-收稿日期：2023-10-08；修回日期：2023-11-18；责任编辑：曹佰迪基金项目：江 西 省 科 技 厅 重 点 研 发 计 划 项 目（20212BBG73045），江 西 省 地 质 局 青 年 科 学 技 术 带 头 人 培 养 计 划 项 目（2022JXDZKJRC02），鹰潭市科技计划项目（20233-185656）联合资助。作者简介：林俞亨（1998），男，硕士研究生，主要从事综合信息矿产预测研究。Email：。*通讯作者：欧 阳 永 棚（1988），男，博 士 研 究 生，高 级 工 程 师，主 要 从 事 勘 查 地 质 学 和 区 域 成 矿 学 研 究。Email：。第 57 卷 第 1 期西 北 地 质Vol.57No.12024 年（总 233 期）NORTHWESTERN GEOLOGY2024（Sum233）odiorite porphyry and copper mineralization.In this study,a predictive model for mineralization potential wasestablished by combining factor analysis(FA),concentration-area(C-A)fractal method,and fuzzy weight of ev-idence(FWofE)based on information related to stream sediment and mineralization.factor analysis was ap-plied to a dataset of 255 stream sediment samples containing 32 elements to identify combinations of elements(principal factors)indicative of copper mineralization.the principal factor scores were interpolated using themultiple inverse distance weighted(MIDW)method to create a raster map,and the C-A fractal model was em-ployed to extract geochemical anomalies associated with copper mineralization.the geochemical anomaly maprelated to copper mineralization was integrated with geological and remote sensing interpretation data,and a pre-dictive model was established using the fuzzy weight of evidence method.The results indicated that:known cop-per deposits are located within high-probability zones defined by the model and are influenced by the distribu-tion of granodiorite porphyry and faults;in addition to the known copper deposit areas,three primary prospec-tive areas identified within the defined regions also exhibit a high probability,meriting further exploration ef-forts for copper prospecting.Keywords：mineral resources evaluation；copper mine；fuzzy weight of evidence；CA fractal model；Jiu-rui region铜具有优良的导热、导电性能和较强的延展性，因其优良的特性，被广泛应用于各个领域并在日常基础设施和技术中发挥着不可替代的作用（Ayres et al.，2003）。世界铜产量在过去 70 年稳步增长，其中，中国是对世界铜产量做出重大贡献的国家之一，具有丰富的铜矿资源（Mudd et al.，2018）。此外，长江中下游成矿带是中国最具悠久历史的多金属成矿带之一（王兆强等，2023），位于该成矿带西部的江西省九瑞地区是中国重要的铜产区，该区域有较好的铜找矿潜力（徐耀明，2014）。另一方面，如地球化学测绘、证据权重法等矿床资源定量评价的方法，能够结合地质、地球物理、地球化学、遥感等数据进行信息综合，进行预测远景区的圈定，对指导找矿勘查具有重要意义（欧阳渊等，2023）。因此，可以利用已有的地球化学数据和地质数据等其他信息来进行处理，对九瑞地区指导铜矿产勘查（薛琮一等，2020；李文明等，2021；Liet al.，2022）。地球化学微量元素作为成矿的直接物质来源，一直是地质工作者研究的重点，如何快速、高效的识别地球化学异常则是处理地球化学数据的一大难题。其中，分形/多重分形模型是识别地球化学异常的强有力的工具（Zuo et al.，2011，2012，2013）。自从 Man-delbrot（1982）提出分形理论以来，已经有许多成果利用分形/多重分析理论从背景中提取出了异常。例如Cheng 等（1994）提出了浓度面积（CA）分形模型，并成功应用在米切尔硫酸盐区金铜异常的提取。各种分形方法，包括能谱面积分形模型（SA）、浓度数量分形模型（CN）、分形奇异值分解法（MSVD）等都能较好的用于不同矿化系统，并成功区分出地球化学异常和背景（Cheng，2004；Ali et al.，2007；Hassan-pour et al.，2013）。源自于没有空间意义的医学诊断上发展而来的证据权模型（Weight of Evidence，简称 WofE）被国内外学者广泛应用于多元信息综合和空间决策支持系统，它可以整合多源地质信息作为诸多证据因子，从而形成潜在矿产分布的后验概率图（Agterberg，1989；Agter-berg et al.，1994；Keller，1995；张津瑞等，2023）。Agter-berg 等（1990）首次将证据权方法用于成矿预测中。在此基础上，Cheng 等（1999）发展出了模糊证据权重法。模糊证据权重（FWofE）是一种用于支持决策预测目的的数据集成方法,是一种空间决策支持方法。相比于普通证据权方法，模糊证据权重法通过引入模糊集合理论和隶属度函数，可以有效地减少由于遗漏或错误标记的数据引起的不确定性和提高地质数据的可信度，在确定模糊隶属度函数以代替后验概率计算的缺失数据方面更加灵活，为地质分析和决策提供更可靠的基础，其已被广泛应用于矿产潜力预测中（Cheng et al.，2002）。成秋明等（2007）采用模糊证据权等技术对老王寨及其邻区的金矿资源进行潜力评价，成功圈定有利成矿地段。Xiao 等（2020）综合普通证据权、模糊证据权和逻辑回归 3 种方法，利用构造、岩性、地球物理和地球化学等数据，绘制新疆东天山造山带的斑岩型铜钼多金属矿床前景并进行定量评价。Zhang 等（2016）运用模糊证据权法得到预测区域166西北地质NORTHWESTERN GEOLOGY2024 年仅占总研究区的 5%且成功识别了 65%的研究区内的铁矿床。Porwal 等（2006）开发的一种基于模糊证据权重方法和模糊逻辑方法相结合的子模型，该子模型在确定高矿产潜力区域方面是有效的，可以作为矿产勘探的工具。在基于多元信息和信息综合技术的矿产资源定量评价与成矿预测领域，模糊证据权重法是行之有效且应用广泛的方法（Cheng，2008a；黄秀等，2010；Zhang et al.，2014；陈风河等，2015；武进，2017；邓军等，2021；Huo et al.，2022；黄鑫怀等，2023；霍雨佳，2023；王佳营等，2023）。江西九瑞地区位于长江中下游铜多金属成矿带中，该地区长期的构造作用、岩浆活动和成矿作用形成了丰富的铜等多金属矿床，具有勘探铜矿的潜力。但区域成矿地质条件复杂，致使传统的地质数据提取找矿信息难度较大。笔者基于九瑞地区 120 万水系沉积物地球化学采样数据、地质构造和控矿岩性等数据，在 ArcGIS、GeoDAS 等软件平台上，综合利用因子分析、多重分形反距离权重插值、CA 分形模型和模糊证据权等方法，构建九瑞地区模糊证据权模型，研究九瑞地区铜矿成矿潜力靶区，圈定可靠的成矿预测远景区，给予未来矿产勘探提供一定的参考意义。1地质背景长江中下游成矿带位于扬子克拉通北缘，秦岭大别造山带和华北克拉通以南。成矿带周边发育有3 条大型断裂，分别为东西向横贯成矿带南侧的阳新常州断裂，西北缘呈 NW 向的襄樊广济断裂以及东北缘的郯庐断裂。研究区九瑞地区位于江西省北部，属长江中下游成矿带之九瑞铜金矿集区（图 1）。九瑞矿集区庐枞宁芜与晚中生代铜-金-钼矿床相关的侵入岩与白垩纪磁铁矿矿床相关的火山岩-次火山岩斑岩-矽卡岩-块状硫化物铜-金-钼-铁矿床磁铁矿-磷灰石矿床扬子板块安庆-贵池大别造山带华北克拉通宁镇鄂东南襄樊-广济断裂九瑞阳新-常州断裂铜陵郯庐断裂长江0N4080 km图1长江中下游成矿带简易地质图（据 Pan et al.，1999 修）Fig.1Simplified geological map of the Mid-Lower Yangtze metallogenic belt九瑞地区地层发育较齐全，以古远古界为结晶基底，中远古界为褶皱基底，晚震旦世以来的地层为盖层，由老到新为奥陶系至第四系，其中主要出露奥陶系至三叠系，研究区内矿床的形成、分布与地层关系密切，区内铜矿床的主要赋矿地层为泥盆系五通组、石炭系黄龙组和二叠系三叠系大冶组（孔凡斌，2014；徐耀明，2014）。区域内曾经历过强烈的岩浆活动，主要活动时期为燕山期，岩石类型以中酸性侵入岩为主（蒋少涌等，2013），产状多为小型岩株、岩脉和岩枝。研究区内岩体种类繁多，以花岗闪长斑岩和石英闪长玢岩为主，另分布有石英斑岩、花岗细晶岩、闪长岩、辉绿岩等。其中，燕山期花岗闪长斑岩是研究区内最主要的铜成矿母岩。研究区内褶皱、断裂构造发育，构造系统分级明显，严格控制了岩浆岩空间展布与矿床的分布。区域第 1 期林俞亨等：基于浓度面积分形模型和模糊证据权的铜矿资源潜力评价：以江西九瑞地区为例167 构造主要为 NW 向与 NEE 向，其次为 NNW 与 NNE向，断裂及其构成的网格状构造节点，控制和影响着岩浆活动与成矿作用（图 2）。九瑞地区成矿岩体及其围岩普遍经受热液蚀变作用，矿床岩体多以岩体中心向外形成环状蚀变分带，与铜矿化相关的蚀变主要有矽卡岩化、硅化、绿泥石化。其中，据徐耀明(2014)制作的遥感蚀变光谱图可以得知，在空间上绿泥石化与已知矿床表现的较为密切。矿体以层状、似层状、透镜状为主，铜矿主要类型有矽卡岩型、层控硫化物型、斑岩型等，区内截至2010 年累计探明的铜储量为 0.3 万 t（所颖萍，2013）。在模糊证据权模型的构建的过程中，需要与成矿相关的地质要素转化为证据图层参与模型的建立。根据九瑞地区的地质背景，选择铜矿化的主要赋存地层、铜成矿相关的花岗闪长斑岩、影响铜矿空间分布区域断裂以及与铜矿空间分布的有关的绿泥石化蚀变作为证据图层参与模型的构建。这些地质要素与铜矿的分布和潜力密切相关，有助于模糊证据权模型更精确地评估矿产资源的分布和潜力。此外，通过对水系沉积物中提取的地球化学数据进行分析，可以识别出与铜矿化相关的地球化学异常，通常这些异常指示潜在的铜矿化区域。因此，地球化学异常也可以作为一类证据图层用于模型的建立，以进一步提高对矿产资源潜力的评估准确性。2研究方法 2.1因子分析因子分析是常用的多元统计分析方法之一，主要用于降维，并可以通过确定用于解释所观察到的相关矩阵所需的公共因子的适当数量和性质来计算出被测变量的更简约结构（Fabrigar et al.，1999）。同时也可以用于解决复杂的地质成因及矿化叠加问题，对叠加地球化学场也有较好的分解作用（董庆吉等，2008）。因子分析步骤包括确定因子载荷、计算变量共同度、因子旋转、方差贡献率和计算因子得分。其中，因子得分是观测变量的线性组合，它可以代替原来的变量用于后续的分析。2.2多重分形反距离权重法在地球化学异常提取分析的过程中常有将点样本的属性转换为连续面的数值插值的处理过程。传统的插值方法（如普通克里金和反距离权重法）是空间插值分析中常用的方法。Cheng（2000）所提出的多重分形反距离权重插值（MIDW）在对地球化学数据处理在局部奇异性中表现出不错的优势。相比于传统的插值方法，多重分形方法能增强地球化学特征的预测能力，能更好的确定矿产勘探目标，其插值结果更 第四系燕山期花岗闪长斑岩三叠系燕山期石英闪长玢岩二叠系铜矿石炭系断裂城门山丁家山N02.55.0 km瑞昌市武山东雷湾宝山邓家山泥盆系志留系奥陶系图2九瑞铜矿区地质图（据 Yang et al.，2011 修）Fig.2Geological map of the Jiurui region ore copper district168西北地质NORTHWESTERN GEOLOGY2024 年能对弱异常进行精细识别（Parsa et al.，2017）。有关局部奇点指数和 MIDW 的更多详细信息见 Cheng（2007，2008b）和 Zheng 等（2014）。MIDW 方法广泛应用于地球化学异常找矿中，可以增强地球化学异常区域，有助于地球化学勘探（Cheng，2004，2008b，2015）。MIDW 的方程式可以表示为公式（1）：(0)=(0)2(0)（1）式中：是多重分形权重；是单位单元大小的尺度（1）；是奇异性指数。2.3CA 分形模型CA 模型是地球化学异常识别的强有力的工具，是人们常用的“地球化学异常建模的基本技术”（Carranza，2009）。浓度面积分形法由 Cheng 等（1994）首先提出，是用来分选地球化学异常的。浓度-面积分形法（CA）是基于浓度值（量值的频率）的从背景中分离异常值的方法，这些量值也可以是地球化学指标的空间和几何特征。双对数坐标下的 CA 图可以用来建立浓度值大于 s 的面积 A（s）和浓度值 s 之间的幂律定律关系。从而可绘制出高、中、低异常区域地图。CA 模型表达方式见公式（2）：A(v)1；A(v)2（2）A()12A()式中：表示浓度大于或等于轮廓值 的区域；是阈值；和是大于零的分形维数。这 2 个分形参数可以从对 的对数-对数图中最佳拟合直线的斜率来估计。2.4模糊证据权重FWofE 将证据层定义为具有多值隶属函数（0（A）1）而不是普通证据权重法中通常涉及的二元或三元集合。二元和三元模式成为模糊集的特例，例如（A）=1 或 0 对应于二元模式，（A）=0、0.5 和 1对应于三元模式。具有隶属函数 1 和 0 的模糊证据的 2 个子集在 FWofE 的实现中特别重要：A1=A:（A）=1或 A2=A:（A）=0。第一组 A1是证据的一部分，通常是勘探良好的区域，大部分矿床已被发现。第二组是发现的矿床很少，发现更多矿床的机会较小。假设与 A1和 A2相关的点事件（E）发生的条件概率为 PE|A1、PE|1、P|A1、P|1、PE|A2、PE|2、P|A2 和 P|2。E 与任何 A 相关的模糊条件概率计算过程见公式（3）与公式（4）：PE|(A)=(B)PE|A1+1-(A)PEA2（3）PE|(A)=(B)P|A1+1-(A)P|A2（4）其中带“bar”的字母表示没有发生事件或没有证据。模糊证据的权重计算过程见公式（5）：W(A)=logPE|(B)/P|(A)（5）式中：W+和 W-是普通证据权重法中涉及的 2 个普通权重；W0是先验权重。公式（5）具有以下特性：W（A）=W+，假设（A）=1 和 W（A）=W-，假设（A）=0 和 W（A）=W0，假设（A）=1/2 和 A1+A2=（整个研究区域）。在条件独立于不同证据 A 和 B 的假设下关于事件E 的发生，给定 A 和 B 的 E 的 logit 计算过程见公式（6）：logitE|(A),(B)=W0+W(A)+W(B)（6）以 应 用 反 演 变 换将 logit 转 换 回 后 验 概 率PE|（A）,（B），也可以计算后验概率的标准偏差，以测试结果的显著性。模糊证据权重方法的实施步骤概括为 6 个步骤：确定研究目标，如预测给定类型但未发现矿产的空间位置。确定与目标相关的空间层（属性）,以便他们可以用于评估目标。提取与目标相关的证据层，用隶属度函数表示证据层的可信度。计算模糊证据层的权重（衡量各个证据层对评估目标的重要性）。综合多个模糊证据层，计算后验概率形成矿产资源潜力分布图。模型检验和预测区解释。2.5接受者操作特征曲线（ROC）接收者操作特征曲线（Receiver Operating Charac-teristic Curve，简称 ROC）自被创建以来已被用于各个学科，是一种可以用来评价分类器性能的可视化方法（Egan，1975；Zou et al.，2007；Swets，2014）。此 外，ROC 也 可 以 用 于 测 试 空 间 预 测 模 型 的 有 效 性（Nyknen，2008；Chauhan et al.，2010；Rossi et al.，2010）。其主要原理是将异常空间分布模式与成矿有利要素或与已知矿床（点）进行比较，并通过相关统计指标来度量两者间的空间对应关系（周光锋，2021）。ROC 曲线通过反映真阳性（True Positive，TP，模型将阳性分为阳性）和假阳性（False Positive，FP，模型将阴性分类为阳性）以及真阴性（True Negative，TP，模型将阴性分类为阴性）和假阴性（False Negative，FN，模型将阳性分类为阴性）的权衡关系作为一种评估分类器性能的方法，ROC 曲线是纵轴上的灵敏度（真阳性率：TP/（TP+FN）与横轴上的 1-特异性（假阳性率：FP/（FP+TN）的比较图，其曲线下面积（Area Under the Curve，简称 AUC）可用作诊断测试准确性的度量，也可用于测量空间预测模型的性能。AUC 值可能在 0 到 1 之第 1 期林俞亨等：基于浓度面积分形模型和模糊证据权的铜矿资源潜力评价：以江西九瑞地区为例169 间变化，如果 AUC 值等于 0.5，则分类性能相当于完全随机猜测；而如果 AUC 值等于 1，则分类性能完美，即分类器可以正确分类所有样本（图 3）。分类器的AUC 值通常介于 0.5 和 1 之间，其值越高代表的预测效果越好（Nyknen et al.，2015；Chen et al.，2016）。最优分类ROC 曲线00.51.00.501.0优差1-特异性敏感度随机分类图3ROC 曲线示例图Fig.3Example plot of ROC curve 3数据与模糊证据权模型构建 3.1数据来源本研究用于建立预测模型的数据集有地质数据、化探异常数据和遥感解译数据，数据包括铜矿床空间分布位置数据、九瑞地区断裂数据、主要赋矿地层数据、花岗闪长斑岩分布数据、绿泥石化蚀变数据、水系沉积物化探异常空间分布数据，这些数据经过处理后作为证据图层输入模糊证据权模型中。其中，断裂、地层、花岗闪长斑岩的数据提取自 1:25 万的九瑞地区构造图。绿泥石化蚀变数据来自徐耀明（2014）对九瑞地区遥感蚀变光谱图的解译结果。化探异常数据来自 1:20 万水系沉积物地球化学数据的分析结果。此外，1:20 万水系沉积物地球化学采样数据由江西省地质局第十地质大队按照采样密度每 4 km2一个样品点位采集所得，且 225 件样品点基本均匀覆盖整个研究区域（图 4）。各微量元素检出限的测定包括：电感耦合等离子体光谱法测定 Be、Co、Cr、Cu、Ba、K2O、La、Li、Mn、Nb、Ni、P、Sr、Th、V、Y、Ti、Zn、Zr 的浓度；掩蔽光谱法测定 As、Bi、Cd、Sb 和 Pb 的浓度；用比色法测定 U 的浓度；催化极谱法测定 Mo 和W 的浓度；Ag 和 Au 的浓度通过原子吸收光谱法测定 B、F 和 Hg 的浓度分别通过发射光谱法、离子选择电极法和冷原子法测定。所有 32 个元素的检测限制如表 1 所示。02.55.0 km采样点N图4水系沉积物样品采样点分布图Fig.4Steam sediment sample location map 3.2数据预处理对水系沉积物样品的 32 个元素进行分析处理，用以提取九瑞地区的铜矿化异常分布图作为构建模糊证据权模型的证据图层。首先对标准化后的数据进行因子分析。因子分析方差最大正交旋转因子载荷表 2 显示了各种因素的组成，其中 F1 因子中的 Cr、170西北地质NORTHWESTERN GEOLOGY2024 年Co、Ni、Ti、V 为亲铁元素，主要存在于基性岩、超基性岩中，可能表明了铁的矿化。因子 F2 含有 Cu、Au、Bi、Mo、Pb、Sb、W 是研究区主要的成矿元素，相应的 Cu、Au、Mo 元素属于矽卡岩-斑岩型铜钼金硫成矿亚系列，而 Pb、Sb、Bi 为亲铜元素，W 则常产出于矽卡岩型铜矿中，因此 F2 因子可以作为研究区主要铜成矿元素空间聚集的综合指标。因子 F3-F7 与区域成矿无明显关系，对于成矿的指示性不强。通过 GeoDAS 软件，使用多重分形 IDW 方法对F2 因子得分（Fac2）的空间分布图进行插值（图 5），发现其现有的矿点与异常区域吻合程度较差，说明只进行插值不足以从背景中提取异常，需要进一步分析来提高异常区域提取的准确性。将浓度面积（CA）模型应用于因子 F2 多重分形 IDW 的结果进行地球化学异常与背景分离。为了研究元素含量与对应面积之间的分形特征值，通过GeoDAS 软件处理，做出对应的浓度-面积双对数散点图。在双对数图上背景值和异常值可以对应于不同的斜率值，不同线段所对应的分界值可作为区分背景和异常的临界值（图 6）。经过分析，利用最小二乘拟合将图上的离散点拟合成 4 条直线，从左往右分别表示为背景、弱异常、异常和高异常，最终得到研究区Cu 矿化异常分级图（图 7）。此外，地质特征一定范围缓冲区与成矿是有关联的，并可以提供重要的勘探指南，所以可以对这些地质要素“延伸”一定的距离作为研究的证据图层（Zhang et al.，2016）。对于地层、断裂、花岗闪长斑岩和绿泥石化蚀变等地质要素，可以构建相应的缓冲区作为证据图层。根据已知矿床与各地质要素的空间 表 1 32 种元素的检出限表Tab.1Detection limits of 32 elements序号元素检出限序号元素检出限1Ag0.0117Mo0.52As2.8218Nb53Au0.000 319Ni54B520P305Ba1021Pb5.46Be0.522Sb0.27Bi0.1623Sn0.148Cd0.124Sr59Co125Th5.110Cr7.226Ti3011Cu127U112F1328V213Hg0.0129W0.514La1030Y1015Li531Zn1016Mn3032Zr10注：元素含量为10-6。表 2 R 型因子分析的正交旋转因子载荷矩阵表Tab.2Orthometric rotating factor loading matrix forR-factor analysis变量因子载荷F1F2F3F4F5F6F7Ag0.1010.0310.9540.0650.0300.0390.062As0.0590.3290.0990.8540.0390.0270.203Au0.0120.7170.0410.0070.0260.0050.069B0.3030.0720.0210.0720.7160.1630.104Ba0.1010.0310.9540.0650.0300.0390.062Be0.9100.0230.1750.0680.0150.1700.016Bi0.0250.9410.0060.0620.0170.0320.007Cd0.0550.2490.0980.9170.0990.0140.102Co0.7670.0060.0930.4350.2240.0210.112Cr0.8310.0010.0070.0610.0130.0820.036Cu0.0890.7860.0220.3240.0220.0140.288F0.8350.0430.0570.0510.0870.1020.123Hg0.0630.2210.0230.1960.1000.1730.443La0.2350.0930.0680.1250.1530.7600.047Li0.8900.0190.1240.0130.0590.1770.038Mn0.4890.0050.0650.3280.1980.0160.212Mo0.0090.4890.0150.2290.0150.0440.679Nb0.2540.1040.1180.0450.6110.1230.029Ni0.8750.0690.0490.0700.0260.1810.052P0.4640.0490.4170.0250.2260.2180.100Pb0.0750.9180.0140.0500.0020.0260.157Sb0.0480.9260.0620.1690.0220.0590.053Sn0.2740.2870.0640.0530.0940.5720.068Sr0.2090.0580.7360.0770.4110.0930.068Th0.3320.0150.2410.0100.6030.1160.094Ti0.5690.0110.5860.0060.2520.0020.110U0.1700.0220.1930.0580.0730.0890.289V0.9310.0420.1460.0490.0360.1340.032W0.0770.6570.1710.1440.0640.0220.444Y0.1080.2160.1810.0310.1520.7390.141Zn0.0260.0490.0170.9590.0550.0230.041Zr0.5700.0700.1660.0980.5460.0970.007注：该因子分析采用的提取方法为主成分分析法，旋转方法为Kaiser标准化最大方差法，旋转在七次迭代后已经收敛。第 1 期林俞亨等：基于浓度面积分形模型和模糊证据权的铜矿资源潜力评价：以江西九瑞地区为例171 关系，当主要赋矿地层（泥盆系五通组、石炭系黄龙组和二叠系三叠系大冶组）的多环缓冲区为 1 km、花岗闪长斑岩的多环缓冲区为 0.5 km、断裂的多环缓冲区为 1.3 km、绿泥石化蚀变的多环缓冲区为 0.8 km 时，平均有 85%以上的矿床落入各个地质要素的多环缓冲区内。在这些缓冲区范围外，随着距离的增加，所包含的矿床的数量增加极少。因此，这些距离可以视为最有利距离，即在这些缓冲距离内与已知矿床的空间相关性最大。经过数据预处理后，共得到了 5 个建 0510 km因子 F2_MIDW9.270铜矿床N图5因子 2 得分的多重分形反距离权重插值结果图Fig.5Multifractal IDW result for the factor 2 scores 53.84.04.24.44.64.85.043Log(浓度)C-ALog(面积)210图6浓度与面积的双对数图Fig.6Log-log plot of concentration versus area 00.01 背景铜矿床0.010.22 弱异常0.220.56 异常0.569.27 高异常0510 kmN图7CA 分形模型识别的异常图Fig.7Anomaly map identified by the CA fractal model172西北地质NORTHWESTERN GEOLOGY2024 年立预测模型所需的证据层：花岗闪长斑岩多环缓冲区证据层。断裂多环缓冲区。主要赋矿地层多环缓冲区证据层。地球化学异常证据层。绿泥石化蚀变多环缓冲区（图 8）。0 2.5 5.0 km0 2.5 5.0 km弱异常异常高异常背景0 2.5 5.0 km0 2.55.0 km0 2.5 5.0 km100 m200 m300 m400 m500 m100 m600 m1 000 m200 m700 m1 100 m300 m800 m1 200 m400 m900 m1 300 m500 m100 m600 m200 m700 m300 m800 m900 m400 m500 m1 000 m100 m600 m200 m700 m300 m800 m400 m500 mabcde铜矿床铜矿床铜矿床铜矿床铜矿床NNNNN图8用于九瑞地区成矿预测的证据图层Fig.8Evidence layers for mineral prediction in the Jiurui region 3.3模糊证据权模型在构建模糊证据权模型前，需要在 GeoDAS 软件对研究范围和相关参数进行定义，研究范围可以分为训练区域和预测区域进行定义，一般而言训练区域和预测区域为同一块区域。此外，为了确保每个栅格单元大小包含并仅包含一个矿床，设置研究区的面积单元大小为 400 000 地图单元，选择 11 个已知矿床作为训练图层，其先验概率为 0.003 6。在设置的训练参数基础上，对上述 5 种证据图层分别进行模糊权重计算来构建模糊证据图层。普通证据权方法一般是以对比度 C 与 C 的标准偏差 S（C）的比值 t（分布 t 值，t=C/S（C）的极大值作为各证据层二值化的阈值，该阈值将各证据图层分成截然不同的二值图，即阈值之前的分类部分对成矿有利，赋值为 1，阈值之后的分类部分对成矿不利，赋值为 0。而对于模糊证据权，证据层由隶属度函数（MSF）进行重分类，隶属度函数是描述某个元素或时间对于某个特定集合（或模糊集合）的隶属程度或归属度，在地质中可以描述地质数据中的不确定性和模糊性并对地质数据进行度量，因此可以使用 MSF 在闭区间 0,1 上任意取值对证据图层的确定性进行度量，根据已知矿点与不同分类的关系赋予不同 MSF，一般的根据矿床在不同分类上的分布情况确定各类的 MSF 取值，各证据图层的 MSF 取值如图 9 所示，其计算的模糊权重值如表 3 所示。4结果与讨论 4.1预测结果每一个模糊证据图层刻画了与已知矿点的空间关系，构建模糊证据图层的过程是训练和学习的过程，综合 5 个模糊证据图层创建后验概率图是达到预测和评价目标的必须过程。根据各证据图层的模糊证据权重，综合 5 个证据图层，进行后验概率的计算，得到后验概率图（图 10）。从后验概率图的结果来看，已知的 11 个矿床的位置都具有较高的后验概率，城门第 1 期林俞亨等：基于浓度面积分形模型和模糊证据权的铜矿资源潜力评价：以江西九瑞地区为例173 1.00.80.60.40.2隶属度函数01.00.80.60.40.2隶属度函数01.00.80.60.40.2隶属度函数01.00.80.60.40.2隶属度函数01234缓冲区分类值567812345671