Analyses on failure modes and effectiveness of the prevention measures of Jianchuandong dangerous rock mass in the Three Gorges Reservoir area
-
摘要: 本文以三峡库区箭穿洞危岩体为例,对涉水厚层危岩体的变形破坏模式和防护措施进行了研究。根据现场调查和长期监测数据可知,干湿循环作用下基座岩体的劣化是加速箭穿洞危岩体变形破坏的主导因素,并判定其破坏模式为基座滑移式崩塌。在此基础上,将危岩体的防护治理定为两部分,分别是基座软弱岩体的补强加固,以及中上部危岩体的锚索加固。通过数值模拟对防护前后危岩体的位移场以及应力场进行了分析,结果表明危岩体上部的锚索加固能够有效控制岩体的变形,基座补强能够有效控制危岩体的最大剪应力,综合防护可以显著提高箭穿洞危岩体的稳定性。该防护措施的理念及方法,可以为库区涉水危岩体的治理提供重要的参考价值。Abstract: In this paper, the Jianchuandong dangerous rock mass in the Three Gorges Reservoir area is taken as an example to study the failure mode and protection of submerged thick-layer dangerous rock mass. Based on field investigation and the monitoring data, it was found that the degradation of the base rock mass under the dry-wet cycles is the main factor to accelerate the deformation of the dangerous rock mass. Its failure mode is determined as slip collapse. According to the deformation and failure characteristics of the dangerous rock mass, the corresponding treatment measure is divided into two parts, inculding the reinforcement of the weak base rock mass, and anchor cable in the upper part of the dangerous rock mass. By comparing the displacement field and the stress field before and after protection, it is found that the prestressed reinforcement of the upper dangerous rock mass can effectively control the deformation of the rock mass, and the reinforcement of base rock mass can effectively control the maximum shear stress of the dangerous rock mass. The comprehensive protection effect is significant to ensure the stability of dangerous rock mass. The protection method can provide important reference for the treatment of submerged dangerous rock mass in the reservoir area.
-
0. 引言
对于内部不存在潜在滑移面和控制性结构面的厚层危岩体,其破坏机理十分复杂[1-2]。这种危岩体具有分布区域广、发生频率高、突发性强、破坏范围大等特点,是一类典型的山区地质灾害,对人民生命财产安全和城镇建设造成严重威胁[3-6]。对于涉水厚层危岩体,除了崩塌的直接危害以外,产生的涌浪次生灾害将进一步扩大威胁范围[7-9]。
自2008年三峡库区正式蓄水以来,由于三峡库区水位的周期性涨落,在库岸形成了高达30 m的劣化带[10-12]。库区涉水危岩体的基座部分位于劣化带区域,基座岩体长期处于上部岩体自重下,并在干湿循环作用下逐渐劣化,导致其变形破坏机理更为复杂,进一步加大了危岩体的防治难度[13-15]。
文章在现场调查、监测数据以及室内试验的基础上,分析了三峡库区箭穿洞危岩体的变形破坏特征,并对其破坏模式进行了判定。根据其变形破坏特征,提出了危岩体的治理方案,并采用数值模拟对治理方案进行了定量评价。该危岩体的防护理念对于岩质库岸的防治具有重要的参考价值。
1. 危岩体概况
箭穿洞危岩体位于重庆市巫山县望霞村。危岩体的上游侧边界为纵张裂缝(LF1: 150~226 m);下游侧边界裂缝在陡崖面上清晰可见,上宽下窄,充填或局部充填碎石土或溶蚀、残积碎石土,并向下逐渐收敛至153 m高程尖灭(LF2)。危岩体后缘边界为卸荷裂缝(LF3)张开度可达3.15 m,裂隙在高程226 m以下底部均被碎石所填充。箭穿洞危岩体的正视图见图1。
箭穿洞危岩体的三维切割边界清楚,其几何形态呈不规则的六面体,后缘高程为278~305 m,基座高程为155 m,平均高差为135 m,危岩体平均横宽约55 m,平均厚度约50 m,危岩体体积约36×104 m3,主崩方向为260°。箭穿洞危岩主要为三叠系下统大冶组第四段(T1d4)、高程280 m以上为嘉陵江组第一段(T1j1),基座以下为大冶组第三段(T1d3)。基座岩体位置处有一平硐,其中布有压力传感器。箭穿洞危岩体的典型剖面图见图2。
2. 危岩体的变形特征
由区域地质构造可知,箭穿洞危岩体是在斜坡岩体不断卸荷,长江不断侵蚀切割、构造应力释放等条件下形成的。三维边界基本形成后,重力成为主导危岩体变形的主要因素。另外,干湿循环作用下基座岩体的持续劣化进一步加速了危岩体的变形。根据现场监测资料可知(图3),危岩体的边界裂缝及基座压力随着库水位周期性升降次数的增加而持续增大。
针对基座的泥质条带灰岩,在室内完成了30次干湿循环试验,得到了初始状态、5次、15次、20次以及30次干湿循环后基座岩体的力学参数(表1)。
表 1 泥质条带灰岩力学参数Table 1. Mechanical parameters of marlstone under dry-wet cycles in the Three Gorges Reservoir area干湿循环次数 单轴抗压强度 抗拉强度 /MPa 抗剪强度 变形参数 天然状态/MPa 饱和状态/MPa 内摩擦角/(°) 黏聚力/MPa 弹性模量/(104 MPa) 泊松比μ 0 19.07 13.24 1.10 32.6 3.36 0.405 0.30 5 18.39 12.25 0.99 32.2 3.16 0.373 0.31 15 16.96 11.10 0.89 31.8 2.91 0.351 0.32 20 16.22 10.44 0.83 30.1 2.83 0.272 0.33 30 14.98 9.67 0.79 28.4 2.76 0.238 0.35 经过30次干湿循环后,基座岩体的单轴抗压强度下降约21%~26%,随着循环次数的增加,其强度的劣化率有所下降,但尚未收敛;抗拉强度和黏聚力下降约28%,随着循环次数增加,其劣化率有所下降,趋于收敛;岩石的内摩擦角下降约17%,在15次循环后劣化率有增大趋势,表明岩体的抗剪强度将持续降低;变形参数下降约40%,变形模量趋于减小,泊松比趋于增大,并且在15次循环后劣化率有增大趋势,表明岩体质量将持续降低。
3. 危岩体的失稳破坏模式
根据危岩体的变形特征可知,箭穿洞危岩体以基座的变形破坏为主导,内部不存在潜在的剪切面或导致倾倒变形的控制性结构面。对于这种类型的危岩体,其变形破坏发展一般有2种趋势,分别是基座压裂型崩塌和基座滑移型崩塌[1-2](图4)。
基座压裂型崩塌见图4(a)。如果缓坡岩体较坚硬,基座底部岩体受压集中,会导致基座岩体和接触岩体出现压致拉裂现象;基座破坏时,大量的拉裂缝和剪裂缝会出现,导致岩体整体失稳。基座滑移型崩塌见图4(b)。如果缓坡岩体较为软弱,在上部压力作用下,基座软弱岩体可能会出现剪切破坏,上覆岩体压力将软弱基座挤出,从而发生后靠滑移式的整体破坏。
基座碎裂型崩塌和基座滑移型崩塌的判定与基座岩体强度有关,根据HUNGR等[1]提出的判定方法[1],可采用应力比值(Ns)来界定危岩的破坏模式,Ns的建议值为0.25,其判定公式如下:
(1) 式中:
——危岩体重度,此处取27.2 kN/m3;H——上部危岩体的高度,此处取135 m;
——基座的抗压强度,此处取5.775 MPa;该抗压强度根据规范对室内试验数据进行了折减[16-17]。计算可知,Ns为0.63,大于0.25,因而确定箭穿洞危岩体将发生滑移破坏。
4. 危岩体的防护设计
基于箭穿洞危岩体的破坏模式,将危岩体的治理定为两部分,分别是基座软弱岩体的补强加固,以及危岩体中上部的锚索加固(图5)。其中,基座岩体的补强是为了阻断危岩体的滑移剪出;中上部锚索加固是为了控制危岩体的变形。防护治理所涉及的力学改善措施如下:
(1)基座软弱岩体补强加固工程
基座平硐采用C30钢筋混凝土键体充填支撑;基底设置3排锚桩,锚桩间距为1.75 m、2.25 m,锚桩孔径为150 mm,锚固段长度为6.00 m,基座涉水岩体的表面采用板肋式锚杆挡墙。
(2)防护工程(锚索、主动防护网、被动防护网、水下柔性防护垫)
在危岩体中上部布置6排2 000 kN级锚索,水平夹角为15°,水平及竖向间距均为6.00 m,锚索为16 φs15.2 mm,锚固段总长度为17.00 m(3.00,3.00,2.50,3.00,3.00,2.50 m分6段设置)。
5. 防护效果分析
针对防护方案,将提高岩体稳定性的防护措施进行简化后,进行有限元数值计算(所采用数值软件为MIDAS GTS),涵盖上部危岩的预应力锚索、消落带区域砂浆锚杆、板肋式锚杆挡墙及平硐充填。未进行防护加固时,平硐区域作隧洞处理;防护加固后,平硐区域采用C30钢筋混凝土的强度参数(参考值)进行分析。此外,砂浆锚杆、预应力锚杆及板肋式锚杆挡墙相关参数均为参照值[18-19],数值分析过程中的计算参数见表2。以初始状态下的危岩防护为例,对防护措施的有效性进行评价。根据相关规范要求[20-21],泥质条带灰岩的岩体黏聚力由岩石黏聚力乘以折减系数,取0.20;岩体内摩擦角由岩石内摩擦角乘以折减系数,取0.80;岩体变形参数由岩石变形参数乘以折减系数,取0.70。数值计算时,对数值模型右侧边界和左侧边界的水平方向进行了约束,底部边界采用固定约束,将危岩体的自重设定为诱发失稳的关键因素。
表 2 有限元数值计算岩体参数Table 2. Mechanical parameters of the marlstone used in the numerical simulation岩性 本构模型 弹性模量/MPa 泊松比 容重/(kN·m−3) 黏聚力/MPa 内摩擦角/(°) 灰岩(基岩) 摩尔库伦 47800 0.26 27.20 5.21 44.4 灰岩(消落带) 摩尔库伦 42000 0.24 24.50 4.82 40.2 泥质条带灰岩(消落带) 摩尔库伦 27200 0.33 26.50 1.79 37.6 泥质条带灰岩(水上基岩) 摩尔库伦 40500 0.30 26.60 2.36 37.6 水上灰岩(基岩) 摩尔库伦 50400 0.28 27.10 5.48 44.4 平硐(充填) 摩尔库伦 27000 0.20 24.20 3.18 48.6 砂浆锚杆 弹性 196000 0.28 78.50 − − 预应力锚杆 弹性 196000 0.28 78.50 − − 板肋式锚杆挡墙 弹性 27000 0.25 23.00 − − 通过对箭穿洞危岩典型剖面的有限元计算,得到该剖面加固前后的位移云图见图6。分析可知,上部岩体的锚索加固是控制危岩体变形的关键措施。危岩体的最大位移位于后缘部分,这是因为危岩体形状不规则,存在偏压应力,导致其具有沿基座滑移的变形趋势,与之前的破坏机制分析相符。在防护加固前,危岩体的最大位移为0.3235 m,综合防护加固后其最大位移为0.1313 m,降低了59.41%,危岩体的变形趋势得以控制。
箭穿洞剖面最大剪应力云图见图7,分析可知,岩体基座加固是控制剪切应力集中的关键措施。防护加固前,剪应力的最大值为38.085 MPa,且在裂隙尖端出现应力集中现象。综合防护加固后,裂隙尖端的最大剪应力为11.117 MPa,降幅可达70.81%。
对应力场及位移场进行分析可知,预应力锚索可有效控制危岩体由于偏压而引发的变形趋势,而基座加固在保证基座岩体完整性的同时,可以有效降低基座岩体的最大剪应力。
通过强度折减法对危岩体的稳定性进行了分析(图8)。根据防护前危岩体的塑性变形可知,其破坏模式为基座滑移式破坏,与前文滑移破坏模式的判定是一致的,其塑性变形区由后缘裂缝尖端向平硐位置延伸,此时危岩体稳定系数为1.04,处于临界失稳状态。在平硐充填的基础上,进行砂浆锚杆以及格构梁的支护,提升基座岩体的完整性之后,其塑性变形区下移见图8(b),危岩体的稳定性大幅度提升,稳定系数可达1.82,其提高幅度为75%。当上部岩体采用预应力锚索进行加固时,见图8(c),可有效控制危岩体的变形,与防护前的危岩体相比,其稳定性提高了17.78%。当完成综合支护后,见图8(d),其稳定性可达2.451,与防护前的稳定性相比提高了135.67%。在综合防护下,基座补强尤其是砂浆锚杆的施工阻断了塑性变形区的连续性,危岩体塑性变形区的剪出口下移到破碎带下方的消落带区域,且上部预应力锚索控制住了危岩体的整体变形,从而大幅提升了危岩体的稳定性。
通过数值模拟可知,在综合防护之后,危岩体的剪出口将下移至145 m高程处。根据原有设计方案,在145 m水位处会设置防水层以及竖向锚杆,因而,能够在之后的劣化中进一步提升危岩体的长期稳定性。由于145 m高程处的防护并非主体设计,本文在计算时并未考虑相关防护措施。箭穿洞危岩体的防护工程已经按照文中所陈述的措施完成了施工,相应的演化趋势将在之后做进一步的研究。
6. 结论
在现场调查、实时监测以及室内试验的基础上,本文采用公式判定和数值模拟等方法对三峡库区箭穿洞危岩体的破坏机理和防护效果进行了分析和研究,得到以下结论:
(1)由于箭穿洞危岩体为内部不存在潜在滑移面和控制性结构面的涉水厚层危岩体,其变形破坏机理较为复杂。箭穿洞危岩体基座部分的软弱岩层不仅承担着上覆岩层的自重,同时在库区水位的周期性升降下持续劣化,加速了危岩体的变形破坏。
(2)通过公式判定可知,在上部压力作用下,箭穿洞危岩的基座软弱岩体可能会出现剪切破坏,上覆岩体压力将软弱基座挤出,并最终发生基座滑移式的整体破坏。
(3)针对该危岩体的变形破坏模式,提出了基座加固及上部岩体固定的防护措施,其中,上部岩体的锚索加固用于控制危岩体的变形,基座加固用于控制危岩体的滑移。
(4)数值模拟结果表明,本文涉及的综合防护措施效果显著,能够有效的控制危岩体的变形,使得危岩体的塑性变形区域下移,并最终提高危岩体的整体稳定性。
-
表 1 泥质条带灰岩力学参数
Table 1 Mechanical parameters of marlstone under dry-wet cycles in the Three Gorges Reservoir area
干湿循环次数 单轴抗压强度 抗拉强度 /MPa 抗剪强度 变形参数 天然状态/MPa 饱和状态/MPa 内摩擦角/(°) 黏聚力/MPa 弹性模量/(104 MPa) 泊松比μ 0 19.07 13.24 1.10 32.6 3.36 0.405 0.30 5 18.39 12.25 0.99 32.2 3.16 0.373 0.31 15 16.96 11.10 0.89 31.8 2.91 0.351 0.32 20 16.22 10.44 0.83 30.1 2.83 0.272 0.33 30 14.98 9.67 0.79 28.4 2.76 0.238 0.35 表 2 有限元数值计算岩体参数
Table 2 Mechanical parameters of the marlstone used in the numerical simulation
岩性 本构模型 弹性模量/MPa 泊松比 容重/(kN·m−3) 黏聚力/MPa 内摩擦角/(°) 灰岩(基岩) 摩尔库伦 47800 0.26 27.20 5.21 44.4 灰岩(消落带) 摩尔库伦 42000 0.24 24.50 4.82 40.2 泥质条带灰岩(消落带) 摩尔库伦 27200 0.33 26.50 1.79 37.6 泥质条带灰岩(水上基岩) 摩尔库伦 40500 0.30 26.60 2.36 37.6 水上灰岩(基岩) 摩尔库伦 50400 0.28 27.10 5.48 44.4 平硐(充填) 摩尔库伦 27000 0.20 24.20 3.18 48.6 砂浆锚杆 弹性 196000 0.28 78.50 − − 预应力锚杆 弹性 196000 0.28 78.50 − − 板肋式锚杆挡墙 弹性 27000 0.25 23.00 − − -
[1] HUNGR O, EVANS SG. The occurrence and classification of massive rock slope failure[J]. Felsbau—Rock and Soil Engineering,2004(22):16 − 23.
[2] HUNGR O, LEROUEIL S, PICARELLI L. The Varnes classification of landslide types, an update[J]. Landslides,2014,11(2):167 − 194. DOI: 10.1007/s10346-013-0436-y
[3] YIN Y P, SUN P, ZHU J L, et al. Research on catastrophic rock avalanche at Guanling, Guizhou, China[J]. Landslides,2011,8(4):517 − 525. DOI: 10.1007/s10346-011-0266-8
[4] XING A G, WANG G H, YIN Y P, et al. Investigation and dynamic analysis of a catastrophic rock avalanche on September 23, 1991, Zhaotong, China[J]. Landslides,2016,13(5):1035 − 1047. DOI: 10.1007/s10346-015-0617-y
[5] FENG Z, LI B, HE K. Rock collapse mechanism on high-steep slopefailure in sub-horizontal thick-bedded mountains[J]. Journal of Geomechanics,2014(20):123 − 131.
[6] FENG Z, LI B, YINY P, HE K. Rockslides on limestone cliffs withsubhorizontal bedding in the southwestern calcareous area of China[J]. NatureHazards and Earth System Sciences,2014(14):2627 − 2635.
[7] YIN Y P, HUANG B L, LIU G N, et al. Potential risk analysis on a Jianchuandong dangerous rockmass-generated impulse wave in the Three Gorges Reservoir, China[J]. Environmental Earth Sciences,2015,74(3):2595 − 2607. DOI: 10.1007/s12665-015-4278-x
[8] YIN Y P, HUANG B L, WANG W P, etal. Reservoir-induced landslides and risk control in Three Gorges Project on Yangtze River, China[J]. Journal of Rock Mechanics and Geotechnical Engineering,2016,8(5):577 − 595. DOI: 10.1016/j.jrmge.2016.08.001
[9] 赵海林, 黄波林, 张全, 等. 柱状危岩体崩塌产生涌浪的物理试验与数值模拟分析[J]. 地质力学学报,2020,26(4):500 − 509. [ZHAO Hailin, HUANG Bolin, ZHANG Quan, et al. Physical experiment and numerical model analysis of surge caused by collapse of columnar dangerous rock mass[J]. Journal of Geomechanics,2020,26(4):500 − 509. (in Chinese with English abstract) DOI: 10.12090/j.issn.1006-6616.2020.26.04.043 [10] 邓华锋, 肖瑶, 方景成, 等. 干湿循环作用下岸坡消落带土体抗剪强度劣化规律及其对岸坡稳定性影响研究[J]. 岩土力学,2017,38(9):2629 − 2638. [DENG Huafeng, XIAO Yao, FANG Jingcheng, et al. Shear strength degradation and slope stability of soils at hydro-fluctuation belt of river bank slope during drying-wetting cycle[J]. Rock and Soil Mechanics,2017,38(9):2629 − 2638. (in Chinese with English abstract) [11] WANG X G, WANG J D, GU T F, et al. A modified Hoek-Brown failure criterion considering the damage to reservoir bank slope rocks under water saturation-dehydration circulation[J]. Journal of Mountain Science,2017,14(4):771 − 781. DOI: 10.1007/s11629-016-4206-x
[12] 余姝, 张枝华, 黄波林, 等. 三峡库区青石-抱龙段顺层灰岩库岸坡变形破坏机理[J]. 中国地质灾害与防治学报,2019,30(3):18 − 23. [YU Shu, ZHANG Zhihua, HUANG Bolin, et al. Deformation mechanism of bedding limestone bank slope from Qingshi to Baolong, the Three Gorges Area[J]. The Chinese Journal of Geological Hazard and Control,2019,30(3):18 − 23. (in Chinese with English abstract) [13] 唐红梅. 群发性崩塌灾害形成机制与减灾技术[D]. 重庆: 重庆大学, 2011. TANG Hongmei. Developing mechanism and mitigation techniques of massive collapse[D]. Chongqing: Chongqing University, 2011. (in Chinese with English abstract)
[14] 张枝华, 杜春兰, 余姝, 等. 三峡库区巫峡箭穿洞危岩体稳定性分析及防治工程设计[J]. 中国地质灾害与防治学报,2018,29(2):48 − 54. [ZHANG Zhihua, DU Chunlan, YU Shu, et al. Stability analysis and design of control works on Jianchuandong dangerous rockmass in Wuxia Gorge, the Three Gorges Reservoir[J]. The Chinese Journal of Geological Hazard and Control,2018,29(2):48 − 54. (in Chinese with English abstract) [15] 陈小婷, 王健, 黄波林, 等. 库水位变动条件下柱状危岩体变形破坏机理[J]. 中国地质灾害与防治学报,2019,30(2):9 − 18. [CHEN Xiaoting, WANG Jian, HUANG Bolin, et al. Deformation and failure mechanism of pillar-shaped dangerous rock mass under resenoir's water level fluctuation[J]. The Chinese Journal of Geological Hazard and Control,2019,30(2):9 − 18. (in Chinese with English abstract) [16] 贺凯. 塔柱状岩体崩塌机理研究[D]. 西安: 长安大学, 2015. HE Kai. Research on collapse mechanism of tower rock[D]. Xi'an: Chang'an University, 2015. (in Chinese with English abstract)
[17] HUANG B L, YIN Y P, ZHANG Z H, etal. A case study of pillar-shaped rock mass failure in the Three Gorges Reservoir Area, China[J]. Quarterly Journal of Engineering Geology and Hydrogeology,2016,49(3):195 − 202. DOI: 10.1144/qjegh2015-034
[18] ULUSAY R. The isrm suggested methods for rock characterization, testing and monitoring: 2007-2014 [M]. Switzerland: Springer International Publishing, 2014.
[19] 闫玉平, 肖世国. 双排抗滑桩加固滑坡的前桩后侧推力算法[J]. 水文地质工程地质,2021,48(4):55 − 63. [YAN Yuping, XIAO Shiguo. A calculation method for thrust on the fore piles of double-row stabilizing piles used to reinforce landslides[J]. Hydrogeology & Engineering Geology,2021,48(4):55 − 63. (in Chinese with English abstract) [20] 三峡库区地质灾害防治工作指挥部. 三峡库区三期地质灾害防治工程地质勘查技术要求[S]. 2012. Headquarters for Geological Hazard Prevention of the Three Gorges Reservoir Region. The technicalrequirements of the 3rd phase geological exploration inthe Three Gorges Reservoir Region[S]. 2012. (in Chinese)
[21] 闫国强, 黄波林, 代贞伟, 等. 三峡库区巫峡段典型岩体劣化特征研究[J]. 水文地质工程地质,2020,47(4):62 − 72. [YAN Guoqiang, HUANG Bolin, DAI Zhenwei, etal. A study of the deterioration effect of limestone bank slope rock mass at the Wuxia section of the Three Gorges Reservoir area[J]. Hydrogeology & Engineering Geology,2020,47(4):62 − 72. (in Chinese with English abstract) -
期刊类型引用(11)
1. 杨皓,魏涛. 三峡库区米仓口危岩体稳定性分析. 科技与创新. 2025(09): 121-124 . 百度学术
2. 殷跃平. 新三峡库区长期地质安全战略研究. 中国水利. 2024(22): 26-35 . 百度学术
3. 李伟,董远峰,邓玖林,靳鹏,高玮阳,李海洋. 基于两期机载LiDAR数据的危岩变形识别方法研究. 人民长江. 2024(S2): 121-124 . 百度学术
4. 柴乃杰,周文梁. 基于优化组合权-模糊可变集的坝基岩体质量分级. 吉林大学学报(地球科学版). 2023(02): 514-525 . 百度学术
5. 张燕,王庆兵,邢文超,王元新,张君,于桑,张建芝,葛江琨. 新疆某山区公路边坡危岩体影响区划分及防治建议. 安全与环境工程. 2023(04): 149-158 . 百度学术
6. 江俊杰,刘东泽,卢应发. 库水位降落与降雨耦合作用下鸡脑壳包滑坡变形分析. 中国农村水利水电. 2023(09): 236-243 . 百度学术
7. 蒋明,李伟,黎景. 西部某水电站枢纽区边坡危岩体防治设计研究. 小水电. 2023(05): 43-47 . 百度学术
8. 檀梦皎,殷坤龙,付智勇,朱春芳,陶小虎,朱延辉. 降雨及库水位影响下麻地湾滑坡地下水响应特征分析. 中国地质灾害与防治学报. 2022(01): 45-57 . 本站查看
9. 陈佳亮. 引水隧洞进口上部危岩体稳定性研究. 水利科学与寒区工程. 2022(05): 70-73 . 百度学术
10. 杨光明,罗垚,张帆,陈也. 三峡库区生态环境三方协同治理演化博弈及系统仿真研究. 重庆理工大学学报(社会科学). 2021(12): 154-166 . 百度学术
11. 庄明水. 厦门岛内孤(滚)石破坏模式及分布规律研究. 地质灾害与环境保护. 2021(04): 34-38+44 . 百度学术
其他类型引用(5)