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

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

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

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

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

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

0 引言

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

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

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

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

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

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

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

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

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

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

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

1.2 主要誤差源

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

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

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

2.1 觀測地點(diǎn)、儀器

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

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

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

2.2.1 經(jīng)緯度解算

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

2.2.2 UT0解算

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

3.1 一致性分析

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

表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é)果的變化，但不會(huì)使組內(nèi)單次觀測精度（限本文中所規(guī)定單次觀測相對組均值的標(biāo)準(zhǔn)差）出現(xiàn)大的波動(dòng)；而觀測環(huán)境的影響比較復(fù)雜，如大氣視寧度的變化，影響CCD拍攝底片的星像質(zhì)量，使隨機(jī)誤差增大。溫度和氣壓的變化，當(dāng)溫度及氣壓梯度在東西方向變化時(shí)主要對經(jīng)度方向產(chǎn)生影響，在南北方向變化時(shí)主要對緯度方向產(chǎn)生影響[16-17]。

3.3 單次觀測精度分析

3.4 單組觀測精度分析

圖7 UT0解算結(jié)果

4 結(jié)論

[1] 韓延本, 王紅旗, 尹志強(qiáng), 等. 建立華北及西北地區(qū)垂線變化試驗(yàn)監(jiān)測網(wǎng)的意義和可行性[J]. 地球物理學(xué)進(jìn)展. 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] 田立麗, 郭金運(yùn), 韓延本, 等. 我國的數(shù)字化天頂望遠(yuǎn)鏡樣機(jī)[J]. 科學(xué)通報(bào), 2014, 59(12): 1094-1099.

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

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

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

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

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

[10] 葉叔華. 徐家匯觀象臺(tái)天文測時(shí)的季節(jié)性變化[J]. 天文學(xué)報(bào), 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] 楊廷高, 馬雨剛. 由陜臺(tái)光電等高儀(I)觀測得到的820顆恒星的位置改正[J]. 陜西天文臺(tái)臺(tái)刊, 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ù)字化天頂望遠(yuǎn)鏡觀測圖像及數(shù)據(jù)處理[J]. 科學(xué)通報(bào), 2014, 59(2): 1100-1107

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

[16] 李德河. 大氣反常折射對星表和參考系的影響[J]. 陜西天文臺(tái)臺(tái)刊, 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ù)字天頂望遠(yuǎn)鏡的初步觀測結(jié)果分析[J]. 時(shí)間頻率學(xué)報(bào), 2021, 44(1): 33-44.

2020-04-25；

2020-06-19

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