溶质在多孔介质中的有效扩散系数分度结构

《溶质在多孔介质中的有效扩散系数分度结构》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、环境科学--13卷第二期、70-172页、2001年====================================================================论文ID：1001-0742（2001）02-0170-03CLC号：X34文件代码：A溶质在多孔介质中的有效扩散系数分度结构摘要：分形方法用来辅助推敲在多孔介质中的有效扩散系数和多孔介质中几何参数表征的幂律关系。该结果和Archie法则中的经验方程相似并预计将会被应用于预测有效扩散系数。关键词：扩散、有效扩散系数、分形、多孔介质引言对于低渗透性的多孔介质（作为自然或人工屏障

2、被广泛应用于废物处置场所的污染控制中），扩散传质机理可能起到主导地位。污染物的扩散传质通过粘土、淤泥或页岩结构以及矿石层和合成聚合物这些因素都被作为选址、风险评估和垃圾填埋场工程设计的审议。此外，废弃污染物的过滤取决于污染物的分子扩散。模拟溶质在多孔介质中的扩散需要有效扩散系数的信息，有效扩散的数值又与驱动力的变化有关。测量有效扩散系数是费时又充满变数的，特别是在地渗透性的多孔介质，所以参数的预测是可行的选择。最经典模拟多孔介质扩散的模型是基于一种假设，即通传递过程是不变的。例如，介质看起来在随机的不同的地点可能是由一个有限的样本大小或固定的统计处理技术。现代

3、分形模型，在另一方面，假定介质的膨胀不变，即看起来放大了相同的倍率。已发现天然的多孔介质具有分形层次结构，并遵循比例和规模之间的上限和下限幂律分布的特点（Gimenez，1997）。Gimenez回顾了经典的半经验法则如Archie法则和Campbell法则，假定分形标度的多孔介质的各种物理特性，应用这些法则去预测土壤中保留水和水力传导性。本篇论文的目的是用分形方法推导一个幂律方程用于预测稳定状态下的有效扩散系数。1.有效扩散系数的定义溶质在多孔介质中扩散被迂回曲折的毛孔结构所阻碍，有效横截面积（孔隙度）和孔径的分布。在稳定状态下，质量通量取决于集中梯度并被F

4、ick第一定律描述：在水饱和的条件下C是与指溶质在水中的浓度有关，De，有效扩散系数被定义为：Daq是水相扩散系数，ε是多孔介质有效空隙率用于说明在扩散中可能减少的有效截面积，当扩散只发生在空隙间时。自然的多孔介质经常还含有溶质无法通过的更小的孔隙且不影响整体通知的移动如死结点或盲孔，ε小于多孔介质的整体孔隙度（Φ）（Lever,1985）.是曲折因子，是说明多孔介质结构并定义为孔长度Le(有效扩散路径)相对与多孔介质的沿主要流程或扩散轴方向的长度L之比的平方。（Epstein,1989）τ是曲折因子。一般来说Le>1,τ>1.实际上，是一个多孔介质一次元的毛

5、细管的参数模型而不是介质的性质（Grathwohl,1997）。毛细管模型（图1）假设多孔介质是蜿蜒的孔束且平行的排列着。图1多孔介质的曲折毛细管束模型2.有效扩散系数的缩放多数情况下，只有多孔介质的整体孔隙度（Φ）可以同时包含有效孔隙度ε和曲折系数τf这两个难以确定的量。类似于Archie法则——被广泛应用于描述多孔岩石电导率的经验法则，这个公式也成为多孔介质有效扩散系数De和整体孔隙率相关的经验幂定律。这里m是一个经验值。如果有效孔隙率ε近似于整体扩散系数Φ，那么ε≈Φ，那么曲折因子τf可与有效孔隙率ε的方程（2）及（5）：图2多孔介质的分形束模型由于方程

6、（5）和（6）是经验的、没有坚实的物理基础和数学基础，它们在已了解多孔介质的扩散现象中有最低的影响，这个限制了他们在数据组以外推断结果的可能性。就像Archie法则，幂律法则方程（5）和（6）表明了几何分形。在蜿蜒的毛细管模型，“曲折”性质反映了孔隙介质界面的粗糙度。深入探究（Gimenez，1997）曾建议，多孔介质的气孔界面具有一定规模的分形特性的限制。假设捆绑如此孔隙束的“曲折”就是一个分形（自相似性），例如VonKoch曲线（图2；Mandelbrot，1982），幂定律方程（5）和（6）可推导如下。鉴于一个长度为L2体积元素，孔隙分形意味着固界面的分

7、形维数和介质的整体孔隙度之间的一个简单关系（Katz，1985）：其中Ds是孔隙固体界面分形维数，L1和L2是分区的上限和下限。且1≤Ds≤2。该分形束的几何元素L2可通过如下方程确定：方程（7）和（8）代入方程（2）：同样假设ε近似等于Φ。方程（9）等价于方程（10）和（11）与经验公式(5)和（6）一致，方程（10）和（11）由方程（5）和（6）推导出。方程（12）的经验显示，指数m（%1）是一个描述多孔介质中界面粗糙度的参数。因此孔界面的粗糙度越大，m值越大。根据方程（11）和（12），Ds=1等价于m=1和τ=1；因此方程（10）写成它对应于多孔介质毛

8、细管束模型。一般来说，Daq可以在文献