台风引发输电线路故障的概率评估方法

1.台风引发输电线路故障的概率评估方法，包括下列步骤：
1）在控制中心汇总实时的台风预报和实况信息、一般气象预报和实况信息、实时的电网工况信息；
2）采用数学建模技术和专家系统技术描述沿线路的地理特征和周边环境特征及其对于风速、风向的影响，根据地理特征和周边环境特征将线路分段，具有相同地理特征和周边环境特征的线路为一段；
3）对于每段线路，将气象、台风与电网工况信息结合各线路的静态参数，评估其各类动态变量，包括：处于大风圈的线路的长度、持续时间、所受平均风速、所受垂直风速，线路平均降雨量系数；
4）根据每段线路的地理特征和周边环境特征结合风的方向计算线路和杆塔上所受实际风速、风向值；
5）对于每段线路，将处于大风圈线路的长度、持续时间、所受平均风速相乘，得到台风量能指标，将地形和周边环境系数、雨量系数相乘，得到线路环境系数，采用分段非线性函数表示的异物挂线概率模型计算该段线路发生异物挂线的概率，根据各段线路的异物挂线概率计算整条线路的异物挂线概率；
6）以杆塔为单位，计算杆塔两侧线路段上以杆塔处地形系数修正后的垂直风速值与设计风速的比值；杆塔雨量系数、档距系数与绝缘子相关系数相乘，得到辅助系数，采用分段非线性函数表示的杆塔风偏闪络概率模型计算杆塔发生风偏闪络的概率，综合各杆塔的风偏概率计算出整条线路的风偏闪络概率；
7）对于每段线路，计算该段线路垂直风速修正值与设计风速的比值，将处于大风圈线路的长度、持续时间、所受平均风速规格化后与雨量系数相乘得到历史环境系数，采用分段非线性函数表示的断线概率模型计算该段线路发生断线的概率，根据各段线路的断线概率计算整条线路发生断线的概率；
8）以杆塔为单位，计算杆塔两侧线路垂直风速修正值与设计风速的比值，将杆塔类型系数、线路服役时间系数相乘得到杆塔动态脆弱度指标，采用分段非线性函数表示的倒塔概率模型计算发生倒塔的概率，基于各杆塔的倒塔概率推算整条线路的倒塔概率；
9）基于四种故障概率推算得到整条线路总的故障概率；
10）根据采集装置记录的相关信息，保存运行过程中线路跳闸信息，宏观的气象信息，自动气象装置采集的局部风速、风力、风向值及其它微气象信息，导线、杆塔或绝缘子的动态数据；
11）按线路故障概率的大小，对所有线路进行排序；
筛选出线路故障概率比较大的线路，为广域测量分析保护控制系统提供风险设备集。

台风引发输电线路故障的概率评估方法\n技术领域\n[0001] 本发明属电力系统及其自动化技术领域。\n背景技术\n[0002] 输电线路的安全可靠运行与气象、环境密切相关。首先，由于受到线路走廊环境制约等原因影响，线路设计日益紧凑，精细化，易受到气象环境的影响；其次，近年来气候时有恶化，对线路运行预防恶劣气象条件提出新的要求；另外，智能电网对输电线路的运行风险管理提出了新的需要。因此，有必要研究台风等外部灾害对于输电线路的影响，评估线路故障概率，加强输电线路的气象环境风险管理，为电网调度运行提供必要的技术支撑。\n[0003] 对于台风引发输电线路故障概率的评估涉及线路内外部各种因素。相关因素的数量众多且机理复杂，很多因素无法精确确定(如局部微气象等)甚至无法近似地确定(如杆塔的基础、各构件的强度等)，因而用于评估线路风荷载与线路强度的模型、参数的有效性也无法保证。因此，这类解析方法，如试图通过实时风荷载模型分析线路故障概率的方法，很难应用于实际的输电线路故障概率评估。\n[0004] 基于历史数据的拟合方法，则过于依赖于历史数据及其准确性，而目前难以获取可靠的历史数据，因此该方法在相当长的时间内无法应用于实际工程。\n[0005] 采用本说明书所述的方法进行评估，则不依赖于历史数据的大量积累，也不必进行复杂且可信度不高的风荷载以及杆塔强度计算，能够考虑各种因素(特别是无法精确量化的因素)，具有很强的适应性。\n[0006] 近年来，自然灾害发生的频率逐步增高，对电网安全稳定的威胁越来越大，但目前还没有满足在线安全分析与预警要求的概率量化评估技术，因此急需进行相关技术的研究开发。\n发明内容\n[0007] 本发明的目的是：克服评估风荷载与杆塔强度时模型、参数不可靠的缺点，以及数据拟合方法过于依赖历史数据的缺点，提供一种能够考虑各种无法精确计算因素的，能够体现无法精确描述的故障机理的，台风引发线路故障概率量化评估方法。\n[0008] 本方法考虑台风引发线路故障的四种途径，对每种途径，分别提取影响线路故障概率的指标并量化评估，根据各种途径下故障发生的机理，建立相应的概率模型，确定输入指标与各类故障概率的对应关系。最后综合各途径下的故障概率得到线路总的故障概率。\n[0009] 具体地说，本发明是采取以下的技术方案来实现的，包括下列步骤：\n[0010] 1)在控制中心汇总实时的台风预报和实况信息、一般气象预报和实况信息、实时的电网工况信息；\n[0011] 2)采用数学建模技术和专家系统技术描述沿线路的地理特征和周边环境特征及其对于风速、风向的影响，根据地理特征和周边环境特征将线路分段，具有相同地理特征和周边环境特征的线路为一段。\n[0012] 3)对于每段线路，将气象、台风与电网工况信息结合各线路的静态参数，评估其各类动态变量，包括：\n[0013] ●处于大风圈线路所受平均风速(未必垂直于线路)Va；\n[0014] ●处于大风圈线路所受的垂直风速Vv；\n[0015] ●线路处于大风圈部分的长度Lw；\n[0016] ●线路受大风圈影响的持续时间(须考虑之前时段的数据)Tk；\n[0017] ●线路平均降雨量(系数)αr；雨量系数以影响较小的小雨为基准(值为1)，根据中(1.05)、大雨(1.1)增加系数。\n[0018] 该评估方法中用到的其它静态信息包括：\n[0019] ●线路最大设计风速Vd；\n[0020] ●采用数学建模技术和专家系统技术所评估出的沿线的地形特征和周边环境特征，沿线地形(是否加强风力)、周边环境(线路交叉、森林、居住区、空旷)系数αc；\n[0021] 其中，地形系数以不影响风力的平原地形为基准(值为1)，遇风口等加强风力的地形调高系数(1～1.3)；遇削弱风力的地形调低系数(0.8～1)；周边环境系数以空旷环境为基准(值为1)，并根据沿线工业区、居民区、树木、交叉线路的数量增加系数(1～2)；\n[0022] ●杆塔两侧档距系数αd，取实际档距与基准值的比值的倒数；\n[0023] ●杆塔安全系数βt；\n[0024] ●绝缘子重量(规格化系数)αg，取为绝缘子重量与基准值之比的倒数；\n[0025] ●串长(规格化系数)αl，取为串长与基准值之比；\n[0026] ●受风表面积(规格化系数)αs，取为表面积与基准值之比；\n[0027] ●线路服役时间(系数)αy；\n[0028] ●自立还是拉线塔(系数)αzl；\n[0029] 4)根据每段线路的地理特征和周边环境特征结合风的方向计算线路和杆塔上所受实际风速、风向值。\n[0030] 对平原地形，线路的实际风速依据规程风压高度变化系数和公式计算得到。\n[0031] 对无遮挡的迎风坡、风沿山谷方向吹时山谷内的坡面，首先将吹向山坡的水平风V分解为平行于山坡走向的顺坡风Vs＝Vcosβ(β为水平风与山坡走向的夹角)和垂直于山坡走向的向坡风Vt＝V sinβ，其中Vt受到山坡的抬升作用将增大，根据风速风压高度变化系数和公式计算得到V′t；然后由V′ t得到垂直于斜面的风速Vt2＝V′ tsinα(α为山坡的地形倾角)和平行于斜面的风速Vt1＝V′ tcosα，其中Vt2受山坡的影响，在斜面上形成气流，将Vt1与斜面直接的夹角改为 ( 为气流坡度)，从而得到新的平行于斜面的风速 垂直于斜面的风速忽略；最后计算得到作用于线路的实际风速\n以及Va与山坡走向的夹角θ＝arccos(Va/Vi)。\n[0032] 对风垂直山谷走向吹时山谷中各坡面、无遮挡迎风坡的背风坡，待求点处的风速a b\n利用已知点的风速和权重函数来插值求解。其中，权重函数w(r，h)＝1/rh，r为待求点与已知点之间的距离，指数α和b为非负数，h表示已知点与待求点之间地形高度变化的总量。那么，待求点的风速 (j≠i，i，j＝1，2，Λ，n)，其中 和\n分别是待求点和已知点点风速，n为相关点的总数。\n[0033] 其它地形(峡谷风道口和山谷形状变化等)，地形系数以不影响风力的平原地形为基准(基准值为1)，风口、山谷空间减小等加强风力的地形系数为1～1.3；山谷空间增大等削弱风力的地形系数为0.8～1。\n[0034] 5)评估异物挂线导致故障的概率。对于每段线路，将该段线路的Va、Lw、Tk相乘，得到该段线路台风量能指标Ewl＝VaLwTk。将地形和周边环境系数、雨量系数相乘，成为环境系数：α＝αcαr，以分段非线性函数f(α，Ewl)模拟该段线路的发生异物挂线的概率，根据各段线路的异物挂线概率计算整条线路的异物挂线概率p1。\n[0035] 6)评估发生风偏闪络的概率。以杆塔为单位，计算线路垂直风速修正值与设计风速的比值(简称风速比)Ed＝(Vvαc)/(Vdβt)；档距系数、服役时间系数与环境系数相乘，得到辅助系数α＝αdαgαlαsαr，以分段非线性函数f(α，Ed)模拟杆塔的风偏闪络概率。\n根据各杆塔的风偏闪络概率计算线路的风偏闪络概率p2。\n[0036] 7)评估发生断线的概率。对于每段线路，计算该段线路垂直风速修正值与设计风速的比值(简称风速比)Ed＝(Vvαc)/(Vdβt)；其它历史环境系数α＝αrαWSαWDαWL，其中αWS＝Va/VB，αWL＝Lw/LB，αWD＝Tk/TB，VB、LB、TB是用于将Va、Lw、Tk规格化的基准值。\n以分段非线性函数f(α，Ed)模拟该段线路发生断线故障的概率，根据各段线路的断线概率计算整条线路发生断线的概率p3。\n[0037] 8)评估发生倒塔的概率。以杆塔为单位，计算杆塔两侧线路垂直风速修正值与设计风速的比值(简称风速比)Ed＝(αdαcVv)/(βtVd)，杆塔动态脆弱度指标Es＝αzlαy，以分段非线性函数f(Es，Ed)模拟发生倒塔故障的概率。综合各杆塔的倒塔概率得到整条线路发生倒塔的概率p4。\n[0038] 9)基于四种故障概率推算得到整条线路的故障概率p＝1-(1-p1)(1-p2)(1-p3)(1-p4)。\n[0039] 10)根据采集装置记录的相关信息，保存运行过程中线路跳闸信息，宏观的气象信息，自动气象装置采集的局部风速、风力、风向值及其它微气象信息，导线、杆塔或绝缘子的动态数据等。\n[0040] 11)按线路故障概率的大小，对所有线路进行排序。筛选出线路故障概率比较大的线路，为广域测量分析保护控制系统提供风险设备集。\n[0041] 本发明的有益效果如下：本方法定量评估风致输电线路故障的内部及外部因素，根据台风引发输电线路故障的机理建立合理的概率模型描述相关内、外部因素与线路故障的关系。使得风致线路故障概率评估不必依赖于历史数据的积累，也不必进行复杂且可信度不高的风荷载以及杆塔强度计算，且能够考虑各种相关的影响因素(特别是无法精确量化的因素)，具有很强的适应性。因此，本方法极大地推动了输电线路风险管理，提高了电力系统防御外部灾害的能力。\n附图说明\n[0042] 图1两点间地形起伏示意图。\n[0043] 图2给定值时的线路的异物挂线概率模型。\n[0044] 图3给定值时杆塔的风偏闪络概率模型。\n[0045] 图4给定值时线路的断线概率模型。\n[0046] 图5给定值时杆塔的倒塔概率模型。\n[0047] 图6评估台风引发输电线路故障概率的流程。\n具体实施方式\n[0048] 下面结合附图5，对本发明方法进行详细描述。\n[0049] 图5中步骤1描述的是在控制中心汇总实时的台风预报和实况信息、一般气象信息以及实时的电网工况信息；\n[0050] 图5中步骤2描述的是采用数学建模技术和专家系统技术描述沿线路的地理特征和周边环境特征及其对于风速、风向的影响，根据地理特征和周边环境特征将线路分段，具有相同地理特征和周边环境特征的线路为一段。\n[0051] 图5中步骤3是对于每段线路，将气象、台风与电网工况信息结合各线路的静态参数，评估其各类动态变量，包括：\n[0052] ●处于大风圈线路所受平均风速(未必垂直于线路)Va；\n[0053] ●处于大风圈线路所受的垂直风速Vv；\n[0054] ●线路处于大风圈部分的长度Lw；\n[0055] ●线路受大风圈影响的持续时间(须考虑之前时段的数据)Tk；\n[0056] ●线路平均降雨量(系数)αr；雨量系数以影响较小的小雨为基准(值为1)，根据中(1.05)、大雨(1.1)增加系数。\n[0057] ●线路实时状态信息，包括：灾害期间的跳闸信息、实时的断线倒塔信息。\n[0058] 该评估方法中用到的其它静态信息包括：\n[0059] ●线路最大设计风速Vd；\n[0060] ●采用数学建模技术和专家系统技术所评估出的沿线的地形特征和周边环境特征，沿线地形(是否加强风力)、周边环境(线路交叉、森林、居住区、空旷)系数αc；\n[0061] 其中，地形系数以不影响风力的平原地形为基准(值为1)，遇风口等加强风力的地形调高系数(1～1.3)；遇削弱风力的地形调低系数(0.8～1)；周边环境系数以空旷环境为基准(值为1)，并根据沿线工业区、居民区、树木、交叉线路的数量增加系数(1～2)；\n[0062] ●杆塔两侧档距系数αd，取实际档距与基准值的比值的倒数；\n[0063] ●杆塔安全系数βt；\n[0064] ●绝缘子重量(规格化系数)αg，取为绝缘子重量与基准值之比的倒数；\n[0065] ●串长(规格化系数)αl，取为串长与基准值之比；\n[0066] ●受风表面积(规格化系数)αs，取为表面积与基准值之比；\n[0067] ●线路服役时间αy；\n[0068] ●自立还是拉线塔(系数)αzl；\n[0069] 图5中步骤4根据地理特征和环境特征修正背景风速、风向，得到线路和杆塔上所受实际风速、风向值。\n[0070] 对平原地形，线路的实际风速依据规程风压高度变化系数和公式计算得到。\n[0071] 对无遮挡的迎风坡、风沿山谷方向吹时山谷内的坡面，首先将吹向山坡的水平风V分解为平行于山坡走向的顺坡风Vs＝Vcosβ(β为水平风与山坡走向的夹角)和垂直于山坡走向的向坡风Vt＝Vsinβ，其中Vt受到山坡的抬升作用将增大，根据风速风压高度变化系数和公式计算得到V′t；然后由V′ t得到垂直于斜面的风速Vt2＝V′ tsinα(α为山坡的地形倾角)和平行于斜面的风速Vt1＝V′ tcosα，其中Vt2受山坡的影响，在斜面上形成气流，将Vt1与斜面直接的夹角改为 ( 为气流坡度)，从而得到新的平行于斜面的风速 垂直于斜面的风速忽略；最后计算得到作用于线路的实际风速\n以及Va与山坡走向的夹角θ＝arccos(Va/Vi)。\n[0072] 对风垂直山谷走向吹时山谷中各坡面、无遮挡迎风坡的背风坡，待求点处的风速a b\n根据已知点的风速和权重函数来插值求解。其中，权重函数w(r，h)＝1/rh，r为待求点与已知点之间的距离，指数α和b为非负数，h表示已知点与待求点之间地形高度变化的总量。那么，待求点的风速 (j≠i，i，j＝1，2，Λ，n)，其中 和\n分别是待求点和已知点点风速，n为相关点的总数。假设点1和点6之间的地形起伏变化a b\n如图1所示，在它们之间还有4个点，则点6处的风速根据权重函数w(r，h)＝1/rh求解。\n其中，r为待求点6与已知点1～5之间的距离，h表示已知点6与待求点1～5之间地形高度变化的总量，如h15＝h12+h23+h34+h45+h56，其中hij＝|zi-zj|，zi表示i点的高程。则待求点6的风速 (j≠i，i，j＝1，2，Λ，6)，其中 和 分别是待求\n点和已知点风速。\n[0073] 其它地形(峡谷风道口和山谷形状变化等)，地形系数以不影响风力的平原地形为基准(基准值为1)，风口、山谷空间减小等加强风力的地形系数为1～1.3；山谷空间增大等削弱风力的地形系数为0.8～1。\n[0074] 图5中步骤5评估异物挂线导致故障的概率。对于每段线路，将该段线路的Va、Lw、Tk相乘，得到该段线路台风量能指标Ewl＝VaLwTk。该段线路的地形系数、周边环境系数、雨量系数相乘，成为环境系数：α＝αcαr，以式(1)所示的分段非线性函数评估该段线路的发生异物挂线的概率，根据各段线路的异物挂线概率计算整条线路的异物挂线概率p1。\n[0075] 令i＝1，2，3，4，j＝1，2，ΛΛn，以i＝1，n＝3为例，该分段非线性函数可以表示为：\n[0076] \n[0077] 对某一给定环境系数α1，该式可以图2来表示，当Ewl＜Ewl1时，异物挂线的概率缓慢增加，当Ewl1≤Ewl＜Ewl2时，异物挂线的概率增加较快，当Ewl2≤Ewl时，异物挂线的概率趋于饱和。其中，Ewl1、Ewl2可根据具体情况修正。\n[0078] 当α增加到α2时，异物挂线的概率相应增加。\n[0079] 图6中步骤6评估发生风偏闪络的概率。以杆塔为单位，计算线路垂直风速修正值与设计风速的比值(简称风速比)Ed＝(Vvαc)/(Vdβt)；档距系数、服役时间系数与环境系数相乘，得到辅助系数α＝αdαgαlαsαr，以式(2)所示的分段非线性函数评估杆塔的风偏闪络概率。\n[0080] 令j＝1，2，ΛΛn，以i＝2，n＝3为例，该分段非线性函数可以表示为：\n[0081] \n[0082] 对某一给定辅助系数α1，该式可以图3来表示，当Ed＜Ed1时，杆塔的风偏闪络概率为0，当Ed1≤Ed＜Ed2时，杆塔的风偏闪络概率增加较快，当Ed2≤Ed时，杆塔的风偏闪络概率趋于饱和。其中，Ed1、Ed2可根据具体情况修正。当α增加到α2时，风偏闪络概率相应增加。根据各杆塔的风偏闪络概率计算线路的风偏闪络概率p2。\n[0083] 图6中步骤7评估发生断线的概率。根据地形特征将线路分段，对于每段线路，计算该段线路垂直风速修正值与设计风速的比值(简称风速比)Ed＝(Vvαc)/(Vdβt)；其它历史环境系数α＝αrαWSαWDαWL，其中αWS＝Va/VB，αWL＝Lw/LB，αWD＝Tk/TB，VB、LB、TB是用于将Va、Lw、Tk规格化的基准值，以式(3)所示的分段非线性函数评估线路发生断线故障的概率p3。\n[0084] 令j＝1，2，ΛΛn，以i＝3，n＝2为例，该分段非线性函数可以表示为：\n[0085] \n[0086] 对某一给定历史环境系数α1，该式可以图4来表示，当Ed＜Ed1时，线路的断线概率为0，当Ed1≤Ed时，线路的断线概率缓慢增加。其中，Ed1可根据具体情况修正。当α增加到α2时，断线概率相应增加。\n[0087] 图6中步骤8评估发生倒塔的概率。以杆塔为单位，计算线路垂直风速修正值与设计风速的比值(简称风速比)Ed＝(αdαcVv)/(βtVd)，杆塔动态脆弱度指标Es＝αzlαy，以式(4)所示的分段非线性函数评估发生倒塔故障的概率。\n[0088] 令j＝1，2，ΛΛn，以i＝4，n＝4为例，该分段非线性函数可以表示为：\n[0089] \n[0090] 对某一给定杆塔动态脆弱度指标Es1，该式可以图5来表示，当Ed＜Ed1时，杆塔的倒塔概率为0，当Ed1≤Ed＜Ed2时，杆塔的倒塔概率缓慢增加，当Ed2≤Ed＜Ed3时，杆塔的倒塔概率增加较快，当Ed3≤Ed时，杆塔的倒塔概率趋于饱和。其中，Ed1、Ed2、Ed3可根据具体情况修正。当Es增加到Es2时，倒塔概率相应增加。综合各杆塔的倒塔概率得到整条线路发生倒塔的概率p4。\n[0091] 图6中步骤9基于四种故障概率推算得到整条线路的故障概率p＝1-(1-p1)(1-p2)(1-p3)(1-p4)。\n[0092] 图6中步骤10自动积累运行过程信息，根据采集装置记录的相关信息，保存运行过程中线路跳闸信息，宏观的气象信息，自动气象装置采集的局部风速、风力、风向值及其它微气象信息，导线、杆塔或绝缘子的动态数据等。\n[0093] 图6中步骤11为广域测量分析保护控制系统提供风险设备集。按线路故障概率的大小，对所有线路进行排序。筛选出线路故障概率大于k(k可以人为设定)的线路，将这些线路信息输入广域测量分析保护控制系统的风险设备集中。