黄金科学技术
img

QQ群聊

img

官方微信

高级检索
作者投稿 专家审稿 编辑办公

黄金科学技术 ›› 2019, Vol. 27 ›› Issue (2): 189-198.doi: 10.11872/j.issn.1005-2518.2019.02.189

• 采选技术与矿山管理 • 上一篇    下一篇

基于PFC2D-DFN的自然崩落法数值模拟研究

方传峰1,2,3(),王晋淼1,4,李剡兵5,贾明涛1,4   

  1. 1. 中南大学资源与安全工程学院，湖南 长沙 410083
    2. 长沙矿山研究院有限责任公司，湖南 长沙 410012
    3. 金属矿山安全技术国家重点实验室，湖南 长沙 410012
    4. 中南大学数字矿山研究中心，湖南 长沙 410083
    5. 云南迪庆有色金属有限责任公司，云南 香格里拉 674400
  • 收稿日期:2018-02-01 修回日期:2018-05-30 出版日期:2019-04-30 发布日期:2019-04-30
  • 作者简介:方传峰（1990-），男，山东泰安人，硕士研究生，从事矿山开采数值模拟分析研究工作。15084850237@163.com
  • 基金资助:
    国家重点研发计划项目“深部集约化开采生产过程智能管控技术”(编号：2017YFC0602905)

Numerical Simulation Research of Natural Caving Method Based on PFC2D-DFN

Chuanfeng FANG1,2,3(),Jinmiao WANG1,4,Shanbing LI5,Mingtao JIA1,4   

  1. 1. School of Resources and Safety Engineering，Central South University，Changsha 410083，Hunan，China
    2. Changsha Institute of Mining Research Co. ，Ltd. ，Changsha 410012，Hunan，China
    3. State Key Laboratory of Safety Technology of Metal Mines，Changsha 410012，Hunan，China
    4. Center of Digital Mine Research，Central South University，Changsha 410083，Hunan，China
    5. Yunnan Diqing Nonferrous Metal Co. ，Ltd. ，Shangri-La 674400，Yunnan，China
  • Received:2018-02-01 Revised:2018-05-30 Online:2019-04-30 Published:2019-04-30

RichHTML

9

PDF (PC)

818

摘要：

采用二维颗粒流软件PFC2D构建矿山模型，引入符合矿山实际节理产状的DFN网络，基于室内岩石试验及矿山实际崩落情况，标定模拟所用细观参数。以某工程应用为实例，模拟分析崩落破坏机理与崩落区演化趋势。结果表明：裂纹发起于推进面与顶板交界处，在拱形区域内密集扩展、贯通，矿岩破裂以节理拉伸破坏为主，矿堆以岩石剪切破坏为主；崩落区前期呈拱形演化，崩透地表后引起邻近高陡边坡滑动；部分拉底仅造成新开挖顶板上方局部崩落。该研究结果可对矿山崩落采矿起到指导作用。

关键词: 自然崩落法, 数值模拟, PFC2D软件, DFN, 拉底, 岩石破裂, 力学分析, 崩落演化

Abstract:

Surface subsidence and rock burst may occur during mining in mines using natural caving method，which affect the safety of mine production and the life safety of underground workers.Therefore，it is necessary to research the evolution law of caving area with the advance of undercutting.The research process of this paper is as follows.Firstly:using two-dimensional discrete element software PFC2D to construct a geometric model consistent with the actual topography of a mine.Secondly，according to the complete rock parameters obtained from laboratory mechanical experiments，by means of simulating physical experiments，the meso-parameters are continuously adjusted to approximate the actual physical parameters，the meso-mechanical parameters used in PFC2D are calibrated.In addition，according to the statistic law of joint sets distribution obtained on the spot，choose an appropriate joint generation method，make use of discrete fracture network DFN to reconstructe a joint model consistent with the statistical law.Then，based on the caving height of a certain period obtained by mine drilling television，the simulation results are approximated to the monitoring values by dichotomy method，and the corresponding joint strength parameters are back analyzed.Finally，the whole mine model is coupled with the joint model，and the particles and joints are assign property.After stress balance，the excavation is simulated.Statistic simulation results，the conclusions are as follows:Cracks first develop at the intersection of the undercutting advance working face and the roof，then extend，develop and penetrate in the arch region，resulting in crushing and collapse of the original rock.The main failure types of ore and rock are joint tensile failure，while shear failure often occurs in caved ore and rock.In the early stage，the caving area develops arch upward，which causes the adjacent steep slope to slide to the valley after penetrating the surface.In some undercutting steps，the caving area only occurs in the local area above the newly excavated undercutting region，which is consistent with the change trend of caving quality and roof height.

Key words: natural caving method, numerical simulation, PFC2D software, DFN, undercutting, rock failure, stress analysis, caving evolution

中图分类号: 

  • TD853.36

图1

矿山地形及地形轮廓线"

图2

PFC2D矿山模型"

表1

优势节理组产状统计结果"

组号占比/%倾向/（°）倾角/（°）
最大值最小值均值方差分布规律最大值最小值均值方差分布规律
136.36166120143.17126.89正态分布417971.7821.75正态分布
232.74196239216.2195.13正态分布447971.8619.83正态分布
35.41329306317.0735.58正态分布517569.2221.58正态分布
43.99314839.9923.69正态分布567770.5128.00正态分布

表2

节理间距与连续度统计结果"

项目最小值最大值均值标准差分布规律
节理间距0.001.380.0930.15对数正态分布
节理连续度0.0510.612.3500.37正态分布

表3

DFN产生所用参数"

节理分组节理与坑道单位长度交点数/个DFN节理倾角/（°）
最小值最小值均值标准差
10.36407071.784.67
20.33100136108.144.45
30.05104130110.784.64
40.04567770.525.29

图3

DFN生成离散节理局部放大图"

图4

DFN节理间距统计"

图5

单轴压缩试验曲线图"

图6

三轴伺服压缩试验曲线图"

表4

岩石实际物理力学参数与岩石模拟物理力学参数对照表"

对比参数试验结果模拟值误差/%
弹性模量/GPa54.5854.290.53
泊松比0.270.263.70
抗压强度/MPa127.96125.931.60
黏聚力/MPa22.0624.8312.50
摩擦角/（°）47.849.002.50

表5

模拟采用细观参数"

参数类型参数名数值
密度/（kg·m-3）2 700
线性参数弹性模量/GPa48.0
刚度比1.0
摩擦系数1.0
平行粘结参数弹性模量/GPa9.0
刚度比2.5
抗拉强度/MPa75.0
黏聚力/MPa25.0
摩擦角/（°）75.0
半径系数1.5

图7

节理参数标定流程图"

表6

节理参数"

参数名称数值参数名称数值
法向刚度/GPa6.0摩擦角/（°）25
切向刚度/GPa3.0抗拉强度/MPa0.7
黏聚力/MPa7.0

图8

第一步拉底裂纹扩展"

图9

裂纹扩展与破裂类型统计"

图10

崩落区形态演化特征"

图11

崩落区统计图"

1 冯兴隆.自然崩落法矿岩工程质量数字化评价及模拟技术研究[D].长沙：中南大学，2010.
FengXinglong.Study on Digital Assessment and Simulation Technology at Rock Mass Engineering Quality of Block Caving[D].Changsha：Central South University，2010.
2 安龙，徐帅，李元辉，等.基于多方法联合的崩落法崩矿步距优化[J].岩石力学与工程学报，2013，32（4）：754-759.
AnLong，XuShuai，LiYuanhui，et al.Optimization of rate of advance during ore breaking of caving method based on multi-method joint application[J].Chinese Journal of Rock Mechanics and Engineering，2013，32（4）：754-759.
3 BradyB H G，BrownE T.Rock Mechanics for Underground Mining[M].Beijing：Science Press，2014：440-441.
4 WooK S，EberhardtE，ElmoD，et al.Empirical investigation and characterization of surface subsidence related to block cave mining[J].International Journal of Rock Mechanics and Mining Sciences，2013，61（61）：31-42.
5 解世俊.金属矿床地下开采[M].北京：冶金工业出版社，1979：202-203.
XieShijun.Underground Mining of Metal Deposit[M].Beijing：Metallurgical Industry Press，1979：202-203.
6 赵文.岩石力学[M].长沙：中南大学出版社，2014：191-192.
ZhaoWen.Rock Mechanics[M].Changsha：Central South University Press，2014：191-192.
7 孔令海，姜福兴，杨淑华，等.基于高精度微震监测的特厚煤层综放工作面顶板运动规律[J].北京科技大学学报，2010，32（5）：552-558，588.
KongLinghai，JiangFuxing，YangShuhua，et al.Movement of roof strata in extra-thick coal seams in top-coal caving mining based on a high precision micro-seismic monitoring system[J].Journal of University of Science and Technology Beijing，2010，32（5）：552-558，588.
8 薛东杰，周宏伟，任伟光，等.浅埋煤层超大采高开采柱式崩塌模型及失稳[J].煤炭学报，2015，40（4）：760-765.
XueDongjie，ZhouHongwei，RenWeiguang，et al.Instability of pillar collapse model and generation of cracks in thick coal seam mining at shallow depth[J].Journal of China Coal Society，2015，40（4）：760-765.
9 CaoS，DuC F，MuC P，et al.UDEC-based modelling of mining surface movement due to transforming from block caving to sublevel filling and its law verification[J].Rock and Soil Mechanics，2015，36（6）：1737-1743，1751.
10 王涛，盛谦，熊将.基于颗粒流方法自然崩落法数值模拟研究[J].岩石力学与工程学报，2007，26（增2）：4202-4207.
WangTao，ShengQian，XiongJiang.Research on numerical simulation based on particle flow method[J].Chinese Journal of Rock Mechanics and Engineering，2007，26（Supp.2）：4202-4207.
11 王连庆，高谦，王建国，等.自然崩落采矿法的颗粒流数值模拟[J].北京科技大学学报，2007，29（6）：557-561.
WangLianqing，GaoQian，WangJianguo，et al. Numerical simulation of natural caving method based on particle flow code in two dimensions[J].Journal of University of Science and Technology Beijing，2007，29（6）：557-561.
12 曾庆田，刘科伟，严体，等.基于多数值模拟方法联合的自然崩落法开采研究[J].黄金科学技术，2015，23（1）：66-73.
ZengQingtian，LiuKewei，YanTi，et al.Study on natural caving mining method based on multi-numerical simulation method[J].Gold Science and Technology，2015，23（1）：66-73.
13 VyazmenskyA，ElmoD，SteadD.Role of rock mass fabric and faulting in the development of block caving induced surface subsidence[J].Rock Mechanics and Rock Engineering，2010，43（5）：533-556.
14 WooK S，EberhardtE，RabusB，et al.Integration of field characterisation，mine production and InSAR monitoring data to constrain and calibrate 3-D numerical modelling of block caving-induced subsidence[J].International Journal of Rock Mechanics and Mining Sciences，2012，53（9）：166-178.
15 ElmoD，SteadD.An integrated numerical modelling-discrete fracture network approach applied to the characterisation of rock mass strength of naturally fractured pillars[J].Rock Mechanics and Rock Engineering，2010，43（1）：3-19.
16 黄彦华，杨圣奇，鞠杨，等.岩石巴西劈裂强度与裂纹扩展颗粒尺寸效应研究[J].中南大学学报（自然科学版），2016，47（4）：1272-1281.
HuangYanhua，YangShengqi，JuYang，et al.Study on particle size effects on strength and crack coalescence behavior of rock during Brazilian splitting test[J].Journal of Central South University（Science and Technology），2016,47（4）：1272-1281.
17 刘波，金爱兵，高永涛，等.基于分形几何理论的DFN模型构建方法研究[J].岩土力学，2016，37（增1）：625-630,638.
LiuBo，JinAibing，GaoYongtao，et al.Construction method research on DFN model based on fractal geometry theory[J].Rock and Soil Mechanics，2016，37（Supp.1）：625-630，638.
18 Itasca Consulting Group.PFC2D User’s Manual （Version5.0） [M].Minneapolis：Itasca Consulting Group, 2014.
19 李坤蒙，李元辉，徐帅，等.PFC2D数值计算模型微观参数确定方法[J].东北大学学报（自然科学版），2016，37（4）：562-567.
LiKunmeng，LiYuanhui，XuShuai，et al.Method to determine microscopic parameters of PFC2D numerical model[J].Journal of Northeastern University（Natural Science），2016，37（4）：562-567.
[1] 何祥锐, 邱贤阳, 史秀志, 李小元, 支伟, 刘军, 王远来. 基于非线性弹性地基梁的地下矿山充填开采覆岩移动规律研究[J]. 黄金科学技术, 2024, 32(4): 640-653.
[2] 王史文, 王晟, 凌炜佳, 朱忠华, 杨鑫源, 沙文忠, 杨家宏. 基于无人机航测的自然崩落法采矿地表塌陷规律研究[J]. 黄金科学技术, 2024, 32(4): 654-665.
[3] 虞云林, 侯克鹏, 杨八九, 程涌, 卢泰宏, 张楠楠. 云锡高峰山矿段矿柱回采方案研究[J]. 黄金科学技术, 2024, 32(3): 445-457.
[4] 李波, 温晨, 史秀志. 高应力扇形中深孔采场边帮控制爆破参数优化[J]. 黄金科学技术, 2024, 32(3): 511-522.
[5] 李昌, 张见, 陈资南, 阚忠辉, 赵锐, 王晓军. 数据趋势融合分析方法在岩石破裂判识与预警中的应用[J]. 黄金科学技术, 2024, 32(3): 523-538.
[6] 刘宽, 莫冠旺, 李响, 沈平欢, 万波, 刘建坤. 超大断面扁平结构隧道施工参数优化研究[J]. 黄金科学技术, 2024, 32(2): 330-344.
[7] 王开彬, 刘钦, 王洪涛. 压力型锚索锚固段荷载传递特征及影响因素研究[J]. 黄金科学技术, 2024, 32(1): 123-131.
[8] 徐泽峰, 史秀志, 黄仁东, 丁文智, 陈新. 基于满管输送的充填管路优化研究[J]. 黄金科学技术, 2024, 32(1): 160-169.
[9] 李杰林, 刘一良, 王玉普, 李在利, 周科平, 程春龙. 高温独头巷道压抽混合式通风参数对人工制冷降温效果的影响[J]. 黄金科学技术, 2024, 32(1): 63-74.
[10] 费鸿禄, 纪海楠, 山杰. 露天台阶水介质间隔装药结构优选及对比试验研究[J]. 黄金科学技术, 2023, 31(6): 930-943.
[11] 单文法, 毛先成, 刘占坤, 邓浩, 陈进, 张维, 王海正, 杨鑫. 胶东大尹格庄金矿床成矿过程数值模拟及其找矿意义[J]. 黄金科学技术, 2023, 31(5): 707-720.
[12] 张玉, 王文己, 孙加奇, 肖永刚. 层理结构板岩动态断裂特性[J]. 黄金科学技术, 2023, 31(5): 803-810.
[13] 赵亚楠, 赵一航, 蒋中明, 赵红敏. 基于离散元法的高放核废料储罐静动力稳定性初步研究[J]. 黄金科学技术, 2023, 31(4): 592-604.
[14] 马恒,高嘉毅,李世虎,高科. 双机并联空气幕射流角度对巷道风流的影响[J]. 黄金科学技术, 2022, 30(5): 743-752.
[15] 郭对明,李国清,侯杰,胡乃联. 基于FLUENT的深井掘进巷道局部通风参数优化[J]. 黄金科学技术, 2022, 30(5): 753-763.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!

©2018 黄金科学技术编辑部

电话：0931-8277791　

E-mail: hjkx@lzb.ac.cn 邮编：730000　

陇ICP备17001318号-1    甘公网安备 62010202000672号
访问总数:    当日访问:    当前在线: