CN111703593B - 地球同步轨道卫星表面充电效应的预警方法及预警系统 - Google Patents

地球同步轨道卫星表面充电效应的预警方法及预警系统 Download PDF

Info

Publication number
CN111703593B
CN111703593B CN202010600627.6A CN202010600627A CN111703593B CN 111703593 B CN111703593 B CN 111703593B CN 202010600627 A CN202010600627 A CN 202010600627A CN 111703593 B CN111703593 B CN 111703593B
Authority
CN
China
Prior art keywords
surface charging
satellite
index
time
early warning
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.)
Active
Application number
CN202010600627.6A
Other languages
English (en)
Other versions
CN111703593A (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.)
Beijing Ruanqi Technology Co ltd
National Satellite Meteorological Center
Original Assignee
Beijing Ruanqi Technology Co ltd
National Satellite Meteorological Center
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 Beijing Ruanqi Technology Co ltd, National Satellite Meteorological Center filed Critical Beijing Ruanqi Technology Co ltd
Priority to CN202010600627.6A priority Critical patent/CN111703593B/zh
Publication of CN111703593A publication Critical patent/CN111703593A/zh
Application granted granted Critical
Publication of CN111703593B publication Critical patent/CN111703593B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00—Cosmonautic vehicles
    • B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/52—Protection, safety or emergency devices; Survival aids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B23/00—Alarms responsive to unspecified undesired or abnormal conditions

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Critical Care (AREA)
  • Emergency Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Emergency Management (AREA)
  • Remote Sensing (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Charge And Discharge Circuits For Batteries Or The Like (AREA)

Abstract

本发明公开了一种地球同步轨道卫星表面充电效应的预警方法及预警系统，包括：获取T0时刻的AE指数数据；推算AE指数时刻所对应卫星的地方时；根据AE指数数据，计算所述地球同步轨道卫星的表面充电电位模式计算值U；根据所述表面充电电位模式计算值U，获得表面充电指数ISC；根据所述表面充电指数ISC，对卫星表面充电危险情况进行预警。本发明所述的方法可以有效地对可能的表面充电效应提供可靠数据支撑，对于卫星异常的分析和处置工作也有参考作用。

Description

地球同步轨道卫星表面充电效应的预警方法及预警系统
技术领域
本发明涉及卫星表面充电技术领域，具体涉及一种地球同步轨道卫星表面充电效应的预警方法及预警系统。
背景技术
当太阳耀斑发生时，太阳风与地球磁场相互作用会产生地磁亚暴现象，此时在地球同步轨道会产生能量在高能大通量等离子体，会使运行在地球同步轨道的卫星发生严重的表面带电现象，特别是没有光照的位置，孤立介质的充电电压可能超过10000伏；完全处在地影中的卫星整体表面也会与周围环境形成很高的电位，高电位所形成的放电过程可能干扰卫星的正常运行，严重时还会造成结构的永久损坏。因此，获得地球磁场与卫星表面充电之间的关系，根据二者之间的关系对卫星表面充电进行危险等级预警显得尤为重要。
综上，有必要提供一种获取卫星表面充电情况并对卫星表面充电进行危险预警的方法。
发明内容
为了解决背景技术中的至少一个技术问题，本发明提供一种地球同步轨道卫星表面充电效应的预警方法及预警系统，该方法通过AE指数计算得到卫星的表面充电电位，根据卫星的表面充电电位获得表面充电指数ISC，进而对卫星表面充电危险情况进行预警，能够在指导卫星采取防护措施和故障快速诊断方面提供可靠信息。
本发明所采用的技术方案具体如下：
根据本发明的一个方面，提供一种地球同步轨道卫星表面充电效应的预警方法，包括以下步骤：
获取T0时刻的AE指数数据；
推算AE指数时刻所对应卫星的地方时；
根据AE指数数据，计算卫星的表面充电电位模式计算值U；
根据表面充电电位模式计算值U，获得表面充电指数ISC；
根据表面充电指数ISC，对卫星表面充电危险情况进行预警。
可选地，获取AE指数数据，还包括：根据获取的AE指数数据曲线，计算距离T0时刻最近的1小时内AE指数的平均值。
可选地，推算AE指数时刻所对应的卫星的地方时还包括以下步骤：
获取卫星定点的经度L；
根据如下公式获得AE指数时刻T0所对应的卫星的地方时TL：TL＝T0+L/15；
式中：卫星定点的经度L的单位为度。
可选地，根据AE指数数据，计算地球同步轨道卫星的表面充电电位模式计算值U，包括：
将卫星所在区域的地方时划分为四个时域；
根据不同时域AE指数与表面充电电位的不同的对应关系，获得表面充电电位模式计算值U。
可选地，当地方时处于6:00到21:00时域时，表面充电电位模式计算值U为0。
可选地，当地方时处于21:00到24:00时域时，表面充电电位模式计算值U为：
Figure GDA0003116192820000021
式中：AE为距离T0时刻最近的1小时内AE指数的平均值，单位为nT；
U为表面充电电位模式计算值，单位为V。
可选地，当地方时处于0:00到4:00时域时，表面充电电位模式计算值U为：
Figure GDA0003116192820000022
式中：AE为距离T0时刻最近的1小时内AE指数的平均值，单位为nT；
U为表面充电电位模式计算值，单位为V。
可选地，当地方时处于4:00到6:00时域时，表面充电电位模式计算值U为：
Figure GDA0003116192820000023
式中：AE为距离T0时刻最近的1小时内AE指数的平均值，单位为nT；
U为表面充电电位模式计算值，单位为V。
可选地，根据表面充电电位模式计算值U，获得表面充电指数ISC，还包括：
当U≥-5000时，ISC＝1；
当-10000≤U<-5000时，ISC＝2；
当U≤-10000，ISC＞2。
可选地，根据表面充电指数ISC，对卫星表面充电危险情况进行预警，还包括：根据表面充电指数ISC值判定卫星表面充电危险等级。
可选地，根据表面充电指数ISC值判定卫星表面充电危险等级，还包括：当表面充电指数ISC＝1，则判断卫星表面充电危险等级为低级。
可选地，根据表面充电指数ISC值判定卫星表面充电危险等级，还包括：当表面充电指数ISC＝2，则判断卫星表面充电危险等级为中级。
可选地，根据表面充电指数ISC值判定卫星表面充电危险等级，还包括：当表面充电指数ISC＞2，则判断卫星表面充电危险等级为高级。
可选地，还包括：当卫星表面充电危险等级为高级时，发出表面充电警报。
根据本发明的另一个方面，还提供一种地球同步轨道卫星表面充电效应的预警系统，包括：
数据获取单元，用于获取T0时刻的AE指数数据；
第一数据处理单元，用于推算AE指数时刻所对应的卫星的地方时；
第二数据处理单元，用于根据AE指数数据，计算卫星的表面充电电位模式计算值U；
第三数据处理单元，用于根据表面充电电位模式计算值U，获得表面充电指数ISC；
预警单元，用于根据表面充电指数ISC，对卫星表面充电危险情况进行预警。
可选地，第一数据处理单元还设置有：
数据获取模块，用于获取卫星定点的经度L；
第一数据计算模块，用于将AE指数时刻T0所对应的卫星的地方时TL。
可选地，第二数据处理单元设置有：
时域划分模块，用于将卫星所在区域的地方时划分为四个时域；
第二数据计算模块，用于根据不同时域AE指数与表面充电电位的不同的对应关系，获得表面充电电位模式计算值U。
可选地，预警单元还包括报警模块，用于根据表面充电指数ISC值判定卫星表面充电危险等级，并且卫星表面充电危险等级为高级时，进行报警。
与现有技术相比，本发明所述的地球同步轨道卫星表面充电效应的预警方法及预警系统至少具备如下有益效果：
本发明所述的地球同步轨道卫星表面充电效应的预警方法及预警系统，首先获取AE指数数据；根据所述AE指数与表面充电电位模式计算值U的关系，计算所述表面充电电位模式计算值U；根据所述表面充电电位模式计算值U与表面充电指数ISC的对应关系，获得表面充电指数ISC；根据所述表面充电指数ISC，对卫星表面充电危险情况进行预警；其中根据AE指数获取卫星表面的充电电位计算值U，可以有效地对可能的表面充电效应提供可靠数据支撑，并且，根据卫星表面的充电电位模式的计算值U及表面充电指数ISC的对应关系，获得表面充电指数ISC，可以对卫星表面的充电效应情况进行预警，确定卫星表面的充电效应等级，快速对卫星表面的充电危险情况作出反应，进而能够在指导卫星采取防护措施和故障快速诊断方面提供可靠信息，对于卫星异常的分析和处置工作有参考作用。
附图说明
图1为本发明一实施例提供的地球同步轨道卫星表面充电效应的预警方法流程图；
图2为本发明一实施例中获取的2020年4月20日的AE指数时间分布图；
图3为本发明一实施例中记录的2020年4月20日风云四号卫星表面的充电电压变化曲线；
图4为本发明另一实施例中获取的2020年3月20日的AE指数时间分布图；
图5为本发明另一实施例中记录的2020年3月20日风云四号卫星表面充电电压变化曲线；
图6为本发明实施例表3中表面充电电位模式计算值与风云四号实测值的对比图；
图7为本发明另一实施例中地球同步轨道卫星表面充电效应的预警系统示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
地球同步轨道是各种应用卫星集中的区域，也是地磁亚暴期间高温等离子体注入显著增加的地区，高温等离子体的快速增加是导致卫星表面充电发生的主要因素。而引发等离子体注入的地磁亚暴的重要标志，就是地磁AE指数的显著增长。
AE指数是空间科学所使用的描述地磁亚暴活动的一个参数。它表示极光电集流所造成的总的磁扰动，单位为nT。该指数的生成是把沿极光带上经度均匀分布的台站观测到的磁场北向分量(以静日观测值为基线)的时间变化画在一张图上，并根据图得到曲线簇的包络线。上包络线给出AU指数，它表示东向电集流引起的最大的磁场扰动；下包络线给出AL指数，它表示西向电集流引起的最大磁场扰动；上下包络线之间的距离给出AE指数，它表示总的扰动强度,称为极光电集流指数。
AE指数与地磁亚暴期间等离子体注入相关，因而是反映地球附近等离子体环境的合适参数。另外，由于该参数有地基台站数据生成，而且相关台站的实时性很好，所以AE指数具有实时性强、连续，可靠性高的特征，适合卫星表面充电这种业务性工作。
通过分析风云四号表面充电观测数据发现：表面充电事件均发生在磁尾时段，即22:00至次日6:00时间区间；通过对比风云四号搭载表面充电探测器和等离子体探测器的数据得出的：表面充电事件对应等离子体温度、密度大幅增加；表面充电事件发生于AE较高的时段，但在不同的地方时，二者幅度的对应关系不同；等离子体显著充电通常发生在AE高的时段及其后1小时区间。基于上述风云四号表面充电探测结果分析，卫星的表面充电事件均发生在卫星处于磁尾区域的时间段；表面充电事件的发生与AE指数的增长密切相关，而且显著的AE指数增加通常对应着严重的充电事件。
基于此，本实施例提供一种地球同步轨道卫星表面充电效应的预警方法，如图1所示，包括以下步骤：
S101：获取T0时刻的AE指数数据；
本实施例中的AE指数数据是由世界数据中心获取，网址为：http://wdc.kugi.kyoto-u.ac.jp/ae_realtime/index.html。
通过上述网址，选择T0时刻的AE指数数据，此处T0时刻对应的时间为国际标准时UT。由于网站AE指数文件所提供的数据为快视图，在本实施例中，需要根据获取的AE指数的数据曲线粗略得出T0时刻最近1小时AE指数的平均值。
S102：推算AE指数时刻所对应的卫星的地方时；
首先，获取卫星定点的经度L；
AE指数T0对应的时间为国际标准时UT，若卫星定点的经度为L(单位为：度)，则该AE指数时刻对应的卫星的地方时TL可表示为：
TL＝T0+L/15 (4)
式中：卫星定点的经度L的取值根据东经和西经进行划分，其中东经为正，西经为负。卫星定点的经度L的单位为度。
S103：根据所述AE指数，计算所述卫星的表面充电电位模式计算值U；
由于不同的地方时，地磁AE指数的增大对应的局地等离子体的增幅不同，本发明实施例中将卫星所在区域的地方时划分为4个时域，具体为：6:00到21:00为第一时域，21:00到24:00为第二时域，0:00到4:00为第三时域及4:00到6:00为第四时域。各时域AE指数与表面充电电位的对应关系不同，AE指数与表面充电电位模式计算值U的关系如下：
1)地方时6:00到21:00为安全期，在此期间不会发生表面充电效应，即表面充电电位模式计算值U始终为0；
2)地方时21:00到24:00为过渡区-I。取AE时间点为卫星地方时TL，则TL-1至TL+1时段(时间单位为h)，表面充电电位模式计算值U为：
Figure GDA0003116192820000061
式中：AE值为AE指数1小时平均值，单位为nT；
U为表面充电电位模式计算值，单位为V。
3)地方时0:00到4:00为扰动区。取AE时间点为卫星地方时TL，则TL-1至TL+1时段(时间单位为h)，表面充电电位模式计算值U为：
Figure GDA0003116192820000062
式中AE及U的意义及单位同上。
4)地方时4:00到6:00为过渡区-II。取AE时间点为卫星地方时TL，则TL-1至TL+1时段(时间单位为h)，表面充电电位模式计算值U为：
Figure GDA0003116192820000063
式中AE及U的意义及单位同上。
根据公式(1)、公式(2)和公式(3)，给出基于AE指数的表面充电电位模式计算值U。
S104：根据所述表面充电电位模式计算值U与表面充电指数ISC的对应关系，获得表面充电指数ISC；
为了能够适应地球同步轨道上其他卫星的应用，本实施例中将基于AE指数的表面充电电位模式计算值U用表面充电危险程度用表面充电指数ISC表示。
其中，表面充电指数ISC与计算得到的表面充电电压模式计算值U的关系为：
若U≥-5000，ISC＝1。
若-10000≤U<-5000，ISC＝2。
若U≤-10000，ISC＞2。
式中：U的单位为V。
根据U值与ISC的对应关系，获得表面充电指数ISC。
S105：根据所述表面充电指数ISC，对卫星表面充电危险情况进行预警。
表面充电指数ISC值与卫星表面充电危险等级对应情况如下：
ISC＝1表示卫星表面充电危险等级为低级；
ISC＝2表示卫星表面充电危险等级为中级；
ISC＞2表示卫星表面充电危险等级为高级。
根据上述表面充电指数ISC值判定卫星表面充电危险等级，对卫星表面充电危险进行预警，当所述卫星表面充电危险等级为高级时，发出表面充电警报。
下面以风云四号卫星为例，对本发明所述地球同步轨道卫星表面充电效应的预警方法及预警系统进行详细的说明。
在一可选实施例中，以2020年4月20日风云四号的表面充电探测数据与当天的AE指数的相关特性为例。
在世界数据中心获取2020年4月20日的AE指数数据，下载2020年4月20日的探测数据(当日的12:00UT前后发生地磁暴，AE指数增长明显)，可得到该日的AE指数连续快视图，如图3所示，其中图中上部的包络线表示AE指数的大小。
将AE指数对应的国际标准时T0转化为卫星地方时TL。风云四号定点于东经104度，由公式(4)计算得到其地方时与国际标准时(UT)的对应关系，风云四号卫星的地方时提前世界标准时约7小时。
根据AE指数计算表面充电电位模式计算值U，由于图中的AE指数数据没有具体的数据文件，使用时可大致估算每小时的AE指数平均值，然后带入上述公式(1)-(3)计算表面充电电位模式计算值U；
根据表面充电电位模式计算值U及其与表面充电指数ISC的对应关系，获得风云四号卫星的表面充电指数。
本实施例详细计算结果如下表所示：
表1 2020年4月20日风云四号卫星位置的表面充电计算结果
Figure GDA0003116192820000071
Figure GDA0003116192820000081
图3给出了当日风云四号表面充电探测器的探测结果，分析图3可得，尽管2020年4月20日当天发生了地磁暴，AE指数也非常高，但风云四号位于磁尾对应的时区时(即22:00至次日6:00时间区间)，AE指数已经明显恢复，因此，风云四号没有探测到表面充电效应，这与本实施例1所得数据表1得到的结果一致。
在另一可选实施例中，以2020年3月20日风云四号的表面充电探测数据与当天的AE指数的相关特性为例。
在世界数据中心获取2020年3月20日的AE指数数据，下载2020年3月20日的探测数据(当日的地磁活动为平静到微扰，但该日的AE指数却显示出明显的扰动)，可得到该日的连续快视图，如图4所示，其中图中上部的包络线表示AE指数的大小。
将AE指数对应的国际标准时T0转化为卫星地方时TL。风云四号定点于东经104度，由公式(4)计算得到其地方时与国际标准时(UT)的对应关系，风云四号卫星的地方时提前世界标准时约7小时。
根据AE指数计算表面充电电位模式计算值U，由于图中的AE指数数据没有具体的数据文件，使用时可大致估算每小时的AE指数平均值，然后带入公式(1)-(3)进行计算表面充电电位模式计算值U；
根据表面充电电位模式计算值U及其与表面充电指数ISC的对应关系，获得风云四号卫星的表面充电指数。
本实施例详细计算结果如下表所示：
表2 2020年4月20日风云四号卫星位置的表面充电计算结果
Figure GDA0003116192820000082
Figure GDA0003116192820000091
图5给出了当日风云四号表面充电探测器的探测结果，分析图5可得，2020年3月20日，当风云四号处于地球磁尾时段，AE指数明显偏高。表2中也显示风云四号卫星在地方时3月20日22:00至21日6:00时间段，充电危险性较高，特别是地方时21日0:00至4:00，危险系数达到最高级别3级。可见，本实施例所得数据表2中的计算结果与图5中的风云四号实测数据相近。
为进一步验证本发明提供的地球同步轨道卫星表面充电效应的预警方法的可行性，在另一可选实施例中，将2020年3月1日至2020年4月26日的风云四号的表面充电事件的时间和规模进行了统计，同时给出了相关时段的AE指数及根据AE指数所得到的表面充电电位U模式，如表3所示。
表3 2020年3月1日至2020年4月26日模式计算值与风云四号实测值对比
Figure GDA0003116192820000092
Figure GDA0003116192820000101
Figure GDA0003116192820000111
如上表所示，表中的磁尾区为地球同步轨道卫星在不同地方时所处的磁尾中的位置：1表示过渡区-I，2表示扰动区，3表示过渡区-II；
U模式为表面充电电位模式计算值；
U实测为风云四号表面充电电位实测值；
ISC模式为充电指数模式值；
ISC实测为风云四号充电指数实测值。
图6为根据表3中计算出的表面充电电位模式值U模式、与风云四号探测得到表面充电电位实测值U实测的对比图，通过对比计算，得出了二者的相关系数的值为0.6964，说明U模式和U实测的数据相关性较强，进一步证明了本发明给出由AE指数推算地球同步轨道卫星的表面充电电位的方法具有较好的可行性。
表3中给出了依照本发明给出的方法得到的表面充电危险系数的模式计算值和风云四号实测值，针对卫星影响较明显的第三类指数事件，统计得到的相关结果如表4。
表4 2020年3月1日-2020年4月26日严重表面充电事件预警结果
I<sub>实测</sub>＝3 I<sub>实测</sub>&lt;3
I<sub>模式</sub>＝3 9 1
I<sub>模式</sub>&lt;3 1
由表4可得，本发明实施例对于危险级别的表面充电事件的预警的准确性为82.3％(9/11)，能够用于卫星表面充电危险事件的预警业务。
另一实施例提供一种地球同步轨道卫星表面充电效应的预警系统，如图7所示，包括数据获取单元，用于获取T0时刻的AE指数数据；第一数据处理单元，用于推算AE指数时刻所对应的卫星的地方时；其中第一数据处理单元设置有：数据获取模块，用于获取卫星定点的经度L；第一数据计算模块，用于将AE指数时刻T0所对应的卫星的地方时TL；第二数据处理单元，用于根据所述AE指数与表面充电电位U的关系，计算卫星的表面充电电位模式计算值U；第二数据处理单元设置有：时域划分模块，用于将卫星所在区域的地方时划分为四个时域；第二数据计算模块，用于根据不同时域AE指数与表面充电电位的不同的对应关系，获得所述表面充电电位模式计算值U；第三数据处理单元，用于根据所述表面充电电位模式计算值U，获得表面充电指数ISC；预警单元，用于根据所述表面充电指数ISC，对卫星表面充电危险情况进行预警；其中预警单元还包括报警模块，用于根据所述表面充电指数ISC值判定卫星表面充电危险等级，当所述卫星表面充电危险等级为高级时，发出警报。
综上，本发明所述的地球同步轨道卫星表面充电效应的预警方法及预警系统至少具备如下有益效果：
本发明所述的地球同步轨道卫星表面充电效应的预警方法及预警系统，首先获取AE指数数据；根据所述AE指数与表面充电电位模式计算值U的关系，计算所述表面充电电位模式计算值U；根据所述表面充电电位模式计算值U与表面充电指数ISC的对应关系，获得表面充电指数ISC；根据所述表面充电指数ISC，对卫星表面充电危险情况进行预警；其中根据AE指数获取卫星表面的充电电位计算值U，可以有效地对可能的表面充电效应提供可靠数据支撑，并且，根据卫星表面的充电电位模式的计算值U及表面充电指数ISC的对应关系，获得表面充电指数ISC，可以对卫星表面的充电效应情况进行预警，确定卫星表面的充电效应等级，快速对卫星表面的充电危险情况作出反应，进而能够在指导卫星采取防护措施和故障快速诊断方面提供可靠信息，对于卫星异常的分析和处置工作有参考作用。
本具体的实施例仅仅是对本发明的解释，而并不是对本发明的限制，本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出修改，但只要在本发明的权利要求范围内都受到专利法的保护。

Claims (13)

1.一种地球同步轨道卫星表面充电效应的预警方法，其特征在于，包括以下步骤：
获取T0时刻的AE指数数据；
推算AE指数时刻所对应卫星的地方时；
根据所述AE指数数据，计算所述卫星的表面充电电位模式计算值U；
根据所述表面充电电位模式计算值U，获得表面充电指数ISC；
根据所述表面充电指数ISC，对卫星表面充电危险情况进行预警；
根据所述AE指数数据，计算所述卫星的表面充电电位模式计算值U，包括：
将卫星所在区域的地方时划分为四个时域；
根据不同时域AE指数与表面充电电位的不同的对应关系，获得所述表面充电电位模式计算值U，其中：
当地方时处于6:00到21:00时域时，表面充电电位模式计算值U为0；
当地方时处于21:00到24:00时域时，表面充电电位模式计算值U为：
Figure FDA0003116192810000011
当地方时处于0:00到4:00时域时，表面充电电位模式计算值U为：
Figure FDA0003116192810000012
当地方时处于4:00到6:00时域时，表面充电电位模式计算值U为：
Figure FDA0003116192810000013
式中：AE为距离T0时刻最近的1小时内AE指数的平均值，单位为nT；
U为表面充电电位模式计算值，单位为V。
2.根据权利要求1所述的地球同步轨道卫星表面充电效应的预警方法，其特征在于，获取T0时刻的AE指数数据，还包括：根据获取的AE指数数据曲线，计算距离T0时刻最近的1小时内AE指数的平均值。
3.根据权利要求1所述的地球同步轨道卫星表面充电效应的预警方法，其特征在于，推算AE指数时刻所对应的卫星的地方时还包括以下步骤：
获取卫星定点的经度L；
根据如下公式获得AE指数时刻T0所对应的卫星的地方时TL：TL＝T0+L/15；
式中：卫星定点的经度L的单位为度。
4.根据权利要求1所述的地球同步轨道卫星表面充电效应的预警方法，其特征在于，根据所述表面充电电位模式计算值U，获得表面充电指数ISC，还包括：
当U≥-5000时，ISC＝1；
当-10000≤U<-5000时，ISC＝2；
当U≤-10000，ISC＞2。
5.根据权利要求1或4所述的地球同步轨道卫星表面充电效应的预警方法，其特征在于，根据所述表面充电指数ISC，对卫星表面充电危险情况进行预警，还包括：根据所述表面充电指数ISC值判定卫星表面充电危险等级。
6.根据权利要求5所述的地球同步轨道卫星表面充电效应的预警方法，其特征在于，根据所述表面充电指数ISC值判定卫星表面充电危险等级，还包括：当所述表面充电指数ISC＝1，则判断卫星表面充电危险等级为低级。
7.根据权利要求5所述的地球同步轨道卫星表面充电效应的预警方法，其特征在于，根据所述表面充电指数ISC值判定卫星表面充电危险等级，还包括：当所述表面充电指数ISC＝2，则判断卫星表面充电危险等级为中级。
8.根据权利要求7所述的地球同步轨道卫星表面充电效应的预警方法，其特征在于，根据所述表面充电指数ISC值判定卫星表面充电危险等级，还包括：当所述表面充电指数ISC＞2，则判断卫星表面充电危险等级为高级。
9.根据权利要求8所述的地球同步轨道卫星表面充电效应的预警方法，其特征在于，还包括：
当所述卫星表面充电危险等级为高级时，发出表面充电警报。
10.一种地球同步轨道卫星表面充电效应的预警系统，其特征在于，包括：
数据获取单元，用于获取T0时刻的AE指数数据；
第一数据处理单元，用于推算AE指数时刻所对应的卫星的地方时；
第二数据处理单元，用于根据所述AE指数数据，经如下步骤计算所述卫星的表面充电电位模式计算值U：
将卫星所在区域的地方时划分为四个时域；
根据不同时域AE指数与表面充电电位的不同的对应关系，获得所述表面充电电位模式计算值U，其中：
当地方时处于6:00到21:00时域时，表面充电电位模式计算值U为0；
当地方时处于21:00到24:00时域时，表面充电电位模式计算值U为：
Figure FDA0003116192810000031
当地方时处于0:00到4:00时域时，表面充电电位模式计算值U为：
Figure FDA0003116192810000032
当地方时处于4:00到6:00时域时，表面充电电位模式计算值U为：
Figure FDA0003116192810000033
式中：AE为距离T0时刻最近的1小时内AE指数的平均值，单位为nT；
U为表面充电电位模式计算值，单位为V；
第三数据处理单元，用于根据所述表面充电电位模式计算值U，获得表面充电指数ISC；预警单元，用于根据所述表面充电指数ISC，对卫星表面充电危险情况进行预警。
11.根据权利要求10所述的地球同步轨道卫星表面充电效应的预警系统，其特征在于，所述第一数据处理单元还设置有：
数据获取模块，用于获取卫星定点的经度L；
第一数据计算模块，用于将AE指数时刻T0所对应的卫星的地方时TL。
12.根据权利要求10所述的地球同步轨道卫星表面充电效应的预警系统，其特征在于，所述第二数据处理单元设置有：
时域划分模块，用于将卫星所在区域的地方时划分为四个时域；
第二数据计算模块，用于根据不同时域AE指数与表面充电电位的不同的对应关系，获得所述表面充电电位模式计算值U。
13.根据权利要求10所述的地球同步轨道卫星表面充电效应的预警系统，其特征在于，所述预警单元还包括报警模块，用于根据所述表面充电指数ISC值判定卫星表面充电危险等级，并且当卫星表面充电危险等级为高级时，进行报警。
CN202010600627.6A 2020-06-28 2020-06-28 地球同步轨道卫星表面充电效应的预警方法及预警系统 Active CN111703593B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010600627.6A CN111703593B (zh) 2020-06-28 2020-06-28 地球同步轨道卫星表面充电效应的预警方法及预警系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010600627.6A CN111703593B (zh) 2020-06-28 2020-06-28 地球同步轨道卫星表面充电效应的预警方法及预警系统

Publications (2)

Publication Number Publication Date
CN111703593A CN111703593A (zh) 2020-09-25
CN111703593B true CN111703593B (zh) 2021-08-17

Family

ID=72544232

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010600627.6A Active CN111703593B (zh) 2020-06-28 2020-06-28 地球同步轨道卫星表面充电效应的预警方法及预警系统

Country Status (1)

Country Link
CN (1) CN111703593B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113189617B (zh) * 2021-04-30 2023-04-14 北京软奇科技有限公司 一种太阳质子事件预警方法及系统、存储介质及终端
CN113158533B (zh) * 2021-04-30 2024-01-30 北京软奇科技有限公司 一种高能质子能谱计算方法及计算系统
CN114325129A (zh) * 2021-12-18 2022-04-12 国家卫星气象中心(国家空间天气监测预警中心) 一种适应fy3bvirr卫星的射出长波辐射值衰减订正方法
CN115600070B (zh) * 2022-10-14 2023-03-28 北京空间飞行器总体设计部 一种极轨卫星遭遇极光电子的预示方法
CN115856986B (zh) * 2023-03-03 2023-04-25 数字太空(北京)智能技术研究院有限公司 一种地球轨道航天器粒子辐射环境的计算方法及设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102117345A (zh) * 2009-12-31 2011-07-06 北京卫星环境工程研究所 一种进行卫星表面带电分析的方法
CN103761417A (zh) * 2013-12-24 2014-04-30 兰州空间技术物理研究所 地球同步轨道卫星表面电位计算方法
CN107657079A (zh) * 2017-08-29 2018-02-02 兰州空间技术物理研究所 一种卫星带电效应诱发故障与环境匹配分析方法
CN109325693A (zh) * 2018-09-28 2019-02-12 兰州空间技术物理研究所 一种航天器在轨表面带电效应风险量化评估方法
RU2683113C1 (ru) * 2018-03-26 2019-03-26 Александр Васильевич Тертышников Способ определения характеристик аврорального овала и состояния магнитного поля земли
CN110231643A (zh) * 2019-06-28 2019-09-13 国家卫星气象中心(国家空间天气监测预警中心) 高能电子暴事件的预报方法、装置及一种存储介质和设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102117345A (zh) * 2009-12-31 2011-07-06 北京卫星环境工程研究所 一种进行卫星表面带电分析的方法
CN103761417A (zh) * 2013-12-24 2014-04-30 兰州空间技术物理研究所 地球同步轨道卫星表面电位计算方法
CN107657079A (zh) * 2017-08-29 2018-02-02 兰州空间技术物理研究所 一种卫星带电效应诱发故障与环境匹配分析方法
RU2683113C1 (ru) * 2018-03-26 2019-03-26 Александр Васильевич Тертышников Способ определения характеристик аврорального овала и состояния магнитного поля земли
CN109325693A (zh) * 2018-09-28 2019-02-12 兰州空间技术物理研究所 一种航天器在轨表面带电效应风险量化评估方法
CN110231643A (zh) * 2019-06-28 2019-09-13 国家卫星气象中心(国家空间天气监测预警中心) 高能电子暴事件的预报方法、装置及一种存储介质和设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Calculating Auroral Oval Pattern by AE Index;CHEN Anqin等;《气象学报（英文版）》;20080131;第22卷(第1期);全文 *
FY2G卫星新一代高能带电粒子探测器观测数据分析;王春琴等;《地球物理学报》;20160930;第59卷(第9期);全文 *
基于COMSOL的GEO航天器表面充电仿真模型;武明志等;《科技创新导报》;20171231(第32期);全文 *
基于SPIS地球同步轨道航天器表面充电仿真;顾超超等;《微机与应用》;20171130;第36卷(第11期);全文 *

Also Published As

Publication number Publication date
CN111703593A (zh) 2020-09-25

Similar Documents

Publication Publication Date Title
CN111703593B (zh) 地球同步轨道卫星表面充电效应的预警方法及预警系统
Tozzi et al. A preliminary risk assessment of geomagnetically induced currents over the Italian territory
CN111242404B (zh) 一种强降雨诱发洪灾事件的极端性评估方法及系统
Bhui et al. Application of recurrence quantification analysis to power system dynamic studies
Mo et al. Estimating the extreme wind speed for regions in China using surface wind observations and reanalysis data
CN110231643B (zh) 高能电子暴事件的预报方法、装置及一种存储介质和设备
Naccarato et al. Improvements in the detection efficiency model for the Brazilian lightning detection network (BrasilDAT)
EP2921887A1 (en) Lightning-strike electric charge estimation system and method
Fan et al. A performance evaluation of the world wide lightning location network (WWLLN) over the Tibetan Plateau
CN113077164A (zh) 一种基于熵值法的输电线路落雷风险系数评估方法
Guo et al. Spatial and temporal variability of annual precipitation during 1958-2007 in Loess Plateau, China
Fabbro et al. GNSS positioning error forecasting in the Arctic: ROTI and Precise Point Positioning error forecasting from solar wind measurements
Gong et al. Dynamic-statistics combined forecast scheme based on the abrupt decadal change component of summer precipitation in East Asia
Herodotou et al. Distribution of lightning peak stroke currents in Ontario using an LLP system
JP2007232485A (ja) 架空地線の損傷区間推定方法
CN113359164B (zh) 一种低频率gnss电离层闪烁因子的有效性验证方法
Tofani et al. Lightning Activity Evaluation in South Sumatera Region
Jiang et al. Analysis of ionospheric vertical total electron content before the 1 April 2014 Mw 8.2 Chile earthquake
Hall Complexity signatures in the geomagnetic H component recorded by the Tromsø magnetometer (70° N, 19° E) over the last quarter of a century
Kappenman et al. Location-centered mitigation of lightning-caused disturbances
Sobouti The relationship between unique geomagnetic and auroral events
CN114638510B (zh) 高压直流干扰分析方法、装置、设备、系统及存储介质
Belov et al. Temporal and latitudinal dependence of the temperature effect for neutron component of cosmic ray
Li et al. Variations in extreme precipitation and relation to the Asia summer monsoon over the Qinling-Dabashan Mountains, China
Myint et al. Analysis of local geomagnetic index under the influence of equatorial electrojet (EEJ) at the equatorial Phuket geomagnetic station in Thailand

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