发明内容
本发明的目的在于,为了解决长期不受重视但十分重要且急需解决的配电线路直击雷跳闸率的计算问题,提出了全新的一种35kV及以下配电线路直击雷跳闸率计算方法。
本发明的技术方案是,一种35kV及以下配电线路直击雷跳闸率计算方法包括以下步骤,简明计算步骤示意图参见图1:
步骤一、获取配电线路所在地区雷害状况和线路、杆塔参数,包括落雷密度和回路数、电压等级、导线空间排列情况、塔形、杆塔波阻抗和接地电阻。
步骤二、根据导线空间排列情况,得到每次雷电流击中杆塔顶部与避雷线的概率和直接击中导线的概率。
步骤三、基于线路参数,根据全新的公式,计算得到雷电流直击导线导致绝缘子双相和三相闪络时的雷电流幅值。
步骤四、通过电磁暂态仿真软件建立特定回路和特定塔形情况下的雷电流反击模型,得到绝缘子双相和三相闪络时的雷电流幅值。
步骤五、根据直击、反击概率和对应的双相闪络、三相闪络的雷电流幅值,计算得出直击跳闸率和反击跳闸率,二者相加,即得总的直击雷跳闸率。
在本发明步骤一,收集雷击跳闸率计算所需的数据,为后续建模和计算做准备。
根据GIS系统确定配电线路的精确路径,利用雷电定位系统获得配电线路所在地区的落雷密度统计数据。确定线路电压水平、回路数、避雷线安装情况杆塔类型、导线空间排列情况、接地电阻等。
本发明步骤二,根据导线空间排列情况,得到每次雷击时,击中杆塔和线路的概率p,其中平原地区雷电流直接击中导线的概率p1和击中杆塔、避雷线的概率p2分别为0.5p。
本发明步骤三,计算雷电流直接击中导线时导致线路跳闸的雷电流幅值。雷电流直接击中导线的示意图如图2所示。由于配电线路中性点采取不接地或者经消弧线圈接地的非直接接地方式,因此,只有绝缘子双相闪络或三相闪络时才构成雷电流流通回路,形成稳定电弧,线路才会跳闸,这与主网不同。
绝缘子第i相和第j相先后闪络时雷电流幅值计算公式如下:
式中:Iij为雷电流幅值,kA;U50%为绝缘子50%击穿电压,kV;yj为第j相导线悬挂高度,m;kij为第i相导线对第j相导线的耦合系数,i,j=A,B,C,且i≠j;Lgt为杆塔电感,μH;Rch为杆塔冲击接地电阻,Ω。
分别求出任意两相先后闪络的雷电流幅值Iij,取其最小值作为双相闪络的雷电流I1,2。
绝缘子第i相、第j相和第k相先后闪络时雷电流幅值计算公式如下:
式中:kijk为第i、j相导线对第k相导线的耦合系数,i,j,k=A,B,C,且互不相同;yk第k相导线悬挂高度,m;
分别求出任意三相先后闪络的雷电流幅值Iij,取其最小值作为三相闪络的雷电流I1,2,3。
当应用于多回路配电线路的计算时,分别求出每一回路的双相闪络和三相闪络雷电流,分别取最小值最为整条线路的双相闪络和三相闪络雷电流。
绝缘子闪络后,形成稳定电弧的概率也与主网不同,可按下式计算:
式中:η1,2为绝缘子双相闪络后产生稳定电弧的概率;UN为系统标称(线)电压,kV;ldis为两相闪络时放电路径总长度,m。
本发明步骤四,通过电磁暂态仿真软件建立特定回路特定塔形情况下的雷电流反击模型,得到绝缘子双相闪络和三相闪络的雷电流幅值I'1,2和I'1,2,3。雷电流击中杆塔示意图如图3所示。
雷电流反击模型包括以下几个模块:线路模块、杆塔模块、雷电流模块和绝缘子闪络模块。
线路模块中可以基于现实情况自定义线路回路数和避雷线条数,设置导线电阻率、内外半径、横纵坐标、垂弧、线路长度、电压水平等,可以自定义是否加装避雷器、具体加装方式(全相加装或者单相加装等)以及避雷器的伏安特性等。
杆塔模块用单一波阻抗模拟,可以定义配电线路常用的各种塔形,包括十字型、上字型、门型等,并经接地电阻接地。
雷电流模块可以选择多种雷电流波形,并可自定义雷幅值、拨前波尾等雷电流参数。
本发明步骤五,根据直击、反击概率和对应的双相闪络、三相闪络的雷电流幅值,计算得出直击跳闸率和反击跳闸率。
我国一般地区雷电流幅值超过I的概率可按下式求得:
我国少雷地区雷电流幅值超过I的概率可按下式求得:
通过上式可以分别求得在雷电流直击导线时,导致绝缘子双相闪络和三相闪络的雷电流幅值概率P1,2和P1,2,3,以及雷电流击中杆塔和避雷线时,绝缘子双相闪络和三相闪络的雷电流幅值概率P′1,2和P′1,2,3。
雷电流直击导线导致的直击跳闸率可按下式计算:
n1=N·p1·[(P1,2-P1,2,3)+P1,2,3(2-η1,2)]·η1,2
雷击杆塔和避雷线导致的反击跳闸率可按下式计算:
n2=N·p2·[(P′1,2-P′1,2,3)+P′1,2,3(2-η1,2)]·η1,2
式中:n1为直击跳闸率,次/(100km·a);n2为反击跳闸率,次/(100km·a);N为落雷概率,次/(100km·a);p1为雷击导线的概率;p2为雷击杆塔和避雷线的概率;P1,2为雷击导线时导致绝缘子双相闪络概率;P1,2,3为雷击导线时导致绝缘子三相闪络概率;P`1,2为雷击杆塔、绝缘子时导致绝缘子双相闪络概率;P`1,2,3为雷击杆塔、绝缘子时导致绝缘子三相闪络概率;η1,2为绝缘子双相闪络后产生稳定工频续流的概率,即建弧率。
落雷概率N按下式计算:
N=0.1γT(b+10hgt)
式中:γ为地面落雷密度,次/(km2·a);T为雷电日数,计算时取40;
hgt为杆塔高度,m;b为避雷线间的距离,m,无避雷线时取0,单避雷线时取避雷线到杆塔中心线之间的垂直距离。
总的雷击跳闸率计算公式为:
n=n1+n2
直击雷跳闸率和感应雷跳闸率的计算步骤如图4所示。
本发明与原有的雷击跳闸率计算方法比较的有益效果如下:
(1)由于配电线路中性点采用非直接接地方式,这与主网不同,跳闸率计算方法也与主网不同。非直接接地方式下,只有绝缘子两相或两相以上闪络才可能导致线路跳闸,因此在计算直击导线跳闸率时,采用全新的公式计算出绝缘子两相闪络和三相闪络时的雷电流幅值,在计算反击跳闸率时,也通过仿真计算得到雷电流击中杆塔、避雷线时绝缘子两相闪络和三相闪络时的雷电流幅值。
(2)由于配电线路采用与主网不同的中性点非直接接地方式,绝缘子闪络后电弧轨迹也不同,因此采用新的两相闪络时的建弧率公式。
(3)得到直击和反击情况下的雷电流幅值和建弧率以后,继而采用的新的直击导线跳闸率和反击跳闸率的计算公式,计算更具有配电线路的针对性和特殊性,也为合理准确。
(4)本发明针对的是分布极为广泛,线路和地理环境复杂多变但却重视不够的35kV及以下的配电线路,提出了影响配电线路雷击跳闸率的所有情形的综合计算方法,具有全面和普遍适用性。包括雷击杆塔和避雷线造成的反击跳闸,雷击导线造成的直击跳闸。
(5)本发明对各种情形导致的雷击跳闸分别运用了不同的,但同时最合理准确的计算方法。其中直击跳闸率采用全新的理论公式法,反击跳闸率采用电磁暂态软件仿真建模与公式结合法。
(6)本发明的方法具有灵活性和针对性,可以自定义回路参数、避雷线参数、塔形、避雷器参数等,加之计算方法的差异性和合理性,保证了改计算方法的准确性和权威性。
本发明适用于配电线路不同塔形和回路下的跳闸率计算。
具体实施方式
本发明实施例一种35kV及以下配电线路直击雷跳闸率计算方法,包括以下步骤:
步骤一,获取配电线路所在地区雷害状况和线路、杆塔参数,包括落雷密度和回路数、电压等级、导线空间排列情况、塔形、杆塔波阻抗和接地电阻,为后续建模和计算做准备。
根据GIS系统确定配电线路的精确路径,利用雷电定位系统获得配电线路所在地区的落雷密度统计数据。确定线路电压水平、回路数、避雷线安装情况杆塔类型、导线空间排列情况、接地电阻等。
步骤二,根据导线空间排列情况,得到每次雷击时,击中杆塔和线路的概率p,其中雷电流击中杆塔顶部的概率p1和直接击中导线的概率p2分别为0.5p;
步骤三,基于线路参数,计算得到雷电流直击导线导致绝缘子双相和三相闪络时的雷电流幅值。雷电流直接击中导线示意图如图2所示。
由于配电线路中性点采取不接地或者经消弧线圈接地的非直接接地方式,因此,只有绝缘子双相闪络或三相闪络时才构成雷电流流通回路,形成稳定电弧,线路才会跳闸,这与主网不同。
绝缘子第i相和第j相先后闪络时雷电流幅值计算公式如下:
式中:Iij——雷电流幅值,kA;
U50%——绝缘子50%击穿电压,kV;
yj——第j相导线悬挂高度,m;
kij——第i相导线对第j相导线的耦合系数,i,j=A,B,C,且i≠j;
Lgt——杆塔电感,μH;
Rch——杆塔冲击接地电阻,Ω;
其中第i相导线对第j相导线的耦合系数kij按下式计算:
式中:Zij——第i相导线与第j相导线之间的互阻抗,Ω;
Zii——第i相导线的自阻抗,Ω;
dij'——第i相导线到第j相导线镜像之间的距离,m;
dij——第i相导线到第j相导线之间的距离,m;
hi——第i相导线的高度,m;
ri——第i相导线的半径,m;
分别求出任意两相先后闪络的雷电流幅值Iij,取其最小值作为双相闪络的雷电流I1,2。
绝缘子第i相、第j相和第k相先后闪络时雷电流幅值计算公式如下:
式中:kijk——第i、j相导线对第k相导线的耦合系数,i,j,k=A,B,C,且互不相同;
yk——第k相导线悬挂高度,m;
其中第i相导线和第j相导线对第k相导线的耦合系数kijk按下式计算:
分别求出任意三相先后闪络的雷电流幅值Iij,取其最小值作为三相闪络的雷电流I1,2,3。
当应用于多回路配电线路的计算时,分别求出每一回路的双相闪络和三相闪络雷电流,分别取最小值最为整条线路的双相闪络和三相闪络雷电流。
绝缘子闪络后,形成稳定电弧的概率也与主网不同,可按下式计算:
式中:η1,2——绝缘子双相闪络后产生稳定电弧的概率。
UN——系统标称(线)电压,kV;
ldis——两相闪络时放电路径总长度,m。
步骤四,通过电磁暂态仿真软件建立特定回路特定塔形情况下的雷电流反击模型,得到反击耐雷水平。该模型包括:线路模块、杆塔模块、雷电流模块、绝缘子闪络模块。雷电流击中杆塔的示意图如图3所示。
当雷电流击中杆塔顶端时,通过观察绝缘子状态,来判断是否闪络。当两相绝缘子闪络时的雷电流幅值即为配电线路耐雷水平。在计算耐雷水平时,可初步设定一个最大耐受电流IN和一个最小闪络电流IS,通过二分法,逐步缩小耐受电流和闪络电流的范围,接近真实耐雷水平值,当耐受电流和闪络电流的范围小于某个很小的确定值时,可停止二分,取耐受电流和闪络电流的平均数作为近似耐雷水平。上述过程可通过编程调用电磁暂态仿真模型,并自动实现循环过程,最终输出绝缘子双相闪络和三相闪络的雷电流幅值I'1,2和I'1,2,3。
我国一般地区雷电流幅值超过I的概率可按下式求得:
少雷地区雷电流幅值超过I的概率可按下式求得:
通过上式可以分别求得在雷电流直击导线时,导致绝缘子双相闪络和三相闪络的雷电流幅值概率P1,2和P1,2,3,以及雷电流击中杆塔和避雷线时,绝缘子双相闪络和三相闪络的雷电流幅值概率P′1,2和P′1,2,3。
雷电流直击导线导致的直击跳闸率可按下式计算:
n1=N·p1·[(P1,2-P1,2,3)+P1,2,3(2-η1,2)]·η1,2
雷击杆塔和避雷线导致的反击跳闸率可按下式计算:
n2=N·p2·[(P′1,2-P′1,2,3)+P′1,2,3(2-η1,2)]·η1,2
式中:n1,n2——直击跳闸率和反击跳闸率,次/(100km·a);
N——落雷概率,次/(100km·a);
p1,p2——雷击导线和雷击杆塔、避雷线概率;
P1,2,P1,2,3——雷击导线时导致绝缘子双相闪络和三相闪络概率;
P′1,2,P′1,2,3——雷击杆塔、绝缘子时导致绝缘子双相闪络和三相闪络概率;
η1,2——绝缘子双相闪络后产生稳定电弧的概率;
直击雷跳闸率和感应雷跳闸率的简明计算步骤如图1和图4所示。
落雷概率N按下式计算:
N=0.1γT(b+10hgt)
式中:γ——地面落雷密度,次/(km2·a);
T——雷电日数,计算时取40;
hgt——杆塔高度,m;
b——避雷线间的距离,m,无避雷线时取0,单避雷线时取避雷线到杆塔中心线之间的垂直距离。
总的雷击跳闸率计算公式为:
n=n1+n2。