海底管道局部沖刷研究綜述

閆 莎,劉名名,王浩丞,金瑞佳,李曦寧,王 銘

(1.中國石油大學(xué)(北京),北京 102249;2.聊城大學(xué),聊城 252000;3.成都理工大學(xué),成都 610059;4.交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所,天津 300456)

海洋油氣資源開發(fā)和海洋油氣經(jīng)濟(jì)的發(fā)展,緩解了近年來掠奪性開采陸地油氣資源而引起油氣資源枯竭這一問題[1]。因此,近年來隨著世界各國對海洋油氣資源開發(fā)的重視,海底管道運(yùn)輸?shù)陌l(fā)展產(chǎn)生了明顯飛躍。根據(jù)中海油能源經(jīng)濟(jì)研究院發(fā)布的《中國海洋能源發(fā)展報告2021》,海洋油氣生產(chǎn)已成為關(guān)鍵的能源增長極,海上油氣等化石能源將成為驅(qū)動能源發(fā)展的主要動力。海底管道作為海洋油氣的主要運(yùn)輸通道,尤為關(guān)鍵地影響海洋油氣開發(fā)作業(yè)的優(yōu)劣性,要提高海洋油氣資源的利用率,必須不斷發(fā)展海底管道相關(guān)研究工作,吸收關(guān)鍵技術(shù),穩(wěn)步提高我國海底管道技術(shù)水平[2]。

相較于陸地管道所處的環(huán)境而言,海底管道所處的環(huán)境條件更為極端。海底管道鋪設(shè)于相對復(fù)雜的海洋環(huán)境中,大多管道又需要埋設(shè)于海床面下一定深度,因此檢查維修等作業(yè)很難進(jìn)行,而部分管段又處于潮差區(qū)或波浪破碎區(qū)(特別是立管),更容易遭受波浪、潮流、冰凌等嚴(yán)重沖擊,有時甚至可能遭受海洋中漂浮的堅硬物體或船舶撞擊產(chǎn)生破壞,以及遭受拋錨破壞,從而造成海底管道的嚴(yán)重破壞,導(dǎo)致管道失效破漏,造成油氣泄露,引發(fā)海洋環(huán)境和海洋油氣經(jīng)濟(jì)的雙重危機(jī)。導(dǎo)致海底管道破壞失效的主要因素包括：海底床面運(yùn)動、波浪及水流沖刷作用[3]。海底管道在海洋水動力因素和海底床面運(yùn)動等相互耦合作用的長期影響下,管道周圍的泥沙等沉積物極易受到侵蝕和遷移,使得管道下方的海床周圍發(fā)生大規(guī)模的侵蝕,附近的泥沙等沉積物在不斷沖刷剝蝕下被掏空,管道下方出現(xiàn)懸空現(xiàn)象,從而使管道極易受到損壞,大幅降低管道的使用周期。

綜上所述,要使海洋油氣資源開發(fā)安全有效地進(jìn)行就要保證海底管道的安全和穩(wěn)定,因此開展海底管道局部沖刷研究顯得尤為重要。如何通過開展海底管道局部沖刷研究,盡可能真實地模擬和還原出海底管道局部沖刷過程,獲得管道周圍沖刷坑的形態(tài)和高度,成為保障管道安全運(yùn)行、維護(hù)海洋油氣資源和海洋環(huán)境安全的一項重要命題,對海洋油氣資源開發(fā)戰(zhàn)略具有重要的實際工程技術(shù)意義,也具有重要的科學(xué)研究價值。

1 海底管道局部沖刷過程概述

海底管道鋪設(shè)于海床面時,管道周圍的流場結(jié)構(gòu)會在來流經(jīng)過管道時瞬間產(chǎn)生變化,管道周圍的來流強(qiáng)度會顯著增加甚至在周圍出現(xiàn)旋渦。由于來流強(qiáng)度的增加以及旋渦的出現(xiàn),原本泥沙不發(fā)生運(yùn)移的海床可能開始出現(xiàn)泥沙運(yùn)移的情況,對于本就存在泥沙運(yùn)移現(xiàn)象的海床而言,海底管道的存在會進(jìn)一步加強(qiáng)泥沙運(yùn)移的強(qiáng)度,從而形成沖刷現(xiàn)象,這種沖刷稱之為海底管道局部沖刷(圖1和圖2)[4]。海底管道局部沖刷現(xiàn)象存在時,管道下方就會出現(xiàn)沖刷坑,且沖刷坑的深度隨著沖刷時間的延續(xù)而不斷增加,從而導(dǎo)致管道懸空,嚴(yán)重危害管道安全,因此,開展海底管道局部沖刷方面的研究十分必要。海底管道局部沖刷研究所關(guān)注的問題包括以下3類：泥沙起動機(jī)理研究、泥沙運(yùn)移問題研究、滲流問題[5]。

圖1 海底管道布置圖Fig.1 Schematic diagram of submarine pipeline

圖2 管道周圍沖刷坑形態(tài)Fig.2 Scour profile around the pipe

(1)泥沙起動機(jī)理研究。泥沙起動機(jī)理研究是開展海底管道局部沖刷研究的首要關(guān)鍵點(diǎn)。泥沙起動研究從本質(zhì)上講,就是研究水流作用下泥沙開始起動的條件,受泥沙顆粒自身固有的性質(zhì)、泥沙顆粒所處位置以及泥沙同水流兩者之間相互作用等多種因素影響,該問題的研究一直是該研究領(lǐng)域的前沿課題。非粘性沙起動沖刷時主要受自身有效重力和流體作用力(即切應(yīng)力和上舉力)作用,而粘性泥沙除了受到上述2個力作用外,還受到顆粒間粘結(jié)力的影響。對于非粘性泥沙,考慮泥沙啟動問題需考慮流體作用帶來的流體啟動和重力作用阻止泥沙運(yùn)動這兩者間的關(guān)系,而對粘性泥沙來說,粘性細(xì)泥沙在淤積固結(jié)條件下會形成一定的穩(wěn)定結(jié)構(gòu),顆粒間粘結(jié)力增大,成為抗拒水流沖刷力的主要因素。因此考慮粘性泥沙啟動問題時就不能忽視淤積固結(jié)條件下粘性細(xì)泥沙形成的穩(wěn)定結(jié)構(gòu)和顆粒間粘結(jié)力。迄今,科學(xué)研究和工程應(yīng)用方面主要采用下列3種方式來判定泥沙顆粒是否起動,即起動拖曳力(又稱切應(yīng)力或者剪切力準(zhǔn)則)、起動流速準(zhǔn)則、起動功率準(zhǔn)則[6]。針對非粘性泥沙的啟動,目前研究上應(yīng)用最廣泛的是時空平均臨界希爾茲剪應(yīng)力(拖曳力)理論,即河床來流作用下的剪切力(拖曳力)[7]大于泥沙起動所需要的剪切力時,泥沙就會在水流作用下開始起動,否則泥沙顆粒保持靜止,沖刷現(xiàn)象也不會發(fā)生。θcr為臨界希爾茲參數(shù),由下式定義

θcr=θcr0(cosα*+sinα*/tanφ)

(1)

式中：α*為海底坡度傾角;φ為泥沙的休止角;θcr0為平坦海床面上的臨界希爾茲參數(shù),可以由SOULSBY等提出的相關(guān)公式求解

(2)

式中：D*為無量綱的泥沙粒徑。

(2)泥沙運(yùn)移問題研究。學(xué)者通過研究提出,泥沙顆粒按照其運(yùn)動形式的不同,可以分為接觸質(zhì)、躍移質(zhì)、懸移質(zhì)及層移質(zhì)4個部分。對于均勻、非粘性泥沙來說,運(yùn)移主要可以分為2種方式,即推移和懸移[8]。泥沙啟動后,多以滾動和滑動形式運(yùn)移,成為推移質(zhì)運(yùn)輸;隨著流速的增加,部分泥沙運(yùn)移仍以推移質(zhì)為主,另一部分則開始懸浮,以懸移質(zhì)運(yùn)動。當(dāng)然也可能存在顆粒較細(xì)的海岸底部浮泥運(yùn)動和顆粒較粗的海岸底部推移質(zhì)運(yùn)動。2種運(yùn)移方式可能隨時會發(fā)生轉(zhuǎn)變。泥沙運(yùn)移形式受泥沙自身性質(zhì)影響,具有“粘性”的泥沙運(yùn)移形式和均勻、非粘性泥沙的運(yùn)移形式存在出入,較小粒徑的粘性泥沙主要偏向于懸移質(zhì)運(yùn)輸。

(3)滲流問題。海底管線上下游存在一定的流體壓差,使流體從高水頭到低水頭,因而使管線下部流體流速較大。在流體流動過程中,土體會對流體產(chǎn)生阻力作用,而流體也會給土體一種反作用力,即滲流力。滲流力將會帶走土體,致使位于管線下部的土體發(fā)生管涌現(xiàn)象[9]。實驗結(jié)果表明,對于非粘性底床來說,管涌現(xiàn)象是造成海底管道沖刷的主要原因。管涌和渦流共同破壞了管道附近床面的平衡,這標(biāo)志著沖刷的開始。

值得說明的是,上文所提出的相關(guān)公式均來自于均勻、非粘性泥沙的沖刷實驗,本文后續(xù)所建立的方程也主要是針對于均勻、非粘性泥沙底床而言。粘性泥沙顆粒之間的粘性效應(yīng)難以考慮,因此目前針對粘性泥沙底床沖刷的研究相對較少[10]。然而海底環(huán)境復(fù)雜,管道所經(jīng)過海域的絕大部分海底底質(zhì)以淤泥質(zhì)為主,研究粘性底床的海底管道沖刷特性,對真實海底管道沖刷問題預(yù)測具有重要意義。

此外,管道下方?jīng)_刷坑形態(tài)及最大平衡深度也是研究海底管道局部沖刷的重要問題。迄今,已有諸多學(xué)者對此方面開展了研究,研究發(fā)現(xiàn)平衡沖刷深度對整個沖刷過程及管道后方出現(xiàn)的旋渦脫落現(xiàn)象起著較大的作用。研究同時發(fā)現(xiàn),粘性底床的沖刷深度還會受稠度系數(shù)的影響。管道下方的沖刷深度在上游水流攜帶的泥沙顆粒量與河床面向下游輸送的泥沙顆粒量一致時達(dá)到最大深度,并將其命名為平衡沖刷深度。

基于以上對于海底管道局部沖刷不同問題的研究,有學(xué)者提出可以將海底管道局部沖刷的整個過程分為4個階段,即：(1)初始沖刷階段：由于來流作用,管線上下游存在一個壓力差,壓力差會使流體給床面一個滲流力。土體隨滲流力的增加不斷被帶走,管線下方開始出現(xiàn)管涌現(xiàn)象,從而造成懸空,這一階段稱為初始沖刷階段;(2)孔道沖刷階段：隨壓力差的不斷增加,在上一階段形成的管涌現(xiàn)象將從管道后方延伸于前方,從而在管道下方出現(xiàn)孔道?？椎捞幩髁魉俸芸?泥沙便會通過孔道不斷向管道后方?jīng)_刷,而且在管道后方逐漸集聚,形成小規(guī)模沙脊,在管道下方形成沖刷坑;(3)尾跡沖刷階段：隨著時間的推移,沖刷深度不斷加劇,當(dāng)海床與管道的間距達(dá)到某一值后,管道尾跡區(qū)就會發(fā)生尤為顯著的旋渦脫落過程。由于旋渦脫落現(xiàn)象會進(jìn)一步增加同床面的剪應(yīng)力,下游管道的沖刷現(xiàn)象將更劇烈。尾跡渦具有一定的周期性,所以觀測到的位于尾跡區(qū)內(nèi)的泥沙推移現(xiàn)象也具有周期性;(4)平衡沖刷階段：隨著沖刷的進(jìn)行,尾跡渦逐漸趨于穩(wěn)定,沖刷坑的深度和形態(tài)也趨于穩(wěn)定,床底表面沖刷形態(tài)也大致保持不變,該階段為平衡沖刷階段?，F(xiàn)階段對于海底管道局部沖刷的研究,大多集中在尾跡沖刷和平衡沖刷階段的研究(圖3)。

圖3 管道周圍沖刷過程Fig.3 Scouring process around the pipe

2 國內(nèi)外研究現(xiàn)狀

目前,世界各國的研究學(xué)者針對海底管道局部沖刷開展了很多工作,主要采用的方法大致可以分為以下3類,即理論分析、模型試驗以及數(shù)值模擬。前期階段,由于計算機(jī)發(fā)展的諸多限制,學(xué)者多通過模型試驗的方法來研究問題。關(guān)于海底管道局部沖刷的模型試驗開展了比較豐富的研究,得到了比較豐碩的成果,但模型試驗取得的結(jié)果無法與技術(shù)實踐相結(jié)合并加以驗證[11]。數(shù)值模擬法隨著計算機(jī)發(fā)展的不斷完善也開始進(jìn)一步走向成熟,能夠與理論分析、模型試驗等方法相結(jié)合,直觀有效地再現(xiàn)海底管道局部沖刷的全過程,從而可以更好地預(yù)測沖刷坑深度和形態(tài)等,促進(jìn)了關(guān)于海底管道沖刷的科學(xué)探究和實際應(yīng)用。下文將按照理論分析、模型試驗、數(shù)值模擬等3類不同方法來概述國內(nèi)外研究發(fā)展歷程及現(xiàn)狀。

2.1 理論研究

海底管道局部沖刷通常會涉及到一些物理過程,如來流運(yùn)動、管道動蕩、泥沙輸送等,以及發(fā)生于所涉及到的物理過程之間的強(qiáng)耦合作用[12]。海底管道局部沖刷是多學(xué)科交叉研究的一個前沿課題,受多種復(fù)雜系統(tǒng)所響應(yīng),關(guān)于海底管道局部沖刷的研究由來已久,下列主要以時間發(fā)展為脈絡(luò),對具有代表性的海底管道局部沖刷理論研究進(jìn)行縱向概述。

早在20世紀(jì)60年代,HERBICH等[13]對埋藏于海底的海底管道穩(wěn)定性和管道周圍存在的沖刷現(xiàn)象展開了研究,并提出新的思路預(yù)測海流作用下海底管線下方的沖刷深度。KJELDSEN[14]最先開展了針對海底管道局部沖刷問題的實驗研究,發(fā)現(xiàn)沖刷平衡深度受到管道幾何尺寸和來流速度的顯著影響,并提出了一個預(yù)測沖刷平衡深度、管徑和來流速度三者間關(guān)系的經(jīng)驗公式,從而更加豐富地完善了關(guān)于海底管道局部沖刷這一課題的理論體系。LUCASSEN[15]在KJELDSEN[14]的研究基礎(chǔ)上對該問題進(jìn)行了更深層次的研究,并進(jìn)行了豐富和優(yōu)化,提出了隧道沖刷和尾渦沖刷2個過程,并且考慮了沖刷受泥沙顆粒大小的影響。自20世紀(jì)80年代末以來,世界各國的眾多專家學(xué)者就海底管道局部沖刷這一問題開展了大量深入性的研究工作。其中,MAO[16]提出海底管道上下游之間在海洋來流作用下將會出現(xiàn)一個壓力差,這種壓力差將會致使流體對土體的滲流力產(chǎn)生,雖然這種滲流力不足以使砂體運(yùn)動,卻會加劇來流對砂體的沖刷。SUMER等[17]則發(fā)現(xiàn)了床面與震動管線在洋流作用下的相互作用關(guān)系,以及沖刷如何受管線尾部湍流的影響,并且得到了一個擬合公式可以用于沖刷平衡深度的計算。CHIEW[18]研究討論了單向定常流條件下海底管道局部沖刷的作用原理以及預(yù)防沖刷發(fā)生的方法,通過探討得到局部沖刷受來流對土體滲流作用的影響,但管道周圍存在的渦流以及由壓力差引起的管涌效應(yīng)才是導(dǎo)致局部沖刷發(fā)生的首要因素。此外,只有在管道上下游的壓力梯度大于泥沙顆粒的懸浮重力時才會發(fā)生局部沖刷。CHIEW提出一個與給定流動條件下的函數(shù)可以預(yù)測給定流量和幾何邊界海底管道的最大沖刷深度。閻通等[19]以堤北海域水下三角洲的地質(zhì)條件及現(xiàn)實工況為依據(jù),對埋置管道和裸置管道進(jìn)行研究,對比分析兩者穩(wěn)定性強(qiáng)弱,通過計算預(yù)測埋置管道沖出掩埋所需要的時間以及管道下的最大沖刷平衡深度。MYRHAUG和RUE[20]則根據(jù)SUMER所研究的用于計算最大平衡沖刷深度的經(jīng)驗公式,進(jìn)一步推導(dǎo)出了新的經(jīng)驗公式可以用于計算海底管道沖刷坑在隨機(jī)波浪作用下的寬度和深度,并且可以用于計算位于樁附近的沖刷坑深度。潘冬子等[21]進(jìn)行了海底管道局部沖刷在波載荷作用下的研究,發(fā)現(xiàn)了波載荷作用下沖刷發(fā)生的臨界條件、沖刷機(jī)理及沖刷形態(tài),將局部沖刷過程分為管道周圍產(chǎn)生旋渦階段、管道下方開口階段、管道下方掏空階段、旋渦趨于穩(wěn)定階段、沖刷平衡穩(wěn)定階段5個階段,并推導(dǎo)出了管道平衡沖刷深度與KC數(shù)及相對埋深的經(jīng)驗公式。楊兵等[22]就海洋來流、海底管道和海床三者間的耦合作用通過采用量綱分析法,展開了研究探討,從而建立海底管道周圍局部沖刷的相似準(zhǔn)則。

2.2 模型試驗

早期計算機(jī)水平有限,并且因為沖刷是一個循環(huán)往復(fù)的過程,來流所形成的沖刷坑容易被泥沙顆粒等沉積物重新填充,在工程開發(fā)上想要現(xiàn)場觀測沖刷坑的形態(tài)和最大深度比較困難,因此研究海底管道局部沖刷相關(guān)問題時現(xiàn)場試驗比較少,大部分研究局部沖刷的方法是在實驗室里面得到的。迄今為止,國內(nèi)外有許多學(xué)者通過在實驗室環(huán)境下建立物理模型,對海底管道模型進(jìn)行了一系列成功的局部沖刷試驗,通過這些物理實驗可以更加直觀有效地觀測海底管道局部沖刷的全過程。下文根據(jù)模型試驗研究側(cè)重點(diǎn)的不同以及按照模型研究時間進(jìn)展,從橫向和縱向上分別選取對海底管道局部沖刷的具有代表性的模型試驗進(jìn)行闡述。

首先對來流作用下管道局部沖刷的模型試驗研究大多是根據(jù)已建立的實驗室內(nèi)試驗?zāi)Ｐ蛠硖接懞５坠艿谰植繘_刷特性同其影響因素之間的關(guān)系,從而可以推導(dǎo)出相關(guān)經(jīng)驗公式來表示各個影響因素同最大沖刷平衡深度之間的關(guān)系。代表性人物有：KJELDSEN[14]最先采用模型試驗的方法進(jìn)行有關(guān)海底管道局部泥沙沖刷問題的研究工作,通過實驗發(fā)現(xiàn)了海底管道的沖刷深度受海底管道幾何尺寸、來流速度等因素的影響且該影響不可忽視,并且開創(chuàng)性地提出了揭示海底管道下方?jīng)_刷深度、海底管道幾何尺寸、來流速度三者之間關(guān)系的經(jīng)驗公式。LUCASSEN等[15]在KJELDSEN的實驗基礎(chǔ)上加以改進(jìn),研究發(fā)現(xiàn)了泥沙粒徑大小對沖刷過程的影響,并得出隨埋深比的增加,平衡沖刷深度將會減小的結(jié)論。IBRAHIM和NALLURI[23]進(jìn)行了試驗?zāi)M并對實驗所得結(jié)果進(jìn)行了分析探討,得出了臨界希爾茲參數(shù)對沖刷的影響,根據(jù)臨界希爾茲參數(shù)的大小,沖刷問題應(yīng)該可以分為清水沖刷(θ<θcr)和動床沖刷(θ>θcr)兩類,其中θcr為臨界希爾茲參數(shù)。在清水沖刷條件下,只有管道附近才會發(fā)生沖刷現(xiàn)象,距離管道較遠(yuǎn)處沖刷現(xiàn)象都不會存在;而在動床沖刷條件下,整個管道周圍的床面都將會發(fā)生沖刷現(xiàn)象。MAO[16]進(jìn)行了一系列研究海底管道局部沖刷的模型試驗,特別是極具代表性的水槽試驗,與海底管道局部沖刷實驗研究形成對比。CHIEW[24]則采用試驗?zāi)Ｐ偷姆椒▽瑢?dǎo)流板的管道在海底局部沖刷現(xiàn)象的相關(guān)問題展開了分析探討,研究發(fā)現(xiàn),帶有導(dǎo)流板的管道在海底局部沖刷的程度和深度將會有所增大,管道自埋發(fā)生的速度也更快。SUMER[25]就清水與動床2種不同的沖刷情況展開了試驗對比研究,分析研究結(jié)果得出,底床切應(yīng)力系數(shù)在清水沖刷條件下,對管道下方的最大沖刷深度有顯著影響;但動床沖刷條件下對其無顯著影響。隨著認(rèn)識的不斷深入,三維實驗研究也不斷加以發(fā)展。CHENG等[26]建立了三維試驗沖刷模型,通過試驗研究了管道在不同的埋置深度、來流速度和迎流角情況下對橫向沖刷的影響,研究發(fā)現(xiàn)管道橫向沖刷速度隨著埋置深度的增大而減小,隨著來流速度的增大而增大。國內(nèi)學(xué)者針對這一問題進(jìn)行了很多研究,楊兵等[27]對海底管道周圍海床沖刷問題進(jìn)行了實驗?zāi)M,并展開系統(tǒng)研究,通過實驗研究分析了局部沖刷的作用原理以及沖刷坑最大平衡深度等,并對已有的成果進(jìn)行總結(jié)和分析。吳鈺驊等[28]進(jìn)行了海底管道水槽模型試驗,研究了海床上泥沙的起動機(jī)理并闡述了間隙比對砂床流速的影響。臧志鵬等[29]開展了單向流條件下海底管道沖刷的三維研究工作,研究側(cè)重于探討海底管道縱向沖刷速率如何受不同來流速度、管道最初埋深深度和來流迎角的影響。常留紅[30]通過試驗研究在不同動力條件下,不同作用角時三維管道沖刷的影響因素以及沖刷過程,并給出了沖刷平衡經(jīng)驗公式。楊少鵬等[31]通過控制變量法進(jìn)行了單向流條件下的試驗研究,探討了平衡沖刷深度受海底管線安裝阻流器與否的影響以及受間隙比差異的影響,研究發(fā)現(xiàn),管線附近流場的分布和沖止流速是決定沖刷深度大小的關(guān)鍵所在,并得出了在安裝阻流器且存在間隙比時用于計算平衡沖刷深度的經(jīng)驗公式。

以上實驗?zāi)Ｐ椭饕轻槍蜗蛄髑闆r下進(jìn)行的,而對于鋪設(shè)于海洋環(huán)境中的海底管道,波浪對管道的作用也不可忽視。迄今,也有許多專家學(xué)者就波浪條件下海底管道沖刷的相關(guān)問題進(jìn)行了一系列模型試驗研究。CHIEW[32]通過實驗研究了海底管線上阻流板的安裝對波浪作用下管道下方的沖刷深度和程度以及沖刷速度的影響。SUMER和FREDS?E[33]采用實驗的方法探討了在波浪條件下海底管道下方的最大沖刷深度和KC數(shù)兩者間的關(guān)系。秦崇仁等[34]采用波浪水槽對不同水深、不同管徑、不同泥沙顆粒大小進(jìn)行了模型試驗,得到了波浪作用下沖刷開始發(fā)生的條件以及沖刷的廣度和深度。李玉成等[35]進(jìn)行了模型試驗研究了波浪作用條件下海底管線上的水平力、舉力等,并通過實驗分析探討了波浪作用條件下海底管線附近所處流場以及流場特性是如何影響作用于管線上的力。SUMER等[36]對海流波浪作用下的埋置于海床下的海底管道局部沖刷的臨界條件做了實驗研究,研究發(fā)現(xiàn),海底管線下方出現(xiàn)的滲流速度增加會造成管涌現(xiàn)象,而管涌現(xiàn)象正是造成局部沖刷的主要因素。NEELAMNI和RAO[37]研究了埋置于黏性土中的管線在波浪作用條件下管線附近的沖刷情況。夏令[38]采用波浪水槽物理模型試驗對不同工況下海底管道下方形成的沖刷坑深度展開了系統(tǒng)研究,研究得到了沖刷坑最大深度和KC數(shù)之間的擬合關(guān)系。潘冬子和王立忠等[39]進(jìn)行了波浪水槽物理模型試驗研究,探討了在波浪條件下海底管道產(chǎn)生局部沖刷的作用原理和床面沖刷形態(tài)以及平衡沖刷深度的大小等。ZANG等[40]假定了波浪、穩(wěn)定流作用下管道開始出現(xiàn)局部沖刷現(xiàn)象的臨界條件,并依此建立了波浪和穩(wěn)定流作用下的數(shù)值波浪水槽。孫國民[41]在基于甌江口海域的工程基礎(chǔ)上,建立了考慮波浪作用的二維潮流泥沙數(shù)學(xué)模型,計算得到了灘面最終平衡狀態(tài)下的沖刷深度。

除單向流和波浪情況外,振蕩流對海底管線沖刷的影響仍然不能忽視,一些學(xué)者也做了該方面的實驗研究。浦群等[42]利用抽氣式“U”型振蕩流水槽產(chǎn)生的周期性振蕩水流,用該振蕩流代替波浪從而研究振蕩流對砂質(zhì)、淤泥質(zhì)海床面的侵蝕,并在實驗完成后進(jìn)行海床面沖刷深度大小的測定。羊皓平[43]進(jìn)行了一系列針對海底管道沖刷的實驗研究,通過實驗探討了海底管道沖刷現(xiàn)象造成管線附近流場及水動力特征的改變,并探討了振蕩流條件下管道沖刷的機(jī)理。殷俊等[44]對原型尺寸下管道振動、泥沙顆粒運(yùn)動和底床運(yùn)動的耦合問題展開了相應(yīng)研究。

總體而言,海底管道沖刷的物理模型實驗研究往往側(cè)重于單向流情況下,而對處于波浪和振蕩流條件下的海底管道沖刷現(xiàn)象的模型實驗研究較少。研究該問題的專家學(xué)者希望可以通過建立物理模型開展實驗,從而探尋出預(yù)測管道下最大平衡沖刷深度的經(jīng)驗公式,但現(xiàn)實情況是由于實驗條件和模型選取存在差異性,導(dǎo)致目前現(xiàn)有的經(jīng)驗公式在公式表現(xiàn)形式和最終預(yù)測結(jié)果上大都存在較大出入?，F(xiàn)將學(xué)者所研究的經(jīng)驗公式整理如表1所示[45]。其中,S為管道的平衡沖刷深度;U為水流流速;g為重力加速度;D為管道直徑;d為泥沙粒徑;Ucr為平均流速;h為水深;H為波高;L為波長;θ為希爾茲數(shù);e0為管道距離初始床面的高度。

表1 沖刷深度經(jīng)驗公式表Tab.1 Empirical formula of scour depth

2.3 數(shù)值模擬

近年來,由于計算機(jī)技術(shù)的高速發(fā)展和持續(xù)優(yōu)化,計算流體力學(xué)(CFD)方法在此基礎(chǔ)上也開始有了陡增式的發(fā)展,數(shù)值模擬方法不斷優(yōu)化,在科學(xué)研究領(lǐng)域和工程開發(fā)領(lǐng)域都得到了極大應(yīng)用,也逐漸作為極為有效的工具應(yīng)用于海底管道局部沖刷領(lǐng)域的研究中,預(yù)測模擬海底管道沖刷過程從而更好預(yù)測沖刷深度等一系列關(guān)注的問題,世界各國的專家學(xué)者針對海底管道局部沖刷方面的數(shù)值模擬也做了大量研究工作。目前常用的計算流體力學(xué)軟件包括 Fluent、Simcenter STAR-CCM+、OpenFOAM、COMSOL、FLOW-3D等[46-48]。

預(yù)測管道下數(shù)值泥沙沖刷模型的方法大致分為兩類,一類為平衡泥沙沖刷模型,通過迭代得到平衡沖刷面,這種方法可以有效地預(yù)測平衡沖刷坑,但卻不能模擬海底床面從沖刷開始到平衡的演化過程;另一類為泥沙質(zhì)量守恒法,該方法同樣可以準(zhǔn)確地預(yù)測出平衡沖刷河床形態(tài),并且可以模擬床面從沖刷開始到?jīng)_刷平衡的演變過程。但是,質(zhì)量守恒方法較平衡泥沙沖刷方法而言需要更多的時間。如今,隨著計算機(jī)的不斷發(fā)展,模擬技術(shù)的計算能力顯著提高,泥沙質(zhì)量守恒方法逐漸成為主流的求解方法。因此,研究中常采用泥沙質(zhì)量守恒法來模擬床面的演變過程,控制方程如下所示

(3)

式中：yb為海床高程;p0為泥沙孔隙率大小;qT為總的泥沙輸送率,包括懸移質(zhì)泥沙輸送率和推移質(zhì)泥沙輸送率。

懸移質(zhì)泥沙輸送率由下式計算求得

(4)

式中：qs為懸移質(zhì)泥沙輸送率;ys為波浪自由面的高度;ya為參考水平高度,以下為推移質(zhì)輸沙,以上為懸移質(zhì)輸沙,文中取ya=2.0d50,其中d50為泥沙顆粒的中值粒徑;yb為海床高程;c為懸移質(zhì)泥沙的體積濃度。

推移質(zhì)泥沙輸沙率的計算采用1987年由Van Rijin提出的下列公式

(5)

式中：qb為推移質(zhì)泥沙輸送率;s為流體中的泥沙密度與流體密度的比值,定義為s=ρs/ρ,其中ρs為計算采用的泥沙密度;ρ為流體密度;其中T0為無因次過剩層剪應(yīng)力。

根據(jù)水流的不同形態(tài)采用的不同理論,又可以分為勢流沖刷模型和湍流沖刷模型。海底管道局部沖刷方面的數(shù)值模擬由來已久,早在20世紀(jì)80年代,便有部分專家學(xué)者提出將數(shù)值模擬方法應(yīng)用于海底管道局部沖刷方面的研究并開始付諸實踐。CHAO等[49]最先采用數(shù)值模擬方法對海底管道沖刷以及海床面間存在的間隙等相關(guān)問題展開了研究,利用勢流模型計算出了河床底部剪應(yīng)力以及臨界剪應(yīng)力,并通過比較兩者的大小判斷是否達(dá)到?jīng)_刷平衡。當(dāng)兩者相等時,沖刷達(dá)到平衡。HANSEN[50]在無粘假定的基礎(chǔ)之上提出了一種勢流模型,并利用該模型進(jìn)行了海底管道局部沖刷的數(shù)值模擬分析。如FREDS?E和HANSEN[51]建立了勢流模型進(jìn)行研究,采用修正VON MüLLER法建立流場模型,泥沙輸運(yùn)方程只考慮了沉積質(zhì),最終模擬結(jié)果表明該模型對管道沖刷坑中上游部分預(yù)測較好,但卻未能預(yù)測管道沖刷坑下游。BERNETTR等[52]建立勢流模型進(jìn)行數(shù)值模擬,模擬研究了波浪作用下的海底管道沖刷,并準(zhǔn)確預(yù)測出了沖刷坑深度。實際上管道中的水流主要以湍流模式存在,勢流模型無法反映真實的流體狀況,能夠精確地預(yù)測出管道上游的地形特征以及最大沖刷深度,但卻無法準(zhǔn)確預(yù)測出管道下游后方的地形特征,這是因為勢流模型無法模擬后方形成的尾渦。正如LI和CHENG[53]所提出的觀點(diǎn),勢流模型的該局限性主要是因為勢流模型不能模擬海底管道后方的尾渦脫落現(xiàn)象,而管道下游所形成的地形關(guān)鍵是因為尾渦的脫落,因此從整體上來看,勢流模型的準(zhǔn)確性差強(qiáng)人意。勢流模型開始逐漸被取代。

為克服勢流模型應(yīng)用于數(shù)值模擬時存在的局限性,湍流模型開始逐步發(fā)展,研究人員開始通過求解水流Navier-Stokes方程來構(gòu)建海底管道局部沖刷過程研究的數(shù)值模型,并采用不同方法封閉N-S方程,目前應(yīng)用廣泛的有k-ε、RNGk-ε、標(biāo)準(zhǔn)k-ω、SSTk-ω模型。早期的湍流模型是在k-ε模型基礎(chǔ)上建立的。LEEUWESTEIN等[54]建立了基于k-ε的湍流模型和泥沙輸運(yùn)方程,但發(fā)現(xiàn)通過泥沙輸運(yùn)方程計算得到的結(jié)果與實際沖刷地形存在很大出入,造成該問題的原因是在建立湍流模型和泥沙輸送方程時只考慮了推移質(zhì)的影響,而忽略了懸移質(zhì)。SUMER等[17]使用離散渦方法進(jìn)行了數(shù)值模擬對海底管道附近沖刷流場展開了研究。BEEK和WIND[55]采用k-ε封閉N-S方程進(jìn)行湍流模型下的數(shù)值模擬,并同時考慮推移質(zhì)和懸移質(zhì),數(shù)值模擬得到的結(jié)果和實驗預(yù)測結(jié)果大致吻合,但下游的模擬效果仍然不夠理想。這些研究都是以粘性流體力學(xué)理論為基礎(chǔ),并考慮了實際流體中的湍流效應(yīng),因此以上研究模型同勢流模型相比可以更加準(zhǔn)確預(yù)測管道尾流區(qū)的沖刷,相較于勢流模型的不足之處有了極大改善,但前期的湍流模型尾流區(qū)的模擬結(jié)果同實際情況仍有較大出入。據(jù)此,BR?RS[56]提出采用標(biāo)準(zhǔn)k-ε湍流模型進(jìn)行沖刷模擬,同時考慮了懸移質(zhì)、推移質(zhì)輸運(yùn)以及密度的影響,該模型針對圓柱繞流方面的模擬結(jié)果和實驗得到的結(jié)果可以進(jìn)行很好的契合,但未能模擬出管道后方尾渦的周期性脫落。李玉成等[57]提出一種新的求解方法,即采用三步有限元法和大渦模擬對海底管道周圍流場進(jìn)行離散化,并分別求得二維和三維條件作用于海底管道上的無量綱水平力和無量綱橫向力,該模型預(yù)測的結(jié)果和試驗結(jié)果可以很好地進(jìn)行契合。LI和CHENG[58]用基于大渦模擬(LES)的方法封閉N-S方程,并建立相關(guān)的數(shù)值模型開展了針對于海底管道沖刷的模擬工作,通過大渦模擬的方法可以準(zhǔn)確預(yù)測出管道后方尾流區(qū)的沖刷情況,并且得到的模擬結(jié)果可以和實驗室結(jié)果很好的吻合。該方法不僅可以成功地模擬出海底管道周圍的流場,而且相較于傳統(tǒng)的模型可以得到更加精確的床面剪切力。LIANG和CHENG[59]在學(xué)者所建立的已有沖刷模型基礎(chǔ)之上進(jìn)行了發(fā)展完善,對比分析了k-ε、k-ω、SSTk-ω等湍流模型在管道周圍沖刷流場中的模擬結(jié)果,并且基于模型進(jìn)行數(shù)值模擬,準(zhǔn)確預(yù)測了沖刷剖面。SMITH等[60]對流場進(jìn)行了數(shù)值模擬。孫婭等[61]基于Navier-Stokes方程、k-ω湍流模型以及泥沙輸運(yùn)模型對海底管道局部沖刷進(jìn)行了數(shù)值模擬,計算了穩(wěn)態(tài)邊界流情況下的海底管道局部沖刷。許文兵等[62]基于標(biāo)準(zhǔn)k-ω湍流模型進(jìn)行數(shù)值模擬,研究了往復(fù)流作用下海底管道沖刷懸空機(jī)理,并基于有限差分法進(jìn)行離散化,研究表明海底管線下方的沖刷坑在往復(fù)流情況下大致具有對稱性。文君鋒等[63]研究了沖刷深度隨雷諾數(shù)和來流流速的變化情況,發(fā)現(xiàn)沖刷深度隨著雷諾數(shù)和流速的增大而增大(圖4和圖5)。顧中浩[64]采用物理模型和計算流體力學(xué)相結(jié)合的方法,對管道周圍流場的運(yùn)動特性及沖刷坑的發(fā)展過程進(jìn)行了數(shù)值模擬,利用FLOW-3D建立海底管道局部沖刷在定床和動床情況下的研究。

圖4 相對沖刷深度隨雷諾數(shù)變化圖Fig.4 Variation of scour depth with Reynolds number

圖5 相對沖刷深度隨恒定流速變化圖Fig.5 Variation of scour depth with flow velocity

上述海底管道沖刷的數(shù)值模擬工作大多是在單向流動的情況下進(jìn)行的,由于波浪場的周期性和強(qiáng)湍流性,針對該方面的研究更為復(fù)雜、且難度更大,因此對于波浪條件下海底管道局部沖刷問題所開展的數(shù)值模擬可供借鑒資料很少。具有代表性的為LIANG和CHENG[65]利用有限差分方法求解N-S方程,建立了一個將波浪作用近似轉(zhuǎn)變?yōu)檎袷幜鞯臄?shù)值模型并進(jìn)行數(shù)值模擬,該模型忽略了管道沖刷受波浪自由表面的影響。趙明[66]和ZHAO等[67]為克服振蕩流未能考慮波浪自由表面對沖刷的影響這一局限性,便采用動網(wǎng)格技術(shù)開展了海底管線在考慮波浪自由表面情況下局部沖刷的數(shù)值模擬研究。FUHRMAN等[68]采用了類似于LIANG之前曾采用的方法,研究了振蕩流條件下海底管道線沖刷的系列問題以及海底管道沖刷受KC數(shù)的影響。劉名名等[69]以及LIU等[70]通過求解利用SSTk-ω湍流模型封閉的Navier-Stokes方程組,實現(xiàn)了對海底管道在波浪作用下產(chǎn)生的局部泥沙沖刷的模擬。LIU等[71]實現(xiàn)了流體運(yùn)動—管道振動—推移質(zhì)輸沙—懸移質(zhì)擴(kuò)散過程的全耦合動力分析,揭示了懸跨管道的振動響應(yīng)與局部沖刷之間的相互影響作用,建立了懸跨管道振動幅值和沖刷深度之間的定量關(guān)系,闡明了其中的物理作用機(jī)理。LIU等[72]提出了近岸斜坡地形上海底管道局部泥沙沖刷深度預(yù)測公式,解決了國際知名管線沖刷專家SUMER教授提出的經(jīng)驗公式,僅適用于水下平坦地形的弊端。

總體而言,隨著計算機(jī)技術(shù)的日益發(fā)展和諸多學(xué)者對海底管道沖刷這一問題深入研究,目前針對于海底管道局部沖刷研究的數(shù)值模擬方法比較全面廣泛,研究成果也比較豐碩。但整體來看,針對于該問題的數(shù)值模擬研究大都是側(cè)重于單向流條件下單管道的模擬,對于波浪以及振蕩流作用下管道周圍沖刷的研究,以及對海底背馱式管道、子母管道、海底并行布置兩管道等多管道的研究較為欠缺。

3 結(jié)論與展望

綜上所述,針對水流條件下的海底管道局部泥沙沖刷這一問題,世界各國專家學(xué)者均開展了廣泛的研究實踐,通過采用理論研究、試驗?zāi)Ｐ秃蛿?shù)值模擬等方法或多種方法相互結(jié)合的手段在該領(lǐng)域獲得了豐碩的科研成果,并促使該領(lǐng)域的研究逐漸步入更完善的階段。但是目前大多研究多集中于單向流情況下單管道附近的沖刷現(xiàn)象,針對往復(fù)流等其他情況下沖刷現(xiàn)象的研究還存在一些欠缺：

(1)針對二維管道沖刷研究較為完善,三維沖刷研究相對欠缺;

(2)單管道局部沖刷方面的數(shù)值模擬較完善,雙管道、并行管道、子母管道等多管道研究相對欠缺;

(3)單向均勻流情況下的沖刷研究較完善,渦激振動下的振動管道局部沖刷研究,波浪自由表面影響作用下的管道沖刷研究較為欠缺;

(4)研究主要針對均勻、非粘性底床的海底管道沖刷,但實際底床為淤泥質(zhì),所以研究和實際工況存在一定出入。

總之,海底管道局部沖刷方面的研究還需要學(xué)者們進(jìn)一步探索、豐富完善,對此展望如下：

(1)對神經(jīng)網(wǎng)絡(luò)等智能算法進(jìn)一步優(yōu)化,從而得到相較于經(jīng)驗公式預(yù)測更為精確的沖刷模型;

(2)對淤泥質(zhì)底床沖刷進(jìn)行進(jìn)一步研究,從而更加接近實際海底管道沖刷;

(3)對于試驗?zāi)Ｐ头矫?若能獲得更多實時數(shù)據(jù)來修正數(shù)值模型并應(yīng)用于數(shù)值模擬,將試驗?zāi)Ｐ秃蛿?shù)值模擬進(jìn)行實時有效對應(yīng),便能獲得更精確的預(yù)測模型;

(4)從數(shù)值模擬角度來看,可以完善三維沖刷的研究來進(jìn)一步提高精度,以及在模擬時考慮振動、溫度、粗糙度等多因素對管道局部沖刷的影響。