滇中主要活斷層現(xiàn)今活動(dòng)性研究

李長(zhǎng)軍 甘衛(wèi)軍 秦姍蘭 郝明 宋尚武

摘要：為研究滇中地區(qū)主要斷層的活動(dòng)特征，利用1999—2007年和2011—2017年2期GPS觀測(cè)資料以及地質(zhì)資料，基于Okada位錯(cuò)模型反演了研究區(qū)域主要斷層的滑動(dòng)速率和閉鎖深度。結(jié)果表明：（1）紅河斷裂帶的走滑速率為（1.5±1.6）～（4.7±1.5）mm/a，傾滑速率為（-3.6±1.6）～（1.9±2.4）mm/a，斷裂南段的活動(dòng)性更強(qiáng)；（2）無(wú)量山斷裂和南華—楚雄—建水?dāng)嗔训淖呋俾蕿椋?.0±1.6）～（5.6±1.5）mm/a 和（4.8±1.4）～（6.6±1.6）mm/a，傾滑速率分別為（-0.7±1.5）～（0.2±1.4）mm/a 和（-5.8±1.5）～（1.7±1.8）mm/a；（3）紅河斷裂帶元江—元陽(yáng)段和洱源—彌渡段、無(wú)量山斷裂帶和南華—楚雄—建水?dāng)嗔褞鞫翁幱谡痖g閉鎖狀態(tài)，閉鎖深度分別為6.8 km，7 km和7.2 km。

關(guān)鍵詞：斷層滑動(dòng)速率；閉鎖深度；GPS；斷層活動(dòng)性

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

0 引言

滇中地區(qū)發(fā)育有3條主要活動(dòng)斷裂帶，分別為紅河斷裂帶、無(wú)量山斷裂帶和南華—楚雄—建水?dāng)嗔褞А?duì)這3條斷裂帶之間的相互關(guān)系、運(yùn)動(dòng)學(xué)特征、變形機(jī)制、應(yīng)變分配特征及其地震危險(xiǎn)性的研究，對(duì)于認(rèn)識(shí)青藏高原隆起及周邊地殼和塊體間變形的運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)特征，具有很重要的理論和實(shí)際意義。

近年來(lái)，大地測(cè)量等地球物理技術(shù)與方法的發(fā)展和完善，為進(jìn)一步認(rèn)識(shí)紅河斷裂帶及其周緣區(qū)域的變形特征提供了新的契機(jī)?；诖蟮販y(cè)量資料的研究結(jié)果表明，紅河斷裂帶現(xiàn)今處于弱滑動(dòng)或強(qiáng)閉鎖狀態(tài)，其活動(dòng)速率約為右旋（0.4±1.6）mm/a，正斷分量為（1.5±1.3）mm/a（王閻昭等，2008；Leloup et al，1993，1995）。劉耀輝等（2015）基于塊體模型（McCaffrey et al，2005，2007，2009）得到紅河斷裂帶的右旋走滑速率，北、中、南段分別為（5.9±1.2）mm/a，（4.8±0.6）mm/a和（4.3±0.4）mm/a，且分別存在3.3 mm/a，2.6 mm/a和2.3 mm的速率虧損。王閻昭等（2008）的研究結(jié)果受觀測(cè)資料分辨率所限，很難給出比較確切的結(jié)果；劉耀輝等（2015）受模型的限制，得到的反演結(jié)果與地質(zhì)上2 mm/a的滑動(dòng)速率相差較大。

為了分析探討紅河斷裂帶、無(wú)量山斷裂帶和南華—楚雄—建水?dāng)嗔褞КF(xiàn)今運(yùn)動(dòng)學(xué)特征和應(yīng)變分配狀態(tài)，本文以紅河斷裂帶及其周邊斷裂附近1999—2007年和2011—2017年2期GPS速度場(chǎng)為基礎(chǔ)資料，結(jié)合斷層位錯(cuò)模型（Okada，1985；Gan et al，2000），反演了上述3條主要活動(dòng)斷裂帶現(xiàn)今滑動(dòng)速率和鎖定深度。

1 研究區(qū)域概況

紅河斷裂帶位于青藏高原橫向擠出變形的前緣，是連接青藏高原東南緣和南海擴(kuò)張的一級(jí)構(gòu)造。該斷裂帶在陸-陸碰撞作用的形成和演化過程中，具有舉足輕重的作用。紅河斷裂帶的演化可分為3個(gè)階段：上新世之前以左旋走滑運(yùn)動(dòng)為主，上新世活動(dòng)方式開始從左旋向右旋轉(zhuǎn)變，第四紀(jì)以來(lái)主要是右旋滑動(dòng)（Allen，Han，1984；Zhu et al，2009；Tapponnier et al，1990；Leloup et al，1993，1995；Schoenbohm et al，2006；常祖峰等，2015）。

地球物理研究結(jié)果表明，紅河斷裂兩邊存在顯著的結(jié)構(gòu)差異，斷裂帶下插到地幔以下（Wang et al，2014），是一條重要的深大邊界斷裂，此外，斷裂帶底部存在低速異常，表明介質(zhì)屬性較弱（Sun et al，2014）。該斷裂帶第四紀(jì)以來(lái)的活動(dòng)速率約為2 mm/a（Wang，1998；虢順民等，2013），斷裂NW段包含小部分正斷分量，而在東段有較小的逆沖分量，總體表現(xiàn)為順時(shí)針旋轉(zhuǎn)和緩慢滑動(dòng)。古地震研究結(jié)果表明，紅河斷裂帶自晚第四紀(jì)以來(lái)存在活動(dòng)，并伴隨有古地震發(fā)生，其北段發(fā)生了包括1652年彌渡7級(jí)地震和1925年大理7級(jí)地震在內(nèi)的9次6.0級(jí)以上地震；南段存在全新世中、晚期活動(dòng)的證據(jù)（李西等，2016）；中段500年以內(nèi)缺少?gòu)?qiáng)震，形成了紅河地震空區(qū)?，F(xiàn)代地震活動(dòng)性研究表明，自1970年以來(lái)紅河斷裂帶中部小震活動(dòng)性較弱（Wen et al，2016），處于低b值高應(yīng)力區(qū)域（謝卓娟等，2015）。

在紅河斷裂帶的東西兩側(cè)，發(fā)育2條第四紀(jì)以來(lái)比較活躍的斷裂帶，即無(wú)量山斷裂帶和南華—楚雄—建水?dāng)嗔褞?。無(wú)量山斷裂帶位于紅河斷裂帶以西的蘭坪—思茅盆地內(nèi)，斷裂帶由多條斷裂組成，由東至西主要包括磨黑斷裂、寧洱斷裂、普文斷裂和景谷—云仙斷裂。該斷裂帶表現(xiàn)為明顯的右旋走滑性質(zhì)，歷史上發(fā)生了10次6級(jí)以上強(qiáng)震及20次5～5.9級(jí)地震，構(gòu)成了一條NW向地震活動(dòng)帶（虢順民等，1999；常祖峰等，2016）。其中，磨黑斷裂的滑動(dòng)速率為1.6～2.5 mm/a，寧洱斷裂的滑動(dòng)速率為2.5～2.8 mm/a，普文斷裂的滑動(dòng)速率為1.0～2.0 mm/a，景谷—云仙斷裂的滑動(dòng)速率為1.3～1.4 mm/a。

南華—楚雄—建水?dāng)嗔褞怯沙邸先A斷裂、曲江斷裂和石屏—建水?dāng)嗔褬?gòu)成的一組斜列的右旋走滑為主的活動(dòng)斷裂帶。其中，楚雄—南華斷裂的滑動(dòng)速率為1.6～2.0 mm/a（常祖峰等，2015），石屏—建水?dāng)嗔训钠骄交瑒?dòng)速率為3.0～3.6 mm/a（韓新民等，1982）或2.8～5.3 mm/a（王怡然等，2015），曲江斷裂晚第四紀(jì)以來(lái)的水平滑動(dòng)速率為2.3～4.0 mm/a（王洋等，2015；聞學(xué)澤等，2011）。目前，紅河斷裂帶主斷裂南段地震活動(dòng)性較低，而地震活動(dòng)主要集中在石屏—建水?dāng)嗔押屯ê！瓟嗔堰@一區(qū)域附近（圖1），比如1970年MS7.7通海地震。

2 數(shù)據(jù)和方法

2.1 GPS數(shù)據(jù)及其處理方法

國(guó)家重大科技基礎(chǔ)設(shè)施建設(shè)項(xiàng)目“中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)”（簡(jiǎn)稱“網(wǎng)絡(luò)工程”）和“中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)”（簡(jiǎn)稱“陸態(tài)網(wǎng)絡(luò)”）的實(shí)施，為本文提供了較密集的GPS觀測(cè)資料。由于“陸態(tài)網(wǎng)絡(luò)”項(xiàng)目的觀測(cè)時(shí)間從2009年開始，同時(shí)為了減小汶川地震的影響，將觀測(cè)資料分為2段，第一段為“網(wǎng)絡(luò)工程”項(xiàng)目的40個(gè)點(diǎn)（圖1），觀測(cè)時(shí)間為1999—2007年；第二段為“網(wǎng)絡(luò)工程”和“陸態(tài)網(wǎng)絡(luò)”2個(gè)項(xiàng)目總計(jì)50個(gè)點(diǎn)（圖1），觀測(cè)時(shí)間為2011—2017年。此外，在進(jìn)行數(shù)值模擬時(shí)，剔除了觀測(cè)噪聲較大的點(diǎn)。數(shù)據(jù)處理采用GAMIT/GLOBK軟件（Herring et al，2010），首先利用GAMIT獲得GPS測(cè)站坐標(biāo)、衛(wèi)星軌道及臺(tái)站對(duì)流層天頂延遲等參數(shù)及其方差-協(xié)方差矩陣的單日松弛解。為避免因框架和模型的差異所引起的定位結(jié)果差異混同于真實(shí)的形變信息，同期的全球IGS站及網(wǎng)基準(zhǔn)站數(shù)據(jù)，采用相同的模型和方法進(jìn)行統(tǒng)一的處理。然后利用GLOBK軟件（卡爾曼濾波）進(jìn)行多時(shí)段綜合解算，以獲得網(wǎng)平差結(jié)果。在進(jìn)行GLOBK綜合解算時(shí)，把區(qū)域站單日松弛解和全球IGS站單日松弛解捆綁在一起，并選取全球均勻分布的90多個(gè)穩(wěn)定站（GAMIT/GLOBK軟件提供）為基準(zhǔn)點(diǎn)，將區(qū)域站速度場(chǎng)固定到ITRF2008框架下。最后利用Altamimi等（2011）提供的 ITRF2008框架下歐亞板塊歐拉極，獲取研究區(qū)域相對(duì)穩(wěn)定歐亞板塊的水平運(yùn)動(dòng)速度場(chǎng)（圖4a，圖5a中的紅色箭頭）。

GAMIT軟件處理時(shí)采用的處理策略和模型如下：采用軌道的松弛模式，在估算測(cè)站位置的同時(shí)，允許衛(wèi)星軌道（IGS精密星歷）和地球自轉(zhuǎn)參數(shù)（Bull_A）有微量的調(diào)整；同時(shí)參數(shù)估計(jì)衛(wèi)星天線的Z-偏差。數(shù)據(jù)處理中涉及到的地球重力場(chǎng)、固體潮和極潮分別采用EGM08模型、IERS03固體潮模型和IERS10規(guī)范的平均地極模型。海潮引發(fā)的測(cè)站地殼形變改正采用最新的全球海潮模型FES2004（Lyard et al，2006），此外，還顧及了海潮導(dǎo)致的地球質(zhì)心變化。對(duì)流層干分量的天頂延遲由GPT模型（Boehm et al，2007）計(jì)算獲得，同時(shí)，每個(gè)測(cè)站每1小時(shí)估計(jì)1個(gè)參數(shù)以改正濕分量導(dǎo)致的天頂延遲，映射函數(shù)采用VMF1（Boehm et al，2006）?？紤]了大氣的不均勻性，對(duì)每個(gè)測(cè)站的東西向和南北向各附加1個(gè)大氣水平梯度參數(shù)。數(shù)據(jù)誤差為衛(wèi)星截止高度角的函數(shù)，函數(shù)系數(shù)由數(shù)據(jù)驗(yàn)后殘差擬合確定。

2.2 小震活動(dòng)性資料

為準(zhǔn)確給定各斷層的鎖定深度，在P、S波成像確定上述主要活斷層深部結(jié)構(gòu)的基礎(chǔ)上（Chen et al，2013；Lei et al，2014；Sun et al，2014；Hu et al，2018），筆者還參考了小震精定位的結(jié)果（Fang et al，2015），如圖2所示最終給定的斷層深度參數(shù)為20 km。

2.3 模型及反演方法

跨斷層GPS剖面結(jié)果表明（圖3），橫跨紅河斷裂帶、無(wú)量山斷裂帶和楚雄—南華—建水?dāng)嗔褞俣葓?chǎng)平行于斷層分量整體呈現(xiàn)出5～6 mm/a的右旋滑動(dòng)。其中，橫跨無(wú)量山斷裂帶和建水?dāng)嗔褞Т嬖诿黠@的震間反“S”型曲線（圖3b）。垂直于斷層方向分量的離散度比較大，難以確定速度場(chǎng)整體差異，但在橫跨紅河斷裂帶和建水?dāng)嗔褞Т嬖?～3 mm/a的速度差異（圖3b）。

Okada（1985）提出了各向同性彈性半無(wú)限空間位錯(cuò)理論，即彈性體內(nèi)一矩形幾何面發(fā)生滑動(dòng)（走滑、傾滑）所引起某一點(diǎn)的位移，與滑動(dòng)面的滑動(dòng)量成正比關(guān)系，比例系數(shù)由該點(diǎn)與滑動(dòng)面的相對(duì)位置、滑動(dòng)面的幾何尺度、傾角、深度和彈性介質(zhì)性質(zhì)所確定。若地下有多個(gè)矩形滑動(dòng)面，則在某一點(diǎn)所引起的位移，就是這些矩形幾何面各自滑動(dòng)所引起的該點(diǎn)位移的疊加，具體計(jì)算公式參考Okada（1985）的研究。在實(shí)際應(yīng)用Okada斷裂位錯(cuò)模型時(shí)，通常將實(shí)際斷裂帶分解成一系列上部斷層面閉鎖、下部斷層面“無(wú)限”延伸并自由錯(cuò)動(dòng)的斷層段矩形面的組合。每個(gè)斷層段矩形面不僅尺度和傾角可以互不相同，閉鎖深度、錯(cuò)動(dòng)量和錯(cuò)動(dòng)方式也可各有差異。其中，下部斷層面無(wú)限延伸，即所謂的“深斷裂位錯(cuò)模型”（Savage et al，1999；Gan et al，2000，2007），主要是為了消除有限深度的斷層面可能產(chǎn)生的不真實(shí)的邊界效應(yīng)。“深斷裂位錯(cuò)模型”將錯(cuò)動(dòng)矩形面看做無(wú)限向下延伸，可在某種程度上等價(jià)體現(xiàn)深部大范圍粘彈性差異拖曳的效果。

在反演過程當(dāng)中，將各斷層段的運(yùn)動(dòng)參數(shù)向量作為未知數(shù)，斷層的鎖定深度、走向和傾角給定先驗(yàn)值，以GPS觀測(cè)值作為約束，利用最小二乘原理，使GPS速度場(chǎng)的擬合達(dá)到最佳，即殘差平方和最小，從而確定斷層的滑動(dòng)量。此外，為了確定最優(yōu)的鎖定深度，利用遺傳算法，給定不同的閉鎖深度先驗(yàn)值，利用殘差平方和最小原理，確定最優(yōu)鎖定深度。

本文首先以紅河斷裂帶、無(wú)量山斷裂帶和楚雄—南華—建水?dāng)嗔褞У目臻g幾何展布為基礎(chǔ)，將3條斷裂帶簡(jiǎn)單地進(jìn)行分段處理，其中紅河斷裂分為5段，楚雄—南華—建水?dāng)嗔逊譃?段，無(wú)量山斷裂未分段（圖1）。其次，基于斷層的彈性半空間位錯(cuò)模型（Okada，1985），利用1999—2007年和2011—2017年2期速度場(chǎng)反演了各斷層段的滑動(dòng)速率（表2）。其中，斷層傾角、滑動(dòng)速率和斷層鎖定深度均先驗(yàn)給定（表1）。

3 反演結(jié)果與討論

2期速度場(chǎng)的反演結(jié)果如表2和圖4，5所示，其中，圖4a和圖5a為2期速度場(chǎng)觀測(cè)值和模型值，圖4b和圖5b為擬合殘差。圖4和圖5表明，除了紅河斷裂帶北段接近鶴慶—洱源段的少數(shù)點(diǎn)之外，其它點(diǎn)的速度場(chǎng)資料在1σ誤差范圍內(nèi)均得到很好地?cái)M合。表2為2期速度場(chǎng)反演的各斷層走滑和傾滑速率以及閉鎖深度，結(jié)果表明，無(wú)量山斷裂、南華—楚雄—建水?dāng)嗔押图t河斷裂總體以右旋走滑為主，并伴隨有較小的正傾滑分量。其中，南華—楚雄—建水?dāng)嗔褞Ш蜔o(wú)量山斷裂帶的活動(dòng)性要明顯強(qiáng)于紅河斷裂帶。

3.1 紅河斷裂帶

紅河斷裂帶作為一條深大的邊界斷裂（Wang et al，2014），與金沙江斷裂一起構(gòu)成了川滇菱形塊體的西南邊界。本文研究結(jié)果表明，整體而言，紅河斷裂帶南段（哀牢山山前斷裂和中谷斷裂）的活動(dòng)性要強(qiáng)于中、北段。北段的2期資料反演的右旋走滑和傾滑速率分別為（單位：mm/a）：（2.5±1.9）（2.1±1.4）和（-0.7±1.5）（0.2±1.4）；中段的2期右旋走滑和拉張速率分別為（單位：mm/a）：（1.6±1.5）（3.2±1.8）和（-2.5±1.6）（-3.0±1.6）；南段的2期右旋走滑和拉張速率分別為（單位：mm/a）：（3.3±2.0）（4.5±2.0）和（1.2±2.2）（-0.3±2.2）。紅河斷裂帶新活動(dòng)地貌地質(zhì)證據(jù)表明，早更新世以來(lái)，紅河斷裂帶北段、中段、南段的水平速率分別為（單位：mm/a）：1.6，3.1，2.6；中更新世以來(lái)的水平活動(dòng)速率分別為（單位：mm/a）：1.6和2.6和3.2；晚更新世以來(lái)的速率分別為（單位：mm/a）：2.6，2.4和2.4 。垂直活動(dòng)速率為（單位：mm/a）：早更新世北段、中段、南段的位錯(cuò)速率分別為（單位：mm/a）：1.6，0.17和0.11；中更新世速率分別為（單位：mm/a）：1.08，0.23和0.15；晚更新世速率分別為（單位：mm/a）：1.4，0.6和0.6（虢順民等，2013）。本文的結(jié)果與地質(zhì)上和其它研究成果基本一致（虢順民等，2013；王閻昭等，2008）。GPS垂向速度場(chǎng)結(jié)果表明（Pan，Shen，2017），紅河斷裂帶北邊位于川滇塊體西南部的區(qū)域以約1 mm/a的速率抬升，而紅河斷裂帶南部塊體以大約1 mm/a的速率下沉，上述運(yùn)動(dòng)可能形成了紅河斷裂帶約2 mm/a的拉張變形。

3.2 無(wú)量山斷裂帶

受青藏高原隆起影響，滇西南塊體向南運(yùn)動(dòng)，中下地殼廣泛存在的低速層為塊體運(yùn)動(dòng)提供了有利條件，但剛性的臨滄花崗巖體對(duì)其南西運(yùn)動(dòng)起著頂托作用，使得東、西兩側(cè)塊體運(yùn)動(dòng)出現(xiàn)差異，且塊體運(yùn)動(dòng)方向與無(wú)量山斷裂帶呈小角度相交，在此背景下，無(wú)量山斷裂帶表現(xiàn)為水平右旋走滑運(yùn)動(dòng)，起著滑動(dòng)分解應(yīng)變的作用（常祖峰等，2016）。無(wú)量山斷裂帶晚第四紀(jì)以來(lái)活動(dòng)特別明顯，孕育了2007年普洱6.4級(jí)和2014年景谷MS6.6等強(qiáng)震?；?期速度場(chǎng)資料，本研究獲得的無(wú)量山斷裂帶右旋走滑速率和傾滑速率分別為（單位：mm/a）：（5.6±1.5）（4.0±1.6）和（-0.7±1.5）（0.2±1.4）；與活動(dòng)構(gòu)造的6.4～8.7 mm/a總滑動(dòng)速率（虢順民等，1999；常祖峰等，2016）及其它研究成果得到的（4.3±1.1）mm/a基本一致（王閻昭等，2008）。

3.3 南華—楚雄—建水?dāng)嗔褞?/p>

南華—楚雄—建水?dāng)嗔褞г趲缀握共忌吓c紅河斷裂帶近乎平行，前人將該斷裂帶作為哀牢山—紅河斷裂帶中南段的重要分支（常祖峰等，2015），與紅河斷裂帶共同構(gòu)成了川滇菱形塊體的西南邊界。斷裂帶內(nèi)孕育了1680年楚雄6級(jí)和1970通海MS7.7地震等，表明該斷裂帶現(xiàn)今活動(dòng)性比較強(qiáng)。本文基于2期速度場(chǎng)資料給出其平均右旋滑動(dòng)和傾滑速率分別為（單位：mm/a）：（5.7±1.5）（6.4±1.4）和（-0.2±1.7）（-3.5±1.7）；與地質(zhì)上的2.3～4.0 mm/a的右旋滑動(dòng)和0.6～1.1 mm/a的構(gòu)造抬升（王洋等，2015），

以及約4.5 mm/a的大地測(cè)量資料研究結(jié)果（王閻昭等，2008；聞學(xué)澤等，2011）基本一致。與紅河斷裂帶的活動(dòng)性相比，南華—楚雄—建水?dāng)嗔褞У幕顒?dòng)性明顯更強(qiáng)，這與紅河斷裂帶中南段近期的低地震活動(dòng)性和南華—楚雄—建水?dāng)嗔褞У妮^強(qiáng)地震活動(dòng)性現(xiàn)狀完全吻合。前人研究認(rèn)為，位于小江和普渡河斷裂西側(cè)的地殼在運(yùn)動(dòng)過程中不斷發(fā)生近東西向的沖斷、逆掩，造成地殼的隆升和水平縮短。由于這種構(gòu)造作用的轉(zhuǎn)換，小江斷裂帶西盤（川滇塊體）的主動(dòng)向南運(yùn)動(dòng)對(duì)建水?dāng)嗔褞Ь哂虚L(zhǎng)期強(qiáng)烈的作用，建水?dāng)嗔褞б杂倚呋?剪切-橫向縮短/逆沖變形的方式吸收與轉(zhuǎn)換小江斷裂帶西盤的向南運(yùn)動(dòng)，從而導(dǎo)致了小江斷裂帶靠近紅河段活動(dòng)性的驟減和南華—楚雄—建水?dāng)嗔褞У幕顒?dòng)性增強(qiáng)（Wang，et al，1998；聞學(xué)澤等，2011）。

3.4 斷層鎖定深度

利用2期資料，反演得到紅河斷裂帶的平均閉鎖深度為6.8 km，閉鎖程度最強(qiáng)的區(qū)域?yàn)樵?yáng)段（F6）和洱源—彌渡段（F3），與M7圈定的地震危險(xiǎn)區(qū)和Wen等（2016）的結(jié)果相符。南華—楚雄—建水?dāng)嗔褞У钠骄]鎖深度為7.2 km，斷裂北段（F1）的閉鎖深度要高于斷裂南段（F2）；無(wú)量山斷裂帶的平均閉鎖深度為7 km。

3.5 討論

基于1999—2007年和2011—2017年2期GPS速度場(chǎng)觀測(cè)資料，反演了滇中地區(qū)主要活斷層現(xiàn)今滑動(dòng)速率和平均閉鎖深度。反演結(jié)果中，2期資料的結(jié)果存在一定的差異，主要原因可能來(lái)自于觀測(cè)數(shù)據(jù)和大地震的影響。為減小汶川地震對(duì)反演結(jié)果的影響，筆者將數(shù)據(jù)分成了2段，剔除了汶川地震之后一年的觀測(cè)數(shù)據(jù)。因此，這種差異可能主要來(lái)自于數(shù)據(jù)質(zhì)量和觀測(cè)點(diǎn)的空間分布，圖4a，5a給出了研究區(qū)域的觀測(cè)站點(diǎn)分布，紅河斷裂帶西北—中段100 Km范圍內(nèi)缺少觀測(cè)數(shù)據(jù)，很難準(zhǔn)確地約束斷裂帶的震間和細(xì)部變形特征。作為M7（2016—2025地震危險(xiǎn)區(qū)劃圖）圈定的主要地震危險(xiǎn)區(qū)，要進(jìn)一步研究紅河斷裂帶的地震危險(xiǎn)性及其與周邊斷裂帶之間的應(yīng)變特征，還需要布設(shè)更詳細(xì)的近斷層觀測(cè)剖面。

感謝匿名審稿人對(duì)本文提供的指導(dǎo)，本文使用的畫圖軟件為GMT，在此對(duì)研發(fā)團(tuán)隊(duì)表示感謝。也感謝中國(guó)地震局第二監(jiān)測(cè)中心和地質(zhì)所各位老師和師兄的指導(dǎo)和幫助。

參考文獻(xiàn)：

常祖峰，常昊，李鑒林，等.2015.楚雄—南華斷裂晚第四紀(jì)活動(dòng)與1680年楚雄6級(jí)地震[J].中國(guó)地震，31（3）：492-500.

常祖峰，陳曉利，陳宇軍，等.2016.景谷MS6.6地震同震地表破壞特征與孕震構(gòu)造[J].地球物理學(xué)報(bào)，59（7）：2539-2552.

虢順民，向宏發(fā)，計(jì)鳳桔，等.2013.云南紅河斷裂帶地質(zhì)圖：1∶ 50000[M].北京：地震出版社：1-21.

虢順民，向宏發(fā)，周瑞琦，等.1999.滇西南龍陵—瀾淪斷裂帶——大陸地殼上一條新生的破裂帶[J].科學(xué)通報(bào)，44（19）：2118-2121.

韓新民，柴天俊，肖九安，等.1982.石屏—建水?dāng)嗔阎卸蔚男禄顒?dòng)與地震[J].地震研究，5（2）：220-225.

李西，冉勇康，陳立春，等.2016.紅河斷裂帶南段全新世地震活動(dòng)證據(jù)[J].地震地質(zhì)，38（3）：596-604.

劉耀輝，李金平，王劉偉，等.2015.紅河斷裂帶閉鎖程度和滑動(dòng)虧損分布特征研究[J].測(cè)繪工程，24（8）：20-22.

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

徐錫偉，張培震，聞學(xué)澤，等.2005.川西及其鄰近地區(qū)活動(dòng)構(gòu)造基本特征與強(qiáng)震復(fù)發(fā)模型[J].地震地質(zhì)，27（3）：446-461.

王洋，張波，侯建軍，等.2015.曲江斷裂晚第四紀(jì)活動(dòng)特征及滑動(dòng)速率分析[J].地震地質(zhì)，37（4）：1177-1192.

王怡然，李有利，司蘇沛.2015.石屏—建水?dāng)嗔阉藉e(cuò)斷河流的幾何特征與斷層活動(dòng)性研究[J].第四紀(jì)研究，35（1）：100-108.

聞學(xué)澤，杜方，龍鋒，等.2011.小江和曲江—石屏兩斷裂系統(tǒng)的構(gòu)造動(dòng)力學(xué)與強(qiáng)震序列的關(guān)聯(lián)性[J].中國(guó)科學(xué)：地球科學(xué)，41（5）：713-724.

謝卓娟，李山有，呂悅軍，等.2015.滇西南地區(qū)主要活動(dòng)斷裂的b值空間分布特征[J].地球科學(xué)——中國(guó)地質(zhì)大學(xué)學(xué)報(bào)，40（10）：1755-1766.

Allen C R，Han Y.1984.Study of the quaternary activities of the red river fault（Ⅱ）-its features of activity，slip rate and recurrence intervals of earthquakes[J].Journal of Seismological Research，41（1）：170-177.

Altamimi X，Collilieux X，Metivier L.2011.ITRF2008：An Improved Solution of the International Terrestrial Reference Frame[J].Journal of Geodesy，85（8）：457-473.

Boehm J，Heinkelmann R，Schuh H.2007.Short Note：A global Model of Pressure and Temperature for Geodetic Applications[J].Journal of Geodesy，81（10）：679-683.

Boehm J，Werl B，Schuh H.2006.Troposphere Mapping Functions for GPS and Very Long Baseline Interferometry from European Centre for Medium-Range Weather Forecasts Operational Analysis Data[J].J Geophys Res，111（B2）：1-9.

Chen Y，Zhang Z，Sun C，et al.2013.Crustal anisotropy from Moho converted Ps wave splitting analysis and geodynamic implications beneath the eastern margin of Tibet and surrounding regions[J].Gondwana Research，24（3-4）：946-957.

Fang L，Wu J，Wang W，et al.2015.Aftershock observation and analysis of the 2013 MS7.0 Lushan earthquake[J].Seismological Research Letters，86（4）：1135-1142.

Gan W，Svarc J L，Savage J C，et al.2000.Strain accumulation across the Eastern California Shear Zone at latitude 36°30′N[J].J Geophys Res，105（B7）：16229-16236.

Gan W，Zhang P，Shen Z K，et al.2007.Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements[J].J Geophys Res，112（B8）：582-596.

Herring T A，King R W，McClusky S C.2010.GAMIT Reference Manual，Release 10.4[M].Cambridge：Mass.Inst.of Technol.

Hu J，Badal J，Yang H，et al.2018.Comprehensive crustal structure and seismological evidence for lower crustal flow in the southeastern margin of Tibet revealed by receiver functions[J].Gondwana Research，55：42-59.

Lei J，Li Y，Xie F，et al.2014.Pn anisotropic tomography and dynamics under eastern Tibetan plateau[J].J Geophys Res，119（3）：2174-2198.

Leloup P H，Harrison T M，Ryerson F J，et al.1993.Structural，petrological and thermal evolution of a Tertiary ductile strike-slip shear zone，Diancang Shan，Yunnan[J].J Geophys Res，98（B4）：6715-6743.

Leloup P H，Lacassin R，Tapponnier P，et al.1995.The Ailao Shan-Red River shear zone（Yunnan，China），Tertiary transform boundary of Indochina[J].Tectonophysics，251（1-4）：3-10.

Lyard F，Lefevre F，Letellier T，et al.2006.Modelling the Global Ocean Tides：Modern Insights from FES2004[J].Ocean Dynamics，56（5），394-415.

McCaffrey R，Qamar A I，King R W，et al.2007.Fault Locking，Block Rotation and Crustal Deformation in the Pacific Northwest[J].Geophys J Int，169（3）：1315-1340.

McCaffrey R.2005.Block Kinematics of the Pacific-North America Plate Boundary in the Southwestern US from Inversion of GPS，Seismological，and Geologic Data[J].J Geophys Res，110（B7）：1-27.

McCaffrey R.2009.Time-dependent inversion of three-component continuous GPS for steady and transient sources in northern Cascadia[J].Geophys Res Lett，36（36）：2497-2502.

Okada Y.1985.Surface deformation due to shear and tensile faults in a half-space[J].Bulletin of the Seismological Society of America，75（2）：1018-1040.

Pan Y，Shen W B.2017.Contemporary crustal movement of southeastern Tibet：Constraints from dense GPS measurements[J].scientific Reports，7.

Savage J C，Svarc J L，Prescott W H.1999.Geodetic estimates of fault slip rates in the San Francisco Bay area[J].J Geophys Res，104（B3）：4995-5002.

Schoenbohm L M，Burchfiel B C，Liang Z C.2006.Propagation of surface uplift，lower crustal flow，and Cenozoic tectonics of the southeast margin of the Tibetan Plateau[J].Geology，34（10）：813.

Sun X，Ban X，Xu M，et al.2014.Crustal structure beneath SE Tibet from joint analysis of receiver functions and Rayleigh wave dispersion[J].Geophys Res Lett，41（5）：1479-1484.

Tapponnier P，Lacassin R，Leloup P H，et al.1990.The Ailao Shan/Red River metamorphic belt：Tertiary left-lateral shear between Indochina and South China[J].Nature，343（6257）：431-437.

Wang E.1998.Late cenozoic Xianshuihe/Xiaojiang and Red River fault systems of southwestern Sichuan and central Yunnan，China[J].Special Paper of the Geological Society of America，327：1-108.

Wang W，Wu J，F(xiàn)ang L H，et al.2014.S wave velocity structure in southwest China from surface wave tomography and receiver functions.J Geophys Res，119（2），1061-1078.

Wen X Z，Ma S L，Du F，et al.2016.Fault Activity，Seismicity and GPS Deformation of the Seismic Gap along the Red River Fault Zone（RRFZ）in Yunnan，China[C]//EGU General Assembly Conference Abstracts，18：3202.

Zhu M，Graham S，Mchargue T.2009.The Red River Fault zone in the Yinggehai Basin，South China Sea[J].Tectonophysics，476（3-4）：397-417.

Abstract To analyze characteristics of the main active faults in the central Yunan province，we inverted slip rates and lock depths of the main fault zones according to Okadas fault dislocation model based on GPS observations during 1999-2007，2011-2017 and geologic slip rates.Our results are as follows：（1）the strike-slip rate of the Red River fault is about（1.5±1.6）～（4.7±1.5）mm/a and the dip-slip rate is（-3.6±1.6）～（1.9±2.4）mm/a，while the southern segment is more active than other parts.（2）the strike-slip rates of the Nanhua-Chuxiong-Jianshui fault and the Wuliangshan fault are（4.8±1.4）～（6.6±1.6）mm/a and（4.0±1.6）～（5.6±1.5）mm/a，respectively.And the dip-slip rates are（-5.8±1.5）～（1.7±1.8）mm/a and（-0.7±1.5）～（0.2±1.4）mm/a，respectively.（3）the Yuanjiang-Yuanyang and Eryuan-Midu segments of the Red River fault，the Wuliangshan fault zone，and the western part of the Nanhua-Chuxiong-Jianshui fault are locked，with locking depths of 6.8 km，7 km，and 7.2 km，respectively.

Keywords：fault slip rates；locking depth；GPS；fault activity