蔣建平，閆澤豐

（上海海事大學(xué)海洋環(huán)境與工程學(xué)院，上海 201306）

地震會(huì)對(duì)地表的各類(lèi)工程造成嚴(yán)重危害[1～3]。相關(guān)調(diào)查結(jié)果表明，地震所造成的建（構(gòu)）筑物破壞，除了與地震烈度和結(jié)構(gòu)物特性有關(guān)外，還與建筑物所在場(chǎng)地的地質(zhì)條件有密切關(guān)系[4～10]。

早在1908年Wood在分析1906年舊金山大地震時(shí)，就已認(rèn)識(shí)到場(chǎng)地條件對(duì)地震震害的重要影響。在這次地震中，市內(nèi)軟弱地基上的震害更為顯著；1967年委內(nèi)瑞拉地震時(shí)，場(chǎng)地覆蓋層厚度不同，震害存在明顯差異；1968年和1970年菲律賓馬尼拉地震中，不同高度建筑物破壞程度隨沉積土層厚度變化而明顯不同。

地震波由震源向地表面?zhèn)鞑?，在?jīng)過(guò)不同地層時(shí)，因各地層對(duì)地震波的透射、反射、阻尼等性質(zhì)不同，會(huì)造成地表各點(diǎn)的地震反應(yīng)不一致。有時(shí)由于地層變化大，相鄰場(chǎng)地的烈度可以相差1～2度，或在區(qū)域內(nèi)出現(xiàn)小范圍的震害異常區(qū)。

地層結(jié)構(gòu)是場(chǎng)地地質(zhì)條件的重要內(nèi)容之一。地層結(jié)構(gòu)是地層的疊覆、堆積或空間組合狀態(tài)，地層的結(jié)構(gòu)效應(yīng)是工程地層學(xué)的重要研究?jī)?nèi)容[11]。目前有關(guān)地層結(jié)構(gòu)對(duì)地震響應(yīng)影響方面的研究很少[12～14]，相關(guān)技術(shù)規(guī)程也只是將不同地層按剪切波速作等效化處理[15]。

地層地震反應(yīng)分析是當(dāng)前工程場(chǎng)地地震安全性評(píng)價(jià)工作中的重要組成部分，直接為上部結(jié)構(gòu)提供地震動(dòng)參數(shù)，為建筑物的抗震設(shè)防提供必不可少的技術(shù)支持；地層地震反應(yīng)分析也是地震工程中的重要研究領(lǐng)域，是結(jié)構(gòu)抗震設(shè)計(jì)的基礎(chǔ)。因此，本研究具有一定的理論價(jià)值和實(shí)際意義。

1 地層結(jié)構(gòu)與地震波設(shè)定

本文采用三層地層結(jié)構(gòu)并按5種結(jié)構(gòu)類(lèi)型模擬計(jì)算地震響應(yīng)。選擇軟土地區(qū)常見(jiàn)的淤泥質(zhì)粘土、粉質(zhì)粘土、粉質(zhì)粘土夾粉砂三種土性地層，分別代表土層強(qiáng)度的軟、中、硬三類(lèi)，并概化組合成不同地層結(jié)構(gòu)形式。同時(shí)對(duì)三種地層的剪切波速進(jìn)行等效化處理，作為均質(zhì)土層，用于比較分析。為利于在理論上分析地層結(jié)構(gòu)的影響，各土層設(shè)為同等厚度(2m)。地層基本特性及結(jié)構(gòu)類(lèi)型組合見(jiàn)表1、表2。

表1 地層基本特性Table 1 The basic characteristics of strata

表2 地層結(jié)構(gòu)組合類(lèi)型Table 2 The compositional type of strata structure

地震動(dòng)是在地層下部的基巖面輸入水平地震波，選用頻譜差異較大的El-Centro地震波和Concrete地震波，其加速度時(shí)程曲線(xiàn)如圖1所示。從圖2的F譜曲線(xiàn)上可發(fā)現(xiàn)兩地震波的頻譜特性有較大差異，El-Centro地震波的頻率范圍較寬，而Concrete地震波的頻率范圍較窄，主要集中在4Hz以?xún)?nèi)。本文計(jì)算中抗震設(shè)防烈度采用8度，因此兩地震波的加速度峰值都調(diào)整到0.30g。

具體計(jì)算時(shí)是采用ABAQUS軟件中的動(dòng)力模塊，土體采用摩爾庫(kù)倫模型，動(dòng)力邊界采用吸收邊界。

圖1 地震波加速度時(shí)程曲線(xiàn)Fig.1 The acceleration time history curves of seismic waves

圖2 地震波Fourier譜曲線(xiàn)Fig.2 The Fourier spectrum of seismic waves

2 模擬結(jié)果及其分析

通過(guò)計(jì)算得到地面點(diǎn)有關(guān)位移、加速度、速度的時(shí)程曲線(xiàn)、峰值曲線(xiàn)及Fourier譜、能量譜、反應(yīng)譜曲線(xiàn)如圖3～圖13所示。

圖3 地面點(diǎn)水平位移時(shí)程曲線(xiàn)Fig.3 The displacement-time curves of point on the ground

圖4 地面點(diǎn)水平位移峰值Fig.4 The displacement peak value curves of point on the ground

圖5 地面點(diǎn)水平位移殘余值比較Fig.5 The displacement residual value curves of point on the ground

圖6 地面點(diǎn)加速度時(shí)程曲線(xiàn)Fig.6 The acceleration time history curves of point A on the ground

圖7 地面點(diǎn)加速度峰值Fig.7 The acceleration peak value curves of point on the ground

圖8 地面點(diǎn)速度時(shí)程曲線(xiàn)Fig.8 The velocity time history curves of point on the ground

圖9 地面點(diǎn)速度峰值比較Fig.9 The velocity peak value curves of point on the ground

圖10 地面點(diǎn)Fourier譜曲線(xiàn)Fig.10 The Fourier spectrum of point on the ground

圖11 地面點(diǎn)能量譜曲線(xiàn)比較Fig.11 The Power spectrum of point on the ground

圖12 地面點(diǎn)反應(yīng)譜(頻率)曲線(xiàn)比較Fig.12 The Response spectrum (frequency)of point on the ground

圖13 地表A點(diǎn)的反應(yīng)譜(周期)曲線(xiàn)比較Fig.13 The Response spectrum (Period)of point on the ground

由圖3的地面點(diǎn)水平位移時(shí)程曲線(xiàn)比較圖中可發(fā)現(xiàn)，在地震波作用下，雖然同一地震動(dòng)作用下波形類(lèi)似，但都顯示出不同的地層結(jié)構(gòu)可導(dǎo)致地面點(diǎn)位移時(shí)程的差異。在El-Centro地震波作用下，ΙΙ、ΙΙΙ地層結(jié)構(gòu)類(lèi)型的曲線(xiàn)位于等效均質(zhì)土層結(jié)構(gòu)V曲線(xiàn)上方，Ι、ΙV地層結(jié)構(gòu)類(lèi)型曲線(xiàn)則位于V曲線(xiàn)下方。說(shuō)明淤泥質(zhì)粘土的C1地層與粉質(zhì)粘土夾粉砂的C3地層位于粉質(zhì)粘土的C2地層上時(shí)，其地面點(diǎn)的位移大于地基為均質(zhì)層C0時(shí)的位移；而當(dāng)C1、C2地層位于C3地層上時(shí)，地面點(diǎn)位移小于C0時(shí)的位移。在Concrete地震波作用下，代表均質(zhì)層的V曲線(xiàn)基本上與ΙΙΙ曲線(xiàn)在大部分時(shí)段幾乎重合，ΙΙ曲線(xiàn)位于V曲線(xiàn)上方，Ι、ΙV曲線(xiàn)位于V曲線(xiàn)下方。說(shuō)明從上往下C3、C1、C2排列的地層結(jié)構(gòu)情形下地面點(diǎn)的位移大于地基為均質(zhì)層C0時(shí)的位移；而當(dāng)?shù)貙覥1、C2位于地層C3上時(shí)，則小于C0狀態(tài)下位移。

圖4顯示，各地層結(jié)構(gòu)情形下地面點(diǎn)的水平位移峰值也各不相同，特別是El-Centro地震波作用條件下。具體來(lái)說(shuō)，El-Centro地震波作用下，C1、C2、C3排列的地層結(jié)構(gòu)的水平位移峰值最大；C2、C1、C3排列，C1、C3、C2排列，C3、C1、C2排列地層結(jié)構(gòu)的位移峰值都小于均質(zhì)層C0的位移峰值。在Concrete地震波作用下，C1、C3、C2地層結(jié)構(gòu)的水平位移峰值最大；C1、C2、C3地層結(jié)構(gòu)與C0均質(zhì)層的水平位移峰值相等；C3、C1、C2結(jié)構(gòu)類(lèi)型則與C2、C1、C3排列的水平位移峰值相等。

從圖5上可看出，兩種地震波作用下地面點(diǎn)的水平位移殘余值的形態(tài)類(lèi)似，都顯示C1C2C3、C2C1C3地層結(jié)構(gòu)的位移殘余值大，C3C1C2、C1C3C2地層結(jié)構(gòu)的位移殘余值小，而C0均質(zhì)層的位移殘余值處于中間層次。

從圖6可發(fā)現(xiàn)，地面點(diǎn)的各加速度時(shí)程曲線(xiàn)與原輸入波的波形類(lèi)似；對(duì)兩種地震波，C0均質(zhì)層的波形幅度處在中間層次，C1C2C3與C1C3C2地層結(jié)構(gòu)的波形幅度處在較大層次，而C3C1C2與C2C1C3的地層結(jié)構(gòu)處在較小層次。圖7的加速度峰值曲線(xiàn)也顯示出這一特點(diǎn)。圖8、圖9所示的地面點(diǎn)速度時(shí)程曲線(xiàn)、速度峰值，也顯示了與加速度類(lèi)似的結(jié)果。

圖10反映El-Centro及Concrete地震波作用下地面點(diǎn)的Fourier譜曲線(xiàn)簇，兩種地震波的作用效應(yīng)有相似規(guī)律。圖示表明，El-Centro地震波作用下，C0均質(zhì)層的波形高度處在中間層次，C1C2C3、C1C3C2排列的地層結(jié)構(gòu)的波形幅度處在較高層次，而C3C1C2、C2C1C3排列的地層結(jié)構(gòu)的波形大小處在較低層次；Concrete地震波作用下，C1C2C3、C1C3C2排列的地層結(jié)構(gòu)的主頻率集中在2.5Hz左右，C3C1C2、C2C1C3排列地層結(jié)構(gòu)的主頻率集中在2.0Hz左右，而C0均質(zhì)層的主頻率集中在2.25Hz左右。

同樣，從圖11的地面點(diǎn)能量譜曲線(xiàn)簇，也可得出與之類(lèi)似的結(jié)論。

圖12為兩地震波作用下以頻率為橫坐標(biāo)的地面點(diǎn)反應(yīng)譜曲線(xiàn)簇，且其中還加進(jìn)了輸入波的反應(yīng)譜曲線(xiàn)。從圖上可看出：C0均質(zhì)層的譜高度處在中間層次，C1C2C3、C1C3C2排列的地層結(jié)構(gòu)的譜高度處在較高層次，而C3C1C2、C2C1C3地層結(jié)構(gòu)的譜高度處在較低層次；各情形的主頻值與圖10、圖11一致；輸入波的譜高處于較低層次，其主頻值也偏小。

圖13為兩地震波作用下以周期為橫坐標(biāo)的地面點(diǎn)反應(yīng)譜曲線(xiàn)簇，同樣可得出與圖12類(lèi)似的結(jié)論，只是主頻率較大的其主周期變成了較小。

3 結(jié)論

本文利用軟土地區(qū)常見(jiàn)的三種巖性土層，通過(guò)不同的地層結(jié)構(gòu)組合，模擬分析了地震波作用下地面點(diǎn)的響應(yīng)，得出如下認(rèn)識(shí)：

（1）其它條件相同而地層結(jié)構(gòu)不同，則地表的地震響應(yīng)也不同。所以，工程場(chǎng)地地震反應(yīng)分析，應(yīng)考慮作為地質(zhì)條件重要內(nèi)容之一的地層結(jié)構(gòu)的影響，不可忽視。

（2）在基巖面輸入頻譜特性有較大差異的El-Centro及Concrete地震波時(shí)，地表的地震響應(yīng)雖不相同但規(guī)律性類(lèi)似，說(shuō)明地震作用下地層結(jié)構(gòu)的影響普遍存在。

（3）地震波作用下，地表的加速度時(shí)程與峰值、速度時(shí)程與峰值、位移時(shí)程與峰值、位移殘余值，以及Fourier譜、能量譜、反應(yīng)譜，基本上都存在均質(zhì)地層的值處在中間層次，軟-中-硬地層結(jié)構(gòu)、軟-硬-中地層結(jié)構(gòu)的值較大，而硬-軟-中地層結(jié)構(gòu)、中-軟-硬地層結(jié)構(gòu)的值較小。說(shuō)明從上往下，土層彈性模量依次為高-低-中及中-低-高的地層組合結(jié)構(gòu)的抗震效果較好。

References)

[1]林宇亮,楊果林,趙煉恒. 地震條件下?lián)鯄笳承酝林鲃?dòng)土壓力研究[J]. 巖土力學(xué),2011,32(8):2479-2487.

Lin Y L, Yang G L, Zhao L H. Effect of foundation stiffness and anisotropy on the seismic response of concrete dams[J]. Rock and Soil Mechanics, 2011,32(8):2479-2487.

[2]Kalkan E, Luco N. Special issue on earthquake ground-motion selection and modification for nonlinear dynamic analysis of structures[J]. Journal of Structural Engineering, 2011,137(3):277-286. doi:10.1061/(ASCE)ST.1943-541X.0000355.

[3]Huang Y N, Whittaker A S, Luco N, et al. Scaling earthquake ground motions for performance-based assessment of buildings[J].Journal of Structural Engineering, 2011,137(3):311-320.doi:10.1061/(ASCE)ST.1943-541X.0000155.

[4]劉昌森,呂美麗. 上海地區(qū)地震放大效應(yīng)的初步探討[J]. 上海地質(zhì),1998,19(1):7-13.

Liu C S, Lu M L. Preliminary analysis of the seismic amplif i cation effect in Shanghai[J]. Shanghai Geology, 1998,19(1):7-13.

[5]高廣運(yùn),陳青生,何俊鋒,等. 地下水位上升對(duì)上海軟土場(chǎng)地地震反應(yīng)的影響[J]. 巖土工程學(xué)報(bào), 2011,33(7):989-996.

Gao G Y, Chen Q S, He J F, et al. Effect of rise of groundwater table on seismic ground response of soft soil in Shanghai[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(7):989-996.

[6]蔣建清,楊果林,李麗民. 基于SΙMULΙNK 的高填方加筋路堤地震反應(yīng)分析[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,42(1):227-233.

Jiang J Q, Yang G L, Li L M. Seismic response of high-filled reinforced embankment based on SΙMULΙNK[J]. Journal of Central South University (Science and Technology), 2011,42(1): 227-233.

[7]陳可洋. 多尺度隨機(jī)孔洞介質(zhì)地震波正演及其逆時(shí)成像[J]. 上海國(guó)土資源,2011,32(1):84-88.

Chen K Y. Seismic wave forward simulation and reverse-time migration with multi-scale random carve medium[J]. Shanghai Land & Resources, 2011,32(1):84-88.

[8]Takahashi H, Kitazume M, Noguchi T, et al. Seismic-resistant effect of cement-treated soil for quay wall[J]. Ground Improvement,2011,164(13):151-160. doi:10.1680/grim.2011.164.3.151.

[9]Hwang G S, Chen C H. Analysis of cases of gravity quaywall movement during earthquakes[J]. Geotechnique,2011, 61(3):199-210. doi:10.1680/geot.3625.

[10]杜建國(guó),林皋. 地基剛度和不均勻性對(duì)混凝土大壩地震響應(yīng)的影響[J]. 巖土工程學(xué)報(bào),2005,27(7):819-824.

Du J G, Lin G. Effect of foundation stiffness and anisotropy on the seismic response of concrete dams[J]. Chinese Journal of Geotechnical Engineering, 2005,27(7):819-824.

[11]蔣建平,劉文白. 巖土工程地層結(jié)構(gòu)效應(yīng)與層組參數(shù)[M]. 北京:人民交通出版社,2011.

Jiang J P, Liu W B. The strata structure effect in geotechnical engineering[M]. Beijing: China Communications Press, 2011.

[12]張嘎,張建民. 地震作用下成層地基中結(jié)構(gòu)物的動(dòng)力響應(yīng)[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2004,44(12):1642-1645.

Zhang G, Zhang J M. Effect of foundation stiffness and anisotropy on the seismic response of concrete dams[J]. Journal of Tsinghua University (Science and Technology), 2004,44(12):1642-1645.

[13]陳之毅,沈昊. 軟土地下工程抗震數(shù)值模擬的若干關(guān)鍵問(wèn)題[J].上海國(guó)土資源,2011,32(4):83-87.

Chen Z Y, Shen H. Key issues for seismic analysis of underground structures in soft soils[J]. Shanghai Land & Resources,2011,32(4):83-87.

[14]Okamura M, Tamamura S. Seismic stability of embankment on soft soil deposit[J]. International Journal of Physical Modelling in Geotechnics, 2011,11(2):50-57. doi:10.1680/ijpmg.2011.11.2.50.

[15]GB50011-2010, 中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn):建筑抗震設(shè)計(jì)規(guī)范[S]. 北京:中國(guó)建筑工業(yè)出版社,2010.

National Standard of PRC: Code for seismic design of buildings(GB50011-2010)[S]. Beijing: China Architecture &Building Press, 2011.