地下水非平稳随机模型及空间变异性与非均匀性相互关系研究的.pdfVIP

水 利 学 报 2004 年10 月 SHUILI XUEBAO 第10 期 文章编号: 0559-9350(2004)10-0013-09 地下水非平稳随机模型及空间变异性与非均匀性相互关系研 究的展望 1 1 2 廖华胜，李连侠，LI Shu-guang (1.四川大学 高速水力学国家重点实验室，四川 成都 610065；2. Dept. of Civil Environmental Eng. Michigan State University, East Lansing, MI 48824,USA) 摘要：本文在对地下水空间变异性现场的调查研究、大尺度示踪试验和数值模拟研究的结果进行比较分析的基础 上，指出平稳随机的假定不能真实反映空间小尺度变异性与大尺度非平稳性间的相互作用，而现有的随机方法需 要分辨小尺度的动力特性，因此在实际应用中受到极大限制。本文还从CPU、内存需求和方程求解难易程度等角度 出发，在客观评述现有地下水随机模型基础上，对非平稳谱方法在研究变异性与非平稳性相互作用方面的应用前 景作了初步展望，认为该方法将使得对实际现场尺度下的流动和输运过程的直接模拟成为可能。 关键词：地下水；随机模型；非平稳谱方法；变异性；平稳性 中图分类号：P641.2 文献标识码：A 1 动机和背景 目前都普遍认同自然的地下水环境其空间变异性是非常显著的。特别是介质的变异性会引起地下水含 水层点的水力传导系数的剧变。尽管有时某些岩层其分层趋势良好，但这种变化有时也是随机的。事实上， 水力传导系数的小尺度脉动会产生明显的大尺度效应。主要的是因为水力传导系数与污染物浓度之间存在 着非线性关系。这种小尺度空间变异性对大尺度溶质运移的影响被认为是与宏尺度弥散现象相关联的[1]， 即实际现场实测到的溶质扩散范围比室内土柱试验的更大的现象[2]。 多年来，对地下水空间变异性作用的研究，已有诸多随机方法[3-8]。总体而言，这些方法都假定变异性 的物理参数（如水力传导系数）在空间上是一已知其统计特性的随机函数。由此得出随这一参数变化的环 境变量也是随机的。比如，水力传导系数不确定性将导致孔隙速度的不确定性，因而导致溶质浓度的不确 定性。流动和输运方程正是联系随机因变量和非确定性水文地质参数的基础。可是，在实际应用中，如果 不作出这样或那样的假定，将很难导出这些方程。因此，随机地下水模拟的技术，就在于了解对给定的实 际例子，做什么样的近似才最为合适。 在随机地下水分析中有一个最重要假定—统计平稳性假定。统计平稳性是描述随机函数时间和空间原 点平移时其统计参数的变化特性。一个空间随机函数只有它的均值和协方差与绝对位置无关时，才能被称 为是广义（或弱的）平稳随机函数。而大多数的水文和地质变量并不是真正地具有平稳性。因为，它们的 均值可以在所感兴趣的尺寸范围内变化。如果这些均值的变化与在该均值上下变化的随机脉动量相对为一 收稿日期：2003-10-15 基金项目：本文得到教育部“跨世纪优秀人才培养计划”基金资助 作者简介：廖华胜(1964-)，男，广西田东人，教授，主要研究领域为工程水力学，地下水数值模拟，可视化应用软件 开发等。 1 水 利 学 报 2004 年10 月 SHUILI XUEBAO 第10 期 小量的话（即在几个相关尺度内此均值接近常数），作为近似或“局部”意义而言，统计平稳性假定是成 立的。这种“局部”平稳概念在地下水应用中也经常遇到[9]。 统计平稳性假定在方法和概念上均具有优越性。它使得我们可以利用各种不同的解析方法来求解随机 地下水流动与输运问题[4～6] ；更重要的是，它可让我们借助各态历经假设，在空间变异性与非确定性之间 建立等价关系[10]。采用平稳随机理论进行实际现场尺度下流动与输运问题研