王伶俐 邵德盛 洪 敏

（云南省地震局，昆明 650041）

云南境內(nèi)陸態(tài)網(wǎng)絡(luò)GNSS觀測(cè)資料數(shù)據(jù)處理與初步結(jié)果1

王伶俐 邵德盛 洪 敏

（云南省地震局，昆明 650041）

利用GAMIT/GLOBK軟件對(duì)云南境內(nèi)陸態(tài)網(wǎng)絡(luò)GNSS觀測(cè)資料，包括27個(gè)GNSS基準(zhǔn)站自2011年至今的連續(xù)觀測(cè)資料，以及170多個(gè)GNSS區(qū)域站2009、2011、2013年3期區(qū)域復(fù)測(cè)資料進(jìn)行了處理，基線處理的重復(fù)性在水平方向優(yōu)于1.79mm+1.1×10-9×S，在垂直方向優(yōu)于5.99mm+1.46×10-9×S（S為基線長(zhǎng)度，以mm為單位），計(jì)算得到了基線及點(diǎn)位時(shí)間序列和歐亞框架下的測(cè)站運(yùn)動(dòng)速度場(chǎng)。初步分析了坐標(biāo)時(shí)間序列的變化規(guī)律，并利用剛體運(yùn)動(dòng)模型反演了云南地區(qū)各活動(dòng)塊體的現(xiàn)今地殼運(yùn)動(dòng)，初步結(jié)果表明云南地塊運(yùn)動(dòng)方向逐漸由SSE向SSW變化，運(yùn)動(dòng)幅度由西向東、由北向南逐漸減弱，菱形塊體外各個(gè)塊體運(yùn)動(dòng)速度大幅衰減。該初步結(jié)果可為后期GNSS數(shù)據(jù)的進(jìn)一步運(yùn)用提供參考。

陸態(tài)網(wǎng)絡(luò) 云南 GAMIT/GLOBK 時(shí)間序列 速度場(chǎng)

王伶俐，邵德盛，洪敏，2015.云南境內(nèi)陸態(tài)網(wǎng)絡(luò)GNSS觀測(cè)資料數(shù)據(jù)處理與初步結(jié)果.震災(zāi)防御技術(shù)，10（1）：141—150. doi：10.11899/zzfy20150115

引言

云南處于歐亞、太平洋和印度三大地質(zhì)板塊碰撞匯聚地帶，青藏高原東南邊緣，位居中國(guó)-東南亞次級(jí)板塊中西部，跨揚(yáng)子、華南、印支及滇緬泰4個(gè)亞板塊（申重陽(yáng)等，2002）。其地殼變動(dòng)劇烈，新構(gòu)造活動(dòng)強(qiáng)烈，活動(dòng)斷裂發(fā)育，是中國(guó)最活動(dòng)的地震構(gòu)造區(qū)之一，也是重點(diǎn)監(jiān)視的強(qiáng)震活動(dòng)區(qū)。“十五”期間，“中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)”在全國(guó)范圍內(nèi)建成了25個(gè)連續(xù)跟蹤站點(diǎn)，1000個(gè)區(qū)域站點(diǎn)，地震分析預(yù)報(bào)人員將處理后的時(shí)間序列及速度場(chǎng)成果用于地震年度會(huì)商，在中長(zhǎng)期地震預(yù)報(bào)方面取得了一些有價(jià)值的研究成果（顧國(guó)華等，2000；孟國(guó)杰等，2001；牛之俊等，2002）。但當(dāng)時(shí)由于云南境內(nèi)站點(diǎn)少，點(diǎn)位分散，尚不能滿足地震預(yù)測(cè)預(yù)報(bào)的實(shí)際需求?！笆晃濉逼陂g，由中國(guó)地震局牽頭的國(guó)家重大科學(xué)工程項(xiàng)目“中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)工程”（簡(jiǎn)稱陸態(tài)網(wǎng)絡(luò)）啟動(dòng)，項(xiàng)目在全國(guó)范圍內(nèi)展開了高密度、大范圍的GNSS觀測(cè)，在云南范圍內(nèi)共建成了27個(gè)連續(xù)跟蹤站（GNSS基準(zhǔn)站），170多個(gè)流動(dòng)觀測(cè)站（GNSS區(qū)域站），并陸續(xù)開展了監(jiān)測(cè)工作。GNSS基準(zhǔn)站和GNSS區(qū)域站的建成，在該區(qū)域形成了點(diǎn)位密集且分布均勻的水平形變監(jiān)測(cè)網(wǎng)，獲取了大量的觀測(cè)資料。本文利用GAMIT/GLOBK 軟件對(duì)云南境內(nèi)陸態(tài)網(wǎng)絡(luò)GNSS觀測(cè)資料，包括27個(gè)GNSS基準(zhǔn)站自2011年至今的連續(xù)觀測(cè)資料，以及170多個(gè)GNSS區(qū)域站2009、2011、2013年3期區(qū)域復(fù)測(cè)資料進(jìn)行了處理。初步分析了坐標(biāo)時(shí)間序列的變化規(guī)律，并利用剛體運(yùn)動(dòng)模型反演了云南地區(qū)各活動(dòng)塊體的現(xiàn)今地殼運(yùn)動(dòng)，以期為后期GNSS數(shù)據(jù)的進(jìn)一步運(yùn)用提供參考，為云南省的強(qiáng)震預(yù)測(cè)預(yù)報(bào)提供基礎(chǔ)資料。

1 云南省陸態(tài)網(wǎng)絡(luò)GNSS網(wǎng)點(diǎn)分布與觀測(cè)

陸態(tài)網(wǎng)絡(luò)項(xiàng)目在云南省范圍內(nèi)共建成27個(gè)GNSS基準(zhǔn)站和170多個(gè)GNSS區(qū)域站，再加上其他行業(yè)專項(xiàng)建成的1個(gè)基準(zhǔn)站，目前云南省內(nèi)共有28個(gè)GNSS基準(zhǔn)站，幾乎覆蓋了云南省的大部份地區(qū)，且跨越了瀾滄江、紅河和小江等主要活動(dòng)斷裂帶（圖1）。其中，27個(gè)基準(zhǔn)站基本上都位于穩(wěn)定的基巖上，基準(zhǔn)站采用Trimble NETR8的GPS 接收機(jī)及扼流圈天線，主要進(jìn)行GPS觀測(cè)，也可擴(kuò)充GLONASS、GALILEO和北斗導(dǎo)航系統(tǒng)觀測(cè)。網(wǎng)絡(luò)中所有測(cè)站均為強(qiáng)制歸心的GPS觀測(cè)墩，避免了對(duì)中誤差的影響。自2010年底開始運(yùn)行至今，28個(gè)GNSS站點(diǎn)已經(jīng)積累了4年多的連續(xù)觀測(cè)數(shù)據(jù)；GNSS區(qū)域站點(diǎn)自2009年開始，每2年觀測(cè)1期，到目前為止共有3期的復(fù)測(cè)資料，資料采集的時(shí)間分別為2009年3—10月、2011年4—10月、2013年4—10月，每期每個(gè)測(cè)點(diǎn)的觀測(cè)時(shí)間長(zhǎng)度均不少于90h。

圖1 云南省GNSS基準(zhǔn)站點(diǎn)位分布Fig. 1 Location of the GNSS continuous stations in Yunnan province

2 數(shù)據(jù)解算

2.1 數(shù)據(jù)預(yù)處理

本文采用MIT（Massachusetts Institute of Technology）和SIO（Scripps Institution of Oceanography）研制的GAMIT/GLOBK軟件對(duì)基線進(jìn)行處理，版本為GAMIT10.5。選取HYDE、BJFS、IRKT、PIMO、GUAO、IISC、USUD、TSKB、WUHN、TWTF、LHAZ、URUM等十幾個(gè)中國(guó)及周邊的IGS跟蹤站與省內(nèi)的28個(gè)GNSS基準(zhǔn)站聯(lián)合進(jìn)行數(shù)據(jù)預(yù)處理。其中，IGS核心站強(qiáng)約束（水平方向0.005—0.020m、垂直方向0.020—0.050m）于IERS（International Earth Rotation Service）公布的ITRF2008值；基準(zhǔn)點(diǎn)松弛約束（水平方向1m、垂直方向1m）于其先驗(yàn)值。主要的模型選擇和參數(shù)設(shè)置如下：①數(shù)據(jù)采樣間隔為30s；②根據(jù)軌道松弛模式對(duì)IGS精密星歷給予10-9量級(jí)的約束；③極移和UT1（經(jīng)極移修正的世界時(shí)）采用USNO的Bull_B或Bull_A值一并給予強(qiáng)約束；④衛(wèi)星截止高度角為10°，觀測(cè)值隨高度角定權(quán)；⑤根據(jù)標(biāo)準(zhǔn)大氣模型，用Saastamoinen模型修正對(duì)流層折射；⑥根據(jù)分段線性模型估計(jì)對(duì)流層天頂延遲和水平大氣梯度，每2個(gè)小時(shí)估計(jì)1個(gè)天頂延遲參數(shù)；⑦分別根據(jù)IERS03標(biāo)準(zhǔn)和FES2004標(biāo)準(zhǔn)建立重力場(chǎng)、固體潮的地球物理模型；⑧根據(jù)GMF模型修正大氣層的映射函數(shù)；⑨根據(jù)NGS提供的絕對(duì)模型修正衛(wèi)星和接收機(jī)天線相位中心。

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

為盡快獲得基線時(shí)序資料以滿足云南省地震局地震監(jiān)測(cè)預(yù)報(bào)時(shí)效性的需要，筆者在數(shù)據(jù)處理中對(duì)基準(zhǔn)站數(shù)據(jù)采用的是快速解方式，即在GAMIT數(shù)據(jù)預(yù)處理中用的是快速精密星歷。同時(shí)為了消除或減少全球參考框架基準(zhǔn)點(diǎn)的先驗(yàn)坐標(biāo)誤差和觀測(cè)誤差對(duì)單日解的影響，采用的是局域基準(zhǔn)。下文的基線重復(fù)性評(píng)估均基于該時(shí)間序列（見(jiàn)圖2b）。

為了獲得單日解與ITRF2008參考框架的聯(lián)系，利用GLOBK將GAMIT處理得到的單日松弛解和SOPAC（Scripps Orbital and Permanent Array Center）給出的全球IGS站的單日松弛解合并，最終利用GLOBK通過(guò)7個(gè)參數(shù)的相似變換可得到ITRF2008參考框架下的單日解。綜合幾年來(lái)的單日解，可得到ITRF2008參考框架下點(diǎn)位坐標(biāo)的時(shí)間序列（見(jiàn)圖2a）?？紤]到點(diǎn)位坐標(biāo)的精度，在點(diǎn)位移的求取過(guò)程中對(duì)每4天的數(shù)據(jù)進(jìn)行綜合平差，以求出一個(gè)基于ITRF2000框架的點(diǎn)位靜態(tài)解。綜合2年來(lái)各期基準(zhǔn)站的坐標(biāo)，可得到ITRF2000參考框架下點(diǎn)坐標(biāo)時(shí)間序列，而進(jìn)一步的分析均基于該時(shí)間序列。

對(duì)于GNSS區(qū)域站觀測(cè)數(shù)據(jù)，為了更好地研究川滇菱塊南段的構(gòu)造變形，筆者收集了部分川滇交界以及四川省區(qū)域站點(diǎn)聯(lián)測(cè)數(shù)據(jù)與云南省內(nèi)站點(diǎn)的觀測(cè)資料進(jìn)行聯(lián)合解算；經(jīng)過(guò)GAMIT預(yù)處理獲得了一個(gè)時(shí)段同步觀測(cè)站的數(shù)據(jù)（目前一般一個(gè)時(shí)段的長(zhǎng)度為1天），利用GLOBK將獲取的基線文件和全球的IGS站基線文件聯(lián)合解算，將同1期內(nèi)每天的單天解文件結(jié)合成一個(gè)總的解文件，選擇研究區(qū)外圍比較穩(wěn)定的IGS站在ITRF2008的坐標(biāo)、運(yùn)動(dòng)速率作為約束條件，進(jìn)而可得到本期觀測(cè)內(nèi)所選擇ITRF框架下的測(cè)站的平均坐標(biāo)；將多期結(jié)合后的總的解文件進(jìn)行聯(lián)合解算，則可獲得多期時(shí)間段內(nèi)的測(cè)站速度。

為了更好地反映研究區(qū)域內(nèi)部的構(gòu)造變形，結(jié)合實(shí)際情況在ITRF2008框架基礎(chǔ)上，利用歐亞板塊穩(wěn)定部分作為測(cè)站速度計(jì)算的參考框架，可得到歐亞框架下云南地區(qū)的速度場(chǎng)，如圖3所示。具體操作是選取KOSG、ONSA等10個(gè)分布在歐亞板塊上較為穩(wěn)定的IGS站，構(gòu)建穩(wěn)定的歐亞板塊，在經(jīng)過(guò)相對(duì)歐亞板塊的最佳歐拉旋轉(zhuǎn)后，這10個(gè)站的水平速度殘差均小于1.5mm/a。

圖2 GNSS基準(zhǔn)站坐標(biāo)及基線位置的時(shí)序資料Fig. 2 Time series of the GNSS reference stations

圖3 云南地區(qū)相對(duì)歐亞板塊的GPS水平速度場(chǎng)（2009—2013年）Fig. 3 GPS-derived horizontal velocity filed relative to Eurasia plate in Yunnan area（2009—2013）

3 精度分析

為評(píng)估GPS所測(cè)數(shù)據(jù)的精度，筆者給出了精度分析或統(tǒng)計(jì)檢驗(yàn)，并對(duì)結(jié)果進(jìn)行了如下評(píng)價(jià)。

（1）GAMIT軟件采用的是網(wǎng)解（即全組合解），其同步環(huán)閉合差在基線解算時(shí)已經(jīng)進(jìn)行了分配。對(duì)于GAMIT軟件基線解的同步環(huán)檢核，可以把解的nrms值（標(biāo)準(zhǔn)均方根誤差）作為同步環(huán)質(zhì)量好壞的一個(gè)指標(biāo)，一般要求nrms值小于0.5，而且nrms值越小，基線估算精度越高。實(shí)際單時(shí)段數(shù)據(jù)處理的nrms值均小于0.2。

（2）單日解的重復(fù)性是衡量GPS觀測(cè)質(zhì)量和處理結(jié)果的一個(gè)重要指標(biāo)（王敏等，2005）?；€重復(fù)性是相對(duì)幾個(gè)時(shí)段解而言的，是基線內(nèi)符合精度的反映，其值越小，基線的內(nèi)符合精度越高，基線質(zhì)量越好；反之，內(nèi)符合精度越低，質(zhì)量越差?；€重復(fù)性可表示為：

利用基線長(zhǎng)度及各分量分別對(duì)固定誤差a和比例誤差b進(jìn)行線性擬合，檢驗(yàn)其是否在我國(guó)GPS測(cè)量規(guī)范規(guī)定的等級(jí)范圍內(nèi)。為了得到固定誤差和比例誤差，利用基線重復(fù)性的結(jié)果，采用最小二乘法的一元線性回歸，對(duì)所有的單天解基線向量的重復(fù)性進(jìn)行計(jì)算，將重復(fù)性表示為固定誤差和比例誤差（a+b）ppm的組合形式。表1是在站心地平坐標(biāo)系中，南北方向、東西方向、垂直方向及基線方向的重復(fù)性統(tǒng)計(jì)結(jié)果。

表1 基線向量的重復(fù)性統(tǒng)計(jì)結(jié)果Table 1 Statistics of repeatability of baseline vector

利用GLOBK軟件對(duì)觀測(cè)期間各測(cè)點(diǎn)的點(diǎn)位重復(fù)性進(jìn)行計(jì)算，可得到每個(gè)測(cè)點(diǎn)的時(shí)間序列。筆者的計(jì)算結(jié)果顯示，除個(gè)別站點(diǎn)（下關(guān)XIAG、耿馬YNDM）因?yàn)橛^測(cè)環(huán)境影響，精度較差以外，云南境內(nèi)GNSS基準(zhǔn)站測(cè)點(diǎn)的點(diǎn)位重復(fù)性精度為：南北向分量在5mm以內(nèi)；東西向分量在2mm以內(nèi)；垂直向分量在8mm以內(nèi)。如圖4所示。

從以上統(tǒng)計(jì)結(jié)果來(lái)看，GNSS基準(zhǔn)站基線向量重復(fù)性在水平方向上優(yōu)于1.79mm+ 1.1×10-9×S，垂直方向優(yōu)于5.99mm+1.46×10-9×S（S為基線長(zhǎng)度，以mm為單位），可滿足地震監(jiān)測(cè)的需要。相對(duì)基準(zhǔn)站而言，GNSS區(qū)域站觀測(cè)精度稍低一些，制約定位精度的因素主要是：大氣折射誤差、軌道誤差、天線相位中心標(biāo)定以及測(cè)站季節(jié)性非構(gòu)造變動(dòng)等方面的影響（王琪等，2003）。其中，測(cè)站季節(jié)性非構(gòu)造變動(dòng)對(duì)垂直分量的測(cè)定有直接的影響。研究表明，大氣、降水等負(fù)荷的季節(jié)性變化可引起地殼的周期性變動(dòng)，從YNSD站點(diǎn)時(shí)序資料可以看出，云南境內(nèi)GPS測(cè)站與季節(jié)有關(guān)的變動(dòng)在垂直方向上有近15mm的幅度，但從基線統(tǒng)計(jì)結(jié)果看并不明顯，影響尚不突出（見(jiàn)圖2）。

圖4 站點(diǎn)坐標(biāo)重復(fù)性精度統(tǒng)計(jì)圖Fig. 4 Plot of site coordinates repeatability precision

4 初步結(jié)果的應(yīng)用分析

4.1 GNSS基準(zhǔn)站時(shí)序資料的應(yīng)用分析

（1）利用GLOBK進(jìn)行無(wú)約束平差后的基線成果可以繪制基線的時(shí)間序列圖，但由于單日解存在誤差不利于分析，通過(guò)對(duì)單日解進(jìn)行滑動(dòng)均值處理（窗長(zhǎng)為15天，步長(zhǎng)為1天），可得到平滑的基線時(shí)間序列。以此為基礎(chǔ)并結(jié)合基線與斷層的相對(duì)位置，可以分析斷層的運(yùn)動(dòng)特征。以小江斷裂為例，通過(guò)斜跨小江斷裂的2條基線變化時(shí)序圖（圖5）可以看出，昆明-魚洞、昆明-蒙自基線均表現(xiàn)出小江斷裂帶的左旋走滑特征（洪敏等，2012）。

圖5 跨過(guò)小江斷裂的2條基線時(shí)序圖Fig. 5 Time series of the two baselines which cross the Xiaojiang fault

（2）綜合近年來(lái)云南GNSS基準(zhǔn)站的單日坐標(biāo)解，可以得到ITRF2008參考框架下測(cè)站坐標(biāo)的時(shí)間序列，并由此刻畫出站點(diǎn)的運(yùn)動(dòng)軌跡。對(duì)于坐標(biāo)的時(shí)間序列分析，位移解中描述了相對(duì)運(yùn)動(dòng)的參照點(diǎn)即所謂的位移解基準(zhǔn)，可以有各種不同的選擇，因此位移解具有多解性。通過(guò)大量分析試算發(fā)現(xiàn)（王伶俐等，2012），選擇位于華南塊體的文山GNSS連續(xù)跟蹤站作為云南GPS位移數(shù)據(jù)處理的參考基準(zhǔn)是可行的。以文山GNSS點(diǎn)作為參考，可求出其他點(diǎn)的位移，從而可得到2010年初至2012年3月的點(diǎn)位移軌跡圖像（圖6）（洪敏等，2012）。從圖6可以看出，點(diǎn)位運(yùn)動(dòng)方向從川滇菱形塊體最初的南南東向，到川滇菱形塊體南端轉(zhuǎn)為南向，到滇南至滇西南地區(qū)轉(zhuǎn)為了南西向運(yùn)動(dòng)，相對(duì)于文山GNSS點(diǎn)的各個(gè)點(diǎn)位運(yùn)動(dòng)趨勢(shì)呈現(xiàn)繞阿薩姆構(gòu)造結(jié)順時(shí)針旋轉(zhuǎn)的特征。

圖6 2010年至2012年3月相對(duì)于文山GNSS點(diǎn)的位移圖像Fig. 6 Point displacement from January 2010 to March 2012

4.2 GNSS區(qū)域速度場(chǎng)的應(yīng)用分析

根據(jù)實(shí)測(cè)速度場(chǎng)并按照闞榮舉等（1980，1992）的觀點(diǎn)，將云南地區(qū)劃分為4個(gè)活動(dòng)地塊以及7個(gè)次級(jí)構(gòu)造單元，利用剛體運(yùn)動(dòng)模型可反演云南地區(qū)各活動(dòng)塊體的現(xiàn)今地殼運(yùn)動(dòng)，計(jì)算各個(gè)活動(dòng)塊體的歐拉旋轉(zhuǎn)矢量。采用剛性地塊假設(shè)下的最小二乘擬合方法，可計(jì)算得到各個(gè)活動(dòng)塊體相對(duì)于歐亞板塊的運(yùn)動(dòng)及應(yīng)變特征（如圖7）。經(jīng)過(guò)筆者的計(jì)算，滇東塊體運(yùn)動(dòng)速率為6.22mm/a，運(yùn)動(dòng)方向與華南地塊相近，為E37.53°S；主壓應(yīng)變率為-20.67×10-9/a，其方向近似為W38.02°N；主張應(yīng)變率為13.48×10-9/a；這一計(jì)算結(jié)果與謝富仁等（2004）推測(cè)得到的NWW主壓應(yīng)力優(yōu)勢(shì)方位的結(jié)果一致。滇中地區(qū)整體優(yōu)勢(shì)運(yùn)動(dòng)為10.78mm/a，運(yùn)動(dòng)方向141.44°；主壓應(yīng)變率為-32.71×10-9/a，方向128.20°，主張應(yīng)變率為17.78×10-9/a；最大剪切應(yīng)變率為25.25×10-9/a；這一計(jì)算結(jié)果與程萬(wàn)正等（2003）由大量震源機(jī)制解得到的P軸優(yōu)勢(shì)分布，推測(cè)的NNW壓應(yīng)力場(chǎng)方向一致。滇西北地區(qū)優(yōu)勢(shì)運(yùn)動(dòng)為10.16mm/a，運(yùn)動(dòng)方向163.84°；主壓應(yīng)變率為-15.88×10-9/a，方向163.43°，主張應(yīng)變率為22.94×10-9/a；這一計(jì)算結(jié)果與謝富仁等（1993）描述的滇西北中甸、麗江、大理一帶主要受SE、SSE向應(yīng)力場(chǎng)控制的結(jié)論一致。滇西南地區(qū)的構(gòu)造應(yīng)力場(chǎng)比較復(fù)雜，從筆者得到的速度場(chǎng)的結(jié)果可清楚地看出：靠近瀾滄江斷裂帶站點(diǎn)的活動(dòng)性強(qiáng)，隨著站點(diǎn)到斷裂帶距離的增加，活動(dòng)性逐步減弱，優(yōu)勢(shì)運(yùn)動(dòng)方向大多數(shù)從SN向轉(zhuǎn)至SSW向；從塊體活動(dòng)性來(lái)看，該單元的優(yōu)勢(shì)運(yùn)動(dòng)為6.74mm/a，方向S16.53°W，主壓應(yīng)變率為-17.97×10-9/a，其方位角為19.81°，主張壓應(yīng)變率為29.69×10-9/a，最大剪切應(yīng)變率為23.83×10-9/a；這一結(jié)果與闞榮舉等（1980；1977）和王紹晉等（2005）所獲得的NNE-NE主壓應(yīng)力方向一致。

圖7 云南地區(qū)活動(dòng)塊體運(yùn)動(dòng)速率及主應(yīng)變率圖（2009—2013年）Fig. 7 Horizontal velocity and principal strain rate of active blocks in Yunnan region（2009—2013）

5 結(jié)論

筆者利用GAMIT/GLOBK軟件對(duì)陸態(tài)網(wǎng)絡(luò)云南境內(nèi)觀測(cè)數(shù)據(jù)進(jìn)行了處理，引入周邊10多個(gè)穩(wěn)定IGS站點(diǎn)參與聯(lián)合解算，并將計(jì)算結(jié)果歸算到區(qū)域參考框架下。GAMIT解算得到的“準(zhǔn)觀測(cè)值”均滿足所須精度要求，所有解時(shí)段的標(biāo)準(zhǔn)化均方差均小于0.2；GNSS基準(zhǔn)站水平和垂直方向基線重復(fù)性分別優(yōu)于1.79mm+1.1×10-9×S和5.99mm+1.46×10-9×S（S為基線長(zhǎng)度，以mm為單位）；GNSS區(qū)域網(wǎng)相對(duì)于穩(wěn)定歐亞框架的速度場(chǎng)東西向、南北向的速度估算精度優(yōu)于2.0mm/a。這些指標(biāo)顯示觀測(cè)精度是比較高的，可以滿足地震監(jiān)測(cè)的實(shí)際需要。

通過(guò)GNSS基準(zhǔn)站時(shí)間序列分析，可不斷得到云南省地殼運(yùn)動(dòng)的最新信息；以GPS速度場(chǎng)為約束，建立各個(gè)塊體的運(yùn)動(dòng)-應(yīng)變模型，可得到各個(gè)塊體的運(yùn)動(dòng)及應(yīng)變參數(shù)。初步結(jié)果表明，云南各活動(dòng)塊體由東至西運(yùn)動(dòng)方向，正逐漸向SSE-SN-SSW變化；其表現(xiàn)出以菱形塊體為主順時(shí)針旋轉(zhuǎn)的特征，且菱形塊體外各個(gè)塊體運(yùn)動(dòng)速度明顯衰減。從應(yīng)變率參數(shù)結(jié)果來(lái)看，華南地塊和滇東塊體主要受SE向壓應(yīng)力場(chǎng)控制；滇西北地區(qū)，應(yīng)力場(chǎng)方向?yàn)镾E-SSE；到滇中地塊，應(yīng)力場(chǎng)方向也轉(zhuǎn)為SE-SSE方向；滇西南受NE、NNE向應(yīng)力場(chǎng)控制，結(jié)果與前人的研究一致（王紹晉等，2005；崔效鋒等，2005；2006）。

陸態(tài)網(wǎng)絡(luò)將在未來(lái)繼續(xù)進(jìn)行觀測(cè)，從而完善現(xiàn)有的GNSS觀測(cè)資料，它將提供更為精細(xì)的應(yīng)變場(chǎng)解算結(jié)果，為云南地震預(yù)測(cè)提供更為精細(xì)的基礎(chǔ)資料。

程萬(wàn)正，刁桂苓，呂戈培等，2003. 川滇地塊的震源力學(xué)機(jī)制、運(yùn)動(dòng)速率和活動(dòng)方式. 地震地質(zhì)，25（1）：71—87.

崔效鋒，謝富仁，趙建濤，2005. 中國(guó)及鄰區(qū)震源機(jī)制解的分區(qū)特征. 地震地質(zhì)，27（2）：298—307.

崔效鋒，謝富仁，張紅艷，2006. 川滇地區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)分區(qū)及動(dòng)力學(xué)意義. 地震學(xué)報(bào)，28（5）：451—461.

顧國(guó)華，牛紅葉，孟國(guó)杰，許永江，孫漢榮，2000. 中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)數(shù)據(jù)處理與分析. 地震，20（增刊）：203—210.

洪敏，邵德盛，李春光，王伶俐，2012. 云南省GNSS連續(xù)觀測(cè)資料應(yīng)用研究. 地震研究，（4）：535—539.

闞榮舉，1980. 西南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)與板內(nèi)強(qiáng)震活動(dòng). 地震研究，3（3）：45—59.

闞榮舉，韓源，1992. 云南遮放至馬龍地學(xué)斷面（說(shuō)明書）. 北京：地震出版社.

闞榮舉，張四昌，晏鳳桐等，1977. 我國(guó)西南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)與現(xiàn)代構(gòu)造活動(dòng)特征的探討. 地球物理學(xué)報(bào)，20（2）：96—109.

孟國(guó)杰，趙承坤，顧國(guó)華等，2001. 中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)首期觀測(cè)數(shù)據(jù)處理. 地震，21（1）：63—68.

牛之俊，馬宗晉，陳鑫連等，2002. 中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò). 大地測(cè)量與地球動(dòng)力學(xué)，22（3）：1—7.

申重陽(yáng)，王琪，吳云等，2002. 川滇地區(qū)活動(dòng)地塊現(xiàn)今地殼形變特征. 地球物理學(xué)報(bào)，47（5）：805—811.

王琪，牛之俊，石俊成，2003. 中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)基本站觀測(cè)精度研究. 大地測(cè)量與地球動(dòng)力學(xué)，23（3）：1—9.

王伶俐，邵德盛，洪敏，李春光，施發(fā)奇，2012. 云南省GNSS基準(zhǔn)站數(shù)據(jù)處理不同參考基準(zhǔn)研究. 地震研究，（4）：449—454.

王紹晉，龍曉帆，余慶坤，2005. 昆明地區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)分析. 地震研究，28（2）：178—184.

王敏，張祖勝，許明元等，2005. 2000國(guó)家GPS 大地控制網(wǎng)的數(shù)據(jù)處理和精度評(píng)估. 地球物理學(xué)報(bào)，48（4）：817—823.

謝富仁，崔效鋒，趙建濤等，2004. 中國(guó)大陸及鄰區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)分區(qū). 地球物理學(xué)報(bào)，47（4）：654—662.

謝富仁，祝景忠，梁海天等，1993. 中國(guó)西南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)基本特征. 地震學(xué)報(bào)，15（4）：407—417.

楊國(guó)華，韓月萍，王敏，馬慶尊，2003. 中國(guó)大陸幾個(gè)主要地震活動(dòng)區(qū)的水平形變. 大地測(cè)量與地球動(dòng)力學(xué)，23（3）：42—49.

Data Processing and Preliminary Results of GNSS Observation based on CMONOC in Yunnan Province

Wang Lingli, Shao Desheng and Hong Min
（Earthquake Administration of Yunnan Province, Kunming 650041, China）

The GNSS observation data in Yunnan region, including continuous observation data from 27

tations since 2011, and a total of three regional survey data from more than 170 regional stations between 2009 to 2013, were processed by using GAMIT/GLOBK. It is shown that the horizontal and vertical proportion errors of baseline repetition rates can reach up to 1.79mm+1.1×10-9and 5.99mm+1.46×10-9, and time series of reference stations and velocity field relative to the Eurasian reference frame was obtained. The results indicate that the direction of crustal movement to the area changes from SSE to SSW gradually. The velocities of those blocks outside the rhombic block decrease swiftly. These preliminary results provide a reference for further use of the GNSS observation in Yunnan Province.

CMONOC；Yunnan；GAMIT/GLOBK；Timeseries；Velocity filed

云南省地震局青年基金課題 “基于高密度GPS區(qū)域網(wǎng)復(fù)測(cè)資料分析云南地區(qū)近期地殼形變特征（201301）” 和中國(guó)地震局2015年度震情跟蹤青年課題“云南省GNSS時(shí)序資料異常跟蹤與地震預(yù)測(cè)（2015010219）”聯(lián)合資助

2014-06-17

王伶俐，女，生于1981年。工程師，碩士。主要從事GPS數(shù)據(jù)處理及分析、地震活動(dòng)性分析及預(yù)測(cè)預(yù)報(bào)工作。E-mail：lingli365@163.com