评测复线高架区段电气化铁路雷击跳闸率的方法

1.一种测量电气化铁路复线高架区段接触网雷击跳闸率的方法，其特征在于它包括以下步骤：
第一步，获取电气化铁路线路参数，包括高架桥对地高度，承力索、回流线对高架桥面高度，承力索、回流线到支柱内侧的距离，复线上下行线路中心距离，绝缘子串的50％冲击放电电压，回流线半径，雷电日，落雷密度，绝缘子串平均运行电压梯度，支柱接地电阻，支柱等效电感，雷电流波头时间，电晕校正系数，支柱两侧相邻回流线的电感并联值；
第二步，计算单侧分界点A、B、C的坐标，计算公式如下：



式中：


2
λ3＝λ2+1
λ4＝2λ1λ2-2(hc+H)
λ5＝λ12+(hc+H)2-rc2
式中：H为高架桥对地高度，单位：m；hg为回流线对高架桥面高度，单位：m；hc为承力索对高架桥面高度，单位：m；a为承力索到支柱内侧的距离，单位：m；b为回流线到支柱内侧的距离，单位：m；c为复线上下行线路中心距离，单位：m；rc为雷电先导对承力索击距，单位：m；rg为雷电先导对回流线击距，单位：m；re为雷电先导对大地击距，单位：m；
第三步，划分不同雷击类型的影响区间，感应雷击对应的区间长度为(-∞，xa)区间的2倍，反击对应的区间长度为(xa，xb)的2倍，绕击对应的区间长度为(xb，xc)的2倍；
第四步，按照下式计算复线感应雷击跳闸率：

其中：Υ＝0.023Td0.3为落雷密度，单位：次/km2·天；Td为雷电日，单位：天/年；f(I)为雷电流概率密度，η为建孤率，ΔL为感应雷击区间的有效投影长度；
η计算公式如下：
η＝(4.5E0.75-14)×10-2
式中：E为绝缘子串平均运行电压梯度，单位：kV/m；
雷电流概率密度f(I)计算公式如下：

ΔL计算公式如下：

其中：

式中：H为高架桥对地高度，I为雷电流幅值，单位：kA，hg为回流线对高架桥面高度，hc为承力索对高架桥面高度，k0为回流线与承力索之间的几何耦合系数，U50％为绝缘子串的50％冲击放电电压，单位：kV；
k0按下式计算：

式中：d′为承力索与回流线镜像间的距离，单位：m；d为承力索与回流线间的距离，单位：m；r为回流线半径，单位：m；
计算感应雷击跳闸率公式中的积分下限Ie按下式确定：

其中：


式中：β为支柱分流系数，Ri为支柱接地电阻，单位：Ω；k为回流线与承力索之间的耦合系数，Lt为支柱等效电感，单位：μH；τf为雷电流波头时间，单位：μs；
其中k＝k1k0，k1为电晕校正系数，接触网取1.15；

式中：Lg为支柱两侧相邻回流线的电感并联值，单位：μH；
积分上限Imax根据线路的重要程度或行业的具体要求取值，或按分布概率为90％或
99％时的雷电流幅值进行估算；
第五步，按照下式计算复线反击雷击跳闸率：

积分下限Ig、积分上限Imax取值同上；
第六步，按照下式计算复线绕击跳闸率：

式中积分下限
积分上限Imax按xb＝c/2＝xc时进行求解；
式中：c为复线上下行线路中心距离；
若令rc＝rg＝re，则Imax按下式计算：

式中：

第七步，按照下式计算接触网复线总雷击跳闸率：
n＝nay+ng+nc
式中：ngy为复线感应雷击跳闸率，ng为复线反击雷击跳闸率，nc为复线绕击跳闸率。
2.一种根据权利要求1所述的方法，其中单侧分界点A、B、C坐标计算公式中的击距rc、rg用以下经验公式进行计算：

或
式中：I为雷电流幅值，单位：kA，H为高架桥对地高度，hc为承力索对高架桥面高度，hg为回流线对高架桥面高度；
a0、b0、c0根据现场实验或模拟线路实验结果进行取值，或参照电力系统的经验取以下数值：
a0＝10，b0＝0.65；或a0＝0.67，b0＝0.74，c0＝0.6；或a0＝1.57，b0＝0.69，c0＝0.45；
re按下式计算：
re＝k2rc
其中k2为击距系数，计算公式如下：
k2＝1.066+(hc+H)/216.45
式中：H为高架桥对地高度，hc为承力索对高架桥面高度；
或k2＝22/(h+H)，或k2＝1.94-(h+H)/26，
或k2＝1.08-(h+H)/59，或k2＝1.05-(h+H)/87；
式中：H为高架桥对地高度，h为支柱对高架桥面高度，单位：m；
或为简化计算，令rc＝rg＝re。

评测复线高架区段电气化铁路雷击跳闸率的方法
技术领域
[0001] 本发明涉及一种测量电气化铁路复线高架区段接触网雷击跳闸率的方法，特别是涉及一种基于电气几何模型计算电气化铁路复线高架区段接触网感应雷击、反击、绕击跳闸率的方法，适用于电气化铁路防雷设计和防雷改造，属于铁路系统过电压领域。
背景技术
[0002] 牵引供电系统雷击跳闸严重影响了我国电气化铁路的安全稳定运行。为保证列车运行的可靠性，准确计算接触网供电线路的雷击跳闸率，评估其对牵引供电可靠性的影响，常需要确定接触网的雷击跳闸率。目前，电气化铁路接触网雷击类型分为感应雷击(雷击大地)、反击(雷击回流线、支柱顶部等接地部分)、绕击(雷击承力索、接触线)三种，以往资料缺乏对电气化铁路这三种雷击跳闸率进行准确计算的方法，为我国电气化铁路进行有针对性地防雷设计与防雷改造带来了非常大的困难。
发明内容
[0003] 本发明的目的在于提供一种测量电气化铁路复线高架区段接触网雷击跳闸率的方法，使用该方法可以计算电气化铁路复线高架区段接触网每百公里的年感应雷击跳闸率、反击跳闸率与绕击跳闸率。
[0004] 本发明实现上述目的的技术原理是利用传统电气几何模型分析电气化铁路复线高架区段接触网的感应雷击、反击、绕击情况，其原理如附图1所示。复线接触网区段呈对称结构，假设左侧为上行线，右侧为下行线，O为坐标原点，分别以上行、下行线的承力索、回流线位置为圆心，以雷电先导对承力索击距rc、雷电先导对回流线击距rg为半径作弧线，再以雷电先导对大地击距re作平行于大地的直线，分别相交于A、B、C、A′、B′点，其中H为高架桥对地高度，hg为回流线对高架桥面高度，hc为承力索对高架桥面高度，a为承力索到支柱内侧的距离，b为回流线到支柱内侧的距离，c为复线上下行线路中心距离。雷电落在A点左侧和A′右侧时击中地面，此时接触网高压导线上产生感应过电压，即发生感应雷击；雷电落在A、B两点中间和A′、B′两点中间时，击中接触网的回流线，在接触网高压导线上产生反击过电压，即发生反击；雷电落在B、B′两点中间时，击中接触网高压导线，接触网上产生绕击过电压，即发生绕击。
[0005] 通常情况下接触网上下行结构相同，为了简化计算，可以假设A点左侧与A′点右侧的感应雷击区间长度完全相等，AB与A′B′反击区间长度完全相等，BC与CB′绕击区间长度完全相等，此时感应雷击对应的区间长度为(-∞，xa)区间的2倍，反击对应的区间长度为(xa，xb)的2倍，绕击对应的区间长度为(xb，xc)的2倍。因此，上行线路的感应雷击、反击、绕击跳闸率分别等于下行线路的感应雷击、反击、绕击跳闸率，即复线高架的感应雷击、反击、绕击跳闸率分别是单侧线路(上行或下行)感应雷击、反击、绕击跳闸率的2倍。
[0006] 本发明解决其技术问题所采用的技术方案主要包括以下步骤：
[0007] 第一步，获取电气化铁路线路参数，包括高架桥对地高度，承力索、回流线对高架桥面高度，承力索、回流线到支柱内侧的距离，复线上下行线路中心距离，绝缘子串的50％冲击放电电压，回流线半径，雷电日，落雷密度，绝缘子串平均运行电压梯度，支柱接地电阻，支柱等效电感，雷电流波头时间，电晕校正系数，支柱两侧相邻回流线的电感并联值等。
[0008] 第二步，确定接触网不同雷击类型的影响区域，计算单侧分界点A、B、C的坐标。
[0009] 建立如附图1所示的坐标系，A点对应的坐标为(xa，ya)，B点对应的坐标为(xb，yb)，C点对应的坐标为(xc，yc)，B′点对应的坐标为(xb′，yb′)，A′点对应的坐标为(xa′，ya′)。此时感应雷击对应的区间为(-∞，xa)和(xa′，+∞)，反击对应的区间为(xa，xb)和(xb′，xa′)，绕击对应的区间为(xb，xb′)。
[0010] 根据各点的几何关系，计算单侧分界点A、B、C的坐标，各点坐标按以下算式确定：
[0011]
[0012]
[0013]
[0014] 式中：
[0015]
[0016]
[0017] λ3＝λ22+1
[0018] λ4＝2λ1λ2-2(hc+H)
[0019] λ5＝λ12+(hc+H)2-rc2
[0020] 式中：H为高架桥对地高度(单位：m)，hg为回流线对高架桥面高度(单位：m)，hc为承力索对高架桥面高度(单位：m)，a为承力索到支柱内侧的距离(单位：m)，b为回流线到支柱内侧的距离(单位：m)，c为复线上下行线路中心距离(单位：m)，rc为雷电先导对承力索击距(单位：m)，rg为雷电先导对回流线击距(单位：m)，re为雷电先导对大地击距(单位：m)。
[0021] rc、rg可采用以下经验公式进行计算：
[0022]
[0023] 或
[0024] 式中：I为雷电流幅值(单位：kA)，H为高架桥对地高度，hc为承力索对高架桥面高度，hg为回流线对高架桥面高度。
[0025] a0、b0、c0可根据现场实验或模拟线路实验结果进行取值，也可以参照电力系统的经验取以下数值：
[0026] a0＝10，b0＝0.65；或a0＝0.67，b0＝0.74，c0＝0.6；或a0＝1.57，b0＝0.69，c0＝
0.45。
[0027] re可按下式计算：
[0028] re＝k2rc
[0029] 其中k2为击距系数，计算公式如下：
[0030] k2＝1.066+(hc+H)/216.45
[0031] 式中：H为高架桥对地高度，hc为承力索对高架桥面高度。
[0032] 或k2＝22/(h+H)，或k2＝1.94-(h+H)/26，
[0033] 或k2＝1.08-(h+H)/59，或k2＝1.05-(h+H)/87。
[0034] 式中：H为高架桥对地高度，h为支柱对高架桥面高度(单位：m)。
[0035] 为简化计算，也可令rc＝rg＝re。
[0036] 第三步，计算复线感应雷击跳闸率。
[0037] 首先，根据下式计算单侧感应雷击区间(-∞，xa)的有效投影长度：
[0038]
[0039] 其中：
[0040]
[0041] 式中：H为高架桥对地高度，I为雷电流幅值，hg为回流线对高架桥面高度，hc为承力索对高架桥面高度，k0为回流线与承力索之间的几何耦合系数，U50％为绝缘子串的50％冲击放电电压(单位：kV)。
[0042] k0可按下式计算：
[0043]
[0044] 式中：d′为承力索与回流线镜像间的距离(单位：m)，d为承力索与回流线间的距离(单位：m)，r为回流线半径(单位：m)。
[0045] 然后，按照下式计算单侧感应雷击跳闸率：
[0046]
[0047] 则复线感应雷击跳闸率为：
[0048]
[0049] 其中： 为落雷密度(单位：次/km2·天)，Td为雷电日(单位：天/年)，f
(I)为雷电流概率密度，η为建孤率。
[0050] f(I)可用以下经验公式进行计算：
[0051]
[0052] 或
[0053] 或
[0054] 或
[0055] 或
[0056] 或
[0057] η的计算可按下式进行：
[0058] η＝(4.5E0.75-14)×10-2
[0059] 式中：E为绝缘子串平均运行电压梯度(单位：kV/m)。
[0060] 计算复线感应雷击跳闸率公式中的积分下限Ie按下式确定：
[0061]
[0062] 其中：
[0063]
[0064]
[0065] 式中：β为支柱分流系数，Ri为支柱接地电阻(单位：Ω)，k为回流线与承力索之间的耦合系数，Lt支柱等效电感(单位：μH)，τf为雷电流波头时间(单位：μs)。
[0066] 其中k＝k1k0，k1为电晕校正系数，接触网可取1.15。
[0067]
[0068] 式中：Lg为支柱两侧相邻回流线的电感并联值(单位：μH)。
[0069] 积分上限Imax根据线路的重要程度或行业的具体要求取值，也可以按分布概率为
90％或99％时的雷电流幅值进行估算。
[0070] 第四步，按照下式计算复线反击跳闸率：
[0071]
[0072] 式中： 为单侧反击跳闸率。
[0073] 积分下限Ig、积分上限Imax取值同上。
[0074] 第五步，按照下式计算复线绕击跳闸率：
[0075]
[0076] 式中： 为单侧绕击跳闸率。
[0077] 积分下限
[0078] 积分上限Imax可按xb＝c/2＝xc时进行求解。
[0079] 式中：c为复线上下行线路中心距离。
[0080] 若令rc＝rg＝re，则Imax可按下式计算：
[0081]
[0082] 式中：
[0083]
[0084] 第六步，复线接触网的总雷击跳闸率为感应、反击与绕击跳闸率三者之和，即按照下式计算复线总跳闸率：
[0085] n＝ngy+ng+nc
[0086] 本发明的技术效果是采用电气几何模型，提出了一种电气化铁路复线高架区段接触网的感应雷击、反击、绕击跳闸率计算方法，解决了电气化铁路雷击跳闸率计算困难的问题。
附图说明
[0087] 下面结合附图和实施例进一步说明本发明。
[0088] 图1为复线高架区段接触网电气几何模型示意图
具体实施方式
[0089] 下面通过实例，结合附图1，对本发明的技术方案进行进一步说明。
[0090] 第一步，获取线路参数。某铁路复线高架区段线路，高架桥对地高度10m，回流线对高架桥面高度8m，承力索对高架桥面高度7.8m，承力索距支柱内侧3m，回流线距支柱内侧
0.8m，复线上下行线路的中心距离6m，回流线半径6.25mm，绝缘子U50％放电电压270kV，雷电流波头时间2.6μs，雷电日40天，绝缘子串平均运行电压梯度20.36kV，支柱接地电阻10Ω，支柱等效电感6.72μH，支柱两侧相邻回流线的电感并联值36.85μH，电晕校正系数1.15。
[0091] 计算时令 且a0＝10，b0＝0.65；选取雷电流概率密度为
[0092] 第二步，计算感应雷击、反击、绕击跳闸率积分上限与下限。
[0093] 利用上述公式计算感应雷击积分下限Ie为：
[0094] Ie＝39kA
[0095] Imax按分布概率为99％时的雷电流幅值进行估算：
[0096] Imax＝176kA
[0097] 计算反击积分上下限Ig、Imax为：
[0098] Ig＝25kA，Imax＝176kA
[0099] 计算绕击积分上下限Ic、Imax为：
[0100] Ic＝4kA，Imax＝31kA
[0101] 第三步，计算复线感应雷击跳闸率。
[0102] 利用A、B、C坐标公式在感应雷击雷电流上下限区间范围Ie～Imax内计算A、B、C坐标分布，确定单侧感应雷击有效区间，再利用下式计算复线感应雷击跳闸率：
[0103]
[0104] 计算结果为ngy＝0.8981次/100km·年。
[0105] 第四步，计算复线反击跳闸率。
[0106] 利用A、B、C坐标公式在反击雷电流上下限区间范围Ig～Imax内计算A、B、C坐标分布，确定单侧反击影响区间，再利用下式计算复线反击跳闸率：
[0107]
[0108] 计算结果为ng＝6.2447次/100km·年。
[0109] 第五步，计算复线绕击跳闸率。
[0110] 利用A、B、C坐标公式在绕击雷电流上下限区间范围Ic～Imax内计算A、B、C坐标分布，确定单侧绕击影响区间，再利用下式计算复线绕击跳闸率：
[0111]
[0112] 计算结果为ng＝0.1327次/100km·年。
[0113] 第六步，利用下式计算总跳闸率：
[0114] n＝ngy+ng+nc
[0115] 计算结果为n＝7.2755次/100km·年。