夏 峰

(中國地震局第一監(jiān)測中心, 天津 300180)

0 引言

土層地震反應(yīng)分析是場地條件對地震動影響研究中一個重要研究內(nèi)容,是場地工程地震安全性評價工作中關(guān)鍵技術(shù)手段之一,其合理性和可靠性對工程結(jié)構(gòu)抗震設(shè)計影響很大[1-2]。

目前,國際上考慮場地條件對設(shè)計地震動影響主要的還是一維土層地震反應(yīng)計算方法,其中最具有代表性的是SHAKE2000計算程序方法[3]和時域非線性DEEPSOIL5.0計算程序方法[4]。我國《工程場地地震安全性評價》規(guī)范中推薦使用的方法是由廖振鵬和李小軍提出的等效線性化法[5],該方法在我國地震安全性評價工作中被廣泛采用。以上國內(nèi)外三個土層計算程序在一定程度上滿足了以往工程設(shè)防要求,但隨著應(yīng)用的深入和不斷的工程經(jīng)驗積累,也發(fā)現(xiàn)存在不少問題。由于等效線性化方法是以線彈性分析為基礎(chǔ),是一種對土體非線性的近似處理方法,通過數(shù)值迭代方法來近似考慮土體非線性動力特征的,并沒有真正地反映土體的非線性,在計算中存在過度的算法共振效應(yīng)現(xiàn)象,尤其是厚層軟土場地在強震作用下尤為明顯,計算反應(yīng)譜明顯呈現(xiàn)峰值低,反應(yīng)譜矮、寬[6-7]的現(xiàn)象,不能很好地反映出場地條件真實的地震響應(yīng),為后繼的設(shè)防參數(shù)確定帶來較大的不確定性,進而留下工程隱患。

為此,袁曉銘等[8-9]采用頻率一致方法(FCM,Frequency-Consistent Method),提出了新一代土層地震反應(yīng)分析方法:SOILQUAKE軟件方法,在有限工況的初步檢驗表明該程序表現(xiàn)良好。另外,李兆焱等[10]用該軟件檢驗了巨厚場地計算結(jié)果,對比結(jié)果表明新一代土層地震反應(yīng)分析SOILQUAKE軟件方法結(jié)果與記錄結(jié)果吻合程度最高。本次工作以天津某厚層淤泥場地資料為基礎(chǔ),同時采用新一代土層地震反應(yīng)分析SOILQUAKE軟件方法和等效線性化方法進行計算,對比分析計算結(jié)果,為該地區(qū)類似厚層淤泥場地設(shè)防參數(shù)確定提供參考,進而探索新一代土層地震反應(yīng)分析SOILQUAKE軟件方法對特定場地條件影響、特定地震影響環(huán)境和與當前抗震規(guī)范銜接應(yīng)用問題。

1 場地地震構(gòu)造環(huán)境及場地條件概述

1.1 場地地震構(gòu)造環(huán)境

場地地震構(gòu)造環(huán)境是影響場地建、構(gòu)筑物設(shè)防的關(guān)鍵因素,它通常是客觀存在的,但又是隨著時間和認知水平變化的,地震工作者的任務(wù)是在給定風(fēng)險水準條件下預(yù)估未來一段時間內(nèi)(通常與建構(gòu)、筑物設(shè)計基準期有關(guān))設(shè)防水準。地震構(gòu)造環(huán)境在目前設(shè)防地震動確定過程中最直接的體現(xiàn)就是潛在震源區(qū)分布圖及相關(guān)參數(shù),它代表了當前對某地區(qū)地震及構(gòu)造工程應(yīng)用處理的認知水平。圖1是新版區(qū)劃圖[11]給出的對工程場區(qū)有影響范圍的潛源分布圖。

由圖1可知,場址區(qū)域范圍潛源多數(shù)呈北東和北西向分布,表明該區(qū)域活動構(gòu)造也多呈北東向和北西向展布。場址位于北西向構(gòu)造控制26號大沽7.0級潛源上。此外,19號天津7.0潛在震源區(qū)和18號唐山8.0級潛在區(qū)離場址較近,對場址地震環(huán)境影響明顯。

圖1 工程場址區(qū)域潛源分布圖Fig.1 Potential source distribution map of the project site area

1.2 場地條件概述

2 計算模型建立

等效線性化方法所需的場地模型土層鉆孔剖面參數(shù)和土動力學(xué)參數(shù)完全引用文獻[12],這里不再詳述。為了與等效線性化方法有可比性,SOILQUAKE軟件方法建模完全參照等效線性化方法分層建模,只是按照SOILQUAKE軟件方法輸入略作變換,見表1。其中土動力學(xué)參數(shù)也完全參照文獻[12]中表2的動三軸實驗測試數(shù)據(jù)。

表1 SOILQUAKE軟件方法場地模型土層鉆孔剖面參數(shù)

表2 兩種方法計算的場地地表峰值加速度結(jié)果Table 2 Calculation results of peak acceleration of ground surface by two methods

3 輸入地震動參數(shù)確定

考慮最新版區(qū)劃圖[11]研究成果,本次工作采用新版區(qū)劃圖地震統(tǒng)計區(qū)、潛源參數(shù)和地震動衰減關(guān)系,對場地進行了場地地震危險性概率分析計算?？紤]到新版區(qū)劃圖給出了極罕遇地震動參數(shù),故本次分別計算了場地100年超越概率1%和50年超越概率63%、10%、2%的基巖5%阻尼比的結(jié)果,分別對應(yīng)的場地設(shè)防水準是極罕遇地震動、罕遇地震動、基本地震動和多遇地震動水平,如圖2所示,其對應(yīng)的基巖加速度峰值分別為554.8 gal、369.2 gal、186.1 gal和48.2 gal。

圖2 場地基巖不同超越概率反應(yīng)譜Fig.2 Response spectra of site bedrock for different exceedance probabilities

基巖加速度的工程特性,主要由加速度峰值、頻譜和振動持續(xù)時間這三個要素決定。將地震危險性分析得到的基巖峰值加速度和反應(yīng)譜作為合成的目標函數(shù),結(jié)合適應(yīng)本場址地區(qū)地震活動特征的強度包絡(luò)函數(shù),采用擬合目標函數(shù)的三角級數(shù)迭加法合成基巖地震加速度時程,作為場地地震動反應(yīng)分析的輸入基巖地震動加速度的時程。給出三個以上相互獨立的隨機樣本時程,反應(yīng)譜擬合周期控制點數(shù)不得少于50個,周期控制點應(yīng)大體均勻地分布于周期的對數(shù)坐標上,控制點譜的相對誤差應(yīng)小于5%。本次工作在滿足上述技術(shù)要求前提下,以場地地震危險性分析得到50年超越概率2%、10%、63%和100年超越概率1%的基巖反應(yīng)譜為目標譜各合成三個相互獨立樣本時程,作為土層地震反應(yīng)分析的輸入(圖3)。

圖3 場地基巖不同超越概率人造地震動時程Fig.3 Artificial ground motion time histories of site bedrock for different exceedance probabilities

4 地震反應(yīng)分析

4.1 峰值加速度計算分析

在上述人造地震動時程輸入條件下分別采用等效線性化方法和SOILQUAKE軟件方法進行了土層反應(yīng)計算,分別得到了地表峰值加速度值(表2)和加速度反應(yīng)譜(圖4,黑實線為不同相位計算地表反應(yīng)譜,紅粗線為標定的設(shè)計譜)。

圖4 土層反應(yīng)計算地表反應(yīng)譜及標定設(shè)計譜Fig.4 Surface response spectrum and calibration design spectrum calculated by soil response

從表2計算結(jié)果來看,等效線性化方法(LSSRLI-1程序)計算的峰值加速度除了場地多遇地震(50 a-63%)條件下有所放大外,其他場地概率水準條件輸入下峰值加速度均有所減小,且輸入地震動強度越大峰值加速度減小幅度越大,這與以往等效線性化計算認識一樣;新一代土層地震反應(yīng)分析SOILQUAKE軟件方法計算的峰值加速度則較輸入值均有不同程度增幅,且增幅隨輸入地震動強度減小峰值加速度增加幅度越大,在場地極罕遇地震動(100 a-1%)輸入下仍是放大的,這與以往等效線性化方法計算認識相左,但這一趨勢與美國國家減輕地震災(zāi)害計劃[13]認識相一致。

4.2 設(shè)計反應(yīng)譜計算分析

為了與當前抗震設(shè)計規(guī)范相協(xié)調(diào),本次工作按建筑規(guī)范[14]給出的地震影響系數(shù)形式進行設(shè)計譜標定,其表達式為:

式中:α(T)為地震影響系數(shù),αmax為地震影響系數(shù)最大值;Tg為特征周期;γ為衰減指數(shù);η1為直線下降段的斜率調(diào)整系數(shù),取為0.02;η2為阻尼調(diào)整系數(shù),取為1.0。

為了與新版區(qū)劃圖[11]結(jié)果有一定可比性,本次計算的設(shè)計譜標定動力放大系數(shù)βmax均取為2.5,標定時先確定設(shè)計譜平臺段高度,然后反算出設(shè)計譜峰值加速度。表3列出了不同方法確定的設(shè)計地震動參數(shù)結(jié)果。為了對不同方法不同概率水準設(shè)計譜更直觀的比較,將表3的參數(shù)繪制成圖5所示。由于新版區(qū)劃圖結(jié)果僅給出設(shè)計譜峰值加速度,未給出設(shè)計譜特征周期值,圖5中新版區(qū)劃圖確定的極罕遇地震動設(shè)計譜特征周期值參照罕遇地震動設(shè)計

表3 不同方法確定的場地地表設(shè)計地震動參數(shù)Table 3 The design ground motion parameters of the ground surface determined by different methods

圖5 不同方法確定的場地地表設(shè)計地震動設(shè)計譜Fig.5 The design ground motion spectrum of the ground surface determined by different methods

譜特征周期值的結(jié)果給出。

從表3和圖5中不同方法確定的厚層淤泥場地不同概率水準的設(shè)計地震動參數(shù)來看,不同方法差別較大。與新版區(qū)劃圖結(jié)果相比,等效線性化方法在多遇地震作用下明顯偏低;在基本地震作用下等效線性化方法設(shè)計譜在特征周期大于0.9 s時與新版區(qū)劃圖結(jié)果基本吻合,但在小于0.9 s時與新版區(qū)劃圖結(jié)果差別較大,且可能會低估場地設(shè)防烈度;在罕遇地震作用下等效線性化方法設(shè)計譜在特征周期約大于0.9 s時大于新版區(qū)劃圖結(jié)果,但在小于0.9 s時略小于新版區(qū)劃圖結(jié)果;在極罕遇地震作用下等效線性化方法設(shè)計譜在特征周期約大于1.3 s時大于新版區(qū)劃圖結(jié)果,但在小于1.3 s時與新版區(qū)劃圖結(jié)果差別較大;相較新版區(qū)劃圖方法和等效線性化方法確定的設(shè)計譜,新一代土層地震反應(yīng)分析方法確定的設(shè)計譜明顯偏高,地震動輸入強度越弱越明顯,克服了文獻[6]提到的等效線性化方法在軟弱場地計算時出現(xiàn)的設(shè)計譜明顯矮、寬現(xiàn)象,尤其是在強地震動(極罕遇地震)輸入條件下與新版區(qū)劃圖結(jié)果較為接近,也與當下美國最新版NEHRP 2015(簡稱NEHRP)認識相一致。

5 結(jié)論

本文以天津濱海某厚層淤泥場地為例,采用新一代土層反應(yīng)方法和等效線性化方法進行了分析計算,并與新版區(qū)劃圖結(jié)果進行了比較,得到以下結(jié)論:

(1) 新一代土層反應(yīng)SOILQUAKE軟件方法在軟弱場地設(shè)計地震動參數(shù)確定時能體現(xiàn)一定的放大作用,尤其是強地震動作用下,克服了等效線性化方法在軟弱場地計算時出現(xiàn)的設(shè)計譜明顯矮、寬現(xiàn)象,與當前認識相一致,為軟弱場地重大工程設(shè)防參數(shù)確定提供了參考。

(2) 新一代土層反應(yīng)SOILQUAKE軟件方法在軟弱場地設(shè)計地震動參數(shù)確定較新版區(qū)劃圖結(jié)果設(shè)防標準有大幅度提高,考慮到相關(guān)抗震設(shè)防規(guī)范的協(xié)調(diào)性,還需進一步對軟弱場地進行大量強震記錄檢驗,以便在工程中更好地應(yīng)用。

(3) 本次分析結(jié)果僅以某厚層淤泥場地為例進行不同水準設(shè)防參數(shù)確定探討,對其他類型軟弱場地還需進一步分析不同強度地震作用及其相關(guān)因素對設(shè)防參數(shù)確定的影響。

致謝:本文計算所采用的新一代土層反應(yīng)方法SOILQUAKE程序由中國地震局工程力學(xué)研究所袁曉銘課題組免費提供,在程序使用過程中得到了李瑞山博士悉心指導(dǎo),在此一一致謝。