基于NETWORK模型的研究底泥中重金属离子的迁移规律.pdfVIP

基于NETWORK 模型研究底泥中重金属离子的迁移规律1) * *2) 王振 ，林缅 * （中国科学院力学研究所，北京，100190） 摘要：在浅水区吸附在底泥中的重金属离子在一定的水动力条件下会随着微粒的上升、迁移、释放，而造成上覆水体的 二次污染。本文考虑微粒的筛滤、表面沉积两种俘获机制，应用NETWORK 模型研究了上覆水体中波浪运动影响下多孔 介质中微粒的迁移规律。研究结果表明，波浪运动造成了微粒在底泥多孔介质中的弥散；弥散和俘获造成了微粒的向上 净迁移。 关键词：NETWORK 模型；筛滤；表面沉积；微粒迁移；底泥多孔介质 引言 随着沿海地区的迅猛发展，越来越多的含重金属离子的污染物被排放到海中，这些含重金属污染物， [1] 一部分当下会被海水带走，还有一部分会吸附在泥沙微粒上而慢慢沉积在底泥中 。在一定的水动力条件 下，这些沉积在底泥中的重金属会随着微粒的迁移释放造成上覆水体的二次污染。所谓微粒是指直径小于 [2] 10 m 的小颗粒，其表面带有电荷，具有比较大的比表面积，对重金属离子具有较强的吸附作用 ，所以 重金属离子是随着微粒的迁移而迁移的。 目前关于重金属离子的迁移扩散的研究大多针对的是土壤的污染问题，相关模拟方法大致可分为以下 两类： [2-4] 1）连续介质模型。基于质量守恒原理，连续介质模型采用对流-弥散方程 （CDE ） 研究微粒的迁移， 微粒的释放和俘获以源 （汇）项的形式表述。该方法忽略了多孔介质孔隙空间的拓扑结构和连通性等，难 以给出孔隙尺度上微粒释放和俘获的物理描述。 2 ）离散模型。Rege 和Fogler[5]发展了一个二维三角形NETWORK 模型，对微粒筛滤效应 （微粒直径 比孔隙直径大，从而堵塞孔隙）进行了研究，其模型考虑了流体流动对微粒路径选择的影响。Idakm 和 Sahimi[6]建立了三维正方体NETWORK 模型，较为准确的预测了由微粒筛滤和表面沉积两种效应引起的多 孔介质渗透率下降。 综合来看，NETWORK 模型可以反映孔隙空间的拓扑结构和连通性，也便于考虑多孔介质的非均质性， 因而在研究孔隙尺度上微粒迁移、释放和俘获过程方面比连续介质模型更有优势。本文采用 NETWORK 模型并运用粒子追踪的方法模拟微粒在多孔介质中的迁移，综合考虑了微粒筛滤和表面沉积两种俘获机 制。运用该方法本文对上覆水体波浪作用下底泥多孔介质中微粒和重金属的迁移规律进行了探索性的研 究。 1) 自然科学基金资助项目 2) E-mail ：linmian@imech.ac.cn 1 1 模型建立 1.1 计算模型构建 我们知道，NETWORK 模型是否规则并不影响微粒的迁移性质，只要平均配位数一致，其输运特性在 [7] 本质上都是一样的 ，因此为了便于研究微粒迁移机理，本文采用二维正方形20 40 的NETWORK 模型， 配位数为4，即每个node 都有4 个bond 与之相连。假设键为圆形，其半径 根据Baghdikian[8] 的实验测 R 得的孔隙半径给定，满足对数正态分布，其均值和标准差分别为 18m 和12m ，键长度为定值l ，其 值可以通过下式计算得到[8] Ri