丁龍兵, 胡進(jìn)軍, 李培旭, 劉名吉
(1.中國地震局工程力學(xué)研究所 地震工程與工程振動(dòng)重點(diǎn)實(shí)驗(yàn)室,哈爾濱 150080;2.中國地震局工程力學(xué)研究所 地震災(zāi)害防治應(yīng)急管理部重點(diǎn)實(shí)驗(yàn)室,哈爾濱 150080)
由近斷層效應(yīng)引發(fā)的速度大脈沖現(xiàn)象,在過去的50年內(nèi)已經(jīng)被證實(shí)和分析。在1994年美國Northridge 6.7級(jí)地震、1995年日本Kobe 7.2級(jí)地震、1999年我國臺(tái)灣集集7.6級(jí)等地震中,近斷層獨(dú)特的長(zhǎng)周期脈沖特性對(duì)工程結(jié)構(gòu)造成嚴(yán)重的破壞更是引起廣泛的關(guān)注[1-3]。這種近斷層長(zhǎng)周期脈沖型地震動(dòng)的形成與震源機(jī)制、斷層破裂的方向性和場(chǎng)地特征有較大的關(guān)聯(lián)[4-5]。近斷層脈沖型地震動(dòng)具有更為復(fù)雜的特性以及對(duì)工程結(jié)構(gòu)的嚴(yán)重破壞性[6-7],使其成為地震工程領(lǐng)域的研究熱點(diǎn)之一。
基于數(shù)據(jù)統(tǒng)計(jì)的方法是目前研究脈沖型地震動(dòng)特征的主要分析方法,謝俊舉等[8]利用165條典型速度脈沖記錄分析原始脈沖記錄和提取脈沖后的殘余記錄進(jìn)行反應(yīng)譜特征研究,結(jié)果表明速度大脈沖在其特征周期Tp附近對(duì)反應(yīng)譜有著較強(qiáng)的放大作用。趙曉芬等[9]通過分析317條速度脈沖記錄的速度脈沖放大作用曲線,同樣證實(shí)了速度大脈沖在其特征周期Tp附近對(duì)反應(yīng)譜有很強(qiáng)的放大作用,并提出了速度脈沖放大作用模型。Somerville[10]利用基巖場(chǎng)地和土層場(chǎng)地記錄到的速度脈沖進(jìn)行對(duì)比分析,研究表明土層場(chǎng)地的脈沖周期大于基巖場(chǎng)地(Mw<7.5)。Yang等[11]選取1999臺(tái)灣集集地震和1994年美國Northridge近斷層地震記錄進(jìn)行頻譜特征分析,結(jié)果表明長(zhǎng)周期脈沖對(duì)近斷層地震動(dòng)頻譜特征影響較大,這種地震動(dòng)長(zhǎng)周期分量主要出現(xiàn)在下盤場(chǎng)地。脈沖型地震動(dòng)對(duì)中長(zhǎng)周期工程結(jié)構(gòu)有著重要影響,這體現(xiàn)在脈沖效應(yīng)對(duì)中長(zhǎng)周期段加速度反應(yīng)譜有著明顯的放大作用。江輝等[12]基于長(zhǎng)周期的近斷層脈沖型地震動(dòng)記錄,通過分析豎向分量與水平向分量的加速度反應(yīng)譜比值與脈沖周期和結(jié)構(gòu)自振周期的變化規(guī)律,表明我國現(xiàn)行規(guī)范取值明顯小于反應(yīng)譜比值。
2023年2月6日,土耳其連續(xù)發(fā)生了Mw7.8級(jí),Mw7.5級(jí)和Mw6.7級(jí)地震,2月20日又發(fā)生了Mw6.3級(jí)地震。這是土耳其歷史上最嚴(yán)重的一次地震,截至2月28日,土耳其國內(nèi)累計(jì)死亡人數(shù)超過44 300人,另外16萬座建筑被毀壞,造成150萬人無家可歸以及不可估量的經(jīng)濟(jì)損失。特別是靠近斷層區(qū)域的卡赫拉曼馬拉什、阿德亞曼等地區(qū)建筑物損毀得更為嚴(yán)重。土耳其災(zāi)害和應(yīng)急管理局(Disaster, and Emergency Management Presidency, AFAD) 數(shù)據(jù)庫獲得了大量的近斷層強(qiáng)震動(dòng)加速度記錄。為了揭示本次地震序列中近斷層地震動(dòng)的特征,本文使用能量法對(duì)近斷層脈沖型地震動(dòng)記錄進(jìn)行識(shí)別,分析四次強(qiáng)震中速度脈沖周期和幅值特征,研究震級(jí)、斷層距對(duì)幅值等參數(shù)的影響。
土耳其Mw=7.8、Mw=7.5、Mw=6.7和Mw=6.3地震序列中,土耳其AFAD分別公布了289組,216組,160組和116組三分量加速度記錄,選取上述781組加速度記錄進(jìn)行研究。根據(jù)美國地質(zhì)調(diào)查局(United States Geological Survey, USGS)和土耳其AFAD發(fā)布的震源機(jī)制信息,四次地震的震源信息如表1所示。地震斷層根據(jù)USGS公布的斷層信息繪制,近斷層臺(tái)站分布及四次地震的震源機(jī)制信息如圖1所示。
圖1 土耳其地震的強(qiáng)震動(dòng)臺(tái)站和脈沖型地震動(dòng)分布
表1 土耳其地震序列的震源機(jī)制
識(shí)別脈沖型地震動(dòng)[13-17]的方法有多種,其中能量法是一種簡(jiǎn)單且有效識(shí)別脈沖型地震動(dòng)的方法,對(duì)于本次土耳其地震序列獲得的強(qiáng)震動(dòng)記錄,能量法能夠準(zhǔn)確識(shí)別出更多的速度脈沖型地震動(dòng)尤其是含有多個(gè)脈沖的單條速度脈沖型地震動(dòng),結(jié)果表明其適用于本次地震序列。因此,本文基于能量法對(duì)地震動(dòng)記錄進(jìn)行預(yù)處理和識(shí)別。
考慮到未經(jīng)處理的加速度記錄可能存在基線偏差和噪聲造成積分后的速度和位移時(shí)程不合理偏移的現(xiàn)象,為了消除這一現(xiàn)象,需要對(duì)地震動(dòng)記錄進(jìn)行基線校準(zhǔn)和濾波處理[18-20],具體步驟如下:
步驟1 從連續(xù)波形中將加速度記錄進(jìn)行截?cái)?將獲取的整個(gè)事件記錄減去事件前5 s的均值得到加速度記錄。
步驟2 將加速度記錄前后補(bǔ)零30 s后,再進(jìn)行波形尖滅。
步驟3 進(jìn)行濾波,選取4階Butterworth非因果濾波器進(jìn)行帶通濾波,濾波選取的頻率段選擇0.1~30 Hz。
步驟4 最后將得到的加速度、速度和位移時(shí)程記錄進(jìn)行人工篩選,剔除無法使用的記錄。
圖2為編號(hào)3124臺(tái)站在Mw7.8地震中獲得的未處理和處理過的東西(EW)方向加速度、速度和位移記錄,經(jīng)過處理后的地震動(dòng)記錄明顯消除了基線漂移的現(xiàn)象。
(a) 處理前
為了提高脈沖識(shí)別的效率以及確定地震動(dòng)記錄中最顯著的速度脈沖方向便于研究其脈沖特性,需要對(duì)于地震動(dòng)記錄的旋轉(zhuǎn)方向進(jìn)行選取。Shahi等[21]曾觀察到,如果在任何方向上發(fā)現(xiàn)脈沖,則地面運(yùn)動(dòng)被標(biāo)記為脈沖型地震動(dòng),本文為了考慮最大速度脈沖方向的不確定性,采用將實(shí)際觀測(cè)水平地震動(dòng)記錄在水平面內(nèi)進(jìn)行旋轉(zhuǎn),進(jìn)而識(shí)別地震動(dòng)在水平0°~180°方位上的最大速度脈沖方法。給定任意兩個(gè)正交水平地震動(dòng)分量按式(1)的關(guān)系可以得到水平向任意方位地震動(dòng)分量
v(t,θ)=v1(t)×cos(θ)+v2(t)×sin(θ)
(1)
式中:v1(t)和v2(t)依次代表臺(tái)站記錄東西(EW)方向與南北(NS)方向的地震動(dòng)速度時(shí)程,以逆時(shí)針旋轉(zhuǎn)為正,則v1(t)旋轉(zhuǎn)90°后與v2(t)的方位重合;v(t,θ)為與v1(t)成任意角度θ方位對(duì)應(yīng)的地震動(dòng)速度時(shí)程。對(duì)于所有方向的地震動(dòng)分量依次進(jìn)行脈沖識(shí)別增加了計(jì)算的時(shí)間和降低了脈沖識(shí)別的效率,Li等[22]發(fā)現(xiàn)最大峰值速度(peak ground velocity,PGV)所在方向在統(tǒng)計(jì)意義上代表著最強(qiáng)脈沖所在的方向,因此本節(jié)將781組水平向兩分量加速度記錄積分得到的速度時(shí)程記錄旋轉(zhuǎn)至最大PGV所在方向進(jìn)行預(yù)處理,將旋轉(zhuǎn)后的記錄作為水平向脈沖識(shí)別的對(duì)象,對(duì)于豎向(UD)地震動(dòng)分量不做旋轉(zhuǎn)處理進(jìn)行識(shí)別。
圖3給出了臺(tái)站地震動(dòng)速度時(shí)程記錄旋轉(zhuǎn)前后的對(duì)比,發(fā)現(xiàn)將速度時(shí)程記錄旋轉(zhuǎn)到最大PGV方向后,可以明顯看出高幅值的速度大脈沖,盡管在EW、NS、UD方向均可觀察到不同程度上的速度脈沖,但在脈沖的形狀和幅值大小上均與最大PGV方向上的速度脈沖存在差異,不具備統(tǒng)計(jì)特征。因此選取最大PGV方向的脈沖進(jìn)行脈沖識(shí)別更加高效。
(a) 2717臺(tái)站旋轉(zhuǎn)前后
能量法通過對(duì)地震動(dòng)記錄速度的平方在時(shí)間尺度進(jìn)行積分,定義了顯著半脈沖的概念和統(tǒng)計(jì)回歸得到速度脈沖型地震動(dòng)的判定準(zhǔn)則,研究表明能量法可以更加有效地識(shí)別速度脈沖型地震動(dòng)。以本次土耳其Mw7.8地震中2718臺(tái)站NS方向的地震動(dòng)記錄為例,圖4給出了半脈沖振動(dòng)區(qū)間和速度時(shí)程的關(guān)系。半脈沖指的是連續(xù)兩個(gè)零點(diǎn)之間的速度時(shí)程,每個(gè)速度半波的能量以兩個(gè)連續(xù)零點(diǎn)之間速度平方的時(shí)間積分來表示,則地震動(dòng)的總能量可以表示為
(2)
式中:N為半脈沖的總數(shù);ΔEk為第k個(gè)速度半脈沖的能量,計(jì)算公式如下
(3)
式中:v(t)為速度時(shí)間序列;t1k和t2k為第k個(gè)速度半脈沖的開始時(shí)間和結(jié)束時(shí)間。為了評(píng)價(jià)半脈沖的重要性和顯著性,第k個(gè)速度半脈沖的相對(duì)能量被定義為
Ek=ΔEk/E
(4)
將Ek≥0.1的半脈沖定義為顯著半脈沖。
根據(jù)脈沖型地震動(dòng)的特征,將脈沖能表示為所有顯著半脈沖所含的能量之和,即
(5)
將所有地震動(dòng)按照所含顯著半脈沖的個(gè)數(shù)分為不同的類別,包括類別0、類別1、類別2、類別3、類別4和類別5,分別代表地震動(dòng)速度時(shí)程中無顯著半脈沖,含有1個(gè),2個(gè),3個(gè),4個(gè)和至少5個(gè)顯著半脈沖。當(dāng)?shù)卣饎?dòng)為類別0時(shí)則為非脈沖,表2給出了其余類別的臨界閾值,當(dāng)?shù)卣饎?dòng)的EP大于相應(yīng)的臨界閾值則判別為脈沖型。
表2 不同類別的臨界閾值
采用上述方法及準(zhǔn)則對(duì)四次地震781組三分量加速度記錄的水平向和豎向分別進(jìn)行識(shí)別。考慮到脈沖周期(Tp)定義的不同會(huì)影響Tp的計(jì)算[23],本文對(duì)識(shí)別出的脈沖型地震動(dòng)記錄采用取速度反應(yīng)譜峰值對(duì)應(yīng)的周期作為Tp的方法[24]和峰點(diǎn)法[25]相結(jié)合確定Tp。表3和表4分別為豎向和水平向的識(shí)別結(jié)果及脈沖參數(shù)。
表3 豎向脈沖型地震動(dòng)的識(shí)別結(jié)果
表4 水平向脈沖型地震動(dòng)的識(shí)別結(jié)果
從表3和表4可以看出,水平向上識(shí)別出了45條速度脈沖記錄,其中有30條記錄來自于Mw7.8地震,5條記錄來自于Mw7.5地震,9條記錄來自于Mw6.7地震,剩余的1條記錄來自于Mw6.3地震。在豎向上識(shí)別出了17條速度脈沖記錄,其中14條來自于Mw7.8地震,1條記錄來自于Mw7.5地震,剩余的2條記錄來自于Mw6.3地震。脈沖記錄主要來自于第一次Mw7.8地震。其中豎向的17條速度脈沖記錄有12條記錄的所在臺(tái)站的水平向記錄上同樣也被歸類為速度脈沖,但3316臺(tái)站和4615臺(tái)站所獲得的Mw7.8地震動(dòng)記錄在水平向被歸類為非速度脈沖,而在豎向被歸類為速度脈沖。水平方向上Tp主要分布在1~6 s之間,其中8002臺(tái)站獲得的Mw7.8地震的記錄的脈沖周期最大,為7.57 s。水平方向上PGV主要分布在30~120 cm/s之間,其中3138臺(tái)站獲得的Mw7.8地震的記錄的PGV最大,為231.30 cm/s。豎向Tp主要分布在1~4 s之間,其中3136臺(tái)站獲得的Mw7.8地震記錄的Tp最大,為5.86 s。豎向PGV主要分布在30~60 cm/s之間,其中3138臺(tái)站獲得的Mw7.8地震記錄的PGV最大,為83.59 cm/s。圖5為水平向和豎向部分臺(tái)站的速度脈沖記錄。
(a) 水平向速度脈沖記錄
脈沖型地震動(dòng)通常可以用四個(gè)參數(shù)來準(zhǔn)確描述[26]:脈沖周期(Tp)、脈沖幅值(Vp)、脈沖個(gè)數(shù)以及相位。由于目前脈沖個(gè)數(shù)和相位的確定的研究工作尚不充分,脈沖型地震動(dòng)的主要的特性一般可以通過脈沖周期和脈沖幅值來表征[27]。由于17條豎向脈沖型地震動(dòng)數(shù)量較少且關(guān)于豎向脈沖型地震動(dòng)特性研究較少,缺乏相應(yīng)模型比較,因此本文研究45組水平向脈沖型地震動(dòng)和脈沖周期與矩震級(jí)的關(guān)系和脈沖幅值隨斷層距和矩震級(jí)的變化。
利用土耳其四次地震識(shí)別出的45條水平向速度脈沖型地震動(dòng)記錄對(duì)脈沖周期和矩震級(jí)的關(guān)系進(jìn)行統(tǒng)計(jì)分析,圖6給出了由45條水平向脈沖型記錄數(shù)據(jù)回歸得到的脈沖周期隨矩震級(jí)變化的關(guān)系曲線。本文統(tǒng)計(jì)回歸得到的脈沖周期隨矩震級(jí)變化的預(yù)測(cè)關(guān)系式為
圖6 本文回歸得到的速度脈沖周期隨震級(jí)的變化關(guān)系
lgTp=-2.139+0.338Mw
(6)
圖7將回歸得到的脈沖周期和震級(jí)的預(yù)測(cè)模型與Alavi和Krawinkler[24]、Somerville[10]、Mavroeidis和Papageorgiou[26]、Bray和Rodriguez-Marek[28]、趙國臣等[29]以及 Shahi和Baker[15]給出的經(jīng)驗(yàn)關(guān)系進(jìn)行對(duì)比。由圖7可知,本文得到的預(yù)測(cè)模型與現(xiàn)有的預(yù)測(cè)模型存在明顯差異,在矩震級(jí)大于6.4時(shí),相同矩震級(jí)下,本文模型給出的脈沖周期預(yù)測(cè)值均低于現(xiàn)有模型的脈沖周期預(yù)測(cè)值。在矩震級(jí)小于6.4時(shí),本文模型給出的脈沖周期預(yù)測(cè)值僅僅只高于Somerville[10]模型在巖石場(chǎng)地下的預(yù)測(cè)值。本次土耳其地震的脈沖型地震動(dòng)的脈沖周期較以往多數(shù)脈沖型地震動(dòng)來說顯著偏低。
圖7 本文回歸得到的脈沖周期和震級(jí)的經(jīng)驗(yàn)關(guān)系與現(xiàn)有模型的比較
Cork等的研究表明對(duì)于震級(jí)較小(Mw<7.5)的地震,土壤場(chǎng)地的脈沖周期比巖石場(chǎng)地更長(zhǎng)。隨著地震震級(jí)的增加,巖石和土壤場(chǎng)地的脈沖周期值趨于一致,這是因?yàn)檎鸺?jí)較大的地震的波長(zhǎng)大于30 m。因此在地表以下30 m的上部,將場(chǎng)地分為巖石或土壤,對(duì)震級(jí)較大地震的脈沖周期幾乎沒有影響。本文選用Shahi等模型選取的243條脈沖型地震動(dòng)記錄(https://www.jackwbaker.com/pulse_classification_v2/Pulse-like-records.html)以及Bray等模型選取的54條脈沖型地震動(dòng)記錄與本文選用的45條土耳其地區(qū)脈沖型地震動(dòng)記錄對(duì)比分析場(chǎng)地條件對(duì)脈沖周期的影響。依據(jù)美國抗震設(shè)計(jì)規(guī)范[30]場(chǎng)地劃分規(guī)則,根據(jù)地表以下30 m深度范圍內(nèi)的等效剪切波速(VS30)值進(jìn)行分類,結(jié)果如圖8所示。
圖8 本文所采用的脈沖型地震動(dòng)記錄與前人模型所用的脈沖型地震動(dòng)記錄的場(chǎng)地分類情況
從圖8可以看出,本次土耳其脈沖型地震動(dòng)記錄較選取的前人模型的脈沖型地震動(dòng)記錄來看,場(chǎng)地條件整體偏硬,這可能是小震級(jí)下(Mw<7.5)導(dǎo)致脈沖周期偏小的原因。
另外脈沖周期與地震動(dòng)記錄的頻率成分息息相關(guān),前人模型中矩震級(jí)7.5級(jí)及以上的脈沖型地震動(dòng)記錄絕大多數(shù)來自于1999年臺(tái)灣集集地震,其中Baker提供的脈沖型地震動(dòng)記錄較其余模型的數(shù)據(jù)來說,最為全面。我們以Baker提供的脈沖型地震動(dòng)記錄為例,對(duì)比本文采用的35條土耳其脈沖型地震動(dòng)記錄(Mw≥7.5)和Baker提供的45條脈沖型地震動(dòng)記錄數(shù)據(jù)庫(Mw≥7.5)的頻率成分。
由于脈沖型地震動(dòng)記錄低頻成分較普通地震動(dòng)低頻成分更豐富,希爾波特譜可以給出地震信號(hào)中其他方法可能遺漏的低頻能量[31]。因此這里我們使用希爾伯特能量譜去計(jì)算脈沖型地震動(dòng)記錄低頻(小于0.3 Hz)頻段的能量與總能量的比率,低頻頻段截止頻率選0.3 Hz是因?yàn)榍叭祟A(yù)測(cè)模型在Mw≥7.5時(shí)的脈沖周期均大于3.33 s(對(duì)應(yīng)頻率為0.3 Hz)。
圖9為本文選取土耳其地區(qū)脈沖型地震動(dòng)記錄與Baker脈沖型地震動(dòng)記錄的低頻成分能量占比。我們發(fā)現(xiàn)本文選取的土耳其地區(qū)脈沖型地震動(dòng)記錄的低頻成分能量占比遠(yuǎn)小于前人所用的脈沖型地震動(dòng)記錄,這是導(dǎo)致大震級(jí)(Mw≥7.5)下本文模型顯著小于前人模型的主要原因。
圖9 土耳其地區(qū)脈沖型地震動(dòng)記錄(Mw≥7.5)與Baker提供的脈沖型地震動(dòng)記錄(Mw≥7.5)的低頻成分(小于0.3 Hz)能量占比
PGV在時(shí)域內(nèi)可以合理并準(zhǔn)確地代表脈沖幅值,因此本文使用PGV來表示脈沖幅值。選取土耳其地震的45條脈沖記錄的脈沖幅值數(shù)據(jù),圖10繪出了脈沖幅值隨震級(jí)和斷層距變化的分布,總體來看,矩震級(jí)相同時(shí),脈沖幅值隨斷層距的增大而顯著減小。斷層距相近時(shí),脈沖幅值隨矩震級(jí)增大而增大。統(tǒng)計(jì)回歸得到的脈沖周期隨震級(jí)和斷層距的預(yù)測(cè)關(guān)系式為
圖10 本文回歸得到的速度脈沖幅值隨震級(jí)和斷層距的變化關(guān)系
lgVp=2.056-0.032Mw-0.550lgRrup
(7)
圖11為文獻(xiàn)[23,28,32-33]提出的脈沖幅值的統(tǒng)計(jì)回歸模型,與本文計(jì)算出的土耳其脈沖型地震動(dòng)的脈沖幅值和模型的比較。從圖11可以發(fā)現(xiàn),本文模型明顯高于Bray等巖石場(chǎng)地上的預(yù)測(cè)值。矩震級(jí)相同的情況下,當(dāng)斷層距小于25 km左右時(shí),本文模型預(yù)測(cè)值與Bray等巖石場(chǎng)地模型預(yù)測(cè)值相近且大于其他模型預(yù)測(cè)值。隨著斷層距不斷變大(Rrup>25 km),本文模型小于Tang等土層場(chǎng)地與巖石場(chǎng)地模型預(yù)測(cè)值,隨著矩震級(jí)不斷減小,本文模型與Chen等走滑斷層的預(yù)測(cè)模型逐漸接近。
圖11 本文回歸得到的脈沖幅值隨震級(jí)和斷層距的經(jīng)驗(yàn)關(guān)系與現(xiàn)有模型的比較
基于能量法對(duì)土耳其地震序列的四次強(qiáng)震記錄在水平方向和豎向上進(jìn)行脈沖型地震動(dòng)的分類識(shí)別,研究了水平方向上震級(jí)對(duì)脈沖周期的影響,分析了震級(jí)、斷層距對(duì)脈沖周期和幅值的影響,揭示了土耳其四次地震的速度脈沖特點(diǎn),研究表明:
(1) 對(duì)土耳其地震序列獲得的781組地震動(dòng)記錄進(jìn)行最大PGV方向上的旋轉(zhuǎn)后,利用能量法進(jìn)行速度脈沖識(shí)別,發(fā)現(xiàn)土耳其地震序列中存在水平向和豎向的近斷層脈沖型地震動(dòng)。
(2) 土耳其地震序列脈沖型地震動(dòng)的速度脈沖具有周期短、幅值大的顯著特點(diǎn)。水平向和豎向速度脈沖周期主要分布在1~6 s和1~4 s之間,豎向速度脈沖周期總體上小于水平向,水平向和豎向速度脈沖周期最大值7.57 s和5.86 s分別出現(xiàn)在8002臺(tái)站和3136臺(tái)站在Mw7.8地震中獲得的地震動(dòng)記錄中。水平向和豎向速度脈沖幅值主要分布在30~120 cm/s和30~60 cm/s,豎向速度脈沖幅值明顯小于水平向。在水平向和豎向3138臺(tái)站獲得的Mw7.8地震記錄的脈沖幅值均為最大,分別為231.30 cm/s和83.59 cm/s。
(3) 水平方向上速度脈沖的周期隨矩震級(jí)呈對(duì)數(shù)線性增大且相同矩震級(jí)情況下,土耳其地震序列脈沖型地震動(dòng)的脈沖周期明顯小于現(xiàn)有模型的預(yù)測(cè)值。脈沖幅值隨斷層距的增加存在明顯衰減的趨勢(shì),當(dāng)斷層距小于25 km時(shí),脈沖幅值普遍高于現(xiàn)有模型預(yù)測(cè)值。
致謝
感謝土耳其災(zāi)害和應(yīng)急管理局 (Disaster, and Emergency Management Presidency, www.afad.gov.tr) 提供的地震動(dòng)數(shù)據(jù)和美國地質(zhì)調(diào)查局(United States Geological Survey, https://earthquake.usgs.gov/)提供的震源機(jī)制信息。