一种大规模光伏逆变器系统的分群等值建模方法
【专利摘要】本发明公开了一种大规模光伏逆变器系统的分群等值建模方法，包括：获取各光伏逆变器的预设控制参数的相对灵敏度；将各所述光伏逆变器的所述预设控制参数进行取值标准化；获取所述光伏逆变器的分群指标；根据K均值聚类算法对所述光伏逆变器的分群指标进行聚类分群；根据所述聚类分群的结果，分别将同群的光伏逆变器通过一个等值光伏逆变器来进行等值替换，进行等值建模。由于大规模光伏逆变器系统中的各光伏逆变器的控制参数存在差异，利用K均值聚类算法将光伏逆变器进行分群，使得同群的光伏逆变器具有类似的动态特性，采用一个等值光伏逆变器来等值替换同群中的光伏逆变器来进行建模，减小了光伏逆变器的建模中的误差，减少了建模的规模。
【专利说明】
一种大规模光伏逆变器系统的分群等值建模方法
技术领域
[0001] 本发明涉及光伏电网技术领域，特别是涉及一种大规模光伏逆变器系统的分群等 值建模方法。
【背景技术】
[0002] 随着科学技术的进步，大型并网光伏电站在我国得到迅速发展。但是，光伏发电因 其自身不同于传统热电厂的发电特性，使得光伏电站接入电网的安全稳定性运行面对着严 峻的挑战。
[0003] 光伏电站具有绿色环保等优点，对于大规模光伏电站接入电网后的电网运行特性 开展研究已经非常必要。对于光伏系统而言，光伏阵列为静态元件，光伏系统的动态特性主 要由光伏逆变器决定，因此，首要解决的问题是光伏逆变器的建模，然而，在建模过程中，若 对光伏电站中的每一个光伏逆变器进行建模，无疑会增大电力系统建模的规模，而且会引 起模型的有效性和数据修正等诸多问题，同时会增加仿真时间。因而，目前多采用等值的方 法来描述大规模光伏逆变器。
[0004] 目前对于光伏逆变器的等值主要有两种方法:第一种为单机等值法，即用一个光 伏逆变器来替代原先的光伏逆变器;第二种为多机等值法，即用若干个光伏逆变器代替原 先的光伏逆变器。其中，单机等值法主要用于各个光伏逆变器的控制参数差异不大，动态特 性基本相同的情况，但在实际的光伏电站中，各个光伏逆变器的动态特性存在差异，当光伏 逆变器的动态特性差距较大时，采用单机等值法与实际情况会产生较大的误差；多机等值 法在一定程度上可以克服单机等值法的精度不足的缺陷，但是缺少科学有效的分群指标， 因而影响了其精度。
[0005] 因而，如何提高大规模光伏逆变器等值模型的精度，是本领域技术人员目前需要 解决的技术问题。
【发明内容】
[0006] 本发明的目的是提供一种大规模光伏逆变器系统的分群等值建模方法，可以提高 大规模光伏逆变器等值模型的精度。
[0007] 为解决上述技术问题，本发明提供了如下技术方案：
[0008] -种大规模光伏逆变器系统的分群等值建模方法，包括：
[0009] 获取各光伏逆变器的预设控制参数的相对灵敏度；
[0010] 将各所述光伏逆变器的所述预设控制参数进行取值标准化；
[0011] 获取所述光伏逆变器的分群指标；
[0012] 根据K均值聚类算法对所述光伏逆变器的分群指标进行聚类分群；
[0013] 根据所述聚类分群的结果，分别将同群的光伏逆变器通过一个等值光伏逆变器来 进行等值替换，进行等值建模。
[0014] 优选地，所述获取各光伏逆变器的预设控制参数的相对灵敏度，包括：
[0015] 调节各所述光伏逆变器的所述预设控制参数，至所述光伏逆变器的公共连接点 PCC处的电流波形失真度G达到预设阈值；
[0016] 获取各所述光伏逆变器的各控制参数的绝对灵敏度；
[0017] 获取由所述光伏逆变器的各控制参数的波动而导致的电流波形失真度G的总波动 参数；
[0018] 获取各所述光伏逆变器各控制参数的相对灵敏度。
[0019] 优选地，所述将各所述光伏逆变器的所述预设控制参数进行取值标准化，包括：
[0020] 获取各所述光伏逆变器的各控制参数的取值标准参数：
[0022] 其中，i为不小于1的整数，t为不小于1的整数，p\为第i个光伏逆变器的第t个控制 参数，n为光伏逆变器的总个数，，〇(p t)为n个光伏逆变器第t个控制参数的标准差。
[0023] 优选地，所述获取所述光伏逆变器的分群指标，包括：
[0024] 获取第i个所述光伏逆变器的分群指标Xl;
[0025] 获取n个所述光伏逆变器的分群指标:X={X1，X2, . . .，Xn}，n为光伏逆变器的总个 数。
[0026] 优选地，所述根据K均值聚类算法对所述光伏逆变器的分群指标进行聚类分群，包 括：
[0027]将所述n个所述光伏逆变器的分群指标X分成c个群，l〈c〈n，c为整数；
[0028] 求取所述c个群的初始群中心:{vYv%，. . .，v°j，. . . V。}，其中，v°j表示第j个 初始群中心，K j<c;
[0029] 根据K均值聚类算法确定分群结果。
[0030] 优选地，所述求取所述c个群的初始群中心包括：
[0031 ]步骤SI 1:定义分群指标集合U，并初始化U为空集，j = 1;
[0032]步骤S12:计算任意两个所述光伏逆变器分群指标间的距离d，生成对应的距离矩 阵DnXn，其中，dab= (xa-xb)T(xa-xb)，Ka<n，Kb<n，dab表示第a个光伏逆变器分群指标与 第b个光伏逆变器分群指标间的距离；
[0033]步骤S13:根据所述距离矩阵DnXn，计算第i个所述光伏逆变器的分群指标与其它光 伏逆变器的分群指标间的距离之和U，其中，
[0034]步骤S14:求取n个所述光伏逆变器的分群指标与其它光伏逆变器的分群指标间的 距离之和，将距离之和最大的光伏逆变器的分群指标作为第j个初始群中心v%，并将该光伏 逆变器的分群指标加入已经使用的分群指标集合U中；
[0035]步骤S15:根据所述距离矩阵DnXn，筛选出与所述已经使用的分群指标集合U中的任 意一个光伏逆变器的分群指标间的距离均大于预设阀值a的所有光伏逆变器的分群指标， 并将其中最短距离最大的光伏逆变器的分群指标作为第j+1个初始群中心，将该光伏 逆变器的分群指标加入已经使用的分群指标集合U中，更新所述集合U;
[0036]步骤S16:判断j+l = c是否成立，若是，则判定已获得c个初始群中心，若否，则将j+ 1赋值给j，并返回步骤SI 5。
[0037] 优选地，所述根据K均值聚类算法确定分群结果，包括：
[0038] 步骤S21:定义迭代变量k，初始化k= 1，将所述c个群的初始群中心V13作为第k-1次 迭代的c个群中心V1- 1;
[0039]步骤S22:求取第i个光伏逆变器的分群指标与第k-1次迭代的第j个群中心间的距 离(1匕，￥1^1」），其中，(1匕，￥1^1」）=匕-￥ 1^1」)1^1-￥1^1」），求取第；[个光伏逆变器的分群指标 与第k-1次迭代的C个群中心间的距离，进而获得n个光伏逆变器的分群指标与第k-1次迭代 的C个群中心间的距离，并将每个光伏逆变器的分群指标归入距该光伏逆变器距离最小的 群中心对应的群中；
[0040] 步骤S23:统计第k-1次迭代的第j个群中心对应的群中包含的光伏逆变器的分群 指标数量为Skl，进而获得第k-1次迭代的c个群中心对应的群中包含的光伏逆变器的分群 指标的数量；
[0041] 步骤S24:更新第k次迭代的第j个群中心，其中，第k次迭代的第j个群中心的更新 计算方法为：
{Skl}表示第k-1次迭代的第j个群中心对应的群中包含的光 伏逆变器的分群指标对应的光伏逆变器的标号构成的集合；
[0042] 步骤S25:根据所述第k次迭代的第j个群中心的更新参数，获得第k次迭代的c个群 中心Vk= {vkl，vk2,…，Vkj，? ? ?，Vkc};
[0043]步骤S26:求取第k次迭代的目标函数：
[0045] 步骤S27:将k+1赋值给k，返回步骤S22执行，直到eSe1-1间的变化值小于允许误 差I，记录此时的k为k end，并将此时的c个群中心记
[0046] 优选地，根据所述聚类分群的结果，分别将同群的光伏逆变器通过一个等值光伏 逆变器来进行等值替换，进行等值建模包括：
[0047] 将第j个群中心中的x个参数除以#个光伏逆变器的第x个控制参数的相 对灵敏度的平均值，获得第j个等值光伏逆变器的第x个控制参数标准化后的值^并 计算第j个等值光伏逆变器的第x个控制参数其中，
[0048]获取第j个等值光伏逆变器的NP个控制参数，进而获得c个等值光伏逆变器各gNP 个控制参数，并通过C个等值光伏逆变分别器对C个群中的光伏逆变器进行等值替换，1 <NP，N P为预设控制参数中的控制参数的总个数。
[0049] 优选地，还包括：
[0050]判断所述等值建模的精度是否达到预设精度阈值；
[0051]若否，则重新进行等值建模。
[0052]优选地，所述判断所述等值建模的精度是否达到预设精度阈值包括：
[0053]获取所述等值建模的满意度指标e，其中，
[0055] tjPts。分别为详细模型和等值模型从动态过程开始到第一摆峰值的时间，A4PA S。 分别为详细模型和等值模型动态过程的第一摆峰值。
[0056] 与现有技术相比，上述技术方案具有以下优点：
[0057]本发明实施例所提供的一种大规模光伏逆变器系统的分群等值建模方法，包括： 获取各光伏逆变器的预设控制参数的相对灵敏度;将各所述光伏逆变器的所述预设控制参 数进行取值标准化;获取所述光伏逆变器的分群指标;根据K均值聚类算法对所述光伏逆变 器的分群指标进行聚类分群;根据所述聚类分群的结果，分别将同群的光伏逆变器通过一 个等值光伏逆变器来进行等值替换，进行等值建模。由于大规模光伏逆变器系统中的各光 伏逆变器的控制参数存在差异，利用K均值聚类算法将光伏逆变器进行分群，使得同群的光 伏逆变器具有类似的动态特性，然后采用一个等值光伏逆变器来等值替换同群中的光伏逆 变器来进行建模，减小了光伏逆变器的建模中的误差，同时减少了建模的规模。
【附图说明】
[0058] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明 的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据 这些附图获得其他的附图。
[0059] 图1为本发明一种【具体实施方式】所提供的大规模光伏逆变器系统的分群等值建模 方法流程图。
【具体实施方式】
[0060] 本发明的核心是提供一种大规模光伏逆变器系统的分群等值建模方法，可以提高 大规模光伏逆变器等值模型的精度。
[0061] 为了使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发 明的【具体实施方式】做详细的说明。
[0062] 在以下描述中阐述了具体细节以便于充分理解本发明。但是本发明能够以多种不 同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类 似推广。因此本发明不受下面公开的【具体实施方式】的限制。
[0063] 请参考图1，图1为本发明一种【具体实施方式】所提供的大规模光伏逆变器系统的分 群等值建模方法流程图。
[0064] 本发明的一种【具体实施方式】提供了一种大规模光伏逆变器系统的分群等值建模 方法，包括：
[0065] S1:获取各光伏逆变器的预设控制参数的相对灵敏度。
[0066] 在本实施方式中，该大规模光伏逆变器系统中的各光伏逆变器设置为相同的控制 策略，但是各光伏逆变器的控制参数及动态特性不同，且各光伏逆变器在公共连接点PCC处 汇流并接入电网，将光伏逆变器进行标号，标号依次为1~n，即本系统中包括n个光伏逆变 器，以构成大规模光伏逆变器的详细模型。
[0067]需要说明的是，在本文中所谓的大规模指的是光伏逆变器系统容量达到兆瓦级。
[0068] 在本发明的一种实施方式中，获取各光伏逆变器的预设控制参数的相对灵敏度， 包括：
[0069] S101:调节各光伏逆变器的预设控制参数，至光伏逆变器的公共连接点PCC处的电 流波形失真度G达到预设阈值，即调节n个光伏逆变器的预设控制参数，使得G达到最小。
[0070] S102:获取各光伏逆变器的各控制参数的绝对灵敏度。
[0071] 定义每个光伏逆变器均有Np个控制参数，即各光伏逆变器的预设控制参数包括Np 个控制参数。将第i个光伏逆变器的第t个控制参数pi t增大A pit，以求取第i个光伏逆变器 的第t个控制参数pit的绝对灵敏度，其中，第i个光伏逆变器的第t个控制参数的绝对灵敏度 为：
[0073] 其中，i为不小于1的整数，t为不小于1的整数，p\为第i个光伏逆变器的第t个控制 参数；
[0074] 进而进行不同的赋值，求取第i个光伏逆变器所有NP个控制参数的绝对灵敏度，再 进而根据不同的赋值，求取n个光伏逆变器各自N P个控制参数的绝对灵敏度， <NP〇
[0075] 在本实施方式中，采用取极限的方式获得光伏逆变器的控制参数的绝对灵敏度， 通过使每个控制参数增大相同的且尽可能小的比例求取，当然这个比例可由用户自行控 制。
[0076] S10 3:获取由光伏逆变器的各控制参数的波动而导致的电流波形失真度G的总波 动参数dG:
[0078]其中，n为光伏逆变器的总个数，Np为各光伏逆变器的预设控制参数中的控制参数 的个数。
[0079] S104:获取各光伏逆变器各控制参数的相对灵敏度。
[0080] 首先求取第i个光伏逆变器的第t个控制参数p\的相对灵敏度，进而进行赋值，求 取第i个光伏逆变器所有NP个控制参数的相对灵敏度，然后进行赋值求取n个光伏逆变器各 自N P个控制参数的相对灵敏度。
[0081 ]其中，第i个光伏逆变器的第t个控制参数pi的相对灵敏度为：
[0083]由于光伏逆变器为光伏系统的动态元件，同一控制模式下的光伏逆变器参数对不 同扰动下的动态响应特性有重要影响，而且不同参数的影响程度不同，本发明实施方式通 过上述的参数灵敏度的"微分摄动法"求取了各个光伏逆变器各个控制参数的相对灵敏度， 使得接下来的光伏逆变器的分群指标更加全面。
[0084] S2:将各光伏逆变器的预设控制参数进行取值标准化。
[0085] 在本发明的一种实施方式中，将各光伏逆变器的预设控制参数进行取值标准化， 包括:获取各光伏逆变器的各控制参数的取值标准参数。首先，将第i个光伏逆变器的第t个 控制参数取值标准化，进而将第i个光伏逆变器所有N P个控制参数取值标准化，最后将n 个光伏逆变器的各自NPf控制参数取值标准化。其中，第i个光伏逆变器的第t个控制参数 取值标准化的计算方法为：
[0087] 其中，i为不小于1的整数，t为不小于1的整数，p\为第i个光伏逆变器的第t个控制 参数，n为光伏逆变器的总个数，，〇(p t)为n个光伏逆变器第t个控制参数的标准差。
[0088] 取值标准化可以消除不同控制参数量纲的影响，同时可以避免不同控制参数间的 "大数吃小数"的现象。
[0089] S3:获取光伏逆变器的分群指标。
[0090] 在本发明的一种实施方式中，获取光伏逆变器的分群指标，包括：
[0091] 获取第i个光伏逆变器的分群指标Xl，每个光伏逆变器与它的分群指标是一一对 应的，其中，第i个光伏逆变器的分群指标 X1的计算方法为：
[0093] 即，根据第i个光伏逆变器的各控制参数的相对灵敏度以及该控制参数取值标准 化依据上式得出该光伏逆变器的分群指标。
[0094] 获取n个光伏逆变器的分群指标:X={X1，X2, . . .，Xn}，n为光伏逆变器的总个数。对 i进行赋值，依次使i取1~n，从而得出X。
[0095] S4:根据K均值聚类算法对光伏逆变器的分群指标进行聚类分群。
[0096] 在本发明的一种实施方式中，根据K均值聚类算法对光伏逆变器的分群指标进行 聚类分群，包括：
[0097] 将n个光伏逆变器的分群指标X分成c个群，l〈c〈n，c为整数;求取c个群的初始群中 心:VQ={v Qi，vQ2,. . .，v°i，. . .，VQC}，其中，V°j 表示第 j 个初始群中心，Kj彡C。
[0098] 其中，求取c个群的初始群中心包括：
[0099]步骤SI 1:定义已经使用的分群指标集合U，并初始化U为空集，并初始化j = 1，然后 根据需要依次对j进行赋值，直至j + l=c时，表示这C个初始群中心已经选择完毕；
[0100] 步骤S12:计算任意两个光伏逆变器分群指标间的距离d，从而生成对应的距离矩 阵DnXn，其中，d的计算方法为:dab = (xa-xb)T(xa-xb)，l<a<n，l<b<n，dab表示第a个光伏逆 变器分群指标与第b个光伏逆变器分群指标间的距离；
[0101] 步骤S13:根据距离矩阵DnXn，计算第i个光伏逆变器的分群指标与其它光伏逆变器 的分群指标间的距离之和U，其中，
[0102] 步骤S14:通过对i进行赋值，求取n个光伏逆变器的分群指标与其它光伏逆变器的 分群指标间的距离之和，筛选出距离之和最大的光伏逆变器的分群指标作为第j个初始群 中心V%，并将该光伏逆变器的分群指标加入已经使用的分群指标集合U中；
[0103] 步骤S15:根据距离矩阵DnXn，筛选出与已经使用的分群指标集合U中的任意一个光 伏逆变器的分群指标间的距离均大于预设阀值a的所有光伏逆变器的分群指标，并将其中 最短距离最大的光伏逆变器的分群指标作为第j+1个初始群中心，将该光伏逆变器的 分群指标加入已经使用的分群指标集合U中，从而更新集合U。
[0104] 步骤S16:判断j+l = c是否成立，若是，则判定已获得c个初始群中心，若否，则将j+ 1赋值给j，并返回步骤S15。
[0105] 本实施方式所提供的上述初始群中心选取原则可避免初始群中心选取的过近使 接下来的聚类算法的迭代次数增大，甚至陷入局部最优解的问题;若无法选取满足条件的c 个初始群中心，则可通过减少阈值a来进行解决。
[0106] 根据K均值聚类算法确定分群结果。
[0107] 在本发明的一种实施方式中，根据K均值聚类算法确定分群结果，包括：
[0108] 步骤S21:定义迭代变量k，初始化k=l，将c个群的初始群中心V13作为第k-1次迭代 的c个群中心V k、
[0109] 步骤S22:求取第i个光伏逆变器的分群指标与第k-1次迭代的第j个群中心间的距 离d(Xi，v1^ 1 j)，其中，d(Xi，v1^1 j) = (Xi-v1^1 j)T(Xi-v1^1 j)，进行不同的赋值，求取第i个光伏逆 变器的分群指标与第k_l次迭代的c个群中心间的距离，进而获得n个光伏逆变器的分群指 标与第k-1次迭代的c个群中心间的距离，并将每个光伏逆变器的分群指标归入距该光伏逆 变器距离最小的群中心对应的群中；
[0110] 步骤S23:统计第k-1次迭代的第j个群中心对应的群中包含的光伏逆变器的分群 指标数量为Skl，进而获得第k-1次迭代的c个群中心对应的群中包含的光伏逆变器的分群 指标的数量；
[0111] 步骤S24:更新第k次迭代的第j个群中心，其中，第k次迭代的第j个群中心的更新
计算方法为 {Skl}表示第k-1次迭代的第j个群中心对应的群中包含的光 .，： 伏逆变器的分群指标对应的光伏逆变器的标号构成的集合；
[0112] 步骤S25:根据第k次迭代的第j个群中心的更新参数，获得第k次迭代的c个群中心 Vk= {vkl,vk2, . . . ,Vkj, . . . ,Vkc}；
[0113] 步骤S26:求取第k次迭代的目标函数：
[0115]步骤S27:将k+1赋值给k，返回步骤S22执行，直到Ek与E1-1间的变化值小于允许误 差I，记录此时的k为kend，并将此时的c个群中心记
[0116]其中，聚类算法是一个迭代寻优的过程，就是使得分到同一群的指标间的相似度 较大，不同群的指标间的相似度较小，相似度的大小可用分群指标间的距离来衡量，采用K 均值聚类算法K-means作为本发明实施方式中的聚类算法，过程简单，几何意义明确，本发 明采用K均值聚类算法K-means来对光伏逆变器的分群指标进行聚类分群，提高了得到的结 果的可靠性。
[0117] 在本实施方式中，利用K均值聚类算法将n个逆变器分成c个群，每个群都有一个中 心，因此共c个群中心。在执行完K均值聚类算法后得到的c个群中心是最终要用的群中 心，也就是最终群中心。但是K均值聚类算法在聚类前需要提供c个群的初始群中心，即V% 有了这c个初始群中心V'K均值聚类算法才能开始聚类分群，利用提供的这c个初始群中心 V'得到最终的c个群的中心。
[0118] S5:根据聚类分群的结果，分别将同群的光伏逆变器通过一个等值光伏逆变器来 进行等值替换，进行等值建模。
[0119] 在本发明的一种实施方式中，根据聚类分群的结果，分别将同群的光伏逆变器通 过一个等值光伏逆变器来进行等值替换，进行等值建模包括：
[0120] 将第j个群中心中的x个参数除以个光伏逆变器的第x个控制参数的相 对灵敏度的平均值，获得第j个等值光伏逆变器的第x个控制参数标准化后的值并 计算第j个等值光伏逆变器的第x个控制参数Z
[0121]获取第j个等值光伏逆变器的NP个控制参数，进而获得c个等值光伏逆变器各自NP 个控制参数，并通过C个等值光伏逆变分别器对C个群中的光伏逆变器进行等值替换，1 <NP，N P为预设控制参数中的控制参数的总个数。
[0122] 通过K均值聚类算法K-means得到的各个群中心，可认为是该群中光伏逆变器的分 群指标的典型值，光伏逆变器控制参数取值标准化的逆变换即为等值逆变器的参数。
[0123] 在上述任一实施方式的基础上，本发明一种实施方式所提供的方法，还包括：
[0124] 判断等值建模的精度是否达到预设精度阈值；
[0125] 其中，判断等值建模的精度是否达到预设精度阈值包括：
[0126] 获取等值建模的满意度指标e，其中，
[0128] tjPts。分别为详细模型和等值模型从动态过程开始到第一摆峰值的时间，A4PA S。 分别为详细模型和等值模型动态过程的第一摆峰值。
[0129] 若否，则重新进行等值建模。
[0130] 在本实施方式中，按照上式可以得出等值模型的满意度指标，以对大规模光伏逆 变器的等值模型进行评价，其中，得出的满意度指标越大则表示等值模型的精度越高，其 中，t s。反应了控制环节的时间参数的影响，As。反应了其他参数的影响，满意度达到所要求 的标准，即达到预设精度阈值，则认为等值方案是合理的，通常可取等值方案中等值光伏逆 变器个数最少的方案作为选用的方案。
[0131] 综上所述，本发明实施方式所提供的一种大规模光伏逆变器系统的分群等值建模 方法，根据不同光伏逆变器的控制参数存在差异的问题，利用K均值聚类算法将光伏逆变器 进行分群，使得同群的光伏逆变器具有类似的动态特征，实现了大规模光伏逆变器的等值 化简，从而提高了等值建模的精度。
[0132] 其次，通过聚类算法得到的群中心来方向求取等值逆变器的控制参数，无需繁琐 的公式推导，计算过程更加简单方便，具有较高的通用性。
[0133] 以上对本发明所提供一种大规模光伏逆变器系统的分群等值建模方法进行了详 细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说 明只是用于帮助理解本发明及其核心思想。应当指出，对于本技术领域的普通技术人员来 说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰 也落入本发明权利要求的保护范围内。
【主权项】
1. 一种大规模光伏逆变器系统的分群等值建模方法，其特征在于，包括： 获取各光伏逆变器的预设控制参数的相对灵敏度； 将各所述光伏逆变器的所述预设控制参数进行取值标准化； 获取所述光伏逆变器的分群指标； 根据K均值聚类算法对所述光伏逆变器的分群指标进行聚类分群； 根据所述聚类分群的结果，分别将同群的光伏逆变器通过一个等值光伏逆变器来进行 等值替换，进行等值建模。2. 根据权利要求1所述的方法，其特征在于，所述获取各光伏逆变器的预设控制参数的 相对灵敏度，包括： 调节各所述光伏逆变器的所述预设控制参数，至所述光伏逆变器的公共连接点PCC处 的电流波形失真度G达到预设阈值； 获取各所述光伏逆变器的各控制参数的绝对灵敏度； 获取由所述光伏逆变器的各控制参数的波动而导致的电流波形失真度G的总波动参 数； 获取各所述光伏逆变器各控制参数的相对灵敏度。3. 根据权利要求2所述的方法，其特征在于，所述将各所述光伏逆变器的所述预设控制 参数进行取值标准化，包括： 获取各所述光伏逆变器的各控制参数的取值标准参数：其中，i为不小于1的整数，t为不小于1的整数，p\为第i个光伏逆变器的第t个控制参 数，n为光伏逆变器的总个数，，〇(pt)为n个光伏逆变器第t个控制参数的标准差。4. 根据权利要求3所述的方法，其特征在于，所述获取所述光伏逆变器的分群指标，包 括： 获取第i个所述光伏逆变器的分群指标Xl; 获取n个所述光伏逆变器的分群指标:X={X1，X2,. . .，Xn}，n为光伏逆变器的总个数。5. 根据权利要求4所述的方法，其特征在于，所述根据K均值聚类算法对所述光伏逆变 器的分群指标进行聚类分群，包括： 将所述n个所述光伏逆变器的分群指标X分成c个群，l〈c〈n，c为整数； 求取所述c个群的初始群中心:. . .，v°j，. . . ，其中，v°j表示第j个初始 群中心，K j<c; 根据K均值聚类算法确定分群结果。6. 根据权利要求5所述的方法，其特征在于，所述求取所述c个群的初始群中心包括： 步骤S11:定义分群指标集合U，并初始化U为空集，j = 1; 步骤S12:计算任意两个所述光伏逆变器分群指标间的距离d，生成对应的距离矩阵 DnXn，其中，dab= (xa-xb)T(xa-xb)，Ka<n，Kb<n，dab表示第a个光伏逆变器分群指标与第 b个光伏逆变器分群指标间的距离； 步骤S13:根据所述距离矩阵DnXn，计算第i个所述光伏逆变器的分群指标与其它光伏逆 变器的分群指标间的距离之和U，其中步骤S14:求取n个所述光伏逆变器的分群指标与其它光伏逆变器的分群指标间的距离 之和，将距离之和最大的光伏逆变器的分群指标作为第j个初始群中心v%，并将该光伏逆变 器的分群指标加入已经使用的分群指标集合U中； 步骤S15:根据所述距离矩阵DnXn，筛选出与所述已经使用的分群指标集合U中的任意一 个光伏逆变器的分群指标间的距离均大于预设阀值a的所有光伏逆变器的分群指标，并将 其中最短距离最大的光伏逆变器的分群指标作为第j+1个初始群中心，将该光伏逆变 器的分群指标加入已经使用的分群指标集合U中，更新所述集合U; 步骤S16:判断j+l = c是否成立，若是，则判定已获得c个初始群中心，若否，则将j+1赋 值给j，并返回步骤S15。7. 根据权利要求6所述的方法，其特征在于，所述根据K均值聚类算法确定分群结果，包 括： 步骤S21:定义迭代变量k，初始化k= 1，将所述c个群的初始群中心V13作为第k-1次迭代 的c个群中心Vk、 步骤S22:求取第i个光伏逆变器的分群指标与第k-1次迭代的第j个群中心间的距离d (义1，￥1^1」），其中，(1(11，￥1^1」）=(11-￥ 1^1」)1'(11-￥1^1」），求取第；[个光伏逆变器的分群指标与第 k-1次迭代的C个群中心间的距离，进而获得n个光伏逆变器的分群指标与第k-1次迭代的C 个群中心间的距离，并将每个光伏逆变器的分群指标归入距该光伏逆变器距离最小的群中 心对应的群中； 步骤S23:统计第k-1次迭代的第j个群中心对应的群中包含的光伏逆变器的分群指标 数量为Sk、，进而获得第k-1次迭代的c个群中心对应的群中包含的光伏逆变器的分群指标 的数量； 步骤S24:更新第k次迭代的第j个群中心，其中，第k次迭代的第j个群中心的更新计算 方法为{Skl}表示第k-1次迭代的第j个群中心对应的群中包含的光伏逆 变器的分群指标对应的光伏逆变器的标号构成的集合； 步骤S25:根据所述第k次迭代的第j个群中心的更新参数，获得第k次迭代的c个群中心 Vk= {vkl,vk2, . . . ,Vkj, . . . ,Vkc}； 步骤S26:求取第k次迭代的目标函数：步骤S27:将k+1赋值给k，返回步骤S22执行，直到Ek与E1-1间的变化值小于允许误差记 录此时的k为kend，并将此时的C个群中心记为。8. 根据权利要求7所述的方法，其特征在于，根据所述聚类分群的结果，分别将同群的 光伏逆变器通过一个等值光伏逆变器来进行等值替换，进行等值建模包括： 将第j个群中心V4"%中的x个参数除以没^^1,个光伏逆变器的第x个控制参数的相对灵 敏度的平均值，获得第j个等值光伏逆变器的第x个控制参数标准化后的机并计算 第j个等值光伏逆变器的第X个控制参数其中获取第j个等值光伏逆变器的NP个控制参数，进而获得C个等值光伏逆变器各自NP个控 制参数，并通过c个等值光伏逆变分别器对c个群中的光伏逆变器进行等值替换，l<x<NP， NP为预设控制参数中的控制参数的总个数。9. 根据权利要求1至8任一项所述的方法，其特征在于，还包括： 判断所述等值建模的精度是否达到预设精度阈值； 若否，则重新进行等值建模。10. 根据权利要求9所述的方法，其特征在于，所述判断所述等值建模的精度是否达到 预设精度阈值包括： 获取所述等值建模的满意度指标e，其中，t4Pts。分别为详细模型和等值模型从动态过程开始到第一摆峰值的时间，A4PAS。分别 为详细模型和等值模型动态过程的第一摆峰值。
【文档编号】G05B17/02GK106054665SQ201610367326
【公开日】2016年10月26日
【申请日】2016年5月27日
【发明人】赵波, 李鹏, 吴红斌, 刘众前, 张雪松, 徐琛
【申请人】国网浙江省电力公司电力科学研究院, 国家电网公司