证据权模型与含矿网格单元蒙特卡洛法相结合的资源潜力-地质科学.pdfVIP

蒲路平 1327 地 质 科 学 2017年10月 CHINESE JOURNAL OF GEOLOGY 52（4）：1327—1339 DOI：10.12017/dzkx.2017.082 证据权模型与含矿网格单元蒙特卡洛法 相结合的资源潜力评价新方法* 浦路平 朱国器 林锦富 陈三明1 2 1 1 （1.桂林理工大学 桂林 541004；2.广西壮族自治区地球物理勘察院 广西柳州 545005） 摘 要 本文提出了计算含矿单元数概率分布函数的新方法，在此基础上，把证据权模型矿 产资源评价方法得到的含矿网格单元作为蒙特卡洛模型中的统计对象，得到了与证据权模型 找矿法结合的含矿网格单元蒙特卡洛法矿产资源潜力评价方法。该方法可以把评价工作区的 矿产空间分布信息、成矿模式和找矿模型、成矿—找矿空间信息引入蒙特卡洛方法中，无需 估计矿床个数分布，减少了工作步骤和评价误差，可提高工作效率和自动化。对实际案例的 应用表明该方法是可行的。 关键词 含矿网格单元 蒙特卡洛法 矿产资源潜力评价 证据权模型 丹池锡矿带 中图分类号：P628 文献标识码：A 文章编号：0563-5020（2017）04-1327-13 蒙特卡罗方法也称统计模拟法，是根据统计抽样理论，通过对随机变量函数的概 率模拟和统计试验来近似求解的方法（Ulametal.，1947；Metropolisetal.，1949； Ulam，1950；VonNeumann，1951；Metropolis，1987）。由于地质成矿及数据获取过程 也是随机过程，蒙特卡洛法也应用于矿产资源评价。 现行蒙特卡洛矿产资源潜力评价方法以矿床为统计对象建立统计模型，进行随机 抽样来计算矿产资源总量，需要组织人类专家采用主观概率法估计或者用密度估计方 法统计工作区潜在的矿床数分布函数，导致误差较大，步骤繁琐，难以利用矿产空间 信息和成矿—找矿模型及成矿—找矿信息。 本文在提出按照证据权后验奇比（或后验概率）阈值统计含矿单元频率得到的含矿 单元数概率分布函数计算方法的基础上，根据总和式与非总和式结合矿产资源潜力评 价方法的思想，用证据权模型矿产资源评价方法得到的含矿网格单元代替矿床作为蒙 特卡洛模型中的统计对象，得到了基于证据权模型的含矿网格单元蒙特卡洛矿产资源 潜力评价方法。这种新方法引入评价工作区的矿产空间分布、成矿模式和找矿模型以 及成矿—找矿空间分布信息到蒙特卡洛方法中，可以省略估计矿床个数的工作步骤， 可以提高矿产资源评价工作效率和工作自动化水平，而且能够减少预测误差。并应用 该方法对广西丹池锡多金属成矿带的热液型锡矿进行了矿产资源潜力评价，表明该方 * 国家自然科学基金项目（编号41372339）、广西自然科学基金项目（编号：2013GXNSFAA019276） 和桂林理工大学博士基金项目资助。 浦路平，男，1958年9月生，博士，副教授，矿产普查与勘探学专业。E⁃mail：lupingpu@ 2016-12-10收稿，2017-07-30改回。 蒲路平 1328 1328 地 质 科 学 2017年 法预测结果与已知矿床分布吻合性较好。 1 原理和方法 1.1 原理 矿产资源评价方法有总和式和非总和式两大类（赵鹏大，1983），分别评价预测工 作区矿产资源的总量和各个子空间含矿性。蒙特卡洛法（Drewetal. 1984）、分形法（申 维，1997）、齐波夫律法（RowlandsandSampey，1977）等总和式方法只能获得工作区矿 产资源总量预测值，没有空间定位功能。而回归方法（ChungandAgterberg，1980）、人 工神经网络 （Singer and Kouda，1996）、证据权法 （Bonham⁃Carter et al.，1989； Agterberg，1992）等非总和式方法虽然可以预测工作区空