CN108830013B - 一种基于动能定理的系统扰动下惯性时间常数评估方法 - Google Patents

一种基于动能定理的系统扰动下惯性时间常数评估方法 Download PDF

Info

Publication number
CN108830013B
CN108830013B CN201810696060.XA CN201810696060A CN108830013B CN 108830013 B CN108830013 B CN 108830013B CN 201810696060 A CN201810696060 A CN 201810696060A CN 108830013 B CN108830013 B CN 108830013B
Authority
CN
China
Prior art keywords
kinetic energy
time constant
disturbance
inertia
power system
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
Application number
CN201810696060.XA
Other languages
English (en)
Other versions
CN108830013A (zh
Inventor
李东东
张佳乐
赵耀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai University of Electric Power
Original Assignee
Shanghai University of Electric Power
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai University of Electric Power filed Critical Shanghai University of Electric Power
Priority to CN201810696060.XA priority Critical patent/CN108830013B/zh
Publication of CN108830013A publication Critical patent/CN108830013A/zh
Application granted granted Critical
Publication of CN108830013B publication Critical patent/CN108830013B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Public Health (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种基于动能定理的系统扰动下惯性时间常数评估方法，包括：在电力系统的扰动下，一次调频还未响应时，假设系统中有功功率的缺额完全由发电机旋转机构所释放的动能来弥补，利用动能定理得到电力系统的惯性时间常数。与现有技术相比，本发明将动能定理与传统的惯量评估方式相结合，不需要对频率变化率进行滤波取值即可获得较为准确的系统惯性时间常数值，而且也可以减小由于扰动发生时间的不精确测量带来的误差，降低了对电网惯量评估的复杂程度，具有良好的实用性。

Description

一种基于动能定理的系统扰动下惯性时间常数评估方法
技术领域
本发明涉及电力系统运行分析领域，尤其是涉及一种基于动能定理的系统扰动下惯性时间常数评估方法。
背景技术
为了解决能源危机和气候问题，各国政府加大了对新能源发电的投资力度，风力发电、太阳能发电的渗透率逐年上升，与此同时，新能源发电对电网产生的不良影响也日益凸显。
相对于传统的发电方式，新能源发电有着很大的区别。首先，新能源发电存在着间歇性，这使得电网波动越来越大，对整个电网的稳定性产生了影响。其次，新能源发电一般会经过电力电子器件和大电网进行并网，这造成了发电机部分和电网的解耦，使得发电机部分的惯量无法传递到电网中。研究发现，惯量越低的系统受扰动影响越大，频率跌落越快，电网稳定性更差，因此计算出整个系统的惯性时间常数是很有必要的。
惯性时间常数H定义为发电机仅通过释放旋转机构中存储的动能提供发电机额定功率的时间，一般单位为秒。
目前计算系统惯性时间常数最为常见的方式是在已知的系统扰动下利用发电机摇摆方程计算系统惯性时间常数。公式表达如式(1)所示
Figure BDA0001713525050000011
式(1)中，Hsys为电力系统的惯性时间常数；Ssys为系统发电机总容量；f为系统当前频率值；fn为系统的额定频率值；Pm为系统总的发电机机械功率；Pe为系统里发电机发出的总电磁功率；ΔP为系统总的有功功率缺额值。
其中，
Figure BDA0001713525050000012
式(2)中，HSG为单台发电机的惯性时间常数值，SSG为单台发电机容量。
但是，利用这种方式计算惯性时间常数的复杂程度也是很大的。为了获得系统的惯性时间常数值，就必须要先获得已知扰动下某一个节点的频率曲线，根据频率曲线计算得到频率变化率
Figure BDA0001713525050000021
的曲线，并对其进行滤波取值代入式(1)进行计算。在此基础上，扰动发生时间的不精确测定会给频率变化率曲线带来影响，从而增加惯量评估的误差，无论是在扰动发生时间的精确测定上，还是在频率变化率的滤波取值上，都给计算系统惯性时间常数带来了难度。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于动能定理的系统扰动下惯性时间常数评估方法。
本发明的目的可以通过以下技术方案来实现：
一种基于动能定理的系统扰动下惯性时间常数评估方法，包括：
在电力系统的扰动下，一次调频还未响应时，假设系统中有功功率的缺额完全由发电机旋转机构所释放的动能来弥补，利用动能定理得到电力系统的惯性时间常数。
优选的，所述利用动能定理得到的电力系统的惯性时间常数为：
Figure BDA0001713525050000022
其中，Hsys为电力系统的惯性时间常数，ΔP为电力系统发生扰动后一次调频未切入时的有功功率缺额值，Ssys为电力系统的发电机容量，fn为电力系统的额定频率值，t0为电力系统扰动发生时或之后的某一时刻，t1为电力系统一次调频切入时刻或选取的一次调频响应前的某一时刻，t0<t1，f1、f0分别为t1、t0时刻对应的系统频率值。
优选的，所述t0为电力系统扰动发生时刻。
优选的，所述t1为电力系统一次调频切入时刻。
优选的，所述t1为电力系统扰动发生后1秒时刻。
优选的，所述电力系统的扰动包括发电机故障切除。
与现有技术相比，本发明将动能定理与传统的惯量评估方式相结合，不需要对频率变化率进行滤波取值即可获得较为准确的系统惯性时间常数值，而且也可以减小由于扰动发生时间的不精确测量带来的误差，降低了对电网惯量评估的复杂程度，具有良好的实用性。
附图说明
图1为实施例中新英格兰10机39节点系统结构示意图；
图2为实施例中切除发电机G03时母线BUS01上频率仿真结果图；
图3为实施例中切除发电机G06时母线BUS01上频率仿真结果图；
图4为实施例中切除发电机G09时母线BUS01上频率仿真结果图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
实施例
一种基于动能定理的系统扰动下惯性时间常数评估方法，包括：
在电力系统不同程度的扰动基础上，忽略负荷的频率特性，在一次调频还未响应时，假设系统中有功功率的缺额完全由发电机旋转机构所释放的动能来弥补，利用动能定理得到电力系统的惯性时间常数。
本方法假定扰动发生后一次调频未响应时刻，忽略负荷频率特性的基础上，发电机通过释放旋转机构存储的动能完全弥补了系统的有功功率缺额，利用动能定理得到表达式如下：
Figure BDA0001713525050000031
上式中，Pinertia为发电机在t0到t1时段内释放旋转机构中存储的动能所提供的有功功率值；ΔP为系统发生扰动后一次调频未切入时的有功功率缺额值；Winertia为发电机在t0到t1时段内做的功；ΔE为t0到t1时段内释放的动能值；JSG为整个系统中同步电机的转动惯量值；w1为t1时刻的系统角频率值，w0为t0时刻的系统角频率值；t1为电力系统一次调频切入时刻或选取的一次调频响应前的某一时刻，t0为扰动发生时或之后的某一时刻，应保证t0<t1且t0到t1的时间范围在扰动发生到一次调频切入的时间范围内。
惯性时间常数的定义为：
Figure BDA0001713525050000041
式中，Hsys为电力系统的惯性时间常数。
将式(4)代入式(3)得到式(5)：
Figure BDA0001713525050000042
式中，Ssys为电力系统的发电机容量。
整理上式，得到的电力系统的惯性时间常数为：
Figure BDA0001713525050000043
其中，fn为电力系统的额定频率值，f1、f0分别为t1、t0时刻对应的系统频率值。
从式(6)可以看出惯性时间常数直接与频率联系，不需要考虑频率变化率的滤波取值，同时也减小了扰动发生时刻的不精确测量对惯量评估带来的误差影响，提高了计算惯性时间常数的速度，降低了惯量评估的复杂程度，具有一定的实用性。
本实施例中，采用新英格兰10机39节点系统(其负荷均设置为恒功率负荷模型)在DIgSILENT PowerFactory仿真软件中进行验证，该系统的拓扑结构如图1所示。验证方法为：设置不同程度的发电机扰动事件，按照公式(2)计算出系统惯性时间常数的实际值，然后在母线上获取合适的t0和t1值以及相对应的频率值，结合扰动后有功功率缺额和系统的发电机容量、额定频率，按照本方法中公式(6)得到惯性时间常数的计算值，将惯性时间常数的计算值和实际值进行比对，计算误差绝对值，若误差绝对值很小，则说明本方法计算得到的惯性时间常数值和系统实际惯性时间常数值基本一致，从而验证了本方法的正确性。
为了验证本方法的正确性，本实施例中设置了3种不同的发电机扰动事件：
扰动事件1：在仿真t＝2s时刻对发电机G03进行切机；
扰动事件2：在仿真t＝2s时刻对发电机G06进行切机；
扰动事件3：在仿真t＝2s时刻对发电机G09进行切机。
参数设置：取t0为扰动发生时刻，即t0＝2s时刻，t1为扰动发生后1秒时刻，即t1＝3s，频率值均在母线BUS01上取值。
上述3个扰动事件相应的BUS01上频率曲线仿真图分别如图2～4所示，由图中可以得到不同的发电机扰动事件对应的t＝3s时刻频率值分别为59.738HZ、59.746HZ、59.666HZ，结合不同扰动事件的有功功率缺额值、系统的发电机容量值、系统额定频率值可以得到惯性时间常数的计算值。经过对比发现，本方法得到的扰动下惯性时间常数计算值相比于系统的惯性时间常数实际值误差较小，说明惯性时间常数计算值与实际值基本一致。
通过上述验证过程可以得出，本方法能较为准确的计算出系统的惯性时间常数并且相对于传统的惯量评估方法，本方法复杂程度较低，计算速度较快，具有良好的实用价值。

Claims (2)

1.一种基于动能定理的系统扰动下惯性时间常数评估方法，其特征在于，包括：在电力系统的扰动下，一次调频还未响应时，假设系统中有功功率的缺额完全由发电机旋转机构所释放的动能来弥补，利用动能定理得到电力系统的惯性时间常数；
所述利用动能定理得到的电力系统的惯性时间常数为：
Figure FDA0002974539350000011
其中，Hsys为电力系统的惯性时间常数，ΔP为电力系统发生扰动后一次调频未切入时的有功功率缺额值，Ssys为电力系统的发电机容量，fn为电力系统的额定频率值，t0为电力系统扰动发生时刻，t1为电力系统扰动发生后1秒时刻，t0<t1，f1、f0分别为t1、t0时刻对应的系统频率值。
2.根据权利要求1所述的一种基于动能定理的系统扰动下惯性时间常数评估方法，其特征在于，所述电力系统的扰动包括发电机故障切除。
CN201810696060.XA 2018-06-29 2018-06-29 一种基于动能定理的系统扰动下惯性时间常数评估方法 Active CN108830013B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810696060.XA CN108830013B (zh) 2018-06-29 2018-06-29 一种基于动能定理的系统扰动下惯性时间常数评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810696060.XA CN108830013B (zh) 2018-06-29 2018-06-29 一种基于动能定理的系统扰动下惯性时间常数评估方法

Publications (2)

Publication Number Publication Date
CN108830013A CN108830013A (zh) 2018-11-16
CN108830013B true CN108830013B (zh) 2021-05-04

Family

ID=64134654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810696060.XA Active CN108830013B (zh) 2018-06-29 2018-06-29 一种基于动能定理的系统扰动下惯性时间常数评估方法

Country Status (1)

Country Link
CN (1) CN108830013B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109921472B (zh) * 2019-03-11 2022-11-04 上海电力学院 基于粒子群优化算法的电力系统等效惯量评估方法
CN110474323B (zh) * 2019-07-29 2020-12-22 上海电力大学 一种电力系统惯性时间常数测量方法
CN111245014B (zh) * 2020-02-27 2022-03-29 上海电力大学 一种基于虚拟惯量的电力系统控制方法
CN111293686A (zh) * 2020-02-29 2020-06-16 上海电力大学 基于armax系统辨识的电力系统惯量实时评估方法
CN112421650B (zh) * 2020-09-29 2022-07-01 中国电力科学研究院有限公司 一种评估电网惯性的方法和系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102064548A (zh) * 2010-12-17 2011-05-18 中国电力科学研究院 一种考虑暂态频率偏移的一次调频旋转备用优化方法
CN104538980A (zh) * 2015-01-28 2015-04-22 国家电网公司 一种微电网自平衡快速减负荷控制方法
KR101673527B1 (ko) * 2016-01-15 2016-11-07 연세대학교 산학협력단 인버터 기반 독립형 마이크로그리드에서의 풍력발전 관성 제어 시스템 및 방법
CN106910142A (zh) * 2017-02-17 2017-06-30 三峡大学 一种含风电有功‑频率耦合作用的电力系统频率特性计算方法
CN107579528A (zh) * 2017-07-18 2018-01-12 国家电网公司 基于能量变化的暂态频率安全紧急控制决策方法
CN107968415A (zh) * 2017-11-27 2018-04-27 上海电力学院 一种虚拟同步发电机的自适应虚拟惯量控制方法
CN108199418A (zh) * 2017-12-29 2018-06-22 湖南大学 一种互联电网暂态稳定紧急控制切机方案制定方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102064548A (zh) * 2010-12-17 2011-05-18 中国电力科学研究院 一种考虑暂态频率偏移的一次调频旋转备用优化方法
CN104538980A (zh) * 2015-01-28 2015-04-22 国家电网公司 一种微电网自平衡快速减负荷控制方法
KR101673527B1 (ko) * 2016-01-15 2016-11-07 연세대학교 산학협력단 인버터 기반 독립형 마이크로그리드에서의 풍력발전 관성 제어 시스템 및 방법
CN106910142A (zh) * 2017-02-17 2017-06-30 三峡大学 一种含风电有功‑频率耦合作用的电力系统频率特性计算方法
CN107579528A (zh) * 2017-07-18 2018-01-12 国家电网公司 基于能量变化的暂态频率安全紧急控制决策方法
CN107968415A (zh) * 2017-11-27 2018-04-27 上海电力学院 一种虚拟同步发电机的自适应虚拟惯量控制方法
CN108199418A (zh) * 2017-12-29 2018-06-22 湖南大学 一种互联电网暂态稳定紧急控制切机方案制定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Inertia time constant design in microgrids with multiple paralleled virtual Inertia time constant design in microgrids with multiple paralleled virtual;Fang Zhuo;《EPE"17 ECCE Europe》;20180402;第61-67页 *
基于频率影响因素的低频减载策略;岑炳成;《电气系统自动化》;20160610;第40卷(第11期);第1-9页 *
风电场等效虚拟惯性时间常数计算;李世春;《电气系统自动化》;20160410;第40卷(第7期);第22-29页 *

Also Published As

Publication number Publication date
CN108830013A (zh) 2018-11-16

Similar Documents

Publication Publication Date Title
CN108830013B (zh) 一种基于动能定理的系统扰动下惯性时间常数评估方法
CN103279639B (zh) 基于响应的受端电网电压稳定全过程态势评估及防控方法
CN104794347B (zh) 一种风电机组机舱传递函数分区拟合方法
CN103500269A (zh) 一种双馈异步发电机组暂态短路电流的计算方法
CN105649875B (zh) 风力发电机组的变桨控制方法及装置
CN105678033B (zh) 一种适用于电磁暂态仿真的风电场等值建模方法
US11340274B2 (en) Method and device for locating equipment-level oscillation sources of DFIG grid-connected system
WO2017016020A1 (zh) 电力系统暂态稳定性在线分析的分区合成方法及其装置
CN109995023A (zh) 基于谱聚类区域划分的电力系统稳定性的判断方法和系统
CN110518631A (zh) 一种直驱风电机组的稳定性评估方法及系统
CN106712097A (zh) 一种含大规模风电场的电力系统暂态稳定风险评估方法
CN110350588A (zh) 一种双馈风机并网系统动态能量稳定性评估方法及系统
CN112434936A (zh) 一种电力系统惯量安全域评估方法、系统、电子设备及可读存储介质
CN104504285A (zh) 一种计及撬棒保护的双馈风电场等值建模方法
CN105186502A (zh) 基于安全域的含双馈风机电力系统暂态稳定性分析方法
CN105243187B (zh) 一种基于暂态稳定测度指标的参数影响程度分析方法
CN109921472B (zh) 基于粒子群优化算法的电力系统等效惯量评估方法
Shabani et al. Real-time transient instability detection in the power system with high DFIG-wind turbine penetration via transient energy
CN107579531A (zh) 一种计及风速不确定性的双馈风电附加阻尼控制器的“域”设计方法
Zhang et al. Dynamic state estimation of a grid-connected converter of a renewable generation system using adaptive cubature Kalman filtering
CN104268316B (zh) 基于离线计算和在线匹配的双馈风电场概率等值建模方法
WO2024041590A1 (zh) 一种风电光伏联合发电的功率控制方法
CN105244887A (zh) 电力系统暂态功角失稳的闭环控制方法
CN103647284B (zh) 一种解决单时间断面问题的电压稳定预测方法
CN105701265A (zh) 一种双馈风电机组建模方法及装置

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