吳微微 梁明劍 龍鋒 蘇金蓉 陳學(xué)芬

摘要：近年來(lái)，青川斷裂及鄰區(qū)多次發(fā)生4.0級(jí)以上地震，為進(jìn)一步研究該區(qū)現(xiàn)今地震活動(dòng)性，結(jié)合現(xiàn)場(chǎng)地震地質(zhì)調(diào)查和無(wú)人機(jī)攝影方法測(cè)繪構(gòu)造地貌，揭示青川斷裂的構(gòu)造變形信息；采用CAP方法反演2008年5月17日至2021年12月31日研究區(qū)32次MW≥4.0地震的震源機(jī)制解，進(jìn)一步討論斷裂的構(gòu)造變形特征；使用ω2模型擬合震源譜方法獲得研究區(qū)十年尺度的應(yīng)力降動(dòng)態(tài)演化，分析3次強(qiáng)震前后研究區(qū)的應(yīng)力調(diào)整過(guò)程；通過(guò)三維P波速度結(jié)構(gòu)反演，探討研究區(qū)深部孕震背景；最后綜合地質(zhì)學(xué)和地震學(xué)研究結(jié)果分析青川斷裂及鄰區(qū)的現(xiàn)今活動(dòng)性。結(jié)果表明：①地震地質(zhì)調(diào)查和震源機(jī)制解的結(jié)果均反映青川斷裂的運(yùn)動(dòng)性質(zhì)為走滑兼擠壓。②2008年汶川8.0級(jí)和2013年蘆山7.0級(jí)強(qiáng)震前、后，研究區(qū)的應(yīng)力水平在一段時(shí)間內(nèi)呈現(xiàn)出升高-回落的調(diào)整過(guò)程。③體波層析成像的結(jié)果顯示研究區(qū)的地殼淺層高低速異常表現(xiàn)以龍門山斷裂帶為分界線的空間分布特征。④青川斷裂具有潛在地震發(fā)生的構(gòu)造背景與能力，但目前該斷裂及鄰區(qū)的地震活動(dòng)正處于強(qiáng)震序列衰減階段，M4.0以上地震活動(dòng)呈逐漸減弱的趨勢(shì)，高應(yīng)力降地震事件明顯減少，區(qū)域整體應(yīng)力水平較低。

關(guān)鍵詞：青川斷裂及鄰區(qū)；汶川地震；震源機(jī)制解；應(yīng)力降；速度結(jié)構(gòu)；地震活動(dòng)性

中圖分類號(hào)：?? 文獻(xiàn)標(biāo)識(shí)碼：?? 文章編號(hào)：1000-0666（2023）02-0188-16

doi：10.20015/j.cnki.ISSN1000-0666.2023.0035

0 引言

青川斷裂是龍門山斷裂帶北段的分支斷裂之一，與北川—映秀斷裂呈左階排列，西起平武東部，經(jīng)青川、寧強(qiáng)至勉縣境內(nèi)，全長(zhǎng)約250 km（樊春等，2008；杜建軍等，2013）。2008年5月12日四川汶川8.0級(jí)地震前，沿青川斷裂發(fā)生的地震活動(dòng)頻度低、震級(jí)小，區(qū)域研究程度相對(duì)較低。地震后，主震在北川—映秀斷裂上產(chǎn)生的地表破裂帶越過(guò)了虎牙斷裂帶，向北東終止于青川斷裂南側(cè)（Xu et al，2009），而余震活動(dòng)一直延伸到青川縣城附近，并切過(guò)青川斷裂繼續(xù)向NE延展近40 km。最大余震2008年5月25日青川6.4級(jí)地震就發(fā)生在余震條帶與青川斷裂的交匯處。之后，隨著2013年蘆山、2017年九寨溝兩次7.0級(jí)地震的相繼發(fā)生，青藏高原東緣活動(dòng)斷裂的地震危險(xiǎn)性更受關(guān)注（Toda et al，2008；單斌等，2009；程佳等，2018；盛書中等，2019），青川斷裂開始受到廣泛重視。

然而，地震學(xué)家對(duì)青川斷裂活動(dòng)性及地震危險(xiǎn)性的認(rèn)識(shí)存在分歧。汶川8.0級(jí)地震前的地震地質(zhì)調(diào)查結(jié)果認(rèn)為，岷山隆起對(duì)龍門山斷裂帶北段具有屏障作用，導(dǎo)致龍門山斷裂北段晚第四紀(jì)不再有新活動(dòng)（鄧起東等，1994；李傳友等，2004），而青川斷裂作為龍門山斷裂帶北段的后緣斷裂未來(lái)地震活動(dòng)性不強(qiáng)。汶川地震后的一些研究成果則認(rèn)為青川斷裂是全新世活動(dòng)斷裂，存在多期古地震事件（Sun，2015；Lin et al，2012，2014，2016）。Liang等（2018）通過(guò)古地震研究，認(rèn)為青川斷裂受岷山斷塊隆起的影響，其發(fā)震能力減弱，斷裂的變形主要是對(duì)北川—映秀斷裂強(qiáng)震活動(dòng)的同震響應(yīng)所產(chǎn)生的近地表變形。

震源機(jī)制解結(jié)果顯示青川斷裂附近的地震類型具有復(fù)雜性。發(fā)生在該區(qū)域的余震既有走滑型，也有逆沖型（鄭勇等，2009；Long et al，2012；易桂喜等，2012），汶川8.0級(jí)主震及其44個(gè)強(qiáng)余震則以逆沖型為主（胡幸平等，2008）。2013年蘆山7.0級(jí)地震后，該區(qū)域也多次發(fā)生M≥4.0地震，其地震類型同樣具有多樣性：2017年7月17日青川MS4.9地震為帶逆沖分量的走滑型，發(fā)震空間位置偏離汶川余震密集帶；9月30日的青川MS5.4地震為帶走滑分量的逆沖型，發(fā)生在汶川余震密集帶；10月13日青川ML4.3地震為純走滑型，空間上也偏離余震密集帶。

前人的研究表明，以青川為界其北側(cè)莫霍面深度更接近于四川盆地，而南側(cè)更類似于青藏高原（龍鋒等，2011）。龍門山斷裂帶及鄰區(qū)的地殼各向異性研究結(jié)果表明，青川斷裂周邊區(qū)域應(yīng)力場(chǎng)方向復(fù)雜（石玉濤等，2013）。這也許意味著青川斷裂帶及鄰區(qū)的地震活動(dòng)不僅受區(qū)域構(gòu)造應(yīng)力場(chǎng)的影響，也與其局部復(fù)雜構(gòu)造有一定關(guān)系。

本文以青川斷裂作為研究對(duì)象，通過(guò)詳細(xì)的現(xiàn)場(chǎng)地震地質(zhì)調(diào)查，采用無(wú)人機(jī)攝影技術(shù)對(duì)斷裂沿線的構(gòu)造變形地貌進(jìn)行測(cè)繪，獲得斷裂活動(dòng)的構(gòu)造變形信息；通過(guò)震源機(jī)制解算，獲得對(duì)震源區(qū)局部斷層性質(zhì)的認(rèn)識(shí)，進(jìn)而討論青川斷裂的構(gòu)造變形特征；使用雙差地震層析成像方法（tomoDD）（Zhang，Thurber，2003，2006）對(duì)研究區(qū)域中小地震進(jìn)行重新定位，結(jié)合震源應(yīng)力降繪制汶川余震區(qū)十年尺度的應(yīng)力分布圖像，討論汶川、蘆山和九寨溝地震對(duì)青川斷裂及鄰區(qū)應(yīng)力水平的影響；通過(guò)對(duì)青川斷裂及鄰區(qū)速度結(jié)構(gòu)的反演，探討深部孕震背景；最后通過(guò)地質(zhì)學(xué)和地震學(xué)方法綜合分析青川斷裂及鄰區(qū)的現(xiàn)今活動(dòng)性，為深入研究龍門山斷裂帶北段深部介質(zhì)狀態(tài)和強(qiáng)震孕育、發(fā)生機(jī)理提供基礎(chǔ)支撐。

1 區(qū)域主要斷裂

龍門山斷裂帶主要由3條主干斷裂構(gòu)成：后山斷裂、主中央斷裂、前山斷裂（圖1）。后山斷裂由青川斷裂（F1）、茂汶—汶川斷裂和耿達(dá)—隴東斷裂構(gòu)成；中央斷裂由北川—映秀斷裂（F2）、茶壩—林庵寺斷裂（F3）與鹽井—五龍斷裂構(gòu)成；前山斷裂由灌縣—江油斷裂（F4）、江油—廣元斷裂和雙石—大川斷裂構(gòu)成。

北川—映秀斷裂（F2）線性特征明顯，沿?cái)嗔寻l(fā)育斷錯(cuò)地貌、斷層槽溝和斷塞塘，晚更新世晚期以來(lái)有強(qiáng)烈活動(dòng)（陳國(guó)光等，2007；Gao et al，2016）。探槽結(jié)果揭露該斷裂平均強(qiáng)震復(fù)發(fā)周期為2 960—2 865 a（Ran et al，2010，2013），以兼有右旋走滑分量的逆斷層型破裂為主（徐錫偉等，2008）。

茶壩—林庵寺斷裂（F3）是右旋走滑斷裂，在衛(wèi)星影像上主要表現(xiàn)為陡崖、斷層三角面、斷層槽地和斷層眉脊等構(gòu)造地貌（謝小平等，2019）。汶川地震前的研究認(rèn)為該斷裂沒(méi)有晚更新世以來(lái)的地表活動(dòng)證據(jù)（李傳友等，2004；陳國(guó)光等，2007）。汶川地震中該斷裂發(fā)生破裂，以走滑為主（Xu et al，2009；Shen et al，2009；Feng et al，2017），具有右旋擠壓特征（謝小平等，2019）。汶川地震后，地震學(xué)家通過(guò)開挖探槽揭露出包括汶川地震在內(nèi)的3次地震（Ran et al，2014）。

灌縣—江油斷裂（F4）主要發(fā)育在三疊系和侏羅系地層，汶川地震后沿該斷裂連續(xù)展布的地表破裂帶長(zhǎng)約72 km，最長(zhǎng)可達(dá)90 km，顯示出典型的純逆斷層型地表破裂特征（徐錫偉等，2008）。

2 現(xiàn)場(chǎng)地震地質(zhì)調(diào)查

青川斷裂呈NE走向，幾何結(jié)構(gòu)連貫，線性地貌顯著，沒(méi)有明顯階區(qū)，顯示出良好的連續(xù)性和高成熟度（李傳友等，2004；Lin et al，2012；Sun，2015）。與龍門山斷裂帶的其它斷裂相比，青川斷裂的走滑運(yùn)動(dòng)作用更為顯著，自中新世晚期以來(lái)平均右旋滑動(dòng)速率約為1.5 mm/a（Wang et al，2008；孫浩越，2015）。汶川8.0級(jí)地震之后，人們通過(guò)探槽研究揭露青川斷裂呈張性運(yùn)動(dòng)性質(zhì)，并認(rèn)為這是斷裂全新世活動(dòng)的證據(jù)（Lin et al，2014，2016；Sun，2015），而Liang等（2018）認(rèn)為這種張性變形是青川斷裂對(duì)北川—映秀斷裂強(qiáng)震活動(dòng)的同震響應(yīng)所產(chǎn)生的地表變形，并不是青川斷裂自身最新活動(dòng)的表現(xiàn)。

為了進(jìn)一步確定青川斷裂的運(yùn)動(dòng)性質(zhì)，我們?cè)谇啻〝嗔焉线x了9個(gè)典型測(cè)點(diǎn)開展詳細(xì)的地震地質(zhì)調(diào)查，并采用無(wú)人機(jī)攝影技術(shù)對(duì)斷裂沿線的構(gòu)造變形地貌進(jìn)行測(cè)繪，結(jié)合地貌點(diǎn)附近的地質(zhì)剖面，揭示青川斷裂的運(yùn)動(dòng)特性。測(cè)點(diǎn)位置及信息如圖2、表1所示，各測(cè)點(diǎn)的典型構(gòu)造地貌、無(wú)人機(jī)航拍及出露地質(zhì)剖面如圖3～6所示。

測(cè)點(diǎn)①位于青川縣城西蒿溪鄉(xiāng)大灣里村一帶，無(wú)人機(jī)獲取的構(gòu)造地貌顯示青川斷裂切過(guò)山腰和山前洪積扇，形成反向槽谷和線性良好的坡折坎，呈現(xiàn)出走滑斷層的典型地貌（圖3a、b）。測(cè)點(diǎn)②位于大灣里村東土地廟的斷層剖面，剖面顯示斷裂發(fā)育于志留系千枚巖與前泥盆系片巖之間，斷裂帶由片理化帶、構(gòu)造角礫、斷層泥等組成，斷裂運(yùn)動(dòng)性質(zhì)為走滑兼擠壓（圖3c）。

測(cè)點(diǎn)③、④、⑤分別是青川斷裂在青川縣城西光輝村、小西壩和喬莊河河床一帶的出露剖面（圖4）。其中，西光輝村附近斷裂發(fā)育于志留系千枚巖與前泥盆系片巖之間，斷層破碎帶寬約6.0～8.0 m，由片理化帶、構(gòu)造角礫、斷層泥帶等組成，顯示斷裂運(yùn)動(dòng)性質(zhì)為走滑兼擠壓（圖4a）。小西壩附近的斷裂呈高角度，產(chǎn)狀N55°E/NW∠75°，斷層帶由片理化帶、斷層泥組成，顯示斷裂運(yùn)動(dòng)性質(zhì)為剪切滑動(dòng)（圖4b）。喬莊河河床附近的斷裂破碎帶寬3.0～4.0 m，主要由剪切片理化帶和構(gòu)造角礫組成，同樣顯示斷裂運(yùn)動(dòng)性質(zhì)為走滑兼擠壓（圖4c）。

測(cè)點(diǎn)⑥位于青川縣城黃土包村一帶，無(wú)人機(jī)獲取的地形地貌為線性的坡折坎和斷層槽谷（圖5a）。測(cè)點(diǎn)⑦是黃土包村西側(cè)發(fā)現(xiàn)的斷層剖面（圖5a中白色正方形位置），主斷裂面寬約0.5～1.0 m，主要由片理化帶和構(gòu)造角礫組成（圖5b），顯示斷裂運(yùn)動(dòng)性質(zhì)為走滑兼擠壓。

測(cè)點(diǎn)⑧位于廣元市路家溝村附近，無(wú)人機(jī)獲取的地形地貌為線性良好的斷層槽谷（圖6a、b），是走滑斷層的典型地貌。測(cè)點(diǎn)⑨是青川斷裂在路家溝東約2.5 km的出露剖面（圖6c），斷裂帶寬約7.0～8.0 m，發(fā)育于志留系千枚巖和遠(yuǎn)古界青白口系變質(zhì)砂巖之間，斷層破碎帶由片理化帶、構(gòu)造角礫巖、碎裂巖等組成，顯示斷裂運(yùn)動(dòng)性質(zhì)為走滑兼擠壓。

綜上所述，青川斷裂在地貌上表現(xiàn)為線性良好的斷層槽谷、坡折帶、斷層埡口等典型的走滑斷層特征，這些線性地貌主要發(fā)育于高臺(tái)地或山坡、山腰處，而斷裂沿線的低階地未見(jiàn)發(fā)育。出露的斷裂剖面顯示，斷裂帶主要由片理化帶、斷層泥為主，局部見(jiàn)構(gòu)造角礫，同樣顯示斷裂以走滑運(yùn)動(dòng)為主，兼有一定逆沖分量的運(yùn)動(dòng)特征。

3 地震學(xué)分析

基于現(xiàn)場(chǎng)地震地質(zhì)調(diào)查獲得的研究區(qū)構(gòu)造變形信息，本文結(jié)合震源機(jī)制解所反映的中強(qiáng)地震錯(cuò)動(dòng)特征，進(jìn)一步討論兩種方法所獲得的研究區(qū)構(gòu)造變形差異。再綜合汶川、蘆山、九寨溝等多次強(qiáng)震前后研究區(qū)應(yīng)力降動(dòng)態(tài)分布演化過(guò)程和三維P波速度結(jié)構(gòu)特征，綜合探討青川斷裂及鄰區(qū)的現(xiàn)今地震活動(dòng)性。

3.1 震源機(jī)制解結(jié)果

本文利用四川數(shù)字地震臺(tái)網(wǎng)最新觀測(cè)資料，挑選信噪比較高、初動(dòng)清晰、震中距小于250 km的寬頻帶地震臺(tái)站波形數(shù)據(jù)，采用CAP方法（Zhao，Helmberger，1994；Zhu Helmberger，1996；Tan et al，2006）反演2008年5月17日至2021年12月31日青川斷裂及鄰區(qū)32次MW≥4.0地震的震源機(jī)制解，見(jiàn)表2。表中編號(hào)1～14是2008年的地震，編號(hào)15～23是2009—2013年的地震，編號(hào)24～32是2014—2021年的地震，震中位置為采用tomoDD方法（Zhang，Thurber，2003，2006）重新定位后的結(jié)果。根據(jù)斷層面解和活動(dòng)斷裂分布推測(cè)所得的發(fā)震斷裂列于表2。

計(jì)算時(shí)，將區(qū)域近震波形旋轉(zhuǎn)后分解為體波（Pnl）和面波，Pnl部分取30 s窗長(zhǎng)截取波形、采用0.05～0.20 Hz的濾波頻帶，面波部分取60 s窗長(zhǎng)截取波形、采用0.05～0.1 Hz的濾波頻帶，使用4階巴特沃斯帶通濾波器濾波（Zoback，1992；Pasyanos et al，1996），選用F-K方法（Zhao，Helmberger，1994；Zhu Helmberger，1996；Tan et al，2006）計(jì)算格林函數(shù)，并最終通過(guò)網(wǎng)格搜索獲得每個(gè)地震的最優(yōu)斷層面解和主應(yīng)力軸。經(jīng)過(guò)比測(cè)發(fā)現(xiàn)，編號(hào)1～14的震源機(jī)制解結(jié)果與前人的研究相近（鄭勇等，2009；Long et al，2012；易桂喜等，2012），驗(yàn)證了計(jì)算結(jié)果的可靠性。

根據(jù)應(yīng)力軸仰角參數(shù)劃分?jǐn)鄬有再|(zhì)（Zoback，1992）得出，32次地震中有11次中強(qiáng)地震為走滑型（P軸仰角＜40°、T軸仰角≤20°或P軸仰角≤20°、T軸仰角＜40°），17次地震為逆斷型地震（P軸仰角≤35°且T軸仰角≥52°），1次地震（編號(hào)2）為帶走滑分量的逆斷型地震（P軸仰角≤20°且T軸仰角在40°～50°），另有3次地震（編號(hào)18，19，22）不在Zoback（1992）分類標(biāo)準(zhǔn)中（表2，圖7）。

基于地質(zhì)構(gòu)造、活動(dòng)斷裂幾何展布及地震活動(dòng)分布特征，將青川斷裂及鄰區(qū)劃分成3個(gè)構(gòu)造單元（圖7）：其中，A區(qū)是北川—映秀斷裂（F2）、灌縣—江油斷裂（F4）和茶壩—林庵寺斷裂（F3）交會(huì)的區(qū)域，為汶川8.0級(jí)地震主破裂區(qū)的東段；B區(qū)是青川斷裂（F1）北東段，為龍門山斷裂帶北段的余震密集區(qū)；C區(qū)是汶川余震序列向北東的擴(kuò)展區(qū)。2008年汶川8.0級(jí)地震以來(lái)，已發(fā)生的中強(qiáng)地震主要集中在A區(qū)和B區(qū)，地震活動(dòng)優(yōu)勢(shì)分布與構(gòu)造線的北東向延伸方向比較一致（圖1）。

在汶川8.0級(jí)地震主破裂區(qū)的東段（A區(qū)），震源機(jī)制解類型（表2，圖7）主要為：北川—映秀斷裂（F2）、灌縣—江油斷裂（F4）附近的地震，以逆沖錯(cuò)動(dòng)為主，未在分類標(biāo)準(zhǔn)中的一次地震（編號(hào)19）也具有較大的逆沖分量，震源機(jī)制解的結(jié)果與地質(zhì)考察顯示的北川—映秀斷裂（F2）、灌縣—江油斷裂（F4）的逆斷層性質(zhì)比較吻合（徐錫偉等，2008）；茶壩—林庵寺斷裂（F3）上發(fā)生的5次中強(qiáng)地震既有逆斷型（編號(hào)23、27、31）也有走滑型（編號(hào)26、28），與該斷裂以走滑運(yùn)動(dòng)為主的破裂特征和右旋兼擠壓性質(zhì)比較一致（Xu et al，2009；Shen et al，2009；Feng et al，2017；謝小平等，2019）；在青川斷裂（F1）、北川—映秀斷裂（F2）和虎牙斷裂（F5）圍限的三角形區(qū)域內(nèi)，既有走滑型地震（編號(hào)13、20）、也有逆斷型（編號(hào)25、29）和未在分類標(biāo)準(zhǔn)中的一次地震（編號(hào)18），推測(cè)可能的發(fā)震斷裂是區(qū)域衍生的次級(jí)斷裂。

在青川斷裂（F1）北東段附近（B區(qū)），已發(fā)生的中強(qiáng)地震既有走滑型，也有帶走滑分量的逆斷型和純逆斷型地震（表2，圖7），表現(xiàn)出較多種類的斷層面解和破裂特征。野外調(diào)查發(fā)現(xiàn)青川斷裂的破碎帶物質(zhì)以剪切滑動(dòng)運(yùn)動(dòng)形成的片理化帶、斷層泥為主，局部見(jiàn)擠壓作用形成的構(gòu)造角礫，具有以走滑為主、兼一定擠壓逆沖分量的錯(cuò)動(dòng)特點(diǎn)（表1，圖3），推測(cè)B區(qū)發(fā)生的走滑型余震可能是青川斷裂的剪切運(yùn)動(dòng)特征在地震學(xué)上的???? 反映。而B區(qū)帶走滑分量的逆斷型地震和逆斷型地震也許可以從水壓致裂地應(yīng)力研究中得到解釋。地應(yīng)力研究認(rèn)為汶川地震影響了龍門山斷裂帶的構(gòu)造應(yīng)力場(chǎng)，在近地表至上地殼15 km深度范圍內(nèi)，北川—青川一帶的震后地應(yīng)力狀態(tài)顯示為逆走滑型，在15～20 km深度范圍內(nèi)該區(qū)域的震后應(yīng)力狀態(tài)由逆沖型占主導(dǎo)（豐成君等，2018）。綜上所述，現(xiàn)場(chǎng)地震地質(zhì)調(diào)查和震源機(jī)制解結(jié)果表明青川斷裂的運(yùn)動(dòng)性質(zhì)以走滑兼擠壓為主，沒(méi)有表現(xiàn)出明顯的張性運(yùn)動(dòng)特征。

C區(qū)是汶川余震序列向北東的擴(kuò)展區(qū)（表2，圖7），地震活動(dòng)以逆沖錯(cuò)動(dòng)為主?！吨袊?guó)及鄰近地區(qū)地震構(gòu)造圖》（徐錫偉等，2016）顯示，這里存在一條NE向不知名次級(jí)斷裂，易桂喜等（2012）的研究指出汶川主震后調(diào)整過(guò)程中發(fā)生較多的逆沖型地震是應(yīng)力在斷裂帶垂向方向上的補(bǔ)充和協(xié)調(diào)。

3.2 強(qiáng)震前后的應(yīng)力降變化

汶川地震后，龍門山斷裂帶及鄰區(qū)相繼發(fā)生蘆山和九寨溝2次7.0級(jí)地震，強(qiáng)震活動(dòng)在周邊斷層上引起的庫(kù)侖應(yīng)力調(diào)整和加卸載作用受到廣泛關(guān)注（Toda et al，2008；單斌等，2013；程佳等，2018；盛書中等，2019）。本文采用ω2模型擬合震源譜獲得2009―2020年青川斷裂及鄰區(qū)3 285次1.5≤ML≤5.9地震的應(yīng)力降值Δσ，從10年尺度的應(yīng)力降動(dòng)態(tài)演化上分析汶川、蘆山和九寨溝地震前、后青川斷裂及鄰區(qū)的應(yīng)力調(diào)整過(guò)程和現(xiàn)今活動(dòng)性。為了有效解決數(shù)據(jù)取值范圍過(guò)大導(dǎo)致的配色問(wèn)題，統(tǒng)一將應(yīng)力降數(shù)據(jù)取對(duì)數(shù)后繪于圖8?；谥行〉卣鸬膽?yīng)力降具有隨震級(jí)增大而升高的特征（趙翠萍等，2011；吳微微等，2017；Wu et al，2020）。在分析比較不同地震的應(yīng)力降與區(qū)域應(yīng)力的關(guān)系時(shí)采用如下原則：震級(jí)相當(dāng)時(shí)，應(yīng)力降值相對(duì)高的震源區(qū)應(yīng)力水平高；應(yīng)力降值相近時(shí)，震級(jí)較低的震源區(qū)應(yīng)力水平高。為了更直觀地描述應(yīng)力較高的區(qū)域，取各時(shí)間段中Δσ≥1.0 MPa（即lgΔσ≥0 MPa）的應(yīng)力降繪于圖9，約定0 MPa≤lgΔσ≤0.8 MPa的應(yīng)力降為“中等應(yīng)力降”，lgΔσ＞0.8 MPa的應(yīng)力降為“高應(yīng)力降”。地震事件的應(yīng)力降空間位置采用經(jīng)TomoDD方法（Zhang，Thurber，2003，2006）重新定位后的結(jié)果。

圖8顯示，在龍門山斷裂帶北東段上，應(yīng)力降Δσ≥1.0 MPa的地震事件主要分布在青川斷裂（F1）、北川—映秀斷裂（F2）北東段和茶壩—林庵寺斷裂（F3）之間的區(qū)域，與汶川余震序列北段的優(yōu)勢(shì)分布比較一致。另有一些Δσ≥1.0 MPa的地震事件分布在虎牙斷裂（F5）附近，而灌縣—江油斷裂（F4）東南側(cè)靠近四川盆地的區(qū)域幾乎沒(méi)有Δσ≥1.0 MPa的地震事件分布。

在時(shí)間上，區(qū)域中小地震應(yīng)力降表現(xiàn)出與強(qiáng)主震較為相關(guān)的階段性變化特征（圖9）。階段一為2009—2012年，為汶川8.0級(jí)地震發(fā)生后到蘆山7.0級(jí)地震發(fā)生前的應(yīng)力變化過(guò)程。期間青川斷裂（F1）與茶壩—林庵寺斷裂（F3）之間集中了較多的高應(yīng)力降地震事件。汶川地震前、后地應(yīng)力數(shù)據(jù)也有類似結(jié)果。汶川地震前，青川測(cè)點(diǎn)的水平最大主應(yīng)力為21～22 MPa（郭啟良等，2009），震后復(fù)測(cè)結(jié)果仍高達(dá)15～16 MPa（郭啟良等，2009），遠(yuǎn)高于龍門山斷裂帶南段寶興測(cè)點(diǎn)的應(yīng)力9.8 MPa（吳滿路等，2010）。2012年開始，隨著汶川余震活動(dòng)的持續(xù)衰減，青川斷裂及鄰區(qū)高應(yīng)力降事件幾乎消失、中等應(yīng)力降地震事件明顯減少。

2013—2014年，在蘆山7.0級(jí)主震發(fā)生前（2013年1月至4月19日）出現(xiàn)高應(yīng)力降事件短暫消失的現(xiàn)象。4月20日蘆山7.0級(jí)地震后，高應(yīng)力降事件集中在距離7.0級(jí)主震較近的北川—映秀斷裂（F2）附近，區(qū)域中等應(yīng)力降事件明顯增多、空間分布廣泛，這種現(xiàn)象持續(xù)至2014年，期間在北川、青川分別發(fā)生多次高應(yīng)力降事件，北東段應(yīng)力水平呈兩端高、中間低的“啞鈴狀”分布格局。

從2015年開始，隨著地震序列的持續(xù)衰減，地震活動(dòng)逐漸減弱，幾乎沒(méi)有高應(yīng)力降事件發(fā)生。2017年九寨溝7.0級(jí)地震對(duì)研究區(qū)應(yīng)力水平?jīng)]有明顯影響，該地震距離青川斷裂（F1）約115 km，距離北川—映秀斷裂（F2）約150 km。庫(kù)侖應(yīng)力相關(guān)研究得出九寨溝地震僅在岷江斷裂北段、塔藏?cái)嗔押突⒀罃嗔涯隙卧斐奢^明顯的同震庫(kù)侖應(yīng)力變化（汪建軍，許才軍，2017）。GNSS觀測(cè)的同震位移數(shù)據(jù)表明九寨溝地震對(duì)龍門山斷裂帶影響非常?。ㄍ蹰啽龋?018）與本文的計(jì)算結(jié)果一致。

十年尺度的應(yīng)力降動(dòng)態(tài)演化表明，強(qiáng)震前、后震源區(qū)及鄰區(qū)的應(yīng)力水平存在明顯調(diào)整過(guò)程：余震密集期的地震，其應(yīng)力降普遍較主震前的地震更高；余震衰減期的地震，其應(yīng)力降值則相對(duì)較低。總的來(lái)說(shuō)，在強(qiáng)主震發(fā)生后，余震序列的應(yīng)力降在一段時(shí)間內(nèi)呈現(xiàn)出升高—回落的變化過(guò)程。

3.3 速度結(jié)構(gòu)特征

對(duì)汶川早期余震的研究表明，應(yīng)力積累水平的差異和龍門山斷裂帶不同段落間介質(zhì)物性的差異可能在汶川地震的孕育與控制破裂擴(kuò)展的過(guò)程中起著重要作用（易桂喜等，2011）。為了研究青川斷裂及鄰區(qū)的最新波速結(jié)構(gòu)特征，本文收集了2009—2020年ML≥1.5的觀測(cè)資料，采用TomoDD雙差層析成像方法（Zhang，Thurber，2003，2006）反演青川斷裂及鄰區(qū)（31.5°～33.5°N，104.0°～106.5°E）的三維P波速度結(jié)構(gòu)，并與汶川地震序列早期余震區(qū)波速特征進(jìn)行比較。反演過(guò)程中，挑選出至少有4個(gè)以上臺(tái)站記錄的地震事件共7 897個(gè)，采用規(guī)則網(wǎng)格節(jié)點(diǎn)構(gòu)建初始模型，??? 其中水平網(wǎng)格為0.2°×0.2°，垂直網(wǎng)格節(jié)點(diǎn)為0、5、10、15和20 km。正演走時(shí)計(jì)算選用偽彎曲射線追蹤技術(shù)，反演求解使用阻尼最小二乘分解算法（LSQR）（Paige，Saunders，1982），阻尼因子300，平滑權(quán)重40。依據(jù)前人研究成果（趙珠等，1997；鄧文澤等，2014）建立P波初始速度模型，波速比為1.71。定位過(guò)程中共拾取P波震相對(duì)99萬(wàn)個(gè)、S波震相對(duì)97萬(wàn)個(gè)、事件對(duì)14萬(wàn)個(gè)，最終獲得定位地震數(shù)7 858個(gè)。經(jīng)檢測(cè)板測(cè)試，模型中心大部分區(qū)域能夠分辨0.2°左右的速度異常。

遙感影像解譯顯示（謝小平等，2019），青川斷裂兩側(cè)的斷塊高程差異最為明顯、地形起伏最大，茶壩—林庵寺斷裂及北川—映秀斷裂北東段兩側(cè)斷塊高程差異次之，廣元—江油斷裂兩側(cè)斷塊的高程差異最小、地形起伏程度顯著降低。在淺部上地殼深度范圍內(nèi)，P波速度結(jié)構(gòu)圖像也表現(xiàn)出與上述斷裂的不同地形地貌、地質(zhì)構(gòu)造和巖相巖體特征相關(guān)的橫向不均勻性（圖10）。

在5 km深度層上（圖10a），青川斷裂（F1）西南段、茶壩—林庵寺斷裂（F3）西南段、北川—映秀斷裂（F2）北東段均為高速異常區(qū)。地質(zhì)資料顯示，青川斷裂尾端與北川—映秀斷裂之間存在地貌上的隆升區(qū)（樊春等，2008），P波速度的高速異常反映了平武、北川之間存在局部隆升和堅(jiān)硬介質(zhì)。青川斷裂的東側(cè)分布著廣元低速異常體和寧強(qiáng)高速異常體，兩個(gè)異常體在汶川地震序列余震階段的研究中也有體現(xiàn)（吳建平等，2009）。在10 km深度上（圖10b），高速異常體主要沿虎牙斷裂、龍門山斷裂帶北段一線分布，反映在此深度上主要構(gòu)造帶具有高強(qiáng)度的特征；青川斷裂西側(cè)的高速異常范圍明顯縮小，整體表現(xiàn)為弱高速特征，局部地區(qū)出現(xiàn)低速異常體，形成尺度不同、高低速相間的分塊結(jié)構(gòu)；青川斷裂西南段和北川—映秀斷裂北東段的高速異常體集中在平武—北川一帶，茶壩—林庵寺斷裂除東北端外整體顯示為弱高速異常。在15 km深度上（圖10c），青川斷裂與北川—映秀斷裂之間的高速異常區(qū)逐漸縮小，寧強(qiáng)一帶的高速異常程度減弱。在20 km深度上（圖10d），青川斷裂與北川—映秀斷裂之間的高速異常區(qū)趨于消失或縮小，其它地區(qū)基本保持15 km深度的速度分布特征。

綜上所述，在中、上地殼深度范圍內(nèi)青川斷裂（F1）、北川—映秀斷裂（F2）北東段和茶壩—林庵寺斷裂（F3）之間的區(qū)域?yàn)楦咚佼惓^(qū)，與汶川地震余震早期階段的研究結(jié)果一致（郭飚等，2009；劉啟元等，2009），該區(qū)域既是汶川余震序列在龍門山斷裂帶北東段的優(yōu)勢(shì)分布區(qū)，也是震源應(yīng)力降Δσ≥1.0 MPa的地震事件的集中發(fā)生區(qū)，其高速特征反映了龍門山斷裂帶北東段的強(qiáng)度明顯高于南段，這種高速堅(jiān)硬巖體的發(fā)育更有利于應(yīng)變能的積累與集中釋放（劉啟元等，2009）。

4 討論

青川斷裂在地表跡線可循，具備明顯的多期活動(dòng)特征，說(shuō)明該構(gòu)造具有潛在地震發(fā)生的構(gòu)造背景與能力，這揭示了汶川地震序列早期余震區(qū)北段5、6級(jí)地震頻發(fā)本身所具有的構(gòu)造因素，另一方面也昭示潛在震級(jí)上限可能不止6級(jí)。同時(shí)，青川斷裂的線性構(gòu)造地貌十分發(fā)育，顯示出典型的走滑運(yùn)動(dòng)特性，出露的斷裂剖面也顯示青川斷裂以走滑運(yùn)動(dòng)為主，兼有一定的擠壓運(yùn)動(dòng)分量，與前人探槽開挖揭露的斷層張性運(yùn)動(dòng)特征（Sun，2015；Lin et al，2014，2016；Liang et al，2018）并不一致。震源機(jī)制解結(jié)果支持青川斷裂具有走滑兼擠壓的運(yùn)動(dòng)特征，12次MW≥4.0地震中有11次走滑型或逆斷型地震。沿北川—映秀斷裂（F2）、灌縣—江油斷裂（F4）和茶壩—林庵寺斷裂（F3）向東北延伸至青川斷裂（F1）北東段，逐漸出現(xiàn)較多的走滑型地震，顯示出區(qū)域應(yīng)力場(chǎng)由逆沖過(guò)渡到走滑的變化過(guò)程。地應(yīng)力研究同樣顯示，8～10 km和10～15 km深度范圍內(nèi)，上地殼應(yīng)力狀態(tài)在青川附近從帶逆沖分量的走滑型過(guò)渡為純走滑型，15～20 km深度范圍內(nèi)從逆沖型過(guò)渡到帶逆沖分量的走滑型（豐成君等，2018）。這與龍門山斷裂帶的運(yùn)動(dòng)學(xué)模式（張培震等，2008）是一致的，即巴顏喀拉塊體SE向運(yùn)移受剛性四川盆地阻擋，在龍門山斷裂帶南段形成疊瓦狀大型逆沖斷裂帶，造成地殼縮短（杜方等，2009），部分物質(zhì)運(yùn)動(dòng)余量沿NE-SW向龍門山斷裂帶NE向側(cè)移，逐漸形成由逆沖到走滑的動(dòng)力學(xué)過(guò)程。

從地震活動(dòng)的時(shí)空分布上看，已發(fā)生的中強(qiáng)地震主要集中在2個(gè)區(qū)域：一個(gè)位于青川斷裂南東側(cè)的北川—映秀斷裂、灌縣—江油斷裂和茶壩—林庵寺斷裂的交會(huì)地帶（圖7中A區(qū)）；另一個(gè)沿青川斷裂走向展布（圖7中B區(qū)），地震活動(dòng)優(yōu)勢(shì)分布與構(gòu)造線的NE向延伸方向基本一致（圖1）。上述斷裂呈近平行分布，多為走滑斷裂或兼有走滑分量，一些研究表明發(fā)生在平行分布的走滑型構(gòu)造中的中強(qiáng)地震往往有較強(qiáng)觸發(fā)作用（龍鋒等，2021），這可能是該區(qū)域中強(qiáng)地震頻發(fā)的原因之一。

強(qiáng)震前后應(yīng)力降結(jié)果表明，汶川8.0級(jí)、蘆山7.0級(jí)2次強(qiáng)震改變了青川斷裂附近的應(yīng)力場(chǎng)，但九寨溝7.0級(jí)地震對(duì)該區(qū)域幾乎沒(méi)有影響。汶川、蘆山地震的震后應(yīng)力釋放過(guò)程主要沿構(gòu)造帶走向發(fā)生，高應(yīng)力降地震事件集中分布在主震或強(qiáng)余震震中附近，2014年魯?shù)?.5級(jí)和景谷6.6級(jí)地震有、后也有類似現(xiàn)象（周少輝，蔣海昆，2017）。蘆山7.0級(jí)主震發(fā)生前（2013年1月至4月19日），研究區(qū)出現(xiàn)高應(yīng)力降地震事件短暫消失的現(xiàn)象；在蘆山地震早期余震密集階段，高應(yīng)力降地震事件多，且集中在余震密集區(qū)；在地震序列衰減階段，地震活動(dòng)明顯減少，幾乎沒(méi)有高應(yīng)力降地震事件發(fā)生。盡管研究樣本有限，但仍然可以看出應(yīng)力降的時(shí)空變化對(duì)后續(xù)地震活動(dòng)趨勢(shì)有一定指示意義，同等震級(jí)條件下余震密集期的中小地震應(yīng)力降值往往高于余震衰減期，而余震區(qū)中高應(yīng)力降事件的明顯減少表明區(qū)域的應(yīng)力水平逐漸降低，這也許意味著地震序列開始進(jìn)入衰減期，近期發(fā)生強(qiáng)震的可能性較小。

體波層析成像的結(jié)果顯示，青川斷裂及鄰區(qū)的地殼淺層高低速異常表現(xiàn)出以青川斷裂帶為分界線的空間分布特征，高速異常體與局部隆升和堅(jiān)硬介質(zhì)比較相關(guān)。隨著深度增加，局部地區(qū)出現(xiàn)低速異常體，形成尺度不一的分塊結(jié)構(gòu)。其中，青川斷裂表現(xiàn)為高速異常體，體現(xiàn)了該斷裂具有較高強(qiáng)度。15 km深度處的青川斷裂NE端表現(xiàn)為強(qiáng)度較弱的低波速區(qū)，這個(gè)深度是汶川地震序列在北段的深度下限，較弱的結(jié)構(gòu)阻止了破裂的持續(xù)擴(kuò)展。

5 結(jié)論

本文以青川斷裂為研究對(duì)象，通過(guò)詳細(xì)的地震地質(zhì)調(diào)查，采用無(wú)人機(jī)攝影技術(shù)對(duì)斷裂沿線的構(gòu)造變形地貌進(jìn)行測(cè)繪，獲得斷裂活動(dòng)的構(gòu)造變形信息；通過(guò)震源機(jī)制解反映的地震錯(cuò)動(dòng)特征，討論青川斷裂的構(gòu)造變形特征；采用應(yīng)力降繪制汶川余震區(qū)十年尺度的應(yīng)力降動(dòng)態(tài)分布圖像，分析汶川、蘆山和九寨溝地震對(duì)青川斷裂及鄰區(qū)應(yīng)力水平的影響；通過(guò)對(duì)青川斷裂及鄰區(qū)速度結(jié)構(gòu)進(jìn)行反演，探討深部孕震背景；最后綜合地質(zhì)學(xué)和地震學(xué)研究結(jié)果分析青川斷裂及鄰區(qū)的現(xiàn)今活動(dòng)性。研究表明：

（1）青川斷裂在地貌上呈現(xiàn)出走滑特征，斷層沿線出露的剖面主要以片理化帶、斷層泥為主，局部見(jiàn)擠壓作用形成的構(gòu)造角礫。

（2）青川斷裂NE段為2008年汶川8.0級(jí)地震余震北延段的密集區(qū)，余震錯(cuò)動(dòng)性質(zhì)既有走滑型，也有走滑兼逆斷型。

（3）震源機(jī)制解和地質(zhì)地貌調(diào)查的結(jié)果均反映青川斷裂的運(yùn)動(dòng)性質(zhì)是走滑兼擠壓，而并非前人探槽所揭露的斷層張性活動(dòng)特征。

（4）汶川、蘆山兩次強(qiáng)震發(fā)生前后，青川斷裂及鄰區(qū)的應(yīng)力水平在一段時(shí)間內(nèi)呈現(xiàn)出升高—回落的調(diào)整過(guò)程。

（5）體波層析成像的結(jié)果顯示青川斷裂及鄰區(qū)的地殼淺層高低速異常表現(xiàn)以龍門山斷裂帶為分界線的空間分布特征。

（6）雖然青川斷裂具有潛在地震發(fā)生的構(gòu)造背景與能力，但目前該斷裂及鄰區(qū)的地震活動(dòng)正處于強(qiáng)震序列衰減階段，M4.0以上地震活動(dòng)呈逐漸減弱的趨勢(shì)，高應(yīng)力降地震事件明顯減少，區(qū)域整體應(yīng)力水平較低。

參考文獻(xiàn)：

陳國(guó)光，計(jì)鳳桔，周榮軍，等.2007.龍門山斷裂帶晚第四紀(jì)活動(dòng)性分段的初步研究［J］.地震地質(zhì)，29（3）：657-673.

程佳，姚生海，劉杰，等.2018.2017年九寨溝地震所受歷史地震黏彈性庫(kù)侖應(yīng)力作用及其后續(xù)對(duì)周邊斷層地震危險(xiǎn)性的影響［J］.地球物理學(xué)報(bào)，61（5）：2133-2151.

鄧起東，陳社發(fā)，趙小麟.1994.龍門山及其鄰區(qū)的構(gòu)造和地震活動(dòng)及動(dòng)力學(xué)［J］.地震地質(zhì)，16（4）：389-403.

鄧文澤，陳九輝，郭飚，等.2014.龍門山斷裂帶精細(xì)速度結(jié)構(gòu)的雙差層析成像研究［J］.地球物理學(xué)報(bào)，57（4）：1101-1110.

杜方，聞學(xué)澤，張培震，等.2009.2008年汶川8.0級(jí)地震前橫跨龍門山斷裂帶的震間形變［J］.地球物理學(xué)報(bào)，52（11）：2729-2738.

杜建軍，馬寅生，尹成明，等.2013.龍門山北部陜甘川交界三角構(gòu)造區(qū)斷裂活動(dòng)特征研究［J］.地震學(xué)報(bào)，35（4）：520-533.

樊春，王二七，王剛，等.2008.龍門山斷裂帶北段晚新近紀(jì)以來(lái)的右行走滑運(yùn)動(dòng)及其構(gòu)造變換研究——以青川斷裂為例［J］.地質(zhì)科學(xué)，43（3）：417-433.

豐成君，戚幫申，張鵬，等.2018.汶川MS8.0地震后龍門山斷裂帶地殼應(yīng)力場(chǎng)及其構(gòu)造意義［J］.地質(zhì)力學(xué)學(xué)報(bào)，24（4）：439-451.

郭飚，劉啟元，陳九輝，等.2009.川西龍門山及鄰區(qū)地殼上地幔遠(yuǎn)震P波層析成像［J］.地球物理學(xué)報(bào)，52（2）：346-355.

郭啟良，王成虎，馬洪生，等.2009.汶川MS8.0級(jí)大震前后的水壓致裂原地應(yīng)力測(cè)量［J］.地球物理學(xué)報(bào)，52（5）：1395-1401.

胡幸平，俞春泉，陶開，等.2008.利用P波初動(dòng)資料求解汶川地震及其強(qiáng)余震震源機(jī)制解［J］.地球物理學(xué)報(bào)，51（6）：1711-1718.

李傳友，宋方敏，冉勇康.2004.龍門山斷裂帶北段晚第四紀(jì)活動(dòng)性討論［J］.地震地質(zhì)，26（2）：248-258.

劉啟元，李昱，陳九輝，等.2009.汶川MS8.0地震：地殼上地幔S波速度結(jié)構(gòu)的初步研究［J］.地球物理學(xué)報(bào)，52（2）：309-319.

龍鋒，倪四道，聞學(xué)澤.2011.用遠(yuǎn)震接收函數(shù)研究龍門山斷裂帶與其鄰區(qū)的莫霍面深度及波速比分布［J］.地球?qū)W報(bào)，32（4）：438-446.

龍鋒，祁玉萍，易桂喜，等.2021.2021年5月21日云南漾濞MS6.4地震序列重新定位與發(fā)震構(gòu)造分析［J］.地球物理學(xué)報(bào)，64（8）：2631-2646.

單斌，熊熊，鄭勇，等.2009.2008年5月12日MW7.9汶川地震導(dǎo)致的周邊斷層應(yīng)力變化［J］.中國(guó)科學(xué)：地球科學(xué)，39（5）：537-545.

單斌，熊熊，鄭勇，等.2013.2013年蘆山地震導(dǎo)致的周邊斷層應(yīng)力變化及其與2008年汶川地震的關(guān)系［J］.中國(guó)科學(xué)：地球科學(xué)，43：1002-1009.

盛書中，萬(wàn)永革，徐錫偉，等.2019.用大量地震震源機(jī)制解檢驗(yàn)2008年汶川地震對(duì)后續(xù)地震的觸發(fā)作用［J］.地球物理學(xué)報(bào)，62（12）：4588-4603.

石玉濤，高原，張永久，等.2013.松潘—甘孜地塊東部、川滇地塊北部與四川盆地西部的地殼剪切波分裂［J］.地球物理學(xué)報(bào)，56（2）：481-494.

孫浩越.2015.青川斷裂晚第四紀(jì)活動(dòng)性及其對(duì)區(qū)域構(gòu)造運(yùn)動(dòng)模式的約束［D］. 北京：中國(guó)地震局地質(zhì)研究所，59.

汪建軍，許才軍.2017.2017年MW6.5九寨溝地震激發(fā)的同震庫(kù)侖應(yīng)力變化及其對(duì)周邊斷層的影響［J］.地球物理學(xué)報(bào)，60（11）：4398-4420.

王閱兵，甘衛(wèi)軍，陳為濤，等.2018.GNSS觀測(cè)的九寨溝7.0級(jí)地震同震位移初步結(jié)果［J］.地球物理學(xué)報(bào)，61（1）：161-170.

吳建平，黃媛，張?zhí)熘?，?2009.汶川MS8.0級(jí)地震余震分布及周邊區(qū)域P波三維速度結(jié)構(gòu)研究［J］.地球物理學(xué)報(bào)，52（2）：320-328.

吳滿路，張?jiān)罉颍未煌?，?2010.汶川地震后沿龍門山斷裂帶原地應(yīng)力測(cè)量初步結(jié)果［J］.地質(zhì)學(xué)報(bào)，84（9）：1291-1299.

吳微微，魏婭玲，龍鋒，等.2020.2017年8月8日四川九寨溝M7.0地震及其余震序列的震源參數(shù)［J］.地震地質(zhì)，42（2）：492-512.

吳微微，吳朋，魏婭玲，等.2017.川滇活動(dòng)塊體中-北部主要活動(dòng)斷裂帶現(xiàn)今應(yīng)力狀態(tài)的分區(qū)特征［J］.地球物理學(xué)報(bào)，60（5）：1735-1745.

謝小平，白毛偉，陳芝聰，等.2019.龍門山斷裂帶北東段活動(dòng)斷裂的遙感影像解譯及構(gòu)造活動(dòng)性分析［J］.國(guó)土資源遙感，31（1）：237-246.

徐錫偉，韓竹軍，楊曉平，等.2016.中國(guó)及鄰近地區(qū)地震構(gòu)造圖［M］.北京：地震出版社.

徐錫偉，聞學(xué)澤，葉建青，等.2008.汶川MS8.0地震地表破裂帶及其發(fā)震構(gòu)造［J］.地震地質(zhì)，30（3）：597-629.

易桂喜，龍鋒，張致偉.2012.汶川MS8.0地震余震震源機(jī)制時(shí)空分布特征［J］.地球物理學(xué)報(bào)，55（4）：1213-1227.

易桂喜，聞學(xué)澤，辛華，等.2011.2008年汶川MS8.0地震前龍門山—岷山構(gòu)造帶的地震活動(dòng)性參數(shù)與地震視應(yīng)力分布［J］.地球物理學(xué)報(bào)，54（6）：1490-1500.

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

趙翠萍，陳章立，華衛(wèi)，等.2011.中國(guó)大陸主要地震活動(dòng)區(qū)中小地震震源參數(shù)研究［J］.地球物理學(xué)報(bào)，54（6）：1478-1489.

趙珠，范軍，鄭斯華，等.1997.龍門山斷裂帶地殼速度結(jié)構(gòu)和震源位置的精確修訂［J］.地震學(xué)報(bào)，19（6）：615-622.

鄭勇，馬宏生，呂堅(jiān)，等.2009.汶川地震強(qiáng)余震（MS≥5.6）的震源機(jī)制解及其與發(fā)震構(gòu)造的關(guān)系［J］.中國(guó)科學(xué)：地球科學(xué)，39（4）：413-426.

周少輝，蔣海昆.2017.景谷6.6級(jí)、魯?shù)?.5級(jí)地震序列應(yīng)力降變化對(duì)比研究［J］.中國(guó)地震，33（1）：23-37.

Feng G C，Jónsson S，Klinger Y.2017.Which fault segments ruptured in the 2008 Wenchuan earthquake and which did not？ New evidence from near-fault 3D surface displacements derived from SAR image offsets［J］.Bulletin of the Seismological Society of America，107（3）：1185-1200.

Gao M X，Zeilinger G，Xu X W，et al.2016.Active tectonics evaluation from geomorphic indices for the central and the southern Longmenshan range on the Eastern Tibetan Plateau，China［J］.Tectonics，35（8）：1812-1826.

Liang M J，Ran Y K，Wang H，et al.2018.A possible tectonic response between the Qingchuan Fault and the Beichuan-Yingxiu Fault of the Longmen Shan Fault Zone？ evidence from geologic observations by Paleoseismic trenching and radiocarbon dating［J］.Tectonics，37：4086-4096.

Lin A M，Chen P，Satsukawa T，et al.2016.Millennium recurrence interval of morphogenic earthquakes on the Qingchuan fault，northeastern segment of the Longmen Shan Thrust Belt，China［J］.Journal of Seismology，20（2）：535-553.

Lin A M，Gang R，Bing Y.2012.Field evidence of rupture of the Qingchuan Fault during the 2008 MW7.9 Wenchuan earthquake，northeastern segment of the Longmen Shan Thrust Belt，China［J］.Tectonophysics，522-523：243-252.

Lin A M，Gang R，Bing Y.2014.Structural analysis of the right-lateral strike-slip Qingchuan fault，northeastern segment of the Longmen Shan thrust belt，central China［J］.Journal of Structural Geology，68：227-244.

Long F，Yi G X，Wen X Z.2012.Spatio-temporal variation of the stress field in the Wenchuan aftershock region［J］.Earthquake Science，25（S1）：517-526.

Paige C C，Saunders M A.1982.LSQR：An algorithm for sparse linear equations and sparse least squares［J］.ACM Transactions on Mathematical Software（TOMS），8（1）：43-71.

Pasyanos M E，Dreger D S，Romanowicz B.1996.Towards real-time estimation of regional moment tensors［J］.Bulletin of the Seismological Society of America，86（5）：1255-1269.

Ran Y K，Chen L C，Chen J，et al.2010.Paleoseismic evidence and repeat time of large earthquakes at three sites along the Longmenshan fault zone［J］.Tectonophysics，491（1-4）：141-153.

Ran Y K，Chen W S，Xu X W，et al.2013.Paleoseismic events and recurrence interval along the Beichuan-Yingxiu fault of Longmenshan fault zone，Yingxiu，Sichuan，China［J］.Tectonophysics，584：81-90.

Ran Y K，Chen W S，Xu X W，et al.2014.Late Quaternary paleoseismic behavior and rupture segmentation of the Yingxiu-Beichuan Fault along the Longmenshan fault zone，China［J］.Tectonics，33（11）：2218-2232.

Shen Z K，Sun J，Zhang P，et al.2009.Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake［J］.Nature Geoscience，2（10）：718-724.

Sun H Y.2015.Late quaternary activity of the Qingchuan fault：Implications for the tectonic movement mechanism in reginal area［D］.Beijing：Institute of Geology，China Earthquake Adminiseration.

Tan Y，Zhu L P，Helmberger D V，et al.2006.Locating and modeling regional earthquakes with two stations［J］.Journal of Geophysical Research：Solid Earth，111（B1）：B01306.

Toda S J，Lin J，Meghraoui M，et al.2008.12 May 2008 M=7.9 Wenchuan，China，earthquake calculated to increase failure stress and seismicity rate on three major fault systems［J］.Geophysical Research Letters，35（17）：L17305.

Wang Y Z，Wang E N，Shen Z K，et al.2008.GPS-constrained inversion of present-day slip rates along major faults of the Sichuan-Yunnan region，China［J］.Science in China Series D：Earth Sciences，51（9）：1267-1283.

Wu W W，Long F，Liang M J，et al.2020.Spatial and Temporal Variations in Earthquake Stress Drops between the 2008 Wenchuan and 2013 Lushan Earthquakes［J］.Acta Geologica Sinica（English edition），94（5）：1635-1650.

Xu X W，Wen X Z，Yu G H，et al.2009.Coseismic reverse-and oblique-slip surface faulting generated by the 2008 MW 7.9 Wenchuan earthquake，China［J］.Geology，37（6）：515-518.

Zhang H J，Thurber C H.2003.Double-difference tomography：The method and its application to the Hayward fault，California［J］.Bulletin of the Seismological Society of America，93（5）：1875-1889.

Zhang H J，Thurber C H.2006.Development and applications of double-difference seismic tomography［J］.Pure and Applied Geophysics，163（2）：373-403.

Zhao L S，Helmberger D V.1994.Source estimation from broadband regional seismograms［J］.Bulletin of the Seismological Society of America，84（1）：91-104.

Zhu L P，Helmberger D V.1996.Advancement in source estimation techniques using broadband regional seismograms［J］.Bulletin of the Seismological Society of America，86（5）：1634-1641.

Zoback M L.1992.First-and second-order patterns of stress in the lithosphere：The world stress map project［J］.Journal of Geophysical Research，97（B8）：11703-11728.

Analysis of the Present-day Activity of the Qingchuan Fault and Its Vicinity

WU Weiwei，LIANG Mingjian，LONG Feng，SU Jinrong，CHEN Xuefen

（Sichuan Earthquake Agency，Chengdu 610041，Sichuan，China）

Abstract

In recent years，earthquakes（M≥4.0）have been repeatedly occurring on the Qingchuan fault and in its adjacent region. Firstly，we conduct a geological investigation and carry out a surveying and mapping of landforms area by the unmanned aviation vehicle in the study in order to obtain the tectonic deformation of the Qingchuan fault. Secondly，we use the cut-and-paste method to invert the focal mechanisms of 32 earthquakes（MW≥4.0）occurring during May 17-December 31，2008 in the study area，and further study the tectonic features of the Qingchuan fault. Thirdly，we use the ω2 model to fit the source spectra to obtain the dynamic evolution of the stress drop in 10-year span in the study area. Then we analyze the stress variation before and after the 2008 Wenchuan earthquake，the 2013 Lushan earthquake，the 2017 Jiuzhaigou earthquake. Fourthly，we invert the P-wave velocity structure in the study area to discuss the tectonic settings for the earthquake genesis in the deep. Finally，on the basis of the above the geological and seismological results we analyze the resent-day activity of the Qingchuan fault. We find that：①The geological survey and the focal mechanism solution show that movement property of the Qingchuan fault is mainly strike-slip，with a certain compressive component. ②Before and after the 2008 Wenchuan earthquake and the 2013 Lushan earthquake，the stress in the study area experienced a rise-fall period. ③Spatially，the body-wave tomography shows that the high-speed anomaly and low-speed anomaly in the shallow Crust are divided by the Longmenshan fault. ④The Qingchuan fault is an earthquake-prone structure because of its tectonic settings. So far，the strong-earthquake activity is in the weakening period; the seismicity（M≥4.0）takes on a trend of gradual attenuation. High stress-drop events are significantly lessening and the overall stress level in the study area is relatively low.

Keywords：the Qingchuan fault and its vicinity；the Wenchuan MS8.0 earthquake；focal mechanism；stress drop；velocity structure；earthquake activity

收稿日期：2022-11-22.

基金項(xiàng)目：中國(guó)地震局地震科技星火計(jì)劃項(xiàng)目（XH202302）；震情跟蹤定向工作任務(wù)（2021020102）；國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2018YFC1504501-02）.

第一作者簡(jiǎn)介：吳微微（1982-），正高級(jí)工程師，主要從事地震監(jiān)測(cè)、數(shù)字地震學(xué)應(yīng)用等方面研究.E-mail：673065247@qq.com.

通訊作者簡(jiǎn)介：蘇金蓉（1971-），正高級(jí)工程師，主要從事地震監(jiān)測(cè)預(yù)警、地震學(xué)應(yīng)用研究.E-mail：sujr0816@163.com.

吳微微，梁明劍，龍鋒，等.2023.青川斷裂及鄰區(qū)現(xiàn)今地震活動(dòng)性研究［J］.地震研究，46（2）：188-203，doi：10.20015/j.cnki.ISSN1000-0666.2023.0035.