CN108521346B - 一种基于终端数据的电信承载网异常节点定位方法 - Google Patents

Abstract

本发明公开了一种基于终端数据的电信承载网异常节点定位方法，包括：S1：获取终端设备的时序数据；S2：构建网络拓扑中每个节点的状态属性序列；S3：识别每个节点的状态属性序列中是否存在状态属性值大于或等于指标阈值，若存在，将状态属性值大于或等于指标阈值的节点作为对应时刻的基本警告节点；S4：筛出同一时刻的所有警示节点；S5：计算同一的每个警示节点的设备发生故障的概率和设备可靠性阈值；S6：分别判断同一的每个警示节点的设备发生故障的概率是否大于同一时刻的同一警示节点的设备可靠性阈值，若是，将对应警示节点的设备标记为真故障设备。本发明可以提高异常节点定位的准确性。

Description

一种基于终端数据的电信承载网异常节点定位方法

技术领域

本发明涉及通信技术领域，尤其涉及一种基于终端数据的电信承载网异常节点定位方法。

背景技术

与传统故障定位方法不同，传统电信承载网设备节点的故障定位方法具有结果误差较大，准确度不高，容易产生错误告警的特点，难以利用定位结果进行实际故障运维派单，这就使得传统故障定位方法难以满足电信运营商对承载网节点设备的管理和运维需求，使得电信运营商对电信运营商对承载网节点设备的管理和运维的效率严重下降。这些方法在进行异常定位处理时主要依赖于单台设备指标异常变化，而较少考虑该指标异常变化是否为受到网络拓扑中周边异常设备影响而变化。因此很容易忽略一个设备的异常特征是否是该设备的真正异常特征，也需要通过分析设备间异常干扰来过滤大部分冗余记录。故需要提供一种基于终端数据的电信承载网异常节点定位方法来提高定位结果准确率。

发明内容

本发明的目的是提供一种基于终端数据的电信承载网异常节点定位方法，其以电信承载网网络拓扑、终端数据作为计算数据，通过分析相关性分析对警告数据进行过滤，最后对警示节点进行设备可靠性分析，识别到真实警告设备，有效地减少了误告率，从而提高了整体定位的准确率。

本发明提供一种基于终端数据的电信承载网异常节点定位方法，包括如下步骤：

S1：获取电信承载网的网络拓扑以及所述网络拓扑中每个终端设备的时序数据；

其中，所述时序数据为终端设备在一个预设周期内各个时刻上异常代码的标记值；

S2：根据所述网络拓扑和终端设备的所述时序数据构建网络拓扑中每个节点的状态属性序列S；

其中，每个节点的状态属性序列S如下所示：

S＝[(t₁,V₁),…(t_i,V_i),……,(t_n,V_n)]

式中，t_i表示所述预设周期中的第i个时刻，V_i表示节点在t_i时刻的状态属性值，节点的状态属性值V_i的计算公式如下：

V_i＝V_i(D)＝f(v)或

其中，V_i(D)表示子节点为单个终端设备表示时的状态属性值，f(v)为终端设备的异常代码的标记值，V_i(F)表示节点为父级节点时的状态属性值，V(D_j)表示父级节点下的子节点的状态属性值，N表示父级节点下的子节点数量；

S3：识别每个节点的状态属性序列S中是否存在状态属性值大于或等于指标阈值，若存在，将状态属性值大于或等于指标阈值的节点作为对应时刻的基本警告节点W，并获取同一时刻的所有基本警告节点W；

S4：基于获取的同一时刻的所有基本警告节点W采用节点的相关性系数分析法依次筛出各个时刻的警示节点；

S5：基于失效率函数计算所述各个时刻的每个警示节点的设备发生故障的概率，以及基于累积失效概率函数计算同一时刻的同一警示节点的设备可靠性阈值；

S6：分别判断所述各个时刻的每个警示节点的设备发生故障的概率是否大于同一时刻的同一警示节点的设备可靠性阈值，若是，将大于同一时刻的同一警示节点的设备可靠性阈值的警示节点标记为异常节点且将对应设备标记为真故障设备，否则，不是真故障设备。

S2中单个终端设备必然是叶子节点，先获取到叶子节点的状态属性值，然后再根据网络拓扑的结构获取父级节点的状态属性值，进而从下往上计算得到每个时刻每个节点的状态属性值，故而构建出每个节点的状态属性序列S。网络拓扑中每个节点是单个终端设备或是由终端设备组成的设备组，此时也将其视为一个设备，一个节点即对应一个设备。

应当理解，S3中若节点的状态属性序列S中存在t_i时刻的状态属性值大于或等于指标阈值时，则将状态属性值大于或等于指标阈值的节点作为对应t_i时刻的基本警告节点W。其中，若节点的状态属性序列S存在多个t_i时刻的状态属性值大于或等于指标阈值时，则针对每个t_i时刻都会执行S4-S6。

进一步优选，将S4中获取的同一时刻的所有基本警告节点W记为同一t_i时刻的所有基本警告节点W，S4中每次筛出任意t_i时刻的警示节点的执行过程如下：

Step A：判断在t_i时刻根节点Q是否是基本警告节点W，若不是，执行Step C；若是，计算根节点Q与所述根节点Q的每个子节点P_j的相关性系数R(Q,P_j)并计算所述根节点Q的子节点相似比例P_WD1；

式中，

表示根节点Q的所有子节点P_j中相关性系数R(Q,P_j)中大于或等于相关系数阈值R_th的子节点P_j的个数，

表示根节点Q的所有子节点P_j的总个数；

Step B：判断所述根节点Q的子节点相似比例P_wD1是否大于或等于告警比例阈值P_wth，若是，所述根节点为t_i时刻的警示节点，若不是，执行Step C；

Step C：基于获取的t_i时刻的所有基本警告节点W得到根节点到每个基本警告节点W的节点链路，并针对每个节点链路分别采用节点相关性系数分析法识别出每个节点链路上t_i时刻的警示节点；

其中，每个基本警告节点W的节点链路是以根节点为链路起始节点，以基本警告节点W为链路终节点。

其中，进入Step C则表示根节点无故障。Step C中每个基本警告节点W对应一个节点链路。

进一步优选，Step C中对每个节点链路分别采用节点相关性系数分析法识别出每个节点链路上t_i时刻的警示节点的过程如下：

Step a：判断节点链路上是否存在子节点，若存在，按照从上至下的方向遍历节点链路获取第一个子节点E；若不存在，流程结束；

Step b：计算所述子节点E的父节点F与所述子节点E的相关性系数R(E,F)，以及判断所述相关性系数R(E,F)是否大于或等于相关系数阈值R_th；若是，执行Step g，否则，执行Step c；

Step c：计算所述子节点E的每个兄弟节点B_i与所述子节点E之间的相关性系数R(E,B_i)，并计算所述子节点E的兄弟相似比例P_WB；

式中，

表示子节点E的所有兄弟节点B_i中相关性系数R(E,B_i)中大于或等于相关系数阈值R_th的兄弟节点B_i的个数，

表示子节点E的兄弟节点B_i的总个数；

Step d：判断所述子节点E的兄弟相似比例P_WB是否小于告警比例阈值P_wth，若是，执行Step e；否则，执行Step g；

Step e：将所述子节点E标记为故障节点，并计算所述子节点E的每个子节点D_j与所述子节点E之间的相关性系数(E,D_j)，以及计算所述子节点E的子节点相似比例P_WD2；

式中，

表示子节点E的所有子节点D_j中相关性系数R(E,D_j)中大于或等于相关系数阈值R_th的子节点D_j的个数，

表示子节点E的子节点D_j的总个数；

Step f：判断所述子节点E的子节点相似比例P_WD2是否大于或等于告警比例阈值P_wth，若是，所述子节点E为t_i时刻的警示节点，流程结束；否则，执行Step l；

Step g：将所述子节点E定义为父节点，并返回Step a；

Step l：将所述子节点E定义为父节点，并将所述子节点E中相关性系数(E,D_j)大于或等于相关系数阈值R_th的子节点D_j作为基本告警节点W，再获取作为基本告警节点W的每个子节点D_j的节点链路，并针对每个子节点D_j的节点链路分别执行Step a-Step l，以识别出子节点D_j节点链路上t_i时刻的警示节点；

其中，所述子节点D_j的节点链路是以根节点为链路起始点，作为基本警告节点W的子节点D_j作为链路终节点。

Step b中若相关性系数R(E,F)大于或等于相关系数阈值R_th，即节点E的特征与父节点F的特征保持一致，没有故障，因此Step g中将当前的节点E定义为父节点F，然后进入下一层迭代，其过程可以理解为将当前的节点E定义为父节点F，将当前的节点E的子节点定义为下一层迭代中的字节点E。若相关性系数R(E,F)小于相关系数阈值R_th，即节点E的特征与父节点F的特征不一致，需要进一步判断。

Step d中若所述子节点E的兄弟相似比例P_WB小于告警比例阈值P_wth，则子节点E和兄弟节点B_i也不具备特征相似性，因此，故障来源子节点E本身，进而将子节点E标记为故障节点，对其进行进一步判断。否则P_WB≥P_wth，说明子节点E和其兄弟节点B_i有着较强相关性，则子节点E和父节点F也有着较强相似性，认为子节点E为非故障节点，进而Step g中当前的节点E定义为父节点F，然后进入下一层迭代。

Step f中若子节点E的子节点相似比例P_WD2大于或等于告警比例阈值P_wth，子节点E与各个子节点D_j故障特征表现一致，则进一步确认告警来源为子节点E；否则P_WD2＜P_wth，认为大部分子节点D_j和告警设备的特征表现不一致，这部分子节点D_j可以保持正常工作状态，因此故障源可能来自于那少部分子节点D_j，即对(E,D_j)大于R_th的部分子节点D_j分别作为基本告警节点W进行迭代，迭代过程实质上也是针对每个节点链路执行Step a-Step l。

其中，得到警示节点后，在同一t_i时刻，将警示节点下全部节点相关告警记录均清除掉，当作误告记录处理。

进一步优选，任意两个节点x与y的相关性系数R(x,y)的计算公式如下：

其中，Sx_i和Sy_i分别为节点x与节点y的状态属性序列S中t_i时刻对应的状态属性值V_i，n为状态属性序列S的长度。

需要说明的是，上述步骤中子节点E与父节点F相关性系数R(E,F)、子节点E的子节点D_j与子节点E之间的相关性系数(E,D_j)、子节点E的兄弟节点B_i与子节点E之间的相关性系数R(E,B_i)均是参照上述公式计算而来。本发明所选用的相关性系数R(x,y)为斯皮尔曼级相关系数。

进一步优选，所述告警比例阈值P_wth为0.7。

进一步优选，所述相关系数阈值R_th为0.8。

进一步优选，S5的执行过程如下：

首先，基于历史记录的每个警示节点的设备故障时刻获取t_i时刻的每个警示节点的设备寿命序列S_l；

S_l＝{l₁,…,l_m…,l_M}＝{t₁-t₀,…,t_m-t₀,…,t_M-t₀}，t_m＜t_i,m∈[1,M]

l_m＝t_m-t₀

式中，l_m表示警示节点的设备寿命序列S_l中t_m时刻对应的寿命，t_m表示历史记录中第m+1次发生故障的时间，t₀表示历史记录中第1次发生故障的时间，M为设备寿命序列S_l的长度；

然后，计算t_i时刻的每个警示节点的设备发生故障的概率；

l_i＝t_i-t₀

式中，P_i为t_i时刻的警示节点的设备发生故障的概率，λ(l_i)表示t_i时刻对应的失效率函数，l_i为警示节点在t_i时刻对应的寿命，α和β为失效率函数的参数，其中α和β的计算公式如下：

其中：

式中，

x_m、y_m均为中间参数：

再次，计算t_i时刻的每个警示节点的设备可靠性阈值；

P_thi＝F(l_i)+0.04，

式中，P_thi为t_i时刻的警示节点的设备可靠性阈值，F(l_i)表示t_i时刻对应的累积失效概率函数。

每个警示节点对应存在一个设备寿命序列S_l，其中，获取到t_i时刻的警示节点后，所选取的警示节点的历史记录是t_i时刻之前的历史数据。例如，状态属性序列S中预设周期为一天，以1小时为间距将预设周期分为24个时刻，t_i时刻即是一天中第i小时；设备寿命序列S_l的周期为7天，即t_i时刻前7天的数据，每一天也是对应24个时刻。

进一步优选，节点的指标阈值的计算公式如下：

V_th＝(1+0.06*U，

其中，V_th表示指标阈值，U为在连续7天的时间范围内节点的同一时刻中的所有状态属性值的平均值。

进一步优选，S1获取的每个终端设备的时序数据还包括每个终端设备的终端MOS值，终端设备的异常代码的标记值是根据终端设备的终端MOS值计算的；

其中，终端设备的异常代码的标记值f(v)的计算公式如下：

式中，v表示终端设备的终端MOS值，v_th为预设标准阈值。

进一步优选，所述预设标准阈值_th等于4.83。

有益效果

与现有预测技术相比，本发明的优点有：

本发明以电信承载网网络拓扑、终端数据作为计算数据，通过将终端数据量化汇集后，根据网络拓扑进行汇聚计算统计，形成时序数据；再对各个节点的时序数据进行判断识别出基本警告节点W，再基于节点的相关性系数分析即设备间的时序关联性对告警数据进行过滤进而获取到警示节点，最后通过对过滤后的警示节点对应进行设备可靠性分析，实现告警有效性判断，实现故障精确定位。

为了尽可能使得从告警设备中分辨出真告警设备以及找出真正故障设备。在进行虚假告警清洗和故障节点定位过程中，除使用网元设备本身监控数据和管理系统的告警数据外，还采用终端数据在网络拓扑节点上的分布状态特征来辅助分析。本发明上述过程通过分析故障带来的拓扑网络时序数据变化和分布状态变化，和其短时间内的周期性、关联性的特殊关系，以及通过较长时间设备可靠性因素有效地减少了误告率，从而提高了整体定位的准确率。能直观地通过终端分布状态来反映网络的健康状态，对于提高网络异常节点的定位的准确性具有一定的意义。本发明为电信运营商指导承载网运维、管理、科学定位故障节点和高效派单管理等需求，提供有价值的参考信息。

附图说明

图1是本发明提供一种基于终端数据的电信承载网异常节点定位方法的流程图；

图2是本发明实施例提供的图1的第一个子流程图；

图3是本发明实施例提供的图1的第三个子流程图；

图4是本发明实施例提供的图1的第二个子流程图；

图5是本发明实施例提供的节点链路示意图。

具体实施方式

下面将结合实施例对本发明做进一步的说明。

如图1所示，本发明提供一种基于终端数据的电信承载网异常节点定位方法，包括以下步骤：

S1：获取电信承载网的网络拓扑以及所述网络拓扑中每个终端设备的时序数据；

其中，所述时序数据包括为每个终端设备的终端MOS值(Mean Opinion Score)以及终端设备在一个预设周期内各个时刻上异常代码的标记值；其中，终端设备的异常代码的标记值是根据终端设备的终端MOS值计算的；

终端设备的异常代码的标记值f(v)的计算公式如下：

式中，v表示终端设备的终端MOS值，v_th为预设标准阈值。本实施例中优选预设标准阈值_th等于4.83。其中，f(v)＝0，表示异常代码不存在，为正常数据；f(v)＝1，表示出现异常代码，为异常数据；本实施例中以5分钟为时间间隔，若一个类别的异常代码多次出现，仅记录一次。

S2：根据所述网络拓扑和终端设备的所述时序数据构建网络拓扑中每个节点的状态属性序列S；

其中，每个节点的状态属性序列S如下所示：

S＝[(t₁,V₁),…(t_i,V_i),……,(t_n,V_n)]

式中，t_i表示所述预设周期中的第i个时刻，V_i表示节点在t_i时刻的状态属性值，节点的状态属性值V_i的计算公式如下：

V_i＝V_i(D)＝f(v)或

其中，V_i(D)表示叶子节点为单个终端设备表示时的状态属性值，f(v)为终端设备的异常代码的标记值，V_i(F)表示节点为父级节点时的状态属性值，V(D_j)表示父级节点下的子节点的状态属性值，N表示父级节点下的子节点数量。应当理解，网络拓扑中的节点可以是单个终端设备；还可以是多个终端设备组成的设备组，也将其视为一个设备，设备组必然是网络拓扑中的父级节点。

S3：识别每个节点的状态属性序列S中是否存在状态属性值大于或等于指标阈值，若存在，将状态属性值大于或等于指标阈值的节点作为对应时刻的基本警告节点W，并获取同一时刻的所有基本警告节点W；

例如：节点的状态属性序列S中存在t_i时刻的状态属性值大于或等于指标阈值，则将该节点作为t_i时刻的基本警告节点W，并获取同一t_i时刻的所有基本警告节点W。

如图2所示，即在t_i时刻，如果设备的状态属性值V_i满足V_i≥V_th，则将该节点标记为作为基本告警节点W。其中，节点的指标阈值V_th的计算公式如下：

V_th＝(1+0.06*U，

其中，U为在连续7天的时间范围内节点的同一时刻中的所有状态属性值的平均值。

S4：基于获取的同一时刻的所有基本警告节点W采用节点的相关性系数分析法依次筛出各个时刻的警示节点。

例如：基于获取的同一t_i时刻的所有基本警告节点W采用节点的相关性系数分析法筛出t_i时刻的警示节点。

本发明中节点的相关性系数分为节点与父节点的相关性系数、节点与兄弟节点的相关系数以及节点与子节点的相关性系数三类。其中，S4的执行过程下文将对此进行具体描述。

S5：基于失效率函数计算所述各个时刻的每个警示节点的设备发生故障的概率，以及基于累积失效概率函数计算同一时刻的同一警示节点的设备可靠性阈值。

例如：基于失效率函数计算t_i时刻的每个警示节点的设备发生故障的概率，以及基于累积失效概率函数计算同一t_i时刻的每个警示节点的设备可靠性阈值。

S5的执行过程如下：

首先，基于历史记录的每个警示节点的设备故障时刻获取t_i时刻的每个警示节点的设备寿命序列S_l；

S_l＝{l₁,…,l_m…,l_M}＝{t₁-t₀,…,t_m-t₀,…,t_M-t₀}，t_m＜t_i,m∈[1,M]

l_m＝t_m-t₀

式中，l_m表示警示节点的设备寿命序列S_l中t_m时刻对应的寿命，t_m表示历史记录中第m+1次发生故障的时间，t₀表示历史记录中第1次发生故障的时间，M为设备寿命序列S_l的长度。

然后，计算t_i时刻的每个警示节点的设备发生故障的概率；

l_i＝t_i-t₀

式中，P_i为t_i时刻的警示节点的设备发生故障的概率，λ(l_i)表示t_i时刻对应的失效率函数，l_i为警示节点在t_i时刻对应的寿命，α和β为失效率函数的参数，其中α和β的计算公式如下：

其中：

式中，

x_m、y_m均为中间参数；

再次，计算t_i时刻的每个警示节点的设备可靠性阈值；

P_thi＝F(l_i)+0.04，

式中，P_thi为t_i时刻的警示节点的设备可靠性阈值，F(l_i)表示t_i时刻对应的累积失效概率函数。

上述中累计失效概率函数和失效率函数是使用已知的最小二乘法估计的两参数韦布尔分布计算所得。

S6：分别判断所述各个时刻的每个警示节点的设备发生故障的概率是否大于同一时刻的同一警示节点的设备可靠性阈值，若是，将大于同一时刻的同一警示节点的设备可靠性阈值的警示节点标记为异常节点且将对应设备标记为真故障设备，否则，不是真故障设备。

例如：分别判断t_i时刻的每个警示节点的设备发生故障的概率是否大于同一t_i时刻的同一警示节点的设备可靠性阈值，若是，将t_i时刻大于设备可靠性阈值的警示节点的设备标记为真故障设备，否则，不是真故障设备。

最终可以得到异常节点列表。

具体的，如果在数据清洗过程中标记为告警节点的设备，在告警时间点t_i计算所得发生故障的概率满足：P_i＝P(l_i)＝λ(l_i)＞P_thi，则认为该设备会发生故障，且超出阈值范围越大，则认为故障可能性越大，并标记为真故障设备；否则不会。

相关性分析：

如图4和图5所示，将S4中获取的同一时刻的所有基本警告节点W记为同一t_i时刻的所有基本警告节点W，S4中每次筛出任意t_i时刻的警示节点的执行过程如下：

Step A：判断在t_i时刻根节点Q是否是基本警告节点W，若不是，执行Step C；若是，计算根节点Q与所述根节点Q的每个子节点P_j的相关性系数R(Q,P_j)并计算所述根节点Q的子节点相似比例P_WD1；

式中，

表示根节点Q的所有子节点P_j中相关性系数R(Q,P_j)中大于或等于相关系数阈值R_th的子节点P_j的个数，

表示根节点Q的所有子节点P_j的总个数。即根节点Q的子节点相似比例P_WD1等于

与

之比。

Step B：判断所述根节点Q的子节点相似比例P_WD1是否大于或等于告警比例阈值P_wth，若是，所述根节点为t_i时刻的警示节点，若不是，执行Step C；

Step C：基于获取的t_i时刻的所有基本警告节点W得到根节点到每个基本警告节点W的节点链路，并针对每个节点链路分别采用节点相关性系数分析法识别出每个节点链路上t_i时刻的警示节点；

其中，每个基本警告节点W的节点链路是以根节点为链路起始节点，以基本警告节点W为链路终节点。每个基本警告节点W对应一个节点链路。

其中，针对每个节点链路，Step C中对每个节点链路分别采用节点相关性系数分析法识别出每个节点链路上t_i时刻的警示节点的过程如下：

Step a：判断节点链路上是否存在子节点，若存在，按照从上至下的方向遍历节点链路获取第一个子节点E；若不存在，流程结束；

Step b：计算所述子节点E的父节点F与所述子节点E的相关性系数R(E,F)，以及判断所述相关性系数R(E,F)是否大于或等于相关系数阈值R_th；若是，执行Step g，否则，执行Step c；

其中，若相关性系数R(E,F)大于或等于相关系数阈值R_th，即节点E的特征与父节点F的特征保持一致，没有故障，因此Step g中将当前的节点E定义为父节点F，然后进入下一层迭代。若相关性系数R(E,F)小于相关系数阈值R_th，即节点E的特征与父节点F的特征不一致，需要进一步判断。

Step c：计算所述子节点E的每个兄弟节点B_i与所述子节点E之间的相关性系数R(E,B_i)，并计算所述子节点E的兄弟相似比例P_WB；

式中，

表示子节点E的所有兄弟节点B_i中相关性系数R(E,B_i)中大于或等于相关系数阈值R_th的兄弟节点B_i的个数，

表示子节点E的兄弟节点B_i的总个数；即子节点E的兄弟相似比例P_WB等于

与

之比。

Step d：判断所述子节点E的兄弟相似比例P_WB是否小于告警比例阈值P_wth，若是，执行Step e；否则，执行Step g；

其中，若所述子节点E的兄弟相似比例P_WB小于告警比例阈值P_wth，则子节点E和兄弟节点B_i也不具备特征相似性，因此，故障来源子节点E本身，进而将子节点E标记为故障节点，对其进行进一步判断。否则P_wB≥P_wth，说明子节点E和其兄弟节点B_i有着较强相关性，则子节点E和父节点F也有着较强相似性，认为子节点E为非故障节点，进而Step g中当前的节点E定义为父节点F，然后进入下一层迭代。

Step e：将所述子节点E标记为故障节点，并计算所述子节点E的每个子节点D_j与所述子节点E之间的相关性系数(E,D_j)，以及计算所述子节点E的子节点相似比例P_WD2；

式中，

表示子节点E的所有子节点D_j中相关性系数R(E,D_j)中大于或等于相关系数阈值R_th的子节点D_j的个数，

表示子节点E的子节点D_j的总个数；即子节点E的子节点相似比例P_WD2等于

与

之比。

Step f：判断所述子节点E的子节点相似比例P_WD2是否大于或等于告警比例阈值P_wth，若是，所述子节点E为t_i时刻的警示节点，流程结束；否则，执行Step l；

其中，若子节点E的子节点相似比例P_WD大于或等于告警比例阈值P_wth，子节点E与各个子节点D_j故障特征表现一致，则进一步确认告警来源为子节点E；否则P_WD2＜P_wth，认为大部分子节点D_j和告警设备的特征表现不一致，这部分子节点D_j可以保持正常工作状态，因此故障源可能来自于那少部分子节点D_j，即对(E,D_j)大于R_th的部分子节点D_j分别作为基本告警节点W进行迭代，迭代过程实质上也是针对每个节点链路执行Step a-Step l。

Step g：将所述子节点E定义为父节点，并返回Step a；

Step l：将所述子节点E定义为父节点，并将所述子节点E中相关性系数(E,D_j)大于或等于相关系数阈值R_th的子节点D_j作为基本告警节点W，再获取作为基本告警节点W的每个子节点D_j的节点链路，并针对每个子节点D_j的节点链路分别执行Step a-Step l，以识别出子节点D_j节点链路上t_i时刻的警示节点；

其中，所述子节点D_j的节点链路是以根节点为链路起始点，作为基本警告节点W的子节点D_j作为链路终节点。

其中，确认警示节点后，在相同t_i时刻，将警示节点下全部节点相关告警记录均清除掉，当作误告记录处理。

本实施例中，优选告警比例阈值P_wth为0.7，相关系数阈值R_th为0.8。

本发明中，任意两个节点x与y的相关性系数R(x,y)的计算公式如下：

其中，Sx_i和Sy_i分别为节点x与节点y的状态属性序列S中第i个时刻对应的状态属性值V_i，n为状态属性序列S的长度。

验证本发明的有效性

本发明在对判定模型的推测效果评估时，考虑到不仅需要检测出网络中可能的故障节点，同时也要尽可能避免误告。因此本发明以承载网7455个设备(OLT及其以上层级)的3个月历史监控数据和系统告警数据中随机抽取进行验证，并将所得结果与设备告警工单进行比较。为了更好的统计方法的准确率(Precision)，准确率定义为TP/(TP+FP)，其中TP是真阳性(True Positive)，FP是假阳性(False Positive)。TP定义为设备故障属于真故障设备，且被正确识别为故障设备，FP定义为设备故障属于非故障设备，但被识别为故障设备。

表1描述的是方法分别用于识别中国电信湖南省承载网全网2017年4月至6月份共3个月时间范围内告警的识别结果。

表1：本发明方法在承载网上的计算精度表

其中，表中部分指标含义：

1)准确率(Precision)：TP/(TP+FP)

2)召回率(Recall)：TP/(TP+FN)

3)F₁值：2*Precision*Recall/(Precision+Recall)

4)TN(True Negative)，真阴性：表示设备故障属于非故障设备，并识别为非故障设备；FN(False Negative)，假阴性：表示设备故障属于真故障设备，但被识别为非故障设备。

需要强调的是，本发明所述的实例是说明性的，而不是限定性的，因此本发明不限于具体实施方式中所述的实例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，不脱离本发明宗旨和范围的，不论是修改还是替换，同样属于本发明的保护范围。

Claims (10)

1.一种基于终端数据的电信承载网异常节点定位方法，其特征在于：包括如下步骤：

S1：获取电信承载网的网络拓扑以及所述网络拓扑中每个终端设备的时序数据；

其中，所述时序数据为终端设备在一个预设周期内各个时刻上异常代码的标记值；

S2：根据所述网络拓扑和终端设备的所述时序数据构建网络拓扑中每个节点的状态属性序列S；

其中，每个节点的状态属性序列S如下所示：

S＝[(t₁,V₁),…(t_i,V_i),……,(t_n,V_n)]

式中，t_i表示所述预设周期中的第i个时刻，V_i表示节点在t_i时刻的状态属性值，节点的状态属性值V_i的计算公式如下：

V_i＝V_i(D)＝f(v)或

其中，V_i(D)表示子节点为单个终端设备表示时的状态属性值，f(v)为终端设备的异常代码的标记值，V_i(F)表示节点为父级节点时的状态属性值，V(D_j)表示父级节点下的子节点的状态属性值，N表示父级节点下的子节点数量；

S3：识别每个节点的状态属性序列S中是否状态属性值大于或等于指标阈值，若存在，将状态属性值大于或等于指标阈值的节点作为对应时刻的基本警告节点W，并获取同一时刻的所有基本警告节点W；

S4：基于获取的同一时刻的所有基本警告节点W采用节点的相关性系数分析法依次筛出各个时刻的警示节点；

S5：基于失效率函数计算所述各个时刻的每个警示节点的设备发生故障的概率，以及基于累积失效概率函数计算同一时刻的同一警示节点的设备可靠性阈值；

S6：分别判断所述各个时刻的每个警示节点的设备发生故障的概率是否大于同一时刻的同一警示节点的设备可靠性阈值，若是，将大于同一时刻的同一警示节点的设备可靠性阈值的警示节点标记为异常节点且将对应设备标记为真故障设备，否则，不是真故障设备。

2.根据权利要求1所述的方法，其特征在于：将S4中获取的同一时刻的所有基本警告节点W记为同一t_i时刻的所有基本警告节点W，S4中每次筛出任意t_i时刻的警示节点的执行过程如下：

Step A：判断在t_i时刻根节点Q是否是基本警告节点W，若不是，执行Step C；若是，计算根节点Q与所述根节点Q的每个子节点P_j的相关性系数R(Q,P_j)并计算所述根节点Q的子节点相似比例P_WD1；

式中，

表示根节点Q的所有子节点P_j中相关性系数R(Q,P_j)中大于或等于相关系数阈值R_th的子节点P_j的个数，

表示根节点Q的所有子节点P_j的总个数；

Step B：判断所述根节点Q的子节点相似比例P_WD1是否大于或等于告警比例阈值P_wth，若是，所述根节点为t_i时刻的警示节点，若不是，执行Step C；

Step C：基于获取的t_i时刻的所有基本警告节点W得到根节点到每个基本警告节点W的节点链路，并针对每个节点链路分别采用节点相关性系数分析法识别出每个节点链路上t_i时刻的警示节点；

其中，每个基本警告节点W的节点链路是以根节点为链路起始节点，以基本警告节点W为链路终节点。

3.根据权利要求2所述的方法，其特征在于：Step C中对每个节点链路分别采用节点相关性系数分析法识别出每个节点链路上t_i时刻的警示节点的过程如下：

Step a：判断节点链路上是否存在子节点，若存在，按照从上至下的方向遍历节点链路获取第一个子节点E；若不存在，流程结束；

Step b：计算所述子节点E的父节点F与所述子节点E的相关性系数R(E,F)，以及判断所述相关性系数R(E,F)是否大于或等于相关系数阈值R_th；若是，执行Step g，否则，执行Stepc；

Step c：计算所述子节点E的每个兄弟节点B_i与所述子节点E之间的相关性系数R(E,B_i)，并计算所述子节点E的兄弟相似比例P_WB；

式中，

表示子节点E的所有兄弟节点B_i中相关性系数R(E,B_i)中大于或等于相关系数阈值R_th的兄弟节点B_i的个数，

表示子节点E的兄弟节点B_i的总个数；

Step d：判断所述子节点E的兄弟相似比例P_WB是否小于告警比例阈值P_wth，若是，执行Step e；否则，执行Step g；

Step e：计算所述子节点E的每个子节点D_j与所述子节点E之间的相关性系数R(E,D_j)，以及计算所述子节点E的子节点相似比例P_WD2；

式中，

表示子节点E的所有子节点D_j中相关性系数R(E,D_j)中大于或等于相关系数阈值R_th的子节点D_j的个数，

表示子节点E的子节点D_j的总个数；

Step f：判断所述子节点E的子节点相似比例P_WD2是否大于或等于告警比例阈值P_wth，若是，所述子节点E为t_i时刻的警示节点，流程结束；否则，执行Step l；

Step g：将所述子节点E定义为父节点，并返回Step a；

Step l：将所述子节点E定义为父节点，并将所述子节点E中相关性系数R(E,D_j)大于或等于相关系数阈值R_th的子节点D_j作为基本告警节点W，再获取作为基本告警节点W的每个子节点D_j的节点链路，并针对每个子节点D_j的节点链路分别执行Step a-Step l，以识别出子节点D_j节点链路上t_i时刻的警示节点；

其中，所述子节点D_j的节点链路是以根节点为链路起始点，作为基本警告节点W的子节点D_j作为链路终节点。

4.根据权利要求1-3任一项所述的方法，其特征在于：任意两个节点x与y的相关性系数R(x,y)的计算公式如下：

其中，Sx_i和Sy_i分别为节点x与节点y的状态属性序列S中t_i时刻对应的状态属性值V_i，n为状态属性序列S的长度。

5.根据权利要求3所述的方法，其特征在于：所述告警比例阈值P_wth为0.7。

6.根据权利要求3所述的方法，其特征在于：所述相关系数阈值R_th为0.8。

7.根据权利要求1所述的方法，其特征在于：S5的执行过程如下：

首先，基于历史记录的每个警示节点的设备故障时刻获取t_i时刻的每个警示节点的设备寿命序列S_l；

S_l＝{l₁,…,l_m…,l_M}＝{t₁-t₀,…,t_m-t₀,…,t_M-t₀}，t_m＜t_i,m∈[1,M]

l_m＝t_m-t₀

式中，l_m表示警示节点的设备寿命序列S_l中t_m时刻对应的寿命，t_m表示历史记录中第m+1次发生故障的时间，t₀表示历史记录中第1次发生故障的时间，M为设备寿命序列S_l的长度；

然后，计算t_i时刻的每个警示节点的设备发生故障的概率；

l_i＝t_i-t₀

式中，P_i为t_i时刻的警示节点的设备发生故障的概率，λ(l_i)表示t_i时刻对应的失效率函数，l_i为警示节点在t_i时刻对应的寿命，α和β为失效率函数的参数，其中α和β的计算公式如下：

其中：

式中，

x_m、y_m分别为中间参数；

再次，计算t_i时刻的每个警示节点的设备可靠性阈值；

P_thi＝F(l_i)+0.04，

式中，P_thi为t_i时刻的警示节点的设备可靠性阈值，F(l_i)表示t_i时刻对应的累积失效概率函数。

8.根据权利要求1所述的方法，其特征在于：节点的指标阈值的计算公式如下：

V_th＝(1+0.06)*U，

其中，V_th表示指标阈值，U为在连续7天的时间范围内节点的同一时刻中的所有状态属性值的平均值。

9.根据权利要求1所述的方法，其特征在于：S1获取的每个终端设备的时序数据还包括每个终端设备的终端MOS值，终端设备的异常代码的标记值是根据终端设备的终端MOS值计算的；

其中，终端设备的异常代码的标记值f(v)的计算公式如下：

式中，v表示终端设备的终端MOS值，v_th为预设标准阈值。

10.根据权利要求9所述的方法，其特征在于：所述预设标准阈值v_th等于4.83。

Images