李凱鋒，歐陽永忠，2，陸秀平，任來平

(1.天津海洋測繪研究所，天津300061;2.武漢大學(xué)測繪學(xué)院，湖北武漢430079)

玉樹地震對周圍衛(wèi)星定位連續(xù)運行站的影響分析

李凱鋒1，歐陽永忠1，2，陸秀平1，任來平1

(1.天津海洋測繪研究所，天津300061;2.武漢大學(xué)測繪學(xué)院，湖北武漢430079)

利用相同時段長度且年份相鄰的兩組連續(xù)運行站的數(shù)據(jù)，均勻地選取7個IGS站作為參考基準(zhǔn)，采用GAMIT進行數(shù)據(jù)處理。比較分析7個中國地殼運動觀測網(wǎng)絡(luò)站單日解坐標(biāo)的變化，得出玉樹地震對周圍衛(wèi)星定位連續(xù)運行站的位移有不同程度影響的結(jié)論，尤其是在高程方向，平均偏差達到16 mm，個別站點的最大偏差達到30 mm。通過對GPS速度場的計算分析，得出玉樹周圍各塊體的運動情況。

衛(wèi)星定位連續(xù)運行站;IGS站;GAMIT軟件;位移;速度場

一、引 言

2010年4月14日在玉樹發(fā)生了7.1級強烈地震。地震發(fā)生后如何準(zhǔn)確、快速地對地殼運動進行調(diào)查分析，對科學(xué)救災(zāi)、災(zāi)后重建以及衛(wèi)星定位連續(xù)運行站的建立和維持都具有重要意義。連續(xù)運行站作為高等級控制測量的基準(zhǔn)，其位置的變化及運動速度將影響大地測量成果的精度。因此，本文利用IGS站和中國地殼運動觀測網(wǎng)絡(luò)工程部分站的數(shù)據(jù)，以目前常用的GPS數(shù)據(jù)處理軟件GAMIT進行處理，分析此次玉樹地震對周圍衛(wèi)星定位連續(xù)運行站的影響。

二、周圍站點選取

在進行GPS數(shù)據(jù)處理時，為了得到準(zhǔn)確的國際地球參考框架下的坐標(biāo)，測站在近3年連續(xù)觀測，且網(wǎng)形合理，針對玉樹周圍網(wǎng)絡(luò)工程站點，選取國內(nèi)的BJFS、WUHN、KUNM和URUM 4個站點，以及國外的ULAB、HYDE和CUSV 3個站點，共7個IGS站作為參考基準(zhǔn)，計算中國地殼運動觀測網(wǎng)絡(luò)中的DLHA、XNIN、DXIN、YANC、LHAS、MMMZ和LUZH 7個連續(xù)運行站的單日解坐標(biāo)。為了分析玉樹地震前后周圍地殼的運動，選取沒有發(fā)生地震的2009年年積日自297至305(對應(yīng)為10月24日—11月1日)為期9 d的觀測數(shù)據(jù)，以及2010年年積日自100至108(對應(yīng)為4月10日—4月18日)為期9 d的觀測數(shù)據(jù)，數(shù)據(jù)采樣間隔為30 s，利用IGS最終精密星歷，數(shù)據(jù)處理采用GAMIT/GLOBK軟件，對兩組數(shù)據(jù)解算結(jié)果進行分析。所選站點分布如圖1所示。

圖1 連續(xù)觀測站站點分布

三、數(shù)據(jù)處理方案

為提高GPS定位解算精度，基線解算中對數(shù)據(jù)作如下設(shè)定:觀測數(shù)據(jù)類型(LC_HELP);衛(wèi)星截止高度角(15°);歷元間隔(30 s);對流層延遲(13個);IGS站點坐標(biāo)約束(NER:0.01 m，0.01 m，0.02 m);網(wǎng)絡(luò)工程站的坐標(biāo)約束(NER:50 m，50 m，50 m);衛(wèi)星軌道約束為松弛約束;電離層約束(0.0 m+8.00×10－6D);對流層模型(Saastamoninen模型);慣性系框架(J2000.0);光壓模型(BERNE);處理模式(RELAX解)。利用GLOBK進行時間序列分析和平差。

四、處理結(jié)果分析

GAMIT基線解算結(jié)果，主要通過均方根殘差、基線向量、基線重復(fù)率3項指標(biāo)進行評定?；€重復(fù)率和相對重復(fù)率的計算公式為

式(1)、式(2)中，Rl為基線向量的重復(fù)率;Rr為基線向量的相對重復(fù)率;n為基線單日解數(shù)目;Li為第i日的基線分量(或邊長);為單日解的基線分量(或邊長)的加權(quán)平均值;是基線分量中誤差。的計算公式為

由式(3)計算得出基線的NRMS值都在0.18 mm左右，基線解算質(zhì)量較高，然后計算各個站點的坐標(biāo)中誤差，采用GLOBK加入基線解的H文件進行聯(lián)合平差。

本文數(shù)據(jù)后處理以單天松弛解為準(zhǔn)觀測值完成的，對于玉樹周圍的中國地震觀測網(wǎng)絡(luò)工程基準(zhǔn)站而言，則通過緊約束地球自轉(zhuǎn)參數(shù)，利用7個均勻分布在周圍的參考站計算出網(wǎng)絡(luò)工程基準(zhǔn)站ITRF2005參考框架下的單日坐標(biāo)。各網(wǎng)絡(luò)工程站至玉樹震中的概況如表1所示。經(jīng)過坐標(biāo)轉(zhuǎn)換得到平面和高程方向的單日坐標(biāo)，單日坐標(biāo)解B、L、H方向的偏差值的時間序列如圖2和圖3所示，其中，橫軸代表年積日，縱軸為偏差值。

表1 各網(wǎng)絡(luò)工程站至玉樹震中的概況

圖2 衛(wèi)星定位連續(xù)運行站B、L、H方向的時間序列變化(2009年)

圖3 玉樹地震前后周圍衛(wèi)星定位連續(xù)運行站B、L、H方向的時間序列變化(2010年)

表2、表3和表4中2009年年積日自297至305為沒有發(fā)生地震的時間段，而2010年年積日自100至108中，年積日為104時發(fā)生玉樹地震。計算平均偏差時，本文用單日坐標(biāo)解減去測站的已知坐標(biāo)值。從圖2(a)、圖3(a)和表2可以看出，沒有發(fā)生地震時同一基準(zhǔn)站緯度方向的變化都在4.5 mm之內(nèi)，且變化較平緩;玉樹地震當(dāng)天，震中東北方向的DLHA、XNIN、DXIN、YANC 4個基準(zhǔn)站的緯度方向同時出現(xiàn)最低值，最大偏差達到7.9 mm，表明震中東北方向地殼往南偏移，其他兩個方向的基準(zhǔn)站變化不明顯。

表2 緯度方向位移 cm

表3 經(jīng)度方向位移 cm

從圖2(b)、圖3(b)和表3可以看出，沒有發(fā)生地震時同一基準(zhǔn)站經(jīng)度方向的變化都在4.5 mm之內(nèi)，且變化非常平緩;年積日為102當(dāng)天，震中東北方向的DLHA、XNIN和DXIN以及西南方向的LHAS 4個基準(zhǔn)站的經(jīng)度方向同時出現(xiàn)最大值，最大偏差達到5.5 mm，說明震中東北方向和西南方向地殼往東偏移，其他兩個方向的基準(zhǔn)站變化不明顯。

從圖2(c)、圖3(c)和表4可以看出，沒有發(fā)生地震時，同一基準(zhǔn)站高程方向的變化都在11 mm之內(nèi)，大部分偏差都小于6 mm，且變化較平緩。玉樹地震當(dāng)天，震中東北方向的DLHA、XNIN、DXIN、YANC 4個基準(zhǔn)站的高程方向同時出現(xiàn)最大值，最大偏差達到19.9 mm，表明震中東北方向地殼發(fā)生隆起，其他兩個方向的基準(zhǔn)站高程方向的偏差也在9 mm左右。玉樹地震后的第二天，由于余震的影響，玉樹東南和西南方向的LUZH、MMMZ、LHAS 3個基準(zhǔn)站的高程方向同時出現(xiàn)最低值，最大偏差達到30.0 mm，而玉樹東北方向的DLHA和XNIN兩個基準(zhǔn)站的高程方向同時出現(xiàn)最低值，偏差達到12.3 mm和6 mm，西南方向地殼比東南、東北方向地殼下沉更嚴(yán)重。

表4 高程方向位移 cm

以上分析是基于各個站點的B、L、H方向的變化情況，下面從GPS速度場分析地震對各個站點的影響。設(shè)X0為參考?xì)v元的坐標(biāo)，X為觀測歷元的坐標(biāo)，本文取時段中的周解坐標(biāo)，X˙代表運動速度，分別代表東向、北向和垂向速度;t0為參考?xì)v元，t為觀測歷元，本文取時段中的平均歷元。于是上述兩個時段的運動速度可表示為

式中，t02、t09是指參考?xì)v元為2002.8和2009.8。玉樹地震的GPS速度場計算結(jié)果如表5所示。

表5 玉樹周圍站點的GPS速率 cm/a

從表5可以看出，由兩個參考?xì)v元計算得出的速率有一定的差別，最大差值出現(xiàn)在LHAS站，差值為64.2 mm。圖4是從水平運動速度來分析各站點的運動方向。

圖4 玉樹周圍的GPS水平速度場

從圖4(a)可以看出，參考?xì)v元為2002.8時，整體水平速度向南，而且速率差別不大。從圖4(b)可以看出，當(dāng)參考?xì)v元為2009.8時，玉樹震中的東北塊體往東南方向移動;玉樹震中的西南塊體往西南方向移動，而且速度較大;玉樹震中的東南塊體往西南方向移動，但運動速率比西南塊體小。

五、結(jié)束語

本文利用玉樹周圍衛(wèi)星定位連續(xù)運行站的數(shù)據(jù)，均勻選取國內(nèi)外的7個IGS站作為控制，采用常用的GPS數(shù)據(jù)處理軟件GAMIT對玉樹周圍的中國地殼運動觀測網(wǎng)絡(luò)站的數(shù)據(jù)進行處理。首先計算出玉樹地震前后幾天的單日坐標(biāo)解B、L、H坐標(biāo)值，并與2009年相同時段長度的解算結(jié)果進行比較。結(jié)果表明，玉樹地震當(dāng)天高程方向的平均偏差達到16 mm;最大偏差值達到30 mm;緯度方向的最大偏差值達到8 mm;而對經(jīng)度方向的影響較小。然后通過對GPS速度場的計算分析，得出了玉樹周圍各塊體的運動情況。由于僅采用7個網(wǎng)絡(luò)工程站9 d的數(shù)據(jù)，且震中的西北方向沒有網(wǎng)絡(luò)工程站，若采用更長時間序列的觀測數(shù)據(jù)，有望得到玉樹地震對周圍衛(wèi)星定位連續(xù)運行站影響的更全面、準(zhǔn)確的分析。

［1］ 鄧明莉，孫和平，徐建橋，等.GPS數(shù)據(jù)約束的昆侖Ms8.1地震斷層的分段模型［J］.大地測量與地球動力學(xué)，2008，28(4):31-37.

［2］ 沈強，喬學(xué)軍，王琪，等.中國玉樹Mw6.9地震InSAR地表形態(tài)特征分析［J］.大地測量與地球動力學(xué)，2010，30(3):5-9.

［3］ HERRING T A，KING R W，MCCLVSKY S C.GAIMT Reference Manual Release 10.3［R］.Washington，DC: NASA，2006.

［4］ 張四新，張希，王雙緒，等.汶川8.0級地震前后地殼垂直形變分析［J］.大地測量與地球動力學(xué)，2008，28(6):43-46.

［5］ 李建虎，張勝凱，鄂棟臣，等.南極中山站數(shù)據(jù)處理中IGS框架站的選擇［J］.大地測量與地球動力學(xué)，2010，30(1):61-65.

［6］ 張小紅，劉經(jīng)南，F(xiàn)ORSBERG R.基于精密單點定位技術(shù)的航空測量應(yīng)用實踐［J］.武漢大學(xué)學(xué)報:信息科學(xué)版，2006，31(1):19-22.

［7］ 葉世榕.GPS非差相位精密單點定位理論與實現(xiàn)［D］.武漢:武漢大學(xué)，2002.

［8］ WITCHAYANGKOON B.Elements of GPS Precise Point Positioning［D］.Orono:the University of Maine，2000.

［9］ 楊元喜.自適應(yīng)動態(tài)導(dǎo)航定位［M］.北京:測繪出版社，2006:54-55.

［10］ 歐陽永忠，陸秀平，孫紀(jì)章，等.GPS測高技術(shù)在無驗潮水深測量中的應(yīng)用［J］.海洋測繪，2005，25(1): 6-9.

The Analysis of Yushu Earthquake’s Effect on Its Surrounding Satellite Continuous Stations

LI Kaifeng，OUYANG Yongzhong，LU Xiuping，REN Laiping

0494-0911(2011)12-0005-04

P228.4

B

2010-11-16

國家海洋局海底科學(xué)重點實驗室開放基金資助項目(KLSG1002);國家863計劃資助項目(2006AA06A202)

李凱鋒(1981—)，男，江西修水人，碩士，主要從事GNSS數(shù)據(jù)處理技術(shù)及其應(yīng)用研究的工作。