包彩虹，王俊杰

（中交三航（上海）新能源工程有限公司，上海 200137）

0 引言

近年來(lái)，海上風(fēng)能作為一種新的可再生能源，受到了大多數(shù)國(guó)家的青睞?！讹L(fēng)電發(fā)展“十三五”規(guī)劃》中明確指出：到2020年底，風(fēng)電累計(jì)并網(wǎng)裝機(jī)容量確保達(dá)到2.1億kW以上，其中海上風(fēng)電并網(wǎng)裝機(jī)容量達(dá)到500萬(wàn)kW以上。因此，海上風(fēng)電場(chǎng)的建設(shè)成為國(guó)家能源結(jié)構(gòu)轉(zhuǎn)型過(guò)程中必不或缺的一部分。而自升式風(fēng)電安裝船作為海上風(fēng)電施工的專用船舶，也越來(lái)越多受到人們的關(guān)注。

自升式風(fēng)電安裝船主要通過(guò)將樁腿插入海底持力層，頂升船體出水形成一個(gè)在海平面以上的穩(wěn)固安裝平臺(tái)來(lái)進(jìn)行作業(yè)。與石油鉆井平臺(tái)不同，風(fēng)電安裝船具有“快插、快拔”的特性，而樁靴基礎(chǔ)在不同土質(zhì)中的承載能力也具有很大差異。因此，對(duì)不同性質(zhì)的海底土承載力的計(jì)算以及平臺(tái)基礎(chǔ)入泥深度的預(yù)測(cè)是保障海上風(fēng)電安全施工的重要前提，對(duì)風(fēng)電安裝船“站得住、拔的出”具有重要的指導(dǎo)意義。

自升式風(fēng)電安裝船插樁深度預(yù)測(cè)的關(guān)鍵在于地基土極限承載力的計(jì)算，一般把自升式風(fēng)電安裝船樁腿的最大預(yù)壓載荷等于或小于地基極限承載力時(shí)所達(dá)到的入泥深度定義為插樁深度[1]，通過(guò)選取幾個(gè)土層層位（文中所選層位為3～4層，根據(jù)土質(zhì)特征進(jìn)行判斷選?。?jì)算每個(gè)土層底面的地基承載力，找到樁腿最大預(yù)壓載荷等于地基承載力時(shí)所在的土層深度來(lái)獲得。文中主要采用《SNAME》[2]和《海洋井場(chǎng)調(diào)查規(guī)范》修正公式[1]對(duì)地基土極限承載力進(jìn)行計(jì)算。

1 《海洋井場(chǎng)調(diào)查規(guī)范》修正公式

目前國(guó)內(nèi)對(duì)于自升式船插樁深度的計(jì)算主要參考了SY/T 6707—2016《海洋井場(chǎng)調(diào)查規(guī)范》[3]的計(jì)算方法，通過(guò)假設(shè)插樁時(shí)樁靴基礎(chǔ)上部土地完全回填（砂土）或完全不回填（黏性土）來(lái)計(jì)算該土層的地基極限承載力[1]。

1）對(duì)于不排水條件的黏性土，內(nèi)摩擦角φ=0°，樁靴極限承載力公式為：

式中：qu為地基極限荷載；cu為土體強(qiáng)度指標(biāo)；Nc為承載能力系數(shù)；γ1為樁靴基線以上土的有效重度；D為樁靴頂面沒(méi)有覆土回填時(shí)取海底泥面到樁靴基線的距離，樁靴頂面被覆土回填時(shí)取樁靴平均厚度h。

2）對(duì)于砂質(zhì)土，由于極限垂向承載力大[4]，因此考慮樁靴頂面被完全覆土回填，主要計(jì)算公式見(jiàn)為：

式中：qu為地基極限荷載；Nq，Nγ為承載能力系數(shù)，根據(jù)土體內(nèi)摩擦角按太沙基極限承載力系數(shù)表確定；γ1為樁靴基線以上土的有效重度；γ2為樁靴基線以下深度B內(nèi)土的有效重度；B為樁靴寬度；Sq，Sγ為樁靴底面形狀系數(shù)，對(duì)于長(zhǎng)方形樁靴Sq取1，Sγ取0.8。

3）對(duì)于粉土，考慮砂土和黏粒土的成分占比，一般通過(guò)土工試驗(yàn)進(jìn)行確定，表現(xiàn)為排水性時(shí)，按砂土承載力公式進(jìn)行計(jì)算，表現(xiàn)為不排水性時(shí)，按黏土承載力公式進(jìn)行計(jì)算。

2 《SNAME》修正公式

《SNAME》相關(guān)規(guī)范通過(guò)將傳統(tǒng)樁靴基礎(chǔ)等效成圓盤模型[2]，并進(jìn)一步考慮樁靴基礎(chǔ)的上覆土層壓力，從而計(jì)算基礎(chǔ)在不排水黏性土和排水性砂質(zhì)土中的垂向承載力。等效模型見(jiàn)圖1。

圖1 等效模型Fig.1 The equivalent model

另外，本規(guī)范還對(duì)影響樁靴基礎(chǔ)承載能力因素進(jìn)行了分析，主要包括極限樁坑深度和深度系數(shù)，《SNAME》中采用穩(wěn)定系數(shù)法來(lái)評(píng)估樁靴基礎(chǔ)在黏性土中插樁的臨界樁坑深度，但由于沒(méi)有考慮土體回流效應(yīng)，因此文中采用Hossain提出的極限樁坑深度計(jì)算方法[5]，其計(jì)算公式見(jiàn)式（3），深度系數(shù)計(jì)算方法則在后面不同土質(zhì)中分別討論。

式中：H為極限樁坑深度；B為樁靴等效橫截面直徑；Su為土體不排水抗剪強(qiáng)度；γ′為土體的有效容重。

1）對(duì)于不排水條件的黏性土，其垂向地基極限承載力為：式中：Fv為樁靴基礎(chǔ)垂向地基極限承載力；cu為泥面D+B/4深度內(nèi)土體不排水抗剪強(qiáng)度；Nc為承載能力系數(shù)，取5.14；Sc為樁靴形狀系數(shù)；dc為垂向承載力深度系數(shù)；p0′為有效上覆壓力；A為有效橫截面面積。

在不排水條件的黏性土中，《SNAME》規(guī)范通過(guò)結(jié)合Meyerhof[6]和Brinch Hansen[7]的研究成果，以D/B=1為分界點(diǎn)計(jì)算垂向承載力深度系數(shù)，但戴兵等人通過(guò)相關(guān)實(shí)驗(yàn)發(fā)現(xiàn)現(xiàn)有的深度系數(shù)低估了樁靴基礎(chǔ)的垂向承載力[8]，因此文中采用Gourvenec提出的深度系數(shù)計(jì)算方法[9]，其公式為：

式中：D為樁靴基礎(chǔ)入泥深度；B為樁靴等效橫截面直徑。

2）對(duì)于排水條件的砂土，其垂向地基極限承載力為：

式中：Fv為樁靴基礎(chǔ)垂向地基極限承載力；γ′為土體的有效容重；B為樁靴等效橫截面直徑；Nγ，Nq為承載能力系數(shù)；Sγ，Sq為樁靴底面形狀系數(shù)；dγ，dq為垂向承載力深度系數(shù)；p0′為有效上覆壓力；A為有效橫截面面積。

其中垂向承載力深度系數(shù)dγ取值1，dq計(jì)算公式為：

式中：φ為土體內(nèi)摩擦角。

由于《SNAME》規(guī)范只考慮了帶樁靴基礎(chǔ)的垂直極限承載力，而忽略了樁靴上部土體回填和樁靴基礎(chǔ)排開(kāi)土的有效重力，因此文中在進(jìn)行不同土質(zhì)地基極限承載力計(jì)算時(shí)還需考慮這兩部分的受力情況，從而得到自升式風(fēng)電安裝船在各個(gè)土層的插樁深度，其計(jì)算公式為：

式中：Q為樁靴基礎(chǔ)總的極限承載力；γD為從泥面到插樁深度的土體平均有效容重；γH為從泥面到極限樁坑深度的土體平均有效容重；γ′為土體的有效容重；V′為樁靴下部排開(kāi)土的體積；V為樁靴基礎(chǔ)體積。

3 實(shí)例分析

通過(guò)選取華能如東300 MW海上風(fēng)電場(chǎng)和江蘇響水200 MW海上風(fēng)電場(chǎng)項(xiàng)目，分別利用文中所介紹的兩種公式對(duì)這兩個(gè)風(fēng)場(chǎng)區(qū)域內(nèi)相關(guān)機(jī)位進(jìn)行自升式風(fēng)電安裝船插樁深度的計(jì)算，與實(shí)際插深進(jìn)行比對(duì)，從而分析這兩種公式在不同土質(zhì)中的適用性。

自升式風(fēng)電船—“三航風(fēng)華”號(hào)主要參數(shù)見(jiàn)表1，樁靴基礎(chǔ)見(jiàn)圖2。

表1 “三航風(fēng)華”號(hào)主要船型參數(shù)Table1 Main parametersof"San Hang Feng Hua"ship

圖2 樁靴基礎(chǔ)圖Fig.2 The drawing of pile-shoe foundation

3.1 華能如東300 MW海上風(fēng)電場(chǎng)項(xiàng)目

海上風(fēng)電場(chǎng)項(xiàng)目位于江蘇省如東縣的爛沙海域，海底地形變化較大，風(fēng)電場(chǎng)形狀呈不規(guī)則四邊形，風(fēng)場(chǎng)中心離岸約23 km。風(fēng)電場(chǎng)擬安裝38臺(tái)上海電氣4.0 MW、12臺(tái)遠(yuǎn)景能源4.2 MW以及20臺(tái)重慶海裝5.0 MW風(fēng)電機(jī)組。其中37～50號(hào)風(fēng)機(jī)采用“三航風(fēng)華”號(hào)進(jìn)行安裝。

現(xiàn)選取48號(hào)機(jī)位作為典型機(jī)位，48號(hào)機(jī)位土層相關(guān)參數(shù)見(jiàn)表2。分別采用《海洋井場(chǎng)調(diào)查規(guī)范》和《SNAME》修正公式，計(jì)算“三航風(fēng)華”號(hào)插樁深度，具體數(shù)據(jù)見(jiàn)圖3。

表2 48號(hào)機(jī)位土層參數(shù)Table2 The soil parametersof No.48 wind turbine position

圖3 48號(hào)機(jī)位插樁深度對(duì)比Fig.3 Comparison of the plug-in depth of pile in No.48 wind turbine position

由圖3可知，采用《SNAME》修正公式計(jì)算出的地基承載力為5 500 t時(shí)，樁靴入泥深度為1.52 m；《海洋井場(chǎng)調(diào)查規(guī)范》修正公式計(jì)算出的地基承載力為5 500 t時(shí)，插樁深度為1.7 m，“三航風(fēng)華”號(hào)實(shí)際插深為1.5 m。對(duì)比圖3（a）、（b）可知，從第二層土層（2 m后）開(kāi)始，入泥深度曲線發(fā)生明顯差異，經(jīng)分析發(fā)現(xiàn)《海洋井場(chǎng)調(diào)查規(guī)范》在砂質(zhì)土的計(jì)算中，考慮砂質(zhì)土垂向承載力極大，樁靴頂面被完全覆土回填，公式中D取值為樁靴平均厚度，但本次計(jì)算了3個(gè)砂性土層的地基承載力，而自升式平臺(tái)樁腿的最大預(yù)壓載荷為5 500 t，因此運(yùn)用《海洋井場(chǎng)調(diào)查規(guī)范》所計(jì)算出的第②-③層的入泥深度曲線不具備實(shí)際應(yīng)用價(jià)值。為進(jìn)一步研究，本文又選取了4個(gè)機(jī)位分別對(duì)4個(gè)樁腿的插樁深度進(jìn)行計(jì)算，其插樁深度計(jì)算結(jié)果見(jiàn)表3。

考慮到自升式風(fēng)電安裝船4個(gè)支腿預(yù)壓載荷各不相同，且所處地質(zhì)環(huán)境也可能存在明顯差異，因此根據(jù)“同一機(jī)位，隨著預(yù)壓載的減小，其插樁深度也愈小”的原理，對(duì)存在明顯誤差的數(shù)據(jù)進(jìn)行剔除，進(jìn)而對(duì)計(jì)算結(jié)果進(jìn)行相關(guān)性分析，相關(guān)性計(jì)算公式為（9）-（10），所得計(jì)算結(jié)果見(jiàn)圖4。

表3 兩種公式插樁深度計(jì)算Table 3 Calculation of the plug-in depth of pile about two formulas

式中：x為理論插樁深度；y為實(shí)際插樁深度；a為修正系數(shù)；b為表相關(guān)性系數(shù)；當(dāng)b越接近1且a越接近0時(shí)，x對(duì)y相關(guān)性越好，即理論插樁深度越接近實(shí)際插深。

圖4 如東風(fēng)場(chǎng)兩種公式計(jì)算值與實(shí)際插深相關(guān)性對(duì)比Fig.4 Thecorrelation contrast of calculated valuesand actual plug-in depth of pilewith two formulasin RUDONGwind farm

由計(jì)算結(jié)果可知：《SNAME》修正公式a值等于0.21，b值等于0.92，而《海洋井場(chǎng)調(diào)查規(guī)范》修正公式中a值為0.16，b值為0.83，則認(rèn)為在排水條件的砂質(zhì)土中，《SNAME》修正公式和《海洋井場(chǎng)調(diào)查規(guī)范》修正公式均可適用。

3.2 江蘇響水200 MW海上風(fēng)電場(chǎng)項(xiàng)目

圖5 18號(hào)機(jī)位插樁深度對(duì)比Fig.5 Comparison of the plug-in depth of pile in No.18 wind turbine position

江蘇響水風(fēng)電場(chǎng)項(xiàng)目位于響水縣灌東鹽場(chǎng)，三圩鹽場(chǎng)外側(cè)海域，風(fēng)電場(chǎng)離岸距離10 km，沿海岸線方向長(zhǎng)13.4 km，垂直于海岸線方向?qū)?.6 km，場(chǎng)區(qū)面積34.7 m2，水深8～12 m。本項(xiàng)工程總裝機(jī)容量200 MW，采用西門子4 MW和金風(fēng)3 MW風(fēng)電機(jī)組進(jìn)行安裝。

現(xiàn)選擇18號(hào)機(jī)位作為典型機(jī)位，分別采用《海洋井場(chǎng)調(diào)查規(guī)范》和《SNAME》修正公式進(jìn)行計(jì)算，計(jì)算“三航風(fēng)華”號(hào)插樁深度見(jiàn)圖5。18號(hào)機(jī)位土層相關(guān)參數(shù)見(jiàn)表4。

表4 18號(hào)機(jī)位土層參數(shù)Table4 Thesoil parametersof No.18wind turbineposition

由圖5可知，使用《SNAME》修正公式計(jì)算出的地基承載力為5 500 t時(shí)，樁腿插入泥土深度為11.6 m，《海洋井場(chǎng)調(diào)查規(guī)范》修正公式在地基承載力為5 500 t時(shí)，樁腿入泥深度為18.4 m，“三航風(fēng)華”號(hào)實(shí)際插深為12 m。為進(jìn)一步研究公式的普遍性，本文又在響水風(fēng)電場(chǎng)選取了4個(gè)機(jī)位，并分別對(duì)4個(gè)樁腿的插樁深度進(jìn)行計(jì)算，其插樁深度計(jì)算結(jié)果見(jiàn)表5，通過(guò)選取表5中所計(jì)算出的兩組數(shù)據(jù)，進(jìn)行理論插深與實(shí)際插深的相關(guān)性分分析，得到相關(guān)性曲線見(jiàn)圖6。

表5 兩種公式插樁深度計(jì)算Table5 Calculation of theplug-in depth of pile about two formulas

圖6 響水風(fēng)場(chǎng)兩種公式計(jì)算值與實(shí)際插深相關(guān)性對(duì)比Fig.6 The correlation contrast of calculated values and the actual plug-in depth of pilewith two formulasin XIANGSHUI wind farm

由計(jì)算結(jié)果可知：《SNAME》修正公式a值等于-0.13，b值等于1.02，而《海洋井場(chǎng)調(diào)查規(guī)范》修正公式中a值為4.09，b值為0.5，因此認(rèn)為在不排水條件的黏性土層中，《SNAME》修正公式所計(jì)算出的插樁深度更具有適用性。

4 基于MATLAB建立的插樁深度軟件

為了更為方便、有效地進(jìn)行插樁深度的計(jì)算，本研究利用Matlab的圖形用戶界面（GUI）功能，編譯了一款程序軟件?；谠撥浖M(jìn)行插樁深度計(jì)算，具有速度快、可靠性高等優(yōu)點(diǎn)，很大的縮減了工作量。另外針對(duì)不同船樁腿參數(shù)以及土壤地質(zhì)情況，該軟件可做出相應(yīng)的調(diào)整，具有普適性，進(jìn)一步方便了插樁深度的研究。

5 結(jié)語(yǔ)

文章通過(guò)選取兩種不同規(guī)范的插樁深度修正公式，分別應(yīng)用于“三航風(fēng)華”號(hào)在華能如東和江蘇響水風(fēng)電場(chǎng)施工作業(yè)，計(jì)算結(jié)果表明，在砂性土中，由于地基土承載力較大，砂質(zhì)層基本為樁靴所在持力層，因此兩種公式得到的插樁深度沒(méi)有顯著差異，而在黏性土層中，由于持力層需要到達(dá)較深的砂土層，土體回流影響較大，相關(guān)系數(shù)的選取也具有很大差異，通過(guò)與實(shí)際插深進(jìn)行相關(guān)性分析，《SNAME》修正公式在該性質(zhì)的土層更為適用。

此外，基于MATLAB建立了插樁深度計(jì)算軟件，可適用于不同的風(fēng)電安裝船及工況，對(duì)于以后插樁深度的研究提供了方便。