《某水库库区渗流场及其水文地质参数敏感性分析》由会员分享，可在线阅读，更多相关《某水库库区渗流场及其水文地质参数敏感性分析（3页珍藏版）》请在金锄头文库上搜索。
1、某水库库区渗流场及其水文地质参数敏感性分析摘要：库区渗流场的分布特征，很大程度上决定了其工程特性，尤其是岩溶区 的发育特征，其是评价库区水文地质条件的一项重要指标，同时定量化分析其水 文地质参数的敏感性也是难题之一。根据实际工程应用的需要，研究了库区渗流 场的地下水等水位，库区地下水分水岭与地表水分水岭基本一致，地下水总体上 由上库向四周补给，由于南北两侧断层的阻水作用，地下水在南侧和北侧以泉的 形式排泄；同时，由水文地质参数敏感性分析可知，第IV渗透性分区渗透参数的 变化对地下水渗流场变化有较大的影响，其它区域渗透参数变化时，地下水渗流 场变化不明显，即地下水渗流场对其它区域渗透参数变化不敏
2、感。关键词：水文地质；渗流场；数值计算；敏感性分析1引言库区地下水渗流分析多以线性的达西定律为基础，这是因为在大多数情况下， 地下水渗流是满足或近似满足达西定律的，而且达西定律所表达的线性关系也使 理论分析与数值分析较为简洁，然而修建在岩溶地区的大坝基础中可能存在强透 水地质构造，渗流的雷诺系数很大，渗流速度与水力坡度不再遵从达西定律，此 时必须考虑惯性力的影响。柴军瑞（2001年）1、代群力（2000年）2认为有 必要对岩溶区进行坝基地下水非线性渗流分析。由于非线性渗流的复杂性，至今 仍没有一个统一的公式来更好地描述这种流动，其中两个常用的经验公式模型为 Forchheimer 公式3。19
3、65 年，津克维茨（Zienkiewiz）和张（Cheung）提出有限 单元法适用于所有可按变分形式进行计算的场问题，为该方法在渗流分析中的应 用提供了理论基础。Leiws R W等（1987年）将这一方法用于模拟多孔介质中 的地下水。国内不少研究5-6 认为有限单元法是一种先进有效的数值模拟法，用 于渗流分析计算时可部分替代模型试验，精度相对较高，可方便地模拟多种外部 条件的特点。2工程地质环境东部某水库总库容1676万m3,北邻长江，南望太湖，属宁镇山脉的低山丘 陵区，地势总体是西高东低，地形切割不深，沟谷稍发育，以冲积堆积、剥蚀地 貌为主，区内普遍可见数级剥夷面和34级阶地，以及相应的成
4、层岩溶、深切 河谷。工程区从岩溶发育程度、地层岩性、地质构造及岩溶水赋存运移特征来划 为四类岩溶水文地质单元（区），一类非可溶岩分布区；其中：A区一裸露型纯 碳酸盐岩单斜岩溶发育强度弱中等水文地质单元（区），地表基岩大片裸露， 岩性由白云岩类逐渐变为灰岩类；西侧以F7断层为界，北侧以F4断层为界，东 侧为仑山东面坡脚下水库库区F12-1断层及其伴生的闪长玢岩脉和龙潭组（P2l），南侧以仑山南面坡脚奥陶系上统（O3）为界，位于工程区上水库、输水 发电系统枢纽工程部位，为本文渗流场分析区；B区一覆盖型纯碳酸盐岩与非可 溶岩互层（或夹层）岩溶发育强度中等水文地质单元，分布在仑山坡脚，下水库 西侧库岸
5、；C区一裸露型背斜岩溶发育强度弱中等水文地质单元，位于上水库 的西侧；非可溶岩由粉砂岩、粗面岩、闪长玢岩等组成，主要分布在下水库的北 侧及仑山的南侧坡脚仑山水库一带。3渗流场模型及结果分析3.1数学模型的概化由于枯水期地下水变化不大，因此可以将地下水运动概化为稳定流。该计算 区域内的地下水运动可以概化为在一个非均质各向同性介质中的二维流动。3.2几何模型的建立及参数的选取根据地层岩性、岩溶水文地质单元，结合地形地貌特征，将研究区分为五个 区。北部边界沿着F3断层，南部边界至F9断层，东部以姊妹桥F12断层附近钻 孔25所在一线为边界，西部以F7为边界。I区：F4断层与F3断层之间的志留系坟头组
6、的砂岩、粉砂岩；H区：F7断层 与F8断层之间北至断层F4,南至仑山山脊；B区：由断层F12与断层F8与仑山 山脊围成的三角形区域；IV区：F7断层与F8断层之间，仑山山脊以南的部分； V区：F8断层以东，仑山山脊以南部分。为了充分利用观测资料，将上库区的钻孔地下水水位作为上库第一类边界条 件；下库区泉水出露的点作为流量边界。参数的选取采用钻孔压水试验换算值和经验值相结合的方法。根据钻孔压水 试验资料，渗透I区渗透系数值取0.0175m/d；渗透II区渗透系数取为 0.00864m/d ；渗透B区渗透系数值取0.0195m/d ；渗透V区渗透系数值取 0.021m/d ；渗透V区渗透系数值取0
7、.0187m/d。考虑库盆中降雨入渗补给，根据 仑山水库降雨量资料，一月份平均降雨量为46.1mm，艮口 1.54mm/d。在库盆中地 表水绝大部分渗入地下，因此降雨入渗补给系数取0.55。3.3数值模拟结果及分析根据上述数学模型，采用FORTRAN编写的二维有限元程序进行了模拟计算， 库区渗流场等水位线。在上库的主体部分，其西北、西及西南三侧分水岭均位于 本亚区。从地下水等水位线图中可以看出，地下水分水岭与地表水分水岭基本一 致。地下水总体上由上库向四周补给，由于南北两侧的非可溶岩（微透水层）的 阻水作用，地下水在南侧和北侧以泉的形式排泄。仑山山脊一带存在一个分水岭， 地下水在总体上向东流的
8、同时，还向分水岭的南、北两侧流动。在上库南侧，地 下水向大哨泉以北的冲沟中排泄；东侧地下水向正东排泄至姊妹桥方向；在北侧， 存在一近南北向的分水岭，地下水向北流动的同时分别向东、西两侧排；在西侧， 地下水向钥匙湾流动，这与地表水的流向基本一致。为了更好的表明计算的精度，分别在五个区各设立了两个地下水水位观测点， 共计十个地下水水位观测点，现将计算值与实测值进行对比。计算结果与实际测 量值较为接近，其中最大误差未超过5m，分别出现在I2和B1处；而最小误差 为2m，出现在V2监测点，最大相对误差为2.6%，可见模型具有一定的精度， 同时证明该方法具有一定的正确性和实用性。4模型水文地质参数敏感性
9、分析由于研究区地质及水文地质研究处于进一步研究阶段，许多参数还不确定， 因此可以通过参数敏感性分析来研究不同参数时的地下水运动状态，分析参数对 渗流场的敏感性。分别将5个分区中的某一区渗透系数加倍或者减小一倍，而进 行计算，结果显示：第V渗透性分区渗透参数的变化对地下水渗流场变化有较大 的影响，其它区域渗透参数变化时，地下水渗流场变化不明显，即地下水渗流场 对其它区域渗透参数变化不敏感。通过对研究区渗流场进行计算分析可知，ZK1钻孔以东仑山南北坡之间仍存 在一分水岭，地下水在向分水岭南、北两侧流动的同时，总体上还向东流动。仑 山南北坡由于非可溶岩的阻水作用，地下水在南侧和北侧以泉的形式排泄。本
10、次 计算中没有单独划分这一高渗透性区域，因而计算结构没有反映低水位槽的情况， 渗流计算仅给出了研究区地下水运动的总体趋势。F8断层以东的区域渗透性参数变化对地下水渗流场有较大的影响，其余区域 渗透参数变化对地下水渗流场不敏感。5结论（1）从地下水等水位线图中可以看出，地下水分水岭与地表水分水岭基本一 致。地下水总体上由上库向四周补给，由于南北两侧非可溶岩的阻水作用，地下 水在南侧和北侧以泉的形式排泄；东侧地下水向正东排泄至姊妹桥方向；在西侧， 地下水向钥匙湾流动。（2）通过参数反演可知，第IV渗透性分区渗透参数的变化对地下水渗流场变 化有较大的影响，其它区域渗透参数变化时对地下水渗流场影响不大
11、。参考文献：1 柴军瑞.坝基非达西渗流分析J.水电能源科学，2001，19（4）： 1-32 代群力.地下水非线性流动模拟J.水文地质工程地质，2000（2）： 50-513 Bear J. Hydraulics of Ground waterM. New York： McGraw-Hill， 19794 Leiws R W， Schrefler B A. The finite element method in the deformation and consolidation of porous media M.New York： John Wiley，19875 李康宏，柴军瑞.坝基渗流分析两种数值方法的比较J.红水河，2003, 22 （4）：14-176 闫澍旺，王瑞钢，王学军，李少明.白石水库渗流场有限元分析J.水力发电学报，2003，2： 53-61作者简介：杨丹，女，1980年6月出生，河南周口人，学士，高级工程师，主要从事水 文水资源与水文地质研究。
