改进蒙特卡洛算法及其在高阶段充填采矿中的应用.pdfVIP

第23卷第4期 湖南科技大学学报(自然科学版) V01．23No．4 l JournalofHunan of Science Dec．2()()8 2008年2月 Edition) UniversityScienceTechnology(Natural 改进蒙特卡洛算法及其在高阶段充填采矿中的应用 张义平，，翁春林2 (1．贵州大学矿业学院，贵州贵阳5500032．昆明冶金高等专科学校，云南昆明650033) 摘要：用简便、经济、可靠的方法产生随机数是蒙特卡洛法应用的关键．用Logistic迭代方程产生随机数，不仅简单可靠，而且 符合要求。改进了蒙特卡洛算法，为工程可靠性计算提供了一种新的方法．对高阶段充填体进行了力学分析，采用微元应力分析法， 推导了分层充填力学计算公式．用改进蒙特卡洛算法分析了高阶段充填体稳定性，并应用于某铜矿5号采场高阶段充填设计．研究 结果表明，充填体暴露高度越高，充填体内应力越大，充填体暴露长度越小，充填体向上下盘围岩转移应力的效果越显著，因此缩短 采场长度有刺于充填体稳定性和降低充填成．图6，表2，参14． 关键词：高阶段充填体；蒙特卡洛算法；可靠性；Logistic迭代方程 中图分类号：TD853．34 文献标识码：A 高阶段采矿大都采用尾砂嗣后充填(矿房胶结充 种依据统计抽样并利用计算机研究随机变量的计算 填、矿柱非胶结充填)，是一高效采矿技术【11．采矿过程 方法．蒙特卡洛法分析工程可靠度时，产生的随机数 中，胶结充填体垂直暴露高度达60．140m，侧向暴露 面积在3000．7000 m2．对高阶段充填体稳定性研究 通常采用确定性方法(即定值法)，该方法得到的安全 速度快，节省计算机时间；(5)达到重复(循环)之前， 系数指标体现了充填体稳定性的一个重要方面，但尾 能产生足够使用的随机数．历史上曾采用过平方取中 砂充填体是一种复杂的非线性力学介质【2】，工程实践 法、乘法取中值、常数乘子法、裴波那契法等方法，都 表明㈨，有诸多因素影响胶结充填体强度，充填体力 因其固有的缺点而被淘汰[91，目前广泛使用线性同法 学参数存在不确定性与随机性，因此用可靠性理论研 (可分为乘法线性同余法、混合线性同余法、二次同余 究高阶段充填体稳定性更符合客观实际． 法)也存在计算工作量大的缺点．混沌(Chaos)研究非 Carlo)模 在目前可靠性分析中，蒙特卡洛(Monte 线性动力系统中出现的内在随机性，具有初始值敏感 拟法是一相对精确的方法唧，它受问题条件限制的影 响较小，适应性强．但蒙特卡洛法分析工程可靠度时， 特性可简便地产生大量遍历、不重复分布的随机数， 必须首先形成系列可靠的随机数，能否有简便、经济 Logistic迭代方程为蒙特卡洛模拟计算提供了一种新 的方法产生随机数，是蒙特卡洛法模拟计算的关键， 的研究思路． 为此本文利用Logistic迭代方程改进蒙特卡洛算法， 1．1 Logistic迭代方程产生随机数方法 并应用于高阶段充填采矿工程实践． 1 改进蒙特卡洛算法 (1) 芝。兰三gi哼s。ti，c雩竺j冀：^水(0，．，． 蒙特卡洛法又称随机模拟法或统率试验法，是一 收稿日期：2008—09—10 基金项目：贵州省教育厅重点科研基金项目(黔教科2006203) 作者简介：张义平(1970-)。男，湖南邵东人，博士，副教授，主要从事岩