032 供热管网阻力系数辨识方法比较的研究.docVIP

供热管网阻力系数辨识方法比较研究 哈尔滨工业大学 刘永鑫 邹平华 摘要 近年来，供热管网发展迅速，因此对水力工况计算、运行调节及故障诊断等提出了更高的要求。而这些研究能否得到准确的结果一般与供热管网阻力系数的准确程度有关，以往的研究都是在应用供热管网阻力系数设计值的基础上进行的，这种条件下有时会使水力工况计算、运行调节及故障诊断等研究的准确性降低。因此如果能够获得管网实际运行的阻力系数值，就会使供热管网水力工况计算等理论的实用性将得到很大提升。 本文在分析管网阻力系数辨识方程组的基础上，分别对三种辨识方法进行分析比较，并将其应用于算例供热管网，结果表明：三种辨识方法均可以应用于实际供热管网阻力系数辨识问题，根据实际情况选择更合适的辨识方法可以得到更准确的辨识结果。 关键词 供热管网，管网阻力系数，辨识方法 1 引言 供热管网的发展趋势是由小型化向大型化发展，这对供热管网的水力工况计算、运行调节及故障诊断等提出了更高的要求。其中，供热管网阻力系数的准确程度直接影响管网水力工况计算结果的准确性，同时也关系到调节和调度方法，以及发生故障时采取应急措施的针对性。因此如果能够通过参数辨识的方法获得管网实际运行的阻力系数值，那么已有的供热管网水力工况计算和运行调节等理论的实用性将得到很大提升。 有关供热管网阻力系数辨识研究较少，Liggett和Chen首先提出拟稳态条件下的给水管网水力工况的校正模型，但其目的是为了检测和量化漏失[1,2]。Kapelan等基于最小二乘原理建立测量值与辨识值之差平方和最小为目标函数的数学模型，并利用遗传算法解决该管网阻力系数辨识问题[3]。刘永鑫在该方法基础上提出了一种测点优化布置方案，在一定程度上提高了辨识的精度[4]。 供热管网阻力系数研究方面，秦绪忠等人提出一种供热管网阻力系数辨识方法——RC2法[5]，但是该方法对供热管网水力工况变化要求较高，当工况变化不大时，易发生不收敛的情况。周志刚提出在某一确定工况条件下根据一次观测得到的测量信息辨识得到全网的管段阻力系数的方法，在一定测点布置方案下，该方法可保证管网阻力系数具有唯一解[6]。但该方法要求在一些管段布置流量测点，因此实际工程应用有一定困难。文献[6]中还将优化算法应用于供热管网阻力系数辨识，得到较好的辨识结果但该方法在解决管网阻力系数辨识问题时。为解决管网阻力系数辨识中有时存在辨识方程组是不确定方程组，即辨识不满足能观性的问题刘永鑫提出基于方程组广义逆的供热管网阻力系数辨识方法和基于部分可信管段阻力系数初值的供热管网阻力系数辨识方法，这两种辨识方法在一定程度上满足能观性条件，可以得到唯一的管网阻力系数辨识结果[7]。 基本方程 供热系统运行时，常规供热管网水力工况计算和阻力系数辨识都满足一系列类似方程所组成的几个矩阵方程组，目前应用矩阵方程已经可以很好地描述并解决供热管网水力工况计算的相关问题，但对于供热管网阻力系数辨识来说，相关研究还不很完善。因此，对应用于供热管网阻力系数辨识基本方程进行组合以进一步提高辨识精度和计算效率。.1 供热管网水力工况计算及阻力系数辨识方程 对于具有节点数为n+1，管段数为b的供热管网，当应用图论于进行水力工况计算时，供热管网中的热媒流动可以用下述矩阵方程表示。节点方程即连续性方程为 (1) 式中 A ——基本关联矩阵，A=(aij)n×b，其中元素aij的取值规定如下 aij= G ——管段流量列向量，m3/s； Q ——节点泄流量列向量，m3/s。 另外，考虑管网内流动介质的能量方程，即管段两端压力差表示为 (2) 式中 P ——节点压力列向量，P=(pi)n×1 ，其中pi为第i个节点上相对于参考节点的节点压力，m； S ——管段阻力系数矩阵，b×b阶对角阵，其中sj为第j根管段的阻力系数，s2/m5； ——管段流量绝对值对角矩阵，b×b阶对角阵，其中gj为第j根管段的流量，m3/s； Z——管段平均标高列向量，m。 进行供热管网水力工况计算时，对于式(1)和(2)联立构成的非线性方程组来说，供热管网阻力系数矩阵S是已知的。此时，该方程组共有b+n个未知量(b个gj 和n个pi)，b+n个独立方程，因此，此时方程组有唯一解，可以唯一地描述供热管网的水力工况。 在进行供热管网阻力系数辨识计算时，可以利用式(1)和(2)联立构成的方程组求解时，供热管网阻力系数矩阵S是未知的，而部分节点压力和管段流量是可观测的。在已知节点压力和管段流量数量不同的条件下讨论该方程组中未知量和独立方程数量对于方程组的解分析十分重要。式(1)和(2)联立表示的单一水力工况的条件下的供热管网阻力系数辨识方程组中有n+ b个独立方程，在其点压力和管段流量均未知的条件下，其中有b个管段流量、n个节