非均质土体中渗流场的非平稳随机模拟.pdfVIP

第26卷第10期 长 江 科 学 院 院 报 Vo1．26 No．10 2009年 10月 JournalofYangtzeRiverScientificResearchInstitute Oct．2 0 09 文章编号：1001—5485(2009)10—0049—05 非均质土体中渗流场的非平稳随机模拟 李少龙 。张家发 (长江科学院 水利部岩土力学与工程重点实验室 ，武汉 430010) 摘要：将渗透系数作为对数正态非平稳随机场，分解为确定性趋势和平稳性扰动两部分，采用基于随机场 Kar- hunen—Loeve展开和摄动分析的KLME随机数值方法，探讨参数场的连续型线性趋势和间断型趋势项对渗流场统计 特征的影响。分析表明，水头均值和标准差的空间分布受到参数场趋势性变化的影响，随线性趋势的变化方向而 呈不同分布，线性趋势的斜率增大，水头统计矩的二阶修正项的作用也增大，在参数场间断型趋势条件下，水头统 计矩呈复杂的非规则空间分布。 关 键 词：随机场 ；渗流 ；非平稳性；随机模拟 中图分类号TB115 文献标识码 ：A 实际渗流场受众多随机因素的影响因而具有一 足平稳场的均值为常量的要求，此时需要考虑非平 定程度的不确定性。近年来，渗流场随机分析得到 稳的随机渗流分析。文献[6]采用一阶非平稳谱方 关注，较多的研究是在考虑渗流参数随机性的基础 法研究了渗透系数具有趋势性变化时渗流场的统计 上将渗流场的不确定性定量化，分析随机渗流场的 特征，为了便于推导解析解而没有考虑边界的影响， 分布规律及结果的可靠性。文献 [1]将各渗透分区 并假定趋势变化的方向与流动方向一致。文献 [7] 的参数作为随机变量，采用一阶Taylor展开随机有 用一阶矩方程法分析了对数渗透系数场的均值函数 限元法对地下水承压流进行随机模拟。文献 [2]将 呈线性、二次和周期性趋势变化时水头方差的异同。 裂隙系统的几何参数作为随机变量，用一阶Taylor 本文将渗透系数作为非平稳随机场，考虑均值函数 展开随机有限元法分析了裂隙岩体渗流场的随机 的连续型和间断型两种空间趋势变化，应用 KLME 性。将岩土体参数作为随机变量，实际上忽略了它 方法计算水头统计矩的二阶结果，探讨参数场的非 们的空间相关性，不符合实际情况。文献 [3]将渗 平稳性对渗流场统计特征的影响。 透系数作为平稳随机场，应用小参数摄动法得到稳 定渗流的随机有限元列式，通过算例分析了水力坡 1 渗透系数随机场 降标准差与渗透系数变异性之间的关系。文献 [4] 将渗透系数作为平稳 随机场，采用 Monte—Cado 土体渗透系数在空间上的变化，既具有随机性 (MC)随机有限元法，研究了土体渗透系数的随机 又具有结构性 ，对这种空间变异的一种处理方式就 性对有压渗流场统计特征的影响。文献 [5]提出了 是用随机场进行描述，如文献 [8]提出土性剖面的 可用于高阶分析的KLME方法，将渗透系数作为对 随机场模型，考虑了土性空间变异的两重性，在岩土 数正态平稳场 ，通过有压渗流的算例与ME方法进 工程 中得到较大应用 。 行比较，分析表明该方法具有较高的计算精度和效 应用随机理论分析问题时通常会引入一些假 率 。 定，较常用的一个假定就是将参数随机场作为二阶 土性参数具有随机性和空问相关性 ，用随机场 平稳场，即：均值在空间上为常数 ；协方差函数仅依 描述这种空间变异性是比较合适的，很多随机渗流 赖于空间两点之间的