不同水位下10m水面差下的杨家湾险段渗透系数模拟计算.docxVIP

不同水位下10m水面差下的杨家湾险段渗透系数模拟计算 1 荆江保水开采中的地下渗流场 河流治理是一种危害水库安全的重要危险。研究河流治理的成因对预防和控制河流控制具有重要的现实意义。农业灌溉引起地下水位上升可以导致管涌，工程建设中的开挖也可以导致管涌。但最常见的、且危害性最大的管涌是在江河持续高水位下土层发生渗透变形引起的管涌。对于管涌形成的机制，有许多文献进行了讨论。文认为1998年长江防讯期间堤坝及其附近出现管涌的原因是坝身结构和地基构造存在缺陷。文详细研究了黄河中游的河床覆盖层渗透变形特征。文认为管涌也是一种地貌过程。 了解在高水位条件下大堤附近的地下渗流场状况更具有实际意义。文从理上论讨论了有限元方法在分析地下渗流场的方法。文用有限元方法模拟了长江安庆城市防洪堤典型地段地下水渗流场，并由此提出了恢复减压井效力的工程建议。 长江中游的荆江河段，每逢长江洪水，荆江大堤必出现管涌险情。虽然经过多次勘察设计、加固治理，但是在20世纪90年代的几次大洪水中，荆江大堤仍然出现了大量的管涌等灾害性险情。监利县杨家湾险段(桩号：638+750～642+500)的荆江大堤有一巨大的拐弯(图1)，堤外为高漫滩，1998年和1999年发生管涌处位于人民大垸围堤与荆江大堤交汇处。外滩宽阔而平坦，北略高于南，地面高程28～31 m，相对高差3 m，局部地势低洼，常年积水。大堤外脚处有宽100～200 m平土坑，坑内常年积水，水深1～3 m。堤内为冲积平原，平坦开阔。 本段堤基土质的工程地质条件极差，上覆为近代湖沼相沉积的淤泥或淤泥质土类，厚度0～5 m；其下为粘性土及砂砾石。 由于堤基土质较差，杨家湾险段历史上发生多次管涌和溃口；该堤段经过多次加固，但仍在1998和1999年的大洪水时出现管涌等重大险情。目前，有关部门正在或准备采用工程措施加固提防，防止土层渗透变形。为此，本文采用美国地质调查局(USGS)的Visual Modflow软件，模拟洪水期长江高水位条件下，不同防渗处理模式下地下水渗流场的变化特征，为实际工程处理提供理论及方法上的指导。 2 w软件模拟 采用美国地质调查局(USGS)的Visual Modflow软件进行模拟。搜集荆江大堤钻孔剖面资料，在综合考虑各剖面的工程地质及水文地质性质的基础上，选择了监利县杨家湾的011横剖面作为模拟对象，建立稳定流条件下潜水剖面二维流模拟模型。 2.1 砂性土壤厚度 该剖面(图2)上层粘土和人工填土厚度为5 m左右，下伏2～9 m厚壤土，砂壤土厚度约为3 m；再往下为粉砂层和细砂层，其总厚度为25 m左右，其间夹有一壤土凸镜体(其主轴直径约3 m)，长江深切下部含水层，与细砂层贯通，水力联系密切。堤身为人工填土。各层渗透系数见表1所列。 2.2 ．堤内侧上边界 定水头边界：将江水与其下伏地层的地形交切线处理为定水头边界，下边界处理为零流量边界，横剖面江水位的分水岭处亦处理为零流量边界，在洪水期由于地下水位可达到地表，堤内侧上边界亦处理为定水头边界。 堤内上边界水头的确定：以堤内侧最低点的地面标高作为基准面，赋零水头值，其他各点的水头值同理可用线性插值的方法来估算。如果堤内有渊塘，以渊塘的水位作为基准面，赋零水头值，渊塘水位以上部分的上边界水头值根据地形标高估算。 堤外上边界水头的确定：以江水位(或设计水位)相对于基准面的水头差来定。在模拟过程中可以根据长江水位的变化对水头差调整，模拟不同水头差下(3,7,10,13 m)堤下地下水渗流场的特征。 2.3 地层剖分叠加 011横剖面用软件剖分为44行76列，堤身部分的列进行了手工剖分加密。根据渗流等效原理，通过改变渗透系数的办法来解决纵横比例尺不一致问题。 2.4 参数的确定和输入 渗透系数值和孔隙度值为实测值，有效孔隙度、各向异性系数以及重力给水度等参数根据经验数据确定(表1)。根据岩性分区输入参数。 3 斜拉合压井深度对边界流速的影响 根据该剖面不同模拟模式下边界流速及流网分布图(图3～8)，可以获得以下几个结果： 在改变水头差(3,7,10 m)模式下(图3)，随江水位的升高，边界流速几乎成倍增加，意味着出流带在相同的岩性条件下，水力梯度成倍增加，管涌危险性亦增加。对于同一水头差条件下，模型上界(出流带)垂直流速距堤脚由近到远逐渐减小，其数值较小而且变化不大，从等水头线可以看出水头主要损失于上层人工填土；模型左边界水平流速和岩性密切相关，渗透系数大的其流速也大，而且数值变化较大，最大与最小的差5个数量级，水流主要从下层的细砂层和粉砂层通过。 在高水位(10 m水头差)减压井处理模式下(图4)，在减压井达到下部含水层时，减压井的深度对边界(出流带)流速的影响较小(图3与图4对比)。随减压井深度的增加，减压井左侧渗流场的等水头线向减压井附近