第 30 卷第 5 期 西 南 大 学 学 报 (自然科学版) 2008 年 5 月
Vol畅30 No畅5 Journal of Southwest University (Natural Science Edition) May 2008
文章编号: 1673 9868(2008)05 0145 04
MATLAB 软 件 在 构 建 潜 流 人 工 湿 地
水 力 学 模 型 中 的 应 用
①
戚景南, 黄玉明
西南大学 化学化工学院, 重庆 400715; 三峡库区生态环境教育部重点实验室, 重庆 400715
摘要:根据示踪试验的数据, 由考虑扩散过程的推流模型计算得到的平均水力停留时间 τ和水力停留时间分布的标
准方差 σ
2
θ
的理论值; 利用数学软件 MATLAB7畅1 模拟潜流人工湿地中的示踪剂浓度———时间响应曲线获得水力停
留时间分布(RTD)函数, 根据 MATLAB 软件中不同函数拟合得到的 τ值与理论值间的比较, 发现对数正态概率密
度函数可以较好地拟合实际得到的示踪试验浓度———时间曲线, 本研究的潜流人工湿地中水流状况为有扩散的推
流模式.
关 键 词: 潜流人工湿地; MATLAB 软件; 水力停留时间分布函数; 示踪试验; 模型
中图分类号: X703 文献标识码: A
人工湿地对废水中污染物有良好的去除作用, 已广泛地用于废水处理中
[1 3]
; 很多水流状况通常用
模型来表示
[4 ,5]
, 目前通常采用理想推流模型来描述湿地内部的水流状况
[6]
, 在该模型中假设示踪剂到
达出口的瞬间, 其浓度达到最大, 没有横向扩散和混合; 也有采用完全混和流模型来模拟湿地中水流状
况的研究
[7]
, 在该模型中, 假设示踪剂进入湿地的瞬间即混合均匀, 而且在流向出口的过程中其浓度保
持不变.由于湿地内部结构及植物根系的复杂性、 外部降水和蒸发作用等因素对水流状态有影响
[8,9]
,
加上“短流”现象的存在, 使得实际得到的示踪试验浓度———时间曲线中, 示踪剂浓度达到峰值之后会有
一个“拖尾”, 因此, 纯粹理想推流或者完全混合流模型都不能全面地模拟实际得到的不对称的示踪剂
浓度———时间曲线.
Kadlec
[10]
提出在湿地中利用一系列平行的网络结构模型来模拟水流状况.Coates 等
[11]
在他们提出的
“微分电容”模型中把平流扩散方程和完全混合存储区域结合起来, 以便更好地模拟多孔介质中的水流流
动.Werner 等
[12]
提出了不完全混合区域模型(ZDM)来描述湿地中的水流状况, ZDM 模型包含了无数多个
推流反应器, 每个都和紧邻的 CSTRs 进行质量交换, 为了使模型更加准确地模拟实际的水流状况, 他们在
垂直于流体运动方向的基础上考虑径向扩散, 得到的模型非常类似于“微分电容”模型和“二元区域”模型.
综合以上分析可以看出必须同时考虑理想推流与完全混合流才能更加准确地模拟实际的水流状况.
本研究利用 MATLAB 软件中的几种函数对示踪试验中得到的数据进行拟合分析, 以获得湿地中的水力
停留时间分布函数, 进一步得到平均水力停留时间 τ及水力停留时间分布的标准方差 σ
2
θ
值, 将几种函数拟
合获得的 τ值与理论值进行比较, 发现用对数正态概率密度函数进行拟合能够得到满意结果; 本研究提出
的方法也有可能用于描述其他水处理单元的水力学行为.
①
收稿日期: 2007 09 07
基金 项目: 国家科技攻关资助项 目( 863 项目, 2008AA06Z312) ; 国家自然科学基金资助项目(20477033) ; 教育部留学回国人员科研启动
基金资助项目( 教外司留 No畅2004 527) ; 重庆市自然科学基金重点资助项目(CSTC2004BA7019) .
作者简介: 戚景南(1 980 ) , 女, 河北唐山人, 硕士研究生, 主要从事环境化学的研究.
通讯作者: 黄玉明,教 授, 博士生导师.