張 偉，宋夢(mèng)然，黃 焱，馬君誠(chéng)

(1.海洋石油工程股份有限公司，天津 300451； 2.天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072)

?

八樁腿大跨度導(dǎo)管架平臺(tái)冰激振動(dòng)模型試驗(yàn)

張 偉1，宋夢(mèng)然2，黃 焱2，馬君誠(chéng)1

(1.海洋石油工程股份有限公司，天津 300451； 2.天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072)

針對(duì)目前關(guān)于導(dǎo)管架平臺(tái)大跨度柔性結(jié)構(gòu)冰激振動(dòng)研究不夠充分的問題，對(duì)我國(guó)渤海海域某帶錐八樁腿導(dǎo)管架平臺(tái)開展系統(tǒng)的冰激振動(dòng)模型試驗(yàn)研究。研究發(fā)現(xiàn)：平臺(tái)在不同冰攻角、冰速、冰厚和水位條件下的振動(dòng)類型及響應(yīng)水平差別較大，5年一遇冰厚條件下，結(jié)構(gòu)的最大響應(yīng)出現(xiàn)在22.5°冰攻角、原型冰速0.8 m/s時(shí)；百年一遇冰厚條件下，結(jié)構(gòu)的最大響應(yīng)出現(xiàn)在22.5°冰攻角、原型冰速1.1 m/s時(shí)。為解決具體的工程實(shí)際問題提供參考。

錐型結(jié)構(gòu)；導(dǎo)管架平臺(tái)；冰激振動(dòng)；模型試驗(yàn)

0 引言

大跨度的導(dǎo)管架平臺(tái)在冰排與結(jié)構(gòu)發(fā)生相互作用時(shí)產(chǎn)生的載荷是交變的，即為一種動(dòng)力作用，這將導(dǎo)致結(jié)構(gòu)發(fā)生劇烈的振動(dòng)。隨著油氣開采業(yè)在寒區(qū)的發(fā)展，越來越多的具有顯著水平的冰激振動(dòng)事件威脅著平臺(tái)的正常作業(yè)及結(jié)構(gòu)的安全。例如，1977年，“海四井”烽火臺(tái)被流冰推倒，生活平臺(tái)振動(dòng)劇烈，平臺(tái)棧橋難以行走[1]；1999年-2000年冬季，渤海JZ20-2中南平臺(tái)由于冰激振動(dòng)導(dǎo)致了平臺(tái)管線的斷裂[2]。

隨著冰邊緣區(qū)范圍擴(kuò)大，作用于渤海海域?qū)Ч芗芷脚_(tái)上的冰力極值有所減小。然而隨著海洋動(dòng)力學(xué)條件更為復(fù)雜，導(dǎo)管架平臺(tái)在海冰作用下的疲勞問題日益顯著[3-4]。由于天然海冰的彎曲強(qiáng)度遠(yuǎn)小于擠壓強(qiáng)度，引導(dǎo)冰排由在直立結(jié)構(gòu)前的擠壓破壞轉(zhuǎn)為在錐形結(jié)構(gòu)前的彎曲破壞，可以有效降低冰載荷[5]。因此，許多寒區(qū)的海洋工程結(jié)構(gòu)物都安裝破冰錐體。然而，安裝破冰錐體可能會(huì)導(dǎo)致一些問題，冰排發(fā)生彎曲破壞的同時(shí)使冰力具備了周期特性，這樣一來，錐體冰力的動(dòng)力特征就變得尤為突出。本文針對(duì)我國(guó)渤海海域的某八樁腿帶錐導(dǎo)管架平臺(tái)的冰激振動(dòng)問題進(jìn)行系統(tǒng)的模型試驗(yàn)研究。

1 冰激錐型結(jié)構(gòu)振動(dòng)概述

無(wú)論是從冰與錐體相互作用的現(xiàn)象來看，還是從冰與錐體相互作用的冰力時(shí)程曲線來看，冰力模型試驗(yàn)和現(xiàn)場(chǎng)原型觀測(cè)都得到了同樣的結(jié)論，即冰排與錐體結(jié)構(gòu)的相互作用具有周期性，會(huì)產(chǎn)生間歇的周期載荷。對(duì)冰激錐體振動(dòng)問題，影響確定性冰力函數(shù)的關(guān)鍵因素有3個(gè)：冰排斷裂長(zhǎng)度、冰速和冰力幅值。

同時(shí)，冰排與錐體的作用過程中存在多次斷裂進(jìn)程，導(dǎo)致冰力中含有不止一種激勵(lì)頻率成分，而當(dāng)結(jié)構(gòu)自振頻率與其中一種頻率成分相匹配的時(shí)候，就可能引發(fā)劇烈的結(jié)構(gòu)振動(dòng)。此外，當(dāng)冰排與多錐體結(jié)構(gòu)相互作用時(shí)，冰排在不同錐體結(jié)構(gòu)間的破壞進(jìn)程存在非同時(shí)性，這種非同時(shí)效應(yīng)弱化了冰激振動(dòng)的水平。

因此，冰激振動(dòng)進(jìn)程是受到冰排破壞模式與結(jié)構(gòu)條件的共同影響和控制的。也就是說，在不同的冰條件與結(jié)構(gòu)條件組合情況下，將發(fā)生不同冰激振動(dòng)機(jī)理所控制的振動(dòng)事件。

2 模型試驗(yàn)概況

2.1 試驗(yàn)?zāi)Ｐ?/p>

選擇順應(yīng)式模擬器作為試驗(yàn)的主要裝置，它可以準(zhǔn)確地模擬結(jié)構(gòu)自身的動(dòng)力特性，如圖1所示。試驗(yàn)?zāi)Ｐ桶凑?∶30的縮尺比加工。表1為試驗(yàn)中模型結(jié)構(gòu)動(dòng)力特性參數(shù)與原型參數(shù)的對(duì)比，其中模型結(jié)構(gòu)動(dòng)力特性參數(shù)的設(shè)定依據(jù)試驗(yàn)的相似體系由原型參數(shù)進(jìn)行縮尺得到。模型固定板直接與順應(yīng)式模擬器裝配，如圖2所示。

圖1 順應(yīng)式模擬器 圖2 裝配后的平臺(tái)冰激振動(dòng)試驗(yàn)?zāi)Ｐ?/p>

表1 冰激振動(dòng)結(jié)構(gòu)動(dòng)力特征參數(shù)與原型對(duì)應(yīng)情況

2.2 試驗(yàn)條件

冰厚和冰強(qiáng)度的確定是針對(duì)平臺(tái)所處海域冰區(qū)的5年一遇和百年一遇條件設(shè)定的。按照模型試驗(yàn)相似體系的幾何比尺1∶30，確定試驗(yàn)?zāi)繕?biāo)冰厚、冰強(qiáng)度與原型間的對(duì)比關(guān)系，見表2。模型冰采用第二代尿素模型冰。

表2 目標(biāo)冰厚、冰強(qiáng)度參數(shù)與原型的對(duì)應(yīng)情況

在平臺(tái)所在海域可能出現(xiàn)的冰速范圍內(nèi)選取8種工況進(jìn)行試驗(yàn)。該海域海冰最大漂移速度為1.4 m/s，縮尺后最大試驗(yàn)冰速為255 mm/s，其余7種冰速則按照線性遞增的方式確定為20 mm/s，50 mm/s，80 mm/s，100 mm/s，120 mm/s，150 mm/s和200 mm/s。

不同的冰攻角將產(chǎn)生不同的破壞模式、冰載荷和冰激振動(dòng)現(xiàn)象。根據(jù)平臺(tái)所在海域環(huán)境條件，在模型試驗(yàn)中共設(shè)置2種冰攻角作用條件：0°和22.5°。

3 試驗(yàn)現(xiàn)象

圖3 冰排在平臺(tái)模型前的非同時(shí)破壞(百年遇冰厚、平均水位、冰攻角0°)

在平均水位工況下，冰排的下邊緣與正-倒錐組合體交接面位置發(fā)生接觸。在這種情況下，冰排均被正錐體表面順利挑起，繼而發(fā)生徑向裂紋與環(huán)向裂紋并存的楔形梁式彎曲破壞。然而，這種破壞形式在不同的冰攻角下有不同的體現(xiàn)。

在0°冰攻角下，位于平臺(tái)兩側(cè)的2座導(dǎo)管架結(jié)構(gòu)上各有2個(gè)抗冰錐體首先迎冰。冰排在錐體模型前發(fā)生彎曲破壞，冰排在4個(gè)迎冰錐體前的非同時(shí)破壞現(xiàn)象可通過肉眼觀察，如圖3所示。在22.5°冰攻角下，位于平臺(tái)兩側(cè)的2座導(dǎo)管架結(jié)構(gòu)上各有3個(gè)抗冰錐體迎冰，而這3個(gè)錐體間的位置關(guān)系并不在一條直線上，這就造成冰排在這一工況下的非同時(shí)破壞效應(yīng)更加突顯。

值得關(guān)注的是，上述非同時(shí)破壞現(xiàn)象均在低冰速(20～150 mm/s)試驗(yàn)中表現(xiàn)得較為明顯，但在高冰速(150～255 mm/s)試驗(yàn)中，很難通過肉眼觀察到顯著的非同時(shí)現(xiàn)象。

4 試驗(yàn)結(jié)果分析

4.1 測(cè)試結(jié)果時(shí)程特征分析

以平均水位、百年一遇冰厚、0°冰攻角工況為例，分別針對(duì)不同冰速下結(jié)構(gòu)的振動(dòng)時(shí)程特征進(jìn)行分析。

4.1.1 低冰速(20～150 mm/s)試驗(yàn)

試驗(yàn)現(xiàn)象中描述的非同時(shí)效應(yīng)可以通過對(duì)冰力時(shí)程的觀察得到更好的體現(xiàn)。以冰排作用方向?yàn)榛鶞?zhǔn)，圖4為百年一遇冰厚、0°冰攻角、50 mm/s冰速試驗(yàn)中左、右兩側(cè)結(jié)構(gòu)上的冰力時(shí)程曲線。2條冰力時(shí)程曲線盡管在冰力幅值和波動(dòng)形態(tài)上相差不大，但二者在波動(dòng)相位上的差異卻十分明顯。將左右兩側(cè)的冰力時(shí)程分別作快速傅里葉變換(Fast Fourier Transform,FFT)處理，即可得到各自的頻譜特征，如圖5和圖6所示。

圖4 冰力時(shí)程曲線 圖5 左側(cè)冰力時(shí)程頻譜

圖6 右側(cè)冰力時(shí)程頻譜

觀察2條冰力時(shí)程頻譜曲線可以發(fā)現(xiàn)：波動(dòng)相位上的差異僅僅致使2條冰力時(shí)程的頻譜特征具有細(xì)微的差別，且頻譜分析結(jié)果中均現(xiàn)了3個(gè)冰力特征頻率。冰排在與大傾角錐體結(jié)構(gòu)作用時(shí)將以二次斷裂的模式發(fā)生破壞，進(jìn)而致使冰力中存在2種主要的波動(dòng)頻率成分。在2條冰力頻譜曲線中，前兩階的冰力主頻保持一致，即0.3 Hz和0.6 Hz。相同的特征頻率說明二者由相同的冰條件和結(jié)構(gòu)條件引導(dǎo)了相同的冰排二次斷裂進(jìn)程。2條冰力頻譜曲線在第三階冰力主頻上出現(xiàn)了差異，分別為0.88 Hz和0.96 Hz。

在具備上述穩(wěn)定波動(dòng)特征的冰力載荷作用下，具備穩(wěn)定動(dòng)力特征的結(jié)構(gòu)就必然會(huì)呈現(xiàn)出與之相對(duì)應(yīng)的特有響應(yīng)形態(tài)。圖7為位移響應(yīng)時(shí)程曲線，圖8為頻譜分析結(jié)果，由圖7和圖8可以看到：模型結(jié)構(gòu)在吸收了集中在0～1 Hz范圍內(nèi)的冰力激振能量后，開始向高頻振動(dòng)上傳輸能量。

圖7 位移響應(yīng)時(shí)程曲線 圖8 位移響應(yīng)時(shí)程頻譜

圖9 冰力時(shí)程曲線

4.1.2 中冰速(150～200 mm/s)試驗(yàn)

隨著冰排作用速度的提升，非同時(shí)破壞效應(yīng)的顯著程度以及影響水平也隨之下降。同樣以冰排作用方向?yàn)榛鶞?zhǔn)，圖9為百年一遇冰厚、0°冰攻角、150 mm/s冰速試驗(yàn)中左、右兩側(cè)模型結(jié)構(gòu)上的冰力時(shí)程曲線?？梢钥闯觯?條冰力時(shí)程曲線在波動(dòng)相位上的差異較低冰速試驗(yàn)時(shí)呈現(xiàn)出明顯的弱化跡象。

圖10為結(jié)構(gòu)的位移響應(yīng)時(shí)程頻譜曲線，可以看出：平臺(tái)模型結(jié)構(gòu)的位移時(shí)程開始呈現(xiàn)出較為穩(wěn)定的振動(dòng)跡象，此時(shí)結(jié)構(gòu)的一階模態(tài)已被激發(fā)，這樣結(jié)構(gòu)將形成對(duì)外部輸入能量的一次過濾，將遠(yuǎn)離其一階自振的激振能量排除，而將大部分的振動(dòng)能量集中在一階模態(tài)上。由此可將這種情況下的結(jié)構(gòu)振動(dòng)歸結(jié)為超諧波共振。

平臺(tái)整體模型結(jié)構(gòu)的轉(zhuǎn)角響應(yīng)時(shí)程如圖11所示，可以看出：此時(shí)盡管冰力作用水平得到了提升，但平臺(tái)模型的整體轉(zhuǎn)角響應(yīng)幅值較之低冰速情況時(shí)并無(wú)明顯升高，而其振動(dòng)的平穩(wěn)程度卻明顯高于低冰速工況。

圖10 位移響應(yīng)時(shí)程頻譜

圖11 轉(zhuǎn)角響應(yīng)時(shí)程

圖12 冰速冰力時(shí)程曲線

4.1.3 高冰速(200～255 mm/s)試驗(yàn)

在冰速達(dá)到高速區(qū)后，冰排的破壞速度也達(dá)到了峰值，非同時(shí)破壞效應(yīng)的顯著程度以及影響水平也隨之下降至最低幾乎可以忽略的水平。同樣以冰排作用方向?yàn)榛鶞?zhǔn)，圖12中為百年一遇冰厚、0°冰攻角、255 mm/s冰速試驗(yàn)中左、右兩側(cè)模型結(jié)構(gòu)上的冰力時(shí)程曲線,可以看出：2條冰力時(shí)程曲線在波動(dòng)相位上的差異基本上已經(jīng)可以忽略。

結(jié)構(gòu)在上述冰力作用下的位移響應(yīng)及其頻譜特征分別如圖13和圖14所示，可以看出：此時(shí)結(jié)構(gòu)的位移幅值明顯升高，同時(shí)，頻譜分析結(jié)果顯示此時(shí)的位移響應(yīng)在結(jié)構(gòu)自身的一階模態(tài)上出現(xiàn)了明顯的集中。

圖13 位移響應(yīng)時(shí)程 圖14 位移響應(yīng)時(shí)程頻譜

圖15 轉(zhuǎn)角響應(yīng)時(shí)程

平臺(tái)整體模型結(jié)構(gòu)的轉(zhuǎn)角響應(yīng)時(shí)程如圖15所示，可以看出：結(jié)構(gòu)響應(yīng)在一階模態(tài)頻率上的集中制約了因冰排非同時(shí)破壞效應(yīng)引發(fā)的扭轉(zhuǎn)響應(yīng)，致使其下降至很低的水平。

上述針對(duì)0°冰攻角試驗(yàn)測(cè)試結(jié)果分析得到的結(jié)論同樣存在于22.5°冰攻角系列試驗(yàn)中，因此這里不再贅述。

4.2 測(cè)試結(jié)果統(tǒng)計(jì)分析

為了針對(duì)測(cè)試物理參量的統(tǒng)計(jì)結(jié)果進(jìn)行清晰的分析，這里定義統(tǒng)計(jì)參數(shù)：

位移波動(dòng)范圍Δd為

式中：dmax為每組次測(cè)試記錄中位移響應(yīng)的最大值；dmin為每組次測(cè)試記錄中位移響應(yīng)的最小值。

加速度波動(dòng)范圍Δa為

式中：amax為每組次測(cè)試記錄中加速度響應(yīng)的最大值；amin為每組次測(cè)試記錄中加速度響應(yīng)的最小值。

轉(zhuǎn)角波動(dòng)范圍Δθ為

式中：θmax為每組次測(cè)試記錄中轉(zhuǎn)角響應(yīng)的最大值；θmin為每組次測(cè)試記錄中轉(zhuǎn)角響應(yīng)的最小值。

4.2.1 0°冰攻角試驗(yàn)

圖16和圖17分別為2種(5年一遇和百年一遇)目標(biāo)冰厚試驗(yàn)中位移波動(dòng)范圍在3種水位工況下隨冰速的變化趨勢(shì)。

圖16 位移波動(dòng)范圍隨冰速的變化曲線(0°冰攻角、5年一遇目標(biāo)冰厚) 圖17 位移波動(dòng)范圍隨冰速的變化曲線(0°冰攻角、百年一遇目標(biāo)冰厚)

由圖16和圖17可以看出：總體上，平臺(tái)結(jié)構(gòu)整體的位移響應(yīng)水平在最高天文潮水位下最高，最低天文潮水位下的平臺(tái)結(jié)構(gòu)整體位移響應(yīng)水平低于另外2個(gè)水位情況；5年一遇目標(biāo)冰厚條件下，平臺(tái)結(jié)構(gòu)整體的最大位移響應(yīng)出現(xiàn)在150 mm/s冰速條件下，對(duì)應(yīng)原型冰速條件為0.8 m/s；百年一遇目標(biāo)冰厚條件下，平臺(tái)結(jié)構(gòu)整體的最大位移響應(yīng)出現(xiàn)在200 mm/s冰速條件下，對(duì)應(yīng)原型冰速條件為1.1 m/s。

值得關(guān)注的問題是，平臺(tái)結(jié)構(gòu)整體的最大位移所對(duì)應(yīng)的冰速條件是隨冰厚的變化而變化的。由于冰排的初始彎曲破壞是促使結(jié)構(gòu)形成高頻濾波響應(yīng)形態(tài)的關(guān)鍵，那么，冰排的初始彎曲破壞所形成的冰排破壞頻率距結(jié)構(gòu)的自振頻率越近，則冰載荷在與結(jié)構(gòu)自振頻率相近的低頻范圍內(nèi)形成的能量集中水平也就越高。自然，結(jié)構(gòu)對(duì)高頻激振能量的濾波功能會(huì)在上述低頻能量集中水平最高的時(shí)刻達(dá)到最佳，而結(jié)構(gòu)的動(dòng)力響應(yīng)也會(huì)在這一時(shí)刻達(dá)到最大水平。由此可見，臨界冰速的出現(xiàn)取決于冰排的初始彎曲破壞頻率與結(jié)構(gòu)自振頻率的接近程度。冰排的彎曲破壞頻率取決于冰速與冰厚這2種基本作用條件，當(dāng)冰排厚度發(fā)生變化時(shí)，導(dǎo)致結(jié)構(gòu)最大響應(yīng)的臨界冰速也必將隨之發(fā)生變化。

圖18和圖19分別顯示了2種冰厚試驗(yàn)中加速度波動(dòng)范圍在3種水位工況下隨冰速的變化趨勢(shì)，可以看出：在平均水位條件下，加速度響應(yīng)水平總體上隨冰速的增長(zhǎng)而上升，這與測(cè)試時(shí)程特征分析的結(jié)論是一致的。

圖18 轉(zhuǎn)角波動(dòng)范圍隨冰速的變化分布 (0°冰攻角、5年一遇目標(biāo)冰厚) 圖19 轉(zhuǎn)角波動(dòng)范圍隨冰速的變化分布 (0°冰攻角、百年遇目標(biāo)冰厚)

4.2.2 22.5°冰攻角試驗(yàn)

圖20～圖23分別為22.5°冰攻角系列試驗(yàn)結(jié)果，可以看出：22.5°冰攻角系列試驗(yàn)中3項(xiàng)響應(yīng)波動(dòng)范圍的變化規(guī)律與0°冰攻角系列試驗(yàn)中的規(guī)律相似，區(qū)別僅在于22.5°冰攻角系列試驗(yàn)的響應(yīng)水平要高于0°冰攻角系列試驗(yàn)。由此可見，所面臨的冰排作用方向變化并不會(huì)致使結(jié)構(gòu)振動(dòng)機(jī)理及模式上發(fā)生改變，但結(jié)構(gòu)承受冰載荷的樁腿數(shù)目增多，會(huì)造成平臺(tái)結(jié)構(gòu)整體響應(yīng)水平升高。

圖20 位移波動(dòng)范圍隨冰速的變化分布 (22.5°冰攻角、5年一遇目標(biāo)冰厚) 圖21 位移波動(dòng)范圍隨冰速的變化分布 (22.5°冰攻角、百年一遇目標(biāo)冰厚)

圖22 轉(zhuǎn)角波動(dòng)范圍隨冰速的變化分布 (22.5°冰攻角、5年一遇目標(biāo)冰厚) 圖23 轉(zhuǎn)角波動(dòng)范圍隨冰速的變化分布 (22.5°冰攻角、百年一遇目標(biāo)冰厚)

5 結(jié)論

(1) 非同時(shí)破壞效應(yīng)的顯著程度以及影響水平隨著冰排作用速度的提升逐漸減小，并且位移響應(yīng)在結(jié)構(gòu)自身的一階模態(tài)上逐漸集中。

(2) 當(dāng)冰排厚度發(fā)生變化時(shí)，結(jié)構(gòu)最大響應(yīng)的臨界冰速隨之發(fā)生變化。在5年一遇冰厚條件下，結(jié)構(gòu)的最大響應(yīng)出現(xiàn)在22.5°冰攻角、原型冰速0.8 m/s時(shí)；在百年一遇冰厚條件下，結(jié)構(gòu)的最大響應(yīng)出現(xiàn)在22.5°冰攻角、原型冰速1.1 m/s時(shí)。

[1] 段夢(mèng)蘭,方華燦，陳如恒，等.渤海老二號(hào)平臺(tái)被冰推倒的調(diào)查結(jié)論[J].石油礦場(chǎng)機(jī)械.1994,23(3)：1-4.

[2] 黃焱.冰激海洋平臺(tái)振動(dòng)的動(dòng)力模型試驗(yàn)研究[D].天津：天津大學(xué),2004.

[3] HUANG Y,MA J J,TIAN Y F.Model Tests of Four-Legged Jacket Platforms in Ice:Part 1.Model Tests and Results[J].Cold Region Science &Technology,2013,95(11):74-85.

[4] HUANG Y,YU M,TIAN Y F.Model Tests of Four-Legged Jacket Platforms in Ice:Part 2.Analyses and Discussions[J].Cold Region Science &Technology,2013,95(11):86-101.

[5] SODHI D S.Crushing Failure During Ice-Structure Interaction[J].Engineering Fracture Mechanics,2001,68(17):1889-1921.

Model Test Study on Ice-Induced Vibration of Eight-Legged
Long Span Jacket Platform

ZHANG Wei1,SONG Mengran2,HUANG Yan2,MA Juncheng1

(1.Offshore Oil Engineering Co.,Ltd.,Tianjin 300451,China； 2.State Key Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300072,China)

Because the research on ice induced vibration of flexible structures with large span of jacket platform is not sufficient,a series of model tests on the ice-induced vibration of an eight-legged conical structure jacket platform in the Bohai Sea are researched.It can be found that the platform in different ice angle of attack,ice velocity,ice thickness and water level has big differences of vibration types and response level.Under the condition of ice thickness generated every 5 years,the maximum response of the structure appeared in the 22.5° angle of ice attack with the prototype ice speed of 0.8 m/s.Under the condition of ice thickness generated every 100 years,the maximum response of the structure appeared in the 22.5° angle of ice attack with the prototype ice speed of 1.1 m/s.It has positive significance on solving the engineering specific problem.

conical structure；jacket platform；ice-induced vibration；model test

2016-11-09

張 偉(1982-)，男，工程師

1001-4500(2017)03-0086-07

U65

A