1、建设与管理水利技术监督2023 年第 2 期DOI:10.3969/j.issn.1008-1305.2023.02.025大兴庄水库降水作用下大坝稳定性分析王起芝(山东省临沂市费县水利局，山东 临沂 273400)摘要:库水位变化是当前引发大坝破坏的重要因素之一，频繁的库水位变化会降低坝体内部的抗剪强度，及时开展水库水位变化后大坝的渗流状态的研究十分重要。文章为分析费县大兴庄水库在降水作用下的稳定性，采用GEO-STUDIO 软件并考虑非饱和土特性，分析了坝基土的导水率、水库水位下降率、水库水位下降等因素对大坝上游边坡稳定性的作用。研究表明，库水位下降会降低坝体稳定性，采用高渗透性土筑成的大
2、坝在水库水位较低时更加稳定。关键词:降水;稳定性;渗透率;大兴庄水库中图分类号:TV698.23文献标识码:B文章编号:1008-1305(2023)02-0097-03收稿日期:2022-10-24作者简介:王起芝(1975 年)，女，工程师。E-mail:我国水库数量位居世界榜首，多数修建于 20世纪 5070 年代，但多数水库由于施工处理不当、填筑标准低等问题，长期处于“带病”运行的状态，时常出现坝体裂缝和坝坡滑移、溃坝等事件1-2。库水位变化是当前引发大坝破坏的重要因素之一，频繁的库水位变化会降低坝体内部的抗剪强度，及时开展水库水位变化后大坝渗流状态的研究十分重要3。岑威钧等4 利用刚
3、体极限平衡理论，通过考虑土工膜后渗流对坝坡稳定的作用，对其进行了数值模拟和分析;唐书童等5 根据水文监测数据，对城西水库堤防渗漏的安全性和今后高水位作业情况下的防渗措施进行了研究;贺芳丁等6 对水库的蓄水安全性评价与水库大坝的安全性评价进行了比较;尹芳芳等7 从多方面对磨刀坑水库大坝安全性进行了全面复核和评价，针对大坝的结构渗流稳定等进行了分析和论证;罗江锋8 参考一个除险加固的工程实例，考虑渗流场和应力场间的耦合，对土石坝加固前、后的坝体受力作用进行计算分析;孙玮玮等9 针对黄河源电站的实际情况，采用数值模拟方法，对黄河源发电站 2 种不同的渗流场进行了数值模拟和计算，对其渗透安全性进行了综
4、合评估;简鸿福等10 根据试验资料，对堤体各个岩石和土体的力学性能进行了反演，以确保堤防渗漏的安全性;岳庆河等11 建立老岚水库土石坝渗流典型断面有限元模型，对不同上下游水位工况及假设可能出现的不利工况进行数值模拟，。文章以大兴庄水库为例，采用 GEO-STUDIO 软件并结合非饱和土性质，主要探讨了大坝的坝基土的导水率、水库水位下降率、水库水位下降等因素对大坝的作用，研究成果可为相关工程提供参考。1工程概况大兴庄水库坝址在薛庄镇大兴庄村北 200m，位于沂河水系的薛庄河，该水库的控制流域为2.8km2，库容22.6 万 m3，兴利库容7.3 万 m3，死库容 0.08 万 m3。大兴庄水库坝
5、体是一种均质坝，坝体宽 4.5m，高 10m。上下游坡比均为 12.0。大兴庄水库地处季风型、半湿润地区，是一种具有显著季节性的温带大陆性气候。春天多风，降水稀少，夏天多雨，秋天有丰富的阳光，冬天又冷又干燥。本区年平均降水量为 852.5mm，降水年际差异大，造成了干旱、洪涝年的不均衡;年平均日照时间为 2529.4h，平均日晒百分比为 59%;最大冻结深度为 28cm;年平均蒸发面积达 1605mm。该区年均最高风速为3.4m/s，最大风速为17.0m/s。风向最多的是东风。2数值模型及计算参数本文建立的大坝数值模型结构如图 1 所示。大坝横断面的尺寸、材料区域和土壤性质根据792023 年
6、第 2 期水利技术监督建设与管理图 1大坝数值模型(单位:m)现场试验取得，在土体非饱和状态下，土体含水量特性曲线是土体结构研究中的一个重要环节。因为实测的 SWCC 不能获得大坝的具体用材，因此本文的心墙计算参数根据勘测数据进行估算取得，大坝和沙滤层是由其颗粒大小来估计的。本次采用GEO-STUDIO 软件进行分析，采用瞬态渗透法测定了水库在下游的孔隙水压，并对上游坝体进行了稳定计算。2 种方法均在自动耦合模式下进行。计算参数见表 1。2.1渗流方程通过 SEEP/W 软件进行了渗透数值计算，分析了大坝在水力作用下的二维稳定渗流场和水压作用下的瞬时渗透特性。应用有限单元法，对渗透问题进行了数
7、值模拟，得到了由公式(1)表示的广义的控制微分方程:xkxH()x+ykxH()y+Q=t(1)式中，H总 水 头，m;kxX 方 向 的 导 水 率，cm/s;kyY 方向的导水率，cm/s;Q边界通量，m3/s;体积含水量，无量纲;t时间，s。坝体有限单元网格的剖分情况如图 2 所示。图 2坝体的有限元网格剖分状况网格剖分分别由 4 和 3 个结点组成的四边形与三角形单元。假定坝顶的最高水位(19.2m)是 ABC与 DE 的界限的一个关键条件，即一个常量的总水头，它的数值为 19.2 和 9.2m。在此基础上，通过稳定的渗透计算，得到了初步的孔隙水压力分布。上游的边界线是以总水位的形式来
8、确定的，这样的界限使得使用者可以根据时间来确定总水位的改变。在瞬时渗透计算中，随着水库水位的降低，上游边坡面的水压由 19.2m 直线降低到 9.2m。在水头降低过程中，坝体由饱和向非饱和转变，因而，各物料的饱和导水系数被称为吸力的函数。通过对坝体在降水过程中的瞬态渗透进行数值模拟，得到了全网格内孔隙水压随着水位的降低而发生的动态变化。2.2坝体稳定分析方程为了准确分析出水位下降时，坝体上游边坡的稳定性。本文利用 SLOPE/W 中已有的 GLE 法对边坡失稳进行了数值仿真。通过计算各个斜坡的潜在失效面，得到了在不同的水位变化条件下，最大安全因子的最优解。在这种情况下，由于大坝的水位降低，出现
9、了黏土心墙坝、沙滤层等坝体不饱和的情况，故本文将材料的不饱和问题和剪切强度也考虑在内，并给出了相应的计算公式(2)。利用抗剪强度参数和 SWCC 来估算材料的非饱和抗剪强度:=c+(n ua)tan+(ua uw)w rs()rtan(2)式中，饱和或非饱和土的抗剪强度，kPa;有效摩擦角，();(ua uw)土壤破坏面内的基质吸力，kPa;Ua非饱和土壤中的孔隙空气压力，kPa;Uw孔隙水压力，kPa。表 1数值计算参数材料名称符号单位数值心墙坝体砂滤层基础饱和导水率ksatm/s1.2 10103 1074.9 1033.9 106饱和单位重量satkN/m32020.520.519.2有
10、效黏聚力ckPa718021有效内摩擦角2930.83432.3体积压缩系数Mvm2/kN18 10512 106012 106体积含水量wm3/m30.40.350.310.33液限LL%5130031塑限PL%201601389建设与管理水利技术监督2023 年第 2 期公式(2)可以用于瞬态渗透计算得到孔隙水压力变化，可应用于上游坝坡在各种水位下的稳定性分析。3数值结果分析3.1导水率对边坡稳定性的影响为分析土体饱和导水率(ksat)对上游边坡稳定性造成的影响，本文对不同水位下降情况下出现的4 种不同 ksat边坡安全系数变化情况进行分析，如图 3 所示。图 3不同 ksat下土坝上游边
11、坡的安全系数随水位变化由图 3 可知，对于给定 ksat下，随着大坝的水位降低，大坝的上游边坡面稳定性逐渐降低，在库容降低到 1/2 2/3 时为最低。这主要是因为大坝的蓄水量减少和大坝暂时不排水造成的，而且随着孔隙水压力的增大，土体剪切强度也会随之减弱，从而降低了安全系数。在蓄水过程中，在蓄水量由原来的 1/2 降低到 2/3 的情况下，稳定性由最低值逐渐增大，直到水位不再下降。坝体材料的抗剪强度因其孔隙水的消失而增大。这对斜坡体系的稳定起到了很好的作用。当水库的水位降低时，从图 3中可以看出，土体的饱和导水率(ksat)对上游边坡稳定性有显著影响。在水位下降的早期阶段，堤坝的上游边坡采用低
12、饱和导水率建造边坡比高导水率更稳定。然而，随着水位的进一步降低，边坡稳定性会随着 ksat值升高而改变。值得注意的是，饱和带(潜水表面以下的区域)在稳定状态下，高 ksat坝体的渗流条件大于低 ksat坝体的渗流条件。由于高ksat土壤允许水更深地渗透到大坝中，饱和区会逐渐扩大，这是坝体在水位下降早期阶段，安全系数降低的主要原因，低 ksat土壤中较小饱和区域的孔隙水压力可能比高 ksat土壤中较大饱和区域的孔隙水压力消散得更快，因此，超孔隙水压力较小的边坡比超孔隙水压力较大或饱和区域较大的边坡更稳定。上游坝坡稳定性在饱和与非饱和状态下其抗剪强度的变化情况如图 4 所示。图 4饱和与非饱和抗剪
13、强度下边坡安全系数变化趋势由图 4 可知，在水位下降期间，通过计算抗剪强度，得出的结果比非饱和抗剪强度的结果稍好。因为非饱和物料可以使水库的水流更深入地渗透到堤坝中，从而使堤坝更加湿润，这是导致非饱和物料堤坝安全性下降的重要因素。相对于不饱和物料坝的高湿度区域，在饱和物料坝体的小湿度区域，其孔隙水压力的消散速度要快，所以在低孔隙水压或在低湿区的水坝要比在高水压条件下的水坝更稳定。但是，应该指出，大坝的饱和状态和现实中的不一样，比如黏土心墙、路堤和沙滤等，在大坝的水位高度到达大坝顶部时，也不会完全饱和。3.2水位下降率对边坡稳定性的影响本节对水位下降率对上游边坡稳定性造成的影响进行分析，分别为
14、8、6、2 和 0.5m/d，如图 5所示。图 5水位下降率对上游边坡稳定性的影响由图5 可知，当水库的水位下降率减小时，水库的上坡面稳定性明显提高。由公式(2)可知，由于坝顶的高度下降速率减小，导致(下转第 189 页)99理论研究水利技术监督2023 年第 2 期能够进行有效调节，且在 0.5 5.0a 重现期，对污染物的削减都能够达到 50%。SWMM 模型能够有效模拟降雨与调蓄池工作情况、掌握水质污染以及污染物削减情况，为建设更好的城市生态系统提供了科学依据，但研究只是将 SWMM 模型利用到个别方面，未来还可以探索模型更多的使用价值。参考文献 1房亚军，于川淇，金鑫，等 基于 SWM
15、M-CCHE2D 单向耦合模型的山地海绵城市内涝管控效果研究J 武汉大学学报:工学版，2021，54(10):898-906 2刘碧云，任建民，张进丽，等 基于 SWMM 模型的兰州交通大学校园区城市雨洪模拟及 LID 效果评价J 水电能源科学，2021，39(7):9-12 3席广朋，王建龙，赵梦圆，等 城市雨水调蓄池水质控制效果及其影响因素分析 J 环境工程，2018，36(12):98-102 4 Tuomela C，Sillanpaa N，Koivusalo H Assessment of stormwaterpollutant loads and source area contri
16、butions with storm watermanagement model(SWMM)JJournalofEnvironmentalManagement，2019，233(MA 1):719-727 5焦春蛟，吕谋，张士官，等 不同模式雨水调蓄池与低影响开发组合对雨洪控制的效果 J 科学技术与工程，2019，19(34):336-342 6王生愿，余黎，陈小龙，等 基于在线监测和 SWMM 模型的合流制排水体制系统化提升方案评估研究J 中国环境监测，2021，37(1):103-111 7陈弘，孙书洪，赵鹏 基于“一河一策”的黑臭水体治理措施研究J 水利规划与设计，2019(4):14-18 8 Luan B，Yin，Xu P，et alEvaluating Green StormwaterInfrastructurestrategiesefficienciesinarapidlyurbanizingcatchment using SWMM-based TOPSISJJournal of CleanerProduction，2019，223(JUN 20):680-691 9王东
