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

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

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

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

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

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

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

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

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

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

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

(1)

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

(2)

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

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

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

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

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

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

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

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

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

2.1 理論研究

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

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

2.2 模型試驗(yàn)

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

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

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

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

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

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

2.3 數(shù)值模擬

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

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

(3)

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

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

(4)

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

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

(5)

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

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

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

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

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

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

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

3 結(jié)論與展望

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

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

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

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

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

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

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

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

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

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