CN105140957B - 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 - Google Patents
基于风电场和光伏电站聚合模型的机电振荡模式估算方法 Download PDFInfo
- Publication number
- CN105140957B CN105140957B CN201510513100.9A CN201510513100A CN105140957B CN 105140957 B CN105140957 B CN 105140957B CN 201510513100 A CN201510513100 A CN 201510513100A CN 105140957 B CN105140957 B CN 105140957B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- photovoltaic power
- mover
- msubsup
- 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.)
- Expired - Fee Related
Links
- 230000010355 oscillation Effects 0.000 title claims abstract description 27
- 238000006116 polymerization reaction Methods 0.000 title abstract description 8
- 230000005662 electromechanics Effects 0.000 title 1
- 238000011156 evaluation Methods 0.000 title 1
- 238000012546 transfer Methods 0.000 claims abstract description 68
- 238000000034 method Methods 0.000 claims abstract description 42
- 230000004044 response Effects 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims abstract description 5
- 238000010248 power generation Methods 0.000 claims description 56
- 239000011159 matrix material Substances 0.000 claims description 31
- 230000002776 aggregation Effects 0.000 claims description 26
- 238000004220 aggregation Methods 0.000 claims description 26
- 238000010586 diagram Methods 0.000 claims description 9
- 230000008846 dynamic interplay Effects 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 5
- 239000013598 vector Substances 0.000 claims description 5
- 238000010587 phase diagram Methods 0.000 claims description 4
- 238000005286 illumination Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 18
- 230000005611 electricity Effects 0.000 abstract 1
- 230000003534 oscillatory effect Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 11
- 230000001360 synchronised effect Effects 0.000 description 7
- 238000011217 control strategy Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 239000002131 composite material Substances 0.000 description 3
- 238000009795 derivation Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000013016 damping Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 241000022852 Letis Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- 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/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/16—Energy services, e.g. dispersed generation or demand or load or energy savings aggregation
Landscapes
- Photovoltaic Devices (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明属于新能源电力系统动态稳定性分析技术领域，尤其涉及一种基于大规模风电场和光伏电站聚合模型的电力系统机电振荡模式估算方法，首先建立单台直驱风机或光伏发电单元的传递函数模型，整个风电场和光伏电站的聚合模型可以通过单台发电设备的模型直接倍乘，或相加后降阶的方法得到；对于参数未知的风电场或光伏电站，聚合的传递函数模型通过现场频率响应试验得到，电力系统机电振荡模式由系统的开环特征值加上聚合传递函数值与残差的乘积来估算。本发明能够在保证一定准确度的情况下，大大降低模型阶数和计算复杂度，从而为新能源发电电源接入后的电力系统小干扰振荡稳定性分析提供帮助。
Description
技术领域
本发明属于新能源电力系统动态稳定性分析技术领域，尤其涉及一种基于大规模风电场和光伏电站聚合模型的电力系统机电振荡模式估算方法。
背景技术
近年来，随着风电机组技术和光伏发电技术的不断进步、风电场和光伏电站规模的不断扩大、区域电网互联的飞速发展，新能源发电电源的接入对电网稳定性的影响越来越大。在这种发展趋势下，全国各地区电网公司对大规模新能源发电电源接入电网所带来的影响，尤其是小干扰稳定性问题非常重视。因此，对于大规模风电场和光伏电站接入电力系统后进行小干扰稳定性及阻尼特性的研究成为了当今最迫切的课题。
目前，在分析电力系统小干扰稳定性问题时，特征值分析法得到了广泛应用。特征值分析法是有效解决计算电力系统机电振荡模式的定量信息问题最有效的手段。它是借助线性系统的稳定性理论来研究非线性系统的稳定性。该方法严格建立在现代控制理论基础上，把电力系统视为一般控制系统，用标准的线性化状态方程来描述。系统的振荡模式完全由状态矩阵的特征根的对数来决定。特征根的实部和虚部分别对应振荡模式的阻尼和频率，而特征根对应的特征向量反映了振荡模式在整个系统中的动态行为。因此，不论系统有多么复杂，都可以采用特征值分析法对系统中任一振荡模式做详细分析，而与其它振荡模式无关。特征值分析法不仅可以用来刻画系统的振荡模式，而且可以用来确定控制器最佳安装地点。
当大规模风电场或光伏电站接入电力系统后，使用特征值分析法分析电力系统机电振荡模式时，一个最大的障碍就是状态矩阵阶数的增加大大增加了计算特征值的复杂度和计算时间。因此对可降低系统阶数的风电场和光伏电站的聚合模型的研究成为新能源电力系统小干扰稳定性分析中很重要的工作。
风电场的聚合模型是指将风电场等效成为单台或几台机组。显然，当风电场包含几百台风电机组时，采用详细模型会使计算规模过于庞大，计算能力往往不能满足要求。因此，当研究风电场对电网的影响时，可将风电场看作一个整体进行聚合，以减小仿真规模，提高计算效率和收敛性。目前，建立风电场聚合模型的方法可以归纳为3种：
1)单机法(完全聚合法)。将风电场等值成一台风力发电机组。等值风电机组的容量等于所有风电机组容量的代数和。其余参数都由等值计算获得。
2)复合模型法。将风电场电气部分等值为一台发电机，而空气动力学部分、传动链和最大功率跟踪模型则根据运行条件分区分别聚合成几台并联，构成一台复合的风电机组。
3)分群法(多机法)。按一定分群原则，将风电场等值成为多台机组。
上述三种方法中，单机法简化程度最大，但往往被认为精度较低；而复合模型法和分群法是为了体现风场中机组运行状态的差异性对输出特性带来的影响，等效机组过多时，简化程度不够。而且上述方法都用于新能源接入后电力系统的动态分析中，而用在小干扰分析中时不免“大材小用”，因为它们都提供了一些无用于小干扰分析的状态量与参数。
发明内容
为了在新能源电源接入电力系统后，建立用于新能源电力系统小干扰稳定分析的风电场和光伏电站聚合模型，以降低特征值分析中的模型阶数和计算复杂度，本发明提出了一种基于大规模风电场和光伏电站聚合模型的电力系统机电振荡模式估算方法，包括：
步骤1、建立大规模风电场和光伏电站的聚合传递函数模型；
步骤2、将电力系统中风电场和光伏电站的输出电流作为控制量，建立开环的多机电力系统的线性化模型；
步骤3、把风电场或光伏电站的传递函数模型看作是控制器，作为连接开环系统输入与输出的反馈回路；开环系统和反馈回路共同构成了多机电力系统的闭环状态空间模型；
步骤4、将风电场和光伏电站看作恒电流源，即令风电场或光伏电站输出电流为零，忽略掉风电场和光伏电站与多机电力系统间的动态交互影响，计算开环系统的特征值；通过系统的开环特征值加上聚合传递函数值与残差的乘积来估算电力系统闭环的机电振荡模式。
所述步骤1若在各台直驱风力发电机或光伏发电单元的型号、参数和运行工况都已知的情况下，则建立的各台直驱风机和光伏发电单元的传递函数模型为：
其中，和分别为第j个风电场或光伏电站内第k台直驱风机或光伏发电单元输出电流的d轴和q轴分量，|Vwj|为第j个风电场或光伏电站内低压母线的电压幅值，为第j个风电场或光伏电站内第k个直驱风机或光伏发电单元的外部变量；对于直驱风机，外部变量为机械转矩，对于光伏发电单元，外部变量为光照强度；和分别为所对应变量间的传递函数。
所述步骤1若在一个风电场或光伏电站内所有直驱风机或光伏发电单元的型号、参数和运行工况都相同的情况下，那么该风电场或光伏电站的聚合传递函数模型由单台风机或光伏发电单元的传递函数倍乘直接得到：
若一个风电场或光伏电站内直驱风机或光伏发电单元的型号、参数或运行工况不同，那么整个风电场或光伏电站的聚合传递函数模型由所有风机或光伏发电单元的传递函数相加得到：
其中，Iwdj和Iwqj分别为第j个风电场或光伏电站内总输出电流的d轴和q轴分量，Lj为第j个风电场或光伏电站内所包含直驱风机或光伏发电单元的数量，Gdj(s)和Gqj(s)分别为第j个风电场或光伏电站内输出量ΔIwdj和ΔIwqj对输入量Δ|Vwj|的聚合传递函数；
对于式(25)，传递函数Gdj(s)和Gqj(s)的阶数会很高，通过目前已有的模型降阶方法，将高阶的传递函数模型降阶成低阶的传递函数模型。
所述步骤1中在直驱风机或光伏发电单元的型号、参数或运行工况未知的情况下，通过现场频率响应试验来得到风电场和光伏电站的聚合传递函数模型；通过在低压母线稳态电压幅值上叠加某频率的正弦信号作为输入，检测风电场或光伏电站的输出电流，得到该频率下的输出电流对输入电压的传递函数值；对各个频率下传递函数值的幅值图和相位图，进行不同阶数的拟合，从而近似得到风电场和光伏电站的聚合模型。
所述步骤2中建立开环的多机电力系统的线性化模型为：
其中，Xg为系统状态变量，A0为开环系统的状态矩阵，b0dj和b0qj为控制矩阵，为输出矩阵，d0dj和d0qj为反馈系数，M为多机电力系统所包含风电场和光伏电站的数量，Iwdj和Iwqj分别为第j个风电场或光伏电站总输出电流的d轴和q轴分量，|Vwj|为第j个风电场或光伏电站内低压母线的电压幅值。
所述步骤4中电力系统机电振荡模式为：
残差为：
其中，为开环系统的状态矩阵的特征值，和分别为对应的左右特征向量，为考虑了风电场和光伏电站动态后闭环阵的特征值估计值，Gdj(s)和Gqj(s)分别为第j个风电场或光伏电站的输出量ΔIwdj和ΔIwqj对输入量Δ|Vwj|的聚合传递函数，b0dj和b0qj为控制矩阵，为输出矩阵，d0dj和d0qj为反馈系数，M为多机电力系统所包含风电场和光伏电站的数量。
本发明的有益效果在于：用本发明中提出的方法计算大规模风电场和光伏电站接入电力系统后的机电振荡模式时，不需要推导闭环系统的状态矩阵。只需要得出把风电场和光伏电站等效为恒电流源后的开环状态矩阵的特征值，再计算残差，即可近似估算闭环系统的机电振荡模式。这种方法可以大大降低矩阵的阶数，减小模态计算中计算负担和复杂度。
附图说明
图1为直驱风机的结构图；
图2为直驱风机机侧换流器的控制方案；
图3为直驱风机网侧换流器的控制方案；
图4为光伏发电单元的结构图；
图5为光伏发电单元DC/DC换流器的控制方案；
图6为大规模风电场和光伏电站接入后的多机电力系统的结构图；
图7为x-y坐标系与d-q坐标系的旋转关系；
图8为大规模风电场和光伏电站接入后的多机电力系统的线性化模型示意图；
图9为纽约-新英格兰互联电力系统的结构图；
图10为连接在母线48处的光伏电站的频率响应曲线。
具体实施方式
下面结合附图，对实施例作详细说明。
本发明提出了一种基于大规模风电场和光伏电站聚合模型的电力系统机电振荡模式估算方法，包括：
步骤1、建立大规模风电场和光伏电站的聚合传递函数模型；
步骤2、每台直驱风机或光伏发电单元通过变压器连接到低压母线上，然后再通过升压变压器连接到公共连接点处的高压母线上；建立开环的多机电力系统的线性化模型；
步骤3、把不考虑风电场和光伏电站动态的多机电力系统的状态空间模型，看作是开环的多机电力系统模型；风电场或光伏电站的传递函数模型看作是控制器，作为连接开环系统输入与输出的反馈回路；开环系统和反馈回路共同构成了多机电力系统的闭环状态空间模型；
步骤4、令风电场或光伏电站输出电流的d轴和q轴分量均为零，忽略掉风电场和光伏电站与多机电力系统间的动态交互影响，计算开环系统、反馈回路、闭环系统的特征值，并估算电力系统机电振荡模式。
所述步骤1，建立大规模风电场和光伏电站的聚合传递函数模型，具体包括下面两种情况：
(1)在各台直驱风力发电机或光伏发电单元的型号、参数和运行工况都已知的情况下，首先建立各台直驱风机和光伏发电单元的传递函数模型。
直驱风机的结构如附图1所示。
在d-q坐标系下，直驱风机的定子电压方程、定子磁链方程和转子运动方程分别为：
其中，ω0为同步角速度，λ0是永磁体的恒定磁链，Ra、Xz、ψpmz、Vpmz和Ipmz(z＝d,q)分别是定子电阻、定子漏抗、定子磁通、定子电压和定子电流，Jpm为直驱风机的惯性常数，ωpm为直驱风机的转子转速，Twm和Tpme分别为机械转矩和电磁转矩，其中：
Tpme＝ψpmqIpmd-ψpmdIpmq (4)
忽略定子电阻，同时假设Xd＝Xq＝Xpm，并把式(2)代入式(1)，可得：
直驱风机的机侧换流器的控制策略如附图2所示。忽略其中的电磁暂态快速过程，那么控制策略中的内环控制可以忽略，由附图2的外环控制可得：
在网侧换流器的控制中，将低压母线电压的方向定为d-q坐标系的q轴方向。由附图1可得：
直驱风机的网侧换流器的控制策略如附图3所示。同样忽略其中的电磁暂态快速过程，那么控制策略中的内环控制可以忽略，由附图3的外环控制可得：
电容器的动态方程为：
其中，
Ppme＝IpmdVpmd+IpmqVpmq
Pw＝|Vw|Iwq (10)
Qw＝|Vw|Iwd
对式(1)-(6)进行线性化，可得：
ΔIpmq＝Kpm(s)Δωpm
JpmsΔωpm＝ΔTwm-λ0ΔIpmq (11)
ΔPpme＝λ0ωpm0ΔIpmq+λ0Ipmq0Δωpm
对式(8)-(10)进行线性化，可得：
ΔIwq＝-Kdc(s)ΔVdc
ΔIwd＝-Kq(s)ΔQw
ΔPw＝Iwq0Δ|Vw|+|Vw0|ΔIwq
ΔQw＝Iwd0Δ|Vw|+|Vw0|ΔIwd
由式(11)(12)可得：
ΔPw＝a1(s)Δ|Vw|+a2(s)ΔTwm
ΔQw＝a3(s)Δ|Vw| (13)
ΔVdc＝a4(s)Δ|Vw|+a5(s)ΔTwm
其中，
将式(13)代入式(12)中，可得：
其中，
b1(s)＝-Kq(s)a3(s)
b2(s)＝-Kdc(s)a4(s)
b3(s)＝-Kdc(s)a5(s)
光伏发电单元的结构如附图4所示。
光伏发电单元的电压-电流(Vpv-Ipv)特性和出力分别为：
Ppv＝IpvVpv (16)
其中，T为光伏板工作温度，Ir为光照强度，Ns和Np分别是串并联光伏板的单元数量，n为理想化因子，k为波尔兹曼常数，q为电子的电荷量，Isc为短路电流，I0为饱和电流。
附图4中左侧电感的动态方程为：
光伏发电单元DC/DC换流器的控制策略如附图5所示，由附图5可得：
对式(15)-(18)进行线性化，可得：
ΔVpv＝aIΔIpv+arΔIr
ΔPpv＝Vpv0ΔIpv+Ipv0ΔVpv
ΔVdc1＝Kpv2(s)ΔIpv+ΔVpv
其中，
继而有：
ΔVdc1＝aVIΔIpv+aVrΔIr
其中，
aVI＝Kpv2(s)+Kpv1(s)Kpv2(s)Vpv0+(1+Kpv1(s)Kpv2(s)Vpv0)aI
aVr＝(1+Kpv1(s)Kpv2(s)Ipv0)ar
光伏发电单元的DC/AC换流器的控制策略与直驱风机的网侧换流器的控制策略相同，如附图3所示。其线性化模型的推导见式(12)-(14)，只需要将其中的ΔPpme和ΔTwm替换为ΔPpv和ΔIr即可：
ΔPw＝a1(s)Δ|Vw|+a2(s)ΔIr
ΔQw＝a3(s)Δ|Vw| (21)
ΔVdc＝a4(s)Δ|Vw|+a5(s)ΔIr
其中，a1(s)、a3(s)和a4(s)与式(13)中的相同，其余参数为：
将式(21)代入到式(12)中，可得：
其中，b1(s)、b2(s)和b3(s)与式(14)中相同。
自此，对于如附图6所示的多机电力系统，在第j个风电场或光伏电站内第k个直驱风机或光伏发电单元的输出电流与低压母线电压幅值之间的关系可以统一表示为：
其中，对于直驱风机，外部变量Δu为机械转矩ΔTwm；对于光伏发电单元，外部变量Δu为光照强度ΔIr。
第j个风电场或光伏电站的总输出电流为：
Gdj(s)和Gqj(s)分别为第j个风电场或光伏电站内，输出量ΔIwdj和ΔIwqj对输入量Δ|Vwj|的聚合传递函数：
当风电场或光伏电站内所有直驱风机或光伏发电单元的型号、参数和运行工况完全相同时，式(25)所示的聚合传递函数可以写为：
显然当风电场或光伏电站内有部分直驱风机或光伏发电单元的型号、参数和运行工况近似相同时，式(25)同样可以用来降低聚合传递函数的阶数。
当风电场或光伏电站内直驱风机或光伏发电单元的型号、参数和运行工况不同时，由式(25)所得的聚合传递函数的阶数可能会非常高。此时利用成熟的模型降阶方法，可以将高阶的传递函数模型降阶成低阶的传递函数模型，这样可以大大降低电力系统小干扰稳定分析的计算压力和复杂度。
(2)在直驱风机或光伏发电单元的型号、参数或运行工况未知的情况下，可以通过现场频率响应试验来得到风电场和光伏电站的聚合传递函数模型。通过在低压母线稳态电压幅值上叠加某频率的正弦信号作为输入，检测整个风电场或光伏电站的输出电流，可以得到该频率下的输出电流对输入电压的传递函数值Gdj(jωi)＝|Gdj(jω)|∠Gdj(jω)和Gqj(jωi)＝|Gqj(jω)|∠Gqj(jω)。将一组连续频率下的传递函数值的幅值和相位作图，形成传递函数幅值图和相位图。然后对幅值图和相位图可以进行不同阶数的曲线拟合，从而可以近似得到风电场和光伏电站的聚合模型。模型降阶的程度可以在曲线拟合过程中给予考虑。
所述步骤2具体包括以下步骤：
建立开环的多机电力系统的线性化模型：
如附图6所示的多机电力系统，有N台同步发电机和M个大规模风电场或光伏电站。第j个风电场或光伏电站内有Lj台直驱风机或光伏发电单元。每台直驱风机或光伏发电单元通过变压器连接到低压母线上，然后再通过升压变压器连接到公共连接点(PCC)处的高压母线上。
含有N台同步发电机的多机电力系统中同步发电机的线性化模型为：
其中，Xg为同步机的状态变量，Vg和Ig为N台同步发电机的端电压和输入网络的电流，其每个元素值为：
其中x和y表示在x-y坐标系下相应物理量的x分量和y分量。
用节点导纳矩阵表示的网络方程为：
其中，Vw和Iw为M个风电场或光伏电站的端电压和输入网络的电流，其每个元素值为：
电压相量对其进行线性化可得：
将式(31)写成矩阵形式为：
其中Δ|Vw|和Δθw的第j个元素为Δ|Vwj|和Δθwj。
对于第j个风电场和光伏电站，令的方向为d-q坐标系中q轴的方向，如附图7所示。由附图7可得：
将式(33)写成矩阵形式为：
因此有：
ΔIw＝CdqΔIwdq+CvθΔVw (35)
其中ΔIwdq的第j个元素值为
ΔIwdqj＝[ΔIwdj ΔIwqj]T (36)
由式(27)和式(29)可得：
由式(37)可得：
将式(38)代入式(37)的第二式中，可得：
ΔVw＝(Yww -1-Yww -1YwgCgw)ΔIw-Yww -1YwgCggΔXg (39)
将式(39)代入式(35)中可得：
ΔIw＝BdqΔIwdq+BwgΔXg (40)
其中，
Bdq＝[I-Cvθ(Yww -1-Yww -1YwgCgw)]-1Cdq
Bwg＝-[I-Cvθ(Yww -1-Yww -1YwgCgw)]-1CvθYww -1YwgCgg
将式(39)(40)代入式(32)中，可得：
其中，
C0V＝A1x[(Yww -1-Yww -1YwgCgw)Bwg-Yww -1YwgCgg]
D0V＝A1x(Yww -1-Yww -1YwgCgw)Bdq
C0θ＝A2x[(Yww -1-Yww -1YwgCgw)Bwg-Yww -1YwgCgg]
D0θ＝A2x(Yww -1-Yww -1YwgCgw)Bdq
由式(38)和式(40)可得：
ΔVg＝(Cgg+CgwBwg)ΔXg+CgwBdqΔIwdq (42)
将式(42)代入式(27)的第一式中，可得：
sΔXg＝A0ΔXg+B0ΔIwdq (43)
其中，
A0＝Ag+Bg(Cgg+CgwBwg)
B0＝BgCgwBdq
根据式(36)中的定义，式(43)和式(41)可以写为：
使用上文所述的方法建立各个风电场和光伏电站的聚合的传递函数模型：
ΔIwdj＝Gdj(s)Δ|Vwj|
当风电场或光伏电站的参数或运行工况未知时，采用下面的近似来得到聚合的传递函数模型：
步骤3具体包括：
式(44)和式(45)构成了含M个风电场或光伏电站的多机电力系统的线性化模型，如附图8所示。附图8中的上半部分为不考虑风电场和光伏电站动态的多机电力系统的状态空间模型，它可以看作是开环的多机电力系统模型。由式(45)表示的风电场或光伏电站的传递函数模型可以看作是“控制器”，作为连接开环系统输入与输出的反馈回路。开环系统和反馈回路共同构成了考虑M个风电场或光伏电站动态的多机电力系统的闭环状态空间模型。
大规模风电场和光伏电站的接入对多机电力系统的影响体现在两个方面：一是风电场和光伏电站向原多机电力系统提供一定的注入电流，改变了原多机电力系统的潮流分布，这部分影响体现在式(44)中的A0矩阵中；另一方面是风电场和光伏电站与多机电力系统的动态交互。如附图8所示，风电场和光伏电站向开环系统输入ΔIwdj和ΔIwqj，ΔIwdj和ΔIwqj正是风电场和光伏电站对Δ|Vwj|的响应，而Δ|Vwj|又是开环系统对ΔIwdj和ΔIwqj的响应。显然，开环系统与反馈回路的动态交互过程与外部变量没有任何关系。因此与外部变量相关的动态过程不会影响到电力系统的小干扰功角稳定，在下面的小信号分析中可以直接忽略掉
步骤4具体包括：
令ΔIwdj＝ΔIwqj＝0，此时风电场和光伏电站褪化为恒电流源，此时系统的线性化模型褪化为：
sΔXg＝A0ΔXg (47)
A0矩阵中包含了大规模风电场和光伏电站的接入对原多机电力系统带来的潮流影响，而忽略了风电场和光伏电站与多机电力系统间的动态交互影响。
对于输入为uj、输出为yj的线性系统：
作为反馈回路的控制器的传递函数模型为：
uj＝hj(s)yj (49)
如果i＝1,2,…,N为式(48)中开环矩阵A的第i个特征值，和为其对应的左右特征向量，由控制器造成的特征值的变化量为：
其中是特征值i＝1,2,…,N对应的残差。
如果式(48)表示的开环系统变为：
那么由控制器造成的特征值的变化量为：
闭环系统的特征值为：
对比式(44)(45)和式(51)(49)，可以看出，式(44)中的ΔIwdj和ΔIwqj相当于式(51)中的控制量uj，Δ|Vwj|相当于式(51)中的输出量yj，式(45)中的Gdj(s)和Gqj(s)相当于式(49)中控制器的传递函数hj(s)。根据式(53)，大规模风电场和光伏电站接入后的多机电力系统的机电振荡模式近似为：
残差为
其中，为式(44)中开环阵A0的特征值，和为其对应的左右特征向量，为考虑了风电场和光伏电站动态后闭环阵的特征值估计值。
如附图9所示为纽约-新英格兰互联电力系统，下面通过该系统来说明本发明提出的基于大规模风电场和光伏电站聚合模型的电力系统机电振荡模式估算方法。
该电力系统中网络参数和同步发电机的参数见文献(G.Rogers,Power SystemOscillations. Norwell, MA: Kluwer, 2000.)。系统中有两个风电场(连接在节点9和26处)和两个光伏电站(连接在节点16和48处)。每个风电场内有50台直驱风力发电机，每个光伏电站内有50台光伏发电单元。
假设节点9处风电场内全部直驱风机的参数都相同，为：
ωpm＝0.2,Xd＝0.01,Xq＝0.01,Xw＝0.3,λ0＝1,Jpm＝6s,Cdc＝50,Vdc＝1,
通过式(1)到式(14)的推导，可以得到该风电场内第k台直驱风机的传递函数模型可以表示为：
其中，
则整个风电场的传递函数模型为：
节点26处风电场内，有20台直驱风机的参数相同，为：
ωpm＝0.2,Xd＝0.01,Xq＝0.01,Xw＝0.3,λ0＝1,Jpm＝6s,Cdc＝50,Vdc＝1,
这20台风机中每台风机的传递函数为：
其中，
另30台直驱风机的参数相同，为：
ωpm＝0.2,Xd＝0.01,Xq＝0.01,Xw＝0.3,λ0＝1,Jpm＝6s,Cdc＝50,Vdc＝1,
这30台风机中每台风机的传递函数为：
其中，
那么该风电场的聚合传递函数为：
节点16处光伏电站内，有20台光伏发电单元的参数相同，为：
通过式(15)到式(22)的推导，可以得到这20台光伏发电单元中每台光伏发电单元的传递函数为：
其中，
另30台光伏发电单元的参数相同，为：
这30台光伏发电单元中每台光伏发电单元的传递函数为：
其中，
那么该光伏电站的聚合传递函数为：
节点48处光伏电站内光伏发电单元的参数未知，需要通过现场频率响应试验来得出其聚合模型。要得到其聚合传递函数模型，需要令输入信号为光伏电站内低压母线的电压幅值Δ|Vwj|，输出信号为整个光伏电站输出到网络中的电流ΔIwd和ΔIwq。通过在输入的低压母线稳态电压幅值上叠加某频率的正弦信号，可以检测到带有一定噪声的风电场或光伏电站的输出电流，即得到该频率下的输出电流ΔIwd和ΔIwq对输入电压Δ|Vwj|的传递函数值。将各个频率下传递函数的幅值和相位作图，进行适当降噪处理后，即得到如附图10所示的幅值图和相位图。最后对附图10所示的曲线进行不同阶数的曲线拟合，从而可以近似得到节点48处连接的光伏电站的聚合传递函数模型：
ΔIwd4＝Gwd4(s)|ΔVw|
ΔIwq4＝Gwq4(s)|ΔVw|+Gwe4(s)ΔIr
其中，
下面验证本发明提出的大规模风电场和光伏电站聚合模型的正确性和计算机电振荡模式的方法的正确性。
表1不同方法求得的特征值
(1)当如附图9所示的互联电力系统中所有直驱风机和光伏发电单元全部使用各自的全模型，即直驱风机使用5阶模型，光伏发电单元使用6阶模型时，此时通过求闭环状态矩阵，可以求得特征值如表1的第一列所示。用直驱风机和光伏发电单元的全模型时，计算最为准确，但是状态矩阵阶数和计算复杂度也是最高的，此时由于风电场和光伏电站的接入而增加的状态矩阵阶数为5×100+6×100＝1100。
(2)当如附图9所示的互联电力系统中所有直驱风机和光伏发电单元全部使用本发明中提出的聚合模型时，即节点9处的风电场为3阶模型、节点26处的风电场为6阶模型、节点16处的光伏电站为6阶模型、节点48处的光伏电站为3阶模型，此时由于风电场和光伏电站的接入而增加的状态矩阵阶数为3+6+6+3＝18。此时通过求闭环状态矩阵而得到的系统特征值如表1中的第二列所示。
(3)通过式(27)到式(44)的推导，可以得到开环特征矩阵，开环特征值即可以求得；再利用式(54)和式(55)，可以得到闭环特征矩阵的估算值，如表1的第三列所示。这种方法并没有由于风电场和光伏电站的接入而增加状态矩阵的阶数。
通过对比表1中三列的数值，可以看出，本发明中提出的大规模风电场和光伏电站的聚合模型是正确的，而且本发明中提出的机电振荡模式估算方法也是比较准确的，可以在工程实际中应用。
此实施例仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
Claims (6)
1.一种基于风电场和光伏电站聚合模型的电力系统机电振荡模式估算方法，其特征在于，包括：
步骤1、建立大规模风电场和光伏电站的聚合传递函数模型；
步骤2、将电力系统中风电场和光伏电站的输出电流作为控制量，建立开环的多机电力系统的线性化模型；
步骤3、把风电场或光伏电站的传递函数模型看作是控制器，作为连接开环系统输入与输出的反馈回路；开环系统和反馈回路共同构成了多机电力系统的闭环状态空间模型；
步骤4、将风电场和光伏电站看作恒电流源，即令风电场或光伏电站输出电流为零，忽略掉风电场和光伏电站与多机电力系统间的动态交互影响，计算开环系统的特征值；通过系统的开环特征值加上聚合传递函数值与残差的乘积来估算电力系统闭环的机电振荡模式。
2.根据权利要求1所述方法，其特征在于，所述步骤1若在各台直驱风力发电机或光伏发电单元的型号、参数和运行工况都已知的情况下，则建立的各台直驱风机和光伏发电单元的传递函数模型为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>d</mi>
<mi>j</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>b</mi>
<mn>1</mn>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mo>|</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>q</mi>
<mi>j</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>b</mi>
<mn>2</mn>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mo>|</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>+</mo>
<msubsup>
<mi>b</mi>
<mn>3</mn>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&Delta;u</mi>
<mi>j</mi>
<mi>k</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，和分别为第j个风电场或光伏电站内第k台直驱风机或光伏发电单元输出电流的d轴和q轴分量，|Vwj|为第j个风电场或光伏电站内低压母线的电压幅值，为第j个风电场或光伏电站内第k个直驱风机或光伏发电单元的外部变量；对于直驱风机，外部变量为机械转矩，对于光伏发电单元，外部变量为光照强度；和分别为所对应变量间的传递函数。
3.根据权利要求1所述方法，其特征在于，所述步骤1若在一个风电场或光伏电站内所有直驱风机或光伏发电单元的型号、参数和运行工况都相同的情况下，那么该风电场或光伏电站的聚合传递函数模型由单台风机或光伏发电单元的传递函数倍乘直接得到：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mo>|</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>L</mi>
<mi>j</mi>
</msub>
<msubsup>
<mi>b</mi>
<mn>1</mn>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mo>|</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>L</mi>
<mi>j</mi>
</msub>
<msubsup>
<mi>b</mi>
<mn>2</mn>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>26</mn>
<mo>)</mo>
</mrow>
</mrow>
若一个风电场或光伏电站内直驱风机或光伏发电单元的型号、参数或运行工况不同，那么整个风电场或光伏电站的聚合传递函数模型由所有风机或光伏发电单元的传递函数相加得到：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mo>|</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<msubsup>
<mi>b</mi>
<mn>1</mn>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mo>|</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<msubsup>
<mi>b</mi>
<mn>2</mn>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
1
其中，Iwdj和Iwqj分别为第j个风电场或光伏电站内总输出电流的d轴和q轴分量，Lj为第j个风电场或光伏电站内所包含直驱风机或光伏发电单元的数量，Gdj(s)和Gqj(s)分别为第j个风电场或光伏电站内输出量ΔIwdj和ΔIwqj对输入量Δ|Vwj|的聚合传递函数；|Vwj|为第j个风电场或光伏电站内低压母线的电压幅值；
对于式(25)，传递函数Gdj(s)和Gqj(s)的阶数会很高，通过目前已有的模型降阶方法，将高阶的传递函数模型降阶成低阶的传递函数模型。
4.根据权利要求1所述方法，其特征在于，所述步骤1中在直驱风机或光伏发电单元的型号、参数或运行工况未知的情况下，通过现场频率响应试验来得到风电场和光伏电站的聚合传递函数模型；通过在低压母线稳态电压幅值上叠加某频率的正弦信号作为输入，检测风电场或光伏电站的输出电流，得到该频率下的输出电流对输入电压的传递函数值；对各个频率下传递函数值的幅值图和相位图，进行不同阶数的拟合，从而近似得到风电场和光伏电站的聚合模型。
5.根据权利要求1所述方法，其特征在于，所述步骤2中建立开环的多机电力系统的线性化模型为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>s&Delta;X</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&Delta;X</mi>
<mi>g</mi>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mrow>
<mn>0</mn>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mn>0</mn>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mo>|</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>=</mo>
<msubsup>
<mi>c</mi>
<mrow>
<mn>0</mn>
<mi>V</mi>
<mi>j</mi>
</mrow>
<mi>T</mi>
</msubsup>
<msub>
<mi>&Delta;X</mi>
<mi>g</mi>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mrow>
<mn>0</mn>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mrow>
<mn>0</mn>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;I</mi>
<mrow>
<mi>w</mi>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>44</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，Xg为系统状态变量，A0为开环系统的状态矩阵，b0dj和b0qj为控制矩阵，为输出矩阵，d0dj和d0qj为反馈系数，M为多机电力系统所包含风电场和光伏电站的数量，Iwdj和Iwqj分别为第j个风电场或光伏电站总输出电流的d轴和q轴分量，|Vwj|为第j个风电场或光伏电站内低压母线的电压幅值。
6.根据权利要求1所述方法，其特征在于，所述步骤4中电力系统机电振荡模式为：
<mrow>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&ap;</mo>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>d</mi>
<mi>j</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>q</mi>
<mi>j</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>d</mi>
<mrow>
<mn>0</mn>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>d</mi>
<mrow>
<mn>0</mn>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>54</mn>
<mo>)</mo>
</mrow>
</mrow>
残差为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>d</mi>
<mi>j</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>b</mi>
<mrow>
<mn>0</mn>
<mi>d</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<msub>
<mi>c</mi>
<mrow>
<mi>o</mi>
<mi>V</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>T</mi>
</msup>
<msub>
<mover>
<mi>r</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>q</mi>
<mi>j</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>b</mi>
<mrow>
<mn>0</mn>
<mi>q</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<msub>
<mi>c</mi>
<mrow>
<mi>o</mi>
<mi>V</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>T</mi>
</msup>
<msub>
<mover>
<mi>r</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>55</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，为开环系统的状态矩阵的特征值，和分别为对应的左右特征向量，为考虑了风电场和光伏电站动态后闭环阵的特征值估计值，Gdj(s)和Gqj(s)分别为第j个风电场或光伏电站的输出量ΔIwdj和ΔIwqj对输入量Δ|Vwj|的聚合传递函数，b0dj和b0qj为控制矩阵，为输出矩阵，d0dj和d0qj为反馈系数，M为多机电力系统所包含风电场和光伏电站的数量，Iwdj和Iwqj分别为第j个风电场或光伏电站总输出电流的d轴和q轴分量，|Vwj|为第j个风电场或光伏电站内低压母线的电压幅值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510513100.9A CN105140957B (zh) | 2015-08-20 | 2015-08-20 | 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510513100.9A CN105140957B (zh) | 2015-08-20 | 2015-08-20 | 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105140957A CN105140957A (zh) | 2015-12-09 |
| CN105140957B true CN105140957B (zh) | 2017-10-10 |
Family
ID=54726214
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510513100.9A Expired - Fee Related CN105140957B (zh) | 2015-08-20 | 2015-08-20 | 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105140957B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110417050A (zh) * | 2019-07-23 | 2019-11-05 | 华北电力大学 | 一种并联结构大型风电场等值计算方法 |
| CN110556871B (zh) * | 2019-10-17 | 2022-11-08 | 电子科技大学 | 基于结构保持方法的大规模光伏发电系统聚合等值方法 |
| CN113809779B (zh) * | 2021-10-14 | 2024-02-02 | 国网江苏省电力有限公司电力科学研究院 | 一种风电场接入柔直系统时稳定性检验方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE60318408T2 (de) * | 2003-06-21 | 2009-01-29 | Abb Research Ltd. | Ermittlung elektromechanischer Schwingungen in Stromversorgungssystemen |
| CN102136036B (zh) * | 2011-03-23 | 2013-01-16 | 天津大学 | 适于电力系统小扰动稳定分析的双馈风电场等值建模方法 |
| CN103036252A (zh) * | 2013-01-06 | 2013-04-10 | 辽宁省电力有限公司电力科学研究院 | 基于轨迹灵敏度法的风电场模型校核系统及方法 |
| CN103138259B (zh) * | 2013-01-09 | 2015-06-10 | 中国电力科学研究院 | 一种间歇性大规模风电接入电网安全稳定性分析方法 |
| CN103532167B (zh) * | 2013-09-29 | 2016-03-02 | 天津理工大学 | 基于fpga的并网光伏系统小干扰稳定性判断系统及方法 |
| CN103928943B (zh) * | 2014-04-21 | 2015-12-02 | 国家电网公司 | 双馈风电机组对互联电力系统低频振荡影响的研究方法 |
-
2015
- 2015-08-20 CN CN201510513100.9A patent/CN105140957B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN105140957A (zh) | 2015-12-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Nonlinear model predictive control for DFIG-based wind power generation | |
| Khedher et al. | Wind energy conversion system using DFIG controlled by backstepping and sliding mode strategies | |
| CN108923460B (zh) | 微网虚拟同步机多机并联动态一致响应的参数配置方法 | |
| Zhang et al. | An improved synchronverter model and its dynamic behaviour comparison with synchronous generator | |
| CN102346442B (zh) | 适于节点分析框架的微网暂态仿真并行计算方法 | |
| CN115313419B (zh) | 一种基于构网型变流器的直驱风电系统控制方法及系统 | |
| CN109256803A (zh) | 虚拟同步机孤岛运行小信号建模及参数灵敏度计算方法 | |
| CN105449699A (zh) | 双馈感应风电机组非线性分数阶自抗扰阻尼控制方法 | |
| CN102055205A (zh) | 一种基于电压不平衡数学模型的无刷双馈电机网侧变换器控制策略 | |
| Teng et al. | Voltage build-up analysis of self-excited induction generator with multi-timescale reduced-order model | |
| CN104993494A (zh) | 一种基于四象限电力电子变流器的电机模拟装置及方法 | |
| CN105140957B (zh) | 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 | |
| CN102801180A (zh) | 一种微电网并网等效模型的参数辨识方法 | |
| Altimania | Modeling of doubly-fed induction generators connected to distribution system based on eMEGASim® real-time digital simulator | |
| CN105162159B (zh) | 基于差分进化算法的光伏并网逆变器控制器参数辨识方法 | |
| Wang et al. | Variable speed drive connected motor load emulator for a multi-converter-based hardware testbed system | |
| Joshi et al. | Wireless speed control of an induction motor using PWM technique with GSM | |
| Wilson et al. | Nonlinear power flow control applications to conventional generator swing equations subject to variable generation | |
| CN106055817B (zh) | 用于电磁暂态分析的双馈异步发电机双工况模型建模方法 | |
| BARRA et al. | Sensorless speed and reactive power control of a DFIG-wind turbine | |
| Wang | Versatile three-phase power electronics converter based real-time load emulators | |
| CN102707161A (zh) | 一种基于动态转移阻抗的短路容量在线监测方法及其装置 | |
| Boussairi et al. | State feedback nonlinear control strategy for wind turbine system driven by permanent magnet synchronous generator for maximum power extraction and power factor correction | |
| Abulizi et al. | Study on Current Controller of Doubly-Fed Induction Generator with Variable Parameters | |
| CN116316893B (zh) | 基于鲁棒虚拟同步的分布式电源并网系统控制方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171010 Termination date: 20180820 |