ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

香丽高速公路边坡地质灾害发育特征与易发性区划

廖小平, 徐风光, 蔡旭东, 周文皎, 魏家旭

廖小平, 徐风光, 蔡旭东, 周文皎, 魏家旭. 香丽高速公路边坡地质灾害发育特征与易发性区划[J]. 中国地质灾害与防治学报, 2021, 32(5): 121-129. DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-15
引用本文: 廖小平, 徐风光, 蔡旭东, 周文皎, 魏家旭. 香丽高速公路边坡地质灾害发育特征与易发性区划[J]. 中国地质灾害与防治学报, 2021, 32(5): 121-129. DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-15
Xiaoping LIAO, Fengguang XU, Xudong CAI, Wenjiao ZHOU, Jiaxu WEI. Development characteristics and susceptibality zoning of slope geological hazards in Xiangli expressway[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(5): 121-129. DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-15
Citation: Xiaoping LIAO, Fengguang XU, Xudong CAI, Wenjiao ZHOU, Jiaxu WEI. Development characteristics and susceptibality zoning of slope geological hazards in Xiangli expressway[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32(5): 121-129. DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-15

香丽高速公路边坡地质灾害发育特征与易发性区划

基金项目: 红层地区线性工程穿越切割的地质体失稳机理和成灾模式项目(2018YFC1504901);香丽高速公路复杂边坡灾变特征及防治新技术研究项目(云交科教2018-35)
详细信息
    作者简介:

    廖小平(1965-),男,研究员,博士生导师,主要从事边坡防护加固工程、滑坡灾害整治工程的勘察设计、咨询评估和监测预警技术等方面的的科学研究与工程实践工作。E-mail:xiaoping.liao@126.com

  • 中图分类号: P694

Development characteristics and susceptibality zoning of slope geological hazards in Xiangli expressway

  • 摘要: 针对香丽高速公路边坡地质灾害,在详细的道路工程勘察设计文件的基础上,结合现场踏勘调查,系统地研究了其边坡地质灾害的主要类型、发育特征和分布规律;提出高速公路等线性工程边坡灾害的基本地质条件、自然诱发因素、人类工程活动和灾害历史记录等完整信息评价指标。采用层次分析法、专家评分法与因素累积法相结合的研究方法,建立了线性工程边坡地质灾害易发性评价模型,并基于GIS平台完成了香丽高速公路边坡地质灾害的易发性区划。为指导香丽高速公路地质灾害的危险性评价及科学防治提供了重要依据,对类似山区道路等线性工程边坡地质灾害危险性评价研究与实践具有一定的示范作用和参考价值。
    Abstract: Focus on theslope geological hazard of Xiangli expressway, systematic study was carried out to find the main types, development characteristics and distribution of slope geological hazards on the basis of detailed road engineering survey.A complete information evaluation index of linear engineering of expressway, including basic geological conditions, natural induced factors, human engineering activities and historical records of slope disasters, was put forward. The susceptibility evaluation model of slope geological hazards of linear engineering was established through the combination of Analytic Hierarchy Process (AHP), expert investigation method and factor accumulation method. The susceptibility zoning of the slope geological hazards in Xiangli expressway was completed, in terms of the GIS platform, which provided an important basis for the risk evaluation and scientific prevention of geological hazards in Xiangli expressway.This paper plays a certain demonstration role and provides reference for the research and practice of slope geological hazard assessment of linear engineering such as mountainious roads.
  • 香丽高速公路起于香格里拉县城,止于大丽高速公路白汉场立交,全长140 km。由于沿线工程地质条件复杂,生态环境脆弱,在工程建设过程中,受人工开挖卸荷、大气降雨等影响,边坡地质灾害将严重危害和威胁工程建设及运营安全,开展边坡地质灾害区划研究与应用具有十分重要的现实意义和学术价值。

    目前,国内外关于地质灾害的区划研究已经较为普遍,最初大多是基于区域性的地质灾害区划研究,后来开始对重点公路区段内的边坡地质灾害风险评估进行研究。康钦勇以永川市实际地质灾害为依托,运用层次分析法,建立了永川市地质灾害易发程度分区评价指标体系[1]。宋世鑫[2]以陕西省定边县地质灾害工点为例,运用层次分析法建立易发性评价体系,并基于GIS系统,绘制了定边县区域内的地质灾害易发性分区图[2]。殷坤龙等[3]建立了滑坡灾害分区的信息分析系统,并结合重庆市的滑坡灾害进行了区划研究。聂忠权等[4]以大连市地质灾害工点为例,建立了一套基于GIS的地质灾害易发程度分区评价模型。1986年交通部颁布《公路自然区划标准》(JTJ 003-86),进行了一级和二级区划系统的划分,并编制了全国范围的区划图[5]。林灿阳等[6]基于山区高速公路边坡工程特点,从边坡易发性、危险性、易损性三个方面进行研究,并建立了相应的评价体系。张建伟等[7]根据大量运用高速公路边坡失稳特点,建立了一套运用高速公路边坡安全评价体系。彭小平等[8]以贵州省公路地质灾害为依托,利用GIS技术,进行了地质灾害危险性区划研究。齐洪亮等[9]采用因素叠加法,基于GIS软件,进行了中国公路地质灾害危险性区划。韦威[10]根据中国公路地质灾害的分布情况,对地质灾害影响因素进行了研究,采用专家打分法建立了分区综合评价模型。综上所述,目前对山区道路等线性工程的边坡地质灾害区划研究相对薄弱和滞后。

    香丽高速公路位于青藏高原南东延横断山脉中段,属高海拔、高烈度和地质灾害高发区。沿线地貌主要有高原盆地地貌、构造剥蚀低中山地貌、河谷侵蚀深切高山峡谷地貌、构造侵蚀剥蚀中高山地貌。河谷多为“V”型,谷坡较陡,多为30°~50°。

    沿线地层主要为第四系全新统坡残积层、冲洪积层、堆积层,更新统冰积层,三叠系中上统,二叠系中上统,泥盆系中上统,石炭系上统,等等。岩性以灰岩、板岩为主,局部揭露玄武岩和凝灰岩。

    场区构造复杂,断裂发育,新构造运动十分强烈,是我国大陆现今地壳构造运动最为强烈的地区。区内构造带为“歹”字型构造体系中的网状构造带和经向构造体系中的川滇南北向构造带,主要活动断裂带为大具-丽江断裂、丽江-剑川断裂、剑川-乔后断裂、中甸-龙蟠-乔后断裂、中甸断裂(含冲江河断裂、尼西断裂)等。

    香丽高速公路地处温带和寒温带季风气候区,年平均气温5.4 ℃,年降雨量618~960 mm,雨季大部分集中在5—10月,雨量充沛。

    公路沿线河流密布,澜沧江、怒江、金沙江三江贯穿并流。以金沙江水系为主,金沙江一、二级支流水系较多,多成树枝状发育分布,总体上场区地表和地下水丰富。

    香丽高速公路沿线工程地质条件复杂,崩塌和滑坡等不良地质发育,第四系覆盖堆积厚度大,自然山体岩石破碎高陡,在降雨、地震、人类工程活动等作用下,沿线边坡地质灾害十分突出。根据详细的道路工程勘察、设计和施工文件等资料,结合工点现场调查及试验监测成果,统计共有51处较为严重的边坡地质灾害。其边坡地质灾害主要有滑坡和崩塌等两个基本类型,以及潜在不稳定斜坡变形体。由于不同成因、性质和规模的边坡工程地质灾害可能带来不同的危害程度和灾害后果,也必将采取不同的防护加固工程措施或者不同的整治工程对策。因此,本文又将边坡地质灾害细分为山体滑坡、边坡滑坡、崩塌坍塌、落石滚石、变形开裂等亚次类型(图1)。其中山体滑坡10处,边坡滑坡12处,崩塌坍塌17处,落石滚石5处,高陡斜坡变形开裂7处。

    图  1  香丽高速公路边坡地质灾害类型
    Figure  1.  Types of slope geological hazards in Xiangli expressway

    (1)山体滑坡

    山体滑坡占全线边坡地质灾害19.6%,主要为沿不同时期或不同成因的堆积界面产生滑动的堆积层滑坡和在岩堆、错落和古老滑坡等堆积体的基础上产生滑动的不良地质堆积体滑坡两大类。

    堆积层滑坡4处,自然斜坡多为15°~35°,较平缓,滑坡上部土体主要为含砾黏土、碎石土,土体松散、空隙大,遇水易软化,下伏三叠系中统板岩。受开挖卸荷和强降雨入渗影响,上部土体沿不同时期或不同成因的堆积界面产生滑动,坡脚剪出,后缘拉裂下错。

    不良地质堆积体滑坡6处,受桥梁、隧道、路基等工程影响,在不良地质堆积体中前部穿越切割和强降雨入渗的作用下,沿堆积界面或老滑动面产生滑动,引起不良地质堆积体复活变形和破坏。典型工点有洼里别隧道进口滑坡、昌格洛隧道出口巨型滑坡、K92老滑坡等。

    (2)边坡滑坡

    边坡滑坡占全线边坡地质灾害的23.5%,主要指公路边坡范围内岩土体沿一定滑动面产生单级或多级滑动变形破坏。沿线边坡滑坡主要为沿顺倾土岩界面滑动的二元结构滑坡和沿沉积岩层面、变质岩似层面等产生顺层滑动的顺倾层状滑坡。

    二元结构滑坡主要特征为边坡上覆土层、下伏基岩、土岩界面顺倾,在工程开挖卸荷和降雨入渗的作用下,边坡上部土体沿土岩界面产生单级或多级边坡变形破坏。顺倾层状滑坡主要表现为边坡岩土体受风化和构造作用,岩土体风化强烈,岩体破碎,节理裂隙发育,受开挖卸荷作用和降雨影响,沿顺倾的灰岩层面和板岩似层面等结构面产生顺层滑坡。

    (3)崩塌(坍塌)

    崩塌(坍塌)占全线边坡地质灾害33.3%,为全线边坡地质灾害发生数量最多。沿线自然斜坡较陡,多为陡倾直线坡,土体松散,岩体破碎,节理裂隙发育。破碎岩石边坡在自然或人为卸荷作用下,因降雨入渗软化结构面,使其沿结构面滑移崩解,堆积在坡脚形成松散堆积体,即崩塌灾害;另一部分边坡在坡面一定松弛范围内,受降雨入渗影响,坡面发生逐步由外至内的坍塌变形,坍塌后坡面上部较陡,下部较缓,略带内凹弧形。

    (4)落石(滚石)

    落石(滚石)占全线边坡地质灾害9.8%,相对较少。一方面主要发育在桥梁修筑过程中,桥墩位于自然陡坡附近,为保护桥墩,在山侧开挖边坡时,引起边坡附近的孤石滚落;另一方面,在山谷、河谷附近的自然坡体受爆破开挖施工影响,引起边坡一定范围内岩体松弛破碎和坡顶岩体松动,使岩石脱离母体沿坡面滚落至坡脚,但滚落岩石尺寸较小。

    (5)变形开裂

    变形开裂占全线边坡地质灾害13.8%,主要发育在沿线桥梁墩台边缘高陡边坡和隧道洞口侧坡、仰坡,此类边坡高陡,岩体破碎,节理裂隙发育,在人类工程活动和降雨入渗的作用下,边坡顶部附近发生局部卸荷变形开裂。路基沉降等引起的路堤边坡变形开裂病害较少发生。

    (1)空间分布

    香丽高速公路边坡地质灾害沿线分布相对集中,主要分布在K15~K20、K35~K47、K83~K108等河谷两侧构造剥蚀缓倾、陡倾直线坡地段和洼里别隧道—金沙江开达古隧道河谷深切高山峡谷地段。其中,K15~K20、K35~K47段主要为山体滑坡、边坡滑坡、崩塌坍塌;K83~K108段主要为落石滚石、崩塌坍塌及不良地质体复活引起的山体滑坡;洼里别隧道—金沙江开达古隧道高山峡谷地段主要为不良地质体复活引起的山体滑坡和隧道洞口侧坡、仰坡的变形开裂。

    (2)时间分布

    根据现场调查和边坡施工资料,沿线边坡地质灾害发生时间主要集中在雨季(6—8月)。受连续强降雨影响,沿线边坡松散土体和破碎岩体迅速从非饱和状态转为饱和状态,使坡体内岩土体抗剪强度大幅降低,极易诱发崩塌坍塌、边坡滑坡、不良地质体复活等边坡地质灾害。

    文章提出的边坡地质灾害易发性评价是一种多指标综合评价方法,即根据边坡地质灾害类型、发育特征和分布规律,提出易发性评价指标,并将评价指标从定性到定量转化,再结合层次分析法、专家评分法和因素累积法,建立易发性评价模型。

    高速公路等线性工程跨越了不同的地貌单元和地质环境,边坡地质灾害具有分布广、种类多、成因复杂的特点,是多种内在因素和外在因素共同作用导致的,影响易发性评价因素多而复杂。为提高边坡地质灾害易发性评价的科学性,减少主观性和随意性,本文在选取评价指标时,遵循系统性、主次性、独立性和可操作性原则,从基本工程地质条件、自然诱发因素、人类工程活动和灾害历史等因素综合分析,提出基本地质条件为主控因素,大气降雨、地震作用、人类工程活动和灾害历史为附加因素,并选择了11个评价因子(图2)。

    图  2  边坡地质灾害易发性评价因子树形结构图
    Figure  2.  Tree diagrams of susceptibility evaluation factors for the slope geological hazards

    根据前述评价因子,结合香丽高速公路沿线工程地质特点和山体滑坡、边坡滑坡、崩塌坍塌、落石滚石、变形开裂等5种边坡地质灾害特征,充分考虑评价因子的可操作性,参考国内外既有研究成果,基于专家评分法,将评价因子进行细分,并按1~9的标度进行量化,标度越大,易发性影响越大。

    (1)基本地质条件

    ①地形地貌

    按重力堆积和侵蚀剥蚀作用,将沿线地形地貌特征分为18种,并按不同边坡灾害类型对其量化,见表1表2

    表  1  重力堆积地形地貌易发性评价分值表
    Table  1.  Susceptibility evaluation scores of gravity accumulation topography and geomorphology
    重力堆积地形
    地貌亚类细分
    不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    坡积裙34511
    洪积扇56711
    坡积台地87511
    岩屑坡45696
    岩堆56797
    滑坡99911
    泥石流67811
    错落体67853
    下载: 导出CSV 
    | 显示表格
    表  2  侵蚀剥蚀地形地貌易发性评价分值表
    Table  2.  Susceptibility evaluation scores of erosion topography and geomorphology
    侵蚀剥蚀地形
    地貌亚类细分
    不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    平坦直线坡11111
    缓倾直线坡11454
    陡倾直线坡14797
    凸形坡12452
    凹形坡47976
    台阶形坡97411
    山顶12321
    鞍部43344
    洼地、谷地11111
    陡崖47997
    下载: 导出CSV 
    | 显示表格

    ②地层岩性

    地层岩性按其坚硬程度综合考虑分为土层、软岩、硬岩等3类,见表3

    表  3  地层岩性评价分值
    Table  3.  Evaluation scores of formation lithology
    地层岩性
    坚硬程度细分
    不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    土层99919
    软岩66797
    硬岩22352
    下载: 导出CSV 
    | 显示表格

    ③岩土体结构

    岩土体结构主要考虑节理面、土岩界面、岩层层面、断层或软弱面等4类,按节理面发育程度和土岩界面、岩层层面、断层或软弱面与坡面的相互关系进行量化,见表4表5表6表7

    表  4  节理发育程度评价分值
    Table  4.  Evaluation scores of joint growth level
    节理发育程度不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    0条/10 m00000
    1~2条/10 m22344
    3~4条/10 m44677
    ≥5条/10 m77899
    下载: 导出CSV 
    | 显示表格
    表  5  土岩界面评价分值
    Table  5.  Evaluation scores of soil-rock interfaces
    土岩界面不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    顺倾(0°~10°)54311
    顺倾(10°~35°)98711
    顺倾(35°~60°)78911
    顺倾(60°~90°)47911
    反倾(0°~90°)11111
    下载: 导出CSV 
    | 显示表格
    表  6  岩层层面评价分值
    Table  6.  Evaluation scores of rock interfaces
    岩层层面不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    顺倾(0°~10°)54331
    顺倾(10°~35°)98772
    顺倾(35°~60°)78994
    顺倾(60°~90°)46888
    反倾(60°~90°)57999
    反倾(35°~60°)24888
    反倾(10°~35°)11114
    反倾(0°~10°)11222
    下载: 导出CSV 
    | 显示表格
    表  7  断层或软弱面评价分值
    Table  7.  Evaluation score of fault or weak surface
    断层或软弱面不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    顺倾(0°~10°)54331
    顺倾(10°~35°)98772
    顺倾(35°~60°)78994
    顺倾(60°~90°)47988
    反倾(0°~90°)11111
    下载: 导出CSV 
    | 显示表格

    ④地层结构

    边坡地层结构是指组成边坡地层的结构特征与组合规律。基于地勘钻孔及露头揭示覆盖土层A、全至强风化岩层B、弱至微风化岩层C的相对厚度及其组合关系,本文建立了7种地层结构模式,(图3)。并根据不同灾害类型进行量化(表8)。

    图  3  边坡地层结构模式分类图
    Figure  3.  Classification of slope stratigraphic structure mode
    表  8  岩层结构评价分值
    Table  8.  Evaluation scores of rock structure
    边坡地层
    结构模式
    不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    模式199901
    模式267801
    模式334501
    模式417901
    模式556751
    模式645651
    模式713575
    下载: 导出CSV 
    | 显示表格

    ⑤地下水发育程度

    地下水发育程度评价根据不同地下水埋藏条件和发育程度分级量化,按类型取大值(表9)。

    表  9  地下水发育程度评价分值
    Table  9.  Evaluation scores of groundwater development level
    类型发育程度评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    上层滞水坡面干燥11111
    坡面潮湿57971
    潜水埋深>2 m11111
    埋深<2 m45211
    溢出泉78911
    承压水无承压水11111
    有承压水75211
    上升泉95211
    下载: 导出CSV 
    | 显示表格

    (2)自然诱发因素

    自然诱发因素主要考虑地震作用和大气降雨共2种,其中大气降雨主要考虑沿线地区最大日降雨量,地震按烈度考虑,地震对各种边坡地质灾害作用大小相同,见表10表11

    表  10  地震作用评价分值
    Table  10.  Evaluation scores of seismic effect
    地震作用(烈度)≤Ⅳ≥Ⅸ
    评价分值013579
    下载: 导出CSV 
    | 显示表格
    表  11  大气降雨评价分值
    Table  11.  Evaluation scores of rainfall
    最大日降雨量/
    mm
    评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    <500000
    5~2511111
    25~5023532
    50~10056765
    100~20067876
    >20078987
    下载: 导出CSV 
    | 显示表格

    (3)人类工程活动

    人类工程活动主要按线性工程特点进行细分,分为路基、桥隧、隧道,不同的工程对边坡地质灾害影响不同。本文将路基分为路堤和路堑,按不同填挖特征量化,桥梁主要按跨越的地貌特征量化,隧道穿越按相对位置量化(表12)。

    表  12  人类工程活动评价分值
    Table  12.  Evaluation scores of human engineering activities
    类型发育程度评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    路堤一般路堤11101
    高路堤75507
    软基路堤57709
    陡坡路堤59909
    路堑浅路堑44411
    深路堑77777
    超深路堑99979
    桥梁(墩台)一般斜坡11100
    陡坡78999
    不良地质98799
    隧道洞身11000
    洞口陡坡78999
    不良地质99999
    下载: 导出CSV 
    | 显示表格

    (4)灾害历史

    灾害历史是指高速公路等线性工程路线走廊带历史上曾发生的边坡地质灾害,按沿线每km走廊带发生的同类灾害点数计取(表13)。

    表  13  灾害点密度评价分值
    Table  13.  Evaluation scores of disaster point density
    灾害点密度(个/km)01234>4
    评价分值013579
    下载: 导出CSV 
    | 显示表格

    目前,用于边坡地质灾害易发性评价的主要方法有专家评分法、模糊评判法、神经网络法和层次分析法等。

    本文根据前述11个评价因子,充分考虑评价因子相互逻辑关系及边坡地质灾害作用机理和影响程度,避免人为性、主观性、随意性。首先将评价因子分为基本地质条件主控因素,大气降雨、地震作用、人类工程活动和灾害历史为附加因素,各评价因素原则上相对独立。其中,灾害历史因素反映或体现其它因素的共同作用或综合作用效果,作为一种重要的附加因素是可行的,也是必要的。然后按前述11个评价因子的量化分值,对5类评价因素分别进行评价分值计算。其中基本地质条件和人类工程活动这两类评价因素采用层次分析法计算,其它三类评价因素为定量指标,直接按前述专家评分法赋值。本文的易发性评价是用于边坡地质灾害的灾势区划及其防治决策依据,不考虑防治工程等因素。

    在边坡地质灾害易发性评价分析中,首先需要对各因素进行“归一化”处理,即各因素评分值分别除以最大标度9,得到各因素的相对易发指数;然后将附加因素易发指数进行“1+x”换算,即换算为附加因素易发性增大系数;最后对主控因素和附加因素采用因素乘积法运算,本文称之为“因素累积法”,从而建立灾害易发性评价模型,如式(1)。

    $$P = \frac{{{P_1}}}{9} \times \left( {1 + \frac{{{P_2}}}{9}} \right) \times \left( {1 + \frac{{{P_3}}}{9}} \right) \times \left( {1 + \frac{{{P_4}}}{9}} \right) \times \left( {1 + \frac{{{P_5}}}{9}} \right)$$ (1)

    式中:P——边坡地质灾害易发指数;

    P1——基本地质条件评价分值;

    P2——人类工程活动评价分值;

    P3——大气降雨评价分值;

    P4——地震作用评价分值;

    P5——灾害历史评价分值。

    “归一化”处理是把各因素评分值变为(0,1)之间的小数,也可以理解为是把各因素变为无量纲表达式;“1+x”换算是为了把附加因素易发性增大系数控制在1~2倍,避免因单个附加因素评价分值为0或较低值而不合理地决定易发性指标的情况。“因素累积法”是基于风险评估中常用的“乘积法”,结合本文的“归一化”处理和“1+x”换算,既反映了主控因素的决定性和有效性(0,1),又体现了附加因素的可控性和相对性(1,2),从而保障边坡地质灾害易发性评价模型的科学性、客观性和准确性。

    (1)基本地质条件评价

    基地地质条件是边坡地质灾害发生的自有内在因素,包括地形地貌、地层岩性、地层结构、岩土体结构、地下水发育程度。不同因素对边坡地质灾害影响不同,同时对不同边坡地质灾害的作用大小也不相同,为确定基本地质条件下各评价因子对边坡地质灾害的作用大小,本文采用层次分析法,构建了基本地质条件评价的层次结构模型和判断矩阵。如表14表15为山体滑坡层次分析法判断矩阵,其他灾种类同。

    表  14  山体滑坡基本地质条件层次分析法判断矩阵
    Table  14.  Judgment matrix in AHP for basic geological conditions of landslide
    基本地质条件地形地貌地层岩性地层结构岩土结构地下水
    地形地貌11/21/31/51/6
    地层岩性211/21/31/4
    地层结构3211/21/3
    岩土体结构53211/2
    地下水64321
    下载: 导出CSV 
    | 显示表格
    表  15  山体滑坡岩土体结构层次分析法判断矩阵
    Table  15.  Judgment matrix in AHP for rock-soil structure of landslide
    岩土体结构节理面土岩界面岩层层面断层软层
    节理面11/31/41/5
    土岩界面311/21/3
    岩层层面4211/2
    断层软层5321
    下载: 导出CSV 
    | 显示表格

    根据上述层次结构模型,结合香丽高速公路沿线边坡地质灾害发育特征,采用AHP软件构造了不同灾害类型的判断矩阵,并进行权重计算和一致性检验,基本地质条件数学计算模型如式(2)。

    $${P_1} = \sum\limits_{i = 1}^{10} {{\alpha _i}{A_{\rm{i}}}} $$ (2)

    式中:P1——基本地质条件评价分值;

    αi——因子权重系数;

    Ai——具体评价因子分值,见表16

    表  16  基本地质条件评价因子权重系数
    Table  16.  Weight coefficient of evaluation factor for basic geological conditions
    评价因子权重系数
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    地形地貌0.05580.05780.04240.04280.0548
    节理发育0.01980.02320.02300.09170.0718
    土岩界面0.04600.05990.07610.04060.0322
    岩层层面0.07630.10810.11810.15290.1083
    断层软层0.12640.11400.12380.16310.1430
    地层岩性0.09600.09030.08140.09820.1153
    地层结构0.15710.14460.12750.13980.1420
    地下水0.42260.40200.40780.27100.3327
    下载: 导出CSV 
    | 显示表格

    (2)人类工程活动评价

    山区道路等线性工程修筑过程中,高填深挖、半挖半填、半桥半路、桥隧相连十分普遍,不同工程类型对边坡地质灾害诱发因素不同。一部分只与单一的工程活动相关,如挖方、填方、桥梁、隧道,一部分是几种工程活动共同作用引起的。本文采用层次分析法,构建人类工程活动的层次结构模型,结合香丽高速公路沿线边坡地质灾害发育特征,采用AHP软件构造了不同灾害类型的判断矩阵,并进行权重计算和一致性检验,人类工程活动数学计算模型如式(3)。

    $${P_2} = \sum\limits_{i = 1}^4 {{\alpha _i}{A_i}} $$ (3)

    式中:P2——人类工程活动评价分值;

    αi——因子权重系数;

    Ai——具体评价因子分值,如表17

    表  17  人类工程活动评价因子权重系数
    Table  17.  Weight coefficient of evaluation factor for human engineering activity
    评价因子权重系数
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    路堤0.13470.15580.14290.12460.1078
    路堑0.40420.46740.42860.49860.4312
    桥梁0.16380.13730.14290.13730.1638
    隧道0.29730.23950.28570.23950.2973
    下载: 导出CSV 
    | 显示表格

    在上述边坡地质灾害易发性评价模型基础上,依托ArcGIS平台,对香丽高速公路边坡地质灾害进行了易发性分级区划研究。

    (1)基础资料

    根据上述建立的边坡地质灾害易发性评价数学模型,需要收集的主要基础资料有公路沿线地形图(包含线位、示坡线)、勘察设计及施工资料、现场调查、最大日降雨量数据和地震烈度等。

    (2)单元划分

    山区线性工程跨越了不同的地区、地形地貌单元、地质环境,为提高易发性区划的准确性和便捷性,根据已收集的基础资料,以公路沿线500 m长为一单元的分段原则(根据需要可以更长或更短),对公路进行单元划分。

    (3)因子评价

    在500 m长一单元划分的基础上,结合基础资料,按边坡地质灾害易发性评价模型包含的评价因子,对每个单元的基本地质条件、大气降雨、地震作用、人类工程活动、灾害历史等5类指标进行赋值。

    (4)图件处理

    将公路带状地形图(包含线位、示坡线)导入ArcGIS软件,按500 m一单元的原则在图件上进行单元离散。

    (5)易发性评价指数

    根据上述建立的山体滑坡、边坡滑坡、崩塌坍塌、落石滚石、变形开裂5类灾害易发性评价模型,首先在ArcGIS软件建立计算模型,然后将500 m一单元的评价因子分值导入ArcGIS系统进行计算,即可得每单元不同灾种的易发性评价指数。

    (6)易发性分级区划

    根据计算得到的易发性评价指数结果,采用间断点法,对不同灾害进行易发性分级区划,最后利用ArcGIS空间叠加功能,将山体滑坡、边坡滑坡、崩塌坍塌、落石滚石、变形开裂易发性评价指数结果图进行叠加,即可得到全线边坡地质灾害易发性分级区划图。

    本文按上述建立的易发性评价模型和分级区划方法,对香丽高速K0+560~K48+900段进行了边坡地质灾害区划研究。

    根据不同边坡地质灾害易发性评价指数计算结果,采用间断点法,按4级分级划分。即:(0~4)为不易发区、[4~8)为低易发区、[8~12)为中易发区、[12~16)为高易发区。

    基于GIS软件,按上述划分分值,得到了香丽高速边坡地质灾害易发性区划图,见图6(受图件大小限制,以K0+560~K7+060段为例展示)。

    图  6  香丽高速公路K0+560~K7+060段边坡地质灾害易发性区划图
    Figure  6.  Susceptibility zoning map of geological hazards at K0+560~K7+060 section of Xiangli expressway

    根据上述区划图,为验证易发性评价模型、区划方法和结果的可靠性,本文将易发性区划图和现场调查的灾害工点分布进行了对比分析,得到已发生边坡地质灾害均位于易发性区划图的中易发区和高易发区,与现场调查结果基本吻合,统计结果见表18

    表  18  香丽高速公路边坡地质灾害易发性分级区划结果
    Table  18.  Classification and regionalization results of geological hazard susceptibility on the slope of Xiangli expressway
    危险性分区不危险区低危险区中危险区高危险区
    区段长度/m202401200091007000
    长度占比/%41.924.818.814.5
    山体滑坡/处0012
    边坡滑坡/处0025
    崩塌(坍塌)/处0009
    落石(滚石)/处0000
    变形开裂/处0010
    合计/处00416
    数量占比/%002080
    下载: 导出CSV 
    | 显示表格

    根据上述香丽高速公路边坡地质灾害易发性分区结果,建立不同分级风险管理,采取不同的风险对策,能有效减少和避免边坡地质灾害的发生。

    例如,边坡灾害高易发区主要位于深和超深路堑开挖、不良地质体等段落,建议通过优化线路平纵断面设计,减少路堑边坡高度,绕避不良地质体,降低边坡灾害风险。并且提高边坡防护加固工程等级,列为重点监测预警区段。

    边坡灾害中危险区主要位于陡坡路基、路堑开挖、隧道洞口等段落,建议加强边坡防护加固工程对策设计,严格控制施工工序,确保主体工程质量,并加强边坡安全监测预警。

    边坡灾害低危险区主要位于软基和一般斜坡桥梁跨越段落,要求明确软基处理方案,控制桥梁墩台基础施工土石方开挖,及时进行边坡安全防护。

    边坡灾害不危险区主要位于一般路基和隧道洞身段落,采用普通防护工程设计,正常施工和质量控制,正常运营养护。

    总结提出了香丽高速公路沿线山体滑坡、边坡滑坡、崩塌坍塌、落石滚石和变形开裂五类地质灾害及其发育特征和时空分布规律,基于基本工程地质条件、地质灾害诱发因素和灾害历史,提出了地形地貌、岩土体结构、地层岩性、地层结构、地下水发育程度、大气降雨、地震作用、路基、桥梁、隧道、灾害历史五类11个地质灾害易发性评价指标及其量化方法。根据评价因子相互逻辑关系及其对边坡地质灾害作用和影响程度,按评价因子量化分值对五类评价因素进行评价分值计算,通过评价因素归一化处理、主控因素与附加因素的因素累积法计算,建立了香丽高速公路边坡地质灾害易发性评价模型,依托GIS平台给出了香丽高速公路边坡地质灾害易发性区划图,与香丽高速公路边坡地质灾害实际基本吻合。可为类似山区道路等线性工程提供一种实用的边坡地质灾害危险性评价方法。

  • 图  1   香丽高速公路边坡地质灾害类型

    Figure  1.   Types of slope geological hazards in Xiangli expressway

    图  2   边坡地质灾害易发性评价因子树形结构图

    Figure  2.   Tree diagrams of susceptibility evaluation factors for the slope geological hazards

    图  3   边坡地层结构模式分类图

    Figure  3.   Classification of slope stratigraphic structure mode

    图  6   香丽高速公路K0+560~K7+060段边坡地质灾害易发性区划图

    Figure  6.   Susceptibility zoning map of geological hazards at K0+560~K7+060 section of Xiangli expressway

    表  1   重力堆积地形地貌易发性评价分值表

    Table  1   Susceptibility evaluation scores of gravity accumulation topography and geomorphology

    重力堆积地形
    地貌亚类细分
    不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    坡积裙34511
    洪积扇56711
    坡积台地87511
    岩屑坡45696
    岩堆56797
    滑坡99911
    泥石流67811
    错落体67853
    下载: 导出CSV

    表  2   侵蚀剥蚀地形地貌易发性评价分值表

    Table  2   Susceptibility evaluation scores of erosion topography and geomorphology

    侵蚀剥蚀地形
    地貌亚类细分
    不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    平坦直线坡11111
    缓倾直线坡11454
    陡倾直线坡14797
    凸形坡12452
    凹形坡47976
    台阶形坡97411
    山顶12321
    鞍部43344
    洼地、谷地11111
    陡崖47997
    下载: 导出CSV

    表  3   地层岩性评价分值

    Table  3   Evaluation scores of formation lithology

    地层岩性
    坚硬程度细分
    不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    土层99919
    软岩66797
    硬岩22352
    下载: 导出CSV

    表  4   节理发育程度评价分值

    Table  4   Evaluation scores of joint growth level

    节理发育程度不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    0条/10 m00000
    1~2条/10 m22344
    3~4条/10 m44677
    ≥5条/10 m77899
    下载: 导出CSV

    表  5   土岩界面评价分值

    Table  5   Evaluation scores of soil-rock interfaces

    土岩界面不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    顺倾(0°~10°)54311
    顺倾(10°~35°)98711
    顺倾(35°~60°)78911
    顺倾(60°~90°)47911
    反倾(0°~90°)11111
    下载: 导出CSV

    表  6   岩层层面评价分值

    Table  6   Evaluation scores of rock interfaces

    岩层层面不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    顺倾(0°~10°)54331
    顺倾(10°~35°)98772
    顺倾(35°~60°)78994
    顺倾(60°~90°)46888
    反倾(60°~90°)57999
    反倾(35°~60°)24888
    反倾(10°~35°)11114
    反倾(0°~10°)11222
    下载: 导出CSV

    表  7   断层或软弱面评价分值

    Table  7   Evaluation score of fault or weak surface

    断层或软弱面不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    顺倾(0°~10°)54331
    顺倾(10°~35°)98772
    顺倾(35°~60°)78994
    顺倾(60°~90°)47988
    反倾(0°~90°)11111
    下载: 导出CSV

    表  8   岩层结构评价分值

    Table  8   Evaluation scores of rock structure

    边坡地层
    结构模式
    不同边坡灾害类型评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    模式199901
    模式267801
    模式334501
    模式417901
    模式556751
    模式645651
    模式713575
    下载: 导出CSV

    表  9   地下水发育程度评价分值

    Table  9   Evaluation scores of groundwater development level

    类型发育程度评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    上层滞水坡面干燥11111
    坡面潮湿57971
    潜水埋深>2 m11111
    埋深<2 m45211
    溢出泉78911
    承压水无承压水11111
    有承压水75211
    上升泉95211
    下载: 导出CSV

    表  10   地震作用评价分值

    Table  10   Evaluation scores of seismic effect

    地震作用(烈度)≤Ⅳ≥Ⅸ
    评价分值013579
    下载: 导出CSV

    表  11   大气降雨评价分值

    Table  11   Evaluation scores of rainfall

    最大日降雨量/
    mm
    评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    <500000
    5~2511111
    25~5023532
    50~10056765
    100~20067876
    >20078987
    下载: 导出CSV

    表  12   人类工程活动评价分值

    Table  12   Evaluation scores of human engineering activities

    类型发育程度评价分值
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    路堤一般路堤11101
    高路堤75507
    软基路堤57709
    陡坡路堤59909
    路堑浅路堑44411
    深路堑77777
    超深路堑99979
    桥梁(墩台)一般斜坡11100
    陡坡78999
    不良地质98799
    隧道洞身11000
    洞口陡坡78999
    不良地质99999
    下载: 导出CSV

    表  13   灾害点密度评价分值

    Table  13   Evaluation scores of disaster point density

    灾害点密度(个/km)01234>4
    评价分值013579
    下载: 导出CSV

    表  14   山体滑坡基本地质条件层次分析法判断矩阵

    Table  14   Judgment matrix in AHP for basic geological conditions of landslide

    基本地质条件地形地貌地层岩性地层结构岩土结构地下水
    地形地貌11/21/31/51/6
    地层岩性211/21/31/4
    地层结构3211/21/3
    岩土体结构53211/2
    地下水64321
    下载: 导出CSV

    表  15   山体滑坡岩土体结构层次分析法判断矩阵

    Table  15   Judgment matrix in AHP for rock-soil structure of landslide

    岩土体结构节理面土岩界面岩层层面断层软层
    节理面11/31/41/5
    土岩界面311/21/3
    岩层层面4211/2
    断层软层5321
    下载: 导出CSV

    表  16   基本地质条件评价因子权重系数

    Table  16   Weight coefficient of evaluation factor for basic geological conditions

    评价因子权重系数
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    地形地貌0.05580.05780.04240.04280.0548
    节理发育0.01980.02320.02300.09170.0718
    土岩界面0.04600.05990.07610.04060.0322
    岩层层面0.07630.10810.11810.15290.1083
    断层软层0.12640.11400.12380.16310.1430
    地层岩性0.09600.09030.08140.09820.1153
    地层结构0.15710.14460.12750.13980.1420
    地下水0.42260.40200.40780.27100.3327
    下载: 导出CSV

    表  17   人类工程活动评价因子权重系数

    Table  17   Weight coefficient of evaluation factor for human engineering activity

    评价因子权重系数
    山体滑坡边坡滑坡崩塌坍塌落石滚石变形开裂
    路堤0.13470.15580.14290.12460.1078
    路堑0.40420.46740.42860.49860.4312
    桥梁0.16380.13730.14290.13730.1638
    隧道0.29730.23950.28570.23950.2973
    下载: 导出CSV

    表  18   香丽高速公路边坡地质灾害易发性分级区划结果

    Table  18   Classification and regionalization results of geological hazard susceptibility on the slope of Xiangli expressway

    危险性分区不危险区低危险区中危险区高危险区
    区段长度/m202401200091007000
    长度占比/%41.924.818.814.5
    山体滑坡/处0012
    边坡滑坡/处0025
    崩塌(坍塌)/处0009
    落石(滚石)/处0000
    变形开裂/处0010
    合计/处00416
    数量占比/%002080
    下载: 导出CSV
  • [1] 康钦容. 地质灾害区划及管理信息系统研究[D]. 重庆: 重庆大学, 2007.

    KANG Qinrong. Study on zoning of geological hazard and manage information system[D]. Chongqing: Chongqing University, 2007. (in Chinese with English abstract)

    [2] 宋世鑫. 基于GIS的定边县地质灾害易发性评价[D]. 西安: 长安大学, 2015.

    SONG Shixin. The evaluation of the emergence about geological disaster in Dingbian County based on GIS[D]. Xi'an: Changan University, 2015. (in Chinese with English abstract)

    [3] 殷坤龙, 柳源. 滑坡灾害区划系统研究[J]. 中国地质灾害与防治学报,2000,11(4):28 − 32. [YIN Kunlong, LIU Yuan. Systematic studies on landslide hazard zonation[J]. The Chinese Journal of Geological Hazard and Control,2000,11(4):28 − 32. (in Chinese with English abstract) DOI: 10.3969/j.issn.1003-8035.2000.04.007
    [4] 聂忠权, 盛丽君, 范文. 基于GIS技术的地质灾害易发程度分区评价系统[J]. 公路交通科技,2005,22(增刊 1):156 − 159. [NIE Zhongquan, SHENG Lijun, FAN Wen. The distriction and assessment system of geologic disaster occurrence level based on geographic information system-taken Dalian for example[J]. Journal of Highway and Transportation Research and Development,2005,22(Sup 1):156 − 159. (in Chinese with English abstract)
    [5] 中华人民共和国交通部. 公路自然区划标准 :JTJ 003—1986[S]. 北京: 中国标准出版社, 1987.

    Ministry of Transport of the People’s Republic of China. Standard of climatic zoning for highway: JTJ 003—1986[S]. Beijing: Standards Press of China, 1987. (in Chinese)

    [6] 林灿阳, 廖小平. 高速公路边坡灾害评估与防控决策研究[J]. 路基工程,2014(5):72 − 76. [LIN Canyang, LIAO Xiaoping. Study on disaster evaluation for highway slope and the countermeasures of control[J]. Subgrade Engineering,2014(5):72 − 76. (in Chinese with English abstract)
    [7] 张建伟, 廖小平, 魏土荣. 高速公路运营期高边坡安全风险评估体系研究[J]. 路基工程,2015(6):197 − 203. [ZHANG Jianwei, LIAO Xiaoping, WEI Turong. Study on safety risk assessment system for high slope of highway in operation[J]. Subgrade Engineering,2015(6):197 − 203. (in Chinese with English abstract)
    [8] 彭小平, 陈开圣, 王成华, 等. 贵州省公路地质灾害基本特征及危险性分区[J]. 中国地质灾害与防治学报,2009,20(2):46 − 51. [PENG Xiaoping, CHEN Kaisheng, WANG Chenghua, et al. Basic characteristics and risk zoning of highway geological hazard in Guizhou Province[J]. The Chinese Journal of Geological Hazard and Control,2009,20(2):46 − 51. (in Chinese with English abstract) DOI: 10.3969/j.issn.1003-8035.2009.02.011
    [9] 齐洪亮, 尹超, 田伟平, 等. 基于ArcGIS的中国公路地质灾害危险性区划[J]. 长安大学学报(自然科学版),2015,35(5):22 − 27. [QI Hongliang, YIN Chao, TIAN Weiping, et al. Risk regionalization of highway geo-hazards in China based on ArcGIS[J]. Journal of Chang'an University (Natural Science Edition),2015,35(5):22 − 27. (in Chinese with English abstract)
    [10] 韦威. 基于GIS的公路地质灾害分区研究[D]. 西安: 长安大学, 2008.

    WEI Wei. Research on the regionalization of highway geological disasters based on GIS[D]. Xi'an: Chang'an University, 2008. (in Chinese with English abstract)

  • 期刊类型引用(21)

    1. 李志国,徐涛,刘永杰,赵立春,徐勇超,杨天鸿,郑小斌. 露天矿边坡稳定性的层次分析-模糊综合评价耦合分析. 中国地质灾害与防治学报. 2024(01): 116-123 . 本站查看
    2. 李晓玮,刘翠娜. 基于AHP法山区公路边坡稳定性评价及危害性分析. 河北地质大学学报. 2024(02): 61-65+134 . 百度学术
    3. 李晓玮. 北京北宫镇大灰厂路牵引式岩质滑坡勘查及防治对策. 城市地质. 2024(02): 139-148 . 百度学术
    4. 李勇,卢浩,黄治国,田卿燕,林森,张豹. 广东省营运高速公路边坡病害特征分析. 广东公路交通. 2024(04): 41-47 . 百度学术
    5. 王有涛. 公路灾害特征及处治设计分析. 科技资讯. 2024(17): 190-194 . 百度学术
    6. 洪勃,唐亚明,冯卫,李政国,潘学树,冯凡,周永恒,尹春旺. 基于3D实景模型和AHP边坡危险性评价系统的应用:以沿黄公路吴堡–永和段为例. 西北地质. 2024(06): 218-233 . 百度学术
    7. 支泽民,刘峰贵,周强,夏兴生,陈琼. 基于流域单元的地质灾害易发性评价——以西藏昌都市为例. 中国地质灾害与防治学报. 2023(01): 139-150 . 本站查看
    8. 翟越,谢怡帆,汪铁楠,王铭,张韵生. 基于改进突变理论的岩体崩塌危险性评价. 甘肃科学学报. 2023(02): 63-70 . 百度学术
    9. 孙璐,柯勇,何云勇,杨昌凤,刘道川. 川西高原山区公路典型灾害原因分析及防灾措施. 公路交通技术. 2023(02): 8-13+20 . 百度学术
    10. 潘世佳. 边坡工程地质灾害的成因分析及防治措施. 世界有色金属. 2023(05): 220-222 . 百度学术
    11. 刘文,余天彬,王猛,代力,黄细超,董继红. 四川宜宾市地质灾害隐患与地层岩性-地质构造关系分析. 中国地质灾害与防治学报. 2023(03): 118-126 . 本站查看
    12. 张燕,王庆兵,邢文超,王元新,张君,于桑,张建芝,葛江琨. 新疆某山区公路边坡危岩体影响区划分及防治建议. 安全与环境工程. 2023(04): 149-158 . 百度学术
    13. 王勇,管年汇,刘佳妍,王婉璐,吉雅娟. 基于GIS的上杭县地质灾害发育特征与主控因素分析. 河南城建学院学报. 2023(03): 92-98 . 百度学术
    14. 张宏蓉,杜娟,郑红卫,刘悦. 基于时空概率的公路崩塌灾害危险性评价. 安全与环境工程. 2023(05): 142-151 . 百度学术
    15. 李立青,何世阳,孙强. 基于专家系统的公路承灾体地质灾害分级. 公路与汽运. 2023(06): 95-100 . 百度学术
    16. 黄艳琴,李为乐,许洲,李鹏飞,铁永波. 国道213汶川—松潘段滑坡隐患遥感识别与易发性评价. 中国地质调查. 2022(04): 121-133 . 百度学术
    17. 梁荣伟. 泡沫轻质土在高速公路改拓建工程中的应用研究. 公路工程. 2022(04): 124-130 . 百度学术
    18. 白光顺,杨雪梅,朱杰勇,张世涛,祝传兵,康晓波,孙滨,周琰嵩. 基于证据权法的昆明五华区地质灾害易发性评价. 中国地质灾害与防治学报. 2022(05): 128-138 . 本站查看
    19. 周胜森,李为乐,陈俊伊,蒋瑜阳,王毅. 西昌—香格里拉拟建高速公路廊道不良地质体遥感识别. 中国地质灾害与防治学报. 2022(06): 90-102 . 本站查看
    20. 支泽民,刘峰贵,陈琼,周强. 青海省生态脆弱度分级1-km栅格数据集(2015)的研发与内容. 全球变化数据学报(中英文). 2022(03): 409-415+586-592 . 百度学术
    21. 南鹏军,刘国斌. 甘肃南部大水金矿集区地质灾害发育特征研究. 世界有色金属. 2021(24): 97-99 . 百度学术

    其他类型引用(7)

图(4)  /  表(18)
计量
  • 文章访问数:  276
  • HTML全文浏览量:  144
  • PDF下载量:  222
  • 被引次数: 28
出版历程
  • 收稿日期:  2020-09-15
  • 修回日期:  2020-09-21
  • 网络出版日期:  2021-10-19
  • 刊出日期:  2021-10-19

目录

/

返回文章
返回