遠(yuǎn)場長周期地震動作用下超高層建筑結(jié)構(gòu)模態(tài)識別

摘要：

遠(yuǎn)場地震動長周期成分顯著，與高層、超高層建筑基本周期相近，容易產(chǎn)生放大效應(yīng)，給高層、超高層建筑結(jié)構(gòu)帶來嚴(yán)重震害。以2022年9月17日21時41分至次日17時39分發(fā)生在臺灣省花蓮縣M6.5、M5.7、M6.9和M5.4地震為例，利用Hilbert-Huang Transform（HHT）方法，初步分析距震中超1 000 km的武漢某超高層建筑結(jié)構(gòu)臺陣對這4次地震的結(jié)構(gòu)動力響應(yīng)觀測數(shù)據(jù)，在對不同時間尺度的本征模態(tài)函數(shù)分量進(jìn)行地震動長周期特性分析的基礎(chǔ)上，重構(gòu)各結(jié)構(gòu)層的地震動，通過功率譜法得到結(jié)構(gòu)層相對地面層的頻率響應(yīng)函數(shù)，并利用復(fù)模態(tài)指數(shù)函數(shù)法獲得結(jié)構(gòu)第一階模態(tài)頻率和振型。對4次地震中該超高層建筑結(jié)構(gòu)模態(tài)頻率和振型的初步分析表明：在花蓮長周期地震動作用下結(jié)構(gòu)未受到損傷。這一認(rèn)識將為遠(yuǎn)場長周期地震動影響下超高層建筑的震后安全性評價提供重要參考。

關(guān)鍵詞：

遠(yuǎn)場長周期地震動; HHT方法; 結(jié)構(gòu)地震響應(yīng)時程重構(gòu); 模態(tài)頻率; 模態(tài)識別

中圖分類號： TU311.3""""" 文獻(xiàn)標(biāo)志碼：A"" 文章編號： 1000-0844（2024）05-1160-12

DOI：10.20000/j.1000-0844.20230714001

Modal identification of super high-rise structures

under far-field long-period ground motions

LUO Yuan1， YANG Jiang1，2，3， FAN Tao1，2，3

（1.Institute of Seismology， CEA， Wuhan 430071， Hubei， China;

2.Wuhan Institute of Seismic Scientific Instruments Co.， Ltd.， Wuhan 430071， Hubei， China;

3.Engineering Technology Research Center for Earthquake Monitoring and Early Warning

Disposal of Major Projects in Hubei Province， Xianning 437000， Hubei， China）

Abstract：

This paper examines the impact of long-period components of far-field ground motion， which can resonate with the natural period of high-rises and super high-rises， potentially leading to significant earthquake damage. The study focuses on the M6.5， M5.7， M6.9， and M5.4 earthquakes that occurred in Hualien County， Taiwan， China， from 21：41 on September 17， 2022， to 17：39 the following day. Using the Hilbert-Huang transform method， the paper preliminarily analyzes the structural dynamic response data of a super high-rise in Wuhan， located more than 1 000 kilometers away from the epicenter. By analyzing the long-period characteristics of the intrinsic mode function components across different time scales， the ground motion of each structural layer was reconstructed. The power spectrum method was used to determine the frequency response function of each structural layer relative to the ground layer. The first modal frequency and vibration mode were identified using the complex-mode Exp-function method. The analysis results indicated that the modal frequencies and vibration modes of the super high-rise structure remained intact under the long-period ground motions from the four earthquakes. This finding provides a crucial reference for post-earthquake safety evaluations of super high-rises in Wuhan when exposed to long-period ground motions.

Keywords：

far-field long-period ground motion; HHT method; time-history reconstruction of structural seismic response; modal frequency; modal identification

0 引言

遠(yuǎn)場長周期地震動是一種特殊類型的地震動，因其顯著的長周期成分與高層、超高層建筑結(jié)構(gòu)基本周期相近而容易產(chǎn)生放大效應(yīng)，使得高層、超高層結(jié)構(gòu)易成為地震潛在發(fā)作區(qū)域內(nèi)的主要受害對象，尤其是大地震產(chǎn)生的長周期（2～10 s）地震動在大盆地中被放大，對城市高層建筑產(chǎn)生了嚴(yán)重威脅［1-2］。如早在1954年的美國Dixie valley地震（MS6.8），距震中300 km 的Sacramento市的儲水池由于與貯水產(chǎn)生共振而造成破壞［3］;2008年中國汶川地震（MW7.9），距震中750 km的西安地區(qū)某高層建筑產(chǎn)生了較為嚴(yán)重的破壞，距震中3 000 km的泰國曼谷市中心一些高樓持續(xù)搖晃了7 min［4］;2011年日本東海地震（M9.0）中，離震源較遠(yuǎn)的東京、名古屋與大阪等城市的超高層建筑都發(fā)生了大幅度、長時間晃動，地震動持續(xù)時間約300 s，其中最典型的是距震源770 km的大阪府政府第二廳舍，上部結(jié)構(gòu)地震反應(yīng)持續(xù)了近10 min，建筑物底部加速度為34 gal，但其頂部最大位移達(dá)到37 cm［5］;1985年墨西哥大地震（M8.1）中，離震源400 km的墨西哥城的高層建筑受到雙重共振的影響而大量損毀［6-8］。這些震害表明，遠(yuǎn)場地震對建筑結(jié)構(gòu)尤其是超高層結(jié)構(gòu)造成的災(zāi)害不可忽視［9-10］ 。這些地震記錄以及震害特點，在分析遠(yuǎn)場長周期地震動對超高層建筑結(jié)構(gòu)的地震響應(yīng)中具有寶貴的借鑒作用。

遠(yuǎn)場地震動對高層建筑產(chǎn)生結(jié)構(gòu)響應(yīng)的是其長周期成分，其頻譜特性是影響長周期結(jié)構(gòu)動力響應(yīng)和破壞程度的重要因素［11］。隨著長周期地震動記錄的不斷積累，科研人員開始關(guān)注表征長周期地震動頻譜特性的周期參數(shù)，并取得了一定的成果，如希爾伯特（Hilbert）邊際譜平均周期能夠有效地反映地震動長周期成分的貢獻(xiàn)，穩(wěn)定性好，可作為遠(yuǎn)場長周期地震動的頻譜特征周期參數(shù)［12］;傅里葉譜平均周期、希爾伯特邊際譜特征周期適宜表征近場地震動長周期成分［13-14］;加速度反應(yīng)譜平均周期是反映長周期地震動頻譜特性的最優(yōu)周期參數(shù)［11］等。在分離長周期地震動的長周期分量方面，目前主要有小波理論、滑動平均濾波器SeismoSignal的數(shù)字濾波、地震動放大系數(shù)β譜及希爾伯特黃變換（Hilbert-Huang Transform，HHT）時頻譜頻率變異系數(shù)等方法［15-19］。同時，基于強(qiáng)震動監(jiān)測臺陣進(jìn)行建筑結(jié)構(gòu)模態(tài)參數(shù)識別，已成為大型結(jié)構(gòu)健康監(jiān)測一種有效而實用的故障診斷和安全檢測方法，獲取的結(jié)構(gòu)模態(tài)參數(shù)在結(jié)構(gòu)動力特性評價、結(jié)構(gòu)振動控制［20］、結(jié)構(gòu)損傷診斷與安全性評價［21］等領(lǐng)域獲得重要應(yīng)用。常用方法有隨機(jī)子空間識別法、稀疏時域法、HHT法、正交多項式擬合法、神經(jīng)網(wǎng)絡(luò)模型及小波變換等［15，22-25］。

在以上研究基礎(chǔ)上，本文以2022年9月17—18日臺灣省花蓮縣M6.5、M5.7、M6.9和M5.4地震作用下武漢市某超高層建筑結(jié)構(gòu)動力響應(yīng)為例，采用經(jīng)驗?zāi)B(tài)分解（Empirical Mode Decomposition，EMD）方法獲得了4次地震動本征模態(tài)函數(shù)（Intrinsic Mode Function，IMF）分量，綜合利用IMF各分量傅里葉譜及其平均周期、Hilbert時頻譜及能量占比等方法，在分析IMF分量頻率特性基礎(chǔ)上重構(gòu)結(jié)構(gòu)地震響應(yīng)時程，再采用互功率譜法對建筑結(jié)構(gòu)的一階模態(tài)頻率及振型進(jìn)行識別，并與數(shù)值模擬結(jié)果進(jìn)行對比分析。

1 地震響應(yīng)臺陣簡介

本文所分析的武漢市高層建筑采用鋼筋混凝土框架剪力墻結(jié)構(gòu)。響應(yīng)臺陣所在建筑1～2層為小開間單元式商店，1層層高5.75 m，2層層高5.40 m，3～4層層高4.20 m，5層層高4.35 m，其余辦公樓層高均為4.50 m;在6層、16層、26層設(shè)避難層，層高5.95 m，屋面標(biāo)高163.25 m?？拐鹪O(shè)防烈度為6度，設(shè)計基本地震加速度值為0.05g，設(shè)計地震分組為第一組。結(jié)構(gòu)臺陣分別在第1層（標(biāo)高）、第13層（56.85 m）、第22層（98.80 m）和第35層（158.75 m）各布設(shè)1個三向加速度傳感器。

2 地震動長周期成分及結(jié)構(gòu)地震響應(yīng)時程重構(gòu)

長周期地震動對建筑結(jié)構(gòu)產(chǎn)生影響的是其長周期分量，因此，分析長周期地震動對結(jié)構(gòu)的影響首先要分析其長周期分量。

結(jié)構(gòu)響應(yīng)觀測臺陣記錄的地震動還包含不同頻率的噪聲，如地球內(nèi)部運動引起的噪聲、不同高度結(jié)構(gòu)樓層環(huán)境噪聲，以及觀測設(shè)備本身的噪聲等。結(jié)構(gòu)地震響應(yīng)時程重構(gòu)最大的特點是對地震動有效降噪，從而使得重構(gòu)結(jié)構(gòu)的地震響應(yīng)時程更加清晰地反映出結(jié)構(gòu)的頻響函數(shù)，進(jìn)而更加準(zhǔn)確地獲取結(jié)構(gòu)的模態(tài)頻率、振型幅值等信息。

2.1 HHT方法

HHT方法具有自適應(yīng)特點，由EMD和Hilbert譜分析兩部分構(gòu)成。

2.1.1 EMD過程

EMD方法主要針對非線性且不平穩(wěn)的數(shù)據(jù)，在時頻域根據(jù)信號自身的時間尺度進(jìn)行自適應(yīng)分析，使信號體現(xiàn)出時頻聚集性。其核心是把原始信號分解成有限個具有瞬時頻率與幅度的本征模態(tài)函數(shù)（分量）IMFs及一個殘余分量（residual）的和。每一個IMF代表了原信號不同頻率段的振蕩變化，反映信號的局部特征；殘余分量則反映信號中的緩慢變化量。EMD方法具體計算過程分為四步：

第一步，提取原始信號S（t）全部極值點，通過3次樣條插值法將局部極大值點連接成上包絡(luò)線Emax（t），將局部極小值點連接成下包絡(luò)線Emin（t），并求得上包絡(luò)線和下包絡(luò)線均值M1（t）：

M1（t）=Emin+Emax2 （1）

第二步，將原始信號S（t）減去包絡(luò)線M1（t），得到中間信號C1（t）：

C1（t）=S（t）-M1（t） （2）

這個過程稱為“篩分”，對原始信號S（t）每經(jīng)過一次“篩分”得到一個新的信號C1（t）。

第三步，判斷C1（t）是否滿足整個時程內(nèi)極值點個數(shù)與過0點個數(shù)相等或最多相差1，上、下包絡(luò)線相對于時間軸局部對稱。如果滿足這兩個條件，該信號就是一個IMF分量;否則，以該信號為基礎(chǔ)，繼續(xù)進(jìn)行“篩分”，直至分解k次后得到的信號滿足IMF條件，即獲得原始信號的第一個IMF分量IMF1：

IMF1=Ck（t）-Mk（t） （3）

第四步，在原始信號S（t）中減去IMF1得到第一次剩余成分R1（t），對R1（t）進(jìn)行上述“篩分”分析，可得到第二個IMF分量IMF2以及第二次剩余分量R2（t）。對R2（t）重復(fù)上面的“篩分”過程，直到剩余信號Ri（t）成為單調(diào)函數(shù)或者序列中只有一個極值點，即完成了信號S（t）的EMD過程。

通過EMD分解構(gòu)成的多個IMF分量蘊含著原始信號從高到低不同頻段的成分，它們反映了信號的特征尺度，代表了信號的內(nèi)在模態(tài)特征，對IMF進(jìn)行HHT變換，即可獲得分量的時頻屬性。

圖1是武漢市超高層第22層結(jié)構(gòu)臺陣記錄到的2022年9月18日花蓮M6.9地震NS向EMD分解結(jié)果。據(jù)圖1可以得出，隨著EMD分解的進(jìn)行，所得到的IMF各階分量IMF1～I(xiàn)MF10的頻率逐漸降低，最終的殘余項為一單調(diào)信號。

通常情況下，分量的振幅越大，說明該頻率段的信號能量越強(qiáng)。就地震動而言，高頻率的第一個IMF分量一般情況下能量都不是最大的。上述10個IMF分量中最大振幅出現(xiàn)在IMF5。

2.1.2 HHT變換

Hilbert變換和Hilbert譜：基于原信號自身的局部均值特征以及時間尺度，HHT首先利用EMD方法將信號按頻率由高到低篩分成一系列IMF分量的和，再對每個IMF分量進(jìn)行Hilbert變換獲得各自的瞬時頻率和瞬時幅值，使每個IMF分量的瞬時頻率和瞬時幅值有了真實的物理意義；再將所有IMF分量的瞬時頻譜綜合，從而得到信號的Hilbert譜。Hilbert譜表示了信號完整的時間-頻率分布［26-27］。

2.2 地震動長周期分量提取

本文采用傅里葉譜平均周期法對EMD各階分量IMF進(jìn)行周期判別：

Tm=∑iC2i1fi∑iC2i （4）

式中：Ci為傅里葉譜幅度值;fi為對應(yīng)的傅里葉譜離散頻率。表1為根據(jù)式（4）計算的M6.9地震各結(jié)構(gòu)層NS向IMF分量的傅里葉譜平均周期。在實際應(yīng)用中，有關(guān)文獻(xiàn)對式（4）約定了不同的截止頻率，如0.1 Hz≤fi≤25 Hz［11］，0.25 Hz≤fi≤20 Hz［28］。

為了比較不同截止頻率對計算結(jié)果的影響，本文分別按照文獻(xiàn)［10］和文獻(xiàn)［29］的截止頻率計算了M6.9地震各結(jié)構(gòu)層NS向傅里葉譜平均周期（表1后面兩組）。據(jù)表1可得出，各結(jié)構(gòu)層地震動IMF分量的傅里葉譜平均周期均小于10 s，大部分IMF的平均周期小于5 s，且同一結(jié)構(gòu)層連續(xù)幾個分量的周期值十分接近，這一結(jié)果沒能有效地反映出EMD分解獲得的IMF分量的頻率漸變分布特性。通過多種截止頻率的嘗試，將最低截止頻率調(diào)整為0.02 Hz，獲得了各IMF分量較為理想的傅里葉平均周期，這可能與長周期地震動具有較豐富的低頻成分、加速度峰值較小、長周期段幅值較大及頻率非平穩(wěn)特性顯著等特點有關(guān)［30-31］。因此，本文將截止頻率調(diào)整為0.02 Hz≤fi≤20 Hz，各IMF分量傅里葉譜平均周期計算結(jié)果列于表1第一組第3列。從這一列傅里葉譜平均周期值的分布情況可以得出，根據(jù)本文的截止頻率獲得的各IMF分量的這一組傅里葉譜平均周期較好地反映了地震動頻率特性的復(fù)雜性。

由于地球介質(zhì)的影響，在實際地震記錄形成過程中存在各種成分地震波的疊加［4］，尤其是遠(yuǎn)場長周期地震動，經(jīng)過長距離的傳播，地震波變得十分復(fù)雜。因此，僅憑地震動頻率分量的傅里葉譜平均周期難以反映出地震動不同頻率的周期特性。為進(jìn)一步分析IMF周期特性，首先計算各IMF的傅里葉譜，圖2（a）為圖1中10個IMF分量，圖2（b）為對應(yīng)的各IMF分量的傅里葉譜。圖2（b）顯示出兩個明顯的特征，一是一個單分量信號可能被分解到多個IMF中；二是一個IMF中可能存在多個單分量信號的成分。這進(jìn)一步表明，由于地震波的復(fù)雜性，僅憑一種方法難以較好地反映出地震動的周期分布信息。

為清楚地對比不同頻率段的能量大小，同時對比不同頻段信號的能量隨時間的變化規(guī)律，對圖2（a）各IMF進(jìn)行HHT變換，結(jié)果如圖3所示。圖3（a）反映的是IMF頻率在時間軸上的變化情況，即頻率時程。從圖中可得出，IMF的頻率均在2 Hz以內(nèi)，大于0.7 Hz的頻率主要分布在40～70 s及105～115 s時間段內(nèi)，結(jié)合圖2可以得出，這個頻段主要以噪聲為主。

根據(jù)上文的分析，在圖2（a）這組分量中IMF5的頻率成分振幅最大;進(jìn)一步分析發(fā)現(xiàn)，在90～110 s間有一個主要峰值成分，與之對應(yīng)的圖2（b）傅里葉譜也顯示IMF5是這組譜中幅度值最大的，最大譜峰頻率為0.28 Hz，說明IMF5是這組分量中的主要成分。同時，在Hilbert時頻譜中在80～110 s時間段內(nèi)也有一個頻率約為0.28 Hz、譜幅值約為1的顯著信號［圖3（a）］。圖3（b）顯示，IMF5的能量占比達(dá)到63.3%，進(jìn)一步說明IMF5是這組成分中的主要分量。根據(jù)HHT理論的地震動能量特性分析［18］，IMF3～I(xiàn)MF6的能量之和占總能量的91.9%［圖3（b）］，是這組分量中的主要成分。

2.3 結(jié)構(gòu)地震響應(yīng)時程重構(gòu)

根據(jù)上節(jié)的分析，結(jié)合IMF傅里葉譜平均周期、傅里葉譜、Hilbert時頻譜及其瞬時能量等信息，選取第22結(jié)構(gòu)層NS向地震動IMF的C3～C6分量對該結(jié)構(gòu)地震響應(yīng)時程進(jìn)行了重構(gòu)（表1粗體字），表1中同時列出了M6.9地震其他結(jié)構(gòu)層的重構(gòu)信息。圖4為根據(jù)重構(gòu)與原數(shù)據(jù)繪制的套合圖。圖4顯示，重構(gòu)的結(jié)構(gòu)地震響應(yīng)時程沒有發(fā)生基線偏移情況，與原數(shù)據(jù)具有良好的一致性。

綜合分析M6.5、M5.7、M5.4等各結(jié)構(gòu)地震響應(yīng)時程的重構(gòu)信息得出，C3～C6分量的頻率為原信號的優(yōu)勢頻率段。

為對比分析重構(gòu)效果，以M6.9主震為例，繪制了結(jié)構(gòu)層NS向各樓層頻響函數(shù)和復(fù)模態(tài)指數(shù)函數(shù)（圖5）。圖5（a）顯示，結(jié)構(gòu)層原數(shù)據(jù)的頻響函數(shù)和復(fù)模態(tài)指數(shù)函數(shù)的模態(tài)信息較亂，尤其是結(jié)構(gòu)的一階模態(tài)，無法判讀模態(tài)頻率和振型等信息。而根據(jù)重構(gòu)的結(jié)構(gòu)地震響應(yīng)時程得到的頻響函數(shù)和復(fù)模態(tài)指數(shù)函數(shù)［圖5（b）］十分清晰地顯示出結(jié)構(gòu)的模態(tài)信息，可以準(zhǔn)確地獲得結(jié)構(gòu)的一階模態(tài)頻率和振型幅值等信息。

2.4 頻率響應(yīng)函數(shù)

結(jié)構(gòu)層加速度傅里葉幅度譜包含地震動頻率成分和結(jié)構(gòu)信息，通過譜比可以得到結(jié)構(gòu)的頻率響應(yīng)函數(shù)（簡稱頻響函數(shù)），利用功率譜法可獲得結(jié)構(gòu)樓層相對地面首層的頻響函數(shù)［29］。

為獲得其他結(jié)構(gòu)層地震動加速度相對于結(jié)構(gòu)首層地震動加速度的頻響函數(shù)，首先計算結(jié)構(gòu)首層地震動加速度自相關(guān)函數(shù)及其他結(jié)構(gòu)層地震動加速度與結(jié)構(gòu)首層地震動加速度的互相關(guān)函數(shù);然后利用傅里葉變換分別計算自相關(guān)函數(shù)和各互相關(guān)函數(shù)的傅里葉自功率譜和互功率譜，通過式（5）構(gòu)建頻響函數(shù)：

Hn1（f）=Pa1an（f）Pa1a1（f） （5）

式中：Pa1an（f）為根據(jù)第n結(jié)構(gòu)層地震動加速度an與結(jié)構(gòu)首層地震動加速度a1構(gòu)建的互相關(guān)函數(shù)求得的互功率譜;Pa1a1（f）為根據(jù)結(jié)構(gòu)首層地震動加速度構(gòu)建的自相關(guān)函數(shù)求得的自功率譜;Hn1（f）為二者的頻響函數(shù)。

根據(jù)式（5）構(gòu)建的主震M6.9地震結(jié)構(gòu)樓層NS向頻響函數(shù)如圖6（c）中（Ⅲ）所示。從圖中可看出，該結(jié)構(gòu)頻響函數(shù)一階模態(tài)十分明顯，頻率值為0.242 Hz。

2.5 模態(tài)參數(shù)識別

模態(tài)參數(shù)的識別可以分為頻域法和時域法兩類。頻域法中具有代表性的有峰值檢測法、復(fù)模態(tài)指數(shù)函數(shù)法等;時域法主要有隨機(jī)子空間法、隨機(jī)減量法與ERA相結(jié)合的方法等。兩類方法中，頻域識別方法概念直觀、結(jié)論可靠。本文中采用復(fù)模態(tài)指數(shù)函數(shù)法進(jìn)行模態(tài)頻率識別。

對式（5）頻響函數(shù)矩陣H通過奇異值分解可得到復(fù)模態(tài)指數(shù)函數(shù)［32］：

H=U∑VT （6）

式中：H均為頻響函數(shù)矩陣;U、V均為分解的特征值向量酉矩陣，其中V包含模態(tài)振型信息;∑為對角陣。

對頻響函數(shù)矩陣進(jìn)行奇異值分解，可得到復(fù)模態(tài)指數(shù)函數(shù)，M6.9地震NS向結(jié)構(gòu)復(fù)模態(tài)指數(shù)函數(shù)如圖7（a）所示。從該圖可得出，結(jié)構(gòu)的第一階模態(tài)特征十分突出，對應(yīng)的頻率值為0.242 Hz，與頻響函數(shù)獲得的結(jié)果一致。

根據(jù)以上的分析思路與過程，本文還分析了主震M6.9 EW向、前震M6.5、M5.7以及余震M5.4地震，分析過程在此不再贅述。

3 模態(tài)初步分析

根據(jù)上文介紹的分析思路，利用式（5）得到M6.5、M5.7、M6.9和M5.4地震結(jié)構(gòu)NS向和EW向各樓層頻響函數(shù)，如圖6所示。從圖6可得出，前震、主震和余震時結(jié)構(gòu)反應(yīng)能量均主要以第一階頻率為主。NS和EW兩個方向第一階頻響函數(shù)十分明顯，頻率值分別為0.23～0.30 Hz和0.27～0.30 Hz。

根據(jù)該超高層建筑抗震設(shè)防可行性論證報告的數(shù)據(jù)顯示，該結(jié)構(gòu)一階模態(tài)頻率數(shù)值模擬結(jié)果NS向為0.23 Hz、EW向為0.214 Hz，與本文獲得的結(jié)果基本一致。

利用式（6）對每個頻率點采用復(fù)模態(tài)指數(shù)函數(shù)法得到的結(jié)構(gòu)NS和EW向第一階模態(tài)頻率，如圖7所示。據(jù)圖7可得出：4次地震的復(fù)模態(tài)指數(shù)十分接近，其中NS向和EW向第一階模態(tài)頻率分別為0.24～0.30 Hz和0.27～0.30 Hz，根據(jù)文獻(xiàn)［29］中相對誤差計算方法可得，NS向和EW向相對誤差分別為11.11%和5.26%。

根據(jù)奇異值分解得到的V矩陣提取模態(tài)振型，得到NS向和EW向第一階振型如圖8所示。由圖8顯示，4次地震作用下結(jié)構(gòu)NS向和EW向第一階的模態(tài)振型都很接近。

為進(jìn)一步研判結(jié)構(gòu)是否受到地震影響，結(jié)合武漢某超限高層建筑工程抗震設(shè)防可行性論證報告提供的相關(guān)數(shù)據(jù)計算了4次地震作用下結(jié)構(gòu)的模態(tài)置信度（Modal Assurance Criterion，MAC）（表2）。

表2顯示，4次地震作用下NS向和EW向第一階模態(tài)的MAC分別在0.95和0.94以上，兩個方向MAC最大互差分別為0.024和0.013，表明4次地震下兩個方向的同一振型均具有良好的一致性，同一振型高度相關(guān)。依據(jù)基于模態(tài)振型的結(jié)構(gòu)損傷識別方法，并考慮NS和EW向的第一階模態(tài)頻率，可以判斷結(jié)構(gòu)沒有受到花蓮地震的影響。

4 結(jié)論與討論

本文利用傅里葉譜平均周期，結(jié)合傅里葉譜、Hilbert時頻譜及能量占比等方法，綜合分析2022年9月17—18日臺灣省花蓮縣M6.5、M5.7、M6.9和M5.4地震遠(yuǎn)場長周期地震動的長周期成分，在此基礎(chǔ)上重構(gòu)了結(jié)構(gòu)地震響應(yīng)時程，獲得了武漢某超高層建筑結(jié)構(gòu)4次地震響應(yīng)的一階模態(tài)頻率和模態(tài)振型。綜合分析認(rèn)為：

（1） 不同截止頻率對傅里葉譜平均周期計算結(jié)果的影響較大。本文設(shè)置的0.02～20 Hz的截止頻率獲得的IMF分量的傅里葉譜平均周期較好地反映了地震動頻率分布特性。傅里葉譜平均周期約為5 s的長周期成分是地震動的主要分量，這對分析遠(yuǎn)場地震動的長周期成分，研判其對高層、超高層建筑結(jié)構(gòu)的影響程度起著至關(guān)重要的作用。

（2） 結(jié)合傅里葉譜、Hilbert時頻譜及其瞬時能量譜綜合分析了IMF特性，在此基礎(chǔ)上識別并提取了長周期IMF分量，通過結(jié)構(gòu)地震響應(yīng)時程重構(gòu)，構(gòu)建了十分清晰的結(jié)構(gòu)頻響函數(shù)和復(fù)模態(tài)指數(shù)函數(shù)，獲得了結(jié)構(gòu)的一階模態(tài)頻率和模態(tài)振型，其中NS向第一階模態(tài)頻率相對誤差為11.11%，EW向為5.26%。

（3） 4次地震作用下結(jié)構(gòu)NS向和EW向同一振型高度相關(guān)，綜合考慮一階模態(tài)頻率和模態(tài)振型認(rèn)為，在花蓮4次地震作用下該超高層建筑結(jié)構(gòu)沒有受到損傷。這一認(rèn)識可為遠(yuǎn)場長周期地震動影響下武漢市超高層建筑的震后安全性評價提供參考。

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

［1］ FURUMURA T，OISHI Y.An early forecast of long-period ground motions of large earthquakes based on deep learning［J］.Geophysical Research Letters，2023，50（6）：e2022GL101774.

［2］ 王亞楠，杜永峰，胡高興，等.遠(yuǎn)場長周期地震下基于規(guī)范設(shè)計隔震結(jié)構(gòu)的抗震性能研究［J］.地震工程學(xué)報，2018，40（2）：219-224.

WANG Ya'nan，DU Yongfeng，HU Gaoxing，et al.Seismic performance of base-isolated structures under far-field long-period ground motion based on seismic codes［J］.China Earthquake Engineering Journal，2018，40（2）：219-224.

［3］ 李春鋒，張旸.長周期地震動衰減關(guān)系研究的迫切性［J］.地震地磁觀測與研究，2006，27（3）：1-8.

LI Chunfeng，ZHANG Yang.Urgency of the study of long-period ground motion［J］.Seismological and Geomagnetic Observation and Research，2006，27（3）：1-8.

［4］ 楊偉林，朱升初，洪海春，等.汶川地震遠(yuǎn)場地震動特征及其對長周期結(jié)構(gòu)影響的分析［J］.防災(zāi)減災(zāi)工程學(xué)報，2009，29（4）：473-478.

YANG Weilin，ZHU Shengchu，HONG Haichun，et al.Characteristics of far-field ground motion of Wenchuan earthquake and the effect on long-periodic structures［J］.Journal of Disaster Prevention and Mitigation Engineering，2009，29（4）：473-478.

［5］ 周福霖，崔鴻超，安部重孝，等.東日本大地震災(zāi)害考察報告［J］.建筑結(jié)構(gòu)，2012，42（4）：1-20.

ZHOU Fulin，CUI Hongchao，ABE Shigetaka，et al.Inspection report of the disaster of the East Japan earthquake by Sino-Japanese Joint Mission［J］.Building Structure，2012，42（4）：1-20.

［6］ BECK J L，HALL J F.Factors contributing to the catastrophe in Mexico City during the earthquake of September 19，1985［J］.Geophysical Research Letters，1986，13（6）：593-596.

［7］ AGUILAR J，JUREZ H，ORTEGA R，et al.The Mexico earthquake of September 19，1985：statistics of damage and of retrofitting techniques in reinforced concrete buildings affected by the 1985 earthquake［J］.Earthquake Spectra，1989，5（1）：145-151.

［8］ ANDERSON J G，BODIN P，BRUNE J N，et al.Strong ground motion from the Michoacan，Mexico，earthquake［J］.Science，1986，233（4768）：1043-1049.

［9］ 杜東升，宋寶璽，許偉志，等.高層鋼結(jié)構(gòu)考慮長周期地震動的減震加固研究［J］.工程力學(xué)，2020，37（7）：189-200.

DU Dongsheng，SONG Baoxi，XU Weizhi，et al.Seismic retrofit of a high-rise steel structure considering long-period and long-duration ground motions［J］.Engineering Mechanics，2020，37（7）：189-200.

［10］ 張令心，張繼文.近遠(yuǎn)場地震動及其地震影響分析［J］.土木建筑與環(huán)境工程，2010，32（2）（增刊1）：84-86.

ZHANG Lingxin，ZHANG Jiwen.Analysis of near far field vibration and its seismic effects［J］.Journal of Civil amp; Environmental Engineering，2010，32（2）（Suppl01）：84-86.

［11］ 蔣連接，白國良，黃偉強(qiáng)，等.長周期地震動頻譜特征周期研究［J］.地震工程與工程振動，2022，42（1）：222-230.

JIANG Lianjie，BAI Guoliang，HUANG Weiqiang，et al.Study on spectrum characteristic period of long-period ground motions［J］.Earthquake Engineering and Engineering Dynamics，2022，42（1）：222-230.

［12］ 王博，劉伯權(quán)，吳濤，等.遠(yuǎn)場長周期地震動頻譜特征周期與強(qiáng)度指標(biāo)研究［J］.振動與沖擊，2018，37（7）：211-219.

WANG Bo，LIU Boquan，WU Tao，et al.Characteristic periods and intensity indexes for far-field long-period ground motion frequency spectra［J］.Journal of Vibration and Shock，2018，37（7）：211-219.

［13］ YANG D X，WANG W.Nonlocal period parameters of frequency content characterization for near-fault groundmotions［J］.Earthquake Engineering amp; Structural Dynamics，2012，41（13）：1793-1811.

［14］ 楊迪雄，王偉.近斷層地震動的頻譜周期參數(shù)和非平穩(wěn)特征分析［J］.地震工程與工程振動，2009，29（5）：26-35.

YANG Dixiong，WANG Wei.Analysis of frequency content parameters and non-stationary characteristic of near-fault ground motions［J］.Journal of Earthquake Engineering and Engineering Vibration，2009，29（5）：26-35.

［15］ BAKER J W.Quantitative classification of near-fault ground motions using waveletanalysis［J］.Bulletin of the Seismological Society of America，2007，97（5）：1486-1501.

［16］ FARIDGHAHARI S，JAHANKHAH H，ALI GHANNAD M.Study on elastic response of structures to near-fault ground motions through record decomposition［J］.Soil Dynamics and Earthquake Engineering，2010，30（7）：536-546.

［17］ 徐龍軍，趙國臣，謝禮立.基于分量分離方法的地震動反應(yīng)譜［J］.天津大學(xué)學(xué)報（自然科學(xué)與工程技術(shù)版，2013，46（11）：1003-1011.

XU Longjun，ZHAO Guochen，XIE Lili.Ground motion response spectra through component decomposition method［J］.Journal of Tianjin University （Science and Technology），2013，46（11）：1003-1011.

［18］ 李雪紅，王文科，吳迪，等.長周期地震動的特性分析及界定方法研究［J］.振動工程學(xué)報，2014，27（5）：685-692.

LI Xuehong，WANG Wenke，WU Di，et al.The bounded method and characteristics analysis for long-period ground motions［J］.Journal of Vibration Engineering，2014，27（5）：685-692.

［19］ 唐俊華.采用HHT方法分析近場及遠(yuǎn)場長周期地震動的特點［J］.建筑技術(shù)開發(fā)，2020，47（12）：151-152.

TANG Junhua.Analysis of characteristics of near-field and far-field long-period ground motions using HHT method［J］.Building Technology Development，2020，47（12）：151-152.

［20］ 施袁鋒，朱正言，陳鵬，等.基于EM算法和模態(tài)形式的狀態(tài)空間模型自降階工作模態(tài)分析［J］.工程力學(xué)，2021，38（9）：15-25.

SHI Yuanfeng，ZHU Zhengyan，CHEN Peng，et al.Operational modal analysis using em algorithm and modal-form state-space model with auto model order reduction［J］.Engineering Mechanics，2021，38（9）：15-25.

［21］ 何定橋.建筑結(jié)構(gòu)地震響應(yīng)監(jiān)測和風(fēng)險評估算法研究［D］.北京：清華大學(xué)，2021.

HE Dingqiao.Structural seismic response monitoring and risk assessment algorithm study［D］.Beijing：Tsinghua University，2021.

［22］ 張鳴祥，廖海飛，王建國，等.大跨度橋梁顫振時域的直接分析法［J］.合肥工業(yè)大學(xué)學(xué)報（自然科學(xué)版），2022，45（4）：487-492.

ZHANG Mingxiang，LIAO Haifei，WANG Jianguo，et al.Time domain flutter direct analysis of long-span bridges［J］.Journal of Hefei University of Technology （Natural Science），2022，45（4）：487-492.

［23］ 孫猛猛.非平穩(wěn)環(huán)境激勵下高層建筑結(jié)構(gòu)模態(tài)參數(shù)識別［D］.武漢：武漢理工大學(xué)，2019.

SUN Mengmeng.Modal parameter identification of high-rise building under non-stationary environmental excitation［D］.Wuhan：Wuhan University of Technology，2019.

［24］ 許斌，龔安蘇，賀佳，等.基于神經(jīng)網(wǎng)絡(luò)模型的結(jié)構(gòu)參數(shù)提取新方法［J］.工程力學(xué)，2011，28（4）：35-41.

XU Bin，GONG Ansu，HE Jia，et al.A novel structural parameter extraction methodology with neural network based nonparametric model［J］.Engineering Mechanics，2011，28（4）：35-41.

［25］ 代煜，張建勛.基于小波變換和反向傳播網(wǎng)絡(luò)的模態(tài)參數(shù)辨識［J］.振動與沖擊，2012，31（3）：55-59.

DAI Yu，ZHANG Jianxun.Identification of modal parameters based on wavelet transform and back propagation network［J］.Journal of Vibration and Shock，2012，31（3）：55-59.

［26］ HUANG N E，SHEN Z，LONG S R，et al.The empirical mode decomposition and hilbert spectrum for nonlinear and nonstationary time series analysis［J］.Proceedings of the Royal Society，1998，A454：903- 997.

［27］ 吳琛，周瑞忠.Hilbert-Huang變換在提取地震信號動力特性中的應(yīng)用［J］.地震工程與工程振動，2006，26（5）：41-46.

WU Chen，ZHOU Ruizhong.Application of Hilbert-Huang transform in extracting dynamic properties of seismic signals［J］.Earthquake Engineering and Engineering Vibration，2006，26（5）：41-46.

［28］ RATHJE E M，F(xiàn)ARAJ F，RUSSELL S，et al.Empirical relationships for frequency content parameters of earthquake ground motions［J］.Earthquake Spectra，2004，20（1）：119-144.

［29］ 陶冬旺，劉泉，馬強(qiáng)，等.云南漾濞地震大理某高層建筑結(jié)構(gòu)地震響應(yīng)觀測數(shù)據(jù)初步分析［J］.世界地震工程，2021，37（3）：19-30.

TAO Dongwang，LIU Quan，MA Qiang，et al.Preliminary analysis of structural seismic response of a tall building in Dali recorded in Yunnan Yangbi earthquakes［J］.World Earthquake Engineering，2021，37（3）：19-30.

［30］ 張振炫，陳清軍.高層建筑結(jié)構(gòu)長周期地震反應(yīng)的比較研究［J］.結(jié)構(gòu)工程師，2009，25（4）：78-84.

ZHANG Zhenxuan，CHEN Qingjun.A comparative study on long-period seismic responses for high-rise structures［J］.Structural Engineers，2009，25（4）：78-84.

［31］ 譚潛，李英民，向淵明，等.遠(yuǎn)場長周期地震動特征［J］.地震工程與工程振動，2019，39（5）：177-188.

TAN Qian，LI Yingmin，XIANG Yuanming，et al.Study on characteristics of far-field long period ground motion［J］.Earthquake Engineering and Engineering Dynamics，2019，39（5）：177-188.

［32］ SHIH C Y，TSUEI Y G，ALLEMANG R J，et al.Complex mode indication function and its applications to spatial domain parameter estimation［J］.Mechanical Systems and Signal Processing，1988，2（4）：367-377.

（本文編輯：張向紅）