軟土層幾何特性與剪切波速對(duì)場(chǎng)地峰值加速度的影響1

劉方成 武景芳 陳 斌 趙成齊

（湖南工業(yè)大學(xué)土木工程學(xué)院，株洲 412008）

軟土層幾何特性與剪切波速對(duì)場(chǎng)地峰值加速度的影響1

劉方成 武景芳 陳 斌 趙成齊

（湖南工業(yè)大學(xué)土木工程學(xué)院，株洲 412008）

軟土層對(duì)場(chǎng)地地震動(dòng)的影響一直以來是地震工程學(xué)的研究重點(diǎn)。本文應(yīng)用一維真非線性場(chǎng)地反應(yīng)分析方法，對(duì)某單層勻質(zhì)場(chǎng)地內(nèi)軟土層的幾何特征（厚度與埋深）和剪切波速變化對(duì)地面峰值加速度的影響進(jìn)行了數(shù)值分析。研究結(jié)果表明：軟土層的存在使得場(chǎng)地加速度幅值分布在軟土層處發(fā)生突變，隨著軟土層幾何特性和剪切波速的不同，這種突變可能使得地面峰值加速度增大或減??；存在一個(gè)由軟土層厚度、埋深和相對(duì)剪切波速三個(gè)參數(shù)構(gòu)成的臨界狀態(tài)面，當(dāng)軟土層的狀態(tài)位于臨界面以內(nèi)時(shí)，軟土層對(duì)地面峰值加速度起放大效應(yīng)，且加速度效應(yīng)系數(shù)隨著上述三個(gè)變量的增大表現(xiàn)出先增大后減小的變化規(guī)律；當(dāng)軟土層位于臨界面以外時(shí)，其對(duì)地面峰值加速度起衰減效應(yīng)，加速度效應(yīng)系數(shù)隨著上述三個(gè)變量的增大而減小。根據(jù)算例的參數(shù)分析建立了其臨界面方程，并提出了估計(jì)軟土層加速度效應(yīng)系數(shù)的經(jīng)驗(yàn)公式，可為類似問題提供參考。

軟弱土層 場(chǎng)地效應(yīng) 土動(dòng)力特性 場(chǎng)地反應(yīng)分析 地面峰值加速度

劉方成，武景芳，陳斌，趙成齊，2015.軟土層幾何特性與剪切波速對(duì)場(chǎng)地峰值加速度的影響.震災(zāi)防御技術(shù)，10（4）：840—852. doi：10.11899/zzfy20150403

引言

軟夾層對(duì)場(chǎng)地反應(yīng)的影響一直是工程地震學(xué)的研究重點(diǎn)。已有的地震災(zāi)害調(diào)查和場(chǎng)地分析（謝君斐等，1978；錢勝國(guó)，1994；黃潤(rùn)秋等，2003；陳國(guó)興等，2004；史丙新等，2010；金丹丹等，2012；張建毅等，2013）均表明，軟弱夾層的存在對(duì)地表地震反應(yīng)存在較顯著的影響，主要表現(xiàn)為：①軟夾層對(duì)地表加速度幅值隨夾層埋深的不同起到衰減或放大的效果；②含軟弱夾層場(chǎng)地的地表加速度反應(yīng)譜呈多峰值形狀，特征周期延長(zhǎng)。我國(guó)現(xiàn)行抗震規(guī)范尚不能敏感地反映軟弱土層對(duì)場(chǎng)地分類的影響（盧華喜等，2010），可能給工程結(jié)構(gòu)的抗震設(shè)計(jì)帶來不確定性。為修正規(guī)范的不足，已有學(xué)者在這方面開展了前瞻性的工作，并對(duì)軟土夾層影響場(chǎng)地效應(yīng)的機(jī)理進(jìn)行了探索（薄景山等，2003a；2003b；2003c；2004；黃玉龍等，2000）?；谝延醒芯?，本文采用能真實(shí)考慮場(chǎng)地土滯回阻尼特性的非線性模型進(jìn)行真非線性分析，考查了軟土層幾何構(gòu)造（厚度和埋深）和剪切波速變化對(duì)地表加速度反應(yīng)幅值的影響規(guī)律，并基于參數(shù)分析，提出了軟土層加速度效應(yīng)系數(shù)的概念及其經(jīng)驗(yàn)估計(jì)公式，為后續(xù)研究和規(guī)范修訂提供參考。

1 計(jì)算模型

1.1 含軟弱土層場(chǎng)地的構(gòu)造

鑒于多層土的復(fù)雜性，很難看出某一參數(shù)對(duì)地面加速度反應(yīng)幅值的影響規(guī)律，本文從單層均質(zhì)場(chǎng)地入手探討軟弱土層的地震效應(yīng)。假設(shè)有一單層勻質(zhì)場(chǎng)地（II類場(chǎng)地），覆蓋層總厚度H=40m，場(chǎng)地剪切波速νs=250m/s，場(chǎng)地土質(zhì)量密度ρ=2.0×103kg/m3，其中含有一層軟弱土層，軟弱土層的剪切波速為νss，質(zhì)量密度為ρi=1.6×103kg/m3，軟土層在場(chǎng)地中的幾何分布特性由其厚度hL、頂部埋深hT確定。含軟土層的場(chǎng)地構(gòu)造示意圖及相關(guān)參數(shù)取值如圖1所示。

1.2 分析程序與計(jì)算參數(shù)

計(jì)算分析在MATLAB語(yǔ)言環(huán)境下編制的場(chǎng)地非線性分析程序（FSRAP）中進(jìn)行。圖2給出了程序所采用的一維有限元分析模型（欒茂田等，1992），圖中hi、ρi、Gi、Di分別表示第i層子單元的厚度、密度、剪切模量與阻尼比。一維有限元模型各子單元的剪切模量Gi和滯回阻尼比Di隨加載時(shí)程的變化，由基于阻尼的滯回模型DBM（尚守平等，2007）來模擬。圖3給出了FSRAP程序?qū)HAKE（Idriss等，1992）、NERA（Bardet等，2001）算例的分析結(jié)果與SHAKE、NERA計(jì)算結(jié)果的對(duì)比情況?？梢姡現(xiàn)SRAP程序與SHAKE、NERA計(jì)算結(jié)果吻合較好，且FSRAP與NERA的結(jié)果更為接近，這是由于FSRAP與NERA均為真非線性分析方法，而SHAKE為等效線性方法。而FSRAP與NERA之間的差異則是由于Masing滯回模型不能真實(shí)模擬場(chǎng)地土滯回阻尼所致。算例對(duì)比分析表明，F(xiàn)SRAP程序具有較好的可靠性。

圖1 含軟弱（夾）層場(chǎng)地剖面構(gòu)造示意圖Fig. 1 Construction diagram of the profile containing a soft soil layer

圖2 FSRAP一維有限元分析模型Fig. 2 One dimensional analysis model for layered site in the program FSRAP

為盡量減少其他因素的影響，選取Seed等（1970；1986）建議的各類土的經(jīng)驗(yàn)動(dòng)力特性曲線（如圖4中粗黑線所示）作為場(chǎng)地土及軟土共同的動(dòng)剪模量曲線和阻尼比曲線。分別選取了Taft波、El-centro波、Loma Pretia波作為基巖輸入地震波，輸入加速度幅值根據(jù)不同工況分別調(diào)整為0.05g、0.1g和0.2g。不同輸入地震波分析結(jié)果具有相似的規(guī)律，限于篇幅文中僅給出Taft波的計(jì)算結(jié)果。輸入地震波特征如圖5所示。

圖3 FSRAP與SHAKE、NERA算例分析對(duì)比Fig. 3 Comparison of the results from different programs of FSRAP， SHAKE and NERA

圖4 本文計(jì)算所采用的場(chǎng)地土動(dòng)力特性曲線Fig. 4 Dynamic properties curves for selected soil in this study

圖5 基巖輸入地震波時(shí)程曲線及FFT曲線（Taft波）Fig. 5 The time history and FFT property of the earthquake wave used as input in the rockbed （Taft wave）

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

2.1 軟土層對(duì)加速度幅值分布曲線的影響

軟土層對(duì)加速度幅值分布曲線的影響如圖6所示。其中，圖6（a）、（b）給出了當(dāng)軟土層厚度分別為2m和4m時(shí)，軟土層埋深變化對(duì)加速度幅值沿深度分布曲線的影響?？梢姡孩偌铀俣确磻?yīng)幅值分布曲線在軟土層位置處發(fā)生突變，軟土層埋深越大，加速度幅值分布曲線突變?cè)矫黠@；②當(dāng)軟弱夾層埋深較小時(shí)，軟弱夾層的存在使得場(chǎng)地地面加速度幅值相比于同條件下不含軟土層場(chǎng)地（以下簡(jiǎn)稱“原場(chǎng)地”）的地面加速度幅值增大；③當(dāng)軟弱夾層埋深大于某一界限值后，將使得地面加速度幅值相比于原場(chǎng)地減小，且軟土層埋深越大，衰減越多；④當(dāng)埋深較小時(shí)，軟土層厚度變化對(duì)接近地表的土層加速度幅值存在顯著影響；當(dāng)埋深較大時(shí)，軟土層厚度對(duì)場(chǎng)地加速度幅值分布曲線的影響較小。如當(dāng)埋深為4m時(shí)，厚度為2m的軟土層使得地面加速度幅值相對(duì)于原場(chǎng)地增大；而厚度為4m的軟土層使得地面加速度幅值相對(duì)于原場(chǎng)地減小。

圖6 深度一定時(shí)軟弱夾層埋深變化對(duì)場(chǎng)地加速度幅值分布曲線的影響Fig. 6 Influence of the buried depth of a soft soil layer with defined thickness on profile of acceleration amplitude

圖6（c）給出了在軟土層的剪切波速和埋深一定時(shí)（νss=100m/s，hL=12m），場(chǎng)地內(nèi)加速度幅值分布曲線隨軟土層厚度的變化情況?？梢娫诩羟胁ㄋ俸吐裆钜欢ǖ那闆r下，隨著軟土層厚度的增大，場(chǎng)地加速度幅值分布曲線的突變也越大，由于這種突變使得地面峰值加速度隨著軟土層厚度的增大而減小。

圖6（d）給出了當(dāng)軟土層厚度和埋深一定時(shí)（hL=2m，hT=12m），軟土層剪切波速變化對(duì)場(chǎng)地反應(yīng)的影響?？梢姡孩匐S著軟土層剪切波速的減小，地面加速度反應(yīng)逐漸減?。虎诋?dāng)軟土層剪切波速相比于原場(chǎng)地土剪切波速降低不多時(shí)（如軟土層剪切波速為νss=200m/s，vss=νss/νs=0.8），場(chǎng)地加速度反應(yīng)幅值沿場(chǎng)地深度分布曲線在軟土層處向增大方向發(fā)生突變、地面加速度幅值相對(duì)于原場(chǎng)地略微增大；③當(dāng)軟土層剪切波速相比于原場(chǎng)地土剪切波速降低較多時(shí)（如軟土層剪切波速小于等于150m/s，vss≤0.6），場(chǎng)地加速度反應(yīng)幅值沿場(chǎng)地深度分布曲線在軟土層處向減小方向發(fā)生突變、地面加速度幅值相對(duì)于原場(chǎng)地減小，且軟土層剪切波速越小，地面加速度幅值相對(duì)于原場(chǎng)地減小得越多。

2.2 軟土層對(duì)地面峰值加速度的影響

圖7給出了含軟土層場(chǎng)地的地面峰值加速度與原場(chǎng)地的地面峰值加速度的比值（以下簡(jiǎn)稱軟土層的加速度效應(yīng)系數(shù)）隨軟土層特性變化的規(guī)律。

由圖7可見：①當(dāng)軟土層埋深較小時(shí)，加速度效應(yīng)系數(shù)大于1，且在一定范圍內(nèi)加速度效應(yīng)系數(shù)隨著埋深的增大先增大后減小，當(dāng)軟土層埋深超過某臨界值以后，加速度效應(yīng)系數(shù)小于1，且隨著軟土層埋深的增加而減?。虎谲浲翆雍穸仍酱?，其對(duì)地面峰值加速度的影響越大，表現(xiàn)為厚度大的曲線兩端均離圖中灰色虛線（代表效應(yīng)系數(shù)為1）越遠(yuǎn)；③軟土層剪切波速越小，對(duì)地面加速度的衰減效應(yīng)越明顯；④軟土層對(duì)地面加速度幅值的影響隨著基巖輸入加速度的增大而更顯著。

圖7 軟弱夾層埋深及厚度對(duì)軟土層加速度效應(yīng)系數(shù)的影響Fig. 7 Effect of thickness and depth of the inter soft soil layer on the ground acceleration amplification factor

圖8進(jìn)一步給出了在不同基巖輸入地震波幅值下，軟土層的加速度效應(yīng)系數(shù)隨軟土層剪切波速變化的規(guī)律，圖中軟土層厚度4m，埋深從0m到36m變化。由圖可見：①當(dāng)軟土層埋深較小時(shí)，軟土層加速度效應(yīng)系數(shù)隨著軟土層剪切波速的減小而先增大后減小，加速度效應(yīng)系數(shù)會(huì)出現(xiàn)大于1的情況；②當(dāng)軟土層埋深大于某一界限深度時(shí)，軟土層的加速度效應(yīng)系數(shù)始終小于1，且隨著軟土層剪切波速的減小而減小；③當(dāng)vs/vss在1—2之間時(shí)，較大基巖輸入加速度幅值下的軟土層加速度效應(yīng)系數(shù)更小，說明基巖輸入加速度增大會(huì)使得原本較硬的軟土層變軟，其軟土夾層效應(yīng)越明顯。

圖8 軟土夾層剪切波速對(duì)其加速度效應(yīng)系數(shù)的影響Fig. 8 Effect of the shear wave velocity of soft layer on ground acceleration amplification factor

圖9給出了軟土層加速度效應(yīng)系數(shù)隨軟土層埋深、厚度、相對(duì)剪切波速變化的三維曲面。圖9（a）為hL-hT平面上的情況；圖9（b）為vss-hT平面上的情況；圖9（c）為vss-hL平面上的情況。圖中淺色表示加速度效應(yīng)系數(shù)大于1的區(qū)域，即軟土層對(duì)地面加速度起放大效應(yīng)；深色表示加速度效應(yīng)系數(shù)小于1的區(qū)域，即軟土層對(duì)地面加速度起衰減效應(yīng)。各圖底部平面的投影給出了軟土層的加速度效應(yīng)在相應(yīng)平面上由放大到衰減的臨界曲線。

由圖9（a）可見，軟土層加速度效應(yīng)系數(shù)隨軟土層埋深和厚度的增大表現(xiàn)出先增大后減小的趨勢(shì)。隨著軟土層剪切波速的不同，hL-hT平面上臨界曲線的形狀和位置也發(fā)生變化。軟土層剪切波速越小，代表加速度效應(yīng)系數(shù)大于1的淺灰色面積也就越小。

由圖9（b）可見，軟土層加速度效應(yīng)系數(shù)隨軟土層埋深和相對(duì)剪切波速的增大表現(xiàn)出先增大后減小的趨勢(shì)。隨著軟土層厚度的不同，vss-hT平面上各曲面的臨界線形狀和位置而變化，軟土層厚度越大，臨界線所圍放大區(qū)域越小。

由圖9（c）可見，軟土層加速度效應(yīng)系數(shù)隨軟土層厚度和相對(duì)剪切波速的增大表現(xiàn)出先增大后減小的趨勢(shì)。vss-hL平面上的臨界線形狀和位置隨著軟土層埋深的不同而變化，隨著軟土層埋深的增加，臨界線所圍放大區(qū)域迅速減小。值得注意的是，當(dāng)軟土層埋深等于0時(shí)（即軟土層位于場(chǎng)地表層），軟土層對(duì)地面加速度幅值的影響以放大效應(yīng)為主。

2.3 軟土層對(duì)地面峰值加速度影響效應(yīng)的臨界曲面

由前面的計(jì)算結(jié)果可知，軟土層對(duì)地面峰值加速度的效應(yīng)系數(shù)受軟土層厚度hL、埋深hT和剪切波速νss三個(gè)參數(shù)的共同影響，且可能存在一個(gè)空間臨界面，當(dāng)軟土層的上述參數(shù)位于該臨界面內(nèi)外時(shí)，軟土層對(duì)地面加速度分別起放大和衰減效應(yīng)。要判斷軟土層的場(chǎng)地效應(yīng)，應(yīng)確定該臨界面。

圖10（a）給出了由圖9中各曲面投影所得到的hL-hT平面上的臨界曲線，圖11（a）則給出了由圖10中各曲面投影所得到的hT-vss平面上的臨界曲線?？梢?，hL-hT平面上的臨界曲線接近于直線，而hT-vss平面上的臨界曲線則近似為雙曲線。從應(yīng)用簡(jiǎn)便出發(fā)，提出如下臨界面方程：

式中，hL表示軟土層厚度；hT表示軟土層埋深；vss=νss/νs表示軟土層的相對(duì)剪切波速；νss為軟土層的剪切波速；νs為原場(chǎng)地的剪切波速；A0、B為曲線參數(shù)。當(dāng)Fcr＜0時(shí)，軟土層對(duì)地面峰值加速度起放大效應(yīng)；當(dāng)Fcr＞0時(shí)，地面峰值加速度起衰減效應(yīng)。

圖9 軟土層加速度效應(yīng)系數(shù)隨軟土層厚度、埋深和相對(duì)剪切波速變化的規(guī)律Fig. 9 Effect of the thickness， buried depth and relative shear velocity of a soft layer on ground acceleration amplification factor

用式（1）對(duì)圖10（a）、圖11（a）中的臨界曲線進(jìn)行擬合，得到兩個(gè)曲線參數(shù)的取值為：A0=13.5，B=0.5。擬合曲線分別繪制于圖10（b）、圖11（b）中，可見擬合的臨界曲線可合理地代表計(jì)算臨界線，且偏于安全（擬合曲線所給出的放大效應(yīng)區(qū)域相對(duì)于計(jì)算曲線偏大）。

圖10 軟土層對(duì)地面加速度影響效應(yīng)的臨界曲線（hL-hT平面）Fig. 10 Critical curves representing the influence of soft soil layer on ground acceleration response （in the hL-hTplane）:

圖11 軟土層對(duì)地面加速度影響效應(yīng)的臨界曲線（hT-vss平面）Fig. 11 Critical curves representing the influence of soft soil layer on ground acceleration response （in the hT-vSSplane）

2.4 軟土層對(duì)地面加速度影響效應(yīng)的經(jīng)驗(yàn)估計(jì)

已有對(duì)軟土層場(chǎng)地效應(yīng)的研究往往只限于定性的分析，且只強(qiáng)調(diào)了軟土層對(duì)場(chǎng)地加速度反應(yīng)的隔震效應(yīng)而忽略了在某些情況下的放大效應(yīng)。本文研究表明，隨著軟土層幾何構(gòu)造和剪切波速的不同，軟土層可能對(duì)地面加速度峰值起衰減效應(yīng)，但也可能起放大效應(yīng)，如果忽略后一種情況，就有可能做出偏于不安全的預(yù)測(cè)。同時(shí)為了在場(chǎng)地分類、地震小區(qū)劃等工程實(shí)踐中更真實(shí)地考慮軟土層的影響，有必要對(duì)軟土層的加速度效應(yīng)系數(shù)進(jìn)行定量分析研究。

根據(jù)對(duì)本文計(jì)算結(jié)果的分析，軟土層對(duì)地面峰值加速度的影響效應(yīng)系數(shù)應(yīng)該與臨界面相關(guān)，在臨界面以內(nèi)，效應(yīng)系數(shù)大于1；在臨界面以外，效應(yīng)系數(shù)可以取小于1的值，以考慮軟土層隔震的有利影響。本文提出軟土層對(duì)地面加速度的影響效應(yīng)系數(shù)估計(jì)公式如下：

式中，Ca表示軟土層對(duì)地面加速度的影響效應(yīng)系數(shù)；hL、hT、vss分別表示軟土層的厚度、埋深和相對(duì)剪切波速；hLB、hTB分別為由臨界面方程確定的軟土層臨界厚度和臨界埋深，由式（1）推得：

式中，A0、B為臨界曲面參數(shù)，取值亦與式（1）相同；p1=0.01，p2=0.5為曲線參數(shù)。

圖12給出了當(dāng)軟土層剪切波速νss為100m/s時(shí)，由式（2）所估計(jì)的軟土層加速度效應(yīng)系數(shù)隨軟土層幾何構(gòu)造特性的變化曲線與計(jì)算曲線的對(duì)比。圖12（a）、（c）表示了不同軟土層埋深下的加速度效應(yīng)系數(shù)隨軟土層厚度變化的曲線；圖12（b）、（d）表示了不同軟土層厚度下的加速度效應(yīng)系數(shù)隨軟土層埋深變化的曲線?？梢娛剑?）較好地反映了軟土層加速度效應(yīng)系數(shù)隨軟土層厚度和埋深變化的規(guī)律，即當(dāng)軟土層厚度和埋深均很小時(shí)，效應(yīng)系數(shù)大于1且隨厚度和埋深的增大表現(xiàn)出先增大后減??；當(dāng)軟土層厚度和埋深超過臨界值后，效應(yīng)系數(shù)小于1且隨著厚度和埋深的增大而降低。由圖可見，在相同條件下的估計(jì)曲線比計(jì)算曲線偏高，說明應(yīng)用式（2）估計(jì)軟土層的場(chǎng)地效應(yīng)是偏于安全的。

對(duì)圖13第一、二排四幅圖進(jìn)一步給出了在不同軟土層剪切波速下，由式（2）所估計(jì)的效應(yīng)系數(shù)隨軟土層厚度和埋深變化的三維曲面與計(jì)算曲面的對(duì)比，可見在不同的軟土層剪切波速下，估計(jì)曲面均較好地模擬了計(jì)算曲面的變化規(guī)律，且估計(jì)值比計(jì)算值高。

根據(jù)本文計(jì)算結(jié)果，軟土層位于臨界面之內(nèi)時(shí)的地面加速度效應(yīng)系數(shù)大于1但最大值亦未超過1.2，而在臨界面外的效應(yīng)系數(shù)最小值接近0.1，直接應(yīng)用式（2）估計(jì)軟土層的效應(yīng)系數(shù)偏于保守。為更合理地預(yù)測(cè)軟土層的加速度效應(yīng)，對(duì)式（2）進(jìn)一步修正如下：

式中，F(xiàn)cr表示臨界曲面，其方程由式（2）代表；Ca，max表示軟土層加速度效應(yīng)系數(shù)最大值，根據(jù)本文計(jì)算結(jié)果，可取為1.2；其他各參數(shù)意義同式（2）。上式的物理意義為，當(dāng)軟土層的幾何構(gòu)造和剪切波速參數(shù)位于臨界面以內(nèi)時(shí)，軟土層對(duì)地面加速度起放大效應(yīng)，此時(shí)為保守起見對(duì)軟土層的加速度效應(yīng)系數(shù)取最大值，因此臨界面內(nèi)的效應(yīng)系數(shù)為一平臺(tái)面；當(dāng)軟土層的幾何構(gòu)造和剪切波速參數(shù)位于臨界面以外時(shí)，軟土層對(duì)地面加速度起衰減效應(yīng)，對(duì)軟土層加速度效應(yīng)系數(shù)取一個(gè)不大于1的估計(jì)值。

圖12 軟土層對(duì)地面加速度影響效應(yīng)系數(shù)隨軟土層幾何構(gòu)造參數(shù)的變化關(guān)系曲線Fig. 12 Relationship between the influence factor of soft soil layer on PGA and geometric parameters of the soft layer

作為對(duì)比，式（4）所代表的加速度效應(yīng)系數(shù)隨軟土層幾何構(gòu)造參數(shù)變化的曲線和三維曲面也同時(shí)分別繪制于圖13。由圖可見，相比于式（2），修正后的式（4）能更好地模擬計(jì)算值，由式（4）所得到的加速度效應(yīng)系數(shù)曲面與計(jì)算曲面更為吻合。

圖13 軟土層加速度效應(yīng)系數(shù)估計(jì)值與計(jì)算值對(duì)比Fig. 13 Comparison of the influence factor of acceleration response between estimated and calculated

3 結(jié)論

本文應(yīng)用一維真非線性場(chǎng)地反應(yīng)分析方法，對(duì)某單層勻質(zhì)場(chǎng)地內(nèi)軟土層的幾何特性和剪切波速變化對(duì)地面加速度幅值的影響進(jìn)行了數(shù)值分析，得到了如下結(jié)論：

（1）軟土層的存在使得場(chǎng)地加速度幅值分布在軟土層處發(fā)生突變，隨著軟土層幾何特性和剪切波速的不同，這種突變可能使得地面峰值加速度增大或減小。

（2）軟土層對(duì)地面峰值加速度的影響受幾何特性和剪切波速的共同影響，存在一個(gè)由厚度hL、埋深hT和相對(duì)剪切波速ss三個(gè)參數(shù)構(gòu)成的臨界狀態(tài)函數(shù)Fcr（hL，hT，ss），當(dāng)Fcr＜0時(shí)，軟土層對(duì)地面峰值加速度起放大效應(yīng)，且加速度效應(yīng)系數(shù)隨著上述三個(gè)變量的增大表現(xiàn)出先增大后減小的變化規(guī)律；當(dāng)Fcr＞0時(shí)，軟土層對(duì)地面峰值加速度起衰減效應(yīng)，加速度效應(yīng)系數(shù)隨著軟土層的厚度、埋深和相對(duì)剪切波速的增大而減小。

（3）提出了可定量估計(jì)軟土層加速度效應(yīng)系數(shù)的經(jīng)驗(yàn)公式，與計(jì)算結(jié)果的對(duì)比驗(yàn)證表明，經(jīng)驗(yàn)公式能合理地反映不同幾何特性和剪切波速下，軟土層對(duì)地面峰值加速度的影響，可為相關(guān)應(yīng)用提供參考。

需要指出的是，本文經(jīng)驗(yàn)公式僅針對(duì)勻質(zhì)場(chǎng)地，對(duì)于復(fù)雜土層結(jié)構(gòu)中軟弱夾層的場(chǎng)地效應(yīng)的適用性還有待進(jìn)一步深入研究。

薄景山，李秀領(lǐng)，劉紅帥，2003a.土層結(jié)構(gòu)對(duì)地表加速度峰值的影響.地震工程與工程震動(dòng)，23（3）：35—40.

薄景山，李秀領(lǐng)，劉德東等，2003b.土層結(jié)構(gòu)對(duì)反應(yīng)譜平臺(tái)值的影響.地震工程與工程震動(dòng)，23（4）：29—33.

薄景山，李秀領(lǐng)，劉德東等，2003c.土層結(jié)構(gòu)對(duì)反應(yīng)譜特征周期的影響.地震工程與工程震動(dòng)，23（5）：42—46.

薄景山，翟慶生，吳兆營(yíng)等，2004.基于土層結(jié)構(gòu)的場(chǎng)地分類方法.地震工程與工程震動(dòng)，24（4）：46—49.

陳國(guó)興，陳繼華，2004.軟弱土層的厚度及埋深對(duì)深厚軟弱場(chǎng)地地震效應(yīng)的影響.世界地震工程，20（3）：66—73.

黃潤(rùn)秋，余嘉順，2003.軟弱夾層的地震動(dòng)強(qiáng)度效應(yīng)研究.自然科學(xué)進(jìn)展，13（11）：1177—1181.

黃玉龍，郭迅，袁一凡等，2000.軟泥夾層對(duì)香港軟土場(chǎng)地地震反應(yīng)的影響.自然災(zāi)害學(xué)報(bào)，9（1）：109—116.

金丹丹，陳國(guó)興，董菲蕃，2012.軟夾層土對(duì)福州盆地地表地震動(dòng)特性的影響.武漢理工大學(xué)學(xué)報(bào)，34（12）：83—88.

盧華喜，梁平英，仇早生，2010.軟弱土層對(duì)中歐場(chǎng)地劃分影響的比較.工業(yè)建筑，40（2）：82—85.

欒茂田，林皋，1992.場(chǎng)地地震反應(yīng)一維非線性計(jì)算模型.工程力學(xué)，9（1）：94—103.

錢勝國(guó)，1994.軟土夾層地基場(chǎng)地土層地震反應(yīng)特性的研究.工程抗震，（1）：32—36.

尚守平，劉方成，王海東，2007.基于阻尼的地震循環(huán)荷載作用下粘土非線性模型.土木工程學(xué)報(bào)，40（3）：74—81.

史丙新，呂悅軍，2010.軟弱場(chǎng)地的地震動(dòng)效應(yīng)研究.見：中國(guó)地震局地殼應(yīng)力研究所，地殼構(gòu)造與地殼應(yīng)力文集（22）.北京：地震出版社，24—33.

謝君斐，石兆吉，1978.天津市震害異常的初步探討.哈爾濱：中國(guó)科學(xué)院工程力學(xué)研究所.

張建毅，薄景山，林瑋等，2013.北川縣城土層特性對(duì)地震動(dòng)的影響.自然災(zāi)害學(xué)報(bào)，22（1）：130—139.

Bardet J.P. and Tobita T.， 2001. NERA- A computer program for nonlinear earthquake site response analyses of layered soil deposits （Manual）. Department of Civil Engineering， University of South California.

Idriss I.M.， Sun J.I.， 1992. Shake91-A computer program for conducting equivalent linear seismic response analysis of horizontally layered soil deposits （User’s Manual）. Department of Civil and Environmental Engineering， University of California， Davis， California.

Seed H.B.， Idriss I.M.， 1970. Soil moduli and damping factors dynamic analysis. Report NO. EERC 70—10，Earthquake Engineering Research Center， University of California， Berkley.

Seed H.B.， Wong R.T.， Idriss I.M.， and Tokimatsu K. ， 1986. Moduli and Damping Factors for Dynamic Analyses of Cohesionless Soils. Journal of Geotechnical Engineering， 112（11）: 1016—1032.

Influence of the Geometric Properties and Shear Wave Velocity of a Soft Soil Layer on the PGA Response of Site

Liu Fangcheng， Wu Jingfang， Chen Bin and Zhao Chengqi

（College of Civil Engineering， Hunan University of Technology， Zhuzhou 412008， China）

The earthquake response of soil deposits containing a soft layer has been of great interest in geotechnical earthquake engineering practice. This paper introduces a study on the influence of a soft layer on the ground motions of a homogeneous deposit， with the variation of geometric properties and shear wave velocity of the soft layer. Numerical studies are conducted by method of real nonlinear one dimensional site response analysis. Results indicate that the existence of a soft soil layer in the deposit will cause sudden change of the acceleration response in the profile， and the PGA may be increased or decreased based on the geometric properties and shear wave velocity of the soft soil. There exists a critical range as which being exceeded the PGA will be decreased monotonically with the increase of the soft soil’s thickness， buried depth and relative shear wave velocity. On the contrary， as the critical range is not exceeded， the PGA will be enlarged in a trend of decreasing after increasing with the increase of the soft soil’s thickness， buried depth and relative shear wave velocity. Formula describing the critical range and empirical methods estimating effects of a soft soil layer on PGA are proposed based on parametric analysis of a homogeneous site. Reasonability of proposed method is validated by numerical calculation results. Foundations of this study may provide reference for further research and application.

Soft soil layer； Site effects； Soil dynamic properties； Site response analysis； Peak ground acceleration

國(guó)家自然科學(xué)基金項(xiàng)目（51108177）；湖南工業(yè)大學(xué)自科基金項(xiàng)目（2012HZX20、2013HZX09）

2015-03-05

劉方成，男，生于1978年。博士，副教授。主要從事土動(dòng)力學(xué)、場(chǎng)地效應(yīng)及土-結(jié)構(gòu)動(dòng)力相互作用研究。E-mail: 455211101@qq.com。