库水位下降条件下堆积层滑坡稳定性分析

尹晓萌，晏鄂川，杜 毅，雷玉菊，谭朝瑞

(1.中国地质大学(武汉)工程学院，湖北 武汉 430074；2. 华中师范大学教育技术协同创新中心，湖北 武汉 430079)

库水位下降条件下堆积层滑坡稳定性分析

尹晓萌1，晏鄂川1，杜 毅1，雷玉菊2，谭朝瑞1

(1.中国地质大学(武汉)工程学院，湖北 武汉 430074；2. 华中师范大学教育技术协同创新中心，湖北 武汉 430079)

为研究库水下降过程中不同因素作用下的堆积层滑坡稳定性响应规律，基于块体极限平衡理论用解析法分析地下水浸润线下降时，影响滑坡稳定性变化的主要因素。概化典型堆积层滑坡地质模型，采用数值模拟手段获取在各因素变化时，库水下降条件下不同时刻滑坡稳定系数，统计分析其变化规律。研究发现，库水下降时，滑坡稳定性变化受滑体渗透系数K、地下水浸润线变动区内的下滑带倾角α和阻滑段与下滑段体积比ε、滑带内摩擦角φ等影响，表现出逐步下降型、逐步上升型、“U”字型、阶梯上升型四类变化规律；K、ε控制稳定系数大小及其曲线形态，α和φ仅影响稳定系数值和不同时刻的稳定系数差异；库水下降始末滑坡的稳定系数变化百分比Q与渗透系数负对数-lnk在逐渐下降型曲线中表现出较好的幂函数关系，ε、α对幂函数的系数和指数的影响呈现出较好的规律性。

库水位下降；堆积层滑坡；稳定性变化规律；数值模拟

0 引言

库水位下降对滑坡稳定性的影响是当前滑坡研究领域十分重要的课题[1]，鉴于单体滑坡内部发育特征和赋存环境的差异性，众多学者从不同因素着手对库水下降条件下的滑坡稳定性变化规律进行了研究。刘新喜[2]、向 玲[3]、宋 琨[4]、梁学战[5]等结合具体实例，研究了滑坡的稳定性对库水位下降速率及滑体渗透系数的响应规律。刘晓红[6]研究了滑坡长度、厚度等对库水位下降条件下的稳定性变化规律的影响。武秀文[7]研究了凹形坡、凸形坡及其组合形态的滑坡稳定性变化规律。时卫民等[8]采用解析法分析库水下降速度V、坡体渗透系数K、给水度L以及含水层厚度h在库水位下降过程对稳定性影响。刘礼领[9]从物质结构的角度，研究了底部分布砂卵石层的万州中学滑坡在库水下降条件下的稳定性变化规律。虽然前人对库水下降条件下滑坡的稳定性问题取得了较丰富的成果，但多结合单体工程进行研究，缺乏高度概括性和广泛适用性，此外不同单体滑坡研究获得的稳定性变化趋势各异，导致对稳定性变化规律缺乏系统性的认识。

基于此，本文采用传力系数法，建立库水下降过程中相邻时刻滑块的剩余下滑力变化公式，分析影响稳定性变化的主要因素，在此基础上，利用前人研究成果概化典型堆积层滑坡地质模型，采用数值模拟手段研究滑坡在各因素下的稳定性变化规律。

1 影响滑坡稳定性的主要因素

库水位降落时，滑体中的地下水位面发生变动。针对地下水位变动条件下的单个条块，采用传力系数法计算其剩余下滑力。选取滑块中的土骨架作为研究对象，计算水的作用力时，考虑静水浮力与动水渗透力两种力[10]。地下水作用下的条块受力见图1。

图1 地下水下降条件下滑块的受力图

在i时刻，条块的下滑力和抗滑力分别表示如下：

Ti=Wisinα+Picos(α-βi)

(1)

Ri=Wicosαtanφ+Pisin(α-βi)tanφ+cl

(2)

其中：Wi=Vuiγ+Vdiγ′

(3)

Pi=γwVdisinβi

(4)

式中：Ti、Ri、Wi、Pi——分别为i时刻滑块的下滑力、抗滑力、重力和渗透力；

c、φ——分别为滑面的粘聚力和内摩擦角；

l、α——滑块底面长度和倾角；

βi——i时刻浸润线倾角；

Vui、Vdi——i时刻浸润线以上和以下滑块体积；

γ、γ′、γw——土体天然重度、浮重度和水的重度。

为简化计算，令Vui=h1ib；Vdi=hwib，其中h1i，hwi分别为i时刻浸润线上、下部分滑块的厚度；b为滑块的宽度，则由式(1)-(2)可得i时刻剩余下滑力为：

Fi=b(h1iγ+hwiγ′)sinα+γwhwibsinβicos(α-βi)-b(h1iγ+hwiγ′)cosαtanφ-γwhwibsinβisin(α-βi)tanφ-cl

(3)

同理可得i+1时刻剩余下滑力为：

Fi+1=b(h1i+1γ+hwi+1γ′)sinα+γwhwi+1bsinβi+1cos(α-βi+1)-b(h1i+1γ+hwi+1γ′)cosαtanφ-γwhwi+1bsinβi+1sin(α-βi+1)tanφ-cl

(4)

则(3)-(4)得相邻时刻剩余下滑力变化量：

ΔF=b(γ-γ′)Δh(cosαtanφ-sinα)+γwhwibsinβ[cos(α-βi)- sin(α-βi)tanφ]-γwhwi+1bsinβi+1[cos(α-βi+1)-sin(α-βi+1)tanφ]

(5)

ΔF=bCγw{hwsin(βi-βi+1)cos[(βi+βi+1)-(φ+α)]+0.5Δhcos(βi-βi+1)sin[(βi+βi+1)-(φ+α)]+Δh[0.5sin(φ+α)+0.8sin(φ-α)]}

(6)

式中，hw代表滑块的浸水高度，当水位与滑块相对位置固定时，与滑面倾角α有关，Δh、βi、βi+1与滑体渗透系数K和地表水下降速率ν相关，因此，当简化浸润线上下部滑体体积，并设置土体天然重度、浮重度与水的浮重度呈固定比时，可得影响ΔF的因素有内摩擦角φ、滑面倾角α、滑体渗透系数K、地表水下降速率ν等。因滑块剩余下滑力变化ΔF与滑坡稳定系数变化ΔFs有较强相关性，可利用式(6)间接分析ΔFs的影响因素，则有ΔFs=F(K，ν，α，φ)。

2 滑坡稳定性变化规律

2.1三峡库区堆积层滑坡基本特征

堆积层滑坡通常发生在第四系及近代松散堆积层中，这类滑坡在三峡库区中分布较广泛。按成因不同，主要分为滑坡堆积体滑坡、崩塌堆积体滑坡、崩滑堆积体滑坡、残坡积层滑坡四类。滑体物质往往由崩积物、崩坡积物、残坡积物及冲洪积与崩坡积混合物堆积而成，成分常为(粉质)粘土夹碎块石、含碎石(粉质)粘土、碎块石土。因物质成分及结构、形成时间的不同，滑体渗透系数差异较大，常呈弱-中等透水性。按滑体厚度，库区堆积层滑坡由浅层 (<10 m)到超深层 (>50 m)皆有分布，主要发育中(10～25 m)、深层(25～50 m)滑坡。堆积层滑坡的坡面形态不一，常见凹形、凸形、凹凸复合形、直线形几类[12]，据不完全统计，库区堆积层滑坡中凸形坡最为发育，其次为直线形坡和凹形坡[13]。

堆积层滑坡的滑带位置与滑坡物质组合、结构特征、地下水的分布、地质构造等密切相关，多位于堆积体和基岩接触面附近，滑带土一般为含砾(粉质)粘土，由于该界面处含泥量增加、渗透性减小，常构成斜坡的相对不透水层。滑动面形态多样，典型滑动面由较平缓的前缘、较陡的中部和陡后缘三段组成。

三峡库区堆积层滑坡主要发生在软岩中，这套软岩类地层常见的有泥岩、粉砂质泥岩、页岩、泥页岩、泥质粉砂岩、泥灰岩等。

大量堆积层滑坡的前缘位于175 m最高库水位以下，库水位的周期性升降导致滑坡渗流场发生变化，引起应力场不断调整，进而对滑坡稳定性产生较大影响。

2.2滑坡基本地质模型建立

三峡库区堆积层滑坡发育形态复杂，唐胜传[14]等根据调查，依据下滑段滑带倾角不同，将库区滑坡分为平缓滑坡、中倾角滑坡、陡倾角滑坡，徐平[15]进一步总结出平缓滑坡、中倾角滑坡发育的形态规律，于永贵[16]研究了堆积层滑坡发育的坡角、坡高、滑体厚度等地质规律。基于前述三峡库区堆积层滑坡的基本特征，结合前人的研究成果，选取分布广泛的中倾角滑坡(下滑段坡度16°～45°)建立堆积层滑坡基本地质模型见图2。该模型滑带分为后部陡倾段，坡度约40°；中部下滑段，坡度18°～22°；前部阻滑段，坡度小于5°，地表坡度20°～32°。

图2 典型滑坡地质模型及其与库水相对位置变化示意图

2.3影响滑坡稳定性的因素设定

为了分析各因素对滑坡稳定性的作用规律，采用如下方式变换各影响因素：

(1)地下水作用范围内滑带倾角影响。从若干滑块组成的整个滑坡体而言，对滑带倾角效应的考虑应包括不同滑段处滑块的数量累积效应和倾角的绝对大小影响两种情况。前者指库水影响范围内下滑段和阻滑段的比例，即ε=阻滑段体积/下滑段体积，可通过滑坡体和库水位的相对位置的改变加以反映(图2)，若固定库水高低水位分别为175 m、145 m，则滑坡前后缘高程发生变化。根据文献[13]中的滑坡分布高程统计规律，建立ε=0、0.37、0.61、0.67、1.04、2.37、∞七种基本数值模型。在基本模型二(ε=0.37)基础上，调整库水影响范围内下滑段滑带倾角α，建立α=18°、21°、23°、25°、28°五种数值模型，分析倾角绝对大小对滑坡稳定性的影响。

(2)渗透系数和库水下降速率的影响。按照库水调度数据，固定库水下降速率为1 m/d，即不考虑库水下降速率的变化对滑坡稳定性的影响。已有研究发现，三峡库区堆积层滑体渗透系数普遍介于 1.0×10-4～1.0 × 10-6m / s[2]，为充分研究渗透系数对滑坡稳定性的影响，将滑体渗透系数K划分为7级，由大到小分别为1e-3、5e-4、1e-4、5e-5、1e-5、5e-6、1e-6 m/s。

(3)内摩擦角φ影响。根据前人统计滑带发育的普遍规律，分别设置为φ=20°、25°、30°、35°进行研究。

根据基本地质模型和影响因素的设定方式，建立各数值计算模型。采用Geo-studio渗流分析软件获取库水降落过程中滑体内地下水渗流场的变化特征，基于渗流分析结果，利用Geo-slope稳定性计算软件，分析库水下降条件下滑坡稳定性的响应规律。图3为滑坡数值计算典型模型和渗流计算结果，渗流场模拟前分别设置地表水头边界条件和滑带不透水边界，由渗流场计算结果可见，库水位降落时，受滑体物质透水能力影响，滑体中的地下水位相对地表水滞后，产生一定水头差，形成指向坡体外的动水压力。

图3 滑坡典型数值计算模型和渗流计算结果

2.4数值计算结果分析

2.4.1稳定性曲线类型变化规律

统计数值计算结果，发现随各因素变化，库水下降过程中的滑坡稳定性表现出4种变化类型，分别为逐渐下降型(A型，图4a)、逐渐上升型(C型，图4c)、U型(B型，图4b)、阶梯上升型(D型，图4d)。其中U型分为左肩较高(B1型)、右肩较高两种类型(B2型)。

图4 稳定性变化曲线类型

图5为内摩擦角φ=20°、下滑段滑带倾角α=18°固定，阻滑/下滑段体积比ε不同的7种基本模型，随滑体的渗透系数K变化时的稳定性变化曲线类型的统计，图4中显示，K和ε对稳定性变化曲线类型有重要影响，当ε和K增大时，曲线由下降型逐渐向上升型过渡，曲线类型皆按A-B1-B2-C-D方向变化，但随ε或K在经验值范围内变化，并非都表现出完整的变化过程。

图5 受渗透系数K和阻滑/下滑段体积比ε影响的稳定性变化曲线类型统计

ε=0、∞时，无论K值如何变化，稳定性变化曲线类型皆不发生改变，表明阻滑/下滑段体积比ε比渗透系数K值对稳定系数变化更具控制性。ε=0时，处于地下水变动带内的全部为滑坡的下滑段，稳定性变化曲线表现为逐渐下降型，ε=∞时，处于地下水变动带内的皆为滑坡的阻滑段，稳定性变化曲线呈阶梯上升型。可见，库水影响范围的滑带倾角较大和较小两种情况下，水位下降对滑坡稳定性表现出相反影响：滑面倾角α较大时，水位下降对稳定性产生不利影响，滑面倾角α较小时，水位下降增大了滑坡稳定性。

水对滑体的作用力有两种，一是渗透力，另一个是水的浮力[10]。当库水下降时，地下水渗透力增大而水的浮力减小，滑坡常受渗透力沿滑面方向的下滑分力影响，稳定性有减小趋势，而浮力对滑坡稳定性作用因滑面倾角而异。α较大时，若不考虑水的影响，滑体天然重力引起的下滑力大于抗滑力，对稳定性不利，当叠加库水下降环境时，因水的浮力减小使得重力对滑面有效作用力增大，导致下滑力和抗滑力差异增大，增强了滑坡的不稳定，因此在水的浮力和渗透力综合影响下，滑坡稳定系数减小。很多库水下降型滑坡都遵循这种规律，在文献[17]、文献[18]中的树坪滑坡、白水河滑坡等的变形监测曲线也反映出这类滑坡在库水下降条件下的稳定性变化趋势。α较小时，随地下水下降，滑体有效重力增大，受小角度滑面倾角影响，抗滑力提高，有助于滑坡稳定，且这种对稳定性的增强作用相对渗透力对稳定性的弱化作用更加显著，导致滑坡的稳定系数整体趋于增大。对这类地下水影响范围内滑面倾角平缓的滑坡，殷跃平[19]、赵代鹏等[20]分别采用解析法、数值模拟与物理试验结合法得到其稳定性在库水下降条件下的增强趋势，一些学者也将这类滑坡称之为浮拖减重型滑坡。由此可见，地下水影响范围内的滑带倾角对滑坡稳定性起首控作用，较大的渗透力并非必然引起滑坡稳定性降低，在实际工作中，不可过于夸大渗透力对稳定性的影响。

除ε=0、ε=∞外，K→1e-3时，稳定性变化曲线类型普遍以上升型为主，K→1e-6时，则以下降型为主，该规律与文献[5]的结论有较好的一致性。这是因为K越小，地下水滞后地表水形成的水力梯度增大，产生的渗透力变大，沿滑面的渗透分力引起下滑段内滑块剩余下滑力显著增大，对滑坡整体稳定性产生控制性影响，导致稳定系数减小。K越大，地下水排泄能力强，水头差降低，渗透力减小，此时水的浮托力对滑坡的稳定性起主导作用，随地下水位下降，浮托力减小，导致阻滑段剩余下滑力显著减小，滑坡整体稳定性趋于提高。由各模型的渗流计算结果发现，K=1e-3时，随库水位降落，地下水浸润线近于平行下降，此时渗透力趋于0，滑体仅受水的浮力影响。

由B1型向B2型过渡时，必然经历了库水下降始末稳定系数值相等的变化类型B3，图5中曲线M为对应K和ε下的B3类的连线，以曲线M为界，图框左下部稳定性变化曲线以下降型为主，含A型和B1型，右上部以上升型为主，含B2型、C型和D型。随ε增大，M曲线呈近似对数型递减，ε越大，对应B3类的渗透系数越小，表明地下水影响范围内的阻滑段比例越大时，水的渗透作用效应趋于弱化，由稳定性上升型转化为下降型所需的渗透力越大。由ε=0.37到ε=2.37，B3型最小稳定系数相对最大值的降低百分比分别为1.24%、2.09%、1.22%、0.69%、0.18%。可见B3型最大最小稳定系数差异比率在M曲线上呈现两侧低中间高的特征，最大值位于ε=0.37～0.66，K=5e-4～5e-5 m/s范围内。

2.4.2稳定系数变化规律

2.4.2.1渗透系数和阻/下滑体积比影响

统计发现，各时刻的稳定系数值及相对差异随渗透系数K和阻滑/下滑段体积比ε发生变化，表现为：

(1)ε不变时，随K减小，同一时刻稳定系数值皆减小，不同时刻Fs差异值的变化规律因稳定性变化曲线类型而异：逐渐上升型曲线，库水下降始末稳定系数增大百分比Q随K值增大而增大；逐渐下降型曲线，库水下降始末稳定系数降低百分比Q表现出相反规律。

(2)K不变，同一时刻稳定系数值与ε呈正相关，而不同类型曲线的始末稳定系数差异百分比Q(除U型外)皆与ε呈反相关。

由下降型曲线形态可见，该类曲线稳定系数降低段，在相邻时刻的渗透系数降低幅度随库水位下降逐渐减小。其中逐渐下降型曲线，随渗透系数K、阻滑/下滑段体积比ε变化时，库水下降末的稳定系数相对初始稳定系数的降低百分比Q表现出较好的规律性。

图6为渗透系数负对数-lnk与稳定系数的降低百分比Q的关系，可见：对应不同ε值，Q皆与-lnk呈较好的幂函数关系，即随渗透系数K值减小，稳定系数的降低百分比Q增大，且增大幅度有快速增加趋势；K不变时，随ε减小，稳定系数的降低百分比Q增大。表明，较低的渗透系数和较高的下滑段比例对稳定系数降低作用更加显著。

图6 稳定系数下降百分比Q与渗透系数K的关系

由拟合的幂函数y=axb发现，系数值a与ε呈正相关，而幂指数b随ε变化表现出负相关性，反映出当K和ε同时变化时，相对于高ε条件，低ε条件下的Q对K值敏感性更大，即相同K值范围内，稳定系数下降百分比Q的变化率随ε减小逐渐增大。由图7可见，-lna、b与ε皆呈较好的指数关系。由此，当渗透系数K和阻/下滑体积比ε发生变化时，逐渐下降型曲线的稳定系数降低百分比可用以下函数表示：

Q=e0.9881e2.3374ε-12.7369×(-lnk)4.6260-0.6197e2.2199ε

图7 幂函数系数和指数与阻滑/下滑段体积比ε的关系

2.4.2.2滑带倾角影响

改变基本模型2(ε=0.37)中的地下水影响范围内下滑段的滑带倾角α，建立不同滑带倾角的滑坡数值模型，计算不同模型在设定渗透系数值下的滑坡稳定性变化。结果发现，同一渗透系数不同模型间基本呈现相同的稳定性曲线类型，但稳定系数值表现出明显变化：渗透系数不变时，随α增大，不同模型相同时刻的稳定系数值逐渐减小。

五种数值模型中，在滑带倾角α变化的同时，阻/下滑体积比ε也发生改变，在α、ε和K三因素影响下，逐渐下降型曲线的稳定系数降低百分比Q仍然表现出较好的规律性。图中显示，Q仍与-lnk呈较好的幂函数关系，与ε和K两因素影响下的Q值呈现的规律不同的是，在三因素综合影响下，随ε增加，拟合幂函数y=axb的b值呈递增而非递减变化，由此反映出滑带倾角α的变化效应，即α递增时，稳定系数下降百分比Q在相同K值范围内的变化幅度逐渐增大，与ε的作用效应恰好相反，受该因素影响，当渗透系数K值不变时，不同的ε值对应的Q值间的差异被弱化。基于图7、图8中的函数关系，分析α与y=axb中a、b关系，由图9可见-lna、b与α呈较好的正线性关系。由此，可建立三因素影响下的稳定系数下降百分比经验公式：

Q=e0.9881e2.3374ε-0.5233α-1.6289×(-lnk)-0.6197e2.2199ε+0.2395α+0.1977

图8 稳定系数下降百分比Q与渗透系数K的关系

图9 幂函数系数和指数与滑带倾角α的关系

建立α=30°滑坡模型，将模拟获取的稳定系数百分比值与经验公式计算所得的值对比，由表1可见，不同渗透系数条件下，经验值普遍比模拟所得值大8%～13%，表明该公式有较高的可靠性。

表1 稳定系数下降百分比Q的模拟值和经验拟合值对比Table 1 The contrast of simulation value and experiencefitting value of decline percentage Q

2.4.2.3内摩擦角影响

改变基本数值模型中的滑带内摩擦角φ，根据数值计算结果分析其对滑坡稳定性变化规律的影响。结果显示，φ值并不改变稳定性变化曲线类型，仅影响稳定系数值的大小和不同时刻该值间的相对差异。随φ增大，同一时刻的稳定系数值逐渐增大，不同类型曲线的稳定系数相对差异对φ的响应不同：φ增大时，逐渐下降型曲线，库水下降始末稳定系数下降百分比逐渐减小；而逐渐上升型曲线，始末增大百分比表现出相反规律。

3 结论

(1)影响库水下降滑坡稳定性变化的主要因素有渗透系数、滑带倾角、滑带内摩擦角。在各因素影响下，滑坡稳定性变化曲线类型分为逐渐下降型(A型)、逐渐上升型(C型)、U型(B型)、阶梯上升型(D型)，其中U型曲线根据对称性分为左肩较高的B1、右肩较高的B2和左右等高的B3型。

(2)阻滑体积与下滑体积之比ε是控制稳定性曲线变化类型的首要因素，其次为渗透系数K。随ε和K增大，曲线由下降型逐渐向上升型过渡，曲线类型皆有A-B1-B2-C-D变化的趋势。连接各对应ε和K值下的B3型组成曲线M，则 B3型的最大、最小稳定系数差异百分比在M曲线上呈现两侧低中间高的特征。

(3)各因素影响下，逐渐下降型曲线呈现较好的规律性。库水下降始末稳定性下降百分比Q与渗透系数K的函数关系可表示为Q=a(-lnk)b。其中-lna、b随ε增大呈指数型递减，随滑带倾角α增大呈线性递增。由此，建立了渗透系数K、阻滑体积与下滑体积之比ε、滑带倾角α等三因素影响下的稳定系数下降百分比经验公式。φ增大时，库水下降始末稳定系数下降百分比逐渐减小。

(4)ε不变时，随渗透系数K增大，同一时刻稳定系数值和逐渐上升型曲线中库水下降始末稳定系数增大百分比Q皆增大；K不变时，同一时刻稳定系数值与ε、φ呈正相关，与α呈反相关，逐渐上升型曲线中，始末增大百分比Q随ε增大而减小，随φ增大而增大。

[1] 郭飞, 付调金, 阮荣乾, 等. 基于 ABAQUS 强度折减法分析库水位下降对边坡稳定性影响[J]. 三峡大学学报: 自然科学版, 2012, 24(3): 15-19.

GUO Fei, FU Tiaojin, RUAN Rongqian, et al. Analysis of influence of drawdown of reservoir water level on landslide stability using strength reduction method based on abaqus[J]. J of China Three Gorges Univ.(Natural Sciences), 2012, 24(3): 15-19.

[2] 刘新喜.库水位下降对滑坡稳定性的影响及工程应用研究[D].武汉: 中国地质大学, 2003.

LIU Xinxi. Research on landslide stability in drawdown of reservoir water level and its application[D].Wuhan: China University of Geosciences, 2003.

[3] 向玲, 王世梅, 王力. 动水压力型滑坡对库水位升降作用的响应——以三峡库区树坪滑坡为例[J]. 工程地质学报, 2014, 22(5): 876-882.

XIANG Ling, WANG Shimei, WANG Li. Response of typical hydrodynamic pressure landslide to reservoir water level fluctuation: using landslide in Three Gorges Reservoir as an example[J]. Journal of Engineering Geology, 2014, 22(5): 876-882.

[4] 宋琨, 晏鄂川, 朱大鹏, 等. 基于滑体渗透性与库水变动的滑坡稳定性变化规律研究[J]. 岩土力学, 2011, 32(9): 1510-1516.

SONG Kun, YAN Echuan, ZHU Dapeng, et al. Base on permeability of landslide and reservoir water change to research variational regularity of landslide stability[J]. Rock and Soil Mechanics, 2011, 32(9): 1510-1516.

[5] 梁学战, 陈洪凯. 库水位升降条件下不同渗透性的滑坡体稳定性变化规律[J]. 中国地质灾害与防治学报, 2012, 23(4): 20-26.

LIANG Xuezhan, CHEN Hongkai. Stability variation of landslide with different permeability coefficient following reservoir water level fluctuation [J]. The Chinese Journal of Geological Hazard and Control, 2012, 23(4): 20-26.

[6] 刘晓红.水库型滑坡库水位动力加载与位移响应特征及其稳定性演化规律研究[D].山东：青岛理工大学, 2013.

LIU Xiaohong. Study on the dynamic increment displacement response ratio of the landslide induced by trservoir water and its stability evolution rules[D]. Shandong: Qingdao Technological University, 2013.

[7] 武秀文.库水位变化加卸载动力效应及其对堆积层边坡稳定性釤响规律研究[D].山东：青岛理工大学, 2012.

WU Xiuwen. The study of the load-unload dynamic impacts of reservoir water level changes and its effect on debris landslide stability[D]. Shandong: Qingdao Technological University, 2012.

[8] 时卫民, 郑颖人. 库水位下降情况下滑坡的稳定性分析[J]. 水利学报, 2004, 35(3): 76-80.

SHI Weimin, ZHENG Yingren. Analysis on stability of landslide during reservoir drawdown[J]. Journal of Hydraulic Engineering, 2004, 35(3): 76-80.

[9] 刘礼领. 库水位升降对堆积体滑坡稳定性影响的研究[J]. 西部探矿工程, 2011, 23(2): 35-39.

LIU Liling. Study on the impact of water level fluctuation on deposit landslide stability[J]. West-china Exploration Engineering, 2011, 23(2): 35-39.

[10] 毛昶熙, 段祥宝. 关于渗流的力及其应用[J]. 岩土力学, 2009, 30(6): 1569-1574.

MAO Changxi，DUAN Xiangbao. On seepage forces and application[J]. Rock and Soil Mechanics, 2009, 30(6): 1569-1574.

[11] 汪洋, 刘艺梁. 浸润线缓慢上升对滑坡剩余推力的影响规律[J]. 灾害学, 2012, 27(4): 21-24.

WANG Yang , LIU Yiliang. Influence law of slowly rising of phreatic line on residual sliding force of land slides[J]. Journal of Catastrophology, 2012, 27(4): 21-24.

[12] 乔建平,吴彩燕,田宏岭. 三峡库区云阳-巫山段坡形因素对滑坡发育的贡献率研宄[J]. 工程地质学报,2006,14(1): 18-22.

QIAO Jianping,WU Caiyan,TIAN Hongling. Contributing ratios of the slope shape toward the land-slide development from Yunyang to Wushan in Three Gorges Reservoir area[J]. Journal of Engineering Geology,2006, 14(1): 18-22.

[13] 汤罗圣.三峡库区堆积成滑坡稳定性与预测预报研究[D].武汉: 中国地质大学, 2013.

TANG Luosheng. Research on stability and prediction for the colluvial landslide in the Three Gorges Reservoir[D].Wuhan: China University of Geosciences, 2013.

[14] 唐胜传, 柴贺军, 冯文凯. 三峡库区岸坡类型划分[J]. 公路交通技术, 2005 (5): 36-39.

TANG Shengchuan, CHAI Hejun, FENG Wenkai. Research on embankment slope classification of Three Gorges Reservoir areas[J]. Technology of Highway and Transport, 2005 (5): 36-39.

[15] 徐平. 三峡库区涉水滑坡体稳定性的可靠度研究[D].西安:长安大学,2011.

XU Ping. Reliability study on stability of submerged slide mass in Three-Gorges Reservoir area[D].Xi′an: Chang′an University, 2011.

[16] 于永贵.三峡库区万州-开县段堆积层滑坡发育规律及破坏模式研究[D].北京: 中国地质大学, 2008.

YU Yonggui. The development regularities and failure modes of colluvium landslides in Wanzhou and Kaixian zones of the Three Gorges Reservior[D].Beijing: China University of Geosciences, 2008.

[17] 向 玲,王世梅,王 力.动水压力型滑坡对库水位升降作用的响应——以三峡库区树坪滑坡为例[J]. 工程地质学报, 2014, 22(5): 876-882.

XIANG Ling, WANG Shimei, WANG Li. Response of typical hydrodynamic pressure landslide toreservoir water level fluctuation: shuping landslide in Three Gorges Reservoir as an example[J]. Journal of Engineering Geology, 2014, 22(5): 876-882.

[18] 高连通,易夏玮,李喜,等.三峡库区典型滑坡变形与高水位涨落关系研究[J]. 地质科技情报, 2011, 30(4):132-136.

GAO Liantong, YI Xiawei, LI Xi, et al. Case study of relationship between typical landslide distortion and water changes in Three Gorges area[J]. Geological Science and Technology Information, 2011, 30(4):132-136.

[19] 殷跃平. 三峡库区地下水渗透压力对滑坡稳定性影响研究[J]. 中国地质灾害与防治学报, 2003, 14(3):1-8.

YIN Yueping. Seepage pressure effect on landslide stability at the Three Gorges Reservoir area[J]. The Chinese Journal of Geological Hazard and Control , 2003, 14(3):1-8.

[20] 赵代鹏,王世梅,谈云志,等. 库水升降作用下浮托减重型滑坡稳定性研究[J]. 岩土力学, 2013, 34(4):1017-1024.

ZHAO Daipeng, WANG Shimei, TAN Yunzhi, et al. Stability studies of buoyancy weight loss landslides under reservoir water level fluctuation[J]. Rock and Soil Mechanics, 2013, 34(4):1017-1024.

Stabilityofcolluviallandslidewithdrawdownofreservoirwaterlevel

YIN Xiaomeng1，YAN Echuan1，DU Yi1，LEI Yuju2，TAN Zhaorui1

(1.FacultyofEngineering,ChinaUniversityofGeosciences(Wuhan),Wuhan,Hubei430074,China； 2.Collaborative&InnovativeCenterforEducationalTechnology,CentralChinaNormalUniversity,Wuhan,Hubei430079,China)

In order to study systematically landslide stability response law influenced by different factors in the process of water drop, based on limit equilibrium theory, the main factors when groundwater invasion line declines were analyzed by analytic method. Then, typical colluvial landslide geological model was generalized and numerical simulation was used to obtain different moment stability coefficient of landslide under the condition of reservoir water dropping. Change laws of landslide stability with different factors were analyzed. The results show that when the reservoir water level falls, the change of landslide stability is observably influenced by permeability coefficient K of sliding mass、internal friction angleφof sliding zone as well as dip angleαof slip plane and the ratio ε of stopping slide sect to main sliding sect in the changing area of groundwater infiltration line. The regularity of stability variation shows four types including gradual decline, gradual ascent, U-shaped, stair-step ascent. The permeability coefficientKand the ratioεof stopping slide sect to main sliding sect control the figure and changing shape of stability coefficients. The dip angleαand internal friction angleφof sliding zone only affect stability coefficient value and difference between different times. In gradually decline type, power function relationship lies in negative logarithm of permeability coefficient , that is -lnk, and the difference percentageQbetween stability coefficient of ending and beginning reservoir water. The coefficient and exponent of power function present a good regularity with the changing of the ratio ε of stopping slide sect to main sliding sect and dip angleαof slip plane.

drawdown of reservoir water level; colluvial landslide; regularity of stability variation; numerical simulation

TU 443

A

1003-8035(2017)03-0008-08

10.16031/j.cnki.issn.1003-8035.2017.03.02

2016-09-11;

2016-09-26

国家重点基础研究发展计划(973)项目(2011CB710605)

尹晓萌(1988-)，男，博士生，从事水文地质工程地质研究。E-mail: xiaomengyin123@163.com

晏鄂川(1969-)，男，博士，教授，从事水文地质工程地质研究。E-mail: yec880903@163.com