潘章容,李同林,崔建文,田秀豐,石文兵,張衛(wèi)東,袁 潔

(1.甘肅省地震局,甘肅 蘭州 730000;2.四川省地震局,四川 成都 610041;3.云南省地震局,云南 昆明 650224)

0 引言

2021年5月21日21時(shí)48分34秒,云南省漾濞縣(25.67°N,99.87°E)發(fā)生M6.4地震,震源深度8.0 km。據(jù)中國地震臺網(wǎng)目錄,截至2021年6月15日,漾濞縣共發(fā)生M3.0以上余震40次,其中M3.0～3.9地震29次,M4.0～4.9地震9次,M5.0以上地震2次,最大余震為2021年5月21日22時(shí)31分10秒(震后43分鐘)發(fā)生的M5.2地震(25.59°N,99.97°E;震源深度8.0 km)。在M6.4主震發(fā)生前26分鐘(2021年5月21日21時(shí)21分25秒),漾濞縣發(fā)生了一次M5.6地震(25.63°N,99.92°E;震源深度10.0 km)。

中國數(shù)字強(qiáng)震動(dòng)觀測網(wǎng)絡(luò)在2021年5月21日云南省漾濞縣前后70分鐘內(nèi)發(fā)生的3次5級以上地震(M5.6、M6.4及M5.2)中,捕獲了豐富的強(qiáng)震動(dòng)觀測記錄?；诖?本文擬對3次地震中獲取的加速度記錄進(jìn)行初步處理與分析,繪制PGA等值線分布圖,并對比研究3次地震中加速度幅值衰減特性及加速度反應(yīng)譜特性。

1 強(qiáng)震動(dòng)記錄收集及處理

1.1 觸發(fā)臺站分布

2021年5月21日云南漾濞M5.6、M6.4及M5.2地震中,中國數(shù)字強(qiáng)震動(dòng)觀測網(wǎng)絡(luò)共獲取了41組三分向加速度記錄,其中M5.6地震9組,M6.4地震27組(云南省強(qiáng)震動(dòng)臺網(wǎng)25組、四川省強(qiáng)震動(dòng)臺網(wǎng)2組),M5.2地震5組。3次地震近場觸發(fā)臺站(震中距<150 km)分布見圖1。由于53YPX(永平臺)UD向數(shù)據(jù)異常,未采用該臺站的UD向加速度記錄,因此在3次地震中共獲取波形完整、震相清晰的加速度記錄120條,其中M5.6地震26條、M6.4地震80條、M5.2地震14條。

圖1 近場觸發(fā)臺站分布圖(震中距<150 km)Fig.1 Distribution of near-field triggered stations(epicentral distance <150 km)

1.2 地震記錄分析

對上述120條未校正的加速度記錄進(jìn)行格式轉(zhuǎn)換、基線校正、濾波等常規(guī)數(shù)據(jù)處理[3-7]。首先根據(jù)《中國地震烈度表》(GB/T 17742—2020)[8]中給定方法計(jì)算PGA(合成地震動(dòng)加速度記錄的最大值)和PGV(合成地震動(dòng)速度記錄的最大值),再按式(1)、式(2)計(jì)算應(yīng)用PGA得到的地震烈度計(jì)算值IA及應(yīng)用PGV得到的地震烈度計(jì)算值IV:

IA=3.17lg(PGA)+6.59

(1)

IV=3.00lg(PGV)+9.77

(2)

按式(3)計(jì)算儀器測定的地震烈度II，結(jié)果取小數(shù)點(diǎn)后一位有效數(shù)字；II<1.0時(shí)取1.0，II>12.0時(shí)取12.0:

(3)

最后按照《中國地震烈度表》(GB/T 17742—2020)中地震烈度與儀器測定的地震烈度II的對應(yīng)關(guān)系,計(jì)算地震烈度值。

將云南漾濞M5.6、M6.4及M5.2地震典型加速度記錄的基本信息及相關(guān)參數(shù)分別列于表1、表2和表3。對有效三分量加速度記錄進(jìn)行初步處理及特征分析,為大理州及周邊地區(qū)的地震區(qū)劃、抗震設(shè)防、場地效應(yīng)、烈度衰減關(guān)系等研究提供基礎(chǔ)資料[6-7]。

表1、表2及表3顯示,M5.6地震獲取記錄的PGA介于0.8～338.9 cm/s2間,最大PGA由距離震中最近的53YBX(漾濞臺)EW向獲取,震中距為5.37 km;M6.4地震獲取記錄的PGA介于0.8～719.6 cm/s2間,最大PGA由震中距為7.77 km的53YBX臺NS向獲取;M5.2地震獲取記錄的PGA介于2.0～190.6 cm/s2間,最大PGA由震中距為9.12 km的53YBX臺NS向獲取。

表1 漾濞M5.6地震強(qiáng)震動(dòng)記錄及其相關(guān)參數(shù)Table 1 Strong motion records and related parameters of the M5.6 earthquake

表2 漾濞M6.4地震強(qiáng)震動(dòng)記錄及其相關(guān)參數(shù)Table 2 Strong motion records and related parameters of the M6.4 earthquake

表3 漾濞M5.2地震獲取強(qiáng)震動(dòng)記錄及其相關(guān)參數(shù)Table 3 Strong motion records and related parameters of the M5.2 earthquake

整體上,3次地震的PGA、PGV值均隨震中距的增加而衰減,個(gè)別異常點(diǎn)可能與臺站本身的場地條件及地震傳播特性有關(guān),還有待進(jìn)一步研究。

2 峰值加速度分布圖

峰值加速度分布圖以各觀測點(diǎn)的加速度值為依據(jù)(通常選取三分向中的某一分向或最大分向來進(jìn)行繪制),可真實(shí)反映地震發(fā)生時(shí)觀測點(diǎn)的地面運(yùn)動(dòng)強(qiáng)度,能夠?yàn)檠芯繌?qiáng)地面運(yùn)動(dòng)的特性和工程結(jié)構(gòu)抗震設(shè)計(jì)方法與技術(shù)提供基礎(chǔ)資料,同時(shí)可為快速評估和制定大震應(yīng)急方案提供重要依據(jù)[7]。

以表1、表2及表3中計(jì)算得到的三分向峰值加速度值為基礎(chǔ),選取三分向的最大峰值加速度值為數(shù)據(jù)源,采用克里格插值法分別繪制了云南漾濞M5.6、M6.4、M5.2地震的峰值加速度等值線分布于圖2。由圖2可知,3次地震的峰值加速度等值線均呈NW-SE方向展布,最大值位于震中附近。

圖2 峰值加速度等值線分布圖Fig.2 PGA contour distribution map

云南省地震局發(fā)布的《云南漾濞6.4級地震烈度圖》(圖3)顯示,此次地震的最高烈度為Ⅷ度,主要分布在漾濞縣蒼山西鎮(zhèn)、漾江鎮(zhèn)、太平鄉(xiāng)3個(gè)鄉(xiāng)鎮(zhèn)。根據(jù)《中國地震烈度表》(GB/T 17742—2020)中地震烈度與儀器測定的地震烈度的對應(yīng)關(guān)系,儀器測定的最大值點(diǎn)(63YBX臺)所在位置對應(yīng)地震烈度為Ⅷ度,與《云南漾濞6.4級地震烈度圖》基本一致。這說明儀器測定的地震烈度值與宏觀調(diào)查的地震烈度值基本符合,測量值可為確定地震烈度提供參考。

圖3 云南漾濞6.4級地震烈度圖(云南省地震局2021年5月25日發(fā)布)Fig.3 The seismic intensity map of Yangbi M6.4 earthquake (Published by Yunnan Earthquake Agency on May 25,2021)

3 幅值衰減特征

地震動(dòng)在傳播過程中受到震源和傳播路徑介質(zhì)的影響會(huì)呈現(xiàn)出不同的衰減特性,目前學(xué)者已提出多種衰減關(guān)系[1-2,9],本文選用國內(nèi)兩種常用的衰減模型與本次地震動(dòng)參數(shù)進(jìn)行對比。

結(jié)合地震發(fā)生區(qū)域的地質(zhì)特征和強(qiáng)震動(dòng)臺站位置,選用第五代《中國地震動(dòng)參數(shù)區(qū)劃圖(GB 18306—2015)》青藏區(qū)(包括川滇地區(qū))地震動(dòng)衰減關(guān)系[1](簡稱YU)和汪素云中國西部地震動(dòng)衰減關(guān)系[2](簡稱WANG)分析M5.6、M6.4及M5.2地震強(qiáng)震動(dòng)記錄的幅值衰減特性。為進(jìn)行對比分析,將YU衰減關(guān)系的適用范圍延長至400 km。3次地震觀測PGA與衰減關(guān)系預(yù)測PGA間的對比關(guān)系見圖4,其中實(shí)線和虛線分別表示衰減關(guān)系長軸及短軸曲線。

圖4 青藏地區(qū)與中國西部地震衰減關(guān)系Fig.4 Earthquake attenuation relationship in the Qinghai-Tibet area and western China

由圖4可見,地震動(dòng)PGA觀測值衰減與地震震級、震中距及地震傳播介質(zhì)等因素有關(guān),具體表現(xiàn)為:

(1)與上述兩種衰減關(guān)系相比較,PGA實(shí)際觀測值與YU預(yù)測值更吻合,與WANG預(yù)測值離散性較大。

(2)在近場震中距小于10 km的53YBX臺,3次地震的PGA觀測值與YU預(yù)測值基本一致,位于WANG衰減曲線上方,即該臺站PGA觀測值大于WANG預(yù)測值;在震中距介于10～30 km的53DLY(月溪井臺),3次地震的PGA觀測值與YU預(yù)測值、WANG預(yù)測值均吻合較好;在震中距介于60～65 km的53BTH(太和臺),3次地震的PGA觀測值均位于兩種預(yù)測曲線下方,即觀測值小于預(yù)測值,且從表1、表2和表3中可以看出該臺三分向加速度峰值均小于類似震中距臺站的加速度峰值,其實(shí)際情況有待進(jìn)一步研究。

(3)震中距大于200 km的遠(yuǎn)場PGA觀測值與YU預(yù)測值、WANG預(yù)測值均吻合較好,PGA觀測值較均勻地分布于兩種衰減曲線兩側(cè)。

(4)較之M6.4地震,M5.6和M5.2地震的PGA觀測值與預(yù)測值離散性更大,造成這些現(xiàn)象的主要原因可能與M5.6和M5.2地震震級較小、臺站獲取記錄較少及臺站分布不均勻有關(guān)。地震動(dòng)參數(shù)衰減關(guān)系是利用觀測的強(qiáng)震記錄回歸得到的,臺站PGA觀測值本身存在較大隨機(jī)性,因此整體上PGA觀測值與兩種衰減關(guān)系的衰減特性一致。

4 加速度反應(yīng)譜分析

加速度反應(yīng)譜是目前國內(nèi)外工程抗震設(shè)計(jì)的重要依據(jù),反映了不同結(jié)構(gòu)周期的單自由度體系在地震力作用下的最大反應(yīng),在地震工程研究和地震動(dòng)區(qū)劃應(yīng)用中有重要作用,計(jì)算和分析加速度反應(yīng)譜有助于進(jìn)一步了解強(qiáng)震動(dòng)記錄的頻域特征[10-11]。

4.1 典型臺站反應(yīng)譜分析

選取距離3次地震震中最近的53YBX臺的加速度記錄進(jìn)行加速度反應(yīng)譜特征分析。該臺距離M6.4、M5.6及M5.2地震震中分別為7.77 km、5.37 km和9.12 km。3次地震中53YBX臺的EW向、NS向及UD向加速度反應(yīng)譜(阻尼比為5%)如圖5所示。由圖5可知:

圖5 3次地震53YBX臺加速度反應(yīng)譜Fig.5 Acceleration response spectra of 53YBX station for the three earthquakes

(1)圖5(a)為3次地震中53YBX臺捕獲記錄的加速度反應(yīng)譜曲線,體現(xiàn)了不同震級地震在類似震中距的情況下對同一臺站反應(yīng)譜的影響?？梢钥闯?不同震級對反應(yīng)譜影響主要表現(xiàn)在0.5 s之后的長周期成分。類似震中距的情況下,震級愈大,長周期分量愈明顯;震級愈小,幅值衰減愈快。

(2)從圖5(a)可以看出,同一震級對反應(yīng)譜的影響主要反映在2 s后的反應(yīng)譜幅值衰減,同一臺站、同一震級的3個(gè)方向反應(yīng)譜形態(tài)整體一致,各分量的最大加速度反應(yīng)譜幅值均在0.05～0.2 s間;各分量在10 s以后的反應(yīng)譜幅值基本衰減至0,其中UD向衰減最快,NS向次之,EW向最慢。

(3)圖5(b)～(d)分別對比了53YBX臺在不同震級、同一方向上的加速度反應(yīng)譜。可以看出:不同震級、類似震中距的情況下,震級越大,臺站同一方向的加速度反應(yīng)譜峰值越大;不同震級、類似震中距的情況下,EW向、NS向及UD向反應(yīng)譜譜形相似,反映了震級對同一臺站場地的卓越周期及反應(yīng)譜形狀的影響有限,但出現(xiàn)這種現(xiàn)象也可能是因?yàn)?次地震震級相差較小。

4.2 加速度反應(yīng)譜平均周期

地震動(dòng)會(huì)對遠(yuǎn)場軟土場地上的高層建筑、大跨橋梁及儲(chǔ)油罐等產(chǎn)生嚴(yán)重的震害,國內(nèi)外專家均對如何評價(jià)遠(yuǎn)場長周期地震動(dòng)進(jìn)行了研究。文獻(xiàn)[12-13]建議采用對地震動(dòng)低頻成分更加敏感的加速度反應(yīng)譜平均周期Tra來評價(jià)遠(yuǎn)場長周期地震動(dòng):

(4)

式中:Ti為5%阻尼比加速度反應(yīng)譜等間距離散周期;Sa(Ti)為Ti對應(yīng)的譜加速度。

根據(jù)式(4)計(jì)算M6.4地震中獲得的80條加速度記錄的加速度反應(yīng)譜平均周期Tra,并繪制Tra擬合曲線圖(圖6)。為對比顯示,圖中離散點(diǎn)為實(shí)際計(jì)算的Tra/2。從圖6中可以看出,在近場區(qū)域(100 km范圍內(nèi))Tra相對較小,這是因?yàn)樵诮鼒鰠^(qū)域,地震動(dòng)高頻成分占主要地位,長周期成分表現(xiàn)不明顯[12];隨著震中距增加,NS向和UD向呈上升趨勢,高頻地震分量在波的傳播過程中能量逐漸消耗,長周期成分成為地震動(dòng)頻譜的主要特征。本次M6.4地震由于缺少350 km以外的遠(yuǎn)場記錄,加速度反應(yīng)譜平均周期的后續(xù)變化趨勢無法再計(jì)算。但整體來看,Tra從反應(yīng)譜的角度比較有效地區(qū)分了近場和遠(yuǎn)場的地震動(dòng)長周期特性,可以作為評價(jià)地震動(dòng)長周期特性的參考。

圖6 M6.4地震加速度反應(yīng)譜平均周期Fig.6 Average period of acceleration response spectrum of the M6.4 earthquake

5 結(jié)論

本文對2021年5月21日云南漾濞M5.6、M6.4及M5.2地震中獲取的120條強(qiáng)震動(dòng)觀測記錄進(jìn)行了處理和初步分析,計(jì)算強(qiáng)震動(dòng)記錄相關(guān)參數(shù),分析了幅值衰減特性及典型臺站加速度反應(yīng)譜,并得出以下結(jié)論:

(1)M5.6地震獲取記錄的PGA介于0.8～338.9 cm/s2間,M6.4地震獲取記錄的PGA介于0.8～719.6 cm/s2間,M5.2地震獲取記錄的PGA介于2.0～190.6 cm/s2間,最大PGA均由53YBX臺獲取。整體上隨震中距增加,各強(qiáng)震動(dòng)臺站的PGA逐漸衰減,直至降為0。

(2)繪制了3次地震的峰值加速度分布圖,發(fā)現(xiàn)其峰值加速度等值線均呈NW-SE方向展布,最大值位于震中附近。其中M6.4地震峰值加速度等值線分布圖與云南省地震局發(fā)布的《云南漾濞6.4級地震烈度圖》基本一致。

(3)將PGA觀測值與兩種衰減關(guān)系進(jìn)行比較,發(fā)現(xiàn)實(shí)際觀測值與俞言祥青藏區(qū)(包括川滇地區(qū))地震動(dòng)衰減關(guān)系預(yù)測曲線吻合更好,而與汪素云中國西部地震動(dòng)衰減關(guān)系預(yù)測值離散性較大。

(4)分析典型臺站的加速度反應(yīng)譜,發(fā)現(xiàn)震級愈大,反應(yīng)譜中長周期分量愈明顯;震級愈小,幅值衰減愈快。

(5)加速度反應(yīng)譜平均周期能比較有效地區(qū)分近場和遠(yuǎn)場的地震動(dòng)長周期特性,可以作為評價(jià)地震動(dòng)長周期特性的參考。

致謝:感謝云南、四川強(qiáng)震動(dòng)觀測人員在3次地震觀測、記錄收集工作中的辛勤付出！