CN106953317B - 电力设备的次/超同步耦合阻抗(导纳)模型的辨识方法 - Google Patents
电力设备的次/超同步耦合阻抗(导纳)模型的辨识方法 Download PDFInfo
- Publication number
- CN106953317B CN106953317B CN201710153029.7A CN201710153029A CN106953317B CN 106953317 B CN106953317 B CN 106953317B CN 201710153029 A CN201710153029 A CN 201710153029A CN 106953317 B CN106953317 B CN 106953317B
- Authority
- CN
- China
- Prior art keywords
- supersynchronous
- subsynchronous
- voltage
- sub
- impedance
- 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
- 238000012850 discrimination method Methods 0.000 title abstract 3
- 238000005259 measurement Methods 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 39
- 230000005284 excitation Effects 0.000 claims abstract description 26
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 238000002847 impedance measurement Methods 0.000 claims abstract description 8
- 230000000295 complement effect Effects 0.000 claims abstract description 7
- 230000008878 coupling Effects 0.000 claims description 65
- 238000010168 coupling process Methods 0.000 claims description 65
- 238000005859 coupling reaction Methods 0.000 claims description 65
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000012360 testing method Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 9
- 230000001360 synchronised effect Effects 0.000 claims description 9
- 230000021615 conjugation Effects 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000000691 measurement method Methods 0.000 description 6
- 238000002347 injection Methods 0.000 description 5
- 239000007924 injection Substances 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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]
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种电力设备的次/超同步耦合阻抗(导纳)模型的辨识方法,该阻抗模型的辨识方法包括:提供电力设备的次/超同步耦合阻抗模型;将所述电力设备接入电网中,并向所述电网施加次/超同步激励,并基于所述次/超同步耦合阻抗模型获取所述电力设备的次同步和超同步电压分量和电流分量;改变所述次/超同步激励的参数或外电网的参数,以获取所述电力设备的多组测量结果,并根据所述多组测量结果生成阻抗量测方程;对所述阻抗量测方程进行预设数值优化方法得到所述次/超同步耦合阻抗模型的阻抗参数。本发明具有如下优点:表征了两种相对工频的互补频率动态之间的不可分割性,能更好地分析和解释电力***的次/超同步谐振和振荡问题。
Description
技术领域
本发明涉及电力***分析与控制技术领域,尤其涉及一种电力设备的次/超同步耦合阻抗(导纳)模型的辨识方法。
背景技术
在现代电力***电力电子化程度日益加深的趋势下,电力电子变流技术不仅改变电网一次设备的形态,而且其特有的多尺度控制机制和低惯性宽频带响应特性,正显著改变电力***的动态特征,带来新的稳定性分析与问题。近年来,尤其以风电、光伏等变流电源大规模接入电网引起的新型次/超同步振荡(Sub Synchronous Oscillation,SSO)问题尤为突出。此外,变流器式恒功率负载在次同步频率上体现的负电阻特性,相邻变流器的锁相环(PLL)回路耦合,以及变流器控制参与滤波电路和电网构成的串/并联谐振回路,都可能激发次/超同步振荡风险。这种海量变流器装置(风电、光伏、逆变器等)-电网相互作用引发的新型次/超同步振荡严重威胁现代电网的设备安全、***稳定和用电质量,在我国已成为制约风光等新能源大规模消纳的瓶颈因素。
针对电力***电力电子化带来的稳定性问题的分析方法大多沿用传统的电磁暂态仿真、特征值分析或简化条件下的阻抗分析方法。目前广泛应用的电磁暂态仿真软件(如PSCAD, EMTP)只支持单一步长(固定速率)的建模与分析,难以兼顾精度与效率以实现多变流器、复杂交直流网络、海量工况下的时域仿真;传统的特征值分析方法面临维数灾问题,难以适用于包括数千台变流器、数万个节点的实际大电网。
发明内容
本发明旨在至少解决上述技术问题之一。
为此,本发明的一个目的在于提出一种电力设备的次/超同步耦合阻抗模型的辨识方法。
为了实现上述目的,本发明的实施例公开了一种电力设备的次/超同步耦合阻抗模型的辨识方法,包括以下步骤:提供电力设备的次/超同步耦合阻抗模型,其中,所述次/超同步耦合阻抗模型包括电压向量、电流向量和阻抗参数;将所述电力设备接入电网中,并向所述电网施加次/超同步激励,并基于所述次/超同步耦合阻抗模型获取所述电力设备的次同步和超同步电压分量和电流分量;改变所述次/超同步激励的参数或外电网的参数,以获取所述电力设备的多组测量结果,并根据所述多组测量结果生成阻抗量测方程;对所述阻抗量测方程进行预设数值优化方法得到所述次/超同步耦合阻抗模型的阻抗参数。
进一步地,向所述电网施加次/超同步激励,并基于所述次/超同步耦合阻抗模型获取所述电力设备的次同步和超同步电压分量和电流分量进一步包括:向所述电网注入一定次同步和/或超同步频率电流源和/或电压源,并基于所述次/超同步耦合阻抗模型在所述电力设备的端口电压和端口电流中出现次同步分量和超同步分量,通过信号处理得到所述次同步和超同步电压分量和电流分量。
进一步地,改变所述次/超同步激励的参数包括改变注入所述电流源和所述电流源的幅值。
根据本发明实施例的电力设备的次/超同步耦合阻抗模型的辨识方法,表征了两种相对工频的互补频率动态之间的不可分割性,能更好地分析和解释电力***的次/超同步谐振和振荡问题。
为此,本发明的一个目的在于提出一种电力设备的次/超同步耦合导纳模型的辨识方法。
为了实现上述目的,本发明的实施例公开了一种电力设备的次/超同步耦合导纳模型的辨识方法,包括以下步骤:提供电力设备的次/超同步耦合导纳模型,其中,所述次/超同步耦合导纳模型包括电压向量、电流向量和导纳参数;将所述电力设备接入电网中,并向所述电网施加次/超同步激励,并基于所述次/超同步耦合导纳模型获取所述电力设备的次同步和超同步电压分量和电流分量;改变所述次/超同步激励的参数或外电网的参数,以获取所述电力设备的多组测量结果,并根据所述多组测量结果生成导纳量测方程;对所述导纳量测方程进行预设数值优化方法得到所述次/超同步耦合导纳模型的导纳参数。
进一步地,向所述电网施加次/超同步激励,并基于所述次/超同步耦合导纳模型获取所述电力设备的次同步和超同步电压分量和电流分量进一步包括:向所述电网注入一定次同步和/或超同步频率电流源和/或电压源,并基于所述次/超同步耦合导纳模型在所述电力设备的端口电压和端口电流中出现次同步分量和超同步分量,通过信号处理得到所述次同步和超同步电压分量和电流分量。
进一步地,改变所述次/超同步激励的参数包括改变注入所述电流源和所述电流源的幅值。
根据本发明实施例的电力设备的次/超同步耦合导纳模型的辨识方法,表征了两种相对工频的互补频率动态之间的不可分割性,能更好地分析和解释电力***的次/超同步谐振和振荡问题。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明实施例的电力设备的次/超同步耦合阻抗模型的辨识方法的流程图;
图2是本发明一个实施例的电力设备并联接入电网的示意图;
图3是本发明一个实施例的电力设备串联接入电网的示意图;
图4是本发明一个实施例的次/超同步耦合阻抗(导纳)模型的通用辨识方法的原理图;
图5是本发明一个实施例的次/超同步耦合阻抗(导纳)模型的实用化测量方法的原理图(测量点位于电力设备端口);
图6是本发明一个实施例的次/超同步耦合阻抗(导纳)模型的实用化测量方法的原理图(测量点位于连接阻抗靠近电网侧)
图7是本发明一个实施例的次/超同步耦合阻抗(导纳)模型的简化测量方法的原理图。
图8是本发明实施例的电力设备的次/超同步耦合导纳模型的辨识方法的流程图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
参照下面的描述和附图,将清楚本发明的实施例的这些和其他方面。在这些描述和附图中,具体公开了本发明的实施例中的一些特定实施方式,来表示实施本发明的实施例的原理的一些方式,但是应当理解,本发明的实施例的范围不受此限制。相反,本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
以下结合附图描述本发明。
图1是本发明实施例的电力设备的次/超同步耦合阻抗模型的辨识方法的流程图。如图 1所示,根据本发明实施例的电力设备的次/超同步耦合阻抗模型的辨识方法,包括以下步骤:
S1:提供电力设备的次/超同步耦合阻抗模型。其中,次/超同步耦合阻抗模型包括电压向量、电流向量和阻抗参数。
对于如图2所示的并联接入电网的三相电力设备(如汽轮发电机、风力发电机、光伏发电、静止无功补偿器、三相变流器负载等)和图3所示的串联接入电网的三相电力设备(如输电线路、变压器、串联型柔性输电控制器、直流输电等),其正序次/超同步耦合阻抗模型定义为:
式中,和分别表示电压和电流相量(以A相基频电压分量为相角参考值),z表示阻抗参数,下标sub和sup分别表示次同步和超同步分量,上标*表示相量取共轭,下同。
需要说明的是,本发明的实施例中所提到的耦合的次/超同步分量,其频率均满足互补关系,即
fsub+fsup=2f1 (2)
式中,f1是基波频率。
或,不用共轭的阻抗模型表达方式为:
式中,θu和θi表示以A相基频电压分量为相角参考值时,和的相角值;
或,共轭放在不同地方的阻抗模型表达方式为:
以上是正序次同步和超同步分量的耦合,研究发现,负序的次、超同步分量和正序的高频分量之间也存在耦合关系。按照公式(2),负序分量的频率为负值,与其互补的频率将大于两倍基波频率,即为正序的高频分量,因此次/超同步耦合阻抗模型定义为
式中,和分别表示电压和电流相量(以A相基频电压分量为相角参考值),z表示阻抗参数,下标p和n分别表示正序和负序分量,其频率也满足互补关系:
fp+fn=2f1(fp>0,fn<0) (6)
S2:将电力设备接入电网中,并向电网施加次/超同步激励,并基于次/超同步耦合阻抗模型获取电力设备的次同步和超同步电压分量和电流分量。
在本发明的一个实施例中,向电网施加次/超同步激励,并基于次/超同步耦合阻抗模型获取电力设备的次同步和超同步电压分量和电流分量进一步包括:
向电网注入一定次同步和/或超同步频率电流源和/或电压源,并基于次/超同步耦合阻抗模型在电力设备的端口电压和端口电流中出现次同步分量和超同步分量,通过信号处理得到次同步和超同步电压分量和电流分量。
在本发明的一个实施例中,改变次/超同步激励的参数包括改变注入电流源和电流源的幅值。
具体地,如图4所示,三相电力设备接入电网中,设采取某种试验手段(对电网进行次超同步激励),如向电网并联和/或串联注入一定次同步和/或超同步频率电流源和/或电压源(即图4中的扰动信号,ud1是扰动电压源,id2是扰动电流源),导致被测设备端口电压和端口电流中出现明显的次同步和超同步分量,通过信号处理方式(如滤波)将次同步和超同步电压和电流分量分别“滤取”出来,设为usub,usup,isub,isup(下标sub和sup分别表示次同步和超同步分量下同)。
以A相基频电压为基准(相位角设为0),将三相次/超电压和电流表示为相量形式,求取其正序分量,设为根据式(1)则有:
Itestz=Utest (7)
式中各矩阵或向量的定义如下:
是可以实测的电流构成的矩阵,上标*表示相量取共轭。
是可以实测的电压构成的向量。
z=[z11 z12 z21 z22]T,是将阻抗矩阵展开形成的阻抗参数列向量。
S3:改变次/超同步激励的参数或外电网的参数,以获取电力设备的多组测量结果,并根据多组测量结果生成阻抗量测方程。
具体地,当改变注入电压和/电流源的幅值或外电网的参数,使得待测设备端部电压和电流发生变化时,皆可以得到多组测量结果,将它们都整理成如式(7)所示的形式,并组合起来构成如下量测方程。
式中,上标(1)和(n)分别表示第一组和第n组测量结果。
如果得到了n组测量结果,则式中是2n×4维矩阵,是2n维向量。
S4:对阻抗量测方程进行预设数值优化方法得到次/超同步耦合阻抗模型的阻抗参数。
具体地,可以通过最小二乘方法或其他数值优化方法得到阻抗参数,从而实现对阻抗模型的精确辨识。
实际使用中,可以对前述通用辨识方法进行简化,得到次/超同步耦合阻抗(导纳)模型的实用化测量方法。如图5所示,扰动信号源(即前述的次同步和/或超同步频率电流源和/或电压源)经过一定的连接阻抗,注入被测设备的端口。对被测设备端口的电压和电流进行测量,通过信号处理方式(如滤波)将次同步和超同步电压和电流分量分别“滤取”出来。改变扰动信号源的大小和/或相位,再次测量,共得到两组测量结果。
根据公式(8),可直接求得阻抗矩阵参数如下:
如果不方便在被测设备端口进行测量且已知连接阻抗值,还可以将测量点移至连接阻抗靠近电网侧,如图6所示。计算阻抗时只需对电压测量结果(设为)进行预处理即可,如下式所示:
式中,Ztest是连接阻抗值。
再将电流测量结果和预处理后的电压测量结果代入公式(9)便可得到阻抗/导纳矩阵参数。
进一步地,如果可以将扰动信号直接施加在被测设备的端口,即连接阻抗为零,则可得到次/超同步耦合阻抗(导纳)模型的简化测量方法,如图7所示。
若扰动源为单一的次同步频率电压源,此时机端电压也只含有次同步频率信号,即改变扰动源为单一的超同步频率电压源,再次测量,此时机端电压也只含有超同步频率信号,即
根据公式(8),阻抗矩阵参数公式简化为
图8是本发明实施例的电力设备的次/超同步耦合阻抗模型的辨识方法的流程图。如图 8所示,本发明的实施例还提供了一种电力设备的次/超同步耦合导纳模型的辨识方法,包括以下步骤:
A:提供电力设备的次/超同步耦合导纳模型,其中,次/超同步耦合导纳模型包括电压向量、电流向量和导纳参数。
具体地,正序次/超同步耦合导纳模型表示方式为:
和分别表示电压和电流相量(以A相基频电压分量为相角参考值),y表示导纳参数,下标sub和sup分别表示次同步和超同步分量,上标*表示相量取共轭,下同。
次/超同步耦合阻抗模型定义为:
B:将电力设备接入电网中,并向电网施加次/超同步激励,并基于次/超同步耦合导纳模型获取电力设备的次同步和超同步电压分量和电流分量。
在本发明的一个实施例中,向电网施加次/超同步激励,并基于次/超同步耦合导纳模型获取电力设备的次同步和超同步电压分量和电流分量进一步包括:
向电网注入一定次同步和/或超同步频率电流源和/或电压源,并基于次/超同步耦合导纳模型在电力设备的端口电压和端口电流中出现次同步分量和超同步分量,通过信号处理得到次同步和超同步电压分量和电流分量。
在本发明的一个实施例中,改变次/超同步激励的参数包括改变注入电流源和电流源的幅值。
以A相基频电压为基准(相位角设为0),将三相次/超电压和电流表示为相量形式,求取其正序分量,设为根据式(12)则有:
Itest=Utesty (14)
式中各矩阵或向量的定义如下:
是可以实测的电流构成的矩阵,上标*表示相量取共轭。
是可以实测的电压构成的向量。
y=[y11 y12 y21 y22]T,是将导纳矩阵展开形成的导纳参数列向量。
C:改变次/超同步激励的参数或外电网的参数,以获取电力设备的多组测量结果,并根据多组测量结果生成导纳量测方程。
具体地,导纳量测方程为:
D:对导纳量测方程进行预设数值优化方法得到次/超同步耦合导纳模型的导纳参数。
具体地,实际使用中,可以对前述通用辨识方法进行简化,得到次/超同步耦合阻抗(导纳)模型的实用化测量方法。如图5所示,扰动信号源(即前述的次同步和/或超同步频率电流源和/或电压源)经过一定的连接阻抗,注入被测设备的端口。对被测设备端口的电压和电流进行测量,通过信号处理方式(如滤波)将次同步和超同步电压和电流分量分别“滤取”出来。改变扰动信号源的大小和/或相位,再次测量,共得到两组测量结果。
导纳矩阵参数如下:
计算阻抗和导纳时只需对电压测量结果(设为)进行预处理即可,如下式所示:
再将电流测量结果和预处理后的电压测量结果代入公式(16)便可得到阻抗/导纳矩阵参数。
进一步地,如果可以将扰动信号直接施加在被测设备的端口,即连接阻抗为零,则可得到次/超同步耦合阻抗(导纳)模型的简化测量方法,如图7所示。
若扰动源为单一的次同步频率电压源,此时机端电压也只含有次同步频率信号,即改变扰动源为单一的超同步频率电压源,再次测量,此时机端电压也只含有超同步频率信号,即
根据公式(15),导纳矩阵参数公式简化为:
另外,本发明实施例的电力设备的次/超同步耦合模型的辨识方法的其它构成以及作用对于本领域的技术人员而言都是已知的,为了减少冗余,不做赘述。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同限定。
Claims (6)
1.一种电力设备的次/超同步耦合阻抗模型的辨识方法,其特征在于,包括以下步骤:
提供电力设备的次/超同步耦合阻抗模型,其中,所述次/超同步耦合阻抗模型包括电压向量、电流向量和阻抗参数,所述次/超同步耦合阻抗模型为:
其中,和分别表示电压和电流相量,z表示阻抗参数,下标sub和sup分别表示次同步和超同步分量,上标*表示相量取共轭,次/超同步分量频率均满足互补关系,即
fsub+fsup=2f1
式中,f1是基波频率;
将所述电力设备接入电网中,并向所述电网施加次/超同步激励,并基于所述次/超同步耦合阻抗模型获取所述电力设备的次同步和超同步电压分量和电流分量,以A相基频电压为基准,相位角设为0,将三相次/超同步电压和电流表示为相量形式,求取其正序分量,设为则有:
Itestz=Utest;
式中各矩阵或向量的定义如下:
是实测的电流构成的矩阵,上标*表示相量取共轭;
是实测的电压构成的向量;
z=[z11 z12 z21 z22]T,是将阻抗矩阵展开形成的阻抗参数列向量;
改变所述次/超同步激励的参数或外电网的参数,以获取所述电力设备的多组测量结果,并根据所述多组测量结果生成阻抗量测方程,其中,所述阻抗量测方程为:
其中,上标(1)和(n)分别表示第一组和第n组测量结果,如果得到了n组测量结果,则式中是2n×4维矩阵,是2n维向量;
对所述阻抗量测方程进行预设数值优化方法得到所述次/超同步耦合阻抗模型的阻抗参数;所述预设数值优化方法为通过最小二乘方法得到阻抗参数,从而实现对阻抗模型的精确辨识。
2.根据权利要求1所述的电力设备的次/超同步耦合阻抗模型的辨识方法,其特征在于,向所述电网施加次/超同步激励,并基于所述次/超同步耦合阻抗模型获取所述电力设备的次同步和超同步电压分量和电流分量进一步包括:
向所述电网注入一定次/超同步频率电流源和/或电压源,并基于所述次/超同步耦合阻抗模型在所述电力设备的端口电压和端口电流中出现次同步分量和超同步分量,通过信号处理得到所述次同步和超同步电压分量和电流分量。
3.根据权利要求2所述的电力设备的次/超同步耦合阻抗模型的辨识方法,其特征在于,改变所述次/超同步激励的参数包括改变注入所述电流源的相位和所述电流源的幅值。
4.一种电力设备的次/超同步耦合导纳模型的辨识方法,其特征在于,包括以下步骤:
提供电力设备的次/超同步耦合导纳模型,其中,所述次/超同步耦合导纳模型包括电压向量、电流向量和导纳参数,所述次/超同步耦合导纳模型为:
其中,和分别表示电压和电流相量,y表示导纳参数,下标sub和sup分别表示次同步和超同步分量,上标*表示相量取共轭;
次/超同步耦合导纳模型定义为:
下标p和n分别表示正序和负序分量
将所述电力设备接入电网中,并向所述电网施加次/超同步激励,并基于所述次/超同步耦合导纳模型获取所述电力设备的次同步和超同步电压分量和电流分量,以A相基频电压为基准,相位角设为0,将三相次/超同步电压和电流表示为相量形式,求取其正序分量,设为则有:
Itest=Utesty;
式中各矩阵或向量的定义如下:
是实测的电流构成的矩阵,上标*表示相量取共轭;
是实测的电压构成的向量;
y=[y11 y12 y21 y22]T,是将导纳矩阵展开形成的导纳参数列向量;
改变所述次/超同步激励的参数或外电网的参数,以获取所述电力设备的多组测量结果,并根据所述多组测量结果生成导纳量测方程,其中,所述导纳量测方程为:
其中,上标(1)和(n)分别表示第一组和第n组测量结果,如果得到了n组测量结果,则式中是2n×4维矩阵,是2n维向量;
对所述导纳量测方程进行预设数值优化方法得到所述次/超同步耦合导纳模型的导纳参数。
5.根据权利要求4所述的电力设备的次/超同步耦合导纳模型的辨识方法,其特征在于,向所述电网施加次/超同步激励,并基于所述次/超同步耦合导纳模型获取所述电力设备的次同步和超同步电压分量和电流分量进一步包括:
向所述电网注入一定次/超同步频率电流源和/或电压源,并基于所述次/超同步耦合导纳模型在所述电力设备的端口电压和端口电流中出现次同步分量和超同步分量,通过信号处理得到所述次同步和超同步电压分量和电流分量。
6.根据权利要求5所述的电力设备的次/超同步耦合导纳模型的辨识方法,其特征在于,改变所述次/超同步激励的参数包括改变注入所述电流源的相位和所述电流源的幅值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710153029.7A CN106953317B (zh) | 2017-03-15 | 2017-03-15 | 电力设备的次/超同步耦合阻抗(导纳)模型的辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710153029.7A CN106953317B (zh) | 2017-03-15 | 2017-03-15 | 电力设备的次/超同步耦合阻抗(导纳)模型的辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106953317A CN106953317A (zh) | 2017-07-14 |
CN106953317B true CN106953317B (zh) | 2019-10-29 |
Family
ID=59471907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710153029.7A Active CN106953317B (zh) | 2017-03-15 | 2017-03-15 | 电力设备的次/超同步耦合阻抗(导纳)模型的辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106953317B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107402334B (zh) * | 2017-07-10 | 2019-08-09 | 清华大学 | 电力***次/超同步耦合阻抗模型的频域辨识方法与*** |
CN108134399B (zh) * | 2018-01-12 | 2020-10-13 | 清华大学 | 网侧次同步阻尼控制器全工况优化方法及装置 |
CN110866338B (zh) * | 2019-11-13 | 2021-08-17 | 清华大学 | 基于耦合阻抗模型求取视在阻抗模型的方法及装置 |
CN111679125B (zh) * | 2020-06-04 | 2021-07-20 | 北京交通大学 | 一种电力***振荡辨识的方法和装置 |
CN112014642B (zh) * | 2020-08-03 | 2021-10-26 | 清华大学 | 静止坐标系下电网频率耦合阻抗模型聚合计算方法和装置 |
CN112865181B (zh) * | 2021-03-02 | 2022-07-29 | 国网冀北电力有限公司电力科学研究院 | 基于端口阻抗特性的光伏逆变器参数辨识方法及装置 |
CN113962171B (zh) * | 2021-10-13 | 2024-04-16 | 西安交通大学 | 一种有损地面上输电线路的高频耦合方法 |
CN114237045B (zh) * | 2021-11-29 | 2022-08-02 | 哈尔滨工业大学 | 一种无传感器式压电驱动闭环控制方法 |
CN116264396A (zh) * | 2021-12-15 | 2023-06-16 | 新疆金风科技股份有限公司 | 风力发电机组、风电变流器的控制方法和控制装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5473528A (en) * | 1991-09-18 | 1995-12-05 | Kabushiki Kaisha Toshiba | Parallel connection of different types of AC power supplies of differing capacities |
CN103259472A (zh) * | 2013-04-19 | 2013-08-21 | 上海交通大学 | 大型发电机组次同步振荡的反振荡自适应抑制*** |
CN103986171A (zh) * | 2014-06-03 | 2014-08-13 | 贵州电力试验研究院 | 抑制次同步振荡的ssdc和sedc协调控制优化方法 |
CN106526328A (zh) * | 2016-12-08 | 2017-03-22 | 浙江大学 | 一种适用于电网及联网设备的广义阻抗测量与计算方法 |
-
2017
- 2017-03-15 CN CN201710153029.7A patent/CN106953317B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5473528A (en) * | 1991-09-18 | 1995-12-05 | Kabushiki Kaisha Toshiba | Parallel connection of different types of AC power supplies of differing capacities |
CN103259472A (zh) * | 2013-04-19 | 2013-08-21 | 上海交通大学 | 大型发电机组次同步振荡的反振荡自适应抑制*** |
CN103986171A (zh) * | 2014-06-03 | 2014-08-13 | 贵州电力试验研究院 | 抑制次同步振荡的ssdc和sedc协调控制优化方法 |
CN106526328A (zh) * | 2016-12-08 | 2017-03-22 | 浙江大学 | 一种适用于电网及联网设备的广义阻抗测量与计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106953317A (zh) | 2017-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106953317B (zh) | 电力设备的次/超同步耦合阻抗(导纳)模型的辨识方法 | |
Lyu et al. | Harmonic state-space based small-signal impedance modeling of a modular multilevel converter with consideration of internal harmonic dynamics | |
CN107121609B (zh) | 一种基于prbs扰动注入的电网阻抗在线辨识方法及装置 | |
CN107608933B (zh) | 一种基于二次侧施加扰动的逆变器广义阻抗测量方法 | |
CN107688722B (zh) | 电压源型变流器的导纳模型和阻抗模型获取方法及装置 | |
CN111239491B (zh) | 采用实物控制器扰动注入的广义阻抗实时实验测量方法 | |
CN109802406A (zh) | 一种分析柔性直流输电接入***谐振稳定性的方法 | |
CN110739678A (zh) | 一种并网换流器串联虚拟阻抗的控制方法 | |
CN108918973A (zh) | 基于upqc的高精度阻抗测量装置及其切换控制方法 | |
CN108767899A (zh) | 一种低压微电网小扰动稳定性分析方法 | |
CN109962479A (zh) | 一种基于交替迭代优化的同步调相机电气参数分布联合辨识方法 | |
CN108155643A (zh) | 一种基于滑模观测器的单相电网电压参数的鲁棒估计方法 | |
CN114325097A (zh) | 一种基于二次侧扰动注入的双馈型发电设备阻抗测量方法 | |
Zong et al. | Accurate aggregated modelling of wind farm systems in modified sequence domain for stability analysis | |
CN111157798A (zh) | 一种基于实时仿真机和实物控制器的阻抗测量*** | |
CN109902347B (zh) | 基于基波阻抗辨识的并网***短路比测量方法及装置 | |
Xue et al. | A comprehensive study on impedance models of grid-tied voltage-source converters | |
CN207472983U (zh) | 一种基于prbs扰动注入的电网阻抗在线辨识装置 | |
CN102570950A (zh) | 一种发电机机端次同步阻尼控制***及次同步阻尼控制方法 | |
CN114325113A (zh) | 基于采样信号扰动叠加的逆变器正负序阻抗测量方法 | |
Wang et al. | Harmonic resonance analysis in high‐renewable‐energy‐penetrated power systems considering frequency coupling | |
Shi et al. | Transient response comparison of virtual oscillator controlled and droop controlled three‐phase inverters under load changes | |
CN109839540A (zh) | 基于rpwm的在线单相并网逆变器弱网阻抗测量方法 | |
CN108347055B (zh) | 一种并网滤波电感器参数评测电路及其控制方法 | |
Zong et al. | Three-port impedance model and validation of VSCs for stability analysis |
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 |