CN103675510A - 评测at方式下单线高架段电气化铁路雷击跳闸率的方法 - Google Patents

Abstract

本发明公开了一种测量AT供电方式下电气化铁路单线高架区段接触网雷击跳闸率的方法，该方法包括：第一步，获取电气化铁路线路参数；第二步，计算雷击类型分界点坐标，确定接触网不同雷击类型的影响区域；第三步，计算感应雷击跳闸率；第四步，计算绕击跳闸率；第五步，确定总跳闸率。该方法计算方便，且解决了AT供电方式下电气化铁路雷击跳闸率计算困难的问题。

Description

评测AT方式下单线高架段电气化铁路雷击跳闸率的方法

技术领域

本发明涉及一种测量AT供电方式下电气化铁路单线高架区段接触网雷击跳闸率的方法，特别是涉及一种基于电气几何模型计算AT供电方式下电气化铁路单线高架区段接触网感应雷击、绕击跳闸率的方法，适用于AT供电方式下电气化铁路防雷设计和防雷改造，属于铁路系统过电压领域。 

背景技术

牵引供电系统雷击跳闸严重影响了我国电气化铁路的安全稳定运行。为保证列车运行的可靠性，准确计算接触网供电线路的雷击跳闸率，评估其对牵引供电可靠性的影响，常需要确定接触网的雷击跳闸率。目前，AT供电方式下的电气化铁路接触网雷击类型分为感应雷击(雷击大地)、绕击(雷击正馈线、承力索、接触线)两种，而以往资料缺乏对AT供电方式下电气化铁路这两种雷击跳闸率进行准确计算的方法，为我国AT供电方式下电气化铁路进行有针对性地防雷设计与防雷改造带来了非常大的困难。 

发明内容

本发明的目的在于提供一种测量AT供电方式下电气化铁路单线高架区段接触网雷击跳闸率的方法，使用该方法可以计算AT供电方式下电气化铁路单线高架区段接触网每百公里的年感应雷击与绕击跳闸率。 

本发明实现上述目的的技术原理是利用传统电气几何模型分析AT供电方式下电气化铁路单线高架区段接触网的感应雷击、绕击情况，其原理如附图1所示。O为坐标原点，分别在桥面上以承力索、正馈线位置为圆心，以雷电先导对承力索击距r_c、雷电先导对正馈线击距r_g为半径作弧线，再以雷电先导对大地击距r_e作平行于大地的直线，分别相交于A、B、C点，其中H为高架桥对地高度，h_g为正馈线对高架桥面高度，h_c为承力索对高架桥面高度，a为承力索到支柱内侧的距离，b为正馈线到支柱内侧的距离。雷电落在A点左侧时击中地面，此时接触网高压导线上产生感应过电压，即发生感应雷击；雷电落在A、B两点以及B、 C两点中间时，分别击中正馈线、承力索，将在接触网上产生绕击过电压，即发生绕击；雷电落在C点右侧时击中地面，在接触网高压导线上产生感应过电压，即同样发生感应雷击。 

本发明解决其技术问题所采用的技术方案主要包括以下步骤： 

第一步，获取电气化铁路线路参数，包括高架桥对地高度，承力索、正馈线对高架桥面高度，承力索、正馈线到支柱内侧的距离，绝缘子串的50％冲击放电电压，正馈线半径，雷电日，落雷密度，绝缘子串平均运行电压梯度，支柱接地电阻，支柱等效电感，雷电流波头时间，电晕校正系数，支柱两侧相邻正馈线的电感并联值等。 

第二步，计算分界点A、C的坐标，确定接触网不同雷击类型的影响区域。建立如附图1所示的坐标系，A点对应的坐标为(x_a，y_a)，C点对应的坐标为(x_c，y_c)。此时感应雷击对应的区间为(-∞，x_a)和(x_c，+∞)，绕击对应的区间为(x_a，x_c)。根据各点的几何关系，各点坐标按以下算式确定： 

$\left\{\begin{array}{c}{x}_a=-\sqrt{{r_g}^2-{[r_e-(h_g+H)]}^2}-(a+b)\\ y_a=r_e\end{array}\right.$

$\left\{\begin{array}{c}{x}_c=\sqrt{{r_c}^2-{[r_e-(h_c+H)]}^2}\\ y_c=r_e\end{array}\right.$

式中：H为高架桥对地高度(单位：m)，h_g为正馈线对高架桥面高度(单位：m)，h_c为承力索对高架桥面高度(单位：m)，a为承力索到支柱内侧的距离(单位：m)，b为正馈线到支柱内侧的距离(单位：m)，r_c为雷电先导对承力索击距(单位：m)，r_g为雷电先导对正馈线击距(单位：m)，r_e为雷电先导对大地击距(单位：m)。 

r_c、r_g可采用以下经验公式进行计算： 

$r_c=r_g=a_0{I}^{b_0}$

或 $r_c=a_0{I}^{b_0}{(h_c+H)}^{c_0},$ $r_g=a_0{I}^{b_0}{(h_g+H)}^{c_0}$

式中：I为雷电流幅值(单位：kA)，H为高架桥对地高度，h_c为承力索对高架桥面高度，h_g为正馈线对高架桥面高度。 

a₀、b₀、c₀可根据现场实验或模拟线路实验结果进行取值，也可以参照电力系统的经验取以下数值： 

a₀＝10，b₀＝0.65；或a₀＝0.67，b₀＝0.74，c₀＝0.6；或a₀＝1.57，b₀＝0.69，c₀＝0.45。 

r_e可按下式计算： 

r_e＝k₂r_c

其中k₂为击距系数，计算公式如下： 

k₂＝1.066+(h_c+H)/216.45 

式中：H为高架桥对地高度，h_c为承力索对高架桥面高度。 

或k₂＝22/(h+H)，或k₂＝1.94-(h+H)/26， 

或k₂＝1.08-(h+H)/59，或k₂＝1.05-(h+H)/87。 

式中：H为高架桥对地高度，h为支柱对高架桥面高度(单位：m)。 

为简化计算，也可令r_c＝r_g＝r_e。 

第三步，计算感应雷击跳闸率。 

首先，根据下式计算感应雷击区间(-∞，x_a)和(x_c，+∞)的有效投影长度： 

$\mathrm{\&Delta;}{L}_a=\left\{\begin{array}{c}{x}_a-x_{\mathrm{ea}},x_{\mathrm{ea}}<x_a\\ 0,x_{\mathrm{ea}}\&GreaterEqual;x_a\end{array}\right.,$ $\mathrm{\&Delta;}{L}_c=\left\{\begin{array}{c}{x}_{\mathrm{ec}}-x_c,x_{\mathrm{ec}}<x_c\\ 0,x_{\mathrm{ec}}\&GreaterEqual;x_c\end{array}\right.$

其中： 

$x_{\mathrm{ea}}=\frac{I\&CenterDot;25(h_c+H)(1-k_0\frac{h_g+H}{h_c+H})}{U_{50\%}},$ $x_{\mathrm{ec}}=\frac{I\&CenterDot;25(h_c+H)(1-k_0\frac{h_g+H}{h_c+H})}{U_{50\%}}$

式中：H为高架桥对地高度，I为雷电流幅值，h_g为正馈线对高架桥面高度，h_c为承力索对高架桥面高度，k₀为正馈线与承力索之间的几何耦合系数，U_50％为绝缘子串的50％冲击放电电压(单位：kV)。 

k₀可按下式计算： 

$k_0=\frac{\ln \frac{d^{\&prime;}}{d}}{\ln \frac{2(h_g+H)}{r}}$

式中：d′为承力索与正馈线镜像间的距离(单位：m)，d为承力索与正馈线间的距离(单位：m)，r为正馈线半径(单位：m)。 

然后，按照下式计算感应雷击跳闸率： 

其中：

    为落雷密度(单位：次/km²·天)，T_d为雷电日(单位：天/年)，f(I)为雷电流概率密度，η为建孤率。 

f(I)可用以下经验公式进行计算： 

$f\left(I\right)=0.026\&times;{10}^{-\frac{I}{88}}$

或 $f\left(I\right)=0.052\&times;{10}^{-\frac{I}{44}}$

或 $f\left(I\right)=\frac{2.0}{{25}^{2.0}}\&times;\frac{I^{(2.0-1)}}{{[1+{(I/25)}^{2.0}]}^2}$

或 $f\left(I\right)=\frac{2.6}{{31}^{2.6}}\&times;\frac{I^{(2.6-1)}}{{[1+{(I/31)}^{2.6}]}^2}$

或 $f\left(I\right)=\frac{2.7}{{12}^{2.7}}\&times;\frac{I^{(2.7-1)}}{{[1+{(I/12)}^{2.7}]}^2}$

或 $f\left(I\right)=\frac{2.8}{{15.9}^{2.8}}\&times;\frac{I^{(2.8-1)}}{{[1+{(I/15.9)}^{2.8}]}^2}$

η的计算可按下式进行： 

η＝(4.5E^0.75-14)×10^-2

式中：E为绝缘子串平均运行电压梯度(单位：kV/m)。 

计算感应雷击跳闸率公式中的积分下限I_ea、I_ec按下式确定： 

$I_{\mathrm{ea}}=\frac{U_{50\%}\left|x_{\mathrm{cg}}\right|}{25(h_c+H)(1-k_0\frac{h_g+H}{h_c+H})},$ $I_{\mathrm{ec}}=\frac{U_{50\%}\left|x_{\mathrm{cg}}\right|}{25(h_c+H)(1-k_0\frac{h_g+H}{h_c+H})}$

其中： 

$x_{\mathrm{ag}}=\sqrt{{\left({{10I}_g}^{0.65}\right)}^2-{[{{10I}_g}^{0.65}-(h_g+H)]}^2}-(a+b),$ $x_{\mathrm{cg}}=\sqrt{{\left({{10I}_g}^{0.65}\right)}^2-{[{{10I}_g}^{0.65}-(h_c+H)]}^2}$

$I_g=\frac{U_{50\%}}{(1-k){\mathrm{\&beta;R}}_i+(\frac{h_g+H}{h_c+H}-k)\mathrm{\&beta;}\frac{L_t}{{\mathrm{\&tau;}}_f}+(1-\frac{h_g+H}{h_c+H}{k}_0)\frac{h_c+H}{{\mathrm{\&tau;}}_f}}$

式中：β为支柱分流系数，R_i为支柱接地电阻(单位：Ω)，k为正馈线与承力索之间的耦合系数，L_t支柱等效电感(单位：μH)，τ_f为雷电流波头时间(单位：μs)。 

其中k＝k₁k₀，k₁为电晕校正系数，接触网可取1.15。 

$\mathrm{\&beta;}=1/(1+\frac{L_t}{L_g}+\frac{R_i{\mathrm{\&tau;}}_f}{{2L}_g})$

式中：L_g为支柱两侧相邻正馈线的电感并联值(单位：μH)。 

积分上限I_max根据线路的重要程度或行业的具体要求取值，也可以按分布概 率为90％或99％时的雷电流幅值进行估算。 

第四步，按照下式计算绕击跳闸率： 

式中积分下限

积分上限I_max取值同上。 

第五步，接触网的总雷击跳闸率为感应雷击跳闸率与绕击雷击跳闸率之和，即按照下式计算总跳闸率： 

n＝n_gy+n_c

本发明的技术效果是采用电气几何模型，提出了一种AT供电方式下电气化铁路单线高架区段接触网的感应雷击、绕击跳闸率计算方法，解决了AT供电方式下电气化铁路雷击跳闸率计算困难的问题。 

附图说明

下面结合附图和实施例进一步说明本发明。 

图1为单线高架区段接触网电气几何模型示意图 

具体实施方式

下面通过实例，结合附图1，对本发明的技术方案进行进一步说明。 

第一步，获取线路参数。某铁路单线高架区段线路，高架桥高度10m，正馈线对高架桥面高度8m，承力索对高架桥面高度7.8m，承力索距支柱内侧3m，正馈线距支柱内侧0.8m，正馈线半径6.25mm，绝缘子U50％放电电压270kV，雷电流波头时间2.6μs，雷电日40天，绝缘子串平均运行电压梯度20.36kV，支柱接地电阻10Ω，支柱等效电感6.72μH，支柱两侧相邻正馈线的电感并联值36.85μH，电晕校正系数1.15。 

计算时令

且a₀＝10，b₀＝0.65；选取雷电流概率密度为  $f\left(I\right)=0.026\&times;{10}^{-\frac{I}{88}}.$ 第二步，计算感应雷击、绕击跳闸率积分上限与下限。 

分别利用上述公式计算感应雷击积分下限I_ea、I_ec为： 

I_ea＝39kA，I_ec＝42kA。 

积分上限I_max按分布概率为99％时的雷电流幅值进行估算： 

I_max＝176kA 

计算绕击积分上下限I_c、I_max为： 

I_c＝4kA，I_max＝176kA。 

第三步，计算感应雷击跳闸率。 

利用A、C坐标公式在感应雷击雷电流上下限区间范围I_ea～I_max、I_ec～I_max内计算A、C坐标分布，确定感应雷击有效区间，再利用下式计算感应雷击跳闸率： 

计算结果为n_gy＝0.9685次/100km·年。 

第四步，计算绕击跳闸率。 

利用A、C坐标公式在绕击雷电流上下限区间范围I_c～I_max内计算A、C坐标分布，确定绕击影响区间，再利用下式计算绕击跳闸率： 

计算结果为n_g＝8.1300次/100km·年。 

第五步，利用下式计算总跳闸率： 

nn＝n_gy+n_c

计算结果为n＝9.0985次/100km·年。 

Claims (4)

1.一种划分AT供电方式下电气化铁路单线高架区段接触网不同雷击类型的影响区域的方法，其特征在于它包括以下步骤：

第一步，获取电气化铁路线路参数，包括高架桥对地高度，承力索、正馈线对高架桥面高度，承力索、正馈线到支柱内侧的距离，绝缘子串的50％冲击放电电压，正馈线半径，雷电日，落雷密度，绝缘子串平均运行电压梯度，支柱接地电阻，支柱等效电感，雷电流波头时间，电晕校正系数，支柱两侧相邻正馈线的电感并联值等；

第二步，计算分界点A、C的坐标，计算公式如下：

$\left\{\begin{array}{c}{x}_a=-\sqrt{{r_g}^2-{[r_e-(h_g+H)]}^2}-(a+b)\\ y_a=r_e\end{array}\right.$

$\left\{\begin{array}{c}{x}_c=-\sqrt{{r_c}^2-{[r_e-(h_c+H)]}^2}\\ y_c=r_e\end{array}\right.$

式中：H为高架桥对地高度(单位：m)，h_g为正馈线对高架桥面高度(单位：m)，h_c为承力索对高架桥面高度(单位：m)，a为承力索到支柱内侧的距离(单位：m)，b为正馈线到支柱内侧的距离(单位：m)，r_c为雷电先导对承力索击距(单位：m)，r_g为雷电先导对正馈线击距(单位：m)，r_e为雷电先导对大地击距(单位：m)；

r_c、r_g、r_e计算公式如下：

r_c＝r_g＝r_e＝10I^0.65

式中：I为雷电流幅值(单位：kA)；

第三步，划分不同雷击类型的区域，其中感应雷的区域为(-∞，x_a)和(x_c，+∞)，绕击对应的区域为(x_a，x_c)。

2.一种测量AT供电方式下电气化铁路单线高架区段接触网雷击跳闸率的方法，其特征在于它包括以下步骤：

第一步，基于权利要求1，获取分界点A、C点的坐标；

第二步，按照下式计算感应雷击跳闸率：

其中：其中：

为落雷密度(单位：次/km²·天)，T_d为雷电日(单位：天/年)，f(I)为雷电流概率密度，η为建孤率，ΔL_a、ΔL_c为感应雷击区间的有效投影长度；

η计算公式如下：

η＝(4.5E^0.75-14)×10^-2

式中：E为绝缘子串平均运行电压梯度(单位：kV/m)；

雷电流概率密度f(I)计算公式如下：

$f\left(I\right)=0.026\&times;{10}^{-\frac{I}{88}}$

ΔLa、ΔLc计算公式如下：

$\mathrm{\&Delta;}{L}_a=\left\{\begin{array}{c}{x}_a-x_{\mathrm{ea}},x_{\mathrm{ea}}<x_a\\ 0,x_{\mathrm{ea}}\&GreaterEqual;x_a\end{array}\right.,$ $\mathrm{\&Delta;}{L}_c=\left\{\begin{array}{c}{x}_{\mathrm{ec}}-x_c,x_{\mathrm{ec}}<x_c\\ 0,x_{\mathrm{ec}}\&GreaterEqual;x_c\end{array}\right.$

其中：

$x_{\mathrm{ea}}=\frac{I\&CenterDot;25(h_c+H)(1-k_0\frac{h_g+H}{h_c+H})}{U_{50\%}},$ $x_{\mathrm{ec}}=\frac{I\&CenterDot;25(h_c+H)(1-k_0\frac{h_g+H}{h_c+H})}{U_{50\%}}$

式中：H为高架桥对地高度，I为雷电流幅值，h_g为正馈线对高架桥面高度，h_c为承力索对高架桥面高度，k₀为正馈线与承力索之间的几何耦合系数，U_50％为绝缘子串的50％冲击放电电压(单位：kV)；

k₀可按下式计算：

$k_0=\frac{\ln \frac{d^{\&prime;}}{d}}{\ln \frac{2(h_g+H)}{r}}$

式中：d′为承力索与正馈线镜像间的距离(单位：m)，d为承力索与正馈线间的距离(单位：m)，r为正馈线半径(单位：m)；

计算感应雷击跳闸率公式中的积分下限I_ea、I_ec按下式确定：

$I_{\mathrm{ea}}=\frac{U_{50\%}\left|x_{\mathrm{ag}}\right|}{25(h_c+H)(1-k_0\frac{h_g+H}{h_c+H})},$ $I_{\mathrm{ec}}=\frac{U_{50\%}\left|x_{\mathrm{cg}}\right|}{25(h_c+H)(1-k_0\frac{h_g+H}{h_c+H})}$

其中：

$x_{\mathrm{cg}}=\sqrt{{\left({{10I}_g}^{0.65}\right)}^2-{[{{10I}_g}^{0.65}-(h_c+H)]}^2}-(a+b),$ $x_{\mathrm{cg}}=\sqrt{{\left({{10I}_g}^{0.65}\right)}^2-{[10{I_g}^{0.65}-(h_c+H)]}^2}$

$I_g=\frac{U_{50\%}}{(1-k){\mathrm{\&beta;R}}_i+(\frac{h_g+H}{h_c+H}-k)\mathrm{\&beta;}\frac{L_t}{{\mathrm{\&tau;}}_f}+(1-\frac{h_g+H}{h_c+H}{k}_0)\frac{h_c+H}{{\mathrm{\&tau;}}_f}}$

式中：β为支柱分流系数，R_i为支柱接地电阻(单位：Ω)，k为正馈线与承力索之间的耦合系数，L_t支柱等效电感(单位：μH)，τ_f为雷电流波头时间(单位：μs)；

其中k＝k₁k₀，k₁为电晕校正系数，接触网可取1.15；

$\mathrm{\&beta;}=1/(1+\frac{L_t}{L_g}+\frac{R_i{\mathrm{\&tau;}}_f}{{2L}_g})$

式中：L_g为支柱两侧相邻正馈线的电感并联值(单位：μH)；

积分上限I_max根据线路的重要程度或行业的具体要求取值，也可以按分布概率为90％或99％时的雷电流幅值进行估算；

第三步，按照下式计算绕击跳闸率：

式中积分下限积分上限I_max取值同上；

第四步，按照下式计算接触网总雷击跳闸率：

n＝n_gy+n_c

式中：n_gy为感应雷击跳闸率，n_c为绕击跳闸率。

3.权利要求1中临界点A、C坐标计算公式中的击距r_c、r_g、r_e可用以下经验公式进行计算：

$r_c=r_g=a_0{I}^{b_0}$

或 $r_c=a_0{I}^{b_0}{(h_c+H)}^{c_0},$ $r_g=a_0{I}^{b_0}{(h_g+H)}^{c_0}$

式中：I为雷电流幅值(单位：kA)，H为高架桥对地高度，h_c为承力索对高架桥面高度，h_g为正馈线对高架桥面高度；

a₀、b₀、c₀可根据现场实验或模拟线路实验结果进行取值，也可以参照电力系统的经验取以下数值：

a₀＝10，b₀＝0.65；或a₀＝0.67，b₀＝0.74，c₀＝0.6；或a₀＝1.57，b₀＝0.69，c₀＝0.45；

r_e可按下式计算：

r_e＝k₂r_c

其中k₂为击距系数，计算公式如下：

k₂＝1.066+(h_c+H)/216.45

式中：H为高架桥对地高度，h_c为承力索对高架桥面高度；

或k₂＝22/(h+H)，或k₂＝1.94-(h+H)/26，

或k₂＝1.08-(h+H)/59，或k₂＝1.05-(h+H)/87；

式中：H为高架桥对地高度，h为支柱对高架桥面高度(单位：m)；

为简化计算，也可令r_c＝r_g＝r_e。

4.权利要求1中感应雷击、绕击跳闸率计算公式中的雷电流概率密度f(I)还可用以下经验公式进行计算：

$f\left(I\right)=0.052\&times;{10}^{-\frac{I}{44}}$

或 $f\left(I\right)=\frac{2.0}{{25}^{2.0}}\&times;\frac{I^{(2.0-1)}}{{[1+{(I/25)}^{2.0}]}^2}$

或 $f\left(I\right)=\frac{2.6}{{31}^{2.6}}\&times;\frac{I^{(2.6-1)}}{{[1+{(I/31)}^{2.6}]}^2}$

或 $f\left(I\right)=\frac{2.7}{{12}^{2.7}}\&times;\frac{I^{(2.7-1)}}{{[1+{(I/12)}^{2.7}]}^2}$

或 $f\left(I\right)=\frac{2.8}{{15.9}^{2.8}}\&times;\frac{I^{(2.8-1)}}{{[1+{(I/15.9)}^{2.8}]}^2}$

式中：I为雷电流幅值(单位：kA)。

Images