弱透水层低速非达西流咸水下移过程的模拟研究 Numerical study on saltwater downward migration in aquitard as low velocity non-Darcy flow.pdfVIP

水 利 学 报 2014年7月 SHUILI XUEBAO 第45卷 第7期 文章编号：0559-9350（2014）07-0875-08 弱透水层低速非达西流咸水下移过程的模拟研究 崔莉红，成建梅，路万里，李敏敏 （中国地质大学（武汉）环境学院，湖北 武汉 430074） 摘要：弱透水层以黏性土为主，其渗流通常表现出低速非达西特性。开采条件下，弱透水层中咸水下移问题是典 型的渗流-固结-溶质运移耦合问题。本文以华北平原典型的“上咸下淡”含水层系统为研究对象，采用修正的多 场耦合模型，研究了弱透水层中非达西渗流固结过程对咸水下移的影响。研究结果表明，非达西渗流不仅增加了 黏性土层的固结时间，而且，在一定程度上，减缓了咸水下移速度。同时，非达西渗流条件下流速具有削峰拖尾 特征，天然水力梯度越大，拖尾越不明显，而削峰始终存在。咸水下移前期，不考虑固结过程时其运移速度更 快，而在中后期，情况相反。阻滞因子决定了弱透水层渗流固结过程对咸水下移的影响程度。 关键词：弱透水层；非达西渗流；渗透固结；数值模拟；咸水下移 中图分类号：P641.2 文献标识码：A doi：10.13243/ki.slxb.2014.07.015 1 研究背景 华北平原水文地质条件的差异及历史成因使得相邻含水层普遍存在着“上咸下淡”的水质格局， 即在浅层广泛分布着一层咸水，而下部的多层含水层中是优质的淡水。由于深部地下淡水的不合理 ［1］ 开采，上下含水层形成巨大的水头差，可能促使咸淡水界面下移 。位于含水层之间的弱透水层，是 联系上下含水层的纽带，其岩性常以黏性土层为主，咸淡水界面位于此层。受下部地下水开采影 ［2-5］ 响，上下水头差加大，弱透水层有效应力增加，必定产生固结变形 。土体固结变形又可能会导致 ［6-7］ 孔隙度减小，渗透性变差和单位储水系数减小，进而改变弱透水层中孔隙水的渗流特性 。对于咸 淡水界面运移而言，弱透水层中的非达西渗流过程直接影响溶质迁移的速度，而弱透水层释水，也 ［8］ 可能会导致阻滞因子及弥散系数的改变，进而影响溶质迁移过程 。可见，由地下水开采引起的咸水 下移问题属于典型的渗流-固结-溶质运移耦合问题。 ［9-10］ 目前，变形多孔介质溶质运移的研究主要集中于土壤环境污染领域，如 Smith、Lewis等 对垃 圾填埋场渗滤液的污染问题进行了研究，张志红、李涛等［11-12］研究了湖泊底泥堆场中污染物迁移问 ［13-16］ 题。耦合模型中黏性土的渗流多假设符合达西定律，而实际上其通常表现出低速非达西渗流特性 。 当前，对非达西渗流规律的刻画多采用Hansbo公式［17-20］。本文以华北平原典型的“上咸下淡”含水层 系统为研究对象，采用一维非达西渗流-固结-溶质运移模型，在忽略固结过程导致弱透水层水文地 质参数变化的前提下，重点考虑弱透水层中非达西渗流过程导致的渗透流速分布改变，讨论其对咸 淡水界面下移过程的影响。 收稿日期：2013-11-15；网络出版日期：2014-05-12 网络出版地址：http：///kcms/detail/11.1882.TV1405.003.html 基金项目：国家重点基础研究发展计划973项目（20