河流影响下包气带水力参数对压力水头及溶质运移的影响

河流影响下包气带水力参数对压力水头及溶质运移的影响刘钢;高志鹏;屈吉鸿【摘要】地表河流入渗补给地下水是自然界普遍存在的一种现象.近些年，由于地表河流污染的加剧污染的河流下渗,必然伴随着污染物对影响范围内的非饱和带土壤造成生态破坏.采用Hydrus软件构建二维饱和-非饱和数值模型，为各影响因子构造3种不同的情景模式，探讨河流影响下的包气带水分及溶质运移随各参数变化而改变的过程.研究结果表明,VG模型5个参数er、es、a、n、Ks对压力水头和溶质运移均具有不同程度的影响，与土壤饱和程度具有一定的相关性，其中es、a、n、Ks的影响较为显著，而0r的影响则比较小.％Groundwaterrechargedbytheriverongroundsurfaceisacommonphenomenoninthenaturalworld.Inrecentyears,duetotheintensificationofriverpollution,theinfiltrationofpollutedrivermustbeaccompaniedbypollutantstransporttotheunsaturatedzoneandcauseddamagetotheecology.Inthisstudy,atwodimensionalunsaturated-saturatednumericalmodelbasedonHydrussoftwarewasusedtoinvestigatethechangesofwaterandsolutetransportwithparametervariationintheunsaturatedzoneundertheinfluenceofriver,andthreedifferentscenarioswerebuiltforeveryinfluencefactor.Theresultsofthestudyshowthat:ThefiveparametersofOr,es,a,n,KsinVGmodelhaveadifferentinfluenceonthepressureheadandsolutetransport,andthefiveparametershaveacertaincorrelationwiththedegreeofsaturation,theeffectofOs,Ks,a,nissignificant,andtheeffectof0risrelativelysmall.【期刊名称】《华北水利水电学院学报》【年(卷)，期】2017(038)002【总页数】5页(P72-76)【关键词】河流;包气带;水力参数;Hydrus;压力水头;溶质运移【作者】刘钢;高志鹏;屈吉鸿【作者单位】黄河勘测规划设计有限公司，河南郑州450003;华北水利水电大学资源与环境学院，河南郑州450045;华北水利水电大学资源与环境学院，河南郑州I450045【正文语种】中文【中图分类】TV131.2+9;S152.7+2自20世纪50年代以来，自然因素的激烈变化和人类活动的加剧，对河流沿岸地下水资源的形成与水质和生态环境的演化产生了重要影响。包气带是河流入渗补给地下水的必经之路，包气带的水力参数控制着水流的入渗，对溶质的迁移转化也具有重要的影响［1］。如今，越来越多的学者意识到包气带对地下水水质的重要性，开始研究饱和-非饱和带相耦合的数值模型［2-3］。描述包气带水力特性的模型很多，以vanGenuchten(VG)模型应用最为广泛，该模型可以描述土壤的水分特征曲线，还可以定量求得非饱和导水率［4］。VG模型包括5个参数，对非饱和带的水流及溶质运移具有很大的影响。但究竟其影响的规律如何，影响的机理又如何，目前这方面的研究还比较少。李世锋等［5］采用Monte-Carlo方法模拟了非饱和带渗透系数的空间不确定性对溶质运移的影响，结果表明其不确定性使浓度呈现明显的不确定性；范严伟［6］采用Hydrus-1D模型研究了VG模型参数变化对土壤入渗特性的影响，指出在推求VG模型参数时，应保证0s.Ks、n的准确性；王志涛［7］采用不同的思路研究了土壤垂直一维入渗对VG模型参数的敏感性，与范严伟的研究得到了相同的结果。本文以VG模型为理论基础，采用Hydrus软件构建二维饱和-非饱和流数值模型，通过构建不同的情景模式，研究某断面上压力水头和铵态氮浓度随深度的变化规律，探讨河流影响下的包气带压力水头及溶质运移随各参数变化的改变过程。1.1构建模型采用Hydrus软件构建理想模型，如图1所示。首先建立基准模型(Std):模型长5m、高2m,采用四边形剖分原则，共剖分132个节点，262个单元格；采用伽辽金有限单元进行空间离散，模拟水分及溶质运移，模拟期为200d;以软件中壤土的VG模型参数经验值为基准模型参数，模拟溶质为铵态氮，铵态氮初始浓度为1mg/L,纵向弥散度(DL)为0.5cm，横向弥散度(DT)为0.1cm，溶质分子扩散系数(D)为4.32cm2/d，线性吸附分配系数为0.002cm3/mg，一级硝化速率为0.005d-1。模型剖分及边界条件如图1所示，河流边界(AB)为定水头边界，河流宽50cm，深10cm,BC为大气边界(因本模型未考虑降雨和蒸发的作用，因此BC边界可简化为隔水边界)、下边界(DE)为自由排水边界，左右边界(AE、CD)为零通量边界；溶质运移河流边界为定浓度边界，左右边界为零通量边界，下边界为零浓度梯度边界。模拟结果以压力水头及铵态氮浓度在AE边界上随深度的变化给出。1.2包气带水分运移模型包气带的水流状态是变饱和流。描述饱和-非饱和流水流运动方程为经典的Richards方程，选取河流影响下的某剖面为研究对象，概化为二维饱和-非饱和流方程:其中：式中：h为压力水头，cm;C(h)为容水度，cm-1;0为土壤体积含水率；K(h)为非饱和导水率，cm/d；h0(x,z)为初始压力水头分布，cm；d0为河流边界上压力水头；x为水平方向空间坐标变量，cm；z为垂直方向空间坐标变量，cm，向上为正；t为时间变量，d；kr为相对渗透系数；Ks为饱和导水率；Os为饱和含水率；由为残余含水率；ee为有效含水率；a、n为拟合参数；h<0表示包气带，h>0表示饱和带。式(3)—(7)为vanGenuchten描述非饱和土壤的经验公式［8］。1.3包气带溶质运移模型铵态氮迁移转化模型的建立仅考虑了主要的迁移转化过程，包括对流、弥散、吸附、硝化等，忽略作物根系吸收作用、矿化作用以及解吸等，采用传统的对流-弥散方程(Convection-DispersionEquation,CDE)来描述污染物运移过程［9］。包气带溶质运移模型为：式中：c为铵态氮的浓度，mg/L；Dxx、Dzz、Dxz、Dzx为水动力弥散系数张量D的分量，cm2/d；qx、qz分别为x、z方向上达西流速分量，cm/d；K1为一级硝化反应速率，d-1；c0为初始浓度分布，mg/L。1.4情景设置本文仅探讨VG模型参数的扰动对包气带压力水头及溶质运移的影响，因此在基准情景的基础上对VG模型的5个参数做出一定的扰动，而不改变其他参数(如弥散度、分配系数、硝化速率等)的取值。由式(3)—(7)可知，VG模型的主要参数有Or、Os、a、n、Ks，选取壤土的经验参数值为基准模型参数，各参数均分2种情景模式进行上下扰动，各情景模式设置见表1。2.1扰动指数定义浓度扰动指数为：式中：c1为某情景下铵态氮的浓度;c0为基准情景下的铵态氮浓度；△为各参数的扰动值，△为正值表示正扰动，△为负值表示负扰动。某参数的浓度扰动指数大于100%(或小于-100%),说明该参数的扰动对溶质浓度的影响较为显著；浓度扰动指数介于100%与-100%之间，则说明该参数的扰动对溶质浓度的影响较小。将某一参数的扰动指数求平均值，平均值为正，表明该参数的扰动与模型输出变量(压力水头或溶质浓度)呈正相关；平均值为负，表明该参数的扰动与模型输出变量呈负相关。2.2各参数变化对压力水头和溶质运移的影响将各情景模式代入数值模型中，分别求解，得到不同时期压力水头及铵态氮浓度扰动指数随深度变化的曲线。当t=1d时，VG模型中5个参数er、es、a、n、Ks的扰动引起压力水头和溶质运移的变化情况如图2所示。A1和A2情景分别为0r扰动-10%和10%，从图2(a)可以看出，t=1d时，压力水头随深度变化的曲线基本无变化，溶质浓度的扰动指数值在模拟深度范围内的最大值为80%，小于100%，说明Or对压力水头和溶质运移的影响很小。同时，其平均浓度扰动指数为正，说明溶质浓度随Or的增大而增大。B1和B2情景分别对Os扰动-8%和10%，从图2(b)可以看出在B1情景下的压力水头较基准情景下的增大，B2情景下的压力水头则减小，说明Os对压力水头的影响为负相关；Os对溶质浓度的扰动指数均小于-100%，说明其影响很大，且其平均值为负，说明Os增大将伴随溶质浓度的减小。从图2(a)—2(e)可以看出，a与Os具有相似的规律，但其影响程度与Os的影响程度相比较小；n、Ks与Os、a的影响规律相反，与Or相似，参数值增大将引起压力水头及溶质浓度的增大。Ks对浓度的影响程度最大(图2(e)),浓度的最大扰动指数超过了1000%；而n的影响相对较小(图2(d)),浓度的最大扰动指数小于200%。由图2还可以看出，在0~50cm范围内形成了饱和带，50~200cm为非饱和带，说明在河床下形成了悬挂饱和带-非饱和带系统［10］。当t=200d时，5个参数的变化对压力水头和溶质运移的影响情况如图3所示。从图3中可以看出，t=200d时，研究深度内的包气带最小压力水头仅为-6cm，其饱和程度已非常高。与图2相比，t=1d时5个参数的扰动对压力水头和溶质运移的影响具有不同的规律。对于压力水头，5个参数的扰动使其随深度的变化呈现出一定的波动性，在深度0~100cm范围与100-200cm范围内的规律相反。整体上，因饱和程度的提高，t=200d时，5个参数对溶质运移的影响程度均有所下降。由图3(a)和3(b)可以看出，由和伽对压力水头随深度的变化影响较为明显。对溶质运移来说，其扰动指数大小均介于100%与-100%之间，说明Or和0s的扰动对溶质运移的影响较小°0r对浓度的最大扰动指数仅为-1.5%;0s对浓度的最大扰动从t=1d时的-591%变为t=200d时的-28.6%。由图3(c)、3(d)、3(e)可以看出，压力水头随深度的变化明显，而溶质浓度的扰动指数最大值均超过了100%，说明a、n、Ks的扰动对压力水头和溶质运移的影响依然较大。但与t=1d时相比，其影响程度均有所下降。综上所述，0r、0s、a、n、Ks对压力水头和溶质运移均具有不同程度的影响，其影响规律随研究区域饱和程度的不同而不同。t=1d时，饱和程度低，VG模型5个参数对压力水头和溶质运移的影响呈现出不同的规律。其中0r、n、Ks的影响规律相似，其对压力水头和溶质运移的影响呈正相关，即0r、n、Ks增大，伴随着压力水头和溶质浓度相应增大；相反的，0s、a对压力水头和溶质运移的影响则呈负相关。整体上，0r对压力水头和溶质浓度的影响较小，0s、a、n、Ks的影响较大。t=200d时，饱和程度很高，0r和0s的扰动对压力水头的影响较大，对溶质运移的影响则较小，扰动指数均小于100%；而a、n、Ks的扰动对压力水头和溶质运移的影响均较大。与t=1d时相比，t=200d时研究范围内土壤饱和度高，描述非饱和带水力特性的VG模型的5个参数的扰动对压力水头和溶质运移的影响程度则有所降低。采用Hydrus构建理想模型，通过模拟河流影响下的非饱和土壤中水分及铵态氮运移的规律，探讨了描述非饱和土壤水力特性的VG模型各参数对水分及溶质运移的影响程度。研究结果表明，VG模型5个参数对压力水头和溶质运移均具有不同程度的影响，与饱和程度具有一定的相关性，其中0s.Ks、a、n的影响较为显著，而0r的影响则比较小。【相关文献】梁斌，王超，王沛芳，等.河流污水非饱和入渗对沿岸地下水质的影响[J].水科学进展，2003,14(5):548-553.赵贵章.鄂尔多斯盆地风沙滩地区包气带水-地下水转化机理研究[D].西安：长安大学，2011.林琳，杨金忠，史良胜，等.区域饱和-非饱和多孔介质的溶质运移简化模型[J].水利学报，2007,38(3)：342-348.LIUYi,KUANGXingxing,JIAOJiuJimmy,etal.Numericalstudyofvariable-densityflowandunsaturated-saturatedporousmedia[J].Journal