多源大型水管网的优化调度论文.docVIP

　　多源大型水管网的优化调度论文 摘要：针对多水源的大型供水管网，运用两级递阶优化原理，建立了反映管网各种水力参数的数值仿真模型，运用该数值模型可达如下目的：对不同水源的供水量进行优化配置，并且对某一水源地泵站内水泵的运行进行优化调度，合理利用能源，以取得良好的经济效益和社会效益；通过计算对比，改变旧有管网的管路参数或增删管网内局部支路.freel计，反映了供水的服务质量； 节点流量(供水量或用水量)，单位以m3/s或1/s计； 管段流量，单位以m3／s或1／s计；与管段两端压差呈非线性关系。 约束——约束分为等式约束条件和不等式约束条件。 等式约束条件：据物质守恒定律，管网中每个节点i的流量平衡满足 图1 供水系统两级递阶优化原理 Qi=∑qij (1) 即从节点i流出的流量等于与该节点相连的各管段ij流入流量之和；i为节点号，qij为i节点与j节点之间管段的流量，j＝1，2，…，n，n为与节点相连的管段数。 不等式约束条件： (2) 即泵站i的供水量有上下限，式中U为水源供水量。 b：min｛Fp｝≥Cp (3) 即要求管网中所有节点的最小自由水头必须大于等于某一常数C，以满足服务质量要求。｛Fp｝为管网节点自由水头的集合。 目标函数J——供水成本(水费+电费，以元计)的表达式为 i=1,2,…,m (4) 式(4)中V，为泵站号。 供水系统优化供水量的目的是在满足约束条件等式(1)、(2)和不等式(3)、(4)的前提下，使式(4)的J值，亦即供水成本最小。这属于有约束条件的非线性规划问题。 解决供水优化的计算分两步：第一步是管网平差(水力分析)，求出系统状态变量如各管段流量、各节点压力等。第二步是将水源地供水量(Ui)和扬程(Pi/γ-Zi)代入式(4)，求出目标函数当前值J(k)。按库恩-塔克(Kuhn-Tucker)的最优必要条件，用拉格朗日乘子法求得供水量增量ΔUi,得到水源下一步的供水量。 在迭代过程中，将前后两次求得的J值进行比较，如果J的绝对值不断下降则继续迭代，直到|Jk+1-Jk≤ε1时停止迭代计算。ε1为指定误差，本文选ε＝0.005元。这时的水源供水量U(k+1)i和(Pi/γ-Zi)(k+1)为所求。 1.1 平差——求各管流量、各节点压力 依据公式 Qi=∑qij，每个节点i流量应等于与之所连各管段流量之和qij。 管段ij沿程能量损失公式： (5) 从上式可得： |qij|=Sij|hij|1/2 (6) 其中：Sij为管段ij的水力特性常数。 平差步骤： ① 进行管网平差之前应预先假设各水源供水量，再按式(1)分配各管段流量； ② 从式(5)求得各管段压降hij；可应用线性理论引入一个变量；Cij=|qij|／|hij|将上式代入每个节点的流量平衡式：Qi=∑qij=∑Cij(Pi-Pj)。如果节点1与节点2，3，4有相连管段，节点2与节点3，5，6有相连管段，节点N-1与节点N，N-2，N-3有相连管段，可写出： 化为： ③如果该系统有N个节点，可以组成N-1个方程的线性方程组：C·P=Q，对于Pi,i=1,2,…,n，选择一个已知最小自由水头Hzk的节点k，则Pk/γ=Zk+Hzk。式中Zk为节点K的高程。［C］可求，［Q］已知，用高斯消去法解此方程组，可求得各管段压差； ④通过式(6)求出各管段的qij值，并代入式(5)，重复步骤②、③，经多次迭代，直到各管都满足时迭代才结束。ε2为各管流量迭代允许误差。 1.2 供水优化 为了找到供水成本J的最小值，应用拉格朗日乘子法将目标函数扩展为： i=1,2,…,m，m为水源总数。其中，L为向量q、p、λ的函数，也称拉格朗日函数，［λT］为拉格朗日乘子向量的转置。 因为根据库恩塔克条件，在满足下式： i=1,2,.freel (8) i=1,2,…,n (9) 时，L取得最小值，即J的最小值，n为节点总数，m为水源总数。Lmin(q,p,λ)的最优解，也是Jmin(q,p,λ)的最优解。 因为式(9)，该式通过平差已得到满足。将平差结果qij代入式(7)，有，可求出向量λ，把λ代入式(8)，得到的值，,是L对U的梯度向量。按梯度法原理，在某一小范围内取负梯度方向能使函数值下降最快3。由此可求各水源供水量的下一步迭代值：。其中C为控制常数，为加快计算收敛，可以先选较大值，快接近最优解时选小些的值。 为满足不等式约束条件(2)，应保持在和之间。迭代若干次，将每次求得的J值进行比较，直到|J(k+1)-J(k)|≤ε1成立并且满足如下条件时可认为为所求最优解。 当 ΔuiL=0； (10) ΔuiL≤0； ΔuiL≥0。 这时各水源节点供水量为目标函数最优解，也就是泵站应供水量。所得各水