热水网路水力工况的计算及其图形显示

提要 　热水网路水力工况计算是研究在管网结构已定，运行条件变化时计算热用户的流量分配和热用户的水力失高度。研究偏离设计工况时管网失调的程度以及在采用各种措施后调整热用户失调度的能力。本文给出了研究热水网路水力工况计算的数学模型和求解方法。可用所编制的程序计算和研究简单技术热水网路在不调节循环水泵出口阀门时其性能曲线与网路曲线的自然交汇点；得出循环水泵扬程变化时出口阀门节流、供、回水干管阀门节流、干管或热用户阀门开度增大或减小、干管或用户堵塞、供、回水干线设中继泵、系统或用户设混水泵时对整个热水网路影响和各热用户的流量变化规律；可计算出水力失调度并显示其水压图。

关键词 　热水网路水力工况 阻力数 水泵特性曲线 水力失调度 1　热水网路水力工况分析与计算的数学模型

设计热水网路时是用已知的用户热负荷去确定各管段的管径、阻力损失以及网路的总阻力损失，选择循环水泵的扬程。分析和计算热水网路的水力工况时正好相反，是对已经设计完毕的或需要改扩建的热网，在已知循环水泵的型号以及各管段的管径时，来确定各管段和热用户的流量。将水泵和网路的特性方程联立求解可以定量和定性解决这一问题。

1．1水泵的特性曲线拟合方程

水泵为网路提高能量，是热媒循环的动力。大型网路中可能有循环水泵、中继泵、加压泵等多组水泵。需对其流量-扬程曲线进行拟合，一般可用下式表示：

Hp=f(G)

(1)

式中：Hp--水泵扬程

f(G)--拟合得到的水泵性能特性曲线公式

本文采用最小二乘拟合水泵特性曲线曲线。该方法可使拟合误差达到最小值，并且该解析式给用矩阵方程分析网路水力工况分析提供了基本条件。大多数离心泵的G-H关系曲线如图1所示，若图中1、2两占之间的曲线为水泵的高效段，可用下式来近似描绘：

图1 水泵G－H曲线

Hp= Hx -SxG2

(2)

式中：Hp--水泵的虚总扬程，mH2O;

Sx --水泵的虚阻耗系数，s2/m2;

G--水泵的总流量，m3/S。

对点1、2可写出

（3）

（4）

求出Sx、Hx，式（2）即被确定。按这种方式确定的解析式，其近似性较差。还须在水泵G-H曲线上取多组数据（G1，H1）、（G，H2）……（Gx，Hx），根据最小二乘原理来确定式（2）中的Sx与Hx。由于在研究水力工况时，流量是未知的，而且在非设计工况下去选择热网也不一定工作在高效段，所以所取数据应涵盖其整个工作区。采用最小二乘原理的Sx与Hx计算式如下：

（5）

例如选取型号为12sh-6A的水泵，转速n=1450转/min，其特性曲线如图2所示。在特性曲线工作段内取13组数据，根据式（5）与式（6）可求出Hx =96.3mH2O、Sx=406.1s2/m5，因此该水泵的特性曲线方程为：H=96.3-406.1(G/3600)对采用多泵的复杂管网而言，可写出如下矩阵方程：

（6）

式中：Hp--水泵扬程矩阵。Hpi为管段i的水泵扬程，对无水泵管段Hpi=0。

图2 　水泵12sh-6A特性曲线

1．2 描绘管路的方程

热媒在管网中流动时产生阻力，消耗能量。当其长度局部构件一定时，其特性可用以下三组方程来描述。

1．2．1节点方程

节点方程就是节点流量连续性方程，即连接于任何节点的所有管段流量的代数和为零。可用以下矩阵表示：

AG=Q

（7）

式中：A--管网图的基本关联矩阵；

G--管段的流量矩阵，G=（G1，G2，…Gi…，GN）T，Gi为管段i的流量；

Q--节点的流量矩阵，Q =（Q1，Q2，…Qi…，QN）T，Qi为管段i的流量，本文中取Q=0

1．2．2 回路方程

回路方程就是能量方程或环方程，即每个环的水头损失闭合差为零，写成矩阵的形式即为

BΔH=0

（8）

式中：B--管网图的基本回路矩阵；

ΔH--管段阻力损失，ΔH =（ΔH1，ΔH2，…ΔHi…，ΔHN）T，ΔHi为管段i的阻力损失；

0--0向量，即0=（0，0，……，0）T。

1.2.3阻力方程

管段流量G与阻力损失ΔH之间的关系可用下式表示：

ΔH=SG2-Hp

（9）

式中：S--阻力数，它与管材、管长、管径以及产生局部阻力损失的管路附件有关；

G--管段的流量

Hp--水泵扬程

当管段中有水泵时，水泵作为一个负阻力损失，管路无水泵时Hp=0，将各管段应用上式并写成矩阵形式：

（10）

2 热水网路水力工况计算与分析的算例

2．1用矩阵方程组求水泵与网路自然交汇工作点

如图3为有六个热用户的供热系统，其管段b和节点n编号如图3所示。分支节点编号为1,2,3，…,11；管段编号为b1，b2，……b16，其中管段b1由两部分组成，即br+ b1（br代表热源内管段，b1′代表热源出口到节点1管段）；相应的各管面流量纺号为Gb1，Gb2，……Gb16，各管段的压降编号为ΔHb1，ΔHb2，……，ΔHb16，其中管段b1的压降由两部分组成，即ΔHr+ΔHb1′；系统循环水泵的扬程为Hp，管段阻力系数编号为Sb1，Sb2，…Sb16。各管段的水流方向如图3中所示。此热水网路的设计数据如下：循环；型号为12sh-6A，设计工况运行时各热用户的流量为100m3/h，各热用户的资用压头为10mH2O，各段供、回水干管管段的阻力损失为5 mH2O。其总阻力损失为80 mH2O，由此可得出各管段的阻力数。

图3 热水网路编码示意图

利用上面所述的公式可列出图3中热水网路系统水力工况数学模型。其中数学模型中流量矩阵G及管段阻力数矩阵S为对角阵，关联矩阵A=（aij）、基本回路矩阵B=（bkj）如下：