1. 引言
近十年来,全球自然灾害呈现出发生次数越加频繁化、发生规模出现极端化、发生可能难以预测化、发生后果衍生级联化等新特点,给防灾减灾构成严峻挑战 [1] 。水库在蓄水运营过程中,水流升降对坝体的渗流特性,会一定程度对坝体渗流与应力变形产生较大的影响,并促进坝体内部渗流活动和增大坝体应力集中,造成坝体失稳。在1928年圣佛兰西斯坝突然溃决约70 min内库水全部泻出,水位骤升坝基内孔隙水压增高,导致溃坝的发生,造成重大损失 [2] 。因而水利大坝类渗流安全稳定性是水库运营重要关注方面。
截止2009年底,全国已建成各类水库大坝87,151座 [3] ,大坝在防洪、供水、灌溉和发电等方面具有重要的作用,发挥着巨大的经济社会效益 [4] 。渗透稳定是影响土石坝安全的重要因素 [5] 。水库大量蓄水会造成局部岩体应力积累、地面沉降以及恢复过程也涉及渗流中的流固耦合问题。在流固耦合中,渗流流体计算和机械力学计算相互影响。王芝银 [6] ,郑会耀 [7] 通过建立了渗流场与应力场的动态全耦合流变模型,研究不同时间的位移、应力和渗透系数的分布。随着计算机技术的快速发展,大变形数值模拟方法在岩土工程领域有了广泛应用 [8] 。朱一飞 [9] 用ANSYS对阜新电厂四灰场主坝进行了数值模拟,对其位移、应力应变等结果,对其精密监测提供依据。黄岗水库位于东南沿海地区,沿海极值水位呈显著增长,更容易造成极端天气和气候所带来的地质灾害,而极端降水的发生又会导致人口较密集的东南沿海城市遭到更严重的洪灾影响。利用FLAC3D软件对黄岗水库大坝进行实时监测,及时避免灾害的发生。
近年来,这些问题都会造成巨大的损失,让我们意识到排除这些隐患对社会经济的重大影响,水库的实时监测可以提前预测事故的发生,以达到大大减少人员伤亡和经济损失的目的。
2. 工程概况
黄岗水库是一座以灌溉、供水为主,结合防洪、发电、养殖等综合利用的重要中型水库。
黄岗水库位于福建省龙岩市。受地理环境和人类活动等因素的综合作用的影响,我国东南沿海地区沿海极值水位呈现显著增长趋势,且更容易发生由于极端天气和气候事件等带来的各类地质灾害。
水库正常蓄水位457.00 m,库容2520万 m3。水库由主坝,副坝,溢洪道,放水系统等部分组成。其中主坝为均质土坝,坝顶高程462.00 m,最大高47米,坝顶长115米。副坝为高12米,长50米的土坝,大坝剖面图尺寸如图1所示,右侧建有侧槽式非常溢洪道,遇特大洪水时,最大可下泄384立方米/秒。
3. 模型建立
针对黄岗水库拟建土坝工程渗流场特征计算原则,水在非饱和土中的渗流也服从达西定律,采用水头h作为控制方程的因变量 [10] 利用FLAC3D有限差分软件,建立坝体模型如图2所示。该计算模型包括黄粘土坝身,排水棱体等,基于微单元网格划分,共获得82,584个单元体,67,330个节点,建立模型长为115 m,坝高为47 m,上游、下游的坡比均按照实际尺寸进行建模。计算模型中的x,z分别表示下游河流向,垂直向上方向 [11] 。边界荷载在坝体四周设置多向约束,顶部设置自由面,根据大坝中不同材料进行不同参数的设定。根据尺寸信息在犀牛进行建模,为了精确度,在FLAC3D中进行分组和命名。应用于本次数值模拟计算的所有参数如表1所示。
![](//html.hanspub.org/file/23-2751622x8_hanspub.png?20230705101655418)
Figure 2. Establishment of finite difference method model
图2. 有限差分法模型的建立
![](Images/Table_Tmp.jpg)
Table 1. Physical and mechanical parameters of materials
表1. 材料物理力学参数
运用FLAC3D软件模拟流体在孔隙介质中的流动及其相互的耦合作用,利用显示差分法研究其渗流过程。当流体在多孔介质中流动时,主要引起孔隙水压以及饱和状态的改变 [12] ,这些变化主要由流体质点平衡方程和Darcy定律方程式(1)所示之间相互关系来描述流体的流动。
(1)
其中:v为断面A上的平均流速;Q为流量;k为渗透系数;j为渗透坡降。
将犀牛建立的三维模型导入FLAC3D中,编辑fish语言,设定渗流参数,边界条件设置,进行渗流计算,得到边坡渗流稳定状态下的孔隙水压分布。随后关闭流体进程,将流体压缩模量设置为0,打开力学进程,利用强度折减法计算安全稳定系数。
通过在FLAC3D软件中建立定义方程式(2)来实现水库的实时监测。
(2)
其中:x为监测到的库水位的值减去坝底标高417.7;y为出水边界;k为库水位变化处的斜率。
4. 计算结果与分析
4.1. 孔隙水压的研究
根据渗流场特征计算,获得坝体静水工况下孔隙水压力特征,如图3所示,工况1 (正常蓄水位457 m处)其中孔隙水压力整体分布从坝底至坝顶递减,上游坝底处分布有最大孔隙水压力,达431.08 kPa,坝身各区间孔隙水压力分布呈块状递减式,上游临水侧孔隙水压力高于下游背水侧,且在下游背水侧孔隙水压力出现负值,表明坝身土体处于非饱和状态,最大负孔隙水压力为9.928kPa。
![](//html.hanspub.org/file/23-2751622x11_hanspub.png?20230705101655418)
Figure 3. Pore water pressure in working condition 1 (normal water level 457 m)
图3. 工况1(正常蓄水位457 m)的孔隙水压
![](//html.hanspub.org/file/23-2751622x12_hanspub.png?20230705101655418)
Figure 4. Pore water pressure in working condition 2 (water level rising to 460 m)
图4. 工况2 (水位上升至460 m处)的孔隙水压
1) 水位上升阶段
当水库蓄水位增高后,此时孔隙水压力分布基本与工况1一致,但量值上出现显著差异,工况2 (水位上升为460 m处)如图4所示,上游坝底最大孔隙水压力增大至445.08 kPa,从坝底至坝顶依然为递减分布,局部块状分布形态与工况1一致。根据浸润线可以看到浸润线有所升高但并未到达坝体高度,处于安全状态下。
2) 水位下降阶段
根据渗流场特征计算,获得坝体静水工况下孔隙水压力特征,工况3 (库水位骤降至455 m)当水库蓄水位降低后的孔隙水压,如图5所示,此时孔隙水压力分布基本与工况1一致,但量值上出现显著差异,从坝底至坝顶依然为递减分布。由于水位骤降,孔隙水压在坝体内逐步蔓延至下游背水侧,导致背水侧非饱和土体逐步吸收水分,非饱和土体分布量减小,进而负孔隙水压有所降低,工况3下游坝底最大孔隙水压力降低至7.892 kPa。
![](//html.hanspub.org/file/23-2751622x13_hanspub.png?20230705101655418)
Figure 5. Pore water pressure at working condition 3 (reservoir water level drops 455 m)
图5. 工况3 (库水位骤降455 m)处孔隙水压
4.2. 位移的变化
![](//html.hanspub.org/file/23-2751622x14_hanspub.png?20230705101655418)
Figure 6. Displacement in working condition 1
图6. 工况1位移
根据渗流场特征计算,获得坝体静水工况下位移变化特征,如图6所示,工况1中坝体位移从上到下总体呈逐渐递减的变化规律,最大的位移集中在坝体的顶部,最大位移达到0.024 m。
1) 水位上升阶段
当库水位骤然升高至工况2库水位为460 m时,此时位移分布基本与工况1一致,如图7所示,由于库水位上升,土中的附加应力减小,地基承载力大幅度下降,考虑其中水的浮力影响,上升承载力降低。与工况1位移相比,下滑最大位移有所增大,增大至0.146 m。
![](//html.hanspub.org/file/23-2751622x15_hanspub.png?20230705101655418)
Figure 7. Displacement in working condition 2
图7. 工况2位移
2) 水位下降阶段
水位降落尤其是骤降于工况3至455 m处的位移如图8所示,由于水位骤降,坡体水位高于库水位,由此产生的渗流作用使潜在滑动面的抗滑能力降低,滑动位移有所增大。
![](//html.hanspub.org/file/23-2751622x16_hanspub.png?20230705101655418)
Figure 8. Displacement at working condition 3 (455 m plunge of reservoir water level)
图8. 工况3 (库水位骤降455 m)处位移
4.3. 应力及安全系数研究
根据渗流场特征计算,获得坝体静水工况下垂直有效应力变化特征,如图9所示,整体呈现压应力。在水位为正常水位457 m时安全系数为1.957。根据GB50330-2013《建筑边坡工程设计规范》 [13] 中的规定,如表2所示。经分析大坝处于稳定状态。
![](Images/Table_Tmp.jpg)
Table 2. Classification of slope stability
表2. 边坡稳定性划分
1) 上升阶段
当库水位骤然升高至工况2库水位460 m时,此时竖直方向应力分布基本与工况1一致,如图10所示,但应力有所增大。其安全系数为1.996,根据技术规范中的规定大坝处于稳定状态。
![](//html.hanspub.org/file/23-2751622x17_hanspub.png?20230705101655418)
Figure 9. Stress and safety factor in vertical direction of working condition 1
图9. 工况1垂直方向的应力及安全系数
![](//html.hanspub.org/file/23-2751622x18_hanspub.png?20230705101655418)
Figure 10. Stress and safety factor in vertical direction of working condition 2
图10. 工况2垂直方向的应力及安全系数
2) 下降阶段
根据渗流场特征计算,获得工况3坝体静水工况下垂直有效应力变化特征,如图11所示。相对于工况1压应力有所增大。降后安全系数为1.996,根据技术规范中的规定大坝仍处于稳定状态。
![](//html.hanspub.org/file/23-2751622x19_hanspub.png?20230705101655418)
Figure 11. Vertical displacement and safety factor at 3 working conditions
图11. 工况3处垂直位移及安全系数
4.4. 极端天气对大坝稳定性研究
在全球变暖的大环境下,极端天气发生的频率也逐渐增多,对极端天气与水文事件精准预测预报,需要相关领域的科研工作者进行原理性创新,精细化解析极端天气与水文事件的成因与机理,提升预测预报技术的精准性 [14] 。而大坝在时间的推移下存在一些损伤,遭遇极端天气的影响,很可能发生溃坝等重大的灾害,对大坝的稳定性产生极大的威胁。其中洪水,暴雨的增加也对大坝的防洪能力产生巨大的压力,大大降低了水库大坝的安全等级。
黄冈水库溢洪道堰型为宽顶堰,溢洪道泄流能力公式如下:
(3)
式中Q-流量;B-总净宽;黄冈水库溢洪道总净宽10.0 m;G-取9.8;M-二元水流宽顶堰流量系数,与相对上游堰高P1/H及堰头形式有关;E-闸墩侧收缩系数;
(4)
注:n为闸孔数目。
时仍取1。
根据不同的泄洪量,由上述公式得到H0。
对于此情况,本文将采用数值模拟技术对100年一遇、1000年一遇的洪水进行模拟并进一步分析研究。根据上述式(3),式(4)综合计算所得100年一遇的洪水下泄洪水量为221 m3/s,库水位为457.75 m。1000年一遇的大洪水下泄洪水量为306 m3/s,库水位为459.25 m。
由下图12,图13,图14所示100年一遇的洪水浸润线,下滑位移以及安全系数显示,在安全范围内。
![](//html.hanspub.org/file/23-2751622x23_hanspub.png?20230705101655418)
Figure 12. Infiltration line of a 100-year flood
图12. 100年一遇洪水的浸润线
![](//html.hanspub.org/file/23-2751622x25_hanspub.png?20230705101655418)
Figure 14. Safety factor and vertical pressure of 100 year flood
图14. 100年一遇洪水的安全系数及竖向压力
由下图15,图16,图17所示1000年一遇的洪水浸润线,下滑位移以及安全系数显示,在安全范围内。
![](//html.hanspub.org/file/23-2751622x26_hanspub.png?20230705101655418)
Figure 15. Infiltration line of 1000 year flood
图15. 1000年一遇洪水的浸润线
![](//html.hanspub.org/file/23-2751622x27_hanspub.png?20230705101655418)
Figure 16. Displacement of a 1000 year flood
图16. 1000年一遇洪水的位移
![](//html.hanspub.org/file/23-2751622x28_hanspub.png?20230705101655418)
Figure 17. Safety factor and vertical pressure of 1000-year flood
图17. 1000年一遇洪水的安全系数及竖向压力
5. 结论
1) 稳定渗流工况下蓄水位上升,不改变坝身孔隙水压的分布特征,孔隙水压力均为从坝底至坝顶递减。最大孔隙水压由431.08 kPa增大至445.08 kPa。浸润线又所升高。同时位移分布基本与工况1一致,由于库水位上升,土中的附加应力减小,地基承载力大幅度下降,考虑其中水的浮力影响,上升承载力降低,与工况1位移相比,下滑最大位移有所增大,增大至0.14611 m。安全系数为1.965,经分析均在安全范围内。
2) 稳定渗流工况下蓄水位骤降,不改变坝身孔隙水压的分布特征,孔隙水压力均为从坝底至坝顶递减。下游非饱和的负孔隙水压9.928 kPa减小至7.892 kPa。浸润线有所降低。同时位移分布基本与工况1一致,由于滑坡渗流的影响,滑坡稳定性下降,与工况1位移相比,下滑最大位移有所增大,增大至0.024 m。安全系数为1.996,经分析均在安全范围内。
3) 经过对极端天气100年一遇洪水和1000年一遇洪水的模拟,检验了大坝在极端天气下的稳定性能。验证了该大坝满足渗流稳定要求。