一种直驱风机变流器的控制参数辨识方法及系统与流程

1.本发明涉及风机变流器技术领域，尤其涉及一种直驱风机变流器的控制参数辨识方法及系统。


背景技术：

2.本部分的陈述仅仅是提供了与本发明相关的背景技术信息，不必然构成在先技术。
3.参数辨识作为控制理论的一个重要分支，始于20世纪60年代，并迅速发展起来。目前，参数辨识理论已经相对成熟，研究成果已广泛应用于电力系统中，如控制参数、机械参数和系统等效模型参数的辨识等。
4.参数辨识理论可分为经典辨识方法和现代辨识方法两大类。对于经典的辨识方法，主要包括频域法和时域法。频域法主要包括最大熵法、静止频率响应法和在线频率响应法。时域法是指最小二乘法、卡尔曼滤波法、分段线性多项式函数法等。对于现代辨识方法，主要包括模拟进化算法、模拟退火算法、禁忌搜索算法、人工神经网络算法等。
5.对于经典辨识方法的应用，最小二乘(ls)技术以其表达式简单、计算量小、对白噪声抑制鲁棒性强等优点而成为主流。ls方法最早由高斯提出，用于计算和预测行星轨道。对于ls方法在电力系统中的应用，文章《基于最小二乘法的光伏逆变器模型辨识》中将ls方法应用于光伏并网逆变器的控制参数辨识，其他研究也有将ls方法应用于异步电机转子和变压器绕组的电气参数辨识。
6.目前，风电机组变流器控制参数辨识方面主要使用的算法是卡尔曼滤波及遗传算法。然而，两种算法均对参数的初值要求较高，本发明将ls方法用于识别dpmsg的变流器控制参数中，该算法对初值无要求且有较高的辨识精度。


技术实现要素：

7.本发明为了解决上述问题，提出了一种直驱风机变流器的控制参数辨识方法及系统，本发明能够有效提高dpmsg的控制参数辨识的适用性。
8.根据一些实施例，本发明采用如下技术方案：
9.一种直驱风机变流器的控制参数辨识方法，包括：
10.获取直驱风机的输入变量；
11.根据获取的输入变量，利用dpmsg模型得到动态输出特性；
12.根据得到的动态输出特性，利用lsa算法得到直驱风机变流器的控制参数的辨识结果。
13.进一步地，所述直驱风机变流器的输入变量包括：正常运行状态下，有功功率参考值和测量值、定子电流的q轴分量、直流母线电压的参考值和测量值、输出电流的dq轴分量、无功功率参考值和测量值；低电压穿越状态下，输出电流q轴分量的参考值和测量值。
14.进一步地，所述利用dpmsg模型得到动态输出特性包括：正常运行状态下，定子电
压的q轴分量的测量值，输出电压的dq轴分量的测量值；低电压穿越状态下，输出电压的q轴分量的测量值。
15.进一步地，所述dpmsg模型包括永磁同步发电机模型、机侧及网侧变流器模型、变流器控制模型、直流母线电容模型、升压变压器模型及外部无穷大电网模型。
16.进一步地，变流器控制结构主要为双环pi控制器，所述双环pi控制器包括内环控制和外环控制。
17.进一步地，所述lsa算法获得辨识结果的必要条件是所得观测矩阵伪逆的存在。
18.进一步地，所述直驱风机变流器包括机侧变流器和网侧变流器。
19.进一步地，所述机侧变流器根据转子磁场矢量控制调节永磁同步发电机的有功功率输出和转子转速。
20.进一步地，所述网侧变流器根据电网电压矢量控制维持直流母线电压的稳定，并调节风电机组的功率因数。
21.一种直驱风机变流器的控制参数辨识系统，包括：
22.数据获取模块，被配置为，获取直驱风机的输入变量；
23.动态模块，被配置为，根据获取的输入变量，利用dpmsg模型得到动态输出特性；
24.辨识模块，被配置为，根据得到的动态输出特性，利用lsa算法得到直驱风机变流器的控制参数的辨识结果。
25.一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行所述的一种直驱风机变流器的控制参数辨识方法。
26.一种终端设备，包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行所述的一种直驱风机变流器的控制参数辨识方法。
27.与现有技术相比，本发明的有益效果为：本发明所获取的输入变量及动态输出特性均为实际风电机组中的可测量值；所用ls算法理论完善，计算简单，具有较好的抗噪音能力，计算速度快，计算精度较高；所用ls算法对初值无要求，具有很高的工程适用性；
28.本发明通过将ls方法用于dpmsg控制参数的辨识中，有效提高了识别dpmsg的控制参数的适用性，并具有较高的辨识精度。
附图说明
29.构成本技术的一部分的说明书附图用来提供对本技术的进一步理解，本技术的示意性实施例及其说明用于解释本技术，并不构成对本技术的不当限定。
30.图1是本实施例1的双环控制的典型结构图；
31.图2是本实施例1的测试流程图；
32.图3是本实施例2的dpmsg风力发电系统图；
33.图4中，(a)为msc的p控制框图；(b)为gsc电压控制框图；(c)为gsc在正常运行模式下的q控制框图；(d)为故障模式下gsc的q控制框图；
34.图5中，(a)为风速阶跃变化图；(b)为dpmsg的有功功率输出图；(c)为dpmsg的无功功率输出图；(d)为dpmsg直流母线电压图；
35.图6是本实施例1的故障状态下35kv侧电网电压图。
具体实施方式：
36.下面结合附图与实施例对本发明作进一步说明。
37.应该指出，以下详细说明都是例示性的，旨在对本技术提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本技术所属技术领域的普通技术人员通常理解的相同含义。
38.需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本技术的示例性实施方式。此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作和/或它们的组合。
39.名词解释：dpmsg，直驱永磁同步风力发电机，简称直驱风机。
40.实施例1
41.如图1所示，本发明实施例1提供了一种直驱风机变流器的控制参数辨识方法，包括：获取直驱风机的输入变量；
42.根据获取的输入变量，利用dpmsg模型得到动态输出特性；
43.根据得到的动态输出特性，利用lsa算法得到直驱风机变流器的控制参数的辨识结果。
44.进一步地，所述直驱风机变流器的输入变量包括：正常运行状态下，有功功率参考值和测量值、定子电流的q轴分量、直流母线电压的参考值和测量值、输出电流的dq轴分量、无功功率参考值和测量值；低电压穿越状态下，输出电流q轴分量的参考值和测量值。
45.进一步地，所述利用dpmsg模型得到动态输出特性包括：正常运行状态下，定子电压q轴分量的测量值、输出电压的dq轴分量的测量值；低电压穿越状态下，输出电压的q轴分量的测量值。
46.进一步地，所述dpmsg模型包括永磁同步发电机模型、机侧及网侧变流器模型、变流器控制模型、直流母线电容模型、升压变压器模型及外部无穷大电网模型。
47.进一步地，变流器控制结构主要为双环pi控制器，所述双环pi控制器包括内环控制和外环控制。
48.进一步地，所述lsa算法获得辨识结果的必要条件是所得观测矩阵伪逆的存在。
49.进一步地，所述直驱风机变流器包括机侧变流器和网侧变流器。
50.进一步地，所述机侧变流器根据转子磁场矢量控制调节永磁同步发电机的有功功率输出和转子转速。
51.进一步地，所述网侧变流器根据电网电压矢量控制维持直流母线电压的稳定，并调节风电机组的功率因数。
52.具体的，
53.s1、直驱风机的输入变量，包括：
54.p
_ref
ꢀꢀꢀꢀꢀꢀꢀ
有功功率参考值；
55.u
dc_ref
ꢀꢀꢀꢀꢀ
直流母线电压的参考值；
56.q
_ref
ꢀꢀꢀꢀꢀꢀꢀ
无功功率参考值；
57.i
sq_ref
ꢀꢀꢀꢀꢀ
定子电流q轴分量的参考值；
58.i
d_ref
ꢀꢀꢀꢀꢀꢀ
输出电流d轴分量的参考值；
59.i
q_ref
ꢀꢀꢀꢀꢀꢀ
输出电流q轴分量的参考值；
60.i
q_ref_fault 故障状态下输出电流q轴分量的参考值；
61.p
_m
ꢀꢀꢀꢀꢀꢀꢀꢀ
有功功率测量值；
62.u
dc_m
ꢀꢀꢀꢀꢀꢀꢀ
直流母线电压测量值；
63.q
_m
ꢀꢀꢀꢀꢀꢀꢀꢀ
无功功率测量值；
64.i
sq
ꢀꢀꢀꢀꢀ
定子电流q轴分量的测量值；
65.i
d_m
ꢀꢀꢀꢀꢀꢀꢀ
输出电流d轴分量的测量值；
66.i
q_m
ꢀꢀꢀꢀꢀꢀꢀ
输出电流q轴分量的测量值；
67.u
sq
ꢀꢀꢀꢀꢀꢀꢀꢀ
定子电压的q轴分量
68.udꢀꢀꢀꢀꢀꢀꢀꢀꢀ
d轴电压分量；
69.uqꢀꢀꢀꢀꢀꢀꢀꢀꢀ
q轴电压分量；
70.k
px
ꢀꢀꢀꢀꢀꢀꢀꢀ
pi调节器的比例系数；
71.k
ix
ꢀꢀꢀꢀꢀꢀ
pi调节器的积分系数。
72.s2、利用最小二乘法得到直驱风机变流器控制参数的辨识结果，其中，dpmsg的控制系统主要由双环pi控制器组成，双环控制器包括内环控制和外环控制，如图1所示。
73.在复频域中，双环pi控制的输入和输出之间的关系表示如下：
[0074][0075]
式中，x1(s)＝x
1ref
(s)-x
1m
(s)，k
p1
、k
p2
是比例系数，k
i1
、k
i2
是积分系数。
[0076]
为了辨识双环pi控制器的参数，表达式(1)可以写为离散形式，如下所示：
[0077][0078]
式中，a＝k
p1kp2
，b＝k
p1ki2
k
p2ki1
，c＝k
i1ki2
，d＝-k
p2
，e＝-k
i2
，δt是采样间隔。
[0079]
因此，只要辨识出方程系数a、b、c、d和e，就可以识别双环pi控制器的参数。
[0080]
对于表达式(2)中所示的双环pi控制器的离散表达形式，其可写为：
[0081][0082]
其中，y(k)为输出向量，φ(k)为观测矩阵，θ为参数向量，e(k)为误差向量。
[0083]
如果进行n(n》》2)次观测，k＝1,2,3
…
n，则可获得以下方程式：
[0084][0085]
则有：y＝φθ e
ꢀꢀ
(5)
[0086]
其中，y为n
×
1输出向量，φ为n
×
5观测矩阵，θ为5
×
1参数向量，e为n
×
1误差向量。
[0087]
评价函数：
[0088][0089]
根据最小二乘算法，有：
[0090][0091]
因此：
[0092][0093]
使用最小二乘法辨识关键参数的必要条件是观测矩阵φ的伪逆的存在。
[0094]
s3、dpmsg模型的构建
[0095]
通过对直驱风电机组和无穷大母线系统的动态仿真，对dpmsg变流器控制参数辨识的ls方法进行了测试。测试系统建立在digsilentpowerfactory平台上。dpmsg系统的结构图如图3所示。表1给出了dpmsg的电气参数。
[0096]
表1 dpmsg参数
[0097][0098]
图4给出dpmsg的机侧变流器(msc)和网侧变流器(gsc)的控制框图。msc根据最大功率跟踪(mppt)控制调节有功功率输出和转子转速。在正常运行情况下，gsc控制确保直流母线电压的稳定，并调节馈入电网的无功功率为零(功率因数为1)。在故障情况下，gsc控制为电网提供额外的无功功率，以支持电压稳定。表2给出了变流器的控制参数。
[0099]
表2 dpmsg的控制参数
[0100][0101]
为了激发风机的动态输出响应，给风机施加一组具有确定性风速的阶跃响应，测试风力涡轮机模型的输出响应。启动风速为3m/s，额定风速为10m/s(5mw有功功率输出)。如
图5中的(a)所示，模拟风速呈现阶跃变化，即每一步增加0.5m/s，每一阶段风速保持60秒。图5中的(b)、(c)、(d)显示了有功功率和无功功率输出，以及直流母线的电压。
[0102]
s4、为了测试lsa在双环pi控制器参数辨识中的应用，提出了一种计算机模拟试验流程，试验流程框图如图2所示。第一步使用模拟软件“digsilentpowerfactory”对直驱风力发电系统进行建模。第二步表示通过给digsilent仿真软件中的直驱风力发电系统施加扰动获得的动态信号。第三步根据lsa算法估计的dpmsg的控制参数，第四步保存所得结果，第五步将所得结果与真实值进行比较，第六步分析误差向量
[0103]
s5、在digsilentpowerfactory模拟低电压穿越故障，在3s时模拟三相短路故障并保持0.7s。为了估计故障状态下gsc的q控制参数，将信号iq_ref_fault和iq_m用作观测数据，并将uq用作输出数据。采样频率为10000hz，数据窗大小为0.7s(从3s到3.7s)。辨识结果如表3所示：
[0104]
表3大扰动数据下q控制参数的辨识结果
[0105][0106]
对于msc的p控制参数的辨识，由于在执行低电压穿越策略时，为满足电网支撑电压的需要，首先要保证无功功率的输出，有功功率的输出会相应的受到限制。因此，对于msc的p控制参数、gsc的udc控制参数及q控制的外环控制参数的辨识，应在风电场正常运行工况下进行，通过模拟风电场输出特性，辨识其控制参数。为了模拟风电场的动态响应，对风电场等值模型施加一组风速的阶跃变化信号。启动风速为3m/s，额定风速为10m/s，如图5(a)所示，模拟风速呈现阶跃变化，即每步增加0.5m/s，每一阶段风速保持60秒。
[0107]
对于msc的p控制估计，信号p
_ref
、p
_m
和i
sq
用作观测数据，u
sq
用作输出数据。为了估计gsc的udc控制，信号u
dc_ref
、u
dc_m
和i
d_m
被用作观测数据，ud是输出数据。对于gsc的q控制估计，q
_ref
、q
_m
和i
q_m
作为观察数据，uq作为输出数据。数据采样频率为1000hz，滑动数据窗口大小为60s，因此总共使用了15个数据窗口。控制参数的辨识结果如表4所示。
[0108]
表4小扰动数据下控制参数的辨识结果
[0109][0110]
实施例2
[0111]
本实施例2提供了一种直驱风机变流器的控制参数辨识系统，包括：
[0112]
数据获取模块，被配置为，获取直驱风机的输入变量；
[0113]
动态模块，被配置为，根据获取的输入变量，利用dpmsg模型得到动态输出特性；
[0114]
辨识模块，被配置为，根据得到的动态输出特性，利用lsa算法得到直驱风机变流器的控制参数的辨识结果。
[0115]
实施例3
[0116]
本实施例3提供了一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行本实施例1提供的一种直驱风机变流器的控制参数辨识方法。
[0117]
实施例4
[0118]
本实施例4提供了一种终端设备，包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行本实施例1提供的一种直驱风机变流器的控制参数辨识方法。
[0119]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0120]
本技术是参照根据本技术实施例的方法和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0121]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0122]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0123]
以上所述仅为本技术的优选实施例而已，并不用于限制本技术，对于本领域的技术人员来说，本技术可以有各种更改和变化。凡在本技术的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本技术的保护范围之内。
[0124]
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。