CN109038616A - 一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统 - Google Patents
一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统 Download PDFInfo
- Publication number
- CN109038616A CN109038616A CN201810847521.9A CN201810847521A CN109038616A CN 109038616 A CN109038616 A CN 109038616A CN 201810847521 A CN201810847521 A CN 201810847521A CN 109038616 A CN109038616 A CN 109038616A
- Authority
- CN
- China
- Prior art keywords
- difference
- formula
- generator
- additional
- tune
- 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.)
- Pending
Links
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/24—Arrangements for preventing or reducing oscillations of power in networks
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统,其中方法包括:对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲‑海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K′5、K′6表示的考虑附加调差后的调差模型;根据无调差时单机无穷大系统的状态方程组的特征值,求解无调差时单机无穷大系统的振荡频率;根据考虑附加调差时单机无穷大系统的状态方程组的特征值,求解考虑附加调差时单机无穷大系统的振荡频率;并计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值;比较无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值。
Description
技术领域
本发明涉及电力系统技术领域,更具体地,涉及一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统。
背景技术
随着特高压交直流大型互联电力系统的发展,电力系统的安全稳定运行日益重要。发电机励磁系统对保证系统电压和无功稳定具有十分显著的作用。现在大型电厂大多为发电机变压器组接线方式,且在主变压器的高压侧均并联于同一条母线。所以,同一电厂内一台机组的励磁电压发生改变,不仅会改变本机组的无功及电压,而且还会影响其它并列运行机组的无功,从而引起母线电压的变化。另外,由于主变压器本身存在较大漏抗,也会影响发电机励磁系统对系统电压的作用效果。为了改善发电机励磁系统对系统无功及电压的控制效果,励磁调节器中引入了附加调差,附加调差既可以提高电力系统的电压稳定性,又可以保证机组间无功功率的合理分配。现有技术采用适当的附加调差系数不仅是机组经济运行的要求,也是电网稳定的必要措施,同时也是新投机组的必须试验,有很大的应用价值。
附加调差改善电压的同时也会对发电机的动态阻尼产生影响,然而,现在工程现场还没有衡量附加调差对发电机动态阻尼影响的量化方法,使技术人员设定附加调差系数时能兼顾对阻尼的影响。
因此,需要一种技术,以实现对励磁附加调差对同步发电机动态阻尼影响的计算。
发明内容
本发明技术方案提供一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统,以解决如何计算励磁附加调差对同步发电机动态阻尼影响的问题。
为了解决上述问题,本发明提供了一种计算励磁附加调差对同步发电机动态阻尼影响的方法,所述方法包括:
对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型;
根据无调差时单机无穷大系统的状态方程组的特征值,求解无调差时所述单机无穷大系统的振荡频率;根据考虑附加调差时单机无穷大系统的状态方程组的特征值,求解所述考虑附加调差时单机无穷大系统的振荡频率;
根据所述振荡频率无调差时所述单机无穷大系统的振荡频率,计算无附加调差Xc=0时励磁系统阻尼转矩系数KD;根据所述考虑附加调差时单机无穷大系统的振荡频率,计算附加调差Xc≠0后励磁系统阻尼转矩系数KD的值;
比较所述无附加调差Xc=0时励磁系统阻尼转矩系数KD和所述附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据所述差值,确定励磁附加调差对同步发电机动态阻尼的影响。
优选地,所述对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型,包括:
根据单机无穷大系统中发电机各电气量的向量关系可得到无功功率的表达式:
式(1)中:Qe为发电机无功;utq和utd分别为发电机端电压Ut在q轴和d轴的分量;id和iq分别是机端电流在q轴和d轴的分量;xd′是发电机d轴暂态电抗;xe为发电机外部电抗;x′d∑=x′d+xe;Eq′是发电机q轴暂态电势;Us为无穷大母线电压;δ为发电机q轴与Us的夹角,即发电机功角;xq是发电机q轴电抗;xq∑=xq+xe;
式(1)写成δ与Eq′的偏差方程的形式为:
ΔQe=K11Δδ+K12ΔEq′ (2)
式(2)中:ΔQe为无功变化量;Δδ为功角变化量;ΔEq′为发电机q轴暂态电势变化量;K11、K12的计算公式为:
式(3)中:E′q0为发电机q轴暂态电势初始值;δ0为发电机功角初始值;发电机并网状态下,机端电压保持在额定值附近,近似认为发电机无功电流与发电机无功成同比例变化;则无功调差公式用下式表示:
Ut′=Ut+XcQe (4)
式中Xc为附加调差系数;对式(4)求微分可得ΔUt′的表达式,并设
ΔUt′=K5′Δδ+K6′ΔEq′ (5)
其中:
考虑附加调差后的调差模型就是将原菲利蒲-海佛隆数学模型的K5和K6分别用K5′和K6′代替,原模型的发电机端电压变化量ΔUt用ΔUt′代替。
优选地,所述根据单机无穷大系统的状态方程组的特征值,求解单机无穷大系统的振荡频率,根据所述振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值,包括:
设励磁系统传递函数为无附加调差时励磁系统附加转矩的表达式为:
式(7)中:KA和TA分别为励磁系统静态放大倍数和时间常数;T′d0为发电机开路d轴暂态时间常数;假设系统作小振幅正弦振荡,振荡角频率为ω1,则以s=jω1代入式(7)求出无附加调差时励磁系统阻尼转矩系数
式中ω0=2πf0,中国电网f0=50Hz,ω1为无附加调差时系统振荡频率;
引入附加调差之后,振荡频率会有所变化;同理以s=jω2代入式(7)并分别用K5′和K6′代替K5和K6,便可得到引入附加调差Xc≠0后励磁系统阻尼转矩系数KD的表达式,见式(9);
式中ω0=2πf0,中国电网f0=50Hz,ω2为有附加调差时系统的振荡频率。
优选地,比较所述无附加调差Xc=0时励磁系统阻尼转矩系数KD和所述附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据所述差值,确定励磁附加调差对同步发电机动态阻尼的影响,包括:
单机-无穷大母线系统,给定系统电压Vs,发电机电压Vt、发电机有功P和无功Q,发电机电抗xd、xd'、xq,外部电抗xe,计算出E'q0、sinδ0,进而计算菲利蒲-海佛隆数学模型系数的K1、K2、K3、K4、K5、K6,给定附加调差系数Xc后,由式(3)计算模型系数K11、K12,进而通过式(6)得到系数K5′和K6′;
设发电机采用理想快速励磁系统(TA=0)传递函数为KA,不考虑发电机阻尼绕组,即D=0,根据菲利蒲-海佛隆数学模型框图可写出无附加调差的系统状态方程组:
设(10)式的系数矩阵为A;可根据下式求得矩阵A的特征值:
A的特征值中有两个互为共轭的复数λ1,λ2和一个实数,其中:
λ1,λ2=α1±jω1
α1为衰减系数、ω1为阻尼振荡频率、ζ为阻尼比;
将ω1代入式(8),求得无附加调差Xc=0时的阻尼转矩系数KD;
同理可写出采用附加调差时的单机无穷大系统的状态方程组:
设其系数矩阵为A’;
根据相同步骤求得振荡频率ω2;将ω2代入式(9),求得采用附加调差Xc≠0后的阻尼转矩系数KD;
附加调差对发电机动态阻尼的影响表达式,由附加调差带来的阻尼转矩系数增量ΔKD为
将式(6)代入采用附加调差Xc≠0后的阻尼转矩系数KD的表达式(9),再将式(9)与无附加调差Xc=0时的阻尼转矩系数KD的表达式(8)代入式(14),则ΔKD的表达式为:
设a1、a2、b和c的表达式为:
则式(15)中d、e、f和g的表达式为:
对一定的发电机运行点,在Xc-ΔKD直角坐标平面上,式(15)是一条二次抛物线;e>0时抛物线开口向下,e<0时抛物线开口向上。
基于本发明的另一方面,提供一种计算励磁附加调差对同步发电机动态阻尼影响的系统,所述系统包括:
建立单元,用于对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型;
第一计算单元,用于根据单机无穷大系统的状态方程组的特征值,求解所述单机无穷大系统的振荡频率;
第二计算单元,用于根据所述振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值;
第三计算单元,用于比较所述无附加调差Xc=0时励磁系统阻尼转矩系数KD和所述附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据所述差值,确定励磁附加调差对同步发电机动态阻尼的影响。
优选地,所述建立单元用于对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型,包括:
根据单机无穷大系统中发电机各电气量的向量关系可得到无功功率的表达式:
式(1)中:Qe为发电机无功;utq和utd分别为发电机端电压Ut在q轴和d轴的分量;id和iq分别是机端电流在q轴和d轴的分量;xd′是发电机d轴暂态电抗;xe为发电机外部电抗;x′d∑=x′d+xe;Eq′是发电机q轴暂态电势;Us为无穷大母线电压;δ为发电机q轴与Us的夹角,即发电机功角;xq是发电机q轴电抗;xq∑=xq+xe;
式(1)写成δ与Eq′的偏差方程的形式为:
ΔQe=K11Δδ+K12ΔEq′ (2)
式(2)中:ΔQe为无功变化量;Δδ为功角变化量;ΔEq′为发电机q轴暂态电势变化量;K11、K12的计算公式为:
式(3)中:E′q0为发电机q轴暂态电势初始值;δ0为发电机功角初始值;发电机并网状态下,机端电压保持在额定值附近,近似认为发电机无功电流与发电机无功成同比例变化;则无功调差公式用下式表示:
Ut′=Ut+XcQe (4)
式中Xc为附加调差系数;对式(4)求微分可得ΔUt′的表达式,并设
ΔUt′=K5′Δδ+K6′ΔEq′ (5)
其中:
考虑附加调差后的调差模型就是将原菲利蒲-海佛隆数学模型的K5和K6分别用K5′和K6′代替,原模型的发电机端电压变化量ΔUt用ΔUt′代替。
优选地,所述第一计算单元用于根据单机无穷大系统的状态方程组的特征值,求解单机无穷大系统的振荡频率;所述第二计算单元用于根据所述振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值包括:
设励磁系统传递函数为无附加调差时励磁系统附加转矩的表达式为:
式(7)中:KA和TA分别为励磁系统静态放大倍数和时间常数;Td′0为发电机开路d轴暂态时间常数;假设系统作小振幅正弦振荡,振荡角频率为ω1,则以s=jω1代入式(7)求出无附加调差时励磁系统阻尼转矩系数
式中ω0=2πf0,中国电网f0=50Hz,ω1为无附加调差时系统振荡频率;
引入附加调差之后,振荡频率会有所变化;同理以s=jω2代入式(7)并分别用K5′和K6′代替K5和K6,便可得到引入附加调差Xc≠0后励磁系统阻尼转矩系数KD的表达式,见式(9);
式中ω0=2πf0,中国电网f0=50Hz,ω2为有附加调差时系统的振荡频率。
优选地,所述第三计算单元用于比较所述无附加调差Xc=0时励磁系统阻尼转矩系数KD和所述附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据所述差值,确定励磁附加调差对同步发电机动态阻尼的影响,包括:
单机-无穷大母线系统,给定系统电压Vs,发电机电压Vt、发电机有功P和无功Q,发电机电抗xd、xd'、xq,外部电抗xe,计算出E'q0、sinδ0,进而计算菲利蒲-海佛隆数学模型系数的K1、K2、K3、K4、K5、K6,给定附加调差系数Xc后,由式(3)计算模型系数K11、K12,进而通过式(6)得到系数K5′和K6′;
设发电机采用理想快速励磁系统(TA=0)传递函数为KA,不考虑发电机阻尼绕组,即D=0,根据菲利蒲-海佛隆数学模型框图可写出无附加调差的系统状态方程组:
设(10)式的系数矩阵为A;可根据下式求得矩阵A的特征值:
A的特征值中有两个互为共轭的复数λ1,λ2和一个实数,其中:
λ1,λ2=α1±jω1
α1为衰减系数、ω1为阻尼振荡频率、ζ为阻尼比;
将ω1代入式(8),求得无附加调差Xc=0时的阻尼转矩系数KD;
同理可写出采用附加调差时的单机无穷大系统的状态方程组:
设其系数矩阵为A’;
根据相同步骤求得振荡频率ω2;将ω2代入式(9),求得采用附加调差Xc≠0后的阻尼转矩系数KD;
附加调差对发电机动态阻尼的影响表达式,由附加调差带来的阻尼转矩系数增量ΔKD为
将式(6)代入采用附加调差Xc≠0后的阻尼转矩系数KD的表达式(9),再将式(9)与无附加调差Xc=0时的阻尼转矩系数KD的表达式(8)代入式(14),则ΔKD的表达式为:
设a1、a2、b和c的表达式为:
则式(15)中d、e、f和g的表达式为:
对一定的发电机运行点,在Xc-ΔKD直角坐标平面上,式(15)是一条二次抛物线;e>0时抛物线开口向下,e<0时抛物线开口向上。
为了在工程现场不同工况下更简易地获得准确的附加调差对发电机动态阻尼的影响的具体数据,本发明技术方案提供一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统,其中方法包括:对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型;根据无调差时单机无穷大系统的状态方程组的特征值,求解无调差时单机无穷大系统的振荡频率;根据考虑附加调差时单机无穷大系统的状态方程组的特征值,求解考虑附加调差时单机无穷大系统的振荡频率;根据振荡频率无调差时单机无穷大系统的振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD,根据考虑附加调差时单机无穷大系统的振荡频率计算附加调差Xc≠0后励磁系统阻尼转矩系数KD的值;比较无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据差值,确定励磁附加调差对同步发电机动态阻尼的影响。本发明技术方案实现了励磁附加调差影响同步发电机动态阻尼的量化,实施过程简单有效,为衡量励磁附加调差对同步发电机动态阻尼的影响提供了解决方案。
附图说明
通过参考下面的附图,可以更为完整地理解本发明的示例性实施方式:
图1为根据本发明优选实施方式的一种计算励磁附加调差对同步发电机动态阻尼影响的方法流程图;
图2为根据本发明优选实施方式的单机无穷大系统的菲利蒲-海佛隆模型框图;
图3为根据本发明优选实施方式的计算用单机-无穷大母线系统结构示意图;以及
图4为根据本发明优选实施方式的一种计算励磁附加调差对同步发电机动态阻尼影响的系统结构图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式,然而,本发明可以用许多不同的形式来实施,并且不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中,相同的单元/元件使用相同的附图标记。
除非另有说明,此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外,可以理解的是,以通常使用的词典限定的术语,应当被理解为与其相关领域的语境具有一致的含义,而不应该被理解为理想化的或过于正式的意义。
图1为根据本发明优选实施方式的一种计算励磁附加调差对同步发电机动态阻尼影响的方法流程图。为实现对励磁附加调差影响同步发电机动态阻尼的量化,本发明实施方式提供了一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统,其中方法包括:对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型;根据无调差时单机无穷大系统的状态方程组的特征值,求解无调差时单机无穷大系统的振荡频率;根据考虑附加调差时单机无穷大系统的状态方程组的特征值,求解考虑附加调差时单机无穷大系统的振荡频率;根据振荡频率无调差时单机无穷大系统的振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD,根据考虑附加调差时单机无穷大系统的振荡频率计算附加调差Xc≠0后励磁系统阻尼转矩系数KD的值;比较无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据差值,确定励磁附加调差对同步发电机动态阻尼的影响。如图1所示,一种计算励磁附加调差对同步发电机动态阻尼影响的方法,方法包括:
优选地,在步骤101:对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型。本申请对菲利蒲-海佛隆模型的拓展,在用单机无穷大系统中发电机电抗、功角、电势、机端电压的函数K1—K6表示的单机无穷大系统常用数学模型(即菲利蒲-海佛隆模型)中,未考虑附加调差。图2为单机无穷大系统的利蒲-海佛隆模型框图。
优选地,本申请对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型,包括:
根据单机无穷大系统中发电机各电气量的向量关系可得到无功功率的表达式:
式(1)中:Qe为发电机无功;utq和utd分别为发电机端电压Ut在q轴和d轴的分量;id和iq分别是机端电流在q轴和d轴的分量;xd′是发电机d轴暂态电抗;xe为发电机外部电抗;x′d∑=x′d+xe;Eq′是发电机q轴暂态电势;Us为无穷大母线电压;δ为发电机q轴与Us的夹角,即发电机功角;xq是发电机q轴电抗;xq∑=xq+xe;
式(1)写成δ与Eq′的偏差方程的形式为:
ΔQe=K11Δδ+K12ΔEq′ (2)
式(2)中:ΔQe为无功变化量;Δδ为功角变化量;ΔEq′为发电机q轴暂态电势变化量;K11、K12的计算公式为:
式(3)中:E′q0为发电机q轴暂态电势初始值;δ0为发电机功角初始值;发电机并网状态下,机端电压保持在额定值附近,近似认为发电机无功电流与发电机无功成同比例变化;则无功调差公式用下式表示:
Ut′=Ut+XcQe (4)
式中Xc为附加调差系数;对式(4)求微分可得ΔUt′的表达式,并设
ΔUt′=K5′Δδ+K6′ΔEq′ (5)
其中:
考虑附加调差后的调差模型就是将原菲利蒲-海佛隆数学模型的K5和K6分别用K5′和K6′代替,原模型的发电机端电压变化量ΔUt用ΔUt′代替。因此考虑附加调差后,并没有改变菲利蒲-海佛隆模型的结构,只是以上3个量发生了变化。
优选地,在步骤102:根据无调差时单机无穷大系统的状态方程组的特征值,求解无调差时单机无穷大系统的振荡频率;根据考虑附加调差时单机无穷大系统的状态方程组的特征值,求解考虑附加调差时单机无穷大系统的振荡频率。
优选地,在步骤103:根据振荡频率无调差时单机无穷大系统的振荡频率,计算无附加调差Xc=0时励磁系统阻尼转矩系数KD;根据考虑附加调差时单机无穷大系统的振荡频率,计算附加调差Xc≠0后励磁系统阻尼转矩系数KD的值。
优选地,根据单机无穷大系统的状态方程组的特征值,求解单机无穷大系统的振荡频率,根据振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值,包括:
设励磁系统传递函数为无附加调差时励磁系统附加转矩的表达式为:
式(7)中:KA和TA分别为励磁系统静态放大倍数和时间常数;T′d0为发电机开路d轴暂态时间常数;假设系统作小振幅正弦振荡,振荡角频率为ω1,则以s=jω1代入式(7)求出无附加调差时励磁系统阻尼转矩系数
式中ω0=2πf0,中国电网f0=50Hz,ω1为无附加调差时系统振荡频率;
引入附加调差之后,振荡频率会有所变化;同理以s=jω2代入式(7)并分别用K5′和K6′代替K5和K6,便可得到引入附加调差Xc≠0后励磁系统阻尼转矩系数KD的表达式,见式(9);
式中ω0=2πf0,中国电网f0=50Hz,ω2为有附加调差时系统的振荡频率。
优选地,在步骤104:比较无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据差值,确定励磁附加调差对同步发电机动态阻尼的影响。
如果已知振荡频率,则可以计算、比较KD(Xc=0)与KD(Xc≠0)的大小。但是,振荡频率本身又是阻尼力矩系数和同步力矩系数的函数。本申请将通过求解系统特征方程的特征值来求得振荡频率。
优选地,比较无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据差值,确定励磁附加调差对同步发电机动态阻尼的影响,包括:
如图3为单机-无穷大母线系统,给定系统电压Vs,发电机电压Vt、发电机有功P和无功Q,发电机电抗xd、xd'、xq,外部电抗xe,计算出E'q0、sinδ0,进而计算菲利蒲-海佛隆数学模型系数的K1、K2、K3、K4、K5、K6,给定附加调差系数Xc后,由式(3)计算模型系数K11、K12,进而通过式(6)得到系数K5′和K6′;
设发电机采用理想快速励磁系统(TA=0)传递函数为KA,不考虑发电机阻尼绕组,即D=0,根据菲利蒲-海佛隆数学模型框图可写出无附加调差的系统状态方程组:
设(10)式的系数矩阵为A;可根据下式求得矩阵A的特征值:
A的特征值中有两个互为共轭的复数λ1,λ2和一个实数,其中:
λ1,λ2=α1±jω1
α1为衰减系数、ω1为阻尼振荡频率、ζ为阻尼比;
将ω1代入式(8),求得无附加调差Xc=0时的阻尼转矩系数KD;
同理可写出采用附加调差时的单机无穷大系统的状态方程组:
设其系数矩阵为A’;
根据相同步骤求得振荡频率ω2;将ω2代入式(9),求得采用附加调差Xc≠0后的阻尼转矩系数KD;
附加调差对发电机动态阻尼的影响表达式,由附加调差带来的阻尼转矩系数增量ΔKD为
将式(6)代入采用附加调差Xc≠0后的阻尼转矩系数KD的表达式(9),再将式(9)与无附加调差Xc=0时的阻尼转矩系数KD的表达式(8)代入式(14),则ΔKD的表达式为:
设a1、a2、b和c的表达式为:
则式(15)中d、e、f和g的表达式为:
对一定的发电机运行点,在Xc-ΔKD直角坐标平面上,式(15)是一条二次抛物线;e>0时抛物线开口向下,e<0时抛物线开口向上。
本申请提供的技术方案中,基于扩展的菲利蒲-海佛隆模型及其特征方程的特征值,推导出的阻尼转矩公式简单有效,适合实际工程应用。本申请提供的技术方案对工况无特殊要求,并网情况下适合所有工况,计算结果鲁棒性较好。本申请提供的技术方案中,通过发电机组实例分析结果,验证了衡量励磁附加调差对同步发电机动态阻尼影响的量化方法的准确性,显示出该实测方法具有较强的工程实用性。
以下通过表1所示的某647MW火电机组测试结果为例对本申请实施方式进一步的详细说明,但本申请不限于所给出的实施方式。
表1 不同调差系数下机组1的动态阻尼特性
附加调差Xc(p.u.) | 阻尼转矩系数增量△KD(p.u.) | 阻尼比 |
-0.134 | -6.14623 | 0.00093 |
-0.1206 | -5.49319 | 0.006168 |
-0.1072 | -4.8303 | 0.011351 |
-0.0938 | -4.1642 | 0.016445 |
-0.067 | -2.8509 | 0.026233 |
-0.0536 | -2.2181 | 0.030865 |
-0.0402 | -1.6107 | 0.035282 |
-0.0268 | -1.035 | 0.039458 |
-0.0134 | -0.4966 | 0.043369 |
0 | 0 | 0.046994 |
0.0134 | 0.451 | 0.050318 |
0.0268 | 0.8542 | 0.053329 |
0.0402 | 1.2085 | 0.056019 |
0.067 | 1.7708 | 0.060436 |
0.0938 | 2.149 | 0.063612 |
0.1072 | 2.2757 | 0.064766 |
0.1206 | 2.3651 | 0.065655 |
0.134 | 2.421 | 0.066299 |
使用本申请提出的方法进行励磁附加调差对同步发电机动态阻尼影响的量化计算,步骤如下:
步骤一:将附加调差引入单机无穷大系统的菲利蒲-海佛隆模型,求得扩展后的菲利蒲-海佛隆模型系数;
根据系统电压Vs、发电机电压Vt、发电机有功P和无功Q,发电机电抗xd、xd'、xq,外部电抗xe,可以计算出E'q0、sinδ0,进而计算菲利蒲-海佛隆模型系数K1~K6。
给定附加调差系数Xc后,由下式计算模型系数K11、K12,
进而通过下式得到系数K5′和K6′
步骤二:根据扩展后的菲利蒲-海佛隆模型,求得有、无附加调差时励磁系统所提供阻尼转矩系数KD(Xc≠0)和KD(Xc=0)的表达式,但此时系统振荡频率仍然未知;
两个阻尼转矩系数公式如下:
步骤三:根据菲利蒲-海佛隆模型求得有、无附加调差时的系统特征方程,通过求解系统特征方程的特征值来求得有、无附加调差时的振荡频率ω1和ω2。
振荡频率ω1计算公式如下:
无附加调差的系统状态方程组为
上述状态方程组的系数矩阵A,根据下式求得矩阵A的特征值
A的特征值中有两个互为共轭的复数λ1,λ2和一个实数。其中:
λ1,λ2=α1±jω1
其中,ω1为无附加调差时的系统振荡频率、ζ为阻尼比、α1为衰减系数;
同理,振荡频率ω2的计算公式如下:
有附加调差的系统状态方程组为:
上述状态方程组的系数矩阵A’,根据下式求得矩阵A’的特征值
A’的特征值中有两个互为共轭的复数λ4,λ5和一个实数。其中:
λ4,λ5=α2±jω2
其中,ω2为有附加调差时的系统振荡频率、ζ为阻尼比、α2为衰减系数;
步骤四:根据有、无附加调差时的励磁系统阻尼转矩系数KD(Xc≠0)和KD(Xc=0)的表达式,两者相减,可以求得由附加调差带来的阻尼转矩系数增量ΔKD。
由附加调差带来的阻尼转矩系数增量ΔKD公式如下:
其中,K2、K3是菲利蒲-海佛隆模型的系数,ω0=2πf0,中国电网f0=50Hz,
且:
步骤五:将系统电压、发电机参数、励磁系统参数、线路参数代入ΔKD的表达式,就可以计算由附加调差带来的阻尼转矩系数增量,从而实现了阻尼的量化;
表1计算结果采用的相关参数为:励磁系统动态放大倍数KA=50,励磁系统时间常数TA=0.02s;发电机参数Xd=Xq=2.04,Xd′=0.29,T′d0=6.6s,发电机惯性常数H=3.846s;主变压器电抗XT=0.134;;500kV等级的线路电抗XL=0.005;发电机的额定容量SN=719MVA,发电机运行有功P=590MW,发电机运行无功Q=50MVar;电力系统电压其中,发电机、主变、线路的电抗标幺基值均为发电机的额定容量SN。
由测试结果可知,使用本申请设计的衡量励磁附加调差对同步发电机动态阻尼影响的量化方法,可以快速估算附加调差对发电机系统阻尼的影响,为现场设置附加调差参数值提供参考。且因本申请提供的方法具有多工况适应的特点,在并网运行情况下具有简单方便、计算结果鲁棒性好的效果,可以满足电网运行需求,从而验证了本申请提供的方法在实际系统分析中的有效性。
图4为根据本发明优选实施方式的一种计算励磁附加调差对同步发电机动态阻尼影响的系统结构图。如图4所示,一种计算励磁附加调差对同步发电机动态阻尼影响的系统,系统包括:
建立单元401,用于对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型。
优选地,建立单元401用于对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K5′、K6′表示的考虑附加调差后的调差模型,包括:
根据单机无穷大系统中发电机各电气量的向量关系可得到无功功率的表达式:
式(1)中:Qe为发电机无功;utq和utd分别为发电机端电压Ut在q轴和d轴的分量;id和iq分别是机端电流在q轴和d轴的分量;xd′是发电机d轴暂态电抗;xe为发电机外部电抗;x′d∑=x′d+xe;Eq′是发电机q轴暂态电势;Us为无穷大母线电压;δ为发电机q轴与Us的夹角,即发电机功角;xq是发电机q轴电抗;xq∑=xq+xe;
式(1)写成δ与Eq′的偏差方程的形式为:
ΔQe=K11Δδ+K12ΔEq′ (2)
式(2)中:ΔQe为无功变化量;Δδ为功角变化量;ΔEq′为发电机q轴暂态电势变化量;K11、K12的计算公式为:
式(3)中:E′q0为发电机q轴暂态电势初始值;δ0为发电机功角初始值;发电机并网状态下,机端电压保持在额定值附近,近似认为发电机无功电流与发电机无功成同比例变化;则无功调差公式用下式表示:
Ut′=Ut+XcQe (4)
式中Xc为附加调差系数;对式(4)求微分可得ΔUt′的表达式,并设
ΔUt′=K5′Δδ+K6′ΔEq′ (5)
其中:
考虑附加调差后的调差模型就是将原菲利蒲-海佛隆数学模型的K5和K6分别用K5′和K6′代替,原模型的发电机端电压变化量ΔUt用ΔUt′代替。
第一计算单元402,用于根据单机无穷大系统的状态方程组的特征值,求解单机无穷大系统的振荡频率。
第二计算单元403,用于根据振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值。
第三计算单元404,用于比较无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据差值,确定励磁附加调差对同步发电机动态阻尼的影响。
优选地,第一计算单元402用于根据单机无穷大系统的状态方程组的特征值,求解单机无穷大系统的振荡频率;第二计算单元403用于根据振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值包括:
设励磁系统传递函数为无附加调差时励磁系统附加转矩的表达式为:
式(7)中:KA和TA分别为励磁系统静态放大倍数和时间常数;T′d0为发电机开路d轴暂态时间常数;假设系统作小振幅正弦振荡,振荡角频率为ω1,则以s=jω1代入式(7)求出无附加调差时励磁系统阻尼转矩系数
式中ω0=2πf0,中国电网f0=50Hz,ω1为无附加调差时系统振荡频率;
引入附加调差之后,振荡频率会有所变化;同理以s=jω2代入式(7)并分别用K5′和K6′代替K5和K6,便可得到引入附加调差Xc≠0后励磁系统阻尼转矩系数KD的表达式,见式(9);
式中ω0=2πf0,中国电网f0=50Hz,ω2为有附加调差时系统的振荡频率。
优选地,第三计算单元404用于比较无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据差值,确定励磁附加调差对同步发电机动态阻尼的影响,包括:
单机-无穷大母线系统,给定系统电压Vs,发电机电压Vt、发电机有功P和无功Q,发电机电抗xd、xd'、xq,外部电抗xe,计算出E'q0、sinδ0,进而计算菲利蒲-海佛隆数学模型系数的K1、K2、K3、K4、K5、K6,给定附加调差系数Xc后,由式(3)计算模型系数K11、K12,进而通过式(6)得到系数K5′和K6′;
设发电机采用理想快速励磁系统(TA=0)传递函数为KA,不考虑发电机阻尼绕组,即D=0,根据菲利蒲-海佛隆数学模型框图可写出无附加调差的系统状态方程组:
设(10)式的系数矩阵为A;可根据下式求得矩阵A的特征值:
A的特征值中有两个互为共轭的复数λ1,λ2和一个实数,其中:
λ1,λ2=α1±jω1
α1为衰减系数、ω1为阻尼振荡频率、ζ为阻尼比;
将ω1代入式(8),求得无附加调差Xc=0时的阻尼转矩系数KD;
同理可写出采用附加调差时的单机无穷大系统的状态方程组:
设其系数矩阵为A’;
根据相同步骤求得振荡频率ω2;将ω2代入式(9),求得采用附加调差Xc≠0后的阻尼转矩系数KD;
附加调差对发电机动态阻尼的影响表达式,由附加调差带来的阻尼转矩系数增量ΔKD为
将式(6)代入采用附加调差Xc≠0后的阻尼转矩系数KD的表达式(9),再将式(9)与无附加调差Xc=0时的阻尼转矩系数KD的表达式(8)代入式(14),则ΔKD的表达式为:
设a1、a2、b和c的表达式为:
则式(15)中d、e、f和g的表达式为:
对一定的发电机运行点,在Xc-ΔKD直角坐标平面上,式(15)是一条二次抛物线;e>0时抛物线开口向下,e<0时抛物线开口向上。
本发明优选实施方式的一种计算励磁附加调差对同步发电机动态阻尼影响的系统400与本发明另一优选实施方式的一种计算励磁附加调差对同步发电机动态阻尼影响的方法100相对应,在此不再进行赘述。
已经通过参考少量实施方式描述了本发明。然而,本领域技术人员所公知的,正如附带的专利权利要求所限定的,除了本发明以上公开的其他的实施例等同地落在本发明的范围内。
通常地,在权利要求中使用的所有术语都根据他们在技术领域的通常含义被解释,除非在其中被另外明确地定义。所有的参考“一个/所述/该[装置、组件等]”都被开放地解释为所述装置、组件等中的至少一个实例,除非另外明确地说明。这里公开的任何方法的步骤都没必要以公开的准确的顺序运行,除非明确地说明。
Claims (8)
1.一种计算励磁附加调差对同步发电机动态阻尼影响的方法,所述方法包括:
对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K′5、K′6表示的考虑附加调差后的调差模型;
根据无调差时单机无穷大系统的状态方程组的特征值,求解无调差时所述单机无穷大系统的振荡频率;根据考虑附加调差时单机无穷大系统的状态方程组的特征值,求解所述考虑附加调差时单机无穷大系统的振荡频率;
根据所述振荡频率无调差时所述单机无穷大系统的振荡频率,计算无附加调差Xc=0时励磁系统阻尼转矩系数KD;根据所述考虑附加调差时单机无穷大系统的振荡频率,计算附加调差Xc≠0后励磁系统阻尼转矩系数KD的值;
比较所述无附加调差Xc=0时励磁系统阻尼转矩系数KD和所述附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据所述差值,确定励磁附加调差对同步发电机动态阻尼的影响。
2.根据权利要求1所述的方法,所述对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K′5、K′6表示的考虑附加调差后的调差模型,包括:
根据单机无穷大系统中发电机各电气量的向量关系可得到无功功率的表达式:
式(1)中:Qe为发电机无功;utq和utd分别为发电机端电压Ut在q轴和d轴的分量;id和iq分别是机端电流在q轴和d轴的分量;xd′是发电机d轴暂态电抗;xe为发电机外部电抗;x′d∑=x′d+xe;E′q是发电机q轴暂态电势;Us为无穷大母线电压;δ为发电机q轴与Us的夹角,即发电机功角;xq是发电机q轴电抗;xq∑=xq+xe;
式(1)写成δ与E′q的偏差方程的形式为:
ΔQe=K11Δδ+K12ΔE′q (2)
式(2)中:ΔQe为无功变化量;Δδ为功角变化量;ΔE′q为发电机q轴暂态电势变化量;K11、K12的计算公式为:
式(3)中:E′q0为发电机q轴暂态电势初始值;δ0为发电机功角初始值;发电机并网状态下,机端电压保持在额定值附近,近似认为发电机无功电流与发电机无功成同比例变化;则无功调差公式用下式表示:
U′t=Ut+XcQe (4)
式中Xc为附加调差系数;对式(4)求微分可得ΔU′t的表达式,并设
ΔU′t=K′5Δδ+K′6ΔE′q (5)
其中:
考虑附加调差后的调差模型就是将原菲利蒲-海佛隆数学模型的K5和K6分别用K′5和K′6代替,原模型的发电机端电压变化量ΔUt用ΔU′t代替。
3.根据权利要求2所述的方法,所述根据单机无穷大系统的状态方程组的特征值,求解单机无穷大系统的振荡频率,根据所述振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值,包括:
设励磁系统传递函数为无附加调差时励磁系统附加转矩的表达式为:
式(7)中:KA和TA分别为励磁系统静态放大倍数和时间常数;T′d0为发电机开路d轴暂态时间常数;假设系统作小振幅正弦振荡,振荡角频率为ω1,则以s=jω1代入式(7)求出无附加调差时励磁系统阻尼转矩系数
式中ω0=2πf0,中国电网f0=50Hz,ω1为无附加调差时系统振荡频率;
引入附加调差之后,振荡频率会有所变化;同理以s=jω2代入式(7)并分别用K′5和K′6代替K5和K6,便可得到引入附加调差Xc≠0后励磁系统阻尼转矩系数KD的表达式,见式(9);
式中ω0=2πf0,中国电网f0=50Hz,ω2为有附加调差时系统的振荡频率。
4.根据权利要求3所述的方法,比较所述无附加调差Xc=0时励磁系统阻尼转矩系数KD和所述附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据所述差值,确定励磁附加调差对同步发电机动态阻尼的影响,包括:
单机-无穷大母线系统,给定系统电压Vs,发电机电压Vt、发电机有功P和无功Q,发电机电抗xd、xd'、xq,外部电抗xe,计算出E′q0、sinδ0,进而计算菲利蒲-海佛隆数学模型系数的K1、K2、K3、K4、K5、K6,给定附加调差系数Xc后,由式(3)计算模型系数K11、K12,进而通过式(6)得到系数K′5和K′6;
设发电机采用理想快速励磁系统TA=0时传递函数为KA,不考虑发电机阻尼绕组,即D=0,根据菲利蒲-海佛隆数学模型框图可写出无附加调差的系统状态方程组:
设(10)式的系数矩阵为A;可根据下式求得矩阵A的特征值:
A的特征值中有两个互为共轭的复数λ1,λ2和一个实数,其中:
λ1,λ2=α1±jω1
α1为衰减系数、ω1为阻尼振荡频率、ζ为阻尼比;
将ω1代入式(8),求得无附加调差Xc=0时的阻尼转矩系数KD;
同理可写出采用附加调差时的单机无穷大系统的状态方程组:
设其系数矩阵为A’;
根据相同步骤求得振荡频率ω2;将ω2代入式(9),求得采用附加调差Xc≠0后的阻尼转矩系数KD;
附加调差对发电机动态阻尼的影响表达式,由附加调差带来的阻尼转矩系数增量ΔKD为
将式(6)代入采用附加调差Xc≠0后的阻尼转矩系数KD的表达式(9),再将式(9)与无附加调差Xc=0时的阻尼转矩系数KD的表达式(8)代入式(14),则ΔKD的表达式为:
设a1、a2、b和c的表达式为:
则式(15)中d、e、f和g的表达式为:
对一定的发电机运行点,在Xc-ΔKD直角坐标平面上,式(15)是一条二次抛物线;e>0时抛物线开口向下,e<0时抛物线开口向上。
5.一种计算励磁附加调差对同步发电机动态阻尼影响的系统,所述系统包括:
建立单元,用于对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K′5、K′6表示的考虑附加调差后的调差模型;
第一计算单元,用于根据单机无穷大系统的状态方程组的特征值,求解所述单机无穷大系统的振荡频率;
第二计算单元,用于根据所述振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值;
第三计算单元,用于比较所述无附加调差Xc=0时励磁系统阻尼转矩系数KD和所述附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据所述差值,确定励磁附加调差对同步发电机动态阻尼的影响。
6.根据权利要求5所述的系统,所述建立单元用于对K1、K2、K3、K4、K5、K6表示的单机无穷大系统菲利蒲-海佛隆数学模型进行拓展,建立K1、K2、K3、K4、K′5、K′6表示的考虑附加调差后的调差模型,包括:
根据单机无穷大系统中发电机各电气量的向量关系可得到无功功率的表达式:
式(1)中:Qe为发电机无功;utq和utd分别为发电机端电压Ut在q轴和d轴的分量;id和iq分别是机端电流在q轴和d轴的分量;xd′是发电机d轴暂态电抗;xe为发电机外部电抗;x′d∑=x′d+xe;E′q是发电机q轴暂态电势;Us为无穷大母线电压;δ为发电机q轴与Us的夹角,即发电机功角;xq是发电机q轴电抗;xq∑=xq+xe;
式(1)写成δ与E′q的偏差方程的形式为:
ΔQe=K11Δδ+K12ΔE′q (2)
式(2)中:ΔQe为无功变化量;Δδ为功角变化量;ΔE′q为发电机q轴暂态电势变化量;K11、K12的计算公式为:
式(3)中:E′q0为发电机q轴暂态电势初始值;δ0为发电机功角初始值;发电机并网状态下,机端电压保持在额定值附近,近似认为发电机无功电流与发电机无功成同比例变化;则无功调差公式用下式表示:
U′t=Ut+XcQe (4)
式中Xc为附加调差系数;对式(4)求微分可得ΔU′t的表达式,并设
ΔU′t=K′5Δδ+K′6ΔE′q (5)
其中:
考虑附加调差后的调差模型就是将原菲利蒲-海佛隆数学模型的K5和K6分别用K′5和K′6代替,原模型的发电机端电压变化量ΔUt用ΔU′t代替。
7.根据权利要求6所述的系统,所述第一计算单元用于根据单机无穷大系统的状态方程组的特征值,求解单机无穷大系统的振荡频率;所述第二计算单元用于根据所述振荡频率,分别计算无附加调差Xc=0时励磁系统阻尼转矩系数KD和附加调差Xc≠0后励磁系统阻尼转矩系数KD的值包括:
设励磁系统传递函数为无附加调差时励磁系统附加转矩的表达式为:
式(7)中:KA和TA分别为励磁系统静态放大倍数和时间常数;T′d0为发电机开路d轴暂态时间常数;假设系统作小振幅正弦振荡,振荡角频率为ω1,则以s=jω1代入式(7)求出无附加调差时励磁系统阻尼转矩系数
式中ω0=2πf0,中国电网f0=50Hz,ω1为无附加调差时系统振荡频率;
引入附加调差之后,振荡频率会有所变化;同理以s=jω2代入式(7)并分别用K′5和K′6代替K5和K6,便可得到引入附加调差Xc≠0后励磁系统阻尼转矩系数KD的表达式,见式(9);
式中ω0=2πf0,中国电网f0=50Hz,ω2为有附加调差时系统的振荡频率。
8.根据权利要求7所述的系统,所述第三计算单元用于比较所述无附加调差Xc=0时励磁系统阻尼转矩系数KD和所述附加调差Xc≠0后励磁系统阻尼转矩系数KD的差值,根据所述差值,确定励磁附加调差对同步发电机动态阻尼的影响,包括:
单机-无穷大母线系统,给定系统电压Vs,发电机电压Vt、发电机有功P和无功Q,发电机电抗xd、xd'、xq,外部电抗xe,计算出E′q0、sinδ0,进而计算菲利蒲-海佛隆数学模型系数的K1、K2、K3、K4、K5、K6,给定附加调差系数Xc后,由式(3)计算模型系数K11、K12,进而通过式(6)得到系数K′5和K′6;
设发电机采用理想快速励磁系统(TA=0)传递函数为KA,不考虑发电机阻尼绕组,即D=0,根据菲利蒲-海佛隆数学模型框图可写出无附加调差的系统状态方程组:
设(10)式的系数矩阵为A;可根据下式求得矩阵A的特征值:
A的特征值中有两个互为共轭的复数λ1,λ2和一个实数,其中:
λ1,λ2=α1±jω1
α1为衰减系数、ω1为阻尼振荡频率、ζ为阻尼比;
将ω1代入式(8),求得无附加调差Xc=0时的阻尼转矩系数KD;
同理可写出采用附加调差时的单机无穷大系统的状态方程组:
设其系数矩阵为A’;
根据相同步骤求得振荡频率ω2;将ω2代入式(9),求得采用附加调差Xc≠0后的阻尼转矩系数KD;
附加调差对发电机动态阻尼的影响表达式,由附加调差带来的阻尼转矩系数增量ΔKD为
将式(6)代入采用附加调差Xc≠0后的阻尼转矩系数KD的表达式(9),再将式(9)与无附加调差Xc=0时的阻尼转矩系数KD的表达式(8)代入式(14),则ΔKD的表达式为:
设a1、a2、b和c的表达式为:
则式(15)中d、e、f和g的表达式为:
对一定的发电机运行点,在Xc-ΔKD直角坐标平面上,式(15)是一条二次抛物线;e>0时抛物线开口向下,e<0时抛物线开口向上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810847521.9A CN109038616A (zh) | 2018-07-27 | 2018-07-27 | 一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810847521.9A CN109038616A (zh) | 2018-07-27 | 2018-07-27 | 一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109038616A true CN109038616A (zh) | 2018-12-18 |
Family
ID=64647341
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810847521.9A Pending CN109038616A (zh) | 2018-07-27 | 2018-07-27 | 一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109038616A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110994642A (zh) * | 2019-11-22 | 2020-04-10 | 中国电力科学研究院有限公司 | 励磁增益对凸极发电机同步转矩影响的量化方法及装置 |
CN111555312A (zh) * | 2020-05-27 | 2020-08-18 | 四川大学 | 一种适用于电力系统超低频振荡稳定性评估的方法 |
CN111817312A (zh) * | 2020-05-26 | 2020-10-23 | 中国电力科学研究院有限公司 | 基于励磁动态增益量化调相机暂态稳定影响的方法及系统 |
CN112398390A (zh) * | 2020-10-22 | 2021-02-23 | 国网江西省电力有限公司电力科学研究院 | 一种发电机组励磁系统参数优化的方法 |
CN112834925A (zh) * | 2021-01-29 | 2021-05-25 | 中国电力科学研究院有限公司 | 一种衡量电压测量时间常数对发电机阻尼影响的方法 |
CN112968447A (zh) * | 2021-01-29 | 2021-06-15 | 中国电力科学研究院有限公司 | 一种衡量电压测量时间常数对机组动态稳定影响的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354332A (zh) * | 2011-09-16 | 2012-02-15 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种用于简化柔性交直流输电系统中rga计算的方法 |
JP5330082B2 (ja) * | 2009-05-11 | 2013-10-30 | 株式会社東芝 | 同期発電機の励磁制御装置 |
CN103401497A (zh) * | 2013-07-10 | 2013-11-20 | 国家电网公司 | 基于提高机组功角稳定性的励磁附加调差系数整定方法 |
CN104979841A (zh) * | 2014-04-08 | 2015-10-14 | 国家电网公司 | 一种平抑间歇式能源出力随机波动的方法 |
-
2018
- 2018-07-27 CN CN201810847521.9A patent/CN109038616A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5330082B2 (ja) * | 2009-05-11 | 2013-10-30 | 株式会社東芝 | 同期発電機の励磁制御装置 |
CN102354332A (zh) * | 2011-09-16 | 2012-02-15 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种用于简化柔性交直流输电系统中rga计算的方法 |
CN103401497A (zh) * | 2013-07-10 | 2013-11-20 | 国家电网公司 | 基于提高机组功角稳定性的励磁附加调差系数整定方法 |
CN104979841A (zh) * | 2014-04-08 | 2015-10-14 | 国家电网公司 | 一种平抑间歇式能源出力随机波动的方法 |
Non-Patent Citations (3)
Title |
---|
霍承祥: ""高压侧电压控制的分析和研究"", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
霍承祥等: ""励磁系统附加调差对发电机阻尼特性影响的机制分析及试验"", 《电网技术》 * |
霍承祥等: ""励磁附加调差对电力系统阻尼特性影响的机理分析"", 《大电机技术》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110994642A (zh) * | 2019-11-22 | 2020-04-10 | 中国电力科学研究院有限公司 | 励磁增益对凸极发电机同步转矩影响的量化方法及装置 |
CN111817312A (zh) * | 2020-05-26 | 2020-10-23 | 中国电力科学研究院有限公司 | 基于励磁动态增益量化调相机暂态稳定影响的方法及系统 |
CN111817312B (zh) * | 2020-05-26 | 2022-05-17 | 中国电力科学研究院有限公司 | 基于励磁动态增益量化调相机暂态稳定影响的方法及系统 |
CN111555312A (zh) * | 2020-05-27 | 2020-08-18 | 四川大学 | 一种适用于电力系统超低频振荡稳定性评估的方法 |
CN111555312B (zh) * | 2020-05-27 | 2021-06-22 | 四川大学 | 一种适用于电力系统超低频振荡稳定性评估的方法 |
CN112398390A (zh) * | 2020-10-22 | 2021-02-23 | 国网江西省电力有限公司电力科学研究院 | 一种发电机组励磁系统参数优化的方法 |
CN112834925A (zh) * | 2021-01-29 | 2021-05-25 | 中国电力科学研究院有限公司 | 一种衡量电压测量时间常数对发电机阻尼影响的方法 |
CN112968447A (zh) * | 2021-01-29 | 2021-06-15 | 中国电力科学研究院有限公司 | 一种衡量电压测量时间常数对机组动态稳定影响的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109038616A (zh) | 一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统 | |
CN108964540A (zh) | 一种计算励磁附加调差对同步发电机同步转矩影响的方法及系统 | |
Wehbe et al. | Least squares based estimation of synchronous generator states and parameters with phasor measurement units | |
CN103236702A (zh) | 一种基于阻抗复合控制的微网无功均衡分配方法 | |
CN114342202A (zh) | 用于电网阻抗和动态估计的方法 | |
Roca et al. | System frequency response model considering the influence of power system stabilizers | |
CN105281329B (zh) | 一种提高2机电力系统暂态频率稳定性的单相重合时序整定方法 | |
Ali et al. | A trajectory piecewise-linear approach to nonlinear model order reduction of wind farms | |
Luo et al. | Stability and accuracy considerations in the design and implementation of wind turbine power hardware in the loop platform | |
CN112436520A (zh) | 一种交流电力弹簧反馈线性化解耦控制方法 | |
WO2019174555A1 (zh) | 一种基于直算法的电力系统机电暂态仿真方法 | |
CN106406272A (zh) | 一种风电场中静止无功发生器的控制器性能测试方法 | |
CN109038674B (zh) | 非线性最小二乘曲线拟合的vsg惯量与阻尼系数测量方法 | |
CN105224728B (zh) | 一种含详细发电机模型的电网暂态稳定能量函数分析方法及系统 | |
Miller et al. | Coordinated voltage control for multiple wind plants in Eastern Wyoming: Analysis and field experience | |
Zou et al. | Voltage stability analysis using a complete model of grid-connected voltage-source converters | |
Samal et al. | Damping of power system oscillations by using unified Power Flow Controller with POD and PID controllers | |
CN107465213B (zh) | 一种发电站avc系统及其无功实时调节量计算方法 | |
Elices et al. | Physical interpretation of state feedback controllers to damp power system oscillations | |
CN112968447A (zh) | 一种衡量电压测量时间常数对机组动态稳定影响的方法 | |
Lin et al. | Impact of Different AC Voltage Control Modes of Wind-farm-side MMC on Stability of MMC-HVDC with Offshore Wind Farms | |
Paszek et al. | Influence of subtransient asymmetry of synchronous generators on the power system waveforms under asymmetric operating conditions | |
Zhao et al. | Corrected Complex Torque Coefficient Method for Small-Signal Stability Analysis of Grid-Connected VSC | |
Awadallah et al. | Mathematical model of steam power plant for voltage and reactive power control | |
Li et al. | Transient Fault Current Calculation Method for the PV Grid-connected System considering Dynamic Response of PLL |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20210317 Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant after: STATE GRID ZHEJIANG ELECTRIC POWER Co.,Ltd. Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant before: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. |
|
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181218 |