CN105140935A - 具有高阻尼值的电力***广域时滞pid阻尼控制器设计方法 - Google Patents

具有高阻尼值的电力***广域时滞pid阻尼控制器设计方法 Download PDF

Info

Publication number
CN105140935A
CN105140935A CN201510492671.9A CN201510492671A CN105140935A CN 105140935 A CN105140935 A CN 105140935A CN 201510492671 A CN201510492671 A CN 201510492671A CN 105140935 A CN105140935 A CN 105140935A
Authority
CN
China
Prior art keywords
power system
pid
damping
time lag
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.)
Granted
Application number
CN201510492671.9A
Other languages
English (en)
Other versions
CN105140935B (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201510492671.9A priority Critical patent/CN105140935B/zh
Publication of CN105140935A publication Critical patent/CN105140935A/zh
Application granted granted Critical
Publication of CN105140935B publication Critical patent/CN105140935B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • 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
    • 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
    • Y04S10/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units

Landscapes

  • Feedback Control In General (AREA)

Abstract

电力***低频振荡广域阻尼控制器设计方法,针对广域电力***在区域电网互联过程中出现的低频振荡问题，提出了具有高阻尼值的广域时滞PID阻尼控制器的设计方法，包括：确定电力***低频振荡模式，选取反应区间低频振荡特性的反馈控制信号以及合适的阻尼控制执行器，建立广域电力***的传递函数数学模型，设计适用于时滞PMU信号反馈控制的电力***PID控制器，计算能够确保电力***稳定的可能的PID参数分布范围，选取使电力***具有高阻尼值的PID参数。本发明不仅适用于广域反馈PMU信号存在时滞的情形、抑制区间低频振荡效果突出，而且所设计的控制器具有简洁的结构、易于工程实现。

Description

具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法
技术领域：
本发明涉及电力***低频振荡广域阻尼控制器设计方法，属于电力***稳定控制技术领域。
背景技术：
随着大区电网的互联，电力***的规模不断扩大，低频振荡问题也日益严重。如何确保大区电网特高压互联下的电力***具有较高的阻尼系数，成为现代电力***安全稳定运行亟待解决的控制难题。
由于大区电网互联后，发生的低频振荡可能同时涉及多个区域电网，分布面非常广，影响也很大。传统的采用本地信号作为反馈信号的电力***稳定控制器(PSS)，受控制信号的可观性限制，在抑制区间低频振荡方面效果非常有限。基于GPS授时技术的相量测量单元(PMU)，使电力***运行状态的同步测量不再是个难题，目前基于PMU的电力***广域测量***(WAMS)也正在成形与完善中。利用PMU同步相量数据进行广域电力***低频振荡阻尼控制，较传统的基于本地信号的PSS控制具有显著的优势，在大量文献中已有明确的结论。但是，PMU同步相量的广域远距离传输也带来了控制信号时滞问题，若不妥善处理时滞影响，阻尼控制器不仅不能起到抑制低频振荡的作用，反而可能进一步恶化电力***的稳定性。
在过去的几十年里，PID控制器在工业控制中得到了广泛应用，工业过程控制中95％以上的控制回路都具有PID结构，并且许多高级控制都是以PID控制为基础的，电力***亦不例外。PID控制器结构和算法简单，应用广泛，但参数整定方法复杂，若要考虑时滞的影响，PID控制器的参数整定将难上加难。
发明内容：
本发明针对广域电力***在区域电网互联过程中出现的低频振荡问题，提出了具有高阻尼值的广域时滞PID阻尼控制器的设计方法，包括：确定电力***低频振荡模式，选取反应区间低频振荡特性的反馈控制信号以及合适的阻尼控制执行器，建立广域电力***的传递函数数学模型，设计适用于时滞PMU信号反馈控制的电力***PID控制器，计算能够确保电力***稳定的可能的PID参数分布范围，选取使电力***具有高阻尼值的PID参数。
本发明所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法的技术方案的实施步骤如下：
(1)确定电力***低频振荡模式，包括特征根、振荡频率、阻尼值、参与机组，从中筛选出区间低频振荡模式；
(2)针对区间低频振荡模式，分析同步PMU信号对该模式的可观性，从中筛选出区间低频振荡的广域反馈控制信号；分析电力***现有的调控装置，从中选取对该模式具有较高可控性的调控装置作为阻尼控制的执行器；
(3)确定电力***从步骤(2)中所选定的阻尼调控装置输入端至广域反馈控制信号的局部线性化传递函数模型；
(4)设计电力***时滞PID阻尼控制器的结构；
(5)根据广域反馈控制信号的时滞和电力***传递函数模型，计算能够确保电力***稳定运行的PID参数分布范围；
(6)将步骤(5)中的PID参数分布范围离散化，从中选取X组PID参数；
(7)将步骤(6)中的PID参数分别加入电力***中，计算与之对应的区间低频振荡阻尼值，从中选取阻尼值最大的那组参数作为时滞PID阻尼控制器的参数。
进一步，所述步骤(1)中，既可以采取电力***小干扰稳定特征根分析，也可以采取电力***动态仿真或者测量数据辨识技术，确定电力***的低频振荡模式的特征根、振荡频率、阻尼值和参与机组，其中振荡频率低于1.0Hz且功率振荡参与机组分布在多个区域电网中的振荡模式就是区间低频振荡模式。
进一步，所述步骤(2)中，可选的广域反馈信号来自各厂站、线路上布置的PMU装置，常见的有发电机功角和转速信号、联络线有功功率信号，通过比较这些信号对区间低频振荡模式的可观性指标，可以选出可观性较好的PMU信号作为反馈控制信号Y；可选的低频振荡稳定调控装置有发电机励磁装置、高压直流(HVDC)控制装置、静态无功补偿(SVC)装置等柔***流输电(FACTS)装置，根据这些装置附加控制输入对该振荡模式的可控性指标，选取可控性较好的调控装置作为低频振荡阻尼控制的执行器，执行器的输入信号为U；反馈信号至阻尼控制执行器的时滞设定为τ。
进一步，所述步骤(3)中，既可以采取电力***线性化数学模型、也可以采取辨识算法确定从阻尼调控装置输入信号U至广域反馈控制信号Y的线性化传递函数模型G(s)＝N(s)e-τs/D(s)，其中N(s)＝bmsm+bm-1sm- 1+…+b1s+b0；D(s)＝sn+an-1sn-1+…+a1s+a0；考虑到实际电力***模型的阶数可能非常高，可以进一步采取降阶措施，使降阶后的传递函数模型包含主要的低频振荡模式，传递函数模型的形式同降阶前G(s)＝N(s)e- τs/D(s)。
进一步，所述步骤(4)中，所设计的电力***广域时滞PID阻尼控制器主要包括：广域测量信号预处理模块、比例(P)环节、积分(I)环节、微分(D)环节、输出限幅环节；其中电力***WAMS中的PMU信号输入到该时滞PID阻尼控制器的测量信号预处理模块，剔除错误的数据和时滞过大的数据，并将其与稳态值进行比较，对不同时滞的PMU信号重新进行排队、根据设定的时滞等待后送入PID环节，PID环节输出的阻尼控制信号经过限幅环节后送入步骤(2)中所选定的低频振荡调控装置，作为附加控制信号参与电力***的稳定控制。
进一步，所述步骤(5)中，针对步骤(3)中获取的电力***传递函数模型G(s)以及步骤(4)中设定的广域PMU反馈控制信号的时滞τ，按以下几个步骤计算时滞PID阻尼控制器的参数分布范围(KP，KD，KI)：
(a)选取足够大的l，若n是偶数，则令Z＝2lπ，否则令Z＝2lπ+π/2；令s＝jz/τ，z为实数；假定Q是曲线f2(z)＝-q1(z)/[Nr 2(z)+Ni 2(z)]与直线f1(z)＝KP在区间(0，Z)内的交点数量，其中q1(z)＝[Dr(z)Nr(z)+Di(z)Ni(z)]cos(z)-[Di(z)Nr(z)-Dr(z)Ni(z)]sin(z)，Nr(z)、Ni(z)、Dr(z)、Di(z)分别为N(jz/τ)和D(jz/τ)的实部和虚部；确定KP的分布范围[KPmin，KPmax]，使Q满足下式：
其中，l(N)、r(N)和j(N)分别为N(s)在s左半平面、右半平面和正虚轴上的零点数量；
(b)将KP分布范围[KPmin，KPmax]等间隔分为F段，间隔点分别为KP0、KP1、KP2、……、KPi、KP(i+1)、……、KPF，其中i＝0、1、2、……、F，KP0＝KPmin，KPF＝KPmax；
(c)对于给定的KP＝KPi，其中i＝0、1、2、……、F，计算q(z，KP)＝{q1(z)+KP[Nr 2(z)+Ni 2(z)]}z/τ在区间[0，Z)内不同的实零点，从小到大依次为z0、z1、z2、……、zc-1，且zc＝Z；
(d)对于t＝1、2、……、c，如果N(-jzt/τ)＝0，则it＝0；若N(-s)在原点有个零点，则令i0＝sgn(d[p1(z)]/dz|z＝0)，其中sgn()为符号函数，p1(z)＝-{[Di(z)Nr(z)-Dr(z)Ni(z)]cos(z)+[Dr(z)Nr(z)+Di(z)Ni(z)]sin(z)}z/τ；否则it＝-1或1，具体由下式决定：
其中
假设满足上述条件的集合I有h组；
(e)若根据步骤(d)得到的I是唯一的，计算由不等式组[KI-A(zt)KD+B(zt)]it>0确定的(KDi，KIi)稳定区间的交集Si，其中A(zt)＝zt 2/τ2，B(zt)＝p1(zt)/[Nr 2(zt)+Ni 2(zt)]，t＝0、1、2、……、c，且t满足N(jzt/τ)≠0；若I不是唯一的，(KDi，KIi)则是步骤(d)中的h组I所对应的稳定区间的并集Si；
(f)返回步骤(c)，直至所有的KPi所对应的(KDi，KIi)稳定区间Si计算完毕。
(g)得到能够确保***G(s)稳定的时滞PID参数分布范围(KP，KD，KI)为(KPi，Si)，其中i＝0、1、2、……、F。
进一步，所述步骤(6)中，在KP、KD、KI三维坐标系中，由步骤(5)计算所得时滞PID参数分布范围(KPi，Si)形成一个立体空间，在该立体空间中等间隔或者随机选取X组PID参数。
进一步，所述步骤(7)中，将步骤(6)中所选的X组PID参数分别加入电力***，既可以采取电力***小干扰稳定特征根分析，也可以采取电力***动态仿真辨识技术，确定电力***的低频振荡模式的特征根、振荡频率和阻尼值；将区间低频振荡模式所对应阻尼值最大的那组PID参数作为所设计的具有高阻尼值的时滞PID控制器的最终参数。
本发明的优点是：不仅适用于广域反馈PMU信号存在时滞的情形、抑制区间低频振荡效果突出，而且所设计的控制器具有简洁的结构、易于工程实现。
附图说明：
图1具有高阻尼值的广域时滞PID阻尼控制器设计流程图
图2电力***广域时滞PID阻尼控制结构图
图3四机两区域电力***
图4时滞PID参数分布范围
图5装设时滞PID阻尼控制器前的单回区间联络线振荡功率曲线
图6装设时滞PID阻尼控制器后的单回区间联络线振荡功率曲线
具体实施方式：
下面结合实施例及附图，对本发明作进一步的详细说明，但本发明的实施方式不限于此。
本发明所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法的技术方案的流程图如图1所示，具体实施步骤如下：
(1)采取电力***小干扰稳定特征根分析方法，或者采取电力***动态仿真或者测量数据辨识技术，确定电力***的低频振荡模式的特征根、振荡频率、阻尼值和参与机组，其中振荡频率低于1.0Hz且功率振荡参与机组分布在多个区域电网中的振荡模式就是区间低频振荡模式。
(2)针对区间低频振荡模式，分析同步PMU信号对该模式的可观性，可选的广域反馈信号来自各厂站、线路上布置的PMU装置，常见的有发电机功角和转速信号、联络线有功功率信号，选出可观性较好的PMU信号作为反馈控制信号Y；选取对低频振荡模式可控性较好的调控装置作为低频振荡阻尼控制的执行器，可选的低频振荡稳定调控装置有发电机励磁装置、高压直流(HVDC)控制装置、静态无功补偿(SVC)装置等柔***流输电(FACTS)装置，执行器的输入信号为U；反馈信号至阻尼控制执行器的时滞设定为τ
(3)采取小干扰稳定分析方法，或者采取辨识算法确定从阻尼调控装置输入信号U至广域反馈控制信号Y的局部线性化传递函数模型G(s)＝N(s)e-τs/D(s)，其中N(s)＝bmsm+bm-1sm-1+…+b1s+b0；D(s)＝sn+an- 1sn-1+…+a1s+a0；考虑到实际电力***模型的阶数可能非常高，可以进一步采取降阶措施，使降阶后的传递函数模型包含主要的低频振荡模式，传递函数模型的形式同降阶前G(s)＝N(s)e-τs/D(s)。
(4)设计电力***时滞PID阻尼控制器，所设计的电力***广域时滞PID阻尼控制器结构如图2所示，主要包括：广域测量信号预处理模块、比例(P)环节、积分(I)环节、微分(D)环节、输出限幅环节；其中电力***WAMS中的PMU信号输入到该时滞PID阻尼控制器的测量信号预处理模块，剔除错误的数据和时滞过大的数据，并将其与稳态值进行比较，对不同时滞的PMU信号重新进行排队、根据设定的时滞等待后送入PID环节，PID环节输出的阻尼控制信号经过限幅环节后送入步骤(2)中所选定的低频振荡调控装置，作为附加控制信号参与电力***的稳定控制。
(5)针对步骤(3)中获取的电力***传递函数模型G(s)以及步骤(2)中设定的广域PMU反馈控制信号的时滞τ，时滞PID阻尼控制器的参数分布范围(KP，KD，KI)的计算步骤如下：
a)选取足够大的l，若n是偶数，则令Z＝2lπ，否则令Z＝2lπ+π/2；令s＝jz/τ，z为实数；假定Q是曲线f2(z)＝-q1(z)/[Nr 2(z)+Ni 2(z)]与直线f1(z)＝KP在区间(0，Z)内的交点数量，其中q1(z)＝[Dr(z)Nr(z)+Di(z)Ni(z)]cos(z)-[Di(z)Nr(z)-Dr(z)Ni(z)]sin(z)，Nr(z)、Ni(z)、Dr(z)、Di(z)分别为N(jz/τ)和D(jz/τ)的实部和虚部；确定KP的分布范围[KPmin，KPmax]，使Q满足下式：
其中，l(N)、r(N)和j(N)分别为N(s)在s左半平面、右半平面和正虚轴上的零点数量；
b)将KP分布范围[KPmin，KPmax]等间隔分为F段，间隔点分别为KP0、KP1、KP2、……、KPi、KP(i+1)、……、KPF，其中i＝0、1、2、……、F，KP0＝KPmin，KPF＝KPmax；
c)对于给定的KP＝KPi，其中i＝0、1、2、……、F，计算q(z，KP)＝{q1(z)+KP[Nr 2(z)+Ni 2(z)]}z/τ在区间[0，Z)内不同的实零点，从小到大依次为z0、z1、z2、……、zc-1，且zc＝Z；
d)对于t＝1、2、……、c，如果N(-jzt/τ)＝0，则it＝0；若N(-s)在原点有个零点，则令i0＝sgn(d[p1(z)]/dz|z＝0)，其中sgn()为符号函数，p1(z)＝-{[Di(z)Nr(z)-Dr(z)Ni(z)]cos(z)+[Dr(z)Nr(z)+Di(z)Ni(z)]sin(z)}z/τ；否则it＝-1或1，具体由下式决定：
其中
假设满足上述条件的集合I有h组；
e)若根据步骤(d)得到的I是唯一的，计算由不等式组[KI-A(zt)KD+B(zt)]it>0确定的(KDi，KIi)稳定区间的交集Si，其中A(zt)＝zt 2/τ2，B(zt)＝p1(zt)/[Nr 2(zt)+Ni 2(zt)]，t＝0、1、2、……、c，且t满足N(jzt/τ)≠0；若I不是唯一的，(KDi，KIi)则是步骤(d)中的h组I所对应的稳定区间的并集Si；
f)返回步骤(c)，直至所有的KPi所对应的(KDi，KIi)稳定区间Si计算完毕。
g)得到能够确保***G(s)稳定的时滞PID参数分布范围(KP，KD，KI)为(KPi，Si)，其中i＝0、1、2、……、F。
(6)在KP、KD、KI三维坐标系中，由步骤(5)计算所得时滞PID参数分布范围(KPi，Si)形成一个立体空间，在该立体空间中等间隔或者随机选取X组PID参数。
(7)将步骤(6)中所选的X组PID参数分别加入电力***，既可以采取电力***小干扰稳定特征根分析，也可以采取电力***动态仿真辨识技术，确定电力***的低频振荡模式的特征根、振荡频率和阻尼值；将区间低频振荡模式所对应阻尼值最大的那组PID参数作为所设计的具有高阻尼值的时滞PID控制器的最终参数。
实施案例：四机两区域电力***仿真分析
如图3所示的四机两区域电力***，发电机G1和G2位于区域1，发电机G3和G4位于区域2，两个区域通过长联络线互联，区域1通过长联络线向区域2传输有功功率。通过小干扰稳定特征根分析可知，该***存在3种低频振荡模式，详见表1，其中模式1是区间低频振荡模式，区域1中的G1和G2相对区域2中的G3和G4存在低频功率振荡，振荡频率为0.5389Hz，阻尼值为0.0099。由于阻尼值很小，因此该***一旦受到扰动便会长时间剧烈振荡，其中联络线功率振荡曲线如图5所示。
根据可观性和可控性原则，结合电力***的实际PMU配置和控制调节能力，选择发电机G2与发电机G4之间的相对角速度(ω2-ω4)作为反馈信号Y，ω1、ω2、ω3和ω4分别为发电机G1、G2、G3和G4的角速度，可由PMU测量；选择发电机G2的励磁***作为附加阻尼控制的执行器，时滞PID阻尼控制器的输出端连接G2的励磁控制器的附加控制输入端U；时滞τ设定为300ms。
采取小干扰稳定分析方法，确定从发电机G2励磁控制装置输入端U至广域反馈控制信号Y的线性化传递函数模型为28阶，降为8阶后的传递函数为G(s)＝N(s)e-τs/D(s)，其中
N(s)＝0.02243s4-5.295s3-4.187s2-99.31s-7.837
D(s)＝s5+21.88s4+87.95s3+1109s2+917.5s+9788
对于如图2所示的电力***时滞PID阻尼控制器，计算PID参数的分布范围如图4所示。从中随机选取4638组PID参数，分别将这些PID参数赋值给电力***中的时滞PID阻尼控制器，计算区间低频振荡模式的阻尼值，阻尼值较大的10组PID参数仿真结果列于表2中，阻尼值最大的那组PID参数为(-1.6362，0.0000，-2.6010)，因此所设计的时滞PID阻尼控制器的参数则为：KP＝-2.6010，KI＝0.0000，KD＝-1.6362。
假设图3所示四机两区域电力***0.1s时刻，联络线7-8一条回路在靠近母线7处发生三相接地短路，50ms后连接母线7的故障联络线端断开，再过50ms后连接母线8的故障联络线端断开，此时该***以单回路联络线互联继续运行。在此过程中，受三相接地短路故障的冲击影响，***将经历一个暂态过程。未安装阻尼控制器前，由于***存在阻尼值很小的低频振荡模式，因此如图5所示，区间联络线8-9的功率出现了幅值很大的振荡；安装广域时滞PID阻尼控制器后，由图6可见，区间联络线8-9的功率振荡得到了有效抑制，***的区间低频振荡阻尼值也由0.0099提高至0.0675。
表1低频振荡模式分析结果
模式编号 特征根 阻尼值 频率/Hz 参与机组
1 -0.0334±3.3862i 0.0099 0.5389 [G1、G2]V.S[G3、G4]
2 -0.7648±6.4155i 0.1184 1.0211 G3V.S G4
3 -0.7830±6.4652i 0.1202 1.0290 G1V.S G2
表2PID控制器阻尼效果分析结果

Claims (8)

1.具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法，包括如下步骤：
(1)确定电力***低频振荡模式，包括特征根、振荡频率、阻尼值、参与机组，从中筛选出区间低频振荡模式；
(2)针对区间低频振荡模式，分析同步PMU信号对该模式的可观性，从中筛选出区间低频振荡的广域反馈控制信号；分析电力***现有的调控装置，从中选取对该模式具有较高可控性的调控装置作为阻尼控制的执行器；
(3)确定电力***从步骤(2)中所选定的阻尼调控装置输入端至广域反馈控制信号的局部线性化传递函数模型；
(4)设计电力***时滞PID阻尼控制器的结构；
(5)根据广域反馈控制信号的时滞和电力***传递函数模型，计算能够确保电力***稳定运行的PID参数分布范围；
(6)将步骤(5)中的PID参数分布范围离散化，从中选取X组PID参数；
(7)将步骤(6)中的PID参数分别加入电力***中，计算与之对应的区间低频振荡阻尼值，从中选取阻尼值最大的那组参数作为时滞PID阻尼控制器的参数。
2.如权利要求1所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法，其特征在于：所述步骤(1)中，既可以采取电力***小干扰稳定特征根分析，也可以采取电力***动态仿真或者测量数据辨识技术，确定电力***的低频振荡模式的特征根、振荡频率、阻尼值和参与机组，其中振荡频率低于1.0Hz且功率振荡参与机组分布在多个区域电网中的振荡模式就是区间低频振荡模式。
3.如权利要求1所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法，其特征在于：所述步骤(2)中，可选的广域反馈信号来自各厂站、线路上布置的PMU装置，常见的有发电机功角和转速信号、联络线有功功率信号，通过比较这些信号对区间低频振荡模式的可观性指标，可以选出可观性较好的PMU信号作为反馈控制信号Y；可选的低频振荡稳定调控装置有发电机励磁装置、高压直流(HVDC)控制装置、静态无功补偿(SVC)装置等柔***流输电(FACTS)装置，根据这些装置附加控制输入对该振荡模式的可控性指标，选取可控性较好的调控装置作为低频振荡阻尼控制的执行器，执行器的输入信号为U；反馈信号至阻尼控制执行器的时滞设定为τ。
4.如权利要求1所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法，其特征在于：所述步骤(3)中，既可以采取电力***线性化数学模型、也可以采取辨识算法确定从阻尼调控装置输入信号U至广域反馈控制信号Y的线性化传递函数模型G(s)＝N(s)e-τs/D(s)，其中N(s)＝bmsm+bm-1sm-1+…+b1s+b0，D(s)＝sn+an-1sn-1+…+a1s+a0，s为复数，m和n分别为多项式N(s)和D(s)的阶数，b0、b1、……、bm为N(s)关于s多项式的系数，a0、a1、……、an-1为D(s)关于s多项式的系数；考虑到实际电力***模型的阶数可能非常高，可以进一步采取降阶措施，使降阶后的传递函数模型包含主要的低频振荡模式，传递函数模型的形式同降阶前G(s)＝N(s)e-τs/D(s)。
5.如权利要求1所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法，其特征在于：所述步骤(4)中，所设计的电力***广域时滞PID阻尼控制器主要包括：广域测量信号预处理模块、比例(P)环节、积分(I)环节、微分(D)环节、输出限幅环节；其中电力***WAMS中的PMU信号输入到该时滞PID阻尼控制器的测量信号预处理模块，剔除错误的数据和时滞过大的数据，并将其与稳态值进行比较，对不同时滞的PMU信号重新进行排队、根据设定的时滞等待后送入PID环节，PID环节输出的阻尼控制信号经过限幅环节后送入步骤(2)中所选定的低频振荡调控装置，作为附加控制信号参与电力***的稳定控制。
6.如权利要求1所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法，其特征在于：所述步骤(5)中，针对步骤(3)中获取的电力***传递函数模型G(s)以及步骤(4)中设定的广域PMU反馈控制信号的时滞τ，按以下几个步骤计算时滞PID阻尼控制器的参数分布范围(KP，KD，KI)，其中KP为比例环节系数，KI为积分环节系数，KD为微分环节系数：
(a)选取足够大的实数l，由其确定实数变量z的分布范围[0，Z)，其中Z为分布范围的上限；若n是偶数，则令Z＝2lπ，否则令Z＝2lπ+π/2；令s＝jz/τ，z为实数变量，j为虚数单位；假定Q是曲线f2(z)＝-q1(z)/[Nr 2(z)+Ni 2(z)]与直线f1(z)＝KP在区间(0，Z)内的交点数量，其中q1(z)＝[Dr(z)Nr(z)+Di(z)Ni(z)]cos(z)-[Di(z)Nr(z)-Dr(z)Ni(z)]sin(z)，Nr(z)、Ni(z)、Dr(z)、Di(z)分别为N(jz/τ)和D(jz/τ)的实部和虚部；确定KP的分布范围[KPmin，KPmax]，其中KPmin和KPmax分别为KP的下限和上限，使Q满足下式：
其中，l(N)、r(N)和j(N)分别为N(s)在s左半平面、右半平面和正虚轴上的零点数量；
(b)将KP分布范围[KPmin，KPmax]等间隔分为F段，间隔点分别为KP0、KP1、KP2、……、KPi、KP(i+1)、……、KPF，其中i＝0、1、2、……、F，KP0＝KPmin，KPF＝KPmax；
(c)对于给定的KP＝KPi，其中i＝0、1、2、……、F，计算q(z，KP)＝{q1(z)+KP[Nr 2(z)+Ni 2(z)]}z/τ在区间[0，Z)内不同的实零点，假设共有c个，从小到大依次为z0、z1、z2、……、zc-1，且zc＝Z；
(d)对于t＝1、2、……、c，如果N(-jzt/τ)＝0，则it＝0；若N(-s)在原点有个零点，则令i0＝sgn(d[p1(z)]/dz|z＝0)，其中sgn()为符号函数，p1(z)＝-{[Di(z)Nr(z)-Dr(z)Ni(z)]cos(z)+[Dr(z)Nr(z)+Di(z)Ni(z)]sin(z)}z/τ；否则it＝-1或1，具体由下式决定：
其中
假设满足上述条件的集合I有h组；
(e)若根据步骤(d)得到的I是唯一的，计算由不等式组[KI-A(zt)KD+B(zt)]it>0确定的(KDi，KIi)稳定区间的交集Si，其中A(zt)＝zt 2/τ2，B(zt)＝p1(zt)/[Nr 2(zt)+Ni 2(zt)]，t＝0、1、2、……、c，且t满足N(jzt/τ)≠0；若I不是唯一的，(KDi，KIi)则是步骤(d)中的h组I所对应的稳定区间的并集Si；
(f)返回步骤(c)，直至所有的KPi所对应的(KDi，KIi)稳定区间Si计算完毕。
(g)得到能够确保***G(s)稳定的时滞PID参数分布范围(KP，KD，KI)为(KPi，Si)，其中i＝0、1、2、……、F。
7.如权利要求1所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法，其特征在于：所述步骤(6)中，在KP、KD、KI三维坐标系中，由步骤(5)计算所得时滞PID参数分布范围(KPi，Si)形成一个立体空间，在该立体空间中等间隔或者随机选取X组PID参数。
8.如权利要求1所述的具有高阻尼值的电力***广域时滞PID阻尼控制器设计方法，其特征在于：所述步骤(7)中，将步骤(6)中所选的X组PID参数分别加入电力***，既可以采取电力***小干扰稳定特征根分析，也可以采取电力***动态仿真辨识技术，确定电力***的低频振荡模式的特征根、振荡频率和阻尼值；将区间低频振荡模式所对应阻尼值最大的那组PID参数作为所设计的具有高阻尼值的时滞PID控制器的最终参数。
CN201510492671.9A 2015-08-12 2015-08-12 具有高阻尼值的电力***广域时滞pid阻尼控制器设计方法 Active CN105140935B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510492671.9A CN105140935B (zh) 2015-08-12 2015-08-12 具有高阻尼值的电力***广域时滞pid阻尼控制器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510492671.9A CN105140935B (zh) 2015-08-12 2015-08-12 具有高阻尼值的电力***广域时滞pid阻尼控制器设计方法

Publications (2)

Publication Number Publication Date
CN105140935A true CN105140935A (zh) 2015-12-09
CN105140935B CN105140935B (zh) 2017-11-07

Family

ID=54726193

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510492671.9A Active CN105140935B (zh) 2015-08-12 2015-08-12 具有高阻尼值的电力***广域时滞pid阻尼控制器设计方法

Country Status (1)

Country Link
CN (1) CN105140935B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107065519A (zh) * 2016-12-08 2017-08-18 浙江工业大学 一种pmu反馈控制信号预处理方法
CN107306033A (zh) * 2016-04-20 2017-10-31 华北电力大学（保定） 一种计及广域信号多时滞的电力***附加阻尼控制方法
CN110365009A (zh) * 2019-06-22 2019-10-22 东北电力大学 一种适用于互联电力***的广域阻尼优化控制方法
CN112234628A (zh) * 2020-10-21 2021-01-15 国网黑龙江省电力有限公司电力科学研究院 一种电力***的低频振荡抑制方法
CN112985530A (zh) * 2021-02-01 2021-06-18 南京航空航天大学 一种基于特征方程根轨迹的燃油计量装置设计参数调整方法
CN113595102A (zh) * 2021-06-24 2021-11-02 国网浙江省电力有限公司嘉兴供电公司 一种基于储能电源阻尼电力***低频振荡的控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101615791A (zh) * 2009-07-29 2009-12-30 清华大学 有效抑制次同步谐振和振荡的次同步阻尼控制***
CN101917003A (zh) * 2010-07-29 2010-12-15 东南大学 一种电力***小干扰振荡稳定性分解式模态分析方法
CN102185326A (zh) * 2011-05-17 2011-09-14 华北电力大学 一种具有鲁棒性的广域阻尼控制***设计方法
CN102474100A (zh) * 2009-08-06 2012-05-23 Abb研究有限公司 电力传输***中的功率或电压振荡阻尼
US20150008672A1 (en) * 2012-01-27 2015-01-08 Vestas Wind Systems A/S Method for damping drive train oscillations in a wind turbine generator

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101615791A (zh) * 2009-07-29 2009-12-30 清华大学 有效抑制次同步谐振和振荡的次同步阻尼控制***
CN102474100A (zh) * 2009-08-06 2012-05-23 Abb研究有限公司 电力传输***中的功率或电压振荡阻尼
CN101917003A (zh) * 2010-07-29 2010-12-15 东南大学 一种电力***小干扰振荡稳定性分解式模态分析方法
CN102185326A (zh) * 2011-05-17 2011-09-14 华北电力大学 一种具有鲁棒性的广域阻尼控制***设计方法
US20150008672A1 (en) * 2012-01-27 2015-01-08 Vestas Wind Systems A/S Method for damping drive train oscillations in a wind turbine generator

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
戚军 等: "电力***传递函数的通用Prony辨识算法", 《中国电机工程学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107306033A (zh) * 2016-04-20 2017-10-31 华北电力大学（保定） 一种计及广域信号多时滞的电力***附加阻尼控制方法
CN107306033B (zh) * 2016-04-20 2020-07-03 华北电力大学（保定） 一种计及广域信号多时滞的电力***附加阻尼控制方法
CN107065519A (zh) * 2016-12-08 2017-08-18 浙江工业大学 一种pmu反馈控制信号预处理方法
CN107065519B (zh) * 2016-12-08 2020-01-14 浙江工业大学 一种pmu反馈控制信号预处理方法
CN110365009A (zh) * 2019-06-22 2019-10-22 东北电力大学 一种适用于互联电力***的广域阻尼优化控制方法
CN112234628A (zh) * 2020-10-21 2021-01-15 国网黑龙江省电力有限公司电力科学研究院 一种电力***的低频振荡抑制方法
CN112234628B (zh) * 2020-10-21 2022-10-11 国网黑龙江省电力有限公司电力科学研究院 一种电力***的低频振荡抑制方法
CN112985530A (zh) * 2021-02-01 2021-06-18 南京航空航天大学 一种基于特征方程根轨迹的燃油计量装置设计参数调整方法
CN112985530B (zh) * 2021-02-01 2022-04-22 南京航空航天大学 一种基于特征方程根轨迹的燃油计量装置设计参数调整方法
CN113595102A (zh) * 2021-06-24 2021-11-02 国网浙江省电力有限公司嘉兴供电公司 一种基于储能电源阻尼电力***低频振荡的控制方法

Also Published As

Publication number Publication date
CN105140935B (zh) 2017-11-07

Similar Documents

Publication Publication Date Title
CN105140935A (zh) 具有高阻尼值的电力***广域时滞pid阻尼控制器设计方法
CN105071414A (zh) 电力***广域时滞pid阻尼控制器设计方法
Jan et al. Adaptive and fuzzy PI controllers design for frequency regulation of isolated microgrid integrated with electric vehicles
CN105207235B (zh) 适用于大范围变化时滞的广域pid阻尼控制器及其设计方法
CN105117535A (zh) 适用于随机时滞的电力***广域pid阻尼控制器设计方法
CN107482649B (zh) 一种基于分频控制的两域互联***负荷频率控制方法
CN103311939A (zh) 基于wams的电力***低频振荡协调阻尼控制方法
CN108288858B (zh) 主被动频率响应切换控制方法
Eriksson On the centralized nonlinear control of HVDC systems using Lyapunov theory
CN105186511B (zh) 电池储能***参与电网二次调频控制方法
Saadatmand et al. PMU-based FOPID controller of large-scale wind-PV farms for LFO damping in smart grid
Piao et al. Control strategy of battery energy storage system to participate in the second frequency regulation
CN107346889B (zh) 考虑一二次调频及最小频率偏差的负荷削减优化模型构建方法
CN103887800B (zh) 基于能量裕度约束的电压在线协调控制方法
CN109546677A (zh) 一种规模化海上风电柔性送出***安全控制策略求解方法
CN109449969A (zh) 一种用于孤岛模式下微电网稳定性分析的建模方法
CN108110776B (zh) 主被动频率响应复合控制方法
CN106532785A (zh) 一种计及新型云模型的负荷频率控制方法
CN105470976B (zh) 一种稳态条件下svc和tcsc的协调配置方法
Retnamony et al. Congestion management is to enhance the transient stability in a deregulated power system using FACTS devices
Shahgholian et al. Transient stability promotion by FACTS controller based on adaptive inertia weight particle swarm optimization method
CN103956767A (zh) 一种考虑尾流效应的风电场并网稳定性分析方法
Shahgholian et al. Design of a new IPFC-based damping neurocontrol for enhancing stability of a power system using particle swarm optimization
CN110086192A (zh) 两区域柔性直流互联***频率一致性控制方法
Hu et al. Robust Cooperative Load Frequency Control for Enhancing Wind Energy Integration in Multi-Area Power Systems

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant