基于線性密集臺(tái)陣的2023年甘肅積石山MS6.2地震斷層低速帶研究

吳柏辰 秦滿忠 郭曉 王亞紅 孫點(diǎn)峰 溫淑琳

摘要：2023年12月18日，甘肅省積石山縣發(fā)生MS6.2逆沖型地震。該地震發(fā)生在祁連活動(dòng)地塊拉脊山北緣斷裂帶，在斷裂兩側(cè)曾發(fā)生過(guò)20余次5級(jí)左右中強(qiáng)地震，地震活動(dòng)頻次高，因此研究拉脊山斷裂帶特征具有重要意義。根據(jù)拉脊山斷裂走向和此次地震余震分布，在劉集鄉(xiāng)布設(shè)一條線性密集臺(tái)陣，利用不同方位的余震事件計(jì)算體波的走時(shí)延時(shí)和放大效應(yīng)。結(jié)果表明，P波和S波分別存在6個(gè)采樣點(diǎn)（約0.012 s）和15個(gè)采樣點(diǎn)（約0.03 s）的走時(shí)延時(shí)，推斷出臺(tái)陣西側(cè)存在一個(gè)約150 m寬近似垂直的斷層低速帶，且該結(jié)果與遠(yuǎn)震P波相似性系數(shù)矩陣結(jié)果一致。

關(guān)鍵詞：甘肅積石山MS6.2地震; 走時(shí)延時(shí); 體波放大效應(yīng); 線性密集臺(tái)陣; 斷層低速帶

中圖分類號(hào)： P315????? 文獻(xiàn)標(biāo)志碼：A?? 文章編號(hào)： 1000-0844（2024）03-0734-08

DOI：10.20000/j.1000-0844.20240306003

Low-velocity fault zone of the 2023 Jishishan Gansu MS6.2

earthquake based on linear dense array data

WU Baichen1，2， QIN Manzhong1，2， GUO Xiao1，2， WANG Yahong2， SUN Dianfeng2，WEN Shulin1，2

（1. Lanzhou Institute of Seismology， CEA， Lanzhou 730000， Gansu， China;

2. Gansu Lanzhou Geophysics Observation and Research Station， Lanzhou 730000， Gansu， China）

Abstract：?On December 18， 2023， an MS6.2 reverse-type earthquake struck Jishishan County， Gansu Province， along the northern margin of the Lajishan fault zone， Qilian active block. Given that over 20 M5.0 earthquakes have occurred on both sides of the Lajishan fault， investigating its characteristics is paramount. Thus， based on the strike of the Lajishan fault and its aftershock distribution， we deployed a linear dense array in Liuji Town. Considering aftershock events in different directions， we determined travel time delays and amplification effects of body waves. The results identified six （approximately 0.012 s） and 15 （approximately 0.03 s） sampling points of travel time delays for the P- and S-waves， respectively. This confirmed the presence of a nearly vertical low-velocity zone with a width of approximately 150 m along the west side of the array. The obtained results align with teleseismic P-wave similarity coefficient matrix calculations.

Keywords：Jishishan MS6.2 earthquake in Gansu; travel time delay; amplification effect of body wave; linear dense array; low-velocity zone of the fault

0 引言

地殼斷層是在地殼運(yùn)動(dòng)過(guò)程中普遍形成的一種地質(zhì)構(gòu)造現(xiàn)象。地震通常發(fā)生在斷層帶附近，是斷層的劇烈運(yùn)動(dòng)形式。因此，了解活動(dòng)斷層帶及周圍的破裂特征在地球物理學(xué)研究中有著重要的意義。地震時(shí)產(chǎn)生的強(qiáng)烈破裂、角礫巖作用、流體飽和度，以及斷層附近孔隙流體壓力的升高等因素都有助于斷層低速帶的形成［1］。斷層帶的產(chǎn)生是地震作用的結(jié)果，該區(qū)域主要由高度破碎的物質(zhì)、斷層角礫巖、黏土和巖石碎屑組成［2］。這些物質(zhì)的存在會(huì)顯著影響地震波的傳播路徑和速度，為進(jìn)行詳細(xì)的斷層結(jié)構(gòu)研究奠定了有利條件［3］。

學(xué)術(shù)界對(duì)斷層的研究?jī)?nèi)容比較廣泛，包括斷層的地面幾何分布、地形變化觀測(cè)研究及斷層帶的內(nèi)部形態(tài)結(jié)構(gòu)、地球物理學(xué)特征和物性參數(shù)等，通常利用重力、電磁測(cè)量、走時(shí)層析成像和基于地震波形的反演建模等多種地球物理方法研究斷裂帶的性質(zhì)［4］。與其他的地球物理學(xué)方法相比，密集地震臺(tái)陣方法具有價(jià)格成本低，布設(shè)方便、環(huán)保［5］，且檢測(cè)微震相的能力強(qiáng)、敏感度高［6］ 等優(yōu)勢(shì)，特別是利用測(cè)線的方式布設(shè)臺(tái)站［7］。該方法已經(jīng)成為地殼斷層結(jié)構(gòu)探測(cè)研究中的重要手段，密集臺(tái)陣可以對(duì)近震及遠(yuǎn)震事件進(jìn)行高質(zhì)量觀測(cè)，可以提高定位精度和區(qū)域地震監(jiān)測(cè)能力［8］。

近年來(lái)，許多地震學(xué)家基于密集臺(tái)陣方法研究地球內(nèi)部不同尺度的速度結(jié)構(gòu)［9-10］，了解斷層帶的物性信息，主要包括近震或遠(yuǎn)震事件的體波到時(shí)延時(shí)方法［11］，近震事件的體波放大效應(yīng)［12-13］并反演斷層帶圍陷波［14-16］，利用斷層帶首波計(jì)算斷層內(nèi)外速度差［17-18］，以及利用背景噪聲計(jì)算放大效應(yīng)［19］。在斷層帶可以產(chǎn)生更長(zhǎng)周期、更大振幅的體波和尾波時(shí)，低速帶與高速圍巖的區(qū)域?qū)?huì)產(chǎn)生更明顯的波形放大效應(yīng)［20］。利用體波在低速介質(zhì)中產(chǎn)生的放大效應(yīng)可以有效地研究斷層帶結(jié)構(gòu)，特別是在確定斷層帶寬度方面具有明顯的優(yōu)勢(shì)。Li Y G等［21］利用體波在低速介質(zhì)下的放大效應(yīng)特征有效約束了美國(guó)加州Landers斷層寬度約200 m;在昆侖山口西MS8.1地震斷層區(qū)，Li S L等［22］利用體波在低速介質(zhì)下的放大效應(yīng)特征有效約束了300 m寬的斷層區(qū)域。野外科考發(fā)現(xiàn)，甘肅積石山MS6.2地震并沒(méi)有產(chǎn)生明顯的地表破裂，但在不同區(qū)域觀察到了較明顯的同震地表現(xiàn)象（地裂縫）［23］。項(xiàng)目組根據(jù)拉脊山斷裂帶走向、余震分布及地震科考同震地表現(xiàn)象等信息，結(jié)合野外實(shí)地勘察，在劉集鄉(xiāng)布設(shè)了一條線性密集地震臺(tái)陣。本研究將基于布設(shè)的線性密集臺(tái)陣觀測(cè)到的近震事件及遠(yuǎn)震事件，計(jì)算體波震相的走時(shí)延時(shí)及放大效應(yīng)，并結(jié)合遠(yuǎn)震P波相似性系數(shù)矩陣，推斷斷層低速帶的具體位置及寬度。

1 構(gòu)造概況及臺(tái)陣布設(shè)

甘肅積石山地震發(fā)生在青藏高原東北緣柴達(dá)木—祁連塊體拉脊山斷裂帶南段附近［24］。拉脊山斷裂帶由拉脊山北緣斷裂和拉脊山南緣斷裂兩條向NE方向凸出的弧形逆沖構(gòu)造帶構(gòu)成［25］，運(yùn)動(dòng)性質(zhì)以擠壓逆沖為主，兼左旋走滑，是調(diào)節(jié)NNW向熱水—日月山右旋走滑斷裂帶與NWW向西秦嶺北緣左旋走滑斷裂帶之間的構(gòu)造轉(zhuǎn)換帶。晚第四紀(jì)以來(lái)該斷裂帶有新活動(dòng)，不斷向盆地內(nèi)部擠壓逆沖產(chǎn)生局部應(yīng)力積累，使得該斷裂未來(lái)存在發(fā)生中強(qiáng)地震的構(gòu)造條件。

甘肅積石山發(fā)生MS6.2地震后，許多地震專家的科考結(jié)果表明，該地震沒(méi)有明顯的地表破裂。本文根據(jù)拉脊山斷裂帶走向、余震分布及地震科考同震地表現(xiàn)象等信息，結(jié)合野外實(shí)地勘察［26］，在劉集鄉(xiāng)布設(shè)了一條線性密集臺(tái)陣。該臺(tái)陣位于拉脊山北緣斷裂與拉脊山南緣斷裂之間［圖1（a）］。各臺(tái)站相對(duì)位置及其高程見(jiàn)圖2。

2 數(shù)據(jù)采集與處理

布設(shè)的線性密集臺(tái)陣包括26個(gè)便攜式一體化地震儀，總長(zhǎng)度約700 m［圖1（b）］，采樣率為500 Hz。觀測(cè)時(shí)間從2023年12月28日開(kāi)始，至2024年1月18日結(jié)束。地震儀記錄了地表速度，并觀測(cè)到0.1～100 Hz相對(duì)平緩的波形信號(hào)（圖3）。

臺(tái)陣在布設(shè)期間記錄到了大量的近震與遠(yuǎn)震事件［圖1（a）］。其中，近震序列整體呈NW向展布，大部分近震集中在主震西側(cè)，并且集中分布在拉脊山北緣斷裂和拉脊山南緣斷裂之間。在深度方向上，根據(jù)左可楨等［27］定位結(jié)果，近震主要集中在5～15 km，且絕大部分事件震源深度比主震淺，精定位結(jié)果揭示了主震的破裂特征。近震在頻譜圖中表現(xiàn)為頻率在0.4～20 Hz間，表明其有更高的能量。與近震事件的波形相比，遠(yuǎn)震產(chǎn)生的地震波則具有更持久的低頻能量，通常低于2 Hz（圖3）。這種能量和頻率的差異可以用來(lái)研究特定區(qū)域的速度變化?？紤]到本研究所使用儀器的固有頻率為0.2 Hz，本文將遠(yuǎn)震數(shù)據(jù)引入分析過(guò)程，旨在補(bǔ)充和驗(yàn)證已收集的近震數(shù)據(jù)，從而提供更加全面和準(zhǔn)確的地震活動(dòng)信息，使得本次研究的地震資料既有高頻數(shù)據(jù)（近震），也有低頻數(shù)據(jù)（遠(yuǎn)震）。

根據(jù)美國(guó)地質(zhì)調(diào)查局（United States Geological Survey，USGS）地震目錄，本文共選取了震中距大于15°，M＞5，且信噪比（Signal Noise Ratio，SNR）大于3的9個(gè)遠(yuǎn)震事件進(jìn)行遠(yuǎn)震P波相似性系數(shù)矩陣計(jì)算（表1）。

本文利用中國(guó)地震局地震臺(tái)網(wǎng)目錄，使用ObsPy軟件包［28］中的ar_pick函數(shù)進(jìn)行體波（P波、S波）到時(shí)檢測(cè)，然后針對(duì)所有臺(tái)站波形人工核對(duì)數(shù)據(jù)質(zhì)量。共觀測(cè)到165個(gè)余震事件，并人工篩選出100個(gè)高信噪比的地震事件進(jìn)行近震事件體波延時(shí)分析、放大效應(yīng)等研究。

在利用近震確定斷層位置時(shí)，不僅可以利用近震體波放大效應(yīng)的觀測(cè)數(shù)據(jù)，還可以結(jié)合近震與遠(yuǎn)震的體波走時(shí)延時(shí)分析法，從而有效地確定斷層的位置［29］。本文參考Yang等［30］數(shù)據(jù)處理流程，并根據(jù)實(shí)際觀測(cè)數(shù)據(jù)做相關(guān)改進(jìn)：（1）遠(yuǎn)震事件使用0.2～2 Hz帶通濾波進(jìn)行計(jì)算，并對(duì)數(shù)據(jù)進(jìn)行去均、滅尖處理;近震事件使用1～10 Hz帶通濾波進(jìn)行處理，并對(duì)數(shù)據(jù)進(jìn)行去均，滅尖處理。（2）考慮到地形和高程對(duì)P波到時(shí)的影響，使用地表速度2.7 km/s對(duì)臺(tái)站高程進(jìn)行校正［31］，從而消除地表高程對(duì)P波到時(shí)的影響。（3） 在2 s長(zhǎng)的窗口范圍內(nèi)拾取精確的近震事件P波到時(shí)，從而盡可能避免在拾取P波到時(shí)的過(guò)程中產(chǎn)生誤差；對(duì)于遠(yuǎn)震事件，本文對(duì)每個(gè)遠(yuǎn)震事件根據(jù)P波到時(shí)截取了10 s的觀測(cè)波形，并計(jì)算波形間的相關(guān)系數(shù)，以評(píng)估斷裂帶兩側(cè)的波形相似性。

在近震事件中觀測(cè)到了部分臺(tái)站記錄的體波走時(shí)延時(shí)，這與體波信號(hào)穿過(guò)斷層低速帶有關(guān)。經(jīng)過(guò)低速帶的體波有長(zhǎng)周期和大振幅的特征，特別是體波的尾波會(huì)產(chǎn)生放大效應(yīng)，該特征為尋找斷層低速帶提供了理論依據(jù)［32-33］。

本文以觀測(cè)到的一個(gè)質(zhì)量較好的近震事件為例（圖4），即密集臺(tái)陣于2023年12月29日15：35：13（UTC）觀測(cè)到的震級(jí)為M1.1的近震事件，并對(duì)各臺(tái)站觀測(cè)波形做了1～10 Hz帶通濾波處理，圖中紅星代表自動(dòng)拾取的P波到時(shí);分別對(duì)各臺(tái)站觀測(cè)到的波形做歸一化處理，并計(jì)算其地表峰值速度（Peak Ground Velocity，PGV）和均方根值（Root Mean Square，RMS），黑色曲線為PGV與RMS的乘積［圖4（b）］。該曲線代表地表介質(zhì)對(duì)地震波的放大效應(yīng)，本文在密集臺(tái)陣西側(cè)100～250 m的區(qū)域觀測(cè)到了明顯的放大效應(yīng)。同時(shí)，在該區(qū)域也發(fā)現(xiàn)了明顯的P波走時(shí)延時(shí)，而走時(shí)延時(shí)與地表高程并沒(méi)有一致性對(duì)應(yīng)特征，說(shuō)明延時(shí)現(xiàn)象與臺(tái)陣下方區(qū)域介質(zhì)速度較低有關(guān)。

3 方法與分析

3.1 遠(yuǎn)震相關(guān)系數(shù)矩陣計(jì)算

遠(yuǎn)震P波在穿過(guò)斷層低速帶時(shí)，位于低速帶附近的臺(tái)站觀測(cè)到的波形具有一定相似性?；赯hang等［34］的遠(yuǎn)震相關(guān)系數(shù)矩陣計(jì)算方法，選取了表1中9個(gè)分布在不同方位的遠(yuǎn)震事件，運(yùn)用窗長(zhǎng)為2 s的P波計(jì)算各臺(tái)站之間的波形相似性，并構(gòu)建了相似性矩陣（圖5），其中每個(gè)顏色塊代表臺(tái)站對(duì)之間同一遠(yuǎn)震事件波形的最大相關(guān)系數(shù)值，顏色塊大小與臺(tái)站間距成正比。在相似性矩陣中，對(duì)角線代表每個(gè)臺(tái)站自相關(guān)的結(jié)果（值均為1）；在非對(duì)角線處通常相關(guān)性系數(shù)會(huì)隨著臺(tái)站間距的增加而減少。然而在圖5的左下角存在相關(guān)系數(shù)高值區(qū)域，說(shuō)明該區(qū)域內(nèi)觀測(cè)到的P波具有高相似性，這與布設(shè)的線性密集臺(tái)陣下方存在斷層帶有關(guān)，由于斷層帶中存在的破碎物質(zhì)使得穿過(guò)斷層的波形相似度變高。

3.2 近震延時(shí)分析

臺(tái)陣觀測(cè)到的直達(dá)P波具有較高的頻率，主頻率大于7 Hz （圖3）。通過(guò)分析近震事件直達(dá)P波到達(dá)每個(gè)臺(tái)站的時(shí)間，可以獲得更高分辨率的局部斷裂帶結(jié)構(gòu)信息。本研究選取了近震事件目錄中高信噪比（SNR>10）的垂直分量波形計(jì)算到時(shí)延時(shí)。

本文在2023年12月29日8：46：39發(fā)生在臺(tái)陣東側(cè)的近震事件波形中觀測(cè)到了明顯的體波放大效應(yīng)及走時(shí)延時(shí)，并對(duì)應(yīng)了具有低速特征的斷層帶（圖6中紅色線段）。

基于參考臺(tái)計(jì)算了每個(gè)近震事件P波、S波與各臺(tái)站間走時(shí)延時(shí)曲線，并計(jì)算了每個(gè)臺(tái)站100個(gè)近震事件的走時(shí)延遲平均值（圖7）。使用地表速度2.7 km/s對(duì)臺(tái)站高程進(jìn)行校正，選取從到時(shí)開(kāi)始窗長(zhǎng)為0.1 s 的P波及窗長(zhǎng)為0.2 s的S波波形，利用互相關(guān)法計(jì)算每個(gè)事件臺(tái)站間的P波與S波走時(shí)延時(shí)，誤差條表示100個(gè)近震事件走時(shí)延時(shí)曲線的兩個(gè)標(biāo)準(zhǔn)差。結(jié)果顯示，在臺(tái)陣西側(cè)存在明顯的走時(shí)延時(shí)特征，P波最大平均延時(shí)約6個(gè)采樣點(diǎn)（0.012 s），S波最大平均延時(shí)約15個(gè)采樣點(diǎn)（0.03 s），其對(duì)應(yīng)的斷層低速帶寬度約150 m（圖7中紅色線段），與P波在斷層低速帶的延時(shí)相比，S波具有更為明顯的高延時(shí)特征，反映出S波在計(jì)算到時(shí)延時(shí)的過(guò)程中比P波更敏感。

3.3 放大效應(yīng)分析

由于地震波的幾何擴(kuò)散與高衰減特性，在臺(tái)陣低速區(qū)觀測(cè)到遠(yuǎn)震波形的信噪比相對(duì)較低［35］。本研究用于體波放大效應(yīng)分析的近震事件均在距離臺(tái)陣30 km的范圍內(nèi)，且有較高信噪比，計(jì)算出不同方位近震事件窗長(zhǎng)為2 s的S波放大效應(yīng)。圖8顯示了PGV和RMS的乘積值歸一化曲線，該曲線可以表示臺(tái)陣下方介質(zhì)對(duì)近震事件S波的放大效應(yīng)。由于該方法對(duì)地殼低速斷層帶更敏感，可以代表不同臺(tái)站覆蓋區(qū)域出現(xiàn)斷層低速帶的概率。本文將所選取的近震事件根據(jù)不同方位進(jìn)行了劃分［圖8（a）］，且計(jì)算得到的結(jié)果一致性較高，表明該臺(tái)陣下方的低速介質(zhì)的方向差異性較小，這與斷層低速帶近似垂直分布有關(guān)。得到的斷層低速帶寬度與近震事件P波延時(shí)曲線、遠(yuǎn)震相關(guān)系數(shù)矩陣結(jié)果一致，寬度約150 m。

4 結(jié)論

甘肅積石山地震是近年來(lái)發(fā)生在青藏高原東北緣的又一次中強(qiáng)地震，研究該地震的發(fā)震斷層對(duì)于了解此次地震孕震機(jī)理及青藏高原東北緣地震風(fēng)險(xiǎn)評(píng)估都具有重要意義。發(fā)震斷層引起的巖石孔隙、壓力場(chǎng)和地震波速度的變化會(huì)導(dǎo)致線性密集臺(tái)陣記錄到的波形間產(chǎn)生差異，從而為研究斷層低速帶位置及寬度奠定理論基礎(chǔ)。許多地震學(xué)家通過(guò)布設(shè)密集臺(tái)陣，基于走時(shí)延時(shí)、放大效應(yīng)等方法開(kāi)展了小尺度斷層精細(xì)結(jié)構(gòu)方面的研究，并取得了豐碩成果。本文根據(jù)拉脊山斷裂走向、地表裂縫及余震分布在劉集鄉(xiāng)布設(shè)了一條線性密集臺(tái)陣，通過(guò)臺(tái)陣觀測(cè)到的近震體波走時(shí)延時(shí)分析、放大效應(yīng)，以及遠(yuǎn)震P波相關(guān)系數(shù)矩陣，在臺(tái)陣西側(cè)觀測(cè)到了一個(gè)寬度約150 m的斷層低速帶，其位于拉脊山北緣斷裂與拉脊山南緣斷裂之間的中部區(qū)域，并存在P波平均走時(shí)延時(shí)約6個(gè)采樣點(diǎn)（約為0.012 s），S波平均走時(shí)延時(shí)約15個(gè)采樣點(diǎn)（約為0.03 s）。通過(guò)放大效應(yīng)計(jì)算得到的斷層低速帶寬度與近震事件體波延時(shí)曲線和遠(yuǎn)震相關(guān)系數(shù)矩陣結(jié)果一致，寬度約150 m。

致謝：本研究中的斷層分布數(shù)據(jù)由蘭州大學(xué)袁道陽(yáng)教授提供，在此特別表示感謝。

參考文獻(xiàn)（References）

［1］ MOONEY W D，GINZBURG A.Seismic measurements of the internal properties of fault zones［J］.Pure and Applied Geophysics，1986，124（1）：141-157.

［2］ EVANS J P，CHESTER F M.Fluid-rock interaction in faults of the San Andreas system：inferences from San Gabriel fault rock geochemistry and microstructures［J］.Journal of Geophysical Research：Solid Earth，1995，100（B7）：13007-13020.

［3］ BEN-ZION Y，SAMMIS C G.Characterization of fault zones［J］.Pure and Applied Geophysics，2003，160（3）：677-715.

［4］ 趙莉，劉白云，范兵，等.2022年1月8日青海門源M6.9地震余震精確定位與斷層面擬合［J］.地震工程學(xué)報(bào)，2023，45（6）：1457-1466.

ZHAO Li，LIU Baiyun，F(xiàn)AN Bing，et al.Relocation and fault plane fitting of the Menyuan M6.9 earthquake on January 8，2022［J］.China Earthquake Engineering Journal，2023，45（6）：1457-1466.

［5］ 潘啟安，沈旭章.基于短周期密集地震臺(tái)陣觀測(cè)的空間自相關(guān)法及其在粵港澳大灣區(qū)的應(yīng)用［J］.地震學(xué)報(bào)，2023，45（2）：246-257.

PAN Qi'an，SHEN Xuzhang.Spatial autocorrelation method based on dense short-period seismic array and its application in the Guangdong—Hong Kong—Macao Greater Bay［J］.Acta Seismologica Sinica，2023，45（2）：246-257.

［6］ 周玲.新疆和田小孔徑臺(tái)陣監(jiān)測(cè)能力分析［J］.地震工程學(xué)報(bào)，2017，39（增刊1）：113-116.

ZHOU Ling.Analysis of monitoring ability of Hetian small-aperture seismic array［J］.China Earthquake Engineering Journal，2017，39（Suppl01）：113-116.

［7］ 王順，廖武林，姚運(yùn)生，等.測(cè)線方位對(duì)微動(dòng)線形臺(tái)陣探測(cè)效果影響研究［J］.地震工程學(xué)報(bào)，2023，45（2）：483-490.

WANG Shun，LIAO Wulin，YAO Yunsheng，et al.Influence of line orientation on the detection effect of a microtremor linear array［J］.China Earthquake Engineering Journal，2023，45（2）：483-490.

［8］ 徐順強(qiáng)，劉巧霞，李怡青，等.密集臺(tái)網(wǎng)微震定位技術(shù)在礦山開(kāi)采動(dòng)態(tài)監(jiān)測(cè)中的應(yīng)用研究［J］.地震工程學(xué)報(bào)，2015，37（1）：266-270.

XU Shunqiang，LIU Qiaoxia，LI Yiqing，et al.Application of the dense network micro-seismic location method to dynamic mining monitoring［J］.China Earthquake Engineering Journal，2015，37（1）：266-270.

［9］ 靳炎錫，秦滿忠，李少華.利用前驅(qū)波研究地幔間斷面的若干進(jìn)展［J］.地震工程學(xué)報(bào)，2023，45（1）：209-219，243.

JIN Yanxi，QIN Manzhong，LI Shaohua.Some progress in the study of mantle discontinuities using precursor waves［J］.China Earthquake Engineering Journal，2023，45（1）：209-219，243.

［10］ 李敏娟，沈旭章，張?jiān)?基于密集臺(tái)陣的青藏高原東北緣地殼精細(xì)結(jié)構(gòu)及九寨溝地震震源區(qū)結(jié)構(gòu)特征分析［J］.地球物理學(xué)報(bào)，2018，61（5）：2075-2087.

LI Minjuan，SHEN Xuzhang，ZHANG Yuansheng，et al.Fine crustal structures of northeast margin of the Tibetan Plateau and structural features of Jiuzhaigou earthquake focal area constrained by the data from a high-density seismic array［J］.Chinese Journal of Geophysics，2018，61（5）：2075-2087.

［11］ QIU H R，BEN-ZION Y，CATCHINGS R，et al.Seismic imaging of the MW 7.1 Ridgecrest earthquake rupture zone from data recorded by dense linear arrays［J］.Journal of Geophysical Research （Solid Earth），2021，126（7）：e2021JB022043.

［12］ QIN L，BEN-ZION Y，QIU H，et al.Internal structure of the San Jacinto fault zone in the trifurcation area southeast of Anza，California，from data of dense seismic arrays［J］.Geophysical Journal International，2018，213（1）：98-114.

［13］ MCGUIRE J，BEN-ZION Y.High-resolution imaging of the Bear Valley section of the San Andreas fault at seismogenic depths with fault-zone head waves and relocated seismicity［J］.Geophysical Journal International，2005，163（1）：152-164.

［14］ HILLERS G，CAMPILLO M，BEN-ZION Y，et al.Seismic fault zone trapped noise［J］.Journal of Geophysical Research：Solid Earth，2014，119（7）：5786-5799.

［15］ COCHRAN E S，LI Y G，SHEARER P M，et al.Seismic and geodetic evidence for extensive，long-lived fault damage zones［J］.Geology，2009，37（4）：315-318.

［16］ LEWIS M A，PENG Z，BEN-ZION Y，et al.Shallow seismic trapping structure in the San Jacinto fault zone near Anza，California［J］.Geophysical Journal International，2005，162（3）：867-881.

［17］ YANG W，PENG Z G，WANG B S，et al.Velocity contrast along the rupture zone of the 2010 MW6.9 Yushu，China，earthquake from fault zone head waves［J］.Earth and Planetary Science Letters，2015，416：91-97.

［18］ LEWIS M A，BEN-ZION Y，MCGUIRE J J.Imaging the deep structure of the San Andreas fault south of Hollister with joint analysis of fault zone head and direct P arrivals［J］.Geophysical Journal International，2007，169（3）：1028-1042.

［19］ 余大新，耿冠世.利用背景噪聲估計(jì)華北地區(qū)場(chǎng)地放大效應(yīng)［J］.地震工程學(xué)報(bào)，2017，39（4）：719-724.

YU Daxin，GENG Guanshi.Estimation of site amplification effect in North China from ambient seismic noise［J］.China Earthquake Engineering Journal，2017，39（4）：719-724.

［20］ KURZON I，VERNON F L，BEN-ZION Y，et al.Ground motion prediction equations in the San Jacinto fault zone：significant effects of rupture directivity and fault zone amplification［J］.Pure and Applied Geophysics，2014，171（11）：3045-3081.

［21］ LI Y G，AKI K，ADAMS D，et al.Seismic guided waves trapped in the fault zone of the Landers，California，earthquake of 1992［J］.Journal of Geophysical Research （Solid Earth），1994，99（B6）：11705-11722.

［22］ LI S L，LAI X L，YAO Z X，et al.Fault zone structures of northern and southern portions of the main central fault generated by the 2008 Wenchuan earthquake using fault zone trapped waves［J］.Earthquake Science，2009，22（4）：417-424.

［23］ 楊傳成，李智敏，熊仁偉，等.2023年甘肅積石山6.2級(jí)地震青海災(zāi)區(qū)地震地質(zhì)災(zāi)害初步調(diào)查分析［J］.地震，2024，44（1）：216-225.

YANG Chuancheng，LI Zhimin，XIONG Renwei，et al.Preliminary investigation and analysis of seismic geological hazards in Qinghai disaster area induced by the 2023 MS6.2 Jishishan earthquake［J］.Earthquake，2024，44（1）：216-225.

［24］ 楊攀新，熊仁偉，胡朝忠，等.2023年甘肅積石山6.2級(jí)地震發(fā)震構(gòu)造淺析［J］.地震，2024，44（1）：153-159.

YANG Panxin，XIONG Renwei，HU Chaozhong，et al.Seismogenic structure of the 2023 M6.2 Jishishan earthquake in Gansu Province［J］.Earthquake，2024，44 （1） ：153-159.

［25］ 李智敏，李延京，田勤儉，等.拉脊山斷裂古地震與喇家遺址災(zāi)變事件關(guān)系研究［J］.地震研究，2014，37（增刊1）：109-115.

LI Zhimin，LI Yanjing，TIAN Qinjian，et al.Study on the relationship between paleoseismic on Laji Mountain fault and catastrophic event on Lajiashan site［J］.Journal of Seismological Research，2014，37（Suppl01）：109-115.

［26］ 王遼，謝虹，蒲小武，等.2023年積石山地震誘發(fā)中川鄉(xiāng)滑坡-泥流滑動(dòng)過(guò)程與成因機(jī)理研究［J/OL］.地震工程學(xué)報(bào)，2024［2024-03-06］.https：//doi.org/10.2000/j.1000-0844.20240113003.

WANG Liao，XIE Hong，PU Xiaowu，et al.Sliding process and causative mechanism of the Zhongchuan landslide-mudflow disaster chain induced by the 2023 Jishishan earthquake［J/OL］.China Earthquake Engineering Journal，2024［2024-03-06］.https：//doi.org/10.2000/j.1000-0844.20240113003.

［27］ 左可楨，趙翠萍.2023年甘肅積石山6.2級(jí)地震序列精定位［J］.地震，2024，44（1）：204-208.

ZUO Kezhen，ZHAO Cuiping.Relocation of the 2023 MS6.2 Jishishan earthquake sequence in Gansu Province［J］.Earthquake，2024，44（1）：204-208.

［28］ BEYREUTHER M，BARSCH R，KRISCHER L，et al.ObsPy：a Python toolbox for seismology［J］.Seismological Research Letters，2010，81（3）：530-533.

［29］ SHARE P E，BEN-ZION Y，ROSS Z E，et al.Internal structure of the San Jacinto fault zone at Blackburn Saddle from seismic data of a linear array［J］.Geophysical Journal International，2017，210（2）：819-832.

［30］ YANG H F，LI Z F，PENG Z G，et al.Low-velocity zones along the San Jacinto fault，southern California，from body waves recorded in dense linear arrays［J］.Journal of Geophysical Research （Solid Earth），2014，119（12）：8976-8990.

［31］ 張勇，徐晨雨.2023/12/18積石山MW6.0地震初步破裂模型［R］.北京：北京大學(xué)地球物理學(xué)系張勇研究組，2023.

ZHANG Yong，XU Chenyu.Preliminary rupture model of the Jishishan MW6.0 earthquake on 2023/12/18［R］.Beijing：Zhang Yong Research Group，Department of Geophysics，Peking University，2023.

［32］ JIANG X H，HU S Q，YANG H F.Depth extent and vP/vS ratio of the Chenghai fault zone，Yunnan，China constrained from dense-array-based teleseismic receiver functions［J］.Journal of Geophysical Research （Solid Earth），2021，126（8）：e2021JB022190.

［33］ SONG J H，YANG H F.Seismic site response inferred from records at a dense linear array across the Chenghai fault zone，Binchuan，Yunnan［J］.Journal of Geophysical Research （Solid Earth），2022，127（1）：e2021JB022710.

［34］ ZHANG Z，DENG Y F，QIU H R，et al.High-resolution imaging of fault zone structure along the creeping section of the Haiyuan fault，NE Tibet，from data recorded by dense seismic arrays［J］.Journal of Geophysical Research （Solid Earth），2022，127（9）：e2022JB024468.

［35］ PENG Z G，BEN-ZION Y，MICHAEL A J，et al.Quantitative analysis of seismic fault zone waves in the rupture zone of the 1992 Landers，California，earthquake：evidence for a shallow trapping structure［J］.Geophysical Journal International，2003，155（3）：1021-1041.

（本文編輯：張向紅）

基金項(xiàng)目：國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目（42230305）；甘肅省科技重大專項(xiàng)（21ZD4FA011）；國(guó)家重點(diǎn)研發(fā)計(jì)劃課題（2018YFC1503206）

第一作者簡(jiǎn)介：吳柏辰（1999-），男，碩士研究生，主要從事密集臺(tái)陣與斷層結(jié)構(gòu)研究。E-mail：523402038@qq.com。

通信作者：秦滿忠（1981-），男，碩士，副研究員，主要從事地球內(nèi)部結(jié)構(gòu)、主動(dòng)源重復(fù)探測(cè)技術(shù)與應(yīng)用等方面的研究工作。E-mail：qinmzh@163.com。

吳柏辰，秦滿忠，郭曉，等.基于線性密集臺(tái)陣的2023年甘肅積石山MS6.2地震斷層低速帶研究［J］.地震工程學(xué)報(bào)，2024，46（3）：734-741.DOI：10.20000/j.1000-0844.20240306003

WU Baichen，QIN Manzhong，GUO Xiao，et al.Low-velocity fault zone of the 2023 Jishishan Gansu MS6.2 earthquake based on linear dense array data［J］.China Earthquake Engineering Journal，2024，46（3）：734-741.DOI：10.20000/j.1000-0844.20240306003