02-云南个旧锡矿预测

1、一、个旧东区区域成矿地质背景一、个旧东区区域成矿地质背景 个个 旧旧 东东 区区 地地 质质 图图个旧锡矿区地质图个旧锡矿区地质图（据西南有色地勘局（据西南有色地勘局308308队）队） T2g -中三叠个旧组；中三叠个旧组； T2f -中三叠统法郎组；中三叠统法郎组； 1 -断裂；断裂； 2 -背斜轴；背斜轴； 3 -花岗岩；花岗岩； 4 -二长岩；二长岩； 5 -辉长岩；辉长岩； 6 -碱性岩；碱性岩； 7 -取样地段。取样地段。三、个旧东区锡多金属多元信息综合三、个旧东区锡多金属多元信息综合定量预测与评价定量预测与评价多元信息综合定量预测与评价流程多元信息综合定量预测与评价流程地、物、化

2、、遥、矿产资料地、物、化、遥、矿产资料地层、构造、岩浆地层、构造、岩浆岩信息及与矿床相岩信息及与矿床相关地质信息图关地质信息图主要成矿元素异常主要成矿元素异常与元素综合异常信与元素综合异常信息提取息提取重力、航磁异重力、航磁异常信息提取常信息提取线性、环形构线性、环形构造解译，蚀变造解译，蚀变信息提取信息提取成矿规律剖析成矿规律剖析建立地质概念模型建立地质概念模型建立证据权模型，选择证据因建立证据权模型，选择证据因子子成矿预测、条件独立性检验成矿预测、条件独立性检验后验概率图后验概率图矿区资源预测与评价矿区资源预测与评价研究区赋矿地层与矿体分布图研究区赋矿地层与矿体分布图 %1001lg1lg

3、NNPiPiH地层组合熵与矿体分布图地层组合熵与矿体分布图 岩浆岩区域分布岩浆岩区域分布GIS缓冲区分析与已知矿点分布图缓冲区分析与已知矿点分布图 ijjiipilAd11ijjiipilAd11 经验告诉我们，矿体不经验告诉我们，矿体不一定分布在重要断裂内部，一定分布在重要断裂内部，因为断裂内部往往是应力因为断裂内部往往是应力集中部位，不利于成矿物集中部位，不利于成矿物质的富集，成矿物质会运质的富集，成矿物质会运移到离断裂一定距离的低移到离断裂一定距离的低应力环境中富集成矿。鉴应力环境中富集成矿。鉴于此，对断裂做于此，对断裂做250m 缓缓冲区分析，从图中可以看冲区分析，从图中可以看出，出，

4、矿体的大部基本都落矿体的大部基本都落在断裂缓冲区内部在断裂缓冲区内部。 研究区构造线密度与矿体分布图研究区构造线密度与矿体分布图 断断 裂裂 缓缓 冲冲 区区 与与 矿矿 体体 分分 布布 图图 物物 探探 异异 常常 与与 矿矿 体体 分分 布布 图图 研究区区域地球化学背景值表研究区区域地球化学背景值表 单位单位 ppm 本区的主要成矿元素本区的主要成矿元素Cu、Pb、Zn、Ag、Sn、W、Mo变异系数均大于变异系数均大于1，表明表明这些元素在本区的分布极不均匀这些元素在本区的分布极不均匀。一般认为。一般认为,正常的沉积岩中元素含量的正常的沉积岩中元素含量的分布应该是相对均匀的，其原始分布

5、量变化非常缓慢，不易察觉，如果分布分布应该是相对均匀的，其原始分布量变化非常缓慢，不易察觉，如果分布极不均匀，则极不均匀，则可能是某种地质作用使其在局部地段发生了富集或分散可能是某种地质作用使其在局部地段发生了富集或分散 。处理特异值之后处理特异值之后Ag的分布型式图的分布型式图 对数变换之后对数变换之后Ag的分布型式图的分布型式图 处理特异值（左）和取其常用对数（右）后处理特异值（左）和取其常用对数（右）后Cu的分布型式图的分布型式图 处理特异值（左）和取其常用对数（右）后处理特异值（左）和取其常用对数（右）后Pb的分布型式图的分布型式图 处理特异值（左）和取其常用对数（右）后处理特异值（左

6、）和取其常用对数（右）后Zn的分布型式图的分布型式图 处理特异值（左）和取其常用对数（右）后处理特异值（左）和取其常用对数（右）后W的分布型式图的分布型式图 处理特异值（左）和取其常用对数（右）后处理特异值（左）和取其常用对数（右）后Mo的分布型式图的分布型式图 处理特异值（左）和取其常用对数（右）后处理特异值（左）和取其常用对数（右）后Sn 的分布型式图的分布型式图 xziixxxziimaxmin异异 常常 区区 与与 矿矿 体体 分分 布布 图图 Sn异异 常常 区区 与与 矿矿 体体 分分 布布 图图 Ag异异 常常 区区 与与 矿矿 体体 分分 布布 图图 Cu异异 常常 区区 与与

7、 矿矿 体体 分分 布布 图图 Pb异异 常常 区区 与与 矿矿 体体 分分 布布 图图 Mo异异 常常 区区 与与 矿矿 体体 分分 布布 图图 W异异 常常 区区 与与 矿矿 体体 分分 布布 图图 ZnCu、Ag、Cu、Pb、Zn、W 、Mo 七元素聚类分析谱系图七元素聚类分析谱系图 其中，其中，p1=0.86，p2=0.84，p3=0.81，p4=0.8，p5=0.78，p6=0.6 为方便讨论，先作以下假设。假设各元素水系沉积物化探取样值分别表示为：为方便讨论，先作以下假设。假设各元素水系沉积物化探取样值分别表示为：铜铜DCu，铅，铅DPb，锌，锌DZn，银，银DAg，锡，锡DSn，

8、钨，钨= DW，钼，钼= DMo 化探综合异常与矿体分布图化探综合异常与矿体分布图 研究区研究区ETM742 遥感影像图遥感影像图 研究区三维遥感影像图研究区三维遥感影像图 1-5、7波段主成分分析波段主成分分析PC3 33低通滤波低通滤波7、PC1 、滤波结果假彩色合成、滤波结果假彩色合成PC2 55高通滤波高通滤波7、PC1、滤波结果假彩色合成、滤波结果假彩色合成主干线性构造解译主干线性构造解译次级线性构造解译次级线性构造解译遥感解译线性影像图遥感解译线性影像图线性构造遥感解译流程图线性构造遥感解译流程图 主干线性构造遥感解译图主干线性构造遥感解译图 次级线性构造遥感解译图次级线性构造遥感

9、解译图 环形构造遥感解译图环形构造遥感解译图 铁染蚀变与矿体叠加图铁染蚀变与矿体叠加图 羟基蚀变与矿体叠加图羟基蚀变与矿体叠加图 预测结果及评价预测结果及评价 以所建立的个旧东区证据权以所建立的个旧东区证据权模型，计算各个预测单元的成模型，计算各个预测单元的成矿有利度（以成矿的后验概率矿有利度（以成矿的后验概率值来表示），按照后验概率相值来表示），按照后验概率相对大小，分为不同等级，不同对大小，分为不同等级，不同等级赋予不同颜色。等级赋予不同颜色。 由图中可以看出，图中深蓝由图中可以看出，图中深蓝色（高后验概率）图斑呈现出色（高后验概率）图斑呈现出良好的分带性，故根据后验概良好的分带性，故根据

10、后验概率的大小与分布，划分出四个率的大小与分布，划分出四个成矿区和四个成矿预测区。成矿区和四个成矿预测区。 后验概率分布图后验概率分布图 成成 矿矿 预预 测测 图图 成矿区：成矿区：I、 小马拉格成矿区小马拉格成矿区II、 松树脚成矿区松树脚成矿区III、老厂成矿区、老厂成矿区IV、 卡房成矿区卡房成矿区成矿预测区：成矿预测区：A、乌谷哨成矿预测区、乌谷哨成矿预测区B、半坡成矿预测区、半坡成矿预测区C、大屯镇成矿预测区、大屯镇成矿预测区D、羊坝底成矿预测区、羊坝底成矿预测区E 、马鹿唐预测区、马鹿唐预测区F 、卡田预测区、卡田预测区FE 个个 旧旧 东东 区区 地地 质质 图图13-2凹陷勘

11、查区凹陷勘查区 13-2凹陷勘查区位于个旧老卡花岗岩体东南缘，凹陷勘查区位于个旧老卡花岗岩体东南缘，因此受岩体作用明显，构造复杂，多深大断裂，成因此受岩体作用明显，构造复杂，多深大断裂，成矿作用显著。矿作用显著。（1）、地层岩性）、地层岩性 13-2勘查区位于个旧岩体东南缘，受岩体作用显著。地层变形勘查区位于个旧岩体东南缘，受岩体作用显著。地层变形强烈。强烈。地层组成主要是三叠系中统个旧组卡房段（地层组成主要是三叠系中统个旧组卡房段（T2g1）碳酸盐）碳酸盐岩岩，为夕卡岩型锡多金属矿床形成提供了有利的围岩条件，矿区，为夕卡岩型锡多金属矿床形成提供了有利的围岩条件，矿区内矿床、矿体的富集产出受到

12、其中某些层位和岩性组合的控制。内矿床、矿体的富集产出受到其中某些层位和岩性组合的控制。卡房段地层的主要特征是灰岩、灰岩与白云质灰岩或含泥质灰岩卡房段地层的主要特征是灰岩、灰岩与白云质灰岩或含泥质灰岩互层，并有玄武岩、辉绿岩床分布，厚度互层，并有玄武岩、辉绿岩床分布，厚度500570m。卡房段容矿。卡房段容矿最为丰富，赋存的矿体占个旧东区探明矿体的最为丰富，赋存的矿体占个旧东区探明矿体的82.5%，锡储量为，锡储量为东区探明量的东区探明量的84.7%，铜储量为探明储量的，铜储量为探明储量的95.1%，铅储量占，铅储量占65.1%，银储量占，银储量占90%。 13-2凹陷勘查区断裂构造格架图凹陷勘

13、查区断裂构造格架图 乔巴硐断裂黑明槽断裂方家山断裂老熊硐断裂北西南东工作区范围形态 工作区DTM实体模型 北东西南工作区地层实体模工作区地层实体模型型 - T2g11；- T2g12；- T2g13；- T2g14；- T2g15；-花岗岩；-玄武岩 北西南东13-2凹陷勘查区断裂构造格架图凹陷勘查区断裂构造格架图 乔巴硐断裂黑明槽断裂方家山断裂老熊硐断裂北西南东13-2凹陷勘查区断裂建模统计表凹陷勘查区断裂建模统计表 工作区花岗岩实体模型工作区花岗岩实体模型 工作区岩浆岩实体模型工作区岩浆岩实体模型 花岗岩玄武岩东东南西北北南西工作区矿体模型总图工作区矿体模型总图 北东工作区氧化矿与花岗岩实

14、体模型关系图工作区氧化矿与花岗岩实体模型关系图 北东工作区硫化矿与花岗岩建模实体关系图工作区硫化矿与花岗岩建模实体关系图 北东工作区夕卡岩与花岗岩建模实体关系图工作区夕卡岩与花岗岩建模实体关系图 北东工作区钻孔分布图工作区钻孔分布图 北东工作区三维网格模型表面形态工作区三维网格模型表面形态 断裂三维网格模型断裂三维网格模型 北东夕卡岩硫化矿氧化矿矿体三维网格模型矿体三维网格模型 北东不同岩性地层中含氧化矿单元体数目比较图不同岩性地层中含氧化矿单元体数目比较图 不同岩性地层中含硫化矿不同岩性地层中含硫化矿单元体数目比较图单元体数目比较图 不同岩性地层中含夕卡岩不同岩性地层中含夕卡岩单元体数目比较

15、图单元体数目比较图 不同地层的氧化矿含矿性统计不同地层的氧化矿含矿性统计 不同地层的夕卡岩含矿性统计不同地层的夕卡岩含矿性统计 不同地层的硫化矿含矿性统计不同地层的硫化矿含矿性统计 2.0 5.0研究区研究区Sn的三维异常图的三维异常图 北东已知矿体已知矿体Sn三维异常三维异常北西南东研究区研究区Pb的三维异常图的三维异常图 1 20北东已知矿体已知矿体Pb三维异常三维异常东东北西南1 2.44研究区研究区Cu的三维异常图的三维异常图 北西已知矿体已知矿体Cu三维异常三维异常北西南东“立方体预测模型立方体预测模型”找矿方法体系结构找矿方法体系结构 “立方体预测模型立方体预测模型”预测流程预测流

16、程 北东断裂及其缓冲区关系图断裂及其缓冲区关系图 北东断裂缓冲区单元格模型断裂缓冲区单元格模型 北东花岗岩体及其缓冲区实体模型花岗岩体及其缓冲区实体模型 花岗岩实体模型花岗岩实体模型花岗岩缓冲区实体模型花岗岩缓冲区实体模型北东玄武岩及其缓冲区实体模型玄武岩及其缓冲区实体模型 北东立方体预测变量统计表立方体预测变量统计表 )()/(lg)(APBAPIBASSNNIjjBAlg)(夕卡岩找矿信息量计算结果表夕卡岩找矿信息量计算结果表 硫化矿找矿信息量计算结果表硫化矿找矿信息量计算结果表 氧化矿找矿信息量计算结果表氧化矿找矿信息量计算结果表 信息量大于信息量大于0.8的夕卡岩矿块模型的夕卡岩矿块模

17、型 0.81.092北北东东0.8001.267信息量大于信息量大于0.8的硫化矿矿块模型的硫化矿矿块模型 北北东东0.8001.230信息量大于信息量大于0.8的氧化矿矿块模型的氧化矿矿块模型 北北东东Sn工程验证区工程验证区2 5 北东工程验证区工程验证区A中心坐标：中心坐标：X：574260 Y：124240 Z：2130 找矿概率：找矿概率：45.63工程验证区工程验证区B中心坐标：中心坐标：X：547200 Y：124310 Z：1800 找矿概率：找矿概率：48.25工程验证区工程验证区AB的预测的预测资源量总和：资源量总和：114324.44 吨。吨。 北东已知矿体已知矿体Sn三

18、维异常三维异常北东1 20工程验证区中心坐标：工程验证区中心坐标： X：574230 Y：124260 Z：2000 找矿概率：找矿概率：48.25预测资源量：预测资源量：215560.09吨吨 Pb工程验证区工程验证区北东已知矿体已知矿体Pb三维异常三维异常西北1 2.44Cu工程验证区工程验证区工程验证区坐标：工程验证区坐标： X：574120 Y：124150 Z：1820 找矿概率：找矿概率：39.39预测资源量：预测资源量：65060.69 吨吨北西已知矿体已知矿体Cu三维异常三维异常信息量大于信息量大于1.0的夕卡岩的夕卡岩立方体模型立方体模型断裂实体模型断裂实体模型地层地层T2g13-2实体模型实体模型研究区夕卡岩重点预测区研究区夕卡岩重点预测区 东北信息量大于信息量大于1.1的硫化矿的硫化矿立方体模型立方体模型断裂实体模型断裂实体模型地层地层T2g11实体模型实体模型研究区硫化矿重点预测区研究区硫化矿重点预测区 北东信息量大于信息量大于1.1的氧化矿的氧化矿立方体模型立方体模型断裂实体模型断裂实体模型地层地层T2g13-2实体模型实体模型 研究区氧化矿重点预测区研究区氧化矿重点预测区 北东信息量大于信息量大于0.8的夕卡岩矿块模型与断裂缓冲区及已知夕卡岩实体模型关系图的夕卡岩矿块模型与断裂缓冲区及已知夕卡岩实体模型关系图 信息量大于信息量大于0.8的夕卡的夕