多尺度非均质多孔介质中溶质运移的蒙特卡罗模拟.pdf

第 19 卷 第 3 期 水 科 学 进 展 Vol19 ,No3 2008 年 5 月 　　 ADVANCES IN WATER SCIENCE 　 May , 2008 　 多尺度非均质多孔介质中溶质运移的蒙特卡罗模拟 陆　乐 , 吴吉春 , 王晶晶 (南京大学水科学系 , 江苏 南京 　210093) ( ) 摘要 : 探讨了将蒙特卡罗 Monte Carlo 方法应用于多尺度非均质含水层中溶质运移模拟的方法 。所研究的含水层由 两种具有不同渗透系数统计特征的多孔介质所组成 , 每一种多孔介质是非均质的 , 且其渗透系数场符合平稳假设 , 而整个模拟区的渗透系数是非平稳的。Monte Carlo 方法要求参数是平稳的 , 因此 , 分别对两种多孔介质产生若干随 机渗透系数场后 , 用两种方法进行组合 , 并进行溶质运移的模拟计算 。通过对计算结果的分析 , 综合考虑计算精 度 、计算时间等因素 , 得出了处理多尺度非均质多孔介质中溶质运移问题的较好方法 。 关 　键 　词 : 多尺度 ; 多孔介质 ; 蒙特卡罗 ; 溶质运移 ; 渗透系数 ; 空间变异性 ( ) 中图分类号 : P641 　　　文献标识码 : A 　　　文章编号 : 2008 自1970 年以来 , 数值模拟方法已逐渐被应用到水文地质的各个领域 , 时至今 日, 地下水数值模拟已成为 评价人类活动对地下水水质及水量的影响、合理开发利用地下水以及预测地下水污染发展趋势的最主要的方法 和手段[ 1] 。确定性方法是进行地下水数值模拟的传统方法 , 已经发展到比较成熟的阶段 。对于水文地质参数的 处理 , 传统的确定性方法都是在研究区内根据实际的水文地质条件进行参数分区 , 并且认为同一参数分区内的 水文地质参数是相同的数值 。但实际上 , 由于天然含水层沉积过程的随机性 , 含水层的水文地质特征以及结 构 、构造等具有非常明显的空间变异性 , 而且这种空间变异性往往在不同尺度上表现出来 : 在小块的土样中 , 可以观测到土粒孔隙大小及团粒结构的差异 ; 在大的区域含水层系统中 , 可以观测到不同含水层间以及同一含 水层的不同地域间的变异性 。众多研究结果表明 , 水文地质参数尤其是渗透系数的空间变异性是影响含水层中 溶质运移的决定性因素[2 ] 。因此 , Dagan 等提出并发展了随机理论的研究方法[3～13 ] 。目前随机方法已成为研究 非均质含水层中地下水流和溶质运移 问题的重要手段 。地下水随机模拟方法主要有矩方程法和 Monte Carlo 法[ 14 ,15 ] 。矩方程方法通过求解有关均值和协方差的随机偏微分方程获得随机问题的解 , 而 Monte Carlo 法则是 通过平均一系列反映含水层实际性质的确定性问题来模拟随机过程的一种计算机模拟方法 。 随机运移理论依赖于水文地质参数满足平稳和遍历性的假设 , 这些假设隐含的意思是非均质介质的结构可 以看成是在空间上具有周期性的 , 从而得到有限的相关尺度[ 16 ] 。在应用随机理论进行地下水溶质运移数值模 拟时 , 通常假设整个研究区域内水文地质参数的空间分布可以由单一的相关