任 鵬， 周 智,， 武文華， 歐進(jìn)萍

(1.大連理工大學(xué) 土木工程學(xué)院，遼寧 大連 116024； 2.大連理工大學(xué) 海岸和近海工程國家重點實驗室，遼寧 大連 116023；3. 大連理工大學(xué) 工程力學(xué)系，工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧 大連 116024)

基于應(yīng)變響應(yīng)統(tǒng)計特征的海洋立管損傷診斷方法

任 鵬1， 周 智1,2， 武文華3， 歐進(jìn)萍2

(1.大連理工大學(xué) 土木工程學(xué)院，遼寧 大連 116024； 2.大連理工大學(xué) 海岸和近海工程國家重點實驗室，遼寧 大連 116023；3. 大連理工大學(xué) 工程力學(xué)系，工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室，遼寧 大連 116024)

考慮海洋立管的實際服役環(huán)境和現(xiàn)有監(jiān)測技術(shù)水平，發(fā)展一種基于應(yīng)變響應(yīng)統(tǒng)計特征的損傷診斷方法。推導(dǎo)以應(yīng)變測量單元為基礎(chǔ)的細(xì)長梁式結(jié)構(gòu)應(yīng)變響應(yīng)與結(jié)構(gòu)模態(tài)參數(shù)的映射關(guān)系，提出連續(xù)測試數(shù)據(jù)的均方根應(yīng)變作為海洋立管的損傷特征參量。針對損傷程度較小或量測噪聲影響較大可能引起的損傷誤判，根據(jù)小波多分辨率分析理論，利用小波去噪和分解處理該損傷特征參量的空間域數(shù)據(jù)，進(jìn)而實現(xiàn)準(zhǔn)確地?fù)p傷定位。在此基礎(chǔ)上，進(jìn)行了某張力腿式平臺頂部張緊式立管的數(shù)值模擬；依據(jù)損傷診斷策略對假定布設(shè)應(yīng)變傳感器的上部500個應(yīng)變測量單元分析得到上述損傷特征參量。結(jié)果表明：該方法對立管局部剛度下降進(jìn)行了有效的識別，對單、多處損傷的損傷判別和定位的效果良好，損傷診斷效果受量測噪聲以及服役工況和頂張力等因素變化的影響較小，可為相關(guān)工程結(jié)構(gòu)健康監(jiān)測提供理論和技術(shù)支持。

海洋立管；損傷診斷；應(yīng)變響應(yīng)；統(tǒng)計特征；小波多分辨率分析

海洋立管是海洋油氣鉆探和輸送的管線，是海洋平臺結(jié)構(gòu)的關(guān)鍵功能設(shè)施。出于減少失效風(fēng)險和增加生產(chǎn)效益的目的，國際各大油氣公司均非常重視對海洋立管、尤其是惡劣服役環(huán)境下的深海立管進(jìn)行全壽命服役期間的完整性管理。結(jié)構(gòu)健康監(jiān)測技術(shù)為立管完整性管理實時提供荷載、響應(yīng)和損傷的原位測試數(shù)據(jù)，為立管生產(chǎn)作業(yè)和風(fēng)險分析提供及時反饋[1-2]。近年來，由于經(jīng)濟(jì)性與可靠性的提升，從英國健康與安全管理部門出臺柔性立管完整性監(jiān)測指導(dǎo)[3]，到歐美各大石油公司在墨西哥灣、北海、西非以及巴西臨近海域進(jìn)行的大量實踐，結(jié)構(gòu)健康監(jiān)測技術(shù)已在立管系統(tǒng)設(shè)計驗證和安全預(yù)警方面取得了應(yīng)用。其中，應(yīng)變傳感技術(shù)的(準(zhǔn))分布式、各測點同步、在線連續(xù)測試以及直接關(guān)聯(lián)疲勞累積等優(yōu)勢，使其逐漸成為立管監(jiān)測的首選：圖1(a)為基于線性差動位移計(Linear Variable Differential Transformer,LVDT)的西非Akpo區(qū)塊立管應(yīng)變監(jiān)測系統(tǒng)[4]；圖1(b)為基于準(zhǔn)分布式光纖光柵(Fiber Bragg Grating,FBG)的墨西哥灣Tahiti區(qū)塊應(yīng)變/曲率監(jiān)測系統(tǒng)[5]。然而，作為健康監(jiān)測的目的之一，實現(xiàn)海洋立管與陸上結(jié)構(gòu)目標(biāo)一致的損傷判別與定位還比較困難，難以對其進(jìn)行全壽命的安全評估和壽命預(yù)測[6]。缺乏面向?qū)嶋H服役環(huán)境以及匹配現(xiàn)有監(jiān)測技術(shù)水平的損傷診斷方法是重要的原因，且相關(guān)研究報道比較少見。

圖1 海洋立管應(yīng)變監(jiān)測的工程實例Fig.1 Engineering examples of strain monitoring for risers

應(yīng)變是對結(jié)構(gòu)局部損傷演化十分敏感的物理量，結(jié)合應(yīng)變傳感技術(shù)的損傷診斷方法是結(jié)構(gòu)健康監(jiān)測的重點發(fā)展方向之一。GLISIC等[7]探討了基于長標(biāo)距光纖光柵和布里淵分布式光纖傳感器獲取應(yīng)變分布從而進(jìn)行損傷探測的有效性并且實施了埋地管道監(jiān)測。作為結(jié)構(gòu)振動系統(tǒng)的固有屬性而不依賴于荷載因素，與位移模態(tài)一樣，應(yīng)變模態(tài)是同一變形能平衡狀態(tài)的另一種表現(xiàn)形式，同時還可作為高敏感性的損傷特征參量。李德葆等[8-9]系統(tǒng)地闡明了應(yīng)變模態(tài)與位移模態(tài)的對應(yīng)關(guān)系以及實驗應(yīng)變模態(tài)分析方法。吳智深等[10-11]考慮長標(biāo)距光纖光柵獲取應(yīng)變響應(yīng)并推導(dǎo)出長標(biāo)距應(yīng)變模態(tài)向量及其衍生損傷指標(biāo)，發(fā)展了基于實測應(yīng)變響應(yīng)的細(xì)長梁式結(jié)構(gòu)損傷診斷方法。在油氣管道方面，提出了相關(guān)的損傷診斷方法，主要包括：管道應(yīng)變模態(tài)變化率的無損檢測法[12]、管道腐蝕檢測的應(yīng)變敏度比法[13-14]與基于應(yīng)變模態(tài)差和神經(jīng)網(wǎng)絡(luò)的管道損傷診斷[15]；盡管上述方法均表現(xiàn)出對管道結(jié)構(gòu)局部損傷的高敏感性，但尚未考慮如何以及能否通過實測響應(yīng)獲取上述損傷特征參量。

海洋立管的實際服役環(huán)境比起陸地管道更為復(fù)雜，海洋環(huán)境風(fēng)、浪和海流等形成的外部激勵通常作為平穩(wěn)窄帶隨機(jī)過程，這與一般模態(tài)參數(shù)識別所需的白噪聲激勵條件相矛盾。另外，上部平臺運(yùn)動、頂部張緊力和海床接觸等邊界條件的時變性以及立管作為細(xì)長柔性結(jié)構(gòu)的幾何非線性造成模態(tài)類損傷特征參量獲取的困難。因此，本文從現(xiàn)有監(jiān)測技術(shù)能夠獲取的應(yīng)變響應(yīng)本身出發(fā)，將連續(xù)測量數(shù)據(jù)的統(tǒng)計特征應(yīng)用在損傷診斷方法中。建立均方根應(yīng)變與結(jié)構(gòu)模態(tài)參數(shù)的映射關(guān)系，將其作為海洋立管的損傷特征參量進(jìn)行損傷的判別和定位。該方法亦能夠避免由于立管結(jié)構(gòu)模型不確定因素給損傷診斷帶來的影響。

小損傷與量測噪聲導(dǎo)致不顯著的損傷特征，文中根據(jù)小波多分辨率分析及其去噪原理，提取損傷特征參量的突變成分用以損傷定位。模擬某張力腿式平臺頂部張緊式生產(chǎn)立管的應(yīng)變響應(yīng)，分別針對兩類服役工況進(jìn)行了單處和多處損傷的診斷，驗證了本文提出方法的有效性。

1 基于應(yīng)變響應(yīng)統(tǒng)計特征的損傷診斷方法

1.1 應(yīng)變響應(yīng)的統(tǒng)計特征

對于細(xì)長梁式結(jié)構(gòu)，離散單元的劃分一般與應(yīng)變傳感器的測量標(biāo)距相對應(yīng)。文中稱此類依照幾何標(biāo)距規(guī)定其尺度的單元為應(yīng)變測量單元。

根據(jù)模態(tài)疊加原理，n自由度彈性結(jié)構(gòu)節(jié)點位移u

(1)

由應(yīng)變與位移的幾何關(guān)系，得到

(2)

由模態(tài)解耦后的頻域分析可知

(3)

模態(tài)解耦后的時域分析可得杜哈梅爾卷積積分

(4)

(5)

將式(4)代入式(5)，得到

(6)

式(6)的具體推導(dǎo)過程參考文獻(xiàn)[16]，式中RPrs(τ)為廣義力的自、互相關(guān)函數(shù)，θ1和θ2為虛時間變量。

上述平穩(wěn)隨機(jī)過程可實現(xiàn)兩對傅里葉變換，包括：功率譜密度Sε(ω)與自相關(guān)函數(shù)Rε(τ)、頻率響應(yīng)函數(shù)Hr(ω)與脈沖響應(yīng)函數(shù)hr(t)。將式(6)代入傅里葉變換對(詳見文獻(xiàn)[16])，得到

(7)

針對第j個應(yīng)變測量單元的響應(yīng)自譜密度，式(7)變?yōu)?/p>

(8)

利用上式可將各階模態(tài)的獨(dú)立貢獻(xiàn)和耦合影響分開

(9)

(10)

1.2 損傷特征參量

對于低阻尼體系，根據(jù)模態(tài)坐標(biāo)系統(tǒng)振型互相獨(dú)立的性質(zhì)，各階模態(tài)的響應(yīng)是統(tǒng)計無關(guān)的，耦合影響可以忽略不計，此時式(9)和式(10)等式右邊交叉項為零[16]。

如響應(yīng)為平穩(wěn)隨機(jī)過程，如果滿足以下兩項識別條件，式(9)可用于峰值提取法識別模態(tài)參數(shù)[17]：

二為隨機(jī)激勵為寬頻帶的近似白噪聲；對于激勵卓越周期處存在峰值的輸入功率譜，認(rèn)為其能量峰值不能與結(jié)構(gòu)共振范圍的峰值相比；

驗證上述條件的方法：結(jié)構(gòu)響應(yīng)點自譜密度峰值曲線，與擬識別第r階振型形狀處于同一相位或相差180°，否則該峰值曲線對應(yīng)激勵能量峰值。

(11)

式中:SPk(ω)為結(jié)構(gòu)離散作用點的海洋環(huán)境荷載譜，如荷載譜已知，可用數(shù)值積分求解應(yīng)變響應(yīng)方差。另一方面，應(yīng)變響應(yīng)方差也可通過立管應(yīng)變測量單元實測獲取[18]；則上式可仿照隱式的靜態(tài)問題，即

(12)

基于上述推導(dǎo)過程，定義應(yīng)變測量單元j無量綱均方根應(yīng)變σεj(Root-Mean-Squared Strain, RMSS)為海洋立管的結(jié)構(gòu)損傷特征參量。該損傷特征參量具有如下特性：

1)當(dāng)海洋環(huán)境荷載作為平穩(wěn)(或分段平穩(wěn))隨機(jī)過程處理時，立管振動響應(yīng)中包含較多反映結(jié)構(gòu)狀態(tài)的信息。式(11)建立的RMSS與結(jié)構(gòu)模態(tài)參數(shù)的映射關(guān)系中，應(yīng)變模態(tài)振型系數(shù)最為敏感。

2)除海洋環(huán)境荷載外，RMSS僅與某一應(yīng)變測量單元j的應(yīng)變、位移模態(tài)振型系數(shù)有關(guān)；應(yīng)變傳感器布置于易損區(qū)域即可，不用獲取全局性的模態(tài)振型。

3)上部平臺運(yùn)動、頂部張緊力和海床接觸等邊界條件改變導(dǎo)致立管軸向應(yīng)變等發(fā)生變化，可通過剔除連續(xù)實測應(yīng)變數(shù)據(jù)的統(tǒng)計平均趨勢，保留應(yīng)變響應(yīng)的脈動分量，從而忽略立管結(jié)構(gòu)模型不確定因素給損傷診斷帶來的影響。

4)根據(jù)式(11)，RMSS的獲取與無損傷狀態(tài)的結(jié)構(gòu)基準(zhǔn)模型無關(guān)，且適合在線處理。

1.3 基于小波多分辨率分析的改進(jìn)

1988年，MALLAT在構(gòu)造正交小波基的工作中提出了多分辨率(尺度)分析的概念，其框架為基于同一空間不同尺度分析信號。對滿足一定條件的空間集合{Vl,l∈Z}存在尺度函數(shù)序列αl構(gòu)成標(biāo)準(zhǔn)正交基，通過構(gòu)建正交于尺度空間的小波空間{Wl,l∈Z}形成標(biāo)準(zhǔn)正交小波基函數(shù)βl。由于空間數(shù)目為無限可數(shù)的，因此多分辨率分析使得所感興趣的信號特征即使在某一分辨率檢測不到，在另一個分辨率卻很容易觀察處理，使得小波變換成為真正的時頻局部化方法。

在MALLAT塔式算法中，使用分別等價于低通和高通濾波器的尺度和小波濾波系數(shù)，從最高分辨率的尺度空間開始，根據(jù)雙尺度方程依次構(gòu)造下一尺度的尺度函數(shù)和小波函數(shù)。該算法使得構(gòu)造出的尺度基函數(shù)具有信號平均的成分，其對應(yīng)的小波變換系數(shù)稱為近似系數(shù)Wla；小波基函數(shù)具有信號變化的成分，其對應(yīng)的小波變換系數(shù)稱為細(xì)節(jié)系數(shù)Wld。從而原始信號s經(jīng)l尺度分解為小波近似成分a和小波細(xì)節(jié)成分d

(13)

本文涉及小波多分辨率分析以下三方面的應(yīng)用

1)信號平穩(wěn)處理：當(dāng)獲取應(yīng)變響應(yīng)數(shù)據(jù)的統(tǒng)計平均變化較大時，采用小尺度范圍進(jìn)行小波分解和重構(gòu)，使得響應(yīng)信號為脈動信號而趨于平穩(wěn)。

2)信號去噪：實測應(yīng)變響應(yīng)不可避免地受到量測噪聲的影響，同時RMSS本身作為響應(yīng)統(tǒng)計特征也具有一定的波動性。對于程度較小的損傷來說，可能不具有較為明顯的突變；此時認(rèn)為損傷引起的空間域信號變化是振動系統(tǒng)物理參數(shù)變化的表征，屬于有用信號，可在某一尺度將其成分體現(xiàn)出來，所以采用分層閾值去噪方法，將無用的噪聲信號濾除，保留包含損傷的信號成分。

(14)

3)信號細(xì)節(jié)分析：對損傷特征參量RMSS處理后的空間域信號進(jìn)行小波分解，提取其細(xì)節(jié)成分，從而放大現(xiàn)有信號的突變。

1.4 損傷診斷策略

1)依據(jù)應(yīng)變傳感器布置獲取應(yīng)變測量單元的應(yīng)變響應(yīng)數(shù)據(jù)；由于海洋立管同時存在彎曲和拉伸變形，應(yīng)變傳感器一般沿管壁環(huán)向布置，即采用90°或120°夾角方向的布設(shè)方案，如圖2所示；管壁應(yīng)變還須考慮雙向應(yīng)力狀態(tài)橫向效應(yīng)的修正，具體換算方法詳見文獻(xiàn)[20-21]，不重復(fù)贅述；

圖2 應(yīng)變傳感器沿管道環(huán)向布置示意圖Fig.2 Schematic diagram of hoop-deployed strain sensors

2)提取連續(xù)測量的平穩(wěn)時間段應(yīng)變響應(yīng)；如響應(yīng)的統(tǒng)計平均變化較大，采用小波分析進(jìn)行處理，剔除統(tǒng)計平均趨勢，得到響應(yīng)的脈動分量；

3)直接計算各應(yīng)變測量單元上述應(yīng)變響應(yīng)數(shù)據(jù)的功率譜密度Sεj(ω)；

5)通過RMSS離散點連接曲線的突變處判別損傷的出現(xiàn)和確定損傷的位置；

6)當(dāng)損傷程度較小或量測噪聲影響較大造成難以人為觀察突變時，采用小波去噪的方法，即通過小波分解、對分解后的細(xì)節(jié)系數(shù)進(jìn)行閾值量化處理和重構(gòu)等步驟保留包含損傷特征的空間域信號RMSSnew；此時式(14)中相關(guān)參數(shù)l,m和λ需重點考慮；

7)對小波去噪處理后的RMSS空間域信號分解，利用其細(xì)節(jié)成分進(jìn)行損傷定位。

2 數(shù)值案例

2.1 立管模型與模擬工況

目標(biāo)立管為南海海域某張力腿式平臺(TLP)頂部張緊式生產(chǎn)立管(Top Tensioned Riser,TTR)[22]，立管模型參數(shù)如表1所示。根據(jù)典型TTR立管的設(shè)計分析結(jié)果，目標(biāo)立管上部受波浪和上部浮體運(yùn)動影響較大可視為易損區(qū)域。假設(shè)立管上部500 m布置標(biāo)距為1 m的應(yīng)變傳感器，并通過尺度為1 m的應(yīng)變測量單元獲取得到應(yīng)變響應(yīng)。

由于本文僅從損傷診斷方法的驗證方面考慮，假設(shè)波浪和海流為同向來流，應(yīng)變傳感器布置于來流方向平面內(nèi)，獲取應(yīng)變響應(yīng)ε(t)，實際工況下依據(jù)圖2所示應(yīng)變傳感器布置及其各應(yīng)變換算關(guān)系，即可得到該主方向的ε(t)。TLP平臺與TTR立管的二維簡化模型如圖3所示。

表1 立管模型參數(shù)

圖3 TLP平臺及其TTR立管模型示意圖Fig.3 Schematic diagram of the TLP platform and its TTR riser

考慮有限元數(shù)值模擬所需環(huán)境荷載資料有限，應(yīng)盡量合理地設(shè)置TTR立管的服役工況。本文針對環(huán)境荷載和頂張力改變設(shè)置兩類工況驗證提出的損傷診斷策略。立管頂部張緊力分別設(shè)為400和1 000 kN。工況一中有義波高7.0 m，譜峰周期12.1 s；工況二有義波高13.8 m，譜峰周期16.1 s。拖拽力系數(shù)取1，慣性力系數(shù)取1.5。海流環(huán)境條件近似為表2。另外，假設(shè)上部平臺運(yùn)動xTLP主要由二階慢漂力引起，并作用在平面內(nèi)的縱蕩方向，慢漂周期為150 s。

表2 兩類工況下的海流流速(m/s)

隨機(jī)波浪被視為脈動隨機(jī)過程，具有各態(tài)歷經(jīng)性。不同海況的波浪頻譜(或稱能量譜)計算得到具有統(tǒng)計特性的有義波高和跨零周期(或譜峰周期)；選擇合適的波浪理論，如線性Airy波或Stokes五階波理論等，通過求得的波浪運(yùn)動的速度勢，計算流場中水質(zhì)點的速度和加速度；在此基礎(chǔ)上，通過Morison公式可計算波浪對立管的環(huán)境荷載，如拖拽力和慣性力等。海流荷載計算方法同理。依據(jù)通用有限元軟件Abaqus提供的Aqua模塊進(jìn)行定義加載，進(jìn)而通過直接積分法模擬立管的動力響應(yīng)。

2.2 結(jié)果分析

圖4 前兩階模態(tài)振型和功率譜峰值曲線的對比Fig.4 Comparisons of the first two mode shapes and peaks of PSDs

如圖4為海洋立管服役工況一前兩階模態(tài)振型與應(yīng)變測量單元(采用Abaqus的B21梁元模擬)響應(yīng)功率譜密度前兩階峰值曲線的無量綱對比?？梢钥闯觯捎诤Ｑ蟓h(huán)境荷載為窄帶隨機(jī)過程，立管結(jié)構(gòu)響應(yīng)受到激勵能量驅(qū)動，難以識別出模態(tài)參數(shù)。因此，本文采用基于響應(yīng)統(tǒng)計特性的RMSS進(jìn)行損傷判別與定位，并通過以下?lián)p傷狀態(tài)驗證其損傷診斷效果。

圖5為根據(jù)第300號單元三種損傷狀態(tài)的100 s應(yīng)變響應(yīng)數(shù)據(jù)得到的損傷特征參量RMSS。其中，圖5(a)為服役工況一，單元剛度下降20%；圖5(b)為服役工況改變，單元剛度仍下降20%；圖5(c)為服役工況二，單元質(zhì)量增加20%。可假定單元剛度下降為管壁裂紋造成，單元質(zhì)量增加由海生物附著引起。結(jié)果說明：當(dāng)服役環(huán)境發(fā)生改變時，只需采集一定時間段的平穩(wěn)隨機(jī)過程或脈動分量，即可使用RMSS進(jìn)行損傷診斷；該損傷特征參量對剛度改變非常敏感并且不受到質(zhì)量改變的影響，上述特性可以避免剛度和質(zhì)量同時變化造成的損傷誤判。

圖5 立管單處損傷狀態(tài)的損傷特征參量Fig.5 RMSS of the riser with a single damage

圖6 立管單處損傷10%的RMSSFig.6RMSSoftheriserwithasingledamageof10%圖7 RMSS的小波去噪Fig.7Waveletde-noisingofRMSS圖8 RMSS小波去噪后的細(xì)節(jié)分解Fig.8DetaildecompositionofRMSSafterwaveletde-noising

為進(jìn)一步驗證提出方法的可行性，面向服役工況二，在第200、300和400號單元分別預(yù)設(shè)20%、10%和15%三處損傷；與此同時，對原應(yīng)變響應(yīng)數(shù)據(jù)加入10%高斯白噪聲，作為量測噪聲的模擬。圖9(a)和(b)分別為未加入和加入噪聲后損傷診斷效果，可見不同RMSS值與損傷程度變化趨勢相同，但量測噪聲對損傷特征參量影響較大。因此，根據(jù)上述損傷診斷策略中f和g兩步，進(jìn)行小波去噪和分解，得到的RMSS改進(jìn)如圖10所示。本文提出的方法能夠由量測噪聲影響下的應(yīng)變測試數(shù)據(jù)對海洋立管的多處損傷進(jìn)行定位。

圖9 立管多處損傷的RMSSFig.9RMSSoftheriserwithmultipledamages圖10 量測噪聲影響下立管多處損傷的RMSSFig.10RMSSoftheriserwithmultipledamagesunderGausswhitemeasurementnoises

3 結(jié) 論

本文面向海洋立管實際服役環(huán)境和現(xiàn)有監(jiān)測技術(shù)水平，進(jìn)行了細(xì)長梁式結(jié)構(gòu)應(yīng)變響應(yīng)統(tǒng)計特征的理論推導(dǎo)，建立了應(yīng)變響應(yīng)與結(jié)構(gòu)模態(tài)參數(shù)的映射關(guān)系，提出了應(yīng)變測量單元均方根應(yīng)變RMSS作為海洋立管的損傷特征參量；基于小波多分辨率分析對RMSS進(jìn)行了改進(jìn)，依據(jù)提出的損傷診斷策略實施例了目標(biāo)立管的數(shù)值驗證。得到主要結(jié)論如下：

(1)針對布設(shè)了應(yīng)變傳感器的立管上部500個測量單元，單、多處損傷的損傷判別和定位的效果良好，RMSS的獲取無須考慮未損傷狀態(tài)的測試數(shù)據(jù)且適合在線處理。

(2)依據(jù)損傷診斷策略對模擬的立管局部剛度下降進(jìn)行了有效的識別，RMSS僅與應(yīng)變測量單元處的應(yīng)變響應(yīng)數(shù)據(jù)相關(guān)。

(3)對損傷程度較小或量測噪聲影響較大的情況，可應(yīng)用小波多分辨率分析處理RMSS空間域數(shù)據(jù)，進(jìn)而實現(xiàn)準(zhǔn)確地?fù)p傷定位，避免損傷誤判。

本文方法或可推廣應(yīng)用到承受脈動激勵的柔性結(jié)構(gòu)，如索體、主桁架、高層建筑和輸電線塔等。

[ 1 ] COOK H H，DOUGLAS E. Riser integrity management for deepwater developments [C]//Proceedings of Offshore Technology Conference，2006，Paper No. OTC-7891.

[ 2 ] DNV. Riser Integrity Management [S]. Recommended Practice F206，2008.

[ 3 ] HSE. Guidelines for Integrity Monitoring of Unbond Flexible Pipe [S]. Offshore Technology Report，OTO-98019，1998.

[ 4 ] DOUARON S L，VUATTIER C，PERROMAT V. Akpo riser integrity monitoring system design, deployment, commissioning and start-up [C]//Proceedings of Deep Offshore Technology International Conference, New Orleans, Louisiana, 2009.

[ 5 ] KARAYAKA M, CHEN J, RUF W. Tahiti online monitoring system for steel catenary risers and flowlines [C]//Proceedings of Offshore Technology Conference，2009，Paper No. OTC-19860.

[ 6 ] GOLDSMITH B，F(xiàn)OYT E，HARIHARAN M. The role of offshore monitoring in an effective deepwater riser integrity management program [C]//Proceedings of 26th International Conference on offshore Mechanics and Arctic Engineering，San Diego，2007: 881-886.

[ 7 ] GLISIC B, HUBBELL D, SIGURDARDOTTIR D H,et al. Damage detection and characterization using long-gauge and distributed fiber optic sensors [J]. Optical Engineering, 2013, 52(8): 1-12.

[ 8 ] 李德葆，諸葛鴻程，王波. 實驗應(yīng)變模態(tài)分析原理和方法 [J]. 清華大學(xué)學(xué)報(自然科學(xué)版)，1990，30(2): 105-112. LI Debao，ZHUGE Hongcheng，WANG Bo. On the principle and technique of experimental strain modal analysis [J]. Journal of Tsinghua University(Science and Technology)，1990，30(2): 105-112.

[ 9 ] 李德葆. 實驗應(yīng)變-應(yīng)力模態(tài)分析若干問題的進(jìn)展評述 [J]. 振動與沖擊，1996，15(1): 13-11. LI Debao. A general review on several fundamental points of experimental strain/stress modal analysis [J]. Journal of Vibration and Shock，1996，15(1): 13-11.

[10] ADEWUYI A P，WU Z S，SERKER K. Strain measurements assessment of vibration-based damage identification methods using displacement and distributed [J]. Structural Health Monitoring-An International Journal, 2009, 8(6): 443-461.

[11] 吳智深，李素貞，ADEWUYI A P. 基于應(yīng)變分布響應(yīng)的模態(tài)分析理論與應(yīng)用 [J]. 科技導(dǎo)報，2010，28(8): 94-103. WU Zhishen，LI Suzhen，ADEWUYI A P. Modal analysis based on distributed strain measurements: theory and application [J]. Science and Technology Review，2010，28(8): 94-103.

[12] 任權(quán)，李洪升，郭杏林. 基于應(yīng)變模態(tài)變化率的壓力管道無損檢測[J]. 大連理工大學(xué)學(xué)報，2001，41(6): 648-652. REN Quan，LI Hongsheng，GUO Xinglin. Undamaged detection of pipe based on modal varied quotiety of strain [J]. Journal of Dalian University of Technology，2001，41(6): 648-652.

[13] 李洪升，楊智榮，李亞民,等. 壓力管道腐蝕災(zāi)變檢測的應(yīng)變敏度比法 [J]. 工程力學(xué)，2007，24(3): 165-169. LI Hongsheng，YANG Zhirong，LI Yamin，et al. Detection of the corrosion damage in pressure pipeline using srain sensitivity ratio method [J]. Engineering Mechanics，2007，24(3): 165-169.

[14] 許成祥，賈善坡，涂金釗,等. 基于工作應(yīng)變模態(tài)的管道損傷識別試驗研究 [J]. 西南交通大學(xué)學(xué)報，2013，48(6): 1031-1037. XU Chengxiang, JIA Shanpo, TU Jinzhao, et al. Experimental study on damage identification of pipeline based on operational strain modal shape [J]. Journal of Southwest Jiaotong University，2013，48(6): 1031-1037.

[15] 周邵萍，郝占峰，韓紅飛,等. 基于應(yīng)變模態(tài)差和神經(jīng)網(wǎng)絡(luò)的管道損傷識別 [J]. 振動、測試與診斷，2015，35(2): 334-338. ZHOU Shaoping，HAO Zhanfeng，HAN Hongfei, et al. Damage identification in Straigh pipeline using strain modal difference and network [J]. Journal of Vibration, Measurement, and Diagnosis，2015，35(2): 334-338.

[16] (美)克拉夫R，彭津J. 王光遠(yuǎn)(譯). 結(jié)構(gòu)動力學(xué) [M]. 北京：高等教育出版社，2006: 339-345.

[17] 寶志雯，陳志鷗. 從建筑物的脈動響應(yīng)確定其動力特性 [J]. 深圳大學(xué)學(xué)報(理工版)，1986，1(1): 36-49. BAO Zhiwen, CHEN Zhiou. Experimental determination of dynamic characteristics of a structure by means of its ambient response signals [J]. Journal of Shenzhen University(Science and Engineering), 1986，1(1): 36-49.

[18] WANG D S, XIANG W, ZHU H P. Damage identification in beam type structures based on statistical moment using a two step method [J]. Journal of Sound and Vibration, 2014, 333: 745-760.

[19] DONOHO D L，JOHNSTONE I M. Ideal spatial adaptation by wavelet shrinkage [J]. Biometrika，1994，81(3): 425-455.

[20] 許學(xué)瑞，帥健，肖偉生. 滑坡多發(fā)區(qū)管道應(yīng)變監(jiān)測應(yīng)變計安裝方法 [J]. 油氣儲運(yùn)，2010，29(10): 780-784. XU Xuerui, SHUAI Jian, XIAO Weisheng. The installation method of strain gage used to monitor pipeline strain in landslide-prone areas [J].Oil & Gas Storage and Transportation，2010，29(10): 780-784.

[21] 李艷萍，冼國棟，王維斌,等. 管道應(yīng)力應(yīng)變監(jiān)測試驗規(guī)程研究 [J]. 油氣儲運(yùn)，2008，27(10): 21-23. LI Yanping, XIAN Guodong, WANG Weibin, et al. Experimental specification research on pipeline stress and strain monitoring [J], Oil & Gas Storage and Transportation，2008，27(10): 21-23.

[22] 王世圣，謝彬，李新仲. 在南海環(huán)境條件下深水典型TLP的運(yùn)動響應(yīng)分析 [J]. 中國造船，2011，52(增刊1): 94-101. WANG Shisheng, XIE Bin, LI Xinzhong. The motion response analysis of deep water typical TLP in environment conditions of South China Sea [J]. Ship Building of China，2011，52(Sup1): 94-101.

Damagediagnosis strategy for offshore risers based on statistical features of strain responses

REN Peng1, ZHOU Zhi1,2, WU Wenhua3, OU Jinping2

(1. School of Civil Engineering, Dalian University of Technology, Dalian 116024, China；2. State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116023, China；3.State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China)

Considering the service environment and the monitoring technology level of offshore risers, a damage diagnosis method was developed based on statistical characteristics of strain responses. A mapping relation between modal parameters and strain responses of a strain-measured element was derived for a beam-like structure. A root mean square strain based on a set of test data was proposed as the damage characteristic parameter of offshore risers. Aiming at a plamage misjudging due to a lower level of damage or a larger measurement noise, the data were processed based on the wavelet de-noising and decomposition. Furthermore, a TTR on a TLP was simulated to obtain the damage characteristic parameter through analyzing 500 strain-measured elements at the top of the TTR. The results showed that the method can be used to identify the local stiffness dropping of the TTR effectively althrough the damage diagnosis is interfered by the measurement noise and other uncertain factors; the damage diagnosis strategy can provide a support for structural health monitoring both theoretically and technically.

offshore riser; damage diagnosis; strain response; statistical characteristics; wavelet multi-resolution

國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)(2011CB013705)；國家自然科學(xué)基金自由申請項目(15572072)

2015-09-07 修改稿收到日期：2016-12-31

任鵬 男，博士生，1984年12月生

周智 男，博士，教授，1973年7月生 E-mail:zhouzhi@dlut.edu.cn

TU311.3；TE58

A

10.13465/j.cnki.jvs.2017.01.004