CN103577691A - 一种电网仿真中孤立电网频率计算方法 - Google Patents
一种电网仿真中孤立电网频率计算方法 Download PDFInfo
- Publication number
- CN103577691A CN103577691A CN201310529190.1A CN201310529190A CN103577691A CN 103577691 A CN103577691 A CN 103577691A CN 201310529190 A CN201310529190 A CN 201310529190A CN 103577691 A CN103577691 A CN 103577691A
- Authority
- CN
- China
- Prior art keywords
- frequency
- sys
- island
- simulation
- electric island
- 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
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种电网仿真中孤立电网频率计算方法,属于电力系统控制和仿真领域,该方法包括:将一个仿真周期划分为若干频率积分时段,在每个时段内基于发电机和负荷的频率特性实现仿真频率计算。在每一个频率积分时段计算出电气岛频率变化、电气岛内每个机组出力变化。机组出力变化改变下一个积分周期的电气岛的加速功率。对于每个频率积分时段,利用频率变化和频率积分时段计算结果,驱动频率自动装置动作切除部分机组或负荷设备。在每个频率积分时段,记录频率计算曲线。如果仿真电网频率无法稳定,则按照电气死岛处理。频率计算完成后,根据最终确定的电气岛发电、负荷情况,进行潮流计算,给出孤立电网各母线电压、相角和各支路潮流及损耗。
Description
技术领域
本发明属于电力系统调度自动控制与仿真技术领域,特别涉及一种孤立电网仿真频率计算方法。
背景技术
电网仿真需要模拟各种可能故障下电网情况,某些故障可能导致电网解列,形成孤立电网。传统的电网仿真对于孤立电网,当发电负荷差额大于一定比例时,通常按照电气死岛处理,不作详细计算。随着电网结构的逐步增强和各种分布式新能源的接入,部分解列后的孤立电网也具备短时期内独立运行的能力;同时电网中通常配置一些频率自动装置,在频率低于某个定值超过一段时间后,自动切除部分负荷或机组,保证孤立电网的频率稳定。传统的孤立电网频率计算方法主要通过孤立电网中的发电和负荷总量来计算频率,难以模拟频率自动装置的作用,无法准确模拟孤立电网频率是否能够保持稳定,导致孤立电网潮流计算结果不正确。
发明内容
本发明的目的是为了解决孤立电网频率难以准确计算的问题,提出一种电网仿真中孤立电网频率计算方法。本方法充分考虑了孤立电网中机组、负荷频率特性参数对频率计算的影响,仿真结果能够准确模拟频率自动装置的快速动作行为,为孤立电网仿真提供精确的潮流计算结果。
本发明提出的一种电网仿真中孤立电网频率仿真计算方法,其特征在于,该方法包括以下步骤:
1)将一个仿真周期划分为若干频率积分时段;具体步骤如下:
①定义频率积分时段Tp
Tp=Tc/Nstep (1)
式(1)中,Tc表示一个仿真周期;Nstep表示最大频率积分次数;
2)在一个仿真周期内对每个频率积分时段内基于发电机和负荷的频率特性进行仿真频率计算,包括计算每个孤立电气岛加速功率Pacc,频率变化Δfsys及发电机出力变化ΔPjG;
3)对于每个频率积分时段,利用频率变化Δfsys和频率积分时段Tp计算结果,仿真频率自动装置动作,根据频率判据切除部分机组或负荷设备;
4)在每个频率积分时段,记录计算的仿真频率曲线,如果仿真频率fsys<0.8*fbase,表示电气岛频率失稳,该频率积分时段按照电气死岛处理,频率计算结束;
当频率积分次数<最大频率积分次数Nstep时,循环进行2)-4),实现频率计算;当频率积分次数>最大频率积分次数Nstep时,该仿真计算周期频率计算完成,进入步骤5);
5)频率计算完成后,根据频率积分最终确定的电气岛机组出力、负荷情况,进行潮流计算,给出孤立电网各母线电压、相角和各支路潮流及损耗;等待下一个仿真计算周期到来转步骤2)。本发明的方法的技术特点:
本发明一方面,将仿真周期划分为若干频率积分时段,在每个时段内基于发电机和负荷的频率特性实现仿真频率计算。另一方面,在每个频率积分时段内模拟频率自动装置动作行为,仿真频率自动装置动作切除部分机组或负荷设备。最后记录仿真周期内整个频率计算曲线,给出孤立电网各机组和负荷的稳定结果。如果仿真电网频率无法稳定,则按照电气死岛处理。
本发明的优点如下:
充分考虑了孤立电网中机组、负荷频率特性参数对频率计算的影响,仿真结果逼真。
详细模拟频率自动装置在频率变化时的快速动作。通过频率积分,将频率变化的分辨率提升为0.1秒,能够准确模拟频率自动装置的快速动作行为。
实现了孤立电网的频率仿真计算,实时记录仿真过程中频率变化曲线,准确模拟低频情况下孤立电网的发电和负荷变化,为孤立电网仿真提供精确的潮流计算结果。
附图说明
图1为本发明方法的总体流程框图。
具体实施方式
本发明提出的电网仿真中孤立电网频率计算方法,结合附图及实施例详细说明如下:
本发明提出的电网仿真中孤立电网频率计算方法,如图1所示,包括以下步骤:
1)将一个仿真周期划分为若干频率积分时段;具体步骤如下:
①定义频率积分时段Tp
Tp=Tc/Nstep (1)
式(1)中,Tc表示一个仿真周期,通常是每5秒计算一次为一个仿真计算周期;Nstep表示最大频率积分次数,通常是50次;
2)在一个仿真周期内对每个频率积分时段内基于发电机和负荷的频率特性进行仿真频率计算,包括计算每个孤立电气岛加速功率,频率变化及发电机出力变化,仿真频率计算采用频率积分模式算法;具体步骤如下:
①对于解列后的多电气岛,计算每个电气岛的惯性时间常数TSum:
式(2)中,NG表示电气岛中的所有机组;TjUn表示第j个机组的惯性时间常数,
可以从发电机的额定参数中获取;
②对于解列后的多电气岛,计算每个电气岛的负荷频率特性系数DSum、电气岛功频特性系数Ks、电气岛时间常数Ts
Ks=1/DSum (4)
Ts=TSum/Dsum (5)
式(3)中,ND表示电气岛中的所有负荷对象;Djtyp表示第j个负荷的频率特性系数,可以从负荷参数中获取;WjD表示该负荷的实际值,由潮流计算结果给出;Δfsys=fsys-fbase表示该电气岛频率fsys与额定频率fbase(通常为50Hz)的偏差,fsys由潮流计算结果给出,其初值为50Hz;
③对于解列后的多电气岛,计算每个电气岛的加速功率Pacc
式(6)中,NL表示电气岛中的所有线路和绕组对象;WjG表示第j个机组的实际出力,由潮流计算结果给出;WjD表示第j个负荷的实际值,由潮流计算结果给出;Wjloss表示该电气岛第j个线路或绕组的损耗,由潮流计算结果给出;
④通过电气岛加速功率Pacc计算电气岛频率变化Δfsys、电气岛内每个机组出力变化ΔPjG:
WjG=WjG+ΔPjG (9)
式(8)中,KjUn表示表示第j个机组的功频特性系数,可以从发电机的额定参数中获取;
通过上述步骤,在每一个频率积分时段计算出电气岛频率变化Δfsys、电气岛内第j个机组出力变化ΔPjG、第j个机组出力为WjG=WjG+ΔPjG,并得到仿真频率fsys=fbase+Δfsys;
用当前电气岛的机组出力代入步骤③得到下一个频率积分时段电气岛的加速功率Pacc;进而计算出下一个频率积分时段电气岛的仿真频率fsys;以此类推得到每个频率积分时段电气岛的仿真频率fsys;
3)对于每个频率积分时段,利用频率变化Δfsys和频率积分时段Tp计算结果,仿真频率自动装置动作,根据频率判据切除部分机组或负荷设备;具体步骤如下:
⑤频率自动装置的判据为:仿真频率fsys低于低频定值fdz(从频率自动装置定值单获取)和判据延时tys大于时延定值tdz(从频率自动装置定值单获取);如果当前仿真频率fsys=fbase-Δfsys小于fdz,则频率判据成立,判据延时tys=tys+Tp;否则频率判据不成立,则判据延时tys=0;
⑥如果判据延时tys<时延定值tdz,则在该频率积分时段频率自动装置不动作;
⑦如果判据延时tys>时延定值tdz,则在该频率积分时段频率自动装置动作切除部分发电或者负荷,切除发电或负荷的数值从频率自动装置定值单获取,用以改变下一个积分频率时段的电气岛的加速功率Pacc;
4)在每个频率积分时段,记录计算的仿真频率曲线,如果仿真频率fsys<0.8*fbase,表示电气岛频率失稳,该频率积分时段按照电气死岛处理,频率计算结束;
当频率积分次数<最大频率积分次数Nstep时,2)-4)循环实现频率计算。当频率积分次数>最大频率积分次数Nstep时,该仿真计算周期频率计算完成,进入步骤5);
5)频率计算完成后,根据频率积分最终确定的电气岛机组出力、负荷情况,进行潮流计算,给出孤立电网各母线电压、相角和各支路潮流及损耗。等待下一个仿真计算周期到来转步骤2)。
Claims (3)
1.一种电网仿真中孤立电网频率仿真计算方法,其特征在于,该方法包括以下步骤:
1)将一个仿真周期划分为若干频率积分时段;具体步骤如下:
①定义频率积分时段Tp
Tp=Tc/Nstep (1)
式(1)中,Tc表示一个仿真周期;Nstep表示最大频率积分次数;
2)在一个仿真周期内对每个频率积分时段内基于发电机和负荷的频率特性进行仿真频率计算,包括计算每个孤立电气岛加速功率Pacc,频率变化Δfsys及发电机出力变化ΔPjG;
3)对于每个频率积分时段,利用频率变化Δfsys和频率积分时段Tp计算结果,仿真频率自动装置动作,根据频率判据切除部分机组或负荷设备;
4)在每个频率积分时段,记录计算的仿真频率曲线,如果仿真频率fsys<0.8*fbase,表示电气岛频率失稳,该频率积分时段按照电气死岛处理,频率计算结束;
当频率积分次数<最大频率积分次数Nstep时,循环进行步骤2)-4),实现频率计算;当频率积分次数>最大频率积分次数Nstep时,该仿真计算周期频率计算完成,进入步骤5);
5)频率计算完成后,根据频率积分最终确定的电气岛机组出力、负荷情况,进行潮流计算,给出孤立电网各母线电压、相角和各支路潮流及损耗;等待下一个仿真计算周期到来转步骤2)。
2.如权利要求1所述的方法,其特征在于,其中所述的步骤2)仿真频率计算采用频率积分模式算法;具体步骤如下:
②对于解列后的多电气岛,计算每个电气岛的惯性时间常数TSum:
式(2)中,NG表示电气岛中的所有机组;TjUn表示第j个机组的惯性时间常数,可以从发电机的额定参数中获取;
③对于解列后的多电气岛,计算每个电气岛的负荷频率特性系数DSum、电气岛功频特性系数Ks、电气岛时间常数Ts
Ks=1/DSum (4)
Ts=TSum/Dsum (5)
式(3)中,ND表示电气岛中的所有负荷对象;Djtyp表示第j个负荷的频率特性系数,可以从负荷参数中获取;WjD表示该负荷的实际值,由潮流计算结果给出;Δfsys=fsys-fbase表示该电气岛频率fsys与额定频率fbase的偏差,fsys由潮流计算结果给出;
④对于解列后的多电气岛,计算每个电气岛的加速功率Pacc
式(6)中,NL表示电气岛中的所有线路和绕组对象;WjG表示第j个机组的实际出力,由潮流计算结果给出;WjD表示第j个负荷的实际值,由潮流计算结果给出;Wjloss表示该电气岛第j个线路或绕组的损耗,由潮流计算结果给出;
⑤通过电气岛加速功率Pacc计算电气岛频率变化Δfsys、电气岛内每个机组出力变化ΔPjG:
WjG=WjG+ΔPjG (9)
式(8)中,KjUn表示表示第j个机组的功频特性系数,可以从发电机的额定参数中获取;
并得到仿真频率fsys=fbase+Δfsys;
⑥用当前电气岛的机组出力代入步骤③得到下一个频率积分时段电气岛的加速功率Pacc;进而计算出下一个频率积分时段电气岛的仿真频率fsys;以此类推得到每个频率积分时段电气岛的仿真频率fsys。
3.如权利要求1所述的方法,其特征在于其中所述的步骤3)具体包括以下步骤:
具体步骤如下:
⑦频率自动装置的判据为:仿真频率fsys低于低频定值fdz和判据延时tys大于时延定值tdz;如果当前仿真频率fsys=fbase-Δfsys小于fdz,则频率判据成立,判据延时tys=tys+Tp;否则频率判据不成立,则判据延时tys=0;
⑧如果判据延时tys<时延定值tdz,则在该频率积分时段频率自动装置不动作;如果判据延时tys>时延定值tdz,则在该频率积分时段频率自动装置动作切除部分发电或者负荷,切除发电或负荷的数值从频率自动装置定值单获取,用以改变下一个积分频率时段的电气岛的加速功率Pacc。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310529190.1A CN103577691B (zh) | 2013-10-30 | 2013-10-30 | 一种电网仿真中孤立电网频率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310529190.1A CN103577691B (zh) | 2013-10-30 | 2013-10-30 | 一种电网仿真中孤立电网频率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103577691A true CN103577691A (zh) | 2014-02-12 |
CN103577691B CN103577691B (zh) | 2016-09-28 |
Family
ID=50049455
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310529190.1A Expired - Fee Related CN103577691B (zh) | 2013-10-30 | 2013-10-30 | 一种电网仿真中孤立电网频率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103577691B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2945278A1 (en) * | 2014-05-12 | 2015-11-18 | Alstom Technology Ltd | Static exciter system for generators |
CN106655283A (zh) * | 2016-12-16 | 2017-05-10 | 国网上海市电力公司 | 一种针对孤网类型电网的频率仿真方法 |
CN111257634A (zh) * | 2020-03-26 | 2020-06-09 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种电网频率实时高精度测量新方法 |
CN112462626A (zh) * | 2020-11-25 | 2021-03-09 | 合肥工业大学 | 一种并网逆变器运行参数的仿真优化方法 |
CN113255485A (zh) * | 2021-05-13 | 2021-08-13 | 云南电网有限责任公司 | 一种水力发电机组并网模式的识别方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005151781A (ja) * | 2003-11-19 | 2005-06-09 | Chugoku Electric Power Co Inc:The | 電力需要予測システムおよび電力系統の負荷周波数制御方法 |
CN102930078A (zh) * | 2012-09-28 | 2013-02-13 | 南方电网科学研究院有限责任公司 | 一种孤立微网分布式电源容量与布点优化方法 |
-
2013
- 2013-10-30 CN CN201310529190.1A patent/CN103577691B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005151781A (ja) * | 2003-11-19 | 2005-06-09 | Chugoku Electric Power Co Inc:The | 電力需要予測システムおよび電力系統の負荷周波数制御方法 |
CN102930078A (zh) * | 2012-09-28 | 2013-02-13 | 南方电网科学研究院有限责任公司 | 一种孤立微网分布式电源容量与布点优化方法 |
Non-Patent Citations (1)
Title |
---|
周专等: "负荷频率特性对孤网频率稳定性影响分析", 《四川电力技术》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2945278A1 (en) * | 2014-05-12 | 2015-11-18 | Alstom Technology Ltd | Static exciter system for generators |
CN105099306A (zh) * | 2014-05-12 | 2015-11-25 | 阿尔斯通技术有限公司 | 用于发电机的静态励磁机系统 |
US9735719B2 (en) | 2014-05-12 | 2017-08-15 | General Electric Technology Gmbh | Static exciter system for generators |
CN105099306B (zh) * | 2014-05-12 | 2019-06-21 | 通用电器技术有限公司 | 用于发电机的静态励磁机系统 |
CN106655283A (zh) * | 2016-12-16 | 2017-05-10 | 国网上海市电力公司 | 一种针对孤网类型电网的频率仿真方法 |
CN111257634A (zh) * | 2020-03-26 | 2020-06-09 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种电网频率实时高精度测量新方法 |
CN111257634B (zh) * | 2020-03-26 | 2022-01-07 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种电网频率实时高精度测量新方法 |
CN112462626A (zh) * | 2020-11-25 | 2021-03-09 | 合肥工业大学 | 一种并网逆变器运行参数的仿真优化方法 |
CN112462626B (zh) * | 2020-11-25 | 2022-04-05 | 合肥工业大学 | 一种并网逆变器运行参数的仿真优化方法 |
CN113255485A (zh) * | 2021-05-13 | 2021-08-13 | 云南电网有限责任公司 | 一种水力发电机组并网模式的识别方法及装置 |
CN113255485B (zh) * | 2021-05-13 | 2022-05-17 | 云南电网有限责任公司 | 一种水力发电机组并网模式的识别方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103577691B (zh) | 2016-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103577691A (zh) | 一种电网仿真中孤立电网频率计算方法 | |
Zhang et al. | Synchrophasor measurement-based wind plant inertia estimation | |
Zaker et al. | A novel measurement-based dynamic equivalent model of grid-connected microgrids | |
Zhang et al. | Angle instability detection in power systems with high-wind penetration using synchrophasor measurements | |
CN104201674A (zh) | 一种考虑负荷低电压释放特性的综合负荷模型的建模方法 | |
Athari et al. | Stochastic cascading failure model with uncertain generation using unscented transform | |
Thakallapelli et al. | Real-time frequency based reduced order modeling of large power grid | |
Zhang et al. | Measurement-based power system dynamic model reductions | |
Piesciorovsky et al. | Comparison of non-real-time and real-time simulators with relays in-the-loop for adaptive overcurrent protection | |
CN109728597A (zh) | 一种光伏逆变器低电压穿越特性拟合方法及系统 | |
Sobbouhi et al. | Online synchronous generator out‐of‐step prediction by electrical power curve fitting | |
Suvorov et al. | Comprehensive validation of transient stability calculations in electric power systems and hardware-software tool for its implementation | |
Chiodo et al. | Probabilistic approach to transient stability evaluation | |
Guardado et al. | A PMU model for wide-area protection in ATP/EMTP | |
CN112952838B (zh) | 一种基于眼动设备的潮流智能分析方法及装置 | |
Osipov et al. | Adaptive model reduction for parareal in time method for transient stability simulations | |
Alvarez-Fernandez et al. | Parameterization of aggregated distributed energy resources (der_a) model for transmission planning studies | |
Kou et al. | Dynamic model validation with governor deadband on the eastern interconnection | |
Chin et al. | Adaptive and optimal under frequency load shedding | |
CN105186440A (zh) | 基于机电暂态方法的继电保护定值整定方法 | |
Meliopoulos et al. | Integration & automation: From protection to advanced energy management systems | |
Wilson et al. | Using dynamic simulations to design the wide-area stability and voltage control system (WACS) | |
Liu et al. | Using terminal circuit for power system electromagnetic transient simulation | |
Wang et al. | Case study of an improved cascading outage analysis model using outage checkers | |
CN105262046A (zh) | 基于机电暂态方法的继电保护定值整定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160928 Termination date: 20181030 |
|
CF01 | Termination of patent right due to non-payment of annual fee |