基于低通滤波的混合储能系统的容量优化方法及装置

1.本技术涉及新能源技术领域，尤其涉及一种基于低通滤波的混合储能系统的容量优化方法及装置。


背景技术：

2.目前，在能源需求不断增加而旧能源逐渐消耗的形势下，发展新能源是解决能源危机的主要途径之一。其中，风力发电和光伏发电作为可再生能源的主力，在风资源和光资源丰富的地区得到了大规模的应用。由于容易受到天气等因素的影响，风、光能源具有随机性和不确定性的特点，在实际供电过程中可能会产生波动，而储能装置能够解决风光能源的功率波动问题，提高风光能源的渗透率，因此，风光储联合发电系统等混合储能系统的普及率逐渐提高。
3.然而，储能装置的成本通常较高，需要合理地配置储能设备，既要保证电网系统的正常运行，又要保证风光储联合发电系统的经济性。因此，如何合理的配置风光储联合发电系统的储能容量，兼顾联合系统的经济性与正常运行的可靠性成为目前亟需解决的问题。


技术实现要素：

4.本技术旨在至少在一定程度上解决相关技术中的技术问题之一。
5.为此，本技术的第一个目的在于提出一种基于低通滤波的混合储能系统的容量优化方法，该方法提出了对于储能容量设计效果的评价指标，定量的描述了风光储联合发电系统在设计和运行阶段系统成本及储能效率，可以合理配置混合储能系统的最优储能容量，提高系统的经济性、稳定性和安全性。
6.本技术的第二个目的在于提出一种基于低通滤波的混合储能系统的容量优化装置；
7.本技术的第三个目的在于提出一种非临时性计算机可读存储介质。
8.为达上述目的，本技术的第一方面实施例在于提出一种基于低通滤波的混合储能系统的容量优化方法，该方法包括以下步骤：
9.获取待优化的混合储能系统中的风电场站和光伏电站的原始样本数据；
10.通过快速傅里叶变换获取预设的多个滤波方式的频谱特性，根据所述频谱特性确定截止频率，并通过所述多个滤波方式对所述原始样本数据进行低通滤波以平抑数据波形；
11.基于优化需求、所述原始样本数据和所述混合储能系统的参数，确定遗传算法的适应度函数和约束条件；
12.基于所述适应度函数和所述约束条件，通过非支配排序遗传算法nsga-ii优化所述混合储能系统的储能容量，确定最优储能容量范围；
13.通过预设的评价指标依次评价配置了所述最优储能容量范围中不同储能容量的所述混合储能系统，确定所述混合储能系统的目标储能容量，并基于优化后的适应度参数
进行滤波方式的寻优。
14.可选地，在本技术的一个实施例中，适应度函数包括第一适应度函数和第二适应度函数，将全寿命周期成本lcc最小函数作为所述第一适应度函数，将波动率fr最小函数作为所述第二适应度函数。
15.可选地，在本技术的一个实施例中，通过以下公式计算所述第一适应度函数：
[0016][0017]
其中，
[0018]
其中，icc是混合储能系统的初始资本成本，n是混合储能项目的寿命，n是混合储能项目运行的年份，dn是年折旧，i是利率，tr是税率，an为年维护和运营成本，r是组件的替换次数，r是混合储能项目周期内的总替换数，iccc是待替换的组件的投资成本，lc是待替换的第c个部件的寿命，s是残值，s表示混合储能项目最后一年系统设备的回收价值，penalty是一年中未达到理想并网值的电量，m为惩罚电价，floor()函数表示将一个数字四舍五入到下一个小的整数。
[0019]
可选地，在本技术的一个实施例中，波动率表示预设时间内并网功率波动幅值与所述混合储能系统的额定功率之比，通过以下公式计算所述第二适应度函数：
[0020][0021]
其中，p'
max
为预设时间内理想并网功率和实际并网功率差值的最大值，p'
min
为预设时间内理想并网功率和实际并网功率差值的最小值，pn为混合储能系统的额定功率。
[0022]
可选地，在本技术的一个实施例中，约束条件包括：荷电状态约束条件、功率平衡约束条件和储能容量约束条件，所述荷电状态约束条件用于约束储能系统的荷电状态在允许的范围内，所述功率平衡约束条件用于约束任意时刻的并网功率与所述混合储能系统的出力相等，所述储能容量约束条件用于约束优化过程中的储能容量在所述储能系统的最小配置容量和最大配置容量之间。
[0023]
可选地，在本技术的一个实施例中，预设的评价指标包括稳定性指标和安全性指标，通过以下公式计算所述稳定性指标：
[0024]
stb＝rate_w+non_fit
[0025]
其中，
[0026]
其中，rate_w是弃风弃光率，non_fit是未贴合率，non_fit表示实际并网值与计算并网值的偏差，e为风电场站和光伏电站的总发电量，grid_ac是电网实际吸纳的电量，grid_aim为计划并网值。
[0027]
可选地，在本技术的一个实施例中，安全性指标表示通过不同的滤波方式进行低通滤波后所述混合储能系统的并网安全性，通过以下公式计算所述安全性指标：
[0028][0029]
其中，p
max
为预设时间内实际并网功率的最大值，p
min
为预设时间内实际并网功率的最小值，pn为混合储能系统的额定功率。
[0030]
可选地，在本技术的一个实施例中，基于优化后的适应度参数进行滤波方式的寻优，包括：
[0031]
通过以下公式计算归一化的波动率：
[0032][0033]
其中，f
min
为多个滤波方法中的波动率最小值，f
max
为多个滤波方法中的波动率最大值，f
t
为算例时间段内t时刻的波动率；
[0034]
结合所述归一化的波动率和所述全寿命周期成本确定不同范围内所述混合储能系统的最优滤波方式。
[0035]
为达上述目的，本技术的第二方面实施例还提出了一种基于低通滤波的混合储能系统的容量优化装置，包括以下模块:
[0036]
获取模块，用于获取待优化的混合储能系统中的风电场站和光伏电站的原始样本数据；
[0037]
滤波模块，用于通过快速傅里叶变换获取预设的多个滤波方式的频谱特性，根据所述频谱特性确定截止频率，并通过所述多个滤波方式对所述原始样本数据进行低通滤波以平抑数据波形；
[0038]
确定模块，用于基于优化需求、所述原始样本数据和所述混合储能系统的参数，确定遗传算法的适应度函数和约束条件；
[0039]
优化模块，用于基于所述适应度函数和所述约束条件，通过非支配排序遗传算法nsga-ii优化所述混合储能系统的储能容量，确定最优储能容量范围；
[0040]
评价模块，用于通过预设的评价指标依次评价配置了所述最优储能容量范围中不同储能容量的所述混合储能系统，确定所述混合储能系统的目标储能容量，并基于优化后的适应度参数进行滤波方式的寻优。
[0041]
为了实现上述实施例，本技术第三方面实施例还提出了一种非临时性计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述实施例中的基于低通滤波的混合储能系统的容量优化方法。
[0042]
本技术的实施例提供的技术方案至少带来以下有益效果：本技术先通过低通滤波处理功率数据，并基于nsga-ii进行储能容量的优化设计，确定最优储能容量范围。然后通过对于储能容量设计效果的评价指标，定量的描述风光储联合发电系统在设计和运行阶段的系统成本及储能效率，在最优储能容量范围内的不同储能容量中，确定最适合该混合储能系统的目标储能容量，提高了混合储能系统的储能容量确定的精确性和合理性。同时提出归一化波动率，可以为风光储联合能源系统的滤波方法寻优提供参考。由此，本技术可以节省容量优化的运算时间，合理配置混合储能系统的最优储能容量，提高风光储联合发电系统的经济性、稳定性和安全性，并可达成对于联合发电系统滤波方法的寻优。
[0043]
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
[0044]
本技术上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
[0045]
图1为本技术实施例提出的一种基于低通滤波的混合储能系统的容量优化方法的流程图；
[0046]
图2为本技术实施例提出的一种滤波后的样本数据与原始样本数据的频谱特性图；
[0047]
图3为本技术实施例提出的一种具体的基于非支配排序遗传算法的储能容量优化方法的流程图；
[0048]
图4为本技术实施例提出的一种混合储能系统的最优储能容量的评价指标示意图；
[0049]
图5为本技术实施例提出的一种结合归一化波动率和全寿命周期成本进行滤波方式寻优的原理示意图；
[0050]
图6为本技术实施例提出的一种基于低通滤波的混合储能系统的容量优化装置的结构示意图。
具体实施方式
[0051]
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
[0052]
需要说明的是，对于风光储系统，目前常用的储能容量优化设计，缺少对储能容量设计效果的评价，导致配置的储能容量的经济性较差或者无法保证联合系统稳定的出力。为此，本技术提出了基于低通滤波的混合储能系统的容量优化方法，该方法在风光储系统的设计和运行阶段通过评价指标来定量描述系统成本及储能效率，可以合理配置混合储能系统的最优储能容量。
[0053]
下面参考附图描述本发明实施例所提出的一种基于低通滤波的混合储能系统的容量优化方法和装置。
[0054]
图1为本技术实施例提出的一种基于低通滤波的混合储能系统的容量优化方法的流程图，如图1所示，该方法包括以下步骤：
[0055]
步骤s101：获取待优化的混合储能系统中的风电场站和光伏电站的原始样本数据。
[0056]
其中，混合储能系统是包括多种新能源发电系统和储能系统的联合系统，比如，混合储能系统可以是风光储联合发电系统，风光储是将风力发电系统、光伏发电系统和储能系统组成联合发电系统，风光储联合发电系统可以实现削峰填谷的功能，通过大容量的电池储能系统可以在较长的时间跨度内，改变联合系统出力特性。本技术中的混合储能系统也可以是还包括其他新能源发电系统的联合发电系统，此处不做限制。
[0057]
其中，原始样本数据包括风电场站和光伏电站的出力数据。
[0058]
具体实施时，本技术可以以一定的采样时间间隔在混合储能系统存储的历史数据中进行采样，获取一段时间内风电场站和光伏电站的出力数据。
[0059]
作为一种可能的实现方式，可以预先为混合储能系统配置数据采集与监视控制系统(supervisory control and data acquisition，简称scada)，通过scada可以对混合储能系统中运行现场的风电机组和光伏组件等设备进行监视和控制，实时采集风电场站和光伏电站的相关数据，并将采集到的数据存储在混合储能系统的数据库中，实现风电场和光伏电站的数据采集和数据存储功能。然后，在需要进行容量优化时，调用数据库中的历史数据，将1h作为采样时间间隔，选取一年的风电出力数据和光伏出力数据进行仿真。
[0060]
在本技术实施例中，还可以根据出力数据的时间间隔设置混合储能系统的充放电时间间隔，比如，在上述示例中以1h为出力数据的时间间隔时，则设置储能系统的充放电时间间隔为1h。
[0061]
102：通过快速傅里叶变换获取预设的多个滤波方式的频谱特性，根据频谱特性确定截止频率，并通过多个滤波方式对原始样本数据进行低通滤波以平抑数据波形。
[0062]
其中，本技术预设的滤波方式可以是各种低通滤波方式，比如，预设的多个滤波方式可以是巴特沃斯(butterworth)滤波和贝赛尔(bessel)滤波，通过butterworth滤波使通频带内的频率响应曲线最大限度平坦，减少纹波，通过bessel滤波在通频带上保持被过滤的信号波形。
[0063]
其中，截止频率是指输出信号能量开始大幅下降或者在滤波器中为大幅上升的边界频率，比如，当保持系统的输入信号的幅度不变，改变频率使输出信号降至最大值的0.707倍，用频响特性来表述-3db点处即为截止频率。
[0064]
具体的，通过预设的多个滤波方式对获取到的原始样本数据进行低通滤波，通过快速傅里叶变换获取频谱特性，根据频谱特性得到截止频率进行低通滤波的工作。
[0065]
举例而言，采用butterworth和bessel两种滤波方式对原始样本数据进行低通滤波，利用快速傅里叶变换获取频谱特性，采用butterworth和bessel两种滤波方式进行低通滤波后的频谱特性，与未滤波的原始样本数据的频谱特性如图2所示，根据该根据频谱特性图，可以确定截止频率，进行低通滤波的工作。通过比较图2中三种场景下的波形可知，本技术通过采用butterworth和bessel两种低通滤波可以平抑原始样本波形。
[0066]
步骤s103：基于优化需求、原始样本数据和混合储能系统的参数，确定遗传算法的适应度函数和约束条件。
[0067]
其中，适应度函数是描述个体性能的函数，用来对比个体与个体的相对优劣，个体的种群分级越靠前，则个体越优秀，即越符合适应度函数，个体越容易被选择，适应度的选取影响到后续进行优化的遗传算法的收敛速度以及能否找到最优解。
[0068]
其中，优化需求可以是容量优化过程中保证混合储能系统的相关性能且容量合理配置的需求。混合储能系统的参数可以是系统的成本参数和运行数据参数等相关参数。
[0069]
在本技术一个实施例中，适应度函数可以为多个，考虑保证混合储能系统的经济性和安全性的优化需求，将全寿命周期成本(life cycle cost，简称lcc)最小函数作为第一适应度函数，将波动率(fluctuation ratio，简称fr)最小函数作为第二适应度函数。
[0070]
具体而言，通过以下公式计算第一适应度函数：
[0071][0072]
其中，
[0073]
其中，icc是混合储能系统的初始资本成本，单位为元。n是混合储能项目的寿命，单位为年。n是混合储能项目运行的年份，指项目运行的第n年。dn是年折旧，单位为元。i是利率，tr是税率。an为年维护和运营成本，单位为元。r是组件的替换次数，指的是组件第r次替换。r是混合储能项目周期内的总替换数，iccc是待替换的组件的投资成本。lc是待替换的第c个部件的寿命，单位为年。s是残值，单位为元，s表示混合储能项目最后一年系统设备的回收价值。penalty是一年中未达到理想并网值的电量，单位为wh。m为惩罚电价，floor()函数表示将一个数字四舍五入到下一个小的整数。替换的总次数r是要替换的组件数量的寿命的函数，因此也是组件寿命的函数，可用上述公式表述。
[0074]
进一步的，计算波动率最小即第二个适应度函数。其中，波动率表示预设时间t内并网功率波动幅值与混合储能系统的额定功率之比，通过以下公式计算第二适应度函数：
[0075][0076]
其中，p'
max
为预设时间t内理想并网功率和实际并网功率差值的最大值，p'
min
为预设时间t内理想并网功率和实际并网功率差值的最小值，pn为混合储能系统的额定功率。
[0077]
需要说明的是，由上述公式可知，本技术在计算适应度函数时需要根据混合储能系统的成本参数和运行数据参数等各个参数进行计算，在获取上述参数时，可以通过读取混合储能系统的数据库中预存的数据，或者从步骤s101中获取原始样本数据中读取相关的历史运行数据等多种方式获取参数。作为一种示例，本技术在适应度函数计算过程中获取的关于成本的一些参数如下表1所示：
[0078]
经济性分析参数
[0079]
[0080][0081]
表1
[0082]
更进一步的，确定约束条件。在本技术一个实施例中，确定的约束条件包括：荷电状态约束条件、功率平衡约束条件和储能容量约束条件。其中，荷电状态约束条件用于约束储能系统的荷电状态在允许的范围内，功率平衡约束条件用于约束任意时刻的并网功率与混合储能系统的出力相等，储能容量约束条件用于约束优化过程中的储能容量在储能系统的最小配置容量和最大配置容量之间。
[0083]
具体而言，在本实施例中，对于荷电状态约束条件，由于过充过放会缩短储能设备寿命，增大储能设备更换次率，因此，为了延长储能电池的寿命，本技术通过荷电状态约束条件保证荷电状态在允许的范围内，荷电状态约束条件如以下公式所示：
[0084]
p
min
≤p(t)≤p
max
[0085]
soc
min
≤soc(t)≤soc
max
[0086]
其中，p(t)表示任意时刻储能系统的充放电功率，p
max
和p
min
分别表示储能的充放电功率的上限和下限，soc(t)表示任意采样时刻储能系统的荷电状态，soc
max
和soc
min
分别表示储能最小和最大荷电状态。
[0087]
对于功率平衡约束条件，通过该约束条件保证系统的有功功率平衡，即任意时刻的并网功率应与风光储出力相等，功率平衡约束条件如以下公式所示：
[0088]
p
grid
(t)＝p
pv
(t)+p
wind
(t)+p
ess
(t)
[0089]
其中，p
grid
(t)、p
pv
(t)、p
wind
(t)和p
ess
(t)分别为采样时刻的并网功率、光伏出力、风电出力和储能出力。
[0090]
对于储能容量约束条件，如以下公式所示：
[0091]cmin
≤c
ess
(t)≤c
max
[0092]
其中，c
ess
表示优化过程中的储能容量，c
max
和c
min
分别表示储能电池的最小和最大配置容量。
[0093]
步骤s104：基于适应度函数和约束条件，通过非支配排序遗传算法nsga-ii优化混合储能系统的储能容量，确定最优储能容量范围。
[0094]
其中，非支配排序遗传算法(non-dominated sorting genetic algorithm，简称nsga-ii)是一种多目标遗传算法，具有降低非劣排序遗传算法的复杂性、运行速度快和解集的收敛性好等特性。本技术通过非支配排序遗传算法nsga-ii优化风光储联合发电系统的储能容量，确定出较为符合上述优化需求的储能容量范围，即最优储能容量范围，便于后续在该范围内评选出最终优化出的目标储能容量。其中，通过非支配排序遗传算法nsga-ii进行优化时，相关步骤中通过步骤s103中确定的适应度函数和约束条件进行计算，以实现确定的最优储能容量范围符合优化需求。
[0095]
为了更加清楚的说明本技术基于非支配排序遗传算法进行储能容量优化的具体实现过程，在本技术一个实施例中还提出了一种基于非支配排序遗传算法的储能容量优化方法，如图3所示，该方法包括以下步骤：
[0096]
步骤s301，种群初始化。
[0097]
具体的，根据问题范围和步骤s103中确定的约束条件初始化种群。在本步骤中，设
置遗传算法的相关计算参数，生成相应规模的初始种群，便于后续根据该初始种群生成子代种群。举例而言，设定种群规模为50、交叉概率pe为0.5、变异概率pm为1.2和最大迭代次数设置为500。
[0098]
步骤s302，进行非支配排序。
[0099]
具体的，进行基于已初始化种群的非支配准则的排序过程。
[0100]
步骤s303，确定拥挤距离。
[0101]
具体的，当排序完成后，对拥挤距离值进行前向分配，便于后续种群中的个体根据等级和拥挤距离进行选择。
[0102]
步骤s304，选择个体。
[0103]
具体的，个体的选择采用带有步骤s303中确定的拥挤距离比较算子的二元锦标赛进行。在二元锦标赛中运用上述步骤确定的适应度函数选择个体。
[0104]
在本实施例中，通过二元锦标赛选择策略每次从种群中取出一定数量个体，然后选择其中最好的一个进入子代种群。重复该操作，直到新的种群规模达到原来的种群规模。具体包括以下步骤：第一步，确定每次选择的个体数量，比如，可以选择个体数量为2个。第二步，以相同的入选概率从种群中随机选择相应数量的个体构成组，再通过上述步骤中确定的适应度函数计算组中每个个体的适应度值，根据每个个体的适应度值，选择其中适应度值最好的个体进入子代种群。第三步，重复上述第二步，依次得到的个体构成新一代种群。
[0105]
步骤s305，通过遗传算子进行交叉和变异。
[0106]
具体的，实数编码ga采用模拟二进制进行交叉和多项式变异，交叉概率和变异概率可以为上述示例中设置的参数。由此，初始种群在进行非支配排序后，通过遗传算法的选择、交叉和变异三个基本操作得到第一代子代种群。
[0107]
步骤s306，进行重组与选择。
[0108]
具体的，将子代种群与当代种群相结合，下一代个体通过选择来设定。新一代随后由每一个前沿填充，直到人口规模超过当前人口规模。
[0109]
在本实施例中，从第二代子代种群开始，将父代种群与子代种群合并，进行快速非支配排序，同时对每个非支配层中的个体进行拥挤度计算，根据非支配关系以及个体的拥挤度选取合适的个体组成新的父代种群。然后，通过遗传算法的基本操作产生新的子代种群，按照上述方式进行迭代处理，直至满足结束条件，比如，人口规模超过当前人口规模，子代种群的代数超过预设的最大代数，即到达最大迭代次数500等。
[0110]
由此，利用适应度函数和约束条件通过非支配排序遗传算法nsga-ii优化混合储能系统的储能容量，求解满足联合系统的经济性和安全性下的如何设置储能容量的优化问题，根据该算法求解的结果确定最优储能容量范围。
[0111]
步骤s105：通过预设的评价指标依次评价配置了最优储能容量范围中不同储能容量的混合储能系统，确定混合储能系统的目标储能容量，并基于优化后的适应度参数进行滤波方式的寻优。
[0112]
其中，预设的评价指标是对储能容量的设计效果进行评价的指标，评价指标可以定量描述混合储能系统的成本、储能效率和安全性。
[0113]
具体的，在确定最优储能容量范围后，进行仿真计算，为混合储能系统依次配置该
范围内的各个储能容量，并获取配置不同储能容量下的运行数据来计算预设的评价指标，通过评价指标依次评价各个储能容量，确定出最适合当前混合储能系统的目标储能容量，目标储能容量可以是系统配置该储能容量后经济性和安全性等相关性能最高，具体可以通过比较不同储能容量对应的评价指标值的大小在最优储能容量范围中选取目标储能容量。
[0114]
在本技术一个实施例中，预设的评价指标包括稳定性指标和安全性指标，即提出稳定性和安全性两个指标评价联合发电系统的性能。具体实施时，首先，本技术实施例通过联合系统的弃风弃光率，以及实际并网值与计算并网值的偏差(即未贴合率)来描述系统的稳定性，通过以下公式计算稳定性指标：
[0115]
stb＝rate_w+non_fit
[0116]
其中，
[0117]
其中，rate_w是弃风弃光率，即能源浪费率，non_fit是未贴合率，non_fit表示实际并网值与计算并网值的偏差，e为风电场站和光伏电站的总发电量，grid_ac是电网实际吸纳的电量，grid_aim为计划并网值。
[0118]
然后，通过安全性指标表示通过不同的滤波方式进行低通滤波后混合储能系统的并网安全性，具体是通过以下公式计算安全性指标：
[0119][0120]
其中，p
max
为预设时间内实际并网功率的最大值，p
min
为预设时间内实际并网功率的最小值，pn为混合储能系统的额定功率。
[0121]
由此，根据本实施例提出的稳定性指标和安全性指标依次评价配置不同储能容量的联合发电系统的性能，确定出对当前混合储能系统优化后的最终的目标储能容量。
[0122]
进一步的，本技术还可以基于优化后的适应度参数进行滤波方式的寻优，即在上述预设的多个滤波方式寻找储能容量固定时，混合储能系统选择何种滤波方式具有更好的性能。其中，适应度参数与上述适应度函数相对应，比如，以波动率fr最小函数为适应度函数时，适应度参数为波动率。
[0123]
在本技术一个实施例中，优化后的适应度参数包括归一化波动率(normalized fr，以下简称n(f))，本技术结合归一化的波动率和全寿命周期成本确定混合储能系统的最优滤波方式。具体实施时，先通过以下公式计算归一化的波动率：
[0124][0125]
其中，f
min
为多个滤波方法中的波动率最小值，f
max
为多个滤波方法中的波动率最大值，f
t
为算例时间段内t时刻的波动率，算例时间段是本技术进行容量优化计算实例中预设的一个时间段，t是该算例时间段内的任一时刻，t为属于该时间段内的正数。
[0126]
然后，通过n(f)-lcc图寻优联合发电系统的最优滤波方式，具体可以是在归一化波动率n(f)相同的基础上，寻找联合发电系统的全寿命周期成本lcc更高的滤波方式。由此，在确定更优的滤波方式后，在后续可根据该滤波方式对混合储能系统中的样本数据进行滤波。
[0127]
综上所述，本技术实施例的基于低通滤波的混合储能系统的容量优化方法，先通过低通滤波处理功率数据，并基于nsga-ii进行储能容量的优化设计，确定最优储能容量范围。然后通过对于储能容量设计效果的评价指标，定量的描述风光储联合发电系统在设计和运行阶段的系统成本及储能效率，在最优储能容量范围内的不同储能容量中，确定最适合该混合储能系统的目标储能容量，提高了混合储能系统的储能容量确定的精确性和合理性。同时提出归一化波动率，可以为风光储联合能源系统的滤波方法寻优提供参考。由此，本技术的储能容量优化方法，可以节省容量优化的运算时间，合理配置混合储能系统的最优储能容量，提高风光储联合发电系统的经济性、稳定性和安全性，并可达成对于联合发电系统滤波方法的寻优。
[0128]
基于上述实施例，为了更加清楚的说明本技术的基于低通滤波的混合储能系统的容量优化方法具体实现过程和实现的技术效果，下面以一个的具体实施例进行详细说明。
[0129]
在本实施例中，以某350kw光伏电站与500kw风电厂为例，将1h作为采样时间间隔，选取一年的风电出力数据和光伏出力数据进行仿真。配置最优储能容量的联合发电系统的评价指标图如图4所示。在本实施例中设定基于butterworth滤波的场景为场景一(scenario 1，以下简称s1)，基于bessel滤波的场景命名为场景二(scenario 2，以下简称s2)。
[0130]
选择通过非支配排序遗传算法nsga-ii优化后确定的最优储能容量范围内的十个代表性储能容量，来评估两类滤波方法处理后的风光储联合系统的性能。其中，十个储能容量(energy storage capacity，esc)分别为esc1：0wh；esc2:4000wh；esc3:8000wh；esc4:15000wh；esc5:20000wh；esc6:25000wh；esc7:30000wh；esc8：35000wh；esc9:40000wh；esc10:45000wh。
[0131]
图4中展示了在十种esc下两种滤波处理后的的lcc、fr、non_fit、non_sft以及rate_w，即该图中横坐标为容量capacity，纵坐标为上述各参数parameter。其中，配备储能容量后，两场景的fr、non_fit、non_sft以及rate_w明显下降，表明配备储能后，联合系统的稳定性、安全性、波动性以及弃风弃光随着储能容量的增大而明显减小。然而随着储能容量的增加，lcc逐渐增高。因此，考虑到经济性，本技术实施例避免选择优化范围内的最大值作为目标储能容量，而根据评价指标选择目标储能容量。
[0132]
在两场景下为联合系统选择相同储能容量时，s1的lcc、fr、non_fit、以及rate_w均优于s2，但s1的并网安全性要低于s2，考虑不同的市场需求，可对于数据处理的方法进行优选，节省时间成本。
[0133]
本技术优化了fr，提出归一化波动率n(f)：
[0134][0135]
通过n(f)-lcc图寻优联合发电系统的最优滤波方式的原理如图5所示。
[0136]
根据s1和s2的n(f)-lcc曲线，将n(f)在0-1的范围内分为a-g 7个性能区域，分别为a区：0-0.59；b区：0.60-0.63；c区：0.64-0.69；d区：0.70-0.86；e区：0.87-0.90；f区：0.91-0.95；g区0.96-1。s1和s2在各性能区域的n(f)值相同，因此经济性较佳的滤波方法为最优滤波方法。
[0137]
a、c、e和g区中，s1的lcc低于s2，因此在a、c、e和g区，s1为适用场景，即达到同样稳
定性时，采用butterworth滤波时，联合系统有较好的的性能。同理，b、d和f区应用在s2有较好的性能。
[0138]
为了实现上述实施例，本技术还提出了一种基于低通滤波的混合储能系统的容量优化装置，图6为本技术实施例提出的一种基于低通滤波的混合储能系统的容量优化装置的结构示意图，如图6所示，该系统包括获取模块100、滤波模块200、确定模块300和优化模块400和评价模块500。
[0139]
其中，获取模块100，用于获取待优化的混合储能系统中的风电场站和光伏电站的原始样本数据。
[0140]
滤波模块200，用于通过快速傅里叶变换获取预设的多个滤波方式的频谱特性，根据频谱特性确定截止频率，并通过多个滤波方式对原始样本数据进行低通滤波以平抑数据波形。
[0141]
确定模块300，用于基于优化需求、原始样本数据和混合储能系统的参数，确定遗传算法的适应度函数和约束条件。
[0142]
优化模块400，用于基于适应度函数和约束条件，通过非支配排序遗传算法nsga-ii优化混合储能系统的储能容量，确定最优储能容量范围。
[0143]
评价模块500，用于通过预设的评价指标依次评价配置了最优储能容量范围中不同储能容量的混合储能系统，确定混合储能系统的目标储能容量，并基于优化后的适应度参数进行滤波方式的寻优。
[0144]
可选地，在本技术的一个实施例中，适应度函数包括第一适应度函数和第二适应度函数，将全寿命周期成本lcc最小函数作为所述第一适应度函数，将波动率fr最小函数作为第二适应度函数。
[0145]
可选地，在本技术的一个实施例中，确定模块300，具体用于：通过以下公式计算第一适应度函数：
[0146][0147]
其中，
[0148]
其中，icc是混合储能系统的初始资本成本，n是混合储能项目的寿命，n是混合储能项目运行的年份，dn是年折旧，i是利率，tr是税率，an为年维护和运营成本，r是组件的替换次数，r是混合储能项目周期内的总替换数，iccc是待替换的组件的投资成本，lc是待替换的第c个部件的寿命，s是残值，s表示混合储能项目最后一年系统设备的回收价值，penalty是一年中未达到理想并网值的电量，m为惩罚电价，floor()函数表示将一个数字四舍五入到下一个小的整数。
[0149]
可选地，在本技术的一个实施例中，波动率表示预设时间内并网功率波动幅值与所述混合储能系统的额定功率之比，确定模块300，具体用于：通过以下公式计算第二适应度函数：
[0150][0151]
其中，p'
max
为预设时间内理想并网功率和实际并网功率差值的最大值，p'
min
为预设时间内理想并网功率和实际并网功率差值的最小值，pn为混合储能系统的额定功率。
[0152]
可选地，在本技术的一个实施例中，约束条件包括：荷电状态约束条件、功率平衡约束条件和储能容量约束条件，荷电状态约束条件用于约束储能系统的荷电状态在允许的范围内，功率平衡约束条件用于约束任意时刻的并网功率与所述混合储能系统的出力相等，储能容量约束条件用于约束优化过程中的储能容量在储能系统的最小配置容量和最大配置容量之间。
[0153]
可选地，在本技术的一个实施例中，预设的评价指标包括稳定性指标和安全性指标，评价模块500，具体用于：通过以下公式计算稳定性指标：
[0154]
stb＝rate_w+non_fit
[0155]
其中，
[0156]
其中，rate_w是弃风弃光率，non_fit是未贴合率，non_fit表示实际并网值与计算并网值的偏差，e为风电场站和光伏电站的总发电量，grid_ac是电网实际吸纳的电量，grid_aim为计划并网值。
[0157]
可选地，在本技术的一个实施例中，安全性指标表示通过不同的滤波方式进行低通滤波后混合储能系统的并网安全性，评价模块500，具体用于：通过以下公式计算安全性指标：
[0158][0159]
其中，p
max
为预设时间内实际并网功率的最大值，p
min
为预设时间内实际并网功率的最小值，pn为混合储能系统的额定功率。
[0160]
可选地，在本技术的一个实施例中，评价模块500，还用于：通过以下公式计算归一化的波动率：
[0161][0162]
其中，f
min
为多个滤波方法中的波动率最小值，f
max
为多个滤波方法中的波动率最大值，f
t
为算例时间段内t时刻的波动率；
[0163]
结合归一化的波动率和全寿命周期成本确定不同范围内混合储能系统的最优滤波方式。
[0164]
需要说明的是，前述对基于低通滤波的混合储能系统的容量优化方法的实施例的解释说明也适用于该实施例的装置，此处不再赘述
[0165]
综上所述，本技术实施例的基于低通滤波的混合储能系统的容量优化装置，先通过低通滤波处理功率数据，并基于nsga-ii进行储能容量的优化设计，确定最优储能容量范围。然后通过对于储能容量设计效果的评价指标，定量的描述风光储联合发电系统在设计和运行阶段的系统成本及储能效率，在最优储能容量范围内的不同储能容量中，确定最适
合该混合储能系统的目标储能容量，提高了混合储能系统的储能容量确定的精确性和合理性。同时提出归一化波动率，可以为风光储联合能源系统的滤波方法寻优提供参考。
[0166]
为了实现上述实施例，本技术还提出了一种非临时性计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现如上述实施例中任一所述的基于低通滤波的混合储能系统的容量优化方法。
[0167]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本技术的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0168]
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本技术的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
[0169]
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本技术的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本技术的实施例所属技术领域的技术人员所理解。
[0170]
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(ram)，只读存储器(rom)，可擦除可编辑只读存储器(eprom或闪速存储器)，光纤装置，以及便携式光盘只读存储器(cdrom)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
[0171]
应当理解，本技术的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。如，如果用硬件来实现和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(pga)，现场可编程门阵列(fpga)等。
[0172]
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
[0173]
此外，在本技术各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
[0174]
上述提到的存储介质可以是只读存储器，磁盘或光盘等。尽管上面已经示出和描述了本技术的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本技术的限制，本领域的普通技术人员在本技术的范围内可以对上述实施例进行变化、修改、替换和变型。