徐克科,伍吉倉

(1.河南理工大學(xué) 測繪與國土信息工程學(xué)院，河南 焦作 454000；2.同濟大學(xué) 測繪與地理信息學(xué)院，上海 200092)

?

川滇地塊現(xiàn)今地殼運動速度場解算與分析

徐克科1,伍吉倉2

(1.河南理工大學(xué) 測繪與國土信息工程學(xué)院，河南 焦作 454000；2.同濟大學(xué) 測繪與地理信息學(xué)院，上海 200092)

利用GAMIT10.4和BERNESE5.2軟件，分別對川滇地區(qū)中國地殼運動觀測網(wǎng)絡(luò)51個GNSS基準(zhǔn)站2010—2014年的觀測數(shù)據(jù)進行處理。兩者解算所得單日解基線和坐標(biāo)時間序列離散度優(yōu)于10 mm，且變化趨勢一致，吻合度較高?；€重復(fù)率水平向優(yōu)于1+4×10-9L；垂向優(yōu)于4+6×10-9L(單位mm，其中L表示基線長度，單位m)。利用PL+VW噪聲模型進行了速度估計，獲取了川滇地區(qū)相對歐亞板塊地殼運動速度場。結(jié)果顯示：除滇西南地區(qū)呈現(xiàn)出由西向東的增大趨勢外，川滇塊體其它區(qū)域地殼運動特征由北往南、由西往東呈逐漸衰減趨勢；運動方向持續(xù)圍繞喜馬拉雅東結(jié)點作順時針旋轉(zhuǎn)，且仍然存在著南北向強烈的擠壓特征。以安寧河-則木河斷裂、小江斷裂和紅河斷裂為界，選用了兩側(cè)能反映震間斷裂構(gòu)造特征的GPS測站，分析了斷裂兩側(cè)的速度差異。并利用平行斷層的速度剖面，擬合反正切函數(shù)變化趨勢，獲取了斷裂帶可能的變形寬度。

GPS；川滇地塊；地殼運動；速度場；斷裂帶

隨著中國地殼運動觀測網(wǎng)絡(luò)的開展，連續(xù)和分期GNSS觀測的持續(xù)，積累了時空分辨率越來越高的地表形變觀測數(shù)據(jù)。如何從一手的觀測資料中獲取高信噪比、真實可靠的測站位移序列或速度信息，為后續(xù)的地殼形變特征分析、斷層活動反演提供技術(shù)支撐。高精度的數(shù)據(jù)處理技術(shù)是重要的必備前提。通過對監(jiān)測網(wǎng)絡(luò)基準(zhǔn)站長年連續(xù)觀測得到的 GPS 數(shù)據(jù)進行精密處理和分析，可以獲得基準(zhǔn)站的坐標(biāo)時間序列，進而獲取研究區(qū)域高精度速度場。利用速度場數(shù)據(jù)通過區(qū)域動力學(xué)模型的正、反演數(shù)值模擬或計算，可進一步分析其地殼位移、應(yīng)力、應(yīng)變場等動力學(xué)特征，探討其構(gòu)造運動的驅(qū)動力源和動力學(xué)演化模式。然而，同一軟件處理結(jié)果間的比較，僅說明了在一致的處理方案與原則下，計算結(jié)果的穩(wěn)定性和內(nèi)符合度。同種軟件處理結(jié)果間符合得再好, 不能完全斷定最終結(jié)果的無偏, 更不能估計系統(tǒng)性偏差的大小和誤差范圍。而不同軟件處理結(jié)果間比對,則可在一定程度上彌補以上缺憾,是一種比較客觀的評價方法。美國MIT的GAMIT/GLOBK軟件，瑞士伯爾尼科技大學(xué)的BERNESE軟件作為兩款最常用的GPS高精度數(shù)據(jù)處理軟件，這些軟件數(shù)據(jù)處理主要分為兩部分：一是對GPS原始數(shù)據(jù)進行解算獲得同步觀測網(wǎng)的基線解；二是對各同步網(wǎng)解進行整體平差和分析，獲得GPS網(wǎng)的整體解。盡管數(shù)據(jù)處理基本原理相同, 但在具體做法上和處理理念上差異明顯。兩種軟件結(jié)果的比較是檢驗網(wǎng)絡(luò)工程最終結(jié)果的一種有效手段。無論使用兩者中的哪一個，都可以利用另一個作為外部參考來進行驗證。只有這樣，才能客觀評價解算結(jié)果的可靠性。中國大陸川滇地區(qū)受青藏高原地殼物質(zhì)的東向移動、阿薩姆頂點楔入及緬甸微板塊的共同影響，成為現(xiàn)今地殼運動最活躍的地區(qū)之一。而該地區(qū)活動斷裂的特殊構(gòu)造特性, 又使得該地區(qū)的地殼運動表現(xiàn)為復(fù)雜多變的特征。為此，長期以來，川滇地區(qū)被視為重要的地震監(jiān)測區(qū)域，并取得了大量的研究成果[1-6]。本文分別利用GAMIT/GLOBK10.4和BERNESE5.2兩款軟件對川滇地區(qū)2010—2014年間中國地殼運動觀測網(wǎng)絡(luò)基準(zhǔn)站數(shù)據(jù)進行了解算，獲取了川滇地區(qū)相對歐亞板塊運動速度場，深入分析了現(xiàn)今川滇地區(qū)地殼運動特征和主要斷裂帶活動情況。

1 解算方案

為了能客觀地比較GAMIT/GLOBK10.4和BERNESE5.2兩種軟件的處理結(jié)果，盡量采用相同的處理參數(shù)，設(shè)置相同的參考框架ITRF2008、慣性系J2000.0、衛(wèi)星截止高度角10°、IGS精密星歷和選用一樣的參考框架站等。

GAMIT/GLOBK10.4軟件解算方案：利用雙差觀測值，組成與觀測值和參數(shù)相關(guān)的非線性數(shù)學(xué)模型，整個GNSS網(wǎng)全部基線一起參與解算。解算中采用了最小二乘算法反復(fù)迭代來估計測站的相對位置、軌道、地球自轉(zhuǎn)參數(shù)、對流層天頂延遲參數(shù)和大氣水平梯度參數(shù)等，得到的載波相位整周模糊度分別為實數(shù)和整數(shù)的約束解與松弛解。選取了11個中國及周邊的IGS站聯(lián)合進行處理。本文參考框架點選擇的依據(jù)是：①全球均勻分布；②站點穩(wěn)定，不受大震或大的構(gòu)造活動的影響；③水平向重復(fù)性誤差不超過3 mm；④站點線性趨勢項明顯，一致性較好等。對IGS核心站設(shè)置強約束，水平向0.05 m，垂直向0.1 m?；鶞?zhǔn)站采用BERNESE PPP定位結(jié)果作為先驗值，水平向1 m，垂直向1 m。主要模型選擇與參數(shù)設(shè)置如下：衛(wèi)星截止高度角設(shè)置10°，觀測值隨高度角定權(quán)。根據(jù)VMF1模型修正大氣層的映射函數(shù)。根據(jù)分段線性模型估計對流層天頂延遲和水平大氣梯度，每2 h估計1個天頂延遲參數(shù)?；€處理是利用偽距和載波相位觀測資料的雙差組合求得臺站坐標(biāo)和衛(wèi)星軌道的單日松弛解。由于基線向量無法提供確定點的絕對坐標(biāo)所必需的絕對位置基準(zhǔn)，GPS衛(wèi)星的軌道和測站坐標(biāo)都不在同一個穩(wěn)定的參考框架里，整個GPS網(wǎng)作為一個剛體每次解算都會存在整體平移和旋轉(zhuǎn)。因此，必須引入外部基準(zhǔn)。先利用雙差數(shù)據(jù)進行最小二乘參數(shù)估計，解算出各時段的基線和模糊度，然后將各同步觀測網(wǎng)自由基準(zhǔn)的法方程矩陣進行疊加，再對平差系統(tǒng)給予確定的基準(zhǔn)，獲得最終的平差結(jié)果。采用卡爾曼濾波的模型，對同步網(wǎng)解進行整體處理，獲取的測站坐標(biāo)和速度。利用GLOBK將整個監(jiān)測網(wǎng)每天的單日松弛解和圣地亞哥海洋研究所軌道中心SOPAC計算出的全球IGS跟蹤站的多個單日松弛解(igs1～igs6)合并，得到一個包含GPS 測站全球分布的合并單日松弛解。最后，以全球單日松弛解做為準(zhǔn)觀測值，利用GLOBK進行卡爾曼濾波參數(shù)估計，并通過IGS站作為參考站進行赫爾默特七參數(shù)相似變換，得到ITRF08參考框架下的測站坐標(biāo)和速率。

BERNESE5.2軟件解算方案：采用GPSEST程序計算基線浮動解，在此基礎(chǔ)上，基于最小二乘估計準(zhǔn)則，解算獨立基線。采用QIF(Quias_Ionosphere Free)方式，忽略基線之間的相關(guān)性，逐條確定基線中的相位模糊度。然后，顧及基線之間的相關(guān)性，根據(jù)相位雙差觀測方程，利用GPSEST程序計算獲得GPS測站的坐標(biāo)自由解。最后，選擇合適的框架基準(zhǔn)，選擇ITRF框架或局部框架，采用重心基準(zhǔn)平差的方法，利用ADDNEQ2程序?qū)⒂^測網(wǎng)依附至特定框架系中，獲得觀測網(wǎng)中GPS測站的坐標(biāo)約束解。在對ITRF或IGS先驗坐標(biāo)作為約束時，設(shè)定閾值：N向0.01 m，E向0.01 m，U向0.03 m。根據(jù)約束閾值，從框架站列表中自動剔除重心基準(zhǔn)下殘差超限的測站。

在GPS高精度地殼形變數(shù)據(jù)處理中，為了得到準(zhǔn)確可靠的結(jié)果，必須對結(jié)果進行驗證。因此，采用了GAMIT10.4和BERNESE5.2兩種軟件不同算法對同一數(shù)據(jù)進行了解算，通過對兩者結(jié)果時序趨勢的一致性和離散度的比較，來判斷解算結(jié)果的質(zhì)量與可靠性。

2 解算結(jié)果對比分析

2.1 基線解算結(jié)果

基線處理是GNSS地殼形變高精度數(shù)據(jù)處理中最關(guān)鍵的環(huán)節(jié)之一，其結(jié)果質(zhì)量好壞直接決定著最終測站坐標(biāo)、速度和地殼形變信息提取的成果。在基線解算過程中，由多臺GNSS接收機在野外通過同步觀測所采集到的觀測數(shù)據(jù)，被用來確定接收機間的基線向量及其方差-協(xié)方差陣?；€解算結(jié)果除了被用于檢驗外業(yè)觀測質(zhì)量和后續(xù)的網(wǎng)平差外，更重要的是，因基線向量提供了點與點之間的相對位置關(guān)系，并且與解算時所采用的衛(wèi)星星歷同屬于一個參考系，不受參考框架和共模誤差的影響，可確定GNSS網(wǎng)高精度的幾何形狀和基線變化量，對捕捉地殼微動態(tài)形變信息有著巨大的優(yōu)勢。利用GAMIT10.4和BERNESE5.2兩種軟件分別解算了川滇地區(qū)2010—2014年中國地殼運動觀測網(wǎng)絡(luò)51個GPS連續(xù)參考站基線單日解序列，其中的兩條基線SCLH-SCXJ、SCDF-SCYX的解算結(jié)果見圖1。由圖1，兩者DX、基線長度DL離散度均在5 mm之內(nèi)，DY、DZ離散度在10 mm之內(nèi)。兩種軟件解算結(jié)果整體吻合度較高，時序曲線趨勢變化一致。可見，兩種軟件基線解算結(jié)果基本能夠滿足高精度地殼形變監(jiān)測的要求。

圖1 基線SCLH-SCXJ、SCDF-SCYX三分量與基線長度單日解時間序列

2.2 坐標(biāo)解算結(jié)果

利用GAMIT10.4和BERNESE5.2解算了川滇地區(qū)陸態(tài)網(wǎng)62個GPS連續(xù)參考站2010—2014共5 a的單日解坐標(biāo)時間序列。其中的4個測站SCLH、SCXJ、SCDF和SCYX的解算結(jié)果見圖2。由圖2可見，兩種軟件網(wǎng)平差后單日解坐標(biāo)時間序列離散度接近，基本都在10 mm之內(nèi)，中誤差均優(yōu)于3 mm。兩者坐標(biāo)序列及其中誤差的整體趨勢變化一致，吻合度較高。雖BERNESE解算中誤差整體略低于GAMIT，但這并不說明中誤差小者解算精度就高，因為可能兩者解算中誤差之間存在有系統(tǒng)偏差。且同一軟件之間的相互比較只能說明處理方案內(nèi)部的一致性。而不同軟件，不同方法的解算結(jié)果之間的比較，能更客觀地評價數(shù)據(jù)處理結(jié)果的外部一致性與可靠性。

2.3 精度評定

單日解的重復(fù)率是衡量GPS觀測質(zhì)量與處理結(jié)果的重要指標(biāo)之一，其值越小，說明內(nèi)符合精度越高，解算質(zhì)量越好?；€重復(fù)率計算公式為

(1)

式中：Rc為基線的重復(fù)性統(tǒng)計值；n為同一基線的總的觀測時段數(shù)；Ci為一個時段所求得的基線某一分量或邊長；σci為相應(yīng)于Ci分量的中誤差；Cm為基線結(jié)果的加權(quán)平均值。

按最小二乘擬合出基線重復(fù)率和基線長度之間的線性關(guān)系，得出固定誤差和比例誤差部分，作為衡量基線精度的參考標(biāo)準(zhǔn)，見下式：

圖2 測站SCLH、SCXJ、SCDF和SCYX單日解坐標(biāo)時間序列及中誤差

Rc=a+bLc.

(2)

其中：Lc為基線長度；a為常數(shù)部分，是與基線長度無關(guān)的誤差部分；b為比例部分，是與基線長度成正比的比例誤差。

圖3 基線重復(fù)率

解算川滇地區(qū)各基線，利用單天解分別計算基線重復(fù)率。選取一周內(nèi)的重復(fù)率結(jié)果見圖3。其中基線重復(fù)率的固定部分單位為 mm，比例部分的基線長度L單位為m。由圖3可見，單天解的相對精度達到了10-9量級。選擇網(wǎng)中最長的基線SCXJ-SCMB，長度為264 428.667 m。經(jīng)計算，N、E、U向和長度L的重復(fù)率分別為0.77 mm、1.14 mm、3.54 mm和1.04 mm。這說明對于最長的基線來說，其絕對重復(fù)率也能達到水平向1 mm左右、垂向3 mm左右的水平。

2.4 測站速度解算

對于長時間連續(xù)GPS觀測，為獲取GPS測站精確的速度場信息，對GPS長時間坐標(biāo)時間序列進行趨勢項擬合獲取測站速度[7-8]。由于閃爍噪聲和隨機漫步噪聲均屬于冪律噪聲類型，而且原始時間序列的最優(yōu)噪聲模型與PL+VW噪聲模型的極大似然值相差甚微。而且，冪律噪聲的水平可以反映時間序列的整體噪聲水平。因此，分別利用白噪聲模型和PL+VW噪聲模型進行了速度估計。

結(jié)果表明，所有測站殘差坐標(biāo)時間序列的有色噪聲功率譜指數(shù)基本處于(-1，0)區(qū)間內(nèi)?？梢?，5 a的數(shù)據(jù)主要表現(xiàn)為閃爍噪聲。由于隨機漫步噪聲的長延時相關(guān)性，而川滇地區(qū)連續(xù)站坐標(biāo)時間序列的時間過短，僅5 a時間，很可能還不足以準(zhǔn)確量化隨機漫步噪聲的特征。不考慮時間序列中的有色噪聲會導(dǎo)致擬合速度精度的嚴(yán)重高估，顧及有色噪聲的影響，GPS測站擬合速度值變化不大，但是擬合速度中誤差將擴大約2～8倍。

3 速度場分析

選取歐亞板塊剛性穩(wěn)定臺站作為測站速度計算的參考框架，解算歐亞板塊歐拉參數(shù)，將解算所得的ITRF2008框架下的速度場轉(zhuǎn)換為穩(wěn)定的歐亞板塊參考框架下的相對速度場[9-10]。圖4為2010—2014年間由GPS獲得的川滇地區(qū)相對歐亞板塊地殼運動水平速度場。東西向、南北向速度估算精度優(yōu)于1 mm/a。

圖4 2010—2014年相對歐亞板塊川滇地區(qū)GPS水平運動速度場

由圖4可見，除滇西南地區(qū)呈現(xiàn)出由西向東的增大趨勢外，川滇塊體其它區(qū)域地殼活動由北往南、由西往東呈逐漸衰減趨勢。在鮮水河斷裂帶的東北地區(qū)，測站運動速度較大，平均約為20 mm/a，運動方向為東向和東南向；滇南地區(qū)測站運動速度相對較小，平均約為7 mm/a。滇東南，運動方向為南向和東南向，滇西南運動方向為南向和南西向，在滇南中部紅河斷裂帶的南部地區(qū)，測站運動方向則表現(xiàn)為南向。因測站SCNC、SCSN、LUZH同位于穩(wěn)定的四川盆地內(nèi)部， 從而呈現(xiàn)出三測站運動速度較小， 約5 mm/a， 且大小相同，運動方向一致。顯然處于同一穩(wěn)定塊體。由以上分析可以得出，川滇地塊運動特征主要表現(xiàn)為持續(xù)圍繞喜馬拉雅東結(jié)點作順時針旋轉(zhuǎn)，且仍然存在著南北向的擠壓特征。從地質(zhì)動力學(xué)分析: 印度板塊與歐亞板塊的碰撞、擠壓是川滇地塊巖石層水平形變的主要驅(qū)動力[1]，從而引起了青藏高原的隆升。在重力勢能作用下，造成青藏高原物質(zhì)東向擠出，青藏高原物質(zhì)的東向擠壓造成了川滇地塊的東移[3]。再在甘青地塊向南的擠壓力及相對穩(wěn)定的華南地塊的阻擋下，青藏高原東南部物質(zhì)相對歐亞板塊轉(zhuǎn)向南東方向運動，川滇地塊總體向東南方向順時針旋轉(zhuǎn)運動，區(qū)域斷裂帶呈旋走滑運動特征。由上述川滇地塊內(nèi)中國地殼運動觀測網(wǎng)絡(luò)GPS基準(zhǔn)站水平運動速度得到的川滇地塊整體運動特征與地質(zhì)動力學(xué)分析的運動結(jié)果完全相符。

斷裂活動是地殼運動的主要形式。一般來說，絕大多數(shù)地塊內(nèi)部比較穩(wěn)定，構(gòu)造變形主要發(fā)生在地塊之間的斷裂上。地塊之間的相互作用，主要通過其間的斷裂活動來實現(xiàn)[11]。為此，分別以安寧河-則木河斷裂、小江斷裂和紅河斷裂為界，選用了兩側(cè)能反映震間斷裂構(gòu)造特征的GPS測站，分析了斷層兩側(cè)的運動速度差異。利用實測垂直斷層剖面平行斷層的速度分量分析了斷裂帶可能的變形寬度[12]。圖5給出了2010—2014年期間三斷裂地表GPS站平行斷層的速度剖面，利用反正切函數(shù)擬合其變化趨勢，確定各斷裂帶可能的變形寬度。其中，安寧河-則木河斷裂兩測運動速度差異較大，為11.4 mm/a，閉鎖寬度約為200 km；小江斷裂兩測運動速度差異為9.5 mm/a，閉鎖寬度約為200 km。紅河斷裂兩測運動速度差異較小，約為6.4 mm/a，變形寬度較大。因研究的斷裂帶運動特征以走滑為主，因此，利用平行斷層的速度剖面可以得到斷裂帶可能的變形寬度。從而看出，小江斷裂帶、安寧河-則木河斷裂帶處于明顯的震間閉鎖狀態(tài)。而紅河斷裂的震間期變形特征與其他斷裂帶變形特征相比，閉鎖狀態(tài)不很清晰。

圖5 2010—2014年各斷裂帶平行斷層速度剖面

4 結(jié) 論

兩種軟件不同方法的解算結(jié)果之間的比較，能更客觀地評價數(shù)據(jù)處理結(jié)果的外部一致性與可靠性，是檢驗網(wǎng)絡(luò)工程最終結(jié)果的一種有效手段。通過兩種軟件對2010—2014年間川滇地區(qū)中國地殼運動觀測網(wǎng)絡(luò)GPS基準(zhǔn)站處理結(jié)果表明，兩種軟件網(wǎng)平差后單日解坐標(biāo)時間序列離散度接近，基本都在10 mm之內(nèi)，中誤差均優(yōu)于3 mm。兩者坐標(biāo)序列及其中誤差的整體趨勢變化一致，吻合度較高。由獲取的川滇地區(qū)相對歐亞板塊的GPS速度場得出，川滇地區(qū)的地殼運動除滇西南地區(qū)呈現(xiàn)出由西向東的增大趨勢外，川滇塊體其它區(qū)域地殼活動由北往南、由西往東呈逐漸衰減趨勢。運動方向主要表現(xiàn)為持續(xù)圍繞喜馬拉雅東結(jié)點作順時針旋轉(zhuǎn)，且仍然存在著南北向的擠壓特征。在鮮水河斷裂帶的東北地區(qū)，測站運動速度較大，平均約為20 mm/a，運動方向為東向和東南向；滇南地區(qū)測站運動速度相對較小，平均約為7 mm/a。滇東南，運動方向為南向和東南向，滇西南運動方向為南向和南西向。在滇南中部紅河斷裂帶的南部地區(qū)，測站運動方向則表現(xiàn)為南向。因測站SCNC、SCSN、LUZH，同位于穩(wěn)定的四川盆地內(nèi)部，從而呈現(xiàn)出三測站運動速度較小，約5 mm/a，且大小相同，運動方向一致，顯然處于同一穩(wěn)定塊體。

以安寧河-則木河斷裂、小江斷裂和紅河斷裂為界，選用了兩側(cè)能反映震間斷裂構(gòu)造特征的GPS測站，獲取了2010—2014年期間三斷裂帶地表GPS站平行斷層的速度剖面，利用反正切函數(shù)擬合其變化趨勢，確定了各斷裂帶可能的變形寬度，分析了現(xiàn)今川滇地區(qū)震間主要斷裂帶的活動特征。其中，安寧河-則木河斷裂兩測運動速度差異較大，為11.4 mm/a，閉鎖寬度約為200 km；小江斷裂兩測運動速度差異為9.5 mm/a，閉鎖寬度約為200 km。紅河斷裂兩測運動速度差異較小，約為6.4 mm/a，變形寬度較大。兩測相對速度越大，變形寬度越大，也就是說其閉鎖程度最高，鎖定深度最深，地震危險性也較高。因研究的斷裂帶運動特征以走滑為主。從所得的斷裂帶可能的變形寬度可以看出，小江斷裂帶、安寧河-則木河斷裂帶處于明顯的震間閉鎖狀態(tài)，很好地反映了震間斷裂帶的活動性。而紅河斷裂的震間期變形特征與其他斷裂帶變形特征相比，閉鎖狀態(tài)不很清晰?？赡艿脑驗椋孩偌t河斷裂帶附近的其他斷裂系統(tǒng)的影響導(dǎo)致了斷裂帶兩測變形的特征；②跨斷層測站數(shù)量太少，統(tǒng)計特征不明顯。因此，需要進一步結(jié)合更多的GPS測站，包括區(qū)域站、連續(xù)站，來重點刻畫主要斷裂帶的活動特征及運動規(guī)律。

感謝中國地震局提供的陸態(tài)網(wǎng)觀測數(shù)據(jù)及吳偉偉提供的BERNESE解算結(jié)果。

[1] 李延興,楊國華,李智,等.中國大陸活動地塊的運動與應(yīng)變狀態(tài)[J].中國科學(xué)(D輯:地球科學(xué)),2003,33(增1):65-81.

[2] 喬學(xué)軍,王琪,杜瑞林.川滇地區(qū)活動地塊現(xiàn)今地殼形變特征[J].地球物理學(xué)報,2004(5):806-812.

[3] 徐錫偉,聞學(xué)澤,鄭榮章,等.川滇地區(qū)活動塊體最新構(gòu)造變動樣式及其動力來源[J].中國科學(xué)(D輯:地球科學(xué)),2003,33(增1):151-162.

[4] 魏文薪.川滇塊體東邊界主要斷裂帶運動特性及動力學(xué)機制研究[D].北京：中國地震局地質(zhì)研究所,2012.

[5] 王閻昭,王恩寧,沈正康,等.基于GPS資料約束反演川滇地區(qū)主要斷裂現(xiàn)今活動速率[J].中國科學(xué)(D輯:地球科學(xué)),2008,38(5):582-597.

[6] 徐克科,伍吉倉,吳偉偉.基于GNSS時空數(shù)據(jù)的瞬態(tài)無震蠕滑信息檢測[J].地球物理學(xué)報,2015,58(7):2330-2338.

[7] DAVIS J L,BENNETT R A,WERNICKE B P.Assessment of GPS velocity accuracy for the Basin and Range Geodetic Network (BARGEN)[J].Geophysical research letters,2003,30(7).

[8] WILLIAMS S D P.The effect of coloured noise on the uncertainties of rates estimated from geodetic time series[J].Journal of Geodesy,2003,76(9-10):483-494.

[9] 徐克科,伍吉倉.利用抗差估計進行剛性板塊異常臺站探測及其運動參數(shù)求取[J].大地測量與地球動力學(xué),2014,34(2):95-99.

[10] 徐克科,伍吉倉,呂志鵬.基于假設(shè)檢驗的剛性板塊臺站選取及全球板塊運動模型的建立[J].大地測量與地球動力學(xué),2013,33(6):121-125.

[11] ENGLAND P, MOLNAR P. Right-lateral shear and rotation as the explanation for strike-slip faulting in eastern Tibet [J].Nature,1990, 344:140-142.

[12] SAVAGE J C,SVARC J I,PRESCOTT W H.Geodetic estimates of fault slip rates in the San Francisco Bay area[J].J Geophys Res.1999,104(B3):4995-5002.

[責(zé)任編輯：劉文霞]

Caculation and analysis of present-day crustal movement velocity fieldin Sichuan-Yunnan block

XU Keke1,WU Jicang2

(1.School of Surveying and Land Information Engineering,Henan Polytechnic University,Jiaozuo 454000,China;2.College of Surveying and Geo-Information,Tongji University, Shanghai 200092,China)

Based on GAMIT and BERNESE softwares,the GPS data of 51 base stations from 2010-2014 from Crustal Movement Observation Network of China is processed and analysed.The discrete degree of both baseline and coordinate time series are better than 10 mm.The change tendency is consistent and the coincidence degree is high.The baseline repeatability is better than 1+4×10-9Lin horizontal direction and 4+6×10-9Lin vertical direction.Crustal movement velocity field of Sichuan-Yunnan region relative to Eurasia plate is obtained.The results show that except the southwest of Yunnan which velocity increases from west to east,the crustal movement characteristics of other Sichuan-Yunnan region gradually reduces from north to south, and from west to east.Still,movement direction shows clockwise rotation around the Himalayan east nodes,and strong extrusion characteristics from north and south.Considering Anning River-Zemu River fault and Xiaojiang fault and Honghe fault as the boundary respectively, the GPS station reflecting fault characteristics on both sides are chozen. The velocity difference on both sides of fault is analyzed.By the velocity profile paralleling with the fault,the arctangent function trend is fitted and the deformation width of the fault zone is obtained.

GPS;Sichuan-Yunnan block;crustal movement;velocity field;fault zone

10.19349/j.cnki.issn1006-7949.2017.03.002

2016-02-04

國家973重點基礎(chǔ)研究發(fā)展計劃(2013CB733304)；國家自然科學(xué)基金資助項目(41404023)

徐克科(1979-)，男，副教授，博士.

P228.4

A

1006-7949(2017)03-0008-06

引用著錄：徐克科,伍吉倉.川滇地塊現(xiàn)今地殼運動速度場解算與分析[J].測繪工程，2017，26(3)：8-13,18.