劉志強，劉全海，王解先，謝中華

（1.同濟大學 測量與國土信息工程系，上海 200092；2.河海大學 地球科學與工程學院，江蘇 南京 210098；3.常州市測繪院，江蘇 常州 213002；4.現(xiàn)代工程測量國家測繪局重點實驗室，上海 200092）

隨著城市信息化建設進程的不斷推進，國內(nèi)的一些城市近年來紛紛開展連續(xù)運行衛(wèi)星定位服務系統(tǒng)CORS（Continuously Operating Reference Stations）的建設。城市CORS是在傳統(tǒng)GNSS（Global Navigation Satellite System）差分定位的基礎上快速發(fā)展起來的一種基于多基準站的區(qū)域衛(wèi)星定位增強系統(tǒng)。它不僅是動態(tài)的、連續(xù)的空間數(shù)據(jù)參考框架，同時也是快速、高精度獲取空間數(shù)據(jù)和地理特征的重要手段。由于受板塊運動、人類活動等因素影響，CORS系統(tǒng)建成后各參考站會隨著時間的變化產(chǎn)生不同程度的三維位移變化。對CORS參考站的位移變化規(guī)律進行分析研究，既可解決CORS運營維護者普遍關心的系統(tǒng)穩(wěn)定性問題，又能夠為城市區(qū)域的地表變形提供一種高精度、有效的監(jiān)測手段。近年來，國內(nèi)外一些學者利用區(qū)域CORS開展了一系列有關參考站穩(wěn)定性分析、板塊構(gòu)造運動、地表變形監(jiān)測的研究工作［1－5］。在對CORS數(shù)據(jù)處理的過程中，不同學者對參考框架的選擇、變形監(jiān)測基準合理性等問題的認識仍存在不同程度的差異，針對不同研究對象得到的參考站位移和區(qū)域地表變形速度場的結(jié)果亦不盡相同。本文以常州市CZCORS系統(tǒng)為例，對比分析不同參考框架和變形監(jiān)測基準對參考站位移監(jiān)測結(jié)果的影響，并探討利用城市CORS進行區(qū)域地表變形監(jiān)測與分析的方法。

1 CORS觀測數(shù)據(jù)與基線解算

常州市CZCORS系統(tǒng)由CZSM、JTAN、LYNG、WJIN、XBEI共5個連續(xù)運行GNSS參考站組成，參考站間平均距離43.9km。各參考站均采用高精度雙頻GNSS接收機（LEICA GRX1200GGPRO），并采用統(tǒng)一的雙頻大地型天線（LEIAX1202GG）。各參考站的GNSS天線均安置于基礎相對穩(wěn)定的建筑物頂部，周圍觀測天空開闊，受多路徑效應影響較小。CZCORS各參考站空間分布合理，能較好地覆蓋常州市各行政區(qū)劃，可作為區(qū)域地表變形的監(jiān)測網(wǎng)。

本文選取2009年10月至2011年10月期間24個月的CZCORS參考站的觀測數(shù)據(jù)進行分析處理，且設定監(jiān)測周期為6個月。為增強監(jiān)測網(wǎng)基線處理的可靠性，每期監(jiān)測數(shù)據(jù)分別取連續(xù)7d的數(shù)據(jù)進行基線解算，得到連續(xù)7d的單天解作為本期觀測的重復基線?；€解算采用GAMIT軟件（10.4版），在進行基線處理時，下載BJFS、KUNM、SHAO、WUHN、TWTF共5個IGS站的同步觀測數(shù)據(jù)進行聯(lián)合解算。聯(lián)測的IGS（International GNSS Service）站近似坐標為ITRF2005參考框架、瞬時歷元坐標，并設置各IGS站東西向、南北向及徑向 先 驗 坐 標 約 束 為0.005m、0.005m和0.010m。其中，固定SHAO站作為起算點，其余IGS站作松弛約束，可得到基線單天無約束解。在采用IGS精密星歷和高精度起算點時，其解算長基線的相對精度能達到10－9量級，能滿足高精度變形監(jiān)測網(wǎng)的要求。

2 參考框架與監(jiān)測基準的選擇

2.1 ITRF框架

在進行變形監(jiān)測數(shù)據(jù)的處理和分析時，首先要將GNSS技術獲得的監(jiān)測點坐標納入到合理、統(tǒng)一的參考框架中。國際地球參考框架ITRF（Internatial Terrestrial Reference Frame）是利用GNSS、VLBI、SLR、LLR及DORIS等多種空間技術，由全球測站觀測資料成果推算得到的四維地心坐標參考框架。迄今為止，IERS已發(fā)布11個版本的ITRF框架。目前最新、精度最高的坐標參考框架是ITRF2008，于2010－05－28日 發(fā) 布。相 對 于ITRF2005，ITRF2008在原點、尺度等框架定義參數(shù)方面有較大的改進和修正，其參考站坐標和速率的精度也均有明顯提高［6］。每次ITRF的更新，都意味著對地殼運動認識的不斷深入。ITRF的實現(xiàn)能夠充分顧及到地殼板塊運動規(guī)律，并且具有高精度、穩(wěn)定性強的特點，因而可以作為較為理想的大尺度范圍內(nèi)地表變形分析的動態(tài)坐標參考基準。

由于CORS得到的區(qū)域地表變形監(jiān)測數(shù)據(jù)時間周期一般較長，期間通常會涉及到不同時期的ITRF坐標參考框架。由于不同的ITRF框架的定義不同，這就導致了ITRF框架間會存在系統(tǒng)性差異。在求取各期監(jiān)測網(wǎng)點位坐標成果時，應歸算到統(tǒng)一的ITRF框架下，由此消除因坐標參考基準差異引起的位置、尺度、方位和時間演變的系統(tǒng)性誤差。本文分析的CZCORS數(shù)據(jù)涉及ITRF2005和ITRF2008兩種坐標參考框架。二者之間的轉(zhuǎn)換可通過基于歐氏空間的相似變換實現(xiàn)。

式中：T，˙T為轉(zhuǎn)換參數(shù)及其變化率；D，˙D為尺度因子及其變化率；R，˙R分別為旋轉(zhuǎn)矩陣及其變化率。ITRF2005和ITRF2008在2005.0歷元的14個轉(zhuǎn)換參數(shù)具體值可參見文獻［6］。

2.2 監(jiān)測網(wǎng)平差基準

在對變形監(jiān)測網(wǎng)進行數(shù)據(jù)處理和分析時，一般采用經(jīng)典的所謂固定參考基準。然而，對于大尺度范圍內(nèi)的CORS監(jiān)測網(wǎng)而言，則需將參考基準納入到國際地球參考框架ITRF中，這時就應充分顧及基準的時空變化規(guī)律，將其視為動態(tài)參考基準。從監(jiān)測網(wǎng)的平差方式來看，若僅固定單個IGS站進行三維無約束平差，則不會使CORS網(wǎng)形產(chǎn)生畸變。由此得到的是各期監(jiān)測網(wǎng)點的相對ITRF坐標，其僅相對于所選的IGS測站具有較高的精度。但在由不同期監(jiān)測成果計算點位位移時，要考慮所選IGS站的坐標及速率的誤差影響，選取固定不同的IGS站所得到的位移場結(jié)果可能不同。若固定多個IGS站進行約束平差，則可將監(jiān)測網(wǎng)納入到由所選的IGS站組的坐標及其速率所確定的坐標參考系統(tǒng)中，由此得到的是各期監(jiān)測網(wǎng)點的絕對ITRF坐標。由于約束平差的方式是將CORS網(wǎng)與所選IGS站組進行強制符合，理論上會引起監(jiān)測網(wǎng)形的畸變。這取決于所選IGS站組坐標及速率的精度及相互之間是否有較好的兼容性。

2.3 對比分析

CZCORS基線解算完成后，進行網(wǎng)平差得到的是對應ITRF框架下的三維空間直角坐標（X，Y，Z）。將2009－10的CZCORS網(wǎng)的平差結(jié)果作為首期成果，其后4期的平差結(jié)果與之相減得到各監(jiān)測點的三維位移變化值（dX，dY，dZ）。為直觀地分析點位變化情況，通常需要將地心坐標系中的點位位移結(jié)果轉(zhuǎn)換到站心坐標系中，得到平面（North、East）和高程（Up）分量的位移變化值（N，E，U）。各期位移（dN，dE，dU）除以時間間隔即可得到監(jiān)測點在平面和高程方向上的位移速率。地心坐標系和站心坐標系間的轉(zhuǎn)換關系如式

（3）所示［7］。

式中：B，L分別為測站點的大地緯度和大地經(jīng)度。

為比較采用不同ITRF框架、不同平差基準對監(jiān)測網(wǎng)點位移結(jié)果的影響，分別采用如下兩種方案對CZCORS監(jiān)測網(wǎng)進行處理：

方案1：在ITRF2005和ITRF2008框架下進行無約束平差（僅固定SHAO站），結(jié)果如表1所示；

表1 ITRF2005和ITRF2008框架各期無約束平差位移值mm

方案2：在ITRF2008框架下，分別進行無約束平差（僅固定SHAO站）和約束平差（固定5個IGS站，即BJFS、SHAO、KUNM、TWTF、WUHN），結(jié)果如表2所示。

表1中SHAO站的各期位移值由ITRF公布的已知值計算得到。從表1可以看出，分別采用ITRF2005和ITRF2008框架進行三維無約束平差，在不同框架下各CORS監(jiān)測站點位移的平面分量值與參考基準SHAO站的變化規(guī)律表現(xiàn)出高度一致性；而高程分量值的變化差異較大，這主要是由于城市CORS參考站與全球IGS站的標墩基礎構(gòu)造條件不同所引起的。由此可見，采用不同的ITRF框架進行無約束平差，各期監(jiān)測點的位移結(jié)果主要取決于所選取的IGS站在不同ITRF框架下的坐標及速率差異。在同一參考框架下，固定不同的IGS站作為參考基準，得到的監(jiān)測網(wǎng)點位移結(jié)果會有所不同。無約束平差結(jié)果得到的各期監(jiān)測網(wǎng)點坐標為相對坐標，其相對于作為參考基準的IGS站具有較高精度，而且監(jiān)測網(wǎng)形無畸變。

從表2可以看出，在ITRF2008框架下，僅固定SHAO站進行無約束平差，得到的BJFS、KUNM、TWTF、WUHN參考站的各期位移結(jié)果，與ITRF2008已知值計算出的結(jié)果相比均存在不同程度的差異。這反映出無約束平差結(jié)果所確定的尺度與ITRF2008框架的尺度存在較為明顯的差異。從CZCORS各監(jiān)測點的位移結(jié)果互差值來看，采用無約束和約束平差對各測站的影響值基本一致。這說明約束平差僅對CZCORS網(wǎng)有整體的平移效應，并未改變CZCORS各網(wǎng)點間的相對位置關系。不同期的約束平差結(jié)果對CZCORS網(wǎng)的平移程度略有不同，這種差異的量級取決于不同期基線解算的精度變化，在能保證基線解算結(jié)果精度和可靠性的情況下，不會對分析各期位移結(jié)果變化情況產(chǎn)生明顯影響。采用固定多個IGS站進行約束平差，可以得到各期監(jiān)測網(wǎng)點在ITRF框架下的絕對坐標結(jié)果，其所確定的尺度與ITRF2008框架本身的尺度較為接近。在現(xiàn)有GNSS基線解算所能達到的精度水平下，由無約束和約束平差結(jié)果計算得到的各期監(jiān)測網(wǎng)點間的相對位移變化結(jié)果并無顯著差異。

表2 ITRF2008框架各期無約束與約束平差位移結(jié)果互差 mm

3 位移結(jié)果與地表變形分析

在ITRF2008框架下，固定5個IGS站（BJFS、SHAO、KUNM、TWTF、WUHN）對CZCORS網(wǎng)進行約束平差，可得到各期監(jiān)測網(wǎng)點相對于2009年10月首期坐標成果的位移變化值，如表3所示。

表3 ITRF2008框架各期約束平差位移結(jié)果 mm

表3中5個IGS站的各期位移值由ITRF公布的已知值計算得到。CZCORS各期監(jiān)測點位移值平面分量（North，East）變化與已知的各IGS站變化規(guī)律有較好吻合。在ITRF2008框架下，CZCORS網(wǎng)有向東南方向整體運動的趨勢，其平面位移變化值及其運動方向與各IGS站對應的結(jié)果基本一致。由表3中CZCORS各監(jiān)測網(wǎng)點的位移變化值，可以計算出常州市區(qū)域地表水平位移速度場如圖1所示。

圖1 ITRF2008框架下水平位移速度場

從表3中給出的各測站垂直方向位移變化值可以看出，不同IGS站及各CZCORS參考站間存在較大差異。對IGS站而言，BJFS、TWTF存在明顯的抬升趨勢，而KUNM、SHAO及WUHN站在監(jiān)測周期內(nèi)的變化并不明顯。圖2為CZCORS各期監(jiān)測點垂直分量變化情況。其中，CZSM、JTAN站表現(xiàn)出一定的抬升趨勢，LYNG、XBEI站總體變化量不大，而WJIN站則表現(xiàn)出較為明顯的下沉趨勢。該監(jiān)測結(jié)果與城市二等幾何水準網(wǎng)平差計算結(jié)果基本一致。由于GNSS高程方向的測量精度一般要比水平方向低2～3倍，而且CZCORS各參考站均為屋頂標墩，季節(jié)溫度變化因素會一定程度上影響建筑物垂直方向上的位移，因此，僅由CZCORS得到的各參考站的垂直分量位移結(jié)果可靠性往往難以保證。將ITRF框架下測站點垂直方向上的相對位移計算結(jié)果與城市二等幾何水準成果進行相互驗證，可以為城市CORS各參考站高程穩(wěn)定性變化提供一種較為有效、可靠的檢驗方法。由于受城市CORS網(wǎng)參考站的數(shù)量及空間分布密度限制，CORS站的變形監(jiān)測成果尚不足以精確、可靠地反映城市大面積區(qū)域的地表垂直位移變化狀況，為此還可進一步結(jié)合該區(qū)域的差分干涉測量D－InSAR（Differential－InSAR）技術成果進行分析［8］。

圖2 ITRF2008框架下垂直位移變化值

4 結(jié) 論

本文通過對常州市CZCORS網(wǎng)監(jiān)測數(shù)據(jù)的處理，對比分析了不同ITRF參考框架和變形監(jiān)測基準對城市CORS參考站位移及區(qū)域地表變形監(jiān)測結(jié)果的影響，得到以下結(jié)論和建議：

1）在利用城市CORS網(wǎng)進行參考站穩(wěn)定性監(jiān)測與區(qū)域地表變形分析時，應采用動態(tài)參考基準，并將其納入到國際地球參考框架ITRF中。

2）采用不同的ITRF框架進行無約束平差，主要取決于所選取的IGS站在不同ITRF框架下的坐標及速率差異，而且得到的監(jiān)測成果為相對坐標，僅相對于作為參考基準的IGS站具有較高精度。

3）在統(tǒng)一的ITRF框架下，在現(xiàn)有GNSS基線處理所能達到的精度水平下，由無約束和約束平差結(jié)果計算得到的各期監(jiān)測網(wǎng)點間的相對位移變化結(jié)果并無顯著差異。在最新的ITRF2008框架下固定一定數(shù)量的IGS站進行約束平差，得到的監(jiān)測成果為絕對坐標，其所確定的尺度與ITRF2008框架本身的尺度較為接近。

4）利用城市CORS網(wǎng)的監(jiān)測成果，可以精確地建立城市區(qū)域地表水平位移速度場；若進一步結(jié)合精密水準測量、D－InSAR等監(jiān)測成果進行相互驗證，則可獲得精確的城市區(qū)域地表垂直位移速度場。

［1］黃立人.GPS觀測結(jié)果形變分析的參考框架及其合理性［J］.測繪學報，2001，30（1）：16－20.

［2］李江衛(wèi)，劉經(jīng)南，肖建華，等.CORS系統(tǒng)穩(wěn)定性監(jiān)測數(shù)據(jù)處理與分析研究［J］.武漢大學學報：信息科學版，2010，35（7）：825－829.

［3］薄萬舉，胡新康，董運洪，等.用GPS位移場進行中小區(qū)域變形分析方法探討［J］.大地測量與地球動力學，2010，30（3）：31－34.

［4］ZHANG KEFEI，HU YOUJIAN，LIU GANGJUN，et al.Deformation monitoring and analysis using Victorian regional CORS data［J］.Journal of Global Positioning Systems，2005，4（1）：129－138.

［5］劉經(jīng)南，施闖，許才軍，等.利用局域復測GPS網(wǎng)研究中國大陸現(xiàn)今地殼運動速度場［J］.武漢大學學報：信息科學版，2001，26（3）：189－195.

［6］ALTAMIMI Z，COLLILIEU X，MéTIVIER L.ITRF2008：an improved solution of the international terrestrial reference frame［J］.Journal of Geodesy，2011，85（8）：457－473.

［7］伍吉倉，鄧康偉，陳永奇.地心坐標系與站心坐標系中的速度轉(zhuǎn)換及誤差傳播［J］.大地測量與地球動力學，2005，25（3）：13－18.

［8］何敏，何秀鳳.利用D－InSAR技術監(jiān)測鹽城地區(qū)地表形變［J］.測繪通報，2010（11）：1－3.