易连兴

(中国地质科学院岩溶地质研究所/自然资源部广西岩溶动力学重点实验室，广西 桂林 541004)

2010年6月28日14时，贵州省关岭县岗乌镇大寨村发生特大滑坡灾害[1]。大寨滑坡发生后，有关单位和专家学者对滑坡区进行了调查[1～3]和研究[4～9]。殷跃平等[4]研究了滑坡失稳后的运动特征和堆积特征，从地形、地质结构以及降雨作用等方面分析了滑坡的形成机制；刘传正[5]通过降雨、斜坡岩体结构两个方面的影响对滑坡的启动机理进行了分析，对大寨边坡岩体的稳定性进行力学平衡计算；孔纪名等[6]分析了滑坡区工程地质条件，认为边坡区岩性、构造、地形等构成了滑坡发生的必要条件，久旱强雨是最重要的激发致灾因素；张信宝[7]采用“超压密岩土边坡失稳理论”来解释大寨高速滑坡的启动机理，认为由于强降雨入渗，岩层裂隙与大气贯通，负压消失，“假”稳定状态消失，导致边坡失稳；张建江等[8]认为特殊的地形地貌、“上硬下软”的地质结构以及“前旱后雨”极端气候是大寨滑坡的成灾原因。上述研究侧重于滑坡区工程地质及滑坡体结构特征，提出了地下水对滑坡启动或诱导作用，但仅笼统提出基岩地下水和岩溶地下水的作用。胡光中[9]在大寨滑坡机理研究基础上，讨论了地下水的影响，并建立地下水作用下的滑坡模拟评价模型，但仅考虑了基岩裂隙水，没有考虑滑坡区后缘上覆灰岩区岩溶地下水的影响作用。

根据统计，80%以上滑坡的发生与地下水有密切关系[10～11]，因此，总结滑坡过程中的水-岩耦合影响模式是众多学者重点研究的内容之一。在这方面，非岩溶区开展了较多研究工作，如黄达等[12]、周剑等[13]、徐永强等[14]开展了库区水位升降耦合滑坡体内基岩裂隙水作用下的滑坡体稳定性评价；周平根[15]总结了土质、岩质滑坡8种水文地质结构类型，其中7种属于非岩溶区，岩溶区划分为管道含水体单一水文地质结构；李根等[16]系统总结了水岩耦合变形破坏过程及机理研究进展，并分析了因水岩耦合变形破坏机制认识不足而导致的惨痛教训案例。相对而言，针对岩溶石山区滑坡水-岩耦合作用的研究较少。基于岩溶含水介质结构和水动力场极其复杂，岩溶地下水影响滑坡模式有待深入研究。

本文通过野外调查，结合水化学和降水数据分析、综合研究等方法，分析和提出滑体区基岩裂隙水及后缘岩溶区地下水两种流场耦合对滑坡的控制作用，对研究大寨滑坡机理及岩溶地下水对滑坡的影响模式有一定意义。

1 研究区概况

大寨滑坡区及周边区域(图1)为溶蚀-侵蚀深切中低山河谷地貌。北盘江河床高程590 m，滑坡区后缘为山脊，最大高程1 695 m，切割深度大于1 000 m；相对光照水电站正常蓄水位730 m，高差也近千米。垂直北盘江发育东西向槽谷，槽谷内发育树枝状冲沟。冲沟地形上游为陡崖或陡坡，坡度25°～35°；下游为缓坡，坡度10°～15°。本次滑坡发生于永窝—大寨槽谷之冲沟上游。

研究区处于普安山字型构造脊柱、大田—法郎向斜北翼，为单斜构造，断裂不发育，断层主要由SW—NE向的小型断层组成，断层倾角60°～80°。出露地层有三叠系中统关岭组(T2g)、下统永宁镇组(T1yn)和夜郎组(T1y)，二叠系上统龙潭-长兴组(P3l+c)、峨眉山玄武岩(P3β)及中统茅口组(P2m)等；地层产状总体倾向SE，倾角28°～38°。

滑坡区及后缘土地利用类型为灌木、草地，滑坡区以下主要为旱地、村民房屋用地等。该区域以农业生产为主，没有采矿等工业活动。

图1 研究区水文地质简图Fig.1 Simplified hydrogeological map of the study area

2 复合水动力场及滑坡影响模式

2.1 含水岩组及其与滑坡的关系

(1)含水岩组

研究区分布有碳酸盐岩含水岩组(Ⅰ)和碎屑岩含水岩组(Ⅱ)。碳酸盐岩含水岩组包括茅口组(P2m)，永宁镇组1段(T1yn1)、3段和4段(T1yn3+4)，关岭组2段和3段(T2g2+3)等，岩溶中等至强发育，发育地下河(表1)、天窗和消水洞，受北盘江排泄基准面控制，岩溶地下水向北盘江或光照水电站库区径流排泄。

表1 地下河主要特征表

研究区内碎屑岩含水岩组包括龙潭-长兴组(P3l+c)、夜郎组(T1y)、永宁镇组2段(T1yn2)、关岭组1段(T2g1)等。碎屑岩从地面至下具有强风化、风化、弱风化特征，浅部节理和裂隙相对发育，深部相对不发育。碎屑岩基岩裂隙水主要受地形控制，顺坡朝最近冲沟或槽谷径流排泄。研究区主要含水岩组分布见图1，特征见表2。

(2)分水岭

永宁镇组2段碎屑岩(Ⅱ2区)含水岩组沿东西方向呈连续条带状分布，厚60～80 m；在西部水库旁公路边(图2a)、东侧上寨公路边(图2b)可见典型剖面，在滑坡后缘山体顶部d3处可见出露。该岩层表层风化后形成粉质黏土，如图2(b)左侧部分，右侧为永宁镇组1段(T1yn1)灰岩。该地层厚度稳定、裂隙不发育，为相对隔水层，其分布区域在地形上为山脊，即地表水分水岭，同时也是地下水分水岭，地表水和地下水分别向南侧下陇古槽谷地下河(gw20)、北侧地下河(gw21)以及滑坡区径流排泄。

(3)滑坡区地下水补径排特征

受分水岭隔水作用，分水岭以南的永宁镇组3段、4段(Ⅰ3区)岩溶管道水、关岭组1段基岩裂隙水和关岭组2段、3段岩溶地下水与滑坡区不存在水力联系。

茅口组灰岩(Ⅰ1区)与滑坡区位于分水岭同一侧，地貌上为岩溶槽谷，发育地下河gw22，主管道水位570～800 m。尽管该地层岩溶强发育，且与滑坡区存在250 m以上的水位差，但与滑坡区之间分布有隔水作用较好的玄武岩和龙潭组长兴组中下部，两者合计厚度达数百米，因此，茅口组灰岩(Ⅰ1区)不会对滑坡区地下水形成袭夺，即不存在地下水水力联系。当然，受地形控制，滑坡区排泄出地表的基岩裂隙水以溪流形式流入茅口组灰岩区，这方面与文章所讨论的地下水与滑坡的关系无关。

表2 研究区主要含水岩组

图2 永宁镇组2段典型照片Fig.2 Typical photos of T1yn2 near the Yongning Town

滑坡区(Ⅱ1区)为基岩裂隙水分布区，大气降水垂直入渗为主要补给源。永宁镇组1段岩溶区(Ⅰ2区)位于分水岭北侧，紧邻且上覆于滑坡区，岩溶管道裂隙水总体由东向西朝水库区径流排泄，在径流过程中，部分岩溶地下水侧向补给滑坡区碎屑岩(Ⅱ1区)，因此，Ⅰ2区与滑坡区有密切水力联系，对滑坡将产生直接影响。对类似岩溶山区，张楠等[18]也提出了降雨入渗且无有效的排泄通道，岩溶区的侧向补给对碎屑岩滑坡将起到重要作用。

滑坡区(Ⅱ1区)接受降水垂直入渗补给、侧向补给后，地下水顺地形向北径流，再转向沿永窝—大寨槽谷向水库区径流排泄，在径流过程中部分地下水就近在沟、槽中排泄。

2.2 复合水动力场及其对滑坡的影响

(1)特枯期水动力场

2009年8月15日以来，我国西南地区遭遇历史罕见的夏秋冬春四季连旱[19]，各地表层土体龟裂，泉点断流。研究区2009年9月至2010年4月合计降水量为108.0 mm，仅为同期多年平均月降水量416.8 mm的26%(表5)。通过访问，在此次特枯期，南坡下陇古槽谷中gw03泉未断流，水位还保持在gw03泉口高程1 432 m；北坡gw18和gw19两个泉断流。北坡永宁镇组1段(Ⅰ2区)岩溶水动力场在相对强径流条件作用下形成低水位径流槽，并袭夺两侧基岩裂隙水向水库径流和排泄，使部分夜郎组(Ⅱ1区)基岩裂隙水向岩溶Ⅰ2区径流并形成临时局部地下水分水岭，永宁镇组2段(Ⅱ2区)内水力坡度相对变陡(图2)。

(2)暴雨期水动力场及其对滑坡的影响

灾害发生前，滑坡区经历了关岭县有气象记录以来的最大一次降水过程，暴雨后的强烈补给，不同地段的地下水位急剧上升，形成与枯季不同的水动力场(图3)。永宁镇组2段(Ⅱ2区)内分水岭向北偏移；永宁镇组1段(Ⅰ2区)地下水位也随着强补给而急剧上涨，伴随强烈补给和受gw21排水能力限制，地下水在滑坡区后缘岩溶区富集，低水位期的径流槽消失，并形成高水头，迫使部分岩溶地下水向滑坡区(Ⅱ1区)内径流，相继夜郎组内的局部临时分水岭也消失。由于侧向补给量大，且碎屑岩(Ⅱ1区)基岩裂隙过水能力小，部分地下水在碎屑岩和灰岩接触带溢出地表而形成泉，这也是gw17和gw18等泉点主要成因。

表3 滑坡区水化学特征

表4 打邦河流域相关地层水化学特征

表5 月降水量统计表

图2 滑坡区剖面水动力场Fig.2 Hydrodynamic field of the landslide section

岩溶区(Ⅰ2区)形成地下水富集并向滑坡区径流补给的同时，滑坡区地层夜郎组、龙潭-长兴组碎屑岩(Ⅱ1区)也接收了强大的垂直入渗补给。滑坡区自身形成的裂隙水动力叠加永宁镇组1段岩溶区(Ⅰ2区)的水头压力作用，当滑坡区内水头达到某种高度时，使得潜在软弱结构面以上岩层失稳而形成本次大型滑坡。

假设P为岩溶区水头、p为滑坡区(A区)内自身形成的碎屑岩裂隙水水位、W为滑坡岩层重力、f为滑动面摩擦阻力，两种地下水共同作用水动力应力场σ(P+p)和岩层应力场σW在滑动面上与摩擦阻力f的平衡关系，即滑坡区水岩耦合模式关系：

σ(P+p)+σW=f

(1)

由于岩溶地下水往往具有快速流动特征，在几十分钟或数小时就会影响到下游数千米距离处水位[20-21]，因此，本案例中考虑上部岩溶地下水影响，不仅需要考虑滑坡区直接后缘岩溶区域，也要考虑后缘以东岩溶区域，以及间接影响岩溶区补给条件的碎屑岩区，即图1中ab线段以东至分水岭的永宁镇组1段B区和永宁镇组2段C区，对应面积分别为0.93，0.30 km2。

2.3 滑坡滞后暴雨分析

2010年6月26日至6月28日强降水过程合计降水296.5 mm。暴雨主要发生在28日6时—12时，合计降水191 mm，其中小时降水量最大达52 mm；12时以后为小雨，其中13时—14时2个小时降水量仅为5.9 mm(表6)。滑坡发生于14时，滞后暴雨2 h，即滑坡区内动、静水压力14时达到最大，因此滑坡暴雨滞后问题在水文地质上转化为滑坡区内最大水动力场产生问题。

滞后暴雨问题与滑坡区本身工程地质结构、降水过程、上游补给区大小、含水介质结构等条件密切相关，一些学者统计得出西南石山区滑坡在暴雨当日或滞后1 d产生的可能性最大[22-23]。本次滑坡区范围南北长280 m，东西宽200～250 m，后缘距离灰岩区仅25～50 m，同时又与两侧裂隙水没有联系，如果不考虑上覆岩溶水影响作用，仅考虑滑坡区碎屑岩自身基岩裂隙水，伴随13时—14时降水和补给量减少以及经过2 h的径流排泄，大部分滑坡区内的水位已经处于下降过程，14时的动、静水压力不是最大峰值。

永宁镇组1段岩溶区地下水，从东部分水岭至滑坡区后缘最大径流距离3 450 m，在暴雨2 h后于滑坡区后缘形成最高水头并传递到滑坡区内，符合实地条件和岩溶管道地下水运动规律。如本文作者实测出桂林寨底岩溶管道1 h影响距离2 230 m[20～21]。因此可以推论，本次滑坡是基岩裂隙水叠加岩溶管道水共同作用的结果，且岩溶管道水动力是主要激发和启动因素。

当然，大寨特大型滑坡存在如相关文献提出的存在高陡斜坡、基岩节理裂隙发育等工程地质方面的滑坡失稳基本要素，但独特的地下水作用是启动或触发本次滑坡以及形成泥石流最重要的因素。从另一角度也说明，泥石流中水体不仅仅来源于滑坡区本身碎屑岩内部，更多的来源于后缘岩溶区，这是泥石流沿沟槽运动更远、破坏力更大的原因。

表6 10月28日0—15时小时降水量

2.4 水动力场耦合

解决上述带有滞后特征的复合水动力场问题有多种方法，采用分块独立计算后再进行耦合是其中的有效方法之一。表示B区岩溶地下水(式2)和A区基岩裂隙水(式3)一维(也可采用二维)水流模型为：

(2)

(3)

式中：Tx——x方向导水系数/(m2·h-1)；

P——水头/m；

R——源汇项/(m·h-1)；

S——储水系数；

p——水位/m；

Ky——y方向渗透系数/(m·h-1)；

r——原汇项(1·h-1)；

μ——单位储水系数(1·m-1)；

t——时间/h。

根据经验，通过调整导水系数Tx很难反映暴雨延迟效应，可通过边界条件达到滞后耦合：

p(y,t)=P(x,t+Δt),x或y∈Γ,t、Δt≥0

(4)

式中：x、y——岩溶区与碎屑岩接触边界Γ上的点；

Δt——滞后时间/h，其值为2。

本文重点关注水动力滑坡影响模式或概念模型，关于水流模型解算的初始条件和其他边界条件以及碎屑岩C区对岩溶B区补给条件处理等不进行深入讨论。

3 结论

(1)大寨滑坡区后缘山体顶部存在地下水分水岭，分水岭两侧地下水补径排条件主要受岩性和含水介质控制。接受补给后，岩溶区地下水以管道集中径流方式由东向西朝北盘江排泄，碎屑岩区基岩裂隙水以分散方式向最近的冲沟或槽谷径流排泄。

(2)滑坡区位于分水岭北侧，滑坡区后缘永宁镇组1段0.93 km2的岩溶区地下水与滑坡区碎屑岩裂隙水存在水力联系，其他区域，特别是下陇古岩溶地下河gw20与滑坡区没有水力联系。

(3)大寨滑坡及泥石流的发生，是暴雨后滑坡区碎屑岩体内裂隙渗流叠加上覆永宁镇组1段岩溶管道快速流共同作用的结果，即复合地下水动力场滑坡影响模式；滞后暴雨2 h启动，是后缘岩溶区管道水对暴雨延迟效应的反应。