乳山金矿深部采场结构参数的数值模拟精.docVIP

乳山金矿深部采场结构参数的数值模拟 董卫军 1, 苏永华 2, 孙玉堂 3, 陈国平 3, 宋文志 3 (1. 北京矿冶研究总院 , 北京 　 100044;2. , 3. 山东金洲矿业集团 , 山东 摘 　 要 :。利用采场结构 , 。 关键词 ; ; 尾砂充填 ; 矿房参数 ; 采场稳定性 1　 问题的提出 山东省乳山金矿已有 30多年的开采历史 , 主要 采矿方法为浅孔留矿法 。 根据多年的采矿实践活动 和地压观察 , 在 -385m 水平以上的开采过程中采场 和巷道均是安全的 。但在 -385m 水平开采时出现 片帮和岩柱剥落现象 。随着采深的增加 , 采矿方法 的适应性成为地下结构稳定的主要影响因素之一 。 根据大量的研究和分析 , 从安全和经济角度考虑 , 在 -485m 以下采用分级尾砂充填采矿法是合适的 , 但 采矿单元的具体参数必须用适当的方法优化确定 。 2　 工程概况 2. 1　 矿床开采地质条件 乳山金矿矿区有 10个矿体 , 分为主矿体 、 次矿 体和零星矿体 。其中 Ⅱ 号矿体为主矿体 。 Ⅰ 、 Ⅳ 号 矿体为次要矿体 。 Ⅱ 号矿体呈脉状产出 , 以含金黄 铁矿石英脉为主 。赋存标高在 120m ~-870m 之 间 , 长 140m ~470m ; 厚度为 1. 5m ~2. 5m , 矿体平均 厚度为 2m 。 区内矿体围岩的蚀变较发育 。根据矿 体及围岩地质条件 , Ⅱ 号矿体围岩可分为 3个区 。 Ⅰ 区位于矿区中部 , 厚为 3~5m 。矿体及近矿围岩 受多组节理裂隙的切割 , 是一复杂的软硬相间的复 合岩体 , 由细块状结构岩体组成 。 Ⅱ 区位于 Ⅰ 区的 两侧 , 厚为 10~15m , 主要是强至中等蚀变的绢云母 钾化混合花岗岩 , 系完整中等块状结构岩体 , 故 Ⅱ 区 为较稳定岩体 。 Ⅲ 区位于 Ⅱ 区的两侧 , 厚为 85~ 190m , 主要为混合花岗岩夹透镜状的黑云母斜长片 麻岩 , 系完整的大块状结构岩体 , 故 Ⅲ 区基本为稳定 的岩体 。 根据国际流行的岩体分级方法 , 各区岩体 稳定性分级结果如表 1。 表 1　 矿岩稳定性 CSIR 分类 分类参数 Ⅰ 区岩体 Ⅱ 区岩体 Ⅲ 区岩体 σ c (MPa (A 1 82. 598. 2147. 5 评分值 71012 RQD (A 2 909698 评分值 171819 节理间距 (A 3 0. 280. 300. 65 评分值 121523 节理条件 (A 4 表面微粗糙 表面较粗糙 , 坚硬 表面粗糙 , 坚硬 评分值 101215 地下水 (A 5 潮 较干燥 干燥 评分值 568 总分值 516177 分级级别 Ⅲ Ⅱ Ⅱ 稳定性描述 中等稳固 稳固 稳固 　　 矿区发育较大的断层有 5组 F 1, F 2, F 3, F 4, F 5, 其中对采矿影响较大的是 F 1断层 。 2. 2　 矿区地应力 根据国家地震局地壳应力研究所在 -195m 水 平和 -235m 水平及中科院武汉岩土力学研究所在 -435m 水平 、 -485m 水平 、 -450m 水平及 -535m 三个水平的地应力测试结果分析 , 矿区地应力场为 : σ x =σx 0+λγΗ σ y =σy 0+λγΗ σ z =σz 0+γΗ (1 式中 :σx 0=8. 95MPa , σy 0=5. 54MPa , σz 0=1. 06MPa , 侧压系数 λ取 0. 668, 岩体容重 γ取 2644kg/m 3。 第 2卷 第 2期 Vol. 2　 No. 2 　　　　　　　　　　　　　　　 采 矿 技 术 Mining Technology 　　　　　　　　　　　　　　　 2002年 6月 J une 　 2002 3矿房长度数值模拟 3. 1　 模拟程序选择 数值模拟方法有多种 , 如有限元 、 离散元 、 边界 元 、 无界元 、 流行元等 。 适合于岩土 、 采矿 、 地质工程 的数值试验程序有很多 。 国内外岩土工程界广为应 用的 FLAC 3D 是有限差分程序 。 FLAC 3D 采用的 “显 式拉格朗日” 算法和 “ 混合 -离散分区” 技术 , 能够非 常准确地模拟材料的塑性破坏和流动 。 成刚度矩阵 , 所以 , FLAC 3D 程问题的理想工具 。 , 到广泛应用 , 。 由于岩石是一种脆性材料 , 当外荷载达到岩石 强度后 , 材料发生断裂破坏 , 产生弱化现象 , 应属于 弹塑性体 。 在 FLAC 3D 中 , 一般对于弹塑性材料 , 判 断其破坏与否的基本准则有 Drucker Prager 准则和 Mohr -Coulomb 准则 。根据大量的理论分析和实 验结果 , 其中 Mohr -Cou