一种全国对流层大气折射率剖面的预测方法与流程

技术特征：

1.一种全国大气折射率剖面的预测方法，其特征在于：建立的全国对流层大气折射率剖面预测数据库应当包含以下两个方面的内容：一是对全国大气环境按照一定的经纬度间隔进行栅格划分后各个栅格的经纬度信息；二是每个栅格中各个月份的大气分段剖面模型系数G、C1与地面折射率N0之间的关系系数。

2.根据权利要求1所述的全国大气折射率剖面的预测方法，其特征在于：首先根据我国的地理位置按照经度1°划分，纬度0.5°划分，将我国划分为5200个格栅，对于栅格内没有探空数据的栅格选用径向基函数插值算法求得该栅格内的大气剖面；其次选用精度较高的分段模型进行大气折射率预测，分段模型中的主要参数有：G为近地面1km大气层的折射率垂直梯度，C1为地面以上1km到9km大气层的指数衰减系数，C9为海拔9km高度以上大气层的指数衰减系数，因为C9在全国范围内的基本变化不大，因此直接使用C9＝0.1434，最后根据每一次探空数据统计出模型中的系数：G＝N1-N0，N0为地面折射率，N1为地面以上1km高度处的折射，根据10的探测数据对N0与G，N0与C1进行回归，建立一般情况下N0与G，N0与C1的关系式，再将每一个栅格中的关系式建库得到全国地面大气折射率与剖面模型系数关系数据库。

3.根据权利要求1所述的全国大气折射率剖面的预测方法，其特征在于：全国对流层大气折射率剖面的预测数据库的建立所包含的步骤主要有：分析全国大气折射率的变化特点、地理栅格划分、基础探测数据的处理、大气分段模型参数的计算、野值处理、有效数据的整理、径向基函数插值、数据统计及数据库的建立。

4.根据权利要求1所述的全国大气折射率剖面的预测方法，其特征在于：全国大气折射率的变化特点及地理环境栅格划分为：首先由我国常用雷达所处的位置可以确定我国的大气环境分布范围是：经度70°～135°，纬度15°～55°，另外，由无线电气象环境的变化特点可知，大气折射率随经度变化较小，随纬度变化较大，所以在全国大气环境范围内经度按照3°进行划分，纬度按照0.5°进行划分，这样一共将全国划分为5200个栅格，那就可以用该栅格中的大气折射率变化代表这一区域的变化，我国的地面大气探测站一共有500多个，而高空探测站有100多个，所以大多数栅格中是没有数据的，则这些栅格中参数就可用空间插值算法来得到。

5.根据权利要求1所述的全国大气折射率剖面的预测方法，其特征在于：基础探测数据的处理为：本发明专利使用的是我国100多个高空探空站10年的有效探空数据，这些高空探测站每天零点探测一次，12点探测一次，探空站探测的是气压p、温度t、相对湿度RH等一些气象参数，因为对流层的折射率计算公式为：

式中，N为大气折射率，单位为N；p为压强，单位为hPa；T表示的是热力学温度，其单位为K，它和摄氏温度之间的转换关系可写为T＝273+t；e为大气中水汽分压强，单位为hPa，我们可以看出在无线电波在对流层里面进行传播时受气压、温度及湿度的影响很大，因而，此时的折射率正是关于温度、气压、湿度的函数，所以我们第一步基础探测数据的处理就是要用探测的数据计算出海拔高度h处相对应的折射率N，计算步骤如下：

(1).海拔高度h的计算

h和地面上空的大气压强密切相关，对于同一个地方在不同时间，不同季节气压强都是不一样的，由于受到重力的影响，气压随高度呈现指数递减，令地面的大气压是P0，高度为h处的气压为P，则有

由于在现实环境中大气潮湿的，因而在工程应用里必须将这一因素考虑进去，对上述式子进行修正，由于，空气压强和高度联系紧密，所以，高度值就可按照其与气压之间的联系进行计算，

实际中经常见到的压高公式为：

h＝hs+∑Δhi i＝1,2,3……(3)

式(3-6)中，hs为当地地面的海拔高度，单位为m；式(3-7)中为当地的经纬度，单位为rad；分别为两大气层间的温度、气压、湿度的平均值，pi及pi+1分别为第i、i+1层大气压强，

(2).温度T的计算

由对流层折射率计算公式(1)可以看出里面使用的是绝对温度T(K)，我们日常用的温度一般是℃，它们两者的相互换算如下所示：

T＝273.16+t(5)

(3).湿度e的计算

大气折射率N的计算中使用的是水汽压e，其表示方法为：

另外，当空气温度为露点温度td时，水汽压e为：

一般如果一地区存在河流、湖泊则此区域的空气湿度通常很大；而在温度较高的沙漠地带空气的一般较小，从总体情况看的话，水汽压与高度呈现下式所示的关系：

e＝e0×10^h/6.3 (8)

在利用压高公式和对流层大气折射率公式计算得到各个站点每时对应的高度和折射率(原始探测数据为：温度t、气压p、和相对湿度RH)之后将各个探空站的数据文件按照月份进行整理，整理后文件格式为第i个探测站，第j个月的探测数据。

6.根据权利要求1所述的全国大气折射率剖面的预测方法，其特征在于：大气分段模型参数的计算为：对流层中大气的分段大气模型如下式所示：

其中的参数主要有地面折射率N0、海拔1千米处的折射率N1、海拔9千米处的折射率N9、近地面1千米的大气折射率垂直梯度G、地面以上1千米至海拔9千米大气层的指数衰减系数C1及海拔大于9千米时的大气层的指数衰减系数C9，

计算这些参数时首先将每次探测数据读入到一个数组中，该数组参数含的参量为：探测站编号，年份，月份，探测时间，数据层数m以及相应高度对应的大气折射率(hi，Ni)，如果某次探测数据层数m<6，则忽略此组数据，同时将此次日期对应的全部参数置零，方便后续计算，因为当探测次数小于6时数据太少，此时的参数的计算产生的误差较大，(1).N0、N1、N9的计算方法

N0是地面折射率，N1是地面以上海拔1千米处的折射率，N9是地面以上海拔9千米处的折射率，这三个参数我们都采用线性插值方法来求得，计算步骤如下：

①将实测数据(hi，Ni)读入到数组h(i)，N(i)中，i＝1,2,……m；

②从i＝1开始比较h(i)-1000-hs，当第一次出现h(i)≥1000+hs时，则令s＝i，

此时：

③从r＝s+1开始比较h(i)-9000，当第一次出现h(i)≥9000时，令t＝r，此时：

若h(m)≤9000，则令t＝m，此时：

(2).G的计算方法

G为地面以上1千米之间的大气折射率垂直梯度，单位1/km，从大气分段模型可知：

G＝N1-N0 (13)

(3).C1的计算方法

C1为地面1千米到海拔9km之间的折射率指数衰减衰减系数，单位1/km，同样从大气分段模型可知：

另外由于参数C9在各个区域里面的变化不太大，通常使用时直接采用年平均值造成的误差对结果影响不大，而且现有探空数据也无法延伸到对流层顶，因此一般采用现有的统计结果。

7.根据权利要求1所述的全国大气折射率剖面的预测方法，其特征在于：去野值处理选用的是拉依达准则，其基本思想是，它认为研究数据里仅包含随机误差，然后通过对研究数据进行运算获得其标准偏差，之后根据得到的数据确定一个区间，而不在此区间内的误差就被认定为野值，将其去除即可，结合本发明使用的数据特点，由于原始探测数据主要包括温度、气压、湿度等数据量比较大，野值剔除比较困难，准确率较低，所以先对基础探测数据进行简单处理，即计算N0、G和C1之后，再对整理后的N0，G，C1，进行野值剔除。

8.根据权利要求1所述的全国大气折射率剖面的预测方法，其特征在于：有效数据的整理为有探空数据的栅格中的数据经过去除野值之后得到的就是有效数据，即各个探测站按照各个月份整理好的大气分段模型的参数值，而进行插值的时候要求每次运算使用的各个站点的数据必须是在同一水平面且是同一时刻的，但是在全国范围内地势高低起伏很大，不同地区的海拔高度相差很多，所以为了使每次进行插值时使用的数据在同一水平面，首先我们对已有探空数据的栅格中的数据进行整理可以得到这些栅格中各个月份的n组大气折射率剖面分段模型(n为10年的探测次数去掉野值之后的剩余组数)，根据此分段模型我们可以计算出此栅格中各个月份任意高度上的大气折射率，然后，我们选取海拔高度从0～20000米，Δh＝500m利用这些栅格中的大剖面分段模型进行计算，这样我们就可以得到这些栅格中各个月份的n组(hi,Ni)的值(hi为海拔高度，Ni为相应海拔高度上对应的折射率)，此时这些栅格中相同的hi对应的为统一水平面的折射率值，最后将相同月份中同一时间的(hi,Ni)放到一起，作为后面的插值使用。

9.根据权利要求1所述的全国大气折射率剖面的预测方法，其特征在于：径向基函数插值使用的是其中的复二次函数，它的表达式为：

φ(r)＝(c²+r²)^β (15)

我们在做插值计算时已知的数据点为全国100多个探空站10年各个月份的n组(hi,Ni)i＝1,2,……m，其中m是hi从0间隔500到20000的个数，插值时构造如下函数：

并且使其满足插值条件：

在式(16)中：为了表示方便，令λi∈R,b1∈R,b2∈R,b2∈R,a∈R，是需要计算的参数；x,y是待插值点处的经纬度坐标；h是待插值点处的海拔高度；xi,yi参与插值计算的点的经纬度坐标；hi是参与插值计算的点的海拔高度；N是参与插值计算的点的个数；Ni是相应hi对应的折射率值，

同时其参数λi满足约束条件：

为了求解插值模型中的参数，我们定义一个N*N的矩阵Φ，Φij＝φ(||ri-rj||).i,j＝1,2,……N，如下所示：

令：α＝[b1,b2,b3,a]，

则加上约束条件后的插值方程组为：

即：

那么，λi,b1,b2,b3,a为：