近斷層長周期脈沖型地震動對豎向反應譜的影響研究

江 輝，李新樂，竇慧娟(.北京交通大學 土木建筑工程學院，北京 00044;.大連民族學院 土木建筑工程學院，大連 6600)

近斷層長周期脈沖型地震動對豎向反應譜的影響研究

江 輝1，李新樂2，竇慧娟2
(1.北京交通大學 土木建筑工程學院，北京 100044;2.大連民族學院 土木建筑工程學院，大連 116600)

長周期大幅值脈沖運動是近斷層地震動的重要特征，也是引起結構震害的重要因素。在既有的地震動數據庫基礎上，篩選建立了長周期脈沖型記錄的統計樣本庫，并構建了地震動特征參數的回歸關系式。針對豎向反應譜，深入分析了豎向分量與水平向分量的加速度反應譜比值隨脈沖周期和結構自振周期的變化規律，結果表明反應譜比值顯著大于現行規范取值，進而采用分段線性化方法，建立了豎向反應譜比值隨脈沖周期、自振周期變化的簡化關系式，并對影響豎向反應譜曲線各拐點周期的特征比值進行統計分析，擬合得到了其隨震級和斷層距的衰減模型。

近斷層地震;長周期脈沖;豎向反應譜;反應譜比;特征比值;回歸

工程地震學的研究表明，靠近發震斷層的地震動往往具有兩大特征：一是斷層破裂的方向性效應明顯，主要表現為水平向地面運動兩個分量的幅值和頻譜相差很大;二是脈沖運動效應顯著，即地面運動速度甚至加速度和位移出現顯著的脈沖波形。Bertero等［1］對1971年San Fernando地震記錄進行研究后首次指出，靠近發震斷層觀測到的地震記錄中，含有不尋常的大幅值、長周期的脈沖運動，這種脈沖作用對結構響應影響很大，造成了地表結構的大量損壞。此后，在1979年美國Imperial Valley地震、1992年美國Landers地震、1994年美國Northridge地震和1995年日本Kobe地震中均觀測到這種脈沖型近斷層地震記錄，特別是在1999年土耳其Kocaeli地震、Duzce地震與我國臺灣地區的Chi－Chi地震中，采集到了大量的具有脈沖的近斷層地震記錄。它們的加速度峰值和速度峰值相當大，部分近斷層記錄可以歸結為近乎簡單沖擊型運動，產生顯著的沖擊力和變形，增大結構的破壞程度。如1999年我國臺灣地區Chi－Chi地震的Tcu084－W記錄分量，其加速度峰值達到1.157g，速度峰值為114.7 cm·s－1。這些脈沖型記錄的脈沖周期一般為1～3 s，更有部分記錄甚至超過6 s。

近年來，中外學者［2－4］對近斷層地震問題進行了深入而廣泛的研究，取得了一定的研究成果。但針對長周期脈沖型地面運動尤其是豎向地震動對結構抗震設計譜的影響研究還遠沒有達到實用的水平。由于地震的復雜性，特別是近斷層地震動脈沖問題的復雜性，需深入探討，并對照現行抗震設計規范開展對比研究。

本文中，基于前期建立的地震動數據庫，篩選具有長周期脈沖特性的地震地面運動記錄，分基巖和土場兩類場地，開展SDOF(Single Degree of Freedom)體系的彈性加速度譜計算，討論了不同周期的脈沖運動對豎向反應譜與水平向反應比值的影響特征及規律，基于現行中外抗震設計規范，建立了不同的脈沖周期及結構自振周期下反應譜比的取值關系式，并對決定豎向反應譜拐點周期的特征比值進行了參數回歸。研究成果可供結構抗震設計及規范修編參考。

1 近斷層長周期脈沖地震記錄庫及回歸模型

1.1 地震動記錄數據庫

為了進行近斷層地震動反應譜的研究，首要問題是建立長周期脈沖地震記錄樣本庫。隨著大型地震的不斷發生以及數字化精密地震儀的廣泛應用，大量靠近斷層的地震波被采集并記錄，本研究中采用的近斷層地震記錄均來自于美國伯克利大學強震記錄中心(http：//peer.berkeley.edu/)。選擇世界范圍內震級大于4.5級、斷層距不大于15 km的30余次主要地震中的132組記錄，已在前期工作中建立了統計分析的樣本庫［5］。斷層距的定義采用美國 UBC97規范［6］所給出的距離發震斷層的投影距離，如圖1所示。在前期已建立的近斷層地震記錄樣本庫的基礎上，進一步篩選具有長周期脈沖運動的記錄，場地劃分標準參見文獻［5］，所得記錄統計數據如表1所示，脈沖周期在0～3 s的記錄為56組，3～5 s的記錄為13組，5 s以上為24組，共計93組，構成本次研究的數據基礎。

圖1 斷層距R定義Fig.1 Definition of fault distance R

表1 近斷層長周期地震動記錄統計表Tab.1 Database of near-fault long-period earthquake records

1.2 地震動參數回歸模型

地震動參數回歸模型的確定，理論上應該分析比較各種模型擬合殘差的分布特征，并對不同模型回歸結果的顯著性進行檢驗，從而確定最優的回歸模型。實際上，做到這一點難度較大，并且對結構抗震的應用意義也不明顯，因此，結合工程地震研究中的常用模型，本研究中采用如下形式的近斷層地震動衰減模型：

式中：Y為需預測的地震動參數;M為矩震級(簡稱震級);R為斷層距，單位為 km;c1、c2、c3為回歸系數。

下文中，采用基于最小二乘法的一階段回歸方法對地震動特征值進行擬合回歸分析。

2 長周期脈沖效應下豎向與水平向反應譜比的分布特征及建議取值

目前，各國現行抗震設計規范仍然是基于反應譜理論構建的［6－9］。在抗震設計規范中，一般給出水平向加速度譜或動力放大系數譜，而豎向反應譜則通過基于水平向反應譜采用某一比例系數進行規定。美國UBC97規范［6］考慮近斷層效應明確了該比值的不確定性，在1 631.2條規定指出：通常條件下豎向設計譜通過在相應水平向設計譜的基礎上乘以2/3得到，但在有確定的二者比值的條件下，應選用確切的系數代替2/3;在考慮近斷層影響的情況下，應當用更確切的系數代替之。但UBC97規范并未說明更確切的系數的取值。

我國《建筑抗震設計規范》(GB50011 －2010)［7］規定豎向反應譜一般取為水平向譜值的0.65。我國《公路橋梁抗震設計細則》(JTG/T B 02 －01 －2008)［8］第5.2.5條規定：豎向設計加速度反應譜由水平向設計加速度反應譜乘以以下譜比值函數R：

式中：T為結構自振周期，單位為s。

美國核工業規范NRG(Nuclear Regulatory Guide1.60，R.G.1.60)［9］規定：

(1)結構自振頻率不超過2.5 Hz(T≥0.4 s)譜段的譜比值為2/3，在高于3.5 Hz(T＜0.286 s)頻段為1.0，中間段的比值與頻率有關(非直線)，見圖2所示。

(2)規范并沒有考慮這一比值與震級和斷層距離的關系，且明確此比值(設計譜)僅適用于非近斷層區。

圖2 美國NRG規范反應譜比Fig.2 Response spectra ratio of USA NRG code

在相關研究方面，中外學者［4，10～11］對于豎向與水平向反應譜比值的變化規律也進行了探討，如王國權［4］對Chi－Chi地震的近斷層地震動進行了研究，指出反應譜比隨斷層距增大而減小的規律并不一定成立，隨周期變化的規律也不明顯。耿淑偉等［11］在研究了豎向分量和水平向分量地震影響系數的比值后認為，在0.04～15 s的結構自振周期范圍內，這一比值與結構頻率有關，反應譜比的總平均值為0.65，接近規范取值的上限(2/3)，但不同頻段的譜比差別很大。雖然各研究者的途徑和方法不盡相同，但結果具有一定的相似性：在反應譜的高頻段，通常的取值(2/3)明顯低估了實際的反應譜比，在近斷層區和高頻段豎向與水平向反應譜比值將顯著大于遠場和長周期段。有必要對規范取值進行近斷層脈沖型地震下的修正。

基于以上研究經驗，采用前文所選取的93組記錄，對近斷層長周期脈沖效應下豎向與水平向加速度反應譜比值進行了統計分析，研究結果見圖3和圖4，圖3為不同脈沖周期的記錄組所對應的均值譜，圖4為考慮84%非超越概率的分布曲線。由圖3、圖4可得出以下結論：

(1)在反應譜的高頻段(T≤0.1 s)，脈沖周期為0～3 s記錄的反應譜比值顯著大于脈沖周期3～5 s和5 s以上記錄組，但都存在一峰點(約在0.1 s周期處)，脈沖周期3～5 s和5 s以上記錄組的峰點譜比值相差不大(約為1.1)，而脈沖周期為0～3 s記錄組的峰點譜比值為1.4。

(2)在反應譜的中頻段(0.1 s＜T≤1.0 s)，反應譜比值呈下降趨勢，并存在一谷點，脈沖周期3～5 s和5 s以上記錄組的谷點譜比值較為接近(約為0.45)，而脈沖周期為0～3 s記錄組的谷點譜比值約為0.35。

(3)在反應譜的低頻段(T＞1.0 s)，反應譜比值呈上升趨勢，并表現為脈沖周期越大則譜比值越大，特別是脈沖周期5 s以上記錄組的譜比最大值與高頻段的峰點譜比值基本相當。

(4)與現行《公路橋梁抗震設計細則》(JTG/T B 02 －01 －2008)［8］第5.2.5 條對比可發現，對于基巖場地，具有脈沖效應的近斷層地震動的豎向與水平向反應譜比顯著大于現行規范規定值(0.65);對于土層場地，對中短周期結構(T≤0.1 s)，脈沖型地震動的反應譜比大于規范值的10%～30%，對于長周期結構(T＞1 s)，脈沖型近斷層地震動反應譜比顯著大于規范給定值(0.5)。這說明具有脈沖效應的近斷層地震對豎向分量存在放大效應，豎向分量反應譜幅值顯著超出了現行規范的給定值。如若仍按照現行規范的規定進行脈沖型近斷層地震作用下結構的響應計算，勢必低估近斷層地震的影響，導致結構地震動響應偏小，加大近斷層區結構破壞的危險性。

將圖3(均值圖)和圖4(均值+σ)對比可發現，反應譜比在高頻段(T＜0.5 s)和低頻段(T＞6 s)顯著變化。84%非超越概率的統計結果表明，高頻段反應譜比在1.5～2.0之間分布，中頻段譜比值接近2/3，低頻段則超過2/3，高頻和低頻段差別明顯。

圖3 各脈沖周期段反應譜比的均值圖Fig.3 Mean curves of response spectra ratio for earthquake records with different pulse periods

圖4 均值+σ(84%的非超越概率)Fig.4 Response spectra ratio curves of mean value plus σ(84%exceeding probability)

圖5 分段線性化得到的反應譜比Fig.5 Response spectra ratio curve obtained by using piecewise linearization

表2 脈沖型近斷層地震動反應譜比建議取值Tab.2 Proposed response spectra ratios for near-fault pulse-type earthquake motions

由于脈沖周期3～5 s和5 s以上記錄組的統計結果在高頻段比較接近，按照平均值進行分段線性化，結果示于圖5，并在表2中列出了控制點周期處反應譜比的取值。

3 長周期脈沖效應下豎向加速度反應譜特征比值的參數回歸

按照Newmark等［12］的研究成果，典型加速度設計譜一般可以劃分為三個區段(如圖6所示)：即加速度敏感區(0～TC周期段)、速度敏感區(TC～TD周期段)和位移敏感區(T≥TD周期段)。TB、TC、TD分別為反應譜曲線的拐點周期值，在抗震設計規范中，TC通常稱為特征周期，相當于我國規范中的Tg。在各國抗震規范中，特征周期TC并不統一［6］，我國東部50年超越概率10%的Tc=0.3 s(均值)，西部Tc=0.4 s(均值);美國基巖場地加速度特征周期Tc處于0.25～0.35 s的占44%，0.3～0.45 s的占26%，0.45 s以上的占24%。TB一般采用TC之倍數確定，如 Eurocode8規范［13］采用0.25，而 UBC97 采用 0.2，Chopra等［14］建議拐點周期采用下式確定：

式中：PGA、PGV、PGD分別為地震動時程的峰值加速度、峰值速度和峰值位移。

圖6 典型加速度設計譜曲線Fig.6 Typical acceleration design response spectra curve

以所構建的長周期脈沖近斷層地震記錄庫為基礎(表1)，對豎向地震動的峰值加速度PGA、峰值速度PGV和峰值位移PGD三者的特征比值進行統計，結果匯總于表3。由表3和式3分析可得到以下結論：

(1)PGV/PGA值隨著震級的增大而顯著增大，說明TC值明顯增大，對于大震級，由于PGV/PGA達到0.192 6，TC值接近1.0 s，即加速度敏感區擴大，而TD值則超過6.0 s(PGD/PGV=0.775 2)，對應的結構自振周期范圍擴大，表明近斷層長周期脈沖地震動蘊含更多的能量，結構受到破壞的可能性顯著增大。

(2)對于基巖和土層兩種場地的對比來看，根據PGV/PGA和PGD/PGV，豎向反應譜曲線三個敏感區的寬度在中小震級地震下差別不大，但大震級下二種場地的反應譜特征周期比值差異明顯。因此，大震級的場地效應對脈沖運動具有放大作用，脈沖周期呈放大趨勢，這一特點與原始記錄脈沖周期的分布吻合。

(3)與我國現行《公路橋梁抗震設計細則》相比，大震級近斷層地震動下基巖和土層場地的特征周期TC分別達到0.82 s(0.165 7×5≌0.82)和1.1 s(0.2195×5≌1.1)，而現行規范中IV場地的特征周期僅為0.9 s，由此可見，近斷層脈沖效應增大了反應譜的加速度敏感區，即有更多的中短周期結構受到近斷層脈沖型地震的影響，基本涵蓋了目前常規的中短周期結構，有必要在常規結構抗震設計中考慮脈沖效應的影響。

由表3可知，對于近斷層區結構抗震設計，針對長周期脈沖地震動，由于各特征周期點顯著后延，有必要建立獨立的豎向抗震設計加速度反應譜的譜曲線，這一工作有待進一步深入開展。

在表3統計分析基礎上，由所建立的回歸分析模型(式1)對長周期脈沖地震動豎向分量的PGV/PGA和PGD/PGV特征比值進行回歸分析，得到了二者隨震級和斷層距變化的關系式，其系數取值見表4，并在圖7和圖8中給出了回歸模型與實際記錄的對比，表明吻合度良好。

由表4和圖7～圖8可見，隨著震級增大，各控制周期的特征比值呈增大趨勢。隨著斷層距增大，特征比值減小，但PGV/PGA和PGD/PGV的衰減程度不同，前者大于后者，說明反應譜加速度敏感區隨著斷層距增大而寬度變小，而速度敏感區則增大，這說明近斷層地震動的脈沖運動對中等周期結構的地震響應影響更為顯著。

表3 豎向地震動反應譜特征比值的統計結果Tab.3 Statistical results of response spectra characteristic ratio for vertical acceleration component

表4 豎向地震動特征值回歸結果Tab.4 Regression coefficients of response spectra characteristic ratio for vertical acceleration component

4 結論

以近斷層長周期脈沖型地震記錄為基礎，建立分析樣本庫和回歸模型，對脈沖效應對豎向反應譜與水平向反應譜比值及豎向分量加速度反應譜的特征比值進行了數據統計及回歸分析，可以得到以下結論：

(1)由于長周期脈沖運動的影響，豎向與水平向反應譜比值顯著高于現行抗震設計規范的取值，且并非一常量，隨著脈沖周期及結構自振周期而變化。

(2)豎向與水平向反應譜比值對中等及短周期結構的影響最大，脈沖周期在0～3 s的反應譜比最大，基本涵蓋了絕大多數實際結構。需對近斷層豎向地震作用引起足夠重視。

(3)確定豎向加速度譜各拐點周期的特征比值隨震級增大而增大，隨斷層距增大而減小，對反應譜曲線的加速度敏感區和速度敏感區影響較大。加速度敏感區的增大涵蓋了更多的常規中短周期結構，反應譜速度敏感區的增大將顯著地改變結構的響應特點，使結構產生較大的內力和位移，增強結構的延性要求。

(4)研究結果表明，對于近斷層區的結構抗震設計，由于長周期脈沖地震動導致豎向與水平向反應譜的比值明顯高于現行規范取值，且各特征周期點顯著后延，有必要建立獨立的豎向抗震設計加速度反應譜，用于此類區域結構的豎向抗震設計，這一工作有待深入開展。

［1］Bertero V V，et al.Establishment of design earthquakes：Evaluation of present methods［J］.Proceedings Interntl.Symp.On Earthq.Struct.Engrg.St.Louis，Univ.of Missouri－Rolla，1976，1：551 － 580.

［2］ Niazi M，Bozorgnia Y M.Behavior of near－source peak vertical and horizontal ground motion over SMART－1 array［M］.Taiwan，BSSA，1991，81：715－732.

［3］ Campbell K W. Empirical near－source attenuation relationships for horizontal and vertical components of peak ground acceleration， peak ground velocityand pseudoabsolution response spectra［J］.Seismological Research Letter，1997，68：154－179.

［4］王國權.臺灣集集地震近斷層地面運動特征［D］.北京：中國地震局地質研究所，2001，106－108.

［5］李新樂.近斷層區橋梁結構的設計地震與抗震性能研究［D］.北京：北京交通大學，2005：130－153.

［6］UBC97，International Council of Building Officials(ICBO).Uniform Building code［M］.Whittier CA，1997.

［7］GB50011－2010.建筑抗震設計規范［S］.北京：中國建筑工業出版社，2010：20－80.

［8］JTG/T B 02－01－2008.公路橋梁抗震設計細則［S］.北京：人民交通出版社，2008，17 －18.

［9］ Niazi M.Behavior of peak values and spectral ordinates of near－source strong ground motion over the SMART 1 array［M］.Tenera，L.P.1995 University Avenue Berkeley，CA 94704，1990.

［10］ Ambraseys N N，Douglas J.Near－field horizontal and vertical earthquake ground motions［J］. SoilDynamics and Earthquake Engineering，2003，23：1－18.

［11］耿淑偉，陶夏新，王國新.對抗震設計規范中地震作用規定的三點建議［C］.現代地震工程進展，2003：919－925.

［12］ Newmark N M，Hall W J.Seismic Design Criteria for Nuclear Reactor Facilities［C］.Proceedings of 4th World Conference on Earthquake Engineering， SandiagodeChile， 1969，II(B5)：1－12.

［13］ Eurocode 8，Design provision for earthquake resistance of structure［S］. Brussels： European Committee for Standardization，2003.

［14］ Chopra A K，Chintanapakdee C.Comparing response of SDOF systems to near－fault and far－fault earthquake motions in the context of spectral regions［J］.Earthquake Engineering and Structural Dynamics，2001，30：769－789.

Influence of near-fault long-period pulse-type earthquake ground motions on vertical response spectra

JIANG Hui1，LI Xin－le2，DOU Hui－juan2
(1.School of Civil Engineering，Beijing Jiaotong University，Beijing 100044，China;2.School of Civil Engineering and Architecture，Dalian Nationalities University，Dalian 116600，China)

The long－period large－amplitude pulse is a remarkable characteristic of near－fault earthquake ground motions，which is also an important factor leading the structure to damage or collapse.The database and regression relation for near－fault long－period pulse－type earthquake records were established based on the available data.The influence of pulse period and structural natural vibration period on the acceleration response spectra ratio of vertical to horizontal components was discussed in－depth.The response spectra ratios obtained from near－fault long－period pulse－type records are greater than those given by existing seismic design codes.The simplified relationship of response spectra ratio obtained by using piecewise linearization method varies as a function of pulse period and structural natural vibration period.The statistical results and attenuation relationship of characteristic ratio for vertical response spectra curves were given，which controls the inflection point periods of vertical response spectra and is related to earthquake magnitude and fault distance.

near fault earthquake;long－period pulse;vertical response spectra;response spectra ratio;characteristic ratio;regression

P315.3

A

國家自然科學基金資助項目(50908014);住房與城鄉建設部項目(2010－k3－56);自主科研基金資助項目(DC10040115，DC120101093);中央高校基本科研業務費專項資金資助項目(2011JBM073)

2011－04－08 修改稿收到日期：2012－04－01

江 輝 男，博士，講師，1977年生

李新樂 男，博士，副教授，1973年生