Risk analysis of the geological hazards during urban tunnel construction in mountainous karst areas
-
摘要: 岩溶山区城市地下隧道工程施工过程中存在大量不确定地质灾害风险因素,由于成本和工期等原因无法详细查明所有水文地质和工程地质灾害风险源,导致各种灾害事故频发。目前对于隧道工程灾害评价方法多依靠专家打分,受人为主观因素影响大,且部分评价指标难以量化。为解决上述问题,依托岩溶山区轨道交通工程—贵阳市轨道交通2号线一期工程,统计其26段区间隧道施工过程中灾害发生位置、类型,仅选取明显相关地质灾害风险因素作为评价因素,根据灾害类型和地质灾害风险因素的位置关系反演出灾害发生频率与地质灾害风险因素间的耦合关系,建立岩溶山区城市隧道工程风险灾害评价体系。评价结果表明灾害发生频率与评价体系分段结果吻合程度均在69%以上,评价体系能够较好预测灾害发生频率。Abstract: Construction of urban tunnel engineering in mountainous karst regions involves a plethora of uncertain geological risk factors. Due to cost and schedule constraints, hydrogeological and engineering geological risk sources cannot be identified in detail, leading to frequent disasters. Expert scoring, which is influenced by human subjectivity, is a common method for disaster evaluation in tunnel engineering, and some evaluation indicators are difficult to quantify. In order to address these issues, Guiyang rail transit line 2 Phase I project was selected as a case study to establish a risk disaster evaluation system for urban tunnel engineering in mountainous karst regions. The location and type of the disasters during construction were recorded for 26 running tunnels, and geological risk factors that were significantly correlated were selected as evaluation factors. The coupling relationship between the frequency of disaster occurrence and the geological risk factors was inverted, leading to the establishment of the risk disaster evaluation system of urban tunnel engineering in mountainous karst region. The evaluation results demonstrate that the degree of consistency between the frequency of disaster occurrence and the segmentation results of the evaluation system is more than 69%, indicating that the evaluation system is capable of predicting the frequency of disaster occurrence effectively.
-
Keywords:
- mountainous karst region /
- tunnel engineering /
- risk evaluation /
- disaster frequency /
- back analysis
-
0. 引言
滑坡是我国发育最广泛的地质灾害之一,尤其是21世纪以来,随着人类工程活动的开展,滑坡发生频率上升,造成了巨大的人员伤亡和经济损失。对滑坡的风险评估和预测是目前防灾减灾的重要途径之一。随着计算机技术的发展,目前对滑坡风险评估逐渐由定性向定量转变,在定量评价中,堆积层厚度对滑坡识别具有重要意义,并且是一个重要的控制因素,能否对其进行准确预测是定量评价的关键。
目前,国内外对于土壤厚度的预测研究较为成熟,主要有基于物理机制、基于点采样插值、基于遥感的方法和基于环境相关推理的方法[1]。李岩等[2]采用多种插值方法对草原栗钙土层厚度空间分布进行预测,并选出适合研究区的最优插值方法。王桂林等[3]基于粒子群优化对协同克里金法进行了优化,结果表明该法对提高土层厚度的预测精度有较大作用。Scarphone等[4]采用广义线性模型(GLM)、随机森林(RF)、GLM残差克里格法和RF残差克里格法预测土壤厚度,结果表明GLM残差克里格法是预测浅层土壤厚度的最佳方法。Pradhan等[5]基于相对起伏度开发了一种空间土壤厚度模型,结果表明该模型的预测能力要优于多元线性回归模型和S模型。Li等[6]采用多元线性回归(MLR)、支持向量回归(SVR)、随机森林(RF)和极端梯度推进(XGBoost)4种常用机器学习方法分别对河南省土壤厚度进行预测。
堆积层厚度预测与土壤厚度预测有相似之处,区别在于堆积层厚度预测需要更多地考虑地质成因,国内对堆积层厚度的研究较少,多以空间插值为主,但对外界因素考虑不足,难以满足精度要求。刘磊等[7]根据已知区域堆积层厚度与其标志因子之间的关系,建立了基于随机森林算法的堆积层厚度模型,估算了整个万州主城区的堆积层厚度的空间分布。段贵娟等[8]以水平相关距离为基础数据,采用普通克里金法对南海区域的厚度分布进行了空间预测,预测效果良好。邵晨灿等[9]通过分析各成因类型土层厚度的影响因素,采用决策树模型建立影响因素与土层厚度之间的关系,推测了贵州省瓮安县的土层厚度分布。虽然目前国内外对于堆积层厚度的研究较少,但是对于土壤厚度的预测研究较为成熟,对堆积层厚度的研究可以借鉴土壤学中基于环境相关推理的方法,该方法假设相同的土壤景观单元发育相同的土壤类型,如Wang等[10]通过采用模糊C-均值法和模糊隶属函数构建土壤景观模型,从而预测中国三峡秭归地区土壤深度的空间变化。Akumu等[11]开发了一种基于GIS模糊逻辑建模方法的土壤深度建模和制图技术,制图精度达到94%,优于随机森林模型。Cheng等[12]提出了一种从单个土壤多边形中挖掘土壤-环境关系的方法,对传统土壤图进行了更新。基于环境相关推理法应用到堆积层厚度中为:相同地质环境下的堆积层厚度相似,可通过已知点的堆积层厚度与影响因子之间的关系预测整个区域的堆积层厚度。
三峡库区万州区铁峰乡地质灾害以滑坡为主,探索区内堆积层厚度对于滑坡稳定性评价显得尤为重要。鉴于此,本文以铁峰乡为研究对象,在实地调查的基础上,首先总结了区域堆积层的成因模式,并据此添加区域堆积层厚度控制点;其次定量分析了堆积层厚度与影响因子之间的相关性,并基于Apriori算法挖掘并分析影响因子与堆积层厚度的关联准则;最后,选用三种机器学习模型对研究区堆积层厚度进行空间预测,并对预测结果进行对比分析。研究结果揭示了区域第四系堆积层的成因机制,为机器学习在堆积层厚度预测领域奠定了基础。
1. 数据挖掘与预测模型
1.1 数据挖掘
数据挖掘流程主要由数据采集,数据预处理,数据离散化和关联准则挖掘四步组成。本文采用两步聚类算法及Apriori算法,通过数据挖掘流程处理堆积层厚度数据,得到堆积层厚度与影响因子间的关联准则。
(1)数据采集:将堆积层厚度进行收集归纳,获得基础数据。
(2)数据预处理:对已有数据进行审核、筛选、排序等处理措施。
(3)数据离散化:由于Apriori算法只能处理离散型变量,故采用两步聚类法将堆积层厚度离散化。
(4)关联准则挖掘:应用Apriori算法对离散后的变量进行关联规则集合选择,生成有效的关联规则。
通过两步聚类算法对堆积层厚度类型进行转化,将分类所得的变量采用Apriori算法计算,挖掘堆积层厚度的关联准则。
Apriori算法由Agrawal和Srikant提出[13],该算法首先生成高于最小支持度的频繁项目集,在第一步产生的频繁项目集中生成高于最小可信度的关联准则。
频繁项目集是指对包含项目a的项集T,其支持度大于或等于用户指定的支持度阈值,即:
(1) 包含k个项目的频繁项目集称为频繁k项集,记为
。如图1上层的项目集ab,abc,abcd,当满足最小支持度时均为频繁k项集。Apriori算法实现过程详见图1。首先,搜索出产生长度为1的频繁项集
, 又用于产生长度为2的频繁项目集 ,如此循环,搜索所有的频繁项目集。从频繁项目集中产生简单的关联准则,按置信度大于置信度阈值的条件,选择出有效规则集合。对每个频繁项目集L,计算L所有非空子集
的置信度,如果 大于用户指定的置信度阈值,即(2) 则生成关联准则
。1.2 GWO-SVM
支持向量回归机模型(Support vector regression,SVR)由美国学者Vapnik[14]提出,主要用于解决非线性问题的回归,其基本思想是:将样本数据分为两部分,分别为训练样本和测试样本,通过事先选择的非线性映射将输入向量(训练样本)映射到高维度特征空间,利用其有限的样本数据进行训练计算,求得达到最佳拟合效果时在这个空间中构造出的最优决策函数模型,对比测试项以分析模型的测试效果。该模型以统计训练理论为基础,实现小样本条件下的统计规律和学习,特点为结构简单,且性能比传统的BP神经网络有明显的优势。因此,本模型在堆积层厚度预测方面具有一定的优越性。
在SVM模型的参数优化方法中,灰狼算法(Grey wolf optimization,GWO)较粒子群算法(Particles swarm optimization,PSO)、遗传算法(Genetic algorithm,GA)等传统启发式算法而言,在研究组合式优化问题上已被证明具有显著优势,能够有效地解决多种参数的联合寻优问题。此外,GWO-SVM模型更适合小样本数据,且具有更好的学习效果。
GWO优化SVM的基本思想是将SVM的惩罚因子
、核函数参数 和不敏感参数 作为猎物,通过不断更新狼群位置,对目标猎物的最优解进行定位。利用得到的最优解构建SVM预测模型。本文设置的参数范围为: ; ; 。采用GWO-SVM耦合模型,对堆积层厚度进行预测,具体预测流程如下:
(1)输入训练集和测试集。
(2)初始化狼群数量、迭代次数,设置惩罚因子
、核函数参数 和不敏感参数 的范围。(3)GWO以
、 和 为猎物进行优化,达到最大迭代次数时输出最优值。(4)利用寻优得到的参数构建预测模型,利用该模型进行堆积层厚度预测。
2. 区域地质背景
2.1 研究区概况
铁峰乡位于重庆市万州区北部山区(图2),长江上游,东邻九池乡,南靠石龙乡、石和乡,西连两汇乡,北接开州区,总面积为51.46 km2。属亚热带季风气候,四季分明,雨量充沛。地势总体上南北两侧高,中间海拔低,以构造侵蚀中低山-河谷地貌为主。区内出露地层有侏罗系中统沙溪庙组和新沟组、侏罗系下统自流井组和珍珠冲组、三叠系须家河组和巴东组。第四系覆盖物广,主要为残坡积物、崩滑堆积物和河流冲洪积物。研究区位于铁峰山背斜北西翼核部,背斜全长约114 km,走向50°~70°,轴线呈凸向NNW的弧形,轴部地层为三叠系巴东组,两翼为侏罗系中、下统。
2.2 研究区地质灾害
铁峰乡发育的地质灾害类型包括堆积层滑坡、顺层岩质滑坡和岩质崩塌。根据地质灾害调查资料,研究区共发育历史滑坡灾害14处(堆积层滑坡12处,岩质滑坡2处),新生滑坡灾害12处(堆积层滑坡9处,岩质滑坡3处),堆积层滑坡主要分布于研究区南侧的顺向坡区段,目前多处于蠕滑变形状态;顺层岩质滑坡主要发育于自流井组和珍珠冲组地层内,以软弱夹层为主要的控滑结构,其稳定状态受公路开挖的影响极大。研究区崩塌灾害点主要分布于研究区斜交坡、北侧逆向坡结构的自然高陡边坡,以及公路开挖形成的高切坡之上,总计87处,以楔形体滑移式崩塌为主,占崩塌灾害点总量的56.3%,其次为坠落式崩塌,占比32.2%。
2.3 研究区堆积物成因模式
本次野外调查共确定了192个堆积层厚度数据控制点。对于全区而言,这些控制点不足以完成铁峰乡的堆积层厚度预测。因此,需在野外调查的基础上分析堆积层的成因模式,并据此添加一定量的堆积层控制点。根据铁峰乡坡体结构和岩层组合关系,对堆积层的成因进行分析总结,将研究区第四系堆积层分为残坡积和崩滑堆积两大类,崩滑堆积又可分为顺向坡远程滑动堆积、顺向坡近距离滑动堆积、顺向坡侧向崩塌堆积和软硬相间岩层崩滑堆积四种成因模式。
2.3.1 残坡积
指地表岩石风化后残留在原地的堆积物(图3),可区分为机械风化作用和化学风化作用。
(1)坡顶及原地残积物
由于重力和流水的搬运,宏观上沟谷及山麓地区多为冲、洪积物分布。残积物范围较为局限,山顶及未受扰动的风化原地堆积较为明显。常见的差异为较软的岩性风化堆积厚度较大,较硬的岩性厚度较薄。
(2)斜坡上部残积物
主要为风化堆积的残积物,未受重力作用及流水侵蚀而发生崩滑运动。该处第四系堆积形态与原岩较为相似,层状结构与地层一致。多位于发生过崩滑区域的上部,且相对于中下部厚度较薄。
(3)软硬相间岩性残积物
由于岩性差异形成的风化差异,不同的岩性处风化厚度不同,堆积形态也有所差异。整体上表现为原生堆积物在软岩处风化强烈,地形较缓且堆积较厚;在硬质岩处风化较弱,地形较陡而堆积较薄。
2.3.2 崩滑堆积
指因崩塌滑坡运动造成的堆积物,原始堆积物形态发生改变。可分为顺向坡远程滑动堆积、顺向坡近距离滑动堆积、顺向坡侧向崩塌堆积和软硬相间岩层崩滑堆积四种成因模式。
(1)顺向坡远程滑动堆积
原岩风化形成的堆积物发生远程滑动从而改变坡面原始形态。从纵剖面上来看,斜坡后部地形较缓,且堆积体较薄,主要为残留的滑坡堆积物和较薄的残坡积物;中前部由于滑坡下滑堆积,厚度较大,且易形成局部平台,主要物质为滑坡堆积物。若滑坡运动距离较大,滑坡堆积物可能覆盖坡底冲沟内的冲洪积物,甚至冲入河道,改变河道形态(图4)。
(2)顺向坡近距离滑动堆积
原岩风化形成的堆积物在内外地质作用的影响下沿软弱层面发生失稳破坏,改变原始地貌形态。原始地貌平面上为正常侵蚀堆积的山坡,沟谷之间过渡平衡,而发生滑动后横剖面中部可见明显的非正常侵蚀,堆积物下移而导致中间地形下错明显,平面上可见明显的洼地。剖面上原地形为自然过渡的缓坡,滑动部分导致前缘鼓胀、堆积较厚;后缘发生滑动出现地形下陷。坡顶为基岩,上部为风化后未完全滑动的残积物,中下部为风化后滑动堆积的滑坡堆积物,厚度较大,前部可能存在被滑坡堆积所部分覆盖的冲洪积物(图5)。
(3)顺向坡侧向崩塌堆积
顺向坡两侧冲沟陡崖在侧向侵蚀的作用下发生崩塌,崩积物堆积于陡崖下部,从而改变原始地貌。近崖处地形较陡,堆积层厚度较薄;随陡崖向外地形逐渐变缓,堆积厚度增厚;由于崩积物运动能量有限,远处堆积厚度逐渐变薄直到停止运动。由于崩坡积物粒径差异大,其形成的堆积结构特点明显,粒径较大的崩积物运动距离较大,多堆积于中前部,形成整体上前粗后细的堆积特点(图6)。
(4)软硬相间岩层崩滑堆积
软硬相间地层且临空条件较好,存在明显的差异风化特征,导致地形上出现了多级陡缓相间的差异变化,且不同位置堆积厚度及范围不同。平台堆积物厚度与陡崖的剖面形态相关,可分为三种类型:①坡下坡上陡崖高度一致,侧向的侵蚀范围相当,堆积厚度一定,以崩积物为主;②坡下陡崖高于坡下时,坡下堆积范围较坡上大,且厚度较大,以崩积物为主;③当侧向侵蚀同时存在滑坡时,坡上以崩积物为主,坡下滑动区则下伏滑坡堆积物,上覆新生崩积物(图7)。
3. 堆积物关联准则挖掘与预测
3.1 堆积层厚度影响因子
在研究区野外调查获得的192个堆积层厚度数据点的基础上,结合上述堆积层成因模式,添加包含0点和非0点在内的堆积层厚度控制点,累计共得到1211个数据样本点。根据堆积层成因模式和堆积特点,本文考虑了坡度、坡向、流向、地形起伏、剖面曲率、斜坡结构等13个影响因子。采用ArcGIS软件对研究区的1∶10000地形图进行栅格划分,并对上述13个因子进行提取并分级,统计各调查点的堆积层厚度,分析堆积层厚度的分布规律。本文选取坡度、坡向、流向、高程、地形起伏、剖面曲率、斜坡结构六个因子进行分析(图8)。
(1)坡度
从坡度因子图可以看出,研究区坡度主要集中在10°~40°范围内,该坡度范围内堆积层平均厚度为2.5~3.0 m。从堆积层平均厚度可以看出,堆积层厚度与坡度成反比,在坡度大于50°时,堆积层厚度接近0,而在0~40°范围内,堆积层厚度大于1.0 m,这表明在坡度低的地方,土层更容易堆积保留,从而形成一定的堆积层厚度,见图8(a)。
(2)坡向
从坡向因子图中可以看出,研究区内整体呈现出越接近于南北向的地方土层厚度越大的规律,在坡向近东西的区域内,堆积层的厚度普遍较小,通常小于0.5 m,在坡向为南北向的区域内堆积层的平均厚度达到最大值,通常在1.5~2 m。此规律与研究区内的地形地貌相吻合,研究区内坡向基本都是沿南北向分布,土体的崩滑运动主要也是沿这个方向发展,从而使土体在此方向上的堆积厚度大于其他方向,见图8(b)。
(3)流向
从流向因子图可以看出,堆积层平均厚度在大体上呈现出越容易向北汇水的区域堆积层厚度越厚,这与坡向因子相似,当堆积层点向北汇水时,该点更可能位于顺向坡内,而顺向坡更易发生崩滑运动,从而造成向北汇水点的堆积层厚度更厚,见图8(c)。
(4)地形起伏
从地形起伏因子可以看出,堆积层厚度分布规律与坡度因子完全一致,与地形起伏程度呈反比,地形起伏越小的区域内,土体越容易堆积,堆积层平均厚度在2.5~3 m之内,而当地形起伏达到4 m时,堆积层平均厚度仅为0.5 m,地形起伏超过5 m时,堆积层难以形成,平均厚度接近0。该规律与研究区地形相符合,由于研究区位于山区,地形起伏较大的区域往往也是坡度较大的区域,受地形限制,土体难以堆积,见图8(d)。
(5)剖面曲率
从剖面曲率因子可以看出,堆积层调查点主要集中在0~10、60~70、>70三个分级内。堆积层厚度分布呈现出与剖面曲率呈反比的规律,当剖面曲率位于0~10之内时,堆积层平均厚度达2.5~3 m,而当剖面曲率达到70时,厚度仅有0.5 m。剖面曲率大代表地表高程变化率大,即地势急剧变化,因此土层难以堆积保留,见图8(e)。
(6)斜坡结构
从斜坡结构因子可以看出,研究区主要为斜交坡和顺向坡。从堆积层厚度上看,土体在顺向坡中更容易堆积,从而形成一定的堆积层,顺向坡中的堆积层厚度在2.5~3 m之间,而斜交坡中堆积层厚度仅为1 m左右。这与地形作用息息相关,在顺向坡中,山体上部的土体由于崩滑等一系列作用从而向下移动,在顺向坡中往下堆积,从而在下部容易形成一定厚度的堆积层,见图8(f)。
3.2 关联准则挖掘
采用GRA(Grey relation analysis)和MIC(Maximal information coefficient)两种算法来分析上述13个影响因子与堆积层厚度的相关性,计算结果如表1所示。根据计算结果,各因子与堆积层厚度的GRA值差异不大,但是对于MIC算法,因子1、2、3、6、10、12与堆积层厚度的MIC值均大于0.2,表示具有较强的相关性。因此,本文选取以上6个因子(坡度、坡向、流向、地形起伏、剖面曲率、斜坡结构)作为关联前项,堆积层厚度作为关联后项,根据已知样本点数据挖掘两者的关联准则。
表 1 各因子与堆积层厚度相关性指标Table 1. Correlation index between each factor and accumulation layer thickness因子 GRA MIC 因子1 坡度 0.6923 0.6977 因子2 坡向 0.6878 0.5285 因子3 流向 0.8673 0.4696 因子4 高程 0.7050 0.0752 因子5 汇水条件 0.6349 0.0417 因子6 地形起伏 0.6856 0.6933 因子7 地层 0.7546 0.0486 因子8 距水系距离 0.6366 0.0525 因子9 距陡崖距离 0.6063 0.0637 因子10 剖面曲率 0.6953 0.2141 因子11 平面曲率 0.6783 0.0784 因子12 斜坡结构 0.6364 0.3271 因子13 植被覆盖 0.6440 0.0943 由于数据挖掘中Apriori算法仅能处理离散型变量,因此,需要对影响因子进行聚类。采用两步聚类算法对坡度、坡向、流向、地形起伏、剖面曲率、斜坡结构6个因子进行聚类,如表2所示。堆积层厚度按照(0~2 m)、(2~4 m)、(4~6 m)、(6~8 m)、(>8 m)划分为HD1、HD2、HD3、HD4、HD5五类。
表 2 各因子属性指标Table 2. Attribute indicators of each factor因子 属性 类别 因子1 坡度 0~20° F11 20°~40° F12 >40° F13 因子2 坡向 337.5°~67.5° F21 67.5°~157.5° F22 157.5°~247.5° F23 247.5°~337.5° F24 因子3 流向 1&2 F31 4&8 F32 16&32 F33 64&128 F34 因子6 地形起伏 0~2 F61 2~4 F62 >4 F63 因子10 剖面曲率 0~20 F101 20~40 F102 40~60 F103 >60 F104 因子12 斜坡结构 顺向坡 F121 逆向坡 F122 斜交坡 F123 根据已知样本点数据,以上述六个影响因子作为关联前项,堆积层厚度作为关联后项,设定支持度阈值为0,置信度阈值为100%,以前项最大为3生成堆积层厚度与影响因子的关联准则。共生成关联准则193条(部分准则如表3所示)。从关联准则可以看出,当堆积层厚度大于2 m时,坡度小于40°、地形起伏小于4 m,而当堆积层厚度小于2 m时,坡度大于40°、地形起伏大于4 m,说明坡度和地形起伏是造成堆积层厚度空间分布差异的主要因素。
表 3 堆积层厚度关联准则Table 3. Association criterion of accumulation layer thickness规则 ID 规则 支持度/% 置信度/% 1 YZ3=F31 & YZ6=F62 & YZ2=F21 HD50.08 100.00 2 YZ3=F31 & YZ2=F21 & YZ1=F12 HD50.08 100.00 3 YZ2=F22 & YZ3=F34 & YZ1=F11 HD30.08 100.00 4 YZ2=F22 & YZ3=F34 & YZ6=F61 HD30.08 100.00 5 YZ3=F34 & YZ1=F12 & YZ2=F24 HD20.17 100.00 6 YZ2=F23 & YZ10=F102 & YZ1=F11 HD20.08 100.00 7 YZ10=F102 & YZ2=F22 & YZ1=F11 HD20.08 100.00 8 YZ3=F31 & YZ6=F63 HD118.58 100.00 9 YZ3=F31 & YZ1=F13 & YZ6=F63 HD117.18 100.00 10 YZ2=F22 & YZ1=F13 & YZ6=F63 HD117.01 100.00 3.3 堆积物厚度预测
研究区共划分了486097个栅格点,其中1211个堆积层调查点,484886个预测栅格点。从1211个调查点中,随机选取1011个点作为训练样本,剩余200个作为测试样本,484886个未知堆积层厚度的栅格点作为预测样本,分别采用PSO-SVM、GA-SVM、GWO-SVM 3种机器学习模型对研究区堆积层厚度分布开展训练及预测。采用PSO、GA、GWO三种算法对SVM模型进行三参数寻优,各种优化方法的搜索区间如下:惩罚因子
;核函数参数 ;不敏感参数 ,其中GWO算法的寻优结果为: 。各模型测试结果及精度如图9、表4所示。表 4 各模型预测精度对比Table 4. Comparison of prediction accuracy of each model预测模型 MSE RMSE Q C 平均厚度/m PSO-SVM 401.22 1.42 0.47 0.79 1.49 GA-SVM 198.57 1.00 0.23 0.76 1.76 GWO-SVM 152.24 0.87 0.19 0.83 1.83 采用上述3种模型对区域堆积层厚度进行预测,预测结果如图10所示。根据图9和表4可知,GWO-SVM模型的测试结果与实际最为吻合,预测结果MSE为152.24,RMSE为0.87,精度最佳。另外,采用GWO-SVM的全区堆积层厚度预测结果平均值为1.83 m,与现场实际调查的结果较为相符(根据调查点数据,堆积层平均厚度为1.43 m)。因此,选取GWO-SVM模型得到的区域堆积层厚度分布图为最终结果,如图10(c)所示。
由图10(c)可知,颜色较深区域为堆积层厚度较薄地区,对应于研究区坡度较陡以及陡崖带附近,与上述坡度越大堆积层厚度越小的规律相符;从近陡崖处到远离陡崖的位置,堆积层厚度呈现出薄-厚-薄的规律;总体沿南北方向呈条带状分布,与研究区地形地貌相契合,在坡脚处堆积较厚。
4. 结论
(1)万州区铁峰乡第四系堆积层可分为残坡积和崩滑堆积两类,崩滑堆积有四种成因模式,分别为:顺向坡远程滑动堆积、顺向坡近距离滑动堆积、顺向坡侧向崩塌堆积和软硬相间岩层崩滑堆积。
(2)堆积层厚度与六种因子表现出较强的相关性:坡向、坡度、流向、地形起伏、剖面曲率和斜坡结构。其中坡度和地形起伏是造成堆积层厚度空间分布差异的主要因素。
(3)采用三种机器学习模型对铁峰乡堆积层厚度进行预测,结果显示GWO-SVM模型的预测结果与实际最为吻合,表明该方法可较好地应用于区域第四系堆积层厚度预测,为滑坡灾害识别提供更为可靠的依据。
-
表 1 灾害数量统计
Table 1 Satistics of Disaster quantity
灾害类型 发生次数 占比/% 涌水涌泥(坍塌) 86 47.5 地面塌陷(脱空) 71 39.2 管道破损 22 12.2 建(构)筑物沉降变形 2 1.1 表 2 涌水涌泥(坍塌)灾害次数统计
Table 2 Statistics of gushing water and bursting mud(collapse) disaster occurrences
岩土分界面与拱顶距离/m 次数 地下水位线与
拱顶距离/m次数 高于拱顶0~4 18 高于拱顶0~4 12 高于拱顶4~8 13 高于拱顶4~8 18 高于拱顶8以上 39 高于拱顶8以上 42 低于拱顶0~4 6 低于拱顶0~4 14 低于拱顶4以上 10 低于拱顶4以上 0 表 3 涌水涌泥(坍塌)灾害评价因素评分标准和权重
Table 3 Scoring criteria and weighting table for gushing water and bursting mud(collapse) risk evaluation factors
评价因素分级 分值计算 权重 岩土分界面高于拱顶:0~4 m和4~8 m 10-平均距离 0.50 岩土分界面低于拱顶:0~4 m 5-平均距离 岩土分界面高于拱顶超过8 m或低于拱顶超过4 m 0 地下水位线高于拱顶:0~8 m 10-平均距离 0.35 地下水位线低于拱顶:0~4 m 5-平均距离 地下水位线高于拱顶超过8 m或低于拱顶超过4 m 0 岩溶发育等级强 8.3 0.15 岩溶发育等级中 5.1 岩溶发育等级弱 1.8 注:当岩土分界面高于拱顶超过8 m或低于拱顶超过4 m时,地下水位线与拱顶距离评价因素分值直接取0。 表 4 地面塌陷(脱空)灾害评价因素评分标准和权重
Table 4 Scoring criteria and weighting table for surface collapse(void) risk evaluation factors
评价因素分级 分值 权重 岩溶发育等级强 8.3 0.55 岩溶发育等级中 5.1 岩溶发育等级弱 1.8 涌水涌泥(坍塌)灾害 涌水涌泥(坍塌)评分结果 0.45 表 5 管道破损灾害风险分级标准
Table 5 Pipeline breakage disaster risk classification standards
风险因素分级 风险等级 地面塌陷(脱空)灾害风险等级为一级,
且一级分段长度大于100 m一级 其他情况均为三级 三级 表 6 七机路口站—云峰路站区间隧道分段
Table 6 The running tunnel section of Qiji-Yunfeng station
岩土分界面与拱顶距离/m 里程 里程代号 高于拱顶0~4 YDK10+273.3~YDK10+462.5 1 YDK10+632.4~YDK10+745.3 3 高于拱顶4~8 YDK10+462.5~YDK10+632.4 2 YDK10+745.3~YDK10+809.5 4 低于拱顶0~4 YDK10+809.5~YDK11+018.3 5 低于拱顶超过4 YDK11+018.3~YDK11+182.2 6 表 7 七机路口站—云峰路站区间隧道涌水涌泥(坍塌)灾害风险评价结果
Table 7 The evaluation of gushing water and bursting mud(collapse) on Qiji-Yunfeng station running tunnel
里程
代号岩土分界面与拱顶平均距离 地下水与拱顶平均距离 岩溶发育等级 评价结果 风险等级 平均距离/m 分值 权重 平均距离/m 分值 权重 等级 分值 权重 1 高于拱顶2.6 7.4 0.5 高于拱顶5.9 4.1 0.35 强 8.3 0.15 6.4 一级 2 高于拱顶5.4 4.6 高于拱顶6.3 3.7 强 8.3 4.8 二级 3 高于拱顶3.6 6.4 高于拱顶6.2 3.8 强 8.3 5.8 一级 4 高于拱顶5.4 4.6 高于拱顶6.0 4 强 8.3 4.9 二级 5 低于拱顶2.1 2.9 高于拱顶4.5 5.5 强 8.3 4.6 二级 6 低于拱顶超过4 0 − 0 强 8.3 1.2 三级 表 8 七机路口站—云峰路站区间隧道地面塌陷(脱空)灾害风险评价结果
Table 8 The evaluation of surface collapse(void) on Qiji -Yunfeng station running tunnel
里程代号 岩溶发育等级 涌水涌泥(坍塌) 评价
结果风险
等级等级 分值 权重 分值 权重 1 强 8.3 0.55 6.4 0.45 7.4 一级 2 强 8.3 4.8 6.7 一级 3 强 8.3 5.8 7.2 一级 4 强 8.3 4.9 6.8 一级 5 强 8.3 4.6 6.6 一级 6 强 8.3 1.2 5.1 二级 表 9 七机路口站—云峰路站区间隧道管道破损灾害风险评价结果
Table 9 The evaluation of pipeline damage on Qiji-Yunfeng station running tunnel
里程代号 地面塌陷(脱空)风险等级 长度/m 风险等级 1 一级 189.2 一级 2 一级 170.0 一级 3 一级 112.9 一级 4 一级 64.2 三级 5 一级 208.8 一级 6 二级 163.9 三级 表 10 水井坡站—长岭路站区间隧道分段
Table 10 The geological summary of running tunnel section Shuijing-Changling station
岩土界面与拱顶距离/m 里程 里程代号 高于拱顶0~4 ZDK20+167.8~ZDK20+311.8 1 高于拱顶4~8 ZDK20+683.4~ZDK20+898.3 3 高于拱顶8以上 ZDK20+311.8~ZDK20+683.4 2 表 11 水井坡站—长岭路站区间隧道涌水涌泥(坍塌)灾害风险评价结果
Table 11 The evaluation of gushing water and bursting mud(collapse) on Shuijing-Changling station running tunnel
里程代号 岩土分界面与拱顶平均距离 地下水与拱顶平均距离 岩溶发育等级 评价
结果风险
等级平均距离/m 分值 权重 平均距离/m 分值 权重 等级 分值 权重 1 高于拱顶1.8 8.2 0.5 高于拱顶0.9 9.1 0.35 强 8.3 0.15 8.5 一级 2 高于拱顶8.1 0 高于拱顶2.9 0 强 8.3 1.2 三级 3 高于拱顶7 3 高于拱顶6.6 3.4 强 8.3 3.9 二级 表 12 水井坡站—长岭路站区间隧道地面塌陷(脱空)灾害风险评价结果
Table 12 The evaluation of surface collapse(void) on Shuijing-Changling station running tunnel
里程
代号岩溶发育等级 涌水涌泥(坍塌) 评价
结果风险
等级等级 分值 权重 分值 权重 1 强 8.3 0.55 8.5 0.45 8.4 一级 2 强 8.3 1.2 5.1 二级 3 强 8.3 3.9 6.3 二级 表 13 水井坡站—长岭路站区间隧道管道破损灾害风险评价结果
Table 13 The evaluation of pipeline damage on Shuijing-Changling station running tunnel
里程代号 地面塌陷(脱空)风险等级 长度/m 风险等级 1 一级 144.0 一级 2 二级 371.6 三级 3 二级 214.9 三级 表 14 阳明祠站—省医站区间隧道分段
Table 14 The running tunnel section of Yangmingci-Provincial hospital station running tunnel
岩土界面与拱顶距离/m 里程 里程
代号高于拱顶0~4 ZDK34+897.8~ZDK34+971.2 3 低于拱顶0~4 ZDK34+793.7~ZDK34+897.8 2 高于拱顶8 ZDK34+636.7~ZDK34+793.7 1 ZDK34+971.2~ZDK35+191.4 4 表 15 阳明祠站—省医站区间隧道涌水涌泥(坍塌)灾害风险评价结果
Table 15 The evaluation of gushing water and bursting mud(collapse) on Yangmingci-Provincial hospital station running tunnel
里程
代号岩土分界面与拱顶平均距离 地下水与拱顶平均距离 岩溶发育等级 评价结果 风险
等级平均距离/m 分值 权重 平均距离(m) 分值 权重 等级 分值 权重 1 高于拱顶8 0 0.5 0 0 0.35 中 5.1 0.15 0.8 三级 2 低于拱顶1.9 3.1 低于拱顶2.4 2.6 中 5.1 3.2 二级 3 高于拱顶2.3 7.7 低于拱顶2.8 2.2 中 5.1 5.4 一级 4 高于拱顶8 0 0 0 中 5.1 0.8 三级 表 16 阳明祠站—省医站区间隧道地面塌陷(脱空)灾害风险评价结果
Table 16 The evaluation of surface collapse(void) on Yangmingci -Provincial hospital station
里程
代号岩溶发育等级 涌水涌泥(坍塌) 评价
结果风险
等级等级 分值 权重 分值 权重 1 中 5.1 0.55 0.8 0.45 3.2 三级 2 中 5.1 3.2 4.2 二级 3 中 5.1 5.4 5.2 二级 4 中 5.1 0.8 3.2 三级 表 17 阳明祠站—省医站区间隧道管道破损灾害风险评价结果
Table 17 The evaluation of pipeline damage on Yangmingci-Provincial station
里程代号 地面塌陷(脱空)风险等级 长度/m 风险等级 1 三级 157.0 三级 2 二级 104.1 三级 3 二级 73.4 三级 4 三级 220.2 三级 表 18 灾害风险等级与灾害发生频率
Table 18 Disaster risk level and frequency of disasters
灾害类型 风险等级 灾害发生频率/(m·处−1) 涌水涌泥(坍塌) 一级 0~51 二级 51~200 三级 >200 地面塌陷(脱空) 一级 0~83 二级 83~210 三级 >210 管道破损 一级 0~270 三级 >270 表 19 与实际灾害发生频率相符的区间隧道长度占比
Table 19 Percentage of running tunnel length corresponds to the actual disaster frequency
灾害类型 与实际发生灾害
频率相符的
区间长度/m区间隧道总长度/m 准确率占比/% 涌水涌泥(坍塌) 16978.8 23254.2 73.0 地面塌陷(脱空) 16038.8 23254.2 69.0 管道破损 21167.2 23254.2 91.0 -
[1] 黄亮,赵卫权,吕思思,等. 喀斯特洞穴空间分布与河网密度关系研究—以贵州省为例[J]. 中国岩溶,2020,39(1):127 − 136. [HUANG Liang,ZHAO Weiquan,LYU Sisi,et al. Study on the relationship between spatial distribution of karst caves and drainage density:An example of Guizhou Province[J]. Carsologica Sinica,2020,39(1):127 − 136. (in Chinese with English abstract) HUANG Liang, ZHAO Weiquan, LYU Sisi, et al. Study on the relationship between spatial distribution of Karst caves and drainage density: an example of Guizhou Province[J]. Carsologica Sinica, 2020, 39(1): 127-136. (in Chinese with English abstract)
[2] 李红卫. 溶洞对贵阳轨道交通工程隧道影响的数值分析[J]. 土木工程与管理学报,2017,34(4):24 − 29. [LI Hongwei. Numerical simulation of the influence of cave on Guiyang railway tunnel[J]. Journal of Civil Engineering and Management,2017,34(4):24 − 29. (in Chinese with English abstract) LI Hongwei. Numerical simulation of the influence of cave on Guiyang railway tunnel[J]. Journal of Civil Engineering and Management, 2017, 34(4): 24-29. (in Chinese with English abstract)
[3] 张文,武科,刘国强,等. 城市地铁隧道岩溶地质灾害预警与成因分析[J]. 水利与建筑工程学报,2015,13(2):42 − 46. [ZHANG Wen,WU Ke,LIU Guoqiang,et al. Geological disaster prediction and cause analysis of karst cave in Guiyang subway construction[J]. Journal of Water Resources and Architectural Engineering,2015,13(2):42 − 46. (in Chinese with English abstract) ZHANG Wen, WU Ke, LIU Guoqiang, et al. Geological disaster prediction and cause analysis of Karst cave in Guiyang subway construction[J]. Journal of Water Resources and Architectural Engineering, 2015, 13(2): 42-46. (in Chinese with English abstract)
[4] 莫伟平. 岩溶地区涌水对地铁隧道轨道结构的影响及防治[J]. 城市轨道交通研究,2021,24(7):115 − 119. [MO Weiping. Impact and prevention of environmental water gushing on metro tunnel track structure in Karst areas[J]. Urban Mass Transit,2021,24(7):115 − 119. (in Chinese with English abstract) MO Weiping. Impact and prevention of environmental water gushing on metro tunnel track structure in Karst areas[J]. Urban Mass Transit, 2021, 24(7): 115-119. (in Chinese with English abstract)
[5] LI Liping,LEI Ting,LI Shucai,et al. Risk assessment of water inrush in karst tunnels and software development[J]. Arabian Journal of Geosciences,2015,8(4):1843 − 1854. DOI: 10.1007/s12517-014-1365-3
[6] 毛正君,杨绍战,朱艳艳,等. 基于F-AHP法的隧道突涌水风险等级评价[J]. 铁道科学与工程学报,2017,14(6):1332 − 1339. [MAO Zhengjun,YANG Shaozhan,ZHU Yanyan,et al. Risk grade evaluation of tunnel water inrush based on fuzzy analytic hierarchy process[J]. Journal of Railway Science and Engineering,2017,14(6):1332 − 1339. (in Chinese with English abstract) MAO Zhengjun, YANG Shaozhan, ZHU Yanyan, et al. Risk grade evaluation of tunnel water inrush based on fuzzy analytic hierarchy process[J]. Journal of Railway Science and Engineering, 2017, 14(6): 1332-1339. (in Chinese with English abstract)
[7] 吴全立,王梦恕,殷明伦. 基于网络分析法盾构下穿隧道风险分析[J]. 地下空间与工程学报,2019,15(6):1881 − 1888. [WU Quanli,WANG Mengshu,YIN Minglun. Risk analysis on shield TBM crossing underneath existing tunnel based on analytic network process[J]. Chinese Journal of Underground Space and Engineering,2019,15(6):1881 − 1888. (in Chinese with English abstract) WU Quanli, WANG Mengshu, YIN Minglun. Risk analysis on shield TBM crossing underneath existing tunnel based on analytic network process[J]. Chinese Journal of Underground Space and Engineering, 2019, 15(6): 1881-1888. (in Chinese with English abstract)
[8] NEZARAT H,SERESHKI F,ATAEI M. Ranking of geological risks in mechanized tunneling by using Fuzzy Analytical Hierarchy Process (FAHP)[J]. Tunnelling and Underground Space Technology,2015,50:358 − 364. DOI: 10.1016/j.tust.2015.07.019
[9] 张文泉,刘毓. 岩溶隧道突涌水危险性的多级模糊综合评价[J]. 西安科技大学学报,2016,36(2):187 − 192. [ZHANG Wenquan,LIU Yu. Multilevel fuzzy comprehensive evaluation of water inrush in karst tunnels[J]. Journal of Xi’an University of Science and Technology,2016,36(2):187 − 192. (in Chinese with English abstract) ZHANG Wenquan, LIU Yu. Multilevel fuzzy comprehensive evaluation of water inrush in Karst tunnels[J]. Journal of Xi’an University of Science and Technology, 2016, 36(2): 187-192. (in Chinese with English abstract)
[10] 匡星,白明洲,王成亮,等. 基于模糊评价方法的隧道岩溶突水地质灾害综合预警方法[J]. 公路交通科技,2010,27(11):100 − 103. [KUANG Xing,BAI Mingzhou,WANG Chengliang,et al. Research of comprehensive warning of water inrush hazards in karst tunnel based on fuzzy evaluation method[J]. Journal of Highway and Transportation Research and Development,2010,27(11):100 − 103. (in Chinese with English abstract) KUANG Xing, BAI Mingzhou, WANG Chengliang, et al. Research of comprehensive warning of water inrush hazards in Karst tunnel based on fuzzy evaluation method[J]. Journal of Highway and Transportation Research and Development, 2010, 27(11): 100-103. (in Chinese with English abstract)
[11] 蒋国云. 深埋隧道岩溶突水安全风险评价模型研究[J]. 地下空间与工程学报,2012,8(2):274 − 279. [JIANG Guoyun. Study on the security risk assessment model of karst water burst in deep buried tunnel[J]. Chinese Journal of Underground Space and Engineering,2012,8(2):274 − 279. (in Chinese with English abstract) DOI: 10.3969/j.issn.1673-0836.2012.02.009 JIANG Guoyun. Study on the security risk assessment model of Karst water burst in deep buried tunnel[J]. Chinese Journal of Underground Space and Engineering, 2012, 8(2): 274-279. (in Chinese with English abstract) DOI: 10.3969/j.issn.1673-0836.2012.02.009
[12] 何怡帆,李天斌,曹海洋. 隧道施工期岩爆危险性评价的属性识别模型及工程应用[J]. 水文地质工程地质,2020,47(2):102 − 111. [HE Yifan,LI Tianbin,CAO Haiyang. Attribute recognition model of fatalness assessment of rockburst in tunnel construction and its application[J]. Hydrogeology & Engineering Geology,2020,47(2):102 − 111. (in Chinese with English abstract) HE Yifan, LI Tianbin, CAO Haiyang. Attribute recognition model of fatalness assessment of rockburst in tunnel construction and its application[J]. Hydrogeology & Engineering Geology, 2020, 47(2): 102-111. (in Chinese with English abstract)
[13] SOUSA R L. Risk analysis for tunnelling projects[D]. Cambridge: Massachusetts Institute of Technology, 2010.
[14] 王. 城市地铁隧道事故案例统计分析与风险评价方法研究[D]. 北京: 北京交通大学, 2018 WANG Gong. Study on accidents statistical analysis and risk assessment methods of urban subway tunnel[D]. Beijing: Beijing Jiaotong University, 2018. (in Chinese with English abstract)
[15] 王智伟,王利,黄观文,等. 基于BP神经网络的滑坡监测多源异构数据融合算法研究[J]. 地质力学学报,2020,26(4):575 − 582. [WANG Zhiwei,WANG Li,HUANG Guanwen,et al. Research on multi-source heterogeneous data fusion algorithm of landslide monitoring based on BP neural network[J]. Journal of Geomechanics,2020,26(4):575 − 582. (in Chinese with English abstract) DOI: 10.12090/j.issn.1006-6616.2020.26.04.050 WANG Zhiwei, WANG Li, HUANG Guanwen, et al. Research on multi-source heterogeneous data fusion algorithm of landslide monitoring based on BP neural network[J]. Journal of Geomechanics, 2020, 26(4): 575-582. (in Chinese with English abstract) DOI: 10.12090/j.issn.1006-6616.2020.26.04.050
[16] 董英,张茂省,李宁,等. 城市地下空间开发利用的地质安全评价内容与方法[J]. 水文地质工程地质,2020,47(5):161 − 168. [DONG Ying,ZHANG Maosheng,LI Ning,et al. Methods and contents of geological safety evaluation for urban underground space development and utilization[J]. Hydrogeology & Engineering Geology,2020,47(5):161 − 168. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.201904060 DONG Ying, ZHANG Maosheng, LI Ning, et al. Methods and contents of geological safety evaluation for urban underground space development and utilization[J]. Hydrogeology & Engineering Geology, 2020, 47(5): 161-168. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.201904060
[17] KAMPMANN J, ESKESEN S D, SUMMERS J W. Risk assessment helps select the contractor for the Copenhagen metro system[C]//Proceedings of the World Tunnel Congress 98 on Tunnels and Metropolises, 1998, (1): 123 − 128.
[18] 黄仁东,吴寒,张惕,等. 基于云模型的岩溶隧道涌水灾害危险性评价及其在青岩头隧道的应用[J]. 中国地质灾害与防治学报,2018,29(5):44 − 51. [HUANG Rendong,WU Han,ZHANG Ti,et al. Evaluation of water burst hazard in Karst tunnel based on cloud model and its application in Qingyantou Tunnel[J]. The Chinese Journal of Geological Hazard and Control,2018,29(5):44 − 51. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2018.05.08 HUANG Rendong, WU Han, ZHANG Ti, et al. Evaluation of water burst hazard in Karst tunnel based on cloud model and its application in Qingyantou Tunnel[J]. The Chinese Journal of Geological Hazard and Control, 2018, 29(5): 44-51. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2018.05.08
[19] 李堃. 青岛地铁13号线区间隧道施工风险管理研究[D]. 青岛: 青岛理工大学, 2018 LI Kun. Study on construction risk management of Qingdao metro line 13 tunnel[D]. Qingdao: Qingdao Tehcnology University, 2018. (in Chinese with English abstract)
[20] 张飞. 贵阳市地铁1号线区间隧道施工安全风险评估[D]. 贵阳: 贵州大学, 2016 ZHANG Fei. Guiyang metro system line 1 tunnel construction safety risk assessment[D]. Guiyang: Guizhou University, 2016. (in Chinese with English abstract)
[21] 中华人民共和国住房和城乡建设部. 城市轨道交通地下工程建设风险管理规范: GB 50652—2011[S]. 北京: 中国建筑工业出版社, 2012. Ministry of Housing and Urban-Rural Development of the People's Republic of China. Code for risk management of underground works in urban rail transit: GB 50652—2011[S]. Beijing: China Architecture & Building Press, 2012. (in Chinese)
[22] 冉光炯, 杨翘楚, 王健. 特长隧道隧址区三维地质模型构建方案[J]. 吉林大学学报(地球科学版), 2022, 52(3): 930-940. RAN Guangjiong, YANG Qiaochu, WANG Jian. 3D geological modeling scheme for super long tunnel site[J]. Journal of Jilin University (Earth Science Edition), 2022, 52(3): 930-940.(in Chinese with English abstract)
[23] 潘涛. 软土地区双线区间盾构隧道施工对周边地表以及建筑物沉降的影响[J]. 水文地质工程地质,2022,49(1):101 − 108. [PAN Tao. Influences of double-track shield tunnel construction on settlements of adjacent ground and buildings in a soft soil area[J]. Hydrogeology & Engineering Geology,2022,49(1):101 − 108. (in Chinese with English abstract) PAN Tao. Influences of double-track shield tunnel construction on settlements of adjacent ground and buildings in a soft soil area[J]. Hydrogeology & Engineering Geology, 2022, 49(1): 101 − 108. (in Chinese with English abstract)
[24] 姜封国, 白丽丽, 宋敏, 等. 哈尔滨城市地铁大断面隧道施工稳定性分析[J]. 吉林大学学报(工学版),2020,50(4):1419 − 1427. [JIANG Fengguo, BAI Lili, SONG Min, et al. Analysis on construction stability of large section tunnel of Harbin City Metro[J]. Journal of Jilin University (Engineering and Technology Edition),2020,50(4):1419 − 1427. (in Chinese with English abstract) JIANG Fengguo, BAI Lili, SONG Min, et al. Analysis on construction stability of large section tunnel of Harbin City Metro[J]. Journal of Jilin University (Engineering and Technology Edition), 2020, 50(4): 1419 − 1427. (in Chinese with English abstract)
-
期刊类型引用(2)
1. 徐开慧. 差异性边界约束边坡稳定性分析及治理措施研究. 江苏建筑. 2024(02): 112-115 . 百度学术
2. 陆烨,黄越川. 基于DIC分析的含石量对碎石土边坡稳定性影响. 中国地质灾害与防治学报. 2023(04): 49-57 . 本站查看
其他类型引用(0)