譚潛 李英民 楊永斌

摘 要：在實際地震記錄形成過程中存在高、低頻成分地震波的疊加，地震動高、低頻成分的組成是地震動反應(yīng)譜特征的重要反映。采用數(shù)字濾波技術(shù)，以加速度反應(yīng)譜平均周期2 s為截止周期，結(jié)合經(jīng)驗?zāi)B(tài)分解，提取2 s以上的本征模態(tài)函數(shù)及殘值重構(gòu)并基線校正后作為長周期分量，其余本征模態(tài)函數(shù)重組為短周期分量;通過原地震動與分量的相關(guān)性、反應(yīng)譜的離散性對提取方法的有效性進(jìn)行驗證，并與既有方法進(jìn)行對比。結(jié)果表明：該長周期地震動分量提取方法能較好地提取出在長周期和短周期部分與原地震動具有較好的相關(guān)性、地震反應(yīng)譜離散性較小的長、短周期分量，且較好地反映原地震動的譜特征，在長周期地震動分量提取中具有更好的適應(yīng)性。

關(guān)鍵詞：長周期地震動;經(jīng)驗?zāi)B(tài)分解;本征模態(tài)函數(shù);反應(yīng)譜;頻譜參數(shù);長周期分量

中圖分類號：P315.5;P315.9 文獻(xiàn)標(biāo)志碼：A 文章編號：2096-6717（2021）03-0119-09

Abstract： In the actual seismic record formation process， there is a superposition of high and low frequency components of seismic waves. The composition of high and low frequency components of ground motion is an important reflection of the characteristics of the response spectrum of ground motion. This study is based on the digital filtering technology， taking the average period of acceleration response spectrum （Tavg） 2 s as the cut-off period and combining with empirical mode decomposition（EMD）， and it extracted the intrinsic mode function（IMF） and residual component which the cut-off period above 2 seconds as long-period components of long period ground motions（LPGM） after refactoring and baseline correction. The remaining IMFs were reorganized into short-period components; The validity of the extraction method is verified by correlation between seismic records and their components and the discreteness of the original ground motion and the response spectrum. It also compare this method with existing methods. The results show that the extraction method for LPGM components in this paper can better extract the long and short period components of the seismic response spectrum that are less discrete and have a good correlation with the original ground motions in the long and short periods respectively. It better reflects the spectral characteristics of original ground motions and has better adaptability in the LPGM component separation.

Keywords： long-period ground motions; empirical mode decomposition（EMD）; intrinsic mode function（IMF）; response spectrum; spectral parameter; long period component

在實際地震記錄形成過程中，存在各種成分地震波的疊加（圖1）[1]，故地震記錄的主要特性由其主要成分特性來表達(dá)更具代表性。對地震動原始記錄的分析處理是以往研究地震動的主要方向，其忽視了地震動分量分析，而地震記錄中存在長周期分量的事實不能否認(rèn)[2]。在長周期地震動界定中，徐龍軍等[3]以規(guī)準(zhǔn)反應(yīng)譜的離散性分析為基礎(chǔ)，指出地震動的反應(yīng)譜特征與地震記錄的長周期分量是不可忽視的兩個方面。對地震記錄各類分量的提取和分析是確定地震動類型的可靠方法。地震動分量提取方法與信號分析方法息息相關(guān)。學(xué)者們在信號處理方法的基礎(chǔ)上提出了相關(guān)分量提取方法，Khanse等[4]利用傅里葉變換方法確定頻率，結(jié)合巴特沃斯濾波器對地震記錄進(jìn)行分量分離，提取出長周期地震動脈沖分量，進(jìn)而對相應(yīng)的位移反應(yīng)進(jìn)行分析評估;陳紅等[5]采用分?jǐn)?shù)階伽柏變換對地震信號進(jìn)行時頻分析，其中一項重要的應(yīng)用就是利用分?jǐn)?shù)階伽柏變換來進(jìn)行分?jǐn)?shù)域譜分解，同時，也可以提取頻域中的不同尺度和不同方向上的特征，優(yōu)勢明顯;Baker[6]以小波理論為基礎(chǔ)，提取地震記錄分量，從而對脈沖和非脈沖的地震記錄進(jìn)行識別，具有較好的區(qū)分度，但是，不能對地震動分量的特性進(jìn)一步的描述是此方法的不足;Mallat分解與重構(gòu)是一種基于多尺度分析的小波變換分離分量的關(guān)鍵技術(shù)，趙國臣等[7]將其引入到地震記錄分量提取中，認(rèn)為是一種有效的方法;Farid等[8]應(yīng)用滑動平均濾波方法，結(jié)合截止頻率，將原地震記錄分解，提取脈沖型分量，剩下作為剩余分量。將提取的分量輸入結(jié)構(gòu)，探索結(jié)構(gòu)的彈性反應(yīng)。商業(yè)軟件Seismo Signal中自帶數(shù)字濾波器，可以對地震動進(jìn)行濾波，提取相應(yīng)的分量，徐龍軍等[3]利用軟件的此功能提取脈沖分量，并對地震記錄反應(yīng)譜及提取的分量反應(yīng)譜的離散性進(jìn)行分析和校驗。

在經(jīng)驗?zāi)B(tài)分解（empirical mode decomposition， EMD）等信號處理與分量分離技術(shù)在地震動信號處理中應(yīng)用的基礎(chǔ)上，采用合適評價長周期地震動的截止周期，實現(xiàn)了表征長周期地震動特性的長周期分量的提取，驗證了提取分量的合理性與有效性，提出可應(yīng)用于長周期地震動識別及其規(guī)律特性研究的地震動分量提取方法。

1 EMD理論

EMD作為一種新型自適應(yīng)非平穩(wěn)信號處理方法，不僅對線性平穩(wěn)信號處理效果較好，而且對地震波等非線性、非平穩(wěn)信號的分析處理有很好的適用性。本征模態(tài)函數(shù)（IMF）是EMD的前提。EMD方法認(rèn)為信號都是由若干IMF組成，如果IMFs相互重組疊加，便可以重構(gòu)為復(fù)合信號[9]。所以，EMD的重要目標(biāo)就是為了得到IMFs，實現(xiàn)這個目標(biāo)受到兩個條件的限制：首先，峰值點的個數(shù)與函數(shù)穿零數(shù)必須相等或者最多相差一個;其次，局部極大值構(gòu)成的上包絡(luò)線和極小值構(gòu)成的下包絡(luò)線的平均值為零，以IMFs為基礎(chǔ)可展開EMD，其過程如下：

1）基于給定的地震記錄信號x（t），分別標(biāo)記出其所有局部極大值xmax（t）和局部極小值xmin（t）;

2）將所有的局部極大值（xmax（t））點、局部極小值（xmin（t））點分別用3次樣條曲線連接起來，形成上、下包絡(luò)線;

3）對包絡(luò)線上的局部極大值、局部極小值取平均

4） 移除初始地震記錄信號中的平均趨勢

2 截止周期確定

數(shù)字信號處理領(lǐng)域常用數(shù)字濾波技術(shù)調(diào)制頻率分量，文獻(xiàn)[4]采用頻率1.67 Hz作為識別脈沖特征的界限。Rathje等[10-11]通過對地震記錄的頻譜周期參數(shù)對比后發(fā)現(xiàn)，加速度反應(yīng)譜平均周期Tavg（式（5））在反映長周期地震動的低頻特性方面具有無可比擬的優(yōu)勢，且采用頻譜參數(shù)Tavg能將地震記錄的高、低頻分量分布情況更準(zhǔn)確地區(qū)分，低頻分量更容易被識別。杜東升等[12]研究表明，Tavg的變化與震級的變化具有一致性，Tavg是敏感性更好的頻譜參數(shù)，在評價長周期地震動時更具優(yōu)勢。由Tavg計算公式可以知道其物理意義表示的是周期關(guān)于對應(yīng)譜值平方的一個加權(quán)平均值，反映了地震動頻譜在整個計算頻域內(nèi)的分布情況[11]。

式中：Ti為加速度譜對應(yīng)的離散等間隔周期點;Sa（Ti）為周期點的加速度反應(yīng)譜;PGA為地震動加速度時程對應(yīng)的峰值。

為了更好地體現(xiàn)Tavg對地震動特性區(qū)分的實用性，以Tavg在0.4～5.1 s范圍內(nèi)的10條汶川地震記錄為數(shù)據(jù)基礎(chǔ)，結(jié)合地震動的速度時程、歸一化傅里葉幅值譜和速度譜（圖2），分析地震動頻譜特性隨Tavg變化而變化的規(guī)律：1）從速度時程看，地震動的速度時程曲線隨Tavg的增大而越來越稀疏，穿零點越來越少，表示其長周期成分不斷豐富，表明Tavg的變化能在一定程度上反映地震動周期的變化;2）從傅里葉譜（0.1～20 Hz）看，灰色背景的主頻區(qū)域隨Tavg的增大而向低頻方向平移，低頻成分占比逐漸增加，顯示出地震記錄主要頻譜成分隨Tavg的變化而變化的趨勢，以Tavg=2 s為界，1 Hz以下的低頻成分占比迅速增加，1 Hz以上的高頻部分占比迅速減小;3）從速度譜看，譜幅值相對較大的區(qū)段則隨著Tavg的增大向長周期段平移，同樣以Tavg=2 s為界，越大長周期成分越豐富，速度反應(yīng)譜峰值2 s后基本不斷向右移動。可見，Tavg=2 s界限區(qū)分明顯，可作為長周期分量提取的頻譜指標(biāo)。

3 分量提取方法步驟

1）地震記錄的提取。以汶川地震中的長周期地震動chnua370505為例，其加速度、速度、位移時程曲線見圖3（a），采用EMD方法的步驟對原始地震記錄（Original）加速度時程分解，提取地震記錄中的IMFs和殘余分量（R），提取的IMFs以I1、I2、I3…I12表示，見圖3（b），通過EMD的步驟可以知道，最終提取的IMFs和R均需符合EMD中的兩個限制條件。

2）地震動長、短周期分量的重構(gòu)?；诘卣饎觕hnua370505的IMFs求解對應(yīng)的Tavg值，如表1所示;由于地震記錄周期10 s以后不確定性的存在，導(dǎo)致在限定周期范圍時僅取值10以內(nèi)。以Tavg=2 s為界，將提取的2 s及以上的IMFs與R疊加重構(gòu)為長周期分量（long period components， LPC），小于2 s的IMFs重新組合為短周期分量（short period components， SPC）。表1顯示：地震動的前4個IMFs的Tavg小于2 s，疊加組合為初始短周期分量，后面7個IMFs與R的Tavg大于2 s，重新組合為初始長周期分量，組合后的長、短周期分量時程見圖4。根據(jù)EMD中地震信號重構(gòu)公式（4）可知，重構(gòu)后的加速度LPC與SPC應(yīng)符合式（6）、式（7）的范式，形成的長、短周期分量與IMFs疊加的誤差應(yīng)在可接受范圍內(nèi)，式中k表示Tavg≤2 s的IMFs最大分量的下標(biāo)。進(jìn)一步對加速度分量積分，得到速度分量（圖4），為后續(xù)校正做準(zhǔn)備。

3）分量基線漂移的消除。由圖4的加速度與速度時程曲線可看出，重新組合的加速度分量時程無基線漂移，速度時程偏移明顯。與原始地震動進(jìn)行基線校正同樣的目的，為了消除基線偏移對信號處理結(jié)果造成的干擾，保證分析結(jié)果的有效性，對重構(gòu)后的加速度分量時程進(jìn)行基線校正與濾波，形成校正后最終的加速度長周期分量與短周期分量，積分得到速度分量時程，消除基線漂移后的加速度、速度時程如圖5所示。可以看出，基線校正對基線漂移的控制效果明顯，速度分量時程的偏向大幅度減小，校正后的地震動分量能更好地代表實際地震動，可將地震動中的干擾因素最大程度地排除。

4 長周期分量提取方法的有效性驗證

4.1 地震波分量提取

以Chi-Chi地震中的地震記錄ILA056-NS與TCU052-EW驗證本文提出的提取方法的有效性?；诘卣饎佑涗浀牟ㄐ翁攸c，王博等[13]將長周期地震動分為遠(yuǎn)場長周期地震動與近斷層長周期地震動，從而可將ILA056-NS歸類于前者，而TCU052-EW歸類于后者;圖6展示了兩條地震記錄的加速度與速度時程與其分量?？梢钥闯觯涸谠加涗浿校€漂移均未在加速度時程和速度時程中出現(xiàn);觀察加速度、速度分量的時程，基線偏移在加速度分量中沒有出現(xiàn)，而速度分量偏移明顯，是對上述分量提取方法二次校正的必要性的進(jìn)一步確認(rèn)。

4.2 原地震記錄與其長周期分量相關(guān)性

驗證提取分量的代表性和有效性，可以分析地震記錄分量與原地震動的相關(guān)性。圖4與圖6顯示，未校正的地震記錄分量速度時程存在明顯基線偏移;對比圖5中校正后的分量速度時程，漂移消除效果較好。分析校正后的地震動分量速度時程與原地震動的相關(guān)性更可靠。從圖7可看出相關(guān)性趨勢：從地震記錄的波形判斷，兩條地震記錄校正后的LPC與原始地震動吻合較好，SPC與原地震動吻合較差，說明長周期地震動的LPC在很大程度上能夠代表長周期地震動，而SPC則不能代表，表明上述提取分量的代表性與有效性，同時，表2中的相關(guān)性系數(shù)計算值也更進(jìn)一步說明基于速度的LPC與長周期地震動速度時程的相關(guān)性更好。

4.3 地震地面運動反應(yīng)譜的離散性分析

分量提取方法的適用性已經(jīng)通過相關(guān)性驗證做出說明，進(jìn)一步可以通過地震記錄分解前后反應(yīng)譜的變化來驗證。在設(shè)計反應(yīng)譜準(zhǔn)確估計的影響原因中，基于統(tǒng)計分析的地震反應(yīng)譜的離散性不可忽視，

且加速度反應(yīng)譜長周期段的離散性更加突出[3]。圖8對比了兩條原始地震記錄與其長、短周期分量的反應(yīng)譜的吻合程度，ILA056-NS反應(yīng)譜與長周期分量反應(yīng)譜在長周期段整體吻合較好，TCU052-EW吻合稍差，短周期分量與兩條原地震記錄在短周期段吻合，長周期段偏差顯著。此種差異，是由遠(yuǎn)場長周期地震動與近斷層長周期地震動的長周期分量更豐富的分布特征決定的。徐龍軍等[3]以離散性系數(shù)k來評價地震反應(yīng)譜的離散性，其計算公式見式（8）。借用參數(shù)k來評估LPC反應(yīng)譜與原始地震記錄反應(yīng)譜的離散性。離散系數(shù)k與離散性是正相關(guān)關(guān)系?；趦蓷l地震記錄與其分量反應(yīng)譜，計算得到k，見表3。從表中k的值可判斷，長周期分量與原長周期地震動的離散性較小，短周期分量與原長周期地震動的離散性較大，同圖8結(jié)果一致。

式中：S1T為原地震動在周期T時的反應(yīng)譜值;S2T為各分量在周期T時的反應(yīng)譜值;ΔT為反應(yīng)譜的計算周期間距，取為0.02;Ttol為反應(yīng)譜的計算最大周期，取15 s。

5 不同分量提取方法的結(jié)果對比

Baker[6]對脈沖型地震動進(jìn)行分量分離時采用了基于小波分析的分量提取方法，小波分析過程中，母波選擇多貝西四階小波，將地震記錄分解為15個分量，并計算每個分量的Tavg值，其值大于2.0 s的重構(gòu)為長周期分量，小于2 s的組合為短周期分量;Farid Ghahari等[8]采用滑動平均濾波方法提取脈沖型分量。其步驟：1）用短時傅里葉變換計算速度時程的卓越周期Tp;2）采用式（9）結(jié)合卓越周期Tp得到參數(shù)m的值，式（9）中參數(shù)α采用文獻(xiàn)[8]中的經(jīng)驗取值0.25，dt為地震記錄采樣點的時間間隔;3）以m值代入式（10）計算截止頻率fc，fc取倒數(shù)得到截止周期Tc（式（11））;4）利用Tc，結(jié)合滑動平均濾波來識別長、短周期分量。將上述兩種方法提取的分量反應(yīng)譜與以EMD方法提取的分量轉(zhuǎn)換反應(yīng)譜一起與原地震動的反應(yīng)譜進(jìn)行吻合性對比，由圖9可以看出，基于EMD的長周期分量提取方法提取的長周期分量反應(yīng)譜在長周期段與原地震動反應(yīng)譜整體上吻合較好;與基于小波分析的方法相比，在ILA056-NS中，長周期段具有相近的效果，在TCU052-EW中，本文的方法提取效果更好;而與基于滑動平均濾波的方法比，本文方法提取的長周期分量均能更好地預(yù)測地震記錄和反應(yīng)譜的趨勢。

6 結(jié)論

研究了基于EMD的長周期地震動分量提取方法，對方法的有效性進(jìn)行了分析，并與以往的分量提取方法進(jìn)行對比，主要結(jié)論包括：

1）Tavg大于2 s的長周期成分在長周期地震動中起著至關(guān)重要的作用。由于對于大多數(shù)長周期地震記錄，周期高于2 s的速度譜振幅通常比周期低于2 s的速度譜振幅更具優(yōu)勢，2 s以下基本不再出現(xiàn)峰值點。歸一化傅里葉譜的頻率成分偏向于低頻，通過地震動特性分析確定了分量分離的周期參數(shù)（Tavg）及取值，該結(jié)果被建議作為長周期地震動識別程序中的關(guān)鍵參數(shù)。

2）采用經(jīng)驗?zāi)B(tài)分解方法，將地震動波分解成IMFs;結(jié)合截止周期Tavg=2 s，將相應(yīng)的IMFs分別重組為長周期分量、短周期分量并校正，進(jìn)而建立地震動長周期分量提取方法，并提出了分量提取實現(xiàn)的程序。

3）基于原地震動與長、短周期分量的相關(guān)性，地震動反應(yīng)譜與長、短周期分量的反應(yīng)譜的離散性對分量提取方法的有效性進(jìn)行了驗證，并將基于EMD的分量提取方法與基于小波、滑動平均濾波的方法進(jìn)行了比較，驗證了方法的合理性與有效性。提取的長周期分量是長周期地震動識別模型化的數(shù)據(jù)基礎(chǔ)與參數(shù)指標(biāo)來源。

參考文獻(xiàn)：

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

ZHOU F L， CUI H C， SHIGETAKA A B E， 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. （in Chinese）

[2] 王亞勇. 關(guān)于設(shè)計反應(yīng)譜、時程法和能量方法的探討[J]. 建筑結(jié)構(gòu)學(xué)報， 2000， 21（1）： 21-28.

WANG Y Y. A review of seismic response spectra， time history analysis and energy method [J]. Journal of Building Structures， 2000， 21（1）： 21-28. （in Chinese）

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

XU L J， ZHAO G C， XIE L L. Ground motion response spectra through component decomposition method [J]. Journal of Tianjin University （Science and Technology），2013， 46（11）： 1003-1011.（in Chinese）

[4] KHANSE A C， LUI E M. Pulse extraction and displacement response evaluation for long-period ground motions [J]. The IES Journal Part A： Civil & Structural Engineering， 2010， 3（4）： 211-223.

[5] 陳紅， 彭真明， 王峻， 等. 地震信號分?jǐn)?shù)階Gabor變換譜分解方法及應(yīng)用[J]. 地球物理學(xué)報， 2011， 54（3）： 867-873.

CHEN H， PENG Z M， WANG J， et al. Spectral decomposition of seismic signal based on fractional Gabor transform and its application [J]. Chinese Journal of Geophysics， 2011， 54（3）： 867-873.（in Chinese）

[6] BAKER J W. Quantitative classification of near-fault ground motions using wavelet analysis [J]. Bulletin of the Seismological Society of America， 2007， 97（5）： 1486-1501.

[7] 趙國臣， 徐龍軍， 謝禮立. 基于多尺度分析方法的近斷層地震動特性分析[J]. 地球物理學(xué)報， 2013， 56（12）： 4153-4163.

ZHAO G C， XU L J， XIE L L. On near-fault ground motion characteristics through multi-scale method [J]. Chinese Journal of Geophysics， 2013， 56（12）： 4153-4163. （in Chinese）

[8] FARID G S， JAHANKHAH H， GHANNAD M A. 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.

[9] HUANG N E， SHEN Z， LONG S R， et al. The empirical mode decomposition and the Hilbert spectrum for nonlinearand non-stationary time series analysis [J]. Proceedings of the Royal Society of London Series A： Mathematical， Physical and Engineering Sciences， 1998， 1971 （454）： 903-995.

[10] RATHJE E M， ABRAHAMSON N A， BRAY J D. Simplified frequency content estimates of earthquake ground motions [J]. Journal of Geotechnical and Geoenvironmental Engineering， 1998， 124（2）： 150-159.

[11] RATHJE E M， FARAJ F， RUSSELL S， et al. Empirical relationships for frequency content parameters of earthquake ground motions [J]. Earthquake Spectra， 2004， 20（1）： 119-144.

[12] 杜東升， 王曙光， 劉偉慶， 等. 長周期地震動影響因素及頻譜參數(shù)研究[J]. 建筑結(jié)構(gòu)學(xué)報， 2014， 35（Sup1）： 1-8.

DU D S， WANG S G， LIU W Q， et al. Study on affecting factors and spectral parameters of long period ground motions [J]. Journal of Building Structures， 2014， 35（Sup1）： 1-8. （in Chinese）

[13] 王博， 白國良， 王超群， 等. 基于Hilbert-Huang變換的長周期地震動能量時頻分布比較研究[J]. 地震工程與工程振動， 2013， 33（3）： 71-80.

WANG B， BAI G L， WANG C Q， et al. Comparative study on energy time-frequency distribution of long-period ground motions based on Hilbert-Huang Transform [J]. Journal of Earthquake Engineering and Engineering Vibration， 2013， 33（3）： 71-80. （in Chinese）

（編輯 胡玲）