CN107908597B - 一种基于复用分解的卫星有效载荷功耗分析方法 - Google Patents

Abstract

本发明公开了一种一种基于复用分解的卫星有效载荷功耗分析方法，包括步骤1：获取卫星执行任务时有效载荷的状态参数；步骤2：建立卫星平台维持常规载荷能耗模型；步骤3：计算卫星执行任务时有效载荷工作产生的能耗；步骤4：建立“单种类型任务”有效载荷能耗复用模型；步骤5：建立“多种类型任务”有效载荷能耗复用分解模型并求解；步骤6：根据卫星执行任务i时的实时遥测数据计算执行该任务时有效载荷总的实时功率。具有通过分析卫星遥测历史数据且不增加物理部件就可合理准确地计算卫星各有效载荷实时功率的优点。

Description

一种基于复用分解的卫星有效载荷功耗分析方法

技术领域

本发明涉及一种基于复用分解的卫星有效载荷功耗分析方法，属于卫星测控领域。

背景技术

卫星有效载荷是指直接执行特定卫星任务的仪器、设备或分系统。卫星由于其功能不同所搭载的有效载荷也各不相同，例如返回式卫星返回舱的有效载荷包括多光谱扫描仪、红外扫描仪、合成孔径雷达、微波辐射计、微波散射计、雷达高度计、超光谱成像仪以及遥感信息等设备。随着航天技术的不断发展，卫星任务规划及健康管理也逐步向精细化、智能化转变。而卫星由于工作环境的特殊性，其有限的能量储备成为卫星任务规划与健康管理的关键制约因素。目前我国对卫星能量智能管理的研究还处在起步阶段，对卫星执行任务所需能耗缺乏准确计算方法导致难以对卫星能源系统状态实现有效预测。所以，目前卫星在执行任务时所采取的策略是给定较高储存能量下限，以牺牲卫星的有效工作量来确保卫星运行的安全性。这种策略大大缩减了卫星可执行任务的规划范围，是对卫星资源的一种浪费。例如，在我国某资源卫星的在轨使用原则中规定的有效载荷在轨使用单圈工作时长、工作频率等大大低于实际可承受范围。这是从卫星在轨能源平衡的角度为确保卫星安全运行给出的一个比较严格的条件。如果能准确掌握卫星有效载荷工作与能耗的关系，对卫星执行任务的能量需求有准确的把握，那么就能实时把握卫星能量系统状态，进而更合理地提出任务规划的约束，提高卫星资源的利用率。

实现卫星能源系统的智能管理，关键一步在于对卫星执行任务时有效载荷功率的准确计算。卫星由于制造工艺及制造成本的限制，无法对每一个部件都安装足够的传感器获得准确数据，因此不能直接对有效载荷的功率进行计算。另一方面，卫星在执行某一任务时，往往通过多个有效载荷协调运行共同完成任务，大部分的卫星有效载荷都存在复用的情况。如何从复用的有效载荷中得到每一个具体载荷的功率是一个关键问题。

目前国内外针对卫星能源智能管理和能耗计算的研究较少，几乎没有现成的方法对其进行求解。由于卫星上的各部件主要依靠电力运行，因此本方法研究的卫星有效载荷功耗针对的是电功耗。目前针对电能计算的研究主要可分为两类，第一类是从理论的角度研究电能的计算公式。计算电能的理论公式主要从电路的电压、电流、电阻考虑，公式主要有Q＝UIT、Q＝I²RT、Q＝U²T/R以及三相交流电中的

虽然利用电能公式计算的方法是最简单最直接的，但是需要知道单个有效载荷准确的电功率或者电流，目前大多数卫星无法满足这要求。关于电能计算另一类研究是从工程实践的角度出发，研究各种修正方法减少由于器件测量不准确等因素造成功耗计算误差，使得电能计算结果更精确。有研究基于相位差校正法的虚拟谐波功率与电能计算方法实现数字化变电站中的虚拟谐波功率与电能计量。也有研究从电能计量装置的误差出发，考虑电能表的误差、互感器的合成误差等对电能计量造成的误差，并对引起误差的原因进行深入分析。有研究从电网电量计算受谐波污染的问题出发，针对间谐波污染严重的电网，为准确计量相应负荷实际消耗的电能量，减小电能计量误差，在传统电能计量方法基础上，增加计及间谐波产生的电能损耗，并推导出存在间谐波条件下的电能计算公式，最后通过数值仿真分析，证明了在间谐波污染严重的电网中，采用提出的计及间谐波的电能计量方法，将更加合理、准确、可靠。上述研究都致力于解决单一载荷的功耗问题，但并不适用于复用的载荷功耗计算。

发明内容

针对现有技术存在的问题，本发明要解决的问题在于提出了一种通过分析卫星遥测历史数据且不增加物理部件就可合理准确地计算卫星各有效载荷实时功率的基于复用分解的卫星有效载荷功耗分析方法。

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

一种基于复用分解的卫星有效载荷功耗分析方法，包括以下步骤：

步骤1：获取卫星执行任务时有效载荷的状态参数；

步骤2：建立卫星平台维持常规载荷能耗模型；

步骤3：计算卫星执行任务时有效载荷工作产生的能耗；

步骤4：建立“单种类型任务”有效载荷能耗复用模型；

步骤5：建立“多种类型任务”有效载荷能耗复用分解模型并求解；

步骤6：根据卫星执行任务i时的实时遥测数据计算执行该任务时有效载荷总的实时功率。

作为本发明的进一步改进：

进一步地，所述步骤2的具体方法为：

步骤2.1：根据充电阵输出电流识别卫星运转周期；

步骤2.2：获取卫星未执行任务时的历史遥测数据，利用RK算法挖掘卫星一个周期内的变化模式，利用多项式拟合方法建立卫星维持常规载荷的母线电压和母线电流的时间函数f(t)和g(t)；

步骤2.3：计算卫星维持常规载荷的实时能耗为p(t)＝f(t)×g(t)；

步骤2.4:计算卫星在一段时间内维持常规载荷的能耗为

T_start为卫星执行任务的开始时间，T_stop为卫星执行任务的结束时间。

进一步地，所述步骤3的具体方法为：

步骤3.1：根据卫星上行指令，获取卫星执行每个任务的开始时间T_start和结束时间T_stop以及任务类型M_type，并据此获取卫星执行任务时的开始时间T_start和结束时间T_stop之间母线电压v(t)及母线电流i(t)的遥测数据，计算卫星执行该任务卫星的总能耗W_total；

步骤3.2：计算卫星平台在任务执行期间用于维持常规载荷的能耗W_base

步骤3.3：得到由于有效载荷工作引起的能耗W_payload为：

W_payload＝W_total-W_base (3)。

进一步地，步骤4中所述“单种类型任务”有效载荷能耗复用模型为：

其中，i为M种任务类型中的第i类任务，i∈{1,...,M}，N_i为参与第i类任务的有效载荷数，

为第j∈{1,...,N_i}个载荷的工作开始时间；

为第j个载荷的工作结束时间；

为第j个有效载荷在第i类任务内的平均功率。

进一步地，步骤5中所述“多种类型任务”有效载荷能耗复用分解模型为

所述求解结果为：卫星有效载荷的平均功率

进一步地，所述步骤6的具体方法为：

步骤6.1：将历史遥测数据包含的n个任务分解出多个任务组，每个任务组里包含卫星所执行的所有M种任务类型，对每一个任务组，利用步骤3至步骤5求解出各种任务类型的有效载荷的平均功率，并保存求解结果；

步骤6.2：根据步骤6.1算出的每个任务组下各有效载荷的平均功率，计算所有任务组中相同有效载荷平均功率的平均值，作为相同有效载荷在工作时的总平均功率，将执行任务的各有效载荷的总平均功率比作为执行该任务的有效载荷的实时功率比；

步骤6.3：计算卫星执行任务i时的有效载荷总的实时功率；

P_payload(t)＝v(t)×i(t)-f(t)×g(t) (6)

步骤6.4：根据执行任务i时有效载荷的实时功率比，计算任务i中有效载荷的实时功率。

其中：a:b:c为执行任务i时所参与的三个有效载荷A,B,C的实时功率比，P_A(t)，P_B(t)，P_C(t)为三个有效载荷A,B,C的实时功率。

进一步地，所述步骤6.1中将历史遥测数据包含的n个任务分解出多个任务组的具体方法是指从第一个任务开始，直至第n个任务，如果从第x个任务至第y个任务之间包含所有M种任务类型，则将其中不含重复任务类型的任务视为一个任务组，然后从第y+1个任务开始往后寻找新的任务组，以此类推。

与现有技术相比，本发明的有益效果在于：

本发明一种基于复用分解的卫星有效载荷功耗分析方法，通过分析已有的历史遥测数据，可以不用增加物理部件，也不用增加成本就计算出卫星有效载荷的平均功耗，并通过各有效载荷总平均功率比估算执行某一任务时各有效载荷的实时功率比，进而可以合理估算出执行某一任务时有效载荷的实时功率。

附图说明

图1本发明的系统流程图。

图2为利用RK算法挖掘卫星运行一个周期内长期负载情况下母线电流和母线电压的变化模式。

具体实施方式：

下面结合具体实施例来做进一步的说明。本发明所研究的卫星属于资源探测卫星，选择2015年6月至7月的历史遥测数据进行研究，期间共执行任务114次，在进行项目分析之前，需要根据有关卫星在轨使用描述文件，该卫星任务执行模式与载荷的关系如表1所示，本实施例中卫星有效载荷包括三线阵相机、多光谱相机、激光测距仪(搭载)、数据传输分系统、数传天线分系统以及数据记录分系统。其中由于激光测距仪属于搭载的载荷，且时刻都在工作，所以不将其作为有效载荷考虑。因此，本实施例中研究的任务类型M＝5。

表1任务模式包含的有效载荷

图1和图2示出了本发明一种基于复用分解的卫星有效载荷功耗分析方法，包括以下步骤：步骤1：获取卫星执行任务时有效载荷的状态参数；步骤2：建立卫星平台维持常规载荷能耗模型；步骤3：计算卫星执行任务时有效载荷工作产生的能耗；步骤4：建立“单种类型任务”有效载荷能耗复用模型；步骤5：建立“多种类型任务”有效载荷能耗复用分解模型并求解；步骤6：根据卫星执行任务i时的实时遥测数据计算执行该任务时有效载荷总的实时功率。本发明通过分析已有的历史遥测数据，计算出卫星各有效载荷的平均功率，以及计算出各个有效载荷在工作时总平均功率，由于卫星在执行某一种任务时，其执行过程固定，因此该模式下包含的有效载荷功率也不变，所以可以根据有效载荷总平均功率比估算执行该任务时各有效载荷的实时功率比，进而可以根据实时遥测数据计算出卫星执行某一任务时有效载荷总的平均功率以及功率比，不用增加物理部件，也不用增加成本就可合理估算出执行某一任务时有效载荷的实时功率。

本实施例中，步骤1：获取卫星执行任务时有效载荷的状态参数；各有效载荷的状态参数如表2所示。

表2有效载荷状态监测参数

步骤2：建立卫星平台维持常规载荷能耗模型；

步骤2.1：根据充电阵输出电流识别卫星运转周期；得到卫星运转一圈的平均周期为5676秒。

步骤2.2：获取卫星未执行任务时的历史遥测数据，利用RK算法挖掘卫星一个周期内的变化模式，利用多项式拟合方法建立卫星维持常规载荷的母线电压和母线电流的时间函数f(t)和g(t)；本实施例中根据卫星上行任务列表识别长期负载时间段，以卫星电源系统母线电压与母线电流为研究对象，利用RK算法挖掘卫星一个周期内的变化模式，其结果如附图2所示。然后利用函数拟合方法建立母线电压和母线电流的时间函数f(t)和g(t)，得到母线电压与母线电流变化的拟合结果为：

其中t∈[0,5676]为一周期内的时间，

步骤2.3：计算卫星维持常规载荷的实时能耗为p(t)＝f(t)×g(t)；

计算卫星平台维持常规载荷的实时能耗为：

p(t)＝f(t)×g(t)＝314.94+29×sin(t×2π/306)瓦

步骤2.4:计算卫星在一段时间内维持常规载荷的能耗为

T_start为卫星执行任务的开始时间，T_stop为卫星执行任务的结束时间。

步骤3：计算卫星执行某一类型任务时有效载荷工作产生的总能耗；

步骤3.1：根据卫星上行指令，获取卫星执行每个任务的开始时间T_start和结束时间T_stop以及执行任务的类型M_type，并据此获取卫星执行任务时的开始时间T_start和结束时间T_stop之间母线电压v(t)及母线电流i(t)的遥测数据；本实施例中以表1中的图像测绘实传模式为例，从卫星历史上行指令列表中找到每种模式的开始时间T_start和结束时间T_stop分别为以2016年6月1日10:33至以2016年6月1日10:41，根据公式(1)计算卫星执行该任务卫星的总能耗W_total，

步骤3.2：计算卫星平台在任务执行期间用于维持常规载荷的能耗W_base。根据步骤2，

计算出卫星在任务执行期间用于维持常规载荷的能耗，

步骤3.3：得到由于有效载荷工作引起的能耗W_payload为：

W_payload＝W_total-W_base＝141117.2(W)。

对前100个任务的有效载荷消耗总能耗计算结果如表3所示。

表3前100次任务有效载荷的总能耗及总功率

步骤4：建立“单种类型任务”有效载荷能耗复用模型；以任务类型1图像测绘实传模式为例，图像测绘实传模式在执行时参与工作的有效载荷有三线阵相机、多光谱相机、数传分系统以及数传天线。以2016年6月1日10:33至以2016年6月1日10:41之间执行的图像测绘实传模式为研究对象，选择三线阵相机成像开/关状态(TMY051)、多光谱相机成像开/关状态(TMY061)、数据处理器a(TMR149)、数传天线加电状态(TMD290)为各个有效载荷开关指示参数。得到在执行该模式期间各载荷的执行时间分别如表4所示。可知在执行该模式时，有效载荷是完全复用的。

表4各有效载荷实际执行时间

有效载荷	执行时段	执行时长
			三线阵相机	10:33—10:41	480
多光谱相机	10:33—10:41	480
			数传分系统	10:33—10:41	480
数传天线	10:33—10:41	480

假设

分别表示三线阵相机、三线阵相机(仅正视相机)、多光谱相机、数传分系统、数传天线以及数据记录分系统在执行测绘实传任务时的平均功率。根据公式(4)

得到该图像测绘实传模式执行的有效载荷能耗复用模型为：

步骤5：建立“多种类型任务”有效载荷能耗复用分解模型并求解；

从历史任务执行列表中选择所有不同种类的任务各一例，建立“多种类型任务”有效载荷能耗复用分解模型，其分别属于图像测绘实传模式、图像资源实传模式、图像测绘记录模式、图像资源记录模式、回放模式，各任务执行时间分别为2016/6/1 10:33:50--2016/6/1 10:41:50；2016/6/1 13:42:50--2016/6/1 13:50:50；2016/6/1 16:58:50--2016/6/118:28:45；2016/6/1 18:23:45--2016/6/1 18:28:45；2016/6/1 21:10:55--2016/6/1 21:16:05；组成一个任务组。分别建立各个任务的“单种类型任务”有效载荷能耗复用模型，联立得到复用分解的线性方程组：

求解上述线性方程组时发现该线性方程组无解，原因在于有效的模式类型只有5类，而求解变量有6个，这里需要对载荷进行合并。通过对载荷工作时间的比对，可知三线阵相机与多光谱相机是同时开始工作的，所以将三线阵相机、三线阵相机(仅正视相机)、多光谱相机组合成“三线阵相机+多光谱相机组”，“三线阵相机(仅正视相机)+多光谱相机组”两个变量，并分别用

和

表示。则线性方程组可以转化为：

所求解出的有效载荷平均功率计算结果如表5所示：

表5有效载荷平均功率计算结果

步骤6：根据卫星执行任务i时的实时遥测数据计算执行该任务时有效载荷总实时功率。

步骤6.1：将历史遥测数据包含的n个任务分解出多个任务组，每个任务组里包含卫星所执行的所有M种任务类型，对每一个任务组，利用步骤3至步骤5求解出各种任务类型的有效载荷的平均功率，并保存求解结果。

步骤6.2：根据步骤6.1算出的历史多个任务组下各有效载荷的平均功率，计算所有任务组中相同有效载荷平均功率的平均值，作为相同有效载荷在工作时的总平均功率，将执行任务的各有效载荷的总平均功率比作为执行该任务的有效载荷的实时功率比；

本实例选择的2016年6月1日以至2016年6月15日之间包含114个任务，由于建立复用分解线性模型需要五种模式的能耗情况，因此首先要选出五种不同的任务组成任务组，根据任务列表，在前100个任务中挑选出6组任务。任务组包含的任务如表6所示，根据步骤3-5计算得到6组任务下各个有效载荷的平均功率如表7所示。

表6挑选出的任务组及其包含的任务代号

表7各个任务组中有效载荷平均功率计算结果

各有效载荷总平均功率会随着新的任务组数据的加入而不断更新。

计算执行各种类型任务时参与工作的有效载荷的功率比；得到6个任务组的各载荷平均功率后，求每种载荷平均功率的平均值即得有效载荷总平均功率值。以“三线阵相机+多光谱相机”总平均功率为例，其计算公式为：

得到总平均功率后，根据总平均功率比可以估算在执行某一任务时各有效载荷的实时功率比。以执行测绘实传任务为例，参与任务的有效载荷包括“三线阵相机+多光谱相机”、“数据传输分系统”以及“数传天线分系统”。则在执行测绘实传任务时三种载荷的总平均功率比计算如下：

用相同方法可以得到在不同任务下各有效载荷的总平均功率比，作为执行该任务的有效载荷的实时功率比。

本实施例中，各种任务类型下有效载荷功率比如表8所示。

表8各任务类型下各有效载荷总平均功率比

步骤6.3：计算卫星执行任务i时的有效载荷总实时功率；

P_payload(t)＝v(t)×i(t)-f(t)×g(t)

则：P_payload(t)＝v(t)×i(t)-314.94+29×sin(t×2π/306)

步骤6.4：根据执行任务i时有效载荷的实时功率比，计算任务i中有效载荷的实时功率。举例如下：

其中：a:b:c为执行任务i时所参与的三个有效载荷A、B、C的实时功率比，P_A(t)，P_B(t)，P_C(t)为三个有效载荷A、B、C的实时功率。

在提出解决方法之前，首先要介绍目前卫星运行时的能量系统变化特点、卫星执行任务的一般原则及有效载荷工作的相关特性进行说明，作为提出本方法的基础。

1)卫星各任务独立执行

卫星执行的任务是通过不同的工作模式来完成的，例如一般探测卫星执行的任务包含图像测绘实传模式、图像资源实传模式、图像测绘记录模式、图像资源记录模式、边记边传模式、回放模式等。任务的执行是按照序列顺序执行。卫星在执行任务时，同一时间内最多只能执行一种工作模式，且执行一种模式必定是从模式开始至模式结束完整、连续的执行过程。

2)卫星平台能耗平稳

卫星在未执行任务时，卫星上的能耗用于维持卫星平台的正常运行。卫星平台的能耗比较平稳且随着卫星所处光照区和阴影区呈周期性变化。

3)相同任务参与执行的有效载荷及工作模式一致

一般而言，卫星有效载荷由于执行任务的不同可以选择不同的工作模式，而在执行相同类型的任务时，有效载荷的工作模式是相同的。且有效载荷在相同的工作模式下其功率是一定的。

4)有效载荷的工作状态可由相关的参数识别出来

目前已有的遥测参数中包含有效载荷工作状态的参数，可以通过状态参数来识别有效载荷的运行状态及时间。例如对于资源卫星上的三线阵相机，标明其工作状态的参数就有三线阵相机成像开/关状态、前后视相机成像加断电状态、正视相机成像加断电状态、三线阵相机A总线状态等，都可以标识三线阵相机的工作状态及工作模式。

在此基础之上，提出一种基于复用分解的卫星有效载荷功耗分析方法解决有效载荷功耗计算与实时功率分析问题。

以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。

Claims (6)

1.一种基于复用分解的卫星有效载荷功耗分析方法，其特征在于：包括以下步骤：

步骤1：获取卫星执行任务时有效载荷的状态参数；

步骤2：建立卫星平台维持常规载荷能耗模型；

步骤3：计算卫星执行任务时有效载荷工作产生的能耗；

步骤4：建立“单种类型任务”有效载荷能耗复用模型；

步骤5：建立“多种类型任务”有效载荷能耗复用分解模型并求解；

步骤6：根据卫星执行任务i时的实时遥测数据计算执行该任务时有效载荷总的实时功率；

所述步骤2的具体方法为：

步骤2.1：根据卫星太阳能帆板的输出电流识别卫星运转周期；

步骤2.2：获取卫星未执行任务时的历史遥测数据，利用RK算法挖掘卫星一个周期内的变化模式，利用多项式拟合方法建立卫星维持常规载荷的母线电压和母线电流的时间函数f(t)和g(t)；

步骤2.3：计算卫星维持常规载荷的实时能耗p(t)＝f(t)×g(t)；

步骤2.4:计算卫星在一段时间内维持常规载荷的能耗为

其中：T_start为卫星执行任务的开始时间，T_stop为卫星执行任务的结束时间。

2.根据权利要求1所述的一种基于复用分解的卫星有效载荷功耗分析方法，其特征在于：所述步骤3的具体方法为：

步骤3.1：根据卫星上行指令，获取卫星执行每个任务的开始时间T_start和结束时间T_stop以及任务类型M_type，并据此获取卫星执行任务时的开始时间T_start和结束时间T_stop之间母线电压v(t)及母线电流i(t)的遥测数据，计算卫星执行该任务的总能耗W_total

步骤3.2：计算卫星平台在任务执行期间用于维持常规载荷的能耗W_base

步骤3.3：得到由于有效载荷工作引起的能耗W_payload为：

W_payload＝W_total-W_base (3)。

3.根据权利要求2所述的一种基于复用分解的卫星有效载荷功耗分析方法，其特征在于：步骤4中所述“单种类型任务”有效载荷能耗复用模型为

i为M种任务类型中的第i类任务，i∈{1,...,M}，N_i为参与第i类任务的有效载荷数，

为第j∈{1,...,N_i}个有效载荷的工作开始时间；

为第j个有效载荷的工作结束时间；

为第j个有效载荷在第i类任务内的平均功率。

4.根据权利要求3所述的一种基于复用分解的卫星有效载荷功耗分析方法，其特征在于：步骤5中所述“多种类型任务”有效载荷能耗复用分解模型为

所述求解结果为：卫星有效载荷的平均功率

5.根据权利要求4所述的一种基于复用分解的卫星有效载荷功耗分析方法，其特征在于：所述步骤6的具体方法为：

步骤6.1：将历史遥测数据包含的n个任务分解成多个任务组，每个任务组里包含卫星所执行的所有M种任务类型，对每一个任务组，利用步骤3至步骤5求解出各种任务类型的有效载荷的平均功率，并保存求解结果；

步骤6.2：根据步骤6.1算出的每个任务组下各有效载荷的平均功率，计算所有任务组中相同有效载荷平均功率的平均值，作为相同有效载荷在工作时的总平均功率，将执行任务的各有效载荷的总平均功率比作为执行该任务的有效载荷的实时功率比；

步骤6.3：计算卫星执行任务i时的有效载荷的总实时功率；

P_payload(t)＝v(t)×i(t)-f(t)×g(t) (6)

步骤6.4：根据执行任务i时有效载荷的实时功率比，计算任务i中有效载荷的实时功率。

6.根据权利要求5所述的一种基于复用分解的卫星有效载荷功耗分析方法，其特征在于：所述步骤6.1中将历史遥测数据包含的n个任务分解出多个任务组的具体方法是指从第一个任务开始，直至第n个任务，如果从第x个任务至第y个任务之间包含所有M种任务类型，则将其中不含重复任务类型的任务视为一个任务组，然后从第y+1个任务开始往后寻找新的任务组，以此类推。

Images