CN109617121A - 一种针对次同步振荡的风电并网系统安全运行方法及系统 - Google Patents
一种针对次同步振荡的风电并网系统安全运行方法及系统 Download PDFInfo
- Publication number
- CN109617121A CN109617121A CN201811600781.2A CN201811600781A CN109617121A CN 109617121 A CN109617121 A CN 109617121A CN 201811600781 A CN201811600781 A CN 201811600781A CN 109617121 A CN109617121 A CN 109617121A
- Authority
- CN
- China
- Prior art keywords
- matrix
- model
- wind
- axis component
- output
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000010355 oscillation Effects 0.000 title claims abstract description 135
- 230000010354 integration Effects 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000012360 testing method Methods 0.000 claims abstract description 36
- 230000000644 propagated effect Effects 0.000 claims abstract description 10
- 238000012544 monitoring process Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 104
- 230000001360 synchronised effect Effects 0.000 claims description 42
- 238000012546 transfer Methods 0.000 claims description 40
- 230000004044 response Effects 0.000 claims description 34
- 238000004804 winding Methods 0.000 claims description 33
- 238000002347 injection Methods 0.000 claims description 25
- 239000007924 injection Substances 0.000 claims description 25
- 230000009466 transformation Effects 0.000 claims description 17
- 239000003990 capacitor Substances 0.000 claims description 16
- 238000010276 construction Methods 0.000 claims description 11
- 230000001052 transient effect Effects 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 8
- 230000003321 amplification Effects 0.000 claims description 7
- 230000005284 excitation Effects 0.000 claims description 7
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 7
- 239000003795 chemical substances by application Substances 0.000 claims description 6
- 230000010363 phase shift Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 description 6
- 238000010248 power generation Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 239000004576 sand Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Classifications
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/01—Arrangements for reducing harmonics or ripples
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P9/00—Arrangements for controlling electric generators for the purpose of obtaining a desired output
- H02P9/007—Control circuits for doubly fed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/002—Flicker reduction, e.g. compensation of flicker introduced by non-linear load
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开一种针对次同步振荡的风电并网系统安全运行方法及系统。该方法包括：获取对风电并网系统中的间谐波进行实时监测所得到的间谐波电流，得到间谐波电流向量；以间谐波电流向量为输入计算振荡传播因子；判断振荡传播因子是否大于预设阈值，得到判断结果；当判断结果表示大于时，确定次同步振荡传播到待测发电机，将引起次同步振荡的风机从风电并网系统中切除；当判断结果表示小于或等于时，确定次同步振荡仅在风机组范围内传播，调整风机组的运行参数抑制间谐波。本发明的针对次同步振荡的风电并网系统安全运行方法及系统，能够对次同步振荡进行监测并在次同步振荡发生时保证风电并网系统的安全运行。
Description
技术领域
本发明涉及电网安全运行领域，特别是涉及一种针对次同步振荡的风电并网系统安全运行方法及系统。
背景技术
风电，光伏等新能源电源已经成为电力系统中的重要组成部分。随着新能源比例的不断提高，由于新能源接入而导致的电力系统稳定性问题得到了广泛的关注。风能是新能源发电的最佳解决方案之一。目前双馈感应风力发电机和直驱永磁风机为风电场中的主流机型。为了提高风电并网电力系统的稳定性和可靠性，有必要对含风电电力系统中的次同步振荡现象进行更深入的研究。
双馈风机(DoublyFed Induction Generator,DFIG)转子侧和网侧变流器中含有大功率电力电子器件，在风速的随机变化引起双馈风力发电机转子励磁电流的频率变化时，转换或调节过程中由于开关效应，转子电流频率与电网频率共同作用会导致DFIG电流中存在次同步频率的间谐波。间谐波电流注入到电网后，会引起电压、电流和功率的振荡，引发各种故障以及电力系统的不稳定。次同步振荡在电力网络中的振荡幅值大，持续时间长，可能损害电力设备甚至危害电力系统的安全稳定运行。如何对次同步振荡进行监测并在次同步振荡发生时保证风电并网系统的安全运行成为目前亟待解决的问题。
发明内容
本发明的目的是提供一种针对次同步振荡的风电并网系统安全运行方法及系统，对次同步振荡进行监测并在次同步振荡发生时保证风电并网系统的安全运行。
为实现上述目的，本发明提供了如下方案：
一种针对次同步振荡的风电并网系统安全运行方法，包括：
获取对风电并网系统中的间谐波进行实时监测所得到的间谐波电流，得到间谐波电流向量；
利用公式计算振荡传播因子；ki(s)为待测发电机的输出向量中第i个元素所对应的振荡传播因子；|Δyi|为待测发电机的输出向量中第i个元素的幅值，Δuj为间谐波电流向量的第j个元素；||Δu||2为间谐波电流向量的二范数；Wij(s)为间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
判断所述振荡传播因子是否大于预设阈值，得到判断结果；
当所述判断结果表示大于时，确定次同步振荡传播到所述待测发电机，将引起次同步振荡的风机从所述风电并网系统中切除；
当所述判断结果表示小于或等于时，确定次同步振荡仅在风机组范围内传播，调整风机组的运行参数抑制所述间谐波。
可选的，所述间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应的确定方法为：
建立双馈风机的动态模型；
将所述双馈风机的动态模型线性化，得到风电系统小干扰模型；
对所述风电系统小干扰模型进行拉氏变换，并设定初始条件为零，得到待测发电机的输出向量与间谐波电流向量之间的矩阵传递函数，其中Wij(s)为所述矩阵传递函数中的元素。
可选的，所述建立双馈风机的动态模型，具体包括：
构建异步发电机的动态模型
其中，E′d为定子暂态电势的d轴分量，Udr为转子绕组电压的d轴分量，E′q为定子暂态电势的q轴分量、Uqr为转子绕组电压的q轴分量；Ids和Iqs分别为定子绕组电流的d轴、q轴分量，s为转差率，Rr为转子绕组电阻，Xm为激磁电抗，Xrr为转子电抗，Xrr＝Xm+Xr，其中Xr为转子漏抗，TJ为电机惯性时间常数，Pe为电磁功率，Pm为机械功率，D为阻尼系数，s0为稳态时的转差率，ω0为同步角频率；
考虑直流电压动态特性，构建中间电容器的动态模型
其中，Pr为转子侧变流器输出到中间电容器的有功功率，Pr3为网侧变流器输出到交流侧的有功功率；C为中间电容器的电容值，Udc为电容器的直流电压；
构建滤波电抗的模型
其中，Idr3为网侧变流器交流侧电流的d轴分量，Ud1为网侧变流器交流侧电压的d轴分量，Iqr3为网侧变流器交流侧电流的q轴分量，Uq1为网侧变流器交流侧电压的q轴分量，Xr3为滤波电抗值；Uds和Uqs分别为定子电压d轴和q轴分量；
构建转子侧变流器的动态模型
其中，x1、x2、x3和x4为引入的状态变量，Ki1、Ki2、Ki3和Ki4为相应PI控制器的积分系数，Ps和Qs分别为定子输出的有功功率和无功功率，Idr和Iqr分别为转子绕组电流的d轴、q轴分量；Ps *和Qs *分别为定子输出的有功功率参考值和无功功率参考值，和分别为转子绕组电流的d轴分量参考值和q轴分量参考值；
构建网侧变流器的动态模型
其中，x5、x6和x7为引入的状态变量，Ki5、Ki6和Ki7相应PI控制器的积分系数；为电容器的直流电压参考值，为网侧变流器交流侧电流的d轴分量参考值，为网侧变流器交流侧电流的q轴分量参考值。
可选的，所述将所述双馈风机的动态模型线性化，得到风电系统小干扰模型，具体包括：
对所述异步发电机的动态模型、所述中间电容器的动态模型、所述滤波电抗的模型、所述转子侧变流器的动态模型和网侧变流器的动态模型进行线性化，得到双馈风机的小干扰模型
ΔIW＝CWΔXW+DWΔUW
其中ΔXW＝[ΔE′d,ΔE′q,Δs,ΔUdc,ΔIdr3,ΔIqr3,Δx1,Δx2,Δx3,Δx4,Δx5,Δx6,Δx7]T为双馈风机状态变量，ΔUW＝[ΔUd,ΔUq]T为双馈风机接入节点电压，ΔIW＝[ΔId,ΔIq]T为双馈风机接入节点注入电流；ΔUd和ΔUq分别为双馈风机接入母线电压的d轴和q轴分量；ΔId和ΔIq分别为双馈风机接入母线注入电流的d轴和q轴分量；AW、BW、CW和DW分别为双馈风机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
线性化后的同步发电机状态方程：
ΔIG＝CGΔXG+DGΔUG
其中ΔXG为同步发电机状态变量，ΔUG为同步发电机接入节点电压，ΔIG为同步发电机接入节点注入电流；AG、BG、CG和DG分别为同步发电机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
考虑负荷的网络方程：
ΔI＝YΔU
ΔI为节点注入电流，Y为考虑负荷的节点导纳矩阵，ΔU为节点电压；
将双馈风机的小干扰模型、线性化后的同步发电机状态方程以及考虑负荷的网络方程联立，得到整个系统的小干扰状态空间模型
Δy＝CΔX+DΔu
其中，ΔX＝[ΔX1,ΔX2…ΔXn]T为n维状态向量，Δu＝[Δu1,Δu2…Δur]T为r维输入向量，Δy＝[Δy1,Δy2…Δym]为m维输出向量，A、B、C和D分别为n×n阶状态矩阵、n×r阶输入矩阵、m×n阶输出矩阵和m×r阶直接传递矩阵。
可选的，所述对所述风电系统小干扰模型进行拉氏变换，并设定初始条件为零，得到待测发电机的输出向量与间谐波电流向量之间的矩阵传递函数，具体包括：
对整个系统的小干扰状态空间模型进行拉氏变换，并设定初始条件为零，得到
ΔX(s)＝(sI-A)-1BΔu(s)
Δy(s)＝CΔX(s)+DΔu(s)
则输出向量Δy与输入向量Δu之间的矩阵传递函数为：
W(s)＝C(sI-A)-1B+D
W(s)为m×r阶矩阵传递函数，矩阵元素Wij(s)表示由输入变量第j个元素Δuj到输出向量第i个元素Δyi的正弦响应；对于角频率为ω的输入变量第j个元素Δuj，输出向量第i个元素Δyi的幅值放大倍数为|Wij(jω)|，相移为∠Wij(jω)，Δyi的总响应等于各个输入响应的线性和，即：
本发明还公开一种针对次同步振荡的风电并网系统安全运行系统，包括：
谐波获取模块，用于获取对风电并网系统中的间谐波进行实时监测所得到的间谐波电流，得到间谐波电流向量；
振荡传播因子计算模块，用于利用公式计算振荡传播因子；ki(s)为待测发电机的输出向量中第i个元素所对应的振荡传播因子；|Δyi|为待测发电机的输出向量中第i个元素的幅值，Δuj为间谐波电流向量的第j个元素；||Δu||2为间谐波电流向量的二范数；Wij(s)为间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
判断模块，用于判断所述振荡传播因子是否大于预设阈值，得到判断结果；
切除模块，用于当所述判断结果表示大于时，确定次同步振荡传播到所述待测发电机，将引起次同步振荡的风机从所述风电并网系统中切除；
运行参数调整模块，用于当所述判断结果表示小于或等于时，确定次同步振荡仅在风机组范围内传播，调整风机组的运行参数抑制所述间谐波。
可选的，该系统还包括正弦响应确定模块，用于确定所述间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
所述正弦响应确定模块包括：
动态模型建立单元，用于建立双馈风机的动态模型；
模型线性化单元，用于将所述双馈风机的动态模型线性化，得到风电系统小干扰模型；
矩阵传递函数确定单元，用于对所述风电系统小干扰模型进行拉氏变换，并设定初始条件为零，得到待测发电机的输出向量与间谐波电流向量之间的矩阵传递函数，其中Wij(s)为所述矩阵传递函数中的元素。
可选的，所述动态模型建立单元，具体包括：
异步发电机模型构建子单元，用于构建异步发电机的动态模型
其中，E′d为定子暂态电势的d轴分量，Udr为转子绕组电压的d轴分量，E′q为定子暂态电势的q轴分量、Uqr为转子绕组电压的q轴分量；Ids和Iqs分别为定子绕组电流的d轴、q轴分量，s为转差率，Rr为转子绕组电阻，Xm为激磁电抗，Xrr＝Xm+Xr，其中Xr为转子漏抗，TJ为电机惯性时间常数，Pe为电磁功率，Pm为机械功率，D为阻尼系数，s0为稳态时的转差率，ω0为同步角频率；
电容模型构建子单元，用于考虑直流电压动态特性，构建中间电容器的动态模型
其中，Pr为转子侧变流器输出到中间电容器的有功功率，Pr3为网侧变流器输出到交流侧的有功功率；C为中间电容器的电容值，Udc为电容器的直流电压；
电抗模块构建子单元，用于构建滤波电抗的模型
其中，Idr3为网侧变流器交流侧电流的d轴分量，Ud1为网侧变流器交流侧电压的d轴分量，Iqr3为网侧变流器交流侧电流的q轴分量，Uq1为网侧变流器交流侧电压的q轴分量，Xr3为滤波电抗值；Uds和Uqs分别为定子电压d轴和q轴分量；
转子侧变流器模型构建子单元，用于构建转子侧变流器的动态模型
其中，x1、x2、x3和x4为引入的状态变量，Ki1、Ki2、Ki3和Ki4为相应PI控制器的积分系数，Ps和Qs分别为定子输出的有功功率和无功功率，Idr和Iqr分别为转子绕组电流的d轴、q轴分量；Ps *和Qs *分别为定子输出的有功功率参考值和无功功率参考值，和分别为转子绕组电流的d轴分量参考值和q轴分量参考值；
网侧变流器模型构建子单元，用于构建网侧变流器的动态模型
其中，x5、x6和x7为引入的状态变量，Ki5、Ki6和Ki7相应PI控制器的积分系数；为电容器的直流电压参考值，为网侧变流器交流侧电流的d轴分量参考值，为网侧变流器交流侧电流的q轴分量参考值。
可选的，所述模型线性化单元，具体包括：
线性化子单元，用于对所述异步发电机的动态模型、所述中间电容器的动态模型、所述滤波电抗的模型、所述转子侧变流器的动态模型和网侧变流器的动态模型进行线性化，得到双馈风机的小干扰模型
ΔIW＝CWΔXW+DWΔUW
其中ΔXW＝[ΔE′d,ΔE′q,Δs,ΔUdc,ΔIdr3,ΔIqr3,Δx1,Δx2,Δx3,Δx4,Δx5,Δx6,Δx7]T为双馈风机状态变量，ΔUW＝[ΔUd,ΔUq]T为双馈风机接入节点电压，ΔIW＝[ΔId,ΔIq]T为双馈风机接入节点注入电流；ΔUd和ΔUq分别为双馈风机接入母线电压的d轴和q轴分量；ΔId和ΔIq分别为双馈风机接入母线注入电流的d轴和q轴分量；AW、BW、CW和DW分别为双馈风机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
线性化后的同步发电机状态方程：
ΔIG＝CGΔXG+DGΔUG
其中ΔXG为同步发电机状态变量，ΔUG为同步发电机接入节点电压，ΔIG为同步发电机接入节点注入电流；AG、BG、CG和DG分别为同步发电机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
考虑负荷的网络方程：
ΔI＝YΔU
ΔI为节点注入电流，Y为考虑负荷的节点导纳矩阵，ΔU为节点电压；
多方程联立子单元，用于将双馈风机的小干扰模型、线性化后的同步发电机状态方程以及考虑负荷的网络方程联立，得到整个系统的小干扰状态空间模型
Δy＝CΔX+DΔu
其中，ΔX＝[ΔX1,ΔX2…ΔXn]T为n维状态向量，Δu＝[Δu1,Δu2…Δur]T为r维输入向量，Δy＝[Δy1,Δy2…Δym]为m维输出向量，A、B、C和D分别为n×n阶状态矩阵、n×r阶输入矩阵、m×n阶输出矩阵和m×r阶直接传递矩阵。
可选的，所述矩阵传递函数确定单元，具体包括：
拉氏变换子单元，用于对整个系统的小干扰状态空间模型进行拉氏变换，并设定初始条件为零，得到
ΔX(s)＝(sI-A)-1BΔu(s)
Δy(s)＝CΔX(s)+DΔu(s)
则输出向量Δy与输入向量Δu之间的矩阵传递函数为：
W(s)＝C(sI-A)-1B+D
W(s)为m×r阶矩阵传递函数，矩阵元素Wij(s)表示由输入变量第j个元素Δuj到输出向量第i个元素Δyi的正弦响应；对于角频率为ω的输入变量第j个元素Δuj，输出向量第i个元素Δyi的幅值放大倍数为|Wij(jω)|，相移为∠Wij(jω)，Δyi的总响应等于各个输入响应的线性和，即：
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明所公开的针对次同步振荡的风电并网系统安全运行方法及系统，获取间谐波电流向量，并利用间谐波电流向量计算振荡传播因子，从而确定次同步振荡的传播范围，当次同步振荡传播到待测发电机时切除引起次同步振荡的风机，当次同步振荡未传播到待测发电机时只需要风机组的运行参数即可。本发明通过对振荡传播因子进行计算，实现了对次同步振荡进行监测并在次同步振荡发生时保证风电并网系统的安全运行。同时，当次同步振荡不可避免时切除对应风机，在可避免时调整运行参数，最大程度上在保证安全运行的基础上减少运行调控成本。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明针对次同步振荡的风电并网系统安全运行方法实施例1的方法流程图；
图2为双馈风机的结构图；
图3为实施例2中四机两区系统的系统结构图；
图4为各发电机的振荡传播因子曲线图；
图5为系统模式1-4固有振荡频率附近不同发电机的振荡传播因子；
图6为DFIG2振荡传播因子的曲线图；
图7为DFIG2有功功率仿真结果图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种针对次同步振荡的风电并网系统安全运行方法及系统，对次同步振荡进行监测并在次同步振荡发生时保证风电并网系统的安全运行。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1：
图1为本发明针对次同步振荡的风电并网系统安全运行方法实施例1的方法流程图。
参见图1，该针对次同步振荡的风电并网系统安全运行方法，包括：
步骤101：获取间谐波电流向量：获取对风电并网系统中的间谐波进行实时监测所得到的间谐波电流，得到间谐波电流向量；
步骤102：计算振荡传播因子：利用公式计算振荡传播因子；ki(s)为待测发电机的输出向量中第i个元素所对应的振荡传播因子；|Δyi|为待测发电机的输出向量中第i个元素的幅值，Δuj为间谐波电流向量的第j个元素；||Δu||2为间谐波电流向量的二范数；Wij(s)为间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
步骤103：阈值判断：判断所述振荡传播因子是否大于预设阈值，得到判断结果；
步骤104：风机切除：当所述判断结果表示大于时，确定次同步振荡传播到所述待测发电机，将引起次同步振荡的风机从所述风电并网系统中切除；
步骤105：运行参数调整：当所述判断结果表示小于或等于时，确定次同步振荡仅在风机组范围内传播，调整风机组的运行参数抑制所述间谐波。
所述间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应的确定方法包括：
(1)建立双馈风机的动态模型；
图2为双馈风机的结构图。
参见图2，双馈风机系统主要包括异步发电机及其传动系统，转子侧变流器及其控制系统和网侧变流器及其控制系统。图2中Us为双馈风机的定子电压，Is为定子电流；Ur为转子电压，Ir为转子电流；U1为网侧变流器交流侧电压，Ir3为网侧变流器交流侧电流；Lr3为网侧变流器交流侧与定子连接线路的滤波电抗。C为中间电容器的电容值，Udc为电容器的直流电压。物理量正方向采用发电机惯例。
忽略电机定子磁链暂态过程，构建异步发电机的动态模型
其中，E′d为定子暂态电势的d轴分量，Udr为转子绕组电压的d轴分量，E′q为定子暂态电势的q轴分量、Uqr为转子绕组电压的q轴分量；Ids和Iqs分别为定子绕组电流的d轴、q轴分量，s为转差率，Rr为转子绕组电阻，Xm为激磁电抗，Xrr为转子电抗，Xrr＝Xm+Xr，其中Xr为转子漏抗，TJ为电机惯性时间常数，Pe为电磁功率，Pm为机械功率，D为阻尼系数，s0为稳态时的转差率，ω0为同步角频率；
考虑直流电压动态特性，构建中间电容器的动态模型
其中，Pr为转子侧变流器输出到中间电容器的有功功率，Pr3为网侧变流器输出到交流侧的有功功率；C为中间电容器的电容值，Udc为电容器的直流电压；
构建滤波电抗的模型
其中，Idr3为网侧变流器交流侧电流的d轴分量，Ud1为网侧变流器交流侧电压的d轴分量，Iqr3为网侧变流器交流侧电流的q轴分量，Uq1为网侧变流器交流侧电压的q轴分量，Xr3为滤波电抗值；Uds和Uqs分别为定子电压d轴和q轴分量；
转子侧变流器的控制目标是通过控制励磁电压，使定子侧的有功功率和无功功率实现解耦独立控制。构建转子侧变流器的动态模型
其中，x1、x2、x3和x4为引入的状态变量，Ki1、Ki2、Ki3和Ki4为相应PI控制器的积分系数，Ps和Qs分别为定子输出的有功功率和无功功率，Idr和Iqr分别为转子绕组电流的d轴、q轴分量；Ps *和Qs *分别为定子输出的有功功率参考值和无功功率参考值，和分别为转子绕组电流的d轴分量参考值和q轴分量参考值；
网侧变流器的控制目标是维持中间电容电压稳定，并控制转子侧与电网的无功功率交换为0。构建网侧变流器的动态模型
其中，x5、x6和x7为引入的状态变量，Ki5、Ki6和Ki7相应PI控制器的积分系数；为电容器的直流电压参考值，为网侧变流器交流侧电流的d轴分量参考值，为网侧变流器交流侧电流的q轴分量参考值。
(2)将所述双馈风机的动态模型线性化，得到风电系统小干扰模型：
对所述异步发电机的动态模型、所述中间电容器的动态模型、所述滤波电抗的模型、所述转子侧变流器的动态模型和网侧变流器的动态模型进行线性化，得到双馈风机的小干扰模型
ΔIW＝CWΔXW+DWΔUW
其中ΔXW＝[ΔE′d,ΔE′q,Δs,ΔUdc,ΔIdr3,ΔIqr3,Δx1,Δx2,Δx3,Δx4,Δx5,Δx6,Δx7]T为双馈风机状态变量，ΔUW＝[ΔUd,ΔUq]T为双馈风机接入节点电压，ΔIW＝[ΔId,ΔIq]T为双馈风机接入节点注入电流；ΔUd和ΔUq分别为双馈风机接入母线电压的d轴和q轴分量；ΔId和ΔIq分别为双馈风机接入母线注入电流的d轴和q轴分量；AW、BW、CW和DW分别为双馈风机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
线性化后的同步发电机状态方程：
ΔIG＝CGΔXG+DGΔUG
其中ΔXG为同步发电机状态变量，ΔUG为同步发电机接入节点电压，ΔIG为同步发电机接入节点注入电流；AG、BG、CG和DG分别为同步发电机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
考虑负荷的网络方程：
ΔI＝YΔU
ΔI为节点注入电流，Y为考虑负荷的节点导纳矩阵，ΔU为节点电压；
含双馈风机系统中的同步发电机采用经典模型，负荷采用恒定阻抗模型。将双馈风机的小干扰模型、线性化后的同步发电机状态方程以及考虑负荷的网络方程联立，得到整个系统的小干扰状态空间模型
Δy＝CΔX+DΔu
其中，ΔX＝[ΔX1,ΔX2…ΔXn]T为n维状态向量，Δu＝[Δu1,Δu2…Δur]T为r维输入向量，Δy＝[Δy1,Δy2…Δym]为m维输出向量，A、B、C和D分别为n×n阶状态矩阵、n×r阶输入矩阵、m×n阶输出矩阵和m×r阶直接传递矩阵。
(3)对所述风电系统小干扰模型进行拉氏变换，并设定初始条件为零，得到待测发电机的输出向量与间谐波电流向量之间的矩阵传递函数，其中Wij(s)为所述矩阵传递函数中的元素。
对整个系统的小干扰状态空间模型进行拉氏变换，并设定初始条件为零，得到
ΔX(s)＝(sI-A)-1BΔu(s)
Δy(s)＝CΔX(s)+DΔu(s)
则输出向量Δy与输入向量Δu之间的矩阵传递函数为：
W(s)＝C(sI-A)-1B+D
W(s)为m×r阶矩阵传递函数，矩阵元素Wij(s)表示由输入变量第j个元素Δuj到输出向量第i个元素Δyi的正弦响应；对于角频率为ω的输入变量第j个元素Δuj，输出向量第i个元素Δyi的幅值放大倍数为|Wij(jω)|，相移为∠Wij(jω)，Δyi的总响应等于各个输入响应的线性和，即：
双馈风机间谐波注入电网后，造成的扰动会在电力系统中的不同点得到响应。在多输入多输出系统中，以间谐波电流向量为输出量，利用振荡传播因子来表征不同频率的输入变量扰动经过系统传播后对输出向量元素振荡幅值造成的影响，待测发电机的输出向量中第i个元素所对应的振荡传播因子由下式计算：
式中|Δyi|为待测发电机的输出向量中第i个元素的幅值，||Δu||2为间谐波电流向量的二范数。振荡传播因子越大，表示系统扰动在该点处影响越强，输出向量元素振荡幅值越大。利用振荡传播因子可以对次同步振荡传播范围进行确定。当振荡传播因子大于预设阈值(一般为-3dB)时，认为次同步间谐波传播到输出向量引对应的发电机引起次同步振荡，输出向量所代表的发电机(待测发电机)包括在传播范围内。当振荡传播因子小于预设阈值时，次同步间谐波没有传播到输出向量，输出向量所代表的发电机不包括在传播范围内。
该实施例中的输出向量可以为输出有功功率。
实施例2：
该实施例主要提供算例仿真与分析。
图3为实施例2中四机两区系统的系统结构图。
参见图3，在四机两区系统中接入两台DFIG，对图3系统进行小干扰稳定分析，得到7个固有振荡模式，表1为系统固有振荡模式表。参见表1，系统特征根实部均为负数，系统小干扰稳定。
表1
频域分析：
假设DFIG1受到扰动，引发间谐波电流注入电网，以该间谐波为输入变量，以各台双馈风机和同步发电机的输出有功功率为输出向量，建立系统传递函数，求其振荡传播因子。
图4为各发电机的振荡传播因子曲线图。
图5为系统模式1-4固有振荡频率附近不同发电机的振荡传播因子。
由图4和图5可知：对于不同角频率的输入变量，同一台发电机输出功率的振荡传播因子不同。同一角频率的输入变量下，各台发电机输出功率的振荡传播因子也不同。当输入变量角频率接近系统模式1-4的固有振荡角频率时，各台发电机的振荡传播因子均超过0dB，DFIG注入系统的间谐波经过系统传播起到了放大作用。输入变量频率接近模式3、4的固有频率时，双馈风机输出功率振荡幅值比同步发电机输出功率振荡幅值大，输入变量频率接近模式1、2的固有频率时，双馈风机输出功率振荡幅值比同步发电机输出功率振荡幅值小。
图6为DFIG2振荡传播因子的曲线图。
参见图6，DFIG2振荡传播因子的3个极大值点横坐标分别为90.79rad/s、89.61rad/s和75.12rad/s。这三个角频率分别与表1中模式1、2和3的系统固有振荡角频率几乎相等。当双馈风机注入电力网络的间谐波中出现该频率段的激发源时就会引发共振，出现幅值较大的功率振荡。
时域分析：
为了验证传递函数分析的准确性，在MATLAB中建立图3系统的时域仿真模型，并在母线6处注入间谐波电流，间谐波电流角频率为89.64rad/s。
图7为DFIG2有功功率仿真结果图。
参见图7，原系统母线6处的间谐波电流引起DFIG2有功功率振荡，振荡幅值为0.0123，振荡传播因子为8.70，与图6中89.61rad/s处的振荡传播因子1018.75/20＝8.66基本相同。
传播范围分析
根据图4和图5中求出的振荡传播因子，可以确定由间谐波引起的次同步振荡传播范围。例如，当DFIG1间谐波角频率等于72rad/s时，只有DFIG1的振荡传播因子大于-3dB，其余发电机的振荡传播因子均小于-3dB。该间谐波引起的次同步振荡只在DFIG1内传播，次同步振荡的传播范围内只有DFIG1，表现为风电机组的振荡。当DFIG1间谐波角频率等于89.64rad/s时，所有发电机的振荡传播因子均大于-3dB。该间谐波引起的次同步振荡在2台DFIG和4台同步发电机之间传播，次同步振荡的传播范围为整个系统，包括全部的发电机，表现为风电机组与火电机组的振荡。
实施例3：
本发明还公开一种针对次同步振荡的风电并网系统安全运行系统，包括：
谐波获取模块，用于获取对风电并网系统中的间谐波进行实时监测所得到的间谐波电流，得到间谐波电流向量；
振荡传播因子计算模块，用于利用公式计算振荡传播因子；ki(s)为待测发电机的输出向量中第i个元素所对应的振荡传播因子；|Δyi|为待测发电机的输出向量中第i个元素的幅值，Δuj为间谐波电流向量的第j个元素；||Δu||2为间谐波电流向量的二范数；Wij(s)为间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
判断模块，用于判断所述振荡传播因子是否大于预设阈值，得到判断结果；
切除模块，用于当所述判断结果表示大于时，确定次同步振荡传播到所述待测发电机，将引起次同步振荡的风机从所述风电并网系统中切除；
运行参数调整模块，用于当所述判断结果表示小于或等于时，确定次同步振荡仅在风机组范围内传播，调整风机组的运行参数抑制所述间谐波。
可选的，该系统还包括正弦响应确定模块，用于确定所述间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
所述正弦响应确定模块包括：
动态模型建立单元，用于建立双馈风机的动态模型；
模型线性化单元，用于将所述双馈风机的动态模型线性化，得到风电系统小干扰模型；
矩阵传递函数确定单元，用于对所述风电系统小干扰模型进行拉氏变换，并设定初始条件为零，得到待测发电机的输出向量与间谐波电流向量之间的矩阵传递函数，其中Wij(s)为所述矩阵传递函数中的元素。
可选的，所述动态模型建立单元，具体包括：
异步发电机模型构建子单元，用于构建异步发电机的动态模型
其中，E′d为定子暂态电势的d轴分量，Udr为转子绕组电压的d轴分量，E′q为定子暂态电势的q轴分量、Uqr为转子绕组电压的q轴分量；Ids和Iqs分别为定子绕组电流的d轴、q轴分量，s为转差率，Rr为转子绕组电阻，Xm为激磁电抗，Xrr为转子电抗，Xrr＝Xm+Xr，其中Xr为转子漏抗，TJ为电机惯性时间常数，Pe为电磁功率，Pm为机械功率，D为阻尼系数，s0为稳态时的转差率，ω0为同步角频率；
电容模型构建子单元，用于考虑直流电压动态特性，构建中间电容器的动态模型
其中，Pr为转子侧变流器输出到中间电容器的有功功率，Pr3为网侧变流器输出到交流侧的有功功率；C为中间电容器的电容值，Udc为电容器的直流电压；
电抗模块构建子单元，用于构建滤波电抗的模型
其中，Idr3为网侧变流器交流侧电流的d轴分量，Ud1为网侧变流器交流侧电压的d轴分量，Iqr3为网侧变流器交流侧电流的q轴分量，Uq1为网侧变流器交流侧电压的q轴分量，Xr3为滤波电抗值；Uds和Uqs分别为定子电压d轴和q轴分量；
转子侧变流器模型构建子单元，用于构建转子侧变流器的动态模型
其中，x1、x2、x3和x4为引入的状态变量，Ki1、Ki2、Ki3和Ki4为相应PI控制器的积分系数，Ps和Qs分别为定子输出的有功功率和无功功率，Idr和Iqr分别为转子绕组电流的d轴、q轴分量；Ps *和Qs *分别为定子输出的有功功率参考值和无功功率参考值，和分别为转子绕组电流的d轴分量参考值和q轴分量参考值；
网侧变流器模型构建子单元，用于构建网侧变流器的动态模型
其中，x5、x6和x7为引入的状态变量，Ki5、Ki6和Ki7相应PI控制器的积分系数；为电容器的直流电压参考值，为网侧变流器交流侧电流的d轴分量参考值，为网侧变流器交流侧电流的q轴分量参考值。
可选的，所述模型线性化单元，具体包括：
线性化子单元，用于对所述异步发电机的动态模型、所述中间电容器的动态模型、所述滤波电抗的模型、所述转子侧变流器的动态模型和网侧变流器的动态模型进行线性化，得到双馈风机的小干扰模型
ΔIW＝CWΔXW+DWΔUW
其中ΔXW＝[ΔE′d,ΔE′q,Δs,ΔUdc,ΔIdr3,ΔIqr3,Δx1,Δx2,Δx3,Δx4,Δx5,Δx6,Δx7]T为双馈风机状态变量，ΔUW＝[ΔUd,ΔUq]T为双馈风机接入节点电压，ΔIW＝[ΔId,ΔIq]T为双馈风机接入节点注入电流；ΔUd和ΔUq分别为双馈风机接入母线电压的d轴和q轴分量；ΔId和ΔIq分别为双馈风机接入母线注入电流的d轴和q轴分量；AW、BW、CW和DW分别为双馈风机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
线性化后的同步发电机状态方程：
ΔIG＝CGΔXG+DGΔUG
其中ΔXG为同步发电机状态变量，ΔUG为同步发电机接入节点电压，ΔIG为同步发电机接入节点注入电流；AG、BG、CG和DG分别为同步发电机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
考虑负荷的网络方程：
ΔI＝YΔU
ΔI为节点注入电流，Y为考虑负荷的节点导纳矩阵，ΔU为节点电压；
多方程联立子单元，用于将双馈风机的小干扰模型、线性化后的同步发电机状态方程以及考虑负荷的网络方程联立，得到整个系统的小干扰状态空间模型
Δy＝CΔX+DΔu
其中，ΔX＝[ΔX1,ΔX2…ΔXn]T为n维状态向量，Δu＝[Δu1,Δu2…Δur]T为r维输入向量，Δy＝[Δy1,Δy2…Δym]为m维输出向量，A、B、C和D分别为n×n阶状态矩阵、n×r阶输入矩阵、m×n阶输出矩阵和m×r阶直接传递矩阵。
可选的，所述矩阵传递函数确定单元，具体包括：
拉氏变换子单元，用于对整个系统的小干扰状态空间模型进行拉氏变换，并设定初始条件为零，得到
ΔX(s)＝(sI-A)-1BΔu(s)
Δy(s)＝CΔX(s)+DΔu(s)
则输出向量Δy与输入向量Δu之间的矩阵传递函数为：
W(s)＝C(sI-A)-1B+D
W(s)为m×r阶矩阵传递函数，矩阵元素Wij(s)表示由输入变量第j个元素Δuj到输出向量第i个元素Δyi的正弦响应；对于角频率为ω的输入变量第j个元素Δuj，输出向量第i个元素Δyi的幅值放大倍数为|Wij(jω)|，相移为∠Wij(jω)，Δyi的总响应等于各个输入响应的线性和，即：
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明所公开的针对次同步振荡的风电并网系统安全运行方法及系统，获取间谐波电流向量，并利用间谐波电流向量计算振荡传播因子，从而确定次同步振荡的传播范围，当次同步振荡传播到待测发电机时切除引起次同步振荡的风机，当次同步振荡未传播到待测发电机时只需要风机组的运行参数即可。本发明通过对振荡传播因子进行计算，实现了对次同步振荡进行监测并在次同步振荡发生时保证风电并网系统的安全运行。同时，当次同步振荡不可避免时切除对应风机，在可避免时调整运行参数，最大程度上在保证安全运行的基础上减少运行调控成本。
对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种针对次同步振荡的风电并网系统安全运行方法，其特征在于，包括：
获取对风电并网系统中的间谐波进行实时监测所得到的间谐波电流，得到间谐波电流向量；
利用公式计算振荡传播因子；ki(s)为待测发电机的输出向量中第i个元素所对应的振荡传播因子；|Δyi|为待测发电机的输出向量中第i个元素的幅值，Δuj为间谐波电流向量的第j个元素；||Δu||2为间谐波电流向量的二范数；Wij(s)为间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
判断所述振荡传播因子是否大于预设阈值，得到判断结果；
当所述判断结果表示大于时，确定次同步振荡传播到所述待测发电机，将引起次同步振荡的风机从所述风电并网系统中切除；
当所述判断结果表示小于或等于时，确定次同步振荡仅在风机组范围内传播，调整风机组的运行参数抑制所述间谐波。
2.根据权利要求1所述的针对次同步振荡的风电并网系统安全运行方法，其特征在于，所述间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应的确定方法为：
建立双馈风机的动态模型；
将所述双馈风机的动态模型线性化，得到风电系统小干扰模型；
对所述风电系统小干扰模型进行拉氏变换，并设定初始条件为零，得到待测发电机的输出向量与间谐波电流向量之间的矩阵传递函数，其中Wij(s)为所述矩阵传递函数中的元素。
3.根据权利要求2所述的针对次同步振荡的风电并网系统安全运行方法，其特征在于，所述建立双馈风机的动态模型，具体包括：
构建异步发电机的动态模型
其中，E′d为定子暂态电势的d轴分量，Udr为转子绕组电压的d轴分量，E′q为定子暂态电势的q轴分量、Uqr为转子绕组电压的q轴分量；Ids和Iqs分别为定子绕组电流的d轴、q轴分量，s为转差率，Rr为转子绕组电阻，Xm为激磁电抗，Xrr为转子电抗，Xrr＝Xm+Xr，其中Xr为转子漏抗，TJ为电机惯性时间常数，Pe为电磁功率，Pm为机械功率，D为阻尼系数，s0为稳态时的转差率，ω0为同步角频率；
考虑直流电压动态特性，构建中间电容器的动态模型
其中，Pr为转子侧变流器输出到中间电容器的有功功率，Pr3为网侧变流器输出到交流侧的有功功率；C为中间电容器的电容值，Udc为电容器的直流电压；
构建滤波电抗的模型
其中，Idr3为网侧变流器交流侧电流的d轴分量，Ud1为网侧变流器交流侧电压的d轴分量，Iqr3为网侧变流器交流侧电流的q轴分量，Uq1为网侧变流器交流侧电压的q轴分量，Xr3为滤波电抗值；Uds和Uqs分别为定子电压d轴和q轴分量；
构建转子侧变流器的动态模型
其中，x1、x2、x3和x4为引入的状态变量，Ki1、Ki2、Ki3和Ki4为相应PI控制器的积分系数，Ps和Qs分别为定子输出的有功功率和无功功率，Idr和Iqr分别为转子绕组电流的d轴、q轴分量；Ps *和Qs *分别为定子输出的有功功率参考值和无功功率参考值，和分别为转子绕组电流的d轴分量参考值和q轴分量参考值；
构建网侧变流器的动态模型
其中，x5、x6和x7为引入的状态变量，Ki5、Ki6和Ki7相应PI控制器的积分系数；为电容器的直流电压参考值，为网侧变流器交流侧电流的d轴分量参考值，为网侧变流器交流侧电流的q轴分量参考值。
4.根据权利要求3所述的针对次同步振荡的风电并网系统安全运行方法，其特征在于，所述将所述双馈风机的动态模型线性化，得到风电系统小干扰模型，具体包括：
对所述异步发电机的动态模型、所述中间电容器的动态模型、所述滤波电抗的模型、所述转子侧变流器的动态模型和网侧变流器的动态模型进行线性化，得到双馈风机的小干扰模型
ΔIW＝CWΔXW+DWΔUW
其中ΔXW＝[ΔE′d,ΔE′q,Δs,ΔUdc,ΔIdr3,ΔIqr3,Δx1,Δx2,Δx3,Δx4,Δx5,Δx6,Δx7]T为双馈风机状态变量，ΔUW＝[ΔUd,ΔUq]T为双馈风机接入节点电压，ΔIW＝[ΔId,ΔIq]T为双馈风机接入节点注入电流；ΔUd和ΔUq分别为双馈风机接入母线电压的d轴和q轴分量；ΔId和ΔIq分别为双馈风机接入母线注入电流的d轴和q轴分量；AW、BW、CW和DW分别为双馈风机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
线性化后的同步发电机状态方程：
ΔIG＝CGΔXG+DGΔUG
其中ΔXG为同步发电机状态变量，ΔUG为同步发电机接入节点电压，ΔIG为同步发电机接入节点注入电流；AG、BG、CG和DG分别为同步发电机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
考虑负荷的网络方程：
ΔI＝YΔU
ΔI为节点注入电流，Y为考虑负荷的节点导纳矩阵，ΔU为节点电压；
将双馈风机的小干扰模型、线性化后的同步发电机状态方程以及考虑负荷的网络方程联立，得到整个系统的小干扰状态空间模型
Δy＝CΔX+DΔu
其中，ΔX＝[ΔX1,ΔX2…ΔXn]T为n维状态向量，Δu＝[Δu1,Δu2…Δur]T为r维输入向量，Δy＝[Δy1,Δy2…Δym]为m维输出向量，A、B、C和D分别为n×n阶状态矩阵、n×r阶输入矩阵、m×n阶输出矩阵和m×r阶直接传递矩阵。
5.根据权利要求4所述的针对次同步振荡的风电并网系统安全运行方法，其特征在于，所述对所述风电系统小干扰模型进行拉氏变换，并设定初始条件为零，得到待测发电机的输出向量与间谐波电流向量之间的矩阵传递函数，具体包括：
对整个系统的小干扰状态空间模型进行拉氏变换，并设定初始条件为零，得到
ΔX(s)＝(sI-A)-1BΔu(s)
Δy(s)＝CΔX(s)+DΔu(s)
则输出向量Δy与输入向量Δu之间的矩阵传递函数为：
W(s)＝C(sI-A)-1B+D
W(s)为m×r阶矩阵传递函数，矩阵元素Wij(s)表示由输入变量第j个元素Δuj到输出向量第i个元素Δyi的正弦响应；对于角频率为ω的输入变量第j个元素Δuj，输出向量第i个元素Δyi的幅值放大倍数为|Wij(jω)|，相移为∠Wij(jω)，Δyi的总响应等于各个输入响应的线性和，即：
6.一种针对次同步振荡的风电并网系统安全运行系统，其特征在于，包括：
谐波获取模块，用于获取对风电并网系统中的间谐波进行实时监测所得到的间谐波电流，得到间谐波电流向量；
振荡传播因子计算模块，用于利用公式计算振荡传播因子；ki(s)为待测发电机的输出向量中第i个元素所对应的振荡传播因子；|Δyi|为待测发电机的输出向量中第i个元素的幅值，Δuj为间谐波电流向量的第j个元素；||Δu||2为间谐波电流向量的二范数；Wij(s)为间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
判断模块，用于判断所述振荡传播因子是否大于预设阈值，得到判断结果；
切除模块，用于当所述判断结果表示大于时，确定次同步振荡传播到所述待测发电机，将引起次同步振荡的风机从所述风电并网系统中切除；
运行参数调整模块，用于当所述判断结果表示小于或等于时，确定次同步振荡仅在风机组范围内传播，调整风机组的运行参数抑制所述间谐波。
7.根据权利要求6所述的针对次同步振荡的风电并网系统安全运行系统，其特征在于，还包括正弦响应确定模块，用于确定所述间谐波电流向量的第j个元素到待测发电机的输出向量中第i个元素的正弦响应；
所述正弦响应确定模块包括：
动态模型建立单元，用于建立双馈风机的动态模型；
模型线性化单元，用于将所述双馈风机的动态模型线性化，得到风电系统小干扰模型；
矩阵传递函数确定单元，用于对所述风电系统小干扰模型进行拉氏变换，并设定初始条件为零，得到待测发电机的输出向量与间谐波电流向量之间的矩阵传递函数，其中Wij(s)为所述矩阵传递函数中的元素。
8.根据权利要求7所述的针对次同步振荡的风电并网系统安全运行系统，其特征在于，所述动态模型建立单元，具体包括：
异步发电机模型构建子单元，用于构建异步发电机的动态模型
其中，E′d为定子暂态电势的d轴分量，Udr为转子绕组电压的d轴分量，E′q为定子暂态电势的q轴分量、Uqr为转子绕组电压的q轴分量；Ids和Iqs分别为定子绕组电流的d轴、q轴分量，s为转差率，Rr为转子绕组电阻，Xm为激磁电抗，Xrr＝Xm+Xr，其中Xr为转子漏抗，TJ为电机惯性时间常数，Pe为电磁功率，Pm为机械功率，D为阻尼系数，s0为稳态时的转差率，ω0为同步角频率；
电容模型构建子单元，用于考虑直流电压动态特性，构建中间电容器的动态模型
其中，Pr为转子侧变流器输出到中间电容器的有功功率，Pr3为网侧变流器输出到交流侧的有功功率；C为中间电容器的电容值，Udc为电容器的直流电压；
电抗模块构建子单元，用于构建滤波电抗的模型
其中，Idr3为网侧变流器交流侧电流的d轴分量，Ud1为网侧变流器交流侧电压的d轴分量，Iqr3为网侧变流器交流侧电流的q轴分量，Uq1为网侧变流器交流侧电压的q轴分量，Xr3为滤波电抗值；Uds和Uqs分别为定子电压d轴和q轴分量；
转子侧变流器模型构建子单元，用于构建转子侧变流器的动态模型
其中，x1、x2、x3和x4为引入的状态变量，Ki1、Ki2、Ki3和Ki4为相应PI控制器的积分系数，Ps和Qs分别为定子输出的有功功率和无功功率，Idr和Iqr分别为转子绕组电流的d轴、q轴分量；Ps *和Qs *分别为定子输出的有功功率参考值和无功功率参考值，和分别为转子绕组电流的d轴分量参考值和q轴分量参考值；
网侧变流器模型构建子单元，用于构建网侧变流器的动态模型
其中，x5、x6和x7为引入的状态变量，Ki5、Ki6和Ki7相应PI控制器的积分系数；为电容器的直流电压参考值，为网侧变流器交流侧电流的d轴分量参考值，为网侧变流器交流侧电流的q轴分量参考值。
9.根据权利要求8所述的针对次同步振荡的风电并网系统安全运行系统，其特征在于，所述模型线性化单元，具体包括：
线性化子单元，用于对所述异步发电机的动态模型、所述中间电容器的动态模型、所述滤波电抗的模型、所述转子侧变流器的动态模型和网侧变流器的动态模型进行线性化，得到双馈风机的小干扰模型
ΔIW＝CWΔXW+DWΔUW
其中ΔXW＝[ΔEd′,ΔEq′,Δs,ΔUdc,ΔIdr3,ΔIqr3,Δx1,Δx2,Δx3,Δx4,Δx5,Δx6,Δx7]T为双馈风机状态变量，ΔUW＝[ΔUd,ΔUq]T为双馈风机接入节点电压，ΔIW＝[ΔId,ΔIq]T为双馈风机接入节点注入电流；ΔUd和ΔUq分别为双馈风机接入母线电压的d轴和q轴分量；ΔId和ΔIq分别为双馈风机接入母线注入电流的d轴和q轴分量；AW、BW、CW和DW分别为双馈风机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
线性化后的同步发电机状态方程：
ΔIG＝CGΔXG+DGΔUG
其中ΔXG为同步发电机状态变量，ΔUG为同步发电机接入节点电压，ΔIG为同步发电机接入节点注入电流；AG、BG、CG和DG分别为同步发电机模型的状态矩阵、控制矩阵、输出矩阵和直接传递阵；
考虑负荷的网络方程：
ΔI＝YΔU
ΔI为节点注入电流，Y为考虑负荷的节点导纳矩阵，ΔU为节点电压；
多方程联立子单元，用于将双馈风机的小干扰模型、线性化后的同步发电机状态方程以及考虑负荷的网络方程联立，得到整个系统的小干扰状态空间模型
Δy＝CΔX+DΔu
其中，ΔX＝[ΔX1,ΔX2…ΔXn]T为n维状态向量，Δu＝[Δu1,Δu2…Δur]T为r维输入向量，Δy＝[Δy1,Δy2…Δym]为m维输出向量，A、B、C和D分别为n×n阶状态矩阵、n×r阶输入矩阵、m×n阶输出矩阵和m×r阶直接传递矩阵。
10.根据权利要求9所述的针对次同步振荡的风电并网系统安全运行系统，其特征在于，所述矩阵传递函数确定单元，具体包括：
拉氏变换子单元，用于对整个系统的小干扰状态空间模型进行拉氏变换，并设定初始条件为零，得到
ΔX(s)＝(sI-A)-1BΔu(s)
Δy(s)＝CΔX(s)+DΔu(s)
则输出向量Δy与输入向量Δu之间的矩阵传递函数为：
W(s)＝C(sI-A)-1B+D
W(s)为m×r阶矩阵传递函数，矩阵元素Wij(s)表示由输入变量第j个元素Δuj到输出向量第i个元素Δyi的正弦响应；对于角频率为ω的输入变量第j个元素Δuj，输出向量第i个元素Δyi的幅值放大倍数为|Wij(jω)|，相移为∠Wij(jω)，Δyi的总响应等于各个输入响应的线性和，即：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811600781.2A CN109617121B (zh) | 2018-12-26 | 2018-12-26 | 一种针对次同步振荡的风电并网系统安全运行方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811600781.2A CN109617121B (zh) | 2018-12-26 | 2018-12-26 | 一种针对次同步振荡的风电并网系统安全运行方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109617121A true CN109617121A (zh) | 2019-04-12 |
| CN109617121B CN109617121B (zh) | 2020-07-03 |
Family
ID=66012284
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811600781.2A Active CN109617121B (zh) | 2018-12-26 | 2018-12-26 | 一种针对次同步振荡的风电并网系统安全运行方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109617121B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110011357A (zh) * | 2019-04-29 | 2019-07-12 | 东南大学 | 一种双馈风电场中强迫振荡扰动源定位方法 |
| CN110501610A (zh) * | 2019-08-26 | 2019-11-26 | 国网山东省电力公司莱芜供电公司 | 一种多逆变器并网振荡特征实时提取方法及装置 |
| CN110794220A (zh) * | 2019-10-28 | 2020-02-14 | 云南电网有限责任公司 | 一种用于评估风电场汇集系统振荡风险的方法及装置 |
| CN110797908A (zh) * | 2019-08-16 | 2020-02-14 | 南京理工大学 | 一种基于广义短路比法的风电并网低频振荡抑制方法 |
| CN113346513A (zh) * | 2021-03-09 | 2021-09-03 | 江苏大学 | 一种辨识直驱风机强迫次同步振荡的方法 |
| CN115360726A (zh) * | 2022-07-05 | 2022-11-18 | 华北电力大学 | 一种次同步振荡传播关键影响因素识别方法及系统 |
| US11581833B2 (en) * | 2020-07-15 | 2023-02-14 | North China Electric Power University | Active damping control method and system for sub-synchronous oscillation of DFIG, and storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6326796B1 (en) * | 1997-07-07 | 2001-12-04 | Nissin Electric Co., Ltd. | Harmonic measuring method and a current injection device for harmonic measurement |
| CN107994579A (zh) * | 2017-12-18 | 2018-05-04 | 华北电力大学 | 一种直驱风机变流控制系统的间谐波参数稳定域构建方法 |
| CN108011364A (zh) * | 2017-11-28 | 2018-05-08 | 郑州轻工业学院 | 一种分析dfig动力学特性与电力系统动态交互影响的方法 |
-
2018
- 2018-12-26 CN CN201811600781.2A patent/CN109617121B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6326796B1 (en) * | 1997-07-07 | 2001-12-04 | Nissin Electric Co., Ltd. | Harmonic measuring method and a current injection device for harmonic measurement |
| CN108011364A (zh) * | 2017-11-28 | 2018-05-08 | 郑州轻工业学院 | 一种分析dfig动力学特性与电力系统动态交互影响的方法 |
| CN107994579A (zh) * | 2017-12-18 | 2018-05-04 | 华北电力大学 | 一种直驱风机变流控制系统的间谐波参数稳定域构建方法 |
Non-Patent Citations (2)
| Title |
|---|
| ALIREZA ALIZADE 等: "Evaluating noise and DC offset due to inter-harmonics and supra-harmonics caused by back-to-back converter of (DFIG) in AC distribution network", 《CIRED - OPEN ACCESS PROCEEDINGS JOURNAL》 * |
| 祁鑫 等: "基于广域信息的次同步振荡传播分析研究", 《控制工程》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110011357A (zh) * | 2019-04-29 | 2019-07-12 | 东南大学 | 一种双馈风电场中强迫振荡扰动源定位方法 |
| CN110011357B (zh) * | 2019-04-29 | 2022-06-03 | 东南大学 | 一种双馈风电场中强迫振荡扰动源定位方法 |
| CN110797908A (zh) * | 2019-08-16 | 2020-02-14 | 南京理工大学 | 一种基于广义短路比法的风电并网低频振荡抑制方法 |
| CN110797908B (zh) * | 2019-08-16 | 2022-09-06 | 南京理工大学 | 一种基于广义短路比法的风电并网低频振荡抑制方法 |
| CN110501610A (zh) * | 2019-08-26 | 2019-11-26 | 国网山东省电力公司莱芜供电公司 | 一种多逆变器并网振荡特征实时提取方法及装置 |
| CN110501610B (zh) * | 2019-08-26 | 2021-10-01 | 国网山东省电力公司莱芜供电公司 | 一种多逆变器并网振荡特征实时提取方法及装置 |
| CN110794220A (zh) * | 2019-10-28 | 2020-02-14 | 云南电网有限责任公司 | 一种用于评估风电场汇集系统振荡风险的方法及装置 |
| US11581833B2 (en) * | 2020-07-15 | 2023-02-14 | North China Electric Power University | Active damping control method and system for sub-synchronous oscillation of DFIG, and storage medium |
| CN113346513A (zh) * | 2021-03-09 | 2021-09-03 | 江苏大学 | 一种辨识直驱风机强迫次同步振荡的方法 |
| CN113346513B (zh) * | 2021-03-09 | 2023-10-10 | 江苏大学 | 一种辨识直驱风机强迫次同步振荡的方法 |
| CN115360726A (zh) * | 2022-07-05 | 2022-11-18 | 华北电力大学 | 一种次同步振荡传播关键影响因素识别方法及系统 |
| CN115360726B (zh) * | 2022-07-05 | 2024-06-25 | 华北电力大学 | 一种次同步振荡传播关键影响因素识别方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109617121B (zh) | 2020-07-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109617121A (zh) | 一种针对次同步振荡的风电并网系统安全运行方法及系统 | |
| Li et al. | Stability control for wind in weak grids | |
| Kunjumuhammed et al. | Electrical oscillations in wind farm systems: Analysis and insight based on detailed modeling | |
| CN102074967B (zh) | 一种具有并网特性的储能型风电场控制方法 | |
| CN105958543B (zh) | 一种提升新能源并网稳定性的控制、实验和仿真方法 | |
| CN102624013B (zh) | 基于相位补偿原理的储能阻尼控制器的设计方法 | |
| CN109980691A (zh) | 一种计及频域阻抗稳定性分析的双馈风电机组虚拟同步控制方法 | |
| Du et al. | Stability analysis for DFIG-based wind farm grid-connected system under all wind speed conditions | |
| CN111130136A (zh) | 一种基于附加虚拟阻抗控制的次同步振荡抑制方法 | |
| Li et al. | Structure preserving aggregation method for doubly-fed induction generators in wind power conversion | |
| Zhu et al. | High-frequency resonance of DFIG-based wind generation under weak power network | |
| Lin et al. | A Voltage Control Strategy of VSG Based on Self‐Adaptive Inertia Coefficient and Droop Coefficient | |
| Jazaeri et al. | Eigenvalue Analysis of a Network Connected to a Wind Turbine Implemented with a Doubly-Fed Induction Generator (DFIG) | |
| Samanes et al. | Sub-synchronous resonance damper based on the stator voltage feedback for DFIG wind turbines | |
| Wang et al. | Effect of AC voltage control on the stability of weak AC grid connected DFIG system | |
| Song et al. | Sensitivity analysis of the wind farm high frequency resonance under transmission cable resistance variation | |
| Yang et al. | Analysis of subsynchronous oscillation propagation characteristics caused by interharmonics in wind integrated power system | |
| Wang et al. | Impedance modeling of wind turbine-variable speed pumped storage combined operation system | |
| Du et al. | A virtual synchronous generator adopting dynamic damping without frequency detection | |
| Riouch et al. | A coordinated control for smoothing output power of a DFIG based wind turbine | |
| Charles et al. | Field oriented control of Doubly Fed Induction Generator in wind power system | |
| Yan et al. | Impact of DFIG-based WTs on subsynchronous oscillation damping of SG based on the net electrical damping method | |
| Li et al. | Impedance Analysis and Stability Optimization based on Phase-Optimization PLL of VSC-HVDC Connected to Weak AC Grids | |
| CN110365048A (zh) | 一种基于最优化潮流的黑启动网架动态稳定性评估方法 | |
| Zhu et al. | A DFIG Harmonic Suppression Method Based on Resonance Control in DC Distribution Network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |