本发明属于气象遥感技术领域，特别涉及一种反演大气水汽含量的方法及装置。
背景技术
水汽及其变化是天气和气候的主要驱动力，是预测降雨、中小尺度恶劣天气以及全球气候变化的一个非常重要的物理量。目前大气水汽含量的传统探测技术主要为无线电探空仪、气象卫星观测以及传统地基bds观测值反演。
无线电探空仪是由传感器、转换器和无线电发射机组成的测定自由大气温、压、湿等气象要素的仪器。该仪器搭载于探空气球上升到高空，用敏感元件直接测量大气压力、温度和相对湿度层结曲线，其机械的或电的输出由转换开关依次接入编码器转变成电信号，再由发射机经调幅或调频发送，接收端在地面进行接收、解调和记录。此方法费用昂贵、充足性差。
气象卫星所载各种气象遥感器，接收和测量地球及其大气层的可见光、红外和微波辐射，并将其转换成电信号传送给地面站。地面站将卫星传来的电信号复原，绘制成各种云层、地表和海面图片，再经进一步处理和计算，得出各种气象资料。但由于技术限制，该方法的数据源有限，也存在成本较高的问题，一般个人和机构无法有效利用。
传统的bds观测值反演大气水汽含量的方法，是根据模型反演后的精密星历以及观测站的气象观测数据获得相当于同样水汽含量的水柱高的大气可降水量(pwv)。由于bds精密星历为事后星历，存在很大的延迟性，因此无法应用于大气水汽含量的实时反演。同时，由于成本等原因，观测站气象观测数据大多是通过气压插值获得，传统的气压插值方法无法很好的处理bds站和探空站高差差异较大的情况。
技术实现要素：
本发明的目的在于提供一种反演大气水汽含量的方法及装置，用于解决现有技术中的水汽含量检测方法不能对水汽含量实时反演的问题。
为实现上述目的，本发明提供了一种反演大气水汽含量的方法，包括如下步骤：
1)采集观测站数据，并下载各个时间点的星历数据，从各个时间点数据里选取与观测时刻最近的星力数据作为反演的星历数据；
2)根据观测站数据计算观测站坐标信息，并计算观测站上空对流层天顶总延迟；获取观测站地面气压，根据观测站坐标信息及地面气压计算获得天顶干延迟，将对流层天顶总延迟减去所述天顶干延迟得到天顶湿延迟；
3)根据水汽的气体常数、大气折射常数及地区平均温度及所述天顶湿延迟计算大气降水量。
进一步地，利用分段高差的气压插值模型获取观测站地面气压。利用分段高差的气压插值模型，增加了地面大气数据的精度，减少了气压观测的成本。
进一步地，所述分段高差的气压插值模型表示为：
p＝pr(1-0.000022δh)5.225
其中，pr为观测站气压，△h为观测站高程与参考站高程的高差。
进一步地，所述天顶干延迟的表达式为：
f(λ，h)＝1-0.00266cos(2λ)-0.000281h
其中，zhd为天顶干延迟，p为观测站实测地面气压，λ为观测站纬度，h为观测站高程。
进一步地于，大气降水量的表达式为：
其中，pwv为大气降水量，ρ为液态水密度，k为转换系数，zwd为天顶湿延迟。
本发明还提供了一种反演大气水汽含量的装置，包括存储器、处理器以及存储在存储器上并可在处理器上运行时的计算机程序，所述处理器执行所述程序时实现以下步骤：
1)采集观测站数据，并下载各个时间点的星历数据，从各个时间点数据里选取与观测时刻最近的星力数据作为反演的星历数据；
2)根据观测站数据计算观测站坐标信息，并计算观测站上空对流层天顶总延迟；获取观测站地面气压，根据观测站坐标信息及地面气压计算获得天顶干延迟，将对流层天顶总延迟减去所述天顶干延迟得到天顶湿延迟；
3)根据水汽的气体常数、大气折射常数及地区平均温度及所述天顶湿延迟计算大气降水量。
进一步地，利用分段高差的气压插值模型获取观测站地面气压。利用分段高差的气压插值模型，增加了地面大气数据的精度，减少了气压观测的成本。
进一步地，所述分段高差的气压插值模型表示为：
p＝pr(1-0.000022δh)5.225
其中，pr为观测站气压，△h为观测站高程与参考站高程的高差。
进一步地，所述天顶干延迟的表达式为：
f(λ，h)＝1-0.00266cos(2λ)-0.000281h
其中，zhd为天顶干延迟，p为观测站实测地面气压，λ为观测站纬度，h为观测站高程。
进一步地，大气降水量的表达式为：
其中，pwv为大气降水量，ρ为液态水密度，k为转换系数，zwd为天顶湿延迟。
本发明的有益效果是：
本发明通过采集观测站数据，并下载各个时间点的星历数据，从各个时间点数据里选取与观测时刻最近的星力数据作为反演的星历数据；根据观测站数据计算观测站坐标信息，并计算观测站上空对流层天顶总延迟；获取观测站地面气压，根据观测站坐标信息及地面气压计算获得天顶干延迟，将对流层天顶总延迟减去天顶干延迟得到天顶湿延迟；根据水汽的气体常数、大气折射常数及地区平均温度及天顶湿延迟计算大气降水量。在星历数据满足反演大气水汽含量的精度要求的前提下，实现了观测实时反演大气水汽含量，提高了大气水汽含量测量的精度，且适用于观测站和参考站高差差异较大的情况。
附图说明
图1为igu预报星历动态选择方案示意图；
图2为本发明的bds近实时反演大气水汽含量的方法流程图；
图3为观测站上空对流层天顶总延迟计算流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明：
为了实现对大气水汽含量的实时预测，本发明的方法与中国北斗卫星导航系统bds结合起来，提供了一种bds近实时反演大气水汽含量的方法，如图2所示，该方法包括如下步骤：
步骤1、将bds天线、接收机安装至观测站上，并将bds接收机连接至服务器，系统连接电源启动。
步骤2、服务器实时采集bds测站的观测数据，并从igmas下载超快速星历数据igu，或者采用其他的一些手段获取星历数据。
步骤3、根据动态igmas超快速星历选择方案获得用于模型反演的星历数据，根据星历数据反演大气水汽含量的过程为，bds信号在穿越地球大气层时，要受到大气和电离层的折射影响，而当信号穿越大气层中对流层，由于大气成分，特别是水汽的折射会导致信号传播的速度减弱，和路径弯曲，造成时间上的延迟。研究发现，这种独特的bds信号时延近似正比于沿信号传播路径上的水汽总含量。而通过bds测站的观测数据，以及近实时的igs超快速星历可以精确的获取这种bds信号时延。
步骤4、利用gamit高精度数据处理软件处理相应数据，获得bds测站坐标信息，以及bds测站上空对流层天顶总延迟。
步骤5、利用双频接收机观测技术改正电离层延迟至毫米级。
步骤6、利用基于分段高差的气压插值模型得到bds测站实测地面气压数据。
步骤7、根据步骤3获得的bds测站坐标信息以及saastamoinen模型计算获得天顶干延迟。
步骤8、将步骤3获得的对流层天顶总延迟及天顶静力学延迟后获得天顶湿延迟。
步骤9、利用公式计算获得相当于同样水汽含量的水柱高的大气可降水量(pwv)。
其中，步骤3中的动态igmas超快速星历选择方案实现过程如下：
如图1所示，要实现bds近实时反演大气水汽含量，必须选择精度较高的igmas超快速星历igu。而igu每天03:00、09:00、15:00、21:00(世界时间)更新。且每次发布均包含48h的轨道和钟差信息，前24h是根据igmas跟踪站的实际观测值计算得到，后24h的轨道信息是外推预报的，在利用igu进行实时水汽反演时，主要使用的是外推部分星历，而外推星历的轨道和钟差精度会随着外推时间的增加而显著降低。因此本发明采用如下所示的igu动态选择方案选择igu星历，实现bds近实时反演大气水汽含量：当bds观测数据时段位于01:30—07:30之间时，采用03:00发布的igu，则01:30-03:00的数据处理时间最多会滞后1.5h，而03:00-07:30的数据可以得到实时处理，且星历外推时间为4.5h，随后的观测时段以此类推。如此既可以保证数据的近实时处理，又可以减少星历数据的外推时间，提高外推星历的轨道和钟差精度。
步骤4与步骤5中，由于bds信号穿越大气时因折射而产生的延迟主要包括电离层延迟、对流层延迟、接收机和卫星钟差引起的延迟与测量误差项等，通过双频接收机观测技术可以改正电离层延迟至毫米级，而利用多站点长短基线的差分测量可以消除接收机和卫星钟差以及测量误差项。采用上述技术手段，再通过配置gamit软件的主要固定解算设置以及约束方案，利用gamit软件实时处理bds测站的bds观测数据，可以解算获得测站的对流层天顶总延迟参数ztd。
步骤6的分段高差的气压插值模型如下：
通过在一定范围内特定的bds测站配置气象仪器采集气象数据，并将这些测站作为参考站，那么实测地面气压数据p可根据实测的参考站地面气压数据pr利用标准大气(sa)模型估计：
p＝pr(1-0.000022δh)5.225(1)
其中，δh＝h-hr(2)
上式中，δh为bds测站高程h与参考站高程hr的高差。
由于bds测站上空实际的大气状况与理想状态下的标准大气存在较大的差异，因此上述标准大气模型估计得到的实测地面地面气压数据存在较大的误差，尤其是bds测站与参考站高差较大时，该模型的误差更大，因此在实际工作中，采用的是根据上述标准大气模型，通过大量实验结果分析出的地面气压差值模型。在该模型中当δh，即bds测站高程与参考站高程的高差一定时，bds测站气压与参考站气压关系为：
p＝a×pr(3)
系数a与δh的关系为：
步骤7中saastamoinen模型计算获得天顶静力学延迟具体实现过程如下：
利用气压计精确获得地面气压p，则根据saastamoinen模型，干延迟精度可以达到mm量级：
f(λ，h)＝1-0.00266cos(2λ)-0.000281h(6)
zwd＝ztd-zhd(7)
其中，ztd为对流层天顶总延迟，zhd为天顶静力学延迟干延迟，zwd为天顶湿延迟，p为bds测站的地面大气压，λ为bds测站的纬度，h为bds测站的高程。
本实施例中，对流层天顶总延迟ztd的具体计算过程如图3所示。
步骤9中，利用公式计算获得大气可降水量(pwv)具体实现过程如下：
其中，转换系数k的计算公式如下：
地区加权平均温度t计算公式为：
其中，k′2，k3均为大气折射常数，ρ为液态水密度，rv为水汽的气体常数，tm为地区加权平均温度，pv为垂直分布上某点的水气分压，t为同一点的气温，hs为bds测站的高程。
本发明根据动态igmas超快速星历选择方案，获得近实时的卫星星历，并利用基于分段高差的气压插值模型得到bds测站实测地面气压数据，最终根据bds反演pmv的原理获得大气可降水量。本发明的方法适用于bds测站和参考站高差差异较大的情况，且提高了大气水汽含量的检测精度，bds测站数据反演的一天时间内，有6h的延迟小于1.5h，18h为实时预报。
以上给出了具体的实施方式，但本发明不局限于以上所描述的实施方式。本发明的基本思路在于上述基本方案，对本领域普通技术人员而言，根据本发明的教导，设计出各种变形的模型、公式、参数并不需要花费创造性劳动。在不脱离本发明的原理和精神的情况下对实施方式进行的变化、修改、替换和变形仍落入本发明的保护范围内。