王伶俐，王青華，張 勇，王 巖（云南省地震局，云南昆明 650041）

?

基于GPS的云南地區(qū)主要斷裂帶現(xiàn)今運(yùn)動(dòng)特征分析

王伶俐，王青華，張 勇，王 巖
（云南省地震局，云南昆明 650041）

摘 要：利用GAMIT/ GLOBK軟件對云南境內(nèi)2009—2013年間三期陸態(tài)網(wǎng)絡(luò)區(qū)域網(wǎng)聯(lián)測數(shù)據(jù)進(jìn)行處理，得到歐亞框架下的測站運(yùn)動(dòng)速度場。將云南地區(qū)劃分為4個(gè)活動(dòng)地塊以及7個(gè)次級構(gòu)造單元，基于塊體整體旋轉(zhuǎn)與均勻應(yīng)變模型（REHSM）研究了云南主要塊體以及斷裂帶的運(yùn)動(dòng)與變形特征，并結(jié)合GPS速度剖面，求得云南省主要斷裂帶的平均滑動(dòng)速率。結(jié)果顯示：位于川滇菱塊東側(cè)的則木河斷裂、小江斷裂、蓮峰巧家斷裂、彌勒師宗斷裂均呈左旋走滑運(yùn)動(dòng)特征，左旋速率分別為5. 4mm/ a、5. 3～7. 2mm/ a、0. 5mm/ a、1. 8mm/ a；而其西側(cè)紅河斷裂、無量山斷裂、瑞麗龍陵斷裂、南汀河斷裂、怒江斷裂、瀾滄江斷裂均以右旋走滑活動(dòng)為主，右旋速率分別為0. 2～4. 3mm/ a、2. 5mm/ a、1. 7mm/ a、1. 1mm/ a、5. 4mm/ a、4. 6mm/ a；位于塊體內(nèi)部的麗江寧蒗斷裂、易門斷裂等均呈弱活動(dòng)特征。

關(guān)鍵詞：云南；GAMIT/ GLOBK；活動(dòng)塊體；GPS速度剖面；應(yīng)變模型

0 引言

地震是地殼活動(dòng)的一種表現(xiàn)，斷層活動(dòng)速率是衡量地殼活動(dòng)程度的重要指標(biāo)之一，它反映了斷層上應(yīng)變能量釋放的速度，與地震的孕育產(chǎn)生有一定聯(lián)系，這已被一些大震前后的形變所證實(shí)［1-8］。斷裂帶滑動(dòng)速率值的確定對地殼穩(wěn)定性評估、強(qiáng)震重復(fù)間隔的研究和地震中長期預(yù)報(bào)等都是很有意義的［10-11］。云南地處歐亞、太平洋和印度三大地質(zhì)板塊碰撞匯聚地帶，其地殼變動(dòng)劇烈，活動(dòng)斷裂發(fā)育［22-23］，作為我國強(qiáng)震活動(dòng)重點(diǎn)監(jiān)視區(qū)之一，該區(qū)域的主要斷裂帶的活動(dòng)特征研究對認(rèn)識(shí)構(gòu)造變形和了解地殼構(gòu)造的動(dòng)力作用過程，捕捉與認(rèn)識(shí)孕震信息，以及活躍期內(nèi)的強(qiáng)震預(yù)測具有重要意義。盡管許多研究者給出的該區(qū)域的變形和應(yīng)變積累方式、斷層活動(dòng)性等方面的結(jié)果大體一致，但不同研究者由于研究方法以及基礎(chǔ)研究數(shù)據(jù)的不同得到的結(jié)果存在一些差異。過去利用GPS觀測得到的研究成果［16-26］主要以大尺度為主，研究對象大多是針對整個(gè)川滇菱形塊體，采用的觀測數(shù)據(jù)多是基于2007年以前的GPS數(shù)據(jù)，對云南地區(qū)主要斷裂帶近期的活動(dòng)性研究尚不充分。自2009年起，“中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)”項(xiàng)目的實(shí)施，使得云南地區(qū)擁有了點(diǎn)位分布比較密集且相對均勻的水平形變監(jiān)測網(wǎng)，并陸續(xù)開展了監(jiān)測工作。2013年10月，最新一期的陸態(tài)網(wǎng)絡(luò)區(qū)域網(wǎng)聯(lián)測工作完成。本文對2009—2013年間三期GPS復(fù)測資料進(jìn)行數(shù)據(jù)處理，建立云南地區(qū)現(xiàn)今地殼水平運(yùn)動(dòng)的速度場。綜合現(xiàn)有的云南地區(qū)現(xiàn)今水平運(yùn)動(dòng)的分區(qū)論述，進(jìn)行合理塊體劃分，以GPS速度場為約束，通過塊體變形模型和實(shí)測剖面定量分析云南主要斷裂帶運(yùn)動(dòng)特征。

1 GPS觀測與數(shù)據(jù)處理

“十一五”期間，國家重大科學(xué)工程項(xiàng)目“中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)工程”（簡稱陸態(tài)網(wǎng)絡(luò)）啟動(dòng)，該項(xiàng)目由中國地震局牽頭，在全國范圍內(nèi)進(jìn)行高密度、大范圍的GNSS觀測。該項(xiàng)目在云南省范圍內(nèi)共建成27個(gè)基準(zhǔn)站和170多個(gè)區(qū)域站，覆蓋了云南的大部地區(qū)且跨越了紅河、小江等主要活動(dòng)斷裂帶。2009年起陸續(xù)開展了陸態(tài)網(wǎng)絡(luò)區(qū)域網(wǎng)聯(lián)測工作，到目前為止，共有復(fù)測資料三期，資料采集的時(shí)間分別為2009年3—10月、2011年4—10月和2013年4—10月。每期每個(gè)測點(diǎn)的觀測時(shí)間長度均不少于90h。

由于涉及多期觀測數(shù)據(jù)，為避免由于處理方法和約束條件不同而產(chǎn)生的速度結(jié)果差別，筆者用MIT（Massachusetts Institute of Technology）和SIO（Scripps Institution of Oceanography）研制的GAMIT/ GLOBK軟件，采用統(tǒng)一的數(shù)據(jù)處理方案對所有數(shù)據(jù)進(jìn)行了處理，GAMIT數(shù)據(jù)預(yù)處理時(shí)，通過IGS的全球跟蹤站作“橋梁”，將區(qū)域站點(diǎn)納入全球GPS網(wǎng)絡(luò)中，從而獲得全球統(tǒng)一參考框架中的坐標(biāo)。具體做法是選取中國及周邊部分IGS跟蹤站（TSKB、USUD、KUNM、IISC、HYDE、BJFS、IRKT、PIMO、GUAO、WUHN、TWTF、LHAZ、URUM、SHAO、KIT3）與云南省內(nèi)的測站聯(lián)合解算，獲得測站和衛(wèi)星軌道的單日區(qū)域松弛解；GLOBK數(shù)據(jù)后處理時(shí)，利用GLOBK將GAMIT處理得到的單日松弛解和SOPAC（Scripps Orbital and Permanent Array Center）網(wǎng)站下載的全球IGS站的單日松弛解合并，并對全部單日解進(jìn)行整網(wǎng)平差；為了更好反映研究區(qū)域內(nèi)部的構(gòu)造變形，筆者選擇KOSG、ONSA等位于歐亞板塊的11個(gè)速度穩(wěn)定的IGS站，構(gòu)建穩(wěn)定的歐亞板塊，以ITRF2008提供的站坐標(biāo)和速度作為先驗(yàn)值，以其2倍誤差作為約束，求出各測站的平均速度，最后將求得的速度扣除歐亞板塊旋轉(zhuǎn)后的值作為云南地區(qū)地殼運(yùn)動(dòng)的速度場（圖1）。

2 基于塊體變形模型研究斷裂帶的運(yùn)動(dòng)特征

從大的尺度來講，邊界帶的變形依賴于塊體的差異運(yùn)動(dòng)。近十幾年來，不少學(xué)者利用塊體變形模型對斷裂帶的運(yùn)動(dòng)特征進(jìn)行研究，已經(jīng)取得了許多重要的成果［12-15］。

根據(jù)李延興等塊體整體旋轉(zhuǎn)與均勻應(yīng)變模型（REHSM）［12-15］：

其中，υe和υn分別代表測站的東方向和北方向速率；λ和φ分別代表測站的經(jīng)度和緯度；r表示地球半徑；λ0和φ0分別表示塊體幾何中心的緯度和經(jīng)度；ωx、ωy和ωz分別代表塊體歐拉運(yùn)動(dòng)矢量的三分量；εe、εn和εen分別代表塊體東方向和北方向的正應(yīng)變率以及剪應(yīng)變率。

參照公式（2）可進(jìn)一步計(jì)算最大主應(yīng)變率ε1和ε2及其方位角A、最大剪應(yīng)變率γmax、面應(yīng)變率Δ等參數(shù)。

筆者參照闞榮舉等觀點(diǎn)將云南地區(qū)劃分為滇中、滇東（Dd）、印支（Ls）以及滇緬泰4個(gè)一級構(gòu)造單元，并在上述一級塊體劃分的基礎(chǔ)上將云南劃分騰沖塊體（Tc）、保山塊體（Ba）、蘭坪—思茅弧后盆地、鹽源—麗江陸緣坳陷（Yi）、滇中坳陷（Dz）、康滇古隆起（Kd）、滇東坳褶帶等7個(gè)二級構(gòu)造單元［21］。受觀測點(diǎn)密度以及分布限制，在上述活動(dòng)地塊中選擇滇東、印支等7個(gè)構(gòu)造單元，同時(shí)為了更詳盡地反應(yīng)斷裂帶運(yùn)動(dòng)特征，筆者在此基礎(chǔ)上增加了川滇塊體南段（Cd）和華南地塊（Hu）共計(jì)9個(gè)活動(dòng)地塊為研究對象。采用REHSM模型，以計(jì)算得到的云南區(qū)域測站基于歐亞框架的GPS速度場為約束，采用最小二乘法，計(jì)算得到9個(gè)活動(dòng)塊體的運(yùn)動(dòng)參數(shù)和應(yīng)變參數(shù)（表1）。再由公式（2）得到每個(gè)塊體的最大主張應(yīng)變率、最小主壓應(yīng)變率和主壓應(yīng)變率方位角（圖2）。

表1　活動(dòng)塊體模型參數(shù)計(jì)算結(jié)果Tab. 1 The calculation results of block model parameters

根據(jù)計(jì)算得到各活動(dòng)塊體的運(yùn)動(dòng)參數(shù)和應(yīng)變參數(shù)（表1），在分隔兩塊體的斷裂帶上選取部分點(diǎn)按照公式（1）計(jì)算模擬速度，將計(jì)算得到的兩個(gè)速度矢量相減再投影到斷裂上，即可得到平行斷裂（走滑分量）和垂直斷裂（擠壓/拉張分量）的速度值。計(jì)算時(shí)設(shè)置斷層走向時(shí)規(guī)定斷層起點(diǎn)緯度小于終點(diǎn)緯度，對于投影到斷層的走滑速率值與斷層走向相反為負(fù)值（即為右旋）、與斷層走向相同為正值（即為左旋），對于投影到斷層的垂直方向運(yùn)動(dòng)速率以斷層走向逆時(shí)針旋轉(zhuǎn)90°方向?yàn)檎较颍礊槔瓘垼?。各主要斷裂的平均滑?dòng)速率見表2。

3 利用GPS實(shí)測剖面分析斷裂帶的運(yùn)動(dòng)特征

利用上述塊體模型計(jì)算斷層滑動(dòng)速率時(shí)，得到的是各塊體的整體變形參數(shù)，計(jì)算模型也是建立在兩塊體邊界為非連續(xù)的前提條件下，用于描述斷裂帶兩側(cè)連續(xù)變形時(shí)必然產(chǎn)生模型偏差，因此，筆者在用塊體模型計(jì)算了主要斷裂帶的滑動(dòng)速率的同時(shí)，還對斷裂帶進(jìn)行GPS速度剖面分析，便于分析斷裂帶遠(yuǎn)場及近場的連續(xù)運(yùn)動(dòng)信息。受GPS臺(tái)站分布的限制，文中選取了小江斷裂、紅河斷裂、無量山斷裂和怒江斷裂等15條主要斷裂帶進(jìn)行剖面分析，由于瀾滄江、怒江斷裂北段臺(tái)站數(shù)量稀疏且均集中在2條斷裂附近，缺少中、遠(yuǎn)場數(shù)據(jù)的約束，因此本文選取他們中南段做剖面分析，并根據(jù)地質(zhì)學(xué)的研究成果，將紅河斷裂、小江斷裂帶分成北段、中段、南段來研究。圖2中灰色線段為實(shí)際斷裂分布，采用的是中國地震局GIS的地震預(yù)報(bào)分析系統(tǒng)MapSIS軟件中文件“中國斷層. Tab”為基礎(chǔ)繪制而成，實(shí)際斷裂構(gòu)造并不完全筆直，只能用多段直線段近似描述，各主要斷裂的分段模型以及剖面分布如圖3所示。

表2　基于塊體模型得到主要斷裂帶活動(dòng)速率Tab. 2 Active rates of principal faults based on REHSM model

具體計(jì)算方法：分別以斷裂模型兩盤的點(diǎn)作為參考，將第一節(jié)中計(jì)算得到的云南地區(qū)地殼運(yùn)動(dòng)的GPS速度場根據(jù)斷層剖面分布分別投影到平行及垂直于斷層走向的軸線上，計(jì)算時(shí)先設(shè)置斷層走向：規(guī)定斷層起點(diǎn)緯度小于終點(diǎn)緯度，投影到斷層的走滑速率值與斷層走向相反為負(fù)值、與斷層走向相同為正值；對于投影到斷層的垂直方向運(yùn)動(dòng)速率以斷層走向順時(shí)針旋轉(zhuǎn)90°方向?yàn)檎较?。圖4中先統(tǒng)一設(shè)定斷層左盤為負(fù)、右盤為正，垂線表示斷層的位置；將斷層剖面兩側(cè)的站點(diǎn)展繪到圖上，橫坐標(biāo)表示站點(diǎn)相對斷層垂向距離，縱坐標(biāo)表示投影后的速度分量。為了便于分析，筆者對同一斷層剖面兩側(cè)的測點(diǎn)縱坐標(biāo)進(jìn)行了平移，將所有斷層剖面展布到同坐標(biāo)系統(tǒng)中，圖4 （a）中斜率反映的是斷裂帶的走滑特征，斜率為正為左旋走滑斷裂，反之為右旋走滑斷裂；圖4 （b）中斜率反映的是斷裂的拉張（擠壓）特征，正表示拉張，反之則為擠壓。將各條斷裂兩盤的投影速率分量分別求平均值，兩盤差異運(yùn)動(dòng)速率（右盤減去左盤）即為各段的錯(cuò)動(dòng)速率（表3），斷裂錯(cuò)動(dòng)速率以左旋走滑和拉張為正。

4 討論與結(jié)論

從大的動(dòng)力學(xué)背景看，云南地區(qū)主要受到來自3個(gè)方面力的作用：一是印度板塊與歐亞板塊碰撞，使西藏地塊東移，產(chǎn)生四川及川滇菱塊的SE、SSE向作用力；二是來自華南地塊的NNW、NW向應(yīng)力的作用；三是印度板塊向東經(jīng)緬甸對云南地區(qū)的側(cè)向擠壓力［19］。從上述塊體運(yùn)動(dòng)及應(yīng)變特征來看，在這3個(gè)方面力的作用下，云南主要地塊運(yùn)動(dòng)方向逐漸由SSE向SSW變化，有順時(shí)針旋轉(zhuǎn)的特征，運(yùn)動(dòng)幅度由西向東、由北向南逐漸減弱，菱形塊體外各個(gè)塊體運(yùn)動(dòng)速度大幅衰減。從應(yīng)變率參數(shù)結(jié)果來看，華南地塊、滇東塊體的主要受SE向壓應(yīng)力場控制，到滇中地塊轉(zhuǎn)為SE—SSE方向，滇西北地區(qū)應(yīng)力場方向?yàn)镾SE方向，滇西南的印支地塊為NNW向，滇西南的騰沖—保山地塊主要受NE—NNE向應(yīng)力場控制。正是這幾方面作用力控制了云南及其鄰區(qū)的主要斷層活動(dòng)及強(qiáng)震活動(dòng)。

作為川滇塊體東邊界的則木河、小江斷裂是研究區(qū)內(nèi)活動(dòng)性最強(qiáng)的大型走滑斷裂帶，本文塊體應(yīng)變反演結(jié)果得到則木河平均走滑速率分別為5. 19mm/ a、5. 29mm/ a，呈左旋拉張?zhí)卣鳌t木河滑動(dòng)速率結(jié)果與任金衛(wèi)等根據(jù)地貌的錯(cuò)斷情況得到4. 7mm/ a滑動(dòng)速率值相近。小江斷裂帶又分為東西兩支，由于GPS臺(tái)站分布密度的限制，本文無法獲得它們各自的錯(cuò)動(dòng)速率，只能給出總的活動(dòng)速率，與Wen［28］根據(jù)歷史地震及大地震同震位錯(cuò)的貢獻(xiàn)得到的9. 6mm/ a滑動(dòng)速率以及He等［9］得到的13. 0～16. 5mm/ a整體滑動(dòng)速率結(jié)果差別較大。本文剖面結(jié)果與塊體模型算得的值基本一致，顯示小江斷裂北、中、南三段滑動(dòng)速率略有不同，分別為5.97mm/ a、7.19mm/ a、5.30mm/ a。與Shen［26］等采用GPS剖面所得結(jié)果（7±2）mm/ a較一致。小江斷裂帶北段和中段的滑動(dòng)速率明顯大于南段，筆者推測導(dǎo)致這種變化特征的原因可能是由于位于小江斷裂南端的曲江—石屏斷裂對塊體運(yùn)動(dòng)有一定吸收作用。

表3　基于剖面分析得到主要斷裂活動(dòng)速率Tab. 3 Active rates of major faults based on profiles analysis method

紅河斷裂是具有深部背景的大型活動(dòng)斷裂，與金沙江斷裂一同被認(rèn)為是川滇菱形塊體的西南邊界，從該塊體應(yīng)變結(jié)果來看，紅河斷裂表現(xiàn)出弱右旋走滑特征。從剖面分析的結(jié)果來看，北段（洱源縣福壽場）和中段（彌渡縣至新平縣北大斗門）右旋滑動(dòng)速率分別為4. 33mm/ a、3. 84mm/ a，中南段（新平水塘至河口縣）現(xiàn)今活動(dòng)不強(qiáng)，與申重陽等人得到的紅河斷裂帶的活動(dòng)劇烈程度以北為最、中段次之、南段最小結(jié)論一致［22-23］。

印度板塊向東經(jīng)緬甸對云南地區(qū)的側(cè)向擠壓力，使得云南西部地區(qū)的原本呈南北展布的大斷裂在小灣附近呈弧形展布特征［19］。該區(qū)域三條大型弧形斷裂，包括怒江斷裂、瀾滄江斷裂、無量山斷裂與鮮水河—小江斷裂帶相比，活動(dòng)性明顯降低，塊體模型結(jié)果顯示怒江斷裂、瀾滄江斷裂均呈右旋走滑特征。剖面分析因?yàn)槭茳c(diǎn)位密度以及分布的限制，僅對這幾條斷裂的中南段進(jìn)行了研究。剖面結(jié)果顯示瀾滄江永平至小灣段右旋走滑速率為4. 62mm/ a，拉張速率為0. 06mm/ a，與王閻昭等得到的碧土以南段右旋走滑速率為（2. 4 ±1. 2）mm/ a［27］結(jié)論相近；怒江斷裂帶北起貢山，經(jīng)福貢到瀘水后，分為東西兩支。塊體應(yīng)變結(jié)果顯示怒江—騰沖，龍陵—瑞麗斷裂帶的整體右旋走滑特征及拉張?zhí)卣?；剖面結(jié)果反應(yīng)了瀘水到龍陵段東西兩支的總體右旋走滑速率為5. 42mm/ a，與Shen［26］等所得結(jié)果6mm/ a一致；無量山斷裂為云南普洱地區(qū)的重要斷裂帶，總體走向NW向，呈反“S”形，具有2. 52mm/ a的右旋走滑活動(dòng)。由于斷裂帶南段由若干條走向接近的次級斷裂組成，上述結(jié)果反映的是斷裂帶總體滑動(dòng)速率，與王閻昭（4. 3±1. 1）mm/ a［27］的右旋走滑結(jié)果相近。同樣位于川滇塊體以西NW向展布的瑞麗龍陵斷裂，南汀河斷裂呈右旋走滑兼張性特征，張性明顯，拉張分量分別為2. 83mm/ a、3. 07mm/ a。

川滇塊體以東的蓮峰巧家斷裂，彌勒師宗斷裂位于穩(wěn)定的華南應(yīng)力區(qū)，從塊體應(yīng)變分析的結(jié)果來看，活動(dòng)性均不太強(qiáng)，剖面定量分析結(jié)果顯示左旋走滑速率分別為0. 49mm/ a、1. 84mm/ a，呈壓性。

川滇塊體內(nèi)部的麗江—小金河斷裂作為一條次級構(gòu)造邊界帶，本文塊體應(yīng)變反演結(jié)果給出該斷裂的中段、南西段平均走滑速率為1. 82mm/ a，擠壓速率為2. 26mm/ a。與Shen［26］等利用GPS數(shù)據(jù)得到麗江—小金河斷裂3mm/ a的左旋走滑速率相當(dāng)，但是剖面結(jié)果顯示南西段麗江寧蒗斷裂呈弱活動(dòng)特征；位于菱塊內(nèi)部的易門斷裂從塊體應(yīng)變結(jié)果來看，活動(dòng)性不太明顯。

總體來看：紅河斷裂作為川滇菱形塊體的西南邊界，在川滇菱形塊體SE-SSE方向和滇西南的印支地塊NNW方向應(yīng)力的共同作用下呈右旋運(yùn)動(dòng)特征；則木河、小江斷裂帶作為川滇菱形塊體東邊界，川滇菱形塊體以及滇東地塊的NW方向共同作用下呈左旋運(yùn)動(dòng)，位于川滇菱塊的東側(cè)蓮峰巧家斷裂、彌勒師宗斷裂均呈左旋走滑運(yùn)動(dòng)特征，而其西側(cè)無量山斷裂、瑞麗龍陵斷裂、南汀河斷裂、怒江斷裂和瀾滄江斷裂均已以右旋走滑活動(dòng)為主；位于塊體內(nèi)部的易門斷裂、麗江—寧蒗斷裂等均呈弱活動(dòng)特征。

本文所得斷裂滑動(dòng)速率結(jié)果與前人所得結(jié)果相比略微偏低，尤其是川滇菱塊邊界的則木河、小江斷裂帶。可能是由于兩方面的原因：一是本文采用的GPS觀測資料是2009—2013年陸態(tài)網(wǎng)絡(luò)區(qū)域網(wǎng)聯(lián)測資料，前人采用基本上是2009年以前觀測資料，由于時(shí)間跨度以及GPS點(diǎn)位分布密度不同，加之采用分析方法及模型的差異，導(dǎo)致結(jié)果略微不同；二是2008年發(fā)生汶川8. 0級地震，位于川滇菱塊邊界斷裂在川滇地區(qū)地殼運(yùn)動(dòng)格局中的作用不容忽視，因此推測這些差異或許與汶川地震引起川滇地區(qū)大范圍應(yīng)力調(diào)整有關(guān)。

參考文獻(xiàn)

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

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

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

［4］ 杜平山.則木河斷裂帶的走滑位移及滑動(dòng)速率［J］.四川地震，2000，1（2）：49-64.

［5］ 顧國華，牛紅葉，孟國杰，等.中國地殼運(yùn)動(dòng)觀測網(wǎng)絡(luò)數(shù)據(jù)處理與分析［J］.地震，2000，20（增刊）：203-210.

［6］ 郭良遷.中國大陸地殼運(yùn)動(dòng)與汶川Ms8. 0級地震孕育的關(guān)系［J］.地球物理學(xué)報(bào)，2009，52（2）：531-537.

［7］ 郭良遷，李延興.中蒙地塊的水平運(yùn)動(dòng)及應(yīng)變特征［J］.大地測量與地球動(dòng)力學(xué)，2005，25（4）：17-22.

［8］ 萬永革，沈正康，盛書中，等. 2008年汶川大地震對周圍斷層的影響［J］.地震學(xué)報(bào)，2009，31（2）：128-139.

［9］ He H L，Yasutaka I，Song F M，et al. Late quaternary slip rate of the Xiaojiang fault and its implication ［J］. Seismol，Geol，2002，24（1）：14-26.

［10］ 宋方敏，汪一鵬，俞維賢，等.小江活動(dòng)斷裂帶［M］.北京：地震出版社，1998.

［11］ 馬宗晉.地震與斷裂關(guān)系的討論，中國活動(dòng)斷裂［M］.北京：地震出版社，1982.

［12］ 李延興，黃珹.板內(nèi)塊體的剛性彈塑性運(yùn)動(dòng)模型與中國大陸主要塊體的應(yīng)變狀態(tài)［J］.地震學(xué)報(bào)，2001，23（6）：565-572.

［13］ 李延興，李智，張靜華，等.中國大陸及周邊地區(qū)的水平應(yīng)變場［J］.地球物理學(xué)報(bào)，2004，7（2）：222-231.

［14］ 李延興，楊國華，李智，等.中國大陸活動(dòng)地塊的運(yùn)動(dòng)與應(yīng)變狀態(tài)［J］.中國科學(xué)D輯，2003，33（增刊）：65-81.

［15］ 李延興，張靜華.汶川Ms8. 0地震孕育發(fā)生的機(jī)制與動(dòng)力學(xué)問題［J］.地球物理學(xué)報(bào)，2009，52（2）：519-529.

［16］ 駱佳驥，崔效鋒，胡幸平，等.川滇地區(qū)活動(dòng)塊體劃分與現(xiàn)代構(gòu)造應(yīng)力場分區(qū)研究綜述［J］.地震研究，2012，35（3）：309-317.

［17］ 呂江寧，沈正康，王敏.川滇地區(qū)現(xiàn)代地殼運(yùn)動(dòng)速度場和活動(dòng)塊體模型研究［J］.地震地質(zhì)，2003，25（4）：543-554.

［18］ 劉宇平，陳智粱.青藏高原東部及周邊現(xiàn)時(shí)地殼運(yùn)動(dòng)［J］.沉積與特提斯地質(zhì)，2003，23（4）：1-7. ［19］ 錢曉東.云南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場研究［J］.地震地質(zhì)，2011，33（1）：91-105.

［20］ 喬學(xué)軍，王琪，杜瑞林.川滇地區(qū)活動(dòng)地塊現(xiàn)今地殼形變特征［J］.地球物理學(xué)報(bào)，2004，47（5）：805-811.

［21］ 皇甫崗，陳顒，秦嘉政，等.云南地震活動(dòng)性［M］.云南：云南科技出版社，2010.

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

［23］ 申重陽，吳云，王琪，等.云南地區(qū)主要斷層運(yùn)動(dòng)模型的GPS數(shù)據(jù)反演［J］.大地測量與地球動(dòng)力學(xué)，2002，22（3）：46-51.

［24］ 宋方敏，李如成，徐錫偉.四川大涼山斷裂帶古地震研究初步結(jié)果［J］.地震地質(zhì)，2002，24（1）：27-34.

［25］ 蘇有錦，秦嘉政.川滇地區(qū)強(qiáng)地震活動(dòng)與區(qū)域新構(gòu)造運(yùn)動(dòng)的關(guān)系［J］.中國地震，2001，17（1）：24-34.

［26］ Shen Z K，Lü J，Wang M，et al. Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau［J］. Geophys Res，2005，110：B11409. doi：10. 1029/2004JB003421.

［27］ 王閻昭，王恩寧，沈正康，等.基于GPS資料約束反演川滇地區(qū)主要斷裂現(xiàn)今活動(dòng)速率［J］.中國科學(xué)（D輯），2008，38（5）：582-597.

［28］ Wen X Z. The uniform-slip method for estimating mean slip-rate of strike-slip fault［J］. EarthqPred Res，1998，7：170-182.

Analysis of Current Activity of Main Faults in Yunnan Region based on GPS

Wang Lingli，Wang Qinghua，Zhang Yong，Wang Yan
（Earthquake Administration of Yunnan Province，Kunming 650041，China）

Abstract：GPS data from CMONOC in the period of 2009—2013 in Yunnan region were processed using GAMIT/ GLOBK software，and a velocity field relative to the Eurasian reference frame was obtained. Based on geological division of blocks，Yunnan region can divide into seven active secondary crustal blocks. We studied the movement velocity and deformation characteristics of these blocks and principal faults in Yunnan region by REHSM model and GPS velocity profiles analysis method. The result shows that the Zemuhe，Xiaojiang，Lianfeng-Qiaojia，Mile-Shizong faults on the east side of the Sichuan-Yunnan rhombic block show left-lateral characteristics，with a rate of 5. 4mm/ a，5. 3- 7. 2mm/ a，0. 5mm/ a，1. 8mm/ a respectively. The Red River，Wuliangshan，Ruili-Longling，Nantinghe，Nujiang，Lancangjiang faults on the west side of the Sichuan-Yunnan rhombic block show right-lateral characteristics，with a rate of 0. 2- 4. 3mm/ a，2. 5mm/ a，1. 7mm/ a，1. 1mm/ a，5. 4mm/ a，4. 6mm/ a respectively. Lijiang-Ninglang and Yimen faults that locate inside the block show weak activities.

Keywords：Yunnan；GAMIT/ GLOBK；active block；GPS velocity profiles；strain model

中圖分類號：P315. 2

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號：1673-8047（2016）01-0001-08

收稿日期：2016-01-06

基金項(xiàng)目：云南省地震局科技人員傳幫帶項(xiàng)目（C2- 2014007）；中國地震局2015年度震情跟蹤定向工作任務(wù)（2015010219）

作者簡介：王伶俐（1981—），女，碩士，高級工程師，主要從事GPS數(shù)據(jù)處理及分析、地殼形變觀測、地震活動(dòng)性分析工作。