魏成前，于彥彥，丁海平

（蘇州科技大學(xué)江蘇省結(jié)構(gòu)工程重點(diǎn)實(shí)驗(yàn)室，江蘇蘇州 215011）

引言

國(guó)內(nèi)外大量的地震觀測(cè)記錄、數(shù)值模擬及震害調(diào)查研究結(jié)果表明，沉積盆地作為一種特殊的場(chǎng)地，會(huì)顯著放大其內(nèi)部長(zhǎng)周期地震動(dòng)，并延長(zhǎng)地震動(dòng)的持續(xù)時(shí)間，這對(duì)于盆地內(nèi)的高層建筑、大跨度橋梁和其它自振周期較長(zhǎng)的工程結(jié)構(gòu)非常不利。1985年墨西哥8.1級(jí)地震中，遠(yuǎn)離震中400 km的墨西哥城中高層建筑破壞嚴(yán)重，震后研究認(rèn)為這與該城市坐落于深厚沉積盆地有關(guān)［1-2］。其后在1995年日本神戶地震中的大阪盆地［3］、1999年中國(guó)臺(tái)灣集集地震中臺(tái)北盆地［4］、2008年四川汶川地震中的渭河盆地［5］及2014年魯?shù)榈卣瘕堫^山鎮(zhèn)盆地［6］等地出現(xiàn)的震害，均表明盆地對(duì)震害會(huì)有加重作用。

現(xiàn)有研究表明，入射地震波的角度、幅值和頻譜特性等參數(shù)對(duì)盆地地震效應(yīng)具有顯著的影響。張建經(jīng)等［7］研究了SH波斜入射時(shí)，地震波的入射角、卓越頻率和盆地傾角對(duì)盆地地表地震響應(yīng)的影響。結(jié)果表明，隨著地震波入射角的增大，加速度反應(yīng)譜的峰值有減小的趨勢(shì)；地震波的卓越頻率越大，影響盆地地表加速度峰值的地震波入射角的變化范圍就越小。Khanbabazadehand等［8］通過輸入不同幅值的強(qiáng)震地震波研究了盆地介質(zhì)分層和土層參數(shù)的變化對(duì)盆地動(dòng)力響應(yīng)的影響，發(fā)現(xiàn)雙層介質(zhì)盆地的共振周期普遍大于單層盆地，雙層軟黏土盆地的譜放大倍數(shù)最大。Zhu等［9］采用顯式有限差分法詳細(xì)研究了SH波垂直下入射波主頻和盆地結(jié)構(gòu)幾何形狀之間的相互作用及盆地內(nèi)外介質(zhì)阻抗比對(duì)盆地地震響應(yīng)的影響，結(jié)果表明盆地地面運(yùn)動(dòng)的振幅和分布特點(diǎn)是由入射波與盆地動(dòng)力特性以及盆地幾何結(jié)構(gòu)相互作用決定的。陳學(xué)良等［10］研究了Rayleigh面波作用下盆地介質(zhì)參數(shù)、盆地尺寸對(duì)盆地場(chǎng)地地震動(dòng)特性的影響規(guī)律，發(fā)現(xiàn)當(dāng)盆地內(nèi)介質(zhì)由軟變硬時(shí)，盆地對(duì)地表地震動(dòng)的放大作用由大變小。劉中憲等［11］基于高精度譜元法，研究了不同斷層傾角及不同震源頻率下盆地的地震響應(yīng)規(guī)律，發(fā)現(xiàn)隨著斷層傾角的增大，地表加速度的峰值有減小的趨勢(shì)。震源的主頻較高時(shí)，盆地的地震動(dòng)放大效應(yīng)顯著增強(qiáng)。梁佳利等［12］基于粘彈性邊界和等效線性化方法研究了二維均勻盆地的非線性地震反應(yīng)特征，結(jié)果顯示入射波的頻譜和幅值對(duì)盆地地震響應(yīng)影響非常明顯。但以往對(duì)于成層盆地模型在SV波斜入射下地震動(dòng)的時(shí)、頻域綜合放大特性的研究相對(duì)較少。此外，當(dāng)前對(duì)于盆地場(chǎng)點(diǎn)的工程地震效應(yīng)評(píng)價(jià)中，采用較多的方法是根據(jù)場(chǎng)點(diǎn)處的鉆孔資料建立一維土層模型，進(jìn)而分析地震波垂直入射下地表點(diǎn)的地震動(dòng)響應(yīng)，而忽略了入射角度的影響以及二維盆地模型橫向上的波速不均勻分布所引起的盆地放大作用。

文中以SV波為入射地震波，結(jié)合顯式有限元法和透射邊界技術(shù)，詳細(xì)研究了梯形形狀的二維成層盆地模型，相比盆地各場(chǎng)點(diǎn)處厚度不等的一維水平成層模型在相同輸入波、不同入射角度下地表地震動(dòng)的放大特征，并以地震動(dòng)峰值大小、地震動(dòng)峰值放大系數(shù)和傅里葉譜譜比作為定量分析的指標(biāo)。

1 計(jì)算方法及計(jì)算模型

1.1 計(jì)算方法

模擬地震波斜入射條件下盆地場(chǎng)地地震響應(yīng)時(shí)，首先需要計(jì)算自由波場(chǎng)，以此作為斜入射條件下盆地模型的波動(dòng)輸入。均勻或成層半空間自由場(chǎng)可以根據(jù)波動(dòng)傳播規(guī)律在頻域內(nèi)計(jì)算得到［13］。當(dāng)使用有限元方法分析地震波散射問題時(shí)，有必要從半無限空間介質(zhì)中切出有限的計(jì)算區(qū)域，同時(shí)，必須在邊界處引入適當(dāng)?shù)娜斯み吔鐥l件模擬區(qū)域內(nèi)外行波的傳播過程，以達(dá)到消除截?cái)噙吔缣幦斯し瓷洳ǖ哪康?。文中采用透射邊界?4］，在模型的左、右兩側(cè)邊界和底邊界設(shè)置二階透射邊界。計(jì)算區(qū)域內(nèi)的節(jié)點(diǎn)可以分為內(nèi)節(jié)點(diǎn)和邊界節(jié)點(diǎn)，位于人工邊界上的節(jié)點(diǎn)u0稱為人工邊界節(jié)點(diǎn)，其余的節(jié)點(diǎn)均為內(nèi)節(jié)點(diǎn)u，全部?jī)?nèi)節(jié)點(diǎn)的運(yùn)動(dòng)方程為

式中：M為質(zhì)量矩陣；C為阻尼矩陣；K為剛度矩；P為外力矢量，文中P為0。

人工邊界節(jié)點(diǎn)運(yùn)動(dòng)采用多次透射公式計(jì)算，即

式中：N為透射階數(shù)（文中取2階）；0表示人工邊界點(diǎn)；j為與節(jié)點(diǎn)0相鄰的計(jì)算點(diǎn)；u0p+1為邊界節(jié)點(diǎn)p+1時(shí)刻的位移；u p+1-jj為p+1時(shí)刻j計(jì)算點(diǎn)的位移。

關(guān)于散射波場(chǎng)的輸入問題，其輸入是來自半無限域的彈性波。在人工邊界附近區(qū)域，既有穿過人工邊界向外傳播的外行波，又有向計(jì)算區(qū)內(nèi)部傳播的內(nèi)行波。由于透射邊界只能模擬外行波，所以必須從總波場(chǎng)中分離出外行波。假設(shè)分離出的外行波位移為us，u為總波場(chǎng)位移，ur為參考波場(chǎng)位移（一般取自由場(chǎng)位移），三者存在以下關(guān)系，

將式（4）代入式（2）得到邊界節(jié)點(diǎn)p+1時(shí)刻的總場(chǎng)位移為

利用式（5）即可實(shí)現(xiàn)散射問題的人工邊界節(jié)點(diǎn)運(yùn)動(dòng)方程的求解，此方程也同時(shí)實(shí)現(xiàn)了地震波的輸入。

1.2 計(jì)算模型

文中所采用的盆地計(jì)算模型如圖1（a）所示，模型尺寸為5 000 m×1 000 m，盆地頂部寬4 000 m，底部寬3 600 m，深度H為200 m，盆地傾角α=45°。盆地內(nèi)土層分3層，厚度從上到下依次為30，70，100 m。模型介質(zhì)參數(shù)見表1。在地表x=100～4 900 m范圍內(nèi)選取55個(gè)觀測(cè)點(diǎn)，其中x=500～700 m，x=4 300～4 500 m范圍內(nèi)觀測(cè)點(diǎn)間隔為40 m，其余間隔均為100 m。盆地內(nèi)最小網(wǎng)格尺寸為2 m，基巖最大網(wǎng)格尺寸為10 m。模擬時(shí)間步長(zhǎng)0.001 s。為了與盆地模型對(duì)比來研究盆地的放大效應(yīng)，建立與圖1（a）各觀測(cè)點(diǎn)處的土層分布和厚度完全相同的一維水平成層模型作為參考模型，共計(jì)55個(gè)一維土層模型。以測(cè)點(diǎn)7和12為例，其對(duì)應(yīng)的一維土層模型分別如圖1（b）、（c）所示。各土層參數(shù)與盆地模型相同，進(jìn)而模擬相同地震波輸入下一維土層的反應(yīng)。采用寬度0.4 s的SV脈沖波作為基巖輸入地震波，從模型左側(cè)傾斜入射，設(shè)地震波入射角（入射方向與垂直方向的夾角）為θ，其位移時(shí)程曲線及其傅里葉譜見圖2。

圖2 輸入波時(shí)程及其傅里葉譜Fig.2 Time history of input wave and its spectrum

表1 計(jì)算模型介質(zhì)參數(shù)Table 1 Physical parameters of the calculation model

對(duì)于該盆地模型，采用剪切波速加權(quán)平均法計(jì)算盆地最深處的等效一維土層模型（土層分布如圖1（c））的基階自振頻率f0，計(jì)算公式為：

圖1 計(jì)算模型示意圖Fig.1 Schematic of the calculation model

式中：vsi和hi分別為各土層的剪切波速和厚度；H為土層總厚度（200 m）。根據(jù)式（6）計(jì)算得f0=0.7 Hz。2.2節(jié)中將基于該自振頻率將模擬結(jié)果對(duì)應(yīng)的頻率無量綱化。

2 計(jì)算結(jié)果

為研究地震波入射角對(duì)成層盆地地表地震動(dòng)放大特征的影響，利用圖1所示盆地模型和盆地不同場(chǎng)點(diǎn)處相同土層厚度的一維水平成層模型，在保持其余參數(shù)不變的基礎(chǔ)上，依次改變地震波的入射角θ為0°、10°、20°和30°，然后分別模擬得到相應(yīng)入射角下模型的地震動(dòng)響應(yīng)。需要指出的是，由于文中研究的是平面SV波入射下的盆地場(chǎng)地效應(yīng)，當(dāng)SV波的入射角大于一定角度時(shí)，反射波就變成了非均勻平面波［15］，不在文中討論范圍之內(nèi)。因此，文中入射SV波將具有臨界入射角，即SV波的入射角需小于按以下公式確定的角度：

式中，β為P波的反射角，當(dāng)β為90°時(shí)，入射角θ達(dá)到臨界值。由表1可知基巖半空間介質(zhì)Vs=1 650 m/s，Vp=2 858 m/s，代入式（7）計(jì)算可得SV波臨界入射角θ約為35o，因此文中最大入射角取為30°。

2.1 入射角對(duì)盆地地震動(dòng)時(shí)域放大特征的影響

不同入射角下盆地模型模擬得到的地表觀測(cè)點(diǎn)水平和垂直分量的位移時(shí)程如圖3所示?？梢钥闯?，盆地內(nèi)觀測(cè)點(diǎn)的位移時(shí)程相比基巖更加復(fù)雜、振蕩，地震動(dòng)持時(shí)顯著延長(zhǎng)。隨著入射角θ的增大，地表各觀測(cè)點(diǎn)直達(dá)體波的到時(shí)會(huì)有滯后的現(xiàn)象。水平分量直達(dá)體波的強(qiáng)度逐漸減弱，垂直分量上則逐漸增強(qiáng)。不同入射角度下盆地兩側(cè)邊緣均有明顯的次生面波產(chǎn)生，但隨入射角度的增大，盆地左側(cè)次生面波的強(qiáng)度逐漸大于右側(cè)。不同入射角度下盆地左側(cè)的面波強(qiáng)度有增強(qiáng)的趨勢(shì)，右側(cè)面波強(qiáng)度有減弱趨勢(shì)，垂直分量的次生面波強(qiáng)度受入射角的影響相對(duì)較小。

圖3 不同入射角下盆地地表觀測(cè)點(diǎn)水平（左）和垂直分量（右）的位移時(shí)程Fig.3 Displacement time histories of horizontal（left）and vertical（right）components of surface observation points corresponding to different incidence angles

圖4為不同入射角下盆地地表觀測(cè)點(diǎn)的位移峰值（PGD）分布。對(duì)于水平分量，當(dāng)入射角θ=0°即垂直入射時(shí)，地表位移峰值基本呈對(duì)稱分布。隨著入射角的增大，盆地左側(cè)邊緣觀測(cè)點(diǎn)的位移峰值普遍大于右側(cè)。盆地邊緣處觀測(cè)點(diǎn)（6～10號(hào)，45～50號(hào)觀測(cè)點(diǎn)）的PGD受入射角的影響十分顯著：6～10號(hào)觀測(cè)點(diǎn)范圍內(nèi)，PGD隨入射角的增加而增大；45～50號(hào)觀測(cè)點(diǎn)范圍內(nèi)，PGD隨入射角的增大而減小。盆地中央?yún)^(qū)域PGD隨入射角的增大而減小，但固定入射角下其值始終比較穩(wěn)定，可能和盆地較大的寬度有關(guān)。同時(shí)，入射角對(duì)PGD在盆地內(nèi)的分布具有重要影響，其中主要影響盆地邊緣效應(yīng)出現(xiàn)的位置。以盆地左側(cè)最強(qiáng)烈地震動(dòng)區(qū)域?yàn)槔?，?dāng)入射角θ=0°或10°時(shí)，PGD有先增大后至平穩(wěn)的趨勢(shì)，PGD最大值發(fā)生在x=900 m處；而入射角θ=20°或30°時(shí)，PGD有先增大后減小至平穩(wěn)的趨勢(shì)，PGD最大值發(fā)生在x=620 m處?？傮w而言，PGD最大值出現(xiàn)的位置有隨入射角的增大而向盆地左側(cè)移動(dòng)的趨勢(shì)。

圖4 不同入射角度下盆地地表觀測(cè)點(diǎn)的位移峰值分布Fig.4 PGD distributions of observation points on the surface of the basin under different incidence angles

對(duì)于垂直分量，與盆地角點(diǎn)相同距離的左側(cè)觀測(cè)點(diǎn)的位移峰值明顯大于右側(cè)，但最大PGD值兩側(cè)基本相同。在盆地左側(cè)和右側(cè)邊緣處（6～15號(hào)，40～50號(hào)觀測(cè)點(diǎn)范圍）均表現(xiàn)出明顯的盆地邊緣效應(yīng)，PGD在此處變化劇烈。PGD最大值與入射角之間具有明顯的規(guī)律性，即隨入射角的增大而增大。盆地中央?yún)^(qū)域的PGD也隨入射角度的增加而增大（10°～30°時(shí)增大最顯著），這與水平分量的趨勢(shì)相反。不同入射角下盆地左側(cè)和右側(cè)觀測(cè)點(diǎn)垂直分量的PGD最大值均出現(xiàn)在6號(hào)和50號(hào)觀測(cè)點(diǎn)處。

55個(gè)觀測(cè)點(diǎn)對(duì)應(yīng)的一維水平成層模型在不同角度SV波入射下的位移峰值分布如圖5所示?？梢钥闯?，同一模型（固定觀測(cè)點(diǎn)）相同厚度覆蓋土層下，隨著入射角的增大水平分量PGD呈減小趨勢(shì)，垂直分量PGD則呈增大趨勢(shì)。同一入射角度下，水平分量上覆蓋土層厚度越大，其PGD值也越大；垂直分量上在測(cè)點(diǎn)7處（覆蓋層厚度100 m）地震動(dòng)幅值最大。此外，垂直入射下水平成層模型垂直分量PGD為0，即垂直入射下，模型沒有產(chǎn)生垂直方向的振動(dòng)。由此可見，垂直入射下二維盆地模型產(chǎn)生的垂直方向上的位移是由盆地特殊的幾何形狀引起的。

圖5 不同入射角度下55個(gè)一維模型地表觀測(cè)點(diǎn)的位移峰值分布Fig.5 PGD distributions of surface observation points of 55 1D models under different incidence angles

定義盆地地震動(dòng)放大系數(shù)（AF）為盆地模型模擬得到的地表位移峰值與對(duì)應(yīng)場(chǎng)點(diǎn)處的一維水平成層模型模擬得到的地表位移峰值的比值，其分布如圖6所示，圖7為圖6（a）在盆地邊緣附近區(qū)域的局部放大?？梢钥闯?，對(duì)于水平分量，當(dāng)?shù)卣鸩ù怪比肷鋾r(shí)，盆地邊緣區(qū)域觀測(cè)點(diǎn)的地震動(dòng)放大系數(shù)大多小于1。隨著入射角度的增大，放大系數(shù)也逐漸增大，當(dāng)入射角θ=30°時(shí)，8號(hào)觀測(cè)點(diǎn)放大系數(shù)達(dá)到了最大值1.43。盆地中央?yún)^(qū)域水平分量地震動(dòng)放大系數(shù)穩(wěn)定在1左右。由于地震波從左側(cè)入射，盆地右側(cè)區(qū)域?qū)θ肷浣堑淖兓鄬?duì)不敏感，地震動(dòng)放大系數(shù)均小于1，且入射角越大，放大系數(shù)總體越小。此外，水平分量上盆地外側(cè)4～5號(hào)測(cè)點(diǎn)區(qū)域的放大系數(shù)也大于1，應(yīng)主要受盆地構(gòu)造的影響。對(duì)于垂直分量，放大系數(shù)隨入射角的增加急劇減小，其放大效應(yīng)在盆地邊緣區(qū)域（6號(hào)及50號(hào)測(cè)點(diǎn)處）最顯著，最大放大系數(shù)可達(dá)4.5。相比而言，盆地垂直分量的放大系數(shù)顯著大于水平分量。因此平面波入射下除考慮盆地構(gòu)造的影響外，還需考慮入射角度的影響。

圖6 不同入射角度下盆地地表觀測(cè)點(diǎn)的放大系數(shù)分布Fig.6 Distributions of amplification factor of surface observation points of the basin under different incidence angles

2.2 入射角對(duì)盆地地震動(dòng)頻域放大特征的影響

選取盆地內(nèi)6、8、10、12、14和20號(hào)觀測(cè)點(diǎn)（距盆地角點(diǎn)540～1 700 m范圍，位置如圖1所示），分別定義為點(diǎn)A、B、C、D、E和F。圖8為A～F觀測(cè)點(diǎn)在不同入射角下水平分量位移時(shí)程的傅里葉譜（本小節(jié)只研究水平分量地震動(dòng)），其中頻率f基于盆地最深處一維等效土層模型的自振頻率f0進(jìn)行了無量綱化。可以看出，對(duì)同一觀測(cè)點(diǎn)，不同入射角度下其頻譜特征基本相似，但卓越頻率（傅里葉譜最大幅值對(duì)應(yīng)的頻率）隨入射角度的增大而略有增加。從盆地邊緣到盆地內(nèi)部，觀測(cè)點(diǎn)的卓越頻率有從高頻向低頻移動(dòng)的趨勢(shì)，如入射角θ=30°時(shí)，A～F號(hào)觀測(cè)點(diǎn)的卓越頻率分別等于2.87f0、1.91f0、1.93f0、1.31f0、1.13f0、1.39f0。此外，A～E點(diǎn)傅氏譜幅值隨入射角增大而增大，而位于盆地中央的F號(hào)觀測(cè)點(diǎn)的傅氏譜幅值對(duì)入射角的變化相對(duì)不敏感，且其幅值小于盆地邊緣處觀測(cè)點(diǎn)。

圖8 不同入射角度下A～F號(hào)觀測(cè)點(diǎn)的傅里葉譜Fig.8 Fourier spectra of observation points A～F under different incidence angles

圖8（續(xù)）Fig.8（Continued）

圖9給出了A～F號(hào)觀測(cè)點(diǎn)在不同入射角下位移時(shí)程傅里葉譜比曲線。本小節(jié)定義頻域放大倍數(shù)（AFF）為盆地模型觀測(cè)點(diǎn)的位移時(shí)程傅里葉譜與一維土層模型位移時(shí)程傅里葉譜之比，用來定量描述頻域內(nèi)盆地場(chǎng)地相比一維土層模型的附加放大效應(yīng)。可以看出，地震波從左側(cè)斜入射時(shí)，入射角越大，盆地地表觀測(cè)點(diǎn)的譜比值越大。如觀測(cè)點(diǎn)A在4種入射角下的最大放大倍數(shù)分別為2.12、2.38、2.65和2.97，這與時(shí)域結(jié)果對(duì)應(yīng)。盆地邊緣處的放大倍數(shù)要大于盆地其它區(qū)域，在C號(hào)觀測(cè)點(diǎn)（x=700 m）處，頻域放大倍數(shù)達(dá)到了最大值3.5，而盆地內(nèi)部F號(hào)觀測(cè)點(diǎn)的頻域放大系數(shù)總體最小。值得注意的是，A點(diǎn)在f0～2.5f0頻段內(nèi)的AFF值基本都小于1，主要放大頻段集中分布在3f0附近；B～E觀測(cè)點(diǎn)的譜比曲線整體呈雙峰分布，一個(gè)大致位于f0～2.5f0的范圍（如C點(diǎn)和E點(diǎn)在此頻段內(nèi)的AFF最大值分別達(dá)到了2.26和2.12），另一個(gè)（也是最主要放大頻段）集中在4f0左右；而F點(diǎn)在較寬頻帶范圍內(nèi)的AFF值均較為均衡。表明相比等效的一維土層模型，盆地邊緣土層較薄區(qū)域?qū)ο鄬?duì)“高頻”的地震動(dòng)放大較顯著，且主要放大頻段集中，而靠近盆地內(nèi)部土層較厚，在較寬的頻帶范圍內(nèi)均表現(xiàn)出明顯放大效應(yīng)。此外，固定觀測(cè)點(diǎn)處譜比曲線的分布特征相近，基本不受入射角度的影響。

圖9 不同入射角度下A～F觀測(cè)點(diǎn)的譜比曲線Fig.9 Fourier spectra ratio curves of observation points A～F under different incidence angles

不同入射角下盆地內(nèi)6～50號(hào)觀測(cè)點(diǎn)譜比最大值分布如圖10所示，可以看出譜比最大值A(chǔ)FFmax與入射角之間有很好的相關(guān)性?？傮w來說，盆地邊緣區(qū)域的AFFmax大于盆地其它區(qū)域，從6號(hào)點(diǎn)到35號(hào)觀測(cè)點(diǎn)，AFFmax有隨入射角增大而增大的趨勢(shì)，但35號(hào)觀測(cè)點(diǎn)之后（逐漸靠近盆地右側(cè)邊界），AFFmax卻隨入射角度的增大而減小。

圖10 不同入射角度下盆地內(nèi)觀測(cè)點(diǎn)的譜比最大值分布Fig.10 Distributions of the maximum spectral ratio of observation points in the basin under different incidence angles

3 成層盆地與等效均勻盆地地震動(dòng)及其放大特征比較

考慮盆地內(nèi)介質(zhì)分層對(duì)不同入射角度的影響，建立和圖1（a）相同尺寸的均勻盆地模型，即盆地內(nèi)土層為均一介質(zhì)，其厚度仍為200 m。剪切波速取成層盆地模型的等效剪切波速［16］，即

進(jìn)而模擬相同入射波（圖2）下均勻盆地的地震響應(yīng)與放大效應(yīng)。

圖11給出了均勻盆地模型地表觀測(cè)點(diǎn)的水平分量位移峰值及相比一維等效土層模型的放大系數(shù)分布?？梢钥闯?，相比成層盆地，均勻盆地模型盆地邊緣效應(yīng)遠(yuǎn)小于成層盆地，其PGD最大值為0.374 m，而成層盆地模型PGD最大值為0.619 m，且地震動(dòng)最大值的位置基本不隨入射角的變化而改變，恒定在12號(hào)觀測(cè)點(diǎn)處，這與成層盆地模型入射角≤10°時(shí)的PGD最大值位置相同。此外，盆地內(nèi)PGD最大值隨入射角的增大而一致減小。從峰值放大系數(shù)來看，在給定入射波下，入射角度較小時(shí)（≤10°）成層盆地內(nèi)部的放大系數(shù)（圖6（a））小于等效均勻盆地，入射角大于10°后則大于等效均勻盆地，如均勻盆地模型在30°入射角時(shí)盆地內(nèi)放大系數(shù)最大值為1.14（顯著小于成層盆地模型的1.45），出現(xiàn)在11號(hào)測(cè)點(diǎn)附近。同時(shí)，均勻盆地模型的放大系數(shù)分布曲線相比成層模型結(jié)果簡(jiǎn)單，盆地內(nèi)最大放大系數(shù)相近，基本不隨入射角的變化而改變。

圖11 不同入射角度下均勻盆地模型地表觀測(cè)點(diǎn)的位移峰值及放大系數(shù)分布Fig.11 Distributions of PGD and amplification factor of observation points on the surface of the single layer basin

4 結(jié)論

文中采用有限元法結(jié)合透射邊界模擬SV波斜入射下成層盆地的地震動(dòng)響應(yīng)。通過改變地震波的入射角度，研究了入射角對(duì)成層盆地相比等效一維土層模型的時(shí)、頻域放大特征的影響，得到以下結(jié)論：

（1）入射角θ對(duì)盆地地表地震動(dòng)峰值影響顯著。水平分量上，盆地邊緣區(qū)域的位移峰值隨θ增加而增大，中間及右側(cè)區(qū)域的位移峰值隨之減?。淮怪狈至可?，盆地所有區(qū)域的位移峰值均隨θ增加而增大。

（2）入射角度對(duì)盆地邊緣效應(yīng)區(qū)域的位置有影響。對(duì)水平分量，θ大于20°時(shí)出現(xiàn)在盆地斜邊區(qū)域范圍內(nèi)，小于20°時(shí)則位于盆地斜邊范圍右側(cè)的區(qū)域；對(duì)垂直分量，邊緣效應(yīng)始終出現(xiàn)在緊鄰盆地角點(diǎn)的位置，不受入射角的影響。

（3）不同入射角度下，盆地左側(cè)邊緣區(qū)域水平分量地震動(dòng)放大系數(shù)顯著大于右側(cè)，其值在此區(qū)域達(dá)到最大，且隨θ增大而增大（最大放大系數(shù)1.43），右側(cè)則反之；垂直分量上左右兩側(cè)放大系數(shù)基本相同，其值隨入射角的增大而減小，最大放大系數(shù)4.5左右。兩分量上放大系數(shù)的分布特征基本不受入射角度的影響，最大放大系數(shù)的位置與最大位移峰值位置基本相同。

（4）固定觀測(cè)點(diǎn)在不同入射角度下的頻譜特征類似，但峰值頻率有隨θ的增大略變大的趨勢(shì)。入射角變化時(shí)不同位置觀測(cè)點(diǎn)的譜比曲線分布特征相似，但譜比值均隨入射角度的增大而增大。不同θ角下成層盆地邊緣區(qū)域以相對(duì)“高頻”的放大為主，放大頻段較為集中，向盆地內(nèi)部則在較寬頻帶范圍內(nèi)均有放大。

（5）盆地內(nèi)介質(zhì)波速的縱向變化對(duì)其地震動(dòng)分布影響顯著。盆地內(nèi)分層介質(zhì)與盆地斜邊的相互作用可能是導(dǎo)致大角度入射時(shí)成層盆地地震動(dòng)放大效應(yīng)相比等效均勻盆地顯著的主要原因，因此在斜入射下盆地地震動(dòng)的估計(jì)中，必須同時(shí)考慮盆地內(nèi)介質(zhì)橫向和縱向波速的變化。