CN101699693B - 铁路牵引冲击负荷对电网电压波动计算方法 - Google Patents

铁路牵引冲击负荷对电网电压波动计算方法 Download PDF

Info

Publication number
CN101699693B
CN101699693B CN2009102135962A CN200910213596A CN101699693B CN 101699693 B CN101699693 B CN 101699693B CN 2009102135962 A CN2009102135962 A CN 2009102135962A CN 200910213596 A CN200910213596 A CN 200910213596A CN 101699693 B CN101699693 B CN 101699693B
Authority
CN
China
Prior art keywords
traction
load
calculate
locomotive
power
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
CN2009102135962A
Other languages
English (en)
Other versions
CN101699693A (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.)
JIANGSU KENENG ELECTRIC ENGINEERING CONSULTATION Co Ltd
Original Assignee
JIANGSU KENENG ELECTRIC ENGINEERING CONSULTATION Co Ltd
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 JIANGSU KENENG ELECTRIC ENGINEERING CONSULTATION Co Ltd filed Critical JIANGSU KENENG ELECTRIC ENGINEERING CONSULTATION Co Ltd
Priority to CN2009102135962A priority Critical patent/CN101699693B/zh
Publication of CN101699693A publication Critical patent/CN101699693A/zh
Application granted granted Critical
Publication of CN101699693B publication Critical patent/CN101699693B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

一种铁路牵引冲击负荷对电网电压波动计算方法,本方法是应用电力系统分析综合程序PSASP进行,包括步骤:1)确定机车启动前牵引站高压侧负荷值并输入电网模型中,作潮流计算,以此潮流计算结果作为暂态稳定计算用的基础潮流;2)分别计算牵引变两臂上两列机车启动时,高压侧无功冲击负荷的大小;3)应用PSASP程序,在牵引站高压侧母线模拟无功冲击,计算PPC点的电压波动.。本发明对于减小牵引负荷对电网冲击影响、提高电力系统电能质量、电网运行安全稳定具有重要意义。

Description

铁路牵引冲击负荷对电网电压波动计算方法
技术领域
本发明属于电力工程技术领域,具体是一种铁路牵引冲击负荷对电网电压波动计算方法。
背景技术
电气化铁路的牵引动力是电力机车,机车本身不带能源,所需能源由电力牵引供电系统提供。牵引供电系统主要是指牵引变电所和接触网两大部分。变电所设在铁道附近,它将从发电厂经高压输电线送来的电流,送到铁路上空的接触网上。接触网是向电力机车直接输送电能的设备。沿着铁路线的两旁,架设着一排支柱,上面悬挂着金属线,即为接触网。电力机车利用车顶的受电弓从接触网获得电能,牵引列车运行。牵引供电制式主要是交流制。交流制是将高压、三相电力在变电所降压和变成单相后,向接触网供交流电。目前我国使用单相25千伏50赫兹的工频交流电供应,这样的供电形式比较经济。
目前我国在建的电气化高速铁路(如京沪高铁、沪宁城际铁路、宁杭客运专线铁路)的列车追踪时间间隔按3min设计,按每3min内两辆机车(上、下行线各一辆)各经过无电区一次考虑,则设备需要投(启动)、切(停止)4次(此处投/切是指合上/断开牵引机车变压器的开关),则每小时设备投切次数达到80次,如此频繁的设备投切次数对系统造成的电压波动和闪变是不可忽视的。
目前在三相对称电力系统中计算冲击负荷在PCC点(即公共连接点)引起的电压波动常用的方法是,应用电力系统分析综合程序PSASP,在程序中建立负荷模型(三相对称负荷),通过在母线上模拟一个无功冲击负荷ΔQk及其变化时间Δt,应用暂态稳定计算模块计算PCC点的电压波动。
但是,由于铁路牵引负荷为单相供电负荷,当电力系统向其供电时,将造成电力系统三相负荷不平衡;此外,无功冲击负荷还会在牵引变压器内部产生无功损耗的增量,在计算电压波动时,该增量也应考虑在内。上述两点因素在传统的计算方法均未考虑,因此传统的模拟三相对称系统的方法不能直接应用于铁路牵引负荷对电压冲击影响的模拟计算。
发明内容
为了解决现有技术中存在的上述问题,本发明针对铁路牵引负荷在电力系统造成三相不平衡的特点,提出的一种计算电压波动的方法,具体技术方案如下:
一种铁路牵引冲击负荷对电网电压波动计算方法,本方法是应用电力系统分析综合程序PSASP进行,包括步骤:
1)确定机车启动前牵引站高压侧负荷值并输入电网模型中,作潮流计算,以此潮流计算结果作为暂态稳定计算用的基础潮流;
2)分别计算牵引变两臂上两列机车启动时,高压侧无功冲击负荷的大小:
2.1)根据牵引变压器两供电臂瞬时最大负荷,以及机车功率ΔP0+jΔQ0,计算机车启动前两供电臂的有功和无功负荷P1+jQ1、P2+jQ2
2.2)根据牵引变压器接线形式(由两个单相变压器接成V/v结线变压器)以及接入相序(相序可取-AB、BC或者-BC、CA或者-CA、AB,具体相序与牵引站所处的电力系统有关,对下面的计算结果没有影响),可计算出机车启动前牵引站高压侧三相有功负荷,分别为
由于电力系统为三相对称系统,为校核电力系统对牵引站的供电能力,应用PSASP程序做潮流计算时,应该取最大一相负荷的3倍为暂态稳定计算基础潮流中的有功负荷P0,即 P 0 = 3 × P 1 2 + P 2 2 + P 1 P 2 / 3 = 3 ( P 1 2 + P 2 2 + P 1 P 2 ) ;
2.3)由Q1、Q2计算机车启动前牵引变压器的无功损耗ΔQT1、ΔQT2,则Q′1=Q1+ΔQT1、Q′2=Q2+ΔQT2,按步骤B的方法可算出暂态稳定计算基础潮流中的无功负荷 Q 0 = 3 ( Q 1 ′ 2 + Q 2 ′ 2 + Q 1 ′ Q 2 ′ ) ;
其中Q′1、Q′2为低压侧两臂计及牵引变压器无功损耗的无功负荷。
3)应用PSASP程序,在牵引站高压侧母线模拟无功冲击,计算PPC点的电压波动:
3.1)在PSASP程序中设定牵引变电站高压侧母线(实际牵引站高压侧无母线,此处假设是方便潮流计算),并在母线上增加一个负荷,大小为P0+jQ0
3.2)先计算一臂上的机车启动时产生的无功冲击负荷ΔQ′=αΔQ0,α为励磁涌流倍数;再计算ΔQ在牵引变压器中产生的无功损耗ΔQ′T1,则无功损耗增量δQ1=ΔQ′T1-ΔQT1;最后得到一臂上的机车启动时,牵引站低压侧的无功冲击负荷ΔQ1=(ΔQ+δQ1);
同样方法可计算得到二臂上机车启动时牵引站低压侧的无功冲击负荷ΔQ2=ΔQ+δQ2
3.3)二个臂上的机车都启动后,牵引站高压侧无功冲击最大一相的值为
Figure G2009102135962D00031
取其3倍作为暂态稳定计算中三相系统的无功冲击值,则 ΔQ = 3 ( Δ Q 1 2 + Δ Q 2 2 + Δ Q 1 Δ Q 2 ) ;
3.4)在PSASP程序中利用暂态稳定计算模块(该模块是PSASP程序自带计算模块),并在基础潮流计算结果的基础上,在牵引站高压侧母线上设定节点(母线)扰动,扰动大小为ΔQ,扰动时间取10ms(此时间是机车启动瞬间产生的冲击扰动时间间隔,是不可缺少的必要条件,扰动时间的取值参考《电气化铁路负荷最大电流突变量的估算方法》(详见《电力系统自动化》2008年第9期、P79~81),观察公共接入点的电压变化大小,电压变化量与额定电压的比值即为电压波动。
本发明采用三相系统计算程序PSASP来模拟单相无功冲击负荷,通过计算牵引站高压侧无功冲击最大一相在PCC点引起的电压波动,根据计算的电压波动确定的低压侧无功补偿容量更为科学合理,对于减小牵引负荷对电网冲击影响、提高电力系统电能质量、电网运行安全稳定具有重要意义。
附图说明
图1为本发明实施例的电网示意图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步说明:
一种铁路牵引冲击负荷对电网电压波动计算方法,本方法是应用电力系统分析综合程序PSASP进行,步骤是,先确定机车启动前牵引站高压侧负荷值并输入电网模型中,作潮流计算,以此潮流计算结果作为暂态稳定计算用的基础潮流(本步骤在现有技术中的常用方法,在此就不进一步加以限制了);
再按照下述步骤进行:
A、根据牵引变压器两供电臂瞬时最大负荷,以及机车功率ΔP0+jΔQ0,计算机车启动前两供电臂的负荷P1+jQ1、P2+jQ2
B、取牵引站变压器为V/v接线变压器,接入系统相序为-AB、BC,计算机车启动前牵引站高压侧三相有功负荷,分别为
Figure G2009102135962D00041
取最大一相的3倍做为暂态稳定计算基础潮流中的有功负荷,即 P 0 = 3 × P 1 2 + P 2 2 + P 1 P 2 / 3 = 3 ( P 1 2 + P 2 2 + P 1 P 2 ) ;
C、由Q1、Q2计算机车启动前牵引变的无功损耗ΔQT1、ΔQT2,则Q′1=Q1+ΔQT1、Q′2=Q2+ΔQT2,按步骤B的方法可算出暂态稳定计算基础潮流中的无功负荷 Q 0 = 3 ( Q 1 ′ 2 + Q 2 ′ 2 + Q 1 ′ Q 2 ′ ) ;
D、在PSASP程序中建立牵引负荷模型,基础潮流计算时牵引负荷为P0+jQ0
所述模型是,在PSASP程序中设定牵引变电站高压侧母线(实际牵引站高压侧无母线,此处假设是方便潮流计算),并在母线上增加一个负荷,大小为P0+jQ0
E、计算一臂机车启动时产生的无功冲击负荷ΔQ′=αΔQ0,式中α为励磁涌流倍数;
F、计算ΔQ在牵引变中产生的无功损耗ΔQ′T1,则无功损耗增量δQ1=ΔQ′T1-ΔQT1
G、一臂机车启动时,牵引站低压侧的无功冲击负荷ΔQ1=(ΔQ+δQ1);
H、同样方法可计算二臂机车启动时牵引站低压侧的无功冲击负荷ΔQ2=ΔQ+δQ2
I、二臂机车启动后,牵引站高压侧无功冲击最大一相的值为
Figure G2009102135962D00044
取其3倍做为暂态稳定计算中三相系统的无功冲击值,则 ΔQ = 3 ( Δ Q 1 2 + Δ Q 2 2 + Δ Q 1 Δ Q 2 ) ;
J、在PSASP程序中利用暂态稳定模块(本模块是PSASP程序中自带工具模块)做暂态稳定计算,即可算得公共接入点的电压波动;
例如,某牵引站负荷及机车功率等参数如表1:
表1
  项目   参数
  牵引变压器容量(V/v接线)   (40+40)MVA
  牵引变压器短路阻抗参数   10.5%
  瞬时最大负荷(一臂、二臂)   80.1MVA、86.6MVA
  负荷功率因数   0.97
  机车有功功率   17.6MW
  机车功率因数   0.98
步骤如下:
(1)上表中两臂的瞬时最大负荷也可表示为P1臂+jQ1臂=77.7+j19.47(MVA)、P2臂+jQ2臂=84.0+j21.05(MVA)
其中:P1臂=80.1×0.97=77.7
P2臂=86.6×0.97=84.0
Figure G2009102135962D00052
(2)根据机车有功功率和机车功率因数(0.98)可知机车功率为
ΔP0+jΔQ0=17.6+j3.57(MVA);
其中:ΔP0=17.6
ΔQ0=17.6×tan(arccos0.98)=3.57
(3)机车启动前两臂负荷分别为
P1+jQ1=77.1+j19.47-17.6-j3.57=60.1+j15.9(MVA)
P2+jQ2=84.0+j21.05-17.6-j3.57=66.4+j17.48(MVA)
(4)计算机车启动前两个单相变压器的无功损耗:
Figure G2009102135962D00053
Figure G2009102135962D00054
(5)计算机车启动前两个单相变压器总的无功损耗:
Q′1=Q1+ΔQT1=15.9+10.14=26.04(Mvar)
Q′2=Q2+ΔQT2=17.48+12.38=29.86(Mvar)
(6)计算牵引站高压侧三相负荷(潮流计算用):
P 0 = 3 ( P 1 2 + P 2 2 + P 1 P 2 ) = 3 × ( 60 . 1 2 + 66.4 2 + 60.1 × 66.4 ) = 189.83 ( MW )
Q 0 = 3 ( Q 1 2 + Q 2 2 + Q 1 Q 2 ) = 3 × ( 26.0 4 2 + 29.86 2 + 26.04 × 29.86 ) = 83.92 ( Mvar )
(7)测量得:一臂机车启动时牵引站低压侧无功冲击负荷ΔQ1=55.63Mvar;
(8)测量得:二臂机车启动时牵引站低压侧无功冲击负荷ΔQ2=56.03Mvar;
(9)计算两臂上的机车同时启动后,折算至高压侧的单相最大无功冲击负荷:
ΔQ = 3 × ( ΔQ 1 2 + Δ Q 2 2 + Δ Q 1 Δ Q 2 ) = 167.49 Mvar
(10)在PSASP程序中利用暂态稳定计算模块,在基础潮流计算结果的基础上,在牵引站高压侧母线上设定节点(母线)扰动,扰动大小为ΔQ,扰动时间取10ms,观察公共接入点的电压变化大小,电压变化量与额定电压的比值即为电压波动;
根据计算结果,公共接入点的电压波动为1.97%。根据国家标准(GB12326-2008《电能质量电压波动和闪变》)的相关规定,公共接入点电压波动限值为1%,则需要的改善率 k = 1.97 - 1 1.97 × 100 % = 49.2 % .
本发明采用三相系统计算程序模拟单相无功冲击负荷在公共接入点产生的电压波动,计算过程中考虑了机车启动前牵引变压器的无功损耗,以及启动瞬间无功冲击负荷在牵引变压器中产生的无功损耗的增量。该计算方法针对铁路牵引负荷三相不对称的特点,同时计及牵引变压器无功损耗的增量,计算结果准确,可有效地针对电铁负荷在公共接入点电压波动超标问题进行相应的治理。

Claims (1)

1.一种铁路牵引冲击负荷对电网电压波动计算方法,本方法是应用电力系统分析综合程序PSASP进行,包括步骤:
1)确定机车启动前牵引站高压侧负荷值并输入电网模型中,作潮流计算,以此潮流计算结果作为暂态稳定计算用的基础潮流;
其特征是还包括步骤:
2)分别计算牵引变两臂上两列机车启动时,高压侧无功冲击负荷的大小,步骤包括:
2.1)根据牵引变压器两供电臂瞬时最大负荷,以及机车功率ΔP0+jΔQ0,计算机车启动前两供电臂的有功和无功负荷P1+jQ1和P2+jQ2
2.2)根据牵引变压器接线形式即由两个单相变压器接成V/v结线变压器,以及接入相序,算得机车启动前牵引站高压侧三相有功负荷,分别为
Figure FSB00000514591400012
Figure FSB00000514591400013
由于电力系统为三相对称系统,为校核电力系统对牵引站的供电能力,应用PSASP程序做潮流计算时,取最大一相负荷的3倍为暂态稳定计算基础潮流中的有功负荷P0,即
2.3)由Q1、Q2计算机车启动前牵引变压器的无功损耗ΔQT1、ΔQT2,则Q′1=Q1+ΔQT1、Q′2=Q2+ΔQT2,按步骤2.2)的方法算得暂态稳定计算基础潮流中的无功负荷
Figure FSB00000514591400015
其中,Q′1、Q′2为低压侧两臂计及牵引变压器无功损耗的无功负荷;
3)应用PSASP程序,在牵引站高压侧母线模拟无功冲击,计算公共接入点PPC的电压波动:
3.1)在PSASP程序中设定牵引变电站高压侧母线,并在母线上增加一个负荷,大小为P0+jQ0;此处,实际牵引站高压侧无母线,此处假设是方便潮流计算;
3.2)先计算一臂上的机车启动时产生的无功冲击负荷ΔQ′=αΔQ0,α为励磁涌流倍数;再计算ΔQ在牵引变压器中产生的无功损耗ΔQ′T1,则无功损耗增量δQ1=ΔQ′T1-ΔQT1;最后得到一臂上的机车启动时,牵引站低压侧的无功冲击负荷ΔQ1=(ΔQ+δQ1);
同样方法可计算得到二臂上机车启动时牵引站低压侧的无功冲击负荷ΔQ2=ΔQ+δQ2
3.3)两个臂上的机车都启动后,牵引站高压侧无功冲击最大一相的值为取其3倍作为暂态稳定计算中三相系统的无功冲击值,则 ΔQ = 3 ( Δ Q 1 2 + Δ Q 2 2 + Δ Q 1 Δ Q 2 ) ;
3.4)在PSASP程序中利用暂态稳定计算模块,并在基础潮流计算结果的基础上,在牵引站高压侧母线上设定节点即母线扰动,扰动大小为ΔQ,扰动时间取10ms,观察公共接入点PPC的电压变化大小,电压变化量与额定电压的比值即为公共接入点PPC的电压波动。
CN2009102135962A 2009-11-06 2009-11-06 铁路牵引冲击负荷对电网电压波动计算方法 Active CN101699693B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009102135962A CN101699693B (zh) 2009-11-06 2009-11-06 铁路牵引冲击负荷对电网电压波动计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009102135962A CN101699693B (zh) 2009-11-06 2009-11-06 铁路牵引冲击负荷对电网电压波动计算方法

Publications (2)

Publication Number Publication Date
CN101699693A CN101699693A (zh) 2010-04-28
CN101699693B true CN101699693B (zh) 2011-09-07

Family

ID=42148144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009102135962A Active CN101699693B (zh) 2009-11-06 2009-11-06 铁路牵引冲击负荷对电网电压波动计算方法

Country Status (1)

Country Link
CN (1) CN101699693B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101917000B (zh) * 2010-07-14 2012-12-26 北京交通大学 电气化铁道牵引变电所负荷的计算方法
CN102570460A (zh) * 2010-12-14 2012-07-11 中国电力科学研究院 一种计算电气化铁路的谐波和负序电流的方法
CN102184296B (zh) * 2011-05-13 2012-11-28 江西省电力科学研究院 基于实测数据的电气化铁路冲击负荷建模方法
CN102243477B (zh) * 2011-06-17 2015-08-19 河南省电力公司信阳供电公司 基于等值故障变换的牵引负荷建模方法
CN102904254B (zh) * 2012-10-13 2015-03-04 山东电力集团公司电力科学研究院 电气化铁路注入电力系统谐波的仿真方法
CN103078347B (zh) * 2013-01-12 2015-03-11 华南理工大学 降低电压越限风险的冲击负荷消纳方法
CN103078329B (zh) * 2013-02-05 2014-08-13 江苏科能电力工程咨询有限公司 海上风电场长距离220kV海缆送出无功补偿分析方法
CN105244878A (zh) * 2015-11-13 2016-01-13 贵州电网有限责任公司电力科学研究院 电力系统三相不平衡时的电力系统损耗增量的分配方法
CN106874691A (zh) * 2017-03-09 2017-06-20 国网黑龙江省电力有限公司电力科学研究院 双线电铁中牵引变电站的负荷电压波动频度计算方法
CN108287936B (zh) * 2017-12-18 2023-05-26 国网陕西省电力公司经济技术研究院 一种牵引供电系统的电磁暂态模型
CN113011004B (zh) * 2021-02-24 2024-03-15 南京地铁建设有限责任公司 一种安全校核方法、装置及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101340098A (zh) * 2008-01-03 2009-01-07 中国南车集团株洲电力机车研究所 一种动态无功补偿方法
CN101453171A (zh) * 2008-09-12 2009-06-10 清华大学 基于变压器串联多重化和链式结构的统一电能质量控制器

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101340098A (zh) * 2008-01-03 2009-01-07 中国南车集团株洲电力机车研究所 一种动态无功补偿方法
CN101453171A (zh) * 2008-09-12 2009-06-10 清华大学 基于变压器串联多重化和链式结构的统一电能质量控制器

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李庚银等.含牵引负荷的电力系统三相不对称谐波潮流计算.《电力系统自动化》.1999,第23卷(第17期), *
胡煜等.沪宁电气化铁路对电网的影响研究.《江苏电机工程》.2008,第27卷(第05期), *

Also Published As

Publication number Publication date
CN101699693A (zh) 2010-04-28

Similar Documents

Publication Publication Date Title
CN101699693B (zh) 铁路牵引冲击负荷对电网电压波动计算方法
CN103124072B (zh) 考虑负荷特性的电网动态无功优化系统及方法
CN106240405B (zh) 单相供电装置、电气化铁路牵引供电系统及其控制方法
CN102998980B (zh) 一种电气化铁路牵引供电系统实时仿真建模方法
CN107732930A (zh) 一种适用于地铁牵引供电系统的多功能变流器系统
CN104617576B (zh) 考虑直流控制特性的多直流馈出交流电网故障计算方法
CN110504698A (zh) 一种电气化铁路同相供电综合补偿装置及其综合补偿方法
CN109217330B (zh) 一种电气化铁路同相供变电系统
CN102694386A (zh) 一种电气化铁路负序不平衡补偿方法
CN106911135A (zh) 一种集中补偿、混连式电动汽车充电系统
CN210224973U (zh) 一种电气化铁路同相供电综合补偿装置
CN203774792U (zh) 一种高速铁路电能质量及牵引网电压综合补偿系统
CN206149148U (zh) 一种基于光伏发电的铁路牵引供电系统
CN106849060A (zh) 一种含光伏电源的高速铁路牵引供电系统及控制方法
CN106356876A (zh) 高铁at供电方式下注入电力系统的负序电流计算模型
CN103475028B (zh) 制动回馈装置对供电系统的防逆流控制系统和控制方法
CN201681389U (zh) 电气化铁路供电系统损耗模型
CN206049416U (zh) 单相供电装置、电气化铁路牵引供电系统
CN105490274A (zh) 地铁供电系统电能质量控制系统及方法
CN103036235B (zh) 基于etap的电气化铁路电能质量综合治理方法
CN103419677A (zh) 一种含城市轨道交通的配电网等值电路
CN201685719U (zh) 一种电气化铁路轨道养护车辆的电源车
CN209159471U (zh) 用于电气化铁路牵引供电系统at供电方式的负序治理系统
CN203942277U (zh) 基于spwm技术的电气化铁路牵引变电所用电源
CN203941238U (zh) 一种35kV高压并联电力电容器组介损在线监测装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP02 Change in the address of a patent holder

Address after: 210036 No. 10 Dujiang Road, Gulou District, Nanjing City, Jiangsu Province

Patentee after: Jiangsu Keneng Electric Engineering Consultation Co., Ltd.

Address before: 210036 No. 96 Caochangmen Street, Gulou District, Nanjing City, Jiangsu Province

Patentee before: Jiangsu Keneng Electric Engineering Consultation Co., Ltd.

CP02 Change in the address of a patent holder