含水层渗透系数空间变异性对地下水数值模拟的影响 - 水科学进展.PDFVIP

第 16 卷 第 4 期 水 科 学 进 展 Vol16 ,No4 2005 年 7 月 　　 ADVANCES IN WATER SCIENCE 　 Jul . , 2005 　 含水层渗透系数空间变异性对地下水数值模拟的影响 陈　彦 , 吴吉春 (南京大学地球科学系 , 江苏 南京 　210093) 摘要 : 地下水数值模拟是 目前定量研究地下水水量和水质的重要手段 。使用基于随机理论的Monte Carlo 方法来进行 地下水数值模拟 。这种方法能较好地考虑水文地质参数的空间变异性 。主要将 Monte Carlo 方法和确定性模型模拟方 法的模拟结果在渗透系数场 、水头场 、速度场和浓度场等方面进行了比较 。结果表明: 在模拟三维非均质含水层中 的溶质运移问题时 , 充分考虑了含水层渗透系数空间变异性的Monte Carlo 法比确定性方法更为有效 , 模拟精度提高 了很多 , 且对模拟误差及误差来源有合理的数学解释 。 关 　键 　词 : 随机理论 ; Monte Carlo 模拟方法 ; 空间变异性 ; 地下水数值模拟 ; 含水层 ; 渗透系数 ( ) 中图分类号 : P641 　　　文献标识码 : A 　　　文章编号 : 2005 目前地下水数值模拟已成为评价人类活动对地下水水质及水量的影响、合理开发利用地下水以及预测地下 水污染发展趋势的最主要的方法和手段[ 1 ,2 ] 。由于地下水污染问题在工农业快速发展的今天已经相当严重 , 因 此 , 建立正确合理的数学模型来预测污染物的运移范围 , 进而采取有效的治理手段就成为必然的要求 。 目前构建地下水数学模型面临的最大困难是模型参数的选择 , 例如含水层渗透系数的确定 。野外实际条件 下的含水介质表现出强烈的空间变异性 。含水介质复杂的变异性迫使我们必须把含水介质作为一种随机变量处 理 。但由于含水介质的一些客观性质又表明它们并非只具有随机性 , 而是同时也具有一定程度上的结构性 。例 如 , 含水介质是在特定的历史时期特定的地质条件下沉积形成的 , 具有一定的成层性 。含水介质的上述特性使 得在考虑它的各种参数性质时 , 必须充分认识到它的随机性以及结构性 , 也就是含水介质的空间变异性 。 地下水数值模拟常用的方法包括有限单元法和有限差分法 。这两种方法的共同点是在进行模拟计算时都要 将复杂的研究区离散为很多小的单元或者网格 , 并根据实际的水文地质条件在这些单元和网格的基础上进行参 数分区。对于在同一参数分区内的单元认为其水文地质参数是相同的常数值 。当研究区面积不大或参数分区较 多时 , 这种方法是有一定精度的。但通常实际研究区的面积很大 , 即使参数分区已经很多 , 每个分区的面积相 对来说还是比较大 。把这样一个较大面积区域的水文地质参数设定为一个相同的数值就很难反映出实际含水介 质的空间变异性 。 在地下水数值模拟中通常使用模型校正的手段来确定水文地质参数 , 这种方法需要具备足够的实际观测资 料 。但在实际工作中获得足够的观测资料是比较困难的 , 因此通过模型校正得到的参数只能被认为是曲线拟合 参数 。大量的野外大尺度溶质运移试验的研究成果表明 , 渗透介质的空间变异性是影响溶质运动的决定性因 素[3 ] 。为了解决水文地质参数的这一特性 , Dagan 等提出了随机理论的研究方法[4～12 ] 。目前随机方法已成为研 究非均质含水层中地