空间热网的拓扑建构与非对称结构下的水力工况计算研究.pdfVIP

空间热网的拓扑建构 与非对称结构下的水力工况计算 哈尔滨工业大学 周志刚 邹平华 摘要：在实际工程中受多种因素的影响，热网往往呈空间非对称性，此时若仍按照传统的 平面网方法计算，必然会产生较大的计算误差，甚至出现无法计算的情况。本文在介绍空 间热网拓扑建构与水力工况计算数学模型的基础上，对实际工程中比较典型的空间非对称 结构热网的水力工况计算方法进行了详细的说明，为实际工程热网的水力计算与分析提供 了理论参考。 关键词：空间热网；拓扑建构；非对称结构；水力工况计算 热网在最初设计时一般按照供水管网与回水管网镜像对称考虑，因此传统的热网水力 工况计算是在人为设定供水管网与回水管网镜像对称的前提下，将热网系统简化为平面网 络进行计算。而在实际工程中，受系统结构、敷设条件以及其他因素的影响，大部分热网 往往呈空间非对称性，此时若仍将管网简化为平面镜像对称，采用平面网方法计算，必然 会产生较大的计算误差，甚至出现无法计算的情况。此时若要准确模拟热网运行状态，必 须在非对称空间热网拓扑建构的基础上，采用空间网方法进行水力工况计算。 1 空间网水力工况数学模型 对于空间供热管网，设其节点数为n ，管段数为b ，则其空间网水力工况计算数学模 型为 ⎧ A Q = Q ⎪ k pk dk B ΔH = 0 ⎪ k pk ⎪ ⎨ (1) ⎪ ΔH pk S k Q pk Q pk =− DH k ⎪ DH A =+ A Q + A Q Q ⎪ k 0 1 pk 2 pk pk ⎩ 其中 Ak 为空间管网的基本关联矩阵，Qpk 为空间管网的管段流量列向量，Qdk 为 空间管网的节点入流流量列向量，B 为空间管网的基本回路矩阵，ΔH 为空间管网管 k pk 段水头损失列向量，S 为管段阻力系数对角阵，Q 为管段流量绝对值对角阵，DH 为 k pk k 管段中的水泵扬程列向量，A0 为系数a0 的对角阵，A1 为系数a1 的对角阵，A2 为系数 a2 的对角阵。 式(1)所示的空间热网水力工况计算矩阵方程组是非线性的。方程组中管段的结构是 已知的，即管段的阻力系数s 已知，而管段质量流量和管段压降为待求的量，共有2b 个 mk,i 未知数，基本关联矩阵A 为 n − 1 × b 阶矩阵，基本回路矩阵B 为(b − n + 1)× b 阶矩阵， k ( ) k 方程的个数也为 2b ，因此方程组有唯一解。上述热网水力工况模拟矩阵方程组可采用 ［ ］ MKP 求解，该计算方法简单、实用，收敛速度快 1 。 2 热网拓扑建构 从数学模型的表达形式来看，空间热网的水力工况计算方法与平面网方法类似。两者 最本质的区别体现在管网拓扑建构上。下面分别介绍平面热网与空间网的拓扑建构方法， 并分析两者的不同。 2.1 平面热网拓扑建构