厚大破碎矿体二步矿房开采采场稳定性分析
|
摘要
山东黄金新城金矿深部V#矿体为一厚大破碎缓倾斜矿体,-580m中段矿块划分一步矿房和二步矿房开采,其中一步矿房宽8m,采用机械化盘区上向水平分层充填采矿法,开采即将结束。因矿体埋藏深,应力大,上盘为焦家断层带,二步矿房两侧为一步采场充填体,且一步充填水对二步矿岩侵蚀破坏,致使二步矿房矿岩异常破碎,二步矿房的开采安全性受到了严重的考验。经过分析研究,拟对二步矿房采用预控顶中深孔落矿嗣后充填采矿法开采。为了安全回采二步矿房,有必要对二步矿房进行一步采场充填体自立性效果和二步采场顶板稳固性等为主的采场稳定性影响因素进行分析。为此,本论文主要开展如下研究：
1.通过现场工程地质调查和岩石力学取样,得到深部矿体的岩体力学参数,并对岩体进行工程质量分级；
2.根据一步采场充填体受力模型,分析充填体内部应力,充填体自立所需的强度；
3.计算一步采场充填体自立高度,研究分析当二步矿房采用中深孔开采高度达到一个分段高时,一步矿房胶结充填体的自立性；
4.利用岩土有限元软件Midas/GTS对新城金矿二步矿房开采进行数值模拟计算,探讨按不同分段高度和不同采场开挖长度,对采场顶板稳定性的影响,并与试验采场监测进行比较验证。
The deep V# orebody of Xincheng Gold Mine in Shandong province is a gently inclined large broken orebody. In-580m level of the orebody extraction is divided into two steps mining. First-step room with 8 meters wide is mined by mechanized panel upward horizontal slice cut-and-fill method, because of deep burial depth, high crustal stress, Jiaojia Fault Zone in the roofs, and the orebody of second-step pillar is extremely broken and unstable. On the other hand filling water from the backfilling in first-step room deteriorates both sides of second-step pillars. The mining safety of second-step pillars mining is seriously challenged. In order to extract second-step stope safety and efficiently, mining method called medium-length hole blasting extraction by pre-controlling roof and first-step filling is put forward. In order to gurantee the safety of mining second-step pillars, it is essential to perform the stability analysis on free-standing of backfill in first-step rooms and the roof stability of second-step pillars. Major works in this paper are listed as follows:
1. Investigating mechanical parameters of rock mass at great depth are and classifying rock mass quality is carried out according to the obtained parameters.by field investigations of engineering geological conditions and rock mechanical sampling experiments.
2. Analysis internal stress and required strength of free-standing backfill based on the mechanical model of backfill in first-step rooms.
3. Calculating the height of free-standing backfill in first-step rooms related formulas. And discussing the free-standing cemented backfill in first-step rooms, as the height of medium-length hole blasting is designed as a sublevel height in current second-step pillars of Xincheng Gold Mine.
4. Numerically simulating according to the effect of different sublevel heights and excavation lengths on roof stability of second-step rooms by using geotechnical finite element software Midas/GTS. After that, comparing with field measurement results in the testing stope.
1.通过现场工程地质调查和岩石力学取样,得到深部矿体的岩体力学参数,并对岩体进行工程质量分级；
2.根据一步采场充填体受力模型,分析充填体内部应力,充填体自立所需的强度；
3.计算一步采场充填体自立高度,研究分析当二步矿房采用中深孔开采高度达到一个分段高时,一步矿房胶结充填体的自立性；
4.利用岩土有限元软件Midas/GTS对新城金矿二步矿房开采进行数值模拟计算,探讨按不同分段高度和不同采场开挖长度,对采场顶板稳定性的影响,并与试验采场监测进行比较验证。
The deep V# orebody of Xincheng Gold Mine in Shandong province is a gently inclined large broken orebody. In-580m level of the orebody extraction is divided into two steps mining. First-step room with 8 meters wide is mined by mechanized panel upward horizontal slice cut-and-fill method, because of deep burial depth, high crustal stress, Jiaojia Fault Zone in the roofs, and the orebody of second-step pillar is extremely broken and unstable. On the other hand filling water from the backfilling in first-step room deteriorates both sides of second-step pillars. The mining safety of second-step pillars mining is seriously challenged. In order to extract second-step stope safety and efficiently, mining method called medium-length hole blasting extraction by pre-controlling roof and first-step filling is put forward. In order to gurantee the safety of mining second-step pillars, it is essential to perform the stability analysis on free-standing of backfill in first-step rooms and the roof stability of second-step pillars. Major works in this paper are listed as follows:
1. Investigating mechanical parameters of rock mass at great depth are and classifying rock mass quality is carried out according to the obtained parameters.by field investigations of engineering geological conditions and rock mechanical sampling experiments.
2. Analysis internal stress and required strength of free-standing backfill based on the mechanical model of backfill in first-step rooms.
3. Calculating the height of free-standing backfill in first-step rooms related formulas. And discussing the free-standing cemented backfill in first-step rooms, as the height of medium-length hole blasting is designed as a sublevel height in current second-step pillars of Xincheng Gold Mine.
4. Numerically simulating according to the effect of different sublevel heights and excavation lengths on roof stability of second-step rooms by using geotechnical finite element software Midas/GTS. After that, comparing with field measurement results in the testing stope.
引文
[1]《金属矿山充填采矿法设计参考资料》编写组编.金属矿山充填采矿法设计参考资料[M].北京：冶金工业出版社,1978
[2]孙恒虎,黄玉诚,杨宝贵.当代胶结充填技术[M].北京：冶金工业出版社,2002
[3]古德生,李夕兵,等.现代金属矿床开采科学技术[M].北京：冶金工业出版社,2006
[4]姚银佩,李夕兵,赵国彦,等,新城金矿二步矿房采矿方法优选研究[J].矿业研究与开发,2010,30(5)：6～8,43
[5]李夕兵,赵国彦,李启月等,深部厚大破碎矿体安全高效低成本开采技术研究[R].长沙：中南大学,2010
[6]张绍国.广西高峰公司深部105号矿体开采稳定性及充填工艺与系统研究[硕士学位论文],长沙：中南大学,2005
[7]王新民,肖卫国,张钦礼.深井矿山充填理论与技术[M].长沙：中南大学出版社,2005
[8]孙朝芦.盘江矿区村下采煤技术研究[硕士学位论文],贵阳：贵州大学,2006
[9]邹连军.充填体对采场稳定性的影响[硕士学位论文].昆明理工大学,2003
[10]刘同有.充填采矿技术与应用[M].北京：冶金工业出版社,2001
[11]万海涛,方勇,肖广哲,秦杨.充填采矿法的应用现状及发展方向[J],世界有色金属,2009.(8)：26～28
[12]王永智.上向进路式尾砂胶结充填采矿法采场结构参数优化研究[硕士学位论文],昆明理工大学,2005
[13]彭续承.充填理论及应用.长沙：中南工业大学出版社,1998
[14]Stewart J M, Clark I H, Morris A N. Assessment of fill quality as basis for selecting and developing optimal backfill systems for South African gold mines. South African Inst of Mining&Metallurgy,1986, v:255-270
[15]Ryder JA.Application of numerical stress analyses to the design of deep mine. In:South African Inst of Mining&Metallurgy, eds. GOLD 100, Johannesburg: South African Inst of Mining&Metallurgy Publisher,1986.245-253
[16]杨根祥.焦家金矿全尾砂胶结充填体强度试验研究[J].沈阳黄金学院学报,1994,13(1)：32-37
[17]李庶林,桑玉发.尾砂胶结充填体的破坏机理及其损伤本构方程[J].黄金,1997,18(1)：24-29
[18]张世超,姚中亮.安庆铜矿特大型采场充填体稳定性分析[J],矿业研究与开发,2001,21(4)：12～15.
[19]石长岩.一步矿房胶结充填的实践及相关问题分析[J],采矿技术,2006,6(2)：7-10
[20]朱志彬,刘成平.充填体强度计算及稳定性分析[J],采矿技术,2008,8(3)：15～17
[21]邓建,李夕兵,彭怀生,确定稳定充填体临界高度的极限分析法[J],金属矿山,1999.7
[22]王劫.尾砂胶结充填体力学作用机理与应用[J].中国有色金属学报,1998,8(2)：797-801
[23]于学馥,邓颖人等,地下工程围岩稳定分析[M],煤炭工业出版社,
[24]刘志祥,李夕兵,爆破动载荷下高阶段充填体稳定性研究[J],矿冶工程,2004.24(3)：21～24
[25]刘可任.充填理论基础.北京：冶金工业出版社,1982
[26]山东黄金新城金矿,山东省莱州市新城金矿V号矿体详查地质报告[R],1994.11
[27]李启月,王树海,范作鹏,赵国彦.盘区阶梯式无间柱连续充填采矿法试验研究[J],矿冶工程,2010,30(3)：16～19
[28]崔栋梁,李夕兵,赵国彦.新城金矿难采矿体采场结构参数数值模拟分析[J].江西有色金属,2006,(3)：13～17
[29]刘滨,姬涛,由伟.新城金矿厚大矿体开采方法研究[J],采矿技术,2010,10(3)：31～32,69
[30]汪亦显.节理岩体力学参数测定及处理方法研究[硕士学位论文],长沙：中南大学,2007
[31]马紫娟.基于拉格朗日差分法的露天边坡稳定性研究[硕士学位论文],长沙：中南大学,2009
[32]伍佑伦,王平,钱锋.数值方法中岩体力学参数的选择[J].武钢技术,2002,40(4)：43～45
[33]李胡生,熊文林.岩体力学参数的工程模糊处理[J].水利学报,1994,06
[34]乔春生,张清,黄修云.岩石工程数值分析中选择岩体力学参数的神经元网络方法[J].岩石力学与工程学报,1999,01.
[35]黄修云.隧道工程岩石力学参数综合处理研究及多媒体实现[博士学位论文].北方交通大学.1998.09.
[36]蔡美峰,何满朝,刘东燕.岩石力学与工程[M].北京科学出版社,2002
[37]陈伊清,岩石质量指标R Q D的应用问题[J].水利科技,2002,01.
[38]杜时贵,许四法,杨树峰等.岩石质量指标RQD与工程岩体分类[J].工程地质学报,2000,8(3)：351～356
[39]王思敬,杨志法,刘竹华,地下工程岩体稳定分析[M],科学出版社,1984
[40]孙广忠.岩体结构力学[M].北京科学出版社,1998.
[41]谌文武,梁收运,刘高等.松动岩体的工程特性研究.岩石力学与工程学报,2000.11
[42]邓代强.安庆铜矿特大型采场充填体力学性能、损伤及稳定性研究[硕士学位论文],长沙矿山研究院,2005
[43]Thomas. E.G.,Nantel. J.H.and Notley. K.R., Fill Technology in Underground Metalliferous Mines, Interna, Acadamic Services Ltd., Canada,1979
[44]卢平.确定胶结充填体强度的理论与实践[J],黄金,1992,13(3)
[45]蔡嗣经.矿山充填材料的优化试验研究[J],江西有色金属.1987,1-6
[46]蔡嗣经.胶结充填材料的强度特性与强度设计(Ⅱ)—胶结充填体强度设计的几个理论模型[J],江西冶金学院学报,1985,(4)：12～21
[47]韩斌.金川二矿区充填体可靠度分析与Ⅰ#矿体回采地压控制优化研究[博士学位论文],长沙：中南大学,2004
[48]李夕兵,刘志祥.高阶段尾砂胶结充填体力学研究与博弈树配比优化,安全与环境学报,2004,4(4)：87～90
[49]薛奕忠.对安庆铜矿高阶段矿房尾胶体稳定性的探讨[J],有色金属：矿山部分,2001,4：11～15
[50]邓泽岱,等.矿柱强度估算及稳定性评价[J].武汉工业大学报,1993,15(3)：59～67
[51]骆玉骑.金川二矿区围岩与充填体稳定性分析[硕士学位论文].昆明理工大学,2005
[52]刘志祥.深部开采高阶段尾砂充填体力学与非线性优化设计[博士学位论文],长沙：中南大学,2005
[53]刘建庄.赵春景.浑宝炬.胶结充填体自立的受力分析与研究[J].有色矿冶,2008,24(4)：12～13,16
[54]余斌.影响尾砂胶结充填体强度的若干因素分析[J],河北冶金,2001,11(3)：3～6.
[55]罗建祥.影响胶结充填体强度的因素分析[J],甘肃冶金,2005,27(2)：24～25.
[56]杨宝贵,孙恒虎,庄百宏.高水固结充填体的自立[J].有色金属, 2000,52(2)：7-10
[57]姜凡均,王子庆.影响块石胶结充填体质量的主要因素[J].湖南有色金属,1999,15(5)：1-3.
[58]Application of Rock Mechanics to Cut-And-Fill Mining, Inst. Min. Met, London,1981
[59]Mining with Backfill, Proc Interna. Symp., Lulea, Sweden，1983
[60]邓代强,高永涛,姚中亮.胶结充填材料力学特性影响因素回归分析[J].有色金属,2008.60(4)：120～123,135.
[61]王振环,李政.安庆铜矿充填体特性现场调查研究[J].甘肃科技,2005,21(5)：85～87.
[62]郭利杰,杨小聪.深部采场胶结充填体力学稳定性研究[J].矿冶,2008,17(3)：10～13.
[63]夏永旭.工程结构数值分析方法[M].北京：人民交通出版社,2008
[64]许万忠.节理裂隙边坡稳定性及锚注加固效应研究[博士学位论文].长沙：中南大学,2006
[65]黎宏.时空效应的海底开采安全隔离厚度研究[硕士学位论文].长沙：中南大学,2009
[66]朱卓慧.公路下矿体开采保护地表的采矿方案优化研究[硕士学位论文].长沙：中南大学,2006
[67]刘浴宇,徐长佑.地下采空区矿柱稳定性分析[J],矿冶工程,2000.3
[68]苗立霞.黄金矿山尾砂充填空区的数值模拟研究[硕士学位论文].青岛理工大学,2007
[69]张强勇,李术才,焦玉勇.岩体数值分析方法与地质力学模型试验原理及工程应用[M].北京：中国水利水电出版社,2005
[70]杨慧.双轴压缩下含闭合裂纹有限岩板的断裂研究[硕士学位论文].长沙：中南大学,2007
[71]G.歌德赫.有限元法在岩土力学中的应用[M].张清,张弥译.北京：中国铁道出版社,1983
[72]李围.隧道及地下工程FLAC解析方法[M].北京：中国水利水电出版社,2009
[73]姚寿广等.边界元数值方法及其工程应用[M].北京：国防工业出版社,1995
[74]王泳嘉.边界元法在岩石力学中的应用[J].岩石力学与工程学报,1986,5(2)：205～222
[75]王泳嘉,邢纪波.离散单元法及其在岩土力学中的应用[M].辽宁：东北工学院出版社,1991
[76]王泳嘉.离散元法及其在岩石力学中的应用[J].金属矿山,1986,8
[77]刘锦华,吕祖巧.块体理论在工程岩体稳定分析中的应用.北京：水利电力出版社,1988
[78]朱伯芳.有限单元法原理与应用(第二版)[M].北京：中国水利水电出版社,1998
[79]王勖成,邵敏.有限单元法基本原理和数值方法(第2版)[M].北京：清华大学出版社,2004
[80]北京迈达斯技术有限公司.midas GTS使用手册[M].2010
[81]张保.深井开采矿柱稳定性分析与可视化验证[硕士学位论文].长沙：中南大学,2009
[82]帅红岩,韩文喜,赵晋乾.MidasGTS软件在边坡三维稳定分析中的应用[J].地质灾害与环境保护,2009,20(3)：104～107
[83]孙钧,汪炳槛.地下结构有限元法解析[M],上海：同济大学出版社1988
[84]谢康和,周健.岩土工程有限元分析理论与应用[M].北京：科学出版社,2002
[85]李爱兵,周先明.安庆铜矿高阶段回采充填体-矿体-岩体稳定性的有限元分析[J].矿业研究与开发,2000,20(1)：19～21.
[86]薛奕忠.对安庆铜矿高阶段矿房尾胶体稳定性的探讨[J],有色金属：矿山部分,2001,4：11～15
[2]孙恒虎,黄玉诚,杨宝贵.当代胶结充填技术[M].北京：冶金工业出版社,2002
[3]古德生,李夕兵,等.现代金属矿床开采科学技术[M].北京：冶金工业出版社,2006
[4]姚银佩,李夕兵,赵国彦,等,新城金矿二步矿房采矿方法优选研究[J].矿业研究与开发,2010,30(5)：6～8,43
[5]李夕兵,赵国彦,李启月等,深部厚大破碎矿体安全高效低成本开采技术研究[R].长沙：中南大学,2010
[6]张绍国.广西高峰公司深部105号矿体开采稳定性及充填工艺与系统研究[硕士学位论文],长沙：中南大学,2005
[7]王新民,肖卫国,张钦礼.深井矿山充填理论与技术[M].长沙：中南大学出版社,2005
[8]孙朝芦.盘江矿区村下采煤技术研究[硕士学位论文],贵阳：贵州大学,2006
[9]邹连军.充填体对采场稳定性的影响[硕士学位论文].昆明理工大学,2003
[10]刘同有.充填采矿技术与应用[M].北京：冶金工业出版社,2001
[11]万海涛,方勇,肖广哲,秦杨.充填采矿法的应用现状及发展方向[J],世界有色金属,2009.(8)：26～28
[12]王永智.上向进路式尾砂胶结充填采矿法采场结构参数优化研究[硕士学位论文],昆明理工大学,2005
[13]彭续承.充填理论及应用.长沙：中南工业大学出版社,1998
[14]Stewart J M, Clark I H, Morris A N. Assessment of fill quality as basis for selecting and developing optimal backfill systems for South African gold mines. South African Inst of Mining&Metallurgy,1986, v:255-270
[15]Ryder JA.Application of numerical stress analyses to the design of deep mine. In:South African Inst of Mining&Metallurgy, eds. GOLD 100, Johannesburg: South African Inst of Mining&Metallurgy Publisher,1986.245-253
[16]杨根祥.焦家金矿全尾砂胶结充填体强度试验研究[J].沈阳黄金学院学报,1994,13(1)：32-37
[17]李庶林,桑玉发.尾砂胶结充填体的破坏机理及其损伤本构方程[J].黄金,1997,18(1)：24-29
[18]张世超,姚中亮.安庆铜矿特大型采场充填体稳定性分析[J],矿业研究与开发,2001,21(4)：12～15.
[19]石长岩.一步矿房胶结充填的实践及相关问题分析[J],采矿技术,2006,6(2)：7-10
[20]朱志彬,刘成平.充填体强度计算及稳定性分析[J],采矿技术,2008,8(3)：15～17
[21]邓建,李夕兵,彭怀生,确定稳定充填体临界高度的极限分析法[J],金属矿山,1999.7
[22]王劫.尾砂胶结充填体力学作用机理与应用[J].中国有色金属学报,1998,8(2)：797-801
[23]于学馥,邓颖人等,地下工程围岩稳定分析[M],煤炭工业出版社,
[24]刘志祥,李夕兵,爆破动载荷下高阶段充填体稳定性研究[J],矿冶工程,2004.24(3)：21～24
[25]刘可任.充填理论基础.北京：冶金工业出版社,1982
[26]山东黄金新城金矿,山东省莱州市新城金矿V号矿体详查地质报告[R],1994.11
[27]李启月,王树海,范作鹏,赵国彦.盘区阶梯式无间柱连续充填采矿法试验研究[J],矿冶工程,2010,30(3)：16～19
[28]崔栋梁,李夕兵,赵国彦.新城金矿难采矿体采场结构参数数值模拟分析[J].江西有色金属,2006,(3)：13～17
[29]刘滨,姬涛,由伟.新城金矿厚大矿体开采方法研究[J],采矿技术,2010,10(3)：31～32,69
[30]汪亦显.节理岩体力学参数测定及处理方法研究[硕士学位论文],长沙：中南大学,2007
[31]马紫娟.基于拉格朗日差分法的露天边坡稳定性研究[硕士学位论文],长沙：中南大学,2009
[32]伍佑伦,王平,钱锋.数值方法中岩体力学参数的选择[J].武钢技术,2002,40(4)：43～45
[33]李胡生,熊文林.岩体力学参数的工程模糊处理[J].水利学报,1994,06
[34]乔春生,张清,黄修云.岩石工程数值分析中选择岩体力学参数的神经元网络方法[J].岩石力学与工程学报,1999,01.
[35]黄修云.隧道工程岩石力学参数综合处理研究及多媒体实现[博士学位论文].北方交通大学.1998.09.
[36]蔡美峰,何满朝,刘东燕.岩石力学与工程[M].北京科学出版社,2002
[37]陈伊清,岩石质量指标R Q D的应用问题[J].水利科技,2002,01.
[38]杜时贵,许四法,杨树峰等.岩石质量指标RQD与工程岩体分类[J].工程地质学报,2000,8(3)：351～356
[39]王思敬,杨志法,刘竹华,地下工程岩体稳定分析[M],科学出版社,1984
[40]孙广忠.岩体结构力学[M].北京科学出版社,1998.
[41]谌文武,梁收运,刘高等.松动岩体的工程特性研究.岩石力学与工程学报,2000.11
[42]邓代强.安庆铜矿特大型采场充填体力学性能、损伤及稳定性研究[硕士学位论文],长沙矿山研究院,2005
[43]Thomas. E.G.,Nantel. J.H.and Notley. K.R., Fill Technology in Underground Metalliferous Mines, Interna, Acadamic Services Ltd., Canada,1979
[44]卢平.确定胶结充填体强度的理论与实践[J],黄金,1992,13(3)
[45]蔡嗣经.矿山充填材料的优化试验研究[J],江西有色金属.1987,1-6
[46]蔡嗣经.胶结充填材料的强度特性与强度设计(Ⅱ)—胶结充填体强度设计的几个理论模型[J],江西冶金学院学报,1985,(4)：12～21
[47]韩斌.金川二矿区充填体可靠度分析与Ⅰ#矿体回采地压控制优化研究[博士学位论文],长沙：中南大学,2004
[48]李夕兵,刘志祥.高阶段尾砂胶结充填体力学研究与博弈树配比优化,安全与环境学报,2004,4(4)：87～90
[49]薛奕忠.对安庆铜矿高阶段矿房尾胶体稳定性的探讨[J],有色金属：矿山部分,2001,4：11～15
[50]邓泽岱,等.矿柱强度估算及稳定性评价[J].武汉工业大学报,1993,15(3)：59～67
[51]骆玉骑.金川二矿区围岩与充填体稳定性分析[硕士学位论文].昆明理工大学,2005
[52]刘志祥.深部开采高阶段尾砂充填体力学与非线性优化设计[博士学位论文],长沙：中南大学,2005
[53]刘建庄.赵春景.浑宝炬.胶结充填体自立的受力分析与研究[J].有色矿冶,2008,24(4)：12～13,16
[54]余斌.影响尾砂胶结充填体强度的若干因素分析[J],河北冶金,2001,11(3)：3～6.
[55]罗建祥.影响胶结充填体强度的因素分析[J],甘肃冶金,2005,27(2)：24～25.
[56]杨宝贵,孙恒虎,庄百宏.高水固结充填体的自立[J].有色金属, 2000,52(2)：7-10
[57]姜凡均,王子庆.影响块石胶结充填体质量的主要因素[J].湖南有色金属,1999,15(5)：1-3.
[58]Application of Rock Mechanics to Cut-And-Fill Mining, Inst. Min. Met, London,1981
[59]Mining with Backfill, Proc Interna. Symp., Lulea, Sweden，1983
[60]邓代强,高永涛,姚中亮.胶结充填材料力学特性影响因素回归分析[J].有色金属,2008.60(4)：120～123,135.
[61]王振环,李政.安庆铜矿充填体特性现场调查研究[J].甘肃科技,2005,21(5)：85～87.
[62]郭利杰,杨小聪.深部采场胶结充填体力学稳定性研究[J].矿冶,2008,17(3)：10～13.
[63]夏永旭.工程结构数值分析方法[M].北京：人民交通出版社,2008
[64]许万忠.节理裂隙边坡稳定性及锚注加固效应研究[博士学位论文].长沙：中南大学,2006
[65]黎宏.时空效应的海底开采安全隔离厚度研究[硕士学位论文].长沙：中南大学,2009
[66]朱卓慧.公路下矿体开采保护地表的采矿方案优化研究[硕士学位论文].长沙：中南大学,2006
[67]刘浴宇,徐长佑.地下采空区矿柱稳定性分析[J],矿冶工程,2000.3
[68]苗立霞.黄金矿山尾砂充填空区的数值模拟研究[硕士学位论文].青岛理工大学,2007
[69]张强勇,李术才,焦玉勇.岩体数值分析方法与地质力学模型试验原理及工程应用[M].北京：中国水利水电出版社,2005
[70]杨慧.双轴压缩下含闭合裂纹有限岩板的断裂研究[硕士学位论文].长沙：中南大学,2007
[71]G.歌德赫.有限元法在岩土力学中的应用[M].张清,张弥译.北京：中国铁道出版社,1983
[72]李围.隧道及地下工程FLAC解析方法[M].北京：中国水利水电出版社,2009
[73]姚寿广等.边界元数值方法及其工程应用[M].北京：国防工业出版社,1995
[74]王泳嘉.边界元法在岩石力学中的应用[J].岩石力学与工程学报,1986,5(2)：205～222
[75]王泳嘉,邢纪波.离散单元法及其在岩土力学中的应用[M].辽宁：东北工学院出版社,1991
[76]王泳嘉.离散元法及其在岩石力学中的应用[J].金属矿山,1986,8
[77]刘锦华,吕祖巧.块体理论在工程岩体稳定分析中的应用.北京：水利电力出版社,1988
[78]朱伯芳.有限单元法原理与应用(第二版)[M].北京：中国水利水电出版社,1998
[79]王勖成,邵敏.有限单元法基本原理和数值方法(第2版)[M].北京：清华大学出版社,2004
[80]北京迈达斯技术有限公司.midas GTS使用手册[M].2010
[81]张保.深井开采矿柱稳定性分析与可视化验证[硕士学位论文].长沙：中南大学,2009
[82]帅红岩,韩文喜,赵晋乾.MidasGTS软件在边坡三维稳定分析中的应用[J].地质灾害与环境保护,2009,20(3)：104～107
[83]孙钧,汪炳槛.地下结构有限元法解析[M],上海：同济大学出版社1988
[84]谢康和,周健.岩土工程有限元分析理论与应用[M].北京：科学出版社,2002
[85]李爱兵,周先明.安庆铜矿高阶段回采充填体-矿体-岩体稳定性的有限元分析[J].矿业研究与开发,2000,20(1)：19～21.
[86]薛奕忠.对安庆铜矿高阶段矿房尾胶体稳定性的探讨[J],有色金属：矿山部分,2001,4：11～15