摘要： 隨著海洋平臺以及海上風(fēng)電規(guī)模的擴(kuò)大，樁基可打入性分析的精度越來越受到關(guān)注。打樁過程中土阻力的合理評估是提高可打入性分析準(zhǔn)確性的前提，而土阻力分析中最重要的是端阻和側(cè)阻折減系數(shù)的取值。以11個(gè)海上平臺的99個(gè)打樁記錄為數(shù)據(jù)庫，首先，通過打樁記錄數(shù)據(jù)證明弗雷歇距離算法在海上平臺樁基打樁記錄之間相似度計(jì)算的適用性，確定該海域弗雷歇距離的閾值，并對閾值結(jié)果的合理性進(jìn)行驗(yàn)證;其次，通過計(jì)算現(xiàn)場打樁記錄反算的土阻力與不同折減系數(shù)下計(jì)算得到的土阻力之間的弗雷歇距離，對該海域的折減系數(shù)進(jìn)行探討，并給出有土塞和無土塞工況下折減系數(shù)的推薦值。研究結(jié)果可為樁基可打入性分析中的土阻力的計(jì)算提供借鑒。

關(guān)鍵詞： 土阻力; 弗雷歇距離; 樁基可打入性分析; 海上平臺; 大直徑管樁

中圖分類號： TU43 文獻(xiàn)標(biāo)志碼：A 文章編號： 1000-0844（2024）06-1345-11

DOI：10.20000/j.1000-0844.20230522003

Method for calculating the soil resistance during offshore

pile driving based on Fréchet distance

YIN Jiangsong， LI Sa

（School of Civil Engineering， Tianjin University， Tianjin 300072， China）

Abstract： With the development of offshore platforms and offshore wind power， the accuracy of pile foundation drivability analysis， which is greatly impacted by soil resistance to driving （SRD）， has attracted more and more attention. In SRD analysis， the key is determining the reduction coefficients for end-bearing resistance and shaft resistance. Based on a database that includes 99 pile driving records from 11 offshore platforms， the applicability of the Fréchet distance algorithm in evaluating the similarity between pile driving records was verified， and the threshold of Fréchet distance in the sea area was determined and verified. Then， the Fréchet distance between the calculated SRDs under different reduction coefficients and SRD from the back analysis of pile driving records was calculated， and the recommended values of reduction coefficient were proposed in the cases of soil plug and no soil plug. The results of this study could be helpful for the calculation of SRD in pile drivability analysis.

Keywords： soil resistance; Fréchet distance; pile drivability analysis; offshore platform; large-diameter pipe pile

0 引言

樁基的可打入性分析對于樁基設(shè)計(jì)和現(xiàn)場施工十分重要，其分析結(jié)果可以提供錘擊數(shù)、樁身應(yīng)力以及打樁時(shí)間，這三個(gè)參數(shù)能夠幫助工程人員合理選擇樁錘，以便在最短的時(shí)間內(nèi)安裝樁基到指定設(shè)計(jì)深度，同時(shí)避免溜樁和拒錘事故的發(fā)生^［1-4］。

在進(jìn)行樁基可打入性分析時(shí)，合理確定打樁過程中的土阻力是提高樁基可打入性分析準(zhǔn)確度的關(guān)鍵。由于打樁是一個(gè)動態(tài)過程，打樁過程中的土阻力與打樁結(jié)束后的樁基長期承載力有明顯差異。近年來，不少學(xué)者在樁基承載力計(jì)算方法的基礎(chǔ)上提出了不同工況下的土阻力的計(jì)算方法。Semple等^［5］提出在砂土層中依據(jù)《美國石油協(xié)會（American Petroleum Institute，API）規(guī)范》^［6］中的樁基承載力計(jì)算方法來計(jì)算土阻力，在黏土層中根據(jù)應(yīng)力歷史理論基于黏土的超固結(jié)比引入承載力折減因子來計(jì)算土阻力。Stevens等^［7］提出在計(jì)算土阻力時(shí)對于有土塞和無土塞工況應(yīng)分別考慮上限結(jié)果和下限結(jié)果。有土塞工況下，土阻力包括樁外壁阻力以及樁端全面積阻力，下限結(jié)果為依據(jù)API規(guī)范計(jì)算的樁外壁阻力與樁端阻力之和，上限結(jié)果為提高了30%的樁外壁阻力和提高了50%的樁端阻力之和;無土塞工況下，土阻力包括樁內(nèi)壁阻力、樁外壁阻力以及樁端環(huán)面積阻力，無土塞的下限結(jié)果中假定樁內(nèi)壁阻力為外壁阻力的一半，上限結(jié)果則假定內(nèi)壁阻力等于外壁阻力。此外，Stevens還提出了對于低質(zhì)量巖石采用砂土參數(shù)計(jì)算土阻力，對于高質(zhì)量巖石依據(jù)無側(cè)限抗壓強(qiáng)度計(jì)算樁端土阻力，樁側(cè)土阻力依然依據(jù)砂土參數(shù)計(jì)算。

API規(guī)范中的承載力依據(jù)各個(gè)土層的室內(nèi)試驗(yàn)參數(shù)計(jì)算，諸如黏性土的不排水抗剪強(qiáng)度或者砂土的摩擦角。不少學(xué)者也提出了依據(jù)靜力觸探試驗(yàn)（Cone Penetration Test，CPT）參數(shù)計(jì)算土阻力的方法?；诒焙５拇笾睆焦軜冬F(xiàn)場打樁記錄，Alm等^［8］提出了考慮側(cè)阻退化的，依據(jù)CPT試驗(yàn)的錐尖土阻力計(jì)算方法。在計(jì)算無土塞工況下的土阻力時(shí)，Alm建議樁內(nèi)壁阻力和外壁阻力都減小50%，由于大直徑管樁直徑與壁厚比值很大，導(dǎo)致內(nèi)壁阻力和外壁阻力差異很小，因此可以近似為無土塞工況下只考慮外壁阻力，不考慮內(nèi)壁阻力。Prendergast等^［9］在3個(gè)基于CPT的樁基承載力計(jì)算方法（ICP-05^［10］，UWA-05^［11］，F(xiàn)urgo-05^［12］）的基礎(chǔ)上，提出了考慮樁基恢復(fù)和基底殘余應(yīng)力的土阻力計(jì)算方法。

上述成果都是根據(jù)現(xiàn)場的數(shù)據(jù)采用半理論半經(jīng)驗(yàn)的方法獲得的。樁基沉樁過程中土阻力的不同會導(dǎo)致樁基可打入性結(jié)果的不同^［13-14］，同時(shí)打樁過程中的土阻力隨深度變化曲線可以利用現(xiàn)場打樁記錄的反分析結(jié)果獲得。為充分利用現(xiàn)場測試數(shù)據(jù)，本文提出采用曲線相似度計(jì)算方法，獲得打樁過程中土阻力與樁基承載力之間的關(guān)系。

曲線相似度計(jì)算方法分為等長曲線相似度計(jì)算方法和不等長曲線相似度計(jì)算方法。在進(jìn)行等長曲線相似度計(jì)算時(shí)，可以將具有n個(gè)點(diǎn)的曲線看作n維歐氏空間上的一個(gè)點(diǎn)，利用空間點(diǎn)之間的距離公式計(jì)算相似度。常見的距離公式包括歐式距離^［15］、曼哈頓距離^［16］、切比雪夫距離^［17］等;不等長曲線相似度計(jì)算包括豪斯多夫距離^［18］、弗雷歇距離^［19-20］、動態(tài)時(shí)間彎曲距離^［21］等。其中，弗雷歇距離在工程應(yīng)用上的適用性已得到證明，諸如通過對各種鐵路道岔故障曲線的分析進(jìn)行道岔故障診斷^［19］，通過對海上目標(biāo)軌跡歷史曲線和研究曲線的相似度計(jì)算進(jìn)行目標(biāo)軌跡的預(yù)測^［20］。

為了驗(yàn)證弗雷歇距離算法在土阻力推薦方法中的適用性，本文基于某海域11個(gè)平臺99根樁的打樁數(shù)據(jù)庫，計(jì)算了不同折減系數(shù)條件下的土阻力與根據(jù)打樁記錄反算得到的土阻力之間的弗雷歇距離，對利用弗雷歇距離獲得土阻力的方法進(jìn)行了探討，并給出了所研究海域的土阻力折減系數(shù)的推薦值。該研究可為工程實(shí)際土阻力的計(jì)算提供借鑒。

1 樁基可打入性分析方法

樁基的可打入性分析通常是指在現(xiàn)場的沉樁工作之前，基于已有的地質(zhì)資料和樁錘資料采用波動方程的方法，對于未來實(shí)際打樁過程的模擬工作。輸出的內(nèi)容主要包括錘擊數(shù)隨深度的變化情況，打樁時(shí)間，樁身的拉、壓應(yīng)力等^［22-23］。

在采用波動方程的方法進(jìn)行分析時(shí)，需要確定樁的參數(shù)（樁身尺寸，樁身材料），錘的參數(shù)（錘型，沖程，錘效，錘墊和替打）以及土的參數(shù)。其中樁的參數(shù)是已知的，錘的參數(shù)除了錘效也都已知，海上打樁常常采用液壓錘進(jìn)行施工，可打入性分析時(shí)一般取90%錘效。不確定性最大的是土的參數(shù)，主要包括土阻力、彈限和阻尼。計(jì)算土阻力時(shí)往往需要通過前期的鉆孔取樣或者現(xiàn)場原位試驗(yàn)（例如CPT試驗(yàn)或者標(biāo)準(zhǔn)貫入試驗(yàn)）確定，計(jì)算結(jié)果存在很大的不確定性。

1.1 樁基承載力計(jì)算方法

在進(jìn)行樁基可打入性動力分析之前，需要先根據(jù)土力學(xué)理論計(jì)算其靜阻力，即樁基的長期靜承載力，以下簡稱樁基承載力。樁基承載力一般通過API規(guī)范中的計(jì)算方法獲得^［6］。計(jì)算如下所示：

Q_u=Q_s+Q_b=πD∫τ_fdz+A_bq_b （1）

式中：Q_u為樁基長期承載力（kN）;Q_s、Q_b分別為樁側(cè)和樁端承載力（kN）;D為樁端直徑（m）;τ_f為沿樁身分布的單位面積側(cè)摩阻（kPa）;z為沉樁深度（m）;A_b為樁端面積（m²）;q_b為單位面積端阻（kPa）。

對于黏性土，樁基的單位面積側(cè)摩阻τ_f的計(jì)算公式為：

式中：α為無量綱系數(shù)，且α≤1.0;S_u為黏性土的不排水抗剪強(qiáng)度（kPa）;σ^′_v0為土體的有效上覆壓力（kPa）。

黏性土中的樁基單位面積端阻q_b的計(jì)算公式為：

q_b=9S_u （5）

非黏性土中的樁基單位面積側(cè)摩阻τ_f的計(jì)算公式為：

τ_f=Kσ^′_v0tanδ≤τ_f，max （6）

式中：K為橫向地基壓力系數(shù)，也就是水平與豎向有效正應(yīng)力之比，對于開口打入管樁，K的取值為0.8，對于閉口打入樁，K取值為1.0;δ為樁壁和非黏性土之間的摩擦角;τ_f，max是單位面積側(cè)摩阻極限值（kPa），該值與非黏性土類別和密實(shí)度有關(guān)，可根據(jù)樁土摩擦角按經(jīng)驗(yàn)取值。

非黏性土的樁基單位面積端阻q_b的計(jì)算公式為：

q_b=N_qσ^′_v0≤q_b，max （7）

式中：N_q是無量綱支撐系數(shù);q_b，max是單位面積端阻極限值（kPa），兩者可根據(jù)樁土摩擦角按經(jīng)驗(yàn)取值。

對于打樁土阻力的計(jì)算，存在無土塞和有土塞兩種工況，本文采取Alm等^［8］的建議方法：對于無土塞工況，土阻力為樁外壁阻力和樁端環(huán)面積阻力之和；對于有土塞工況，土阻力為樁外壁阻力和樁端全面積阻力之和。

1.2 土阻力計(jì)算方法

打樁過程中的土阻力與長時(shí)間停錘之后的樁基長期承載力Q_u之間可以用一個(gè)損失系數(shù)f_RD來表示：

SRD=f_RD·Q_u （8）

式中：SRD是打樁過程的土阻力（kN）;Q_u是樁基的長期承載力（kN）;f_RD是土阻力相對承載力的損失系數(shù)，無量綱。

基于波動方程的樁基可打入性分析方法中，需要確定損失系數(shù)f_RD，其由兩部分組成：一部分是控制土阻力相對承載力絕對變化的土阻力折減系數(shù)f_GL;另一部分是控制不同土層之間土阻力損失相對變化的土層恢復(fù)系數(shù)f_s。對于土阻力折減系數(shù)f_GL，其代表整個(gè)打樁過程中土阻力相對承載力的整體折減。但是由于施工現(xiàn)場沉樁往往不止貫入一個(gè)土層，不同土層之間土體的土阻力折減比例并不相同，因此引入了土層恢復(fù)系數(shù)f_s，代表不同土層未擾動情況下的強(qiáng)度與完全擾動情況下的強(qiáng)度的比值。對于黏性土，該值與土的靈敏度相關(guān)。

損失系數(shù)f_RD與土阻力折減系數(shù)f_GL、土層恢復(fù)系數(shù)f_s的關(guān)系如下：

f_RD=（1-f^*_s+f^*_s·f_GL） （9）

式中：f^*_s為該土層相對最敏感土層的相對土層恢復(fù)系數(shù)，計(jì)算公式如下所示：

f^*_s=1-1/f_s/1-1/f_sx （10）

式中：f_s為該土層的土層恢復(fù)系數(shù);f_sx為最敏感土層的土層恢復(fù)系數(shù)。

以上公式同樣適用于單一土層情況，即f_sx=f_s的情況，此時(shí)f^*_s=1，f_RD=f_GL，與前文所說的單土層結(jié)果一致。

在計(jì)算過程中，土阻力折減系數(shù)f_GL目前沒有推薦值，需要由工程師根據(jù)經(jīng)驗(yàn)進(jìn)行取值。

由于打樁過程中的土阻力由側(cè)阻和端阻兩部分組成，因此土阻力折減系數(shù)f_GL涉及側(cè)阻折減和端阻折減兩個(gè)部分，而側(cè)阻和端阻這兩個(gè)折減系數(shù)并不相同。例如，Stevens等^［7］對于黏土側(cè)摩推薦的折減系數(shù)為0.5，而端阻不折減，即端阻的折減系數(shù)為1.0。因此一個(gè)土阻力結(jié)果往往需同時(shí)獲得兩個(gè)不同的折減系數(shù)，這就使得確定f_GL的難度很大。

2 弗雷歇距離方法原理

弗雷歇距離也被稱為狗繩距離^［17］，即若人的行走路徑為A，狗的行走路徑為B，兩者走完各自路徑所需的最短狗繩長度即為弗雷歇距離。如圖1所示，情境一和情境二分別代表兩組不同的人、狗速度，A和B分別代表人和狗的路徑，虛線代表不同時(shí)刻的狗繩長度，D_A-D_B和E_A-E_B線段長度分別代表圖1的兩種情境下人、狗走完各自路徑所需的最短狗繩長度。通過改變?nèi)撕凸犯髯缘乃俣?，可以得到不同情境下的人狗走完全部路徑所需的最短狗繩長度，在所有情境下的最短狗繩長度的最小值即為人狗路徑的弗雷歇距離。圖1中的D_A-D_B線段相比E_A-E_B線段更短，更接近弗雷歇距離。

弗雷歇距離的數(shù)學(xué)表達(dá)形式為：

式中：A和B是度量空間S中的兩條曲線;α（t）和β（t）分別是運(yùn)動位置的描述函數(shù);A［α（t）］和B［β（t）］分別是在t時(shí)刻A曲線上人和B曲線上狗的位置;d是S上的度量函數(shù)，這里表示狗繩長度，那么A和B之間的弗雷歇距離F（A，B）就被定義為在t∈［0，1］區(qū)間內(nèi)最長狗繩長度的下確界，也就是任意一對 α（t）和β（t），總能找到最長的狗繩長度，而改變α（t）和β（t）使得最長的狗繩長度最小，這時(shí)的狗繩長度就是弗雷歇距離。利用弗雷歇距離可以對兩條曲線的相似度進(jìn)行判別。

3 弗雷歇距離方法適用性

為了說明弗雷歇距離算法在樁基打樁記錄曲線相似度識別上的適用性，首先運(yùn)用該算法計(jì)算了某海域相同海上平臺和不同平臺的打樁記錄的曲線相似度。選擇的兩個(gè)平臺的工程概況如下所示。

平臺9和平臺3屬于同一海域，兩平臺之間的直線距離為68 km。其中，平臺9的樁基由8根樁組成，樁徑為2.438 m，樁長為124.1 m，設(shè)計(jì)入泥深度為101 m，由MHU 800S液壓沖擊錘貫入;平臺3的樁基由12根樁組成，樁徑為2.438 m，樁長為121.3 m，設(shè)計(jì)入泥深度為98 m，由MHU 1200S液壓沖擊錘貫入。兩平臺現(xiàn)場均進(jìn)行了CPT試驗(yàn)，兩平臺通過CPT試驗(yàn)得到的錐尖阻力隨深度變化如圖2所示。

由圖2可以看出，兩平臺除0～20 m以及80～100 m處土體的錐尖阻力較為接近之外，其他深度處土層的錐尖阻力均存在較大差異，由此可以推斷出兩平臺樁基下的土阻力隨深度曲線也存在較大差異。從平臺9和平臺3的樁基中各抽取兩根樁，根據(jù)打樁記錄反分析結(jié)果得到上述4根樁的土阻力隨深度變化曲線如圖3所示。

從圖3中可以得到，平臺9的樁1和樁2之間，平臺3的樁1和樁2之間曲線的走向趨勢存在很大的相似性。相反，平臺9的兩根樁和平臺3的兩根樁之間相似度不高。采用上一章節(jié)提到的弗雷歇距離計(jì)算方法得到上述4根樁之間的弗雷歇距離結(jié)果如圖4所示。

其中樁1、2取自同一平臺，樁3、4取自另一平臺。4根樁的打樁記錄相互之間的弗雷歇距離計(jì)算結(jié)果在0～52之間，其中同一根樁的打樁記錄的相似度是100%，體現(xiàn)在圖4中就是1-1，2-2，3-3，4-4的弗雷歇距離均為0。弗雷歇距離越大，表明兩個(gè)曲線的相似度越低，圖4中來自同一平臺的樁 （1-2，3-4）弗雷歇距離均在25以下，而不同平臺的樁的弗雷歇距離均在40以上，十分符合實(shí)際的工程情況，證明了弗雷歇距離算法在對這4根樁歸屬平臺上的判別的準(zhǔn)確性。

以上4根樁的結(jié)果說明設(shè)定弗雷歇距離閾值為25～40均可合理劃分不同樁的平臺歸屬，為了明確該海域11個(gè)平臺99根樁的弗雷歇距離閾值，對這99根樁進(jìn)行計(jì)算。11個(gè)平臺的工程背景匯總概況如表1所列。從表1可以看出，11個(gè)平臺的樁基樁徑只有2.134 m和2.438 m兩種情況;樁長匯總為99～124.1 m;設(shè)計(jì)入泥深度匯總為77～101 m;錘型包括三種，分別為IHC S-500液壓沖擊錘、MHU 800S液壓沖擊錘、MHU 1200S液壓沖擊錘;樁基數(shù)量為8或12（平臺11為特殊情況只有3根樁）;錘效分布為42%～98%不等;總錘擊能量分布為659～6 451 MJ不等。

為了驗(yàn)證上述閾值在新數(shù)據(jù)上的適用性，采用機(jī)器學(xué)習(xí)中的劃分訓(xùn)練集和測試集的方法，將平臺10和平臺11的數(shù)據(jù)作為測試集，沒有用于確定閾值。平臺1～9的同一平臺打樁記錄數(shù)據(jù)進(jìn)行了404次計(jì)算，不同平臺打樁記錄數(shù)據(jù)進(jìn)行了3 424次計(jì)算。隸屬于同一平臺和不同平臺的樁基弗雷歇距離統(tǒng)計(jì)情況分別如圖5、6所示。

從圖5、6中可以看出，同一平臺樁基的弗雷歇距離主要集中在20～60之間，不同平臺樁基的弗雷歇距離集中在30～80之間。由于兩者之間存在較大的重合區(qū)域，無法像圖4一樣可以直觀地確定閾值，本文采取了通過分析概率分布的方法確定閾值。對圖5、6中的數(shù)據(jù)進(jìn)行了累積概率分布的計(jì)算，計(jì)算結(jié)果如圖7所示。

為了確定合理的閾值從而確保依據(jù)閾值劃分的同平臺、不同平臺樁基劃分結(jié)果最佳，這里閾值的確定需要同時(shí)考慮同平臺、不同平臺的累積分布曲線。具體來說，就是該閾值能使該弗雷歇距離對應(yīng)的同一平臺正累積概率和不同平臺逆累積概率之和最大。依據(jù)圖7的結(jié)果，閾值的位置確定在32處，該處的同一平臺正累積概率是56.4%，不同平臺逆累積概率是75.8%。也就是將32作為弗雷歇距離的閾值，同一平臺的準(zhǔn)確率為56.4%，不同平臺的準(zhǔn)確率為75.8%，平均準(zhǔn)確率為66.1%。

上述閾值結(jié)果是基于已有的數(shù)據(jù)統(tǒng)計(jì)出來的，為了驗(yàn)證該數(shù)據(jù)在新平臺上的適用性，對平臺10、11的打樁記錄數(shù)據(jù)進(jìn)行了弗雷歇距離計(jì)算，并依據(jù)閾值32進(jìn)行了分類，與實(shí)際結(jié)果進(jìn)行了比較，計(jì)算結(jié)果如表2所列。從表2中可以看出，同一平臺的樁基隸屬劃分準(zhǔn)確率高達(dá)100%，相應(yīng)的不同平臺劃分準(zhǔn)確率為70.8%，總體準(zhǔn)確率為87.3%。弗雷歇距離閾值在新樣本集上的表現(xiàn)是可以接受的，甚至比訓(xùn)練集上的更好。

4 弗雷歇距離用于土阻力折減系數(shù)推薦

弗雷歇距離方法可以用來確定樁基的相似度，進(jìn)而獲得樁基承載力的折減系數(shù)。在確定樁基的相似度中，具體的方法為：

（1） 有土塞和無土塞工況下都對樁基的樁側(cè)折減系數(shù)和樁端折減系數(shù)從0.1到1.0，按0.1增量計(jì)算10×10共計(jì)100次可打入性結(jié)果，兩個(gè)工況共得到200條樁基的土阻力隨深度變化曲線。

（2） 各個(gè)工況下比較實(shí)際的打樁記錄和可打入性結(jié)果，計(jì)算兩者之間的弗雷歇距離，弗雷歇距離最小者對應(yīng)的可打入性結(jié)果即為最接近實(shí)際打樁記錄的結(jié)果，其對應(yīng)的折減系數(shù)也就是該工況下推薦的折減系數(shù)。

以該海域平臺9樁1為例，樁1的100個(gè)有土塞和100個(gè)無土塞可打入性結(jié)果與樁1的實(shí)際打樁記錄之間的弗雷歇距離結(jié)果如圖8所示。如前所述，為了表述清楚，按折減系數(shù)（1.0，1.0）至（1.0，0.1）從1到10編號，以此類推直到（0.1，0.1）編號為100，折減系數(shù)編號和對應(yīng)的側(cè)、端部折減系數(shù)對應(yīng)情況如表3所列。

從圖8中可以看出，基于無土塞假定下的可打入性結(jié)果與實(shí)際打樁記錄的弗雷歇距離隨編號的變大（表現(xiàn)為隨整體土阻力的減小）而減小，變化趨勢相對較緩和，其中最小的弗雷歇距離為49.02，對應(yīng)的折減系數(shù)編號是55，對應(yīng)的折減系數(shù)情況是側(cè)部折減系數(shù)0.5，端部折減系數(shù)0.6;基于有土塞假定的可打入性結(jié)果與實(shí)際打樁記錄之間的弗雷歇距離隨編號的變化存在周期性的波動，其原因在于有土塞的端部土阻力占比相對無土塞更大，從而有土塞的總土阻力也會隨著編號存在周期性的變化。如圖9所示，土阻力和弗雷歇距離隨折減系數(shù)編號的變化趨勢是完全一致的。圖9中有土塞弗雷歇距離最小值為50.06，對應(yīng)的折減系數(shù)編號是70，表示側(cè)部折減系數(shù)為0.4，端阻折減系數(shù)為0.1。具體來說，無土塞中折減系數(shù)編號55對應(yīng)的側(cè)部土阻力為34.1 MN，端部土阻力為2.8 MN，總土阻力為36.9 MN。有土塞折減編號70對應(yīng)的側(cè)部土阻力為27.3 MN，端部土阻力為4.5 MN，總土阻力為31.8 MN，比無土塞的總土阻力小13.8%，兩者總土阻力差異不大。

為了研究弗雷歇距離與土阻力關(guān)系，得到弗雷歇距離隨土阻力的變化情況如圖10所示。根據(jù)反分析結(jié)果得知，該算例在打樁結(jié)束時(shí)的土阻力為31.7 MN，從圖10可以看出，當(dāng)土阻力與現(xiàn)場土阻力最為吻合時(shí)，弗雷歇距離最小。

依據(jù)對平臺9樁1的處理方式，對于平臺9剩下的7根樁做了相同計(jì)算處理，計(jì)算結(jié)果如圖11所示。考慮到有土塞和無土塞是人為確定的兩種情況，且對于同一位置，打樁的土阻力只存在一個(gè)數(shù)值，因此選擇采用兩種工況得到的土阻力結(jié)果接近且出現(xiàn)頻率較高的情況作為該平臺的折減系數(shù)，統(tǒng)計(jì)結(jié)果見表4。

依據(jù)上述方法，對其余10個(gè)平臺的現(xiàn)場實(shí)測打樁記錄和可打入性結(jié)果之間的弗雷歇距離進(jìn)行計(jì)算，無土塞和有土塞的匯總結(jié)果分別如圖12、13所示。從圖12可以看出，無土塞工況下，除了平臺5、7，其余平臺的折減系數(shù)編號集中在55附近;圖13中，有土塞工況下的結(jié)果相對無土塞更集中，除了平臺7以外，其他10個(gè)平臺的數(shù)據(jù)主要都集中在折減系數(shù)編號85附近。即對于該海域，在無土塞工況下，側(cè)摩阻力的折減系數(shù)為0.4～0.6，端阻折減系數(shù)范圍為0.4～1.0;有土塞工況下側(cè)摩阻力的折減系數(shù)范圍為0.1～0.4，端阻折減系數(shù)范圍為0.1～0.6。

得到上述推薦折減系數(shù)之后，利用平臺9的打樁記錄反算得到的土阻力隨深度曲線對推薦的折減系數(shù)進(jìn)行驗(yàn)證，驗(yàn)證方法如下所示：

（1）選擇推薦的折減系數(shù)上限和下限得到4種工況，分別計(jì)算4種工況下的土阻力隨深度變化結(jié)果，如表5所列（以平臺9為例）。

（2） 再將每個(gè)平臺的實(shí)際打樁記錄反算的土阻力結(jié)果與依據(jù)上述計(jì)算條件得到的土阻力結(jié)果相比較，驗(yàn)證上述折減參數(shù)是否合理。

平臺9的8根樁的反算土阻力結(jié)果如圖14所示。從圖14可以看出，盡管個(gè)別深度部分樁的土阻力在數(shù)值上不完全相同，但是整體上平臺9上8根樁的反算土阻力結(jié)果具有很好的一致性。

將表5中的4種工況的計(jì)算結(jié)果與圖14中反算土阻力結(jié)果比較，比較結(jié)果如圖14所示。從圖14可以看出，無土塞工況下，側(cè)部和端部折減系數(shù)分別為0.4/0.4的土阻力結(jié)果與平臺9上8根樁的反算土阻力結(jié)果均值更接近;有土塞工況下，側(cè)部和端部折減系數(shù)分別為0.1/0.1的土阻力結(jié)果整體小于實(shí)際反算土阻力，也就是計(jì)算結(jié)果偏危險(xiǎn)，該種工況更適合用于預(yù)測溜樁等不能高估土阻力的工況;無土塞工況下，側(cè)部和端部折減系數(shù)分別為0.6/1.0和有土塞工況下，側(cè)部和端部折減系數(shù)分別為0.4/0.6的土阻力結(jié)果整體大于實(shí)際反算土阻力，也就是計(jì)算結(jié)果偏安全。

綜合圖14的結(jié)果可以看出，表5推薦的折減系數(shù)可為實(shí)際工程的折減系數(shù)的選取提供借鑒。

5 結(jié)語

本文基于弗雷歇距離方法的曲線相似度算法，提出了海上平臺沉樁過程中的土阻力折減系數(shù)推薦方法，并得出以下結(jié)論：

（1） 11個(gè)平臺上99根樁之間的打樁記錄曲線弗雷歇距離的計(jì)算結(jié)果顯示，同平臺、不同平臺之間的弗雷歇距離結(jié)果差異明顯。當(dāng)以弗雷歇距離32作為判斷是否同一平臺的閾值時(shí)，9個(gè)平臺結(jié)果中同一平臺的準(zhǔn)確率為56.4%，不同平臺的準(zhǔn)確率為75.8%。將該閾值結(jié)果應(yīng)用于剩下2個(gè)平臺時(shí)，同一平臺的準(zhǔn)確率為100%，不同平臺的準(zhǔn)確率為70.8%。利用弗雷歇距離進(jìn)行曲線相似度判別具有可行性。

（2） 在進(jìn)行樁基可打入性分析中采用有土塞和無土塞工況，均可獲得與現(xiàn)場土阻力接近的結(jié)果。盡管在不同的工況下，所采用的折減系數(shù)不同，但兩者計(jì)算得到的打樁過程中的總土阻力數(shù)值接近。

（3） 綜合所有計(jì)算結(jié)果，某海域在進(jìn)行可打入性分析時(shí)，無土塞工況下，側(cè)部折減系數(shù)范圍為0.4～0.6，端阻折減系數(shù)范圍為0.4～1.0;有土塞工況下側(cè)摩阻力的折減系數(shù)范圍為0.1～0.4，端阻折減系數(shù)范圍為0.1～0.6。

參考文獻(xiàn)（References）

［1］閆澍旺，李嘉，賈沼霖，等.海洋石油平臺超長樁拒錘分析及工程實(shí)例［J］.巖土力學(xué)，2015，36（增刊2）：559-564.

YAN Shuwang，LI Jia，JIA Zhaolin，et al.Analysis and engineering case of super-long pile refused hammer in offshore oil platform［J］.Rock and Soil Mechanics，2015，36（Suppl02）：559-564.

［2］李颯，王耀存，蒲玉成，等.海洋平臺打樁過程中溜樁對樁基影響的研究［J］.地震工程學(xué)報(bào)，2014，36（3）：462-467.

LI Sa，WANG Yaocun，PU Yucheng，et al.Influence of pile sinking on pile capacity during pile driving on offshore platforms［J］.China Earthquake Engineering Journal，2014，36（3）：462-467.

［3］BYRNE T，GAVIN K，PRENDERGAST L J，et al.Performance of CPT-based methods to assess monopile driveability in North Sea sands［J］.Ocean Engineering，2018，166：76-91.

［4］LI S，TIAN Y H，ZHOU Y R，et al.Premature refusal of large-diameter，deep-penetration piles on an offshore platform［J］.Applied Ocean Research，2013，42：55-59.

［5］SEMPLE R M，GEMEINHARDT J P.Stress history approach to analysis of soil resistance to pile driving［C］//Proceedings of the Offshore Technology Conference.Houston，Texas：＼[s.n.＼]，1981：165-172.

［6］American Petroleum Institute.WSD Recommended practice for planning，designing and constructing fixed offshore platforms-working stress design （21st ed.）［S］.USA：American Petroleum Institute，2005.

［7］STEVENS R S，WILTSIE E A，TURTON T H.Evaluating drivability for hard clay，very dense sand，and rock［C］//Proceedings of the Offshore Technology Conference.Houston，Texas：＼[s.n.＼]，1982：465-481.

［8］ALM，T，HAMRE L.Soil model for pile driveability predictions based on CPT interpretations［C］// Proceedings of 15th International Conference on Soil Mechanics and Geotechnical Engineering.Istanbul：＼[s.n.＼]，2001：27-31.

［9］PRENDERGAST L J，GANDINA P，GAVIN K.Factors influencing the prediction of pile driveability using CPT-based approaches［J］.Energies，2020，13（12）：3128.

［10］JARDINE R，CHOW F，OVERY R，et al.ICP design methods for driven piles in sands and clays［M］.London：Thomas Telford，2005.

［11］LEHANE B M，SCHNEIDER J A，XU X.Development of the UWA-05 design method for open and closed ended driven piles in siliceous sand［C］//Proceedings of the Contemporary Issues in Deep Foundations.Denver，Colorado，USA：American Society of Civil Engineers，2007：1-10.

［12］KOLK H J，BAAIJENS A E，SENDERS M.Design criteria for pipe piles in silica sands［C］// Proceedings of the First International Symposium on Frontiers in Offshore Geotechnics.Taylor，F(xiàn)rancis，London：Press of University of Western Australia，2005：711-716.

［13］李颯，吳興州，王耀存，等.打樁過程中考慮溜樁影響的土阻力研究［J］.巖土工程學(xué)報(bào)，2015，37（6）：1150-1157.

LI Sa，WU Xingzhou，WANG Yaocun，et al.Soil resistance to driving considering effect of pile running on pile installation［J］.Chinese Journal of Geotechnical Engineering，2015，37（6）：1150-1157.

［14］QIN W，CAI S Y，DAI G L，et al.Soil resistance during driving of offshore large-diameter open-ended thin-wall pipe piles driven into clay by impact hammers［J］.Computers and Geotechnics，2023，153：105085.

［15］黨倩，崔阿軍，尚聞博，等.采用歐式形態(tài)距離的負(fù)荷曲線近鄰傳播聚類方法［J］.西安交通大學(xué)學(xué)報(bào)，2022，56（1）：165-176.

DANG Qian，CUI Ajun，SHANG Wenbo，et al.Affinity propagation clustering method of typical load curve with euclidean morphological distance［J］.Journal of Xi'an Jiaotong University，2022，56（1）：165-176.

［16］胡濤，任仙怡，張基宏.高分辨率PCB裸板圖像輪廓快速DP曲線抽?。跩］.光電工程，2013，40（7）：83-88.

HU Tao，REN Xianyi，ZHANG Jihong.A fast DP polyline simplification algorithm for high-resolution PCB bare board image contours［J］.Opto-Electronic Engineering，2013，40（7）：83-88.

［17］黎佳，陳元坤，陳曉強(qiáng)，等.基于切比雪夫網(wǎng)格的自由曲面網(wǎng)格自適應(yīng)及桿長優(yōu)化［J］.土木工程與管理學(xué)報(bào)，2020，37（3）：91-95，102.

LI Jia，CHEN Yuankun，CHEN Xiaoqiang，et al.Adaptive and length optimization research on free-form surface grid generation based on Chebyshev net［J］.Journal of Civil Engineering and Management，2020，37（3）：91-95，102.

［18］劉冠蘭，徐亞明，柏文峰.地鐵隧道斷面變形分析方法研究［J］.測繪通報(bào)，2016（2）：103-106，135.

LIU Guanlan，XU Yaming，BAI Wenfeng.A method to determine the deformation of tunnel section［J］.Bulletin of Surveying and Mapping，2016（2）：103-106，135.

［19］黃世澤，陳威，張帆，等.基于弗雷歇距離的道岔故障診斷方法［J］.同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版），2018，46（12）：1690-1695.

HUANG Shize，CHEN Wei，ZHANG Fan，et al.Method of turnout fault diagnosis based on Fréchet distance［J］.Journal of Tongji University （Natural Science），2018，46（12）：1690-1695.

［20］趙一鑒，林利，王茜蒨，等.基于大地距離計(jì)算相似度的海上目標(biāo)軌跡預(yù)測［J］.計(jì)算機(jī)應(yīng)用，2023，43（11）：3594-3598.

ZHAO Yijian，LIN Li，WANG Xiqian，et al.Trajectory prediction of sea targets based on geodetic distance similarity calculation［J］.Journal of Computer Applications，2023，43（11）：3594-3598.

［21］李映輝，鐘平安，錢睿智，等.基于動態(tài)時(shí)間彎曲算法的相似洪水識別方法［J］.水電能源科學(xué)，2018，36（11）：51-55.

LI Yinghui，ZHONG Ping'an，QIAN Ruizh，et al.Similar flood recognition method based on dynamic time warping algorithm［J］.Water Resources and Power，2018，36（11）：51-55.

［22］郭全生，王向平，付永剛，等.大直徑樁波動方程可打性分析及應(yīng)用［J］.地震工程學(xué)報(bào)，2017，39（增刊1）：232-235.

GUO Quansheng，WANG Xiangping，F(xiàn)U Yonggang，et al.Drivability analysis of offshore large diameter piles using a wave equation analysis program and its application［J］.China Earthquake Engineering Journal，2017，39（Suppl01）：232-235.

［23］王耀存，李颯，張培獻(xiàn)，等.打樁對樁周土體阻尼系數(shù)的影響［J］.地震工程學(xué)報(bào)，2014，36（4）：946-951.

WANG Yaocun，LI Sa，ZHANG Peixian，et al.Influence of piling on damping coefficients of soil surrounding piles［J］.China Earthquake Engineering Journal，2014，36（4）：946-951.

（本文編輯：賈源源）