[go: up one dir, main page]

CN111817313B - Optical storage capacity optimal configuration method and system based on sub-band mixed energy storage - Google Patents

Optical storage capacity optimal configuration method and system based on sub-band mixed energy storage Download PDF

Info

Publication number
CN111817313B
CN111817313B CN202010672667.1A CN202010672667A CN111817313B CN 111817313 B CN111817313 B CN 111817313B CN 202010672667 A CN202010672667 A CN 202010672667A CN 111817313 B CN111817313 B CN 111817313B
Authority
CN
China
Prior art keywords
energy storage
storage device
power
capacity
frequency
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
CN202010672667.1A
Other languages
Chinese (zh)
Other versions
CN111817313A (en
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.)
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
State Grid Shandong Electric Power Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
State Grid Shandong Electric Power Co Ltd
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 Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd, State Grid Shandong Electric Power Co Ltd filed Critical Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Priority to CN202010672667.1A priority Critical patent/CN111817313B/en
Publication of CN111817313A publication Critical patent/CN111817313A/en
Application granted granted Critical
Publication of CN111817313B publication Critical patent/CN111817313B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for AC mains or AC distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for AC mains or AC distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for AC mains or AC distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/22The renewable source being solar energy
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E70/00Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30Systems combining energy storage with energy generation of non-fossil origin

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Abstract

一种基于分频段混合储能的光储容量优化配置方法及系统,方法包括:基于希尔伯特黄变换分解光伏出力数据,将其分解为高频分量和低频分量;随机生成储能装置容量作为变量;设置变量的上下限;将生成的初始储能装置容量输入到适应度函数中,计算适应度函数中的目标函数;对当前的种群进行遗传变异,从而形成下一代种群;检测遗传代数是否达到设置好的最大遗传代数,根据检测结果选择继续计算或者结束计算,将最后一代适应度最高的个体作为最终计算结果;利用高频组分和低频组分的储能优化结果分别确定超级电容和蓄电池的最优容量,有利于克服光伏分散性、能源密度低、间歇性的缺点,实现成本大幅降低。

Figure 202010672667

A method and system for optimal configuration of optical storage capacity based on sub-band hybrid energy storage, the method includes: decomposing photovoltaic output data based on Hilbert-Huang transform, and decomposing it into high-frequency components and low-frequency components; randomly generating the capacity of an energy storage device as a variable; set the upper and lower limits of the variable; input the generated initial energy storage device capacity into the fitness function to calculate the objective function in the fitness function; genetically mutate the current population to form the next generation population; detect the genetic algebra Whether the set maximum genetic algebra is reached, choose to continue the calculation or end the calculation according to the test results, and take the individual with the highest fitness in the last generation as the final calculation result; use the energy storage optimization results of the high-frequency components and low-frequency components to determine the supercapacitor respectively. And the optimal capacity of the battery is beneficial to overcome the shortcomings of photovoltaic dispersion, low energy density, and intermittent, and achieve a significant reduction in cost.

Figure 202010672667

Description

一种基于分频段混合储能的光储容量优化配置方法及系统A method and system for optimal configuration of optical storage capacity based on sub-band hybrid energy storage

技术领域technical field

本发明属于光伏发电领域,更具体地,涉及一种基于分频段混合储能的光储容量优化配置方法及系统。The invention belongs to the field of photovoltaic power generation, and more particularly, relates to a method and a system for optimizing the configuration of optical storage capacity based on sub-band hybrid energy storage.

背景技术Background technique

光伏发电是太阳能利用的最重要的方式之一,也是世界上能源开发利用中技术最成熟、最具有商业化前景的领域之一。太阳能天然存在着许多不足,例如具有波动性强、间歇性强、季节性变化显著等特征,导致光伏出力波动性强、稳定性低、难以控制等问题,因此发电量受到了限制。为光伏配置合适的储能是光伏长久发展的最有效的途径之一,其中超级电容已经得到了广泛而成功的应用。针对光伏发电的特点,对储能技术的要求应该是响应速度快、存储量大,综合考虑,本发明提出一种混合储能方式,与光伏出力配合运行以达到优化超短期并网计划的效果。Photovoltaic power generation is one of the most important ways to utilize solar energy, and it is also one of the fields with the most mature technology and the most commercial prospects in energy development and utilization in the world. Solar energy naturally has many shortcomings, such as strong volatility, strong intermittency, and significant seasonal changes, which lead to problems such as strong photovoltaic output volatility, low stability, and difficulty in control, so power generation is limited. Configuring suitable energy storage for photovoltaics is one of the most effective ways for the long-term development of photovoltaics, among which supercapacitors have been widely and successfully applied. According to the characteristics of photovoltaic power generation, the requirements for energy storage technology should be fast response speed and large storage capacity. Considering comprehensively, the present invention proposes a hybrid energy storage method, which cooperates with photovoltaic output to achieve the effect of optimizing the ultra-short-term grid connection plan. .

对于光伏-蓄能联合优化运行及容量配置,国内外已经有很多研究,但对于如何选择及设计混合储能方式的研究仍然较少。对于装机容量的优化设计涉及了很多因素,比如超级电容的地理位置、开发条件、单位投资等等。在配置超级电容器容量时,如果配置的容量较大,虽然可以更好的消除光伏的波动性的影响,但是相应的容量和土地的投资就更多,如果配置较小的容量,初始投资自然会降低,但是对于解决光伏的波动性的作用也小了。There have been many studies at home and abroad on the optimal operation and capacity allocation of photovoltaic-energy storage, but there are still few studies on how to select and design hybrid energy storage methods. The optimal design of installed capacity involves many factors, such as the geographical location of supercapacitors, development conditions, unit investment and so on. When configuring the supercapacitor capacity, if the configured capacity is larger, although the impact of photovoltaic fluctuations can be better eliminated, the corresponding capacity and land investment will be more. If the configuration is smaller, the initial investment will naturally be Reduced, but the effect on solving the volatility of photovoltaics is also small.

但对于光伏迅速发展的今天,配置最佳容量的超级电容器和蓄电池容量配置才是最合理的,不论是从建设投资还是电网运行方面都具有重要意义。However, for the rapid development of photovoltaics today, the configuration of supercapacitors and battery capacity with the best capacity is the most reasonable, which is of great significance in terms of construction investment and power grid operation.

发明内容SUMMARY OF THE INVENTION

为解决现有技术中存在的不足,本发明的目的在于,提供一种基于希尔伯特黄变换和遗传算法的光储容量配置方法,能够给出一定容量的光伏系统对应的储能装置的最佳容量,而且计算方法简单,实用性强。In order to solve the deficiencies in the prior art, the purpose of the present invention is to provide an optical storage capacity configuration method based on Hilbert-Huang transform and genetic algorithm, which can provide the energy storage device corresponding to a photovoltaic system with a certain capacity. The best capacity, and the calculation method is simple and practical.

本发明采用如下的技术方案。一种基于分频段混合储能的光储容量优化配置方法,其特征在于,包括以下步骤:The present invention adopts the following technical solutions. A method for optimizing the configuration of optical storage capacity based on frequency-frequency hybrid energy storage, characterized in that it includes the following steps:

步骤1,获取设定时间跨度内的历史光伏发电数据;Step 1: Obtain historical photovoltaic power generation data within a set time span;

步骤2,使用HHT分解处理步骤1获得的历史光伏发电数据,获得高频发电数据和低频发电数据;Step 2, using the HHT decomposition to process the historical photovoltaic power generation data obtained in step 1 to obtain high-frequency power generation data and low-frequency power generation data;

步骤3,初始化遗传算法参数,以超级电容容量和蓄电池容量为自变量,随机生成一个初始种群,种群总数为常数N,遗传代数为GEN,最大遗传代数为M,初始化GEN=0;Step 3, initialize the genetic algorithm parameters, take the supercapacitor capacity and the battery capacity as independent variables, randomly generate an initial population, the total population is a constant N, the genetic algebra is GEN, the maximum genetic algebra is M, and the initialization GEN=0;

步骤4,分别输入低频发电数据和高频发电数据,计算步骤3中种群个体的适应度函数,以储能装置生命周期总花费为第一个目标函数,以净功率大于零的个数为第二个目标函数;Step 4: Input the low-frequency power generation data and the high-frequency power generation data respectively, calculate the fitness function of the individual population in step 3, take the total life cycle cost of the energy storage device as the first objective function, and take the number of net power greater than zero as the first objective function. two objective functions;

步骤5,对当前的种群进行遗传变异形成下一代种群;利用当前种群中的不同个体的适应度函数值,对种群进行运算;对选择运算产生的可行种群进行交叉运算产生新的个体;把新个体带入到适应度函数中进行计算,并比较与上一代目标函数的优异;Step 5: Perform genetic variation on the current population to form the next generation population; use the fitness function values of different individuals in the current population to perform operations on the population; perform crossover operations on the feasible populations generated by the selection operation to generate new individuals; The individual is brought into the fitness function for calculation, and compared with the excellence of the previous generation objective function;

步骤6,判断遗传变异是否达到最大遗传代数,如果没有达到最大遗传代数,返回执行步骤5,如果达到最大遗传代数M,则计算结束,将最后一代种群适应度最高的作为最终结果,并输出最佳储能装置容量,其中输入低频发电数据的输出结果为最佳蓄电池容量,输入高频发电数据的输出结果为最佳超级电容容量;Step 6: Judge whether the genetic variation has reached the maximum genetic algebra. If the genetic variation has not reached the maximum genetic algebra, return to step 5. If the maximum genetic algebra M is reached, the calculation is over, and the last generation of population with the highest fitness is the final result, and output the maximum genetic algebra. The best energy storage device capacity, in which the output result of inputting low-frequency power generation data is the best battery capacity, and the output result of inputting high-frequency power generation data is the best supercapacitor capacity;

其中HHT是指希尔伯特黄变换。where HHT refers to the Hilbert-Huang transform.

优选地,步骤1中,所述时间跨度为一年,数据采样周期为15min。Preferably, in step 1, the time span is one year, and the data sampling period is 15 minutes.

优选地,步骤1中,获取设定时间跨度内的历史光伏发电功率Ppv和历史光伏计划功率Pplan;其中,Ppv={Ppv(i)}i=1,2,…,n,Pplan={Pplan(i)}i=1,2,…n,序号i对应表示历史时刻。Preferably, in step 1, the historical photovoltaic power generation power P pv and the historical photovoltaic planned power P plan within the set time span are obtained; wherein, P pv ={P pv (i)} i=1,2,...,n , P plan ={P plan (i)} i=1, 2, . . . n , and the sequence number i corresponds to a historical moment.

优选地,步骤2中,对步骤1获取设定时间跨度内的历史光伏发电功率Ppv和历史光伏计划功率Pplan数据使用HHT进行分解处理,将光伏发电数据分解为高频历史功率Phf_pv={Phf_pv(i)}i=1,2,…,n、低频历史功率Plf_pv={Plf_pv(i)}i=1,2,…,n、高频计划功率Phf_plan={Phf_plan(i)}i=1,2,…,n和低频计划功率Plf_plan={Plf_plan(i)}i=1,2,…,nPreferably, in step 2, the historical photovoltaic power generation power P pv and historical photovoltaic plan power P plan data obtained in step 1 within the set time span are decomposed using HHT, and the photovoltaic power generation data is decomposed into high-frequency historical power P hf_pv = {P hf_pv (i)} i=1,2,...,n , low frequency historical power P lf_pv ={P lf_pv (i)} i=1,2,...,n , high frequency plan power P hf_plan ={P hf_plan (i)} i=1,2,...,n and low frequency plan power P lf_plan ={P lf_plan (i)} i=1,2,...,n .

优选地,步骤2具体包括:Preferably, step 2 specifically includes:

步骤2.1,使用EMD分解历史光伏发电功率Ppv和历史光伏计划功率PplanStep 2.1, using EMD to decompose historical photovoltaic power P pv and historical photovoltaic plan power P plan ;

Figure BDA0002582903170000031
Figure BDA0002582903170000031

式中:where:

s(t)表示光伏发电历史数据,即解步骤1获取的Ppv和Pplans(t) represents the historical data of photovoltaic power generation, namely P pv and P plan obtained in step 1,

ck(t)表示IMF分量,即c1(t)表示高频分量,c2(t)表示低频分量,c k (t) represents the IMF component, that is, c 1 (t) represents the high-frequency component, c 2 (t) represents the low-frequency component,

r(t)表示残余函数;r(t) represents the residual function;

步骤2.2,使用HSA获得光伏出力历史数据时频谱;Step 2.2, use HSA to obtain the frequency spectrum of the historical PV output data;

Figure BDA0002582903170000032
Figure BDA0002582903170000032

Figure BDA0002582903170000033
Figure BDA0002582903170000033

Figure BDA0002582903170000034
Figure BDA0002582903170000034

式中:where:

Re表示取实部;Re means to take the real part;

ak(t)表示每一个IMF分量的瞬时振幅;a k (t) represents the instantaneous amplitude of each IMF component;

ω(t)表示瞬时频率;ω(t) represents the instantaneous frequency;

θ(t)表示瞬时相位,

Figure BDA0002582903170000035
θ(t) represents the instantaneous phase,
Figure BDA0002582903170000035

H(ω,t)表示瞬时振幅在时间、频率平面的分布;H(ω,t) represents the distribution of the instantaneous amplitude in the time and frequency planes;

H(ω)表示瞬时振幅在频率平面的分布;H(ω) represents the distribution of the instantaneous amplitude in the frequency plane;

优选地,其中,EMD是指经验模态分解;HSA是指希尔伯特谱分析,IMF是指固有模态函数。Preferably, wherein, EMD refers to empirical mode decomposition; HSA refers to Hilbert spectral analysis, and IMF refers to intrinsic mode function.

优选地,步骤3具体包括:Preferably, step 3 specifically includes:

定义自变量超级电容容量x1和蓄电池容量x2,x1_min≤x1≤x1_max,x1_min和x1_max分别为超级电容容量下限约束和上限约束,

Figure BDA0002582903170000041
和x2_max分别为蓄电池容量下限约束和上限约束;Define the independent variables supercapacitor capacity x 1 and battery capacity x 2 , x 1_min ≤x 1 ≤x 1_max , x 1_min and x 1_max are the lower limit and upper limit constraints of the super capacitor capacity, respectively,
Figure BDA0002582903170000041
and x 2_max are the lower limit and upper limit constraints of the battery capacity, respectively;

定义储能装置的充电功率Pp={Pp(i)i=1,2,…,n},Ppmin≤Pp(i)≤Ppmax,Ppmin和Ppmax分别为充电功率下限约束和上限约束; Define the charging power P p = { P p (i) i = 1, 2, . and upper bound;

定义储能装置的放电功率Ph={Ph(i)i=1,2,…,n},Phmin和Phmax分别为放电功率下限约束和上限约束;Define the discharge power of the energy storage device P h = {P h (i) i = 1, 2, ..., n }, and P hmin and P hmax are the lower and upper limit constraints of the discharge power, respectively;

定义荷电状态SOC={SOC(i)i=1,2,…,n}、荷电状态上限SOCmax、荷电状态下限SOCmin,SOC(i)=SOCmax时储能装置不能继续充电,SOC(i)=SOCmin时,储能装置不能继续放电,初始化SOC(1)=100%;Define the state of charge SOC={SOC(i) i=1,2,...,n }, the upper limit of the state of charge SOC max , the lower limit of the state of charge SOC min , the energy storage device cannot continue to charge when SOC(i)=SOC max , when SOC(i)=SOC min , the energy storage device cannot continue to discharge, and the initialization SOC(1)=100%;

定义单位容量储能装置投资成本。Define the investment cost of energy storage device per unit capacity.

优选地,步骤4中的适应度函数包括:运行策略计算、储能装置状态计算和目标函数计算。Preferably, the fitness function in step 4 includes: operation strategy calculation, energy storage device state calculation and objective function calculation.

优选地,运行策略计算包括:Preferably, running the policy calculation includes:

以如下公式计算各个时刻储能装置的Pextro(i)i=1,2,…n,其中Calculate the P extro (i) i=1,2,...n of the energy storage device at each moment with the following formula, where

Pextro(i)=Pplan(i)-Ppv(i),P extro (i)=P plan (i)-P pv (i),

如果Pextro(i)≤0,储能装置在i时刻处于充电状态,储能装置放电功率为零,即Ph(i)=0;If P extro (i)≤0, the energy storage device is in the charging state at time i, and the discharge power of the energy storage device is zero, that is, P h (i)=0;

如果Pextro(i)>0,储能装置在i时刻处于放电状态,储能装置放电功率为零,即Pp(i)=0。If P extro (i)>0, the energy storage device is in a discharge state at time i, and the discharge power of the energy storage device is zero, that is, P p (i)=0.

优选地,储能装置状态计算包括:以如下公式计算储能装置的荷电状态SOC(i)i=1,2,…,n,其中SOC(1)=100%,Preferably, the calculation of the state of the energy storage device includes: calculating the state of charge of the energy storage device SOC(i) i=1, 2, . . .

如果Pextro(i)≤0,储能装置在i时刻处于充电状态,If P extro (i)≤0, the energy storage device is in the charging state at time i,

如果Pextro(i)>0,表明储能装置在i时刻处于放电状态,If P extro (i)>0, it means that the energy storage device is in the discharge state at time i,

式中:where:

x表示自变量超级电容容量x1和蓄电池容量x2x represents the independent variable supercapacitor capacity x 1 and battery capacity x 2 ,

eta表示逆变效率,eta represents the inverter efficiency,

σ表示储能装置自放电损失。σ represents the self-discharge loss of the energy storage device.

优选地,储能装置状态计算包括:以如下公式计算充放电功率Pp(i)i=1,2,…,n和Ph(i)i=1,2,…,nPreferably, the calculation of the state of the energy storage device includes: calculating the charging and discharging power P p (i) i=1,2,...,n and P h (i) i=1,2,...,n by the following formulas,

如果Pextro(i)≤0,表明储能装置在i时刻处于充电状态,Ph(i)=0,以如下公式计算Pp(i),If P extro (i)≤0, it means that the energy storage device is in the charging state at time i, and P h (i)=0, and P p (i) is calculated by the following formula,

如果Pextro(i)>0,表明储能装置在i时刻处于放电状态,Pp(i)=0,以如下公式计算Ph(i),If P extro (i)>0, it means that the energy storage device is in the discharge state at time i, P p (i)=0, and P h (i) is calculated by the following formula,

式中:where:

x表示自变量超级电容容量x1和蓄电池容量x2x represents the independent variable supercapacitor capacity x 1 and battery capacity x 2 ,

eta表示逆变效率,eta represents the inverter efficiency,

σ表示储能装置自放电损失。σ represents the self-discharge loss of the energy storage device.

优选地,目标函数计算包括:计算储能装置生命周期总花费objective1和达到计划功率时刻的个数objective2;定义目标函数1:Preferably, the calculation of the objective function includes: calculating the total life cycle cost of the energy storage device objective 1 and the number of times when the planned power is reached objective 2 ; defining the objective function 1:

objective1=F1+F2+Fpenalty·nobjective 1 =F 1 +F 2 +F penalty n

式中:where:

F1表示超级电容在生命周期内的总花费,F 1 represents the total cost of the supercapacitor in its life cycle,

F2表示蓄电池在生命周期内的总花费,F 2 represents the total cost of the battery in the life cycle,

Fpenalty表示一年内光伏发电功率未达到计划值的惩罚费用,F penalty represents the penalty fee for photovoltaic power generation that does not reach the planned value within one year,

n表示储能装置生命周期;n represents the life cycle of the energy storage device;

定义目标函数2:Define objective function 2:

objective2=Reliabilityobjective 2 = Reliability

Reliability是净功率power_grid大于零的数量,即满足Pextro(i)≤0,且power_grid=Ppv(i)-Pplan(i)-Ph(i)>0的时刻i的数量。Reliability is the number of times that the net power power_grid is greater than zero, that is, the number of times i when P extro (i)≤0, and power_grid=P pv (i)-P plan (i)-P h (i)>0 is satisfied.

优选地,步骤4,将低频历史功率Plf_pv={Plf_pv(i)}i=1,2,…,n和低频计划功率Plf_plan={Plf_plan(i)}i=1,2,…,n作为低频组数据,将高频历史功率Phf_pv={Phf_pv(i)}i=1,2,…,n和高频计划功率Phf_plan={Phf_plan(i)}i=1,2,…,n作为高频组数据,分别输入适应度函数中进行计算,以储能装置生命周期总花费objective1为第一个目标函数,以净功率power_grid大于零的个数objective2为第二个目标函数。Preferably, in step 4, the low frequency historical power P lf_pv ={P lf_pv (i)} i=1,2,...,n and the low frequency plan power P lf_plan ={P lf_plan (i)} i=1,2,... ,n as the low frequency group data, the high frequency historical power P hf_pv ={P hf_pv (i)} i=1,2,...,n and the high frequency plan power P hf_plan ={P hf_plan (i)} i=1, 2 , . two objective functions.

优选地,步骤1-6中,以常规水电站替代蓄电池,以抽水蓄能电站替代超级电容。Preferably, in steps 1-6, the battery is replaced by a conventional hydropower station, and the super capacitor is replaced by a pumped storage power station.

本发明还提供了一种基于所述光储容量优化配置方法的基于分频段混合储能的光储容量优化配置系统,其特征在于,光储容量优化配置系统包括:The present invention also provides an optimal configuration system for optical storage capacity based on the hybrid energy storage in sub-frequency bands based on the method for optimal configuration of optical storage capacity, characterized in that the optimal configuration system for optical storage capacity includes:

数据采集模块,用于获取设定时间跨度内的历史光伏发电数据;The data acquisition module is used to obtain the historical photovoltaic power generation data within the set time span;

数据变换模块,内置HHT分解单元,将数据采集模块获得的历史光伏发电数据分解为高频发电数据和低频发电数据;The data conversion module has a built-in HHT decomposition unit, which decomposes the historical photovoltaic power generation data obtained by the data acquisition module into high-frequency power generation data and low-frequency power generation data;

优化模块,内置遗传算法单元,输入数据变换模块获得的高频发电数据和低频发电数据,随机生成一个初始种群,以超级电容容量和蓄电池容量为自变量,并获得最佳储能装置容量,The optimization module, with built-in genetic algorithm unit, inputs the high-frequency power generation data and low-frequency power generation data obtained by the data conversion module, randomly generates an initial population, takes the supercapacitor capacity and battery capacity as independent variables, and obtains the optimal energy storage device capacity,

显示输出模块,用于可视化显示数据采集模块、数据变换模块和优化模块使用和获得的数据。The display output module is used to visually display the data used and obtained by the data acquisition module, the data transformation module and the optimization module.

优选地,储能装置为超级电容和蓄电池,或者常规水电站和抽水蓄电站。Preferably, the energy storage devices are supercapacitors and batteries, or conventional hydropower stations and pumped storage power stations.

本发明的有益效果在于,与现有技术相比,本发明基于希尔伯特黄变换和遗传算法,通过将上一代的结果进行遗传变异,寻找最优解,最终得到最佳配置储能容量,从而更好地解决光伏发电波动的缺陷。在容量设计时采用分段处理,分别针对低频功率组分和高频功率组分进行储能设备容量的设计和发电计划的优化确定。The beneficial effect of the present invention is that, compared with the prior art, the present invention is based on Hilbert-Huang transform and genetic algorithm, by genetically mutating the results of the previous generation to find the optimal solution, and finally obtain the optimal configuration energy storage capacity , so as to better solve the defects of photovoltaic power generation fluctuations. In the capacity design, segmental processing is adopted, and the capacity design of the energy storage equipment and the optimization of the power generation plan are carried out for the low-frequency power components and the high-frequency power components respectively.

本发明基于光储联合系统经济性和稳定性进行分析计算,用于解决光伏大规模集中并网困难的问题,为储能系统在电力中起到更好地作用奠定了基础,有显著的社会价值和经济价值。The invention analyzes and calculates based on the economy and stability of the combined photovoltaic and storage system, is used to solve the problem of the difficulty of large-scale centralized photovoltaic grid connection, lays a foundation for the energy storage system to play a better role in electric power, and has significant social benefits. value and economic value.

附图说明Description of drawings

图1为本发明总体流程图;Fig. 1 is the overall flow chart of the present invention;

图2为适应度函数工作流程图;Fig. 2 is the work flow chart of fitness function;

图3为低稳定性联合系统输出功率与计划功率;Figure 3 shows the low stability combined system output power and planned power;

图4为高稳定性联合系统输出功率与计划功率;Figure 4 shows the output power and planned power of the high-stability combined system;

图5为高稳定性联合系统与低温定性联合系统储能装置充放电功率对比图。Figure 5 is a comparison chart of the charging and discharging power of the energy storage device of the high-stability combined system and the low-temperature qualitative combined system.

具体实施方式Detailed ways

下面结合附图对本申请作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本申请的保护范围。The present application will be further described below with reference to the accompanying drawings. The following examples are only used to more clearly illustrate the technical solutions of the present invention, and cannot be used to limit the protection scope of the present application.

如图1所示,本发明提供了一种基于分频段混合储能的光储容量优化配置方法,包括以下步骤:As shown in FIG. 1 , the present invention provides an optimal configuration method for optical storage capacity based on sub-band hybrid energy storage, including the following steps:

步骤1,获取设定时间跨度内的历史光伏发电功率Ppv和历史光伏计划功率Pplan;其中,Ppv={Ppv(i)}i=1,2,…,n,Pplan={Pplan(i)}i=1,2,…n,序号i对应表示历史时刻。Step 1: Obtain the historical photovoltaic power generation power P pv and the historical photovoltaic planned power P plan within the set time span; wherein, P pv ={P pv (i)} i = 1,2,...,n , P plan ={ P plan (i)} i=1, 2, . . . n , and the serial number i corresponds to a historical moment.

所属领域技术人员可以任意设定时间跨度和数据采样的间隔时间,作为一种非限制性的优选,数据时间间隔15分钟,数据时间跨度至少为1年,采用一年以上的运行数据进行时频变换可以更科学直观地捕捉低频和高频的变化规律。Those skilled in the art can arbitrarily set the time span and the data sampling interval. As a non-limiting preference, the data time interval is 15 minutes, the data time span is at least 1 year, and the time-frequency data is carried out using the operation data of more than one year. The transformation can capture the variation law of low frequency and high frequency more scientifically and intuitively.

步骤2,对获取设定时间跨度内的历史光伏发电功率Ppv和历史光伏计划功率Pplan数据进行分解处理,将光伏发电数据分为高频与低频数据,高频出力数据上下波动较大,低频出力数据波动较为缓慢,对于低频出力,由于波动缓慢,对应的需要储能装置充放电的变换就没那么迅速,蓄电池就可以处理;相应的对于高频出力数据,所需的充放电变化较为明显,因此需要超级电容进行处理。或者使用常规水电站替代蓄电池处理低频出力数据,以抽水蓄电站代替超级电容处理高频出力数据。Step 2, decompose the acquired historical photovoltaic power generation power P pv and historical photovoltaic planned power P plan data within the set time span, and divide the photovoltaic power generation data into high-frequency and low-frequency data, and the high-frequency output data fluctuates greatly. The low-frequency output data fluctuates slowly. For low-frequency output, due to the slow fluctuation, the corresponding change of charge and discharge of the energy storage device is not so fast, and the battery can handle it; correspondingly, for the high-frequency output data, the required charge and discharge changes are relatively Obviously, supercapacitors are therefore needed for processing. Or use conventional hydropower stations instead of batteries to process low-frequency output data, and use pumped-storage power stations instead of super capacitors to process high-frequency output data.

具体地,使用EMD(Empirical Mode Decomposition,经验模态分解)分解步骤1获取的Ppv和Pplan,能使非平稳数据进行平稳化处理,然后使用HAS(Hilbert SpectrumAnalysis,希尔伯特谱分析)获得Hilbert时频谱图,即使用HHT处理步骤1获取的Ppv和Pplan,获得Specifically, using EMD (Empirical Mode Decomposition, Empirical Mode Decomposition) to decompose P pv and P plan obtained in step 1, the non-stationary data can be stabilized, and then HAS (Hilbert Spectrum Analysis, Hilbert Spectrum Analysis) Obtain the spectrogram when Hilbert, that is, use the HHT to process the P pv and P plan obtained in step 1, and obtain

高频历史功率Phf_pv,Phf_pv={Phf_pv(i)}i=1,2,…,nHigh-frequency historical power P hf_pv , P hf_pv ={P hf_pv (i)} i=1,2,...,n ,

低频历史功率Plf_pv,Plf_pv={Plf_pv(i)}i=1,2,…,nLow frequency historical power P lf_pv , P lf_pv ={P lf_pv (i)} i=1,2,...,n ,

高频计划功率Phf_plan,Phf_plan={Phf_plan(i)}i=1,2,…,nHigh frequency plan power P hf_plan , P hf_plan ={P hf_plan (i)} i=1,2,...,n ,

低频频计划功率Plf_plan,Plf_plan={Plf_plan(i)}i=1,2,…,nLow frequency plan power P lf_plan , P lf_plan ={P lf_plan (i)} i=1,2,...,n .

更具体地,使用如下的HHT公式处理步骤1获取的Ppv和PplanMore specifically, the P pv and P plan obtained in step 1 are processed using the following HHT formula:

Figure BDA0002582903170000081
Figure BDA0002582903170000081

Figure BDA0002582903170000082
Figure BDA0002582903170000082

Figure BDA0002582903170000083
Figure BDA0002582903170000083

Figure BDA0002582903170000084
Figure BDA0002582903170000084

式中:where:

s(t)表示光伏发电历史数据,即解步骤1获取的Ppv和Pplans(t) represents the historical data of photovoltaic power generation, namely P pv and P plan obtained in step 1;

ck(t)表示IMF(Intrinsic Mode Function,固有模态函数)分量,即c1(t)表示高频分量,c2(t)表示低频分量;分别包含了信号不同时刻特征尺度大小的成分。c k (t) represents the IMF (Intrinsic Mode Function, intrinsic mode function) component, that is, c 1 (t) represents the high-frequency component, and c 2 (t) represents the low-frequency component; they respectively contain the components of the characteristic scale at different times of the signal. .

r(t)表示残余函数,代表信号的平均趋势;r(t) represents the residual function, which represents the average trend of the signal;

Re表示取实部;Re means to take the real part;

ak(t)表示每一个IMF分量的振幅;a k (t) represents the amplitude of each IMF component;

ω(t)表示瞬时频率;ω(t) represents the instantaneous frequency;

θ(t)表示瞬时相位,

Figure BDA0002582903170000085
θ(t) represents the instantaneous phase,
Figure BDA0002582903170000085

H(ω,t)表示瞬时的振幅在时间、频率平面的分布;H(ω,t) represents the distribution of the instantaneous amplitude in the time and frequency planes;

H(ω)表示瞬时的振幅在频率平面的分布。H(ω) represents the distribution of the instantaneous amplitude in the frequency plane.

步骤3,初始化遗传算法参数,具体包括:Step 3: Initialize genetic algorithm parameters, including:

定义自变量超级电容容量x1和蓄电池容量x2,x1_min≤x1≤x1_max,x1_min和x1_max分别为超级电容容量下限约束和上限约束,

Figure BDA0002582903170000086
x2_min和x2_max分别为蓄电池容量下限约束和上限约束;Define the independent variables supercapacitor capacity x 1 and battery capacity x 2 , x 1_min ≤x 1 ≤x 1_max , x 1_min and x 1_max are the lower limit and upper limit constraints of the super capacitor capacity, respectively,
Figure BDA0002582903170000086
x 2_min and x 2_max are the lower limit and upper limit constraints of the battery capacity, respectively;

定义储能装置的充电功率Pp={Pp(i)i=1,2,…,n},Ppmin≤Pp(i)≤Ppmax,Ppmin和Ppmax分别为充电功率下限约束和上限约束; Define the charging power P p = { P p (i) i = 1, 2, . and upper bound;

定义储能装置的放电功率Ph={Ph(i)i=1,2,…,n},Phmin和Phmax分别为放电功率下限约束和上限约束;Define the discharge power of the energy storage device P h = {P h (i) i = 1, 2, ..., n }, and P hmin and P hmax are the lower and upper limit constraints of the discharge power, respectively;

定义荷电状态SOC={SOC(i)i=1,2,…,n}、荷电状态上限SOCmax、荷电状态下限SOCmin,SOC(i)=SOCmax时储能装置不能继续充电,SOC(i)=SOCmin时,储能装置不能继续放电,初始化SOC(1)=100%。Define the state of charge SOC={SOC(i) i=1,2,...,n }, the upper limit of the state of charge SOC max , the lower limit of the state of charge SOC min , the energy storage device cannot continue to charge when SOC(i)=SOC max , when SOC(i)=SOC min , the energy storage device cannot continue to discharge, and initializes SOC(1)=100%.

步骤4,分别输入高频发电数据和低频发电数据,计算步骤3中种群个体的适应度函数,以储能装置生命周期总花费为第一个目标函数,以净功率大于零的个数为第二个目标函数。Step 4: Input the high-frequency power generation data and the low-frequency power generation data respectively, calculate the fitness function of the individual population in step 3, take the total life cycle cost of the energy storage device as the first objective function, and take the number of net power greater than zero as the first objective function. two objective functions.

如图2所以,适应度函数中包括系统运行策略、稳定性计算、生命周期总花费计算,定义两个目标函数。As shown in Figure 2, the fitness function includes system operation strategy, stability calculation, and total life cycle cost calculation, and two objective functions are defined.

运行策略:输入光伏出力低频功率、光伏低频计划功率和储能装置容量,输出每个时刻的净功率、此储能装置容量下生命周期总花费和储能装置每个时刻的充放电功率。将没有达到计划功率时刻的净功率与高频功率合并在一起,通过高频数据的处理达到使低频发电功率完美匹配计划功率的结果。Operation strategy: input the low-frequency power of photovoltaic output, the low-frequency planned power of photovoltaics and the capacity of the energy storage device, and output the net power at each moment, the total life cycle cost under the capacity of the energy storage device, and the charging and discharging power of the energy storage device at each moment. The net power that has not reached the planned power moment is combined with the high-frequency power, and the low-frequency power generation power perfectly matches the planned power through the processing of high-frequency data.

本发明适应度函数具体步骤如下:The specific steps of the fitness function of the present invention are as follows:

1)输入种群个体,输入光伏功率和计划功率。1) Input the individual population, input the photovoltaic power and the planned power.

2)通过比较光伏功率与计划功率的大小计算需要的储能装置充放电功率。2) Calculate the required charging and discharging power of the energy storage device by comparing the photovoltaic power and the planned power.

3)将所需的充放电功率输入到电池模块中。具体地,在通过利用储能装置上一时刻的荷电状态进行计算,此时需要利用储能装置的容量约束,考虑上一时刻剩余电量、自损失和充放电损失,计算出下一时刻储能装置实际充放电功率和荷电状态。3) Input the required charging and discharging power into the battery module. Specifically, when calculating by using the state of charge of the energy storage device at the previous moment, it is necessary to use the capacity constraint of the energy storage device at this time, considering the remaining power, self-loss and charge-discharge loss at the previous moment, to calculate the storage capacity at the next moment. The actual charging and discharging power and state of charge of the device can be measured.

4)将上一步骤得出的结果输入到运行模块中,通过其他参数计算目标函数。4) Input the result obtained in the previous step into the running module, and calculate the objective function through other parameters.

5)下一组个体种群输入后,重复步骤1。5) After the next group of individual populations are input, repeat step 1.

本发明采用的遗传算法有益效果至少包括:在遗传算法中选择运算产生的种群名为可行种群,在本方法中由于设定了自变量取值范围,所以选择运算后不需要采用判定种群个体是否可行;在比较中采用两个目标函数都更优秀才能保留这个结果的方式。The beneficial effects of the genetic algorithm adopted in the present invention at least include: the population generated by the selection operation in the genetic algorithm is called the feasible population, and since the value range of the independent variable is set in this method, after the selection operation, it is not necessary to use to determine whether the individual population is Feasible; adopt a way in which both objective functions are better in comparison to preserve this result.

具体地,将低频历史功率Plf_pv={Plf_pv(i)}i=1,2,…,n和低频计划功率Plf_plan={Plf_plan(i)}i=1,2,…,n作为低频组数据,将高频历史功率Phf_pv={Phf_pv(i)}i=1,2,…,n和高频计划功率Phf_plan={Phf_plan(i)}i=1,2,…,n作为高频组数据,分别输入适应度函数中进行计算,以储能装置生命周期总花费objective1为第一个目标函数,以净功率power_grid大于零的个数objective2为第二个目标函数;Specifically, the low frequency historical power P lf_pv ={P lf_pv (i)} i=1,2,...,n and the low frequency plan power P lf_plan ={P lf_plan (i)} i=1,2,...,n are taken as For the low frequency group data, the high frequency historical power P hf_pv = {P hf_pv (i)} i=1,2,...,n and the high frequency plan power P hf_plan ={P hf_plan (i)} i=1,2,... ,n is used as high-frequency group data, respectively input into the fitness function for calculation, take the total life cycle cost of the energy storage device objective 1 as the first objective function, and take the number of net power power_grid greater than zero objective 2 as the second objective function;

以如下公式计算各个时刻储能装置的Pextro(i)i=1,2,…n,其中Calculate the P extro (i) i=1,2,...n of the energy storage device at each moment with the following formula, where

Pextro(i)=Pplan(i)-Ppv(i);P extro (i)=P plan (i)-P pv (i);

以如下公式计算储能装置的荷电状态SOC(i)i=1,2,…,n,其中SOC(1)=100%,Calculate the state of charge SOC(i) i=1,2,...,n of the energy storage device with the following formula, where SOC(1)=100%,

如果Pextro(i)≤0,表明储能装置在i时刻处于充电状态,If P extro (i) ≤ 0, it means that the energy storage device is in the charging state at time i,

Figure BDA0002582903170000101
Figure BDA0002582903170000101

式中:where:

eta表示逆变效率;eta represents the inverter efficiency;

σ表示储能装置自放电损失;σ represents the self-discharge loss of the energy storage device;

其中

Figure BDA0002582903170000102
表示如果储能装置把电量全部吸收,电量会溢出,所以没办法全部吸收电能。in
Figure BDA0002582903170000102
It means that if the energy storage device absorbs all the electricity, the electricity will overflow, so there is no way to absorb all the electricity.

如果Pextro(i)>0,表明储能装置在i时刻处于放电状态,If P extro (i)>0, it means that the energy storage device is in the discharge state at time i,

Figure BDA0002582903170000103
Figure BDA0002582903170000103

以如下公式计算充放电功率Pp(i)i=1,2,…,n和Ph(i)i=1,2,…,n Calculate the charging and discharging power P p (i) i=1,2,...,n and P h (i) i=1,2,...,n with the following formulas

如果Pextro(i)≤0,表明储能装置在i时刻处于充电状态,Ph(i)=0,以如下公式计算Pp(i),If P extro (i)≤0, it means that the energy storage device is in the charging state at time i, and P h (i)=0, and P p (i) is calculated by the following formula,

Figure BDA0002582903170000111
Figure BDA0002582903170000111

如果Pextro(i)>0,表明储能装置在i时刻处于放电状态,Pp(i)=0,以如下公式计算Ph(i),If P extro (i)>0, it means that the energy storage device is in the discharge state at time i, P p (i)=0, and P h (i) is calculated by the following formula,

Figure BDA0002582903170000112
Figure BDA0002582903170000112

计算储能装置生命周期总花费objective1和达到计划功率时刻的个数objective2;定义目标函数1:Calculate the total life cycle cost of the energy storage device objective 1 and the number of times the planned power is reached objective 2 ; define the objective function 1:

objective1=F1+F2+Fpenalty·nobjective 1 =F 1 +F 2 +F penalty n

式中:where:

F1表示超级电容在生命周期内的总花费,F 1 represents the total cost of the supercapacitor in its life cycle,

F2表示蓄电池在生命周期内的总花费,F 2 represents the total cost of the battery in the life cycle,

Fpenalty表示一年内光伏发电功率未达到计划值的惩罚费用,F penalty represents the penalty fee for photovoltaic power generation that does not reach the planned value within one year,

n表示储能装置生命周期;n represents the life cycle of the energy storage device;

定义目标函数2:Define objective function 2:

objective2=Reliabilityobjective 2 = Reliability

Reliability是净功率power_grid大于零的数量,即满足Pextro(i)≤0,且power_grid=Ppv(i)-Pplan(i)-Ph(i)>0的时刻i的数量。Reliability is the number of times that the net power power_grid is greater than zero, that is, the number of times i when P extro (i)≤0, and power_grid=P pv (i)-P plan (i)-P h (i)>0 is satisfied.

步骤5,对当前的种群进行遗传变异形成下一代种群;利用当前种群中的不同个体的适应度函数值,对种群进行运算;对选择运算产生的可行种群进行交叉运算产生新的个体;把新个体带入到适应度函数中进行计算,并比较与上一代目标函数的优异;Step 5: Perform genetic variation on the current population to form the next generation population; use the fitness function values of different individuals in the current population to perform operations on the population; perform crossover operations on the feasible populations generated by the selection operation to generate new individuals; The individual is brought into the fitness function for calculation, and compared with the excellence of the previous generation objective function;

步骤6,判断遗传变异是否达到最大遗传代数,如果没有达到最大遗传代数,返回执行步骤5,如果达到最大遗传代数M,则计算结束,将最后一代种群适应度最高的作为最终结果,并输出最佳储能装置容量,其中输入高频发电数据的输出结果为最佳超级电容容量,输入低频发电数据的输出结果为最佳蓄电池容量;Step 6: Judge whether the genetic variation has reached the maximum genetic algebra. If the genetic variation has not reached the maximum genetic algebra, return to step 5. If the maximum genetic algebra M is reached, the calculation is over, and the last generation of population with the highest fitness is the final result, and output the maximum genetic algebra. The optimal energy storage device capacity, in which the output result of inputting high-frequency power generation data is the best super capacitor capacity, and the output result of inputting low-frequency power generation data is the best battery capacity;

步骤6中,以常规水电站替代蓄电池,以抽水蓄能电站替代超级电容。In step 6, the battery is replaced by a conventional hydropower station, and the super capacitor is replaced by a pumped storage power station.

本发明还提供了一种基于所述光储容量优化配置方法的基于分频段混合储能的光储容量优化配置系统,光储容量优化配置系统包括:The present invention also provides an optimal configuration system for optical storage capacity based on the hybrid energy storage in sub-frequency bands based on the method for optimal configuration of optical storage capacity. The optimal configuration system for optical storage capacity includes:

数据采集模块,用于获取设定时间跨度内的历史光伏发电数据;The data acquisition module is used to obtain the historical photovoltaic power generation data within the set time span;

数据变换模块,内置HHT分解单元,将数据采集模块获得的历史光伏发电数据分解为高频发电数据和低频发电数据;The data conversion module has a built-in HHT decomposition unit, which decomposes the historical photovoltaic power generation data obtained by the data acquisition module into high-frequency power generation data and low-frequency power generation data;

优化模块,内置遗传算法单元,输入数据变换模块获得的高频发电数据和低频发电数据,随机生成一个初始种群,以超级电容容量和蓄电池容量为自变量,并获得最佳储能装置容量,The optimization module, with built-in genetic algorithm unit, inputs the high-frequency power generation data and low-frequency power generation data obtained by the data conversion module, randomly generates an initial population, takes the supercapacitor capacity and battery capacity as independent variables, and obtains the optimal energy storage device capacity,

显示输出模块,用于可视化显示数据采集模块、数据变换模块和优化模块使用和获得的数据。The display output module is used to visually display the data used and obtained by the data acquisition module, the data transformation module and the optimization module.

储能装置既可以是超级电容和蓄电池的组合,也可以是常规水电站和抽水蓄电站的组合。The energy storage device can be a combination of a super capacitor and a battery, or a combination of a conventional hydropower station and a pumped storage power station.

为更清楚地介绍本发明的技术方案和有益效果,以下介绍算例分析:In order to introduce the technical solutions and beneficial effects of the present invention more clearly, the following calculation example analysis is introduced:

输入某地区一年内每个15min的光伏出力数据,通过希尔伯特黄变换将数据分为高频与低频发电功率,设置种群个数为50,设置蓄电池的容量上限为3000kWh,设置超级电容器的容量上限为1000kWh,通过建立的数学模型和遗传算法,最终保留下来50个最优解。Enter the photovoltaic output data of each 15min in a year in a certain area, and divide the data into high-frequency and low-frequency power generation by Hilbert-Huang transformation. The upper limit of the capacity is 1000kWh. Through the established mathematical model and genetic algorithm, 50 optimal solutions are finally retained.

将最优解中的系统总花费最高和最低对应的两种储能装置的容量进行对比,其中总花费最高的对应蓄电池的容量是2872.52kWh,对应的超级电容器的容量为340.38kWh。总花费最低对应的蓄电池的容量是425.659kWh,对应的超级电容器的容量为890.9kWh。在运行策略板块中设置的是将光伏低频发电功率不满足计划功率的时刻所需的功率由高频部分来补充,因为在低频发电功率比高频发电功率要大的多,所以尽量使低频发电功率输入给电网的电更稳定。花费较高的联合系统30个小时输出的功率与计划功率对比如图3所示,花费较低的联合系统30个小时输出的功率与计划功率对比如图4所示,图5为花费较高的储能装置与花费较低的储能装置充放电对比图。The capacity of the two energy storage devices corresponding to the highest and lowest total system cost in the optimal solution is compared. The capacity of the corresponding battery with the highest total cost is 2872.52kWh, and the corresponding supercapacitor capacity is 340.38kWh. The battery capacity corresponding to the lowest total cost is 425.659kWh, and the corresponding supercapacitor capacity is 890.9kWh. In the operation strategy section, the power required by the low-frequency photovoltaic power generation does not meet the planned power to be supplemented by the high-frequency part. Because the low-frequency power generation is much larger than the high-frequency power generation, try to make the low-frequency power generation as much as possible. The power input to the grid is more stable. Figure 3 shows the power output of the combined system with higher cost and the planned power in 30 hours. Figure 4 shows the power output of the combined system with low cost and the planned power for 30 hours. Figure 5 shows the higher cost. The comparison chart of charging and discharging of energy storage device with lower cost energy storage device.

对于相同参数下的光储联合系统,使用混合能源的方法与不使用混合能源的系统总花费对比如下表所示,明显使用混合能源的系统总花费更低。For the combined solar-storage system with the same parameters, the total cost of the method using mixed energy is compared with that of the system that does not use mixed energy, as shown in the table below. Obviously, the total cost of the system using mixed energy is lower.

表1用混合能源的方法与不使用混合能源的系统总花费对比Table 1 Comparison of the total cost of the system with the method of mixed energy and the system without the use of mixed energy

Figure BDA0002582903170000131
Figure BDA0002582903170000131

本发明的有益效果在于,与现有技术相比,本发明基于希尔伯特黄变换和遗传算法,通过将上一代的结果进行遗传变异,寻找最优解,最终得到最佳配置储能容量,从而更好地解决光伏发电波动的缺陷。在容量设计时采用分段处理,分别针对低频功率组分和高频功率组分进行储能设备容量的设计和发电计划的优化确定。The beneficial effect of the present invention is that, compared with the prior art, the present invention is based on Hilbert-Huang transform and genetic algorithm, by genetically mutating the results of the previous generation to find the optimal solution, and finally obtain the optimal configuration energy storage capacity , so as to better solve the defects of photovoltaic power generation fluctuations. In the capacity design, segmental processing is adopted, and the capacity design of the energy storage equipment and the optimization of the power generation plan are carried out for the low-frequency power components and the high-frequency power components respectively.

本发明基于光储联合系统经济性和稳定性进行分析计算,用于解决光伏大规模集中并网困难的问题,为储能系统在电力中起到更好地作用奠定了基础,有显著的社会价值和经济价值。The invention analyzes and calculates based on the economy and stability of the combined photovoltaic and storage system, is used to solve the problem of the difficulty of large-scale centralized photovoltaic grid connection, lays a foundation for the energy storage system to play a better role in electric power, and has significant social benefits. value and economic value.

名词释义:Definition of noun:

HHT:Hilbert-Huang Transform,希尔伯特-黄变换;HHT: Hilbert-Huang Transform, Hilbert-Huang Transform;

EMD:Empirical Mode Decomposition,经验模态分解EMD: Empirical Mode Decomposition, empirical mode decomposition

IMF:Intrinsic Mode Function,固有模态函数;IMF: Intrinsic Mode Function, intrinsic mode function;

HSA:Hilbert Spectrum Analysis,希尔伯特谱分析。HSA: Hilbert Spectrum Analysis, Hilbert Spectrum Analysis.

本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述,但是本领域技术人员应该理解,以上实施示例仅为本发明的优选实施方案,详尽的说明只是为了帮助读者更好地理解本发明精神,而并非对本发明保护范围的限制,相反,任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。The applicant of the present invention has described and described the embodiments of the present invention in detail with reference to the accompanying drawings, but those skilled in the art should understand that the above embodiments are only preferred embodiments of the present invention, and the detailed description is only to help readers better It should be understood that the spirit of the present invention is not limited to the protection scope of the present invention. On the contrary, any improvement or modification made based on the spirit of the present invention should fall within the protection scope of the present invention.

Claims (11)

1. An optical storage capacity optimal configuration method based on sub-band mixed energy storage is characterized by comprising the following steps:
step 1, obtaining historical photovoltaic power generation power P in a set time spanpvAnd historical photovoltaic planned power Pplan(ii) a Wherein, Ppv={Ppv(i)}i=1,2,…,n,Pplan={Pplan(i)}i=1,2,…nThe serial number i correspondingly represents historical time;
step 2, obtaining the historical photovoltaic power generation power P in the set time span from the step 1pvAnd historical photovoltaic planned power PplanDecomposing the data by using Hilbert-Huang transform to decompose the photovoltaic power generation data into high-frequency historical power Phf_pv={Phf_pv(i)}i=1,2,…,nLow frequency historical power Plf_pv={Plf_pv(i)}i=1,2,…,nHigh frequency planned power Phf_plan={Phf_plan(i)}i=1,2,…,nAnd low frequency planned power Plf_plan={Plf_plan(i)}i=1,2,…,n(ii) a The method specifically comprises the following steps: step 2.1, decomposing historical photovoltaic power generation power P by using EMDpvAnd historical photovoltaic planned power Pplan
Figure FDA0003560812910000011
In the formula:
s (t) represents photovoltaic power generation historical data, namely P acquired in the step 1pvAnd Pplan
ck(t) denotes the IMF component, i.e. c1(t) represents a high frequency component, c2(t) represents a low-frequency component,
r (t) represents a residual function;
step 2.2, obtaining a photovoltaic output historical data time frequency spectrum by using HSA;
Figure FDA0003560812910000012
Figure FDA0003560812910000013
Figure FDA0003560812910000014
in the formula:
re represents a real part;
ak(t) represents the instantaneous amplitude of each IMF component;
ω (t) represents the instantaneous frequency;
theta (t) represents the instantaneous phase,
Figure FDA0003560812910000021
h (ω, t) represents the distribution of instantaneous amplitude in the time, frequency plane;
h (ω) represents the distribution of instantaneous amplitude in the frequency plane;
wherein, EMD refers to empirical mode decomposition; HSA refers to Hilbert spectrum analysis, IMF refers to intrinsic mode function;
step 3, initializing genetic algorithm parameters, randomly generating an initial population by taking the capacity of the super capacitor and the capacity of the storage battery as independent variables, wherein the total number of the population is a constant N, the genetic algebra is GEN, the maximum genetic algebra is M, and the initial GEN is 0; the method specifically comprises the following steps:
defining independent variable super capacitor capacity x1And battery capacity x2,x1_min≤x1≤x1_max,x1_minAnd x1_maxRespectively a super capacitor capacity lower limit constraint and an upper limit constraint, x2_min≤x2≤x2max,x2_minAnd x2_maxRespectively a lower limit constraint and an upper limit constraint of the capacity of the storage battery; defining the charging power P of an energy storage devicep={Pp(i)i=1,2,…,n},Ppmin≤Pp(i)≤Ppmax,PpminAnd PpmaxRespectively a charging power lower limit constraint and an upper limit constraint; defining the discharge power P of an energy storage deviceh={Ph(i)i=1,2,…,n},PhminAnd PhmaxRespectively a lower limit constraint and an upper limit constraint of the discharge power; define State of Charge SOC ═ { SOC (i)i=1,2,…,n}, upper limit of state of charge SOCmaxLower limit of state of charge SOCmin,SOC(i)=SOCmaxThe energy storage device can not be charged continuously, and SOC (i) ═ SOCminWhen the battery is charged, the energy storage device can not continue to discharge, and the initial SOC (1) is 100%; defining investment cost of an energy storage device with unit capacity;
step 4, respectively inputting low-frequency power generation data and high-frequency power generation data, calculating a fitness function of population individuals in the step 3, taking the total cost of the life cycle of the energy storage device as a first objective function, and taking the number of net powers greater than zero as a second objective function;
step 5, carrying out genetic variation on the current population to form a next generation population; calculating the population by using fitness function values of different individuals in the current population; performing cross operation on feasible populations generated by the selection operation to generate new individuals; the new individual is brought into the fitness function to be calculated, and the superiority of the fitness function is compared with that of the previous generation of objective function;
and 6, judging whether the genetic variation reaches the maximum genetic algebra, if not, returning to the step 5, if so, ending the calculation, taking the individual with the highest population fitness of the last generation as a final result, and outputting the optimal energy storage device capacity, wherein the output result of the input low-frequency power generation data is the optimal storage battery capacity, and the output result of the input high-frequency power generation data is the optimal super-capacitor capacity.
2. The optimal configuration method for optical storage capacity based on sub-band hybrid energy storage according to claim 1, wherein:
in step 1, the time span is one year, and the data sampling period is 15 min.
3. The optimal configuration method for optical storage capacity based on sub-band hybrid energy storage according to claim 1, wherein:
the fitness function in step 4 comprises: the method comprises the steps of operation strategy calculation, energy storage device state calculation and objective function calculation.
4. The optimal configuration method for optical storage capacity based on sub-band hybrid energy storage according to claim 2, wherein:
the operation strategy calculation comprises the following steps:
calculating P of the energy storage device at each moment by the following formulaextro(i)i=1,2,…nWherein
Pextro(i)=Pplan(i)-Ppv(i),
If P isextro(i) Less than or equal to 0, the energy storage device is in a charging state at the moment i, and the discharge power of the energy storage device is zero, namely Ph(i)=0;
If P isextro(i) When the voltage is higher than 0, the energy storage device is in a discharge state at the moment i, and the charging power of the energy storage device is zero, namely Pp(i)=0。
5. The optimal configuration method for optical storage capacity based on sub-band hybrid energy storage according to claim 3, wherein:
the energy storage device state calculation comprises: calculating the state of charge (SOC (i)) of the energy storage device according to the formulai=1,2,…,nWherein SOC (1) ═ 100%,
if P isextro(i) Less than or equal to 0, the energy storage device is in a charging state at the moment i,
Figure FDA0003560812910000031
if P isextro(i) Is greater than 0, indicating that the energy storage device is in a discharge state at time i,
Figure FDA0003560812910000041
in the formula:
x represents the independent variable super capacitor capacity x1And battery capacity x2
eta represents the efficiency of the inversion,
σ represents the energy storage device self-discharge loss.
6. The method of claim 4 or 5, wherein the method further comprises:
the energy storage device state calculation comprises: calculating the charge-discharge power P according to the following formulap(i)i=1,2,…,nAnd Ph(i)i=1,2,…,n
If P isextro(i) Less than or equal to 0, indicating that the energy storage device is in a charging state at time i, Ph(i) P is calculated as follows when P is 0p(i),
Figure FDA0003560812910000042
If P isextro(i) Greater than 0, indicating that the energy storage device is in a discharge state at time i, Pp(i) P is calculated as follows when P is 0h(i),
Figure FDA0003560812910000043
In the formula:
x represents the independent variable super capacitor capacity x1And battery capacity x2
eta represents the efficiency of the inversion,
σ represents the energy storage device self-discharge loss.
7. The optimal configuration method for optical storage capacity based on sub-band hybrid energy storage according to claim 6, wherein:
the objective function calculation includes: calculating total cost objective of life cycle of energy storage device1And the number of objective times to reach the planned power2(ii) a Defining an objective function 1:
objective1=F1+F2+Fpenalty·n
in the formula:
F1representing the total cost of the supercapacitor over the life cycle,
F2represents the total cost of the battery over the life cycle,
Fpenaltyrepresents the penalty cost of the photovoltaic power generation power not reaching the planned value within one year,
n represents the energy storage device lifecycle;
defining an objective function 2:
objective2=Reliability
reliabilitity is the amount of net power _ grid greater than zero, i.e., P is satisfiedextro(i) Not more than 0, and power _ grid ═ Ppv(i)-Pplan(i)-Ph(i) Number of time instants i > 0.
8. The optimal configuration method for optical storage capacity based on sub-band hybrid energy storage according to claim 7, wherein:
step 4, low-frequency historical power Plf_pv={Plf_pv(i)}i=1,2,…,nAnd low frequency planned power Plf_plan={Plf_plan(i)}i=1,2,…,nAs low frequency group data, high frequency historical power Phf_pv={Phf_pv(i)}i=1,2,…,nAnd high frequency planning power Phf_plan={Phf_plan(i)}i=1,2,…,nAs high-frequency group data, the high-frequency group data are respectively input into a fitness function for calculation so as to obtain the total cost objective of the life cycle of the energy storage device1As the first objective function, the number objective with net power _ grid larger than zero2Is the second objective function.
9. The optimal configuration method for optical storage capacity based on sub-band hybrid energy storage according to any one of claims 1 to 5, wherein:
in the steps 1-6, a conventional hydropower station is used for replacing a storage battery, and a pumped storage power station is used for replacing a super capacitor.
10. An optical storage capacity optimal configuration system based on sub-band hybrid energy storage based on the optical storage capacity optimal configuration method of any one of claims 1 to 9, wherein the optical storage capacity optimal configuration system comprises:
the data acquisition module is used for acquiring historical photovoltaic power generation data within a set time span;
the data transformation module is internally provided with an HHT decomposition unit and is used for decomposing the historical photovoltaic power generation data obtained by the data acquisition module into high-frequency power generation data and low-frequency power generation data;
an optimization module, a built-in genetic algorithm unit, which inputs the high-frequency power generation data and the low-frequency power generation data obtained by the data transformation module, randomly generates an initial population, takes the capacity of the super capacitor and the capacity of the storage battery as independent variables, and obtains the optimal capacity of the energy storage device,
and the display output module is used for visually displaying the data used and obtained by the data acquisition module, the data transformation module and the optimization module.
11. The system for optimized configuration of optical storage capacity based on sub-band hybrid energy storage according to claim 10,
the energy storage device is a super capacitor and a storage battery, or a conventional hydropower station and a water pumping and power storage station.
CN202010672667.1A 2020-07-14 2020-07-14 Optical storage capacity optimal configuration method and system based on sub-band mixed energy storage Active CN111817313B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010672667.1A CN111817313B (en) 2020-07-14 2020-07-14 Optical storage capacity optimal configuration method and system based on sub-band mixed energy storage

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010672667.1A CN111817313B (en) 2020-07-14 2020-07-14 Optical storage capacity optimal configuration method and system based on sub-band mixed energy storage

Publications (2)

Publication Number Publication Date
CN111817313A CN111817313A (en) 2020-10-23
CN111817313B true CN111817313B (en) 2022-05-06

Family

ID=72843398

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010672667.1A Active CN111817313B (en) 2020-07-14 2020-07-14 Optical storage capacity optimal configuration method and system based on sub-band mixed energy storage

Country Status (1)

Country Link
CN (1) CN111817313B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104795828A (en) * 2015-04-24 2015-07-22 清华大学 Wind storage capacity configuration method based on genetic algorithm
CN106707021A (en) * 2016-08-14 2017-05-24 华北电力大学(保定) Resonance detection method of distributed photovoltaic power generation cluster grid-connected system
CN109146131A (en) * 2018-07-11 2019-01-04 国网河北省电力有限公司 A kind of wind-power electricity generation prediction technique a few days ago
CN110797896A (en) * 2019-11-28 2020-02-14 国网江苏电力设计咨询有限公司 An optimal configuration method for hybrid energy storage capacity in a microgrid wind-solar storage system

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2805747A1 (en) * 2010-07-16 2012-01-19 The Trustees Of Columiba University In The City Of New York Machine learning for a power grid
CN106600022B (en) * 2015-10-20 2020-02-21 上海交通大学 Capacity Optimum Configuration Method of Wind-Solar-Gas-Sea Water Pumped Storage Isolated Power System Based on Multi-objective Optimization
CN109103941B (en) * 2018-10-26 2019-11-19 贵州电网有限责任公司 It is a kind of that method is stabilized based on the fluctuation of the wind-powered electricity generation of Hilbert-Huang transform and filtering algorithm
CN110676870B (en) * 2019-04-30 2021-09-21 国网新疆电力有限公司经济技术研究院 Hybrid energy storage capacity configuration method suitable for wind power grid connection

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104795828A (en) * 2015-04-24 2015-07-22 清华大学 Wind storage capacity configuration method based on genetic algorithm
CN106707021A (en) * 2016-08-14 2017-05-24 华北电力大学(保定) Resonance detection method of distributed photovoltaic power generation cluster grid-connected system
CN109146131A (en) * 2018-07-11 2019-01-04 国网河北省电力有限公司 A kind of wind-power electricity generation prediction technique a few days ago
CN110797896A (en) * 2019-11-28 2020-02-14 国网江苏电力设计咨询有限公司 An optimal configuration method for hybrid energy storage capacity in a microgrid wind-solar storage system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
独立海岛微网分布式电源容量优化设计;黎嘉明 等;《电工技术学报》;20160531;全文 *

Also Published As

Publication number Publication date
CN111817313A (en) 2020-10-23

Similar Documents

Publication Publication Date Title
CN104734166B (en) hybrid energy storage system and wind power generation power smooth control method
CN110797896A (en) An optimal configuration method for hybrid energy storage capacity in a microgrid wind-solar storage system
CN112736973A (en) Battery energy storage capacity configuration method and system for stabilizing output fluctuation of wind power and photovoltaic power station
Fang et al. Research on demand management of hybrid energy storage system in industrial park based on variational mode decomposition and Wigner–Ville distribution
CN103944177B (en) Hybrid energy-storing micro-grid system filter cutoff frequency defining method and output power control method
CN117081128A (en) Micro-grid hybrid energy storage capacity optimal configuration method based on variation modal decomposition
CN109787239B (en) Alternate iterative power flow calculation method and system for power flow model and multi-frequency AC hybrid system
CN105140942B (en) Consider the hybrid energy-storing optimization power distribution method of state-of-charge deviation
CN114201916A (en) A residual power decomposition method of microgrid hybrid energy storage system based on parameter optimization
CN104377724A (en) Coordinated optimization control method for improving economical efficiency of wind power/photovoltaic hybrid energy storage system
CN113644674B (en) Hydrogen hybrid energy storage capacity configuration system, method and application based on quantum particle swarm
CN107994593B (en) Capacity configuration method of composite energy storage system and terminal equipment
CN105720573A (en) Actually-measured data based modeling method for active power and reactive power control system of wind-light power storage station
CN111146793B (en) Photovoltaic-energy storage system capacity optimization design method and system based on power feature extraction
CN111313478A (en) Renewable energy storage optimization configuration method based on power smoothing
CN111817313B (en) Optical storage capacity optimal configuration method and system based on sub-band mixed energy storage
CN116911695B (en) Flexible resource adequacy evaluation method and device
CN113555901A (en) Hybrid energy storage capacity optimization method based on improved S-shaped function particle swarm optimization algorithm
CN118920432A (en) Light hydrogen storage isolated direct current micro-grid power balance control strategy
Das et al. Energy Management for a RES-Powered DC Microgrid Under Variable Load
CN117578957A (en) Cluster division method and system suitable for distributed light Fu Qunkong cluster modulation
Yuan et al. Research on the capacity configuration of the “flywheel+ lithium battery” hybrid energy storage system that assists the wind farm to perform a frequency modulation
Yang et al. Optimal allocation of energy storage capacity of high-permeability photovoltaic power generation system based on elastic neural network
CN110932312B (en) Reliability evaluation method for wind-solar storage micro-grid system
CN106849189A (en) Consider the micro-capacitance sensor economy and method for optimizing stability of regenerative resource randomness

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