一种卫星稳态可用度计算方法

摘要

一种卫星稳态可用度计算方法，根据卫星在轨安全模式的种类，并结合对每个安全模式的故障树分析结果和对底事件发生概率的统计，计算出卫星的进入每个安全模式的概率，进而得到卫星进入每个安全模式的相对发生概率。再统计在轨相似型号的进入每个安全模式的平均持续时间和两次进入安全模式的平均间隔时间。所有安全模式的相对发生概率乘以平均持续时间然后累加，得到卫星进入一次安全模式的平均中断时间，平均两次进安全模式的间隔时间减去卫星进入一次安全模式的平均中断时间再除以平均两次进安全模式的间隔时间就得到了卫星的稳态可用度。

说明书

技术领域

本发明涉及一种卫星稳态可用度计算方法，属于卫星总体设计技术领域。

背景技术

对于卫星来说，一般在轨不可修，属于不可修产品，而不可修产品一般强调可靠性而忽视可用性。近年来，随着对外出口型号的增多，对外交流时发现国外(特别是欧空局)比较重视卫星的可用性，比国内更加重视卫星的用户体验。而国内卫星目前还没有可用性指标和相应的计算方法，这在承担对外出口型号时就遇到了很多问题。

可用性是指产品在任一随机时刻需要开始使用时，处于工作状态或可使用状态的程度，用概率来度量这种程度就叫可用度。可用度一般有稳态可用度和瞬时可用度之分，瞬时可用度用来表示特定时刻产品处于工作或可使用状态的概率，而稳态可用度用来表示产品在一段时间内处于工作或可用状态的时间占总时间的比例，具体的计算方法有很多种。本发明提出了一种基于安全模式的卫星稳态可用度计算方法。

卫星在轨运行时，不可避免地会出现各种各样的问题。发生较严重的问题，而且卫星在轨又无法自主解决时，卫星就会进入安全模式，中断业务从而保证整星的安全。待地面故障排查，确定应对措施和遥控操作后，业务恢复才能正常。可以认为，卫星在轨进入安全模式即业务中断，从安全模式恢复即业务正常。可以通过计算卫星的平均每次进安全模式的中断时间占平均两次进安全模式的间隔时间的比例来计算卫星的在轨稳态可用度。

发明内容

本发明要解决的技术问题是：解决国内卫星没有可用性相关指标，没有适用的卫星可用性指标的分析计算方法。提供一种卫星稳态可用度计算方法，根据卫星在轨安全模式的种类，并结合对每个安全模式的故障树分析结果和对底事件发生概率的统计，计算出卫星的进入每个安全模式的概率，进而得到卫星进入每个安全模式的相对发生概率。再统计在轨相似型号的进入每个安全模式的平均持续时间和进入安全模式的平均间隔时间。所有安全模式的相对发生概率乘以平均持续时间然后累加，得到卫星进入一次安全模式的平均中断时间，卫星进入安全模式的平均间隔时间减去卫星进入一次安全模式的平均中断时间再除以进入安全模式的平均间隔时间就得到了卫星的稳态可用度。本方法针对目前国内卫星尚无可用性指标和可用性指标计算方法的情况，解决了对外出口型号用户提出的可用性要求无对应指标和无法分析验证的问题。

本发明目的通过以下技术方案予以实现：

一种卫星稳态可用度计算方法，包括如下步骤：

S1、以卫星在轨安全模式中每个安全模式作为顶事件建立故障树；

S2、根据故障树，统计各底事件的发生概率，计算出每个顶事件的发生概率，即每个安全模式的发生概率；

S3、根据所有安全模式的发生概率，计算每个安全模式的相对发生概率；

S4、基于现有数据中卫星在轨进入安全模式的时间和次数进行统计，获得卫星进入安全模式的平均间隔时间和每个安全模式的平均持续时间；

S5、根据每个安全模式的相对发生概率、每个安全模式的平均持续时间，计算卫星进入一次安全模式的平均中断时间；

S6、利用卫星进入安全模式的平均间隔时间、卫星进入一次安全模式的平均中断时间，确定卫星的稳态可用度。

上述卫星稳态可用度计算方法，优选的，S3中，计算每个安全模式的相对发生概率的方法为：

上述卫星稳态可用度计算方法，优选的，S4中，计算卫星进入安全模式的平均间隔时间的方法为：

上述卫星稳态可用度计算方法，优选的，S4中，计算卫星进入每个安全模式的平均持续时间的方法为：

上述卫星稳态可用度计算方法，优选的，S5中，计算卫星进入一次安全模式的平均中断时间的方法为：

MDT＝∑(MDT

上述卫星稳态可用度计算方法，优选的，S6中，卫星稳态可用度的计算方法为：

一种卫星稳态可用度计算装置，包括概率计算模块1、概率计算模块2、统计模块、平均中断时间计算模块、稳态可用度计算模块；

所述概率计算模块1根据故障树统计各底事件的发生概率，计算出每个顶事件的发生概率，即每个安全模式的发生概率；

所述概率计算模块2根据所有安全模式的发生概率，计算每个安全模式的相对发生概率；

所述统计模块统计现有数据中卫星在轨进入安全模式的时间和次数，获得卫星进入安全模式的平均间隔时间和每个安全模式的平均持续时间；

所述平均中断时间计算模块根据每个安全模式的相对发生概率、每个安全模式的平均持续时间，计算卫星进入一次安全模式的平均中断时间；

所述稳态可用度计算模块利用卫星进入安全模式的平均间隔时间、卫星进入一次安全模式的平均中断时间，确定卫星的稳态可用度。

上述卫星稳态可用度计算装置，优选的，计算卫星进入安全模式的平均间隔时间的方法为：

计算卫星进入每个安全模式的平均持续时间的方法为：

上述卫星稳态可用度计算装置，优选的，计算卫星进入一次安全模式的平均中断时间的方法为：

MDT＝∑(MDT

上述卫星稳态可用度计算装置，优选的，卫星稳态可用度的计算方法为：

本发明相比于现有技术具有如下有益效果：

(1)提出了使用稳态可用度作为卫星在轨运行时的可用性指标，解决了国内卫星目前还没有可用性指标的问题。

(2)通过卫星进出安全模式来统计卫星的在轨平均中断时间，较为方便统计数据和分析计算。

(3)通过计算卫星进入每个安全模式的相对发生概率来表征卫星在轨每次进入安全模式是特定安全模式的概率，从而可以计算卫星进入一次安全模式的平均中断时间。

(4)通过计算卫星在平均两次进安全模式之间的间隔时间内的稳态可用度来表征卫星整个在轨期间的稳态可用度，降低了分析计算的难度，使分析计算具有很好的实际可操作性。

(5)利用大量已经在轨运行的相似型号的统计数据来计算卫星平均进入一次安全模式的中断时间，进入安全模式的平均间隔时间，使分析计算结果更接近真实情况，更准确。

附图说明

图1为本发明的方法流程图。

图2为实施例第一个安全模式的故障树分析图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明的实施方式作进一步详细描述。

一种卫星稳态可用度计算方法，包括如下步骤：

S1、以卫星在轨安全模式中每个安全模式作为顶事件建立故障树；

S2、根据故障树，统计各底事件的发生概率，计算出每个顶事件的发生概率，即每个安全模式的发生概率；

S3、根据所有安全模式的发生概率，计算每个安全模式的相对发生概率；

S4、基于现有数据中卫星在轨进入安全模式的时间和次数进行统计，获得卫星进入安全模式的平均间隔时间和每个安全模式的平均持续时间；

S5、根据每个安全模式的相对发生概率、每个安全模式的平均持续时间，计算卫星进入一次安全模式的平均中断时间；

S6、利用卫星进入安全模式的平均间隔时间、卫星进入一次安全模式的平均中断时间，确定卫星的稳态可用度。

一种卫星稳态可用度计算装置，包括概率计算模块1、概率计算模块2、统计模块、平均中断时间计算模块、稳态可用度计算模块；

所述概率计算模块1根据故障树统计各底事件的发生概率，计算出每个顶事件的发生概率，即每个安全模式的发生概率；

所述概率计算模块2根据所有安全模式的发生概率，计算每个安全模式的相对发生概率；

所述统计模块统计现有数据中卫星在轨进入安全模式的时间和次数，获得卫星进入安全模式的平均间隔时间和每个安全模式的平均持续时间；

所述平均中断时间计算模块根据每个安全模式的相对发生概率、每个安全模式的平均持续时间，计算卫星进入一次安全模式的平均中断时间；

所述稳态可用度计算模块利用卫星进入安全模式的平均间隔时间、卫星进入一次安全模式的平均中断时间，确定卫星的稳态可用度。

作为本发明的一种优选方案，计算每个安全模式的相对发生概率的方法为：

作为本发明的一种优选方案，计算卫星进入安全模式的平均间隔时间的方法为：

作为本发明的一种优选方案，计算卫星进入每个安全模式的平均持续时间的方法为：

作为本发明的一种优选方案，计算卫星进入一次安全模式的平均中断时间的方法为：

MDT＝∑(MDT

作为本发明的一种优选方案，卫星稳态可用度的计算方法为：

实施例：

图1所示为本发明的方法流程图，本发明提供的一种卫星稳态可用度计算方法，具体步骤如下：

(1)在卫星在轨安全模式中依次选择一个安全模式作为顶事件进行故障树分析，底事件分析到可以统计其发生概率为止；

(2)根据故障树分析结果，统计各底事件的发生概率，计算顶事件的发生概率P

(3)依次计算所有安全模式的发生概率，按照下面的公式计算每个安全模式的相对发生概率P

(4)统计相似型号在轨进入安全模式的时间和次数，按照下面的公式计算卫星进入安全模式的平均间隔时间MTBSM：

式中T

按照下面的公式计算卫星进入每个安全模式的平均持续时间MDT

式中T

(5)根据每个安全模式的相对发生概率P

MDT＝∑(MDT

S6、按照下面的公式计算卫星的稳态可用度A：

具体的，以航天东方红卫星有限公司基于CAST-2000平台的某光学遥感卫星为例：

(1)在该卫星在轨安全模式中依次选择一个安全模式作为顶事件进行故障树分析，底事件分析到可以统计发生概率为止。该卫星的全部在轨安全模式如表1所示。

表1

以表1中第一个安全模式为例，故障树分析图如图2所示。

(2)根据故障树分析结果，统计各底事件的发生概率，计算顶事件的发生概率P

表2

根据图2的故障树逻辑关系，假设各底事件相互独立，各底事件之间是一种并联关系，则顶事件的发生概率P

P

(3)依次计算所有安全模式的发生概率。按照下面的公式计算每个安全模式的相对发生概率P

其他各安全模式的分析计算过程略，所有安全模式分析完成后，各安全模式的发生概率P

表3

(4)统计相似型号在轨进入安全模式的时间和次数，按照下面的公式计算卫星进入安全模式的平均间隔时间MTBSM：

根据统计，目前CAST-2000平台光学遥感卫星的MTBSM为75.147天，合1803.529小时。

按照下面的公式计算卫星进入每种安全模式的平均持续时间MDT

根据统计，得到CAST-2000平台光学遥感卫星在轨卫星进入每个安全模式的平均持续时间如表4所示。

表4

(5)根据每个安全模式的相对发生概率P

MDT＝∑(MDT

表5

(6)按照下面的公式计算卫星的稳态可用度A：

经过计算，该卫星的稳态可用度为：

A＝1-45.90799/1803.529＝0.974545≈97.45％。

本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。

本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。