CN107872328B - 基于通信链路多种状态计算iec61850通信系统失效率的方法 - Google Patents

基于通信链路多种状态计算iec61850通信系统失效率的方法 Download PDF

Info

Publication number
CN107872328B
CN107872328B CN201610840108.0A CN201610840108A CN107872328B CN 107872328 B CN107872328 B CN 107872328B CN 201610840108 A CN201610840108 A CN 201610840108A CN 107872328 B CN107872328 B CN 107872328B
Authority
CN
China
Prior art keywords
state
communication system
calculating
link
communication
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.)
Expired - Fee Related
Application number
CN201610840108.0A
Other languages
English (en)
Other versions
CN107872328A (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.)
State Grid Information and Telecommunication Co Ltd
North China Electric Power University
Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Information and Telecommunication Co Ltd
North China Electric Power University
Nanjing Power Supply Co of State Grid Jiangsu 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 State Grid Information and Telecommunication Co Ltd, North China Electric Power University, Nanjing Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Information and Telecommunication Co Ltd
Priority to CN201610840108.0A priority Critical patent/CN107872328B/zh
Publication of CN107872328A publication Critical patent/CN107872328A/zh
Application granted granted Critical
Publication of CN107872328B publication Critical patent/CN107872328B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
    • H04L41/14—Network analysis or design
    • H04L41/142—Network analysis or design using statistical or mathematical methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00—Computing arrangements based on specific mathematical models
    • G06N7/01—Probabilistic graphical models, e.g. probabilistic networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Probability & Statistics with Applications (AREA)
  • Algebra (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Maintenance And Management Of Digital Transmission (AREA)

Abstract

基于通信链路多种状态计算IEC61850通信系统失效率的方法，创新点在于形成通信系统故障状态空间和提出通信系统状态序列类的概念。所述方法包括步骤：A计算单条通信链路的最大允许连续数据集丢失数量；B定义单条通信链路的状态，计算通信链路的状态数量；C计算通信系统故障状态空间数量r，并将状态G(1)~G(r)对应于具体链路状态，且规定通信系统各链路都正常的状态记为G(1)；D计算通信系统状态一步转移概率矩阵P (1) ；E按照首次出现系统传送有效的时刻特征对通信系统状态分类；F计算各状态序列类的概率；G计算通信系统规定时间内的失效率。本发明所述方法能够对通信系统有效性简化计算，实现容易，精度高。

Description

基于通信链路多种状态计算IEC61850通信系统失效率的方法
技术领域
本发明涉及电力系统通信及系统可靠性计算技术领域，特别是涉及IEC61850通信系统领域的可靠性计算方法，属IEC61850通信系统可靠性领域。
背景技术
为了适应智能变电站的建设，智能电子设备在变电站得到大规模应用。智能化一次设备可以实现状态信息的数字化采集和网络传输，从而建立基于IEC 61850的统一通信系统和通信平台。IEC 61850通信系统的状态监测是变电站中智能一次设备的核心组成部分，IEC61850的通信系统是一个分布式实时系统，其功能需要多个IED通过交换信息共同协作，并在规定时间内完成功能。对IEC61850信息系统可靠性的计算，涉及到信息系统在规定时间内是否顺利完成特定功能，其计算和评估方法在当前阶段研究较少。
IEC61850信息系统可靠性，是指对特定功能在规定时间内完成的可能性的考察。针对此思想，有研究者提出了基于离散时间马尔科夫链的可靠性评估方法。此类方法的基本思想如下：（1）以IED为顶点或通信链路为基本单元构成通信系统，根据失效的基本单元不同，列举出所有的通信系统状态（包括正常和故障状态）；（2）根据各基本单元的失效率经验数据，可计算出各通信系统状态之间的概率转移矩阵；（3）再以概率转移矩阵为基础，计算出在规定的时间内仍处于正常状态的概率，即可得到通信系统在规定时间内的可用性。
由上可知，此类方法计算结果的合理性，取决于各基本单元失效率的确定是否合理。在已提出的方法中，相关的失效率数据一般取自于二次设备厂商的维护数据。
这样的失效率数据选取方法，对于通信链路失效率的确定可能存在一定问题。因为维护数据所给出的失效率，对应的都是不可修复故障。而实际上，通信链路在传输两逻辑节点间的数据集时，发生的故障是多种多样的。既有不可修复的故障，也有随机发生的零散数据集丢失（协议、缓存等原因）。而当后者发生时，IEC61850通信系统仍有可能在规定时间完成特定功能（当零散数据集丢失情况不严重时）。
因此，通信链路的故障情况应该是多样的。在形成通信系统故障状态空间时，应充分考虑这一因素。
发明内容
鉴于以上情况，本专利在基于离散时间马尔科夫链的可靠性计算方法的基础上，提出了一种充分考虑通信链路各种故障状态的IEC61850通信系统失效率计算方法。本发明定义了通信链路的各种状态，在链路状态的基础上构成通信系统的状态空间，通过构成通信系统的各链路的状态空间结构，即可形成通信系统状态变化的转移概率矩阵；其次，本发明提出了通信系统状态序列类的概念，将庞大的状态转换数量分为较少的几类，每一类状态序列对应一个通信系统有效的序列集合，可更简化的计算在规定时间内的通信系统失效率，使得本发明涉及的方法对于状态空间较大的情景仍能够有很好的通用性。
由于本技术方案的改进主要体现在对通信链路失效率的确定，因而在计算系统的失效率时，为计算简便，认为各IED不发生硬件失效。
为实现此目的，本专利采用的技术方案如下。
基于通信链路多种状态计算IEC61850通信系统失效率的方法，所述方法包括步骤：（A）计算单条通信链路的最大允许连续丢包数量n；（B）定义单条通信链路的状态，计算通信链路可能存在的状态数量；（C）计算出通信系统故障状态空间的状态数量r，并规定通信系统各条链路都正常的状态记为G(1)，其余的状态分别记为G(2)~G(r) ，并将状态G(1)~G(r)一一对应于具体链路状态；（D）根据通信链路状态之间的转移概率为p j 计算通信系统状态一步转移概率矩阵P (1) ；（E）根据通信系统所有可能的状态序列中，首次出现系统传送有效的时刻特征对通信系统状态进行分类；（F）按照首中概率计算方法，计算各状态序列类的概率；（G）计算在规定时间内通信系统的失效率。
在所述的步骤B中，关于信息链路状态的定义如下：
将事件“某固定时刻通信链路i将要发生连续j个时间间隔（包含当前时刻）的传送故障”记为状态L i (j)
该通信链路i在某一时刻可能存在的状态数量为m+1个，其中m为自该时刻至规定时间结束剩余的采样次数（1≤m≤n+1），故事件“通信链路i在当前时刻正常”的状态应记为L i (0)
步骤E中对状态序列分为若干个状态序列类，其步骤包括：
E1、定义初始时刻（0时刻）系统为正常，然后按照发送数据的时间顺序形成通信系统的状态，每一个按照时间顺序排列的固定顺序称之为一个“状态序列”。每一序列均拥有n +2个状态，其中第一个（0时刻的）状态为G(1)。故通信系统拥有rn+1个状态序列。
E2、对rn+1个状态序列进行分类，为使通信系统有效，故要求在所有序列类中，出现至少一个正常状态G(1)的即为系统有效，故按照在0时刻后首次出现状态G(1)状态的顺序，对所有含状态G(1)的状态序列进行分类，可分为n+1个状态序列。
而对于不包含G(1)的状态序列，由于该状态序列对应于系统无效，不在状态序列类的分类之中。
E3、对每一状态序列类，将其中的G(1)状态表示为1，确定为非G(1)的状态1表示为0，其余不确定的状态用“*”表示。
在以上规则基础上，将通信链路发生故障时的数据集丢失的时间长度，作为区分通信链路各类故障的依据，各类故障都有一定的发生概率。以此为基础，形成通信系统的全部故障状态空间，并计算各状态之间的一步转移概率矩阵。在计算规定时间内通信系统的失效率时，假设只要有一个时刻各条链路均处于正常状态，就认为通信系统没有失效。基于以上链路通信可用的实际情况，定义了状态序列类的概念，并将其用于通信系统在规定时间内的失效率计算中。
附图说明
图1为本发明实施方式中基于通信链路多种状态计算IEC61850通信系统失效率的流程图。
图2为本发明中一个单条链路的可能状态示意图。
图3为本发明中一个包含3条链路的通信系统示意图，且该图中每条通信链路的可能状态均如图2所示。
图4为对应于图3通信系统在某一时刻形成的状态空间。
具体实施方式
下面结合附图，对本发明作详细说明。
以下公开详细的示范实施例。然而，此处公开的具体步骤和功能细节仅仅是出于描述示范实施例的目的。
然而，应该理解，本发明不局限于公开的具体示范实施例，而是覆盖落入本公开范围内的所有修改、等同物和替换物。在对全部附图的描述中，相同的附图标记表示相同的元件。
同时应该理解，如在此所用的术语“链路”表示一条单独的通信传送渠道，“通信系统”表示多条链路通过一定方式连接成的一个通信单元。“和/或”包括一个或多个相关的列出项的任意和所有组合。此外，用来描述部件或单元之间关系的其他词语应该按照相同的方式理解（例如，“之间”对“直接之间”、“相邻”对“直接相邻”等）。
另外，在本发明中计算通信系统可靠性的步骤中，涉及到针对首中概率的求取。下面给出首中概率的基本概念和计算方法，以便对应后续内容的说明。
假设已知齐次有限马尔科夫链的一步转概阵P，状态i经过n(n>=2)步首次到状态j的概率为
Figure DEST_PATH_IMAGE001
Figure DEST_PATH_IMAGE002
称为状态i到j的n步首中概率。特别地，当i=j时，
Figure 264675DEST_PATH_IMAGE002
称为状态i的n步首返概率。
上式中，r为P的阶数，
Figure DEST_PATH_IMAGE003
为k到m的不经过j的n-2步转移概率。写成矩阵形式为
Figure DEST_PATH_IMAGE004
行矩阵R i 为P中状态i对应的行去掉p ij ，列矩阵C j 为P中状态j对应的列去掉p ij ，
Figure DEST_PATH_IMAGE005
为其他状态间不经过j的n-2步转移概率矩阵。
在一步转移概率矩阵P中，去掉状态j所对应的行和列，得到余子阵H j ，则
Figure DEST_PATH_IMAGE007
即
Figure DEST_PATH_IMAGE008
通过上式即可求取各种情形的首中（返）概率。
如图1所示，本发明实施方式中，基于通信链路多种状态计算IEC61850通信系统失效率方法包括步骤：（A）计算单条通信链路的最大允许连续丢包数量n；（B）定义单条通信链路的状态，计算通信链路可能存在的状态数量；（C）计算出通信系统故障状态空间的状态数量r，并规定通信系统各条链路都正常的状态记为G(1)，其余的状态分别记为G(2)~G(r) ，并将状态G(1)~G(r)一一对应于具体链路状态；（D）根据通信链路状态之间的转移概率为p j 计算通信系统状态一步转移概率矩阵P (1) ；（E）根据通信系统所有可能的状态序列中，首次出现系统传送有效（首中）的时刻特征对通信系统状态进行分类；（F）计算各状态序列类的概率（首中概率）；（G）计算在规定时间内通信系统的失效率。
对于步骤A，最大允许丢包数量指在规定的时间间隔内，除有且仅有一次传送成功外，其余传送均连续失败的最大传送次数。
根据通信链路传送的规定时间范围、采样间隔时间、一次传送成功所需要时间等参数，若通信链路i在规定时间t 1 内传送若干次数据，其中有一次传送成功即为该链路传送成功（链路有效），若通信链路相邻两次数据传送的时间间隔为t 2，一次传送成功所需要的时间为t 0 。则最大允许丢包数量的计算公式为
Figure DEST_PATH_IMAGE009
以图2所示的通信链路为例，t 1 取3ms，t 0 取1.5ms，t 2取1ms，则可计算得出
Figure DEST_PATH_IMAGE010
，即通信链路最大允许丢包数量为1。
对于步骤B，本专利所提出的关于信息链路状态的定义：
将事件“某固定时刻通信链路i将要发生连续j个时间间隔（包含当前时刻）的传送故障”记为状态
Figure DEST_PATH_IMAGE011
。
由于通信链路在规定时间内的最大传送次数为最大允许丢包数+1，考虑到每一时刻通信链路可能出现传送“正常”和“失败”两种情况，故通信链路的状态数为n+2。
图2所示的通信链路在第一次传送的对应时刻（当前时刻）可能存在的状态数量为2个，故其有可能在此时刻发生连续2个、1个传送失败、以及当前正常传送等3个状态，分别可记为
Figure DEST_PATH_IMAGE012
、
Figure DEST_PATH_IMAGE013
、
Figure DEST_PATH_IMAGE014
。
对于步骤C，考虑到单条通信链路的状态数为n+2，对含有m条链路的通信系统，其状态空间所包含的状态数量为
Figure DEST_PATH_IMAGE015
。分别对应状态空间G(e 1 ,e 2 ,…,e m )。其中m为通信系统包含的链路数，e 1 ,e 2 ,…,e m 对应通信系统各条链路的状态，e i =j表示通信链路i的状态为
Figure 273137DEST_PATH_IMAGE011
。
以图3为例，图中所示通信系统状态包含3条如图2所示的通信链路，每一链路均有3个状态，故图3所对应通信系统应有27个状态空间。分别对应于状态G(e 1 ,e 2 , e 3 )，其中e 1 , e 2 , e 3 均为0,1,2之中间的任一整数。
规定通信系统各条链路都正常的状态为G(1)，其余状态对应于G(2)~ G(r)。故图3中通信系统G(1)状态对应的状态空间为G(0,0,0)，而其余状态所对应状态空间中至少有一个链路的状态不为0。具体可进行一一对应如下表
图3所示通信系统状态与对应状态空间
通信系统状态 对应状态空间
<i>G</i>(1) <i>G</i>(0,0,0)
<i>G</i>(2) <i>G</i>(0,0,1)
<i>G</i>(3) <i>G</i>(0,0,2)
<i>G</i>(4) <i>G</i>(0,1,0)
<i>G</i>(5) <i>G</i>(0,2,0)
<i>G</i>(6) <i>G</i>(0,1,1)
<i>G</i>(7) <i>G</i>(0,1,2)
<i>G(8)</i> <i>G</i>(0,2,1)
<i>G</i>(9) <i>G</i>(0,2,2)
<i>G</i>(10) <i>G</i>(1,0,0)
<i>G</i>(11) <i>G</i>(1,0,1)
<i>G</i>(12) <i>G</i>(1,0,2)
<i>G</i>(13) <i>G</i>(1,1,0)
<i>G</i>(14) <i>G</i>(1,2,0)
<i>G</i>(15) <i>G</i>(1,1,1)
<i>G</i>(16) <i>G</i>(1,1,2)
<i>G</i>(17) <i>G</i>(1,2,1)
<i>G(18)</i> <i>G</i>(1,2,2)
<i>G</i>(19) <i>G</i>(2,0,0)
<i>G</i>(20) <i>G</i>(2,0,1)
<i>G</i>(21) <i>G</i>(2,0,2)
<i>G</i>(22) <i>G</i>(2,1,0)
<i>G</i>(23) <i>G</i>(2,2,0)
<i>G</i>(24) <i>G</i>(2,1,1)
<i>G</i>(25) <i>G</i>(2,1,2)
<i>G</i>(26) <i>G</i>(2,2,1)
<i>G</i>(27) <i>G</i>(2,2,2)
对于步骤D，通信链路状态之间的转移概率指某链路i经过一个时间间隔的转移概率，其值可通过仿真实验或者经验数据得到。据此得知，若某时刻链路i处于状态
Figure 730664DEST_PATH_IMAGE011
（1≤j ≤n+1），那么在经过一个时间间隔后，该链路的状态一定是L i (j-1)。故其转移概率一定为1，即p{L i (j)→L i (j-1)}=1，只有在
Figure 545036DEST_PATH_IMAGE014
时，其转移概率才为介于0~1之间的实数。
以图2为例，图中通信链路的状态转移概率具有如下特点：
p{L i (2)→L i (1)}=1， p{L i (1)→L i (0)}=1
p{L i (0)→L i (2)}=x p{L i (0)→L i (1)}=y， p{L i (0)→L i (0)}=z
其中x,y,z∈[0,1] 且x+y+z=1。其余转移概率均为0。
将通信系统由状态i转移到状态j的概率定义为p ij ，若通信系统状态G(i)的状态空间对应于G(a 1 ,a 2 ,…,a m )，通信系统状态G(j)的状态空间对应于G(b 1 ,b 2 ,…,b m ) 则
Figure DEST_PATH_IMAGE016
从而，通信系统的一步转移概率矩阵可表示为
Figure DEST_PATH_IMAGE017
对角元素在实际中只有在a 1 =a 2 =…=a m =0时才不为0，若其中
Figure DEST_PATH_IMAGE018
Figure DEST_PATH_IMAGE019
则对应的
Figure DEST_PATH_IMAGE020
，故矩阵中p ii =0（i=2,3,…,r）。
以图3为例，图中通信系统有27个状态，故其一步转移概率矩阵可表示为
Figure DEST_PATH_IMAGE021
为了更清楚的理解，矩阵中第一行和第一列的元素计算结果如下表所示
图3所示信息系统一步转移矩阵中第一行元素计算结果
矩阵元素 计算结果 矩阵元素 计算结果 矩阵元素 计算结果
<i>P</i><sub><i>1,1</i></sub> <i>z</i><sup><i>3</i></sup> <i>P</i><sub><i>1,10</i></sub> <i>yz</i><sup><i>2</i></sup> <i>P</i><sub><i>1,19</i></sub> <i>xz</i><sup><i>2</i></sup>
<i>P</i><sub><i>1,2</i></sub> <i>yz</i><sup><i>2</i></sup> <i>P</i><sub><i>1,11</i></sub> <i>y</i><sup><i>2</i></sup><i>z</i> <i>P</i><sub><i>1,20</i></sub> <i>xyz</i>
<i>P</i><sub><i>1,3</i></sub> <i>xz</i><sup><i>2</i></sup> <i>P</i><sub><i>1,12</i></sub> <i>xyz</i> <i>P</i><sub><i>1,21</i></sub> <i>x</i><sup><i>2</i></sup><i>z</i>
<i>P</i><sub><i>1,4</i></sub> <i>yz</i><sup><i>2</i></sup> <i>P</i><sub><i>1,13</i></sub> <i>y</i><sup><i>2</i></sup><i>z</i> <i>P</i><sub><i>1,22</i></sub> <i>xyz</i>
<i>P</i><sub><i>1,5</i></sub> <i>xz</i><sup><i>2</i></sup> <i>P</i><sub><i>1,14</i></sub> <i>xyz</i> <i>P</i><sub><i>1,23</i></sub> <i>x</i><sup><i>2</i></sup><i>z</i>
<i>P</i><sub><i>1,6</i></sub> <i>y</i><sup><i>2</i></sup><i>z</i> <i>P</i><sub><i>1,15</i></sub> <i>y</i><sup><i>3</i></sup> <i>P</i><sub><i>1,24</i></sub> <i>xy</i><sup><i>2</i></sup>
<i>P</i><sub><i>1,7</i></sub> <i>xyz</i> <i>P</i><sub><i>1,16</i></sub> <i>xy</i><sup><i>2</i></sup> <i>P</i><sub><i>1,25</i></sub> <i>x</i><sup><i>2</i></sup><i>y</i>
<i>P</i><sub><i>1,8</i></sub> <i>xyz</i> <i>P</i><sub><i>1,17</i></sub> <i>xy</i><sup><i>2</i></sup> <i>P</i><sub><i>1,26</i></sub> <i>x</i><sup><i>2</i></sup><i>y</i>
<i>P</i><sub><i>1,9</i></sub> <i>x</i><sup><i>2</i></sup><i>z</i> <i>P</i><sub><i>1,18</i></sub> <i>x</i><sup><i>2</i></sup><i>y</i> <i>P</i><sub><i>1,27</i></sub> <i>x</i><sup><i>3</i></sup>
图3所示信息系统一步转移矩阵中第一列元素计算结果
矩阵元素 计算结果 矩阵元素 计算结果 矩阵元素 计算结果
<i>P</i><sub><i>1,1</i></sub> <i>z</i><sup><i>3</i></sup> <i>P</i><sub><i>10,1</i></sub> <i>z</i><sup><i>2</i></sup> <i>P</i><sub><i>19,1</i></sub> <i>0</i>
<i>P</i><sub><i>2,1</i></sub> <i>z</i><sup><i>2</i></sup> <i>P</i><sub><i>11,1</i></sub> <i>Z</i> <i>P</i><sub><i>20,1</i></sub> <i>0</i>
<i>P</i><sub><i>3,1</i></sub> <i>0</i> <i>P</i><sub><i>12,1</i></sub> <i>0</i> <i>P</i><sub><i>21,1</i></sub> <i>0</i>
<i>P</i><sub><i>4,1</i></sub> <i>z</i><sup><i>2</i></sup> <i>P</i><sub><i>13,1</i></sub> <i>Z</i> <i>P</i><sub><i>22,1</i></sub> <i>0</i>
<i>P</i><sub><i>5,1</i></sub> <i>0</i> <i>P</i><sub><i>14,1</i></sub> <i>0</i> <i>P</i><sub><i>23,1</i></sub> <i>0</i>
<i>P</i><sub><i>6,1</i></sub> <i>z</i> <i>P</i><sub><i>15,1</i></sub> <i>1</i> <i>P</i><sub><i>24,1</i></sub> <i>0</i>
<i>P</i><sub><i>7,1</i></sub> <i>0</i> <i>P</i><sub><i>16,1</i></sub> <i>0</i> <i>P</i><sub><i>25,1</i></sub> <i>0</i>
<i>P</i><sub><i>8,1</i></sub> <i>0</i> <i>P</i><sub><i>17,1</i></sub> <i>0</i> <i>P</i><sub><i>26,1</i></sub> <i>0</i>
<i>P</i><sub><i>9,1</i></sub> <i>0</i> <i>P</i><sub><i>18,1</i></sub> <i>0</i> <i>P</i><sub><i>27,1</i></sub> <i>0</i>
按照所述的步骤E，按照0时刻后首次出现1状态的时刻分类，可得到n+1个状态序列类，具体如下表所示：
状态序列类的形成
Figure DEST_PATH_IMAGE022
以图3为例，图中通信系统有27个状态，可形成272=729个状态序列，在摒弃了系统无效的状态序列后，仅有27+26=52个有效状态序列，这些状态序列可分为两类，分别对应于首中在时刻1或时刻2。具体如下表
图3系统对应的状态序列类
Figure DEST_PATH_IMAGE023
注：表中“1”表示状态G(1)，“0”表示G(2)~G(27)中的任一状态，“*”表示G(1)~G (27)中的任一状态。
在所述的步骤F中，若状态序列k的概率表示为S k ，则
Figure DEST_PATH_IMAGE024
其中，
Figure DEST_PATH_IMAGE025
表示k步转移概率，指通信系统在时刻0以状态1为起点，经过k个时间间隔后在时刻k第一次变为状态1，对时刻k之后的各时刻状态没有要求。
以图3所示通信系统为例，
Figure DEST_PATH_IMAGE026
表示从0时刻起，从状态1出发经过1个时间间隔首次变为状态1的概率，即一步转移概率矩阵中的元素p 11 ，也对应一步首中概率。
Figure DEST_PATH_IMAGE027
表示从0时刻起，从状态1出发经过2个时间间隔首次变为状态1的概率，即状态1的两步首中概率。
在所述的步骤G中，通信系统在规定时间内的可用率为
Figure DEST_PATH_IMAGE028
规定时间内失效率为
Figure DEST_PATH_IMAGE029
以图3所示通信系统为例，
Figure DEST_PATH_IMAGE030
Figure DEST_PATH_IMAGE031
故图3通信系统在3ms内的可用率为
Figure DEST_PATH_IMAGE032
所以，图3通信系统在3ms内的失效率为
Figure DEST_PATH_IMAGE033
由此，通过采用本发明的基于通信链路多种状态计算IEC61850通信系统失效率的方法，针对通信链路的状态，计算通信链路的状态数量，符合通信链路实际数据采集和传送模式；依据通信系统数据传送时各不同时刻状态空间的变换所表现的时间特征，提出状态序列类的概念，并将其应用于通信系统概率计算，体现出了通信系统传送有效所具有的特性。另外，本发明基于通信系统的状态序列类的转移概率的计算方法，可以很方便的处理含多条链路、多个传送时刻的大型通信系统的传送概率计算问题，通过一步转移概率矩阵的较少次计算，实现了大型通信系统失效率计算。
需要说明的是，上述实施方式仅为本发明较佳的实施方案，不能将其理解为对本发明保护范围的限制，在未脱离本发明构思前提下，对本发明所做的任何微小变化与修饰均属于本发明的保护范围。

Claims (2)

1.基于通信链路多种状态计算IEC61850通信系统失效率的方法，所述方法步骤包括：
(A)计算单条通信链路的最大允许连续丢包数量n；
(B)定义单条通信链路的状态，计算通信链路可能存在的状态数量；其中，单条通信链路状态定义为某固定时刻通信链路i将要发生连续j个时间间隔的传送故障；通信链路的状态数量为最大允许丢包数+2；
(C)计算出通信系统故障状态空间的状态数量r，并规定通信系统各条链路都正常的状态记为G(1)，其余的状态分别记为G(2)～G(r)，并将状态G(1)G(r)一一对应于具体链路状态；
(D)根据通信链路状态之间的转移概率为pj计算通信系统状态一步转移概率矩阵P(1)；
(E)根据通信系统所有可能的状态序列中，首次出现系统传送有效的时刻特征对通信系统状态进行分类；所述步骤E中，对通信系统状态进行分类的具体方法：按照在0时刻后首次出现状态G(1)状态的顺序，对所有含状态G(1)的状态序列进行分类；不对0时刻后不包含G(1)的状态序列分类，其余所有状态序列分为n+1个状态序列类；其中，G(1)为通信系统各条链路都正常的状态；(F)按照首中概率计算方法，计算各状态序列类的概率；
计算在规定时间内通信系统的失效率。
2.根据权利要求1中所述的基于通信链路多种状态计算IEC61850通信系统失效率的方法，其特征在于，所述步骤E中，对通信系统状态序列的排列方式，其中对初始时刻的通信系统规定为正常，标记其状态为G(1)，然后提取发送数据对应的通信系统状态，将所提取的通信系统状态按照时间顺序排列形成“状态序列”，由于共计有n+1个发送时刻，故可形成若干由n+2个状态元素构成的状态序列，且每一个状态序列的第一个元素一定为G(1)。
CN201610840108.0A 2016-09-22 2016-09-22 基于通信链路多种状态计算iec61850通信系统失效率的方法 Expired - Fee Related CN107872328B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610840108.0A CN107872328B (zh) 2016-09-22 2016-09-22 基于通信链路多种状态计算iec61850通信系统失效率的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610840108.0A CN107872328B (zh) 2016-09-22 2016-09-22 基于通信链路多种状态计算iec61850通信系统失效率的方法

Publications (2)

Publication Number Publication Date
CN107872328A CN107872328A (zh) 2018-04-03
CN107872328B true CN107872328B (zh) 2021-04-02

Family

ID=61750937

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610840108.0A Expired - Fee Related CN107872328B (zh) 2016-09-22 2016-09-22 基于通信链路多种状态计算iec61850通信系统失效率的方法

Country Status (1)

Country Link
CN (1) CN107872328B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110177029B (zh) * 2019-05-30 2022-05-27 国网甘肃省电力公司经济技术研究院 一种电力通信网业务运行质量评估方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130053261A (ko) * 2011-11-15 2013-05-23 한양대학교 산학협력단 Ｉｅｃ 61850 규격에 따른 패킷 분석 장치
CN103323707A (zh) * 2013-06-05 2013-09-25 清华大学 基于半马尔科夫过程的变压器故障率预测方法
CN104156515A (zh) * 2014-07-28 2014-11-19 湖北文理学院 一种基于马尔可夫过程的ied形式化描述性能分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130053261A (ko) * 2011-11-15 2013-05-23 한양대학교 산학협력단 Ｉｅｃ 61850 규격에 따른 패킷 분석 장치
CN103323707A (zh) * 2013-06-05 2013-09-25 清华大学 基于半马尔科夫过程的变压器故障率预测方法
CN104156515A (zh) * 2014-07-28 2014-11-19 湖北文理学院 一种基于马尔可夫过程的ied形式化描述性能分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
分布式环境下软件可靠性模型的扩展与改进;韦群等;《第三届全国软件测试会议与移动计算、栅格、智能化高级论坛论文集》;20091231;第22卷;第83-87页 *

Also Published As

Publication number Publication date
CN107872328A (zh) 2018-04-03

Similar Documents

Publication Publication Date Title
DE102012102770B9 (de) System und Verfahren zur Fehlereingrenzung und Fehlerabschwächung basierend auf einer Netzwerkmodellierung
CN106056212B (zh) 一种人工神经网络计算核
CN107968840A (zh) 一种大规模电力设备监测报警数据实时处理方法及系统
CN105278460A (zh) 基于级联故障分析的数控机床系统组件可靠性评价方法
CN108848146A (zh) 一种基于时间触发通信业务的调度优化方法
CN101013955B (zh) 用于业务量矩阵估计的快速模拟退火
CN109408424A (zh) 一种基于PCIe接口的SpaceFibre总线数据采集方法
CN106779269A (zh) 基于变异系数法的智能变电站二次系统可靠性评估方法
CN107872328B (zh) 基于通信链路多种状态计算iec61850通信系统失效率的方法
CN116032845A (zh) 数据中心网络开销管理方法及系统
CN107454009B (zh) 面向数据中心的离线场景低带宽开销流量调度方案
CN102946443B (zh) 一种实现大规模数据传输的多任务调度方法
CN106130929A (zh) 基于图论算法的互联网保险领域的业务报文自动处理方法和系统
CN107832920A (zh) 海上风电场集群输电网可靠性评估方法和系统
CN106709593B (zh) 一种电力系统可靠性预测方法
CN102253978B (zh) 基于混合蚁群算法的公交车载多媒体文件更新方法
CN111210060B (zh) 一种工作日期间机房温度预测方法
Zhang et al. Transfer stability of urban subway network with passenger flow: Evidence in Beijing
CN105874488B (zh) 用于被动分布式定量集群知识的平台设备
Huang et al. Multiobjective multistate system preventive maintenance model with human reliability
WO2019141313A1 (de) Steuerungssystem, verwendung des steuerungssystems und steuerungsverfahren
Zhang et al. Design of CAN bus application layer protocol for aerocraft control system
CN116260128B (zh) 一种区域电网分析方法和装置
CN115062771B (zh) 一种分布式机器学习梯度汇聚方法、装置及模型训练方法
CN115660324B (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: 20210402

CF01 Termination of patent right due to non-payment of annual fee