驪山觀測站數(shù)字天頂望遠鏡的初步觀測結(jié)果分析

劉娜，高玉平，蔡宏兵，張鵬飛，王平利，李琳

驪山觀測站數(shù)字天頂望遠鏡的初步觀測結(jié)果分析

劉娜1,2,3，高玉平1,2,3，蔡宏兵1,2，張鵬飛1,2,3，王平利1,2,3，李琳1,2,3

（1. 中國科學(xué)院 國家授時中心，西安 710600；2. 中國科學(xué)院 時間頻率基準重點實驗室，西安 710600；3. 中國科學(xué)院大學(xué)，北京 100049）

數(shù)字天頂望遠鏡（DZT）是一種新型的光學(xué)測量裝置，測站大氣擾動的影響是主要的誤差源之一。通過驪山觀測站上兩臺儀器的并址觀測避免了觀測環(huán)境差異的影響，對儀器的初步觀測結(jié)果進行了分析。同步觀測試驗結(jié)果表明兩臺儀器的一致性較好，當測站觀測環(huán)境較理想時DZT單次觀測精度緯度方向為0.16″，經(jīng)度方向為0.18″；緯度及經(jīng)度方向20 min一組的單組觀測標準差為0.05″~0.06″;兩儀器觀測UT0組間標準差分別為3.8 ms和3.5 ms。

數(shù)字天頂望遠鏡；天文經(jīng)緯度；世界時；誤差分析

0 引言

世界時（universal time，UT1）是中國科學(xué)院國家國家授時中心時間服務(wù)內(nèi)容之一。國家授時中心擬采用數(shù)字天頂望遠鏡（digital zenith telescope，DZT）開展世界時的測量。DZT是應(yīng)國家授時中心的應(yīng)用需求由國家天文臺研制的天文光學(xué)測量裝置，儀器設(shè)計參考傳統(tǒng)照相天頂筒（photographic zenith tube，PZT）的觀測原理，通過觀測天頂附近的恒星實現(xiàn)地方恒星時（或天文經(jīng)緯度）的測量。DZT采用CCD（charge coupled device）技術(shù)，以高精度傾斜儀與平面鏡代替水銀盤，利用GPS時間信號控制曝光時間、標定曝光歷元，實現(xiàn)了自動觀測，消除人儀差的影響[1]，提高了觀測精度。同時采用天文照相較差測量代替觀測固定星組，最快34 s可完成一個單次觀測，大大提高了觀測效率。與傳統(tǒng)光學(xué)儀器相同的是，儀器誤差及觀測環(huán)境的影響是制約觀測精度提高的主要因素。

本文介紹了數(shù)字天頂望遠鏡的觀測原理，觀測系統(tǒng)組成及觀測過程，利用兩臺儀器并址、同步觀測數(shù)據(jù)，采用分組處理，對比分析方法，避免了觀測環(huán)境差異的影響，對數(shù)字天頂望遠鏡的初步觀測結(jié)果進行了分析。

1 數(shù)字天頂望遠鏡觀測原理

式（3）中，，可由觀測的UTC時刻結(jié)合IERS發(fā)布的參數(shù)計算，因此已知UT1可根據(jù)式（1）、（2）解算測站緯度及經(jīng)度坐標。

多臺站觀測UT1的過程中，對測站在緯度方向與經(jīng)度方向可建立如下觀測方程：

對比式（5）可得世界時UT1與UT0之間的關(guān)系：

1.1 觀測系統(tǒng)組成與觀測過程

DZT自動觀測系統(tǒng)由上位機、控制箱、望遠鏡組成，其中望遠鏡主要由鏡筒及精密轉(zhuǎn)臺組成[4]，系統(tǒng)構(gòu)成如圖2所示。其中控制箱讀取DZT的狀態(tài)并控制望遠鏡工作，同時上位機與控制箱通過串口進行通訊，實時控制望遠鏡工作，可進行儀器調(diào)平、旋轉(zhuǎn)、觀測等操作。

轉(zhuǎn)身觀測過程中由傾斜儀與2個置平電機組成高精度置平系統(tǒng)，可快速將儀器水平狀態(tài)調(diào)整至2″以內(nèi)；置平后由GPS時間信號控制CCD快門進行曝光，系統(tǒng)自動記錄曝光時刻、傾斜儀讀數(shù)及CCD圖像；再由定位電機驅(qū)動望遠鏡進行180°轉(zhuǎn)身，轉(zhuǎn)身后讀取傾斜讀數(shù)，若傾斜儀相對測量結(jié)果在2″以內(nèi)則重復(fù)上述觀測與數(shù)據(jù)采集過程，若超過2″則由精置平系統(tǒng)進行水平調(diào)整然后觀測，至此0°和180°位的兩幅CCD底片、曝光時刻、相對傾斜測量結(jié)果構(gòu)成1個單次觀測，為了提高觀測效率，采用180°—180°—0°—0°的觀測模式，觀測過程如圖3所示。實際觀測中，由于DZT的視場較小云層遮擋易使拍攝的底片上無星或觀測星數(shù)過少而無法解算，加之觀測環(huán)境、數(shù)據(jù)處理等因素的影響，解算獲取的觀測值個數(shù)要小于相機拍攝的圖像對數(shù)，觀測值會出現(xiàn)中斷、不連續(xù)的情況。本文中的單次觀測數(shù)是指經(jīng)數(shù)據(jù)解算獲得的觀測值的個數(shù)。

圖2 數(shù)字天頂望遠鏡觀測系統(tǒng)組成

圖3 數(shù)字天頂望遠鏡觀測流程

1.2 主要誤差源

DZT為地面光學(xué)測地儀器，通過夜晚觀測恒星實現(xiàn)鉛垂線在天球背景中的方向測量，其特殊的180°轉(zhuǎn)身觀測消除儀器準直差、傾斜儀零點差的影響。除此之外，影響觀測精度的因素有儀器誤差、星表位置及自行誤差、數(shù)據(jù)處理過程中的誤差以及觀測環(huán)境的影響。其中儀器誤差主要有：儀器旋轉(zhuǎn)平臺及軸系誤差，CCD相機與高精度傾斜儀安裝方位角誤差，系統(tǒng)時延、傾斜儀測量誤差等；數(shù)據(jù)處理過程中的誤差主要是星像量度誤差；此外觀測環(huán)境中溫度、氣壓、風[5]、光等因素的變化，不僅會對大氣視寧度、天光背景產(chǎn)生影響，使星像產(chǎn)生閃爍、抖動，降低CCD拍攝底片的星像質(zhì)量[6-8]，還會引起儀器誤差改變[9-10]，如振動影響高精度電子傾斜儀觀測精度[11]，溫度變化會影響望遠鏡的焦距。

2 試觀測與數(shù)據(jù)處理

DZT的誤差源中觀測環(huán)境的影響不僅與測站環(huán)境有關(guān)還和觀測時段有關(guān)。星表位置誤差[12]也與測站緯度、觀測時段（DZT為非跟蹤觀測，一天內(nèi)觀測時間不同，所觀測的赤經(jīng)區(qū)不同）有關(guān)。因此并址、同步觀測，對觀測數(shù)據(jù)分組處理，能夠避免上述兩項誤差對觀測結(jié)果一致性分析的影響。

2.1 觀測地點、儀器

為了避免環(huán)境差異及星表誤差的影響，2017年夏季在驪山測站（原陜西天文臺子午環(huán)觀測室，緯度34.356 186 1°，經(jīng)度109.206 909 2°）進行了并址觀測實驗。觀測室內(nèi)有固建于同一觀測基墩的4個80 cm見方的小型觀測基墩，儀器DZT-I、DZT-II分別安裝在南北向相距約1.5 m的兩個基墩上，由觀測位置引起的緯度差約0.05″。觀測所用兩臺儀器均為折射式天頂望遠鏡，望遠鏡口徑200 mm，其中DZT-I安裝CCD為ALTA F9000，視場約1.3°，DZT-II安裝 CCD為ALTA U9，視場約1.73°，由此會引起兩臺儀器單次觀測星數(shù)的差異。兩儀器CCD快門時延均已校準，觀測中同時開機預(yù)熱，同步觀測。

2.2 數(shù)據(jù)處理

由數(shù)字天頂望遠鏡的觀測原理可知，已知ERP參數(shù)可實現(xiàn)測站經(jīng)緯度坐標的測量，已知測站坐標單臺站可解算世界時的初值UT0，具體解算方法如下。

2.2.1 經(jīng)緯度解算

為了分析觀測過程中兩臺儀器觀測精度的變化，以20 min為窗口進行滑動平均，20 min內(nèi)單次觀測值相對于組內(nèi)觀測平均值的均方差為該組單次觀測精度，如圖5中散點所示。

2.2.2 UT0解算

3 觀測結(jié)果與分析

3.1 一致性分析

兩臺儀器觀測的一致性不僅受到儀器誤差的影響、還受觀測環(huán)境的影響。圖4所示為DZT-I、DZT-II緯度及經(jīng)度觀測一致性及觀測的平均單次觀測星數(shù)，可見DZT-I每晚的平均單次觀測星數(shù)明顯小于DZT-II。經(jīng)度方向觀測值的一致性與平均單次觀測星數(shù)相關(guān)性較大，當DZT-I平均單次觀測星數(shù)明顯少于50顆時，兩臺儀器經(jīng)度方向觀測的一致性變差。其中經(jīng)度觀測值差異較大的6月1日、6月16日、6月26日以及7月20日，觀測值一致性分別為0.11″、0.12″、0.09″以及0.09″，對應(yīng)的DZT-1平均單次觀測星數(shù)分別為34顆、43顆、42顆與48顆，遠小于同期DZT-II的平均單次觀測星數(shù)95顆、122顆、122顆與130顆，也低于其8晚平均單次觀測星數(shù)53顆。當DZT-I單次平均觀測星數(shù)達到50顆以上時，如7月9日（58顆），7月11日（68顆），7月19日（84顆）兩儀器經(jīng)度一致性分別為0.06″，0.02″，0.03″。數(shù)字天頂望遠鏡測定經(jīng)緯度坐標是由CCD底片上拍攝星像的圖像位置與利用星表計算的參考星的視位置進行匹配，建立底片模型經(jīng)內(nèi)插解算得到，因此底片上拍攝星像的個數(shù)及分布會影響觀測精度，DZT-I在6月1日、6月16日、6月26日以及7月20日平均單次觀測星數(shù)偏少是影響兩臺儀器觀測一致性的一項重要因素。

表1 DZT-I、DZT-II同步觀測結(jié)果（扣除基墩位置差異的影響）

圖4 DZT-I、DZT-II觀測一致性與平均單次觀測星數(shù)

表2 DZT-I、DZT-II 7月9日觀測星數(shù)均大于50顆，連續(xù)同步觀測結(jié)果（扣除基墩位置差異的影響）

3.2 誤差分析

并址觀測的公共誤差源中星表位置誤差為系統(tǒng)誤差，引起觀測結(jié)果的變化，但不會使組內(nèi)單次觀測精度（限本文中所規(guī)定單次觀測相對組均值的標準差）出現(xiàn)大的波動；而觀測環(huán)境的影響比較復(fù)雜，如大氣視寧度的變化，影響CCD拍攝底片的星像質(zhì)量，使隨機誤差增大。溫度和氣壓的變化，當溫度及氣壓梯度在東西方向變化時主要對經(jīng)度方向產(chǎn)生影響，在南北方向變化時主要對緯度方向產(chǎn)生影響[16-17]。

3.3 單次觀測精度分析

3.4 單組觀測精度分析

圖7 UT0解算結(jié)果

4 結(jié)論

[1] 韓延本, 王紅旗, 尹志強, 等. 建立華北及西北地區(qū)垂線變化試驗監(jiān)測網(wǎng)的意義和可行性[J]. 地球物理學(xué)進展. 2012, 27(4): 1287-1293

[2] PETIT G, LUZUM B. The IERS conventions (2010)[J]. IERS Technical Note, 2010(36): 1-95.

[3] 趙銘. 天體測量學(xué)導(dǎo)論[M]. 北京: 中國科學(xué)技術(shù)出版社, 2006.

[4] 田立麗, 郭金運, 韓延本, 等. 我國的數(shù)字化天頂望遠鏡樣機[J]. 科學(xué)通報, 2014, 59(12): 1094-1099.

[5] 劉次沅. 陜西天文臺光電中星儀1972—1976年測時結(jié)果的風效應(yīng)[J]. 陜西天文臺臺刊, 1978, 1(1): 10-13.

[6] 譚徽松. 選址中的大氣視寧度測量[J]. 天文學(xué)進展, 1992, 10(1): 48-55.

[7] 高布錫. 反常折射引起的地面天文觀測殘差的周年變化[J]. 天文學(xué)報, 1994(1): 93-97.

[8] 冒蔚, 彭青玉, 王曉彬. 研究大氣反常折射的迫切性[J]. 天文研究與技術(shù), 1997(s1): 101-105.

[9] 葉叔華. 徐家匯觀象臺中星儀方位差的變化[J]. 天文學(xué)報, 1955, 3(2): 197-216.

[10] 葉叔華. 徐家匯觀象臺天文測時的季節(jié)性變化[J]. 天文學(xué)報, 1957, 5(2): 281-291.

[11] HIRT C. Entwicklung und erprobung eines digitalen zenitkamerasystems für die hochpr?zise lotabweichungsbestimmung[D]. Munich, Bayern: Technische Universit?t München, 2004.

[12] 楊廷高, 馬雨剛. 由陜臺光電等高儀(I)觀測得到的820顆恒星的位置改正[J]. 陜西天文臺臺刊, 1983, 6(1): 75-109.

[13] PENG Q Y, HAN Y B, ZHANG C L, et al. Image-processing techniques in precisely measuring positions of Jupiter and its Galilean satellites[J]. Astronomy Astrophysics, 2003, 401(2): 773-779

[14] 王博, 田立麗, 王政, 等. 數(shù)字化天頂望遠鏡觀測圖像及數(shù)據(jù)處理[J]. 科學(xué)通報, 2014, 59(2): 1100-1107

[15] 季凱帆, 王鋒. CCD圖像的二維修正矩定心方法[J]. 天文學(xué)報, 1996, 37(1): 85-90.

[16] 李德河. 大氣反常折射對星表和參考系的影響[J]. 陜西天文臺臺刊, 1998, 21(1): 59-64.

[17] HIRT C. Monitoring and analysis of anomalous refraction using a digital zenith camera system[J]. Astronomy & Astrophysics, 2006, 459(1): 283-290.

Preliminary observation results of digital zenith telescope at Li-Shan observatory

LIU Na1,2,3, GAO Yu-ping1,2,3, CAI Hong-bing1,2, ZHANG Peng-fei1,2,3,WANG Ping-li1,2,3, LI Lin1,2,3

(1. National Time Service Center, Chinese Academy of Sciences, Xi’an 710600, China;2. Key Laboratory of Time and Frequency Primary Standards, Chinese Academy of Sciences, Xi’an 710600, China;3. University of Chinese Academy of Sciences, Beijing 100049, China)

Digital Zenith Telescope (DZT) is a new type of optical instrument. One of the main errors of DZT is the atmospheric turbulence at the observational station. In this paper, the effects of observation environments can be effectively cancelled by co-location observations of two instruments at the Li-Shan station. The performance of the two instruments are evaluated. The simultaneous observations of two instruments show good consistency. When the observation environment of the station is well, the standard deviation of the derived station position using one scan data is about 0.16″ in the latitude direction and 0.18″ in the longitude direction. Combining 20-min observations, yields precision of 0.05″ to 0.06″ for measurements of station’s latitude and longitude. The standard deviations of UT0 derived by the two instruments are 3.8 ms and 3.5 ms, respectively.

digital zenith telescope; astronomical latitude and longitude; universal time; error analysis

10.13875/j.issn.1674-0637.2021-01-0033-12

劉娜, 高玉平, 蔡宏兵, 等. 驪山觀測站數(shù)字天頂望遠鏡的初步觀測結(jié)果分析[J]. 時間頻率學(xué)報, 2021, 44(1): 33-44.

2020-04-25；

2020-06-19

國家自然科學(xué)基金資助項目（91736207）；中國科學(xué)院國家授時中心“青年創(chuàng)新人才”資助項目（國授發(fā)字〔2017〕48號）；科技大數(shù)據(jù)湖北省重點實驗室開放基金資助項目（20KF011015）