地下水流数值模型的分解与叠加.pdfVIP

第 37卷 第 1期 水文地质工程地质 Vo1．37 NO．1 2010年 1月 HYDR0GEOLOGY ENGINEERING GEOL0GY Jan．20l0 地下水流数值模型的分解与叠加 王晓明，黄 新，阁黎 东 (西安科技大学地质与环境学院，西安 710054) 摘要 ：叠加原理在水文地质计算 中的应用相 当普遍 ，但一般局 限于地下水流模型 的解析解 ，本文扼要介绍 了线性定解 问 题 的叠加原理 ，在此基础上提 出了线性离散模 型的叠加原理 ，进而讨论 了线性离散模型的叠加原理在地下水流数值模型 中的应用一一地下水流数值模型的分解与叠加。模型分解与叠加 既可用于解析分析也可用于数值计算 ；同时不仅可 以 简化复杂模型的求解问题 ，更有助于理解地下水变化过程 。 关键词 ：地下水流数值模型 ；模 型分解 ；模型叠加 ；离散模型 中图分类号 ：P641．2 文献标识码 ：A 文章编号 ：1000—3665(2010)01-003l-03 在水文地质计算 中，如解决群井干扰 、边界附近井 线性方程组 。为便于说 明问题，这里讨论一般情况 流、变流量等问题时 ，都用到了叠加原理。但叠加原理 下线性方程组的分解与解 的叠加问题。考虑如下 的n 的上述应用仅局限于地下水流模型的解析解 ，在数值 ×n阶非齐次线性方程组 ： all 1 + a12 2 + … + 0l = bl+c1 计算中却容易被忽略。实际上 ，叠加原理不但可以用 Ⅱ a21 t + a212gt72 + … + 02 = b2+ c2 于解析解 ，也可以用于数值模 型。合理地应用叠加原 n 理不仅有助于复杂模型的求解 ，更有助于对水文地质 a Ⅱl 1 + a 2 2 + … + a = bn+cn 条件 ，地下水变化过程的理解 。本文试 图对离散模 型 n n n 叠加原理的基本概念以及地下水流数值模型分解与叠 写成矩阵形式为 ： [A]{}= {b}+ {c} (1) 加的方法进行简单 的阐述 ，并着重讨论离散模 型的叠 将非齐次线性方程组 (1)分解为如下 的两个非齐 加原理在地下水流数值模 型中的应用。 次线性方程组 1 离散模型的叠加原理 [A]{}= {b} (2) [A]{}= {c} (3) 如果泛定方程和定解条件都是线性 的，则可把定 注意线性 方程组 (1)，(2)，(3)的系数矩阵都一 解问题看作是几个部分的线性叠加 ，只要这些部分各 样 ，所不 同的是常数项 ，且方程组 (1)的常数项为方程 自所满足的泛定方程和定解条件线性叠加 ，满足原来