首页 / 专利库 / 燃气和蒸汽发生器 / 热电联产 / 一种应用于热电耦合网络的异步分布式状态估计方法

一种应用于热电耦合网络的异步分布式状态估计方法

阅读:1014发布:2020-05-20

专利汇可以提供一种应用于热电耦合网络的异步分布式状态估计方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种应用于热电耦合网络的异步分布式状态估计方法,包括步骤:1)获取热电耦合网络的量测数据;2)用 节点 法建立热网拟动态模型,使用热网量测数据计算热网拟动态模型系数;3)基于量测数据与拟动态模型系数构造热电耦合网络的状态估计目标函数;4)构造目标函数的增广拉格朗日函数,通过交替方向乘子法将增广拉格朗日函数进行分布式求解,并使用异步计算 框架 对各子问题的求解时长进行控制,最终求解得到热电联合网络的运行状态数据;5)输出热电耦合网络的运行状态。本发明将热电耦合网络的状态估计问题分解并将各子问题进行独立求解,保证解的可行性以及子系统的相互独立;用异步计算框架调节子问题计算时间,改善计算效率。,下面是一种应用于热电耦合网络的异步分布式状态估计方法专利的具体信息内容。

1.一种应用于热电耦合网络的异步分布式状态估计方法,其特征在于,所述的异步分布式状态估计方法包括以下步骤:
S1、获取热电耦合网络的量测数据;
S2、用节点法建立热网拟动态模型,使用热网量测数据计算热网拟动态模型中的流速相关系数以及温度时滞系数;
S3、基于量测数据与热网拟动态模型系数构造热电耦合网络的状态估计目标函数;
S4、基于热电耦合网络模型,构造热电耦合网络的状态估计目标函数的增广拉格朗日函数,通过交替方向乘子法将增广拉格朗日函数分解为对应热电耦合网络各子系统的子问题进行求解,并使用异步计算框架对各子问题的求解时长进行控制,最终求解得到热电耦合网络的运行状态数据;
S5、输出热电耦合网络的运行状态。
2.根据权利要求1所述的一种应用于热电耦合网络的异步分布式状态估计方法,其特征在于,所述的热电耦合网络的量测数据包括当前时刻电网的节点有功注入功率、节点无功注入功率、线路传输有功潮流、线路传输无功潮流及节点电压幅值,以及热网的节点压强、管道流速、节点温度、热源/热负荷的供应/消耗热功率、耗电量;所述的热电耦合网络的量测数据还包括历史时刻的节点温度和管道流速数据。
3.根据权利要求1所述的一种应用于热电耦合网络的异步分布式状态估计方法,其特征在于,所述的步骤S2中,用节点法建立热网拟动态模型,使用热网量测数据计算热网拟动态模型中的流速相关系数与温度时滞系数,从而建立当前时刻温度与历史时刻的温度、流量之间的关系,并在管道热量损耗方面对温度进行校正,包括以下步骤:
S2.1)计算拟动态模型中的流速相关系数:
式中,Rb,t、Sb,t、γb,t和φb,t为流速相关系数,b为管道编号,t和k分别为当前时刻与历史时刻的编号,Δt为每个时刻的时间长度,mP为管道内热水流速,上标P表示管道,ρ为水的密度,A为管道横截面积,L为管道长度;
S2.2)根据流速相关系数计算当前时刻管道b中不同段的温度时滞系数 上标TD表示时滞:
S2.3)基于热网拟动态模型的系数,将当前时刻的未考虑温度损耗的管道流出端温度与多个时刻的管道流入端温度的关系用下式表示:
式中, 为当前时刻的管道b未考虑温度损耗的管道流出端温度,上标out表示流出端,上标*表示不考虑热损耗, 表示时刻管道的流入端温度,上标in表示流入端;
S2.4)基于未考虑温度损耗的管道流出端温度以及热网拟动态模型系数,计算考虑温度损耗的管道流出端温度
式中,TA为外界温度,ch为管道热传导系数,Cp为水的比热容。
4.根据权利要求3所述的一种应用于热电耦合网络的异步分布式状态估计方法,其特征在于,所述的步骤S3中,基于量测数据与热网拟动态模型系数,构造热电耦合网络的状态估计目标函数:
式中,i和j分别为热网、电网子系统的编号,SDHN和SEPS分别为热网、电网子系统的编号集,项 为第i个热网子系统的量测残差加权和,项 为第j个电
网子系统的量测残差加权和,以上两项分别表示为:
式中,ps、Ppump、m、Φ、T分别代表热网中的节点压强、水泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度变量, 分别代表节点压强、水
泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度的量测值,分别代表由量测方程得到的节点压强、水泵耗电量、管道流速、
热源/热负荷的供应/消耗热功率以及节点温度计算值,Wips、Wipump、Wim、WiΦ、WiT分别代表节点压强、水泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度对应的权重值, 分别代表节点压强、水泵耗电量、管道流速、热源/热
负荷的供应/消耗热功率以及节点温度量测项的编号集,i为量测项的编号;
式中,V、θ、Pinj、Qinj、Pflow、Qflow分别代表节点电压幅值、节点电压相、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流变量,
分别代表节点电压幅值、节点有功注入功率、节点无功注入
功率、线路有功传输潮流、线路无功传输潮流的量测值, 分别
代表由量测方程得到的节点电压幅值、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流计算值,WiV、 分别代表节点电压幅
值、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流对应的权重值, 分别代表节点电压幅值、节点有功注入功率、
节点无功注入功率、线路有功传输潮流、线路无功传输潮流量测项的编号集;
在式(9)中,除了满足管道拟动态特性约束式(1)-式(6)和温度损耗式(7)外,热网各变量需满足节点水流流入流出量平衡约束、温度混合约束、压强-流量约束、节点供热/耗热量约束、以及水泵耗能约束,公式分别如下:
式中, 和 分别为节点处对应流入、流出管道的编号集合, 和 分别表示
管道流入端、流出端的压强,Kf为压强损耗系数, 为节点温度, 和 分别为热电联产机组供热量和节点热负荷, 和 分别为流经热电联产机组、负荷的热水流速,iS和iLD分别为热电联产机组、负荷节点的编号,Cp为比热容系数, 和
分别表示热电联产机组与负荷处流入端、流出端水流温度, 为流经水泵的水流速率, 和 分别为水泵流入端、流出端压强, 为水泵的耗电量,iP为水泵的
编号,ρ为水的密度, 为水泵效率系数,cm为热电联产机组热-电转换比例, 为热电联产机组发电量;
在式(10)中,电网各变量需满足节点功率平衡约束以及线路传输功率约束:
Pij,t=Vi2gij-ViVj(gijcosθij+bijsinθij)    (20)
Qij,t=-Vi2(bsi+bij)-ViVj(gijsinθij+bijcosθij)     (21)
式中,i和j均为电网中节点的编号,t为所在时刻编号,Pi,t和Qi,t分别为节点有功、无功注入功率,Pij,t和Qij,t分别为线路有功、无功传输功率,Vi和θij为节点电压相量和两节点之间电压相角差,Gij和Bij分别表示导纳矩阵中电导、电纳元素,gij和bij分别表示支路电导和电纳,bsi表示对地电纳;
热网、电网的耦合约束由水泵耗电公式(16)、热电联产机组产电量-供热量关系式(17)组成。
5.根据权利要求4所述的一种应用于热电耦合网络的异步分布式状态估计方法,其特征在于,所述的步骤S4包括:
S4.1)基于热电耦合网络模型中的热电耦合关系,构造目标函数的增广拉格朗日函数:
式中, 为热网子系统变量, 为电网子系统变量,λ为拉格朗日乘子,β为增广拉格朗日函数的惩罚因子,z为辅助变量,项 代表热网子系统内部约束以及与
电网子系统之间的耦合约束,即式(1)-(8)和式(11)-(17),项 代表电网子
系统内部约束以及与热网子系统之间的耦合约束,即式(17)-(21);
S4.2)对应热电耦合网络中的每一个子系统,将增广拉格朗日函数划分为各个子问题,热网、电网子问题对应的目标函数分别如下:
S4.3)对热网、电网子问题分别进行求解,使用异步计算框架对热网子问题的求解时长进行控制。
6.根据权利要求5所述的一种应用于热电耦合网络的异步分布式状态估计方法,其特征在于,所述的步骤S4.3)对热网、电网子问题分别进行求解,使用异步计算框架对热网子问题的求解时长进行控制具体包括:
S4.3.1)电网子问题直接求解,
热网子问题在求解过程中,将热网变量分为两类进行计算,第一类为拟动态模型涉及的系数Rb,t、Sb,t、γb,t、φb,t和 称为复杂变量,第二类为热网其余变量,称为普通变量,在已知时滞系数的基础上,首先对热网子问题进行求解,并用所得的普通变量解更新拟动态模型系数,即复杂变量;之后,使用更新后的拟动态模型系数构造新的热网子问题目标函数并再一次进行求解,更新拟动态模型系数;复杂变量、普通变量的交替计算将一直进行,直至前后两次解的差值低于某一特定值,视为热网子问题求解完成;
S4.3.2)若电网、热网子问题同时求解完成,则视为该轮迭代完成,通过子问题边界变量计算辅助变量与拉格朗日乘子,并进入下一轮迭代;
若热网子问题的求解早于电网子问题完成,则热网部分暂时不进行计算,待电网子问题完成求解后,更新辅助变量与拉格朗日乘子,并进入下一轮迭代;
若电网子问题的求解早于热网子问题完成,则热网在完成当前时刻的热网普通变量计算后停止计算,直接将边界节点数据与电网交换,更新辅助变量与拉格朗日乘子,进入下一轮迭代;
S4.3.3)当系统变量满足以下要求时,迭代满足收敛条件:
res
‖rk‖<∈     (25)
‖sk‖<∈dual    (26)
式中,k为迭代次数,‖rk‖和‖sk‖分别代表本次迭代计算的原始残差、对偶残差,∈res和∈dual为分别对应原始残差、对偶残差的收敛标准,原始残差、对偶残差的计算方法如下:
7.根据权利要求1所述的一种应用于热电耦合网络的异步分布式状态估计方法,其特征在于,所述的步骤S5中,输出热电耦合网络的运行状态数据包括:电网节点电压的相角与幅值、节点有功注入功率、节点无功注入功率、线路传输有功潮流、线路传输无功潮流;热网节点压强、管道流速、节点温度、热源/热负荷的供应/消耗热功率、水泵耗电量。

说明书全文

一种应用于热电耦合网络的异步分布式状态估计方法

技术领域

[0001] 本发明涉及多能流系统的状态估计技术领域,具体涉及一种应用于热电耦合网络的异步分布式状态估计方法。

背景技术

[0002] 由于能源局势日趋紧张,近年来,高效率、低排放的能源技术愈发受到能源行业的重视。其中,热电耦合网络具有较大的优势,研究表明:热电联产机组的能源利用效率很高;热电耦合技术能够有效地减少污染物与温室气体的排放;热电耦合网络中含有多类能源转换装置,网络内部电、热的耦合性强,提高了系统中能源供应的可靠性及能源使用的灵活性,故热电耦合网络得到的越来越多的使用与推广。对热电耦合网络进行合理的分析和控制,需要首先获取准确的系统运行状态。通常运行状态信息来自于系统各处的量测设备的数据,为了减小误差、排除故障设备数据,状态估计常被用于处理量测数据,从而获取最接近系统真实运行状态。因此,如何对热电耦合网络进行快速、准确的状态估计是一个重要的研究课题。
[0003] 热电耦合网络的特性之一即为热网与电网的耦合关系,故对于热电耦合网络进行整体的状态估计是有必要的,而不是单独对电网、供热网络进行状态估计。联合状态估计可以保证热网、电网的耦合关系,避免独立网络状态估计出现的耦合处数据不匹配的问题。同时,热电耦合约束使状态估计中的约束增加,相应减少了变量,这使得量测数据的冗余度增加,可以提高状态估计的准确性。因此,必须对热电耦合网络进行联合的状态估计。
[0004] 除此之外,对于多个供能系统耦合的网络而言,各子系统的运行、控制实际由该子系统负责,各子系统之间的决策相互独立且仅有部分数据共享,故需要分布式的计算方法保证各系统计算的独立性。对于热电耦合网络,总优化问题可分解为针对各子系统的问题,各子问题同时进行独立求解,并通过边界信息的交互保证热电耦合关系。考虑到不同子问题的求解耗时不同,也可以采用异步计算等思路对各子问题的计算顺序进行调整,进一步改善计算速度。

发明内容

[0005] 本发明的目的是为了解决现有技术中的上述缺陷,提供一种应用于热电耦合网络的异步分布式状态估计方法,该状态估计方法对热电耦合网络进行联合状态估计,与已有研究不同,该状态估计方法在保留了电热子系统之间的耦合关系的同时,将状态估计问题划分为子问题并实现各子问题的独立求解,保证子系统之间的决策独立和信息私密性;同时该估计方法采用了异步计算框架,提高了计算效率,能够较快地为热电耦合网络的运行、控制、分析提供数据。
[0006] 本发明的目的可以通过采取如下技术方案达到:
[0007] 一种应用于热电耦合网络的异步分布式状态估计方法,包括以下步骤:
[0008] S1、获取热电耦合网络的量测数据;
[0009] S2、用节点法建立热网拟动态模型,使用热网量测数据计算热网拟动态模型中的流速相关系数以及温度时滞系数;
[0010] S3、基于量测数据与拟动态模型系数构造热电耦合网络的状态估计目标函数;
[0011] S4、基于热电耦合网络模型,构造目标函数的增广拉格朗日函数,通过交替方向乘子法将增广拉格朗日函数分解为对应热电耦合网络各子系统的子问题进行求解,并使用异步计算框架对各子问题的求解时长进行控制,最终求解得到热电联合网络的运行状态数据;
[0012] S5、输出热电耦合网络的运行状态。
[0013] 进一步地,在步骤S1中,所述的热电耦合网络的量测数据包括当前时刻电网的节点有功注入功率、节点无功注入功率、线路传输有功潮流、线路传输无功潮流及节点电压幅值,以及热网的节点压强、管道流速、节点温度、热源/热负荷的供应/消耗热功率、耗电量;此外,热电耦合网络的量测数据还包括历史时刻的节点温度和管道流速数据。
[0014] 进一步地,在步骤S2中,用节点法建立热网拟动态模型,使用热网量测数据计算热网拟动态模型中的流速相关系数与温度时滞系数,从而建立当前时刻温度与历史时刻的温度、流量之间的关系,并在管道热量损耗方面对温度进行校正,包括以下步骤:
[0015] S2.1)计算拟动态模型中的流速相关系数:
[0016]
[0017]
[0018]
[0019]
[0020] 式中,Rb,t、Sb,t、γb,t和φb,t为流速相关系数,b为管道编号,t和k分别为当前时刻与历史时刻的编号,Δt为每个时刻的时间长度,mP为管道内热水流速(上标P表示管道),ρ为水的密度,A为管道横截面积,L为管道长度;
[0021] S2.2)根据流速相关系数计算当前时刻管道b中不同段的温度时滞系数 (上标TD表示时滞):
[0022]
[0023] S2.3)基于热网拟动态模型的系数,将当前时刻的未考虑温度损耗的管道流出端温度与多个时刻的管道流入端温度的关系用下式表示:
[0024]
[0025] 式中, 为当前时刻的管道b未考虑温度损耗的管道流出端温度(上标out表示流出端,上标*表示不考虑热损耗), 表示时刻管道的流入端温度(上标in表示流入端);
[0026] S2.4)基于未考虑温度损耗的管道流出端温度以及热网拟动态模型系数,计算考虑温度损耗的管道流出端温度
[0027]
[0028] 式中,TA为外界温度,ch为管道热传导系数,Cp为水的比热容。
[0029] 进一步地,在步骤S3中,基于量测数据与热网拟动态模型系数,构造热电耦合网络的状态估计目标函数:
[0030]
[0031] 式中,i和j分别为热网、电网子系统的编号,SDHN和SEPS分别为热网、电网子系统的编号集,项 为某个热网子系统的量测残差加权和,项 为某个电网子系统的量测残差加权和,以上两项分别表示为:
[0032]
[0033]
[0034] 式中,ps、Ppump、m、Φ、T分别代表热网中的节点压强、水泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度变量, 分别代表节点压强、水泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度的量测值,分别代表由量测方程得到的节点压强、水泵耗电量、管道流速、
热源/热负荷的供应/消耗热功率以及节点温度计算值,Wips、Wipump、Wim、WiΦ、WiT分别代表节点压强、水泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度对应的权重值, 分别代表节点压强、水泵耗电量、管道流速、热源/热
负荷的供应/消耗热功率以及节点温度量测项的编号集,i为量测项的编号。
[0035] 式中,V、θ、Pinj、Qinj、Pflow、Qflow分别代表节点电压幅值、节点电压相、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流变量,分别代表节点电压幅值、节点有功注入功率、节点无功注入
功率、线路有功传输潮流、线路无功传输潮流的量测值, 分别
代表由量测方程得到的节点电压幅值、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流计算值,WiV、 分别代表节点电压幅
值、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流对应的权重值, 分别代表节点电压幅值、节点有功注入功率、
节点无功注入功率、线路有功传输潮流、线路无功传输潮流量测项的编号集。
[0036] 在式(9)中,除了满足管道拟动态特性约束(式(1)-(6))和温度损耗(式(7))外,热网各变量需满足节点水流流入流出量平衡约束、温度混合约束、压强-流量约束、节点供热/耗热量约束、以及水泵耗能约束,公式如下:
[0037]
[0038]
[0039]
[0040]
[0041]
[0042]
[0043]
[0044] 式中, 和 分别为节点处对应流入、流出管道的编号集合, 和 分别表示管道流入端、流出端的压强,Kf为压强损耗系数, 为节点温度, 和 分别为热电联产机组供热量和节点热负荷, 和 分别为流经热电联产机组、负荷的热水流速,iS和iLD分别为热电联产机组、负荷节点的编号,Cp为比热容系数,
和 分别表示热电联产机组与负荷处流入端、流出端水流温度, 为流经水泵的水
流速率, 和 分别为水泵流入端、流出端压强, 为水泵的耗电量,iP为水泵
的编号,ρ为水的密度, 为水泵效率系数,cm为热电联产机组热-电转换比例, 为热电联产机组发电量;
[0045] 在式(10)中,电网各变量需满足节点功率平衡约束以及线路传输功率约束:
[0046]
[0047]
[0048] Pij,t=Vi2gij-ViVj(gij cosθij+bij sinθij)       (20)
[0049] Qij,t=-Vi2(bsi+bij)-ViVj(gij sinθij+bij cosθij)       (21)
[0050] 式中,i和j均为电网中节点的编号,t为所在时刻编号,Pi,t和Qi,t分别为节点有功、无功注入功率,Pij,t和Qij,t分别为线路有功、无功传输功率,Vi和θij为节点电压相量和两节点之间电压相角差,Gij和Bij分别表示导纳矩阵中电导、电纳元素,gij和bij分别表示支路电导和电纳,bsi表示对地电纳。
[0051] 热网、电网的耦合约束由水泵耗电公式(式(16))、热电联产机组产电量-供热量关系式(式(17))组成。
[0052] 在步骤S4中,基于热电耦合网络模型,构造目标函数的增广拉格朗日函数,通过交替方向乘子法将增广拉格朗日函数分解为对应热电耦合网络各子系统的子问题进行求解,并使用异步计算框架对各子问题的求解时长进行控制,最终求解得到热电联合网络的运行状态数据,包括以下步骤:
[0053] S4.1)基于热电耦合网络模型中的热电耦合关系,构造目标函数的增广拉格朗日函数:
[0054]
[0055] 式中, 为热网子系统变量, 为电网子系统变量,λ为拉格朗日乘子,β为增广拉格朗日函数的惩罚因子,z为辅助变量,项 代表热网子系统内部约束以及与电网子系统之间的耦合约束(式(1)-(8)、式(11)-(17)),项 代表电网
子系统内部约束以及与热网子系统之间的耦合约束(式(17)-(21))。
[0056] S4.2)对应热电耦合网络中的每一个子系统,将增广拉格朗日函数划分为各个子问题,热网、电网子问题对应的目标函数分别如下:
[0057]
[0058]
[0059] S4.3)对热网、电网子问题分别进行求解,使用异步计算框架对热网子问题的求解时长进行控制。包括以下内容:
[0060] S4.3.1)电网子问题直接求解。
[0061] 热网子问题在求解过程中,需将热网变量分为两类进行计算,第一类为拟动态模型涉及的系数Rb,t、Sb,t、γb,t、φb,t和 (称为复杂变量),第二类为热网其余变量(称为普通变量)。在已知时滞系数的基础上,首先对热网子问题进行求解,并用所得的普通变量解更新拟动态模型系数,即复杂变量;之后,使用更新后的拟动态模型系数构造新的热网子问题目标函数并再一次进行求解,更新拟动态模型系数;复杂变量、普通变量的交替计算将一直进行,直至前后两次解的差值低于某一特定值,视为热网子问题求解完成。
[0062] S4.3.2)若电网、热网子问题同时求解完成,则视为该轮迭代完成,通过子问题边界变量计算辅助变量与拉格朗日乘子,并进入下一轮迭代。
[0063] 若热网子问题的求解早于电网子问题完成,则热网部分暂时不进行计算,待电网子问题完成求解后,更新辅助变量与拉格朗日乘子,并进入下一轮迭代。
[0064] 若电网子问题的求解早于热网子问题完成,则热网在完成当前时刻的热网普通变量计算后停止计算,直接将边界节点数据与电网交换,更新辅助变量与拉格朗日乘子,进入下一轮迭代。
[0065] S4.3.3)当系统变量满足以下要求时,迭代满足收敛条件:
[0066] ‖rk‖<∈res          (25)
[0067] ‖sk‖<∈dual          (26)
[0068] 式中,k为迭代次数,‖rk‖和‖sk‖分别代表本次迭代计算的原始残差、对偶残差,∈res和∈dual为分别对应原始残差、对偶残差的收敛标准。原始残差、对偶残差的计算方法如下:
[0069]
[0070]
[0071] 在步骤S5中,所述输出热电耦合网络的运行状态包括:电网节点电压的相角与幅值、节点有功注入功率、节点无功注入功率、线路传输有功潮流、线路传输无功潮流;热网节点压强、管道流速、节点温度、热源/热负荷的供应/消耗热功率、水泵耗电量。
[0072] 本发明相对于现有技术具有如下的优点及效果:
[0073] 1、本发明首次考虑了热电耦合网络的分布式状态估计,在保留状态估计结果符合热电耦合关系的同时,保证热电耦合网络中各子系统的决策独立性与数据隐私性。
[0074] 2、本发明首次在热电耦合网络的状态估计问题求解过程中使用了异步计算框架,有效地改善了计算效率,缩短了计算时间。
[0075] 3、本发明使用节点法建立了热网的拟动态模型,更贴近实际热网管道内部温度变化动态,减小计算误差;拟动态模型的使用增加了系统的量测量个数,提高系统量测冗余度,进一步改善状态估计准确度。
[0076] 4、本发明将热网变量分为两类进行交替计算,有效降低了计算难度,缩短计算时间。
[0077] 5、本发明可以为热电耦合网络的规划、分析、优化与控制等提供直观、准确的系统运行状态数据,该估计方法计算复杂度低、计算效率高,有较为广阔的应用前景。附图说明
[0078] 图1是本发明公开的种应用于热电耦合网络的异步分布式状态估计方法的流程步骤图。

具体实施方式

[0079] 为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
[0080] 实施例
[0081] 如图1所示,本实施例提供一种应用于热电耦合网络的异步分布式状态估计,基于当前时刻与历史时刻的量测数据,构造各子系统的状态估计问题目标函数与约束,并通过交替方向乘子法实现分布式求解,同时该方法采用了异步计算框架安排子问题求解顺序,最终实现热电耦合网络的异步分布式状态估计。该状态估计方法包括以下步骤:
[0082] S1、获取热电耦合网络的量测数据,包括电网当前时刻的节点有功注入功率、节点无功注入功率、线路传输有功潮流、线路传输无功潮流及节点电压幅值,热网当前时刻的节点压强、管道流速、节点温度、热源/热负荷的供应/消耗热功率、水泵耗电量,以及热网历史时刻的节点温度和管道流速数据。您
[0083] S2、用节点法建立热网拟动态模型,使用热网量测数据计算热网拟动态模型中的流速相关系数以及温度时滞系数,包括以下步骤:
[0084] S2.1)计算拟动态模型中的流速相关系数:
[0085]
[0086]
[0087]
[0088]
[0089] 式中,Rb,t、Sb,t、γb,t和φb,t为流速相关系数,b为管道编号,t和k分别为当前时刻与历史时刻的编号,Δt为每个时刻的时间长度,mP为管道内热水流速(上标P表示管道),ρ为水的密度,A为管道横截面积,L为管道长度。
[0090] S2.2)根据流速相关系数计算温度时滞系数 (上标TD表示时滞):
[0091]
[0092] S2.3)基于拟动态模型的系数,当前时刻的管道流出端温度(未考虑温度损耗)与多个时刻的管道流入端温度的关系用下式表示:
[0093]
[0094] 式中, 为当前时刻的管道b未考虑温度损耗的管道流出端温度(上标out表示流出端,上标*表示不考虑热损耗), 表示时刻管道的流入端温度(上标in表示流入端);
[0095] S2.4)基于未考虑温度损耗的管道流出端温度以及拟动态模型系数,计算考虑温度损耗的管道流出端温度:
[0096]
[0097] 式中,TA为外界温度,ch为管道热传导系数,Cp为水的比热容。
[0098] S3、基于量测数据与拟动态模型系数构造热电耦合网络的状态估计目标函数:
[0099]
[0100] 式中,i和j分别为热网、电网子系统的编号,SDHN和SEPS分别为热网、电网子系统的编号集。项 为某个热网子系统的量测残差加权和,项 为某个电网子系统的量测残差加权和,以上两项可分别表示为:
[0101]
[0102]
[0103] 式中,ps、Ppump、m、Φ、T分别代表热网中的节点压强、水泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度变量, 分别代表节点压强、水泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度的量测值,分别代表由量测方程得到的节点压强、水泵耗电量、管道流速、
热源/热负荷的供应/消耗热功率以及节点温度计算值,Wips、Wipump、Wim、WiΦ、WiT分别代表节点压强、水泵耗电量、管道流速、热源/热负荷的供应/消耗热功率以及节点温度对应的权重值, 分别代表节点压强、水泵耗电量、管道流速、热源/热
负荷的供应/消耗热功率以及节点温度量测项的编号集,i为量测项的编号。
[0104] 式中,V、θ、Pinj、Qinj、Pflow、Qflow分别代表节点电压幅值、节点电压相角、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流变量,分别代表节点电压幅值、节点有功注入功率、节点无功注入
功率、线路有功传输潮流、线路无功传输潮流的量测值, 分别
代表由量测方程得到的节点电压幅值、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流计算值,WiV、 分别代表节点电压幅
值、节点有功注入功率、节点无功注入功率、线路有功传输潮流、线路无功传输潮流对应的权重值, 分别代表节点电压幅值、节点有功注入功率、
节点无功注入功率、线路有功传输潮流、线路无功传输潮流量测项的编号集。
[0105] S4、基于热电耦合网络模型,构造目标函数的增广拉格朗日函数,通过交替方向乘子法将增广拉格朗日函数分解为对应热电耦合网络各子系统的子问题进行求解,并使用异步计算框架对各子问题的求解时长进行控制,最终求解得到热电联合网络的运行状态数据,包括以下步骤:
[0106] S4.1)基于热电耦合网络模型中的热电耦合关系,构造目标函数的增广拉格朗日函数:
[0107]
[0108] 式中, 为热网子系统变量, 为电网子系统变量,λ为拉格朗日乘子,β为增广拉格朗日函数的惩罚因子,z为辅助变量,项 与项 分别代表热网子系统、电网子系统与其他子系统之间的耦合约束。
[0109] S4.2)对应热电耦合网络中的每一个子系统,将增广拉格朗日函数划分为各个子问题,热网、电网子问题对应的目标函数分别如下:
[0110]
[0111]
[0112] S4.3)对热网、电网子问题分别进行求解,使用异步计算框架对热网子问题的求解时长进行控制。包括以下内容:
[0113] S4.3.1)电网子问题直接求解。
[0114] 热网子问题在求解过程中,需将热网变量分为两类进行计算,第一类为拟动态模型涉及的系数Rb,t、Sb,t、γb,t、φb,t和 (称为复杂变量),第二类为热网其余变量(称为普通变量)。在已知时滞系数的基础上,首先对热网子问题进行求解,并用所得的普通变量解更新拟动态模型系数,即复杂变量;之后,使用更新后的拟动态模型系数构造新的热网子问题目标函数并再一次进行求解,更新拟动态模型系数;复杂变量、普通变量的交替计算将一直进行,直至前后两次解的差值低于某一特定值,视为热网子问题求解完成。
[0115] S4.3.2)若电网、热网子问题同时求解完成,则视为该轮迭代完成,通过子问题边界变量计算辅助变量与拉格朗日乘子,并进入下一轮迭代。
[0116] 若热网子问题的求解早于电网子问题完成,则热网部分暂时不进行计算,待电网子问题完成求解后,更新辅助变量与拉格朗日乘子,并进入下一轮迭代。
[0117] 若电网子问题的求解早于热网子问题完成,则热网在完成当前时刻的热网普通变量计算后停止计算,直接将边界节点数据与电网交换,更新辅助变量与拉格朗日乘子,进入下一轮迭代。
[0118] S4.3.3)当系统变量满足以下要求时,迭代满足收敛条件:
[0119] ‖rk‖<∈res
[0120] ‖sk‖<∈dual
[0121] 式中,k为迭代次数,‖rk‖和‖sk‖分别代表本次迭代计算的原始残差、对偶残差,∈res和∈dual为分别对应原始残差、对偶残差的收敛标准。原始残差、对偶残差的计算方法如下:
[0122]
[0123]
[0124] S5、输出热电耦合网络的运行状态,计算完成。输出的运行状态数据包括:电网节点电压的相角与幅值、节点有功注入功率、节点无功注入功率、线路传输有功潮流、线路传输无功潮流;热网节点压强、管道流速、节点温度、热源/热负荷的供应/消耗热功率、水泵耗电量。
[0125] 综上所述,在采用以上方案后,本发明为热电耦合网络的状态估计提供了新的方法,不同于已有的非分布式热电耦合网络状态估计,该方法通过交替方向乘子法对状态估计问题进行分布式求解,并采用异步计算框架对各子问题的求解时长进行控制。该方法能在保证数据准确性的同时,保留热电耦合网络各子系统的决策独立性和数据隐私性,并有效地提高了分布式计算的效率,具有较高的实际应用价值,值得推广。
[0126] 上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈