[发明专利]一种基于双侧等效模型的热网稳态运行状态估计方法

摘要

本发明涉及一种基于双侧等效模型的热网稳态运行状态估计方法，属于综合能源系统的运行和控制技术领域。该方法包括：建立热网双侧等效模型、基于热网双侧等效模型对稳态运行的热网进行状态估计、建立状态估计的目标函数、建立热网稳态运行的约束条件、利用拉格朗日乘数法和牛顿‑拉夫逊法进行状态估计计算并对状态估计结果进行收敛性判断。本发明同时考虑了供回水网络，并且考虑了传热过程，可以很好地解决热网稳态状态估计问题，在量测配置不全的情下补全量测，存在坏数据时辨识坏数据。

主权项

一种基于热网双侧等效模型的热网稳态状态估计方法，其特征在于该方法包括以下步骤：(1)建立热网双侧等效模型的节点‑支路关联矩阵，包括：在热网双侧等效模型中同时考虑热网供水支路和回水支路，并将热源和热负荷等效为连接支路，对于由N个节点、B条支路构成的热网，形成以下表示节点和支路关系的矩阵：(1‑1)一个节点‑支路关联矩阵A，节点‑支路关联矩阵A表示网络中节点和支路的拓扑关系，矩阵A由0、1、‑1三个元素组成，A中元素定义如下：其中，i为热网中的任意一个节点，j为热网中任意一条支路；(1‑2)一个正节点‑支路关联矩阵Af：正节点‑支路关联矩阵Af表示各支路的首端节点与支路的关系，Af＝{A│Aij＞0}，Af中元素定义如下：(1‑3)一个负节点‑支路关联矩阵At：负节点‑支路关联矩阵At表示各支路的末端节点与支路的关系，At＝{‑A|Aij＜0}，At中元素定义如下式：(2)基于热网双侧等效模型对稳态运行的热网进行状态估计：(2‑1)设定热网状态估计的收敛精度δ和最大循环次数d，初始化时设循环次数a为0；(2‑2)从热网的数据采集与监视控制系统中获取实时测量的t时刻热网的运行数据，包括热网中各节点压力H，任意两个节点间支路的流量m，两个节点间支路的首端温度Tf和末端温度Tt，热源和热负荷等效成的连接支路，热功率φq，其中上标q表示连接支路，上述运行数据构成一个测量值列向量zh；(2‑3)将热网所有待估计的状态量构成一个列向量xh，其中包括热网中各节点压力以及任两个节点间支路的首端温度和末端温度(2‑4)建立一个描述热网状态量与测量值之间关系的量测函数f(x)，f(x)＝f(xh)，f(xh)为热力系统潮流方程组，热力系统潮流方程组包括以下方程：(2‑4‑1)一个支路压力损失方程：支路压力损失方程表示一条支路两端节点的压力差，支路压力损失方程的矩阵形式如下：ATH＝ΔH‑Hp其中H为上述步骤(2‑2)中的热网各节点压力组成的列向量，AT为上述步骤(1‑1)中节点‑支路关联矩阵A的转置，Hp为支路上泵的扬程组成的列向量，a、b、c为泵参数，从泵的产品铭牌上获取，mp为泵所在支路的流量，ΔH为热网中每条支路压力损失组成的列向量，支路压力损失ΔH通过下式计算得到：ΔH＝K·m·|m|其中，K为热网中支路的摩阻系数，取值为10～500帕/(千克/秒)2，m为热网中任一支路流量；(2‑4‑2)一个连接支路热功率方程：热功率方程表示连接支路q的首末端温度关系，表示为如下形式：φq=Cpmq(Tfq-Ttq)]]>其中，上标q表示连接支路，φq为连接支路的用热功率，热负荷处的用热功率为正，热源处的用热功率为负，Cp为供热介质的比热容，由流体的物性参数表获取，mq为连接支路流量，为连接支路首端温度，为连接支路末端温度；(3)根据上述步骤(2‑2)的测量值，建立一个热网稳态运行状态估计的目标函数如下：minJ(xh)＝min{[zh‑f(xh)]TW[zh‑f(xh)]}其中W为测量值的协方差矩阵，上标T表示矩阵转置，J(xh)表示目标函数表达式；(4)建立热网稳态运行的约束条件c(xh)，包括：(4‑1)对所有节点建立流量连续性约束，流量连续性约束表示为如下矩阵形式：AM＝0其中，M为热网每条支路流量构成的列向量，将供水支路和回水支路统一等效为普通支路，用上标p表示，连接支路等效为特殊支路，用上标q表示，则M表示为：M=MpMq]]>其中，Mp表示供水支路和回水支路即普通支路流量组成的子向量，Mq表示连接支路流量组成的子向量；(4‑2)对热网中的所有节点建立温度混合约束：(∑mout)Tn＝∑(minTin)其中，mout为供热介质流出节点的支路流量，min为供热介质流入节点的支路流量，Tn为节点处供热介质混合后的温度，Tin为不同支路供热介质在节点处混合前的温度；用不同支路末端温度Tt代替支路供热介质在节点处混合前的温度Tin，则节点温度混合约束表示为如下矩阵形式：diag(AfM)Tn＝Atdiag(M)Tt其中，Af、At分别为上述步骤(1‑2)中的正节点‑支路关联矩阵和上述步骤(1‑3)中的负节点‑支路关联矩阵，diag(·)表示对角阵；(4‑3)对热网中的所有普通支路建立支路温降约束，支路温降约束的矩阵形式如下：Ttp-Ta=(Tfp-Ta)e-λLCpMp]]>其中，Ta为环境温度，为普通支路末端温度，为普通支路首端温度，L为普通支路长度，λ为热网中普通支路的散热系数，从相应数据手册中获取，e是自然对数，Cp为供热介质的比热容，Mp表示供水支路和回水支路即普通支路流量组成的子向量，；(5)利用拉格朗日乘数法，将上述步骤(3)的目标函数和上述步骤(4)的约束条件构成一个拉格朗日函数如下：L(xh,ω)＝J(xh)+ωTc(xh)其中，J(xh)为上述步骤(3)的目标函数，ω为拉格朗日乘子，c(xh)为上述步骤(4)建立的热网稳态运行的约束条件，上标T是矩阵转置；利用最优化理论中的牛顿‑拉夫逊法，求解上述热网稳态运行时的拉格朗日函数，得到热网稳态运行时的状态估计结果；(6)对上述步骤5的状态估计结果进行收敛性判断：若循环次数a达到预设循环次数d，即a≥d，则将本次状态估计结果作为t时刻基于双侧等效模型的热网稳态状态估计结果；若循环次数a未达到预设循环次数d，即a＜d，则进一步根据热网状态估计的精度δ对状态估计结果收敛性进行判断：若最近相邻两次状态估计结果中的状态变xa和xa‑1的差值小于状态估计精度δ，即max|xa‑xa‑1|＜δ，则将本次状态估计结果作为t时刻基于双侧等效模型的热网稳态状态估计结果，若最近两次状态估计结果中的状态变量估计值xa和xa‑1的差值大于或等于状态估计精度δ，即max|xa‑xa‑1|≥δ，则更新状态变量，并根据本次状态估计所得的温度值更新热网中节点压力和支路首末端温度，同时使a＝a+1，并返回步骤4，继续本次状态估计过程。