CN111049184A - 一种多能互补园区响应电力系统风电消纳需求的计算方法 - Google Patents
一种多能互补园区响应电力系统风电消纳需求的计算方法 Download PDFInfo
- Publication number
- CN111049184A CN111049184A CN201911362865.1A CN201911362865A CN111049184A CN 111049184 A CN111049184 A CN 111049184A CN 201911362865 A CN201911362865 A CN 201911362865A CN 111049184 A CN111049184 A CN 111049184A
- Authority
- CN
- China
- Prior art keywords
- power
- energy complementary
- ith
- energy
- complementary park
- 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
- 230000000295 complement effect Effects 0.000 title claims abstract description 155
- 238000004364 calculation method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 15
- 230000004044 response Effects 0.000 claims abstract description 6
- 230000009194 climbing Effects 0.000 claims description 72
- 238000005457 optimization Methods 0.000 claims description 62
- 238000003860 storage Methods 0.000 claims description 45
- 238000010521 absorption reaction Methods 0.000 claims description 28
- 238000001816 cooling Methods 0.000 claims description 26
- 238000007599 discharging Methods 0.000 claims description 18
- 238000010438 heat treatment Methods 0.000 claims description 16
- 238000010248 power generation Methods 0.000 claims description 10
- 238000004134 energy conservation Methods 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 230000008859 change Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
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/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
-
- 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)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种多能互补园区响应电力系统风电消纳需求的计算方法,属于综合能源系统的运行控制技术领域。本发明计算方法,在日前计算多能互补园区次日的联络线初始计划即联络线功率基线;在电力系统存在风电消纳需求时,在日前计算多能互补园区对电力系统发布的风电消纳需求的最大响应程度;在电力系统存在风电消纳需求时,在日前计算不同程度风电消纳时的联络线功率计划曲线、多能互补园区内部的能源设备运行计划和最小调节成本。电力系统可根据这些计算结果向多能互补园区下发联络线功率计划。本发明能够充分利用多能互补园区内部的能源设备响应电力系统风电消纳需求,从而提升风电利用率,减少弃风的产生。
Description
技术领域
本发明涉及一种多能互补园区响应电力系统风电消纳需求的计算方法,属于综合能源系统的运行控制技术领域。
背景技术
近年来,电力系统中风电的比例逐年上升。风电的随机性、间歇性、波动性和逆调峰性等特点使得弃风问题极为突出。
食品、材料加工、制药、机械制造等工业用户同时需要冷、热、电等多种能量形式进行生产作业,这些工厂聚集形成多能互补园区。其中冷、热能量在园区内平衡,由园区内的供热、供冷设备供给,这些供冷、供热设备常常与电能存在耦合关系。电能既可利用园区内的发电机产生,又可通过联络线从电力系统获取。由于不同的能流形式之间存在互补性,园区可通过调节不同特性的供冷、供热设备间的功率分配,在保证冷、热功率平衡的前提下调节联络线功率,同时园区内各种发电设备的发电功率也可实时调节,故多能互补园区联络线的电功率存在一定的灵活性,可作为电力系统的灵活性资源参与辅助服务。
发明内容
本发明涉及一种多能互补园区响应电力系统风电消纳需求的计算方法,以解决已有技术中存在的问题,计算多能互补园区可响应的不同风电消纳比例下的联络线功率计划曲线、多能互补园区内部的能源设备运行计划和最小调节成本,通过合理安排热电联产机组、燃气锅炉、吸收式制冷机、电锅炉和电制冷机能源设备的出力,以及合理控制蓄电池的充/放电时间和充/放电功率,可以改变园区联络线功率,从而将多能互补园区作为灵活性资源,响应电力系统风电消纳需求。
本发明提出的多能互补园区响应电力系统风电消纳需求的计算方法,包括以下步骤:
(1)设定多能互补园区中的能源设备包括发电机、热电联产机组、燃气锅炉、吸收式制冷机组、电锅炉、电制冷机和蓄电池,该多能互补园区与电力系统通过联络线相连,建立一个第一优化模型,第一优化模型用于计算多能互补园区向电力系统上报的联络线功率基线,作为电力系统日前调度的功率基线,第一优化模型以不加调控时多能互补园区自身总运行成本c0最小为目标,表达式如下:
其中,x为该第一优化模型中所有待求解变量构成的列向量,即:
为第iG台发电机在调度时刻t发出的有功功率;和分别为第iCHP台热电联产机组在调度时刻t的有功功率和供热功率;为第iGB台燃气锅炉在调度时刻t的供热功率;为第iAC台吸收式制冷机组在调度时刻t的供冷功率;和分别为第iEB台电锅炉在调度时刻t的耗电功率和供热功率;和分别为第iEC台电制冷机在调度时刻t的耗电功率和供冷功率;和分别为第iES台蓄电池在调度时刻t的充电功率和放电功率;为第iES台蓄电池在调度时刻t的充电状态的0-1变量,代表蓄电池在调度时刻t处在充电状态,代表蓄电池在调度时刻t未处在充电状态;为描述第iES台蓄电池在调度时刻t的放电状态的0-1变量,代表蓄电池在调度时刻t处在放电状态,代表蓄电池在调度时刻t未处在放电状态;为第iES台蓄电池在调度时刻t的充电转换状态的0-1变量,即代表蓄电池在调度时刻t-1未在充电、在调度时刻t处在充电状态,代表蓄电池在调度时刻t-1在充电、调度时刻t处在未充电状态;为第iES台蓄电池在调度时刻t的放电转换状态的0-1变量,即代表蓄电池在调度时刻t-1未在放电、调度时刻t处在放电状态,代表蓄电池在调度时刻t-1处在放电、调度时刻t未在放电状态;为调度时刻t的联络线功率,即多能互补园区与电力系统交换的电功率,以电功率流入多能互补园区为正方向;上标T为向量转置;γ为电力系统所有调度时刻t构成的集合;SG为多能互补园区内所有发电机构成的集合;为第iG台发电机在一个调度时刻发出单位有功功率所需成本,对于分布式光伏发电装置和分布式风电机组该值可取为0;SGB为所有燃气锅炉构成的集合;为第iGB台燃气锅炉在一个调度时刻内维持单位供热功率所需成本;SAC为所有吸收式制冷机构成的集合;为第iAC台吸收式制冷机组在一个调度时刻内维持单位供冷功率所需成本;SCHP为所有热电联产机组构成的集合;为第iCHP台热电联产机组在一个调度时刻内维持单位发电功率所需成本;为第iCHP台热电联产机组在一个调度时刻内维持单位供热功率所需成本;为调度时刻t的联络线电价;ΔT为相邻两个调度时刻的时间间隔;
第一优化模型的约束条件包括:
(1-1)多能互补园区中发电机的有功功率范围和爬坡约束:
(1-2)多能互补园区中热电联产机组的有功功率范围、供热功率范围和热电联产机组的有功功率爬坡约束:
(1-3)多能互补园区中燃气锅炉供热功率范围及爬坡约束:
(1-4)多能互补园区中吸收式制冷机的供冷功率范围及爬坡约束:
(1-5)多能互补园区中电锅炉的供热功率范围及爬坡约束:
(1-6)多能互补园区中电制冷机的供冷功率范围及爬坡约束:
(1-7)多能互补园区中蓄电池的运行约束:
其中,和为第iES台蓄电池的充电功率上限和下限;和为第iES台蓄电池的放电功率上限和下限;为一个调度日中第i台蓄电池的最大充放电转换次数;和为第iES台蓄电池中可储存的最大能量和最小能量;上述所有参数均在对应设备的说明书或操作规程中得到;为第iES台蓄电池在该调度日开始时储存的能量,由蓄电池前一日的调度结果得到;
(1-8)多能互补园区中的电能、热能、冷能的能量守恒约束:
(1-9)多能互补园区中冷、热负荷的惯性及室内温度约束:
其中,为热负荷的热容,和分别为热负荷t和t-1时刻的室内温度,为热负荷的热导,SHL为所有热负荷构成的集合,为t时刻的环境温度,和分别为热负荷室内温度最小值和最大值,为冷负荷的热容,和分别为冷负荷t和t-1时刻的室内温度,为冷负荷的热导,SCL为所有冷负荷构成的集合,和分别为冷负荷室内温度最小值和最大值。
采用分支定界法,求解由上述目标函数和约束条件组成的第一优化模型,得到多能互补园区与电力系统之间的联络线在调度时刻t的有功功率的最优值即多能互补园区联络线功率基线,以及目标函数最优值c0,即多能互补园区最小运行成本,多能互补园区将联络线功率基线上报电力系统;
(2)多能互补园区接收电力系统下发的调度时段t风电消纳需求 为电力系统期望联络线在调度时段t的有功功率高于联络线功率基线的有功功率值;根据该风电消纳需求计算多能互补园区计算可提供的风电消纳能力,具体方法如下:
(2-1)选取计算步长ε,ε取值为0~1,置循环次数k初始值1;
(2-2)求解如下的第二优化模型,该第二优化模型用于在调度时段t,当联络线功率为时,求多能互补园区运行费用最小时多能互补园区中的能源设备的运行计划;第二优化模型的目标函数与第一优化模型的目标函数相同;
第二优化模型的约束条件包括:
a、多能互补园区中发电机的有功功率范围和爬坡约束:与第一优化模型中约束条件(1-1)相同;
b、多能互补园区中热电联产机组的有功功率和供热功率范围约束和热电联产机组的有功功率爬坡约束:与第一优化模型中约束条件(1-2)相同;
c、多能互补园区中燃气锅炉供热功率范围及爬坡约束:与第一优化模型中约束条件(1-3)相同;
d、多能互补园区中吸收式制冷机的供冷功率范围及爬坡约束:与第一优化模型中约束条件(1-4)相同;
e、多能互补园区中电锅炉的供热功率范围及爬坡约束:与第一优化模型中约束条件(1-5)相同;
f、多能互补园区中电制冷机的供冷功率范围及爬坡约束:与第一优化模型中约束条件(1-6)相同;
g、多能互补园区中蓄电池的运行约束:与第一优化模型中约束条件(1-7)相同;
h、多能互补园区中的电能、热能、冷能的能量守恒约束:
i、多能互补园区中冷、热负荷的惯性及室内温度约束:与第一优化模型中约束条件(1-9)相同;
采用分支定界法,求解由上述目标函数和约束条件组成的第二优化模型,若第二优化模型可解,则得到目标函数最优值ck,以及多能互补园区与电力系统之间的联络线在调度时刻t的有功功率的最优值即多能互补园区向电力系统上报的功率基线,转步骤(2-3);若该模型无解,则结束循环,进行步骤(3);
(2-3)计算风电消纳比例为kε时的多能互补园区的调节成本ck*,ck*=ck-c0,记录ck*,并置k=k+1,对kε进行判断,若kε>1,转步骤(3),若kε≤1,返回步骤(2-2),直到kε>1,得到多个与不同风电消纳比例kε相对应的多能互补园区的调节成本ck*;
(3)将步骤(2-3)的多个与不同风电消纳比例kε相对应的多能互补园区的调节成本ck*上报电力系统,电力系统下发最终采用的多能互补园区的风电消纳比例,实现由多能互补园区对电力系统风电消纳需求的响应。
本发明提出的多能互补园区响应电力系统风电消纳需求的计算方法,其优点是:
本发明多能互补园区响应电力系统风电消纳需求的计算方法,在日前计算多能互补园区次日的联络线初始计划即联络线功率基线;在电力系统存在风电消纳需求时,在日前计算多能互补园区对电力系统发布的风电消纳需求的最大响应程度;在电力系统存在风电消纳需求时,在日前计算不同程度风电消纳时的联络线功率计划曲线、多能互补园区内部的能源设备运行计划和最小调节成本。电力系统可根据这些计算结果向多能互补园区下发联络线功率计划。本发明能够充分利用多能互补园区内部的能源设备响应电力系统风电消纳需求,从而提升风电利用率,减少弃风的产生。
具体实施方式
本发明提出的多能互补园区响应电力系统风电消纳需求的计算方法,包括以下步骤:
(1)计算多能互补园区次日的联络线初始计划即联络线功率基线,一般以园区总供能成本最小化为目标,求解如下优化问题:设定多能互补园区中的能源设备包括发电机、热电联产机组、燃气锅炉、吸收式制冷机组、电锅炉、电制冷机和蓄电池,其中的发电机可以为常规发电机、分布式光伏发电装置或分布式风电机组,该多能互补园区与电力系统通过联络线相连,建立一个第一优化模型,第一优化模型用于计算多能互补园区向电力系统上报的联络线功率基线,作为电力系统日前调度的功率基线,第一优化模型以不加调控时多能互补园区自身总运行成本c0最小为目标,第一优化模型的目标函数为使c0为最小:
其中,x为该第一优化模型中所有待求解变量构成的列向量,即:
为第iG台发电机在调度时刻t发出的有功功率;和分别为第iCHP台热电联产机组在调度时刻t的有功功率和供热功率;为第iGB台燃气锅炉在调度时刻t的供热功率;为第iAC台吸收式制冷机组在调度时刻t的供冷功率;和分别为第iEB台电锅炉在调度时刻t的耗电功率和供热功率;和分别为第iEC台电制冷机在调度时刻t的耗电功率和供冷功率;和分别为第iES台蓄电池在调度时刻t的充电功率和放电功率;为第iES台蓄电池在调度时刻t的充电状态的0-1变量,代表蓄电池在调度时刻t处在充电状态,代表蓄电池在调度时刻t未处在充电状态;为描述第iES台蓄电池在调度时刻t的放电状态的0-1变量,代表蓄电池在调度时刻t处在放电状态,代表蓄电池在调度时刻t未处在放电状态;为第iES台蓄电池在调度时刻t的充电转换状态的0-1变量,即代表蓄电池在调度时刻t-1未在充电、在调度时刻t处在充电状态,代表蓄电池在调度时刻t-1在充电、调度时刻t处在未充电状态;为第iES台蓄电池在调度时刻t的放电转换状态的0-1变量,即代表蓄电池在调度时刻t-1未在放电、调度时刻t处在放电状态,代表蓄电池在调度时刻t-1处在放电、调度时刻t未在放电状态;为调度时刻t的联络线功率,即多能互补园区与电力系统交换的电功率,以电功率流入多能互补园区为正方向;上标T为向量转置;γ为电力系统所有调度时刻t构成的集合;SG为多能互补园区内所有发电机构成的集合;为第iG台发电机在一个调度时刻发出单位有功功率所需成本,对于分布式光伏发电装置和分布式风电机组该值可取为0;SGB为所有燃气锅炉构成的集合;为第iGB台燃气锅炉在一个调度时刻内维持单位供热功率所需成本;SAC为所有吸收式制冷机构成的集合;为第iAC台吸收式制冷机组在一个调度时刻内维持单位供冷功率所需成本;SCHP为所有热电联产机组构成的集合;为第iCHP台热电联产机组在一个调度时刻内维持单位发电功率所需成本;为第iCHP台热电联产机组在一个调度时刻内维持单位供热功率所需成本;为调度时刻t的联络线电价;ΔT为相邻两个调度时刻的时间间隔;
第一优化模型的约束条件包括:
(1-1)多能互补园区中发电机的有功功率范围和爬坡约束:
(1-2)多能互补园区中热电联产机组的有功功率范围、供热功率范围和热电联产机组的有功功率爬坡约束:
(1-3)多能互补园区中燃气锅炉供热功率范围及爬坡约束:
(1-4)多能互补园区中吸收式制冷机的供冷功率范围及爬坡约束:
(1-5)多能互补园区中电锅炉的供热功率范围及爬坡约束:
(1-6)多能互补园区中电制冷机的供冷功率范围及爬坡约束:
(1-7)多能互补园区中蓄电池的运行约束:
其中,和为第iES台蓄电池的充电功率上限和下限;和为第iES台蓄电池的放电功率上限和下限;为一个调度日中第i台蓄电池的最大充放电转换次数;和为第iES台蓄电池中可储存的最大能量和最小能量;上述所有参数均在对应设备的说明书或操作规程中得到;为第iES台蓄电池在该调度日开始时储存的能量,由蓄电池前一日的调度结果得到;
(1-8)多能互补园区中的电能、热能、冷能的能量守恒约束:
(1-9)多能互补园区中冷、热负荷的惯性及室内温度约束:
其中,为热负荷的热容,和分别为热负荷t和t-1时刻的室内温度,为热负荷的热导,SHL为所有热负荷构成的集合,为t时刻的环境温度,和分别为热负荷室内温度最小值和最大值,为冷负荷的热容,和分别为冷负荷t和t-1时刻的室内温度,为冷负荷的热导,SCL为所有冷负荷构成的集合,和分别为冷负荷室内温度最小值和最大值。
采用分支定界法,求解由上述目标函数和约束条件组成的第一优化模型,得到多能互补园区与电力系统之间的联络线在调度时刻t的有功功率的最优值即多能互补园区联络线功率基线,以及目标函数最优值c0,即多能互补园区最小运行成本,多能互补园区将联络线功率基线上报电力系统;
(2)多能互补园区接收电力系统下发的调度时段t风电消纳需求 为电力系统期望联络线在调度时段t的有功功率高于联络线功率基线的有功功率值;根据该风电消纳需求计算多能互补园区计算可提供的风电消纳能力,具体方法如下:
(2-1)选取计算步长ε,ε取值为0~1,本发明的一个实施例中取0.1;置循环次数k初始值1;
(2-2)求解如下的第二优化模型,该第二优化模型用于在调度时段t,当联络线功率为时,求多能互补园区运行费用最小时多能互补园区中的能源设备的运行计划;第二优化模型的目标函数与第一优化模型的目标函数相同;
第二优化模型的约束条件包括:
a、多能互补园区中发电机的有功功率范围和爬坡约束:与第一优化模型中约束条件(1-1)相同;
b、多能互补园区中热电联产机组的有功功率和供热功率范围约束和热电联产机组的有功功率爬坡约束:与第一优化模型中约束条件(1-2)相同;
c、多能互补园区中燃气锅炉供热功率范围及爬坡约束:与第一优化模型中约束条件(1-3)相同;
d、多能互补园区中吸收式制冷机的供冷功率范围及爬坡约束:与第一优化模型中约束条件(1-4)相同;
e、多能互补园区中电锅炉的供热功率范围及爬坡约束:与第一优化模型中约束条件(1-5)相同;
f、多能互补园区中电制冷机的供冷功率范围及爬坡约束:与第一优化模型中约束条件(1-6)相同;
g、多能互补园区中蓄电池的运行约束:与第一优化模型中约束条件(1-7)相同;
h、多能互补园区中的电能、热能、冷能的能量守恒约束:
i、多能互补园区中冷、热负荷的惯性及室内温度约束:与第一优化模型中约束条件(1-9)相同;
采用分支定界法,求解由上述目标函数和约束条件组成的第二优化模型,若第二优化模型可解,则得到目标函数最优值ck,以及多能互补园区与电力系统之间的联络线在调度时刻t的有功功率的最优值即多能互补园区向电力系统上报的功率基线,转步骤(2-3);若该模型无解,则结束循环,进行步骤(3);
(2-3)计算风电消纳比例为kε时的多能互补园区的调节成本ck*,ck*=ck-c0,记录ck*,并置k=k+1,对kε进行判断,若kε>1,转步骤(3),若kε≤1,返回步骤(2-2),直到kε>1,得到多个与不同风电消纳比例kε相对应的多能互补园区的调节成本ck*;
(3)将步骤(2-3)的多个与不同风电消纳比例kε相对应的多能互补园区的调节成本ck*上报电力系统,电力系统下发最终采用的多能互补园区的风电消纳比例,实现由多能互补园区对电力系统风电消纳需求的响应。
Claims (1)
1.一种多能互补园区响应电力系统风电消纳需求的计算方法,其特征在于该方法包括以下步骤:
(1)设定多能互补园区中的能源设备包括发电机、热电联产机组、燃气锅炉、吸收式制冷机组、电锅炉、电制冷机和蓄电池,该多能互补园区与电力系统通过联络线相连,建立一个第一优化模型,第一优化模型用于计算多能互补园区向电力系统上报的联络线功率基线,作为电力系统日前调度的功率基线,第一优化模型以不加调控时多能互补园区自身总运行成本c0最小为目标,表达式如下:
其中,x为该第一优化模型中所有待求解变量构成的列向量,即:
为第iG台发电机在调度时刻t发出的有功功率;和分别为第iCHP台热电联产机组在调度时刻t的有功功率和供热功率;为第iGB台燃气锅炉在调度时刻t的供热功率;为第iAC台吸收式制冷机组在调度时刻t的供冷功率;和分别为第iEB台电锅炉在调度时刻t的耗电功率和供热功率;和分别为第iEC台电制冷机在调度时刻t的耗电功率和供冷功率;和分别为第iES台蓄电池在调度时刻t的充电功率和放电功率;为第iES台蓄电池在调度时刻t的充电状态的0-1变量,代表蓄电池在调度时刻t处在充电状态,代表蓄电池在调度时刻t未处在充电状态;为描述第iES台蓄电池在调度时刻t的放电状态的0-1变量,代表蓄电池在调度时刻t处在放电状态,代表蓄电池在调度时刻t未处在放电状态;为第iES台蓄电池在调度时刻t的充电转换状态的0-1变量,即代表蓄电池在调度时刻t-1未在充电、在调度时刻t处在充电状态,代表蓄电池在调度时刻t-1在充电、调度时刻t处在未充电状态;为第iES台蓄电池在调度时刻t的放电转换状态的0-1变量,即代表蓄电池在调度时刻t-1未在放电、调度时刻t处在放电状态,代表蓄电池在调度时刻t-1处在放电、调度时刻t未在放电状态;为调度时刻t的联络线功率,即多能互补园区与电力系统交换的电功率,以电功率流入多能互补园区为正方向;上标T为向量转置;γ为电力系统所有调度时刻t构成的集合;SG为多能互补园区内所有发电机构成的集合;为第iG台发电机在一个调度时刻发出单位有功功率所需成本,对于分布式光伏发电装置和分布式风电机组该值可取为0;SGB为所有燃气锅炉构成的集合;为第iGB台燃气锅炉在一个调度时刻内维持单位供热功率所需成本;SAC为所有吸收式制冷机构成的集合;为第iAC台吸收式制冷机组在一个调度时刻内维持单位供冷功率所需成本;SCHP为所有热电联产机组构成的集合;为第iCHP台热电联产机组在一个调度时刻内维持单位发电功率所需成本;为第iCHP台热电联产机组在一个调度时刻内维持单位供热功率所需成本;为调度时刻t的联络线电价;ΔT为相邻两个调度时刻的时间间隔;
第一优化模型的约束条件包括:
(1-1)多能互补园区中发电机的有功功率范围和爬坡约束:
(1-2)多能互补园区中热电联产机组的有功功率范围、供热功率范围和热电联产机组的有功功率爬坡约束:
(1-3)多能互补园区中燃气锅炉供热功率范围及爬坡约束:
(1-4)多能互补园区中吸收式制冷机的供冷功率范围及爬坡约束:
(1-5)多能互补园区中电锅炉的供热功率范围及爬坡约束:
(1-6)多能互补园区中电制冷机的供冷功率范围及爬坡约束:
(1-7)多能互补园区中蓄电池的运行约束:
其中,和为第iES台蓄电池的充电功率上限和下限;和为第iES台蓄电池的放电功率上限和下限;为一个调度日中第i台蓄电池的最大充放电转换次数;和为第iES台蓄电池中可储存的最大能量和最小能量;上述所有参数均在对应设备的说明书或操作规程中得到;为第iES台蓄电池在该调度日开始时储存的能量,由蓄电池前一日的调度结果得到;
(1-8)多能互补园区中的电能、热能、冷能的能量守恒约束:
(1-9)多能互补园区中冷、热负荷的惯性及室内温度约束:
其中,为热负荷的热容,和分别为热负荷t和t-1时刻的室内温度,为热负荷的热导,SHL为所有热负荷构成的集合,为t时刻的环境温度,和分别为热负荷室内温度最小值和最大值,为冷负荷的热容,和分别为冷负荷t和t-1时刻的室内温度,为冷负荷的热导,SCL为所有冷负荷构成的集合,和分别为冷负荷室内温度最小值和最大值。
采用分支定界法,求解由上述目标函数和约束条件组成的第一优化模型,得到多能互补园区与电力系统之间的联络线在调度时刻t的有功功率的最优值即多能互补园区联络线功率基线,以及目标函数最优值c0,即多能互补园区最小运行成本,多能互补园区将联络线功率基线上报电力系统;
(2)多能互补园区接收电力系统下发的调度时段t风电消纳需求为电力系统期望联络线在调度时段t的有功功率高于联络线功率基线的有功功率值;根据该风电消纳需求计算多能互补园区计算可提供的风电消纳能力,具体方法如下:
(2-1)选取计算步长ε,ε取值为0~1,置循环次数k初始值1;
(2-2)求解如下的第二优化模型,该第二优化模型用于在调度时段t,当联络线功率为时,求多能互补园区运行费用最小时多能互补园区中的能源设备的运行计划;第二优化模型的目标函数与第一优化模型的目标函数相同;
第二优化模型的约束条件包括:
a、多能互补园区中发电机的有功功率范围和爬坡约束:与第一优化模型中约束条件(1-1)相同;
b、多能互补园区中热电联产机组的有功功率和供热功率范围约束和热电联产机组的有功功率爬坡约束:与第一优化模型中约束条件(1-2)相同;
c、多能互补园区中燃气锅炉供热功率范围及爬坡约束:与第一优化模型中约束条件(1-3)相同;
d、多能互补园区中吸收式制冷机的供冷功率范围及爬坡约束:与第一优化模型中约束条件(1-4)相同;
e、多能互补园区中电锅炉的供热功率范围及爬坡约束:与第一优化模型中约束条件(1-5)相同;
f、多能互补园区中电制冷机的供冷功率范围及爬坡约束:与第一优化模型中约束条件(1-6)相同;
g、多能互补园区中蓄电池的运行约束:与第一优化模型中约束条件(1-7)相同;
h、多能互补园区中的电能、热能、冷能的能量守恒约束:
i、多能互补园区中冷、热负荷的惯性及室内温度约束:与第一优化模型中约束条件(1-9)相同;
采用分支定界法,求解由上述目标函数和约束条件组成的第二优化模型,若第二优化模型可解,则得到目标函数最优值ck,以及多能互补园区与电力系统之间的联络线在调度时刻t的有功功率的最优值即多能互补园区向电力系统上报的功率基线,转步骤(2-3);若该模型无解,则结束循环,进行步骤(3);
(2-3)计算风电消纳比例为kε时的多能互补园区的调节成本ck*,ck*=ck-c0,记录ck*,并置k=k+1,对kε进行判断,若kε>1,转步骤(3),若kε≤1,返回步骤(2-2),直到kε>1,得到多个与不同风电消纳比例kε相对应的多能互补园区的调节成本ck*;
(3)将步骤(2-3)的多个与不同风电消纳比例kε相对应的多能互补园区的调节成本ck*上报电力系统,电力系统下发最终采用的多能互补园区的风电消纳比例,实现由多能互补园区对电力系统风电消纳需求的响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911362865.1A CN111049184B (zh) | 2019-12-26 | 2019-12-26 | 一种多能互补园区响应电力系统风电消纳需求的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911362865.1A CN111049184B (zh) | 2019-12-26 | 2019-12-26 | 一种多能互补园区响应电力系统风电消纳需求的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111049184A true CN111049184A (zh) | 2020-04-21 |
CN111049184B CN111049184B (zh) | 2023-04-18 |
Family
ID=70239950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911362865.1A Active CN111049184B (zh) | 2019-12-26 | 2019-12-26 | 一种多能互补园区响应电力系统风电消纳需求的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111049184B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111861551A (zh) * | 2020-07-08 | 2020-10-30 | 海宁同兴智慧园区管理有限公司 | 一种园区用电智慧监测管理系统 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB912762A (en) * | 1960-07-12 | 1962-12-12 | Rolls Royce | Heat exchange apparatus |
US20160154397A1 (en) * | 2014-11-25 | 2016-06-02 | Acciona Energia, S.A. | Method for the control of power ramp-rates minimizing power storage requirements in intermittent power generation plants |
US20160273518A1 (en) * | 2014-09-02 | 2016-09-22 | Siemens Industry, Inc. | Systems, methods and apparatus for improved energy management systems with security constrained dynamic dispatch for wind power management |
CN106372742A (zh) * | 2016-08-19 | 2017-02-01 | 天津大学 | 考虑电转气多源储能型微网日前最优经济调度方法 |
WO2017107246A1 (zh) * | 2015-12-25 | 2017-06-29 | 北京四方继保自动化股份有限公司 | 一种含介质储能的智能微电网发电功率控制系统及控制方法 |
CN107341593A (zh) * | 2017-06-19 | 2017-11-10 | 东北电力大学 | 一种基于场景划分的电热综合系统弃风消纳协调调度模型 |
CN107634518A (zh) * | 2017-09-21 | 2018-01-26 | 国网福建省电力有限公司 | 一种“源‑网‑荷”相协调的主动配电网经济调度方法 |
CN108063451A (zh) * | 2017-11-17 | 2018-05-22 | 国网辽宁省电力有限公司 | 一种提高风电接纳能力的电网多元优化调度方法 |
CN109245183A (zh) * | 2018-05-21 | 2019-01-18 | 国网河南省电力公司安阳供电公司 | 一种基于负荷控制的高风光渗透率地区电网调峰方法 |
CN109409595A (zh) * | 2018-10-19 | 2019-03-01 | 南京南瑞继保电气有限公司 | 一种园区多能互补系统日前调度方法 |
CN109449993A (zh) * | 2018-12-26 | 2019-03-08 | 西安交通大学 | 一种面向多能互补电力系统的多时间尺度生产模拟方法 |
CN109638818A (zh) * | 2018-12-13 | 2019-04-16 | 中国电力科学研究院有限公司 | 一种储热式电采暖系统容量优化配置方法和装置 |
CN111049134A (zh) * | 2019-12-26 | 2020-04-21 | 清华大学 | 一种多能互补园区响应电力系统削峰需求的计算方法 |
-
2019
- 2019-12-26 CN CN201911362865.1A patent/CN111049184B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB912762A (en) * | 1960-07-12 | 1962-12-12 | Rolls Royce | Heat exchange apparatus |
US20160273518A1 (en) * | 2014-09-02 | 2016-09-22 | Siemens Industry, Inc. | Systems, methods and apparatus for improved energy management systems with security constrained dynamic dispatch for wind power management |
US20160154397A1 (en) * | 2014-11-25 | 2016-06-02 | Acciona Energia, S.A. | Method for the control of power ramp-rates minimizing power storage requirements in intermittent power generation plants |
WO2017107246A1 (zh) * | 2015-12-25 | 2017-06-29 | 北京四方继保自动化股份有限公司 | 一种含介质储能的智能微电网发电功率控制系统及控制方法 |
CN106372742A (zh) * | 2016-08-19 | 2017-02-01 | 天津大学 | 考虑电转气多源储能型微网日前最优经济调度方法 |
CN107341593A (zh) * | 2017-06-19 | 2017-11-10 | 东北电力大学 | 一种基于场景划分的电热综合系统弃风消纳协调调度模型 |
CN107634518A (zh) * | 2017-09-21 | 2018-01-26 | 国网福建省电力有限公司 | 一种“源‑网‑荷”相协调的主动配电网经济调度方法 |
CN108063451A (zh) * | 2017-11-17 | 2018-05-22 | 国网辽宁省电力有限公司 | 一种提高风电接纳能力的电网多元优化调度方法 |
CN109245183A (zh) * | 2018-05-21 | 2019-01-18 | 国网河南省电力公司安阳供电公司 | 一种基于负荷控制的高风光渗透率地区电网调峰方法 |
CN109409595A (zh) * | 2018-10-19 | 2019-03-01 | 南京南瑞继保电气有限公司 | 一种园区多能互补系统日前调度方法 |
CN109638818A (zh) * | 2018-12-13 | 2019-04-16 | 中国电力科学研究院有限公司 | 一种储热式电采暖系统容量优化配置方法和装置 |
CN109449993A (zh) * | 2018-12-26 | 2019-03-08 | 西安交通大学 | 一种面向多能互补电力系统的多时间尺度生产模拟方法 |
CN111049134A (zh) * | 2019-12-26 | 2020-04-21 | 清华大学 | 一种多能互补园区响应电力系统削峰需求的计算方法 |
Non-Patent Citations (3)
Title |
---|
BOMING ZHANG ET AL: "Real-Time Power Dispatch Method under Wind Power Curtailed Conditions for Large Scale Power Grid" * |
丁涛等: "考虑风电不确定性的区间经济调度模型及空间分支定界法" * |
陈沼宇;王丹;贾宏杰;王伟亮;郭炳庆;屈博;范孟华;: "考虑P2G多源储能型微网日前最优经济调度策略研究" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111861551A (zh) * | 2020-07-08 | 2020-10-30 | 海宁同兴智慧园区管理有限公司 | 一种园区用电智慧监测管理系统 |
CN111861551B (zh) * | 2020-07-08 | 2023-12-05 | 海宁同兴智慧园区管理有限公司 | 一种园区用电智慧监测管理系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111049184B (zh) | 2023-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107832979B (zh) | 一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 | |
CN110689189B (zh) | 考虑供能侧和需求侧的冷热电联合供需平衡优化调度方法 | |
US11095127B2 (en) | Method for real-time scheduling of multi-energy complementary micro-grids based on rollout algorithm | |
CN106950840B (zh) | 面向电网削峰的综合能源系统分层分布式协调控制方法 | |
CN106058942B (zh) | 考虑风电不确定性的含电转气和cchp的能量枢纽优化模型 | |
CN109004686A (zh) | 一种考虑冰蓄冷空调多模式的冷热电联供型微网系统 | |
CN109193815A (zh) | 一种提高风电消纳的热电联合调度方法 | |
CN110311421A (zh) | 基于需求侧响应的微电网多时间尺度能量管理方法 | |
CN111400641A (zh) | 一种含蓄热式电采暖的综合能源系统日前优化调度方法 | |
CN108521132B (zh) | 多能互补支撑电网频率调节的多时间尺度优化控制方法 | |
CN107294212B (zh) | 考虑不同空调负荷特性的微网双层优化调度方法及系统 | |
CN111313400B (zh) | 一种基于鲁棒修正的多能虚拟电厂运行参数聚合方法 | |
CN107749645B (zh) | 一种控制高压大容量储热供暖装置的方法 | |
CN110807588A (zh) | 一种多能源耦合综合能源系统的优化调度方法 | |
CN107957681A (zh) | 一种研究冷热电联供微网环保经济调度的方法 | |
CN113240166A (zh) | 考虑高比例新能源消纳的微能源网日前经济调度方法 | |
CN115173470A (zh) | 一种基于电网调峰的综合能源系统调度方法及系统 | |
CN116308881A (zh) | 利用供热管网储热的综合能源系统多时间尺度调度方法 | |
CN108199369A (zh) | 面向多种电力市场准入关联的地域性微电网低成本优化运营方法 | |
CN111049184B (zh) | 一种多能互补园区响应电力系统风电消纳需求的计算方法 | |
CN204407890U (zh) | 一种可再生能源冷热电微网系统 | |
CN111049134B (zh) | 一种多能互补园区响应电力系统削峰需求的计算方法 | |
CN113537618A (zh) | 一种考虑居民用户需求响应的综合能源系统优化调度方法 | |
CN112072645A (zh) | 一种日前聚合与日内滚动修正衔接的虚拟电厂等值方法 | |
CN109995030B (zh) | 一种考虑脱网风险的储能装置soc下限值最优设置方法 |
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 |