多热源环状管网定压点个数数学推论.docVIP

多热源环状管网定压点个数数学推论 马琳1, 2官燕玲1 1 长安大学环境科学与工程学院 2 中国建筑西北设计研究院有限公司摘要：本文针对多热源环状管网允许设几个定压点的问题，利用图论的方法，将管网物理模型转化为数学模型，并利用线性代数的方法对该数学模型进行理论推导，从而得到多热源环状管网只允许设一个定压点的结论。 关键词：多热源环状管网定压点个数图论线性代数随着城市集中供热的快速发展，多热源环状管网将成为城市供热发展的必然趋势。多热源环状管网的供热以其节能、高效以及供热保障的优点备受关注。但由于管网水力工况的复杂性，有些内在的关联性和特性尚不清楚，因此发展受到一定的限制。环状管网定压点个数的问题是确保环状管网安全运行的关键，该问题业内持有不同的认识[1、2]。本文将利用图论的方法将管网模型转化成数学模型，并利用线性代数的方法对该数学模型进行理论推导，从而确定多热源环状管网允许定压点的个数。1 环状管网数学模型的建立 利用节点流量的概念，可以将具体的管网简化为只包含管段和节点两类元素的管网模型。管段是管网中流量和管径不发生变化的一段管道；节点是管段的端点，也是一些管段的交点。如果只考虑管段和节点之间的关联关系，管网即被抽象为管网图。与管段或节点有关的其它属性参数作为“权”值，赋予管段和节点，以便进行计算和研究。1.1 流体输配管网的矩阵表示 在流体输配管网图中，各管段的交汇点称为节点，各交汇点间的管段称为分支，由节点和分支构成的集合称为图，用表示。为了方便进行管网工况分析和水力计算，提出用矩阵代数来表示管网结构。 假设管网有J个节点，N条分支，可以用矩阵来表示分支和节点之间的关系，为图的关联矩阵，定义为[3]： 在关联矩阵中，每一行表示一个节点，每一列表示一个管段。的特点是：一列中总有一个数是1，一个数是 -1，其余数字都是0。是由0，1，-1组成的矩阵。关联矩阵的秩为J-1，去掉节点k所对应的一行，得到基本关联矩阵。 以图1为例，其中表示管段，表示节点，其关联矩阵为： 其基本关联矩阵为： 图1 流体管网图 1.2 环状管网的特性方程组 环状管网应满足基尔霍夫第一、第二定律，即节点流量平衡方程、管段能量平衡方程，表现如下： 节点流量平衡方程 其中为J行N列的关联矩阵，表示环网中J个节点与N个管段之间的拓扑关系，它的秩为J-1， 已知量；N维列向量，表示N个管段的流量，未知量；J维列向量，表示J个节点的流量，已知量。该方程的物理意义为：与任何节点关联的所有分支流量其代数和等于该节点的流量。管段能量平衡方程 其中为的转置矩阵；J维列向量，表示J个节点的压力，未知量； N维列向量，表示N个管段的压力降，未知量。该方程的物理意义为：管段两端节点水头之差等于该管段的压降。 分枝阻力定律及泵的扬程曲线代数方程 其中是式3）中的第i个分量是式2）中的第i个分量，i=1，…，N。因为流动能量损失公式为，为了表示与方向的一致性，故表示为。当管段中有动力装置泵时，为泵的扬程曲线函数关系式，通常表现为一元二次式，其中表示泵所在管段的管段流量，的值为该泵的扬程。该管段的压力降应该等于阻力损失与泵的扬程的差值：，即，和分别为第i段管段的压力降和流量，二者为一一映射关系。N个管段的压力降和流量均满足式，管段中有动力装置时，取；管段中没有动力装置时，取0。 目前的已知量为 ，， ，未知量为，，。环状管网的水力计算就是通过已知的基本关联矩阵和所有节点的节点流量，求出所有管段的压力损失、管段流量以及所有节点的节点压力。借助以下的数学表示法来分析环状管网单点、多点定压问题。1.3 环状管网的数学模型 环状管网应满足节点流量平衡方程、管段能量平衡方程、分枝阻力定律及泵的扬程曲线代数方程。通过以上转化，可以用如下方程组将环状管网表示出来，将式(2)、(3)、(4)组成的方程组称为方程组。 一个拥有J个节点，N个管段的环状管网系统，其数学表示如下： 其中为J行N列的关联矩阵，它的秩为J-1，已知量；N维列向量，未知量；J维列向量，已知量；J维列向量，未知量为N维列向量，未知量；泵的扬程-流量函数关系式，是的二次函数，函数关系已知可以表示为其中，式中的是式中的第i个分量，是式中的第i个分量，i=1，…，N。2 确定环状管网定压点个数的数学推导过程 首先在本节将引入一个数学概念：矩阵A的所有非零阶子式的最高阶数称为A的秩，记为并规定零矩阵的秩等于0[4]。将式、、组成的方程组称为方程组I）。本小节需要对方程组I）进行一系列的推导及分析，证明环状管网只能进行单点定压，其具体推导过程如下。