CN110829465B - 计及多调速器死区的电力系统超低频振荡风险评估方法 - Google Patents
计及多调速器死区的电力系统超低频振荡风险评估方法 Download PDFInfo
- Publication number
- CN110829465B CN110829465B CN201911201185.1A CN201911201185A CN110829465B CN 110829465 B CN110829465 B CN 110829465B CN 201911201185 A CN201911201185 A CN 201911201185A CN 110829465 B CN110829465 B CN 110829465B
- Authority
- CN
- China
- Prior art keywords
- frequency oscillation
- power system
- dead zone
- speed regulator
- generator
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
Abstract
本发明公开了计及多调速器死区的电力系统超低频振荡风险评估方法:采集联网发电机及其调速器模型和运行参数;确定等值发电机的传递函数,根据调速器模型确定各调速器线性部分的传递函数,计算各调速器死区环节的描述函数;结合各调速器及其死区环节,计算各调速器扩展描述函数的幅值和相角,将多个扩展描述函数相加组合成系统的扩展描述函数;建立统一频率模型;在统一频率模型的基础上,计算产生超低频频振荡的临界振幅,通过评估电力系统超低频振荡的风险。本发明用以解决现有技术中无法分析包含不同死区的多机系统的频率振荡的问题,实现对含有多个调速器死区环节且死区大小不同的电力系统超低频振荡风险进行评估的目的。
Description
技术领域
本发明涉及电力系统超低频振荡领域,具体涉及计及多调速器死区的电力系统超低频振荡风险评估方法。
背景技术
电力系统的超低频振荡是近年来最受关注的频率稳定问题之一。调速器死区对于频率振荡的影响不可忽略,调速器死区在实际电网中广泛存在并且各不相同,但是目前的研究方法仅能分析单一死区环节,还没有能够分析包含不同死区的多机系统的频率振荡的方法。
发明内容
本发明的目的在于提供计及多调速器死区的电力系统超低频振荡风险评估方法,以解决现有技术中无法分析包含不同死区的多机系统的频率振荡的问题,实现对含有多个调速器死区环节且死区大小不同的电力系统超低频振荡风险进行评估的目的。
本发明通过下述技术方案实现:
计及多调速器死区的电力系统超低频振荡风险评估方法,包括:
S1.采集联网发电机及其调速器模型和运行参数;
S2.确定等值发电机的传递函数Ggen,根据调速器模型确定各调速器线性部分的传递函数Ggovi,计算各调速器死区环节的描述函数NLi;
S3.结合各调速器Ggovi及其死区环节NLi,计算各调速器扩展描述函数的幅值|Gei|和相角∠Gei,将多个扩展描述函数相加组合成系统的扩展描述函数Ge;
S4.利用Ge和Ggen,建立统一频率模型;在统一频率模型的基础上,计算产生超低频频振荡的临界振幅Ac,通过Ac评估电力系统超低频振荡的风险。
进一步的,步骤S4中,通过Ac评估电力系统超低频振荡的风险的方法为:临界振幅Ac越高、能够使系统产生振荡的最小扰动幅值越大、系统的稳定性越高。
进一步的,步骤S4中,产生超低频频振荡的临界振幅Ac,根据极限环产生的条件进行计算。
进一步的,所述极限环产生的条件为:
|Ge|·|Ggen|=1,∠Ge+∠Ggen=-π/2。
进一步的,步骤S2中,等值发电机的传递函数Ggen根据发电机的转动惯量、阻尼系数和负荷频率调节效应进行确定。
其中,s为拉普拉斯算子,Ds为发电机阻尼系数与负荷频率调节效应系数之和,TJ为等值发电机的转动惯量;
进一步的,步骤S2中,调速器死区环节的描述函数NLi通过下式确定:
其中,b是死区的宽度;A0为待求参数,表示振荡的幅值大小。
进一步的,步骤S3中,扩展描述函数的幅值|Gei|和相角∠Gei通过下式确定:
|Gei(A0,ω)|=|NLi||Ggovi|
∠Gei(A0,ω)=∠NLi+∠Ggovi。
进一步的,步骤S2中,各调速器死区环节的描述函数NLi通过描述函数法计算得到。
本发明与现有技术相比,具有如下的优点和有益效果:
本发明计及多调速器死区的电力系统超低频振荡风险评估方法,能够分析含有不同调速器死区的多机电力系统的超低频振荡问题,计算结果能够有效反应调速器死区、模型和参数都不同的情况下的多机电力系统发生超低频振荡的风险。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明具体实施例的统一频率模型示意图;
图2为本发明具体实施例2中电力系统的结构模型;
图3a为本发明具体实施例2在场景1下的示例电力系统频率响应曲线;
图3b为本发明具体实施例2在场景2下的示例电力系统频率响应曲线;
图3c为本发明具体实施例2在场景3下的示例电力系统频率响应曲线。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1:
计及多调速器死区的电力系统超低频振荡风险评估方法,
包括如下步骤:
S1.采集联网发电机及其调速器模型和运行参数;
S2.根据发电机的转动惯量、阻尼系数以及负荷频率调节效应确定等值发电机的传递函数Ggen,根据调速器模型确定各调速器线性部分的传递函数Ggovi,利用描述函数法计算各调速器死区环节的描述函数NLi;
S3.结合各调速器Ggovi及其死区环节NLi,计算其扩展描述函数的幅值|Gei|和相角∠Gei,将多个扩展描述函数相加组合成系统的扩展描述函数Ge;
S4.利用表示调速器及死区的扩展描述函数Ge和等值发电机模型Ggen,建立统一频率模型,如图1所示,图1中,fref为电力系统的额定运行频率,Δf为电力系统的运行频率偏差。
再根据极限环产生的条件,计算产生超低频频振荡的临界振幅Ac,临界振幅Ac越高,表明能够使系统产生振荡的最小扰动幅值越大,系统的稳定性越高,所以临界振幅Ac的大小能够衡量系统稳定性的高低。
上式中,Si表示第i台发电机的容量,TJi表示第i台发电机的转动惯量,m表示发电机的总数量。
在步骤S2中,调速器死区环节的描述函数NLi通过如下方法确定:
其中,b是死区的宽度,A0是振荡的幅值大小,在此作为待求参数。
在步骤S3中,扩展描述函数的幅值|Gei|和相角∠Gei通过如下方法确定:
|Gei(A0,ω)|=|NLi||Ggovi|
∠Gei(A0,ω)=∠NLi+∠Ggovi
幅值|Gei|和相角∠Gei是振荡幅值A0和振荡频率ω的函数。
然后,在步骤S4中,利用表示调速器及死区的扩展描述函数Ge和等值发电机模型Ggen,建立统一频率模型,如图1所示。
联立上述方程组,以振荡幅值A0和振荡频率ω为未知数进行求解,求解出的A0值即为临界幅值Ac,即使系统产生超低频振荡的最小扰动幅值。
实施例2:
本实施例所提供的电力系统结构模型如图2所示,在图2中,忽略了网络结构,发电机组由最常用的水电机组(水轮机+调速器)和火电机组(汽轮机+调速器)两个机组构成,发电机为等值发电机。其中,水轮机Ght及其调速器模型Ghgov分别为:
式中,TW为水锤效应时间常数,KD、KP、KI分别为调速器的微分、比例、积分系数,BP为调差系数,TG为伺服系统时间常数。
汽轮机及其调速器模型分别为:
式中,FHP为高压缸产生的功率在总轮机功率中占比,TRH为再热时间常数,TCH为主进气容积和气室时间常数。
将本申请提出的计及多调速器死区环节的电力系统超低频振荡风险评估方法运用至本实施例中:
本实施例在步骤S2中,水轮机死区描述函数NLh通过如下方法确定:
其中,b1是水轮机调速器死区的宽度,A0是振荡的幅值大小。
汽轮机死区描述函数NLs通过如下方法确定:
其中,b2是汽轮机调速器死区的宽度,A0是振荡的幅值大小。
本实施例在步骤S3中,含死区的水轮机及其调速器的扩展描述函数的幅值|Ghe|和相角∠Ghe分别为:
|Ghe(A0,ω)|=|NLh||Ghgov·Ght|
∠Ghe(A0,ω)=∠NLh+∠Ghgov·Ght
含死区的汽轮机及其调速器的扩展描述函数的幅值|Gse|和相角∠Gse分别为:
本实施例在步骤S3中的扩展描述函数Ge通过如下方法确定:
Ge=Ghe+Gse,Ge是输入扰动幅值A0和扰动频率ω的函数。
具体的,本实施例,即图2系统中,各模型参数为如下所示:
发电机Ggen参数:TJ=10.0s,DS=0.4;
水轮机Ght及其调速器Ghgov参数:KP=0.5,KD=0.7,KI=1,TW=1,Tgh=0.2,bp=0.04;
汽轮机Gst及其调速器Gsgov参数:R=0.0303,Tgg=0.2,Fhp=1,Trh=10,Tch=12。
另外,b1和b2分别是水轮机调速器和汽轮机调速器的死区。通过配置不同死区b1和b2的大小,利用本文方法计算出使系统产生振荡的临界振幅,结果如表1所示。
表1
对于每个不同死区配置的算例,都有一个临界振幅值,当扰动的幅值大于此临界值时,产生频率振荡,系统不稳定;当扰动的幅值在临界值内时,振荡无法持续,系统稳定。
从表1中可以看出,场景1、2、3的临界振幅值均大于对应的死区,所以分别以大于和小于临界振幅的两个扰动输入系统,得到如图3a~图3c所示的系统频率响应曲线。图3a~图3c分别为在不同死区配置场景下的示例电力系统频率响应曲线,从图中可以看出,在场景1中,临界振幅为0.0116pu,幅值Ad1=0.0113pu的扰动小于临界振幅,激发的振荡逐渐衰减;幅值Ad2=0.0122pu的扰动大于临界振幅,则产生了频率振荡,系统不稳定。场景2、3中也是类似的情况。因此通过上述计算结果表明,本发明计算出的临界振幅大小能够反应系统发生超低频振荡的风险高低。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.计及多调速器死区的电力系统超低频振荡风险评估方法,其特征在于,包括:
S1.采集联网发电机及其调速器模型和运行参数;
S2.确定等值发电机的传递函数Ggen,根据调速器模型确定各调速器线性部分的传递函数Ggovi,计算各调速器死区环节的描述函数NLi;
S3.结合各调速器Ggovi及其死区环节NLi,计算各调速器扩展描述函数的幅值|Gei|和相角∠Gei,将多个扩展描述函数相加组合成系统的扩展描述函数Ge;
S4.利用Ge和Ggen,建立统一频率模型;在统一频率模型的基础上,计算产生超低频振荡的临界振幅Ac,通过Ac评估电力系统超低频振荡的风险;
步骤S3中,扩展描述函数的幅值|Gei|和相角∠Gei通过下式确定:
|Gei(A0,ω)|=|NLi||Ggovi|
∠Gei(A0,ω)=∠NLi+∠Ggovi,
其中:A0为待求参数,表示振荡的幅值大小;ω为振荡频率;
步骤S4中,产生超低频振荡的临界振幅Ac,根据极限环产生的条件进行计算;
所述极限环产生的条件为:
|Ge|·|Ggen|=1,∠Ge+∠Ggen=-π/2。
2.根据权利要求1所述的计及多调速器死区的电力系统超低频振荡风险评估方法,其特征在于,步骤S4中,通过Ac评估电力系统超低频振荡的风险的方法为:临界振幅Ac越高、能够使系统产生振荡的最小扰动幅值越大、系统的稳定性越高。
3.根据权利要求1所述的计及多调速器死区的电力系统超低频振荡风险评估方法,其特征在于,步骤S2中,等值发电机的传递函数Ggen根据发电机的转动惯量、阻尼系数和负荷频率调节效应进行确定。
7.根据权利要求1所述的计及多调速器死区的电力系统超低频振荡风险评估方法,其特征在于,步骤S2中,各调速器死区环节的描述函数NLi通过描述函数法计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911201185.1A CN110829465B (zh) | 2019-11-29 | 2019-11-29 | 计及多调速器死区的电力系统超低频振荡风险评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911201185.1A CN110829465B (zh) | 2019-11-29 | 2019-11-29 | 计及多调速器死区的电力系统超低频振荡风险评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110829465A CN110829465A (zh) | 2020-02-21 |
CN110829465B true CN110829465B (zh) | 2023-02-24 |
Family
ID=69541912
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911201185.1A Active CN110829465B (zh) | 2019-11-29 | 2019-11-29 | 计及多调速器死区的电力系统超低频振荡风险评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110829465B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114398761B (zh) * | 2021-12-20 | 2023-04-07 | 华北电力大学 | 一种考虑调速器强化型死区的同步机系统频率响应解析计算方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2733607A1 (fr) * | 1995-04-28 | 1996-10-31 | Centre Nat Rech Scient | Procede et disposittif d'ajustement d'un regulateur pid |
CN1146532A (zh) * | 1996-01-10 | 1997-04-02 | 贾元丁 | 自动式重力动力发电系统 |
CN104682764A (zh) * | 2013-11-28 | 2015-06-03 | 哈尔滨功成科技创业投资有限公司 | 一种应用于飞轮储能系统电网侧变流装置控制器 |
CN105138739A (zh) * | 2015-07-30 | 2015-12-09 | 南方电网科学研究院有限责任公司 | 一种计及死区效应的电力系统最低频率快速计算方法 |
CN106357112A (zh) * | 2016-10-18 | 2017-01-25 | 东南大学 | 一种减小dc/dc变换器adc量化效应造成非线性的补偿方法 |
CN106849078A (zh) * | 2017-04-25 | 2017-06-13 | 国网四川省电力公司电力科学研究院 | 一种计及死区非线性的电力系统稳态频率预测方法 |
CN107070338A (zh) * | 2017-04-13 | 2017-08-18 | 李颖玉 | 电动车辆驱动系统 |
CN110224429A (zh) * | 2019-06-04 | 2019-09-10 | 中国南方电网有限责任公司 | 基于增强型调速器对电网超低频振荡稳定性分析方法 |
JP2019168777A (ja) * | 2018-03-22 | 2019-10-03 | 三菱重工業株式会社 | 制御系設計支援装置、制御系設計支援方法、及び制御系設計支援プログラム |
CN110492483A (zh) * | 2019-09-17 | 2019-11-22 | 华北电力大学 | 一种一次调频反馈通道非线性环节参数配置方法和系统 |
-
2019
- 2019-11-29 CN CN201911201185.1A patent/CN110829465B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2733607A1 (fr) * | 1995-04-28 | 1996-10-31 | Centre Nat Rech Scient | Procede et disposittif d'ajustement d'un regulateur pid |
CN1146532A (zh) * | 1996-01-10 | 1997-04-02 | 贾元丁 | 自动式重力动力发电系统 |
CN104682764A (zh) * | 2013-11-28 | 2015-06-03 | 哈尔滨功成科技创业投资有限公司 | 一种应用于飞轮储能系统电网侧变流装置控制器 |
CN105138739A (zh) * | 2015-07-30 | 2015-12-09 | 南方电网科学研究院有限责任公司 | 一种计及死区效应的电力系统最低频率快速计算方法 |
CN106357112A (zh) * | 2016-10-18 | 2017-01-25 | 东南大学 | 一种减小dc/dc变换器adc量化效应造成非线性的补偿方法 |
CN107070338A (zh) * | 2017-04-13 | 2017-08-18 | 李颖玉 | 电动车辆驱动系统 |
CN106849078A (zh) * | 2017-04-25 | 2017-06-13 | 国网四川省电力公司电力科学研究院 | 一种计及死区非线性的电力系统稳态频率预测方法 |
JP2019168777A (ja) * | 2018-03-22 | 2019-10-03 | 三菱重工業株式会社 | 制御系設計支援装置、制御系設計支援方法、及び制御系設計支援プログラム |
CN110224429A (zh) * | 2019-06-04 | 2019-09-10 | 中国南方电网有限责任公司 | 基于增强型调速器对电网超低频振荡稳定性分析方法 |
CN110492483A (zh) * | 2019-09-17 | 2019-11-22 | 华北电力大学 | 一种一次调频反馈通道非线性环节参数配置方法和系统 |
Non-Patent Citations (4)
Title |
---|
Analysis of power system frequency oscillations with intentional governor deadbands using describing functions;Lei Chen,等;《Electrical Power and Energy Systems》;20191031;第390–397页 * |
Estimation of Inherent Governor Dead-Band and Regulation Using Unscented Kalman Filter;Pratyasa Bhui,等;《IEEE TRANSACTIONS ON POWER SYSTEMS》;20180731;第33卷(第4期);第3546-3558页 * |
单机等值系统调速器死区对频率振荡的影响;陈磊,等;《电力系统自动化》;20190410;第43卷(第7期);第107-112页 * |
电力系统一次调频过程的超低频振荡分析;路晓敏,等;《电力系统自动化》;20170825;第41卷(第16期);第64-70页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110829465A (zh) | 2020-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Okou et al. | Power systems stability enhancement using a wide-area signals based hierarchical controller | |
Zenelis et al. | Online PMU-based wide-area damping control for multiple inter-area modes | |
Elkington et al. | On the stability of power systems containing doubly fed induction generator-based generation | |
JP2014168351A (ja) | 系統連系する電力変換装置 | |
CN109256803A (zh) | 虚拟同步机孤岛运行小信号建模及参数灵敏度计算方法 | |
CN112332459B (zh) | 一种基于灵敏度分析的多机系统调差系数优化方法 | |
CN110875600B (zh) | 一种两机等值电力系统动态频率响应近似解析模型 | |
Silva et al. | Assessment of power swings in hydropower plants through high-order modeling and eigenanalysis | |
Luo et al. | An optimal modal coordination strategy based on modal superposition theory to mitigate low frequency oscillation in FCWG penetrated power systems | |
CN110829465B (zh) | 计及多调速器死区的电力系统超低频振荡风险评估方法 | |
CN106684866A (zh) | 一种静态开环综合增益的计算方法及装置 | |
CN110649596B (zh) | 一种考虑系统初始状态的频率全响应解析模型 | |
CN114759618A (zh) | 基于动态能量确定双馈风电并网系统稳定性的方法及系统 | |
Wang et al. | Linear parameter varying control of a doubly fed induction generator based wind turbine with primary grid frequency support | |
Cheng et al. | Damping analysis and parameter adjustment of hydraulic turbine governor under ultra-low frequency oscillation | |
Bento et al. | Design of Power Systems Stabilizers for distributed synchronous generators using Linear Matrix Inequality solvers | |
CN116436042B (zh) | 计及风电机组调频死区的风水火系统稳定性分析方法 | |
CN116937698A (zh) | 一种含异构电源的电力系统小信号建模方法 | |
CN114243762A (zh) | 一种风机并网的分析与控制方法 | |
Calgan et al. | RSM-based optimization of excitation capacitance and speed for a self-excited induction generator | |
Zarei et al. | A robust direct current control of DFIG wind turbine with low current THD based predictive approach | |
Saxena et al. | Investigation of STATCOM performance with different participation of dynamic and static load in isolated hybrid power system | |
Silva et al. | A fully modular tool for small-signal stability analysis of hydroelectric systems | |
Wang et al. | Modelling analysis in power system small signal stability with grid-connected wind farms of DFIG Type | |
Kumar et al. | A novel pid-fopd controller for lfc including ipfc and rfb |
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 |