一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法技术

本发明专利技术公开了一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，包括以下步骤，步骤1，计算各双馈风力发电机组的撬棒保护动作新判据，进而判断各机组撬棒保护是否动作，以动作与否为依据将风电场中的双馈机组进行一次分群；步骤2，计算各机组的分群特征量，并基于密度峰值聚类算法对一次分群结果中的群组进行二次分群，得到风电场的分群结果；步骤3，分别将同一群组的机组等值为一台机，并计算相应等值参数；步骤4，根据所得双馈风电机组风电场的等值模型，计算双馈风电场的短路电流。本发明专利技术更加精确地表达了风电场并网点短路电流的大小，计算精度高。

A New Criterion for Wind Farms Short Circuit Current Calculation Based on Doubly-fed Machine Prying Bar Protection

The invention discloses a method for calculating short-circuit current of wind farm based on a new criterion of pry bar protection of doubly-fed generators, which includes the following steps: step 1, calculating a new criterion for pry bar protection action of doubly-fed wind turbines, and then judging whether pry bar protection of each generating unit operates or not, dividing the doubly-fed generators in wind farm into groups once according to the action or not; step 2, calculating each generating unit; Based on the clustering characteristic, the cluster in the first clustering result is subdivided into two groups based on the peak density clustering algorithm to get the clustering result of the wind farm. Step 3, the same group of units is equivalent to one machine, and the corresponding equivalent parameters are calculated. Step 4, the short-circuit current of the doubly-fed wind farm is calculated according to the equivalent model of the doubly-fed wind farm. The invention can more accurately express the short-circuit current of the wind farm grid-connected point, and has high calculation accuracy.

【技术实现步骤摘要】
一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法
本专利技术涉及一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，属于风电场电气经济安全分析与设计领域。
技术介绍
大多数装设撬棒保护的双馈风力发电机作为现阶段风电场的主流机型，撬棒的动作情况将直接影响双馈电机的短路电流特性，同时不同型号、不同的参数的双馈电机内阻也不同，且风力发电机的组成型式和繁多的集电线路使得风电场拓扑结构更加复杂，如果忽略这几点，按现有方法计算并网点(PCC端)三相短路电流，可能得到过大的结果，造成不必要的经济浪费。
技术实现思路
为了解决上述技术问题，本专利技术提供了一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，更加精确地表达了风电场并网点短路电流的大小，计算精度高。为了达到上述目的，本专利技术所采用的技术方案是：本专利技术的一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，包括以下步骤，步骤1，计算各双馈风力发电机组的撬棒保护动作新判据(判据是专业术语，撬棒保护动作新判据是指在扰动发生后，判断撬棒保护是否动作的新的依据，有别于旧的依据)，进而判断各机组撬棒保护是否动作，以动作与否为依据将风电场中的双馈机组进行一次分群；步骤2，计算各机组的分群特征量，并基于密度峰值聚类算法对一次分群结果中的群组进行二次分群，得到风电场的分群结果；步骤3，分别将同一群组的机组等值为一台机，并计算相应等值参数；步骤4，根据所得双馈风电机组风电场的等值模型，计算双馈风电场的短路电流。上述双馈风力发电机组的撬棒保护动作新判据的计算方法为，其中，T表示周期，通常T＝1/f＝1/50＝0.02s，f表示工频，为50Hz；t0为短路故障发生时刻，ir为dq同步旋转坐标系下的转子电流空间矢量；为转子电流参考值，j指复数域中的虚部部分，Us0为故障前机端电压幅值；Ps,ref和Qs,ref分别为有功功率参考值和无功功率参考值，Kd为双馈电机机端电压跌落系数，ψsm为定子磁链幅值，Ls为定子绕组等效电感，Lm为励磁电感；为故障后转子感应电动势e直流分量的衰减时间常数；ωp＝ωs–ωr为转差角速度，ωr为转子角速度，ωs为同步转速，τs＝L's/Rs为定子回路衰减时间常数；Kd为机端电压跌落系数；为故障后转子感应电动势e的直流分量；s＝ωp/ωs为转差率，为转子暂态等效电感，Lr为转子电感，μ＝(Rr+kp+jωpL'r)/L'r和λ＝ki/L'r为参数变量，Rr为转子电阻，kp和ki分别为转子电流内环PI控制器的比例系数和积分时间常数，ωp为转差角速度；为常微分方程的特征根；为转子电流初值，和分别代表定子、转子电感耦合系数，为定子磁链初值，为转子磁链初值，Rs为定子电阻，Is0为故障后定子电流初值；Irth为转子电流允许极限值；比较临界动作电压跌落系数Kd的大小与机端电压幅值跌落量Kdcr大小，若满足Kdcr≥Kd，则该台机组撬棒保护动作；否则撬棒保护不动作。上述一次分群的结果为：撬棒保护动作的机组为一群，撬棒保护不动作的机组为一群。上述各机组的分群特征量包括故障前机端电压初始值us0、故障后机端电压稳态值usf、转子转速ωr、风力发电机的有功功率Ps0和无功功率Qs0。上述各机组的分群特征量具体的计算方法如下：在潮流计算中将双馈风力发电机视为PQ节点，按照牛顿拉夫逊法求得各机组故障前机端电压Us0；按照电力系统短路计算的方法，将双馈风力发电机等效为次暂态电势1和次暂态电抗x'，求解各机组故障后机端电压Usf；根据机型厂家提供手册中的标准风速-功率曲线，由风速的大小查出各风力发电机组对应的故障前有功功率Ps0大小；双馈风力发电机组均以恒功率因数方式运行，其无功功率Qs0根据功率因数和有功功率Ps0求解；双馈风力发电机组的转子转速ωr根据对应型号机组的技术手册中已知的转速-功率曲线求得。将同一群组的机组等值为一台机，其相应等值参数的计算方法为，若等值后风电场中共有h台等值机，其中等值机p是由等值前风电场中的m台机组组成，则将这m台机组的每一个的下标编号构成的有序数对集合记为Mij＝{(i,j)}；1)等值发电机参数所述等值发电机参数计算包括等值发电机的容量Seq、电阻Req、电抗Leq、电容Ceq和撬棒电阻Rcb_eq，应有以下关系成立：式中，Si为等值前第i台风力发电机的额定容量，Ri为等值前第i台风力发电机的电阻；Li为等值前第i台风力发电机的电抗；Ci为等值前第i台风力发电机的电容；Rcbi为等值前第i台风力发电机的撬棒电阻；2)等值变压器参数每一台风力发电机出口都经过一台小容量的升压变压器升压后接入系统，对于同群机组等值后的等值机出口同样需要一台等值升压变压器；变压器参数计算即相当于对等值变压器的容量ST_eq、阻抗ZT_eq进行计算，它们应该满足如下关系：式中，STi为等值前第i台变压器的容量；ZTi为等值前第i台变压器的阻抗；3)等值集电线路参数第p台等值机对应的等值集电线路阻抗公式为Slj为第l行第j列双馈电机的视在功率；Zkj为第k行第j列双馈电机所对应的集电线路阻抗；Sij为第i行第j列双馈电机的视在功率；Zeq_p为第p台等值机所对应的等值集电线路阻抗。上述双馈风电场的短路电流的计算方法如下：其中，对于撬棒保护不动作的等值机的A相短路电流计算公式为，其中，Ts和Tr分别为定子和转子回路暂态时间常数，Ts＝Ls’/Rs、Tr＝Lr’/Rr；Asa、Bsa、Csa分别为各频率分量系数，Rs为定子侧电阻，is为定子侧电流，ωs为同步角速度，ωp为转差角速度，Lm为激磁电感，ψs为定子侧磁链，Rr为转子侧电阻，ir为转子侧电流，ψr为定子侧磁链，Ls和Lr分别为定子自感和转子自感，Ls＝Lsσ+Lm，Lr＝Lrσ+Lm，Lsσ和Lrσ分别为定子漏感和转子漏感，Ls’和Lr’分别代表定子绕组暂态电感和转子绕组暂态电感，并且Ls’＝Ls-Lm2/Lr、Lr’＝Lr-Lm2/Ls，kr为转子电感耦合系数，kr＝Lm/Lr，ks为定子电感耦合系数，ks＝Lm/Ls；Kd为定子侧电压跌落系数；对于撬棒保护动作的等值机的A相短路电流计算公式为，其中，Tr”为转子回路暂态时间常数，Tr”＝L’r/(Rr+Rcb)；Asa2、Bsa2、Csa2分别为Crowbar保护动作下双馈电机短路电流各频率分量系数，R1为Crowbar保护等效阻抗，R＝Rr+Rcb，Rcb为Crowbar保护电阻。本专利技术的有益效果：本专利技术在计算风电场短路电流时，基于撬棒保护动作新判据和密度峰值聚类算法对风电场进行一次和二次分群，同时计及了集电线路转移阻抗和不同参数双馈风力发电机组成对风电场短路电流的影响，更加精确地表达了风电场并网点(PCC端)短路电流的大小，与现有计算公式相比，提高了计算的精度，对于风电场的电气设备经济安全设计具有重要意义。附图说明图1为本专利技术的工作流程图；图2为风电场系统接线图；图3为风电场的等值机组接线图。具体实施方式下面结合附图对本专利技术作进一步描述。以下实施例仅用于更加清楚地说明本专利技术的技术方案，而不能以此来限制本专利技术的保护范围。如图1所示，一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，包括以下步骤：步骤1，计算各双馈风力发电机组的撬棒保护动作本文档来自技高网...

【技术保护点】
1.一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，其特征在于，包括以下步骤，步骤1，计算各双馈风力发电机组的撬棒保护动作新判据，进而判断各机组撬棒保护是否动作，以动作与否为依据将风电场中的双馈机组进行一次分群；步骤2，计算各机组的分群特征量，并基于密度峰值聚类算法对一次分群结果中的群组进行二次分群，得到风电场的分群结果；步骤3，分别将同一群组的机组等值为一台机，并计算相应等值参数；步骤4，根据所得双馈风电机组风电场的等值模型，计算双馈风电场的短路电流。

【技术特征摘要】
1.一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，其特征在于，包括以下步骤，步骤1，计算各双馈风力发电机组的撬棒保护动作新判据，进而判断各机组撬棒保护是否动作，以动作与否为依据将风电场中的双馈机组进行一次分群；步骤2，计算各机组的分群特征量，并基于密度峰值聚类算法对一次分群结果中的群组进行二次分群，得到风电场的分群结果；步骤3，分别将同一群组的机组等值为一台机，并计算相应等值参数；步骤4，根据所得双馈风电机组风电场的等值模型，计算双馈风电场的短路电流。2.根据权利要求1所述的一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，其特征在于，所述双馈风力发电机组的撬棒保护动作新判据的计算方法为，其中，T表示周期，f表示工频；t0为短路故障发生时刻，ir为dq同步旋转坐标系下的转子电流空间矢量；为转子电流参考值，j指复数域中的虚部部分，Us0为故障前机端电压幅值；Ps,ref和Qs,ref分别为有功功率参考值和无功功率参考值，Kd为双馈电机机端电压跌落系数，ψsm为定子磁链幅值，Ls为定子绕组等效电感，Lm为励磁电感；为故障后转子感应电动势e直流分量的衰减时间常数；ωp＝ωs–ωr为转差角速度，ωr为转子角速度，ωs为同步转速，τs＝L's/Rs为定子回路衰减时间常数；Kd为机端电压跌落系数；为故障后转子感应电动势e的直流分量；s＝ωp/ωs为转差率，为转子暂态等效电感，Lr为转子电感，μ＝(Rr+kp+jωpL'r)/L'r和λ＝ki/L'r为参数变量，Rr为转子电阻，kp和ki分别为转子电流内环PI控制器的比例系数和积分时间常数，ωp为转差角速度；为常微分方程的特征根；为转子电流初值，和分别代表定子、转子电感耦合系数，为定子磁链初值，为转子磁链初值，Rs为定子电阻，Is0为故障后定子电流初值；Irth为转子电流允许极限值；比较临界动作电压跌落系数Kd的大小与机端电压幅值跌落量Kdcr大小，若满足Kdcr≥Kd，则该台机组撬棒保护动作；否则撬棒保护不动作。3.根据权利要求1所述的一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，其特征在于，所述一次分群的结果为：撬棒保护动作的机组为一群，撬棒保护不动作的机组为一群。4.根据权利要求1所述的一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，其特征在于，所述各机组的分群特征量包括故障前机端电压初始值us0、故障后机端电压稳态值usf、转子转速ωr、风力发电机的有功功率Ps0和无功功率Qs0。5.根据权利要求4所述的一种基于双馈电机撬棒保护新判据的风电场短路电流计算方法，其特征在于，所述各机组的分群特征量具体的计算方法如下：在潮流计算中将双馈风力发电机视为PQ节点，按照牛顿拉夫逊法求得各机组故障前机端电压Us0；按照电力系统短路计算的方法，将双馈风力发电机等效为次暂态电势1和次暂态电抗x'，求解各机组故障后机端电压Usf；根据机型厂家提供手册中的标准风速-功率曲线，由风速的大小查出各风力发电机组对应的故障前有功...

【专利技术属性】
技术研发人员：潘文霞，刘明洋，赵坤，张艺博，郭家圣，杨刚，李建成，陈昕，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32