CN110425092B - 一种基于能量频谱的风电并网系统振荡源定位系统及方法 - Google Patents
一种基于能量频谱的风电并网系统振荡源定位系统及方法 Download PDFInfo
- Publication number
- CN110425092B CN110425092B CN201910694996.3A CN201910694996A CN110425092B CN 110425092 B CN110425092 B CN 110425092B CN 201910694996 A CN201910694996 A CN 201910694996A CN 110425092 B CN110425092 B CN 110425092B
- Authority
- CN
- China
- Prior art keywords
- unit
- energy
- energy spectrum
- oscillation source
- 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 222
- 238000001228 spectrum Methods 0.000 title claims abstract description 166
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000001360 synchronised effect Effects 0.000 claims abstract description 29
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000010183 spectrum analysis Methods 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 12
- 230000001174 ascending effect Effects 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000010248 power generation Methods 0.000 abstract description 2
- 230000000630 rising effect Effects 0.000 abstract description 2
- 238000013016 damping Methods 0.000 description 10
- 230000000737 periodic effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- 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
-
- 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/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
- H02J3/0012—Contingency detection
-
- 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
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
- H03L7/00—Automatic control of frequency or phase; Synchronisation
- H03L7/06—Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
- H03L7/08—Details of the phase-locked loop
- H03L7/099—Details of the phase-locked loop concerning mainly the controlled oscillator of the loop
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
- H03L7/00—Automatic control of frequency or phase; Synchronisation
- H03L7/06—Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
- H03L7/16—Indirect frequency synthesis, i.e. generating a desired one of a number of predetermined frequencies using a frequency- or phase-locked loop
- H03L7/18—Indirect frequency synthesis, i.e. generating a desired one of a number of predetermined frequencies using a frequency- or phase-locked loop using a frequency divider or counter in the loop
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明涉及一种基于能量频谱的风电并网系统振荡源定位系统及方法,属于风力发电系统技术领域,解决了现有针对系统受扰后激发多振荡模态的情况下难以准确定位振荡源机组的难题。该方法包括以下步骤:采集系统各机组输出端口的电压和电流信息,获得各发电机组的动态能量随时间变化曲线;将呈上升趋势的曲线对应的发电机组作为备选机组;获取并处理备选机组中同步机组、含锁相环双馈风电机组的能量频谱和/或含虚拟惯量的双馈风电机组的类能量频谱,将主导振荡模态最大的机组作为振荡源参考机组;计算其余每一备选机组与振荡源参考机组的能量频谱相似度系数,确定振荡源机组。本发明实现了系统受扰后激发多振荡模态的情况下振荡源机组的准确定位。
Description
技术领域
本发明涉及风力发电系统技术领域,尤其涉及一种基于能量频谱的风电并网系统振荡源定位系统及方法。
背景技术
虚拟惯量作为双馈风电机组的附加控制环节,虽可在一定程度上改善系统的惯量和频率响应特性,但也打开了双馈风电机组在扰动下与系统功率互通的通道,加强了机网间的耦合作用:虚拟惯量会降低风电并网系统的阻尼水平,诱使双馈风电机组参与的振荡模态呈现弱阻尼,系统受扰后易振荡发散失稳。另一方面,系统振荡过程中,除了主导振荡模态还可能激发其余弱阻尼振荡模态,多振荡模态相耦合使得系统振荡波形中的频率分量更为复杂,主导振荡源机组难以被准确定位。因此,为保障风电并网系统受扰后的稳定水平,对主导模态下的振荡源机组进行准确定位已成为亟待解决的关键问题。
目前,电力系统低频振荡扰动源定位已逐渐成为国内外众多专家和学者研究的热点。然而,现有振荡源定位的研究主要针对传统电力系统,均未考虑含虚拟惯量的风电并网系统受扰产生多振荡模态耦合的影响,难以区分主导振荡模态和非主导振荡模态的振荡源机组,无法准确定位系统振荡源。
发明内容
鉴于上述的分析,本发明旨在提供一种基于能量频谱的风电并网系统振荡源定位系统及方法,用以解决针对系统受扰后激发多振荡模态的情况下难以准确定位振荡源机组的难题。
本发明的目的主要是通过以下技术方案实现的:
一方面,提供一种基于能量频谱的风电并网系统振荡源定位方法,包括以下步骤:
步骤1:采集系统各机组输出端口的电压和电流信息,获得各发电机组的动态能量随时间变化曲线;
步骤2:将所述动态能量随时间变化曲线中呈上升趋势的曲线对应的发电机组作为备选机组;
步骤3:获取备选机组中同步机组、含锁相环双馈风电机组的能量频谱和/或含虚拟惯量的双馈风电机组的类能量频谱;
步骤4:对所述能量频谱和/或类能量频谱的幅值进行标幺化处理,将主导振荡模态最大的机组作为振荡源参考机组;
步骤5:在振荡源参考机组的主导振荡模态±0.3Hz频段下,计算其余每一备选机组与所述振荡源参考机组的能量频谱相似度系数,根据相似度系数大小确定振荡源机组。
进一步地,所述步骤1中,通过下式得到各发电机组的动态能量,以便获得各发电机组的动态能量随时间变化曲线:
其中:Pe为各机组输出有功功率;δ为xy轴与dq轴的夹角,当该发电机组为同步机组时,δ为机组的功角;当该机组为含锁相环双馈风电机组时,δ为锁相环的锁相角δpll;当该机组为含虚拟惯量的双馈风电机组时,δ为锁相环的锁相角δpll;Ud、Uq、Id和Iq分别为各机组输出端口的电压和电流的d、q轴分量。
进一步地,所述步骤3包括,将各发电机组的动态能量在时域上进行傅里叶分解,得到频率成分f(k)对应的能量幅值A(k),k=1,2,3,4,…,N/2,N为离散信号的长度;
同步机组、含锁相环双馈风电机组的能量频谱均为A(k);
含虚拟惯量的双馈风电机组的类能量频谱为A(k)/f2(k)。
进一步地,还包括对所述能量频谱的幅值和/或类能量频谱的幅值进行归一化处理后,再对归一化后的所述能量频谱和/或类能量频谱的幅值进行标幺化处理;所述归一化处理的计算公式为:
其中,S(k)为各频率对应的能量频谱幅值,同步机组、含锁相环双馈风电机组的能量频谱幅值S(k)=A(k),含虚拟惯量的双馈风电机组的类能量频谱幅值S(k)=A(k)/f2(k),为能量频谱幅值的平均值,N为频率个数,S'(k)为归一化处理后的能量频谱幅值。
进一步地,步骤4中,对0.1Hz-2.5Hz频段上所有归一化后的能量频谱幅值和/或类能量频谱的幅值进行标幺化处理,得到每台机组各振荡模态能量频谱幅值和/或类能量频谱的幅值对相应机组所有能量频谱幅值之和的占比,所述占比最大的机组为振荡源参考机组。
进一步地,步骤5中,通过下式计算其余每一备选机组与所述振荡源参考机组的能量频谱相似度系数:
其中,x,y分别为主导振荡模态±0.3Hz频段下两个机组的能量频谱数据,r为计算得出的相似度系数。
进一步地,步骤5中,当与所述振荡源参考机组的能量频谱相似度系数大于等于0.7±0.05时,将该机组确定为振荡源机组。
另一方面,提供一种基于能量频谱的风电并网系统振荡源定位系统,其特征在于,包括数据采集模块、动态能量计算模块、能量频谱分析模块和结果输出模块;
所述数据采集模块用于采集系统各机组输出端口的电压和电流信息,并将其传递给动态能量计算模块;
所述动态能量计算模块接收数据采集模块传递的数据,计算各发电机组的动态能量随时间变化趋势,并将其传递给能量频谱分析模块;
所述能量频谱分析模块接收各机组输出端口的动态能量数据,进行能量频谱计算,并通过对比各机组与参考振荡源机组能量频谱的相似度,定位系统振荡源;
所述结果输出模块用于输出主导振荡模态及对应的系统振荡源机组。
进一步地,所述动态能量计算模块利用下述公式计算各发电机组的动态能量随时间变化曲线:
其中:Pe为各机组输出有功功率;δ为xy轴与dq轴的夹角,当该发电机组为同步机组时,δ为机组的功角;当该机组为含锁相环双馈风电机组时,δ为锁相环的锁相角δpll;当该机组为含虚拟惯量的双馈风电机组时,δ为锁相环的锁相角δpll;Ud、Uq、Id和Iq分别为各机组输出端口的电压和电流的d、q轴分量。
进一步地,所述能量频谱分析模块的工作过程为:
将所述动态能量随时间变化曲线中呈上升趋势的曲线对应的发电机组作为备选机组;
获取备选机组中同步机组、含锁相环双馈风电机组的能量频谱和/或含虚拟惯量的双馈风电机组的类能量频谱;
对所述能量频谱和/或类能量频谱的幅值进行标幺化处理,将主导振荡模态最大的机组作为振荡源参考机组;
在振荡源参考机组的主导振荡模态±0.3Hz频段下,计算其余每一备选机组与所述振荡源参考机组的能量频谱相似度系数,根据相似度系数大小确定振荡源机组。
与现有技术相比,本发明至少具有如下有益效果之一:
1)本发明提供的基于能量频谱的风电并网系统振荡源定位方法,通过测量机组端口的电压、电流信息,实时监测振荡过程中各发电机组的动态能量随时间变化曲线,初步确定备选机组,实现系统阻尼的在线评估以及振荡分量传导路径的在线追踪;通过比较振荡模态临近频段下,不同机组端口的动态能量频谱相似度,能够避免非主导振荡模态对振荡源定位的影响,实现复杂系统工况下主导振荡源机组的快速、准确识别,有助于及时隔离故障机组,抑制振荡的发展和发散,保障系统的安全稳定运行。
2)本发明提供的基于能量频谱的风电并网系统振荡源定位系统,所述能量频谱分析模块接收各机组输出端口的动态能量数据,进行能量频谱计算,并通过对比各机组与参考振荡源机组能量频谱的相似度,能够快速、准确的识别并定位出系统振荡源,实现了在线识别振荡源,为系统稳定安全运行提供有力的保障。
本发明中,上述各技术方案之间还可以相互组合,以实现更多的优选组合方案。本发明的其他特征和优点将在随后的说明书中阐述,并且,部分优点可从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过说明书、权利要求书以及附图中所特别指出的内容中来实现和获得。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1为实施例中基于能量频谱的风电并网系统振荡源定位方法流程图;
图2为实施例中锁相环模型;
图3为实施例中典型虚拟惯量控制框图;
图4为实施例中含双馈风电机组的互联系统模型;
图5为实施例中新英格兰10机39节点系统;
图6为实施例中各发电机组动态能量变化曲线;
图7为实施例中备选机组能量频谱;
图8为实施例中各机组能量频谱相似度柱状图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例,其中,附图构成本申请一部分,并与本发明的实施例一起用于阐释本发明的原理,并非用于限定本发明的范围。
本发明涉及一种基于能量频谱的风电并网系统振荡源定位方法,该方法是基于振荡源机组的频谱特征是主导振荡模态远大于非主导振荡模态这一理论支撑实现的。为了便于理解本申请的技术方案,先通过以下理论分析,详细阐述振荡源机组的频谱特征是主导振荡模态远大于非主导振荡模态,具体分析过程如下:
步骤S1、建立含虚拟惯量及锁相环控制环节的双馈风电机组动态能量模型。
锁相环模型动态方程可以表示为:
ωs_PLL=KI_PLLxPLL-KP_PLLusd
式中:xPLL为跟踪定子d轴电压的误差积累,usd为双馈风机定子d轴电压,δPLL表示观测到的定子电压矢量领先xy坐标中x轴的角度,ωs_PLL为锁相环测得的dq标系旋转角速度,ωn为电网频率,Kp和Ki为锁相环PI控制器参数,如图2所示。
虚拟惯量控制模型可以表示为:
根据系统频率偏差,通过虚拟惯量控制释放出的有功功率变化量Pvir的数学表达式为:
其中,Kω为虚拟惯量控制的比例参数,Kω>0。
引入虚拟惯量控制后,转子变流器的有功参考值Pref随系统频率的变化进行相应调整。由虚拟惯量引起的双馈风电机组定子侧有功功率变化值ΔPs为:
基于节点电流方程的系统能量函数构建方法,对节点电压UB、发电机注入电流IG及负荷电流IL沿系统轨迹积分得到系统的能量:
W=∫Im(((YUB-IG+IL)*)TdUB) (1)
式中,Y为系统导纳矩阵,UB为母线电压列向量,IG和IL分别为发电机节点注入电流和负荷节点注入电流列向量。Im表示取复数的虚部。
取系统能量中的发电机部分能量,并将节点电压、节点电流值进行dq变换,结合发电机转子运动方程,得到机端动态能量表达式:
其中:Pe为各机组输出有功功率;δ为xy轴与dq轴的夹角,当该发电机组为同步机组时,δ为机组的功角;当该机组为含锁相环双馈风电机组时,δ为锁相环的锁相角δpll;当该机组为含虚拟惯量的双馈风电机组时,δ为锁相环的锁相角δpll;Ud、Uq、Id和Iq分别为各机组输出端口的电压和电流的dq轴分量。
考虑到直接测量的电气量中包含与振荡无关的稳态分量,由式(1)得到的动态能量无法准确反映振荡中能量的累积过程。为滤除能量稳态分量的影响,利用电气量的变化量计算振荡过程中产生的动态能量为:
式(3)中各机组风机端口动态能量分为两部分:
ΔW1=-∫ΔPedΔδpll (4)
其中:ΔW1部分为锁相环产生的能量,ΔW2部分为网侧线路产生的能量。
虚拟惯量将锁相环输出的频率值作为输入,通过改变定子输出有功功率实现系统惯量支撑。因此,虚拟惯量主要对式(4)中ΔW1产生影响,将含虚拟惯量的有功功率变化量代入式(4)中,可得含虚拟惯量的ΔW1表达式为:
式(5)中ΔW2主要由线路参数决定。设风电并网系统中出现ωd频率的低频振荡,电流振荡分量幅值为In,双馈风机并网点电流dq轴分量可分别表示为:
并网点电压和电流的dq轴分量之间的关系可表示为:
将式(7)和式(8)代入式(5)中可得:
由式(9)可知,线路产生的动态能量为恒定值,不包含周期分量,且不影响动态能量的上升和下降趋势。因此,低频振荡过程中,双馈风电机组的动态能量主要受锁相环和虚拟惯量影响。
步骤S2、根据式(6)得出的含虚拟惯量的双馈风电机组动态能量表达式,推导含虚拟惯量及锁相环的双馈风电机组作为振荡源和非振荡源时的能量频谱,并分析其能量频谱特征,具体过程如下:
步骤S201、构建含双馈风电机组的互联电力系统模型。
对互联系统而言,低频振荡存在两种振荡模式,分别为局部振荡模式和区域间振荡模式。这两种振荡模式均可用两个机群间的相对振荡来表示。因此,本发明以等值的三机系统为例,进行能量频谱的理论推导,系统结构图如图4所示。其中,G3为参考机群。系统受扰后,G1和G3间产生的振荡模态为ω1,且为主导振荡模态,G1为振荡源机群;G2和G3间产生的振荡模态为ω2,为非主导振荡模态,G2为非振荡源机群;母线4的电压相角产生的振荡分量为Δθ4=asinω1t+bsinω2t,且a>b。
步骤S202、推导同步机组、含锁相环双馈风电机组以及含虚拟惯量的双馈风电机组作为振荡源和非振荡源时的能量频谱。
(1)同步机作为振荡源时的能量频谱表达式为:
其中,C1、C2、M1、N1、M2和N2分别为G1功角振荡的幅值分量系数,
其中,KD、Ks、KN分别为:
(2)同步机作为非振荡源时的能量频谱表达式为:
其中,G0、H0、G1、H1、G2和H2分别为G2功角振荡的幅值分量系数。
其中,KD2、Ks2、KN2分别为:
(3)含锁相环双馈风电机组作为振荡源时的能量频谱表达式为:
其中,A0、B0、A1、B1、A2和B2的表达式分别为:
(4)含虚拟惯量双馈风电机组作为振荡源时的能量频谱表达式为:
步骤S203、分析同步机组、含锁相环双馈风电机组以及含虚拟惯量的双馈风电机组作为振荡源和非振荡源时的能量频谱特征。
含锁相环的双馈风电机组为振荡源时,其能量频谱与作为振荡源的同步机能量频谱类似,即主导频率分量远大于其他频率分量。当双馈风电机组为非振荡源时,其能量频谱与同步机非振荡源能量频谱类似,与振荡源机组能量频谱具有明显区别。
虚拟惯量的引入会改变风电机组的能量频谱,即在各能量分量幅值上引入频率系数无法直接判断主导振荡模态和非主导振荡模态的大小。因此,本发明构建类能量频谱,以消除虚拟惯量引入的频率系数的影响。双馈风电机组的类能量频谱可表示为:
含虚拟惯量的双馈风电机组类能量频谱与含锁相环的双馈风电机组能量频谱特征一致,因此,含锁相环和虚拟惯量的双馈风电机组振荡源与非振荡源之间的类能量频谱有明显区别。
根据上述能量频谱分析可知,振荡源机组的能量频谱特征为:主导振荡模态能量幅值分量远大于非主导振荡模态能量幅值分量,而非振荡源机组的能量频谱特征为:主导振荡模态能量幅值分量与非主导振荡模态能量幅值分量接近,甚至小于非主导振荡模态幅值分量。
实施例一
基于前述理论分析,本实施例中的基于能量频谱的风电并网系统振荡源定位方法,如图1所示,具体包括以下步骤:
步骤1:采集系统各机组输出端口的电压U和电流IG信息,代入式(2)中,以便获得各发电机组的动态能量随时间变化曲线。
步骤2:将动态能量随时间变化曲线中呈现上升趋势的曲线对应的发电机组划分为备选机组,即将负阻尼机组划分为备选机组。
步骤3:获取备选机组中同步机组、含锁相环双馈风电机组的能量频谱和/或含虚拟惯量的双馈风电机组的类能量频谱。
将各发电机组的动态能量在时域上进行傅里叶分解,得到频率成分f(k)对应能量幅值为A(k),k=1,2,3,4,…,N/2,N为离散信号的长度;同步机组、含锁相环双馈风电机组的能量频谱均为A(k);含虚拟惯量的双馈风电机组的类能量频谱为A(k)/f2(k)。
由于能量频谱幅值数据中存在奇异样本数据,会影响相似度分析结果的准确性,因此为了消除数据中噪声影响,利用下述公式对能量频谱的幅值和/或类能量频谱的幅值数据进行归一化处理:
其中,A(k)为各频率对应的能量频谱幅值,同步机组、含锁相环双馈风电机组的能量频谱幅值S(k)=A(k),含虚拟惯量的双馈风电机组的类能量频谱幅值为S(k)=A(k)/f2(k),为能量频谱幅值的平均值,N为频率个数,S'(k)为归一化处理后的能量频谱幅值。
步骤4:对能量频谱的幅值和/或类能量频谱的幅值进行归一化处理后,对能量频谱和/或类能量频谱的幅值进行标幺化处理,得到每台机组各振荡模态能量频谱幅值和/或类能量频谱的幅值对相应机组所有能量频谱幅值之和的占比,所述占比最大的机组为振荡源参考机组。
动态能量中包含非周期分量,可能影响各振荡模态的能量频谱幅值大小。由于每台机组都有多个振荡模态能量频谱幅值,为消除非周期分量的影响,便于统一比较不同振荡模态频谱幅值大小,对低频振荡的频段如0.1Hz-2.5Hz频段上所有归一化后能量频谱幅值和/或类能量频谱的幅值进行标幺化处理,即求取各振荡模态能量频谱幅值和/或类能量频谱的幅值对相应机组所有能量频谱幅值之和的占比,并从备选机组中选取主导振荡模态占比最大的机组,作为振荡源参考机组。
步骤5:在振荡源参考机组的主导振荡模态±0.3Hz频段下,计算其余每一备选机组与振荡源参考机组的能量频谱相似度系数,根据相似度系数大小确定振荡源机组。
其中,x,y分别为主导振荡模态±0.3Hz频段下两个机组的能量频谱数据,r为计算得出的相似度系数。
通过比较相似度系数r的大小,当与振荡源参考机组的能量频谱相似度系数大于等于0.7±0.05时,该机组为振荡源机组,当与振荡源参考机组的能量频谱相似度系数小于0.7±0.05时,该机组为非振荡源机组。
与现有技术相比,本实施例提供的基于能量频谱的风电并网系统振荡源定位方法,根据同步机能量频谱与双馈风电机组的类能量频谱特征,提出了一种基于能量频谱相似度的低频振荡源追溯方法,通过测量机组端口的电压、电流信息,实时监测振荡过程中各发电机组的动态能量随时间变化曲线,初步确定备选机组,实现系统阻尼的在线评估以及振荡分量传导路径的在线追踪;通过比较振荡模态临近频段下,不同机组端口的动态能量频谱相似度,能够避免非主导振荡模态对振荡源定位的影响,实现复杂系统工况下主导振荡源机组的快速、准确识别,有助于及时隔离故障机组,抑制振荡的发展和发散,保障系统的安全稳定运行。与传统方法相比,本发明在系统受扰后激发多振荡模态的情况下,能够避免非主导振荡模态对振荡源定位的影响,实现复杂系统工况下振荡源机组的准确识别。
实施例二
本实施例提供一种基于能量频谱的风电并网系统振荡源定位系统,包括数据采集模块、动态能量计算模块、能量频谱分析模块和结果输出模块;
所述数据采集模块用于采集系统各机组输出端口的电压和电流信息,并将其传递给动态能量计算模块;
所述动态能量计算模块接收数据采集模块传递的数据,计算各发电机组的动态能量随时间变化趋势,并将其传递给能量频谱分析模块;
所述能量频谱分析模块接收各机组输出端口的动态能量数据,进行能量频谱计算,并通过对比各机组与参考振荡源机组能量频谱的相似度,定位系统振荡源;
所述结果输出模块用于输出主导振荡模态及对应的系统振荡源机组。
本实施例中,所述动态能量计算模块利用下述公式计算各发电机组的动态能量随时间变化曲线:
其中:Pe为各机组输出有功功率;δ为xy轴与dq轴的夹角,当该发电机组为同步机组时,δ为机组的功角;当该机组为含锁相环双馈风电机组时,δ为锁相环的锁相角δpll;当该机组为含虚拟惯量的双馈风电机组时,δ为锁相角;Ud、Uq、Id和Iq分别为各机组输出端口的电压和电流的d、q轴分量。
本实施例中,所述能量频谱分析模块的工作过程为:
将所述动态能量随时间变化曲线中呈上升趋势的曲线对应的发电机组作为备选机组;
获取备选机组中同步机组、含锁相环双馈风电机组的能量频谱和/或含虚拟惯量的双馈风电机组的类能量频谱;
对所述能量频谱和/或类能量频谱的幅值进行标幺化处理,将主导振荡模态最大的机组作为振荡源参考机组;
在振荡源参考机组的主导振荡模态±0.3Hz频段下,计算其余每一备选机组与所述振荡源参考机组的能量频谱相似度系数,根据相似度系数大小确定振荡源机组。
本实施例供的基于能量频谱的风电并网系统振荡源定位系统,能量频谱分析模块接收各机组输出端口的动态能量数据,进行能量频谱计算,并通过对比各机组与参考振荡源机组能量频谱的相似度,能够快速、准确的识别并定位出系统振荡源,实现了在线识别振荡源,为系统稳定安全运行提供有力的保障。由于本实施例提供的系统与实施例一提供的方法实施例原理相同,所以本系统也具有上述方法实施例一相应的技术效果。
实施例三
本实施例采用双馈风机并联等效模型来模拟整个风电场,进一步验证实施例一的定位方法,如图3所示的新英格兰10机39节点系统图,风电场由1000台完全相同的双馈异步风机并联组成,风机出口并联在母线39处向电网供电,每台风机的额定容量为1.5MW,并且按照额定出力运行;风电场首先经过0.69/20kV场内变压器,再经过20/230kV变压器连接到母线5处进行输电,其余发电机为同步发电机。
本实施例设置场景为:t=2s时在风机机械转矩设置一个频率为0.48Hz的周期性小扰动引发系统振荡,风机的有功功率发生近似幅值为43kW的等幅振荡。
运行动态能量计算模块,结果为:
测量系统中各发电机组端口有功功率、相角及电压电流值,计算各发电机组振荡过程中动态能量变化曲线,如图6所示。其中,DFIG、G2和G4的动态能量均呈现上升趋势,即为负阻尼作用。
运行能量频谱分析模块,结果为:
将三台负阻尼机组划分为备选机组,进行FFT变换,求取动态能量频谱并进行标准化如图7所示。其中DFIG采取类能量频谱进行分析。由图7可知,四台机组的最大能量频谱均出现在2倍振荡频率0.96Hz处。同时,对比该频段下所有机组的能量频谱可知,DFIG主导振荡模态占比最大,根据振荡源能量频谱特征判断此时DFIG为振荡源机组之一。
为了进一步验证其他机组是否为振荡源,本实施例取主导振荡模态±0.3Hz的频段,根据相似度公式,对各机组能量频谱相似度进行计算,由于相似度系数矩阵为实对称阵,故取上三角阵数据,同时为了对比效果更加显著,将各机组与其自身相似度(应为1)置0作柱状图,如图8所示。根据图8可知,在该频段下,其余各机组与DFIG的能量相似度均较在0.6以下,相关性较低,因此主导振荡模态下振荡源仅为DFIG,G2、G3和G4均不是主导振荡模态下的振荡源。
为探究G2和G4被主导振荡模态激发的新振荡模态,对比G2和G4的能量频谱可知,除了主导振荡频率分量,这两台机组中还存在1.64Hz左右的非主导能量频率分量。本发明进一步针对1.64±0.3Hz的频段进行能量频谱分析。对比G2、G3和G4中1.64Hz频率分量的占比可知,G2中该频率分量占比最大,因此G2为1.64Hz非主导振荡模态对应的振荡源之一。计算G3、G4分别与G2在非主导频段下的能量频谱相似度如表1所示。
表1 相关机组能量频谱相似度系数
由表1可知,此时G4与G2的能量频谱具有较高相似性为0.878,而G3在该频段下与G2的能量频谱相似度仅为0.521,相关性较弱。由此可以判断得出,G2、G4为非主导模态下的振荡源,是由强迫振荡激发产生的区间振荡。虽然这两台机组的动态能量呈现增长趋势,机组阻尼为负,但不是主导模态下的振荡源。由此可以验证,仅依赖能量变化趋势无法准确辨识振荡源,而本发明所提方法能够有效准确地定位主导振荡模态下的振荡源机组。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种基于能量频谱的风电并网系统振荡源定位方法,其特征在于,包括以下步骤:
步骤1:采集系统各机组输出端口的电压和电流信息,获得各发电机组的动态能量随时间变化曲线;
步骤2:将所述动态能量随时间变化曲线中呈上升趋势的曲线对应的发电机组作为备选机组;
步骤3:获取备选机组中同步机组、含锁相环双馈风电机组的能量频谱和/或含虚拟惯量的双馈风电机组的类能量频谱;
步骤4:对所述能量频谱的幅值和/或类能量频谱的幅值进行归一化处理,对低频振荡频段上所有归一化后的所述能量频谱和/或类能量频谱的幅值进行标幺化处理,将主导振荡模态最大的机组作为振荡源参考机组;
步骤5:在振荡源参考机组的主导振荡模态±0.3Hz频段下,计算其余每一备选机组与所述振荡源参考机组的能量频谱相似度系数,根据相似度系数大小确定振荡源机组。
3.根据权利要求1所述的振荡源定位方法,其特征在于,所述步骤3包括,将各发电机组的动态能量在时域上进行傅里叶分解,得到频率成分f(k)对应的能量幅值A(k),k=1,2,3,4,…,N/2,N为离散信号的长度;
同步机组、含锁相环双馈风电机组的能量频谱均为A(k);
含虚拟惯量的双馈风电机组的类能量频谱为A(k)/f2(k)。
5.根据权利要求4所述的振荡源定位方法,其特征在于,步骤4中,对0.1Hz-2.5Hz频段上所有归一化后的能量频谱幅值和/或类能量频谱的幅值进行标幺化处理,得到每台机组各振荡模态能量频谱幅值和/或类能量频谱的幅值对相应机组所有能量频谱幅值之和的占比,所述占比最大的机组为振荡源参考机组。
7.根据权利要求1所述的振荡源定位方法,其特征在于,步骤5中,当与所述振荡源参考机组的能量频谱相似度系数大于等于0.7±0.05时,将该机组确定为振荡源机组。
8.一种基于能量频谱的风电并网系统振荡源定位系统,其特征在于,包括数据采集模块、动态能量计算模块、能量频谱分析模块和结果输出模块;
所述数据采集模块用于采集系统各机组输出端口的电压和电流信息,并将其传递给动态能量计算模块;
所述动态能量计算模块接收数据采集模块传递的数据,计算各发电机组的动态能量随时间变化趋势,并将其传递给能量频谱分析模块;
所述能量频谱分析模块接收各机组输出端口的动态能量数据,进行能量频谱计算,并通过对比各机组与参考振荡源机组能量频谱的相似度,定位系统振荡源;
所述结果输出模块用于输出主导振荡模态及对应的系统振荡源机组。
10.根据权利要求8所述的基于能量频谱的风电并网系统振荡源定位系统,其特征在于,所述能量频谱分析模块的工作过程为:
将所述动态能量随时间变化曲线中呈上升趋势的曲线对应的发电机组作为备选机组;
获取备选机组中同步机组、含锁相环双馈风电机组的能量频谱和/或含虚拟惯量的双馈风电机组的类能量频谱;
对所述能量频谱和/或类能量频谱的幅值进行标幺化处理,将主导振荡模态最大的机组作为振荡源参考机组;
在振荡源参考机组的主导振荡模态±0.3Hz频段下,计算其余每一备选机组与所述振荡源参考机组的能量频谱相似度系数,根据相似度系数大小确定振荡源机组。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910694996.3A CN110425092B (zh) | 2019-07-30 | 2019-07-30 | 一种基于能量频谱的风电并网系统振荡源定位系统及方法 |
US16/891,943 US11245263B2 (en) | 2019-07-30 | 2020-06-03 | System and method of locating oscillation sources of wind power integrated system based on energy spectrums |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910694996.3A CN110425092B (zh) | 2019-07-30 | 2019-07-30 | 一种基于能量频谱的风电并网系统振荡源定位系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110425092A CN110425092A (zh) | 2019-11-08 |
CN110425092B true CN110425092B (zh) | 2020-06-16 |
Family
ID=68411298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910694996.3A Active CN110425092B (zh) | 2019-07-30 | 2019-07-30 | 一种基于能量频谱的风电并网系统振荡源定位系统及方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11245263B2 (zh) |
CN (1) | CN110425092B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110425092B (zh) * | 2019-07-30 | 2020-06-16 | 华北电力大学 | 一种基于能量频谱的风电并网系统振荡源定位系统及方法 |
CN110417056B (zh) * | 2019-07-30 | 2020-12-18 | 华北电力大学 | 双馈风机并网系统的设备级振荡源定位方法及其装置 |
CN110826927B (zh) * | 2019-11-11 | 2020-07-10 | 国网山东省电力公司日照供电公司 | 一种风电场可用惯量评估方法 |
CN111769596B (zh) * | 2020-07-15 | 2022-03-01 | 华北电力大学 | 一种基于附加能量支路的双馈风电场控制方法及系统 |
CN111969622B (zh) * | 2020-08-13 | 2022-07-05 | 清华大学 | 电力系统等幅振荡的频率和幅值参数解析计算方法及装置 |
CN112308470B (zh) * | 2020-12-28 | 2021-04-20 | 北京隆普智能科技有限公司 | 一种风电并网频率响应异常监测方法及系统 |
CN113866565B (zh) * | 2021-10-22 | 2024-03-29 | 福州大学 | 一种基于svmd的风能渗透型配电网事件检测方法 |
CN114611987B (zh) * | 2022-03-25 | 2022-12-09 | 南通大学 | 一种并网型风电场的局域振荡的辨识方法 |
CN114792984B (zh) * | 2022-04-18 | 2024-07-05 | 国网山东省电力公司莱芜供电公司 | 一种次/超同步振荡源的快速定位方法 |
CN115498652B (zh) * | 2022-05-30 | 2024-04-30 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 基于cps的综合能源园区电能质量控制方法及系统 |
CN114966296B (zh) * | 2022-08-01 | 2022-12-02 | 锦浪科技股份有限公司 | 一种振荡源定位方法及装置 |
CN115015697B (zh) * | 2022-08-09 | 2022-10-11 | 中国电力科学研究院有限公司 | 一种新能源的宽频振荡扰动源定位方法及系统 |
CN115483691B (zh) * | 2022-11-15 | 2023-02-28 | 华北电力大学 | 储能协调漂浮式风电机组的控制方法及装置 |
CN116203350B (zh) * | 2023-03-07 | 2024-02-27 | 西南交通大学 | 一种高水电占比电网中频率振荡源定位方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103901882A (zh) * | 2014-04-15 | 2014-07-02 | 北京交通大学 | 一种列车动力系统在线监控故障诊断的系统及其方法 |
CN109217362A (zh) * | 2018-08-22 | 2019-01-15 | 华北电力大学 | 一种双馈风机并网系统低频振荡扰动源定位系统及方法 |
CN109283427A (zh) * | 2018-08-21 | 2019-01-29 | 华北电力大学 | 含双馈风电机组的互联电网低频扰动源定位系统及方法 |
CN110011357A (zh) * | 2019-04-29 | 2019-07-12 | 东南大学 | 一种双馈风电场中强迫振荡扰动源定位方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8471610B2 (en) * | 2011-09-29 | 2013-06-25 | Electric Power Research Institute, Inc. | Application of phase-locked loop (PLL) in oscillation monitoring for interconnected power systems |
JP6348862B2 (ja) * | 2015-03-30 | 2018-06-27 | 株式会社日立製作所 | 系統安定化制御装置および電力系統制御システム |
CN106286129B (zh) | 2016-10-12 | 2021-04-06 | 北京金风科创风电设备有限公司 | 风力发电机组及其控制方法 |
US10063161B2 (en) * | 2016-10-18 | 2018-08-28 | Abb Schweiz Ag | Active neutral point clamped converter control system and method |
CN110425092B (zh) * | 2019-07-30 | 2020-06-16 | 华北电力大学 | 一种基于能量频谱的风电并网系统振荡源定位系统及方法 |
-
2019
- 2019-07-30 CN CN201910694996.3A patent/CN110425092B/zh active Active
-
2020
- 2020-06-03 US US16/891,943 patent/US11245263B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103901882A (zh) * | 2014-04-15 | 2014-07-02 | 北京交通大学 | 一种列车动力系统在线监控故障诊断的系统及其方法 |
CN109283427A (zh) * | 2018-08-21 | 2019-01-29 | 华北电力大学 | 含双馈风电机组的互联电网低频扰动源定位系统及方法 |
CN109217362A (zh) * | 2018-08-22 | 2019-01-15 | 华北电力大学 | 一种双馈风机并网系统低频振荡扰动源定位系统及方法 |
CN110011357A (zh) * | 2019-04-29 | 2019-07-12 | 东南大学 | 一种双馈风电场中强迫振荡扰动源定位方法 |
Non-Patent Citations (1)
Title |
---|
项晓强.双馈风机并网次同步振荡机理与控制措施研究.《华北电力大学硕士学位论文》.2019,第7-36页. * |
Also Published As
Publication number | Publication date |
---|---|
US11245263B2 (en) | 2022-02-08 |
US20210036518A1 (en) | 2021-02-04 |
CN110425092A (zh) | 2019-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110425092B (zh) | 一种基于能量频谱的风电并网系统振荡源定位系统及方法 | |
CN109217362B (zh) | 一种双馈风机并网系统低频振荡扰动源定位系统及方法 | |
CN110943451B (zh) | 一种基于扰动数据的系统分区域惯量评估的方法 | |
CN110518631B (zh) | 一种直驱风电机组的稳定性评估方法及系统 | |
CN101119095A (zh) | 基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统 | |
CN110417054B (zh) | 一种双馈风电并网系统稳定性调整系统及方法 | |
CN110120686A (zh) | 一种基于电力系统在线惯量估计的新能源承载力预警方法 | |
CN111769573B (zh) | 基于能量补偿的直驱风机次/超频振荡抑制方法及系统 | |
CN109861296B (zh) | 一种孤岛虚拟同步发电机系统的惯性和阻尼辨识方法 | |
CN114006383B (zh) | 含vsg的双馈风机并网系统建模与次同步振荡分析的方法 | |
CN111009921A (zh) | 基于奈奎斯特稳定判据的双馈风机并网系统振荡分析方法 | |
CN105259414A (zh) | 一种基于逆变器的电网阻抗在线检测方法 | |
Fanglei et al. | Estimating maximum penetration level of renewable energy based on frequency stability constrains in power grid | |
CN111082436B (zh) | 一种基于锁相一致的直驱风电场振荡抑制方法及系统 | |
CN108258725A (zh) | 一种基于等效功角同调的双馈风电机动态等值方法 | |
Xiang et al. | Equivalent modeling of wind farms considering wake effects and analysis of probability of sub-synchronous oscillations | |
Long et al. | Analysis of sub-synchronous oscillations on hybrid wind farms with DFIG and PMSG | |
CN112861326A (zh) | 一种基于量测的新能源电网发电机阻尼评估装置及方法 | |
Melhem et al. | Frequency support and stability analysis for an integrated power system with wind farms | |
CN112186767A (zh) | 含高比例可再生能源的海岛微电网频率稳定的优化控制方法 | |
CN116937671B (zh) | 一种孤岛交流微电网的恒频分布式控制方法 | |
Li et al. | An improved quantitative analysis method of oscillation mode for DFIG-based wind power base with LCC-HVDC transmission considering frequency coupling characteristic | |
CN117154760A (zh) | 计及电压特性的大电网频率响应模型的构建方法与系统 | |
Liu et al. | The influence of wind/PV generation uncertainties on power system low frequency oscillation | |
CN116780502B (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 |