0 引言

大量事实表明降雨将对坝坡稳定性产生重要影响^[1]，实际降雨入渗降低坝坡安全系数情况可分为2类：前者为坝体本身材料性质的改变如黏聚力、内摩擦角及基质吸力等； 后者为坝体综合环境因素的变化，如降雨导致浸润线抬高，孔隙水压力升高使潜在滑移面的抗滑能力降低等^[2].目前，国内外大都致力于降雨对边坡渗流稳定性影响研究，如陈浩、黄静等^[3]利用Geostudio 软件对降雨条件下的土质边坡内部渗流场进行模拟，分析了地下水位线、坡体饱和度、压力水头随降雨时间的变化规律；谢瑾荣等^[4]运用二维渗流数值计算方法，对降雨条件下的边坡孔隙水压力大小及暂态饱和区面积在空间及时间上的分布进行了模拟；李兆平等^[5]以土壤体积含水率作为控制变量，改变边界条件模拟了降雨过程中土壤入渗能力的变化，然而，针对尾矿坝在降雨条件下渗流稳定性研究较少，此外，类似上述研究未能建立不同降雨强度与尾矿坝可能失稳时间一一对应的关系，离解决工程实际问题仍有较大差距，要想深入了解降雨入渗条件下坝坡稳定性问题，首先要对其进行正确的模拟，本文以新余良山太平尾矿坝为工程对象，应用Geostudio软件中的SEEP/W 模块模拟不同降雨强度入渗工况下坝坡渗流场发展趋势，并结合SLOPE/W 模块在相应工况下对其稳定性进行数值计算分析，研究成果对太平尾矿坝在实际降雨条件下安全稳定运行具有重大参考和借鉴意义^[6-8].

1 太平尾矿坝区域概况

太平尾矿坝坝址位于新余良山铁矿.初期坝为透水石坝，是在水库大坝基础上改建而成的.原水库大坝坝基标高130.60 m, 坝顶标高151.30 m, 坝顶宽10.00 m, 内外坡比均为1∶1.1，为混凝土心墙堆石坝，改建后初期坝为透水堆石坝，坝顶高程160.00 m, 坝顶宽度4.00 m, 上游边坡坡率1∶1.5，下游边坡坡率1∶1.65，坝高29.40 m.尾矿库采用上游式尾矿堆坝法筑坝，目前尾矿堆积坝坝顶标高约206.5 m, 总坝高约75.90 m, 属Ⅲ级尾矿库和三等堆积坝.调查未见明显裂缝、变形位移等迹象.

据当地气象资料，良山镇雨量充沛，但时空分布不均，年实际变化较大，多年平均降雨量1 403 mm, 降雨量主要集中在4～6 月，占全年降雨量的47 %.常有阴雨天气.年降雨量变化最大值达510 mm, 降雨成为可能诱发坝坡失稳的突出因素.

2 建模及计划方案的确定 2.1 模型网格划分

根据良山太平尾矿坝相关工程地质条件，针对选取的典型剖面，利用SEEP/W 模块建立坝坡渗流场模型.具体有限元网格剖分见图 1，网格中共有3 058 个节点和3 020 个单元.

2.2 模型参数

模型参数选取主要包括岩土体的力学参数以及渗流相关参数的选取.相关物理力学参数取值见表 1.

渗流参数的正确选取直接关系到最终模拟效果的准确性.根据各层岩土体岩性不同确定体积含水量样本函数，结合土体饱和含水率计算出岩土体土水特征曲线见图 2，即基质吸力与体积含水量关系图，而渗透系数K 与基质吸力、体积含水量都有直接关系^[9-11]，因此，可采用Van Genuchten 法结合资料中每层土体的饱和渗透性系数，推导出的渗透性函数特征曲线见图 3.

2.3 边界条件

根据实际降雨情况，整个坝坡表面都应有同一强度降雨，具体设为510 mm/d、154 mm/d、120 mm/d、70 mm/d 4 个等级流量边界；模型底面定义为不透水边界^[12-14].

2.4 计算工况

为了综合分析降雨强度、降雨持续时间对滑坡稳定性的影响，通过查算江西省暴雨洪水查算手册，结合当地实际降雨情况可得出相应重现期的不同降雨强度值，以暴雨重现期为依据得到的降雨强度值分别为510 mm/d、154 mm/d、120 mm/d、70 mm/d 4 个等级见表 2.

重现期代表当地出现同等强度降雨的时间周期，降雨强度频率表示单位时间内某种事件出现的次数或比例，水文统计上，用频率反映水文时间出现的频繁情况.

为了更好地突出降雨历时对坝坡稳定影响程度，得出不同降雨对应的预计失稳时间关系图，降雨持续时间设置为60 d.每种降雨强度历经20 个时间步长，计算初始状态对应初始地下水位情形.

3 渗流场分析及稳定性计算 3.1 渗流场结果分析

通过对太平尾矿坝进行瞬态非饱和渗流模拟分析，可得到不同时刻不同降雨等级工况下孔隙水压力等值线图，限于篇幅，仅列出初始状态孔隙水压力等值线图及第②级降雨强度（50 年一遇）工况下分别持续3 d、9 d、30 d、60 d 时刻坝体内孔隙水压力等值线图见图 4~图 8.

由图 4~图 8 可知，降雨开始后，坝坡顶部和坝坡坡面以下一定区域范围内孔隙水压力升高，基质吸力相应降低；随着降雨持续进行，负孔隙水压力分布区进一步缩小，基质吸力逐渐降低，当上部局部饱和带与坡体内的饱和带汇合贯通后，浸润线开始整体向坡体内侧后移，15 d 以后，坝体基本处于饱和状态.同样，当降雨强度分别处于第①、③、④级工况时，渗流场分布情形与上述情况类似，只是在时间上分别相应提前9 d、推后3 d、推后18 d.

3.2 稳定性分析

采用非饱和土抗剪强度理论，将SEEP/W 的分析结果代入SLOPE/W 中进行耦合计算，用Morgenstern-Price 法分析降雨时坝体的瞬态稳定性见图 9~图 12，并进一步统计降雨强度、降雨持继时间与滑坡稳定性的关系见图 13.

结果表明：随着降雨持续天数的增加，坡内土体抗剪强度降低，坡体的稳定性逐渐下降；降雨强度越大，坝体安全系数下降越明显.深入分析各级降雨强度工况下降雨持时与安全系数之间关系，进一步通过拟合得到上述因素之间关系公式及对应相关系数如下：

①级：

(1)

相关系数R=0.992

②级：

(2)

相关系数R=0.99

③级：

(3)

相关系数R=0.992

④级：

(4)

相关系数R=0.988

式（1）~式（4）中：Y为安全系数；X 为降雨持续天数.

太平尾矿坝为三级尾矿坝，其对应的最小安全系数为1.10 见表 3，若太平尾矿坝处于第④级（10 年一遇） 工况时，60 d 内持续降雨坝坡安全系数均大于1.10，满足尾矿库安全规程要求，若处于第③级（25年一遇）工况时降雨工况时，第27 天前，尾矿坝尚处于基本稳定状态，接近第36 天时，坝体处于欠失稳状态，36 d 后，坝体整体失稳；若尾矿坝处于第①、②级降雨工况时，基本情况与前述情况类似，只是坝体达到失稳状态较之前分别提前31 d、16 d 左右.

进一步分析可知，当降雨强度值大于100 mm /d降雨持时60 d 范围内时，从理论分析，坝体将失稳，且通过计算可得出降雨强度值与可能发生失稳时对应降雨持时之间关系拟合曲线图见图 14，降雨强度范围在100 mm 到1 000 mm 之间.通过拟合得出二者关系如下所示：

(5)

相关系数R=0.979

式（5）中：Y 为尾矿库预计失稳时间；X 为降雨强度值.

4 结论

以良山太平尾矿坝为研究对象，借助Geostudio软件对不同降雨量和降雨持续时间下矿坝的稳定性进行研究，得出了以下结论：

1） Geostudio 的SLOPE/W 模块与SEEP/W 模块进行耦合计算，为降雨条件下良山太平尾矿坝的稳定性分析提供了基础；

2） 同一强度降雨工况下，随着降雨持续时间延长，尾矿库坝的稳定性逐渐降低，在坝脚处渗流集聚，易造成局部失稳； 同一时刻不同降雨强度工况下，随着降雨强度的增强，尾矿坝安全系数逐渐降低；

3） 经模拟和统计，获取了降雨量、降雨持续时间与坝坡稳定性的关系曲线，针对不同的降雨强度，均可参照关系曲线，对尾矿库的运行状态进行安全预测和警戒，工程意义深远；

4） 影响尾矿库稳定性的因素众多，文章只对降雨进行了讨论，只能为尾矿库的安全性提供参考，而不能作为决断性依据.