CN109713686B - 一种高渗透率风电电力系统风电支路模式能量解析方法 - Google Patents
一种高渗透率风电电力系统风电支路模式能量解析方法 Download PDFInfo
- Publication number
- CN109713686B CN109713686B CN201811502670.8A CN201811502670A CN109713686B CN 109713686 B CN109713686 B CN 109713686B CN 201811502670 A CN201811502670 A CN 201811502670A CN 109713686 B CN109713686 B CN 109713686B
- Authority
- CN
- China
- Prior art keywords
- wind power
- mode
- potential energy
- power branch
- branch 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.)
- Expired - Fee Related
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 18
- 238000005381 potential energy Methods 0.000 claims abstract description 64
- 230000035945 sensitivity Effects 0.000 claims abstract description 24
- 238000000034 method Methods 0.000 claims abstract description 14
- 230000010355 oscillation Effects 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 5
- 230000000007 visual effect Effects 0.000 abstract description 4
- 238000012423 maintenance Methods 0.000 abstract description 2
- 230000035699 permeability Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
Images
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/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
一种高渗透率风电电力系统风电支路模式能量解析方法,属于电力系统运维技术领域,本发明通过对风电电力系统风电支路模式能量建模,构建了简单直观、物理意义清晰的风电电力系统风电支路模式势能函数,可以以此分析风电支路模式势能的构成及变化规律,研究支路模式势能变化对系统小干扰稳定性的影响;还可以分析风电支路模式势能对状态变量的灵敏度,分析不同状态变量的敏感程度,从而通过调整状态变量来控制和抑制低频振荡;构造仿真算例,将扰动后量测计算的风电支路模式势能和频域模式势能对比,验证了所提方法的有效性,具有较高的实际应用价值。
Description
技术领域
本发明属于电力系统运维技术领域,特别是涉及到一种应用于高渗透率风电电力系统低频振荡定量分析和抑制的能量解析。
背景技术
高比例可再生能源并网将成为未来电力系统的重要特征,高渗透率风电出力的随机性和不确定性,加剧了电力系统小干扰稳定的复杂程度,使得电力系统安全稳定运行面临着更大的挑战。
能量函数能描述电力系统能量变化关系,已有学者将能量函数引入到低频振荡问题,从能量角度分析低频振荡的作用机理,有效识别网络中的薄弱环节、实现振荡度评估,为解决低频振荡的分析与抑制提供切实可行的新思路,对研究小干扰稳定性具有重要意义。但对高渗透率风电电力系统能量角度分析并未深入。因此现有技术当中亟需要一种新型的技术方案来解决这一问题。
发明内容
本发明所要解决的技术问题是:提供一种高渗透率风电电力系统风电支路模式能量解析方法,通过对风电电力系统风电支路模式能量建模,构建了简单直观、物理意义清晰的风电电力系统风电支路模式势能函数,可以以此分析风电支路模式势能的构成及变化规律,研究支路模式势能变化对系统小干扰稳定性的影响;还可以分析风电支路模式势能对状态变量的灵敏度,分析不同状态变量的敏感程度,从而通过调整状态变量来控制和抑制低频振荡;构造仿真算例,将扰动后量测计算的风电支路模式势能和频域模式势能对比,验证了所提方法的有效性,具有较高的实际应用价值。
一种高渗透率风电电力系统风电支路模式能量解析方法,其特征是:包括以下步骤,且以下步骤顺次进行,
步骤一、获得风电支路模式势能表达式
其中:△P为风机功率增量,△θ为风电支路母线电压相角增量;
风电支路有功表达式为
DFIG为双馈风机,xm为DFIG的电磁电抗,xs为DFIG的定子电抗,iqr为DFIG的q轴转子电流,V为节点电压;
式中,ci=ΨiΔx(0),Δx(0)为状态变量增量的初值,λi(i=1,…,7)为特征值;Ψi为左特征向量;“Φ7i为右特征向量;
通过网络变量线性化增量函数获得
式中,F=-D-1C=[F1,F2,F3,F4,F5,F6,F7]T
步骤二、获得风电支路模式势能总表达式
根据步骤一获得的各状态变量各自对应的特征值,获得风电支路模式势能总表达式,
风电支路势能由自相关模式和互相关模式构成,自相关模式由机电自相关模式和风电状态量自相关模式构成,互相关模式由机电互相关模式、风电状态量互相关模式和机电-风电状态量互相关模式构成,
其中:△Vpz为风电支路模式势能自相关部分,
ΔVph为风电支路模式势能互相关部分;
式中Φ7i、Φki、Φkq、Φkp、Φ7p、F1k、F2k、Fqk,下角标第一个符号代表行,第二个符号代表列;
步骤三、分析风电支路模式势能对状态变量的灵敏度
对风机出口支路模式势能状态变量灵敏度的状态变量求偏导数,表达式为:
式中:y1~y16为常系数,△δi为发电机的功角增量,△ωm为DFIG的转子角速度增量,△idr为DFIG的d轴转子电流增量;
获得风电支路模式势能对状态变量的灵敏度与发电机功角、风机转子角速度、风机转子d轴和q轴电流有关,与发电机转子角速度、风速和桨距角无关,并通过调整状态变量控制和抑制低频振荡。
所述步骤一中的系统小干扰下状态方程解的一般形式为
式中:ci=ΨiΔx(0),“λi(i=1,,7)为特征值,Ψi为左特征向量,Φi为右特征向量。
所述步骤一中的网络变量线性化增量函数为
式中:F=-D-1C=[F1,F2,F3,F4,F5,F6,F7]T,
式中A、B、C、D是系数矩阵。
通过上述设计方案,本发明可以带来如下有益效果:一种高渗透率风电电力系统风电支路模式能量解析方法,通过对风电电力系统风电支路模式能量建模,构建了简单直观、物理意义清晰的风电电力系统风电支路模式势能函数,可以以此分析风电支路模式势能的构成及变化规律,研究支路模式势能变化对系统小干扰稳定性的影响;还可以分析风电支路模式势能对状态变量的灵敏度,分析不同状态变量的敏感程度,从而通过调整状态变量来控制和抑制低频振荡;构造仿真算例,将扰动后量测计算的风电支路模式势能和频域模式势能对比,验证了所提方法的有效性,具有较高的实际应用价值。
附图说明
以下结合附图和具体实施方式对本发明作进一步的说明:
图1为本发明发电机和双馈风机接入无穷大系统示意图。
图2为本发明风电渗透率为20%时风电支路总模式势能分析示意图。
图3为本发明风电渗透率为20%时风电支路自相关模式势能分析示意图。
图4为本发明风电渗透率为20%时风电支路互相关模式势能分析示意图。
图5为本发明风电渗透率为20%时风电支路机电风电状态量分量分析示意图。
图6为本发明不同风电渗透率时风电支路主导模式分析示意图。
图7为本发明20%风电渗透率时风电支路模式势能对状态变量灵敏度示意图。
图8为本发明80%风电渗透率时风电支路模式势能对状态变量灵敏度示意图。
图9为本发明理想的风电支路机电自相关模式势能示意图。
图10为本发明量测的风电支路机电自相关模式势能示意图。
图中1-发电机、2-双馈风机、3-变压器、4-输电线路、5-无穷大系统母线。
具体实施方式
一种高渗透率风电电力系统风电支路模式能量解析方法,首先,对含有风电电力系统小干扰稳定性进行理论分析,
同步发电机采用经典二阶模型,忽略励磁系统和调速动态,双馈风机DFIG采用五阶模型,同步发电机和双馈风机构成微分方程为
其中:δ为发电机的功角,ω为发电机的转速,νw为风速,ωm为DFIG的转子角速度,θp为DFIG的桨距角,idr为DFIG的d轴转子电流,iqr为DFIG的q轴转子电流,M为发电机的惯性时间常数,Pm为发电机的机械功率,PE为发电机的电磁功率,τ为时间常数,Tm为DFIG的机械转矩,Te为DFIG的电磁转矩,Hm为DFIG的转子惯量,Kp为桨距控制增益,φ为随机相角,Tp为DFIG的桨距角时间常数,KV为DFIG的电压控制增益,xm为DFIG的电磁电抗,xs为DFIG的定子电抗,Pw*(ωm)为功率转速特性;
将发电机和双馈风机构成微分方程简化成矩阵形式为
式中:xo=[νw,ωm,θp,idr,iqr]T为风机状态向量,
u=[θi,Vi,vf,pm,p,q,ωs,pwa,vref]T为输入向量;
系统小干扰下状态方程解的一般形式为
式中:ci=ΨiΔx(0),“λi(i=1,,7)为特征值,Ψi为左特征向量,Φi为右特征向量;
网络变量线性化增量函数为
式中:F=-D-1C=[F1,F2,F3,F4,F5,F6,F7]T
风电支路模式势能为
其中:△P为风机功率增量,△θ为风电支路母线电压相角增量;
风电支路有功表达式为
线性化得:
由(3)式得:
由(4)式得:
风电支路模式势能总表达式为
风电支路势能由自相关模式和互相关模式构成,自相关模式由机电自相关模式和风电状态量自相关模式构成,互相关模式由机电互相关模式、风电状态量互相关模式和机电-风电状态量互相关模式构成,
其中:ΔVpz为风电支路模式势能自相关部分,ΔVph为风电支路模式势能互相关部分;
风机出口支路模式势能对状态变量灵敏度即对状态变量求偏导数,表达式为
式中:y1~y16为常系数;
由(12)式可知风电支路模式势能对状态变量的灵敏度与发电机功角、风机转子角速度、风机转子d轴和q轴电流有关,与发电机转子角速度、风速和桨距角无关。分析不同状态变量的敏感程度,从而通过调整状态变量来控制和抑制低频振荡。
以发电机和双馈风机接入无穷大系统为对象进行详细分析,用电力系统分析软件PSAT(Power System Analysis Toolbox)仿真分析如图1所示,对20%渗透率风电时风电支路模式势能进行分析,如图2~图5所示。机电自相关模式是风电支路模式势能的主导模式,因此对不同渗透率风电机电自相关模式进行分析,如图6所示.再分析20%和80%风电渗透率时风电支路模式势能对各状态变量的灵敏度,如图7和图8所示,已知风机转子角速度与机电尺度下的转速环有关,风机转子d轴和q轴电流与非机电尺度下的电流环有关。本算例风机转子角速度灵敏度大于风机转子d轴和q轴电流灵敏度,即机电尺度下的转速环灵敏度大于非机电尺度下的电流环灵敏度,随着风电渗透率的增加,机电尺度下的转速环始终占主导地位。考虑风电渗透率为20%时,当母线3设置0.05秒小扰动,提取扰动后的风电支路有功功率和母线电压角度信号,计算风电支路势能,进行模态辨识提取机电自相关模式势能。将扰动后计算的风机模式势能和频域模式势能对比,如图9和图10所示两者基本一致,从而证明所提方法的有效性。
本发明通过对风电电力系统风电支路模式能量建模,构建了简单直观、物理意义清晰的风电电力系统风电支路模式势能函数,可以以此分析风电支路模式势能的构成及变化规律,研究支路模式势能变化对系统小干扰稳定性的影响;还可以分析风电支路模式势能对状态变量的灵敏度,分析不同状态变量的敏感程度,从而通过调整状态变量来控制和抑制低频振荡;构造仿真算例,将扰动后量测计算的风电支路模式势能和频域模式势能对比,验证了所提方法的有效性,具有较高的实际应用价值。
Claims (3)
1.一种高渗透率风电电力系统风电支路模式能量解析方法,其特征是:包括以下步骤,且以下步骤顺次进行,
步骤一、获得风电支路模式势能表达式
其中:△P为风机功率增量,△θ为风电支路母线电压相角增量;
风电支路有功表达式为
DFIG为双馈风机,xm为DFIG的电磁电抗,xs为DFIG的定子电抗,iqr为DFIG的q轴转子电流,V为节点电压;
式中,ci=ΨiΔx(0),Δx(0)为状态变量增量的初值,λi(i=1,…,7)为特征值;Ψi为左特征向量;Φ7i为右特征向量;
通过网络变量线性化增量函数获得
式中,F=-D-1C=[F1,F2,F3,F4,F5,F6,F7]T
步骤二、获得风电支路模式势能总表达式
根据步骤一获得的各状态变量各自对应的特征值,获得风电支路模式势能总表达式,
风电支路势能由自相关模式和互相关模式构成,自相关模式由机电自相关模式和风电状态量自相关模式构成,互相关模式由机电互相关模式、风电状态量互相关模式和机电-风电状态量互相关模式构成,
其中:△Vpz为风电支路模式势能自相关部分,
ΔVph为风电支路模式势能互相关部分;
式中,Φ7i、Φki、Φkq、Φkp、Φ7p、F1k、F2k、Fqk下角标第一个符号代表行,第二个符号代表列;
步骤三、分析风电支路模式势能对状态变量的灵敏度
对风机出口支路模式势能状态变量灵敏度的状态变量求偏导数,表达式为:
式中:y1~y16为常系数,△δi为发电机的功角增量,△ωm为DFIG的转子角速度增量,△idr为DFIG的d轴转子电流增量;
获得风电支路模式势能对状态变量的灵敏度与发电机功角、风机转子角速度、风机转子d轴和q轴电流有关,与发电机转子角速度、风速和桨距角无关,并通过调整状态变量控制和抑制低频振荡。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811502670.8A CN109713686B (zh) | 2018-12-10 | 2018-12-10 | 一种高渗透率风电电力系统风电支路模式能量解析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811502670.8A CN109713686B (zh) | 2018-12-10 | 2018-12-10 | 一种高渗透率风电电力系统风电支路模式能量解析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109713686A CN109713686A (zh) | 2019-05-03 |
CN109713686B true CN109713686B (zh) | 2020-12-22 |
Family
ID=66255564
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811502670.8A Expired - Fee Related CN109713686B (zh) | 2018-12-10 | 2018-12-10 | 一种高渗透率风电电力系统风电支路模式能量解析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109713686B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110707694B (zh) * | 2019-10-29 | 2022-09-20 | 东北电力大学 | 风电并网电力系统通过网络模式能量识别振荡路径的方法 |
CN110707730B (zh) * | 2019-11-01 | 2022-06-28 | 东北电力大学 | 一种双馈风机并网电力系统低频振荡调控方法 |
CN113162113B (zh) * | 2021-04-19 | 2022-10-28 | 东北电力大学 | 确定含永磁风机电力系统机电振荡源的节点模式能量方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102185325A (zh) * | 2011-04-25 | 2011-09-14 | 东北电力大学 | 基于网络量测信息的电力系统暂态稳定性定量评价方法 |
CN102354997A (zh) * | 2011-10-22 | 2012-02-15 | 东北电力大学 | 电力系统低频振荡关键支路及振荡割集辨识方法 |
CN102411118B (zh) * | 2011-12-01 | 2013-06-26 | 武汉华中电力电网技术有限公司 | 一种区域互联电网强迫功率振荡扰动源位置判断方法 |
CN103178535B (zh) * | 2013-02-27 | 2015-04-29 | 中国电力科学研究院 | 一种兼顾两类机制的电力系统低频振荡在线防控方法 |
CN103475010B (zh) * | 2013-08-09 | 2015-05-13 | 国家电网公司 | 一种冲击扰动后互联电网的暂态稳定控制方法 |
CN104967135B (zh) * | 2015-07-08 | 2017-03-01 | 上海电力学院 | 含风电的多域时滞互联电力系统滑模负荷频率控制方法 |
GB2557220B (en) * | 2016-11-30 | 2019-05-15 | Reactive Tech Limited | Determining a characteristic of an inertial contribution to an electric power grid |
CN108306287A (zh) * | 2018-02-09 | 2018-07-20 | 清华大学 | 一种含风电电力系统的小干扰失稳风险识别方法及系统 |
CN108964125A (zh) * | 2018-07-19 | 2018-12-07 | 长沙理工大学 | 提高主网功角安全鲁棒性的大型风电接入点选择方法 |
-
2018
- 2018-12-10 CN CN201811502670.8A patent/CN109713686B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN109713686A (zh) | 2019-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109713686B (zh) | 一种高渗透率风电电力系统风电支路模式能量解析方法 | |
CN105156271B (zh) | 受扰风力发电机系统的抗干扰最优功率跟踪控制方法 | |
Cidras et al. | A linear dynamic model for asynchronous wind turbines with mechanical fluctuations | |
CN103336909A (zh) | 一种风电接入电网的低频振荡辨识方法 | |
CN107451364A (zh) | 一种dfig风电场等值参数的辨识方法 | |
De Kooning et al. | Shaft speed ripples in wind turbines caused by tower shadow and wind shear | |
CN110212575A (zh) | 双馈风机单机无穷大系统的小干扰稳定性模态分析方法 | |
CN104218864A (zh) | 一种双馈风电机组转子侧控制器参数辨识的频域方法 | |
CN116402017B (zh) | 一种用于电力系统频率动态分析的双馈风机简化建模方法 | |
Benaaouinate et al. | Development of a useful wind turbine emulator based on permanent magnet DC motor | |
Corradini et al. | An aerodynamic torque observer for the robust control of variable-speed wind turbines | |
CN103485978A (zh) | 一种电磁转矩补偿实现最大风能快速平稳跟踪的控制方法 | |
Syskakis et al. | MPPT for small wind turbines: Zero-oscillation sensorless strategy | |
Pan et al. | Parameter estimation of drive system in a fixed‐speed wind turbine by utilising turbulence excitations | |
Mishra et al. | Design of non-linear controller for ocean wave energy plant | |
Wu et al. | Backstepping terminal sliding mode control of DFIG for maximal wind energy captured | |
Tshiloz et al. | Sensorless speed estimation in wound rotor induction machines by spectral search of the stator phase power signal | |
Amêndola et al. | Fuzzy-logic control system of a variable-speed variable-pitch wind-turbine and a double-fed induction generator | |
CN116722533A (zh) | 一种基于随机风电预测区间的参数化暂态稳定预防控制求解方法 | |
Ortega et al. | A globally convergent wind speed estimator for windmill systems | |
Xing et al. | On-line Identification of Equivalent Inertia for DFIG Wind Turbines Based on Extended Kalman Filters | |
Minka et al. | Power Control of a DFIG Driving by Wind Turbine: Comparison Study Between ADRC and PI Controller | |
Benaaouinate et al. | Emulation of Wind Turbine for Standalone Wind Energy Conversion Systems | |
CN104343628A (zh) | 一种含有死区补偿器的风电机组变桨距控制方法 | |
Budisan et al. | Rotation speed and wind speed indirect measurement methods for the control of windmills with fixed blades turbine |
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 | ||
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: 20201222 |