連月勇,張 超,謝宗特
(1.信息工程大學 導航與空天目標工程學院,鄭州 450001;2.75711部隊,廣州 510000)
地球定向參數(shù)對天文定位定向的影響分析
連月勇1,張 超1,謝宗特2
(1.信息工程大學 導航與空天目標工程學院,鄭州 450001;2.75711部隊,廣州 510000)
在天文定位定向中需要用到地球定向參數(shù)的預報值,實現(xiàn)地心地固坐標系和地心慣性坐標系之間的相互轉換,但在實時解算時,無法獲得最新的地球定向參數(shù)值,研究地球定向參數(shù)預報誤差將有助于提高天文解算精度。介紹了天文定位定向原理,對地球定向參數(shù)在其中的應用進行了研究,并分別分析了地球定向參數(shù)短期和長期預報值的精度,以及參數(shù)預報誤差對天文定位定向的影響。結果表明,一等天文測量中緯度和方位角的解算可采用365d地球定向參數(shù)預報值,經(jīng)度的解算需要進行每季度預報數(shù)據(jù)的更新。以上結論在天文導航和天文測量中有較好的借鑒意義。
地球定向參數(shù);天文定位;天文定向;極移;預報誤差
國際地球自轉服務(International Earth Rotation Service,IERS)于1988年由國際大地測量學和地球物理學聯(lián)合會與國際天文學聯(lián)合會共同建立,綜合全球各個新技術觀測處理中心的結果開展地球定向參數(shù)服務。地球定向參數(shù)(Earth orientation parameters,EOP)是一組反應地球自轉變化的參數(shù),包括歲差、章動、極移、UT1-UTC和地球自轉角(Earth rotation angle,ERA)等。EOP參數(shù)描述了地心地球參考系(geocentric terrestrial reference system,GTRS)和地心天球坐標系(geocentric celestial reference system,GCRS)之間的轉換關系,應用于天文測量、天文導航和衛(wèi)星導航等各個領域中[1]。
目前所獲取的EOP參數(shù),由于解算過程復雜,并不能提供參數(shù)的實時數(shù)據(jù),因此進行天文定位定向時經(jīng)常使用的是EOP參數(shù)的預報值。公報A文件是快速服務和預報產(chǎn)品,每周更新一次,數(shù)據(jù)文件包含各項地球定向參數(shù)及其計算誤差,文件格式為ASCII,包含了跳秒和地球自轉參數(shù)等信息,其中EOP參數(shù)的最長預報時間為1 a。
1991年至2000年,由于天體測量精度的不斷提高及甚長基線干涉技術(very long baseline interferometry,VLBI)的發(fā)展等極大的提高了天體位置的觀測精度,國際天文學協(xié)會(International Astronomical Union,IAU)通過了一系列關于天文參考系理論、地球自轉模型和時間尺度的決議[2-4]。其中,IAU2000A歲差章動模型的預報精度可達0.2 mas,IAU2000B歲差章動模型預報精度達1 mas,在進行天文定位定向解算時可忽略歲差章動參數(shù)的誤差不計,因此本文將重點針對EOP參數(shù)中的極移分量xp、yp和時間尺度關系值UT1-UTC進行具體的分析。
2.1 天文定位原理
如圖1所示,QWQ′是赤道,NWS為地平圈,P為北天極,Z為天頂,σ1、σ2及σ3為天球上的等高圈上的三顆恒星。由球面公式得[5-6]
cosz=sinφsinδ+cosφcosδcost
(1)
式(1)中,z為天體的天頂距,φ為測站天文緯度,δ天體的赤緯,t是測瞬天體的時角。
t=S-α+λ
(2)
式(2)中,α為天體的赤經(jīng),λ為測站的天文經(jīng)度,S為觀測瞬間的格林尼治真恒星時,由觀測瞬間的UTC時刻T觀解算得到。采用多星觀測數(shù)據(jù)聯(lián)立式(1)即可解得天文經(jīng)緯度值。
圖1 天球上球面三角形示意圖
2.2 天文定向原理
天文定向方法包括時角法和高度法兩種,下面僅介紹時角法天文定向的原理。如圖2所示,AB是地面目標B的天文方位角,Aσ是觀測時刻自然天體σ的方位角。設在格林尼治世界時Tutc觀測得到自然天體的水平度盤值為Lσ, 地面目標的水平度盤值為LB, 則地面目標的方位角AB可通過式(3)解得[7-8]
AB=LB-(Lσ-Aσ)
(3)
式(3)中,天體方位角計算公式
(4)
式(4)中
t=S0+Tut1+Tut1×u+λ-α
(5)
式(5)中,α為自然天體在測瞬時刻的視赤經(jīng),δ為自然天體在測瞬時刻的視赤緯,t為天體的時角,λ和φ分別為測站的天文經(jīng)度和緯度。u=0.002 737 91,S0為世界時0時對應的格林尼治恒星時。
圖2 天文定向原理圖
2.3 恒星視位置計算
由2.1節(jié)及2.2節(jié)可知,觀測瞬間恒星的真赤道坐標α、δ是天文定位定向解算必不可少的參數(shù),恒星視位置計算的精確度將直接影響定位定向的解算精度。恒星的星表歷元平位置換算到觀測歷元的視位置,即視位置的計算,其計算流程如圖3所示。
圖3 恒星視位置計算流程圖
此時的觀測位置所解算的天文經(jīng)緯度(λ、φ)和方位角(a)是相對于測瞬真北極的,由于地極移動使地面各點的子午圈發(fā)生變動,需要利用EOP的極移參數(shù)將天文成果歸算到平北極(λ0、φ0、a0)。改正方法為[9]
(1)經(jīng)緯度的極移改正為
(6)
(2)方位角的極移改正為
(7)
利用2010年到2013年6月這三年半期間所公布的182個公報A文件,對IERS發(fā)布的EOP參數(shù)的預報誤差從短期預報和長期預報[10]兩個方面進行詳細分析。
3.1 短期預報誤差分析
1)1 d預報誤差
把IERS網(wǎng)站所發(fā)布的2010年到2014年的CO4文件所提供的EOP參數(shù)作為真值,將以上的182個公報A文件中的第一天預報值同對應真值進行比較,得到了182個1 d預報數(shù)值的誤差,其分布如圖4所示。
從圖4可以看出,1 d預報誤差精度較高,且穩(wěn)定性好,浮動較小,并以零值為中心呈現(xiàn)不規(guī)則分布。極移分量xp的誤差絕對值基本上在1mas以內(nèi),極移分量yp的誤差絕對值大都在0.5mas以內(nèi),UT1-UTC誤差絕對值都在0.2mas內(nèi)。
圖4 1 d的EOP參數(shù)預報誤差
表1給出了1d預報誤差值的最大絕對值、均值和標準差,作為預報精度的參考。
表1 1 d的EOP預報誤差指標/mas
2)7 d預報誤差
因為IERS網(wǎng)站每周進行一次公報A的更新,研究前7 d的預報值精度對天文測量和天文導航具有很好的借鑒意義[11]。將以上的182個公報A文件中的第一周預報數(shù)據(jù)同對應真值進行比較,得到參數(shù)的誤差值如圖5所示。
從圖5可以看出,7 d預報誤差隨著天數(shù)增加而增大,每周的誤差情況相似,誤差值比較穩(wěn)定,整體上呈現(xiàn)以零值為中心不規(guī)則分布。極移分量xp的誤差絕對值在10mas以內(nèi),極移分量yp的誤差絕對值都在5mas以內(nèi),UT1-UTC誤差絕對值大都在1mas。
表2給出了3.5(a)內(nèi)7d預報誤差的最大絕對值、均值和標準差。
3.2 長期預報誤差分析
選取2010年到2012年發(fā)布的公報A文件和之后發(fā)布的相應的CO4文件,兩者比對求出1 a預報值和真值的差值,結果如圖6所示。
圖5 7 d的EOP參數(shù)預報誤差
表2 7 d的EOP預報誤差指標/mas
因為極移中包含周年周期和張德勒等分量,故而其預測有周期性,由圖6可看出,1 a的EOP誤差值存在一定的周期性。表3給出了三年總的及各個年份的EOP參數(shù)的各項指標。
由表3可知,xp全年的預報誤差在67mas以內(nèi),yp全年的預報誤差在80mas以內(nèi),UT1-UTC全年的預報誤差在62mas以內(nèi),2011年的各個參賽預報誤差最小,隨著年份推移,整體預報準確性有所提高,但不穩(wěn)定。
圖6 長期EOP參數(shù)預報誤差
表3 1 a的EOP預報誤差指標/mas
分量最大誤差的絕對值均值標準差xp66.5521.1150329.02378yp79.559-0.1416529.04795UT1-UTC62.0046-5.0125123.3792764.32320.8790926.17982010年79.559-11.980336.7674819.7479-0.7657311.8919341.862-8.3256512.839992011年28.263-3.5533413.5034162.0046-28.711822.3210866.552-9.2083433.069512012年50.84315.1087224.8268522.679514.439976.511917
根據(jù)文獻[12]規(guī)定的各等天文測量要求(如表4所示),分析以上三個EOP參數(shù)不同長度預報的誤差影響大小。
表4 天文等級測量標準
4.1 日常預報誤差
EOP參數(shù)中的UT1-UTC稱為日長預報,下面利用第二節(jié)中EOP預報誤差的分析結果研究日長預報誤差對天文定位定向的影響。
由天文定向原理可知,日長預報對于緯度的測定不產(chǎn)生影響。UT1-UTC主要用于恒星時角的計算,其預報誤差主要體現(xiàn)為天文經(jīng)度的誤差。天文經(jīng)度觀測值計算如為
λ=α+n+m{+δu+δr}+δa-Sg
(8)
式(8)中,Sg是表面時刻(XE+XW)/2時的格林尼治恒星時,微分式(8)得
dλ=dSg
(9)
即日常預報誤差對天文經(jīng)度的影響為
(10)
微分式(4)得:
(11)
根據(jù)目前天文定向常采用觀測北極星(δ≈90°)來實現(xiàn),則tanδ≈+∞, 即dAN≈0, 因此日常預報誤差對天文方位角基本上也不產(chǎn)生影響。
4.2 極移預報誤差
根據(jù)式(6)及式(7)并結合短期預報誤差情況,可得短期EOP預報值能夠滿足一等天文測量精度要求。我國地域遼闊,緯度范圍3°51′N~53°33′N,經(jīng)度范圍73°33′E~135°05′E,由式(6)及式(7)知極移誤差引起的天文解算影響關于緯度成增函數(shù),故而取漠河(121°10′E,53°33′N)為例,分析2010年、2011年和2012年極移的365 d預報誤差對天文定位定向的影響。
表5中給出了 1 a預報對天文定位定向影響的最大值、平均值和標準差,從表中可以看出,2010年和2012年的經(jīng)度誤差約為0.045″,超出了一等天文測量中經(jīng)度±0.02″的要求。
針對1 a預報無法滿足一等天文定位精度要求的情況,分別分析了1 a預報、半年預報以及一個季度預報對天文定位的影響,如表6所示。由表可知,1 a預報引起的經(jīng)度誤差的標準差在0.045″左右,可以滿足三等天文測量經(jīng)度±0.08″的要求;半年預報引起的經(jīng)度誤差的標準差在0.02″左右,可以滿足二等天文測量經(jīng)度±0.04″的要求;一個季度預報引起的經(jīng)度誤差的標準差不超過0.02″,可以滿足一等天文測量經(jīng)度±0.02″的要求。
表5 1 a預報對天文測量影響
表6 經(jīng)度誤差分析
本文綜合各參考文獻所述關于EOP參數(shù)預報分析的方法和結論,借鑒其中提出的比較好的分析方法并結合自身所做工作,對EOP參數(shù)的預報從長期和短期兩方面進行了詳細深入的分析,分別得到了1 d、7 d和1 a的預報分析結果。結合平時工程運用情況,將EOP預報誤差對天文定位定向的影響按照國家天文測量標準進行分門別類的分析,得到以下有益結論:
(1)1 d預報誤差精度較高,且穩(wěn)定性好,浮動較小,并以零值為中心呈現(xiàn)不規(guī)則分布。極移分量xp的誤差絕對值基本上在1mas以內(nèi),極移分量yp的誤差絕對值大都在0.5mas以內(nèi),UT1-UTC誤差絕對值都在0.2mas內(nèi)。7d預報誤差隨著天數(shù)增加而增大,極移分量xp的誤差絕對值在10mas以內(nèi),極移分量yp的誤差絕對值都在5mas以內(nèi),UT1-UTC誤差絕對值大都在1mas內(nèi)。
(2)1a的EOP誤差值存在一定的周期性。xp全年的預報誤差在67mas以內(nèi),yp全年的預報誤差在80mas以內(nèi),UT1-UTC全年的預報誤差在62mas以內(nèi),隨著年份推移,整體預報準確性有所提高,但不穩(wěn)定。
(3)日常預報(UT1-UTC)誤差對天文緯度沒有影響,對天文方位角基本不產(chǎn)生影響。1d日常預報誤差對經(jīng)度的最大影響約為0.207 5mas,一周日常預報誤差對經(jīng)度的最大影響約為1.223 8mas,1a日常預報誤差對經(jīng)度的最大影響約為62.046mas,都能夠滿足一等天文測量精度要求,可在解算時直接使用日常預報1a的預報值。
(4)極移預報誤差對緯度誤差不超過±0.03″,方位角誤差的標準差在0.05″左右,能滿足一等天文測量精度要求。1a極移預報引起的經(jīng)度誤差的標準差在0.045″左右,可以滿足三等天文測量要求;0.5a極移預報引起的經(jīng)度誤差的標準差在0.02″左右,可以滿足二等天文測量要求;一個季度極移預報引起的經(jīng)度誤差的標準差不超過0.02″,可以滿足一等天文測量要求。因此,進行一等天文定向解算時可直接使用365d的極移預報值,進行一等天文定位解算需對極移預報值進行季度更新。
綜合本文的分析,在進行天文定向定位時,各EOP參數(shù)預報誤差所產(chǎn)生的影響不可忽略,進行天文解算時必須加以考慮,必要時及時進行預報數(shù)據(jù)的更新。
[1]MCCARTHYDD,PETITC.IERSConventions(2003)[M].FrankfurtamMain,Germany:VerlagdesBundesamtesfürKartographieundGeod?sie,2004.
[2] 劉佳成,朱紫.2000年以來國際天文學聯(lián)合會(IAU)關于基本天文學的決議及其應用[J].天文學進展,2012,30(4):411-437.
[3] 劉佳成.新天文參考系若干問題的研究[D].南京:南京大學,2012.
[4] 夏一飛,金文敬.新參考系的引入對天體測量學的影響[J].天文學進展,2004,22(3):200-208.
[5] 張超,鄭勇,夏治國.多星近似等高法同時測定天文經(jīng)緯度的研究[J].測繪通報,2006(11):3-5+15.
[6] 張捍衛(wèi),許厚澤,王愛生.天文經(jīng)緯度和天文方位角測定的基本原理[J].測繪科學,2006,31(4):157-160+10.
[7] 張超.基于電子經(jīng)緯儀的天文測量系統(tǒng)及應用研究[D].鄭州:解放軍信息工程大學,2009.
[8] 詹銀虎.基于自然天體的快速定向理論及技術研究[D].鄭州:信息工程大學,2012.
[9] 王永海.CGCS2000的啟用對天文定位定向的影響研究[D].鄭州:解放軍信息工程大學,2013.
[10]張衛(wèi)星,劉萬科,龔曉穎.EOP預報誤差對自主定軌結果影響分析[J].大地測量與地球動力學,2011,31(5):106-110.
[11]柳文明,李崢嶸,劉文祥,等.EOP預報誤差對導航衛(wèi)星軌道預報的影響分析[J].全球定位系統(tǒng),2009(6):17-22+27.
[12]GB/T17943—2000,大地天文測量規(guī)范[S].
[13]劉彩璋,黃杰安,戴四軍.新天文常數(shù)、系統(tǒng)對天文經(jīng)緯度和方位角的影響[J].武漢測繪科技大學學報,1988,13(3):47-55.
[14]張捍衛(wèi),鄭勇,馬高峰.GCRS與ITRS之間的坐標轉換研究[J].大地測量與地球動力學,2011,31(1):63-67.
Influence of Earth Orientation Parameter on Celestial Fix and Astronomical Orientation
LIAN Yue-yong1,ZHANG Chao1,XIE Zong-te2
(1.College of Navigation and Aerospace Engineering,Information Engineering University,Zhengzhou 450001,China;2.Troops 75711,Guangzhou 510000,China)
The prediction of earth orientation parameters (EOP) can be used in the celestial fix and astronomical orientation,it can achieve the translation of conventional terrestial system and geocentric celestial reference system,but the latest EOP parameters value can not reach real time when calculation,so study on the prediction error of EOP is helpful to improve the calculation accuracy of astronomy.Celestial fix and astronomical orientation principle was introduced in the paper,and the application of EOP parameters was studied,the accuracy of forecasted EOP parameters on short-term and long-term were analyzed respectively,studied the effects of EOP parameters error on astronomical positioning and orientation.The results show that,365_day daily forecasting value can be used in the calculation of latitude and azimuth for an astronomical measurement,but the calculation of longitude should update the prediction value each quarter.The astronomical navigation and astronomical measurement can use these conclusions for reference.
EOP;celestial fix;astronomical orientation;polar motion;prediction error
2014-05-08
國家自然科學基金(41174025;41174026;41174027)。
連月勇(1990),男,福建泉州人,碩士生,現(xiàn)主要從事天文大地測量和天文導航研究。
P
A
2095-4999(2015)-01-0041-06