基于GPS資料研究云南地區(qū)地殼形變動(dòng)態(tài)特征

王巖 洪敏 邵德勝 汪志民

摘要：基于1999—2007年和2011—2013年2期云南地區(qū)的GPS速度場(chǎng)資料，結(jié)合最小二乘配置方法給出了2個(gè)時(shí)段應(yīng)變參數(shù)和GPS速度剖面結(jié)果，綜合分析了汶川8.0級(jí)地震前后云南地區(qū)地殼變形動(dòng)態(tài)演化特征。GPS 速度場(chǎng)動(dòng)態(tài)演化結(jié)果表明：整體上GPS速度場(chǎng)方向從川滇塊體南部的南向逐漸轉(zhuǎn)為滇南塊體的南東向，塊體邊界斷裂帶附近存在明顯的相對(duì)運(yùn)動(dòng)特征。相比于1999—2007年，點(diǎn)位更加密集的2011—2013年GPS速度場(chǎng)結(jié)果還反映了斷裂帶分段變形特征；在滇南塊體及周邊地區(qū)，南東向的運(yùn)動(dòng)有顯著的增加。應(yīng)變參數(shù)動(dòng)態(tài)演化結(jié)果表明：整體上呈現(xiàn)“中部拉張、兩端擠壓”的特點(diǎn)；區(qū)域變形（剪切、拉張或擠壓）的分布特征與區(qū)域主要斷裂帶的背景運(yùn)動(dòng)及變形特征密切相關(guān)；相比于1999—2007年，2011—2013年的拉張變形區(qū)有所擴(kuò)大，并且有向東和向南擴(kuò)展的趨勢(shì)；川滇塊體東邊界（小江斷裂帶）整體的擠壓變形已經(jīng)不再顯著，主要集中在小江斷裂帶北段靠北的區(qū)域?？鐢嗔褞PS速度剖面結(jié)果顯示：小江斷裂帶左旋滑動(dòng)速率從北向南逐漸減?。ㄓ?0 mm/a減至5 mm/a），北段變形寬度較寬；紅河斷裂帶右旋滑動(dòng)速率約為4 mm/a，變形寬度較寬；動(dòng)態(tài)結(jié)果表明，汶川地震對(duì)2條斷裂帶變形模式影響較小。

關(guān)鍵詞：最小二乘配置；GPS速度場(chǎng)；應(yīng)變參數(shù)；GPS速度剖面

中圖分類(lèi)號(hào)：P315.725 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1000-0666（2018）03-0368-07

0 引言

云南地區(qū)地處青藏高原東南緣地區(qū)，是我國(guó)地震活動(dòng)最為顯著的地區(qū)之一，該區(qū)域活動(dòng)構(gòu)造豐富且復(fù)雜（張俊昌等，1980；王義昭等，1988；李西，2015），這種特殊的構(gòu)造環(huán)境也造就了該區(qū)域特殊的地殼變形特征。本區(qū)域位于南北地震南段，多位研究學(xué)者針對(duì)南北地震帶地區(qū)地殼變形的特征開(kāi)展了大量研究（江在森等，2001；方穎等，2008；楊國(guó)華等，2009；丁開(kāi)華等，2013；張希等，2003；武艷強(qiáng)等，2012；魏文薪等，2012；季靈運(yùn)等，2015），另外，也有學(xué)者針對(duì)云南地區(qū)的地殼變形特征開(kāi)展了相關(guān)的研究（楊國(guó)華等，2003；洪敏等，2014；王伶俐等，2015，周海濤等，2017），但針對(duì)云南地區(qū)的地殼變形的動(dòng)態(tài)演化特征研究還略顯不足。2008年汶川8.0級(jí)地震發(fā)生對(duì)云南地區(qū)地殼變形的影響如何，該地震前后區(qū)域的變形是否有顯著變化等問(wèn)題值得進(jìn)一步研究。因此，本文基于1999—2007年和2011—2013年云南地區(qū)的GPS速度場(chǎng)資料，結(jié)合最小二乘配置方法給出2個(gè)時(shí)段應(yīng)變參數(shù)結(jié)果和GPS速度剖面結(jié)果，綜合分析了汶川地震前后云南地區(qū)地殼變形動(dòng)態(tài)演化特征。

1 GPS資料概況及數(shù)據(jù)處理

我國(guó)從20世紀(jì)80年代末逐步開(kāi)始GPS觀測(cè)，較大規(guī)模、高精度的GPS觀測(cè)從 “中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)”建設(shè)完成后開(kāi)始，該網(wǎng)絡(luò)分布于中國(guó)大陸的GPS區(qū)域站約1 000個(gè)。在中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)基礎(chǔ)上，2007年進(jìn)行二期建設(shè)的“中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)”（簡(jiǎn)稱(chēng)陸態(tài)網(wǎng)絡(luò)）把區(qū)域網(wǎng)點(diǎn)增加到2 000個(gè)。為研究云南地區(qū)地殼變形的動(dòng)態(tài)演化特征，本文主要利用中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)建成后的GPS復(fù)測(cè)資料進(jìn)行研究，這些資料包括：中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)GPS區(qū)域網(wǎng)1999，2001，2004，2007年復(fù)測(cè)資料，中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)GPS區(qū)域網(wǎng)2011，2013年復(fù)測(cè)資料。

GPS區(qū)域網(wǎng)速度場(chǎng)精密解算工作與GPS時(shí)間序列解算類(lèi)似，分2步進(jìn)行。第一步為GAMIT解算，由于同步觀測(cè)的測(cè)點(diǎn)較多，考慮到GAMIT單次解算的測(cè)點(diǎn)有限，一般采用分區(qū)解算方式，此步的主要產(chǎn)出為單日松弛解hfile。第二步解算以單日松弛解hfile作為輸入，融合中國(guó)大陸GPS連續(xù)站和IGS單日松弛解結(jié)果，采用GLOBK軟件或QOCA軟件（Dong et al，1998）進(jìn)行GPS速度解算。由于GPS速度解算至少需要2期GPS觀測(cè)，因此統(tǒng)一2期（多期）解算結(jié)果的參考基準(zhǔn)就顯得尤為重要，此時(shí)參考框架點(diǎn)的選取原則與時(shí)間序列類(lèi)似。另外一個(gè)需要關(guān)注的問(wèn)題是同震位移的扣除，如果在地震后進(jìn)行過(guò)GPS觀測(cè)，則可以利用觀測(cè)數(shù)據(jù)建立同震位移分布模型，對(duì)同震影響進(jìn)行扣除。

基于上述處理原則，采用GAMIT/GLOBK軟件解算，得出了1999—2007年共4期和2011—2013年共2期云南地區(qū)GPS速度場(chǎng)結(jié)果?？紤]到南北地震帶中南部西側(cè)為中國(guó)大陸最大的相對(duì)穩(wěn)定的華南地塊，在云南地區(qū)的GPS站實(shí)測(cè)速度中，扣除華南地塊的歐拉剛體運(yùn)動(dòng)參數(shù)確定的理論值，獲得該區(qū)域相對(duì)于華南地塊的地殼運(yùn)動(dòng)速度圖，既能夠反映云南地區(qū)屬于青藏地塊區(qū)的強(qiáng)烈變形與穩(wěn)定的華南地塊的差異，也能夠兼顧反映云南地區(qū)主要塊體邊界帶的相對(duì)運(yùn)動(dòng)。筆者還針對(duì)制約多期GPS速度場(chǎng)動(dòng)態(tài)演化特征分析的關(guān)鍵問(wèn)題，利用改進(jìn)擬準(zhǔn)檢定（QUAD）初選指標(biāo)的方法，有效解決多期速度場(chǎng)的基準(zhǔn)偏移對(duì)粗差剔除的影響，給出了相對(duì)華南基準(zhǔn)的速度場(chǎng)結(jié)果（鄒鎮(zhèn)宇等，2015）。圖1給出了由上述GPS觀測(cè)資料解算獲得的1999—2007年及2011—2013年云南地區(qū)相對(duì)于華南地塊的GPS速度場(chǎng)結(jié)果。

2 地殼相對(duì)運(yùn)動(dòng)及其動(dòng)態(tài)變化特征 從云南地區(qū)總體地殼相對(duì)運(yùn)動(dòng)與變形特征來(lái)看（圖1），青藏高原東部物質(zhì)向邊緣流動(dòng)受到邊界較穩(wěn)定塊體的阻擋，造成在接近垂直青藏地塊邊界帶方向地殼縮短與沿地塊邊界帶方向地殼伸展，兼有扭動(dòng)與剪切變形的特征。根據(jù)速度衰減或方向改變，大致為南東向或南向。作為本區(qū)域重要的塊體邊界帶的小江斷裂帶GPS站速度矢量大小方向在空間上的不斷變化，顯示有較強(qiáng)的相對(duì)運(yùn)動(dòng)與變形。這表明由于其特殊的構(gòu)造環(huán)境及其東側(cè)被穩(wěn)定的地塊阻擋，其對(duì)青藏地塊向東運(yùn)動(dòng)及物質(zhì)東移的吸收十分顯著。云南地區(qū)以近南北向地殼縮短變形、近東西向地殼伸展為主。在最南部南北向800 km以?xún)?nèi)最大縮短量約8 mm/a，而在最南端近東西向地殼伸張量為近10 mm/a。該區(qū)地殼相對(duì)運(yùn)動(dòng)與變形以川滇地塊向南擠出滑移最為顯著，且川滇地塊中部運(yùn)動(dòng)速度最大。南部整體上呈現(xiàn)順時(shí)針扭轉(zhuǎn)運(yùn)動(dòng)。GPS站速度從西到東及從北到南不斷向南、南部偏轉(zhuǎn)，使站速度從近東西向運(yùn)動(dòng)偏轉(zhuǎn)到川滇地塊南側(cè)地區(qū)為向南或南西運(yùn)動(dòng)。在南北地震帶南部東邊緣，即川滇地塊與華南地塊邊界帶的GPS站速度在平行地塊邊界帶的分量自西向東衰減顯著，表現(xiàn)為明顯的左旋變形特征，總體左旋扭動(dòng)量約為10～12 mm/a。

由圖1a可見(jiàn)，1999—2007年 GPS速度場(chǎng)特征主要體現(xiàn)在：（1）GPS速度場(chǎng)方向從川滇塊體南部南向逐漸轉(zhuǎn)為在滇南塊體的南西向，反映出了繞喜馬拉雅東構(gòu)造結(jié)做順時(shí)針旋轉(zhuǎn)的變形特征；（2）在川滇塊體東邊界的南段存在顯著的左旋剪切運(yùn)動(dòng)，滇西北地區(qū)存在顯著的拉張變形；（3）青藏高原物質(zhì)東流，受穩(wěn)定的華南塊體的阻擋，使得云南地區(qū)形成現(xiàn)有的變形特征，雖然區(qū)域的構(gòu)造十分復(fù)雜，但主要的塊體邊界斷裂帶的相對(duì)運(yùn)動(dòng)特征仍能在速度場(chǎng)中清晰的顯示出來(lái)。由圖1b可見(jiàn)，2011—2013年GPS點(diǎn)位的密度顯著增加，但是由于觀測(cè)時(shí)間較短，誤差大于1999—2007年的結(jié)果；速度場(chǎng)結(jié)果不僅反映了塊體邊界帶的相對(duì)運(yùn)動(dòng)特征，也可以反映斷裂帶分段變形

特征（例如小江斷裂帶的北段和中段的左旋滑動(dòng)速率顯著的大于南段的滑動(dòng)速率）；并且，相比于1999—2007年速度場(chǎng)結(jié)果，在滇南塊體及周邊地區(qū)，2011—2013年南東向的運(yùn)動(dòng)有顯著的增加。

3 應(yīng)變率場(chǎng)動(dòng)態(tài)變化特征

雖然采用相對(duì)運(yùn)動(dòng)速度場(chǎng)可以較直觀地反映地殼變形特征，但就對(duì)地殼變形的定量描述方面來(lái)看，應(yīng)變場(chǎng)的分析更有優(yōu)勢(shì)。應(yīng)變場(chǎng)各種參數(shù)能夠全面表達(dá)變形的不同性質(zhì)與強(qiáng)度，從一般概念講，由于應(yīng)變量是由觀測(cè)點(diǎn)之間的相對(duì)運(yùn)動(dòng)來(lái)決定的，它不依賴(lài)于參考基準(zhǔn)。但目前應(yīng)變計(jì)算有多種方法，不同研究者采用不同的方法用相同來(lái)源的GPS資料所給出的應(yīng)變場(chǎng)圖像也有較明顯差異（石耀霖，朱守彪，2006）。Wu 等（2011）借助模擬實(shí)驗(yàn)數(shù)據(jù)對(duì)最小二乘配置應(yīng)變場(chǎng)解算方法結(jié)果與三角單元法、多面函數(shù)擬合和球諧函數(shù)展開(kāi)這幾種應(yīng)變場(chǎng)方法解算結(jié)果進(jìn)行了比較分析，根據(jù)其與理論值的逼近程度，客觀估計(jì)不同方法模型的優(yōu)缺點(diǎn)及其抗差性與實(shí)用性。結(jié)果表明，最小二乘方法最優(yōu)。因此本文利用最小二乘配置球面應(yīng)變率計(jì)算方法對(duì)研究區(qū)域的應(yīng)變率場(chǎng)特征進(jìn)行分析（江在森，劉經(jīng)南，2010）。

圖2給出了2個(gè)時(shí)段的面應(yīng)變率及最大主應(yīng)變率的分布圖，1999—2007年的結(jié)果顯示：（1）拉張變形主要集中在紅河斷裂帶、程海斷裂帶、麗江小金河斷裂帶南段等區(qū)域，擠壓變形主要集中在滇西及川滇塊體東邊界南段的區(qū)域，并且整體上呈現(xiàn)“中部拉張、兩端擠壓”的變形特征；（2）拉張變形的區(qū)域與地質(zhì)結(jié)果給出的斷裂帶變形特征具有一致性（劉光勛等，1986；黃小巾等，2014），擠壓變形特征可能與此區(qū)域的動(dòng)力源有關(guān)，即青藏高原物質(zhì)東流受到穩(wěn)定的華南塊體阻擋，在川滇塊體東邊界南段（小江斷裂帶）形成擠壓區(qū)；（3）在滇南的邊界區(qū)域形成了擠壓變形的高值區(qū)，由于最小二乘配置方法在計(jì)算的過(guò)程中存在邊界效應(yīng)，因此，高擠壓變形區(qū)可信度不高。2011—2013年的結(jié)果一方面繼承了上一期整體變形特征，另一方面還存在著顯著的差異性，主要表現(xiàn)為：拉張變形區(qū)有所擴(kuò)大，并且有向東和向南擴(kuò)展的趨勢(shì)；在靠近則木河斷裂帶附近存在顯著的擠壓變形的高值區(qū)，川滇塊體東邊界（小江斷裂帶）整體的擠壓變形已經(jīng)不再顯著，主要集中在小江斷裂帶北段靠北的區(qū)域。

圖3給出了最大剪應(yīng)變率的結(jié)果，1999—2007年的空間分布最顯著的特征是應(yīng)變率高值主要沿川滇地塊東邊界帶南段（小江斷裂帶）分布。并且以程海斷裂帶為界，形成了最大剪應(yīng)變率的高值與低值的過(guò)渡區(qū)域，可能預(yù)示著該區(qū)域剪切應(yīng)力積累水平較高，需關(guān)注該區(qū)域的強(qiáng)震危險(xiǎn)性。另外，區(qū)域剪切變形的分布特征與區(qū)域主要斷裂帶的背景運(yùn)動(dòng)及變形特征密切相關(guān)（例如左旋走滑運(yùn)動(dòng)為主的小江斷裂帶處于最大剪應(yīng)變率的高值區(qū)，右旋走滑運(yùn)動(dòng)為主的紅河斷裂帶處于最大剪應(yīng)變率的低值區(qū)）。從2011—2013年的空間分布特征來(lái)看，與上一時(shí)段的一致性較好，這也表明了地殼變形隨著時(shí)間的推移，具有較好的繼承性特征；另外，相比于上一時(shí)段，程海斷裂帶附近的最大剪應(yīng)變的高值區(qū)有向北擴(kuò)展的趨勢(shì)；與其同步的還有紅河斷裂帶的最大剪應(yīng)變的高值區(qū)也有向北擴(kuò)展的趨勢(shì)。

4 主要斷裂帶運(yùn)動(dòng)特性研究

由于最小二乘配置應(yīng)變率計(jì)算方法基于連續(xù)變形假設(shè)，因此其結(jié)果在識(shí)別具體斷裂帶的具體變形模式方面不具優(yōu)勢(shì)。GPS速度剖面圖像可直觀展現(xiàn)位移分布與斷層的關(guān)系，經(jīng)常被應(yīng)用于斷層變形特征的分析中（張培震等，2005，2008；Maurin et al，2010；Wu et al，2015）。因此，本文以小江斷裂帶和紅河斷裂帶為例分析斷裂帶運(yùn)動(dòng)變形的動(dòng)態(tài)演化特征。

4.1 小江斷裂帶

小江斷裂帶北起巧家以北，南至建水東南，全長(zhǎng)約400 km，總體呈南北走向。根據(jù)其內(nèi)部結(jié)構(gòu)分成3段，分別為北段、中段和南段。在長(zhǎng)期的活動(dòng)過(guò)程中，曾經(jīng)歷壓、張、扭不同力學(xué)性質(zhì)的轉(zhuǎn)化，晚第四紀(jì)以來(lái)，斷裂以強(qiáng)烈左旋走滑為特征（宋方敏，汪一鵬，1998）。按照上述的分段原則，以斷裂帶為中心，左右兩側(cè)分別圈定一定范圍，利用該范圍內(nèi)的GPS實(shí)測(cè)速度場(chǎng)數(shù)據(jù)進(jìn)行高斯投影，求得各GPS點(diǎn)的沿?cái)嗔褞ё呦蚍较蚣按怪睌嗔褞ё呦蚍较虻乃俾?；并按照走滑斷層的震間變形模式對(duì)剖面進(jìn)行擬合，最后給出了2個(gè)時(shí)間段平行斷裂帶GPS速度剖面和擬合結(jié)果，如圖4所示。

由圖4可見(jiàn)：（1）小江斷裂帶北段2個(gè)時(shí)段的滑動(dòng)速率變化不大，均為10 mm/a，且表現(xiàn)為一定的應(yīng)變積累的變形特征，變形寬度較寬；（2）小江斷裂帶中段：相比于1999—2007年，2011—2013年的遠(yuǎn)場(chǎng)滑動(dòng)速率有所增加（從7 mm/a增加到10 mm/a），小江斷裂帶從中段開(kāi)始分為2支，現(xiàn)有的GPS點(diǎn)位密度還不能識(shí)別出東西2支的變形特征，所得到的滑動(dòng)速率結(jié)果應(yīng)為2支滑動(dòng)速率之和。雖然滑動(dòng)速率有所增加，并且從剖面的形態(tài)來(lái)看仍然表現(xiàn)為一定的應(yīng)變積累特征，但2個(gè)時(shí)段的變形寬度變化不大；相比于北段，中段的變形寬度較窄；（3）小江斷裂帶南段：相比于北段，該段的滑動(dòng)速率有了明顯的減少，約為5 mm/a，這可能與曲江石屏斷裂帶吸收了部分變形有關(guān)，從2個(gè)時(shí)段的動(dòng)態(tài)演化特征來(lái)看，變形寬度沒(méi)有顯著的變化，且變形寬度比小江斷裂帶北段窄。綜合上述分析可以看出，小江斷裂帶以左旋走滑運(yùn)動(dòng)為主，不同段落變形有著明顯的差異，滑動(dòng)速率也有所不同；根據(jù)走滑斷層的震間期的變形模型推測(cè)小江斷裂帶的應(yīng)變積累水平較高；另外，從各段落的動(dòng)態(tài)演化特征來(lái)看，小江斷裂帶不同段落在不同時(shí)段變化并不明顯，汶川地震同震及震后調(diào)整對(duì)此斷裂帶影響較小。

4.2 紅河斷裂帶中段

根據(jù)紅河斷裂帶幾何結(jié)構(gòu)特征，可以分為3個(gè)變形區(qū)，并且北西段和南東段存在多個(gè)次級(jí)斷裂（向宏發(fā)等，2004），顧及到實(shí)際GPS點(diǎn)位分布特征，本文僅對(duì)紅河斷裂的中段進(jìn)行GPS剖面分析。由于紅河斷裂帶是以右旋走滑為主，因此，圖5給出了平行斷裂帶的GPS速度剖面和擬合結(jié)果。

從GPS速度剖面及擬合結(jié)果可以看出，紅河斷裂帶中段主要以右旋走滑運(yùn)動(dòng)為主，滑動(dòng)速率約為4 mm/a，2個(gè)時(shí)段的滑動(dòng)速率變化不大；從斷裂帶的變形形態(tài)來(lái)看，存在著一定的應(yīng)變積累特征，以現(xiàn)有的GPS點(diǎn)位做剖面擬合可以看出，2個(gè)時(shí)段的變形寬度變化不大，但變形寬度較寬，可能預(yù)示著該段落可能處于孕震晚期；從以上分析也可看出，汶川地震的發(fā)生對(duì)此斷裂的變形影響較小。

5 討論與結(jié)論

基于1999—2007年和2011—2013年的2期GPS速度場(chǎng)資料，并結(jié)合最小二乘配置方法給出的2個(gè)時(shí)段應(yīng)變參數(shù)結(jié)果，綜合分析了汶川地震前后云南地區(qū)地殼變形動(dòng)態(tài)演化特征。獲得了以下主要認(rèn)識(shí)：

（1）2個(gè)時(shí)段GPS速度場(chǎng)的動(dòng)態(tài)結(jié)果表明：整體上GPS速度場(chǎng)方向從川滇塊體南部南向逐漸轉(zhuǎn)為在滇南塊體的南西向，反映出了繞喜馬拉雅東構(gòu)造結(jié)做順時(shí)針旋轉(zhuǎn)的變形特征。塊體邊界斷裂帶附近存在明顯的相對(duì)運(yùn)動(dòng)特征。2011—2013年點(diǎn)位的密度顯著增加，速度場(chǎng)結(jié)果不僅反映了塊體邊界帶的相對(duì)運(yùn)動(dòng)特征，也可以反映斷裂帶分段變形特征；相比于1999—2007年的GPS速度場(chǎng)結(jié)果，在滇南塊體及周邊地區(qū)，南東向的運(yùn)動(dòng)有顯著的增加。

（2）2個(gè)時(shí)段的應(yīng)變率場(chǎng)動(dòng)態(tài)結(jié)果表明，整體上“中部拉張、兩端擠壓”；區(qū)域變形（剪切、拉張或擠壓）的分布特征與區(qū)域主要斷裂帶的背景運(yùn)動(dòng)及變形特征密切相關(guān)；相比于1999—2007年，2011—2013年給出的拉張變形區(qū)有所擴(kuò)大，并且有向東和向南擴(kuò)展的趨勢(shì)；川滇塊體東邊界（小江斷裂帶）整體的擠壓變形已經(jīng)不再顯著，主要集中在小江斷裂帶北段靠北的區(qū)域；程海斷裂帶附近的最大剪應(yīng)變的高值區(qū)有向北擴(kuò)展的趨勢(shì)；與其同步的還有紅河斷裂帶的最大剪應(yīng)變的低值區(qū)也有向北擴(kuò)展的趨勢(shì)。

（3）跨斷裂帶GPS速度剖面結(jié)果顯示：小江斷裂帶左旋滑動(dòng)速率從北向南逐漸減?。ㄓ?0 mm/a減至5 mm/a），北段變形寬度較寬；紅河斷裂帶右旋滑動(dòng)速率約為4 mm/a，變形寬度較寬；動(dòng)態(tài)結(jié)果表明，汶川地震對(duì)2條斷裂帶變形模式影響較小。

GPS 數(shù)據(jù)處理得到中國(guó)地震局第一監(jiān)測(cè)中心武艷強(qiáng)研究員數(shù)據(jù)解算團(tuán)隊(duì)的技術(shù)支持，GPS的變形分析得到了中國(guó)地震局地震預(yù)測(cè)研究所魏文薪博士的有益討論，在此一并表示感謝。

參考文獻(xiàn)：

丁開(kāi)華，許才軍，鄒蓉，等.2013.利用GPS分析川滇地區(qū)活動(dòng)地塊運(yùn)動(dòng)與應(yīng)變模型[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），38（7）：822-827.

方穎，江在森，顧國(guó)華.2008.昆侖山MS8.1地震對(duì)川滇地區(qū)地殼運(yùn)動(dòng)影響分析[J].大地測(cè)量與地球動(dòng)力學(xué)，28（4）：25-30.

洪敏，張勇，邵德盛，等.2014.云南地區(qū)近期地殼活動(dòng)特征[J].地震研究，37（3）：367-372.

黃小巾，吳中海，李家存，等.2014.滇西北裂陷帶的構(gòu)造地貌特征與第四紀(jì)構(gòu)造活動(dòng)性[J].地質(zhì)通報(bào)，33（4）：578-593.

季靈運(yùn)，劉立煒，郝明.2015.利用InSAR技術(shù)研究滇西南鎮(zhèn)康—永德地區(qū)現(xiàn)今地殼形變特征[J].地震研究，38（1）：84-89.

江在森，劉經(jīng)南.2010.應(yīng)用最小二乘配置建立地殼運(yùn)動(dòng)速度場(chǎng)與應(yīng)變場(chǎng)的方法[J].地球物理學(xué)報(bào)，53（5）：1109-1117.

江在森，張希，崔篤信，等.2001.青藏塊體東北緣近期水平運(yùn)動(dòng)與變形[J].地球物理學(xué)報(bào)，44（5）：636-644.

李西.2015.川滇地塊云南地區(qū)不同發(fā)育階段邊界斷裂破裂特征研究[D].北京：中國(guó)地震局地質(zhì)研究所.

劉光勛，李方全，李桂榮.1986.我國(guó)滇西北地震活動(dòng)區(qū)的活動(dòng)構(gòu)造與應(yīng)力狀態(tài)[J].地震地質(zhì)，8（1）：1-14.

石耀霖，朱守彪.2006.用GPS位移資料計(jì)算應(yīng)變方法的討論[J].大地測(cè)量與地球動(dòng)力學(xué)，26（1）：1-8.

宋方敏，汪一鵬.1998.中國(guó)活斷層研究專(zhuān)輯——小江活動(dòng)斷裂帶[M].北京：地震出版社，257.

王伶俐，邵德盛，施發(fā)奇.2015.基于GPS的云南地區(qū)活動(dòng)地塊現(xiàn)今運(yùn)動(dòng)及應(yīng)變特征分析[J].地震地磁觀測(cè)與研究，36（1）：28-36.

王義昭，熊家鏞，林堯明.1988.云南地質(zhì)構(gòu)造的若干特點(diǎn)[J].云南地質(zhì)，7（2）：3-9.

魏文薪，江在森，武艷強(qiáng)，等.2012.利用GPS數(shù)據(jù)研究川滇塊體東邊界主要斷裂帶運(yùn)動(dòng)特性[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），37（9）：1041-1044.

武艷強(qiáng)，江在森，楊國(guó)華，等.2012.南北地震帶北段近期地殼變形特征研究[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），37（9）：1045-1048.

向宏發(fā)，韓竹軍，虢順民，等.2004.紅河斷裂帶大型右旋走滑運(yùn)動(dòng)與伴生構(gòu)造地貌變形[J].地震地質(zhì)，26（4）：597-610.

楊國(guó)華，韓月萍，楊博.2009.川滇地區(qū)地殼水平運(yùn)動(dòng)與變形場(chǎng)的演化特征及其機(jī)制討論[J].地震研究，32（3）：275-282.

楊國(guó)華，王琪，王敏，等.2003.云南地區(qū)現(xiàn)今地殼水平運(yùn)動(dòng)與變形特征[J].大地測(cè)量與地球動(dòng)力學(xué)，23（2）：7-14.

張俊昌，姚孝執(zhí)，吳伯黔.1980.云南主要構(gòu)造體系與強(qiáng)震分布規(guī)律的初步分析[J].地震研究，3（1）：47-57.

張培震，甘衛(wèi)軍，沈正康，等.2005.中國(guó)大陸現(xiàn)今構(gòu)造作用的地塊運(yùn)動(dòng)和連續(xù)變形耦合模型[J].地質(zhì)學(xué)報(bào)，79（6）：748-756.

張培震，徐錫偉，聞學(xué)澤，等.2008.2008年汶川8.0級(jí)地震發(fā)震斷裂的滑動(dòng)速率、復(fù)發(fā)周期和構(gòu)造成因[J].地球物理學(xué)報(bào)，51（4）：1066-1073.

張希，江在森，王琪，等.2003.川滇地區(qū)地殼水平運(yùn)動(dòng)特征與強(qiáng)震關(guān)系研究[J].大地測(cè)量與地球動(dòng)力學(xué)，23（3）：35-41.

鄒鎮(zhèn)宇，江在森，武艷強(qiáng)，等.2015.基于GPS速度場(chǎng)變化結(jié)果研究汶川地震前后南北地震帶地殼運(yùn)動(dòng)動(dòng)態(tài)特征[J].地球物理學(xué)報(bào)，58（5）：1597-1609.

周海濤，陳兆輝，朱爽，等.2017.川滇地區(qū)部分連續(xù)GPS測(cè)點(diǎn)形變異常分析[J].地震研究，40（1）：94-100.

Dong D，Hering T A，King R W.1998.Estimating regional deformation from a combination of space and terrestrial geodetic data[J].J Geophys Res，72（4）：200-214.

Maurin T，Masson F，Rangin C，et al.2010.First global positioning system results in northern Myanmar：Constant and localized slip rate along the Sagaing fault[J].Geology，38（7）：591-594.

Wu Y Q，Jiang Z S，Yang G H.et al.2011.Comparison of GPS strain rate computing methods and their reliability[J].Geophys J Int，185（2）：703-717，doi：10.1111/j.1365-246X.2011.04976.x.

Wu Y Q，Jiang Z S，Zhao J，et al.2015.Crustal Deformation before the 2008 Wenchuan MS8.0 Earthquake studied using GPS Data[J].J Geodyn，（85）：11-23，doi：10.1016/j.jog.2014.12.002.

Abstract Based on the GPS velocity field data of 1999—2007 and 2011—2013，we used the least squares configuration method and GPS velocity profile results to synthetically analyze the dynamic evolution characteristics of crustal deformation in Yunnan area before and after the Wenchuan earthquake. The dynamic evolution of GPS velocity field shows that the direction is gradually changed from the south of the southern part of the Sichuan-Yunnan block to the south-west of the southern Yunnan block and there is a clear relative motion characteristic near the block boundary fault zone. Compared with the GPS velocity of 1999—2007，the results of 2011—2013 also reflect segmental deformation characteristics of the block boundary fault zone. South-east movement has a significant increase. This may be related to crustal deformation adjustment after the Wenchuan earthquake. The dynamic evolution of strain parameters shows that ‘middle tension，squeeze at both ends is in the whole area and the distribution of deformation（shear，tension or extrusion）is closely related to the background motion and deformation characteristics of the main fault zone. Compared with results during the period of 1999—2007，the tensile zone in 2011—2013 has been expanded in the eastward and southward direction. The extrusion deformation of the eastern boundary（Xiaojiang fault zone）of the Sichuan-Yunnan block is no longer significant，mainly concentrated in the northern section of the Xiaojiang fault zone，which may be related to the post-seismic deformation adjustment of the Wenchuan earthquake. The GPS velocity profile results show that the left-lateral slip velocity of the Xiaojiang fault zone reduced gradually from north to south（from 10 mm/a to 5 mm/a），and the width of the northern section is wider. The right-hand slip rate of the Honghe fault zone is about 4 mm/a，and the deformation width is wider. The dynamic results show that the Wenchuan earthquake has little effect on the deformation modes of the two fault zones.

Keywords：least squares configuration；GPS velocity field；strain parameter；GPS velocity profile