川贫矿体三维节理网络的数值模拟.pdf
金川贫矿体三维节理网络的数值模拟 董卫军 (北京矿冶研究总院, 北京 “ “ “ ) 摘要运用计算机模拟技术, 采用蒙特卡罗模拟方法, 结合现场实测节理几何参数, 综合考虑组优势节理和占一定比例 的零散节理, 模拟金川铜镍矿贫矿体内的三维节理分布。由模拟得出节理频度、 节理半径、 节理面与剖面的交线, 节理圆盘中心到 直线的距离, 生成节理网络。模拟结果可作为自然崩落法设计与实施的依据。 关键词自然崩落法; 节理; 计算机模拟 中图分类号 “ ( ’ * ) 式中 .服从负指数分布随机数;2负指数分布 的参数。 () 正态分布复合抽样方法。 ,,.. 5 ; ( ’ * , ) , (. ’ * .) .,.. 5 ; ( ’ * , ) , 6 ; (. ’ * .) 式中 ,, ..服从正态分布的随机数; ’ * ,, ’ * ..[1,,] 区间均匀分布的随机数。 () 对数正态分布抽样方法。由正态分布与对 数正态分布的关系, 可以得到对数正态分布的随机 数为 ,/ ,; ./ . 第 卷第期 .11.年月 有色金属 3 - 3 A B B - C DE A “ 0 F D G 5 7 7 / 7 6 7 A B 7 9 9 C D E 9 7 对于某一区域内岩体, 由节理频度可以求得在 该区域内节理的总体数目。再由各组节理占总体的 百分比和抽样函数, 得到该区域内岩体节理裂隙的 真实分布情况。节理裂隙在岩体中的真实分布状况 是节理网络生成和三维块度生成的基础。 根据蒙特卡罗法生成的节理系统, 运用空间几 何知识, 得到节理系统与空间某平面的交线段, 平面 内的这些线段就组成了该平面内的节理网。 例如, 对于某个剖面* 6(* 6为常数) , 节理 系统中的某条半径为, 中心坐标为 ( “,““, “)的 节理, 节理面方程可表示为 7 8 9 “ 8 80“ (/ “) , 3 . 1 3 . 0 , 4 A ’B , C , . 6 3 4 , * 0 6 3 * 0 1 . , , C , * 5 * 3 4 . 6 /D 4 6 /5 E 5 , C 5 5 * 3 7 4 6 3 AB 4 7 4 7 89 7 B 3 9 6 - 8 “,* 3 4 5 4 7 89 8 8,, 4 7 9) 1 3 - 1 5 Q * * , - . / 0 1 , D 3 . . , * . , D E 5 5 C 4 3 , - 1 6 C F C 4 3 *A , 1 ’S 1 , Q * * , - . / 5 D 3 * , D E 5 * 1 , Q * , * 5 , 5,. 3 5, * , . 5 , 6 , 4 * , 0 1 , Q * F 4 3 * , 3 * 5 , 6 *F 4 3 * ,,3 * 1 , 5 3 * 6 , 0 . CQ * 5 6 6 , * , . 1 , 4 * , , . 7 , 0 . C5 C 4 3 * ’S 1 , . , 5 4 6 3 *D , 5 , 3 5 1 ,D 3 5 5 0 . , 5 *3 * F . 3 6 6 , 0 6 3 7 * ’ , 4* -6 3 7 * C , 1 ;Q * ;6 C F , . 5 C 4 3 * OL 有色金属第P 8卷 万方数据