陳 亮

(華東電力設計院，上海 200331)

深厚軟土地基上新建煤場與海堤的最小安全距離分析

陳 亮

(華東電力設計院，上海 200331)

軟土在我國廣泛分布，其工程性質極差。本文結合實際工程，運用傳統(tǒng)極限平衡法和FLAC3D數值計算方法，綜合分析了深厚軟土地基上新建煤場大面積堆載對鄰近海堤地基穩(wěn)定性的影響，提出了煤場與海堤之間最小的安全距離。

軟土地基；最小安全距離；極限平衡法；FLAC3D。

1 概述

軟土在我國沿海一帶分布較廣，如渤海灣、天津、長江三角洲以及浙江、福建等地的沿海地區(qū)都存在海相沉積的深厚軟土層，普遍具有高含水量、高孔隙比、高壓縮性、高靈敏度和低強度等特性，工程性質極差。

近年來我國沿海軟土地區(qū)建設了數量眾多的燃煤火力發(fā)電廠，為便于運輸，其新建煤場布置往往緊靠廠區(qū)外側的海堤。由于軟土地基的性質較差，當煤場與海堤的間距過小時，煤場的大面積堆載作用可能會對海堤的地基穩(wěn)定性造成明顯的不利影響，而兩者間距過大時，則又會造成電廠建設用地的浪費。針對這種現(xiàn)狀，對軟土地基進行專門的穩(wěn)定性分析，提出新建煤場與海堤之間最小的安全距離是十分必要的。

本文以我國沿海某大型火力發(fā)電廠為例，通過對深厚軟土地基上新建煤場的大面積堆載對海堤地基穩(wěn)定性影響的分析，論述相應的計算方法，供類似工程參考。

2 工程概況

本工程采用斜坡式土石混合堤結構，南北向布置，軸線長887m，按1級堤防工程設計，堤頂高程為6.5 m，堤頂寬度為8.0 m，海堤外坡采用1∶2.5復式斷面，內坡采用1∶5，外海側采用細石混凝土灌砌塊石護面，結構斷面后背采用土體閉氣并用無紡反濾布進行反濾。該海堤已建成，建設時采用了分級堆載+插設塑料排水板加快排水固結的地基處理方法，現(xiàn)處于正常運行期。

電廠新建煤場長約580 m，寬約220 m，南北向布置。設計最大堆煤高度12 m，為保證煤場地基的穩(wěn)定性，擬采用分級堆載，煤堆休止角按46°考慮。海堤及新建煤場的剖面示意圖見圖1。

圖1 海堤、煤場剖面示意圖

3 場地工程地質條件

工程場地所處地貌類型為低山丘陵和海灣灘涂，場地內第四系海相沉積軟土層厚度變化較大，靠近山腳處約0 m～15.0 m，向大海方向逐漸變厚，海堤處軟土厚度普遍大于40 m。

場地內主要地層的分布情況描述如下：

①填土：主要為電廠建設時回填的開山土石，松散～稍密，厚度約3.0 m～6.0 m。

②淤泥：流塑，夾少量腐植質和貝殼碎屑，層厚約10.6 m～28.5 m。

④淤泥質粘土：軟塑，含少量腐植質及貝殼碎片，具片狀層理，層厚約4.4 m～23.0 m。

⑤2碎石土：中密，碎石粒徑4 cm～10 cm，少量大于10 cm，碎石成分為凝灰熔巖，含量50%～60%，碎石間充滿泥質和砂礫，未膠結，層頂埋深約43.5 cm～53.0 m，層厚約3.8 m～17.0 m。

⑥2中等風化凝灰熔巖：堅硬，節(jié)理發(fā)育，巖芯呈短柱和碎塊狀。

根據本工程場地的地震安全性評價報告，在50年超越概率10%情況下，場地地表水平峰值加速度為0.057 g，對應的地震基本烈度為6°。

4 影響海堤地基穩(wěn)定性因素分析

地基土體結構類型是穩(wěn)定性分析的基礎。從本工程的地質條件來看，場地內分布有深厚的軟土層，其滲透性小，抗剪強度較低，工程性質不良，在外部荷載及特定條件下，當土體內部某個面上的剪應力達到或超過了它的抗剪強度時，原有的穩(wěn)定平衡遭到破壞，土體發(fā)生剪切變形，極有可能會引起軟土地基的大規(guī)模整體滑移。根據以往經驗，軟土地基破壞模式以圓弧滑動為主。

顯然，影響本工程海堤地基穩(wěn)定性的內在因素是深厚軟土及其不良工程性質，而外在因素則是煤場的大面積堆載、潮水、地形及地震作用等。

5 海堤穩(wěn)定性計算與分析

5.1 計算工況

由于本工程海堤已處于正常運行期，故不需進行海堤施工期的整體抗滑穩(wěn)定計算。按《海堤工程設計規(guī)范》(SL435-2008)10.2.2條規(guī)定，僅需進行正常運用情況及非常運用情況Ⅱ下的穩(wěn)定計算。另外，煤場的堆載能起到反壓海堤背海坡側坡腳的效果，對海堤內側的穩(wěn)定性是有利的，故計算時只需考慮海堤臨海坡側的穩(wěn)定性即可，即可能滑動面的剪出口位于海堤外側。

具體的計算工況及其水位組合見表1。

表1 計算工況及水位組合

5.2 安全系數

本海堤工程的級別為1級，按《海堤工程設計規(guī)范》(SL435-2008)10.2.3條規(guī)定，采用瑞典圓弧滑動法進行海堤整體抗滑穩(wěn)定計算時，正常運用條件下安全系數不應小于1.30，非常運用條件Ⅱ下安全系數不應小于1.10。

5.3 計算參數

做抗滑穩(wěn)定性分析時，土的抗剪強度指標可采用三軸抗剪強度、直剪指標等，具體應根據海堤的工作狀態(tài)和采用的計算方法選用不同的強度指標。本文中各土層的物理力學計算參數是綜合地勘報告的實測結果和當地相關工程經驗確定的，并考慮到海堤下土體經長期排水固結，物理力學性質有所改善，適當提高了海堤下土體的計算參數，但總體上計算參數仍偏于保守。具體取值見表2。

表2 計算參數

5.4 極限平衡法計算

⑴計算方法和原理

極限平衡法是邊坡穩(wěn)定計算中最廣泛采用的方法，它只需提供少量力學參數就能計算出邊坡設計需要的穩(wěn)定性指標。本文針對地基的可能滑動形式，采用瑞典圓弧滑動法[2]進行計算，其計算簡圖見圖2。

圖2 圓弧滑動法計算簡圖

抗滑穩(wěn)定安全系數K按式1計算。

式中： K為抗滑穩(wěn)定安全系數；W1i、W2i、W＇

2i、W＇3i為第i個土條浸潤線以上的土體的天然重量、浸潤線與外坡水位線之間的土體的飽和重量、浸潤線與外坡水位線之間的土體的浮重量、外坡水位線以下的土體浮重量，kN；αi為第i個土條底面中點的徑向與豎直方向的夾角，(°)；φi、Ci為第i個土條底部土體的總抗剪強度指標，(°、kPa)。

⑵計算結果及分析

綜合分析場地的地層分布和海堤外側的地形情況，本文選取3個最不利的海堤斷面，按新建煤場與海堤的不同間距條件，對海堤臨海坡的穩(wěn)定性進行反復試算。

計算結果表明，當煤堆坡腳距海堤內側堤腳9 m時，在水位降落工況下，3個剖面計算得到的最小安全系數均超過1.50，滿足不小于1.30的要求；而在地震工況下，3個剖面的計算得到的最小安全系數值均略大于1.10，較接近1.10，這說明此時煤場堆載對海堤穩(wěn)定性的影響達到規(guī)范要求的安全臨界狀態(tài)。3個剖面的計算結果見表3，典型1—1′剖面的計算滑面圖見圖3。

表3 穩(wěn)定性計算結果

圖3 1—1′剖面大圓弧滑動計算滑面圖(水位降落工況、地震工況)

從滑面圖可以看到，當煤堆坡腳距海堤內側堤腳9 m時，在水位降落工況下，滑動面剪出口位置在-115 m～-160 m(以海堤臨海坡坡腳作為原點，靠海側取負值，下同)，滑動面剪入口位置在130 m～150 m，即距煤堆坡腳30 m～50 m處；在地震工況下，滑動面剪出口位置在-130 m～-180 m，滑動面剪入口位置在130 m～150 m，即距煤堆坡腳30 m～50 m處。兩種工況下的滑動面弧底埋深約-40 m～-50 m，位于基巖面以上的軟土中。

5.5 FLAC3D數值分析

⑴ 算方法和原理

考慮到本文探討的是巖土大變形問題，而FLAC3D程序在計算過程中允許材料發(fā)生屈服及流變，可以摸擬巖土的力學性能，尤其在彈塑性分析、大變形分析方面有其獨到的優(yōu)點，因此本文采用FLAC3D進行數值模擬，對海堤的變形和穩(wěn)定性做進一步分析。

由于上述1—1′剖面涉及海堤外側水下地形最為復雜，相對2—2′、3—3′剖面而言，發(fā)生整體滑動的可能性更大，故本次FLAC3D分析只針對1—1′剖面進行，重點關注當煤堆坡腳距海堤內側堤腳9 m時，煤場堆載是否會造成超大規(guī)模的圓弧滑動。

⑵計算模型

FLAC3D建模時，若模型嚴格與原始地形匹配，則單元剖分會過于精細，計算耗時太久，而對整體圓弧的穩(wěn)定性計算結果則影響甚微，故本文建模時，適當忽略一些細節(jié)的地形變化。具體的計算模型見圖4，計算模型中土層分組圖見圖5。

圖4 FLAC3D模型圖

圖5 FLAC3D土層分組圖

對圖4、圖5表示的計算模型說明如下：

①模型參數：x方向，海堤外側距離取150 m，海堤寬度取90 m，海堤內側預留9 m(即煤堆坡腳到海堤內側坡腳距離9m)，煤堆寬度取72 m；z方向，模型底部定義為0 m(坐標原點在模型左下角)，模型最左側高48 m，模型局部最低處高46 m，海堤底面標高60 m，海堤頂面標高66.5 m，煤堆底面標高63 m(煤堆堆高12 m，自設計場坪標高3 m起算)，煤堆頂面標高75 m；y方向，取3 m，以保證FLAC3D差分計算的精度。

②網格剖分：x方向3m，z方向2 m，y方向3 m，整個模型共劃分為3131個區(qū)域單元，6520個節(jié)點。

③邊界條件：模型左右側x方向約束，模型底部z=0處x、y、z三個方向約束，垂直紙面y方向約束。

④煤場堆載的模擬：煤場堆載是按分級加載進行的，F(xiàn)LAC3D針對該種情形有兩種處理方式：分層堆載，一層層建實體單元(煤堆)的模型，或者用等量的荷載替代實體單元(煤堆)分批次作用于地基上，本文選用后一種方式模擬煤場的堆載。

⑤土體模型：土體本構模型取摩爾-庫侖彈塑性模型，其物理力學性質指標選取同極限平衡法。

⑶計算結果及分析

FLAC3D中判定潛在滑動面存在的依據，有三類圖：位移矢量圖或位移云圖、剪切應變率云圖、塊狀態(tài)圖。本文將煤場堆載前后的模型位移云圖和剪切應變率云圖對比，以分析煤場堆載對海堤的影響。

煤場堆載前的模型位移云圖和剪切應變率云圖見圖6和圖7。從圖6可以看出，煤場堆煤前模型位移最大值出現(xiàn)在海堤處，其主要原因是地基在海堤自重作用下產生了變形。從圖7可以看出，在屏蔽海底陡降坡處剪切變形較大的區(qū)域后，整個模型顯示出潛在的從海堤背水坡到海底陡降處的圓弧滑動面，但其剪切應變率在1E-9數量級，非常微小，說明煤場堆載前海堤是相當穩(wěn)定的。

圖6 未堆煤前模型位移云圖

圖7 未堆煤前模型剪切應變率云圖

煤場堆載后，由堆煤荷載造成的新增位移云圖見圖8?？梢钥闯觯m然位移云圖顯示出模型存在超大圓弧滑動面的可能，但新增的位移值量級在厘米級，說明當煤堆坡腳距海堤內側堤腳9 m時，海堤整體上仍是穩(wěn)定的。

圖8 堆煤后模型位移云圖

5.6 綜合分析

分析極限平衡法和FLAC3D數值模擬的結果，可以看出，雖然這兩種計算方法反映的滑動面位置略有差異，但極限平衡法與FLAC3D中的潛在破壞面形式非常相似，F(xiàn)LAC3D數值分析結果對由極限平衡法得到的穩(wěn)定性計算結果也做了較好的驗證，因此綜合判定本工程新建煤場與海堤之間最小的安全距離取為9 m是合理的。

6 結語

在充分利用建設用地的基礎上，確保海堤的穩(wěn)定性對于沿海軟土地區(qū)的電力工程而言，是一個極其重要的課題。本文結合實際工程，運用傳統(tǒng)極限平衡法和FLAC3D數值計算方法，分析了深厚軟土地基上新建煤場大面積堆載對鄰近海堤地基穩(wěn)定性的影響，提出了煤場與海堤之間最小的安全距離，對本工程建設起到了較好的指導作用，其計算方法對類似工程也有一定的借鑒意義。由于巖土工程普遍存在一定程度的不確定性，在類似工程的穩(wěn)定性評價中，應盡量綜合運用多種計算方法進行定性和定量分析，以提高結論的可靠性。

[1]SL435-2008，海堤工程設計規(guī)范[S]．

[2]黃求順，張四平，胡岱文．邊坡工程[M]．重慶：重慶大學出版社，2003．

[3]盧延浩．巖土數值分析[M]．北京：中國水利水電出版社，2008．

Stability Analysis of Sea Dike In fl uenced by the Nearby Coal Yard Stowage on Deep Soft Clay Ground

CHEN Liang
(East China Electric Power Design Institute, Shanghai 200331, China)

Soft clay is widely distributed in our country and its engineering properties are very bad. According the project sample, in this paper the traditional limit equilibrium method and FLAC3D are comprehensively used to analyze the stability of the sea dike on deep soft clay ground in fl uenced by the nearby new-bulit widespread coal yard stowage, the minimum safe distance between the coal yard and the sea dike is raised.

soft clay ground; minimum safe distance; limit equilibrium method; FLAC3D.

TU4

B

1671-9913(2011)01-0010-05

2010-09-07

陳亮(1979- )，男，江蘇啟東人，碩士，工程師，注冊土木工程師(巖土)。