陳龍偉，馮浩，袁曉銘

（中國地震局地震工程與工程振動重點實驗室；中國地震局工程力學研究所，哈爾濱 150080）

場地地震反應分析是地震安全評價的關鍵環(huán)節(jié)，目前主要采用的等效線性化程序，如Shake2000及LSSRLI-1程序，進行場地地震反應數(shù)值計算，計算結果的可靠性難以得到實際地震數(shù)據(jù)的檢驗，且存在很大的不確定性。實際上，同一場地在不同地震中的反應亦不同，體現(xiàn)在地震動幅值、頻率成分的差異，這種差異主要源于震源機制、地震波傳播路徑以及場地土層動力特性的差異。準確評價場地地震反應的不確定性，可用于驗證現(xiàn)有場地地震反應分析結果的合理性。通過單臺地震數(shù)據(jù)探討場地反應不確定性及其影響因素，可以彌補數(shù)值結果可靠性的不足。但是，以往相關研究受到了單臺強震記錄缺乏這一客觀條件的限制，進而采用數(shù)值實驗手段。隨著全球地震臺網(wǎng)建設的加強，在地震活動性較為頻繁的地區(qū)可以獲取較為豐富的單臺地震數(shù)據(jù)。通過單臺地震數(shù)據(jù)研究場地反應是目前研究的熱門課題［1－4］。

地震中軟厚場地響應不同于其它類別場地，體現(xiàn)在對地震動長周期成分的放大作用。中國的建筑物抗震設計規(guī)范中明確規(guī)定對于軟弱土、可液化場地抗震設計要進行專門研究，抗震設計時需采取相應措施［5］。軟厚場地地基設計的關鍵問題之一就是地震動強度估計，現(xiàn)有場地反應分析程序在計算軟厚場地地震反應時基本失效，或者計算結果差異性很大。原因主要包括：場地基巖深度空間分布不均勻，即場地土層厚度分布空間不均，目前，雖然有較成熟面波測試法以及先進的地震測試設備，但是軟厚場地基巖面的確定仍然存在很大誤差［6］，這種誤差可能導致地震動強度估計誤差及規(guī)范中抗震設計譜中地震影響系數(shù)估計誤差。第二，對于軟厚場地尚無恰當特征參數(shù)能夠描述其地震反應，已有研究結果顯示場地特征參數(shù)都不能作為一般參數(shù)，用于描述軟厚場地地震反應的不確定性［7－8］。第三，中、強地震中軟厚場地土層表現(xiàn)出明顯的非線性特征，對于飽和砂土場地甚至會液化，所以，在選取場地模型描述土的粘彈性及耗能特性進行非線性分析時會產(chǎn)生很大的誤差，勢必導致分析結果的不確定性［9］。對軟厚場地地震反應不確定性分析最有效的方法就是采用實際地震數(shù)據(jù)，隨著單個地震臺記錄的累積，采用單臺地震記錄研究場地反應的不確定性。

通過單臺地震記錄研究場地反應的不確定性，首先需要選擇恰當?shù)膱龅貐?shù)。筆者采用場地放大函數(shù)（Site Amplification Function，SAF），場地放大函數(shù)綜合反應了地震場地效應，是衡量場地反應的綜合指標?？朔酝芯咳鄙賹崪y地震數(shù)據(jù)驗證的限制，特別是軟厚場地地震記錄，選取新西蘭Canterbury平原地區(qū)軟厚場地臺站記錄，以及作為參考臺站的B類場地LPCC和作為對比分析的C類場地HVSC。通過實測地震數(shù)據(jù)分析軟厚場地地震反應不確定性，同時剖析地震反應不確定性的影響因素（即地震動強度PGA、震級、震源距、震源方位角）及影響規(guī)律。

1 地震數(shù)據(jù)及場地地震反應

1.1 地震臺站及強震數(shù)據(jù)的選取

Christchurch市位于Canterbury平原地區(qū)，為多條源于南Alps山脈的河流沖積而成的沖積扇形平原。Christchurch市坐落在 Waimakiriri古河道上，Avon和Heathcote兩條河流穿城而過。附近地區(qū)地質構造復雜，地層主要由中間夾雜淤泥、黏土、泥炭、貝殼砂等的一系列卵石層組成，勘察結果顯示Christchurch地區(qū)的卵石層厚度多超過500 m，甚至超過2 km［10］。Christchurch市地區(qū)多為可液化的軟弱土層，2011年ChristchurchMw6.3地震中市區(qū)發(fā)生嚴重砂土液化，是自有液化調(diào)查歷史以來以液化震害為主要震害的一次地震，直接導致部分城區(qū)的廢棄［11］。

新西蘭擁有世界上較為完備的地震臺網(wǎng)（Get Net），擁有100多地震儀和180多個強震臺站，在地震中可以獲得高質量的地震數(shù)字記錄。在2010—2012年Canterbury地區(qū)一系列地震及余震中，Get Net獲得了豐富的地震數(shù)據(jù)，為地震研究提供了必要的寶貴資源。

選取Christchurch地區(qū)12個強震臺站（圖1），臺站主要位于Christchurch市區(qū)軟弱場地（主要是D、E類場地），搜集臺站2010—2012年系列強震記錄584條，地震震級Mw4.0～7.1，震源距100 km以內(nèi)。圖2顯示選取臺站記錄的Mw與震源距的分布圖。根據(jù)新西蘭抗震設計規(guī)范NZS1170.5（2004），Christchurch地區(qū)多為D類場地，然而 Wood等人根據(jù)現(xiàn)場實測剪切波速數(shù)據(jù)及場地土層厚度對Christchurch市區(qū)臺站場地類別進行修正，修正結果多為E類場地［12］。表1列出選取臺站場地信息，圖3給出臺站場地的剪切波速剖面。剪切波速剖面顯示，Christchurch市區(qū)覆蓋厚度超過20 m的軟弱土層，即軟厚場地。

表1 新西蘭臺站信息［12］Table 1 Information for strong motion stations in Christchurch［12］

圖1 選取的Christchurch地區(qū)強震臺站位置Fig.1 Locations of selected strong motion stations in Christchurch area

圖2 選取臺站地震記錄Mw－震源距分布圖Fig.2 Mw－Rhypo distribution of selected earthquake data

圖3 選取的強震臺站剪切波速剖面Fig.3 Shear wave velocity profiles of the selected stations

1.2 不確定性分析

選取反映場地地震反應指標，即場地地震放大函數(shù)（SAF），定義為場地地表加速度反應譜值與基巖輸入地震動反應譜的比值。由于新西蘭強震臺網(wǎng)沒有安裝井下強震儀，不能直接確定基巖輸入。為克服這一難題，選取參考臺站LPCC為出露基巖輸入。臺站LPCC位于Christchurch市東南部Lyttleton港口，臺站位于B類場地，場地為覆蓋于火山基巖上的淺薄土層（小于5 m）。選取參考臺站計算傳遞函數(shù)時，場地與參臺站地之間的距離相對于其震源距應很小，可以忽略傳播路徑的影響，否則進行傳播路徑幾何擴散的修正，即反應譜比值乘以震源距的比值近似進行修正［13］，計算式如式（1），采用對數(shù)值以便于進行分析。

式中：PSARhypo為目標臺站加速度譜值及臺站地震源距，PSAr、Rrhypo為參考臺站加速度譜值及臺站地震源距。

基于式（1），圖4給出選取臺站記錄計算得到的lg（SAF）曲線。由于選取的地震震級主要集中在4.0～4.8，且長周期信號噪音干擾影響較大，避免噪音影響，選取周期選取0～3.0 s，涵蓋了工程結構普遍的周期范圍。

由圖4可以得到：1）周期0～0.2 s，除CMHS、HVSC及PRPC臺站，其余臺站lg（SAF）均值小于0，即這些場地相對于參考LPCC對地震動具有減震效應；2）周期0.2～0.4 s，lg（SAF）均值隨著周期的增大而增大，且均大于0，即場地相對于參考LPCC對地震動具有放大效應；3）周期大于0.4 s，C類場地臺站HVSC，相對較硬，剪切波速vs30為422 m／s，lg（SAF）均值最小，且隨著周期的增大趨于減小，即硬土層對地震動中長周期的放大效應較弱；剪切波速最小的臺站REHS，vs30為141 m／s，lg（SAF）均值最大，說明軟弱土層對地震動長周期放大效應明顯，但隨著周期增大趨于降低；其余位于市區(qū)臺站的lg（SAF）均值基本相當，且隨著周期的增長保持穩(wěn)定，這些臺站中，RHSC場地剪切波速vs30為293 m／s，場地相對較硬，lg（SAF）均值在長周期段相對小于其他臺站的lg（SAF）均值。

圖4 臺站場地lg（SAF）均值Fig.4 Means lg（SAF）for selected stations

圖5給出所選臺站基于地震記錄的lg（SAF）標準差隨著周期的變化曲線?？梢钥闯觯琹g（SAF）標準差隨著周期的增長趨于降低。而較硬場地HVSC的lg（SAF）標準差最小，而對于軟弱場地D和E類場地，lg（SAF）標準差基本相當。圖6分別給出C、D和E類場地lg（SAF）均值和標準差均值。軟厚D和E類場地的場地放大函數(shù)較為接近與C類場地，HVSC有明顯區(qū)別；而E類場地離散性比D類場地小，比C類場地大。原因可能是D類場地空間分布為分散，而E類場地則集中在市區(qū)東部近海地區(qū)。要說明的是C、D和E類場地是指沒有進行修正的場地類別，即表1中的第4列。

圖5 臺站場地lg（SAF）標準差Fig.5 Standard deviations of lg（SAF）for selected stations

圖6 不同類別場地lg（SAF）均值（左）和標準差（右）Fig.6 Means（left）and standard deviations（right）of lg（SAF）for different site classes

2 不確定性影響因素分析

LPCC作為參照場地，那么僅剩HVSC為C類場地，記錄數(shù)量有限。由圖3的剪切波速剖面可見，位于Christchurch市區(qū)的D與E類場地均是厚度超過20 m軟厚場地，剪切波速140～300 m／s，多數(shù)200 m／s左右。通過上一節(jié)分析，D和E類場地對地震動響應較為接近（圖4和圖5），所以市區(qū)的D和E類場地記錄統(tǒng)一分析，即CBGS、CCCC、CHHC、CMHS、HPSC、NNBS、PRPC、REHS、RHSC和SHLC臺站記錄編為統(tǒng)一數(shù)據(jù)庫，計356條地震記錄，進行場地地震響應不確定性影響因素分析。

2.1 地震動強度（PGA）

首先，考慮地震動強度與場地SAF不確定性的關聯(lián)性。根據(jù)地表加速度PGA值，將記錄分為5組，即PGA≤0.2 m／s2、0.2 m／s2＜PGA≤0.4 m／s2、0.4 m／s2＜PGA≤0.6 m／s2、0.6 m／s2＜PGA≤1.0 m／s2和1.0 m／s2＜PGA，同時，保證每組擁有足夠多的地震記錄。圖7給出每個分組的lg（SAF）曲線，圖8給出每個分組lg（SAF）均值及標準差曲線。PGA水平不同，lg（SAF）均值基本一致；僅在周期大于1.0 s時，lg（SAF）均值差異明顯，而且PGA越大lg（SAF）均值越大，即說明周期大于1.0 s時地震動強度越大，場地對地震動放大效應越大。對于lg（SAF）標準差而言，隨著周期的增大，每個分組的lg（SAF）標準差均趨于降低。在PGA大于0.6 m／s2的2個分組，lg（SAF）標準差差別不大，但顯著大于其他分組；在PGA小于0.6 m／s2的3個分組，lg（SAF）標準差較為接近。

通過分析可知，場地地震動強度，雖僅選擇PGA作為地震動強度指標，對場地的放大函數(shù)lg（SAF）均值影響很小，僅在長周期階段有增大放大函數(shù)的趨勢；而地震動強度對場地放大函數(shù)不確定性影響明顯，當?shù)卣饎訌姸却笥谝欢〝?shù)值時顯著增大放大函數(shù)的不確定性。原因可能與場地土層強震下非線性有關，即場地土非線性提高了場地對地震動放大效應的不確定性。

圖7 不同地震動強度分組下場地放大函數(shù)lg（SAF）曲線Fig.7 lg（SAF）curves with respect to ground shaking intensity PGA bins

圖8 不同地震動強度分組下場地放大函數(shù)均值及標準差Fig.8 Means and standard deviations of lg（SAF）with respect to ground shaking intensity PGA bins

2.2 地震震級（Mw）

2010—2012年新西蘭一系列地震中，特別是2010年9月4日（Mw7.2）、2011年2月22日（Mw6.2），2012年6月14日（Mw6.0），新西蘭南島Canterbury地區(qū)遭遇了嚴重破壞。同時，此系列地震中也為地震研究提供寶貴的數(shù)據(jù)資源。所搜集的地震數(shù)據(jù)中Mw＞4.0的地震73次。分析場地地震動放大函數(shù)不確定性與地震震級的關聯(lián)性，同時，平衡每個分組中數(shù)據(jù)量，將地震數(shù)據(jù)根據(jù)震級Mw進行分組，即4.0≤Mw≤4.4、4.5≤Mw≤4.9和5.0≤Mw≤7.2。

圖9顯示不同地震震級分組下場地的放大函數(shù)lg（SAF）曲線。圖10給出不同地震震級分組放大函數(shù)lg（SAF）均值及標準差。由圖得到，地震震級越大，場地lg（SAF）均值趨向減小，但變化范圍不大；lg（SAF）的標準差基本相當，且與地震震級無關。

圖9 不同地震震級分組下場地放大函數(shù)lg（SAF）曲線Fig.9 lg（SAF）curves with respect to earthquake magnitude bins

圖10 不同地震震級分組下場地放大函數(shù)均值及標準差Fig.10 Means and standard deviations of lg（SAF）with respect to earthquake magnitude bins

2.3 地震震源距

將搜集的地震記錄按照震源距Rhypo進行分組，即5 km＜Rhypo≤15 km、15 km＜Rhypo≤25 km和25 km＜Rhypo。

圖11 不同地震震源距分組下場地放大函數(shù)lg（SAF）曲線Fig.11 lg（SAF）curves with respect to hypocentral distance bins

圖12 不同地震震源距分組下場地放大函數(shù)均值及標準差Fig.12 Means and standard deviations of lg（SAF）with respect to hypocentral distance bins

圖11顯示不同震源距分組下，場地放大函數(shù)lg（SAF）曲線。圖12給出lg（SAF）均值及標準差?？芍?，當頻率小于0.3 s時，lg（SAF）與距離無關；而當周期大于0.3 s時，震源距越大lg（SAF）越小，但是不同震源距分組lg（SAF）差別不顯著。對于lg（SAF）標準差而言，5 km＜Rhypo≤15 km分組的標準差明顯大于其它2個分組，而且其它2個分組lg（SAF）標準差基本重合。所以，距離震源較近區(qū)域場地放大函數(shù)的不確定性明顯大于距離震源較遠的區(qū)域，而且震源距超過一定數(shù)值的場地，其放大函數(shù)的不確定性與距離基本無關。

2.4 震源方位角

震源方位角是影響場地地震反應不確定性的另一因素，其體現(xiàn)了地震波傳播路徑以及不同震源區(qū)對場地放大函數(shù)不確定性的影響。采用的地震數(shù)據(jù)的震源與臺站的方位角主要分布在0～220°范圍內(nèi)，所以按照20°、30°、40°、60°、90°和120°震源方位角進行分組。同一方位角分組中，場地lg（SAF）標準差取平均值。這樣得到的不同震源方位角分組下，場地lg（SAF）標準差（圖13）以及給定周期下lg（SAF）標準差隨著方位角的變化規(guī)律（圖14）。由圖13可以看出，不同震源方位角分組場地lg（SAF）標準差基本接近，即場地放大函數(shù)不確定性與震源方位角選取基本無關。另一方面，給定周期條件下，圖14直觀顯示場地lg（SAF）標準差存在隨著方位角的增長的趨勢，但其增長幅值甚微，基本保持穩(wěn)定。

圖13 不同震源方位角分組場地放大函數(shù)lg（SAF）標準差Fig.13 Standard deviations of lg（SAF）with respect to azimuth bins

圖14 給定周期下震源方位角對lg（SAF）標準差的影響Fig.14 Standard deviations of lg（SAF）to azimuth under different periods

3 結 論

選取新西蘭Christchurch市區(qū)12個強震臺站2010—2012年間的臺站記錄584條，以淺硬場地LPCC作為參考輸入場地，對比分析選取場地的場地放大函數(shù)的不確定性。結果發(fā)現(xiàn)，較硬C類場地場地放大函數(shù)的不確定性最小。為探究場地放大函數(shù)不確定性源因，克服地震記錄數(shù)據(jù)量上的限制將Christchurch市區(qū)軟厚D類和E類場地臺站作統(tǒng)一分析，分析了地震動強度PGA、地震震級Mw及震源距Rhypo以及震源方位角等因素對場地不確定性的影響，分析結果發(fā)現(xiàn)：1）地震動強度水平對場地放大函數(shù)均值影響較小，僅在周期大于1.0 s段有差別，但場地放大函數(shù)標準差在PGA大于0.6 m／s2分組下顯著高于PGA小于0.6 m／s2分組；2）不同地震震級下，不同震級分組的場地放大函數(shù)標準差表現(xiàn)基本相當，即地震震級對場地放大函數(shù)不確定性影響?。?）地震震源距是影響場地放大函數(shù)不確定性的一個重要因素，距離震源越近的場地，其放大函數(shù)不確定性越高；4）不同震源方位角分組下場地lg（SAF）標準差差別很小，即場地放大函數(shù)不確定性與震源方位角的選取基本無關。

震源距（近場）、強地震動是場地放大函數(shù)不確定性的主要來源，而地震震級以及震源方位角的影響可忽略；另外，場地土層非線性趨于提高場地反應不確定性，其中機理有待深入探討。

軟厚土層場地反應分析的關鍵是場地放大函數(shù)模型的構建以及不確定性評估，采用合適的參數(shù)量化場地效應的不確定性是目前研究的重點，也是軟厚場地地震反應分析需要解決的問題。

［1］Kokusho T，Sato K.Surface-to-base amplification evaluated from Kik-net vertical array strong motion records ［J］. Soil Dynamics and Earthquake Engineering，2008，28：707-716.

［2］Rodriguez-Marek T，Montalva G A，Cotton F，et al.Analysis of single-station standard deviation using the Kik-net data［J］.Bulletin of the Seismological Society of America，2011，101（2）：1242-1258.

［3］Abrahamson N A，Hollenback J C.Application of single-station sigma ground motion prediction equations in practice ［C］／／Proceedings of the 15thWorld Conference on Earthquake Engineering，2012，Lisbon.

［4］Assimaki D，Li W，Steidl J H，et al.Site amplification and attenuation via downhole array seismogram inversion：a comparative study of the 2003 Miyagi-Oki aftershock sequence［J］.Bulletin of the Seismological Society of America，2008，98（1）：301-330.

［5］GB 50011—2010建筑物抗震設計規(guī)范 ［S］.北京：中國建筑工業(yè)出版社，2010.GB 50011—2010 Code for seismic design of building［S］.Beijing：China Architecture ＆ Building Press，2010.（in Chinese）

［6］Arai H，Tokimatsu K.S-wave velocity profiling by joint inversion of microtremor dispersion curve and horizontal-to vertical（H／V）spectrum ［J］.Bulletin of the Seismological Society of America，2005，95（5）：1766-1778.

［7］Boore D M，Thompson E M，Cadet H.Regional correlation of Vs30 and velocities averaged over depths less than and greater than 30 m ［J］.Bulletin of the Seismological Society of America，2011，101（6）：3046-3059.

［8］Choi Y，Stewart J P，Graves R W.Empirical model for basin effects accounts for basin depth and source location［J］.Bulletin of the Seismological Society of America，2005，95（4）：1412-1427.

［9］Bazzurro P，Cornell C A.Ground-motion amplification in nonlinear soil sites with uncertain properties ［J］.Bulletin of the Seismological Society of America，2011，94（6）：2090-2109.

［10］Tonkin ＆ Taylor Ltd.Christchurch central city geological interpretative report（Version 1.1）［R］.Report for Christchurch City Council，Christchurch，New Zealand，2012.

［11］陳龍偉，袁曉銘，孫銳.2011年新西蘭Mw6.3地震液化及巖土震害述評［J］.世界地震工程，2013，29（3）：1-9.

Chen L W，Yuan X M，Sun R.Review of liquefaction phenomena and geotechnical damage in the 2011 New Zealand Mw6.3 earthquake ［J］.World Earthquake Engineering，2013，29（3）：1-9.（in Chinese）

［12］Wood C M，Cox B R，Wotherspoon L M，et al.Dynamic site characterization of Christchurch strong motion stations ［J］.Bulletin of the New Zealand Society for Earthquake Engineering，2011，44（4）：195-204.

［13］Steidl J H，Tumarkin A G，Archuleta R J.What is a reference site？ ［J］.Bulletin of the Seismological Society of America，1996，86（6）：1733-1748.

（編輯胡英奎）