供热管网水力仿真模型辨识修正方法及系统、操作方法技术方案

本发明专利技术涉及一种供热管网水力仿真模型辨识修正方法及系统，本供热管网水力仿真模型辨识修正方法包括：步骤S1，选取供热系统在稳态工况下的多组实测数据；步骤S2，计算各管段的阻力系数理论值；以及步骤S3，建立供热管网水力仿真模型，利用粒子群算法结合多组实测数据辨识计算各阻力系数修正量，并获得修正过的阻力系数。发明专利技术的供热管网水力仿真模型辨识修正方法及系统，采用粒子群算法在每一代产生不同的阻力系数修正量向量，并通过供热管网水力仿真模型中不同阻力系数修正量向量在多个工况下管网的流动状态与相应实测数据进行对比，辨识出最优的阻力系数修正量。修正后的模型能够更好地描述实际物理系统，提高供热管网水力计算仿真的计算精度。采用并行计算架构进行粒子群算法阻力系数修正量的辨识，辨识效率高。

Hydraulic network simulation model identification correction method and system and operation method of heating pipe network

The invention relates to a heating pipe network hydraulic simulation model identification correction method and system, including the identification of the hydraulic simulation model of heating pipe network correction method: step S1, multi group data under steady condition selection of heating system; step S2 theory to calculate the resistance coefficient of each pipe; and step S3, a hydraulic simulation of heating the network model, the correction calculation of the resistance coefficient of the particle swarm algorithm combined with the measured data of multi group identification, and the modified drag coefficient obtained. The invention of the heating pipe network hydraulic simulation model identification correction method and system using particle swarm algorithm in each generation resistance coefficient of different correction vectors, and through different resistance coefficient of hydraulic simulation model in the heating network correction contrast flow state vector network in a number of conditions with the corresponding measured data, identify the resistance coefficient optimal correction. The revised model can better describe the actual physical system, and improve the calculation accuracy of the hydraulic network simulation of the heating pipe network. The parallel computing architecture is used to identify the resistance coefficient correction of particle swarm optimization algorithm, and the identification efficiency is high.

【技术实现步骤摘要】
供热管网水力仿真模型辨识修正方法及系统、操作方法
本专利技术属于供热领域，具体涉及一种供热管网水力仿真模型辨识修正方法及系统、操作方法。
技术介绍
建立符合实际供热管网运行特性的水力仿真模型是实现热网智慧化管理的重要基础。随着供热规模的日益扩大和网络结构的日趋复杂，建立严格精确的供热管网的水力仿真模型变得愈加困难。通常管网的连接结构、节点标高等参数较为容易得到，而管道阻力特性受多种物理参数影响，且随时间变化，通过严格机理计算获得的理论阻力特性往往与实际特性存在偏差，导致建立的管网水力仿真模型不能很好地模拟实际运行工况。为此，需要充分利用实际运行历史数据，采用辨识的方法修正管段阻力特性。目前大多数供热管网在热力站处有供、回水压力测点和流量测点，而其它管段和节点处并不存在测点，因而难以通过直接计算的方法获得整个管网管段的阻力系数。
技术实现思路
本专利技术的目的是提供一种供热管网水力仿真模型辨识修正方法及系统，以解决供热管网水力仿真模型不能准确描述实际管网运行的技术问题。为了解决上述技术问题，本专利技术提供了一种供热管网水力仿真模型辨识修正方法，包括如下步骤：步骤S1，选取供热系统在稳态工况下的多组实测数据；步骤S2，计算各管段的阻力系数理论值；以及步骤S3，建立供热管网水力仿真模型，利用粒子群算法结合多组实测数据辨识计算各阻力系数修正量，并获得修正过的阻力系数。进一步，所述步骤S1中选取供热系统在稳态工况下的多组实测数据，即识别供热管网系统的多个稳态运行工况，选取多组稳态运行工况的实测数据，其方法包括：步骤S11，工况数据预处理；步骤S12，求各热力站的稳态工况时间区间；步骤S13，求整个供热管网系统的稳态工况时间区间；以及步骤S14，选取所述多组实测数据。进一步，所述步骤S2中计算各管段的阻力系数理论值的方法包括：通过如下机理公式计算出各管段的阻力系数理论值；即δi为管段i的阻力系数(Pa/(t/h)2)；K为管壁的当量绝对粗糙度；di为管道内径，单位为m；li为管段长度，单位为m；ldi为管段局部阻力当量长度，单位为m；ρ为管内热水的平均密度，单位为kg/m3；D为管段数量。进一步，所述步骤S3中建立供热管网水力仿真模型的方法包括：将供热管网转换为由节点和区段构成的有向图模型，即管网图；其中节点表示存在流量进出的点，用集合V表示，V＝{V1,V2,…,Vn}，式中，n为管网中的节点个数；区段表示节点间的连接管段，用集合E表示，E＝{E1,E2,…,Em}，式中D为区段数即管段数量；转换为有向图，则表示为G＝<V,E>；根据网络图论得到供热管网的关联矩阵A和基本回路矩阵B，其中A为n×D阶矩阵，B为s×D阶矩阵，s为基本回路个数s＝D-n+1；构建供热管网水力仿真模型，即AGT＝0；BΔHT＝0；式中，G为记录管网图中各管段内体积流量的行向量，即G＝[G1,G2,…,GD]；ΔH为记录管网图中各管段总阻力损失的列向量，即ΔH＝ζ*|G|*G+Z-Hb；式中，ζ为管网中各管段的阻力系数矩阵；ζ＝diag{δ1+Δδ1,δ2+Δδ2,…,δD+ΔδD}；其中δ1、δ2、……、δD为阻力系数理论值，Δδ1、Δδ2、……、ΔδD为阻力系数修正量，即需要辨识的参数；D表示待辨识的阻力系数修正量的维数；Z为管段两节点位能差的列向量Z＝[Z1,Z2,…,ZD]T；Hb为热网中的水泵扬程列向量Hb＝[Hb1,Hb2,…,HbD]T；在供热管网水力仿真模型应用于辨识时，将其抽象表示为如下形式：Y＝Func(X,δ)式中X表示输入的多组实测数据；Y表示输出向量；δ表示待辨识管段的阻力系数向量,δ＝[δ1+Δδ1,δ2+Δδ2,…,δD+ΔδD]。进一步，所述步骤S3中利用粒子群算法结合多组实测数据辨识计算各阻力系数修正量，即采用粒子群算法在每一代产生不同的阻力系数修正量，并通过供热管网水力仿真模型中不同阻力系数修正量向量在多个工况下管网的流动状态与相应实测数据进行对比，辨识出最优的阻力系数修正量。进一步，采用粒子群算法结合多组实测数据辨识计算各阻力系数修正量向量的方法包括如下步骤：步骤S31，粒子群算法参数设置，即微粒个数N、微粒维数D，以及决策变量上界U＝[Δδ1u,Δδ2u,…,ΔδDu]和下界L＝[Δδ1l,Δδ2l,…,ΔδDl]；步骤S32，初始化，即读取阻力系数理论值，生成N个微粒，即N组阻力系数修正量向量；步骤S33，建立用于阻力系数修正量辨识的目标函数，以计算各微粒的目标函数值；所述目标函数如下：当阻力系数修正量Δδ1,Δδ2,…,ΔδD作为决策变量给定值后，通过供热管网水力仿真模型计算得到相应输出数据；对每一个微粒k∈{1,2,…,N}，将其与阻力系数理论值向量[δ1,δ2,…,δD]相加，得到该微粒代表的阻力系数向量将δk输入供热管网水力仿真模型，并依次利用第1、2，……，第NC组实测数据对其进行水力计算，得到相应工况下各热力站的压力计算值；式中，PCij：稳态工况j下由仿真模型计算得到的热力站i的压力计算值；POij：稳态工况j下热力站i对应的压力测量值；NP：压力测点数量；步骤S34，微粒飞行，产生下一代微粒；步骤S35，达到设定的运行代数，算法终止，否则回到步骤S33继续执行。又一方面，本专利技术还提供了一种供热管网水力仿真模型辨识修正系统，包括：粒子群算法单元，产生每一代的阻力系数修正量对应的微粒，并将微粒写入数据库，并且读取各微粒的目标函数值；目标函数计算单元，建立用于阻力系数修正量辨识的目标函数，驱动多个供热管网水力仿真模型，即产生N个仿真模型进行并行的模拟计算，并从数据库获取微粒，通过仿真模型计算微粒的目标函数值并写回数据库。数据库单元，存储阻力系数修正量向量的群体和目标函数值，以及目标函数的计算状态，并且还将稳态工况的实测数据以及每一代优化计算的结果也存入数据库中。进一步，所述目标函数如下：当阻力系数修正量Δδ1,Δδ2,…,ΔδD作为决策变量给定值后，通过供热管网水力仿真模型计算得到相应输出数据；对每一个微粒k∈{1,2,…,N}，将其与阻力系数理论值向量[δ1,δ2,…,δD]相加，得到该微粒代表的阻力系数向量将δk输入供热管网水力仿真模型，并依次利用第1、2，……，第NC组实测数据对其进行水力计算，得到相应工况下各热力站的压力计算值；式中，PCij：稳态工况j下由仿真模型计算得到的热力站i的压力计算值；POij：稳态工况j下热力站i对应的压力测量值；NP：压力测点数量。第三方面，本专利技术还提供了一种供热管网水力仿真模型辨识修正系统的操作方法，包括如下步骤：步骤Sa，建立待辨识管网的拓扑结构模型，以生成管网各部件的唯一标识ID；步骤Sb，获得工况实测数据和计算阻力系数理论值的准备工作；步骤Sc，设置粒子群算法单元中的算法参数，启动该单元；步骤Sd，启动目标函数计算单元。进一步，所述步骤Sb中所述准备工作包括：将待辨识管段的ID和阻力系数理论值对应地整理到配置文件中，供粒子群算法单元读取；以及将筛选的多组实测数据写入数据库中，供目标函数计算单元查询，每一组工况数据包括工况编号、热源或热力站ID、热源或热力站对应的流量和供回水压力。本专利技术的有益效果是，本专利技术的供热管网水力仿真模型辨本文档来自技高网...

【技术保护点】
一种供热管网水力仿真模型辨识修正方法，其特征在于，包括如下步骤：步骤S1，选取供热系统在稳态工况下的多组实测数据；步骤S2，计算各管段的阻力系数理论值；以及步骤S3，建立供热管网水力仿真模型，利用粒子群算法结合多组实测数据辨识计算各阻力系数修正量，并获得修正过的阻力系数。

【技术特征摘要】
1.一种供热管网水力仿真模型辨识修正方法，其特征在于，包括如下步骤：步骤S1，选取供热系统在稳态工况下的多组实测数据；步骤S2，计算各管段的阻力系数理论值；以及步骤S3，建立供热管网水力仿真模型，利用粒子群算法结合多组实测数据辨识计算各阻力系数修正量，并获得修正过的阻力系数。2.根据权利要求1所述的供热管网水力仿真模型辨识修正方法，其特征在于，所述步骤S1中选取供热系统在稳态工况下的多组实测数据，即识别供热管网系统的多个稳态运行工况，选取多组稳态运行工况的实测数据，其方法包括：步骤S11，工况数据预处理；步骤S12，求各热力站的稳态工况时间区间；步骤S13，求整个供热管网系统的稳态工况时间区间；以及步骤S14，选取所述多组实测数据。3.根据权利要求2所述的供热管网水力仿真模型辨识修正方法，其特征在于，所述步骤S2中计算各管段的阻力系数理论值的方法包括：通过如下机理公式计算出各管段的阻力系数理论值；即δi为管段i的阻力系数(Pa/(t/h)2)；K为管壁的当量绝对粗糙度；di为管道内径，单位为m；li为管段长度，单位为m；ldi为管段局部阻力当量长度，单位为m；ρ为管内热水的平均密度，单位为kg/m3；D为管段数量。4.根据权利要求3所述的供热管网水力仿真模型辨识修正方法，其特征在于，所述步骤S3中建立供热管网水力仿真模型的方法包括：将供热管网转换为由节点和区段构成的有向图模型，即管网图；其中节点表示存在流量进出的点，用集合V表示，V＝{V1,V2,…,Vn}，式中，n为管网中的节点个数；区段表示节点间的连接管段，用集合E表示，E＝{E1,E2,…,Em}，式中D为区段数即管段数量；转换为有向图，则表示为G＝<V,E>；根据网络图论得到供热管网的关联矩阵A和基本回路矩阵B，其中A为n×D阶矩阵，B为s×D阶矩阵，s为基本回路个数s＝D-n+1；构建供热管网水力仿真模型，即AGT＝0；BΔHT＝0；式中，G为记录管网图中各管段内体积流量的行向量，即G＝[G1,G2,…,GD]；ΔH为记录管网图中各管段总阻力损失的列向量，即ΔH＝ζ*|G|*G+Z-Hb；式中，ζ为管网中各管段的阻力系数矩阵；ζ＝diag{δ1+Δδ1,δ2+Δδ2,…,δD+ΔδD}；其中δ1、δ2、……、δD为阻力系数理论值，Δδ1、Δδ2、……、ΔδD为阻力系数修正量，即需要辨识的参数；D表示待辨识的阻力系数修正量的维数；Z为管段两节点位能差的列向量Z＝[Z1,Z2,…,ZD]T；Hb为热网中的水泵扬程列向量Hb＝[Hb1,Hb2,…,HbD]T；在供热管网水力仿真模型应用于辨识时，将其抽象表示为如下形式：Y＝Func(X,δ)式中X表示输入的多组实测数据；Y表示输出向量；δ表示待辨识管段的阻力系数向量,δ＝[δ1+Δδ1,δ2+Δδ2,…,δD+ΔδD]。5.根据权利要求4所述的供热管网水力仿真模型辨识修正方法，其特征在于，所述步骤S3中利用粒子群算法结合多组实测数据辨识计算各阻力系数修正量，即采用粒子群算法在每一代产生不同的阻力系数修正量，并通过供热管网水力仿真模型中不同阻力系数修正量向量在多个工况下管网的流动状态与相应实测数据进行对比，辨识出最优的阻力系数修正量。6.根据权利要求5所述的供热管网水力仿真模型辨识修正方法，其特征在于，采用粒子群算法结合多组实测数据辨识计算各阻力系数修正量向量的方法包括如下步骤：步骤S31，粒子群算法参数设置...

【专利技术属性】
技术研发人员：于春娣，方大俊，瞿广峰，
申请(专利权)人：常州英集动力科技有限公司，
类型：发明
国别省市：江苏,32