首页 / 专利库 / 电脑编程 / 算法 / 一种多区域电-气综合能源系统分布式鲁棒优化方法

一种多区域电-气综合能源系统分布式鲁棒优化方法

阅读:196发布:2024-01-04

专利汇可以提供一种多区域电-气综合能源系统分布式鲁棒优化方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种多区域电-气综合 能源 系统分布式鲁棒优化方法,首先建立多区域电-气综合能源系统模型,并设定相应的约束条件;根据不确定变量的上下界对所设定的约束条件作等价转换,消除所述多区域电-气综合能源系统模型中的不确定变量;针对系统模型非凸非线性的特点,采用一个包络Weymouth曲线的四边形区域代替Weymouth曲线,将非凸优化问题转化为凸优化问题,以降低求解复杂度并保证 算法 的收敛性。上述方法提出了基于凸优化的多区域电-气综合能源系统分布式鲁棒优化手段,保证了算法的收敛性并降低了求解复杂度。,下面是一种多区域电-气综合能源系统分布式鲁棒优化方法专利的具体信息内容。

1.一种多区域电-气综合能源系统分布式鲁棒优化方法,其特征在于,所述方法包括:
步骤1、建立多区域电-气综合能源系统模型,并设定相应的约束条件;
步骤2、根据不确定变量的上下界对所设定的约束条件作等价转换,消除所述多区域电-气综合能源系统模型中的不确定变量;
步骤3、针对系统模型非凸非线性的特点,采用一个包络Weymouth曲线的四边形区域代替Weymouth曲线,将非凸优化问题转化为凸优化问题,以降低求解复杂度并保证算法的收敛性。
2.根据权利要求1所述多区域电-气综合能源系统分布式鲁棒优化方法,其特征在于,在步骤1中,所建立的多区域电-气综合能源系统模型表示为:
式中, 表示时段数、常规机组、天然气井、天然气储能的集合;
表示常规机组出、上备用容量、下备用容量;fwt表示气井的产气量; 表示考虑不确定性情况下天然气储能所应对的不确定区间,以流出为正; 表示天然气储能的放气、储气效率; 表示常规机组常数项、一次项、二次项成本系数; 表示常规机组的上备用、下备用成本系数;cw表示气井产气量的成本系数; 表示天然气储能的上备用、下备用成本系数。
3.根据权利要求1所述多区域电-气综合能源系统分布式鲁棒优化方法,其特征在于,在步骤1中,所设定的约束条件包括电力网络约束条件、天然气网络约束条件和区域间耦合约束条件,其中:
电力网络约束条件表示为:
该式2表示不考虑不确定性情况下的全网功率平衡方程约束,其中, 表示燃气机组、可再生能源机组、电力负荷集合;put, 表示常规机组和燃气机组出力、可再生能源出力预测值、电力负荷预测值; 表示与区域i通过电联络线相连的区域集合; 表示区域i向区域j输送的功率;
该式3表示不考虑不确定性情况下的线路功率方程约束,其中, 表示线路功率;
表示线路传输容量; 表示常规机组和燃气机组、可再生能源机组、与区域j间的联络线、电力负荷对支路的功率转移分布因子;
在不确定场景下,可再生能源、负荷的实际功率分别为 为保持
全网功率平衡,基于线性决策规则,常规机组、燃气机组按照一定比例响应可再生能源出力、电力负荷的不确定性,如下式所示:
式中,αut,Δput表示常规机组和燃气机组的参与因子、出力调整量;
αut应满足式4的约束:
在不确定场景下,常规机组和燃气机组的功率约束、爬坡约束如式5和6所示,式中,表示常规机组和燃气机组的出力上限、下限; 表示常规机组和燃气机
组的上、下爬坡率:
在不确定场景下,随着节点注入功率的改变,线路功率也发生变化,其变化量如式7和8所示,式中,Δplet表示线路功率变化量,plet表示考虑不确定性后的线路功率:
4.根据权利要求3所述多区域电-气综合能源系统分布式鲁棒优化方法,其特征在于,所述天然气网络约束条件表示为:
式9和10表示气井产气量约束、节点压强约束,式中,fwt, 表示气井产气量、节点压强平方; 表示气井产量上限、下限; 表示节点压强平方的上限、下限;
该式表示天然气管道的Weymouth方程约束,式中, 表示管道流量, 表示管道常数,表示管道首端、末端压强平方;
式12表示压缩机入口、出口压强约束,式13表示压缩机入口、出口流量约束,压缩机通过消耗天然气实现加压,耗气量与入口流量成正比;式中, 表示压缩机首端、末端压强平方, 表示压缩机入口、出口流量,Zc,λc表示压缩机最大压缩比、耗气量系数;
式14-17表示不考虑不确定性情况下天然气储能约束,其中,式14和15表示天然气储能的储气、放气流量约束,式16表示天然气储能的储气量约束,式17表示天然气储能的初始、终止状态下的储气量约束;式中, sst表示天然气储能的输入流量、输出流量、储气量; 表示天然气储能的储气、放气流量上限; 表示天然气储能的储
气量上限、下限、初始值;T表示终止时段数;
式18表示节点流量平衡约束,式中, 表示与节点ng相连的气井、天
然气储能、天然气负荷; 表示与天然气负荷dg相连的燃气机组集合;
表示以节点ng为首端节点、末端节点的管道集合; 表示以节
点ng为首端节点、末端节点的压缩机集合; 表示区域i通过节点ng相连的区域集合;
表示天然气负荷的预测值; 表示燃气机组的效率;K表示天然气的热值;
在不确定场景下,天然气负荷的偏差量可以表示为:
因此,天然气负荷的实际为 其偏离预测值的部分 由负荷节点配置的
储能装置承担;天然气储能在不确定场景下的约束条件表示为:
式中, 表示与天然气储能s相连的天然气负荷集合; 表示负荷实际值高于预测值, 表示负荷实际值低于预测值。
5.根据权利要求3所述多区域电-气综合能源系统分布式鲁棒优化方法,其特征在于,所述区域间耦合约束条件表示为:
其中,式22表示区域间联络线功率约束,式23区域间联络管道流量约束,式中,表示联络线、联络管道的传输容量。
6.根据权利要求1所述多区域电-气综合能源系统分布式鲁棒优化方法,其特征在于,在步骤2中,根据不确定变量的上下界对所设定的约束条件作等价转换的过程具体为:
将上述式5等价为
式中, 表示常规机组和燃气机组的向上、向下旋转备用容量;
将上述式6等价为
上述式7做如下转化
其中, 为辅助变量;
因此,将上述式7和8等价为
将上述式20等价为:
将上述式21等价为:
式中, 的计算公式如下:
由此得到多区域电-气综合能源系统的鲁棒模型。
7.根据权利要求1所述多区域电-气综合能源系统分布式鲁棒优化方法,其特征在于,所述步骤3的过程具体为:
所采用的四边形区域由四条线段构成,其中,线段AB、CD为Weymouth曲线以A、C点为切点的切线段;线段AD、BC为Weymouth曲线过点A、C的切线段;该四边形区域的数学表达式如下所示:
上述式11由式36代替后,得到的模型记为模型II,原模型记为模型I;
然后根据松弛解得到管道流向,进而消除Weymouth方程中的绝对值符号,即将式11转化为式37,进一步,式37等价为式38和39:
式38为凸的二次约束,式39看做两个凸函数之差的形式;
设 为管道流量初值,由松弛解得到,对式39中的二次项 在初值点处做一阶泰勒展开,由此式39转化为式40,式中 为松弛变量,经过上述处理后得到的模型记为模型III,
然后求解模型III,直至CCP算法收敛,得到单个区域子问题的求解结果。

说明书全文

一种多区域电-气综合能源系统分布式鲁棒优化方法

技术领域

[0001] 本发明涉及电气系统技术领域,尤其涉及一种多区域电-气综合能源系统分布式鲁棒优化方法。

背景技术

[0002] 随着新能源机组的大规模接入,其随机性和波动性对电系统运行带来了严峻的挑战,电力系统的运行灵活性亟待提高,鲁棒优化能够应对不确定区间的最坏情况,保证系统的安全稳定运行,因而广泛用于电力系统的优化运行。燃气轮机具有调节速度快、调节范围广的优点,常作为调频机组使用,以应对电力系统不确定性。由于燃气轮机和电转气装置的广泛应用,电力系统与天然气系统的耦合程度进一步加深,进而形成电-气综合能源系统。近年来,电-气综合能源系统成为能源领域的研究热点。
[0003] 区域电-气综合能源系统往往独立运行,彼此间缺乏协调,不利于系统整体的经济性和安全性。多区域电-气综合能源系统的协同运行能够整合各个区域的资源,增强系统的经济性和安全性,以及系统应对不确定性的能力。多区域系统的管理方式分为集中式和分布式,集中式管理采用单一决策单元对所有区域统一管理,分布式管理即每个决策单元仅负责各自区域内的能量管理,通过彼此间交换信息实现区域间的协同。与集中式管理方法相比,多区域综合能源系统的分布式管理能够降低计算和通信成本,并保障用户的隐私信息。但现有技术方案中,分布式优化算法对模型的要求较高,当子问题均为凸优化问题时,其收敛性方可保证,然而两阶段鲁棒模型、天然气系统模型呈现非凸非线性的特点,不利于分布式算法的实现。

发明内容

[0004] 本发明的目的是提供一种多区域电-气综合能源系统分布式鲁棒优化方法,该方法提出了基于凸优化的多区域电-气综合能源系统分布式鲁棒优化手段,保证了算法的收敛性并降低了求解复杂度。
[0005] 本发明的目的是通过以下技术方案实现的:
[0006] 一种多区域电-气综合能源系统分布式鲁棒优化方法,所述方法包括:
[0007] 步骤1、建立多区域电-气综合能源系统模型,并设定相应的约束条件;
[0008] 步骤2、根据不确定变量的上下界对所设定的约束条件作等价转换,消除所述多区域电-气综合能源系统模型中的不确定变量;
[0009] 步骤3、针对系统模型非凸非线性的特点,采用一个包络Weymouth曲线的四边形区域代替Weymouth曲线,将非凸优化问题转化为凸优化问题,以降低求解复杂度并保证算法的收敛性。
[0010] 由上述本发明提供的技术方案可以看出,上述方法提出了基于凸优化的多区域电-气综合能源系统分布式鲁棒优化手段,保证了算法的收敛性并降低了求解复杂度。附图说明
[0011] 为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
[0012] 图1为本发明实施例提供的多区域电-气综合能源系统分布式鲁棒优化方法流程示意图;
[0013] 图2为本发明实施例所述四边形松弛的示意图;
[0014] 图3为本发明所举实例基于两区域电-气综合能源系统的拓扑结构示意图;
[0015] 图4为本发明所举实例收敛过程曲线示意图;
[0016] 图5为本发明所举实例不同运行方式下总成本随不确定区间的变化曲线示意图。

具体实施方式

[0017] 下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
[0018] 下面将结合附图对本发明实施例作进一步地详细描述,如图1所示为本发明实施例提供的多区域电-气综合能源系统分布式鲁棒优化方法流程示意图,所述方法包括:
[0019] 步骤1、建立多区域电-气综合能源系统模型,并设定相应的约束条件;
[0020] 在该步骤中,所建立的多区域电-气综合能源系统模型表示为:
[0021]
[0022] 式中, 表示时段数、常规机组、天然气井、天然气储能的集合;表示常规机组出力、上备用容量、下备用容量;fwt表示气井的产气量;
表示考虑不确定性情况下天然气储能所应对的不确定区间,以流出为正;
表示天然气储能的放气、储气效率; 表示常规机组常数项、一次项、二次项成本系数; 表示常规机组的上备用、下备用成本系数;cw表示气井产气量的成本系数; 表示天然气储能的上备用、下备用成本系数。
[0023] 所设定的约束条件包括电力网络约束条件、天然气网络约束条件和区域间耦合约束条件,具体为:
[0024] 1、电力网络约束条件
[0025] 首先考虑可再生能源出力、电力负荷的不确定性,其偏离预测值的偏差量表示如下:
[0026]
[0027]
[0028] 式中, 表示可再生能源机组、负荷的集合; 表示可再生能源出力、电力负荷的偏差量; 表示可再生能源出力偏差量上下限; 表示电力负荷偏差量上下限。
[0029] 则电力网络约束条件表示为:
[0030]
[0031] 该式2表示不考虑不确定性情况下的全网功率平衡方程约束,其中, 表示燃气机组、可再生能源机组、电力负荷集合;put,PuRt,Pdet表示常规机组和燃气机组出力、可再生能源出力预测值、电力负荷预测值; 表示与区域i通过电联络线相连的区域集合;表示区域i向区域j输送的功率;
[0032]
[0033] 该式3表示不考虑不确定性情况下的线路功率方程约束,其中, 表示线路功率;表示线路传输容量; 表示常规机组和燃气机组、可再生能源机组、与
区域j间的联络线、电力负荷对支路的功率转移分布因子(GSDF);
[0034] 在不确定场景下,可再生能源、负荷的实际功率分别为 为保持全网功率平衡,基于线性决策规则,常规机组、燃气机组按照一定比例响应可再生能源出力、电力负荷的不确定性,如下式所示:
[0035]
[0036] 式中,αut,Δput表示常规机组和燃气机组的参与因子、出力调整量;
[0037] αut应满足式4的约束:
[0038]
[0039] 在不确定场景下,常规机组和燃气机组的功率约束、爬坡约束如式5和6所示,式中, 表示常规机组和燃气机组的出力上限、下限; 表示常规机组和燃气机组的上、下爬坡率:
[0040]
[0041]
[0042] 在不确定场景下,随着节点注入功率的改变,线路功率也发生变化,其变化量如式7和8所示,式中, 表示线路功率变化量, 表示考虑不确定性后的线路功率:
[0043]
[0044]
[0045] 2、天然气网络约束条件:
[0046] 首先考虑天然气负荷的不确定性,其偏离预测值的偏差量表示如下:
[0047]
[0048] 式中, 表示天然气负荷的集合, 表示天然气负荷的偏差量, 表示天然气负荷偏差量上、下限。
[0049] 该天然气网络约束条件表示为:
[0050]
[0051] 式9和10表示气井产气量约束、节点压强约束,式中,fwt, 表示气井产气量、节点压强平方; 表示气井产量上限、下限; 表示节点压强平方的上限、下限;
[0052]
[0053] 该式表示天然气管道的Weymouth方程约束,式中, 表示管道流量, 表示管道常数, 表示管道首端、末端压强平方;
[0054]
[0055]
[0056] 式12表示压缩机入口、出口压强约束,式13表示压缩机入口、出口流量约束,压缩机通过消耗天然气实现加压,耗气量与入口流量成正比;式中, 表示压缩机首端、末端压强平方, 表示压缩机入口、出口流量,Zc,λc表示压缩机最大压缩比、耗气量系数;
[0057]
[0058]
[0059]
[0060]
[0061] 式14-17表示不考虑不确定性情况下天然气储能约束,其中,式14和15表示天然气储能的储气、放气流量约束,式16表示天然气储能的储气量约束,式17表示天然气储能的初始、终止状态下的储气量约束;式中, sst表示天然气储能的输入流量、输出流量、储气量; 表示天然气储能的储气、放气流量上限; 表示天然气储能的储气量上限、下限、初始值;T表示终止时段数;
[0062]
[0063] 式18表示节点流量平衡约束,式中, 表示与节点ng相连的气井、天然气储能、天然气负荷; 表示与天然气负荷dg相连的燃气机组集合;
[0064] 表示以节点ng为首端节点、末端节点的管道集合; 表示以节点ng为首端节点、末端节点的压缩机集合; 表示区域i通过节点ng相连的区域集合; 表示天然气负荷的预测值; 表示燃气机组的效率;K表示天然气的热值;
[0065] 在不确定场景下,由于燃气机组参与应对电力负荷的不确定性,其耗气量将跟随其出力的变化而变化,进而将电力负荷的不确定性转化为天然气负荷的不确定性,因此,天然气负荷的偏差量可以表示为:
[0066]
[0067] 因此,天然气负荷的实际为 其偏离预测值的部分 由负荷节点配置的储能装置承担;天然气储能在不确定场景下的约束条件表示为:
[0068]
[0069]
[0070] 式中, 表示与天然气储能s相连的天然气负荷集合; 表示负荷实际值高于预测值, 表示负荷实际值低于预测值。
[0071] 3、所述区域间耦合约束条件表示为:
[0072]
[0073]
[0074] 其中,式22表示区域间联络线功率约束,式23区域间联络管道流量约束,式中,表示联络线、联络管道的传输容量。
[0075] 步骤2、根据不确定变量的上下界对所设定的约束条件作等价转换,消除所述多区域电-气综合能源系统模型中的不确定变量;
[0076] 在该步骤中,上述式(44)-(47)以及(58)-(60)中均含有不确定变量,无法直接进行求解,按照鲁棒优化的思想,在最恶劣情况下满足约束条件,则在任意不确定条件下均满足约束条件。
[0077] 上述根据不确定变量的上下界对所设定的约束条件作等价转换的过程具体为:将上述式5等价为
[0078]
[0079]
[0080]
[0081] 式中, 表示常规机组和燃气机组的向上、向下旋转备用容量;将上述式6等价为
[0082]
[0083]
[0084] 上述式7做如下转化
[0085]
[0086] 其中, 为辅助变量;
[0087]
[0088] 因此,将上述式7和8等价为
[0089]
[0090] 将上述式20等价为:
[0091]
[0092]
[0093] 将上述式21等价为:
[0094]
[0095]
[0096] 式中, 的计算公式如下:
[0097]
[0098]
[0099] 综上所示,式(44)-(47)以及(58)-(60)可以由式(63)-(74)等价替换,即转化为确定性的优化问题,由此得到多区域电-气综合能源系统的鲁棒模型。
[0100] 步骤3、针对系统模型非凸非线性的特点,采用一个包络Weymouth曲线的四边形区域代替Weymouth曲线,将非凸优化问题转化为凸优化问题,以降低求解复杂度并保证算法的收敛性。
[0101] 在该步骤中,所采用的四边形区域由四条线段构成,如图2所示为本发明实施例所述四边形松弛的示意图,图2中:线段AB、CD为Weymouth曲线以A、C点为切点的切线段;线段AD、BC为Weymouth曲线过点A、C的切线段;该四边形区域的数学表达式如下所示:
[0102]
[0103] 上述式11由式36代替后,得到的模型记为模型II,原模型记为模型I;
[0104] 然后根据松弛解得到管道流向,进而消除Weymouth方程中的绝对值符号,即将式11转化为式37,进一步,式37等价为式38和39:
[0105]
[0106]
[0107]
[0108] 式38为凸的二次约束,式39看做两个凸函数之差的形式;
[0109] 设 为管道流量初值,由松弛解得到,对式39中的二次项 在初值点处做一阶泰勒展开,由此式39转化为式40,式中 为松弛变量,经过上述处理后得到的模型记为模型III,
[0110]
[0111] 然后求解模型III,直至CCP算法收敛,得到单个区域子问题的求解结果。
[0112] 另外,在上述步骤的基础上,还可以调用ADMM算法迭代求解各个区域子问题,当ADMM算法收敛时,即可实现各区域间的协同。
[0113] 本实施例中,ADMM算法嵌套调用CCP算法,ADMM算法的过程为,上述多区域电-气能源系统分布式鲁棒优化问题可以抽象为如下问题:
[0114]
[0115] 式中,xi为区域i的决策变量;Ωi表示决策变量xi的可行域;Cost(xi)表示区域i的目标函数;Ae,Βe,ce,Ag,Βg,cg为常数矩阵,可以由式(61)和(62)得到;εpri,εdual表示原残差和对偶残差,可以由εabs,εrel计算得到。
[0116] 上述技术方案针对多区域电-气综合能源系统,各个区域之间通过联络线、联络管道相连,其功率、流量均可控。每个区域电-气综合能源系统由电力网络、天然气网络构成,并通过燃气机组耦合。实际运行中,各区域可以协同运行,亦可以独立运行。
[0117] 下面以具体的实例对上述方法的过程及结果进行详细说明,本实例基于两区域电-气综合能源系统,区域间通过联络线和联络管道相连,拓扑结构如图3所示。其中区域1为5节点电力系统和7节点天然气系统,区域2为6节点电力系统和6节点天然气系统,测试工具采用Matlab2015b编程软件和GUROBI8.1商用求解器。
[0118] 在本实例中,设置εabs=0.0001、εrel=0.004、de=5×103、dg=1.5×104,得到如图4所示的收敛过程曲线示意图。其中,(a)表示联络线功率、联络管道流量偏差量的变化曲线,(b)表示联络线功率、联络管道流量的变化曲线,(c)表示原残差、对偶残差的变化曲线,(d)表示目标函数的变化曲线。由图4可知,随着迭代次数的增加,联络线功率、联络管道流量逐渐趋于稳定,其相应偏差量逐渐减小,目标函数也逐渐趋于稳定。在前20次迭代中,联络线功率、联络管道流量、目标函数变化较为明显,经过20次迭代后基本趋于稳定;同时联络线功率、联络管道流量偏差量经过20次迭代后趋于稳定,且数量级较小,作为ADMM算法收敛指标的原残差与对偶残差,分别在第68、15次迭代时满足收敛判据。
[0119] 另外,本实例的方案还探究了多区域协同运行对于提升系统应对不确定性能力的影响,在不同的不确定性区间下,比较了协同运行方式(CS)和独立运行方式(IS)下的运行成本,如图5所示为0不同运行方式下总成本随不确定区间的变化曲线示意图,图5中:左边的(a)表示可再生能源出力不确定性为5%、电负荷不确定性变化时的运行成本曲线,右边的(b)表示电负荷不确定性保持5%、可再生能源出力不确定性变化时的运行成本曲线,(a)、(b)两图中天然气负荷的不确定性均为1.5%。由图可知,协同运行方式下的成本低于独立运行方式下的成本,且成本随不确定区间的扩大而增大。此外,运行成本变化趋势并非线性,主要由于不同机组的上下备用成本不同。
[0120] 由于系统应对不确定的能力有限,当不确定性增大到一定值,优化模型无可行解,此时必须采取弃弃光或切负荷的措施,如下表1所示,当两个区域在独立运行方式下均不存在可行解时,通过协同运行方式仍可以找到可行解,由此可以得出多区域系统协同运行可以提升系统应对不确定性的能力。
[0121] 表1不同不确定区间下优化模型的可行性(一)
[0122]
[0123] 另外,由于所提出的四边形松弛法与现有技术仅在天然气网络模型上有所差别,故以单一天然气网络为测试对象测试不同方法的计算性能。测试结果如表2所示,Gas7、Gas20分别表示7节点天然气网络和20节点天然气网络,误差表示Weymouth方程两端的差值,表征方法的精确度。
[0124] 表2所提方法与现有技术的计算性能对比
[0125]
[0126]
[0127] 值得注意的是,本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
[0128] 以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。
高效检索全球专利

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

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

申请试用

分析报告

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

申请试用

QQ群二维码
意见反馈