首页 > 专利 > 杭州电子科技大学 > 一种考虑配电系统负荷波动的日前鲁棒调度方法专利详情

一种考虑配电系统负荷波动的日前鲁棒调度方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2019-09-30
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2020-01-24
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2021-04-13
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2039-09-30
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201910942849.3 申请日 2019-09-30
公开/公告号 CN110635496B 公开/公告日 2021-04-13
授权日 2021-04-13 预估到期日 2039-09-30
申请年 2019年 公开/公告年 2021年
缴费截止日
分类号 H02J3/28H02J15/00 主分类号 H02J3/28
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 0
权利要求数量 1 非专利引证数量 1
引用专利数量 0 被引证专利数量 0
非专利引证 1、CN 103904664 A,2014.07.02CN 103259314 A,2013.08.21CN 110084465 A,2019.08.02CN 107742901 A,2018.02.27CN 105870950 A,2016.08.17杨震 等.带有储能的可再生能源调度问题研究《.能源互联网专栏》.2017,第45卷(第7期),吴霜 等.含分布式储能的配电网多目标运行优化策略研究《.电力工程技术》.2018,第37卷(第2期),吴晨曦 等.考虑电价影响的电动汽车削峰填谷水平评价《.电力系统保护与控制》.2019,第47卷(第17期),;
引用专利 被引证专利
专利权维持 3 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 吴晨曦、何章露、陈泽昊 第一发明人 吴晨曦
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 3
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
杨舟涛
摘要
本发明公开了一种考虑配电系统负荷波动的日前鲁棒调度方法,包括:步骤一:建立先进绝热压缩空气储能AA‑CAES电站运行与风力发电系统发电功率模型,确定模型的约束条件。步骤二:建立目标函数,并采用鲁棒动态规划方法求解模型。步骤三:根据步骤一和步骤二将所提方法在IEEE‑RBTS测试系统进行仿真验证,并对只考虑风力发电机的不确定性和考虑风机及负荷不确定性两种场景进行仿真。本发明提出基于风机和负荷预测误差的CAES鲁棒调度方法,使整个系统的负荷波动更小。本发明将功率不足量与不确定度集的不同权重进行比较,分析系统的预测误差对可靠性的影响。
  • 摘要附图
    一种考虑配电系统负荷波动的日前鲁棒调度方法
  • 说明书附图:图1
    一种考虑配电系统负荷波动的日前鲁棒调度方法
  • 说明书附图:图2
    一种考虑配电系统负荷波动的日前鲁棒调度方法
  • 说明书附图:图3
    一种考虑配电系统负荷波动的日前鲁棒调度方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2021-04-13 授权
2 2020-01-24 实质审查的生效 IPC(主分类): H02J 3/28 专利申请号: 201910942849.3 申请日: 2019.09.30
3 2019-12-31 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种考虑配电系统负荷波动的日前鲁棒调度方法,其特征在于该方法的具体步骤是:
步骤一:建立先进绝热压缩空气储能AA‑CAES电站运行与风力发电系统发电功率模型,确定其模型的约束条件为:
式(1)表示AA‑CAES电站压缩功率约束;其中,Pc,t表示时段t的压缩功率;ηc表示压缩过程效率; 表示时段t流入压缩机的流量;γ表示空气比热容比;Rg表示理想气体常数;nc表示压缩机级数;Tc,l,in和 分别表示进入第l级压缩机和末级压缩机的空气温度;βc,l和分别表示第l级压缩机的额定压缩比和末级压缩机在时段t的压缩比,其与储气室气压有关;
Pc,minνc,t≤Pc,t≤Pc,maxνc,t                     (2)
式(2)表示压缩功率上下限约束;其中,Pc,min和Pc,max分别表示压缩功率的上下限;νc,t为二进制决策变量,用于表示AA‑CAES电站是否处于压缩工况,当AA‑CAES电站处于压缩工况时,νc,t=1,反之νc,t=0;
式(3)表示膨胀机发电功率约束;其中,Pg,t表示时段t的发电功率;ηg表示发电过程效率; 表示时段t流入膨胀机的流量;Tg,j,in,t表示t时段进入第j级膨胀机的空气温度,其与储气室气温相关;ng表示膨胀级数;βg,j表示第j级膨胀机额定膨胀比;
Pg,minνg,t≤Pg,t≤Pg,maxνg,t                     (4)
式(4)表示发电功率上下限约束;其中,Pg,min和Pg,max分别表示发电功率的上下限;νg,t为二进制决策变量,用于表示AA‑CAES电站是否处于发电工况,当AA‑CAES电站处于发电工况时,νg,t=1,反之νg,t=0;
νc,tνg,t=0                          (5)
式(5)表示AA‑CAES电站运行工况约束;该约束用于保证AA‑CAES电站不同时工作在压缩工况和发电工况;
AFR(SOCt)=‑0.2985(SOCt‑0.33)+1.85                 (7)
AFR(Pg,t)=‑1.3143(Pg,t/Pg,max‑0.3)+2.3                 (8)
式(6)表示CAES储气室约束;其中,SOCt和SOCt+1分别表示储气室在t和t+1时刻的状态;
CAmax表示储气室总容量;式(7)、(8)分别为空气流量,AFR与SOCt、Pg,t之间的线性函数;
SOCmin≤SOCt≤SOCmax                     (9)
式(9)表示储气室的储气室上下限约束;其中,SOCmin和SOCmax分别表示储气室SOC的最小值和最大值;
νc,t(Pc,t+1‑Pc,t)+νg,t(Pg,t+1‑Pg,t)≥‑νc,tRc,down‑νg,tRg,downνc,t(Pc,t+1‑Pc,t)+νg,t(Pg,t+1‑Pg,t)≤νc,tRc,up+νg,tRg,up   (10)式(10)表示CAES储能电站的爬坡率约束;其中,Rc,up、Rc,down分别表示CAES在压缩工况下功率的上升、下降速率;Rg,up、Rg,down分别表示CAES在膨胀工况下功率的上升、下降速率;
式(11)表示忽略风机的机械损耗及风机间发电功率的相关性,双馈风力发电机的发电
3
功率;其中,PW表示风机捕获的风能功率;ρ表示空气密度,单位为kg/m ;S表示风轮扫风面
2
积,单位为m;ν表示风速,Cp表示风能利用系数;
式(12)表示在负荷水平较低时削减后的双馈风力发电机的发电功率;其中,PW(t)表示t时刻风机的实际值,ζ表示削减系数,取值范围为10%‑20%;
0≤αW(t)≤1
0≤κP(t)≤1
ΓW+ΓL≤Γ
式(13)表示风力发电机发电功率和负荷的不确定集;其中,PL(t)表示t时刻负荷的实际值; 分别表示t时刻风机、负荷的预测值; 分别表示t时刻风机、负
荷的下限值; 分别表示t时刻风机、负荷的上限值;αW(t)、κP(t)分别表示t时刻风机、负荷的权重;ΓW、ΓL分别表示风机、负荷的不确定集;ΓT表示时间集;Γ表示总的不确定集;
步骤二:建立模型的目标函数,并采用鲁棒动态规划方法求解模型;
式(14)表示在T为时间周期以负荷波动最小为目标函数;其中,Pav是等效平均荷载;
式(15)表示等效平均荷载的定义;
式(14)中的优化函数是二次函数,利用鲁棒等效模型将式(14)转换为确定性二次规划问题;CAES的最优控制是一个多步骤决策过程的优化问题,通过利用不同阶段之间的关系,利用动态规划DP求解;从储气室的初始状态来看,每个时间步的CAES控制过程都是一个DP状态;初始状态是没有CAES的初始等效负载;在CAES压缩/膨胀控制的初始过程中,负载向量由[PL]变为[PL]+[PCAES(1)],第t小时后的递推方程为:
其中PCAES(t)为第t小时AA‑CAES储能电站的充放电功率,充电为正,且PCAES(t)=Pc,t,发电为负,且PCAES(t)=Pg,t;υn为第t小时的决策变量;u(t)是由PCAES(t)和容量约束确定的允许决策集合;υn(PCAES(t),u(t))为等效负荷平方和;f(P(t‑1))为第t‑1阶段式(14)的目标函数,f(P(t))为式(14)第t阶段的目标函数;
步骤三:将所提方法在IEEE‑RBTS测试系统进行对系统鲁棒性的最优日前调度方案进行了仿真验证,对只考虑风力发电机的不确定性和考虑风机及负荷不确定性两种场景进行仿真。
说明书

技术领域

[0001] 本发明设计电力系统运行领域,涉及一种考虑配电系统负荷波动的日前鲁棒调度方法。技术背景
[0002] 近年来,节能和减少碳排放导致可再生能源的快速发展。随着风力发电(wind power generation,WPG)等可再生能源(renewable energy resource,RES)的渗透率逐步提高,能量存储系统(energy storage system,ESS)已用于应对RES引入的挑战。压缩空气储能(compressed air energy storage,CAES)是一种大规模的物理储能技术,为电力负荷的调峰和削峰填谷提供了新的解决方案,并缓解了可再生能源过剩的问题。相关研究已经[1‑2]进行 。CAES使用高压,等温空气压缩机/膨胀机,避免了碳排放,相比使用化石燃料的传[3]
统系统效率更高 。CAES具有大规模存储容量,高爬坡率和短启动时间,可以减轻可再生能[4‑6]
源发电的波动。目前,已成为可再生能源储能的良好选择 。
[0003] 不确定性是WPG的主要问题,这对WPG参与电力市场带来了巨大挑战。WPG的不确定[7]性随着预测时间的增加而增加 。在具有不确定性的调度问题中采用随机优化和鲁棒优[8]
化。随机优化方法无法得到精确解。随着场景的增加,计算量也将大幅增加 。鲁棒优化已成为在不确定环境下实现具有合理实用性,经济性和可靠性的解决方案的有效决策工具。
文献[9]提出基于风险条件值(conditional value at risk,CVaR)的分布式鲁棒方程来获得风电储备需求。文献[10]通过将随机模型转化为可以有效求解的确定性双线性矩阵不等式问题,消除了调度模型中风电预测误差的随机变量。鲁棒最优调度(ROD)已被用于可再生[11‑12]
发电系统 。
[0004] 文献[13]结果表明,混合风力发电系统可以在波动的风速条件下提供平稳的功率输出。采用CAES的小型风力涡轮机的建模和实验研究在[13]中提出。该系统在具有1kW [14‑15]WPG‑CAES系统的实验室中进行测试。数百兆瓦的CAES系统也已开展研究 。文献[16]中使用'PLEXOS'对电力系统运行的经济效益以及WPG‑CAES系统的收益和总发电成本进行了评估。文献[17]为参与能源旋转和非旋转备用市场的CAES设施提出了一种自我调度方法。
文献[18]通过双层规划以确定性方法优化了分布式发电和CAES在孤岛微网的发电规模。此外,所提出的模型考虑了旋转储备以响应负载和可再生能源输出的不确定性。遗传算[19] [20] [21]
法 、粒子群算法 、神经网络 等智能算法常用于经济调度问题。智能算法的一个缺点是很难收敛到全局优化解决方案。动态规划(Dynamic programming,DP)方法具有获得全局[22]
最优控制策略的能力。DP可以处理不连续和非线性约束 。分布式动态规划在文献[23]中用于传统发电系统的经济调度。因此,本发明采用鲁棒动态规划解决WPG‑CAES系统中的日前调度问题。
[0005] 参考文献
[0006] [1]Luo,X.,Wang,J.,Krupke,C.,et al.:‘Modelling study,efficiency analysis and optimisation of large‑scale Adiabatic Compressed Air Energy Storage systems with low‑temperature thermal storage’,Applied Energy,2016,162:589‑600.
[0007] [2]Colthrope,A.,‘Canadian firms NRStor and Hydrostor partner up on utility‑scale adiabatic compressed air storage’,https://www.energy‑storage.news/news/canadian‑firms‑nrstor‑and‑hyd rostor‑partner‑up‑on‑utility‑scale‑adiabatic‑c,accessed 12December 2017.
[0008] [3]Perry,Y.L.,Mohsen,S.:‘An approach to reduce the flow requirement for a liquid piston near‑isothermal air compressor/expander in a compressed air energy storage system’,IET Renewable Power Generation,2015,10,(10),pp.1506‑1514.
[0009] [4]National Renewable Energy Laboratory.:‘Cost and performance data for power generation technologies’,NREL technology report,2012.
[0010] [5]Klumpp,F.:‘Comparison of pumped hydro,hydrogen storage and compressed air energy storage for integrating high shares of renewable energies‑Potential,cost‑comparison and ranking’,J.Energy Storage,2016,8,pp.119‑128.
[0011] [6]Shafiee,S.,Zareipour,H.,Knight,A.M.,Amjady,N.,and Mohammadi‑Ivatloo,B.:‘Risk‑constrained bidding and offering strategy for a merchant compressed air energy storage plant’,IEEE Transactions on Power System,2017,32,(2),pp.946‑957.
[0012] [7]Meng,K.,Yang,H.,Dong,Z.,Guo,W.,Wen,F.,Xu,Z.:‘Flexible operational planning framework considering multiple wind energy forecasting service providers’,IEEE Transactions on Sustainable Energy,2016,7,(2),pp.708‑717[0013] [8]CONEJO,A.J., M.,MORALES,J.M.:‘Decision making under uncertainty in electricity markets’,USA:Springer,2010.
[0014] [9]Wang,Z.,Bian,Q.,Xin,H.and Gan,D.:‘A distributionally robust co‑ordinated reserve scheduling model considering CVaR‑based wind power reserve requirements’,IEEE Transactions on Sustainable Energy,2016,7,(2),pp.625‑636.[0015] [10]Li,Y.,Li,Z.,Wen,F.,Mohammad,S.:‘Minimax‑regret robust co‑optimization for enhancing the resilience of integrated power distribution and natural gas system’,accepted for publication by IEEE Transactions on Sustainable Energy.
[0016] [11]Peng,C.,Xie,P.,Pan,L.and Yu,R.:‘Flexible Robust Optimization Dispatch for HybridWind/Photovoltaic/Hydro/Thermal Power System’,IEEE Transactions on Smart Grid,2016,7,(2),pp.751‑762.
[0017] [12]Wei,W.,Liu,F.and Mei,S.:‘Distributionally Robust Co‑Optimization of Energy and Reserve Dispatch’,IEEE Transactions on Sustainable Energy,2016,7,(1),pp.289‑300.
[0018] [13]Christopher,K.,Wang,J.,Jonathan,C.and Lou,X.:‘Modeling and Experimental Study of a Wind Turbine System in Hybrid Connection With Compressed Air Energy Storage’,IEEE Transactions on Energy Conversion,2017,32,(1),pp.137‑145.
[0019] [14]Evans,A.,Strezov,V.and Evans,T.J.:‘Assessment of utility energy storage options for increased renewable energy penetration’,Renewable&Sustainable Energy Reviews,2012,16,(6),pp.4141–4147.
[0020] [15]‘A.CAES.(2013).APEX CAES’,http://www.apexcaes.com,accessed 31 October 2013.
[0021] [16]Brendan,C.,Aidan,D.,Alan,O.,Michael,C.and Vasilis F.:‘Assessing the Economic Benefits of Compressed Air Energy Storage for Mitigating Wind Curtailment’,IEEE Transactions on Sustainable Energy.2015,6,(3),pp.1021‑1028.[0022] [17]Soroush,S.,Hamidreza,Z.and Andrew,M.K.:‘Considering Thermodynamic Characteristics of a CAES Facility in Self‑Scheduling in Energy and Reserve Markets’,IEEE Transactions on Smart Grid,2018,9,(4),pp.3476‑3485.[0023] [18]Zhang,J.,Li,K.,Wang,M.,et al.:‘A Bi‑Level Program for the Planning of an Islanded Microgrid Including CAES’,IEEE Transactions on Industry Applications.2016,52,(4),pp.2768‑2777.
[0024] [19]Damousis,I.G.,Bakirtzis,A.G.and Dokopoulos.P.S.:‘Network constrained economic dispatch using real‑coded genetic algorithm’,IEEE Transactions on Power System,2003,18,(1),pp.198‑205.
[0025] [20]Selvakumar,A.I.and Thanushkodi,K.:‘A new particle swarm optimization solution to nonconvex economic dispatch problems’,IEEE Transactions on Power System,2007,22,(1),pp.42–51.
[0026] [21]Umamaheswari,K.and Devarajan,N.:‘Hybrid weighted probabilistic neural network and biogeography based optimization for dynamic economic dispatch of integrated multiple‑fuel and wind power plants’,International Journal of Electrical Power&Energy Systems,2016,(77),pp.385‑394.
[0027] [22]Zhang,S.,Xiong,R.:‘Adaptive energy management of a plug‑in hybrid electric vehicle  based  on driving  pattern recognition and dynamic programming’,Applied Energy,2015,(155),pp.68‑78.
[0028] [23]Xu,Y.,Zhang,W.and Liu,W.:‘Distributed Dynamic Programming‑Based Approach for Economic Dispatch in Smart Grids’,IEEE Transactions on Industrial Informatics,2015,11,(1),pp.166‑175.

发明内容

[0029] 本发明所要解决的问题是考虑CAES运行特性、风力发电运行特性及不确定集的前提下,提出一种考虑配电系统负荷波动的日前鲁棒调度方法。
[0030] 本发明解决其技术问题所采用的技术方案是:考虑配电系统负荷波动的压缩空气储能的日前鲁棒调度方法。包括如下步骤:
[0031] 步骤一:建立先进绝热压缩空气储能AA‑CAES电站运行与风力发电系统发电功率模型,确定其模型的约束条件为:
[0032]
[0033] 式(1)表示AA‑CAES电站压缩功率约束。其中,Pc,t表示时段t的压缩功率;ηc表示压缩过程效率; 表示时段t流入压缩机的流量;γ表示空气比热容比;Rg表示理想气体常数;nc表示压缩机级数;Tc,l,in和Tc,nc,in分别表示进入第l级压缩机和末级压缩机的空气温度;
βc,l和βc,nc,t分别表示第l级压缩机的额定压缩比和末级压缩机在时段t的压缩比,其与储气室气压有关。
[0034] Pc,minνc,t≤Pc,t≤Pc,maxνc,t   (2)
[0035] 式(2)表示压缩功率上下限约束。其中,Pc,min和Pc,max分别表示压缩功率的上下限;νc,t为二进制决策变量,用于表示AA‑CAES电站是否处于压缩工况,当AA‑CAES电站处于压缩工况时,νc,t=1,反之νc,t=0。
[0036]
[0037] 式(3)表示膨胀机发电功率约束。其中,Pg,t表示时段t的发电功率;ηg表示发电过程效率; 表示时段t流入膨胀机的流量;Tg,j,in,t表示t时段进入第j级膨胀机的空气温度,其与储气室气温相关;ng表示膨胀级数;βg,j表示第j级膨胀机额定膨胀比。
[0038] Pg,minνg,t≤Pg,t≤Pg,maxνg,t   (4)
[0039] 式(4)表示发电功率上下限约束。其中,Pg,min和Pg,max分别表示发电功率的上下限;νg,t为二进制决策变量,用于表示AA‑CAES电站是否处于发电工况,当AA‑CAES电站处于发电工况时,νg,t=1,反之νg,t=0。
[0040] νc,tνg,t=0   (5)
[0041] 式(5)表示AA‑CAES电站运行工况约束。该约束用于保证AA‑CAES电站不同时工作在压缩工况和发电工况。
[0042]
[0043] AFR(SOCt)=‑0.2985(SOCt‑0.33)+1.85   (7)
[0044] AFR(Pg,t)=‑1.3143(Pg,t/Pg,max‑0.3)+2.3   (8)
[0045] 式(6)表示CAES储气室约束。其中,SOCt和SOCt+1分别表示储气室在t和t+1时刻的状态;CAmax表示储气室总容量;式(7)、(8)分别为空气流量,AFR与SOCt、Pg,t之间的线性函数。
[0046] SOCmin≤SOCt≤SOCmax   (9)
[0047] 式(9)表示储气室的储气室上下限约束。其中,SOCmin和SOCmax分别表示储气室SOC的最小值和最大值。
[0048] νc,t(Pc,t+1‑Pc,t)+νg,t(Pg,t+1‑Pg,t)≥‑νc,tRc,down‑νg,tRg,down[0049] νc,t(Pc,t+1‑Pc,t)+νg,t(Pg,t+1‑Pg,t)≤νc,tRc,up+νg,tRg,up   (10)[0050] 式(10)表示CAES储能电站的爬坡率约束。其中,Rc,up、Rc,down分别表示CAES在压缩工况下功率的上升、下降速率;Rg,up、Rg,down分别表示CAES在膨胀工况下功率的上升、下降速率。
[0051]
[0052] 式(11)表示忽略风机的机械损耗及风机间发电功率的相关性,双馈风力发电机的3
发电功率。其中,PW表示风机捕获的风能功率;ρ表示空气密度,单位为kg/m ;S表示风轮扫风
2
面积,单位为m;ν表示风速;Cp表示风能利用系数。
[0053] PW(t)=(1‑ζ)PW(t)   (12)
[0054] 式(12)表示在负荷水平较低时削减后的双馈风力发电机的发电功率。其中,PW(t)表示t时刻风机的实际值,ζ表示削减系数,取值范围为10%‑20%。
[0055]
[0056] 0≤αW(t)≤1
[0057]
[0058]
[0059] 0≤κP(t)≤1
[0060]
[0061] ΓW+ΓL≤Γ
[0062] 式(13)表示风力发电机发电功率和负荷的不确定集。其中,PL(t)表示t时刻负荷的实际值; 分别表示t时刻风机、负荷的预测值; 分别表示t时刻风机、负荷的下限值; 分别表示t时刻风机、负荷的上限值;αW(t)、κP(t)分别表示t时刻风机、负荷的权重;ΓW、ΓL分别表示风机、负荷的不确定集;ΓT表示时间集;Γ表示总的不确定集。
[0063] 步骤二:建立模型的目标函数,并采用鲁棒动态规划方法求解模型。
[0064]
[0065] 式(14)表示在T为时间周期以负荷波动最小为目标函数。其中,Pav是等效平均荷载。
[0066]
[0067] 式(15)表示等效平均荷载的定义。
[0068] 式(14)中的优化函数是二次函数,利用鲁棒等效模型将式(14)转换为确定性二次规划问题。CAES的最优控制是一个多步骤决策过程的优化问题,通过利用不同阶段之间的关系,利用动态规划DP求解。从储气室的初始状态来看,每个时间步的CAES控制过程都是一个DP状态。初始状态是没有CAES的初始等效负载。在CAES压缩/膨胀控制的初始过程中,负载向量由[PL]变为[PL]+[PCAES(1)],第t小时后的递推方程为:
[0069]
[0070] 其中PCAES(t)为第t小时AA‑CAES储能电站的充放电功率,充电为正,且PCAES(t)=Pc,t,发电为负,且PCAES(t)=Pg,t;υn为第t小时的决策变量;u(t)是由PCAES(t)和容量约束等确定的允许决策集合;υn(PCAES(t),u(t))为等效负荷平方和。f(P(t‑1))为第t‑1阶段式(14)的目标函数,f(P(t))为式(14)第t阶段的目标函数。
[0071] 步骤三:将所提方法在IEEE‑RBTS测试系统进行对系统鲁棒性的最优日前调度方案进行了仿真验证,对只考虑风力发电机的不确定性和考虑风机及负荷不确定性两种场景进行仿真。
[0072] 本发明专利区别于现有的研究工作有如下特点。
[0073] 1)与传统的MPPT风电跟踪方法不同,本文采用的是WPG缩减法;
[0074] 2)本发明将CAES作为风力发电配电网的储能系统。针对风力发电和负荷的预测误差,提出了一种CAES的鲁棒日前最优调度方案;
[0075] 3)本发明对系统鲁棒性的最优日前调度方案进行了仿真,并将功率不足量(energy not supplied,ENS)与不确定度集的不同权重进行比较与分析。

实施方案

[0079] 以下结合附图具体说明本发明方法:
[0080] 步骤一:建立先进绝热压缩空气储能AA‑CAES电站运行与风力发电系统发电功率模型,如图1所示,确定其模型的约束条件为:
[0081]
[0082] 式(1)表示AA‑CAES电站压缩功率约束。其中,Pc,t表示时段t的压缩功率;ηc表示压缩过程效率; 表示时段t流入压缩机的流量;γ表示空气比热容比;Rg表示理想气体常数;nc表示压缩机级数;Tc,l,in和Tc,nc,in分别表示进入第l级压缩机和末级压缩机的空气温度;
βc,l和βc,nc,t分别表示第l级压缩机的额定压缩比和末级压缩机在时段t的压缩比,其与储气室气压有关。
[0083] Pc,minνc,t≤Pc,t≤Pc,maxνc,t   (2)
[0084] 式(2)表示压缩功率上下限约束。其中,Pc,min和Pc,max分别表示压缩功率的上下限;νc,t为二进决策制变量,用于表示AA‑CAES电站是否处于压缩工况,当AA‑CAES电站处于压缩工况时,νc,t=1,反之νc,t=0。
[0085]
[0086] 式(3)表示膨胀机发电功率约束。其中,Pg,t表示时段t的发电功率;ηg表示发电过程效率; 表示时段t流入膨胀机的流量;Tg,j,in,t表示t时段进入第j级膨胀机的空气温度,其与储气室气温相关;ng表示膨胀级数;βg,j表示第j级膨胀机额定膨胀比。
[0087] Pg,minνg,t≤Pg,t≤Pg,maxνg,t   (4)
[0088] 式(4)表示发电功率上下限约束。其中,Pg,min和Pg,max分别表示发电功率的上下限;νg,t为二进制变量,用于表示AA‑CAES电站是否处于发电工况,当AA‑CAES电站处于发电工况时,νg,t=1,反之νg,t=0。
[0089] νc,tνg,t=0   (5)
[0090] 式(5)表示AA‑CAES电站运行工况约束。该约束用于保证AA‑CAES电站不同时工作在压缩工况和发电工况。
[0091]
[0092] AFR(SOCt)=‑0.2985(SOCt‑0.33)+1.85   (7)
[0093] AFR(Pg,t)=‑1.3143(Pg,t/Pg,max‑0.3)+2.3   (8)
[0094] 式(6)表示CAES储气室约束。其中,SOCt和SOCt+1分别表示储气室在t和t+1时刻的状态;CAmax表示储气室总容量;式(7)、(8)分别为空气流量,AFR与SOCt、Pg,t之间的线性函数。
[0095] SOCmin≤SOCt≤SOCmax   (9)
[0096] 式(9)表示储气室的储气室上下限约束。其中,SOCmin和SOCmax分别表示储气室SOC的最小值和最大值。
[0097]
[0098] 式(10)表示CAES储能电站的爬坡率约束。其中,Rc,up、Rc,down分别表示CAES在压缩工况下功率的上升、下降速率;Rg,up、Rg,down分别表示CAES在膨胀工况下功率的上升、下降速率。
[0099]
[0100] 式(11)表示忽略风机的机械损耗及风机间发电功率的相关性,双馈风力发电机的3
发电功率。其中,PW表示风机捕获的风能功率;ρ表示空气密度,单位为kg/m ;S表示风轮扫风
2
面积,单位为m;ν表示风速;Cp表示风能利用系数。
[0101] PW(t)=(1‑ζ)PW(t)
[0102] (12)
[0103] 式(12)表示在负荷水平较低时削减后的双馈风力发电机的发电功率。其中,PW(t)表示t时刻风机的实际值,ζ表示削减系数,取值范围为10%‑20%。
[0104]
[0105] 式(13)表示风力发电机发电功率和负荷的不确定集。其中,PL(t)分别表示t时刻负荷的实际值; 分别表示t时刻风机、负荷的预测值; 分别表示t时刻风机、负荷的下限值; 分别表示t时刻风机、负荷的上限值;αW、κP(t)分别表示t时刻风机、负荷的权重;ΓW、ΓL分别表示风机、负荷的不确定集;ΓT表示时间集;Γ表示总不确定集。
[0106] 步骤二:建立模型的目标函数,并采用鲁棒动态规划方法,如图2所示,求解模型。
[0107]
[0108] 式(14)表示在T为时间周期以负荷波动最小为目标函数。其中,Pav是等效平均荷载。
[0109]
[0110] 式(15)表示等效平均荷载的定义。
[0111] 其中在建立目标函数时,需要对压缩空气储能电站进行优化调度,其表达式为:
[0112]
[0113] 式(16)表示CAES的爬坡率在t时段随风机和负荷的波动而变化。其中 分别表示预测风机功率和预测负荷功率。
[0114] 根据式(2)、式(4)及式(9)可以得到:
[0115] Pg,t(t)+ΔPg,t(t)≤Pg,max   (17‑1)
[0116] Pg,t(t)‑ΔPg,t(t)≥Pg,min   (17‑2)
[0117] Pc,t(t)+ΔPc,t(t)≤Pc,max   (17‑3)
[0118] Pc,t(t)‑ΔPc,t(t)≥Pc,min   (17‑4)
[0119] SOCmin≤SOC(t)≤SOCmax   (17‑5)其中ΔPg,t(t),ΔPc,t(t)分别表示一段时间内的膨胀功率和压缩功率,式(17‑1)和(17‑2)分别表示发电模式下的最大和最小发电功率。式(17‑3)和式(17‑4)分别表示压缩模式下的最大压缩功率和最小压缩功率。
[0120] 将式(16)、式(17)合并后可重写为式(18):
[0121]
[0122]
[0123]
[0124]
[0125]
[0126]
[0127] 式(14)中的优化函数是二次函数,含不确定性变量难以求解。将式(18)转化为鲁棒等效模型进行求解,以(18‑1)为例说明等效模型:
[0128]
[0129] 式(19)表示将不等式约束转化为鲁棒对应等价模型,如图3所示。其中,zt,m,pt,m表示辅助变量,m=1表示与风机相关,m=2表示与负荷相关。
[0130] 鲁棒等效模型将式(14)转换为确定性二次规划问题。CAES的最优控制是一个多步骤决策过程的优化问题。通过利用不同阶段k之间的关系,动态规划适合于解决这一优化问题。从储气室的初始状态来看,每个时间步的CAES控制过程都是一个DP状态。初始状态是没有CAES的初始等效负载。在CAES压缩/膨胀控制的初始状态中,负载向量由[PL]变为[PL]+[PCAES(1)],第t小时后的递推方程为:
[0131]
[0132] 其中PCAES(t)为第t小时AA‑CAES储能电站的充放电功率,充电为正,且PCAES(t)=Pc,t,发电为负,且PCAES(t)=Pg,t;υn为第t小时的决策变量;u(t)是由PCAES(t)和容量约束等确定的允许决策集合;υn(PCAES(t),u(t))为等效负荷平方和。f(P(t‑1))为第t‑1阶段式(14)的目标函数。f(P(t))为式(14)第t阶段的目标函数。
[0133]
[0134] 式(21)表示功率不足量的指标。
[0135] 步骤三:将所提方法在IEEE‑RBTS测试系统进行对系统鲁棒性的最优日前调度方案进行了仿真验证,并对只考虑风力发电机的不确定性和考虑风机及负荷不确定性两种场景进行仿真。
[0136] RBTS的总负载为111.5MW。每个用户类型的小时平均负载被视为采用时间系统小时负荷曲线作为数值试验的预测负荷。每台风力涡轮机的额定发电功率为1.5MW。风力涡轮机的总数是100台。每小时的平均功率作为预测的风机功率。额定膨胀功率100MW,额定压缩功率60MW,压缩空气总量640万千克,储气能力8小时满排能力。最小膨胀、压缩和SOC水平分别为30MW、10MW和33%。假定一天开始时SOC的初始水平为50%。为了采用本发明专利提出的鲁棒动态规划方法,膨胀/压缩功率在最大和最小膨胀/压缩功率之间离散。
[0137] 场景一:只考虑风力发电机的不确定性
[0138] 表1风力发电机预测误差为40%时的ΥENS
[0139]
[0140] 表2风力发电机预测误差为20%时的ΥENS
[0141]
[0142] 表3风力发电机预测误差为10%时的ΥENS
[0143]
[0144] 通过表1、2、3可以看出ROD变化的权重将影响优化结果的可靠性。权重越接近预测误差,结果越可靠。
[0145] 场景2:同时考虑风机及负荷不确定性
[0146] 表4风力发电机权重为1时的ΥENS
[0147]
[0148] 当同时风机及考虑负荷预测的不确定性时,表4中列出了具有不同κP(t)和不确定性集的ΥENS。风力发电和负荷的预测误差分别为30%、5%。得出与上述表1‑3相同的结论。

附图说明

[0076] 图1为典型的两级压缩两级膨胀绝热CAES电站结构图;
[0077] 图2为动态规划CAES控制优化路径示意图;
[0078] 图3为鲁棒DP流程图;
专利联系人(活跃度排行)
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号