1.一种基于中性原子成像探测的极光亚暴模拟方法，包括以下步骤：
建立以离子投掷角扩散引导的环电流离子通量分布动态演化的中性原子成像模拟方程；
利用中性原子成像模拟方程，正演地磁宁静期间和极光亚暴期间不同离子投掷角分布函数对应的黄道面内的中性原子仿真图像，以重现极光亚暴过程中增长相阶段和恢复相阶段的环电流离子通量分布形态。
2.根据权利要求1所述的基于中性原子成像探测的极光亚暴模拟方法，其特征在于，所述建立以离子投掷角扩散引导的环电流离子通量分布动态演化的中性原子图像模拟方程，具体包括以下步骤：将中性原子图像的每个像素中记录的计数C(δ，ε)、仰角δ和方位角ε，表示在中性原子图像模拟方程中：
式中，ΔE为离子能道宽度；ΔT为像素的积分时间；ΔΩ是体积元素指向仰角δ和方位角ε像素的立体角； 为积分体积元处的离子微分通量；A(δ,ε)为探测器的响应函数；σ(E)为电荷交换截面； 为外逸层中性原子密度，其中r为地心距离，为地方时，θ为磁纬度；dV是沿探测器视线方向的体积元积分，L为磁壳指数，E为能量，α为磁层中性原子图像发射体元内的离子投掷角。
3.根据权利要求2所述的基于中性原子成像探测的极光亚暴模拟方法，其特征在于，所述积分体积元处的离子微分通量 为：
其中，L为磁壳参数，为地方时，θ为磁纬度，E为能量，α为磁层中性原子图像发射体元内的离子投掷角，αeq为赤道离子的投掷角；κ＝5.5；Emax0为环电流离子通量最大的能量初始值；
κ+1
e＝(1+1/κ) ；
其中, 为赤道离子通量；
为环电流离子通量的地方时分布函数：
其中, ξ＝0.73，且环电流离子通量的地方时 分布对称；
fL为环电流离子通量的磁壳参数分布函数：
其中,L为磁壳参数，L11为环电流第一边界值，L22为环电流第二边界值，γ1为第一模型参数，γ2为第二模型参数，γ为第三模型参数；其中，
其中，Kp11为第一地磁活动指数，Kp22为第二地磁活动指数，Kp11＝5.5，
Kp22＝0；
fα为环电流离子的投掷角分布函数；
其中，磁层中性原子图像发射体元内的离子投掷角α用赤道离子的投掷角αeq表示。
4.根据权利要求3所述的基于中性原子成像探测的极光亚暴模拟方法，其特征在于，所述第一模型参数γ1取值为0.53，第二模型参数γ2取值为0.88，第三模型参数γ取值为
1.16。
5.根据权利要求2所述的基于中性原子成像探测的极光亚暴模拟方法，其特征在于，所述外逸层中性原子密度 为：
其中，设定中性原子的径向通量在强地磁风暴期间是守恒的，n0为中性密度常数且取值‑3
为1600cm ，exp表示指数函数，r为地心距离，地心距离r以外逸层高度RE为单位，θ为磁纬度，为地方时，中间值a0＝1.78。
6.根据权利要求2所述的基于中性原子成像探测的极光亚暴模拟方法，其特征在于，所述中性原子图像由高能中性原子成像仪记录；
所述高能中性原子成像仪，在赤道环电流沉降到极光区的能量离子与从外逸层蒸发的低能中性原子进行电荷交换后，记录中性原子图像；
所述高能中性原子成像仪的最小能量通道为4keV。
7.根据权利要求1所述的基于中性原子成像探测的极光亚暴模拟方法，其特征在于，在正演地磁宁静期间和极光亚暴期间不同离子投掷角分布函数对应的黄道面内的中性原子仿真图像后，还包括以下步骤：
在地磁宁静期间，利用黄道面内的中性原子仿真图像，对环电流区域储存的能量离子通量分布、质谱和能谱进行测量评估；
基于地磁宁静期间，黄道面内的中性原子仿真图像对环电流区域储存的能量离子通量分布、质谱和能谱测量评估结果，监测太阳风弓激波扰动增长相阶段，高能中性原子通量在环电流区域的减少情况和在极光带区域的增强情况，以及对应的AE指数变化，以评估极光亚暴沉降离子的加速过程，并分析对应的加速机制。
8.根据权利要求1所述的基于中性原子成像探测的极光亚暴模拟方法，其特征在于，在正演地磁宁静期间和极光亚暴期间不同离子投掷角分布函数对应的黄道面内的中性原子仿真图像后，还包括以下步骤：
在极光亚暴过程中的恢复相阶段，通过中性原子图像监测极光带区域和环电流区域的能量中性原子通量变化，监测磁场偶极化过程中注入环电流的能量离子成分和能谱，以评估极光亚暴过程中，环电流能量离子的沉降损失和注入更新情况。
9.根据权利要求1所述的基于中性原子成像探测的极光亚暴模拟方法，其特征在于，在正演地磁宁静期间和极光亚暴期间不同离子投掷角分布函数对应的黄道面内的中性原子仿真图像后，还包括以下步骤：
基于极光亚暴过程中恢复相阶段的环电流离子通量分布形态，评估连续极光亚暴产生的几率。
10.一种基于中性原子成像探测的极光亚暴宏观模型，其特征在于，包括：模拟模块和仿真模块，其中，
所述模拟模块，用于建立以离子投掷角扩散引导的环电流离子通量分布动态演化的中性原子成像模拟方程；
所述仿真模块，用于利用中性原子成像模拟方程，正演地磁宁静期间和极光亚暴期间不同离子投掷角分布函数对应的黄道面内的中性原子仿真图像，以重现极光亚暴过程中增长相阶段和恢复相阶段的环电流离子通量分布形态。