MODFLOW在磁县地下水数值模拟中的应用.pdfVIP

第33卷第 2期 华 北 水 利 水 电 学 院 学 报 VoL33No．2 2012年 4月 JournalofNoAhChinaInstituteofWaterConservancyandHydroelectricPower Apr．2012 文章编号 ：1002—5634(2012)02—0053—03 MODFLOW 在磁县地下水数值模拟中的应用 李树文 ，于虎广 (1．河北工程大学，河北 邯郸 056038；2．河北省资源勘测研究重点实验室，河北 邯郸 056038) 摘 要 ：磁县东部平原区水环境 由于水资源利用不合理、地下水超采等原因发生恶化．通过分析研究区水文 地质条件 ，运用 MODFLOW进行浅层地下水数值模拟．结果表明，在保持开采条件不变 的情况下，至2020年 ， 漏斗面积 由l8．59km 增大为 23．31km ．提出合理调整和优化地下水开采方案可以有效地避免降落漏斗影 响范围的增大． 关键词 ：地下水 ；数值模拟 ；MODFLOW 研究区位于磁县东部平原 ，属暖温带半湿润半干 受补给和排泄，补给方式主要为地表水体入渗，排泄 旱大陆性季风气候区，多年平均降水量 531。4mill，每 方式为垂直蒸发 ，故将上部平面边界概化为二类流 年汛期 (6—9月份)的降水量约 占全年降水量的 量边界．下部边界主要为几乎不透水的泥岩 ，故概化 70％ 一80％ ． 为隔水边界． 研究区是 由漳河 、滏 阳河和!}亡牛河 的山前冲积 综上所述 ，研究区可以概化为非均质各向同性 物组成的太行山前冲积扇．其中南部为漳河冲积扇 ， 三维非稳定地下水流系统． 北部是沙沼河和漳河冲积、洪积扇之间的交接带，多 1．2 数学模拟模型 是滏阳河和．牛亡牛河 出山口的冲积和洪积物． 数学模型是在研究区水文地质条件的基础上建 地下水赋存于第 四系松散岩层 中，主要受大气 立的，如下式所示 ： 降水和地表水体渗入补给及山前潜水 的侧 向补给． 毫( +( +去( +W=SsOH， 含水砂层分布显示 了典型的山前冲积扇规律 ，主流 部位，含水层砂粒粗、厚度大 ，向主流两翼渐变为 中 It(x，Y，=，￡)Is,= (，)，，，f)， (，Y，)∈s，￡≥0， 细砂 ，地下水流向与地形基本一致． KO _ a_ (，y,z)∈s2，t>I0， 近年来该地 区地下水供水量约占区域总供水量 nIs=g(，，，，z，)， 的70％左右 ，长期超量开采致使地 下水资源衰 n(x，Y，z，￡)l；。=Ho(x，Y，)， (，y，z)∈力， 竭 ，造成地面沉降，当地的供水安全受到严重影响． 式中：K为含水层渗透系数 ，m／d； 为单位时间、单 位体积上从垂直方 向上流人或流 出含水层 的水量 1 三维地下水流模拟 (流人为正，流出为负)，In／d；S 为贮水率；力为渗 1．1 水文地质概念模型 流区域；H为水位 ，in；t为时间，d；Ho为初始时刻水 研究区西部为山区，接受山前侧 向径流补给 ，将 位值 ，In． 研究 区西部边界概化为二类流量边界 ，流量 以注水 模型采用 VisualMODFLOW 4．2进行求解．利 井的形式赋入模型，各段流量值利用达西定律分别 用 MapGIS对研究区 1：50000第 四系水文地质图做 计算．