kl展开在渗流场随机分析中的初步应用∫d-数学进展-长江科学院.pdfVIP

　　第２６卷第１０期 长　江　科　学　院　院　报 Ｖｏｌ．２６　Ｎｏ．１０　　 　　２００９年１０月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ｏｃｔ．２００９　 　 　　文章编号：１００１－５４８５（２００９）１０－００３９－０５ ＫＬ展开在渗流场随机分析中的初步应用 １ ２ １ 李少龙 ，杨金忠 ，张家发 （１．长江科学院水利部岩土力学与工程重点实验室，武汉　４３００１０； ２．武汉大学 水资源与水电工程科学国家重点实验室，武汉　４３００７２） 摘要：在水力特性参数空间变异性的影响下，渗流场具有一定程度的不确定性。分析了随机场ＫＬ展开的特点，阐 述ＫＬＭＥ渗流随机分析方法的基本原理，发现利用ＫＬ展开的均方收敛性和摄动方法的渐进收敛性，并将随机偏微 分方程转化为递推的确定性偏微分方程组求解，提高了计算效率；对非饱和渗流场进行随机数值模拟，考虑了非饱 和渗透函数的参数场之间不相关、正相关和负相关３种情形，分析表明参数场之间的互相关性对渗流场的统计特 征具有较大影响。 关　键　词：随机场；渗流；ＫａｒｈｕｎｅｎＬｏｅｖｅ展开；空间变异性 中图分类号：ＴＶ１３９．１　　　文献标识码：Ａ 　　为了准确反映土性参数的空间变异性，一般将 非饱和渗流场的统计特征进行分析。 渗透系数等参数作为空间随机场。将参数随机场和 反映物理规律的微分方程相结合并发展相应的数值 １　随机场ＫＬ展开 方法，成为随机分析的有力工具。文献［１］将渗透 系数作为随机场，采用蒙特卡罗随机有限元法研究 随机场可以表示为级数展开，即用确定性函数 了土体渗透系数的随机性对有压渗流场统计特征的 的一个完备集与相应的随机系数来表征。在随机场 影响。文献［２］将渗透系数作为随机场，应用小参 的级数展开中，较为常用的是 ＫＬ展开。随机场 Ｕ（Ｘ，）的ＫＬ展开式为 数摄动法得到稳定渗流的随机有限元列式，通过算 ω 例分析水力坡降标准差与渗透系数变异性之间的关 珚 ∞ Ｕ（ｘ，）＝Ｕ（ｘ）＋ （）（ｘ）　，（１） ω Σχξ ω 槡ｉ ｉ ｉ 系。以上分析仅考虑了一个参数的空间变异性，实 ｉ＝１ 珚 式中：Ｕ（ｘ）为随机场的均值函数；｛ （ ）｝为零均 ξ ω 际问题中通常是多个参数共同作用。 ｉ 值互不相关的随机变量序列；（ｘ）和 分别为随 χ 为了将参数随机场引入数值分析方法中，一般