CN116203350A - 一种高水电占比电网中频率振荡源定位方法 - Google Patents

一种高水电占比电网中频率振荡源定位方法 Download PDF

Info

Publication number
CN116203350A
CN116203350A CN202310210299.2A CN202310210299A CN116203350A CN 116203350 A CN116203350 A CN 116203350A CN 202310210299 A CN202310210299 A CN 202310210299A CN 116203350 A CN116203350 A CN 116203350A
Authority
CN
China
Prior art keywords
frequency
generator
oscillation
time
positioning
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
Application number
CN202310210299.2A
Other languages
English (en)
Other versions
CN116203350B (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202310210299.2A priority Critical patent/CN116203350B/zh
Publication of CN116203350A publication Critical patent/CN116203350A/zh
Application granted granted Critical
Publication of CN116203350B publication Critical patent/CN116203350B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241—The oscillation concerning frequency
    • 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
    • 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/08—Locating faults in cables, transmission lines, or networks
    • G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • 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/08—Locating faults in cables, transmission lines, or networks
    • G01R31/088—Aspects of digital computing
    • 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]

Abstract

本发明公开了一种高水电占比电网中频率振荡源定位方法，属于电力系统稳定分析技术领域；本发明基于发电机频率信号和有功功率信号计算发电机转速偏差信号和机械功率偏差信号，并计算得到两者的同步压缩小波变换系数矩阵，后建立频率振荡时域能量耗散模型，利用交叉同步压缩小波变换系数矩阵计算频域耗散能量谱，后根据振荡源定位判据和耗散能量谱辨识频率振荡主导振荡频率，定位提供负阻尼转矩的发电机调速系统，评估各发电机调速系统对振荡的贡献程度，输出定位结果；相较于现有算法，本发明实现了超低频段范围内的频率振荡源定位，提高了振荡频率辨识精度和振荡源定位效率。

Description

一种高水电占比电网中频率振荡源定位方法
技术领域
本发明专利属于电力系统稳定分析技术领域，具体涉及一种高水电占比电网中频率振荡源定位方法。
背景技术
为解决电能生产和消耗的空间逆向分布问题，我国通过大容量、远距离输电实现跨区域电网互联。异步联网可以减弱区域间同步发电机组的耦合，增强区域电网暂态功角稳定性，如云南电网、西南电网均先后实现异步联网运行，传统功角稳定问题得到有效解决。但在此类高水电占比电网中，出现了多次电力系统一次调频过程中小扰动失稳的超低频（≤0.1Hz）频率振荡问题，其振荡持续时间长、影响范围广、频率失稳风险大。为找到诱发频率振荡的机组或设备以便于振荡后紧急抑制，亟需提出一种高水电占比电网中频率振荡源定位方法。
目前较为实用的振荡源定位分析方法主要为：基于模型动态仿真法与基于广域测量系统（WAMS）的暂态能量法。但基于模型动态仿真法依赖于电网精确建模和准确数据，对于实际电系统的复杂性往往无能为力。随着WAMS的快速发展，基于WAMS的能量流法在复杂电网中的振荡源定位得到广泛应用。但该类方法受限于定位单模式振荡源，而频段为超低频的频率振荡经常伴随低频振荡模式或其他频段分量，使用此类方法则需要对WAMS测量信号进行模式辨识、信号滤波、重构等预处理操作，过程较为繁琐，计算效率较低。
发明内容
有鉴于此，本发明专利基于SSWT（同步压缩小波变换），提出了一种高水电占比电网中频率振荡源定位方法。该方法在兼顾频率振荡信号含低频振荡分量或其他振荡分量时，改进使用进能量法时对WAMS测量信号预处理过程，实现高水电占比电网中频率振荡源准确、高效定位，为振荡发生后紧急抑制提供前提和参考，以保障电网频率稳定性能，具有较高的科学意义及实用价值。
本发明专利采用的技术方案如下：
一种高水电占比电网中频率振荡源定位方法，包括以下步骤：
步骤1：获取高水电占比电网发生频率振荡后各发电机频率信号和有功功率信号，基于各发电机频率信号和有功功率信号计算各发电机转速偏差信号和机械功率偏差信号，其中各发电机转速偏差信号和机械功率偏差信号的计算方法如下：
在频率振荡过程中，各发电机转速和频率同调振荡，其转速偏差信号
Figure SMS_1
可以通过机端频率信息计算得到：
Figure SMS_2
（1）
式中：
Figure SMS_3
为发电机机端频率；/>
Figure SMS_4
为频率参考值；/>
Figure SMS_5
单位为/>
Figure SMS_6
，/>
Figure SMS_7
单位为Hz；/>
由发电机的转子运动方程，忽略发电机阻尼系数D，可计算得到发电机机械功率偏差信号
Figure SMS_8
：
Figure SMS_9
（2）
式中，
Figure SMS_10
为发电机电磁功率，t为时间，/>
Figure SMS_11
为发电机惯性时间常数。
步骤2：基于常规信号同步压缩小波变换系数矩阵计算公式，分别计算得到各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵，具体包括以下步骤：
步骤2.1：常规信号同步压缩小波变换系数矩阵计算公式的获取方法如下：
对于常规信号
Figure SMS_12
，其中A为幅值，/>
Figure SMS_13
为角频率，t为时间，选择频率集中在正半轴的小波基函数/>
Figure SMS_14
，其傅里叶变换/>
Figure SMS_15
满足频率/>
Figure SMS_16
时，/>
Figure SMS_17
，该常规信号的连续小波变换系数为：
Figure SMS_18
（3）
式中：
Figure SMS_19
为s(t)的连续小波变换的系数，/>
Figure SMS_20
为尺度因子；b为时间平移因子，/>
Figure SMS_21
表示取共轭；
时间尺度平面上满足s(t)的连续小波变换系数
Figure SMS_22
的任意一点/>
Figure SMS_23
的瞬时频率为：
Figure SMS_24
(4)
通过对瞬时频率
Figure SMS_26
的计算，实现尺度-时间/>
Figure SMS_28
到频率-时间
Figure SMS_31
的映射，即同步压缩的过程，离散情况下，取尺度/>
Figure SMS_27
的第/>
Figure SMS_29
个步长为
Figure SMS_32
，频率/>
Figure SMS_34
第/>
Figure SMS_25
个步长为/>
Figure SMS_30
，定义在
Figure SMS_33
的同步压缩小波变换系数矩阵/>
Figure SMS_35
为：
Figure SMS_36
(5)
式中，
Figure SMS_39
为常规信号s(t)在尺度/>
Figure SMS_40
和时间b处的连续小波变换系数，
Figure SMS_42
为经同步压缩后常规信号s(t)在频率/>
Figure SMS_38
和时间b处的同步压缩小波变换系数矩阵，/>
Figure SMS_41
提高了能量的时频聚集性；需要说明的是，符号/>
Figure SMS_43
既可以表示发电机转速，也可以表示频率信号。这里符号重了，但是概念不一样：前面步骤1中/>
Figure SMS_44
是转速偏差信号，后面步骤2中加了个下标1,特指第/>
Figure SMS_37
个频率步长。
步骤2.2：基于公式（5）分别计算得到各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵
Figure SMS_45
和/>
Figure SMS_46
。
步骤3：建立频率振荡时域能量耗散模型，基于频率振荡时域能量耗散模型并结合各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵计算各发电机调速系统频率振荡的频域耗散能量谱，具体包括以下步骤：
步骤a：建立发电机调速系统频率振荡时域能量耗散模型：
Figure SMS_47
(6)
式中，
Figure SMS_48
为发电机调速系统所消耗的能量，/>
Figure SMS_49
为发电机机械功率偏差信号，/>
Figure SMS_50
为发电机转子角偏差信号，/>
Figure SMS_51
为发电机转速偏差信号，t为时间，/>
Figure SMS_52
，/>
Figure SMS_53
为频率参考值；
步骤b：对（6）式中两个实信号
Figure SMS_54
和/>
Figure SMS_55
，利用CWT和其逆变换可以得其Parseval恒等式：
Figure SMS_56
(7)
式中：
Figure SMS_57
为小波变换的容许常数，t为时间，/>
Figure SMS_58
和/>
Figure SMS_59
分别为转速偏差信号/>
Figure SMS_60
和机械功率偏差信号/>
Figure SMS_61
的在尺度/>
Figure SMS_62
和时间b处的CWT系数；
步骤c：定义发电机调速系统基于CWT的频域耗散能量为：
Figure SMS_63
(8)
式中，
Figure SMS_64
即为调速系统CWT时频域耗散能量，经过SSWT后，得到/>
Figure SMS_65
和/>
Figure SMS_66
的映射关系，则基于SSWT的时频耗散能量可表示为：
Figure SMS_67
(9）
式中，
Figure SMS_68
和/>
Figure SMS_69
分别为转速偏差信号/>
Figure SMS_70
和机械功率偏差信号
Figure SMS_71
的在频率/>
Figure SMS_72
和时间b处的SSWT变换系数矩阵，/>
Figure SMS_73
为调速系统的SSWT时频域耗散能量；
步骤d：对式（9）在时间所对应的位移系数b上积分后，可得关于频率
Figure SMS_74
的频域耗散能量表达式：
Figure SMS_75
（10）
式中，
Figure SMS_76
为调速系统频域耗散能量在频率/>
Figure SMS_77
点处的值，/>
Figure SMS_78
为
Figure SMS_79
和/>
Figure SMS_80
的交叉同步压缩小波变换系数矩阵对位移系数b积分后在频率/>
Figure SMS_81
处的值；
由互能量密度谱的概念可知：
Figure SMS_82
（11）/>
式中，Re表示取实部，提取式（11）中的被积项，定义为调速系统的耗散能量谱：
Figure SMS_83
（12）
式中，
Figure SMS_84
为发电机调速系统在频率/>
Figure SMS_85
处的SSWT耗散能量谱。
步骤4：基于频域耗散能量谱并结合振荡源定位依据辨识频率振荡主导振荡频率，定位提供负阻尼转矩的发电机调速系统，评估各发电机调速系统对振荡的贡献程度，输出定位结果，其中振荡源定位依据具体为：
Figure SMS_86
为波谷时，调速系统提供负阻尼转矩，为振荡源；/>
Figure SMS_87
为波峰时，调速系统提供正阻尼转矩，为非振荡源；此外，/>
Figure SMS_88
峰/谷值所在频率F点为频率振荡的主导振荡频率，/>
Figure SMS_89
为波谷时，其绝对值越大，对诱发振荡的贡献越大，其发电机调速系统是紧急抑制时的首要控制对象。
综上所述，由于采用了上述技术方案，本发明的有益效果是：
本发明在连续小波变换(CWT)的基础上使用SSWT（同步压缩小波变换），提高了振荡频率辨识分辨率，将时域暂态能量法推广至时频域；兼顾了频率振荡信号中混含低频信号或其他模式分量的问题，实现了超低频段范围内的频率振荡源定位，保证了振荡源定位精度；同时改进了对WAMS测量信息进行预处理过程，避免了对信号模式辨识、滤波、重构等过程，提高了计算和定位效率。
附图说明
本发明将通过例子并参照附图的方式说明，其中：
图1是本发明的高水电占比电网频率振荡源定位方法流程图；
图2是4机2区系统结构图；
图3是发电机转速偏差信号；
图4是发电机机械功率偏差信号；
图5是各发电机调速系统SSWT耗散能量谱；
图6是退出G1、G2调速系统一次调频后各发电机转速偏差信号；
图7是振荡发生后各发电机调速系统时域耗散能量流。
实施方式
为使本发明实施方式的目的、技术方案和优点更加清楚，下面将结合本发明实施方式中的附图，对本发明实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式是本发明一部分实施方式，而不是全部的实施方式。通常在此处附图中描述和示出的本发明实施方式的组件可以以各种不同的配置来布置和设计。
因此，以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。
需要说明的是，在不冲突的情况下，本发明中的实施方式及实施方式中的特征可以相互组合。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
需要说明的是，在不冲突的情况下，本发明中的实施方式及实施方式中的特征可以相互组合。
实施例
本发明实施例中公开了一种高水电占比电网中频率振荡源定位方法，包括以下步骤：
步骤1：获取高水电占比电网发生频率振荡后各发电机频率信号和有功功率信号，基于各发电机频率信号和有功功率信号计算各发电机转速偏差信号和机械功率偏差信号，其中各发电机转速偏差信号和机械功率偏差信号的计算方法如下：
在频率振荡过程中，各发电机转速和频率同调振荡，其转速偏差信号
Figure SMS_90
可以通过机端频率信息计算得到：
Figure SMS_91
（1）
式中：
Figure SMS_92
为发电机机端频率；/>
Figure SMS_93
为频率参考值；/>
Figure SMS_94
单位为/>
Figure SMS_95
，/>
Figure SMS_96
单位为Hz；
由发电机的转子运动方程，忽略发电机阻尼系数D，可计算得到发电机机械功率偏差信号
Figure SMS_97
：
Figure SMS_98
（2）
式中，
Figure SMS_99
为发电机电磁功率，t为时间，/>
Figure SMS_100
为发电机惯性时间常数。
步骤2：基于常规信号同步压缩小波变换系数矩阵计算公式，分别计算得到各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵，具体包括以下步骤：
步骤2.1：常规信号同步压缩小波变换系数矩阵计算公式的获取方法如下：
对于常规信号
Figure SMS_101
，其中A为幅值，/>
Figure SMS_102
为角频率，t为时间，选择频率集中在正半轴的小波基函数/>
Figure SMS_103
，其傅里叶变换/>
Figure SMS_104
满足频率/>
Figure SMS_105
时，/>
Figure SMS_106
，该常规信号的连续小波变换系数为：
Figure SMS_107
（3）
式中：
Figure SMS_108
为s(t)的连续小波变换的系数；/>
Figure SMS_109
为尺度因子；b为时间平移因子，/>
Figure SMS_110
表示取共轭；
时间尺度平面上满足s(t)的连续小波变换系数
Figure SMS_111
的任意一点/>
Figure SMS_112
的瞬时频率为：/>
Figure SMS_113
（4）
通过对瞬时频率
Figure SMS_116
的计算，实现尺度-时间/>
Figure SMS_119
到频率-时间
Figure SMS_122
的映射，即同步压缩的过程，离散情况下，取尺度/>
Figure SMS_115
的第/>
Figure SMS_118
个步长为
Figure SMS_121
，频率/>
Figure SMS_124
第/>
Figure SMS_114
个步长为/>
Figure SMS_117
，定义在
Figure SMS_120
的同步压缩小波变换系数矩阵/>
Figure SMS_123
为：
Figure SMS_125
（5）
式中，
Figure SMS_126
为常规信号s(t)在尺度/>
Figure SMS_127
和时间b处的连续小波变换系数，
Figure SMS_128
为经同步压缩后常规信号s(t)在频率/>
Figure SMS_129
和时间b处的同步压缩小波变换系数矩阵，/>
Figure SMS_130
提高了能量的时频聚集性；
步骤2.2：基于公式（5）分别计算得到各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵
Figure SMS_131
和/>
Figure SMS_132
。
步骤3：建立频率振荡时域能量耗散模型，基于频率振荡时域能量耗散模型并结合各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵计算各发电机调速系统频率振荡的频域耗散能量谱，具体包括以下步骤：
步骤a：建立发电机调速系统频率振荡时域能量耗散模型：
Figure SMS_133
（6）
式中，
Figure SMS_134
为发电机调速系统所消耗的能量，/>
Figure SMS_135
为发电机机械功率偏差信号，/>
Figure SMS_136
为发电机转子角偏差信号，/>
Figure SMS_137
为发电机转速偏差信号，t为时间，/>
Figure SMS_138
，/>
Figure SMS_139
为频率参考值；
步骤b：对（6）式中两个实信号
Figure SMS_140
和/>
Figure SMS_141
，利用CWT和其逆变换可以得其Parseval恒等式：
Figure SMS_142
（7）
式中：
Figure SMS_143
为小波变换的容许常数，t为时间，/>
Figure SMS_144
和/>
Figure SMS_145
分别为转速偏差信号/>
Figure SMS_146
和机械功率偏差信号/>
Figure SMS_147
的在尺度/>
Figure SMS_148
和时间b处的CWT系数；
步骤c：定义发电机调速系统基于CWT的频域耗散能量为：
Figure SMS_149
（8）
式中，
Figure SMS_150
即为调速系统CWT时频域耗散能量，经过SSWT后，得到/>
Figure SMS_151
和
Figure SMS_152
的映射关系，则基于SSWT的时频耗散能量可表示为：/>
Figure SMS_153
（9）
式中，
Figure SMS_154
和/>
Figure SMS_155
分别为转速偏差信号/>
Figure SMS_156
和机械功率偏差信号
Figure SMS_157
的在频率/>
Figure SMS_158
和时间b处的SSWT变换系数矩阵，/>
Figure SMS_159
为调速系统的SSWT时频域耗散能量；
步骤d：对式（9）在时间所对应的位移系数b上积分后，可得关于频率
Figure SMS_160
的频域耗散能量表达式：
Figure SMS_161
（10）
式中，
Figure SMS_162
调速系统频域耗散能量在频率/>
Figure SMS_163
点处的值，/>
Figure SMS_164
为/>
Figure SMS_165
和/>
Figure SMS_166
的交叉同步压缩小波变换系数矩阵对位移系数b积分后在频率/>
Figure SMS_167
处的值；
由互能量密度谱的概念可知：
Figure SMS_168
（11）
式中，Re表示取实部，提取式（11）中的被积项，定义为调速系统的耗散能量谱：
Figure SMS_169
（12）
式中，
Figure SMS_170
为发电机调速系统在频率/>
Figure SMS_171
处的SSWT耗散能量谱。
步骤4：基于频域耗散能量谱并结合振荡源定位依据辨识频率振荡主导振荡频率，定位提供负阻尼转矩的发电机调速系统，评估各发电机调速系统对振荡的贡献程度，输出定位结果，其中振荡源定位依据具体为：
Figure SMS_172
为波谷时，调速系统提供负阻尼转矩，为振荡源；/>
Figure SMS_173
为波峰时，调速系统提供正阻尼转矩，为非振荡源；此外，/>
Figure SMS_174
峰/谷值所在频率/>
Figure SMS_175
点为频率振荡的主导振荡频率，/>
Figure SMS_176
为波谷时，其绝对值越大，对诱发振荡的贡献越大，其发电机调速系统是紧急抑制时的首要控制对象。
为验证本发明方法的有效性，建立如图2所示典型4机2区系统，4台发电机均为水电，10s时在母线7-9之间220km线路上加入0.05s三相短路故障，设置各发电机调速系统参数，此时产生了全网频率同调的频率振荡，根据步骤1可得到如图3所示的各发电机转速偏差信号
Figure SMS_177
和图4所示的机械功率偏差信号/>
Figure SMS_178
。
根据步骤2、步骤3计算各发电机调速系统SSWT耗散能量谱，结果如图5所示，从图5可知，各发电机的SSWT耗散能量谱均在频率为0.0603Hz处出现峰/谷值，说明系统产生了的主导振荡频率为0.0603Hz的超低频频率振荡，为验证所提方法辨识系统主导振荡频率的准确性，对频率振荡数据分别利用传统Prony分析法、TLS-ESPRIT算法进行分析和对比，结果如表1所示：
Figure SMS_179
/>
根据图5的SSWT耗散能量谱可以看出，发电机G1、G2调速系统提供负阻尼转矩，是振荡源，且G1绝对值最大，对振荡贡献最多，其次为G2，为验证本发明所提定位方法的有效性，分别通过在110s时退出G1、G2一次调频和与传统时域能量流法对比，结果如图6和图7所示，从图6可以看出，退出振荡源G1、G2调速系统一次调频后，振荡得到有效抑制，且G1效果比G2更好，从图7可以看出G1、G2发出能量，提供负阻尼，总的来说，三者从定位结果上是一致的，本发明所提频率振荡源定位方法是有效的。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (6)

1.一种高水电占比电网中频率振荡源定位方法，其特征在于，包括以下步骤：
步骤1：获取高水电占比电网发生频率振荡后各发电机频率信号和有功功率信号，基于各发电机频率信号和有功功率信号计算各发电机转速偏差信号和机械功率偏差信号；
步骤2：基于常规信号同步压缩小波变换系数矩阵计算公式，分别计算得到各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵；
步骤3：建立频率振荡时域能量耗散模型，基于频率振荡时域能量耗散模型并结合各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵，计算各发电机调速系统频率振荡的频域耗散能量谱；
步骤4：基于频域耗散能量谱并结合振荡源定位依据辨识频率振荡主导振荡频率，定位提供负阻尼转矩的发电机调速系统，评估各发电机调速系统对振荡的贡献程度，输出定位结果。
2.根据权利要求1所述的一种高水电占比电网中频率振荡源定位方法，其特征在于，步骤1中各发电机转速偏差信号和机械功率偏差信号的计算方法如下：
Figure QLYQS_1
（1)
通过发电机频率计算得到转速偏差信号
Figure QLYQS_2
：
式中：
Figure QLYQS_3
为发电机机端频率；/>
Figure QLYQS_4
为频率参考值；/>
Figure QLYQS_5
单位为/>
Figure QLYQS_6
，/>
Figure QLYQS_7
单位为Hz；
基于发电机的转子运动方程，忽略发电机阻尼系数D，计算得到发电机机械功率偏差信号
Figure QLYQS_8
：
Figure QLYQS_9
(2)
式中，
Figure QLYQS_10
为发电机电磁功率，t为时间，/>
Figure QLYQS_11
为发电机惯性时间常数。
3.根据权利要求2所述的一种高水电占比电网中频率振荡源定位方法，其特征在于，步骤2具体包括以下步骤：
步骤2.1：常规信号同步压缩小波变换系数矩阵计算公式的获取方法如下：
对于常规信号
Figure QLYQS_12
，其中A为幅值，/>
Figure QLYQS_13
为角频率，t为时间，选择频率集中在正半轴的小波基函数/>
Figure QLYQS_14
，其傅里叶变换/>
Figure QLYQS_15
满足频率/>
Figure QLYQS_16
时，/>
Figure QLYQS_17
，该常规信号的连续小波变换系数为：
Figure QLYQS_18
(3)
式中：
Figure QLYQS_19
为s(t)的连续小波变换的系数，/>
Figure QLYQS_20
为尺度因子；b为时间平移因子，/>
Figure QLYQS_21
表示取共轭；
时间尺度平面上满足s(t)的连续小波变换系数
Figure QLYQS_22
的任意一点/>
Figure QLYQS_23
的瞬时频率为：
Figure QLYQS_24
(4)/>
通过对瞬时频率
Figure QLYQS_26
的计算，实现尺度-时间/>
Figure QLYQS_28
到频率-时间/>
Figure QLYQS_31
的映射，即同步压缩的过程，离散情况下，取尺度/>
Figure QLYQS_25
的第/>
Figure QLYQS_30
个步长为
Figure QLYQS_33
，频率/>
Figure QLYQS_35
第/>
Figure QLYQS_27
个步长为/>
Figure QLYQS_29
，定义在
Figure QLYQS_32
的同步压缩小波变换系数矩阵/>
Figure QLYQS_34
为：
Figure QLYQS_36
(5)
式中，
Figure QLYQS_37
为常规信号s(t)在尺度/>
Figure QLYQS_38
和时间b处的连续小波变换系数，/>
Figure QLYQS_39
为经同步压缩后常规信号s(t)在频率/>
Figure QLYQS_40
和时间b处的同步压缩小波变换系数矩阵；
步骤2.2：基于公式（5）分别计算得到各发电机转速偏差信号和机械功率偏差信号的同步压缩小波变换系数矩阵
Figure QLYQS_41
和/>
Figure QLYQS_42
。
4.根据权利要求3所述的一种高水电占比电网中频率振荡源定位方法，其特征在于，步骤3中所述建立频率振荡时域能量耗散模型具体包括以下步骤：
步骤a：建立发电机调速系统频率振荡时域能量耗散模型：
Figure QLYQS_43
(6)
式中，
Figure QLYQS_44
为发电机调速系统所消耗的能量，/>
Figure QLYQS_45
为发电机机械功率偏差信号，/>
Figure QLYQS_46
为发电机转子角偏差信号，/>
Figure QLYQS_47
为发电机转速偏差信号，t为时间，/>
Figure QLYQS_48
，/>
Figure QLYQS_49
为频率参考值；
步骤b：对（6）式中两个实信号
Figure QLYQS_50
和/>
Figure QLYQS_51
，利用连续小波变换和其逆变换可以得其Parseval恒等式：
Figure QLYQS_52
(7)
式中：
Figure QLYQS_53
为小波变换的容许常数，t为时间，/>
Figure QLYQS_54
和/>
Figure QLYQS_55
分别为转速偏差信号/>
Figure QLYQS_56
和机械功率偏差信号/>
Figure QLYQS_57
的在尺度/>
Figure QLYQS_58
和时间b处的连续小波变换系数。
5.根据权利要求4所述的一种高水电占比电网中频率振荡源定位方法，其特征在于，步骤3中所述利用同步压缩小波变换系数矩阵计算各发电机调速系统频率振荡的频域耗散能量谱具体包括以下步骤：
步骤c：基于公式（6）和公式（7），定义发电机调速系统基于连续小波变换的频域耗散能量为：
Figure QLYQS_59
(8)
式中，
Figure QLYQS_60
即为调速系统连续小波变换时频域耗散能量，经过同步压缩小波变换后，得到/>
Figure QLYQS_61
和/>
Figure QLYQS_62
的映射关系，则基于同步压缩小波变换的时频域耗散能量可表示为：
Figure QLYQS_63
(9)
式中，
Figure QLYQS_64
和/>
Figure QLYQS_65
分别为转速偏差信号/>
Figure QLYQS_66
和机械功率偏差信号
Figure QLYQS_67
的在频率/>
Figure QLYQS_68
和时间b处的同步压缩小波变换系数矩阵，/>
Figure QLYQS_69
为调速系统的同步压缩小波变换时频域耗散能量；
步骤d：对式（9）在时间所对应的位移系数b上积分后，可得关于频率
Figure QLYQS_70
的频域耗散能量表达式：
Figure QLYQS_71
(10)
式中，
Figure QLYQS_72
为调速系统频域耗散能量在频率/>
Figure QLYQS_73
点处的值，/>
Figure QLYQS_74
为/>
Figure QLYQS_75
和
Figure QLYQS_76
的交叉同步压缩小波变换系数矩阵对位移系数b积分后在频率/>
Figure QLYQS_77
处的值；
由互能量密度谱的概念可知：
Figure QLYQS_78
(11)
式中，Re表示取实部，提取式（11）中的被积项，定义为调速系统的耗散能量谱：
Figure QLYQS_79
(12)
式中，
Figure QLYQS_80
为发电机调速系统在频率/>
Figure QLYQS_81
处的同步压缩小波变换耗散能量谱。
6.根据权利要求5所述的一种高水电占比电网中频率振荡源定位方法，其特征在于，步骤4中所述振荡源定位依据具体为：
Figure QLYQS_82
为波谷时，调速系统提供负阻尼转矩，为振荡源；/>
Figure QLYQS_83
为波峰时，调速系统提供正阻尼转矩，为非振荡源；此外，/>
Figure QLYQS_84
峰/谷值所在频率/>
Figure QLYQS_85
点为频率振荡的主导振荡频率，/>
Figure QLYQS_86
为波谷时，其绝对值越大，对诱发振荡的贡献越大，其发电机调速系统是紧急抑制时的首要控制对象。/>
CN202310210299.2A 2023-03-07 2023-03-07 一种高水电占比电网中频率振荡源定位方法 Active CN116203350B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310210299.2A CN116203350B (zh) 2023-03-07 2023-03-07 一种高水电占比电网中频率振荡源定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310210299.2A CN116203350B (zh) 2023-03-07 2023-03-07 一种高水电占比电网中频率振荡源定位方法

Publications (2)

Publication Number Publication Date
CN116203350A true CN116203350A (zh) 2023-06-02
CN116203350B CN116203350B (zh) 2024-02-27

Family

ID=86515617

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310210299.2A Active CN116203350B (zh) 2023-03-07 2023-03-07 一种高水电占比电网中频率振荡源定位方法

Country Status (1)

Country Link
CN (1) CN116203350B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104111405A (zh) * 2014-06-23 2014-10-22 华中科技大学 一种基于阻尼转矩分析法的电力系统低频振荡源定位方法
CN106526384A (zh) * 2016-12-01 2017-03-22 西南交通大学 一种用于大规模电力系统的振荡源定位方法
WO2020007970A1 (de) * 2018-07-06 2020-01-09 Wobben Properties Gmbh Verfahren zum erkennen niederfrequenter schwingungen und erfassungseinrichtung dafür
CN110879355A (zh) * 2019-12-03 2020-03-13 国网湖南省电力有限公司 一种基于fft的发电机调速系统引发低频振荡辨识方法
KR20200081591A (ko) * 2018-12-27 2020-07-08 현대자동차주식회사 웨이블릿변환을 이용한 동기전동기의 고장진단방법
CN112069453A (zh) * 2020-09-08 2020-12-11 东北电力大学 一种基于小波变换的电力系统强迫振荡源时频域定位方法
US20210036518A1 (en) * 2019-07-30 2021-02-04 North China Electric Power University System and method of locating oscillation sources of wind power integrated system based on energy spectrums
CN114690034A (zh) * 2022-03-14 2022-07-01 东北电力大学 基于小波耗散能量谱的强迫振荡源定位方法
CN114707538A (zh) * 2022-03-14 2022-07-05 东北电力大学 一种电力系统强迫振荡源定位的耗散能量谱方法
CN115327374A (zh) * 2022-07-21 2022-11-11 东北电力大学 基于多通道同步压缩变换的强迫振荡源定位方法及装置
CN115411721A (zh) * 2022-07-21 2022-11-29 东北电力大学 多元变分模态分解的电力系统强迫振荡源定位方法及装置
US20220399719A1 (en) * 2021-06-02 2022-12-15 North China Electric Power University Oscillation active damping control method and system for grid-tied type-4 wind turbine generator

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104111405A (zh) * 2014-06-23 2014-10-22 华中科技大学 一种基于阻尼转矩分析法的电力系统低频振荡源定位方法
CN106526384A (zh) * 2016-12-01 2017-03-22 西南交通大学 一种用于大规模电力系统的振荡源定位方法
WO2020007970A1 (de) * 2018-07-06 2020-01-09 Wobben Properties Gmbh Verfahren zum erkennen niederfrequenter schwingungen und erfassungseinrichtung dafür
KR20200081591A (ko) * 2018-12-27 2020-07-08 현대자동차주식회사 웨이블릿변환을 이용한 동기전동기의 고장진단방법
US20210036518A1 (en) * 2019-07-30 2021-02-04 North China Electric Power University System and method of locating oscillation sources of wind power integrated system based on energy spectrums
CN110879355A (zh) * 2019-12-03 2020-03-13 国网湖南省电力有限公司 一种基于fft的发电机调速系统引发低频振荡辨识方法
CN112069453A (zh) * 2020-09-08 2020-12-11 东北电力大学 一种基于小波变换的电力系统强迫振荡源时频域定位方法
US20220399719A1 (en) * 2021-06-02 2022-12-15 North China Electric Power University Oscillation active damping control method and system for grid-tied type-4 wind turbine generator
CN114690034A (zh) * 2022-03-14 2022-07-01 东北电力大学 基于小波耗散能量谱的强迫振荡源定位方法
CN114707538A (zh) * 2022-03-14 2022-07-05 东北电力大学 一种电力系统强迫振荡源定位的耗散能量谱方法
CN115327374A (zh) * 2022-07-21 2022-11-11 东北电力大学 基于多通道同步压缩变换的强迫振荡源定位方法及装置
CN115411721A (zh) * 2022-07-21 2022-11-29 东北电力大学 多元变分模态分解的电力系统强迫振荡源定位方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
RUI YU等: "《Mechanism and Damping Torque Analysis of Power System Stabilizer for Suppressing Ultra-low Frequency Oscillation》", 《 2021 IEEE 2ND CHINA INTERNATIONAL YOUTH CONFERENCE ON ELECTRICAL ENGINEERING》 *
TENG, YF等: "《Mechanism and Characteristics analysis of Ultra-low Frequency Oscillation phenomenon in a Power Grid with a High Proportion of Hydropower》", 《 INTERNATIONAL CONFERENCE ON POWER SYSTEM TECHNOLOGY》 *
吴熙;陈曦;吕万;袁超;杨宏宇;: "电力系统次同步振荡检测与在线定位技术综述", 电力自动化设备, no. 09 *
栾某德;刘涤尘;廖清芬;董超;欧阳利平;: "基于改进小波系数奇异值分解和小波去噪的低频振荡时变模式辨识", 电网技术, no. 06, 5 June 2012 (2012-06-05) *
薛安成;付潇宇;乔登科;王永杰;王嘉伟;: "风电参与的电力系统次同步振荡机理研究综述和展望", 电力自动化设备, no. 09 *

Also Published As

Publication number Publication date
CN116203350B (zh) 2024-02-27

Similar Documents

Publication Publication Date Title
Yang et al. Integrated mechanical and electrical DFIG wind turbine model development
Liu et al. Oscillation analysis and wide-area damping control of DFIGs for renewable energy power systems using line modal potential energy
Xie et al. Fundamental study of common mode small-signal frequency oscillations in power systems
CN110412462A (zh) 一种航空用永磁同步电动机瞬态温度预测方法
CN109871577A (zh) 一种基于叠加原理的永磁游标电机转矩分析方法
CN105762816A (zh) 一种抑制次同步相互作用的方法
Sarma et al. Modeling, analysis, and validation of controller signal interharmonic effects in DFIG drives
Sheng et al. Blade aerodynamic asymmetry fault analysis and diagnosis of wind turbines with doubly fed induction generator
Bashir et al. Frequency limited & weighted model reduction algorithm with error bound: Application to discrete-time doubly fed induction generator based wind turbines for power system
Ferrari et al. Open-Code, Real-Time Emulation Testbed Of Grid-Connected Type-3 Wind Turbine System With Hardware Validation
CN116203350B (zh) 一种高水电占比电网中频率振荡源定位方法
CN106294959A (zh) 模型参考自适应控制与双馈风机降阶模型相结合的建模仿真方法
Li et al. Improved identification method of doubly-fed induction generator based on trajectory sensitivity analysis
Schumann et al. Axial-flux-machine with an ironless fractional slot winding
CN113482852B (zh) 永磁直驱风力发电变换器控制方法、终端及存储介质
CN114254540A (zh) 基于数字栾生的永磁同步发电机的模拟方法以及模拟系统
Vargas et al. Experimental validation of a hybrid TD/FEHD model of a wind turbine generator for harmonic transient analysis
Almozayen et al. Dynamic phasor finite element modeling of grid-connected DFIG considering winding space harmonics
Dong et al. Comprehensive wave-to-wire model and control strategy design for wave energy conversion system
CN105589985A (zh) 风力发电机组参数对并网特性影响的确定方法、装置
Ding et al. Comparison investigation of fluid rheological characteristics effect on stator temperature field for large hydro-generator based on experiment
CN111740449A (zh) 基于附加励磁信号注入法的风电机群对同步发电机组阻尼特性分析方法
Singh et al. Time Domain Simulation of DFIG-Based Wind Power System using Differential Transform Method
Jovanović et al. Comparisons of vector control algorithms for doubly-fed reluctance wind generators
Yonggang et al. Blade mass imbalance fault diagnosis using rotor and stator current based on coordinate transformation and HED

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