地下水数值模拟07反求水文地质参数数值方法.pptVIP

第六章 反求水文地质参数的数值方法 、基本概念 ·正演问题与反演问题 正演问题(正问题): 在已知地下水流动的微分方程及水文地质参数 K、M、W、卩、μ*和边界条件的前提下,求解渗流区域 内的水头分布规律和流量。 反演问题(逆问题): 根据地下水的天然动态或抽水试验的观测资料 研究所选用的方程类型是否适当、确定方程的参数和检 验定解条件。 、基本概念 反求参数问题的适定性圖可以诫 解的存在性——解是否存在? 根据实际资料反求渗流区的水文地质参数(有时包括边界流 量),这样的解是否存在? 很难满 解的唯一性——解是否唯一?。。 根据实际资料反求的水文地质参数是否唯 解的稳定性—解是否稳定?。0很难满 当实测资料有微小误差时,反求的水文地质参数的误差是 否也微小,即水文地质参数是否连续依赖于实测资料 解的噍一性若单从观测数据反求参数,可能存在多种解!! OHOH 例1 均质等厚二维承压含 水层稳定流动问题 HR=H 在同一区域同一边界条件下,不管导水系数T取何值,水头分布相同 a(aH ◆正问题:H为未知量,T为已知量 例2 可以求得: C1→H=-x+C2 x=0 x=0→C2=H1 x=L→C1= T(H2-H1) ◆逆问题:T为未知量 若已知r 可以求得 OH OH →T C→T=C/ 需要补充相应条件以确定参数 解的稳定性水头较小的误差,可能会引起所求参数较大的误差! aH aH arcy2+W=0→7= 02H 若实测水头存在误差ε,即:H=H+E OHOH 88 88 E OHOH a8 O8YOH a2 从数学上讲:虽然很小,但其导数何能很大,若a+→H ax ay a-ha-h a-e ox Oy 数值模型反演问题解的适定性 ·以非均质承压二维非稳定流为例 1) 81)8 +1 如图所示的有限差分网格,其差分方程为 于1) H +2(4x 7=,+7,) H,-H T,:+T 1, H++I H ◆正问题:H为未知量 H H H [B{H}=f} T (4y) ◆逆问题:T、μ*、W为 未知量 A}={f 数值模型反演问题解的适定性 Akk= 根据方程个数N和未知参数个数M(m1+m2+m3)之间的关系: (1)NM: 方程组可能存在无穷多组解—不满足唯一性 (2)N=M 当系数矩阵行列式不等于0,能求出唯一解(需要足够多资料) 当系数矩阵行列式的值很小,则水头观测的微小变化可能造成参 数的很大变化—不满足稳定性 (3)NM: 一般不存在任何一组参数使所有方程同时满足。然而,根据问题本身的 物理特性,应该存在一组参数使这些方程基本得到满足 反求参数问题的适定性 ·反求参数问题本身不一定是唯一的,也不 定是稳定的。 但是对于实际问题,我们可以根据对水文 地质条件的初步认识以及通过其它手段得 到一些辅助的资料和参数的约束条件从而 使得反求参数问题在一定程度上是唯一的 和稳定的。 、反求参数的直接解法 求解思路: 指在地下水流动微分方程(或描述地下水流动的数值模型)中 ,将水头值作为已知量,将待求的参数(往往包括源汇项及 边界流量)等作为未知量,直接求解未知参数的方法。 常用方法: 局部直接求逆法 数学规划法 应用现状: 对数据误差十分敏感,对观测资料有过高的要求,因而目 前还难以应用 局部直接求逆法 例下→ aHaH aH 已知不同时刻局部区域水头H和源汇项W的实测值,需要求解T和μ 设,/02 H aH b (w) 22h ah H b2(w) 1 T+W,=u b a, T+w,=u Z r=b2W-b吗 2b-a, b 优点:不需要初始条件,也不需要边界条件。对均质各向同性、等厚的渗流区来说 只需要知道某个局部区域在两个不同时刻的水头值和垂直方向的水量交换的实测资料, 因此,这个方法称为“局部直接求逆法” 难点:如何求出方程中水头对时间和空间的导数。实际计算发现系数矩阵对应的行列 式之值通常很小,因此求参数住往会产生很大的误差