云南个旧矿区某隐伏矿床大比例尺三维预测

　 2009年 1月 地 　质 　科 　学 CH INESE JOURNAL OF GEOLOGY 　 44 (1) : 324—337 云南个旧矿区某隐伏矿床大比例尺三维预测3 陈建平　尚北川　吕　鹏　赵　洁　胡　青 (中国地质大学国土资源信息开发研究重点实验室 北京　100083) 摘 　要 　使用钻孔数据提取云南个旧某矿区的元素异常分布 ,应用三维地质建模软件对之进 行地下隐伏矿体的三维模拟预测。将三维可视化技术与找矿信息量法结合提出“立方体预测 模型 ”找矿 ,选取 5个统计分析变量 (地层、断裂、花岗岩、Sn异常及 Pb异常 )的地质标志 进行了找矿信息量计算 ,进而估算其金属资源储量。确定地层和断裂是研究区主要的控矿因 素 ,以断裂周围 150 m区域 ,个旧组卡房段的第五和第六亚段地层为远景预测区 ,预测氧化矿 体远景区的单元数约占总单元数的四分之一。这表明“立方体预测模型 ”法在该区隐伏矿体 定位、定量和定概率预测方面的有效性 ,实现了传统的二维找矿向三维找矿的新突破。 关键词 　三维地质建模 　立方体预测模型 　成矿预测 　隐伏 Pb2Sn矿体 　云南个旧 中图分类号 : P618 文献标识码 : A 文章编号 : 0563 - 5020 (2009) 01 - 324 - 14 　3 北京市教育委员会共建项目建划 (JD104910556)的部分成果。 陈建平 ,男 , 1959年 6月生 ,教授 ,经济地质学专业。e2mail: 3 s@ cugb1edu1cn 2007 - 11 - 26收稿 , 2008 - 03 - 28改回。 随着地表和浅部矿及易识别矿的日益减少 ,找矿难度日益增大 ,找矿效果日益降低。 进入 20世纪 80年代以来 ,重点寻找深部隐伏矿、新类型矿及新领域矿已成为世界各国所 关注的找矿方向 (胡慧民等 , 1995)。隐伏矿床的大比例尺三维预测成为成矿预测领域的 前沿课。本研究是应用“立方体预测模型 ”找矿方法并利用 M icrom ine公司的 KANTAN 3D软件 (M icrom ine国际矿业软件有限公司 , 1999～2004) ,对云南个旧某矿区进行的地 下隐伏矿体三维建模及预测 ,突破二维找矿的传统 ,实现了对该区矿体的定位、定量和定 概率的预测 (赵鹏大 , 2003; 邓明国 , 2005)。 1　“立方体预测模型”找矿方法 找矿信息量法是区域矿产预测中一种简便有效的方法。它是以地质异常理论为指 导 ,以地质、物探、化探和遥感等找矿信息为基础 ,通过计算各种地质因素和找矿标志所 提供的找矿信息量 ,定量地评价各地质因素和标志对指示找矿作用大小 ,确定有利成矿 部位 ,预测隐伏矿体 (王勇等 , 2002; 温礼琴等 , 2004)。 为实现传统的二维找矿向三维找矿的转换 ,开展研究区深部矿体定位、定量和定概 率一体化的三维预测 ,本研究以找矿信息量法为基础结合三维可视化技术提出了“立方 体预测模型 ”找矿方法。 “立方体预测模型”找矿方法首先通过研究矿区控矿地质条件和找矿标志在空间上、 特别是在深部的变化规律 ,综合分析处理各种深部找矿评价的定量化信息异常的三维实 体模型 ,再根据实体模型进行研究区三维立方体提取 ,并将找矿定量化信息赋予每一个 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 1期 陈建平等 :云南个旧矿区某隐伏矿床大比例尺三维预测 立方体 ,最后使用地质统计学等预测理论实现研究区深部矿体的三维预测 (侯景儒等 , 1998)。其基本如下。 1)对已有地质研究资料和矿区生产资料进行收集与整理 ,对矿区的地质、构造、物 探、化探和钻孔等数据建立矢量化数据库 ,特别要收集矿区地质剖面图、中段平面图和钻 孔编录资料。2)使用合适三维建模软件 ,利用收集到的矿区地质剖面图、中段地质图和钻 孔编录数据 ,建立研究区地质体三维实体模型。3)建立研究区的找矿模型 ,列出相关可利 用的找矿标志作为矿床三维预测变量并对其进行立方体预测单元的提取。4)将找矿定量 化信息赋予每一个立方体 ,对立方体单元所包含的数据进行统计处理 ,应用统计方法进 行找矿信息运算 ,进而确定找矿靶区、重点工作区和成矿有利区。 2　研究区地质背景及三维地质实体模型 (1)研究区地质背景及研究现状概述 个旧超大型锡 2铜 2多金属矿区位于云南省东南部 ,其中心地理位置 ,东经 103°09′ 26″,北纬 23°22′40″。它处于欧亚板块与太平洋板块、印度洋板块俯冲碰撞相接的部位 , 属于中国东南微板块西南缘右江陆缘盆地中的南盘江断褶带之西南隅。个旧矿区具有 优越的成矿地质背景及控矿条件 ,因而在地质矿床勘查中已取得了显著成效。目前个旧 已探明 Sn、Cu、Pb、Zn、W、B i、Mo、Au和 Ag等有色金属及贵金属矿产 20余种 ,矿体 10300 多个 ,累计探明有色金属储量逾 8 ×106 t,其中 Sn的探明储量约占世界 Sn储量的 24166% ; 占全国 Sn储量的 39114%。南北走向的个旧断裂纵贯个旧矿区 ,将其分为东区 和西区两部分。近东西向断裂控制了东区的马格拉、松树脚、高松、老厂和卡房 5大矿田 (图 1)。 个旧矿区锡 2多金属矿床的成因 ,长期以来以“花岗岩成因论 ”占主导地位 ,并作为矿 区找矿勘探的主要指导思想 (冶金工业部西南冶金地质勘探公司 , 1984; 杨世瑜 , 1990; 杨世瑜等 , 1994)。随着矿区开采深部硫化物型矿床 ,尤其是陆续发现新的层状型矿体 , 金祖德 (1991)明确指出 :个旧地区存在非热液成因的层间赤铁矿型锡矿 (即层间氧化 矿 ) ,是典型的三叠纪同生沉积矿床 ; 彭张翔 (1992)认为个旧锡矿为一成因比较复杂的 层控型锡 2多金属矿床 ,已有的花岗岩成矿模式不符合个旧矿区的实际情况 ; 谭允谦①认 为 : 个旧矿区高松矿田大箐东、马吃水和芦塘坝地区的层间氧化矿床是燕山晚期岩浆岩 期后含矿热液充填至古岩溶内形成的。秦德先等 (2006)基于一系列同位素年龄值提出 3 个系列的成矿作用 ,即 :印支期海底基性火山 —沉积成矿、海底喷流 —沉积成矿和燕山期 的花岗岩叠加改造成矿。 本次工作的研究区属高松矿田 ,位于大菁 —阿西寨向斜北翼 ,区内地表无岩浆岩体 露头 ,所出露地层为中三叠统个旧组碳酸盐岩 ,断裂构造十分发育 ,矿体的定位及产出形 态受地层、岩性及地质构造的控制十分明显 (图 2,图 3)。该区找矿应从研究向斜和断裂 构造入手 ,不仅注重锡矿化成矿 ,而且还应利用矿区金属元素分带规律进行多种矿种的 综合找矿。已知此区锡矿主要分布在芦塘坝断裂以西 ,矿体平均 Sn品位一般都可达到 523 ①谭允谦. 1993. 个旧层间氧化矿成矿空间特征. 个旧地质 , (1) : 9—21. © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 地 　质 　科 　学 2009年 图 1　个旧东矿区地质图 (据西南有色地勘局 308队 , 1984年未刊资料修改 ) 1. 全新统 ; 2. 上三叠统火把冲组砂砾岩和页岩 ; 3. 上三叠统鸟格组砂页岩 ; 4. 中三叠统法郎组砂页岩、灰岩夹玄 武岩 ; 5. 中三叠统个旧组白泥洞段灰岩、白云岩夹玄武岩 ; 6. 中三叠统个旧组马拉格段白云岩、钙质白云岩夹灰 岩 ; 7. 中三叠统个旧组卡房段灰岩、白云岩 ; 8. 下三叠统永宁镇组砂页岩夹灰岩 ; 9. 大箐—阿西寨向斜 ; 10. (中 2)晚燕山期花岗岩 ; 11. 印支期玄武岩 ; 12. 地层界线 ; 13. 断裂及产状 ; 断裂名称 : F1. 白沙冲断裂 ; F2. 岬界 山断裂 ; F3. 芦塘坝断裂 ; F4. 个松断裂 ; F5. 背阴山断裂 ; F6. 黄泥洞断裂 ; F7. 老熊洞断裂 ; F8. 仙人洞断裂 ; 　　　　　　　 F9. 白龙断裂 ; F10. 个旧断裂 ; 为研究区范围 Fig. 1　Geological map of the eastern m ining area of Gejiu 623 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 1期 陈建平等 :云南个旧矿区某隐伏矿床大比例尺三维预测 图 2　研究区断裂及勘探线分布图 F1. 芦塘坝断裂 ; F2. 麒麟山断裂 ; F3. 大箐东断裂 ; F4. 1号断裂 ; F5. 608断裂 ; F6. 黑蚂石断裂 ; F7. 坝西断裂 ; 　　　　　　　F8. 131断裂 ; A1Sn的 A验证区 ; B1Sn的 B工程验证区 ; C1Pb的 C工程验证区 ; D1Pb的 D工程验证区。 Fig. 2　D istribution of the fractures and p rospecting lines in the studied area 图 3　研究区地层实体模型 T2 g125. 中三叠统个旧组卡房段第五亚段灰岩、钙质白云岩 ,具最好的含矿性 ; T2 g126. 中三叠统个旧组卡房段第六亚 段钙质白云岩、白云质灰岩 ,具良好的含矿性 ; T2 g221. 中三叠统个旧组马拉格段第一亚段钙质白云岩 ; T2 g222. 中三 叠统个旧组马拉格段第二亚段钙质白云岩 ; T2 g223. 中三叠统个旧组马拉格段第三亚段钙质白云岩 (庄永秋等 , 1996) Fig. 3 　A model for ore2bearing strata of the studied area (1～2) %以上 ,伴生有 Ag、In、Fe和 A s,少数矿体伴生 Pb和 Zn, 8个矿体的 Sn储量达 50000 t以上 ; Ag、Pb和 Sn矿主要分布于芦塘坝断裂以东。 723 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 地 　质 　科 　学 2009年 (2)地层 该研究区内分布的地层均属中三叠统个旧组 (图 3) ,只发育下段卡房段 ( T2 g1 )和中 段马拉格段 ( T2 g2 ) ,缺失上段白泥洞段 ( T2 g3 )。 卡房段以中 2厚层灰岩为主夹薄层灰岩、白云岩。岩性变化频繁 ,互层带厚达数百米 , 层理发育 ,具去膏化及富藻层 ,有基性火山岩及凝灰岩。多含层间缓倾斜矿及接触带矿。 卡房段是个旧矿区最主要的容矿层位 ,矿区 19个容矿层中有 13个容矿层属卡房段 ,占矿 区 Sn储量的 90% , Cu储量的 96% , Pb储量的 44%。该段中的第五亚段和第六亚段 ( T2 g125和 T2 g126 )是本区的主要含矿层位。 马拉格段主要为中至厚层白云岩、钙质白云岩夹灰岩 ,互层带较卡房段发育差 ,互层 带厚度小 ,是矿区内第二个重要容矿层位。以发育层间陡倾斜管状矿为特征 ,也有缓倾 斜层间矿 ,接触带矿不发育。此段可分为 3个亚段 ,从下至上分别称第一亚段 ( T2 g221 )、 第二亚段 ( T2 g222 )和第三亚段 ( T2 g223 ) ,它们对本研究区矿体产出的贡献较少。 (3)构造 该研究区处于北西西向的大菁 —阿西寨向斜北翼 ,小的挠曲较明显 ,但对矿体的产 出起控制作用者以断裂构造 (图 4)及其派生裂隙为突出。 东西走向断裂主要有两条 ,均有较强矿化 ,并产出含锡氧化矿体。一是大菁东断裂 的近东西走向段 ,横贯整个研究区 ; 二是麒麟山断裂 ,为芦塘坝矿段南界 ,南侧有数条平 行破碎带产出。 北东走向断裂以芦塘坝断裂为主干构造 ,但本区更多体现的为其派生的一系列断 裂 ,它们通过个旧组卡房段第六亚段 ( T2 g126 )地层时产生频繁的层间剥离和破碎 ,形成了 多层次的层间脉与层间氧化矿床相交。南北向断裂主要为隐伏的小断裂 ,多见于芦塘坝 断裂的下盘 ,一般未见工业矿化。 (4)岩体 地表无岩浆岩出露 ,但深部有隐伏花岗岩体分布。据钻探工程揭露 ,该区隐伏花岗 岩体标高在海拔 900～10200 m ±,距地表约 10200～10400 m。花岗岩由北西向南东逐渐 倾伏 ,形成凹槽 ,在凹陷之中有数个岩脉或岩株小突起。花岗岩属燕山 (中 2)晚期 , 为壳 源重熔型黑云母斑状花岗岩 (赖大信 , 2005)。本区花岗岩对矿床形成的影响较小。 (5)已知矿体实体模型 研究区内与花岗岩密切相关的硫化矿及夕卡岩型矿体较少 ,以氧化矿为主。氧化矿 体属锡石 2硫化物矿床经氧化后形成的层间氧化矿床 ,按产状及形态可分为两种类型 : 1) 陡倾斜脉状矿体是成矿热液沿芦塘坝断裂及其派生的断裂破碎带充填交代而形成的 ,矿 体大小不一 ,一般含 Sn、Pb和 Ag品位较高 ; 2)缓倾斜似层状矿体主要分布在北东向芦 塘坝断裂旁侧的卡房段第六亚段 ( T2 g126 )地层中 ,受主断裂旁侧的次级层间剥离带或层 间破碎带控制 ,产状与地层一致 ,具有多层性。 (6)研究区钻孔模型 钻孔模型 (图 5)由原始数据库的 3个文件组成 ,分别为钻孔孔口三维坐标文件、钻孔 侧斜文件和化验数据文件。该模型共收集和使用钻孔 164个 ,化验数据 90610 个 ,为分 析该区元素异常三维分布提供基础。 823 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 1期 陈建平等 :云南个旧矿区某隐伏矿床大比例尺三维预测 3　研究区地质体三维立方体提取 根据现有地质资料对矿体的揭示 ,特别是勘探线的分布 ,结合矿体的形态、走向、倾 向和空间分布特征确定了建模的范围和基本参数。建模范围由南端的 206 /902103号钻 孔至北段的 205 /8924号钻孔 ,从东端的 214 /9722号钻孔至西端的 203 /91225号钻孔。模 型区形态是一个南北长 10908 m,东西长 10599 m,垂高为 10909 m 的立方体 ,总体积 3181 ×109 m3 ,单元块行 ×列 ×层为 20 m ×20 m ×20 m,次分块的规格为 10 m ×10 m × 10 m ,模型包括次分块总共有 308130678个单元块 (陈爱兵等 , 2004)。 在建立三维网格模型后 ,可以将找矿数字模型所确定的预测参数作为属性赋给每一 个单元块。如用地层实体模型对三维网格进行限定 ,划分出不同地层所包含的网格单 元 ,作为矿床预测中的岩性变量 ; 使用断裂实体模型对三维网格进行限定 ,划分出不同断 裂所包含的网格单元 ,作为矿床预测中的断裂变量。使用已知矿体实体模型对三维网格 进行限定 ,划分出不同矿体所包含的网格单元 (图 6) ,作为矿床预测中的先验条件。 同时可利用钻孔数据来分析这些单元块的元素三维异常分布。本研究选用距离反 比加权的方法对未知单元块进行了 Sn和 Pb含量的插值分析 (勒国栋 , 2003; 邓明国 , 2005)。最终确定 Sn异常下限取 413% , Pb的异常下限取 315% ,为工程验证区的确立提 供依据。图 7示出了研究区 Sn及 Pb异常与已知矿体的相对位置。 4　研究区三维成矿预测结果 (1)研究区找矿标志单元统计 研究区有效立方体总数为 303270577。在仔细研究主要控矿因素和有利找矿标志的 基础上 ,选取 5个统计分析变量的地质标志 (地层、断裂、花岗岩、Sn异常及 Pb异常 ) ,分 别对这些变量进行立方体单元的划分与提取如下 (刘春学等 , 2003)。 1)地层 　研究区地层有个旧组卡房段第五亚段 ( T2 g125 )和第六亚段 ( T2 g126 ) ,个旧组 马拉格段第一亚段 ( T2 g221 )、第二亚段 ( T2 g222 )和第三亚段 ( T2 g223 ) , 5个单元的地层分别 作为 5个变量 ,实际选取 T2 g125、T2 g126和 T2 g221。 2)断裂 　实际选取断裂面两侧 150 m的区域作为变量并进行立方体提取。 3)花岗岩 　以岩体顶面及四周向外延伸 150 m的区域作为变量进行立方体提取。 4) Sn异常和 Pb异常 　钻孔模型中已经对样品的 Sn和 Pb的含量进行了分析 ,得到 每个立方体单元的元素含量 ,将 Sn异常和 Pb异常作为两个预测变量。 具体预测所用变量及其统计结果见表 1。 (2)找矿信息量计算 以研究区找矿标志单元的统计结果为基础 ,使用找矿信息量法进行找矿有利度的计 算 ,具体步骤如下。 1)某找矿标志的找矿信息量 ( I)可用条件概率计算 ,但在实际应用中由于概率估计 上的困难 ,故采用频率值来估计概率值 (赵鹏大 , 2003) ,即 : IA (B) = lg N j N S j S (1) 923 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 地 　质 　科 　学 2009年 表 1　立方体预测变量统计表 Table 1　Statistical table of the cubic p redicting variates 找矿标志 标志所占立方体数 标志内氧化矿 立方体数 标志内夕卡岩 立方体数 标志内硫化矿 立方体数 研究区有效 立方体总数 T2 g 223 335042 0 0 0 3327577 T2 g222 310070 0 0 0 3327577 T2 g221 475850 18 0 0 3327577 T2 g 126 515465 4484 0 59 3327577 T2 g125 1691150 1710 250 384 3327577 断裂 2344020 6212 238 443 3327577 花岗岩 475692 0 250 343 3327577 Sn异常 788806 5741 17 90 3327577 Pb异常 685919 4530 17 84 3327577 矿块总数 — 6212 250 443 — 式中 , IA (B)为 A标志有 B矿的信息量 ; N j 为研究区内具有标志 A的含矿单元数 ; N 为研 究区内的含矿单元总数 ; S j 为研究区内具有标志 A的单元数 ; S为研究区的单元总数。 2)信息量 ( I)为正值的标志状态 ,虽都可提供找矿信息 ,但并非均为有利找矿标志 , 故利用式 (2)确定找矿有利标志的临界值①: ΔI+ = K∑ N j =1 Ij (2) 式中 ,ΔI+ 为有利找矿标志的临界值 , Ij 为选取有利找矿标志的界限 ,一般取有用信息水 平 K = 0175, ∑Nj =1 Ij 为正信息量的总和 ( j为正信息量的序号 , N 为正信息量的个数 )。由 于本研究区涉及变量数量较少 ,就不利用此式计算有用信息的累积临界值 ,而直接使用 所有正值信息量。此外 ,本区夕卡岩矿和硫化矿相对较少 ,且主要集中在花岗岩区域内 , 故略去 ,只针对氧化矿进行计算。计算各标志的氧化矿找矿信息量结果见表 2。 3)利用各找矿标志信息量分别计算每个基本单元的信息量总和 ( I) ,对于包含了找 矿信息量总和的预测单元 ,可以使用主观选择法、统计法和作图法确定找矿靶区包含信 息量的临界值 (据中国有色金属工业总公司北京矿产地质研究所① )。本文根据主观选择 法确定氧化矿找矿信息量临界值为 ≥1141,确定的含矿远景单元包含已知矿体 40393个 单元 ,占已知矿体的 70172%。 图 8为氧化矿综合信息量单元数统计 ,可得氧化矿找矿信息量 I≥0174的单元数有 9160138个 ,占总单元数的 27153% ; 氧化矿找矿信息量 I≥1129的单元数有 6620155个 , 占总单元数的 1919% ; 氧化矿找矿信息量 I≥1141 的单元数有 2880497个 ,占总单元数 的 8167% ; 氧化矿找矿信息量 I≥1196的单元数有 2290065个 ,占总单元数的 6188%。 (3)预测结果 1)工程验证区金属资源储量估算 033 ①中国有色金属工业总公司北京矿产地质研究所. 1987. 大比例尺矿床统计预测工作手册. 1—372. © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 1期 陈建平等 :云南个旧矿区某隐伏矿床大比例尺三维预测 图 4　研究区断裂构造三维格架图 F1. 麒麟山断裂 ; F2. 大菁东断裂 ; F3. 131断裂 ; F4. 608断裂 ; F5. 黑蚂石断裂 ; F6. 坝西断裂 ; F7. 104断裂 Fig. 4　A model for 3D fractures in the studied area 图 5　建模范围内钻孔与勘探线分布 Fig. 5　D istribution of drilling wells and p rospecting lines in the modelling area 根据对 Sn和 Pb含量进行的插值分析 ,确定 Sn异常下限取 413% , Pb的异常下限取 315% ,由此为工程验证区的确立提供依据 ,得到 A和 B两个 Sn的工程验证区 (图 9)及 133 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 地 　质 　科 　学 2009年 图 6　矿体三维网格模型 Fig. 6　3D gridding modle of the ore bodies 图 7　研究区 Sn、Pb的元素异常叠加图 (从北向南看 ) Fig. 7　Sn and Pb geochem ical anomalypis superposition of the studied area ( looking towards the south) C、D两个 Pb的工程验证区 (图 10) ,预测结果示于图 11。 233 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 1期 陈建平等 :云南个旧矿区某隐伏矿床大比例尺三维预测 图 8　氧化矿综合信息量单元数统计 Fig. 8　Statistics of the oxidized ore information variable 图 9　研究区 Sn的三维工程验证区 (A和 B区 ) Fig. 9　3D p roject verifying zone of Sn (A and B) 333 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 地 　质 　科 　学 2009年 图 10　研究区 Pb的三维工程验证区 (C区和 D区 ) Fig. 10　3D p roject verifying zone of Pb (C and D) 图 11　研究区重点预测区 Fig. 11　Key p redicting area of m ineralization in the studied area 矿体内部的不符合工业要求的岩石叫做夹石 ,杂石率是夹石在矿体中所占的百分 比。本研究中取杂石率为 20% (经验值 ) ,氧化矿验证区其比重平均值为 21451 t/m3。研 433 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 1期 陈建平等 :云南个旧矿区某隐伏矿床大比例尺三维预测 究区 Sn的工程验证区 A有立方体块 498个 , B区的立方体块数为 612个 ,得到 Sn的工程 验证区 (A区与 B区 )共 101100000 m3 ,预测金属 Sn的资源总量为 390569117 t; 研究区 Pb的工程验证区 C有立方体块 752个 , D区有立方体块 279个 , 得到 Pb的工程验证区 (C区与 D区 )共 100310000 m3 ,预测金属 Pb的资源总量为 1050099184 t。 表 2　氧化矿找矿信息量计算结果表 Table 2　Results of the information variable of oxidized ore in the studies area 找矿标志 具有标志 A的 含矿单元数 (N j ) 研究区内的 含矿单元总数 (N ) 具有标志 A的 单元数 (S j ) 研究区的 单元总数 S 标志 A存在已知 矿体的信息量 (N j /N ) / (S j /S ) 选取有利找矿 标志的界限 ( Ij ) T2 g221 18 6212 475850 3327577 01020262783 - 11693300902 T2 g 126 4484 6212 515465 3327577 11420300176 01668363088 T2 g125 1710 6212 1691150 3327577 01541639897 - 01266289352 断裂 6212 6212 2344020 3327577 11419602648 011521668 花岗岩 — 6212 475692 3327577 — — Sn异常 5741 6212 788806 3327577 31898647845 01590914008 Pb异常 4530 6212 685919 3327577 31537708424 01548722035 2)远景预测结果 根据信息量计算结果及预测结果 ,确定以找矿信息量 I≥1141的立方体作为研究区 重点预测区 (图 11) ,该区的位置受到地层 ( T2 g126 )和断裂的严格控制。 综上所述 ,地层和断裂是研究区主要的控矿因素 ,因此以断裂周围 150 m区域、T2 g125 和 T2 g126地层范围作为远景预测区 ,单元总数共计 301830905个。以信息量 ≥1141为综 合地质标志 ,通过计算得到研究区预测远景区预测氧化矿矿体单元总数为 70908个 ,而 研究区共有已知氧化矿矿体单元共计 60212个 ,因此还有 10696个氧化矿矿体单元未被 发现。 5　结束语 (1)立方体预测模型找矿方法从矿区尺度的定性找矿向定量找矿迈出了新步伐 ,实 现了传统的二维找矿向三维找矿的新突破 ,同时为当前寻找隐伏矿床的研究领域提供了 新思路 ,具有广阔的应用前景。 (2)本研究在建立了个旧矿区某勘察区三维矿床模型的基础上 ,使用钻孔数据提取 了该区的元素异常分布 ,据此确定了 Sn和 Pb的工程验证区。此外 ,应用“立方体预测模 型 ”找矿方法 ,计算了该研究区立方体单元块的综合信息量 ,并通过选取适当的信息量下 限确定了远景预测区。 (3)预测结果表明 :氧化矿体受到地层和断裂的严格控制 ,水平方向呈层状 ,垂向上 局限于海拔 10400 m 至 20000 m 之间。工程验证区预测金属 Sn 的资源总量为 390569117 t,金属 Pb的资源总量为 1050099184 t。确定的含矿远景单元中找到矿的评价 概率为 70172%。 533 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 地 　质 　科 　学 2009年 参 考 文 献 陈爱兵 ,秦德先 ,张学书 ,贾国相 ,谈树成 ,范柱国. 2004. 基于 M ICROM INE矿床三维立体模型的应用. 地质与勘探 , 40 (5) : 77—80. Chen A ibing, Q in Dexian, Zhang Xueshu, J ia Guoxiang, Tan Shucheng and Fan Zhuguo. 2004. 3D model for deposit based on M ICROM INE technology. Geology and Prospecting, 40 (5) : 77—80. 邓明国. 2005. 个旧矿区芦塘坝 10号矿群矿床模型研究 [硕士学位论文 ]. 昆明 : 昆明理工大学. 1—101. DengM ingguo. 2005. The Research of Lutongba No. 10 O rebody Group in Gejiu O re Field [Master Degree Thesis ]. Kun2 m ing: Kunm ing University of Science and Technology. 1—101. 侯景儒 ,尹镇南 ,李维明 ,向永生 ,黄竞先 ,胡平昭. 1998. 实用地质统计学 (空间信息统计学 ) . 北京 :地质出版社. 1— 200. Hou J ingru, Yin Zhennan, L iW eim ing, Xiang Yongsheng, Huang J ingxian and Hu Pingzhao. 1998. The App lied Geostatistics ( Statistics of Spatial Information) . Beijing: Geological Publ. House. 1—200. 胡慧民 ,赵鹏大 ,李紫金. 1995. 大比例尺成矿预测方法. 北京 :地质出版社. 1—170. Hu Huim in, Zhao Pengda and L i Zijin. 1995. The W ays of Large2Scale Metallogenetic Prognosis. Beijing: Geological Publ. House. 1—170. 金祖德. 1991. 个旧层状赤铁矿型锡矿热液成因之否定. 地质与勘探 , 27 (1) : 19—20. J in Zhude. 1991. Negation of hydrothermal origin of bedded hematite type tin ore deposit in Gejiu. Geology and Prospecting, 27 (1) : 19—20. 赖大信. 2005. 个旧锡矿大马芦层间氧化矿床地质特征及找矿方向研究. 矿产与地质 , 19 (1) : 25—28. Lai Daxin. 2005. Geological characteristics and ore p rospecting orientation of Damalu interlayer oxidized deposit in Gejiu tin m ine. M ineral R esources and Geology, 19 (1) : 25—28. 勒国栋 ,刘衍聪 ,牛文杰. 2003. 距离加权反比插值法和克拉金插值法的比较. 长春工业大学学报 , 24 (3) : 53—57. Lei Guodong, L iu Yancong and N iu W enjie. 2003. Comparison between Inverse D istance W eightingMethod and Kriging. Jour2 nal of Changchun U niversity of Technology, 24 (3) : 53—57. 刘春学 ,秦德先 ,党玉涛 ,谈树成. 2003. 个旧锡矿高松矿田综合信息矿产预测. 地球科学进展 , 18 (6) : 921—927. L iu Chunxue, Q in Dexian, Dang Yutao and Tan Shucheng. 2003. Synthesis information based m ineral resource p rediction of Gaosong field in Gejiu tin deposit. Advances in Earth Science, 18 (6) : 921—927. 彭张翔. 1992. 个旧锡矿成矿模式商榷. 云南地质 , 11 (4) : 362—368. Peng Zhangxiang. 1992. The m ineralizing mode of Gejiu tin deposit. Yunnan Geology, 11 (4) : 362—368. 秦德先 ,黎应书 ,谈树成 ,陈爱兵 ,薛传东 ,范柱国 ,党玉涛 ,童 　祥 ,武俊德 ,李玉新 ,王海云. 2006. 云南个旧锡矿的成 矿时代. 地质科学 , 41 (1) : 122—132. Q in Dexian, L i Yingshu, Tan Shucheng, Chen A ibing, Xue Chuandong, Fan Zhuguo, Dang Yutao, Tong Xiang, W u Junde, L i Yuxin and W ang Haiyun. 2006. Metallogenic ages of Gejiu tin ore deposit in Yunnan p rovince. Chinese Jour. Geol. , 41 (1) : 122—132. 王　勇 ,陈邦国. 2002. 找矿信息量法在铜山铜矿成矿预测中的应用. 地质与勘探 , 38 (3) : 49—51. W ang Yong and Chen Bangguo. 2002. The app lication of p rospecting2information contents in metallogenic p rognosis at Tongs2 han copper deposit. Geology and P rospecting, 38 (3) : 49—51. 温礼琴 ,游远航 ,杨　斌. 2004. 信息量计算法在武当西缘银多金属矿成矿预测中的应用. 资源环境与工程 , 18 ( 3) : 74—77. W en L iqin, You Yuanhang and Yang B in. 2004. The app lication of information content computation in metallogenic p rognosis in western W udang silver2polymetalic deposit. Resources Environm ent & Engineering, 18 (3) : 74—77. 杨世瑜. 1990. 滇东南锡矿带矿床类型及其组合特征. 矿床地质 , 9 (1) : 35—48. Yang Shiyu. 1990. The type and combinatorial characteristic of stannary in south2east of Yunnan. O re D eposit Geology, 9 (1) : 633 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net 1期 陈建平等 :云南个旧矿区某隐伏矿床大比例尺三维预测 35—48. 杨世瑜 ,颜以彬. 1994. 云南的锡矿床与花岗岩类在时空分布上的关系. 云南地质 , 13 (2) : 149—157. Yang Shiyu and Yan Yibin. 1994. The spatio2temporal relation between tin ore deposit and granite. Yunnan Geology, 13 (2) : 149—157. 冶金工业部西南冶金地质勘探公司. 1984. 个旧锡矿地质. 北京 :冶金工业出版社. 1—256. Southwest Metallurgical and Geological Prospect Corp. ,M inistry ofMetallurgical Industry. 1984. Geology of Gejiu Tin O re De2 posit. Bejing:Metallurgical Industry Press. 1—256. 赵鹏大. 2003. 定量地学方法及应用. 北京 :高等教育出版社. 1—459. Zhao Pengda. 2003. Quantitative Geoscience Methodes and Its App lications. Beijing: H igher Education Press. 1—459. 庄永秋 ,王任重 ,杨树培 ,尹金明. 1996. 云南个旧锡铜多金属矿床. 北京 :地震出版社. 1—145. Zhuang Yongqiu,W ang Renzhong, Yang Shupei and Yin J inm ing. 1996. Geology of Gejiu Tin2Copper Polymetallic Deposit. Beijing: Seismological Press. 1—145. Large2sca le 3D m eta llogen ic pred iction of concea led orebody in Gejiu , Y unnan prov ince Chen J ianp ing　Shang Beichuan　LüPeng　Zhao J ie　Hu Q ing (Laboratory of Land Resources and Inform ation D evelopm ent, China U niversity of Geosciences, B eijing　100083) Abstract Based on the drilling data of the m ining area in Gejiu, Yunnan p rovince, the 3D p redic2 tion of a conceal orebody was p rocessed whith the p roper 3D geological modeling software1One new method, “the cubic p redicting model”method, was brought forward by combining the 3D visualization technology and the traditional quantitative p redicting method. Five statistical vari2 ates, namely, the strata, faults, granites, anomalies of Sn and Pb, were selected to be geolog2 ical marks, and we established the p rospective zone and tested p roject zones of Sn and Pb, cal2 culated the ore2containing p robabilities and the quantities of the metal resource. The key con2 trolling factors in the studied area should be strata and faults, so we thought that the 5 th and 6 th Submembers, Kafang Member of M iddle Triassic Gejiu Formation, and the region within 150 m of the both sides of fault should be potential metallurgic regions. This app lication showed the validity of“the cubic p redicting model”method in orebodypis p rediction for determ2 ing quantity, position and p robability. Fourthermore, this study made a great breakthrough of changing from traditional 2D to the 3D p rospecting. Key W ords　3D geological modeling, Cubic p redicting model, Metallogenic p rediction, Concealed Pb2Sn orebody, Gejiu in Yunnan 733 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net