老挝铜资源成矿规律与基于机器学习的远景预测

老挝处于特提斯成矿域南东段，具有丰富的矿产资源，但其地质工作基础薄弱，厘定矿产资源成矿规律并开展远景区预测是老挝在重点区实现找矿突破的有效途径。老挝 1∶100 万国家尺度地球化学填图由中老双方合作完成，为其矿产资源和环境评价提供了高质量的地球化学基础数据和图件。本文主要利用国家尺度地球化学填图数据，结合老挝已发现矿产成矿规律，利用机器学习技术，开展铜资源远景区预测。研究结果表明，（1）老挝铜矿床的形成明显受到构造-岩浆-沉积作用控制，铜矿床主要类型有斑岩型、矽卡岩型、热液型以及砂岩型。（2） 老挝全国水系沉积物中铜含量在 1.20~459 μg/g 之间，平均值为 21.96 μg/g，中位值为 16.50 μg/g，在 7 个三级大地构造单元中, 长山地块和哀牢山—马江等三个缝合带的平均值高于其他几个构造单元，地球化学图显示铜在老挝分布不均匀，存在多个大面积分布的高背景区和异常区。（3）构建了包括单元素异常、矿化元素组合异常、指示中酸性岩体元素组合、控矿构造分布、碳酸盐岩和碎屑岩分布等要素的老挝铜矿多源信息定量信息预测模型。（4）利用随机森林成矿预测方法，共圈定 9 个成矿远景区，具有寻找斑岩型、矽卡岩型等类型铜矿找矿前景。

关键词：远景区预测；机器学习；铜成矿规律；地球化学填图；老挝

老挝地处东南亚中南半岛中北部，在大地构造上位于印支板块的中东部，处于特提斯成矿域南东段，是我国三江成矿带向南延伸的部分。老挝具有丰富的矿产资源，目前已发现有铜、金、铅、锌、铁、钼、锑、锡、锰、 稀土、铝土矿、钾盐、石膏、煤等 20 余种矿产。已经发现的诸如富开（Phu Kham）、塞奔（Sepon）等超大型斑岩型铜金矿床， 甘蒙、万象盆地超大型钾盐矿床，表明老挝具有巨大的找矿潜力。

老挝是与中国陆地相连的友好国家，也是中国“一带一路”倡议的重要合作伙伴之一。由于老挝矿业投资环境较好，地理位置优越，资源潜力巨大，加之地质勘查程度低，矿业开发活动不多，具有较大勘探开发前景，并能为中国在铜、金、钾盐等战略性矿产资源供给起到互补性作用，已成为中国地勘单位和企业境外矿业开发与投资的首选之地。老挝的地质工作基础薄弱， 至今仍处于较低水平。在中国、 越南、 日本等国援助下， 累计完成的 1∶20万地质矿产填图面积不到全国总国土面积的 70%。地球化学填图工作更是开展的少，1∶20万区域化探仅在近十年完成了零星几个图幅。更大比例尺的地质矿产调查和化探工作仅由一些矿业公司围绕勘查需要完成了极小范围区域。

2014 年，中老双方在寻求加强地质调查合作背景下，本团队依托联合国教科文组织全球尺度地球化学国际研究中心与老挝能源矿产部地矿司共同开展了老挝全境 1∶100 万国家尺度地球化学填图，覆盖面积约 23 万平方千米，采集水系沉积物样品 1905 件，分析元素69 种。该合作项目作为“化学地球”国际大科学计划的一部分， 为老挝矿产资源和环境评价提供了高质量地球化学基础数据，制作完成老挝首套覆盖全境的地球化学图，获老方高度肯定， 取得了重大合作成果。1∶100 万超低密度地球化学填图方法可以达到迅速有效圈定找矿远景区的目的， 是勘查地球化学 “迅速掌握全局， 逐步缩小靶区”优势和特点的重要体现。本文主要利用中老合作项目取得的数据，结合老挝已发现矿产的成矿规律，利用机器学习技术，开展铜资源的远景区预测，为老挝下一步开展更大比例尺的地质和地球化学填图工作提供选区，为其矿产资源勘查与开发形成有力支撑。

1.1 大地构造与区域地质概况

老挝位于冈瓦那大陆与劳亚大陆的交接部位，重要断裂可划分为北东向、北西向两组。北东向断裂主要位于老挝西北部， 包含澜沧江—班南坎断裂、 南本河断裂、 琅勃拉邦断裂、普雷山断裂；北西向断裂位于老挝南东部， 包含长山—岘港断裂、蓝江断裂及马江断裂。可划分为 4 个一级构造单元，北西部地区为中缅马陆块、 昌都—思茅—南邦陆块， 属我国云南“三江”造山带的南延部分，中南部地区属印支陆块，北东部（华潘省东部） 为华南陆块。自西向东可划分为 7 个三级大地构造单元， 即景洪—素可泰火山弧、 思茅—彭世洛地块、奠边府—黎府缝合带、 万象—昆嵩地块、 色潘—三岐缝合带、 长山地块、 哀牢山—马江缝合带（图 1）。

老挝出露地层主要包括元古宙高级变质片麻岩和片岩; 古生代变质海相火山沉积岩系和少量片岩、灰岩等; 中生界砂岩、灰岩; 新生界砾岩、砂页岩、粉砂岩等。其中老挝的大部地区均有古生界和中生界出露，前中生界分布于北部及长山山脉，中、新生界分布于湄公河流域与西南部，元古宙地层零星分布于西北部和东部。老挝共有西北、北中、南部 3 个褶皱带。西北部北北东向印支期褶皱带主要由泥盆系—三叠系组成，北中部北西—北北西向海西期褶皱带主要由奥陶系—石炭系组成，南部北北西向印支期褶皱带主要由变形较弱的中生代陆相沉积岩组成。老挝断裂带较为发育，主要包括北东向、北北东向和北西向 3 组深大断裂及其次级断裂。

老挝岩浆活动从晋宁期、海西期、印支期、燕山期和喜山期均有分布，其中最为发育的为印支期。全国范围来看，超基性岩、基性岩、中性岩及酸性岩均有出现。晋宁期侵入岩主要零星分布于老挝东部，最主要岩性为片麻状黑云母和白云母花岗岩。海西期岩浆活动以喷发为主，中-基性火山岩是其主要岩性，分布于中北部和部分南部地区，其次还有部分花岗闪长岩和花岗岩在中部地区零星分布。印支期岩浆活动规模最大，主要岩性为英云闪长岩、花岗闪长岩等，分布于西北部、中北部和部分南部地区。喜山期岩浆活动总体上规模较小，主要是更新世玄武岩，岩性以碱性玄武岩类为主，分布于南部地区。

1.2 矿产资源概况

据老挝能矿部最新资料，目前已发现各类金属矿床和矿点（矿化点）近 500 处，其中铜矿 70 余处，金、银矿 160 余处，铅、锌、锑矿近 100 处，铁、锰、铬矿 60 多处，钨、锡、钼等矿 70 多处，此外还有铝土、钾盐、稀土等其他矿种 30 处以上。老挝的金属矿产主要集中分布在万象、沙耶武里、琅勃拉邦、川圹、华潘等中北部省区，以及甘蒙、沙湾拿吉和阿速波等中南部省区。

铜、金、铝土、锡、钾盐、铅锌、铁、煤等 8 种是老挝传统优势矿产资源，锰、铬、钨、钼、铋、镍、锑、汞、铂、铀、钍、铍、独居石等 13 种为其一般性矿产资源，迄今只部分开采了铜、金、锡、钾盐、煤等资源。近几年，随着全球对稀土等关键矿产资源的关注，老挝在稀土勘查工作方面取得了较大进展，寻找到了大量具有开发潜力的资源，将是该国下一步重点开采的资源之一。

老挝矿产资源的区域成矿时代主要集中在海西期、 印支期和喜山期， 其中海西—印支成矿期集中了老挝主要铁、有色金属和稀有金属矿产，并以热液型和矽卡岩型为主，喜山成矿期集中了老挝主要非金属岩盐、钾盐、石膏、宝石以及部分砂金、砂锡等矿产，矿床类型主要为沉积型、砂矿型和风化淋滤型。热液型、矽卡岩型、沉积型、砂矿型和风化淋滤型是老挝主要矿床类型。

矿产资源预测通过整合多源地球科学数据，进而量化某一位置矿床存在的可能性。机器学习是人工智能的重要分支，凭借其强大的数据处理和模式识别能力，成为矿产资源预测领域的一种主流技术。

随着大数据技术和人工智能在地球科学中的引入，各种机器学习方法用于矿产资源的定量评估和建模，例如逻辑回归（LR），支持向量机（SVM）、随机森林（RF）、卷积神经网络（CNN）、最大熵、玻尔兹曼机和极端学习机回归。这些方法通过处理多源地理空间数据集之间的复杂和非线性关系，为矿产前景建模提供了强大的技术。相较于传统统计学方法，其具备更强的泛化性与适用性。

随机森林作为机器学习中的一种的集成算法，通过对多源地球科学数据特征进行自动化提取，来完成不同的分类或预测任务，已被广泛应用到矿产资源定量预测中。相较于其他的机器学习方法，随机森林在矿产资源预测中表现突出，能够深入剖析数据分布规律，挖掘数据间的非线性特征，不仅对噪声和异常值有良好的鲁棒性，并且预测准确率也较高，能够在复杂多变的地质环境中提供可靠的预测结果。

老挝 1:100 万国家尺度地球化学填图采样网格大小约为 10 km×10 km。根据地球化学制图要求, 按图纸上每 1 cm² 至少有 1 个数据的原则, 换算相应的采样密度。为便于网格划分和点位部署， 实际采用经度 7.5′、 纬度 5′为一个采样单元格进行布点。该项目按中国地质调查局发布的标准开展工作。采样点布置主要借助遥感影像图进行部署。布点时主要考虑采样点宜均匀分布在整个调查区域，控制采样网格内的最大汇水域，当单一采样点所代表的汇水域无法控制 2/3 网格面积时，在采样网格内增加采样点。老挝地理位置处于中南半岛，境内多山，低山丘陵次之，平原低谷较少，海拔总体较高，采样介质以水系沉积物为主。野外采样时，以手持式全球卫星导航定位系统实时实地进行定位，每个样品由 3 个子样组成，原则上以横切河床断面， 左、 中、 右采集 3 个子样组合成一个样品。样品粒级小于 2 mm(– 10 目)。样品采集重量过筛后大于 1 kg。重复样采集比例达到总样品数量 3%。包括重复样，老挝全境共采集样品 1905 件。

老挝国家尺度地球化学填图项目总共分析 69 种元素。样品测试分析在河南岩石矿物测试中心完成。分析方法以电感耦合等离子体质谱（ICP-MS）、电感耦合等离子体发射光谱仪（ICP-OES） 和 X 射线荧光光谱分析法（XRF）为主。所有元素分析检出限均低于国家尺度地球化学填图对沉积物或土壤样品分析测试的最低要求。分析方法的准确度和精密度均采用分析国家一级地球化学标准物质或国际地球化学标准物质（根据样品类型选择水系沉积物或土壤）的方法进行检验。所有元素的报出率均≥ 90%。所有元素准确度和精密度合格率均达到 100%。铜的检测方法为 XRF， 检出限为 1μg/g。

4.1 典型铜矿床

（1）富开（Phu Khan）铜金矿

富开铜金矿位于万象省赛宋奔(文)县中部。大地构造位置位于印支地块东部的长山构造岩浆岩带之长山成矿带北段，区域内次级断裂和皱褶比较发育，早古生代至新生带中酸性岩浆侵入活动强烈。长山构造岩浆岩带，为铜金铁成矿作用提供了热源和成矿物质，是长山多金属矿带最重要成矿因素之一。上石炭统—二叠系以碳酸盐岩为主， 夹火山岩和碎屑岩， 是该矿区铜金矿主要赋矿围岩。该矿床类型复杂，斑岩型、矽卡岩型、热液型铜矿均有发现。其中，斑岩型铜矿体产于花岗闪长斑岩、石英二长斑岩及流纹英安质火山杂岩中，火山杂岩具同心环状构造（图 2）。含矿母岩发生强烈蚀变，具钾化、硅化、绢云母化。矿区已圈出 2 个主要矿体和 10 多个次要矿体。斑岩型铜矿含铜矿物为黄铜矿、蓝铜矿、斑铜矿及孔雀石。矽卡岩型铜矿体产于花岗斑岩及细粒花岗岩的外接触带内，矿体受构造控制明显，呈透镜状、似层状产出。热液型铜矿体产于构造破碎带中，含矿母岩为火山碎屑岩。矿体呈透镜状、脉状产出。从矿区岩浆岩的活动时期、区域地质资料分析，富开铜金矿的成矿时期可能为晚二叠世—中三叠世。

（ 2）色（塞）奔（ Sepon）铜金矿

色奔铜金矿位于老挝中南部沙湾拿吉省色奔镇。大地构造位置位于印支地块长山皱褶带南东端色奔盆地内。地质历史上遭受了多期强烈的构造-岩浆作用，岩浆活动主要发生在早二叠世，以发育花岗闪长斑岩和少量岩脉为特征。区内出露地层以古生界碎屑岩、碳酸盐岩建造为主（图 3）。北西向和东西向断裂构造发育，交汇部位是该区矿化最有利部位。矿区内矿体总体呈东西向分布，铜矿化主要分布在矿区西部，金矿化分布于矿区东部。铜矿化主要与斑岩铜矿体系有关，存在三种矿化类型。铜钼铋金矿化主要发育于花岗闪长斑岩体内部，以石英网脉状矿化产出。类矽卡岩型铜矿化主要沿斑岩体与围岩接触带发育，形成块状硫化物铜矿化层和类矽卡型铜矿化层。次生铜矿化主要为后期风化淋滤作用形成，原生矿发生氧化次生富集成矿。类矽卡岩型矿化和次生铜矿化是该矿主要矿化类型。

（ 3）班康姆（ Pangkuam）铜金矿

班康姆铜金矿位于老挝沙耶武里省巴莱县。大地构造位置位于印支地块黎府褶皱带之琅勃拉邦-黎府铜金多金属成矿带的中南段。古生代—早中生代酸性火山岩系及同成分的浅成侵入岩体控制着该成矿带的矿床，岩浆岩主要为安山岩、少量斜长斑岩及花岗闪长岩。地层主要为一套由灰岩、 泥质粉砂岩、 大理岩等组成的石炭系—二叠系沉积地层。区内主要构造为北北东向断裂和轴向为北东向的背斜褶皱，其中北北东向断裂 F3 贯穿整个研究区，其为琅勃拉邦深大断裂的次级断裂，控制矿区的成矿作用并提供成矿空间，还有 F1、 F2 两条近平行断裂与 F3 组成了北北东向断裂组（图 4）。班康姆铜金矿总体呈现上金下铜的空间分布特征，属于中低温浅层热液铜金矿床，矿体围岩蚀变主要有矽卡岩化、青磐岩化、硅化及碳酸盐化。该矿成矿时期主要为晚二叠世至早三叠世。

（4）银水山（Pha Yinshui）铜矿

银水山铜矿位于南塔省孟龙县。大地构造位置位于特提斯造山区南东段。矿区内未发现侵入岩体，主要出露形成于岛弧环境的中—基性火山岩。地层主要由二叠系陆相—浅海相沉积与火山喷发沉积，上部为陆相沉积，下部为火山沉积序列。矿区构造较发育，主要表现一个背斜构造，主要发育近南北向 F1 断层及北西、北东向的 F2、 F3 断层。矿体主要赋存于南北向产出的中-基性火山岩中，受地层岩性（细碧岩）及构造（南北向与东西向断裂）的控制，矿体呈似层状、透镜状沿断裂产出。该矿床成因上属于海底中一基性火山热液喷流沉积、后期热液叠加富集改造型铜矿床。

（5）新寨（Xinzhai）铜矿

新寨铜矿位于丰沙里省北部。大地构造位置位于兰坪-思茅盆地南部。矿区为沉积岩分布区，以中生界侏罗系、白垩系为主，属典型的陆相碎屑沉积。矿区内构造发育，主要表现为一个向斜构造，发育有大量断裂，其中北西向断裂贯穿整个矿区。矿体由两部分组成，包括层状矿体和脉状矿体。层状矿体主要赋存在侏罗系灰白色砂岩和紫色砂岩带中，呈层状或似层状，矿石矿物主要为黄铜矿、孔雀石等。脉状矿体沿断层分布，产于白垩系景星组石英砂岩中。该砂岩型铜矿床在形成过程中存在强烈的还原流体蚀变作用，成矿时期可能是喜马拉雅早期。

（6）班纳莫（Ban Namo）铜矿

班纳莫铜矿位于乌多姆赛省纳莫县。区域构造上位于思茅-丰沙里微陆块。矿区出露地层为上三叠统，属红色陆源碎屑沉积。矿区内断层发育，地层多呈断块状产出。沿断层破碎带，发育了数米至数十米宽的构造破碎带。矿体呈脉状、透镜状产于构造破碎带中。班纳莫铜矿类型为热液型，含矿热液源于沉积盆地中的卤水，成矿时期可能为古近纪晚期。

（7）南永（Nanyong）铜矿

南永铜矿位于老挝丰沙里省西部，大地构造位置处于昌都-思茅-南邦陆块中之思茅-丰沙里微陆块，属于云南“三江”铜多金属巨型成矿带的南延部分。矿区出露地层下白垩统景星组，属灰白色中细粒石英砂岩，属于河流-湖泊沉积环境。矿区断裂构造较为发育，其中北东向断裂是本区的导矿构造和容矿构造。构造破碎带对矿体的控制作用十分明显，矿体一般呈脉状、似层状、透镜状沿碎裂岩构造角砾岩带产出沿破碎带发育，赋存于下白垩统景星组中。矿床为沉积改造型之中低温热液矿床。

4.2 主要成因类型和时空分布特征

老挝主要铜矿床受岩浆热液及构造条件所控制，类型有斑岩型、矽卡岩型、热液型铜矿，砂岩型铜矿也是一个重要类型。斑岩型铜金矿，主要分布在沙耶武里、万象省和川圹-沙湾拿吉-阿速坡一带，矿床规模很大，一般与金共生，有较大的找矿前景。热液型和矽卡岩型铜矿分布范围广，主要分布在川圹省、沙耶武里省、万象省，矿体分布于中酸性侵入岩的外接触带，通常呈透镜状和脉状产出，矿石品位较高，常与铅锌伴生。砂岩型铜矿主要分布在老挝北部的南塔、乌多姆赛和丰沙里等省以及南部的占巴色、阿速坡一带，南部普遍含金，矿体呈层状、似层状产出，矿石品位较低，但矿床规模较大。

老挝铜矿成矿时代主要集中在晚石炭世-早三叠世和侏罗纪-早白垩世两个主要成矿期。老挝自早古生代开始进入板块活动演化阶段，主要表现为区内几个小洋盆的发生、发展和消亡及与其间地块的洋陆相互作用，相应的成矿作用也受制于洋陆转换体制的演化。在演化过程中，在不同构造部位发育不同的构造-岩浆-沉积-成矿作用，出现多成矿期、多矿源、多成因和叠加成矿等有利条件，使这个时期的矿产种类、成因类型不断增多，尤其是晚古生代-早中生代出现了许多大型-超大型矿床。老挝铜矿在该时期主要分布于长山成矿带和琅勃拉邦—黎府成矿带， 形成了一批大中型的斑岩-矽卡岩型矿床， 典型代表如富开铜金矿。晚石炭—中二叠世形成北西—南东向的火山弧， 洋盆在中二叠世闭合。俯冲阶段岩浆活动在长山带南西侧的响应较早，首先形成富开埃达克质岩石及一系列流纹质英安斑岩、 英安岩、 安山岩， 以及花岗岩、 花岗闪长岩等。在晚二叠世—早三叠世， 幔源岩浆出溶作用形成初始富硫富铜金矿质的岩浆流体，在其上侵过程中与安山质围岩发生普遍的水岩交互作用，然后与雨水型地下水混合形成混合成矿流体，成矿流体沿断裂构造体系上升，广泛发生成矿作用，典型代表如班康姆铜矿。在经历板块俯冲、伸展作用、陆陆碰撞作用和陆内造山作用后，中生代中期老挝境内铜成矿作用主要与陆内沉积盆地有关。在这一时期，盆地内主要发育下侏罗统-下白垩统河流-湖泊相陆源碎屑红色砂岩建造，是砂岩型铜矿床的主要赋矿地层，矿床成因为同生沉积（-后期热液叠加改造）成矿，铜矿化多产于红色碎屑岩所夹的浅色砂页岩中，典型代表如新寨铜矿。

4.3 成矿要素

老挝铜矿床的形成明显受到构造-岩浆-沉积作用控制。综合老挝境内区域地质背景及典型铜矿床的成矿规律和矿化特征等要素，按其在成矿过程中的不同作用，对老挝境内的主要成矿要素进行总结（表 1）。

5.1 铜地球化学特征

对老挝全境 1905 件样品铜元素分析结果进行参数统计（表 2）。老挝全国水系沉积物中铜的含量在 1.20~459 μg/g 之间，平均值为 21.96 μg/g，中位值为 16.50 μg/g。平均值与我国南方湿润低山丘陵区水系沉积物平均值 21 μg/g 和热带雨林区水系沉积物平均值 20 μg/g[86]相当。老挝境内 7 个三级大地构造单元中, 长山地块和哀牢山—马江等三个缝合带的平均值大于 23 μg/g，中位值大于 19 μg/g，均高于其他几个构造单元。哀牢山—马江缝合带平均值最大，达 25.53μg/g，其次为长山地块和奠边府—黎府缝合带，最小的为景洪—素可泰火山弧带。奠边府—黎府缝合带中位值最大，达 19.72 μg/g，其次为长山地块和色潘—三岐缝合带，最小的为万象—昆嵩地块。

采用 Geochem Studio 软件制作地球化学图。网格化处理时网格单元间距 Dx, Dy 均为10 km，计算模型采用指数加权，数据搜索模式采用圆域搜索，半径为网格单元间距 2.5 倍，采用累积频率的分级方法按累积频率 2.5%, 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%,60%, 70%, 75%, 80%, 85%, 90%, 95%, 97.5%所对应的含量划分成 18 级，制作形成铜地球化学图。

从老挝铜地球化学图（图 5）中可以看出，铜在老挝分布并不均匀，存在多个大面积分布的高背景区和异常区。这些区域主要位于长山地块琅勃拉邦省—川圹省—华潘省，思茅—彭世洛地块丰沙里省和沙耶武里省南部区域，万象—昆嵩地块占巴塞省—塞公省—阿速坡省，此外在这三个地块局部地区零星存在小面积的异常区。铜的低背景区和低值区主要分布于万象—昆嵩地块沙拉湾省和沙湾拿吉省，以及占巴塞省的西南部地区。从老挝铜地球化学图和老挝已知铜矿床（点）的分布可以看出，部分铜矿床（点）位于高背景区和异常区范围内，如色奔铜金矿。部分铜矿床（点）位于背景区范围内，如班纳莫铜矿。此外，受地质高背景影响，有些区域铜元素含量高，面积大，但尚未发现成规模的铜矿床，如位于巴塞省—阿速坡省内的波罗芬高原玄武岩区。

5.2 地球化学异常信息提取

地球化学数据具有成分数据特性，通过对数比变换可以“打开”原始数据，使每个元素的含量可以独立分析，而不受总量限制的影响。这种方法有助于更深入地理解成分数据中各元素之间的关系，为后续的数据处理和解释奠定基础。

建立在对数比变化方法之上的连续二值分解（Sequential Binary Partition, SBP）方法是一种在单形空间中利用正交基构建对抗关系的方法。通过将需要研究的变量分成两组（一组为“+”， 另一组为“-”）， 然后对这两组变量分别构建对抗关系， 不断重复这个过程， 直到只剩下两个单独的变量为止。这种方法作为一种数据降维技术，能够帮助聚焦于感兴趣的组合，并研究它们所携带的相对信息。

在第 i次划分中， 标记为“+”的变量个数 r_i+与标记为“-”的变量个数si-
的几何均值分别为gi（x₊）与gi（x_-） ，通过相应的数学运算得到相应的成分矢量Z_i。

将经过成分数据变换后的地球化学数据进行主成分分析，得到第一个主成分（PC1）和第二个主成分（PC2）的 Biplot 图（图 6）。

在 biplot 图中，第一象限的元素与中酸性岩浆岩有关。第二象限出现的元素可能与地层分布有关。第三象限出现多种铁族元素，推测与基性、超基性有关。此外，第三象限出现的Zn、 Hg、 Mo 等元素、第四象限的 As、 Sb、 Pb 等元素以及 W、 Mo、 Bi 作为典型的高温元素，与K₂O 出现在第四象限可能暗示了与这些元素成矿与中酸性岩浆热液活动有关。从数据驱动结果来看， U、 Th 元素的相关性较高，符合中酸性岩浆岩中 U、 Th 分布相似的事实。Co、 Ni、 V、 Fe₂O₃、 Cu、 Cr 元素相关性较高，表明存在部分玄武岩出露。As、 Sb、 Pb、 W、Mo、 Bi、 Sn 等元素存在一定的相关性，一定程度上表明老挝严格受岩浆热液控制。

本文通过对老挝铜矿成矿要素的详细梳理，认为中酸性岩浆活动对于区域铜矿化具有重要影响，为铜多金属矿的形成提供了物源，特别是斑岩型铜矿体主要产于花岗闪长岩、石英二长斑岩中，矽卡岩型和热液型铜矿产于中酸性侵入体的外接触带中。老挝深部断裂和岛弧的发育导致岩浆侵入和火山活动的多阶段性，分布广泛。依据上述认识，通过 SBP 技术，应用成分数据分析中的连续二值分解功能划分出 2 种元素组合（表 3）。

老挝岩浆活动分布比较广泛，中酸性侵入岩的控制范围对矿床具有一定的指示意义。根据 B1（图 7）元素组合图发现，中酸性岩体分布的区域基本都对应着图中的高值区域。但在未有中酸性岩体分布的沙湾拿吉省存在大面积异常，根据地形和水系分布，初步分析是因该处地势低洼，区域物源主体来自于东部特别是越南境内的花岗岩体所致。整体来说， B1能够有效的识别中酸性岩体的分布。

与铜矿相关的成矿元素主要包括 Cu、 Mo、 Pb、 Zn、 Au、 Ag、 W、 Sn 等，从 B2 元素组合图来看（图 7），正值 0.72 以上的异常范围与区内已发现的矿床都比较吻合，无论是斑岩矿床、矽卡岩矿床还是热液脉型矿床。这说明了利用成矿元素的组合异常来圈定矿集区矿化的范围是有效的。

6.1 定量矿产预测模型建立

以地质成矿模式为核心，结合地质、地球化学等技术开展成矿有利信息定量提取是实现成矿预测工作的基础。前文通过对地球化学数据的处理，提取了两种不同的元素组合，构建了两个不同的成矿预测指标。基于此，本节进一步采用 GIS 空间分析对区域地层、构造等数据开展预测指标定量化提取，从而构建包含地质、地球化学的定量矿产预测模型。

① 断裂构造

断裂构造作为含矿热液流体运输的通道，是岩浆矿床、矽卡岩矿床、热液矿床等在内的内生矿床成矿的关键影响因素。依据矿体在断裂构造旁侧就位的规律，断裂构造距离分析可以作为一种定量提取构造信息的有效方法，经 GIS 空间分析，形成断裂构造距离分析图（图8a）。

② 碳酸盐岩和碎屑岩

根老挝铜成矿规律，很多铜矿控矿沉积建造均与碳酸盐岩、碎屑岩有关，特别是矽卡岩型的铜矿。因此，将碳酸盐岩和碎屑岩主要分布区距离作为成矿有利指标之一，经 GIS空间分析，形成碳酸盐岩和碎屑岩距离分析图（图 8b）。

根据老挝区域地质背景，通过成分数据分析完成了地球化学数据成矿有利信息定量提取，并结合上述工作，依托断裂、碳酸盐岩和碎屑岩分布等地质信息共同构建老挝铜矿多源信息定量信息预测模型。如表 4 所示：

6.2 基于随机森林的成矿预测

6.2.1 理论概述

随机森林方法是一种非参数集成机器学习方法，由一系列相互独立的决策树组成。该方法通过结合多个决策树的预测结果，利用投票或平均方式来获得最优结果。随机森林具有较高的准确性和鲁棒性，尤其在处理高维数据和防止过拟合方面表现突出。

随机森林的构建过程主要包括两个步骤。首先，采用 bootstrap 方法从原始数据集中随机选择 67%的样本作为训练样本用于构建每个决策树，其余 33%的样本（袋外数据）用于内部交叉验证，以估计模型的泛化误差。其次，在构建每棵决策树时，随机森林方法在每个节点随机选择一部分特征。具体而言，对于每个节点，从全部特征中随机选择一个子集，然后在这个子集中找到最佳分裂条件，从而确定该节点的分裂方式（图 9）。在节点处，使用基尼指数（Gini impurity）来评估特征的不纯度。基尼指数是一种衡量分类问题中节点纯度的指标，该值越低，表示节点中的样本大多属于同一类别。具体公式如下：

6.2.2 训练过程

采用 10km×10km 的网格将研究区划分为 2432 个预测单元，将 12 个预测变量的数据赋值到预测单元中，将数据网格化的同时保留数据的地质属性，最终形成的数据表包含 2432行数据（即预测单元个数）， 12 列属性（即定量化指标个数）。

基于随机森林的成矿预测目标变量是一个二值变量，选择合适的训练样本是成功预测关键之一。研究区已知矿床/点为正样本，数量为 131 个。相对于矿化来说，非矿化是一种普遍的现象。从确定的、远离矿化的区域，远离构造区域、在主成矿元素含量较低的区域，随机选取 119 个负样本。因此，构建了一个正负样本比例约为 1:1 的数据集，将 70%的数据为训练集， 30%的数据为测试集。

随机森林模型是基于 Python 中 sklearn 包实现的。在模型训练过程中主要对以下参数进行参设置。

①树的数量（n_estimators）：基评估器的数量，对模型的精确性有单调影响。通常，随着 n_estimators 的增加，模型的效果会不断提升。然而，每个模型都有其决策边界，当n_estimators 达到某个临界点后，随机森林的精确性可能不再提高，甚至会开始波动。此外， n_estimators 越大，所需的计算资源和内存也越多， 训练时间也会相应增加。因此，需要在训练难度和模型效果之间找到平衡。

② 树的最大深度（max_depth）。一般来说，当数据量或特征较少时，可以使用默认设置忽略该参数。但在数据量大且特征多的情况下，需要限制最大深度，其具体取值应根据数据的分布特征来确定。

③ 树的特征个数（max_features）。当特征数量较多（超过 50 个）时，通常选择特征总数 n 的平方根作为树的特征数。反之，如果特征数量较少，则通常会使用全部特征。

本文使用 Python 中 sklearn 库的 GridSearchCV 模块进行参数优化，该模块能够遍历所有可能的参数组合，以确定最佳参数值。通过分析三个主要参数，并考虑到本文使用的样本量较少且研究的是分类问题，因此基于网格搜索法设定最优的决策树数量为 150，树的特征数量和最大深度保持默认值，采样方法选择有放回的采样法。训练后发现，虽然其他参数对模型训练速度有一定影响，但对准确率影响不大，因此这些参数采用默认设置。最终，所构建的随机森林模型在测试集上的预测准确率为 0.85， AUC 值为 0.9（图 10）。

6.3 成矿远景区的圈定

根据训练集数据所建立的随机森林模型，对老挝进行了全区域预测，并对预测结果进行了可视化。结合区域地质背景，成矿地质条件，共圈定了 9 个成矿远景区（图 11）。

远景区 I 和 II 位于丰沙里-会嘎金、铜、铅锌成矿带上，与我国三江成矿带相连。该区域断裂构造发育，出露地层主要包括中新生代的海相和陆相沉积岩以及碳酸盐岩，富含金、铜、铅锌等矿床。新寨铜矿、南永铜矿等均位于该远景区内。基于区域地质情况和已发现矿床的类型，该远景区具有找寻砂岩型铜矿或沉积改造型铜矿的潜力。

远景区Ⅲ位于琅勃拉邦成矿带， 其周边典型矿床为帕奔金矿， 区内断裂以 NE 向为主，韧性剪切带较为发育，岩浆活动频繁，这些对铜金多金属矿的形成极为有利。元素组合 B1、B2 均在此处有异常分布，结合随机森林的预测结果，此成矿远景区具有寻找斑岩型、矽卡岩型铜矿的潜力。

远景区Ⅳ横跨万象省、 沙耶武里省、 琅勃拉邦省和川圹省， 位于琅勃拉邦-黎府成矿带、川圹-长山成矿带与万象-昆嵩成矿带的交汇处。该区拥有多个典型矿床，如富开斑岩型铜金矿床和纳勐铜矿。区内出露地层主要包括石炭系、二叠系碳酸盐岩和侏罗系海陆交互相碎屑。岩浆活动、构造特征和岩性条件共同为斑岩型铜金矿床的成矿提供了有利环境。主要的蚀变类型包括矽卡岩化、硅化和碳酸盐化，这些蚀变与铜多金属矿化相关。元素组合B2 在此处呈现出高异常，结合随机森林的预测结果认为该区域具有较大找矿潜力。

远景区Ⅴ位于川圹省北部， 主要由 NW-SE 走向的控矿断裂构造控制， 区内发育有花岗闪长斑岩和石英二长斑岩等岩体。中酸性岩浆的侵入和喷发活动为本区金属矿床的形成提供了丰富的物源，复杂的断裂构造则为成矿过程提供了必要的构造条件，利于斑岩型铜矿的形成。Cu 元素在此处呈现出高异常，与矿化点对应较好。同时，随机森林模型将此处预测为高概率区域，是斑岩型铜矿勘查的极为有利区域。

远景区Ⅵ和Ⅶ位于华藩省，属于桑怒成矿带，为岛弧成矿带。该区内中酸性岩浆活动频繁，构造复杂，为铁、铜、铅、锌等多金属矿床的形成提供了丰富的成矿热液。远景区内马江断裂具有切割深度大、多期活动等特征，由该断裂构造运动派生的一系列 NW 向构造构成了区域性重要的控岩控矿构造，区域内发育镁铁质-超镁铁质岩体，已发现多个铜镍矿点，显示出找寻铜镍硫化物矿床的找矿前景。

远景区Ⅷ位于沙湾拿吉省， 典型矿床包括色奔铜金矿。区内主要出露中石炭系至二叠系地层，中酸性岩浆岩发育， Cu 元素与元素组合 B2 都在此处呈现出高异常，基于随机森林预测结果，该区应为找寻斑岩-矽卡岩铜多金属矿床的潜在区域。

远景区Ⅸ位于老挝最南部， 塞公省和阿速坡省东部区域。区内深大断裂发育， 岩浆活动频繁， Au、 Pb、 Zn 等元素异常强度高。区内爆发的大规模的中酸性岩浆岩侵入与喷发活动为形成斑岩型铜金矿床提供了良好的成矿条件。该区存在大量砂金矿点，为寻找原生铜金矿床提供了线索。

本文主要利用老挝 1∶100 万国家尺度地球化学填图项目取得的高质量基础数据， 在总结老挝铜矿成矿规律基础上，利用机器学习技术，开展铜资源远景区预测，结论如下：

（1）老挝铜矿床主要类型有斑岩型、矽卡岩型、热液型以及砂岩型。铜矿床的形成明显受到构造-岩浆-沉积作用控制。矿床主要成矿期为晚石炭世-早三叠世、侏罗纪-早白垩世；控矿构造主要为区域性南北向、北东向、北西向断裂；控矿侵入岩主要为古生代-早中生代酸性侵入岩、中性闪长玢岩、花岗闪长斑岩；控矿沉积建造包括二叠系陆相-浅海相碎屑沉积和碳酸盐岩、白垩系陆源碎屑沉积。

（2） 老挝全国水系沉积物中铜的含量在 1.20~459 μg/g 之间，平均值为 21.96 μg/g，中位值为 16.50 μg/g。境内 7 个三级大地构造单元中, 长山地块和哀牢山—马江等三个缝合带的平均值大于 23 μg/g，中位值大于 19 μg/g，高于其他几个构造单元。地球化学图表明铜在老挝分布并不均匀，存在多个大面积分布的高背景区和异常区。

（3） 通过成分数据分析完成了地球化学数据成矿有利信息定量提取，并结合该工作，依托断裂、碳酸盐岩分布等地质信息共同构建了老挝铜矿多源信息定量信息预测模型。模型预测因子包含了单元素异常、矿化元素组合异常、指示中酸性岩体元素组合、控矿构造分布、碳酸盐岩和碎屑岩分布。

（4）利用随机森林成矿预测方法，共圈定 9 个成矿远景区。远景区 I 和 II 位于丰沙里-会嘎金、 铜、 铅锌成矿带上， 具有找寻砂岩型铜矿或沉积改造型铜矿的潜力。远景区Ⅲ位于琅勃拉邦成矿带， 远景区Ⅳ位于琅勃拉邦-黎府成矿带、 川圹-长山成矿带与万象-昆嵩成矿带的交汇处， 远景区Ⅴ位于景洪-会晒成矿带北部， 远景区Ⅷ位于沙湾拿吉省东部， 远景区Ⅸ位于塞公省和阿速坡省东部区域， 以上远景区具有寻找斑岩型、 矽卡岩型铜矿的潜力。远景区Ⅵ和Ⅶ位于华藩省， 属于桑怒成矿带， 显示出找寻铜镍硫化物矿床找矿前景。