基于随机理论的渗透参数反演分析方法及应用.pdfVIP

基于随机理论的渗透参数反演分析方法及应用 魏海沈振中 江沆聂琴 (河海大学水利水电工程学院，江苏南京210098) 摘要：首先利用试验资料求出坝体和坝基材料渗透参数的概率分布，并结合钻孔资料对坝体和坝基进行分区， 建立有限元模型，根据坝体材料渗透参数的概率分布随机生成单元的渗透参数，然后采用DFP方法进行渗透 系数反分析，反演出坝体和坝基材料渗透参数的均值和标准差。该法应用于某平原地区水库大坝渗透参数的反 演分析，得到了较为符合实际的成果。算例表明，该法具有精度高、操作简便的优点，是一种切实可行的反演 方法。 关键词：随机理论渗透参数反演分析有限元 1 引言 岩土渗透参数的选择在岩土渗流分析中十分重要，参数选择的合理与否直接影响计算的合理性和可靠性。 岩土渗透参数通常是通过室内试验或现场试验测得，但室内试验由于岩土试样的扰动测出的参数并不能代表 真实的情况；而现场测试的费用高又不能大量进行，得到的参数有时就缺乏代表性，因此岩土渗透参数反演 分析被工程界普遍青睐。反演分析已是岩土工程中的一个重要分支。通常的反演分析是通过沉降量反演岩土 体的变形参数，通过渗流量和测压管水位反演岩土体渗透系数、贮水率和给水度等。通常反演的参数更能代 表真实情况，反映实际。但由于岩土的复杂性，有时同一种土的参数变化很大，而传统的方法是用一个参数 代表一类土层的性质，这不能反映岩土参数的一些变化，因此有时分析结果很难和实测数据吻合。本文基于 随机理论把岩土渗透参数视为某种概率分布，把坝体划分为许多单元，根据渗透参数的概率分布随机生成单 元的渗透参数，然后采用DFP方法进行岩土渗透系数反分析，反演出岩土体渗透参数的均值和标准差，并反 演出坝体的渗流场。 2 随机理论 随机理论的数学基础主要是概率论和数理统计，它是利用不确定的数学方法去研究工程中的确定性问题 和随机问题。随机理论包括随机变量函数、参数估计、随机变量分布模型、回归分析与相关分析、抽样检验、 质的孑L隙在空间上的分布通常都是随机的，于是土层的有效孔隙度、贮水系数、渗透参数等在空间上也呈随 机分布，致使地下水在多孔介质中的运动具有很大的不确定性，因此大坝的渗流量和渗流场分布是许多随机 因素共同作用的结果。传统的确定性数学方法很难反映这些因素的随机性，而随机理论方法却能很好体现这 些因素的随机性【21。随机理论求解随机问题时。一般首先根据问题的物理性质建立随机模型，然后再根据模 型中各个随机变量的分布，求出问题的大量统计试验值，最后由这些试验结果求出它的统计特征值【3】。 3参数随机反演分析方法 ． 随机反演分析通常将材料参数、结构的响应及外载视为随机变量，然后根据结构响应的概率分布建立目 标函数，再采用确定化优化反演方法求出需要反演的参数。一般随机反演方法有极大似然参数反演和Bayes 参数反演法。 极大似然参数反演法用于分析大坝渗透参数时，是把实测值(水头h’)视为随机变量，则幽’(水头变化) 也具有随机性，于是幽’的概率分布是大坝材料渗透参数k的条件分布，可用f(AhTk)来描述。极大似然法 的思想，是把待求参数k取某一确定值时，计算值△J，恰好为实测值幽’的概率最大，即找某一k使得条件概 率厂(幽．／七)取最大值14J。 ·162· 极大似然参数反演时。没有考虑坝体渗透参数k的不确定性。而Bayes参数反演法可以将渗透参数视为 随机变量，充分利用工程经验大致确定的参数先验信息。在利用Bayes参数反演时，先假设渗透参数k的概 度函数f(k／Ah9，当f(k／劬f)取最大值时的k就是待求材料的渗透参数【5，6】。 Bayes参数反演法认为大坝的水头是随机变化的，考虑了坝体渗透参数的变化，但它把坝体视为绝对均 匀，没有考虑坝体渗透参数在空间上的变化。本文首先根据钻孔资料对坝体进行分区，这种分区称为物理分 区，它与实际的地质情况一致。在物理分区的基础上再进行数学分区，数学分区是把各物理层(即地层)进行 再次划分。各分区的形状可以是三角形、四边形或其他多边形(本文采用三角形和四边形)，分区的大小可以 根据问题的复杂情况和计算精度要求确定。在数学分区后，就可以根据各物理分区渗透参数的概率分布类型， 随机生成各数学分区的渗透参数。在随机生成渗透参数时，通常根据经验或试验确定的渗透参数的均值和标 准差。在确定坝体的渗透参数后，就可以利用稳定渗流方程或非稳定流方程进行渗流分析，采用常规参数反 演方法进行参数反演【J71。由于各数学分区