img

Wechat

Adv. Search

Gold Science and Technology ›› 2019, Vol. 27 ›› Issue (2): 189-198.doi: 10.11872/j.issn.1005-2518.2019.02.189

• Mining Technology and Mine Management • Previous Articles    

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

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

CLC Number: 

  • TD853.36

Fig.1

Mine topography and terrain contour"

Fig.2

PFC2D mine model"

Table 1

Statistical results of dominant joint set occurrence"

组号占比/%倾向/(°)倾角/(°)
最大值最小值均值方差分布规律最大值最小值均值方差分布规律
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正态分布

Table 2

Statistical results of joint spacing and continuity(m)"

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

Table 3

Parameters of DFN generated"

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

Fig.3

Local enlarged diagram of discrete joints generated with DFN"

Fig.4

Statistical chart of DFN spacing"

Fig.5

Plots of uniaxial compression test"

Fig.6

Plots of triaxial servo compression test"

Table 4

Comparison of actual physical mechanics parameters of rock with physical and mechanical parameters of rock simulation"

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

Table 5

Microscopic parameters used in the sumulation"

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

Fig.7

Flow chart of joint parameters calibration"

Table 6

Parameters of joint"

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

Fig.8

Crack propagation in the first step"

Fig.9

Crack propagation and statistics of fracture types"

Fig.10

Evolution characteristics of caving area shape"

Fig.11

Statistics of caving area"

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] Xiang LI,Zhen HUAI,Xibing LI,Zhuoyao ZHANG. Study on Fracture Characteristics and Mechanical Properties of Brittle Rock Based on Crack Propagation Model [J]. Gold Science and Technology, 2019, 27(1): 41-51.
[2] Chong CHEN,Xibing LI,Fan FENG. Numerical Study on Damage Zones of the Induced Roadway Surrounding Rock [J]. Gold Science and Technology, 2018, 26(6): 771-779.
[3] Liansheng LIU,Qingliang ZHONG,Lei YAN,Wen ZHONG,Longhua LIANG. Measurement of Sound Waves to Study Cumulative Damage Effect on Saturated Rock Under Frequent Production Blast Loading [J]. Gold Science and Technology, 2018, 26(6): 750-760.
[4] Bo PENG. Numerical Simulation Study of Seepage Stability of Mudstone Slope with Dif-ferent Rainfall Types Coupled with Sudden Drop of Reservoir Water Level [J]. Gold Science and Technology, 2018, 26(5): 647-655.
[5] Guang LI,Fengshan MA,Gang LIU,Jie GUO,Huigao GUO,Yongyuan KOU. Study on Supporting Parametric Optimizing Design and Evaluate Supporting Effect of Deep Roadway in Jinchuan Mine [J]. Gold Science and Technology, 2018, 26(5): 605-614.
[6] Jianhua HU,Qifan REN,Zhonghua QI,Jiwei ZHANG. Regress Optimize Model of Limit Exposure Area to Stope in Wohushan Iron Mine [J]. Gold Science and Technology, 2018, 26(4): 503-510.
[7] HU Guiying, LIU Kewei, DU Xin, LI Xiaohan. Research on Smooth-cutting Method and Its Application in Tunnel Excavation [J]. Gold Science and Technology, 2018, 26(3): 349-356.
[8] WANG Haibo, WEI Guoli, ZONG Qi, XU Ying. Numerical Simulation and Application Research on Joint Development Rock Roadway Blasting Excavation [J]. Gold Science and Technology, 2018, 26(3): 342-348.
[9] YANG Mingcai, SHENG Jianlong, YE Zuyang, DONG Shu. Analysis of Sensitivity Factors of Open-pit Mine Slope Stability and Impact Based on FlAC3D [J]. Gold Science and Technology, 2018, 26(2): 179-186.
[10] CHENG Hao, XU Tao, ZHOU Guanglei, FANG Ke. Numerical Simulation of Fracture and Acoustic Emission Evolution of Hetero-geneous Rocks [J]. Gold Science and Technology, 2018, 26(2): 170-178.
[11] REN Weicheng,QIAO Dengpan,ZHOU Zhiwei. Research on the Change Law of Volume Fraction of Tailings with Tailings Silo Height in Vertical Tailings Silo [J]. Gold Science and Technology, 2018, 26(1): 64-73.
[12] MA Fengshan,LI Kepeng,DU Yunlong,HOU Chenglu,LI Wei2,ZHANG Guodong. Analysis on the Possible Failure Modes of Water Burst Prevention Structures of F1 Fault Caused by Undersea Mining in Sanshandao Gold Mine [J]. Gold Science and Technology, 2017, 25(5): 47-56.
[13] LI Guang,MA Fengshan,GUO Jie,LU Rong,CHEN Zhongjie,KOU Yongyuan. Application of Steel Tube Confined Concrete Support in Jinchuan Mine [J]. Gold Science and Technology, 2017, 25(5): 57-66.
[14] BAI Zhaoyang,WANG Guowei,ZHANG Peng,LIU Shuanping. Research on Pillaring Blasting Location Under Level Gob Group Conditions [J]. Gold Science and Technology, 2017, 25(4): 81-86.
[15] HU Jianhua,YANG Chun,ZHOU Bingren,ZHOU Keping,ZHANG Shaoguo. Simulation of Fracture Propagation and Optimization of Parameters for Smooth Blasting of Coping in Roadway [J]. Gold Science and Technology, 2017, 25(2): 45-53.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!