采用高斯混合模型的风电场最优聚类等值方法
1.本发明属于电力系统新能源并网技术领域，具体涉及采用高斯混合模型的风电场最优聚类等值方法。
背景技术：
2.随着风电渗透率在电力系统中比例的不断攀升，双馈风电场的建设和实时仿真技术也日趋重要。由于实际风电场中风机数过多，采用详细的风电场模型进行仿真不仅耗费大量时间，甚至可能造成维数灾的发生。因此综合考虑风电场的实时仿真精度和仿真时间，有必要建立合适的双馈风电场等值模型。
3.由于风电场等值研究主要包括单机等值和多机等值两大类。由于风电场中各机组之间距离较长，线路参数也不尽相同，单机等值不能准确地描述整个风电场的动态特性。通常的多机等值可分为以下三个步骤：1)选取风机的实时参数或者某个时间序列电气量的特征信息得出不同的机划分指标；2)其次选取合适的聚类方法对机进行划分；3)分别计算不同聚类中的等值参数得到最终的聚类结果。
4.故障切除时风电机组机械特性变量、状态变量和转子电流受扰曲线等已经被用于风电机组分。较多文献研究了dfig输出有功功率的指标，采用模糊c均值算法对与机组相关的特征向量进行聚类分析。从避免出现聚类中心偏移方面来改进传统k-means算法，来实现风电场的动态等值多机表征。aic准则和bic准则被广泛应用于指标评价方面，可以考虑与风电场聚类结合。综上，具有良好等值效果且达到最优等值分数的等值方法亟需提出。研究风电场等值方法对于新能源系统经济稳定运行具有重要意义，并具有一定的工程与科研意义。
技术实现要素：
5.本发明的目的在于提供采用高斯混合模型的风电场最优聚类等值方法，解决了现有技术中针对于大型风电场难以得出最优等值分数的问题。
6.本发明所采用的技术方案是：采用高斯混合模型的风电场最优聚类等值方法，具体包括如下步骤：
7.步骤1，通过风电场总控获取场内不同位置风机的风速；
8.步骤2，通过matlab/simulink搭建多台含风电机组、变压器、集电线路物理参数的风电场模型，分别得到稳态有功功率和故障引起的电压跌落下的无功功率；
9.步骤3，选定高斯分布聚类分采用的风速、稳态有功功率和步骤2故障引起的电压跌落下的无功功率，计算每一个样本点数据由其对应簇生成的概率，多步迭代得到更新簇的概率直到高斯收敛；分别将每一个样本划分到概率最大的簇中，最终得到n类聚类分结果；
10.步骤4，根据步骤3的聚类分结果计算似然函数值，然后采用赤池信息准则aic和贝叶斯信息准则bic评判选取聚类中心，所述聚类中心对应的风速为内等效风速，由此获
得n个聚类；
11.步骤5，利用步骤4的聚类结果，在风电场等值计算中采用基于等值前后功率和电流相等的原则，分别将n个聚类内的风电机组、变压器和集电线路进行简化等值，最终得到具有最优分数的风电场简化等值结果。
12.本发明的特点还在于，
13.步骤2搭建的风电场模型包括多台1.5mw且具有低电压穿越能力的双馈风电机组，由0.69kv/35kv箱式变压器汇集至集电线路，经35kv/110kv升压变压器接入大电网系统。
14.步骤2中的风电场中集电线路参数复杂，并网点处电压因故障发生变化时反应到各单台风电机组上的电压变化是不同的。由于每台风机的风速不同，风机变桨距系统的桨距角和风机工作于启动区和恒转速恒功率区的状态也都有差异。双馈风机在电压跌落引起的故障中，由于定转子电流变化会引起电磁转矩的改变，进而会影响风机转子转速和输出的有功功率和无功功率。
15.高斯混合模型聚类是一种基于概率模型的方法，假设输入的样本服从k个参数未知的高斯分布，每一个高斯分布都对应不同的均值μj和协方差矩阵∑j(1≤j≤k)。那么服从同一分布的样本则被聚为一类。将步骤2的风速、稳态有功功率和无功功率作为样本集d＝{x1,x2,...,xm}，初始化高斯混合分布的参数混合系数αi、均值相量μi和协方差矩阵∑i，并合理假设每个混合元素都具有对角矩阵。其次计算每一个样本点xi(i＝1,2,3,
…
m)属于第j个高斯分布的概率：
[0016][0017]
根据式(2)、式(3)和式(4)计算更新迭代的模型参数值均值相量μi'、协方差矩阵∑i'和混合系数αi'。
[0018][0019][0020][0021]
根据将每一个样本划分到概率最大的簇中，依据三维的风速、稳态有功功率和无功功率样本集d＝{x1,x2,...,xm}进行三维聚类。得到最后的簇划分c＝{c1,c2,...,cn},最终得到n类聚类结果。
[0022]
步骤4中利用aic和bic准则进行判别是否得到最优分数：aic和bic定义分别如下所示：
[0023]
aic＝-2ln(l)+2k
ꢀꢀꢀ
(5)
[0024]
bic＝kln(n)-2ln(l)
ꢀꢀꢀ
(6)
[0025]
式中,k为模型的参数个数，n为样本数量，为似然函数；其中xj为第j个不同样本点，μ
l
为不同高斯分布对应的均值，∑l为不同高斯分布对应的协方差矩阵；结合风机三维聚类数据，在不同的聚类数下对应不同aic值和bic值，当两者都取最小，此时的风电场聚类数达到最优解。
[0026]
步骤5的计算过程如下：
[0027]
等值后第a风机的等值功率和等值电流可表示为：
[0028][0029]
等效后的第a风电机组的阻抗可表示为：
[0030][0031]
式(7)和式(8)中s
eq_a
(t)为内风电场的等值功率，i
eq_a
(t)为内风电场出口处的等值电流，z
eq_a
(t)为内风电场的等值阻抗，s
i(t)
为单独每一台风机对应的功率，i
i(t)
为单独每一台风机对应的出口处电流，v
i(t)
为单独每一台风机对应的出口处电压，v
eq_t
(t)为内风电场出口处电压，v
pcc
(t)为风电场并网点电压，i
ceq_a
(t)为内风电场的对地等效电容电流。
[0032]
考虑内每台风电机组的电容，第a风电机组对地电容可以近似表示为单台风机对地电容之和。
[0033][0034]
式(9)ci为第a内所有风机的对地电容。
[0035]
变压器参数容量采用内单台容量的倍乘，发电机内部参数和阻抗依据并联原则采用内单台阻抗的并联。第a内等值升压变压器容量和阻抗可表示为：
[0036]seq_a
＝as
t
ꢀꢀꢀ
(10)
[0037][0038]
式(10)和式(11)中，a为第a内风机的个数，s
t
为单台变压器的容量，z
t
为单台变压器的阻抗值。
[0039]
最后根据风电机组的聚类结果，计算出等值机的等效风速，以等效风速作为内各等值机的风速，最终得到等值机的各类参数。
[0040]
本发明的有益效果是，本发明的采用高斯混合模型的风电场最优聚类等值方法采用高斯混合模型聚类方法以及aic准则和bic准则进行判别得出最优的分数，解决了传统k-means算法需要给定初始聚类中心和聚类数的缺点，为风电场等值方法研究提出新的评价指标。计算原理简单且易于实现，使得风电场等值分结果更加符合实际，具有一定的工
程实用价值。
附图说明
[0041]
图1是本发明采用高斯混合模型的风电场最优聚类等值方法的风电场详细模型示意图；
[0042]
图2(a)是本发明采用高斯混合模型的风电场最优聚类等值方法的有功功率分指标；
[0043]
图2(b)是本发明采用高斯混合模型的风电场最优聚类等值方法的无功功率分指标；
[0044]
图3是本发明采用高斯混合模型的风电场最优聚类等值方法的等值效果流程图；
[0045]
图4是本发明采用高斯混合模型的风电场最优聚类等值方法的聚类效果图；xyz坐标分别代表风速、有功功率和无功功率，不同的聚类效果图采用不同的图形来表示；
[0046]
图5(a)是本发明采用高斯混合模型的风电场最优聚类等值方法的aic准则判别图；
[0047]
图5(b)是本发明采用高斯混合模型的风电场最优聚类等值方法的bic准则判别图；
[0048]
图6(a)是本发明采用高斯混合模型的风电场最优聚类等值方法在不同电压跌落度下的有功功率波形；
[0049]
图6(b)是本发明采用高斯混合模型的风电场最优聚类等值方法在不同电压跌落度下的无功功率波形。
具体实施方式
[0050]
下面结合附图和具体实施方式对本发明进行详细说明。
[0051]
步骤1，通过风电场总控获取场内不同位置风机的风速；
[0052]
步骤2，通过matlab/simulink搭建多台含风电机组、变压器、集电线路物理参数的风电场模型，分别得到稳态有功功率和故障引起的电压跌落下的无功功率；
[0053]
步骤2中的风电场中集电线路参数复杂，并网点处电压因故障发生变化时反应到各单台风电机组上的电压变化是不同的。由于每台风机的风速不同，风机变桨距系统的桨距角和风机工作于启动区和恒转速恒功率区的状态也都有差异。双馈风机在电压跌落引起的故障中，由于定转子电流变化会引起电磁转矩的改变，进而会影响风机转子转速和输出的有功功率和无功功率
[0054]
步骤3，选定高斯分布聚类分采用的风速、稳态有功功率和步骤2故障引起的电压跌落下的无功功率，计算每一个样本点数据由其对应簇生成的概率，多步迭代得到更新簇的概率直到高斯收敛；分别将每一个样本划分到概率最大的簇中，最终得到n类聚类分结果；
[0055]
步骤3中采用三维高斯混合模型聚类方法：
[0056]
高斯混合模型聚类是一种基于概率模型的方法，假设输入的样本服从k个参数未知的高斯分布，每一个高斯分布都对应不同的均值μj和协方差矩阵∑j(1≤j≤k)。那么服从同一分布的样本则被聚为一类。将步骤2的风速、稳态有功功率和无功功率作为样本集d
＝{x1,x2,...,xm}，初始化高斯混合分布的参数混合系数αi、均值相量μi和协方差矩阵∑i，并合理假设每个混合元素都具有对角矩阵。其次计算每一个样本点xi(i＝1,2,3,
…
m)属于第j个高斯分布的概率：
[0057][0058]
根据式(2)、式(3)和式(4)计算更新迭代的模型参数值均值相量μi'、协方差矩阵∑i'和混合系数αi'。
[0059][0060][0061][0062]
根据将每一个样本划分到概率最大的簇中，依据三维的风速、稳态有功功率和无功功率样本集d＝{x1,x2,...,xm}进行三维聚类。得到最后的簇划分c＝{c1,c2,...,cn},最终得到n类聚类结果。
[0063]
步骤4，根据步骤3的聚类分结果计算似然函数值，然后采用赤池信息准则aic和贝叶斯信息准则bic评判选取聚类中心，所述聚类中心对应的风速为内等效风速，由此获得n个聚类；
[0064]
步骤4中利用aic和bic准则进行判别是否得到最优分数：aic和bic定义分别如下所示：
[0065]
aic＝-2ln(l)+2k
ꢀꢀꢀ
(5)
[0066]
bic＝kln(n)-2ln(l)
ꢀꢀꢀ
(6)
[0067]
式中，k为模型参数个数，n为样本数量，为似然函数；其中xj为第j个不同样本点，μ
l
为不同高斯分布对应的均值，∑l为不同高斯分布对应的协方差矩阵；结合风机三维聚类数据，在不同的聚类数下对应不同aic值和bic值，当两者都取最小，此时的风电场聚类数达到最优解。
[0068]
步骤5，利用步骤4的聚类结果，在风电场等值计算中采用基于等值前后功率和电流相等的原则，分别将n个聚类内的风电机组、变压器和集电线路进行简化等值，最终得到具有最优分数的风电场简化等值结果。
[0069]
步骤5的计算过程如下：
[0070]
等值后第a风机的等值功率和等值电流可表示为：
[0071][0072]
等效后的第a风电机组的阻抗可表示为：
[0073][0074]
式(7)和式(8)中s
eq_a
(t)为内风电场的等值功率，i
eq_a
(t)为内风电场出口处的等值电流，z
eq_a
(t)为内风电场的等值阻抗，s
i(t)
为单独每一台风机对应的功率，i
i(t)
为单独每一台风机对应的出口处电流，v
i(t)
为单独每一台风机对应的出口处电压，v
eq_t
(t)为内风电场出口处电压，v
pcc
(t)为风电场并网点电压，i
ceq_a
(t)为内风电场的对地等效电容电流。
[0075]
考虑内每台风电机组的电容，第a风电机组对地电容可以近似表示为单台风机对地电容之和。
[0076][0077]
式(9)ci为第a内所有风机的对地电容。
[0078]
变压器参数容量采用内单台容量的倍乘，发电机内部参数和阻抗依据并联原则采用内单台阻抗的并联。第a内等值升压变压器容量和阻抗可表示为：
[0079]seq_a
＝as
t
ꢀꢀꢀ
(10)
[0080][0081]
式(10)和式(11)中，a为第a内风机的个数，s
t
为单台变压器的容量，z
t
为单台变压器的阻抗值。
[0082]
最后根据风电机组的聚类结果，计算出等值机的等效风速，以等效风速作为内各等值机的风速，最终得到等值机的各类参数。
[0083]
本发明采用高斯混合模型的风电场最优聚类等值方法原理为：以图1双馈风电场详细模型为例说明。
[0084]
双馈风电场由多台风机、变压器以及复杂集电线路组成，风力发电机发出电能经690v/35kv变压器输出汇集到35kv母线上，再经过35kv/110kv变压器升压并入大电网。若采用详细的风电场模型进行仿真不仅会耗费大量时间，甚至可能造成维数灾现象的发生，因此综合考虑风电场的实时仿真精度和仿真时间，进而建立合适的双馈风电场等值模型。
[0085]
(1)对风速、稳态期间有功功率和故障期间无功功率采用高斯混合模型聚类算法进行三维聚类，得到聚类结果图；
[0086]
(2)对于所得聚类分采用aic准则和bic准则进行判别，观测判别结果是否与聚类结果相对应，得出最终聚类结果；
[0087]
(3)采用等值前后功率和电流相等原则对变压器和集电参数进行等值计算，得出最终风电场等值模型。根据风电机组的聚类结果，计算出等值机的等效风速，以等效风速作
为内各等值机的风速，最终得到等值机的各类参数。
[0088]
实施例1
[0089]
以图1所述的双馈风电场详细模型为例说明，为了验证所提出的采用高斯混合模型聚类等值方法的性能，在matlab/simulink中建立了30台风机的详细模型。该模型中风机出口处690v电压经变压器升压至35kv，汇集到35kv交流母线上再经35kv/110kv变压器送至大电网，完成新能源风电场的并网。详细模型参考某含30台风机的风电场的变压器线路物理参数并合理修改以满足仿真要求。
[0090]
双馈风电场的等值步骤如下：
[0091]
1)通过风电场总控获取场内不同位置风机的不同风速；
[0092]
2)通过matlab/simulink搭建多台含风电机组、变压器、集电线路物理参数的风电场模型，分别得到稳态有功功率和故障引起的电压跌落下的无功功率；
[0093]
3)选定高斯分布聚类分采用的风速、稳态有功功率和步骤2故障引起的电压跌落下的无功功率，计算每一个样本点数据由其对应簇生成的概率，多步迭代得到更新簇的概率直到高斯收敛；分别将每一个样本划分到概率最大的簇中，最终得到n类聚类分结果；
[0094]
4)根据步骤3的聚类分结果计算似然函数值，然后采用赤池信息准则aic和贝叶斯信息准则bic评判选取聚类中心，所述聚类中心对应的风速为内等效风速，由此获得n个聚类；
[0095]
5)利用步骤4的聚类结果，在风电场等值计算中采用基于等值前后功率和电流相等的原则，分别将n个聚类内的风电机组、变压器和集电线路进行简化等值，最终得到具有最优分数的风电场简化等值结果。
[0096]
利用matlab/simulink进行仿真测试，针对不同电压跌落度情况进行等值精度验证。
[0097]
图6(a)是电压跌落度为0.5pu下详细模型、单机等值模型与本发明所提等值模型的有功功率和无功功率对比情况，图6(b)是电压跌落度为0.85pu下详细模型、单机等值模型与本发明所提等值模型的有功功率和无功功率对比情况，
[0098]
上述仿真结果皆证明本发明所提仿真精度高。
