一种基于垂直探测数据修正等效风场的方法技术

本发明专利技术公开了一种基于垂直探测数据修正等效风场的方法，包括如下步骤：步骤1，调试电离层物理模型和中性风模型；步骤2，下载t时刻电离层垂直探测电离图，反演得到t时刻电离层参数；步骤3，下载物理模型输入参数；步骤4，按照电离层物理模型输入参数格式，计算电离层；步骤5，计算系数α：步骤6，计算差值，步骤7，推算t+Δt时刻输入电离层物理模型所需等效风场U（t+Δt），步骤8，计算t+Δt时刻电离层参数。本发明专利技术所公开基于垂直探测数据修正等效风场的方法，经过对电离层物理模型驱动因素中性风场的修正，能够使输出结果更加接近实际，显著提高模型计算精度。高模型计算精度。高模型计算精度。

【技术实现步骤摘要】
一种基于垂直探测数据修正等效风场的方法


[0001]本专利技术属于电离层环境预报领域，特别涉及该领域中的一种基于垂直探测数据修正等效风场的方法。

技术介绍

[0002]电离层是空间环境探测与研究中极为重要的区域，在空间工程、空间应用及空间环境的利用与开发中一直受到特别关注。电离层中存在足够的由于太阳辐射而离化的电离气体成分，以至于它能影响绝大部分波段的无线电波的传播特性，是影响空间信息链路畅通的重要环节。因此，对电离层状态的精确现报预报意义重大。
[0003]基于数值方法和计算技术，考虑大量物理、化学过程，通过连续性方程、动量方程及能量方程，物理模型可以模拟复杂的物理、化学过程，获得电离层状态参量如密度、速度、温度的时空分布。相比于经验模型，物理模型不依赖观测数据，能够从理论上洞察何种因素影响电离层的变化。然而，模型模拟结果受到多种因素影响，包括太阳辐射通量、中性成分密度、温度、风、电场等等，这些驱动因素均由经验模型给定，必然对计算结果带来误差。另外，理论模型不可能全面准确地考虑电离层当中的物理化学过程及扰动情况，再加上边界条件和初始条件的不确定性，及数值方法带来的误差，理论模型给出的结果往往与实际情况有很大偏差。

技术实现思路

[0004]本专利技术所要解决的技术问题就是针对电离层物理模型受驱动因素影响，模拟精度不高的问题，提供一种基于垂直探测数据修正等效风场的方法。
[0005]本专利技术采用如下技术方案：
[0006]一种基于垂直探测数据修正等效风场的方法，其改进之处在于，包括如下步骤：
[0007]步骤1，调试电离层物理模型和中性风模型HWM14；
[0008]步骤2，下载t时刻电离层垂直探测电离图，反演得到t时刻电离层参数：F2层峰值高度hmF2
ionosonde
(t)；
[0009]步骤3，下载物理模型输入参数，包括太阳活动指数F107、地磁活动指数AP；
[0010]步骤4，按照电离层物理模型输入参数格式，输入日期、时间、经纬度、太阳活动指数、地磁活动指数，利用HWM14计算的中性风场U
HWM14
(t)信息，计算电离层hmF2
model
(t)；
[0011]步骤5，其他输入参数不变，改变电离层物理模型输入参数中性风场，输入U1(t)计算得到hmF21(t)；输入U2(t)计算得到hmF22(t)，根据以下公式计算系数α：
[0012][0013]步骤6，计算HWM14风场与电离层物理模型计算所需等效风场的差值ΔU，计算方法如下：
[0014][0015]步骤7，推算t+Δt时刻输入电离层物理模型所需等效风场U(t+Δt)，计算方法如下：
[0016]U(t+Δt)＝U
HWM14
(t+Δt)+ΔU
[0017]步骤8，将步骤7中计算得到的U(t+Δt)作为物理模型输入参数，计算t+Δt时刻电离层参数。
[0018]本专利技术的有益效果是：
[0019]电离层中存在较多复杂的物理化学过程，而物理模型难以全面准确地描述这些过程，对电离层状态起至关重要作用的中性风场是导致计算结果不准确的重要原因。依照本专利技术所公开基于垂直探测数据修正等效风场的方法，经过对电离层物理模型驱动因素中性风场的修正，能够使输出结果更加接近实际，显著提高模型计算精度。作为电离层数值预报的核心，电离层物理模型模拟精度的提升，可为发展电离层数值现报、预报技术提供技术支撑。
附图说明
[0020]图1是本专利技术方法的流程示意图；
[0021]图2是用本专利技术方法修正前后模拟值与实际观测值的对比图。
具体实施方式
[0022]为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。
[0023]实施例1，本实施例公开了一种基于垂直探测数据修正等效风场的方法，基于垂直探测获取的电离层F2层峰值高度hmF2，修正等效风场(物理模型驱动因素)，提高模型模拟精度。
[0024]如图1所示，包括如下步骤：
[0025]步骤1，调试电离层物理模型(高度范围：90
‑
550km，时间步长：6分钟)和中性风模型HWM14(Horizontal Wind Model)；
[0026]步骤2，原始数据获取：
[0027]互联网下载t时刻电离层垂直探测电离图，本实施例中具体为：2020年1月1日，12LT，北京站(40.11
°
N,116.28
°
E)，反演(或经验模型)得到t时刻电离层参数：F2层峰值高度hmF2
ionosonde
(t)；
[0028]步骤3，互联网下载物理模型输入参数，包括太阳活动指数F107、地磁活动指数AP；
[0029]步骤4，hmF2获取，从电离图中获取hmF2参数信息，利用电离层物理模型，计算观测数据相应时段及站点电离层hmF2：
[0030]按照电离层物理模型输入参数格式，输入日期、时间、经纬度、太阳活动指数、地磁活动指数，利用HWM14计算的中性风场U
HWM14
(t)等信息，计算电离层hmF2
model
(t)；
[0031]步骤5，hmF2与中性风场U的关系求解，通过改变电离层物理模型输入参数(风场)计算相应hmF2，求得当前条件下hmF2与中性风场U的关系：
[0032]其他输入参数不变，改变电离层物理模型输入参数中性风场，输入U1(t)计算得到hmF21(t)；输入U2(t)计算得到hmF22(t)，根据以下公式计算系数α：
[0033][0034]步骤6，模型等效风场与经验计算值关系获取，根据等效风场与模型输出参数hmF2关系，获取模型等效风场与经验计算值关系：
[0035]计算HWM14风场与电离层物理模型计算所需等效风场的差值ΔU，计算方法如下：
[0036][0037]步骤7，计算下一时刻等效风场，并作为电离层物理模型输入：
[0038]推算t+Δt时刻输入电离层物理模型所需等效风场U(t+Δt)，计算方法如下：
[0039]U(t+Δt)＝U
HWM14
(t+Δt)+ΔU
[0040]步骤8，计算下一时刻电离层参数：
[0041]将步骤7中计算得到的U(t+Δt)作为物理模型输入参数，计算t+Δt时刻电离层参数。
[0042]图2是将本专利技术方法应用于电离层物理模型前后的计算结果与实际观测值的对比图。
[0043]综上所述，本专利技术提出了一种基于垂直探测数据修正等效风场的方法，基于单站电离层F2层峰高hmF2垂直探测数据，以及电离层本身的物理过程(即电离层峰高与风场存在一定相关性)，修正电离层物理模型关键输入参数：等效风场，使模型计算结果更加本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于垂直探测数据修正等效风场的方法，其特征在于，包括如下步骤：步骤1，调试电离层物理模型和中性风模型HWM14；步骤2，下载t时刻电离层垂直探测电离图，反演得到t时刻电离层参数：F2层峰值高度hmF2
ionosonde
(t)；步骤3，下载物理模型输入参数，包括太阳活动指数F107、地磁活动指数AP；步骤4，按照电离层物理模型输入参数格式，输入日期、时间、经纬度、太阳活动指数、地磁活动指数，利用HWM14计算的中性风场U
HWM14
(t)信息，计算电离层hmF2
m...

【专利技术属性】
技术研发人员：胡艳莉，徐彤，孙树计，邓忠新，朱梦言，
申请(专利权)人：中国电波传播研究所中国电子科技集团公司第二十二研究所，
类型：发明
国别省市：