唐紅濤 張學(xué)剛 高子讓 丁國麗 李瑞莎

1）中國地震局第二監(jiān)測中心，西安 710054

2）青海省測繪質(zhì)量監(jiān)督檢驗中心，西寧 810001

3）湖南省第二測繪院，長沙 410118

4）國家測繪地理信息局第一航測遙感院，西安 710054

0 引言

我國大陸地殼運動的主要動力源自印度板塊北推，當遇到堅硬的阿拉善、鄂爾多斯地塊，地殼物質(zhì)便由鄂爾多斯地塊西南緣向南運移，從而構(gòu)成了我國中西部地震多發(fā)的南北地震帶這一地殼活動的重要組成部分（張培震等，2013；王雙緒等，2005、2012；張希等，2004；李瑞莎等，2012；蘇琴，2005）。該構(gòu)造過程中的青藏塊體東北緣處于青藏高原、阿拉善地塊以及鄂爾多斯地塊等相互碰撞、擠壓的物質(zhì)遷移交匯區(qū)，對該地區(qū)地殼運動的監(jiān)測有助于分析研究我國地殼整體運動特征。2008年汶川8.0級地震后，南北地震帶進入了地震高發(fā)期，斷層活動頻繁，2017年該地區(qū)發(fā)生了門源6.4級地震，因此，對青藏塊體東北緣的觀測研究顯得尤為重要。

中國地震局第二監(jiān)測中心自20世紀80年代起，通過跨斷層水準、基線等手段對該地區(qū)進行了斷層活動性監(jiān)測，至今已累積了長達30多年的觀測資料。本文綜合運用合成速率法與趨勢累積率，對祁連山、海原、六盤山斷裂帶跨斷層水準資料進行了處理分析，進而揭示了門源地震前及鄰近期內(nèi)青藏塊體東北緣的地殼運動特征與斷層活動情況。

1 區(qū)域構(gòu)造運動特征

青藏塊體東北緣內(nèi)部自北向南發(fā)育著NWW-NW向的新構(gòu)造時期以來的活動斷裂帶，主要包括：祁連山北緣斷裂帶、昌馬-俄博斷裂、冷龍嶺斷裂、陰凹槽斷裂、香山-天景山斷裂、海原斷裂及六盤山斷裂等（圖1）。受印度板塊北推碰撞歐亞大陸的主動力驅(qū)使，加之受到相對穩(wěn)定的阿拉善地塊、鄂爾多斯地塊相當力度的阻擋（丁國瑜，1991；虢順民等，2000；李清河等，1999；張國民等，2004；魯如魁等，2008；吳滿路等，2008；張先康等，2008，方穎等，2009；姜衛(wèi)平等，2008；吳云等，2002；朱文耀等，1999；張強等，2000；國家地震局地質(zhì)研究所等，1993），青藏塊體東北緣活動斷裂帶地質(zhì)構(gòu)造活動以NE向擠壓-逆沖和走滑-旋轉(zhuǎn)（左旋走滑）為主要特征。

整個青藏塊體東北緣地區(qū)布設(shè)了40個跨斷層監(jiān)測點（圖1），監(jiān)測場地自西向東分布較為均勻，每年觀測3個周期（分別在3、7、11月）。圖1中數(shù)據(jù)的截止時間為2016年11月，紅色標注的跨斷層場地為中-長期異常監(jiān)測點，藍色場地為無異常監(jiān)測點。

圖1 青藏塊體東北緣監(jiān)測點分布

2 方法與處理分析

2.1 合成速率法

本文使用的合成速率法是將同一斷裂帶上的異常監(jiān)測點進行對比，從整體上分析該斷裂活動的總體特征，把多個場地數(shù)據(jù)經(jīng)該方法處理后繪制在同一個圖上，觀察其活動性質(zhì)是否具有一致性，從而獲取斷裂帶附近區(qū)域應(yīng)變積累狀況。若研究區(qū)域再擴大，可將多條斷裂進行速率合成以對比分析。監(jiān)測點第j年的平均年變化速率為（孫軍等，2014；陳兵等，2000；王永安等，2011）

把同一斷裂上的跨斷層監(jiān)測點的年平均速率的均值視為該斷裂運動的合成速率，則斷裂第j年合成速率為

2.2 解算分析

根據(jù)合成速率理論，本文對青藏塊體東北緣地區(qū)祁連山斷裂帶、海原斷裂帶及六盤山斷裂帶進行了合成速率解算分析，并繪制了每個測點的年變化速率（圖2～5）及每個斷裂的合成速率時序圖（圖6），以研究同一斷層不同時段與不同斷層之間的運動差異特征。為了較好地反映斷層的活動特性，本文所選資料均為以往觀測穩(wěn)定的跨斷層短水準數(shù)據(jù)。

圖2 祁連山斷裂帶西段測點年平均活動速率

圖3 祁連山斷裂帶東段測點年平均活動速率

2.2.1 祁連山斷裂帶西段

由于祁連山斷裂帶EW向空間跨度較大，本文計算與分析時將其分為祁連山斷裂帶西段與東段。由圖2可知，紅柳峽場地開始觀測時間較晚，4個場地年平均活動速率在時間軸上表現(xiàn)較為一致（1997、2001、2004、2008、2012年等時間節(jié)點變化規(guī)律明顯一致），玉門5.9級、汶川8.0級地震前，曲線顯示了斷裂由張—壓的活動特征，2010～2016年間，多個測點呈現(xiàn)張性運動狀態(tài)。

圖4 海原斷裂帶測點年平均活動速率

圖5 六盤山斷裂帶測點年平均活動速率

2.2.2 祁連山斷裂帶東段

祁連山斷裂帶東段扁都口場地B-水A2測段表現(xiàn)為運動速率起伏較大，但其起伏態(tài)勢與其他測段有多點的時間一致性（圖3）。從該5處場地曲線整體來看，在2010～2016年間，多表現(xiàn)為斷裂由負轉(zhuǎn)正的變化規(guī)律，尤其從2016年門源6.4級地震前2年開始，統(tǒng)一呈現(xiàn)斷裂張—壓的活動特征。

2.2.3 海原斷裂帶

由圖4可知，該斷裂帶活動在時間上有顯著的分段性變化規(guī)律，1981～1999年斷層活動較弱，1999～2005年斷層活動較為活躍，其間發(fā)生了景泰5.9級、民樂6.1級地震，2010～2016年，異常測段均表現(xiàn)為由張—壓的運動態(tài)勢，曲線線性趨勢走向顯著，此區(qū)間發(fā)生了岷縣漳縣6.6級、門源6.4級地震，目前海原斷裂帶仍趨向壓性活動變化。

2.2.4 六盤山斷裂帶

該斷裂帶中長期異常監(jiān)測點分布相對集中（圖1），其年平均活動速率時序曲線有顯著的一致性（圖5）。自1981年開始觀測至2016年，該斷裂帶多處監(jiān)測點在許多時點呈現(xiàn)斷裂活動特征的明顯一致性（如1982、1983、1995、1996、2000、2001、2004、2010年等）。從圖5可知，永登5.8級、景泰5.9級、玉樹7.1級地震同震反應(yīng)較大，均出現(xiàn)了震后張、拉交替變化現(xiàn)象。但在2016年1月門源6.4級地震前，多數(shù)測段出現(xiàn)了壓—張的運動特征。門源地震后，六盤山斷裂帶大多數(shù)異常測段年平均活動速率轉(zhuǎn)折回升，這可能是此次地震同震反應(yīng)的結(jié)果。從量值上可以看出，六盤山斷裂帶監(jiān)測點年平均活動速率變化幅值大于另外2條斷裂，壓性最大約4.8mm／a，張性最大約2.8mm／a。

圖6 各斷裂帶合成速率

由多個監(jiān)測點運動速率合成為同一條斷裂的速率時序圖（圖6），從圖6可知，2000～2010年是祁連山斷裂帶、海原斷裂帶、六盤山斷裂帶，亦即整個青藏塊體東北緣地區(qū)的斷層活躍期。2010～2016年，祁連山斷裂帶整體偏向微弱的張性運動；海原斷裂帶由負轉(zhuǎn)正，張性變壓性；六盤山斷裂帶在門源地震前顯示張性運動，震后轉(zhuǎn)折。

2.3 趨勢累積率

趨勢累積率可用以分析斷裂的活動能量和應(yīng)變的聚集狀況，它反映了自有資料記錄以后監(jiān)測點的趨勢形變累積情況以及趨勢累積量歸一化后的整個研究區(qū)域地殼應(yīng)變強度比。大范圍的斷層形變趨勢累積率的時空演化過程可以代表區(qū)域地殼形變的分布狀況，它實際上反映的是一種廣義應(yīng)變的累積速率趨勢強度比（李媛等，2016），其計算公式為

式中，hi為斷層兩盤測點之間的相對高差觀測值；n為觀測期數(shù)。根據(jù)定義，Dc絕對值為0～1，當Dc絕對值為0時，形變不存在趨勢性變化；當Dc絕對值為1時，形變記錄的趨勢性最強。Dc絕對值具有自然歸一性，這為建立不同演化程度的數(shù)據(jù)序列之間的可比性標度提供了基礎(chǔ)。

2.4 數(shù)據(jù)處理與分析

圖7 祁連山斷裂帶（a）和海原-六盤山斷裂帶（b）跨斷層水準形變趨勢累積率等值線分布

利用中長期趨勢性異常的跨斷層資料（自有觀測以來的全部數(shù)據(jù)），計算了各場地的Dc值，并繪制了祁連山斷裂帶與海原-六盤山斷裂帶的Dc等值線圖（圖7）。由圖7可見，祁連山斷裂帶（圖7（a））垂直形變累積率高值區(qū)與高梯度帶主要分布在中、西段，東段相對較穩(wěn)定；海原-六盤山斷裂帶（圖7（b））Dc值高梯度帶主要分布于六盤山地區(qū)，整體來看，海原-六盤山斷裂帶NE方向應(yīng)力積累效應(yīng)明顯強于其WS方向。這是長期地殼應(yīng)力積累的結(jié)果，其間必然經(jīng)歷了多次地震能量的釋放與再累積的交替過程。

3 結(jié)語

通過將速率合成、趨勢累積率等方法與青藏塊體東北緣跨斷層短水準資料相結(jié)合，可獲得以下初步結(jié)論：

（1）由單個斷裂多個異常監(jiān)測點的年平均活動速率時序曲線可以獲知，用速率合成方法分析跨斷層的監(jiān)測數(shù)據(jù)，可以較好地找到某一條斷裂在時空分布上的具體活動規(guī)律（活躍期、張、拉等變化過程）；而由多個斷裂集中對比分析與跨斷層形變趨勢累積率，則可獲得更廣時間與空間的地殼形變運動特征。

（2）從合成速率結(jié)果可知，2010～2016年，祁連山斷裂帶呈微弱拉張的運動態(tài)勢；海原斷裂帶由張轉(zhuǎn)壓；六盤山斷裂帶在門源地震前顯示張性運動，震后發(fā)生轉(zhuǎn)折。趨勢累積率則從長期觀測數(shù)據(jù)的角度顯示了祁連山斷裂帶中、西段累積應(yīng)變略高于東段，從量值上看整個青藏塊體東北緣Dc高值區(qū)位于六盤山斷裂帶附近。