CN113162113B - 确定含永磁风机电力***机电振荡源的节点模式能量方法 - Google Patents
确定含永磁风机电力***机电振荡源的节点模式能量方法 Download PDFInfo
- Publication number
- CN113162113B CN113162113B CN202110416760.0A CN202110416760A CN113162113B CN 113162113 B CN113162113 B CN 113162113B CN 202110416760 A CN202110416760 A CN 202110416760A CN 113162113 B CN113162113 B CN 113162113B
- Authority
- CN
- China
- Prior art keywords
- generator
- power
- permanent magnet
- energy
- node mode
- 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.)
- Active
Links
Images
Classifications
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
-
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明的一种确定含永磁风机电力***机电振荡源的节点模式能量方法,其特点是,通过采用建立电力***的小干扰稳定分析微分代数方程、对***状态方程在小扰动下求解、对发电机电磁功率计算、对发电机电磁功率变化量ΔPE进行线性化计算、对小扰动下同步发电机内节点模式能量计算和对小扰动下同步发电机内节点模式能量角度计算步骤,对含高渗透率永磁风机电力***的振荡机理进行分析,能够准确确定***中的振荡源,通过确定振荡源探究永磁风机对***低频振荡的影响。具有方法科学合理,适用性强,效果佳等优点。
Description
技术领域
本发明涉及电力***规划与运行领域,是一种确定含永磁风机电力***机电振荡源的节点模式能量方法,特别适用于对含高渗透率永磁风机并网电力***进行低频振荡机理的分析。
背景技术
随着互联电网的急速发展,越来越多的新能源发电装置接入电网。风电作为一种高效、便捷的可再生能源,正逐步的扩大其在电网中的渗透率。其中永磁风机(Permanentmagnet synchronous generator,PMSG)以结构简单、价格低廉、维修便捷、可实现有功功率和无功功率完全解耦的特点被广泛的应用在风力发电中。低频振荡是电力***面临的重大威胁之一,它的振荡频率低、周期长,当电网发生低频振荡事故时,可能会导致大规模的停电,严重的威胁着电力***的稳定。而永磁风机的大规模接入会加剧低频振荡问题对电网的威胁,此时电力***的低频振荡问题变得更加复杂。故对于高渗透率永磁风机并网电力***来说,及时找到振荡源十分重要。迄今未见有关确定含永磁风机电力***机电振荡源的节点模式能量方法的文献报导和实际使用。
发明内容
本发明所要解决的问题是:针对高渗透率永磁风机并网的电力***存在的问题,提供一种科学合理,适用性强,效果佳的确定含永磁风机电力***机电振荡源的节点模式能量方法,用所述的方法对含高渗透率永磁风机电力***的振荡机理进行分析,能够准确确定振荡源。
采用的技术方案是:一种确定含永磁风机电力***机电振荡源的节点模式能量方法,其特征是,它包括以下步骤:
1)建立电力***的小干扰稳定分析微分代数方程
其中:ΔM—n维状态变量;k1、k2、k3、k4—系数矩阵;Δu—节点电压的幅值和角度增量;
2)对***状态方程在小扰动下求解
其中:Δmk(t)—***状态方程在小扰动下的解;β1,β2,…,βz—特征值;z—特征值数量;k—节点数量;c1,c2,…,c3—系数;
3)对发电机电磁功率计算
PE=Eq′uksinδ/x (3)
其中:PE—发电机电磁功率;Eq′—发电机q轴暂态电势;uk—发电机内节点k电压;δ—发电机功角;x—发电机d轴暂态电抗;
4)对发电机电磁功率变化量ΔPE进行线性化计算
ΔPE=Eq′(Δu1sinδ0+Δδu10cosδ0)/x (4)
其中:ΔPE—发电机电磁功率变化量;Eq′—发电机q轴暂态电势;Δu1—发电机内节点电压变化量;u10—发电机内节点电压初值;Δδ—发电机功角变化量;δ0—发电机功角初值;
5)对小扰动下同步发电机内节点模式能量计算
ΔLG=∫ΔPE·Δωdt (5)
其中:ΔLG—小干扰下的同步发电机内节点模式能量;ΔPE—小干扰下发电机电磁功率增量;Δω—发电机转子角速度增量;
6)对小扰动下同步发电机内节点模式能量角度计算
其中:γG—小扰动下同步发电机内节点模式能量角度;ωG—发电机转子角速度;A—代数系数。
本发明的一种确定含永磁风机电力***机电振荡源的节点模式能量方法,是基于高渗透率永磁风机的接入电力***后,与传统的电力***相比更不稳定、更易受到小干扰的影响,同时高渗透率永磁风机接入电力***后的振荡特性也更加复杂的实际情况,通过采用建立电力***的小干扰稳定分析微分代数方程、对***状态方程在小扰动下求解、对发电机电磁功率计算、对发电机电磁功率变化量ΔPE进行线性化计算、对小扰动下同步发电机内节点模式能量计算和对小扰动下同步发电机内节点模式能量角度计算步骤,对含高渗透率永磁风机电力***的振荡机理进行分析,能够准确确定***中的振荡源,通过确定振荡源探究永磁风机对***低频振荡的影响;并且通过二机***和四机两区***的仿真算例验证所提出方法的有效性;***能量是由多种模式能量耦合而成,计算各模式下发电机的内节点模式能量及能量角度能够确定各模式下主要参与振荡的同步机。具有方法科学合理,适用性强,效果佳等优点。
附图说明
图1为本发明所使用的二机永磁风机并网***示意图;
图2为本发明所使用的二机永磁风机并网***区间振荡同步机内节点模式势能示意图;
图3为本发明所使用的二机永磁风机并网***区间振荡同步机内节点模式势能角度示意图;
图4为本发明所使用的四机两区永磁风机并网***示意图;
图5为本发明所使用的四机两区永磁风机并网***区间振荡同步机内节点模式势能示意图;
图6为本发明所使用的四机两区永磁风机并网***区间振荡同步机内节点模式势能角度示意图。
具体实施方式
以下结合附图和具体实施方式对本发明作进一步的说明。
本发明的一种基于模式能量确定永磁风机并网电力***机电振荡源的方法,包括以下步骤:
1)建立电力***的小干扰稳定分析微分代数方程
其中:ΔM—n维状态变量;k1、k2、k3、k4—系数矩阵;Δu—节点电压的幅值和角度增量;
2)对***状态方程在小扰动下求解
其中:Δmk(t)—***状态方程在小扰动下的解;β1,β2,…,βz—特征值;z—特征值数量;k—节点数量;c1,c2,…,c3—系数;
3)对发电机电磁功率计算
PE=Eq′uksinδ/x (3)
其中:PE—发电机电磁功率;Eq′—发电机q轴暂态电势;uk—发电机内节点k电压;δ—发电机功角;x—发电机d轴暂态电抗;
4)对发电机电磁功率变化量ΔPE进行线性化计算
ΔPE=Eq′(Δu1sinδ0+Δδu10cosδ0)/x (4)
其中:ΔPE—发电机电磁功率变化量;Eq′—发电机q轴暂态电势;Δu1—发电机内节点电压变化量;u10—发电机内节点电压初值;Δδ—发电机功角变化量;δ0—发电机功角初值;
5)对小扰动下同步发电机内节点模式能量计算
ΔLG=∫ΔPE·Δωdt (5)
其中:ΔLG—小干扰下的同步发电机内节点模式能量;ΔPE—小干扰下发电机电磁功率增量;Δω—发电机转子角速度增量;
6)对小扰动下同步发电机内节点模式能量角度计算
其中:γG—小扰动下同步发电机内节点模式能量角度;ωG—发电机转子角速度;A—代数系数。
本发明对永磁风机并入二机***为例,此时***中风电渗透率达40%,***如图1所示。在母线L3-4处设置三相短路故障,对该***进行传统小干扰计算,得到结果如表1所示。其中主导振荡模式时频率为1.6412Hz的振荡模式,此时最大的参与因子为发电机G1,但与G1产生相对振荡的机组不明确,对该***进行发电机内节点模式能量及能量角度计算,结果如图2、图3所示。根据计算结果进行振荡源识别,图2中可见G2的内节点模式能量振荡方向与G1相反,幅值相近,说明G2是除G1外参与程度最大的机组。图3中可见G1、G2的内节点模式能量角度呈反向,说明该模式下发生相对振荡的机组是{G1-G2}。
再选用四机两区算例,在母线2和母线4处同时等容量将同步机替换为永磁风机,此时***风电渗透率为43%,如图4所示。在母线7处设置0.01秒三相短路故障,对该***进行传统小干扰计算,得到结果如表2所示。其中主导振荡模式时频率为0.72305Hz的振荡模式,此时最大的参与因子为发电机G1,但与G1产生相对振荡的机组不明确,对该***进行发电机内节点模式能量及能量角度计算,结果如图5、图6所示。根据计算结果进行振荡源识别,图5中可见G2的内节点模式能量振荡方向与G1相反,幅值相近,说明G2是除G1外参与程度最大的机组。图6中可见G1、G2的内节点模式能量角度呈反向,说明该模式下发生相对振荡的机组是{G1-G2}。
表1二机***接入PMSG小干扰计算结果
表2四机两区***并入PMSG小干扰计算结果
上述实施例仅为清楚地说明本发明技术方案所作的举例,而并非是对本发明的实施方式的限定。在不改变本发明基本构思和实质的情况下,任何其它等同技术特征的变换或修改,都应属于本发明权利要求的保护范围。
Claims (1)
1.一种确定含永磁风机电力***机电振荡源的节点模式能量方法,其特征是,它包括以下步骤:
1)建立电力***的小干扰稳定分析微分代数方程作为***状态方程
其中:ΔM—n维状态变量;k1、k2、k3、k4—系数矩阵;Δu—节点电压的幅值和角度增量;
2)对***状态方程在小扰动下求解
其中:Δmk(t)—***状态方程在小扰动下的解;β1,β2,…,βz—特征值;z—特征值数量;k—节点数量;c1,c2,…,cz—系数;
3)对发电机电磁功率计算
PE=E′quksinδ/x (3)
其中:PE—发电机电磁功率;E′q—发电机q轴暂态电势;uk—发电机内节点k电压;δ—发电机功角;x—发电机d轴暂态电抗;
4)对发电机电磁功率变化量ΔPE进行线性化计算
ΔPE=E′q(Δu1sinδ0+Δδu10cosδ0)/x (4)
其中:ΔPE—发电机电磁功率变化量;Eq′—发电机q轴暂态电势;Δu1—发电机内节点电压变化量;u10—发电机内节点电压初值;Δδ—发电机功角变化量;δ0—发电机功角初值;
5)对小扰动下同步发电机内节点模式能量计算
ΔLG=∫ΔPE·Δωdt (5)
其中:ΔLG—小干扰下的同步发电机内节点模式能量;Δω—发电机转子角速度增量;
6)对小扰动下同步发电机内节点模式能量角度计算
其中:γG—小扰动下同步发电机内节点模式能量角度;ωG—发电机转子角速度;A—代数系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110416760.0A CN113162113B (zh) | 2021-04-19 | 2021-04-19 | 确定含永磁风机电力***机电振荡源的节点模式能量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110416760.0A CN113162113B (zh) | 2021-04-19 | 2021-04-19 | 确定含永磁风机电力***机电振荡源的节点模式能量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113162113A CN113162113A (zh) | 2021-07-23 |
CN113162113B true CN113162113B (zh) | 2022-10-28 |
Family
ID=76868393
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110416760.0A Active CN113162113B (zh) | 2021-04-19 | 2021-04-19 | 确定含永磁风机电力***机电振荡源的节点模式能量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113162113B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178535A (zh) * | 2013-02-27 | 2013-06-26 | 中国电力科学研究院 | 一种兼顾两类机制的电力***低频振荡在线防控方法 |
CN106526384A (zh) * | 2016-12-01 | 2017-03-22 | 西南交通大学 | 一种用于大规模电力***的振荡源定位方法 |
CN109217362A (zh) * | 2018-08-22 | 2019-01-15 | 华北电力大学 | 一种双馈风机并网***低频振荡扰动源定位***及方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105243214B (zh) * | 2015-10-12 | 2018-12-07 | 国家电网公司 | 一种互联电网强迫功率振荡节点敏感性评估方法 |
CN109713686B (zh) * | 2018-12-10 | 2020-12-22 | 国家电网有限公司 | 一种高渗透率风电电力***风电支路模式能量解析方法 |
CN110795840A (zh) * | 2019-10-22 | 2020-02-14 | 海南电网有限责任公司电力科学研究院 | 基于dmd的电力***主导振荡模式及参数识别方法 |
-
2021
- 2021-04-19 CN CN202110416760.0A patent/CN113162113B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178535A (zh) * | 2013-02-27 | 2013-06-26 | 中国电力科学研究院 | 一种兼顾两类机制的电力***低频振荡在线防控方法 |
CN106526384A (zh) * | 2016-12-01 | 2017-03-22 | 西南交通大学 | 一种用于大规模电力***的振荡源定位方法 |
CN109217362A (zh) * | 2018-08-22 | 2019-01-15 | 华北电力大学 | 一种双馈风机并网***低频振荡扰动源定位***及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113162113A (zh) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ju et al. | Analytical assessment for transient stability under stochastic continuous disturbances | |
Zhang et al. | Responses and stability of power system under small Gauss type random excitation | |
CN104362622B (zh) | 一种带小水电机组地区电网的负荷模型和建模方法 | |
CN113034002B (zh) | 一种地磁暴对电力***小扰动电压稳定性的分析方法 | |
CN103267897A (zh) | 一种基于反Park变换的三相锁相环 | |
Morgan et al. | Wind farm dynamic models assessment under weak grid conditions | |
Benbouhenni et al. | Backstepping control for multi-rotor wind power systems | |
Ren et al. | Component damping evaluation in sub‐synchronous oscillation based on transient energy flow method | |
Feng et al. | Influences of DC bus voltage dynamics in modulation algorithm on power oscillations in PMSG-based wind farms | |
Deng et al. | Analysis of electromagnetic transient characteristics of doubly-fed induction generator under grid voltage swell | |
Wang et al. | Fault diagnosis based on current signature analysis for stator winding of Doubly Fed Induction Generator in wind turbine | |
CN113162113B (zh) | 确定含永磁风机电力***机电振荡源的节点模式能量方法 | |
He et al. | Investigation of the effects of various types of wind turbine generators on power-system stability | |
CN105302999B (zh) | 发现汽轮发电机内负序分量对转矩与谐波影响规律的方法 | |
Li et al. | DFIG-based wind farm electromagnetic dynamic model and impact on protection relay of transmission network | |
Yu et al. | Simultaneous suppression of torque ripple and flexible load vibration for PMSM under stator current vector orientation | |
CN101340121A (zh) | 一种发电机内部故障暂态仿真简化模型及其建模方法和应用 | |
El Azzaoui et al. | Comparative study of the sliding mode and backstepping control in power control of a doubly fed induction generator | |
CN114465280A (zh) | 一种新能源并网***动态等效建模方法 | |
CN113612238A (zh) | 一种分析风电并网引发火电机组轴系小干扰振荡的改进阻尼转矩方法 | |
Yan et al. | Impact of DFIG-based WTs on subsynchronous oscillation damping of SG based on the net electrical damping method | |
Cui et al. | Synchronous Condenser Parameter Estimation Based on Lifting Wavelet Transform and Particle Swarm Optimization | |
Zhang et al. | The Load Equivalence Method for Dynamic Frequency Analysis of large Power Grid | |
Yang et al. | European Journal of Electrical Engineering | |
Zhou et al. | Calculation and analysis of unbalanced magnetic pulls of different stator winding setups in static eccentric induction motor |
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 |