CN103761414A - 输电线路的雷击停运概率分析方法 - Google Patents
输电线路的雷击停运概率分析方法 Download PDFInfo
- Publication number
- CN103761414A CN103761414A CN201310700689.4A CN201310700689A CN103761414A CN 103761414 A CN103761414 A CN 103761414A CN 201310700689 A CN201310700689 A CN 201310700689A CN 103761414 A CN103761414 A CN 103761414A
- Authority
- CN
- China
- Prior art keywords
- lightning
- thunder
- thunderbolt
- time
- probability
- 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.)
- Granted
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及输电线路的雷击停运概率分析方法,本发明方法基于三时次雷区信息对电网雷击跳闸率进行预测,采用三时次雷区信息对雷电进行短时预测,识别出预测雷电发生范围内的杆塔,计算杆塔在预测雷电强度下的雷击停运概率,从而获得线路的预测雷击停运概率。该方法以实时雷电监测信息为基础,考虑了不同雷电天气间的差异性,使线路雷击停运概率的预测更具有实际意义。
Description
技术领域
本发明涉及电网气象风险分析方法。
背景技术
近年来,随着电网的快速发展和强对流天气的增多,雷害故障频繁发生,雷击成为我国输电线路安全可靠运行的主要危害。如2009年7月的某大型水电站近区电网500kV双回线路遭受雷击相继跳闸,造成水电厂切机甩负荷;2009年8月浙江某变电站220kV线路遭遇多重雷击,开关断口被击穿,造成设备损坏、220kV母线失电。根据电网故障分类统计,在高压线路的总跳闸次数中,由雷击所引起的跳闸次数占40%~70%[1],尤其是在多雷、土壤电阻率高、地形复杂的地区,雷击输电线路所引起的故障率更高。因而若能实现电网因雷击所引起的故障概率预测、分析,并进而形成相应的预警系统,对增强电网抵抗因雷击所引起的故障具有重要意义。若要进行电网的雷击故障概率预测和分析,首先需要实现对雷电发生范围的预报,进而挑选出可能受雷电影响的输电线路,并计算线路雷击跳闸概率。对于雷电发生范围的预测,涉及到许多其它气象及自然环境因素,且其发生及消失过程非常快,如何进行准确其准确预报一直是气象预报领域的热点课题之一。文献[2]提出了基于雷达回波信息的TITAN系统,对雷暴进行识别、追踪、分析和临近预报。文献[3]在此基础上,进一步提出了雷暴单体识别追踪算法(SCIT),为下一代气象雷达(WSR-88D)的雷达回波信息识别设置了多个阈值,增强了对雷暴单体的识别能力。文献[4]考虑到雷暴气象生消时间短、变化快的特点,将临近预报点的时空分辨率设定在0.5km-1min,以提高预报的准确度。文献[5]则提出利用连续三时间段的雷电统计数据代替连续两时间段的统计数据进行分析,获得雷电发展的变化量的同时,也能获得雷电发展的变化率,以得到更准确的雷区预测。雷击对于电网运行的影响,以及雷击引起电网故障率预测和分析,不少专家和学者对此进行了有意义的研究和探索。文献[6]基于有限的数个严重雷暴事件及相关的其它 气象数据和电网故障数据,采用数理统计分析,而建立了雷暴引起的电网故障率回归模型;但由于雷暴事件比较少,且事件之间间隔比较长,故基于此所建立的模型的推广性有很大的局限性。文献[7]对比了实测落雷密度下的雷击跳闸率与规程落雷密度下的理论雷击跳闸率,发现线路的雷击故障情况时空差异性大,认为根据实际雷电活动分布进行线路雷击跳闸率的预测更为准确。文献[8]利用时变矩阵分析两相邻时段雷电区域的关联,从而预测未来雷电发生区域,通过以预测雷电区域的预测落雷密度来动态评估线路雷击跳闸率,但关联分析法并不十分适用于临近预报点的时空分辨率较小的情况,另外,线路跳闸率动态评估方法中只考虑了线路的实时落雷情况,而忽略了雷电强度的影响。可以看出,目前线路雷电停运率的计算方法主要借鉴了线路防雷设计的近似公式,但防雷公式建立在模拟实验结果统计的基础之上,反映的是线路年平均雷击风险,而选择采用统一的雷电流幅值分布也无法体现出地域间雷电天气的差异性。各种雷电停运率计算方法各有优点,均能判别出输电线路间受雷电影响的差异,但不适用于描述线路的实时雷电停运故障情况。
参考文献:
[1]张田乾,张希峰.架空电力线路的防雷保护[J].电源技术应用,2005,8(7):41-42.
[2]Dixon M,Wiener G.TITAN:Thunderstorm identification,tracking,analysis,and nowcasting-A radar-based methodology[J].Journal of Atmospheric and Oceanic Technology,1993,10(6):785-797.
[3]Johnson J T,MacKeen P L,Witt A,et al.The storm cell identification and tracking algorithm:An enhanced WSR-88D algorithm[J].Weather and Forecasting,1998,13(2):263-276.
[4]Ruzanski E,Chandrasekar V,Wang Y.The CASA nowcasting system[J].Journal of Atmospheric and Oceanic Technology,2011,28(5):640-655.
[5]Peters S,Meng L.Visual Analysis for Nowcasting of Multidimensional Lightning Data[J].ISPRS International Journal of Geo-Information,2013,2(3):817-836.
[6]Bonelli P,Marcacci P,Bertolotti E,et al.Nowcasting and assessing thunderstorm risk on the Lombardy region(Italy)[J].Atmospheric Research,2011,100(4):503-510.
[7]黄伟超,何俊佳,陆佳政,等.实际雷电活动分布下的线路雷击跳闸率计算[J].高电压技术,2008,34(7):1368-1373.
[8]谢云云,薛禹胜,王昊昊,等.电网雷击故障概率的时空在线预警[J].电力系统自动化,2013,37(17):44-51.
发明内容
本发明要解决的技术问题是提供一种电网雷击停运概率的预测方法,该方法该方法以实时雷电监测信息为基础,考虑了不同雷电天气间的差异性,使线路雷击停运概率的预测更具有实际意义。
本发明的具体实现方法包括以下步骤:
1.设定统计时间步长,称为一个时次,统计连续三时次的雷电监测实时信息,识别每个时次的雷电分区。
2.对相邻两时次的雷电分区进行关联,得到最优轨迹,连续三时次雷电信息可获得两段相连的雷电分区最优轨迹。
3.根据连续三时次的雷电分区的最优轨迹,预测下一时次的雷电发生情况,包括雷电分区的位置、面积、雷电流强度和落雷密度。
4.根据杆塔坐标信息筛选出预测雷电分区内的杆塔,根据杆塔在雷电分区中的位置获得杆塔处落雷概率。
5.在预测雷电流强度下,根据受雷电影响杆塔的杆塔参数,采用改进电气几何模型计算杆塔绕击跳闸概率,利用蒙特卡洛模拟方法计算杆塔反击跳闸概率,得到杆塔雷电跳闸概率。
6.综合杆塔处落雷概率和杆塔雷电跳闸概率,获得杆塔预测雷击停运概率,从而获得线路预测雷击停运概率。
以上各步骤具体为以下方法实现:
输电线路的雷击停运概率分析方法,设定统计时间段步长,称为一个时次,统计连续三时次t-2,t-1,t的雷电监测实时信息,对每一个时次进行雷电分区的识别,每一时次的雷电分区识别具体步骤如下:
步骤1:对每个时次的雷电统计信息进行预处理,首先将监控地域按经纬度方向 划分成面积相等的小方格,将其等价于一个2D网格图像,以经纬度方向为横纵坐标建立新坐标系,小方格的边长为单位长度,以小方格右下角点的坐标作为小方格的标记,取方格大小为0.05°×0.05°,然后将统计的落雷点根据落雷经纬度信息累计到对应的小方格;
步骤2:根据小方格的累计落雷个数n,按公式(1)进行全局阈值分割,对表示监控区域的网格图像进行二值化处理,以二值化落雷网格表示监测地域的落雷情况,
步骤3:在落雷网格中,利用八邻域边界跟踪算法对落雷区域进行轮廓提取,并统计区域内落雷信息,落雷网络轮廓提取采用3×3网格的八邻域分析;
步骤4:以圆形代替各个雷区的不规则轮廓,作为雷电发生区域,并求取各雷电发生区域的特性参数,根据公式—计算各落雷区域中心点经纬度L(x,y)、各雷区所有落雷点距离相应中心点的最大距离r及落雷密度p;以L(x,y)为中心r为半径画圆,即为实际雷电发生区域,各个圆的中心坐标L(x,y)、半径r和雷电分区落雷密度p为相应雷电分区特性参数,也为当前时次雷电分区的识别结果;
r=Max{[(xi-x)2+(yi-y)2]1/2} (3)
式中,n为落雷区域内的落雷总数;(xi,yi)为第i个落雷的经纬度,i=1,…,n;N为当前时次的落雷总数。
输电线路的雷击停运概率分析方法,连续三时次t-2,t-1,t中,分别对连续两时次的雷电分区进行追踪,即寻求t-2时次雷电分区到t-1时次雷电分区的最优轨迹和t-1时次雷电分区到t时次雷电分区的最优轨迹,对于连续的 两个时次t1和t2,t1时次有n1个雷电分区,t2时次有n2个雷电分区,(x1i,y1i)表示t1时次第i个雷电分区的位置,(x2j,y2j)表示t2时次第j个雷电分区的位置,为了求取两个时次所有雷区之间的最有可能轨迹,建立如下优化模型:
约束条件:
Cij=ω1dp+ω2dA (6)
dp=[(x1i-x2j)2+(y1i-y2j)2]1/2 (7)
若n1<n2,则 否则
其中:Xij为相邻时次雷电分区关联矩阵的第i行第j列元素,表示t1时次的i雷电分区到t2时次的j雷电分区的关联情况;Cij为价值矩阵第i行第j列元素,表示t1时次的i雷电分区到t2时次的j雷电分区的变化价值,价值函数值越小表示两雷电分区关联性越强;dp表示雷电分区位移的变化;dA表示雷电分区面积的变化;ω1和ω2为权重;基于统计分析,雷区的移动速度可设定阈值vmax=60km/h;两时次的时间间隔为Δt,若dL/Δt>vmax,设定Cij为一个较大的值,使得该轨迹不会成为最佳匹配结果。
输电线路的雷击停运概率分析方法,基于连续三时次t-2,t-1,t的雷电分区信息及相邻两时次的雷电分区最优轨迹进行下一时次的雷电分区预测,具体预测步骤如下:
其中,(xt-2,yt-2),(xt-1,yt-1),(xt,yt)分别表示t-2,t-1,t三个时次的雷电分区位置,rt-2,rt-1,rt分别表示t-2,t-1,t三个时次的雷电分区的半径;
步骤3:按式(16)-(19)预测t+1时次的每一雷电分区位置(xt+1,yt+1)和雷电分区的特性(包括雷电流It+1、半径rt+1和落雷密度pt+1),由于雷电流I具有很强的随机性,本文在预测中使t+1时次各雷电分区的雷电流预测值等于与之相关的t时次雷电分区雷电流,以保证:
rt+1=βtrt (17)
It+1=max{It,1,It,2,…It,m} (18)
其中,对于t+1时次的每一雷电分区,m为t时次与之相关联的雷电分区数目,It,1,It,2,…It,m分别为t时次与之相关联的m个雷电分区的雷电流,pt,k为t时次与之相关联的m个雷电分区中第k个雷电分区的落雷密度。
输电线路的雷击停运概率分析方法,根据杆塔的经纬度信息识别出位于预测雷电分区内的杆塔,并根据杆塔在雷电分区内的位置,计算杆塔的落雷概率,对于每个杆塔,以杆塔为中心取垂直线路方向±350m、线路方向1/2档距的区域为杆塔落雷概率计算的有效区域,杆塔的落雷概率等于该区域的落雷概率,利用公式(20)计算杆塔的落雷概率,
式中,A为杆塔所在雷电分区的面积;A′为杆塔落雷概率计算的有效区域与杆塔所在雷电分区的重合面积;p为杆塔所在雷电分区的单落雷密度。
输电线路的雷击停运概率分析方法,计算杆塔的雷击跳闸概率,雷击跳闸概率为绕击跳闸概率和反击跳闸概率的总和,如式(21)所示,
P=Ps+Pc (21)
步骤1:采用电气几何模型计算绕击跳闸概率;通过几何分析,考虑避雷线路、三相线路间的相互屏蔽作用,以屏蔽弧的比率计算得到线路的绕击概率和绕击跳闸概率,电气几何模型分析过程中,引入地面击距rg和导线(避雷线)击距rc(rs),考虑雷电先导对大地和线路击距的差异性;引入地面倾斜角θg和雷电入射角来考虑地形与雷电入射方向对输电线路的绕击率的影响;通过旋转坐标系的方法按雷电垂直入射平原地区的情形进行分析,以暴露弧段的比值计算绕击概率Pα;最后,如式(22)所示,得到绕击跳闸概率Ps,
Ps=ηPα (22)
式中,η为建弧率,
步骤2:应用蒙特卡洛模拟方法计算反击跳闸概率,设定模拟次数为N,以yk表 示第k次模拟的结果,若反击引起闪络,则yk=1,否则yk=0,每一次雷击模拟过程中,首先计算出杆塔的绕击跳闸概率,产生一个[0,1]均匀分布的随机数r1,判断是否发生绕击,若r1>Pα,则认为没有发生绕击;若没有发生绕击,则再产生一个[0,1]均匀分布的随机数r2,判断是否发生反击,若r2<g,则发生反击;若发生反击,则继续产生一个[0,1]均匀分布的随机数r3,得到工频相角φ及三相工频电压瞬时值,如式(23)-(26)所示;最后根据闪络判据判断是否发生闪络,N次模拟结束后,如式(27)、(27)所示,得反击跳闸率的渐进统计估计值P,并计算反击跳闸概率Pc,
φ=2π·r3 (23)
UA=Umsinφ (24)
UB=Umsin(φ-2π/3) (25)
UC=Umsin(φ+2π/3) (26)
Pc=ηP (28)
输电线路的雷击停运概率分析方法,基于杆塔的预测落雷概率及预测雷击跳闸概率,计算线路的预测雷击停运概率,同一线路上的杆塔之间为可靠性逻辑串联关系,对于有n个杆塔的线路,如式(29)所示,计算其雷击停运概率:
式中,Pline为线路的雷击停运概率;n为该线路在雷区内的杆塔数;PiT为第i基杆塔的落雷概率;Pi为第i基杆塔的雷击跳闸概率,为第i基杆塔的反击跳闸率Pic和绕击跳闸概率Pis的和;PiTPi称为第i基杆塔的雷击停运概率。
本发明基于三时次雷区信息对电网雷击跳闸率进行预测,采用三时次雷区信息对雷电进行短时预测,识别出预测雷电发生范围内的杆塔,计算杆塔在预测雷 电强度下的雷击停运概率,从而获得线路的预测雷击停运概率。该方法以实时雷电监测信息为基础,考虑了不同雷电天气间的差异性,使线路雷击停运概率的预测更具有实际意义。
附图说明
图1输电线路预测雷击停运概率计算流程
图22D形式的雷电监测数据识别示例
图3区域识别流程
图4落雷区域识别结果
图5相邻时次雷电分区移动的所有可能轨迹
图6推广匈牙利算法求解相邻雷电分区间最优轨迹的流程图
图7雷电区域的位置预测
图8杆塔落雷概率计算的有效区域
图9电气几何模型图
具体实施方式:
下面结合附图对本发明方法进行详细描述,附图1为本发明方法的流程图。1、设定统计时间段步长,称为一个时次,统计连续三时次t-2,t-1,t的雷电监测实时信息,对每一个时次进行雷电分区的识别。每一时次的雷电分区识别具体步骤如下:
步骤1:对每个时次的雷电统计信息进行预处理。首先将监控地域按经纬度方向划分成面积相等的小方格,将其等价于一个2D网格图像。以经纬度方向为横纵坐标建立新坐标系,小方格的边长为单位长度,以小方格右下角点的坐标作为小方格的标记。取方格大小为0.05°×0.05°。然后将统计的落雷点根据落雷经纬度信息累计到对应的小方格,如附图2所示。
步骤2:根据小方格的累计落雷个数n,按公式(1)进行全局阈值分割,对表示监控区域的网格图像进行二值化处理。以二值化落雷网格表示监测地域的落雷情况。
步骤3:在落雷网格中,利用八邻域边界跟踪算法对落雷区域进行轮廓提取,并统计区域内落雷信息。落雷网络轮廓提取采用3×3网格的八邻域分析,流程见附图3。
步骤4:以圆形代替各个雷区的不规则轮廓,作为雷电发生区域,并求取各雷电发生区域的特性参数。根据公式-计算各落雷区域中心点经纬度L(x,y)、各雷区所有落雷点距离相应中心点的最大距离r及落雷密度p;以L(x,y)为中心r为半径画圆,即为实际雷电发生区域。各个圆的中心坐标L(x,y)、半径r和雷电分区落雷密度p为相应雷电分区特性参数,也为当前时次雷电分区的识别结果。附图4中的虚圆型即为辨识之后雷区边界。
r=Max{[(xi-x)2+(yi-y)2]1/2} (3)
式中,n为落雷区域内的落雷总数;(xi,yi)为第i个落雷的经纬度,i=1,…,n;N为当前时次的落雷总数。
2、连续三时次t-2,t-1,t中,分别对连续两时次的雷电分区进行追踪,即寻求t-2时次雷电分区到t-1时次雷电分区的最优轨迹和t-1时次雷电分区到t时次雷电分区的最优轨迹。
对于连续的两个时次t1和t2,如附图5所示,t1时次有n1个雷电分区,t2时次有n2个雷电分区。(x1i,y1i)表示t1时次第i个雷电分区的位置,(x2j,y2j)表示t2时次第j个雷电分区的位置。为了求取两个时次所有雷区之间的最有可能轨迹,建立如下优化模型:
约束条件:
Cij=ω1dp+ω2dA (6)
dp=[(x1i-x2j)2+(y1i-y2j)2]1/2 (7)
若n1<n2,则 否则
其中:Xij为相邻时次雷电分区关联矩阵的第i行第j列元素,表示t1时次的i雷电分区到t2时次的j雷电分区的关联情况;Cij为价值矩阵第i行第j列元素,表示t1时次的i雷电分区到t2时次的j雷电分区的变化价值,价值函数值越小表示两雷电分区关联性越强;dp表示雷电分区位移的变化;dA表示雷电分区面积的变化;ω1和ω2为权重;基于统计分析,雷区的移动速度可设定阈值vmax=60km/h;两时次的时间间隔为Δt,若dL/Δt>vmax,设定Cij为一个较大的值,使得该轨迹不会成为最佳匹配结果;约束条件保证了t2时次的每个雷电分区只由一个t1时次的雷电分区发展而来或t1时次的每个雷电分区只会发展到t2时次的一个雷电分区,体现了雷区的分裂和融合。
对于由式(9)所组成的优化模型,本文采用推广匈牙利算法求解最优轨迹,流程如附图6所示,具体求解思路如下:
2)设定m=Min{n1,n2};
3)利用标准匈牙利算法进行求解:
①变换价值方阵Cn×n,使方阵的各行各列均出现0元素,变换原则为:①从价值方阵的每行元素中减去该行的最小元素,②再从所得价值方阵的每列元素中减去该列的最小元素;
④反复进行步骤①②,直到所有0元素都被标记和划掉为止;
⑤判断标记出的数目是否等于m,若相等,转到步骤4),若不等,转到步骤⑥;
⑥调整价值方阵:取出未被直线穿过的数据区中最小元素,该数据区所有元素减去该最小元素,同时直线相交处每个元素加上该最小元素,去掉直线和圈,更新价值方阵,转到步骤①进行计算;
4)获取新价值方阵:当n1>n2时,取剩余列构成新价值方阵,当n1<n2时,取剩余行构成新价值方阵,并取m为新价值方阵的行(列)数,转到步骤3);
3、基于连续三时次t-2,t-1,t的雷电分区信息及相邻两时次的雷电分区最优轨迹进行下一时次的雷电分区预测。具体预测步骤如下:
其中,(xt-2,yt-2),(xt-1,yt-1),(xt,yt)分别表示t-2,t-1,t三个时次的雷电分区位置,rt-2,rt-1,rt分别表示t-2,t-1,t三个时次的雷电分区的半径;
步骤3:按式(17)-(20)预测t+1时次的每一雷电分区位置(xt+1,yt+1)和雷电分区的特性(包括雷电流It+1、半径rt+1和落雷密度pt+1),由于雷电流I具有很强的随机性,本文在预测中使t+1时次各雷电分区的雷电流预测值等于与之相关的t时次雷电分区雷电流,以保证:
rt+1=βtrt (18)
It+1=max{It,1,It,2,…It,m} (19)
其中,对于t+1时次的每一雷电分区,m为t时次与之相关联的雷电分区数目,It,1,It,2,…It,m分别为t时次与之相关联的m个雷电分区的雷电流,pt,k为t时次与之相关联的m个雷电分区中第k个雷电分区的落雷密度。
4、根据杆塔的经纬度信息识别出位于预测雷电分区内的杆塔,并根据杆塔在雷电分区内的位置,计算杆塔的落雷概率。如图8所示,对于每个杆塔,以杆塔为中心取垂直线路方向±350m、线路方向1/2档距的区域为杆塔落雷概率计算的有效区域,杆塔的落雷概率等于该区域的落雷概率。利用公式(21)计算杆塔的落雷概率。
式中,A为杆塔所在雷电分区的面积;A′为杆塔落雷概率计算的有效区域与杆塔所在雷电分区的重合面积;p为杆塔所在雷电分区的单落雷密度。
5、计算杆塔的雷击跳闸概率。雷击跳闸概率为绕击跳闸概率和反击跳闸概率的总和,如式(22)所示。
P=Ps+Pc (22)
步骤1:采用电气几何模型计算绕击跳闸概率。电气几何模型如附图9所示,通过几何分析,考虑避雷线路、三相线路间的相互屏蔽作用,以屏蔽弧的比率计算得到线路的绕击概率和绕击跳闸概率,该方法综合联系了雷电放电特性、线路结构尺寸、线路所处地形等实际因素,可较准确描述线路的绕击特性。电气几何模型分析过程中,通过引入地面击距rg和导线(避雷线)击距rc(rs),考虑雷电先导对大地和线路击距的差异性;引入地面倾斜角θg和雷电入射角来考虑地形与雷电入射方向对输电线路的绕击率的影响;通过旋转坐标系的方法按雷电垂直入射平原地区的情形进行分析,以暴露弧段的比值计算绕击概率Pα;最后,如式(23)所示,得到绕击跳闸概率Ps。
Ps=ηPα (23)
式中,η为建弧率。
步骤2:应用蒙特卡洛模拟方法计算反击跳闸概率。蒙特卡罗方法是一种统计计算方法,能够模拟雷击的随机过程,考虑雷电对三相输电线路的不同影响,更加符合实际。设定模拟次数为N,以yk表示第k次模拟的结果,若反 击引起闪络,则yk=1,否则yk=0。每一次雷击模拟过程中,首先计算出杆塔的绕击跳闸概率,产生一个[0,1]均匀分布的随机数r1,判断是否发生绕击,若r1>Pα,则认为没有发生绕击;若没有发生绕击,则再产生一个[0,1]均匀分布的随机数r2,判断是否发生反击,若r2<g,则发生反击;若发生反击,则继续产生一个[0,1]均匀分布的随机数r3,得到工频相角φ及三相工频电压瞬时值,如式(24)-(27)所示;最后根据闪络判据判断是否发生闪络。N次模拟结束后,如式(28)、(29)所示,得反击跳闸率的渐进统计估计值P,并计算反击跳闸概率Pc。
φ=2π·r3 (24)
UA=Umsinφ (25)
UB=Umsin(φ-2π/3) (26)
UC=Umsin(φ+2π/3) (27)
Pc=ηP (29)
6、基于杆塔的预测落雷概率及预测雷击跳闸概率,计算线路的预测雷击停运概率。同一线路上的杆塔之间为可靠性逻辑串联关系,对于有n个杆塔的线路,如式(30)所示,计算其雷击停运概率:
式中,Pline为线路的雷击停运概率;n为该线路在雷区内的杆塔数;PiT为第i基杆塔的落雷概率;Pi为第i基杆塔的雷击跳闸概率,为第i基杆塔的反击跳闸率Pic和绕击跳闸概率Pis的和;PiTPi称为第i基杆塔的雷击停运概率。
Claims (7)
1.输电线路的雷击停运概率分析方法,其特征在于包括以下步骤:
(1)设定统计时间步长,称为一个时次,统计连续三时次的雷电监测实时信息,识别每个时次的雷电分区;
(2)对相邻两时次的雷电分区进行关联,得到最优轨迹,连续三时次雷电信息可获得两段相连的雷电分区最优轨迹;
(3)根据连续三时次的雷电分区的最优轨迹,预测下一时次的雷电发生情况,包括雷电分区的位置、面积、雷电流强度和落雷密度;
(4)根据杆塔坐标信息筛选出预测雷电分区内的杆塔,根据杆塔在雷电分区中的位置获得杆塔处落雷概率;
(5)在预测雷电流强度下,根据受雷电影响杆塔的杆塔参数,采用改进电气几何模型计算杆塔绕击跳闸概率,利用蒙特卡洛模拟方法计算杆塔反击跳闸概率,得到杆塔雷电跳闸概率;
(6)综合杆塔处落雷概率和杆塔雷电跳闸概率,获得杆塔预测雷击停运概率,从而获得线路预测雷击停运概率。
2.根据权利要求1所述的输电线路的雷击停运概率分析方法,其特征在于设定统计时间段步长,称为一个时次,统计连续三时次t-2,t-1,t的雷电监测实时信息,对每一个时次进行雷电分区的识别,每一时次的雷电分区识别具体步骤如下:步骤1:对每个时次的雷电统计信息进行预处理,首先将监控地域按经纬度方向划分成面积相等的小方格,将其等价于一个2D网格图像,以经纬度方向为横纵坐标建立新坐标系,小方格的边长为单位长度,以小方格右下角点的坐标作为小方格的标记,取方格大小为0.05°×0.05°,然后将统计的落雷点根据落雷经纬度信息累计到对应的小方格;
步骤2:根据小方格的累计落雷个数n,按公式(1)进行全局阈值分割,对表示监控区域的网格图像进行二值化处理,以二值化落雷网格表示监测地域的落雷情况,
步骤3:在落雷网格中,利用八邻域边界跟踪算法对落雷区域进行轮廓提取,并 统计区域内落雷信息,落雷网络轮廓提取采用3×3网格的八邻域分析;
步骤4:以圆形代替各个雷区的不规则轮廓,作为雷电发生区域,并求取各雷电发生区域的特性参数,根据公式(2)—(4)计算各落雷区域中心点经纬度L(x,y)、各雷区所有落雷点距离相应中心点的最大距离r及落雷密度p;以L(x,y)为中心r为半径画圆,即为实际雷电发生区域,各个圆的中心坐标L(x,y)、半径r和雷电分区落雷密度p为相应雷电分区特性参数,也为当前时次雷电分区的识别结果;
r=Max{[(xi-x)2+(yi-y)2]1/2} (3)
式中,n为落雷区域内的落雷总数;(xi,yi)为第i个落雷的经纬度,i=1,…,n;N为当前时次的落雷总数。
3.根据权利要求1所述的输电线路的雷击停运概率分析方法,其特征在于连续三时次t-2,t-1,t中,分别对连续两时次的雷电分区进行追踪,即寻求t-2时次雷电分区到t-1时次雷电分区的最优轨迹和t-1时次雷电分区到t时次雷电分区的最优轨迹,对于连续的两个时次t1和t2,t1时次有n1个雷电分区,t2时次有n2个雷电分区,(x1i,y1i)表示t1时次第i个雷电分区的位置,(x2j,y2j)表示t2时次第j个雷电分区的位置,为了求取两个时次所有雷区之间的最有可能轨迹,建立如下优化模型:
约束条件:
Cij=ω1dp+ω2dA (6)
dp=[(x1i-x2j)2+(y1i-y2j)2]1/2 (7)
4.根据权利要求1所述的输电线路的雷击停运概率分析方法,其特征在于基于连续三时次t-2,t-1,t的雷电分区信息及相邻两时次的雷电分区最优轨迹进行下一时次的雷电分区预测,具体预测步骤如下:
其中,(xt-2,yt-2),(xt-1,yt-1),(xt,yt)分别表示t-2,t-1,t三个时次的雷电分区位置, rt-2,rt-1,rt分别表示t-2,t-1,t三个时次的雷电分区的半径;
步骤3:按式(16)-(19)预测t+1时次的每一雷电分区位置(xt+1,yt+1)和雷电分区的特性(包括雷电流It+1、半径rt+1和落雷密度pt+1),由于雷电流I具有很强的随机性,本文在预测中使t+1时次各雷电分区的雷电流预测值等于与之相关的t时次雷电分区雷电流,以保证:
rt+1=βtrt (17)
It+1=max{It,1,It,2,…It,m} (18)
其中,对于t+1时次的每一雷电分区,m为t时次与之相关联的雷电分区数目,It,1,It,2,…It,m分别为t时次与之相关联的m个雷电分区的雷电流,pt,k为t时次与之相关联的m个雷电分区中第k个雷电分区的落雷密度。
6.根据权利要求1所述的输电线路的雷击停运概率分析方法,其特征在于计算杆塔的雷击跳闸概率,雷击跳闸概率为绕击跳闸概率和反击跳闸概率的总和,如式(21)所示,
P=Ps+Pc (21)
步骤1:采用电气几何模型计算绕击跳闸概率;通过几何分析,考虑避雷线路、三相线路间的相互屏蔽作用,以屏蔽弧的比率计算得到线路的绕击概率和绕击跳闸概率,电气几何模型分析过程中,引入地面击距rg和导线(避雷线)击距rc(rs),考虑雷电先导对大地和线路击距的差异性;引入地面倾斜角θg和雷电入射角来考虑地形与雷电入射方向对输电线路的绕击率的影响;通过旋转坐标系的方法按雷电垂直入射平原地区的情形进行分析,以暴露弧段的比值计算绕击概率Pα;最后,如式(22)所示,得到绕击跳闸概率Ps,
Ps=ηPα (22)
式中,η为建弧率,
步骤2:应用蒙特卡洛模拟方法计算反击跳闸概率,设定模拟次数为N,以yk表示第k次模拟的结果,若反击引起闪络,则yk=1,否则yk=0,每一次雷击模拟过程中,首先计算出杆塔的绕击跳闸概率,产生一个[0,1]均匀分布的随机数r1,判断是否发生绕击,若r1>Pα,则认为没有发生绕击;若没有发生绕击,则再产生一个[0,1]均匀分布的随机数r2,判断是否发生反击,若r2<g(g为击杆率),则发生反击;若发生反击,则继续产生一个[0,1]均匀分布的随机数r3,得到工频相角φ及三相工频电压瞬时值,如式所示(23)-(26)所示;最后根据闪络判据判断是否发生闪络,N次模拟结束后,如式(27)、(28)所示,得反击跳闸率的渐进统计估计值P,并计算反击跳闸概率Pc,
φ=2π·r3 (23)
UA=Umsinφ (24)
UB=Umsin(φ-2π/3) (25)
UC=Umsin(φ+2π/3) (26)
Pc=ηP (28) 。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310700689.4A CN103761414B (zh) | 2013-12-18 | 2013-12-18 | 输电线路的雷击停运概率分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310700689.4A CN103761414B (zh) | 2013-12-18 | 2013-12-18 | 输电线路的雷击停运概率分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103761414A true CN103761414A (zh) | 2014-04-30 |
CN103761414B CN103761414B (zh) | 2017-05-24 |
Family
ID=50528650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310700689.4A Expired - Fee Related CN103761414B (zh) | 2013-12-18 | 2013-12-18 | 输电线路的雷击停运概率分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103761414B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105067904A (zh) * | 2015-07-17 | 2015-11-18 | 云南电力试验研究院(集团)有限公司 | 一种基于历史雷电定位数据和电网跳闸数据的雷电预警数据评估方法 |
CN105868872A (zh) * | 2016-05-30 | 2016-08-17 | 东北大学 | 一种配电网雷电灾害故障预测方法 |
CN106443342A (zh) * | 2016-09-29 | 2017-02-22 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 高压直流输电线路雷击故障点的识别方法 |
CN107633320A (zh) * | 2017-08-17 | 2018-01-26 | 广东电网有限责任公司惠州供电局 | 一种基于气象预测和风险评估的电网线路重要度评估方法 |
CN108693450A (zh) * | 2018-04-09 | 2018-10-23 | 国网电力科学研究院武汉南瑞有限责任公司 | 用于大跨越高杆塔的输电线路绕击闪络耐雷性能分析方法 |
CN108898258A (zh) * | 2018-07-06 | 2018-11-27 | 广州供电局有限公司 | 雷电灾害天气下电力系统连锁故障风险的分析方法与系统 |
CN109541409A (zh) * | 2018-11-30 | 2019-03-29 | 国家电网有限公司 | 基于电气几何模型的配电线路闪络率改进算法 |
CN110851791A (zh) * | 2019-11-05 | 2020-02-28 | 云南电网有限责任公司 | 一种雷电灾害下线路故障概率计算方法及系统 |
CN112001070A (zh) * | 2020-08-10 | 2020-11-27 | 国网山东省电力公司滨州供电公司 | 一种输电线受外部环境影响停运概率的建模方法 |
CN114154766A (zh) * | 2022-02-07 | 2022-03-08 | 国网江苏省电力有限公司苏州供电分公司 | 雷电动态预测下电网的易损性预警方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009264943A (ja) * | 2008-04-25 | 2009-11-12 | Chubu Electric Power Co Inc | 落雷予報装置および落雷予報方法 |
CN102590651A (zh) * | 2011-12-31 | 2012-07-18 | 国网电力科学研究院 | 基于雷电实测数据的输电线路故障概率评估方法 |
-
2013
- 2013-12-18 CN CN201310700689.4A patent/CN103761414B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009264943A (ja) * | 2008-04-25 | 2009-11-12 | Chubu Electric Power Co Inc | 落雷予報装置および落雷予報方法 |
CN102590651A (zh) * | 2011-12-31 | 2012-07-18 | 国网电力科学研究院 | 基于雷电实测数据的输电线路故障概率评估方法 |
Non-Patent Citations (1)
Title |
---|
谢云云,薛禹胜,王昊昊: "《电网雷击故障概率的时空在线预警》", 《电力系统自动化》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105067904B (zh) * | 2015-07-17 | 2018-02-02 | 云南电力试验研究院(集团)有限公司 | 一种基于历史雷电定位数据和电网跳闸数据的雷电预警数据评估方法 |
CN105067904A (zh) * | 2015-07-17 | 2015-11-18 | 云南电力试验研究院(集团)有限公司 | 一种基于历史雷电定位数据和电网跳闸数据的雷电预警数据评估方法 |
CN105868872A (zh) * | 2016-05-30 | 2016-08-17 | 东北大学 | 一种配电网雷电灾害故障预测方法 |
CN105868872B (zh) * | 2016-05-30 | 2019-08-06 | 东北大学 | 一种配电网雷电灾害故障预测方法 |
CN106443342A (zh) * | 2016-09-29 | 2017-02-22 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 高压直流输电线路雷击故障点的识别方法 |
CN106443342B (zh) * | 2016-09-29 | 2018-11-23 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 高压直流输电线路雷击故障点的识别方法 |
CN107633320B (zh) * | 2017-08-17 | 2021-03-02 | 广东电网有限责任公司惠州供电局 | 一种基于气象预测和风险评估的电网线路重要度评估方法 |
CN107633320A (zh) * | 2017-08-17 | 2018-01-26 | 广东电网有限责任公司惠州供电局 | 一种基于气象预测和风险评估的电网线路重要度评估方法 |
CN108693450A (zh) * | 2018-04-09 | 2018-10-23 | 国网电力科学研究院武汉南瑞有限责任公司 | 用于大跨越高杆塔的输电线路绕击闪络耐雷性能分析方法 |
CN108898258A (zh) * | 2018-07-06 | 2018-11-27 | 广州供电局有限公司 | 雷电灾害天气下电力系统连锁故障风险的分析方法与系统 |
CN109541409A (zh) * | 2018-11-30 | 2019-03-29 | 国家电网有限公司 | 基于电气几何模型的配电线路闪络率改进算法 |
CN109541409B (zh) * | 2018-11-30 | 2020-05-26 | 国家电网有限公司 | 基于电气几何模型的配电线路闪络率改进算法 |
CN110851791A (zh) * | 2019-11-05 | 2020-02-28 | 云南电网有限责任公司 | 一种雷电灾害下线路故障概率计算方法及系统 |
CN112001070A (zh) * | 2020-08-10 | 2020-11-27 | 国网山东省电力公司滨州供电公司 | 一种输电线受外部环境影响停运概率的建模方法 |
CN112001070B (zh) * | 2020-08-10 | 2023-10-03 | 国网山东省电力公司滨州供电公司 | 一种输电线受外部环境影响停运概率的建模方法 |
CN114154766A (zh) * | 2022-02-07 | 2022-03-08 | 国网江苏省电力有限公司苏州供电分公司 | 雷电动态预测下电网的易损性预警方法及系统 |
CN114154766B (zh) * | 2022-02-07 | 2022-05-10 | 国网江苏省电力有限公司苏州供电分公司 | 雷电动态预测下电网的易损性预警方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103761414B (zh) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103761414A (zh) | 输电线路的雷击停运概率分析方法 | |
CN103837769B (zh) | 一种输电线路雷害预警方法及其系统 | |
CN102435921B (zh) | 同塔双回输电线路绝缘及耐雷电冲击性能的判定方法 | |
CN103207340B (zh) | 一种输电线路雷电绕击跳闸在线预警方法 | |
CN101315400B (zh) | 基于雷电参数统计的输电线路防雷性能评估方法 | |
CN106410792B (zh) | 一种电网防雷方法及系统 | |
CN105242133B (zh) | 一种改进配电线路雷电跳闸率计算方法 | |
CN105868872B (zh) | 一种配电网雷电灾害故障预测方法 | |
CN105224797B (zh) | 一种特高压大跨越线路反击跳闸率计算方法 | |
CN102967785A (zh) | 一种高速铁路牵引网防雷性能评估方法 | |
CN105426671B (zh) | 一种雷暴天气下架空配电线路的可靠性评测方法 | |
CN110866693B (zh) | 一种基于gis模型的架空输电线路覆冰风险评估方法 | |
CN104182830A (zh) | 一种基于多维分析的配电系统可靠性薄弱环节挖掘方法 | |
CN103093097A (zh) | 基于规范切的电力系统脆弱断面识别方法 | |
CN105095668B (zh) | 基于亚洲极涡因子的电网覆冰长期预报方法 | |
CN103777094A (zh) | 一种基于层次分析法的高速铁路牵引网雷害风险评估方法 | |
CN105809574A (zh) | 一种结合地形参数判定线路走廊云地闪高危区段的方法 | |
CN105067904A (zh) | 一种基于历史雷电定位数据和电网跳闸数据的雷电预警数据评估方法 | |
CN112365038A (zh) | 一种配电线路受雷击跳闸概率的评估方法及系统 | |
CN106295871A (zh) | 气象因子与输电线路山火发生概率的关联性计算方法 | |
Yuan et al. | Research on fusion algorithm of lightning strike trip warning for mountain transmission lines | |
CN106159940A (zh) | 基于电网负荷特性分析的pmu最优布点方法 | |
Xie et al. | Extensions of power system early-warning defense schemes by integrating wide area meteorological information | |
Ciapessoni et al. | A Methodology for Resilience-oriented Planning in the Italian Transmission System | |
CN106372392A (zh) | 基于核密度估计的输电线路走廊雷害分布统计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170524 Termination date: 20171218 |
|
CF01 | Termination of patent right due to non-payment of annual fee |