榮 澤，年廷凱，張 浩，鄭德鳳，付崔偉
(1.大連理工大學 海岸和近海工程國家重點實驗室，遼寧 大連 116024；
2.遼寧師范大學 自然地理與空間信息科學重點實驗室，大連 116029；
3.中鐵第一勘察設計院集團有限公司，陜西 西安 710043)
島礁岸坡是一種特殊的海岸形式，通常由陡峭的礁前斜坡和相對平坦的礁坪組成[1]。當遭遇臺風海況時會出現巨浪作用而發生滲流，在坡體內產生不利的孔隙水壓力或滲透壓力，會引發岸坡失穩，造成吹填體滑塌、上部基礎設施破壞等鏈式災害問題；
而地震荷載會導致坡體內產生動應力和超孔壓，同樣會危及坡體穩定性甚至誘發滑坡。因此開展巨浪和強震等極端海況下島礁岸坡的穩定性分析十分重要。
島礁岸坡穩定性分析，整體上與陸上各類邊坡穩定性分析方法類似，但由于島礁地形和水下地質環境的復雜性，導致海上工程地質勘測工作有些難度，要獲得準確的地層分布和土性參數難度更大，特別是要考慮島礁岸坡的滲流和動力特性，使得其穩定性分析變得更加困難，故起步較晚。郭見揚[2]對珊瑚島礁的地層、巖石成因和區域穩定性等進行了定性化描述；
詹文歡等[3]通過對珊瑚島礁群的調查研究，總結歸納了場地內的活動斷裂和多發性災害地質問題；
崔永圣等[4-5]基于各類先進的物探方法，進一步調查了珊瑚島礁的特殊巖土工程地質特征，為后續島礁岸坡的穩定性定量分析提供了基礎數據。胡進軍等[6]針對南海島礁場地的吹填土開展了物理力學性質方面的試驗研究；
劉清君等[7]通過模型試驗研究了規則波作用下島礁地形上波浪破碎點相對距離與Irribarren數、礁坪相對水深以及前坡坡度的密切關系，推導了經驗公式；
并進一步開展了島礁地形在拋石護岸條件下穩定性的試驗研究，發現拋石護岸邊坡的穩定性與岸坡寬度有關[1]；
葉劍紅等[8]采用大型物理模型水槽試驗，研究了波浪荷載作用下南海島礁護岸防波堤的穩定性，探究了極端波浪條件下防波堤的位移、越浪量以及地基內的孔隙壓力是否存在累計上升等現象。對于地震作用下島礁地形的動力響應方面，胡進軍等[9]探究了不同脈沖型地震動情況下珊瑚島礁場地的放大倍數，對珊瑚島礁場地地震反應分析進行了初步探索；
陳國興等[10]進一步研究了珊瑚島礁場地的地震反應特征，填補了南海島礁場地地震效應研究的空白，但未考慮海水-島礁體相互作用；
張巍等[11]則建立了海水-島礁-地震耦合的島礁數值模型，考慮流固耦合作用，深入分析了島礁場地的地震反應規律。
綜上所述，前人對島礁岸坡在波浪、地震荷載作用下的動態反應進行了不少研究，且開展了實際島礁岸坡的物理模型試驗，但采用數值模擬手段進行極端海況下島礁岸坡的滲流和動力響應方面的研究還相對較少，也缺乏針對具體案例的深入穩定性分析。
基于此，本文以南海吹填島礁典型護岸邊坡工程為例，采用數值分析方法與極限平衡法相結合的研究手段，對極端海浪條件下島礁岸坡進行有限元滲流穩定性分析，并開展地震荷載作用下的動力穩定性分析，以期豐富島礁岸坡的穩定性分析方法，為同類工程建設和安全運營提供技術參考。
1.1 滲流有限元控制方程
基于達西定律和質量守恒方程，各向異性土中非穩定滲流的微分方程可表示為:
(1)
式中:ρ為水的密度;β為水壓縮模量;α為土壓縮模量;n為土體孔隙率[12]。穩定滲流的微分方程即等式(1)右側為零。
采用Seep/W有限元分析模塊進行求解，獲得島礁岸坡的滲流場、變形場和應力場以及動態反饋。
1.2 動力有限元控制方程
對于地震荷載作用下的島礁岸坡動應力場求解，其動態響應的有限元控制方程為式(2)，其中荷載{F}由體荷載、邊界荷載、集中荷載和地震荷載共同組成，見式(3)。由牛頓第二定律可知，地震動力荷載為節點加速度向量與總體質量矩陣之積，見式(4)。通過平方根法求解動力有限元方程，可得到每一單元高斯點的應力應變結果[13]。
(2)
{F}={Fg}+{Fb}+{Fs}+{Fn}
(3)
(4)
1.3 基于有限元法的邊坡極限平衡穩定性分析
極限平衡法通過靜力平衡原理分析岸坡的受力特點，在滑移體范圍內將土體劃分成若干個豎向的條塊，針對每個小條塊建立力或力矩的平衡方程，經由抗滑力與下滑力之間的相對大小關系來評價其穩定性。在運用極限平衡法進行穩定性分析時，邊坡破壞的滑移面遵循摩爾-庫侖強度準則，其表達式為：
τf=αtanφ+c
(5)
式中:τf為土體的抗剪強度;σ為剪切面上的法向應力；
φ為土的內摩擦角;c為土的黏聚力。當剪應力τ≥τf時，滑移面被“剪壞”，即當某一平面上的剪應力和法向應力滿足式(5)時，該平面處于極限平衡狀態。
對于滲流及動力問題的分析，有限元極限平衡法很好地結合了極限平衡法和有限元應力分析法的特點[14]，通過彈塑性有限元分析，精準把握邊坡整體應力場的影響，進而確定當前應力場中最危險滑動面的位置及其對應的安全系數。滑動面安全系數Fs為滑動面上土體的抗剪強度和實際剪應力的比值：
(6)
式中:l為滑移面滑弧長度；
c、φ分別為土體的黏聚力和內摩擦角；
本文有限元極限平衡法是通過應力變形分析實現的[15]。
2.1 典型島礁岸坡工程地質條件
島礁護岸防波堤是保證礁坪上吹填體穩定的工程構筑物，以我國南海某一吹填體島礁的典型護岸邊坡工程地質剖面為例，南海典型吹填體岸坡剖面各部分具體尺寸參照圖1(a)，各分區材料的物理力學參數見表1。針對極端海況下吹填島礁岸坡無防波堤(圖1(b))和有防波堤(圖1(c))的兩種岸坡形式，分別開展滲流和動力穩定性分析。其中，島礁護岸防波堤由胸墻及S型護坡兩大部分組成，由混凝土澆筑而成[8]。島礁護岸防波堤的地基材料根據南海海域南沙群島珊瑚礁淺地層的地質特征劃分為上下兩層：上層為鈣質砂，下層為礁灰巖(由砂與碎石的混合物等效)。由于現場島礁礁坪坡面地形變化顯著，因此在建立模型時采用統一坡度近似表示[8]。
圖1 南海典型島礁護岸邊坡剖面圖[8](單位：m)
表1 島礁護岸邊坡各分區材料參數取值
根據邊坡地質原型以及工程地質條件，有護岸防波堤的模型分為四部分：①碎石，②砂碎石混合物，③鈣質砂，④預制混凝土；
無護岸防波堤的模型只有三部分，不包含④預制混凝土。滲流計算應當考慮吹填體有無護岸防波堤保護，因此滲流分析按圖1(b)和圖1(c)兩種不同分區情況考慮(吹填體坡度均為1∶2.5)，網格劃分無護岸防波堤的剖面共計1 534個節點，1 370個單元(見圖2(a))，有護岸防波堤的剖面共計1 532個節點，1 364個單元(見圖2(b))，均采用混合單元(三角形&四邊形)模式。
圖2 南海典型島礁護岸邊坡有限元網格劃分圖
滲流建模應當考慮海水漲退潮，尤其是潮位可能出現的極端情況，因此荷載程度分類如下：
(1) 正常海平面位置，即如圖1所示7.1 m處。
(2) 極端海平面位置，與鈣質砂土層上位面相平，即9.5 m。
2.2 考慮退潮過程的島礁岸坡滲流穩定性分析
實際沿海環境，海水總處于漲潮與退潮的動態變化，本節考慮水位從初始極端海平面9.5 m降落至正常海平面7.1 m，這一動態過程對島礁岸坡穩定性的影響。假定整個退潮過程在12 h內完成，取水位下降結束瞬間進行對比，滑移體及總水頭如圖3所示。按動態水位模擬，無防波堤的情況下最小安全系數由4.506降至2.792(見圖3(a))，有防波堤的情況下由4.959降至3.287(見圖3(b))。瞬態水位降落，引發順坡滲流，使得斜坡穩定性有所降低，而對于有防波堤的岸坡，其安全系數均高于無防波堤的岸坡，這說明了護岸防波堤保證了在大潮位變化情形下島礁岸坡的穩定性。
圖3 水位驟降時島礁岸坡臨界滑移面位置及孔壓分布
2.3 極端波浪荷載下島礁岸坡穩定性分析
風暴潮產生的波浪荷載是南海島礁岸坡所受到的最頻繁動荷載，砂土質邊坡失穩更是與波浪作用密切相關。波浪荷載導致島礁岸坡失穩的主要原因為：(1)土體的剪應力由于波浪作用下發生突變；
(2)土體的抗剪強度由于有效應力下降而減小；
(3)土體內部孔壓急劇變化。隨著剪應力突增和土體抗剪強度下降，滑坡更易失穩[16]。本節將簡諧波浪荷載轉化為等效均布荷載，利用擬靜力方法，探究波浪荷載對島礁岸坡穩定性的影響。
島礁岸坡屬于淺水區，應用艾里波理論定義等效海底波壓力，波浪對海底波壓力的表達式為[17]：
(7)
式中:γw為水重度，取10 kN/m3;H為波高;h為水深。根據前人的研究可得，島礁岸坡防波堤的胸墻受到設防極端風浪沖擊時，最大沖擊力能夠達到50 kPa，護坡所受沖擊力峰值在10 kPa～15 kPa之間[8]。同時，水深越深，波浪作用越弱。對波浪荷載按最不利因素進行擬靜力等效，即認為波浪力均布于坡面，在宏觀上對斜坡表現為垂直坡面向外的拉力，如圖4所示。
圖4 波浪作用下島礁岸坡臨界滑移面及孔壓分布
通過計算可得，無防波堤時在極端水位和退潮至正常水位條件下受波浪作用的岸坡安全系數分別為0.917和0.431，此時岸坡處于失穩狀態；
相對地，有防波堤時在極端水位下最小安全系數為1.784，在退潮后正常水位下為1.654，岸坡處于穩定狀態。
對比3.2節考慮退潮過程的穩定性分析，可知當波浪力對岸坡作用時，有防波堤和無防波堤的岸坡最小安全系數分別降低了約50%、85%；
相對地，防波堤使得極端水位和退潮至正常水位條件下的岸坡最小安全系數相比受波浪作用的無護坡情形分別提高了約49%、74%。以上結果表明風浪沖擊進一步降低了岸坡的穩定性，對岸坡安全穩定產生更大威脅，同時護岸防波堤在應對風浪沖擊時，弱化了波浪沖擊影響，顯著提高了岸坡的安全系數，保證了島礁岸坡的穩定性和耐久性。
3.1 島礁岸坡地震動力響應分析
地震產生的加速度增加了海底斜坡的下滑力，同時土體內部孔隙水壓力增大，形成流體超壓，從而降低斜坡土體的抗剪強度，致使海底斜坡發生失穩破壞[18]。南海吹填島礁場地地震頻發，因此進行島礁岸坡在地震荷載作用下的動力響應分析具有重要意義[19]。本節探究島礁岸坡在20 s實際水平地震加速度作用下(如圖5所示唐山地震區地震波數據)作用下的動力響應，以0.01 s為時間步長，岸坡坡頂前端頂點為位移響應點，觀察岸坡穩定性及變形情況。
圖5 地震時程曲線(20 s)
由地震時程曲線可知，某些振動時刻(加速度為負值，積分位移逆滑移方向增加)可以抑制滑移體滑動，某些振動時刻(加速度為正值，積分位移順滑移方向增加)則促使滑移體滑動。因此，地震動力作用在某一時刻安全系數可能降低，也可能升高。對超出屈服加速度的加速度量進行積分運算，可得出地震動力作用下的岸坡累積變形[20]，位移響應特征點結果如圖6所示。無護岸防波堤時，該砂質邊坡無法保持穩定，變形不斷積累，最終剪切破壞(見圖6(a))；
有護岸防波堤時，邊坡累積變形小于0.1 m(見圖6(b))，相比前者形變較小，穩定性在一定程度上得到保證。
圖6 岸坡累積變形示意圖
3.2 島礁岸坡擬靜力地震穩定性分析
擬靜力法也稱為等效荷載法：加速度可以產生慣性力，引發地震動力效應，這些力以水平力Fh和豎向力Fv的形式作用在每一個條塊的質心上，即：
(8)
(9)
式中：ah和av為水平和豎向加速度;g為重力加速度;W為條塊重量。慣性效應由上述兩個無量綱常數kh和kv來表征。
其中，豎向地震系數kv對邊坡穩定性影響很小，相較而言，水平地震系數kh對島礁岸坡穩定性有很大影響，即便是較小的水平系數也能大大降低島礁岸坡的穩定性，甚至導致其失穩[21]。因此在穩定水面的條件下，控制豎向系數kv≡0，水平地震系數kh從零逐漸遞增，探究安全系數關于水平地震系數變化的敏感性，計算所得結果列于表2，島礁岸坡安全系數的變化如圖7所示。可以看出，隨著水平地震系數增大，最小安全系數平滑地逐漸減小，發生滑坡的風險逐漸升高。對比無護岸防波堤的岸坡可知，防波堤有效提高了各級地震荷載作用下岸坡的安全系數，其作為一種風險儲備而存在，有利于島礁工程的長期安全穩定運營。
表2 最小安全系數統計表
圖7 安全系數Fs關于水平系數kh的變化圖
3.3 島礁岸坡地震動力穩定性分析
擬靜力分析無法表達安全系數在地震活動這一動態過程中的變化，實際上在某個時刻安全系數可能已經低于臨界值。因此，進行島礁岸坡的地震動力穩定性分析是對擬靜力穩定性分析結果的補充和對照；
同時，在3.1節，通過累積位移初步判斷了島礁岸坡兩種剖面形式的穩定性，直觀結論是否正確，應經由動力穩定性分析予以驗證。在考慮地震影響時，由于垂直作用對穩定性的影響較小，往往忽略垂直地震作用對海底邊坡穩定性的影響[22]，基于此，本文動力分析只施加水平方向的地震動力作用。
在地震動力作用下，臨界(或任一非臨界)滑動面在每個時間步上的安全系數計算結果如圖8所示。當遇到地震時，無護岸防波堤的吹填島礁岸坡極易發生滑坡，安全系數降低了約60%。故在南海地震多發海域，護岸防波堤具有十分重要的工程意義，保障了島礁岸坡的安全穩定。
圖8 安全系數Fs隨時間變化圖(10 s)
另外，由圖8可得，地震過程中安全系數變化波動很大，動態安全系數難以全面判斷島礁岸坡是否穩定。因此，要針對各單元進行應力狀態計算，評價岸坡穩定性。定義新的Slope/W模塊分析類型為“地震應力”，以Quake/W模塊為上級目錄，計算可得到滑移體形態如圖9所示，進一步分析島礁岸坡在地震作用下的穩定性(區別于擬靜力法)。無護岸的島礁岸坡最小安全系數為0.419(見圖9(a))，有護岸的島礁岸坡則為1.677(見圖9(b))；
安全系數提高了約75%，與擬靜力法所得結論一致，即護岸防波堤可以顯著提升島礁工程的穩定性。
圖9 地震作用下島礁岸坡總水頭及穩定性示意圖
針對南海典型島礁岸坡，考慮穩態滲流及非穩態退潮情況下的滲流、波浪力和地震荷載等不同海況開展了島礁岸坡的穩定性分析，得出以下結論：
(1) 基于海水退潮這一非穩態滲流的水動力過程，島礁岸坡穩定性會有明顯的降低，上部吹填體整體性受到沖擊，而護岸防波堤可以有效降低滲流影響，防止吹填體在滲流作用下發生滑塌破壞。
(2) 風暴潮對島礁岸坡穩定性威脅較大，建設島礁防波堤有利于降低風暴潮對岸坡的威脅。通過擬靜力法等效波浪荷載，對某一極端全表面拉力狀態進行穩定性分析，發現無護岸防波堤的島礁岸坡在極端水位和退潮后正常水位下均無法承受設防波浪的沖擊，處于極不穩定狀態。
(3) 強震作用對島礁岸坡穩定性影響很大。通過岸坡隨時程變化的震動，對加速度積分所得累積變形可知，未防護的吹填體在強震作用下位移顯著，喪失穩定性。擬靜力和動力穩定性分析進一步說明，累積變形較大的無護岸防波堤處于失穩狀態，護岸防波堤可大幅提升島礁岸坡在各類極端海況條件下的穩定性。
猜你喜歡防波堤護岸島礁淺談水利工程中堤防護岸工程施工技術應用建材發展導向(2022年10期)2022-07-28多孔結構防波堤消浪特性實驗研究港工技術(2022年2期)2022-05-12預制裝配式波浪樁生態護岸研究中國水運(2022年4期)2022-04-27高強度塑鋼板樁在河道護岸工程中的應用研究房地產導刊(2022年4期)2022-04-19生態護岸技術在深圳防洪治理中的應用黑龍江水利科技(2020年8期)2021-01-21寬肩臺防波堤穩定性數值模擬方法研究海洋通報(2020年6期)2020-03-19關于浮式防波堤消能效果及透射系數的研究建材發展導向(2019年10期)2019-08-24體系作戰條件下島礁作戰中輔助決策問題研究航天電子對抗(2019年4期)2019-06-02島礁多能源軍事供電系統容量優化配置研究軍事運籌與系統工程(2018年1期)2018-11-10頂升平臺在強涌浪海域深水防波堤地基處理中的應用中國港灣建設(2017年11期)2017-12-19