复杂层状土层地质对堤岸渗流特性的影响研究_王海龙.pdf

：复杂层状土层地质对堤岸渗流特性的影响研究收稿日期：作者简介：王海龙（），男，工程师，从事资源勘查工程工作 通信作者：王一朵（），女，在读硕士，从事岩土工程研究王海龙，杜宁宇，王一朵（江苏省地质局，江苏 南京；河海大学河海里尔学院，江苏 南京；河海大学土木与交通学院，江苏 南京）摘 要：随着大坝构筑物的发展，安全问题也时有发生。渗流问题是关系到水库大坝安全运行的主要因素。而现场监测仅仅是对断面局部测压管水位的监测，具有一定的局限性，根据合肥市大房郢水库大坝工程的现场监测结果，采用数值模拟对现场监测数据的扩展，可以模拟现场不同工况下大坝内部渗流场的变化，分析不同地质条件对大坝渗流及稳定性的影响。结果表明，与均质土层坝基相比，层状土层坝基测线处的孔隙水压力减小，浸润线高程降低 ，上下游坝坡稳定安全系数分别提高了 ，；在坝基较上层下卧阻水层内水力梯度随着高程的增大而增大，且下卧阻水层厚度越大，水力梯度增大幅度越大；在非饱和区域内，流速均随高程的增大而降低，在饱和区域内，流速的变化趋势主要受水力梯度及不同土层渗透系数的影响。关键词：层状土层；水库大坝；二维；数值模拟；渗流中图分类号：文献标识码：文章编号：（）引言随着大坝构筑物的发展，大坝渗流及稳定性研究得到越来越多的重视。水库建成蓄水后，在上、下游水位作用下，坝体和坝基均会出现渗流现象，进而影响坝坡稳定性。韩文娟通过建立均质坝体模型，研究了坝体土料不同渗透系数在水位下降条件下的坝体浸润线和上游逸出点渗透坡降的变化规律，岳庆河等深入分析了在水位下降时坝体材料渗透系数大小以及渗透系数的各向异性对大坝渗流及稳定的影响，李子阳基于饱和 非饱和渗流理论，考虑了渗透系数与基质吸力的非线性关系，研究大坝非稳定渗流场和渗透稳定性的变化规律，坝体材料的渗透系数大小及渗透系数的各向异性会不同程度地影响大坝的渗流及稳定性，为此钟学梅对两个影响因素进行敏感性分析，在此基础上，王宁研究了坝体土体参数对大坝渗流及稳定的影响，包括土体渗透系数大小以及强度参数。在实际工程中，坝基的地质条件并不是单一均质的，以上众多学者 在研究大坝渗流分析模型中，仅将坝基作为同一种材料来考虑，忽略了坝基的复杂地质情况。以合肥市大房郢水库为例，该工程沿坝轴线方向的长度较长，沿坝轴线方向的土层分布也发生变化，不同断面位置的地质条件存在差异，如分层土层种类、土层厚度以及土层起伏等。由于各土层渗透性的差异，不同土层分布将对大坝渗流产生不同的影响。因此，建立复杂的坝基地质条件模型更加符合实际情况。本文以合肥市大房郢水库大坝为依托，基于现场监测结果，通过数值模拟研究不同库水位下复杂地质条件对大坝渗流及稳定性的影响。本文采用饱和 非饱和渗透模拟分析，根据饱和 非饱和渗流理论，开展二维有限元数值模拟，通过二维数值模拟，得出复杂地质对大坝渗流影响的一般规律。模拟几何模型图 展示了合肥市大房郢水库工程堤岸的计算模型。203932.520119.58迎水坡 13cdbahg34.00背水坡 12.5e1粉质黏土（21 m）21.00f119.5坝身填土203932.520834.00迎水坡 13背水坡 12.51粉质黏土2砂壤土1全风化砂岩21.00216105203932.5208119.534.00迎水坡 13背水坡 12.5坝身填土1粉质黏土2砂壤土1全风化砂岩21610521.00（a）粉质黏土 21 m（b）粉质黏土 5 m，砂壤土 10 m（c）粉质黏土 10 m，砂壤土 5 m单位：m图 1二维几何模型坝身填土第 卷 第 期 年 月 山西建筑 图 中坝高为 ，上游坝坡比为 ，下游坝坡比为 ，上下游坝坡脚外延伸的长度为 ，将下游坝脚外延伸的长度分别取，以研究模型尺寸对渗流结果的影响，结果表明，当下游坝脚外延伸的长度为，时的渗流场变化一致，因此在不影响大坝渗流场的情况下，取最小下游坝脚外延伸长度为 。土体参数大坝为碾压式土坝，坝体材料的各向异性比值为，即坝身填土水平向渗透系数与竖直向渗透系数比值为，坝基老填土的水平向渗透系数与竖直向渗透系数比值为 ，坝基其余土层水平向渗透系数与竖直向渗透系数相同。大房郢水库大坝各土层的土体物理力学性质指标见表。表 大房郢水库土层分区及主要参数土层名称土层性质黏聚力 内摩擦角（）天然密度（）孔隙率 渗透系数（）坝身填土 重粉质壤土 粉质黏土 砂壤土（含轻粉质壤土）全风化砂岩 由主坝纵剖图及各土层物理力学性质指标可知，大房郢水库大坝地面线以上即坝体部分土层较为一致，均为渗透系数较小的填土层，地面线以下土层均含有老填土层，其渗透系数大于填土层的渗透系数，而老填土层以下的坝基土层较为复杂，坝基地质条件复杂性体现在坝基含有 土层不等，且各土层的位置和厚度不同，它们的渗透系数差异较大，而坝基的最底层即相对不透水层均为层全风化砂岩，该土层的渗透系数很小。综上可将此复杂土层简化为三个断面进行计算。边界条件二维边界条件如图（）所示，图中 为水头边界条件，为下游坡面，设置为自由渗出面，坝基底部 及坝基两侧设置为不透水边界条件，此边界的流量为零，下游边界 为自由渗出面，设置为流量边界条件。模拟结果分析 坝体浸润线变化坝基不同工况下坝体浸润线如图 所示，图 中虚线表示大坝中心线即测线的所在位置。模拟编号 坝基为均质粉质黏土 ，模拟编号 坝基为粉质黏土 砂壤土 的层状土层，模拟编号 坝基为粉质黏土 砂壤土 的层状土层。坝基为均质土层与层状土层的坝体浸润线有一定差异，坝基为均质粉质黏土的坝体浸润线最高，在层状土层坝基中，坝基较上层粉质黏土厚 的坝体浸润线高于粉质黏土厚 的坝体浸润线。由图 测线处孔隙水压力变化可知，测线处均质粉质黏土坝基的浸润线高程较高，坝基为粉质黏土厚 的浸润线高程次之，坝基为粉质黏土厚 的浸润线高程最低。具体表现为：测线处均质粉质黏土的坝体浸润线高程比层状土层坝基中粉质黏土厚 的坝体浸润线高 ；层状土层坝基中，粉质黏土厚 的坝体浸润线比粉质黏土厚 的坝体浸润线高 。表明与均质土层坝基相比，层状土层坝基的浸润线高程低 左右，在层状土层坝基中，坝基较上层的下卧阻水层厚度增大、下卧导水层厚度减小，坝体浸润线升高。当下卧阻水层厚度增大 、下卧导水层厚度减小 时，坝体浸润线升高 。302826242220-30-20-10010203040水位高程/m上游水位模拟编号 1（粉质黏土 21 m）模拟编号 2（粉质黏土 5 m，砂壤土 10 m）模拟编号 3（粉质黏土 10 m，砂壤土 5 m）轴距/m图 2坝基不同工况下坝体浸润线位置 水力梯度变化图 分别展示了均质粉质黏土坝基、粉质黏土厚 的层状土层坝基以及粉质黏土厚 的层状土层坝基的水力梯度的分布规律，其水力梯度最大值均位于坡脚处，小于临界水力梯度，故都不会发生渗透破坏。XY-水梯度00.050.050.10.10.150.150.20.20.250.250.30.30.350.350.40.40.4500.10.10.20.20.30.30.40.40.50.50.60.60.7XY-水梯度00.10.10.20.20.30.30.40.40.50.50.60.60.70.70.80.80.9XY-水梯度（a）模拟编号 1（粉质黏土 21 m）（b）模拟编号 2（粉质黏土 5 m，砂壤土 10 m）（c）模拟编号 3（粉质黏土 10 m，砂壤土 5 m）图 3坝基不同工况下水力梯度云图坝基不同工况下测线处水力梯度随高程的变化如图 所示，图 中虚线为高程 位置处，即坝体和坝基的分界线。均质土层、层状土层及其不同厚度的坝基，其水力梯度的分布有一定的差异。坝基不同工况下测线处的水力梯度均未超过各土层的临界水力梯度，均质粉质黏土坝基测线处最大水力梯度为，层状土层坝基中，粉质黏土厚 的坝基，其测线处最大水力梯度为 ，粉质黏土厚 的坝基测线处最大水力梯度为 。表明均质粉质黏土坝基测线处最大水力梯度最大，坝基较上层下卧阻水层厚度越大，水力梯度越大，当下卧阻水层厚度增大 时，测线处最大水力梯度增大。对于层状土层坝基，在坝体土层范围内，水力梯度随高程的增大呈现先增大后减小的趋势，具体表现为：坝基为粉质黏土厚 的水力梯度由 增大至 后逐渐降低至；坝基为粉质黏土厚 的水力梯度由 增大至 后逐渐降低至 。在坝基土层范围 第 卷 第 期 年 月 王海龙等：复杂层状土层地质对堤岸渗流特性的影响研究内，层状土层坝基的水力梯度的变化趋势与均质土层坝基有所不同，其水力梯度的变化趋势呈现两段式，在坝基较下层即全风化砂岩和砂壤土层内水力梯度不变，而在坝基较上层下卧阻水层即粉质黏土层内，水力梯度随着高程的增大而增大，下卧阻水层厚度为 的坝基，其水力梯度增大了 ，下卧阻水层厚度为 的坝基，其水力梯度增大了 ，表明下卧阻水层厚度越大，水力梯度增大幅度越大。图 4测线水力梯度随高程的变化4540353025201510500.020.040.060.080.100.12水力梯度高程/m坝体坝基模拟编号 1（粉质黏土 21 m）；模拟编号 2（粉质黏土 5 m，砂壤土 10 m）；模拟编号 3（粉质黏土 10 m，砂壤土：5 m）流速变化坝基不同工况下测线处流速随高程的变化如图 所示，图 中虚线为饱和区与非饱和区的分界线。对于层状土层坝基，在坝基较上层粉质黏土层内，土体渗透系数不变，由于水力梯度随高程的增大而增大，因而流速也随之增大。而在坝基较下层即砂壤土以及全风化砂岩土层内，水力梯度不变，渗透系数不变，水在各土层内的流速保持不变，但由于受到土层分层渗透系数发生变化的影响，在不同土层交界处流速会发生突变。具体表现为：对于粉质黏土厚 的层状土层坝基，在 ，饱和区域内，流速保持不变，在 处由于土层渗透系数变大，流速发生突变，在 粉质黏土层内，流速随水力梯度的增大而呈折线型上升，在 处由于土层渗透系数变大，流速发生突变。对于粉质黏土厚 的层状土层坝基，在 ，饱和区域内，流速保持不变，在 处由于土层渗透系数变大，流速发生突变，在 粉质黏土层内，流速随水力梯度的增大而呈折线型上升，在 处由于土层渗透系数变大，流速发生突变。对于不同厚度的层状土层坝基，坝基较上层粉质黏土厚 的流速略大于粉质黏土厚 的流速，表明坝基较上层下卧阻水层的厚度越厚，其测线处流速也越大，当下卧阻水层的厚度增大为 时，测线处最大流速为 。结论本文根据饱和 非饱和渗流理论，开展二维有限元数值模拟，研究库水位不变条件下复杂地质二维大坝的坝体浸润线、水力梯度以及流速矢量的变化，同时将数值模拟结果与现场测压管水位数据进行比较，分析不同地图 5测线流速随高程的变化454035302520151050高程/m非饱和区饱和区1E-111E-101E-091E-81E-7流速/（m s-1）模拟编号 2（粉质黏土 5 m，砂壤土 10 m）；模拟编号 3（粉质黏土 10 m，砂壤土 5 m）模拟编号 1（粉质黏土 21 m）；质条件对坝坡稳定性的影响。主要结论如下：）与均质土层坝基相比，层状土层坝基测线处的孔隙水压力值小 ，浸润线高程降低 ，上下游坝坡稳定安全系数分别提高了 ，；对于不同厚度的层状土层坝基，坝基较上层下卧阻水层厚度越厚，坝体浸润线越高，坝坡稳定性越差。当下卧阻水层厚度增大 、下卧导水层厚度减小 时，坝体浸润线高程升高 ，上下游坝坡安全系数分别降低了，。）在坝基土层范围内水力梯度有明显差异，具体层状坝基水力梯度的变化趋势呈现两段式，在坝基较下层水力梯度不变，在坝基较上层下卧阻水层内水力梯度随着高程的增大而增大，且下卧阻水层厚度越大，水力梯度增大幅度越大，当下卧阻水层厚度增大 时，测线处最大水力梯度增大 。）对于层状土层坝基，在坝基较上层粉质黏土层内，土体渗透系数不变，由于水力梯度随高程的增大而增大，因而流速也随之增大。在饱和区域内，流速的变化趋势主要受水力梯度及不同土层渗透系数的影响，当土层的渗透系数发生变化时，流速在土层交界处发生突变。对于不同厚度的层状土层坝基，坝基较上层下卧阻水层的厚度增大，其测线处流速也变大，当下卧阻水层的厚度增大为 时，测线处最大流速为 。参考文献：原有芬 双龙水电站心墙土石坝浸润线分析东北水利水电，（）：，（）：韩文娟 库水位降落时均质土坝非稳定渗流分析 扬州：扬州大学，岳庆河 基于 的土石坝渗流与稳定分析研究 泰安：山东农业大学，李子阳，马福恒，张湛，等 水位骤变下土石坝非稳定渗流及稳定分析 人民黄河，（）：钟学梅 库水位升降条件下土石坝渗流稳定分析及应用研究 成都：西华大学，第 卷 第 期 年 月 山西建筑 王宁 土石坝渗流影响因素及坝坡稳定性研究 成都：西华大学，何秉顺，丁留谦 各向异性渗流对堤坝稳定性的影响 中国水利水电科学研究院学报，（）：孙旭曙，梁佳灯，陈超，等 均质土坝正交各向异性渗流场数值模拟 水科学与工程技术，（）：王亚月 水位下降对均质土石坝坝坡稳定的影响分析 石河子：石河子大学，（，；，；，）：，：；（上接第 页）李磊，李寿君 基于 交替法求解分析浅埋近距离煤层开采覆岩结构演化规律 采矿与岩层控制工程学报，（）：皇甫鹏鹏，张路青，伍法权 基于高效保角变换的单个复杂洞室围岩应力场解析解研究 岩石力学与工程学报，（增刊）：王润富 一种保角映射法及其微机实现 河海大学学报，（）：吕爱钟，覃 媛，陈虹宇 马蹄形隧洞考虑支护滞后过程的应力分析 岩土力学，（增刊）：，（）：邓涛，魏雯，关振长，等 公路隧道开挖问题的复变函数解及应用 中国公路学报，（）：吕爱钟 地下洞室最优开挖形状的确定方法岩石力学与工程学报，（）：刘建友，吕刚，赵勇，等 应力流守恒原理及地下洞室群支护结构设计方法 隧道建设（中英文），（增刊）：，（，）：，：，（），（），：；第 卷 第 期 年 月 王海龙等：复杂层状土层地质对堤岸渗流特性的影响研究