基于一維波動理論的特厚覆蓋層壩基動力特性研究

余 挺，邵 磊

（中國電建集團成都勘測設計研究院有限公司，四川省成都市 610072）

0 引言

在覆蓋層上建高土石壩，若地基在地震工況下保持穩定或經工程措施處理后能夠保持穩定，則建壩就具有可行性。天然覆蓋層地基地震動力反應分析是判斷其動力穩定性的基礎。

運用傳統地震動一致性輸入方法（剛性邊界+慣性力法）進行大壩動力反應分析時，需考慮壩體近域地基和遠域地基的影響，一般將臨近壩體地基的河床與山體劃分為近域地基，與壩體一起作為計算分析模型。由于外行波能量具有向遠域地基傳播的逸散效應，理論上近域地基取的范圍越大計算精度越高[1]，但存在網格數量龐大、模型復雜，計算工作量不經濟等局限，導致工程實用性差。若近域地基范圍取的偏小，且采用固定邊界，則遠域地基耗散的外行波能量將從邊界反射回壩體和近域地基內，導致大壩地震響應計算結果失真。

在目前的土石壩動力分析中，對與壩體直接建在基巖上或壩基覆蓋層較淺的情況，一般選取基巖面做為模型邊界；對于壩基為深厚覆蓋層或基巖條件較為復雜的情況，則取向下一定深度和兩岸一定寬度區域的覆蓋層或基巖作為大壩計算模型的近域地基，并做固定邊界處理，對壩體施加慣性力，在邊界處輸入地震加速度[2]。這種計算方法就像是將整個計算模型放在剛度很大的振動臺上，整個系統是能量封閉的。計算分析研究表明[3]，固定邊界條件的計算結果在量值上與較合理的計算結果（波動輸入方法）相差高達20%，且與壩體體積和地基變形模量密切相關，已不能很好地適應厚達數百米級覆蓋層上高土石壩的抗震計算要求。因此，研究將建在深厚覆蓋層上的土石壩作為開放系統的波動問題來求解。

一維波動理論是地震工程界較為常用的波動問題求解方法。該該理論假設土體為水平分層，每一水平子層為線黏彈性體，其動剪切模量與阻尼比均為剪應變的函數，通過計算任一子層對地震波輸入層的傳遞函數，來推求土層地表的地震反應，常用于計算覆蓋層等天然地基的動力響應[4-8]。巖性地基條件下，將基巖考慮為具有自由邊界的彈性半空間，假設上覆土層對基巖內的入射波無影響，而基巖面以下的入射波在基巖內被全部吸收，在計算中取由下臥基巖處向上傳播的輸入加速度峰值等于露頭基巖場地表面設計地震加速度的二分之一，則基巖面的加速度即可達到設計值[9-11]。但針對存在天然覆蓋層的非巖性地基，考慮覆蓋層與基巖相互作用以及覆蓋層土體本身的強非線性動力特性，使得在覆蓋層與基巖的接觸面處，基巖內豎直向上的入射波并不等于向下返回基巖內的反射波，即基巖內入射波在覆蓋層土體與基巖的接觸面處不可能等于覆蓋層不存在時的輸入地震波的一半[12]。故在計算地震波輸入層的傳遞函數時，還應考慮覆蓋層與基巖的相互作用。

有研究表明[13]，是否考慮覆蓋層與基巖的相互作用，對于動力反應時程、加速度峰值及反應譜的影響均不可忽視。若不考慮覆蓋層與基巖與的相互作用，計算得到的動力反應偏小，使得抗震設計趨于不安全。故此，本文根據文獻[13]考慮基巖與上覆土層的相互作用的方法，計算了500m級特厚覆蓋層壩基系統的一維土層地震反應。

1 計算方法與參數

1.1 基本理論[13]

在一維波動理論中，將自由場地地震運動過程簡化為垂直向上傳播的一維剪切波，產生的水平位移可以表示為：

式中：ρ——質量密度；

G——剪切模量；

η——黏滯阻尼系數。

黏滯阻尼系數η可由粘滯阻尼比β表示為：η2Gβ/ω=。為第j層土中Zj處的水平位移。

在頻域中，可以將地震波化為若干簡諧波，其中一條的運動方程可表示為：

將式（2）代入式（1），可得：

式中：

式中：kj——第j層的復波數；

vj——復剪切波速；

——復剪切模量。

G和β是頻率的函數，工程上通常取為常數，則可定義復剪切模量：

式中：jζ——臨界阻尼比。

由式（3）可以看出，介質中某一點的運動Uj由兩部分組成。第一部分是向上傳播的波動，表示的是向下傳播的波動。其中Ej和Fj分別表示振動幅度的大小，稱為波幅系數。

在相鄰界面上的剪應力和位移連續條件下，可得如下遞推關系：

式中：hj——第j層土的厚度；

aj——第j層土的復阻抗比，與頻率無關。

如果從第j層的頂面輸入位移幅值為的地震動，第m層頂面的位移反應為則第m層對第j層的位移傳遞函數為：

當地震動從基巖頂面輸入時，假定基巖以上土層對基巖內的入射波無影響，基巖面以下的反射波在基巖內全部吸收。即基巖內入射波和反射波幅值之間，在基巖頂面處存在關系：（帶′表示不考慮基巖和上覆土層相互作用），從而第m層對于基巖的位移傳遞函數為：

上述假定是建立在把基巖以上全部土層去掉，基巖面變成半無限空間的自由邊界，基巖內的入射波到達自由邊界時被完全反射。現行工程上通常采用的一維波動理論即是采用該簡化原理進行計算。

實際上由于上覆土層與基巖的相互作用，土層與基巖的入射波和反射波形成了復雜的系統，在基覆界面處，基巖內的反射波和入射波并不相等。故此，有學者[13]給出了傳遞函數公式為：

1.2 大壩與河床覆蓋層參數

規劃中的某水電站擋水建筑物采用土質心墻堆石壩，建基面最低高程為2930m，正常蓄水位3070m，死水位3015m，壩頂高程為3080m，最大壩高150m。壩址區河床覆蓋層最大厚度達550m，按其物質組成、層次結構、成因、形成時代和分布情況等，自上而下（由新至老）可劃為6層，其中，③～⑤層為軟弱土層，厚度約240m。第①層由塊碎石土及含砂的塊碎礫石土組成，厚約130m；第②層由含（塊）碎（卵）礫石砂層，粒徑大小懸殊，分選性差，級配不良，粗顆粒未形成骨架，厚約160m；第③層以含礫中粗砂、中細砂為主，分布連續；第④層為粉質黏土層；第⑤層以含礫粗砂、中砂、細砂為主，分布連續，厚度較大；第⑥層以砂卵礫石層為主，局部分布有灰黃色粉土、細砂層，厚約20m。各土層剪切波速、動剪切模量、動剪切模量和阻尼比等參數見表1。

表1 計算參數表Table1 Table of calculation parameters

1.3 地震等級

中國地震局地球物理研究所根據相關標準[14]與壩址區工程地質條件，得到壩址區50年超越概率10%地震動峰值加速度為0.20～0.30g，對應地震基本烈度為Ⅷ度。設計加速度代表值的概率水準，對壅水建筑物應取基準期100年內超越概率2%。根據工程場地所處的區域工作背景與地震地質環境，結合工程經驗判斷，此設計概率下的場地基巖地震動參數水平一般在0.5g以上。

1.4 輸入地震波

在場址區附近對工程場地地震危險性貢獻最大的主要潛源區見表2所示。

表2 主要潛源對壩址峰值加速度貢獻率近似值表 （50年超越概率10%）Table 2 An approximate table of the contribution of main potential sources to the peak acceleration of the dam site

通過地震危險性概率分析，壩址50年超越概率10%的基巖水平峰值加速度為258.8gal，100年超越概率2%的基巖水平峰值加速度為528.9gal，100年超越概率1%的基巖水平峰值加速度為624.0gal。

輸入地震波采用100年超越概率2%的場地譜人工地震波，水平向加速度峰值為0.55g，豎向加速度峰值為水平向的2/3，場地譜人工地震波如圖1所示。

圖1 場地譜順河向水平加速度時程Figure 1 The site spectrum along the river horizontal acceleration time history

2 計算結果分析

2.1 加速度放大倍數

在場地譜地震波輸入作用下加速度放大倍數沿覆蓋層深度分布如圖2所示。可見在場地譜地震波作用下，覆蓋層放大倍數隨高程呈現出先減小后增大的趨勢。考慮和不考慮基巖-上覆覆蓋層相互作用條件下，頂部自由表面放大倍數分別達到0.625和0.536，相差約16%。在軟弱土層為主的③～⑤層，地震動出現較大的衰減。

為了對比分析考慮基巖-上覆覆蓋層相互作用的影響，根據式（14）和式（15）分別計算了具自由邊界基巖工況與基巖-土層相互作用工況下的位移傳遞函數。圖3顯示了兩種情形的位移傳遞函數的幅值比，可以看出前一種情形的傳遞函數之幅值，總是小于（不超過）后一種情況下傳遞函數的幅值。相對應的，覆蓋層頂部的峰值加速度小于考慮到基巖與上覆土層存在相互作用時的結果。

圖2 加速度放大倍數沿覆蓋層深度的分布圖Figure 2 The distribution of acceleration magnification along the depth of the overburden

圖3 具自由邊界基巖工況與基巖-土層相互作用工況之傳遞函數比Figure 3 Transfer function ratio between free boundary bedrock and bedrock-soil interaction

2.2 考慮覆蓋層天然地基上新加土層

為了考慮天然覆蓋層地基上筑壩的情況，在原有概化一維土層模型的基礎上，新加一層50m厚的堆石材料。新加層參數如下：剪切波速根據文獻[15]表7.2建議的土石壩平均剪切波速Vs=300m/s，最大動剪切模量Gdmax=0.89GPa，密度=2180 kN/m3，厚度=50m，λmax=0.254，k1=26。計算結果與天然覆蓋層地基對比見圖4。

圖4 加速度放大倍數沿覆蓋層深度的分布Figure 4 The distribution of acceleration magnification along the depth of the overburden

考慮基巖-上覆覆蓋層相互作用的一維波動理論計算結果表明，地震波從基巖輸入，在覆蓋層中的傳播模式為土層中減小，至頂部自由面增大，且頂部自由表面的放大倍數小于1。新加層情況下（550m原有天然覆蓋層地基+50m新加上覆土層），相當于增加了覆蓋層厚度，頂部加速度放大倍數由天然覆蓋層地基（550m）的0.625減小為0.536，減少了14%，可見隨著高度增加（針對高度500m以上情況），建基面地震動的放大效應呈現逐漸減小的趨勢。取同一水平截面的建基面（550m高度處）對比來看，新加層情況下，放大倍數由0.625減少為0.418，減少約33%，說明上覆結構對建基面影響較大。這一結論不一定適用于覆蓋層上筑壩，因為壩體結構和覆蓋層相比很小，類似于一根筷子，地震波順著壩體傳播上去，壩頂會產生鞭梢效應。另外，以上兩種情況，軟弱土層③～⑤的濾波效應均較為明顯。

2.3 覆蓋層自由表面加速度

圖5示出在基覆界面處輸入場地譜人工波時，覆蓋層頂部表面的地震響應，其對應的加速度反應譜示于圖6。

圖5 覆蓋層頂部自由面加速度時程Figure 5 Time history of overburden surface acceleration

圖6 覆蓋層頂部自由面加速度反應譜Figure 6 Acceleration response spectrum of overburden surface

可見，通過對比計算是否考慮基巖-覆蓋層相互作用兩種工況算出的覆蓋層頂部表面的動力響應，無論加速度峰值還是其譜特征均有差異。與圖4對應的位移傳遞函數相比較，在不考慮基巖-覆蓋層相互作用工況下，覆蓋層頂部表面的動力反應的峰值加速度和反應譜峰值均偏小。

3 結論

本文采用考慮基巖-覆蓋層相互作用的一維波動理論，對強震區特厚覆蓋層壩基進行動力反應分析，得到以下結論：

（1）峰值加速度為0.55g的地震波從基巖輸入，在含軟弱土層的500m級特厚覆蓋層中的傳播模式為入射后先衰減，至頂部自由面增大，且頂部自由表面的放大倍數小于1。傳播過程中覆蓋層所含的軟弱土層濾波效應明顯。

（2）考慮和不考慮基巖-覆蓋層相互作用，頂部自由表面放大倍數分別達到0.625和0.536，相差約16%，后者縮減了土層的放大效應。

（3）上覆結構對建基面影響較大，隨著建筑物高度的增加，覆蓋層建基面地震動的放大倍數呈現逐漸減小的趨勢。