2. 中国地质环境监测院, 北京 100081;
3. 长安大学, 陕西 西安 710054
2. China Institute of Geo-Environment Monitoring, Beijing 100081, China;
3. Chang'an University, Xi'an 710054, Shaanxi, China
山区城镇地质环境容量有限,随着山区城镇化建设的加快与建设规模的持续扩张,“向山发展、向沟发展”成为山区城镇化发展的必然趋势,削山造地、开挖坡脚势必诱发大量地质灾害[1],不仅威胁城镇居民的生命财产安全,同时,严重阻碍城镇化建设与发展。针对吕梁山区城镇边坡安全问题,首先要解决的是边坡风险分级问题,继而基于边坡的风险级别,采取相应的风险防控措施。因此,山区城镇边坡风险分级问题成为吕梁山区城镇化过程中亟待开展的一项研究工作。边坡风险分级是基于城镇边坡单元的危险性与危害性多个评价指标的综合性风险评价。边坡分级系统的概念最早由Pierson等[2]提出,将地质条件、危岩体规模、灾害历史、降雨等因素作为评价因子,开展公路边坡落石风险评估。在此基础上,一些学者[3-5]开展了公路边坡岩崩识别与分级系统研究,基于公路、铁路特点,建立并改进了边坡分级系统;2005年加拿大温哥华召开的“国际滑坡风险管理会议”提出了在土地利用方面地质灾害风险管理的理论框架[6]。依据地质灾害风险理念,Mölk等[7]针对土地利用的边坡开展了岩崩风险分级研究。相对于公路、铁路等线状工程开展的边坡分级研究,基于城镇的边坡风险分级研究则较少。香港曾分别于20世纪70年代末、90年代中期开展了定性的土质边坡分级研究与定量的边坡风险分级研究[8]。在以上边坡风险分级研究的基础上,随着山区城镇化建设与发展的需要,以城镇斜坡或边坡为单元的边坡风险分级研究成果不断涌现。唐亚明等[9]开展了城镇边坡风险分级研究,建立了基于层次分析法的半定量黄土边坡风险分级系统;易靖松等[10]基于城镇斜坡单元,开展了城镇边坡地质灾害风险分级评价。本文以吕梁山区吉县城区为研究区,开展吕梁山区城镇边坡地质灾害风险分级研究。结合研究区地质灾害风险调查与孕灾条件分析,优化吉县城区边坡风险分级系统评分指标;基于数理统计学方法,依据统计趋势线的自然拐点法,对边坡的危险性与危害性分级标准进行改进与厘定,建立城镇边坡风险分级表,并确定吉县城区边坡的风险。在此基础上,探索吕梁山区城镇边坡地质灾害风险分级方法与技术,为吕梁山区城镇边坡地质灾害风险分级与防控、土地利用规划及城镇化建设发展提供技术支撑。
1 研究区概况吉县位于山西吕梁山区的南部黄土地区(图 1),隶属于山西省临汾市。吉县县城属于黄土梁塬型河谷地貌,受地形条件的制约,城区地质环境容量有限,随着城镇化发展的需要,城区范围不断向清水河两侧的黄土支沟拓展,切坡、开挖等工程活动引发了大量的地质灾害及隐患,尤其是黄土崩塌灾害及隐患尤为突出(图 2)。吉县城区共发育地质灾害37处,其中,黄土崩塌33处,黄土滑坡4处,直接威胁人民群众的生命财产安全。
![]() |
图 1 研究区地理位置 Fig.1 Geographic location of the study area |
![]() |
图 2 吉县城区公路边坡滑塌(a)和土窑边坡崩塌(b) Fig.2 Collapses of highway slope(a) and cave-dwelling slope(b) in urban area of Jixian |
吉县城区属暖温带大陆性气候,四季分明,年平均气温10℃,极端最低气温零下21.3℃,每年一般于11月开始封冻,解冻始于次年3月上旬,最大冻土深度0.82 m,昼夜温差变化与冻融循环作用导致斜坡土体强度不断降低,是区内春季解冻期黄土崩塌、滑塌灾害频发的主要诱发因素。城区年降水量为560~580 mm,降雨时间集中,主要集中于每年的6—9月,秋季以连阴雨为主,集中降雨与长历时降雨的特征是区内汛期崩塌、滑坡灾害高发的主要因素。
吉县城区位于清水河河谷区,整体地形平缓,但城区南北两侧与黄土梁塬、黄土斜坡相接,受人类开挖、切坡修路、建房等工程活动影响,形成大量高陡边坡崩塌隐患。吉县主城区位于河谷阶地之上,全新统冲洪积层在研究区广泛分布,城区南北两侧为黄土斜坡,表部为上更新统黄土披覆,斜坡主体为中更新统黄土,沿清水河河谷及公路两侧零星出露中生界三叠系中统延长组(T2y)与二马营组(T2e)及下统和尚沟组(T1h)与刘家沟组(T1l),新生界新近系上新统(N1b)等(图 3)。
![]() |
图 3 黄土斜坡(a)和黄土-基岩斜坡(b)剖面 Fig.3 Profile of the loess slope(a) and loess-bedrock slope(b) |
城区边坡风险分级系统包括边坡危险性分级与危害性分级,评分指标亦分为危险性评分指标与危害性评分指标。评分指标是根据研究区边坡地质灾害及隐患的特征、类型与孕灾条件、承灾体类型等因素确定的评价指标因子,不同研究区应根据地质灾害类型、发育特征、孕灾条件、承灾体等,制定的符合实际城镇边坡的风险评价指标因子。
2.1 评分指标的选取在前人研究[9-10]的基础上,考虑山区城镇尺度地貌单元、地质构造、降雨条件等具有共性地质背景的条件下,结合吉县城区地质灾害的类型、特征,以及主控诱发因素、承灾体类型、特征等,认为吉县城区边坡危险性评分指标包括边坡的几何特征、灾害潜在破坏类型与强度、灾害发育的主控诱发因素等,具体包括边坡的地形条件、边坡结构类型、岩土体特性、地表水与地下水、人类工程活动、历史灾害、变形迹象、潜在地质灾害的类型及强度等9项、21个指标;危害性评分指标包括受城镇边坡地质灾害及隐患威胁的承灾体类型、数量、财产的价值、承灾体的易损性等,具体包括6项、15个指标。
2.2 评分指标权重的确定评分指标权重可以反映评价指标的相对重要程度及对评分结果的贡献程度。评分指标权重的计算方法有熵权法[11]、模糊数学法[12-13]、层次分析法[14-15]等。层次分析法因兼具定性经验判断与定量计算的优势[16],计算结果更加客观、准确,在地质灾害评价中应用广泛。前人[14, 17]通过建立层次分析模型、构造判断矩阵、合成权重计算、一致性检验等步骤得到危险性与危害性评分指标的权重。因前人对获取评分指标权重的方法、步骤有详细的阐述,本文不再对指标合成权重的获取过程与一致性检验进行详细阐述。利用该方法、步骤,可得到吉县城区边坡通过一致性检验的危险性指标合成权重,WHi=(0.03, 0.06, 0.02;0.01, 0.03;0.02, 0.07;0.03, 0.08;0.06;0.02, 0.04;0.06, 0.03, 0.02;0.06, 0.13;0.03, 0.05, 0.10, 0.05)(表 1)。同理,可得到吉县城区边坡危害性评分指标合成权重,WCi=(0.128, 0.256;0.074, 0.074;0.025, 0.012, 0.005, 0.002;0.039, 0.008, 0.023, 0.013;0.03;0.206, 0.103)。
![]() |
表 1 危险性指标合成权重 Table 1 Integrated weights of hazards factors |
在充分考虑地质灾害发育特征、孕灾背景条件、潜在灾害类型与强度等因素的基础上,结合前人评分标准的划分方法[9],依据评价指标对地质灾害危险性的贡献程度,吉县城区边坡危险性评分标准划分为5档,即1分、3分、5分、7分及10分(满分10分)5档分值(表 2)。单个指标的评分分值乘以其对应的权重,可得到该项指标的评分,各指标评分累加可得到边坡单元的危险性评分。危害性评分标准则是基于承灾体的类型与数量、承灾体的易损程度等评价指标对危害性的贡献程度,也划分为5个档次,即1分、3分、5分、7分及10分(满分10分)5档分值(表 3)。危害性的单项评价指标分值乘以其对应的权重,可得到该项指标的评分,并累加得到边坡单元的危害性评分。该风险评价方法相比常用的R=H*V评价方法,评价指标更全面、客观,基于研究区边坡的危险性与危害性分级标准建立了风险分级表,进而确定研究区边坡的风险,评价结果更贴近实际情况,也更具合理性。
![]() |
表 2 吉县城区边坡风险分级危险性评分 Table 2 Grading table for hazards of risk classification system of town slope in Jixian |
![]() |
表 3 吉县城区边坡风险分级危害性评分系统 Table 3 Grading table for consequence of risk classification system of town slope in Jixian |
边坡风险调查首先是对城镇边坡单元进行划分,依据划分的边坡单元,采用边坡风险分级表进行风险调查与评分。边坡单元的划分基于DEM数据,利用ArcGIS水文学方法实现边坡单元的自动划分,并在此基础上,依据城区边坡的自然形态,即坡顶以斜坡顶部的山脊线或坡折线为界,以冲沟或山梁作为斜坡单元的两侧天然分界线,结合斜坡的结构类型、坡向等因素对自动划分的单元进行人工修正,最终得到斜坡单元面数据集。吉县城区边坡单元总面积约5.9 km2,划分为135个边坡单元(图 4),平均边坡单元面积4.37×104m2。在此基础上,依据设计的吉县城区边坡风险分级调查表格逐坡逐段开展城区边坡风险调查,填表人依据边坡的危险性与危害性调查要素指标逐一进行调查、分析与判断,并给出各项指标的评分,得到边坡单元的危险性与危害性评分结果。
![]() |
图 4 吉县城区边坡单元划分 Fig.4 Unit partition of slopes in urban area of Jixian |
城区边坡风险分级是在危险性与危害性分级的基础上,建立城镇边坡风险分级表,并据该表判定边坡的风险级别。关于边坡危险性与危害性分级标准,有学者[9]按照边坡单元数的百分比作为划分标准,也有学者[10]依据边坡单元评分的自然间断法进行级别标准的划分,划分标准均具有较强的主观性,难以客观地反映边坡的危险性、危害性及风险级别。
在以上研究的基础上,本文对边坡的危险性与危害性分级标准进行了改进、优化。基于数理统计学理论[18],分别进行城区边坡单元危险性与危害性评分排序,依据统计趋势线的自然拐点对应的得分作为划分危险性与危害性级别的标准,并对拐点对应的边坡单元危险性与危害性进行核查,确保分级标准符合城镇边坡的实际情况,并将危险性与危害性分别划分为高、中、低、很低4个等级。根据危险性与危害性评分分布的趋势线自然拐点(图 5、图 6),结合对边坡的危险性与危害性实地核查,确定危险性评分在6.04及以上的地质灾害隐患为高危险,评分在4.49~6.04之间的为中危险,评分在3.38~4.49之间为低危险,评分低于3.38的为危险性很低;确定危害性评分为2.284及以上的为高危害性,在1.109~2.284之间的为中危害性,在0.954~1.109之间的为低危害,0.954以下的为危害性很低。依据边坡的危险性与危害性评分标准,判定吉县城区高危险边坡16处,中危险边坡78处,低危险边坡38处,很低危险边坡3处(图 7);吉县城区高危害边坡17处,中危害边坡82处,低危害边坡27处,很低危害边坡9处(图 8)。
![]() |
图 5 吉县城区及新区边坡危险性评分趋势图 Fig.5 Tendency chart of geohazards hazard scores of slopes in urban area and new urban area of Jixian |
![]() |
图 6 吉县城区及新区边坡危害性评分趋势图 Fig.6 Tendency chart of geohazards consequence scores of slopes in urban area and new urban area of Jixian |
![]() |
图 7 吉县城区边坡危险性分布图 Fig.7 Distribution of slope hazard in urban area of Jixian |
![]() |
图 8 吉县城区边坡危害性分布图 Fig.8 Distribution of slope consequence in urban area of Jixian |
依据边坡的危险性与危害性分级标准,建立吉县城区边坡风险分级表(表 4)。依据该表,判定吉县城区边坡的风险,即高风险边坡22处,中风险84处,低风险21处,风险很低8处(图 9)。同时,依据城区边坡风险分级结果,对城区边坡的风险进行了核验。经核验,采用该城镇边坡风险分级系统进行城区边坡风险评价的结果基本符合吉县城区边坡的实际情况。
![]() |
表 4 吉县城区边坡风险分级 Table 4 Risk classification for slopes in urban area of JixianCounty |
![]() |
图 9 吉县城区边坡风险分布图 Fig.9 Distribution of slope risk in urban area of Jixian |
(1) 基于数理统计学理论,依据统计趋势线的自然拐点法对城镇边坡的危险性与危害性分级标准进行了改进与厘定,并将危险性与危害性分别划分为高、中、低、很低4个等级。
(2) 依据边坡危险性与危害性分级标准,吉县城区高危险边坡16处,中危险边坡78处,低危险边坡38处,很低危险边坡3处;吉县城区高危害边坡17处,中危害边坡82处,低危害边坡27处,很低危害边坡9处。
(3) 根据边坡的危险性与危害性分级标准,建立了边坡风险分级表,判定吉县城区高风险边坡22处,中风险边坡84处,低风险21处,风险很低8处。高风险边坡建议采取搬迁避让,辅以工程治理措施;中风险边坡以工程治理措施为主,加强绿化与排水,有必要可采取避让搬迁;低风险边坡与风险很低边坡建议定期巡查,规范工程活动,加强地质灾害知识教育与宣传。
致谢: 西安科技大学叶万军教授对论文提出了宝贵的修改意见与建议,在此表示衷心感谢。
[1] |
张茂省, 薛强, 贾俊, 等. 山区城镇地质灾害调查与风险评价方法及实践[J]. 西北地质, 2019, 52(2): 125-135. |
[2] |
Pierson L A, Vickle R V. Rock-fall hazard rating system participants' manual[J]. U.S.Department of Transportation, Publication No.FHWA SA-93-057, 1993, 11: 1-99. |
[3] |
Stover B K. Highway rock-fall[J]. Colorado Geological Survey, Department of Natural Resources, Special Publication, 1992, 37: 27-27. |
[4] |
Hungr O, Fletcher L, Jakob M, et al.A system of rock fall and rock slide hazard rating for a railway-geo-hazard[M].Edmonton, Alberta, Canada, 2003: 277-283.
|
[5] |
Vandewater C J, Dunne W M. Classifying and assessing the geologic contribution to rock fall hazard[J]. Environmental & Engineering Geo-science, 2005, 11(2): 141-154. |
[6] |
Fell R, Corominas J, Bonnard C, et al. Guidelines for landslide susceptibility, hazard and risk zoning for land use planning[J]. Engineering Geology, 2008, 102(3/4): 85-98. |
[7] |
Mölk M, Poisel R, Weilbold J, et al. Rock fall rating systems:is there a comprehensive method for hazard zoning in populated areas?[J]. Interpraevent 2008-Conference Proceedings, 2008, 2: 207-218. |
[8] |
Wong H N.Landslide risk assessment for individual facilities[C]//Hunger O, Fell R, Couture R, et al.Proceedings of the International Conference on Landslide Risk Management, London: Taylor and Francis, 2005: 237-296.
|
[9] |
唐亚明, 薛强, 李清, 等. 黄土滑塌灾害风险分级系统研究[J]. 工程地质学报, 2012, 20(3): 378-386. DOI:10.3969/j.issn.1004-9665.2012.03.011 |
[10] |
易靖松, 张勇, 石胜伟, 等. 基于斜坡单元的山区城镇地质灾害风险快速评价研究——以江口镇为例[J]. 探矿工程(岩土钻掘工程), 2018, 45(8): 78-72. |
[11] |
王佳运, 何红前. 基于墒权法滑坡稳定性模糊综合评判[J]. 工程地质学报, 2007, 15(增刊): 129-132. |
[12] |
王秀英, 聂高众, 王松. 利用模糊数学方法建立汶川地震滑坡灾害评判标准[J]. 岩土力学, 2011, 32(2): 403-410. DOI:10.3969/j.issn.1000-7598.2011.02.014 |
[13] |
吴云, 张锐, 乐琪浪, 等. 基于模糊数学和Geo-Studio的滑坡稳定性综合评价[J]. 科学技术与工程, 2017, 17(23): 168-173. DOI:10.3969/j.issn.1671-1815.2017.23.027 |
[14] |
唐亚明, 程秀娟, 薛强, 等. 基于层次分析法的黄土滑塌风险评价指标权重分析[J]. 中国地质灾害与防治学报, 2012, 23(4): 40-46. |
[15] |
邱曼, 魏云杰. 基于层次分析法的叶城二牧场地质灾害危险性评价[J]. 防灾减灾学报, 2019, 35(3): 9-14. |
[16] |
张征. 环境评价学[M]. 北京: 高等教育出版社, 2004.
|
[17] |
许树柏. 实用决策方法——层次分析法原理[M]. 天津: 天津大学出版社, 1988.
|
[18] |
汪荣鑫. 数理统计[M]. 西安: 西安交通大学出版社, 1986.
|
[19] |
李瑞敏, 殷志强, 李小磊, 等. 资源环境承载协调理论与评价方法[J]. 地质通报, 2020, 39(1): 80-87. |
[20] |
曾毅, 项广鑫, 蒋星祥, 等. 基于国土空间自然适宜性的三类空间划分方法[J]. 地质通报, 2020, 39(1): 138-145. |