浅析嵌入法和响应矩阵法建立水资源管理模型各自特点与联系.docVIP

浅析嵌入法和响应矩阵法建立水资源管理模型各自特点与联系 　　[摘要]建立地下水资源管理模型通常有两种方法：嵌入法和响应矩阵法，二者适用于解决不同类型的地下水资源管理问题，本文介绍了嵌入法和响应矩阵法的原理方法，通过三个实例模型分别使用两种建模方法求解，对两种方法所建模型的目标函数、约束条件和变量类型、数量、计算速度及输出结果进行对比分析，总结出两种建模方法的适用对象范围、各自特点和联系。 　　[关键词]嵌入法 响应矩阵法 水资源管理模型 特点 联系 　　[中图分类号] P641.83 [文献码] B [文章编号] 1000-405X（2015）-2-200-4 　　0引言 　　嵌入法和响应矩阵法是建立地下水资源管理模型的两种基本方法，经过国内外众多学者的研究探索与实践，嵌入法和响应矩阵法在地下水资源管理领域得到广泛的应用，目前已经发展到使用大型软件建立模型求解的阶段。在理论方法研究方面，关于两种建模方法的适用对象领域、限制条件、应用范围等，已经形成了一定的共识，同时也存在对两种建模方法的认识误区。本文首先从理论方法方面对这两种建模方法进行简述，然后通过三个实例分别采用两种方法建立模型并计算求解，从建模方法步骤、公式、计算结果、运算速度等方面进行对比分析，得出对两种建模方法的总结性认识。 　　1建模理论方法简述 　　1.1嵌入法 　　嵌入法是把地下水分布参数模拟模型以“嵌入”的方式作为管理模型的等式约束条件，建立地下水管理模型的方法。它把地下水水流模拟模型用一个线性方程组表示，然后把这个方程组作为优化管理模型的约束条件，实现了水流模拟模型与管理模型的耦合[1]。 　　以承压二维地下水流为例，其偏微分方程的离散矩阵形式如下： 　　式中： [A]为线性方程组的系数矩阵；H为由状态变量水头或水头降深构成的列向量；[B]为一个单位对角矩阵，存储各结点进出含水层水量的系数；[Q]是由可控脉冲变量和不可控脉冲变量构成向量；[D]为方程右端常数项。 　　将式（1）方程组作为约束条件“嵌入”到优化管理模型中，形成地下水管理模型： 　　目标函数： 　　约束条件： 　　式中：max（min）Z为目标管理项，取极大值或极小值；[C]为决策变量的价值系数向量；F（H，Q）为决策变量H或Q构成的函数向量；Gi（H，Q）为约束条件表示式向量；G0i为为右端常数项向量；H，Q≥0为非负约束。 　　1.2响应矩阵法 　　响应矩阵法是以线性系统的叠加原理和地下水位响应函数为基础，以矩阵形式反映地下水系统本身特征的一系列函数（单位脉冲响应函数），它表示地下水流系统中，水位与抽（注）水量之间的关系，由地下水流的模拟模型求得，作为优化管理模型中的一组优化约束条件，来实现水流模拟模型与优化模型的耦合[2]。 　　线性地下水系统可用单位脉冲响应函数来表征其特征，并通过它获得系统的输入输出关系表达式。线性地下水系统通常描述为： 　　抽水作为系统输入，地下水水位作为系统输出，其输入输出关系可用卷积的形式表示为： 　　目前，主要采用数值法建立地下水流模拟模型，在齐次初始边界渗流场条件下，通过抽水后观测剩余降深值，可求得单位脉冲响应函数。对于初始边界值不为零的非齐次条件，或存在天然补给（排泄）时，可用一个附加降深场S0（i，n）来模拟初始流场、边界条件及天然补给排泄对实际流场的影响[3]，此时降深场满足如下叠加关系： 　　式中：S（i，n）为i点n时刻的总降深；S0（i，n）为i点在非齐次条件或有天然补给排泄条件下，n时段实际流场的降深。β（i，j，n-k+1）为j点在k时段抽水，在n时段末在i点引起的水位降深值；Q（j，k）为j点k阶段相对于脉冲抽水量的抽水量，它以脉冲抽水量（I）为单位； 　　利用响应矩阵法来建立线性规划水资源管理模型，有如下水位与降深关系式[3]： 　　式中： h（i，n）为i观测点n时段末水位；h0i为i观测点初始水位；S（i，n）为i观测点n时段末的水位降深值。 　　以追求地下水系统最大开采量的管理模型为例，目标函数为规划期内L个规划阶段N口抽水井的地下水总开采量最大： 　　目标函数： 　　约束条件包括每口井各规划阶段出水能力限制、地下水总开采量指标约束、水位降深极限约束、地下水流状态约束和非负性约束，其中地下水流状态约束以式（6）的矩阵形式来表述： 　　通过式（6）、式（7）将所有的降深约束转换为水量约束，去掉了水流方程约束，减少了变量和约束条件数量。 　　2潜水含水层两单元稳定抽水量最优分配实例建模对比 　　2.1实例模型条件简述 　　如图1所示模型1，某地有一近似矩形的均质含水层，三面为隔水边界，一面傍河，河水位视为定水头，以此为基准面，地下水接受大气降水补给，向河流排泄，已知含水层几何尺