我国西北盆地地下水流建模的若干问题——以河西走廊黑河干流中游盆地为例研究.pdf

我国西北盆地地下水流建模的若干问题 ——以河西走廊黑河干流中游盆地为例 陈崇希’，胡立堂2 摘要：地下水数值模拟是地下水资源评价的有效方法，它涉及的是井孔一含水系统的问题。当今在建立地下水流的 概念模型、数学模型上往往忽视了对井孔的描述和井流特征的分析，或者缺乏符合水流机理的模拟技术。本文以河 西走廊黑河干流中游盆地地下水流问题为例，分析了建立三维流基本模型的必要性，重点剖析了目前国际上应用最 广的MODFLOW软件在模拟混合抽水井的不妥，也涉及到了混合观测孔和泉等的模拟问题以及河流、渠系、降雨滞 后补给潜水的刻画问题，这些问题在提高我国西北盆地地下水流概念模型和数值模型的仿真性上十分重要。 关键词：地下水混合抽水井混合观测孔泉滞后补给数值模拟渗流一管流耦合模型 黑河 地下水流数值模拟是地下水资源评价和管理的有效方法，它涉及的是井孔一含水系统的问题…。 然而当今在建立地下水流的概念模型、数学模型上往往忽视了对井孔的描述和井流特征的分析，尽 管研究区存在混合抽水井、混合观测孔等，但在建模中却往往没有考虑这些源汇项或边界条件的刻 画；泉流量动态的模拟和河流、渠系、降雨滞后补给潜水的刻画也是建模过程中易被忽视或刻画不 佳的问题，以上这些问题都会导致模拟失真。本文针对这些当前较普遍的问题，以河西走廊黑河干 流中游盆地为例，作适当的分析及相应的改进。 1水文地质概况 研究区为南东一北西走向的盆地(张液盆地及酒泉东盆地)，长度约200km，最大宽度处约80kin， 面积约8700 km2。南、北、东侧为山区，西侧以地下水分水岭为界；地形：总体上为南高北低，东 高西低。该区为典型的干旱、半干旱大陆性气候特征，多年平均气温为7C，多年平均降水量约为 129 mm，多年平均蒸发量为2048mm。区内大小河流都发源于祁连山，从祁连山进入到黑河中 过张液市后转向酲北，通过正义峡流出本盆地，进入额济纳盆地(内陆盆地)。 盆地内以大型洪积扇为骨架，形成陡倾斜砾石平原和缓倾斜细土平原，前者主要分布于区内南 北两侧山前，后者分布于区内中部地带。研究对象第四系孔隙含水系统的厚度以盆地中部为最大， 研究区南部洪积扇上部为比较均匀的砂卵石、砂砾石含水层，向扇的前缘及黑河冲积平原逐渐转变 为多层含水层系统。 中的数值为甘肃二水1984年的研究成果(下同)。 作者简介：陈崇希，中国地质大学教授、博导，主要从事地下水流数值模拟研究。E—mail：—ccx33@1—26．corn 936 (3．60亿m3／年)的排泄和人工开采地下水(1．08亿m3／年)。 图l河西走廊黑河干流中游盆地白然地理图 LocationofthemiddlereachesofHeiheRiverBasinatHexi Fig．1 map Corridor 2概念模型与模拟技术中的几个问题 2．1基本模型 模型设计第一个问题是采用什么基本模型，是平面二维、准三维还是三维模型，它们在资料上、 运算分析工作量上及模拟仿真性上都有很大的差异。基于本区的地下水流动的垂直分量在许多因 素作用下不可忽略不计，甚至在某些条件下是主导分量，例如：洪积扇顶部出山口河水渗漏补给后，地下 水将以垂向为主地向下游流动；洪积扇前缘地下水向泉口溢出也是向上垂直流动；含水层之间弱透水 层中的地下水也以垂向流动为主；黑河排泄段(张液黑河大桥以下)非完整河底部的地下水也以垂直向 上流动为主；非完整抽水井附近的地下水流存在垂直流速分量，等等。因此，本区地下水流属于三维流 类型。那么采用准三维流模型还是三维流模型呢?考虑到具体条件和建模目标一区域地下水资源评价 与管理，实际含水系统的层数很多而模型分层又不宜过多，因此概化后的“弱透水层”中含有一定数量 的砂层，即“弱透水层”是各向异性的。陈崇希等(2001)旺耶1在研究苏州市地面沉降的模型时指出：“准 三维流模型‘误差小于5％’的结论仅适用于弱透水层为各向同性介质；对于弱透水层的各向异性比为 10时，误差已超过27．7％，而且随着弱透水层单位储水系数∥：的减小和模拟时间的延长，误差还要 增大，已存在超过30．5％的情况。基于此研究结果，本专题对苏州市采用三维不稳定流动模型”。薛禹 群(2003)【4l也提出相似的定