CN103413048A - 基于三参数威布尔分布确定电网设备最佳退役时间的方法 - Google Patents
基于三参数威布尔分布确定电网设备最佳退役时间的方法 Download PDFInfo
- Publication number
- CN103413048A CN103413048A CN2013103561501A CN201310356150A CN103413048A CN 103413048 A CN103413048 A CN 103413048A CN 2013103561501 A CN2013103561501 A CN 2013103561501A CN 201310356150 A CN201310356150 A CN 201310356150A CN 103413048 A CN103413048 A CN 103413048A
- Authority
- CN
- China
- Prior art keywords
- equipment
- parameter
- retired
- weibull distribution
- failure
- 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.)
- Pending
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于三参数威布尔分布确定电网设备最佳退役时间的方法,包括:基于威布尔分布规律建立设备失效率模型;对失效率模型进行参数估计以得到失效率函数;根据失效率函数确定设备最佳退役时间。与现有技术相比,由于本发明的方法具备以下优点:1.具有普适性;2.相对于二参数威布尔分布而言,三参数威布尔分布对于各种类型的实验数据拟合能力更强,尤其是对那些以损耗失效为特征的电网设备拟合精度更高,更能反映设备失效的实际情况;3.能够提供一个相对准确的设备最佳退役时间,从而为制定资产全寿命周期管理策略(尤其是运维策略和退役策略)提供了科学、定量的界定标准,充分发挥了设备自身价值,减少了能源浪费和经济损失。
Description
技术领域
本发明涉及电力企业设备全寿命周期管理技术领域,更具体地涉及一种基于三参数威布尔分布确定电网设备最佳退役时间的方法。
背景技术
对于制定资产全寿命周期管理策略而言,电网设备最佳退役时间是一个非常重要的参考值。然而当前在实际应用中,对电网设备最佳退役时间的界定一直停留在定性分析的层面,多以设备的设计寿命、使用寿命参考值、在运年限等现场经验作为判断依据,缺乏科学、定量的标准和界限,因此容易导致设备提前退役,其价值没有得到充分发挥,造成能源浪费和经济损失。
因此,有必要提供一种改进的确定电网设备最佳退役时间的方法来克服上述缺陷。
发明内容
本发明的目的是提供一种基于三参数威布尔分布确定电网设备最佳退役时间的方法,为制定资产全寿命周期管理策略(尤其是运维策略和退役策略)提供科学、定量的界定标准,充分发挥设备自身价值,减少能源浪费和经济损失。
为实现上述目的,本发明提供了一种基于三参数威布尔分布确定电网设备最佳退役时间的方法,包括:
基于威布尔分布规律建立设备失效率模型;
对所述失效率模型进行参数估计以得到失效率函数,所述参数包括形状参数、尺度参数和位置参数;
根据所述失效率函数确定设备最佳退役时间。
与现有技术相比,由于本发明的方法先基于威布尔分布规律建立设备失效模型,再对失效率模型进行参数估计以得到失效率函数,最后根据失效率函数确定设备最佳退役时间,使得该方法具备以下优点:
1.历史经验表明,电网设备的寿命分布类型中指数分布、威布尔分布、正态分布比例占80%以上,而指数分布和正态分布是威布尔分布的特例,所以本发明以威布尔分布规律建立设备失效率模型,最终确定设备最佳退役时间具有普适性;
2.相对于二参数威布尔分布而言,三参数威布尔分布对于各种类型的实验数据拟合能力更强,尤其是对那些以损耗失效为特征的电网设备拟合精度更高,更能反映设备失效的实际情况;
3.能够提供一个相对准确的设备最佳退役时间,从而为制定资产全寿命周期管理策略(尤其是运维策略和退役策略)提供了科学、定量的界定标准,充分发挥了设备自身价值,减少了能源浪费和经济损失。
通过以下的描述并结合附图,本发明将变得更加清晰,这些附图用于解释本发明的实施例。
附图说明
图1为本发明方法一实施例的流程图。
图2为本发明运维成本曲线和退役处置成本折现曲线图。
具体实施方式
现在参考附图描述本发明的实施例,附图中类似的元件标号代表类似的元件。
请参考图1,以SF6断路器为例,详细描述本发明方法的具体步骤。如图1所示,该方法包括:
S101,基于威布尔分布规律建立设备失效率模型;具体为,假设设备在全寿命周期内的累计失效概率F(t)服从具有形状参数m、尺度参数η、位置参数γ三个可变参数的威布尔分布规律,表示为由此得到断路器失效率模型为
需要说明的是,电网设备在全寿命周期内的失效率(故障率)是时间的函数,因其分布规律呈浴盆形状而称为浴盆曲线,分为早期失效期、偶发失效期和耗损失效期三个阶段。威布尔分布能够很好地拟合浴盆曲线,上述三个阶段分别对应形状参数m小于1(大于零)、等于1、大于1三种情况;
S102,对设备样本组进行预处理;具体为,设备样本预处理是对设备寿命试验数据进行筛选,选取设计/制造相同、运行环境相似的10台SF6断路器作为一个样本组,寿命统计区间以5年为一段;本实施例样本组的寿命试验数据为15.6,16.2,16.8,17.3,17.9,18.1,18.7,19.5,19.7,20.1;
S103,根据经预处理后的设备样本组对失效率模型进行参数估计以得到失效率函数;具体为,参数估计是利用经过预处理的一组样本数据对三参数威布尔分布中的形状参数、尺度参数、位置参数进行估算,从而得到失效率函数。在估算时,样本数据中第i个设备对应的累计失效概率可利用中位秩法近似求得,其中N为样本数量;本实施例中N=10,对应的累计失效概率为0.07,0.16,0.26,0.36,0.45,0.55,0.64,0.74,0.84,0.93;形状参数m和尺度参数η均由最小二乘法进行回归求出;位置参数γ用最大相关系数优化法求出;具体地,最小二乘法回归计算是对累计失效概率F(t)两侧同取两次自然对数,将得到的方程转化为线性方程y=mx-B,其中x=ln(t-γ),B=lnηm;本实施例中x,y的值如下表所示:
假设已知位置参数γ,则可以利用最小二乘法进行回归求出形状参数m、尺度参数η;本实施例中,形状参数和尺度参数分别为2.81,4.47;
具体地,最大相关系数优化法是将线性方程y=mx-B中自变量x和因变量y之间的相关系数对位置参数γ一阶求导,令导数为零,得到的γ值即为所求;本实施例中,γ=13.65。因此与所选样本组设计/制造相同、运行环境相似的SF6断路器失效率函数为
S104,根据失效率函数确定设备最佳退役时间;具体为,假设设备运维成本与失效率函数成正比,比例系数为设备初期投资费用。本实施例中SF6断路器的投资成本为270000元,与上述失效率的乘积得到运维成本曲线如图2实线所示,随时间增加呈现非线性递增的规律;假设该断路器退役处置费用(包括进行退役处置时的人工、设备费用以及运输费)为110000元,设备残值(设备在计算周期末的剩余价值)为2800元,则退役处置成本终值为110000–2800=107200元;若假设利率为8%,则退役处置成本折现曲线如图2虚线所示,随时间增加呈现非线性递减的规律。两条曲线的交点——约16.9年——即为与符合该样本组累计失效概率分布的SF6断路器的最佳退役时间。
从以上描述可以看出,本发明的方法具有以下优点:
1.历史经验表明,电网设备的寿命分布类型中指数分布、威布尔分布、正态分布比例占80%以上,而指数分布和正态分布是威布尔分布的特例,所以本发明以威布尔分布规律建立设备失效率模型,最终确定设备最佳退役时间具有普适性;
2.相对于二参数威布尔分布而言,三参数威布尔分布对于各种类型的实验数据拟合能力更强,尤其是对那些以损耗失效为特征的电网设备拟合精度更高,更能反映设备失效的实际情况;
3.利用最小二乘法和最大相关系数优化法确定威布尔分布系数,计算量不大,且能保证一定的精度;
4.能够提供一个相对准确的设备最佳退役时间,从而为制定资产全寿命周期管理策略(尤其是运维策略和退役策略)提供了科学、定量的界定标准,充分发挥了设备自身价值,减少了能源浪费和经济损失。
以上结合最佳实施例对本发明进行了描述,但本发明并不局限于以上揭示的实施例,而应当涵盖各种根据本发明的本质进行的修改、等效组合。
Claims (7)
1.一种基于三参数威布尔分布确定电网设备最佳退役时间的方法,其特征在于,包括:
基于威布尔分布规律建立设备失效率模型;
对所述失效率模型进行参数估计以得到失效率函数;
根据所述失效率函数确定设备最佳退役时间。
2.如权利要求1所述的方法,其特征在于,“基于威布尔分布规律建立设备失效率模型”之后还包括:
对设备样本组进行预处理。
3.如权利要求2所述的方法,其特征在于,“对设备样本进行预处理”具体包括:
选取设计/制造相同、运行环境相似的多个设备作为一个所述设备样本组;
统计预设寿命统计区间内所述设备样本组内的多个寿命试验数据。
5.如权利要求4所述的方法,其特征在于,所述“采用最大相关系数法计算出所述位置参数”具体包括:
将线性方程y=mx-B中自变量x和因变量y之间的相关系数对位所述置参数γ一阶求导,令导数为零以得到所述置参数γ。
6.如权利要求4所述的方法,其特征在于,“根据所述失效率函数确定设备最佳退役时间”具体包括:
根据设备初期投资费用和所述失效率函数绘制运维成本曲线;
根据设备退役处置成本和设备残值绘制退役处置成本折现曲线;
根据所述运维成本曲线和退役处置成本折现曲线确定所述设备最佳退役时间。
7.如权利要求1所述的方法,其特征在于,所述参数包括形状参数、尺度参数和位置参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013103561501A CN103413048A (zh) | 2013-08-15 | 2013-08-15 | 基于三参数威布尔分布确定电网设备最佳退役时间的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013103561501A CN103413048A (zh) | 2013-08-15 | 2013-08-15 | 基于三参数威布尔分布确定电网设备最佳退役时间的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103413048A true CN103413048A (zh) | 2013-11-27 |
Family
ID=49606059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013103561501A Pending CN103413048A (zh) | 2013-08-15 | 2013-08-15 | 基于三参数威布尔分布确定电网设备最佳退役时间的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103413048A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483192A (zh) * | 2014-11-27 | 2015-04-01 | 福达合金材料股份有限公司 | 一种基于威布尔分布的电触头材料静熔焊力数据处理方法 |
CN105891645A (zh) * | 2016-05-31 | 2016-08-24 | 西安航空制动科技有限公司 | 确定防滑刹车控制装置振动故障分布的方法 |
CN105930632A (zh) * | 2016-03-02 | 2016-09-07 | 航天科工防御技术研究试验中心 | 机电整机产品贮存寿命建模方法 |
CN106054105A (zh) * | 2016-05-20 | 2016-10-26 | 国网新疆电力公司电力科学研究院 | 一种智能电表的可靠性预计修正模型建立方法 |
CN106647273A (zh) * | 2016-12-26 | 2017-05-10 | 北京天源科创风电技术有限责任公司 | 预测部件的预防性更换时间的方法和设备 |
CN108108542A (zh) * | 2017-12-14 | 2018-06-01 | 河北工业大学 | 低压成套开关设备的寿命预测方法 |
CN109101466A (zh) * | 2018-11-22 | 2018-12-28 | 中国人民解放军国防科技大学 | 基于分布函数取对数变换的威布尔分布参数估计方法 |
CN113168597A (zh) * | 2018-11-08 | 2021-07-23 | 施乐百有限公司 | 用于预测风扇组的故障的方法和***以及对应的风扇组 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101038638A (zh) * | 2007-04-25 | 2007-09-19 | 上海发电设备成套设计研究院 | 一种发电机组自动控制***电子元器件剩余寿命预测方法 |
US20110144924A1 (en) * | 2006-12-22 | 2011-06-16 | The Boeing Company | Computational method for load enhancement factors |
CN102945316A (zh) * | 2012-10-25 | 2013-02-27 | 华北电力大学 | 一种考虑协变量的继电保护设备失效率计算方法 |
-
2013
- 2013-08-15 CN CN2013103561501A patent/CN103413048A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110144924A1 (en) * | 2006-12-22 | 2011-06-16 | The Boeing Company | Computational method for load enhancement factors |
CN101038638A (zh) * | 2007-04-25 | 2007-09-19 | 上海发电设备成套设计研究院 | 一种发电机组自动控制***电子元器件剩余寿命预测方法 |
CN102945316A (zh) * | 2012-10-25 | 2013-02-27 | 华北电力大学 | 一种考虑协变量的继电保护设备失效率计算方法 |
Non-Patent Citations (1)
Title |
---|
周瑞祥等: "随机截尾机电设备三参数威布尔分布模型及其应用研究", 《现代制造工程》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483192B (zh) * | 2014-11-27 | 2017-01-25 | 福达合金材料股份有限公司 | 一种基于威布尔分布的电触头材料静熔焊力数据处理方法 |
CN104483192A (zh) * | 2014-11-27 | 2015-04-01 | 福达合金材料股份有限公司 | 一种基于威布尔分布的电触头材料静熔焊力数据处理方法 |
CN105930632A (zh) * | 2016-03-02 | 2016-09-07 | 航天科工防御技术研究试验中心 | 机电整机产品贮存寿命建模方法 |
CN105930632B (zh) * | 2016-03-02 | 2018-10-02 | 航天科工防御技术研究试验中心 | 机电整机产品贮存寿命建模方法 |
CN106054105B (zh) * | 2016-05-20 | 2019-01-15 | 国网新疆电力公司电力科学研究院 | 一种智能电表的可靠性预计修正模型建立方法 |
CN106054105A (zh) * | 2016-05-20 | 2016-10-26 | 国网新疆电力公司电力科学研究院 | 一种智能电表的可靠性预计修正模型建立方法 |
CN105891645A (zh) * | 2016-05-31 | 2016-08-24 | 西安航空制动科技有限公司 | 确定防滑刹车控制装置振动故障分布的方法 |
CN105891645B (zh) * | 2016-05-31 | 2018-10-09 | 西安航空制动科技有限公司 | 确定防滑刹车控制装置振动故障分布的方法 |
CN106647273A (zh) * | 2016-12-26 | 2017-05-10 | 北京天源科创风电技术有限责任公司 | 预测部件的预防性更换时间的方法和设备 |
CN108108542A (zh) * | 2017-12-14 | 2018-06-01 | 河北工业大学 | 低压成套开关设备的寿命预测方法 |
CN113168597A (zh) * | 2018-11-08 | 2021-07-23 | 施乐百有限公司 | 用于预测风扇组的故障的方法和***以及对应的风扇组 |
CN109101466A (zh) * | 2018-11-22 | 2018-12-28 | 中国人民解放军国防科技大学 | 基于分布函数取对数变换的威布尔分布参数估计方法 |
CN109101466B (zh) * | 2018-11-22 | 2019-03-22 | 中国人民解放军国防科技大学 | 基于分布函数取对数变换的威布尔分布参数估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103413048A (zh) | 基于三参数威布尔分布确定电网设备最佳退役时间的方法 | |
Zhou et al. | Framework for capacity credit assessment of electrical energy storage and demand response | |
Warren | A review of demand-side management policy in the UK | |
CN105069535A (zh) | 一种基于arima模型的配电网运行可靠性的预测方法 | |
CN102509173B (zh) | 一种基于马尔科夫链的电力***负荷准确预测方法 | |
CN103473438B (zh) | 风功率预测模型优选及修正方法 | |
CN105225010A (zh) | 一种基于可靠性的变压器设备寿命评估方法 | |
CN103383445A (zh) | 一种智能电表寿命和可靠性预测***及其方法 | |
CN103457279A (zh) | 大规模电网的分区d-statcom综合优化配置方案的确定方法 | |
CN105740975A (zh) | 一种基于数据关联关系的设备缺陷评估与预测方法 | |
CN102306346A (zh) | 基于可靠性规划的中压配电网目标网架结构优化决策方法 | |
Shayesteh et al. | Scenario reduction, network aggregation, and DC linearisation: which simplifications matter most in operations and planning optimisation? | |
CN103996147A (zh) | 配电网综合评估方法 | |
CN105610192A (zh) | 考虑大规模风电接入的在线风险评估方法 | |
Wang et al. | Optimal sizing of energy storage considering the spatial‐temporal correlation of wind power forecast errors | |
Qiu et al. | Probabilistic transmission expansion planning for increasing wind power penetration | |
CN102751725B (zh) | 一种配电网过负荷风险状态辨识方法 | |
CN105825295A (zh) | 一种计及元胞发展程度的空间负荷预测方法 | |
Zhou et al. | Dynamic battery loss evaluation and its application for optimal online wind‐storage integrated scheduling | |
CN107332239B (zh) | 一种基于配电网等值的输配电网协调规划方法 | |
Guo et al. | The Forecasting Procedure for Long‐Term Wind Speed in the Zhangye Area | |
CN103700034A (zh) | 一种基于边际成本理论的配电网网架规划方法 | |
Wang et al. | A trend‐based approach for situation awareness in power systems | |
Uppal et al. | Weather biased optimal delta model for short‐term load forecast | |
Yang | Electricity Price Forecast Based on Metabolic GM (1, 1) Model under Wind Power Development Prospect |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20131127 |
|
RJ01 | Rejection of invention patent application after publication |