在井流影响下对流—弥散问题研究.pdfVIP

维普资讯 水 利 学 报 1995年 5月 $HUiLI XUEBAO 第5期 在井流影响下的对流一弥散 问题研究 、J薛禹群‘ ． 张志 吴 吉春 ＼ 1 南京 种 z 提 要 。 本文针对在对流一弥散同题的数值计算 中，井点附近容易出现不合理的现象， 理论上推导 出井流项(源汇项)在对流一弥散方程中的具体表达式，探讨 了问题出现的原因，并提出了求解对流 散 问题的4种计算方案 最后通过一个算铡，验证 rl本文处理源汇项的可行性 ，并对 4个计 算方案进行 了比较．计算结果表 明，在实际应用中，可 以忽略弹性储(释)水对溶质运移的影响， 也可把注水井作为浓度第一类边界处理，经这样处理后所引入的误差不超过o．5％． 糊 踯 ’ 蛳 水 ． 耐l浇 ． 对哦 撇 ．． 哪 一 、 前 言 井流对污染物和热量的运移有重要作用，但在许多文献中对考虑开采井或注水井影 响的描述溶质运移的对流一弥散方程一般写为 … 考(Dc善) 十 、‘ ㈩ 式中 为单位体积含水层的抽(注)水强度(抽水为负，注水为正)， 为有效孔隙度，而 rc，为注入水的浓度(W0) c． 为抽出水的浓度 ( 0) 当利用这一方程计算时，常发现井点及其附近浓度计算值会出现令人无法相信的情 况．脓明佑 曾以从盛糖水的缸里取水为例把源汇项修改为 一c’)／ ，但没有给 出严格的推导，或明确指出上述问题的真正所在，也没能用算例验证，不能令人信服． 实际上，G．Galeati和 G．Gambolati 曾经得出类似的结果，但同时又对其正确性 表示怀疑，为此他们建议把源汇项作为适当的外边界条件处理．其实，上式只有在满足 0的条件下才成立，即_只有在稳定流且水流区内投有井等源汇项( =0)的条件 下，或稳定流且井作为边界条件处理时才是正确的，即把它应用于非稳定流或水流区有 井的稳定流情况就不正确了，计算结果必然出现错误． 本文从她下水水流方程、溶质运移方程及达西定律出发，从数学上严密推导出源汇 项在对流一 蔽秀褪币丽真诔毳i； 并给出存在井流影响时溶质运移问题的4种计算 +本文于1992年loft13日收到．丰研究得到国家 自然科学基金的赞助 一 6。一 维普资讯 方案及其有限差解．最后通过一个实例计算，对几种计算方案进行了比较 二、存在井流时的对流一弥散方程 对流一弥散方程 中考虑溶质随水流一起运移的对流项有形式a(V，c)／ ．因此对流 一 弥散方程的一般形式为 “ (D )一 c) ac一 (1) ax ax ax ‘ aI ‘ 当有抽(注)水井等垂 向水量交换时，根据质量守恒原理