CN106772042B - 一种发电机原动系统阻尼在线检测方法 - Google Patents
一种发电机原动系统阻尼在线检测方法 Download PDFInfo
- Publication number
- CN106772042B CN106772042B CN201611209173.XA CN201611209173A CN106772042B CN 106772042 B CN106772042 B CN 106772042B CN 201611209173 A CN201611209173 A CN 201611209173A CN 106772042 B CN106772042 B CN 106772042B
- Authority
- CN
- China
- Prior art keywords
- generator
- energy stream
- dissipation energy
- excitation
- damping
- 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
Links
- 238000013016 damping Methods 0.000 title claims abstract description 56
- 238000010998 test method Methods 0.000 title description 5
- 230000005284 excitation Effects 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 18
- 230000010355 oscillation Effects 0.000 claims description 16
- 230000001419 dependent effect Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
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
- G01R31/34—Testing dynamo-electric machines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明提供一种发电机原动系统阻尼在线测量方法，包括以下步骤：S1、计算电网流入发电机的发电机耗散能量流和单位阻尼耗散能量流；S2、计算发电机流入励磁系统的励磁耗散能量流；S3、根据电网流入发电机的耗散能量流和励磁耗散能量流之差计算出发电机流入原动系统的原动耗散能量流；S4、对原动耗散能量流和单位阻尼耗散能量流分别进行线性拟合，并计算两者斜率比，得到原动系统的阻尼转矩系数。通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的原动耗散能量流，并在线计算发电机原动系统的阻尼转矩系数。
Description
技术领域
本发明涉及电力系统分析技术领域，更具体地，涉及一种发电机原动系统阻尼在线检测方法。
背景技术
超低频振荡是威胁电力系统安全运行的一个突出问题，当发生超低频振荡时，系统内所有机组共同振荡，频率和功率也发生持续振荡，严重威胁电网安全。超低频振荡是一次调频过程中出现的稳定问题，与发电机的原动机和调速系统密切相关，原动系统产生负阻尼是导致超低频振荡的重要原因之一。目前没有对原动系统阻尼在线评估方法，使得超低频振荡的快速有效控制成为难题。本发明通过在线评估计算发电机原动系统的阻尼转矩系数，能够在线得到各发电机原动系统的阻尼情况，当系统发生超低频振荡时可以快速确定提供负阻尼的原动系统并将其调速器退出，使振荡快速平息，以满足电网安全需求。
发明内容
本发明提供一种克服上述问题或者至少部分地解决上述问题的发电机原动系统阻尼在线测量方法，通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的原动耗散能量流，并在线计算发电机原动系统的阻尼转矩系数。
根据本发明的一个方面，提供一种发电机原动系统阻尼在线测量方法，包括以下步骤：
S1、计算电网流入发电机的发电机耗散能量流和单位阻尼耗散能量流；
S2、计算发电机流入励磁系统的励磁耗散能量流；
S3、根据电网流入发电机的耗散能量流和励磁耗散能量流之差计算发电机流入原动系统的原动耗散能量流；
S4、对原动耗散能量流和单位阻尼耗散能量流分别进行线性拟合，并计算两者斜率比，得到原动系统的阻尼转矩系数。
作为优选的，所述步骤S1具体包括：
S11、监测电力系统振荡后发电机的有功功率振幅，确定有功功率振幅稳定的一段时间窗口；
S12、获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f，并进行滤波预处理；
S13、计算该时间窗口内每个采样点电网流入发电机的耗散能量流和每个采样点单位阻尼耗散能量流。
作为优选的，所述步骤S12具体包括：获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f；计算电压U的自然对数lnU，对P、Q、lnU、f进行高通滤波，去除直流分量和低频分量，得到相应偏差△P、△Q、△lnU、△f。
作为优选的，所述步骤S13具体包括：
对时间窗口内的每个采集样点t，计算电网流入发电机的发电机耗散能量流：
式中，为时刻t电网流入发电机的的耗散能量流，ts为计算的起始时刻，初始值
对时间窗口内每个采样点t，计算单位阻尼耗散能量流：
式中，f0为额定频率，为时刻t单位阻尼消耗能量流，ts为计算的起始时刻，初始值
作为优选的，所述步骤S13进一步包括，当采样数据为离散数据时，采样周期为h，则：
电网流入发电机的发电机耗散能量流为：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
单位阻尼消耗能量流：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
作为优选的，所述步骤S2包括：
S21、采集该段时间窗口内发电机的励磁电压Ufd、励磁电流Ifd并进行滤波预处理，得到励磁电压、励磁电流的偏差量△Ufd、△Ifd；
S22、并计算计算机流入励磁系统的励磁耗散能量流。
作为优选的，所述步骤S22具体包括：
对于每个时间窗口内的每个采样点t，计算发电机流入励磁系统的励磁耗散能量流：
式中，为时刻t发电机流入励磁系统的励磁耗散能量流，ts为计算的起始时刻，初始值
当采集数据为离散数据时，采集周期为h，发电机流入励磁系统的励磁耗散能量流为：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
作为优选的，所述步骤S3具体包括：
根据电网流入发电机的耗散能量流和励磁耗散能量流得到原动耗散能量流：
当采集数据为离散数据时，原动耗散能量流为：
本申请提出一种发电机原动系统阻尼在线检测方法，通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的原动耗散能量流，并在线计算发电机原动系统的阻尼转矩系数，获取各个发电机原动系统的阻尼情况，当系统发生超低频振荡时可以快速确定提供负阻尼的原动系统，并将其一次调频退出，使振荡可以快速平息，保证电网的安全性。
附图说明
图1为本发明实施例中的原动系统阻尼在线测量方法流程图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
图1示出了一种发电机原动系统阻尼在线检测方法，利用能量流方法计算发电机流入原动系统的耗散能量流，进而通过耗散能量流分析原动系统的阻尼作用。由于计算发电机流入原动系统的耗散能量流公式中很多变量在实际电网中无法测量或测量精度较低，因此本发明采取了间接计算的方法，通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的耗散能量流。将得到的能量流进行线性拟合得到斜率即为流入原动系统的耗散功率，该功率与单位阻尼耗散功率的比值即为原动系统阻尼转矩系数，具体包括以下步骤：
S1、计算电网流入发电机的发电机耗散能量流和单位阻尼耗散能量流；
S2、计算发电机流入励磁系统的励磁耗散能量流；
S3、根据发电机耗散能量流和励磁耗散能量流之差计算发电机流入原动系统的原动耗散能量流；
S4、对原动耗散能量流和单位阻尼耗散能量流分别进行线性拟合，并计算两者斜率比，得到原动系统的阻尼转矩系数,该值为正代表原动系统提供正阻尼，否则为负阻尼。
作为优选的，所述步骤S1具体包括：
S11、监测电力系统振荡后发电机的有功功率振幅，确定有功功率振幅稳定的一段时间窗口；
S12、获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f，并进行滤波预处理；其中P、Q分别为发电机发出的有功功率和无功功率，f为发电机所连母线处的频率，U为发电机所连母线处电压。
S13、计算该时间窗口内每个采样点电网流入发电机的发电机耗散能量流和每个采样点单位阻尼耗散能量流。
作为优选的，所述步骤S12具体包括：获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f；计算电压U的自然对数lnU，对P、Q、lnU、f进行高通滤波，去除直流分量和低频分量，得到相应偏差△P、△Q、△lnU、△f。
作为优选的，所述步骤S13具体包括：
对时间窗口内的每个采集样点t，计算电网流入发电机的发电机耗散能量流：
式中，为时刻t电网流入发电机的的耗散能量流，ts为计算的起始时刻，初始值
当采样数据为离散数据时，采样周期为h，则：
电网流入发电机的发电机耗散能量流为：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
对时间窗口内每个采样点t，计算单位阻尼耗散能量流：
式中，f0为额定频率，为时刻t单位阻尼消耗能量流，ts为计算的起始时刻，初始值
当数据为离散数据时，采样周期为h时，单位阻尼消耗能量流：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
作为优选的，所述步骤S2包括：
S21、采集该段时间窗口内发电机的励磁电压Ufd、励磁电流Ifd并进行滤波预处理，得到励磁电压、励磁电流的偏差量△Ufd、△Ifd；Ufd、Ifd分别为励磁电压、流入发电机励磁绕组的励磁电流。
S22、并计算计算机流入励磁系统的励磁耗散能量流。
作为优选的，所述步骤S22具体包括：
对于每个时间窗口内的每个采样点t，计算发电机流入励磁系统的励磁耗散能量流：
式中，为时刻t发电机流入励磁系统的励磁耗散能量流，ts为计算的起始时刻，初始值
当采集数据为离散数据时，采集周期为h，发电机流入励磁系统的励磁耗散能量流为：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
作为优选的，所述步骤S3具体包括：
根据电网流入发电机的耗散能量流和励磁耗散能量流的得到原动耗散能量流：
当采集数据为离散数据时，原动耗散能量流为：
最后，对发电机流入原动系统的耗散能量流和单位阻尼耗散能量流分别进行线性拟合，得到拟合结果Wpm D≈kpmt+cpm，使Wpm D＝kpmt+cpm，计算两者斜率比值即为发电机原动系统的阻尼转矩系数。
本申请提出一种发电机原动系统阻尼在线检测方法，通过计算电网流入发电机和发电机流入励磁系统的耗散能量流来间接计算流入原动系统的原动耗散能量流，并在线计算发电机原动系统的阻尼转矩系数，获取各个发电机原动系统的阻尼情况，当系统发生超低频振荡时可以快速确定提供负阻尼的原动系统，并将其一次调频退出，使振荡可以快速平息，保证电网的安全性。
最后，本申请的方法仅为较佳的实施方案，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种发电机原动系统阻尼在线测量方法，其特征在于，包括以下步骤：
S1、计算电网流入发电机的发电机耗散能量流和单位阻尼耗散能量流；
S2、计算发电机流入励磁系统的励磁耗散能量流；
S3、根据电网流入发电机的发电机耗散能量流和励磁耗散能量流之差计算出发电机流入原动系统的原动耗散能量流；
S4、对原动耗散能量流和单位阻尼耗散能量流分别进行线性拟合，并计算两者斜率比，得到原动系统的阻尼转矩系数。
2.根据权利要求1所述的发电机原动系统阻尼在线测量方法，其特征在于，所述步骤S1具体包括：
S11、监测电力系统振荡后发电机的有功功率振幅，确定有功功率振幅稳定的一段时间窗口；
S12、获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f，并进行滤波预处理；
S13、计算该时间窗口内每个采样点电网流入发电机的发电机耗散能量流和每个采样点单位阻尼耗散能量流。
3.根据权利要求2所述的发电机原动系统阻尼在线测量方法，其特征在于，所述步骤S12具体包括：
获取该时间窗口内发电机的有功功率P、无功功率Q、电压U、频率f；
计算电压U的自然对数lnU，对P、Q、lnU、f进行高通滤波，去除直流分量和低频分量，得到相应偏差△P、△Q、△lnU、△f。
4.根据权利要求3所述的发电机原动系统阻尼在线测量方法，其特征在于，所述步骤S13具体包括：
对时间窗口内的每个采集样点t，计算电网流入发电机的发电机耗散能量流：
式中，为时刻t电网流入发电机的耗散能量流，ts为计算的起始时刻，初始值
对时间窗口内每个采样点t，计算单位阻尼耗散能量流：
式中，f0为额定频率，为时刻t单位阻尼消耗能量流，ts为计算的起始时刻，初始值
5.根据权利要求4所述的发电机原动系统阻尼在线测量方法，其特征在于，所述步骤S13进一步包括，当采样数据为离散数据时，采样周期为h，则：
电网流入发电机的发电机耗散能量流为：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
单位阻尼消耗能量流：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
6.根据权利要求5所述的发电机原动系统阻尼在线测量方法，其特征在于，所述步骤S2包括：
S21、采集该段时间窗口内发电机的励磁电压Ufd、励磁电流Ifd并进行滤波预处理，得到励磁电压、励磁电流的偏差量△Ufd、△Ifd；
S22、并计算计算机流入励磁系统的励磁耗散能量流。
7.根据权利要求6所述的发电机原动系统阻尼在线测量方法，其特征在于，所述步骤S22具体包括：
对于每个时间窗口内的每个采样点t，计算发电机流入励磁系统的励磁耗散能量流：
式中，为时刻t发电机流入励磁系统的励磁耗散能量流，ts为计算的起始时刻，初始值
当采集数据为离散数据时，采集周期为h，发电机流入励磁系统的励磁耗散能量流为：
式中，下标i表示对应变量在t＝ts+ih处的值，初始值
8.根据权利要求7所述的发电机原动系统阻尼在线测量方法，其特征在于，所述步骤S3具体包括：
根据电网流入发电机的发电机耗散能量流和励磁耗散能量流的得到原动耗散能量流：
当采集数据为离散数据时，原动耗散能量流为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611209173.XA CN106772042B (zh) | 2016-12-23 | 2016-12-23 | 一种发电机原动系统阻尼在线检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611209173.XA CN106772042B (zh) | 2016-12-23 | 2016-12-23 | 一种发电机原动系统阻尼在线检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106772042A CN106772042A (zh) | 2017-05-31 |
| CN106772042B true CN106772042B (zh) | 2019-06-18 |
Family
ID=58920232
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611209173.XA Active CN106772042B (zh) | 2016-12-23 | 2016-12-23 | 一种发电机原动系统阻尼在线检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106772042B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108896908B (zh) * | 2018-04-03 | 2021-07-09 | 上海电器科学研究所（集团）有限公司 | 电动机有限元振动计算模型中阻尼参数的精确识别方法 |
| CN110417054B (zh) * | 2019-07-30 | 2020-10-16 | 华北电力大学 | 一种双馈风电并网系统稳定性调整系统及方法 |
| CN110875599A (zh) * | 2019-11-04 | 2020-03-10 | 云南电网有限责任公司 | 一种电网频率振荡的控制方法及系统 |
| CN110879324B (zh) * | 2019-12-02 | 2022-05-17 | 国网湖南省电力有限公司 | 一种励磁系统在功率振荡中阻尼特性的辨别方法 |
| CN111555312B (zh) * | 2020-05-27 | 2021-06-22 | 四川大学 | 一种适用于电力系统超低频振荡稳定性评估的方法 |
| CN112861328B (zh) * | 2021-01-22 | 2022-08-30 | 东北电力大学 | 一种基于随机响应信号的发电机阻尼评估装置及方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201303308Y (zh) * | 2008-11-19 | 2009-09-02 | 湖南大学 | 一种发电机原动系统动态仿真的装置 |
| CN102944798B (zh) * | 2012-11-29 | 2015-05-20 | 武汉华中电力电网技术有限公司 | 一种负阻尼低频振荡与强迫功率振荡判别方法 |
| CN103117547A (zh) * | 2013-03-04 | 2013-05-22 | 中国南方电网有限责任公司 | 一种电力系统低频振荡的振荡源定位方法 |
| CN104111405B (zh) * | 2014-06-23 | 2016-10-05 | 华中科技大学 | 一种基于阻尼转矩分析法的电力系统低频振荡源定位方法 |
| CN106058897B (zh) * | 2016-07-28 | 2018-09-21 | 东南大学 | 一种基于相量的发电机强迫扰动源定位方法 |
-
2016
- 2016-12-23 CN CN201611209173.XA patent/CN106772042B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106772042A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106772042B (zh) | 一种发电机原动系统阻尼在线检测方法 | |
| CN102288283B (zh) | 汽轮发电机组振动故障的监测方法及装置 | |
| CN103036498B (zh) | 一种基于pmu的同步发电机实用模型参数校核与辨识方法 | |
| CN105606895A (zh) | 电力系统次同步振荡成份的在线检测及滤除方法 | |
| CN104865474B (zh) | 一种基于pmu数据实时监测低频振荡源的方法 | |
| EP3824523B1 (en) | Grid impedence estimation for power converter controller | |
| CN106948796A (zh) | 一种抽油机不停机间抽控制装置及方法 | |
| CN103344914B (zh) | 基于归一化的风电机组故障预警方法 | |
| CN103117547A (zh) | 一种电力系统低频振荡的振荡源定位方法 | |
| CN102928697B (zh) | 励磁调节器pss2a模型低频段阻尼检测方法及系统 | |
| CN113991702B (zh) | 一种基于准稳态数据的电力系统惯量评估方法 | |
| CN101915601B (zh) | 一种1000mw汽轮发电机组轴系模态阻尼的求取方法 | |
| CN109642542A (zh) | 用于风力涡轮机的诊断系统和方法 | |
| Jinbang et al. | Detection and reduction of middle frequency resonance for an industrial servo | |
| CN106058897B (zh) | 一种基于相量的发电机强迫扰动源定位方法 | |
| CN110927442A (zh) | 一种基于边缘计算的多模态振荡在线监测预警系统 | |
| CN105891625B (zh) | 基于能量流的电力系统次同步振荡扰动源辨别方法 | |
| CN107104448A (zh) | 电力系统稳定器及方法 | |
| CN104777396A (zh) | 一种综合考虑电机电流和电压的转子故障诊断方法 | |
| CN104101805B (zh) | 一种基于振荡能量注入的励磁系统负阻尼检测方法 | |
| Dan et al. | Monitoring, intelligent perception, and early warning of vortex‐induced vibration of suspension bridge | |
| CN110365028B (zh) | 一种一次调频影响功率振荡的判断方法 | |
| Guan et al. | Adaptive surge detection of magnetic suspension centrifugal blower based on rotor radial displacement signal and SOGI-FLL with prefilter | |
| DE102015218934A1 (de) | Überwachen einer sensorlosen Läuferlageerfassung | |
| CN116014760B (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 |