陳西宏 吳文溢 劉 贊

?

基于改進(jìn)射線描跡法的對流層斜延遲估計(jì)

陳西宏 吳文溢*劉 贊

(空軍工程大學(xué)防空反導(dǎo)學(xué)院 西安 710051)

該文針對傳統(tǒng)對流層延遲模型和射線描跡法在估計(jì)對流層延遲方面的局限性，如效率低、成本高、精度受地表參數(shù)和探空數(shù)據(jù)限制等不足，提出一種基于改進(jìn)射線描跡法的對流層斜延遲估計(jì)方法。該方法結(jié)合中緯度大氣模式氣象參數(shù)公式和UNB3m氣象參數(shù)模型，改進(jìn)了射線描跡法中折射率剖面的計(jì)算，克服了氣象數(shù)據(jù)對射線描跡法的限制。選取亞洲地區(qū)10個(gè)站點(diǎn)2012年的氣象數(shù)據(jù)，分別采用改進(jìn)射線描跡法和傳統(tǒng)對流層延遲模型估計(jì)各個(gè)站天頂方向至高度角區(qū)間15個(gè)方向的對流層斜延遲，并與基于探空數(shù)據(jù)獲取的對流層斜延遲真值進(jìn)行比較，計(jì)算結(jié)果表明該方法的估計(jì)精度優(yōu)于傳統(tǒng)對流層延遲模型，為非氣象數(shù)據(jù)情況下對流層斜延遲實(shí)時(shí)估計(jì)提供了新的思路。

無線電探空；對流層斜延遲；射線描跡法；UNB3m模型；折射率

1 引言

針對現(xiàn)有對流層斜延遲估計(jì)方法的不足，本文提出一種利用改進(jìn)射線描跡模型估計(jì)對流層斜延遲的方法。此方法結(jié)合中緯度大氣模型和UNB3m模型中氣象參數(shù)插值公式，改進(jìn)了射線描跡法中折射率剖面的計(jì)算，在確保精度的同時(shí)，克服了描跡法對探空技術(shù)的依賴，具有成本低、效率高、操作簡單和實(shí)時(shí)性好等特點(diǎn)。通過對2012年日本地區(qū)內(nèi)4個(gè)氣象站數(shù)據(jù)的計(jì)算和分析，驗(yàn)證了此方法的可行性。

2 對流層斜延遲

電磁波信號在對流層內(nèi)因大氣折射而引起的傳播延遲稱為對流層延遲。在衛(wèi)星導(dǎo)航定位中，通常將任意高度角方向的對流層延遲表達(dá)為天頂延遲和映射函數(shù)的乘積。

3 對流層斜延遲估計(jì)方法

3.1基于探空數(shù)據(jù)的射線描跡法

射線描跡法在計(jì)算電磁波大氣折射誤差時(shí)改正精度高且適用于各種高度角，而基于探空數(shù)據(jù)的射線描跡法獲取的對流層斜延遲可以視為真值[19]。其簡化示意圖如圖1所示。

圖1 射線描跡法示意圖

3.2 UNB3m氣象參數(shù)模型

UNB3m模型[20,21]中的氣象參數(shù)是從海平面算起的，包括溫度(K)，大氣壓(mbar)，水汽壓(mbar)，大氣相對濕度(%)，溫度變化率(K/m)，水汽壓變化率(mbar/m)等，這些參數(shù)的解算根據(jù)測站點(diǎn)的地理位置和年積日，按照UNB3m對流層氣象參數(shù)格網(wǎng)值(如表1)進(jìn)行內(nèi)插獲得。內(nèi)插公式為

表1 UNB3m氣象參數(shù)格網(wǎng)值

3.3 改進(jìn)射線描跡法

通過對射線描跡法模型的分析，不難發(fā)現(xiàn)，如果能不依賴氣象觀測而建立測站點(diǎn)比較準(zhǔn)確的大氣折射率剖面模型，則能克服探空數(shù)據(jù)對描跡法的限制。根據(jù)球面大氣分層定律，大氣折射率可由各層的氣象參數(shù)計(jì)算得到，具體公式[11]為

在中緯度大氣模式中，不同海拔的氣象參數(shù)可以通過地表氣象參數(shù)計(jì)算獲取，其公式為

根據(jù)測站點(diǎn)的地理位置和年積日等信息，利用UNB3m氣象參數(shù)模型獲取氣象參數(shù)剖面后，再將其代入式(11)中可得對流層大氣折射率剖面：

結(jié)合射線描跡法中折射指數(shù)公式，令

通過利用UNB3m氣象參數(shù)模型和中緯度大氣模式氣象參數(shù)公式，建立測站的大氣折射率剖面模型，并基于此對射線描跡法進(jìn)行改進(jìn)，能有效地克服探空數(shù)據(jù)對描跡法的限制，快速精確地估計(jì)不同仰角的對流層斜延遲。

3.4對流層斜延遲估計(jì)方法

本文以基于探空數(shù)據(jù)的射線描跡法計(jì)算的對流層延遲作為參考值，分析改進(jìn)射線描跡法的估計(jì)精度，并與傳統(tǒng)對流層延遲模型的計(jì)算結(jié)果進(jìn)行比較分析。其研究步驟如下：

步驟2 根據(jù)亞洲地區(qū)10個(gè)站點(diǎn)2012年全年的探空和氣象數(shù)據(jù)(間隔時(shí)長為6 h)，利用式(11)獲得大氣折射率參考剖面；然后，分別在= 90, 80, 70, 60, 50, 30, 25, 20, 10, 9, 8, 7, 6, 5, 4(=0,1,,14)共計(jì)15個(gè)高度角方向進(jìn)行射線描跡，對流層層頂?shù)母叨瓤捎蒆opfield模型確定，其計(jì)算見式(16)。將,以及等參數(shù)代入式(4)~式(6)，可得到每站各高度角方向全年的對流層延遲參考值序列。

步驟3 根據(jù)3.3節(jié)中式(13)~式(15)，獲取折射率剖面；然后，重復(fù)步驟2的過程，將,以及等參數(shù)代入式(4)~式(6)，利用3.3節(jié)的改進(jìn)射線描跡法同樣估算得到每站上15個(gè)高度角全年的對流層延遲值序列；

步驟4 根據(jù)式(1)，利用Neil映射函數(shù)，將由對流層延遲模型Hopfield模型和Saastamoinen模型解算的對流層天頂延遲映射到上述的15個(gè)高度角方向獲取對流層斜延遲值序列，其中=1,2分別表示為模型和模型，計(jì)算公式為

步驟5 按式(18)計(jì)算得到3種方法在各高度角方向的對流層斜延遲偏差序列，下標(biāo)= 0,1,2分別代表模型，模型和改進(jìn)射線描跡法。

步驟6 對3種方法估算的各站點(diǎn)的延遲時(shí)間序列進(jìn)行統(tǒng)計(jì)分析，得出結(jié)論。

4 算例與結(jié)果分析

根據(jù)國際GPS服務(wù)站(International GPS Service, IGS)提供的亞洲地區(qū)10個(gè)IGS站點(diǎn)2012年的氣象和探空數(shù)據(jù)進(jìn)行對流層斜延遲分析。所選測站的信息如表2所示。

表2測站情況表

表2中，“Lat.”代表測站緯度(Latitude)，“Lon.”代表測站經(jīng)度(Longitude)，“Hgt.”代表測站海拔高度(Height)。

按照3.4節(jié)中的步驟進(jìn)行計(jì)算。圖2表示3種方法2012年全年對流層斜延遲估算結(jié)果的3維柱狀圖，為便于繪圖站點(diǎn)名用ID表示，并將偏差均值和均取絕對值。

圖2 3種方法2012年平均斜延遲偏差及RMS

圖3 各站高度角方向的偏差

5 結(jié)論

本文研究并設(shè)計(jì)了一種在缺少氣象數(shù)據(jù)情況下的對流層斜延遲估計(jì)方法，該方法利用中緯度大氣模式氣象參數(shù)公式和UNB3m氣象參數(shù)模型建立了對流層折射剖面模型，并利用該模型改進(jìn)了射線描跡法中折射指數(shù)的計(jì)算，使其克服了對探空數(shù)據(jù)的依賴。通過對IGS站實(shí)測氣象數(shù)據(jù)的計(jì)算和分析，結(jié)果表明此方法能在缺少氣象數(shù)據(jù)的情況下，較精確地估算出對流層斜延遲，且其精度優(yōu)于傳統(tǒng)對流層延遲模型，驗(yàn)證了此方法的可靠性和可行性。該方法對我國衛(wèi)星導(dǎo)航定位系統(tǒng)中對流層延遲實(shí)時(shí)監(jiān)測應(yīng)用具有重要的參考意義。

參考文獻(xiàn)

[1] 姚宜斌,曹娜, 許超鈐, 等. GPT2模型的精度檢驗(yàn)與分析[J].測繪學(xué)報(bào), 2015, 44(7): 726-733. doi: 10.11947/j.AGCS.2015. 20140356.

YAO Yibin, CAO Na, XU Chaoqian,. Accuracy assessment and analysis for GPT2[J]., 2015, 44(7): 726-733. doi: 10.11947/j. AGCS.2015.20140356.

[2] QU Weijing, ZHU Wenyao, SONG Shuli,. Evaluation of the precision of three tropospheric delay correction models[J]., 2008, 32(4): 429-438. doi: 10.1016/j.chinastron.2008.10.010.

[3] JOHANNES B, GREGOR M, MICHAEL S,. Development of an improved empirical model for slant delays in the troposphere (GPT2w)[J]., 2015, 19(3): 433-441. doi: 10.1007/s10291-014-0403-7.

[4] HEIKKI J, REIMA E, HENRIK V,A variational data assimilation system for ground-based GPS slant delays[J]., 2007, 133(625):969-980. doi: 10.1002/qj.79.

[5] 姚宜斌,張豹, 嚴(yán)鳳,等. 兩種精化的對流層延遲改正模型[J].地球物理學(xué)報(bào), 2015, 58(5): 1492-1501. doi: 10.6038/ cjg20150503.

YAO Yibin, ZHANG Bao, YANFeng,. Two new sophisticated models for tropospheric delay corrections[J]., 2015, 58(5): 1492-1501. doi: 10.6038/cjg20150503.

[6] HOPFIELD H S. Two-quartic tropospheric refractivity profile for correction satellite data[J]., 1969, 74(18): 4487-4499.

[7] SAASTAMOINEN J. Contribution to the theory of atmospheric refraction[J]., 1973, 107(1): 13-34.

[8] MARINI J W. Correction of satellite tracking data for an arbitrary tropospheric profile[J]., 1972, 7(2):223-231.

[9] NEIL A E. Global mapping functions for the atmospheric delay at radio wavelengths[J]., 1996, 101:3227-3246.

[10] 郭際明, 章迪, 史俊波, 等. 利用射線追蹤法分析三種典型對流層映射函數(shù)在中國區(qū)域的精度[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2015, 40(2): 182-187. doi: 10.13203/j.whugis20130096.

GUO Jiming, ZHANG Di, SHI Junbo,. Using ray-tracing to analyse the precision of three classical tropospheric mapping functions in china[J]., 2015, 40(2): 182-187. doi: 10.13203/j.whugis20130096.

[11] DAVID E, MACMILLAN D S, and JOHN M G. Tropospheric delay ray tracing applied in VLBI analysis[J].:, 2014, 119(12): 9156-9170. doi: 10.1002/2014JB011552.

[12] LANDON U, FELIPE G N, and MARCELO C S. Ray-traced slant factors for mitigating the tropospheric delay at the observation level[J]., 2012, 86(2): 149-160. doi: 10.1007/s00190-011-0503-x.

[13] HONG Zhenjie J, YE Fan, and GUO Peng. The ray-tracing mapping operator in an asymmetric atmosphere[J]., 2007, 18(1): 123-127. doi: 10.1016/j.pnsc. 2007.06.004.

[14] REZA GF, PETER D, and RICHARD B L. A web-based package for ray tracing the netural atmosphere radiomeric path delay[J]., 2008, 35(6): 1113-1124. doi: 10.1016/j.cageo.2008.02.027.

[15] 朱慶林,吳振森, 趙振維, 等. 單臺地基衛(wèi)星導(dǎo)航接收機(jī)測量對流層斜延遲[J]. 電波科學(xué)學(xué)報(bào), 2010,25(1):37-41. doi:10.13443/j.cjors.2010.01.008.

ZHU Qinglin, WU Zhensen, ZHAO Zhenwei,Tropospheric slant delay measured by singular ground-based satellite navigation receiver[J]., 2010,25(1):37-41. doi:10.13443/j.cjors.2010.01.008.

[16] 趙靜旸,宋淑麗, 陳欽明, 等.基于垂直剖面函數(shù)式的全球?qū)α鲗犹祉斞舆t模型的建立[J].地球物理學(xué)報(bào),2014,57(10):3140-3153. doi:10.6038/cjg20141005.

ZHAO Jingyang, SONG Shuli, CHEN Qinming,. Establishment of a new global model for zenith troposphere delay based on functional expression for its vertical profile[J]., 2014,57(10):3140-3153. doi:10.6038/cjg20141005.

[17] 黃良珂, 劉立龍, 文鴻雁, 等. 亞洲地區(qū)EGNOS天頂對流層延遲模型單站修正與精度分析[J].測繪學(xué)報(bào), 2014,43(8):807-817.doi:10.13485/j.cnki.11-2089.2014.0126.

HUANG Liangke, LIU Lilong, WEN Hongyan,Single-site improvement and accuracy analysis for zenith tropospheric delay of EGNOS model over asia area[J]., 2014,43(4): 807-817. doi:10.13485/j.cnki.11-2089.2014.0126.

[18] 毛健, 朱長青, 郭繼發(fā). 一種新的全球?qū)α鲗犹祉斞舆t模型[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2013, 38 (6): 684-688.

MAO Jian, ZHU Changqing, and GUO Jifa. A new global zenith tropospheric delay model[J]., 2013,38(6):684-688.

[19] 孫方, 王紅光, 康士峰, 等.大氣波導(dǎo)環(huán)境下的射線跟蹤算法[J].電波科學(xué)學(xué)報(bào),2008,23(1):179-183. doi:10.13443/j.cjors.2008.01.019.

SUN Fang, WANG Hongguang, KANG Shifeng,A ray tracing algorithm for duct environment [J]., 2008,23(1): 179-183. doi:10.13443/j.cjors.2008.01.019.

[20] LEANDRO R, SANTO M, and LANGLY R B. UNB neutral atmosphere model: Development and performance[C]. Proceeding of ION NTM 2006C. Monterey, USA, 2006: 564-573.

[21] LEANDRO R, LANGLY R B, and SANTO M. UNB3m_ pack: A neutral atmosphere delay package for radiometric space techniques[J]., 2008, 12(1): 65-70. doi: 10.1007/s10291-007-0077-5.

[22] SAKIDIN H. Nonhydrostatics UNBab mapping function simplification for tropospheric delay[J]., 2013, 53(12): 710-714. doi: 10.1016/j.proeng.2013.02.093.

Estimation of Tropospheric Slant Delay Based on Improved Ray Tracing Method

CHEN Xihong WU Wenyi LIU Zan

(,,’710051,)

The tradition tropospheric delay models and the method of ray-tracing have some limits like inefficiency, high cost and the the restrictions of radiosonde data and surface parameters. An improved method based on ray tracing is proposed. In this method, meteorological parameter formulas in the middle latitude model and meteorological parameter models in UNB3m model are combined to modify the arithmetic of refractive index, which gets rid of the restrictions of radiosonde data. Meteorological data from 10 Asia stations in 2012 are analyzed based on ray tracing technology and tradition models like Hopfield and Saastamoinen. Slant delays for 15 directions, from the zenith toelevations, are computed. The results are compared with ray-tracing tropospheric slant delays from nearby radiosonde measurements, which demonstrates that the accuracy of improved ray tracing method is superior to tradition models, and the proposed method provides a new real time way for estimating tropospheric slant delay in the case of lack of meteorological data.

Radiosonde; Tropospheric slant delay; Ray tracing; UNB3m model; Refractive index

P228

A

1009-5896(2016)10-2468-07

10.11999/JEIT160023

2016-01-11；改回日期：2016-06-06；網(wǎng)絡(luò)出版：2016-07-19

吳文溢 1440524558@qq.com

國家自然科學(xué)基金(61172169)

The National Natural Science Foundation of China (61172169)

陳西宏： 男，1961年生，教授，博士生導(dǎo)師，主要從事信息與控制技術(shù)的研究.

吳文溢： 男，1993年生，碩士生，研究方向?yàn)楦呔葧r(shí)間同步技術(shù).