岩体渗透系数反演的数值方法及适定性.pdfVIP

岩体渗透系数反演的数值方法及其适定性① 刘迎曦李守巨 周承芳 李正国 林龙泽 (大连理工大学大连 116024) (丰满发电厂吉林132108) 路振刚 (东北电业管理局沈阳110006) 摘要本文分析了渗透系数反演中的适定性向题，包括唯一性和稳定性同题。分别对直接反演 法和间接反演法进行丁论述。侧重研究了同接法中的正则化最小二乘法和逐个修正法，建立了 两套反演数值计算方法，主要内容包括叠代格式的建立，Jacobian矩阵的计算．高斯一牛顿方向 的修正。计算结果表明，此两种方法有较快的收敛速度和较强的稳定性。 关键词反分析 渗流最小二乘法适定性 高斯一牛顿法 1渗流问题反演的适定性 根据水头的观测资料反求岩土水文地质参数问题被称为解渗流反问题渖f究表明，地下 水渗流模型的正问题(Iowardproblem)都是适定的，而反问题(inverseproblem)往往是不适 定的”]。不适定一般可由被识别参数的不唯一和不稳定来表征。反问题的不稳定是指求头的 少量误差会引起被识别参数的严重误差。不唯一性是指被识别参数的数值会随着初始参数 的改变而变化。参数识别中的唯一性与参数的可识别性密切相关。可识别性提出了这样一 个问题，根据所测得的水头数据总是能够获得反问题未知参数的唯一解。如果在研究范围的 所有点，在使用系统的输出一输入关系和输出一输入数据时，未知参数能够被唯一确定，则称 未知参数是可识别的。地下水渗流问题的基本偏微分方程属于抛物线型的热传导方程，这类 方程的正问题是适定的，其解有很好的稳定性。含水层参数的少量变化，对求出的水头影响 很小。但是，渗流的反问题是非适定的，即水头观测上的少量误差及不可避免的截断误差会 使反求出的参数有较大幅度的跳动，甚至出现不应有的负值。无论有限差分法还是有限元 式。’可以看出，对于正演问题，水位h为待求未知数。^的系数矩阵主对角线约为其它非零 元素的4倍多，这样的系数矩阵行列式值较大．导水系数71的少量变化对解影响很小。将 题时，T成了未知数，系数矩阵主对角优势没有了，解的稳定性很差。反问题的经典定义为一 ④国家自然科学基金资助项目，项目批准号；5977903 397 丁(j)，即铒存在；对于任何给定的h(z，￡)，解唯一}解连续地依赖于h(x，f)，即解是稳定的 的。如果反问题不能满足此三条的任何一条，则称反声题是不适定的(ill—posed)。 2渗透系数反演的直接方法 按照Neuman的分类方法。’，渗流参数识别反问题可分为两类：直接法和间接法．如果 水头的导数在整个渗流区域已知(常常采用数值近似，如有限差分和有限元法)，初始的支配 方程(Darcy定律)就成为包含有待识别参数的线性一阶偏微分方程，借助边界条件，可以直 接解出未知的参数(渗透系数了1)。该种方法被称为直接法。间接法是先给定待定水文地质参 数71假定初值，通过解正演模型计算相应的水头分布，然后将计算水头值与实测水头值进 行比较，看二者拟合程度如何，如果不满足要求，则修改参数，重新进行计算。直到满足要求 为止。无论是直接法还是间接法，反问题最终都归结为数学规划中的求极小化平方和函数问 题或极小化残差问题。在过去描述参数系统时，一直使用两种误差准则：输出误差准则(Out． Error Error put Criterion)(间接法)和方程误差准则(EquationCriterion)(直接法)。如果水头 变化和导数(常常靠估计)在整个渗流区域已知，并且测量误差和模型误差可以忽略不计，原 来的控制方程就变为双曲型线性一阶偏微分土‘程。借助于边界条件和流量数据，直接求解未 知参数是可能的。实际上，水头观测点是随意散步在渗流区域的，并且观{贝l点是有限的。为 j。采用方程误差准则进行反分析，某些点的观测数据(已丢失或无效)不得不采用插值法进 行估计。在插值过程中，所得到的数据含有误差。如果伴随着观测的插值数据被代人控制方 程，并且观测数据带有噪音。这样．在计算时误差项就出现了。这样的误差被称为方程误差。 然而，通过采用台适的参数使该误差降至最小。如果不考虑样士的统计