CN111679124B - 一种电力***的振荡辨识的方法和装置 - Google Patents
一种电力***的振荡辨识的方法和装置 Download PDFInfo
- Publication number
- CN111679124B CN111679124B CN202010499195.4A CN202010499195A CN111679124B CN 111679124 B CN111679124 B CN 111679124B CN 202010499195 A CN202010499195 A CN 202010499195A CN 111679124 B CN111679124 B CN 111679124B
- Authority
- CN
- China
- Prior art keywords
- component
- subsynchronous
- supersynchronous
- frequency
- oscillation
- 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
- 230000010355 oscillation Effects 0.000 title claims abstract description 118
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000001228 spectrum Methods 0.000 claims abstract description 82
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 80
- 230000001360 synchronised effect Effects 0.000 claims abstract description 64
- 238000013016 damping Methods 0.000 claims abstract description 58
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000004458 analytical method Methods 0.000 claims description 18
- 230000003595 spectral effect Effects 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 3
- 241001505100 Succisa pratensis Species 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009615 fourier-transform spectroscopy Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
- G01R23/06—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage by converting frequency into an amplitude of current or voltage
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Complex Calculations (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明实施例提供了一种电力***的振荡辨识的方法和装置。所述方法,包括:在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗;对所述同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果;针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组;求解所述次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比;基于所述多对次同步和超同步振荡频谱信息,通过多次计算,获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比。
Description
技术领域
本发明涉及电力***领域,尤其涉及一种电力***的振荡辨识的方法和装置。
背景技术
电力***在正常运行过程中受到扰动或运行状态改变,可能产生各种振荡,振荡往往持续十几秒甚至几分钟。电力***的运行状态与其产生的振荡特性紧密相关。跨区域互联电网的长输电线在重载条件下易发生功率振荡,大量双馈风力发电机为主的大规模风力发电与输电通道串联补偿器同时工作的条件下易发生仅包含次同步分量的次同步振荡,大量永磁风力发电机为主的大规模风力发电在通过弱输电通道连接电网的条件下易发生次同步分量与超同步分量耦合的次同步振荡。
电力***广域测量***为电力***的动态监测提供了有效手段。广域测量***的测量数据是同步相量数据,同步相量数据的采样率通常为二倍电力***额定频率,可实现电力***动态过程的有效监测。但是,由于同步相量数据仅包含对应基波频率的基波相量,在电力***中出现振荡时,振荡分量仅由于频谱泄露而出现在基波相量中,因此同步相量数据中虽然包含振荡信息,但不能根据同步相量数据的频谱分析结果直接计算出电力***的振荡分量。
目前已有的基于同步相量的电力***振荡辨识方法,主要存在的技术问题表现在:一方面,针对同步相量数据进行频谱分析时,为保证足够高的频谱分辨率而不得不采用较大的数据窗,导致振荡辨识的实时性较差,例如为保证0.1Hz的频谱分辨率,则至少需要10秒长度的数据窗;进一步的,振荡辨识实时性差将导致辨识结果的变差,因为数据窗长度越大则电力***的振荡模式越可能变化。另一方面,已有的基于同步相量的电力***振荡辨识方法仅能辨识出频率和幅值,而不能辨识出相位,尤其是当振荡过程包含正的或负的阻尼比,即衰减因子时,更加难以辨识出衰减因子的影响而进一步严重影响辨识的准确性。
发明内容
本发明的实施例提供了一种电力***的振荡辨识的方法和装置,具有更好的实时性。
一种电力***的振荡辨识的方法,包括:
在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗;对所述同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果;
针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组;
求解所述次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比;
基于所述多对次同步和超同步振荡频谱信息,通过多次计算,获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比;
所述电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比为所述电力***的振荡进行辨识的结果。
一种电力***的振荡辨识的装置,包括:
第一分析单元,在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗;对所述同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果;
构建单元,针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组;
第一求解单元,求解所述次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比;
第二求解单元,基于所述多对次同步和超同步振荡频谱信息,通过多次计算,获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比;
辨识结果单元,所述电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比为所述电力***的振荡进行辨识的结果。
由上述本发明的实施例提供的技术方案可以看出,本发明实施例中,使用的傅里叶频谱分析数据窗长度小具有更好的实时性,频谱结果差分化克服了数据窗小导致的频谱混叠,可辨识出二倍额定频率以内的多对次同步分量和超同步分量的幅值、相位、频率和阻尼比。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的电力***振荡辨识方法的一种实施例的处理流程图;
图2为本发明提供的基于同步相量频谱差分分析的电力***振荡辨识方法的一种实施例的处理流程图;
图3为本发明提供的基于同步相量频谱差分分析的电力***振荡辨识方法的第二种实施例的同步相量数据序列示意图;
图4为本发明提供的基于同步相量频谱差分分析的电力***振荡辨识方法处理流程中的复数域傅里叶变换频谱分析结果示意图;
图5为本发明提供的基于同步相量频谱差分分析的电力***振荡辨识方法处理流程中的复数域傅里叶变换频谱分析结果第一对次同步与超同步分量的结果示意图;
图6为本发明提供的基于同步相量频谱差分分析的电力***振荡辨识方法处理流程中的复数域傅里叶变换频谱分析结果第二对次同步与超同步分量的结果示意图。
图7为本发明提供的电力***振荡辨识装置的一种实施例的连接示意图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
为便于对本发明实施例的理解,下面将结合附图以几个具体实施例为例做进一步的解释说明,且各个实施例并不构成对本发明实施例的限定。
如图1所示,为本发明所述的一种电力***的振荡辨识的方法,包括:
步骤11,在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗;对所述同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果;
步骤12,针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组;
步骤13,求解所述次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比;
步骤14,基于所述多对次同步和超同步振荡频谱信息,通过多次计算,获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比;
步骤15,将所述电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比,作为对所述电力***的振荡进行辨识的结果。
其中,所述在给定的电力***的同步相量数据序列中选取用于振荡辨识的同步相量数据窗,对所述同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果包括:
在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗k=1,…,K,其中K为复数域傅里叶变换频谱分析的数据窗的数据点数,同步相量数据的采样频率为fS,通过如下公式进行复数域傅里叶变换频谱分析,获得同步相量数据的复数域频谱分析结果m=1,…,M且M=K,
所述的复数域傅里叶变换频谱分析的数据窗的时间长度为K/fS。
其中,所述针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组的步骤包括:
构建次同步分量和超同步分量与频谱差分结果联立方程组,如下公式所示:
其中,“*”标记表示复数共轭,函数Q和D为如下公式所示
N为同步相量计算中傅里叶变换数据窗的数据点数,K为复数域傅里叶变换频谱分析的数据窗的数据点数,f0为电力***基波分量频率,fS为同步相量数据的采样频率且fS=2f0,M为同步相量数据的复数域频谱分析结果的个数,m1、m2、m3、m4为同步相量数据的复数域频谱分析结果的索引;
m1、m2、m3、m4的取值如下公式,
其中,所述的求解所述次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比包括:
根据以下公式所述的关系,获得次同步分量的频率fsub、幅值xsub、相位φsub和阻尼比δsub,超同步分量的频率fsup、幅值xsup、相位φsup和阻尼比δsup;
其中,电力***振荡中的瞬时值包括基波分量、多对次同步分量和超同步分量,瞬时值x(t)的表示为如下公式,
其中,f0、x0、φ0为基波分量的频率、幅值和相位,且f0等于电力***额定频率,共有I对次同步分量和超同步分量,xsub,i、φsub,i、fsub,i、δsub,i为第i个次同步分量的幅值、相位、频率和阻尼比,xsup,i、φsup,i、fsup,i、δsup,i为第i个超同步分量的幅值、相位、频率和阻尼比,且fsub,i+fsup,i=2f0。
相比较现有技术中同步相量频谱分析方法,本发明提供的方法使用的傅里叶频谱分析数据窗长度小具有更好的实时性,频谱结果差分化克服了数据窗小导致的频谱混叠,可辨识出二倍额定频率以内的多对次同步分量和超同步分量的幅值、相位、频率和阻尼比。
本发明提供了一种基于同步相量频谱差分分析的电力***振荡辨识方法,对同步相量数据窗进行复数域傅里叶变换频谱分析并获得频谱分析结果,构建次同步分量和超同步分量与频谱差分结果联立方程组;求解该方程组获得次同步分量和超同步分量的频率、幅值、相位和阻尼比,通过多次计算获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比。
以下描述本发明的应用场景。所述方法包括:
步骤1,在给定的电力***的同步相量数据序列中选取用于振荡辨识的同步相量数据窗,对同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果;
步骤2,针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组;
步骤3,求解次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比;
步骤4,基于该多对次同步和超同步振荡频谱信息,通过多次计算获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比。
其中,所述的电力***振荡包括基波分量、次同步分量和超同步分量,电力***振荡中的瞬时值包括基波分量、多对次同步分量和超同步分量,瞬时值x(t)的表示为如下公式,
其中,f0、x0、φ0为基波分量的频率、幅值和相位,且f0等于电力***额定频率,共有I对次同步分量和超同步分量,xsub,i、φsub,i、fsub,i、δsub,i为第i个次同步分量的幅值、相位、频率和阻尼比,xsup,i、φsup,i、fsup,i、δsup,i为第i个超同步分量的幅值、相位、频率和阻尼比,且fsub,i+fsup,i=2f0。
进一步的,在一些优选实施例中,第一个步骤所述的在给定的电力***的同步相量数据序列中选取用于振荡辨识的同步相量数据窗,对同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果包括如下子步骤:
在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗k=1,…,K,其中K为复数域傅里叶变换频谱分析的数据窗的数据点数,同步相量数据的采样频率为fS,复数域傅里叶变换频谱分析的数据窗的时间长度为K/fS,为提升所述方法的实时性,所述的复数域傅里叶变换频谱分析的数据窗的时间长度通常可选为2秒,对应的同步相量数据的复数域频谱分析结果的频率分辨率为0.5Hz,该频率分辨率较低,将导致复数域频谱分析结果的混叠。
复数域傅里叶变换频谱分析的数据窗的时间长度选择需兼顾方法的实时性和频谱分辨率,当复数域傅里叶变换频谱分析的数据窗的时间长度选为1秒以下或5秒以上时,虽然该方法的实时性更好或频谱分辨率更高,但将牺牲更多的频谱分辨率或实时性。
复数域傅里叶变换频谱分析的公式为
和
进一步的,所述的针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组包括以下子步骤:
在公式(14)和(15)中,频谱泄露分量和对复数域频谱分析结果产生影响,为消除频谱泄露分量和的影响,对复数域频谱分析结果进行二阶差分化处理,进而由公式(14)和(15)构建次同步分量和超同步分量与频谱差分结果联立方程组,如下公式所示
其中,“*”标记表示复数共轭,N为同步相量计算中傅里叶变换数据窗的数据点数,K为复数域傅里叶变换频谱分析的数据窗的数据点数,f0为电力***基波分量频率,fS为同步相量数据的采样频率且fS=2f0,M为同步相量数据的复数域频谱分析结果的个数。
m1、m2、m3、m4为同步相量数据的复数域频谱分析结果的索引,m1、m2、m3、m4的取值如下公式,
进一步的,所述的求解次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比包括:
求解公式(16)所述的次同步分量和超同步分量与频谱差分结果联立方程组,求解结果ωsub、ωsup、和均为复数,进一步根据公式(19)所述的关系,获得次同步分量的频率fsub、幅值xsub、相位φsub和阻尼比δsub,超同步分量的频率fsup、幅值xsup、相位φsup和阻尼比δsup,
进一步的,所述的基于该多对次同步和超同步振荡频谱信息,通过多次计算获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比包括:
电力***振荡过程中可能包含不止一对次同步和超同步振荡分量,当电力***振荡中包含多对次同步和超同步分量时,多次重复上述过程可计算得到其他对次同步和超同步分量;
本发明还提供一个实施例,用于示例性地显示利用本发明提供的方法进行电力***振荡辨识的过程,步骤如图2所示。
设电力***振荡中的瞬时值x(t)中包含5个分量,分别为:
1、基波分量,频率50Hz,阻尼比0,幅值100,相位为0rad;
2、第一个次同步分量,频率10.2Hz,阻尼比0.1,幅值40,相位2rad;
3、第一个超同步分量,频率89.8Hz,阻尼比0.05,幅值60,相位0.5rad;
4、第二个次同步分量,频率30.6Hz,阻尼比-0.2,幅值20,相位-1rad;
5、第二个超同步分量,频率69.4Hz,阻尼比-0.1,幅值15,相位1rad。
该瞬时值对应的同步相量如图3所示。
执行所述的第一步,图3所示的同步相量数据序列共有200个数据点,对数据序列进行复数域离散傅里叶变换频谱分析,复数域离散傅里叶变换频谱分析的数据窗长度有K=M=200个数据点,同步相量数据序列的采样率为fS=100Hz,数据窗时间长度为K/fS=2秒,因此复数域离散傅里叶变换频谱分析的频率分辨率为0.5Hz。
复数域离散傅里叶变换频谱分析的结果如图4所示,图4所示结果可见,共有两对次同步与超同步分量存在,第一对次同步分量与超同步分量出现在对应m=80和m=120附近,次同步分量的频率约为10Hz,超同步分量的频率约为90Hz;第二对次同步分量与超同步分量出现在对应m=39和m=161附近,次同步分量的频率约为30.5Hz,超同步分量的频率约为69.5Hz。
执行所述的第二步,针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组。
N为同步相量计算中傅里叶变换数据窗的数据点数,N=256;K为复数域傅里叶变换频谱分析的数据窗的数据点数,K=200;f0为电力***基波分量频率,f0=50Hz;fS为同步相量数据的采样频率且fS=2f0,fS=100Hz;M为同步相量数据的复数域频谱分析结果的个数,M=200。
[m1,m2,m3,m4]=[mk-2,mk-1,mk,mk+1]=[78,79,80,81],进而,
即,
执行所述的第三步,求解次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比。
进而根据计算获得第一对次同步分量的频率fsub=10.20Hz、幅值xsub=40.00、相位φsub=2.000rad和阻尼比δsub=0.1000,超同步分量的频率fsup=89.80Hz、幅值xsup=60.00、相位φsup=0.5000rad和阻尼比δsup=0.0500。
上述结果与实际结果相符。
执行所述的第四步,基于该多对次同步和超同步振荡频谱信息,通过多次计算获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比。
第二对次同步分量与超同步分量出现在对应m=39和m=161附近,次同步分量的频率约为30.5Hz,超同步分量的频率约为69.5Hz。第二对次同步和超同步振荡频谱信息对应的频谱结果如图6所示,由于复数域频谱分析结果的模在mk=39处的取得最大值且则[m1,m2,m3,m4]=[37,38,39,40]。进而重复上述步骤,计算得到第二对次同步分量的频率fsub=30.60Hz、幅值xsub=20.00、相位φsub=-1.000rad和阻尼比δsub=-0.2000,超同步分量的频率fsup=69.40Hz、幅值xsup=15.00、相位φsup=1.000rad和阻尼比δsup=-0.1000。
如上,可见本发明专利方法的正确性和实用性。
如图7所示,为本发明所述的一种电力***的振荡辨识的装置,包括:
第一分析单元,在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗;对所述同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果;
构建单元,针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组;
第一求解单元,求解所述次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比;
第二求解单元,基于所述多对次同步和超同步振荡频谱信息,通过多次计算,获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比;
辨识结果单元,所述电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比为所述电力***的振荡进行辨识的结果。
本领域普通技术人员可以理解:附图只是一个实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置或***实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的装置及***实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (5)
1.一种电力***的振荡辨识的方法,其特征在于,包括:
在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗;对所述同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果;
针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组,包括:构建次同步分量和超同步分量与频谱差分结果联立方程组,如下公式所示:
其中,“*”标记表示复数共轭,函数Q和D为如下公式所示
N为同步相量计算中傅里叶变换数据窗的数据点数,K为复数域傅里叶变换频谱分析的数据窗的数据点数,f0为电力***基波分量频率,fS为同步相量数据的采样频率且fS=2f0,M为同步相量数据的复数域频谱分析结果的个数,m1、m2、m3、m4为同步相量数据的复数域频谱分析结果的索引;
m1、m2、m3、m4的取值如下公式,
式中,次同步分量的频率fsub、幅值xsub、相位φsub和阻尼比δsub,超同步分量的频率fsup、幅值xsup、相位φsup和阻尼比δsup;
求解所述次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比;
基于所述多对次同步和超同步振荡频谱信息,通过多次计算,获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比;
所述电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比为对所述电力***的振荡进行辨识的结果。
5.一种电力***的振荡辨识的装置,其特征在于,包括:
第一分析单元,在给定的电力***的同步相量数据序列中,选取用于振荡辨识的同步相量数据窗;对所述同步相量数据窗进行复数域傅里叶变换频谱分析,获得包含多对次同步和超同步振荡频谱信息的同步相量数据的复数域频谱分析结果;
构建单元,针对一对次同步和超同步振荡频谱信息,构建次同步分量和超同步分量与频谱差分结果联立方程组,包括:构建次同步分量和超同步分量与频谱差分结果联立方程组,如下公式所示:
其中,“*”标记表示复数共轭,函数Q和D为如下公式所示
N为同步相量计算中傅里叶变换数据窗的数据点数,K为复数域傅里叶变换频谱分析的数据窗的数据点数,f0为电力***基波分量频率,fS为同步相量数据的采样频率且fS=2f0,M为同步相量数据的复数域频谱分析结果的个数,m1、m2、m3、m4为同步相量数据的复数域频谱分析结果的索引;
m1、m2、m3、m4的取值如下公式,
式中,次同步分量的频率fsub、幅值xsub、相位φsub和阻尼比δsub,超同步分量的频率fsup、幅值xsup、相位φsup和阻尼比δsup;
第一求解单元,求解所述次同步分量和超同步分量与频谱差分结果联立的方程组,获得次同步分量和超同步分量的频率、幅值、相位和阻尼比;
第二求解单元,基于所述多对次同步和超同步振荡频谱信息,通过多次计算,获得电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比;
辨识结果单元,所述电力***振荡中的各次同步分量和超同步分量的频率、幅值、相位和阻尼比为所述电力***的振荡进行辨识的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010499195.4A CN111679124B (zh) | 2020-06-04 | 2020-06-04 | 一种电力***的振荡辨识的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010499195.4A CN111679124B (zh) | 2020-06-04 | 2020-06-04 | 一种电力***的振荡辨识的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111679124A CN111679124A (zh) | 2020-09-18 |
CN111679124B true CN111679124B (zh) | 2021-07-23 |
Family
ID=72453799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010499195.4A Active CN111679124B (zh) | 2020-06-04 | 2020-06-04 | 一种电力***的振荡辨识的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111679124B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101068079A (zh) * | 2007-06-18 | 2007-11-07 | 北京交通大学 | 一种基于导纳原理实现输电线路保护的方法 |
WO2012134458A1 (en) * | 2011-03-30 | 2012-10-04 | Amsc Windtec Gmbh | Dynamic braking and low voltage ride through |
CN105223418A (zh) * | 2015-09-22 | 2016-01-06 | 清华大学 | 次同步和超同步谐波相量的测量方法及测量装置 |
CN106841778A (zh) * | 2016-12-28 | 2017-06-13 | 国电南瑞科技股份有限公司 | 基于pmu实现的次同步和超同步谐波参数的处理方法 |
CN107086588A (zh) * | 2017-05-27 | 2017-08-22 | 清华大学 | 电力电子***的次同步/超同步振荡分析方法及*** |
CN107247182A (zh) * | 2017-06-23 | 2017-10-13 | 华北电力大学 | 一种基于量测相量数据的间谐波分量还原方法 |
CN108957129A (zh) * | 2018-07-27 | 2018-12-07 | 中国南方电网有限责任公司 | 一种新能源并网***次/超同步振荡自适应监测方法 |
CN109473978A (zh) * | 2018-11-26 | 2019-03-15 | 北京金风科创风电设备有限公司 | 可再生能源场站振荡监测方法、控制方法和装置 |
CN110569577A (zh) * | 2019-08-26 | 2019-12-13 | 清华大学 | 基于频域模式法的新能源***次同步振荡溯源方法及*** |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE502008000783D1 (de) * | 2007-02-14 | 2010-07-29 | Semikron Elektronik Gmbh | Umrichterschaltung für einen doppeltgespeisten Asynchrongenerator mit variabler Leistungsabgabe und Verfahren zu deren Betrieb |
CN106655123A (zh) * | 2016-12-29 | 2017-05-10 | 北京四方继保自动化股份有限公司 | 一种次同步/超同步振荡监测保护装置及振荡监测保护方法 |
CN108303609B (zh) * | 2018-01-19 | 2024-02-02 | 国家电网公司 | 一种次同步振荡在线监测及扰动源定位方法及*** |
-
2020
- 2020-06-04 CN CN202010499195.4A patent/CN111679124B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101068079A (zh) * | 2007-06-18 | 2007-11-07 | 北京交通大学 | 一种基于导纳原理实现输电线路保护的方法 |
WO2012134458A1 (en) * | 2011-03-30 | 2012-10-04 | Amsc Windtec Gmbh | Dynamic braking and low voltage ride through |
CN105223418A (zh) * | 2015-09-22 | 2016-01-06 | 清华大学 | 次同步和超同步谐波相量的测量方法及测量装置 |
CN106841778A (zh) * | 2016-12-28 | 2017-06-13 | 国电南瑞科技股份有限公司 | 基于pmu实现的次同步和超同步谐波参数的处理方法 |
CN107086588A (zh) * | 2017-05-27 | 2017-08-22 | 清华大学 | 电力电子***的次同步/超同步振荡分析方法及*** |
CN107247182A (zh) * | 2017-06-23 | 2017-10-13 | 华北电力大学 | 一种基于量测相量数据的间谐波分量还原方法 |
CN108957129A (zh) * | 2018-07-27 | 2018-12-07 | 中国南方电网有限责任公司 | 一种新能源并网***次/超同步振荡自适应监测方法 |
CN109473978A (zh) * | 2018-11-26 | 2019-03-15 | 北京金风科创风电设备有限公司 | 可再生能源场站振荡监测方法、控制方法和装置 |
CN110569577A (zh) * | 2019-08-26 | 2019-12-13 | 清华大学 | 基于频域模式法的新能源***次同步振荡溯源方法及*** |
Non-Patent Citations (3)
Title |
---|
High-accurate estimation method of sub-synchronous and super-synchronous harmonic parameters on PMU;Tan Y. et.al;《The Journal of Engineering》;20171231(第13期);第1275-1279页 * |
基于PMU的分布式次同步振荡在线辨识方法;陈剑 等;《电气技术》;20171231(第10期);第92-95页 * |
风电并网***次/超同步振荡的动态监测方法研究;李金 等;《现代电力》;20190630;第36卷(第3期);第19-26页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111679124A (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ayhan et al. | On the use of a lower sampling rate for broken rotor bar detection with DTFT and AR-based spectrum methods | |
US5488281A (en) | Method and apparatus for predicting winding failure using zero crossing times | |
Karimi-Ghartemani et al. | Processing of symmetrical components in time-domain | |
Rauhala et al. | Detection of subsynchronous torsional oscillation frequencies using phasor measurement | |
JP2016214064A (ja) | 発電機の制動トルクの非侵入的な推定のためのシステムおよび方法 | |
CN106645919B (zh) | 基于三相瞬时功率的电网全谱功率振荡相量同步测量方法 | |
CN109861296B (zh) | 一种孤岛虚拟同步发电机***的惯性和阻尼辨识方法 | |
CN107144734B (zh) | 一种适用于pmu的配电网高精度相量测量方法 | |
Khelfi et al. | Induction motor rotor fault diagnosis using three-phase current intersection signal | |
Taskin et al. | Determination of the spectral properties and harmonic levels for driving an induction motor by an inverter driver under the different load conditions | |
Khalique et al. | A novel detection method for diagnosis of rotor eccentricity in three-phase induction motor | |
CN114374210A (zh) | 一种采用改进Prony算法的电力***低频振荡监控方法 | |
CN111679124B (zh) | 一种电力***的振荡辨识的方法和装置 | |
Bossa et al. | A field test to determine PSS effectiveness at multigenerator power plants | |
Cablea et al. | Method for computing efficient electrical indicators for offshore wind turbine monitoring | |
Amin et al. | Experimental validation of high-performance HIL-based real-time PMU model for WAMS | |
CN111679125B (zh) | 一种电力***振荡辨识的方法和装置 | |
Gomez et al. | Motor current signature analysis apply for external mechanical fault and cage asymmetry in induction motors | |
CN111273103B (zh) | 基于同步相量复数域频谱分析的电力***振荡辨识方法 | |
Wu et al. | Estimation of power system inertia using traditional swing equation, polynomial approximation and RV methods | |
CN111257636B (zh) | 一种电力***次同步和超同步振荡频率辨识方法 | |
Nayeripour et al. | An advanced analytical calculation and modeling of the electrical and mechanical harmonics behavior of Doubly Fed Induction Generator in wind turbine | |
CN114548149A (zh) | 电力***次同步振荡与超同步振荡的辨识方法和装置 | |
Garcia-Calva et al. | Rotor speed estimation for half-broken bar detection in induction motors using Kalman filtering | |
DK3047291T3 (en) | PROCEDURE AND DEVICE FOR DETERMINING THE ELECTRICAL PROPERTIES IN A WINDMILL |
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 |