地球化学信息挖掘在云南易门铜矿区找矿预测中的应用.doc

文章编号: 1674 9057 ( 2011) 03 0328 11地球化学信息挖掘在云南易门铜矿区找矿预测中的应用刘洪滔1 ，陈三明2，3 ，朱家群1( 1. 云南省有色地质局 地质地球物理化学勘查院，昆明 650216; 2. 中国地质大学 地球科学与资源学院，北京 100083;3. 桂林理工大学 a. 地球科学学院; b. 广西地质工程中心重点实验室，广西 桂林 541004)摘 要: 以云南易门铜矿区及其外围为研究区，对不同时代地层中 12 种微量元素的地球化学分布特征、水系沉积物中微量元素的空间分布特征以及地球化学指标与成矿的关系进行了分析。依据地球化 学信息挖掘理论，采用微量元素主成分因子 泛克里格模型、三维化探异常可视化分析，以及基于 GIS 的多元信息预测等方法，对地学数据进行了信息的挖掘与提取，建立了基于 GIS 的地球化学综合 找矿模型。通过对地球化学弱异常信息增强和基于地球化学综合指标的成矿预测，圈定了 9 个找矿远 景区及 17 个找矿靶区，取得了较好的找矿效果。关键词: 地球化学; 信息挖掘; 找矿预测; 主成分因子 泛克里格模型; 易门铜矿; 云南中图分类号: P628文献标志码: A地 球空间数据挖掘 ( spatialdata mining，以突现，本文综合运用地球空间数据挖掘的理论和方法，在寻找元素地球化学场与成矿的空间关 系上取得了新的认识。SDM ) ， 又 称 知 识 发 现( knowledge discovery，KD) ，是指从地球空间数据库中发现隐含的、未知的、潜在有用的信息，提取感兴趣的空间模式 与特征、空间与非空间数据的普遍关系及其他一 些隐含在数据库中的数据特征，目的是把大量的 原始数据转换成有价值的知识1 2。地学可视化 是科学计算可视化与地球科学结合而形成的概念， 是关于地学数据的视觉表达与分析3，微量元素 的主成分分析和三维可视化及基于 GIS 的多元信 息预测均是地学数据挖掘的方 法4。 自 阮 天 健 等5提出地球化学场的概念以来，如何将同一点 上各种地质、地球化学参数之间的动态关系有效区域地质成矿背景云南易门地区大地构造位置位于扬子准地台 西南缘，川滇台背斜南部7; 昆阳裂谷武定 易 门 元江裂陷槽的中部8 ( 图 1) 。云南易门地区 构造复杂，变质作用广泛，成矿流体活跃，成矿 作用复杂多样，是云南省著名的铜矿集中区之一。 迄今已探 ( 查) 明大型铜矿床 1 处 ( 易门三家厂 铜矿) ，中型铜矿床 2 处 ( 易门狮子山铜矿、易门 铜厂铜矿) ，小型规模 10 余处。1地提取挖掘并直观地表达出来运用于找矿工作中， 1. 1成矿构造一直以来均是化探工作者探索的热点6。结合近年来在云南易门矿区开展的1 5万水系沉积物地球 化学测量工作，运用传统的异常提取方法难以准 确把握异常的分带规律，隐伏矿床的弱异常也难区内的构造环境主要是扬子地台西缘昆阳裂谷系，受裂谷作用的影响，区内的铜铁矿床多受 近南北向与东西向两组断裂控制，尤其是褶皱及 交 叉或网状断裂处往往是成矿的有利地段。而绿收稿日期: 2010 12 23基金项目: 中国地质调查局战略性矿产资源远景调查项目 ( 1212010880803) ; 广西地质工程中心重点实验室主任基金 ( 桂科 能 07109011 Z013) ; 广西研究生教育创新计划项目 ( 2010105960818M30)作者简介: 刘洪滔 ( 1969) ，男，硕士，高级工程师，研究方向: 勘查地球化学，lym0813 tom. com。引文格式: 刘洪滔，陈三明，朱家群 地球化学信息挖掘在云南易门铜矿区找矿预测中的应用 J 桂林理工大学学报，图 1 云南易门铜矿区区域构造纲要图8Fig. 1 Sketch structure map of Yimen copper ore district in Yunnan汁江断裂是区内最主要的铜矿导矿构造。区内断裂主要有北北东向、北西向和北东向近东西向 3 组。北东向断裂组控制了全区地层的展布形态和 矿体的产出形态，属区内主要的控矿容矿构造; 北西向组断裂多错断北北东向组断裂，该断裂组 主要对北东东向近东西向组构造和地质体起到 破坏作用。受绿汁江导矿断裂长期活动及东西向挤压造 山作用等影响，在南北向和北北西向构造带内或 其交汇部位，分布着大量的紧密褶皱，还有与褶 皱方向基本一致的相配套的北东向、东西向断裂 构造，以北西走向及北东走向最发育，将本区切 割成大小不等的地质块体 ( 图 2) 。区内与铜矿有面积最广，铁铜矿赋存的层位有中元古界绿汁江组 ( Pt2 lz) 、鹅 头 厂 组 ( Pt2 e ) 、落 雪 组( Pt2 l ) 、因民 组 ( Pt2 y ) 、 大 龙 口 组 ( Pt2 d ) 、 富 良 棚 组( Pt2 f) 、黑山头组 ( Pt2 hs) 。其中: 铜矿赋矿层位 主要有落雪组、鹅头厂组、因民组，以及绿汁江 组下部与因民组角砾岩断裂接触带; 铁矿主要赋 存于大龙口组、富良棚组、黑山头组 ( 表 1) 。1. 3 岩浆岩裂谷早期的岩浆活动以基性为主，仅有少量 中酸性，且多位于含矿层之下。“万宝厂式” 铜钴 矿明显受海相火山 喷流沉积建造所控制。铜钴 矿层下盘为中基性 基性为主的钠质喷出岩 复 杂角砾岩带，复杂角砾岩主体具一定层位，呈多 层产出，为一套由岩浆活动引起的火山泥石流、 同生断裂生长发育及其相伴的水下地震 海底崩 裂 陡壁塌陷形成的同生至准同生角砾岩带。研 究表明，昆阳裂谷中铜铁矿床多产于封闭较好、 水体较深的断陷盆地内。其聚矿场所，铁矿主要在关的构造带主要有狮子山 绿宝冲 ( 青龙山)皱带、铜厂褶皱带和南部的甸中褶皱带。褶赋矿地层9 12区内地层新生界、中生界、新元古界、中元 古界均有出露 ( 图 2) ，以中元古界、中生界出露1. 2变化特征。根据地表铁锰帽 ( 矿石) 元素组合特征估计矿体剥蚀深度，评价深部矿体远景。如果 地表样品的地球化学元素 As、Ag 相对较高，而 Pb、Zn 较低，说明矿体剥蚀深度不大，深部可能 存在工业矿体。反之，矿体剥蚀深度较大，深部 有矿的可能性较小13。为系统研究原生地球化学场中元素在不同时 代地层中的分配特征，化探普查工作中系统采集 了易门矿区内出露地层中未矿化地段 ( 背景区) 不同岩性岩石样品，分析测试了 Cu、Pb、Zn、Ag 等 12 个元素，统计结果( 表 2) 表明: ( 1) 不同时 代不同岩性的岩石微量元素含量差异较大。( 2) 中 元古界绿汁江组白云岩，鹅头厂组碳质板岩、粉 砂岩、板岩，落雪组白云岩，因民组 ( 角) 砾岩， 美党组板岩、( 粉) 砂岩，黑山头组石英砂岩、板 岩中 Cu 含量远高于地壳同类岩石中的平均含量， 其中尤以落雪组白云岩 Cu 含量最为显著。反映这 些层位在成岩时期伴随沉淀了丰富的 Cu 元素，为 后期成矿提供了稳定的物质来源，故 Cu 是区内的 优势矿种。 ( 3) 与地壳同类岩石中的平均含量相 比较，Co、Ni 在不同时代地层中均明显高度富集， 其中以 Co 富集程度最高，这一特征与区内存在多 个钴矿，且铜矿中多伴生钴相吻合，Co 是区内的 优势矿种也说明了这一点。 ( 4) 区内碳酸盐岩中 的 Mn 含量普遍较高，与区内铁锰矿化在碳酸盐岩 中广泛发育现象相吻合，反映了 Fe、Mn 也属区内 的优势矿种。 ( 5) 与地壳同类岩石中的平均含量 相比较，Pb、Zn 在区内多个层位中均相对富集，图 2 易门铜矿区地质矿产简图Fig. 2 Geological and mineral sketch of Yimen copper ore district反映区内具有一定的铅锌矿的找矿前景。 ( 6)Ag在上禄丰群泥岩、澄江组砂岩、因民组 ( 角) 砾岩中相对富集。( 7) Au 除在因民组板岩中高于地 壳丰度外，其余地层中含量均不高，区内至今也火山活动形成的穹丘内及其附近的中基性熔岩相中，而铜铁矿多在火山穹丘间的沉积洼地内的火 山喷流 热水沉积岩相组合中8。未发现有工业价值的金矿。 ( 8)反映了取样地段构造发育程度。As、Sb、Hg、Bi微量元素在不同年代地层中的分2布特征地球化学元素的组合是隐伏矿体预测的地球 化学标志之一。研究表明，热液矿床成矿元素及 成矿伴生元素沿轴向有明显的分带，分带序列自水系沉积物中微量元素的分布特征33. 1元素地球化学的主成分因子 泛克里格模型在元素地球化学空间分析中，一个非常重要 的问题是区域异常识别提取，而空间求异是一个 极富价值的理论与实用性问题，求异过程的基本 原则是: 减少空间数据相关性、离散空间数据结 构性、寻找空间目标受控条件多元性差异14 。其下 而 上 为:MnPbZnCu ( As、 Ag ) 。 As、Ag 主要富集在矿体的顶部，Pb、Zn 富集在矿体的底部。而且元素特征比值如( As + Ag) / ( Pb +Zn) 、 ( AsAg) / ( PbZn) 等具有明显的单向表 1 易门铜矿区昆阳群主要地层特征9 12Table 1 Major stratigraphic characteristics of Kunyang Group in Yimen copper ore district地 层 单 元厚度 /m主要岩石类型火山岩典型金属矿床界系统组代号柳坝塘组Pt2 lb砾岩、铁质砾岩、板岩、硅质板岩256脉状铅锌 矿、三 家厂( 狮、凤山) 式铜 矿、峨腊厂铜矿上部泥质白云岩夹板岩; 中部藻礁白云岩; 下部角砾状白云岩、白云岩夹板岩绿汁江组Pt2 lz1 321 2 030上部黑色板岩夹灰岩; 中部板岩夹砂岩;下部碳质板岩，局部有铜、铁矿体鹅头厂组Pt2 e993 1 796顶部碳泥质白云岩夹粗晶灰岩; 中部柱状叠层石白云岩; 下部藻席白云岩、同生角砾 岩; 底部火山质泥质硅质白云岩、硅质岩流纹岩、粗面火山岩东川式铜矿( 狮子山、铜厂、老厂)Pt2 l落雪组50 540中前昆中上部砂泥质白云岩夹紫色板岩、砂岩;中下部火山角砾岩、凝灰岩、细碧岩。产 铁铜矿细碧 角 斑 岩及火山碎屑岩稀 矿 山 式 铁 铜 矿( 万宝厂)元震因民组Pt y30 8182阳古旦灰绿色板岩，顶部夹粉砂岩及钙质砾岩;底部夹泥灰岩钙质板岩美党组Pt2 m1 448 5 000界系群薄层灰岩蠕条灰岩叠层石白云岩，下部砂屑灰岩夹钙质石英砂岩、塌积岩，下部及 中上部为含铁层位鲁奎山式 铁 矿、大六龙 他 达 式 铁 矿、化念铁矿Pt2 d大龙口组551 3 538枕状细碧岩顶部为气孔杏仁状安山岩、安山凝灰岩、沉凝灰岩与泥灰岩及灰岩互层，中下部为 变质石英砂岩、石英岩夹板岩气 孔 杏 仁 状 安山岩、凝灰岩贡山铁矿、小 假 佐铁矿富良棚组Pt2 f105 240上部为含铁变质石英砂岩、石英岩夹板岩; 中下部夹铁质砂板岩尖山箐铁 矿、天 星厂铁矿黑山头组Pt2 hs765 4 211中的因子分析是一种去相关分析，即将原始地球化学数据看作为 N 维空间中的向量，通过坐标变 换并取得各主成分在变换后的坐标方向上的方差 极大化15。因此，地球化学元素进行主成分分析 可将特定的信息集中分布到某几个主成分因子分 量中。在云南易门铜矿区元素地球化学建模过程 中综合借鉴了因子分析与泛克里格法的原理，构 建了如下的主成分因子 泛克里格模型。地球化学勘查通常将测点数据表达成规则网 格 ( 栅格) 点阵形式，设在某区域获取网格数据 X，表示为数据集X = ( xij ) n m ，n 为样品数，m 为地球化学元素变量数( x1 ，x2 ，xm ) 。按因子分析原始模型p的 m 个变量简化成少数的因子，即 p m; 系数 akj和 j ( j = 1，2，m; k = 1，2，p) 叫做因子载荷，akj 是公因子负载，j 是单因子负载，由于本文关心的只是公因子，通常所说的因子载荷仅指 a。kj求其旋转变换解( 方差极大旋转或四次幂极小旋转)x1x2 xmf1b11b21b12b22b1m b2m f2B = bp mfbbp1p2pmpF = ( f ，f ，f ) ，则 F = XR 1 B( R 1 为 X 之令1 2p相关阵逆阵) 。即F = ( f )，F 为因子计量阵。kj n p= akj fk + j j ，j= 1，2，m;xj然后对原始网格数据 X 重新构造，目的在于对变量进行必要的线性组合，即构造所谓主成分，于 是有mk = 1其中: f1 ，f2 ，fp 和 1 ，2 ，m 都是理论上的待求新变量( 因子) ，f1 ，f2 ，fp 叫做公因子，反映各 变量间的相关信息; 1 ，2 ，m 因子仅出现于一个变量的表达式中，因此称单因子，反映相应变量 的特征信息。正整数 p 是公因子个数，意思是把原来( 1)ikXBT= X( 1) 。 ij kj或x b= xj = 1将 X( 1)中的元素 x( 1) 乘以相应因子计量，得ik表 2 微量元素在不同时代地层 (岩石) 中的分配特征Table 2Distribution features of microelements in different stratums ( rock)Z( x) = m( x)+ R( x) 。mm x( 1) x( 1) 1 bpj f1p1j b1j f11通过求数学期望 m( x) 可得 R( x) ，即为区域异常。当然这里假定 Z( x) 存在变异函数，其增量具有 一阶矩与二阶矩，在无偏性与最优性条件下通过解 克里格方程组可求出 m( x) ，这里的 R( x) 分布形态 进一步可由类似于 DEM 数字高程模型以图像形式 表达出来。按云南省地球化学单元划分，易门矿区属滇中Cu、( Fe) 地球化学省( ) 之东亚区( B ) 罗茨 易 j = 1j = 1 X*=mmx( 1)x( 1)nj bnj fn11nj bpj fnp n pj = 1j = 1X*其中:中第 j 列可解释为以主因子 f 为主的空间j组合加权新变量，共 p 个。将 X* 中元素列按原始样品坐标顺序排列，构成 p 个点阵数据，每一点阵都相 当于新的网格数据，于是可依空间分布的顺序进行 求异。关于求异方法，将上述点阵数据理解成具有 漂移的区域化变量 Z( x) ，则有16门 青龙厂 Cu 异常带( B 3 )。2009 年，云南省有色地质局地质地球物理化学勘查院对云南易门铜矿区及其外围 4 个 1 5万标准分幅区域开展了1 5万地球化学水系沉积物测量，采样合计 8 459 件 样品，采样密度 4 5 点 /km2 ，测试 Cu、Pb、Mn、Ni、 Ag、Zn、Co、As、Sb、Bi、Hg、Au 等 12 个元素，数据网 格化方法采用泛克里格法，从地球化学元素的因子 分析正交因子解( 表 3) 可见，测区内不同地层单 元、构造格局、岩浆活动具有不同的地球化学分布 特征。其中: F1 为 Cu、Bi、Ag 的组合，反映区内 Cu 与 Ag、Bi 密切相关。从 F1 正交值因子异常图( 图3) 上可以明显地看到，老吾街断裂以北，Cu 元素呈 区域性富集，断裂以南则相对贫化，根据 F1 分值的 级别不同将本区划分为、 2 个一级地质 地球 化学单元( 图 3、表 4) 。F2 为 Co、Mn、Ni 的组合，该 因子反映了成岩作用中金属元素的富集及基性岩 浆岩、断裂带、第四系的分布概况( 图 4) ，利用 F2 因子将本区划分为 7 个二级地质 地球化学单元( 表 4) 。F3 为 Pb、Hg、Zn 组合，反映了受裂隙控制 的热液矿化元素分布特征( 图 5) 。F4 为 Au、As、Sb 组合，反映区内构造裂隙发育程度及 Au 异常的分 布特征( 图 6) 。其中 As、Sb 在其 F1、F3、F4 等因子 上均有较高的因子负荷，规律性不明显，反映出区 内的热液矿化与构造活动密切相关并且具有多期 次叠加的地球化学特征。从图 3、图 4 可见，区内的铜铁矿产分布于特定 的地球化学单元中，形成了多个矿化集中区。表 3 元素因子荷载矩阵(正规方差极大法)Table 3 Load matrix of elements factor(normalize-varimax method)图 3 易门铜矿区 F 1 地球化学因子异常Fig. 3 Factor F1 score in Yimen copper ore district元素F1F2F3F4CuPb Mn Ni Ag Zn Co As Sb Bi HgAu0. 758 920. 033 2960. 093 019 0. 005 970. 461 1660. 035 4050. 198 9480. 420 6200. 424 9290. 695 2710. 110 7130. 008 2360. 248 4250. 080 9050. 801 3500. 680 4270. 295 7380. 369 9630. 832 9040. 106 8580. 030 2930. 016 7660. 230 4480. 141 035 0. 123 380. 820 3370. 118 6810. 310 0530. 304 5300. 623 6380. 060 5460. 213 8350. 507 2240. 257 1720. 648 243 0. 036 34 0. 026 94 0. 094 040. 074 3490. 249 2560. 163 7470. 158 3840. 055 5340. 644 2790. 518 9480. 319 9240. 297 6930. 866 492表 4 易门铜矿区地球化学单元划分Table 4 Geochemical unit in Yimen copper ore district一 级二 级编号名称编号名 称青龙山 狮子山铜( 铁) 地球化学亚区铜厂 老厂铜( 铁) 地球化学亚区 土官乡 铁( 铜) 地球化学亚区 小尖山 三元村铁( 锰) 地球化学亚区营盘山 白家泥地球化学亚区123狮子山 易门县铜 ( 铁 ) 多 金 属 地球化学区45老吾街 峨腊铜( 铁) 地球化学区绿房 汞山铜( 铁) 地球化学亚区老吾街 甸中地球化学亚区12图 4 易门铜矿区 F 2 地球化学因子异常Fig. 4 Factor F 2 score in Yimen copper ore district图 5 易门铜矿区 F 3 地球化学因子异常Fig. 5 Factor F 3 score in Yimen copper ore district会因矿体埋深增大而逐渐减弱17 18。 虽 然 基 于GIS 的地球化学找矿模型能快速有效地圈定找矿远 景区及找矿靶区，但往往会漏掉许多有一定寻找3. 2地球化学异常分带与成矿的空间位置关系全区可分为 10 个地球化学异常带 ( 图 7) 。其中:H1 、H2 、 H3是 区 内 主要的铜矿化 带;H4 、19隐伏矿意义的低缓弱异常 。在云南易门地区，地形相对高差较大，因此 DEM 高程的突变掩盖了 低缓异常的识别，难以发现隐伏矿床表达的异常 信息。因此需要将异常信息叠加在 DEM 地形图上 再进行视觉上的三维地球化学图像分析，其目的 在于寻找到因地形突然抬升而形成的低缓弱异常。 三维地球化学图像包括: 元素含量曲面图( 图 8) 、 三维元素组合异常图( 地形图上叠覆元素组合异常H5 、H6 是区内主要的铁矿化带，H7 为铁 ( 锰)矿化带; H8 、H9 、H10 为受裂隙控制的铁矿化带。 除 H5 异常带因受北西向米茂断裂及近东西向老吾 街断裂的影响而呈东西向转北东向弧形弯曲外， 其他均为北北东或近南北向展布。三维地球化学图像可视化分析从深穿透地球化学的理念模型来看，异常均4图 7 易门铜矿区地球化学异常分带Fig. 7 Geochemical anomaly zonality in Yimen copper ore district图 6 易门铜矿区 F 4 地球化学因子异常Fig. 6 Factor F 4 score in Yimen copper ore district2. 0 ( mineral ore resources perspective and assess-ment system) 系统证据权重法模块对区内的成矿 有利区段进行预测21。云南易门矿区预测的目标 矿种是以铜铁为主的多金属矿，在对区域地质背 景深入分析的基础上，针对所预测的矿种与地质 背景的关系，把预测区划分为 1 km 1 km 的网格 基础单元，保证一个网格中至多出现一个已知矿 点，采用人工方式剔除同一网格内的多余矿点， 对不同图幅的地层分组方案进行相对统一。预测采用证据权重法预测模型，首先要根据 已知矿床 ( 点) 与各种控矿成矿条件之间的条件 概率来确定每种条件的权重值 ( 表 5 ) ， 其 中， Cu 、Ag 、Bi 衬 度 累 加 异常组合主要据表3 中 正 交图，图 9) 、三维地形地质图 ( 地形图上叠覆地质图) 等系列图像。在云南易门矿区的化探普查工 作中，通过三维地球化学图像视觉分析发现了多 个有找矿意义的低缓弱异常，例如在 A、B 两矿区 之间的 Cu 低缓异常区即是通过三维地球化学图像 分析重新圈定的新靶区，经工程验证取得了较好 的找矿效果 ( 图 8、图 9) 。5基于 GIS 的地球化学找矿预测近年来，基于 GIS 的地球化学找矿预测方法日渐普及20，本次工作采用中国地质大学数学地质遥感地质研究所基于 MAPGIS 开发的 MORPAS图 8 铜矿找矿靶区预测与低缓异常区识别Prediction of prospecting target for copper ore and characterization of weak anomaliesFig. 8图 9 地形、地质、地球化学异常三维组合Fig. 9 3D information of landform，geology and geochemical anomaly因子解 F1 确定，从图 3、表 3 可见，F1 较好地反映了区内 Cu 元素的地球化学场特征，且与已知矿 床 ( 点) 分布情况吻合，Mn 元素衬度反映区内 Fe、Mn 矿化强度。绿汁江组 ( Pt2 lz) 、鹅头厂组( Pt2 e) 、落雪组 ( Pt2 l) 、因民组 ( Pt2 y) 、美党组( Pt2 m) 、大龙口组 ( Pt2 d) 、富良棚组 ( Pt2 f) 、黑山头组( Pt2 hs)是区内 Cu、Mn 地球化学场高值层位和铁铜矿主要赋存层位; 然后推广到全区，并表 5 易门铜矿区证据权法预测变量赋值条件Table 5 Variable assignment of evidence weight model in Yimen copper ore district裂隙的交汇部位，从图 5、图 6 可见，S 5 3、S 7 1，S 7 2 地段存在多期次的构造活动 ( 同 时出现 F3、F4 高因子得分) ，与其处于米茂断裂 与老吾街断裂部位的地质情况相吻合，同时也反 映了老吾街断裂东段的隐伏矿床特征。S 5 1， S 5 2，S 7 1，S 7 2 除了寻找 Fe、Cu 矿 产外，应 注 意 寻 找 Au 矿 化 信 息 ( F4 高 因 子 得 分) 。变量名称及含义赋值的条件2 2 累加异常值，小于全区衬度平均值的置为 02 2 衬度异常值，小于全区衬 度平均值的置为 0大于 20% 保留原值，否则为 0 大于 10% 保留原值，否则为 0 大于 5% 保留原值，否则为 0 大于 10% 保留原值，否则为 0 大于 20% 保留原值，否则为 0 大于 20% 保留原值，否则为 0 大于 20% 保留原值，否则为 0 大于 20% 保留原值，否则为 0 大于 20% 保留原值，否则为 0 实际值实际值 实际值( Cu + Ag + Bi) 衬度Mn 元素衬度绿汁江组地层百分比鹅头厂组地层百分比 落雪组地层百分比 因民组地层百分比 美党组地层百分比 大龙口地层百分比 富良棚组地层百分比 黑山头组地层百分比 地层组合熵 断裂等密度 断裂交点数断裂条数结 论( 1) 在云南易门地区，不同时代地层中的微 量元素分布特征是地球化学找矿预测的基础，它 是分析成矿物质来源的依据之一。( 2) 成矿关系中综合指标的地球化学因子更 能表达其整体性，因此从地学信息的挖掘角度来 看，利用改进的主成分因子 泛克里格模型成功 地提取易门铜矿区的地球化学示矿指标因子作为 成矿预测变量是找矿的关键。6表 6易门铜矿区找矿预测区及分级编号( 3)在综合预测过程中，采用基于 GIS 的地Table 6Prospective area of prospecting and classificationin Yimen copper ore district球化学证据权重找矿预测模型法，并结合三维地球化学多元信息图像进行人机交互式的综合分析， 对寻找因地形因素而掩盖的低缓弱地球化学异常 具有重要意义。预 测目标区级别靶 区 编 号 S 2 1，S 4 2S 2 2，S 4 1，S 4 3，S 8 2S 1 1，S 8 1Cu 矿找矿靶区( 4)在面临着地表矿产基本调查清楚，找矿难度愈趋艰难的情况下，常规地质方法难以直接发现隐伏矿床，地质地球化学综合预测方案的优 势日显其效，同时地球化学信息的充分挖掘对有 效、快速、准确地圈定找矿靶区的作用日趋显著。 工作中得到了云南省有色地质局地质地球物 理化学勘 查院晏建国总工程师、张 兴 林 工 程 师、张豫华工程 师、陈 永东高级工程师的大力支持和 帮助，在此表示衷心感谢。 S 9 1S 7 1，S 7 2，S 5 2S 3 1，S 3 2，S 6 1，S 5 1，S 5 3Fe 矿找矿靶区利用该模型计算得到包括已知矿 ( 床) 点在内的后验概率等值线，本次预测工作把后验概率大于0. 3 的预测单元作为 I 级远景区，共有 4 个; 大于0. 2 的预测单元作为级远景区，共有 7 个，大于0. 05 的预测单元作为远景区，共有 9 个。如图 5、图 6 所示，S 2、S 4、S 8 为 级 Cu 矿找矿远景区，S 1 为级 Cu 矿找矿远景 区，S 9 为 级 Fe 矿找矿远景区，S 5、S 7参考文献:1 Daskalaki S，Kopanas I，Goudara M，et al Data mining fordecision support on customer insolvency in telecommunications business J European Journal of Operational Research，2003，145 ( 2) : 239 255为级 Fe( Cu、Au)矿找矿远景区，S 3、S 6矿化找矿远景区。 2 王燕，李睿，李明 数据挖掘技术应用研究 J 甘肃为级受裂隙控制的 Fe( Mn)科技，2001，17 ( 1) :49 50在此基础上，结合地质情况，进一步圈定了 17 个找矿靶区 ( 表 6) 。由于 F3、F4 因子得分图反映了区内构造裂隙 发育程度，其因子得分较高的地段也是多组构造3 龚建华，林珲 虚拟地理环境 M 北京: 高等教育出版社，2002:100 1014 Ananthanarayana V S，Narasimha Murty M，Subramanian DK Tree structure for efficient data mining using rough setsJ Pattern Recognition Letters，2003，24( 6) : 851 8625 阮天健，朱有光 地球化学找矿 M 北京: 地质出版社，19856 方维萱 试论地球化学场的研究方法 J 西北地质，1999，32 ( 1) : 28 3512715 Agterberg F P，Bonham-Carter G F，Wright D F Statistical pattern integration for mineral exploration C / /Gaal G， Merriam D F Computer Applications for Mineral Exploration in Resource Exploration， Oxford: Pergamon Press，1990:1 2116 陈元坤，吴上龙 云南省区域物化探资料综合研究报告R 昆明: 云南省地质矿产局地球物理地球化学勘查 队，19907 云南省地质矿产局 云南省区域地质志 M质出版社，1990: 572 610北京: 地8 龚琳， 何 毅 特， 陈 天 佑 云南东川元古宙裂谷型铜矿 M 北京: 冶金工业出版社，1996: 4 99 云南省地质局第二区域地质测量大队 1 20万昆明幅区域地质调查报告 R 昆明: 云南省地质矿产勘查开发局， 17 谢学锦，王学求 深穿透地球化学新进展 J 地学前197110 云南省地质局第二区域地质测量大队 1 20万玉溪幅区 域地质调查报告 R 昆明: 云南省地质矿产勘查开发 局，1969缘，2003，10 ( 1) :225 23818 陈三明，钱建平，陈宏毅 桂东南植被覆盖区的抗干扰遥感蚀变信息优化提取与找矿预测 J 桂林理工大学学报，2010，30 ( 1) :33 4011 张远志 云南省岩石地层 M版社，1996武汉: 中国地质大学出19 陈三明，吴虹，谭泛，等 基于 ASTER 的遥感地球化学统计预测模型及应用以金川铜镍矿床外围找矿为例12 中国有色金属工业总公司西南地质勘查局 昆阳裂谷铜矿地质及成矿预测云南省东川 元江区带科研找矿 总结报告 R 昆明 云南省有色地质局，199413 潘龙驹，付金德，孙恩守，等 满洲里 新右旗铜银多 金属矿带大型矿床地质特征 R 北京: 中国有色金属 工业总公司地质勘查总局14 赵鹏大，王京贵，饶明辉，等 中国地质异常 J 地 球科学中国地质大学学报，1995，20 ( 2) : 117 J桂林理工大学学报，2010，30 ( 4) :480 48920 刘洪滔，李家盛，晏建国 “三江” 中南段综合信息找矿模型 及 找 矿 方 向 J 中 国 地 质，2008，35 ( 1 ) :101 11021 刘星，胡光道 应用 MORPAS 系统证据权重法进行多源 信息成矿预测以澜沧江南段地区为例 J 地质与 勘探，2004，40 ( 4) : 66 69Prospecting Prediction by Geochemical Information Miningin Yimen Copper Area，YunnanLIU Hong-tao1 ，CHEN San-ming2，3 ，ZHU Jia-qun1( 1. Institute of Geological Geophysical and Geochemical Exploration，Yunnan Non-ferrous Metals Geological