CN107480337B - 多因素驱动架空线路故障率建模方法 - Google Patents

多因素驱动架空线路故障率建模方法 Download PDF

Info

Publication number
CN107480337B
CN107480337B CN201710571200.6A CN201710571200A CN107480337B CN 107480337 B CN107480337 B CN 107480337B CN 201710571200 A CN201710571200 A CN 201710571200A CN 107480337 B CN107480337 B CN 107480337B
Authority
CN
China
Prior art keywords
overhead line
fault rate
factor
comprehensive
calculating
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.)
Active
Application number
CN201710571200.6A
Other languages
English (en)
Other versions
CN107480337A (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.)
Zhejiang University ZJU
State Grid Corp of China SGCC
State Grid Zhejiang Electric Power Co Ltd
Shaoxing Power Supply Co of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Zhejiang University ZJU
State Grid Corp of China SGCC
State Grid Zhejiang Electric Power Co Ltd
Shaoxing Power Supply Co of State Grid Zhejiang Electric Power Co 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 Zhejiang University ZJU, State Grid Corp of China SGCC, State Grid Zhejiang Electric Power Co Ltd, Shaoxing Power Supply Co of State Grid Zhejiang Electric Power Co Ltd filed Critical Zhejiang University ZJU
Priority to CN201710571200.6A priority Critical patent/CN107480337B/zh
Publication of CN107480337A publication Critical patent/CN107480337A/zh
Application granted granted Critical
Publication of CN107480337B publication Critical patent/CN107480337B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Monitoring And Testing Of Exchanges (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Abstract

本发明公开了一种多因素驱动架空线路故障率建模方法。包括如下步骤:步骤1、根据架空线路的材料,计算得到架空线路期望寿命;步骤2、根据架空线路实际运行时间及其所处环境温度统计数据,计算得到架空线路等效服役时间;步骤3、根据所研究地区架空线路故障率与天气情况的统计数据,通过假设检验的方式得到各种天气因素与架空线路故障率的相关性,并通过加权的方式得到天气综合状况评分值;步骤4、根据架空线路的基准故障率函数、天气综合状况评分值、线路健康状态值和负载率得到综合故障率函数;步骤5、得到适用于该地区的多因素驱动架空线路故障率模型。本发明可以科学准确地给出架空线路的故障率。

Description

多因素驱动架空线路故障率建模方法
技术领域
本发明涉及电力分析技术,尤其是电力系统架空线路故障率建模领域。
背景技术
近年来,随着我国电网的快速发展,特别是特高压建设的不断推进落地,与之配套的大量新建、改扩建电力工程的不断实施,以及各级电网规模的不断扩大,电网运行的风险隐患逐渐增多,电网安全面临着更为严峻的考验。如何有效提高电网风险管控能力、提升电网安全可靠供电水平,成为电网安全管理工作的重中之重。架空线路作为电力传输的重要通道,在电网安全中占据着重要地位,由于其运行环境复杂多变,极易受到内外部因素的共同影响而发生故障停运,严重影响人民的正常生活秩序。因此,建立一个可以全面量化内外部影响因素的架空线故障率模型,对于了解电力系统运行风险、提高输电可靠性具有重要意义。
目前存在的针对电力系统架空线路故障率建模的方法,或是考虑单一的影响因素,不能全面反映导致架空线路故障的各影响因素之间的耦合关系;或是虽然整合多种影响因素建模的方法,但这些模型中对于影响因素的选择更多依靠直觉判断,客观性不足。
发明内容
本发明所要解决的技术问题就是提供一种多因素驱动架空线路故障率建模方法,科学准确地给出架空线路的故障率。
为解决上述技术问题,本发明采用如下技术方案:多因素驱动架空线路故障率建模方法,包括下列步骤:
步骤1、根据架空线路的材料,计算得到架空线路期望寿命;
步骤2、根据架空线路实际运行时间及其所处环境温度测量数据,计算得到架空线路等效服役时间;
步骤3、根据所研究地区架空线路故障率与天气情况的统计数据,通过假设检验的方式得到各种天气因素与架空线路故障率的相关性,并通过加权的方式得到天气综合状况评分值;
步骤4、根据架空线路的基准故障率函数、天气综合状况评分值、线路健康状态值和负载率得到综合故障率函数;
步骤5、通过所研究地区的历史统计数据,对得到的综合故障率函数进行参数估计,得到适用于该地区的多因素驱动架空线路故障率模型。
作为优选,步骤1中架空线路期望寿命的计算方法为:
Figure BDA0001349715020000021
式中:W为导线抗拉强度损失百分比;Wa为导线完全退火后抗拉强度损失百分比;Wmax为导线完成服役退出运行时抗拉强度损失百分比;θ为导线实际测量温度,K、A、B、C、D和R为对应导线的属性常数。
作为优选,步骤2中架空线路等效服役时间的计算方法为:
将架空线路实际运行时间T划分为n个小的时间区间t1,t2,…,tn,全区间恒定运行温度θ(ti)取时间区间ti初始时刻的温度,架空线路等效服役时间Teq为:
Figure BDA0001349715020000022
式中,
Figure BDA0001349715020000023
作为优选,步骤3中通过假设检验的方式得到各种天气因素与架空线路故障率的相关性,其方法为:
首先采用Pearson相关系数的概念计算影响因素与架空线故障之间的相关系数:
Figure BDA0001349715020000031
式中,xi为第i个观测样本对应影响因素归一化后数值,yi为第i个观测样本的故障率;
然后采用t分布检验法对与r相应的总体相关系数ρ进行假设检验,
假设H0:ρ=0;H1:ρ≠0,t检验法的检验统计量为:
Figure BDA0001349715020000032
根据假设检验p-的定义,计算公式如下:
Figure BDA0001349715020000033
定义影响因素的重要度水平为(1-p-),根据上述算法,给出影响因素的重要度,通过设定门槛值,选择出重要度水平高于门槛值的影响因素,将其纳入整体天气状况评价的指标。
作为优选,步骤3所述的通过加权的方式得到天气综合状况评分值,其具体方法为:
根据重要度计算各影响因素权重:
Figure BDA0001349715020000034
式中,qj为第j个影响因素权重,(1-p-j)为对应因素重要度,
得到天气综合状况评分值Z1,即:
Figure BDA0001349715020000041
作为优选,步骤4中的综合故障率函数为:
Figure BDA0001349715020000042
式中,L为架空线路期望寿命;Teq为架空线路等效服役时间;Z1为天气综合状况评分值;Z2、Z3分别代表设备健康状态评分值和架空线路负载率,β,γ123为待估计参数。
作为优选,步骤5中采用最小二乘法、极大似然估计法或Levenberg-Marquardt法对步骤4得到的模型中β,γ123参数进行估计,得到参数的估计值,带回到式(9)中,得到适用于该地区的多因素驱动架空线路故障率模型。
本发明根据所研究地区的相关历史统计数据,构建一种多因素驱动架空线路故障率模型,能够综合考虑内外部等多种因素对架空线路故障率的影响,并通过假设检验的方式避免人为主观盲目性,科学准确地给出架空线路的故障率。使得调度运行人员能够更加直观地掌握当前风险水平,辅助调度运行人员做出快速而准确的决策,提高调度运行人员对电网的管控能力,保障电网的安全稳定运行。本发明的方法可靠、易行,便于推广。
附图说明
下面结合附图和具体实施方式对本发明作进一步描述:
图1为多因素驱动架空线路故障率建模方法的流程图。
具体实施方式
请参阅图1,本发明多因素驱动架空线路故障率建模方法,包括下列步骤:
步骤1、根据架空线路的材料,计算得到架空线路期望寿命;
所述的架空线路期望寿命,其计算方法为:
Figure BDA0001349715020000051
式中:W为导线抗拉强度损失百分比,即导线损失强度与初始强度的比值;Wa为导线完全退火后抗拉强度损失百分比;Wmax为导线完成服役退出运行时抗拉强度损失百分比;θ为导线实际测量温度,K;A,B,C,D和R为对应导线的属性常数。Wmax、Wa、K、A、B、C、D和R可以根据导线的材料参照表1得到。
对应导线的属性常数均可以查表中获得,不同材质和温度下的输电线路材料参数表如表1所示。
Figure BDA0001349715020000052
表1不同材质和温度下的输电线路材料参数表
步骤2、根据架空线路实际运行时间及其所处环境温度测量数据,计算得到架空线路等效服役时间;
所述的架空线路等效服役时间,其具体计算方法为:
将架空线路实际运行时间T划分为n个小的时间区间t1,t2,…,tn,全区间恒定运行温度θ(ti)取时间区间ti初始时刻的温度(将架空线路运行时间分为了T段,每一段都有起始时刻和终止时刻,起始时刻对应的温度就是区间初始时刻温度,该温度可以通过温度传感器测量),由于沿线穿越多个不同气候区域,区间温度取全线温度最高值。架空线路等效服役时间Teq为:
Figure BDA0001349715020000061
式中,
Figure BDA0001349715020000062
步骤3、根据所研究地区架空线路故障率与天气情况的统计数据,通过假设检验的方式得到各种天气因素与架空线路故障率的相关性,并通过加权的方式得到天气综合状况评分值;
所述的通过假设检验的方式得到各种天气因素与架空线路故障率的相关性,其具体方法为:
首先采用Pearson相关系数的概念计算影响因素与架空线故障之间的相关系数:
Figure BDA0001349715020000063
式中,xi为第i个观测样本对应影响因素归一化后数值,yi为第i个观测样本的故障率。
观测样本指的是电网公司历史统计的数据,包括当时的各种天气状况和对应的架空线路故障率。各种天气因素指的是电网公司历史统计数据中统计的各种天气因素,可以包括但不限于最大风速、平均相对湿度、瞬时降雨量、平均温度等。
然后采用t分布检验法对与r相应的总体相关系数ρ进行假设检验。假设H0:ρ=0;H1:ρ≠0,t检验法的检验统计量为:
Figure BDA0001349715020000064
根据假设检验p-的定义,计算公式如下:
Figure BDA0001349715020000071
定义影响因素的重要度水平为(1-p-),根据上述算法,给出影响因素的重要度,通过设定门槛值,选择出重要度水平高于门槛值的影响因素,将其纳入整体天气状况评价的指标。
所述的通过加权的方式得到天气综合状况评分值,其具体方法为:
根据重要度计算各影响因素权重:
Figure BDA0001349715020000072
式中,qj为第j个影响因素权重,(1-p-j)为对应因素重要度。
得到天气综合状况评分值Z1,即:
Figure BDA0001349715020000073
步骤4、根据架空线路的基准故障率函数、天气综合状况评分值、线路健康状态值和负载率得到综合故障率函数;
所述的综合故障率函数,其具体形式为:
Figure BDA0001349715020000074
式中,架空线路期望寿命L为架空线路预期寿命;Teq为架空线路等效服役时间;Z1为天气综合状况评分值;Z2、Z3分别代表设备健康状态评分值和架空线路负载率,可通过查询所研究地区电力系统统计数据计算得到。β,γ123为待估计参数。
步骤5、通过所研究地区的历史统计数据,对得到的综合故障率函数进行参数估计,得到适用于该地区的多因素驱动架空线路故障率模型。
所述的对得到的综合故障率函数进行参数估计,得到适用于该地区的多因素驱动架空线路故障率模型。其具体方法为:
根据所研究地区的相关历史统计数据,采用最小二乘法、极大似然估计法或Levenberg-Marquardt法等参数估计方法,对步骤4得到的模型中β,γ123等参数进行估计,得到参数的估计值,带回到综合故障率函数中,得到适用于该地区的多因素驱动架空线路故障率模型。
除上述优选实施例外,本发明还有其他的实施方式,本领域技术人员可以根据本发明作出各种改变和变形,只要不脱离本发明的精神,均应属于本发明权利要求书中所定义的范围。

Claims (4)

1.多因素驱动架空线路故障率建模方法,其特征在于,包括下列步骤:
步骤1、根据架空线路的材料,计算得到架空线路期望寿命;
步骤2、根据架空线路实际运行时间及其所处环境温度测量数据,计算得到架空线路等效服役时间;
步骤3、根据所研究地区架空线路故障率与天气情况的统计数据,通过假设检验的方式得到各种天气因素与架空线路故障率的相关性,并通过加权的方式得到天气综合状况评分值;
步骤4、根据架空线路的基准故障率函数、天气综合状况评分值、线路健康状态值和负载率得到综合故障率函数:
Figure FDA0002536646220000011
式中,L为架空线路期望寿命;Teq为架空线路等效服役时间;Z1为天气综合状况评分值;Z2、Z3分别代表设备健康状态评分值和架空线路负载率,β,γ123为待估计参数;
步骤5、通过所研究地区的历史统计数据,对得到的综合故障率函数进行参数估计,得到适用于该地区的多因素驱动架空线路故障率模型;
步骤1中架空线路期望寿命的计算方法为:
Figure FDA0002536646220000012
式中:W为导线抗拉强度损失百分比;Wa为导线完全退火后抗拉强度损失百分比;Wmax为导线完成服役退出运行时抗拉强度损失百分比;θ为导线实际测量温度,A、B、C、D和R为对应导线的属性常数;
步骤2中架空线路等效服役时间的计算方法为:
将架空线路实际运行时间T划分为n个小的时间区间t1,t2,…,tn,全区间恒定运行温度θ(ti)取时间区间ti初始时刻的温度,架空线路等效服役时间Teq为:
Figure FDA0002536646220000021
式中,
Figure FDA0002536646220000022
2.根据权利要求1所述的多因素驱动架空线路故障率建模方法,其特征在于,步骤3中通过假设检验的方式得到各种天气因素与架空线路故障率的相关性,其方法为:
首先采用Pearson相关系数的概念计算影响因素与架空线故障之间的相关系数:
Figure FDA0002536646220000023
式中,xi为第i个观测样本对应影响因素归一化后数值,yi为第i个观测样本的故障率;
然后采用t分布检验法对与r相应的总体相关系数ρ进行假设检验,
假设H0:ρ=0;H1:ρ≠0,t检验法的检验统计量为:
Figure FDA0002536646220000024
根据假设检验p-的定义,计算公式如下:
Figure FDA0002536646220000025
定义影响因素的重要度水平为(1-p-),根据上述算法,给出影响因素的重要度,通过设定门槛值,选择出重要度水平高于门槛值的影响因素,将其纳入整体天气状况评价的指标。
3.根据权利要求2所述的多因素驱动架空线路故障率建模方法,其特征在于,步骤3所述的通过加权的方式得到天气综合状况评分值,其具体方法为:
根据重要度计算各影响因素权重:
Figure FDA0002536646220000031
式中,qj为第j个影响因素权重,(1-p-j)为对应因素重要度,
得到天气综合状况评分值Z1,即:
Figure FDA0002536646220000032
4.根据权利要求1所述的多因素驱动架空线路故障率建模方法,其特征在于,步骤5中采用最小二乘法、极大似然估计法或Levenberg-Marquardt法对步骤4得到的模型中β,γ123参数进行估计,得到参数的估计值,带回到式(9)中,得到适用于该地区的多因素驱动架空线路故障率模型。
CN201710571200.6A 2017-07-13 2017-07-13 多因素驱动架空线路故障率建模方法 Active CN107480337B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710571200.6A CN107480337B (zh) 2017-07-13 2017-07-13 多因素驱动架空线路故障率建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710571200.6A CN107480337B (zh) 2017-07-13 2017-07-13 多因素驱动架空线路故障率建模方法

Publications (2)

Publication Number Publication Date
CN107480337A CN107480337A (zh) 2017-12-15
CN107480337B true CN107480337B (zh) 2020-10-20

Family

ID=60596582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710571200.6A Active CN107480337B (zh) 2017-07-13 2017-07-13 多因素驱动架空线路故障率建模方法

Country Status (1)

Country Link
CN (1) CN107480337B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110533213A (zh) * 2019-07-08 2019-12-03 广东工业大学 基于支持向量机的输电线路缺陷风险建模及其预测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462718A (zh) * 2014-12-25 2015-03-25 国家电网公司 一种变电站经济运行年限区间评估方法
CN106096827A (zh) * 2016-06-02 2016-11-09 国网山东省电力公司济南供电公司 基于混合抽样与最小切负荷计算的可靠性计算方法及平台

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8682623B1 (en) * 2007-04-24 2014-03-25 University Of South Florida Electric power distribution interruption risk assessment calculator
US8165811B2 (en) * 2011-07-25 2012-04-24 Clean Power Research, L.L.C. Computer-implemented system and method for determining point-to-point correlation of sky clearness for photovoltaic power generation fleet output estimation
CN104166788B (zh) * 2014-07-22 2017-02-22 国家电网公司 一种架空输电线路最佳经济寿命区间评估方法
CN105631578A (zh) * 2015-12-10 2016-06-01 浙江大学 一种面向风险评估的输变电设备故障概率模型的建模方法
CN105826954B (zh) * 2016-04-28 2018-12-04 四川瑞霆电力科技有限公司 架空线路故障指示器的感应取电系统及方法
CN106056468B (zh) * 2016-05-27 2019-11-26 国网陕西省电力公司经济技术研究院 一种基于可靠性的中压线路典型接线模式选择方法
CN105868872B (zh) * 2016-05-30 2019-08-06 东北大学 一种配电网雷电灾害故障预测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462718A (zh) * 2014-12-25 2015-03-25 国家电网公司 一种变电站经济运行年限区间评估方法
CN106096827A (zh) * 2016-06-02 2016-11-09 国网山东省电力公司济南供电公司 基于混合抽样与最小切负荷计算的可靠性计算方法及平台

Also Published As

Publication number Publication date
CN107480337A (zh) 2017-12-15

Similar Documents

Publication Publication Date Title
CN107315138B (zh) 功率mosfet的故障预测和健康处理方法及测试系统
WO2020215721A1 (zh) 一种激光雷达的寿命预测方法
CN112098876B (zh) 蓄电池内单体电池的异常检测方法
CN104794206A (zh) 一种变电站数据质量评价系统及方法
CN108092622B (zh) 一种基于电阻计算的光伏组串故障诊断方法
CN110502777A (zh) 基于神经网络预测的igbt模块状态检测系统及方法
CN103530511A (zh) 紊流激励条件下的风洞颤振试验颤振边界预测方法
Liu et al. Health assessment and prediction of overhead line based on health index
CN105956752A (zh) 用于绿色公共建筑运营阶段的能效评估与诊断分析方法
CN113759225B (zh) Igbt剩余寿命预测和状态评估实现方法
CN111669123B (zh) 用于对光伏组串进行故障诊断的方法和装置
CN106776480B (zh) 一种无线电干扰现场测量异常值的剔除方法
CN116111727A (zh) 基于动态温度阈值的综合配电箱异常监测方法
CN105719094A (zh) 一种输电设备状态评估方法
CN115358726A (zh) 一种基于物联网的工程施工进度监控系统及方法
CN107480337B (zh) 多因素驱动架空线路故障率建模方法
CN108120932A (zh) 对充电电池的电池健康状态进行估算的方法和装置
CN116090821A (zh) 一种考虑极端天气的配电网线路安全风险评估方法
CN107103425B (zh) 发电设备运行状态计算机智能量化评价系统
CN117368644A (zh) 一种传感器电缆检测方法
CN105488572A (zh) 一种配电设备健康状态评价方法
CN109285331B (zh) 一种基于数据分析和温度预测的电力电缆温度预警系统
CN117080624B (zh) 一种电池加热控制方法及系统
CN106290129A (zh) 一种计算在用钢芯铝绞线的腐蚀寿命的方法及系统
CN112507290B (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