基于运行数据的风电场参数辨识及动态等值方法.pdf

《基于运行数据的风电场参数辨识及动态等值方法.pdf》由会员分享，可在线阅读，更多相关《基于运行数据的风电场参数辨识及动态等值方法.pdf（19页完整版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 103887815 A (43)申请公布日 2014.06.25 CN 103887815 A (21)申请号 201410060760.1 (22)申请日 2014.02.21 H02J 3/38(2006.01) G06F 17/50(2006.01) (71)申请人 华南理工大学 地址 511458 广东省广州市南沙区环市大道 南 25 号华工大广州产研院 (72)发明人 赵卓立 杨苹 尹旭 许志荣 (74)专利代理机构 广州粤高专利商标代理有限 公司 44102 代理人 何淑珍 (54) 发明名称 基于运行数据的风电场参数辨识及动态等值 方法 (57) 摘要 。

2、本发明公开了基于运行数据的风电场参数辨 识及动态等值方法， 包括以下步骤 ： a ： 基于实测 数据对风电机组控制模型参数进行辨识 ； b ： 结合 运行数据， 选取能反映扰动下的风电机组和风电 场的特征变量， 应用改进的模糊 K 均值动态聚类 算法进行机群划分 ； c ： 基于全局最优位置变异粒 子群算法对风电场进行网络简化、 参数寻优， 获得 风电场动态等值模型 ； d ： 在扰动输入下， 比较风 电场动态等值模型与详细模型动态响应， 验证等 值模型的有效性。本发明建立的风场动态等值模 型能准确地反映风场并网点的动态特性， 具有重 要的工程应用价值， 可用于双馈型风电场接入电 力系统稳定性。

3、分析中， 为含风电场电力系统的规 划与运行调度提供理论支撑。 (51)Int.Cl. 权利要求书 2 页 说明书 10 页 附图 6 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书10页 附图6页 (10)申请公布号 CN 103887815 A CN 103887815 A 1/2 页 2 1. 基于运行数据的风电场参数辨识及动态等值方法， 其特征在于， 包括以下步骤 ： 步骤 a ： 基于实测数据对风电机组控制模型的功率和电流控制器参数进行参数辨识， 并对辨识结果进行通用性验证 ； 步骤 b ： 结合运行数据， 选取能反映扰动下的风电机组和风电场的特。

4、征变量， 应用改进 的模糊 K 均值动态聚类算法进行机群划分 ； 步骤 c ： 基于全局最优位置变异粒子群算法对风电场进行网络简化、 参数聚合， 获得风 电场动态等值模型 ； 步骤 d ： 在相同扰动输入下， 比较风电场动态等值模型与详细模型动态响应， 验证风电 场动态等值模型的有效性。 2. 如权利要求 1 所述的基于运行数据的风电场参数辨识及动态等值方法， 其特征在 于， 所述步骤 a 包括以下具体步骤 ： a1） 选择风机模型中功率控制器的控制参数 Kp、 Tp、 Kq、 Tq 和电流控制器的控制参数 Kd、 Td、 Kq、 Tq 作为需要重点辨识的控制器模型参数 ； 其中， Kp、 T。

5、p 分别是有功功率控制器 的比例常数和积分常数 ； Kq、 Tq 分别是无功功率控制器的比例常数和积分常数 ； Kd、 Td 分别 是有功电流控制器的比例常数和积分常数 ； Kq1、 Tq1 无功电流控制器的比例常数和积分常 数 ； a2） 搭建参数辨识模块， 所述参数辨识模块包括测量文件， 要辨识的模型元件、 参数辨 识比较器、 优化器和目标函数 ； a3） 设置参数初始值， 所述参数初始值为双馈风力发电机组经典控制模型参数 ； a4） 实测数据处理及导入， 所述实测数据为风电机组的输入风速 ； a5） 获得控制器参数辨识结果 ； a6） 基于控制器参数辨识结果， 进行模型验证。 3. 如权。

6、利要求 1 所述的基于运行数据的风电场参数辨识及动态等值方法， 其特征在 于， 所述步骤 b 包括以下具体步骤 ： b1） 运行数据选取风电场在扰动期间设定时间段内各风力发电机组的实测有功功率数 据， 并作为聚类分群的特征变量， 建立实测有功功率样本矩阵 ； b2） 确定聚类数目 K， 将全部数据随机分成 K 类， 其中， 聚类的簇的个数 K 值根据实际风 电场规模和模型精度要求选取 ； b3） 定义隶属度函数， 用隶属度函数定义的聚类损失函数 ； b4） 重复计算隶属度函数， 用当前的隶属度函数更新计算各类聚类中心直至算法收敛。 4. 如权利要求 1 所述的基于运行数据的风电场参数辨识及动态。

7、等值方法， 其特征在 于， 所述步骤 c 包括以下具体步骤 ： c1） 搭建双馈风电场等值模型， 所述风电场等值模型的功率控制器和电流控制器参数 采用双馈风机的单机模型经过辨识后的参数值， 机组参数按风机容量的变化作等比例增 加， 风轮半径保持不变 ； 定义目标函数 ：将大小扰动下实测有功功率数据作为功率参考， 将双 馈电机参数辨识问题转化为参数优化问题 ； c2） 选择待辨识参数 Rs、 Rr、 Ls、 Lr、 Xm， 初始化粒子群体的位置和速度， 设定种群大小， 设 权 利 要 求 书 CN 103887815 A 2 2/2 页 3 置粒子位置和速度上下限， 设定迭代次数、 收敛精度和相。

8、关权系数， 其中， Rs和 Ls分别为定 子电阻和定子漏抗 ； Rr和 Lr分别为归算后转子电阻和转子电抗 ； Xm为励磁电抗 ； c3） 计算粒子适应值， 将每个粒子适应值与个体极值进行比较， 如果较优， 则更新当前 的个体极值 ； 将每个粒子适应值与全局极值进行比较， 如果较优， 则更新当前的全局极值 ； c4） 更新每个粒子的位置和飞行速度 ； c5） 计算变异时机， 确定是否需要变异 ； c6） 重复步骤 c3 到 c5， 通过不断调整解向量， 将等值模型输出的功率曲线 P 与 P0进行 拟合， 直到达到迭代次数或收敛精度， 确定最优解向量， 从而获得等值风电机组各项参数的 最优值。 。

9、5. 如权利要求 1 所述的基于运行数据的风电场参数辨识及动态等值方法， 其特征在 于， 步骤 d 中验证等值模型的有效性方法为 ： 结合搜索辨识后的参数完善风电场等值模型， 对比采用详细风电场模型建立的仿真系统在相同输入条件下的动态响应特性， 以验证基于 运行数据双馈型风电场动态等值模型的有效性。 权 利 要 求 书 CN 103887815 A 3 1/10 页 4 基于运行数据的风电场参数辨识及动态等值方法 技术领域 0001 本发明属于智能电网中电力系统建模领域， 具体涉及一种基于运行数据的风电场 参数辨识及动态等值方法。 背景技术 0002 风力发电作为目前可再生能源开发利用中技术最。

10、成熟、 最具规模开发和商业化发 展前景的发电方式， 由于其在减轻环境污染、 调整能源结构、 解决偏远地区居民用电问题等 方面的突出作用， 越来越受到世界各国的重视并得到了广泛的开发和利用。 近几年， 风电已 成为增长速度最快的发电方式， 世界风力发电的发展速度每年都在 30% 以上。大型风电机 组的制造水平及其控制技术不断提高， 商业化风力机的容量已由原来的十几千瓦增长到了 目前的310MW， 接入电网后大型风电场的运行、 控制技术也日趋成熟。 风电并网已明显呈 现大规模、 区域化集中接入电网与更高电压等级接入电网的态势， 多个风电基地即将实现 百万甚至千万千瓦风电集中送出的目标。 0003 。

11、与常规能源不同， 由于风力发电出力的强波动性、 出力动态特性的功率弱支撑性 / 同步强随动性 / 暂态弱生存性等特点， 大规模风电并网将对电力系统的安全稳定运行带来 重要影响。因此， 在风电场建设前， 需要对风电场接入电网的可行性进行论证 ； 在运行调度 中， 也需要分析风电场接入的动态响应及其对电网的影响。 0004 为了定量分析风电场并网对局部地区电网的影响， 这就首先要求建立风电机组和 风电场精确的动态和稳态模型数学模型， 以便能够对系统的动态行为进行模拟仿真， 进而 为风电场的规划和运行调度提供依据。 0005 现有的风电并网问题研究中， 大家过多的把注意力集中于风电接入后对系统的各 。

12、种影响及防范措施上， 并未深入细致的思考所用风电模型是否可靠， 而模型的准确性将直 接影响分析结果的可靠性。 0006 目前风电机组模型参数大都由厂家给定或采用经验值， 然而在风电机组的实际运 行中， 其实际参数值并不是保持不变的， 而是随着工作环境和状况的变化在一定范围内变 动， 如频率不同、 温度变化引起的集肤效应， 会使电机的定转子电阻值发生变化， 磁场饱和 程度不同会影响到电机的绕组电感参数等。 模型参数给定值与实际值的偏离将直接影响仿 真结果的可信度， 因此有必要基于实测运行数据进行风电机组的参数辨识和模型验证。我 国 2009 年国家电网公司颁布的 风电并网运行控制技术规定 中对风。

13、电并网分析工作中应 采用风电机组的详细数学模型， 模型的参数应由风电场提供实测参数。对没有实测参数的 风电机组， 暂时可以采用同类机型的典型模型和参数， 风电机组模型和参数确定后应重新 校核。 0007 国内使用较多的大型仿真计算软件有 BPA（Bonneville PowerAdministration） ， BPA 中风电机组模型以 GE 机组为原型， 采用封装模型， 用户仅能设置几个外部参数， 而用户 使用时一般直接采用中国电科院提供的 GE 模型典型参数， 或稍作修改， 并未检验其模型的 有效性。国际电工委员会 （IEC） 风电机组电气仿真模型工作组 （TC88-WG27） 主要致力于。

14、建 说 明 书 CN 103887815 A 4 2/10 页 5 立标准的通用风电机组 / 风电场动态仿真模型并制定相关的风电机组模型方面的标准， 用 于电力系统分析与电网稳定计算研究， 但在双馈风电机组模型结构、 参数估计、 模型验证及 标准化方面仍处于探索阶段。 0008 在风电场动态建模方面， 大型 / 巨型风电场由几百台甚至上千台风电机组组成， 若对风电场中每台风电机都进行详细建模， 整个风电场的模型将变得十分复杂， 运算量巨 大， 可能导致电力系统分析与计算的维数灾， 且潮流计算极可能不收敛， 甚至无法进行仿 真。风电场动态等值则涉及更复杂的问题， 如地面粗糙度、 地势海拔、 尾流。

15、效应、 塔筒效应、 机组运行特性分散等等， 国内外尚未给出较有权威性的风电机组和风电场等值的通用模 型， 能够结合风电场实际运行数据进行建模和验证的模型更少。 0009 经对现有技术文献的检索发现， 双馈感应风力发电机的参数辨识分析 （金宇清， 赵 泽， 鞠平等.双馈感应风力发电机的参数辨识分析.高电压技术,2011,37(7):1700-1705)， 基于轨迹灵敏度和频域灵敏度研究了双馈感应风力发电机的模型参数的可辨识性， 选取发 电机参数 Rs、 Xs、 Xm、 Rr作为重点待辨识参数， 设计了基于蚁群优化算法的双馈感应风力发电 机参数辨识方法， 并验证了该辨识方法的有效性 ； 中国专利申。

16、请号为 ： 201210060923.7， 名 称为 ： 一种双馈电机参数在线辨识方法， 该申请案运用递推最小二乘法反复采样和迭代计 算， 得到双馈电机参数。然而， 对于风机模型来说， 除了风机本身固有的机械参数和发电机 本体参数之外， 风电机组的控制器参数也是很重要的， 它将直接影响风机模型输出特性。 据 检索， 目前对风力发电机组控制系统参数辨识的研究还没有相关报道。 0010 另经检索发现， 中国专利申请号为 ： 20110162399.X， 名称为 ： 一种含双馈风电机 组的大型风电场动态等值方法， 该申请案根据风电场输入风速的特性对风机进行分群， 而 同群的机组将被等值为一台机， 之。

17、后计算输入等值机的等值风速， 并得到分群机组的等值 参数 ； 然而， 申请案并没有对风力发电机及控制系统的参数进行辨识和模型验证， 分群机组 等值参数的有效性有待商榷 ； 中国专利申请号为 ： 200910238099.8， 名称为 ： 一种大规模风 电集中接入电网的风电场等值建模方法， 该申请案基于接入电网的各风电场并网点之间的 电气距离、 场内风电机组拓扑分布、 机型和控制模式， 采用同调机群识别聚类方法来确定同 调风电机群， 采用加权等效参数聚合方法求出等值风电机组各项动态电气参数， 根据静态 等值理论近似等效为等值机串接等值阻抗形式， 并求出相关的静态参数， 得到风电场等值 模型。 然。

18、而， 申请案采用加权等效参数聚合方法求出等值风电机组各项动态电气参数， 单台 发电机参数来源于厂家给定或采用经验值， 等值参数的有效性同样有待商榷。 0011 据检索， 目前基于运行数据对含双馈风电机组的风电场进行发电机和控制系统的 参数辨识， 并在此基础上采用优化算法进行动态等值建模研究， 还没有相关报道， 均属于空 白。 发明内容 0012 本发明的目的是为了解决上述现有技术中存在的不足之处， 提供一种基于运行数 据的风电场参数辨识及动态等值方法， 首先基于实测运行数据对风力发电机组控制模型进 行参数辨识， 并对控制器参数辨识结果进行模型验证， 然后使用实测数据对风电场进行聚 类分析， 基。

19、于全局最优位置变异粒子群算法获得含双馈风电机组的风电场动态等值模型。 建立的动态等值模型能准确地反映双馈感应发电机型风电场的动态特性， 可用于含双馈风 说 明 书 CN 103887815 A 5 3/10 页 6 电机组风电场接入电力系统稳定性分析中。 0013 为实现上述目的， 本发明通过以下技术方案实现 ： 0014 基于运行数据的风电场参数辨识及动态等值方法， 具体包括以下步骤 ： 0015 （1） 首先基于实测运行数据对风力发电机组控制模型进行参数辨识， 并对控制器 参数辨识结果进行模型验证 ； 0016 所述基于实测运行数据对风力发电机组控制模型进行参数辨识方法包括以下步 骤 ： 。

20、0017 a） 选择需要重点辨识的控制器模型参数。需要进行辨识的必须是变化时能引起 风电机输出功率明显变化的控制器参数， 经测试， 模型中功率控制器的控制参数 Kp、 Tp、 Kq、 Tq 和电流控制器的控制参数 Kd、 Td、 Kq1、 Tq1 对发电机输出功率的影响较大 ； 由此， 选用上 述 8 个参数进行辨识 ； 其中， Kp、 Tp 分别是有功功率控制器的比例常数和积分常数 ； Kq、 Tq 分别是无功功率控制器的比例常数和积分常数 ； Kd、 Td 分别是有功电流控制器的比例常数 和积分常数 ； Kq1、 Tq1 无功电流控制器的比例常数和积分常数 ； 0018 b） 搭建参数辨识。

21、模块。所述参数辨识模块包括测量文件， 要辨识的模型元件 1 n、 参数辨识比较器、 优化器， 目标函数默认是实测数据和仿真结果的差值的平方和 ： 0019 0020 其中， Mi是测量文件中的数据 ； Si是仿真文件中的数据 ； i是权重因子， 可根据实 际情况进行调整 ； p 是指数值， 默认为 2 ； 目标函数越小指示参数辨识效果越好 ； 0021 测量文件通过读取实际的发电机有功功率给比较器模块， 再通过功率测量环节将 仿真中发电机有功功率特性提供给比较器模块， 根据实际的数据和仿真结果计算目标函数 值 ； 0022 c） 设置参数初始值。所述参数初始值为双馈风力发电机组经典控制模型参数。

22、 ； 0023 d） 实测数据处理及导入。所述实测数据为风电机组的输入风速 ； 0024 e） 获得控制器参数辨识结果 ； 0025 f） 基于控制器参数辨识结果， 进行模型验证。将获得的控制器参数辨识结果作为 该时段风电机组控制器模型参数， 与对应采样功率数据进行比较， 以验证控制器参数辨识 的可用性。 （2） 结合运行数据， 选取能反映扰动下的风电机组和风电场的特征变量， 应用改 进的模糊 K 均值动态聚类算法进行机群划分 ： 结合风电机组和风电场实测数据， 分析风力 发电机组的有功、 无功、 电压、 转速、 桨距角等参量对风力发电机组和风电场暂态特性的影 响， 并提炼出能反映大扰动下的风。

23、电机组和风电场的特征变量。以大扰动发生前一刻、 大 扰动期间和大扰动消失后风力发电机组具有相同或相近的运行点为风电机组动态机群划 分为准则， 综合利用各风电机组在大扰动期间实测风速、 有功、 无功等反映发电机组暂态运 行点的特征变量的测量数据为数据样本， 应用改进的模糊 K 均值动态聚类算法进行机群划 分， 使同一机群内的风电机组的动态响应尽可能相似， 不同机群的风电机组的动态响应差 别尽可能大。 借鉴同调等值的思想， 认为同一风电组群具有相同或相近的动态过程， 并将其 等效为一台等值风电机组。 0026 本发明的优选方案 ： 采用扰动期间实测风力发电机组有功功率作为聚类分群的特 征变量， 具。

24、体包括以下步骤 ： 说 明 书 CN 103887815 A 6 4/10 页 7 0027 a） 首先选取风电场在扰动期间设定时间段内各风力发电机组的实测有功功率数 据作为运行数据样本， 样本的测量数据点数为 L， 将 N 台风力发电机组的有功功率数据记录 按行建立实测有功功率样本矩阵 ： 0028 0029 其中第 i 行表示第 i 台风力发电机组实际测得的有功功率记录， 第 j 列表示 n 台 风力发电机组在第 j 个时刻测得的有功功率。把实测有功功率样本矩阵 P 的每一行看作为 空间中的一个数据点 ( 样本数量为 n， 样本维数为 k) ； 0030 b） 确定聚类数目 K， 将全部数。

25、据随机分成 K 类。其中， 聚类的簇的个数 K 值根据实 际风电场规模和模型精度要求选取 ； 0031 c） 定义 j（xi） 是第 i 个样本对于第 j 类的隶属度函数， 用隶属度函数定义相应 的聚类损失函数为 ： 公式一 0032 其中， b1 是一个可以控制聚类结果的模糊程度的常数 ； xi是聚类样本。 0033 改进模糊 K 均值算法要求所有样本对于各个聚类的隶属度总和为 n， 即 0034 公式二 0035 在公式一条件下， 令 Jf对 mi和 j的偏导数为 0， 可得必要条件 ： 0036 公式三 0037 公式四 0038 d） 重复下面的运算， 直到各个样本的隶属度值稳定 ： 。

26、0039 用当前的聚类中心根据公式四计算隶属度函数 ； 0040 用当前的隶属度函数按公式三更新计算各类聚类中心 ； 0041 当算法收敛时， 就得到了各类的聚类中心和各个样本对与各类的隶属度值， 从而 完成了基于改进模糊的 K 均值动态聚类算法的风力发电机组机群划分， 风电场的风力发电 机组按相似动态特性被划分成 K 群。 0042 （3） 基于全局最优位置变异粒子群算法对风电场进行网络简化、 参数聚合， 获得风 电场动态等值模型 : 基于全局最优位置变异粒子群算法以等值前后风电场有功、 无功、 电 压偏差的平方和作为优化目标， 在整个解空间上进行寻优， 按概率意义找到满足优化目标 函数的等。

27、值风电机组模型的参数。其中， 风电场等值模型的功率控制器和电流控制器参数 说 明 书 CN 103887815 A 7 5/10 页 8 采用双馈风机的单机模型经过辨识后的参数值。这样， 整个风电场就能用几台等值风电机 组来代替。据此基于运行数据的风电场内机群划分、 网络简化和参数聚合三个步骤实现了 风电场动态等值建模。 0043 所述基于全局最优位置变异粒子群算法对风电场进行动态等值方法包括以下步 骤 ： 0044 a） 搭建双馈风电场等值模型。 其中， 电机参数未知， 风电场等值模型的功率控制器 和电流控制器参数采用双馈风机的单机模型经过辨识后的参数值， 机组参数按风机容量的 变化作等比例。

28、增加， 风轮半径保持不变。 0045 定义目标函数 ：将大小扰动下实测有功功率数据作为功率参考， 将有功功率偏差平方和送入MATLAB中workplace区， 以供适应度函数采集调用。 式中， M为 等值前后有功功率的采样点数， Poi为样本的有功功率曲线中第 i 个采样点的有功功率值， Pi为等值后第 i 点的有功功率值。经过处理就将双馈电机参数辨识问题转化为参数优化问 题 ； 0046 b） 选择待辨识参数 Rs、 Rr、 Ls、 Lr、 Xm， 初始化粒子群体的位置和速度， 设定种群大小， 设置粒子位置和速度上下限， 设定迭代次数、 收敛精度和相关权系数。其中， Rs和 Ls分别为 定子。

29、电阻和定子漏抗 ； Rr和 Lr分别为归算后转子电阻和转子电抗 ； Xm为励磁电抗 ； 0047 c） 计算粒子适应值， 将每个粒子适应值与个体极值进行比较， 如果较优， 则更新 当前的个体极值 ； 将每个粒子适应值与全局极值进行比较， 如果较优， 则更新当前的全局极 值 ； 0048 d） 更新每个粒子的位置和飞行速度， 在第 t+1 次迭代计算时， 粒子 i 更新规则如 下 ： 0049 vik(t+1) wvik(t)+c1rand1(0,1)(pik(t)-xik(t)+c2rand2(0,1)(pgk(t) -xik(t) 0050 xik(t+1) xik(t)+vik(t+1) 。

30、0051 式中 为惯性权重 ,c1、 c2为两个学习因子 ； rand1(0,1) 和 rand2(0,1) 为两个均 匀分布在(0， 1)之间的随机数 ； i1,2,.,m;k1,2,.,d ； 另外， 粒子在每一维的速度 vi都被一个最大速度 Vmax所限制。vik(t+1) 和 xik(t+1) 分别为粒子 i 的第 t+1 次飞行速度 和位置 ； 0052 e） 计算变异时机， 确定是否需要变异 ； 0053 f） 重复步骤 c 到 e， 通过不断调整解向量， 将等值模型输出的功率曲线 P 与 P0进 行拟合， 直到达到迭代次数或收敛精度， 确定最优解向量， 从而获得等值风电机组各项参。

31、数 的最优值。 0054 （4） 模型验证， 在相同扰动输入下， 比较风电场动态等值模型与详细模型动态响 应， 验证风电场动态等值模型的有效性 ： 结合搜索辨识后的参数完善风电场等值模型， 对比 采用详细风电场模型建立的仿真系统输入输出特性， 以验证基于运行数据双馈型风电场动 态等值模型的有效性。 0055 本发明具有以下的优点和有益效果 ： 0056 1、 本发明基于实测运行数据对含双馈风电机组的风电场进行动态等值建模， 首先 说 明 书 CN 103887815 A 8 6/10 页 9 结合经过校正处理后的采样数据， 对风电机组的控制器参数进行直接辨识， 使风电机组的 功率输入输出特性与。

32、实际相符合， 得到一套具有一定实用性的双馈风电机组模型。然后在 聚类分群结果的基础上， 对风电场进行动态等值， 建立的风电场动态等值模型能够准确地 反映风电场并网点的动态特性， 具有重要的工程应用价值， 可用于含双馈风电机组风电场 接入电力系统稳定性分析中。 0057 2、 本发明通过对风电场内相关机群划分、 网络简化和参数聚合三个步骤， 大大降 低了模型的阶数， 显著减少仿真时间， 降低阶数， 从而避免维数灾， 从而解决大型 / 巨型风 电场的建模难题， 结合运行数据的风电场模型将为含风电场电力系统的规划与运行调度提 供理论支撑。 0058 3、 本发明采用的全局最优位置变异的粒子群算法设计。

33、了一种不同于传统变异的 新的变异方法以避免陷入 “早熟” ， 这样整个粒子群也不至于因全局最优位置处于某个局部 最优点时的强大吸引而快速聚集到局部最优点。这也确保了动态等值结果的有效性。 附图说明 0059 图 1 是本发明基于运行数据含双馈风电机组的风电场参数辨识及动态等值建模 方法总流程图 ； 0060 图 2 是本发明实施例国内某地双馈风电机组构建的风电场结构示意图 ； 0061 图 4 是本发明基于实测运行数据对风力发电机组控制模型进行参数辨识的流程 图 ； 0062 图5是本发明基于改进的模糊K均值动态聚类算法对风电场中风电机组进行聚类 分群的流程图 ； 0063 图 6 是本发明基。

34、于全局最优位置变异粒子群算法对风电场进行动态等值方法的 流程图 ； 0064 图 7 是本发明基于全局最优位置变异粒子群算法寻优结束前中风电场的并网点 电压、 电流、 有功功率、 无功功率、 直流侧电压和转子转速波形图 ； 0065 图8a、 图8b、 图8c分别是本发明实施例风电场详细模型和动态等值模型在输电线 路三相短路故障过程中并网点电压、 有功功率和无功功率的暂态响应过程波形图。 具体实施方式 0066 下面结合实施例及附图， 对本发明做进一步地详细说明， 但本发明的实施方式不 限于此。 0067 图 1 是本发明基于运行数据含双馈风电机组的风电场参数辨识及动态等值建模 方法总流程图。。

35、 方法包括以下步骤 ： 步骤a ： 基于实测数据对风电机组控制模型的功率和电 流控制器参数进行参数辨识， 并对辨识结果进行通用性验证 ； 步骤 b ： 结合运行数据， 选取 能反映扰动下的风电机组和风电场的特征变量， 应用改进的模糊 K 均值动态聚类算法进行 机群划分 ； 步骤c： 基于全局最优位置变异粒子群算法对风电场进行网络简化、 参数聚合， 获 得风电场动态等值模型 ； 步骤 d ： 在相同扰动输入下， 比较风电场动态等值模型与详细模型 动态响应， 验证等值模型的有效性。 0068 图 2 是本发明实施例国内某地 33 台双馈风电机组构建的风电场结构示意图。所 说 明 书 CN 1038。

36、87815 A 9 7/10 页 10 述风电场每台双馈风机均经由一台变比为 0.69kV/35kV 的理想变压器连接到每组集总母 线上， 每 11 台双馈风机为一组， 再经由风电场集总变压器 Trf_Grid（变比为 35kV/110kV） 连接到外部系统， External Grid 是外部等值网络。 0069 图 3 是本发明基于实测运行数据对风力发电机组控制模型进行参数辨识的系统 框图。所述参数辨识模块包括测量文件、 要辨识的模型元件 1 n、 参数辨识比较器、 优化 器。 0070 图 4 是本发明基于实测运行数据对风力发电机组控制模型进行参数辨识的流程 图。所述步骤 a， 具体包括。

37、以下步骤 ： 0071 a） 选择需要重点辨识的控制器模型参数。 需要进行辨识的必须是变化时能引起风 电机输出功率明显变化的控制器参数， 经测试， 模型中功率控制器的控制参数 Kp、 Tp、 Kq、 Tq 和电流控制器的控制参数 Kd、 Td、 Kq、 Tq 对发电机输出功率的影响较大。由此， 选用上述 8 个参数进行辨识 ； 其中， Kp、 Tp 分别是有功功率控制器的比例常数和积分常数 ； Kq、 Tq 分别 是无功功率控制器的比例常数和积分常数 ； Kd、 Td 分别是有功电流控制器的比例常数和积 分常数 ； Kq1、 Tq1 无功电流控制器的比例常数和积分常数 ； 0072 b） 搭建。

38、参数辨识模块。目标函数默认是实测数据和仿真结果的差值的平方和 ： 0073 0074 其中， Mi是测量文件中的数据 ； Si是仿真文件中的数据 ； i是权重因子， 可根据实 际情况进行调整 ； p 是指数值， 默认为 2。目标函数越小指示参数辨识效果越好。 0075 测量文件通过读取实际的发电机有功功率给比较器模块， 再通过功率测量环节将 仿真中发电机有功功率特性提供给比较器模块， 根据实际的数据和仿真结果计算目标函数 值 ； 0076 c） 设置参数初始值。所述参数初始值为双馈风力发电机组经典控制模型参数 ； 0077 d） 实测数据处理及导入。所述实测数据为风电机组的输入风速。 0078。

39、 作为优选实施方式， 对实际的采样数据经过数据过滤与还原处理后作为输入风速 数据 ； 0079 作为优选实施方式， 为了能更好地反映双馈感应电机的动态性能， 输入数据选择 风速波动较大的数据进行辨识， 设定标准方差在 0.75 以上的即为风速波动较大的数据 ； 0080 e） 获得控制器参数辨识结果 ； 0081 f） 模型验证。将获得的控制器参数辨识结果作为该时段风电机组控制器模型参 数， 与对应采样功率数据进行比较， 以验证控制器参数辨识的可用性。 0082 图5是本发明基于改进的模糊K均值动态聚类算法对风电场中风电机组进行聚类 分群的流程图。 结合风电机组和风电场实测数据， 分析风力发电。

40、机组的有功、 无功、 电压、 转 速、 桨距角等参量对风力发电机组和风电场暂态特性的影响， 并提炼出能反映大扰动下的 风电机组和风电场的特征变量。以大扰动发生前一刻、 大扰动期间和大扰动消失后风力发 电机组具有相同或相近的运行点为风电机组动态机群划分为准则， 综合利用各风电机组在 大扰动期间实测风速、 有功、 无功等反映发电机组暂态运行点的特征变量的测量数据为数 据样本， 应用改进的模糊 K 均值动态聚类算法进行机群划分， 使同一机群内的风电机组的 动态响应尽可能相似， 不同机群的风电机组的动态响应差别尽可能大。借鉴同调等值的思 说 明 书 CN 103887815 A 10 8/10 页 1。

41、1 想， 认为同一风电组群具有相同或相近的动态过程， 并将其等效为一台等值风电机组。 0083 本发明的优选方案 ： 采用扰动期间实测风力发电机组有功功率作为聚类分群的特 征变量， 所述步骤 b， 具体包括以下步骤 ： 0084 a） 首先选取风电场在扰动期间设定时间段内各风力发电机组的实测有功功率数 据作为运行数据样本， 样本的测量数据点数为 L， 将 N 台风力发电机组的有功功率数据记录 按行建立实测有功功率样本矩阵 ： 0085 0086 其中第 i 行表示第 i 台风力发电机组实际测得的有功功率记录， 第 j 列表示 n 台 风力发电机组在第 j 个时刻测得的有功功率。把实测有功功率样。

42、本矩阵 P 的每一行看作为 空间中的一个数据点 ( 样本数量为 n， 样本维数为 k) ； 0087 b） 确定聚类数目 K=3， 将全部数据随机分成 K 类。其中， 聚类的簇的个数 K 值根据 实际风电场规模和模型精度要求选取。计算每类的重心， 并将重心作为其代表点。各代表 点自成一类， 将样本依顺序归入与其距离最近的代表点一类， 并重新计算该类的重心以代 替原来的代表点， 然后再计算下一个样本的归类， 直至所有的样本都归到相应的类中位置， 初始化各个聚类中心 mi完成 ； 0088 c） 定义 j（xi） 是第 i 个样本对于第 j 类的隶属度函数， 用隶属度函数定义相应 的聚类损失函数为。

43、 ： 公式一 0089 其中， b1 是一个可以控制聚类结果的模糊程度的常数 ； xi是聚类样本。 0090 改进模糊 K 均值算法要求所有样本对于各个聚类的隶属度总和为 n， 即 0091 公式二在公式一条件下， 令 Jf对 mi和 j的偏导数为 0， 可 得必要条件 ： 0092 公式三 0093 公式四 0094 d） 重复下面的运算， 直到各个样本的隶属度值稳定 ： 0095 用当前的聚类中心根据公式四计算隶属度函数 ； 0096 用当前的隶属度函数按公式三更新计算各类聚类中心。 0097 当算法收敛时， 就得到了各类的聚类中心和各个样本对与各类的隶属度值， 从而 说 明 书 CN 1。

44、03887815 A 11 9/10 页 12 完成了基于改进模糊的 K 均值动态聚类算法的风力发电机组机群划分， 风电场的风力发电 机组按相似动态特性被划分成 K=3 群。 0098 图 6 是本发明基于全局最优位置变异粒子群算法对风电场进行动态等值方法的 流程图。 基于全局最优位置变异粒子群算法以等值前后风电场有功、 无功、 电压偏差的平方 和作为优化目标， 在整个解空间上进行寻优， 按概率意义找到满足优化目标函数的等值风 电机组模型的参数。其中， 风电场等值模型的功率控制器和电流控制器参数采用双馈风机 的单机模型经过辨识后的参数值。这样， 整个风电场就能用几台等值风电机组来代替。据 此基。

45、于运行数据的风电场内机群划分、 网络简化和参数聚合三个步骤实现了风电场动态等 值建模。 0099 所述步骤 c， 具体包括以下步骤 ： 0100 a） 搭建双馈风电场等值模型。 其中， 电机参数未知， 风电场等值模型的功率控制器 和电流控制器参数采用双馈风机的单机模型经过辨识后的参数值， 机组参数按风机容量的 变化作等比例增加， 风轮半径保持不变。 0101 定义目标函数 ：将大小扰动下实测有功功率数据作为功率参考， 将有功功率偏差平方和送入MATLAB中workplace区， 以供适应度函数采集调用。 式中， M为 等值前后有功功率的采样点数， Poi为样本的有功功率曲线中第 i 个采样点的。

46、有功功率值， Pi为等值后第 i 点的有功功率值。经过处理就将双馈电机参数辨识问题转化为参数优化问 题 ； 0102 b） 选择待辨识参数 Rs、 Rr、 Ls、 Lr、 Lm， 初始化粒子群体的位置和速度， 设定种群大小， 设置粒子位置和速度上下限， 设定迭代次数、 收敛精度和相关权系数。其中， Rs和 Ls分别为 定子电阻和定子漏抗 ； Rr和 Lr分别为归算后转子电阻和转子电抗 ； Xm为励磁电抗 ； 0103 作为优选实施方式， 采用容量加权法求取等值风电机组的发电机参数， 然后分别 乘以 2 和 0.3 作为搜索范围的上限与下限。 0104 容量加权法中， 等值机的参数为 ： 010。

47、5 0106 其中， Zs=Rs+jXs表示定子阻抗 ； Zr=Rr+jXr表示转子漏抗 ； Zm=jXm表示励磁电抗。AS 表示参与等值的全部双馈发电机的集合， 下标 A 表示等值后的双馈发电机， wi为容量权系 数 ； i 表示第 i 台发电机 ； 0107 c） 计算粒子适应值， 将每个粒子适应值与个体极值进行比较， 如果较优， 则更新 当前的个体极值 ； 将每个粒子适应值与全局极值进行比较， 如果较优， 则更新当前的全局极 值 ； 0108 d） 更新每个粒子的位置和飞行速度， 在第 t+1 次迭代计算时， 粒子 i 更新规则如 下 ： 0109 vik(t+1) wvik(t)+c1。

48、rand1(0,1)(pik(t)-xik(t)+c2rand2(0,1)(pgk(t) -xik(t) 0110 xik(t+1) xik(t)+vik(t+1) 说 明 书 CN 103887815 A 12 10/10 页 13 0111 式中 为惯性权重 ,c1、 c2为两个学习因子 ； rand1(0,1) 和 rand2(0,1) 为两个均 匀分布在(0， 1)之间的随机数 ； i1,2,.,m;k1,2,.,d ； 另外， 粒子在每一维的速度 vi都被一个最大速度 Vmax所限制。vik(t+1) 和 xik(t+1) 分别为粒子 i 的第 t+1 次飞行速度 和位置 ； 011。

49、2 e） 计算变异时机， 确定是否需要变异 ； 0113 f） 重复步骤 c 到 e， 通过不断调整解向量， 将等值模型输出的功率曲线 P 与 P0进 行拟合， 直到达到迭代次数或收敛精度， 确定最优解向量， 从而获得等值风电机组各项参数 的最优值。 0114 图 7 是本发明基于全局最优位置变异粒子群算法寻优结束前风电场的并网点电 压Vabc_B、 电流Iabc_B、 有功功率P、 无功功率Q、 直流侧电压Vdc和转子转速wr波形图。 由 图可知算法结束时， 风电场各参量已达到稳态。 0115 图 8 是本发明实施例风电场详细模型和动态等值模型在输电线路三相短路故障 过程中并网点电压、 有功功率和无功功率的暂。