王銘基，孫珊珊，謝孜源，王延林

（1.大連理工大學(xué)運載工程與力學(xué)學(xué)部，遼寧大連 116024；2.大連理工大學(xué)海洋科學(xué)與技術(shù)學(xué)院，遼寧盤錦 124221；3.大連理工大學(xué)物理學(xué)院，遼寧盤錦 124221）

0 引 言

海洋石油平臺遠(yuǎn)離陸地，其一旦發(fā)生意外，不但會造成巨大的經(jīng)濟(jì)損失及人員傷亡，而且還將對海洋環(huán)境造成災(zāi)難性的破壞，因此海洋石油平臺的安全問題尤為重要。寒區(qū)海上油氣開發(fā)中，海冰對油氣開發(fā)設(shè)施的破壞能力要遠(yuǎn)遠(yuǎn)高于波浪、海流、風(fēng)等其它常規(guī)環(huán)境荷載。冰荷載是冰區(qū)海洋結(jié)構(gòu)設(shè)計與使用過程中需要重點考慮的外荷載，是冰與結(jié)構(gòu)作用發(fā)生破碎而產(chǎn)生的[1]。海冰與柔性導(dǎo)管架平臺作用會引起導(dǎo)管架結(jié)構(gòu)的冰激振動問題，冰振進(jìn)而會引起平臺主體結(jié)構(gòu)以及上部設(shè)備與管線的疲勞破壞，同時會影響平臺上部工作人員的舒適度，降低工作效率[2]。寒區(qū)抗冰導(dǎo)管架平臺在服役期間受到海冰的長期作業(yè)，動冰荷載會使結(jié)構(gòu)產(chǎn)生疲勞損傷問題。同時波浪荷載、風(fēng)荷載、環(huán)境的腐蝕等其他環(huán)境荷載長期積累也會引起結(jié)構(gòu)不同程度的損傷問題。這種損傷往往是隱性的，需要定期對結(jié)構(gòu)進(jìn)行檢測，以保證結(jié)構(gòu)的安全。傳統(tǒng)的損傷檢測方法主要有聲波法、超聲波法、磁場檢測法、X 射線檢測法和電渦流檢測法等[3]。但導(dǎo)管架平臺是一個復(fù)雜的大尺度的結(jié)構(gòu)物，利用傳統(tǒng)的檢測方法不但效率低且實時性差。當(dāng)結(jié)構(gòu)發(fā)生損傷時，結(jié)構(gòu)的剛度、強(qiáng)度、阻尼等結(jié)構(gòu)特征參數(shù)會發(fā)生改變，從而使結(jié)構(gòu)響應(yīng)的時域、頻域、時頻域特征發(fā)生改變，所以可以通過監(jiān)測到的結(jié)構(gòu)冰激振動響應(yīng)實現(xiàn)對結(jié)構(gòu)的損傷識別、定量和定位。冰荷載引起的結(jié)構(gòu)振動響應(yīng)相對于風(fēng)、波浪等常規(guī)載荷要大一個甚至幾個量級，是十分容易監(jiān)測到的，因此更易通過響應(yīng)進(jìn)行損傷識別。這種方法相對于傳統(tǒng)的結(jié)構(gòu)檢測方法效率高，且可以實現(xiàn)對損傷的長期實時監(jiān)測。

海冰與導(dǎo)管架平臺作用是一個隨機(jī)過程[4]，得到的結(jié)構(gòu)響應(yīng)是一個非線性非平穩(wěn)的信號。常見的振動響應(yīng)信號處理方法包括傅里葉變換和小波變換，但是這兩種信號處理方法都具有一定的局限性。傅里葉變換是把信號展開為一系列三角函數(shù)的和，代表信號的全局性質(zhì)，難以捕捉突變的信號。小波變換的信號分析結(jié)果依賴于小波基函數(shù)的選取。HHT（Hilbert-Huang Transform）是一種針對非線性非平穩(wěn)信號的處理方法[5]，其可以根據(jù)信號的特點自適應(yīng)地產(chǎn)生固有模態(tài)函數(shù)，并不受測不準(zhǔn)原理制約，適合突變信號，可以較好地分析信號的瞬時特性，被廣泛應(yīng)用于地震信號分析[6]、海洋結(jié)構(gòu)損傷檢測[7]中。

本文應(yīng)用HHT 方法對結(jié)構(gòu)振動響應(yīng)信號進(jìn)行分析和處理，通過結(jié)構(gòu)損傷前后Hilbert邊際能量譜的變化構(gòu)建損傷指標(biāo)，利用該指標(biāo)對JZ20-2MUQ 平臺模型進(jìn)行損傷識別，通過對有限元仿真得到的響應(yīng)數(shù)據(jù)的分析處理，實現(xiàn)冰荷載作用下的抗冰導(dǎo)管架平臺模型的損傷識別，并探討海冰的冰速、冰厚對該損傷指標(biāo)的影響，探究利用HHT對抗冰導(dǎo)管架平臺進(jìn)行損傷識別的可行性。

1 基于Hilbert-Huang 變換的響應(yīng)信號的損傷特征量的構(gòu)造

1.1 經(jīng)驗?zāi)B(tài)分解（EMD）

Hilbert-Huang 變換包含兩部分，第一部分為經(jīng)驗?zāi)B(tài)分解（EMD），它是由Huang 提出的；第二部分為Hilbert譜分析（HAS）。HHT處理非線性非平穩(wěn)信號的基本過程是：首先利用EMD 方法將給定的信號分解為若干固有模態(tài)函（IMF），這些IMF是滿足一定條件的分量；然后，對每一個IMF進(jìn)行Hilbert變換，并表示在時頻域中，聯(lián)合這些變換結(jié)果得到信號的時頻特征；最后，利用信號的時頻特征進(jìn)行損傷特征量的構(gòu)造。

經(jīng)驗?zāi)B(tài)分解（EMD）方法的實質(zhì)是通過特征時間尺度來識別信號中所內(nèi)含的所有振動模態(tài)。其分解出的固有模態(tài)函數(shù)（IMF）是根據(jù)信號本身特性決定的。EMD 分解方法是基于以下假設(shè)條件：數(shù)據(jù)至少有兩個極值，一個最大值和一個最小值；數(shù)據(jù)的局部時域特性是由極值點間的時間尺度唯一確定；如果數(shù)據(jù)沒有極值點但有拐點，則可以通過對數(shù)據(jù)微分一次或多次求得極值，然后再通過積分來獲得分解結(jié)果。

EMD 方法的分解過程其實是一個數(shù)據(jù)選擇的過程，這個數(shù)據(jù)選擇過程依據(jù)信號特點自適應(yīng)地把一個復(fù)雜信號分解為一列固有模態(tài)函數(shù)（IMF）。它滿足如下兩個條件：信號極值點的數(shù)量與零點數(shù)相等或相差是一；信號由極大值定義的上包絡(luò)和由極小值定義的下包絡(luò)的局部均值為零。

其數(shù)據(jù)選擇過程如下：

1.2 Hilbert變換

對信號經(jīng)過經(jīng)驗?zāi)B(tài)分解（EMD）得到的每一階固有模態(tài)函數(shù)（IMF）ci(t)進(jìn)行Hilbert變換，對信號進(jìn)行卷積，即

基于此變換，建立解析函數(shù)zi，并得到信號的瞬時幅值ai和瞬時相位函數(shù)φi：

瞬時頻率ωi定義為瞬時相位φi對時間t的導(dǎo)數(shù)：

信號x(t)被表示為

在這里忽略了原始信號的殘差部分rn，Re代表級數(shù)展開的實部部分，把上述公式的右邊部分叫做Hilbert譜，即

將式（9）與時間積分得到的希爾伯特邊際譜為

將H(w,t)的平方對時間積分得到的希爾伯特能量邊際譜為

希爾伯特邊際譜代表對應(yīng)信號的瞬時頻率波動的幅值分布，希爾伯特能量邊際譜代表對應(yīng)信號的瞬時頻率波動的能量分布，兩種譜均具有很高的時頻分辨率，對損傷前后信號的變化十分敏感。

1.3 基于響應(yīng)信號的損傷特征量的構(gòu)造

當(dāng)結(jié)構(gòu)產(chǎn)生損傷時，結(jié)構(gòu)的特征參數(shù)（剛度、阻尼）會發(fā)生改變，同時這些改變會影響結(jié)構(gòu)的響應(yīng)，海洋結(jié)構(gòu)的響應(yīng)信息主要包括位移、傾角、振動加速度等?？贡鶎?dǎo)管架平臺的樁腿在振動時可以看作一個多自由度系統(tǒng)，當(dāng)樁腿發(fā)生微小損傷時，特別是在海冰荷載這種較大的激勵時，其損傷位置周圍測點的結(jié)構(gòu)響應(yīng)會比平時劇烈，這時的結(jié)構(gòu)響應(yīng)信息肯定會比平時更有利于損傷的檢測。所以我們可以利用結(jié)構(gòu)的響應(yīng)信息對抗冰導(dǎo)管架平臺進(jìn)行損傷識別，利用Hilbert-Huang 變換對結(jié)構(gòu)損傷前后的希爾伯特能量邊際譜進(jìn)行比較。相比于邊際譜中的信號振幅，希爾伯特能量邊際譜可以將因結(jié)構(gòu)損傷產(chǎn)生的信號能量變化放大。對于像海洋平臺這種大尺度結(jié)構(gòu)，有利于對結(jié)構(gòu)很微小損傷的識別。因此基于希爾伯特能量邊際譜，定義信號的能量變化DI為

式中，Eh(ωi)和Ed(ωi)分別為健康的海洋結(jié)構(gòu)和損傷的海洋結(jié)構(gòu)的信號在各時刻的希爾伯特能量邊際譜值。通過式（12）可以看出，當(dāng)DI值為0 時，結(jié)構(gòu)響應(yīng)信號希爾伯特能量邊際譜前后沒有變化，表明結(jié)構(gòu)沒有損傷；而當(dāng)DI值大于0時，說明結(jié)構(gòu)響應(yīng)信號在各頻帶上會發(fā)生變化，DI值越大，說明和健康結(jié)構(gòu)的響應(yīng)信號能量差異越大，說明結(jié)構(gòu)損傷越大。因此，利用響應(yīng)信號的希爾伯特能量邊際譜變化量可以判斷結(jié)構(gòu)是否有損傷發(fā)生。

2 損傷指標(biāo)的仿真分析

2.1 結(jié)構(gòu)參數(shù)

本文選取的平臺原型是位于渤海海域遼東灣北部的JZ20-2MUQ 平臺，其為4 樁腿加錐的抗冰導(dǎo)管架結(jié)構(gòu)。該平臺所處海域設(shè)計水深為15.6 m，平臺水上部分高為39 m，水下部分為17 m，4 個樁腿間距為14 m，樁腿本身有5°左右斜度，樁腿直徑約1.5 m。平臺整體質(zhì)量為2 700 t，其中自重約1 000 t，上部結(jié)構(gòu)等效質(zhì)量約1 700 t。樁腿水面位置裝有抗冰錐體，錐體斜面為60°，錐體交界處的最大錐徑為4 m。本文通過ANSYS 建立有限元模型，平臺為鋼結(jié)構(gòu)，楊氏模量為200 GPa，材料密度為7 800 kg/m3，泊松比為0.3，利用MASS21單元模擬上部質(zhì)量，利用Pipe20 單元模擬樁腿的水上和水下部分，利用Beam188單元來模擬導(dǎo)管架中的斜撐部分，導(dǎo)管架底部約束采用6 倍樁徑法[8]，不考慮樁土的相互作用，結(jié)構(gòu)阻尼比為0.025?？贡鶎?dǎo)管架平臺有限元模型如圖1所示。

圖1 抗冰導(dǎo)管架平臺的有限元模型Fig.1 Finite element model of ice-resistant jacket platform

2.2 結(jié)構(gòu)冰荷載

寒區(qū)海上油氣開發(fā)過程中，海冰對油氣開發(fā)設(shè)施的破壞能力要遠(yuǎn)遠(yuǎn)高于其它常規(guī)環(huán)境荷載。海冰與錐體的作用是一個周期性過程，由于海冰與錐體結(jié)構(gòu)作用的特殊性，海冰會在加載一段時間后斷裂，會發(fā)生加載和卸載循環(huán)的狀態(tài)。針對錐體結(jié)構(gòu)的冰荷載，可以采用岳前進(jìn)等[9]根據(jù)大量現(xiàn)場實測數(shù)據(jù)得到的錐體的確定性冰力函數(shù)來計算，如圖2所示。用公式可以表達(dá)為

圖2 確定性冰力函數(shù)Fig.2 Determined ice force function

式中：F0為錐體結(jié)構(gòu)極值靜冰力，可近似等于錐體水平冰力FH；T為冰力周期，τ取冰力周期T的1/3；LB為冰板的破碎長度，一般取7倍的冰厚[10]；v為海冰速度，單位為m/s。

針對錐體結(jié)構(gòu)的極值靜冰力計算，我們采用Hirayama-Obara 冰力計算模型[11]的彎曲破壞項對現(xiàn)場試驗得到的冰力數(shù)據(jù)進(jìn)行修正[12]。選取這一公式的主要原因是在Hirayama-Obara 的室內(nèi)試驗中，結(jié)構(gòu)形式與渤海的錐體比較接近。由于我們分析的錐體較窄，冰力公式中只包括冰的彎曲破碎分量。Hirayama-Obara給出的冰力彎曲破碎分量計算公式為

式中：B為待定系數(shù)，渤海取3.7；σf為冰的彎曲強(qiáng)度，取渤海遼東灣設(shè)計值700 kPa；h是冰的厚度，單位為m；D為冰作用處錐體直徑，取4 m；E為海冰彈性模量，單位為Pa；ρw為海水密度，單位為kg/m3；g代表重力加速度，單位為m/s2，ν為海冰的泊松比。

在仿真分析中，把對應(yīng)工況下冰力時程進(jìn)行離散化處理并導(dǎo)入ANSYS 中形成插值形式的冰力時程曲線，冰力正向（沿著X方向）作用于平臺兩個樁腿的入水點處（節(jié)點8、9），由于考慮屏蔽效應(yīng)，因此只作用于兩個樁腿。

2.3 損傷指標(biāo)分析

為了研究本文構(gòu)建的損傷指標(biāo)對損傷識別的效果，本文在荷載和損傷位置不變的情況下，比較不同損傷程度時損傷指標(biāo)的變化情況。常用的模擬損傷的方法有：改變壁厚和減小楊氏模量。本文采用減小壁厚來模擬損傷，損傷位置為圖1 中的點6，是長為0.1 m 的單元，設(shè)置七種損傷情況來對損傷指標(biāo)進(jìn)行分析，詳見表1。首先對導(dǎo)管架平臺模型進(jìn)行模態(tài)分析，選取結(jié)構(gòu)一階平動和一階扭轉(zhuǎn)自振頻率計算得到結(jié)構(gòu)瑞利阻尼系數(shù)為：α=0.143 0，β=0.004 4。因為結(jié)構(gòu)的損傷相對結(jié)構(gòu)整體而言很微小，所以可以認(rèn)為結(jié)構(gòu)的瑞利阻尼系數(shù)是不變的。利用2.2 節(jié)中的理論，把冰荷載輸入到ANSYS 中，其中海冰的厚度為20 cm，速度為0.5 m/s，冰荷載正向作用于平臺兩個樁腿（沿著X方向），加載的時程為10 個冰力周期，利用ANSYS 瞬態(tài)動力學(xué)分析對結(jié)構(gòu)模型進(jìn)行動力響應(yīng)分析，由于冰荷載作用于X方向，所以結(jié)構(gòu)的主要動力響應(yīng)都發(fā)生在X方向上。在真實的結(jié)構(gòu)監(jiān)測過程中，傳感器往往布置在海洋平臺的每一層甲板上，同時為了研究損傷的定位，我們選定平臺的1、2、3、4、5、7點作為測點，上述測點測量的均是結(jié)構(gòu)在X方向上的振動加速度，采樣頻率為100 Hz，分析損傷指標(biāo)的敏感程度，以及縱向、橫向的定位能力。分析得到結(jié)構(gòu)各測點的振動加速度響應(yīng)信號后，選取平穩(wěn)數(shù)據(jù)段對其進(jìn)行Hilbert-Huang 變換，并計算各個工況下?lián)p傷指標(biāo)DI值。下面以6 點處未發(fā)生損傷和6 點處損傷35%（損傷情況7）兩種狀態(tài)下測點2 的振動加速度響應(yīng)值為例，說明對信號進(jìn)行Hilbert-Huang 變換和計算DI值的過程。

表1 損傷程度Tab.1 Damage degree

把測點2 輸出的加速度響應(yīng)時程信號導(dǎo)入MATLAB 對其進(jìn)行經(jīng)驗?zāi)B(tài)分解。圖3 和圖4 分別顯示了健康和損傷35%（損傷情況7）狀態(tài)下，測點2處原始響應(yīng)信號及其經(jīng)驗?zāi)B(tài)分解（EMD）的結(jié)果。

圖3 健康狀況下原始信號與EMD分解的結(jié)果Fig.3 Results of original signal and EMD decomposition in healthy condition

圖4 損傷情況7下原始信號與EMD分解的結(jié)果Fig.4 Results of original signal and EMD decomposition under Damage Condition 7

通過對各階IMF 分量進(jìn)行Hilbert 變換得到信號的Hilbert 譜，其頻率分辨率為0.01 Hz。對Hilbert譜在時域內(nèi)先平方再積分得到信號的Hilbert能量邊際譜，Hilbert能量邊際譜顯示了信號在0～2.5 Hz頻率范圍內(nèi)的能量分布。圖5（a）和圖5（b）分別為點6處未發(fā)生損傷和點6處損傷35%（損傷情況7）兩種狀態(tài)下測點2 的Hilbert 能量邊際譜，可以看出信號的頻率成分主要集中在1 Hz 左右，與結(jié)構(gòu)水平方向低階自振頻率相符。

圖5 兩種情況下Hilbert能量邊際譜Fig.5 Hilbert energy marginal spectra in two conditions

由式（12）計算得到在損傷情況7 下測點2 處損傷指標(biāo)DI值為0.006 6，與健康狀態(tài)相比有很大變化。對各種損傷情況下每個測點的響應(yīng)信號均按上述過程分析處理，得到DI值如圖6～7 所示。從圖中可以看出，損傷指標(biāo)與損傷位置和損傷程度成正相關(guān)關(guān)系，其表現(xiàn)出的規(guī)律如下：

圖6 同層甲板測點損傷指標(biāo)隨損傷情況變化Fig.6 Damage index of the monitoring points changing with damage degree on same deck

（1）如圖6（a）所示，在冰荷載作用下，頂層甲板隨著損傷位置的損傷程度的增加，當(dāng)損傷比較小的時候，靠近損傷位置6點的測點2損傷指標(biāo)DI值比較大，當(dāng)損傷比較大時，靠近冰荷載沖擊的測點1 的損傷指標(biāo)DI值變化更明顯。如圖6（b）所示，下層甲板測點7的損傷指標(biāo)DI值始終遠(yuǎn)小于靠近損傷位置6 點的測點4，在結(jié)構(gòu)監(jiān)測的過程中可以選擇海冰漂流方向頂層甲板布置傳感器，去評估損傷程度，選擇下層甲板布置傳感器，可以識別出損傷橫向位置。

（2）如圖7 所示，在冰荷載作用下，不同層甲板隨著損傷位置的損傷程度的增加，遠(yuǎn)離損傷位置的測點5的損傷指標(biāo)DI值始終小于靠近損傷位置6點的測點2和4，可以識別出損傷縱向位置。

圖7 不同層甲板測點損傷指標(biāo)隨損傷程度變化Fig.7 Damage index of monitoring points changing with damage degree on different decks

（3）由損傷情況下各個測點的DI值可以看出，損傷位置一定（點6）時，隨著損傷程度的增加，DI值變大，說明DI值與模型結(jié)構(gòu)的損傷程度存在著正相關(guān)關(guān)系。

2.4 冰厚對損傷指標(biāo)的影響

冰厚影響著冰荷載的極值大小，是最重要的海冰參數(shù)之一。保持冰速0.5 m/s、損傷位置（點6）、損傷程度（10%）不變，根據(jù)中國海海冰條件及應(yīng)用規(guī)定[13]，此海洋平臺所在海域最大冰厚為35 cm，所以改變冰厚分別為20 cm、25 cm、30 cm、35 cm，取10個周期計算各個測點處的損傷指標(biāo)DI值，計算結(jié)果如圖8所示。從圖中可以看出，當(dāng)冰層較厚時，由于激勵頻率降低，極值冰力的增大，上層甲板靠近海冰激勵的測點1的DI值會有較大的變化，下層甲板靠近損傷位置的測點4的DI值會有較大的變化，對于不同層甲板測點2、4、5，靠近損傷位置的測點DI值始終大于遠(yuǎn)離損傷位置的點。所以，損傷指標(biāo)在不同冰厚下均是有效的，可以明顯識別出損傷位置。

圖8 測點損傷指標(biāo)隨冰厚變化Fig.8 Damage index of monitoring points changing with ice thickness

2.5 冰速對損傷指標(biāo)的影響

冰速影響著冰荷載的頻率，是最重要的海冰參數(shù)之一。保持冰厚20 cm、損傷位置（點6）、損傷程度（10%）不變，根據(jù)中國海海冰條件及應(yīng)用規(guī)定[13]，此海洋平臺所在海域平均冰速為0.4 m/s、最大冰速為1 m/s，改變冰速分別為0.4 m/s、0.6 m/s、0.8 m/s 和1 m/s，取10 個周期計算各個測點處的損傷指標(biāo)DI值，計算結(jié)果如圖9所示。從圖中可以看出，隨著冰速變快，激勵頻率的加快，其他測點損傷指標(biāo)DI值影響不是很大，但是靠近損傷位置的測點4會有巨大變化，且對于不同層甲板測點2、4、5，靠近損傷位置的測點DI值始終大于遠(yuǎn)離損傷位置的點，所以損傷指標(biāo)在不同冰速下均是有效的，可以明顯識別出損傷位置。

圖9 測點損傷指標(biāo)隨冰速變化Fig.9 Damage index of monitoring points changing with ice velocity

3 結(jié) 論

本文應(yīng)用Hilbert-Huang 變換對抗冰導(dǎo)管架平臺有限元模型損傷前后的振動加速度響應(yīng)信號進(jìn)行了分析，探討了損傷程度、海冰參數(shù)對該損傷指標(biāo)的影響，得到了以下結(jié)論：

（1）當(dāng)結(jié)構(gòu)某部位發(fā)生損傷時，激勵力作用下得到的結(jié)構(gòu)各測點的振動響應(yīng)信號的模態(tài)成分都會發(fā)生變化，頂層甲板靠近海冰激勵的位置測點隨損傷程度的變化損傷指標(biāo)變化更加明顯，從而實現(xiàn)損傷程度的識別。損傷位置下層甲板橫向測點和不同層甲板測點之間距離損傷部位較近的測點DI值比較大，DI值可以表征靠近損傷位置和遠(yuǎn)離損傷部位的差異，從而實現(xiàn)損傷位置的識別。不同損傷程度、海冰厚度和海冰速度下該指標(biāo)的有效性不受影響。

（2）在同一損傷位置，隨著損傷程度的增加，損傷指標(biāo)DI值相應(yīng)變大，DI值與模型的損傷程度存在著正相關(guān)關(guān)系。

（3）損傷指標(biāo)在不同冰厚、冰速下均是有效的，可以明顯識別出損傷位置。海冰厚度的改變對損傷前后能量邊際譜的影響主要集中在損傷位置周圍和受到海冰激勵方向的上層甲板測點，會使損傷指標(biāo)有較大的波動。海冰速度對損傷位置下層甲板的DI值的影響十分明顯，當(dāng)冰速比較大時，損傷位置下層甲板的DI值遠(yuǎn)遠(yuǎn)高于其他測點的，當(dāng)冰速較大時，更有利于此方法對損傷的識別定位。

基于響應(yīng)信號Hilbert能量邊際譜構(gòu)造的損傷指標(biāo)DI值能夠?qū)崿F(xiàn)結(jié)構(gòu)損傷的定位識別，并能根據(jù)其大小對損傷嚴(yán)重程度進(jìn)行定性判斷，相比于傳統(tǒng)方法更加方便簡單，但對于損傷程度和位置的更準(zhǔn)確的識別，還需進(jìn)一步改進(jìn)優(yōu)化算法。