CN104392286A - 考虑冷热电联供和储能运行策略的微电网运行优化方法 - Google Patents
考虑冷热电联供和储能运行策略的微电网运行优化方法 Download PDFInfo
- Publication number
- CN104392286A CN104392286A CN201410724835.1A CN201410724835A CN104392286A CN 104392286 A CN104392286 A CN 104392286A CN 201410724835 A CN201410724835 A CN 201410724835A CN 104392286 A CN104392286 A CN 104392286A
- Authority
- CN
- China
- Prior art keywords
- power
- energy
- micro
- capacitance sensor
- fuel cell
- 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
- 238000010438 heat treatment Methods 0.000 title claims abstract description 52
- 238000001816 cooling Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000013486 operation strategy Methods 0.000 title claims abstract description 18
- 238000004146 energy storage Methods 0.000 claims abstract description 174
- 239000000446 fuel Substances 0.000 claims abstract description 125
- 230000005611 electricity Effects 0.000 claims abstract description 97
- 238000005457 optimization Methods 0.000 claims abstract description 26
- 239000007789 gas Substances 0.000 claims description 77
- 238000007599 discharging Methods 0.000 claims description 27
- 238000010248 power generation Methods 0.000 claims description 25
- 238000005057 refrigeration Methods 0.000 claims description 25
- 238000003860 storage Methods 0.000 claims description 16
- WKBOTKDWSSQWDR-UHFFFAOYSA-N Bromine atom Chemical compound [Br] WKBOTKDWSSQWDR-UHFFFAOYSA-N 0.000 claims description 12
- GDTBXPJZTBHREO-UHFFFAOYSA-N bromine Substances BrBr GDTBXPJZTBHREO-UHFFFAOYSA-N 0.000 claims description 12
- 229910052794 bromium Inorganic materials 0.000 claims description 12
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 238000004378 air conditioning Methods 0.000 claims description 9
- 238000012423 maintenance Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 239000011435 rock Substances 0.000 claims description 9
- 238000007667 floating Methods 0.000 claims description 8
- 238000002922 simulated annealing Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 229910052799 carbon Inorganic materials 0.000 claims description 5
- 238000007710 freezing Methods 0.000 claims description 5
- 230000008014 freezing Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000002485 combustion reaction Methods 0.000 claims description 4
- 238000007620 mathematical function Methods 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 4
- 239000002918 waste heat Substances 0.000 claims description 4
- 235000019504 cigarettes Nutrition 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 15
- 230000007613 environmental effect Effects 0.000 abstract description 6
- 239000000809 air pollutant Substances 0.000 abstract description 2
- 231100001243 air pollutant Toxicity 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 210000004027 cell Anatomy 0.000 description 73
- 230000006870 function Effects 0.000 description 19
- 239000002245 particle Substances 0.000 description 15
- 238000010586 diagram Methods 0.000 description 7
- 230000009467 reduction Effects 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- AMXOYNBUYSYVKV-UHFFFAOYSA-M lithium bromide Chemical compound [Li+].[Br-] AMXOYNBUYSYVKV-UHFFFAOYSA-M 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000003344 environmental pollutant Substances 0.000 description 3
- 231100000719 pollutant Toxicity 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 229940059936 lithium bromide Drugs 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 230000005619 thermoelectricity Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000003850 cellular structure Anatomy 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 244000144992 flock Species 0.000 description 1
- 239000003546 flue gas Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Fuel Cell (AREA)
Abstract
本发明公开了一种考虑冷热电联供和储能运行策略的微电网运行优化方法,包括:建立CCHP型微电网优化模型;构建以微电网并网运行状态下发电成本最低为目标的目标优化函数;根据建立的目标优化函数,综合考虑前后时段的负荷水平以及各微电源的出力趋势,确定储能单元以及燃料电池的运行原则,进而确定调度周期内整个CCHP系统的最佳运行方式。本发明有益效果:对储能优化能够有效地增加微电网售电收益,对降低微电网的运行成本具有明显优势;相对于传统的火电厂、热电厂能够提高系统能源利用率并实现能源梯级利用,为电力系统节能提供参考;能够节能减排,极大的减少空气污染物的排放,对于提高环境效益具有重大意义。
Description
技术领域
本发明涉及微电网运行优化领域,具体涉及一种考虑冷热电联供和储能运行策略的微电网运行优化方法。
背景技术
冷热电联供(combined cooling heating and power,CCHP)系统将制冷、制热和供电三者相结合,能够实现能量梯级利用并减少污染物排放,是一种节能、环保的供能方式,具有良好的社会和经济效益。而微电网则解决了分布式电源大规模接入电网的问题,能够有效、灵活、智能地利用各种分布式能源,在缓解能源短缺现状、保护环境和提高电能质量方面具有巨大潜力。在发展低碳绿色经济的背景下,两者结合的CCHP型微电网的发展受到了极大地重视,但目前技术还不成熟,应用尚处于起步阶段。
CCHP型微电网存在冷、热、电三种能量之间的平衡关系,而微电网具有并网和孤岛运行两种方式。基于能源利用率和环保效益的考虑,在满足系统冷热负荷的前提下,如何根据各微电源特性制定最佳的运行策略以提高微电网的经济效益这一问题受到越来越多的关注。
目前,关于CCHP型微电网的研究多以优化电源配置和储能容量、最佳备用储能等为目标,主要包括以下方法:
1.建立热电联产型微电网,在模糊优化理论的基础上,运用改进遗传算法优化并网运行方式下各微源的有功、无功出力和多目标优化的满意度。
2.应用机会约束规划理论建立了经济运行优化模型,并采用基于随机模拟技术的粒子群优化算法求解模型,根据不同的微源配置,对系统的运行方案进行了优化。
3.采用外推近似法对冷电联供系统进行优化求解,以日综合运行经济效益最高为目标,在满足给定冷负荷的基础上,确定优化的运行模式,并对蓄冷器的使用进行了研究。
以上研究对微电网的经济运行和储能容量的优化配置等方面做了大量的工作,但在提高微电网能源利用率以及优化储能策略的方面上仍需要进一步研究。
发明内容
本发明的目的就是为了解决上述问题,提出了一种考虑冷热电联供和储能运行策略的微电网运行优化方法,该方法以能源梯级利用为基础,以一个包含光伏(photovoltaic,PV)、风机(wind turbine,WT)、微型燃气轮机(micro turbine,MT)、空调机(air conditioner,AC)、燃料电池(fuel cell,FC)、储能单元(energy storage,ES)的CCHP型微电网为研究对象,在满足电能供需平衡、储能容量限制等约束下,确定最佳的储能运行策略并采用基于模拟退火的粒子群(SimuAPSO)动态规划法求得调度周期内各微电源的最佳运行方式,并通过能源利用率、污染物减排量及减排费用和运行成本对比等过程评价指标对结果进行量化分析。
为了实现上述目的,本发明采用如下技术方案:
一种考虑冷热电联供和储能运行策略的微电网运行优化方法,包括以下步骤:
(1)假设在调度周期单位时间间隔Δt内各微电源的出力恒定并满足冷、热、电负荷的平衡关系,微电网与主网之间功率交互功率恒定,交互电价跟随实时电价并在Δt内保持不变的情况下,建立CCHP型微电网优化模型;
(2)根据功率平衡约束条件、分布式电源有功出力约束条件、空调机出力约束条件以及储能单元运行约束条件,构建以微电网并网运行状态下发电成本最低为目标的目标优化函数;
(3)根据建立的目标优化函数,利用基于模拟退火粒子群的动态规划法求解输出储能单元、交互功率和燃料电池的最佳出力;综合考虑前后时段的负荷水平以及各微电源的出力趋势,确定储能单元以及燃料电池的运行原则,进而确定调度周期内整个CCHP系统的最佳运行方式。
所述步骤(1)中建立的CCHP型微电网优化模型包括:微型燃气轮机模型、空调机模型、燃料电池模型、风力发电模型、光伏发电模型和储能单元模型;
假设微型燃气轮机排烟温度与溴冷机的进烟温度保持不变,并忽略外界环境的变化对微型燃气轮机供电、燃料燃烧效率的影响,所述微型燃气轮机模型的数学模型为:
QMT-h(t)=QMT(t)ηhCOPh (2)
QMT-c(t)=QMT(t)ηcCOPc (3)
微型燃气轮机的燃料成本为:
式中,令t=t0+rΔt,t0为调度初始时刻,r为非负整数,调度周期T=NTΔt,NT为正整数;QMT(t)为燃气轮机排气余热量;PMT(t)为t时刻燃气轮机输出的电功率;ηMT为燃气轮机发电效率;ηL为散热损失系数;QMT-h(t)、QMT-c(t)为溴冷机t时刻提供的制热量和制冷量;ηh、ηc和COPh、COPc分别为溴冷机的制热和制冷时的烟气回收率及制热和制冷系数;;微型燃气轮机不同功率下的发电效率由生产商提供的相关参数拟合得到;CMT为微型燃气轮机的燃料成本;Cf为单位燃料的价格;LHVNG为天然气低热热值;
所述空调机模型的数学函数为:
Qair-h(t)=Pair(t)(1-ηLh)ηahCOPair-h (5)
Qair-c(t)=Pair(t)(1-ηLc)ηacCOPair-c (6)
式中,Qair-h(t)、Qair-c(t)分别为t时刻空调机制热、制冷功率;Pair(t)为空调机用电功率;ηLh、ηLc分别为空调机制热、制冷时损失系数,ηah、ηac分别为空调机的制热、制冷效率;COPair-h、COPair-c分别为空调制热和制冷时的能效比;空调机所耗电能作为系统用电负荷统一调度。
所述燃料电池模型的数学函数为:
不考虑燃料电池的热能利用,燃料电池的燃料消耗-功率输出特性为:
式中,CFC为燃料电池的燃料成本;PFC(t)为t时刻燃料电池发电功率;ηFC为燃料电池发电效率;LHVNG为天然气低热热值,Cf为单位燃料价格、r为调度时间段、NT为总的调度时间段数目,;燃料电池不同功率下的发电效率由生产商提供的相关参数拟合得到;
所述储能单元模型的状态有充电、放电、浮充三种状态,储能在t时刻的剩余电量与t-1时刻的剩余电量、t-1到t时刻的充放电量以及自放电量有关;储能的输入、输出功率及充放电状态满足下列关系:
式中,Csoc(t)为t时刻储能剩余电量;τ为储能自放电率;Pes为储能t时刻的充放电功率;ηch、ηdis分别为储能充电和放电效率;ufl(t)、uch(t)、udis(t)∈﹛0,1﹜,为储能单元的状态标识,分别代表浮充,充电,放电三种状态,对应Pes(t)取值分别等于、小于、大于零,Δt为单位调度时间间隔。
所述步骤(2)中构建以微电网并网运行状态下发电成本最低为目标的目标优化函数具体为:
minF=F1-F2 (9)
式中,F为微电网的调度周期内总发电成本;F1为微电源运行成本函数;F2为微电网经济收益函数;
微电源运行成本函数为:
F1=Ct+Cfu+Cm (10)
式中,Ct、Cfu、Cm分别代表微电网折旧成本、燃料成本和维护成本;N为微电源总数;Ci为第i台微电源的初始成本;l为利息率;ki为第i台微电源的寿命周期;CBAT为储能初始成本;v为储能寿命周期,微电源与储能折旧的平摊成本为期初成本;Cmi为第i台微电源的单位出力维护成本;Pi(t)为第i台微电源单位时间内有功出力;Cmes为储能的单位出力维护成本,CMT为微型燃气轮机的燃料成本、CFC为燃料电池的成本、NT为调度周期间隔数、Pes(t)为储能t时刻的充放电功率。
所述微电网经济收益函数为:
式中,Prs(t)、Prb(t)分别为t时刻的售电和购电电价;Ps(t)、Pb(t)分别为t时刻的售电和购电功率;Chc为制热(冷)收益;Che、Cco分别为单位制热和制冷量的售价;a=1表示冬季制热,a=0表示夏季制冷;QMT-h(t)为微型燃气轮机的产热量、QMT-c(t)为微型燃气轮机的制冷量。
所述步骤(2)中的功率平衡约束条件具体为:
式中,Pload(t)、Qhe(t)、Qco(t)分别为t时刻网内普通电负荷和热、冷负荷;Ploss(t)为t时刻微电网功率损耗,QMT-h(t)为微型燃气轮机的产热量、QMT-c(t)为微型燃气轮机的制冷量,PMT(t)为微型燃气轮机的发电功率,Pair(t)为空调机用电功率;Pb(t)、Ps(t)分别为微电网从大电网的购电功率和售电功率。,Qair-h(t)、Qair-c(t)分别为t时刻空调机制热、制冷功率。Pi(t)为第i个微电源的出力功率。
所述步骤(2)中的分布式电源有功出力约束条件具体为:
Pimin≤Pi(t)≤Pimax (17)
式中,Pimax,Pimin分别为第i台微电源有功出力的上下限;
所述空调机出力约束条件具体为:
Pairmin≤Pair(t)≤Pairmax (18)
式中,Pairmin、Pairmax为空调机出力上下限。
所述步骤(2)中的储能单元运行约束条件具体为:
①储能单元的出力功率约束条件:
Pesmin≤Pes(t)≤Pesmax (19)
式中:Pesmax,Pesmin分别为储能出力上下限;
②储能单元的能量约束条件:
Csocmin≤Csoc(t)≤Csocmax (20)
式中:Csocmax,Csocmin分别为储能容量的最大、最小值;其中,SOC表示储能的核电状态,表征其能量范围;
③单位调度周期始末状态储能能量约束条件:
Csoc(t0)=Csoc(t0+NrΔt) (21)
考虑储能单元对微电网的周期性优化,其能量状态满足在调度周期始末相等;Csoc(t)为t时刻储能剩余电量;Δt为单位调度时间间隔,NT为总的调度时间段数目;
单位调度周期储能的充放电状态和次数约束:
uch(t)+udis(t)+ufl(t)=1 (22)
式中:λ1、λ2分别表示储能的充放电次数,具体的数值由系统负荷、储能寿命和优化策略来确定;uch(t)、udis(t)、ufl(t)分别为代表t时刻储能的充电、放电、浮充三种状态,Δt为单位调度时间间隔。
所述步骤(3)中,调度周期内整个CCHP系统的最佳运行方式的确定原则为:
空调机与微型燃气轮机单独使用;
定义Pe(t)为t时刻储能和燃料电池均不启动时微电网广义负荷功率,正值表示微电网输出电能,负值表示微电网缺电并假设大电网能完全消纳与弥补微电网功率差额;
依据Pe(t)与分时电价制定的储能单元交换功率和燃料电池的运行原则为:
1)若微电网售电价格较低,储能单元优先充电:若微电网电量剩余,则向大电网售出电能,燃料电池不工作;若微电网电量不足,比较购电与燃料电池发电成本,选择成本较低者;
2)若微电网售电价格较高,储能单元优先放电:若微电网电量剩余,则向大电网售出电能,燃料电池不工作;若微电网电量不足,比较购电与燃料电池发电成本,选择价格较低者;
3)若购电成本低于燃料电池发电成本,交互功率为储能单元充放电后微电网功率缺额;
4)若燃料电池发电成本低于购电成本,其发电功率为储能充电之前或放电之后系统功率缺额。
定义仅考虑储能单元交换功率和燃料电池运行的期望成本函数如下:
式中,Fin_t1(t)、Fin_t2(t)分别表示燃料电池不启动下售电与购电的成本函数;
Ptex(t)=Pe(t)+Pes(t),表示在燃料电池不启动时系统广义交互功率,正值表示微电网售电功率,负值表示从大电网购电功率;Pe(t)、为t时刻为t时刻储能和燃料电池均不启动时微电网广义负荷功率,Pes(t)为储能的充放电功率;
式中,Prsmax为调度周期内最高售电电价;Fa(t)、Fc(t)、Fe(t)为各情况下非峰时放电期望损失;Fb(t)、Fd(t)、Ff(t)分别为各情况下非峰时充电期望收益;Ptex(t)表示在燃料电池不启动时系统广义交互功率、Prs(t)、Prb(t)分别为t时刻的售电和购电电价、Δt为单位调度时间间隔、CFC(t)为燃料电池的发电成本、ηch和ηdis分别为储能的充、放电功率;
燃料电池的出力为
其中:当PFC(t)=-Pe(t)时,Pes(t)=0;
微电网与大电网的交互功率为
由式(26)、(27)所示,在计算电能交互收益时,在广义负荷功率为正或负的情况下,对储能单元充放电过程分别考虑了充电的期望收益与放电过程中期望损失,即若非峰时电价放电时,售电收益减去的非峰时放电期望损失,若非峰时充电时,售电收益减去非峰时充电期望收益;
其中,Fa(t)、Fb(t)、Ff(t)表示系统功率过剩时,储能系统充放电与替代供电方式之间的期望损失或收益,而式Fc(t)、Fd(t)、Fe(t)表示系统功率不足时,储能系统充放电与替代供电方式之间的期望损失或收益;
因此,储能单元在动态寻优的任一时刻,将尽量避开非峰时放电与非谷时充电状态,并依据微电网前后时刻广义负荷功率Pe(t)、储能单元的剩余电量、实时电价和燃料电池发电成本决策最佳充放电出力,以保证系统期望收益函数取最大,满足运行原则1)和2);
由式(29)、(30)表示,当系统广义交互功率Ptex(t)为负时,比较从大电网购电与燃料电池成本,从而确定是否启动燃料电池及其出力大小,同时由式(29)可知,燃料电池仅满足系统的功率缺额,不向储能单元充电及大电网售电,满足运行原则3)和4);
储能单元以及燃料电池的运行原则确定后,调度周期内整个CCHP型微电网的最佳运行方式随之确定。
本发明的有益效果是:
本发明基于小型微电网建立了考虑冷热电联产和储能优化策略并计及制冷与制热收益的经济调度模型,在能源梯级利用和储能优化策略的基础上,采用基于模拟退火粒子群的动态规划法进行求解,分析对比了多种策略,同时采用了量化指标对结果进行评估,得出以下结论:
1)储能系统采取文中策略能够有效地增加微电网售电收益,对降低微电网的运行成本具有明显优势。
2)CCHP系统相对于传统的火电厂、热电厂能够提高系统能源利用率并实现能源梯级利用,为电力系统节能提供参考。
3)CCHP型微电网能够节能减排,极大的减少空气污染物的排放,对于提高环境效益具有重要参考意义。
附图说明
图1为本发明实施例CCHP型微电网结构示意图;
图2为本发明实施例冬季典型日热电负荷需求示意图;
图3为本发明实施例夏季典型日冷电负荷需求示意图;
图4为本发明实施例冬季典型日三种策略下微网与主网交互功率示意图;
图5为本发明实施例冬季典型日三种策略下燃料电池的出力示意图;
图6为本发明实施例策略2和3中储能单元的出力示意图;
图7为本发明实施例冬季典型日策略3的电负荷平衡情况示意图;
图8为本发明实施例夏季典型日策略3的电负荷平衡情况示意图。
具体实施方式:
下面结合附图与实施例对本发明做进一步说明:
一种考虑冷热电联供和储能运行策略的微电网运行优化方法,包括以下步骤:
(1)假设在调度周期单位时间间隔Δt内各微电源的出力恒定并满足冷、热、电负荷的平衡关系,微电网与主网之间功率交互功率恒定,交互电价跟随实时电价并在Δt内保持不变的情况下,建立CCHP型微电网优化模型;
(2)根据功率平衡约束条件、分布式电源有功出力约束条件、空调机出力约束条件以及储能单元运行约束条件,构建以微电网并网运行状态下发电成本最低为目标的目标优化函数;并对CCHP系统能源利用率指标和节能减排指标进行分析;
(3)根据建立的目标优化函数,综合考虑前后时段的负荷水平以及各微电源的出力趋势,确定储能单元以及燃料电池的运行原则,进而确定调度周期内整个CCHP系统的最佳运行方式。
(4)针对实际的微电网系统进行算例仿真,分析所提模型和策略的合理性,并对优化结果进行CCHP系统能源利用率指标和节能减排指标分析。
1、CCHP型微电网优化模型
假设在调度周期单位时间间隔Δt内各微电源的出力恒定并满足冷、热、电负荷的平衡关系,微电网与主网之间功率交互功率恒定,交互电价跟随实时电价并在Δt内保持不变。
CCHP型微电网优化模型包括:微型燃气轮机模型、空调机模型、燃料电池模型、风力发电模型、光伏发电模型和储能单元模型;
1.1微型燃气轮机模型
CCHP系统核心装置为微型燃气轮机,微型燃气轮机发电同时排出的高温烟气通过溴化锂吸收式制冷机用于制冷、取暖和供应生活热水等。假设微型燃气轮机排烟温度与溴冷机的进烟温度保持不变,并忽略外界环境的变化对微型燃气轮机供电、燃料燃烧效率等的影响,溴化锂吸收式制冷机采用双效或多效制冷机。含微型燃气轮机的CCHP系统的数学模型为
QMT-h(t)=QMT(t)ηhCOPh (32)
QMT-c(t)=QMT(t)ηcCOPc (33)
式中:令t=t0+rΔt,t0为调度初始时刻,r为非负整数,调度周期T=NTΔt,NT为正整数;QMT(t)为燃气轮机排气余热量;PMT(t)为t时刻燃气轮机输出的电功率;ηMT为燃气轮机发电效率;ηL为散热损失系数;QMT-h(t)、QMT-c(t)为溴冷机t时刻提供的制热量和制冷量;ηh、ηc和COPh、COPc分别为溴冷机的制热和制冷时的烟气回收率及制热和制冷系数;VMT为调度周期内消耗的天然气量;微型燃气轮机不同功率下的发电效率可由生产商提供的相关参数拟合得到。
微型燃气轮机的燃料成本为
式中:CMT为微型燃气轮机的燃料成本;Cf为天然气的价格,取2.54元/m3;LHVNG为天然气低热热值,取9.7kWh/m3。
1.2空调机模型
关于空调出力,当电压高于额定值时,空调出力随着电压下降较快,而低于额定电压时,空调出力下降比较缓慢,因此文中假设空调负荷在制冷或制热时为恒功率出力。
空调机出力的数学模型为
Qair-h(t)=Pair(t)(1-ηLh)ηahCOPair-h (35)
Qair-c(t)=Pair(t)(1-ηLc)ηacCOPair-c (36)
式中:Qair-h(t)、Qair-c(t)分别为t时刻空调机制热、制冷功率;Pair(t)为空调机用电功率;ηLh、ηLc分别为空调机制热、制冷时损失系数,ηah、ηac分别为空调机的制热、制冷效率;COPair-h、COPair-c分别为空调制热和制冷时的能效比。空调机所耗电能作为系统用电负荷统一调度。
1.3燃料电池模型
由于燃料电池发电效率较高,且作为微电网主要电能调度型电源,因此不考虑燃料电池的热能利用,燃料电池的燃料消耗-功率输出特性为
式中:CFC为燃料电池的燃料成本;PFC(t)为t时刻燃料电池发电功率;ηFC为燃料电池发电效率;燃料电池不同功率下的发电效率可由生产商提供的相关参数拟合得到。
1.4储能单元模型
储能单元的状态有充电、放电、浮充三种状态,储能在t时刻的剩余电量与t-1时刻的剩余电量、t-1到t时刻的充放电量以及自放电量有关。储能的输入、输出功率及充放电状态应满足下列关系:
式中:Csoc(t)为t时刻储能剩余电量;τ为储能自放电率,取3.6%/月;Pes为储能t时刻的充放电功率;ηch、ηdis分别为储能充电和放电效率;ufl(t)、uch(t)、udis(t)∈﹛0,1﹜,为储能单元的状态标识,分别代表浮充,充电,放电三种状态,对应Pes(t)取值分别等于、小于、大于零。
风力发电模型、光伏发电模型
1.5风力发电模型
风力发电的输出与风速关系成非线性,理想的功率曲线数学表达式为:
式中:Pwt为风电出力;v、vci、vco、vrate分别为风机的当前风速、切入风速、切除风速和额定风速;Prate为额定功率;λ1、λ2、λ3分别为风机的出力系数,可以通过曲线拟合得到。
1.6光伏电池模型
对于玻璃-玻璃封装的光伏电池组件,由如下经验计算式,可以计算环境温度从而估算出组件温度:
Tmod=Tamd+30G/1000
式中:Tmod、Tamd分别为组件温度和环境温度;G为组件接受到的太阳能辐射值。
光伏电池的输出功率计算式为:
PPV=PSTCGAC[1+k(Tc-Tt)]/GSTC
式中:PPV、PSTC分别为当前和标准测试条件下(环境温度为25℃)光伏电池的出力;GAC、GSTC分别为当前和标准测试下的光照强度,其中标准测试下光照强度取1000W/m2;k为功率温度系数,其值取-0.47%/K,Tc、Tt分别为电池板的工作温度和参考温度。
风机出力和光伏电池年出力、月出力、日出力等变化较大,为更好地描述其出力变化特性,采用时序模拟仿真法在给定负荷条件和预测风光出力的情况下将风电光伏出力等看做随时间变化的时间序列,并采用最大功率跟踪模式。
2、CCHP型微电网运行优化
CCHP型微电网主要是在能源梯级利用基础上,以最低运行成本实现微电网的运行效益最优化。因此,微电网中冷热负荷主要由CCHP系统供应,通过跟踪系统热(冷)负荷确定微型燃气轮机的出力,对于不满足微型燃气轮机出力的情况下启用空调,并对冷热电联产系统进行能源利用率、CO2减排量及单位节能减排成本进行评估。对网内的电负荷供应,考虑环境效益及调度可控性等因素,优先使用可再生清洁能源并制定合理的储能单元和燃料电池的运行策略,从而实现对负荷的削峰填谷基础上充分利用分时电价的差异以使得微电网更好的参与到电力市场中去。
2.1经济效益分析
微电网经济优化的主要内容为并网运行状态下发电成本最低。忽略突发事故导致的损失和新能源发电补贴等经济因素,微电网发电成本主要有运行成本和经济收益两部分,假设分布式电源发电仅与有功出力相关,其发电成本函数为:
minF=F1-F2 (39)
式中:F为微电网的调度周内总发电成本;F1为微电源运行成本;F2为微电网经济收益。
微电源运行成本函数为:
F1=Ct+Cfu+Cm (40)
式(43)、(43)中:Ct、Cfu、Cm分别代表微电网折旧成本、燃料成本和维护成本;N为微电源总数;Ci为第i台微电源的初始成本;l为利息率;ki为第i台微电源的寿命周期;CBAT为储能初始成本;v为储能寿命周期,微电源与储能折旧的平摊成本为期初成本;Cmi为第i台微电源的单位出力维护成本;Pi(t)为第i台微电源单位时间内有功出力;Cmes为储能的单位出力维护成本。
微电网经济收益为:
式(45)、(46)中:Prs(t)、Prb(t)分别为t时刻的售电和购电电价;Ps(t)、Pb(t)分别为t时刻的售电和购电功率;Chc为制热(冷)收益;Che、Cco分别为单位制热和制冷量的售价;本文冬季不考虑冷负荷且夏季不考虑热负荷,因此,a=1表示冬季制热,a=0表示夏季制冷。
2.2 CCHP系统节能减排分析
1)能源利用率指标。为评价CCHP系统的能源利用效率,利用一次能源利用效率指标(primary energy ratio,PER)VPER对系统的运行状况进行评估,一般联产系统要求VPER≥75%,并与传统分产系统能源利用率进行比较;其中一次能源利用率为有效利用能量与总耗能之比:
式中:ΣQf(t)为调度周期内为联供系统消耗天然气量,且ΣQf(t)=VMT×LHVNG,LHVNG为天然气低热热值,PMT(t)为微型燃气轮机的发电功率、Qhe(t)为系统的热负荷、Qco(t)为系统的冷负荷。
2)节能减排指标。由于微型燃气轮机的主要排放物为CO2、SO2、NOx等,因此评价CCHP系统的节能减排效益,以相对传统的分产火电机组的日平均排放物减排率、减排节约成本等为指标,对其减排运行状况进行评估。
式中,j表示排放物的数目,由于仅考虑CO2,SO2,NOx三种排放物,因此MT=3;ηdj、Cdcj为CCHP系统污染物减排率、减排节约成本;Fcoalj、FMTj分别为分产火电机组和微型燃气轮机的污染物排放因子;Peq(t)为产生与CCHP系统同等冷热电能的等效火电机组出力;ξeh、ξec分别为热负荷和冷负荷与火力机组出力的转化系数,具体数值由转换机组的效率决定;ζgasj为排放物惩罚因子;QMT-h(t)为微型燃气轮机的产热量、QMT-c(t)为微型燃气轮机的制冷量,PMT(t)为微型燃气轮机的发电功率。
2.3约束条件
若忽略燃气轮机与溴冷机之间的热能损失,微电网在运行期间主要满足系统冷热电功率平衡与分布式电源、空调机及储能单元运行约束。
1)功率平衡约束
式(50)-(52)中:Pload(t)、Qhe(t)、Qco(t)分别为t时刻网内普通电负荷和热、冷负荷;Ploss(t)为t时刻微电网功率损耗。
2)分布式电源有功出力约束
Pimin≤Pi(t)≤Pimax (50)
式中:Pimax,Pimin分别为第i台微电源有功出力的上下限。
3)空调机出力约束
Pairmin≤Pair(t)≤Pairmax (51)
式中:Pairmin、Pairmax为空调机出力上下限;
4)储能单元运行约束
④储能单元的出力功率约束
Pesmin≤Pes(t)≤Pesmax (52)
式中:Pesmax,Pesmin分别为储能出力上下限;
⑤储能单元的能量约束
Csocmin≤Csoc(t)≤Csocmax (53)
式中:Csocmax,Csocmin分别为储能容量的最大、最小值。其中,SOC表示储能的核电状态,表征其能量范围,文中取值为0.2-0.8;
⑥单位调度周期始末状态储能能量约束:
Csoc(t0)=Csoc(t0+NrΔt) (54)考虑储能单元对微电网的周期性优化,其能量状态满足在调度周期始末相等;
⑦单位调度周期储能的充放电状态和次数约束:
uch(t)+udis(t)+ufl(t)=1(55)
式中:λ1、λ2分别表示储能的充放电次数,具体的数值由系统负荷、储能寿命和优化策略来确定。
3、储能系统最佳运行策略
CCHP系统中微型燃气轮机出力低于最大出力50%时,发电效率大幅度降低。对于不满足微型燃气轮机出力限制的热(冷)负荷采用空调机补燃。虽然空调机同溴冷机一样为风机盘管设计,但是2个系统的水泵压力不同,可能出现一方无法将冷水压入管道的情况。此外,将联供系统的冷却水与空调机冷却水混合,不易分配两者之间合理的制热(冷)量。因此,空调机与微型燃气轮机单独使用。
综上,CCHP系统与可再生能源等运行策略已定,根据系统各时段的冷、热电负荷的需求,可确定储能单元及燃料电池的运行策略。由于研究重点为在满足系统负荷需求下采用合理策略调度储能,以充分利用峰谷电价达使微电网运行成本最低,因此在满足储能充放电次数约束下忽略储能非线性充放电损耗等经济损失。储能系统的优化属于多阶段非线性动态规划问题,需要对储能的运行在调度周期内进行动态规划处理。
粒子群优化算法(PSO)模拟了鸟群飞行觅食的行为,通过鸟之间的集体协作使群体达到最优。PSO求解优化问题时,每个粒子都有自己的位置和速度,还有一个由目标函数决定的适应值。粒子通过跟踪自身最好和全局最好值来寻优。该算法概念简明、实现方便、参数设置少,是一种高效的搜索算法,已经广泛应用在电力系统的求解无功优化、动态安全边界辨识、配电系统状态估计等方面。但是,由于粒子的移动没有选择性,粒子容易跳出最优解附近的某一邻域,使其局部搜索能力较弱,易收敛为局部最优,因此将其与模拟退火算法相结合,使粒子的飞行无记忆性并结合模拟退火算法重新生成停止进化粒子的位置,增强了全局搜索能力。此外SimuAPSO算法采用双群体搜索机制,一个群体保存具有可行解的粒子,使其逐步搜索到最优可行解;另一群体保存具有不可行解的粒子,且可行解群体以一定概率接受具有不可行解的粒子,有效地维持了群体的多样性。具体解题步骤如下:
Step1:输入原始数据(风速、温度、冷热电负荷等)并确定Pwt(t)、Ppv(t)、Pmt(t)和Pair(t),同时确定Pe(t)。
Step2:置粒子群维数K,粒子个数M,最大迭代次数N,并根据约束条件初始化粒子群位置与速度。
Step3:根据适应度函数f评价每个微粒的适应度,将记录各微子历史最佳位置pi和群体最优位置pg。
Step4:设置初始温度t0=f(pg)/ln5,根据式(下)确定当前温度下各pi的适配值TF(pi)。
Step5:根据适配值TF(pi)并采用轮盘赌策略从所有pi中确定全局最优的某个替代值pg′,然后根据Step1、Step2、Step3更新各微粒的速度和位置:
xi,j(t+1)=xi,j(t)+vi,j(t+1),j=1,2,...K (63)
C=c1+c2 (65)
式中:c1、c2为粒子权重系数;为加速因子;r1、r2为(0,1)之内均匀分布的随机数;xij(t)、vij(t)分别为粒子t时刻的位置和速度。
Step6:计算各微粒新的适应值,更新各微粒的pi值及群体的pg值,并按式进行退温操作。
tn=λtn-1,n=1,2……N (66)
式中:λ为退温常数。
Step7:若满足最大迭代次数,搜索停止,输出储能单元、交互功率和燃料电池的最佳出力,否则返回Step4继续搜索。
在调度周期的NT个时段对储能单元进行动态规划,以储能单元剩余电量作为各阶段的决策变量,为模拟储能单元充放电过程,首先将储能单元出力变化的连续过程离散化,同时为让储能单元更好的参与到经济调度,必须要综合考虑前后时段的负荷水平以及各微电源的出力趋势,在满足系统的负荷需求下对储能动态充放电。而考虑环境、能源利用率等因素,燃料电池出力应该只满足系统缺电功率,不应向储能充电及大电网售电。
定义Pe(t)为t时刻储能和燃料电池均不启动时微电网广义负荷功率,正值表示微电网输出电能,负值表示微电网缺电并假设大电网能完全消纳与弥补微电网功率差额。依据Pe(t)与分时电价制定的储能单元交换功率和燃料电池的运行原则为:
1)若微电网售电价格较低,储能单元优先充电:若微电网电量剩余,则向大电网售出电能,燃料电池不工作;若微电网电量不足,比较购电与燃料电池发电成本,选择成本较低者。
2)若微电网售电价格较高,储能单元优先放电:若微电网电量剩余,则向大电网售出电能,燃料电池不工作;若微电网电量不足,比较购电与燃料电池发电成本,选择价格较低者。
3)若购电成本低于燃料电池发电成本,交互功率为储能单元充放电后微电网功率缺额。
4)若燃料电池发电成本低于购电成本,其发电功率为储能充电之前或放电之后系统功率缺额。
为动态地表示上述原则,定义仅考储能、交互功率和燃料电池运行的期望成本函数如下:
式中:Fin_t1(t)、Fin_t2(t)分别表示燃料电池不启动下售电与购电的成本函数;Ptex(t)=Pe(t)+Pes(t),表示在燃料电池不启动时系统广义交互功率,正值表示微电网售电功率,负值表示从大电网购电功率。
式中:Prsmax为调度周期内最高售电电价;Fa(t)、Fc(t)、Fe(t)为各情况下非峰时放电期望损失;Fb(t)、Fd(t)、Fd(t)分别为各情况下非峰时充电期望收益;
燃料电池的出力为
其中:当PFC(t)=-Pe(t)时,Pes(t)=0。
微电网与大电网的交互功率为
由式(71)、(72)所示,在计算电能交互收益时,在广义负荷功率为正或负的情况下,对储能单元充放电过程分别考虑了充电的期望收益与放电过程中期望损失,即若非峰时电价放电时,售电收益减去的非峰时放电期望损失,若非峰时充电时,售电收益减去非峰时充电期望收益。其中,Fa(t)、Fb(t)、Ff(t)表示系统功率过剩时,储能系统充放电与替代供电方式之间的期望损失或收益,而式Fc(t)、Fd(t)、Fe(t)表示系统功率不足时,储能系统充放电与替代供电方式之间的期望损失或收益。因此,储能单元在动态寻优的任一时刻,将尽量避开非峰时放电与非谷时充电状态,并依据微电网前后时刻广义负荷功率Pe(t)、储能单元的剩余电量、实时电价和燃料电池发电成本决策最佳充放电出力,以保证系统期望收益函数取最大,满足运行原则1),2)。由(71)、(72)表示,当系统广义交互功率Ptex(t)为负时,比较从大电网购电与燃料电池成本,从而确定是否启动燃料电池及其出力大小,同时由式(29)可知,燃料电池仅满足系统的功率缺额,不向储能单元充电及大电网售电,满足运行原则3)、4)。
储能单元、交互功率及燃料电池的运行原则确定后,调度周期内整个CCHP型微电网的最佳运行方式随之确定。
4、算例分析
4.1 CCHP型微电网的结构
本文在保证系统安全、综合考虑系统各部分比例的基础上,以北方某地区实际的小型CCHP型微电网为例对其进行运行优化。微电网结构、联供方式及能量流向等如图1所示。
图1中,主网电压为10kv,微电网电压为0.38kv,初始时刻储能剩余电量值对应于电池最低核电状态;网内电负荷包括普通电负荷和空调电负荷,当微燃机不满足运行条件时,需启动空调机填补冷热负荷功率缺额。微电网内部分元件基础参数及排放物污染惩罚参数见如表1。
表1微电网内各元件参数
表2污染物排放和惩罚因子
选取冬季和夏季典型日为研究对象,时间间隔Δt=1h,NT=24,交互电价谷时段为23:00—07:00,平时段为07:00—08:00、11:00—17:00和22:00—23:00,峰时段为08:00—11:00和17:00—22:00。各时段售电、购电电价如表2所示。
表3调度周期实时电价
在研究冬夏季CCHP型微电网的电能调度时,制定以下三种调度策略:
策略1:储能单元停用,由交互功率和燃料电池出力满足微电网负荷需求。
策略2:储能单元不参与电能交易,仅为满足微电网的内部负荷需求而启用。
策略3:采用3.2节的方法,即储能单元参与电能交易,并动态规划储能充放电。
4.2典型日CCHP型微电网运行分析
本文冷热电负荷选取该微电网内冬夏季典型日负荷,并采取基于模拟退火的粒子群算法对模型优化求解。该微电网内冬夏季典型日冷热电负荷需求功率分别如图2、图3所示。
风电及光伏发电出力采用最大功率点追踪模式,冷热电联产系统运行在以热(冷)定电模式下并按3.1节策略运行,取微电网制热和制冷收益均为0.1元/(kWh),微电网内固定成本与收益如表4所示。
表4微电网的固定成本与收益
以冬季典型日为例,当CCHP型微电网与电网间交互取实时电价时,系统采取不同控制策略,交互功率与燃料电池的出力如图4、图5所示。
图4中,相对于策略1,策略2和3的储能能够在一定程度上平抑微电网与大电网之间的交互功率的波动,减少交互功率为负的时段,从而减少微电网购电量,为电网层的调度提供了更好的环境。相对策略2,策略3交互功率在电价较低时,售电功率低而购电功率高;电价较高时,售电功率高而购电功率低,增加了微电网的交互收益。图5中,相对于策略1和2,策略3由于采取了文中的储能优化策略,有效减少了燃料电池的启停次数并能极大发挥峰时电价下燃料电池代替交互功率满足系统功率缺额的作用。
策略2和3中储能的出力如图6所示,图6中,策略2在微电网功率不足时,储能放电,在微电网功率过剩时,储能充电,储能起到了削峰填谷的作用。相对策略2,策略3中储能除了削峰填谷外,在分时电价下,储能参与微电网电能经济调度并尽量在电价较低的时充电,而在电价较高的时放电:1-4时,电网电价较低,储能以最大功率充电,购进低价电能;9-11时,电网电价较高,储能以最大功率放电,获得最大收益。12-17时,微电网功率不足,此时电网电价相对较低,微电网优先购进大电网低价电能,储能单元不放电。18-20时,微电网能量不足,此时电价相对较高,微电网优先利用储能放电。由于对储能单元的充放电依据分时电价进行了经济调度,因此相对策略2,策略3还能有效减少储能的充放电次数,降低了充放电损失,从而整体提高微电网经济效益,降低微电网运行成本。
冬夏季典型日策略3的电负荷平衡情况如图7、图8所示,图7、图8中,调度周期内风电及光伏发电置于最大功率点追踪模式,CCHP系统运行在以热(冷)定电模式且微型燃气轮机与空调互补出力。冬季典型日,0-10和16-23时,热负荷功率较大,热电联产系统满足最低出力限制,微型燃气轮机启动;而在11-15时,系统热负荷较小,微型燃气轮机出力低于最低出力限制,微型燃气轮机停用,空调机启动,空调机出力为网内空调电负荷。由于夏季典型日冷负荷较大,调度周期内冷电联产系统均满足最低出力限制,微型燃气轮机启用,空调及停用,网内空调电负荷为0。风光出力及热电联产系统出力确定后,微电网的功率缺额主要由储能单元、交互功率及燃料电池提供。
三种策略下CCHP型微电网发电成本如表5所示,由于考虑了制冷制热收益,对系统的热能经济效益充分利用,相对于直接由电能供应冷热负荷的分产型微电网而言,运行成本大大降低。策略2、3在引入储能单元后,降低了微电网的运行成本,证明了储能系统的经济性。策略2中购电成本和售电收益均较低,表明储能单元对电网功率进行了平抑,但是由于没有考虑储能单元的经济调度策略,优化结果仍然不够经济。在分时电价下,采用本文的策略能够合理的优化储能单元的运行方式,提高微电网的售电收益,从而降低微电网的运行成本。
表5三种策略下微电网发电成本
4.3 CCHP系统指标评估
1)能源利用率指标。CCHP系统能源利用率如表6所示。冬夏季联产系统能源利用率均满足要求,且相对于传统火电厂、热电厂能源利用率(分别为35%和45%)有较大提高,实现了能源梯级高效利用。同时,由于溴冷机夏季制冷系数比冬季制热系数大,在相同的燃料产热量下,能够提供更多制冷量,因此冷电联产系统能源利用率高于热电联产系统,该结果可以为微电网的冷热调度及节能效益提供参考。
表6联产系统能源利用率
2)节能减排指标。CCHP系统的减排量、减排节约成本如表7所示。热电分产的火电机组制热量由余热锅炉供热提供,制冷量由压缩式制冷机组提供。联产系统由于主要燃料为天然气,单位含碳量相对较高,因此对提高CO2的减排率作用较小,同时由于热电联产中锅炉产热效率低于电压缩机的制冷效率,因此冷电联产系统的CO2排放较传统分产系统较高,而热电联产系统低于分产系统。但由于天然气所含SO2、NOx等污染物较少,因此能够极大地提高SO2和NOx的减排率,并减少排污惩罚费用,该结果可以为能源利用的环境代价活效益的定量化评价提供可操作性方法。
表7联产系统减排率和减排节约成本
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,包括以下步骤:
(1)假设在调度周期单位时间间隔Δt内各微电源的出力恒定并满足冷、热、电负荷的平衡关系,微电网与主网之间功率交互功率恒定,交互电价跟随实时电价并在Δt内保持不变的情况下,建立CCHP型微电网优化模型;
(2)根据功率平衡约束条件、分布式电源有功出力约束条件、空调机出力约束条件以及储能单元运行约束条件,构建以微电网并网运行状态下发电成本最低为目标的目标优化函数;
(3)根据建立的目标优化函数,利用基于模拟退火粒子群的动态规划法求解输出储能单元、交互功率和燃料电池的最佳出力;综合考虑前后时段的负荷水平以及各微电源的出力趋势,确定储能单元以及燃料电池的运行原则,进而确定调度周期内整个CCHP系统的最佳运行方式。
2.如权利要求1所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,所述步骤(1)中建立的CCHP型微电网优化模型包括:微型燃气轮机模型、空调机模型、燃料电池模型、风力发电模型、光伏发电模型和储能单元模型;
假设微型燃气轮机排烟温度与溴冷机的进烟温度保持不变,并忽略外界环境的变化对微型燃气轮机供电、燃料燃烧效率的影响,所述微型燃气轮机模型的数学模型为:
QMT-h(t)=QMT(t)ηhCOPh (2)
QMT-c(t)=QMT(t)ηcCOPc (3)
微型燃气轮机的燃料成本为:
式中,令t=t0+rΔt,t0为调度初始时刻,r为非负整数,调度周期T=NTΔt,NT为正整数;QMT(t)为燃气轮机排气余热量;PMT(t)为t时刻燃气轮机输出的电功率;ηMT为燃气轮机发电效率;ηL为散热损失系数;QMT-h(t)、QMT-c(t)为溴冷机t时刻提供的制热量和制冷量;ηh、ηc和COPh、COPc分别为溴冷机的制热和制冷时的烟气回收率及制热和制冷系数;;微型燃气轮机不同功率下的发电效率由生产商提供的相关参数拟合得到;CMT为微型燃气轮机的燃料成本;Cf为单位燃料的价格;LHVNG为天然气低热热值;
所述空调机模型的数学函数为:
Qair-h(t)=Pair(t)(1-ηLh)ηahCOPair-h (5)
Qair-c(t)=Pair(t)(1-ηLc)ηacCOPair-c (6)
式中,Qair-h(t)、Qair-c(t)分别为t时刻空调机制热、制冷功率;Pair(t)为空调机用电功率;ηLh、ηLc分别为空调机制热、制冷时损失系数,ηah、ηac分别为空调机的制热、制冷效率;COPair-h、COPair-c分别为空调制热和制冷时的能效比;空调机所耗电能作为系统用电负荷统一调度。
3.如权利要求2所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,所述燃料电池模型的数学函数为:
不考虑燃料电池的热能利用,燃料电池的燃料消耗-功率输出特性为:
式中,CFC为燃料电池的燃料成本;PFC(t)为t时刻燃料电池发电功率;ηFC为燃料电池发电效率;LHVNG为天然气低热热值,Cf为单位燃料价格、r为调度时间段、NT为总的调度时间段数目,;燃料电池不同功率下的发电效率由生产商提供的相关参数拟合得到;
所述储能单元模型的状态有充电、放电、浮充三种状态,储能在t时刻的剩余电量与t-1时刻的剩余电量、t-1到t时刻的充放电量以及自放电量有关;储能的输入、输出功率及充放电状态满足下列关系:
式中,Csoc(t)为t时刻储能剩余电量;τ为储能自放电率;Pes为储能t时刻的充放电功率;ηch、ηdis分别为储能充电和放电效率;ufl(t)、uch(t)、udis(t)∈﹛0,1﹜,为储能单元的状态标识,分别代表浮充,充电,放电三种状态,对应Pes(t)取值分别等于、小于、大于零,Δt为单位调度时间间隔。
4.如权利要求1所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,所述步骤(2)中构建以微电网并网运行状态下发电成本最低为目标的目标优化函数具体为:
minF=F1-F2 (9)
式中,F为微电网的调度周期内总发电成本;F1为微电源运行成本函数;F2为微电网经济收益函数;
微电源运行成本函数为:
F1=Ct+Cfu+Cm (10)
式中,Ct、Cfu、Cm分别代表微电网折旧成本、燃料成本和维护成本;N为微电源总数;Ci为第i台微电源的初始成本;l为利息率;ki为第i台微电源的寿命周期;CBAT为储能初始成本;v为储能寿命周期,微电源与储能折旧的平摊成本为期初成本;Cmi为第i台微电源的单位出力维护成本;Pi(t)为第i台微电源单位时间内有功出力;Cmes为储能的单位出力维护成本,CMT为微型燃气轮机的燃料成本、CFC为燃料电池的成本、NT为调度周期间隔数、Pes(t)为储能t时刻的充放电功率。
5.如权利要求4所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,所述微电网经济收益函数为:
式中,Prs(t)、Prb(t)分别为t时刻的售电和购电电价;Ps(t)、Pb(t)分别为t时刻的售电和购电功率;Chc为制热(冷)收益;Che、Cco分别为单位制热和制冷量的售价;a=1表示冬季制热,a=0表示夏季制冷;QMT-h(t)为微型燃气轮机的产热量、QMT-c(t)为微型燃气轮机的制冷量。
6.如权利要求1所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,所述步骤(2)中的功率平衡约束条件具体为:
式中,Pload(t)、Qhe(t)、Qco(t)分别为t时刻网内普通电负荷和热、冷负荷;Ploss(t)为t时刻微电网功率损耗,QMT-h(t)为微型燃气轮机的产热量、QMT-c(t)为微型燃气轮机的制冷量,PMT(t)为微型燃气轮机的发电功率,Pair(t)为空调机用电功率;Pb(t)、Ps(t)分别为微电网从大电网的购电功率和售电功率。,Qair-h(t)、Qair-c(t)分别为t时刻空调机制热、制冷功率。Pi(t)为第i个微电源的出力功率。
7.如权利要求1所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,所述步骤(2)中的分布式电源有功出力约束条件具体为:
Pi min≤Pi(t)≤Pi max (17)
式中,Pi max,Pi min分别为第i台微电源有功出力的上下限;
所述空调机出力约束条件具体为:
Pair min≤Pair(t)≤Pair max (18)
式中,Pair min、Pair max为空调机出力上下限。
8.如权利要求1所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,所述步骤(2)中的储能单元运行约束条件具体为:
①储能单元的出力功率约束条件:
Pes min≤Pes(t)≤Pes max (19)
式中:Pes max,Pes min分别为储能出力上下限;
②储能单元的能量约束条件:
Csoc min≤Csoc(t)≤Csoc max (20)
式中:Csoc max,Csoc min分别为储能容量的最大、最小值;其中,SOC表示储能的核电状态,表征其能量范围;
③单位调度周期始末状态储能能量约束条件:
Csoc(t0)=Csoc(t0+NrΔt) (21)
考虑储能单元对微电网的周期性优化,其能量状态满足在调度周期始末相等;Csoc(t)为t时刻储能剩余电量;Δt为单位调度时间间隔,NT为总的调度时间段数目;
单位调度周期储能的充放电状态和次数约束:
uch(t)+udis(t)+ufl(t)=1 (22)
式中:λ1、λ2分别表示储能的充放电次数,具体的数值由系统负荷、储能寿命和优化策略来确定;uch(t)、udis(t)、ufl(t)分别为代表t时刻储能的充电、放电、浮充三种状态,Δt为单位调度时间间隔。
9.如权利要求1所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,所述步骤(3)中,调度周期内整个CCHP系统的最佳运行方式的确定原则为:
空调机与微型燃气轮机单独使用;
定义Pe(t)为t时刻储能和燃料电池均不启动时微电网广义负荷功率,正值表示微电网输出电能,负值表示微电网缺电并假设大电网能完全消纳与弥补微电网功率差额;
依据Pe(t)与分时电价制定的储能单元交换功率和燃料电池的运行原则为:
1)若微电网售电价格较低,储能单元优先充电:若微电网电量剩余,则向大电网售出电能,燃料电池不工作;若微电网电量不足,比较购电与燃料电池发电成本,选择成本较低者;
2)若微电网售电价格较高,储能单元优先放电:若微电网电量剩余,则向大电网售出电能,燃料电池不工作;若微电网电量不足,比较购电与燃料电池发电成本,选择价格较低者;
3)若购电成本低于燃料电池发电成本,交互功率为储能单元充放电后微电网功率缺额;
4)若燃料电池发电成本低于购电成本,其发电功率为储能充电之前或放电之后系统功率缺额。
10.如权利要求9所述的一种考虑冷热电联供和储能运行策略的微电网运行优化方法,其特征是,定义仅考虑储能单元交换功率和燃料电池运行的期望成本函数如下:
式中,Fin_t1(t)、Fin_t2(t)分别表示燃料电池不启动下售电与购电的成本函数;Ptex(t)=Pe(t)+Pes(t),表示在燃料电池不启动时系统广义交互功率,正值表示微电网售电功率,负值表示从大电网购电功率;Pe(t)、为t时刻为t时刻储能和燃料电池均不启动时微电网广义负荷功率,Pes(t)为储能的充放电功率;
式中,Prs max为调度周期内最高售电电价;Fa(t)、Fc(t)、Fe(t)为各情况下非峰时放电期望损失;Fb(t)、Fd(t)、Ff(t)分别为各情况下非峰时充电期望收益;Ptex(t)表示在燃料电池不启动时系统广义交互功率、Prs(t)、Prb(t)分别为t时刻的售电和购电电价、Δt为单位调度时间间隔、CFC(t)为燃料电池的发电成泵、ηch和ηdis分别为储能的充、放电功率;
燃料电池的出力为
其中:当PFC(t)=-Pe(t)时,Pes(t)=0;
微电网与大电网的交互功率为
由式(26)、(27)所示,在计算电能交互收益时,在广义负荷功率为正或负的情况下,对储能单元充放电过程分别考虑了充电的期望收益与放电过程中期望损失,即若非峰时电价放电时,售电收益减去的非峰时放电期望损失,若非峰时充电时,售电收益减去非峰时充电期望收益;
其中,Fa(t)、Fb(t)、Ff(t)表示系统功率过剩时,储能系统充放电与替代供电方式之间的期望损失或收益,而式Fc(t)、Fd(t)、Fe(t)表示系统功率不足时,储能系统充放电与替代供电方式之间的期望损失或收益;
因此,储能单元在动态寻优的任一时刻,将尽量避开非峰时放电与非谷时充电状态,并依据微电网前后时刻广义负荷功率Pe(t)、储能单元的剩余电量、实时电价和燃料电池发电成本决策最佳充放电出力,以保证系统期望收益函数取最大,满足运行原则1)和2);
由(29)、(30)表示,当系统广义交互功率Ptex(t)为负时,比较从大电网购电与燃料电池成本,从而确定是否启动燃料电池及其出力大小,同时由式(29)可知,燃料电池仅满足系统的功率缺额,不向储能单元充电及大电网售电,满足运行原则3)和4);
储能单元以及燃料电池的运行原则确定后,调度周期内整个CCHP型微电网的最佳运行方式随之确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410724835.1A CN104392286B (zh) | 2014-12-02 | 2014-12-02 | 考虑冷热电联供和储能运行策略的微电网运行优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410724835.1A CN104392286B (zh) | 2014-12-02 | 2014-12-02 | 考虑冷热电联供和储能运行策略的微电网运行优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104392286A true CN104392286A (zh) | 2015-03-04 |
CN104392286B CN104392286B (zh) | 2017-07-21 |
Family
ID=52610186
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410724835.1A Expired - Fee Related CN104392286B (zh) | 2014-12-02 | 2014-12-02 | 考虑冷热电联供和储能运行策略的微电网运行优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104392286B (zh) |
Cited By (65)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104734168A (zh) * | 2015-03-13 | 2015-06-24 | 山东大学 | 一种基于电热联合调度的微电网运行优化系统及方法 |
CN104820363A (zh) * | 2015-03-09 | 2015-08-05 | 山东大学 | 基于生物质气的冷热电联供系统及能量优化控制方法 |
CN105391061A (zh) * | 2015-12-01 | 2016-03-09 | 贵州电网有限责任公司 | 冷热电风光储联合优化的主动配电网优化系统和优化方法 |
CN105931136A (zh) * | 2016-04-25 | 2016-09-07 | 天津大学 | 一种融合需求侧虚拟储能系统的楼宇微网优化调度方法 |
CN105930980A (zh) * | 2016-06-08 | 2016-09-07 | 河海大学 | 一种电转气的综合能源系统多点线性概率能量流方法 |
CN106058942A (zh) * | 2016-08-03 | 2016-10-26 | 长沙理工大学 | 考虑风电不确定性的含电转气和cchp的能量枢纽优化模型 |
CN106712086A (zh) * | 2017-01-17 | 2017-05-24 | 无锡协鑫分布式能源开发有限公司 | 微网优化控制模式 |
CN107023933A (zh) * | 2015-10-27 | 2017-08-08 | 阿自倍尔株式会社 | 热源运行辅助装置及方法 |
CN107065523A (zh) * | 2017-01-17 | 2017-08-18 | 无锡协鑫分布式能源开发有限公司 | 一种多能互补微能源的能效运行策略 |
CN107203136A (zh) * | 2017-06-08 | 2017-09-26 | 国网甘肃省电力公司电力科学研究院 | 一种智慧农业大棚微型能源网的优化调度方法及装置 |
CN107294201A (zh) * | 2017-06-02 | 2017-10-24 | 北京科诺伟业科技股份有限公司 | 一种同步发电机和储能变流器并联模拟平台 |
CN107769246A (zh) * | 2017-10-20 | 2018-03-06 | 海南电网有限责任公司 | 一种孤立岛屿的微电网控制方法及系统 |
CN107807523A (zh) * | 2017-10-18 | 2018-03-16 | 国网天津市电力公司电力科学研究院 | 考虑分时电价的区域能源互联网多源协调优化运行策略 |
CN107844055A (zh) * | 2017-11-03 | 2018-03-27 | 南京国电南自电网自动化有限公司 | 一种基于博弈论的冷热电三联供微网系统优化运行方法 |
CN107863773A (zh) * | 2017-11-24 | 2018-03-30 | 清华-伯克利深圳学院筹备办公室 | 一种多能流供能方法及系统 |
CN107918290A (zh) * | 2017-10-19 | 2018-04-17 | 上海科梁信息工程股份有限公司 | 一种仿真系统及仿真方法 |
CN108133285A (zh) * | 2017-12-21 | 2018-06-08 | 南京工程学院 | 一种接入大规模可再生能源的混合能源系统实时调度方法 |
CN108171384A (zh) * | 2017-12-30 | 2018-06-15 | 国网天津市电力公司电力科学研究院 | 一种基于复合粒子群算法微网能量管理方法 |
CN108387013A (zh) * | 2018-01-09 | 2018-08-10 | 湘潭大学 | 天然气-地热能-太阳能耦合的多能源互补系统 |
CN108491973A (zh) * | 2018-03-21 | 2018-09-04 | 浙江大学 | 一种考虑需求侧响应的电-热联合系统优化调度方法 |
CN108564230A (zh) * | 2018-04-28 | 2018-09-21 | 湖南红太阳新能源科技有限公司 | 一种户用分布式能源管理方法及系统 |
CN108564242A (zh) * | 2018-01-09 | 2018-09-21 | 湖南大学 | 微能源网系统、微能源网配置方法及装置 |
CN108631343A (zh) * | 2018-06-12 | 2018-10-09 | 上海电力学院 | 一种多能互补能源互联网优化调度方法 |
CN108649553A (zh) * | 2018-03-23 | 2018-10-12 | 中国电力科学研究院有限公司 | 一种微能源网供能系统的配置选优方法及装置 |
CN108808653A (zh) * | 2017-05-02 | 2018-11-13 | 南京理工大学 | 一种考虑可控负荷的风光储微电网储能容量优化配置方法 |
CN108832656A (zh) * | 2018-06-20 | 2018-11-16 | 天津大学 | 基于电转气和可再生能源利用微能源网多目标规划方法 |
CN108964048A (zh) * | 2018-08-15 | 2018-12-07 | 国家电网有限公司 | 微电网容量配置与运行调度方法及装置 |
CN108960564A (zh) * | 2018-05-23 | 2018-12-07 | 国网天津市电力公司电力科学研究院 | 基于集中式建模的冷热电联供型微网配电系统调度方法 |
CN108984894A (zh) * | 2018-07-09 | 2018-12-11 | 四川大学 | 一种基于ccp的ehs气电购置优化方法 |
CN109066750A (zh) * | 2018-09-11 | 2018-12-21 | 重庆大学 | 基于需求侧响应的光伏-电池微电网混合能量调度管理方法 |
CN109217349A (zh) * | 2017-06-30 | 2019-01-15 | 国网山西省电力公司经济技术研究院 | 基于多微网并网场景下的计及运行策略的微电网外特性求解方法 |
CN109214561A (zh) * | 2018-08-21 | 2019-01-15 | 上海电力学院 | 考虑主动配电系统动态路径优化的分布式电源配置方法 |
CN109255471A (zh) * | 2018-08-17 | 2019-01-22 | 国网山东省电力公司电力科学研究院 | 一种含风电的电-气-热综合能源系统扩展规划优化方法 |
CN109255560A (zh) * | 2018-11-20 | 2019-01-22 | 成都大学 | 一种基于冷热电负荷比例的cchp系统评估优化方法 |
CN109286187A (zh) * | 2018-10-19 | 2019-01-29 | 国网宁夏电力有限公司经济技术研究院 | 一种面向多主体利益均衡的微网日前经济调度方法 |
CN109345019A (zh) * | 2018-10-10 | 2019-02-15 | 南京邮电大学 | 一种基于改进粒子群算法的微电网经济调度优化策略 |
CN109473976A (zh) * | 2018-10-22 | 2019-03-15 | 华润智慧能源有限公司 | 一种冷热电联供型微网能量调度方法及系统 |
CN109472716A (zh) * | 2018-11-09 | 2019-03-15 | 浙江大学 | 一种基于遗传算法的考虑多因素的产能端功率分配方法 |
CN109659927A (zh) * | 2018-10-24 | 2019-04-19 | 国网天津市电力公司电力科学研究院 | 一种考虑储能参与度的综合能源微网储能容量配置方法 |
CN109687504A (zh) * | 2018-11-14 | 2019-04-26 | 杭州赫智电子科技有限公司 | 基于源荷平衡的微电网能源优化分配方 |
CN109687518A (zh) * | 2018-12-29 | 2019-04-26 | 南京工程学院 | 一种家庭微电网系统的优化调度方法 |
CN109802447A (zh) * | 2019-01-30 | 2019-05-24 | 广东工业大学 | 一种cchp系统调度评价方法 |
CN109861260A (zh) * | 2018-11-30 | 2019-06-07 | 中国电力科学研究院有限公司 | 一种混合能源系统中电能存储设备的充放电控制方法及装置 |
CN109919478A (zh) * | 2019-02-28 | 2019-06-21 | 天津大学 | 一种考虑综合供能可靠性的综合能源微网规划方法 |
CN110190598A (zh) * | 2019-06-06 | 2019-08-30 | 新奥数能科技有限公司 | 获取冷热电联供多微网系统的优化调度方案的方法及系统 |
CN110263966A (zh) * | 2019-05-06 | 2019-09-20 | 天津大学 | 考虑动态传热过程的电-热综合能源系统优化调度方法 |
CN110264051A (zh) * | 2019-06-06 | 2019-09-20 | 新奥数能科技有限公司 | 确定微网的优化调度方案的方法及装置 |
CN110391655A (zh) * | 2018-04-20 | 2019-10-29 | 中国电力科学研究院有限公司 | 一种含多能源耦合的微能源网经济优化调度方法及装置 |
CN110504684A (zh) * | 2019-08-21 | 2019-11-26 | 东北大学 | 一种区域多微网系统日前优化调度方法 |
CN110689189A (zh) * | 2019-09-24 | 2020-01-14 | 国网天津市电力公司 | 考虑供能侧和需求侧的冷热电联合供需平衡优化调度方法 |
CN111049138A (zh) * | 2019-12-31 | 2020-04-21 | 四川大学 | 一种基于云储能系统的微电网多源协调优化的方法和装置 |
CN111275251A (zh) * | 2020-01-16 | 2020-06-12 | 广东工业大学 | 一种含污水源热泵的cchp系统冷热电联供优化方法 |
CN111313400A (zh) * | 2019-11-11 | 2020-06-19 | 国网吉林省电力有限公司 | 一种基于鲁棒修正的多能虚拟电厂运行参数聚合方法 |
CN111339474A (zh) * | 2020-02-17 | 2020-06-26 | 山东大学 | 一种基于趋势预测分析法的综合能源系统预测运行方法 |
CN112103955A (zh) * | 2020-09-16 | 2020-12-18 | 湖南大学 | 一种综合能源系统的电储能事故备用容量优化利用方法 |
CN112103946A (zh) * | 2020-08-20 | 2020-12-18 | 西安理工大学 | 一种基于粒子群算法的微电网储能优化配置方法 |
CN112103940A (zh) * | 2020-07-03 | 2020-12-18 | 国网江西省电力有限公司电力科学研究院 | 一种含温控设备的多能源微电网最优经济运行方法 |
CN112488363A (zh) * | 2020-11-09 | 2021-03-12 | 国网湖北省电力有限公司电力科学研究院 | 基于广义储能的多能源电力系统优化调度方法 |
CN112713590A (zh) * | 2020-12-22 | 2021-04-27 | 南昌大学 | 计及idr的冷热电联供微网与主动配电网联合优化调度方法 |
CN113065981A (zh) * | 2021-03-30 | 2021-07-02 | 东北电力大学 | 一种考虑储能综合收益的源-储-网协调规划方法 |
CN113098036A (zh) * | 2021-03-23 | 2021-07-09 | 国网安徽省电力有限公司电力科学研究院 | 基于氢燃料电池的综合能源系统运行方法 |
CN113779757A (zh) * | 2021-08-05 | 2021-12-10 | 成都信息工程大学 | 一种电气热联供能源集线器日前优化调度方法 |
CN113883682A (zh) * | 2021-09-28 | 2022-01-04 | 青岛海尔中央空调有限公司 | 空调计费方法、装置及计算机可读存储介质 |
CN115062869A (zh) * | 2022-08-04 | 2022-09-16 | 国网山东省电力公司东营供电公司 | 一种考虑碳排放量的综合能源调度方法和系统 |
CN116646908A (zh) * | 2023-04-06 | 2023-08-25 | 中国质量认证中心 | 基于光储直柔配电系统的节能用电评估方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050209906A1 (en) * | 2002-04-30 | 2005-09-22 | Bayoumi Deia S | Distribution/power transformers customer support, tracking problems and recalls |
CN103617460A (zh) * | 2013-12-06 | 2014-03-05 | 天津大学 | 冷热电联供微网系统的双层优化规划设计方法 |
-
2014
- 2014-12-02 CN CN201410724835.1A patent/CN104392286B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050209906A1 (en) * | 2002-04-30 | 2005-09-22 | Bayoumi Deia S | Distribution/power transformers customer support, tracking problems and recalls |
CN103617460A (zh) * | 2013-12-06 | 2014-03-05 | 天津大学 | 冷热电联供微网系统的双层优化规划设计方法 |
Non-Patent Citations (3)
Title |
---|
万术来: "基于改进粒子群算法的微网环保经济运行的优化", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 * |
刘天琪等: "基于储能单元运行方式优化的微电网经济运行", 《电网技术》 * |
闫晓娟等: "冷热电联供系统能耗与CO2排放特性分析", 《建筑热能通风空调》 * |
Cited By (91)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104820363A (zh) * | 2015-03-09 | 2015-08-05 | 山东大学 | 基于生物质气的冷热电联供系统及能量优化控制方法 |
CN104820363B (zh) * | 2015-03-09 | 2017-12-12 | 山东大学 | 基于生物质气的冷热电联供系统及能量优化控制方法 |
CN104734168A (zh) * | 2015-03-13 | 2015-06-24 | 山东大学 | 一种基于电热联合调度的微电网运行优化系统及方法 |
CN107023933A (zh) * | 2015-10-27 | 2017-08-08 | 阿自倍尔株式会社 | 热源运行辅助装置及方法 |
CN107023933B (zh) * | 2015-10-27 | 2019-07-23 | 阿自倍尔株式会社 | 热源运行辅助装置及方法 |
CN105391061A (zh) * | 2015-12-01 | 2016-03-09 | 贵州电网有限责任公司 | 冷热电风光储联合优化的主动配电网优化系统和优化方法 |
CN105931136A (zh) * | 2016-04-25 | 2016-09-07 | 天津大学 | 一种融合需求侧虚拟储能系统的楼宇微网优化调度方法 |
CN105930980A (zh) * | 2016-06-08 | 2016-09-07 | 河海大学 | 一种电转气的综合能源系统多点线性概率能量流方法 |
CN105930980B (zh) * | 2016-06-08 | 2019-10-15 | 河海大学 | 一种电转气的综合能源系统多点线性概率能量流方法 |
CN106058942A (zh) * | 2016-08-03 | 2016-10-26 | 长沙理工大学 | 考虑风电不确定性的含电转气和cchp的能量枢纽优化模型 |
CN106058942B (zh) * | 2016-08-03 | 2020-01-10 | 长沙理工大学 | 考虑风电不确定性的含电转气和cchp的能量枢纽优化模型 |
CN107065523A (zh) * | 2017-01-17 | 2017-08-18 | 无锡协鑫分布式能源开发有限公司 | 一种多能互补微能源的能效运行策略 |
CN106712086A (zh) * | 2017-01-17 | 2017-05-24 | 无锡协鑫分布式能源开发有限公司 | 微网优化控制模式 |
CN108808653A (zh) * | 2017-05-02 | 2018-11-13 | 南京理工大学 | 一种考虑可控负荷的风光储微电网储能容量优化配置方法 |
CN107294201A (zh) * | 2017-06-02 | 2017-10-24 | 北京科诺伟业科技股份有限公司 | 一种同步发电机和储能变流器并联模拟平台 |
CN107203136A (zh) * | 2017-06-08 | 2017-09-26 | 国网甘肃省电力公司电力科学研究院 | 一种智慧农业大棚微型能源网的优化调度方法及装置 |
CN109217349A (zh) * | 2017-06-30 | 2019-01-15 | 国网山西省电力公司经济技术研究院 | 基于多微网并网场景下的计及运行策略的微电网外特性求解方法 |
CN107807523A (zh) * | 2017-10-18 | 2018-03-16 | 国网天津市电力公司电力科学研究院 | 考虑分时电价的区域能源互联网多源协调优化运行策略 |
CN107918290A (zh) * | 2017-10-19 | 2018-04-17 | 上海科梁信息工程股份有限公司 | 一种仿真系统及仿真方法 |
CN107769246A (zh) * | 2017-10-20 | 2018-03-06 | 海南电网有限责任公司 | 一种孤立岛屿的微电网控制方法及系统 |
CN107844055A (zh) * | 2017-11-03 | 2018-03-27 | 南京国电南自电网自动化有限公司 | 一种基于博弈论的冷热电三联供微网系统优化运行方法 |
CN107863773B (zh) * | 2017-11-24 | 2020-04-24 | 清华-伯克利深圳学院筹备办公室 | 一种多能流供能方法及系统 |
CN107863773A (zh) * | 2017-11-24 | 2018-03-30 | 清华-伯克利深圳学院筹备办公室 | 一种多能流供能方法及系统 |
CN108133285B (zh) * | 2017-12-21 | 2021-05-11 | 南京工程学院 | 一种接入大规模可再生能源的混合能源系统实时调度方法 |
CN108133285A (zh) * | 2017-12-21 | 2018-06-08 | 南京工程学院 | 一种接入大规模可再生能源的混合能源系统实时调度方法 |
CN108171384A (zh) * | 2017-12-30 | 2018-06-15 | 国网天津市电力公司电力科学研究院 | 一种基于复合粒子群算法微网能量管理方法 |
CN108564242A (zh) * | 2018-01-09 | 2018-09-21 | 湖南大学 | 微能源网系统、微能源网配置方法及装置 |
CN108387013A (zh) * | 2018-01-09 | 2018-08-10 | 湘潭大学 | 天然气-地热能-太阳能耦合的多能源互补系统 |
CN108491973A (zh) * | 2018-03-21 | 2018-09-04 | 浙江大学 | 一种考虑需求侧响应的电-热联合系统优化调度方法 |
CN108649553A (zh) * | 2018-03-23 | 2018-10-12 | 中国电力科学研究院有限公司 | 一种微能源网供能系统的配置选优方法及装置 |
CN110391655A (zh) * | 2018-04-20 | 2019-10-29 | 中国电力科学研究院有限公司 | 一种含多能源耦合的微能源网经济优化调度方法及装置 |
CN108564230A (zh) * | 2018-04-28 | 2018-09-21 | 湖南红太阳新能源科技有限公司 | 一种户用分布式能源管理方法及系统 |
CN108960564A (zh) * | 2018-05-23 | 2018-12-07 | 国网天津市电力公司电力科学研究院 | 基于集中式建模的冷热电联供型微网配电系统调度方法 |
CN108631343A (zh) * | 2018-06-12 | 2018-10-09 | 上海电力学院 | 一种多能互补能源互联网优化调度方法 |
CN108832656A (zh) * | 2018-06-20 | 2018-11-16 | 天津大学 | 基于电转气和可再生能源利用微能源网多目标规划方法 |
CN108832656B (zh) * | 2018-06-20 | 2023-08-22 | 天津大学 | 基于电转气和可再生能源利用微能源网多目标规划方法 |
CN108984894A (zh) * | 2018-07-09 | 2018-12-11 | 四川大学 | 一种基于ccp的ehs气电购置优化方法 |
CN108964048A (zh) * | 2018-08-15 | 2018-12-07 | 国家电网有限公司 | 微电网容量配置与运行调度方法及装置 |
CN108964048B (zh) * | 2018-08-15 | 2020-07-28 | 国家电网有限公司 | 微电网容量配置与运行调度方法及装置 |
CN109255471A (zh) * | 2018-08-17 | 2019-01-22 | 国网山东省电力公司电力科学研究院 | 一种含风电的电-气-热综合能源系统扩展规划优化方法 |
CN109214561A (zh) * | 2018-08-21 | 2019-01-15 | 上海电力学院 | 考虑主动配电系统动态路径优化的分布式电源配置方法 |
CN109066750A (zh) * | 2018-09-11 | 2018-12-21 | 重庆大学 | 基于需求侧响应的光伏-电池微电网混合能量调度管理方法 |
CN109345019B (zh) * | 2018-10-10 | 2021-08-31 | 南京邮电大学 | 一种基于改进粒子群算法的微电网经济调度优化策略 |
CN109345019A (zh) * | 2018-10-10 | 2019-02-15 | 南京邮电大学 | 一种基于改进粒子群算法的微电网经济调度优化策略 |
CN109286187A (zh) * | 2018-10-19 | 2019-01-29 | 国网宁夏电力有限公司经济技术研究院 | 一种面向多主体利益均衡的微网日前经济调度方法 |
CN109286187B (zh) * | 2018-10-19 | 2022-01-04 | 国网宁夏电力有限公司经济技术研究院 | 一种面向多主体利益均衡的微网日前经济调度方法 |
CN109473976A (zh) * | 2018-10-22 | 2019-03-15 | 华润智慧能源有限公司 | 一种冷热电联供型微网能量调度方法及系统 |
CN109659927A (zh) * | 2018-10-24 | 2019-04-19 | 国网天津市电力公司电力科学研究院 | 一种考虑储能参与度的综合能源微网储能容量配置方法 |
CN109659927B (zh) * | 2018-10-24 | 2022-06-17 | 国网天津市电力公司电力科学研究院 | 一种考虑储能参与度的综合能源微网储能容量配置方法 |
CN109472716B (zh) * | 2018-11-09 | 2022-05-17 | 浙江大学 | 一种基于遗传算法的考虑多因素的产能端功率分配方法 |
CN109472716A (zh) * | 2018-11-09 | 2019-03-15 | 浙江大学 | 一种基于遗传算法的考虑多因素的产能端功率分配方法 |
CN109687504A (zh) * | 2018-11-14 | 2019-04-26 | 杭州赫智电子科技有限公司 | 基于源荷平衡的微电网能源优化分配方 |
CN109255560A (zh) * | 2018-11-20 | 2019-01-22 | 成都大学 | 一种基于冷热电负荷比例的cchp系统评估优化方法 |
CN109861260A (zh) * | 2018-11-30 | 2019-06-07 | 中国电力科学研究院有限公司 | 一种混合能源系统中电能存储设备的充放电控制方法及装置 |
CN109861260B (zh) * | 2018-11-30 | 2022-09-02 | 中国电力科学研究院有限公司 | 一种混合能源系统中电能存储设备的充放电控制方法及装置 |
CN109687518B (zh) * | 2018-12-29 | 2022-06-17 | 南京工程学院 | 一种家庭微电网系统的优化调度方法 |
CN109687518A (zh) * | 2018-12-29 | 2019-04-26 | 南京工程学院 | 一种家庭微电网系统的优化调度方法 |
CN109802447A (zh) * | 2019-01-30 | 2019-05-24 | 广东工业大学 | 一种cchp系统调度评价方法 |
CN109919478A (zh) * | 2019-02-28 | 2019-06-21 | 天津大学 | 一种考虑综合供能可靠性的综合能源微网规划方法 |
CN110263966B (zh) * | 2019-05-06 | 2023-04-18 | 天津大学 | 考虑动态传热过程的电-热综合能源系统优化调度方法 |
CN110263966A (zh) * | 2019-05-06 | 2019-09-20 | 天津大学 | 考虑动态传热过程的电-热综合能源系统优化调度方法 |
CN110190598A (zh) * | 2019-06-06 | 2019-08-30 | 新奥数能科技有限公司 | 获取冷热电联供多微网系统的优化调度方案的方法及系统 |
CN110264051A (zh) * | 2019-06-06 | 2019-09-20 | 新奥数能科技有限公司 | 确定微网的优化调度方案的方法及装置 |
CN110504684A (zh) * | 2019-08-21 | 2019-11-26 | 东北大学 | 一种区域多微网系统日前优化调度方法 |
CN110504684B (zh) * | 2019-08-21 | 2023-06-20 | 东北大学 | 一种区域多微网系统日前优化调度方法 |
CN110689189B (zh) * | 2019-09-24 | 2023-05-09 | 国网天津市电力公司 | 考虑供能侧和需求侧的冷热电联合供需平衡优化调度方法 |
CN110689189A (zh) * | 2019-09-24 | 2020-01-14 | 国网天津市电力公司 | 考虑供能侧和需求侧的冷热电联合供需平衡优化调度方法 |
CN111313400B (zh) * | 2019-11-11 | 2022-07-12 | 国网吉林省电力有限公司 | 一种基于鲁棒修正的多能虚拟电厂运行参数聚合方法 |
CN111313400A (zh) * | 2019-11-11 | 2020-06-19 | 国网吉林省电力有限公司 | 一种基于鲁棒修正的多能虚拟电厂运行参数聚合方法 |
CN111049138B (zh) * | 2019-12-31 | 2021-04-20 | 四川大学 | 一种基于云储能系统的微电网多源协调优化的方法和装置 |
CN111049138A (zh) * | 2019-12-31 | 2020-04-21 | 四川大学 | 一种基于云储能系统的微电网多源协调优化的方法和装置 |
CN111275251B (zh) * | 2020-01-16 | 2022-10-04 | 广东工业大学 | 一种含污水源热泵的cchp系统冷热电联供优化方法 |
CN111275251A (zh) * | 2020-01-16 | 2020-06-12 | 广东工业大学 | 一种含污水源热泵的cchp系统冷热电联供优化方法 |
CN111339474A (zh) * | 2020-02-17 | 2020-06-26 | 山东大学 | 一种基于趋势预测分析法的综合能源系统预测运行方法 |
CN112103940A (zh) * | 2020-07-03 | 2020-12-18 | 国网江西省电力有限公司电力科学研究院 | 一种含温控设备的多能源微电网最优经济运行方法 |
CN112103946A (zh) * | 2020-08-20 | 2020-12-18 | 西安理工大学 | 一种基于粒子群算法的微电网储能优化配置方法 |
CN112103955B (zh) * | 2020-09-16 | 2022-02-08 | 湖南大学 | 一种综合能源系统的电储能事故备用容量优化利用方法 |
CN112103955A (zh) * | 2020-09-16 | 2020-12-18 | 湖南大学 | 一种综合能源系统的电储能事故备用容量优化利用方法 |
CN112488363B (zh) * | 2020-11-09 | 2022-05-31 | 国网湖北省电力有限公司电力科学研究院 | 基于广义储能的多能源电力系统优化调度方法 |
CN112488363A (zh) * | 2020-11-09 | 2021-03-12 | 国网湖北省电力有限公司电力科学研究院 | 基于广义储能的多能源电力系统优化调度方法 |
CN112713590A (zh) * | 2020-12-22 | 2021-04-27 | 南昌大学 | 计及idr的冷热电联供微网与主动配电网联合优化调度方法 |
CN112713590B (zh) * | 2020-12-22 | 2022-11-08 | 南昌大学 | 计及idr的冷热电联供微网与主动配电网联合优化调度方法 |
CN113098036A (zh) * | 2021-03-23 | 2021-07-09 | 国网安徽省电力有限公司电力科学研究院 | 基于氢燃料电池的综合能源系统运行方法 |
CN113098036B (zh) * | 2021-03-23 | 2024-06-07 | 国网安徽省电力有限公司电力科学研究院 | 基于氢燃料电池的综合能源系统运行方法 |
CN113065981B (zh) * | 2021-03-30 | 2022-04-26 | 东北电力大学 | 一种考虑储能综合收益的源-储-网协调规划方法 |
CN113065981A (zh) * | 2021-03-30 | 2021-07-02 | 东北电力大学 | 一种考虑储能综合收益的源-储-网协调规划方法 |
CN113779757A (zh) * | 2021-08-05 | 2021-12-10 | 成都信息工程大学 | 一种电气热联供能源集线器日前优化调度方法 |
CN113883682A (zh) * | 2021-09-28 | 2022-01-04 | 青岛海尔中央空调有限公司 | 空调计费方法、装置及计算机可读存储介质 |
CN115062869A (zh) * | 2022-08-04 | 2022-09-16 | 国网山东省电力公司东营供电公司 | 一种考虑碳排放量的综合能源调度方法和系统 |
CN116646908A (zh) * | 2023-04-06 | 2023-08-25 | 中国质量认证中心 | 基于光储直柔配电系统的节能用电评估方法 |
CN116646908B (zh) * | 2023-04-06 | 2023-12-19 | 中国质量认证中心 | 基于光储直柔配电系统的节能用电评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104392286B (zh) | 2017-07-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104392286B (zh) | 考虑冷热电联供和储能运行策略的微电网运行优化方法 | |
Luo et al. | Multi-objective capacity optimization of a distributed energy system considering economy, environment and energy | |
CN113344736B (zh) | 一种园区级综合能源系统及其控制方法 | |
Pan et al. | Multi-objective and two-stage optimization study of integrated energy systems considering P2G and integrated demand responses | |
Shen et al. | Multi-objective capacity configuration optimization of an integrated energy system considering economy and environment with harvest heat | |
CN109784569A (zh) | 一种区域综合能源系统优化控制方法 | |
CN106022503A (zh) | 面向具有耦合型电冷热需求的微电网容量规划方法 | |
CN103997062B (zh) | 一种微网能量管理控制方法 | |
CN112600253B (zh) | 基于用能效率最优的园区综合能源协同优化方法及设备 | |
CN109659927A (zh) | 一种考虑储能参与度的综合能源微网储能容量配置方法 | |
CN111463836A (zh) | 一种综合能源系统优化调度方法 | |
CN108197768A (zh) | 一种能源系统与管网布局联合优化方法 | |
CN110110897A (zh) | 一种考虑不同储能运行策略的综合能源系统优化方法 | |
CN111737884B (zh) | 一种含多种清洁能源微能源网多目标随机规划方法 | |
CN110110904A (zh) | 考虑经济性、独立性和碳排放的综合能源系统优化方法 | |
CN109510241A (zh) | 工业园区风光燃储能能源的并网模式优化配置系统及方法 | |
CN113344249B (zh) | 基于区块链的冷热电联供多微网优化调度方法和系统 | |
CN109685332A (zh) | 一种综合能源多主体利益均衡优化调度方法及设备 | |
CN112131712B (zh) | 客户侧多能源系统多目标优化方法和系统 | |
Zhang et al. | Energy scheduling optimization of the integrated energy system with ground source heat pumps | |
CN111969657A (zh) | 一种基于谈判博弈的多能源站系统经济优化运行方法 | |
Wu et al. | Optimal design method and benefits research for a regional integrated energy system | |
Li et al. | Role of different energy storage methods in decarbonizing urban distributed energy systems: A case study of thermal and electricity storage | |
Ren et al. | Life-cycle-based multi-objective optimal design and analysis of distributed multi-energy systems for data centers | |
CN111882452A (zh) | 综合能源系统参与需求侧响应的边际成本计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170721 Termination date: 20211202 |
|
CF01 | Termination of patent right due to non-payment of annual fee |