地下水溶质运移数值模型.docx

《地下水溶质运移数值模型.docx》由会员分享，可在线阅读，更多相关《地下水溶质运移数值模型.docx（3页珍藏版）》请在第一文库网上搜索。

1、地下水溶质运移数值模型(资料性附录)水是溶质运移的载体，地下水溶质运移数值模拟宜在地下水流场模拟基础上，因此地下水溶质运移数值模型包括水流模型和溶质运移模型两部分。DJ地下水水流模型非均质、各向异性、空间三维结构、非稳定地下水流系统:1)控制方程hhy3(“h(hS,=Kv+Ky+K+/t3xIx)为(yJAzIz)式中：SS一一给水度I/；h水位1；KfKy,Kz分别为X,y,Z方向上的渗透系数E;T一一时间T；Qs一一源汇项m注：方括号中的符号为量纲，以下同。2)初始条件h(xyy9zyt)=Zz0(x,y,z)(x,y,z),/=O式中：4*，y，z)已知水位分布：Q模型模拟区。3)边界

2、条件：第一类边界：(x,y,z)=(x,y,z)(x,y,z)1,rO式中：r一一类边界；h(x,y,z,t)一一类边界上的己知水位函数。第二类边界：nq(x,y,Z)(x,y,z)2式中：2二类边界;三维空间上的渗透系数张量;nn边界r2的外法线方向；q(x,ytz)二类边界上已知流量函数。第三类边界：r(k(h-z)-+ah)=q(x,y,z)加r3式中：0一一系数；3一二类边界；k一一三维空间上的渗透系数张量；n边界G的外法线方向;q(x,yfz)三类边界上已知流量函数。D.2地下水水质模型1)控制方程R啜喘吗(他C)Fe/元式中：R迟滞系数，无量纲PbSCCph介质密度IM1-3；介质

3、孔隙度，无量纲；C组分的浓度M1，；亍一一介质骨架吸附的溶质浓度M1，；t时间T；X,y,Z一空间位置坐标1；Dij水动力弥散系数张量1?T；Vi地下水渗流速度张量E;qs源和汇T;CJ一一源或汇水流中组分的浓度M1；4一一溶解相一级反应速率T；4一一吸附相一级反应速率Tj。2)初始条件C(x,y,z)=C0(x,y,z)(x,y,z)=O式中：C0*,y,z)一已知浓度分布；Q模型模拟区域。3)定解条件第一类边界Dirich1et边界：C(x,y9ZJ)=Co(X，Xz)(x,y,z)1,rO式中：r一定浓度边界；C(X,y,z,f)一定浓度边界上的浓度分布。第二类边界Neumann边界：dCODij=y,zj)(x,y,z)2,OOXj式中：2一通量边界；E(X,y,zj)一一边界2上已知的弥散通量函数。第三类边界Cauchy边界：ecODii-qiC=gi(x,y,z)(x,y,z)3,OOXj式中：3一一混合边界；g,(,y,z,z)一一边界3上己知的对流弥散总的通量函数。