基于modflow的地下水硝酸盐运移模型.docxVIP

基于modflow的地下水硝酸盐运移模型 地下水通常受亚硝酸盐（no-3）污染（康熙等人，1993；朱ang等人，1996；诺兰等人，1997；张云等人，2006）。影响地下水NO-3浓度的因素复杂,而数值模型提供了地下水管理的有力工具,因为这种模型综合了气象、水文、地质、水文地质和土地利用等多方面数据,直接或间接地建立了相互间的联系,使得人们能够定量地分析和预测含水层中NO-3浓度的变化。但是分析和预测的精度直接受边界条件、含水层参数和源汇项等有关数据精度的影响。为了有效地校正这些估计值,近年来国际上采用的方法是延长模拟时间(Delwyn, 1997;Liu et al.,2004;Nishikawa et al.,2003),通常是从地下水开发之初到最近某一年。这样,估计值的精度可以通过长期的累积效应表现出来,而估计值的校正可以通过长期的模拟值和实测值的拟合获得。 石家庄市以地下水作为主要供水水源。随着人类活动的加剧,地下水位持续下降,水质不断恶化,尤其是NO-3浓度不断升高,从1959年的4.38 mg/l上升到2000年的均值54.07 mg/l,最大值达126.90 mg/l。出于不同的目的,该区已建立多个水流模型(林学钰,1987;许广明等,2001),而NO-3运移模型少见报导。本次研究的目的是评价各种管理方案对石家庄市地下水NO-3浓度的影响。采用的方法是利用Visual Modflow软件(Waterloo Hydrogeologic Inc.,2000)进行地下水水流和硝酸盐运移模拟,建立长达42 a的数值模型。 1 地下水上、副三大区域 石家庄市位于河北省中南部,西依太行山,东接河北低平原。研究区范围介于东经114°16′～114°41′,北纬 37°56′～38°17′,包括市区、郊区和周围县市部分地区,面积约1500 km2(图1)。地下水主要赋存在第四系松散沉积物中,其中上部中、上更新统和全新统沉积物主要由卵砾石和砂砾石组成,构成统一的非承压含水层。地下水流场概化为非均值各向同性二维非稳定流。 1.1 模拟时间和网格大小 水流模型和NO-3运移模型采用相同的有限差分网格,在水平方向上按210行185列将研究区剖分成38850个网格, 每个网格大小约为200 m×200 m, 面积40000 m2。模拟时间为1959-2000年,共42 a。根据获得的水位和硝酸盐浓度观测数据,模拟时间分两部分,第一部分是1959-1995年,以1 a为一个应力期,每个应力期分10个时间步长。第二部分是1996-2000年,以0.5 a为1个应力期,每个应力期也分10个时间步长。 1.2 坝体边界质区 模拟区范围是: 西边界以山区与平原交界线为边界;西北边界以黄壁庄水库坝体渗漏段为边界,北边界以滹沱河北岸的一条流线为边界,东南边界以地下水水位下降漏斗分水岭线为边界。模拟区面积约895 km2。所有侧向边界均为二类边界,其中北边界为0流量边界,西北水库坝体边界的渗漏量为(0.75 ～1.60)×108m3/a。 西边界为透水—不透水边界。东南边界的性质随水位下降漏斗的发展而变化,1980年之前为流出边界, 之后为0流量边界。西边界和东南流出边界流量通过收集以往调查成果或根据历史时期的等水位线图计算得到。 底边界以潜水含水层的底板为边界,西部以基岩为底板,中、东部为一厚度大于15 m的粘土层,均为隔水边界。 顶边界为潜水面,受补排量的影响,为一随时间而变化的界面。 所有流入边界均赋浓度值,其中西北边界坝下渗漏水浓度在1990年以前为背景浓度4.43 mg/l,以后为水库水现测平均浓度9.31 mg/l,西边界流入水浓度也为背景浓度。 1.3 no-3具有生物生态学意义 水流模型参数主要包括渗透系数K和给水度μ。溶质运移参数包括迟滞因子Rf、衰变系数λ和水动力弥散系数D。本次模拟对象是NO-3,它是一种易迁移的化学组分,不被孔隙介质的固相吸附,因此,设定Rf等于1。使NO-3浓度减少的生物化学作用是反硝化作用。然而,含水层为潜水含水层,岩性颗粒粗,溶解氧DO浓度高,均不利于厌氧反硝化作用的发生(张翠云等,2005),因此,假定NO-3没有衰减, λ等于0。水动力弥散系数D是描述由机械弥散和分子扩散引起溶质迁移的因子。分子扩散较机械弥散小得多,因此忽略分子扩散。根据地形、地貌和岩性,对渗透系数、给水度、纵向弥散度、孔隙度、入渗系数等参数进行分区。参数初值估计来源于前人成果(林学钰,1987)以及国外总结的不同岩性参考值(Gelhar et al., 1992;Zheng et al., 2002)。 1.4 污水的补给源 模型补给包括大气降水入渗、河渠渗漏、农田灌溉水入渗和污水入渗。模型排泄包括农业地下水开采和城市地下水开采