含内源污染平面二维水流-水质耦合模型 2-D plane flow-water quality-bed sediment pollutant coupled model.pdfVIP

水 利 学 报 2006年2月 SHUIUXUEBAO 第37卷第2期 文章编号：0559．9350(2006)02．0205．05 含内源污染平面二维水流一水质耦合模型 龚春生1，姚琪1，范成新2，赵棣华3 (1．河海大学环境科学与环境工程学院，江苏南京210098；2．中国科学院南京地理与湖泊研究所，江苏南京210008； 3．南京水利科学研究院，江苏南京210024) 摘要：本文从紊动扩散理论出发根据浅水湖泊紊流扰动特性，推导出包含底泥污染的水质控制方程，并与浅水方程 耦合，建立相应的平面二维水流一水质一底泥污染的数值模型。模型在有限体积法框架下应用黎曼近似解求得耦 合方程的数值解。并以玄武湖为例进行模型验证。结果表明，数值解与实测值较吻合，证明预测结果可信。 关键词：底泥污染；有限体积法；黎曼近似解；对流扩散方程 中图分类号：TVl31．2 文献标识码：A 国内外学者对浅水湖泊的水流一水质模型作了大量研究，但大多数研究没能将内源污染与水质模 型综合起来¨以3。本文从紊动扩散理论的观点出发、根据湖底处沉积物一水界面质量守恒建立了浅水 湖泊平面二维水流一水质一底泥污染耦合模型。在有限体积法框架下应用黎曼近似解格式解算模型中 各跨单元边界的数值通量，进而求得耦合方程的数值解。根据玄武湖调水实验的实测资料对模型进行 了验证，数据检验表明，结果合理，模型可靠。 1 包含底泥污染的水质方程的建立 根据紊动扩散理论H“。，包含内源污染的三维对流一扩散方程为 瓦3C+未(cⅡ)+未(踟)+未(c埘)+毫(c5二) (1) 间隙水浓度；二为底泥在外力扰动作用下的综合再悬浮速度；D，、D，、D：分别为污染物在水体中的扩散 系数；D：为沉积物一水界面的扩散系数；d为泥一水界面的厚度；K为污染物在水体中的综合降解系 数；S。为源、汇项。 对式(1)沿水深^进行积分 d彳 ㈦ 』：昌(见蓦)¨f芝㈦誓卜+r未(D：筹卜一．f。h删z+卜h 由于水浅，湖泊的水流和水质在垂直方向混合比较充分，因此可假设：(1)水力变量及污染物浓度沿 水深不变．乖向流速为0：(2)通过紊动扩散及污染物再悬浮产生的底泥污染发生在泥一水界面后迅速 收稿日期：2004—10—15 基金项目：国家自然科学基金重点资助项目；中国科学院知识创新工程项目(KZCX3一SW．348) -·-——205·-——— 万方数据 在垂向混合。因此，包含底泥污染的水平二维水质基本方程为 a(hC)．a(hCu)．a(hCv)，一。。+ o一 ——犷+——百焉—一+——否。i—一+L叫I =未(D，掣)+苦(珥掣)+D，等呼一勋c+s 式中：D：嘉I：：为沉积物一水界面物质交换通量；Cw～l：：表示底泥表层污染物颗粒再悬浮引起的通 量；S为垂向平均后的源、汇项。 1．1沉积物一水界面物质交换通量的确定沉积物一水界面物质交换主要由表层沉积物的分解、底泥 间隙水的扩散引起，考虑到上覆水与间隙水的浓差较大，再加上厌氧条件下表层沉积物的分解速率较 慢，沉积物一水界面物质交换主要归因于底泥间隙水的扩散。因此扩散通量D：嘉l：：=一D：