CN108537413A - 基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 - Google Patents
基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 Download PDFInfo
- Publication number
- CN108537413A CN108537413A CN201810225114.4A CN201810225114A CN108537413A CN 108537413 A CN108537413 A CN 108537413A CN 201810225114 A CN201810225114 A CN 201810225114A CN 108537413 A CN108537413 A CN 108537413A
- Authority
- CN
- China
- Prior art keywords
- state
- time
- typhoon
- failure rate
- probability
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000012512 characterization method Methods 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 27
- 230000015556 catabolic process Effects 0.000 claims abstract description 14
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 10
- 238000005070 sampling Methods 0.000 claims abstract description 9
- 230000007257 malfunction Effects 0.000 claims abstract description 5
- 230000009897 systematic effect Effects 0.000 claims abstract description 5
- 238000012546 transfer Methods 0.000 claims abstract description 5
- 238000011156 evaluation Methods 0.000 claims abstract description 3
- 230000007704 transition Effects 0.000 claims description 19
- 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
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000004048 modification Effects 0.000 claims description 3
- 238000012986 modification Methods 0.000 claims description 3
- 230000008707 rearrangement Effects 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 230000007306 turnover 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/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于马尔科夫链的考虑台风影响时空特性的电网韧性评估方法,包括:步骤一、基于传统的元件脆弱性曲线,并考虑台风天气对电网影响的时空特性,计算在台风天气下的时变故障率;步骤二、采用抽样方法,根据元件故障率的时变曲线,确定***中受台风影响元件的可能故障时间和相应的故障概率;步骤三、基于马尔科夫链,根据各元件的可能故障时间及相应的故障概率,计算各元件的可能故障时间所对应的***状态转移概率矩阵,从而得出各时段下的***状态概率向量;步骤四、根据各时段下的***状态概率向量,对各时段***的可能状态进行韧性水平评价,并应用蒙特卡洛方法对***的韧性水平进行期望评估。
Description
技术领域
本发明属于电力***分析领域,尤其涉及一种电网韧性评估方法。
背景技术
电力***在正常运行条件下,需要对预想事故集满足相应的安全可靠准则。现代电力***一般能够满足相应的可靠性指标。然而,当遇到极端天气时,电力***却显得极为脆弱。究其原因就在于,电力***的可靠性一般关注的是经常发生的、对电网影响较小的事故,而极端天气发生的概率往往很低,对电网造成的影响却极为严重。目前用韧性衡量电力***的这种准备应对并适应不断改变条件的能力以及承受扰动并从中快速恢复的能力。对电力***的韧性进行评估,可以得知电力***承受自然威胁的能力,以便及时采取适当的韧性提升措施,减少灾害来临时对电网造成的威胁。韧性评估存在于极端天气对***产生影响的整个过程中,而处于***中不同位置的元件,在不同时刻所经历的极端天气强度均是不同的,因此元件故障率的变化是受台风影响的具有时空特性的曲线。针对实时变化的故障率,***在不同时刻也可能会有不同的状态,在韧性评估时需要综合计及***状态发生改变的各可能时间和对应阶段存在的各种状态,以分析各个可能阶段***的韧性水平。综上所述,建立基于马尔科夫链的考虑台风影响时空特性的韧性评估方法,反映电力***在台风天气影响的整个过程中的韧性水平,是亟待解决的实际问题,具有良好的理论价值和应用价值。
发明内容
针对台风天气对于***影响的时空特性所带来的元件时变故障率和台风天气对***元件产生影响的整个过程中的韧性评估问题,本发明提供一种基于马尔科夫链的考虑台风时空特性的电网韧性评估方法。
为了解决上述技术问题,本发明提出的一种基于马尔科夫链的考虑台风时空特性的电网韧性评估方法,包括以下步骤:
步骤一、基于传统的元件脆弱性曲线,并考虑台风天气对电网影响的时空特性,计算受台风影响的元件在台风天气下的时变故障率;
步骤二、根据步骤一得到的时变故障率绘制出时变曲线,并采用抽样方法确定受台风影响元件的可能故障时间和相应的故障率;
步骤三、基于马尔科夫链,根据各元件的可能故障时间及相应的故障率,计算各元件的可能故障时间所对应的***状态转移概率矩阵,从而得出各时段下的***状态概率向量;
步骤四、根据各时段下的***状态概率向量,对各时段***的可能状态进行韧性水平评价,并应用蒙特卡洛方法对***的韧性水平进行期望评估,最终得到***可满足的期望负荷量百分比。
进一步讲,步骤一的具体内容包括:
步骤1)根据元件的脆弱性曲线建立元件故障率与台风强度的关系模型:
式(1)中,p为元件故障率,v为台风风速,vmax为元件可以承受的最大台风风速,vmin表示台风风速大于此风速时,元件的故障率开始大于元件在正常运行时的故障率;
步骤2)不同地理位置的元件在不同时刻所经历的台风强度为:
v=v(t) (2)
式(2)中,v(t)表示不同元件在不同时间经历的台风风速,由公式(1)和(2)可得,元件的时变故障率为:
p=p_time(t) (3)
步骤二的具体内容包括:
设***中受台风影响的元件按照受台风影响的先后顺序依次为元件1,元件2,…,元件i,…,元件m,其总数为m个;其中,元件i的故障率随时间的变化曲线为:
pi=p_timei(t) t∈(ts,te) (4)
式(4)中,ts和te分别为元件i受影响的起止时间;
在ts至te时间内对于元件i进行故障时刻抽样:在元件i的受影响时间内,即ts至te,首先,产生随机时间t,并产生与t对应的0-1之间的随机数r;t时刻元件i对应的故障率为pi,t,如果r<pi,t,那么元件i可能在此时刻故障,且故障率为pi,t;如果r>pi,t,则再次产生新的随机数,并进行相同的比较,直到产生相应的故障时刻和故障率。
步骤三的具体内容包括:
步骤1)记元件1,2,…,m的抽样时间分别为T1′,T2′,…,T′m,对所述抽样时间按照前后顺序重新排列为T1,T2,…,Tm,同时按重新排序的对应顺序调整各元件的顺序;
步骤2)不考虑故障时刻,计入所有元件抽样得到的故障率,按照下述方法求解得到状态转移概率矩阵P:
记元件的状态转移率矩阵为:
式(5)中,λ为元件故障率;
设***中受影响的元件数为m,定义***的状态s是由每个元件的状态组成的m维数组,数组中各元件的状态仍然用0和1表示;***的所有状态集合为S,s∈S;状态集合S中共包含2m个状态;
设***的状态转移矩阵为P,且状态x,y∈S,其元素P[x,y]表示状态x向状态y转移的概率,则P[x,y]的计算公式有下述三种情形之一:
(1)当状态x和状态y对应的m维数组存在两个或两个以上不同元素时:
P[x,y]=0 (6)
(2)当状态x和状态y对应的m维数组有一个元素相异时:
式(7)中,i表示元件编号,Tixy表示状态x和状态y对应的m维数组中,元件i由状态x向状态y转移的概率;
(3)当状态x和状态y对应的m维数组没有元素相异时,即状态x向自身转移时:
步骤3)通过对状态转移概率矩阵P进行如下修改得到Ti时刻的状态转移概率矩阵PTi,
(1)对于PTi的非对角元素,即不同状态间的转换概率:
当状态x向状态y转移时,若数组中存在元件i的状态由0变为1,即由正常运行状态变为故障状态,且数组中其余位置元素不变,则元素PTi[x,y]=P[x,y];位于矩阵PTi第x行的其余元素均置0;
(2)对于PTi的对角元素,即状态向自身转换的概率,
由此求出Ti时刻的状态转移概率矩阵PTi;
步骤4)Ti时刻的状态概率向量为:
式(10)中,State0为2m维初始状态概率向量,其元素代表了相应状态的初始概率,STi中各元素代表各状态存在的概率。
步骤四的具体内容包括:
步骤1)采用蒙特卡洛仿真,在第n次仿真中,共包含m个可能的故障时间,每个可能故障时间Ti下,对应的该时刻状态概率向量为STi;每个故障时间后,***的期望失负荷量为:
RTi=STi×L i=1,2,...,m (11)
式(11)中,L为每个***状态对应的失负荷量,RTi为每个阶段的期望失负荷量;
步骤2)对于第n次仿真,其最终期望失负荷量为:
Rn=RTm n=1,2,...,N (12)
步骤3)蒙特卡洛仿真最终结果为:
式(13)中RMC表示***最终的期望失负荷量;
步骤4)***的最终韧性水平为:
式(14)中,R0表示***受影响的负荷量,RE即为***可满足的期望负荷量百分比。
与现有技术相比,本发明的有益效果是:
本发明提出一种基于马尔科夫链的考虑台风影响时空特性的电网韧性评估方法,本发明方法可以计及在台风对电网产生影响的整个过程中,***中多个受影响元件的时空特性故障率变化曲线,通过提出的抽样方法,对元件可能的故障时间进行抽样,在元件故障率实时变化的情况下生成***状态变化的可能时间和相应的***状态,并基于马尔科夫模型对各状态概率求解,得到大量情景下的***各阶段期望失负荷量,用于韧性评估,所得结果为***各可能阶段的失负荷量和多次仿真下的最终期望失负荷量,可以反映各场景下***各阶段的韧性期望水平和综合的韧性期望水平。
附图说明
图1元件在台风天气下的脆弱性曲线;
图2元件故障时刻随机抽样流程图;
图3韧性水平评估流程图;
图4IEEE33节点***及受影响线路示意图;
图5元件故障率示意图;
图6 N=10时不同仿真中各时段的期望失负荷量;
图7 N=2000时***的期望失负荷量结果。
具体实施方式
下面结合附图和具体实施算例对本发明技术方案作进一步详细描述。
本发明提出的一种基于马尔科夫链的考虑台风时空特性的电网韧性评估方法,下面采用如图4所示的IEEE33节点测试***验证本发明所提方法的有效性与正确性。图中X符号示出了台风经过路径。台风依次经过线路2-19,线路5-23,线路26-27和线路6-7。
借助matlab编程实现对该电网韧性评估方法中的相关算法,计算机的配置为:Intel Core i3-2120处理器,6G内存。
基于马尔科夫链的考虑台风时空特性对本发明中实施例IEEE33节点测试***进行电网韧性评估方法,包括以下步骤:
步骤一、基于传统的元件脆弱性曲线,并考虑台风天气对电网影响的时空特性,计算受台风影响的元件在台风天气下的时变故障率;具体内容包括:
步骤1)根据元件的脆弱性曲线,如图1所示,建立元件故障率与台风强度的关系模型:
式(1)中,p为元件故障率,v为台风风速,vmax为元件可以承受的最大台风风速,vmin表示台风风速大于此风速时,元件的故障率开始大于元件在正常运行时的故障率;
步骤2)不同地理位置的元件在不同时刻所经历的台风强度为:
v=v(t) (2)
式(2)中,v(t)表示不同元件在不同时间经历的台风风速,由公式(1)和(2)可得,元件的时变故障率为:
p=p_time(t) (3)
针对实施例,台风对元件的时间空间影响表现为图5中不同元件的故障率随时间的变化曲线(自左至右依次为线路2-19,线路5-23,线路26-27和线路6-7)。
步骤二、根据步骤一得到的时变故障率绘制出时变曲线,并采用抽样方法确定受台风影响元件的可能故障时间和相应的故障率;如图2所示,具体内容包括:
设***中受台风影响的元件按照受台风影响的先后顺序依次为元件1,元件2,…,元件i,…,元件m,其总数为m个;其中,元件i的故障率随时间的变化曲线为:
pi=p_timei(t)t∈(ts,te) (4)
式(4)中,ts和te分别为元件i受影响的起止时间;
在ts至te时间内对于元件i进行故障时刻抽样:在元件i的受影响时间内,即ts至te,首先,产生随机时间t,并产生与t对应的0-1之间的随机数r;t时刻元件i对应的故障率为pi,t,如果r<pi,t,那么元件i可能在此时刻故障,且故障率为pi,t;如果r>pi,t,则再次产生新的随机数,并进行相同的比较,直到产生相应的故障时刻和故障率。
步骤三、基于马尔科夫链,根据各元件的可能故障时间及相应的故障率,计算各元件的可能故障时间所对应的***状态转移概率矩阵,从而得出各时段下的***状态概率向量;如图3所示,具体内容包括:
步骤1)记元件1,2,…,m的抽样时间分别为T1′,T2′,…,T′m,对所述抽样时间按照前后顺序重新排列为T1,T2,…,Tm,同时按重新排序的对应顺序调整各元件的顺序;
步骤2)不考虑故障时刻,计入所有元件抽样得到的故障率,按照下述方法求解得到状态转移概率矩阵P:
记元件的状态转移率矩阵为:
式(5)中,λ为元件故障率;
设***中受影响的元件数为m,定义***的状态s是由每个元件的状态组成的m维数组,数组中各元件的状态仍然用0和1表示;***的所有状态集合为S,s∈S;状态集合S中共包含2m个状态;
设***的状态转移矩阵为P,且状态x,y∈S,其元素P[x,y]表示状态x向状态y转移的概率,则P[x,y]的计算公式有下述三种情形之一:
(4)当状态x和状态y对应的m维数组存在两个或两个以上不同元素时:
P[x,y]=0 (6)
(5)当状态x和状态y对应的m维数组有一个元素相异时:
式(7)中,i表示元件编号,Tixy表示状态x和状态y对应的m维数组中,元件i由状态x向状态y转移的概率;
(6)当状态x和状态y对应的m维数组没有元素相异时,即状态x向自身转移时:
步骤3)通过对状态转移概率矩阵P进行如下修改得到Ti时刻的状态转移概率矩阵PTi,
(1)对于PTi的非对角元素,即不同状态间的转换概率:
当状态x向状态y转移时,若数组中存在元件i的状态由0变为1,即由正常运行状态变为故障状态,且数组中其余位置元素不变,则元素PTi[x,y]=P[x,y];位于矩阵PTi第x行的其余元素均置0;
(2)对于PTi的对角元素,即状态向自身转换的概率:
由此求出Ti时刻的状态转移概率矩阵PTi;
步骤4)Ti时刻的状态概率向量为:
式(10)中,State0为2m维初始状态概率向量,其元素代表了相应状态的初始概率,STi中各元素代表各状态存在的概率。
针对实施例,设定***初始状态为State0=[1 0 ... 0]1×16,即所有线路初始正常运行,计算每个时刻的状态转移概率矩阵,然后,计算各个时段的状态概率矩阵,其中各个阶段的期望失负荷量如图6所示,不同的曲线代表不同仿真中所得的各阶段期望失负荷量,该期望失负荷量为该阶段所有可能***状态及其后果的期望值,同一曲线不同的转折点表示一次仿真中不同线路所对应的不同可能故障时间。
步骤四、根据各时段下的***状态概率向量,对各时段***的可能状态进行韧性水平评价,并应用蒙特卡洛方法对***的韧性水平进行期望评估,最终得到***可满足的期望负荷量百分比。如图3所示,具体内容包括:
步骤1)采用蒙特卡洛仿真,在第n次仿真中,共包含m个可能的故障时间,每个可能故障时间Ti下,对应的该时刻状态概率向量为STi;每个故障时间后,***的期望失负荷量为:
RTi=STi×L i=1,2,...,m (11)
式(11)中,L为每个***状态对应的失负荷量,RTi为每个阶段的期望失负荷量;
步骤2)对于第n次仿真,其最终期望失负荷量为:
Rn=RTm n=1,2,...,N (12)
步骤3)蒙特卡洛仿真最终结果为:
式(13)中RMC表示***最终的期望失负荷量;
步骤4)***的最终韧性水平为:
式(14)中,R0表示***受影响的负荷量,RE即为***可满足的期望负荷量百分比。
针对实施例,进行多次仿真,得到仿真结果最终如图7所示,由图7可知,仿真结果最终趋***直,为239.5MW。因此***的韧性水平以预计失负荷量表示为239.5MW。用百分比表示为:
尽管上面结合附图对本发明进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨的情况下,还可以做出很多变形,这些均属于本发明的保护之内。
Claims (5)
1.一种基于马尔科夫链的考虑台风时空特性的电网韧性评估方法,其特征在于,所述方法包括以下步骤:
步骤一、基于传统的元件脆弱性曲线,并考虑台风天气对电网影响的时空特性,计算受台风影响的元件在台风天气下的时变故障率;
步骤二、根据步骤一得到的时变故障率绘制出时变曲线,并采用抽样方法确定受台风影响元件的可能故障时间和相应的故障率;
步骤三、基于马尔科夫链,根据各元件的可能故障时间及相应的故障率,计算各元件的可能故障时间所对应的***状态转移概率矩阵,从而得出各时段下的***状态概率向量;
步骤四、根据各时段下的***状态概率向量,对各时段***的可能状态进行韧性水平评价,并应用蒙特卡洛方法对***的韧性水平进行期望评估,最终得到***可满足的期望负荷量百分比。
2.根据权利要求1所述的一种考虑台风影响时空特性的电网韧性评估方法,其中,步骤一的具体内容包括:
步骤1)根据元件的脆弱性曲线建立元件故障率与台风强度的关系模型:
式(1)中,p为元件故障率,v为台风风速,vmax为元件可以承受的最大台风风速,vmin表示台风风速大于此风速时,元件的故障率开始大于元件在正常运行时的故障率;
步骤2)不同地理位置的元件在不同时刻所经历的台风强度为:
v=v(t) (2)
式(2)中,v(t)表示不同元件在不同时间经历的台风风速,由公式(1)和(2)可得,元件的时变故障率为:
p=p_time(t) (3)。
3.根据权利要求1所述的一种考虑台风影响时空特性的电网韧性评估方法,其中,步骤二的具体内容包括:
设***中受台风影响的元件按照受台风影响的先后顺序依次为元件1,元件2,…,元件i,…,元件m,其总数为m个;其中,元件i的故障率随时间的变化曲线为:
pi=p_timei(t) t∈(ts,te) (4)
式(4)中,ts和te分别为元件i受影响的起止时间;
在ts至te时间内对于元件i进行故障时刻抽样:在元件i的受影响时间内,即ts至te,首先,产生随机时间t,并产生与t对应的0-1之间的随机数r;t时刻元件i对应的故障率为pi,t,如果r<pi,t,那么元件i可能在此时刻故障,且故障率为pi,t;如果r>pi,t,则再次产生新的随机数,并进行相同的比较,直到产生相应的故障时刻和故障率。
4.根据权利要求1所述的一种考虑台风影响时空特性的电网韧性评估方法,其中,步骤三的具体内容包括:
步骤1)记元件1,2,…,m的抽样时间分别为T1′,T′2,…,T′m,对所述抽样时间按照前后顺序重新排列为T1,T2,…,Tm,同时按重新排序的对应顺序调整各元件的顺序;
步骤2)不考虑故障时刻,计入所有元件抽样得到的故障率,按照下述方法求解得到状态转移概率矩阵P:
记元件的状态转移率矩阵为:
式(5)中,λ为元件故障率;
设***中受影响的元件数为m,定义***的状态s是由每个元件的状态组成的m维数组,数组中各元件的状态仍然用0和1表示;***的所有状态集合为S,s∈S;状态集合S中共包含2m个状态;
设***的状态转移矩阵为P,且状态x,y∈S,其元素P[x,y]表示状态x向状态y转移的概率,则P[x,y]的计算公式有下述三种情形之一:
(1)当状态x和状态y对应的m维数组存在两个或两个以上不同元素时:
P[x,y]=0 (6)
(2)当状态x和状态y对应的m维数组有一个元素相异时:
式(7)中,i表示元件编号,Tixy表示状态x和状态y对应的m维数组中,元件i由状态x向状态y转移的概率;
(3)当状态x和状态y对应的m维数组没有元素相异时,即状态x向自身转移时:
步骤3)通过对状态转移概率矩阵P进行如下修改得到Ti时刻的状态转移概率矩阵PTi,
(1)对于PTi的非对角元素,即不同状态间的转换概率:
当状态x向状态y转移时,若数组中存在元件i的状态由0变为1,即由正常运行状态变为故障状态,且数组中其余位置元素不变,则元素PTi[x,y]=P[x,y];位于矩阵PTi第x行的其余元素均置0;
(2)对于PTi的对角元素,即状态向自身转换的概率,
由此求出Ti时刻的状态转移概率矩阵PTi;
步骤4)Ti时刻的状态概率向量为:
式(10)中,State0为2m维初始状态概率向量,其元素代表了相应状态的初始概率,STi中各元素代表各状态存在的概率。
5.根据权利要求1所述的一种考虑台风影响时空特性的电网韧性评估方法,其中,步骤四的具体内容包括:
步骤1)采用蒙特卡洛仿真,在第n次仿真中,共包含m个可能的故障时间,每个可能故障时间Ti下,对应的该时刻状态概率向量为STi;每个故障时间后,***的期望失负荷量为:
RTi=STi×L i=1,2,...,m (11)
式(11)中,L为每个***状态对应的失负荷量,RTi为每个阶段的期望失负荷量;
步骤2)对于第n次仿真,其最终期望失负荷量为:
Rn=RTm n=1,2,...,N (12)
步骤3)蒙特卡洛仿真最终结果为:
式(13)中RMC表示***最终的期望失负荷量;
步骤4)***的最终韧性水平为:
式(14)中,R0表示***受影响的负荷量,RE即为***可满足的期望负荷量百分比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810225114.4A CN108537413B (zh) | 2018-03-19 | 2018-03-19 | 基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810225114.4A CN108537413B (zh) | 2018-03-19 | 2018-03-19 | 基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108537413A true CN108537413A (zh) | 2018-09-14 |
CN108537413B CN108537413B (zh) | 2021-09-10 |
Family
ID=63483896
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810225114.4A Expired - Fee Related CN108537413B (zh) | 2018-03-19 | 2018-03-19 | 基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108537413B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109274124A (zh) * | 2018-11-22 | 2019-01-25 | 国网黑龙江省电力有限公司电力科学研究院 | 基于场景马尔科夫法的风电就地消纳能力预测方法 |
CN109784780A (zh) * | 2019-03-04 | 2019-05-21 | 广东工业大学 | 一种电力***韧性的评估方法、装置及设备 |
CN110233832A (zh) * | 2019-05-22 | 2019-09-13 | 中国电子科技集团公司第二十八研究所 | 韧性信息服务云环境试验评估***及其使用方法 |
CN110535144A (zh) * | 2019-09-27 | 2019-12-03 | 国网甘肃省电力公司经济技术研究院 | 风沙天气下含多类型负荷的智能配电网韧性定量分析方法 |
CN110570122A (zh) * | 2019-09-10 | 2019-12-13 | 重庆大学 | 一种计及风速季节特性及集电***元件故障的海上风电场可靠性评估方法 |
CN111507606A (zh) * | 2020-04-13 | 2020-08-07 | 华北电力大学 | 用于复杂能源互联***的韧性评估方法 |
CN111582512A (zh) * | 2020-03-31 | 2020-08-25 | 清华大学深圳国际研究生院 | 一种配电网韧性恢复方法及计算机可读存储介质 |
CN112348697A (zh) * | 2020-10-21 | 2021-02-09 | 国网天津市电力公司 | 基于大数据的电网运行状态综合评估方法及装置 |
CN113872191A (zh) * | 2021-09-24 | 2021-12-31 | 武汉大学 | 一种极端天气下电力***弹性提升方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110231028A1 (en) * | 2009-01-14 | 2011-09-22 | Ozog Michael T | Optimization of microgrid energy use and distribution |
CN105260952A (zh) * | 2015-11-19 | 2016-01-20 | 华北电力大学(保定) | 基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 |
CN107681691A (zh) * | 2017-09-30 | 2018-02-09 | 太原理工大学 | 计及不确定性因素的风电并网***运行可靠性评估方法 |
CN107769202A (zh) * | 2017-11-28 | 2018-03-06 | 西华大学 | 计及天气模型的基于模糊Petri网的配电网可靠性评估方法 |
-
2018
- 2018-03-19 CN CN201810225114.4A patent/CN108537413B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110231028A1 (en) * | 2009-01-14 | 2011-09-22 | Ozog Michael T | Optimization of microgrid energy use and distribution |
CN105260952A (zh) * | 2015-11-19 | 2016-01-20 | 华北电力大学(保定) | 基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 |
CN107681691A (zh) * | 2017-09-30 | 2018-02-09 | 太原理工大学 | 计及不确定性因素的风电并网***运行可靠性评估方法 |
CN107769202A (zh) * | 2017-11-28 | 2018-03-06 | 西华大学 | 计及天气模型的基于模糊Petri网的配电网可靠性评估方法 |
Non-Patent Citations (1)
Title |
---|
陈小青: "基于蒙特卡洛模拟的电网调度运行风险评估研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109274124B (zh) * | 2018-11-22 | 2022-08-12 | 国网黑龙江省电力有限公司电力科学研究院 | 基于场景马尔科夫法的风电就地消纳能力预测方法 |
CN109274124A (zh) * | 2018-11-22 | 2019-01-25 | 国网黑龙江省电力有限公司电力科学研究院 | 基于场景马尔科夫法的风电就地消纳能力预测方法 |
CN109784780A (zh) * | 2019-03-04 | 2019-05-21 | 广东工业大学 | 一种电力***韧性的评估方法、装置及设备 |
CN109784780B (zh) * | 2019-03-04 | 2022-12-30 | 广东工业大学 | 一种电力***韧性的评估方法、装置及设备 |
CN110233832A (zh) * | 2019-05-22 | 2019-09-13 | 中国电子科技集团公司第二十八研究所 | 韧性信息服务云环境试验评估***及其使用方法 |
CN110570122A (zh) * | 2019-09-10 | 2019-12-13 | 重庆大学 | 一种计及风速季节特性及集电***元件故障的海上风电场可靠性评估方法 |
CN110570122B (zh) * | 2019-09-10 | 2023-03-14 | 重庆大学 | 一种计及风速季节特性及集电***元件故障的海上风电场可靠性评估方法 |
CN110535144A (zh) * | 2019-09-27 | 2019-12-03 | 国网甘肃省电力公司经济技术研究院 | 风沙天气下含多类型负荷的智能配电网韧性定量分析方法 |
CN111582512A (zh) * | 2020-03-31 | 2020-08-25 | 清华大学深圳国际研究生院 | 一种配电网韧性恢复方法及计算机可读存储介质 |
CN111582512B (zh) * | 2020-03-31 | 2023-04-28 | 清华大学深圳国际研究生院 | 一种配电网韧性恢复方法及计算机可读存储介质 |
CN111507606B (zh) * | 2020-04-13 | 2022-05-20 | 华北电力大学 | 用于复杂能源互联***的韧性评估方法 |
CN111507606A (zh) * | 2020-04-13 | 2020-08-07 | 华北电力大学 | 用于复杂能源互联***的韧性评估方法 |
CN112348697A (zh) * | 2020-10-21 | 2021-02-09 | 国网天津市电力公司 | 基于大数据的电网运行状态综合评估方法及装置 |
CN113872191A (zh) * | 2021-09-24 | 2021-12-31 | 武汉大学 | 一种极端天气下电力***弹性提升方法 |
CN113872191B (zh) * | 2021-09-24 | 2023-02-24 | 武汉大学 | 一种极端天气下电力***弹性提升方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108537413B (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108537413A (zh) | 基于马尔科夫链的考虑台风时空特性的电网韧性评估方法 | |
CN104809722B (zh) | 一种基于红外热像的电气设备故障诊断方法 | |
CN103746370B (zh) | 一种风电场可靠性建模方法 | |
CN106101252A (zh) | 基于大数据和可信计算的信息安全风险防护*** | |
CN108693771A (zh) | 一种基于多种群遗传算法的配电网故障区段定位算法 | |
CN111507509A (zh) | 一种电力***极端事件的风险评估方法 | |
CN104133143B (zh) | 一种基于Hadoop云计算平台的电网线路故障诊断***及方法 | |
CN104331572A (zh) | 考虑风速和风电机组故障相关性的风电场可靠性建模方法 | |
CN105278460A (zh) | 基于级联故障分析的数控机床***组件可靠性评价方法 | |
CN102768701A (zh) | 基于量子遗传算法的高压开关柜绝缘子电场优化方法 | |
CN112381282B (zh) | 基于宽度学习***的光伏发电功率预测方法 | |
CN111191832A (zh) | 一种台风灾害配电网杆塔故障预测方法及*** | |
CN103887795A (zh) | 电力***实时概率潮流在线的计算方法 | |
CN106684885A (zh) | 基于多场景分析的风电机组配电网无功优化方法 | |
CN113569411A (zh) | 一种面向灾害天气的电网运行风险态势感知方法 | |
CN106655201A (zh) | 一种基于安全域的电力***热稳定安全优化控制方法 | |
CN104779608A (zh) | 基于离线策略的电力***安全稳定在线附加紧急控制方法 | |
CN115526333A (zh) | 边缘场景下动态权重的联邦学习方法 | |
Ma et al. | A measurement error prediction framework for smart meters under extreme natural environment stresses | |
CN116167609A (zh) | 基于神经网络模型的电力***风险评估方法 | |
CN115935407A (zh) | 一种面向联邦学习的自适应差分隐私保护方法 | |
CN103347028A (zh) | 云架构下基于贝叶斯的对等网络信任度量模型 | |
CN113065218B (zh) | 考虑lr攻击的电力***可靠性评估方法、装置及*** | |
Rouhani et al. | Soft clustering based probabilistic power flow with correlated inter temporal events | |
CN106056305A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210910 |
|
CF01 | Termination of patent right due to non-payment of annual fee |