基于HHT地震動分量分離的長周期地震動界定方法

李英民， 趙晨曉， 譚 潛

(1.重慶大學 土木工程學院，重慶 400045；2.山地城鎮建設與新技術教育部重點實驗室,重慶 400045)

隨著大型工程結構體系的快速發展，長周期地震動的潛在破壞能力已不容忽視。如因為超高層建筑物、大型儲液罐、大型橋梁等對長周期地震動都較為敏感，導致長周期地震動對此類大型工程結構體系的影響可能超出了設計時的估計水平。然而，從Hanks[1]提出長周期地震動的概念，到墨西哥地震中長周期震害的爆發，再到目前長周期地震動已經引起廣泛關注為止， 雖然已經認識到了存在有遠場長周期和近斷層脈沖型兩類特殊的長周期地震動[2-3]，但是基于分量分離方法的長周期地震動的區分和界定方法的相關研究卻很少。Baker[4]基于小波理論對地震動進行分解，對脈沖與非脈沖地震動進行了界定；而Farid等[5]采用滑動平均濾波器對地震動進行分解；分析各分量對結構響應的影響，指出小波分析法不能較好反映實際地震動的各分量特征;徐龍軍等[6]則采用了基于seismo signal的數字濾波對地震動進行分解，分析反應譜的相關特點并對規范反應譜進行修正，卻沒有對地震動的區分深入探討；李雪紅等[7]基于β譜的概念提出了長周期地震動的界定指標，此界定方法只停留在長周期地震動表面，不能合理的過濾地震動的高頻分量的影響，并沒有深入到影響長周期地震動的本質長周期分量上來。長周期地震動與地震動長周期分量是兩個不同的概念，前者是一種特殊的地震波形式，后者是地震動中的長周期成分。對于任何地震動，都會有其長周期分量，而長周期地震動中的長周期分量在長周期結構的響應中起著決定性作用。對于長周期分量的研究能更清晰地得到長周期地震動的特性規律，并有利于人造長周期地震動的合成。

本文通過對三種頻譜平均周期的對比評價，選取加速度反應譜平均周期作為長周期分量劃分的指標，采用HHT方法(Hilbert-Huang Transform，HHT)將符合條件的固有模態函數(Intrinsic Mode Funtion，IMF)和殘余分量進行重構形成新的地震動并對重構的地震動進行基線校正后作為長周期分量；在此基礎上，通過對240條地震動的長周期分量進行分析，提出長周期地震動的界定方法以及兩類特殊長周期地震動的劃分指標。

1 HHT基本原理

HHT是由Huang等[8]提出的一種非平穩信號的時頻域分析方法。相比常用模擬信號的處理方法，其徹底擺脫了線性和平穩性的束縛，具有完全的自適應性，并不受Heisenberg測不準原理的制約。更適合用于地震動的分量分離研究。HHT變換主要分為兩個步驟：首先，對信號進行經驗模態分解(Empirical Mode Decomposition，EMD)；然后，對分解后得到的固有模態函數進行Hilbert譜分析(Hilbert Spectral Analysis，HSA)[9]。

EMD分解過程采用的是“篩分”算法，主要分為三個步驟：

步驟1確定地震動信號x(t)所有的局部極值點，并用三次樣條曲線連接所有的局部極大值和局部極小值，分別形成上包絡線xmax(t)和下包絡線xmin(t)

步驟2求出上下包絡線的平均值，記為m11(t)，將原地震動時程曲線x(t)減去該平均值后得到新數據序列h11(t)

m11(t)=[xmax(t)+xmin(t)]/2

(1)

h11(t)=x(t)-m11(t)

(2)

判斷h11(t)是否具備IMF的兩個條件，若不滿足，則將h11(t)作為原始數據重復上述過程，直到得到的新數據序列

h1k(t)=h1(k-1)(t)-m1k(t)

(3)

滿足IMF的兩個條件，這樣就得到了第一個IMF分量c1(t)

c1(t)=h1k(t)

(4)

步驟3從原地震動x(t)中去掉c1(t)得到剩余序列r1(t)

r1(t)=x(t)-c1(t)

(5)

將r1(t)作為新信號序列，按照上述步驟，依次得出第2、第3、…、直至第n個固有模態函數cn(t)，當殘量rn(t)為時間的單調函數或者小于某一預定值時，則EMD分解結束。原始地震動時程曲線則可以表達為各IMF分量與最終殘量的總和，即

(6)

對每一個IMF分量進行HSA，得到相應的Hilbert譜，匯總后得到原始信號的Hilbert譜H(ω,t)

(7)

將H(ω,t)對時間積分，得到Hilbert幅值邊際譜

(8)

將H(ω,t)的平方對頻率積分，得到瞬時能量曲線

(ω,t)2dω

(9)

再將IE(t)對時間積分，得到Hilbert累積能量曲線

(10)

2 地震動的選取

本文選取30條文獻[10-14]中所認為的遠場長周期地震動，以及30條文獻[15-17]中所認為的近斷層長周期地震動來進行研究，所選地震動信息見表1。由于所選60條長周期地震動，基本上都來自于1999年集集地震，且集集地震場地變化多，峰值、震中距涉及范圍廣，具有普遍性。為便于更好的對比分析，從集集地震中另外隨機選取180條水平地震動作為區分長周期地震動的對照組數據，由于篇幅原因文中便不一一列出。

3 長周期分量提取方法

3.1 長周期地震動各頻譜平均周期的分析

平均周期根據其計算所依據的頻譜不同，可以分為傅里葉譜平均周期Tm和阻尼比為5%的加速度反應譜的平均周期Tr，以及Hilbert幅值邊際譜的平均周期Tmh。其計算公式分別見式(11)～式(13)。平均周期參數能夠較好的反應地震動的頻譜成分[18]，而哪一個周期參數能更好的評價長周期地震動還需要進一步探索。

表1 長周期地震動基本信息

(11)

(12)

(13)

式中：fi為經快速傅里葉變換后對應的頻率點；Ci為傅里葉幅值譜中fi所對應的幅值；Ti為阻尼比為5%的加速度反應譜的周期點；Sa為Ti所對應的加速度幅值；fhi為HHT變換后所得幅值邊際譜的頻率點；Hi為Hilbert幅值邊際譜中fhi所對應的幅值。

離散程度和相關程度是判斷一個指標是否合理的重要依據，通過對長周期地震動各頻譜的平均周期的相關性及離散性分析，確定長周期分量的劃分指標。并通過對地震動EMD分解后得到的各IMF分量進行分析，提取出地震動的長周期分量。

離散性和相關性分別用各頻譜平均周期的變異系數及平均相關系數表示，其計算公式分別見式(14)和式(15)。

(14)

式中:Ti為第i條地震動的頻譜的平均周期。

(15)

本文中使用240條地震記錄的各頻譜平均周期的變異系數及平均相關系數結果,如表2、表3所示。由表2可知,對于長周期地震動Tm變異系數最大，離散性最大；對于對照組地震動Tmh離散性最大；而Tr離散性使終保持在中下，具有更好的穩定性。從表3可知，Tm、Tr、Tmh之間都具有較大的相關性，并且Tr更具有代表性。總的來說Tr能夠更為準確地反映地震動高低頻分量的分布情況并且對低頻分量具有更好的辨別能力，這與以往的研究成果也是一致的[19-20]。因此，平均周期Tr更適合作為長周期分量的劃分指標。

3.2 基于EMD的地震動長周期分量提取

EMD能把復雜非平穩信號分解為有限基本模式分量之和，且基于EMD的濾波方法具有良好的頻率選擇性和自適應性的多分辨力[21]。以TCU052-EW加速度時程為例，對其進行EMD過程處理。該地震動分解得到11個IMF分量和1個殘余分量，如圖1所示。隨著EMD分解的進行，所得到的分量頻率也在逐漸降低。根據式(12)計算出各IMF分量和殘余分量的加速度反應譜平均周期的變化也驗證了此結論，如表4所示。

表2地震動平均周期的變異系數

Tab.2Variationcoefficientofaverageperiodsofeachgroup’sgroundmotions

變異系數TmTrTmh遠場長周期地震動0．3120．2980．251近斷層長周期地震動0．3750．2840．310遠、近場長周期地震動0．3750．3120．284對照組地震動0．4450．4490．459

表3地震動平均周期的平均相關系數

Tab.3Averagecorrelationcoefficientofaverageperiodsofeachgroup’sgroundmotions

平均相關系數TmTrTmh遠場長周期地震動0．9190．9350．897近斷層長周期地震動0．8890．8810．853遠、近長周期地震動0．9080．9130．880對照組地震動0．9250．9360．903全部240條地震動0．9370．9530．936

圖1 TCU052-EW加速度時程的EMD分解(gal)

邱立珊的研究表明，在Tr≥2 s以后，傅里葉譜1 Hz以上的高頻部分幅值大幅衰減，速度譜幅值較大的區域分布在較長周期段。同時，加速度反應譜與結構響應和解決工程問題密切相關，考慮到超高層結構自振周期>2.5 s[22]且共振區間的范圍0.75～1.25，Tr≥2 s也十分合理，并具有一定的工程意義。因此，將EMD分解后得到的分量(包括殘余分量)中Tr≥2 s的分量篩選出來并重構。但重構得到的分量，其速度時程具有較大的基線漂移，如圖2(a)所示。進行基線校正后,如圖2(b)所示。將基線校正后的信號作為該地震動的長周期分量。以TCU052-EW為例，將其C6～C11及Res疊加重構并進行基線校正后作為其長周期分量，與原地震動信號的對比,如圖3所示。

表4 TCU052-EW各IMF分量及殘余分量的Tr

(a)

(b)

圖3 TCU052-EW及其長周期分量

3.3 長周期分量與原地震動相關性分析

地震動的長周期分量與地震動的相關性越好，說明長周期分量能夠更好的代表這種地震動的特點。在反映地震動長周期分量的參數里，李雪紅等提出譜指標βl和加速度反應譜平均周期(Tr)可以較好地反映地震動中低頻成分的含量，計算所選240條地震動的長周期分量與原地震動的相關性，相關性系數隨著長周期分量參數的變化關系,如圖4所示。從圖4可知，隨著地震動中低頻含量的增加，其長周期成分與原地震動的相關性越好；并且對于長周期地震動，其長周期分量與原地震動的相關性處在一個相對較高的水平，采用地震動的長周期分量來界定地震動是合適的。

3.4 長周期地震動與長周期分量對結構響應影響的相關性驗證

以地震動ILA056-NS和TCU052-EW為例，兩條地震動及分解后的各分量加速度時程及反應譜,如圖5所示。根據圖形特點認為其分別代表類和諧長周期地震動和脈沖型長周期地震動，其中ORG(Original Ground motions)代表原地震動，LPC(Long Period Component)代表長周分量，SPC(Short Period Component)代表短周期分量。由反應譜可以看出，兩類地震動在長周期段，長周期分量與原地震動的吻合效果是較好的，能反映原地震動對超高層結構的影響。通過Perform3D軟件，對比地震動及其長周期分量與短周期分量對超高層建筑物地震響應的影響，驗證其長周期分量與原地震動的相關性。

圖4 長周期分量與原地震動相關系數

Fig.4 Correlation coefficient of long period component and original ground motion

模型為處于7度0.1g地區的層高為4.5 m，總高247.5 m,共55層的框架-核心筒結構，其自振周期為5.447 s。模型中荷載標準值為：樓面恒載取5 kN/m2，樓面活載選取2.5 kN/m2，基本風壓選取0.4 kN/m2。重力荷載代表值按照100%恒荷載加50%活荷載計算。模型基本信息和平面布置,見表5和圖6。

(a) ILA056-NS

(b) TCU052-EW

圖5 地震動及各分量加速度時程和加速度反應譜

圖6 結構平面布置圖

將兩條地震動加速度峰值均調幅至100 gal，進行超高層結構響應分析。其內力和變形,如表6所示，從中可以看出，兩類長周期地震動及對應的長周期分量對結構基底剪力以及位移響應的影響是相似的，且最大位移發生的時刻也相近，反觀短周期分量與原地震動對結構響應的影響差異極大。特別是脈沖型長周期地震動，其短周期分量對于位移和基底剪力的貢獻遠不及類和諧長周期地震動。產生這種差異的原因是脈沖型地長周期地震動的主要能量集中在低頻脈沖上，高頻分量影響很小，而類諧和地震動能量分布比較均勻。因此長周期分量在長周期地震動對超高層結構的影響中起著決定性的作用，能夠反映甚至代表長周期地震動在長周期結構中的響應特征。進一步驗證了基于長周期分量對長周期地震動界定的合理性。

表6結構頂層最大位移及基底剪力

Tab.6Themaximumvertexdisplacementandbaseshearofthestructure

地震動及分量PGA/gal發生時刻/s頂層最大位移/m發生時刻/s基底剪力/kNILA056?NSORGLPCSPC100．0058．7886．8744．4783．1341．981．3851．3890．32077．4277．4247．07587395299153804TCU052?EWORGLPCSPC100．0068．8055．3133．3434．3333．861．0541．0970．08341．3341．3548．35571225371425716

4 基于長周期成分的多類型長周期地震動的界定

4.1 長周期與常規地震動的界定

基于上述分量的提取方法可以發現，長周期地震動的長周期分量在自振周期較大的結構中起著決定性作用，因此長周期分量的一些特性指標可以作為描述長周期地震動的依據且過濾掉高頻分量的影響，對地震動的界定有更好的區分性。根據以往的研究結果顯示，地震動的頻譜和能量是區分各類型地震動的最具代表性的特性。累積能量比可以很好地從能量角度反映出長周期分量的含量。圖7為地震動及其長周期分量累積能量曲線。根據式(10)計算其長周期分量與原地震動的Hilbert總累積能量ELPC和EORG的比值作為一個參數來區分的地震動，如式(16)所示

(16)

圖7 長周期分量累積能量曲線

同時，PGV/PGA是能反映地震波頻譜特性，Liao等[23]的統計分析表明，當地震動的PGV/PGA的比值≥0.2時地震動的脈沖特性明顯。然而李雪紅等的研究發現不僅僅脈沖型的地震動有此特點，遠場長周期地震動同樣具有相同的特性。根據地震動單一的PGV/PGA區分脈沖型長周期地震動存在明顯的不足。基于上述特點，本文將基于240條地震動的長周期分量的PGV/PGA和rE界定長周期地震動與常規地震動，將其值表示在以rE為橫坐標，以PGV/PGA為縱坐標的圖8中。從圖8可知，長周期地震動與常規地震動具有較為明顯的界限，對數據進行基于MATLAB的邏輯回歸分析，得到長周期地震動界定指標LPGI(Long Period Ground motions Indicator)，如式(17)，且IPGI>0.6的地震動為長周期地震動，IPGI<0.6的地震動為常規地震動

(17)

由上述界定指標反過來觀察根據文獻所選出的60條長周期地震動中，一些由于其長周期分量含量過少或強度較低，筆者認為將其認作長周期地震動欠妥。

圖8 長周期地震動界定的變量散點圖

4.2 多類型的長周期地震動的區分

目前被認為特殊的兩類長周期地震動，分別是近斷層脈沖型長周期地震動和遠場類和諧長周期地震動。兩類長周期地震動及其長周期分量的能量分布特性有明顯的差異。而長周期分量的能量分布更能準確地反映地震動在長周期階段的脈沖特性。脈沖型地震動的能量主要集中在脈沖處，其積累能量的時間較短。因此結合Hilbert累積能量曲線，借助Trifunac等[24]的90%的能量持時的概念，定義長周期分量的90%的能量持時t90%，其含義為地震動長周期分量的歸一化Hilbert累積能量曲線從0.05增長至0.95所經歷的時間。圖9展示了三類長周期地震動長周期分量的累積能量曲線與速度時程的能量持時的特點。

Fig.9 Cumulative energy curve and velocity time history of long period component of TCU052-EW(left)、KAU044-EW(middle)、ILA056-NS(right)

基于“4.1”中LPGI>0.6區分出來的長周期地震動，統計其90%的能量持時和震中距，來區分遠場長周期地震動和近斷層脈沖型長周期地震動，發現t90%≤38 s時長周期分量呈現出較為明顯的脈沖特性，且時間越短脈沖特性越明顯；t90%≥45 s則呈現出類和諧性；而38 s

5 結 論

(1)所選的三類平均周期(Tm、Tr、Tmh)能較準確地描述高低頻分量的分布情況。由相關性與離散性分析發現，Tr具有更好的穩定性和相關性，更適合作為長周期分量的提取依據。

圖10 長周期地震動震中距與t90%關系

Fig.10 Relation between epicenter distance andt90%of long period ground motion

(2)將地震動時程信號進行EMD分解，將Tr≥2 s的分量篩選出來并進行疊加重構基線校正后得到地震動的長周期分量。長周期地震動的長周期分量與原地震動具有良好的相關性。能夠很好的反應長周期地震動對長周期結構的影響，長周期地震動的特性主要決定于長周期分量的特性。

表7兩類長周期地震動劃分指標

Tab.7Classificationindicatoroftwokindsoflongperiodgroundmotions

長周期地震動類型LPGI震中距/kmt90%/s遠場類和諧長周期地震動近斷層脈沖型長周期地震動>0．6≥100≥45≤50≤38

(3)長周期分量的強度和能量供獻決定著地震動的長周期特性是否顯著。分析240條地震動的長周期分量，根據長周期分量的PGV/PGA比值及其與原地震動的積累能量的比值rE，進行邏輯回歸，提出長周期地震動界定指標LPGI，發現LPGI>0.6時為長周期地震動，且區分效果明顯。

(4)不同類型的地震動的能量分布有明顯的差異，提出用長周期地震動長周期分量的歸一化Hilbert累積能量曲線從0.05增長至0.95所經歷的時間t90%來表征長周期地震動能量分布特性。t90%≤38 s時具有較明顯的脈沖性，t90%≥45 s則表現出類和諧性，而38 s

[1] HANKS T C. Strong ground motion of the San Fernando, California, earthquake: ground displacements[J]. Bulletin of the Seismological Society of America, 1975, 65(1): 193-225.

[2] KOKETSU K, MIYAKE H. A seismological overview of long-period ground motion[J]. Journal of Seismology, 2008, 12(2): 133-143.

[3] 徐龍軍, 胡進軍, 謝禮立. 特殊長周期地震動的參數特征研究[J]. 地震工程與工程振動，2008，28(6)：20-27.

XU Longjun, HU Jinjun, XIE Lili. On characteristics of ground motion parameters for special long-period ground motions[J]. Journal of Earthquake Engineering and Engineering Vibration，2008，28(6)：20-27.

[4] 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.

[5] 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.

[6] 徐龍軍, 趙國臣, 謝禮立. 基于分量分離方法的地震動反應譜[J]. 天津大學學報(自然科學與工程技術版)，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.

[7] 李雪紅, 王文科, 吳迪, 等. 長周期地震動的特性分析及界定方法研究[J]. 振動工程學報，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.

[8] HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proceedings of the Royal Society, 1998, A(454): 903-995.

[9] 張郁山. 希爾伯特-黃變換(HHT)與地震動時程的希爾伯特譜[D].北京：中國地震局地球物理研究所，2003.

[10] 張振炫, 陳清軍. 高層建筑結構長周期地震反應的比較研究[J]. 結構工程師，2009，225(4)：78-84.

ZHANG Zhenxuan, CHEN Qingjun. A comparative study on long-period seismic responses for high-rise structures[J]. Structural Engineers，2009，225(4)：78-84.

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

WANG Bo, BAI Guoliang, WANG Chaoqun, 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.

[12] 王博. 多類型長周期地震動作用下高層結構的動力行為與設計方法研究[D].西安：西安建筑科技大學，2013.

[13] 韓建平, 蘇彥兵. 長周期地震動輸入下某超限高層結構反應分析——以汶川地震和日本東北地震記錄為例[J]. 建筑結構，2014，44(18)：59-62.

HAN Jianping, SU Yanbing. Seismic response analysis of an out-of-code high-rise building under long-period ground motion: case studies of Wenchuan earthquake and Japan Tohoku earthquake records[J]. Building Structure，2014，44(18)：59-62.

[14] 陳清軍, 袁偉澤. 基于長周期地震動記錄的SDOF體系能量譜探討[J]. 振動與沖擊，2013，32(10)：36-42.

CHEN Qingjun, YUAN Weize. Energy spectrum of SDOF systerm based on long-period ground motion records[J]. Journal of Vibration and Shock，2013，32(10)：36-42.

[15] 王博, 白國良, 代慧娟. 典型地震動作用下長周期單自由度體系地震反應分析[J]. 振動與沖擊，2013，32(15)：190-196.

WANG Bo, BAI Guoliang, DAI Huijuan. Seismic response analysis of long-period SDOF system under typical ground motions[J]. Journal of Vibration and Shock，2013，32(15)：190-196.

[16] 楊迪雄, 王偉. 近斷層地震動的頻譜周期參數和非平穩特征分析[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.

[17] 楊迪雄, 趙巖, 李剛. 近斷層地震動運動特征對長周期結構地震響應的影響分析[J]. 防災減災工程學報，2007，27(2)：133-140.

YANG Dixiong, ZHAO Yan, LI Gang. Influence analysis of motion characteristics of near-fault ground motions on seismic responses of long-period structures[J]. Journal of Disaster Prevention and Mitigation Engineering，2007，27(2)：133-140.

[18] 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.

[19] 邱立珊. 基于汶川地震記錄的遠場長周期地震動設計反應譜研究[D].重慶：重慶大學，2016.

[20] 杜東升, 王曙光, 劉偉慶, 等. 長周期地震動影響因素及頻譜參數研究[J]. 建筑結構學報，2014，35(增刊1)：1-8.

DU Dongsheng, WANG Shuguang, LIU Weiqing, et al. Study on affection factors and spectral parameters of long period ground motions[J]. Journal of Building Structures，2014，35(Sup1)：1-8.

[21] 程軍圣,于德介,楊宇,等.基于EMD 的濾波方法研究及應用[C]//第八屆全國振動理論及應用學術會議論文集.上海,2003:105-111.

[22] 徐培福, 肖從真, 李建輝. 高層建筑結構自振周期與結構高度關系及合理范圍研究[J]. 土木工程學報，2014，47(2)：1-11.

XU Peifu, XIAO Congzhen, LI Jianhui. Study on relationship between natural vibration periods and heights of structures for high-rise buildings and its reference range[J]. China Civil Engineering Journal, 2014, 47(2): 1-11.

[23] LIAO W I, LOH C H, WAN S. Earthquake responses of RC moment frames subjected to near-fault ground motions[J]. Structural Design of Tall Buildings，2001，10(3)：219-229.

[24] TRIFUNAC M D, BRADY A G. A study on the duration of strong earthquake ground motion[J]. Bulletin of the Seismological Society of America，1975，65(3)：581-626.