海洋立管渦激振動特性研究綜述

殷布澤，胡其會，李玉星，王武昌，朱建魯，王冬旭

（中國石油大學（華東），山東省油氣儲運安全省級重點實驗室，山東 青島 266580）

0 引 言

隨著長期大規(guī)模開采，陸地油氣資源日漸枯竭，開采難度及成本日益增加，海洋油氣資源開發(fā)已成為國際油氣生產的重要增長點。我國海洋油氣資源蘊藏豐富，近年來先后開展了荔灣3-1 和番禺34-1等南海深水油氣田開發(fā)工程項目，但深水區(qū)域油氣勘探開發(fā)仍處于早期階段，工程中仍存在許多技術難題亟待解決。

在海洋油氣開發(fā)中，海洋立管作為連接海上平臺與海底油氣生產系統(tǒng)的關鍵構件，常年承受風、浪、洋流等多流場載荷作用，不可避免地會產生渦激振動，引起管道疲勞損傷、管束碰撞、磨損等管道破壞問題。一旦發(fā)生事故，必將引起嚴重的海洋環(huán)境污染和次生災害，這也使得海洋立管的VIV研究成為海洋工程中的重要研究課題。

本文從半經驗預測模型、VIV 的多模態(tài)響應及三維特性、外部波浪作用的影響、內外流共同作用的影響、立管群的VIV 特性以及VIV 抑制等方面闡述管道VIV 的研究現(xiàn)狀，最后對海洋立管尚待研究的重點進行展望。

1 渦激振動（VIV）原理

渦激振動是一種典型的流固耦合現(xiàn)象，包含復雜的多頻多自由度振動響應。圖1 描繪了在內部和外部激勵下立管系統(tǒng)的渦激振動。簡單的渦激振動形式是由立管在平面內兩個正交自由度的彈性約束下和周圍的均勻流之間的相互作用產生的，如圖2所示。

圖1 內外流場共同激勵作用下海洋立管系統(tǒng)VIV示意圖Fig.1 Schematic for VIV of the riser system under internal and external excitation

圖2 VIV原理示意圖Fig.2 Schematic of VIV

當流體通過圓柱體時，會產生以脫落旋渦形式存在的非穩(wěn)定尾跡。漩渦交替地從圓柱的兩側脫落，產生沿橫流和順流兩個方向的周期性流體力。假設圓柱靜止不動時，圓柱受到的阻力和升力與x和y坐標軸重合，如圖3（a）所示，分別表示為

式中，F(xiàn)x、Fy分別為圓柱受到的阻力和升力，CD、CL分別為阻力系數(shù)和升力系數(shù)，ρ為流體密度，D為圓柱直徑，U為流體速度。

但實際上，圓柱受到流體力的作用會產生位移和振動。此時流體力不再與x和y坐標軸重合，而是存在一定的夾角，如圖3（b）所示。此時阻力和升力可分別表示為

圖3 圓柱在均勻流下受到流體力圖解[1]Fig.3 Diagram of fluid force under uniform flow of cylinder[1]

式中，fD為圓柱受到的平均阻力，fD'、fL分別為渦誘導產生的阻力和升力，θ是x軸與瞬時速度的夾角。

二維立管振動可以用簡化的歐拉-伯努利梁方程描述：

式中，M為質量，C為阻尼，EI為抗彎剛度，T為有效軸向拉力。

從方程的形式可以看出，立管的渦激振動是一種復雜多維的非線性、非對稱流固耦合問題，這使其理論求解、數(shù)值模擬及實驗研究均具有較大難度。

2 海洋立管VIV特性

海洋立管是連接海洋平臺與海底生產系統(tǒng)的關鍵部件，其長徑比大，且同時承受管外海洋環(huán)境荷載作用和管內流體流動作用，易誘發(fā)立管振動引起管道疲勞損傷及破壞。現(xiàn)分別從立管VIV 特性半經驗預測模型、大長徑比條件下的立管VIV 多模態(tài)響應及三維特性、外部流動以及內外流場共同作用下的立管VIV特性、立管VIV抑制等方面進行綜述。

2.1 立管VIV半經驗預測模型

立管長期振動容易產生疲勞損傷，建立完善精準的立管振動響應模型、及時準確地預測立管振動和壽命非常有必要。然而，由于立管長期處于非定常外流和內流共同作用之下，復雜的結構響應和高度的流固耦合作用對理論分析預測工作帶來了很大挑戰(zhàn)，所以立管VIV特性預測常用半經驗法。半經驗法是根據(jù)實驗獲得的立管水動力參數(shù)并與力學結構模型進行耦合求解的一種半經驗半理論方法。

尾流振子模型是最經典的半經驗模型，最早由Bishop 和Hassan（1964）[2]提出，用尾跡振蕩器來模擬尾渦的運動。在此基礎上，Hartlen 和Currie（1970）[3]提出一種用非線性VanderPol 方程來模擬彈性支撐圓柱的橫向振動以及流體力的經驗模型，但其方程系數(shù)嚴重依賴模型實驗的經驗參數(shù)。之后許多學者對其進行了改進，如Skop和Griffin（1973）[4]引入修正的VanderPol方程并與圓柱振動方程耦合，適當?shù)卣{整了方程的經驗系數(shù)，并與雷諾數(shù)范圍在400～20 000內的實驗結果進行對比，成功預測了圓柱的振幅、頻率及相位角等參數(shù)。Skop 和Balasubramanian（1997）[5]又提出了一種改進模型，該模型以局部橫向振動作為升力的分量驅動VanderPol方程，另一分量由與圓柱局部橫向速度呈線性比例的失速項控制，該改進模型具備之前模型沒有的自激性和自我限制性。Facchinetti（2004）等[6]提出了一種加速耦合形式的尾流振子模型，并探究了該模型的適用范圍。Guo（2008）等[7]在Facchinetti 提出的加速耦合尾流振子模型的基礎上，又通過功能原理考慮內流和頂張力的作用，建立了一種更為全面的立管渦激振動流固耦合預測模型。Chaplin（2005）等[8]通過將11種半經驗模型與階梯流立管渦激振動實驗數(shù)據(jù)進行對比，結果發(fā)現(xiàn)半經驗模型關于橫向振幅的預測比較準確，但關于曲率的預測存在很大誤差，說明半經驗模型的精度和求解方法還有待提高。

近年來，許多學者基于半經驗方法開發(fā)了大量的非線性VIV模型，以解決線性模型中存在的渦激振幅過高問題，并考慮多模態(tài)相互作用的影響。Maincon（2011）[9]開發(fā)了一種基于神經網絡的VIV 時域預測方法，該方法具備一定的應用前景，但欠缺穩(wěn)定性，算法和神經網絡的學習方法也有待改進；Thorsen（2017）等[10]提出了一種非線性橫流渦激振動時域分析模型，并與莫里森方程進行耦合，該模型可以考慮大位移和時變接觸條件等非線性結構效應，可以很好地預測VIV 橫流中的能量傳遞和附加質量以及阻力放大等特征；Xue（2014）等[11]基于立管強迫振動實驗數(shù)據(jù)和能量平衡原理，提出一種考慮內流和頂張力影響的順流和橫流振動耦合模型，與實驗數(shù)據(jù)對比發(fā)現(xiàn)，該模型對于振幅預測結果較為準確，而對于曲率的預測卻有較大誤差。目前大部分模型都存在這一缺點，這都是由于立管振動是多模態(tài)參與，且模態(tài)激勵作用存在重疊區(qū)域引起的，而大部分模型都忽略或簡化這一現(xiàn)象從而導致對曲率的預測存在很大誤差。

通過半經驗模型研究和實驗結果的對比發(fā)現(xiàn)，對于長徑比較大的立管，順流方向的振動對于整體VIV 的貢獻也很大，不可忽略。除此之外，內部流動尤其是段塞流對立管VIV 響應也有耦合作用，它可以降低立管系統(tǒng)的固有頻率，使之更容易產生“鎖定”現(xiàn)象加劇疲勞損傷。關于管道振動的描述多采用模態(tài)疊加原理，對各模態(tài)的貢獻做了線性處理，使得預測結果精度不高。

2.2 立管VIV多模態(tài)響應及三維特性

表1 匯總了國內外諸多立管VIV 實驗，但大部分的實驗中立管尺寸較小。實際生產中，海洋立管往往有上千米長，長徑比非常大，振動形式非常復雜，所以一般的小型實驗無法完全展現(xiàn)其高模態(tài)主導、多模態(tài)組合振動的特性。Chaplin（2005）[12]指出立管渦激振動實驗設計的參數(shù)中，雷諾數(shù)大于104量級、長徑比高于500 才能較好地保留高模態(tài)對立管振動的貢獻。而以往的大部分實驗或理論模型中的立管都是基于剛性圓柱且長徑比和雷諾數(shù)均較小的情況，所以其振動模態(tài)水平較低且單一，反映規(guī)律與實際有較大差距。

表1 VIV實驗參數(shù)Tab.1 VIV experimental parameters

Huse（1998）等[13]在挪威一個海灣進行了著名的Han?ytangen 實驗，立管的長徑比高達3000，該實驗的成功為后續(xù)研究提供了可信的數(shù)據(jù)支持；Lie 和Kaasen（2006）[14]通過分析Huse（1998）的實驗發(fā)現(xiàn)，該實驗的VIV 振動響應具有高模態(tài)主導和多模態(tài)參與的特點；Vandiver（2005）等[15]在紐約塞納卡湖做了高模態(tài)的渦激振動實驗，立管的振動模態(tài)高達24；Trim（2005）等[16]開展了長徑比高達1400的立管振動實驗，發(fā)現(xiàn)了多模態(tài)振動響應現(xiàn)象，并通過模態(tài)分析發(fā)現(xiàn)VIV 振動最大位移通常高于主導模態(tài)，隨著折算速度的增大，參與振動的模態(tài)增多，而且不同形式的外流對立管的振動模態(tài)也有影響；Gao（2016）等[17]也發(fā)現(xiàn)隨著折算速度的增大，固有頻率開始發(fā)散，即參與的模態(tài)變多，且在剪切流中行波出現(xiàn)比在均勻流中快，即剪切流下參與的模態(tài)更多；Violette（2010）等[18]使用線性穩(wěn)定方法研究了圓柱的渦激振動特性，并與Chaplin（2005）[12]的實驗數(shù)據(jù)進行對比，發(fā)現(xiàn)計算結果并不能很好地驗證實驗數(shù)據(jù)，這也證明了VIV是多模態(tài)參與，具有很大的隨機性和非線性，從而使預測難度增大。

VIV 的三維效應在實驗中很難被可視化，一般都采用數(shù)值模擬軟件進行研究。Han（2020）等[19]使用LES 方法模擬三維立管尾跡的泄渦模式發(fā)現(xiàn)，在低雷諾數(shù)下立管的三維效應較小，而在高雷諾數(shù)下，立管的三維效應較強。實際立管工況中的雷諾數(shù)往往都很高，三維效應顯著，水動力不穩(wěn)定。但由于計算資源的限制，對于大長徑比，高雷諾數(shù)的三維流固耦合很難進行直接模擬。Willden 和Graham（2001，2004）[20-21]提出了切片法理論，如圖4（a）所示。該方法假設立管VIV 在長度方向為若干帶條，且在展開方向無相關性，此種處理方法大大減少了計算資源，計算速度快，且具有較高的精度，該方法已被Meneghini（2004）[22]和Yamamoto（2004）[23]驗證。Bao（2019）等[24]考慮了展開方向的關聯(lián)效應、三維效應和計算資源、計算效率等因素，在傳統(tǒng)切片法的基礎上提出了改良的viv3D-FOAM-SJTU求解器，如圖4（b）所示，并將振動響應、水動力和尾跡模式等的模擬結果與實驗結果進行了比較。Deng（2020）等[25]應用viv3D-FOAM-SJTU 求解器結合RANS 法模擬了長徑比為481.5 的立管在均勻流下的VIV 響應的三維特性，并與實驗結果對比驗證了該方法的有效性，結果發(fā)現(xiàn)切片的厚度越厚，立管的三維效應越明顯。

圖4 兩種切片法對比圖[25]Fig.4 Comparison chart of two strip theories[25]

由此可見，立管VIV的多模態(tài)響應特性和尾跡的三維效應是目前研究的重點和難點，其作用機理的復雜性和隨機性使得預測模型的精度很難得到實質性的提高。亟需提出一個更準確的多模態(tài)參與振動的描述方法，或對現(xiàn)有的歐拉—伯努利梁方程和模態(tài)疊加原理進行修正和改進，以研究尾跡三維效應如何影響水動力的穩(wěn)定性，進而更深入地研究多因素共同作用下的VIV機理。

2.3 外部波浪作用下的立管VIV特性

目前室內立管渦激振動實驗大多只考慮簡單穩(wěn)定的外部流動，而事實上海洋立管所處環(huán)境里還包含更為復雜的風浪。平臺由于受到波浪的周期作用也會產生周期性的震蕩，如圖5 所示。立管在振蕩流作用下會產生間歇性渦激振動，過程可分為成長、鎖定、衰弱三個階段。隨著折算速度和KC數(shù)的不同，表現(xiàn)出不同的振動特性[32]：“幅值調制”、“分時特性”、“遲滯效應”和“高頻諧振”等，這也是振蕩流與穩(wěn)定來流引起立管振動的本質區(qū)別[31]。

圖5 考慮波浪作用下海洋平臺振蕩的立管VIV原理圖[32]Fig.5 Schematic of VIV under oscillating flow induced by surging of offshore platform caused by wave[32]

Bishop和Hassan（1964）[26]第一次發(fā)現(xiàn)振蕩流下圓柱受迫振動產生的渦脫落發(fā)生“遲滯現(xiàn)象”，即在某個頻率的附近，相位和幅值的變化發(fā)生了階躍，當振蕩的頻率接近渦的泄放頻率的整數(shù)倍時會出現(xiàn)升力和阻力同步的現(xiàn)象；Carberry（2005）等[27]發(fā)現(xiàn)尾渦在低頻和高頻之間轉換時，立管的升力和阻力在振幅和頻率上發(fā)生了階躍，這種變化的敏感性隨振蕩頻率和泄渦頻率的比值而變化，且是一種非線性變化；Fu（2013）[28]、Thorsen（2016）[29]等發(fā)現(xiàn)VIV 主頻在整體上與KC數(shù)呈正相關，與折算速度呈負相關，隨著KC數(shù)的增加，鎖定域趨于穩(wěn)定，隨著脫落頻率的增加，鎖定區(qū)域變得更寬；Wang（2015）[30]發(fā)現(xiàn)當約化速度較大時會出現(xiàn)明顯的“模態(tài)轉換”現(xiàn)象，當KC數(shù)較小時引起的渦激振動更為明顯，甚至出現(xiàn)了“完全鎖定”現(xiàn)象，隨著KC數(shù)的增大，響應又逐漸分化成三個階段；王俊高（2015）[33]針對平臺垂蕩誘發(fā)的懸鏈線立管渦激振動開展了實驗研究，并建立了頻率響應模型；袁昱超（2019）[32]建立了時域數(shù)值模型，并針對定常洋流、垂蕩和縱蕩之間的組合流動對懸鏈線立管VIV 響應的影響進行了分析，結果發(fā)現(xiàn)縱蕩引起的非線性時變特性強于垂蕩，組合工況下的VIV響應更為激烈。

目前關于波浪引起海洋平臺振蕩下的立管VIV 研究都簡化為定常的正弦簡單振蕩激勵，振動自由度低，且未考慮來流的多角度和隨機性。未來需要開展更符合實際工況復雜流場的實驗研究，并建立更為準確的時域預測模型以適應真實的非定常流場。

2.4 內外流動耦合作用下的立管VIV特性

海洋油氣開采中，海底管道通常采用混輸技術，管道內輸送介質包括油、氣、水甚至還有固體沙粒，屬于多相流范疇，在超長立管和水平管段組成的系統(tǒng)中很容易生成具有明顯周期性和壓力波動特性的嚴重段塞流，引起立管振動[34-35]，從而使得立管的振動規(guī)律更為復雜。由于開展考慮內流尤其是含段塞流的渦激振實驗比較困難，近年來部分學者在半經驗半理論模型中考慮了內流的影響。Blanco 和Casanova（2010）[36]在Facchinetti[6]改進模型基礎之上考慮了內流是段塞流的情況，引入了Cooper（2009）[37]的段塞流計算方法，結果發(fā)現(xiàn)段塞流可以引起立管振動模態(tài)的變化；Ortega和Rivera（2013）[38]將一種非線性VIV模型與拉格朗日段塞流追蹤模型結合，將外流處理為艾里波，分析了外部波浪對內部段塞流的影響以及內部段塞流對頂張力的影響；Liu（2013）等[39]利用哈密頓原理建立了柔性立管隨時間變化的內部流體和海洋擾動作用下的VIV響應模型，并利用Lyapunov直接法進行邊界控制，保證了立管系統(tǒng)在內外干擾下的閉環(huán)穩(wěn)定性，并通過有限差分數(shù)值仿真驗證了該模型的有效性。

然而到目前為止，很少有學者在實驗中考慮內部流體對立管振動的影響，部分學者也只是將內部流體簡化成單質穩(wěn)定流動。Lou（2007）等[40]在考慮內部為單質流和外部均勻流的情況下，測量了柔性立管的振動響應，并利用尾流振子模型進行了相應的數(shù)值模擬，結果發(fā)現(xiàn)隨著內流速的增加，振動響應幅值增大，振動頻率減小，頻率的降低可以使立管在較低的流速下產生渦激振動，而振幅的增加可以使渦激振動響應更強烈，因此應重視立管內部多相流的影響；Guo（2008）[41]和Zhang（2009）[42]發(fā)現(xiàn)在外流靜止時，內流的存在會降低立管固有頻率，在均勻外流下，內流流速越大立管振幅越大、頻率越小；Zhu（2019）等[43]研究了氣液兩相流懸鏈線立管在剪切流作用下的渦激振動規(guī)律，實驗裝置如圖6所示，結果發(fā)現(xiàn)渦激振動會減小管內液塞長度，平面內兩種振動的耦合與兩者之間的相位差和水深有關，耦合作用對段塞流的氣液比很敏感。

圖6 Zhu（2019）實驗裝置原理圖[43]Fig.6 Schematic of Zhu's experiment(2019)[43]

目前開展含內流特別是段塞流與立管渦激振動耦合的實驗極少，內流對立管振動的影響規(guī)律尚不清晰，亟需開展內部多相流尤其是段塞流與外部流場共同作用下的立管VIV 實驗，為立管VIV 理論研究、立管設計及安全生產提供數(shù)據(jù)支撐。

2.5 立管群VIV特性

隨著海洋資源開采規(guī)模的不斷增加，單根立管生產能力已不能滿足生產需要，因此多立管群逐漸被廣泛應用。相比單根立管的VIV 響應，多立管群的VIV 響應更為復雜。多立管群的組合形式包括：串列雙立管、并列雙立管、錯列雙立管（如圖7）、不同來流角度的三根立管、管群等。Sumner（2010）[44]針對均勻流中兩個“無限長”等徑圓柱不同交錯形式排列的近尾渦流型、過渡渦結構行為、雷諾數(shù)影響、升力阻力系數(shù)以及斯特勞哈爾數(shù)等進行了詳細綜述。對于串聯(lián)兩圓柱的泄渦模式，可以分為三大類（如圖8所示）：（i）加長模式、（ii）再附著模式和（iii）共同脫落模式。并列圓柱泄渦模式可分為（如圖9所示）：（i）“單體”模式、（ii）偏向流模式和（iii）平行流模式。

圖7 橫向流中兩個等直徑的圓柱體流動排列形式[44]Fig.7 Arrangement of two equal diameter cylinders under cross flow[44]

圖8 橫向流動中兩個串列的圓柱體的泄渦模式示意圖[44]Fig.8 Schematic of vortex flow patterns of two tandem cylinders under cross flow[44]

圖9 橫向流動中兩個并排的圓柱體的泄渦模式示意圖[44]Fig.9 Schematic of vortex flow patterns of two side-by-side cylinders under cross flow[44]

Hu(2008)等[45]開展了Re=7000、不同角度錯列的雙立管振動特性實驗，采用了熱線、流動可視化和粒子圖像測速技術，將不同角度、不同間距下兩圓柱的泄渦模式分為四種，并繪制了模式圖；Huera-Huarte（2011a）等[46]開展了并列雙立管渦激振動實驗，主要分析了折算速度和立管間距對系統(tǒng)振動幅值、頻率和相位差的影響；Assi（2010）等[47]研究了串列雙立管尾流誘導振動(WIV)機理，結果表明WIV并不是共振現(xiàn)象，其產生是依賴于上游立管尾渦的能量激勵，這種能量激勵只有在兩根立管位移存在大于0°小于180°的相位差時才會產生；Huera-Huarte（2011b）等[48]之后又做了關于串列立管渦流誘導振動和尾流誘導振動機理研究的實驗，并利用上游立管的橫向位移幅值和主頻率分析了下游立管的動態(tài)響應以及與下游立管瞬時橫向位移運動的同步性，得到了上游立管的動態(tài)響應與下游立管動態(tài)響應的關系；Han（2018）[49]采用RANS 模型模擬了三角陣列剛性圓柱不同攻角的渦激振動，結果發(fā)現(xiàn)來流攻角對圓柱陣列間的相互作用以及鎖定頻率范圍有很大影響。

目前立管群的研究大多是關于并列或串列雙立管，而三立管及以上管群的理論和實驗研究很少，只有少量的數(shù)值模擬研究，理論上還未能定量描述立管間干涉的作用力。后續(xù)需要進一步開展多立管群VIV實驗和理論研究，為探究不同工況下立管間的安全距離以及防止立管間碰撞提供理論支持。

2.6 VIV抑制

長期的渦激振動會使立管疲勞受損，如何防止和減少VIV 對海洋立管的損傷已成為國內外學者研究的熱點。首先在設計階段可以使結構的約化速度避開渦激振動或“鎖定”的范圍，也可以通過增加結構的質量來增大其自振頻率使其遠離渦脫落的頻率，顯然這樣會大大增加成本，而且適應性差，一旦海況變化則會失效。因此，大部分學者從改變立管周圍流場、影響渦的泄放入手，主要采用兩種方法：主動控制和被動控制。

2.6.1 主動控制

主動控制是通過計算機實時監(jiān)測立管受力和周圍流場變化，再根據(jù)PID 閉環(huán)控制理論添加干擾或直接添加開環(huán)控制性的干擾到立管或周圍的流場從而達到抑制渦激振動的目的。如抽吸與吹噴海水或氣體、聲激勵系統(tǒng)、敲擊振動等。主動控制法對技術要求高，且成本也很高，故不受大多數(shù)學者的青睞，但其具有適應性好、精確、高效的優(yōu)點。

選取本院2011年5月—2017年6月收治的術后病理確診為陰莖鱗狀細胞癌的50例患者為研究對象，患者均于術后1年內（中位時間為8個月）進行18F-FDG PET/CT檢查。中位年齡為64歲（年齡范圍：30～76歲）。臨床分期采用Jackson分期法，Ⅰ期25例，Ⅱ期16例，Ⅲ期6例，Ⅳ期3例。

Zhang（2005）等[50]將壓電陶瓷振動器嵌入結構體內部，使結構表面產生微小振動從而改變流體與結構之間的作用關系，影響渦的泄放；Shaharuddin（2017）等[51]通過MATLAB 將模糊迭代原理與PID 控制原理結合，將測量的流速作為輸入參數(shù)，經過計算得到立管的模擬橫向振幅，通過立管兩端相連的彈簧實現(xiàn)立管與系統(tǒng)的共振，進而達到抑制效果；Li（2020）[52]通過在立管上游設置水平噴氣裝置來干擾流場，抑制下游立管VIV 的產生；隨后Liu（2020）[53]開展了在立管下游不同間距放置不同噴氣級數(shù)干擾的實驗，確定了最優(yōu)抑制效果的布置間距。

其中Zhang（2005）和Shaharuddin（2017）的研究屬于閉環(huán)控制，而Li（2020）和Liu（2020）的研究屬于開環(huán)控制。開環(huán)控制具有單向控制的局限性，沒有信號反饋，不能及時調整，且具有一定的適用范圍，如果超出這個適用范圍則會加劇VIV。現(xiàn)實中的海洋環(huán)境多變，流動條件很不穩(wěn)定，所以開環(huán)控制具有一定的風險性。相反，閉環(huán)控制沒有這樣的缺點，未來的應用前景也更好。

2.6.2 被動控制

被動控制是指改變立管剖面形狀或在立管表面添加附件來擾亂流場結構、干擾漩渦的形成與泄放過程，達到抑制或削弱VIV的目的。通常VIV抑制器分為三類：

（1）影響分離線和剪切層的凸起裝置如螺旋列板等。Trim（2005）等[16]測評了不同螺距/鰭高比的螺旋列板在均勻流和剪切流動下的抑制效果，實驗立管的長徑比高達1400，呈現(xiàn)了立管的多模態(tài)組合動力響應。結果表明螺距/鰭高比較大的螺旋列板抑制效果較好，且覆蓋率在82%以上時，兩種抑制器都能100%抑制剪切流下的振動和90%均勻流下的振動。同時發(fā)現(xiàn)覆蓋率和抑制效果不是正比例關系，若覆蓋率低于82%以下，抑制效率會迅速減弱。Marcollo（2018）等[54]對一種截面為“仙人掌型”的抑制器進行了實驗評價，給出了不同的截面參數(shù)對抑制效果的分析。Li（2019）等[55]對不同幾何參數(shù)的離散螺旋列板和離散板條進行了抑制實驗，結果表明離散的螺旋列板可將渦引導至多個方向泄放，更好地破壞了渦的空間相關性和壓力場分布，而離散板條的抑制效果卻表現(xiàn)一般。

（2）影響卷吸層的裹覆裝置，如控制桿、穿孔絲網、軸向板條等。Cicolin（2014）等[56]開展了三種不同形式的網格型抑制器對VIV 抑制效果的測試，結果表明網格型抑制器的抑制效果與其布置形式有很大關系，但其本質是增加了整個系統(tǒng)的流動阻尼。Huera-Huarte（2017）[57]開展了不同參數(shù)的包裹鋼絲網對立管渦激振動抑制效果的測試實驗。結果發(fā)現(xiàn)，立管的邊界層分離和剪切層的形成發(fā)生了改變，產生了較窄的尾流。在流動方向上渦的排列更緊密，剪切層相互作用區(qū)域更加遠離立管。鋼絲網不僅修正了剪切層的形成過程，而且阻止了剪切層在圓柱兩側的相互作用。

總的來說，幾何參數(shù)合適的各種形式的抑制器對VIV 振幅抑制效果都不錯，但也存在缺點，如螺旋列板會造成拖拽力增大，增加頂張力負荷，而且海洋生物的吸附對抑制效果會大打折扣[60]。整流罩對振幅抑制效果好，還可以減弱拖拽力，但設計難度高，制作安裝成本高，且存在動態(tài)不穩(wěn)定性的問題。所以未來還需要做更多創(chuàng)新嘗試（組合兩種控制方法）來開發(fā)低成本，抑制性能優(yōu)異，適應性更強的抑制手段。

3 總結與展望

本文對近年來關于海洋立管渦激振動的研究做了較為全面的綜述，包括了之前未被深入研究的方面，如半經驗預測模型的發(fā)展、高長徑比帶來的高模態(tài)、多模態(tài)響應特性、海洋波浪引起平臺和立管系統(tǒng)的VIV 特性、立管干涉或管群之間尾流的耦合振動和內部多相流與外部流場的耦合振動。盡管有大量學者投身到VIV 的研究之中，但因為其復雜特性和存在諸多影響因素，目前關于VIV 的研究還遠遠不足，未來需要做的工作和開展的研究方向有：

（1）關于半經驗預測模型：由于做了一些簡化假設（只考慮橫流，忽略順流振動、忽略或線性處理模態(tài)之間的重疊區(qū)域等），導致大多數(shù)半經驗模型存在振幅預測偏大，曲率預測不準確等缺點。目前缺少一個更具完整性和適用性的模型，相比較理論上的忽略或簡化，可以考慮利用模糊數(shù)學和神經網絡等手段分析描述VIV的三維非線性和隨機性特征，以得到更好的預測效果。

（2）關于大長徑比立管多模態(tài)響應特性：目前大部分實驗立管的長徑比與現(xiàn)場立管相差很大，得到的結果很難應用到實際生產。亟需開展大量的大長徑比立管VIV 實驗，為深入研究大長徑比立管VIV 的多模態(tài)響應規(guī)律、三維特性以及建立更準確的預測模型提供數(shù)據(jù)基礎。在數(shù)值模擬方面，傳統(tǒng)的CFD 方法直接模擬全尺寸復雜流場下立管VIV 的計算成本極高，需要繼續(xù)深入探索其它模擬手段在立管VIV方面的應用。

（3）關于外部波浪作用的影響：有學者研究發(fā)現(xiàn)，波浪引起的振蕩流下的立管VIV 出現(xiàn)了許多定常流條件下沒有的特性，而實際只有少數(shù)預測模型和實驗考慮了外部波浪載荷對VIV 的影響。實際立管工況變化劇烈且不定常性顯著，未來的VIV實驗和預測模型的建立也應該考慮到波浪、內波等不定常流場的作用，以適應現(xiàn)場的復雜工況，得到更符合實際的結果。

（4）關于內外流共同作用的影響：目前很少有學者考慮管內流體對立管VIV的影響。有限的研究證明，管內多相流和管外流場具有相互影響作用，但二者的耦合作用機理尚不清晰，仍需要開展大量實驗和理論研究。

（5）關于多立管群的VIV 特性：由于多立管系統(tǒng)下，尾流之間的干涉情況復雜，導致定量描述困難，國內外大多數(shù)關于立管群的研究都局限在實驗研究，理論研究匱乏。而實際生產中管群的應用越來越廣泛，立管之間經常會出現(xiàn)碰撞、摩擦等破壞現(xiàn)象而影響立管的安全。因此，更深入開展多立管群的VIV特性理論研究對保證生產安全具有重要指導意義。

（6）關于VIV 抑制：主動控制存在技術難度高、成本高的限制，被動控制的VIV 抑制器也普遍存在適應性、穩(wěn)定性和附著海洋生物效率降低等問題。未來仍需開展大量工作，探究新的抑制方式，盡量減少上述缺點，在降低成本的同時達到更好的抑制效果。