一种at单线供电系统牵引网谐振频率的确定方法
【技术领域】
[0001] 本发明属电气化铁道技术领域，涉及一种AT单线供电系统牵引网谐振频率的确 定方法。
【背景技术】
[0002] 随着交-直-交机车在电气化铁路中的大量应用，牵引供电系统的谐波特性发生 变化，原交-直型机车富含的3、5、7次等低次谐波含量大大降低，而高次谐波含量明显增 加，频谱变宽。当电力机车注入的谐波与牵引网的谐振频率匹配时，将造成谐波放大、激发 高次谐振，产生过电压、过电流，造成保护动作、设备损失。牵引网谐振已成为各级部门广泛 关注的问题，为了保证电气化铁路安全运营，通过简便、有效的方法直接计算牵引网谐振频 率，对工程设计和实际运营均有重要意义。
[0003]目前关于谐振的研究大多是采用仿真工具建立牵引网分布参数模型，通过仿真扫 描电流、电压或阻抗频谱，确定谐振是否发生及谐振频率，分析过程较为复杂。AT供电是我 国高速铁路的主流供电方式，对AT供电牵引网，即便经合并化简也存在接触网、正馈线、钢 轨至少三根导线，直接应用必然导致偏差。因此针对AT单线供电系统，急需发明简便有效 的牵引网谐振频率的确定方法。
【发明内容】
[0004]针对现有技术的状况，本发明的目的是提供一种AT单线供电系统牵引网谐振频 率的确定方法，使之克服现有技术的以上不足，为后续谐振抑制处理提供指导和参考。
[0005]本发明解决上述技术问题所采用的技术方案是：
[0006] -种AT单线供电系统牵引网谐振频率的确定方法，通过对单线AT供电方式下的 牵引网进行合并化简和去耦，得到接触网对地单位长度等效电容，进而获得牵引网谐振频 率，包含如下主要步骤
[0007]第一步：由牵引网横向截面图获得导线i对地高度h、半径ri，导线i与导线j之 间的距离、导线i与导线j的镜象之间的距离，并确定空气介电常数e。;可由下式计 算出导线i的自电位系数Pn、导线i与导线j之间的互电位系数P,，从而得到牵引网的电 位系数矩阵P:
[0008]
12 式中，下标1、]_可分别取1(：、1?1、1?2、？、？，其中1代表接触线、(：代表承力索、1?1 代表钢轨一、R2代表钢轨二、P代表保护线、F代表正馈线； 2 第二步：对6X6的电位系数矩阵P求逆，得到6X6的电容系数矩阵B，B的各行、 列的元素记为blS，下标i、j可分别取J、C、R1、R2、P、F;
[0011] 第三步：将牵引网进行等值合并，接触线（J)和承力索（C)合并成接触网等值导线 (T)，两条钢轨（Ri、R2)和保护线（P)合并成钢轨等值导线（R)，从而由接触网等值导线（T)、 钢轨等值导线（R)和正馈线（F)构成三导线的牵引网，其电容系数矩阵B'为3X3矩阵，各 行、列元素为：
[0012]
[0013] 第四步：得到单位长度的接触网对地等值电容C。：
[0014] C〇=bTT-bTF
[0015] 第五步：AT单线供电系统牵引网谐振频率f确定为：
[0016]
[0017] 式中，LT为牵引变压器折算到二次侧的等效电感，1为牵引网长度；
[0018] 第六步：输出牵引网谐振频率f至后续谐振抑制处理
[0019] 本发明涉及的原理如下：
[0020] 设牵引网多导体传输线数目为m，则对地电压Ui与其线电荷密度Qi之间存在下列 关系
[0021]
[0022] 忒由
[0023]
线i的自电位系数（km/F)
[0024]
导线i与导线j之间的互电位系数（km/F)
[0025] Qi-导线i的线电荷密度（C/km)
[0026] e。一空间介电系数，e。= 8. 8542X10 1F/km
[0027] r;一导线半径（m)
[0028] h-导线i的对地高度（m)
[0029] 屯一导线i与导线j之间的距离（m)
[0030] Dij-导线i与导线j之间的镜像距离（m)
[0031] 上述各导线及其镜像分布图见图1。
[0032] 将式（1)写成
[0033] U=PQ(2)
[0034] 所以
[0035] 〇 =P:11 =RIJ⑶
[0036]
(4.)
[0037] P矩阵为电位系数矩阵，B矩阵为电容系数矩阵，通过P矩阵求逆即可得到。
[0038] 由于牵引网的平行多导体传输线数目较多，单线AT供电方式牵引网有6根导线： 接触线（J)、承力索（C)、钢轨汛）、钢轨（R2)、保护线（P)、正馈线（F)，其横截面图见图2。 可根据纵向连接点对多导线进行等值合并：由于接触线和承力索每隔6_7m用一根吊弦连 接，并且承力索和接触线在短距离内可认为是平行传输的，所以可以将同一条线路的接触 线（J)和承力索（C)合并成接触网等值导线（T);同时可以将同一条线路的两条钢轨（&、 R2)和保护线（P)合并成钢轨等值导线（R)，导线合并图见图3。
[0039] 等值合并应该满足下列条件：
[0040] ①合并后的等值导线流过的电流等于需要合并的各根导线电流之和
[0041]
(5)
[0042] ②需要合并的各根导线单位长压降相等
[0043]
(6、
[0044] ③合并后的等值导线电荷等于需要合并的各根导线电荷之和
[0045]
C7)
[0046] ④需要合并的各根导线对地电压相等
[0047]
(8)
[0048] 导线等值合并后的电容系数矩阵
[0049] 单线AT供电方式牵引网6根导线电容系数矩阵方程为
[0050]
[0051] 将第一列、第二列按矩阵乘法展开并相加得：
[0052]Qc+Q】=[(bcc+bjc)Uc+(bcj+bjj)Uj] + (10)
[0053] [ (bCR1+bjR1)UR1+ (bCR2+bjR2)UR2+ (bcp+bjp)Up] + [bCF+bjF]UF
[0054] 将合并条件式（6) - (8)代入上式，得：
[0055]Qt=(bcc+2bcj+bjj)UT+ (bCR1+bCR2+bCP+bjR1+bjR2+bJP)UR+ (bCF+bJF)UF (11)
[0056] 同理，可将第三、四、五列按矩阵乘法展开并相加，结合合并条件，可得到：
[0057]
[0058]
[0059] 合并后的接触网等值导线T、钢轨等值导线R以及正馈线F三导体的部分电容分布 如图4所示，其中CTT、CFF、CRR为导体对地的部分电容，CTF、CFR、CTR为导体之间的部分电容。 T、R、F三导体系统的各部分电容取值为：
[0060]
[0061] 由图4可知：
[0062]QT=CTTUT+CTR (UT_UR)+CTF (UT_UF) = (CTT+CTR+CTF)UT_CTRUR _CTFUF (14)
[0063] 对应式（12)可知
[0064]
[0065] 且在AT供电系统中，UT=-UF，U产0,故式（14)可改写为
[0066]Qt=(Ctt+Ctr+2Ctf)Ut (16)
[0067] 所以得到单位长度的接触网对地等值电容C。：
[0068] C〇=Ctt+Ctr+2Ctf=bTT-bTF (17)
[0069] 理论上可将牵引网线路进行简化，可采用如图5所示的分布参数电路来等效。图 中，
[0070]
〖网对地等值线路特征阻抗
[0071]
接触网对地等值线路传播常数
[0072]Z。一一单位长度接触网对地等值阻抗
[0073]C。--单位长度接触网对地等值电容
[0074]Ls一一牵引变压器折算到二次侧的等值电感之和
[0075] 1--牵引网长度
[0076] x一一机车到牵引变压器的距离
[0077] 从负荷处向变电所看的等值阻抗ZiS:
[0078]
(18):
[0079] 从负荷处向分区所看的等值阻抗厶为：
[0080]
'^9)
[0081] 负荷处看到的阻抗Z为：
[0082]
(2Q)
[0083] 当Z=m时，牵引网将发生谐振，于是谐振条件是为
[0084]j〇Lsshy1+Zcchy1 = 0 (21)
[0085]
(22)
[0086]由干Y1 <<〗，thY1~y1，代入上式得
[0087]
(乃）
[0088] 于是，牵引网谐振频率
[0089]
_
[0090] 本发明方法计算方便、精度较高，计算所得的牵引网谐振频率可作为牵引供电系 统设计阶段谐振抑制的依据。
【附图说明】
[0091] 图1为导线及其镜像图
[0092] 图2为单线AT供电方式牵引网悬挂导线图
[0093] 图3为单线AT供电方式牵引网导线合并图
[0094] 图4为合并后牵引网三导体电容分布图
[0095] 图5为分布参数电路图
[0096] 图6为本发明方法的步骤框图。
【具体实施方式】
[0097] 下面结合附图和实施例，对本发明作进一步的详细说明。
[0098] 如图2所示的牵引网配置如下：
[0099] 表1导线型号及位置
[0102] 表2各导线等效半径与电阻参数
[0103]
[0104] 电容系数矩及电容系数矩阵的求解：
[0105] 由上文知可利用公5
ft求解多导线系统的电 位系数P，对P求逆矩阵即可得到电容系数矩阵B，其具体计算方法上文已做详细说明，这里 直接给出电位系数矩阵及电容系数矩阵的计算结果：
[0106] 表3电位系数矩阵P (F/km) 1
[0107]
[0108](注：以上各参数均X10°)
[0109] 表4电容系数矩阵B (nF/km)
[0110]
[0111]
[0112] 将接触线、承力索等值成一根导线（T线），两条钢轨与保护线合成一根（R线），由 式（12)得到合并后的电容系数矩阵为
[0113] 表5合并后的电容系数矩阵B'（nF/km)
[0114]
[0115] 由式（13)，将表5再次变换为分布电容矩阵C，如表6。
[0116] 表6合并后的分布电容矩阵（nF/km)
[0117]
[0118] 因此由式（17)得单位长度的接触网对地等值电容为12. 59InF。
[0119] 系统及牵引变压器折算到二次侧等值电感为
[0120]
[0121] 因此由式（24)得牵引网谐振频率为
[0122]
[0123]对应为28. 21次谐波，由于牵引网富含谐波为奇次谐波，实际取29次谐波。
【主权项】
1. 一种AT单线供电系统牵引网谐振频率的确定方法，通过对单线AT供电方式下的 牵引网进行合并化简和去耦，得到接触网对地单位长度等效电容，进而获得牵引网谐振频 率，包含如下主要步骤 第一步：由牵引网横向截面图获得导线i对地高度Ii1、半径Γι，导线i与导线j之间的 距离Cl1,、导线i与导线j的镜象之间的距离Du，并确定空气介电常数ε。;可由下式计算出 导线i的自电位系数Pn、导线i与导线j之间的互电位系数Ρ,，从而得到牵引网的电位系 数矩阵P :式中，下标i、j可分别取J、C、札、R2、P、F，其中J代表接触线、C代表承力索、Rl代表 钢轨一、R2代表钢轨二、P代表保护线、F代表正馈线； 第二步：对6X6的电位系数矩阵P求逆，得到6X6的电容系数矩阵B，B的各行、列的 元素记为blS，下标i、j可分别取J、C、R1、R2、P、F ; 第三步：将牵引网进行等值合并，接触线（J)和承力索（C)合并成接触网等值导线 (T)，两条钢轨（R^R2)和保护线（P)合并成钢轨等值导线（R)，从而由接触网等值导线（T)、 钢轨等值导线（R)和正馈线（F)构成三导线的牵引网，其电容系数矩阵B'为3X3矩阵，各 行、列元素为：第四步：得到单位长度的接触网对地等值电容C。： Cq t) TT bff 第五步：AT单线供电系统牵引网谐振频率f确定为：式中，Lt为牵引变压器折算到二次侧的等效电感，1为牵引网长度； 第六步：输出牵引网谐振频率f至后续谐振抑制处理。
【专利摘要】本发明公开了一种AT单线供电系统的牵引网谐振频率的确定方法。该方法通过对单线AT供电方式下的牵引网进行合并化简和去耦，得到接触网对地单位长度等效电容，获得牵引网谐振频率。本发明方法计算方便、精度较高，计算所得的牵引网谐振频率可作为牵引供电系统设计阶段谐振抑制的依据。
【IPC分类】G06F19/00, H02J3/00
【公开号】CN105186505
【申请号】CN201510624842
【发明人】郭蕾, 周利军, 杨佳
【申请人】西南交通大学
【公开日】2015年12月23日
【申请日】2015年9月28日