冰川冻土, 2023, 45(5): 1564-1574 DOI: 10.7522/j.issn.1000-0240.2023.0120

寒区工程与灾害

降雨对多年冻土区浅层滑坡失稳的影响研究

周保,1, 魏刚,2, 魏赛拉加1, 沈凌铠3, 张明哲3, 常文斌3, 臧佳园3

1.青海省地质环境监测总站，青海 西宁 810007

2.青海省环境地质勘查局，青海 西宁 810007

3.上海交通大学，上海 200240

Study on influence of rainfall on active layer detachment failure in the permafrost region

Bao ZHOU,1, Gang WEI,2, Sailajia WEI1, Lingkai SHEN3, Mingzhe ZHANG3, Wenbin CHANG3, Jiayuan ZANG3

1.Qinghai Geological Environment Monitoring Station，Xining 810007，China

2.Qinghai Environmental Geology Prospecting Bureau，Xining 810007，China

3.Shanghai Jiao Tong University，Shanghai 200240，China

通讯作者: 魏刚，正高级工程师，主要从事工程地质灾害调查、勘查及研究. E-mail：158213989@qq.com

收稿日期: 2022-01-25   修回日期: 2022-05-03  

基金资助: 青海省创新平台建设专项.  2021-ZJ-T08
青海省基础研究计划.  2019-ZJ-7053

Received: 2022-01-25   Revised: 2022-05-03  

作者简介 About authors

周保，正高级工程师，主要从事工程地质、地质灾害调查研究.E-mail：zhb820322@163.com

摘要

近年来青藏高原多年冻土地区降雨量呈增大趋势，导致活动层沿多年冻土层滑脱，诱发的冻土浅层滑坡灾害严重影响区域生态环境和人类活动。冻土浅层滑坡失稳是渗流、温度和应力复杂耦合的过程，明确降雨条件下多年冻土斜坡水热力响应机制，揭示降雨诱发冻土浅层滑坡失稳的机理十分关键。基于冻土水热力耦合数值模拟方法，建立了仅施加气温变化的模型一和在气温变化基础上施加强度为9 mm·d-1、持续降雨18 d的模型二，探讨了低强度、长时间降雨对多年冻土斜坡水热力演化的影响。结果表明：夏季雨水入渗对斜坡浅层温度场产生扰动，进而影响土体冻融过程，活动层以下有形成富水层的可能。雨水入渗导致融土饱和度大幅增加，水分渗流方向由竖直向下逐渐转变为顺坡方向。极限状态下斜坡位移分布在活动层，符合冻土浅层滑坡变形特征，降雨入渗数天后活动层位移有显著增大的趋势，最大位移所在位置向坡脚转移。降雨对斜坡稳定性影响显著，雨水入渗对活动层水热力产生持续影响，斜坡安全系数最小值出现明显滞后。研究结果为青藏地区冻土浅层滑坡灾害防治提供了科学指导。

关键词： 水热力耦合 ; 冻土融化 ; 顺坡渗流 ; 安全系数 ; 冻土浅层滑坡

Abstract

Under the influence of global climate change， rainfall in the permafrost region of the Qinghai-Tibet Plateau has tended to increase in recent years. The active layer of slope in the permafrost region is melting， a lot of water accumulates at the freezing-thawing interface， and the shear strength of soil decreases rapidly， which caused active layer detachments to occur frequently. Such natural hazards seriously affect the regional ecological environment， and human activities， and restrict social and economic development. Failure of active layer detachment is a complex coupling process of seepage， temperature， and stress. Driven by climate factors such as rainfall and temperature change， the hydro-thermal-mechanical evolution of the active layer is significant. Most scholars at home and abroad have researched the deformation of permafrost slopes in the process of slope failure and the variation of slope stability under rainfall conditions using model tests， numerical analysis， field monitoring， etc. Therefore， there is a lack of research on the mechanism of the hydro-thermal-mechanical response of permafrost slopes under rainfall conditions. This paper takes a typical active layer detachment disaster in Qilian County， Qinghai Province as the research object. Based on the numerical simulation method of hydro-thermal-mechanical coupling numerical of frozen soil， two models to discuss the effect of low intense and longtime rainfall on the hydro-thermal-mechanical evolution of permafrost slopes were constructed by COMSOL Multiphysics. Model 1 only considered the variation of temperature， model 2 applied continuous rainfall for 18 days with the intensity of 9 mm·d-1 under the condition of air temperature variation. By solving the partial differential equations of multiple physical fields， the numerical simulation of the hydro-thermal-mechanical evolution in the temperature change and rainfall process of permafrost slopes is realized， and the slope stability is calculated by the strength reduction method. The results are as follows： （1） the infiltration of rainwater in summer disturbed the temperature field of the active layer， Furthermore， the freezing and thawing process of soil is affected， even promoting the melting of fine soil with high ice content below the active layer to form a water-rich layer， and the pore water pressure in the basement of the active layer accumulates， which increases the risk of slope instability. （2） Rainwater infiltration causes a large increase in the saturation of the thawed soil， which increases the bulk density of soil and causes the mechanical property degradation of soil. The direction of the water seepage gradually changes from downward to parallel to the slope. It produces a seepage force parallel to the slope， which results in increased deformation of the active layer. （3） The displacement of the slope is distributed in the active layer in the limit state， which accords with the deformation characteristics of active layer detachments. The evolution of soil displacement is closely related to the seepage field formed by rainwater infiltration. The displacement of the active layer tends to increase significantly after several days of rainfall infiltration， and low intensity and longtime rainfall caused the location of maximum displacement moves towards slope toe. （4） Rainfall has a significant adverse effect on slope stability， the factor of safety of the slope decreases rapidly after rainwater infiltration for a while. Rainwater infiltration has a continuous influence on the thermal-hydro-mechanical evolution process of the active layer， the minimum factor of safety of the slope appeared to be several days lag behind the end of rainfall. The above results of the study provide some theoretical basis for revealing the mechanism of active layer detachments induced by rainfall and provide scientific guidance for the prevention of active layer detachments in the Qinghai-Tibet Plateau.

Keywords： coupling of hydro-thermal-mechanical ; thawing of frozen soil ; seepage ; factor of safety ; active layer detachment

PDF (5394KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

周保, 魏刚, 魏赛拉加, 沈凌铠, 张明哲, 常文斌, 臧佳园. 降雨对多年冻土区浅层滑坡失稳的影响研究. 冰川冻土[J], 2023, 45(5): 1564-1574 DOI:10.7522/j.issn.1000-0240.2023.0120

Bao ZHOU, Gang WEI, Sailajia WEI, Lingkai SHEN, Mingzhe ZHANG, Wenbin CHANG, Jiayuan ZANG. Study on influence of rainfall on active layer detachment failure in the permafrost region. Journal of Glaciology and Geocryology[J], 2023, 45(5): 1564-1574 DOI:10.7522/j.issn.1000-0240.2023.0120

0 引言

活动层沿多年冻土层滑脱引起的冻土浅层滑坡广泛分布于青藏高原多年冻土地区，是该区域热融灾害的主要类型之一［1-2］。近年来青藏地区气候变暖、降雨量增大趋势显著［3-4］，直接影响多年冻土稳定性，导致冻土浅层滑坡灾害频发［5-6］，造成高寒生态系统破坏，对社会经济发展产生威胁［7］。鉴于降雨是诱发冻土浅层滑坡的主要外部因素之一［8］，研究降雨对冻土浅层滑坡失稳的影响对青藏地区热融灾害的预防工作具有重要意义。

随着高寒地区重点工程的修建，降雨引发的多年冻土地区斜坡稳定性问题受到了广泛关注，已有学者通过模型试验、数值分析和现场监测等手段开展了降雨对冻土斜坡稳定性影响的研究。Vidie等［9］通过模型试验研究不同降雨强度对冻土斜坡失稳的影响，提出中等降雨会引起冻胀蠕变和小规模滑坡，强降雨会诱发规模大、速度快的滑坡和泥流；段东明等［10］采用数值模拟对多年冻土斜坡路基进行降雨入渗分析，得出降雨导致斜坡稳定性的减小存在一定滞后；Niu等［11］基于无限斜坡模型评估多年冻土斜坡稳定性，计算夏季降雨触发冻土浅层滑坡的临界地下水位；赵文等［12］开展模型试验研究冻土斜坡在反复降雨-日晒-冻融循环条件下的演化特征，认为冻融循环引起坡面土体密实度降低，松散土体在降雨过程中流失导致斜坡逐年后退。

冻土浅层滑坡失稳是一个水热力动态耦合过程，雨水入渗和气温变化均会影响冻土斜坡内部的水热迁移［13-14］，但现有的研究主要关注降雨条件下冻土斜坡失稳过程及稳定性变化规律，气温变化作用下降雨对多年冻土斜坡水热力影响机制尚不明确。本文通过对青海省祁连地区冻土浅层滑坡开展现场调查，总结了当地冻土浅层滑坡发育的地质环境特征与气候条件，以典型冻土浅层滑坡灾害为研究对象，基于现场勘察和原位钻孔资料，根据当地夏季降雨特征［15］，采用冻土水热力耦合数值方法，系统研究了气温变化条件下低强度连续降雨对多年冻土斜坡水热力演化的影响，研究结果为青藏地区冻土浅层滑坡灾害防治提供了理论依据。

1 青海省祁连县冻土浅层滑坡

1.1　地质环境特征

基于多源遥感数据调查青海省多年冻土区热融灾害分布，共解译冻土浅层滑坡灾害290处，在青海省祁连县、治多县和曲麻莱县广泛分布，结合现场调查验证，祁连县共发育有冻土浅层滑坡灾害54处，占青海省该类型灾害总数的18.6%。采用无人机航拍、动力触探和钻探等方式对祁连县冻土浅层滑坡开展了现场调查。

青海省祁连县位于祁连山中段腹地，地形呈北高南低、西高东低的特点，地貌类型较丰富。经统计，祁连县冻土浅层滑坡发育的海拔范围为3 550~3 850 m，主要集中在3 750 m；发育位置多为汇水沟中上部，斜坡坡度集中在10°~20°，整体上陡下缓，发育坡向多为阴坡；斜坡表层植被以异针茅和高山蒿草为主，上部植被较下部稀疏，植被总体覆盖率约50%~80%，滑体表层植被破坏严重；滑体表面由第四纪上更新统坡（洪）积物（Q3dl+pl）组成，下伏基岩主要为三叠系（T）砂岩，活动层土体主要为细粒黏土和泥炭，受冻融循环影响土体力学性质较差；后缘往往发育小型不规则陡坎，底部常有地下水汇集或地下冰出露；滑体向沟口流动堆积，前缘多发育有鱼鳞状草皮。

1.2　气候条件

青海省祁连县属于典型的高原大陆性气候，地形条件复杂和海拔梯度悬殊导致水热条件呈显著的区域性差异和垂直差异。东部降雨多、气温较高，西部降雨少、气温较低，低海拔地区的气温较高海拔地区高；年内气温变化符合正弦函数分布，1月平均气温达到最低值-11 ℃，7月平均气温达到最高值15 ℃；研究表明当地多年平均气温呈波动上升趋势，气温倾向率达0.33 ℃·（10a）-1［16］。

根据当地气象站资料可得：从祁连县降雨量月际分配来看，5—9月多年平均降雨量达到324 mm，占全年总量的84.6%，7月往往达到降雨量最大值；从祁连县降雨量的季节分配来看，夏季多年平均降雨量为234 mm，占全年降雨量的57.8%，秋季和春季平均降雨量分别占20.9%和17.8%，冬季仅占3.5%。研究表明近年来青海省祁连山地区降雨量呈波动上升趋势，尤其进入21世纪后，降雨量呈快速增大趋势，降雨量倾向率为46.3 mm·（10a）-1［17］。

1.3　滑坡概况

2018年8月末青海省祁连县默勒镇海浪村（100°49′58″ E， 37°43′47″ N）发生冻土浅层滑坡（图1）。现场调查表明，滑坡高程范围为3 763~3 800 m，平面呈不规则舌状，位于阴坡中下部，斜坡整体坡度约15°。根据运动及堆积特征将滑坡分为滑源区和流动堆积区，滑源区长32 m，流动堆积区长58 m。根据2020年10月10日钻孔资料获取滑坡地层信息（图2），地表以下0~1.6 m为活动层，由粉质黏土、砂质黏土和泥炭组成，其中0~1.4 m土层较松散，1.4~1.6 m土层相对致密；1.6~12.7 m为多年冻土层，由黏土和泥岩组成，含有大量冰晶；12.7~15.0 m为基岩层，主要为砂砾岩。结合现场调查和钻孔资料获取主滑方向地质剖面如图3所示。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   祁连县冻土浅层滑坡

Fig. 1   Scene photograph of the active layer detachment in Qilian


图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   钻孔地层信息

Fig. 2   Stratum information based on borehole


图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   滑坡体主滑剖面

Fig. 3   Section of the main sliding mass


图4为根据气象站观测资料所得祁连县降雨概况。根据祁连县2008—2020年的年降雨量统计图［图4（a）］，2018年降雨量达574 mm，为2008—2020年间降雨量峰值，2018年6—8月降雨量占全年总降雨量67.6%。根据2018年6—8月的日降雨量统计图［图4（b）］，当地6月22日至7月9日（冻土浅层滑坡发生前）连续18 d总降雨量为162 mm，日均降雨量为9 mm·d-1，最大日降雨量为33 mm·d-1。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   青海省祁连县降雨概况

Fig. 4   Rainfall situation of Qilian， Qinghai Province： annual rainfall in Qilian from 2008 to 2020 （a）； daily rainfall in Qilian from June to August in 2018 （b）


2 水热力耦合数值模拟

2.1　计算原理

基于有限元软件COMSOL Multiphysics求解多物理场偏微分方程，实现多年冻土斜坡气温变化和降雨过程中的水热力演化数值仿真，并用强度折减法计算斜坡稳定性［18］。鉴于气温变化和降雨对冻土浅层滑坡的影响是复杂的水热力动态耦合过程，为简化这一过程，假设入渗雨水温度与地表温度一致，温度场变化受热传递、水分运移和冰水相变影响；水分场变化由基质吸力驱动，并由孔隙冰阻隔作用控制；水热过程单向影响土体应力和位移。相应物理场控制方程如下：

（1）温度场方程

冻土二维传热中，忽略对流传热、考虑相变的温度场方程可写为［19］

ρC(θ)∂T∂t=λ(θ)∇2T+LρI∂θI∂t

式中：ρ为土体的密度（kg·m-3）；C(θ)为土体的热容量（J·kg-1·℃-1）；λ(θ)为土体导热系数（W·m-1·℃-1）；T为土体温度（℃）；L为水冻结成冰的相变潜热，取335 kJ·kg-1，ρI为冰的密度（kg·m-3）；θI为体积含冰量。

（2）水分场方程

冻土中水分迁移规律与非饱和土中水的运移相似，用Richards方程描述［19］

∂θu∂t+∇(-D(θu)∇θu-kg(θu))+ρIρw∂θI∂t=0

式中：D为未冻土的扩散率，未冻土的扩散率可写为D=kc，（m2·s-1），冻土中孔隙冰的存在削弱了水分的迁移效应，因此，引入阻抗系数考虑孔隙冰的隔水作用［20］，冻土扩散率写为D=kc⋅I，（m2·s-1），I=1010θI；θu为体积含水率；ρw为水的密度（kg·m-3）；k为导水率（m·s-1）；c为比水容量（1·m-1）；kg为重力方向导水率（m·s-1）。

（3）应力场方程

土体任意微元的平衡微分方程和土体位移-应变关系表达式可写为［21］

∂Tσ-F=0
ε=-∂u

式中：∂T为微分算子矩阵；σ为应力张量，σ=σx  σy  σxyT；F为单位土体体积力，F=Fx  FyT；ε为应变，ε=εx  εy  εzT；u为位移，u=ux  uyT。

2.2　网格模型和计算参数

根据地层信息（图1）和主滑剖面（图2）建立二维有限元网格计算模型如图5所示，斜坡模型长90 m，坡度为15°，地表以下0~1.6 m为活动层，1.6~12.7 m为多年冻土层，12.7 m以下为基岩层；采用自由三角形网格，由于活动层水热力演化剧烈，故将活动层网格细化。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   有限元计算模型

Fig. 5   Finite element model


冻土比热容和导热系数与未冻水含量和含冰量的关系［22-23］、抗剪强度参数与土温的关系［24］参考相关研究设置，土水特征曲线和渗透系数采用Van Genuchten提出的模型［25-26］。土骨架的比热容和导热系数分别取1.4×106 J·m-3·℃-1和1.3 W·m-1·℃-1；活动层深度0~1.4 m土层的饱和渗透系数取2×10-6 m·s-1，1.4~1.6 m土层饱和渗透系数取4×10-9 m·s-1［27］；其余数值模拟所需参数根据钻孔取样进行室内试验的结果设置如表1所示。

表1   地层物理力学参数

Table 1  Physical and mechanical parameters of the formation

参数活动层多年冻土基岩
重度/（kN·m-3）182025
弹性模量/MPa40305 000
泊松比0.250.30.15
饱和渗透系数/（m·s-1）2×10-6/4×10-98×10-100
黏聚力/kPa1235/
内摩擦角/（°）2220/
初始体积含水率0.220.430.1

新窗口打开| 下载CSV


2.3　计算方案及约束条件

建立不考虑降雨的模型一和考虑降雨的模型二，对照分析降雨对冻土浅层滑坡的影响：模型一仅施加气温变化工况模拟斜坡10月10日至次年12月31日的水热力演化；模型二以模型一7月1日的水热力计算结果为初始状态，在气温变化的基础上施加强度9 mm·d-1、历时18 d的降雨工况，模拟斜坡7月1日至7月18日降雨过程中和7月19日至7月30日降雨停止后12 d内的水热力演化。

根据当地气温变化情况和气候变暖特征［28］，结合附面层理论［29］将地表温度视为平均气温并简化成三角函数形式作为模型上表面温度边界条件，表达式如下：

T=2+13sin(2πt3.1104×107+π2)+0.037t3.1104×107

式中：t为时间（s），模型下表面温度边界条件为热流边界条件［30］，根据深层地热梯度与土体导热系数的乘积进行换算，取0.06 W·m-2；模型两侧为绝热边界；将钻孔实测地温设为初始值，下表面采取热流边界条件，计算时间设为100年，使地温分布基本稳定，结果作为模型一2020年10月10日温度场初始状态。模型两侧和下表面为不透水边界条件；不考虑降雨时上表面不透水边界条件，考虑降雨时上表面为流量边界；如表1所示，土体初始体积含水率作为水分场初始状态。模型上表面为自由位移边界条件，两侧水平位移为0，下表面为固定位移边界条件。

3 模拟结果与分析

3.1　降雨对温度场的影响

2018年7月1日至18日模型一和模型二地表以下0.4 m和1.6 m处地温变化分别如图6（a）和图6（b）所示：仅施加气温变化时，7月18日0.4 m深度地温升高0.90 ℃，1.6 m深度地温升高0.087 ℃；在气温变化基础上施加降雨工况后，7月18日0.4 m深度地温升高1.84 ℃，1.6m深度地温升高0.18 ℃。降雨导致0.4 m深度地温增长速率随时间减小，1.6 m深度地温增长速率随时间增大。活动层基底存在0.2 m厚的相对致密层，阻碍雨水入渗，降雨0~13 d对1.6 m深度地温不产生影响。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   降雨导致不同深度地温随时间变化

Fig. 6   Ground temperature at different depths varies with time caused by rainfall： ground temperature at a depth of 0.4 m varies with time （a）； ground temperature at a depth of 1.6 m varies with time （b）


3.2　降雨对水分场的影响

图7为7月18日和7月30日模型一和模型二的含冰量云图，结果表明：7月18日模型一融深位于1.35 m，模型二融深位于1.49 m，降雨18 d导致模型二融深下降14 cm，多年冻土上限位置含冰量较模型一减小2.94%；7月30日模型一融深位于1.43 m，模型二活动层已完全融化，融深位于1.60 m，可见雨水渗流对土体含冰量产生持续影响，降雨停止后12 d模型二融深下降17 cm，模型二多年冻土上限位置含冰量较模型一减小40.74%。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   7月18日和7月30日模型一和模型二含冰量

Fig. 7   Ice content of model 1 and model 2 on July 18th and July 30th


图8为7月18日和7月30日模型一和模型二的饱和度云图。7月18日模型一0~1.4 m为饱和度约0.58的融土层，活动层底部20 cm为饱和度约0.45的冻土层；7月30日模型一融深增大，冻土层厚度减小为10 cm，饱和度增至0.48。雨水入渗导致活动层饱和度大大增加，7月18日模型二融土层的饱和度约0.88，活动层底部13 cm为饱和度约0.45的冻土层；由于水分下渗和冻土进一步融化，7月30日活动层完全融化，饱和度约0.85，活动层基底以下形成10 cm厚饱和度约0.90的富水层。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   7月18日和7月30日模型一和模型二饱和度

Fig. 8   Saturation of model 1 and model 2 on July 18th and July 30th


根据不同深度渗流速度随时间变化曲线［图9（a）］和降雨过程中流速分布［图9（b）］，分析降雨条件下渗流速度变化规律如下：7月1日1.2 m深度内土体已完全融化，渗流速度随深度减小，此时融水沿坡面法向迁移，降雨0~3 d雨水入渗导致流速迅速增加，此时雨水渗透深度有限，流速增大速率随深度增大而减小；随着雨水继续下渗，自降雨3 d起，流速增大速率随深度增大而增大；1.4 m深度土体的孔隙冰在降雨9 d融化，流速开始增大；降雨9 d雨水沿重力方向竖直下渗，在孔隙冰的阻隔作用下，降雨18 d随土层深度增加，渗流方向愈发接近顺坡；降雨停止后，1.4 m深度范围内土体流速下降，此时流速随深度增大而增大；由于1.4~1.6 m深度为渗透性较差的相对致密层，1.6 m深度处流速在降雨停止后2 d才逐渐增大，且流速增量很小；降雨停止后12 d，渗流方向完全转变为顺坡。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9   降雨条件下渗流速度变化规律

Fig. 9   Variation of seepage velocity under rainfall condition： seepage velocity at different depths varies with time （a）； distribution of seepage velocity at different stages （b）


3.3　降雨对斜坡稳定性的影响

图10（a）为模型一年内的安全系数变化曲线，可以得出斜坡稳定性与气温均呈正弦函数变化，且两者呈负相关，稳定性相对气温具有一定滞后，最高气温和最低气温出现在7月和1月，而8月和9月斜坡稳定性最差，2月稳定性最好，年内气温升高导致安全系数下降约12.1%。

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10   模型一和模型二稳定性变化规律

Fig. 10   Variation of slope stability of model 1 and model 2： variation of safety factor of model 1 （a）； variation of safety factor of model 2 （b）


图10（b）为模型二7月1日至7月30日安全系数变化曲线，结果表明：降雨导致安全系数下降约49.1%。降雨0~6 d安全系数缓慢下降，6~15 d安全系数下降趋势显著增大，安全系数从4.30下降至2.61，降雨15 d安全系数下降趋势逐渐减小，于降雨停止后10 d达到最小值2.22，此时斜坡仍处于稳定状态。

图11（a）直观反映了不同时刻斜坡极限状态下的位移。降雨条件下斜坡位移主要集中于活动层，降雨0、9、18 d活动层最大位移分别为5.2 mm、8.4 mm和35.9 mm，降雨停止后第12 d达到19.1 cm，最大位移所在位置从坡顶逐渐向坡脚移动，降雨停止后12 d最大位移位于斜坡中部和坡脚间。根据斜坡极限状态下不同位置的位移随时间变化曲线［图11（b）］，降雨期间位移缓慢增大，坡顶的位移较大，坡脚位移较小；降雨入渗至第15 d，坡表变形迅速发展，其中坡中和坡脚位移增大趋势尤为显著，坡顶位移增幅相对较小。

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11   极限状态下模型二位移变化规律

Fig. 11   Variation of displacement of model 2 under limit state： distribution of slope displacement under limit state （a）； variation of slope displacement at different positions under limit state （b）


4 讨论

降雨18 d导致0.4 m深度地温升高0.94 ℃、1.6 m深度地温升高0.093 ℃，可见雨水入渗导致斜坡浅层地温升高，降雨对地温的影响随深度增加而减小。地温升高促进冻土融化，同一时刻降雨条件下活动层融深较大，降雨停止后12 d多年冻土上限附近含冰量减小40.74%。雨水入渗和冻土融化导致活动层饱和度提高77%~95%，降雨停止后12 d活动层以下出现了10 cm厚富水层。分析认为，降雨导致活动层融土含水量增加、容重增大，水分聚集对土颗粒产生动态浮托力［31］，土体力学性质下降；7月1日至18日地表温度达到最大值，地温随深度增加而降低，雨水下渗过程中内能向温度较低的土体转移，且活动层饱和度上升提高了土的导热性，有利于大气与土体之间的热传导；降雨条件下融深可能进一步增大，活动层以下细粒富冰冻土层融化产生富水层，孔隙水压难以消散，土体有液化的可能［32］，且水分的润滑作用削弱了冻融交界面的抗滑力［33］，发生冻土浅层滑坡的风险大大增加。

极限状态下斜坡的位移集中在活动层，潜在滑面位于活动层和多年冻土交界面，符合冻土浅层滑坡变形特征。降雨初期水分沿重力方向下渗，此时位移缓慢增大；当雨水入渗量较大时，地下冰起隔水作用，导致渗流方向逐渐转变为顺坡，流速随深度增大而增大，同一深度斜坡中部流速最大，此时位移显著增大，最大位移所在位置由坡顶附近转移至坡中和坡脚之间。由此可见，降雨条件下融土达到一定饱和度后，活动层形成顺坡方向的渗流场，产生顺坡方向的渗透力［34］，增大了上覆融土的下滑力，此时土体饱和度较高，力学性质降低，导致斜表变形加剧，增大了活动层沿多年冻土层滑脱的可能性。

斜坡稳定性与气温均以正弦函数形式变化，两者呈负相关，相关研究表明，气温对多年冻土斜坡稳定性的影响需要进一步考虑短时间尺度下的极端高温［35-36］和长时间尺度下的冻融循环作用［8］。施加降雨工况后，斜坡安全系数从4.30下降至2.22，下降幅度达48.4%，降雨入渗前3 d安全系数下降速率缓慢，3~15 d安全系数迅速下降，第15 d起安全系数下降速率逐渐降低，于降雨停止后10 d达到最小值，此时斜坡仍处于稳定状态，推测后续的降雨将进一步降低斜坡安全系数导致斜坡在8月末失稳。综上所述，降雨对多年冻土斜坡浅层稳定的影响显著，雨水入渗对多年冻土斜坡浅层水热力产生持续性的影响，在低强度、长时间的降雨工况下，斜坡最小安全系数出现时间滞后于降雨数天，斜坡可能在降雨停止后的数天失稳，推测此时斜坡的破坏形式为坡脚附近活动层最先沿冻融交界面薄弱层破坏，从而牵引上部土体下滑引起冻土浅层滑坡。

5 结论

通过COMSOL Multiphysics建立了多年冻土斜坡模型，在气温变化的基础上模拟了9 mm·d-1、持续18 d降雨条件下多年冻土斜坡水热力演化过程，基于模拟结果探讨了气温变化条件下青藏地区低强度、长持时降雨对冻土浅层滑坡失稳的影响，得出以下结论：

（1）夏季降雨扰动斜坡浅层温度场，影响了活动层的冻融过程，可能导致活动层以下含冰量较高的细粒土融化形成富水层，增大了斜坡因活动层基底孔隙水压力积累而失稳的风险。

（2）雨水入渗导致活动层饱和度大幅提高，增大了土体容重、削弱了土体力学性质，在多年冻土层隔水作用下水分沿顺坡方向迁移，产生顺坡方向的渗透力，加剧活动层变形，对多年冻土斜坡浅层稳定性产生威胁。

（3）降雨对多年冻土斜坡浅层稳定性产生的不利影响较明显，雨水入渗一段时间后斜坡稳定性迅速下降，且雨水入渗对活动层水热力的演化产生持续影响，安全系数最小值滞后于降雨数天。

（4）极限状态下斜坡位移演化规律印证了冻土浅层滑坡是活动层沿多年冻土层滑脱的过程，土体位移变化与雨水入渗形成的渗流场关系密切，低强度的持续降雨导致坡体中下部变形最显著。

参考文献

Wu Qingbai， Zhang Zhongqiong， Liu Ge.

Relationships between climate warming and engineering stability of permafrost on Qinghai-Tibet plateau

［J］. Journal of Engineering Geology， 2021， 29（2）： 342-352.

[本文引用: 1]

吴青柏， 张中琼， 刘戈.

青藏高原气候转暖与冻土工程的关系

［J］. 工程地质学报， 2021， 29（2）： 342-352.

[本文引用: 1]

Luo Jing， Niu Fujun， Lin Zhanju， et al.

Development of thawing hazards and thermal influence on permafrost along Qinghai-Tibet engineering corridor

［J］. Journal of Engineering Geology， 2014， 22（2）： 326-333.

[本文引用: 1]

罗京， 牛富俊， 林战举， 等.

青藏工程走廊典型热融灾害现象及其热影响研究

［J］. 工程地质学报， 2014， 22（2）： 326-333.

[本文引用: 1]

Kuang Xingxing， Jiao J J.

Review on climate change on the Tibetan Plateau during the last half century

［J］. Journal of Geophysical Research： Atmospheres， 2016， 121（8）： 3979-4007.

[本文引用: 1]

Yubonchit S， Chinkulkijniwat A， Horpibulsuk S， et al.

Influence factors involving rainfall-induced shallow slope failure： Numerical study

［J］. International Journal of Geomechanics， 2017， 17（7）： 04016158.

[本文引用: 1]

Goodfellow B W， Boelhouwers J.

7.31 hillslope processes in cold environments： An illustration of high-latitude mountain and hillslope processes and forms

［M］//Treatise on Geomorphology. Amsterdam： Elsevier， 2013： 320-336.

[本文引用: 1]

Huscroft C， Lipovsky P， Bond J.

Permafrost and landslide activity： Case studies from southwestern Yukon territory

［J］. Yukon exploration and geology， 2004.

[本文引用: 1]

Yin Guo’an， Niu Fujun， Lin Zhanju， et al.

The distribution characteristics of permafrost along the Qinghai-Tibet Railway and their response to environmental change

［J］. Journal of Glaciology and Geocryology， 2014， 36（4）： 772-781.

[本文引用: 1]

尹国安， 牛富俊， 林战举， 等.

青藏铁路沿线多年冻土分布特征及其对环境变化的响应

［J］. 冰川冻土， 2014， 36（4）： 772-781.

[本文引用: 1]

Lewkowicz A G， Harris C.

Frequency and magnitude of active-layer detachment failures in discontinuous and continuous permafrost， northern Canada

［J］. Permafrost and Periglacial Processes， 2005， 16（1）： 115-130.

[本文引用: 2]

Vedie E， Lagarde J L， Font M.

Physical modelling of rainfall- and snowmelt-induced erosion of stony slope underlain by permafrost

［J］. Earth Surface Processes and Landforms， 2011， 36（3）： 395-407.

[本文引用: 1]

Duan Dongming， Shen Yupeng， Zhang Luxin.

Stability analysis of Qinghai-Tibet Railway slope embankment in the event of rainstorm at the southern boundary region of the permafrost zone

［J］. Journal of Engineering Geology， 2008， 16（2）： 250-254.

[本文引用: 1]

段东明， 沈宇鹏， 张鲁新.

暖季强降雨对多年冻土南界斜坡路基稳定性影响分析

［J］. 工程地质学报， 2008， 16（2）： 250-254.

[本文引用: 1]

Niu Fujun， Luo Jing， Lin Zhanju， et al.

Thaw-induced slope failures and stability analyses in permafrost regions of the Qinghai-Tibet Plateau， China

［J］. Landslides， 2016， 13（1）： 55-65.

[本文引用: 1]

Zhao Wen， Wang Xiaojing， Xu Zhengxuan， et al.

Experimental study on natural evolution characteristics of coarse-grained soil slope in seasonal frozen soil region

［J］. Journal of Engineering Geology， 2021， 29（5）： 1497-1506.

[本文引用: 1]

赵文， 汪小静， 徐正宣， 等.

季节性冻土区粗颗粒土边坡自然演化特征试验研究

［J］. 工程地质学报， 2021， 29（5）： 1497-1506.

[本文引用: 1]

Zhang Mingli， Wen Zhi， Dong Jianhua， et al.

Mechanism of climate warming on thermal-moisture dynamics of active permafrost layer considering effect of rainfall

［J］. Rock and Soil Mechanics， 2019， 40（5）： 1983-1992.

[本文引用: 1]

张明礼， 温智， 董建华， 等.

考虑降雨作用的气温升高对多年冻土活动层水热影响机制

［J］. 岩土力学， 2019， 40（5）： 1983-1992.

[本文引用: 1]

Zhang Mingli， Wen Zhi， Xue Ke， et al.

A coupled model for liquid water， water vapor and heat transport of saturated–unsaturated soil in cold regions： Model formulation and verification

［J］. Environmental Earth Sciences， 2016， 75（8）： 1-19.

[本文引用: 1]

Li Desheng， Zhang Mingli， Wen Zhi， et al.

Quantitative evaluation of characteristic precipitation event （CPE）’s effect on hydro-thermal state of active layer in Beiluhe area

［J］. Natural Hazards， 2021， 105（2）： 2307-2326.

[本文引用: 1]

Fu Jianxin， Cao Guangchao， Li Lingqin， et al.

Spatiotemporal variation of air temperature in the middle and eastern parts of the Qilian Mountains and the nearby regions during the period of 1960—2014

［J］. Arid Zone Research， 2018， 35（3）： 549-561.

[本文引用: 1]

付建新， 曹广超， 李玲琴， 等.

1960—2014年祁连山中东段及其附近地区气温时空变化特征

［J］. 干旱区研究， 2018， 35（3）： 549-561.

[本文引用: 1]

Xiao Liangui， Shi Mingzhang， Yuxian La， et al.

Characteristics of precipitation change in Qilian Mountains during 1961—2019

［J］. Qinghai Science and Technology， 2020， 27（5）： 76-80.

[本文引用: 1]

肖莲桂， 石明章， 喇玉先， 等.

1961—2019年祁连山地区降水变化特征分析

［J］. 青海科技， 2020， 27（5）： 76-80.

[本文引用: 1]

Chen Guoqing， Huang Runqiu， Shi Yuchuan， et al.

Stability analysis of slope based on dynamic and whole strength reduction methods

［J］. Chinese Journal of Rock Mechanics and Engineering， 2014， 33（2）： 243-256.

[本文引用: 1]

陈国庆， 黄润秋， 石豫川， 等.

基于动态和整体强度折减法的边坡稳定性分析

［J］. 岩石力学与工程学报， 2014， 33（2）： 243-256.

[本文引用: 1]

Chen Xiaobai， Liu Jiankun， Liu Hongxu. Frost action of soil and foundation engineering［M］. Beijing： Science Press， 2006.

[本文引用: 2]

陈肖柏， 刘建坤， 刘鸿绪. 土的冻结作用与地基［M］. 北京： 科学出版社， 2006.

[本文引用: 2]

Ma Wei， Wang Dayan. Mechanics of frozen soil［M］. Beijing： Science Press， 2014.

[本文引用: 1]

马巍， 王大雁. 冻土力学［M］. 北京： 科学出版社， 2014.

[本文引用: 1]

Li Shuangyang， Niu Fujun， Sun Zhizhong， et al.

Numerical study on moisture-heat-mechanics coupled process of Jiangluling Tunnel construction in permafrost region

［J］. Journal of Glaciology and Geocryology， 2018， 40（5）： 966-973.

[本文引用: 1]

李双洋， 牛富俊， 孙志忠， 等.

多年冻土区姜路岭隧道施工水热力数值研究

［J］. 冰川冻土， 2018， 40（5）： 966-973.

[本文引用: 1]

Zhang Yuzhi， Du Yanliang， Sun Baochen.

Temperature distribution in roadbed of high-speed railway in seasonally frozen regions

［J］. Chinese Journal of Rock Mechanics and Engineering， 2014， 33（6）： 1286-1296.

[本文引用: 1]

张玉芝， 杜彦良， 孙宝臣.

季节性冻土地区高速铁路路基地温分布规律研究

［J］. 岩石力学与工程学报， 2014， 33（6）： 1286-1296.

[本文引用: 1]

Liu Weimin， He Ping， Zhang Zhao.

A calculation method of thermal conductivity of soils

［J］. Journal of Glaciolgy and Geocryology， 2002， 24（6）： 770-773.

[本文引用: 1]

刘为民， 何平， 张钊.

土体导热系数的评价与计算

［J］. 冰川冻土， 2002， 24（6）： 770-773.

[本文引用: 1]

Li Shuangyang， Lai Yuanming， Zhang Minyi， et al.

Study on long-term stability of Qinghai-Tibet Railway embankment

［J］. Cold Regions Science and Technology， 2009， 57（2/3）： 139-147.

[本文引用: 1]

van Genuchten M T.

A closed-form equation for predicting the hydraulic conductivity of unsaturated soils

［J］. Soil Science Society of America Journal， 1980， 44（5）： 892-898.

[本文引用: 1]

Lu Ning， Likos W J. Unsaturated soil mechanics［M］. Wei Changfu， Hou Long， Jian Wenxing， trans. Beijing： Higher Education Press， 2012.

[本文引用: 1]

卢宁， Likos William J.. 非饱和土力学［M］. 韦昌富， 侯龙， 简文星， 译. 北京： 高等教育出版社， 2012.

[本文引用: 1]

Zhang Mingli， Wen Zhi， Xue Ke， et al.

Influence of slope seepage on thermal regime of Qinghai-Tibet highway in permafrost regions

［J］. Highway， 2016， 61（7）： 13-19.

[本文引用: 1]

张明礼， 温智， 薛珂， 等.

斜坡渗流对青藏公路多年冻土路基温度场的影响

［J］. 公路， 2016， 61（7）： 13-19.

[本文引用: 1]

Bai Ming， Niu Lijuan， Wei Rongni.

Temporal and spatial variation of air temperature in the middle Qilian Mountains in recent 58 years

［J］. Journal of Qinghai Meteorology， 2019（2）： 12-14.

[本文引用: 1]

白明， 牛丽娟， 魏荣妮.

近58年来祁连山中段气温时空变化

［J］. 青海气象， 2019（2）： 12-14.

[本文引用: 1]

Zhu Linnan.

Study of the adherent layer on different types of ground in permafrost regions on the Qinghai-Xizang plateau

［J］. Journal of Glaciology and Geocryology， 1988， 10（1）： 8-14.

[本文引用: 1]

朱林楠.

高原冻土区不同下垫面的附面层研究

［J］. 冰川冻土， 1988， 10（1）： 8-14.

[本文引用: 1]

Sun Zhe， Zhao Lin， Hu Guojie， et al.

Influence of lower boundary conditions on the numerical simulation of permafrost temperature field changes

［J］. Journal of Glaciology and Geocryology， 2021， 43（2）： 357-369.

[本文引用: 1]

孙哲， 赵林， 胡国杰， 等.

下边界条件对多年冻土温度场变化数值模拟的影响

［J］. 冰川冻土， 2021， 43（2）： 357-369.

[本文引用: 1]

Wu Huozhen， Feng Meiguo， Jiao Yuyong， et al.

Analysis of sliding mechanism of accumulation horizon landslide under rainfall condition

［J］. Rock and Soil Mechanics， 2010， 31（）： 324-329.

[本文引用: 1]

吴火珍， 冯美果， 焦玉勇， 等.

降雨条件下堆积层滑坡体滑动机制分析

［J］. 岩土力学， 2010， 31（）： 324-329.

[本文引用: 1]

Wang Baolin， Nichol S， Su Xueqing.

Geotechnical field observations of landslides in fine-grained permafrost soils in the Mackenzie Valley

［J］. CanadaLandslides， 2005： 203-212.

[本文引用: 1]

Jin Dewu， Sun Jianfeng， Fu Shaolan.

Discussion on landslides hazard mechanism of two kinds of low angle slope in permafrost region of Qinghai-Tibet Plateau

［J］. Rock and Soil Mechanics， 2005， 26（5）： 774-778.

[本文引用: 1]

靳德武， 孙剑锋， 付少兰.

青藏高原多年冻土区两类低角度滑坡灾害形成机理探讨

［J］. 岩土力学， 2005， 26（5）： 774-778.

[本文引用: 1]

Li Guangxin.

Reply to discussion on “On soil skeleton and seepage force”

［J］. Chinese Journal of Geotechnical Engineering， 2017， 39（2）： 387-388.

[本文引用: 1]

李广信.

对“论土骨架与渗透力”讨论的答复

［J］. 岩土工程学报， 2017， 39（2）： 387-388.

[本文引用: 1]

Lamoureux S F， Lafrenière M J.

Fluvial impact of extensive active layer detachments， cape bounty， Melville Island， Canada

［J］. Arctic， Antarctic， and Alpine Research， 2009， 41（1）： 59-68.

[本文引用: 1]

Harris C， Lewkowicz A G.

An analysis of the stability of thawing slopes， Ellesmere Island， Nunavut， Canada

［J］. Canadian Geotechnical Journal， 2000， 37（2）： 449-462.

[本文引用: 1]

/

〈 〉

陇ICP备2021001824号-10
　 甘公网安备 62010202000676号
网站版权 © 《冰川冻土》编辑部
地址：兰州市天水中路8号　
邮编：730000　电话：0931-8260767　
E-mail：edjgg@lzb.ac.cn
总访问量 　 今日访问 　 在线人数
技术支持：北京玛格泰克科技发展有限公司