譚潛 李英民 楊永斌

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

關鍵詞：長周期地震動;經驗模態(tài)分解;本征模態(tài)函數(shù);反應譜;頻譜參數(shù);長周期分量

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

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

1 EMD理論

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

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

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

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

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

2 截止周期確定

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

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

為了更好地體現(xiàn)Tavg對地震動特性區(qū)分的實用性，以Tavg在0.4～5.1 s范圍內的10條汶川地震記錄為數(shù)據(jù)基礎，結合地震動的速度時程、歸一化傅里葉幅值譜和速度譜（圖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為界，越大長周期成分越豐富，速度反應譜峰值2 s后基本不斷向右移動?？梢?，Tavg=2 s界限區(qū)分明顯，可作為長周期分量提取的頻譜指標。

3 分量提取方法步驟

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

2）地震動長、短周期分量的重構?；诘卣饎觕hnua370505的IMFs求解對應的Tavg值，如表1所示;由于地震記錄周期10 s以后不確定性的存在，導致在限定周期范圍時僅取值10以內。以Tavg=2 s為界，將提取的2 s及以上的IMFs與R疊加重構為長周期分量（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中地震信號重構公式（4）可知，重構后的加速度LPC與SPC應符合式（6）、式（7）的范式，形成的長、短周期分量與IMFs疊加的誤差應在可接受范圍內，式中k表示Tavg≤2 s的IMFs最大分量的下標。進一步對加速度分量積分，得到速度分量（圖4），為后續(xù)校正做準備。

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

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

4.1 地震波分量提取

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

4.2 原地震記錄與其長周期分量相關性

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

4.3 地震地面運動反應譜的離散性分析

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

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

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

5 不同分量提取方法的結果對比

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

6 結論

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

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

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

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

參考文獻：

[1] 周福霖， 崔鴻超， 安部重孝， 等. 東日本大地震災害考察報告[J]. 建筑結構， 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] 王亞勇. 關于設計反應譜、時程法和能量方法的探討[J]. 建筑結構學報， 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] 徐龍軍， 趙國臣， 謝禮立. 基于分量分離方法的地震動反應譜[J]. 天津大學學報（自然科學與工程技術版）， 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] 陳紅， 彭真明， 王峻， 等. 地震信號分數(shù)階Gabor變換譜分解方法及應用[J]. 地球物理學報， 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]. 地球物理學報， 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]. 建筑結構學報， 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）

（編輯 胡玲）