CN111244991A - 一种考虑电池寿命的储能电站运行控制方法及系统 - Google Patents
一种考虑电池寿命的储能电站运行控制方法及系统 Download PDFInfo
- Publication number
- CN111244991A CN111244991A CN202010047125.5A CN202010047125A CN111244991A CN 111244991 A CN111244991 A CN 111244991A CN 202010047125 A CN202010047125 A CN 202010047125A CN 111244991 A CN111244991 A CN 111244991A
- Authority
- CN
- China
- Prior art keywords
- power
- energy storage
- station
- soc
- power station
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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/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
-
- 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)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑电池寿命的储能电站运行控制方法及系统,包括:获取含储能的快充电站区域的参数数据;根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率;根据所述含储能的快充电站区域的参数数据确定考虑电池寿命的含储能的快充电站区域储能电站的运行控制目标函数;通过遗传算法确定储能电站的充放电目标功率;以及通过预设的控制策略确定储能电站的最优充放电功率。本发明使得储能电站在参与含储能的快充电站区域削峰填谷控制和提高含储能的快充电站区域的经济收益的同时,兼顾电池使用寿命,使储能电池的放电深度控制在合理范围内,达到延缓电池寿命衰减的目的,提高了含储能的快充电站区域经济性。
Description
技术领域
本发明涉及储能与新能源领域,并且更具体地,涉及一种考虑电池寿命的含储能的快充电站区域储能电站运行控制方法及系统。
背景技术
含储能的快充电站区域主要指城市配电网或其他典型配电网中广域分散式接入的充(含快充或快/慢充混合)储电站区域、光储充(含快充或快/慢充混合)电站区域等。近年来城市中含储能的快充电站区域的数量和规模有不断扩大的趋势,随着我国城镇化水平的加速,含储能的快充电站区域也成为“智慧城市”的标志之一。含储能的快充电站区域作为城市的标志,其用能需求越来越大,负荷也呈现多类型、差异化和个性化的特点。现代含储能的快充电站区域中安装电池储能电站,有利于解决快充电站区域中负荷峰谷差大、能源利用率低等问题,随着能源互联网的逐步推进,未来的含储能的快充电站区域将是多种能源互联的智能化控制节点。
含储能的快充电站区域中的储能电站具有快速充放电的特性,储能电站在负荷低谷时段进行充电,在负荷高峰时段放电,从而缓解含储能的快充电站区域中峰时段用电压力并保障电网可靠运行,提高含储能的快充电站区域经济性。国内外已有储能系统或储能电站工程中,其控制策略专门针对含储能的快充电站区域的较少,含储能的快充电站区域作为城市能源互联网的重要组成部分,将成为一个小型的综合能源系统,储能系统作为其重要组成部分,其运行控制策略也尤为重要。
因此,需要一种对含储能的快充电站区域储能电站充放电进行有效控制的方法。
发明内容
本发明提出一种考虑电池寿命的储能电站运行控制方法及系统,以解决如何实现含储能的快充电站区域储能电站的经济运行的问题。
为了解决上述问题,根据本发明的一个方面,提供了一种考虑电池寿命的储能电站运行控制方法,其特征在于,所述方法包括:
获取含储能的快充电站区域的参数数据;所述参数数据包括:负荷、光伏、日前负荷预测功率、日前光伏预测功率、光伏装机容量和储能电站装机容量;
根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率;
根据所述含储能的快充电站区域的参数数据确定考虑电池寿命的含储能的快充电站区域储能电站的运行控制目标函数;
通过遗传算法确定储能电站的充放电目标功率;
通过预设的控制策略确定储能电站的最优充放电功率。
优选地,其中所述根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率,包括:
步骤1,根据所述日前负荷预测功率和日前光伏预测功率确定日前预测最大净负荷曲线的净负荷最大功率和净负荷最小功率;
步骤2,确定迭代步长;
步骤3,计算所述净负荷最大功率和迭代步长的差值为当前的削峰功率,计算所述净负荷最小功率和迭代步长的和作为当前的填谷功率;
步骤4,根据迭代步长反复调整削峰功率和填谷功率的值,直至统计的储能电站的总充电量和总放电量均大于储能电站的额定容量的修订值的差值;
步骤5,若当前的储能电站的总充电量与总放电量的差值的绝对值在预设范围内,则确定当前的削峰功率和填谷功率即为出储能电站的削峰功率和填谷功率;反之,返回步骤2,重新确定迭代步长。
优选地,其中所述运行控制目标函数为:
其中,F1为储能电站提高含储能的快充电站区域经济效益的评价函数;F2为提高储能电站寿命的评价函数;F3为SOC波动评价目标函数;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;p(t)为t时刻含储能的快充电站区域的商用电实时电价;Cprofit_max为储能电站SOC经历一个预设倍数的循环周期,且在谷时段充电、峰时段放电时的盈利值的最大值;Nmax[(DOD(i)]表示第i次放电深度对应的电池循环寿命(次数);SOCmax(T)为该时段内SOC最大值,SOCav(T)为T+Δt时段内SOC平均值。
优选地,其中所述方法还包括:
采用雨流计数法确定电池的放电深度。
优选地,其中在通过遗传算法确定储能电站的充放电目标功率时的约束条件,包括:
电网功率平衡约束为:
Pg=Pload-Ppv-PESS,
储能电站SOC及功率约束为:
SOCmin≤SOC(t)≤SOCmax,
0≤Pc≤Pc_e,
0≤Pd≤Pd_e,
并且,在储能电站充电过程中,其SOC还需满足:
在储能电站放电过程中,其SOC还需满足:
其中,Pg为电网向含储能的快充电站区域注入的功率;Ppv为含储能的快充电站区域的分布式光伏发电功率;Pload为含储能的快充电站区域的各类负荷的总功率;PESS为储能系统吸收或发出的功率;SOC(t)为t时刻储能电站电池的SOC值;SOCmin为储能电站电池SOC的最小约束值;SOCmax为储能电站电池SOC的最大约束值;Pc_e为储能电站的额定充电功率;Pd_e为储能电站的额定放电功率;γ为储能电站自放电率,SOC(t)为t时刻的储能电站的SOC状态值;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;Ee为储能电站额定容量。
优选地,其中所述通过预设的控制策略确定储能电站的最优充放电功率,包括:
当负荷高于削峰功率时,储能电站放电,放电功率为:
Pd=min(Pload-P1,Pd_obj),
当负荷低于填谷功率时,储能系统充电,充电功率为:
Pc=min(P2-Pload,Pc_obj),
其中,P1为削峰功率;Pd为放电功率;Pload为实际负荷功率;Pd_obj为通过遗传算法确定的储能电站的放电目标功率;P2为填谷功率;Pc为充电功率或放电功率;Pload为实际负荷功率;Pc_obj为基于考虑电池寿命的储能电站运行控制方法计算得出的充电功率值。
根据本发明的另一个方面,提供了一种考虑电池寿命的储能电站运行控制系统,其特征在于,所述系统包括:
参数数据获取单元,用于获取含储能的快充电站区域的参数数据;所述参数数据包括:负荷、光伏、日前负荷预测功率、日前光伏预测功率、光伏装机容量和储能电站装机容量;
削峰功率和填谷功率确定单元,用于根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率;
运行控制目标函数确定单元,用于根据所述含储能的快充电站区域的参数数据确定考虑电池寿命的含储能的快充电站区域储能电站的运行控制目标函数;
充放电目标功率确定单元,用于通过遗传算法确定储能电站的充放电目标功率;
最优充放电功率确定单元,用于通过预设的控制策略确定储能电站的最优充放电功率。
优选地,其中所述削峰功率和填谷功率确定单元,根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率,包括:
净负荷功率确定模块,用于根据所述日前负荷预测功率和日前光伏预测功率确定日前预测最大净负荷曲线的净负荷最大功率和净负荷最小功率;
迭代步长确定模块,用于确定迭代步长;
计算模块,用于计算所述净负荷最大功率和迭代步长的差值为当前的削峰功率,计算所述净负荷最小功率和迭代步长的和作为当前的填谷功率;
迭代模块,用于根据迭代步长反复调整削峰功率和填谷功率的值,直至统计的储能电站的总充电量和总放电量均大于储能电站的额定容量的修订值的差值;
削峰功率和填谷功率确定模块,若当前的储能电站的总充电量与总放电量的差值的绝对值在预设范围内,则确定当前的削峰功率和填谷功率即为出储能电站的削峰功率和填谷功率;反之,进入迭代步长确定模块,重新确定迭代步长。
优选地,其中所述运行控制目标函数为:
其中,F1为储能电站提高含储能的快充电站区域经济效益的评价函数;F2为提高储能电站寿命的评价函数;F3为SOC波动评价目标函数;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;p(t)为t时刻含储能的快充电站区域的商用电实时电价;Cprofit_max为储能电站SOC经历一个预设倍数的循环周期,且在谷时段充电、峰时段放电时的盈利值的最大值;Nmax[(DOD(i)]表示第i次放电深度对应的电池循环寿命(次数);SOCmax(T)为该时段内SOC最大值,SOCav(T)为T+Δt时段内SOC平均值。
优选地,其中所述运行控制目标函数确定单元,还包括:
放电深度计算模块,用于采用雨流计数法确定电池的放电深度。
优选地,其中在所述运行控制目标函数确定单元,通过遗传算法确定储能电站的充放电目标功率时的约束条件,包括:
电网功率平衡约束为:
Pg=Pload-Ppv-PESS,
储能电站SOC及功率约束为:
SOCmin≤SOC(t)≤SOCmax,
0≤Pc≤Pc_e,
0≤Pd≤Pd_e,
并且,在储能电站充电过程中,其SOC还需满足:
在储能电站放电过程中,其SOC还需满足:
其中,Pg为电网向含储能的快充电站区域注入的功率;Ppv为含储能的快充电站区域的分布式光伏发电功率;Pload为含储能的快充电站区域的各类负荷的总功率;PESS为储能系统吸收或发出的功率;SOC(t)为t时刻储能电站电池的SOC值;SOCmin为储能电站电池SOC的最小约束值;SOCmax为储能电站电池SOC的最大约束值;Pc_e为储能电站的额定充电功率;Pd_e为储能电站的额定放电功率;γ为储能电站自放电率,SOC(t)为t时刻的储能电站的SOC状态值;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;Ee为储能电站额定容量。
优选地,其中所述最优放电功率确定单元,通过预设的控制策略确定储能电站的最优充放电功率,包括:
当负荷高于削峰功率时,储能电站放电,放电功率为:
Pd=min(Pload-P1,Pd_obj),
当负荷低于填谷功率时,储能系统充电,充电功率为:
Pc=min(P2-Pload,Pc_obj),
其中,P1为削峰功率;Pd为放电功率;Pload为实际负荷功率;Pd_obj为通过遗传算法确定的储能电站的放电目标功率;P2为填谷功率;Pc为充电功率或放电功率;Pload为实际负荷功率;Pc_obj为基于考虑电池寿命的储能电站运行控制系统计算得出的充电功率值。
本发明提供了一种考虑电池寿命的储能电站运行控制方法及系统,包括:获取含储能的快充电站区域的参数数据;根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率;根据所述含储能的快充电站区域的参数数据确定考虑电池寿命的含储能的快充电站区域储能电站的运行控制目标函数;通过遗传算法确定储能电站的充放电目标功率;以及通过预设的控制策略确定储能电站的最优充放电功率。本发明的控制方法具有一定全局性和预见性,使得储能电站在参与含储能的快充电站区域削峰填谷控制和提高含储能的快充电站区域的经济收益的同时,兼顾电池使用寿命,使储能电池的放电深度控制在合理范围内,达到延缓电池寿命衰减的目的,本发明能够实现对含储能的快充电站区域储能电站充放电的有效控制,达到了延缓电池寿命衰减,提高了含储能的快充电站区域经济性。
附图说明
通过参考下面的附图,可以更为完整地理解本发明的示例性实施方式:
图1为根据本发明实施方式的考虑电池寿命的含储能的快充电站区域储能电站运行控制方法100的流程图;
图2为根据本发明实施方式的削峰功率和填谷功率的示意图;
图3为根据本发明实施方式的确定削峰功率和填谷功率的流程图;
图4为根据本发明实施方式的SOC变化曲线图;
图5为根据本发明实施方式的雨流技术法中雨滴下落的示意图;
图6为根据本发明实施方式的循环次数的示意图;
图7为根据本发明实施方式的利用遗传算法确定储能电站充放电功率最优个体的流程图;以及
图8为根据本发明实施方式的考虑电池寿命的含储能的快充电站区域储能电站运行控制系统800的结构示意图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式,然而,本发明可以用许多不同的形式来实施,并且不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中,相同的单元/元件使用相同的附图标记。
除非另有说明,此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外,可以理解的是,以通常使用的词典限定的术语,应当被理解为与其相关领域的语境具有一致的含义,而不应该被理解为理想化的或过于正式的意义。
图1为根据本发明实施方式的考虑电池寿命的含储能的快充电站区域储能电站运行控制方法100的流程图。如图1所示,本发明的实施方式提供的考虑电池寿命的含储能的快充电站区域储能电站运行控制方法具有一定全局性和预见性,使得储能电站在参与含储能的快充电站区域削峰填谷控制和提高含储能的快充电站区域的经济收益的同时,兼顾电池使用寿命,使储能电池的放电深度控制在合理范围内,达到延缓电池寿命衰减的目的,本发明能够实现对含储能的快充电站区域储能电站充放电的有效控制,达到了延缓电池寿命衰减,提高了含储能的快充电站区域经济性。本发明的实施方式提供的考虑电池寿命的含储能的快充电站区域储能电站运行控制方法100从步骤101处开始,在步骤101获取含储能的快充电站区域的参数数据;所述参数数据包括:负荷、光伏、日前负荷预测功率、日前光伏预测功率、光伏装机容量和储能电站装机容量。
在步骤102,根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率。
优选地,其中所述根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率,包括:
步骤1021,根据所述日前负荷预测功率和日前光伏预测功率确定日前预测最大净负荷曲线的净负荷最大功率和净负荷最小功率;
步骤1022,确定迭代步长;
步骤1023,计算所述净负荷最大功率和迭代步长的差值为当前的削峰功率,计算所述净负荷最小功率和迭代步长的和作为当前的填谷功率;
步骤1024,根据迭代步长反复调整削峰功率和填谷功率的值,直至统计的储能电站的总充电量和总放电量均大于储能电站的额定容量的修订值的差值;
步骤1025,若当前的储能电站的总充电量与总放电量的差值的绝对值在预设范围内,则确定当前的削峰功率和填谷功率即为出储能电站的削峰功率和填谷功率;反之,返回步骤1022,重新确定迭代步长。
本发明的实施方式的削峰功率和填谷功率的示意图如图2所示。
本发明的实施方式的计算削峰功率和填谷功率的流程如图3所示,包括:根据所述日前负荷预测功率Pload和日前光伏预测功率Ppv确定日前预测最大净负荷曲线的净负荷最大功率Pmax=max(Pload-Ppv)和净负荷最小功率Pmin=min(Pload-Ppv);确定迭代步长△P;计算当前的削峰功率为P1=Pmax-ΔP,计算当前的填谷功率为P2=Pmax+ΔP;根据迭代步长反复调整削峰功率和填谷功率的值,直至统计的储能电站的总充电量Ec和总放电量Ed均大于储能电站的额定容量的修订值ζEN的差值;若当前的储能电站的总充电量与总放电量的差值的绝对值在预设范围内,即|Ec-Ed|<ε,则确定当前的削峰功率和填谷功率即为出储能电站的削峰功率和填谷功率;反之,则重新确定迭代不长以重新确定削峰功率和填谷功率。
在本发明的实施方式中,当统计的储能总充电量和储能总放电量与储能电站额定容量修订值ζEN相近时输出P1和P2,其中EN为储能电站额定容量。但实际工程中储能电站SOC需限定在一定范围内,不允许满充满放运行,因此引入变量ζ,ζ通常取值0.6~0.8,计算完成后,可根据含储能的快充电站区域负荷及储能电站实时工况适当调整P1和P2的值。
在步骤103,根据所述含储能的快充电站区域的参数数据确定考虑电池寿命的含储能的快充电站区域储能电站的运行控制目标函数。
优选地,其中所述运行控制目标函数为:
其中,F1为储能电站提高含储能的快充电站区域经济效益的评价函数;F2为提高储能电站寿命的评价函数;F3为SOC波动评价目标函数;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;p(t)为t时刻含储能的快充电站区域的商用电实时电价;Cprofit_max为储能电站SOC经历一个预设倍数的循环周期,且在谷时段充电、峰时段放电时的盈利值的最大值,预设倍数的取值范围为0.2-0.8;Nmax[(DOD(i)]表示第i次放电深度对应的电池循环寿命(次数);SOCmax(T)为该时段内SOC最大值,SOCav(T)为T+Δt时段内SOC平均值。
优选地,其中所述方法还包括:采用雨流计数法确定电池的放电深度。
在本发明的实施方式中,采用雨流计数法确定电池的放电深度,计算规则如下:
(1)将SOC曲线顺时针旋转90°,并标注出所有极值点,极值点如图5所示;其中,SOC曲线如图4所示。
(2)雨流计数的起点依次从SOC变化曲线每一个极值的内边开始,即从1,2,3···等尖点开始。
(3)雨流在流到峰值处(即屋檐处)竖直下滴,直到对面处有一个比开始时最大值(最小值)更正的最大值(或更负的最小值)为止。
(4)当雨流遇到来自上面屋顶留下的雨时,就停止流动,并计1个循环。
(5)根据雨每个流的起始点记录循环,并单独取出,取出结果如图6所示。
(6)确定每个雨流对应的纵轴长度即为为该循环的放电深度。
在运用雨流计数法时,可根据电池SOC曲线得到电池工作DOD的n次循环,记为DOD(1),DOD(2),···,DOD(n),则电池寿命衰减率可表示为:
若电池经历了N次循环,则电池剩余寿命可表示为:
例如,通过雨流计数法计算电池DOD时,计算某磷酸铁锂电池循环寿命与DOD的对应关系如下表1所示。
表1磷酸铁锂循环寿命与电池DOD关系表
根据表1中的数据,采用5阶函数可以拟合电池DOD与其循环寿命(次数)的关系为:
N=-1302DOD5+4427DOD3
-8925DOD+10500。
在本发明的实施方式中,运用主观赋值法确定各子目标函数权重,目标函数F1权重λ1值越大,储能电站每日带来的经济效益越明显,但相应的电池寿命衰减越严重;目标函数F2、F3的权重λ2、λ3的值越大,越有利于减缓电池寿命衰减,但储能电站带来的日收益也会相应减少。
在步骤104,通过遗传算法确定储能电站的充放电目标功率。
优选地,其中在通过遗传算法确定储能电站的充放电目标功率时的约束条件,包括:
电网功率平衡约束为:
Pg=Pload-Ppv-PESS,
储能电站SOC及功率约束为:
SOCmin≤SOC(t)≤SOCmax,
0≤Pc≤Pc_e,
0≤Pd≤Pd_e,
并且,在储能电站充电过程中,其SOC还需满足:
在储能电站放电过程中,其SOC还需满足:
其中,Pg为电网向含储能的快充电站区域注入的功率;Ppv为含储能的快充电站区域的分布式光伏发电功率;Pload为含储能的快充电站区域的各类负荷的总功率;PESS为储能系统吸收或发出的功率;SOC(t)为t时刻储能电站电池的SOC值;SOCmin为储能电站电池SOC的最小约束值;SOCmax为储能电站电池SOC的最大约束值;Pc_e为储能电站的额定充电功率;Pd_e为储能电站的额定放电功率;γ为储能电站自放电率,SOC(t)为t时刻的储能电站的SOC状态值;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;Ee为储能电站额定容量。
在本发明的实施方式中,SOCmin取0.2,SOCmax取0.8。
在步骤105,通过预设的控制策略确定储能电站的最优充放电功率。
优选地,其中所述通过预设的控制策略确定储能电站的最优充放电功率,包括:
当负荷高于削峰功率时,储能电站放电,放电功率为:
Pd=min(Pload-P1,Pd_obj),
当负荷低于填谷功率时,储能系统充电,充电功率为:
Pc=min(P2-Pload,Pc_obj),
其中,P1为削峰功率;Pd为放电功率;Pload为实际负荷功率;Pd_obj为通过遗传算法确定的储能电站的放电目标功率;P2为填谷功率;Pc为充电功率或放电功率;Pload为实际负荷功率;Pc_obj为基于考虑电池寿命的储能电站运行控制方法计算得出的充电功率值。
在本发明的实施方式中,储能电站的最优充放电功率Pc表示如下,包括:
(1)负荷高于削峰功率P1时,储能系统放电,放电功率可表示为:
Pd=min(Pload-P1,Pd_obj),
(2)当负荷低于填谷功率P2时,储能系统充电,充电功率可表示为:
Pc=min(P2-Pload,Pc_obj),
其中,P1为削峰功率;Pd为放电功率;Pload为实际负荷功率;Pd_obj为通过遗传算法确定的储能电站的放电目标功率;P2为填谷功率;Pc为充电功率或放电功率;Pload为实际负荷功率;Pc_obj为基于考虑电池寿命的储能电站运行控制方法计算得出的充电功率值。
(3)当负荷处于削峰功率P1和填谷功率P2区间时,储能无需完成削峰填谷任务,但考虑到削峰功率P1和填谷功率P2根据预测负荷得出,如果削峰功率P1过小,后期会出现储能放电功率不足的问题,显然这会增重电力系统负担,造成经济损失。因此,当含储能的快充电站区域负荷处于P1和P2之间时,如果后期储能功率不足,可适当对储能系统进行充电,但要综合考虑储能系统SOC状态以及电池寿命等因素的影响,不能出现储能系统过度充电的情况,此时储能充电功率由本发明所建立的数学模型求得,无法严格解析表达,可简单表示为:
PC=f(SOC,DOD,R,L),
其中,其中SOC为储能电站荷电状态状态、DOD为储能电站放电深度R为储能电站剩余电池寿命,除此之外还有时序功率、负荷及光伏功率等其他因素的影响。
图7为根据本发明实施方式的利用遗传算法确定储能电站充放电功率最优个体的流程图,首先根据典型日负荷预测曲线确定削峰功率P1和填谷功率P2,然后产生储能电站充放电功能率初始群体,通过执行遗传算法的选择、交叉、变异操作,输出通过遗传算法得出的储能电站充放电功率计算结果,然后再根据步骤105所述的方法确定储能电站的最终充放电功率。
图8为根据本发明实施方式的考虑电池寿命的含储能的快充电站区域储能电站运行控制系统800的结构示意图。如图8所示,本发明的实施方式提供的考虑电池寿命的含储能的快充电站区域储能电站运行控制系统800,包括:参数数据获取单元801、削峰功率和填谷功率确定单元802、运行控制目标函数确定单元803、充放电目标功率确定单元804和最优充放电功率确定单元805。
优选地,所述参数数据获取单元801,用于获取含储能的快充电站区域的参数数据;所述参数数据包括:负荷、光伏、日前负荷预测功率、日前光伏预测功率、光伏装机容量和储能电站装机容量。
优选地,所述削峰功率和填谷功率确定单元802,用于根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率。
优选地,其中所述削峰功率和填谷功率确定单元802,包括:净负荷功率确定模块8021、迭代步长确定模块8022、计算模块8023、迭代模块8024和削峰功率和填谷功率确定模块8025。
优选地,所述净负荷功率确定模块8021,用于根据所述日前负荷预测功率和日前光伏预测功率确定日前预测最大净负荷曲线的净负荷最大功率和净负荷最小功率。
优选地,所述迭代步长确定模块8022,用于确定迭代步长。
优选地,所述计算模块8023,用于计算所述净负荷最大功率和迭代步长的差值为当前的削峰功率,计算所述净负荷最小功率和迭代步长的和作为当前的填谷功率。
优选地,所述迭代模块8023,用于根据迭代步长反复调整削峰功率和填谷功率的值,直至统计的储能电站的总充电量和总放电量均大于储能电站的额定容量的修订值的差值。
优选地,所述削峰功率和填谷功率确定模块8025,若当前的储能电站的总充电量与总放电量的差值的绝对值在预设范围内,则确定当前的削峰功率和填谷功率即为出储能电站的削峰功率和填谷功率;反之,进入迭代步长确定模块8022,重新确定迭代步长。
优选地,所述运行控制目标函数确定单元803,用于根据所述含储能的快充电站区域的参数数据确定考虑电池寿命的含储能的快充电站区域储能电站的运行控制目标函数。
优选地,其中所述运行控制目标函数为:
其中,F1为储能电站提高含储能的快充电站区域经济效益的评价函数;F2为提高储能电站寿命的评价函数;F3为SOC波动评价目标函数;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;p(t)为t时刻含储能的快充电站区域的商用电实时电价;Cprofit_max为储能电站SOC经历一个预设倍数的循环周期,且在谷时段充电、峰时段放电时的盈利值的最大值,预设倍数的取值范围为0.2-0.8;Nmax[(DOD(i)]表示第i次放电深度对应的电池循环寿命(次数);SOCmax(T)为该时段内SOC最大值,SOCav(T)为T+Δt时段内SOC平均值。
优选地,其中所述运行控制目标函数确定单元,还包括:放电深度计算模块,用于采用雨流计数法确定电池的放电深度。
优选地,所述充放电目标功率确定单元804,用于通过遗传算法确定储能电站的充放电目标功率。
优选地,其中在所述运行控制目标函数确定单元,通过遗传算法确定储能电站的充放电目标功率时的约束条件,包括:
电网功率平衡约束为:
Pg=Pload-Ppv-PESS,
储能电站SOC及功率约束为:
SOCmin≤SOC(t)≤SOCmax,
0≤Pc≤Pc_e,
0≤Pd≤Pd_e,
并且,在储能电站充电过程中,其SOC还需满足:
在储能电站放电过程中,其SOC还需满足:
其中,Pg为电网向含储能的快充电站区域注入的功率;Ppv为含储能的快充电站区域的分布式光伏发电功率;Pload为含储能的快充电站区域的各类负荷的总功率;PESS为储能系统吸收或发出的功率;SOC(t)为t时刻储能电站电池的SOC值;SOCmin为储能电站电池SOC的最小约束值;SOCmax为储能电站电池SOC的最大约束值;Pc_e为储能电站的额定充电功率;Pd_e为储能电站的额定放电功率;γ为储能电站自放电率,SOC(t)为t时刻的储能电站的SOC状态值;ηd和ηc(t)分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;Ee为储能电站额定容量。
优选地,所述最优充放电功率确定单元805,用于通过预设的控制策略确定储能电站的最优充放电功率。
优选地,其中所述最优放电功率确定单元,通过预设的控制策略确定储能电站的最优充放电功率,包括:
当负荷高于削峰功率时,储能电站放电,放电功率为:
Pd=min(Pload-P1,Pd_obj),
当负荷低于填谷功率时,储能系统充电,充电功率为:
Pc=min(P2-Pload,Pc_obj),
其中,P1为削峰功率;Pd为放电功率;Pload为实际负荷功率;Pd_obj为通过遗传算法确定的储能电站的放电目标功率;P2为填谷功率;Pc为充电功率或放电功率;Pload为实际负荷功率;Pc_obj为基于考虑电池寿命的储能电站运行控制系统计算得出的充电功率值。
本发明的实施例的考虑电池寿命的含储能的快充电站区域储能电站运行控制系统800与本发明的另一个实施例的考虑电池寿命的含储能的快充电站区域储能电站运行控制方法100相对应,在此不再赘述。
已经通过参考少量实施方式描述了本发明。然而,本领域技术人员所公知的,正如附带的专利权利要求所限定的,除了本发明以上公开的其他的实施例等同地落在本发明的范围内。
通常地,在权利要求中使用的所有术语都根据他们在技术领域的通常含义被解释,除非在其中被另外明确地定义。所有的参考“一个/所述/该[装置、组件等]”都被开放地解释为所述装置、组件等中的至少一个实例,除非另外明确地说明。这里公开的任何方法的步骤都没必要以公开的准确的顺序运行,除非明确地说明。
Claims (12)
1.一种考虑电池寿命的储能电站运行控制方法,其特征在于,所述方法包括:
获取含储能的快充电站区域的参数数据;所述参数数据包括:负荷、光伏、日前负荷预测功率、日前光伏预测功率、光伏装机容量和储能电站装机容量;
根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率;
根据所述含储能的快充电站区域的参数数据确定考虑电池寿命的含储能的快充电站区域储能电站的运行控制目标函数;
通过遗传算法确定储能电站的充放电目标功率;
通过预设的控制策略确定储能电站的最优充放电功率。
2.根据权利要求1所述的方法,其特征在于,所述根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率,包括:
步骤1,根据所述日前负荷预测功率和日前光伏预测功率确定日前预测最大净负荷曲线的净负荷最大功率和净负荷最小功率;
步骤2,确定迭代步长;
步骤3,计算所述净负荷最大功率和迭代步长的差值为当前的削峰功率,计算所述净负荷最小功率和迭代步长的和作为当前的填谷功率;
步骤4,根据迭代步长反复调整削峰功率和填谷功率的值,直至统计的储能电站的总充电量和总放电量均大于储能电站的额定容量的修订值的差值;
步骤5,若当前的储能电站的总充电量与总放电量的差值的绝对值在预设范围内,则确定当前的削峰功率和填谷功率即为储能电站的削峰功率和填谷功率;反之,返回步骤2,重新确定迭代步长。
3.根据权利要求1所述的方法,其特征在于,所述运行控制目标函数为:
其中,F1为储能电站提高含储能的快充电站区域经济效益的评价函数;F2为提高储能电站寿命的评价函数;F3为SOC波动评价目标函数;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能系电站放电功率和充电功率;p(t)为t时刻含储能的快充电站区域的商用电实时电价;Cprofit_max为储能电站SOC经历一个预设倍数的循环周期,且在谷时段充电、峰时段放电时的盈利值的最大值;Nmax[(DOD(i)]表示第i次放电深度对应的电池循环寿命(次数);SOCmax(T)为该时段内SOC最大值,SOCav(T)为T+Δt时段内SOC平均值。
4.根据权利要求3所述的方法,其特征在于,所述方法还包括:
采用雨流计数法确定电池的放电深度。
5.根据权利要求1所述的方法,其特征在于,在通过遗传算法确定储能电站的充放电目标功率时的约束条件,包括:
电网功率平衡约束为:
Pg=Pload-Ppv-PESS,
储能电站SOC及功率约束为:
SOCmin≤SOC(t)≤SOCmax,
0≤Pc≤Pc_e,
0≤Pd≤Pd_e,
并且,在储能电站充电过程中,其SOC还需满足:
在储能电站放电过程中,其SOC还需满足:
其中,Pg为电网向含储能的快充电站区域注入的功率;Ppv为含储能的快充电站区域的分布式光伏发电功率;Pload为含储能的快充电站区域的各类负荷的总功率;PESS为储能系统吸收或发出的功率;SOC(t)为t时刻储能电站电池的SOC值;SOCmin为储能电站电池SOC的最小约束值;SOCmax为储能电站电池SOC的最大约束值;Pc_e为储能电站的额定充电功率;Pd_e为储能电站的额定放电功率;γ为储能电站自放电率,SOC(t)为t时刻的储能电站的SOC状态值;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;Ee为储能电站额定容量。
6.根据权利要求1所述的方法,其特征在于,所述通过预设的控制策略确定储能电站的最优充放电功率,包括:
当负荷高于削峰功率时,储能电站放电,放电功率为:
Pd=min(Pload-P1,Pd_obj),
当负荷低于填谷功率时,储能系统充电,充电功率为:
Pc=min(P2-Pload,Pc_obj),
其中,P1为削峰功率;Pd为放电功率;Pload为实际负荷功率;Pd_obj为通过遗传算法确定的储能电站的放电目标功率;P2为填谷功率;Pc为充电功率或放电功率;Pload为实际负荷功率;Pc_obj为基于通过遗传算法确定储能电站的充电目标功率值。
7.一种考虑电池寿命的储能电站运行控制系统,其特征在于,所述系统包括:
参数数据获取单元,用于获取含储能的快充电站区域的参数数据;所述参数数据包括:负荷、光伏、日前负荷预测功率、日前光伏预测功率、光伏装机容量和储能电站装机容量;
削峰功率和填谷功率确定单元,用于根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率;
运行控制目标函数确定单元,用于根据所述含储能的快充电站区域的参数数据确定考虑电池寿命的含储能的快充电站区域储能电站的运行控制目标函数;
充放电目标功率确定单元,用于通过遗传算法确定储能电站的充放电目标功率;
最优充放电功率确定单元,用于通过预设的控制策略确定储能电站的最优充放电功率。
8.根据权利要求7所述的系统,其特征在于,所述削峰功率和填谷功率确定单元,根据所述含储能的快充电站区域的参数数据确定储能电站的削峰功率和填谷功率,包括:
净负荷功率确定模块,用于根据所述日前负荷预测功率和日前光伏预测功率确定日前预测最大净负荷曲线的净负荷最大功率和净负荷最小功率;
迭代步长确定模块,用于确定迭代步长;
计算模块,用于计算所述净负荷最大功率和迭代步长的差值为当前的削峰功率,计算所述净负荷最小功率和迭代步长的和作为当前的填谷功率;
迭代模块,用于根据迭代步长反复调整削峰功率和填谷功率的值,直至统计的储能电站的总充电量和总放电量均大于储能电站的额定容量的修订值的差值;
削峰功率和填谷功率确定模块,若当前的储能电站的总充电量与总放电量的差值的绝对值在预设范围内,则确定当前的削峰功率和填谷功率即为出储能电站的削峰功率和填谷功率;反之,进入迭代步长确定模块,重新确定迭代步长。
9.根据权利要求7所述的系统,其特征在于,所述运行控制目标函数为:
其中,F1为储能电站提高含储能的快充电站区域经济效益的评价函数;F2为提高储能电站寿命的评价函数;F3为SOC波动评价目标函数;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;p(t)为t时刻含储能的快充电站区域的商用电实时电价;Cprofit_max为储能电站SOC经历一个预设倍数的循环周期,且在谷时段充电、峰时段放电时的盈利值的最大值;Nmax[(DOD(i)]表示第i次放电深度对应的电池循环寿命(次数);SOCmax(T)为该时段内SOC最大值,SOCav(T)为T+Δt时段内SOC平均值。
10.根据权利要求9所述的系统,其特征在于,所述系统还包括:
放电深度计算模块,用于采用雨流计数法确定电池的放电深度。
11.根据权利要求7所述的系统,其特征在于,在所述运行控制目标函数确定单元,通过遗传算法确定储能电站的充放电目标功率时的约束条件,包括:
电网功率平衡约束为:
Pg=Pload-Ppv-PESS,
储能电站SOC及功率约束为:
SOCmin≤SOC(t)≤SOCmax,
0≤Pc≤Pc_e,
0≤Pd≤Pd_e,
并且,在储能电站充电过程中,其SOC还需满足:
在储能电站放电过程中,其SOC还需满足:
其中,Pg为电网向含储能的快充电站区域注入的功率;Ppv为含储能的快充电站区域的分布式光伏发电功率;Pload为含储能的快充电站区域的各类负荷的总功率;PESS为储能系统吸收或发出的功率;SOC(t)为t时刻储能电站电池的SOC值;SOCmin为储能电站电池SOC的最小约束值;SOCmax为储能电站电池SOC的最大约束值;Pc_e为储能电站的额定充电功率;Pd_e为储能电站的额定放电功率;γ为储能电站自放电率,SOC(t)为t时刻的储能电站的SOC状态值;ηd和ηc分别为电池的放电效率和充电效率;Pd(t)和Pc(t)分别为t时刻储能电站的放电功率和充电功率;Ee为储能电站额定容量。
12.根据权利要求7所述的系统,其特征在于,所述最优放电功率确定单元,通过预设的控制策略确定储能电站的最优充放电功率,包括:
当负荷高于削峰功率时,储能电站放电,放电功率为:
Pd=min(Pload-P1,Pd_obj),
当负荷低于填谷功率时,储能系统充电,充电功率为:
Pc=min(P2-Pload,Pc_obj),
其中,P1为削峰功率;Pd为放电功率;Pload为实际负荷功率;为通过遗传算法确定的储能电站的放电目标功率;P2为填谷功率;Pc为充电功率或放电功率;Pload为实际负荷功率;Pc_obj为基于通过遗传算法确定储能电站的充电目标功率值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010047125.5A CN111244991B (zh) | 2020-01-16 | 2020-01-16 | 一种考虑电池寿命的储能电站运行控制方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010047125.5A CN111244991B (zh) | 2020-01-16 | 2020-01-16 | 一种考虑电池寿命的储能电站运行控制方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111244991A true CN111244991A (zh) | 2020-06-05 |
CN111244991B CN111244991B (zh) | 2023-01-17 |
Family
ID=70874617
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010047125.5A Active CN111244991B (zh) | 2020-01-16 | 2020-01-16 | 一种考虑电池寿命的储能电站运行控制方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111244991B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112260300A (zh) * | 2020-10-10 | 2021-01-22 | 国网能源研究院有限公司 | 储能配置及最佳延缓年限的确定方法及装置 |
CN112803464A (zh) * | 2021-03-16 | 2021-05-14 | 中国电力科学研究院有限公司 | 一种储能系统充放电控制方法、系统、设备和存储介质 |
CN112858942A (zh) * | 2020-11-30 | 2021-05-28 | 浙江大有实业有限公司杭州科技发展分公司 | 一种退役电池健康状态均衡方法 |
CN113131503A (zh) * | 2021-04-25 | 2021-07-16 | 山东电工电气集团有限公司 | 基于多电池组soc一致性的储能电站能量管理方法 |
CN114006442A (zh) * | 2021-11-26 | 2022-02-01 | 国网甘肃省电力公司电力科学研究院 | 一种考虑荷电状态一致的电池储能电站能量管理方法 |
CN115579923A (zh) * | 2022-10-18 | 2023-01-06 | 国网上海能源互联网研究院有限公司 | 一种用于多场景多目标的储能优化运行方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102361327A (zh) * | 2011-10-17 | 2012-02-22 | 张家港智电可再生能源与储能技术研究所有限公司 | 一种考虑电池寿命的电池储能系统削峰填谷方法 |
CN102664401A (zh) * | 2012-04-16 | 2012-09-12 | 中国电力科学研究院 | 一种基于电池寿命模型的微电网控制方法 |
US20160197514A1 (en) * | 2013-07-30 | 2016-07-07 | Zhangjiakou Wind And Solar Power Energy Demonstration Station Co. Ltd., State Grid Xin | Energy management method of multi-type battery energy storage power station considering charge and discharge rates |
CN109245143A (zh) * | 2018-11-26 | 2019-01-18 | 江苏云上电力科技有限公司 | 一种考虑锂离子电池寿命的储能调峰电站优化运行方法 |
-
2020
- 2020-01-16 CN CN202010047125.5A patent/CN111244991B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102361327A (zh) * | 2011-10-17 | 2012-02-22 | 张家港智电可再生能源与储能技术研究所有限公司 | 一种考虑电池寿命的电池储能系统削峰填谷方法 |
CN102664401A (zh) * | 2012-04-16 | 2012-09-12 | 中国电力科学研究院 | 一种基于电池寿命模型的微电网控制方法 |
US20160197514A1 (en) * | 2013-07-30 | 2016-07-07 | Zhangjiakou Wind And Solar Power Energy Demonstration Station Co. Ltd., State Grid Xin | Energy management method of multi-type battery energy storage power station considering charge and discharge rates |
CN109245143A (zh) * | 2018-11-26 | 2019-01-18 | 江苏云上电力科技有限公司 | 一种考虑锂离子电池寿命的储能调峰电站优化运行方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112260300A (zh) * | 2020-10-10 | 2021-01-22 | 国网能源研究院有限公司 | 储能配置及最佳延缓年限的确定方法及装置 |
CN112858942A (zh) * | 2020-11-30 | 2021-05-28 | 浙江大有实业有限公司杭州科技发展分公司 | 一种退役电池健康状态均衡方法 |
CN112803464A (zh) * | 2021-03-16 | 2021-05-14 | 中国电力科学研究院有限公司 | 一种储能系统充放电控制方法、系统、设备和存储介质 |
CN112803464B (zh) * | 2021-03-16 | 2022-05-31 | 中国电力科学研究院有限公司 | 一种储能系统充放电控制方法、系统、设备和存储介质 |
CN113131503A (zh) * | 2021-04-25 | 2021-07-16 | 山东电工电气集团有限公司 | 基于多电池组soc一致性的储能电站能量管理方法 |
CN114006442A (zh) * | 2021-11-26 | 2022-02-01 | 国网甘肃省电力公司电力科学研究院 | 一种考虑荷电状态一致的电池储能电站能量管理方法 |
CN115579923A (zh) * | 2022-10-18 | 2023-01-06 | 国网上海能源互联网研究院有限公司 | 一种用于多场景多目标的储能优化运行方法及装置 |
CN115579923B (zh) * | 2022-10-18 | 2024-02-23 | 国网上海能源互联网研究院有限公司 | 一种用于多场景多目标的储能优化运行方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111244991B (zh) | 2023-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111244991B (zh) | 一种考虑电池寿命的储能电站运行控制方法及系统 | |
Liu et al. | Optimal sizing of a wind-energy storage system considering battery life | |
CN108667052B (zh) | 一种面向虚拟电厂优化运行的多类型储能系统规划配置方法及系统 | |
CN109599856B (zh) | 一种微网多楼宇中电动汽车充放电管理优化方法及装置 | |
Li et al. | On the determination of battery energy storage capacity and short-term power dispatch of a wind farm | |
CN106972516A (zh) | 一种适用于微网的多类型储能多级控制方法 | |
CN107248751A (zh) | 一种实现配电网负荷功率削峰填谷的储能站调度控制方法 | |
CN103485977B (zh) | 风力发电系统实时功率预测的修正方法 | |
CN102005771A (zh) | 风、光、储微电网系统的储能容量选取方法 | |
CN103311942A (zh) | 一种用于配电网削峰填谷的电池储能系统的控制方法 | |
CN109378846B (zh) | 储能电站中储能变流器的电池模块充放电控制方法与装置 | |
CN109102185B (zh) | 一种评估光储充一体化电站全寿命周期经济性的方法和系统 | |
CN110086187A (zh) | 计及负荷特性的储能调峰日前优化调度方法 | |
CN113627762B (zh) | 一种基于激励电价的虚拟电厂调峰方法 | |
CN111509745B (zh) | 一种解决变压器超容问题的储能系统ems放电控制方法 | |
CN111626527A (zh) | 计及可调度电动汽车快/慢充放电形式的智能电网深度学习调度方法 | |
CN113659623B (zh) | 一种基于布林线理论的风储联合系统的优化方法及系统 | |
CN114154790A (zh) | 基于需量管理和柔性负荷的工业园区光储容量配置方法 | |
CN114899856A (zh) | 一种电动汽车充电桩功率调节方法、系统、设备及介质 | |
CN107565587B (zh) | 考虑爬坡功率有限度平抑的光伏电网储能配置方法及系统 | |
CN114944661A (zh) | 一种基于储能系统滚动优化的微电网三阶段优化控制方法 | |
Boyouk et al. | Peak shaving of a grid connected-photovoltaic battery system at helmholtz institute ulm (hiu) | |
Martinsen | A business model for an EV charging station with battery energy storage | |
CN117856261A (zh) | 基于时序模拟仿真计算电网调节能力缺口的方法和装置 | |
CN117575185A (zh) | 一种面向生产计划的储能分时梯级控制经济调度策略方法 |
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 |