CN106802425B - 一种估算天顶对流层延迟的积分方法 - Google Patents
一种估算天顶对流层延迟的积分方法 Download PDFInfo
- Publication number
- CN106802425B CN106802425B CN201710047564.4A CN201710047564A CN106802425B CN 106802425 B CN106802425 B CN 106802425B CN 201710047564 A CN201710047564 A CN 201710047564A CN 106802425 B CN106802425 B CN 106802425B
- Authority
- CN
- China
- Prior art keywords
- atmospheric
- delay
- refractive index
- hydrostatic
- zenith
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000010354 integration Effects 0.000 title claims abstract description 9
- 230000002706 hydrostatic effect Effects 0.000 claims abstract description 45
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 25
- 239000005436 troposphere Substances 0.000 claims description 23
- 239000000203 mixture Substances 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims description 2
- 238000012986 modification Methods 0.000 abstract description 3
- 230000004048 modification Effects 0.000 abstract description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提出一种估算对流层天顶延迟的积分方法,在计算天顶对流层静水延迟和天顶对流层湿度延迟时分别考虑了相邻两高度层上静水折射指数和大气压强的改变量以及湿度折射指数与大气压强的改变量,并分别将静水折射指数和湿度折射指数转为大气压强随高度变换的函数,然后对其进行积分后分别估算出天顶对流层静水延迟和天顶对流层湿度延迟。在空间分布上,当存在水汽压“逆增层”和大气温度“逆温”的情况下,新的积分模型能够更加精确的估算ZTD。
Description
技术领域
本发明属于全球卫星导航系统气象学技术领域,特别涉及一种估算天顶对流层延迟的积分方法。
背景技术
全球卫星导航系统(后简称为GNSS)卫星信号经过对流层时,对流层对它产生传播速度延迟和传播路径弯曲延迟两部分影响。由此产生的路径延迟被称作对流层延迟。在GNSS数据处理的过程中,通常将信号传播路径的对流层延迟投影到天顶方向,并将其分为天顶对流层静水延迟和天顶对流层湿度延迟两部分。再已知大气静水折射指数Nd h和大气湿度折射指数Nw h后,通常采用线性离散化的方法估算出天顶静水延迟ZHD及天顶湿度延迟ZWD。事实上,大气静水折射指数Nd h和大气湿度折射指数Nw h在垂直方向上呈现指数形式的变换,故线性离散化估算ZHD和ZWD会进入较大的误差。
发明内容
针对上述问题,本发明提出一种新的估算对流层天顶延迟的积分方法,用于解决目前卫星定位系统气象学中,利用探空产品、无线掩星产品、ERA-Interim产品等采用线性离散化的积分模型估算天顶对流层延迟时存在误差较大的问题。
为解决上述技术问题,本发明采用如下技术方案:
一种估算对流层天顶延迟的积分方法,其特征在于包括如下步骤:
首先根据垂直与地平面的高度方向上,大气压强、大气水汽压和大气温度的空间分布信息,分别估算出不同高度层上的大气静水折射指数Nd h和大气湿度折射指数Nw h;
然后分别得出大气静水折射指数与大气压强在空间分布中存在的关系、大气湿度折射指数与大气水汽压强在空间分布中存在的关系、以及大气压强与大气水汽压强在空间分布中存在的关系;
之后,在估算天顶对流层静水延迟时,考虑相邻两高度层上大气静水折射指数和大气压强的改变量,并将大气静水折射指数转化为大气压强随高度变换的函数;在估算天顶对流层湿度延迟时,在相邻两高度层上,考虑大气湿度折射指数、大气水汽压强、大气压强和大气水汽压强的改变量,并最终将大气湿度折射指数转换为大气压强随高度变换的函数;
最后,对两个函数分别进行积分后既可以估算出天顶对流层静水延迟和天顶对流层湿度延迟。
进一步的,在垂直与地平面的高度方向上,大气静水折射指数Nd h与大气压强的关系、大气湿度折射指数Nw h与水汽压的关系、以及水汽压强和大气压强的关系表示如下:
其中,Nd h表示大气静水折射指数;Nw h表示大气湿度折射指数;τd为静水延迟混合比指数;τw为湿度延迟混合比指数;e代表为水汽压强;P为大气压强;e0和P0表示地表处的水汽压强和大气压强;γ定义为大气混合比指数;
大气压强随高度的分布通常用指数形式表示:
Pi=P0exp(H·hi) (4)
其中,H为大气层等效高度,且h为高度坐标。
进一步的,选取地表向上不同高度层hi(i=0,1,2,3,…,n)的高度区间[h1,h2],对大气静水折射指数Nd h和大气湿度折射指数Nw h分别积分,估算天顶静水延迟ZHD及天顶湿度延迟ZWD:
由公式(4)可知,在区间[h1,h2]中:
P=P1exp(H1,2·(h-h1)) (6)
再根据公式(6),公式(5)可以变换为:
其中,
之后,在全部的高度层上,通过累积相加所有高度层估算的ZHD来求取
整个大气层的天顶对流层静水延迟ZHD为:
其中,
同理,
根据公式(9),天顶对流层湿度延迟可以表示为:
其中,
最后,得到天顶对流层延迟ZTD为:
本发明针对水汽压存在“逆增层”时、或者温度存在“逆温”的情况下,采用无线探空产品、无线掩星产品或ERA-Interim产品等通过线性的离散化估算天顶对流层延迟会引入较大的误差的问题,提出顾及大气折射指数与大气压强在空间分布上呈现近似指数变换的形式,并对原有估算方法进行重组,在计算天顶对流层静水延迟和天顶对流层湿度延迟时分别考虑了相邻两高度层上静水折射指数和大气压强的改变量以及湿度折射指数与大气压强的改变量,并分别将静水折射指数和湿度折射指数转为大气压强随高度变换的函数,然后对其进行积分后分别估算出天顶对流层静水延迟和天顶对流层湿度延迟。相对于现有的误差估算方式,本发明只需要引入相邻两高度层上静水折射指数和大气压强的改变量以及湿度折射指数与大气压强的改变量,就使得最后的结果更加客观科学,参数易于获取、整个估算方法更加简单便捷。在空间分布上,存在水汽压“逆增层”现象或者温度存在“逆温”的情况下,新的积分估算方法能够提高估算对流层延迟的精度。
附图说明
图1为本发明实施例的流程图。
具体实施方式
以下结合附图和实施例详细说明本发明技术方案。
首先根据无线探空产品、无线掩星产品和ERA-Interim产品提供的高垂直分辨率的大气压强、大气水汽压和大气温度的空间分布信息,分别估算出不同高度层上的静水折射指数和湿度折射指数。并分别得出静水折射指数与大气压强、湿度折射指数与大气水汽压以及大气压强与大气水汽压在空间分布中存在的关系。在估算天顶对流层静水延迟时,顾及了相邻两高度层上静水折射指数和大气压强的改变量,并将静水折射指数转化为大气压强随高度变换的函数。在估算天顶对流层湿度延迟时,在相邻两高度层上,顾及了湿度折射指数与大气水汽压、大气压强和大气水汽压的改变量,并最终将湿度折射指数转换为大气压强随高度变换的函数。
在已知大气静水折射指数Nd h和大气湿度折射指数Nw h后,估算天顶静水延迟ZHD及天顶湿度延迟ZWD的函数模型如公式(1.1)和(2.1)所示:
由公式(1.1)和(2.1)可以看出ZHD和ZWD都是一个积分量,分别与不同高度上的静水折射率Nd h和湿度折射率Nh w相关。实际计算中,可采用高垂直分辨率的探空资料、无线掩星观测资料、ERA-Interim产品等。这些产品提供了从地表向上不同高度层hi(i=0,1,2,3,…,n)的温度Ti与水汽压ei。根据这些气象产品可以事先估算出每个高度层上的静水折射指数和湿度折射指数,然后再分别对公式(1.1)和公式(2.1)中的积分进行线性离散化可以得到:
通常情况下,Nd h和Nw h随高度呈现出近似指数的形式变化,故采用线性的模型对公式(1.1)进行积分会存在较大的误差。为了减弱线性离散化对估算ZHD和ZWD的影响,我们分别对公式(1.1)和(2.1)进行了重新处理。根据无线探空、无线掩星等高垂直分辨率、高精度的气象观测值进行研究发现:Nd h与P,Nw h与e之间的关系可近似表示为:
在垂直方向上,大气水汽压和大气压强之间的关系可以近似的表示为:
其中,e0和P0表示地球表面处的水汽压强和大气压强;γ定义为大气混合比指数,P表示大气压强,e表示水汽压强。
大气压强随高度的分布通常可以用指数形式表示:
Pi=P0exp(H·hi) (4)
P0为地球表面的大气压强;H为大气层等效高度,且h为高度坐标。
根据公式(1),在区间[h1,h2]中,Nd h的积分可以表示为:
再利用公式(3)给出的压强随高度变换的关系式,区间[h1,h2]中,任意高度处的压强可以表示为:
P=P1exp(H1,2·(h-h1)) (6)
再根据公式(6),公式(5)可以变换为:
其中,
公式(10)仅仅给出了其中一个高度层内ZHD的离散化方法。在全部的高度层上,可以通过累积相加所有高度层估算的ZHD来求取整个大气层的ZHD:
其中,
利用公式(8)的模型离散化公式(1.1)来估算ZHD,考虑了不同高度处大气静水折射指数Nd h与大气压强P之间的近似指数变换的关系。并将静水折射指数转化为大气压强随高度变换的函数,然后对其进行积分。
同理,联合公式(2)、(6)和(3),在区间[h1,h2]中,Nw h的积分可以表示为:
其中,
公式(13)仅给出了某一个高度层内估算ZWD的积分算法。在整个大气层内,将所有单个高度层内利用公式(13)估算出来,然后求和即可获得整个大气层的ZWD。
其中,
联合公式(8)和公式(10),估算天顶对流层延迟ZTD如下式如示:
公式(11)提供的积分方法估算ZTD时,顾及了单个高度层内,静水折射指数与大气压强之间的关系及湿度折射指数与大气压强之间的关系,并分别将静水折射指数和湿度折射指数转换为大气压强随高度变化的函数。在空间分布上,当存在水汽压强“逆增层”现象和大气温度存在“逆温”的情况下,新的积分方法可以提高估算ZTD的精度。
以上实施例仅供说明本发明之用,而非对本发明的限制,有关技术领域的技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变换或变型,因此所有等同的技术方案,都落入本发明的保护范围。
Claims (2)
1.一种估算对流层天顶延迟的积分方法,其特征在于包括如下步骤:
首先根据垂直于地平面的高度方向上的大气压强、大气水汽压强和大气温度的空间分布信息,分别估算出不同高度层上的大气静水折射指数Nd h和大气湿度折射指数Nw h;
然后分别得出大气静水折射指数与大气压强在空间分布中存在的关系、大气湿度折射指数与大气水汽压强在空间分布中存在的关系、以及大气压强与大气水汽压强在空间分布中存在的关系;
之后,在估算天顶对流层静水延迟时,考虑相邻两高度层上大气静水折射指数和大气压强的改变量,并将大气静水折射指数转化为大气压强随高度变换的函数;在估算天顶对流层湿度延迟时,在相邻两高度层上,考虑大气湿度折射指数、大气水汽压强、大气压强和大气水汽压强的改变量,并最终将大气湿度折射指数转换为大气压强随高度变换的函数;
最后,对两个函数分别进行积分后即可以估算出天顶对流层静水延迟和天顶对流层湿度延迟;
其中,在垂直于地平面的高度方向上,大气静水折射指数Nd h与大气压强P的函数关系、大气湿度折射指数Nw h与水汽压强e的函数关系、以及水汽压强e和大气压强P的函数关系表示如下:
其中,τd为静水延迟混合比指数;τw为湿度延迟混合比指数;e0和P0表示地表处的水汽压强和大气压强;γ定义为大气混合比指数;
大气压强随高度的分布用指数函数形式表示:
Pi=P0 exp(H·hi) (4)
其中,H为大气层等效高度,且h为高度坐标。
2.根据权利要求1所述的估算对流层天顶延迟的积分方法,其特征在于选取地表向上不同高度层hi的高度区间[h1,h2],其中,i=0,1,2,3,…,n;对大气静水折射指数Nd h和大气湿度折射指数Nw h分别积分,估算天顶对流层静水延迟ZHD及天顶对流层湿度延迟ZWD:
由公式(4)得到公式(6)P=P1 exp(H1,2·(h-h1))并代入公式(5)后:
其中,
之后,在全部的高度层上,通过累积相加所有高度层估算的ZHD来求取整个大气层的天顶对流层静水延迟ZHD为:
其中,
同理,
根据公式(9),天顶对流层湿度延迟可以表示为:
其中,
最后,得到天顶对流层延迟ZTD为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710047564.4A CN106802425B (zh) | 2017-01-22 | 2017-01-22 | 一种估算天顶对流层延迟的积分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710047564.4A CN106802425B (zh) | 2017-01-22 | 2017-01-22 | 一种估算天顶对流层延迟的积分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106802425A CN106802425A (zh) | 2017-06-06 |
CN106802425B true CN106802425B (zh) | 2019-07-23 |
Family
ID=58987069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710047564.4A Active CN106802425B (zh) | 2017-01-22 | 2017-01-22 | 一种估算天顶对流层延迟的积分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106802425B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108008416B (zh) * | 2017-12-04 | 2020-01-14 | 武汉大学 | 一种估算斜路径对流层延迟的积分方法 |
CN108920414B (zh) * | 2018-05-18 | 2022-03-29 | 中国人民解放军61540部队 | 一种利用气象资料计算局部天顶对流层湿延迟的新方法 |
CN111679346B (zh) * | 2019-12-27 | 2022-06-21 | 广东电网有限责任公司电力科学研究院 | 一种大气层可降水量预估方法及装置 |
CN111538943B (zh) * | 2020-04-24 | 2021-05-18 | 桂林理工大学 | 新的高时空分辨率全球ztd垂直剖面格网模型构建方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5731786A (en) * | 1994-12-29 | 1998-03-24 | Trimble Navigation Limited | Compaction of SATPS information for subsequent signal processing |
CN103558609A (zh) * | 2013-11-08 | 2014-02-05 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统对流层天顶延迟修正的方法 |
CN104965207A (zh) * | 2015-05-19 | 2015-10-07 | 同济大学 | 一种区域对流层天顶延迟的获取方法 |
CN105116468A (zh) * | 2015-06-24 | 2015-12-02 | 中国人民解放军63655部队 | 一种边界层大气湍流强度高度分布的检测方法 |
EP2995972A1 (en) * | 2014-09-15 | 2016-03-16 | Fugro N.V. | Integer ambiguity-fixed precise point positioning method and system |
CN105676239A (zh) * | 2014-11-19 | 2016-06-15 | 香港理工大学 | 一种gnss大气折射延迟估值方法和系统 |
CN105785407A (zh) * | 2016-02-23 | 2016-07-20 | 东南大学 | 一种适用于中国地区的无气象参数对流层延迟改正方法 |
CN106324620A (zh) * | 2016-08-02 | 2017-01-11 | 中国人民解放军空军工程大学 | 一种不依赖地表气象数据实时测量的对流层天顶延迟方法 |
-
2017
- 2017-01-22 CN CN201710047564.4A patent/CN106802425B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5731786A (en) * | 1994-12-29 | 1998-03-24 | Trimble Navigation Limited | Compaction of SATPS information for subsequent signal processing |
CN103558609A (zh) * | 2013-11-08 | 2014-02-05 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统对流层天顶延迟修正的方法 |
EP2995972A1 (en) * | 2014-09-15 | 2016-03-16 | Fugro N.V. | Integer ambiguity-fixed precise point positioning method and system |
CN105676239A (zh) * | 2014-11-19 | 2016-06-15 | 香港理工大学 | 一种gnss大气折射延迟估值方法和系统 |
CN104965207A (zh) * | 2015-05-19 | 2015-10-07 | 同济大学 | 一种区域对流层天顶延迟的获取方法 |
CN105116468A (zh) * | 2015-06-24 | 2015-12-02 | 中国人民解放军63655部队 | 一种边界层大气湍流强度高度分布的检测方法 |
CN105785407A (zh) * | 2016-02-23 | 2016-07-20 | 东南大学 | 一种适用于中国地区的无气象参数对流层延迟改正方法 |
CN106324620A (zh) * | 2016-08-02 | 2017-01-11 | 中国人民解放军空军工程大学 | 一种不依赖地表气象数据实时测量的对流层天顶延迟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106802425A (zh) | 2017-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106802425B (zh) | 一种估算天顶对流层延迟的积分方法 | |
Rocken et al. | Improved mapping of tropospheric delays | |
CN104965207B (zh) | 一种区域对流层天顶延迟的获取方法 | |
Wilgan et al. | Multi-observation meteorological and GNSS data comparison with Numerical Weather Prediction model | |
CN110031877B (zh) | 一种基于grnn模型的区域nwp对流层延迟改正方法 | |
CN106814373B (zh) | 大气加权平均温度估算及对流层延迟积分方法 | |
Stephan et al. | Ship-and island-based atmospheric soundings from the 2020 EUREC 4 A field campaign | |
CN109543353A (zh) | 三维水汽反演方法、装置、设备和计算机可读存储介质 | |
Benevides et al. | 4D wet refractivity estimation in the atmosphere using GNSS tomography initialized by radiosonde and AIRS measurements: results from a 1-week intensive campaign | |
CN107037439A (zh) | 针对陆地目标的激光测高仪大气延迟测距误差修正方法 | |
CN108008416B (zh) | 一种估算斜路径对流层延迟的积分方法 | |
CN113009531A (zh) | 一种小尺度高精度的低空对流层大气折射率模型 | |
CN101702037B (zh) | 基于遥控飞艇的近地层气象要素廓线剖面数据获取方法 | |
Liangke et al. | A zenith tropospheric delay correction model based on the regional CORS network | |
TWI424183B (zh) | 定位方法與裝置 | |
CN115061167A (zh) | 一种适用于短距离大高差rtk的对流层延迟改正方法 | |
CN113093241B (zh) | 一种考虑高度角的单测站对流层斜延迟计算方法 | |
CN112711022B (zh) | 一种GNSS层析技术辅助的InSAR大气延迟改正方法 | |
CN113960642A (zh) | 青藏高原区域光学遥感卫星影像大气折射定位误差修正方法 | |
Wilgan et al. | Comparison of tropospheric parameters from Meteodrone measurements with GNSS estimates from ground-based stations | |
CN111273319B (zh) | 一种基于余弦函数的区域对流层湿延迟计算方法 | |
De Haan et al. | Comparison of GPS slant delay measurements to a numerical model: case study of a cold front passage | |
CN111199092A (zh) | 太阳辐射遥感估算方法、系统及数据处理装置 | |
Kim et al. | Estimation of monthly precipitation in North Korea using PRISM and digital elevation model | |
Rózsa et al. | The role of GNSS infrastructure in the monitoring of atmospheric water vapor |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |