一种反演大气水汽含量的方法及装置制造方法及图纸

本发明专利技术提供了一种反演大气水汽含量的方法及装置，通过采集BDS观测数据，并从IGMAS下载星历数据IGU，根据星历数据IGU获取用于反演的星历数据；根据上述信息获取观测站坐标信息，及观测站上空对流层天顶总延迟，根据观测站坐标信息及saastamoinen模型计算获得天顶干延迟，将对流层天顶总延迟减去所顶干延迟得到天顶湿延迟；根据转换系数、液态水密度及天顶湿延迟计算大气降水量。在星历数据满足BDS反演大气水汽含量的精度要求的前提下，实现了BDS观测实时反演大气水汽含量，提高了大气水汽含量测量的精度，且适用于观测站和参考站高差差异较大的情况。

【技术实现步骤摘要】
一种反演大气水汽含量的方法及装置
本专利技术属于气象遥感
，特别涉及一种反演大气水汽含量的方法及装置。
技术介绍
水汽及其变化是天气和气候的主要驱动力，是预测降雨、中小尺度恶劣天气以及全球气候变化的一个非常重要的物理量。目前大气水汽含量的传统探测技术主要为无线电探空仪、气象卫星观测以及传统地基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的轨道信息是外推预报的，在利用本文档来自技高网...

【技术保护点】
1.一种反演大气水汽含量的方法，其特征在于，包括如下步骤：1)采集观测站数据，并下载各个时间点的星历数据，从各个时间点数据里选取与观测时刻最近的星力数据作为反演的星历数据；2)根据观测站数据计算观测站坐标信息，并计算观测站上空对流层天顶总延迟；获取观测站地面气压，根据观测站坐标信息及地面气压计算获得天顶干延迟，将对流层天顶总延迟减去所述天顶干延迟得到天顶湿延迟；3)根据水汽的气体常数、大气折射常数及地区平均温度及所述天顶湿延迟计算大气降水量。

【技术特征摘要】
1.一种反演大气水汽含量的方法，其特征在于，包括如下步骤：1)采集观测站数据，并下载各个时间点的星历数据，从各个时间点数据里选取与观测时刻最近的星力数据作为反演的星历数据；2)根据观测站数据计算观测站坐标信息，并计算观测站上空对流层天顶总延迟；获取观测站地面气压，根据观测站坐标信息及地面气压计算获得天顶干延迟，将对流层天顶总延迟减去所述天顶干延迟得到天顶湿延迟；3)根据水汽的气体常数、大气折射常数及地区平均温度及所述天顶湿延迟计算大气降水量。2.根据权利要求1所述的反演大气水汽含量的方法，其特征在于，利用分段高差的气压插值模型获取观测站地面气压。3.根据权利要求2所述的反演大气水汽含量的方法，其特征在于，所述分段高差的气压插值模型表示为：P＝Pr(1-0.000022ΔH)5.225其中，Pr为观测站气压，△H为观测站高程与参考站高程的高差。4.根据权利要求3所述的反演大气水汽含量的方法，其特征在于，所述天顶干延迟的表达式为：f(λ，H)＝1-0.00266cos(2λ)-0.000281H其中，ZHD为天顶干延迟，P为观测站实测地面气压，λ为观测站纬度，H为观测站高程。5.根据权利要求1所述的反演大气水汽含量的方法，其特征在于，大气降水量的表达式为：其中，PWV为大气降水量，ρ为液态水密度，K为转换系数，ZWD为天顶湿延迟。6.一种反演大气水汽含量的装置，包括存储器、处理...

【专利技术属性】
技术研发人员：陈明剑，闫欢欢，李俊毅，李万里，李广云，
申请(专利权)人：中国人民解放军战略支援部队信息工程大学，
类型：发明
国别省市：河南,41