CN114548149A - 电力***次同步振荡与超同步振荡的辨识方法和装置 - Google Patents
电力***次同步振荡与超同步振荡的辨识方法和装置 Download PDFInfo
- Publication number
- CN114548149A CN114548149A CN202210027959.9A CN202210027959A CN114548149A CN 114548149 A CN114548149 A CN 114548149A CN 202210027959 A CN202210027959 A CN 202210027959A CN 114548149 A CN114548149 A CN 114548149A
- Authority
- CN
- China
- Prior art keywords
- frequency
- component
- sinusoidal component
- oscillation
- fundamental wave
- 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.)
- Pending
Links
- 230000010355 oscillation Effects 0.000 title claims abstract description 148
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000001228 spectrum Methods 0.000 claims abstract description 71
- 230000001360 synchronised effect Effects 0.000 claims abstract description 51
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 15
- 150000001875 compounds Chemical class 0.000 claims description 12
- 230000003595 spectral effect Effects 0.000 claims description 11
- 238000010183 spectrum analysis Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Human Resources & Organizations (AREA)
- General Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Discrete Mathematics (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Signal Processing (AREA)
- Algebra (AREA)
- Artificial Intelligence (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种电力***次同步振荡与超同步振荡的辨识方法及装置,方法包括:基于电力***的振荡信号建立分析模型,包括一个基波正弦分量、一对频率耦合的次同步正弦分量和超同步正弦分量;将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量,将次同步正弦分量与超同步正弦分量分别合成为正频率振荡分量和负频率振荡分量;构建同步相量频谱拟合方程组,基于傅里叶变换后的同步相量序列,对频谱误差进行加权;以加权后的频谱误差的范数最小为目标,求解同步相量频谱拟合方程组,得到正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量,确定基波正弦分量、次同步正弦分量和超同步正弦分量的频率、幅值和相位。
Description
技术领域
本发明涉及电力***广域监测和次同步振荡技术领域,尤其涉及一种电力***次同步振荡与超同步振荡的辨识方法和装置。
背景技术
振荡过程是电力***中最常见的扰动过程之一。由于严重的振荡将导致***失稳、线路过载等问题而严重威胁电网安全,因此当电力***出现振荡过程时,需要有效监测并辨识振荡的参数。由于电力***是运行在额定频率下的交流电力***,当发生振荡时,振荡信号最主要包含偏移额定频率的基波正弦分量、频率之和为二倍额定频率的次同步正弦分量和超同步正弦分量,除此以外还有少量的其他频率分量。对于电力***运行管理人员而言,基波正弦分量、次同步正弦分量和超同步正弦分量的频率、幅值和相位是最重要的振荡参数。
随着监测手段的逐渐丰富,电力***广域测量***成为了新一代电力***动态监测***,成为了电力***动态监测与控制的有效手段。广域测量***的测量终端是相量测量终端PMU,相量测量终端将以最高二倍额定频率上传实时测量的同步相量数据。但是同步相量数据仅为计算出的对应额定频率的基波相量,在电力***振荡过程中的频率偏移基波分量、次同步分量和超同步分量将受频谱泄露影响而出现在同步相量数据的基波相量中,因此,同步相量数据中将包含振荡信息而可用于振荡辨识,需要基于同步相量的电力***振荡辨识方法以有效辨识出基波正弦分量、次同步正弦分量和超同步正弦分量的频率、幅值和相位。
目前已有的基于同步相量的电力***振荡辨识方法均为基于傅里叶变换频谱分析的方法,主要存在的技术问题表现在,针对同步相量数据进行频谱分析时,为保证足够高的频谱分辨率而不得不采用较大的数据窗以避免频谱混叠,导致振荡辨识的实时性较差;进而,振荡辨识实时性差将导致辨识结果的变差,因为数据窗长度越大则电力***的振荡模式越可能变化,因此有学者提出采用非线性超定方程组计算振荡参数,该方法能克服频谱分析方法中数据窗小而频谱混叠导致的频率分辨率差的缺点,可有效辨识出二倍额定频率以内的频率偏移基波正弦分量、次同步正弦分量和超同步正弦分量的幅值、相位和频率,但是该方法建立的是时域的同步相量轨迹拟合方程组,在以时域的误差最小为目标进行非线性曲线拟合求解时,难以控制振荡参数的求解精度,导致振荡参数的精度较低。
发明内容
本发明的实施例提供了一种电力***次同步振荡与超同步振荡的辨识方法和装置,用于解决现有技术中存在的问题。
为了实现上述目的,本发明采取了如下技术方案。
一种电力***次同步振荡与超同步振荡的辨识方法,包括:
S1基于获得的电力***的振荡信号,建立分析模型,该分析模型包括一个基波正弦分量、一对频率耦合的次同步正弦分量和超同步正弦分量;
S2将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量,将次同步正弦分量分解为正频率次同步正弦分量和负频率次同步正弦分量,将超同步正弦分量分解为正频率超同步正弦分量和负频率超同步正弦分量,将正频率次同步正弦分量和负频率超同步正弦分量合成为正频率振荡分量,将负频率次同步正弦分量和正频率超同步正弦分量合成为负频率振荡分量;
S3构建同步相量频谱拟合方程组,该同步相量频谱拟合方程组等号的左侧包括正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量,以及频谱误差,该同步相量频谱拟合方程组等号的右侧为傅里叶变换后的同步相量序列;
S4基于傅里叶变换后的同步相量序列,对频谱误差进行加权;
S5以加权后的频谱误差的最小范数为目标,通过加权最小二乘法求解同步相量频谱拟合方程组,获得正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量的计算值;
S6基于正频率基波正弦分量或负频率基波正弦分量的计算值,获得基波正弦分量的频率、幅值和相位;基于正频率振荡分量和负频率振荡分量的计算值,获得次同步正弦分量的频率、幅值和相位,以及超同步正弦分量的频率、幅值和相位。
优选地,步骤S1包括:
基于获得的电力***的振荡信号,通过式
x(t)=x0cos(2πf0t+φ0)+xsubcos(2πfsubt+φsub)+xsupcos(2πfsupt+φsup) (1)
建立分析模型;式中,x(t)为振荡信号瞬时值,f0、x0、φ0分别为瞬时值的基波正弦分量的频率、幅值和相位,xsub、φsub、fsub分别为次同步正弦分量的幅值、相位和频率,xsup、φsup、fsup分别为超同步正弦分量的幅值、相位和频率。
优选地,步骤S2包括:
通过式
将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量;式中,和分别为第k个同步相量的正频率基波正弦分量和负频率基波正弦分量,N为同步相量计算数据窗中的数据点数,fN为电力***额定频率,fs为同步相量数据上传频率,*为复数共轭标记;
通过式
通过式
优选地,步骤S3包括:
通过式
优选地,步骤S4包括:
通过式
优选地,步骤S6包括:
通过式
计算获得基波正弦分量的频率、幅值和相位;
通过式
fsup=2fN-fsub (11)
计算获得次同步正弦分量的频率、幅值和相位,以及超同步正弦分量的频率、幅值和相位。
第二方面,本发明提供一种电力***次同步振荡与超同步振荡的辨识装置,包括信号获取模块、分析计算模块和辨识模块;
信号获取模块用于:基于获得的电力***的振荡信号,建立分析模型,该分析模型包括一个基波正弦分量、一对频率耦合的次同步正弦分量和超同步正弦分量;
分析计算模块用于:将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量,将次同步正弦分量分解为正频率次同步正弦分量和负频率次同步正弦分量,将超同步正弦分量分解为正频率超同步正弦分量和负频率超同步正弦分量,将正频率次同步正弦分量和负频率超同步正弦分量合成为正频率振荡分量,将负频率次同步正弦分量和正频率超同步正弦分量合成为负频率振荡分量;
构建同步相量频谱拟合方程组,该同步相量频谱拟合方程组等号的左侧包括正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量,以及频谱误差,该同步相量频谱拟合方程组等号的右侧为傅里叶变换后的同步相量序列;
基于傅里叶变换后的同步相量序列,对频谱误差进行加权;
以加权后的频谱误差的范数最小为目标,采用加权最小二乘法求解同步相量频谱拟合方程组,得到正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量的计算值;
辨识模块用于:基于正频率基波正弦分量或负频率基波正弦分量的计算值,获得基波正弦分量的频率、幅值和相位;基于正频率振荡分量和负频率振荡分量的计算值,获得次同步正弦分量的频率、幅值和相位,以及超同步正弦分量的频率、幅值和相位。
第三方面,本发明提供一种计算机设备,包括存储器和处理器,该存储器存储有程序,该程序能够执行上述的方法,该处理器用于执行该程序。
第四方面,本发明提供一种计算机可读存储介质,该计算机可读存储介质存储有程序,该程序能够执行上述的方法。
由上述本发明的实施例提供的技术方案可以看出,本发明实施例提供的电力***次同步振荡与超同步振荡的辨识方法及装置具有以下优点:
一方面采用非线性超定方程组求解方法解决了频谱泄露不同分量互相影响的问题,能够有效辨识出次同步振荡参数和超同步振荡参数,而传统频谱分析方法是选取频谱的尖峰点作为频率点,认为频率尖峰点不受其他频谱的影响,或可通过差分方法减小这种影响;另一方面在求解非线性超定方程组时,将频谱幅值作为频谱误差的权重,增加了尖峰的影响,同时削弱了非尖峰的影响,即噪声的影响,进一步提升了在频谱混叠和噪声影响下的辨识结果精度。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的电力***次同步振荡与超同步振荡的辨识方法的处理流程图;
图2为本发明提供的电力***次同步振荡与超同步振荡的辨识装置的逻辑框图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元件被“连接”或“耦接”到另一元件时,它可以直接连接或耦接到其他元件,或者也可以存在中间元件。此外,这里使用的“连接”或“耦接”可以包括无线连接或耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的任一单元和全部组合。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
为便于对本发明实施例的理解,下面将结合附图以几个具体实施例为例做进一步的解释说明,且各个实施例并不构成对本发明实施例的限定。
本发明提供一种电力***次同步振荡与超同步振荡的辨识方法,用以提高次同步振荡与超同步振荡的辨识精度,图1为本发明实施例中电力***次同步振荡与超同步振荡的辨识方法的处理流程图,如图1所示,该方法包括:
S1:基于获得的电力***的振荡信号,建立分析模型,该分析模型包括一个频率可能偏移的基波正弦分量、一对频率耦合的次同步正弦分量和超同步正弦分量;
S2:将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量,将次同步正弦分量分解为正频率次同步正弦分量和负频率次同步正弦分量,将超同步正弦分量分解为正频率超同步正弦分量和负频率超同步正弦分量,将正频率次同步正弦分量和负频率超同步正弦分量合成为正频率振荡分量,将负频率次同步正弦分量和正频率超同步正弦分量合成为负频率振荡分量;
S3:构建同步相量频谱拟合方程组,该同步相量频谱拟合方程组等号的左侧包括正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量,以及频谱误差,该同步相量频谱拟合方程组等号的右侧为傅里叶变换后的同步相量序列;
S4:基于傅里叶变换后的同步相量序列,对频谱误差进行加权;
S5:以加权后的频谱误差的范数最小为目标,采用加权最小二乘法求解同步相量频谱拟合方程组,得到正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量的计算值;
S6:基于正频率基波正弦分量或负频率基波正弦分量的计算值,获得基波正弦分量的频率、幅值和相位;基于正频率振荡分量和负频率振荡分量的计算值,获得次同步正弦分量的频率、幅值和相位,以及超同步正弦分量的频率、幅值和相位。
在本发明提供的优选实施例中,步骤S1具体包括:
基于获得的电力***的振荡信号,通过式
x(t)=x0cos(2πf0t+φ0)+xsubcos(2πfsubt+φsub)+xsupcos(2πfsupt+φsup) (1)
建立分析模型;式中,x(t)为振荡信号瞬时值,f0、x0、φ0分别为瞬时值的基波正弦分量的频率、幅值和相位,xsub、φsub、fsub分别为次同步正弦分量的幅值、相位和频率,xsup、φsup、fsup分别为超同步正弦分量的幅值、相位和频率。
步骤S2具体包括:
通过式
将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量;式中,和分别为第k个同步相量的正频率基波正弦分量和负频率基波正弦分量,N为同步相量计算数据窗中的数据点数,fN为电力***额定频率,fs为同步相量数据上传频率,“*”为复数共轭标记,e是常数;
通过式
通过式
进一步的,在一个实施例中,具体实施步骤S2的过程如下。
快速傅里叶变换同步相量计算数据窗长度为一个电力***额定频率对应的一个周波,N为快速傅里叶变换同步相量计算数据窗中的数据点数,fS为同步相量数据上传频率且fS=2fN,由于fS=2fN且快速傅里叶变换同步相量计算数据窗为一个额定频率周波,故存在(ejkπ)一项。
由式(14)可见,由两部分组成,考虑到0<fsub<fN,第一部分对应的频率为为负频率;第二部分对应的频率为为正频率。相似地,由式(15)可见,由两部分组成,第一部分对应的频率为正频率且与的第二部分一致;第二部分为对应的频率为负频率且与的第一部分一致。
由式(16)可见,由两部分组成,两部分对应的频率分别为和考虑到可能存在f0<fN、f0=fN和f0>fN的情况,两部分对应频率分别为正频率和负频率或负频率和正频率。由于电力***的频率不会严格存在f0=fN的关系,因此不失一般性的先假设f0<fN并进行后续计算,并根据后续计算结果判断最终结果是f0<fN或f0>fN。
根据以上推导,将正频率次同步正弦分量和负频率超同步正弦分量合成为正频率振荡分量将负频率次同步正弦分量和正频率超同步正弦分量合成为负频率振荡分量如公式(4)和公式(5)所示,正频率基波正弦分量和负频率基波正弦分量分别为和如公式(2)和公式(3)所示。
在本发明提供的优选实施例中,步骤S3具体包括:
通过式
考虑到公式(1)的建模误差而得到:
现有技术根据公式(2)至(6)改写公式(15)构建同步相量轨迹拟合方程组,但是难以确定同步相量轨迹拟合方程组中的误差与振荡参数之间的关联关系,在以误差的范数最小为目标进行非线性曲线拟合求解时,难以控制振荡参数的求解精度,导致计算结果精度较低,而本发明实施例步骤S3中建立的是如公式(6)所示的同步相量频谱拟合方程组。
在本发明提供的优选实施例中,步骤S4包括:
通过式
具体实施时,在建立了同步相量频谱拟合方程组之后,步骤S4中,可以按照公式(7)对频谱误差进行加权处理,Fr(m)为同步相量的傅里叶频谱分析结果,能够反映振荡参数的频谱幅值特性,因此,加权后的频谱误差中考虑了振荡参数的频谱幅值特性,从而建立了频谱误差与振荡参数之间的关联关系。
可以采用加权最小二乘法,以加权后的频谱误差的范数最小为目标求解非线性超定方程组,利用解析或数值方法求得迭代雅克比矩阵,进而求得变量α和β,以及正频率基波正弦分量负频率基波正弦分量正频率振荡分量和负频率振荡分量
在本发明提供的优选实施例中,步骤S6包括:
通过式
计算获得基波正弦分量的频率、幅值和相位;
通过式
fsup=2fN-fsub (11)
计算获得次同步正弦分量的频率、幅值和相位,以及超同步正弦分量的频率、幅值和相位。
在一个具体的实施例中,根据求得的α和β及式可计算出公式(10)、(11)所示的次同步正弦分量的频率超同步正弦分量的频率fsup=2fN-fsub,以及假设f0<fN条件下的基波正弦分量的频率当f0<fN时,当f0>fN时,因此,根据求得的正频率基波正弦分量和负频率基波正弦分量可以按照公式(8)最终判断基波正弦分量的频率f0。
可以根据公式(8)计算得到的基波正弦分量的频率f0以及公式(4)计算Q*(f0,+1),进而根据公式(12)计算得到基波正弦分量的幅值x0和相位φ0。
相似地,根据次同步正弦分量的频率fsub和超同步正弦分量的频率fsup,以及公式(4)计算Q(fsub,-1)、Q(fsup,-1)、Q*(fsub,+1)和Q*(fsup,+1),进而根据公式(5)和公式(6)联立方程组,由于设定k为参考基准点,故可令k=0代入公式(5)和公式(6),得到公式(12),求解公式(12),得到次同步正弦分量的幅值和相位为xsub和φsub,超同步正弦分量的幅值和相位xsup和φsup。这样就求得了公式(1)中待定的基波正弦分量、次同步正弦分量和超同步正弦分量的幅值、相位和频率。
第二方面,本发明提供一种电力***次同步振荡与超同步振荡的辨识装置,图2为本发明提供的电力***次同步振荡与超同步振荡的辨识装置的逻辑框图,如图2所示,该装置包括信号获取模块1021、分析计算模块1022和辨识模块1023;
信号获取模块1021用于:基于获得的电力***的振荡信号,建立分析模型,该分析模型包括一个基波正弦分量、一对频率耦合的次同步正弦分量和超同步正弦分量;
分析计算模块1022用于:将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量,将次同步正弦分量分解为正频率次同步正弦分量和负频率次同步正弦分量,将超同步正弦分量分解为正频率超同步正弦分量和负频率超同步正弦分量,将正频率次同步正弦分量和负频率超同步正弦分量合成为正频率振荡分量,将负频率次同步正弦分量和正频率超同步正弦分量合成为负频率振荡分量;
构建同步相量频谱拟合方程组,该同步相量频谱拟合方程组等号的左侧包括正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量,以及频谱误差,该同步相量频谱拟合方程组等号的右侧为傅里叶变换后的同步相量序列;
基于傅里叶变换后的同步相量序列,对频谱误差进行加权;
以加权后的频谱误差的范数最小为目标,采用加权最小二乘法求解同步相量频谱拟合方程组,得到正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量的计算值;
辨识模块1023用于:基于正频率基波正弦分量或负频率基波正弦分量的计算值,获得基波正弦分量的频率、幅值和相位;基于正频率振荡分量和负频率振荡分量的计算值,获得次同步正弦分量的频率、幅值和相位,以及超同步正弦分量的频率、幅值和相位。
第三方面,本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述电力***次同步振荡与超同步振荡的辨识方法。
第四方面,本发明实施例还提供一种计算机可读存储介质,计算机可读存储介质存储有执行上述电力***次同步振荡与超同步振荡的辨识方法的计算机程序。
综上所述,本发明提供的一种电力***次同步振荡与超同步振荡的辨识方法及装置,如下优点:一方面采用非线性超定方程组求解方法解决了频谱泄露不同分量互相影响的问题,能够有效辨识出次同步振荡参数和超同步振荡参数,而传统频谱分析方法是选取频谱的尖峰点作为频率点,认为频率尖峰点不受其他频谱的影响,或可通过差分方法减小这种影响;另一方面在求解非线性超定方程组时,将频谱幅值作为频谱误差的权重,增加了尖峰的影响,同时削弱了非尖峰的影响,即噪声的影响,进一步提升了在频谱混叠和噪声影响下的辨识结果精度。
本领域普通技术人员可以理解:附图只是一个实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置或***实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的装置及***实施例仅仅是示意性的,其中作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (9)
1.电力***次同步振荡与超同步振荡的辨识方法,其特征在于,包括:
S1基于获得的电力***的振荡信号,建立分析模型,该分析模型包括一个基波正弦分量、一对频率耦合的次同步正弦分量和超同步正弦分量;
S2将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量,将次同步正弦分量分解为正频率次同步正弦分量和负频率次同步正弦分量,将超同步正弦分量分解为正频率超同步正弦分量和负频率超同步正弦分量,将正频率次同步正弦分量和负频率超同步正弦分量合成为正频率振荡分量,将负频率次同步正弦分量和正频率超同步正弦分量合成为负频率振荡分量;
S3构建同步相量频谱拟合方程组,该同步相量频谱拟合方程组等号的左侧包括正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量,以及频谱误差,该同步相量频谱拟合方程组等号的右侧为傅里叶变换后的同步相量序列;
S4基于傅里叶变换后的同步相量序列,对频谱误差进行加权;
S5以加权后的频谱误差的最小范数为目标,通过加权最小二乘法求解同步相量频谱拟合方程组,获得正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量的计算值;
S6基于正频率基波正弦分量或负频率基波正弦分量的计算值,获得基波正弦分量的频率、幅值和相位;基于正频率振荡分量和负频率振荡分量的计算值,获得次同步正弦分量的频率、幅值和相位,以及超同步正弦分量的频率、幅值和相位。
2.根据权利要求1所述的方法,其特征在于,步骤S1包括:
基于获得的电力***的振荡信号,通过式
x(t)=x0cos(2πf0t+φ0)+xsubcos(2πfsubt+φsub)+xsupcos(2πfsupt+φsup) (1)
建立分析模型;式中,x(t)为振荡信号瞬时值,f0、x0、φ0分别为瞬时值的基波正弦分量的频率、幅值和相位,xsub、φsub、fsub分别为次同步正弦分量的幅值、相位和频率,xsup、φsup、fsup分别为超同步正弦分量的幅值、相位和频率。
7.电力***次同步振荡与超同步振荡的辨识装置,其特征在于,包括信号获取模块、分析计算模块和辨识模块;
信号获取模块用于:基于获得的电力***的振荡信号,建立分析模型,该分析模型包括一个基波正弦分量、一对频率耦合的次同步正弦分量和超同步正弦分量;
分析计算模块用于:将基波正弦分量分解为正频率基波正弦分量和负频率基波正弦分量,将次同步正弦分量分解为正频率次同步正弦分量和负频率次同步正弦分量,将超同步正弦分量分解为正频率超同步正弦分量和负频率超同步正弦分量,将正频率次同步正弦分量和负频率超同步正弦分量合成为正频率振荡分量,将负频率次同步正弦分量和正频率超同步正弦分量合成为负频率振荡分量;
构建同步相量频谱拟合方程组,该同步相量频谱拟合方程组等号的左侧包括正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量,以及频谱误差,该同步相量频谱拟合方程组等号的右侧为傅里叶变换后的同步相量序列;
基于傅里叶变换后的同步相量序列,对频谱误差进行加权;
以加权后的频谱误差的范数最小为目标,采用加权最小二乘法求解同步相量频谱拟合方程组,得到正频率基波正弦分量、负频率基波正弦分量、正频率振荡分量和负频率振荡分量的计算值;
辨识模块用于:基于正频率基波正弦分量或负频率基波正弦分量的计算值,获得基波正弦分量的频率、幅值和相位;基于正频率振荡分量和负频率振荡分量的计算值,获得次同步正弦分量的频率、幅值和相位,以及超同步正弦分量的频率、幅值和相位。
8.一种计算机设备,包括存储器和处理器,其特征在于,该存储器存储有程序,该程序能够执行如权利要求1至6任一所述的方法,该处理器用于执行该程序。
9.一种计算机可读存储介质,其特征在于,该计算机可读存储介质存储有程序,该程序能够执行如权利要求1至6任一所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210027959.9A CN114548149A (zh) | 2022-01-11 | 2022-01-11 | 电力***次同步振荡与超同步振荡的辨识方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210027959.9A CN114548149A (zh) | 2022-01-11 | 2022-01-11 | 电力***次同步振荡与超同步振荡的辨识方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114548149A true CN114548149A (zh) | 2022-05-27 |
Family
ID=81668950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210027959.9A Pending CN114548149A (zh) | 2022-01-11 | 2022-01-11 | 电力***次同步振荡与超同步振荡的辨识方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114548149A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115792347A (zh) * | 2023-02-10 | 2023-03-14 | 安徽省国盛量子科技有限公司 | 基于微波移频的交流电测算方法 |
-
2022
- 2022-01-11 CN CN202210027959.9A patent/CN114548149A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115792347A (zh) * | 2023-02-10 | 2023-03-14 | 安徽省国盛量子科技有限公司 | 基于微波移频的交流电测算方法 |
CN115792347B (zh) * | 2023-02-10 | 2023-06-09 | 国网安徽省电力有限公司电力科学研究院 | 基于微波移频的交流电测算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPWO2008126240A1 (ja) | 同期フェーザ測定装置及びこれを用いた母線間位相角差測定装置 | |
CN1349102A (zh) | 检测交流电动机定子匝故障的方法和装置 | |
CN113285471B (zh) | 海上风电场次超同步振荡源感知与定位方法、装置及设备 | |
CN107851995A (zh) | 用于检测供给电网中的电压的方法和设备 | |
CN108155643B (zh) | 一种基于滑模观测器的单相电网电压参数的鲁棒估计方法 | |
CN111487561A (zh) | 一种分析电能计量二次接线错误及计算其正确电量的方法 | |
CN110333404B (zh) | 一种非侵入式的负荷监测方法、装置、设备及存储介质 | |
CN114548149A (zh) | 电力***次同步振荡与超同步振荡的辨识方法和装置 | |
Samal et al. | New signal subspace approach to estimate the inter‐area oscillatory modes in power system using TLS‐ESPRIT algorithm | |
CN114884094A (zh) | 一种风电场阻抗特性监测方法、***、设备和介质 | |
Znidi et al. | Coherency detection and network partitioning based on hierarchical DBSCAN | |
CN113156247B (zh) | 一种电力***低频振荡的预警方法及装置 | |
CN112557751B (zh) | 一种基于dft迭代法的谐波参数估计方法 | |
CN110334087A (zh) | 一种基于改进鲁棒容积卡尔曼滤波的数据清洗方法 | |
CN116595395B (zh) | 一种基于深度学习的逆变器输出电流预测方法及*** | |
CN111679125B (zh) | 一种电力***振荡辨识的方法和装置 | |
Tan et al. | High‐accurate estimation method of sub‐synchronous and super‐synchronous harmonic parameters on PMU | |
Shanmugapriya et al. | A fast Broyden's approximation‐based weighted least square state estimation for power systems | |
Chen et al. | Probabilistic load flow calculation method based on polynomial normal transformation and extended Latin Hypercube | |
Jia et al. | A novel adaptive linear prediction-based parameter estimation method for monitoring sub-/inter-harmonics during SSI events | |
CN111257636A (zh) | 一种电力***次同步和超同步振荡频率辨识方法 | |
Hassak et al. | Implementing a pmu based linear state estimator on a single board computer | |
Antonios et al. | Approximation of the time alignment error for measurements in electricity grids | |
CN111679124B (zh) | 一种电力***的振荡辨识的方法和装置 | |
CN118091526B (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 |