一种地磁Dst指数计算方法及系统、存储介质及终端

一种地磁dst指数计算方法及系统、存储介质及终端
技术领域
1.本发明涉及地磁的技术领域，特别是涉及一种地磁dst指数计算方法及系统、存储介质及终端。


背景技术：

2.地磁暴是高速等离子体云到达地球空间后，引发的最具代表性的全球空间环境扰动事件。
3.地磁暴期间，高能粒子沉降和焦耳加热等过程使低层大气受热膨胀，引起高层大气密度增加；高层大气密度、成分和风场的变化，会引起电离层暴；磁层剧烈扰动时，磁尾中的热等离子体被加速向地球方向运动，形成热等离子体注入；带电粒子沿磁力线沉降，轰击高层大气，形成绚烂多彩的极光；磁层扰动期间，磁层中的电子可能被加速至很高的能量，引起全球范围的高能电子增强现象—高能电子暴。
4.当地磁场扰动时，磁场方向和大小的改变会影响它们之间的力矩，致使卫星的姿态发生变化。卫星的姿态发生变化后，通信卫星将无法正常通信，甚至有时可能会中断通信；气象卫星、军事卫星也无法监测地球。
5.强磁暴时，地磁场会发生剧烈扰动变化，变化的地磁场会在土壤电阻率高的地区产生每公里几伏特到十几伏特，持续时间从几分钟到几小时的地面电势(earth surface potential，esp)。而在高压、超高压输电系统中，由于电网变压器中性点直接接地，所以esp会在东西走向、长距离输电线路与大地构成的回路中产生地磁感应电流(geomagnetically induced currents，gic)。容易引起大型变压器半波饱和而缩短其使用寿命，极端情况下会使其烧毁而造成永久损坏。
6.地磁暴的强度可以表征太阳风暴中高速等离子体云的影响大小。地磁暴的强度等级一般用kp(planetary k，k为德文kennziffer)指数和dst(disturbance storm time)指数这两类地磁指数来划分。在研究中通常采用dst指数分级，而在预警应用中采用kp指数。它们大小和磁层、电离层等区域的活动有关，因此可以用于指示近地空间和磁层空间中的活动状况。
7.传统的地磁算法需要利用坐落在地球上一定纬度范围内、不同经度的地磁台站的磁场测量值。如dst指数利用了赤道附近南北磁纬度20
°
到35
°
之间的四个台站(honolulu,san juan,hermanus,kakioka)的测量值。在这个纬度范围内，地磁场主要受地球磁暴环电流的影响。因此，磁场的信号可以反映磁暴活动的发生以及强弱。
8.具体来说，计算dst指数需要先选取过往地磁平静日的测量数据，确定磁场的长期变化趋势h
base
(t)，再通过数据的拟合获取以天为周期的日变化量hsq(t)，这两种变化产生的原因可能是地磁极的漂移、电离层的地方时差异等，和磁暴活动关联不大。将观测磁场h
obs
(t)减去上述两种变化分量后，就可以得到主要有磁暴环电流贡献的残余磁场分量δh(t)。最后将各个台站的残余磁场分量归一化求和，得到dst指数。
9.然而，上述方法需要设置遍布全球的地磁台站，建立成本和维护成本较高，计算复
杂度较高。


技术实现要素：

10.鉴于以上所述现有技术的缺点，本发明的目的在于提供一种地磁dst指数计算方法及系统、存储介质及终端，基于太阳同步轨道晨昏卫星磁场来实现地磁指数的计算，快速精准，响应迅速，实用性强。太阳同步轨道晨昏卫星(如风云3e卫星)在环绕期间高度较高，受电离层电流影响较地面和更低轨卫星更小；且覆盖的地方时固定，不会引入由地磁扰动在地方时上的不均匀性导致的偏差；特别地，晨昏星覆盖了地方时06点和18点，这两个地方时可以分别反映了磁层环电流强度的最小值和最大值，因此晨昏星算法以及算法导出的指数更能反映全球性的地磁变化。
11.为实现上述目的及其他相关目的，本发明提供一种地磁dst指数计算方法，包括以下步骤：获取太阳同步轨道晨昏卫星的卫星位置和所采集的观测磁场；基于所述卫星位置，计算国际参考地磁场模型的模型磁场；计算所述模型磁场与所述观测磁场之间的残差磁场；基于地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型；基于所述残差磁场分布模型，对所述观测磁场进行去趋势化处理；在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场，作为地磁dst指数。
12.于本发明一实施例中，计算所述模型磁场与所述观测磁场之间的残差磁场采用以下任一方式：
13.1)δb＝b
obs-b
igrf
，其中，b
obs
表示所述观测磁场，b
igrf
表示国际参考地磁场模型，δb表示所述残差磁场；
14.2)δb＝b
obs-b
igrf-symh，其中，symh表示sym-h指数；
15.于本发明一实施例中，基于地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型包括以下步骤：
16.获取地磁平静日的残差磁场；
17.基于所述残差磁场构建磁场残差分布模型δb
ascend
(glongi，glati)和δb
descend
(glongi，glati)；其中δb
ascend
(glongi，glati)表示卫星处于升轨时期的残差磁场分布模型，δb
descend
(glongi，glati)表示卫星处于降轨时期的残差磁场分布模型；
18.将所述磁场残差分布模型δb
ascend
(glongi，glati)和δb
descend
(glongi，glati)插值到经度范围[-180
°
，+180
°
]之间，获取沿经纬度连续分布的残差磁场分布模型δb
ascend
(glong，glat)和δb
descend
(glong，glat)；其中δb
ascend
(glong，glat)表示卫星处于升轨时期的沿经纬度连续分布的残差磁场分布模型，δb
descend
(glong，glat)表示卫星处于降轨时期的沿经纬度连续分布的残差磁场分布模型。
[0019]
于本发明一实施例中，根据δb(t)＝b
obs
(t)-b
igrf
(glong，glat)-δb(glong，glat)对所述观测磁场进行去趋势化处理，其中，b
obs
(t)表示所述观测磁场，b
igrf
(glong，glat)表示国际参考地模型磁场，δb(glong，glat)表示残差磁场分布模型，在卫星处于升轨时期为δb
ascend
(glong，glat)，在卫星处于降轨时期为δb
descend
(glong，glat)。
[0020]
于本发明一实施例中，在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场包括以下步骤：
[0021]
在geo坐标系下，分别计算每个轨道升轨和降轨期间南北磁纬30度以内的残余磁
场z向平均值；
[0022]
计算升轨和降轨的残余磁场z向平均值的均值，作为所述dst指数。
[0023]
于本发明一实施例中，还包括当检测到地磁平静日时，将所述地磁平静日的观测磁场进行插值处理，并基于插值处理后的观测磁场更新所述残差磁场分布模型。
[0024]
本发明提供一种地磁dst指数计算系统，包括获取模块、磁场计算模块、残差计算模块、构建模块、去趋势化模块和dst指数计算模块；
[0025]
所述获取模块用于获取太阳同步轨道晨昏卫星的卫星位置和所采集的观测磁场；
[0026]
所述磁场计算模块用于基于所述卫星位置，计算国际参考地磁场模型的模型磁场；
[0027]
所述残差计算模块用于计算所述模型磁场与所述观测磁场之间的残差磁场；
[0028]
所述构建模块用于基于地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型；
[0029]
所述去趋势化模块用于基于所述残差磁场分布模型，对所述观测磁场进行去趋势化处理；
[0030]
所述dst指数计算模块用于在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场，作为地磁dst指数。
[0031]
本发明提供一种存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述的地磁dst指数计算方法。
[0032]
本发明提供一种地磁dst指数计算终端，包括：处理器及存储器；
[0033]
所述存储器用于存储计算机程序；
[0034]
所述处理器用于执行所述存储器存储的计算机程序，以使所述地磁dst指数计算终端执行上述的地磁dst指数计算方法。
[0035]
如上所述，本发明的地磁dst指数计算方法及系统、存储介质及终端，具有以下有益效果：
[0036]
(1)基于太阳同步轨道卫星磁场来实现地磁指数的计算，快速精准；
[0037]
(2)只需要一颗低轨卫星，不需要遍布全球的地磁台站，建立成本、维护成本较低；
[0038]
(3)能够通过调整参数(如纬度范围)快速建立新版本的地磁指数；
[0039]
(4)相比于地基观测，天基观测具有良好的南北对称性，更适于建立衡量全球地磁活动的指数；
[0040]
(5)特别地，太阳同步轨道晨昏卫星固定地围绕着06点和18点两个地方时，更能反映地磁扰动的上下限。
附图说明
[0041]
图1显示为本发明的地磁dst指数计算方法于一实施例中的流程图；
[0042]
图2显示为地磁sym-h指数于一实施例中的示意图；
[0043]
图3显示为卫星观测磁场和国际参考地磁场模型的残差磁场的z分量于一实施例中的示意图；
[0044]
图4显示为选取的地磁平静日内的升降轨残差磁场分布以及利用插值算法构建的残差磁场分布模型于一实施例中的示意图；
[0045]
图5所示为选取的地磁平静日之后的一天的原始残差磁场和经残差磁场分布模型修正后的磁场，以及卫星所在的纬度和磁纬度的示意图；
[0046]
图6显示为本发明的地磁dst指数计算方法获取的地磁dst指数序列和kyoto sym-h指数的对比示意图；
[0047]
图7显示为风云3号e星地磁dst指数序列和kyoto sym-h指数的散点图；
[0048]
图8显示为利用最初的地磁平静日建立残差磁场分布模型后得到的地磁dst指数以及使用了更新的残差磁场分布模型计算出的地磁dst指数的比对图；
[0049]
图9显示为图8中的两版本dst指数的散点图；
[0050]
图10显示为本发明的地磁dst指数计算系统于一实施例中的结构示意图；
[0051]
图11显示为本发明的地磁dst指数计算终端于一实施例中的结构示意图。
具体实施方式
[0052]
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需说明的是，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。
[0053]
需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其组件布局型态也可能更为复杂。
[0054]
本发明的地磁dst指数计算方法及系统、存储介质及终端利用太阳同步轨道卫星磁场计算地磁dst指数。太阳同步轨道晨昏卫星周期短，约100分钟，每天可以对经度进行扫描，对扫描磁场进行插值或平滑等处理，就可以得到总磁场经纬度的分布，通过限制纬度范围、国际参考地磁场模型的校正，能够尽可能得将磁层环电流体系以外的电流系统所带来的磁场贡献缩小，其系统建立成本、维护成本较低，计算精度高，实用性强。
[0055]
如图1所示，于一实施例中，本发明的地磁dst指数计算方法包括以下步骤：
[0056]
步骤s1、获取太阳同步轨道晨昏卫星的卫星位置和所采集的观测磁场b
obs
。
[0057]
具体地，选取太阳同步轨道晨昏卫星的卫星位置和卫星上的磁强计所采集的观测磁场。优选地，所述卫星可采用风云3号-e卫星。该卫星携带有4组巨磁阻磁强计，是一颗太阳同步轨道卫星，轨道倾角小，主要经历晨侧6点地方时以及昏测18点地方时。
[0058]
步骤s2、基于所述卫星位置，计算国际参考地磁场模型的模型磁场。
[0059]
具体地，计算国际参考地磁场模型在所述卫星位置处的模型磁场b
igrf
。
[0060]
步骤s3、计算所述模型磁场与所述观测磁场之间的残差磁场。
[0061]
于本发明一实施例中，根据δb＝b
obs-b
igrf
计算所述模型磁场与所述观测磁场之间的残差磁场δb。其中，残差磁场根据卫星的升降轨道分为两类：δb
ascend
和δb
descend
，以避免地方时或卫星本身的效应，其中下标ascend表示升轨，下标descend表示降轨。
[0062]
于本发明另一实施例中，为了进一步提升精度，尽可能减小固有地磁活动对残差的影响，利用kyoto sym-h指数进一步计算所述残差磁场，即δb＝b
obs-b
igrf-symh，其中symh
表示sym-h指数。具体地，sym-h指数与dst指数相似，采用了更多的地磁台站观测，时间精度更高，和dst指数相关性极高。
[0063]
步骤s4、基于地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型。
[0064]
具体地，地磁平静日是指连续24小时dst指数绝对值不大于10的时间区间。获取地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型。首先，基于所述残差磁场构建磁场残差随地理经度(geographic longitude，glong)和地理纬度(geographic latitude，glat)的分布模型δb
descend
(glongi，glati)和δb
ascend
(glongi，glati)，其中i为输入量在输入事件序列中的序号；然后，将所述磁场残差分布δb
descend
(glongi，glati)和δb
ascend
(glongi，glati)插值到经度范围[-180
°
，+180
°
]之间，获取沿经纬度连续分布的残差磁场分布模型δb
ascend
(glong，glat)和δb
descend
(glong，glat)。其中δb
ascend
(glong，glat)表示卫星处于升轨时期的沿经纬度连续分布的残差磁场分布模型，δb
descend
(glong，glat)表示卫星处于降轨时期的沿经纬度连续分布的残差磁场分布模型。
[0065]
步骤s5、基于所述残差磁场分布模型，对所述观测磁场进行去趋势化处理。
[0066]
具体地，根据δb(t)＝b
obs
(t)-b
igrf
(glong，glat)-δb(glong，glat)对所述观测磁场进行去趋势化处理，其中，b
obs
(t)表示所述观测磁场，b
igrf
(glong，glat)表示国际参考地磁场，δb(glong，glat)表示系统性残差，在卫星处于升轨时期为δb
ascend
(glong，glat)，在卫星处于降轨时期为δb
descend
(glong，glat)。
[0067]
步骤s6、在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场，作为地磁dst指数。
[0068]
具体地，在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场包括以下步骤：
[0069]
61)在geo坐标系下，分别计算每个轨道升轨和降轨期间南北磁纬30度以内的残余磁场δb的z向平均值。其中，需要对升轨和降轨分别计算南北磁纬30度以内的残余磁场δb的z向平均值。
[0070]
62)计算升轨和降轨的残余磁场z向平均值的均值，作为所述z向平均残余磁场
[0071]
于本发明一实施例中，为了进一步提升精度，本发明的地磁dst指数计算方法还包括当检测到地磁平静日时，将所述地磁平静日的观测磁场进行插值处理，并基于插值处理后的观测磁场更新所述残差磁场分布模型，以用于后续的地磁dst指数计算。
[0072]
下面通过具体实施例来进一步阐述本发明的地磁dst指数计算方法。
[0073]
图2所示为地磁sym-h指数，图中灰色阴影区域为选取的地磁平静日。图3所示为卫星观测磁场和国际参考地磁场模型的残差磁场(δb＝b
obs-b
igrf
)的z分量。图4中(a)部分和(b)部分所示为选取的地磁平静日内的升降轨残差分布，颜色深浅代表δbz的大小。图(c)部分和(d)部分为根据(a)部分和(b)部分中展示的残差磁场分布利用插值算法构建的残差磁场分布模型。
[0074]
图5的上半部分所示为选取的地磁平静日之后的一天的原始残差磁场和经残差磁场分布模型修正后的磁场。图5的下半部分所示为卫星所在的磁纬度，图中阴影区域为磁纬度大于30度的区域。由图可知，大幅的磁场残差被残差磁场分布模型修正掉了，剩余残差比
较大的区域集中在高纬(磁纬度大于30度)区域，对赤道附近的影响不大。
[0075]
图6所示为本发明的地磁dst指数计算方法获取的地磁dst指数序列和国际公认的kyoto (kyoto world data center for geomagnetism，京都地磁数据中心)sym h指数的对比。图7所示为风云3号e星地磁dst指数序列和kyoto dst指数的散点图，横轴为地面台站观测得出的dst指数(kyoto发布)，纵轴为风云3号e星推算得出的dst指数，二者相关系数可达0.86。
[0076]
图8显示为利用最初的地磁平静日建立残差磁场分布模型后得到的地磁dst指数以及使用了更新的残差磁场分布模型计算出的地磁dst指数。图9是图8中的两版本dst指数的散点图。由图可知，更新后的地磁dst指数和kyoto dst指数间均方根偏移为3.2，小于更新前残差磁场分布模型推算出的均方根偏移3.4。
[0077]
如图10所示，于一实施例中，本发明的地磁dst指数计算系统包括获取模块101、磁场计算模块102、残差计算模块103、构建模块104、去趋势化模块105和dst指数计算模块106。
[0078]
所述获取模块101用于获取太阳同步轨道晨昏卫星的卫星位置和所采集的观测磁场。
[0079]
所述磁场计算模块102与所述获取模块101相连，用于基于所述卫星位置，计算国际参考地磁场模型的模型磁场。
[0080]
所述残差计算模块103与所述获取模块101和所述磁场计算模块102相连，用于计算所述模型磁场与所述观测磁场之间的残差磁场。
[0081]
所述构建模块104与所述残差计算模块103相连，用于基于地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型。
[0082]
所述去趋势化模块105与所述构建模块104相连，用于基于所述残差磁场分布模型，对所述观测磁场进行去趋势化处理。
[0083]
所述dst指数计算模块106与所述去趋势化模块105相连，用于在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场，作为地磁dst指数。
[0084]
其中，获取模块101、磁场计算模块102、残差计算模块103、构建模块104、去趋势化模块105和dst指数计算模块106的结构和原理与上述地磁dst指数计算方法中的步骤一一对应，故在此不再赘述。
[0085]
需要说明的是，应理解以上装置的各个模块的划分仅仅是一种逻辑功能的划分，实际实现时可以全部或部分集成到一个物理实体上，也可以物理上分开。且这些模块可以全部以软件通过处理元件调用的形式实现；也可以全部以硬件的形式实现；还可以部分模块通过处理元件调用软件的形式实现，部分模块通过硬件的形式实现。例如，x模块可以为单独设立的处理元件，也可以集成在上述装置的某一个芯片中实现，此外，也可以以程序代码的形式存储于上述装置的存储器中，由上述装置的某一个处理元件调用并执行以上x模块的功能。其它模块的实现与之类似。此外这些模块全部或部分可以集成在一起，也可以独立实现。这里所述的处理元件可以是一种集成电路，具有信号的处理能力。在实现过程中，上述方法的各步骤或以上各个模块可以通过处理器元件中的硬件的集成逻辑电路或者软件形式的指令完成。
[0086]
例如，以上这些模块可以是被配置成实施以上方法的一个或多个集成电路，例如：
一个或多个特定集成电路(application specific integrated circuit，简称asic)，或，一个或多个微处理器(digital signal processor，简称dsp)，或，一个或者多个现场可编程门阵列(field programmable gate array，简称fpga)等。再如，当以上某个模块通过处理元件调度程序代码的形式实现时，该处理元件可以是通用处理器，例如中央处理器(central processing unit，简称cpu)或其它可以调用程序代码的处理器。再如，这些模块可以集成在一起，以片上系统(system-on-a-chip，简称soc)的形式实现。
[0087]
本发明的存储介质上存储有计算机程序，该程序被处理器执行时实现上述的地磁dst指数计算方法。所述存储介质包括：rom、ram、磁碟、u盘、存储卡或者光盘等各种可以存储程序代码的介质。
[0088]
如图11所示，于一实施例中，本发明的地磁dst指数计算终端包括：处理器111及存储器112。
[0089]
所述存储器112用于存储计算机程序。
[0090]
所述存储器112包括：rom、ram、磁碟、u盘、存储卡或者光盘等各种可以存储程序代码的介质。
[0091]
所述处理器111与所述存储器112相连，用于执行所述存储器112存储的计算机程序，以使所述地磁dst指数计算终端执行上述的地磁dst指数计算方法。
[0092]
优选地，所述处理器111可以是通用处理器，包括中央处理器(central processing unit，简称cpu)、网络处理器(network processor，简称np)等；还可以是数字信号处理器(digital signal processor，简称dsp)、专用集成电路(application specific integrated circuit，简称asic)、现场可编程门阵列(field programmable gate array，简称fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
[0093]
综上所述，本发明的地磁dst指数计算方法及系统、存储介质及终端基于太阳同步轨道卫星晨昏磁场来实现地磁指数的计算，快速精准；只需要一颗低轨卫星，不需要遍布全球的地磁台站，建立成本、维护成本较低；能够通过调整参数(如纬度范围)快速建立新版本的地磁指数；相比于地基观测，天基观测具有良好的南北对称性，特别地，太阳同步轨道晨昏卫星固定地围绕着06点和18点两个地方时，更能反映地磁扰动的上下限，更适于建立衡量全球地磁活动的指数。所以，本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
[0094]
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。

技术特征：
1.一种地磁dst指数计算方法，其特征在于：包括以下步骤：获取太阳同步轨道晨昏卫星的卫星位置和所采集的观测磁场；基于所述卫星位置，计算国际参考地磁场模型的参考磁场；计算所述国际参考地磁场模型的参考磁场与所述观测磁场之间的残差磁场，并基于地磁平静时期的残差磁场构建残差磁场经验数值模型；基于所述残差磁场分布模型，对所述观测磁场进行去趋势化处理；在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场，作为地磁dst指数。2.根据权利要求1所述的地磁dst指数计算方法，其特征在于：计算所述模型磁场与所述观测磁场之间的残差磁场采用以下任一方式：1)δb＝b
obs-b
igrf
，其中，b
obs
表示所述观测磁场，b
igrf
表示所述国际参考地磁场模型磁场，δb表示所述残差磁场；2)δb＝b
obs-b
igrf-symh，其中，symh表示sym-h指数；3.根据权利要求1所述的地磁dst指数计算方法，其特征在于：基于地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型包括以下步骤：获取地磁平静日的残差磁场；基于所述残差磁场构建磁场残差分布模型δb
ascend
(glong
i
，glat
i
)和δb
descend
(glong
i
，glat
i
)；其中δb
ascend
(glong
i
，glat
i
)表示卫星处于升轨时期的残差磁场分布模型，δb
descend
(glong
i
，glat
i
)表示卫星处于降轨时期的残差磁场分布模型；将所述磁场残差分布模型δb
ascend
(glong
i
，glat
i
)和δb
descend
(glong
i
，glat
i
)插值到经度范围[-180
°
，+180
°
]之间，获取沿经纬度连续分布的残差磁场分布模型δb
ascend
(glong，glat)和δb
descend
(glong，glat)；其中δb
ascend
(glong，glat)表示卫星处于升轨时期的沿经纬度连续分布的残差磁场分布模型，δb
descend
(glong，glat)表示卫星处于降轨时期的沿经纬度连续分布的残差磁场分布模型。4.根据权利要求1所述的地磁dst指数计算方法，其特征在于：根据δb(t)＝b
obs
(t)-b
igrf
(glong，glat)-δb(glong，glat)对所述观测磁场进行去趋势化处理，其中，b
obs
(t)表示所述观测磁场，b
igrf
(glong，glat)表示国际参考地模型磁场，δb(glong，glat)表示残差磁场分布模型，在卫星处于升轨时期为δb
ascend
(glong，glat)，在卫星处于降轨时期为δb
descend
(glong，glat)。5.根据权利要求1所述的地磁dst指数计算方法，其特征在于：在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场包括以下步骤：在geo坐标系下，分别计算每个轨道升轨和降轨期间南北磁纬30度以内的残余磁场z向平均值；计算升轨和降轨的残余磁场z向平均值的均值，作为所述z向平均残余磁场。6.根据权利要求1所述的地磁dst指数计算方法，其特征在于：还包括当检测到地磁平静日时，将所述地磁平静日的观测磁场进行插值处理，并基于插值处理后的观测磁场更新所述残差磁场分布模型。7.一种地磁dst指数计算系统，其特征在于：包括获取模块、磁场计算模块、残差计算模块、构建模块、去趋势化模块和dst指数计算模块；
所述获取模块用于获取太阳同步轨道晨昏卫星的卫星位置和所采集的观测磁场；所述磁场计算模块用于基于所述卫星位置，计算国际参考地磁场模型的模型磁场；所述残差计算模块用于计算所述模型磁场与所述观测磁场之间的残差磁场；所述构建模块用于基于地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型；所述去趋势化模块用于基于所述残差磁场分布模型，对所述观测磁场进行去趋势化处理；所述dst指数计算模块用于在geo坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的z向平均残余磁场，作为地磁dst指数。8.一种存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1至6中任一项所述的地磁dst指数计算方法。9.一种地磁dst指数计算终端，其特征在于，包括：处理器及存储器；所述存储器用于存储计算机程序；所述处理器用于执行所述存储器存储的计算机程序，以使所述地磁dst指数计算终端执行权利要求1至6中任一项所述的地磁dst指数计算方法。

技术总结
本发明提供一种地磁Dst指数计算方法及系统、存储介质及终端，包括获取太阳同步轨道晨昏卫星的卫星位置和所采集的观测磁场；基于所述卫星位置，计算国际参考地磁场模型的模型磁场；计算所述模型磁场与所述观测磁场之间的残差磁场；基于地磁平静日的残差磁场，构建沿经纬度连续分布的残差磁场分布模型；基于所述残差磁场分布模型，对所述观测磁场进行去趋势化处理；在GEO坐标系下，计算每个轨道南北磁纬30度以内的残余磁场的Z向平均残余磁场，作为地磁Dst指数。本发明的地磁Dst指数计算方法及系统、存储介质及终端基于太阳同步轨道卫星磁场来实现地磁指数的计算，快速精准，实用性强。实用性强。实用性强。


技术研发人员：赵玖桐 宗秋刚 王劲松 任杰
受保护的技术使用者：北京大学
技术研发日：2021.11.29
技术公布日：2022/3/8