CN112529398A - 一种高海拔山区风电场集电线路雷击跳闸率的估算方法 - Google Patents

一种高海拔山区风电场集电线路雷击跳闸率的估算方法 Download PDF

Info

Publication number
CN112529398A
CN112529398A CN202011419812.1A CN202011419812A CN112529398A CN 112529398 A CN112529398 A CN 112529398A CN 202011419812 A CN202011419812 A CN 202011419812A CN 112529398 A CN112529398 A CN 112529398A
Authority
CN
China
Prior art keywords
lightning
wind power
conductor
current
point
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
Application number
CN202011419812.1A
Other languages
English (en)
Other versions
CN112529398B (zh
Inventor
李国庆
屠劲林
刘庭
王森
张晓朝
杨英
叶林
王超群
蒋贲
向念文
刘吉辰
杨劲
谢伟
张伟
王启江
王寿福
陈志才
杜雨晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yunnan Branch Of Huaneng New Energy Co ltd
Hefei University of Technology
Huaneng Group Technology Innovation Center Co Ltd
Huaneng Renewables Corp Ltd
Original Assignee
Yunnan Branch Of Huaneng New Energy Co ltd
Hefei University of Technology
Huaneng Group Technology Innovation Center Co Ltd
Huaneng Renewables Corp Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yunnan Branch Of Huaneng New Energy Co ltd, Hefei University of Technology, Huaneng Group Technology Innovation Center Co Ltd, Huaneng Renewables Corp Ltd filed Critical Yunnan Branch Of Huaneng New Energy Co ltd
Priority to CN202011419812.1A priority Critical patent/CN112529398B/zh
Publication of CN112529398A publication Critical patent/CN112529398A/zh
Application granted granted Critical
Publication of CN112529398B publication Critical patent/CN112529398B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/20Administration of product repair or maintenance
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了一种高海拔山区风电场集电线路雷击跳闸率的估算方法,属于风电场集电线路雷电防护领域。所述方法首先获取风电场集电线路区域内部沿线的地形高度数据,然后采用三维电气几何模型建立风电场集电线路的绕击跳闸率计算模型和的规程法建立集电线路的反击跳闸率计算模型,最后分析高海拔山区风电场集电线路直击雷的跳闸率。本发明在计算集电线路直击雷的绕击跳闸率时,考虑了山区高海拔地形因素以及高海拔环境对绝缘子闪络的影响;此外,通过考虑地形因素下导地线高度差异建立地形对集电线路的反击跳闸率关系,更加准确的计算高风电场集电线路直击雷跳闸率的评估方法。

Description

一种高海拔山区风电场集电线路雷击跳闸率的估算方法
技术领域
本发明属于风电场集电线路雷电防护领域,具体涉及一种高海拔山区风电场集电线路直击雷跳闸率的估算方法。
背景技术
基于风力发电对风能资源利用的特殊性,山区高海拔风电场常建于高海拔、强风速、高湿度、高土壤电阻率、雷电活动频繁的自然环境中。风电场35kV集电线路作为连接风机和风电场升压站的重要输电通道,保障其安全可靠运行具有至关重要的作用。但据某风电场雷击跳闸故障信息统计报告显示,在两年内该风电场集电线路因雷击故障造成跳闸次数多达20次,给风电场造成了巨大经济损失,并对电网稳定性产生了不良影响。造成山区高海拔风电场集电线路雷击跳闸率较为严重的原因有以下两点:首先由于高海拔山区复杂的地理环境,风电场集电线路采用的杆塔高度达到20-30m,基本接近或高于110kV输电线路,其引雷能力远远大于电网35kV配电网络,然而集电线路绝缘配置水平依然按照35kV线路设计;其次则是山区高海拔风电场雷电活动较平原地区更频繁,发生雷击事故的概率更高。因此,开展高海拔山区风电场集电线路雷击跳闸率计算方法分析研究具有重要工程应用价值。
雷电直击集电线路过电压形式分为反击和绕击两种,所谓反击是指当雷电击中杆塔顶部或避雷线时,使得雷电流通过雷击点的对地电位被大大的抬升,当电位差超过了杆塔绝缘子两端的耐受电压时,便引发绝缘子闪络,发生击穿,导致导线发生接地故障,引起跳闸;所谓绕击指的是有避雷线的线路,雷电绕过避雷线(屏蔽失效)而击中导线,直接在导线上引起过电压,且该过电压超过绝缘子闪络电压,从而引起集电线路发生雷击跳闸。因此,风电场集电线路直击雷跳闸率包括反击闪络率和绕击闪络率两部分。而目前对于风电场集电线路雷击跳闸率计算分析的研究现状有以下两点:
1)目前,大多数关于风电场集电线路雷击跳闸的研究主要集中在雷击故障分析和传统防雷措施优化方面,对雷击跳闸原因多为定性讨论;同时受到高海拔山区风电场特殊的自然和地理环境,对于直击雷跳闸的计算分析多采用传统的公式法进行讨论,未能考虑到高海拔及特殊的地形因素对计算结果带来的影响。
2)建于高海拔山区的风电场集电线路由于地理环境的特殊性,在计算其雷击跳闸率时,如果采用传统的电气几何模型及公式法时,未能考虑到沿线路的地线和海拔影响因素,势必对计算结果产生影响,从而无法准确定量分析集电线路雷击跳闸率。
文献“风电场集电线路杆塔雷击原因和防雷接地优化保护措施”、“风电场防雷接地系统的相关技术”、“.风电场集电线路防雷保护的研究”和专利“一种电力系统雷击闪络风险评估方法”、“一种输电线路区段雷害风险的评估方法及系统”等大多数关于风电场集电线路雷击跳闸的研究集中在雷击故障分析和传统防雷措施优化方面,对雷击跳闸原因分析多为定性讨论;对于直击雷跳闸的计算分析多采用传统的公式法进行讨论,未能考虑到山区风电场高海拔及特殊的地形因素对集电线路雷击跳闸率计算带来的影响。
综上所述,现有技术还存在以下问题:
1、目前,未考虑整条线路的绕击闪络率随地形变化情况。
2、未考虑风电场在实际的高海拔高湿度下得绝缘子闪络电压修正情况。
因此,需要提出一种可考虑地形和海拔因素计算集电线路雷击跳闸率的方法,为风电场集电线路雷电防护及评估提供科学分析手段。
发明内容
本发明要解决的技术问题是针对我国风电场雷击跳闸率相对较高,采用传统的电气几何模型计算绕击闪络率以及采用规程法计算反击闪络率的计算方法未考虑山区地形因素对风电场直击雷跳闸率影响这一现状,从而导致计算结果与实际的历史雷击跳闸次数存在一定偏差的问题。本文以考虑地形、海拔高度和杆塔结构,并基于杆塔不同位置分别计算集电线路反击闪络率,完成对集电线路直击雷跳闸率影响因素的定量分析,并与实际风电场集电线路雷击跳闸故障数据进行对比分析,验证本文提出的雷击跳闸率计算分析方法的准确性。
本发明的目的是这样实现的,本发明提出了一种高海拔山区风电场集电线路雷击跳闸率的估算方法,包括以下步骤:
步骤1,对风电场集电线路内的杆塔进行编号、进行基本参数的采集和雷电流幅值范围的设定;
设需要估算的任一个风电场集电线路所在的区域为U,在区域U内对杆塔进行编号,设区域U中的任一个杆塔为杆塔Tm,m为区域U内的杆塔按照排列进行编号的顺序号,m为正整数;
采集区域U内的风电场集电线路的基本参数,所述基本参数包括杆塔高度尺寸参数、导线参数、避雷线参数、雷电参数和绝缘子闪络电压参数;
根据已知的观测数据,首先将需要估计的雷电流幅值范围记为雷电流幅值范围A,A=0~300kA;然后将雷电流幅值范围A等分为n个区间,令每个区间起始处的雷电流幅值作为该区间评估点的雷电流幅值,得到n个区间评估点的雷电流幅值,在n个区间中任取一个区间,并将该区间评估点的雷电流幅值记为雷电流幅值Ij,j=1,2...n;最后将n个区间评估点的雷电流幅值组成一个雷电流幅值的集合,记为集合I,I={I1,I2.......Ij......In};
步骤2,进行二级区段的设定,并进行导线对地高度和避雷线对地高度的计算;
步骤2.1,进行二级区段的设定;
将杆塔Tm与杆塔Tm+1之间的区间定义为后一级区段,在后一级区段内,将杆塔Tm与杆塔Tm+1之间的水平距离按长度1米等分得到K1个后二级区段,将该其中的任意一个记为后二级区段Bh,h=1,2…..K1;将K1个后二级区段的集合记成集合B,集合
Figure BDA0002820249150000031
令集合B内每一个后二级区段终点处的海拔高度为该后二级区段的海拔高度,对集合B内的每一个后二级区段的海拔高度进行测量,得到K1个后二级区段海拔高度的集合HB
Figure BDA0002820249150000032
其中,
Figure BDA0002820249150000033
为后二级区段Bh的海拔高度;
将杆塔Tm与杆塔Tm-1之间的区间定义为前一级区段,在前级一区段内,将杆塔Tm与杆塔Tm-1杆塔之间的水平距离按长度1米等分得到K2个前二级区段,并将其中的任意一个记为前二级区段Cγ,γ=1,2…..K2,且序号γ的增大方向为自杆塔Tm向杆塔Tm-1的方向;将K2个前二级区段的集合记成集合C,集合
Figure BDA0002820249150000041
K2≠K1;令集合C内每一个前二级区段终点处的海拔高度为该前二级区段的海拔高度,对集合C内的每一个前二级区段的海拔高度进行测量,得到前一级区段内K2个前二级区段海拔高度的集合HC
Figure BDA0002820249150000042
其中,
Figure BDA0002820249150000043
为前二级区段Cγ的海拔高度;
步骤2.2,计算后导线对地高度LBh和前导线对地高度
Figure BDA0002820249150000044
令后二级区段Bh中的导线的中点为点0,点0处导线的对地垂直高度记为后导线对地高度LBh,后导线对地高度LBh的计算式如下:
Figure BDA0002820249150000045
式中,Hm为杆塔Tm的海拔高度,HC为杆塔Tm上导线悬挂点的对地垂直距离,fO为点0处的导线弧垂,D1为点0到杆塔Tm的直线水平距离,α1为杆塔Tm和杆塔Tm+1之间的高差角,HO为点0处对应的地面海拔高度;
其中,点0处的导线弧垂fo的计算公式如下:
Figure BDA0002820249150000046
式中:g1为比载,σ1为导线最大应力;
令前二级区段Cγ中的导线的中点为点Q,点Q处导线的对地垂直高度记为前导线对地高度
Figure BDA0002820249150000047
前导线对地高度
Figure BDA0002820249150000048
的计算式如下:
Figure BDA0002820249150000049
式中,fQ为点Q处的导线弧垂,D2为点Q到杆塔Tm的直线水平距离,α2为杆塔Tm和杆塔编号Tm-1之间的高差角,HQ为点Q对应的地面海拔高度;
其中,点Q处的导线弧垂fQ的计算公式如下:
Figure BDA00028202491500000410
步骤2.3,计算后避雷线对地高度
Figure BDA00028202491500000411
和前避雷线对地高度
Figure BDA00028202491500000412
令后二级区段Bh中的避雷线的中点为点S,点S处避雷线的对地垂直高度记为后避雷线对地高度
Figure BDA00028202491500000413
后避雷线对地高度
Figure BDA00028202491500000414
的计算公式如下:
Figure BDA00028202491500000415
式中,HG为杆塔Tm上避雷线悬挂点的对地垂直距离,fS为在点S处的避雷线弧垂,E1为点S到杆塔Tm的直线水平距离,HS为点S处对应的地面海拔高度;
其中,点S处的避雷线弧垂fS的计算公式如下:
Figure BDA0002820249150000051
式中:g2为比载,σ2为避雷线最大应力;
令前二级区段Cγ中的避雷线的中点为点X,点X处避雷线的对地垂直高度记为前避雷线对地高度
Figure BDA0002820249150000052
前避雷线对地高度
Figure BDA0002820249150000053
的计算式如下:
G=Hm+HG-fX-E2×tanα2-HX
式中,fX为点X处的避雷线弧垂,E2为点X到杆塔Tm的直线水平距离,HX为点X对应的地面海拔高度;
其中,点X处的避雷线弧垂fX的计算公式如下:
Figure BDA0002820249150000054
步骤3,考虑海拔高度的影响,对杆塔Tm的绝缘子雷电冲击50%放电电压U50%进行修正,得到修正过的绝缘子雷电冲击50%放电电压U150%
Figure BDA0002820249150000055
式中,Ka为修正系数,
Figure BDA0002820249150000056
q为海拔高度修正因子,本发明中,取修正因子为1.0;
步骤4,计算后二级区段Bh的风电场集电线路最小反击耐雷水平
Figure BDA0002820249150000057
和风电场集电线路最小绕击耐雷水平
Figure BDA0002820249150000058
和风电场集电线路最大绕击耐雷水平
Figure BDA0002820249150000059
计算前二级区段Cγ的风电场集电线路最小反击耐雷水平
Figure BDA00028202491500000510
和风电场集电线路最小绕击耐雷水平
Figure BDA00028202491500000511
和风电场集电线路最大绕击耐雷水平
Figure BDA00028202491500000512
步骤4.1,计算后二级区段Bh的风电场集电线路最小反击耐雷水平
Figure BDA00028202491500000513
和前二级区段Cγ中的风电场集电线路最小反击耐雷水平
Figure BDA00028202491500000514
计算式分别如下:
Figure BDA0002820249150000061
Figure BDA0002820249150000062
式中,km为考虑冲击电晕影响的耦合系数,βm为杆塔Tm的杆塔分流系数;Rm为杆塔Tm的冲击接地电阻;
Figure BDA0002820249150000063
为杆塔Tm横担高度;Vm为杆塔总电感;
Figure BDA0002820249150000064
为导线、地线间的几何耦合系数;
步骤4.2,计算后二级区段Bh的风电场集电线路最小绕击耐雷水平
Figure BDA0002820249150000065
和风电场集电线路最大绕击耐雷水平
Figure BDA0002820249150000066
计算式如下:
Figure BDA0002820249150000067
Figure BDA0002820249150000068
其中,
Figure BDA0002820249150000069
为后二级区段的风电场集电线路最大击矩,计算式如下:
Figure BDA00028202491500000610
式中:
θ01是后二级区段Bh的保护角,
Figure BDA00028202491500000611
L1为从杆塔Tm的导线到避雷线的水平横担距离;
F是第一中间变量,F=kh 2-sin201),其中,kh是在后二级区段Bh的击距系数,计算式如下:
Figure BDA00028202491500000612
Figure BDA00028202491500000613
kh=0.55;
G1是第二中间变量,
Figure BDA00028202491500000614
步骤4.3,计算前二级区段Cγ的风电场集电线路最小绕击耐雷水平
Figure BDA00028202491500000615
和风电场集电线路最大绕击耐雷水平
Figure BDA00028202491500000616
计算式如下:
Figure BDA00028202491500000617
Figure BDA0002820249150000071
其中,
Figure BDA0002820249150000072
为前二级区段的风电场集电线路最大击矩,计算式如下:
Figure BDA0002820249150000073
式中:
θ02是在前二级区段Cγ的保护角,
Figure BDA0002820249150000074
F2是第三中间变量,F2=kγ 2-sin202),kγ是在前二级区段Cγ的击距系数,,计算式如下:
Figure BDA0002820249150000075
Figure BDA0002820249150000076
kγ=0.55;
G2是第四中间变量,
Figure BDA0002820249150000077
步骤5,计算风电场集电线路直击雷反击闪络率SBFOR
步骤5.1,计算与后二级区段Bh、雷电流幅值范围A3对应的风电场集电线路发生雷电反击闪络次数BFORBhA3
从雷电流幅值范围A中分解出雷电流幅值为
Figure BDA0002820249150000078
的雷电流幅值范围A3,设雷电流幅值范围A3内共包含n3个区间,n3≤n;对雷电流幅值范围A3,逐个区间计算与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh,得到n3个与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh,对n3个与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh进行累加得到与后二级区段Bh、雷电流幅值范围A3对应的风电场集电线路发生雷电反击闪络次数BFORBhA3,计算式分别如下:
ΔBFORBh=ηgNg[P(Ij)-P(Ij+1)]
Figure BDA0002820249150000079
式中:
g为击杆率,Ng为地面落雷密度,
Figure BDA00028202491500000710
Td为年雷暴日数,η为建弧率;
P()为雷电流幅值分布概率函数,
Figure BDA0002820249150000081
即雷电流幅值超过Ij的雷电流出现的概率,P(Ij+1)为雷电流幅值超过(Ij+1)的雷电流出现的概率;
步骤5.2,计算与前二级区段Cγ、雷电流幅值范围A4对应的风电场集电系统发生雷电反击闪络次数BFORCγA4
从雷电流幅值范围A中分解出雷电流幅值为
Figure BDA0002820249150000082
的雷电流幅值范围A4,设雷电流幅值范围A4内共包含n4个区间,n4≤n;对雷电流幅值范围A4,逐个区间计算与前二级区段Cγ对应的集电线路发生雷电反击闪络次数ΔBFOR,得到n4个与前二级区段Cγ对应的集电线路发生雷电反击闪络次数ΔBFOR,对n4个与前二级区段Cγ对应的集电线路发生雷电反击闪络次数ΔBFOR进行累加得到与前二级区段Cγ、雷电流幅值范围A4对应的风电场集电线路发生雷电反击闪络次数BFORCγA4,计算式分别如下:
ΔBFOR=ηgNg[P(Ij)-P(Ij+1)]
Figure BDA0002820249150000083
步骤5.3,计算风电场集电线路直击雷反击闪络率SBFOR
在后一级区段的集合
Figure BDA0002820249150000084
中截取包括
Figure BDA0002820249150000085
的前半部分组成一个新的集合B’,
Figure BDA0002820249150000086
对集合B’中的每一个后二级区段,按照步骤5.1的方法进行计算,得到
Figure BDA0002820249150000087
个风电场集电线路发生雷电反击闪络次数BFORBhA3
在前一级区段的集合
Figure BDA0002820249150000088
中截取包括
Figure BDA0002820249150000089
的前半部分组成一个新的集合C’,
Figure BDA00028202491500000810
对C’中的每一个前二级区段,按照步骤5.2的方法进行计算,得到
Figure BDA00028202491500000811
个风电场集电线路发生雷电反击闪络次数BFORCγA4
风电场集电线路直击雷的反击闪络率SBFOR的计算式如下:
Figure BDA0002820249150000091
步骤6,计算风电场集电线路直击雷绕击闪络率SSFFOR
步骤6.1,计算与后二级区段Bh、雷电流幅值范围A1对应的风电场集电线路发生雷电绕击闪络次数SFFORBhA1
从雷电流幅值范围A中分解出雷电流幅值为
Figure BDA0002820249150000092
的雷电流幅值范围A1,设雷电流幅值范围A1内共包含n1个区间,n1≤n;对雷电流幅值范围A1,逐个区间计算与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh,得到n1个与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh,对n1个与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh进行累加得到与后二级区段Bh、雷电流幅值范围A1对应的风电场集电系统发生雷电绕击闪络次数SFFORBhA1,计算式分别如下:
ΔSFFORBh=0.1NglG1[P(Ij)-P(Ij+1)]
Figure BDA0002820249150000093
式中:
lG1为后二级区段Bh的集电线路导线的暴露弧
Figure BDA0002820249150000094
在大地上的水平投影;
步骤6.2,计算与前二级区段Cγ、雷电流幅值范围A2对应的风电场集电线路发生雷电绕击闪络次数SFFORCγA2
从雷电流幅值范围A中分解出
Figure BDA0002820249150000095
的雷电流幅值范围A2,设雷电流幅值范围A2内共包含n2个区间,n2≤n;对雷电流幅值范围A2,逐个区间计算与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数ΔSFFOR,得到n2个与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数ΔSFFOR,对n2个与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数ΔSFFOR进行累加得到与前二级区段Cγ、雷电流幅值范围A2对应的风电场集电线路发生雷电绕击闪络次数SFFORCγA2
ΔSFFOR=0.1NglG2[P(Ij)-P(Ij+1)]
Figure BDA0002820249150000101
式中:
lG2为前二级区段Cγ的集电线路导线的暴露弧
Figure BDA0002820249150000102
在大地上的水平投影;
步骤6.3,计算风电场集电线路直击雷的绕击闪络率SSFFOR
对步骤5.3中定义的集合B’中的每一个后二级区段,按照步骤6.1的方法进行计算,得到
Figure BDA0002820249150000103
个风电场集电线路发生雷电绕击闪络次数SFFORBhA1
对步骤5.3中定义的集合C’中的每一个前二级区段,按照步骤6.2的方法进行计算,得到
Figure BDA0002820249150000104
个风电场集电线路发生雷电绕击闪络次数SFFORCγA2
风电场集电线路直击雷的绕击闪络率SSFFOR的计算式如下:
Figure BDA0002820249150000105
步骤7,记考虑地形因素影响下的风电场集电线路雷击跳闸率为S,计算式如下;
S=SSFFOR+SBFOR
优选地,后二级区段Bh的集电线路导线的暴露弧
Figure BDA0002820249150000106
在大地上的水平投影lG1的计算式如下;
Figure BDA0002820249150000107
其中,
rc1为后二级区段Bh雷电流幅值Ij逐渐增加时的集电线路导线击距,rc1=10Ij 0.65
rg1为后二级区段Bh雷电流幅值Ij逐渐增加时的集电线路大地击距,rg1=khrc1
所述暴露弧
Figure BDA0002820249150000108
的获取方式为:以避雷线的中点为圆心、以导线击距rc1为半径做避雷线圆弧,以导线的中点为圆心、以导线击距rc1为半径做导线圆弧,以大地击距rg1为高度做水平线,使导线圆弧和水平线相交得到点a、导线圆弧与避雷线圆弧相交得到点b,点a和点b形成暴露弧
Figure BDA0002820249150000109
优选地,前二级区段Cγ的集电线路导线的暴露弧
Figure BDA0002820249150000111
在大地上的水平投影lG2的计算式如下:
Figure BDA0002820249150000112
其中,
rc2为前二级区段Cγ雷电流幅值Ij逐渐增加时的集电线路导线击距,rc2=10Ij 0.65
rg2为前二级区段Cγ雷电流幅值Ij逐渐增加时的集电线路大地击距rg2,rg2=kγrc2
所述暴露弧
Figure BDA0002820249150000113
的获取方式为:以避雷线的中点为圆心、以导线击距rc2为半径做避雷线圆弧,以导线的中点为圆心、以导线击距rc2为半径做导线圆弧,以大地击距rg2为高度做水平线,使导线圆弧和水平线相交得到点c、导线圆弧与避雷线圆弧相交得到点d,点c和点d形成暴露弧
Figure BDA0002820249150000114
与现有的估算方法相比,本发明的有益效果如下:
1、考虑了风电场所处的高海拔山区地形高度差异性,采用三维电气几何模型计算方法对风电场集电线路绕击闪络率的计算更加准确。
2、考虑了风电场所处的高海拔山区地形高度差异性,对传统的规程法计算反击跳闸率采用考虑杆塔两侧沿档距分布的集电线路实际高度进行积分再取平均值的计算方法,能够充分考虑集电线路杆塔附近的地形情况对反击闪络率的影响,从而使得集电线路反击闪络率的计算更加准确。
附图说明
图1为本发明的流程图。
图2考虑地形作用下沿档距方向的海拔高度及任一点导线对地高度示意图。
图3为本发明中考虑地形情况下的三维电气几何模型图。
具体实施方式
下面结合附图和实施例对对本发明进行详细的说明
图1是一种高海拔山区风电场集电线路雷击跳闸率的估算方法的流程图。由该图可见,本发明的评估方法包括以下步骤:
步骤1,对风电场集电线路内的杆塔进行编号、进行基本参数的采集和雷电流幅值范围的设定。
设需要估算的任一个风电场集电线路所在的区域为U,在区域U内对杆塔进行编号,设区域U中的任一个杆塔为杆塔Tm,m为区域U内的杆塔按照排列进行编号的顺序号,m为正整数。
采集区域U内的风电场集电线路的基本参数,所述基本参数包括杆塔高度尺寸参数、导线参数、避雷线参数、雷电参数和绝缘子闪络电压参数。
根据已知的观测数据,首先将需要估计的雷电流幅值范围记为雷电流幅值范围A,A=0~300kA;然后将雷电流幅值范围A等分为n个区间,令每个区间起始处的雷电流幅值作为该区间评估点的雷电流幅值,得到n个区间评估点的雷电流幅值,在n个区间中任取一个区间,并将该区间评估点的雷电流幅值记为雷电流幅值Ij,j=1,2...n;最后将n个区间评估点的雷电流幅值组成一个雷电流幅值的集合,记为集合I,I={I1,I2.......Ij....In}。在本实施例中n=1000~2000。
步骤2,进行二级区段的设定,并进行导线对地高度和避雷线对地高度的计算。即建立考虑地形情况下的三维电气模型,图3为本发明中考虑地形情况下的三维电气几何模型图。
步骤2.1,进行二级区段的设定。
将杆塔Tm与杆塔Tm+1之间的区间定义为后一级区段,在后一级区段内,将杆塔Tm与杆塔Tm+1之间的水平距离按长度1米等分得到K1个后二级区段,将该其中的任意一个记为后二级区段Bh,h=1,2…..K1;将K1个后二级区段的集合记成集合B,集合
Figure BDA0002820249150000121
令集合B内每一个后二级区段终点处的海拔高度为该后二级区段的海拔高度,对集合B内的每一个后二级区段的海拔高度进行测量,得到K1个后二级区段海拔高度的集合HB
Figure BDA0002820249150000122
其中,
Figure BDA0002820249150000123
为后二级区段Bh的海拔高度。
将杆塔Tm与杆塔Tm-1之间的区间定义为前一级区段,在前级一区段内,将杆塔Tm与杆塔Tm-1杆塔之间的水平距离按长度1米等分得到K2个前二级区段,并将其中的任意一个记为前二级区段Cγ,γ=1,2…..K2,且序号γ的增大方向为自杆塔Tm向杆塔Tm-1的方向。将K2个前二级区段的集合记成集合C,集合
Figure BDA0002820249150000124
K2≠K1;令集合C内每一个前二级区段终点处的海拔高度为该前二级区段的海拔高度,对集合C内的每一个前二级区段的海拔高度进行测量,得到前一级区段内K2个前二级区段海拔高度的集合HC
Figure BDA0002820249150000131
其中,
Figure BDA0002820249150000132
为前二级区段Cγ的海拔高度。
步骤2.2,计算后导线对地高度LBh和前导线对地高度
Figure BDA0002820249150000133
令后二级区段Bh中的导线的中点为点0,点0处导线的对地垂直高度记为后导线对地高度LBh,后导线对地高度LBh的计算式如下:
Figure BDA0002820249150000134
式中,Hm为杆塔Tm的海拔高度,HC为杆塔Tm上导线悬挂点的对地垂直距离,fo为点0处的导线弧垂,D1为点0到杆塔Tm的直线水平距离,α1为杆塔Tm和杆塔Tm+1之间的高差角,H0为点O处对应的地面海拔高度。
其中,点0处的导线弧垂fO的计算公式如下:
Figure BDA0002820249150000135
式中:g1为比载,σ1为导线最大应力。
令前二级区段Cγ中的导线的中点为点Q,点Q处导线的对地垂直高度记为前导线对地高度
Figure BDA0002820249150000136
前导线对地高度
Figure BDA0002820249150000137
的计算式如下:
Figure BDA0002820249150000138
式中,fQ为点Q处的导线弧垂,D2为点Q到杆塔Tm的直线水平距离,α2为杆塔Tm和杆塔编号Tm-1之间的高差角,HQ为点Q对应的地面海拔高度。
其中,点Q处的导线弧垂fQ的计算公式如下:
Figure BDA0002820249150000139
步骤2.3,计算后避雷线对地高度
Figure BDA00028202491500001310
和前避雷线对地高度
Figure BDA00028202491500001311
令后二级区段Bh中的避雷线的中点为点S,点S处避雷线的对地垂直高度记为后避雷线对地高度
Figure BDA00028202491500001312
后避雷线对地高度
Figure BDA00028202491500001313
的计算公式如下:
Figure BDA00028202491500001314
式中,HG为杆塔Tm上避雷线悬挂点的对地垂直距离,fS为在点S处的避雷线弧垂,E1为点S到杆塔Tm的直线水平距离,HS为点S处对应的地面海拔高度。
其中,点S处的避雷线弧垂fS的计算公式如下:
Figure BDA0002820249150000141
式中:g2为比载,σ2为避雷线最大应力。
令前二级区段Cγ中的避雷线中点为点X,点X处避雷线的对地垂直高度记为前避雷线对地高度
Figure BDA0002820249150000142
前避雷线对地高度
Figure BDA0002820249150000143
的计算式如下:
G=Hm+HG-fX-E2×tanα2-HX
式中,fX为点X处的避雷线弧垂,E2为点X到杆塔Tm的直线水平距离,HX为点X对应的地面海拔高度。
其中,点X处的避雷线弧垂fX的计算公式如下:
Figure BDA0002820249150000144
图2考虑地形作用下沿档距方向的海拔高度及任一点导线对地高度示意图。由图2、图3可见前一级区段和后一级区段中二级区段的分段及导线、避雷线的状态。
步骤3,考虑海拔高度的影响,对杆塔Tm的绝缘子雷电冲击50%放电电压U50%进行修正,得到修正过的绝缘子雷电冲击50%放电电压U150%
Figure BDA0002820249150000145
式中,Ka为修正系数,
Figure BDA0002820249150000146
q为海拔高度修正因子,本发明中,取q为1.0。
在本实施例中取U150%=380kV。
步骤4,计算后二级区段Bh的风电场集电线路最小反击耐雷水平
Figure BDA0002820249150000147
和风电场集电线路最小绕击耐雷水平
Figure BDA0002820249150000148
和风电场集电线路最大绕击耐雷水平
Figure BDA0002820249150000149
计算前二级区段Cγ的风电场集电线路最小反击耐雷水平
Figure BDA00028202491500001410
和风电场集电线路最小绕击耐雷水平
Figure BDA00028202491500001411
和风电场集电线路最大绕击耐雷水平
Figure BDA00028202491500001412
步骤4.1,计算后二级区段Bh的风电场集电线路最小反击耐雷水平
Figure BDA00028202491500001413
和前二级区段Cγ中的风电场集电线路最小反击耐雷水平
Figure BDA00028202491500001414
计算式分别如下:
Figure BDA0002820249150000151
Figure BDA0002820249150000152
式中,km为考虑冲击电晕影响的耦合系数,βm为杆塔Tm的杆塔分流系数;Rm为杆塔Tm的冲击接地电阻;
Figure BDA0002820249150000153
为杆塔Tm横担高度;Vm为杆塔总电感;
Figure BDA0002820249150000154
为导线、地线间的几何耦合系数。
步骤4.2,计算后二级区段Bh的风电场集电线路最小绕击耐雷水平
Figure BDA0002820249150000155
和风电场集电线路最大绕击耐雷水平
Figure BDA0002820249150000156
计算式如下:
Figure BDA0002820249150000157
Figure BDA0002820249150000158
其中,
Figure BDA0002820249150000159
为后二级区段的风电场集电线路最大击矩,计算式如下:
Figure BDA00028202491500001510
式中:
θ01是在后二级区段Bh的保护角,
Figure BDA00028202491500001511
L1为从杆塔Tm的导线到避雷线的水平横担距离。
F是第一中间变量,F=kh 2-sin20),其中,kh是在后二级区段Bh的击距系数,计算式如下:
Figure BDA00028202491500001512
Figure BDA00028202491500001513
kh=0.55;
G1是第二中间变量,
Figure BDA00028202491500001514
步骤4.3,计算前二级区段Cγ的风电场集电线路最小绕击耐雷水平
Figure BDA00028202491500001515
和风电场集电线路最大绕击耐雷水平
Figure BDA00028202491500001516
计算式如下:
Figure BDA00028202491500001517
Figure BDA0002820249150000161
其中,
Figure BDA0002820249150000162
为前二级区段的风电场集电线路最大击矩,计算式如下:
Figure BDA0002820249150000163
式中:
θ02是在前二级区段Cγ的保护角,
Figure BDA0002820249150000164
F2是第三中间变量,F2=kγ 2-sin202),kγ是在前二级区段Cγ的击距系数,,计算式如下:
Figure BDA0002820249150000165
Figure BDA0002820249150000166
kγ=0.55;
G2是第四中间变量,
Figure BDA0002820249150000167
步骤5,计算风电场集电线路直击雷反击闪络率SBFOR
步骤5.1,计算与后二级区段Bh、雷电流幅值范围A3对应的风电场集电线路发生雷电反击闪络次数BFORBhA3
从雷电流幅值范围A中分解出雷电流幅值为
Figure BDA0002820249150000168
的雷电流幅值范围A3,设雷电流幅值范围A3内共包含n3个区间,n3≤n;对雷电流幅值范围A3,逐个区间计算与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh,得到n3个与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh,对n3个与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh进行累加得到与后二级区段Bh、雷电流幅值范围A3对应的风电场集电线路发生雷电反击闪络次数BFORBhA3,计算式分别如下:
ΔBFORBh=ηgNg[P(Ij)-P(Ij+1)]
Figure BDA0002820249150000169
式中:
g为击杆率,Ng为地面落雷密度,
Figure BDA00028202491500001610
Td为年雷暴日数,η为建弧率。在本实施例中,Td为40天,η=0.8,击杆率g的选择如下:平原为1/6,山区为1/4,本发明为高海拔山区,取g=1/4。
P()为雷电流幅值分布概率函数,
Figure BDA0002820249150000171
即雷电流幅值超过Ij的雷电流出现的概率,P(Ij+1)为雷电流幅值超过(Ij+1)的雷电流出现的概率;
步骤5.2,计算与前二级区段Cγ、雷电流幅值范围A4对应的风电场集电系统发生雷电反击闪络次数BFORCγA4
从雷电流幅值范围A中分解出雷电流幅值为
Figure BDA0002820249150000172
的雷电流幅值范围A4,设雷电流幅值范围A4内共包含n4个区间,n4≤n;对雷电流幅值范围A4,逐个区间计算与前二级区段Cγ对应的集电线路发生雷电反击闪络次数ΔBFOR,得到n4个与前二级区段Cγ对应的集电线路发生雷电反击闪络次数ΔBFOR,对n4个与前二级区段Cγ对应的集电线路发生雷电反击闪络次数ΔBFOR进行累加得到与前二级区段Cγ、雷电流幅值范围A4对应的风电场集电线路发生雷电反击闪络次数BFORCγA4,计算式分别如下:
ΔBFOR=ηgNg[P(Ij)-P(Ij+1)]
Figure BDA0002820249150000173
步骤5.3,计算风电场集电线路直击雷反击闪络率SBFOR
在后一级区段的集合
Figure BDA0002820249150000174
中截取包括
Figure BDA0002820249150000175
的前半部分组成一个新的集合B’,
Figure BDA0002820249150000176
对集合B’中的每一个后二级区段,按照步骤5.1的方法进行计算,得到
Figure BDA0002820249150000177
个风电场集电线路发生雷电反击闪络次数BFORBhA3
在前一级区段的集合
Figure BDA0002820249150000178
中截取包括
Figure BDA0002820249150000179
的前半部分组成一个新的集合C’,
Figure BDA00028202491500001710
对C’中的每一个前二级区段,按照步骤5.2的方法进行计算,得到
Figure BDA00028202491500001711
个风电场集电线路发生雷电反击闪络次数BFORCγA4
风电场集电线路直击雷的反击闪络率SBFOR的计算式如下:
Figure BDA0002820249150000181
步骤6,计算风电场集电线路直击雷绕击闪络率SSFFOR
步骤6.1,计算与后二级区段Bh、雷电流幅值范围A1对应的风电场集电线路发生雷电绕击闪络次数SFFORBhA1
从雷电流幅值范围A中分解出雷电流幅值为
Figure BDA0002820249150000182
的雷电流幅值范围A1,设雷电流幅值范围A1内共包含n1个区间,n1≤n;对雷电流幅值范围A1,逐个区间计算与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh,得到n1个与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh,对n1个与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh进行累加得到与后二级区段Bh、雷电流幅值范围A1对应的风电场集电系统发生雷电绕击闪络次数SFFORBhA1,计算式分别如下:
ΔSFFORBh=0.1NglG1[P(Ij)-P(Ij+1)]
Figure BDA0002820249150000183
式中:
lG1为后二级区段Bh的集电线路导线的暴露弧
Figure BDA0002820249150000184
在大地上的水平投影。
步骤6.2,计算与前二级区段Cγ、雷电流幅值范围A2对应的风电场集电线路发生雷电绕击闪络次数SFFORCγA2
从雷电流幅值范围A中分解出
Figure BDA0002820249150000185
的雷电流幅值范围A2,设雷电流幅值范围A2内共包含n2个区间,n2≤n;对雷电流幅值范围A2,逐个区间计算与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数ΔSFFOR,得到n2个与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数ΔSFFOR,对n2个与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数ΔSFFOR进行累加得到与前二级区段Cγ、雷电流幅值范围A2对应的风电场集电线路发生雷电绕击闪络次数SFFORCγA2
ΔSFFOR=0.1NglG2[P(Ij)-P(Ij+1)]
Figure BDA0002820249150000191
式中:
lG2为前二级区段Cγ的集电线路导线的暴露弧
Figure BDA0002820249150000192
在大地上的水平投影。
步骤6.3,计算风电场集电线路直击雷的绕击闪络率SSFFOR
对步骤5.3中定义的集合B’中的每一个后二级区段,按照步骤6.1的方法进行计算,得到
Figure BDA0002820249150000193
个风电场集电线路发生雷电绕击闪络次数SFFORBhA1
对步骤5.3中定义的集合C’中的每一个前二级区段,按照步骤6.2的方法进行计算,得到
Figure BDA0002820249150000194
个风电场集电线路发生雷电绕击闪络次数SFFORCγA2
风电场集电线路直击雷的绕击闪络率SSFFOR的计算式如下:
Figure BDA0002820249150000195
步骤7,记考虑地形因素影响下的风电场集电线路雷击跳闸率为S,计算式如下;
S=SSFFOR+SBFOR
在步骤6.1中的后二级区段Bh的集电线路导线的暴露弧
Figure BDA0002820249150000196
在大地上的水平投影lG1的计算式如下;
Figure BDA0002820249150000197
其中,
rc1为后二级区段Bh雷电流幅值Ij逐渐增加时的集电线路导线击距,rc1=10Ij 0.65
rg1为后二级区段Bh雷电流幅值Ij逐渐增加时的集电线路大地击距,rg1=khrc1
所述暴露弧
Figure BDA0002820249150000198
的获取方式为:以避雷线的中点为圆心、以导线击距rc1为半径做避雷线圆弧,以导线的中点为圆心、以导线击距rc1为半径做导线圆弧,以大地击距rg1为高度做水平线,使导线圆弧和水平线相交得到点a、导线圆弧与避雷线圆弧相交得到点b,点a和点b形成暴露弧
Figure BDA0002820249150000201
在步骤6.2中的前二级区段Cγ的集电线路导线的暴露弧
Figure BDA0002820249150000202
在大地上的水平投影lG2的计算式如下:
Figure BDA0002820249150000203
其中,
rc2为前二级区段Cγ雷电流幅值Ij逐渐增加时的集电线路导线击距,rc2=10Ij 0.65
rg2为前二级区段Cγ雷电流幅值Ij逐渐增加时的集电线路大地击距rg2,rg2=kγrc2
所述暴露弧
Figure BDA0002820249150000204
的获取方式为:以避雷线的中点为圆心、以导线击距rc2为半径做避雷线圆弧,以导线的中点为圆心、以导线击距rc2为半径做导线圆弧,以大地击距rg2为高度做水平线,使导线圆弧和水平线相交得到点c、导线圆弧与避雷线圆弧相交得到点d,点c和点d形成暴露弧
Figure BDA0002820249150000205
暴露弧
Figure BDA0002820249150000206
的位置见图3。
在本实施例中的计算结果为:SBFOR=4.1674次/(百公里·年),SSFFOR=1.2387次/(百公里·年),直击雷跳闸率为S=5.4061次/(百公里·年),根据风电场集电线路雷击跳闸故障信息统计数据整理并折合后,该风电场集电线路雷击跳闸率为5.0158次/百公里·年。该计算结果与实际数据相比,较为吻合
由以上可见,本发明利用在输电线路绕击闪络率计算中通过平断面定位图提取地形参数,从集电线路采用的塔形结构和所处的地形入手,编制可考虑杆塔结构及地形因素的计算程序对绕击跳闸率计算分析;然后通过公式法计算集电线路反击闪络率;最后,对集电线路面临的直击雷击故障风险进行计算分析,提出了一种高海拔山区风电场集电线路雷击跳闸率的估算方法。

Claims (3)

1.一种高海拔山区风电场集电线路雷击跳闸率的估算方法,其特征在于,包括以下步骤:
步骤1,对风电场集电线路内的杆塔进行编号、进行基本参数的采集和雷电流幅值范围的设定;
设需要估算的任一个风电场集电线路所在的区域为U,在区域U内对杆塔进行编号,设区域U中的任一个杆塔为杆塔Tm,m为区域U内的杆塔按照排列进行编号的顺序号,m为正整数;
采集区域U内的风电场集电线路的基本参数,所述基本参数包括杆塔高度尺寸参数、导线参数、避雷线参数、雷电参数和绝缘子闪络电压参数;
根据已知的观测数据,首先将需要估计的雷电流幅值范围记为雷电流幅值范围A,A=0~300kA;然后将雷电流幅值范围A等分为n个区间,令每个区间起始处的雷电流幅值作为该区间评估点的雷电流幅值,得到n个区间评估点的雷电流幅值,在n个区间中任取一个区间,并将该区间评估点的雷电流幅值记为雷电流幅值Ij,j=1,2...n;最后将n个区间评估点的雷电流幅值组成一个雷电流幅值的集合,记为集合I,I={I1,I2.......Ij.......In};
步骤2,进行二级区段的设定,并进行导线对地高度和避雷线对地高度的计算;
步骤2.1,进行二级区段的设定;
将杆塔Tm与杆塔Tm+1之间的区间定义为后一级区段,在后一级区段内,将杆塔Tm与杆塔Tm+1之间的水平距离按长度1米等分得到K1个后二级区段,将该其中的任意一个记为后二级区段Bh,h=1,2…..K1;将K1个后二级区段的集合记成集合B,集合
Figure FDA0002820249140000011
令集合B内每一个后二级区段终点处的海拔高度为该后二级区段的海拔高度,对集合B内的每一个后二级区段的海拔高度进行测量,得到K1个后二级区段海拔高度的集合HB
Figure FDA0002820249140000012
其中,
Figure FDA0002820249140000013
为后二级区段Bh的海拔高度;
将杆塔Tm与杆塔fm-1之间的区间定义为前一级区段,在前级一区段内,将杆塔Tm与杆塔Tm-1杆塔之间的水平距离按长度1米等分得到K2个前二级区段,并将其中的任意一个记为前二级区段Cγ,γ=1,2…..K2,且序号γ的增大方向为自杆塔Tm向杆塔Tm-1的方向;将K2个前二级区段的集合记成集合C,集合
Figure FDA0002820249140000021
K2≠K1;令集合C内每一个前二级区段终点处的海拔高度为该前二级区段的海拔高度,对集合C内的每一个前二级区段的海拔高度进行测量,得到前一级区段内K2个前二级区段海拔高度的集合HC
Figure FDA0002820249140000022
其中,
Figure FDA0002820249140000023
为前二级区段Cγ的海拔高度;
步骤2.2,计算后导线对地高度LBh和前导线对地高度
Figure FDA0002820249140000024
令后二级区段Bh中的导线的中点为点0,点0处导线的对地垂直高度记为后导线对地高度LBh,后导线对地高度LBh的计算式如下:
Figure FDA0002820249140000025
式中,Hm为杆塔Tm的海拔高度,HC为杆塔Tm上导线悬挂点的对地垂直距离,f0为点0处的导线弧垂,D1为点0到杆塔Tm的直线水平距离,α1为杆塔Tm和杆塔Tm+1之间的高差角,H0为点0处对应的地面海拔高度;
其中,点0处的导线弧垂f0的计算公式如下:
Figure FDA0002820249140000026
式中:g1为比载,σ1为导线最大应力;
令前二级区段Cγ中的导线的中点为点Q,点Q处导线的对地垂直高度记为前导线对地高度
Figure FDA0002820249140000027
前导线对地高度
Figure FDA0002820249140000028
的计算式如下:
Figure FDA0002820249140000029
式中,fQ为点Q处的导线弧垂,D2为点Q到杆塔Tm的直线水平距离,α2为杆塔Tm和杆塔编号Tm-1之间的高差角,HQ为点Q对应的地面海拔高度;
其中,点Q处的导线弧垂fQ的计算公式如下:
Figure FDA00028202491400000210
步骤2.3,计算后避雷线对地高度
Figure FDA00028202491400000211
和前避雷线对地高度
Figure FDA00028202491400000212
令后二级区段Bh中的避雷线的中点为点S,点S处避雷线的对地垂直高度记为后避雷线对地高度
Figure FDA00028202491400000213
后避雷线对地高度
Figure FDA00028202491400000214
的计算公式如下:
Figure FDA00028202491400000215
式中,HG为杆塔Tm上避雷线悬挂点的对地垂直距离,fS为在点S处的避雷线弧垂,E1为点S到杆塔Tm的直线水平距离,HS为点S处对应的地面海拔高度;
其中,点S处的避雷线弧垂fS的计算公式如下:
Figure FDA0002820249140000031
式中:g2为比载,σ2为避雷线最大应力;
令前二级区段Cγ中的避雷线的中点为点X,点X处避雷线的对地垂直高度记为前避雷线对地高度
Figure FDA0002820249140000032
前避雷线对地高度
Figure FDA0002820249140000033
的计算式如下:
Figure FDA0002820249140000034
式中,fX为点X处的避雷线弧垂,E2为点X到杆塔Tm的直线水平距离,HX为点X对应的地面海拔高度;
其中,点X处的避雷线弧垂fX的计算公式如下:
Figure FDA0002820249140000035
步骤3,考虑海拔高度的影响,对杆塔Tm的绝缘子雷电冲击50%放电电压U50%进行修正,得到修正过的绝缘子雷电冲击50%放电电压U150%
Figure FDA0002820249140000036
式中,Ka为修正系数,
Figure FDA0002820249140000037
q为海拔高度修正因子,本发明中,取修正因子为1.0;
步骤4,计算后二级区段Bh的风电场集电线路最小反击耐雷水平
Figure FDA0002820249140000038
和风电场集电线路最小绕击耐雷水平
Figure FDA0002820249140000039
和风电场集电线路最大绕击耐雷水平
Figure FDA00028202491400000310
计算前二级区段Cγ的风电场集电线路最小反击耐雷水平
Figure FDA00028202491400000311
和风电场集电线路最小绕击耐雷水平
Figure FDA00028202491400000312
和风电场集电线路最大绕击耐雷水平
Figure FDA00028202491400000313
步骤4.1,计算后二级区段Bh的风电场集电线路最小反击耐雷水平
Figure FDA00028202491400000314
和前二级区段Cγ中的风电场集电线路最小反击耐雷水平
Figure FDA00028202491400000315
计算式分别如下:
Figure FDA00028202491400000316
Figure FDA00028202491400000317
式中,km为考虑冲击电晕影响的耦合系数,βm为杆塔Tm的杆塔分流系数;Rm为杆塔Tm的冲击接地电阻;
Figure FDA0002820249140000041
为杆塔Tm横担高度;Vm为杆塔总电感;
Figure FDA0002820249140000042
为导线、地线间的几何耦合系数;
步骤4.2,计算后二级区段Bh的风电场集电线路最小绕击耐雷水平
Figure FDA0002820249140000043
和风电场集电线路最大绕击耐雷水平
Figure FDA0002820249140000044
计算式如下:
Figure FDA0002820249140000045
Figure FDA0002820249140000046
其中,
Figure FDA0002820249140000047
为后二级区段的风电场集电线路最大击矩,计算式如下:
Figure FDA0002820249140000048
式中:
θ01是后二级区段Bh的保护角,
Figure FDA0002820249140000049
L1为从杆塔Tm的导线到避雷线的水平横担距离;
F是第一中间变量,F=kh 2-sin201),其中,kh是在后二级区段Bh的击距系数,计算式如下:
Figure FDA00028202491400000410
Figure FDA00028202491400000411
G1是第二中间变量,
Figure FDA00028202491400000412
步骤4.3,计算前二级区段Cγ的风电场集电线路最小绕击耐雷水平
Figure FDA00028202491400000413
和风电场集电线路最大绕击耐雷水平
Figure FDA00028202491400000414
计算式如下:
Figure FDA00028202491400000415
Figure FDA00028202491400000416
其中,
Figure FDA00028202491400000417
为前二级区段的风电场集电线路最大击矩,计算式如下:
Figure FDA0002820249140000051
式中:
θ02是在前二级区段Cγ的保护角,
Figure FDA0002820249140000052
F2是第三中间变量,F2=kγ 2-sin202),kγ是在前二级区段Cγ的击距系数,,计算式如下:
Figure FDA0002820249140000053
Figure FDA0002820249140000054
kγ=0.55;
G2是第四中间变量,
Figure FDA0002820249140000055
步骤5,计算风电场集电线路直击雷反击闪络率SBFOR
步骤5.1,计算与后二级区段Bh、雷电流幅值范围A3对应的风电场集电线路发生雷电反击闪络次数BFORBhA3
从雷电流幅值范围A中分解出雷电流幅值为
Figure FDA0002820249140000059
的雷电流幅值范围A3,设雷电流幅值范围A3内共包含n3个区间,n3≤n;对雷电流幅值范围A3,逐个区间计算与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh,得到n3个与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh,对n3个与后二级区段Bh对应的集电线路发生雷电反击闪络次数ΔBFORBh进行累加得到与后二级区段Bh、雷电流幅值范围A3对应的风电场集电线路发生雷电反击闪络次数BFORBhA3,计算式分别如下:
ΔBFORBh=ηgNg[P(Ij)-P(Ij+1)]
Figure FDA0002820249140000056
式中:
g为击杆率,Ng为地面落雷密度,
Figure FDA0002820249140000057
Td为年雷暴日数,η为建弧率;
P()为雷电流幅值分布概率函数,
Figure FDA0002820249140000058
即雷电流幅值超过Ij的雷电流出现的概率,P(Ij+1)为雷电流幅值超过(Ij+1)的雷电流出现的概率;
步骤5.2,计算与前二级区段Cγ、雷电流幅值范围A4对应的风电场集电系统发生雷电反击闪络次数
Figure FDA00028202491400000614
从雷电流幅值范围A中分解出雷电流幅值为
Figure FDA0002820249140000061
的雷电流幅值范围A4,设雷电流幅值范围A4内共包含n4个区间,n4≤n;对雷电流幅值范围A4,逐个区间计算与前二级区段Cγ对应的集电线路发生雷电反击闪络次数
Figure FDA00028202491400000615
得到n4个与前二级区段Cγ对应的集电线路发生雷电反击闪络次数
Figure FDA00028202491400000616
对n4个与前二级区段Cγ对应的集电线路发生雷电反击闪络次数
Figure FDA00028202491400000617
进行累加得到与前二级区段Cγ、雷电流幅值范围A4对应的风电场集电线路发生雷电反击闪络次数
Figure FDA00028202491400000613
计算式分别如下:
Figure FDA0002820249140000062
Figure FDA0002820249140000063
步骤5.3,计算风电场集电线路直击雷反击闪络率SBFOR
在后一级区段的集合
Figure FDA0002820249140000064
中截取包括
Figure FDA0002820249140000065
的前半部分组成一个新的集合B’,
Figure FDA0002820249140000066
对集合B’中的每一个后二级区段,按照步骤5.1的方法进行计算,得到
Figure FDA0002820249140000067
个风电场集电线路发生雷电反击闪络次数BFORBhA3
在前一级区段的集合
Figure FDA0002820249140000068
中截取包括
Figure FDA00028202491400000618
的前半部分组成一个新的集合C’,
Figure FDA0002820249140000069
对C’中的每一个前二级区段,按照步骤5.2的方法进行计算,得到
Figure FDA00028202491400000610
个风电场集电线路发生雷电反击闪络次数
Figure FDA00028202491400000611
风电场集电线路直击雷的反击闪络率SBFOR的计算式如下:
Figure FDA00028202491400000612
步骤6,计算风电场集电线路直击雷绕击闪络率SSFFOR
步骤6.1,计算与后二级区段Bh、雷电流幅值范围A1对应的风电场集电线路发生雷电绕击闪络次数SFFORBhA1
从雷电流幅值范围A中分解出雷电流幅值为
Figure FDA0002820249140000071
的雷电流幅值范围A1,设雷电流幅值范围A1内共包含n1个区间,n1≤n;对雷电流幅值范围A1,逐个区间计算与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh,得到n1个与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh,对n1个与后二级区段Bh对应的集电线路发生雷电绕击闪络次数ΔSFFORBh进行累加得到与后二级区段Bh、雷电流幅值范围A1对应的风电场集电系统发生雷电绕击闪络次数SFFORBhA1,计算式分别如下:
ΔSFFORBh=0.1NglG1[P(Ij)-P(Ij+1)]
Figure FDA0002820249140000072
式中:
lG1为后二级区段Bh的集电线路导线的暴露弧
Figure FDA0002820249140000073
在大地上的水平投影;
步骤6.2,计算与前二级区段Cγ、雷电流幅值范围A2对应的风电场集电线路发生雷电绕击闪络次数SFFORCγA2
从雷电流幅值范围A中分解出
Figure FDA0002820249140000074
的雷电流幅值范围A2,设雷电流幅值范围A2内共包含n2个区间,n2≤n;对雷电流幅值范围A2,逐个区间计算与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数
Figure FDA0002820249140000075
得到n2个与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数
Figure FDA0002820249140000076
对n2个与前二级区段Cγ对应的集电线路发生雷电绕击闪络次数
Figure FDA0002820249140000077
进行累加得到与前二级区段Cγ、雷电流幅值范围A2对应的风电场集电线路发生雷电绕击闪络次数
Figure FDA0002820249140000078
Figure FDA0002820249140000079
Figure FDA00028202491400000710
式中:
lG2为前二级区段Cγ的集电线路导线的暴露弧
Figure FDA00028202491400000711
在大地上的水平投影;
步骤6.3,计算风电场集电线路直击雷的绕击闪络率SSFFOR
对步骤5.3中定义的集合B’中的每一个后二级区段,按照步骤6.1的方法进行计算,得到
Figure FDA0002820249140000081
个风电场集电线路发生雷电绕击闪络次数SFFORBhA1
对步骤5.3中定义的集合C’中的每一个前二级区段,按照步骤6.2的方法进行计算,得到
Figure FDA0002820249140000082
个风电场集电线路发生雷电绕击闪络次数
Figure FDA0002820249140000083
风电场集电线路直击雷的绕击闪络率SSFFOR的计算式如下:
Figure FDA0002820249140000084
步骤7,记考虑地形因素影响下的风电场集电线路雷击跳闸率为S,计算式如下;
S=SSFFOR+SBFOR
2.根据权利要求1所述的一种高海拔山区风电场集电线路雷击跳闸率的估算方法,其特征在于,后二级区段Bh的集电线路导线的暴露弧
Figure FDA0002820249140000085
在大地上的水平投影lG1的计算式如下;
Figure FDA0002820249140000086
其中,
rc1为后二级区段Bh雷电流幅值Ij逐渐增加时的集电线路导线击距,rc1=10Ij 0.65
rh1为后二级区段Bh雷电流幅值Ij逐渐增加时的集电线路大地击距,rg1=khrc1
所述暴露弧
Figure FDA0002820249140000087
的获取方式为:以避雷线的中点为圆心、以导线击距rc1为半径做避雷线圆弧,以导线的中点为圆心、以导线击距rc1为半径做导线圆弧,以大地击距rg1为高度做水平线,使导线圆弧和水平线相交得到点a、导线圆弧与避雷线圆弧相交得到点b,点a和点b形成暴露弧
Figure FDA0002820249140000088
3.根据权利要求1所述的一种高海拔山区风电场集电线路雷击跳闸率的估算方法,其特征在于,前二级区段Cγ的集电线路导线的暴露弧
Figure FDA0002820249140000089
在大地上的水平投影lG2的计算式如下:
Figure FDA0002820249140000091
其中,
rc2为前二级区段Cγ雷电流幅值Ij逐渐增加时的集电线路导线击距,rc2=10Ij 0.65
rg2为前二级区段Cγ雷电流幅值Ij逐渐增加时的集电线路大地击距rg2,rg2=kγrc2
所述暴露弧
Figure FDA0002820249140000092
的获取方式为:以避雷线的中点为圆心、以导线击距rc2为半径做避雷线圆弧,以导线的中点为圆心、以导线击距rc2为半径做导线圆弧,以大地击距rg2为高度做水平线,使导线圆弧和水平线相交得到点c、导线圆弧与避雷线圆弧相交得到点d,点c和点d形成暴露弧
Figure FDA0002820249140000093
CN202011419812.1A 2020-12-07 2020-12-07 一种高海拔山区风电场集电线路雷击跳闸率的估算方法 Active CN112529398B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011419812.1A CN112529398B (zh) 2020-12-07 2020-12-07 一种高海拔山区风电场集电线路雷击跳闸率的估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011419812.1A CN112529398B (zh) 2020-12-07 2020-12-07 一种高海拔山区风电场集电线路雷击跳闸率的估算方法

Publications (2)

Publication Number Publication Date
CN112529398A true CN112529398A (zh) 2021-03-19
CN112529398B CN112529398B (zh) 2023-10-03

Family

ID=74998042

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011419812.1A Active CN112529398B (zh) 2020-12-07 2020-12-07 一种高海拔山区风电场集电线路雷击跳闸率的估算方法

Country Status (1)

Country Link
CN (1) CN112529398B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116070794A (zh) * 2023-03-29 2023-05-05 合肥工业大学 一种集电线路反击跳闸概率预测、告警方法和系统
CN116151516A (zh) * 2023-01-03 2023-05-23 湖北能源集团新能源发展有限公司 一种基于混合决策模型的新能源外送线路防雷性能评估方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109460602A (zh) * 2018-08-03 2019-03-12 山东农业大学 基于大数据和神经网络的架空线路雷击跳闸率计算方法
CN110309527A (zh) * 2019-04-08 2019-10-08 国网电力科学研究院武汉南瑞有限责任公司 一种基于电气几何模型的架空配电线路雷害风险评估方法
CN111751638A (zh) * 2020-05-12 2020-10-09 南方电网科学研究院有限责任公司 输电线路雷击跳闸率的分析方法、装置、设备及存储介质

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109460602A (zh) * 2018-08-03 2019-03-12 山东农业大学 基于大数据和神经网络的架空线路雷击跳闸率计算方法
CN110309527A (zh) * 2019-04-08 2019-10-08 国网电力科学研究院武汉南瑞有限责任公司 一种基于电气几何模型的架空配电线路雷害风险评估方法
CN111751638A (zh) * 2020-05-12 2020-10-09 南方电网科学研究院有限责任公司 输电线路雷击跳闸率的分析方法、装置、设备及存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周唯;陈兴福;任伟;: "超高压及特高压直流输电线路雷击闪络研究", 四川电力技术, no. 02 *
朱雪凌;张娟;平增;张伟;: "山区架空输电线路绕击跳闸率的计算研究", 科学技术与工程, no. 35 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116151516A (zh) * 2023-01-03 2023-05-23 湖北能源集团新能源发展有限公司 一种基于混合决策模型的新能源外送线路防雷性能评估方法
CN116151516B (zh) * 2023-01-03 2024-04-02 湖北能源集团荆门象河风电有限公司 一种基于混合决策模型的新能源外送线路防雷性能评估方法
CN116070794A (zh) * 2023-03-29 2023-05-05 合肥工业大学 一种集电线路反击跳闸概率预测、告警方法和系统
CN116070794B (zh) * 2023-03-29 2023-06-27 合肥工业大学 一种集电线路反击跳闸概率预测、告警方法和系统

Also Published As

Publication number Publication date
CN112529398B (zh) 2023-10-03

Similar Documents

Publication Publication Date Title
Chowdhuri et al. Parameters of lightning strokes: A review
WO2021109633A1 (zh) 一种基于粒子群算法的输电线路绕击跳闸率测评方法
CN105242133B (zh) 一种改进配电线路雷电跳闸率计算方法
CN110309527B (zh) 一种基于电气几何模型的架空配电线路雷害风险评估方法
WO2022127616A1 (zh) 一种输电线路全线雷电绕击跳闸风险多维分析方法
CN112529398B (zh) 一种高海拔山区风电场集电线路雷击跳闸率的估算方法
CN109753703A (zh) 一种配电线路防雷等级评估方法
CN110879333B (zh) 一种高压直流输电线路雷击重启率计算和评估方法
CN108009351A (zh) 雷击跳闸风险的分布图绘制方法
CN110865266A (zh) 一种十字型接地装置的输电线路耐雷水平试验方法
CN113011099A (zh) 一种对输电线路雷击跳闸率进行计算修正的方法
WO2021109632A1 (zh) 一种110kV输电线路绕击跳闸率的评估方法
CN105929264B (zh) 一种750kV输电线路雷电绕击性能评估方法
CN117291418A (zh) 一种基于分布式监测数据的线路雷害风险评估方法及装置
CN110865268B (zh) 低土壤电阻率地区输电杆塔雷击跳闸率测试方法
CN110865270B (zh) 雷击下220kV输电线路反击跳闸率测试方法
CN113312682B (zh) 考虑风力发电机屏蔽作用的集电线路绕击闪络率计算方法
CN110687371A (zh) 一种用于确定同塔多回线路的雷电绕击性能的方法及系统
CN110865265B (zh) 一种山地地区输电线路反击跳闸率测试方法
CN112149946A (zh) 一种输电线路避雷器位置规划方法及系统
Lihong et al. Parameters for wind caused overhead transmission line swing and fault
CN110879329A (zh) 一种适用于高海拔景观塔的防雷优化与故障定位方法
CN109888711B (zh) 500kV输电线路杆塔接地状态的分析运维方法
CN111062609A (zh) 一种输电线路杆塔雷害风险评估方法
CN109460602A (zh) 基于大数据和神经网络的架空线路雷击跳闸率计算方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant