一种基于非相干散射雷达的热层温度与氧原子密度反演方法与流程

本发明属于地球空间环境探测与研究领域，特别涉及该领域中的一种基于非相干散射雷达的热层大气温度与氧原子密度反演方法。

背景技术：

热层大气是指中层顶(约90km)到约500km的地球中性大气层，它是地球电离层的背景大气层。受太阳、地磁活动以及低层大气等因素的影响，热层大气温度、密度、成分等具有复杂的时空变化特性。热层大气温度与密度的观测反演具有重要的科学研究与工程应用价值。热层大气是国际空间站、神舟飞船、航天飞机等低轨航天器驻留和运行的空间环境，热层大气变化对低轨航天器运行具有重要影响。比如：热层大气密度直接决定了低轨道航天器运行的大气阻力，实验表明热层大气密度变化引起的低轨航天器轨道变化和摄动可达数百m；热层大气温度对低轨航天器的材料选择、环境适应性设计等有重要影响。

目前的热层大气观测手段非常有限，主要观测手段是通过卫星搭载光谱仪、大气成分探测器、宽带辐射计等实现热层探测，此外通过卫星等低轨航天器运行轨道也可推算热层大气密度等。在地基探测方面，通过发射高空探空火箭可探测火箭飞行区域的热层大气密度等参数，地面被动光学(气辉、fpi光谱仪)等可探测局地(如250km)的热层大气密度。此外，非相干散射雷达作为目前地面探测电离层的最强大手段，它可以直接探测电离层电子密度、电子与离子温度，间接推导热层大气温度、密度、热层成分、运动(风)等。目前全世界科学研究与工程应用领域普遍使用的热层大气模型msis就是基于卫星质谱仪、地基探空火箭与非相干散射雷达的观测数据而开发的。

由于观测手段的不足，我国的热层大气监测手段与数据严重不足，在热层大气反演与建模方面尚属于空白，急需加强。在子午工程(国家“十一五”大科学工程)支持下，我国在云南曲靖建设了国内首套非相干散射雷达，为我国开展热层大气温度与密度的反演提供了条件。

技术实现要素：

本发明所要解决的技术问题就是提供一种基于非相干散射雷达的热层温度与氧原子密度反演方法。

本发明采用如下技术方案：

一种基于非相干散射雷达的热层温度与氧原子密度反演方法，其改进之处在于，包括如下步骤：

(1)热层温度与氧原子密度模型、氧原子离子热量平衡方程：

研究表明，热层大气温度随高度的变化特点是从较低的某一高度(约120km)出发，开始温度梯度很大，随后随高度增加而减小。在很高的高度上，热层大气温度趋近于恒定值(通常称为逃逸层温度tn∞)，逃逸层温度与太阳活动密切相关，从太阳活动低年的600k左右到高年的2000k左右。

bates(空间物理学家的姓名)的热层温度经验模型如下：

其中tn120为120km的热层温度，z为相对120km的高度，s表示热层温度形态因子；

热层氧原子密度的经验模型如下：

其中m(o)为氧原子质量，g为重力加速度，k是波尔兹曼常数，z0是参考高度，[o]z0是参考高度处的氧原子密度，可见中性氧原子密度主要由热层温度、参考高度氧原子密度与高度确定；

在高度150-500km的电离层f层，主要离子成分为氧原子离子，忽略氮分子与氧分子的热损失、输运过程，那么电子-氧原子离子碰撞能量损失lei与氧原子离子-氧原子碰撞能量损失lin满足能量平衡，即lei＝lin；

根据banks(空间物理学家的姓名)的经验公式：

lin＝2.1×10^-15×[o]×(tn+ti)^1/2×(ti-tn)×ne(3-2)

这里ne为电子密度，te是电子温度，ti是离子温度，tn是热层中性大气温度；

(2)反演初值与收敛条件的确定：

根据(1)式，未知变量为逃逸层温度tn∞、120km热层温度tn120和温度梯度s，其中tn120为常数，其变化范围为350-450k，由于120km位于电离层e区，在夜间热层温度与离子温度处于热平衡，因此通过对非相干散射雷达夜间的离子温度观测数据进行统计分析可得出；在高度大于300km时，热层温度接近逃逸层温度，温度梯度很小，为0.02k/km，近似为零，难以有效反演温度梯度；根据(2)式，未知变量为热层温度与参考高度处的氧原子密度；根据(3)式，未知变量为氧原子密度与热层温度；利用msis模型(massspectrometerandincoherentscatterradar，质谱与非相干散射雷达热层大气模型，是国外开发的且目前比较广泛使用的模型)可计算逃逸层温度与参考高度氧原子密度初值；

非相干散射雷达可提供100-500km内电离层电子密度、电子温度、离子温度剖面的观测数据，即观测方程个数多于反演参数个数，属于超定方程组，可采用最小二乘法迭代求解；即假设tn∞、s与[o]z0的初值，通过方程(1)-(3)，求解离子温度tiinv，并与非相干散射雷达实测离子温度timea比较，在满足迭代次数与反演精度条件下，即可得到一组tn∞、s与[o]z0的反演值；

迭代次数与反演精度有关，也影响了计算量与反演速度，一般迭代次数数百至数千次。离子温度反演精度一般小于10％。迭代收敛条件可表示如下：

n为与高度范围、高度步长有关的高度数目，δtimea,i为相对于非相干散射雷达观测值的第i个高度的离子温度偏差；

从(1)式可推导出热层温度梯度与热层温度形态参数的表达式如下：

在300km以上，热层温度趋近于常数，温度梯度接近于零，此时无法有效反演热层温度形态因子；

从(3)式可知，未知量为氧原子密度[o]、离子温度与热层温度差δtin，从(3)式可推导出氧原子密度随着离子温度与热层温度差而变化的表达式如下：

可见δtin较小时，氧原子密度变化剧烈，误差增大；当离子温度与热层温度相近时，较小的温度误差可能导致很大的氧原子密度变化或误差(见图1)，甚至得不到有效的氧原子密度。比如，夜间时，离子温度与热层温度几乎相等，电子温度与离子温度接近，此时(3)式能量损失项为小项，无法反演得到有效的氧原子密度。

基于以上考虑，若利用非相干散射雷达观测的160-500km的电子密度、电子温度、离子温度数据同时反演tn∞、s与[o]z0三个参数将存在很大误差，或者可信度大大降低，因此我们采用以下分步反演法。

(3)分步反演法：

(31)对于300km以上高度，设温度梯度为常数值s0，利用非相干散射雷达350-500km的电子密度、电子温度、离子温度观测数据，可以初步反演tn∞与[o]z0；

(32)对于160-500km，将以上步骤的[o]z0作为准确值，进一步反演tn∞和s。

进一步的，参考高度需要满足扩散平衡，取150-250km。

本发明的有益效果是：

本发明所公开的基于非相干散射雷达的热层温度与氧原子密度反演方法，对提升我国的热层大气探测与研究水平具有重要意义，与现有其它的热层探测反演方法相比，具有简单、快捷、有效的突出优点。该方法的研究提出，有助于充分发挥已建大型设备(非相干散射雷达)的作用，具有重要的科学与应用价值。

附图说明

图1是氧原子密度变化量与离子温度—热层温度差的变化关系图；

图2是曲靖非相干散射雷达实测的电子密度、电子温度、离子温度剖面的日变化图；

图3是本方法基于曲靖非相干散射雷达观测数据反演的热层温度图；

图4是本方法基于曲靖非相干散射雷达观测数据反演的氧原子密度图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

实施例1，本发明所公开的方法已经基于子午工程曲靖非相干散射雷达进行了非公开保密应用，该雷达的主要参数包括：工作频率为500mhz、发射峰值功率为2mw、最大占空比为5％、脉冲重复周期为1-20ms、天线增益为41db、系统噪声温度为150k。其具体实施步骤如下：

1)首先，获取非相干散射雷达的电子密度、电子温度、离子温度实测数据，逃逸层温度(1000k)与氧原子密度(10¹⁵/m³)初值，离子温度误差值(取实测离子温度的10％)，最大迭代次数1000次。曲靖非相干散射雷达实测的电子密度、电子温度、离子温度剖面日变化见图2。

2)假设温度梯度为常数(0.02k/km)，通过迭代反演，可得到tn∞与[o]z0的初步反演值，分别为1100k与1.76x10¹⁵/m³。

3)根据2)获取的[o]z0反演值，进一步反演得到tn∞与s修正值，作为最终反演结果。

4)根据公式(1)和(2)计算得到热层温度、氧原子密度剖面，见图3与图4。