gps信號對流層延遲映射函數(shù)的研究

gps信號對流層延遲映射函數(shù)的研究

gps信號從衛(wèi)星傳輸?shù)浇邮諜C的過程應(yīng)該受到層的影響。根據(jù)大氣層對信號傳播的不同影響,可將其分為電離層和對流層。電離層屬于彌散性介質(zhì),故其引起的信號傳播誤差可以通過不同頻率的相位組合消除;而對流層屬于非彌散性介質(zhì),故無法通過相位組合消除其對信號傳播的延遲。目前改正對流層延遲的方法是建立大氣模型,但是因?qū)α鲗拥拇髿獬煞趾軓?fù)雜,很難建立一個十分精確的模型,所以對流層延遲修正不夠精確,這已經(jīng)成為阻礙GPS定位精度提高的一個重要因素。1觀測衛(wèi)星高度角及以最大濕分量的集合式氣象機色模型vmf和集氣數(shù)集,假設(shè)有以下幾種形式,即對流層延遲映射函數(shù)在數(shù)學(xué)理論上可以分為兩類:第一類是以Saastamoine改正公式為代表的,首先把大氣折射積分中的被積函數(shù)按高度角(或天頂距)三角函數(shù)進行級數(shù)展開,然后逐項積分;第二類是以Marini連分式映射函數(shù)和在此基礎(chǔ)上發(fā)展起來的Chao、CfA2.2、Ifadis、MTT和NMF模型為代表,后來又出現(xiàn)了VMF1模型和GMF模型。1)NMF模型。是Niell在1996年利用一個全球分布的無線電探空氣球網(wǎng)的結(jié)果,是只與測站的地理緯度、高度和觀測日期有關(guān)的連分式映射函數(shù)系數(shù),并且把這些系數(shù)用一個與緯度相關(guān)的表列了出來。其加入了一個干延遲映射函數(shù)改正量Δmf(EL),并且其使用的列連分式系數(shù)表使它具有很廣的使用范圍。干分量的計算。ΜF(xiàn)d(EL)=mfd(EL)+Δmfd(EL)(1)mfd(EL)=1+ad/(1+bd/(1+cd))/(sin(EL)+ad/(sin(EL)+bd/(sin(EL)+cd)))(2)MFd(EL)=mfd(EL)+Δmfd(EL)(1)mfd(EL)=1+ad/(1+bd/(1+cd))/(sin(EL)+ad/(sin(EL)+bd/(sin(EL)+cd)))(2)Δmfd(EL)=[1/sinEL-f(EL?aht?bht?cht)]Η(3)ad(φ?t)=aavg(φ)-aamp(φ)cos(2π(t-28)/365.25)(4)式中,EL代表觀測衛(wèi)星的高度角,H是正高,單位為km,f(EL,aht,bht,cht)為式(2)表達的函數(shù),φ是測站緯度,t是觀測年積日,aavg,aamp在文獻中按測站緯度用列表的形式給出,bd,cd計算公式與ad一樣。aht,bht,cht分別為2.53×10-5,5.49×10-3,1.14×10-3。濕分量計算。ΜF(xiàn)w(EL)=1+aw/(1+bw/(1+cw))/(sinEL+aw/(sinEL+bw/(sinEL+cw)))(5)式中,aw,bw,cw在文獻中給出。2)VMF1模型。是全球范圍內(nèi)精度最高的對流層映射函數(shù)模型,可以在全球范圍內(nèi)依照2.5°×2.0°劃分的各區(qū)域內(nèi)每6h提供一套對流層映射函數(shù)。VMF1模型將連分式映射函數(shù)中的系數(shù)b,c固定,通過NWM(NumercialWeatherModels)確定初始高度角為3.3°時的映射函數(shù)值,再代入式(1)或式(5)反求系數(shù)a。VMF1模型重新定義了干延遲映射函數(shù)的系數(shù)b,c。在VMF1模型中bd=0.0029。cd=c0+[c11/2(cos(2π(doy-28)/365+ph)+1)+c10](1-cosφ)(6)當(dāng)測站位于北半球時,計算cd用到的系數(shù)c0,c10,c11和ph分別為0.062、0.001、0.005和0;當(dāng)測站位于南半球時,分別為0.062、0.002、0.007和π。VMF1模型中用到的參數(shù)每61)和式(5)計算。3)GMF模型。是由J.Boehm(2006)等人提出的一種新的全球映射函數(shù),其投影函數(shù)是基于全球ECMWF數(shù)字氣象模型數(shù)據(jù)建立的。干延遲包括與地理測站高程有關(guān)的改正,反映了大氣密度隨高度增加而減少的變化率。在GMF模型中輸入的參數(shù)是年積日、測站經(jīng)緯度和測站高程。這也是GAMIT軟件中默認(rèn)使用GMF模型的原因。計算干分量。計算方法與NMF模型干分量計算方法一樣,參照式(1)和式(2)的算法,但式中bd=0.0029。ad=55∑i=1{[ahm(i)ap(i)+bhm(i)bp(i)]×10-5+[aha(i)ap(i)+bha(i)bp(i)]×10-5}式中,ap(i)=9∑n=0n∑m=0Ρn+1?m+1cos(mλ);Ρn+1?m+1=9∑i=0min(j?m)∑j=0{(1/2)i+√(1-sin2φ)jS}?S=int(i+j/2)∑k=0(-1)k(2i-2k+1)！(i-j-2k+1)！(sinφ)i-j-2k/(k+1)！(i-k+1)！。當(dāng)測站處于北半球時,φ≥0,ph=0,ch1=0.005,cho=0.001;當(dāng)測站處于南半球時,φ≤0,ph=π,ch1=0.007,ch0=0.002。干延遲映射函數(shù)的改正量計算方法與NMF模型相同,同式(3)。計算濕分量。濕延遲映射函數(shù)同式(5),但系數(shù)計算方法不同。GMF模型中用到的系數(shù)為:aw=55∑i=1{[awm(i)ap(i)+bwm(i)bp(i)]×10-5+[awa(i)ap(i)+bwa(i)bp(i)]×10-5}×cos((doy-28)2π/365.25)bw=0.00146?cw=0.04391干、濕映射函數(shù)中用到的ahm,bhm,aha,bha,awm,bwm,awa,bwa均可按模型中提供的表查詢。2種模型的對比利用德國境內(nèi)的SAPOS網(wǎng)(Satellitenpositionierungsdienst)的10個測站,連續(xù)7天(2009年042-048天),每天24h,采樣間隔為30s的觀測數(shù)據(jù),采用IGS精密星歷,使用GAMIT10.34版軟件解算基線。天頂延遲模型采用Saastamoinen模型,氣象數(shù)據(jù)全部采用標(biāo)準(zhǔn)氣象數(shù)據(jù)(1013.25mb,20℃,50%RH),對流層天頂延遲參數(shù)約束設(shè)置為0.5m,這樣的設(shè)置保證天頂延遲參數(shù)包含所有的濕延遲。天頂延遲變化設(shè)置為0.02m/√h,相關(guān)時間值為100h,對流層參數(shù)個數(shù)設(shè)為13個。在截止高度角為7°和15°時分別用三種模型計算觀測數(shù)據(jù),比較數(shù)據(jù)處理的結(jié)果。評價基線解算結(jié)果通常有兩種標(biāo)準(zhǔn)。1)在GAMIT解算結(jié)果中,NRMS表示單時段解算出的基線值偏離其加權(quán)平均值的程度,是從歷元模糊度解算中得出的殘差,是衡量GAMIT解算結(jié)果的一個重要指標(biāo),其計算公式為:ΝRΜS=(1/Νn∑i=1(Yi-Y)2/σ2i)1/2一般來說,NRMS值越小,基線估算精度越高。根據(jù)國內(nèi)外GPS數(shù)據(jù)處理經(jīng)驗,其值一般應(yīng)小于0.3。截止高度角為7°時的NRMS值如表1所示。截止高度角為15°時的NRMS值如表2所示。從表1、表2可以看出,不同模型計算的基線NRMS都小于0.3,滿足數(shù)據(jù)處理的要求,三個模型相互之間的NRMS差值在0～0.0002之間,說明各個模型解算基線的精度相近,沒有較大差異。2)基線分量的重復(fù)性反映了單天解之間的內(nèi)符合精度,是衡量GPS基線解算結(jié)果的重要質(zhì)量指標(biāo)之一。計算基線向量的重復(fù)性的公式為:Rl=[n/(n-1)n∑i=1(Li-ˉL)2/δ2i/n∑i=11/δ2i]12式中,Rl為基線向量的重復(fù)性,n為基線單天解的數(shù)目,Li為第i天的基線分量(或邊長),ˉL為單天解基線分量的加權(quán)平均值。在衛(wèi)星截止高度角為7°時,基線重復(fù)性統(tǒng)計如表3所示。從表3可以看出,當(dāng)衛(wèi)星截止高度角設(shè)為7°時,模型之間的基線重復(fù)性在各個方向上都沒有明顯的差異,差值最大為0.0001。衛(wèi)星截止高度角為15°時,基線0577-0597的重復(fù)性如圖1所示?；€0577-0601的重復(fù)性如圖2所示。從圖1和圖2可以看出:衛(wèi)星截止高度角設(shè)為15°時,三個模型之間的精度一致,差異最大為2mm。對同一組數(shù)據(jù),衛(wèi)星截止高度角為15°時基線的重復(fù)性,優(yōu)于衛(wèi)星截止高度角為7°的情況。大氣延遲主要影響垂直方向的定位精度,從圖1圖2也可以清楚地看出,垂直方向的基線重復(fù)性比其他三個方向大。通過以上的數(shù)據(jù)處理實例可以看出,在中高緯度地區(qū)各種映射函數(shù)模型的連分式之間,僅在形