一种面向省内新能源消纳的储能电站容量规划方法技术

一种面向省内新能源消纳的储能电站容量规划方法，涉及储能应用技术领域，该方法包括如下步骤：步骤一：建立新能源发电出力的概率分布函数集合；步骤二：建立鲁棒机会约束规划模型；步骤三：求解分布鲁棒规划模型。本发明专利技术的有益效果：本发明专利技术为提高新能源消纳提供了一种储能电站容量规划方法，该方法以降低系统年新能源弃电率为目标，提出一种采用分布鲁棒优化方法的新能源储能装置容量优化配置方案，与传统的处理不确定因素的随机规划和鲁棒优化方法相比具有显著的优势。法相比具有显著的优势。法相比具有显著的优势。

【技术实现步骤摘要】
一种面向省内新能源消纳的储能电站容量规划方法


[0001]本专利技术涉及储能应用
，特别是涉及一种面向省内新能源消纳的储能电站容量规划方法。

技术介绍

[0002]新能源出力呈现出间歇性和波动性的特点，给电力系统的功率平衡带来巨大压力，成为高比例新能源消纳的主要障碍。储能系统可突破传统电力在时间与空间上的供需约束，具备了参与电力系统运行调控与安全稳定控制的能力。对于新能源比例较高的电力系统而言，配置一定的储能设备具有显著的必要性。为了保证储能灵活调节作用的有效发挥，提高新能源消纳率，有必要对新能源储能装置容量优化配置进行深入研究。

技术实现思路

[0003]本专利技术通过提出一种面向省内新能源消纳的储能电站容量规划方法，采用分布鲁棒优化方法研究新能源储能装置容量优化配置达到降低系统年新能源弃电率的目的。
[0004]针对以上现有技术的不足，本专利技术提供一种面向省内新能源消纳的储能电站容量规划方法，该方法包括如下步骤：
[0005]步骤一：建立新能源发电出力的概率分布函数集合。根据历史数据构造新能源发电出力的经验分布，以Kullback
‑
Leibler(KL)散度作为分布函数距离测度建立新能源发电出力的概率分布函数集合。
[0006]步骤二：建立鲁棒机会约束规划模型。以储能投资成本最小为目标，以年新能源弃电率为约束建立鲁棒机会约束规划模型。
[0007]步骤三：求解分布鲁棒规划模型。通过矫正机会约束中的风险阈值将鲁棒机会约束转化为传统机会约束，并借助凸近似和抽样平均构建线性规划进行高效求解。
[0008]进一步地，采用分布鲁棒优化的方法研究新能源的储能容量规划降低系统年新能源弃电率。
[0009]该专利技术的具体方法如下：
[0010]步骤一：建立新能源发电出力的概率分布函数集合
[0011]1)直流潮流模型
[0012]由于储能装置主要对有功进行调控，故采用直流潮流对网络进行建模。
[0013][0014]P
ij,t
＝(θ
i,t
‑
θ
j,t
)/x
ij
ꢀꢀꢀ
(2)
[0015]其中，P
i,t
表示节点i的注入有功功率之和(机组出力减去负荷)；P
ij,t
表示线路(i,j)的有功潮流；θ
i,t
表示节点i的电压相角；x
ij
表示线路(i,j)的电抗值。
[0016]2)储能模型
[0017]储能装置在电网中发挥削峰填谷的作用，当新能源发电出力过剩时，储能装置充
电；当系统新能源发电出力不足时，储能装置放电。储能装置充放电过程的数学模型如下：
[0018][0019]其中，表示储能装置在时刻t的储能量；表示储能装置的损失率；和表示储能装置的充电和放电功率；和表示充电和放电效率；P
cmin
和P
cmax
表示储能装置的最小和最大充电功率；和表示储能装置的最小和最大放电功率；和表示储能装置的最小和最大储能量。
[0020]3)确定性储能容量规划模型
[0021]考虑储能装置的规划问题，网络架构、传统火电机组容量位置和新能源容量位置均为给定值。
[0022]在系统中无储能装置时，由于系统中只含有传统火电机组和新能源发电机组，调节能力不足，造成了大量新能源弃电。为了降低新能源弃电率，采用在新能源节点配置储能装置的手段，依靠储能装置削峰填谷的能力以提升系统灵活性，减小新能源弃电率。因此，在不考虑新能源发电不确定性情况下的储能规划模型如下：
[0023][0024]其中，目标函数旨在最小化投资成本，I
E
是储能装置单位投资造价，是节点i储能装置的容量；Cons
‑
PF表示电网潮流方程；Cons
‑
EES表示储能装置充放电约束(3)；Cons
‑
BD表示潮流变量上下界约束；D
curt
≤R
curt
则规定实际新能源弃电率必须小于给定阈值R
curt
。
[0025]实际新能源弃电率D
curt
定义如下：从春、夏、秋、冬四个季节中各取一个典型日，调度时间间隔为1小时，则以此96个点代表全年的新能源发电出力情况，则D
curt
可表述为
[0026][0027]其中，表示t时刻新能源i的最大可发电量；表示t时刻新能源i的实际发电量。
[0028]在不考虑新能源发电不确定性的情况下，新能源储能规划模型是一个线性规划模型，可以高效求解。然而，新能源发电出力具有随机性和波动性，在系统规划这种长期优化决策问题中如果忽略不确定性的影响，可能会导致最终的结果无法满足实际系统对新能源弃电率的要求。因此，采用了基于KL散度的分布鲁棒优化方法处理新能源发电的不确定性。
[0029]4)不确定集合建模
[0030]①
生成参考分布P0：
[0031]目前最常用的生成参考分布的方法是利用历史数据进行估计。例如，假设有M个抽样可以分类到N个区间中，则在每个区间中有M1,M2,
…
,MN个抽样。每个区间中代表样本是区间中样本的期望值，对应的概率是π1＝M
i
/M,i＝1,...,N，则参考分布P0为{π1,...,π
N
}。此外，也可以假设不确定因素符合某一特定分布，如高斯分布等，从而利用参数估计的方法确定分布函数。
[0032]②
建立不确定集合
[0033]首先，采用KL散度描述两个分布函数之间的距离，距离越小，表明两个分布越相似。对于连续型的分布其定义如下：
[0034][0035]对于离散型的分布其定义如下：
[0036][0037]基于KL散度的描述，考虑了与参考分布P0的KL距离不超过d
KL
的所有分布函数，从而构建如下的不确定集合(集合中的元素为分布函数)：
[0038]W＝{P|D
KL
(P||P0)≤d
KL
}
ꢀꢀꢀ
(8)
[0039]当d
KL
＞0时，不确定集合W中含有无穷多个分布函数；随着d
KL
趋近于0，W变成单元素P0，后续描述的分布鲁棒规划模型也转变为一个传统的随机规划模型。
[0040]③
选择集合距离d
KL
[0041]在实际决策中，决策制定者需要根据风险偏好决定d
KL
的大小。显然，历史数据越多，则估计出来的参考分布与真实分布越近，可以设定更小的d
KL
值。根据分布鲁棒优化理论，d
KL
可以采用如下的选取方法：
[0042][0043]其中，代表N
‑...

【技术保护点】

【技术特征摘要】
1.一种面向省内新能源消纳的储能电站容量规划方法，该方法包括如下步骤：步骤一：建立新能源发电出力的概率分布函数集合，根据历史数据构造新能源发电出力的经验分布，以Kullback
‑
Leibler(KL)散度作为分布函数距离测度建立新能源发电出力的概率分布函数集合；(一)直流潮流模型由于储能装置主要对有功进行调控，故采用直流潮流对网络进行建模，P
ij,t
＝(θ
i,t
‑
θ
j,t
)/x
ij (2)其中，P
i,t
表示节点i的注入有功功率之和(机组出力减去负荷)；P
ij,t
表示线路(i,j)的有功潮流；θ
i,t
表示节点i的电压相角；x
ij
表示线路(i,j)的电抗值；(二)储能模型储能装置在电网中发挥削峰填谷的作用，当新能源发电出力过剩时，储能装置充电；当系统新能源发电出力不足时，储能装置放电，储能装置充放电过程的数学模型如下：其中，表示储能装置在时刻t的储能量；表示储能装置的损失率；和表示储能装置的充电和放电功率；和表示充电和放电效率；和表示储能装置的最小和最大充电功率；和表示储能装置的最小和最大放电功率；和表示储能装置的最小和最大储能量；(三)确定性储能容量规划模型考虑储能装置的规划问题，网络架构、传统火电机组容量位置和新能源容量位置均为给定值，在系统中无储能装置时，由于系统中只含有传统火电机组和新能源发电机组，调节能力不足，造成了大量新能源弃电，为了降低新能源弃电率，采用在新能源节点配置储能装置的手段，依靠储能装置削峰填谷的能力以提升系统灵活性，减小新能源弃电率，因此，在不考虑新能源发电不确定性情况下的储能规划模型如下：其中，目标函数旨在最小化投资成本，I
E
是储能装置单位投资造价，是节点i储能装置的容量；Cons
‑
PF表示电网潮流方程；Cons
‑
EES表示储能装置充放电约束(3)；Cons
‑
BD表示潮流变量上下界约束；D
curt
≤R
curt
则规定实际新能源弃电率必须小于给定阈值R
curt
，实际
新能源弃电率D
curt
定义如下：从春、夏、秋、冬四个季节中各取一个典型日，调度时间间隔为1小时，则以此96个点代表全年的新能源发电出力情况，则D
curt
可表述为其中，表示t时刻新能源i的最大可发电量；表示t时刻新能源i的实际发电量，在不考虑新能源发电不确定性的情况下，新能源储能规划模型是一个线性规划模型，可以高效求解，然而，新能源发电出力具有随机性和波动性，在系统规划这种长期优化决策问题中如果忽略不确定性的影响，可能会导致最终的结果无法满足实际系统对新能源弃电率的要求，因此，采用了基于KL散度的分布鲁棒优化方法处理新能源发电的不确定性；(四)不确定集合建模1、生成参考分布P0：目前最常用的生成参考分布的方法是利用历史数据进行估计。例如，假设有M个抽样可以分类到N个区间中，则在每个区间中有M1,M2,
…
,MN个抽样。每个区间中代表样本是区间中样本的期望值，对应的概率是π1＝M
i
/M,i＝1,...,N，则参考分布P0为{π1,...,π
N
}。此外，也可以假设不确定因素符合某一特定分布，如高斯分布等，从而利用参数估计的方法确定分布函数；2、建立不确定集合首先，采用KL散度描述两个分布函数之间的距离，距离越小，表明两个分布越相似，对于连续型的分布其定义如下：对于离散型的分布其定义如下：基于KL散度的描述，考虑了与参考分布P0的KL距离不超过d
KL
的所有分布函数，从而构建如下的不确定集合(集合中的元素为分布函数)：W＝{...

【专利技术属性】
技术研发人员：刘庭响，李春来，杨立滨，李正曦，安娜，周万鹏，马俊雄，王恺，高金，
申请(专利权)人：国网青海省电力公司经济技术研究院国网青海省电力公司国家电网有限公司，
类型：发明
国别省市：