深井矿山盘区隔离矿柱回采方案及数值模拟研究
|
摘要
长期以来,在岩体稳定性问题的研究中多是依靠经验、半经验的方法分析,一些潜在的不安全因素难以得到揭示。随着工程规模和复杂性的日益增大,这种依靠经验的方法越来越不能适应发展的工程实际需要,难以有效解决岩体工程中的稳定性问题。在岩土和地下工程中,经常遇到的岩石力学问题是需要确定开挖前后的应力和位移及其在开挖过程中的变化,以便为设计和施工提供可靠的定量依据。然而已有的理论往往只能解决圆形或椭圆形等简单孔周应力场和位移场,对于几何形状复杂的孔洞如采场问题只能用数值模拟方法来近似求解。随着计算分析软件的发展数值模拟方法有了显著的进步,已成为岩石力学研究和工程计算的重要手段。
冬瓜山铜矿第三步骤回采隔离矿柱的时候,由于经过前面大规模的第一步骤矿房回采和第二步骤矿柱回采,采矿技术条件已经发生了更本性的变化。首先,隔离矿柱周围的介质由强度很好的矿岩体基本变成了强度较低的尾胶充填体或基本没有强度的尾砂充填体；其次,由于受地压转移效应的影响,隔离矿柱内的应力会明显增大,矿岩体次生节理裂隙可能大量发育(矿岩体原生节理均不发育,平均1条/m)；总之,到第三步骤回采隔离矿柱的时候,采矿技术条件比第一步骤回采矿房和第二步骤回采矿柱的时候严重恶化。
本文针对冬瓜山铜矿盘区隔离矿柱回采时围岩力学性质的复杂性,主要研究内容为受大规模采动影响后隔离矿柱开采技术条件的研究,以及隔离矿柱采场合理结构尺寸参数的数值模拟研究。具体内容包括以下几点：
(1)对深井开采岩体破坏机理进行研究,明确岩体破坏的基本方式,然后对岩体破坏判据进行讨论和分析。最后,通过对岩体稳定性的分析,总结出综合考虑极限应力和应变条件,来对地下工程岩体稳定性进行综合评判。
(2)分析冬瓜山铜矿暂留隔离矿柱阶段空场嗣后充填采矿方法的技术与经济条件,分析隔离矿柱回采技术与工艺的复杂性。
(3)通过分析冬瓜山铜矿隔离矿柱回采的各项技术条件,选取影响隔离矿柱回采稳定性的6项因素,按照正交试验的原理,安排8种不同的数值模拟试验方案。
(4)应用FLAC3D软件数值模拟,分析隔离矿柱回采前的开采技术条件,完成对8种方案的模拟试验,并比较各方案条件下的最大应力应变数值和规律,进而分析围岩稳定性状态。然后通过正交分析,得出6项因素对隔离矿柱回采稳定性影响的权重。
(5)根据数值模拟结果正交分析得出的各因素影响权重,结合冬瓜山铜矿开采的技术与经济条件,建议一个隔离矿柱回采结构参数较优的方案。
Over the years, it mainly relies on the method of multiple, semi-empirical experience to analyze the rock stability problems, however, some potential insecure elements are difficult to be revealed. This method can not meet the needs of the actual engineering with the increasing project size and complexity, it can not effectively solve the stability problems in rock engineering. The need to determine the stress and strain around the excavation and changes in the excavation process is often encountered in rock mechanics problems in geotechnical and underground engineering, it can provide reliable quantitative basis for the design and construction. However, the existed theories can only solve such a simple round or oval holes on stress and strain fields, the complex geometry of the holes, such as mining field problems, can only be approximated by numerical simulation solution. With the development of analysis software, numerical simulation method has made substantial progress, and has become an important means in rock mechanics research and engineering calculations.
Before the excavation of the third step, the barrier pillar in Dongguashan Copper Mine, we have excavated a large scale of the first steps and the second steps, so the mining technology conditions have essentially changed. First, the medium around the barrier pillar are turned into tailings backfill, which have low strength or little plastic strength. Secondly, due to the affect of the pressure transfer, the stress within the barrier pillar have been significantly increased, the mine fractured rock mass may be a lot of secondary development. In short, when excavating the third steps, the technical conditions of mining is more seriously deteriorated than it in the first steps and the second steps.
The main research content of this paper are the mining technical conditions of the barrier pillar after the large-scale mining effect, and the numerical simulation about the reasonable structure size of the barrier pillar. Contents include the following:
(1) Ascertaining the basic form of mass failure after the research of rock failure mechanism in deep mine. Then discussing and analyzing the rock mass failure.
(2) Analyzing the technical and economic conditions of the stay barrier pillar level room subsequent filling mining method, analyzing the technical complexity of the excavation of barrier pillar.
(3) Selecting six factors of affecting the barrier pillar stability, in accordance with the principle of orthogonal design, arranging 8 kinds of numerical simulation test program.
(4) Analyzing the barrier pillar's mining technical condition by applying FLAC3D numerical simulation software. Complete the 8 kinds of numerical simulation test program, and compare the maximum stress values and discipline under the various programs, and then analyzes the stability of surrounding rock state. At last, obtain the weight of the six factors on the barrier pillar's stability by orthogonal analysis.
(5) According to the weight of 6 factors obtained in the orthogonal numerical simulation, and combined with the technical and economic conditions in Dongguashan Copper Mine, propose an optimum scheme of the structure parameters in the excavation of barrier pillar.
冬瓜山铜矿第三步骤回采隔离矿柱的时候,由于经过前面大规模的第一步骤矿房回采和第二步骤矿柱回采,采矿技术条件已经发生了更本性的变化。首先,隔离矿柱周围的介质由强度很好的矿岩体基本变成了强度较低的尾胶充填体或基本没有强度的尾砂充填体；其次,由于受地压转移效应的影响,隔离矿柱内的应力会明显增大,矿岩体次生节理裂隙可能大量发育(矿岩体原生节理均不发育,平均1条/m)；总之,到第三步骤回采隔离矿柱的时候,采矿技术条件比第一步骤回采矿房和第二步骤回采矿柱的时候严重恶化。
本文针对冬瓜山铜矿盘区隔离矿柱回采时围岩力学性质的复杂性,主要研究内容为受大规模采动影响后隔离矿柱开采技术条件的研究,以及隔离矿柱采场合理结构尺寸参数的数值模拟研究。具体内容包括以下几点：
(1)对深井开采岩体破坏机理进行研究,明确岩体破坏的基本方式,然后对岩体破坏判据进行讨论和分析。最后,通过对岩体稳定性的分析,总结出综合考虑极限应力和应变条件,来对地下工程岩体稳定性进行综合评判。
(2)分析冬瓜山铜矿暂留隔离矿柱阶段空场嗣后充填采矿方法的技术与经济条件,分析隔离矿柱回采技术与工艺的复杂性。
(3)通过分析冬瓜山铜矿隔离矿柱回采的各项技术条件,选取影响隔离矿柱回采稳定性的6项因素,按照正交试验的原理,安排8种不同的数值模拟试验方案。
(4)应用FLAC3D软件数值模拟,分析隔离矿柱回采前的开采技术条件,完成对8种方案的模拟试验,并比较各方案条件下的最大应力应变数值和规律,进而分析围岩稳定性状态。然后通过正交分析,得出6项因素对隔离矿柱回采稳定性影响的权重。
(5)根据数值模拟结果正交分析得出的各因素影响权重,结合冬瓜山铜矿开采的技术与经济条件,建议一个隔离矿柱回采结构参数较优的方案。
Over the years, it mainly relies on the method of multiple, semi-empirical experience to analyze the rock stability problems, however, some potential insecure elements are difficult to be revealed. This method can not meet the needs of the actual engineering with the increasing project size and complexity, it can not effectively solve the stability problems in rock engineering. The need to determine the stress and strain around the excavation and changes in the excavation process is often encountered in rock mechanics problems in geotechnical and underground engineering, it can provide reliable quantitative basis for the design and construction. However, the existed theories can only solve such a simple round or oval holes on stress and strain fields, the complex geometry of the holes, such as mining field problems, can only be approximated by numerical simulation solution. With the development of analysis software, numerical simulation method has made substantial progress, and has become an important means in rock mechanics research and engineering calculations.
Before the excavation of the third step, the barrier pillar in Dongguashan Copper Mine, we have excavated a large scale of the first steps and the second steps, so the mining technology conditions have essentially changed. First, the medium around the barrier pillar are turned into tailings backfill, which have low strength or little plastic strength. Secondly, due to the affect of the pressure transfer, the stress within the barrier pillar have been significantly increased, the mine fractured rock mass may be a lot of secondary development. In short, when excavating the third steps, the technical conditions of mining is more seriously deteriorated than it in the first steps and the second steps.
The main research content of this paper are the mining technical conditions of the barrier pillar after the large-scale mining effect, and the numerical simulation about the reasonable structure size of the barrier pillar. Contents include the following:
(1) Ascertaining the basic form of mass failure after the research of rock failure mechanism in deep mine. Then discussing and analyzing the rock mass failure.
(2) Analyzing the technical and economic conditions of the stay barrier pillar level room subsequent filling mining method, analyzing the technical complexity of the excavation of barrier pillar.
(3) Selecting six factors of affecting the barrier pillar stability, in accordance with the principle of orthogonal design, arranging 8 kinds of numerical simulation test program.
(4) Analyzing the barrier pillar's mining technical condition by applying FLAC3D numerical simulation software. Complete the 8 kinds of numerical simulation test program, and compare the maximum stress values and discipline under the various programs, and then analyzes the stability of surrounding rock state. At last, obtain the weight of the six factors on the barrier pillar's stability by orthogonal analysis.
(5) According to the weight of 6 factors obtained in the orthogonal numerical simulation, and combined with the technical and economic conditions in Dongguashan Copper Mine, propose an optimum scheme of the structure parameters in the excavation of barrier pillar.
引文
[1]古德生,李夕兵等.现代金属矿床开采科学技术[M].北京：冶金工业出版社,2006
[2]何满潮,谢和平,彭苏萍等.深部开采岩体力学研究[J].岩石力学与工程学报,2005,24(16)：2803-2813
[3]孟稳权.冬瓜山矿床深部开采岩层控制的措施[J].矿业快报,2006,452(12)： 58-59
[4]杨承祥.深井金属矿床高效开采及地压监控技术研究[博士学位论文].长沙：中南大学,2007
[5]朱维申,何满潮.复杂条件下围岩稳定性与岩体动态施工力学[M].北京：科学出版社,1996
[6]李华晔.地下洞室围岩稳定性分析[M].北京：中国水利电力出版社,1999
[7]许传华,任青文,李瑞.地下工程围岩稳定性分析方法研究进展[J].金属矿山,2003,320(2)：34-37
[8]邵国建,卓家寿,章青.岩体稳定性分析与评判准则研究[J].岩石力学与工程学报,2003,22(5)：691-696
[9]王生全,夏玉成,樊怀仁.邢台显德汪井田煤层顶板稳定性评价与预测[J].西北地质,1997,18(4)：46-50
[10]杨仕教,丁德馨,古德生等.丰山铜矿采场稳定性的FUZZY综合评判[J].有色金属,2002,54(3)：99-102
[11]郭立,吴爱祥,周克平等.基于RES理论的岩体失稳模式判别及其智能实现[J].中国有色金属学报,2003,13(3)：749-753
[12]Hudson J A. Rock Engineering Systems:Theory and Practice[M]. Horwood Chichester,1992
[13]张保.深井开采矿柱稳定性分析与可视化验证[硕士学位论文].长沙：中南大学,2009
[14]宫凤强,李夕兵,高科.地下工程围岩稳定性分类的突变级数法研究[J].中南大学学报(自然科学版),2008,39(5)：1081-1086
[15]蔡美峰,孔广亚,贾立宏.岩体工程系统失稳的能量突变判断准则及其应用[J].北京科技大学学报,1997,19(4)：325-328
[16]朱永全.隧道稳定性位移判别准则[J].中国铁道科学,2001,22(6)： 80-83
[17]唐礼忠,潘长良,谢学斌等.冬瓜山铜矿深井开采岩爆危险区分析与预测[J].中南工业大学学报,2002,33(4)：335-338
[18]解联库,杨小聪,刘庆林等.基于霍克-布朗破坏准则的采空区稳定性分析[J].有色金属,2008,60(3)：109-111
[19]高峰,周科平,胡建华.采场稳定性的模糊物元评价模型及应用研究[J].采矿与安全工程学报,2006,23(2)：164-168
[20]李增学,吕大炜,王怀洪等.龙固井田主采煤层顶板稳定型特征研究[J].岩土工程学报,2008,30(4)：473-479
[21]冯玉国.灰色优化理论模型在地下工程围岩稳定性分类中的应用[J].岩土工程学报,1996,18(3)：62-66
[22]胡全舟,吴超,陈沅江.风险评价方法在地下工程围岩稳定性分级中的应用[J].地下空间与工程学报,2005,1(6)：874-877
[23]曹刚.康家湾矿深部难采矿体采场稳定性及安全开采技术研究[硕士学位论文].中南大学,2008
[24]刘培慧.基于应力边界法厚大矿体采场结构参数数值模拟优化研究[硕士学位论文].中南大学,2009
[25]布朗.工程岩体力学中的解析与数值计算方法(佘诗刚,王可均译)[M].北京：科学技术出版社,1997
[26]Chandrakant S.Dcsai,John T.Chritian.岩土工程数值方法(潘善德译)[M].北京：科学技术出版社
[27]姚振汉.边界元法[M].北京：高等教育出版社,2010
[28]Peter Hunter,Andrew Pul lan. FEM/BEM NOTES.Department of Engineering Science,The University of Auckland,New Zealand,2003
[29]江见鲸,何放龙,何益斌等.有限元法及其应用[M].北京：机械工业出版社,2006
[30]王泳嘉,邢纪波.离散元法及其在岩土力学中的应用[M].沈阳：东北大学出版社,1993
[31]孙广忠.岩体力学基础[M].北京：科学出版社,1983
[32]陈子光.岩石力学性质与构造应力场[M].北京：地质出版社,1986
[33]朱德仁.岩石工程破坏准则[J].煤炭学报,1994,19(1)：15-19
[34]王文星.岩体力学[M].长沙：中南大学出版社,2004
[35]吴爱祥,郭立,张卫锋.深井开采岩体破坏机理及工程控制方法综述[J].矿业研究与开发,2001,21(2)：4-7
[36]张孟喜,邓荣基,朱文喜.岩石破坏机理与岩石强度准则新探[J].甘肃科学学报,1996,8(1)56-61
[37]俞茂宏,昝月稳,范文等.20世纪岩石强度理论的发展-纪念Mohr-Coulomb强度理论100周年[J].岩石力学与工程学报,2000,19(5)：545-549
[38]Wittke W. Rock Mechanics:Theory and Applications with Case Histories [M]. Berlin:Springer-Verlag,1990,中译本：岩石力学(曾国熙译).杭州：浙江大学出版,1996
[39]Yu M H, Zhao J, Guan L W. Strength theory for rock and concrete: history, present situation and development [J]. Progress in Natural Science (English Edition),1998,8 (4):394-402
[40]白国良,梁冰.巷道岩体能量破坏机理探讨[J].矿山压力与顶板管理,2004,(1)：86-88
[41]杨官涛.地下采场结构参数优化及稳定性的能量突变分析[硕士学文论文].长沙：中南大学,2007
[42]付成华,陈胜宏.基于突变理论的地下工程洞室围岩失稳判据研究[J].岩土力学,2008,29(1)：167-172
[43]李忠华,官福海,潘一山.基于损伤理论的圆形巷道围岩应力场分析[J].岩土力学,2004,25(增刊)：160-163
[44]马青.基于损伤理论的圆形水工压力隧洞应力场分析[硕士学位论文].西安：长安大学,2009
[45]黄梦宏.岩石损伤演化模型和损伤临界值理论与实验研究[硕士学位论文].长沙：中南大学,2003
[46]凌建明,孙钧.建立在损伤应变空间的岩体破坏准则[J].同济大学学报,1995,23(5)：483-487
[47]徐卫亚,韦立德.岩石损伤统计本构模型的研究[J].岩石力学与工程学报,2002,21(6)：787-791
[48]田运生,张舒畅,李战军等.岩巷爆破威严损伤演化分析及损伤判据的建立[J].黄金,2005,26(8)：22-24
[49]李树茂,齐伟,刘红帅.岩体损伤力学理论进展[J].世界地质,2001,20(1)：72-78
[50]仵彦卿,丁卫华,蒲毅彬等.压缩条件下岩石密度损伤增量的CT动态观测[J].自然科学进展,2000,10(9)：830-835
[51]许传华.岩体破坏的非线性理论研究及应用[博士学位论文].南京：河海大学,2004
[52]孔广亚,蔡美峰,贾立宏.岩体工程系统破坏失稳的能量突变准则[J].有色矿冶,1996,(3)：1-3
[53]许传华,任青文,李瑞.围岩稳定的熵突变理论研究[J].岩石力学与工程学报,2004,23(12)：1992-1995
[54]杨官涛,李夕兵,王其胜等.地下采场失稳的能量突变判断准则及其应用[J].采矿与安全工程学报,2008,25(3)：268-276
[55]朱和玲.基于开采环境再造人工结构的稳定性与可开度研究[硕士学位论文].长沙：中南大学,2008
[56]北京有色冶金设计研究总院.铜陵有色金属(集团)公司冬瓜山铜矿可行性研究报告[R].1998
[57]北京有色冶金设计研究总院.安徽铜都铜业股份有限公司狮子山铜矿千米深井300万t级矿山开采条件研究成果[R].2001
[58]李冬青,王李管,等.深井硬岩大规模开采理论与技术——冬瓜山铜矿床开采研究与实践[M].北京：冶金工业出版社,2009
[59]北京有色冶金设计研究总院.狮子山铜矿冬瓜山矿床岩石力学试验报告[R].1998
[60]贾明涛.基于过程模拟的回采方案优化技术及其在深井开采中的应用[博士学位论文].长沙：中南大学,2007
[61]马春德,汪为巍,陈枫.铜陵冬瓜山铜矿蛇纹岩力学性质实验报告[R].2005
[62]中国有色工程设计研究总院.冬瓜山初步设计[R].2003
[63]李仲奎,戴荣,姜逸明.FLAC3D分析的初始应力场生成及在大型地下洞室群计算中的应用[J].岩石力学与工程学报,2002,21(增2)：2876-2392
[64]唐辉明,晏鄂川,胡新丽.工程地质数值模拟的理论与方法[M].武汉：中国地质大学出版社,2001
[65]汪吉林,丁陈建,吴圣林.基于FLAC3D的复杂地貌三维地质建模[J].地质力学学报,2008,14(2)：149-157
[66]艾志雄,罗先启,刘波等.FLAC基本原理及其在边坡稳定性分析中的应用[J].灾害与防治工程,2006,60(1)：19-24
[67]魏学勇,武强,赵树贤.FLAC3D在矿井防治水中的应用[J].煤炭学报,2004,29(6)：704-707
[68]刘波,韩彦辉.FLAC原理、实例与应用指南[M].北京：人民交通出版社,2006
[69]张天宝.地下洞室群围岩稳定性分析[硕士学位论文].杭州：浙江大学, 2001
[70]赵兴东,杨天鸿,唐春安等.冲击地压下煤巷锚杆支护机理数值模拟研究[J].中国矿业,2004,13(3)：58-59
[71]梅松华,李文秀,盛谦.FLAC在岩土工程参数反演中的应用[J].矿冶工程,2000,20(4)：23-26
[72]肖红飞,何学秋,冯涛等.FLAC在煤炭变形破裂电耦合规律研究中的应用[J].煤炭学报,2004,29(6)：694-698
[73]朱建明,徐秉业,朱峰等.FLAC有限差分程序及其在矿山工程中的应用[J].中国矿业,20009(4)：78-81
[74]吴洪词,胡兴,包太.采场围岩稳定性的FLAC算法分析[J].矿山压力与顶板管理,2002,(4)：96-98
[75]张志扬.粘土中扩底桩承载性状及承载力预测研究[硕士学位论文].长沙：中南大学,2006
[76]彭文斌.FLAC3D实用教程[M].北京：机械工业出版社,2007
[77]陈育民,徐鼎平.FLAC/FLAC3D基础与工程实例[M].北京：中国水利水电出版社,2008
[78]王国富,王志忠.应用统计[M].长沙：中南大学出版社,2003
[2]何满潮,谢和平,彭苏萍等.深部开采岩体力学研究[J].岩石力学与工程学报,2005,24(16)：2803-2813
[3]孟稳权.冬瓜山矿床深部开采岩层控制的措施[J].矿业快报,2006,452(12)： 58-59
[4]杨承祥.深井金属矿床高效开采及地压监控技术研究[博士学位论文].长沙：中南大学,2007
[5]朱维申,何满潮.复杂条件下围岩稳定性与岩体动态施工力学[M].北京：科学出版社,1996
[6]李华晔.地下洞室围岩稳定性分析[M].北京：中国水利电力出版社,1999
[7]许传华,任青文,李瑞.地下工程围岩稳定性分析方法研究进展[J].金属矿山,2003,320(2)：34-37
[8]邵国建,卓家寿,章青.岩体稳定性分析与评判准则研究[J].岩石力学与工程学报,2003,22(5)：691-696
[9]王生全,夏玉成,樊怀仁.邢台显德汪井田煤层顶板稳定性评价与预测[J].西北地质,1997,18(4)：46-50
[10]杨仕教,丁德馨,古德生等.丰山铜矿采场稳定性的FUZZY综合评判[J].有色金属,2002,54(3)：99-102
[11]郭立,吴爱祥,周克平等.基于RES理论的岩体失稳模式判别及其智能实现[J].中国有色金属学报,2003,13(3)：749-753
[12]Hudson J A. Rock Engineering Systems:Theory and Practice[M]. Horwood Chichester,1992
[13]张保.深井开采矿柱稳定性分析与可视化验证[硕士学位论文].长沙：中南大学,2009
[14]宫凤强,李夕兵,高科.地下工程围岩稳定性分类的突变级数法研究[J].中南大学学报(自然科学版),2008,39(5)：1081-1086
[15]蔡美峰,孔广亚,贾立宏.岩体工程系统失稳的能量突变判断准则及其应用[J].北京科技大学学报,1997,19(4)：325-328
[16]朱永全.隧道稳定性位移判别准则[J].中国铁道科学,2001,22(6)： 80-83
[17]唐礼忠,潘长良,谢学斌等.冬瓜山铜矿深井开采岩爆危险区分析与预测[J].中南工业大学学报,2002,33(4)：335-338
[18]解联库,杨小聪,刘庆林等.基于霍克-布朗破坏准则的采空区稳定性分析[J].有色金属,2008,60(3)：109-111
[19]高峰,周科平,胡建华.采场稳定性的模糊物元评价模型及应用研究[J].采矿与安全工程学报,2006,23(2)：164-168
[20]李增学,吕大炜,王怀洪等.龙固井田主采煤层顶板稳定型特征研究[J].岩土工程学报,2008,30(4)：473-479
[21]冯玉国.灰色优化理论模型在地下工程围岩稳定性分类中的应用[J].岩土工程学报,1996,18(3)：62-66
[22]胡全舟,吴超,陈沅江.风险评价方法在地下工程围岩稳定性分级中的应用[J].地下空间与工程学报,2005,1(6)：874-877
[23]曹刚.康家湾矿深部难采矿体采场稳定性及安全开采技术研究[硕士学位论文].中南大学,2008
[24]刘培慧.基于应力边界法厚大矿体采场结构参数数值模拟优化研究[硕士学位论文].中南大学,2009
[25]布朗.工程岩体力学中的解析与数值计算方法(佘诗刚,王可均译)[M].北京：科学技术出版社,1997
[26]Chandrakant S.Dcsai,John T.Chritian.岩土工程数值方法(潘善德译)[M].北京：科学技术出版社
[27]姚振汉.边界元法[M].北京：高等教育出版社,2010
[28]Peter Hunter,Andrew Pul lan. FEM/BEM NOTES.Department of Engineering Science,The University of Auckland,New Zealand,2003
[29]江见鲸,何放龙,何益斌等.有限元法及其应用[M].北京：机械工业出版社,2006
[30]王泳嘉,邢纪波.离散元法及其在岩土力学中的应用[M].沈阳：东北大学出版社,1993
[31]孙广忠.岩体力学基础[M].北京：科学出版社,1983
[32]陈子光.岩石力学性质与构造应力场[M].北京：地质出版社,1986
[33]朱德仁.岩石工程破坏准则[J].煤炭学报,1994,19(1)：15-19
[34]王文星.岩体力学[M].长沙：中南大学出版社,2004
[35]吴爱祥,郭立,张卫锋.深井开采岩体破坏机理及工程控制方法综述[J].矿业研究与开发,2001,21(2)：4-7
[36]张孟喜,邓荣基,朱文喜.岩石破坏机理与岩石强度准则新探[J].甘肃科学学报,1996,8(1)56-61
[37]俞茂宏,昝月稳,范文等.20世纪岩石强度理论的发展-纪念Mohr-Coulomb强度理论100周年[J].岩石力学与工程学报,2000,19(5)：545-549
[38]Wittke W. Rock Mechanics:Theory and Applications with Case Histories [M]. Berlin:Springer-Verlag,1990,中译本：岩石力学(曾国熙译).杭州：浙江大学出版,1996
[39]Yu M H, Zhao J, Guan L W. Strength theory for rock and concrete: history, present situation and development [J]. Progress in Natural Science (English Edition),1998,8 (4):394-402
[40]白国良,梁冰.巷道岩体能量破坏机理探讨[J].矿山压力与顶板管理,2004,(1)：86-88
[41]杨官涛.地下采场结构参数优化及稳定性的能量突变分析[硕士学文论文].长沙：中南大学,2007
[42]付成华,陈胜宏.基于突变理论的地下工程洞室围岩失稳判据研究[J].岩土力学,2008,29(1)：167-172
[43]李忠华,官福海,潘一山.基于损伤理论的圆形巷道围岩应力场分析[J].岩土力学,2004,25(增刊)：160-163
[44]马青.基于损伤理论的圆形水工压力隧洞应力场分析[硕士学位论文].西安：长安大学,2009
[45]黄梦宏.岩石损伤演化模型和损伤临界值理论与实验研究[硕士学位论文].长沙：中南大学,2003
[46]凌建明,孙钧.建立在损伤应变空间的岩体破坏准则[J].同济大学学报,1995,23(5)：483-487
[47]徐卫亚,韦立德.岩石损伤统计本构模型的研究[J].岩石力学与工程学报,2002,21(6)：787-791
[48]田运生,张舒畅,李战军等.岩巷爆破威严损伤演化分析及损伤判据的建立[J].黄金,2005,26(8)：22-24
[49]李树茂,齐伟,刘红帅.岩体损伤力学理论进展[J].世界地质,2001,20(1)：72-78
[50]仵彦卿,丁卫华,蒲毅彬等.压缩条件下岩石密度损伤增量的CT动态观测[J].自然科学进展,2000,10(9)：830-835
[51]许传华.岩体破坏的非线性理论研究及应用[博士学位论文].南京：河海大学,2004
[52]孔广亚,蔡美峰,贾立宏.岩体工程系统破坏失稳的能量突变准则[J].有色矿冶,1996,(3)：1-3
[53]许传华,任青文,李瑞.围岩稳定的熵突变理论研究[J].岩石力学与工程学报,2004,23(12)：1992-1995
[54]杨官涛,李夕兵,王其胜等.地下采场失稳的能量突变判断准则及其应用[J].采矿与安全工程学报,2008,25(3)：268-276
[55]朱和玲.基于开采环境再造人工结构的稳定性与可开度研究[硕士学位论文].长沙：中南大学,2008
[56]北京有色冶金设计研究总院.铜陵有色金属(集团)公司冬瓜山铜矿可行性研究报告[R].1998
[57]北京有色冶金设计研究总院.安徽铜都铜业股份有限公司狮子山铜矿千米深井300万t级矿山开采条件研究成果[R].2001
[58]李冬青,王李管,等.深井硬岩大规模开采理论与技术——冬瓜山铜矿床开采研究与实践[M].北京：冶金工业出版社,2009
[59]北京有色冶金设计研究总院.狮子山铜矿冬瓜山矿床岩石力学试验报告[R].1998
[60]贾明涛.基于过程模拟的回采方案优化技术及其在深井开采中的应用[博士学位论文].长沙：中南大学,2007
[61]马春德,汪为巍,陈枫.铜陵冬瓜山铜矿蛇纹岩力学性质实验报告[R].2005
[62]中国有色工程设计研究总院.冬瓜山初步设计[R].2003
[63]李仲奎,戴荣,姜逸明.FLAC3D分析的初始应力场生成及在大型地下洞室群计算中的应用[J].岩石力学与工程学报,2002,21(增2)：2876-2392
[64]唐辉明,晏鄂川,胡新丽.工程地质数值模拟的理论与方法[M].武汉：中国地质大学出版社,2001
[65]汪吉林,丁陈建,吴圣林.基于FLAC3D的复杂地貌三维地质建模[J].地质力学学报,2008,14(2)：149-157
[66]艾志雄,罗先启,刘波等.FLAC基本原理及其在边坡稳定性分析中的应用[J].灾害与防治工程,2006,60(1)：19-24
[67]魏学勇,武强,赵树贤.FLAC3D在矿井防治水中的应用[J].煤炭学报,2004,29(6)：704-707
[68]刘波,韩彦辉.FLAC原理、实例与应用指南[M].北京：人民交通出版社,2006
[69]张天宝.地下洞室群围岩稳定性分析[硕士学位论文].杭州：浙江大学, 2001
[70]赵兴东,杨天鸿,唐春安等.冲击地压下煤巷锚杆支护机理数值模拟研究[J].中国矿业,2004,13(3)：58-59
[71]梅松华,李文秀,盛谦.FLAC在岩土工程参数反演中的应用[J].矿冶工程,2000,20(4)：23-26
[72]肖红飞,何学秋,冯涛等.FLAC在煤炭变形破裂电耦合规律研究中的应用[J].煤炭学报,2004,29(6)：694-698
[73]朱建明,徐秉业,朱峰等.FLAC有限差分程序及其在矿山工程中的应用[J].中国矿业,20009(4)：78-81
[74]吴洪词,胡兴,包太.采场围岩稳定性的FLAC算法分析[J].矿山压力与顶板管理,2002,(4)：96-98
[75]张志扬.粘土中扩底桩承载性状及承载力预测研究[硕士学位论文].长沙：中南大学,2006
[76]彭文斌.FLAC3D实用教程[M].北京：机械工业出版社,2007
[77]陈育民,徐鼎平.FLAC/FLAC3D基础与工程实例[M].北京：中国水利水电出版社,2008
[78]王国富,王志忠.应用统计[M].长沙：中南大学出版社,2003