CN105891625B - 基于能量流的电力系统次同步振荡扰动源辨别方法 - Google Patents

基于能量流的电力系统次同步振荡扰动源辨别方法 Download PDF

Info

Publication number
CN105891625B
CN105891625B CN201610053096.7A CN201610053096A CN105891625B CN 105891625 B CN105891625 B CN 105891625B CN 201610053096 A CN201610053096 A CN 201610053096A CN 105891625 B CN105891625 B CN 105891625B
Authority
CN
China
Prior art keywords
value
frequency
energy flow
current
voltage
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
CN201610053096.7A
Other languages
English (en)
Other versions
CN105891625A (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201610053096.7A priority Critical patent/CN105891625B/zh
Publication of CN105891625A publication Critical patent/CN105891625A/zh
Application granted granted Critical
Publication of CN105891625B publication Critical patent/CN105891625B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

本发明提供一种基于能量流的电力系统次同步振荡扰动源辨别方法，包括：在待测设备工作状态下，在预设时间段内采集待测设备的三相电压值和电流值；对采集的电流值中的任一相电流进行识别获取基频、次同步频率和超同步频率，以及扭振频率；对采集的电压值和电流值进行变换获得XY坐标下的电压值和电流值；根据XY坐标下的电压值、电流值和扭振频率获得扭振频率的电压分量值和电流分量值；根据电压分量值和电流分量值获取能量流并计算能流功率；根据能流功率值的正负值确定扰动源和设备阻尼状态。本发明利用流入设备的能量流辨别扰动源并评估设备的阻尼特性。

Description

基于能量流的电力系统次同步振荡扰动源辨别方法
技术领域
本发明涉及电力系统分析技术领域，尤其涉及一种基于能量流的电力系统次同步振荡扰动源辨别方法。
背景技术
次同步振荡是威胁电力系统安全运行的一个突出问题。次同步振荡是指电力系统中频率低于工频但高于低频振荡频率(0.1～2.5Hz)的振荡。早期的次同步振荡主要由串联电容引发。近年来，由电力电子装置如高压直流、风电场引发的次同步振荡屡有发生，严重影响发电机和电网的安全运行。系统发生次同步振荡后，需要快速准确的辨别出扰动源，以采取针对性的措施平息振荡，而且振荡事件平息后的事后分析和设备缺陷检查，也需要查找扰动源并评估各个设备的阻尼特性。现有的次同步振荡监测分析方法还不能满足上述需求。
发明内容
本发明提供一种基于能量流的电力系统次同步振荡扰动源辨别方法，用于解决现有技术中不能准确找到设备扰动源和评估阻尼状态的问题。
为了解决上述问题，本发明提供一种基于能量流的电力系统次同步振荡扰动源辨别与阻尼评估方法，包括：
在待测设备工作状态下，在预设时间段内采集待测设备的三相电压值和三相电流值；
对采集的三相电流值中的任一相电流采用傅里叶变换进行识别获取基频、次同步频率和超同步频率，并根据所述基频、次同步频率和超同步频率获取扭振频率；
对采集的三相电压值和三相电流值采用变换矩阵进行变换获得XY坐标下的电压值和电流值；
根据XY坐标下的电压值、电流值和扭振频率采用滤波算法获得扭振频率的电压分量值和电流分量值；
根据电压分量值和电流分量值采用能量流计算公式获取能量流，并根据能量流获得能流功率；
根据能流功率值的正负值确定扰动源。
优选地，还包括根据能流功率值的正负值确定待测设备的阻尼状态。
优选地，所述扭振频率为所述基频与所述次同步频率之差，或为所述超同步频率与所述基频之差。
优选地，所述变换矩阵为：
其中，
f0为基频，θ为预设系数。
优选地，利用变换矩阵获得XY坐标下的电压值和电流值具体为：
其中，Ux,Uy,Ix,Iy为XY坐标下的电压值和电流值，ua, ub,uc为三相电压值，ia,ib,ic为三相电流值。
优选地，所述能量流计算公式为：其中，W(t)为时刻t的能量流，ts为初始时刻，初始值W(ts)＝0。
优选地，当在预设时间段内采集的数据为离散数据且采样周期为h时，所述能流流计算公式为：
其中，Wk为能 量流，i表示对应变量在t＝ts+ih处的值，W0＝W(ts)＝0。
优选地，通过线性拟合获取能量流与时间坐标上线性关系为W(t)＝αt+β，其中，α为能流功率值，β为系数。
优选地，当能流功率值为负值时，确定待测设备为扰动源，具有负阻尼状态；当能流功率值为正值时，确定待测设备具有正阻尼状态。
由上述技术方案可知，本发明的电力系统次同步振荡扰动源辨别方法，利用流入设备的能量流辨别扰动源并评估设备的阻尼特性。当能流功率为正，设备具有正阻尼，当能量功率为负，设备具有负阻尼，且是次同步振荡的扰动源。本方法提高了电网的次同步振荡的分析技术和控制技术水平，一旦电网发生次同步振荡，电网运行人员能够及时的辨别出扰动源，从而采取有效的控制措施快速平息振荡，此外对于事后分析，需要重点检查阻尼为负的设备，提高了设备缺陷检查的针对性。
附图说明
图1为本发明实施例提供的基于能量流的电力系统次同步振荡扰动源辨别方法的流程示意图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
图1示出了本发明实施例提供的一种基于能量流的电力系统次同步振荡扰动源辨别方法，包括：
S1、在待测设备工作状态下，在预设时间段内采集待测设备的三相电压值和三相电流值。在本步骤中，在预设时间段内采集的三相电压值和三相电流值实际上是关于预设时间t的函数。所述待测设备为发电机、串联电容、风电场等多种设备。在本发明实施例中，可以对所有待测的设备进行一一检测，以达到辨别目的。在本实施例中，所述预设时间过长会影响本方法的实时性，过短会影响分析结果的可靠性。因此所述预设时间可为1-5s，但本发明所述方案中不局限于该时间长度。当采用离散数据方式对电压值和电流值进行采集，在预设时间段内可确定一时间点ts，并以该时间点为起点按预设周期h确定其他时间点。即t＝ts+ih，其中，t为后续采样点处的时间点，i为采样周期的个数。
S2、对采集的三相电流值中的任一相电流采用傅里叶变换进行识别获取基频f0、次同步频率f-和超同步频率f+，并根据所述基频、次同步频率和超同步频率获取扭振频率f1。在本步骤中，对每个时间点采集的三相电压值和三相电流值定为一组三相电压值和三相电流值，即每组电压值和电流值与各自的时间点一一对应。所述扭振频率f1为所述基频与所述次同步频率之差，或为所述超同步频率与所述基频之差，即f1＝f0-f-或f1＝f+-f0。
S3、对采集的三相电压值和三相电流值采用变换矩阵进行变换获得XY坐标下的电压值和电流值。在本步骤中，需要将每个时间点采集的三相电压值和三相电流值转换成XY坐标下的电压值和电流值，以便后续步骤使用。因此，需要建立变换矩阵，如下：
其中，
f0为基频，θ为预设系数。
利用上述变换矩阵D获得XY坐标下的电压值和电流值具体为：
其中，Ux,Uy,Ix,Iy为XY坐标下的电压值和电流值，ua, ub,uc为三相电压值，ia,ib,ic为三相电流值。
S4、根据XY坐标下的电压值、电流值和扭振频率采用滤波算法获得扭振频率的电压分量值和电流分量值。在本步骤中，根据步骤S3获得的电压值(Ux,Uy)和电流值(Ix,Iy)，以及步骤S2获得的扭振频率采用滤波算法获得扭振频率f1的电压分量值Ux1,Uy1和电流分量值Ix1,Iy1。需要说明的是，所述滤波算法可采用任何可适应且能够达到本步骤目的的算法。其中，优选采用傅里叶变换算法。
S5、根据电压分量值和电流分量值采用能量流计算公式获取能量流，并根据能量流获得能流功率。在本步骤中，所述能量流计算公式如下：
其中，W(t)为时刻t的能量流，ts为初始时刻，初始值W(ts)＝0。
对于上述能量流计算公式中的Ux,Uy,Ix,Iy实际上是关于预设时间t的函数。为此，需要说明的是，对于步骤S1-S4中提到的三相电压值(Ux,Uy)、三相电流值(Ix,Iy)、基频f0、次同步频率f-和超同步频率f+均是建立在预设时间段内有关时间t的函数。为此，依据上述能量流计算公式可以得到在预设时间段内的能量流，然后根据能量流与时间的关系计算获得能流功率。
因为电力系统中一般数据采集是离散的。所述能量流计算公式如下：
其中，Wk为能 量流，Ux1,Uy1,Ix1,Iy1为扭振频率在XY坐标下的电压分量值和电流分量值。
计算得到的能量流是时间t的函数，对该函数进行线性拟合，得到能量流和时间的近似线性关系：W(t)＝αt+β，其中，α为待测设备的能流功率值，β为系数。
S6、据能流功率值的正负值确定扰动源和阻尼状态。在本步骤中，当能流功率值为负值时，表示该待测设备产生能流，可确定待测设备为扰动源，具有负阻尼状态。当能流功率值为正值时，表示该待测设备消耗能流，确定待测设备具有正阻尼状态。
本发明的电力系统次同步振荡扰动源辨别方法，利用流入设备的能量流辨别扰动源并评估设备的阻尼特性。当能流功率为正，设备具有正阻尼，当能量功率为负，设备具有负阻尼，且是次同步振荡的扰动源。本方法提高了电网的次同步振荡的分析技术和控制技术水平，一旦电网发生次同步振荡，电网运行人员能够及时的辨别出扰动源，从而采取有效的控制措施快速平息振荡，此外对于事后分析，需要重点检查阻尼为负的设备，提高了设备缺陷检查的针对性。
此外，本领域的技术人员能够理解，尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如，在下面的权利要求书中，所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制，并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中，不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。
本领域普通技术人员可以理解：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。

Claims (8)

1.一种基于能量流的电力系统次同步振荡扰动源辨别方法，其特征在于，包括：
在待测设备工作状态下，在预设时间段内采集待测设备的三相电压值和三相电流值；
对采集的三相电流值中的任一相电流采用傅里叶变换进行识别获取基频、次同步频率和超同步频率，并根据所述基频、次同步频率和超同步频率获取扭振频率；
对采集的三相电压值和三相电流值采用变换矩阵进行变换获得XY坐标下的电压值和电流值；
根据XY坐标下的电压值、电流值和扭振频率采用滤波算法获得扭振频率的电压分量值和电流分量值；
根据电压分量值和电流分量值采用能量流计算公式获取能量流，并根据能量流获得能流功率；
根据能流功率值的正负值确定扰动源；
其中，所述能量流计算公式为：其中，W(t)为时刻t的能量流，ts为初始时刻，初始值W(ts)＝0，Ux1,Uy1为扭振频率f1的电压分量值，Ix1,Iy1为扭振频率f1的电流分量值Ix1,Iy1，t为预设时间段内采样的时间点。
2.根据权利要求1所述的方法，其特征在于，还包括根据能流功率值的正负值确定待测设备的阻尼状态。
3.根据权利要求1所述的方法，其特征在于，所述扭振频率为所述基频与所述次同步频率之差，或为所述超同步频率与所述基频之差。
4.根据权利要求1所述的方法，其特征在于，所述变换矩阵为：其中，
f0为基频，θ为预设系数，t为预设时间段内采样的时间点。
5.根据权利要求4所述的方法，其特征在于，利用变换矩阵获得XY坐标下的电压值和电流值具体为：
其中，Ux,Uy,Ix,Iy为XY坐标下的电压值和电流值，ua,ub,uc为三相电压值，ia,ib,ic为三相电流值。
6.根据权利要求1所述的方法，其特征在于，当在预设时间段内采集的数据为离散数据且采样周期为h时，所述能量流计算公式为：
其中，Wk为能量流，i表示对应变量在t＝ts+ih处的值，W0＝W(ts)＝0，Ux1,Uy1,Ix1,Iy1为扭振频率在XY坐标下的电压分量值和电流分量值，t为预设时间段内采样的时间点，i为采样周期的个数。
7.根据权利要求1或6所述的方法，其特征在于，通过线性拟合获取能量流与时间坐标上线性关系为W(t)＝αt+β，其中，α为能流功率值，β为系数，t为预设时间段内采样的时间点。
8.根据权利要求1所述的方法，其特征在于，当能流功率值为负值时，确定待测设备为扰动源，具有负阻尼状态；当能流功率值为正值时，确定待测设备具有正阻尼状态。
CN201610053096.7A 2016-01-26 2016-01-26 基于能量流的电力系统次同步振荡扰动源辨别方法 Active CN105891625B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610053096.7A CN105891625B (zh) 2016-01-26 2016-01-26 基于能量流的电力系统次同步振荡扰动源辨别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610053096.7A CN105891625B (zh) 2016-01-26 2016-01-26 基于能量流的电力系统次同步振荡扰动源辨别方法

Publications (2)

Publication Number Publication Date
CN105891625A CN105891625A (zh) 2016-08-24
CN105891625B true CN105891625B (zh) 2018-10-30

Family

ID=57013713

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610053096.7A Active CN105891625B (zh) 2016-01-26 2016-01-26 基于能量流的电力系统次同步振荡扰动源辨别方法

Country Status (1)

Country Link
CN (1) CN105891625B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106526384B (zh) * 2016-12-01 2018-12-25 西南交通大学 一种用于大规模电力系统的振荡源定位方法
CN108303609B (zh) * 2018-01-19 2024-02-02 国家电网公司 一种次同步振荡在线监测及扰动源定位方法及系统
CN111679125B (zh) * 2020-06-04 2021-07-20 北京交通大学 一种电力系统振荡辨识的方法和装置
CN113866501B (zh) * 2021-09-22 2024-01-02 南京南瑞继保电气有限公司 一种基于能流功率的次同步振荡扰动源定位方法
CN117543627B (zh) * 2024-01-08 2024-04-02 华北电力大学 一种双馈风机振荡扰动源定位方法、系统及电子设备

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4355241A (en) * 1980-08-15 1982-10-19 Electric Power Research Institute, Inc. Means for damping subsynchronous oscillations in an AC power system including overload protection
CN102122823B (zh) * 2011-02-23 2013-08-14 南京南瑞继保电气有限公司 电力系统振荡扰动源定位的方法
CA2844731A1 (en) * 2011-09-12 2013-03-21 Alstom Technology Ltd Sub-synchronous oscillation damping by shunt facts apparatus
CN103117547A (zh) * 2013-03-04 2013-05-22 中国南方电网有限责任公司 一种电力系统低频振荡的振荡源定位方法
CN103208808B (zh) * 2013-03-07 2015-01-14 武汉大学 一种电力系统次同步振荡模态辨识方法
CN105223418B (zh) * 2015-09-22 2018-12-11 清华大学 次同步和超同步谐波相量的测量方法及测量装置
CN105226679A (zh) * 2015-09-29 2016-01-06 特变电工新疆新能源股份有限公司 一种在线动态抑制电网低频振荡和次同步振荡的控制方法
CN105119286B (zh) * 2015-10-08 2017-05-17 南京南瑞继保电气有限公司 一种次同步振荡源定位方法及装置

Also Published As

Publication number Publication date
CN105891625A (zh) 2016-08-24

Similar Documents

Publication Publication Date Title
CN105891625B (zh) 基于能量流的电力系统次同步振荡扰动源辨别方法
Li et al. Application of bandwidth EMD and adaptive multiscale morphology analysis for incipient fault diagnosis of rolling bearings
Gargoom et al. Automatic classification and characterization of power quality events
JP5607772B2 (ja) 太陽電池パネル監視プログラム、太陽電池パネル監視装置及び太陽電池パネル監視方法
CN108802535B (zh) 筛选方法、主干扰源识别方法及装置、服务器及存储介质
US11658490B2 (en) Wind energy system and method for identifying low-frequency oscillations in an electrical supply network
JP2022515935A (ja) 誘導電動機の軸受故障重大度を推定する方法
CN113285471B (zh) 海上风电场次超同步振荡源感知与定位方法、装置及设备
Ma et al. Application of multisynchrosqueezing transform for subsynchronous oscillation detection using PMU data
US20210281071A1 (en) Method and apparatus for detecting low-frequency oscillations
CN107545253B (zh) 一种电力系统低频振荡监测方法
CN113238142A (zh) 用于集成电路的方法和系统
CN117471227B (zh) 汽车线束参数性能测试方法及测试系统
CN104091056A (zh) 一种实验室设备使用寿命预测系统和方法
CN103926494B (zh) 一种干扰源的确定方法及装置
CN107844670A (zh) 一种谐波统计所需样本数量的计算方法
CN105954695A (zh) 一种基于同步的同质传感器突变参数识别方法与装置
CN105068035B (zh) 一种电压互感器误差水平动态检测方法及系统
CN107798149B (zh) 一种飞机维修性评估方法
CN109613372B (zh) 一种基于多元电网数据库的电网故障诊断方法
CN105203865B (zh) 分布式电源并网逆变器及电容的在线测试系统的工作方法
JP2011202956A (ja) 避雷装置の故障判定方法
Abidullah et al. Real-time power quality disturbances detection and classification system
CN103928923B (zh) 一种基于灵敏度分析的网络稳态电能质量预警方法
Eidson et al. An evaluation of the extent of correlation between interharmonic and voltage fluctuation measurements

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant