一种基于最小互信息的光伏出力分解方法技术

本发明专利技术公开了一种基于最小互信息的光伏出力分解方法，包括以下步骤：建立光伏负荷混合模型：将台区内的负荷分为光伏负荷与非光伏负荷两类；构建光伏出力分解算法；基于卡尔曼滤波的出力估计：即采用卡尔曼滤波算法对光伏倍率进行最优无偏估计，对光伏倍率进行状态估计，光伏负荷混合模型采用一维卡尔曼滤波，并将其分为预测与校正两个过程；引入光伏倍率动态变量a，计算得出t时刻实时出力表达式。本发明专利技术仅基于台区总表的量测与光伏出力曲线轮廓即可实现算法分解，无需安装额外量测，算法易于推广；算法从负荷与光伏之间的物理特性出发，当负荷特性发生变化时仍具有较好的精度，具有一定的普适性。具有一定的普适性。具有一定的普适性。

【技术实现步骤摘要】
一种基于最小互信息的光伏出力分解方法


[0001]本专利技术涉及一种基于最小互信息的光伏出力分解方法，属于分布式光伏监测


技术介绍

[0002]在分布式光伏项目推进的背景下，台区内光伏渗透率将不断提高，仅2021年的分布式光伏容量就新增2928万千瓦，占全国新增装机的16.6％。出于电量平衡、配电网安全运行、容量管理的需求，电网迫切需要对台区内分布式光伏进行监测。但是，一方面台区内光伏信息标准缺失，对信息采集与通信的成本敏感，另一方面维护水平低导致监测数据质量较差，不能满足光伏监测需求。
[0003]盲信号分离是指在未知源信号与源信号混合方式的情况下，以观测的混合信号去估计源信号的技术，在信号处理、生物医学、数字水印、语音处理领域有众多应用。从分解原理上主要分为两类，一类以负熵、累计张量等高阶统计量为分离标准，一类采用深度学习技术实现在复杂混合环境下的分离。
[0004]在现有的关于台区内分布式光伏出力感知的相关研究方面，国内相关研究较少，而是研究分离光伏、负荷后进行更高精度的负荷预测，国外由于能源聚合商对负荷监测与窃电检测的需求，已经开展了相关研究，从实现手段上有两种：
[0005]1)根据台区内负荷特性进行光伏出力分解。有文献假设台区内各用户负荷特性一致，以少量量测完全的负荷特性采用聚类等方式进行分离。有文献采用非侵入式负荷监测的方式实现分解，但是需要安装额外的设备。文献采用相关性分析结合网格搜索方法，但是需要光伏规模大，无法在台区内有效应用。
[0006]2)采用深度学习的方式分离光伏出力。存在方法采用数据驱动的方式进行识别，方法需要积累一定的量测数据，实际应用困难。存在方法结合了聚类与组合神经网络的方法，但对于负荷的处理仍是静态的，面对台区负荷特性变化的情况仍存在限制。
[0007]那么，如何实现台区内分布式光伏出力的分解就成为本领域人员急需要解决的技术问题。

技术实现思路

[0008]目的：为了克服现有技术中存在的不足，本专利技术提供一种基于最小互信息的光伏出力分解方法。
[0009]技术方案：为解决上述技术问题，本专利技术采用的技术方案为：
[0010]一种基于最小互信息的光伏出力分解方法，包括如下步骤：
[0011]构建光伏负荷混合模型。
[0012]对光伏负荷混合模型进行分解获得评价函数，采用粒子群算法求取评价函数为最小值时的倍率λ
′
。
[0013]利用倍率λ
′
对光伏负荷混合模型进行卡尔曼滤波，得到估计值λ
t
。
[0014]根据λ
′
、λ
t
求取光伏负荷电度总曲线。
[0015]作为优选方案，所述光伏负荷混合模型计算公式如下：
[0016]A
T
(t)＝A
L
(t)+A
pv
(t)
[0017][0018]式中，A
T
(t)表示台区变压器总表的电度曲线，A
L
(t)、A
pv
(t)分别表示非光伏负荷电度总曲线与光伏负荷电度总曲线，表示光伏出力曲线轮廓，λ表示实际光伏总负荷与光伏出力曲线轮廓之间的倍率。
[0019]作为优选方案，对光伏负荷混合模型进行分解获得评价函数，采用粒子群算法求取评价函数为最小值时的倍率λ
′
的方法，包括如下步骤：
[0020]根据光伏负荷混合模型得到I(x,y)标准互信息表达式，表达式如下：
[0021][0022]式中，H(x)表示随机变量x的信息熵，随机变量x为光伏负荷电度总曲线，H(y)表示随机变量y的信息熵，随机变量y为非光伏负荷电度总曲线。H(x,y)表示随机变量x、y之间的联合熵。
[0023]将x、y用A
L
(t)、A
pv
(t)进行替换，得到评价函数的表达式如下：
[0024][0025]采用粒子群算法求取评价函数为最小值时的倍率λ
′
。
[0026]作为优选方案，采用粒子群算法求取评价函数为最小值时的倍率λ
′
的方法，包括如下步骤：
[0027]初始化多组粒子(e
i
,f
i
),e
i
表示待选λ，f
i
表示该粒子的速度，i代表粒子群中第i个粒子；
[0028]迭代更新粒子速度与位置，得出更新后的e
i
′
；
[0029]f
i
′
＝f
i
+c1*rand()*(pbest
‑
e
i
)+c2*rand()*(gbest
‑
e
i
)
ꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0030]e
i
′
＝e
i
+f
i
′ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0031]式中，c1、c2表示学习因子，pbest为本次迭代中使评价函数最小对应的e
i
，gbest表示全局迭代中使评价函数最小对应的e
i
，rand()表示取值为0
‑
1的随机数；*代表乘号，f
i
′
为更新后的粒子速度，e
i
′
为更新后的粒子位置。
[0032]达到最大迭代次数，e
i
′
输出最优倍率λ
′
。
[0033]作为优选方案，利用倍率λ
′
对光伏负荷混合模型进行卡尔曼滤波，得到估计值λ
t
的方法，包括如下步骤：
[0034]获取t时刻的倍率预测值
[0035][0036]其中，λ
t
‑1为t
‑
1时刻经过卡尔曼滤波的倍率。
[0037]获取t时刻的误差矩阵预测值P
t
'。
[0038]P
t
'＝P
t
‑1+Q
[0039]其中，P
t
‑1为t
‑
1时刻的误差矩阵，Q为预测过程的噪声方差。
[0040]获取卡尔曼增益K
t
。
[0041][0042]其中，R为倍率λ
′
的噪声方差。
[0043]根据K
t
、λ
′
计算估计值λ
t
。
[0044][0045]作为优选方案，根据λ
′
、λ
t
求取光伏负荷电度总曲线的方法，包括如下步骤：
[0046]将λ
′...

【技术保护点】

【技术特征摘要】
1.一种基于最小互信息的光伏出力分解方法，其特征在于：包括如下步骤：构建光伏负荷混合模型；对光伏负荷混合模型进行分解获得评价函数，采用粒子群算法求取评价函数为最小值时的倍率λ
′
；利用倍率λ
′
对光伏负荷混合模型进行卡尔曼滤波，得到估计值λ
t
；根据λ
′
、λ
t
求取光伏负荷电度总曲线。2.根据权利要求1所述的一种基于最小互信息的光伏出力分解方法，其特征在于：所述光伏负荷混合模型计算公式如下：A
T
(t)＝A
L
(t)+A
pv
(t)式中，A
T
(t)表示台区变压器总表的电度曲线，A
L
(t)、A
pv
(t)分别表示非光伏负荷电度总曲线与光伏负荷电度总曲线，表示光伏出力曲线轮廓，λ表示实际光伏总负荷与光伏出力曲线轮廓之间的倍率。3.根据权利要求1所述的一种基于最小互信息的光伏出力分解方法，其特征在于：对光伏负荷混合模型进行分解获得评价函数，采用粒子群算法求取评价函数为最小值时的倍率λ
′
的方法，包括如下步骤：根据光伏负荷混合模型得到I(x,y)标准互信息表达式，表达式如下：式中，H(x)表示随机变量x的信息熵，随机变量x为光伏负荷电度总曲线，H(y)表示随机变量y的信息熵，随机变量y为非光伏负荷电度总曲线；H(x,y)表示随机变量x、y之间的联合熵；将x、y用A
L
(t)、A
pv
(t)进行替换，得到评价函数的表达式如下：采用粒子群算法求取评价函数为最小值时的倍率λ
′
。4.根据权利要求3所述的一种基于最小互信息的光伏出力分解方法，其特征在于：采用粒子群算法求取评价函数为最小值时的倍率λ
′
的方法，包括如下步骤：初始化多组粒子(e
i
,f
i
),e
i
表示待选λ，f
i
表示该粒子的速度，i代表粒子群中第i个粒子；迭代更新粒子速度与位置，得出更新后的e
i
′
；f
i
′
＝f
i
+c1*rand()*(pbest

【专利技术属性】
技术研发人员：卜强生，吕朋蓬，窦晓波，罗飞，周旭峰，彭志强，张潼，郑明忠，叶志刚，张岩，袁洲茂，
申请(专利权)人：东南大学国网江苏省电力有限公司江苏省电力试验研究院有限公司，
类型：发明
国别省市：