基于马尔可夫理论的能量管理系统可靠性分析方法与流程

本发明涉及软件可靠性评估技术领域，具体涉及基于马尔可夫理论的能量管理系统可靠性分析方法。

背景技术：

目前能量管理系统的可靠性关乎电力系统的安全运行。传统的能量管理系统可靠性分析方法大致分为两种：其一为根据软件开发过程中测试结果，对软件进行修改和可靠性评估；其二为根据应用的相关关系，视整个系统的应用为串并联方式复合使用的软件。然而，测试的方式计算可靠性在需要极高可靠性和实际运行的电力系统显然不适用；而视为串并联的式复合使用的方式分析显现不出应用具体的内部关系，也不能很好结合软件实际应用频率等基本情况，更不能反映各个应用在系统中的重要程度。

技术实现要素：

为了克服现在技术的缺陷，基于离散马尔可夫理论，对能量管理系统的应用结构进行建模，结合应用实际使用频率确定概率转移概率矩阵计算管理系统的可靠性。

为解决上述问题，本发明提供以下技术方案：

一种基于马尔可夫理论的能量管理系统可靠性分析方法，包括以下步骤：

1)基于离散马尔可夫理论，根据多个应用的调用关系，对能量管理系统的应用结构进行建模；

2)结合应用实际使用频率确定转移概率矩阵；

3)计算单个应用的可靠性；

4)根据应用的转移概率和可靠性计算出系统的可靠性。

进一步的，所述步骤1)建模方法为：把能量管理系统的若干个应用等效为马尔可夫模型中的若干个节点；再根据每个应用的实际调用关系表示应用的调用转移关系。

进一步的，所述步骤3)单个应用可靠性计算公式为：

mtbf为平均无故障时间,mttr为平均故障恢复时间；

当一个应用中有n个服务，若采用串联方式，应用可靠性若为两服务并联，则：s＝1-(1-s1)(1-s2)。

进一步的，所述步骤2的使用频率的确定方式为：根据实际开发的能量管理系统的应用所设定的刷新频率，计算出应用在单位时间内使用的次数。

进一步的，架构中每个节点都是马尔可夫过程的一种状态，一个节点存在一个或多个分支，有向分支(ni，nj)表示ni到nj的可能转移，单一构件ni的控制转移概率之和为1，即：

所述节点转移概率计算方法为：

其中：pij表示能量管理系统中应用nj调用ni的概率，kij为应用构件nj在单位时间内调用ni的次数。

进一步的，根据转移概率和应用的可靠性，可得到随机转移矩阵：

包含了系统的入口和出口的两个吸收态，去除吸收态后得到应用间转移概率矩阵为：

系统可靠性可表示为：

rsys＝s(1,n)rn＝(i-q)^-1(1,n)rn。

本发明还提供一种基于马尔可夫理论的能量管理系统可靠性分析方法用于评估应用敏感性：

根据实际运行经验可知，不同应用在系统中的使用情况和重要性不同，定义应用ni的敏感度si为：

令y＝i-q，可得

式中|y|是y的行列式，值为1，|x|是y(n，1)对应的余子式；求取构件ni的敏感度si，则对|x|的第i行进行展开：

|x|＝x(i,1)αi1+x(i,2)αi2+…+x(i,n)αin,i＝1,…,n-1

式中αij是x(i，j)对应的代数余子式，αij和pij与ri不相关，因此可得

x＝βiri

其中βi和ri不相关，系统结构和其他已知构件的可靠性相关，具体值由pij和α决定；因此系统构件敏感度si可表示为：

本发明具有以下有益效果：

1.本发明根据电力能量管理系统的实际应用使用频率确定转移概率矩阵，可靠性结果更符合系统的实际运行情况；根据可靠性计算结果进行应用的敏感度分析，分析结果能够为制定提升系统可靠性方案提供参考，对进行日常系统维护工作具有指导意义。在计算出系统可靠性后，对系统应用进行敏感度分析，计算出各个应用的敏感度，反映各个应用在系统中的重要程度，为系统开发制定应用可靠性指标提供参考。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。

图1本发明基于马尔可夫理论的能量管理系统可靠性分析方法流程图

图2为能量管理系统中应用调用控制图；

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例1

请参见附图1，一种基于马尔可夫理论的能量管理系统可靠性分析方法，包括以下步骤：

1)基于离散马尔可夫理论，根据多个应用的调用关系，对能量管理系统的应用结构进行建模。本发明通过将电力能量管理系统的应用等效为马尔可夫模型中的一个节点n，模型中用圆圈表示，如附图2所示：n1-n14分别表示能量管理系统中的实时数据库、历史数据库、实时报警、断路器状态、主接线、状态估计、变压器状态、实时线损监测、负荷预测、电能质量检测、无功优化、潮流计算、故障选线、调度主界面应用；把用节点代表的系统应用的调用关系用有向线段连接起来，表示应用的调用转移关系，如n1分别被n2-n6调用，n2分别被n7-n10调用。

2)结合应用实际使用频率确定转移概率矩阵；根据实际开发的能量管理系统的应用所设定的刷新频率，计算出应用在单位时间内使用的次数。

架构中每个节点都是马尔可夫过程的一种状态，一个节点存在一个或多个分支，有向分支(ni，nj)表示ni到nj的可能转移，单一构件ni的控制转移概率之和为1，即：

其中，pij表示能量管理系统中应用nj调用ni的概率。

结合步骤2)系统应用的使用频率，计算各个节点的转移概率，计算公式如下：

其中：kij为应用构件nj在单位时间内调用ni的次数。

3)计算单个应用的可靠性；

单个应用可靠性计算公式为：

mtbf为平均无故障时间,mttr为平均故障恢复时间；

当一个应用中有n个服务，若采用串联方式，应用可靠性计算公式如下：

若为两服务并联，应用可靠性计算公式如下：

s＝1-(1-s1)(1-s2)。

4)根据应用的转移概率和可靠性计算出系统的可靠性；

根据转移概率和应用的可靠性，可得到随机转移矩阵：

上述转移矩阵包含了系统的入口和出口的两个吸收态，去除吸收态后得到应用间转移概率矩阵为：

系统可靠性可表示为：

rsys＝s(1,n)rn＝(i-q)^-1(1,n)rn

实施例2

本发明还提供一种基于马尔可夫理论的能量管理系统可靠性分析方法用于评估应用敏感性：

根据实际运行经验可知，不同应用在系统中的使用情况和重要性不同，定义应用ni的敏感度si为：

需要寻找ri和si的函数关系即可得到si，令y＝i-q，可得

式中|y|是y的行列式，值为1，|x|是y(n，1)对应的余子式；求取构件ni的敏感度si，则对|x|的第i行进行展开：

|x|＝x(i,1)αi1+x(i,2)αi2+…+x(i,n)αin,i＝1,…,n-1

式中αij是x(i，j)对应的代数余子式，αij和pij与ri不相关，因此可得

x＝βiri

其中βi和ri不相关，系统结构和其他已知构件的可靠性相关，具体值由pij和α决定；因此系统构件敏感度si表示为：

通过计算出系统各个应用的敏感性，可以为在实际开发中优先保证搞灵敏度应用的可靠性以保证系统的可靠性。

对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。