CN107611999B - 大区交流联络线功率振荡响应的计算方法 - Google Patents

大区交流联络线功率振荡响应的计算方法 Download PDF

Info

Publication number
CN107611999B
CN107611999B CN201710894105.XA CN201710894105A CN107611999B CN 107611999 B CN107611999 B CN 107611999B CN 201710894105 A CN201710894105 A CN 201710894105A CN 107611999 B CN107611999 B CN 107611999B
Authority
CN
China
Prior art keywords
power
alternating current
generator
oscillation response
power 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
Application number
CN201710894105.XA
Other languages
English (en)
Other versions
CN107611999A (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201710894105.XA priority Critical patent/CN107611999B/zh
Publication of CN107611999A publication Critical patent/CN107611999A/zh
Application granted granted Critical
Publication of CN107611999B publication Critical patent/CN107611999B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种大区交流联络线功率振荡响应的计算方法,包括构建简化交流系统模型;基于简化交流系统模型,推导出脉冲功率冲击下交流联络线功率振荡响应公式和阶跃功率冲击下交流联络线功率振荡响应公式;将特高压直流故障下的功率波动分解为若干个脉冲函数和阶跃函数的叠加;计算各脉冲函数下的交流联络线功率振荡响应以及各阶跃函数下的交流联络线功率振荡响应;将所有交流联络线功率振荡响应叠加即为特高压直流故障下的大区交流联络线功率振荡响应。本发明能够比较准确地得出大区交流联络线的功率振荡响应,同时计算方法简单易用,适合于如今对交直流混联电网的交流联络线功率振荡研究。

Description

大区交流联络线功率振荡响应的计算方法
技术领域
本发明涉及一种大区交流联络线功率振荡响应的计算方法,具体涉及一种特高压直流故障下大区交流联络线功率振荡响应的计算方法,属于电力系统分析领域。
背景技术
特高压直流输电(UHVDC)具有输送容量大,输电损耗小,传输距离长等优点,有助于我国实现能源资源的优化配置,因而发展十分迅速,其在我国“西电东送”的战略中发挥了重要的作用。
对交流系统而言,特高压直流系统可看为大容量的电源,当其发生故障后,会对交流系统产生大功率冲击,而大区域之间交流联络线的功率振荡与特高压直流功率波动的冲击密切相关。
所以,若能有一种计算方法,可以离线计算出在特高压直流故障下大区交流联络线功率振荡响应,对于研究交直流电网的交互影响将会十分有帮助。
发明内容
为了解决上述技术问题,本发明提供了一种大区交流联络线功率振荡响应的计算方法。
为了达到上述目的,本发明所采用的技术方案是:
大区交流联络线功率振荡响应的计算方法,包括,
构建简化交流系统模型;
基于简化交流系统模型,推导出脉冲功率冲击下交流联络线功率振荡响应公式和阶跃功率冲击下交流联络线功率振荡响应公式;
将特高压直流故障下的功率波动分解为若干个脉冲函数和阶跃函数的叠加;
计算各脉冲函数下的交流联络线功率振荡响应以及各阶跃函数下的交流联络线功率振荡响应;
将所有交流联络线功率振荡响应叠加即为特高压直流故障下的大区交流联络线功率振荡响应。
简化交流系统模型包括通过交流联络线连接的送端系统A和受端系统B,特高压直流往送端系统A输送功率,送端系统A发电机组和受端系统B发电机为两等值发电机。
脉冲功率冲击下交流联络线功率振荡响应公式为,
Figure BDA0001421803690000021
其中,ΔPtie(t)为t时刻的交流联络线功率振荡响应,TJ2为受端系统B的总惯性常数,TJ1为送端系统A的总惯性常数,
Figure BDA0001421803690000022
ω0=2πf0,f0为基准频率,KS为同步转矩系数,
Figure BDA0001421803690000023
D1TJ2=D2TJ1=D,D1为送端系统A发电机阻尼系数,D2为受端系统B发电机阻尼系数,
Figure BDA0001421803690000024
F为功率冲击的幅值。
脉冲功率冲击下交流联络线功率振荡响应公式推导过程为,
假设交流系统在受扰期间发电机机械功率和系统负荷保持恒定,发电机转子运动方程在工作点附近线性化,
Figure BDA0001421803690000031
其中,Δδ1为送端系统A发电机转子角变化,ΔPm1为送端系统A原动机机械功率变化,ΔPL为功率冲击,ΔPtie为交流联络线功率振荡响应,Δδ2为受端系统B发电机转子角变化,ΔPm2为受端系统原动机机械功率变化;
对交流联络线功率线性化有,
Figure BDA0001421803690000032
其中,E1',E2'分别为送端系统A和受端系统B发电机内电势的幅值,δ120表示稳态点的等值发电机功角差,X=X1+X2+XL,X1为送端系统A发电机等值电抗,X2为受端系统B发电机等值电抗,XL为交流联络线电抗,Δδ12为两个等值发电机功角差变化;
忽略ΔPm1和ΔPm2,有,
Figure BDA0001421803690000033
由线性均匀性原理得,
Figure BDA0001421803690000034
Figure BDA0001421803690000035
功率冲击ΔPL为脉冲功率冲击,ΔPL=-Fδ(t),δ(t)为脉冲函数,令x=ΔPtie,有,
Figure BDA0001421803690000041
系统为欠阻尼系统时,系统输出的拉式变换式为,
Figure BDA0001421803690000042
分解为部分分式为,
Figure BDA0001421803690000043
系统输出响应为,
Figure BDA0001421803690000044
Figure BDA0001421803690000045
阶跃功率冲击下交流联络线功率振荡响应公式为,
Figure BDA0001421803690000046
其中,ΔPtie(t)为t时刻的交流联络线功率振荡响应,TJ2为受端系统B的总惯性常数,TJ1为送端系统A的总惯性常数,
Figure BDA0001421803690000047
ω0=2πf0,f0为基准频率,KS为同步转矩系数,
Figure BDA0001421803690000048
D1TJ2=D2TJ1=D,D1为送端系统A发电机阻尼系数,D2为受端系统B发电机阻尼系数,
Figure BDA0001421803690000049
阶跃功率冲击下交流联络线功率振荡响应公式推导过程为,
在受扰期间发电机机械功率和系统负荷保持恒定,功率冲击ΔPL为阶跃功率冲击,ΔPL=-Fε(t),ε(t)为阶跃函数,令x=ΔPtie,有,
Figure BDA0001421803690000051
系统为欠阻尼系统时,系统输出的拉式变换式为,
Figure BDA0001421803690000052
分解为部分分式为,
Figure BDA0001421803690000053
系统输出响应为,
Figure BDA0001421803690000054
Figure BDA0001421803690000055
本发明所达到的有益效果:本发明针对特高压直流发生故障后对大区交流联络线造成的功率振荡问题,通过在简化交流系统中的理论分析和公式推导,给出了对功率振荡响应的计算方法,能够比较准确地得出大区交流联络线的功率振荡响应,同时,计算方法简单易用,适合于如今对交直流混联电网的交流联络线功率振荡研究。
附图说明
图1为本发明的流程图;
图2为简化交流系统的结构图;
图3为有功功率实测录波曲线;
图4为脉冲函数曲线;
图5为第一个阶跃函数曲线;
图6为第二个阶跃函数曲线;
图7为第三个阶跃函数曲线;
图8为功率振荡曲线。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示,大区交流联络线功率振荡响应的计算方法,包括以下步骤:
步骤1,构建简化交流系统模型。
如图2所示,简化交流系统模型包括通过交流联络线连接的送端系统A和受端系统B,送端系统A发电机组和受端系统B发电机为两等值发电机,特高压直流往送端系统A输送功率,其可以等效为一个功率为负的负荷,因此在图中将特高压直流与送端系统A的总负荷合并为PL1,特高压直流发生故障时产生的功率冲击,定义为ΔPL,PL2为受端系统B的总负荷,Pm1和Pm2分别为送端系统A和受端系统B原动机机械功率,TJ1和TJ2分别为送端系统A和受端系统B的总惯性常数,Ptie为交流联络线功率,X1为送端系统A发电机等值电抗,X2为受端系统B发电机等值电抗,XL为交流联络线电抗,两等值发电机内电势分别用E1'∠δ1和E2'∠δ2表示,E1',E2'分别为送端系统A和受端系统B发电机内电势的幅值,δ1为送端系统A发电机转子角,δ2为受端系统B发电机转子角,。
步骤2,基于简化交流系统模型,推导出脉冲功率冲击下交流联络线功率振荡响应公式和阶跃功率冲击下交流联络线功率振荡响应公式。
脉冲功率冲击下交流联络线功率振荡响应公式推导过程如下:
假设交流系统在受扰期间发电机机械功率和系统负荷保持恒定,发电机转子运动方程在工作点附近线性化,
Figure BDA0001421803690000071
其中,Δδ1为送端系统A发电机转子角变化,ΔPm1为送端系统A原动机机械功率变化,ΔPtie为交流联络线功率振荡响应,Δδ2为受端系统B发电机转子角变化,ΔPm2为受端系统原动机机械功率变化,ω0=2πf0,f0为基准频率,D1为送端系统A发电机阻尼系数,D2为受端系统B发电机阻尼系数;
对交流联络线功率线性化有,
Figure BDA0001421803690000072
其中,δ120表示稳态点的等值发电机功角差,X=X1+X2+XL,Δδ12为两个等值发电机功角差变化,δ12为两个等值发电机的功角差,KS为同步转矩系数,
Figure BDA0001421803690000073
忽略ΔPm1和ΔPm2,有,
Figure BDA0001421803690000081
由线性均匀性原理,假设D1TJ2=D2TJ1=D,
Figure BDA0001421803690000082
得,
Figure BDA0001421803690000083
Figure BDA0001421803690000084
整理上式得,
Figure BDA0001421803690000085
功率冲击ΔPL为脉冲功率冲击,ΔPL=-Fδ(t),F为功率冲击的幅值,δ(t)为脉冲函数,令x=ΔPtie,有,
Figure BDA0001421803690000086
系统为欠阻尼系统时,系统输出的拉式变换式为,
Figure BDA0001421803690000087
分解为部分分式为,
Figure BDA0001421803690000088
系统输出响应为,
Figure BDA0001421803690000089
即,
Figure BDA0001421803690000091
其中,ΔPtie(t)为t时刻的交流联络线功率振荡响应,
Figure BDA0001421803690000092
单个脉冲功率冲击下交流联络线功率振荡时域响应为一个自由振荡过程,振荡幅值受脉冲幅值、脉冲宽度、子系统惯性时间常数等因素影响,振荡频率由系统固有振荡频率决定。
阶跃功率冲击下交流联络线功率振荡响应公式推导过程如下:
在受扰期间发电机机械功率和系统负荷保持恒定,功率冲击ΔPL为阶跃功率冲击,ΔPL=-Fε(t),ε(t)为阶跃函数,令x=ΔPtie
代入
Figure BDA0001421803690000093
有,
Figure BDA0001421803690000094
系统为欠阻尼系统时,系统输出的拉式变换式为,
Figure BDA0001421803690000095
分解为部分分式为,
Figure BDA0001421803690000096
系统输出响应为,
Figure BDA0001421803690000097
即,
Figure BDA0001421803690000101
其中,
Figure BDA0001421803690000102
阶跃功率冲击下交流联络线功率振荡时域响应,等式右边第一项为一个与阶跃幅值有关的稳态分量,等式右边第二项可分解为由初始状态引起的自由振荡分量和由阶跃扰动引起的伴随自由振荡分量。
步骤3,将特高压直流故障下的功率波动分解为若干个脉冲函数和阶跃函数的叠加。
步骤4,计算各脉冲函数下的交流联络线功率振荡响应以及各阶跃函数下的交流联络线功率振荡响应;
步骤5,将所有交流联络线功率振荡响应叠加即为特高压直流故障下的大区交流联络线功率振荡响应。
以华北、华中电网为例,华北、华中电网通过1000kV交流联络线相联,华中电网有特高压直流工程天中直流馈入,直流落点位于河南电网。华北—华中区间惯性时间常数比约为1.1:1,分析计算在天中直流故障下华中-华北交流联络线功率振荡响应。
计算过程如下:
1、将特高压直流故障下天中直流的有功功率波动分解为多个脉冲函数与阶跃函数的叠加。
图3的有功功率冲击可以分解为图2所示的1个脉冲函数以及图5~图7所示的3个阶跃函数的叠加,对应的表达式为:
脉冲功率冲击ΔPL1=-98δ(t-t0);
第一个阶跃功率冲击ΔPL2=-1730ε(t-t0);
第二个阶跃功率冲击ΔPL3=270ε(t-t0-0.18);
第三个阶跃功率冲击ΔPL4=-255ε(t-t0-3.18);
其中,t0为特高压直流故障开始时间。
2、将ΔPL1代入脉冲功率冲击下交流联络线功率振荡响应公式,得到功率振荡响应ΔPtie1(t);
将ΔPL2、ΔPL3和ΔPL4代入阶跃功率冲击下交流联络线功率振荡响应公式,分别得到功率振荡响应ΔPtie2(t)、ΔPtie3(t)和ΔPtie4(t);
3、将所有交流联络线功率振荡响应叠加ΔPtie1(t)+ΔPtie2(t)+ΔPtie3(t)+ΔPtie4(t),即为天中特高压直流故障下华中-华北交流联络线功率振荡响应,如图8所示。
针对特高压直流发生故障后对大区交流联络线造成的功率振荡问题,通过在简化交流系统中的理论分析和公式推导,给出了对功率振荡响应的计算方法,该方法能够比较准确地得出大区交流联络线的功率振荡响应,同时,计算方法简单易用,适合于如今对交直流混联电网的交流联络线功率振荡研究。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。

Claims (4)

1.大区交流联络线功率振荡响应的计算方法,其特征在于:包括,
构建简化交流系统模型;简化交流系统模型包括通过交流联络线连接的送端系统A和受端系统B,特高压直流往送端系统A输送功率,送端系统A发电机组和受端系统B发电机为两等值发电机;
基于简化交流系统模型,推导出脉冲功率冲击下交流联络线功率振荡响应公式和阶跃功率冲击下交流联络线功率振荡响应公式;
脉冲功率冲击下交流联络线功率振荡响应公式为,
Figure FDA0002377827150000011
其中,ΔPtie(t)为t时刻的交流联络线功率振荡响应,TJ2为受端系统B的总惯性常数,TJ1为送端系统A的总惯性常数,
Figure FDA0002377827150000012
ω0=2πf0,f0为基准频率,KS为同步转矩系数,
Figure FDA0002377827150000013
D1TJ2=D2TJ1=D,D1为送端系统A发电机阻尼系数,D2为受端系统B发电机阻尼系数,
Figure FDA0002377827150000014
F为功率冲击的幅值;
将特高压直流故障下的功率波动分解为若干个脉冲函数和阶跃函数的叠加;
计算各脉冲函数下的交流联络线功率振荡响应以及各阶跃函数下的交流联络线功率振荡响应;
将所有交流联络线功率振荡响应叠加即为特高压直流故障下的大区交流联络线功率振荡响应。
2.根据权利要求1所述的大区交流联络线功率振荡响应的计算方法,其特征在于:脉冲功率冲击下交流联络线功率振荡响应公式推导过程为,
假设交流系统在受扰期间发电机机械功率和系统负荷保持恒定,发电机转子运动方程在工作点附近线性化,
Figure FDA0002377827150000021
其中,Δδ1为送端系统A发电机转子角变化,ΔPm1为送端系统A原动机机械功率变化,ΔPL为功率冲击,ΔPtie为交流联络线功率振荡响应,Δδ2为受端系统B发电机转子角变化,ΔPm2为受端系统原动机机械功率变化;
对交流联络线功率线性化有,
Figure FDA0002377827150000022
其中,E1',E2'分别为送端系统A和受端系统B发电机内电势的幅值,δ120表示稳态点的等值发电机功角差,X=X1+X2+XL,X1为送端系统A发电机等值电抗,X2为受端系统B发电机等值电抗,XL为交流联络线电抗,Δδ12为两个等值发电机功角差变化;
忽略ΔPm1和ΔPm2,有,
Figure FDA0002377827150000023
由线性均匀性原理得,
Figure FDA0002377827150000024
Figure FDA0002377827150000031
功率冲击ΔPL为脉冲功率冲击,ΔPL=-Fδ(t),δ(t)为脉冲函数,令x=ΔPtie,有,
Figure FDA0002377827150000032
系统为欠阻尼系统时,系统输出的拉式变换式为,
Figure FDA0002377827150000033
分解为部分分式为,
Figure FDA0002377827150000034
系统输出响应为,
Figure FDA0002377827150000035
Figure FDA0002377827150000036
3.根据权利要求1所述的大区交流联络线功率振荡响应的计算方法,其特征在于:阶跃功率冲击下交流联络线功率振荡响应公式为,
Figure FDA0002377827150000037
其中,ΔPtie(t)为t时刻的交流联络线功率振荡响应,TJ2为受端系统B的总惯性常数,TJ1为送端系统A的总惯性常数,
Figure FDA0002377827150000038
ω0=2πf0,f0为基准频率,KS为同步转矩系数,
Figure FDA0002377827150000041
D1TJ2=D2TJ1=D,D1为送端系统A发电机阻尼系数,D2为受端系统B发电机阻尼系数,
Figure FDA0002377827150000042
F为功率冲击的幅值。
4.根据权利要求3所述的大区交流联络线功率振荡响应的计算方法,其特征在于:阶跃功率冲击下交流联络线功率振荡响应公式推导过程为,
在受扰期间发电机机械功率和系统负荷保持恒定,功率冲击ΔPL为阶跃功率冲击,ΔPL=-Fε(t),ε(t)为阶跃函数,令x=ΔPtie,有,
Figure FDA0002377827150000043
系统为欠阻尼系统时,系统输出的拉式变换式为,
Figure FDA0002377827150000044
分解为部分分式为,
Figure FDA0002377827150000045
系统输出响应为,
Figure FDA0002377827150000046
Figure FDA0002377827150000047
CN201710894105.XA 2017-09-28 2017-09-28 大区交流联络线功率振荡响应的计算方法 Active CN107611999B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710894105.XA CN107611999B (zh) 2017-09-28 2017-09-28 大区交流联络线功率振荡响应的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710894105.XA CN107611999B (zh) 2017-09-28 2017-09-28 大区交流联络线功率振荡响应的计算方法

Publications (2)

Publication Number Publication Date
CN107611999A CN107611999A (zh) 2018-01-19
CN107611999B true CN107611999B (zh) 2020-05-22

Family

ID=61059032

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710894105.XA Active CN107611999B (zh) 2017-09-28 2017-09-28 大区交流联络线功率振荡响应的计算方法

Country Status (1)

Country Link
CN (1) CN107611999B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101847872A (zh) * 2010-06-01 2010-09-29 中国电力科学研究院 两大区互联电力系统交流联络线功率波动峰值的计算方法
CN103928939A (zh) * 2014-04-14 2014-07-16 国家电网公司 直流换相失败冲击下交流联络线功率波动峰值的计算方法
CN104122819A (zh) * 2014-07-18 2014-10-29 东北电力大学 基于用户习惯的家庭智能用电方法及其系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101847872A (zh) * 2010-06-01 2010-09-29 中国电力科学研究院 两大区互联电力系统交流联络线功率波动峰值的计算方法
CN103928939A (zh) * 2014-04-14 2014-07-16 国家电网公司 直流换相失败冲击下交流联络线功率波动峰值的计算方法
CN104122819A (zh) * 2014-07-18 2014-10-29 东北电力大学 基于用户习惯的家庭智能用电方法及其系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"电力系统‘广域同步相量测量’中几个关键问题的研究";杨贵玉;《中国博士学位论文全文数据库 工程科技II辑》;20051015(第6期);C042-25 *
"送端电网稳定性若干问题及控制策略研究";吕清洁;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20160715(第7期);C042-349 *

Also Published As

Publication number Publication date
CN107611999A (zh) 2018-01-19

Similar Documents

Publication Publication Date Title
CN103559347B (zh) 一种大规模交直流电力系统电磁暂态仿真模型的构建方法
CN104362668B (zh) 一种电压不平衡/谐波畸变时双馈风力发电机的控制方法
CN106159974A (zh) 一种输配协调的分布式无功电压优化方法
CN103269088B (zh) 一种基于电磁暂态算法的双馈型风力发电机组低电压穿越控制方法
CN104993767B (zh) 考虑谐波耦合的多相电机缺相故障容错运行电流重构方法
CN103647467B (zh) 一种基于粒子群算法的不平衡电网下vsc多目标优化直接功率控制方法
CN102592003A (zh) 一种用于机电与电磁暂态混合仿真的数据交换方法
CN102281029B (zh) 一种无轴承同步磁阻电机悬浮系统构造方法
CN105226679A (zh) 一种在线动态抑制电网低频振荡和次同步振荡的控制方法
CN105281350A (zh) 微电网频率控制方法和系统
CN104377697A (zh) 模拟同步发电机动态特性的电压源换流站的控制方法
CN102769297A (zh) 一种基于全控型逆变器的次同步振荡抑制方法
CN105790306B (zh) 基于逆变器电流单传感器的lcl型并网逆变器控制方法
CN105811470A (zh) 一种电网不平衡条件下虚拟同步机控制系统及控制方法
CN112186758B (zh) 一种自适应捕捉频点的直驱风电机组次同步振荡抑制方法
CN103326399A (zh) 一种不平衡及谐波电网下的并网逆变器控制方法
CN102832643A (zh) 一种基于逆系统的三相光伏并网逆变器的控制方法
CN103259472A (zh) 大型发电机组次同步振荡的反振荡自适应抑制系统
CN103346728A (zh) 双馈风力发电机定转子电流传感器故障判断方法
CN112165106A (zh) 基于双二阶广义积分器的dfig次同步振荡抑制方法
CN107134962A (zh) 一种三级式同步电机转子位置估算方法
CN103944186B (zh) 一种三相光伏并网逆变器控制装置
CN107611999B (zh) 大区交流联络线功率振荡响应的计算方法
CN101719678A (zh) 双馈型风力发电机空载并网建模及实验的方法
CN106655270B (zh) 抑制三相lcl型并网变换器整流启动冲击电流的控制方法

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