工程科学学报.第41卷.第11期:1374-1386.2019年11月 Chinese Journal of Engineering,Vol.41,No.11:1374-1386,November 2019 D0L:10.13374/.issn2095-9389.2019.05.07.004,http://journals.ustb.edu.cn 高海拔寒区岩质边坡变形破坏机制研究现状及趋势 李长洪12,肖永刚1,2区,王宇1,2),卜磊2),侯志强12) 1)北京科技大学金属矿山高效开采与安全教育部重点实验室.北京1000832)北京科技大学土木与资源工程学院,北京100083 通信作者,E-mail:b20170013@xs.ustb.edu.cn 摘要高海拔寒区矿山岩质边坡变形破坏机制研究已取得一定的研究成果,但基于现行理论与技术还难以全面解决未来 高寒边坡失稳机理和灾害防控的所有问题,至今尚未建立起完善的高寒边坡开采研究体系和边坡稳定性判别标准.本文对高 寒岩质边坡变形破坏的室内岩石力学试验、边坡物理相似模拟、多场多相耦合数值模拟、变形破坏原位监测、高海拔寒区岩 质边坡失稳机理五个方面开展了大量的文献调研,总结高寒岩体变形破坏有关的研究成果,继而对存在的问题进行探讨并分 析当前研究的不足,总结出高寒岩质边坡变形破坏研究领域亟待解决的关键问题:一是开采扰动条件下高海拔寒区矿山边坡 岩体结构损伤劣化机制,二是冻融循环条件下流-固-气多相多场耦合边坡失稳时效特征与评价方法:并就未来高寒边坡变形 和破坏研究方向及发展趋势予以分析,指出开展不同应力路径冻融循环耦合作用下岩体结构损伤劣化机理研究,开展爆破采 动条件下高海拔寒区岩质边坡结构面致遗机制及边坡失稳破坏研究,开展地震荷载作用下高海拔寒区节理岩质边坡地震动 力响应及致灾规律研究,研究多场多相耦合条件下节理岩体损伤劣化机理,开展高海拔寒区矿山边坡抗寒多参量实时安全监 测及失稳预警技术研究五个方面是未来研究的趋势 关键词高寒边坡:冻融循环:爆破开采:动力响应;多相耦合;失稳机理 分类号P642.3 Review and prospects for understanding deformation and failure of rock slopes in cold regions with high altitude LI Chang-hong 2),XIAO Yong-gang2.WANG Yu2),BU Lei2),HOU Zhi-giang2) 1)Key Laboratory of Ministry of Education for High-Efficient Mining and Safety of Metal Mines,University of Science and Technology Beijing,Beijing 100083,China 2)School of Civil and Resource Engineering,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:b20170013@xs.ustb.edu.cn ABSTRACT Research on the deformation and failure mechanism of rock slopes in high-altitude cold areas has obtained certain results, but based on the current theory and technology,it is difficult to comprehensively solve problems related to the cold-slope instability mechanisms and disaster prevention and control.As yet,no overall research system for high-altitude slope mining or criteria for slope stability have been established.In this paper,based on an extensive literature review,five measures of the deformation and failure of alpine rock slopes were presented,including the indoor rock mechanics test,simulation of physically similar slopes,multi-field multi- phase coupled numerical simulation,in-situ monitoring of deformation and damage,and the instability mechanism of rock slopes in high-altitude cold areas.After summarizing the research results related to the deformation and failure of alpine rock masses,existing problems were discussed and current research deficiencies were analyzed.The key problems that require urgent solutions in the research of the deformation and failure of alpine rock slopes were summarized.The first problem is the damage mechanism of a rock mass in a 收稿日期:2019-05-07 基金项目:国家重点研发计划课题资助项目(2018YFC0808402)
高海拔寒区岩质边坡变形破坏机制研究现状及趋势 李长洪1,2),肖永刚1,2) 苣,王 宇1,2),卜 磊1,2),侯志强1,2) 1) 北京科技大学金属矿山高效开采与安全教育部重点实验室,北京 100083 2) 北京科技大学土木与资源工程学院,北京 100083 苣通信作者,E-mail:b20170013@xs.ustb.edu.cn 摘 要 高海拔寒区矿山岩质边坡变形破坏机制研究已取得一定的研究成果,但基于现行理论与技术还难以全面解决未来 高寒边坡失稳机理和灾害防控的所有问题,至今尚未建立起完善的高寒边坡开采研究体系和边坡稳定性判别标准. 本文对高 寒岩质边坡变形破坏的室内岩石力学试验、边坡物理相似模拟、多场多相耦合数值模拟、变形破坏原位监测、高海拔寒区岩 质边坡失稳机理五个方面开展了大量的文献调研,总结高寒岩体变形破坏有关的研究成果,继而对存在的问题进行探讨并分 析当前研究的不足,总结出高寒岩质边坡变形破坏研究领域亟待解决的关键问题:一是开采扰动条件下高海拔寒区矿山边坡 岩体结构损伤劣化机制,二是冻融循环条件下流−固−气多相多场耦合边坡失稳时效特征与评价方法;并就未来高寒边坡变形 和破坏研究方向及发展趋势予以分析,指出开展不同应力路径冻融循环耦合作用下岩体结构损伤劣化机理研究,开展爆破采 动条件下高海拔寒区岩质边坡结构面致溃机制及边坡失稳破坏研究,开展地震荷载作用下高海拔寒区节理岩质边坡地震动 力响应及致灾规律研究,研究多场多相耦合条件下节理岩体损伤劣化机理,开展高海拔寒区矿山边坡抗寒多参量实时安全监 测及失稳预警技术研究五个方面是未来研究的趋势. 关键词 高寒边坡;冻融循环;爆破开采;动力响应;多相耦合;失稳机理 分类号 P642.3 Review and prospects for understanding deformation and failure of rock slopes in cold regions with high altitude LI Chang-hong1,2) ,XIAO Yong-gang1,2) 苣 ,WANG Yu1,2) ,BU Lei1,2) ,HOU Zhi-qiang1,2) 1) Key Laboratory of Ministry of Education for High-Efficient Mining and Safety of Metal Mines, University of Science and Technology Beijing, Beijing 100083, China 2) School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China 苣 Corresponding author, E-mail: b20170013@xs.ustb.edu.cn ABSTRACT Research on the deformation and failure mechanism of rock slopes in high-altitude cold areas has obtained certain results, but based on the current theory and technology, it is difficult to comprehensively solve problems related to the cold-slope instability mechanisms and disaster prevention and control. As yet, no overall research system for high-altitude slope mining or criteria for slope stability have been established. In this paper, based on an extensive literature review, five measures of the deformation and failure of alpine rock slopes were presented, including the indoor rock mechanics test, simulation of physically similar slopes, multi-field multiphase coupled numerical simulation, in-situ monitoring of deformation and damage, and the instability mechanism of rock slopes in high-altitude cold areas. After summarizing the research results related to the deformation and failure of alpine rock masses, existing problems were discussed and current research deficiencies were analyzed. The key problems that require urgent solutions in the research of the deformation and failure of alpine rock slopes were summarized. The first problem is the damage mechanism of a rock mass in a 收稿日期: 2019−05−07 基金项目: 国家重点研发计划课题资助项目(2018YFC0808402) 工程科学学报,第 41 卷,第 11 期:1374−1386,2019 年 11 月 Chinese Journal of Engineering, Vol. 41, No. 11: 1374−1386, November 2019 DOI:10.13374/j.issn2095-9389.2019.05.07.004; http://journals.ustb.edu.cn
李长洪等:高海拔寒区岩质边坡变形破坏机制研究现状及趋势 ·1375 high-altitude cold area under mining disturbance conditions.The second problem is the aging characteristics and evaluation methods of the instability of a flow-solid-gas multi-phase multi-field coupled slope under freeze-thaw cycles.The future research direction and development trends in the deformation and failure of rock slopes in cold regions with high altitude were also analyzed.Research should be conducted on the following:(1)the damage degradation mechanism of rock masses with different stress paths coupled with freeze-thaw cycles,(2)the structural plane collapse mechanism and instability of rock slopes in high-altitude cold areas under the condition of blast mining,(3)the dynamic response and disaster occurrence law of jointed rock slopes in high-altitude cold areas under earthquake loading.(4)the mechanism of damage deterioration of jointed rock masses under multi-field and multi-phase coupling conditions,and(5)real-time safety monitoring and early-warning technology regarding the instability of multi-parameter cold resistance of mine slopes in high-altitude cold areas.These five research areas constitute the trends of future research. KEY WORDS alpine slope;freeze-thaw cycle;blast mining;dynamic response;multi-phase coupling;failure mechanism 我国低海拔金属矿资源量减少,而其需求量却 岩体边坡失稳问题亟待解决 不断增加,《全国矿产资源规划(2016一2020)》四 1高寒岩质边坡变形和破坏研究现状调研 指出金属矿产资源开发应该逐步向西部地区进 军,将高海拔寒区金属矿产资源开发上升到国家 高寒岩质边坡多数为岩性组合变化多、高陡 战略高度.高海拔高寒地区低温、低气压、冻融循 坡、强卸荷、多裂隙、冻融作用明显的边坡,一旦 环、生态脆弱等特点,导致矿山开采条件差、灾变 出现了变形和失稳,会严重影响工程各方面的进 机理复杂、人机功效低、救援难度大,而现有的开 展,本文从高寒岩体变形破坏室内岩石力学试 采理论与技术不能适用于高海拔寒区金属矿产资 验、相似模型试验、数值模拟、高寒边坡原位监测 源安全开采,迫切需要建立新理论,开发新技术与 以及失稳机理5个方面进行总结 装备 1.1高海拔寒区岩体室内岩石力学试验研究现状 爆破开采、冻融循环、地震等作用引起岩质边 高寒岩体变形破坏室内试验包括冻融循环条 坡地质灾害,致使高寒地区露天矿山岩质边坡的 件下的静态试验、爆破开采等应力扰动下的动态 破坏问题尤为突出.由于露天金属矿的开采致使 试验等 边坡长期暴露在自然环境中,受温度、渗流、采 1.1.1静态试验研究 矿、冻融等因素的影响,随着采动后暴露时间的延 ince和Fenerl)利用卡帕多西亚l0个火山碎 长,岩质边坡结构完整性和强度将逐渐下降,一旦 屑岩样品进行冻融(freeze-thaw简写为FT)循环 诱发灾害,将会给国家经济建设和人民生命财产 试验,建立强度损失百分比统计模型预测出火山 带来巨大损失.例如在2013年,我国西藏甲玛矿 碎屑岩单轴抗压强度损失百分比,具有较高的可 由于冻融循环和爆破开采等动载作用导致岩体劣 靠性.Jia等通过不同饱和度的砂岩试样的冻胀 化失稳,引起200多万方边坡塌方量,83名现场作 变形与单轴拉伸变形对比(如图1),不仅说明砂岩 业人员被埋,“3.29”塌方给人民的生命财产造成巨 冻胀变形与饱和度有关,而且证明超过一定饱和 大损失) 度之后的最大冻胀变形将会存在不可逆的残余变 目前,有关寒区煤矿露天边坡工程和金属矿 形.Han等阿研究砂岩试样在不同化学溶液和快 露天矿山边坡工程岩体在冻融循环条件下岩土体 速FT循环耦合作用下的损伤机理,结果表明砂岩 介质物理力学性质的研究多有报道.露天煤矿边 试样的断裂韧性和强度随着FT循环次数的增加 坡岩体强度较低,主要由沉积岩构成,层理明显且 出现不同程度的劣化.Chen等对软岩材料试件 夹层多,岩体较破碎,而露天金属矿边坡岩体强度 进行了T循环试验,建立冻融边坡稳定性的计算 较高,主要由岩浆岩、变质岩构成,断层、节理发 模型,结果表明含水率和冻融循环次数对软岩的 育,地质条件复杂.有关高寒高海拔地区金属矿露 力学参数有显着影响.Lu0等m采集某寒区铁矿边 天矿山边坡节理岩体在开采扰动条件下的力学响 坡辉绿岩,对比辉绿岩和变质辉绿岩的冻融系数 应、结构面灾变扩展机理的研究较少.随着“一带 (K)和岩石强度损失率(R)等指标,以此研究两 一路”建设的不断深入,确保露天矿山边坡开采边 种岩样冻融损伤劣化规律. 坡安全与提高经济效益之间的矛盾日渐明显,有 Ghobadi和Babazadeh对伊朗西部某地区砂 关高寒高海拔地区露天金属矿山开采相关的节理 岩试样进行了60次T循环试验,结果表明同一
high-altitude cold area under mining disturbance conditions. The second problem is the aging characteristics and evaluation methods of the instability of a flow –solid –gas multi-phase multi-field coupled slope under freeze –thaw cycles. The future research direction and development trends in the deformation and failure of rock slopes in cold regions with high altitude were also analyzed. Research should be conducted on the following: (1) the damage degradation mechanism of rock masses with different stress paths coupled with freeze –thaw cycles, (2) the structural plane collapse mechanism and instability of rock slopes in high-altitude cold areas under the condition of blast mining, (3) the dynamic response and disaster occurrence law of jointed rock slopes in high-altitude cold areas under earthquake loading, (4) the mechanism of damage deterioration of jointed rock masses under multi-field and multi-phase coupling conditions, and (5) real-time safety monitoring and early-warning technology regarding the instability of multi-parameter cold resistance of mine slopes in high-altitude cold areas. These five research areas constitute the trends of future research. KEY WORDS alpine slope;freeze-thaw cycle;blast mining;dynamic response;multi-phase coupling;failure mechanism 我国低海拔金属矿资源量减少,而其需求量却 不断增加,《全国矿产资源规划(2016—2020)》[1] 指出金属矿产资源开发应该逐步向西部地区进 军,将高海拔寒区金属矿产资源开发上升到国家 战略高度. 高海拔高寒地区低温、低气压、冻融循 环、生态脆弱等特点,导致矿山开采条件差、灾变 机理复杂、人机功效低、救援难度大,而现有的开 采理论与技术不能适用于高海拔寒区金属矿产资 源安全开采,迫切需要建立新理论,开发新技术与 装备. 爆破开采、冻融循环、地震等作用引起岩质边 坡地质灾害,致使高寒地区露天矿山岩质边坡的 破坏问题尤为突出. 由于露天金属矿的开采致使 边坡长期暴露在自然环境中,受温度、渗流、采 矿、冻融等因素的影响,随着采动后暴露时间的延 长,岩质边坡结构完整性和强度将逐渐下降,一旦 诱发灾害,将会给国家经济建设和人民生命财产 带来巨大损失. 例如在 2013 年,我国西藏甲玛矿 由于冻融循环和爆破开采等动载作用导致岩体劣 化失稳,引起 200 多万方边坡塌方量,83 名现场作 业人员被埋,“3.29”塌方给人民的生命财产造成巨 大损失[2] . 目前,有关寒区煤矿露天边坡工程和金属矿 露天矿山边坡工程岩体在冻融循环条件下岩土体 介质物理力学性质的研究多有报道. 露天煤矿边 坡岩体强度较低,主要由沉积岩构成,层理明显且 夹层多,岩体较破碎,而露天金属矿边坡岩体强度 较高,主要由岩浆岩、变质岩构成,断层、节理发 育,地质条件复杂. 有关高寒高海拔地区金属矿露 天矿山边坡节理岩体在开采扰动条件下的力学响 应、结构面灾变扩展机理的研究较少. 随着“一带 一路”建设的不断深入,确保露天矿山边坡开采边 坡安全与提高经济效益之间的矛盾日渐明显,有 关高寒高海拔地区露天金属矿山开釆相关的节理 岩体边坡失稳问题亟待解决. 1 高寒岩质边坡变形和破坏研究现状调研 高寒岩质边坡多数为岩性组合变化多、高陡 坡、强卸荷、多裂隙、冻融作用明显的边坡,一旦 出现了变形和失稳,会严重影响工程各方面的进 展. 本文从高寒岩体变形破坏室内岩石力学试 验、相似模型试验、数值模拟、高寒边坡原位监测 以及失稳机理 5 个方面进行总结. 1.1 高海拔寒区岩体室内岩石力学试验研究现状 高寒岩体变形破坏室内试验包括冻融循环条 件下的静态试验、爆破开采等应力扰动下的动态 试验等. 1.1.1 静态试验研究 İnce 和 Fener[3] 利用卡帕多西亚 10 个火山碎 屑岩样品进行冻融(freeze–thaw 简写为 FT)循环 试验,建立强度损失百分比统计模型预测出火山 碎屑岩单轴抗压强度损失百分比,具有较高的可 靠性. Jia 等[4] 通过不同饱和度的砂岩试样的冻胀 变形与单轴拉伸变形对比(如图 1),不仅说明砂岩 冻胀变形与饱和度有关,而且证明超过一定饱和 度之后的最大冻胀变形将会存在不可逆的残余变 形. Han 等[5] 研究砂岩试样在不同化学溶液和快 速 FT 循环耦合作用下的损伤机理,结果表明砂岩 试样的断裂韧性和强度随着 FT 循环次数的增加 出现不同程度的劣化. Chen 等[6] 对软岩材料试件 进行了 FT 循环试验,建立冻融边坡稳定性的计算 模型,结果表明含水率和冻融循环次数对软岩的 力学参数有显着影响. Luo 等[7] 采集某寒区铁矿边 坡辉绿岩,对比辉绿岩和变质辉绿岩的冻融系数 (Kf)和岩石强度损失率(Rσ)等指标,以此研究两 种岩样冻融损伤劣化规律. Ghobadi 和 Babazadeh[8] 对伊朗西部某地区砂 岩试样进行了 60 次 FT 循环试验,结果表明同一 李长洪等: 高海拔寒区岩质边坡变形破坏机制研究现状及趋势 · 1375 ·
1376 工程科学学报.第41卷,第11期 砂岩单轴直接拉伸应力-应变曲线 0.8 弹性变形阶段 塑性变形阶段 0.6 0.4 50%饱和试样的 最大冻胀变形 90%饱和试样的 0.2 最大冻胀变形 0 0.001 0.002 0.003 应变 图1砂岩在不同饱和度下的冻胀变形与单轴拉伸变形对比州 Fig.I Comparison of frost heave deformation and uniaxial tensile deformation of sandstone at different saturation levels! 地层的岩石风化特性存在较大差异,而且干湿条 明,冻结状态下的峰值抗压强度、黏聚力和内摩擦 件对砂岩特性有较大影响.Fang等9指出随着腐 角均有不同程度的提高:②5、10、20和40次冻融 蚀溶液FT最低温度和pH值的降低,黄砂岩抗压 后红砂岩力学试验结果表明,随着FT循环试验次 强度和弹性模量均减小,说明承载力和抗变形能 数的增加,红砂岩的纵波波速、弹性模量和抗压强 力随着延性的增大而减小.De Kock等o探讨了 度逐渐减小;③采用0°、45°和90°三个角度制备三 微米CT技术在研究冻融过程中孔隙尺度动态变 种单裂隙类砂岩岩样,研究冻融循环和周期荷载 化的可能性,试验表明多孔石灰岩冷却至-9.7℃ 耦合下不同冻融循环的破坏特征,表明岩体的损 时(如图2),裂缝网络发育,此时冰结晶并释放潜 伤及强度与岩体裂隙、冻融作用、周期荷载均有 热,吸水路径以及岩石纹理结构影响着裂缝生成 关联,并使用缺陷面积来描述裂隙损伤区 的位置. 综上所述,冻融循环作用使岩体的力学特性 杨更社等叫,张慧梅等1和申艳军等3对 发生了显著的变化,但是关于冻融循环条件下岩 饱和砂岩、红砂岩、单裂隙类砂岩岩样进行了不 体结构劣化的研究不足,尤其是针对高寒边坡采 同冻融次数的FT循环试验,得出的结论主要有: 动条件下岩体结构劣化研究不足,应开展有针对 ①不同冻结温度下的饱和砂岩三轴压缩试验表 性的室内力学试验,特别是接近工程实际和环境 20 (a) 15 10 0 潜热释放引起的放热峰 -5 冻结温度为-5℃ -10 冻结温度为-15℃ -15 -20 0 50b 100015002000250030003500 时间俗 (b) 中点 (c) 中点A 5 mm 图2冻结温度设定在-5℃和-15℃时读出数据四(a)岩样中传感器温度随时间的变化:(b)放热峰前CT切片及裂缝半高全宽值曲线:(©)放热 峰后CT切片及裂缝半高全宽值曲线 Fig.2 Readouts of freezing temperature with set points-5Cand-15:(a)the temperature of the sensor in the rock sampler as a function of time; (b)CT slice and Full Width Half Maximum (fwhm)before the exothermal peak;(c)CT slice and Full Width Half Maximum(fwhm)after the exothermal peak
地层的岩石风化特性存在较大差异,而且干湿条 件对砂岩特性有较大影响. Fang 等[9] 指出随着腐 蚀溶液 FT 最低温度和 pH 值的降低,黄砂岩抗压 强度和弹性模量均减小,说明承载力和抗变形能 力随着延性的增大而减小. De Kock 等[10] 探讨了 微米 CT 技术在研究冻融过程中孔隙尺度动态变 化的可能性,试验表明多孔石灰岩冷却至−9.7 ℃ 时(如图 2),裂缝网络发育,此时冰结晶并释放潜 热,吸水路径以及岩石纹理结构影响着裂缝生成 的位置. 杨更社等[11] ,张慧梅等[12] 和申艳军等[13−14] 对 饱和砂岩、红砂岩、单裂隙类砂岩岩样进行了不 同冻融次数的 FT 循环试验,得出的结论主要有: ①不同冻结温度下的饱和砂岩三轴压缩试验表 明,冻结状态下的峰值抗压强度、黏聚力和内摩擦 角均有不同程度的提高;②5、10、20 和 40 次冻融 后红砂岩力学试验结果表明,随着 FT 循环试验次 数的增加,红砂岩的纵波波速、弹性模量和抗压强 度逐渐减小;③采用 0°、45°和 90°三个角度制备三 种单裂隙类砂岩岩样,研究冻融循环和周期荷载 耦合下不同冻融循环的破坏特征,表明岩体的损 伤及强度与岩体裂隙、冻融作用、周期荷载均有 关联,并使用缺陷面积来描述裂隙损伤区. 综上所述,冻融循环作用使岩体的力学特性 发生了显著的变化,但是关于冻融循环条件下岩 体结构劣化的研究不足,尤其是针对高寒边坡采 动条件下岩体结构劣化研究不足,应开展有针对 性的室内力学试验,特别是接近工程实际和环境 0 0.001 0.002 0.003 0.8 0.4 0.2 0 0.6 应变 应力/MPa 50%饱和试样的 最大冻胀变形 90%饱和试样的 最大冻胀变形 砂岩单轴直接拉伸应力−应变曲线 弹性变形阶段 塑性变形阶段 图 1 砂岩在不同饱和度下的冻胀变形与单轴拉伸变形对比[4] Fig.1 Comparison of frost heave deformation and uniaxial tensile deformation of sandstone at different saturation levels[4] 0 500 1000 1500 2000 2500 3000 3500 时间/s 温度/°C −20 −5 −10 −15 20 15 10 0 5 冻结温度为−5 ℃ 冻结温度为−15 ℃ (a) (b) 5 mm 潜热释放引起的放热峰 中点 (c) 中点 边 缘 点 边 缘 点 图 2 冻结温度设定在−5 ℃ 和−15 ℃ 时读出数据[10] . (a)岩样中传感器温度随时间的变化;(b)放热峰前 CT 切片及裂缝半高全宽值曲线;(c)放热 峰后 CT 切片及裂缝半高全宽值曲线 Fig.2 Readouts of freezing temperature with set points −5 ℃ and −15 ℃[10] : (a) the temperature of the sensor in the rock sampler as a function of time; (b) CT slice and Full Width Half Maximum (fwhm) before the exothermal peak; (c) CT slice and Full Width Half Maximum (fwhm) after the exothermal peak · 1376 · 工程科学学报,第 41 卷,第 11 期
李长洪等:高海拔寒区岩质边坡变形破坏机制研究现状及趋势 1377 下的冻融循环扰动试验 环作用以外,还会频繁受到爆破开采、大型机械化 1.1.2动态试验研究 施工以及地震等影响(如图3所示),将产生较大 高海拔寒区露天金属矿山边坡除受到冻融循 的应力扰动或冲击荷载,使节理裂隙扩展发育 动力波源: 渗流场 岩体A岩体B 温度场 爆破开采、 冻融循环 大型机械施工, 地震、 应力波 爆炸冲击等 波的特性: 44 纵波和横被、 传播 面波、 频率 振幅. 微观结构影响 接触面效应 岩体 单一节理效应 边界效应等 节理 尺度 效 ,节理组效应 →图 节理岩体效应 图3岩体工程受环境因素及内在因素影响) Fig.3 Effects of environmental and internal factors on rock mass engineering 为研究高海拔寒区岩体工程受爆破或地震等 并不是一种耦合关系,所以可分开单独考虑.李维 动力荷载的影响,通过使用分离式Hopkinson压杆 光2设计了可以模拟不同爆破方向的爆破振动试 冲击(SHPB)测试和传统的单轴抗压强度(UCS)测 验台,以某露天层理岩质边坡为工程背景系统研 试,Ke等16、Wang等n7-ly、Zhou等2o、Li等和 究了爆破采动对顺层边坡失稳的影响.徐拴海2 Zhang等21得出的主要结论如下: 采集某岩质边坡多种砂岩试样进行FT循环试验, ①冲击加载试验结果表明,应力-应变曲线均 考虑爆破开挖对冻岩边坡的影响,定义了岩体的 经历了致密化、弹性、屈服和破坏四个阶段,随着 冻融扰动系数D,分析了冻融扰动和爆破采动对 冻融循环次数的增加,动态强度、杨氏模量和动态 岩体的劣化影响. 峰值应力减小,能量吸收随着冻融循环次数的增 综上所述,高寒岩体实验室动态试验主要采 加而增加 用的方法是SHPB,且主要针对砂岩等多孔隙、易 ②随循环次数的增加,砂岩孔隙结构分布趋 变形岩性,专门针对高寒露天金属矿山边坡岩体 势呈现多尺度特征,其孔隙率随着冻融作用的持 在爆破采动下的动态力学试验开展较少,应针对 续影响而增大,弹性模量逐渐减小,峰值应变不断 高寒边坡研发可以模拟不同爆破荷载大小、不同 增大. 爆破方向、不同爆破荷载加载方式的爆破振动试 ③砂岩的饱和质量、孔隙度和宏观损伤随着 验台,采用相似原理设计并制作符合工程实际的 冻融循环次数的增加而增大 边坡岩体复杂结构模型,进行爆破振动台模拟试 ④FT循环作用后,红砂岩的力学性能有明显 验,全面模拟爆破振动效应 的劣化,其UCS和模量的降低比静载荷作用下更 1.2边坡物理相似模型试验研究现状 为显著 国内外许多学者在研究边坡动载作用下的响 Liu等21采用分离式Hopkinson压杆冲击 应时采用物理模型振动台试验,该方法能够较为 (SHPB)系统进行了花岗岩巴西劈裂试验,提出了 准确的再现整个地质体的变形破坏过程,由于其 考虑加载速率和冻融循环的动态拉伸强度预测公 成本低、周期短、成果形象直观且对影响因素能 式.Ma等24对经过不同次数的FT周期后的软岩 进行重复分析,因此,国内外学者广泛使用边坡模 和砂质泥岩进行了SHPB动态压缩试验,结果表明 型试验来研究边坡的变形和破坏 两种岩样单轴动态抗压强度随着FT周期次数的 Yang等7以汶川地震为背景通过边坡模型 增加呈对数下降趋势.由于岩体爆破和冻融循环 试验研究地震滑坡发生机理及其灾害防控.Huang 不一定同时发生,在时间上存在先后的叠加作用, 等2!通过大型振动台模拟地震荷载作用下边坡变
下的冻融循环扰动试验. 1.1.2 动态试验研究 高海拔寒区露天金属矿山边坡除受到冻融循 环作用以外,还会频繁受到爆破开采、大型机械化 施工以及地震等影响(如图 3 所示),将产生较大 的应力扰动或冲击荷载,使节理裂隙扩展发育. 为研究高海拔寒区岩体工程受爆破或地震等 动力荷载的影响,通过使用分离式 Hopkinson 压杆 冲击(SHPB)测试和传统的单轴抗压强度(UCS)测 试,Ke 等[2, 16]、Wang 等 [17−19]、Zhou 等[20]、Li 等[21] 和 Zhang 等[22] 得出的主要结论如下: ① 冲击加载试验结果表明,应力−应变曲线均 经历了致密化、弹性、屈服和破坏四个阶段,随着 冻融循环次数的增加,动态强度、杨氏模量和动态 峰值应力减小,能量吸收随着冻融循环次数的增 加而增加. ② 随循环次数的增加,砂岩孔隙结构分布趋 势呈现多尺度特征,其孔隙率随着冻融作用的持 续影响而增大,弹性模量逐渐减小,峰值应变不断 增大. ③ 砂岩的饱和质量、孔隙度和宏观损伤随着 冻融循环次数的增加而增大. ④ FT 循环作用后,红砂岩的力学性能有明显 的劣化,其 UCS 和模量的降低比静载荷作用下更 为显著. Liu 等 [23] 采 用 分 离 式 Hopkinson 压 杆 冲 击 (SHPB)系统进行了花岗岩巴西劈裂试验,提出了 考虑加载速率和冻融循环的动态拉伸强度预测公 式. Ma 等[24] 对经过不同次数的 FT 周期后的软岩 和砂质泥岩进行了 SHPB 动态压缩试验,结果表明 两种岩样单轴动态抗压强度随着 FT 周期次数的 增加呈对数下降趋势. 由于岩体爆破和冻融循环 不一定同时发生,在时间上存在先后的叠加作用, 并不是一种耦合关系,所以可分开单独考虑. 李维 光[25] 设计了可以模拟不同爆破方向的爆破振动试 验台,以某露天层理岩质边坡为工程背景系统研 究了爆破采动对顺层边坡失稳的影响. 徐拴海[26] 采集某岩质边坡多种砂岩试样进行 FT 循环试验, 考虑爆破开挖对冻岩边坡的影响,定义了岩体的 冻融扰动系数 Dn,分析了冻融扰动和爆破采动对 岩体的劣化影响. 综上所述,高寒岩体实验室动态试验主要采 用的方法是 SHPB,且主要针对砂岩等多孔隙、易 变形岩性,专门针对高寒露天金属矿山边坡岩体 在爆破采动下的动态力学试验开展较少,应针对 高寒边坡研发可以模拟不同爆破荷载大小、不同 爆破方向、不同爆破荷载加载方式的爆破振动试 验台,采用相似原理设计并制作符合工程实际的 边坡岩体复杂结构模型,进行爆破振动台模拟试 验,全面模拟爆破振动效应. 1.2 边坡物理相似模型试验研究现状 国内外许多学者在研究边坡动载作用下的响 应时采用物理模型振动台试验,该方法能够较为 准确的再现整个地质体的变形破坏过程,由于其 成本低、周期短、成果形象直观且对影响因素能 进行重复分析,因此,国内外学者广泛使用边坡模 型试验来研究边坡的变形和破坏. Yang 等[27] 以汶川地震为背景通过边坡模型 试验研究地震滑坡发生机理及其灾害防控. Huang 等[28] 通过大型振动台模拟地震荷载作用下边坡变 动力波源: 岩体A 岩体B 爆破开采、 大型机械施工、 地震、 爆炸冲击等. 波的特性: 纵波和横波、 面波、 频率、 振幅. 渗流场 应力波 传播 岩体 节理 尺度 效应 微观结构影响 单一节理效应 节理组效应 节理岩体效应 接触面效应 边界效应等 温度场 最大主 冻融循环 应力σ1 最小主 应力σ3 图 3 岩体工程受环境因素及内在因素影响[15] Fig.3 Effects of environmental and internal factors on rock mass engineering[15] 李长洪等: 高海拔寒区岩质边坡变形破坏机制研究现状及趋势 · 1377 ·
·1378 工程科学学报,第41卷,第11期 形破坏来研究硬岩边坡的失稳机理.Che等采 用振动台试验分析了高陡边坡在地震作用下非连 渗透压力 岩 自重应力 续节理对边坡的影响规律.Fan等Bo认为边坡的 体应变 动力响应与地震能量的不均匀分布有重大关联 渗流场 渗透压力 应力场 Zhang等B川通过进行泥岩边坡离心振动台物理模 拟试验,研究了地震荷载作用下泥岩边坡的动力 透 孔隙用力 响应规律.Lⅰ等四通过两个相似的边坡物理模型 变化 渗透压力 渗透系数 的振动台试验对比,研究了孔隙水压力、加速度、 频谱等特征与边坡场地动力响应特征之间的关联 水相变 冰 规律.黄震配制了流固耦合相似材料,开展了 冰水相变 任移 围岩渗透破坏灾变演化规律的物理模型试验,研 机械能 究了开挖扰动对围岩渗流特征的影响.杨忠平等 正温 负温 温度场 采用振动台架试验并结合通用离散单元法程序 图4岩-水-冰系统与多场耦合关系B8 (UDEC)数值分析研究反复微小地震对两种岩质 边坡(顺层和反倾)的动力响应.王来贵等通过 Fig.4 Relationship of coupled rock-water-ice system 边坡振动台模拟试验揭示边坡内部应力状态的调 定律并考虑水冰相变的静力平衡原理,推求低温 整变化规律.冯细霞等通过振动台试验研究具 冻结岩体热-水-力(thermo-hydro-mechanical简写 有两组正交节理的边坡动力响应,揭示了反倾岩 为THM)完全耦合控制方程,并应用于寒区隧道的 质边坡在不同波形、不同振幅和频率下的动力响 稳定性分析.刘乃飞2建立了高海拔寒区裂隙岩 应规律 体的应力-渗流-温度-化学四场耦合数学模型,基 综上所述,岩质边坡在动载作用下的振动台 于自行开发的4G2017岩土工程仿真软件建立典 物理模型试验取得了较多的研究成果,但针对高 型寒区岩体工程模型.李国峰等推导寒区岩体 海拔寒区岩质边坡在地震荷载、爆破采动等动载 含相变的温度场-渗流场-应力场三场耦合简化算 作用下的动力响应特征研究较少,存在的主要问 法,通过FLAC3D对室内岩样力学试验仿真验证了 题有:①传统的物理模型较大,而高寒边坡存在冻 该算法的正确性.Shen等对准砂岩试样进行了 融环境,如何将较大尺寸的边坡模型进行冻融之 从20℃到-20℃冻融循环试验,建立了考虑潜热 后再加载,这是进行高寒边坡物理模型试验的一 效应的表观热容模型(AHCM),并将该模型嵌入 大难点;②传统的边坡模型振动台试验采用接触 COMSOL程序中,模拟岩体内部温度场的变化规 式手段测量变形,当边坡模型内发生较大变形时, 律.李滨锷建立非饱和土固-液-气三相耦合模 接触测量装置便与监测对象分离,无法继续获取 型,揭示了非饱和边坡渗流、气体运移对边坡稳定 该点监测数据,需要采用新的监测手段或研发新 系数的影响.陈益峰等啊建立了岩土固-液-气三 装备来解决边坡模型振动台试验中这一问题; 相THM完全耦合模型.班改革和戴剑勇6建立 ③如何在边坡物理模型中构建真实复杂的节理, 了非饱和岩体基于固-液-气三相系统守恒定律的 如何模拟控制高寒边坡模型的含水量和含冰量也 热一流-固完全耦合数学模型 是边坡振动台物理模拟需要研究的难点之一 从以上研究成果来看,国内外学者针对寒区 13高海拔寒区岩体结构数值模拟研究现状 岩体建立了多个耦合模型,但目前关于高寒岩体 高海拔寒区岩体工程破坏是由各种化学和力 的多场耦合模型多是在冻土力学的理论基础上, 学机制多场耦合作用造成的,影响高海拔寒区岩 将岩体看作各向同性的孔隙介质,对于寒区露天 体的因素众多,边坡稳定性研究非常复杂,但是边 矿山岩体各向异性、裂隙内水冰相变、流-固-气 坡失稳的根源是水-热-力多场耦合作用刃,寒区 多相多场耦合模型的研究较少,边坡多场多相耦 金属矿山开采过程中冻融温度场、应力场、渗流 合致灾过程方面的研究尚不完善 场耦合作用影响岩质边坡稳定性,与常温多场耦 1.4高海拔寒区岩质边坡变形破坏原位监测研究 合不同的是寒区岩体内会出现水冰相变,寒区岩 现状 体工程岩-水-冰与多场耦合关系如图4所示 为了更好地研究高寒边坡的变形破坏机理,国 Huang等9-)根据能量守恒定律、质量守恒 内外专家学者对寒区边坡使用多种手段进行了监测
形破坏来研究硬岩边坡的失稳机理. Che 等[29] 采 用振动台试验分析了高陡边坡在地震作用下非连 续节理对边坡的影响规律. Fan 等[30] 认为边坡的 动力响应与地震能量的不均匀分布有重大关联. Zhang 等[31] 通过进行泥岩边坡离心振动台物理模 拟试验,研究了地震荷载作用下泥岩边坡的动力 响应规律. Li 等[32] 通过两个相似的边坡物理模型 的振动台试验对比,研究了孔隙水压力、加速度、 频谱等特征与边坡场地动力响应特征之间的关联 规律. 黄震[33] 配制了流固耦合相似材料,开展了 围岩渗透破坏灾变演化规律的物理模型试验,研 究了开挖扰动对围岩渗流特征的影响. 杨忠平等[34] 采用振动台架试验并结合通用离散单元法程序 (UDEC)数值分析研究反复微小地震对两种岩质 边坡(顺层和反倾)的动力响应. 王来贵等[35] 通过 边坡振动台模拟试验揭示边坡内部应力状态的调 整变化规律. 冯细霞等[36] 通过振动台试验研究具 有两组正交节理的边坡动力响应,揭示了反倾岩 质边坡在不同波形、不同振幅和频率下的动力响 应规律. 综上所述,岩质边坡在动载作用下的振动台 物理模型试验取得了较多的研究成果,但针对高 海拔寒区岩质边坡在地震荷载、爆破采动等动载 作用下的动力响应特征研究较少,存在的主要问 题有:①传统的物理模型较大,而高寒边坡存在冻 融环境,如何将较大尺寸的边坡模型进行冻融之 后再加载,这是进行高寒边坡物理模型试验的一 大难点;②传统的边坡模型振动台试验采用接触 式手段测量变形,当边坡模型内发生较大变形时, 接触测量装置便与监测对象分离,无法继续获取 该点监测数据,需要采用新的监测手段或研发新 装备来解决边坡模型振动台试验中这一问题 ; ③如何在边坡物理模型中构建真实复杂的节理, 如何模拟控制高寒边坡模型的含水量和含冰量也 是边坡振动台物理模拟需要研究的难点之一. 1.3 高海拔寒区岩体结构数值模拟研究现状 高海拔寒区岩体工程破坏是由各种化学和力 学机制多场耦合作用造成的,影响高海拔寒区岩 体的因素众多,边坡稳定性研究非常复杂,但是边 坡失稳的根源是水−热−力多场耦合作用[37] ,寒区 金属矿山开采过程中冻融温度场、应力场、渗流 场耦合作用影响岩质边坡稳定性,与常温多场耦 合不同的是寒区岩体内会出现水冰相变,寒区岩 体工程岩−水−冰与多场耦合关系如图 4 所示. Huang 等[39−41] 根据能量守恒定律、质量守恒 定律并考虑水/冰相变的静力平衡原理,推求低温 冻结岩体热−水−力 (thermo-hydro-mechanical 简写 为 THM)完全耦合控制方程,并应用于寒区隧道的 稳定性分析. 刘乃飞[42] 建立了高海拔寒区裂隙岩 体的应力−渗流−温度−化学四场耦合数学模型,基 于自行开发的 4G2017 岩土工程仿真软件建立典 型寒区岩体工程模型. 李国峰等[38] 推导寒区岩体 含相变的温度场−渗流场−应力场三场耦合简化算 法,通过 FLAC3D 对室内岩样力学试验仿真验证了 该算法的正确性. Shen 等[43] 对准砂岩试样进行了 从 20 ℃ 到−20 ℃ 冻融循环试验,建立了考虑潜热 效应的表观热容模型(AHCM),并将该模型嵌入 COMSOL 程序中,模拟岩体内部温度场的变化规 律. 李滨锷[44] 建立非饱和土固−液−气三相耦合模 型,揭示了非饱和边坡渗流、气体运移对边坡稳定 系数的影响. 陈益峰等[45] 建立了岩土固−液−气三 相 THM 完全耦合模型. 班改革和戴剑勇[46] 建立 了非饱和岩体基于固−液−气三相系统守恒定律的 热−流−固完全耦合数学模型. 从以上研究成果来看,国内外学者针对寒区 岩体建立了多个耦合模型,但目前关于高寒岩体 的多场耦合模型多是在冻土力学的理论基础上, 将岩体看作各向同性的孔隙介质,对于寒区露天 矿山岩体各向异性、裂隙内水冰相变、流−固−气 多相多场耦合模型的研究较少,边坡多场多相耦 合致灾过程方面的研究尚不完善. 1.4 高海拔寒区岩质边坡变形破坏原位监测研究 现状 为了更好地研究高寒边坡的变形破坏机理,国 内外专家学者对寒区边坡使用多种手段进行了监测. 渗透压力 冻胀力 自重应力 机械能 水热迁移 正温 负温 岩 温度场 渗流场 应力场 水 冰 水冰相变 冰水相变 体应变 渗透压力 孔隙压力 附加应力 孔隙压力 热胀冷缩 渗透压力 渗透系数降低 热应变 渗透系数变化 图 4 岩−水−冰系统与多场耦合关系[38] Fig.4 Relationship of coupled rock−water−ice system[38] · 1378 · 工程科学学报,第 41 卷,第 11 期