申请/专利权人：昆明理工大学

申请日：2018-06-20

公开（公告）日：2020-09-25

公开（公告）号：CN108919029B

主分类号：G01R31/52(20200101)

分类号：G01R31/52(20200101);G01R29/16(20060101);G01R31/34(20060101);F03D17/00(20160101);G06F17/13(20060101)

优先权：

专利状态码：有效-授权

法律状态：2020.09.25#授权;2018.12.25#实质审查的生效;2018.11.30#公开

摘要：本发明涉及一种计及RSC控制的双馈风机三相短路电流的解析方法，属风力发电系统故障分析领域。本方法为：列写DFIG定子和转子电压方程、定子和转子磁链方程、RSC控制的转子电压方程；确定故障发生后定子磁链解析式；求取转子电流二阶微分方程；再求解转子电流二阶微分方程，得到计及RSC控制的DFIG转子短路电流解析式；再求得定子电流方程，将定子磁链解析式以及转子电流解析式代入定子电流方程，即求得计及RSC控制的DFIG定子短路电流解析式。本发明能够准确计算计及转子侧变流器的双馈感应风力发电机三相短路电流的解析式，对含双馈感应风力发电机的电力系统设备选型和保护动作特性分析具有重要意义。

主权项：1.一种计及RSC控制的双馈风机三相短路电流的解析方法，其特征在于：所述方法包括如下步骤：步骤一：根据双馈感应风力发电机等效电路图列写定子和转子电压方程，根据双馈感应风力发电机定子和转子的自感磁链和互感磁链列写定子和转子磁链方程，根据RSC控制原理列写转子电压方程，以上所列写的方程均在dq同步旋转坐标系下；步骤二：根据步骤一中所列写的定子电压方程以及磁链守恒原则，求取故障后定子磁链解析式；步骤三：根据步骤一中所列写的定子和转子磁链方程，即可求得转子磁链与转子电流和定子磁链相关的方程；再将步骤二中所求的定子磁链解析式代入所求得的转子磁链与转子电流和定子磁链相关的方程中，求得与转子电流相关的转子磁链方程；将所求得的与转子电流相关的转子磁链方程代入步骤一中根据双馈感应风力发电机等效电路图列写的转子电压方程中，可得与转子电流相关的转子电压方程；最后将所求得的与转子电流相关的转子电压方程代入步骤一中根据RSC控制原理列写的转子电压方程中，即可求得转子电流二阶微分方程；步骤四：首先求解转子电流二阶微分方程对应的齐次方程的特征方程的根，得到转子电流二阶微分方程对应的齐次方程通解的形式，再根据转子电流二阶微分方程的初值求得转子电流二阶微分方程对应的齐次方程解的表达式；接着根据待定系数法求解转子电流二阶微分方程的特解；最后，将所求的转子电流二阶微分方程对应的齐次方程解的表达式与转子电流二阶微分方程的特解相加，所得的和即为dq同步旋转坐标系下转子短路电流解析式；步骤五：根据步骤一中所列写的定子磁链方程，可求得定子电流方程；将步骤二中所求得的定子磁链解析式与步骤四中所求得的转子短路电流解析式代入定子电流方程中，即可求得dq同步旋转坐标系下定子短路电流解析式；步骤六：将步骤四中所求得的dq同步旋转坐标系下转子短路电流和步骤五中的所求得的dq同步旋转坐标系下定子短路电流变换到三相静止坐标系下，得到三相静止坐标系下定子短路电流和转子短路电流的解析式。

全文数据：一种计及RSC控制的双馈风机三相短路电流的解析方法技术领域[0001]本发明涉及一种计及RSC控制的双馈风机三相短路电流的解析方法，属于风力发电系统故障分析技术领域。背景技术[0002]随着全球经济的发展和能源消耗量的大幅度增长，能源的储量、生产和使用之间的矛盾日益突出，成为目前世界各国急待解决的重要问题之一。因此，为解决能源危机、环境污染等问题，风能、太阳能等新能源的研究开发已成当前人类十分迫切的需求。其中，风能是一种清洁永续的能源，与传统能源相比，风力发电具有不依赖外部能源、没有燃料价格风险、发电成本稳定、没有碳排放等环境成本特点；与太阳能、潮汐能相比，风能的产业基础最好，经济优势最为明显，没有大的环境影响；而且，全球范围内可利用的风能分布十分广泛。由于风力发电具有的这些独特优势，使其逐渐成为许多国家可持续发展战略的重要组成部分，并在世界各国得到迅速发展。[0003]风力发电机组种类较多，双馈感应风力发电机通过使用双PWM变流器控制其励磁电流实现了发电机组与风力系统良好的柔性连接，这种连接方式便于并网操作，具有有功、无功功率独立控制，可变速运行及励磁变流器容量小等优点，因此成为了风电场MW级风力发电机的主要机型。但另一方面并网型双馈感应风力发电机组在并网电压突降时的暂态特性相当复杂，不同于传统的同步和异步电机。[0004]当风电大规模接入系统后，变压器，线路阻抗器以及断路器等电气设备的动、热稳定性校验，以及线路、变压器等各元件的保护动作特性主要依靠系统的短路电流计算整定，因此随着风机大规模的并网，确定双馈感应发电机在故障过程中的短路电流特性是目前双馈风机并网需解决的重要问题。发明内容[0005]本发明提供了一种计及RSC控制的双馈风机三相短路电流的解析方法，利用DFIG定子和转子电压方程、定子和转子磁链方程、RSC控制的转子电压方程，计算了计及RSC转子侧变流器控制后的DFIG三相短路电流。[0006]本发明的技术方案是:一种计及RSC控制的双馈风机三相短路电流的解析方法，所述方法包括如下步骤：[0007]步骤一:根据双馈感应风力发电机简称双馈风机等效电路图（图1列写定子和转子电压方程，根据双馈感应风力发电机定子和转子的自感磁链和互感磁链列写定子和转子磁链方程，根据RSC控制原理列写转子电压方程，以上所列写的方程均在dq同步旋转坐标系下；[0008]步骤二:根据步骤一中所列写的定子电压方程以及磁链守恒原则，求取故障后定子磁链解析式；[0009]步骤三:根据步骤一中所列写的定子和转子磁链方程，即可求得转子磁链与转子电流和定子磁链相关的方程;再将步骤二中所求的定子磁链解析式代入所求得的转子磁链与转子电流和定子磁链相关的方程中，求得与转子电流相关的转子磁链方程;将所求得的与转子电流相关的转子磁链方程代入步骤一中根据图1列写的转子电压方程中，可得与转子电流相关的转子电压方程;最后将所求得的与转子电流相关的转子电压方程代入步骤一中根据RSC控制原理列写的转子电压方程中，即可求得转子电流二阶微分方程；[0010]步骤四：首先求解转子电流二阶微分方程对应的齐次方程的特征方程的根，得到转子电流二阶微分方程对应的齐次方程通解的形式，再根据转子电流二阶微分方程的初值求得转子电流二阶微分方程对应的齐次方程解的表达式;接着根据待定系数法求解转子电流二阶微分方程的特解;最后，将所求的转子电流二阶微分方程对应的齐次方程解的表达式与转子电流二阶微分方程的特解相加，所得的和即为dq同步旋转坐标系下转子短路电流解析式；[0011]步骤五:根据步骤一中所列写的定子磁链方程，可求得定子电流方程;将步骤二中所求得的定子磁链解析式与步骤四中所求得的转子电流电流解析式代入定子电流方程中，即可求得dq同步旋转坐标系下定子短路电流解析式。[0012]步骤六:将步骤四中所求得的dq同步旋转坐标系下转子短路电流和步骤五中的所求得的dq同步旋转坐标系下定子短路电流变换到三相静止坐标系下，得到三相静止坐标系下定子短路电流和转子短路电流的解析式。[0013]1、所述步骤一的具体步骤为：[0014]根据双馈感应风力发电机等效电路图（图1列写系统发生三相短路时，dq同步旋转坐标系下定子和转子电压方程为：[0017]式中，U^DFIG机端的稳态电压;k为机端电压跌落系数;Ur为转子电压；is为DFIG定子三相短路电流；ir为DFIG转子三相短路电流;队为定子磁链;知为转子磁链;Rs为定子绕组等效电阻;Rr为转子绕组等效电阻；O1为同步转速;S=CO1-COrCO1为转差率，COr为转子转速；为虚单位;t为时间；[0018]根据双馈感应风力发电机定子和转子的自感磁链和互感磁链列写dq同步旋转坐标系下定子和转子磁链方程为：[0021]式中，Ls=Lm+Ljs为定子绕组等效电感，Lm为定转子绕组互感，Lcjs为定子漏感;Lr=Lm+Lw为转子绕组等效电感，Lw为转子漏感；[0022]根据转子侧变流器RSC的控制原理，可得dq同步旋转坐标系下转子电压方程为：[0024]式中，kp、ki分别为PI控制器的比例参数和积分参数;ir.rrf为转子电流参考值;ir为DFIG转子电流；=Μ。.-总，为定转子等效电感。[0025]2、求定子磁链解析式，所述步骤二的具体步骤为：[0026]假设Os时刻系统发生三相短路对称故障，根据磁链守恒原则，发生故障瞬间，DFIG机端电压由us骤降至（1-kus，但定子磁链不会发生突变，因此故障后的定子磁链有两个分量:一是与机端残压（l_kUs相对应的定子磁链稳态分量itsf;二是与电压跌落部分相对应的定子磁链暂态分量Φπ，该分量以定子时间常数衰减；[0027]1、求取发生故障后定子磁链稳态分量；[0028]系统发生三相短路对称故障后DFIG的机端电压由Us跌落至（1-kUs，稳态时忽略定子电阻以及磁链变换率的影响，由式⑴可得在短路发生后与机端残压（1-kUs相对应的定子磁链稳态分量为：[0030]式中，为定子磁链稳态分量;k为系统发生三相短路故障后DFIG机端电压的跌落率；[0031]2、求取发生故障后定子磁链暂态分量；[0032]根据磁链守恒原则，定子磁链暂态分量为：[0034]式中，I为定子磁链暂态分量;is=RsLrLD为定子衰减时间常数；[0035]3、由以上分析可知系统发生三相短路对称故障后，dq旋转坐标下定子磁链^为：[0037]3、求取转子电流二阶微分方程，所述步骤三的具体步骤为：[0038]1、联立式（3、式4所示的定子和转子磁链方程，消去定子电流is可得转子磁链与转子电流和定子磁链相关的方程为：[0040]2、将式8定子磁链解析式代入式9转子电流方程，可进一步得到转子磁链的方程为：[0042]3、将式10代入式⑵转子电压方程，可得转子电压与转子电流之间的关系为：[0044]⑷、对式11方程两边同时微分可得：[0045][0046]5、对式⑸RSC控制的转子电压方程两边同时微分可得：[0048]⑹、将式13代入式12，可得转子电流的二阶微分方程为：[0050]式中，ft=Rr+kpLsLd[0051]02=kiLsLD[0052]βι·=jsωi_jωi-Tsjωi+tsLmLs。[0053]4、解转子电流二阶微分方程，求得转子短路电流解析式，所述步骤四的具体步骤为：[0054]通过求解转子电流二阶微分方程式（14，其解即为转子电流解析式，转子电流二阶微分方程的解为其对应的齐次方程的通解与转子电流二阶微分方程的特解之和，具体求解过程如下：[0055]1求转子二阶微分方程对应的齐次方程的通解ireai的表达式；[0056]1求转子二阶微分方程对应的齐次方程的通解ireai的形式；[0057]转子二阶微分方程式14对应的齐次方程为：[0059]转子二阶微分方程式对应的齐次方程式15的特征方程为：[0061]其中，ar为转子特征方程的特征值；[0062]对式（16进行求解，可求得特征方程式（16的特征根为Ct1，a2两个不相等的实数根：[0065]式中，ai、a2为转子特征方程两个不相等的特征根；[0066]根据二阶常系数齐次微分方程的通解的定义可知，当所求的特征方程有两个不相等的实数根时，则是转子二阶微分方程式（15两个线性无关的特解，故所求转子二阶微分方程式15的通解iresii的表不为：[0068]式中，ir纖D为转子二阶微分方程的通解;C1、C2为待求常数；[0069]2求转子二阶微分方程对应的齐次方程的通解ireai的表达式；[0070]转子二阶微分方程对应的齐次方程的通解ireai的形式中常数C1X2的值是根据故障瞬间转子电流的初值确定的，故障瞬间Os时刻转子电流初值为稳态运行时转子电流的指令值，即irQIt=Q=ir.ref;故障瞬间转子电流初值的一阶导为iIt=Q=0;[0071]对所求转子二阶微分方程对应的齐次方程通解式17求一阶导得：[0073]令式（17中t=0、ir_D=ir.ref代入式（17可得：[0074]ir.ref=Cl+C219[0075]令式18中ΐ=0、;ΐ：Γ_=0代入式18可得：[0076]〇=aiCi+a2C220[0077]联立方程式19、式20后求解可得常数C1X2的值分别为：[0079]将式（21代入式（17，可求得转子二阶微分方程对应的齐次方程通解的表达式为：[0081]2求转子二阶微分方程的特解；[0082]令二阶微分方程为：[0085]设匕、4分别为二阶微分方程式23、式24的特解，则转子二阶微分方程的特解为：[0087]其中，ΐ；：为转子二阶微分方程的特解；[0088]1求二阶微分方程式23的特解;[0089]转子电流二阶微分方程对应的齐次方程的特征根为两个不相等的实数根，故可设二阶微分方程式23的特解‘为：[0091]式中，bi为待求常数；[0092]对式26求一阶导和二阶导得：[0095]将式26—式28代入式21中，可得：[0097]对式29进行求解，可得常数h的值为In=Lef;将所求的匕的值代入式26，即可求得二阶微分方程式23的特解匕的值为：[0099]2求二阶微分方程式24的特解2:;[0100]转子电流二阶微分方程对应的齐次方程的特征根为两个不相等的实数根，故可设二阶微分方程式24的特解:2为：[0102]式中，b2为待求常数；[0103]对式31求一阶导和二阶导得：[0106]将式31—式33代入式24，可得：[0108]对式34进行求解，可得b2的值为：[0110]将式35代入式31所设的特解匕中，可得二阶微分方程式24特解4的值为：[0112]3转子二阶微分方程的特解[0113]将所求的二阶微分方程式23的特解C1和式24的特解:2代入式25中，即可求得转子二阶微分方程的特解为：[0115]⑷求取dq同步旋转坐标系下转子短路电流的解析式[0116]转子短路电流的解析式为转子二阶微分方程对应的齐次方程的通解ireai和转子二阶微分方程的特解之和，即为：[0118]式中，Ar=ir.ref[0120]Cr=ir.refCX2-CX1。[0121]5、求取dq同步旋转坐标系下定子电流解析式，所述步骤五的具体步骤为：[0122]由定子磁链方程式⑶可得定子电流方程为：[0124]将定子磁链解析式式8和转子电流解析式式38代入定子电流方程式39，可得dq同步旋转坐标系下定子短路电流解析式为：[0126]式中，[0129]6、求取三相静止坐标系下定子和转子短路电流[0130]所述步骤六的具体步骤为：[0131]将所求得的dq旋转坐标系下DFIG定子和转子短路电流转换到三相静止坐标系下为：[0134]式中，ira、irb、irc分别为三相静止坐标系下转子a、b、c各相电流；isa、isb、isc分别为三相静止坐标系下定子A、B、C各相电流;Re为取实部符号；[0135]将式38代入式41可得计及转子侧变流器的双馈感应风力发电机的转子三相短路电流的解析式为：[0137]将式39代入式42可得计及转子侧变流器的双馈感应风力发电机的定子三相短路电流的解析式为：[0139]本发明的有益效果是：[0140]1、本发明在dq同步旋转坐标系下求解计及RSC控制的DFIG三相短路电流的解析式，再通过反Park变换得到三相静止坐标系下计及RSC控制的DFIG三相短路电流的解析式。在三相静止坐标系下DFIG磁链方程是非线性的代数方程，电压方程是时变系数的微分方程，对于求解三相短路电流带来了很大的困难;dq旋转坐标系下磁链方程变成线性代数方程，电压方程为常微分方程，Park变换简化了复杂的电压和磁链方程。[0141]2、本发明采用全量法求解计及RSC控制的DFIG三相短路电流的解析式，全量法与故障分量法相比较，故障分量法求得的短路电流故障分量需要与故障之前的稳态量进行叠加才能得到短路电流全解，但全量法能够一次求解就得到结果，不需要进行叠加。附图说明[0142]图1为双馈感应风力发电机等效电路图；[0143]图2为DFIG转子三相短路电流随时间的变化曲线图；[0144]图3为DFIG定子三相短路电流随时间的变化曲线图；[0145]图4为双馈感应风力发电机的工作原理图；[0146]图5为DFIG三相静止坐标系下的物理模型；[0147]图6为DFIG转子侧变频器控制框图。具体实施方式[0148]实施例1:如图1-6所示，一种计及RSC控制的双馈风机三相短路电流的解析方法，包括如下步骤：[0149]步骤一:根据双馈感应风力发电机等效电路图（图1列写定子和转子电压方程，根据双馈感应风力发电机定子和转子的自感磁链和互感磁链列写定子和转子磁链方程，根据RSC控制原理列写转子电压方程，以上所列写的方程均在dq同步旋转坐标系下；[0150]步骤二:根据步骤一中所列写的定子电压方程以及磁链守恒原则，求取故障后定子磁链解析式；[0151]步骤三:根据步骤一中所列写的定子和转子磁链方程，即可求得转子磁链与转子电流和定子磁链相关的方程;再将步骤二中所求的定子磁链解析式代入所求得的转子磁链与转子电流和定子磁链相关的方程中，可求得与转子电流相关的转子磁链方程;将所求得的与转子电流相关的转子磁链方程代入步骤⑴中根据图1列写的转子电压方程中，可得与转子电流相关的转子电压方程;最后将所求得的与转子电流相关的转子电压方程代入步骤一中根据RSC控制原理列写的转子电压方程中，即可求得转子电流二阶微分方程；[0152]步骤四：首先求解转子电流二阶微分方程对应的齐次方程的特征方程的根，得到转子电流二阶微分方程对应的齐次方程通解的形式，再根据转子电流二阶微分方程的初值求得转子电流二阶微分方程对应的齐次方程解的表达式;接着根据待定系数法求解转子电流二阶微分方程的特解;最后，将所求的转子电流二阶微分方程对应的齐次方程解的表达式与转子电流二阶微分方程的特解相加，所得的和即为dq同步旋转坐标系下转子短路电流解析式；[0153]步骤五:根据步骤一中所列写的定子磁链方程，可求得定子电流方程;将步骤二中所求得的定子磁链解析式与步骤四中所求得的转子电流解析式代入定子电流方程中，即可求得dq同步旋转坐标系下定子短路电流解析式。[0154]步骤六:将步骤四中所求得的dq同步旋转坐标系下转子短路电流和步骤五中的所求得的dq同步旋转坐标系下定子短路电流变换到三相静止坐标系下，得到三相静止坐标系下定子短路电流和转子短路电流的解析式。[0155]本发明为一种计及RSC控制的双馈风机三相短路电流的解析方法，按以下步骤进行：[0156]以一台并网双馈感应风力发电机为例，假设在Os时刻系统发生三相短路故障，双馈感应风力发电机的具体参数如表1所示：[0157]表1双馈感应发电机的参数设置[0158][0159]1、根据双馈感应风力发电机等效电路图（图1列写dq同步旋转坐标系下定子和转子电压方程为：[0162]式中，uADFIG机端的稳态电压;k为机端电压跌落系数;ur为转子电压；iADFIG定子三相短路电流；ir为DFIG转子三相短路电流;队为定子磁链;也为转子磁链;私为定子绕组等效电阻;Rr为转子绕组等效电阻；W1为同步转速；COr为转子转速;S=CO1-COrZiCO1S转差率；为虚单位;t为时间。[0163]根据双馈感应风力发电机定子和转子的自感磁链和互感磁链列写dq同步旋转坐标系下定子和转子磁链方程为：[0166]式中，Lm为定转子绕组互感;Ls=Lm+Ljs为定子绕组等效电感，Ljs为定子漏感;Lr=Lm+Lw为转子绕组等效电感，Lw为转子漏感。[0167]根据转子侧变流器RSC的控制原理，可得dq同步旋转坐标系下转子电压方程为：[0169]式中，kp、ki分别为PI控制器的比例参数和积分参数；ir.ref为转子电流参考值矢量;ir为DFIG转子电流；G=I.0148,为定转子等效电感。[0170]2、求定子磁链解析式[0171]假设Os时刻系统发生三相短路对称故障，根据磁链守恒原则，发生故障瞬间，DFIG机端电压由us骤降至（1-kus，但定子磁链不会发生突变。因此故障后的定子磁链有两个分量:一是与机端残压（l_kUs相对应的定子磁链稳态分量itsf;二是与电压跌落部分相对应的定子磁链暂态分量Φπ，该分量以定子时间常数衰减。[0172]1求取发生故障后定子磁链稳态分量[0173]系统发生三相短路对称故障后DFIG的机端电压由Us跌落至（1-kUs。稳态时忽略定子电阻以及磁链变换率的影响，由式⑴可得在短路发生后与机端残压（1-kUs相对应的定子磁链稳态分量为：[0175]式中，k为系统发生三相短路故障后DFIG机端电压的跌落率。[0176]2求取发生故障后定子磁链暂态分量[0177]根据磁链守恒原则，定子磁链暂态分量为：[0179]式中，φ3η为定子磁链暂态分量;Ts=RsLrLD=0.0694为定子衰减时间常数。[0180]3由以上分析可知系统发生三相短路对称故障后，dq旋转坐标下定子磁链也为：[0182]3、求取转子电流二阶微分方程[0183]1联立式3、式4所示的定子和转子磁链方程，消去定子电流is可得转子磁链与转子电流和定子磁链相关的方程为：[0185]2将式8定子磁链解析式代入式9转子电流方程，可进一步得到转子磁链的[0187]3将式10代入式⑵转子电压方程，可得转子电压与转子电流之间的关系为：[0188][0189]⑷对式（11方程两边同时微分可得：[0191]5对式⑸RSC控制的转子电压方程两边同时微分可得：[0193]⑹将式13代入式12，可得转子电流的二阶微分方程为：[0195]式中，β1=Rr+kpLsLD=0.3521[0196]02=kiLsLD=5.O595[0197]βΓ=jsQ1-J'ωι-τ3jo1+TsLmLs=3.4115-j0.436[0198]4、解转子电流二阶微分方程，求得转子短路电流解析式[0199]通过求解转子电流二阶微分方程式（14，其解即为转子电流解析式。转子电流二阶微分方程的通解为其对应的齐次方程的通解与转子电流二阶微分方程的特解之和，具体求解过程如下：[0200]1求转子二阶微分方程对应的齐次方程的通解ireai的解析式[0201]1求转子二阶微分方程对应的齐次方程的通解ireai*的形式[0202]转子二阶微分方程式14对应的齐次方程为：[0204]转子二阶微分方程式对应的齐次方程式15的特征方程为：[0206]对式（16进行求解，可求得特征方程式（16的特征根为Ct1，α2两个不相等的实数根：[0209]根据二阶常系数齐次微分方程的通解的定义可知，当所求的特征方程有两个不相等的实数根时，则L=Pi、是方程（15的两个线性无关的特解，故所求方程式（15的通解ir纖®的形式可以表示为：[0211]式中，C1X2为常数。[0212]2求转子二阶微分方程对应的齐次方程的通解ireai的表达式[0213]转子二阶微分方程对应的齐次方程的通解ireai的形式中常数Cl、C2的值是根据故障瞬间转子电流的初值确定的。故障瞬间Os时刻转子电流初值为稳态运行时转子电流的指令值，即ilflIt=Q=ir.rrf;故障瞬间转子电流初值的一阶导为iIt=Q=0。[0214]对所求转子二阶微分方程对应的齐次方程通解式17求一阶导得：[0216]令式（17中t=0、ir_D=ir.ref代入式（17可得：[0217]ir.ref=Cl+C219[0218]令式（18中ΐ=0、;ΐ：Γ_=0代入式（18可得：[0219]〇=aiCi+a2C220[0220]联立方程式19、式20后求解可得常数C1X2的值分别为：[0222]将式（21代入式（17，可求得转子二阶微分方程对应的齐次方程通解的表达式为：[0224]2求转子二阶微分方程的特解[0225]令二阶微分方程为：[0228]设匕、匕分别为二阶微分方程式23、式24的特解，则转子二阶微分方程的特解为：[0230]其中，i,!为转子二阶微分方程的特解[0231]1求二阶微分方程式23特解匕[0232]转子电流二阶微分方程对应的齐次方程的特征根为两个不相等的实数根，故可设二阶微分方程式23的特解为：26[0234]式中，bi为待求常数。[0235]对式26求一阶导和二阶导得：[0238]将式26—式28代入式23中，可得：[0240]对式29进行求解，可得常数h的值为h=ir.ref=0.8-j0.25;将所求的brl的值代入式26，即可求得二阶微分方程式23特解匕的值为：[0242]2求二阶微分方程式24特解4[0243]转子电流二阶微分方程对应的齐次方程的特征根为两个不相等的实数根，故可设二阶微分方程式24的特解4为：[0245]式中，132为常数。[0246]对式31求一阶导和二阶导得：[0249]将式31—式33代入式24，可得：[0251]对式34进行求解，可得b2的值为：[0253]将式35代入式31所设的特解匕中，可得二阶微分方程式24特解4的值为：[0255]3转子二阶微分方程的特解[0256]将所求的二阶微分方程式23的特解匕和式24的特解C2代入式25中，即可求得转子二阶微分方程的特解为；[0258]3求取dq同步旋转坐标系下转子短路电流的解析式[0259]转子短路电流的解析式为转子二阶微分方程对应的齐次方程的通解ireai和转子二阶微分方程的特解之和，即为：[0261]式中，Ar=ir.ref=〇.8_j0.25[0263]Cr=ir.refa2-〇i=O.0557+j0.1784[0264]5、求取dq同步旋转坐标系下定子电流解析式[0265]由定子磁链方程式⑶可得定子电流方程为：[0267]将定子磁链解析式式8和转子电流解析式式38代入定子电流方程式39，可得dq同步旋转坐标系下定子短路电流解析式为：[0272]6、求取三相静止坐标系下定子和转子短路电流[0273]将所求得的dq旋转坐标系下DFIG定子和转子短路电流转换到三相静止坐标系下为：[0276]式中，ira、irb、ire分别为转子a、b、c各相电流；isa、isb、isc分别为定子A、B、C各相电流;Re为取实部符号。[0277]将式38代入式41可得计及转子侧变流器的双馈感应风力发电机的转子三相短路电流的解析式为：[0279]将式39代入式42可得计及转子侧变流器的双馈感应风力发电机的定子三相短路电流的解析式为：[0281]根据式43可以绘制出转子三相短路电流随时间的变化趋势图，如图2所示；[0282]根据式44可以绘制出定子三相短路电流随时间的变化趋势图，如图3所示。本发明原理：[0283]1、双馈感应风力发电机的工作原理[0284]双馈风力发电机由风力机，齿轮箱，双馈电机，背靠背变流器及控制系统四部分组成。如图4所示，定子绕组直接与电网相连;转子绕组则通过背靠背变频器与电网相连;齿轮箱是将在风力作用下的风轮所产生的动力传递给发电机并使其得到相应的转速;背靠背变流器PWM是由两个电压源变流器构成的交-直-交电压源变流器，与转子侧相连的称为转子侧变流器Rotor-sideConverter，RSC，与定子侧相连的称为定子侧变流器Grid-sideC〇nverter，GSC，转子侧变流器RSC和定子侧变流器GSC之间通过一个直流电压相连；控制系统主要是控制定子和转子侧变流器。[0285]双馈感应风力发电机工作原理可以简述为:风速发生变化时，即转子的转速ω1随之变化，通过控制转子励磁电流的频率f2,即控制转子励磁电流的转速ω2,使得气隙合成磁场相对于定子转速保持不变，保证定子转速为同步转速，DFIG实现变速恒频运行。[0286]2、双馈感应风力发电机数学模型的建立[0287]双馈风力发电机的数学模型是高阶、非线性、强耦合的多变量系统。为建立数学模型，一般做以下假设：[0288]1定子和转子绕组三相均对称，在空间上角度互差120°，定子和转子电流中只考虑基波分量，谐波分量忽略不计，空间磁动势沿气隙圆周呈正弦分布；[0289]2忽略磁路饱和、涡流损耗和铁耗；[0290]3忽略电机参数受温度和频率变化的影响；[0291]⑷定子和转子绕组的自感和互感恒定不变[0292]基于以上分析假设，同时规定DFIG定转子侧均采用电动机惯例，则可建立DFIG等效物理模型，如图5的所示，A、B、C表示定子三相绕组轴线，在空间上是静止的;a、b、c表示转子三相绕组轴线，在空间上以转子转速ωr旋转的；0r=ωrt表示转子a轴与定子A轴之间的角·[0293]根据图5双馈电机的物理模型，可以列写出三相静止坐标系下DFIG的数学模型：[0296]式中：ua、ub、uc为定子A、B、C三相电压，iA、iB、ic为定子A、B、C三相电流;ua、ub、uc为转子a、b、c三相电压，ia、ib、ic为转子a、b、c三相电流;Hik为定子A、B、C三相绕组的全磁链，为转子a、b、c三相绕组的全磁链;Rs、Rr分别为定子绕组和转子绕组的电阻;D表示微分算子[0297]定、转子各相绕组的全磁链是由各相绕组的自电感磁链和其它绕组的互感磁链组成，定、转子磁链方程如式47所示：[0299]式中：定子各相自感为:LM=LBB=Lcc=Lms+Ljs[0300]转子各相自感为:Laa=Lbb=Lcc=Lmr+Lor[0301]Us、!^分别表示定子绕组每相的主电感和漏感；[0302]Lmr、L。^别表示转子绕组每相的主电感和漏感；[0303]定子三相彼此之间的互感为：[0305]转子三相彼此之间的互感为：[0307]转子与定子之间的互感：[0309]由于双馈感应发电机定子绕组与转子饶组之间的互感为时变函数，从而磁链方程式47将是非线性的代数方程。将磁链方程代入电压方程后，电压方程式45和式46将成为时变系数的微分方程，给分析带来了很大的困难。所以一般会对三相静止坐标系下的电压、磁链方程进行派克Park变换，从而得到dq旋转坐标系下的电压、磁链方程，在dq旋转坐标系下的磁链方程变成线性代数方程，电压方程变成了常微分方程。[0310]对三相静止坐标系下DFIG电压和磁链方程进行Park变换，得到dq旋转坐标系下DFIG电压和磁链方程为：[0315]式中，Us为DFIG机端的稳态电压;k为机端电压跌落系数;Ur为转子电压；is为DFIG定子三相短路电流；ir为DFIG转子三相短路电流;队为定子磁链;知为转子磁链;Rs为定子绕组等效电阻;Rr为转子绕组等效电阻；W1为同步转速；COr为转子转速;S=ωι-ωτωι*转差率；为虚单位;t为时间。[0316]2、转子侧变流器控制原理[0317]目前转子侧变流器采用的控制方式主要有矢量控制和直接功率控制，不同的控制方式产生的控制效果有一定的差别，本文转子侧变流器控制采用了定子磁链定向控制原理。根据定子磁链定向原理，转子电流的参考值为：[0320]式中，ird.ref、irq.ref分别为转子d、q轴电流分量的参考值;Φιη为定子磁链幅值；PS.rrf、QS.rrf分别为定子有功、无功的参考值;Us为定子电压。[0321]DFIG的控制原理如图6所示。图6中转子电流的参考值ird.ref、irq.ref与转子电流ird、irq的反馈值比较后的差值送入PI控制器，输出电压IZrtKl^rq与电压补偿分量AUrd、AUrq叠加，就可以获得转子电压指令值Urd.rAUrq.rrf，经过坐标变换后得到的转子电压进行脉宽调制后输出对转子侧变频器的驱动信号，实现对DFIG的控制。[0322]图6中，Urd.ref、Urq.ref分另Ij为跟踪转子电流戶斤需要的转子电压参考值。转子侧输出的转子电压方程为：[0324]式中，kp、kl分别为PI控制器的比例参数和积分参数;ir.ref为转子电流参考值;ir为DFIG转子电流;4,为定转子等效电感。[0325]上面结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。

权利要求：1.一种计及RSC控制的双馈风机三相短路电流的解析方法，其特征在于:所述方法包括如下步骤：步骤一:根据双馈感应风力发电机等效电路图列写定子和转子电压方程，根据双馈感应风力发电机定子和转子的自感磁链和互感磁链列写定子和转子磁链方程，根据RSC控制原理列写转子电压方程，以上所列写的方程均在dq同步旋转坐标系下；步骤二:根据步骤一中所列写的定子电压方程以及磁链守恒原则，求取故障后定子磁链解析式；步骤三:根据步骤一中所列写的定子和转子磁链方程，即可求得转子磁链与转子电流和定子磁链相关的方程;再将步骤二中所求的定子磁链解析式代入所求得的转子磁链与转子电流和定子磁链相关的方程中，求得与转子电流相关的转子磁链方程;将所求得的与转子电流相关的转子磁链方程代入步骤一中根据列写的转子电压方程中，可得与转子电流相关的转子电压方程；最后将所求得的与转子电流相关的转子电压方程代入步骤一中根据RSC控制原理列写的转子电压方程中，即可求得转子电流二阶微分方程；步骤四：首先求解转子电流二阶微分方程对应的齐次方程的特征方程的根，得到转子电流二阶微分方程对应的齐次方程通解的形式，再根据转子电流二阶微分方程的初值求得转子电流二阶微分方程对应的齐次方程解的表达式;接着根据待定系数法求解转子电流二阶微分方程的特解;最后，将所求的转子电流二阶微分方程对应的齐次方程解的表达式与转子电流二阶微分方程的特解相加，所得的和即为dq同步旋转坐标系下转子短路电流解析式；步骤五:根据步骤一中所列写的定子磁链方程，可求得定子电流方程;将步骤二中所求得的定子磁链解析式与步骤四中所求得的转子电流电流解析式代入定子电流方程中，即可求得dq同步旋转坐标系下定子短路电流解析式。步骤六:将步骤四中所求得的dq同步旋转坐标系下转子短路电流和步骤五中的所求得的dq同步旋转坐标系下定子短路电流变换到三相静止坐标系下，得到三相静止坐标系下定子短路电流和转子短路电流的解析式。2.根据权利要求1所述的计及RSC控制的双馈风机三相短路电流的解析方法，其特征在于:所述步骤一的具体步骤为：根据双馈感应风力发电机等效电路图列写系统发生三相短路时，dq同步旋转坐标系下定子和转子电压方程为：式中，Us为DFIG机端的稳态电压;k为机端电压跌落系数;Ur为转子电压；is为DFIG定子三相短路电流；ir为DFIG转子三相短路电流;队为定子磁链;屯为转子磁链;Rs为定子绕组等效电阻;Rr为转子绕组等效电阻；CO1为同步转速；s=CO1-COrCO1为转差率，COr为转子转速；％虚单位;t为时间；根据双馈感应风力发电机定子和转子的自感磁链和互感磁链列写dq同步旋转坐标系下定子和转子磁链方程为：式中，Ls=Lm+Ljs为定子绕组等效电感，Lm为定转子绕组互感，Lcjs为定子漏感;Lr=Lm+L。!·为转子绕组等效电感，Lw为转子漏感；根据转子侧变流器RSC的控制原理，可得dq同步旋转坐标系下转子电压方程为：式中，kP、k4V别为PI控制器的比例参数和积分参数；ir.ref为转子电流参考值；ir为DFIG转子电流，为定转子等效电感。3.根据权利要求2所述的计及RSC控制的双馈风机三相短路电流的解析方法，其特征在于:所述步骤二的具体步骤为：假设Os时刻系统发生三相短路对称故障，根据磁链守恒原则，发生故障瞬间，DFIG机端电压由us骤降至（1-kus，但定子磁链不会发生突变，因此故障后的定子磁链有两个分量:一是与机端残压（1-kUs相对应的定子磁链稳态分量itsf;二是与电压跌落部分相对应的定子磁链暂态分量Φπ，该分量以定子时间常数衰减；1、求取发生故障后定子磁链稳态分量；系统发生三相短路对称故障后DFIG的机端电压由Us跌落至（1-kus，稳态时忽略定子电阻以及磁链变换率的影响，由式⑴可得在短路发生后与机端残压（1-kUs相对应的定子磁链稳态分量为：式中，为定子磁链稳态分量;k为系统发生三相短路故障后DFG机端电压的跌落率；2、求取发生故障后定子磁链暂态分量；根据磁链守恒原则，定子磁链暂态分量为：式中，I为定子磁链暂态分量;ts=RsLrLD为定子衰减时间常数；3、由以上分析可知系统发生三相短路对称故障后，dq旋转坐标下定子磁链也为：4.根据权利要求3所述的计及RSC控制的双馈风机三相短路电流的解析方法，其特征在于:所述步骤三的具体步骤为：1、联立式3、式4所示的定子和转子磁链方程，消去定子电流is可得转子磁链与转子电流和定子磁链相关的方程为：2、将式8定子磁链解析式代入式9转子电流方程，可进一步得到转子磁链的方程为：3、将式10代入式⑵转子电压方程，可得转子电压与转子电流之间的关系为：⑷、对式11方程两边同时微分可得：5、对式⑸RSC控制的转子电压方程两边同时微分可得：6、将式13代入式12，可得转子电流的二阶微分方程为：5.根据权利要求4所述的计及RSC控制的双馈风机三相短路电流的解析方法，其特征在于:所述步骤四的具体步骤为：通过求解转子电流二阶微分方程式（14，其解即为转子电流解析式，转子电流二阶微分方程的解为其对应的齐次方程的通解与转子电流二阶微分方程的特解之和，具体求解过程如下：1求转子二阶微分方程对应的齐次方程的通解ir纖®的表达式；1求转子二阶微分方程对应的齐次方程的通解ir纖®的形式；转子二阶微分方程式14对应的齐次方程为：转子二阶微分方程式对应的齐次方程式15的特征方程为：其中，ar为转子特征方程的特征值；对式16进行求解，可求得特征方程式16的特征根为αι，α2两个不相等的实数根：式中，αι、α2为转子特征方程两个不相等的特征根；根据二阶常系数齐次微分方程的通解的定义可知，当所求的特征方程有两个不相等的实数根时，则是转子二阶微分方程式（15两个线性无关的特解，故所求转子二阶微分方程式（15的通解ir_D的表不为：式中，irea®为转子二阶微分方程的通解;C1、C2为待求常数；2求转子二阶微分方程对应的齐次方程的通解irel»的表达式；转子二阶微分方程对应的齐次方程的通解irea®的形式中常数C1X2的值是根据故障瞬间转子电流的初值确定的，故障瞬间Os时刻转子电流初值为稳态运行时转子电流的指令值，即irOIt=Q=ir.:ref;故障瞬间转子电流初值的一阶导为i’ιΌIt=Q=0;对所求转子二阶微分方程对应的齐次方程通解式17求一阶导得：令式（17中t=0、ir_D=ir.ref代入式（17可得：令式18中t=0、i’r_=0代入式18可得：联立方程式19、式20后求解可得常数、：2的值分别为：将式21代入式17，可求得转子二阶微分方程对应的齐次方程通解的表达式为：⑵求转子二阶微分方程的特解；令二阶微分方程为：设G、2分别为二阶微分方程式23、式24的特解，则转子二阶微分方程的特解为：其中，为转子二阶微分方程的特解；1求二阶微分方程式23的特解¢,;转子电流二阶微分方程对应的齐次方程的特征根为两个不相等的实数根，故可设二阶微分方程式23的特解4为：式中，bi为待求常数；对式26求一阶导和二阶导得：将式26—式28代入式21中，可得：对式29进行求解，可得常数h的值为bde;将所求的匕的值代入式26，即可求得二阶微分方程式23的特解匕的值为：2求二阶微分方程式24的特解；转子电流二阶微分方程对应的齐次方程的特征根为两个不相等的实数根，故可设二阶微分方程式24的特解4为：式中，b2为待求常数；对式31求一阶导和二阶导得：将式31—式33代入式24，可得：对式34进行求解，可得b2的值为：将式35代入式31所设的特解2中，可得二阶微分方程式24特解4的值为：3转子二阶微分方程的特解将所求的二阶微分方程式23的特解和式24的特解2代入式25中，即可求得转子二阶微分方程的特解为：⑶求取dq同步旋转坐标系下转子短路电流的解析式转子短路电流的解析式为转子二阶微分方程对应的齐次方程的通解irea®和转子二阶微分方程的特解之和，即为：6.根据权利要求5所述的计及RSC控制的双馈风机三相短路电流的解析方法，其特征在于:所述步骤五的具体步骤为：由定子磁链方程式⑶可得定子电流方程为：将定子磁链解析式式8和转子电流解析式式38代入定子电流方程式39，可得dq同步旋转坐标系下定子短路电流解析式为：7.根据权利要求6所述的计及RSC控制的双馈风机三相短路电流的解析方法，其特征在于:所述步骤六的具体步骤为：将所求得的dq旋转坐标系下DFIG定子和转子短路电流转换到三相静止坐标系下为：式中，ira、irb、ire分别为三相静止坐标系下转子a、b、C各相电流；isa、isb、isc分别为三相静止坐标系下定子A、B、C各相电流;Re为取实部符号；将式38代入式41可得计及转子侧变流器的双馈感应风力发电机的转子三相短路电流的解析式为：将式39代入式42可得计及转子侧变流器的双馈感应风力发电机的定子三相短路电流的解析式为：

百度查询： 昆明理工大学 一种计及RSC控制的双馈风机三相短路电流的解析方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD