CN116599086A - 一种考虑新能源发电与负荷不确定性的储能调频控制方法 - Google Patents
一种考虑新能源发电与负荷不确定性的储能调频控制方法 Download PDFInfo
- Publication number
- CN116599086A CN116599086A CN202310675994.6A CN202310675994A CN116599086A CN 116599086 A CN116599086 A CN 116599086A CN 202310675994 A CN202310675994 A CN 202310675994A CN 116599086 A CN116599086 A CN 116599086A
- Authority
- CN
- China
- Prior art keywords
- energy storage
- power
- representing
- frequency modulation
- battery
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004146 energy storage Methods 0.000 title claims abstract description 161
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000010248 power generation Methods 0.000 title claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 38
- 238000009826 distribution Methods 0.000 claims abstract description 27
- 230000009194 climbing Effects 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000007599 discharging Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 238000013528 artificial neural network Methods 0.000 claims description 6
- 230000015556 catabolic process Effects 0.000 claims description 6
- 238000006731 degradation reaction Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 229910018503 SF6 Inorganic materials 0.000 claims description 3
- 108010074506 Transfer Factor Proteins 0.000 claims description 3
- 230000003828 downregulation Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- SFZCNBIFKDRMGX-UHFFFAOYSA-N sulfur hexafluoride Chemical compound FS(F)(F)(F)(F)F SFZCNBIFKDRMGX-UHFFFAOYSA-N 0.000 claims description 3
- 229960000909 sulfur hexafluoride Drugs 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 claims 2
- 238000005315 distribution function Methods 0.000 abstract description 2
- 230000008901 benefit Effects 0.000 description 5
- 229910052799 carbon Inorganic materials 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 241000588653 Neisseria Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 241000039077 Copula Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- 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/002—Flicker reduction, e.g. compensation of flicker introduced by non-linear load
-
- 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/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Nonlinear Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种考虑新能源发电与负荷不确定性的储能调频控制方法。首先,建立考虑AGC的综合成本最小化的经济调度模型;其次得到AGC信号预测数据的特征并建立其概率分布函数;然后构建AGC信号不确定量的模糊集并用分布式鲁棒优化模型进行求解,得到调频资源的基本调度点;随后根据AGC功率需求的高频分量和储能运行约束,考虑电池储能SOC均衡及调频成本最优化建立电池储能AGC控制模型;最后用分布式一致性算法对AGC控制模型进行求解,得到储能集群的实时功率输出,维持电网频率稳定。本发明不仅通过综合考虑经济调度与频率控制,实现始终以经济成本最小化来调度和控制配电网内的调频资源,同时考虑了储能的SOC均衡,实现调频控制的可持续性。
Description
技术领域
本发明涉及一种考虑新能源发电与负荷不确定性的储能调频控制方法,属于电力系统技术领域。
背景技术
近年来,为应对能源供应、环境保护以及气候变化带来的巨大挑战,能源革命在世界范围内逐渐兴起,世界能源结构发生了重大变化。“双碳”战略的提出就要求我国能源结构需要加速向多元化、清洁化转变,同时发展重点由传统能源向新能源转变。以光伏、风力为代表的新能源发电,不仅优化了能源结构,同时还具有清洁、低碳、可再生等诸多优点。因此,构建以新能源为主体的新型电力系统对实现“双碳”目标、推动我国高质量发展意义重大。
然而,随着大规模新能源的并网,骨干电网消纳新能源的压力开始逐渐显现,弃风弃光现象频发。为了维持骨干电网的稳定,配电网作为能源生产、转换、消费的关键环节,更是可再生能源消纳的支撑平台,在双碳战略的牵引下,配电网逐渐成为电力系统核心。传统配电网通过经济调度和频率控制维持有功功率的平衡,但随着新能源在配网中渗透率的增加,频率调节由原本仅应对负荷不确定性,转而同时需要新能源发电和负荷的双重不确定性,不仅对调频资源的容量和爬坡率的要求变得更高,同时频率调节辅助服务的成本也增加,尤其是类似于储能等快速调频资源的加入,使得对调频资源的控制变得更加复杂。
当前的多数研究从调频资源运行的可持续性入手来提高调度的经济性能,忽略了一个合理且考虑全面的优化目标函数对于频率控制经济性的影响。
申请号为201910862624.7、名称为《一种考虑风光储提供柔性爬坡容量的实时经济调度方法》的发明专利申请(下文简称“专利一”)提出一种考虑风光储提供柔性爬坡容量,大大减轻了常规机组提供FRC的压力,同时通过提供多种提供FRC的“灵活源”组合,执行实时经济调度,灵活性高。
申请号为201910294381.1、名称为《一种基于储能装置SOC动态约束的微网经济调度方法》的发明专利申请(下文简称“专利二”)提出一种考虑储能SOC动态约束的经济调度方法,解决了现有为保障重要负荷的储能装置预留容量不足问题,确保在任意时刻脱网都能保障其供电不中断,减少了损失,提高了经济效益。
上述专利一以及专利二从经济调度的约束条件入手,确保运行时的持续性来提高经济调度的经济效益,但忽视了新能源高渗透率下频繁的频率调节导致AGC里程成本的显著增加,没有在目标函数中考虑调频成本,进而导致经济效益的下降。
申请号为202010333512.5、名称为《一种配电网侧规模化储能参与调频的双层控制方法》的发明专利申请(下文简称“专利三”)提出一种储能调频的双层控制方法,其中上层控制考虑经济性,下层控制考虑调频性能,在维持系统频率稳定的同时降低成本。
专利三考虑了经济性能但忽视了储能的SOC均衡,没有考虑到储能运行的可持续性。
综上所述,这些方法没有综合考虑经济调度与频率控制的协同,在提高经济性的同时维持储能的SOC均衡。
发明内容
本发明要解决的技术问题是:现有的。储能调频控制方法没有综合考虑经济调度与频率控制的协同,在提高经济性的同时维持储能的SOC均衡
为了解决上述技术问题,本发明的技术方案是提供了一种考虑新能源与负荷双重不确定性的储能调频控制方法,其特征在于,包括以下步骤:
步骤S1、获取经济调度模型的约束条件,建立考虑自动发电控制成本的综合成本最小的经济调度优化模型;
步骤S2、根据历史AGC信号数据,建立AGC信号特征值的概率分布并构建模糊集,将随机优化模型转化为分布式鲁棒优化模型并求解,得到调频资源的基本功率调度点,其中,AGC表示自动发电控制成本;
步骤S3、根据AGC功率需求的高频分量和储能运行约束,考虑电池储能SOC均衡及调频成本最优化建立电池储能AGC控制模型,如下式所示:
式中:I为储能有效运行范围集合;PSref(t)表示t时刻电池储能需要满足的总功率需求;DS,i(t)表示t时刻第i台电池储能的调频耗量;Pi S(t)表示第i台电池储能的实时充放电功率;Pi SL(t)和Pi SU(t)表示储能满足其运行约束的充放电功率上下限;
步骤S4、对电池储能AGC控制模型用一致性分布式算法求解,不断迭代一致性变量直至得到最优解作为储能的功率输出值,实现电网频率控制。
优选地,步骤S1中,所述考虑自动发电控制成本的综合成本最小的经济调度优化模型的目标函数,如下式所示:
式中:表示电池储能的退化成本;/>和/>分别表示第i台电池储能在第n时刻的充电功率和放电功率;ΕP表示求数学期望;/>表示第j台传统电机的发电成本;/>为第j台电机的发电功率,/>表示第j台传统电机在n时刻的平均输出功率值;/>表示调频资源的里程成本;Mn表示调节里程;PFk,n表示第n时刻第k台设备的参与系数;S表示电池储能;G表示传统电机;{G,S}表示所有调频资源,即传统电机和电池储能。
优选地,步骤S1中,所述约束条件包括实时经济调度的约束条件、考虑AGC功率需求的机会约束以及储能的SOC约束。
优选地,其特征在于,所述实时经济调度的约束条件包括:
功率平衡约束:
功率上下限约束:
功率爬坡约束:
线路容量约束:
式中,是第i台电池储能在第n时刻的输出功率;/>是节点b在n时刻的负荷功率,b表示一个节点;Pk,min和Pk,max分别是传统电机或储能输出功率的上下限;/>表示第k个传统机组或者电池储能在n时刻的输出功率;rrk是调频资源的爬坡率;T表示一个经济调度周期;fl是线路l的容量;SFl,b是线路l上节点b的转移因子;/>分别表示节点b中传统机组、电池储能以及新能源发电的向量值;/>分别表示传统电机、电池储能以及新能源发电的输出功率;
考虑AGC功率需求的机会约束有:
响应AGC信号需求下,调频资源的功率输出在其自身约束范围之内:
式中:PFk,n是标号为k的调频资源的参与系数;P(·)表示该事件发生的概率;是一个经济调度周期中AGC信号特征值中的上分位数,/>是一个经济调度周期中AGC信号特征值中的下分位数;β表示规定的功率输出机会约束满足概率值;
同时,调频资源的爬坡率高于AGC功率需求的调节率:
式中,表示调节率上分位数,/>表示调节率下分位数;P(·)表示该事件发生的概率,β表示规定的爬坡率机会约束满足概率值
所有调频资源的参与系数应介于0-1之间,且参与系数之和应为1,即
储能的SOC约束有:
式中,SOCi,n-1表示第i台电池储能在第n-1时刻的电池储能SOC值;是一个周期内AGC信号向上调节的累计值;/>是一个周期内AGC信号向下调节的累计值;CEi表示第i台电池储能的容量;/>和/>分别表示储能的充电效率和放电效率;/>和/>分别表示电池储能SOC值的上限和下限。
优选地,所述步骤S2包括:
步骤S201、对AGC信号历史数据进行处理,提取累计调节能量累计调节里程Mn、调节幅度分位数/>以及调节率分位数/>然后用LSTM神经网络根据历史数据进行预测,并得到其概率分布;
步骤S202、采用Wasserstein度量从预测的AGC信号样本集中构建模糊集,从而找到一个最小化最坏情况下的决策;
步骤S203、根据步骤S1确立的目标函数及其约束条件,重新表述优化模型如下:
式中:表示为参与经济调度的所有储能的总退化成本;表示所有传统电机的总发电成本,EQ表示对于满足模糊集的最差分布求期望,Βki表示调频参与系数矩阵,/>表示AGC信号的特征值向量,Dki表示传统电机的发电成本;/>表示所有调频资源的调频里程成本,表示调频里程成本的参与系数矩阵;Mm表示机会约束中的系数矩阵,Nm表示和调频资源输出功率相关的函数矩阵,βm表示规定的第m个机会约束式子满足的概率值,J表示所有机会约束的集合;
步骤S204、将AGC功率需求的机会约束转化为条件风险价值约束,随后将条件风险价值约束表述的约束条件乘以风险规避参数并添加到目标函数中;
步骤S205、用MATLAB中带的求解器进行求解,得到调频资源的基本调度点。
优选地,所述步骤S3包括:。
建立单个储能调频耗量函数如下:
式中:DS,i(t)表示t时刻第i台电池储能的调频耗量,其为实时充放电功率和容量偏差的加权平方和,Pi S(t)表示第i台电池储能的实时充放电功率,表示第i台电池储能的容量偏差,CEi(t)表示t时刻第i个电池储能容量,/>表示第i个电池储能额定容量,aS,i和bS,i分别为系数,Δt表示调频控制间隔时间;
储能运行受到额定功率、额定容量和爬坡速率及功率平衡约束的限制,如下所示:
Pi,min≤Pi S≤Pi,max
Pi S(t)∈ΩP,i(t)∈{Pi S(t)|Pi SL(t)≤Pi S(t)≤Pi SU(t)}
式中:和/>表示储能的充放电爬坡速率范围;ΩP,i(t)表示电池储能输出功率的可行区间;Pi SL(t)和Pi SU(t)表示储能满足其运行约束的充放电功率上下限;
考虑储能实时容量受其充放电功率影响,将储能调频耗量函数变换为以储能充放电功率为单一变量的函数形式,如下所示:
DS,i(t)=αS,i[Pi S(t)]2+βS,iPi S(t)+γS,i
式中:αS,i、βS,i和γS,i分别是调频耗量函数二次项、一次项和常数项系数;
则建立AGC的优化模型如下:
式中:I为储能有效运行范围集合;PSref(t)表示t时刻电池储能需要满足的总功率需求。
优选地,所述步骤S4包括:
步骤S401、引入拉格朗日乘子λ,将原优化问题描述为如下所示的新优化问题:
步骤S402、当储能充放电功率在其上、下界内取值时,集群内各储能调频边际成本趋于一致且等于λ*时,优化问题获得最优解,相应公式如下:
步骤S403、不断在相邻储能集群之间迭代一致性变量,直至达到一致,得到最优边际成本和最佳输出功率如下:
式中,βS,i(t)表示电池储能调频损耗函数一次项系数,αS,i表示电池储能调频损耗函数二次项系数,Ωi(t)表示电池储能输出功率可行区间。
与现有技术相比,本发明具有以下有益效果:
(1)本发明综合考虑经济调度与频率控制,实现始终以经济成本最优化实现储能的功率分配并在实时经济调度模型中考虑自动发电控制(Automatic GenerationControl,下文简称AGC)的成本,提高了实时经济调度的经济性;
(2)本发明采用基于Wasserstein度量的分布式鲁棒优化模型,其相对于鲁棒优化模型在具有更优的经济效益的同时兼顾鲁棒性。
(3)本发明在AGC控制模型中考虑了储能的SOC均衡,在降低储能调频成本的同时实现储能的SOC均衡,实现储能调频的可持续性运行。
附图说明
图1为本发明实施流程示意图;
图2为实时经济调度模型流程图;
图3为频率跟踪曲线图。
具体实施方式
下面结合具体实施例进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
本发明公开了一种考虑新能源与负荷双重不确定性的储能调频控制方法,可适用于分布式的电池储能集群。首先,本发明建立考虑AGC成本的综合成本最小化的经济调度模型;其次,本发明根据AGC信号历史数据的几个特征值,得到AGC信号预测数据的特征并建立其概率分布函数;然后,本发明衡量AGC信号特征量预测值的真实概率分布和参考概率分布之间的距离,构建AGC信号不确定量的模糊集并用分布式鲁棒优化模型进行求解,得到调频资源的基本调度点;随后,本发明根据AGC功率需求的高频分量和储能运行约束,考虑电池储能SOC均衡及调频成本最优化建立电池储能AGC控制模型;最后,本发明用分布式一致性算法对AGC控制模型进行求解,得到储能集群的实时功率输出,维持电网频率稳定。
具体而言,如图1所示,本发明公开的一种考虑新能源与负荷双重不确定性的储能调频控制方法包括以下步骤:
步骤S1、获取经济调度模型的约束条件,建立考虑AGC成本的综合成本最小的经济调度优化模型。
作为一种优选实施方式,步骤S1进一步包括以下内容:
建立经济调度模型的实时经济调度目标函数,如下式所示:
式中:表示电池储能的退化成本;/>和/>分别表示第i台电池储能在第n时刻的充电功率和放电功率;ΕP表示求数学期望;/>表示第j台传统电机的发电成本;/>为第j台电机的发电功率,/>表示第j台传统电机在n时刻的平均输出功率;/>表示调频资源的里程成本;Mn表示调节里程;PFk,n表示第n时刻第k台设备的参与系数;S表示电池储能;G表示传统电机;{G,S}表示所有调频资源,即传统电机和电池储能。
经济调度模型的约束条件包括实时经济调度的约束条件,考虑AGC功率需求的机会约束以及储能的SOC约束,其中:
实时经济调度的约束条件有功率平衡约束、功率上下限约束、功率爬坡约束、线路容量约束,分别表示为
功率平衡约束:
功率上下限约束:
功率爬坡约束:
n线路容量约束:
式中,是第i台电池储能在第n时刻的输出功率;/>是节点b在n时刻的负荷功率,b表示一个节点;Pk,min和Pk,max分别是传统电机或储能输出功率的上下限;/>表示第k个传统机组或者电池储能在n时刻的输出功率;rrk是调频资源的爬坡率;T表示一个经济调度周期;fl是线路l的容量;SFl,b是线路l上节点b的转移因子;/>分别表示节点b中传统机组、电池储能以及新能源发电的向量值;/>分别表示传统电机、电池储能以及新能源发电的输出功率。
考虑AGC功率需求的机会约束有:
响应AGC信号需求下,调频资源的功率输出应在其自身约束范围之内:
式中:PFk,n是标号为k的调频资源的参与系数;P(·)表示该事件发生的概率;是一个经济调度周期中AGC信号特征值中的上分位数,/>是一个经济调度周期中AGC信号特征值中的下分位数;β表示规定的功率机会约束应满足概率值。
同时,调频资源的爬坡率应高于AGC功率需求的调节率:
式中,表示调节率上分位数,/>表示调节率下分位数,P(·)表示该事件发生的概率,β表示规定的爬坡率机会约束应满足概率值。
所有调频资源的参与系数应介于0-1之间,且参与系数之和应为1,即
储能的SOC约束有:
式中,SOCi,n-1表示第i台电池储能在第n-1时刻的电池储能SOC值;是一个周期内AGC信号向上调节的累计值;/>是一个周期内AGC信号向下调节的累计值;CEi表示第i台电池储能的容量;/>和/>分别表示储能的充电效率和放电效率;/>和/>分别表示电池储能SOC值的上限和下限。
步骤S2、根据历史AGC信号数据,建立AGC信号特征值的概率分布并构建模糊集,将随机优化模型转化为分布式鲁棒优化模型并求解,得到调频资源的基本功率调度点。
作为一种优选实施方式,步骤S2进一步包括以下内容:
步骤S201、对AGC信号历史数据进行处理,提取以下几个特征值分别为:累计调节能量累计调节里程Mn、调节幅度分位数/>以及调节率分位数/>然后用LSTM神经网络根据历史数据进行预测,并得到其概率分布。
步骤S202、采用Wasserstein度量从预测的AGC信号样本集中构建模糊集,从而找到一个最小化最坏情况下的决策。基于Wasserstein距离的鲁棒优化算法需要先粗略设定随机变量所在集合,一般取多面集Ξ={ξ|Cξ≤d}。同时,瓦瑟斯坦球提供的半径量化置信上限相对于其他方法来说可以实现更优异的性能,其具体定义如下:
基于瓦瑟斯坦度量,确定AGC信号数据的模糊集如下:
式中:Q表示多面集上所有支持的AGC概率分布;是AGC信号的基于Copula函数预测的经验分布;/>表示两个分布之间的Wasserstein距离;S(Ξ)表示满足半径要求的Wasserstein球;ε是瓦瑟斯坦球半径;/>则是以预测的分布为中心、以ε为半径的分布集,即AGC信号分布的模糊集。
步骤S203、根据步骤S1确立的目标函数及其约束条件,重新表述优化模型如下(除此之外,还需满足先前所述的实时经济调度约束):
式中:表示为参与经济调度的所有储能的总退化成本;表示所有传统电机的总发电成本,EQ表示满足模糊集的最差分布的期望值,Βki表示调频参与系数矩阵,/>表示AGC信号的特征值向量,Dki表示传统电机的发电成本矩阵形式;/>表示所有调频资源的调频里程成本,/>表示调频里程成本的参与系数矩阵;约束条件则是对上述AGC功率需求以及储能SOC机会约束的简写,Mm表示机会约束的系数矩阵,Nm表示调频资源输出功率相关的函数矩阵,βm表示第m个机会约束式子满足的概率值,J表示所有机会约束的集合。
步骤S204、将AGC功率需求的机会约束转化为条件风险价值(CVaR)约束,其约束条件表示如下:
式中,表示第m个机会约束的风险值。随后将CvaR表述的约束条件乘以风险规避参数并添加到目标函数中,其表达如下:
Q表示多面集上所有支持的AGC概率分布;是AGC信号的经验分布;ρ为风险规避参数;/>表示第m个机会约束的风险值。
步骤S205、用MATLAB中带的求解器进行求解,得到调频资源的基本调度点。
步骤S3、根据AGC功率需求的高频分量和储能运行约束,考虑电池储能SOC均衡及调频成本最优化建立电池储能AGC控制模型。
作为一种优选实施方式,步骤S3进一步包括以下内容:
获取AGC信号的高频分量,确立面向调频需求确立储能集群均衡控制的目标函数,包括以下步骤:
步骤S301、建立单个储能调频耗量函数如下:
式中:DS,i(t)表示t时刻第i台电池储能的调频耗量,其为实时充放电功率和容量偏差的加权平方和,Pi S(t)表示第i台电池储能的实时充放电功率,表示第i台电池储能的容量偏差,CEi(t)表示t时刻第i个电池储能容量,/>表示第i个电池储能额定容量,aS,i和bS,i分别为系数,Δt表示调频控制间隔时间。
进一步地,储能运行受到额定功率、额定容量和爬坡速率及功率平衡约束的限制,如下所示:
Pi,min≤Pi S≤Pi,max
Pi S(t)∈ΩP,i(t)∈{Pi S(t)|Pi SL(t)≤Pi S(t)≤Pi SU(t)}
式中:和/>表示储能的充放电爬坡速率范围;ΩP,i(t)表示电池储能输出功率的可行域;Pi SL(t)和Pi SU(t)表示储能满足其运行约束的充放电功率上下限。
进一步地,考虑储能实时容量受其充放电功率影响,可将储能调频耗量函数变换为以储能充放电功率为单一变量的函数形式,如下所示:
DS,i(t)=αS,i[Pi S(t)]2+βS,iPi S(t)+γS,i
式中:αS,i、βS,i和γS,i分别是调频耗量函数二次项、一次项和常数项系数。
则综合上述式子建立AGC的优化模型如下:
式中:I为储能有效运行范围集合;PSref(t)表示t时刻电池储能需要满足的总功率需求。
步骤S4、对电池储能AGC优化模型用一致性分布式算法求解,不断迭代一致性变量直至得到最优解作为储能的功率输出值,实现电网频率控制。
作为一种优选实施方式,步骤S4进一步包括以下内容:
用一致性分布式算法求解包括以下步骤:
步骤S401、引入拉格朗日乘子λ,可将原优化问题描述为如下所示的新优化问题,两者具有相同的最优解。
步骤S402、根据Lagrange极值条件可推算得到,当储能充放电功率在其上、下界内取值时,集群内各储能调频边际成本趋于一致且等于λ*时,优化问题获得最优解,相应公式如下:
步骤S403、不断在相邻储能集群之间迭代一致性变量,直至达到一致,得到最优边际成本和最佳输出功率如下:
式中,βS,i(t)表示电池储能调频损耗函数一次项系数,αS,i表示电池储能调频损耗函数二次项系数,Ωi(t)表示电池储能输出功率可行域。
综上,得到储能的输出功率分配方案。
本发明公开的方法不仅通过综合考虑经济调度与频率控制,实现始终以经济成本最小化来调度和控制配电网内的调频资源,同时考虑了储能的SOC均衡,实现调频控制的可持续性。
Claims (7)
1.一种考虑新能源与负荷双重不确定性的储能调频控制方法,其特征在于,包括以下步骤:
步骤S1、获取经济调度模型的约束条件,建立考虑自动发电控制成本的综合成本最小的经济调度优化模型;
步骤S2、根据历史AGC信号数据,建立AGC信号特征值的概率分布并构建模糊集,将随机优化模型转化为分布式鲁棒优化模型并求解,得到调频资源的基本功率调度点,其中,AGC表示自动发电控制成本;
步骤S3、根据AGC功率需求的高频分量和储能运行约束,考虑电池储能SOC均衡及调频成本最优化建立电池储能AGC控制模型,如下式所示:
式中:I为储能有效运行范围集合;PSref(t)表示电池储能需要满足的总功率需求;DS,i(t)表示t时刻第i台电池储能的调频耗量;Pi S(t)表示第i台电池储能的实时充放电功率;Pi SL(t)和Pi SU(t)表示储能满足其运行约束的充放电功率上下限;
步骤S4、对电池储能AGC控制模型用一致性分布式算法求解,不断迭代一致性变量直至得到最优解作为储能的功率输出值,实现电网频率控制。
2.如权利要求1所述的一种考虑新能源与负荷双重不确定性的储能调频控制方法,其特征在于,步骤S1中,所述考虑自动发电控制成本的综合成本最小的经济调度优化模型的目标函数,如下式所示:
式中:表示电池储能的退化成本;/>和/>分别表示第i台电池储能在第n时刻的充电功率和放电功率;ΕP表示求数学期望;/>表示传统电机的发电成本;/>为第j台电机的发电功率,/>表示第j台传统电机在n时刻的平均输出功率值;/>表示调频资源的里程成本;Mn表示调节里程;PFk,n表示第n时刻第k台设备的参与系数;S表示电池储能;G表示传统电机;{,S}表示所有调频资源,即传统电机和电池储能。
3.如权利要求1所述的一种考虑新能源与负荷双重不确定性的储能调频控制方法,其特征在于,步骤S1中,所述约束条件包括实时经济调度的约束条件、考虑AGC功率需求的机会约束以及储能的SOC约束。
4.如权利要求3所述的一种考虑新能源与负荷双重不确定性的储能调频控制方法,其特征在于,所述实时经济调度的约束条件包括:
功率平衡约束:
功率上下限约束:
功率爬坡约束:
线路容量约束:
式中,是第i台电池储能在第n时刻的输出功率;/>是节点b在n时刻负荷功率,b表示一个节点;Pk,min和Pk,max分别是传统电机或储能输出功率的上下限;/>表示第k个传统机组或者电池储能在n时刻的输出功率;rrk是调频资源的爬坡率;T表示一个经济调度周期;fl是线路l的容量;SFl,b是线路l上节点b的转移因子;/>分别表示节点b中传统机组、电池储能以及新能源发电的向量值;/>分别表示传统电机、电池储能以及新能源发电的输出功率;
考虑AGC功率需求的机会约束有:
响应AGC信号需求下,调频资源的功率输出在其自身约束范围之内:
式中:PFk,n是标号为k的调频资源的参与系数;P(·)表示该事件发生的概率;是一个经济调度周期中AGC信号特征值中的上分位数,/>是一个经济调度周期中AGC信号特征值中的下分位数;β表示规定的功率输出机会约束应满足的概率值;
同时,调频资源的爬坡率高于AGC功率需求的调节率:
式中,表示调节率上分位数,/>表示调节率下分位数,P(·)表示该事件发生的概率,β表示规定的爬坡率机会约束应满足的概率值;
所有调频资源的参与系数应介于0-1之间,且参与系数之和应为1,即
储能的SOC约束有:
式中,SOCi,n-1表示第i台电池储能在第n-1时刻的电池储能SOC值;是一个周期内AGC信号向上调节的累计值;/>是一个周期内AGC信号向下调节的累计值;CEi表示第i台电池储能的容量;/>和/>分别表示储能的充电效率和放电效率;SOCmax i和SOCmin i分别表示电池储能SOC值的上限和下限。
5.如权利要求4所述的一种考虑新能源与负荷双重不确定性的储能调频控制方法,其特征在于,所述步骤S2包括:
步骤S201、对AGC信号历史数据进行处理,提取累计调节能量累计调节里程Mn、调节幅度分位数/>以及调节率分位数/>然后用LSTM神经网络根据历史数据进行预测,并得到其概率分布;
步骤S202、采用Wasserstein度量从预测的AGC信号样本集中构建模糊集,从而找到一个最小化最坏情况下的决策;
步骤S203、根据步骤S1确立的目标函数及其约束条件,重新表述优化模型如下:
式中:表示为参与经济调度的所有储能的总退化成本;表示所有传统电机的总发电成本,EQ表示对于满足模糊集最差分布的求期望,Βki表示调频参与系数矩阵,/>表示AGC信号的特征值向量,Dki表示传统电机的发电成本;/>表示所有调频资源的调频里程成本,表示调频里程成本的参与系数矩阵;Mm表示机会约束系数矩阵,Nm表示调频资源输出功率相关的函数矩阵,βm表示第m个机会约束式子应满足的概率值,J表示机会约束的集合;
步骤S204、将AGC功率需求的机会约束转化为条件风险价值约束,随后将条件风险价值约束表述的约束条件乘以风险规避参数并添加到目标函数中;
步骤S205、用MATLAB中带的求解器进行求解,得到调频资源的基本调度点。
6.如权利要求4所述的一种考虑新能源与负荷双重不确定性的储能调频控制方法,其特征在于,所述步骤S3包括:。
建立单个储能调频耗量函数如下:
式中:DS,i(t)表示t时刻第i台电池储能的调频耗量,其为实时充放电功率和容量偏差的加权平方和,Pi S(t)表示第i台电池储能的实时充放电功率,表示第i台电池储能的容量偏差,CEi(t)表示t时刻电池储能的容量,/>表示电池储能的额定容量,aS,i和bS,i分别为系数,Δt表示调频控制间隔时间;
储能运行受到额定功率、额定容量和爬坡速率及功率平衡约束的限制,如下所示:
Pi,min≤Pi S≤Pi,max
Pi S(t)∈ΩP,i(t)∈{Pi S(t)|Pi SL(t)≤Pi S(t)≤Pi SU(t)}
式中:和/>表示储能的充放电爬坡速率范围;ΩP,i(t)表示电池储能输出的可行区间;Pi SL(t)和Pi SU(t)表示储能满足其运行约束的充放电功率上下限;
考虑储能实时容量受其充放电功率影响,将储能调频耗量函数变换为以储能充放电功率为单一变量的函数形式,如下所示:
DS,i(t)=αS,i[Pi S(t)]2+βS,iPi S(t)+γS,i
式中:αS,i、βS,i和γS,i分别是调频耗量函数二次项、一次项和常数项系数;
则建立AGC的优化模型如下:
式中:I为储能有效运行范围集合;PSref(t)表示t时刻电池储能需要满足的总功率需求。
7.如权利要求6所述的一种考虑新能源与负荷双重不确定性的储能调频控制方法,其特征在于,所述步骤S4包括:
步骤S401、引入拉格朗日乘子λ,将原优化问题描述为如下所示的新优化问题:
步骤S402、当储能充放电功率在其上、下界内取值时,集群内各储能调频边际成本趋于一致且等于λ*时,优化问题获得最优解,相应公式如下:
步骤S403、不断在相邻储能集群之间迭代一致性变量,直至达到一致,得到最优边际成本和最佳输出功率如下:
式中,βS,i(t)表示电池储能调频损耗函数一次项系数,αS,i表示电池储能调频损耗函数二次项系数,Ωi(t)表示电池储能输出功率可行域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310675994.6A CN116599086A (zh) | 2023-06-08 | 2023-06-08 | 一种考虑新能源发电与负荷不确定性的储能调频控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310675994.6A CN116599086A (zh) | 2023-06-08 | 2023-06-08 | 一种考虑新能源发电与负荷不确定性的储能调频控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116599086A true CN116599086A (zh) | 2023-08-15 |
Family
ID=87611658
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310675994.6A Pending CN116599086A (zh) | 2023-06-08 | 2023-06-08 | 一种考虑新能源发电与负荷不确定性的储能调频控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116599086A (zh) |
-
2023
- 2023-06-08 CN CN202310675994.6A patent/CN116599086A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109559035B (zh) | 一种考虑灵活性的城市配电网双层规划方法 | |
CN110676870B (zh) | 一种适用于风电并网的混合储能容量配置方法 | |
CN105680474B (zh) | 一种储能抑制光伏电站快速功率变化的控制方法 | |
CN109103912B (zh) | 考虑电网调峰需求的工业园区主动配电系统调度优化方法 | |
CN107069776B (zh) | 一种平滑微网联络线功率的储能前瞻分布式控制方法 | |
CN109599856B (zh) | 一种微网多楼宇中电动汽车充放电管理优化方法及装置 | |
CN113515884A (zh) | 分散式电动汽车实时优化调度方法、系统、终端及介质 | |
CN110311397B (zh) | 基于自适应识别的储能电站多场景协同控制方法 | |
CN112688347B (zh) | 一种电网负荷波动的平缓系统和方法 | |
CN112381375B (zh) | 一种基于潮流分配矩阵的电网经济运行域快速生成方法 | |
CN108009672B (zh) | 基于双层优化模型的水光互补电站日发电计划编制方法 | |
CN112993979A (zh) | 配电网无功优化方法、装置、电子设备和存储介质 | |
CN112053035B (zh) | 考虑经济性与灵活性的输电通道与储能联合规划方法 | |
CN113991661B (zh) | 一种考虑集群电动汽车内部充放电的微电网预测控制方法 | |
CN111641233A (zh) | 一种考虑新能源及负荷不确定性的电力系统日内灵活调峰方法 | |
CN114819373A (zh) | 一种基于合作博弈的共享混合储能电站的储能规划方法 | |
Li et al. | Fuzzy logic-based coordinated control method for multi-type battery energy storage systems | |
CN109149658B (zh) | 基于一致性理论的独立微电网分布式动态经济调度方法 | |
CN112202207A (zh) | 基于新能源出力随机特性的电网深度调峰运行优化方法 | |
CN111915161A (zh) | 基于可调鲁棒算法的区域光储充系统短期调度方法和系统 | |
CN115021295A (zh) | 针对一次调频的风电场混合储能容量优化配置方法及系统 | |
CN109768626B (zh) | 一种储能电站即插即用的能量流实现方法 | |
CN110808600A (zh) | 一种计算变电站内接入电池储能站的方法及系统 | |
CN108321916B (zh) | 一种具有能量协作功能的基站及能量协作方法 | |
CN114285086A (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 |