Influences of antecedent rainfall and change of reservoir water level on seepage and stability characteristics of earth rock dams
-
摘要: 气候变化给水库大坝安全运行管理带来前所未有的挑战，前期降雨、库水位变化是影响坝坡失稳的关键因素，且碾压式土石坝设计规范未考虑降雨对大坝坝坡或浅层坝坡稳定性的影响，水库大坝滑坡或溃坝严重威胁人民的生命安全并造成巨大经济损失。以某水库心墙坝为研究对象，建立了有限元计算模型，研究了前期降雨、库水位变化对心墙坝渗流和稳定特性的影响，揭示了孔隙水压力、渗透比降和坝坡抗滑稳定安全系数的变化规律，并结合高密度电法和安全监测资料验证了计算模型的准确性。研究结果表明：前期降雨对大坝坝坡浅层范围孔隙水压力和抗滑稳定安全系数影响较大，对下游坝坡孔隙水压力影响程度由大到小依次为坡顶 > 坡肩（坡脚），随着降雨量增加下游坝坡10 m范围内孔隙水压力逐渐增大，坝坡表层土体达到饱和，主降雨后下游坝坡抗滑稳定安全系数减小较大；计算模型的浸润线高程与高密度电法、测压管水位高程均吻合，验证了计算模型的准确性；库水位是影响上游坝坡孔隙水压力和抗滑稳定安全系数的主要因素，对上游坝坡孔隙水压力影响程度由大到小依次为坡脚 > 坡肩 > 坡顶，降雨是导致下游坝坡孔隙水压力和抗滑稳定安全系数变化的主要因素；综合计算模型、高密度电法和安全监测资料，较好地分析了心墙坝孔隙水压力、渗透比降和抗滑稳定安全系数的变化规律，分析前期降雨、库水位变化为完善土石坝设计规范和水库大坝安全评价提供了科学依据。Abstract: The climate change has brought unprecedented challenges to the safe operation and management of reservoir dams. The antecedent rainfall and change of reservoir water level are the key factors affecting the instability of reservoir dam slopes, and the impact of rainfall on the stability of dam slopes or shallow dam slopes is not considered in the existing design code for rolled earth-rock fill dams. The landslide or dam break of reservoir dams seriously threatens human lives and causes huge economic losses. Taking the core wall dam of a reservoir as the research object, a finite element model is established, the influences of the antecedent rainfall and change of reservoir water level on the seepage and stability characteristics of the core wall dam are investigated, the variation rules of pore water pressure, hydraulic gradient and safety factor of anti-sliding stability of the dam slope are released, and the accuracy of the model is verified by using the high-density electrical method and the safety monitoring data. The research results show that the antecedent rainfall has a great impact on the pore water pressure and safety factor of anti-sliding stability in the shallow depth of the dam slope. The influence degree of the antecedent rainfall on the pore water pressure of the downstream dam slope has the descending order: slope top, slope shoulder and slope toe. With the increase of the antecedent rainfall, the pore water pressure within 10 m of the downstream dam slope gradually increases, and the surface soil of the downstream dam slope reaches the saturated state. After the main rainfall, the safety factor of anti-sliding stability of the downstream dam slope greatly decreases. The phreatic line elevation of the proposed model is consistent with the measured water level using the high-density electrical method and the piezometric tube, which verifies its accuracy. The reservoir water level is the main factor affecting the pore water pressure and safety factor of anti-sliding stability of the upstream slope. The influence degree on the pore water pressure of the upstream slope of the dam is in the descending order of slope toe, slope shoulder and slope top. The rainfall is the main factor inducing the change of pore water pressure and affecting the safety factor and anti-sliding stability of the downstream slope. The comprehensive model, high-density electrical method and safety monitoring data can be used to better analyze the change laws of pore water pressure, seepage gradient and safety factor of core wall dam. The analyses of the effects of the antecedent rainfall and change of reservoir water level provide a scientific basis for improving the design specifications of earth-rock dams and the safety evaluation of reservoir dams.
-
0. 引言
土拱效应是岩土工程领域重要的现象，显著影响地层应力传递和变形传播规律，常见于隧道、桩承式路堤、挡土墙等工程[1-5]。基坑工程中，土体开挖会导致坑底地层应力释放并产生不均匀回弹变形，引起土拱效应[6]。
Meng等[7]基于现场实测发现，当基坑开挖深度达到一定值时下方既有隧道上浮量会急剧增大。该现象与坑底土拱效应有关[6]，隧道残余覆土厚度随上方基坑开挖的进行逐渐减小，当土拱松动区发展至隧道范围时，隧道顶部地层竖向应力快速下降且土体松动弱化。此外，为控制既有地下结构变形，上方基坑工程通常采用坑底地层加固、分区开挖等措施[8-9]。这些措施实质上限制了坑底土拱松动区扩展，进而达到了减小地层不均匀变形、降低应力释放的作用。然而，上述措施设计参数依赖经验，包括加固范围、分区开挖尺寸等，因而难以实现变形控制要求[10-11]。以隧道为例，据统计，超过64%的基坑上跨运营隧道工程出现了隧道上浮超限的现象[12]。因此，有必要关注坑底地层的应力转移和变形传播机制，明确坑底土拱效应边界条件、影响范围和演化规律。
实际基坑工程中，因机械土方开挖作业，坑底地层应力和变形测试极为困难，坑底土拱效应的相关研究鲜见报道。离心模型试验可构建均质地层并使用高精度传感器，获取规律性结论，是研究此类问题的重要手段[13-14]。Ng等[15]、Shi等[16]和陈仁朋等[17]在基坑离心模型试验中测试了坑底回弹，但隧道的存在会显著影响结果。刘炀镔等[18]基于离心试验讨论了围护墙插入深度对基坑稳定性的影响，但并未关注坑底土体的应力场和变形场。综上所述，坑底土拱效应的研究尚缺少必要的试验支撑。
本文开展了干砂地层基坑开挖坑底土拱效应离心模型试验和不同开挖深度和宽度的多工况数值模拟。试验和数值分析重点关注坑底竖向应力和回弹变形。本文研究确定了坑底土拱效应边界迹线，揭示了基坑开挖坑底土拱效应形成-演化机制，明确了不同参数的影响规律。研究成果可为预测基坑坑底地层响应和下方既有地下结构变形控制提供支撑。
1. 基坑开挖坑底逆土拱效应
坑底土拱效应如图 1（a）所示，和活动门上方土拱效应（图 1（b））的主要区别在于土拱发展方向和力-位移边界条件。根据发展方向，可将活动门上方自下而上发展的土拱效应定义为“正土拱效应”。坑底土拱效应由上方卸载所致，其发展方向自上而下，与正土拱效应相反，可被定义为“逆土拱效应”。相应的，对于围护墙外侧土体，若存在因墙体变形水平卸载而产生的土拱效应，可被定义为“侧土拱效应”。
![]()
图 1 基坑坑底逆土拱效应和活动门上方正土拱效应示意图Figure 1. Schematic of downward soil arching below excavation base and upward soil arching above trapdoor随活动门下移，正土拱向地表发展演化，其土拱区和松动区边界上均存在摩擦力Ff（土-土界面）作用，阻碍土体向下滑动。与之相对的，坑底逆土拱会随开挖深度加深而逐渐向下发展，且基坑围护结构的存在会显著影响地层中的应力传递机制。逆土拱的侧部（围护墙边）和下部边界（围护墙底）上分别存在摩擦力Ff（土-结构界面）和支承力Fs，二者共同作用，限制了地层应力释放和回弹变形。本文将这种自上而下发展且围护结构影响应力传递的土拱效应定义为坑底逆土拱效应，并将显著受到开挖扰动影响的坑底浅部地层定义为松动区。
2. 离心试验方案
2.1 试验方案
本文开展的干砂地层基坑开挖坑底土拱效应离心试验在浙江大学ZJU-400土工离心机上完成。离心机有效旋转半径为4.5 m，试验过程中离心加速度保持为60g。本模型为左右对称的平面应变模型，图 2给出模型平面布置的示意图。试验用刚性模型箱内部尺寸为1000 mm×400 mm×1000 mm（长×宽×高）。根据相似定律，开挖基坑的模型和原型尺寸如表 1所示，原型基坑深度和宽度分别为24 m和15 m。若无特殊说明，后文中所给数据均为原型尺寸。
表 1 基坑模型与原型尺寸Table 1. Model and prototype dimensions of excavation变量名 模型尺寸/mm 原型尺寸/m 基坑开挖深度 400 24 基坑开挖宽度 250 15 围护墙嵌入深度 200 12 围护墙厚度 25 1.5 模型深度 1000 60 模型宽度 1000 60 2.2 试验材料
福建标准砂被广泛用于土工离心试验[16, 19]，其颗粒组成主要为细粒（图 3）。本试验通过“砂雨法”制备地基土，控制相对密实度为75%。基坑内土体开挖采用“排液法”[15-17]模拟，待开挖土体由ZnCl2重液替代，其实际密度与地基土差别约为2%。试验准备过程中先在开挖范围内贴设厚度为0.5 mm的柔性薄膜，开挖范围底部薄膜设有开孔并通过管路与模型箱外重液收集容器连接，然后在开挖范围内灌入重液。
本试验主要关注基坑坑底地层响应，为确保“排液法”顺利开展并限制围护墙变形，在围护墙外加支撑固定。围护墙和支撑模型均由铝合金材料制成，其杨氏模量约为72 GPa。试验过程围护墙水平位移可以忽略不计。
试验使用LVDT和微型土压计测试坑底回弹和土压力，具体布置如图 2所示。两种传感器精度分别为5×10-3 mm和1 kPa，可以保证测试结果的准确性。此外，其他断面布置了3对弯曲元，具体位置未在图中画出。弯曲元用于测试土体剪切波速，进而获得土体小应变剪切模量G0，结果在表 2给出。
表 2 地基土HSS模型参数Table 2. HSS model parameters for soil参数 数值 参数 数值 Erefoed/MPa 22.8 G0/MPa 86.4 Eref50/MPa 28.6 γ0.7 2×10-4 Erefur/MPa 103.0 einit 0.693 γ/kPa 15.53 K0 0.45 c'/kPa 0 νur 0.2 φ'/(°) 34.4 Rf 0.9 ψ/(°) 4.4 m 0.5 2.3 试验过程及结果
模型制备完成后吊入离心机吊篮。模型箱与重液收集箱通过管道连接，管道中设有电磁控制阀。开启离心机使模型重力加速度达到60g。待传感器示数稳定后，打开电磁阀开关开始排液，排液过程中保持重液流速恒定。试验过程中实时采集所有传感器读数。上机实物图像如图 4所示。
3. 数值分析方案
3.1 基于离心试验的数值模型
为进一步研究坑底逆土拱效应形成-演化机制，采用PLAXIS 2D开展基坑开挖坑底逆土拱效应数值分析。有限元模型基于离心试验原型尺寸建立，地层水力条件为干，模型边界设置为100 m（长）×100 m（宽），远大于2倍开挖范围，边界效应可以忽略。
基坑开挖岩土体的应变量级主要为小应变范围（10-6~10-3）[20]。为考虑土体小应变阶段模量的非线性、应力相关的特点[21]，地基土使用小应变硬化土（HSS）模型。土体本构参数主要由室内单元体试验获得[22]，其中剪切模量取试验时弯曲元测试结果的平均值，具体如表 2所示。因试验中围护墙变形很小，计算时限制围护墙位移为0。
3.2 多工况基坑开挖数值模型
基坑开挖导致的坑底逆土拱效应示意图如图 5（a）所示，Hl和Ha为基坑坑底到松动区和土拱区最大影响范围的竖向距离。基坑开挖深度He和宽度B是坑底逆土拱效应形成-演化的关键影响因素，显著影响拱效应的范围与发挥程度。为研究不同开挖尺度下坑底逆土拱效应的演化规律，共开展了25组模拟计算，各算例地层参数均与前文一致。计算工况如表 3所示，其中围护墙高度Hw均为36 m。
![]()
图 5 坑底逆土拱分区及其演化Figure 5. Division of downward soil arching below excavation base and its evolution表 3 基坑开挖计算工况Table 3. Simulation conditions of excavation变量 数值/m 开挖宽度B 5，15，30，45，60 开挖深度He 6，12，18，21，24 围护墙高度Hw 36 4. 离心试验及数值分析结果
4.1 离心试验结果及数值模型验证
试验前后地层竖向应力的比值是评价应力变化、反映土拱效应发挥程度的重要指标[23-24]。针对坑底逆土拱效应，定义竖向应力比为基坑开挖后某点竖向土压力与开挖前之比。坑内中心线上和围护墙底部的竖向应力比实测结果在图 6，7用实心点表示。图 6所示，坑内竖向应力比小于1，表明开挖后竖向土压力下降，坑底发生应力释放。然而，围护墙底部区域测试结果显示（图 7），开挖后土压力出现明显增大（竖向应力比＞1），这表明坑底存在应力转移现象，证实了坑底逆土拱效应存在。坑底回弹如图 8所示，呈“倒锅底形”，即拱形分布，最大回弹位于基坑中心。
![]()
图 6 基坑中心线上竖向应力比和归一化回弹Figure 6. Vertical stress ratios and normalized ground heave ratios in centerline of excavation为验证模型参数的准确性，将试验实测数据和数值模型的计算结果进行对比，结果显示开挖后坑底地层的响应基本吻合。基坑中心线上的竖向应力比与归一化回弹δ/δmax如图 6所示，竖向应力比的模拟结果略小于实测值。随深度增加，竖向应力比和δ/δmax逐渐增大和减小，最终趋于开挖前的状态。
在土拱效应的相关研究中，应力的拐点标志着地层应力状态的改变，常常被用于确定土拱边界[25-26]。但对于坑底逆土拱效应，竖向应力比随深度持续增加，曲线在数学上不存在拐点。因此，本文基于卸荷单元体试验中临界卸荷比和极限卸荷比以确定土拱区和松动区的高度，即Ha和Hl。既有试验结论显示，干砂地层中极限卸荷比约为0.8[27-28]，当卸荷比小于该值时，扰动效应明显且土体回弹急剧增加。该值对应的
地层竖向应力比为0.2，由此可确定松动区的边界。砂土地层中针对临界卸荷比的研究鲜有报道，因此考虑通过其他方法确定土拱区高度。分层总和法中，砂土地层的计算深度可根据地层应力变化小于20%确定。该深度以下不存在显著的地层相对位移且土体变形忽略不计，可视为坑底逆土拱效应的下边界，对应的地层竖向应力比为0.8。根据上述方法，可确定本模型中松动区和土拱区的最大高度分别为4.6 m和21.6 m。土拱区边界对应的δ/δmax约为20%。
围护墙底部竖向应力比如图 7所示。开挖后，试验和模拟均观测到该处竖向应力比大于1，实测最大竖向应力比为1.2。值得注意的是，PLAXIS 2D中板单元没有厚度，因此仅对围护墙外侧区域进行对比。距基坑较远处，实测竖向应力比减小至小于1，表明该区域出现竖向卸载，而有限元模型不能模拟这一现象。一般认为应力变化小于5%的区域内土拱效应的影响可以忽略不计[25]。此处将竖向应力比大于1.05的部分视为拱脚区，可确定该模型拱脚宽度Ba约为4.7 m，此时拱脚的边界也约为实测竖向应力比连线与竖向应力比为1.0直线的交点。
基坑坑底回弹如图 8所示，模拟结果与实测结果较吻合，均呈“倒锅底”形。坑底地层存在明显差异变形，最大和最小回弹量分别位于基坑中心和靠近围护墙处，差值为31 mm。数值模型中最大回弹量为44.8 mm，略小于实测结果。根据《建筑地基基础设计规范：GB50007—2011》推荐的坑底回弹计算方法[29]，计算的坑底最大回弹量为31.5 mm。该值小于实测值和模拟值，结果偏于危险，这可能与该方法未考虑卸载导致的地基土模量衰减有关。
4.2 坑底土拱分区
图 9给出基坑数值模型中基坑底土体偏应变云图，基坑坑底上方的无关部分已经略去。由前文分析可知，围护墙是坑底逆土拱效应的边界条件之一。基坑开挖后，土-结构界面产生切应力并向坑内传递，导致围护墙侧壁附近出现了显著的偏应变集中带。此外，围护墙底部以下也存在剪切带，其与水平方向的夹角约为61.5°，与主动朗肯状态下滑裂面和水平面的夹角相近（45°+φ/2=61.5°）。
图 10给出基坑数值模型中基坑底土体相对密实度云图。地基土初始相对密实度为75%，开挖后坑底绝大部分区域相对密实度下降且越靠近坑底土体越松散。其中松动区内相对密实度变化最为明显，土拱区其余范围内略有减小，而稳定区基本不变。基于图 10可知，松动区土体显著受到上方卸载和扰动的影响，表现为结构性减弱、力学性能劣化。此时土拱效应难以充分发挥，可能导致下方既有地下结构承受过大附加荷载进而出现变形超限。因此实际工程中有必要重点关注松动区影响范围，控制地下结构间的安全距离。值得注意的是，在围护墙底部，开挖后相对密实度反而增加，该现象由坑底逆土拱效应发挥和地层应力转移所致，与前文中该区域的实测和模拟竖向应力比结论相印证。
![]()
图 10 地层相对密实度分布Figure 10. Distribution of relative density of soil below excavation base图 11给出基坑数值模型中基坑底土体在不同深度处的归一化回弹结果。各剖面的间隔深度为4 m，同一深度处的最大回弹量在右侧给出。随深度增加，开挖后的地层响应减弱，最大回弹量减少，回弹曲线也更加平缓。因围护墙的阻隔作用，回弹主要发生在围护墙内部。值得注意的是，在围护墙内侧及围护墙底部地层（即图 9中偏应变集中带）回弹分布发生明显变化。在活动门上方正土拱效应的研究中，有学者将这种地层应变发生急剧变化的区域视为土拱的边界[31]。
根据前文所述方法，可同时取多个纵向剖面确定松动区和土拱区下边界迹线，如图 12所示。各纵向剖面距离基坑中心的水平距离xp为0，3，6 m，对应确定的松动区高度为4.6，4.1，3.2 m，土拱区高度为21.6，20.0，17.1 m。松动区和土拱区的下边界迹线函数均接近于开口向上的抛物线，其中松动区边界变化更加平缓。可以推测，松动区和土拱区下边界的形状参数与开挖深度He和宽度B密切相关。
4.3 坑底逆土拱效应演化规律
坑底逆土拱效应的演化规律类似于活动门上方正土拱效应。基坑开挖的深度He和宽度B类比于活动门试验中门板的下移量和宽度，二者的大小决定了卸载程度，与土拱效应的范围成正相关[30]。此外，区别于正土拱效应，围护墙深度Hw也会影响土拱的演化机制。坑底逆土拱效应随开挖深度He演化如图 5（b）所示。当开挖深度较浅时，松动区和土拱区范围较小，此时逆土拱效应仅存在于围护墙深度以内。随开挖进一步增加，逆土拱效应向下发展并超过围护墙底部，此时墙底竖向应力增加并产生剪切带。随后，逆土拱效应进一步发展，松动区和拱脚范围随开挖深度增大而逐渐增加。
图 13给出松动区高度随开挖深度变化关系。不同开挖宽度下，松动区高度与开挖深度基本呈线性正相关。此外，开挖宽度对松动区高度也略有影响。减小基坑开挖的深度和宽度有利于控制松动区范围。对各案例结果线性拟合可以得到松动区高度与开挖深度的关系（图 13），根据拟合关系可简单估算坑底的松动区高度约为0.2倍开挖深度。松动区内土体松散、变形显著，实际工程中应尽可能控制该范围与下方既有地下结构的安全距离。需要注意的是，该结论仅适用于干砂地层，不同地层条件下的规律仍有待进一步验证。
![]()
图 13 坑底松动区高度随开挖深度演化规律Figure 13. Evolution of height of loosened zone with increasing excavation depth图 14给出坑底逆土拱高度随开挖深度的关系，图中横纵坐标均基于围护墙高度归一化处理。随开挖深度增加，坑底逆土拱效应产生并向下发展，且其影响范围逐渐增加。不同开挖宽度的基坑坑底均呈现相似的规律。当开挖深度增加至接近围护墙底部时，明显可见土拱区高度的增长趋势变缓。此外，基坑开挖宽度也显著影响了坑底逆土拱效应的发展。相同开挖深度下，逆土拱高度随开挖宽度的增加而增加。值得注意的是，开挖宽度大于30 m时，这种变化不再明显。
![]()
图 14 坑底逆土拱高度随开挖深度演化规律Figure 14. Evolution of height of arching zone with increasing excavation depth由图 12可知，土拱区的下边界迹线函数接近于抛物线。为描述坑底逆土拱效应下边界的形状特征，定义下边界高跨比α=[Ha-(Hw-He)]/(B+2Ba)。如图 15，下边界高跨比α与He/Hw基本呈线性正相关。由此可知，坑底逆土拱效应随开挖深度的增加主要向下部发展而非向两边扩张。当He/Hw为0.374时，拟合直线与横坐标相交，该点表示逆土拱效应发展到围护墙以下的临界状态。为减少基坑开挖对下方地下结构的影响，实际工程中可适当增加围护墙插入深度以减小土拱下边界高跨比，避免土拱区发展至既有地下结构范围。
![]()
图 15 坑底逆土拱下边界高跨比Figure 15. Depth-span ratios of lower boundary of downward soil arching and its evolution with increasing excavation depth5. 结论
本文开展了干砂地层基坑开挖土拱效应离心模型试验及多工况下坑底逆土拱效应数值分析，主要得到以下3点结论。
（1）基坑开挖导致坑底产生土拱效应，地层中存在应力转移和差异变形现象。坑内土压力减小，而围护墙底部土压力增大，最大竖向应力比为1.2。坑底回弹呈拱形分布，最大和最小回弹分别位于基坑中心和靠近围护墙处，差值为31 mm。
（2）坑底逆土拱效应与活动门上方正土拱效应的主要区别在于土拱发展方向和力-位移边界条件。基于土拱效应可将坑底划分为松动区、土拱区和稳定区，松动区内土体相对密实度明显下降。
（3）基坑深度和宽度是影响坑底逆土拱效应形成-演化的关键因素。坑底松动区和土拱区范围与基坑宽度和深度呈正相关，实际工程中应尽量避免松动区发展到坑底下方既有地下结构。减小开挖深度和宽度、增加围护墙插入深度均有利于控制基坑开挖对下方既有地下结构的影响。
-
![]()
图 4 前期降雨作用下心墙坝孔隙水压力随时间变化图
Figure 4. Variation of pore water pressure of core wall dam with time under action of antecedent rainfall
![]()
图 5 前期降雨作用下上下游坝坡孔隙水压力随时间变化图
Figure 5. Variation of pore water pressure of upstream and downstream dam slopes with time under action of antecedent rainfall
![]()
图 6 库水位变化作用下上游坝坡孔隙水压力随时间变化图
Figure 6. Variation of pore water pressure of upstream dam slope with time under change of reservoir water level
![]()
图 7 库水位变化作用下下游坝坡孔隙水压力随时间变化图
Figure 7. Variation of pore water pressure of downstream dam slope with time under change of reservoir water level
![]()
图 8 前期降雨和库水位变化共同作用下上游坝坡孔隙水压力
随时间变化图
Figure 8. Variation of pore water pressure of upstream dam slope with time under action of antecedent rainfall and change of reservoir water level
![]()
图 9 前期降雨和库水位变化共同作用下下游坝坡孔隙水压力
随时间变化图
Figure 9. Variation of pore water pressure of downstream dam slope with time under action of antecedent rainfall and change of reservoir water level
![]()
图 11 #1测线高密度电法视电阻率反演剖面图
Figure 11. Apparent resistivity inversion profile of high-density electrical method by detection line 1
![]()
图 12 #2测线高密度电法视电阻率反演剖面云图
Figure 12. Apparent resistivity inversion profile of high-density electrical method by detection line 2
![]()
图 13 工况一水头等值线分布图
Figure 13. Contour distribution of water head under working condition 1
![]()
图 14 前期降雨作用下心墙坝抗滑稳定安全系数随时间变化图
Figure 14. Variation of safety factor of anti-sliding stability of core wall dam with time under action of antecedent rainfall
![]()
图 15 库水位变化作用下心墙坝抗滑稳定安全系数随时间
变化图
Figure 15. Variation of safety factor of anti-sliding stability of core wall dam with time under action of change of reservoir water level
![]()
图 16 前期降雨和库水位变化共同作用下心墙坝抗滑稳定安全
系数随时间变化图
Figure 16. Variation of safety factor of anti-sliding stability of core wall dam with time under joint action of antecedent rainfall and change of reservoir water level
表 1 坝体分区材料物理力学参数表
Table 1 Physical and mechanical parameters of dam partition
坝体分区 渗透系数/(cm·s-1) 饱和含水率/% 重度/(kN·m-3) 有效黏聚力/kPa 有效内摩擦角/(°) 坝前抛石 1.00×10-3 52.0 21.0 0 35.0 坝体代料 1.43×10-4 40.0 18.5 20.0 16.8 心墙 9.82×10-5 40.0 18.4 23.0 16.5 下游反滤 1.00×10-3 52.0 19.0 0 35.0 坝基 9.62×10-5 36.0 28.7 100.0 45.0 防渗墙 1.00×10-7 2.0 24.0 500.0 35.0 防渗帷幕 1.00×10-7 2.0 24.0 500.0 35.0 
下载: 导出CSV
表 2 计算工况表
Table 2 Working conditions
计算工况 库水位/m 骤变速率/(m·d-1) 降雨类型 工况一 100.00 — 最不利降雨 工况二 100.00~78.00 1 — 工况三 2 — 工况四 4 — 工况五 100.00~78.00 1 最不利降雨 工况六 2 最不利降雨 工况七 4 最不利降雨 注：2020年7月降雨开始前的初始库水位为100.00 m，死水位为78.00 m。
下载: 导出CSV
表 3 测压管实测值与计算模型计算值对比表
Table 3 Comparison between measured values by pressure tube and calculated values
(m) 日期 U1 U2 实测值 模拟值 差值 实测值 模拟值 差值 07-01 95.29 96.20 0.91 79.50 80.31 0.81 07-20 97.89 98.55 0.66 79.45 80.37 0.92
下载: 导出CSV
表 4 渗透比降统计表
Table 4 Statistical list of infiltration ratio
工况 渗透比降 防渗墙 心墙 坝体代料 下游逸出处 初始状态 24.53 0.11 0.13 0.14 一 37.38 0.22 0.14 0.11 二 16.53 0.26 0.13 0.11 三 18.58 0.24 0.13 0.10 四 19.02 0.21 0.13 0.11 五 37.38 0.25 0.13 0.11 六 38.20 0.23 0.13 0.10 七 39.13 0.21 0.13 0.10
下载: 导出CSV
表 5 心墙坝抗滑稳定安全系数统计表
Table 5 Statistical list of safety factor of anti-sliding stability of core dam
工况 上游坡最小抗滑稳定安全系数 降幅/% 下游坡最小抗滑稳定安全系数 增幅/% 降幅/% 工况一 2.14 6.56 1.54 — 21.08 工况二 1.53 33.19 1.97 0.96 — 工况三 1.54 32.76 1.97 0.96 — 工况四 1.55 32.32 1.97 0.96 — 工况五 1.39 39.31 1.59 — 18.52 工况六 1.41 38.43 1.58 — 19.03 工况七 1.40 38.87 1.57 — 19.54
下载: 导出CSV
-
[1] 张建云. 气候变化与水利工程安全[J]. 岩土工程学报, 2009, 31(3): 326–330. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract13166.shtml ZHANG Jian-yun. Climate change and safety of water conservancy projects[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(3): 326–330. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract13166.shtml
[2] 张建云, 向衍. 气候变化对水利工程安全影响分析[J]. 中国科学: 技术科学, 2018, 48(10): 1031–1039. https://www.cnki.com.cn/Article/CJFDTOTAL-JEXK201810002.htm ZHANG Jian-yun, XIANG Yan. Analysis on the impact of climate change on the water conservancy project safety[J]. Scientia Sinica (Technologica), 2018, 48(10): 1031–1039. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JEXK201810002.htm
[3] 徐旺敏, 傅琼华, 马秀峰, 等. 滨田水库大坝滑坡机理分析及防治对策[J]. 人民长江, 2012, 43(20): 36–38. https://www.cnki.com.cn/Article/CJFDTOTAL-RIVE201220012.htm XU Wang-min, FU Qiong-hua, MA Xiu-feng, et al. Analysis of landslide mechanism of Bintian reservoir dam and countermeasures[J]. Yangtze River, 2012, 43(20): 36–38. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-RIVE201220012.htm
[4] HUANG Q X, WANG J L, XUE X. Interpreting the influence of rainfall and reservoir infilling on a landslide[J]. Landslides, 2016, 13(5): 1139–1149. doi: 10.1007/s10346-015-0644-8
[5] 黄润秋. 20世纪以来中国的大型滑坡及其发生机制[J]. 岩石力学与工程学报, 2007, 26(3): 433–454. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200703000.htm HUANG Run-qiu. Large-scale landslides and their sliding mechanisms in China since the 20th century[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(3): 433–454. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200703000.htm
[6] RAHARDJO H, ONG T H, REZAUR R B, et al. Factors controlling instability of homogeneous soil slopes under rainfall[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2007, 133(12): 1532–1543. doi: 10.1061/(ASCE)1090-0241(2007)133:12(1532)
[7] 唐栋, 李典庆, 周创兵, 等. 考虑前期降雨过程的边坡稳定性分析[J]. 岩土力学, 2013, 34(11): 3239–3248. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201311032.htm TANG Dong, LI Dian-qing, ZHOU Chuang-bing, et al. Slope stability analysis considering antecedent rainfall process[J]. Rock and Soil Mechanics, 2013, 34(11): 3239–3248. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201311032.htm
[8] 李卓, 何勇军, 李宏恩, 等. 前期降雨作用下边坡滑坡模型试验[J]. 河海大学学报(自然科学版), 2016, 44(5): 400–405. https://www.cnki.com.cn/Article/CJFDTOTAL-HHDX201605004.htm LI Zhuo, HE Yong-jun, LI Hong-en, et al. Model test on slope landslides under antecedent rainfall[J]. Journal of Hohai University (Natural Sciences), 2016, 44(5): 400–405. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HHDX201605004.htm
[9] JIA G W, ZHAN T L T, CHEN Y M, et al. Performance of a large-scale slope model subjected to rising and lowering water levels[J]. Engineering Geology, 2009, 106(1/2): 92–103.
[10] JIAN W X, XU Q, YANG H F, et al. Mechanism and failure process of Qianjiangping landslide in the Three Gorges Reservoir, China[J]. Environmental Earth Sciences, 2014, 72(8): 2999–3013.
[11] 李卓, 何勇军, 盛金保, 等. 降雨与库水位共同作用下近坝库岸边坡滑坡模型试验研究[J]. 岩土工程学报, 2017, 39(3): 452–459. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract16853.shtml LI Zhuo, HE Yong-jun, SHENG Jin-bao, et al. Landslide model for slope of reservoir bank under combined effects of rainfall and reservoir water level[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(3): 452–459. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract16853.shtml
[12] SUN G H, YANG Y T, CHENG S G, et al. Phreatic line calculation and stability analysis of slopes under the combined effect of reservoir water level fluctuations and rainfall[J]. Canadian Geotechnical Journal, 2017, 54(5): 631–645.
[13] 碾压式土石坝设计规范: SL274—2020[S]. 2020. Design Code for Rolled Earth-Rock Fill Dams: SL274—2020[S]. 2020. (in Chinese)
[14] Geo-Slope International Ltd. Seepage Modeling with SEEP/W 2007[R]. Calgary: Geo-Slope International Ltd, 2007.
[15] FREDLUND D G, MORGENSTERN N R, WIDGER R A. The shear strength of unsaturated soils[J]. Canadian Geotechnical Journal, 1978, 15(3): 313–321.
[16] Geo-Slope International Ltd. Stability Modeling with SLOPE/W 2007[R]. Calgary: Geo-Slope International Ltd, 2007.
-
期刊类型引用(8)
1. 李珍玉，单世杰，林坤，朱航. 微生物促进矿物相变及其改良膨胀土胀缩特性试验研究. 岩石力学与工程学报. 2025(01): 209-220 . 
百度学术
2. 朱文羲，邓华锋，李建林，马林建，李锦瑞，陈勇琪，陈向阳. 干湿循环作用下微生物改良花岗岩残积土劣化规律研究. 岩石力学与工程学报. 2025(02): 482-491 . 百度学术
3. 史金权，付贵永，刘汉龙，肖杨. 微生物加固模拟月壤强度特性试验研究. 土木与环境工程学报(中英文). 2025(02): 20-29 . 百度学术
4. 胡雯璐，刘鹏. 脲酶喷洒工艺对赤泥矿化胶结及抑尘效果试验. 林业工程学报. 2025(02): 173-179 . 百度学术
5. 汤禹，付俊，陈安，周罕，张宇，罗磊. 改进灌注方式下MICP固化尾矿中重金属形态特征及风险评价. 中国环境科学. 2025(03): 1385-1394 . 百度学术
6. 柳嘉豪，林文彬，卓祖磊，吴杉颖，高玉朋，张佳源，罗承浩. 微生物诱导碳酸钙沉淀技术固化海相淤泥试验. 福建理工大学学报. 2025(01): 33-39 . 百度学术
7. 王东星，许凤丽，泮晓华，商武锋，吴章平，郭克诚. GGBS-MICP协同固化淤泥质砂土工程特性研究. 岩石力学与工程学报. 2025(05): 1349-1362 . 百度学术
8. 何文杰，郑文杰，谢毅鑫，薛中飞，秦鹏，吕鑫江. 基于纳米羟基磷灰石的矿化技术修复铅污染水和一维土柱的试验研究. 土木工程学报. 2024(11): 45-56 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 232
- HTML全文浏览量: 79
- PDF下载量: 64
- 被引次数: 11
