CN108280302A - 一种热网热流计算方法 - Google Patents
一种热网热流计算方法 Download PDFInfo
- Publication number
- CN108280302A CN108280302A CN201810079288.4A CN201810079288A CN108280302A CN 108280302 A CN108280302 A CN 108280302A CN 201810079288 A CN201810079288 A CN 201810079288A CN 108280302 A CN108280302 A CN 108280302A
- Authority
- CN
- China
- Prior art keywords
- heat
- water
- pipeline
- supply network
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000012530 fluid Substances 0.000 title abstract 5
- 238000000205 computational method Methods 0.000 title abstract 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 485
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000012546 transfer Methods 0.000 claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims description 24
- 238000003860 storage Methods 0.000 claims description 6
- 238000010438 heat treatment Methods 0.000 abstract 4
- 238000012797 qualification Methods 0.000 abstract 1
- 238000004519 manufacturing process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005485 electric heating Methods 0.000 description 2
- 238000012887 quadratic function Methods 0.000 description 2
- 239000008400 supply water Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Heat-Pump Type And Storage Water Heaters (AREA)
Abstract
本发明提供一种热网热流计算方法，包括：基于热负荷模型，应用各热负荷的消耗量和各管道的入口温度最小值计算各热负荷的循环水质量流量；基于热源模型，应用各热源的供热量和各管道的入口温度最小值计算各热源的循环水质量流量；其中，各热源的供热量根据供热量最小值对各热源进行分配；供热量最小值为各热负荷的消耗量与热网热损耗最小值之和；热网热损耗最小值根据各管道的入口温度最小值计算获取；基于热网流量模型，根据各热负荷和热源的循环水质量流量，计算各管道的循环水质量流量；将各管道的循环水质量流量和热网传热模型作为限定条件，求解热流方程，将热流作为热网最优热流。本发明提供的方法，实现了最优热流模型的高效计算。
Description
技术领域
本发明涉及能源经济技术领域，尤其涉及一种热网热流计算方法。
背景技术
随着能源危机和环境污染问题的日益严重，通过多能源系统的互补和耦合实现能源利用效率的提高，受到了越来越多研究人员的关注。
其中，基于电力系统和热力系统的耦合运行，能够实现能量在电力系统和热力系统间的灵活转换和能量的梯级利用与开发，有利于提高整体能源利用效率。
在目前的电热耦合系统的研究中，大多数研究仅仅将热力系统建模为一个供热与用热平衡的系统，并未考虑到热网的网络拓扑、管道参数及热量损耗等等一系列的因素，难以全面地反映出热网的运行特性。而实际上，城市供热系统是以供回水管道为线路、以水为传热工质构建起来的能量网络模型。因此，对于热力系统的分析必须考虑详细的网络模型，获得能量流动、损耗、温度变化等具体信息，才能对热网进行安全经济的分析和调度。
为了进行热网的经济运行调度分析，必须根据热网的详细网络模型建立起相应的最优热流模型。然而，热网的最优热流模型是一个非凸非线性的高度复杂的优化模型，现有的研究中并没有提出有效的求解方法，难以投入到工程实践中，给电热耦合系统的运行决策带来极大的不便。
发明内容
本发明为解决现有技术中存在的最优热流模型难以求解的问题，提供了一种热网热流计算方法。
一方面，本发明提出一种热网热流计算方法，包括：S1，基于热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算各热负荷的循环水质量流量；基于热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算各热源的循环水质量流量；其中，所述各热源的供热量根据供热量最小值对各热源进行分配获取；所述供热量最小值为所述各热负荷的消耗量与热网热损耗最小值之和；所述热网热损耗最小值根据所述各供水管道和回水管道的入口温度最小值计算获取；S2，基于热网流量模型，根据所述各热负荷和各热源的循环水质量流量，计算各供水管道和回水管道的循环水质量流量；S3，将所述各供水管道和回水管道的循环水质量流量和热网传热模型作为限定条件，求解热流方程对应的热流，将所述热流作为热网最优热流。
优选地，所述步骤S1中，基于热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算各热负荷的循环水质量流量，进一步包括：将所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值代入下式：
式中，为热负荷k的消耗量，cp为水的比热容，为热负荷k的循环水质量流量，为热负荷k的供水温度，为热负荷k的回水温度，SHL为热负荷集合；分别将与热负荷k连接的供水管道和回水管道的入口温度最小值作为和
优选地，所述步骤S1中，基于热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算各热源的循环水质量流量，进一步包括：将所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值代入下式：
式中，为热源j的供热量，cp为水的比热容，为热源j的循环水质量流量，为热源j的供水温度，为热源j的回水温度，SHS为热源集合；分别将与热源j连接的供水管道和回水管道的入口温度最小值作为和
优选地，所述步骤S1中，热网热损耗最小值根据所述各供水管道和回水管道的入口温度最小值计算获取，进一步包括：将任一管道的入口温度最小值代入下式，获取该管道的热损耗：
式中，ΔQb为管道b的热损耗，λb为管道b的温度损耗系数，为管道b的入口温度，Lb为管道b的长度，Ta为环境温度，Spipe为供水管道和回水管道集合；
将各管道的热损耗之和作为热网热损耗最小值。
优选地，所述步骤S2中，所述热网流量模型包括流量连续性方程和压力降落平衡方程；所述流量连续性方程包括：
其中，和分别为供水管道a和回水管道d的循环水质量流量，和分别为热源j和热负荷k的循环水质量流量，和分别为以节点i为末端的供水管道集合和回水管道集合，和分别为节点i为首端的供水管道集合回水管道集合，和分别为连接节点i的热源集合和热负荷集合，Snode为热网节点集合；所述节点包括热源和热负荷；所述压力降落平衡方程包括：
式中，和分别为供水管道a和回水管道d的水压降落，和分别为供水管道a和回水管道d的压力损耗系数，和分别为供水管道a和回水管道d的循环水质量流量，为供水管道集合，为回水管道集合，为供水管道环路集合，为回管道环路集合。
优选地，所述步骤S3中，所述热网传热模型包括热负荷模型、热源模型、温度损耗方程和温度混合方程；所述温度损耗方程包括：
式中，和分别为管道b的出口温度和入口温度，Ta为环境温度，λb为管道b的温度损耗系数，Lb为管道b的长度，cp为水的比热容，为管道b的循环水质量流量，Spipe为供水管道和回水管道集合；所述温度混合方程包括：
式中，和分别为供水管道a和回水管道d的出口温度，和分别为供水管道a和回水管道d的入口温度，和分别为供水管道a和回水管道d的循环水质量流量,和分别为热源j的供水温度和回水温度，和分别为热负荷k的供水温度和回水温度，和分别为热源j和热负荷k的循环水质量流量，和分别为节点i的供水混合温度和回水混合温度，和分别为以节点i为末端的供水管道集合和回水管道集合，和分别为节点i为首端的供水管道集合回水管道集合，和分别为连接节点i的热源集合和热负荷集合，Snode为热网节点集合。
优选地，所述步骤S3中，当热源为燃煤产热机组时，所述热流方程如下：
min CCU
式中，
为热网成本，为热流，b0,i、b1,i和b2,i分别为成本系数的常数项、一次项和二次项，SCU为燃煤产热机组集合。
另一方面，本发明提出一种热网热流计算设备，包括：至少一个处理器；以及与所述处理器通信连接的至少一个存储器，其中：所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行如前所述的方法。
又一方面，本发明提出一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行如前所述的方法。
本发明提供的一种热网热流计算方法，应用固定的热网管道的循环水质量流量将热流方程约束退化为线性等式约束，简化了热网的最优热流模型，实现了最优热流模型的高效计算，有助于电网详细模型的推广和热网运行特性的全面反映。
附图说明
图1为本发明具体实施例的一种热网热流计算方法的流程示意图；
图2为本发明具体实施例的一种热网热流计算设备的结构示意图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
图1为本发明具体实施例的一种热网热流计算方法的流程示意图，如图1所示，一种热网热流计算方法，包括：S1，基于热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算各热负荷的循环水质量流量；基于热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算各热源的循环水质量流量；其中，所述各热源的供热量根据供热量最小值对各热源进行分配获取；所述供热量最小值为所述各热负荷的消耗量与热网热损耗最小值之和；所述热网热损耗最小值根据所述各供水管道和回水管道的入口温度最小值计算获取；S2，基于热网流量模型，根据所述各热负荷和各热源的循环水质量流量，计算各供水管道和回水管道的循环水质量流量；S3，将所述各供水管道和回水管道的循环水质量流量和热网传热模型作为限定条件，求解热流方程对应的热流，将所述热流作为热网最优热流。
具体地，热网模型是一个包含供回水管道、热源和热负荷等的复杂能量网络，根据其建模出来的最优热流模型是一个高度非凸非线性的优化模型，现有的商业优化软件无法进行可靠求解，难以投入到工程实践当中。本发明具体实施例中，求解热网管道的循环水质量流量，并应用固定的热网管道的循环水质量流量将热流方程约束退化为线性等式约束，以便于现有商业优化软件能够高效求解。其具体求解方法如下：
首先，基于热网模型中的热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算热网中各热负荷的循环水质量流量。此处，所述各供水管道和回水管道的入口温度最小值为所述热网中预先设定的各供水管道和回水管道的入口温度的下限。所述热负荷的循环水质量流量是指热负荷从供水管道流出到回水管道的循环水质量流量。
与此同时，基于热网模型中的热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算所述热网中各热源的循环水质量流量。此处，所述热源的循环水质量流量是指热源从回水管道注入供水管道的循环水质量流量。
其中，根据热网中各供水管道和回水管道的入口温度最小值计算可以获取所述热网的热网热损耗最小值。所述热网的热网热损耗最小值与预先设定的各热负荷的消耗量相加，得到所述热网的供热量最小值。随后，根据能量守恒定律，在不考虑热流方程约束的条件下，以热源生产成本为目标进行热源出力的经济分配，获取各热源的供热量，并将所述各热源的供热量应用于各热源循环水质量流量的求解，即本发明具体实施例中的各热源的供热量，是依据热网的整体热损耗最小原则获取的。
其次，将上一步骤中计算获取的热网中各热负荷的循环水质量流量和各热源的循环水质量流量代入热网模型中的热网流量模型，计算热网中各供水管道和回水管道的循环水质量流量。
最后，应用上一步骤中计算获取的热网中各供水管道和回水管道的循环水质量流量替代原最优热流模型中应用热网流量模型构成的限定条件，将固定的各供水管道和回水管道的循环水质量流量和热网传热模型作为限定条件，求解热流方程中对应的热流，并将上述计算过程中获取的热流，即热流方程中生产成本最小时对应的热流作为所述热网的最优热流。
本发明具体实施例中，应用固定的热网管道的循环水质量流量将热流方程约束退化为线性等式约束，简化了热网的最优热流模型，实现了最优热流模型的高效计算，有助于电网详细模型的推广和热网运行特性的全面反映。
基于上述具体实施例，一种热网热流计算方法，所述步骤S1中，基于热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算各热负荷的循环水质量流量，进一步包括：将所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值代入下式：
式中，为热负荷k的消耗量，cp为水的比热容，为热负荷k的循环水质量流量，为热负荷k的供水温度，为热负荷k的回水温度，SHL为热负荷集合。
具体地，在热网中的各热负荷节点，循环水从供水管道经过换热器到达回水管道，将热量释放给各热负荷节点。上式所示的热负荷模型表示了热负荷的消耗量与热网中热负荷的循环水质量流量以及供水温度和回水温度的关系。
在应用各热负荷的消耗量和各供水管道和回水管道的入口温度最小值对各热负荷的循环水质量流量进行计算时，将与任一热负荷连接的回水管道的入口温度最小值作为式中的该热负荷的回水温度，将与任一热负荷连接的供水管道的入口温度最小值作为式中的该热负荷的供水温度。
由此计算获取的各热负荷的循环水质量流量是基于热网的整体热损耗最小原则的热负荷循环水质量流量。
本发明具体实施例中，通过对热网中各热负荷的循环水质量流量的计算，为后续获取整体热损耗最小原则下的固定的管道循环水质量流量提供了条件。
基于上述任一具体实施例，一种热网热流计算方法，所述步骤S1中，基于热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算各热源的循环水质量流量，进一步包括：将所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值代入下式：
式中，为热源j的供热量，cp为水的比热容，为热源j的循环水质量流量，为热源j的供水温度，为热源j的回水温度，SHS为热源集合。
具体地，在热网中，热源通过换热器加热循环水，为热网提供热量。上式所示的热源模型表示了热源的供热量与热网中热源的循环水质量流量以及供水温度和回水温度的关系。
在应用各热源的供热量和各供水管道和回水管道的入口温度最小值对各热源的循环水质量流量进行计算时，将与任一热源连接的回水管道的入口温度最小值作为式中的该热源的回水温度，将与任一热源连接的供水管道的入口温度最小值作为式中的该热源的供水温度。
由此计算获取的各热源的循环水质量流量是基于热网的整体热损耗最小原则的热源循环水质量流量。
本发明具体实施例中，通过对热网中各热源的循环水质量流量的计算，为后续获取整体热损耗最小原则下的固定的管道循环水质量流量提供了条件。
基于上述任一具体实施例，一种热网热流计算方法，所述步骤S1中，热网热损耗最小值根据所述各供水管道和回水管道的入口温度最小值计算获取，进一步包括：将任一管道的入口温度最小值代入下式，获取该管道的热损耗：
式中，ΔQb为管道b的热损耗，λb为管道b的温度损耗系数，为管道b的入口温度，Lb为管道b的长度，Ta为环境温度，Spipe为供水管道和回水管道集合；将各管道的热损耗之和作为热网热损耗最小值。
具体地，在影响管道热损耗的相关因素中，温度损耗系数、管道长度和环境温度均为常数，仅管道入口温度为可变变量。由此可知，管道入口温度越低，管道热损耗越小，热网的整体热损耗也越小。
管道热损耗与管道入口温度的关系式推导过程如下：
根据能量守恒定律，热网中各热源的供热量等于各热负荷的消耗量和热损耗之和。对于热网中的任一管道，其热损耗反映在其出口温度和入口温度之差，其具体关系描述如下：
式中，ΔQb为管道b的热损耗，cp为水的比热容，为管道b的循环水质量流量，和分别为管道b的出口温度和入口温度，Spipe为供水管道和回水管道集合。
随后，将温度损耗方程代入上式，所述温度损耗方程如下：
式中，Ta为环境温度，Lb为管道b的长度，Ta为环境温度。
由上得到如下公式：
通常，数值小于0.1。因而，对上式中的指数想进行泰勒展开并忽略高阶项：
本发明具体实施例中，验证了管道热损耗仅与管道入口温度有关，并通过对管道热损耗的计算，获取了通过对热网中各热源的循环水质量流量的计算，为后续获取整体热损耗最小原则下的固定的管道循环水质量流量提供了条件。
基于上述任一具体实施例，一种热网热流计算方法，所述步骤S2中，所述热网流量模型包括流量连续性方程和压力降落平衡方程；所述流量连续性方程包括：
其中，和分别为供水管道a和回水管道d的循环水质量流量，和分别为热源j和热负荷k的循环水质量流量，和分别为以节点i为末端的供水管道集合和回水管道集合，和分别为节点i为首端的供水管道集合回水管道集合，和分别为连接节点i的热源集合和热负荷集合，Snode为热网节点集合；所述节点包括热源和热负荷。
具体地，所述流量连续性方程是基于质量守恒定律建立的。对于热网中的任一节点，注入该节点的循环水质量流量之和与流出该节点的循环水质量流量之和相等。
由上述公式可知，热网中以各节点为末端的各供水管道的循环水质量流量与各热源从回水管道注入供水管道的循环水质量流量之和等于以各节点为首端的供水管道的循环水质量流量与各热负荷从供水管道流出到回水管道的循环水质量流量之和；
热网中以各节点为末端的各回水管道的循环水质量流量与各热负荷从供水管道流出到回水管道的循环水质量流量之和等于以各节点为首端的回水管道的循环水质量流量与各热源从回水管道注入供水管道的循环水质量流量之和。
所述压力降落平衡方程包括：
式中，和分别为供水管道a和回水管道d的水压降落，和分别为供水管道a和回水管道d的压力损耗系数，和分别为供水管道a和回水管道d的循环水质量流量，为供水管道集合，为回水管道集合，为供水管道环路集合，为回管道环路集合。
具体地，由于管道摩擦的存在，循环水压在管道中发生变化，所述压力降落平衡方程分别表示了供水管道的水压降落和回水管道的水压降落，以及闭合环路中供水管道的压降平衡方程和回水管道的压降平衡方程。
其中，闭合环路中，为了满足压降平衡，其环路水压降落之和为0。
本发明具体实施例中，通过热网流量模型中的流量连续性方程和压力降落平衡方程获取了各供水管道和回水管道的循环水质量流量，为将热流方程约束退化为线性等式约束，简化热网的最优热流模型提供了条件。
基于上述任一具体实施例，一种热网热流计算方法，所述步骤S3中，所述热网传热模型包括热负荷模型、热源模型、温度损耗方程和温度混合方程；所述温度损耗方程包括：
式中，和分别为管道b的出口温度和入口温度，Ta为环境温度，λb为管道b的温度损耗系数，Lb为管道b的长度，cp为水的比热容，为管道b的循环水质量流量，Spipe为供水管道和回水管道集合。
具体地，由于管道与环境的温度差等因素，循环水在管道的流动过程中会产生热损耗，从而导致温度的变化。所述温度损耗方程用于表示管道入口温度与管道出口温度之间的关系。
所述温度混合方程包括：
式中，和分别为供水管道a和回水管道d的出口温度，和分别为供水管道a和回水管道d的入口温度，和分别为供水管道a和回水管道d的循环水质量流量，和分别为热源j的供水温度和回水温度，和分别为热负荷k的供水温度和回水温度，和分别为热源j和热负荷k的循环水质量流量，和分别为节点i的供水混合温度和回水混合温度，和分别为以节点i为末端的供水管道集合和回水管道集合，和分别为节点i为首端的供水管道集合回水管道集合，和分别为连接节点i的热源集合和热负荷集合，Snode为热网节点集合。
具体地，所述温度混合方程是根据能量守恒定律建立的。上述四个方程中，后两个方程表示从任一节点流程的循环水温度等于该节点的混合温度。
基于上述任一具体实施例，一种热网热流计算方法，所述步骤S3中，当热源为燃煤产热机组时，所述热流方程如下：
min CCU
式中，
其中，为热网成本，为热流，b0,i、b1,i和b2,i分别为成本系数的常数项、一次项和二次项，SCU为燃煤产热机组集合。
具体地，燃煤产热机组是热网中常见的热源，其生产成本可由自变量为热流的二次函数表示。
根据上述基于热网成本最小化的热流方程，应用上一步骤中求取的各供水管道和回水管道的循环水质量流量和热网传热模型作为限定条件，构建最优热流模型：
min CCU
s.t.热网传热模型
和为定值
式中，和分别为供水管道a和回水管道d的循环水质量流量。
在和为定值的情况下，热网传热模型约束为线性约束，同时目标函数为二次形式，则上式表示的最优热流模型为一个线性约束下的凸优化模型，利用商业优化软件可进行高效求解，得到各热源出力、管道温度分布等最优热流结果。
为了更好地理解与应用本发明提出的一种热网热流计算方法，本发明进行以下示例，且本发明不仅局限于以下示例。
首先，建立热网的最优热流模型：
本示例中，热网模型包括供回水管道、热源和热负荷。水作为热网中的传热工质在供水管道和回水管道中循环流动，热源生产热量通过换热器加热循环水，同时热负荷通过换热器从循环水中获取热量。
热网模型分为两个部分：热网流量模型和热网传热模型。热网流量模型描述了热网中的管道循环水的质量流量、管道压力损耗的情况；热网传热模型描述了热网中热量传递过程和管道温度变化的情况。
其中，所述热网流量模型包括流量连续性方程和压力降落平衡方程：
所述流量连续性方程是基于质量守恒定律建立的。对于热网中的任一节点，注入该节点的循环水质量流量之和与流出该节点的循环水质量流量之和相等。所述流量连续性方程如下：
其中，和分别为供水管道a和回水管道d的循环水质量流量，和分别为热源j和热负荷k的循环水质量流量，和分别为以节点i为末端的供水管道集合和回水管道集合，和分别为节点i为首端的供水管道集合回水管道集合，和分别为连接节点i的热源集合和热负荷集合，Snode为热网节点集合；所述节点包括热源和热负荷。
由于管道摩擦的存在，循环水压在管道中发生变化，所述压力降落平衡方程分别表示了供水管道的水压降落和回水管道的水压降落，以及闭合环路中供水管道的压降平衡方程和回水管道的压降平衡方程：
式中，和分别为供水管道a和回水管道d的水压降落，和分别为供水管道a和回水管道d的压力损耗系数，和分别为供水管道a和回水管道d的循环水质量流量，为供水管道集合，为回水管道集合，为供水管道环路集合，为回管道环路集合。
所述热网传热模型包括热负荷模型、热源模型、温度损耗方程和温度混合方程：
热负荷模型用于表征热负荷的消耗量与热网中热负荷的循环水质量流量以及供水温度和回水温度的关系：
式中，为热负荷k的消耗量，cp为水的比热容，为热负荷k的循环水质量流量，为热负荷k的供水温度，为热负荷k的回水温度，SHL为热负荷集合。
热源模型用于表征热源的供热量与热网中热源的循环水质量流量以及供水温度和回水温度的关系：
式中，为热源j的供热量，cp为水的比热容，为热源j的循环水质量流量，为热源j的供水温度，为热源j的回水温度，SHS为热源集合。
所述温度损耗方程用于表示管道入口温度与管道出口温度之间的关系：
式中，和分别为管道b的出口温度和入口温度，Ta为环境温度，λb为管道b的温度损耗系数，Lb为管道b的长度，cp为水的比热容，为管道b的循环水质量流量，Spipe为供水管道和回水管道集合。
所述温度混合方程用于在不同管道的温度不同的循环水汇合在某一节点时，节点混合温度的计算：
式中，和分别为供水管道a和回水管道d的出口温度，和分别为供水管道a和回水管道d的入口温度，和分别为供水管道a和回水管道d的循环水质量流量，和分别为热源j的供水温度和回水温度，和分别为热负荷k的供水温度和回水温度，和分别为热源j和热负荷k的循环水质量流量，和分别为节点i的供水混合温度和回水混合温度，和分别为以节点i为末端的供水管道集合和回水管道集合，和分别为节点i为首端的供水管道集合回水管道集合，和 点集合。
基于上述热网流量模型和热网传热模型，建立仅含燃煤产热机组的热网最优热流模型如下：
min CCU
s.t.热网传热模型
热网流量模型
式中，
为热网成本，为热流，b0,i、b1,i和b2,i分别为成本系数的常数项、一次项和二次项，SCU为燃煤产热机组集合。
目前，上述最优热流模型的目标函数是二次函数形式，属于凸函数；若将热网传热模型和热网流量模型直接联合求解，由于模型中的等式约束的非凸性较强，特别是温度损耗方程中含有复杂的指数形式约束，给模型的求解造成较大的困难，目前的商业优化软件往往无法直接求得其最优热流的结果。
在此基础上，对热网热损耗进行估算：
根据能量守恒定律，热网中各热源的供热量等于各热负荷的消耗量和热损耗之和。对于热网中的任一管道，其热损耗反映在其出口温度和入口温度之差，其具体关系描述如下：
式中，ΔQb为管道b的热损耗，cp为水的比热容，为管道b的循环水质量流量，和分别为管道b的出口温度和入口温度，Spipe为供水管道和回水管道集合。
随后，将温度损耗方程代入上式，所述温度损耗方程如下：
式中，Ta为环境温度，Lb为管道b的长度，Ta为环境温度。
由上得到如下公式：
通常，数值小于0.1。因而，对上式中的指数想进行泰勒展开并忽略高阶项：
由上式可知，在影响管道热损耗的相关因素中，温度损耗系数、管道长度和环境温度均为常数，仅管道入口温度为可变变量。由此可知，管道入口温度越低，管道热损耗越小，热网的整体热损耗也越小。
因而，基于上式，应用任一管道的入口温度最小值，获取该管道的热损耗；随后，将各管道的热损耗之和作为热网热损耗最小值。
将热网热损耗最小值与预先设定的各热负荷的消耗量相加，得到所述热网的供热量最小值。随后，根据能量守恒定律，在不考虑热流方程约束的条件下，以热源生产成本为目标进行热源出力的经济分配，获取各热源的供热量。基于热网模型中的热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算所述热网中各热源的循环水质量流量。
与此同时，基于热网模型中的热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算热网中各热负荷的循环水质量流量。
其次，将上一步骤中计算获取的热网中各热负荷的循环水质量流量和各热源的循环水质量流量代入热网模型中的热网流量模型，计算热网中各供水管道和回水管道的循环水质量流量。
最后，应用上一步骤中计算获取的热网中各供水管道和回水管道的循环水质量流量替代原最优热流模型中应用热网流量模型构成的限定条件，将固定的各供水管道和回水管道的循环水质量流量和热网传热模型作为限定条件，根据下式求解热最优热流：
min CCU
s.t.热网传热模型
和为定值
式中，和分别为供水管道a和回水管道d的循环水质量流量。
在和为定值的情况下，热网传热模型约束为线性约束，同时目标函数为二次形式，则上式表示的最优热流模型为一个线性约束下的凸优化模型，利用商业优化软件可进行高效求解，得到各热源出力、管道温度分布等最优热流结果。
本示例中，应用固定的热网管道的循环水质量流量将热流方程约束退化为线性等式约束，简化了热网的最优热流模型，实现了最优热流模型的高效计算，有助于电网详细模型的推广和热网运行特性的全面反映。
图2为本发明具体实施例的一种热网热流计算设备的结构示意图，如图2所示，该设备包括：至少一个处理器201；以及与所述处理器201通信连接的至少一个存储器202，其中：所述存储器202存储有可被所述处理器201执行的程序指令，所述处理器201调用所述程序指令能够执行上述各实施例所提供的热网热流计算方法，例如包括：S1，基于热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算各热负荷的循环水质量流量；基于热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算各热源的循环水质量流量；其中，所述各热源的供热量根据供热量最小值对各热源进行分配获取；所述供热量最小值为所述各热负荷的消耗量与热网热损耗最小值之和；所述热网热损耗最小值根据所述各供水管道和回水管道的入口温度最小值计算获取；S2，基于热网流量模型，根据所述各热负荷和各热源的循环水质量流量，计算各供水管道和回水管道的循环水质量流量；S3，将所述各供水管道和回水管道的循环水质量流量和热网传热模型作为限定条件，求解热流方程对应的热流，将所述热流作为热网最优热流。
本发明实施例还提供一种非暂态计算机可读存储介质，该非暂态计算机可读存储介质存储计算机指令，该计算机指令使计算机执行对应实施例所提供的热网热流计算方法，例如包括：S1，基于热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算各热负荷的循环水质量流量；基于热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算各热源的循环水质量流量；其中，所述各热源的供热量根据供热量最小值对各热源进行分配获取；所述供热量最小值为所述各热负荷的消耗量与热网热损耗最小值之和；所述热网热损耗最小值根据所述各供水管道和回水管道的入口温度最小值计算获取；S2，基于热网流量模型，根据所述各热负荷和各热源的循环水质量流量，计算各供水管道和回水管道的循环水质量流量；S3，将所述各供水管道和回水管道的循环水质量流量和热网传热模型作为限定条件，求解热流方程对应的热流，将所述热流作为热网最优热流。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后，本申请的方法仅为较佳的实施方案，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (9)
1.一种热网热流计算方法，其特征在于，包括：
S1，基于热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算各热负荷的循环水质量流量；
基于热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算各热源的循环水质量流量；
其中，所述各热源的供热量根据供热量最小值对各热源进行分配获取；所述供热量最小值为所述各热负荷的消耗量与热网热损耗最小值之和；所述热网热损耗最小值根据所述各供水管道和回水管道的入口温度最小值计算获取；
S2，基于热网流量模型，根据所述各热负荷和各热源的循环水质量流量，计算各供水管道和回水管道的循环水质量流量；
S3，将所述各供水管道和回水管道的循环水质量流量和热网传热模型作为限定条件，求解热流方程对应的热流，将所述热流作为热网最优热流。
2.根据权利要求1所述的方法，其特征在于，所述步骤S1中，基于热负荷模型，应用所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值计算各热负荷的循环水质量流量，进一步包括：
将所述热网中预先设定的各热负荷的消耗量和各供水管道和回水管道的入口温度最小值代入下式：
式中，为热负荷k的消耗量，cp为水的比热容，为热负荷k的循环水质量流量，为热负荷k的供水温度，为热负荷k的回水温度，SHL为热负荷集合；分别将与热负荷k连接的供水管道和回水管道的入口温度最小值作为和
3.根据权利要求1所述的方法，其特征在于，所述步骤S1中，基于热源模型，应用所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值计算各热源的循环水质量流量，进一步包括：
将所述热网中各热源的供热量和各供水管道和回水管道的入口温度最小值代入下式：
式中，为热源j的供热量，cp为水的比热容，为热源j的循环水质量流量，为热源j的供水温度，为热源j的回水温度，SHS为热源集合；分别将与热源j连接的供水管道和回水管道的入口温度最小值作为和
4.根据权利要求1所述的方法，其特征在于，所述步骤S1中，热网热损耗最小值根据所述各供水管道和回水管道的入口温度最小值计算获取，进一步包括：
将任一管道的入口温度最小值代入下式，获取该管道的热损耗：
式中，ΔQb为管道b的热损耗，λb为管道b的温度损耗系数，为管道b的入口温度，Lb为管道b的长度，Ta为环境温度，Spipe为供水管道和回水管道集合；
将各管道的热损耗之和作为热网热损耗最小值。
5.根据权利要求1所述的方法，其特征在于，所述步骤S2中，所述热网流量模型包括流量连续性方程和压力降落平衡方程；
所述流量连续性方程包括：
其中，和分别为供水管道a和回水管道d的循环水质量流量，和分别为热源j和热负荷k的循环水质量流量，和分别为以节点i为末端的供水管道集合和回水管道集合，和分别为节点i为首端的供水管道集合回水管道集合，和分别为连接节点i的热源集合和热负荷集合，Snode为热网节点集合；所述节点包括热源和热负荷；
所述压力降落平衡方程包括：
式中，和分别为供水管道a和回水管道d的水压降落，和分别为供水管道a和回水管道d的压力损耗系数，和分别为供水管道a和回水管道d的循环水质量流量，为供水管道集合，为回水管道集合，为供水管道环路集合，为回管道环路集合。
6.根据权利要求1至5中任一权利要求所述的方法，其特征在于，所述步骤S3中，所述热网传热模型包括热负荷模型、热源模型、温度损耗方程和温度混合方程；
所述温度损耗方程包括：
式中，和分别为管道b的出口温度和入口温度，Ta为环境温度，λb为管道b的温度损耗系数，Lb为管道b的长度，cp为水的比热容，为管道b的循环水质量流量，Spipe为供水管道和回水管道集合；
所述温度混合方程包括：
式中，和分别为供水管道a和回水管道d的出口温度，和分别为供水管道a和回水管道d的入口温度，和分别为供水管道a和回水管道d的循环水质量流量,和分别为热源j的供水温度和回水温度，和分别为热负荷k的供水温度和回水温度，和分别为热源j和热负荷k的循环水质量流量，和分别为节点i的供水混合温度和回水混合温度，和分别为以节点i为末端的供水管道集合和回水管道集合，和分别为节点i为首端的供水管道集合回水管道集合，和分别为连接节点i的热源集合和热负荷集合，Snode为热网节点集合。
7.根据权利要求1至5中任一权利要求所述的方法，其特征在于，所述步骤S3中，当热源为燃煤产热机组时，所述热流方程如下：
min CCU
式中，
为热网成本，为热流，b0,i、b1,i和b2,i分别为成本系数的常数项、一次项和二次项，SCU为燃煤产热机组集合。
8.一种热网热流计算设备，其特征在于，包括：
至少一个处理器；以及与所述处理器通信连接的至少一个存储器，其中：所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行如权利要求1至7任一所述的方法。
9.一种非暂态计算机可读存储介质，其特征在于，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行如权利要求1至7任一所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810079288.4A CN108280302A (zh) | 2018-01-26 | 2018-01-26 | 一种热网热流计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810079288.4A CN108280302A (zh) | 2018-01-26 | 2018-01-26 | 一种热网热流计算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108280302A true CN108280302A (zh) | 2018-07-13 |
Family
ID=62805477
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810079288.4A Pending CN108280302A (zh) | 2018-01-26 | 2018-01-26 | 一种热网热流计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108280302A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110762602A (zh) * | 2019-10-08 | 2020-02-07 | 华南理工大学 | 一种分布式能源供热管网储能的量化计算方法 |
| CN113108353A (zh) * | 2021-04-26 | 2021-07-13 | 杭州鸿晟电力设计咨询有限公司 | 考虑热损非线性的变流量变温度区域供热系统调度方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11305646A (ja) * | 1998-04-23 | 1999-11-05 | Mitsubishi Electric Corp | 初期値生成方法および装置 |
| CN106294961A (zh) * | 2016-08-03 | 2017-01-04 | 清华大学 | 一种计及管网储热效益的集中供热系统建模方法 |
| CN106339772A (zh) * | 2016-08-11 | 2017-01-18 | 清华大学 | 基于供热管网储热效益的热‑电联合优化调度方法 |
| CN106845701A (zh) * | 2017-01-11 | 2017-06-13 | 东南大学 | 一种基于热网和房屋热惯性的综合能源系统优化方法 |
| CN109579100A (zh) * | 2018-10-26 | 2019-04-05 | 南京艾科美热能科技有限公司 | 一种数据中心余热利用装置 |
-
2018
- 2018-01-26 CN CN201810079288.4A patent/CN108280302A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11305646A (ja) * | 1998-04-23 | 1999-11-05 | Mitsubishi Electric Corp | 初期値生成方法および装置 |
| CN106294961A (zh) * | 2016-08-03 | 2017-01-04 | 清华大学 | 一种计及管网储热效益的集中供热系统建模方法 |
| CN106339772A (zh) * | 2016-08-11 | 2017-01-18 | 清华大学 | 基于供热管网储热效益的热‑电联合优化调度方法 |
| CN106845701A (zh) * | 2017-01-11 | 2017-06-13 | 东南大学 | 一种基于热网和房屋热惯性的综合能源系统优化方法 |
| CN109579100A (zh) * | 2018-10-26 | 2019-04-05 | 南京艾科美热能科技有限公司 | 一种数据中心余热利用装置 |
Non-Patent Citations (1)
| Title |
|---|
| 李平等: "《利用建筑物与热网热动态特性提高热电联产机组调峰能力》", 《电力系统自动化》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110762602A (zh) * | 2019-10-08 | 2020-02-07 | 华南理工大学 | 一种分布式能源供热管网储能的量化计算方法 |
| CN110762602B (zh) * | 2019-10-08 | 2021-02-12 | 华南理工大学 | 一种分布式能源供热管网储能的量化计算方法 |
| CN113108353A (zh) * | 2021-04-26 | 2021-07-13 | 杭州鸿晟电力设计咨询有限公司 | 考虑热损非线性的变流量变温度区域供热系统调度方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110263966B (zh) | 考虑动态传热过程的电-热综合能源系统优化调度方法 | |
| Xu et al. | Quantification of flexibility of a district heating system for the power grid | |
| Awad et al. | Integrated optimal power flow for electric power and heat in a microgrid | |
| CN110175311B (zh) | 一种基于多能耦合模型的优化潮流计算方法 | |
| CN110991092A (zh) | 一种余能一体化利用计算方法及系统 | |
| CN108229753A (zh) | 一种热电耦合系统的分布式调度方法 | |
| CN112330493B (zh) | 一种能源系统的建模及综合分析的方法、装置及存储介质 | |
| CN108280302A (zh) | 一种热网热流计算方法 | |
| CN110020506B (zh) | 基于电热型综合能源系统运行优化的差分格式选择方法 | |
| CN115859686A (zh) | 考虑扩展碳排放流的综合能源系统低碳调度方法及系统 | |
| Man et al. | State estimation for integrated energy system containing electricity, heat and gas | |
| CN112686447B (zh) | 海上油气田开发多能流耦合负荷预测方法 | |
| CN113886761A (zh) | 一种综合能源系统能效分析评估方法 | |
| CN112182905B (zh) | 一种用于综合能源系统的供热管网仿真方法和装置 | |
| CN211044253U (zh) | 一种余能一体化利用计算系统 | |
| Tang et al. | Optimal dispatch of integrated energy campus microgrids considering the time-delay of pipelines | |
| Loukarakis et al. | A sequential programming method for multi-energy districts optimal power flow | |
| CN111047069B (zh) | 一种考虑热网储能特性的热电联合调度灵活性评估方法 | |
| CN116739404A (zh) | 一种综合能源系统熵态模型中网络化特征模块的构建方法 | |
| CN117391249A (zh) | 一种电气互联综合能源系统凸优化方法 | |
| Dai et al. | Modeling and analysis of electrical heating system based on entransy dissipation-based thermal resistance theory | |
| CN111062598A (zh) | 综合能源系统的分布式优化调度方法及系统 | |
| CN111523792B (zh) | 综合能源系统调度参数计算、设备控制方法及相关装置 | |
| CN114825325A (zh) | 一种配电网运行参数优化方法及装置 | |
| Chen et al. | Multi-energy flow calculation considering the convexification network constraints for the integrated energy system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180713 |
|
| RJ01 | Rejection of invention patent application after publication |