CN107122878A - 一种基于大数据与人为因素的电力系统风险评估方法 - Google Patents
一种基于大数据与人为因素的电力系统风险评估方法 Download PDFInfo
- Publication number
- CN107122878A CN107122878A CN201710124484.4A CN201710124484A CN107122878A CN 107122878 A CN107122878 A CN 107122878A CN 201710124484 A CN201710124484 A CN 201710124484A CN 107122878 A CN107122878 A CN 107122878A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- probability
- equipment
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004458 analytical method Methods 0.000 title claims abstract description 25
- 238000011156 evaluation Methods 0.000 title claims abstract description 15
- 230000005540 biological transmission Effects 0.000 claims abstract description 15
- 230000001131 transforming effect Effects 0.000 claims abstract description 15
- 238000005516 engineering process Methods 0.000 claims abstract description 13
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000003542 behavioural effect Effects 0.000 claims description 16
- 238000004891 communication Methods 0.000 claims description 12
- 238000012423 maintenance Methods 0.000 claims description 12
- 238000004445 quantitative analysis Methods 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000008439 repair process Effects 0.000 claims description 3
- 230000000717 retained effect Effects 0.000 claims description 3
- 230000004888 barrier function Effects 0.000 claims 3
- 238000011160 research Methods 0.000 abstract description 4
- 238000012502 risk assessment Methods 0.000 abstract description 4
- 230000007257 malfunction Effects 0.000 abstract description 3
- 230000009183 running Effects 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于大数据与人为因素的电力系统风险评估方法。本方法分为四个步骤,第一,利用AHP‑SLIM方法评估系统运行过程中操作人员发生人为失误的概率;第二,采用大数据挖掘技术得出输变电设备、二次系统设备及软件的失效概率;第三,基于IEC‑61850的功能分解对二次设备或软件的各种功能的失效风险进行评估;最后,两者相结合对传统一次设备的故障概率进行修正,完成系统风险评估。本发明在现有研究的基础上,将操作人员与二次系统的失效概率考虑进一次系统的风险评估中,使得电力系统的风险评估的结果更为准确,相较于传统评估方法能更为客观地表征当前电网所面临的风险。
Description
技术领域
本发明属于电力系统风险评估领域,具体涉及了一种基于大数据与 人为因素的电力系统风险评估方法。
背景技术
目前,针对电力系统风险评估的研究较多,但主要集中在对一次设 备的故障概率与故障后果的研究上,而少有将人因可靠性与二次设备停 运引发的系统故障考虑进系统的失效机理中,因此,评估得到的风险结 果往往不准确,有时甚至会存在较大出入。
随着电力信息技术的不断发展,电力系统的数据呈爆炸式增长,具 有海量、异构的特点,采用常规软件、工具分析费时、耗力,而大数据 挖掘技术具有快速处理此类数据的优势,因此本发明将采用大数据挖掘 技术对输变电设备、二次系统设备、软件的可用率与其历史数据、当前 数据、人为因素等进行关联分析,得到当前的设备可用率。
本发明提出了一种基于大数据与人为因素的电力系统风险评估方法。 首先利用AHP-SLIM方法评估系统运行过程中操作人员发生人为失误的 概率;其次采用大数据挖掘技术得出输变电设备、二次系统设备及软件 的失效概率;再次基于IEC-61850的功能分解对二次设备或软件的各种 功能的失效风险进行评估,两者相结合对传统一次设备的故障概率进行 修正,实现电力系统风险评估的精确化,可弥补目前电力系统风险评估 存在的局限。
发明内容
针对现有电力系统风险评估存在的上述问题,本发明的目的在于提 出一种基于大数据与人为因素的电力系统风险评估方法。本方法分为四 个步骤,第一,利用AHP-SLIM方法评估系统运行过程中操作人员发生 人为失误的概率;第二,采用大数据挖掘技术得出输变电设备、二次系 统设备及软件的失效概率;第三,基于IEC-61850的功能分解对二次设备或软件的各种功能的失效风险进行评估;最后,两者相结合对传统一 次设备的故障概率进行修正,完成系统风险评估。
本方法为解决其技术问题所采用的技术方案包括如下步骤:
步骤1、构建行为影响体系,选择情景下的影响因素。
针对电力系统工作人员的操作环境,结合常规人因可靠性分析理论, 建立了相应的行为影响体系。该体系共分为5类一级影响因素,包括: 个人因素、组织因素、团队因素、环境因素以及信息因素。每种一级影 响因素下又包含多种二级因素,详细情况如表1所示。每种影响因素代 表了电力系统工作人员在工作时其操作可能受到影响的方面,该体系不仅适用于人为失误根源的定性分析,还可作为各因素对人员行为影响的 定量分析的基础。为了定量分析每种因素的影响程度,一般采用面向专 家的问卷调查方法。每类二级因素均具有1至10的等级划分,数字越小 说明情况越恶劣。
表1行为影响体系
步骤2、基于步骤1中的行为影响体系,利用AHP-SLIM方法评估 人为失效概率。
步骤2-1、利用层次分析法(AHP)计算各项行为影响因素的权重系 数。
步骤2-1-1、假设分析场景下有n个二级行为影响因素,则征询本领 域专家意见,依据表2对这些影响因素进行重要度评分(两两打分,评 分范围为1-10的整数)并进行比较,构建判别矩阵A如下:
该矩阵满足aji=1/aij,且aji(j<=i)的取值遵循表2。
步骤2-1-2、对矩阵A中的各列进行归一化处理,得到A1矩阵。
步骤2-1-3、A1矩阵各行求和,得到列向量W。
步骤2-1-4、对列向量W进行归一化处理,即得权重系数向量W1。
表2AHP判别矩阵取值说明
步骤2-2、利用式(2)计算成功似然指数SLI。
SLI=Σωivi,0≤SLI≤10 (2)
式中,ωi为第i项影响因素的重要度权重,vi为第i项影响因素的评 定等级,a和b均为常数。
步骤2-3、利用式(3)计算人为失误概率,
HEP=exp(a×SIL+b) (3)
式中,a与b均为常数,依据历史数据进行拟合得到。
步骤3、采用大数据挖掘技术,结合HEP及其他各项因素,得出输 变电设备、二次系统设备及软件的失效概率。
步骤3-1、获取当前数据、历史数据,包括负载率K、天气W、检 修R、历史故障数据F、设备故障率λ、设备修复率μ,以及步骤2所 得的人为失误概率HEP。对各项数据进行离散化处理,在历史数据的各 个时间点上,负载率K分为轻载1,重载2以及过载3;天气W依据是 否有气象局的预警分为恶劣天气1与正常天气2;检修R分为最近6个 月内有检修1和无检修2;历史故障数据F分为最近6个月内发生故障1 和未发生故障2;设备故障率λ分为低故障率1,一般故障率2,高故障 率3;人为失误概率HEP分为低概率1、一般概率2与高概率3。负载率K、设备故障率λ以及人为失误概率HEP的离散化对应关系如表3所示, 其中λ6为近6个月来的设备平均故障率。
表3离散化对应关系
步骤3-2、枚举关联规则
1)对于输变电设备,考虑K、W、R与F对λ的影响,共321条规 则;
2)对于二次系统设备、软件,考虑W、R、F与HEP对λ的影响, 共321条规则。
规则形如A→B,A为条件(可包含任意项),B为结果,例如
K1,W1,R2、F2→λ3
K2,W2,R2→λ2
F1,HEP2→λ1
步骤3-3、计算各条关联规则的支持度S和置信度C,N(A∪B)表示 A和B同时出现的次数,N表示总数据记录数,N(A)表示A出现的次数。
舍弃支持度小于20%或置信度小于50%的规则。
步骤3-4、对于被保留的所有规则,利用极大似然方法拟合设备故障 率关于规则条件中影响因素的关系函数。
步骤3-5、根据步骤3-4中得出的关系函数,代入当前各台设备或软 件的负载率、天气、检修、历史故障情况以及人为失误概率,得到当前 设备故障率。若没有符合的规则,则故障率取为历史平均故障率。计算 输变电设备的失效概率P1,计算二次系统设备或软件的失效概率P2。
步骤4、对系统二次设备或软件进行功能分解。
基于IEC-61850的功能分解,将系统各类功能分解为逻辑节点、逻 辑连接,通信信息片经逻辑连接在逻辑节点间传递。
步骤5、计算各功能失效概率。
5-1.在计算各功能失效概率前,作以下规定:
1)某二次设备或软件中的所有逻辑节点与逻辑连接具有相同的失效 概率,与该设备失效概率相一致;
2)功能、逻辑节点及逻辑连接仅存在两种状态,工作状态与失效状 态;
3)当某逻辑节点或逻辑连接处于失效状态时,其所属功能也处于失 效状态;
4)不同功能间的失效事件相互独立;
5)忽略所有通信时间延迟。
5-2.基于步骤5-1的相关规定,可得到:
1)设备D(失效概率为PD)中的逻辑节点LN,其失效概率为 pLN=pD。
2)连接设备D中两个逻辑节点的逻辑连接LC,其失效概率为 pLC=pD
3)对于连接不同设备(D1,D2)中的两个逻辑节点(LN1,LN2) 的逻辑连接LC,其失效概率为:
式中,m为两台设备间的通信连接开关数,pswi是开关i的失效概率, n为两台设备间的通信线路数,pclj为通信线路j的失效概率。
依据功能分解原理,可将功能视作一系列逻辑节点与逻辑连接的集 合,因此功能F的失效概率为:
式中,NLN和NLC分别表示本功能中所包含的逻辑节点数与逻辑连接 数,PLNi和PLCj则分别表示第i个逻辑节点与第j个逻辑连接的失效概率。
步骤6、对电力系统风险进行评估。
步骤6-1、采用非序贯蒙特卡洛采样技术,获取每个场景下一次系统 设备与二次系统功能的状态,统计得到一次系统设备与二次系统功能的 失效概率。
步骤6-2、利用二次系统功能的失效概率修正一次设备的失效概率:
式中,pE和pE’分别表示不考虑与考虑二次系统功能失效概率的一次 系统设备失效概率,SE表示影响一次设备的功能序号集合,pFi表示第i 个功能的失效概率。
步骤6-3、在修正后的一次系统设备失效概率基础上,进一步采用非 序贯蒙特卡洛采样技术,对电网场景采样,并计算每个场景下的失负荷 量。
步骤6-4、利用式(10)计算得到最终的失负荷期望值。
pi为系统在状态i的概率;S为系统所有切负荷状态的集合,Ci为系 统在状态i的失负荷量。
本发明的有益效果是:
本发明在现有研究的基础上,将操作人员与二次系统的失效概率考 虑进一次系统的风险评估中,使得电力系统的风险评估的结果更为准确, 相较于传统评估方法能更为客观地表征当前电网所面临的风险。
附图说明
图1是本发明的基于大数据与人为因素的电力系统风险评估方法的 总体流程图。
图2是本发明的基于大数据挖掘技术的输变电设备、二次系统设备 及软件失效概率计算的流程图。
图3是本发明的电力系统风险评估的流程图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
如图1所示,本方法为解决其技术问题所采用的技术方案包括如下 步骤:
步骤1、构建行为影响体系,选择情景下的影响因素。
针对电力系统工作人员的操作环境,结合常规人因可靠性分析理论, 建立了相应的行为影响体系。该体系共分为5类一级影响因素,包括: 个人因素、组织因素、团队因素、环境因素以及信息因素。每种一级影 响因素下又包含多种二级因素,详细情况如表1所示。每种影响因素代 表了电力系统工作人员在工作时其操作可能受到影响的方面,该体系不仅适用于人为失误根源的定性分析,还可作为各因素对人员行为影响的 定量分析的基础。为了定量分析每种因素的影响程度,一般采用面向专 家的问卷调查方法。每类二级因素均具有1至10的等级划分,数字越小 说明情况越恶劣。
表1行为影响体系
步骤2、基于步骤1中的行为影响体系,利用AHP-SLIM方法评估 人为失效概率。
步骤2-1、利用层次分析法(AHP)计算各项行为影响因素的权重系 数。
步骤2-1-1、假设分析场景下有n个二级行为影响因素,则征询本领 域专家意见,依据表2对这些影响因素进行重要度评分(两两打分,评 分范围为1-10的整数)并进行比较,构建判别矩阵A如下:
该矩阵满足aji=1/aij,且aji(j<=i)的取值遵循表2。
步骤2-1-2、对矩阵A中的各列进行归一化处理,得到A1矩阵。
步骤2-1-3、A1矩阵各行求和,得到列向量W。
步骤2-1-4、对列向量W进行归一化处理,即得权重系数向量W1。
表2AHP判别矩阵取值说明
步骤2-2、利用式(2)计算成功似然指数SLI。
SLI=∑ωivi,0≤SLI≤10 (2)
式中,ωi为第i项影响因素的重要度权重,vi为第i项影响因素的评 定等级,a和b均为常数。
步骤2-3、利用式(3)计算人为失误概率,
HEP=exp(a×SIL+b) (3)
式中,a与b均为常数,依据历史数据进行拟合得到。
步骤3、采用大数据挖掘技术,结合HEP及其他各项因素,得出输 变电设备、二次系统设备及软件的失效概率。
如图2所示,本发明的基于大数据挖掘技术的输变电设备、二次系 统设备及软件失效概率计算,包括如下步骤:
步骤3-1、获取当前数据、历史数据,包括负载率K、天气W、检 修R、历史故障数据F、设备故障率λ、设备修复率μ,以及步骤2所 得的人为失误概率HEP。对各项数据进行离散化处理,在历史数据的各 个时间点上,负载率K分为轻载1,重载2以及过载3;天气W依据是 否有气象局的预警分为恶劣天气1与正常天气2;检修R分为最近6个 月内有检修1和无检修2;历史故障数据F分为最近6个月内发生故障1 和未发生故障2;设备故障率λ分为低故障率1,一般故障率2,高故障 率3;人为失误概率HEP分为低概率1、一般概率2与高概率3。负载率K、设备故障率λ以及人为失误概率HEP的离散化对应关系如表3所示, 其中λ6为近6个月来的设备平均故障率。
表3离散化对应关系
步骤3-2、枚举关联规则
1)对于输变电设备,考虑K、W、R与F对λ的影响,共321条规 则;
2)对于二次系统设备、软件,考虑W、R、F与HEP对λ的影响, 共321条规则。
规则形如A→B,A为条件(可包含任意项),B为结果,例如
K1,W1,R2、F2→λ3
K2,W2,R2→λ2
F1,HEP2→λ1
步骤3-3、计算各条关联规则的支持度S和置信度C,N(A∪B)表示 A和B同时出现的次数,N表示总数据记录数,N(A)表示A出现的次数。
舍弃支持度小于20%或置信度小于50%的规则。
步骤3-4、对于被保留的所有规则,利用极大似然方法拟合设备故障 率关于规则条件中影响因素的关系函数。
步骤3-5、根据步骤3-4中得出的关系函数,代入当前各台设备或软 件的负载率、天气、检修、历史故障情况以及人为失误概率,得到当前 设备故障率。若没有符合的规则,则故障率取为历史平均故障率。计算 输变电设备的失效概率P1,计算二次系统设备或软件的失效概率P2。
步骤4、对系统二次设备或软件进行功能分解。
基于IEC-61850的功能分解,将系统各类功能分解为逻辑节点、逻辑 连接,通信信息片经逻辑连接在逻辑节点间传递。
步骤5、计算各功能失效概率。
5-1.作以下规定:
1)某二次设备或软件中的所有逻辑节点与逻辑连接具有相同的失效 概率,与该设备失效概率相一致;
2)功能、逻辑节点及逻辑连接仅存在两种状态,工作状态与失效状 态;
3)当某逻辑节点或逻辑连接处于失效状态时,其所属功能也处于失 效状态;
4)不同功能间的失效事件相互独立;
5)忽略所有通信时间延迟。
5-2.基于5-1的规定,可得到:
1)设备D(失效概率为PD)中的逻辑节点LN,其失效概率为 pLN=pD。
2)连接设备D中两个逻辑节点的逻辑连接LC,其失效概率为 pLC=pD
3)对于连接不同设备(D1,D2)中的两个逻辑节点(LN1,LN2) 的逻辑连接LC,其失效概率为:
式中,m为两台设备间的通信连接开关数,pswi是开关i的失效概率, n为两台设备间的通信线路数,pclj为通信线路j的失效概率。
依据功能分解原理,可将功能视作一系列逻辑节点与逻辑连接的集 合,因此功能F的失效概率为:
式中,NLN和NLC分别表示本功能中所包含的逻辑节点数与逻辑连接 数,PLNi和PLCj则分别表示第i个逻辑节点与第j个逻辑连接的失效概率。
步骤6、对电力系统风险进行评估。
如图3所示,本发明的电力系统风险评估,包括如下步骤:
步骤6-1、采用非序贯蒙特卡洛采样技术,获取每个场景下一次系统 设备与二次系统功能的状态,统计得到一次系统设备与二次系统功能的 失效概率。
步骤6-2、利用二次系统功能的失效概率修正一次设备的失效概率:
式中,pE和pE’分别表示不考虑与考虑二次系统功能失效概率的一次 系统设备失效概率,SE表示影响一次设备的功能序号集合,pFi表示第i 个功能的失效概率。
步骤6-3、在修正后的一次系统设备失效概率基础上,进一步采用非 序贯蒙特卡洛采样技术,对电网场景采样,并计算每个场景下的失负荷 量。
步骤6-4、利用式(10)计算得到最终的失负荷期望值。
pi为系统在状态i的概率;S为系统所有切负荷状态的集合,Ci为系 统在状态i的失负荷量。
Claims (1)
1.一种基于大数据与人为因素的电力系统风险评估方法,其特征在于包括如下步骤:
步骤1、构建行为影响体系,选择情景下的影响因素;
针对电力系统工作人员的操作环境,结合常规人因可靠性分析理论,建立了相应的行为影响体系;该体系共分为5类一级影响因素,包括:个人因素、组织因素、团队因素、环境因素以及信息因素;每种一级影响因素下又包含多种二级因素;每种影响因素代表了电力系统工作人员在工作时其操作可能受到影响的方面,该体系不仅适用于人为失误根源的定性分析,还可作为各因素对人员行为影响的定量分析的基础;为了定量分析每种因素的影响程度,一般采用面向专家的问卷调查方法;每类二级因素均具有1至10的等级划分,数字越小说明情况越恶劣;
步骤2、基于步骤1中的行为影响体系,利用AHP-SLIM方法评估人为失效概率;
步骤2-1、利用层次分析法计算各项行为影响因素的权重系数;
步骤2-1-1、假设分析场景下有n个二级行为影响因素,则征询本领域专家意见,依据表2对这些影响因素进行重要度评分并进行比较,构建判别矩阵A如下:
该矩阵满足aji=1/aij,且aji(j<=i)的取值遵循表2;
步骤2-1-2、对矩阵A中的各列进行归一化处理,得到A1矩阵;
步骤2-1-3、A1矩阵各行求和,得到列向量W;
步骤2-1-4、对列向量W进行归一化处理,即得权重系数向量W1;
表2 AHP判别矩阵取值说明
步骤2-2、利用式(2)计算成功似然指数SLI;
SLI=∑ωivi,0≤SLI≤10 (2)
式中,ωi为第i项影响因素的重要度权重,vi为第i项影响因素的评定等级,a和b均为常数;
步骤2-3、利用式(3)计算人为失误概率,
HEP=exp(a×SIL+b) (3)
式中,a与b均为常数,依据历史数据进行拟合得到;
步骤3、采用大数据挖掘技术,结合HEP及其他各项因素,得出输变电设备、二次系统设备及软件的失效概率;
步骤3-1、获取当前数据、历史数据,包括负载率K、天气W、检修R、历史故障数据F、设备故障率λ、设备修复率μ,以及步骤2所得的人为失误概率HEP;对各项数据进行离散化处理,在历史数据的各个时间点上,负载率K分为轻载1,重载2以及过载3;天气W依据是否有气象局的预警分为恶劣天气1与正常天气2;检修R分为最近6个月内有检修1和无检修2;历史故障数据F分为最近6个月内发生故障1和未发生故障2;设备故障率λ分为低故障率1,一般故障率2,高故障率3;人为失误概率HEP分为低概率1、一般概率2与高概率3;负载率K、设备故障率λ以及人为失误概率HEP的离散化对应关系如表3所示,其中λ6为近6个月来的设备平均故障率;
表3 离散化对应关系
步骤3-2、枚举关联规则
1)对于输变电设备,考虑K、W、R与F对λ的影响,共321条规则;
2)对于二次系统设备、软件,考虑W、R、F与HEP对λ的影响,共321条规则;
规则形如A→B,A为条件(可包含任意项),B为结果,例如
K1,W1,R2、F2→λ3
K2,W2,R2→λ2
F1,HEP2→λ1
步骤3-3、计算各条关联规则的支持度S和置信度C,N(A∪B)表示A和B同时出现的次数,N表示总数据记录数,N(A)表示A出现的次数;
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>&RightArrow;</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>&cup;</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>N</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>&RightArrow;</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>&cup;</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
舍弃支持度小于20%或置信度小于50%的规则;
步骤3-4、对于被保留的所有规则,利用极大似然方法拟合设备故障率关于规则条件中影响因素的关系函数;
步骤3-5、根据步骤3-4中得出的关系函数,代入当前各台设备或软件的负载率、天气、检修、历史故障情况以及人为失误概率,得到当前设备故障率;若没有符合的规则,则故障率取为历史平均故障率;计算输变电设备的失效概率P1,计算二次系统设备或软件的失效概率P2;
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>,</mo>
<mi>W</mi>
<mo>,</mo>
<mi>R</mi>
<mo>,</mo>
<mi>F</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>,</mo>
<mi>W</mi>
<mo>,</mo>
<mi>R</mi>
<mo>,</mo>
<mi>F</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&mu;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>W</mi>
<mo>,</mo>
<mi>R</mi>
<mo>,</mo>
<mi>F</mi>
<mo>,</mo>
<mi>H</mi>
<mi>E</mi>
<mi>P</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>W</mi>
<mo>,</mo>
<mi>R</mi>
<mo>,</mo>
<mi>F</mi>
<mo>,</mo>
<mi>H</mi>
<mi>E</mi>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&mu;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤4、对系统二次设备或软件进行功能分解;
基于IEC-61850的功能分解,将系统各类功能分解为逻辑节点、逻辑连接,通信信息片经逻辑连接在逻辑节点间传递;
步骤5、计算各功能失效概率;
5-1.在计算各功能失效概率前,作以下规定:
1)某二次设备或软件中的所有逻辑节点与逻辑连接具有相同的失效概率,与该设备失效概率相一致;
2)功能、逻辑节点及逻辑连接仅存在两种状态,工作状态与失效状态;
3)当某逻辑节点或逻辑连接处于失效状态时,其所属功能也处于失效状态;
4)不同功能间的失效事件相互独立;
5)忽略所有通信时间延迟;
5-2.基于步骤5-1的相关规定,可得到:
1)设备D(失效概率为PD)中的逻辑节点LN,其失效概率为pLN=pD;
2)连接设备D中两个逻辑节点的逻辑连接LC,其失效概率为pLC=pD
3)对于连接不同设备(D1,D2)中的两个逻辑节点(LN1,LN2)的逻辑连接LC,其失效概率为:
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>L</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<munderover>
<mi>&Pi;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>sw</mi>
<mi>i</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<munderover>
<mi>&Pi;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>cl</mi>
<mi>j</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,m为两台设备间的通信连接开关数,pswi是开关i的失效概率,n为两台设备间的通信线路数,pclj为通信线路j的失效概率;
依据功能分解原理,可将功能视作一系列逻辑节点与逻辑连接的集合,因此功能F的失效概率为:
<mrow>
<msub>
<mi>p</mi>
<mi>F</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>L</mi>
<mi>N</mi>
</mrow>
</msub>
</munderover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>LM</mi>
<mi>i</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>L</mi>
<mi>C</mi>
</mrow>
</msub>
</munderover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>LC</mi>
<mi>j</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,NLN和NLC分别表示本功能中所包含的逻辑节点数与逻辑连接数,PLNi和PLCj则分别表示第i个逻辑节点与第j个逻辑连接的失效概率;
步骤6、对电力系统风险进行评估;
步骤6-1、采用非序贯蒙特卡洛采样技术,获取每个场景下一次系统设备与二次系统功能的状态,统计得到一次系统设备与二次系统功能的失效概率;
步骤6-2、利用二次系统功能的失效概率修正一次设备的失效概率:
<mrow>
<msubsup>
<mi>p</mi>
<mi>E</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>p</mi>
<mi>E</mi>
</msub>
<mo>)</mo>
</mrow>
<munder>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mi>E</mi>
</msub>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>F</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,pE和pE’分别表示不考虑与考虑二次系统功能失效概率的一次系统设备失效概率,SE表示影响一次设备的功能序号集合,pFi表示第i个功能的失效概率;
步骤6-3、在修正后的一次系统设备失效概率基础上,进一步采用非序贯蒙特卡洛采样技术,对电网场景采样,并计算每个场景下的失负荷量;
步骤6-4、利用式(10)计算得到最终的失负荷期望值;
<mrow>
<mi>E</mi>
<mi>E</mi>
<mi>N</mi>
<mi>S</mi>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>S</mi>
</mrow>
</munder>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
pi为系统在状态i的概率;S为系统所有切负荷状态的集合,Ci为系统在状态i的失负荷量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710124484.4A CN107122878A (zh) | 2017-03-03 | 2017-03-03 | 一种基于大数据与人为因素的电力系统风险评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710124484.4A CN107122878A (zh) | 2017-03-03 | 2017-03-03 | 一种基于大数据与人为因素的电力系统风险评估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107122878A true CN107122878A (zh) | 2017-09-01 |
Family
ID=59718196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710124484.4A Pending CN107122878A (zh) | 2017-03-03 | 2017-03-03 | 一种基于大数据与人为因素的电力系统风险评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107122878A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108133083A (zh) * | 2017-12-08 | 2018-06-08 | 中国航天标准化研究所 | 一种空间人机系统安全性影响因素筛选方法 |
CN108985627A (zh) * | 2018-07-16 | 2018-12-11 | 浙江大学 | 一种综合考虑灾害和人因的电力系统可靠性评估算法 |
CN110019407A (zh) * | 2017-12-29 | 2019-07-16 | 北京金风科创风电设备有限公司 | 风电场中枚举数据的处理方法、装置、设备和存储介质 |
CN110633906A (zh) * | 2019-09-11 | 2019-12-31 | 绍兴建元电力集团有限公司 | 一种计及人因可靠性的变电站自动化方法失效风险评估方法 |
CN112395200A (zh) * | 2020-11-25 | 2021-02-23 | 电信科学技术第十研究所有限公司 | 一种软件可靠性操作剖面概率的确定方法及装置 |
CN112580916A (zh) * | 2019-09-30 | 2021-03-30 | 深圳无域科技技术有限公司 | 数据评估方法、装置、计算机设备和存储介质 |
TWI752273B (zh) * | 2018-03-28 | 2022-01-11 | 日商三菱電機股份有限公司 | 更新計畫作成支援裝置及更新計畫作成支援方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012109548A2 (en) * | 2011-02-10 | 2012-08-16 | Eaton Corporation | Equipment-related risk assessment from electrical event analysis |
CN102819813A (zh) * | 2012-08-20 | 2012-12-12 | 浙江大学 | 一种智能变电站自动化系统的安全风险评估方法 |
CN103679544A (zh) * | 2012-11-27 | 2014-03-26 | 江苏省电力公司南京供电公司 | 一种智能配电网运行综合评估方法 |
CN105470951A (zh) * | 2015-12-10 | 2016-04-06 | 浙江大学 | 一种基于大数据挖掘技术的配电网可靠供电能力评估方法 |
-
2017
- 2017-03-03 CN CN201710124484.4A patent/CN107122878A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012109548A2 (en) * | 2011-02-10 | 2012-08-16 | Eaton Corporation | Equipment-related risk assessment from electrical event analysis |
CN102819813A (zh) * | 2012-08-20 | 2012-12-12 | 浙江大学 | 一种智能变电站自动化系统的安全风险评估方法 |
CN103679544A (zh) * | 2012-11-27 | 2014-03-26 | 江苏省电力公司南京供电公司 | 一种智能配电网运行综合评估方法 |
CN105470951A (zh) * | 2015-12-10 | 2016-04-06 | 浙江大学 | 一种基于大数据挖掘技术的配电网可靠供电能力评估方法 |
Non-Patent Citations (2)
Title |
---|
俞斌: "基于功能分解的电力二次系统风险评估方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
唐俊熙: "电网操作人因可靠性及事故预控研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108133083A (zh) * | 2017-12-08 | 2018-06-08 | 中国航天标准化研究所 | 一种空间人机系统安全性影响因素筛选方法 |
CN108133083B (zh) * | 2017-12-08 | 2021-09-03 | 中国航天标准化研究所 | 一种空间人机系统安全性影响因素筛选方法 |
CN110019407A (zh) * | 2017-12-29 | 2019-07-16 | 北京金风科创风电设备有限公司 | 风电场中枚举数据的处理方法、装置、设备和存储介质 |
TWI752273B (zh) * | 2018-03-28 | 2022-01-11 | 日商三菱電機股份有限公司 | 更新計畫作成支援裝置及更新計畫作成支援方法 |
CN108985627A (zh) * | 2018-07-16 | 2018-12-11 | 浙江大学 | 一种综合考虑灾害和人因的电力系统可靠性评估算法 |
CN110633906A (zh) * | 2019-09-11 | 2019-12-31 | 绍兴建元电力集团有限公司 | 一种计及人因可靠性的变电站自动化方法失效风险评估方法 |
CN112580916A (zh) * | 2019-09-30 | 2021-03-30 | 深圳无域科技技术有限公司 | 数据评估方法、装置、计算机设备和存储介质 |
CN112580916B (zh) * | 2019-09-30 | 2024-05-28 | 深圳无域科技技术有限公司 | 数据评估方法、装置、计算机设备和存储介质 |
CN112395200A (zh) * | 2020-11-25 | 2021-02-23 | 电信科学技术第十研究所有限公司 | 一种软件可靠性操作剖面概率的确定方法及装置 |
CN112395200B (zh) * | 2020-11-25 | 2024-04-02 | 电信科学技术第十研究所有限公司 | 一种软件可靠性操作剖面概率的确定方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107122878A (zh) | 一种基于大数据与人为因素的电力系统风险评估方法 | |
Englin et al. | An empirical examination of optimal rotations in a multiple-use forest in the presence of fire risk | |
CN106602548B (zh) | 一种基于eeac的暂态稳定关键输电断面识别方法 | |
CN103618638B (zh) | 评估电力通信网检修方案的方法 | |
CN105139095A (zh) | 基于属性区间模型的配电网运行状态评估方法 | |
CN104408549A (zh) | 一种城市配电网运行状态评估方法 | |
CN104779608B (zh) | 基于离线策略的电力系统安全稳定在线附加紧急控制方法 | |
CN110619386A (zh) | 一种tmr运行监测及故障智能研判方法及系统 | |
CN103065057B (zh) | 一种电力系统差异化规划经济性评估方法 | |
CN105069692A (zh) | 准确的电网安全风险评估方法 | |
CN102915515A (zh) | 一种电网连锁性故障识别和风险评估方法 | |
CN107132762A (zh) | 基于预想故障集自动筛选的在线静态安全评估方法 | |
CN109633370A (zh) | 一种基于故障信息编码和融合方法的电网故障诊断方法 | |
CN106527394A (zh) | 一种考虑多重信息因素的ecps连锁故障风险评估方法 | |
CN103246939A (zh) | 基于安全稳定裕度的电网运行安全风险事件在线辨识方法 | |
CN112308438B (zh) | 基于状态耦合矩阵的能源互联系统可靠性评估方法及装置 | |
CN109001596A (zh) | 一种电网故障诊断与输电线路参数辨识系统 | |
CN106972517A (zh) | 基于双极对称特征的特高压直流输电系统可靠性计算方法 | |
CN105335819A (zh) | 一种基于大数据的信息系统风险预警模型构建方法 | |
CN113569411A (zh) | 一种面向灾害天气的电网运行风险态势感知方法 | |
CN112712205A (zh) | 一种基于长短期记忆神经网络的配电网故障预防方法 | |
EP3955012A1 (de) | Verfahren und einrichtung zum ermitteln des fehlerortes eines fehlers auf einer leitung eines elektrischen energieversorgungsnetzes | |
Zhou et al. | Robustness analysis of power system dynamic process and repair strategy | |
CN106355510A (zh) | 一种电力系统安全性判定方法 | |
CN104680431A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170901 |