CN110134908A - 一种电力系统稳定器pss增益计算方法以及计算装置 - Google Patents
一种电力系统稳定器pss增益计算方法以及计算装置 Download PDFInfo
- Publication number
- CN110134908A CN110134908A CN201910418603.6A CN201910418603A CN110134908A CN 110134908 A CN110134908 A CN 110134908A CN 201910418603 A CN201910418603 A CN 201910418603A CN 110134908 A CN110134908 A CN 110134908A
- Authority
- CN
- China
- Prior art keywords
- pss
- unit
- coefficient
- generator
- modulus value
- 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
Links
Classifications
-
- 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本申请提供一种电力系统稳定器PSS增益计算方法以及计算装置,该方法包括:计算机组参与强相关振荡模式的角频率ωd;获取机组的PSS的传递函数在强相关振荡模式下的第一模值获取机组的励磁系统的传递函数在强相关振荡模式下的第二模值线性化机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2;根据第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算机组的本机振荡角频率ωn;根据第一模值、第二模值、第二系数K2和本机振荡角频率ωn,计算PSS的增益。这样,可以根据第一模值、第二模值、第二系数K2和本机振荡角频率ωn,计算PSS的增益。节省时间,效率较高,且操作简单。
Description
技术领域
本申请涉及电力系统控制技术领域,尤其涉及一种电力系统稳定器PSS增益计算方法以及计算装置。
背景技术
系统如果发生低频振荡,会导致电网无法安全稳定运行。目前主要采用电力系统稳定器(Power System Stabilizer,PSS)来提高动态阻尼,抑制低频振荡发生。PSS是励磁系统的一个附加装置,它通过提取与振荡有关的信号在PSS、励磁系统以及发电机电磁回路上产生一个附加阻尼转矩,来抑制低频振荡的发生。
相关技术中,PSS增益的整定方法,是通过临界增益试验来整定的:逐渐增加PSS增益,观察机组电压阶跃试验的稳定性,直到出现不稳定现象为止。此时对应的增益就是PSS的临界增益。该机组PSS增益的整定值,整定为PSS的临界增益的1/3至1/2之间的一个值。
通过试验整定方法确定PSS增益的整定值,试验过程较长,效率较低。
发明内容
本申请提供了一种电力系统稳定器PSS增益计算方法以及计算装置,以解决相关技术中,通过试验整定方法确定PSS增益的整定值,试验过程较长,效率较低的问题。
一方面,本申请提供一种电力系统稳定器PSS增益计算方法,包括:
计算机组参与强相关振荡模式的角频率ωd;
获取所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值
获取所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值
线性化所述机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2;
根据所述第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算所述机组的本机振荡角频率ωn;
根据所述第一模值、所述第二模值、所述第二系数K2和所述本机振荡角频率ωn,计算所述PSS的增益;
其中,所述PSS的增益通过以下公式表示:
所述KPSS为所述PSS的增益,所述ε为阻尼比,D为所述发电机的转子的机械阻尼转矩系数,De为无PSS时发电机的电气阻尼转矩系数。
可选的,所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值通过以下公式表示:
其中,T1、T2、T3、T4、T5和T6为所述PSS的补偿环节参数。
可选的,在PSS信号接入所述励磁系统的第一位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
在PSS信号接入所述励磁系统的第二位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
其中,K、KA、KF、KV、TA、TF、T1、T2、T3和T4为所述励磁系统的参数。
可选的,所述第一系数K1通过以下公式表示:
K1=Fd(xq-x′d)Iq0+Fq[E′q0+(xq-x′d)Id0]
所述第二系数K2通过以下公式表示:
K2=Iq0+Yd(xq-x′d)Iq0+Yq[E′q0+(xq-x′d)Id0]
Z2=R2+(X+xq)(X+x′d)
其中,xq为所述发电机的q轴同步电抗,x′d为所述发电机的d轴暂态电抗,Iq0为所述发电机的定子电流q轴分量的初值,Id0为所述发电机的定子电流d轴分量的初值,E′q0为所述发电机的暂态电势的初值,U为电网电压,R为电厂至所述电网的电阻,X为所述电厂至所述电网的电抗,δ0为所述发电机的功角的初值。
可选的,所述机组的本机振荡角频率ωn通过以下公式表示:
另一方面,本申请还提供一种电力系统稳定器PSS增益计算装置,包括:
第一计算模块,用于计算机组参与强相关振荡模式的角频率ωd;
第一获取模块,用于获取所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值
第二获取模块,用于获取所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值
线性化模块,用于线性化所述机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2;
第二计算模块,用于根据所述第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算所述机组的本机振荡角频率ωn;
第三计算模块,用于根据所述第一模值、所述第二模值、所述第二系数K2和所述本机振荡角频率ωn,计算所述PSS的增益;
其中,所述PSS的增益通过以下公式表示:
所述KPSS为所述PSS的增益,所述ε为阻尼比,D为所述发电机的转子的机械阻尼转矩系数,De为无PSS时发电机的电气阻尼转矩系数。
可选的,所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值通过以下公式表示:
其中,T1、T2、T3、T4、T5和T6为所述PSS的补偿环节参数。
可选的,在PSS信号接入所述励磁系统的第一位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
在PSS信号接入所述励磁系统的第二位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
其中,K、KA、KF、KV、TA、TF、T1、T2、T3和T4为所述励磁系统的参数。
可选的,所述第一系数K1通过以下公式表示:
K1=Fd(xq-x′d)Iq0+Fq[E′q0+(xq-x′d)Id0]
所述第二系数K2通过以下公式表示:
K2=Iq0+Yd(xq-x′d)Iq0+Yq[E′q0+(xq-x′d)Id0]
Z2=R2+(X+xq)(X+x′d)
其中,xq为所述发电机的q轴同步电抗,x′d为所述发电机的d轴暂态电抗,Iq0为所述发电机的定子电流q轴分量的初值,Id0为所述发电机的定子电流d轴分量的初值,E′q0为所述发电机的暂态电势的初值,U为电网电压,R为电厂至所述电网的电阻,X为所述电厂至所述电网的电抗,δ0为所述发电机的功角的初值。
可选的,所述机组的本机振荡角频率ωn通过以下公式表示:
由以上技术方案可知,本申请提供一种电力系统稳定器PSS增益计算方法以及计算装置,所述方法包括:计算机组参与强相关振荡模式的角频率ωd;获取所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值获取所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值线性化所述机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2;根据所述第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算所述机组的本机振荡角频率ωn;根据所述第一模值、所述第二模值、所述第二系数K2和所述本机振荡角频率ωn,计算所述PSS的增益;其中,所述PSS的增益通过以下公式表示:
所述KPSS为所述PSS的增益,所述ε为阻尼比,D为所述发电机的转子的机械阻尼转矩系数,De为无PSS时发电机的电气阻尼转矩系数。这样,可以根据第一模值、第二模值、第二系数K2和本机振荡角频率ωn,计算PSS的增益。节省时间,效率较高,且操作简单。可以降低试验对机组的不利影响,减少试验人员的安全风险。
附图说明
为了更清楚地说明本申请的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请提供的一种电力系统稳定器PSS增益计算方法的流程图;
图2为本申请提供的PSS2B模型的示意图;
图3为本申请提供的励磁系统模型的示意图;
图4为本申请提供的一种计算电力系统稳定器PSS增益的示例的流程图;
图5为本申请提供的励磁系统参数表;
图6为本申请提供的PSS参数表;
图7为本申请提供的一种电力系统稳定器PSS增益计算装置的结构图。
具体实施方式
下面将详细地对实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下实施例中描述的实施方式并不代表与本申请相一致的所有实施方式。仅是与权利要求书中所详述的、本申请的一些方面相一致的系统和方法的示例。
参见图1,图1是本申请提供的一种电力系统稳定器PSS增益计算方法的流程图。
如图1所示,包括以下步骤:
步骤101、计算机组参与强相关振荡模式的角频率ωd。
在步骤101中,可以计算机组参与强相关振荡模式的角频率ωd。例如,可以利用电力系统分析软件的小扰动功能,计算机组参与强相关振荡模式的频率fd,则角频率为ωd=2πfd。电力系统分析软件可以为电力系统分析综合程序(Power System AnalysisSynthesis Program,PSASP)。
步骤102、获取所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值
在步骤102中,可以获取机组的PSS的传递函数在强相关振荡模式下的第一模值GPSS(s)为PSS的传递函数,表示求其复数的模。
可选的,所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值通过以下公式表示:
其中,T1、T2、T3、T4、T5和T6为所述PSS的补偿环节参数。
如图2所示,为PSS2B模型的示意图。在图2中,示出了测量环节、隔离环节、比较环节、滤波环节、增益、补偿环节和限幅环节。机组的PSS的传递函数在强相关振荡模式下的第一模值可以通过以下公式表示:
其中,T1、T2、T3、T4、T5和T6为PSS的补偿环节参数。例如,假设步骤101中计算出机组参与强相关振荡模式的频率fd为0.5Hz,则角频率ωd=2πfd=π。
s=jωd=j2πfd=jπ
在PSS的补偿环节参数T1、T2、T3、T4、T5和T6已知的条件下,就可以依照上述公式计算出机组的PSS的传递函数在强相关振荡模式下的第一模值
步骤103、获取所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值
在步骤103中,可以获取机组的励磁系统的传递函数在强相关振荡模式下的第二模值GE(s)为励磁系统的传递函数,表示求其复数的模。
可选的,在PSS信号接入所述励磁系统的第一位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
在PSS信号接入所述励磁系统的第二位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
其中,K、KA、KF、KV、TA、TF、T1、T2、T3和T4为所述励磁系统的参数。
如图3所示,为励磁系统模型的示意图。需要说明的是,对于实际机组的PSS信号接入励磁系统的位置是确定的。例如,在图3所示的励磁系统中,PSS信号的接入位置可以为位置1,也可以为位置2。
在PSS信号接入励磁系统的第一位置的情况下,也即在PSS信号接入励磁系统的位置1的情况下,机组的励磁系统的传递函数在强相关振荡模式下的第二模值通过以下公式表示:
在PSS信号接入励磁系统的第二位置的情况下,也即在PSS信号接入励磁系统的位置2的情况下,机组的励磁系统的传递函数在强相关振荡模式下的第二模值通过以下公式表示:
其中,K、KA、KF、KV、TA、TF、T1、T2、T3和T4为励磁系统的参数。这样,在励磁系统的参数K、KA、KF、KV、TA、TF、T1、T2、T3和T4已知的条件下,就可以依照上述公式计算出机组的励磁系统的传递函数在强相关振荡模式下的第二模值
步骤104、线性化所述机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2。
在步骤104中,可以线性化机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2。
可选的,所述第一系数K1通过以下公式表示:
K1=Fd(xq-x′d)Iq0+Fq[E′q0+(xq-x′d)Id0]
所述第二系数K2通过以下公式表示:
K2=Iq0+Yd(xq-x′d)Iq0+Yq[E′q0+(xq-x′d)Id0]
Z2=R2+(X+xq)(X+x′d)
其中,xq为所述发电机的q轴同步电抗,x′d为所述发电机的d轴暂态电抗,Iq0为所述发电机的定子电流q轴分量的初值,Id0为所述发电机的定子电流d轴分量的初值,E′q0为所述发电机的暂态电势的初值,U为电网电压,R为电厂至所述电网的电阻,X为所述电厂至所述电网的电抗,δ0为所述发电机的功角的初值。
第一系数K1可以通过以下公式表示:
K1=Fd(xq-x′d)Iq0+Fq[E′q0+(xq-x′d)Id0]
第二系数K2可以通过以下公式表示:
K2=Iq0+Yd(xq-x′d)Iq0+Yq[E′q0+(xq-x′d)Id0]
Z2=R2+(X+xq)(x+x′d)
其中,xq为发电机的q轴同步电抗,x′d为发电机的d轴暂态电抗,Iq0为发电机的定子电流q轴分量的初值,Id0为发电机的定子电流d轴分量的初值,E′q0为发电机的暂态电势的初值,U为电网电压,R为电厂至电网的电阻,X为电厂至电网的电抗,δ0为发电机的功角的初值。
步骤105、根据所述第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算所述机组的本机振荡角频率ωn。
在步骤105中,可以根据第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算机组的本机振荡角频率ωn。其中,ω0=2πf0。f0一般为50赫兹。
可选的,所述机组的本机振荡角频率ωn通过以下公式表示:
机组的本机振荡角频率ωn可以通过以下公式表示:
步骤106、根据所述第一模值、所述第二模值、所述第二系数K2和所述本机振荡角频率ωn,计算所述PSS的增益;
其中,所述PSS的增益通过以下公式表示:
所述KPSS为所述PSS的增益,所述ε为阻尼比,D为所述发电机的转子的机械阻尼转矩系数,De为无PSS时发电机的电气阻尼转矩系数。
在步骤106中,可以根据第一模值第二模值第二系数K2和本机振荡角频率ωn,计算PSS的增益。
其中,PSS的增益可以通过以下公式表示:
KPSS为PSS的增益,ε为有PSS时提供给系统的阻尼比,要求ε=0.1~0.3,一般取ε=0.15即可。D为发电机的转子的机械阻尼转矩系数,水轮机取D=0.5~1,汽轮机取D≈0。De为无PSS时发电机的电气阻尼转矩系数,De≈0。
本申请提供的电力系统稳定器PSS增益计算方法,计算机组参与强相关振荡模式的角频率ωd;获取所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值获取所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值线性化所述机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2;根据所述第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算所述机组的本机振荡角频率ωn;根据所述第一模值、所述第二模值、所述第二系数K2和所述本机振荡角频率ωn,计算所述PSS的增益。这样,可以根据第一模值、第二模值、第二系数K2和本机振荡角频率ωn,计算PSS的增益。节省时间,效率较高,且操作简单。可以降低试验对机组的不利影响,减少试验人员的安全风险。
参见图4,图4是本申请提供的一种计算电力系统稳定器PSS增益的示例的流程图。
如图4所示,包括以下步骤:
步骤401、利用电力系统分析软件的小扰动功能,计算机组参与强相关振荡模式的频率fd,则角频率为ωd=2πfd。假设机组参与强相关振荡模式的频率fd=0.64Hz,则机组参与强相关振荡模式的角频率ωd=2πfd=4.021弧度/秒。
1、发电机参数:发电机额定功率为700MW,额定功率因数为0.9,发电机d轴同步电抗xd=1.015pu,发电机q轴同步电抗xq=0.73pu,发电机d轴暂态电抗x′d=0.305pu,机组轴系的惯性时间常数TJ=8.72秒。
2、发电厂升压变压器(等效双绕组):阻抗电压12.8%,额定容量为800MVA。
3、送出线路:电阻=0,电抗=0.02461pu(100MVA基准容量)。
如图5所示,为励磁系统参数表。即图5为图3所示的励磁系统模型的参数表。
如图6所示,为PSS参数表。即图6为图2所示的PSS2B模型的参数表。需要说明的是,在图6的状态下,PSS信号接入励磁系统的第一位置,即PSS信号接入励磁系统的位置1。
步骤402、在强相关振荡模式下,计算机组的PSS的传递函数在强相关振荡模式下的第一模值计算机组的励磁系统的传递函数在强相关振荡模式下的第二模值
步骤403、利用第一系数K1和第二系数K2的表达式,计算得到第一系数K1=1.3056,第二系数K2=1.3377。
步骤404、根据本机振荡角频率ωn的计算公式计算得到ωn=6.858弧度/秒。
步骤405、将上述各量代入PSS的增益的计算公式取D=0.5、De=0、ε=0.15。计算得到PSS的增益KPSS=20.7。且ε=0.1时,计算得到KPSS=13.6;ε=0.3时,计算得到KPSS=42.1。
PSS的增益的计算值所属的区间为13.6~42.1,中间值为20.7。由图6所示的PSS参数表可知,通过试验整定方法确定PSS的增益的整定值为14。虽然满足阻尼比的要求,但PSS的增益的整定值偏低,造成附加阻尼不足的缺陷。
而采用本申请的电力系统稳定器PSS增益计算方法,可以根据第一模值、第二模值、第二系数K2和本机振荡角频率ωn,计算PSS的增益。节省时间,效率较高,且操作简单。可以降低试验对机组的不利影响,减少试验人员的安全风险,避免了PSS的增益的整定值可能偏小的问题。
参见图7,图7是本申请提供的一种电力系统稳定器PSS增益计算装置的结构图。如图7所示,电力系统稳定器PSS增益计算装置700包括第一计算模块701、第一获取模块702、第二获取模块703、线性化模块704、第二计算模块705和第三计算模块706,其中:
第一计算模块701,用于计算机组参与强相关振荡模式的角频率ωd;
第一获取模块702,用于获取所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值
第二获取模块703,用于获取所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值
线性化模块704,用于线性化所述机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2;
第二计算模块705,用于根据所述第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算所述机组的本机振荡角频率ωn;
第三计算模块706,用于根据所述第一模值、所述第二模值、所述第二系数K2和所述本机振荡角频率ωn,计算所述PSS的增益;
其中,所述PSS的增益通过以下公式表示:
所述KPSS为所述PSS的增益,所述ε为阻尼比,D为所述发电机的转子的机械阻尼转矩系数,De为无PSS时发电机的电气阻尼转矩系数。
可选的,所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值通过以下公式表示:
其中,T1、T2、T3、T4、T5和T6为所述PSS的补偿环节参数。
可选的,在PSS信号接入所述励磁系统的第一位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
在PSS信号接入所述励磁系统的第二位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
其中,K、KA、KF、KV、TA、TF、T1、T2、T3和T4为所述励磁系统的参数。
可选的,所述第一系数K1通过以下公式表示:
K1=Fd(xq-x′d)Iq0+Fq[E′q0+(xq-x′d)Id0]
所述第二系数K2通过以下公式表示:
K2=Iq0+Yd(xq-x′d)Iq0+Yq[E′q0+(xq-x′d)Id0]
Z2=R2+(X+xq)(X+x′d)
其中,xq为所述发电机的q轴同步电抗,x′d为所述发电机的d轴暂态电抗,Iq0为所述发电机的定子电流q轴分量的初值,Iq0为所述发电机的定子电流d轴分量的初值,E′q0为所述发电机的暂态电势的初值,U为电网电压,R为电厂至所述电网的电阻,X为所述电厂至所述电网的电抗,δ0为所述发电机的功角的初值。
可选的,所述机组的本机振荡角频率ωn通过以下公式表示:
电力系统稳定器PSS增益计算装置700能够实现图1和图4的方法实施例中电力系统稳定器PSS增益计算装置实现的各个过程,为避免重复,这里不再赘述。且电力系统稳定器PSS增益计算装置700可以实现根据第一模值、第二模值、第二系数K2和本机振荡角频率ωn,计算PSS的增益。节省时间,效率较高,且操作简单。可以降低试验对机组的不利影响,减少试验人员的安全风险,避免了PSS的增益的整定值可能偏小的问题。
本申请提供的实施例之间的相似部分相互参见即可,以上提供的具体实施方式只是本申请总的构思下的几个示例,并不构成本申请保护范围的限定。对于本领域的技术人员而言,在不付出创造性劳动的前提下依据本申请方案所扩展出的任何其他实施方式都属于本申请的保护范围。
Claims (10)
1.一种电力系统稳定器PSS增益计算方法,其特征在于,包括:
计算机组参与强相关振荡模式的角频率ωd;
获取所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值
获取所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值
线性化所述机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2;
根据所述第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算所述机组的本机振荡角频率ωn;
根据所述第一模值、所述第二模值、所述第二系数K2和所述本机振荡角频率ωn,计算所述PSS的增益;
其中,所述PSS的增益通过以下公式表示:
所述KpSS为所述PSS的增益,所述ε为阻尼比,D为所述发电机的转子的机械阻尼转矩系数,De为无PSS时发电机的电气阻尼转矩系数。
2.如权利要求1所述的方法,其特征在于,所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值通过以下公式表示:
其中,T1、T2、T3、T4、T5和T6为所述PSS的补偿环节参数。
3.如权利要求1或2所述的方法,其特征在于,在PSS信号接入所述励磁系统的第一位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
在PSS信号接入所述励磁系统的第二位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
其中,K、KA、KF、KV、TA、TF、T1、T2、T3和T4为所述励磁系统的参数。
4.如权利要求3所述的方法,其特征在于,所述第一系数K1通过以下公式表示:
K1=Fd(xq-x′d)Iq0+Fq[E′q0+(xq-x′d)Id0]
所述第二系数K2通过以下公式表示:
K2=Iq0+Yd(xq-x′d)Iq0+Yq[E′q0+(xq-x′d)Id0]
Z2=R2+(X+xq)(X+x′d)
其中,xq为所述发电机的q轴同步电抗,x′d为所述发电机的d轴暂态电抗,Iq0为所述发电机的定子电流q轴分量的初值,Id0为所述发电机的定子电流d轴分量的初值,E′q0为所述发电机的暂态电势的初值,U为电网电压,R为电厂至所述电网的电阻,X为所述电厂至所述电网的电抗,δ0为所述发电机的功角的初值。
5.如权利要求4所述的方法,其特征在于,所述机组的本机振荡角频率ωn通过以下公式表示:
6.一种电力系统稳定器PSS增益计算装置,其特征在于,包括:
第一计算模块,用于计算机组参与强相关振荡模式的角频率ωd;
第一获取模块,用于获取所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值
第二获取模块,用于获取所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值
线性化模块,用于线性化所述机组的发电机对应的电磁转矩方程,得到电磁转矩增量表达式的第一系数K1和第二系数K2;
第二计算模块,用于根据所述第一系数K1、工频角速度ω0和机组轴系的惯性时间常数TJ计算所述机组的本机振荡角频率ωn;
第三计算模块,用于根据所述第一模值、所述第二模值、所述第二系数K2和所述本机振荡角频率ωn,计算所述PSS的增益;
其中,所述PSS的增益通过以下公式表示:
所述KpSS为所述PSS的增益,所述ε为阻尼比,D为所述发电机的转子的机械阻尼转矩系数,De为无PSS时发电机的电气阻尼转矩系数。
7.如权利要求6所述的电力系统稳定器PSS增益计算装置,其特征在于,所述机组的PSS的传递函数在所述强相关振荡模式下的第一模值通过以下公式表示:
其中,T1、T2、T3、T4、T5和T6为所述PSS的补偿环节参数。
8.如权利要求6或7所述的电力系统稳定器PSS增益计算装置,其特征在于,在PSS信号接入所述励磁系统的第一位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
在PSS信号接入所述励磁系统的第二位置的情况下,所述机组的励磁系统的传递函数在所述强相关振荡模式下的第二模值通过以下公式表示:
其中,K、KA、KF、KV、TA、TF、T1、T2、T3和T4为所述励磁系统的参数。
9.如权利要求8所述的电力系统稳定器PSS增益计算装置,其特征在于,所述第一系数K1通过以下公式表示:
K1=Fd(xq-x′d)Iq0+Fq[E′q0+(xq-x′d)Id0]
所述第二系数K2通过以下公式表示:
K2=Iq0+Yd(xq-x′d)Iq0+Yq[E′q0+(xq-x′d)Id0]
Z2=R2+(X+xq)(X+x′d)
其中,xq为所述发电机的q轴同步电抗,x′d为所述发电机的d轴暂态电抗,Iq0为所述发电机的定子电流q轴分量的初值,Id0为所述发电机的定子电流d轴分量的初值,E′q0为所述发电机的暂态电势的初值,U为电网电压,R为电厂至所述电网的电阻,X为所述电厂至所述电网的电抗,δ0为所述发电机的功角的初值。
10.如权利要求9所述的电力系统稳定器PSS增益计算装置,其特征在于,所述机组的本机振荡角频率ωn通过以下公式表示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418603.6A CN110134908B (zh) | 2019-05-20 | 2019-05-20 | 一种电力系统稳定器pss增益计算方法以及计算装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910418603.6A CN110134908B (zh) | 2019-05-20 | 2019-05-20 | 一种电力系统稳定器pss增益计算方法以及计算装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110134908A true CN110134908A (zh) | 2019-08-16 |
CN110134908B CN110134908B (zh) | 2023-04-11 |
Family
ID=67571576
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910418603.6A Active CN110134908B (zh) | 2019-05-20 | 2019-05-20 | 一种电力系统稳定器pss增益计算方法以及计算装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110134908B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103795069A (zh) * | 2013-12-19 | 2014-05-14 | 国网上海市电力公司 | 基于单机无穷大系统的全网pss协调整定方法 |
CN103812122A (zh) * | 2013-11-19 | 2014-05-21 | 云南电网公司电网规划研究中心 | 基于约束变尺度pss参数优化整定方法 |
CN105826937A (zh) * | 2016-05-23 | 2016-08-03 | 国网浙江省电力公司电力科学研究院 | 一种超前相位自适应型电力系统稳定器 |
CN106026138A (zh) * | 2016-06-22 | 2016-10-12 | 国网江苏省电力公司 | 一种即插即用式的电力系统稳定器设计算法 |
US20170331289A1 (en) * | 2015-06-19 | 2017-11-16 | Zhejiang University | A parameter tuning approach for bypass damping filter to suppress subsynchronous resonance in power systems |
CN107732942A (zh) * | 2017-11-10 | 2018-02-23 | 云南电力试验研究院(集团)有限公司 | 电力系统稳定器整定参数的筛选方法 |
-
2019
- 2019-05-20 CN CN201910418603.6A patent/CN110134908B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103812122A (zh) * | 2013-11-19 | 2014-05-21 | 云南电网公司电网规划研究中心 | 基于约束变尺度pss参数优化整定方法 |
CN103795069A (zh) * | 2013-12-19 | 2014-05-14 | 国网上海市电力公司 | 基于单机无穷大系统的全网pss协调整定方法 |
US20170331289A1 (en) * | 2015-06-19 | 2017-11-16 | Zhejiang University | A parameter tuning approach for bypass damping filter to suppress subsynchronous resonance in power systems |
CN105826937A (zh) * | 2016-05-23 | 2016-08-03 | 国网浙江省电力公司电力科学研究院 | 一种超前相位自适应型电力系统稳定器 |
CN106026138A (zh) * | 2016-06-22 | 2016-10-12 | 国网江苏省电力公司 | 一种即插即用式的电力系统稳定器设计算法 |
CN107732942A (zh) * | 2017-11-10 | 2018-02-23 | 云南电力试验研究院(集团)有限公司 | 电力系统稳定器整定参数的筛选方法 |
Non-Patent Citations (2)
Title |
---|
余骏: "励磁系统的电力系统稳定器研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技II辑》 * |
邓若梅等: "采用PSS-AVR并联模型对励磁系统进行相位补偿的理论与实践", 《电气开关》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110134908B (zh) | 2023-04-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fan et al. | On active/reactive power modulation of DFIG-based wind generation for interarea oscillation damping | |
Wang et al. | Virtual synchronous control for grid-connected DFIG-based wind turbines | |
Fan et al. | Mitigating SSR using DFIG-based wind generation | |
Suriyaarachchi et al. | A procedure to study sub-synchronous interactions in wind integrated power systems | |
Wei et al. | Synchronous motor-generator pair to enhance small signal and transient stability of power system with high penetration of renewable energy | |
CN101834446B (zh) | 一种基于可控串补的次同步振荡抑制方法 | |
CN108631338A (zh) | 一种用于抑制双馈风电场并网次同步振荡的方法 | |
CN109120001A (zh) | 基于虚拟电阻的双馈风电场并网系统次同步振荡抑制方法 | |
CN108270240A (zh) | 一种海上风电场的次同步源-网联合阻尼抑制方法 | |
CN104329220B (zh) | 用于抑制风电机组扭振的扭转载荷控制器及控制方法 | |
CN101119095A (zh) | 基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统 | |
CN105826937B (zh) | 一种超前相位自适应型电力系统稳定器 | |
Gao et al. | Impact of DFIG‐based wind farm integration on sub‐synchronous torsional interaction between HVDC and thermal generators | |
CN109038649A (zh) | 一种基于虚拟惯量控制的dfig轴系振荡的控制方法 | |
CN108321841A (zh) | 宽频转子附加阻尼的双馈风电机组次同步振荡抑制方法 | |
CN109217365A (zh) | 一种无刷双馈电机虚拟同步控制方法 | |
CN113779697A (zh) | 一种汽轮发电机组扭振联合仿真方法、装置及存储介质 | |
CN106374524A (zh) | 一种upfc对异步风力发电机系统次同步谐振的抑制方法 | |
Tabesh et al. | Small-signal model and dynamic analysis of variable speed induction machine wind farms | |
Limebeer et al. | Research Paper No 8: Subsynchronous resonance of the Koeberg turbo-generators and of a laboratory micro-alternator system | |
Tsourakis et al. | A power system stabilizer for variable-speed wind generators | |
Singh et al. | Using generic wind turbine models to compare inertial response of wind turbine technologies | |
CN111740449A (zh) | 基于附加励磁信号注入法的风电机群对同步发电机组阻尼特性分析方法 | |
CN109755968A (zh) | 一种双馈风电机组的神经网络保性能虚拟同步控制方法 | |
CN107623336A (zh) | 基于hvdc整流侧变换器优化同步机次同步振荡的方法和装置 |
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 |