山東基準站網實現區域參考框架維持的方法

張海平，徐彥田，趙 碩

山東基準站網實現區域參考框架維持的方法

張海平1，徐彥田2，趙 碩2

（1. 山東省國土測繪院，濟南 250102；2. 中國測繪科學研究院，北京 100036）

為了實現山東省級區域大地坐標參考框架動態維持，通過SDCORS速度場模型改正板塊運行，在此基礎上通過ITRF框架間Helmert轉換的方法轉換到CGCS2000，試驗證明山東區域的坐標參考框架維持精度優于2 cm，考慮到SDCORS高精度服務精度指標需要每年更新坐標。

2000國家大地坐標系；山東基準站網；國際地球參考框架；區域參考框架；速度場模型

0 引言

為了描述地球空間信息的幾何形態和時空分布需要一個與地球固連在一起，與地球共同旋轉且在空間隨地球運行的非慣性參考系統即地球參考系統。地球參考系統是理論體系，定義了原點、尺度和定向及其實現等系列理論算法等；坐標參考框架是參考系統的實現，由一系列固連在地球表面的具有任意歷元點位坐標和速度場的基準點組成[1-2]。

地球是一個運動的系統，自轉軸內部不斷變化，地殼也存在復雜的運動，無法找到相對靜止的基準點；并且隨著大地測量技術的不斷發展，地球參考系統和地球參考框架不斷演變精化。目前，以國際地球自轉與參考系統服務（international Earth rotation and reference systems service,IERS）為代表的國際組織致力于地球參考框架的不斷完善，國際地球參考框架（international terrestrial reference frame,ITRF）是當前理論最完善、實現精度最高的全球參考框架，為其他全球和區域參考框架提供基準。1998年后，IERS先后發布了ITRF1997、ITRF2000、ITRF2005、ITRF2008和ITRF2014。

隨著山東省衛星導航定位基準服務系統（Shandong continuously operating reference stations，SDCORS）的運行，建立了覆蓋全省區域的間隔70~90 km的基準站119個，可以連續實時的獲取高精度的觀測數據，根據基準站坐標的時間序列建立和維持山東省區域坐標參考框架。山東基準站坐標周期性解算是基于當前的ITRF框架和觀測歷元，為了保持一致的坐標基準需要建立不同空間不同歷元坐標的轉換關系，并且由于山東CORS系統采用統一的國家規定的2000國家大地坐標系（China geodetic coordinate system 2000, CGCS2000），由于CGCS2000基于ITRF97框架、參考時刻為2000.0，需要將基準站站坐標轉換到2000國家大地坐標系[3-5]。

1 ITRF框架轉換方法

不同的ITRF框架對應著不同的坐標參考系統的定義，2個ITRF框架間關系可以通過赫爾默特（Helmert）相似變換關系描述，涉及14個轉換參數：3個平移參數、3個旋轉參數和1個比例因子及其變化率7參數；假定基準站坐標是線性變化、且原點差異在幾百米、尺度和旋轉在1×10-5量級時，模型表示為

由于坐標框架的動態特性，2、1、、、為時間的函數，求偏導得

考慮坐標向量歷元時刻和0，轉換整理得

1.1 板塊運動改正

由于不同ITRF框架對應不同的歷元時刻，需要計算歷元間板塊運動引起的坐標變化值，即

1.2 轉換參數歷元歸算

表1 ITRF2014和ITRF2008到ITRF1997轉換參數

表2 ITRF2014和ITRF2008到ITRF1997轉換參數變化率

1.3 框架點轉換

經過板塊運動改正和轉換參數歷元歸算獲得了轉換歷元的坐標向量和轉換參數，通過Helmert變換關系計算框架點坐標，將不同框架轉換到ITRF1997公式為

式中ITRFyy 為某一框架，若轉到2000國家大地坐標系則=2000.0。

2 山東區域速度場模型

SDCORS系統2012年正式運行，省內和鄰省共享基準站共150個，采用相同的策略每年定期解算基準站精確坐標[9]，坐標解算中誤差優于1 mm。從中選擇符合條件的基準站110個[10]，根據基準站坐標（2012—2018年）7 a的時間序列忽略非線性運行計算線性速度（圖1）[11]，并聯合使用國家陸態網絡2001—2010年觀測數據建立的中國大陸板塊運動模型（CPM-CGCS2000）[12]，采用反距離加權法建立山東省區域的15′×15′的格網速度場（圖2），模型精度中誤差東西（）方向0.98 mm/a，南北（）方向0.78 mm/a。

圖1 SDCORS速度場

圖2 山東格網速度場

3 試驗分析

采用山東省范圍內11個基準站2019年2月18日至2019年2月24日7 d的觀測數據，基準站分布如圖3。觀測數據采樣率30 Hz，采用相同的數據處理策略計算基準站觀測歷元時刻為2019.135 a ITRF2014框架的精確坐標，采用上述的ITRF框架間轉換方法轉化到CGCS2000下，并與已有的CGCS2000差值比較分析。

圖3 測站分布圖

表2統計了11個基準站觀測時間為2000.0~2019.135 a期間的板塊運動改正值，可以看出南北（）方向約為22 cm，每年向南緩慢移動約1.15 cm；東西（）方向約為57 cm，每年向東移動約為3 cm，與文獻[13]板塊趨勢一致。山東區域板塊運動引起的平面變化量約為3.25 cm，因此高精度坐標框架維持需要考慮板塊運動，并且考慮SDCORS基準的服務的精度指標需要1 a更新一次坐標值[14]。

表3 歷元時間為2 000.0~2 019.135 a期間的板塊運動改正值 　　 cm

表3統計了11個基準站在在第2000.0個歷元時，從ITRF2014轉換到ITRF1997的改正值。從表3可以看出：南北（）方向約為2.6 cm；東西（）方向約為0.4 cm。由于ITRF框架實現越來越完善精確，框架間的差異越來越小，考慮山東坐標框架的高精度維持，需要計算ITRF框架間差異量。

表3 第2000.0個歷元時ITRF2014和ITRF1997間差值 　　 cm

圖4顯示了11個基準站轉換到CGCS2000的坐標值和已有的坐標值間的差異，可以看出南北（）方向和東西（）方向方向偏差基本都在2 cm內，大部分在1 cm以內，顯然山東省區域的坐標框架維護精度在2 cm內，實現了全省范圍坐標基準空間上和時間上的統一。

圖4 基準站轉換精度

4 結束語

山東省區域板塊運動導致坐標每年位移約3 cm，是框架維護的主要改正量，綜合SDCORS服務的精度要求，SDCORS坐標需要每年更新一次滿足高精度和大比例尺測圖的需求。SDCORS衛星觀測數據經過相同的處理策略，經過板塊運動的改正和ITRF框架的轉換能夠實現全省坐標基準的統一，建立不同空間上和時間上的轉換關系，維護精度小于2 cm，高精度的坐標框架為自然資源監測調查分析提供時間上和空間上一致的坐標基準。

[1] 程鵬飛, 成英燕, 秘金鐘, 等. 國家大地坐標系建立的理論與實踐[M]. 北京: 測繪出版社, 2017: 1-2.

[2] 姜衛平, 馬一方, 鄧連生, 等. 毫米級地球參考框架的建立方法與展望[J]. 測繪地理信息, 2016, 41(4): 1-6.

[3] 成英燕, 黨亞民, 秘金鐘, 等. CGCS2000框架維持方法分析[J]. 武漢大學學報(信息科學版), 2017, 42(4): 543-549.

[4] 張西光, 呂志平. 論地球參考框架的維持[J]. 測繪通報, 2009(5): 1-4.

[5] 夏俊晨. 區域性參考框架的穩定性分析方法研究[D]. 上海: 華東師范大學, 2015.

[6] 彭小強, 高井祥, 王堅. WGS84和CGCS2000坐標轉換研究[J]. 大地測量與地球動力學, 2015, 35(2): 219-221.

[7] 徐世依. 北斗參考框架建立的理論與實驗研究[D]. 鄭州: 信息工程大學, 2017.

[8] 施闖, 魏娜, 李敏, 等. 利用北斗系統建立和維持國家大地坐標參考框架的方法研究[J]. 武漢大學學報(信息科學版), 2017, 42(11): 1635-1643.

[9] 程鵬飛, 成英燕. 基于GNSS的CGCS2000數據處理技術綜述[J]. 武漢大學學報(信息科學版), 2018, 43(12): 2071-2078.

[10] 國家測繪地理信息局. 大地測量控制點坐標轉換技術規范: CH/T 2014-2016[S]. 北京: 測繪出版社, 2017: 5-6.

[11] 孟國杰, 顧國華, 王武星. GPS測站位移速率的解算方法[J]. 大地測量與地球動力學, 2003, 23(3): 21-25.

[12] 程鵬飛, 成英燕, 秘金鐘, 等. CGCS2000板塊模型構建[J]. 測繪學報, 2013, 42(2): 159-167.

[13] 王東振, 趙斌, 余建勝, 等. 大規模GPS揭示中國現今地殼構造形變特征[J]. 地震工程學報, 2017, 39(3): 521-526, 544.

[14] 劉文建, 彭冬林. CORS坐標參考框架更新探討[J]. 測繪通報, 2017(2): 70-74.

An algorithm for maintenance of regional terrestrial reference frame base on SDCORS

ZHANG Haiping1, XU Yantian2, ZHAO Shuo2

（1. Shandong Provincial Institute of Land Surveying and Mapping, Jinan 250102, China；2. Chinese Academy of Surveying & Mapping, Beijing 100830, China）

In order to maintain of Shandong regional terrestrial coordinate frame based on SDCORS, the coordinate transformation model with space and time should be established. An algorithm of ITRF framework transformation is introduced and high-precision velocity model of Shandong Province is established which is about 3 cm per year; The Shandong regional terrestrial reference frame is maintained by velocity model and ITRF transformation parameters, the accuracy is better than 2 cm. So the coordinates of SDCORS need to be updated every year in view of the Accuracy of SDCORS Service.

CGCS2000; SDCORS; ITRF; regional terrestrial reference frame; velocity model

P228

A

2095-4999(2019)03-0131-04

2018-12-28

中國測繪科學研究院基本科研業務費（AR1903）。

張海平（1977—），男，山東五蓮人，碩士，高級工程師，研究方向為衛星導航定位基準站數據高精度處理。

張海平，徐彥田，趙碩.山東基準站網實現區域參考框架維持的方法[J].導航定位學報,2019,7(3):131-134.（ZHANG Haiping, XU Yantian, ZHAO Shuo.An algorithm for maintenance of regional terrestrial reference frame base on SDCORS[J].Journal of Navigation and Positioning,2019,7(3):131-134.）

10.16547/j.cnki.10-1096.20190322.