CN113256045A - 考虑风光不确定性的园区综合能源***日前经济调度方法 - Google Patents
考虑风光不确定性的园区综合能源***日前经济调度方法 Download PDFInfo
- Publication number
- CN113256045A CN113256045A CN202010773898.1A CN202010773898A CN113256045A CN 113256045 A CN113256045 A CN 113256045A CN 202010773898 A CN202010773898 A CN 202010773898A CN 113256045 A CN113256045 A CN 113256045A
- Authority
- CN
- China
- Prior art keywords
- power
- gas
- energy
- wind
- maximum
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000004044 response Effects 0.000 claims abstract description 53
- 238000005457 optimization Methods 0.000 claims abstract description 37
- 239000007789 gas Substances 0.000 claims description 150
- 238000003860 storage Methods 0.000 claims description 46
- 238000005338 heat storage Methods 0.000 claims description 40
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 34
- 238000004146 energy storage Methods 0.000 claims description 24
- 238000006243 chemical reaction Methods 0.000 claims description 23
- 238000010248 power generation Methods 0.000 claims description 22
- 230000009194 climbing Effects 0.000 claims description 19
- 238000010438 heat treatment Methods 0.000 claims description 18
- 239000003345 natural gas Substances 0.000 claims description 17
- 230000005611 electricity Effects 0.000 claims description 15
- 230000033228 biological regulation Effects 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 8
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 230000009977 dual effect Effects 0.000 claims description 5
- WKBOTKDWSSQWDR-UHFFFAOYSA-N Bromine atom Chemical compound [Br] WKBOTKDWSSQWDR-UHFFFAOYSA-N 0.000 claims description 4
- GDTBXPJZTBHREO-UHFFFAOYSA-N bromine Substances BrBr GDTBXPJZTBHREO-UHFFFAOYSA-N 0.000 claims description 4
- 229910052794 bromium Inorganic materials 0.000 claims description 4
- 238000002485 combustion reaction Methods 0.000 claims description 4
- 239000003546 flue gas Substances 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 4
- 238000012546 transfer Methods 0.000 claims description 4
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 3
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 claims description 3
- 238000005265 energy consumption Methods 0.000 claims description 3
- 230000017525 heat dissipation Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 238000009826 distribution Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 abstract 1
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 230000001965 increasing effect Effects 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000002356 single layer Substances 0.000 description 2
- 241000711969 Chandipura virus Species 0.000 description 1
- 208000015951 Cytophagic histiocytic panniculitis Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000002355 dual-layer Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
- G06Q10/06375—Prediction of business process outcome or impact based on a proposed change
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑风光不确定性的园区综合能源***日前经济调度方法,首先以基础场景运行成本最小为目标函数并结合园区***相关设备及运行约束建立了考虑联合热电需求响应的两阶段可调鲁棒优化日前调度模型。由于CCG法在处理鲁棒优化模型上具有更好的计算效率,本发明应用CCG法进行求解。最后应用商业求解器Gurobi进行算例仿真,得到考虑联合热电需求响应的两阶段鲁棒优化调度策略,验证了该方法可以较好的处理***中出现的不确定性,同时可以降低园区运行成本、促进新能源的消纳。
Description
技术领域
本发明属于综合能源***优化运行技术领域,特别涉及一种考虑风光不 确定性的园区综合能源***日前经济调度方法。
背景技术
风电实际出力与预测出力方差的大小远大于传统的负荷方差,传统的电 力***运行方式不足以保证***的可靠性。由于电力***的物理限制(如常 规发电机和输电线路容量的爬坡限制),风电功率削减频繁发生,导致风电 利用率低,从长远来看抑制了风电投资的积极性。当大规模可再生能源接入 电力***,可再生能源的强随机性、间歇性使得其出力难以被精确预测,在 电力***规划调度中,可再生能源出力预测值的误差将造成不可忽略的影响, 换言之,可再生能源发电的随机波动降低了能源***的供应灵活性。因此,传统确定性优化调度方法不能满足要求,须在原有确定性优化调度方法的基 础上考虑可再生能源的不确定性。
随机优化和鲁棒优化是解决含不确定性优化问题的典型方法。与随机优 化相比,鲁棒优化具有获取数据简单、求解速度快,适用于求解大型不确定 性问题的求解等优点,在处理不确定性问题中得到了广泛的应用。目前,已 有许多学者对鲁棒优化在综合能源***中的应用做了研究,但现有文献在利 用鲁棒优化处理综合能源***风光不确定性中未充分考虑***中各类设备的 能源耦合关系和负荷侧需求侧响应对可再生能源不确定性的影响。
较为全面的综合能源***模型可通过各设备间的能量转换关系,自适应 地调整机组出力以适应可再生能源发电量的变化,保证***运行的安全性。 联合热电需求响应可以充分利用各种设备之间的耦合关系,进一步提高*** 应对可再生能源不确定性的能力,减少弃风、弃光,加大可再生能源的渗透 性。因此,在现有研究基础上,进一步考虑联合热电需求响应的园区综合能 源***两阶段可调鲁棒优化模型对综合能源***的优化运行具有重要意义。
发明内容
本发明所要解决的技术问题是提出一种考虑风光不确定性的园区综合 能源***日前经济调度方法,建立一种考虑联合热电需求响应的两阶段可调 鲁棒优化模型来处理园区内风光出力造成的不确定性问题,通过对偶理论将 模型中的双层max-min问题转换为单层max问题,并采用CCG(column and constraint generation)法进行求解,提高求解速度的同时使得调度结果更 加符合实际,进一步促进了新能源的消纳。
为解决上述技术问题,本发明采用的技术方案是:为考虑可再生能源不 确定性,针对园区综合能源***的日前经济调度问题,提出一种考虑风光不 确定性的园区综合能源***日前经济调度方法,建立一种考虑联合热电需求 响应的两阶段可调鲁棒模型。园区综合能源***日前经济调度是针对以风电、 光伏出力值为预测值的基础场景,当运行中出现不确定性时,它可以自适应、 安全地重新分配发电机组、供热机组、P2G设备、储能设备及园区与上级网 络的能量交换。这种优化调度决策完全符合两阶段可调鲁棒模型的思想。换言之,第一阶段在基础场景下决定燃气轮机和热电联产机组的机组启停状态, 第二阶段在不确定性出现时搜索造成***不安全最恶劣的场景。本发明所提 出的两阶段可调鲁棒调度模型假设机组启停状态是第一阶段变量,即当*** 不确定变量在其波动区间内变化时,机组启停状态将保持不变,这是因为大 多数发电机组的物理特性限制了它们在不确定情况下无法快速改变机组启停 状态。
一种考虑风光不确定性的园区综合能源***日前经济调度方法,包括以 下步骤:
(1)确定多能源园区的具体组成,包含引入的新能源形式和具体的设备组成;
(2)建立园区内各能源转换设备模型;
(3)建立需求响应模型;
(4)在满足***安全约束的前提下,以基础场景的运行成本最小为目标函数, 建立考虑联合热电需求响应的两阶段可调鲁棒模型;
(5)得到园区综合能源***日前经济调度的两阶段可调鲁棒优化模型的抽象 表达式;
(6)建立最恶劣场景识别的最大最小子问题;
(7)利用CCG法求解考虑联合热电需求响应的两阶段可调鲁棒模型;
(8)输入园区综合能源***能源接入、新能源出力数据、各设备参数、运行 参数,采用商业求解器Gurobi对考虑风光不确定性的园区综合能源***日前经济 调度两阶段鲁棒优化模型进行求解,得到其调度策略。
进一步的,在步骤(1)所述园区综合能源***的具体组成具体如下:
(1.1)接入园区综合能源***的新能源形式为:风电和光伏发电;
(1.2)引入园区综合能源***的能源转换设备有:电转气设备、电锅炉、燃 气轮机、热电联产机组、储气/储热设备。
步骤(2)所述园区内各能源转换设备模型如下;
(2.1)电转气设备模型
(2.2)电锅炉模型
式中:t为调度时间;n为电锅炉索引;和分别为第n台电锅炉在t时 段的耗电功率和产热功率;为第n台电锅炉的电热转换效率,分 为第n台电锅炉最小、最大制热功率;为第n台电锅炉在t时段的启停状态(1 代表开机,0代表停机)。
(2.3)燃气轮机模型
式中:t为调度时间;q为燃气轮机索引;和分别表示燃气轮机的发电 功率和耗气功率;F(·)表示燃气轮机能耗曲线;和分别表示燃气轮机的 开机和关机所需天然气消耗量;aq、bq和cq表示F(·)的燃气系数;LHANG为天然气低 热值,取9.7kWh/m3;分为第q台燃气轮机最小、最大发电功率。表示第q台燃气轮机在t时段的发电功率;为第q台燃气轮机在t时段的启停状 态(1代表开机,0代表停机),为第q台燃气轮机在t时段的启停状态; 为第q台燃气轮机的上爬坡率、下爬坡率,为第q台燃气轮机在t-1 时段内连续开机、停机时间,为第q台燃气轮机在时段t内的最小开机、 停机时间。
(2.4)热电联产机组模型
式中:t为调度时间;p为热电联产机组索引;和分别表示热电联产 机组(CHP)的产热功率和耗气功率;为微燃机t时段的发电功率、发电 效率,为散热损失率,和分别为溴冷机的制热系数和烟气回收率;LHANG为 天然气低热值,取9.7kWh/m3;分为第p台CHP机组最小、最大发 电功率;表示第p台热电联产机组在t时段的发电功率;为第p台CHP机组 在t时段的启停状态(1代表开机,0代表停机),为第P台热电联产机组在t 时段的启停状态;为第p台CHP机组的上爬坡率、下爬坡率; 为第p台CHP机组在t-1时段内连续开机、停机时间,为第p台 CHP机组在时段t内的最小开机、停机时间。
(2.5)储能设备模型
ESS(t)=(1-σS)·ESS(t-1)+ESin(t)·ηin-ESout(t)/ηout
式中:t为调度时间;ESS(t)、ESin(t)、ESout(t)分别为t时段储能设备的储能量、 存储能量功率和释放能量功率;ESS(t-1)表示t-1时段储能设备的储能量;σS为储 能***的自耗率;ηin、ηout分别为储能设备储能、放能效率;为t时段储 气、放气功率,GSin,max、GSout,max分别为储气设备的最大储气、放气功率;为t时 段储气设备的储气量,为t-1时段储气设备的储气量;CGS,min、CGS,max分别为储 气设备的最小、最大储气容量,ηCGS、ηGS,in、ηGS,out为储气设备自耗率、储气效率、 放气效率。为t时段储热、放热功率,HSin,max、HSout,max分别为储热设备 的最大储热、放热功率;为t时段储热设备的储气量,为t-1时段储热设备 的储气量;CHS,min、CHS,max分别为储热设备的最小、最大储热容量,ηCHS、ηHS,in、ηHS,out为储热设备自耗率、储热效率、放热效率。
步骤(3)所述需求响应模型具体如下:
Pt DR=Pt DR,inte+Pt DR,shif
式中:t为调度时间;Pt DR、Pt DR,inte、Pt DR,shif、Pt LD,max为t时段电负荷预测 值、需求侧响应电负荷、需求侧转移电负荷、***允许的最大电负荷,Pinte,max为调 度时间内最大可中断电负荷。为t时段最大可中断和可转出的电负荷 比例。Pt DR,shif为正代表转出可平移负荷,反之,代表转入可平移负荷。
式中:t为调度时间;Ht DR、和分别为时间段t热负荷 预测值、需求侧响应热负荷、需求侧可响应热负荷比例、***允许的最大热负荷 和考虑需求响应后的电负荷;HDR,max为调度时间内最大可中断热负荷;Nt为整个 调度时间段。
步骤(4)所述考虑联合热电需求响应的两阶段可调鲁棒模型具体如下:
(4.1)目标函数
式中:t为调度时间;分别为购电、气费用;分别为 第p台CHP的开机、关机费用;为弃风惩罚费用,为弃光惩罚费用; 为售电收益,CE,cc为需求侧响应中断电负荷的补偿成本;Pt in为购电功率,为 购气功率,分别为第i台风机、第j组光伏电池在t时段的弃风、弃 光功率;Pt out为售电功率,Pt DR,inte为需求侧中断电负荷。Δt为调度时间间隔,Nt为 调度时段数。分别为单位购电、购气、售电价格,φt WT、φt PV分别 为单位弃风、弃光的惩罚价格,为单位需求侧响应中断电负荷的补偿价格。为第p台CHP在t时段的启停状态(1代表开机,0代表停机),分别 为第p台CHP的开机、停机一次的费用。NWT、NPV、NCHP分别为风机、光伏电池、 电锅炉、CHP的数量。
(4.2)约束条件
v1t,v2t,v3t,v4t≥0
Pt u,out≤Pout,min,Pt u,out≤Pout,max:(λ9_3,t,λ9_4,t)
Pt u,DR=Pt u,DR,inte+Pt u,DR,shif:(λ10_6,t)
式中:t为调度时间;(·)u为在风光出力实时变化下对应的调整后变量;v1,t和v2,t为电力平衡约束松弛变量;v3,t和v4,t为热量平衡约束松弛变量;ΩWT、ΩPV分别为风 电、光伏出力不确定合集;分别为风电、光伏出力与预测值的偏差, 为风电、光伏出力偏差值与预测值的比例;为不确定合集中的 0-1变量;Δi、Δj分别为风电、光伏出力不确定性预算值;为第p 台CHP机组上爬坡和下爬坡纠错能力;为第q台燃气轮机机组上爬坡 和下爬坡纠错能力;λ(·)为约束条件对应的对偶变量。
步骤(5)所述园区综合能源***日前经济调度的两阶段可调鲁棒优化模 型的抽象表达式具体如下:
s.t. Ax+By≤b
式中:x代表与CHP、燃气轮机相关的机组启停状态,y、z分别代表基础场 景和根据风光出力变换调整的***其余机组调度出力,u为与风电、光伏出力不确 定性相关的不确定变量,cb、cg、A、B、b、f、h、C、D、E、F、G可通 过(4)中的目标函数和约束条件得出。
步骤(6)所述最恶劣场景识别的最大最小子问题具体如下:
s.t. λ9_1,t≤0,λ9_2,t≤0,λ9_3,t≤0,λ9_4,t≤0,λ9_5,t≤0,λ9_6,t≤0,t∈NT
λ6_2,q,t≤0,λ6_3,q,t≤0,λ6_4,q,t≤0,λ6_6,q,t≤0,λ6_7,q,t≤0,q∈NGT,t∈NT
λ7_2,p,t≤0,λ7_3,p,t≤0,λ7_4,p,t≤0,λ7_6,p,t≤0,λ7_7,p,t≤0,p∈NCHP,t∈NT
λ8_1,t≤0,λ8_2,t≤0,λ8_3,t≤0,λ8_4,t≤0,t∈NT
λ13_1,t≤0,λ13_2,t≤0,λ13_3,t≤0,λ13_4,t≤0,t∈NT
λ4_2,m,t≤0,λ4_3,m,t≤0,t∈NT,m∈NP2G,t∈NT
λ5_2,n,t≤0,λ5_3,n,t≤0,t∈NT,n∈NEB,t∈NT
λ10_1,t≤0,λ10_2,t≤0,λ10_3,t≤0,λ10_4,t≤0,t∈NT
λ11_3,t≤0,λ11_3,t≤0,λ12_1,t≤0,λ12_2,t≤0,t∈NT
λ10_7≤0,λ10_8≤0,λ11_4≤0,λ11_5≤0,t∈NT
s.t. -λ1,t-λ6_1,q,t·bq-λ6_2,q,t+λ6_3,q,t-λ6_4,q,t+1+
λ6_4,q,t+λ6_5,p,t+1-λ6_5,p,t-λ6_6,q,t+λ6_7,q,t≤0
-λ1,t-λ9_1,t+λ9_2,t≤0,t∈NT
λ1,t-λ9_3,t+λ9_4,t≤0,t∈NT
λ1,t+λ10_5,t≤0,t∈NT
-λ10_1,t-λ10_2,t+λ10_3,t+λ10_5,t-λ10_7+λ10_8=0,t∈NT
-λ10_1,t+λ10_4,t+λ10_5,t+λ10_5,t+λ10_6=0,t∈NT
λ1,t+λ12_1,i,t≤0,t∈NT,i∈NWT
λ1,t+λ12_2,j,t≤0,t∈NT,j∈NPV
-λ2,t+λ7_1,p,t≤0,t∈NT,p∈NCHP
-λ2,t+λ6_1,q,t≤0,t∈NT,q∈NGT
-λ2,t+λ4_1,m,t-λ4_2,m,t+λ4_3,m,t≤0,t∈NT,m∈NP2G
λ2,t-λ9_5,t+λ9_6,t≤0,t∈NT
λ2,t+λ8_2,t+λ8_5,t/ηGS,out≤0,t∈NT
-λ2,t+λ8_1,t-λ8_5,t·ηGS,in≤0,t∈NT
-λ8_3,t+λ8_4,t-λ8_5,t+1·ηGS,in+λ8_5,t≤0,t∈NT-λ3,t·ηh+λ13_2,t+λ13_5,t/ηHS,out≤0,t∈NT
λ3,t·ηh+λ13_1,t-λ13_5,t·ηHS,in≤0,t∈NT-λ13_3,t+λ13_4,t-λ13_5,t+1·ηHS,in+λ13_5,t≤0,t∈NT
-λ3,t·ηh+λ5_1,n,t-λ5_2,n,t+λ5_3,n,t≤0,t∈NT,n∈NEB-λ3,t·ηh+λ7_8,n,t≤0,t∈NT,n∈NEB
λ3,t+λ11_1,t≤0,t∈NTλ11_1,t-λ11_2,t+λ11_3,t-λ11_4+λ11_5,t≤0,t∈NT
-1≤λ1,t≤1,-1≤λ3,t≤1,t∈NT
步骤(7)所述利用CCG法求解考虑联合热电需求响应的两阶段可调鲁棒 模型聚具体如下:
(7.1)主问题
Ax+By≤b
(7.2)子问题
(7.3)CCG求解步骤
步骤1:令迭代计数器s=0,设置***允许的违反安全规定最大值εRO。
步骤2:求解主问题,若有解,得到***机组启停状态x和机组出力安排 y,进行步骤3;反之,停止迭代并输出无解。
步骤3:根据步骤2中求解得到的x和y,求解最大最小子问题,找到导 致最大可能违反安全规定值的最恶劣场景下的风力和光伏出力大小。
步骤4:如果步骤3中求解出的最大可能违反安全规定值小于εRO,则x和 y是最终优化方案并停止迭代;反之,令s=s+1,根据步骤3中求解出来的 最恶劣场景下风电和光伏出力值向主问题中增加以下所 示的CCG约束,返回步骤2。
fTvs≤εRO
步骤(8)所述园区综合能源***数据还包括园区综合能源***的具体组 成、能源价格、各能源转换设备的设备参数及取值、需求响应比例、基础场 景和最恶劣场景下的新能源出力波动情况、最大违反安全规定值和新能源出 力与负荷的预测值。
与现有技术相比,本发明的有益效果是:
(1)采用鲁棒优化解决含不确定性的优化问题,并建立了两阶段可调鲁 棒优化模型,使得园区综合能源***调度过程中只需获取简单数据就能得到 贴近实际的结果,且求解速度快。改进的两阶段可调鲁棒优化模型在考虑综 合能源***风光不确定性的同时充分考虑***中各类设备的能源耦合关系和 负荷侧需求侧响应对可再生能源不确定性的影响。园区综合能源***日前经 济调度是针对以风电、光伏出力值为预测值的基础场景,当运行中出现不确 定性时,它可以自适应、安全地重新分配发电机组、供热机组、P2G设备、储能设备及园区与上级网络的能量交换。两阶段鲁棒优化保证了***在任何情 况下都能满足安全约束,最小化了基础场景运行成本,同时满足***安全性 和经济性要求。
(2)考虑热电联合需求响应后,充分发挥气储能、热储能以及电热负荷 需求响应的作用,联合热电需求响应可以充分利用各种设备之间的耦合关系, 进一步提高***应对可再生能源不确定性的能力,减少弃风、弃光,加大可 再生能源的渗透性。较为全面的综合能源***模型可通过各设备间的能量转 换关系,自适应地调整机组出力以适应可再生能源发电量的变化,保证*** 运行的安全性,促进可再生能源的消纳、提升***的鲁棒性,提高园区的经 济效益。
附图说明
图1是本发明所述方法的步骤流程图;
图2是园区综合能源***的具体组成图;
图3是CCG法求解流程图;
图4是园区多能源综合***冬季典型日的风机、光伏出力、电负荷和热 负荷的预测值。
具体实施方式
为了详尽说明本发明所公开的技术方案,下面结合附图和具体实施例对 本发明作进一步说明。
本发明公开的是一种考虑风光不确定性的园区综合能源***日前经济调 度方法。具体实施步骤流程如图1所示,本发明技术方案包括以下步骤:
步骤1:确定园区综合能源***的具体组成,包含引入的新能源形式和具 体的设备组成。
(1.1)接入园区综合能源***的新能源形式为:风电和光伏发电;
(1.2)引入园区综合能源***的能源转换设备有:电转气设备、电锅炉、 燃气轮机、热电联产机组、储气/储热设备。
步骤2:建立园区内各能源转换设备模型。
(2.1)电转气设备模型P2G技术可实现电能向天然气的转换,天然气通 过天然气管道进行大容量、长时间和远距离传输,为可再生能源消纳提供了 有力的技术支持,能实现风电的大范围、长距离的时空转移,同时P2G响应 迅速,具备较强的应用前景。P2G设备制气功率与耗电功率间的关系和制气功 率限制如下:
(2.2)电锅炉模型
电锅炉的引入能打破CHP机组的电热耦合硬性约束,改变传统“以热定 电”调度方式。电锅炉可协调电热负荷的峰谷,其制热量与耗电量之间的关 系和制热量约束如下:
式中:t为调度时间;n为电锅炉索引;和分别为第n台电锅炉在t 时段的耗电功率和产热功率;为第n台电锅炉的电热转换效率, 分为第n台电锅炉最小、最大制热功率;为第n台电锅炉在t时段的 启停状态(1代表开机,0代表停机)。
(2.3)燃气轮机模型
燃气轮机将天然气中化学能转换为电能,其消耗天然气功率与发电功率 间的关系、发电功率限制、爬坡约束与最小开关机时间约束如下:
式中:t为调度时间;q为燃气轮机索引;和分别表示燃气轮机的 发电功率和耗气功率;F(·)表示燃气轮机能耗曲线;和分别表示燃气 轮机的开机和关机所需天然气消耗量;aq、bq和cq表示F(·)的燃气系数;LHANG为天然气低热值,取9.7kWh/m3;分为第q台燃气轮机最小、最 大发电功率。表示第q台燃气轮机在t时段的发电功率;为第q台燃气 轮机在t时段的启停状态(1代表开机,0代表停机),为第q台燃气轮机 在t时段的启停状态;为第q台燃气轮机的上爬坡率、下爬坡率, 为第q台燃气轮机在t-1时段内连续开机、停机时间,为第q台燃气轮机在时段t内的最小开机、停机时间。
(2.4)热电联产机组模型
忽略外界环境变化对发电、燃料燃烧效率的影响,其热电关系数学、气 电关系、发电功率间的关系、发电功率限制、爬坡约束与最小开关机时间约 束如下:
式中:t为调度时间;p为热电联产机组索引;和分别表示热电 联产机组(CHP)的产热功率和耗气功率;为微燃机t时段的发电 功率、发电效率,为散热损失率,和分别为溴冷机的制热系数和烟气 回收率;LHANG为天然气低热值,取9.7kWh/m3;分为第p台CHP 机组最小、最大发电功率;表示第p台热电联产机组在t时段的发电功率; 为第p台CHP机组在t时段的启停状态(1代表开机,0代表停机),为第P台热电联产机组在t时段的启停状态;为第p台CHP机 组的上爬坡率、下爬坡率;为第p台CHP机组在t-1时段内连 续开机、停机时间,为第p台CHP机组在时段t内的最小开机、 停机时间。
(2.5)储能设备模型
储能设备模型、储气、储热设备的储/放气、储/放热功率约束和t时刻储 气、储热设备的储气、储热容量约束以及t时刻储气、储热设备的储气、储热 容量与t-1时刻储气、储热量和t时刻储/放气、储/放热功率之间的关系如下:
ESS(t)=(1-σS)·ESS(t-1)+ESin(t)·ηin-ESout(t)/ηout
式中:t为调度时间;ESS(t)、ESin(t)、ESout(t)分别为t时段储能设备的储 能量、存储能量功率和释放能量功率;ESS(t-1)表示t-1时段储能设备的储能 量;σS为储能***的自耗率;ηin、ηout分别为储能设备储能、放能效率; 为t时段储气、放气功率,GSin ,max、GSout,max分别为储气设备的最大储气、 放气功率;为t时段储气设备的储气量,为t-1时段储气设备的储气量; CGS,min、CGS,max分别为储气设备的最小、最大储气容量,ηCGS、ηGS,in、ηGS,out为储 气设备自耗率、储气效率、放气效率。为t时段储热、放热功率, HSin,max、HSout,max分别为储热设备的最大储热、放热功率;为t时段储热设 备的储气量,为t-1时段储热设备的储气量;CHS,min、CHS,max分别为储热设 备的最小、最大储热容量,ηCHS、ηHS,in、ηHS,out为储热设备自耗率、储热效率、 放热效率。
步骤3:建立需求响应模型。
将需求侧可响应电负荷分为可平移电负荷和可中断电负荷,对于可中断 负荷考虑其中断成本,需求侧响应电负荷与当前时段电负荷需求之和必须小 于当前时段允许的最大电负荷,电负荷预测值与需求响应后电负荷的关系, 即:
Pt DR=Pt DR,inte+Pt DR,shif
式中:t为调度时间;Pt DR、Pt DR,inte、Pt DR,shif、Pt LD,max为t时段电负荷预 测值、需求侧响应电负荷、需求侧转移电负荷、***允许的最大电负荷,Pinte,max为调度时间内最大可中断电负荷。为t时段最大可中断和可转出的 电负荷比例。Pt DR,shif为正代表转出可平移负荷,反之,代表转入可平移负荷。
需求侧可响应热负荷约束如下:
式中:t为调度时间;Ht DR、αt HLD、和分别为时间段t热 负荷预测值、需求侧响应热负荷、需求侧可响应热负荷比例、***允许的最 大热负荷和考虑需求响应后的电负荷;HDR,max为调度时间内最大可中断热负 荷;Nt为整个调度时间段。
步骤4:在满足***安全约束的前提下,以基础场景的运行成本最小为目 标函数,建立考虑联合热电需求响应的两阶段可调鲁棒模型。
(4.1)目标函数
本发明所提出的两阶段鲁棒优化模型的目标为在满足***安全约束的前 提下,最小化基础场景的运行成本,由于基础场景下不允许失负荷,因此目 标函数中不含有失负荷相关变量,其目标函数和相关等式约束如下:
式中:t为调度时间;分别为购电、气费用;分 别为第p台CHP的开机、关机费用;为弃风惩罚费用,为弃光惩 罚费用;为售电收益,CE,cc为需求侧响应中断电负荷的补偿成本;为购 电功率,为购气功率,分别为第i台风机、第j组光伏电池在 t时段的弃风、弃光功率;Pt out为售电功率,Pt DR,inte为需求侧中断电负荷。Δt为 调度时间间隔,Nt为调度时段数。分别为单位购电、购气、 售电价格,分别为单位弃风、弃光的惩罚价格,为单位需求侧响 应中断电负荷的补偿价格。为第p台CHP在t时段的启停状态(1代表开机, 0代表停机),分别为第p台CHP的开机、停机一次的费用。NWT、 NPV、NCHP分别为风机、光伏电池、电锅炉、CHP的数量。
(4.2)约束条件
为保证***安全运行,本节采用一个双层max-min模型识别造成***运 行最不安全即最大违反安全规定值(security violation)的最恶劣场景, 最恶劣场景下***的最大违反安全规定值需小于预先设定值εRO,εRO值的设 定与预先确定的***安全级别有关,以确保园区综合能源***的安全运行。 在不确定情况下的最恶劣场景造成违反***安全规定的最大值、风力出力和 光伏出力的不确定性集合、***能量平衡约束、松弛变量恒大于零约束、机 组出力约束及可再生能源出力不确定下的机组纠错爬坡约束、储能约束、与上级网络功率交换约束、需求侧响应约束和弃风弃光约束如下:
v1t,v2t,v3t,v4t≥0
-Pt u,in≤-Pin,min,Pt u,in≤Pin,max:(λ9_1,t,λ9_2,t)
Pt u,DR=Pt u,DR,inte+Pt u,DR,shif:(λ10_6,t)
式中:t为调度时间;Pt u,E,out、分别为弃风、弃光、需求响应后电负荷、P2G消耗电功率、电锅炉消耗电功率、园区售电功率、CHP机组出力、燃气轮机出力根据不同风电出力和光伏出力调整后的实时值;(·)u为在风光出力实时变化下对应的调整后变量;v1,t和v2,t为电力 平衡约束松弛变量;v3,t和v4,t为热量平衡约束松弛变量;ΩWT、ΩPV分别为风电、光 伏出力不确定合集;分别为风电、光伏出力与预测值的偏差,为 风电、光伏出力偏差值与预测值的比例;为不确定合集中的0-1 变量;Δi、Δj分别为风电、光伏出力不确定性预算值;为第p台 CHP机组上爬坡和下爬坡纠错能力;为第q台燃气轮机机组上爬坡和 下爬坡纠错能力;λ(·)为约束条件对应的对偶变量。
步骤5:得到园区综合能源***日前经济调度的两阶段可调鲁棒优化模型的抽 象表达式。
提出的两阶段鲁棒优化模型,涉及电、气、热三个***,涉及的约束较多, 并且含有不确定参数,同时是一个非线性的混合整数规划问题。为便于讨论,本 章所提出的两阶段鲁棒优化调度模型可采用抽象形式的鲁棒优化模型式如下:
s.t. Ax+By≤b
式中:x代表与CHP、燃气轮机相关的机组启停状态,y、z分别代表基础场 景和根据风光出力变换调整的***其余机组调度出力,u为与风电、光伏出力不确 定性相关的不确定变量,cb、cg、A、B、b、f、h、C、D、E、F、G可通 过步骤4中的目标函数和约束条件得出。
步骤6:建立最恶劣场景识别的最大最小子问题。
子问题为识别最恶劣场景的问题,通过步骤4所示的双层的max-min问题, 寻找到造成***最大违反安全规定值的场景,即确定最恶劣场景中不确定量的具 体取值,再通过对偶理论将双层max-min问题转化为如下所示的单层双线性的极 大值优化子问题。
ΔD=ΔD1+ΔD2
s.t. λ9_1,t≤0,λ9_2,t≤0,λ9_3,t≤0,λ9_4,t≤0,λ9_5,t≤0,λ9_6,t≤0,t∈NT
λ6_2,q,t≤0,λ6_3,q,t≤0,λ6_4,q,t≤0,λ6_6,q,t≤0,λ6_7,q,t≤0,q∈NGT,t∈NT
λ7_2,p,t≤0,λ7_3,p,t≤0,λ7_4,p,t≤0,λ7_6,p,t≤0,λ7_7,p,t≤0,p∈NCHP,t∈NT
λ8_1,t≤0,λ8_2,t≤0,λ8_3,t≤0,λ8_4,t≤0,t∈NT
λ13_1,t≤0,λ13_2,t≤0,λ13_3,t≤0,λ13_4,t≤0,t∈NT
λ4_2,m,t≤0,λ4_3,m,t≤0,t∈NT,m∈NP2G,t∈NT
λ5_2,n,t≤0,λ5_3,n,t≤0,t∈NT,n∈NEB,t∈NT
λ10_1,t≤0,λ10_2,t≤0,λ10_3,t≤0,λ10_4,t≤0,t∈NT
λ11_3,t≤0,λ11_3,t≤0,λ12_1,t≤0,λ12_2,t≤0,t∈NT
λ10_7≤0,λ10_8≤0,λ11_4≤0,λ11_5≤0,t∈NT
s.t. -λ1,t-λ6_1,q,t·bq-λ6_2,q,t+λ6_3,q,t-λ6_4,q,t+1+
λ6_4,q,t+λ6_5,p,t+1-λ6_5,p,t-λ6_6,q,t+λ6_7,q,t≤0
-λ1,t-λ9_1,t+λ9_2,t≤0,t∈NT
λ1,t-λ9_3,t+λ9_4,t≤0,t∈NT
λ1,t+λ10_5,t≤0,t∈NT
-λ10_1,t-λ10_2,t+λ10_3,t+λ10_5,t-λ10_7+λ10_8=0,t∈NT
-λ10_1,t+λ10_4,t+λ10_5,t+λ10_5,t+λ10_6=0,t∈NT
λ1,t+λ12_1,i,t≤0,t∈NT,i∈NWT
λ1,t+λ12_2,j,t≤0,t∈NT,j∈NPV
-λ2,t+λ7_1,p,t≤0,t∈NT,p∈NCHP
-λ2,t+λ6_1,q,t≤0,t∈NT,q∈NGT
-λ2,t+λ4_1,m,t-λ4_2,m,t+λ4_3,m,t≤0,t∈NT,m∈NP2G
λ2,t-λ9_5,t+λ9_6,t≤0,t∈NT
λ2,t+λ8_2,t+λ8_5,t/ηGS,out≤0,t∈NT
-λ2,t+λ8_1,t-λ8_5,t·ηGS,in≤0,t∈NT
-λ8_3,t+λ8_4,t-λ8_5,t+1·ηGS,in+λ8_5,t≤0,t∈NT-λ3,t·ηh+λ13_2,t+λ13_5,t/ηHS,out≤0,t∈NT
λ3,t·ηh+λ13_1,t-λ13_5,t·ηHS,in≤0,t∈NT-λ13_3,t+λ13_4,t-λ13_5,t+1·ηHS,in+λ13_5,t≤0,t∈NT
-λ3,t·ηh+λ5_1,n,t-λ5_2,n,t+λ5_3,n,t≤0,t∈NT,n∈NEB-λ3,t·ηh+λ7_8,n,t≤0,t∈NT,n∈NEB
λ3,t+λ11_1,t≤0,t∈NTλ11_1,t-λ11_2,t+λ11_3,t-λ11_4+λ11_5,t≤0,t∈NT
-1≤λ1,t≤1,-1≤λ3,t≤1,t∈NT
步骤7:利用CCG(column and constraint generation)法求解考虑联合热电 需求响应的两阶段可调鲁棒模型。
(7.1)主问题
Ax+By≤b
(7.2)子问题
(7.3)CCG求解步骤
步骤(7.3.1):令迭代计数器s=0,设置***允许的违反安全规定最大 值εRO。
步骤(7.3.2):求解主问题,若有解,得到***机组启停状态x和机组 出力安排y,进行步骤(7.3.3);反之,停止迭代并输出无解。
步骤(7.3.3):根据步骤(7.3.2)中求解得到的x和y,求解最大最小 子问题,找到导致最大可能违反安全规定值的最恶劣场景下的风力和光伏出 力大小。
步骤(7.3.4):如果步骤(7.3.3)中求解出的最大可能违反安全规定 值小于εRO,则x和y是最终优化方案并停止迭代;反之,令s=s+1,根据步 骤(7.3.3)中求解出来的最恶劣场景下风电和光伏出力值向主问题中增加以下所示的CCG约束,返回步骤(7.3.2)。
fTvs≤εRO
步骤8输入园区综合能源***数据还包括园区综合能源***的具体组成、 能源价格、各能源转换设备的设备参数及取值、需求响应比例、基础场景和 最恶劣场景下的新能源出力波动情况、最大违反安全规定值和新能源出力与 负荷的预测值等,采用商业求解器Gurobi对园区综合能源***两阶段可调鲁 棒优化运行模型进行求解,得出鲁棒优化调度结果。
下面通过具体实施例详细说明本发明效果。
(1)算例介绍。
如图2所示为园区综合能源***的具体组成。仿真选取包含风、光、气、 储和考虑电转气、电转热技术的多能互补园区***模型。该***包含燃气轮 机、风力发电机、电锅炉、P2G设备、储热设备、储气设备各一台,CHP两台 和一组光伏电池。溴冷机的制热系数和烟气回收率分别为0.9、1.2,燃 气轮机、CHP、电锅炉的一次开机和停机成本分别为:3.5、1.94、2.74元。 假设CHP、燃气轮机初始状态为停运状态,储气设备初始储气量为10m3,储热 设备初始储热量为100kW·h,储气、储热设备自耗率为0.01。园区多能源综 合***冬季典型日的风机、光伏出力、电负荷和热负荷的预测值如图4所示。
(2)实施例场景描述。
为验证提出的考虑热电联合需求响应的两阶段鲁棒可调优化模型的有效 性,设置表1所示7种调度运行方式。
表1 7种调度运行方式
将风力、光伏预测出力调整为原来2倍,与此同时考虑风光出力的不确 定性,得到新的运行方案1-7。违反安全规定值阈值设置为0,即任何场景下 ***都不允许失负荷和过负荷,假设风电出力和光伏出力的不确定预算为24。
将电负荷增加为原来的2.4倍,与此同时考虑风光出力的不确定性,得到 新的运行方案8-14,违反安全规定值设置为0。
(3)实施例结果分析。
表2给出了不同不确定性比例下运行方案1-7的运行成本和违反安全规定 值,从中可以得到:随着风电、光伏出力不确定比例的加大,园区的运行成 本增加,园区可能出现的最大违反安全规定的值也越大。但是随着不断加入 储热设备、P2G及其储气设备等能源转换设备,园区能不断实时调整储能状 态来应对风电光伏出力的实时变化,降低***运行风险,考虑需求响应可以 灵活应对***的不确定性,考虑联合热电需求响应,进一步降低了园区基础 场景下的运行成本。
表2不同不确定性比例下运行方案1-7的运行成本和违反安全规定值
表3给出了不同风电、光伏不确定性预算下运行方案1的迭代次数和违 反安全规定值,从中可以得到:可再生能源不确定性预算值越大,得到满足 ***需求的优化方案所需的迭代次数越少,同时***可能出现的违反安全规 定值也越大。通过调节不确定性预算值能够调节***的鲁棒性。
表3不同风电、光伏不确定性预算下运行方案1的迭代次数和违反安全规定值
表4给出了不同不确定性比例下运行方案8-14的运行成本和违反安全规 定值,从中可以得到:当***负荷需求明显大于可再生能源出力时,可再生 能源出力能被***完全消纳,气转电优于电转气,P2G设备不工作,储能设 备的引入不能提升***的鲁棒性。当可再生能源出力不确定比例为0.2时,对 比方案8-10,可知***含有的设备越全面,***在基础场景下的运行成本越 低,对比方案10-14可知,需求侧响应的引入能进一步降低***的运行成本。 对比方案11-14在不同可再生能源不确定比例的运行成本和违反安全规定值可知,联合热电需求响应可以通过降低热负荷来减少电力负荷损失,极大地 增强了园区综合能源***的鲁棒性、灵活性和经济性。
表4不同不确定性比例下运行方案8-14的运行成本和违反安全规定值
以上所述,仅为本发明的具体实施例,但并不因此限值本发明的专利保 护范围,凡是利用本发明说明书以及附图内容进行等效变化或替换,直接或 间接运用到其他相关技术领域,都应包括在本发明的保护范围之内。
Claims (9)
1.提出一种考虑风光不确定性的园区综合能源***日前经济调度方法,其特征在于,包括以下步骤:
步骤1:确定园区综合能源***的具体组成,包含引入的新能源形式和具体的设备组成;
步骤2:建立园区内各能源转换设备模型;
步骤3:建立需求响应模型;
步骤4:在满足***安全约束的前提下,以基础场景的运行成本最小为目标函数,建立考虑联合热电需求响应的两阶段可调鲁棒模型;
步骤5:得到园区综合能源***日前经济调度的两阶段可调鲁棒优化模型的抽象表达式;
步骤6:建立最恶劣场景识别的最大最小子问题;
步骤7:利用CCG法求解考虑联合热电需求响应的两阶段可调鲁棒模型;
步骤8:输入园区综合能源***能源接入、新能源出力数据、各设备参数、运行参数,采用商业求解器Gurobi对考虑风光不确定性的园区综合能源***日前经济调度两阶段鲁棒优化模型进行求解,得到其调度策略。
2.根据权利要求1所述的考虑风光不确定性的园区综合能源***日前经济调度方法,其特征在于,步骤1所述园区综合能源***的具体组成如下:
(1)接入园区综合能源***的新能源形式为:风电和光伏发电;
(2)引入园区综合能源***的能源转换设备有:电转气设备、电锅炉、燃气轮机、热电联产机组、储气/储热设备。
3.根据权利要求1所述的考虑风光不确定性的园区综合能源***日前经济调度方法,其特征在于,步骤2所述各能源转换设备模型如下;
(1)电转气设备模型
(2)电锅炉模型
式中:t为调度时间;n为电锅炉索引;和分别为第n台电锅炉在t时段的耗电功率和产热功率;为第n台电锅炉的电热转换效率, 分为第n台电锅炉最小、最大制热功率;为第n台电锅炉在t时段的启停状态(1代表开机,0代表停机)。
(3)燃气轮机模型
式中:t为调度时间;q为燃气轮机索引;和分别表示燃气轮机的发电功率和耗气功率;F(·)表示燃气轮机能耗曲线;和分别表示燃气轮机的开机和关机所需天然气消耗量;aq、bq和cq表示F(·)的燃气系数;LHANG为天然气低热值,取9.7kWh/m3;分为第q台燃气轮机最小、最大发电功率。表示第q台燃气轮机在t时段的发电功率;为第q台燃气轮机在t时段的启停状态(1代表开机,0代表停机),为第q台燃气轮机在t时段的启停状态;为第q台燃气轮机的上爬坡率、下爬坡率,为第q台燃气轮机在t-1时段内连续开机、停机时间, 为第q台燃气轮机在时段t内的最小开机、停机时间。
(4)热电联产机组模型
式中:t为调度时间;p为热电联产机组索引;和分别表示热电联产机组(CHP)的产热功率和耗气功率;为微燃机t时段的发电功率、发电效率,为散热损失率,和分别为溴冷机的制热系数和烟气回收率;LHANG为天然气低热值,取9.7kWh/m3;分为第p台CHP机组最小、最大发电功率;表示第p台热电联产机组在t时段的发电功率;为第p台CHP机组在t时段的启停状态(1代表开机,0代表停机),为第P台热电联产机组在t时段的启停状态;为第p台CHP机组的上爬坡率、下爬坡率;为第p台CHP机组在t-1时段内连续开机、停机时间,为第p台CHP机组在时段t内的最小开机、停机时间。
(5)储能设备模型
ESS(t)=(1-σS)·ESS(t-1)+ESin(t)·ηin-ESout(t)/ηout
式中:t为调度时间;ESS(t)、ESin(t)、ESout(t)分别为t时段储能设备的储能量、存储能量功率和释放能量功率;ESS(t-1)表示t-1时段储能设备的储能量;σS为储能***的自耗率;ηin、ηout分别为储能设备储能、放能效率;
为t时段储气、放气功率,GSin,max、GSout,max分别为储气设备的最大储气、放气功率;为t时段储气设备的储气量,为t-1时段储气设备的储气量;CGS,min、CGS,max分别为储气设备的最小、最大储气容量,ηCGS、ηGS,in、ηGS,out为储气设备自耗率、储气效率、放气效率。为t时段储热、放热功率,HSin,max、HSout,max分别为储热设备的最大储热、放热功率;为t时段储热设备的储气量,为t-1时段储热设备的储气量;
CHS,min、CHS,max分别为储热设备的最小、最大储热容量,ηCHS、ηHS,in、ηHS,out为储热设备自耗率、储热效率、放热效率。
4.根据权利要求1所述的考虑风光不确定性的园区综合能源***日前经济调度方法,其特征在于,步骤3所述需求响应模型如下:
Pt DR=Pt DR,inte+Pt DR,shif
式中:t为调度时间;Pt DR、Pt DR,inte、Pt DR,shif、Pt LD,max为t时段电负荷预测值、需求侧响应电负荷、需求侧转移电负荷、***允许的最大电负荷,Pinte,max为调度时间内最大可中断电负荷。为t时段最大可中断和可转出的电负荷比例。Pt DR,shif为正代表转出可平移负荷,反之,代表转入可平移负荷。
5.根据权利要求1所述的考虑风光不确定性的园区综合能源***日前经济调度方法,其特征在于,步骤4所述考虑联合热电需求响应的两阶段可调鲁棒模型如下:
(1)目标函数
式中:t为调度时间;分别为购电、气费用;分别为第p台CHP的开机、关机费用;为弃风惩罚费用,为弃光惩罚费用;为售电收益,CE,cc为需求侧响应中断电负荷的补偿成本;为购电功率,为购气功率,分别为第i台风机、第j组光伏电池在t时段的弃风、弃光功率;Pt out为售电功率,Pt DR,inte为需求侧中断电负荷。△t为调度时间间隔,Nt为调度时段数。分别为单位购电、购气、售电价格,分别为单位弃风、弃光的惩罚价格,为单位需求侧响应中断电负荷的补偿价格。为第p台CHP在t时段的启停状态(1代表开机,0代表停机),分别为第p台CHP的开机、停机一次的费用。NWT、NPV、NCHP分别为风机、光伏电池、电锅炉、CHP的数量。
(2)约束条件
式中:t为调度时间;(·)u为在风光出力实时变化下对应的调整后变量;v1,t和v2,t为电力平衡约束松弛变量;v3,t和v4,t为热量平衡约束松弛变量;ΩWT、ΩPV分别为风电、光伏出力不确定合集;分别为风电、光伏出力与预测值的偏差,为风电、光伏出力偏差值与预测值的比例; 为不确定合集中的0-1变量;△i、△j分别为风电、光伏出力不确定性预算值;为第p台CHP机组上爬坡和下爬坡纠错能力;
8.根据权利要求1所述的考虑风光不确定性的园区综合能源***日前经济调度方法,其特征在于,步骤7所述利用CCG法求解考虑联合热电需求响应的两阶段可调鲁棒模型如下:
(1)主问题
Ax+By≤b
(2)子问题
(3)CCG求解步骤
步骤1:令迭代计数器s=0,设置***允许的违反安全规定最大值εRO。
步骤2:求解主问题,若有解,得到***机组启停状态x和机组出力安排y,进行步骤3;反之,停止迭代并输出无解。
步骤3:根据步骤2中求解得到的x和y,求解最大最小子问题,找到导致最大可能违反安全规定值的最恶劣场景下的风力和光伏出力大小。
步骤4:如果步骤3中求解出的最大可能违反安全规定值小于εRO,则x和y是最终优化方案并停止迭代;反之,令s=s+1,根据步骤3中求解出来的最恶劣场景下风电和光伏出力值 向主问题中增加以下所示的CCG约束,返回步骤2。
fTvs≤εRO
9.根据权利要求1所述的考虑气电联合需求响应的气电综合能源***配网优化运行方法,其特征在于,步骤8所述园区综合能源***数据还包括园区综合能源***的具体组成、能源价格、各能源转换设备的设备参数及取值、需求响应比例、基础场景和最恶劣场景下的新能源出力波动情况、最大违反安全规定值和新能源出力与负荷的预测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010773898.1A CN113256045B (zh) | 2020-08-04 | 2020-08-04 | 考虑风光不确定性的园区综合能源***日前经济调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010773898.1A CN113256045B (zh) | 2020-08-04 | 2020-08-04 | 考虑风光不确定性的园区综合能源***日前经济调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113256045A true CN113256045A (zh) | 2021-08-13 |
CN113256045B CN113256045B (zh) | 2023-04-07 |
Family
ID=77220096
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010773898.1A Active CN113256045B (zh) | 2020-08-04 | 2020-08-04 | 考虑风光不确定性的园区综合能源***日前经济调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113256045B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113705991A (zh) * | 2021-08-17 | 2021-11-26 | 国网四川省电力公司技能培训中心 | 一种多能源园区的建立及低碳调度方法 |
CN114219221A (zh) * | 2021-11-19 | 2022-03-22 | 贵州电网有限责任公司 | 一种考虑多重不确定性的综合能源***日前经济协调与调度方法 |
CN114519449A (zh) * | 2021-12-01 | 2022-05-20 | 中国华能集团有限公司河北雄安分公司 | 一种园区能源***运行优化方法 |
CN114784797A (zh) * | 2022-04-25 | 2022-07-22 | 东南大学溧阳研究院 | 一种计及多重不确定性的住宅综合能源***的热电优化日前调度方法 |
CN114912714A (zh) * | 2022-06-17 | 2022-08-16 | 国网江苏省电力有限公司苏州供电分公司 | 一种考虑雷电气候下风光出力不确定性的低碳经济调度方法和*** |
CN117039901A (zh) * | 2023-08-10 | 2023-11-10 | 国网山东省电力公司威海供电公司 | 基于低压柔性直流互联的台区不确定性就地管理调度方法 |
WO2024103450A1 (zh) * | 2022-11-15 | 2024-05-23 | 云南电网有限责任公司电力科学研究院 | 多能***的调度方法、计算机设备和计算机可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108832665A (zh) * | 2018-07-04 | 2018-11-16 | 四川大学 | 一种考虑风电不确定性的电热综合***分布式鲁棒协调优化调度模型 |
CN110009152A (zh) * | 2019-04-03 | 2019-07-12 | 东南大学 | 一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法 |
CN110555595A (zh) * | 2019-08-09 | 2019-12-10 | 四川大学 | 一种基于能量枢纽的沼-风-光全可再生能源***及其方法 |
CN111355265A (zh) * | 2020-04-10 | 2020-06-30 | 中国人民解放军国防科技大学 | 一种微电网能量两阶段鲁棒优化方法及*** |
-
2020
- 2020-08-04 CN CN202010773898.1A patent/CN113256045B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108832665A (zh) * | 2018-07-04 | 2018-11-16 | 四川大学 | 一种考虑风电不确定性的电热综合***分布式鲁棒协调优化调度模型 |
CN110009152A (zh) * | 2019-04-03 | 2019-07-12 | 东南大学 | 一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法 |
CN110555595A (zh) * | 2019-08-09 | 2019-12-10 | 四川大学 | 一种基于能量枢纽的沼-风-光全可再生能源***及其方法 |
CN111355265A (zh) * | 2020-04-10 | 2020-06-30 | 中国人民解放军国防科技大学 | 一种微电网能量两阶段鲁棒优化方法及*** |
Non-Patent Citations (4)
Title |
---|
CHUAN HE 等: "Co-optimization scheduling of interdependent power and gas systems with electricity and gas uncertainties" * |
JUN LU 等: "Robust day-ahead coordinated scheduling of multi-energy systems with integrated heat-electricity demand response and high penetration of renewable energy" * |
单福州 等: "基于改进两阶段鲁棒优化的区域综合能源***经济调度" * |
税月 等: "考虑风电不确定性的电气能源***两阶段分布鲁棒协同调度" * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113705991A (zh) * | 2021-08-17 | 2021-11-26 | 国网四川省电力公司技能培训中心 | 一种多能源园区的建立及低碳调度方法 |
CN113705991B (zh) * | 2021-08-17 | 2023-08-29 | 国网四川省电力公司技能培训中心 | 一种多能源园区的建立及低碳调度方法 |
CN114219221A (zh) * | 2021-11-19 | 2022-03-22 | 贵州电网有限责任公司 | 一种考虑多重不确定性的综合能源***日前经济协调与调度方法 |
CN114519449A (zh) * | 2021-12-01 | 2022-05-20 | 中国华能集团有限公司河北雄安分公司 | 一种园区能源***运行优化方法 |
CN114784797A (zh) * | 2022-04-25 | 2022-07-22 | 东南大学溧阳研究院 | 一种计及多重不确定性的住宅综合能源***的热电优化日前调度方法 |
CN114784797B (zh) * | 2022-04-25 | 2024-01-19 | 东南大学溧阳研究院 | 一种计及多重不确定性的住宅综合能源***的热电优化日前调度方法 |
CN114912714A (zh) * | 2022-06-17 | 2022-08-16 | 国网江苏省电力有限公司苏州供电分公司 | 一种考虑雷电气候下风光出力不确定性的低碳经济调度方法和*** |
CN114912714B (zh) * | 2022-06-17 | 2023-11-07 | 国网江苏省电力有限公司苏州供电分公司 | 一种考虑雷电气候下风光出力不确定性的低碳经济调度方法和*** |
WO2024103450A1 (zh) * | 2022-11-15 | 2024-05-23 | 云南电网有限责任公司电力科学研究院 | 多能***的调度方法、计算机设备和计算机可读存储介质 |
CN117039901A (zh) * | 2023-08-10 | 2023-11-10 | 国网山东省电力公司威海供电公司 | 基于低压柔性直流互联的台区不确定性就地管理调度方法 |
CN117039901B (zh) * | 2023-08-10 | 2024-07-05 | 国网山东省电力公司威海供电公司 | 基于低压柔性直流互联的台区不确定性就地管理调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113256045B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113256045B (zh) | 考虑风光不确定性的园区综合能源***日前经济调度方法 | |
Zhou et al. | Optimized operation method of small and medium-sized integrated energy system for P2G equipment under strong uncertainty | |
CN109004686B (zh) | 一种考虑冰蓄冷空调多模式的冷热电联供型微网*** | |
CN108717594B (zh) | 一种冷热电联供型多微网***经济优化调度方法 | |
CN104734168B (zh) | 一种基于电热联合调度的微电网运行优化***及方法 | |
CN111950807B (zh) | 计及不确定性与需求响应的综合能源***优化运行方法 | |
CN112464477A (zh) | 计及需求响应的多能耦合综合能源运行仿真方法 | |
CN111614121A (zh) | 考虑需求响应的含电动汽车的多能源园区日前经济调度方法 | |
CN106786793A (zh) | 一种基于鲁棒优化的冷热电联供型微网运行方法 | |
CN112701687B (zh) | 考虑价格型联合需求响应的气电配网***鲁棒优化运行方法 | |
CN107358345B (zh) | 计及需求侧管理的分布式冷热电联供***优化运行方法 | |
CN111737884B (zh) | 一种含多种清洁能源微能源网多目标随机规划方法 | |
CN104537443A (zh) | 一种热电联供型微网经济协调优化调度方法 | |
CN113193602A (zh) | 含低热值发电和分布式电源的配电网优化运行***及方法 | |
CN110620403A (zh) | 计及可再生能源的能源***协同运行日前调度方法及*** | |
CN112836882B (zh) | 考虑设备负载率变化的区域综合能源***运行优化方法 | |
CN113779783A (zh) | 计及多重不确定性的区域综合能源***规划与运行联合优化方法 | |
CN112446546B (zh) | 一种考虑能源可靠性的综合能源***两阶段优化配置方法 | |
CN111160636A (zh) | 一种cchp型微电网调度优化方法 | |
CN115099007B (zh) | 基于综合成本-能耗曲线的综合能源***优化运行方法 | |
CN115860413A (zh) | 基于负荷需求响应与双层可调鲁棒优化的并网冷热电联产微网***经济调度方法 | |
CN117081143A (zh) | 促进分布式光伏就地消纳的园区综合能源***协调优化运行方法 | |
CN111030101A (zh) | 一种基于多元化大数据清洁能源消纳联动调控方法及*** | |
CN114362152B (zh) | 综合能源***多时间尺度调度方法 | |
CN114091917A (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 |