矿岩崩落块度预测的节理面三维模拟技术及其应用.pdf

矿岩崩落块度预测的节理面三维模拟技术及其应用 ① 彭平安１，２， 王李管１，２ （１．中南大学 资源与安全工程学院，湖南 长沙 ４１００８３； ２．中南大学 数字矿山研究中心，湖南 长沙 ４１００８３） 摘 要 基于传统二维状态下块度预测存在的弊端，提出以构造面的空间展布规律为基础，利用 Ｍｏｎｔｅ Ｃａｒｌｏ 技术模拟节理系统，采 用三维实体切割技术建立块体三维截取模型的矿岩块度预测新模型，实现了矿岩崩落块度的三维预测，并开发了相应软件ＯＦＰ １．０ （Ｏｒｅ Ｆｒａｇｍｅｎｔｓ Ｐｒｅｄｉｃｔｉｏｎ）。 研究成果应用于某特大型斑岩铜矿，预测结果表明该矿山等效尺寸大于 １．６５ ｍ 的块体筛上累积百分 比为 ４０％，大于 １．２６ ｍ 的块体筛上累积百分比为 ５９．８５％。 根据等效尺寸与块体体积的转换关系，建议放矿口尺寸设置为 ５ ｍ ５ ｍ， 此时放矿口堵塞几率较小。 研究成果为自然崩落法矿山采矿方法设计和工程中的放矿控制提供了重要依据。 关键词 自然崩落法； 块度预测； Ｍｏｎｔｅ Ｃａｒｌｏ 法； 节理面模拟 中图分类号 ＴＤ８５３文献标识码 Ａｄｏｉ１０．３９６９／ ｊ．ｉｓｓｎ．０２５３－６０９９．２０１５．０３．００６ 文章编号 ０２５３－６０９９（２０１５）０３－００２２－０５ Ｔｅｃｈｎｉｑｕｅ ａｎｄ Ａｐｐｌｉｃａｔｉｏｎ ｏｆ Ｔｈｒｅｅ⁃ｄｉｍｅｎｓｉｏｎａｌ Ｓｉｍｕｌａｔｉｏｎ ｏｆ Ｊｏｉｎｔ Ｐｌａｎｅ ｆｏｒ Ｏｒｅ Ｆｒａｇｍｅｎｔａｔｉｏｎ Ｐｒｅｄｉｃｔｉｏｎ ｉｎ Ｂｌｏｃｋ Ｃａｖｉｎｇ ＰＥＮＧ Ｐｉｎｇ⁃ａｎ１，２， ＷＡＮＧ Ｌｉ⁃ｇｕａｎ１，２ （１．Ｓｃｈｏｏｌ ｏｆ Ｒｅｓｏｕｒｃｅｓ ａｎｄ Ｓａｆｅｔｙ Ｅｎｇｉｎｅｅｒｉｎｇ， Ｃｅｎｔｒａｌ Ｓｏｕｔｈ Ｕｎｉｖｅｒｓｉｔｙ， Ｃｈａｎｇｓｈａ ４１００８３， Ｈｕｎａｎ， Ｃｈｉｎａ； ２．Ｄｉｇｉｔａｌ Ｍｉｎｅ Ｒｅｓｅａｒｃｈ Ｃｅｎｔｅｒ， Ｃｅｎｔｒａｌ Ｓｏｕｔｈ Ｕｎｉｖｅｒｓｉｔｙ， Ｃｈａｎｇｓｈａ ４１００８３， Ｈｕｎａｎ， Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ Ｉｎ ｖｉｅｗ ｏｆ ｄｒａｗｂａｃｋｓ ｉｎ ｔｈｅ ｔｒａｄｉｔｉｏｎａｌ ｔｗｏ⁃ｄｉｍｅｎｓｉｏｎａｌ ｍｏｄｅｌ ｆｏｒ ｆｒａｇｍｅｎｔ ｓｉｚｅ ｐｒｅｄｉｃｔｉｏｎ ｉｎ ｂｌｏｃｋ ｃａｖｉｎｇ， ｂａｓｅｄ ｏｎ ｔｈｅ ｓｐａｃｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｒｏｃｋ ｓｔｒｕｃｔｕｒｅ， ｉｔ ｗａｓ ｐｒｏｐｏｓｅｄ ｔｏ ｃｏｎｓｔｒｕｃｔ ａ ｎｅｗ ３Ｄ ｍｏｄｅｌ ｂｙ ｃｕｔｔｉｎｇ ｔｈｅ ｔｈｒｅｅ⁃ ｄｉｍｅｎｓｉｏｎａｌ ｅｎｔｉｔｙ ａｆｔｅｒ ｒｏｃｋ ｊｏｉｎｔｓ ｓｉｍｕｌａｔｅｄ ｗｉｔｈ Ｍｏｎｔｅ Ｃａｒｌｏ ｍｅｔｈｏｄ， ｔｈｕｓ ｒｅｓｕｌｔｉｎｇ ｉｎ ｔｈｅ ｆｒａｇｍｅｎｔ ｓｉｚｅ ｐｒｅｄｉｃｔｉｏｎ． Ｔｈｅｎ ａ ｓｏｆｔｗａｒｅ， ＯＦＰ １．０ （Ｏｒｅ Ｆｒａｇｍｅｎｔｓ Ｐｒｅｄｉｃｔｉｏｎ）， ｗａｓ ｄｅｖｅｌｏｐｅｄ． Ｔｈｅ ａｐｐｌｉｃａｔｉｏｎ ｏｆ ｓｕｃｈ ｔｅｃｈｎｉｑｕｅ ｉｎ ａ ｌａｒｇｅ ｐｏｒｐｈｙｒｙ ｃｏｐｐｅｒ ｍｉｎｅ ｇａｖｅ ａ ｐｒｅｄｉｃｔｉｏｎ ｔｈａｔ ｔｈｅ ｃｕｍｕｌａｔｉｖｅ ｐｅｒｃｅｎｔａｇｅ ｏｆ ｏｒｉｇｉｎａｌ ｂｌｏｃｋ ｗｉｔｈ ｅｑｕｉｖａｌｅｎｔ ｓｉｚｅ ｌａｒｇｅｒ ｔｈａｎ １．６５ ｍ ｗａｓ ４０％， ａｎｄ ｔｈａｔ ｗｉｔｈ ｅｑｕｉｖａｌｅｎｔ ｓｉｚｅ ｌａｇｅｒ ｔｈａｎ １．２６ ｍ ｗａｓ ５９．８５％． Ｂａｓｅｄ ｏｎ ｔｈｅ ｃｏｎｖｅｒｓｉｏｎ ｂｅｔｗｅｅｎ ｅｑｕｉｖａｌｅｎｔ ｓｉｚｅ ａｎｄ ｂｌｏｃｋ ｓｉｚｅ， ｉｔ ｉｓ ｒｅｃｏｍｍｅｎｄｅｄ ｔｈａｔ ｔｈｅ ｏｒｅ ｐａｓｓ ｓｈｏｕｌｄ ｂｅ５ ｍ ５ ｍ ｉｎ ｄｉｍｅｎｓｉｏｎ ｓｏ ａｓ ｔｏ ｐｒｅｖｅｎｔ ｔｈｅ ｏｃｃｕｒｒｅｎｃｅ ｏｆ ｃｌｏｇｇｉｎｇ． Ｓｕｃｈ ｆｉｎｄｉｎｇｓ ｃａｎ ｂｅ ｏｆ ｒｅｆｅｒｅｎｃｅ ａｎｄ ｇｕｉｄａｎｃｅ ｉｎ ｂｌｏｃｋ ｃａｖｉｎｇ ｄｅｓｉｇｎ ｆｏｒ ｕｎｄｅｒｇｒｏｕｎｄ ｍｉｎｉｎｇ ａｎｄ ｐｒａｃｔｉｃａｌ ｃｏｎｔｒｏｌ ｏｆ ｏｒｅ ｄｒａｗｉｎｇ ｉｎ ｍｉｎｅ． Ｋｅｙ ｗｏｒｄｓ ｂｌｏｃｋ ｃａｖｉｎｇ ｍｅｔｈｏｄ； ｆｒａｇｍｅｎｔ ｓｉｚｅ ｐｒｅｄｉｃｔｉｏｎ； Ｍｏｎｔｅ Ｃａｒｌｏ ｍｅｔｈｏｄ；ｊｏｉｎｔ ｐｌａｎｅ ｓｉｍｕｌａｔｉｏｎ 自然崩落法是一种借助矿石自重和地压作用引起 矿石崩落的采矿方法，具有低成本、高效率、作业安全 等特点。 与其他采矿方法相比，其缺点也十分明显，如 技术要求严格、实施风险高等［１－２］。 在该采矿方法的 设计阶段，底部结构尺寸设计、设备选型等与矿岩崩落 块度的分布存在着密不可分的联系［３］。 块度预测的 可靠程度直接影响矿石的回收率、贫化率、采矿系统的 工作效率和工作安全性。 因此，在采用自然崩落法的 矿山对矿岩崩落块度进行预测显得尤为重要。 国内外学者在块度预测方面进行了诸多研究［４－７］。 由于矿石块度形成过程极其复杂，这些研究多是基于 二维状态下进行，不能准确反映真实情况，给采矿设计 和实际生产带来一定风险。 基于此，本文提出基于岩体构造空间分布规律的 统计分析，利用 Ｍｏｎｔｅ Ｃａｒｌｏ 模拟方法，实现了矿岩崩 落块度的三维预测，并以某特大型斑岩铜矿为例，验证 了该系统的有效性。 １ 三维矿岩崩落块度预测基本原理 岩体的破坏将首先沿其体内原有的结构面发展， 通过对岩体结构进行调查，获得有关岩石强度、结构面 产状、节理间距等的统计特征与统计参数，因此就不难 应用计算机来推估未知区域的岩体和结构面特征以及 再现符合这种统计特征的岩体和结构面网络，从而研 ①收稿日期 ２０１４－１２－１６ 基金项目 国家高技术研究发展计划（８６３ 计划）项目（２０１１ＡＡ０６０４０７） 作者简介 彭平安（１９８９－），男，江西宜春人，博士研究生，主要研究方向为数字矿山。 第 ３５ 卷第 ３ 期 ２０１５ 年 ０６ 月 矿矿 冶冶 工工 程程 ＭＩＮＩＮＧ ＡＮＤ ＭＥＴＡＬＬＵＲＧＩＣＡＬ ＥＮＧＩＮＥＥＲＩＮＧ Ｖｏｌ．３５ №３ Ｊｕｎｅ ２０１５ 究矿岩的崩落块度。 １．１ 崩落块度影响因素分析 随着人们对自然崩落法的深入研究，普遍认识到 影响矿岩崩落块度的主要因素有三类一是节理的空 间展布状态，即岩体结构面的空间几何参数；二是地应 力状态，包括原始地应力状态和拉底崩落过程中产生 的次应力场；三是研磨作用，崩落下来的矿石在放矿过 程中由于速度场不同，必然造成矿岩块之间的相互研 磨与碰撞，从而使崩落块体再次破碎［８－９］。 因此，通过 对影响崩落矿岩块度的各主要因素的考察可得出以下 结论① 由于应力作用，岩体破坏的第一步将沿已有 弱面产生；② 结构面的切割将导致岩体形成不同大 小、不同形状的小块集合；③ 结构面的几何形状将决 定切割后的块体分布。 因此，将结构面的倾向、倾角和间距分布规律反演 到待评价岩体中，即可实现对真实岩体构造面的模拟。 １．２ 三维节理面的 Ｍｏｎｔｅ Ｃａｒｌｏ 随机模拟 Ｍｏｎｔｅ Ｃａｒｌｏ 法源于二战时期美国对裂变物质的 中子随机扩散进行模拟。 其基本思想是为解决工程 应用中的各类复杂问题，通过建立与待解问题相似或 相关联的概率模型或概率过程，对概率模型进行随机 模拟或统计抽样，用所得样本得到这些特征量的估计 值，从而近似推估问题的解［１０－１１］。 节理面的 Ｍｏｎｔｅ Ｃａｒｌｏ 模拟则根据构造面的参数分布规律，计算各样本 的参数估计值，以此代替母体参数，然后采用线性同余 法产生随机数进行随机抽样，用以模拟原系统参数。 根据国内外的研究，结构面一般服从正态分布、对 数正态分布和指数分布 ３ 种形式［１２］，对应的伪随机数 发生器分别为 正态分布 Ｎｉ ＝ σ － ２ｌｎｒｉ ｃｏｓ（２πｒｉ＋１） ＋ μ Ｎｉ＋１ ＝ σ － ２ｌｎｒｉ ｓｉｎ（２πｒｉ＋１） ＋ μ { （１） 式中 σ 为正态分布随机变量的均方差；μ 为正态分布 随机变量的均值；Ｎｉ为正态分布随机数；ｒｉ为［０，１］均 匀分布随机数。 对数正态分布 Ｓｉ＝ ｅｘｐ［σ′ｉ（ － ２ｌｎｒｉ） ｃｏｓ（２πｒｉ＋１） ＋ μ′ｉ］ Ｓｉ＋１＝ ｅｘｐ［σ′ｉ（ － ２ｌｎｒｉ） ｓｉｎ（２πｒｉ＋１） ＋ μ′ｉ］ { （２） 式中 Ｓｉ为对数正态分布随机数；μ′和 σ′分别为对数转 换后的分布方差和均值。 指数分布 Ｅｒ＝－ １ Ｅ ｌｎｒｉ（３） 式中 Ｅｒ为服从指数分布的随机数；Ｅ 为指数分布的随 机参数。 ２ 三维矿岩崩落块度预测系统的实现 ２．１ 算法流程 基于上述分析，提出矿岩崩落块度预测系统的基 本算法，其流程如图 １ 所示 7,3C 0;ε ,9330,5MD9 038,κ ;-* ;2 0360 0300 0 0** A3MD BD,0 D30, 9,090 3;,, 090 30 05-136 EBD023 Monte Carlo066 *11136, ,;14 8- 0D3 0,5 **9 N N N N N N Y Y 136;-D/ Y B;3,B53,B C1302 YY Y .B10136CA4 -D0,* 0; 01360C;088, 0; 图 １ 矿岩崩落块度预测系统算法流程 １） 输入结构面统计分布特征参数，执行 Ｍｏｎｔｅ Ｃａｒｌｏ 模拟，形成节理面随机素描文件。 ３２第 ３ 期彭平安等 矿岩崩落块度预测的节理面三维模拟技术及其应用 ２） 确定模拟位置及范围，生成模拟岩体的实体模 型，用三维节理面切割实体模型。 ３） 通过多边形三角剖分方法把切割产生的不封 闭多边形封闭起来。 ４） 删除外表面分割形成的块，剩余块即为模拟产 生的崩落块体。 ５） 计算岩块的体积、面积，最大、中间和最小弦 长，并判别岩块的形状特征。 ６） 绘制图形结果。 ２．２ 系统功能 根据所提算法，开发了块度预测软件 ＯＦＰ １．０ （Ｏｒｅ Ｆｒａｇｍｅｎｔｓ Ｐｒｅｄｉｃｔｉｏｎ １．０）。 系统采用模块化管 理，包括模拟参数输入、三维节理系统模拟、原始和崩 落块度模拟、块体形状特征分析和块体分布统计共 ５ 个模块，各模块的功能为 １） 模拟参数输入用于输入模拟参数，包括节理 参数、节理面力学参数、完整岩石力学参数、地应力参 数等。 ２） 三维节理系统模拟根据输入的模拟参数，利 用 Ｍｏｎｔｅ Ｃａｒｌｏ 方法产生随机素描文件。 ３） 原始和崩落块度模拟主要用于产生原始和崩 落块度文件，并将结果保存。 ４） 块体形状特征分析对上一步生成的块体按特 征规律进行分析。 ５） 块体分布统计统计块体形状和特征规律，绘 制筛上累计百分率曲线。 ３ 工程应用 ３．１ 工程概况 某特大型斑岩铜矿位于我国西南部，远景铜金属 储量 ４３７ 万吨。 该矿构造裂隙发育，岩石蚀变强烈，具 典型的“斑岩型”蚀变分带。 岩浆岩、岩浆侵位的地 层、热液蚀变作用、热液运移和矿质沉淀的构造空间控 制了本矿床的产出。 １＃矿体是其主要矿体，矿体长 ２ ２４０ ｍ，垂深 ２．２３～７６６．３５ ｍ。 矿体呈大透镜状，北西 向展布，空间上呈一北西向展布“马鞍”状，在平面上 为一不规则的“多节葫芦”形。 考虑到矿床厚大但品 位低，同时岩层破碎、完整性差，工程地质条件恶劣，拟 采用自然崩落法开采。 ３．２ 参数设置 软件输入参数主要包括① 节理面倾角、倾向、间 距的分布规律及样本参数估计值；② 现场主应力的大 小及方向；③ 完整岩石的强度指标；④ 不连续相关因 子和节理面持续性系数。 ３．２．１ 结构面参数 采用详细线观测法进行了大量 岩体原位构造调查工作，获得有效节理产状数据 １０ ４３４ 条，根据调查结果，作 Ｓｃｈｍｉｄｔ 极点等密图及玫 瑰花图，见图 ２。 图 ２ 节理产状的极点等密图及玫瑰花图 （ａ） 极点等密图； （ｂ） 玫瑰花图 从图 ２ 可以看出，整体可分为 ３ 组优势节理组，其 倾向范围分别为１２５～１９５，１９６～２４０，３２０～３５９。 优势节理组中倾向倾角的分布情况如表 １ 所示。 表 １ 矿区优势节理分组情况 组 号 样本 数据 百分比 ／ ％ 倾向／ （）倾角／ （） 最小 最大均值方差最小 最大均值方差 １４ ６１７４４．２１２５１９５１６３．８５３５６．４９１９０６５．８１２１６．８０ ２１ ８８０１８．０１９６２４０２１４．２７１２０．６９３９０６９．９４１８０．７５ ３１ ４５３１３．９３２０３５９３４１．６６１０３．００２８９５１．１９５７３．１９ 在利用 Ｐ⁃Ｐ 图对优势节理组倾向倾角进行分布 检验时，正态分布检验 Ｐ⁃Ｐ 图中散点基本位于 ｙ＝ｘ 直 线附近，正态趋势 Ｐ⁃Ｐ 图中散点基本在 ｙ ＝ ０ 上下波 动，波动范围不大于０．１，拟合性较好，因此可认为 ３ 组优势节理组倾向倾角均服从正态分布。 ３．２．２ 地应力测量 在工程地质调查基础上，结合现 场施工条件，在 ３ ７２０ ｍ、３ ４５０ ｍ 两个中段选择了 ６ 个 测点（其中以首采中段 ３ ７２０ ｍ 作为研究重点，布置了 ４ 个测点），利用国际先进的瑞典 ＬＵＴ 地应力测定仪 对矿区的地应力场的变化及分布规律进行测试研究， 测量结果列于表 ２ 中。 ４２矿 冶 工 程第 ３５ 卷 表 ２ 地应力测量结果 测点 编号 埋深 ／ ｍ 最大主应力中间主应力最小主应力 大小 ／ ＭＰａ 方位 ／ （） 倾角 ／ （） 大小 ／ ＭＰａ 方位 ／ （） 倾角 ／ （） 大小 ／ ＭＰａ 方位 ／ （） 倾角 ／ （） １２０３１７．０９２５６．５１０．８９ １１．８５ ３４７．４５４．９５７．７８１０１．５４ ７８．０１ ２２３８１５．９５ ２８５．１９５．７３１０．３１ １９０．５６ ３８．８２７．０６２２．２０５０．５９ ３１９３１１．６０ ２４７．９１９．７７６．３７３４４．７１ ３７．４８５．４２１４６．８７ ５１．１４ ４２０７１５．７９９８．０９１３．２２９．５０１９８．６２ ３７．８８２．６４３５２．３７ ４９．０７ ５２１４１３．３３９４．１７２９．４１５．６９１９１．７７ １３．２１２．１９３０３．１３５７．２ ６２０９１２．６３８６．５９２２．２１７．４０１８１．９４ １２．８６０．０３２９９．８６ ６３．９８ ３．２．３ 岩石物理力学参数测试 为准确预测崩落块 度大小，还需掌握主要岩体的物理力学性质。 在本次 工作中，制作了与采矿工程稳定性关系较为密切的 ４ 种主要岩石（石英二长斑岩、闪长玢岩、大理岩和角 岩）试样共计 ７０ 个，对该矿岩石物理力学参数进行了 测试，结果如表 ３ 所示。 表 ３ 岩石物理力学参数测试结果 岩石 名称 单轴抗压强度 ／ ＭＰａ 抗拉强度 ／ ＭＰａ 内聚力 ／ ＭＰａ 内摩擦角 ／ （） 石英二长斑岩１２７．９６７．０７２２．０６４７．３１ 闪长玢岩１８５．６７１２．２９２２．７０４１．１８ 大理岩１２６．３８８．６５２７．１７４４．５１ 角 岩１９２．５１１４．７２３２．４９４０．０３ ３．３ 矿区整体块度分布预测 以上述参数为输入，对该矿矿岩崩落块度进行预 测。 软件进行切割完成后产生的模拟块体形状如图 ３ 所示，图中 ε 值用于区分物体的细长性，κ 用来表述物 体的容积系数，Ｄ２５，Ｄ５０，Ｄ７５分别表示占总体积 ２５％、 ５０％和 ７５％时块体体积的大小。 对产生的切割块体结 果进行统计，得到的矿岩块度体积区间分布如图 ４ 所 示，图中 Ｃ 代表立方体，ＣＥ 代表立方⁃细长体，ＰＣ 代表 扁平⁃立方体，Ｐ 代表扁平体，ＥＰ 代表细长⁃扁平体，Ｅ 代表细长体。 从图 ４（ｂ）可以看出，块体筛下体积的四 分位点分别为０．１６４，０．４０３ 和０．８５３ ｍ３。 从图４（ｃ）可 以看出，当等效尺寸在 ０．５～３ ｍ 之间时，等效尺寸曲 线下降较快，这说明块体等效尺寸落在这个区间体积 百分比较大。 等效尺寸大于 １．６５ ｍ 的块体筛上累积 百分比为 ４０％，大于 １．２６ ｍ 的块体筛上累积百分比为 ５９．８５％。 根据等效尺寸与块体体积的转换关系，建议 12 345 678 910 κ 0.0010.010.07 0.110.20.37 0.5 0.6 0.8 ε D75 ab 图 ３ 产生的块体切割结果及形状 （ａ）块体切割结果； （ｂ） 不同块体体积的形状分布 18 16 14 12 10 8 6 4 2 0 D75D25D50D50D75 100 80 60 40 20 0 100 80 60 40 20 00 1 2 3 4 5 6 0.0001 0.001 0.01 0.1 1 10 3090 30m3 3**m CD30,530,5 D25 D250.164 D500.403 D750.853 D50 D75 1.65 1C 3E 5P 2CE 4EP 6PC 12 3 4 56 2 3 4 5 6 12 3 4 56 123 4 6 b c a 图 ４ 矿岩崩落块度分布曲线 放矿口的尺寸设置为５ ｍ ５ ｍ， 此时放矿时堵塞的几 率较小。 ３．４ 矿体、上盘和下盘区域块度分布预测 应用所开发系统对矿体、上盘和下盘区域崩落块 度分别进行了预测，等效尺寸分布曲线如图 ５ 所示。 100 80 60 40 20 00 1 2 3 4 5 6 3**m 30,5 1.51 1.61 1.91 37 图 ５ 崩落块体等效尺寸分布曲线 ５２第 ３ 期彭平安等 矿岩崩落块度预测的节理面三维模拟技术及其应用 从图 ５ 可见，矿体崩落产生的块度较小，下盘次 之，上盘较大。 下盘块度较小，不利于开拓工程的稳定 性，但上盘块度较大，不易引起矿石贫化。 块体筛上累 积百分比为 ４０％时，矿体、上盘和下盘岩体的等效尺 寸分别为１．５１，１．９１ 和１．６１ ｍ。 因此，设计时可根据放 矿口的具体位置及预测结果对放矿口尺寸进行调整。 ４ 结 语 １） 基于岩体构造空间分布规律的统计分析，利用 Ｍｏｎｔｅ Ｃａｒｌｏ 模拟方法，提出了矿岩崩落块度预测算 法，开发了相应软件，实现了矿岩崩落块度的三维 预测。 ２） 应用所开发系统对某特大型斑岩铜矿自然崩 落法矿岩崩落块度进行了预测。 结果表明，等效尺寸 大于１．６５ ｍ 的块体筛上累积百分比为４０％，大于１．２６ ｍ 的块体筛上累积百分比为 ５９．８５％。 因此建议放矿口 的尺寸为５ ｍ ５ ｍ， 此时放矿时堵塞的几率较小。 按 矿体、上盘和下盘区域分别进行预测的结果表明矿体 崩落产生的块度较小，下盘次之，上盘较大。 ３） 研究成果为矿山采矿方法设计和工程中的放 矿控制提供了重要依据，也为类似自然崩落法矿山实 现安全大规模开采提供了参考。 参考文献 ［１］ 袁海平，曹 平． 我国自然崩落法发展现状与应用展望［Ｊ］． 金属 矿山，２００４（８）２５－２８． ［２］ 于少峰，吴爱祥，韩 斌． 自然崩落法在厚大破碎矿体中的应用 ［Ｊ］． 金属矿山，２０１２（９）１－４． ［３］ 褚洪涛，章 光，李爱兵． 某矿区矿岩可崩性评价及块度分布预测 研究［Ｊ］． 矿业研究与开发，２０１０（３）１５－１７． ［４］ Ｗｈｉｔｅ Ｄ Ｈ． Ｐｒｅｄｉｃｔｉｎｇ ｆｒａｇｍｅｎｔａｔｉｏｎ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ａ ｂｌｏｃｋ ｃａｖｉｎｇ ｏｒｅｂｏｄｙ［Ｄ］． Ａｒｉｚｏｎａ Ｔｈｅ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ａｒｉｚｏｎａ，１９７７． ［５］ Ｍｕｋｈｅｒｊｅｅ Ａ， Ｍａｈｔａｂ Ａ． Ｓｉｚｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｏｒｅ ｆｒａｇｍｅｎｔｓ ｉｎ ｂｌｏｃｋ ｃａｖｉｎｇ［Ｃ］∥Ｐｒｏｃ １３ｔｈ Ｗｏｒｌｄ Ｍｉｎｉｎｇ Ｃｏｎｇｒｅｓｓ， Ｓｔｏｃｋｈｏｌｍ，１９８７． ［６］ 董卫军． 矿石崩落块度的三维模型与块度预测［Ｊ］． 矿冶，２００２， １１（２）１－３． ［７］ 王家臣，熊道慧，方君实． 矿石自然崩落块度的拓扑研究［Ｊ］． 岩 石力学与工程学报，２００１，２０（４）４４３－４４７． ［８］ 王李管，潘长良． 矿石块度影响因素的定量分析研究［Ｊ］． 中南矿 冶学院学报，１９９２，２３（５）５１４－５１９． ［９］ 冯兴隆，李 德，王李管，等． 矿岩块度预测影响因素的分析研究 ［Ｊ］． 矿冶工程，２０１１，３１（３）８－１２． ［１０］ 曲双石，王会娟． Ｍｏｎｔｅ Ｃａｒｌｏ 方法及其应用［Ｊ］． 统计教育，２００９ （１）４５－５５． ［１１］ 邵 伟． 蒙特卡洛方法及在一些统计模型中的应用［Ｄ］． 济南 山东大学数学学院，２０１２． ［１２］ 吴 琼． 复杂节理岩体力学参数尺寸效应及工程应用研究［Ｄ］． 武汉中国地质大学（武汉）工程学院，２０１２． 􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎􀪎 􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃 􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃 􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃 􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀥃􀧃 􀧃 􀧃 􀧃 矿冶工程杂志 ２０１５ 年征订启事 矿冶工程（双月刊）创刊于 １９８１ 年，由中国金属学会、长沙矿冶研究院有限责任公司主办。 本刊 是中国期刊方阵“双效期刊”、全国中文核心期刊、中国科学引文数据库（ＣＳＣＤ）及中国学术期刊综 合评价数据库来源期刊、中国核心学术期刊（ＲＣＣＳＥ），是集学术性和技术性于一体的综合性刊物，已 被中国知网（ＣＮＫＩ）、万方数据库、重庆维普资讯、台湾华艺数据库等全文收录，是国外多家知名检索机 构的检索对象。 矿冶工程杂志主要栏目为采矿、选矿、冶金、材料、矿冶行业企业管理等。 读者对象是采矿、选 矿、冶金、材料、地质、煤炭、化工等系统的有关生产人员、院校师生和管理人员。 矿冶工程面向国内外公开发行，读者分布在 １７ 个国家与地区，机构用户 ４０００ 多个。 矿冶工 程为国际标准大 １６Ｋ 本，１２０ 版。 国内统一刊号ＣＮ４３－１１０４／ ＴＤ，国际标准刊号ＩＳＳＮ０２５３－６０９９，邮 发代号４２－５８。 真诚欢迎新、老订户向全国各地邮局订阅本刊，也可直接向编辑部订阅。 定价 １０ 元， 全年 ６０ 元。 矿冶工程编辑部承接彩色、黑白及文字广告业务，欢迎各企事业单位来电来函联络。 地 址 湖南省长沙市麓山南路 ９６６ 号矿冶工程编辑部联 系 人 黄小芳 邮 编 ４１００１２开户名称 长沙矿冶研究院有限责任公司 电 话 （０７３１）８８６５７０７０／８８６５７１７６／８８６５７１７３开户银行 长沙市工商银行左家垅分行 传 真 （０７３１）８８６５７１８６帐号 １９０１０１３００９０１４４５０２７６ Ｅ－ｍａｉｌ ｋｕａｎｇｙｅｇｏｎｇｃｈｅｎｇｚｚ＠ １６３．ｃｏｍ ６２矿 冶 工 程第 ３５ 卷