李經(jīng)緯 陳長(zhǎng)云 占 偉 武艷強(qiáng)

1)中國(guó)地震局第一監(jiān)測(cè)中心， 天津 300180 2)中國(guó)科學(xué)院精密測(cè)量科學(xué)與技術(shù)創(chuàng)新研究院， 武漢 430071

0 引言

據(jù)中國(guó)地震臺(tái)網(wǎng)正式測(cè)定， 2021年5月22日2時(shí)4分， 在青海省果洛藏族自治州瑪多縣(34.59°N， 98.34°E)發(fā)生7.4級(jí)地震， 震源深度17km， 此后發(fā)生了多次4～5級(jí)余震(圖1b)。經(jīng)專(zhuān)家學(xué)者現(xiàn)場(chǎng)考察及地震序列的精定位初步確定， 此次地震發(fā)生昆侖山口-江錯(cuò)斷裂(李智敏等， 2021； 潘家偉等， 2021； 王未來(lái)等， 2021)(圖 1)。昆侖山口-江錯(cuò)斷裂位于巴顏喀拉塊體北邊界——東昆侖斷裂帶以南， 向E與東昆侖斷裂銜接， 向W延伸至甘德附近， 全長(zhǎng)>500km， 位于瑪多-甘德斷裂與甘德南緣斷裂之間(鄧起東等， 2007； 徐錫偉等， 2016； 李智敏等， 2021)。該斷裂所在地區(qū)發(fā)育一系列大型活動(dòng)斷層， 包括西藏大溝-昌馬河斷裂、 瑪多-甘德斷裂、 甘德南緣斷裂及達(dá)日斷裂等。第四紀(jì)以來(lái)該地區(qū)地震活動(dòng)性顯著， 近百年來(lái)發(fā)生了多次大地震。例如， 此次地震震中NW的東昆侖斷裂在1937年和1963年發(fā)生了2次7級(jí)以上地震、 SE的達(dá)日斷裂中段在1947年發(fā)生了M7.7地震(中國(guó)地震局震害防御司， 1999)。該地區(qū)的GNSS在ITRF參考框架下呈現(xiàn)N偏E的運(yùn)動(dòng)趨勢(shì)， EW向運(yùn)動(dòng)分量由西向東逐漸增大， 由南向北逐漸減小。研究表明， 北部東昆侖斷裂帶的走滑速率自西向東由13～14mm/a逐步減弱至3mm/a， 其中托索湖段約為5mm/a(高翔等， 2013； 梁明劍等， 2020)。Zhu等(2021)通過(guò)2005—2010年的GNSS觀測(cè)數(shù)據(jù)給出昆侖山口-江錯(cuò)的滑動(dòng)速率約為(1.2±0.8)mm/a。雖然該斷層的運(yùn)動(dòng)速率小于附近的東昆侖托索湖段運(yùn)動(dòng)速率， 但這次地震的發(fā)生證明該地區(qū)的構(gòu)造形變發(fā)生在整個(gè)區(qū)域而不是少數(shù)幾條主要斷層上(張培震等， 2003)。

圖 1 瑪多7.4級(jí)地震的構(gòu)造背景及余震分布Fig. 1 The tectonic setting of the Maduo MS7.4 earthquake and the distribution of aftershocks.a 青海瑪多7.4級(jí)地震震中區(qū)歷史地震、 活動(dòng)斷裂及GNSS速度場(chǎng)； b 青?，敹?.4級(jí)地震的主震及余震分布。斷層數(shù)據(jù)來(lái)自于鄧起東等(2007)和徐錫偉等(2016)， 同震和余震數(shù)據(jù)來(lái)自于USGS(1)https: ∥earthquake.usgs.gov/earthquakes。， GNSS速度場(chǎng)結(jié)果由中國(guó)地震局第一監(jiān)測(cè)中心計(jì)算(2) ftp: ∥60.30.77.19。

同震近、 遠(yuǎn)場(chǎng)位移場(chǎng)是認(rèn)識(shí)發(fā)震斷裂運(yùn)動(dòng)學(xué)、 動(dòng)力學(xué)特征必不可少的資料(Banerjeeetal.， 2005； Ganetal.， 2007； Wuetal.， 2016； Wangetal.， 2020)， 快速產(chǎn)出該資料對(duì)于震后快速判定未來(lái)地震趨勢(shì)、 研究區(qū)域地殼變形特征及地震動(dòng)力學(xué)特征具有極其重要的意義?？焖佾@取GNSS測(cè)站的坐標(biāo)常采用動(dòng)態(tài)數(shù)據(jù)解算的方法， 可在發(fā)震后利用15min的高頻GNSS觀測(cè)數(shù)據(jù)獲取， 常用的軟件包括各類(lèi)精密單點(diǎn)定位(PPP)軟件及GAMIT/GLOBK的TRACK動(dòng)態(tài)差分處理模塊(3)Herring T， King R， Mcclusky S， 等， 2018， GAMIT/GLOBK reference manuals， Release 10.7， MIT Technical Reports。(Psimoulisetal.， 2016)。由于受到星歷、 鐘差等產(chǎn)品精度的限制， PPP的精度常為cm級(jí)(Guoetal.， 2015； Suetal.， 2020)。TRACK模塊采用雙差卡爾曼濾波方法可消除部分誤差， 但其主要適用于中短距離的基線(xiàn)， 定位精度也在cm級(jí)(殷海濤等， 2012； Tiryakiogluetal.， 2017)。除了采用GNSS高頻數(shù)據(jù)動(dòng)態(tài)定位外， 可利用4～8h觀測(cè)數(shù)據(jù)獲取mm級(jí)精度的定位結(jié)果， 這在測(cè)繪工程的高等級(jí)控制網(wǎng)觀測(cè)中已有眾多案例(周命端等， 2015)。因此， 可利用地震前、 后數(shù)h的靜態(tài)觀測(cè)數(shù)據(jù)采用雙差的靜態(tài)解算模式快速獲取同震位移結(jié)果。

本文利用震中所在的巴顏喀拉塊體內(nèi)部及周邊18個(gè)連續(xù)GNSS站點(diǎn)的數(shù)據(jù)(其中7個(gè)測(cè)站記錄了1Hz高頻數(shù)據(jù))， 采用震前、 后15min1Hz高頻數(shù)據(jù)動(dòng)態(tài)解算， 震后5h靜態(tài)解算和多天靜態(tài)解算3種方法獲取了同震位移并進(jìn)行對(duì)比分析， 探索和測(cè)試了采用震前、 震后數(shù)h的觀測(cè)資料靜態(tài)解算獲取形變資料的可行性， 該測(cè)試成果快速揭示了青?，敹?.4級(jí)地震產(chǎn)生的近、 遠(yuǎn)場(chǎng)同震位移， 為災(zāi)情評(píng)估和地震趨勢(shì)的快速判定提供了基礎(chǔ)資料， 有助于加深對(duì)青?，敹?.4級(jí)地震發(fā)生機(jī)理的認(rèn)識(shí)。

1 GNSS數(shù)據(jù)及處理方法

地震發(fā)生后， 中國(guó)地震局臺(tái)網(wǎng)中心及時(shí)發(fā)布了周邊地區(qū)的 “陸態(tài)網(wǎng)絡(luò)”GNSS連續(xù)站的原始觀測(cè)數(shù)據(jù)。本文挑選了震中附近的18個(gè)GNSS連續(xù)觀測(cè)站的數(shù)據(jù)。觀測(cè)數(shù)據(jù)可分為2類(lèi)， 一類(lèi)是采樣間隔為30s的觀測(cè)數(shù)據(jù)， 包含同震當(dāng)天、 震前3d和震后1d的數(shù)據(jù)； 另一類(lèi)是其中7個(gè)測(cè)站記錄的地震當(dāng)天1Hz的高頻觀測(cè)數(shù)據(jù)。針對(duì)這2類(lèi)數(shù)據(jù)， 分別采用靜態(tài)解算和動(dòng)態(tài)解算方法進(jìn)行計(jì)算。在進(jìn)行靜態(tài)解算時(shí)， 采用GAMIT/GLOBK網(wǎng)解模式。 為滿(mǎn)足網(wǎng)解需求，本文選取了18個(gè)中國(guó)大陸周邊的IGS站作為框架參考站。首先利用GAMIT/GLOBK軟件獲取GNSS基線(xiàn)松弛解， 然后以IGS站點(diǎn)作為框架點(diǎn)采用GLOBK進(jìn)行約束平差獲得ITRF14參考框架下的測(cè)站坐標(biāo)①。這種雙差模式的優(yōu)勢(shì)在于能夠在雙差組合方程中消除鐘差、 電離層延遲等帶來(lái)的誤差影響。在動(dòng)態(tài)解算中， 本文采用GAMIT/GLOBK的TRACK模塊進(jìn)行計(jì)算， 其可利用震前、 震后約10min的高頻數(shù)據(jù)獲取地震發(fā)生瞬間的動(dòng)態(tài)變化過(guò)程。在解算中首先需要找一個(gè)基準(zhǔn)點(diǎn)， 然后進(jìn)行載波相位動(dòng)態(tài)差分解算， 這種方法獲取結(jié)果的速度較快， 但基線(xiàn)>100km時(shí)的解算精度為cm級(jí)(范士杰等， 2013)。

無(wú)論是進(jìn)行靜態(tài)解算還是動(dòng)態(tài)解算， 都是為了快速而準(zhǔn)確地獲取地震前后的坐標(biāo)， 從而獲取同震位移。在獲取震前、 震后坐標(biāo)時(shí)， 根據(jù)不同的解算模式及所用數(shù)據(jù)的時(shí)長(zhǎng)， 將其劃分為3種方法： 1)利用1Hz高頻觀測(cè)數(shù)據(jù)進(jìn)行動(dòng)態(tài)差分定位獲取逐歷元解， 然后分別在震前、 震后取平均作差獲得最終同震位移， 這種方法可以在獲取高頻數(shù)據(jù)后快速產(chǎn)出結(jié)果； 2)僅采用震后數(shù)h的觀測(cè)值進(jìn)行靜態(tài)解算， 解算前需要將單天數(shù)據(jù)中的震前、 震后數(shù)據(jù)分割為獨(dú)立的文件進(jìn)行靜態(tài)解算， 針對(duì)此次青?，敹?.4級(jí)地震， 本文將同震當(dāng)天的數(shù)據(jù)去除UTC18:00—18:30時(shí)刻的數(shù)據(jù)后分割為2個(gè)觀測(cè)文件， 同時(shí)為了驗(yàn)證5h結(jié)果的可靠性， 對(duì)震前18h的文件也進(jìn)行了分割處理； 3)利用多天的數(shù)據(jù)獲取精確解， 本文選取了震前4d和震后1d獲取的穩(wěn)定結(jié)果， 由于該地區(qū)在ITRF參考框架下， 變化最大的方向量值≤50mm/a， 即單天變化約為0.1mm， 因此認(rèn)為在震前3～4d測(cè)站是固定的， 不存在較大的線(xiàn)性速率問(wèn)題， 在獲取震前坐標(biāo)時(shí)直接采用4d的平均值。

2 結(jié)果與討論

2.1 1Hz高頻數(shù)據(jù)的動(dòng)態(tài)差分方法結(jié)果

依據(jù)TRACK模塊的解算特點(diǎn)， 其獲取的結(jié)果為2點(diǎn)之間基線(xiàn)長(zhǎng)度的動(dòng)態(tài)變化值， 因此對(duì)應(yīng)地選取UTC 17:45—18:15時(shí)段的數(shù)據(jù)進(jìn)行動(dòng)態(tài)差分解算， 衛(wèi)星星歷為IGS提供的快速星歷(igr)。圖2a 展示了GSMA-QHMQ測(cè)站的結(jié)果， 基線(xiàn)向W和向S分別發(fā)生了近1cm和近2cm的位移。圖2b 展示了XZCD-QHMD測(cè)站的結(jié)果， 從結(jié)果中可以明顯地看出地震發(fā)生過(guò)程中在近斷層處的位移過(guò)程， 基線(xiàn)向E發(fā)生了近22cm的位移， 向N發(fā)生了近6cm的位移。圖2c 展示了GSMA-QHBM測(cè)站結(jié)果， 基線(xiàn)顯示向E發(fā)生了近1.4cm位移。利用此結(jié)果結(jié)合地質(zhì)構(gòu)造背景可以基本判斷此次地震發(fā)生于左旋走滑型斷層上。

圖 2 利用高頻GNSS數(shù)據(jù)獲取的同震位移及分布(參考站為GSMA)Fig. 2 Coseismic displacement and distribution gained from high-frequency GNSS data(reference station GSMA).a 利用高頻GNSS數(shù)據(jù)通過(guò)動(dòng)態(tài)差分獲取的同震位移(相對(duì)于GSMA站)分布； b GSMA-QHMD基線(xiàn)的長(zhǎng)度變化； cGSMA-QHBM基線(xiàn)的長(zhǎng)度變化。圖b、 c中的灰色虛線(xiàn)表示瑪多地震發(fā)生時(shí)刻所在歷元， 藍(lán)色實(shí)線(xiàn)表示基線(xiàn)坐標(biāo)差值序列， 紅色實(shí)線(xiàn)表示震前、 震后各自的中位值線(xiàn)

從圖 2 中我們發(fā)現(xiàn)， 由于TRACK模塊采用卡爾曼濾波方法， 解算中歷元間值存在關(guān)聯(lián)， 在地震瞬時(shí)歷元附近的序列會(huì)抖動(dòng)， 作為離震中較遠(yuǎn)的基準(zhǔn)站GSMA受到地震影響的時(shí)間晚于其他站， 因此在基線(xiàn)時(shí)間序列中存在2次抖動(dòng)的特征。由于GSMA站與QHBM站距離較近， 2次抖動(dòng)存在重疊， 同時(shí)基線(xiàn)N向位移量值較小， 在圖2c 中顯示不明顯。在解算精度方面， 經(jīng)統(tǒng)計(jì)基線(xiàn)解算收斂后的整體解算精度為1～2cm， 除上述幾個(gè)測(cè)站外其他基線(xiàn)的變化量值較小， 離散度已經(jīng)超過(guò)了變化量值， 較難采用平均值或中值的辦法來(lái)獲取準(zhǔn)確的位移量。除此之外， 從地震發(fā)生時(shí)刻到地震波傳播到測(cè)站位置時(shí)存在時(shí)間差， 由于數(shù)據(jù)離散程度較大， 當(dāng)位移量較小時(shí)較難從數(shù)據(jù)的角度確定準(zhǔn)確的到達(dá)時(shí)刻。因此， TRACK模塊受制于精度的影響， 對(duì)于本次解算中位移量較大的QHMD、 QHMQ站可以采用這種方法合理地估算出位移量， 而對(duì)于其他測(cè)站則較難獲取合理的位移量值， 但從位移的趨勢(shì)與分布可以準(zhǔn)確地呈現(xiàn)斷層兩側(cè)的運(yùn)動(dòng)現(xiàn)象， 在快速揭示破裂性質(zhì)方面仍有著重要作用。

2.2 震后5h數(shù)據(jù)的靜態(tài)解算方法結(jié)果

此次地震發(fā)生在5月21日UTC18:04:13時(shí)， 在震后第2天(北京時(shí)間5月22日)獲取了相關(guān)數(shù)據(jù)， 地震當(dāng)天的數(shù)據(jù)中包括震后5h的觀測(cè)值， 滿(mǎn)足GNSS靜態(tài)雙差處理的時(shí)長(zhǎng)需求。衛(wèi)星星歷采用IGS中心發(fā)布的快速星歷(igr)。采用GFZRNX軟件(4)Nischan T， 2016， GFZRNX-RINEX GNSS Data Conversion and Manipulation Toolbox， GFZData Services， http: ∥dx.doi.org/10.5880/GFZ.1.1.2016.002。對(duì)當(dāng)天的數(shù)據(jù)進(jìn)行分時(shí)段提取， 提取時(shí)段分為第1時(shí)段(UTC00:00—18:00)、 第2時(shí)段(UTC00:00—09:00)、 第3時(shí)段(UTC09:00—18:00)及第4時(shí)段(UTC18:30—24:00)。同時(shí)假定UTC00:00—18:00時(shí)段的結(jié)果為真， 以此為基準(zhǔn)經(jīng)處理后獲得的每個(gè)測(cè)站的震前水平方向坐標(biāo)差值見(jiàn)表1， 第4段與第1時(shí)段的差值表示GNSS的同震位移結(jié)果， 如圖 3 所示。

表1 利用不同時(shí)間段的數(shù)據(jù)獲取的坐標(biāo)值差異(以UTC00:00—18:00為基準(zhǔn))(單位： mm)Table1 The difference of coordinate values obtained from data of different time periods (UTC00:00—18:00 as the benchmark)(unit: mm)

圖 3 青?，敹?.4級(jí)地震GNSS震后5h解和多天單日解獲取的同震位移結(jié)果Fig. 3 Coseismic displacement obtained respectively from the GNSS station’s five-hour and multi-day solutions of the MS7.4 Maduo earthquake in Qinghai Province.箭頭表示測(cè)站的同震位移， 紅色和藍(lán)色箭頭表示震后5h靜態(tài)解結(jié)果， 金色和綠色箭頭表示多天單日解結(jié)果

表1 中， 第1時(shí)段和第2時(shí)段結(jié)果與第4時(shí)段坐標(biāo)結(jié)果差異的中位值為0.3mm。除了青海的QHYS站和西藏的XZCD站及震中距約為40km的QHMD站(青?，敹?在震前的結(jié)果差異>2mm外， 其他測(cè)站的差異均<2mm。經(jīng)統(tǒng)計(jì)， 第2、 3時(shí)段水平方向解算精度的中位值為5.4mm， 第4時(shí)段為6.7mm， 第1時(shí)段為2.2mm。結(jié)果表明， 觀測(cè)時(shí)長(zhǎng)縮短會(huì)降低解算中的內(nèi)符合精度， 雖然表1 中展現(xiàn)的結(jié)果差異大部分<2mm， 但第4時(shí)段解算的內(nèi)符合精度約為6.7mm， 同時(shí)由于在震后5h受到諸多余震的影響， 其結(jié)果也可能存在偏差。我們認(rèn)為將位移量值小于第3時(shí)段內(nèi)符合精度值的結(jié)果作為參考值是可行的， 但作為最終結(jié)果并不可靠。圖 3中的藍(lán)色箭頭展示了5h靜態(tài)解算獲得的位移場(chǎng)， 結(jié)果表現(xiàn)出明顯的左旋走滑特性， 這與該地區(qū)的地質(zhì)構(gòu)造背景運(yùn)動(dòng)性質(zhì)一致。在距震中200～300km的QHDL、 QHMQ站發(fā)生了約10mm位移， 在QHBM、 QHYS發(fā)生了約5mm的位移， 超過(guò)該范圍的測(cè)站的位移<5mm。綜上所述， 分析認(rèn)為采用震后數(shù)h獲取的雙差靜態(tài)結(jié)果作為快速產(chǎn)出的同震位移場(chǎng)資料在近場(chǎng)基本可靠， 在遠(yuǎn)場(chǎng)位移量值較小的測(cè)站位移可作為初步的參考值。

2.3 多天數(shù)據(jù)的靜態(tài)解算方法結(jié)果

為了獲取更為準(zhǔn)確的同震形變結(jié)果， 本文對(duì)震前4d及震后1d的數(shù)據(jù)進(jìn)行了分析處理， 在處理前去除了地震當(dāng)天數(shù)據(jù)文件中震后時(shí)段(UTC 18:30—24:00)的數(shù)據(jù)。衛(wèi)星星歷仍然采用IGS提供的快速星歷(igr)， 同時(shí)添加高階電離層改正。經(jīng)解算統(tǒng)計(jì)獲取的多天單日解基線(xiàn)NRMS<0.25， 約束平差后的測(cè)站坐標(biāo)解算精度為1～2mm， 獲取的GNSS同震位移結(jié)果如圖 3 中的綠色箭頭所示。

圖 3 結(jié)果顯示， 多天靜態(tài)解算與震后5h靜態(tài)解算獲取的同震位移結(jié)果基本保持一致， 從表2 發(fā)現(xiàn)兩者的數(shù)值相差1～4mm。從數(shù)值和地理分布分析， 近場(chǎng)測(cè)站的EW向和向SN向位移都呈現(xiàn)繼續(xù)在原有方向上增加的趨勢(shì)， 且離震中越近這種增量越大， 如QHMD站增加了約4mm位移、 QHMQ增加了約1.5mm位移。造成這種現(xiàn)象主要是由于在斷層面瞬間發(fā)生破裂后斷層仍將保持快速滑動(dòng)而使得形變累積增大， 這種現(xiàn)象也是早期震后形變的重要特征(Liu-Zengetal.， 2020； Lietal.， 2021)。當(dāng)我們使用震后1d的單日解作為震后測(cè)站坐標(biāo)時(shí)， 其代表震后1d的平均坐標(biāo)。因此， 如果考慮震后的持續(xù)變形過(guò)程， 其坐標(biāo)值必將大于震后5h的靜態(tài)解結(jié)果， 因此存在這種差值是合理的， 在分析dm級(jí)的同震位移時(shí)可以忽略。但如果作為震后形變分析， 表明震后十幾h內(nèi)近場(chǎng)發(fā)生了顯著變形， 例如QHMD站出現(xiàn)了W向3.7mm、 N向3.3mm的位移， 因此該地區(qū)后續(xù)的震情仍需持續(xù)關(guān)注。

表2 利用震后5h數(shù)據(jù)與多天數(shù)據(jù)通過(guò)靜態(tài)解方法獲取的同震位移量及差值(單位： mm)Table2 Coseismic displacement and difference obtained by postseismic 5-hour and multi-day static solutions， respectively(unit： mm)

除此之外， 較遠(yuǎn)區(qū)域測(cè)站的多天靜態(tài)結(jié)果與5h靜態(tài)解結(jié)果相比量值更小。造成這種現(xiàn)象的原因主要是由于震后5h數(shù)據(jù)受到觀測(cè)時(shí)長(zhǎng)及余震的影響， 其解算結(jié)果的內(nèi)符合精度較低， 多天靜態(tài)解算的精度能夠提高1～2mm， 在位移量值較小的情況下這種改正非常顯著。另外， 青海玉樹(shù)測(cè)站(QHYS)的運(yùn)動(dòng)速率呈現(xiàn)與其他測(cè)站運(yùn)動(dòng)不一致的現(xiàn)象， 該測(cè)站的解算精度較其他測(cè)站低， 具體還需后續(xù)與測(cè)站環(huán)境等因素結(jié)合進(jìn)行討論。

2.4 同震形變特征

根據(jù)本文所獲得的GNSS位移場(chǎng)結(jié)果， 結(jié)合華俊等(2021)的InSAR結(jié)果、 王未來(lái)等(2021)的余震精定位結(jié)果及李智敏等(2021)的野外考察結(jié)果， 綜合分析認(rèn)為此次同震地表形變的影響范圍主要集中在斷裂帶兩側(cè)97°～101°E的范圍內(nèi)。距離震中約40km的QHMD測(cè)站的W向位移約為240mm， N向位移約為90mm， 運(yùn)動(dòng)方向與斷層的走向保持一致。在此方向延伸處的QHGE站(格爾木)存在同向的微小位移。震中北側(cè)的QHDL和DLHA站呈現(xiàn)NW向的位移， 但與QHMD站的運(yùn)動(dòng)方向相比明顯更偏于N向。位于青?，斍腝HMQ測(cè)站的S向運(yùn)動(dòng)約為8mm， W向運(yùn)動(dòng)約為11mm， S向運(yùn)動(dòng)的比例相對(duì)于同區(qū)域QHMD、 QHDL站的位移方向存在明顯差異， 造成這種現(xiàn)象可能是QHMQ站位于破碎帶上。位于斷層另一側(cè)的QHBM站(班瑪)， 其S向位移約為6mm， E向位移約為4mm， 呈現(xiàn)SE向的運(yùn)動(dòng)位移。同時(shí)甘孜、 爐霍等地區(qū)受到較小的波動(dòng)， SCGZ、 SCLH和SCXJ測(cè)站的E向位移不明顯， 呈現(xiàn)S向位移特征， 位移量約為4mm。整體而言， 斷裂北側(cè)呈現(xiàn)NW向位移， 南側(cè)呈現(xiàn)SE向位移， 在遠(yuǎn)離斷層處垂直于斷層方向存在背離斷層面的運(yùn)動(dòng)趨勢(shì)。

3 結(jié)論

本文對(duì)通過(guò)1Hz高頻數(shù)據(jù)動(dòng)態(tài)差分、 震后5h靜態(tài)解算及多天靜態(tài)解算3種方法獲取的青?，敹?.4級(jí)地震地表同震位移開(kāi)展了對(duì)比分析， 得到以下結(jié)論：

(1)TRACK模塊具有快速解算的能力， 其高頻動(dòng)態(tài)差分結(jié)果在近場(chǎng)位移量較大的站點(diǎn)上同震位移顯著， 但在遠(yuǎn)場(chǎng)受測(cè)站距離所限解算精度為cm級(jí)， 較難顯示<1cm的形變量。雖然TRACK的結(jié)果不精確， 但可以作為快速獲取近場(chǎng)測(cè)站位移量及初步判定斷層運(yùn)動(dòng)特征的手段；

(2)震后5h靜態(tài)解算獲得的同震位移與多天靜態(tài)解算結(jié)果相比， 兩者的位移大小、 方向以及影響范圍都保持較強(qiáng)的一致性。多天靜態(tài)解的同震結(jié)果與震后5h靜態(tài)解的結(jié)果相差1～4mm。因此， 考慮震后5h數(shù)據(jù)的時(shí)效性， 認(rèn)為震后5h靜態(tài)解算方法的結(jié)果作為同震位移快速結(jié)果是可信的。

(3)綜合多方面結(jié)果開(kāi)展分析， 可初步判斷此次地震是以左旋走滑為主的地震。目前測(cè)站中離震中最近的約40km的QHMD站的最大位移約為240mm。地震的影響范圍波及了震中瑪多縣周邊區(qū)域， 包括北部西寧、 都蘭、 德令哈地區(qū)及南部斑馬、 甘孜地區(qū)。由于多天靜態(tài)解包含了震后1d的形變， 其形變方向與震后5h形變保持方向一致， 且量值增加， 這種現(xiàn)象與震后余滑較吻合， 因此認(rèn)為近場(chǎng)的震后形變現(xiàn)象顯著， 后期仍需持續(xù)關(guān)注。

致謝中國(guó)地震局臺(tái)網(wǎng)中心提供了GNSS高頻數(shù)據(jù)， 在此表示衷心感謝!