CN109001382B - 一种基于cors的区域大气水汽实时监测方法及系统 - Google Patents

一种基于cors的区域大气水汽实时监测方法及系统 Download PDF

Info

Publication number
CN109001382B
CN109001382B CN201811098226.4A CN201811098226A CN109001382B CN 109001382 B CN109001382 B CN 109001382B CN 201811098226 A CN201811098226 A CN 201811098226A CN 109001382 B CN109001382 B CN 109001382B
Authority
CN
China
Prior art keywords
real
time
water vapor
precise
atmospheric
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811098226.4A
Other languages
English (en)
Other versions
CN109001382A (zh
Inventor
姚宜斌
孔建
刘邢巍
许超钤
张良
彭文杰
张泽烈
刘磊
蒲德祥
吴国梁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing Geographic Information And Remote Sensing Application Center (chongqing Surveying And Mapping Product Quality Inspection And Testing Center)
Wuhan University WHU
Original Assignee
Chongqing Geographic Information And Remote Sensing Application Center (chongqing Surveying And Mapping Product Quality Inspection And Testing Center)
Wuhan University WHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing Geographic Information And Remote Sensing Application Center (chongqing Surveying And Mapping Product Quality Inspection And Testing Center), Wuhan University WHU filed Critical Chongqing Geographic Information And Remote Sensing Application Center (chongqing Surveying And Mapping Product Quality Inspection And Testing Center)
Priority to CN201811098226.4A priority Critical patent/CN109001382B/zh
Publication of CN109001382A publication Critical patent/CN109001382A/zh
Application granted granted Critical
Publication of CN109001382B publication Critical patent/CN109001382B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004Gaseous mixtures, e.g. polluted air
    • G01N33/0009General constructional details of gas analysers, e.g. portable test equipment
    • G01N33/0027General constructional details of gas analysers, e.g. portable test equipment concerning the detector
    • G01N33/0036Specially adapted to detect a particular component

Abstract

本发明公开了一种基于CORS的区域大气水汽实时监测方法及系统,其中方法包括步骤:将各CORS站的观测值实时汇集到数据中心形成GNSS数据;获取IGS精密轨道和精密钟差实时改正产品,并读入GNSS数据,采用精密单点定位技术估测各CORS站上空的大气总延迟;利用全球加权平均温度模型将大气总延迟中的湿延迟转换为大气水汽含量,得到所有CORS站天顶上空的大气水汽含量;采用克里金插值法将所述大气水汽含量插值获得特定时间分辨率与特定空间分辨率的区域大气水汽含量;利用区域大气水汽含量,实现区域上空水汽含量的实时监测。其显著效果是:实现了实时、稳定、高精度、高时空分辨率、全天候和全天时的水汽监测。

Description

一种基于CORS的区域大气水汽实时监测方法及系统
技术领域
本发明涉及到大气水汽探测技术领域,具体涉及一种基于CORS的区域大气水汽实时监测方法。
背景技术
水汽是地球大气中的微量气体,在大气中所占比例很小,仅0.1%~3%,却是大气中最活跃的组分,水汽相比于其他微量气体更具有重要性。天气现象大多是大气水汽变化的结果,水汽在相变过程中吸收和释放大量潜热,直接影响地面和空气温度,进而影响大气垂直稳定度和对流天气系统的形成。大气水汽的三维分布、水汽垂直输送和相变是制约中尺度天气系统发展的动力机制之一。
目前,常规大气水汽探测手段有无线电探空、微波辐射计和卫星遥感。然而这些探测手段有如下不足:无线电探空仪不能重新利用,成本高昂;同时气象站一般会采用每6小时或12小时进行一次探空观测,其时间分辨率太低,无法做到实时大气水汽监测;微波辐射计是被动地接收各个高度传来的温度辐射的微波信号来测量大气参数,设备成本更高,其单台设备的价格达到了数十万元人民币;卫星遥感方式中由于卫星重访同一地点的时间间隔较长,以MODIS为例,其卫星的重访周期为一天四次,远远达不到实时水汽监测的目的,此外,MODIS进行水汽探测还受到云雨的影响。
随着GNSS应用领域的不断发展,GNSS气象学成为一种较新的水汽探测方式,利用放置在地面上的GNSS接收机,接收测量GNSS卫星的信号纵向穿过大气层到达地面所引起的延迟量,进而反演出天顶方向整层大气层或信号斜路径上的大气水汽的累积量,具有精度高、时空分辨率高、连续、不受云雨影响、成本低等特点,是一种新型廉价的大气水汽探测方法。
常用的基于GNSS进行大气水汽探测方法是精密单点定位(Precision PointPositioning,PPP)技术。该方法利用精密轨道和精密钟差,解算GNSS测站的天顶大气湿延迟,再利用加权平均温度将大气湿延迟转换为大气水汽。国际GNSS组织向全球提供精密星历和钟差产品,成为精密单点定位的必要条件。但是,事后的精密星历和钟差产品,如最终星历或快速星历存在十几天至几天的延迟,只能进行事后大气水汽的计算,无法满足实时大气水汽的监测。超快星历延迟3小时发布,并且预报的轨道精度满足精密单点定位的要求,但钟差精度有限,无法进行毫米级大气水汽监测。
综上而言,虽然各种基于GNSS的大气水汽探测方法被提出,但目前基于GNSS尚无真正实时的低成本廉价的大气水汽监测手段。
发明内容
针对现有技术的不足,本发明的目的是提供一种基于CORS的区域大气水汽实时监测方法及系统,通过收集CORS站观测数据和IGS精密轨道和钟差实时改正产品,并利用精密单点定位技术估计测站天顶方向的对流层延迟,最后转换得到区域大气水汽含量产品,达到对区域上空水汽含量的实时监测。
为达到上述目的,本发明采用的技术方案如下:
一种基于CORS的区域大气水汽实时监测方法,其关键在于,包括以下步骤:
步骤1、将GNSS系统中各CORS站的观测值,以特定频率实时汇集到数据中心形成GNSS数据,并存入观测值数据库;
步骤2、获取IGS精密轨道和精密钟差实时改正产品,并读入GNSS数据,采用精密单点定位技术估测各CORS站上空的大气总延迟;
步骤3、利用全球加权平均温度模型将大气总延迟中的湿延迟转换为大气水汽含量,得到所有CORS站天顶上空的大气水汽含量;
步骤4、采用克里金插值法将所述大气水汽含量进行插值处理,获得特定时间分辨率与特定空间分辨率的区域大气水汽含量;
步骤5、利用区域大气水汽含量,实现区域上空水汽含量的实时监测。
进一步的,步骤1中所述特定频率为1HZ。
进一步的,所述IGS精密轨道和精密钟差实时改正产品的获取过程为:将IGS实时服务产品提供的精密轨道改正数与精密钟差改正数,应用于广播星历的轨道和钟差,生成IGS精密轨道和精密钟差实时改正产品,并存入精密轨道和钟差数据库。
进一步的,所述精密单点定位技术进行大气总延迟估测时,所采用的数据为每隔时间段T1从观测值数据库提取的包含最临近时间段T2数据的RINEX文件,以及每隔时间段T3从精密轨道和钟差数据库提取的包含最临近时间段T4数据的精密轨道和精密钟差文件。
进一步的,所述时间段T1与T3的值均为5min,所述最临近时间段T2与T4的值均为2h。
进一步的,所述全球加权平均温度模型采用GTm-III模型,公式如下:
PWV=Π×ZWD,
其中,PWV为大气水汽含量,
Figure BDA0001806047620000041
为转换因子,ρw为液态水密度,Rv为水汽气体常数,k'2、k3为大气折射常数,Tm为加权平均温度,ZWD=ZTD-ZHD为大气湿延迟,ZTD为大气天顶总延迟,ZHD为静力学延迟。
进一步的,所述克里金插值法的数学表达式为:
Figure BDA0001806047620000042
其中,Z(x0)为未采样点x0处的大气水汽含量,Z(xi)为采样点xi处的大气水汽含量,λi为权系数,i=1,2,…,n。
根据上述的基于CORS的区域大气水汽实时监测方法,本发明还提出了一种区域大气水汽实时监测系统,其关键在于:包括CORS数据实时采集模块、RTS数据实时采集模块、实时精密星历生成模块、实时数据管理模块与数据处理模块;
所述CORS数据实时采集模块,用于以每秒的时间间隔接收RTCM格式的CORS站观测数据;
所述RTS数据实时采集模块,用于实时接收RTS轨道改正数和RTS钟差改正数,并将数据传入实时精密星历生成模块;
所述实时精密星历生成模块,用于将RTS轨道改正数和RTS钟差改正数应用于广播星历,生成精密轨道和精密钟差实时改正产品;
所述实时数据管理模块,用于以数据库的形式管理CORS站观测数据、精密轨道和精密钟差实时改正产品,并将CORS站观测数据每隔时间段T1生成一次包含最临近时间段T2内数据的RINEX文件,以及将精密轨道和精密钟差产品每隔时间段T3生成一次包含最临近时间段T4内数据的精密轨道和精密钟差文件;
所述数据处理模块,用于利用精密单点定位技术对RINEX文件、精密轨道和精密钟差产品进行精密单点定位,生成各CORS站上空的大气总延迟,并分离出天顶湿延迟;利用全球加权平均温度模型将大气总延迟中的湿延迟转换为各CORS站上空的大气水汽含量;采用克里金差值法插值得到特定时间分辨率与特定空间分辨率的区域大气水汽含量数据;依据区域大气水汽含量数据实现区域上空水汽含量的实时监测。
所述的基于CORS的区域大气水汽实时监测系统还包括显示模块和产品服务模块,所述显示模块用于实时显示区域大气水汽含量随时间的变化情况;所述产品服务模块用于查询、统计和下载大气水汽含量数据。
进一步的,所述时间段T1与T3的值均为5min,所述最临近时间段T2与T4的值均为2h。
显著效果是:
本发明通过收集GNSS站观测数据和IGS精密轨道和钟差产品,并利用精密单点定位技术估计测站天顶方向的对流层延迟,最后利用区域大气水汽含量产品,达到对区域上空水汽含量的实时监测,从而实现了实时稳定地提供基于GPS、GLONASS和北斗系统的高时空分辨率区域对流层延迟产品,以及毫米级精度的大气水汽产品,为改善短期临近的精细化气象预报提供了高精度可靠的数据依据;
区域CORS为服务于区域经济建设的基础设施,一经建设完成,即可长期无间断跟踪观测GNSS卫星,而基于CORS的实时水汽产品,是CORS的增值服务,在不增加建设成本的情况下,利用该技术能够提供稳定、高精度、实时、高分辨率、全天候和全天时的水汽监测,在有限的费用成本的前提下,是任何其他水汽探测技术都无法做到的。
附图说明
图1是本发明的方法流程图;
图2是本发明的原理框图。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
如图1所示,一种基于CORS的区域大气水汽实时监测方法,其处理步骤如下:
步骤1、将GNSS系统中第1~第n个CORS站的观测值,以1HZ的特定频率实时汇集到数据中心形成GNSS数据,并存入观测值数据库;
步骤2、每间隔5秒或1分钟接收RTS精密轨道改正数,同时每5秒或10接收RTS精密钟差改正数,并将获得的精密轨道改正数和精密钟差改正数应用于广播星历的轨道和钟差,生成IGS精密轨道和精密钟差实时改正产品,同时存入精密轨道和钟差数据库;
然后,每隔特定时间T1(本例T1=5min)从观测值数据库中取出最临近时间段T2(本例T2=2小时)的观测值生成为RINEX文件,同时每隔时间段T3(本例T3=5min)从精密轨道和钟差数据库提取包含最临近时间段T4(本例T4=2小时)的精密轨道和精密钟差文件,生成待处理任务;
之后,采用精密单点定位技术,每秒等待解算任务,当有新的解算任务出现时,及时解算生成所有CORS站上空的大气总延迟;
在精密单点定位中,主要通过处理测码伪距P和载波相位观测值L来得到定位结果,原始观测如下:
Figure BDA0001806047620000071
Figure BDA0001806047620000072
其中,i=1,2,表示观测值频率f1,f2,本文采用双频观测数据,f1=1575.42MHZ,f2=1227.60MHZ;ρ表示卫星与测站之间的实际几何距离;dtrcv和dtsat分别表示接收机(rcv)和卫星(sat)的钟差;Tr为对流层延迟
Figure BDA0001806047620000073
为电离层延迟系数,与频率相关,αiI为电离层延迟;
Figure BDA0001806047620000074
Figure BDA0001806047620000075
分别为接收机和卫星在fi上的伪距硬件延迟;
Figure BDA0001806047620000081
Figure BDA0001806047620000082
分别为接收机和卫星在fi上的相位硬件延迟;λi为载波相位Li的波长;Ni为载波相位Li的整周模糊度项;w为天线相位缠绕;M和m分别为伪距和相位的多路径误差;εPi和εLi分别为伪距和相位噪声。
需要说明的是,由于精密单点定位需要一段时间的收敛过程,因此在进行实时精密单点定位时,需要累积一段时间的观测数据和精密轨道、精密钟差产品。也即是选用最相近时间段的数据进行处理的原因所在。
步骤3、利用全球加权平均温度模型将大气总延迟中的湿延迟转换为大气水汽含量,得到所有CORS站天顶上空的大气水汽含量;
当时间分辨率要求更高时,适当将生成RINEX文件的频率提高,即可得到更高时间分辨率的区域大气水汽含量产品。
全球大气加权平均温度系列模型包括GTm-I、GTm-II和GTm-III三个模型。优选的,本发明中所述全球加权平均温度模型采用GTm-III模型,大气加权平均温度的是将大气湿延迟转换为大气水汽的必要参数。大气天顶总延迟、大气湿延迟与大气水汽的关系可以表达为:
PWV=Π×ZWD,
其中,PWV为大气水汽含量,
Figure BDA0001806047620000083
为转换因子,ρw为液态水密度,Rv为水汽气体常数,k'2、k3为大气折射常数,Tm为加权平均温度,ZWD=ZTD-ZHD为大气湿延迟,ZTD为大气天顶总延迟,ZHD为静力学延迟。
Π为转换因子,是大气加权平均温度Tm的函数,其他参数为经验常数。由于Tm的计算需要采用长期累积的气象资料计算得到,在没有气象资料的地区无法计算得到该参数。而所述全球大气加权平均温度模型,则提供了球谐函数表达的全球范围内的大气加权平均温度。只要输入待计算数据点的经纬度和年积日,即可得到大气加权平均温度Tm。进而可以计算转换因子Π,并根据解算出的天顶总延迟提取出大气湿延迟,即可将大气湿延迟转换为大气水汽含量。
步骤4、采用普通克里金插值法将所述大气水汽含量进行插值处理,获得特定时间分辨率与特定空间分辨率的区域大气水汽含量;所述普通克里金插值法的数学表达式为:
Figure BDA0001806047620000091
其中,Z(x0)为未采样点x0处的大气水汽含量,Z(xi)为采样点xi处的大气水汽含量,λi为权系数,i=1,2,…,n。其中,权系数λi的求解方程组为
Figure BDA0001806047620000092
其中,j=1,2,…,nμ为拉格朗日乘子。
步骤5、利用区域大气水汽含量,实现区域上空水汽含量的实时监测。
在具体实施时,利用精密单点定位技术进行大气水汽监测的精度取决于精密单点定位的精度。目前,精密单点定位收敛后的静态定位精度可达厘米级,相应的天顶总延迟含量的精度也可达厘米级。将天顶总延迟分离出湿延迟后,转换成大气水汽含量的精度在2-3mm。
如图2所示,一种基于CORS的区域大气水汽实时监测系统,包括CORS数据实时采集模块、RTS数据实时采集模块、实时精密星历生成模块、实时数据管理模块、数据处理模块、显示模块和产品服务模块,其中:
所述CORS数据实时采集模块,用于以每秒的时间间隔接收RTCM格式的CORS站观测数据;
所述RTS数据实时采集模块,用于实时接收RTS轨道改正数和RTS钟差改正数,并将数据传入实时精密星历生成模块;
所述实时精密星历生成模块,用于将RTS轨道改正数和RTS钟差改正数应用于广播星历,生成精密轨道和精密钟差实时改正产品;
所述实时数据管理模块,用于以数据库的形式管理CORS站观测数据、精密轨道和精密钟差实时改正产品,并将CORS站观测数据每隔时间段T1生成一次包含最临近时间段T2内数据的RINEX文件,以及将精密轨道和精密钟差产品每隔时间段T3生成一次包含最临近时间段T4内数据的精密轨道和精密钟差文件;
所述数据处理模块,用于利用精密单点定位技术对RINEX文件、精密轨道和精密钟差实时改正产品进行精密单点定位,生成各CORS站上空的大气总延迟,并分离出天顶湿延迟;利用全球加权平均温度模型将大气总延迟中的湿延迟转换为各CORS站上空的大气水汽含量;采用克里金差值法插值得到特定时间分辨率与特定空间分辨率的区域大气水汽含量数据;依据区域大气水汽含量数据实现区域上空水汽含量的实时监测;
所述显示模块用于实时显示区域大气水汽含量随时间的变化情况;
所述产品服务模块用于查询、统计和下载大气水汽含量数据。
本发明相较于传统的大气水汽检测,有如下优势:
(1)基于实时服务产品(RTS)的实时大气水汽监测方法
目前,尚无真正实时的GNSS水汽探测手段。见诸报道的利用GNSS进行水汽探测方法,主要是通过事后精密星历进行处理,并在水汽反演的可行性和精度上做出论证。而本发明采用实时服务产品(RTS)结合精密单点定位技术反演大气水汽含量,并提出行之有效的基于RTS的实时大气水汽监测系统结构。虽然从理论上,相对定位模型也可以实现实时大气水汽监测,但基于RTS和精密单点定位技术在实现上较相对定位模式更容易。
(2)基于实时CORS的区域大气水汽监测方法
当前见诸报道的GNSS大气水汽探测方法均是通过事后存储的观测文件进行。只有基于实时数据流的区域大气水汽反演,才能真正做到实时水汽监测。本发明基于实时CORS网络,实时接收RTCM3数据流格式并解码,存入观测值数据库,进而进行大气水汽解算,实现了数据的实时性,数据的时效性更好。
(3)可用于事后精密单点定位程序的实时处理模式
目前大多数的精密单点定位程序均是基于文件模式的事后处理方式,为了使这种精密单点定位程序可以用于实时精密单点定位的处理流程中来,本发明提出将实时数据流存入观测值数据库和精密轨道和钟差数据库,并每隔特定时间间隔生成一定时长的观测文件以及精密轨道和钟差文件,这样的处理模式既能满足实时性要求,又可以不改动已有精密单点定位程序。
(4)基于全球大气加权平均温度的实时大气水汽转换方法
将精密单点定位得到的大气总延迟分离出湿延迟后,可通过大气加权平均温度转换为大气水汽含量。然而,常用的大气加权平均温度模型需要长期气象观测资料计算求出,对于无实测气象资料的地区,无法将湿延迟有效转换为大气水汽含量,也就无法进行实时水汽监测。本发明采用全球大气加权平均温度模型进行转换,该模型是利用了全球气象观测资料计算出的经验模型,可以在全球任何地区,包括没有大气观测资料的地区使用。

Claims (8)

1.一种基于CORS的区域大气水汽实时监测方法,其特征在于,包括以下步骤:
步骤1、将GNSS系统中各CORS站的观测值,以特定频率实时汇集到数据中心形成GNSS数据,并存入观测值数据库;
步骤2、获取IGS精密轨道和精密钟差实时改正产品,并读入GNSS数据,采用精密单点定位技术估测各CORS站上空的大气总延迟;
所述IGS精密轨道和精密钟差实时改正产品的获取过程为:将IGS实时服务产品提供的精密轨道改正数与精密钟差改正数,应用于广播星历的轨道和钟差,生成IGS精密轨道和精密钟差实时改正产品,并存入精密轨道和钟差数据库;
所述精密单点定位技术进行大气总延迟估测时,所采用的数据为每隔时间段T1从观测值数据库提取的包含最临近时间段T2数据的RINEX文件,以及每隔时间段T3从精密轨道和钟差数据库提取的包含最临近时间段T4数据的精密轨道和精密钟差文件;
采用精密单点定位技术,每秒等待解算任务,当有新的解算任务出现时,及时解算生成所有CORS站上空的大气总延迟,观测方程为:
Figure FDA0002326597970000011
Figure FDA0002326597970000012
其中,Pi表示观测值频率i上的测码伪距,Li表示观测值频率i上的载波相位观测值,i=1,2表示观测值频率f1,f2,ρ表示卫星与测站之间的实际几何距离;dtrcv和dtsat分别表示接收机和卫星的钟差;Tr为对流层延迟,
Figure FDA0002326597970000021
为电离层延迟系数,与频率相关,αiI为电离层延迟;
Figure FDA0002326597970000022
Figure FDA0002326597970000023
分别为接收机和卫星在fi上的伪距硬件延迟;
Figure FDA0002326597970000024
Figure FDA0002326597970000025
分别为接收机和卫星在fi上的相位硬件延迟;λi为载波相位Li的波长;Ni为载波相位Li的整周模糊度项;w为天线相位缠绕;M和m分别为伪距和相位的多路径误差;εPi和εLi分别为伪距和相位噪声;
步骤3、利用全球加权平均温度模型将大气总延迟中的湿延迟转换为大气水汽含量,得到所有CORS站天顶上空的大气水汽含量;
步骤4、采用克里金插值法将所述大气水汽含量进行插值处理,获得特定时间分辨率与特定空间分辨率的区域大气水汽含量;
步骤5、利用区域大气水汽含量,实现区域上空水汽含量的实时监测。
2.根据权利要求1所述的基于CORS的区域大气水汽实时监测方法,其特征在于:步骤1中所述特定频率为1HZ。
3.根据权利要求1所述的基于CORS的区域大气水汽实时监测方法,其特征在于:所述时间段T1与T3的值均为5min,所述最临近时间段T2与T4的值均为2h。
4.根据权利要求1所述的基于CORS的区域大气水汽实时监测方法,其特征在于:所述全球加权平均温度模型采用GTm-III模型,公式如下:
PWV=Π×ZWD,
其中,PWV为大气水汽含量,
Figure FDA0002326597970000031
为转换因子,ρw为液态水密度,Rv为水汽气体常数,k'2、k3为大气折射常数,Tm为加权平均温度,ZWD=ZTD-ZHD为大气湿延迟,ZTD为大气天顶总延迟,ZHD为静力学延迟。
5.根据权利要求1所述的基于CORS的区域大气水汽实时监测方法,其特征在于:所述克里金插值法的数学表达式为:
Figure FDA0002326597970000032
其中,Z(x0)为未采样点x0处的大气水汽含量,Z(xi)为采样点xi处的大气水汽含量,λi为权系数,i=1,2,…,n。
6.一种基于CORS的区域大气水汽实时监测系统,其特征在于:包括CORS数据实时采集模块、RTS数据实时采集模块、实时精密星历生成模块、实时数据管理模块与数据处理模块;
所述CORS数据实时采集模块,用于以每秒的时间间隔接收RTCM格式的CORS站观测数据;
所述RTS数据实时采集模块,用于实时接收RTS轨道改正数和RTS钟差改正数,并将数据传入实时精密星历生成模块;
所述实时精密星历生成模块,用于将RTS轨道改正数和RTS钟差改正数应用于广播星历,生成精密轨道和精密钟差实时改正产品;
所述实时数据管理模块,用于以数据库的形式管理CORS站观测数据、精密轨道和精密钟差实时改正产品,并将CORS站观测数据每隔时间段T1生成一次包含最临近时间段T2内数据的RINEX文件,以及将精密轨道和精密钟差产品每隔时间段T3生成一次包含最临近时间段T4内数据的精密轨道和精密钟差文件;
所述数据处理模块,用于利用精密单点定位技术对RINEX文件、精密轨道和精密钟差产品进行精密单点定位,生成各CORS站上空的大气总延迟,并分离出天顶湿延迟;利用全球加权平均温度模型将大气总延迟中的湿延迟转换为各CORS站上空的大气水汽含量;采用克里金差值法插值得到特定时间分辨率与特定空间分辨率的区域大气水汽含量数据;依据区域大气水汽含量数据实现区域上空水汽含量的实时监测;
所述精密单点定位技术解算生成所有CORS站上空的大气总延迟时的观测方程为:
Figure FDA0002326597970000041
Figure FDA0002326597970000042
其中,Pi表示观测值频率i上的测码伪距,Li表示观测值频率i上的载波相位观测值,i=1,2表示观测值频率f1,f2,ρ表示卫星与测站之间的实际几何距离;dtrcv和dtsat分别表示接收机和卫星的钟差;Tr为对流层延迟,
Figure FDA0002326597970000043
为电离层延迟系数,与频率相关,αiI为电离层延迟;
Figure FDA0002326597970000044
Figure FDA0002326597970000045
分别为接收机和卫星在fi上的伪距硬件延迟;
Figure FDA0002326597970000046
Figure FDA0002326597970000047
分别为接收机和卫星在fi上的相位硬件延迟;λi为载波相位Li的波长;Ni为载波相位Li的整周模糊度项;w为天线相位缠绕;M和m分别为伪距和相位的多路径误差;εPi和εLi分别为伪距和相位噪声。
7.根据权利要求6所述的基于CORS的区域大气水汽实时监测系统,其特征在于:还包括显示模块和产品服务模块,所述显示模块用于实时显示区域大气水汽含量随时间的变化情况;所述产品服务模块用于查询、统计和下载大气水汽含量数据。
8.根据权利要求6所述的基于CORS的区域大气水汽实时监测系统,其特征在于:所述时间段T1与T3的值均为5min,所述最临近时间段T2与T4的值均为2h。
CN201811098226.4A 2018-09-20 2018-09-20 一种基于cors的区域大气水汽实时监测方法及系统 Active CN109001382B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811098226.4A CN109001382B (zh) 2018-09-20 2018-09-20 一种基于cors的区域大气水汽实时监测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811098226.4A CN109001382B (zh) 2018-09-20 2018-09-20 一种基于cors的区域大气水汽实时监测方法及系统

Publications (2)

Publication Number Publication Date
CN109001382A CN109001382A (zh) 2018-12-14
CN109001382B true CN109001382B (zh) 2020-05-29

Family

ID=64592491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811098226.4A Active CN109001382B (zh) 2018-09-20 2018-09-20 一种基于cors的区域大气水汽实时监测方法及系统

Country Status (1)

Country Link
CN (1) CN109001382B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110031877B (zh) * 2019-04-17 2020-05-26 山东科技大学 一种基于grnn模型的区域nwp对流层延迟改正方法
JP7280947B2 (ja) * 2019-05-13 2023-05-24 古野電気株式会社 水蒸気観測システム及び水蒸気観測方法
CN110244387B (zh) * 2019-07-30 2021-06-22 成都润联科技开发有限公司 一种基于大气可降水量预测降雨天气的方法、装置、设备以及存储介质
CN110580453B (zh) * 2019-08-20 2020-06-30 清华大学 大气水汽含量监测方法、装置、计算机设备和存储介质
CN110531395B (zh) * 2019-09-05 2021-10-01 北京百度网讯科技有限公司 用于无人车定位的方法、装置和设备
CN110995381A (zh) * 2019-11-26 2020-04-10 天津大学 一种短波无线传输可用频率预测方法
CN111474606B (zh) * 2020-03-30 2022-04-05 中国人民解放军国防科技大学 基于双频双极化微波链路的降雨与水汽综合测量装置及方法
CN114019585A (zh) * 2021-10-11 2022-02-08 武汉大学 一种大高差地区高精度定位cors网fkp解算方法
CN114019584A (zh) * 2021-10-11 2022-02-08 武汉大学 一种大高差地区高精度cors网vrs解算方法
CN115144878B (zh) * 2022-07-13 2023-05-30 武汉大学 基于ppp的短距离大高差nrtk对流层延迟改正方法
CN115099159B (zh) * 2022-07-20 2023-03-07 武汉大学 基于神经网络且顾及地表差异的modis水汽反演方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007096466A1 (en) * 2006-02-27 2007-08-30 Antti Aarne Llmari Lange A method for calibrating teh carrier-phases of radio signals from satellites and other transmitters by using fast kalman filtering
KR101221931B1 (ko) * 2012-10-17 2013-01-15 한국해양과학기술원 위성신호 미약 환경에서의 관성센서를 이용한 선박의 위성측정치 생성방법 및 장치

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101526516B (zh) * 2009-04-10 2013-05-01 中国科学院测量与地球物理研究所 利用超导重力连续观测确定大气水汽垂直分布的方法
CN104020252A (zh) * 2014-05-15 2014-09-03 中国农业科学院农业资源与农业区划研究所 全球和区域水汽变化监测方法
CN107389029B (zh) * 2017-08-24 2019-10-29 北京市水文地质工程地质大队 一种基于多源监测技术融合的地面沉降集成监测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007096466A1 (en) * 2006-02-27 2007-08-30 Antti Aarne Llmari Lange A method for calibrating teh carrier-phases of radio signals from satellites and other transmitters by using fast kalman filtering
KR101221931B1 (ko) * 2012-10-17 2013-01-15 한국해양과학기술원 위성신호 미약 환경에서의 관성센서를 이용한 선박의 위성측정치 생성방법 및 장치

Also Published As

Publication number Publication date
CN109001382A (zh) 2018-12-14

Similar Documents

Publication Publication Date Title
CN109001382B (zh) 一种基于cors的区域大气水汽实时监测方法及系统
Haase et al. Accuracy and variability of GPS tropospheric delay measurements of water vapor in the western Mediterranean
Lee et al. Monitoring precipitable water vapor in real-time using global navigation satellite systems
Lu et al. Estimation and evaluation of real-time precipitable water vapor from GLONASS and GPS
Jin et al. Variability and climatology of PWV from global 13-year GPS observations
US6067852A (en) Method and apparatus using slant-path water delay estimates to correct global positioning satellite survey error
Sun et al. The FengYun-3C radio occultation sounder GNOS: a review of the mission and its early results and science applications
CN101893714B (zh) 全球卫星导航系统广播电离层时延修正方法
Deblonde et al. Evaluation of GPS precipitable water over Canada and the IGS network
Bianchi et al. Multi-year GNSS monitoring of atmospheric IWV over Central and South America for climate studies
Shoji A study of near real-time water vapor analysis using a nationwide dense GPS network of Japan
Wu et al. Validating HY-2A CMR precipitable water vapor using ground-based and shipborne GNSS observations
Webb et al. Kinematic GNSS estimation of zenith wet delay over a range of altitudes
Lu et al. Real-time retrieval of precipitable water vapor from Galileo observations by using the MGEX network
Liu et al. Evaluation of HY-2A satellite-borne water vapor radiometer with shipborne GPS and GLONASS observations over the Indian Ocean
Satirapod et al. Performance assessment of GPS-sensed precipitable water vapor using IGS ultra-rapid orbits: a preliminary study in Thailand
Gurbuz et al. Assessment of precipitable water vapor over Turkey using GLONASS and GPS
CN116068670B (zh) 一种无地面网络区适用的北斗水汽场实时重构方法及装置
Liu et al. The impact of different mapping function models and meteorological parameter calculation methods on the calculation results of single-frequency precise point positioning with increased tropospheric gradient
Wu et al. Retrieval of vertical distribution of tropospheric refractivity through ground-based GPS observation
CN115980317A (zh) 基于修正相位的地基gnss-r数据土壤水分估算方法
Wielgosz et al. Research on GNSS positioning and applications in Poland in 2015–2018
CN112528213B (zh) 基于低地球轨道卫星全球电离层总电子含量多层解析方法
Tunalı Water vapor monitoring with IGS RTS and GPT3/VMF3 functions over Turkey
Klügel et al. Atmospheric data set from the Geodetic Observatory Wettzell during the CONT-17 VLBI campaign

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 400020 Jiangbei District, Chongqing electric measuring Village No. 231

Applicant after: Chongqing geographic information and Remote Sensing Application Center (Chongqing surveying and mapping product quality inspection and testing center)

Applicant after: WuHan University

Address before: 430000 Hubei city of Wuhan province Wuchang Luojiashan

Applicant before: WuHan University

Applicant before: Chongqing Geographical Information Center

CB02 Change of applicant information
CB02 Change of applicant information

Address after: 430000 Hubei city of Wuhan province Wuchang Luojiashan

Applicant after: WuHan University

Applicant after: Chongqing geographic information and Remote Sensing Application Center (Chongqing surveying and mapping product quality inspection and testing center)

Address before: 400020 Jiangbei District, Chongqing electric measuring Village No. 231

Applicant before: Chongqing geographic information and Remote Sensing Application Center (Chongqing surveying and mapping product quality inspection and testing center)

Applicant before: WuHan University

GR01 Patent grant
GR01 Patent grant