CN110061528B - 一种气电联合系统日前鲁棒调度方法 - Google Patents
一种气电联合系统日前鲁棒调度方法 Download PDFInfo
- Publication number
- CN110061528B CN110061528B CN201910287669.6A CN201910287669A CN110061528B CN 110061528 B CN110061528 B CN 110061528B CN 201910287669 A CN201910287669 A CN 201910287669A CN 110061528 B CN110061528 B CN 110061528B
- Authority
- CN
- China
- Prior art keywords
- gas
- unit
- thermal power
- generating unit
- constraint
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 62
- 239000003345 natural gas Substances 0.000 claims abstract description 33
- 238000006243 chemical reaction Methods 0.000 claims abstract description 28
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 12
- 238000004519 manufacturing process Methods 0.000 claims abstract description 11
- 238000010977 unit operation Methods 0.000 claims abstract description 11
- 238000012795 verification Methods 0.000 claims abstract description 6
- 239000007789 gas Substances 0.000 claims description 142
- 239000003245 coal Substances 0.000 claims description 18
- 238000005457 optimization Methods 0.000 claims description 13
- 230000009194 climbing Effects 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000012886 linear function Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 2
- 230000002146 bilateral effect Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 230000009977 dual effect Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 102000003712 Complement factor B Human genes 0.000 description 2
- 108090000056 Complement factor B Proteins 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Wind Motors (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种气电联合系统日前鲁棒调度方法。该方法收集各个元件的技术参数,找到风电出力在其波动范围内的都能满足系统运行安全性和经济性的变量决策方案。本发明通过日前制定火电机组开机方案,并进行可行性校验以保证该机组开机方案能够应对风电不确定性;采用C&CG算法求解该气电联合系统日前鲁棒调度问题:求解决策日前开机方式,校验所得到的开机方式,若不满足风电场景则返回风电场景生产新的相关约束,迭代得到所有风电场景下均能满足机组运行约束的开机方式。本发明考虑高比例风电接入后由风电的不确定性和波动性对日前调度产生的影响,考虑天然气的动态过程和气电联合系统双边能量转换关系,提高了系统运行灵活性、安全性和经济性。
Description
技术领域
本发明涉及电气工程领域,更具体地,涉及一种气电联合系统日前鲁棒调度方法。
背景技术
随着煤炭、石油等化石能源的日益减少,风能、光能等可再生能源受到越来越多的关注。世界各国都在大力推进风电的开发,但是随着系统中接入风电比例越来越大,风电出力的预测误差和波动性无疑对日前机组组合带来了巨大的挑战,具体体现在风电出力的预测误差较大,而在当前电力系统的背景下,具有出力水平和起停状态调整缓慢特征的火电机组仍在电网电源中占有显著比重,当前火电机组的出力状态会很大限度地制约未来长时间尺度的出力调整范围。如何有效地克服风电场出力的不确定性,保证系统安全运行的同时,保证系统高效经济运行成为调度人员亟需解决的难题。
发明内容
针对现有技术的缺陷,本发明的目的在于解决由于风电出力波动性和随机性影响导致的日前机组组合不准确,系统不能及时调整的技术问题。
为实现上述目的,本发明提供一种气电联合系统日前鲁棒调度方法,包括以下步骤:
确定气电联合系统中各元件的参数,所述气电联合系统的元件包括:火电机组、燃气机组、电转气机组、风电机组、交流电网以及天然气网络;
确定气电联合系统日前鲁棒调度模型主问题的目标函数,所述主问题目标函数与火电机组运行费用、燃气机组运行费用、电转气机组运行费用以及电转气机组收益有关;
确定所述调度模型主问题的约束条件,所述约束条件包括:电功率平衡约束、燃气机组运行约束、电转气机组运行约束以及天然气系统状态约束,所述电功率平衡约束考虑火电机组的功率、燃气机组的功率、风电机组的功率、电转气机组功率以及负荷;
基于所述主问题的约束条件求解调度模型的主问题,得到火电机组开机方式,对求解得到的开机方式基于所有风电场景出力进行可行性子问题校验,确定所述调度模型子问题的目标函数,所述调度模型子问题需满足添加火电机组松弛变量的火电运行约束;
根据可行性子问题的目标函数判断由主问题求解得到的火电机组开机方式下,所有风电场景是否均可满足气电联合系统约束,若不满足,则存在风电场景不满足系统约束,重新求解主问题,得到排除前一种开机方式下的新的火电机组开机方式,重新进行可行性子问题校验,直到得到满足所有风电场景的火电机组开机方式。
可选地,所述气电联合系统日前鲁棒调度模型的目标函数为气电联合系统运行费用,包括火电机组运行费用、燃气机组运行费用及电转气机组运行费用,并扣除电转气机组收益;采用分段线性函数对火电机组的费用进行线性化,得以下目标函数:
式中,Ng为火电机组台数,Ngas为燃气机组台数,Np2g为电转气机组台数;为火电机组i在t时刻的煤耗费用,ccoal为标准煤当月实时价格,Csui,t为火电机组i在t时刻的开机费用,Csdi,t为火电机组i在t时刻的停机费用;为燃气机组i在t时刻的耗气量,为燃气机组i在t时刻对应的气价,为燃气机组i在t时刻的发电功率,为燃气机组i在t时刻的二次项成本系数;为电转气机组i在t时刻的产气量,为电转气机组i在t时刻对应的气价,为电转气机组i在t时刻的用电功率,为电转气机组i在t时刻的二次项成本系数;T表示系统的优化时段。
可选地,所述电功率平衡约束为:
可选地,所述燃气机组运行约束为:
可选地,所述电转气机组运行约束为:
可选地,所述天然气系统状态约束为:
式中MDk,t为t时刻节点k的天然气负荷质量流量,包括燃气机组电站耗气与电转气机组电站产气,MEik,t为t时刻末端为节点k的管道的天然气质量流量,MFkj,t为t时刻始端为节点k的管道的天然气质量流量;
源节点气压约束为:pi,t=psi,0,i∈源节点;
式中,pi,t为t时刻节点i的气压,psi,0为源节点恒定气压。
可选地,所述约束条件还包括:电力系统支路容量约束、火电机组运行约束、最小启停时间约束、火电煤耗松弛约束以及火电启停费用松弛约束。
可选地,可行性子问题校验需满足如下公式:
s.t.Gx+Mw+Eu*≤h
式中,x为可行性子问题优化变量,包括火电机组出力、燃气机组出力、电转气机组出力以及辅助变量s,w为随机变量,在此模型中代表风电机组随机出力,w为随机变量的下限,为随机变量的上限,b为目标函数系数列向量,对应辅助变量s的系数为1,其余均为0,h为常数,G、M、E分别表示x、w、u*的系数。
可选地,所述火电运行约束为:
式中,为松弛变量,为对应随机风电场景下的火电机组i在t时刻的出力,为火电机组i在t时刻第k种松弛变量,为火电机组i的最大上爬坡功,为火电机组i的最大下爬坡功率,Pi g,max为火电机组出力上限,Pi g,min为火电机组出力下限,为火电机组i的最大开机功率,为火电机组i的最大停机功率,为火电机组i在t时刻的最优开机方式,0代表停机,1代表开机;
若所有松弛变量s为0,说明在所得开机方式下,所有风电场景均可满足气电联合系统约束;若存在松弛变量s不为0,说明存在风电场景不满足气电联合系统约束。
可选地,采用列C&CG算法求解气电联合系统日前鲁棒调度问题,具体包括如下步骤:
步骤1,在预测场景下,求解调度模型目标函数,得到火电机组的初始开机方式u*0,并令r=0,r表示算法迭代次数;
步骤3,如果则火电机组开机方式下,所有风电场景均满足系统运行约束,终止循坏得到系统调度结果;否则,设置r=r+1,且得到可行性子问题中不满足系统约束的风电场景wr,生成新的风电场景出力约束:Gp+Mwr+Eu≤h,p表示下一次迭代计算的可行性子问题优化变量,u表示下一次迭代计算的火电机组开机方式;
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
本发明充分考虑了高比例风电接入后由于风电的不确定性和波动性对日前调度产生的影响,同时加入燃气机组和P2G(电转气)电站,考虑天然气的动态过程和气电联合系统双边能量转换关系,大大提高了系统运行的灵活性,有利于解决风电出力随机性和波动性带来的系统运行安全性隐患,有助于减小系统运行成本,提高电力系统运行的安全性和经济性。
采用本发明提出的气电联合系统日前鲁棒调度方法对气电联系系统进行调度,当风电实际出力与预测出力相差较大时,我们的系统依然能够正常运行,大大提高了系统运行的安全性和经济性。
本发明采用C&CG算法调度,大大加快计算速度,减少计算复杂性。因为风电的实际值是和预测值会有出入的,如果直接计算,将所有风电场景考虑,约束会非常多,求解非常困难。因此本发明先将主问题的风电约束限定为只考虑我们预测场景的风电,这样风电场景就一个,得到一个日前的开机方式。然后代入子问题验证是否可行,子问题中是所有的风电场景。若子问题验证不成功,说明至少有一条风电场景不满足,将这个风电场景的约束生成,添加到主问题约束中,重新求解主问题,这样求解出来的开机方式就排除了上一次不满足要求的开机方式。依次循环,直到子问题验证成功,这样就说明主问题求出来的开机方式满足所有风电场景。
附图说明
图1本发明提供的气电联合系统日前鲁棒调度方法流程图;
图2本发明提供的具有39节点的电力系统;
图3本发明提供的具有27节点的天然气系统;
图4本发明提供的基于C&CG算法调度的流程图;
图5本发明提供的基于气电联合系统日前鲁棒调度方法调度的风电鲁棒区间示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
需要说明的是,天然气作为一种清洁能源,可减少煤和石油的用量,并且天然气燃烧后产生CO2和H2O,不产生污染环境的SO2,可以改善环境,更重要的是,燃气机组具有输出范围广、调整迅速的特点,若在系统中加入燃气机组和P2G(电转气)电站,利用燃气机组调整迅速的特点,可以大大提高系统运行的灵活性,有助于改善风电出力的不确定性和波动性对系统的冲击对日前调度进行合理的决策产生重要而深远的影响。
图1为本发明实施例提供的一种气电联合系统日前鲁棒调度方法流程图,本实例以39节点电力系统与27节点天然气系统组成的气电联合系统进行分析,如图2、图3所示。该系统共有6台燃煤机组,2台燃气机组,2个风电机组,2台P2G机组。电力系统节点32、33通过燃气机组与天然气系统节点13、22相连,电力系统节点17、22通过P2G机组与天然气系统节点27、8相连。
如图1所示,本实例的主要步骤包括如下:
步骤1:收集所研模型系统中各元件经济技术参数
所述模型系统各元件包含交流电网、联络线、火电机组、燃气机组、电转气机组以及风电机组。
各元件参数含:
1)交流电网节点个数Nb,节点有功负荷Pd;
2)交流电网线路条数Nl,线路首末端节点编号,线路电抗标幺值x;交流电网节点m与节点n之间传输功率上限Cmn,系统基准容量Sb;
3)火电机组所在节点编号,火电机组出力上下限Pg,max和Pg,max,最大上下爬坡率RU和RD,火电机组最小开停机时间TU和TD,火电运行成本Fg,火电机组开停机费用系数,(见表1),标准煤实时价格ccoal(此处设为500);
表1火电机组参数
4)风电机组所在节点编号,风电机组出力Pw;
5)燃气机组所在节点编号,燃气机组的出力上限Pgas,max,燃气机组的气-电转化效率ηgas,(见表2),燃气机组的气价cgas(此处为500),燃气机组的二次项成本系数Bgas(此处为1);
6)电转气机组所在节点编号,电转气机组的功率上限Pp2g-max,电转气机组的转化效率ηp2g,(见表2),电转气机组对应的气价cp2g(此处为500),电转气机组的二次项成本系数Bp2g(此处为1);
表2燃气机组与p2g机组参数
节点 | 出力下限/MW | 出力上限/MW | 效率 |
32 | 0 | 200 | 1.8MW/(kg/s) |
33 | 0 | 200 | 1.6MW/(kg/s) |
17 | 0 | 120 | 0.10kg/s/MW |
22 | 0 | 120 | 0.15kg/s/MW |
7)天然气网络参数。
步骤2:确定气电联合系统日前鲁棒调度模型的目标函数
所提气电联合系统日前鲁棒调度模型的目标函数为气电联合系统运行费用,包括传统火电机组运行费用、燃气机组运行费用及电转气机组运行费用,并扣除电转气收益,以上费用均为在期望风电场景下的机组运行费用。采用分段线性函数对火电机组的费用进行线性化,得以下目标函数:
式中,Ng为传统火电机组台数,Ngas为燃气机组台数,Np2g为电转气机组台数;为火电机组i在t时刻的煤耗费用,Csui,t为火电机组i在t时刻的开机费用,Csdi,t为火电机组i在t时刻的停机费用;为燃气机组i在t时刻的发电功率,为燃气机组i在t时刻的耗气量,为燃气机组i在t时刻对应的气价,为燃气机组i在t时刻的二次项成本系数;为电转气机组i在t时刻的用电功率,为电转气机组i在t时刻的产气量,为电转气机组i在t时刻对应的气价,为电转气机组i在t时刻的二次项成本系数,ccoal为标准煤当月实时价格。
步骤3:确定气电联合系统日前鲁棒调度方法主问题的约束条件
1)电功率平衡约束:
2)电力系统支路容量约束:
式中Cmn为节点m与节点n之间传输功率上限,为电力系统火电机组节点注入功率与支路潮流功率关系矩阵,为电力系统燃气机组节点注入功率与支路潮流功率关系矩阵,为电力系统电转气机组节点注入功率与支路潮流功率关系矩阵,为电力系统风电场节点注入功率与支路潮流功率关系矩阵,为电力系统负荷节点注入功率与支路潮流功率关系矩阵。
3)火电机组运行约束:
式中ui,t为火电机组i在t时刻的开停机状态,为0-1变量,其中0代表停机,1代表开机,为火电机组i的最大上爬坡功率,为火电机组i的最大下爬坡功率,为火电机组i的最大开机功率,为火电机组i的最大停机功率。式(4)为火电机组i出力上下限约束,式(5)为机组i上爬坡约束,式(6)为机组i下爬坡约束,式(7)为火电机组i开机时刻最大功率约束,式(8)为火电机组i停机时刻最大功率约束。M为较大常数,为保证涉及机组停机时刻时机组爬坡约束恒成立,以及火电机组开机状态下机组最大启停功率约束恒成立。
4)最小启停时间约束:
式中Ti U为机组i需要维持的最小开机时间,Ti D为机组i需要维持的最小停机时间。仿真时段初期的最小启停时间约束需要结合机组初始开停机状态进行约束,若机组在初始时刻t=0为开机状态,则最小启停约束如下:
ui,t≤ui,t-1 t≤Ti D (12)
式中Ti U0为机组i在仿真时段之前已经启动的时长,式(12)为退化后的最小停机时间约束,保证在前Ti D内机组不能重新开机。
若机组在初始时刻t=0为停机状态,则最小启停约束如下:
ui,t-1≤ui,t t≤Ti D (14)
式中Ti D0为机组i在仿真时段之前已经停机的时长,式(14)为退化后的最小开机时间约束,保证在前Ti U内机组不能重新停机。
5)燃气机组运行约束:
6)电转气(P2G)机组运行约束:
式中Pi p2g,max为电转气机组i的功率上限,ηi p2g为电转气机组i的转化效率,为电转气机组i在t时刻的产气量。式(17)为电转气机组i功率上下限约束,式(18)为电转气机组i的耗电-产气对应关系约束。
7)火电煤耗松弛约束:
式中ai,k、bi,k为火电机组i煤耗曲线分段线性化后第k段一次项系数与纵轴截距系数,kmax为火电机组煤耗曲线分段数。
8)火电启停费用松弛约束:
9)天然气系统状态约束:
基于流体动力学模型,在一定假设情况下,我们简化得到天然气系统的物质平衡方程与动量方程,其可分别由式(22)和式(23)表达:
式中:c为声速;MEij,t/MFij,t为t时刻管道ij的流入/流出天然气质量流量;pi,t为t时刻节点i的气压;Lij为管道ij的长度;Aij为管道ij横截面积;dij为管道ij的直径;ωij为表示天然气在管道ij的流速;t为时间步长,本申请取t=1h,λ表示摩擦系数。
天然气节点处流出与流入的天然气质量流量守恒,如下式所示:
式中MDk,t为t时刻节点k的天然气负荷质量流量,包括燃气电站耗气与P2G电站产气。
源节点气压约束,如下式所示:
pi,t=psi,0,i∈源节点 (25)
式中psi,0为源节点恒定气压。
节点气压上下限余数:
步骤4:确定气电联合系统日前鲁棒调度方法的子问题
为了保证所有风电场景下均满足气电联合系统运行约束,需要对由主问题求解得到的机组组合(开机方式)进行可行性校验。可行性子问题可以简写成如下形式:
式中x为可行性子问题优化变量,包括火电出力、燃气机组出力、P2G机组出力以及辅助变量s,w为随机变量,即风电随机出力,b为目标函数系数列向量,对应辅助变量s的系数为1,其余均为0。s.t.代表约束,G,M,E分别代表x,w,u的系数,h代表“常数”,他们泛指约束。例如在电功率平衡约束中,x代表火电出力、燃机出力等等,w代表风电出力,u代表机组开机方式,h代表负荷;在电功率平衡约束下,系数G全为1,系数M全为1,系数E全为0。
若可行性子问题优化目标函数为0,即所有松弛变量为0,说明在所得开机方式下,所有风电场景均可满足气电联合系统约束;若可行性子问题优化目标函数不为0,即存在松弛变量不为0,说明存在风电场景不满足系统约束。
步骤5:确定求解气电联合系统日前鲁棒调度的方法
本发明采用列C&CG算法求解该气电联合系统日前鲁棒调度问题:主问题用于求解决策日前开机方式,子问题用于校验所得到的开机方式能够满足所有风电场景,若不满足则返回风电场景生产新的相关约束,最后迭代得到所有风电场景下均能满足机组运行约束的开机方式。
1)可行性子问题求解方法
式(27)为max-min问题,为内外层优化目标函数相同的双层优化问题,可先采用对偶理论将可行性子问题中的min问题转化为其对偶形式,对于min问题而来,随机变量w为已知参数(非优化变量),因此其式(27)问题可转化为如下形式:
其中π为原问题中约束Gx+Mw+Eu*≤h的对偶变量。对偶问题(30)的目标函数中出现了非线性项(Mw)Tπ,本文中采用二进制表示与大M法将该非线性项简化。
其中(MΔwv)Tπ为整数变量与连续变量相乘的双线性项,可用大M法将其线性化,令xtem为中间变量表示(MΔwv)Tπ,其满足如下约束:
式中Nc为足够大常数,当v=0时,(MΔwv)Tπ=xtem=0;当v=1时,(MΔwv)Tπ=xtem=(MΔw)Tπ。至此,原可行性子问题转化为可直接求解的混合整数优化问题:
2)C&CG算法迭代
C&CG算法通过求解可行性子问题获得不满足系统运行约束的风电场景,再将该风电场景带入主问题求解同时满足该风电场景与预测风电场景的开机方式。求解步骤可以阐述为:
1.在预测场景下,求解主问题(1),得到初始开机方式u*0,并令r=0;
Gp+Mwr+Eu≤h (35)
整体算法流程图如图4所示。
仿真结果:风电的波动区间如图5,分别在风电预测值的基础上向下波动50%与向上波动30%。通过图5我们可以看到,通过我们提出的气电联合系统日前鲁棒调度方法,当风电实际出力与预测出力相差较大时,我们的系统依然能够正常运行,我们所提出的方法完全能够应用到实际算例当中,这将会大大提高了电力系统运行的安全性和经济性。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种气电联合系统日前鲁棒调度方法,其特征在于,包括以下步骤:
确定气电联合系统中各元件的参数,所述气电联合系统的元件包括:火电机组、燃气机组、电转气机组、风电机组、交流电网以及天然气网络;
确定气电联合系统日前鲁棒调度模型主问题的目标函数,所述主问题目标函数与火电机组运行费用、燃气机组运行费用、电转气机组运行费用以及电转气机组收益有关;
确定所述调度模型主问题的约束条件,所述约束条件包括:电功率平衡约束、燃气机组运行约束、电转气机组运行约束以及天然气系统状态约束,所述电功率平衡约束考虑火电机组的功率、燃气机组的功率、风电机组的功率、电转气机组功率以及负荷;
基于所述主问题的约束条件求解调度模型的主问题,得到火电机组开机方式,对求解得到的开机方式基于所有风电场景出力进行可行性子问题校验,确定所述调度模型子问题的目标函数,所述调度模型子问题需满足添加火电机组松弛变量的火电运行约束;
根据可行性子问题的目标函数判断由主问题求解得到的火电机组开机方式下,所有风电场景是否均可满足气电联合系统约束,若不满足,则存在风电场景不满足系统约束,重新求解主问题,得到排除前一种开机方式下的新的火电机组开机方式,重新进行可行性子问题校验,直到得到满足所有风电场景的火电机组开机方式;
所述气电联合系统日前鲁棒调度模型的目标函数为气电联合系统运行费用,包括火电机组运行费用、燃气机组运行费用及电转气机组运行费用,并扣除电转气机组收益;采用分段线性函数对火电机组的费用进行线性化,得以下目标函数:
式中,Ng为火电机组台数,Ngas为燃气机组台数,Np2g为电转气机组台数;为火电机组i在t时刻的煤耗费用,ccoal为标准煤当月实时价格,Csui,t为火电机组i在t时刻的开机费用,Csdi,t为火电机组i在t时刻的停机费用;为燃气机组i在t时刻的耗气量,为燃气机组i在t时刻对应的气价,为燃气机组i在t时刻的发电功率,为燃气机组i在t时刻的二次项成本系数;为电转气机组i在t时刻的产气量,为电转气机组i在t时刻对应的气价,为电转气机组i在t时刻的用电功率,为电转气机组i在t时刻的二次项成本系数;T表示系统的优化时段;
所述天然气系统状态约束中,天然气系统的物质平衡方程与动量方程,分别由式(22)和式(23)表达:
式中:c为声速;MEij,t/MFij,t为t时刻管道ij的流入/流出天然气质量流量;pi,t为t时刻节点i的气压;Lij为管道ij的长度;Aij为管道ij横截面积;dij为管道ij的直径;ωij为表示天然气在管道ij的流速;t为时间步长,λ表示摩擦系数,Δt表示优化的时段。
6.根据权利要求1至5任一项所述的气电联合系统日前鲁棒调度方法,其特征在于,所述约束条件还包括:电力系统支路容量约束、火电机组运行约束、最小启停时间约束、火电煤耗松弛约束以及火电启停费用松弛约束。
9.根据权利要求8所述的气电联合系统日前鲁棒调度方法,其特征在于,采用列C&CG算法求解气电联合系统日前鲁棒调度问题,具体包括如下步骤:
步骤1,在预测场景下,求解调度模型目标函数,得到火电机组的初始开机方式u*0,并令r=0,r表示算法迭代次数;
步骤3,如果则火电机组开机方式u*r下,所有风电场景均满足系统运行约束,终止循坏得到系统调度结果;否则,设置r=r+1,且得到可行性子问题中不满足系统约束的风电场景wr,生成新的风电场景出力约束:Gp+Mwr+Eu≤h,p表示下一次迭代计算的可行性子问题优化变量,u表示下一次迭代计算的火电机组开机方式;
步骤4,将新的风电场景出力约束代入主问题求解,得到新的火电开机方式u*r,并返回步骤2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910287669.6A CN110061528B (zh) | 2019-04-11 | 2019-04-11 | 一种气电联合系统日前鲁棒调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910287669.6A CN110061528B (zh) | 2019-04-11 | 2019-04-11 | 一种气电联合系统日前鲁棒调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110061528A CN110061528A (zh) | 2019-07-26 |
CN110061528B true CN110061528B (zh) | 2021-03-26 |
Family
ID=67318632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910287669.6A Expired - Fee Related CN110061528B (zh) | 2019-04-11 | 2019-04-11 | 一种气电联合系统日前鲁棒调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110061528B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110854929B (zh) * | 2019-11-12 | 2021-05-18 | 华中科技大学 | 一种计及时段内不确定性的日前调度方法 |
CN111401647B (zh) * | 2020-03-23 | 2022-04-08 | 清华大学 | 考虑不确定性传递的电气耦合系统分布式优化调度方法 |
CN111769553A (zh) * | 2020-07-06 | 2020-10-13 | 四川大学 | 考虑风电不确定性的气电互联系统优化运行方法 |
CN112018756A (zh) * | 2020-08-04 | 2020-12-01 | 四川大学 | 气电联合系统日前鲁棒协调优化调度方法 |
CN112508236A (zh) * | 2020-11-20 | 2021-03-16 | 华北电力大学 | 一种基于电转气提供灵活调节服务的日前调度方法及系统 |
CN112564182A (zh) * | 2020-11-25 | 2021-03-26 | 云南电网有限责任公司 | 一种基于分布鲁棒机会约束的多能联合系统统一调度方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9785130B2 (en) * | 2014-04-10 | 2017-10-10 | Nec Corporation | Decentralized energy management platform |
US9926852B2 (en) * | 2015-03-03 | 2018-03-27 | General Electric Company | Methods and systems for enhancing control of power plant generating units |
CN107977744B (zh) * | 2017-11-30 | 2021-07-13 | 三峡大学 | 一种基于传统Benders分解法的电力系统日前鲁棒调度方法 |
-
2019
- 2019-04-11 CN CN201910287669.6A patent/CN110061528B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110061528A (zh) | 2019-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110061528B (zh) | 一种气电联合系统日前鲁棒调度方法 | |
CN110571789B (zh) | 数据驱动下基于风电不确定性的电热气网三阶段调度方法 | |
CN107944757A (zh) | 电力互联系统可再生能源消纳能力评估分析方法 | |
CN109031952B (zh) | 一种电-气互联综合能源系统混合控制方法 | |
CN108123492A (zh) | 一种考虑火电电蓄热联合调峰的日前发电计划优化方法 | |
CN109325621B (zh) | 一种园区能源互联网两阶段优化调度控制方法 | |
CN103683326A (zh) | 一种区域电网风电多点接入最佳接纳能力的计算方法 | |
CN109787227B (zh) | 一种多能互补系统多时间尺度优化调度方法 | |
CN112701687B (zh) | 考虑价格型联合需求响应的气电配网系统鲁棒优化运行方法 | |
CN109818347B (zh) | 一种电力系统风电消纳能力的评估方法 | |
CN104299173B (zh) | 一种适用于多种能源接入的鲁棒优化日前调度方法 | |
CN112053035B (zh) | 考虑经济性与灵活性的输电通道与储能联合规划方法 | |
CN104795828A (zh) | 基于遗传算法的风储容量配置方法 | |
Zhang et al. | Mid-long term optimal dispatching method of power system with large-scale wind-photovoltaic-hydro power generation | |
CN104935017A (zh) | 基于改进轻鲁棒优化模型的风电与火电机组组合方法 | |
CN105244870A (zh) | 一种快速计算电网风电场弃风率及机组发电量的方法 | |
CN114266382A (zh) | 一种考虑热惯性的热电联产系统两阶段优化调度方法 | |
CN106684928A (zh) | 一种基于调峰成本的电网调峰裕度计算方法 | |
CN107769266A (zh) | 一种多时间尺度发电与备用联合优化方法 | |
CN113837589B (zh) | 一种综合能源系统n-1调度方法及装置 | |
CN109510238B (zh) | 一种高效求解水电火电风电的协调调度机组组合方法 | |
CN106655153A (zh) | 二类状态变量的电力系统机组组合模型及其建立方法 | |
Su et al. | An energy efficiency index formation and analysis of integrated energy system based on exergy efficiency | |
CN105515062B (zh) | 基于风电和水电的电力系统调峰方法 | |
CN110429591B (zh) | 一种基于电力系统时序耦合性的输电网利用率评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210326 |