CN112636367B - 一种平抑风电波动的混合储能容量优化配置方法 - Google Patents

一种平抑风电波动的混合储能容量优化配置方法 Download PDF

Info

Publication number
CN112636367B
CN112636367B CN202011432158.8A CN202011432158A CN112636367B CN 112636367 B CN112636367 B CN 112636367B CN 202011432158 A CN202011432158 A CN 202011432158A CN 112636367 B CN112636367 B CN 112636367B
Authority
CN
China
Prior art keywords
power
energy storage
super capacitor
hybrid energy
storage battery
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.)
Active
Application number
CN202011432158.8A
Other languages
English (en)
Other versions
CN112636367A (zh
Inventor
吴瑶
孟高军
孙玉坤
常青青
李建林
程哲明
刘海涛
夏涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Institute of Technology
Original Assignee
Nanjing Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Institute of Technology filed Critical Nanjing Institute of Technology
Priority to CN202011432158.8A priority Critical patent/CN112636367B/zh
Publication of CN112636367A publication Critical patent/CN112636367A/zh
Application granted granted Critical
Publication of CN112636367B publication Critical patent/CN112636367B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • 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/28—Arrangements for balancing of the load in a network by storage of energy
    • 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/28—Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • 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
    • 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
    • Y02E70/00—Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Abstract

本发明公开了一种平抑风电波动的混合储能容量优化配置方法，首先，通过滑动平均滤波获得混合储能系统参考功率PHESS，采用合成平均经验模态分解将PHESS分解为多个固有模态函数分量；随后，通过递归TF变换得到各IMF的瞬时频率‑时间曲线，在分界点临近分量处，以全寿命周期年化成本为目标函数构建HESS年均总成本优化模型，并结合智能IPQ算法求解最优频率分界点；最后，重构高、低频信号，考虑HESS充放电效率、换流器的能量转换效率及荷电状态来确定蓄电池和超级电容器的额定功率和容量。本发明所述平抑风电波动的混合储能容量优化配置方法，可以提高风储联合发电系统的电能质量，降低风力发电并网功率的波动，提高系统运行的稳定性和经济性。

Description

一种平抑风电波动的混合储能容量优化配置方法
技术领域
本发明涉及一种储能系统的容量配置方法，尤其涉及一种平抑风电波动的混合储能容量优化配置方法。
背景技术
近年来，为了转变以化石能源为主的能源格局，一场以推动可再生能源和新能源发展为核心的新能源革命已风起云涌。但是，风电的波动和间歇行为都具有强烈的不确定性，其对微电网的电力可靠性、电能质量、经济性的影响随着渗透率的增加而越发突出。采用储能技术能够提高风储联合发电系统的电能质量，降低风力发电并网功率的波动，提高系统运行的稳定性。随着储能技术的发展，目前广泛采用混合储能系统(Hybrid EnergyStorage System，HESS)来平抑风电功率波动，可分为能量型储能和功率型储能。
然而，在混合储能系统容量配置时，储能系统容量配置过大会增加微电网的投资成本，容量配置过小则会导致储能系统频繁切换工作状态，缩短电池的使用寿命。因此，如何合理配置HESS的容量，在满足系统平抑风功率波动要求的同时，降低建设成本成为HESS应用于风储联合发电系统的研究热点之一。
发明内容
发明目的：针对以上问题，本发明提出一种平抑风电波动的混合储能容量优化配置方法，能够提高风储联合发电系统的电能质量，降低风力发电并网功率的波动，提高系统运行的稳定性和经济性。
技术方案：本发明所采用的技术方案是一种平抑风电波动的混合储能容量优化配置方法，包括以下步骤：
(1)通过滑动平均滤波获得混合储能系统的参考功率，采用CAMD分解法将所述参考功率分解为多个IMFs分量，重复上述步骤若干次，将每次得到的IMF做合成平均处理后作为最终结果递归TF变换的输入。
其中，所述的通过滑动平均滤波获得混合储能系统的参考功率，是将风电输出功率与系统并网参考功率做差，得到所述参考功率，所述系统并网参考功率是风电输出功率经过滑动平均滤波处理后所得。
所述的CAMD分解法是将正态分布的白噪声信号叠加给所述参考功率，随后再进行经验模态分解，得到各IMF分量。
(2)通过递归TF变换得到各IMF分量的瞬时频率-时间曲线，在分界点临近分量处，以全寿命周期年化成本为目标函数，以最小容量、最小功率为约束条件，构建混合储能系统年均总成本优化模型，并结合智能IPQ算法求解最优频率分界点。
其中所述的递归TF变换，递归公式为：
Figure SMS_1
式中，ci(t)表示步骤(1)结束时所得的IMF分量中第i个模态分量，ci+1(t)为下一个第i+1个模态分量，T[ci(t)]为模态分量ci(t)的TF变换；
每经过一次递归计算相应的幅值函数ai(t)和相位函数φi(t)为：
Figure SMS_2
经所述递归TF变换后，瞬时频率-时间关系表示为：
Figure SMS_3
式中，φm(t)为递归TF变换计算到幅值函数趋于1时对应的相位函数。
步骤(2)中所述的混合储能系统年均总成本优化模型，目标函数为：
Csum＝(1+q1)αpbPB+[(1+q1)βeb+N′k1]EB+(1+q2)αpcPC+[(1+q2)βec+N′k2]EC
式中Csum表示全寿命周期的年化成本，PB、PC分别为蓄电池、超级电容器的额定功率；αpb、βeb分别代表了蓄电池的功率成本系数、容量成本系数；αpc、βec分别代表了超级电容器的功率成本系数、容量成本系数；q1、q2分别代表了蓄电池、超级电容器的更换次数；k1、k2分别为蓄电池、超级电容器的维护单价；N’为设备使用期数；EB为蓄电池的额定容量，EC为超级电容器的额定容量；
所述目标函数的约束条件为：
Figure SMS_4
Figure SMS_5
EB≥PBΔt
EC≥PCΔt
式中
Figure SMS_6
为需要平抑的最大功率，/>
Figure SMS_7
为需要平抑的风电功率容量最大值，Δt为采样时段。
采用智能IPQ算法求解所述的混合储能系统年均总成本优化模型。所述智能IPQ算法是将I算法所产生的解作为Q算法的初始点进行求解，其中I算法指的是以分界点临近分量作为变量的路径优化算法，Q算法指的是通过二次规划算法。
(3)根据最优频率分界点将高、低频信号重构，然后根据混合储能系统充放电效率、换流器的能量转换效率及荷电状态来确定蓄电池和超级电容器的额定功率和容量，采用功率型储能和能量型储能分别对风电机组的输出功率进行平抑。
其中所述的根据最优频率分界点将高、低频信号重构，是指将低于最优频率分界点的各个IMF重构，分配给蓄电池，高于最优频率分界点的各个IMF重构，分配给超级电容器。
所述重构的计算公式为：
Figure SMS_8
式中，PB(t)、PC(t)分别为蓄电池、超级电容器的额定功率，ci(t)表示步骤(1)结束时所得的IMF分量中第i个模态分量，rn(t)为分解筛除n个IMF分量后的信号残余分量，n为所有的IMF分量的个数，k为分配给超级电容器的IMF分量的个数。
确定蓄电池和超级电容器的额定功率和容量，包括以下步骤：
(31)确定超级电容器的实际补偿功率，计算式为：
Figure SMS_9
式中，ηDC-DC为换流器的能量转换效率，ηC，out为超级电容器的放电效率，PC(t)为超级电容器的功率，ηC,in为超级电容器的充电效率；
(32)确定蓄电池的充放电功率PESS(t)，计算式如下：
Figure SMS_10
式中，PB(t)为蓄电池的功率，ηDC-DC为换流器的能量转换效率，ηB,in和ηB,out为蓄电池的充放电效率；
(33)蓄电池的额定功率为：
Figure SMS_11
式中PB(t)为蓄电池的功率，ηDC-DC为换流器的能量转换效率，ηB，in和ηB，out为蓄电池的充放电效率；t0为起始时间，T为不平衡功率时间周期；
蓄电池的额定容量为：
Figure SMS_12
式中，SOC0为混合储能系统荷电状态的初始值；SOCup、SOClow分别为荷电状态的上、下限；ΔT为蓄电池功率指令间隔，i′为ΔT的个数。
超级电容器的额定功率大小为实际补偿功率
Figure SMS_13
绝对值的最大值；
超级电容器的额定容量为：
Figure SMS_14
式中，
Figure SMS_15
其中T为不平衡功率时间周期；SOCup、SOClow分别为荷电状态的上、下限。
有益效果：相比于现有技术，本发明具有以下优点：(1)蓄电池作为能量型储能，主要特征为能量密度大，功率密度小，可用来平抑低频部分；而超级电容器作为功率型储能元件能量密度较小，功率密度大，可进行高倍率充放电，负责平抑高频部分。混合储能系统充分利用蓄电池与超级电容器优势互补的特性，可增加储能系统的运行寿命，提高储能系统的整体性能。(2)本发明所述的CAMD分解法将白噪声扰动加入待分解的原始信号中，白噪声的频谱均匀分布，使得信号会自动分布到合适的参考尺度上。再对添加了白噪声的原始信号进行多次EMD分解，噪声经过多次的平均计算后会相互抵消，将多次分解所得到的IMF做合成平均处理后直接视作最终的IMF。再将合成平均处理后的各IMF经过递归TF变换，得到对应的各个IMF分量的瞬时频率-时间曲线，由于零均值噪声的特性，合成均值的计算结果与原始信号的差值随着集成平均的次数增加而减少，有效避免了EMD模态混叠的现象，提高HESS的经济性。经过递归TF变换计算得到的各个IMF分量的瞬时频率-时间曲线还可以观察到瞬时幅值、频率，并构建HESS年均总成本优化模型，结合智能IPQ算法求解最优频率分界点。
附图说明
图1是风储联合发电系统结构图；
图2是本发明所述平抑风电波动的混合储能容量优化配置方法的流程图。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步的说明。
在风电机组并网处安装储能装置不仅可以平滑功率波动，还可以获得良好的经济效益。为取得更好的平滑功率波动效果，降低储能系统投入成本，本专利提出了合成平均经验模态分解-递归TF变换的时频分析来优化配置储能系统的容量。本发明所述的平抑风电波动的混合储能容量优化配置方法具体步骤为：首先，通过滑动平均滤波获得HESS参考功率PHESS，采用CAMD(Compound AverageMode Decomposition，CAMD)分解法将参考功率PHESS分解为多个IMFs分量；随后，通过递归TF变换得到各IMF(intrinsic mode function，IMF)分量的瞬时频率-时间曲线，在分界点临近分量处，以全寿命周期年化成本为目标函数，以最小容量、最小功率为约束条件，构建HESS年均总成本优化模型，并结合智能IPQ算法求解最优频率分界点；最后，将高、低频信号重构后，采用功率型储能和能量型储能分别对其进行平抑，考虑HESS充放电效率、换流器的能量转换效率及荷电状态SOC(State Of Charge，SOC)来确定蓄电池和超级电容器的额定功率和容量。本发明所述平抑风电波动的混合储能容量优化配置方法的流程图如图2所示，具体技术方案如下：
(1)通过滑动平均滤波获得HESS参考功率PHESS，采用CAMD分解法将PHESS分解为多个IMFs分量。
图1为风储联合发电系统结构图，主要由风电机组、逆变器、混合储能系统(超级电容器和蓄电池)、直流母线组成。图中：PWG(t)为风力发电机组t时刻的输出功率，PHESS(t)为混合储能系统t时刻的输出功率，PG(t)为系统并网参考功率，PB(t)和PC(t)分别为蓄电池和超级电容器t时刻的充放电功率，放电时为正，充电时为负。为简化计算，本文暂不考虑负载对输出功率的影响。
将风力发电机组的输出功率通过为保证风电并网后电力系统的安全稳定运行、提升并网功率的电能质量，须对风电输出功率进行平滑处理后才能并入电网。本专利选择滑动平均法分解风电输出功率，以满足并网要求。我国2012年在《风电场接入电力系统技术规定》中划分了1min与10min两个时间尺度风电并网标准，如表1所示。将风电输出功率分为N个采样点，对非平稳的风电输出功率PwG(t)在适当的小区间上视为接近平稳的。将N个采样数据分别以窗口大小为S的最小区间内进行局部滑动平均滤波，得到较平滑的t时刻系统并网参考功率Pc(t)。
表1我国风电并网标准
Figure SMS_16
将获得的系统并网参考功率PG(t)与风电输出功率PWG(t)做差后得到HESS参考功率：
PG(t)＝PWG(t)-PHESS(t) (1)
PHESS(t)＝PB(t)+PC(t) (2)
当风电输出功率大于由滑动平均滤波得到的系统并网参考功率时，HESS充电；当风电输出功率小于由滑动平均滤波得到的系统并网参考功率时，HESS放电。
对HESS参考功率PHESS进行CAMD分解步骤如下：
(11)对HESS参考功率PHESS(t)叠加一组白噪声信号w(t)得到一个总体信号P(t)：
P(t)＝PHESS(t)+w(t) (3)
(12)对P(t)进行EMD分解，得到各IMF分量：
Figure SMS_17
式中，ci(t)为EMD分解得到的第i个IMF分量；rn(t)为分解筛除n个IMF分量后的信号残余分量，常常代表信号的直流分量或信号的趋势。
(13)将第j(j＝2，3，...，N)组白噪声信号叠加给HESS参考功率PHESS(t)，重复步骤(11)和(12)：
Figure SMS_18
式中，wj(t)为第j组相互独立的白噪声信号；Pj(t)为第j次加入白噪声信号后的不平衡功率信号；cji(t)为第j次EMD分解后得到的第i个IMF分量；rjn(t)为第j次分解得到的剩余部分。
(14)经过N次EMD分解，得到N组分解结果：
Figure SMS_19
(15)由于零均值噪声的特性，将多次叠加的白噪声信号对IMF分量的影响相互抵消，HESS参考功率PHESS(t)对应的IMF分量ci(t)及剩余分量rn(t)可表示为：
Figure SMS_20
(16)经过CAMD分解后重构的PHESS(t)可表示为：
Figure SMS_21
(2)通过递归TF变换得到各IMF分量的瞬时频率-时间曲线，在分界点临近分量处，以全寿命周期年化成本为目标函数，以最小容量、最小功率为约束条件，构建HESS年均总成本优化模型，并结合智能IPQ算法求解最优频率分界点。该步骤包括以下内容：
(21)HESS参考功率信号如式(8)所示，对每个模态分量ci(t)进行递归TF变换后为T[ci(t)]，则：
Figure SMS_22
构造各模态分量的解析信号zi(t)，即：
zi(t)＝ci(t)+jT[ci(t)]＝ai(t)exp[jφi(t)] (10)
Figure SMS_23
Figure SMS_24
式中，ai(t)为ci(t)的瞬间幅值；φi(t)为ci(t)的瞬时相位。
将c0(t)表示为瞬时幅a0(t)和纯调频信号cosφ0(t)的乘积。
c0(t)＝a0(t)cosφ0(t) (13)
以纯调频函数c1(t)＝cosφ0(t)作为新的信号，并对其进行TF变换，得到新的瞬时幅值和纯调频信号。不断进行递归计算，其递归公式为
Figure SMS_25
每经过一次递归计算便可得到相应的幅值函数和相位函数。
Figure SMS_26
不断进行递归计算直到幅值函数ai(t)趋向于1，则cm+1(t)＝cosφm(t)，综合上述递归计算步骤，信号ci(t)可表示为
Figure SMS_27
瞬时频率可表示为
Figure SMS_28
(22)在分界点临近分量的区间内，以全寿命周期年化成本为目标函数，以最小容量、最小功率为约束条件，构建HESS年均总成本优化模型：
在充分考虑HESS整个建设及运行过程中的经济效应，本专利以全寿命周期的年化成本Csum最小，即HESS的安装设备成本Cinvest、设备寿命终止产生的更换成本Cchange、运行维护成本Cmaintain之和最小为目标函数，计算公式如下：
Csum＝Cinvest+Cchange+Cmaintain (18)
Cinvest＝αpbPB+βebEB+αpcPC+βecEC (19)
Cchange＝q1(αpbPB+βebEB)+q2(αpcPC+βecEC) (20)
Cmaintain＝N′(k1EB+k2Ec) (21)
式中，Cinvest、Cchange、Cmaintain分别为HESS的安装设备成本、设备寿命终止产生的更换成本、运行维护成本；PB、PC分别为蓄电池、超级电容器的额定功率；αpb、βeb分别代表了蓄电池的功率成本系数、容量成本系数；αpc、βec分别代表了超级电容器的功率成本系数、容量成本系数；q1、q2分别代表了蓄电池、超级电容器的更换次数；k1、k2为蓄电池、超级电容器的维护单价；N’为设备使用期数。
将式(19)～式(21)代入式(18)中可得：
Csum＝(1+q1)αpbPB+[(1+q1)βeb+N′k1]EB+(1+q2)αpcPC+[(1+q2)βec+N′k2]EC (22)
对于目标函数，最小功率约束为：在任一时段t内，应满足PB+PC≥PHESS(t)，即要求蓄电池和超级电容器的配置功率之和不小于需要平抑的最大功率
Figure SMS_29
Figure SMS_30
同时，上述目标函数的最小容量约束为：对PHESS(t)按时段积分可得到需要平抑的风电功率容量，将该容量的最大值设为
Figure SMS_31
则HESS的最小容量至少要不低于/>
Figure SMS_32
考虑为了最大限制地平抑功率波动，HESS的初始SOC设定为0.5，即有：
Figure SMS_33
在采样时段Δt内，蓄电池和超级电容器的容量至少满足该时段内的充放电容量，即：
EB≥PBΔt (25)
EC≥PCΔt (26)
(23)采用智能IPQ算法求解最优分频频率，其主要原理是将I算法所产生的解作为Q算法的初始点进行求解，Q算法对I算法所产生的所有解进行微调，最终得到较好的解。其中I算法指的是变量(即分界点临近分量)根据环境以及自身属性选择一条较优的路径的优化算法。Q算法指的是通过二次规划来求解非线性问题。其算法步骤为：
(231)确定I数目NI，每条路径上的初始面积量、最大迭代次数、速度更新参数αv、bv和cv，面积量更新参数as、bs和cs，根据分界点临近分量的区间内设置参数的范围。
(232)设置I的初始速度以及携带的初始面积量，I访问过的节点集。对每个I的初始位置进行随机初始化。
(233)I根据路径中的面积量大小，按照式(27)概率大小选择下一个节点。I根据被选择的节点更新局部信息，按式(28)～(31)局部更新I的信息。
Figure SMS_34
式中：f(s)为面积量的函数，JI(i)为I在区域i可访问区域集合；s(i，j)表示在区域i时选择j作为下一个节点的面积量，Pi I(j)表示在区域i选择j作为下一个节点的概率。
s(I)＝s(I)+Δs(i，j) (28)
Figure SMS_35
s(i，j)＝(1-ρn)·s(i，j)-ρn·Δs(i，j) (30)
Figure SMS_36
式中：time(i，j)为时间函数，与I速度成反比，s(i，j)。
(234)重复步骤(233)，直到所有的I都完成移动，即所有的I都有一条路径形成完整的解。
(235)根据每个I所形成的移动路径，将其分别按式(32)、(33)计算出各个参数的值，得到本次迭代I的解。
Figure SMS_37
Pi＝Pi，min+(Pi，max-Pi，min)·xi/xi，max (33)
式中：Pi为参数xi的有效位数，Pi，min为有效位数的最小值，Pi，max为有效位数的最大值，mi为第i个参数的整数位个数，yj为第j个有效位的取值。
(236)将步骤(235)中本次迭代I所产生的解作为Q的初始点进行计算，得到本次迭代的解与目标值。
(237)确定本次迭代最优解并根据本次迭代最优解更新全局最优解，同时根据式(34)更新全局面积量。
Figure SMS_38
式中：n为参数有效位个数，
Figure SMS_39
为构造最优路径I携带的面积量，TIB指的是本次迭代的较优路径。
(238)当所有的I都收敛到同一路径上或者算法达到设定的最大迭代次数时转到步骤(239)，否则转到步骤(232)。
(239)将最优解与目标值保存，算法结束。
(3)将高、低频信号重构后，采用功率型储能和能量型储能分别对风机组的输出功率进行平抑，考虑HESS充放电效率、换流器的能量转换效率及SOC来确定蓄电池和超级电容器的额定功率和容量。该步骤包括以下内容：
经过步骤(2)的递归TF变化得到各IMF分量的瞬时频率-时间曲线，并在分界点临近分量区间内构建HESS年均总成本优化模型，结合智能IPQ算法求解出最优频率分界点ωk。将低于ωk的各个IMF分量重构，分配给蓄电池，高于ωk的各个IMF分量重构，分配给超级电容器，即：
Figure SMS_40
蓄电池的额定容量为EB，设混合储能系统荷电状态的初始值为SOC0，在ti时刻的荷电状态为：
Figure SMS_41
式中，ΔT为蓄电池功率指令间隔；SOCup、SOClow为荷电状态的上、下限；PESS(t)为计及换流器能量转换效率和充放电效率后蓄电池的充放电功率，其公式如下：
Figure SMS_42
式中，ηDC-DC为换流器的能量转换效率；ηB，in和ηB，out为蓄电池的充放电效率。
当研究时长为tiΔT时，根据式(36)可得蓄电池额定容量EB的计算公式为：
Figure SMS_43
其中ΔT为蓄电池功率指令间隔，ti为ΔT的个数。
在满足式(38)的前提下，还要满足下式(39)对SOC0的要求，蓄电池额定容量EB才能最小。
Figure SMS_44
蓄电池额定功率为PB，考虑变流器效率和蓄电池充放电效率：
Figure SMS_45
在考虑超级电容器的充电效率ηC,in和放电效率ηC，out的情况下，超级电容器的实际补偿功率式为：
Figure SMS_46
式中，PC(t)＞0为超级电容放电状态，PC(t)＜0为超级电容充电状态。
超级电容器的额定功率大小PC为实际补偿功率
Figure SMS_47
绝对值的最大值：
Figure SMS_48
由此，超级电容器的额定容量即可通过下面两个步骤求得：
①首先通过实际补偿功率
Figure SMS_49
求出超级电容器在T时段累积的能量ec(t)：
Figure SMS_50
式中，T为不平衡功率时间周期。
②获得累积能量后，额定容量可以通过下式确定：
Figure SMS_51
/>

Claims (9)

1.一种平抑风电波动的混合储能容量优化配置方法，其特征在于，包括以下步骤：
(1)通过滑动平均滤波获得混合储能系统的参考功率，采用CAMD分解法将所述参考功率分解为多个IMFs分量，重复上述步骤若干次，将每次得到的IMF做合成平均处理后作为最终结果递归TF变换的输入；
(2)通过递归TF变换得到各IMF分量的瞬时频率-时间曲线，在分界点临近分量处，以全寿命周期年化成本为目标函数，以最小容量、最小功率为约束条件，构建混合储能系统年均总成本优化模型，并结合智能IPQ算法求解最优频率分界点；采用智能IPQ算法求解最优分频频率的主要原理是将I算法所产生的解作为Q算法的初始点进行求解，使用Q算法对I算法所产生的所有解进行微调，得到最终解；其中I算法指的是以分界点临近分量作为变量，根据环境以及自身属性选择一条较优的路径的优化算法；Q算法指的是通过二次规划来求解非线性问题；
(3)根据最优频率分界点将高、低频信号重构，然后根据混合储能系统充放电效率、换流器的能量转换效率及荷电状态来确定蓄电池和超级电容器的额定功率和容量，采用功率型储能和能量型储能分别对风电机组的输出功率进行平抑。
2.根据权利要求1所述的平抑风电波动的混合储能容量优化配置方法，其特征在于：步骤(1)中所述的通过滑动平均滤波获得混合储能系统的参考功率，是将风电输出功率与系统并网参考功率做差，得到所述参考功率，所述系统并网参考功率是风电输出功率经过滑动平均滤波处理后所得。
3.根据权利要求1所述的平抑风电波动的混合储能容量优化配置方法，其特征在于：步骤(1)中所述的CAMD分解法是将正态分布的白噪声信号叠加给所述参考功率，随后再进行经验模态分解，得到各IMF分量。
4.根据权利要求1所述的平抑风电波动的混合储能容量优化配置方法，其特征在于：步骤(2)中所述的递归TF变换，递归公式为：
Figure FDA0004133964190000011
式中，ci(t)表示步骤(1)结束时所得的IMF分量中第i个模态分量，ci+1(t)为下一个第i+1个模态分量，T[ci(t)]为模态分量ci(t)的TF变换；
每经过一次递归计算相应的幅值函数ai(t)和相位函数φi(t)为：
Figure FDA0004133964190000012
5.根据权利要求4所述的平抑风电波动的混合储能容量优化配置方法，其特征在于：经所述递归TF变换后，瞬时频率-时间关系表示为：
Figure FDA0004133964190000021
式中，φm(t)为递归TF变换计算到幅值函数趋于1时对应的相位函数。
6.根据权利要求1所述的平抑风电波动的混合储能容量优化配置方法，其特征在于：步骤(2)中所述的混合储能系统年均总成本优化模型，目标函数为：
Csum＝(1+q1)αpbPB+[(1+q1)βeb+N′k1]EB+(1+q2)αpcPC+[(1+q2)βec+N′k2]EC
式中Csum表示全寿命周期的年化成本，PB、PC分别为蓄电池、超级电容器的额定功率；αpb、βeb分别代表了蓄电池的功率成本系数、容量成本系数；αpc、βec分别代表了超级电容器的功率成本系数、容量成本系数；q1、q2分别代表了蓄电池、超级电容器的更换次数；k1、k2分别为蓄电池、超级电容器的维护单价；N’为设备使用期数；EB为蓄电池的额定容量，EC为超级电容器的额定容量；
所述目标函数的约束条件为：
Figure FDA0004133964190000022
Figure FDA0004133964190000023
EB≥PBΔt
EC≥PCΔt
式中
Figure FDA0004133964190000024
为需要平抑的最大功率，/>
Figure FDA0004133964190000025
为需要平抑的风电功率容量最大值，△t为采样时段。
7.根据权利要求1所述的平抑风电波动的混合储能容量优化配置方法，其特征在于：步骤(3)中所述的根据最优频率分界点将高、低频信号重构，是指将低于最优频率分界点的各个IMF重构，分配给蓄电池，高于最优频率分界点的各个IMF重构，分配给超级电容器。
8.根据权利要求7所述的平抑风电波动的混合储能容量优化配置方法，其特征在于：所述重构的计算公式为：
Figure FDA0004133964190000026
式中，PB(t)、PC(t)分别为蓄电池、超级电容器的额定功率，ci(t)表示步骤(1)结束时所得的IMF分量中第i个模态分量，rn(t)为分解筛除n个IMF分量后的信号残余分量，n为所有的IMF分量的个数，k为分配给超级电容器的IMF分量的个数。
9.根据权利要求1所述的平抑风电波动的混合储能容量优化配置方法，其特征在于：步骤(3)中所述的根据混合储能系统充放电效率、换流器的能量转换效率及荷电状态来确定蓄电池和超级电容器的额定功率和容量，包括以下步骤：
(31)确定超级电容器的实际补偿功率，计算式为：
Figure FDA0004133964190000031
式中，ηDC-DC为换流器的能量转换效率，ηC,out为超级电容器的放电效率，PC(t)为超级电容器的功率，ηC,in为超级电容器的充电效率；
(32)确定蓄电池的充放电功率PEss(t)，计算式如下：
Figure FDA0004133964190000032
式中，PB(t)为蓄电池的功率，ηDC-DC为换流器的能量转换效率，ηB,in和ηB,out为蓄电池的充放电效率；
(33)蓄电池的额定功率为：
Figure FDA0004133964190000033
式中PB(t)为蓄电池的功率，ηDC-DC为换流器的能量转换效率，ηB,in和ηB,out为蓄电池的充放电效率；t0为起始时间，T为不平衡功率时间周期；
蓄电池的额定容量为：
Figure FDA0004133964190000034
式中，SOC0为混合储能系统荷电状态的初始值；SOCup、SOClow分别为荷电状态的上、下限；ΔT为蓄电池功率指令间隔，ti为ΔT的个数；
超级电容器的额定功率大小为实际补偿功率
Figure FDA0004133964190000035
绝对值的最大值；
超级电容器的额定容量为：
Figure FDA0004133964190000036
式中，
Figure FDA0004133964190000037
其中T为不平衡功率时间周期；SOCup、SOClow分别为荷电状态的上、下限。/>
CN202011432158.8A 2020-12-09 2020-12-09 一种平抑风电波动的混合储能容量优化配置方法 Active CN112636367B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011432158.8A CN112636367B (zh) 2020-12-09 2020-12-09 一种平抑风电波动的混合储能容量优化配置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011432158.8A CN112636367B (zh) 2020-12-09 2020-12-09 一种平抑风电波动的混合储能容量优化配置方法

Publications (2)

Publication Number Publication Date
CN112636367A CN112636367A (zh) 2021-04-09
CN112636367B true CN112636367B (zh) 2023-05-23

Family

ID=75309082

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011432158.8A Active CN112636367B (zh) 2020-12-09 2020-12-09 一种平抑风电波动的混合储能容量优化配置方法

Country Status (1)

Country Link
CN (1) CN112636367B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113162080A (zh) * 2021-04-23 2021-07-23 安徽信息工程学院 混合储能系统容量配置方法及系统
CN113162079A (zh) * 2021-04-23 2021-07-23 安徽信息工程学院 一种面向风电平抑的混合储能系统容量配置的方法及系统
CN113690912A (zh) * 2021-08-26 2021-11-23 南京理工大学 一种多类型储能装置分级协调控制方法及系统
CN113629737B (zh) * 2021-08-31 2023-06-27 国网新源控股有限公司 一种风光储蓄系统中化学储能的容量配置方法
CN114912848B (zh) * 2022-06-27 2023-05-26 南通大学 一种基于自适应滤波的全生命周期混合储能容量配置方法
CN116131291B (zh) * 2023-04-14 2023-08-29 武汉理工大学 用于平抑波动功率的氢电耦合系统容量配置方法
CN117713144B (zh) * 2024-02-06 2024-04-30 西安热工研究院有限公司 一种基于熔盐储能的火电机组的调频方法和系统

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5013372B2 (ja) * 2007-09-06 2012-08-29 国立大学法人 琉球大学 風力発電機用蓄電池設備の製作方法
CN104779630B (zh) * 2015-05-08 2017-07-11 武汉大学 一种平抑风电输出功率波动的混合储能系统容量配置方法
CN104933490B (zh) * 2015-07-03 2018-09-07 浙江大学 一种考虑社会车辆的土石坝工程的土石方调配优化方法
CN106998072A (zh) * 2017-05-15 2017-08-01 国网江苏省电力公司电力科学研究院 一种面向配电网优化运行的混合储能系统容量优化配置方法
CN109659970A (zh) * 2018-10-23 2019-04-19 国网甘肃省电力公司电力科学研究院 基于变分模态分解的储能系统平抑风功率波动的控制方法
CN110676870B (zh) * 2019-04-30 2021-09-21 国网新疆电力有限公司经济技术研究院 一种适用于风电并网的混合储能容量配置方法
CN110797896A (zh) * 2019-11-28 2020-02-14 国网江苏电力设计咨询有限公司 一种微网风光储系统中混合储能容量优化配置方法

Also Published As

Publication number Publication date
CN112636367A (zh) 2021-04-09

Similar Documents

Publication Publication Date Title
CN112636367B (zh) 一种平抑风电波动的混合储能容量优化配置方法
CN107104433B (zh) 一种光储系统参与配电网优化运行策略的获取方法
CN105896582B (zh) 一种微网储能容量优化配置方法
CN102684222B (zh) 一种基于储能技术的风力发电功率平滑控制方法
CN110797896A (zh) 一种微网风光储系统中混合储能容量优化配置方法
CN109659970A (zh) 基于变分模态分解的储能系统平抑风功率波动的控制方法
CN109617105B (zh) 一种基于下垂控制的分布式复合储能协同控制方法
CN107947211A (zh) 采用小波包分解并计及频率响应的孤岛型微电网储能优化配置方法
CN104659804A (zh) 含有混合储能的微电网及其控制方法
CN105260797A (zh) 一种计及经济性评估的微网储能电站规划控制方法
Yuming et al. Optimal capacity configuration of hybrid energy storage system considering smoothing wind power fluctuations and economy
CN112769149B (zh) 混联风光微电网快速频率响应分布式协调控制方法及系统
CN112086997A (zh) 一种基于变功率跟踪和超级电容储能的光伏协调调频控制方法
CN112260296A (zh) 一种用于平抑间歇性负荷的混合储能系统优化方法及装置
CN105140942A (zh) 考虑核电状态偏差的混合储能优化功率分配方法
CN112736973A (zh) 平抑风电、光伏电站出力波动的电池储能容量配置方法及系统
CN117081128A (zh) 基于变分模态分解的微电网混合储能容量优化配置方法
CN115811070A (zh) 辅助火电机组调频的飞轮储能自适应容量配置方法及系统
CN114201916A (zh) 一种基于参数优化的微网混合储能系统剩余功率分解方法
CN111654054A (zh) 一种储能基于自适应神经网络(ann)平抑短期风电波动的控制方法
CN113644674A (zh) 基于量子粒子群的氢混合储能容量配置系统、方法及应用
CN117458534A (zh) 一种新型液流储能调峰调频方法及装置
Wang et al. Research on Evaluation ofMulti-Timescale Flexibility and Energy Storage Deployment for the High-Penetration Renewable Energy of Power Systems.
CN104638634A (zh) 主从模式下基于带通滤波器的直流微电网振荡抑制方法
CN113746120A (zh) 基于ga的储能系统优化配置方法

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