CN107947172A - 一种基于广域信息的电力系统惯性水平评估方法 - Google Patents
一种基于广域信息的电力系统惯性水平评估方法 Download PDFInfo
- Publication number
- CN107947172A CN107947172A CN201711364766.8A CN201711364766A CN107947172A CN 107947172 A CN107947172 A CN 107947172A CN 201711364766 A CN201711364766 A CN 201711364766A CN 107947172 A CN107947172 A CN 107947172A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- synchronous generator
- slow
- matrix
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000001360 synchronised effect Effects 0.000 claims abstract description 142
- 230000001427 coherent effect Effects 0.000 claims abstract description 67
- 230000008859 change Effects 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims description 99
- 239000013598 vector Substances 0.000 claims description 24
- 230000009897 systematic effect Effects 0.000 claims description 15
- 230000003416 augmentation Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 claims description 6
- 238000012512 characterization method Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims description 2
- 238000010248 power generation Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
-
- 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
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于广域信息的电力系统惯性水平评估方法,将电力系统同步发电机组分为多个慢同调机群,对每个慢同调机群,计算其包含的同步发电机组的参与因子和相关因子,选择起主导作用的机组;根据慢同调机群主导机组转子运动方程,确定主导机组电角速度、电角速度变化率、惯性时间常数与有功功率差额之间的关系;根据主导机组有功功率差额与惯性时间常数等参数间的关系,利用各慢同调机群主导机组的广域信息,得到电力系统总有功功率差额的估计值,用来评估电力系统整体惯性水平。本发明保证了系统整体惯性水平评估的准确性。根据慢同调理论对电力系统进行分群,减少了惯性水平评估中所需获取的信息量,显著提高了计算效率。
Description
技术领域
本发明涉及电力系统运行与控制技术领域,特别是涉及一种基于广域信息的电力系统惯性水平评估方法。
背景技术
随着以风电为代表的可再生能源发电的开发和利用,电力系统能源结构也随之变化,火力发电机组等传统同步发电机组所占比例逐渐降低。传统上,当电力系统遭受扰动时,同步发电机组通过调节存储在旋转质量块中的动能,补偿由扰动引起的系统总有功功率差额,以此抑制系统频率的波动。可再生能源发电机组通常通过电力电子换流器接入电力系统,实现解耦控制,不能直接响应系统的频率波动。
为了反映可再生能源发电并网对电力系统频率响应能力的影响,需要对系统的惯性水平进行准确的评估。在电力系统出现有功功率差额时,频率波动特性受系统惯性水平的影响。单台同步发电机组的惯性水平可以由其惯性时间常数来表征,但是系统整体惯性水平受多种因素影响,并不是各台同步发电机组惯性水平的简单加和。
在可再生能源发电比例不断增加的趋势下,对电力系统整体惯性水平进行准确评估,增强系统频率响应能力,保证系统安全稳定运行,是亟待解决的技术问题。
发明内容
为了解决现有技术的不足,本发明提供了一种基于广域信息的电力系统惯性水平评估方法,本发明根据慢同调理论对电力系统同步发电机组进行分群,选取各机群中起主导作用的机组,确定主导机组相关参变量与有功功率差额之间的关系,利用各主导机组的广域信息,得到系统总有功功率差额的估计值,评估电力系统整体惯性水平。
一种基于广域信息的电力系统惯性水平评估方法,包括以下步骤:
步骤1:根据慢同调理论,将电力系统同步发电机组分为多个慢同调机群;
步骤2:对每个慢同调机群,计算其包含的同步发电机组的参与因子和相关因子,选择起主导作用的机组,获得该机组的惯性时间常数;
步骤3:根据慢同调机群主导机组转子运动方程,确定主导机组电角速度、电角速度变化率、惯性时间常数与有功功率差额之间的关系;
步骤4:根据主导机组有功功率差额与惯性时间常数等参数间的关系,利用各慢同调机群主导机组的广域信息,包括电角速度、电角速度变化率,得到系统总有功功率差额的估计值,用来评估电力系统整体惯性水平。
进一步的,所述步骤1,根据慢同调理论,将电力系统同步发电机组分为多个慢同调机群的具体步骤为:
步骤1.1:根据同步发电机二阶经典模型,构造收缩到同步发电机内电势节点的导纳矩阵,并在此基础上求取系统状态矩阵;
步骤1.2:根据步骤1.1中所求系统状态矩阵,求其特征值与特征向量,并在此基础上确定慢同调分群数,选取慢同调分群的参考机组,计算分群矩阵,得到慢同调机群。
进一步的,所述步骤1.1中,构造收缩到同步发电机内电势节点的导纳矩阵的具体步骤为:
在不计入同步发电机内电势节点时,设系统共计有m+n个节点,其中,前m个节点为同步发电机端节点,后n个节点为其他节点,系统节点导纳矩阵可以写为:
其中,YA、YB、YC、YD分别为同步发电机端节点的自导纳矩阵、同步发电机端节点与其他节点的互导纳矩阵、其他节点与同步发电机端节点的互导纳矩阵、其他节点的自导纳矩阵。
计入同步发电机内电势节点,并将负荷表示为等值导纳的形式,形成增广节点导纳矩阵为:
其中,Ym为由同步发电机的等值导纳构成的对角子矩阵,Ym的对角线元素为Ymi=1/jx'd,x'd为同步发电机暂态电抗;YLm与YLn分别为由同步发电机端节点所连接负荷和其他节点所连接负荷的等值导纳构成的对角子矩阵,定义为:
其中
Pli+jQli为第i个节点的负荷功率,Uli为第i个节点的电压幅值;
为求出收缩到同步发电机内电势节点的导纳矩阵,将前述增广节点导纳矩阵表示为:
其中,Y1=Ym,Y2=[-Ym 0],
对应的网络方程为:
其中,为同步发电机内电势节点的注入电流向量,为同步发电机内电势节点的电压向量,为网络节点电压向量。
仅保留m个同步发电机内电势节点,由得到等值导纳矩阵:
YR=Y1-Y2Y4 -1Y3
其中,YR即为收缩到同步发电机内电势节点的导纳矩阵,导纳矩阵元素为Yij=Gij+jBij i,j=1,2,…,m,其中,当i=j时,Yij=Yii表示矩阵YR的自导纳元素,当i≠j时,Yij表示矩阵YR的互导纳元素,Gij、Bij分别为相应导纳元素的电导和电纳。
进一步的,所述步骤1.1中,求取系统状态矩阵的步骤为:
忽略阻尼常数,将同步发电机转子运动方程在其功率及角度稳态值附近线性化,可得计入同步发电机惯性时间常数后的系统状态矩阵A的表达式为:
A=T-1K
其中,T为同步发电机惯性时间常数的对角阵,K为系统状态矩阵,矩阵元素为:
其中,Ei'和E'j分别表示第i台和第j台同步发电机暂态电抗后的电势相量和的幅值,δi和δj分别为第i台和第j台同步发电机的转子电角度,Yij=Gij+jBij为收缩到同步发电机内电势节点的导纳矩阵元素,m为同步发电机的台数。
进一步的,所述步骤1.2中,确定慢同调分群数及其参考机组的具体步骤为:
求取状态矩阵A的m个特征值λ1,λ2,…,λm,按照特征值绝对值的大小进行排序,依次为|λ1|,|λ2|,…,|λm|,令
则r即为慢同调分群数目,并可获得慢模式组σr={λ1,λ2,…,λr}以及对应的右特征向量矩阵Ur=[u1,u2,…,ur],其中,ui满足方程Aui=λiui i=1,2,…r。
对右特征向量矩阵Ur采用列主元高斯消去法,可得r组线性无关的向量,对应的同步发电机定为慢同调分群的参考机组,并将Ur中对应的r行提取形成参考机组向量空间U1。
进一步的,所述步骤1.2中,获取分群矩阵并得到慢同调分群的具体步骤为:
将右特征向量矩阵Ur提取参考机组向量空间U1后剩下的m-r行特征向量提取形成矩阵U2,即U=[U1U2]T,定义分群矩阵为:L=U2U1 -1,根据分群矩阵L中每行的最大值,将剩余同步发电机组归类到对应的参考机组中,形成慢同调机群。
进一步的,所述步骤2,计算同步发电机组参与因子和相关因子并选取主导机组的具体步骤为:
vi和ui分别是步骤1.1中状态矩阵的特征值λi的左特征向量和右特征向量,则相对应的参与因子Pi和相关因子Cki定分别定义为:
Pi=vi Tui
Cki=|vkiuki/(vi Tui)|
vki和uki分别为vi和ui的第k个元素,参与因子Pi表征状态模式中变量的相对参与程度,相关因子Cki值的大小反映了状态量对相应特征值的能观性和能控性。
进一步的,所述步骤3,确定主导机组电角速度、电角速度变化率、惯性时间常数与有功功率差额之间的关系的具体步骤为:
慢同调机群中主导机组的转子运动方程为:
其中,J为转子转动惯量,单位为kg·m2;Ω是转子机械运动角速度,单位为rad/s;M为作用在转子上所有转矩的代数和;
主导机组的转子机械惯性时间常数可利用其额定转速下存储的动能与额定视在功率的比值来定义:
Ω0为同步发电机转子额定转速,SN为额定容量。
同步发电机电角速度ω与机械转速Ω之间的关系为:
ω=pΩ
其中,p为同步发电机极对数。
综合上述三式,考虑到同步发电机有功功率等于其转速与转矩的乘积,从而将同步发电机电角速度、电角速度变化率与惯性时间常数的关系表示为:
其中,ω*为同步发电机电角速度标幺值,ΔP为同步发电机有功功率差额。
进一步的,所述步骤4,根据主导机组有功功率差额与惯性时间常数等参数间的关系,利用各慢同调机群主导机组的广域信息,得到电力系统总有功功率差额的估计值,用来评估电力系统整体惯性水平具体步骤为:
采集、汇总各慢同调机群主导机组的电角速度、电角速度变化率信息,以全球定位系统的时钟信号作为统一的时间基准,形成反映系统整体动态特性的广域信息;
根据同步发电机有功功率差额与惯性时间常数等参数间的关系,利用各慢同调机群主导机组的广域信息,包括电角速度、电角速度变化率,将第i个慢同调机群所对应的有功功率差额ΔPi表示为:
其中,r为慢同调机群数,ki为第i个慢同调机群中同步发电机的数目,即TJij为第i个慢同调机群中第j个同步发电机的惯性时间常数,SNij为其额定容量;ωis*为第i个慢同调机群中主导机组的电角速度标幺值;
通过在主导机组上设置广域同步测量装置,测量电力系统中r个慢同调机群中主导机组的电角速度变化率,则利用每个慢同调机群所对应的有功功率差额来求取系统总有功功率差额估计值
电力系统中常用惯性时间常数与额定容量的乘积表征系统的惯性水平G,即:
G=TJSN
则电力系统中所有慢同调机群的总体惯性水平为:
进一步的,在实际电力系统中,如果实际的系统总有功功率差额ΔP已知,系统所有慢同调机群的总体惯性水平Ggen也已知,则计算出系统的整体惯性水平Gtotal:
与现有技术相比,本发明的有益效果是:
本发明提出的一种基于广域信息的电力系统惯性水平评估方法,利用基于统一时间基准的广域信息,能够反映同步发电机之外其他因素对系统整体惯性水平的影响,保证了系统整体惯性水平评估的准确性。根据慢同调理论对电力系统进行分群,减少了惯性水平评估中所需获取的信息量,显著提高了计算效率。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1是本发明方法的流程图;
图2是本发明实施例一New England 39节点系统慢同调分群结果图。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
正如背景技术所介绍的,现有技术中存在电力系统的惯性估计上不精确的不足,为了解决如上的技术问题,本申请提出了一种基于广域信息的电力系统惯性水平评估方法。
本申请的一种典型的实施方式中,如图1所示,提供了一种基于广域信息的电力系统惯性水平评估方法。
在本实施例中,New England 39节点系统包含10台同步发电机、39个节点以及19个负荷。
采用本发明提出的一种基于广域信息的电力系统惯性评估方法,评估NewEngland 39节点系统的惯性水平。该方法根据慢同调理论对电力系统同步发电机组进行分群,选取各机群中起主导作用的机组,确定主导机组相关参变量与有功功率差额之间的关系,利用各主导机组的广域信息,得到系统总有功功率差额的估计值,评估电力系统整体惯性水平。具体步骤包括:
步骤1:根据慢同调理论,将电力系统同步发电机组分为多个慢同调机群;
步骤2:对每个慢同调机群,计算其包含的同步发电机组的参与因子和相关因子,选择起主导作用的机组,获得该机组的惯性时间常数;
步骤3:根据慢同调机群主导机组转子运动方程,确定主导机组电角速度、电角速度变化率、惯性时间常数与有功功率差额之间的关系;
步骤4:利用各慢同调机群主导机组的电角速度、电角速度变化率等广域信息,结合系统中各台同步发电机的惯性时间常数与额定容量等数据,得到系统总有功功率差额的估计值,评估电力系统整体惯性水平。
进一步的,所述步骤1,根据慢同调理论,将电力系统同步发电机组分为多个慢同调机群的具体步骤为:
步骤1.1:根据发电机二阶经典模型,构造收缩到同步发电机内电势节点的导纳矩阵,并在此基础上求取系统状态矩阵;
步骤1.2:根据步骤1.1中所求系统状态矩阵,求其特征值与特征向量,并在此基础上确定慢同调分群数,选取慢同调分群的参考机组,计算分群矩阵,得到慢同调机群。
进一步的,所述步骤1.1中,构造收缩到同步发电机内电势节点的导纳矩阵的具体步骤为:
在不计入同步发电机内电势节点时,设系统共计有m+n个节点,其中,前m个节点为同步发电机端节点,后n个节点为其他节点,系统节点导纳矩阵可以写为:
其中,YA、YB、YC、YD分别为同步发电机端节点的自导纳矩阵、同步发电机端节点与其他节点的互导纳矩阵、其他节点与同步发电机端节点的互导纳矩阵、其他节点的自导纳矩阵。
计入同步发电机内电势节点,并将负荷表示为等值导纳的形式,形成增广节点导纳矩阵为:
其中,Ym为由同步发电机的等值导纳构成的对角子矩阵,Ym的对角线元素为Ymi=1/jx'd,x'd为同步发电机暂态电抗;YLm与YLn分别为由同步发电机端节点所连接负荷和其他节点所连接负荷的等值导纳构成的对角子矩阵,定义为:
其中
Pli+jQli为第i个节点的负荷功率,Uli为第i个节点的电压幅值;
为求出收缩到同步发电机内电势节点的导纳矩阵,将前述增广节点导纳矩阵表示为:
其中,Y1=Ym,Y2=[-Ym 0],
对应的网络方程为:
其中,为同步发电机内电势节点的注入电流向量,为同步发电机内电势节点的电压向量,为网络节点电压向量。
仅保留m个同步发电机内电势节点,由得到等值导纳矩阵:
YR=Y1-Y2Y4 -1Y3
其中,YR即为收缩到同步发电机内电势节点的导纳矩阵,导纳矩阵元素为Yij=Gij+jBij i,j=1,2,…,m,其中,当i=j时,Yij=Yii表示矩阵YR的自导纳元素,当i≠j时,Yij表示矩阵YR的互导纳元素,Gij、Bij分别为相应导纳元素的电导和电纳。
进一步的,所述步骤1.1中,求取系统状态矩阵的步骤为:
忽略阻尼常数,将同步发电机转子运动方程在其功率及角度稳态值附近线性化,可得计入同步发电机惯性时间常数后的系统状态矩阵A的表达式为:
A=T-1K
其中,T为同步发电机惯性时间常数的对角阵,K为系统状态矩阵,矩阵元素为:
其中,Ei'和E'j分别表示第i台和第j台同步发电机暂态电抗后的电势相量和的幅值,δi和δj分别为第i台和第j台同步发电机的转子电角度,Yij=Gij+jBij为收缩到同步发电机内电势节点的导纳矩阵元素,m为同步发电机的台数。
进一步的,所述步骤1.2中,确定慢同调分群数及其参考机组的具体步骤为:
求取状态矩阵A的m个特征值λ1,λ2,…,λm,按照特征值绝对值的大小进行排序,依次为|λ1|,|λ2|,…,|λm|,令
则r即为慢同调分群数目,并可获得慢模式组σr={λ1,λ2,…,λr}以及对应的右特征向量矩阵Ur=[u1,u2,…,ur],其中,ui满足方程Aui=λiui i=1,2,…,r。
进一步的,所述步骤1.2中,获取分群矩阵并得到慢同调分群的具体步骤为:
将右特征向量矩阵Ur提取参考机组向量空间U1后剩下的m-r行特征向量提取形成矩阵U2,即U=[U1U2]T。定义分群矩阵为:L=U2U1-1。根据分群矩阵L中每行的最大值,将剩余同步发电机组归类到对应的参考机组中,形成慢同调机群。
本实施例中,New England 39节点系统包含10个同步发电机节点,通过求取系统状态矩阵,识别出本系统的慢同调分群数为3个,参考机组为1、5、9号同步发电机,进一步获得分群矩阵如表1所示。
表1New England 39节点系统慢同调分群矩阵
根据分群矩阵中每行的最大值,将剩余同步发电机组归类到对应的参考机组中,可知New England 39节点系统中,同步发电机组1、2、3为慢同调机群1;同步发电机组4、5、6、7为慢同调机群2;同步发电机组8、9、10为慢同调机群3。系统慢同调分群结果如图2所示。
进一步的,所述步骤2,计算同步发电机组参与因子和相关因子并选取主导机组的具体步骤为:
vi和ui分别是步骤1.1中状态矩阵的特征值λi的左特征向量和右特征向量,则相对应的参与因子Pi和相关因子Cki定分别定义为:
Pi=vi Tui
Cki=|vkiuki/(vi Tui)|
vki和uki分别为vi和ui的第k个元素。
本实施例中,通过比较3个慢同调机群中各同步发电机组的参与因子和相关因子,选取出各慢同调机群的主导机组,可得慢同调机群1中的主导机组为1号同步发电机组,慢同调机群2中的主导机组为4号同步发电机组,慢同调机群3中的主导机组为8号同步发电机组。
进一步的,所述步骤3,确定主导机组电角速度、电角速度变化率、惯性时间常数与有功功率差额之间的关系的具体步骤为:
慢同调机群中主导机组的转子运动方程为:
其中,J为转子转动惯量,单位为kg·m2;Ω是转子机械运动角速度,单位为rad/s;M为作用在转子上所有转矩的代数和。
主导机组的转子机械惯性时间常数可利用其额定转速下存储的动能与额定视在功率的比值来定义:
Ω0为同步发电机转子额定转速,SN为额定容量。
同步发电机电角速度ω与机械转速Ω之间的关系为:
ω=pΩ
其中,p为同步发电机极对数。
综合上述三式,考虑到同步发电机有功功率等于其转速与转矩的乘积,从而将同步发电机电角速度、电角速度变化率与惯性时间常数的关系表示为:
其中,ω*为同步发电机电角速度标幺值,ΔP为同步发电机有功功率差额。
进一步的,所述步骤4,利用各慢同调机群主导机组的广域信息,结合同步发电机的参数,得到系统有功功率差额的估计值,评估电力系统整体惯性水平的具体步骤为:
采集、汇总各慢同调机群主导机组的电角速度、电角速度变化率信息,以全球定位系统的时钟信号作为统一的时间基准,形成反映系统整体动态特性的广域信息。
根据同步发电机有功功率差额与惯性时间常数等参数间的关系,利用各慢同调机群主导机组的电角速度、电角速度变化率等广域信息,可将第i个慢同调机群所对应的有功功率差额ΔPi表示为:
其中,r为慢同调机群数,ki为第i个慢同调机群中同步发电机的数目,即TJij为第i个慢同调机群中第j个同步发电机的惯性时间常数,SNij为其额定容量;ωis*为第i个慢同调机群中主导机组的电角速度标幺值。
通过在主导机组上设置广域同步测量装置,系统中r个慢同调机群中主导机组的电角速度变化率都可以准确的测得,则可以利用每个慢同调机群所对应的有功功率差额来求取电力系统总有功功率差额估计值
电力系统中常用惯性时间常数与额定容量的乘积表征系统的惯性水平G,即:G=TJSN,则电力系统中所有慢同调机群的总体惯性水平为:
在实际电力系统中,如果实际的系统总有功功率差额ΔP已知,系统所有慢同调机群的总体惯性水平Ggen也已知,则可以计算出系统的整体惯性水平Gtotal:
本实施例中,采集、汇总3个慢同调机群主导机组(1号、4号、8号同步发电机组)的电角速度、电角速度变化率信息,以全球定位系统的时钟信号作为统一的时间基准,形成反映系统整体动态特性的广域信息,并根据上述步骤,分别在2组故障情况下,对NewEngland39节点系统的整体惯性水平进行评估,结果如表2所示。
表2New England 39节点系统惯性水平评估结果
表2中,N表示需要提取电角速度、电角速度变化率等信息的同步发电机组的数目,%Gres表示系统整体的惯性水平与同步发电机组惯性水平之间的偏差百分比,即同步发电机之外其他因素对系统整体惯性水平的影响占比。由表2可知,采用本发明所提出的评估方法能够反映同步发电机之外其他因素对系统整体惯性水平的影响,保证了系统整体惯性水平评估的准确性;根据慢同调理论对电力系统进行分群,减少了惯性水平评估中所需获取的信息量,显著提高了计算效率。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (10)
1.一种基于广域信息的电力系统惯性水平评估方法,其特征是,包括以下步骤:
步骤1:根据慢同调理论,将电力系统同步发电机组分为多个慢同调机群;
步骤2:对每个慢同调机群,计算其包含的同步发电机组的参与因子和相关因子,选择起主导作用的机组,获得该机组的惯性时间常数;
步骤3:根据慢同调机群主导机组转子运动方程,确定主导机组电角速度、电角速度变化率、惯性时间常数与有功功率差额之间的关系;
步骤4:根据主导机组有功功率差额与惯性时间常数参数间的关系,利用各慢同调机群主导机组的广域信息,包括电角速度、电角速度变化率,得到电力系统总有功功率差额的估计值,用来评估电力系统整体惯性水平。
2.如权利要求1所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,所述步骤1,根据慢同调理论,将电力系统同步发电机组分为多个慢同调机群的具体步骤为:
步骤1.1:根据同步发电机二阶经典模型,构造收缩到同步发电机内电势节点的导纳矩阵,并在此基础上求取系统状态矩阵;
步骤1.2:根据步骤1.1中所求系统状态矩阵,求其特征值与特征向量,并在此基础上确定慢同调分群数,选取慢同调分群的参考机组,计算分群矩阵,得到慢同调机群。
3.如权利要求2所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,所述步骤1.1中,构造收缩到同步发电机内电势节点的导纳矩阵的具体步骤为:
在不计入同步发电机内电势节点时,设系统共计有m+n个节点,其中,前m个节点为同步发电机端节点,后n个节点为其他节点,系统节点导纳矩阵可以写为:
其中,YA、YB、YC、YD分别为同步发电机端节点的自导纳矩阵、同步发电机端节点与其他节点的互导纳矩阵、其他节点与同步发电机端节点的互导纳矩阵、其他节点的自导纳矩阵;
计入同步发电机内电势节点,并将负荷表示为等值导纳的形式,形成增广节点导纳矩阵为:
其中,Ym为由同步发电机的等值导纳构成的对角子矩阵,Ym的对角线元素为Ymi=1/jx′d,x′d为同步发电机暂态电抗;YLm与YLn分别为由同步发电机节点所连接负荷和其他节点所连接负荷的等值导纳构成的对角子矩阵,定义为:
其中
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>jQ</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msup>
<msub>
<mi>U</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
Pli+jQli为第i个节点的负荷功率,Uli为第i个节点的电压幅值;
为求出收缩到同步发电机内电势节点的导纳矩阵,将前述增广节点导纳矩阵表示为:
其中,Y1=Ym,Y2=[-Ym 0],
对应的网络方程为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>4</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,为同步发电机内电势节点的注入电流向量,为同步发电机内电势节点的电压向量,为网络节点电压向量;
仅保留m个同步发电机内电势节点,由得到等值导纳矩阵:
YR=Y1-Y2Y4 -1Y3
其中,YR即为收缩到同步发电机内电势节点的导纳矩阵,导纳矩阵元素为Yij=Gij+jBiji,j=1,2,…,m,其中,当i=j时,Yij=Yii表示矩阵YR的自导纳元素,当i≠j时,Yij表示矩阵YR的互导纳元素,Gij、Bij分别是相应的导纳元素的电导和电纳。
4.如权利要求2所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,所述步骤1.1中,求取系统状态矩阵的步骤为:
忽略阻尼常数,将同步发电机转子运动方程在其功率及角度稳态值附近线性化,可得计入同步发电机惯性时间常数后的系统状态矩阵A的表达式为:
A=T-1K
其中,T为同步发电机惯性时间常数的对角阵,K为系统状态矩阵,矩阵元素为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<msubsup>
<mi>E</mi>
<mi>j</mi>
<mo>&prime;</mo>
</msubsup>
<mo>&lsqb;</mo>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&delta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&delta;</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&delta;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&delta;</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>,</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>K</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,E′i和E′j分别表示第i台和第j台同步发电机暂态电抗后的电势相量和的幅值,δi和δj分别为第i台和第j台同步发电机的转子电角度,Yij=Gij+jBij为收缩到同步发电机内电势节点的导纳矩阵元素,m为同步发电机的台数。
5.如权利要求2所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,所述步骤1.2中,确定慢同调分群数及其参考机组的具体步骤为:
求取状态矩阵A的m个特征值λ1,λ2,…,λm,按照特征值绝对值的大小进行排序,依次为|λ1|,|λ2|,…,|λm|,令
<mrow>
<mo>|</mo>
<msub>
<mi>&lambda;</mi>
<mi>r</mi>
</msub>
<mo>|</mo>
<mo>/</mo>
<mo>|</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<mo>=</mo>
<munder>
<mi>min</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munder>
<mo>|</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>/</mo>
<mo>|</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
则r即为慢同调分群数目,并可获得慢模式组σr={λ1,λ2,…,λr}以及对应的右特征向量矩阵Ur=[u1,u2,…,ur],其中,ui满足方程Aui=λiui i=1,2,…,r;
对右特征向量矩阵Ur采用列主元高斯消去法,可得r组线性无关的向量,对应的同步发电机定为慢同调分群的参考机组,并将Ur中对应的r行提取形成参考机组向量空间U1。
6.如权利要求5所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,所述步骤1.2中,获取分群矩阵并得到慢同调分群的具体步骤为:
将右特征向量矩阵Ur提取参考机组向量空间U1后剩下的m-r行特征向量提取形成矩阵U2,即U=[U1 U2]T,定义分群矩阵为:L=U2U1 -1,根据分群矩阵L中每行的最大值,将剩余同步发电机组归类到对应的参考机组中,形成慢同调机群。
7.如权利要求5所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,所述步骤2,计算同步发电机组参与因子和相关因子并选取主导机组的具体步骤为:
vi和ui分别是步骤1.1中状态矩阵的特征值λi的左特征向量和右特征向量,则相对应的参与因子Pi和相关因子Cki定分别定义为:
Pi=vi Tui
Cki=|vkiuki/(vi Tui)|
vki和uki分别为vi和ui的第k个元素,参与因子Pi表征状态模式中变量的相对参与程度,相关因子Cki值的大小反映了状态量对相应特征值的能观性和能控性。
8.如权利要求1所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,所述步骤3,确定主导机组电角速度、电角速度变化率、惯性时间常数与有功功率差额之间的关系的具体步骤为:
慢同调机群中主导机组的转子运动方程为:
<mrow>
<mi>J</mi>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&Omega;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>M</mi>
</mrow>
其中,J为转子转动惯量,单位为kg·m2;Ω是转子机械运动角速度,单位为rad/s;M为作用在转子上所有转矩的代数和;
主导机组的转子机械惯性时间常数可利用其额定转速下存储的动能与额定视在功率的比值来定义:
<mrow>
<msub>
<mi>T</mi>
<mi>J</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>J&Omega;</mi>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
<msub>
<mi>S</mi>
<mi>N</mi>
</msub>
</mfrac>
</mrow>
Ω0为同步发电机转子额定转速,SN为额定容量;
同步发电机电角速度ω与机械转速Ω之间的关系为:
ω=pΩ
其中,p为同步发电机极对数;
综合上述三式,考虑到同步发电机有功功率等于其转速与转矩的乘积,从而将同步发电机电角速度、电角速度变化率与惯性时间常数的关系表示为:
<mrow>
<mfrac>
<mrow>
<msub>
<mi>d&omega;</mi>
<mo>*</mo>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>P</mi>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mi>J</mi>
</msub>
<msub>
<mi>S</mi>
<mi>N</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mo>*</mo>
</msub>
</mrow>
</mfrac>
</mrow>
其中,ω*为同步发电机电角速度标幺值,ΔP为同步发电机有功功率差额。
9.如权利要求1所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,所述步骤4,根据主导机组有功功率差额与惯性时间常数参数间的关系,利用各慢同调机群主导机组的广域信息,得到系统总有功功率差额的估计值,用来评估电力系统整体惯性水平具体步骤为:
采集、汇总各慢同调机群主导机组的电角速度、电角速度变化率信息,以全球定位系统的时钟信号作为统一的时间基准,形成反映系统整体动态特性的广域信息;
根据同步发电机有功功率差额与惯性时间常数参数间的关系,利用各慢同调机群主导机组的广域信息,包括电角速度、电角速度变化率,将第i个慢同调机群所对应的有功功率差额ΔPi表示为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
</munderover>
<msub>
<mi>T</mi>
<mrow>
<mi>J</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>N</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>i</mi>
<mi>s</mi>
<mo>*</mo>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>d&omega;</mi>
<mrow>
<mi>i</mi>
<mi>s</mi>
<mo>*</mo>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>r</mi>
<mo>;</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,r为慢同调机群数,ki为第i个慢同调机群中同步发电机的数目,即TJij为第i个慢同调机群中第j个同步发电机的惯性时间常数,SNij为其额定容量;ωis*为第i个慢同调机群中主导机组的电角速度标幺值;
通过在主导机组上设置广域同步测量装置,测量电力系统中r个慢同调机群中主导机组的电角速度变化率,则利用每个慢同调机群所对应的有功功率差额来求取系统总有功功率差额估计值
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mover>
<mi>P</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>r</mi>
</munderover>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>r</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
电力系统中常用惯性时间常数与额定容量的乘积表征系统的惯性水平G,即:
G=TJSN
则电力系统中所有慢同调机群的总体惯性水平为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>r</mi>
</munderover>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
</munderover>
<msub>
<mi>T</mi>
<mrow>
<mi>J</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>N</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>r</mi>
<mo>;</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>.</mo>
</mrow>
10.如权利要求9所述的一种基于广域信息的电力系统惯性水平评估方法,其特征是,在实际电力系统中,如果实际的系统总有功功率差额ΔP已知,系统所有慢同调机群的总体惯性水平Ggen也已知,则计算出系统的整体惯性水平Gtotal:
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>t</mi>
<mi>o</mi>
<mi>t</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>P</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mover>
<mi>P</mi>
<mo>^</mo>
</mover>
</mrow>
</mfrac>
<mo>&times;</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>g</mi>
<mi>e</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711364766.8A CN107947172B (zh) | 2017-12-18 | 2017-12-18 | 一种基于广域信息的电力系统惯性水平评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711364766.8A CN107947172B (zh) | 2017-12-18 | 2017-12-18 | 一种基于广域信息的电力系统惯性水平评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107947172A true CN107947172A (zh) | 2018-04-20 |
CN107947172B CN107947172B (zh) | 2023-11-21 |
Family
ID=61943716
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711364766.8A Active CN107947172B (zh) | 2017-12-18 | 2017-12-18 | 一种基于广域信息的电力系统惯性水平评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107947172B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108767876A (zh) * | 2018-05-31 | 2018-11-06 | 大连理工大学 | 一种面向大电网的主动频率响应模型预测控制方法 |
CN110943451A (zh) * | 2019-12-12 | 2020-03-31 | 华北电力大学 | 一种基于扰动数据的系统分区域惯量评估的方法 |
CN111293685A (zh) * | 2020-02-28 | 2020-06-16 | 华北电力大学 | 一种基于同调识别的系统分区惯量的评估方法 |
CN111371125A (zh) * | 2020-01-17 | 2020-07-03 | 武汉大学 | 一种考虑风机接入下提升系统同调性的解列分群优化方法 |
CN112329241A (zh) * | 2020-11-06 | 2021-02-05 | 国网新疆电力有限公司 | 一种计及负荷惯量的电力系统惯量评估方法 |
CN112421650A (zh) * | 2020-09-29 | 2021-02-26 | 中国电力科学研究院有限公司 | 一种评估电网惯性的方法和系统 |
CN112865076A (zh) * | 2021-01-12 | 2021-05-28 | 天津大学 | 面向电力电子化柔性配电网的改进增广节点方程建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104795812A (zh) * | 2015-04-07 | 2015-07-22 | 中国南方电网有限责任公司电网技术研究中心 | 一种考虑定子阻尼效应并采用变量代换的慢同调分区法 |
CN105281326A (zh) * | 2015-10-15 | 2016-01-27 | 中国南方电网有限责任公司电网技术研究中心 | 考虑定子阻尼影响的慢同调分区方法 |
CN106450528A (zh) * | 2016-11-23 | 2017-02-22 | 许继电源有限公司 | 储能系统及其功率均衡控制方法和控制装置 |
CN107017635A (zh) * | 2017-06-09 | 2017-08-04 | 广东电网有限责任公司电力调度控制中心 | 一种自适应解列控制方法及系统 |
-
2017
- 2017-12-18 CN CN201711364766.8A patent/CN107947172B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104795812A (zh) * | 2015-04-07 | 2015-07-22 | 中国南方电网有限责任公司电网技术研究中心 | 一种考虑定子阻尼效应并采用变量代换的慢同调分区法 |
CN105281326A (zh) * | 2015-10-15 | 2016-01-27 | 中国南方电网有限责任公司电网技术研究中心 | 考虑定子阻尼影响的慢同调分区方法 |
CN106450528A (zh) * | 2016-11-23 | 2017-02-22 | 许继电源有限公司 | 储能系统及其功率均衡控制方法和控制装置 |
CN107017635A (zh) * | 2017-06-09 | 2017-08-04 | 广东电网有限责任公司电力调度控制中心 | 一种自适应解列控制方法及系统 |
Non-Patent Citations (5)
Title |
---|
夏道止等: "第六章 同步电机的数学模型", 《电力系统分析》 * |
姜涛等: "基于广域量测信息相关性的电力系统同调辨识", 《电工技术学报》 * |
宋洪磊: "基于WAMS信息的大区域互联电网主动解列控制策略研究", 《中国博士学位论文全文数据库 工程科技II辑》 * |
林琳: "基于暂态能量函数法的电力系统安全性评估", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
马振玢: "低惯量电力系统的模式分析与同调识别研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108767876B (zh) * | 2018-05-31 | 2021-09-17 | 大连理工大学 | 一种面向大电网的主动频率响应模型预测控制方法 |
CN108767876A (zh) * | 2018-05-31 | 2018-11-06 | 大连理工大学 | 一种面向大电网的主动频率响应模型预测控制方法 |
CN110943451B (zh) * | 2019-12-12 | 2021-05-25 | 华北电力大学 | 一种基于扰动数据的系统分区域惯量评估的方法 |
CN110943451A (zh) * | 2019-12-12 | 2020-03-31 | 华北电力大学 | 一种基于扰动数据的系统分区域惯量评估的方法 |
CN111371125A (zh) * | 2020-01-17 | 2020-07-03 | 武汉大学 | 一种考虑风机接入下提升系统同调性的解列分群优化方法 |
CN111371125B (zh) * | 2020-01-17 | 2021-12-07 | 武汉大学 | 一种考虑风机接入下提升系统同调性的解列分群优化方法 |
CN111293685A (zh) * | 2020-02-28 | 2020-06-16 | 华北电力大学 | 一种基于同调识别的系统分区惯量的评估方法 |
CN112421650A (zh) * | 2020-09-29 | 2021-02-26 | 中国电力科学研究院有限公司 | 一种评估电网惯性的方法和系统 |
CN112421650B (zh) * | 2020-09-29 | 2022-07-01 | 中国电力科学研究院有限公司 | 一种评估电网惯性的方法和系统 |
CN112329241A (zh) * | 2020-11-06 | 2021-02-05 | 国网新疆电力有限公司 | 一种计及负荷惯量的电力系统惯量评估方法 |
CN112329241B (zh) * | 2020-11-06 | 2023-06-20 | 国网新疆电力有限公司 | 一种计及负荷惯量的电力系统惯量评估方法 |
CN112865076A (zh) * | 2021-01-12 | 2021-05-28 | 天津大学 | 面向电力电子化柔性配电网的改进增广节点方程建模方法 |
CN112865076B (zh) * | 2021-01-12 | 2022-05-20 | 天津大学 | 面向电力电子化柔性配电网的改进增广节点方程建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107947172B (zh) | 2023-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107947172A (zh) | 一种基于广域信息的电力系统惯性水平评估方法 | |
CN109217362B (zh) | 一种双馈风机并网系统低频振荡扰动源定位系统及方法 | |
CN102146812B (zh) | 电力系统原动机及其调速器实测建模方法 | |
Vittal et al. | Analysis of the inter-area mode phenomenon in power systems following large disturbances | |
Huang et al. | Generator dynamic model validation and parameter calibration using phasor measurements at the point of connection | |
CN110943451B (zh) | 一种基于扰动数据的系统分区域惯量评估的方法 | |
CN103618492B (zh) | 一种基于时频变换的同步发电机参数辨识方法 | |
CN103336909B (zh) | 一种风电接入电网的低频振荡辨识方法 | |
US11714934B2 (en) | Aggregated model of large-scale wind farms for power system simulation software tools | |
CN106526384B (zh) | 一种用于大规模电力系统的振荡源定位方法 | |
Zhang et al. | Synchrophasor measurement-based wind plant inertia estimation | |
CN110120686A (zh) | 一种基于电力系统在线惯量估计的新能源承载力预警方法 | |
CN109861258B (zh) | 一种虚拟同步机一次调频性能在线评价方法 | |
CN106058922A (zh) | 一种含虚拟惯量控制的双馈风电机组降阶仿真系统及方法 | |
CN106383270A (zh) | 基于广域测量信息的电力系统次同步振荡监测方法及系统 | |
CN104820741A (zh) | 兼顾风场分散性与机组差异性的风电场动态等值方法 | |
CN106410825A (zh) | 预测多端直流系统接入后电力系统机电振荡模式的方法 | |
CN113629729B (zh) | 基于频率测量点选取的含风电电力系统区域惯量估计方法 | |
CN110518632A (zh) | 一种风电场并网对电网惯量削弱的定量计算方法 | |
CN105244887A (zh) | 电力系统暂态功角失稳的闭环控制方法 | |
CN112448410B (zh) | 应用于含光伏电站电力系统低频振荡的模态综合分析方法 | |
Prévost et al. | An emulator for static and dynamic performance evaluation of small wind turbines | |
Jeffries et al. | Development of a dynamic model for no storage wind/diesel systems | |
CN106026182A (zh) | 考虑参数随机性的双馈风电机组仿真系统及其搭建方法 | |
Liu et al. | Inertia Characteristic Analysis and Inertia Estimation of New Energy Power System |
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 |