电力系统转动惯量在线实用化估计核心方法

1.本发明涉及电力系统运行状态量化评估领域，尤其涉及一种高比例新能源接入和大规模交直流电网结构下的电力系统转动惯量计算评估方法。


背景技术：

2.我国经济正从高速发展转向高质量发展。在能源领域，随着能源生产和消费革命持续推进，生产侧清洁化和消费侧电气化成为当前我国能源体系重要的趋势和特点。2019年国家电网公司提出到2050年我国能源清洁化率达到50％和终端电气化率达到50％的双50％目标
1.，将使我国已形成的世界上规模最大的交直流互联电网动态特性更加复杂、安全稳定运行风险加大。客观上对新形势下电力系统状态，特别是更具时效的转动惯量感知方法、更加精准的转动惯量计算技术和更加可靠的分析评估手段提出更高要求。
3.一个准确快速的系统惯量在线评估方法是在调控运行中时刻监视系统具有惯量水平的核心内容。关于惯量计算方法，文献均以发生切机扰动为计算前提。文献[7]-[11]从模拟电网事故contingencies角度出发，用扰动法将事故后功率和频率差额带入公式计算转动惯量。然而扰动法存在2个问题：1)只能计算切机事件当时的转动惯量，无法计算正常运行方式的转动惯量；2)需人工设置、整理计算数据，过程繁琐、费时，不满足在线计算需求。
[0004]
电力系统的转动惯量包括发电机和电动机及其拖动的转动机械的惯量。电源侧常规同步发电机、新能源机组和负荷侧带负载电动机均可为系统提供惯量支撑。以调度快速掌握系统实时转动惯量为出发点，在调度实时掌握当前运行方式机组类型、投入与否、出力、负荷等数据的基础上，利用统计法统计源、荷侧在线机组惯量。以此为线索，从物理构成要素出发，首先提出了系统整体转动惯量为发电侧发电机与负荷侧电动机惯量之和的物理概念。在理论计算源、荷侧惯量的基础上，推导出一种用统计法累加源、荷侧惯量的系统整体转动惯量在线实用化估计方法。


技术实现要素：

[0005]
本发明针对上述问题，在理论计算源、荷侧惯量的基础上，推导出一种用统计累加源、荷侧惯量的系统整体转动惯量在线实用化估计方法，所提算法可在计划层面对转动惯量进行摸底研究，在运行层面为实时控制策略提供参考，可复制推广性强。
[0006]
为了达到上述目的，本发明提供了一种电力系统转动惯量在线实用化估计核心方法，其包括：
[0007]
s1：定义电力系统转动惯量为电源侧同步发电机、非同步发电机和负荷侧带负载电动机的转动惯量和；
[0008]
s2：获取发电机铭牌参数，计算同步发电机单机的转动惯量jm，其计算公式为：
[0009]
[0010]
式中，hm为国际通用转动惯量时间常数，ss为系统容量基准值，f0为系统频率基准值；
[0011]
s3：计算所有并网的同步发电机的转动惯量j
sm
，其计算公式为：
[0012][0013]
式中，h1、h2、h3...hn为发电侧各台同步发电机转动惯量有名值；
[0014]
s4：计算非同步发电机的转动惯量j
asm
，其计算公式为：
[0015]jasm
＝k
×
p
asmp
[0016]
式中，p
asmp
为非同步发电机出力，k为区域电网非同步发电机的“转动惯量/额定功率”经验拟合系数；
[0017]
s5：计算负荷侧带负载电动机的转动惯量j
pl
，其计算公式为：
[0018]jpl
＝(g 1)j
mm
[0019]
式中，g为电动机惯量比，j
mm
为电动机自身转动惯量；
[0020]
s6：计算电力系统转动惯量js，其计算公式为：
[0021][0022]
优选方式下，所述非同步发电机出力p
asmp
的计算公式为：
[0023]
p
asmp
＝p
sum
×
α
[0024]
式中，p
sum
为某一运行方式下的系统总有功出力，α为某一运行方式下系统新能源出力占比。
[0025]
优选方式下，负荷侧带负载电动机的转动惯量j
pl
包括电动机自身转动惯量j
mm
以及负载折算在电机转子轴上的转动惯量j
lm
；
[0026]
所述电动机惯量比g的计算公式为：
[0027][0028]
本发明的有益效果为：本发明所提算法计算得到的转动惯量可有效评估高比例新能源接入下的电网安全水平，即对电网频率调节能力的量化评估；该结果可指导电力系统电源结构比例，可为新能源接入比例的边界控制提供量化依据。
附图说明
[0029]
图1为本发明的电力系统转动惯量各组成部分示意图。
具体实施方式
[0030]
理论分析与公式推导：
[0031]
1、转动惯量的具体组成
[0032]
转动惯量物理意义是刚体绕轴旋转时惯性的量度，发电机转动惯量是描述发电机或等效发电机所产生的惯性。
[0033]
发电机单台机的惯性时间常数国际通用转动惯量时间常数国际通用转动惯量时间常数与发电机转动惯量jm的关系为
[0034][0035]
实际中，一般用h表示t，二者关系为
1.。
[0036][0037]
sm为发电机额定容量，单位mva；ss为系统基准值，单位mva，单位为s，jm单位为kg.m2。
[0038][0039]
sm为单机容量，单位mva，hm是发电机惯性时间常数的有名值，单位mw.s。由ω＝2πf得到：
[0040][0041]
由公式(5)得：
[0042][0043]
若假设系统ss为100mva，100mva换算成与jm单位kg.m2对应的数量级是105kw，假设f0为50hz，则由式(1)可知，jm约等于hm：
[0044][0045]
根据集总参数电路理论
[2]
，可将n个电机组成的系统等效成一个电机模型，此时系统总转动惯量由式(4)单机转动惯量转变为：
[0046][0047]
δpr为电网发生扰动的瞬间功率变化，其可以是缺额，也可是盈余，jr为扰动发生的瞬间系统总转动惯量，式(7)即为已知功率变化量和频率变化率后计算系统整体转动惯量的“扰动法”基本原理。
[0048]
在考虑发电侧的非同步发电机如新能源、储能设备等发电设备转动惯量前提下，电力系统整体转动惯量由发电侧同步发电机、非同步发电机和恒功率负荷组成。
[0049]jt
＝jg j
l
＝j
sm
j
asm
j
pl
ꢀꢀꢀ
(8)
[0050]
式中，j
t
为统计法计算出的系统整体转动惯量，jg为发电侧转动惯量，j
l
为负荷侧转动惯量，j
sm
为同步发电机转动惯量，j
asm
为非同步发电机转动惯量，j
pl
为负荷侧恒功率负
荷转动惯量。式(8)即为已知发电侧在线发电机参数和负荷侧在线电动机参数后统计计算出的系统整体转动惯量“统计法”基本原理。
[0051]
根据式(5)(6)(7)，扰动法计算得到的系统整体转动惯量应等于统计计算出的发电侧与负荷可提供的转动惯量之和，即：
[0052]
jr＝j
t
ꢀꢀꢀ
(9)
[0053]
2、源侧同步发电机惯量的计算方法
[0054]
在已知发电机铭牌参数后，首先计算发电机单机转动惯量，将所有并网发电机的转动惯量累加，获得发电侧的转动惯量。利用集总参数电路理论，得：
[0055][0056]
式中h1、h2、h3...hn为发电侧单台同步发电机转动惯量有名值。
[0057]
由于国内目前只对同步发电机有强制甩负荷试验规定，进而给出惯性时间常数，因此式(10)仅能统计累加给定方式下在线同步发电机转动惯量。
[0058]
3、非同步发电机惯量估计
[0059]
目前，有若干惯量存在于系统中无法直接计算。如新能源机组，其并网逆变器本身属于不具备转动惯量的静止元件。虚拟同步机可以使并网逆变器具有类似于同步发电机的外特性[引，从而使逆变器参与调频调节并为电力系统增加虚拟惯量。但虚拟同步机技术存在2个问题：1)注重惯量模拟的控制回路，忽略了惯量模拟的能量来源与时候恢复；2)能够参与调频的逆变器许配备足够的储能。由于定向研究的某区域电网实际投运的虚拟同步机技术和储能设备较少，因此无法直接利用虚拟同步机技术计算新能源机组客观存在的转动惯量。
[0060]
由文献
[4]
可知，非同步发电机的虚拟同步机模型中转子运动方程中转动惯量与虚拟同步机输出有功功率正相关，因此区域电网在线非同步发电机总转动惯量j
asm
与其额定功率p
asme
正相关。根据“最大功率点跟踪”即maximum power point tracking原理
[5][6]7.，为了追求经济效益，目前普遍存在风电和光伏为代表的新能源发电机组并网即满发的现象，即
[0061]
p
asme
＝p
asmp
ꢀꢀꢀ
(11)
[0062]
式中p
asme
为非同步发电机额定功率，p
asmp
为非同步发电机出力，假设系统中非同步发电机出力占比为α，则有
[0063]
p
asme
＝p
asmp
＝p
sum
×
α
ꢀꢀꢀ
(12)
[0064]
p
sum
为某一运行方式下，系统总有功出力。α为某一运行方式下系统新能源出力占比。可将区域电网的j
asm
与其额定功率p
asme
的正相关性简化成倍数关系，即设定一个系数k，由式(8)(9)得：
[0065]jasm
＝k
×
p
sum
×
α＝j
t-j
sm-j
pl
ꢀꢀꢀ
(13)
[0066]
k为区域电网非同步发电机的“转动惯量/额定功率”经验拟合系数。在给定运行方式下，k可以通过扰动法得到的j
t
与j
sm
和j
pl
之差后计算求得。即可通过大量计算j
sum
和同方式下的j
sm
、j
pl
和p
asm
，将k利用历史数据拟合的方式估计出来，得出某区域电网的k值。
[0067]
4、荷侧电动机惯量的计算方法
[0068]
负荷是系统中所有电力用户的用电设备所消耗的电功率总和，是所有用户的负荷总加。电力系统综合负荷模型具有区域性和时间性，在没有精确综合负荷模型的情况下，典型区域电网一般按50％恒功率 50％恒阻抗经验值计算。
[0069]
恒功率负荷模型可以简化成一个电动机带负载。电动机自身转动惯量为j
mm
，负载折算在电机转子轴上的转动惯量为j
lm
，电动机与负载整体转动惯量为j
pl
，3者关系为：
[0070]jpl
＝j
mm
j
lm
ꢀꢀꢀ
(14)
[0071]jpl
与j
mm
的比例为电动机惯量比
[8]
g，单位100％，计算公式为：
[0072][0073]
进而得恒功率负荷的转动惯量为：
[0074]jpl
＝(g 1)j
mm
ꢀꢀꢀ
(16)
[0075]
在电机系统选型时，一般要求惯量比g在3-20之间
[8][9][10]
。未来需要对负荷侧进行精确建模和实时分类分析，获取动态的、实时的恒功率恒阻抗模型。
[0076]
5、快速实用化估计方法
[0077]
若假设f0为50hz、s为100mva，根据式(8)(9)(13)(16)得：
[0078][0079]
给定运行方式下，扰动法求解jr时容易受到实际扰动发生次数和仿真软件限制。如图1所示，而快速实用化估计方法求解j
sm
、j
acm
和j
pl
不受外界条件限制，过程简单、快速。因此可采取快速计算系统整体转动惯量js的方式替代扰动法求解j
t
：
[0080]js
＝jr＝j
t
ꢀꢀꢀ
(18)
[0081]
式(17)为电力系统转动惯量在线实用化估计核心方法。
[0082]
具体实施例：
[0083]
1、仿真计算拟合求取k
[0084]
以2019-2020年某区域电网实际数据为基础进行仿真计算，计算转动惯量。
[0085]
设定1％～10％，每档递增1％的10个新能源占比级别，每个级别配比2019-2020年该区域电网的冬大、冬小、夏大、夏小等10个运行方式，共100个运行方式。用扰动法模拟切除某1000mw级最大转动惯量单台发电机，将电力系统分析综合程序psasp即power system analysis software package，仿真得到系统频率变化曲线数据导入matlab程序，计算jr。同时，分别将所有并网的同步发电机厂家给定的甩负荷试验得到的值计算出j
sm
。进而对未知参数进行拟合。
[0086]
每个新能源占比对应的10个运行方式的负荷、j
sm
、j
pl
、j
t
及k的平均值，如表1所示为100个运行方式的转动惯量及相关参数，其中每10个为一组求平均值。求得j
t
值后，根据式(16)，拟合计算出了k，表中100个方式k的平均经验值为2.303。
[0087]
编号占比(％)负荷平均值(mw)j
sm
平均值(mw.s)j
pl
平均值(mw.s)j
t
平均值(mw.s)k平均值1163723.35489787.6254893.4744963.022.012251149.79382966.4204599.16589614.482.883362003.25465895.5248013719074.612.53
4470787519973.4283148810081.752.325563523.35469269.9254093.4729570.391.946659299.89431544.7237199.56677306.642.297758355.33422511.4233421.32664840.262.098855313.85381974.7221255.4613802.682.119963520.24445857254080.96715737.662.47101058272385706.81233088642857.112.39
[0088]
2、国内电网仿真算例验证
[0089]
2019年11月某日，国内某区域电网发生了大容量发电机脱网事件。该机脱网时出力1126mw，造成电网频率由50.05hz下降至49.90hz。故障时该全网发电功率6221万千瓦，火电占比69％，新能源占比21％，核电占比9.2％，负荷约58322mw，其不算厂用负荷。故障发生后2秒内df/dt约0.03478hz/s。扰动法计算得转动惯量jr为809373mw.s。
[0090]
2020年7月某日，国内某区域电网再次发生大型单台发电机跳闸事件。该机脱网前出力1031mw。故障时总功率59890.70mw，同步发电机出力55548mw，新能源占比7.25％，负荷约56538mw。电网频率从49.96hz下降至49.88hz。故障发生后2秒内df/dt约0.0325hz/s。扰动法计算得转动惯量jr为793076mw.s。
[0091]
将2次事件中实测数据带入式(10)求得j
t
，将仿真计算求得的k平均值2.303代入式(17)求js，通过对比j
t
和js，对快速实用化算法进行验证，如表2所示为2019年和2020年系统发生故障时的系统转动惯量及相关参数。2次事故js与j
t
误差在-0.084％～0.056％之间，误差很小。
[0092][0093]
2019、2020年故障中，用扰动法计算jr值后，根据式(17)，计算出了k分别为2.36和2.20，与仿真中计算得到的k值(2.303)相差2.48％和-4.47％，均在
±
5％之内，证明仿真结果较为准确，式(17)(18)可对该区域转动惯量进行快速、准确的在线计算。
[0094]
本发明使用的技术术语：
[0095]
转动惯量：转动惯量是刚体绕轴旋转时惯性的量度，发电机转动惯量是描述发电机或等效发电机所产生的惯性。
[0096]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。