高祥宇，周 偉，李書亮，段子冰，高正榮，俞竹青

(1.南京水利科學(xué)研究院 港口航道泥沙工程交通行業(yè)重點實驗室，江蘇 南京 210024; 2.長江南京航道局，江蘇 南京 210011; 3.港珠澳大橋管理局，廣東 珠海 519060; 4.廣東港珠澳大橋材料腐蝕與工程安全國家野外科學(xué)觀測研究站，廣東 珠海 519060)

在經(jīng)濟建設(shè)和可持續(xù)發(fā)展中，交通起著至關(guān)重要的作用，橋梁建設(shè)是跨越河流、海域和峽谷等天然屏障縮短交通距離的一種重要的手段。20世紀(jì)80年代以來，我國進(jìn)入了橋梁建設(shè)的高速發(fā)展期，每年開工建設(shè)橋梁有1萬多座。進(jìn)入21世紀(jì)，為適應(yīng)我國經(jīng)濟的全球化發(fā)展，促進(jìn)沿海城市、沿海城市與島嶼之間的交流與發(fā)展，跨海通道工程建設(shè)成為一種趨勢，如已建成的港珠澳大橋、杭州灣跨海大橋、青島海灣大橋、東海大橋、舟山跨海大橋、平潭跨海大橋、廈門海灣大橋等。我國大陸淤泥質(zhì)海岸長逾4 000 km，約占我國大陸海岸的22%[1]，已建的跨海大橋多數(shù)處于淤泥質(zhì)海岸。橋梁建設(shè)后橋墩周圍水流環(huán)境發(fā)生改變，會引起橋墩基礎(chǔ)沖刷。我國所建跨河跨海的大型橋梁，多數(shù)為深水群樁基礎(chǔ)，體積和阻水面積相對較大，在循環(huán)荷載和基礎(chǔ)沖刷的情況下，橋梁樁基的承載能力減小，水流的橫向沖擊力增大，由此可能引起橋梁毀壞和基礎(chǔ)裸露等嚴(yán)重安全隱患。國內(nèi)外已發(fā)生的橋梁事故統(tǒng)計表明，橋墩基礎(chǔ)沖刷是導(dǎo)致橋梁毀壞的主要原因之一[2]。橋梁建設(shè)投資一般都比較大，跨海大橋由于建設(shè)所處的地理環(huán)境特殊，建設(shè)投資更加巨大。一旦發(fā)生橋梁毀壞事故，會造成更大的經(jīng)濟損失和人員傷亡。橋梁基礎(chǔ)沖淤后為了保障橋梁安全就需進(jìn)行防護(hù)，防護(hù)經(jīng)費也巨大，東海大橋橋墩樁基沖刷防護(hù)應(yīng)急搶險工程投資經(jīng)費需1個多億。

自20世紀(jì)50年代以來，國內(nèi)外眾多學(xué)者針對橋墩沖刷進(jìn)行了大量的研究，2000年以前研究較多的是單向流沖刷計算公式，美國聯(lián)邦公路局發(fā)布的“水力工程通報”HEC-18[3]認(rèn)為橋下總沖刷由河床長期的沖淤、束窄河床沖刷、局部沖刷組成，推薦了橋墩局部沖刷公式。國內(nèi)規(guī)范[4-5]中橋墩附近沖刷包含河床自然演變沖刷、一般沖刷、局部沖刷3個部分，并按照黏性土和非黏性土給出了相應(yīng)的計算公式。美國在2000年行業(yè)規(guī)范中首次提出了橋墩局部沖刷在潮汐水流中計算方法[6]，還有一些橋墩局部沖刷深度計算公式也相繼被提出[7-11]，但這些公式由于是采用了特定的研究方法和試驗條件，其適用性仍有待進(jìn)一步考證[12]。目前，潮流橋梁局部沖刷研究或認(rèn)識也存在一定的爭議：有的學(xué)者認(rèn)為在足夠長的沖刷周期作用下，其與徑流相同流速條件下局部最大沖刷深度基本相同[13-14]，也有的認(rèn)為小于徑流作用下的局部最大沖刷深度[15-16]?？绾蛄夯A(chǔ)受到水流、波浪和潮汐等綜合作用，橋墩局部沖刷情況比單向水流和往復(fù)流沖刷更加復(fù)雜[17-18]。

港珠澳大橋是連接香港、廣東珠海和澳門的橋隧工程，位于珠江口伶仃洋海域內(nèi)。伶仃洋內(nèi)地形復(fù)雜、島嶼眾多，受地形、潮汐、徑流的綜合作用，水動力條件復(fù)雜。港珠澳大橋橋隧全長55 km，其中主橋29.6 km，設(shè)計壽命120年，項目總投資額達(dá)1 269億元。橋墩局部沖刷深度是橋梁基礎(chǔ)設(shè)計的重要參數(shù)之一，現(xiàn)有的數(shù)值計算方法和公式還無法很好地解決局部沖刷問題。下文采用系列模型，對青州航道橋主橋墩進(jìn)行水流、潮流和波流共同作用下的局部沖刷物理模型試驗，為工程設(shè)計提供依據(jù)。

1 試驗概況

1.1 橋墩結(jié)構(gòu)

港珠澳大橋有3座通航孔橋，分別為青州航道橋、江海直達(dá)船航道橋和九洲航道橋，青州航道橋跨徑最大，主塔最高的通航孔橋，雙塔之間距離458 m，其主橋墩是全橋最大的橋墩。橋墩基礎(chǔ)為樁承臺墩型，平面呈長八邊形，迎水面寬32.50 m，長53.50 m。自頂面標(biāo)高+3.8 m向下的承臺厚6.0 m。承臺底布置36根圓樁，按高程的不同分為兩種樁型：承臺底～-58.5 m之間為鋼護(hù)筒樁，標(biāo)高-58.5 m向下的為混凝土樁，樁徑為2.5 m(圖1)。

圖1 橋墩結(jié)構(gòu)型式

1.2 工程區(qū)水文和地質(zhì)條件

伶仃洋為一個東南向開敞的海域，匯集珠江入海八個口門中的虎門、蕉門、洪奇瀝和橫門四個口門的徑流，由四個口門進(jìn)入伶仃洋的平均流量約5 663 m3/s，徑流總量約1 742億m3/a，約占珠江出?？偭康?3.4%。徑流年內(nèi)分配極不均勻，汛期4月—9月約占全年徑流總量的80%。潮汐類型屬不正規(guī)半日潮混合潮型，日不等現(xiàn)象明顯，平均潮差1.06～1.16 m，屬于弱潮海灣。潮流基本為沿槽線走向的周期性往復(fù)流，落潮流速一般都大于漲潮流速，最大流速可達(dá)2.0 m/s。洪季受徑流影響，流速明顯大于枯季。內(nèi)伶仃島以內(nèi)流向以NNW～SSE向為主，內(nèi)伶仃島以外流向轉(zhuǎn)為S～N向。橋區(qū)潮流漲急時垂線平均流向大致為N向，落急時垂線平均流向大致為S向。外海大浪主要是臺風(fēng)浪，臺風(fēng)大浪的波向主要是向岸的SE～SSW向，根據(jù)橋區(qū)實測波浪資料，常浪向為偏南向，其次為偏北向，在臺風(fēng)“鸚鵡”期間最大有效波高為3.64 m，周期為5.3 s，浪向為SSW。

根據(jù)橋位沿線地質(zhì)勘探資料，橋位區(qū)下伏基巖主要為燕山期花崗巖和震旦系片麻狀混合花崗巖、混合花崗巖、片巖。基巖上部為全新統(tǒng)～更新統(tǒng)海相、陸相、河流相、海陸交互相沉積層。其中全新世海相沉積物巖性為淤泥(厚度4.7～8.6 m)、淤泥質(zhì)黏土(厚度9.3～13.3 m)和淤泥質(zhì)黏土夾砂(厚度7.4～17.4 m)；晚更新世陸相沉積物巖性主要為軟～可塑狀黏土，海相沖積物巖性主要為淤泥質(zhì)粉質(zhì)黏土、淤泥質(zhì)黏土和軟～可塑狀黏土，河流相沖積物主要由中密～密實砂類土組成。

在青州航道橋位附近取了2個原狀土泥樣，泥沙中值粒徑分別為0.023 mm和0.028 mm，根據(jù)原狀土起動試驗，起動流速分別為1.25 m/s和1.59 m/s[19]，為使沖刷試驗得出的結(jié)果偏于安全，選取泥沙中值粒徑0.023 mm作為原型泥沙代表粒徑，容重γs=2.65 t/m3。

1.3 試驗水槽和設(shè)備

模型試驗在寬水槽內(nèi)進(jìn)行，水槽總長42 m，寬4.2 m。水槽一端為往復(fù)流雙向泵，另一端為翻板尾門，并布置了造波板，水槽的外緣為回水廊道和水庫，水槽能夠產(chǎn)生單向流、潮流和波流共同作用動力條件。動床段布置在水槽中間段長度10 m，鋪沙厚度0.6 m，橋墩模型布置在中間位置(圖2)。由計算機、自動調(diào)速裝置、雙向泵頻率控制裝置、搖擺式生波裝置和數(shù)據(jù)采集箱等設(shè)備組成一套自動控制和采集系統(tǒng)。試驗中流速、水位和波浪測量分別采用二維電磁式流速儀、自動跟蹤式水位儀和CBY-Ⅱ型波高儀，橋墩外圍沖淤地形采用直讀式地形測量儀測量，橋墩區(qū)沖淤地形采用測針人工讀取。

圖2 試驗水槽布置示意

2 試驗方法和條件

2.1 系列模型延伸法

橋墩局部沖刷是一個三維問題，垂向水流運動是決定沖刷深度和形態(tài)的重要因素之一，正態(tài)動床模型預(yù)測橋墩局部沖刷形態(tài)和深度是通常采用的方法。動床模型設(shè)計時除滿足幾何相似、動力相似等條件外，模型沙的相似尤為重要。在實際工作中，由于縮尺影響，模型的動力條件遠(yuǎn)弱于原型，模型沙與原型沙達(dá)到完全相似相當(dāng)困難，為了消除模型沙和原型沙運動不完全相似而帶來的試驗結(jié)果偏差常借助系列模型延伸法。

假設(shè)完全滿足相似條件的正態(tài)模型幾何比尺為λh0，系列模型選用的不滿足泥沙起動相似的模型幾何比尺為λh,沖刷深度比尺為λhs=hp/hm(hp、hm分別為原型和模型沖刷深度)，當(dāng)模型完全滿足正態(tài)模型相似條件時，則λhs=λh=λh0；當(dāng)系列模型幾何比尺λh偏離幾何比尺λh0時，λhs≠λh≠λh0，λhs與λh存在下列指數(shù)函數(shù)關(guān)系：

λhs/λh=(λh0/λh)n

(1)

將沖刷深度比尺代入式(1)中，則原型沖刷深度hp為：

hp=hmλh(λh0/λh)n

(2)

采用系列模型試驗結(jié)果推算原型沖刷深度時，常利用雙對數(shù)坐標(biāo)系，以λh為橫軸，hm為縱軸作直線關(guān)系，推求出原型沖深。當(dāng)系列模型采用原型沙進(jìn)行試驗時，λh0=1，將λh延伸到1的hm即為原型沖淤值hp；當(dāng)模型沙與原型沙不同時，λh0等于按泥沙起動相似求得一個數(shù)值，系列模型結(jié)果應(yīng)延伸至λh=λh0，得到延伸的沖淤值，然后再與λh0相乘得到原型沖淤值。也可以將所得試驗結(jié)果進(jìn)行延伸，求得指數(shù)n，建立外插方程，直接求得原型沖刷深度。

為了較好地模擬原型動力條件，綜合考慮場地和設(shè)備等因素，本次試驗分別采用幾何比尺為50、70和100進(jìn)行系列模型試驗。

2.2 模型沙

研究橋墩附近沖刷坑深度及形態(tài)，模型沙選擇主要考慮泥沙起動相似和水下泥沙休止角的相似。采用經(jīng)過相似條件下的比尺計算和模型沙性質(zhì)的比選，采用株洲精煤研磨的煤粉作為試驗?zāi)Ｐ蜕?，試驗前進(jìn)行了脫脂處理，模型沙中值粒徑d50=0.50 mm，容重為γs=1.33 t/m3。采用張紅武泥沙水下休止角計算公式[20]得到原型沙休止角為30.3°，通過模型試驗?zāi)Ｐ蜕承葜菇菫?1.2°，模型沙與原型沙水下休止角基本相當(dāng)。通過試算模型比尺為λh0=40時模型沙與天然沙泥沙起動相似。對于波浪作用下的泥沙起動，采用周益人提出的公式[21]計算波浪泥沙起動相似比尺為λh0=67，模型沙略容易起動，偏于安全。

2.3 試驗動力條件

分別按照單純水流、洪水與潮流組合和波流共同作用三種情況進(jìn)行試驗。單純水流分別采用300年一遇和20年一遇設(shè)計水流；波流共同作用時，水流動力分別采用100年一遇和20年一遇設(shè)計水流，波浪條件采用水流對應(yīng)條件下的H4%規(guī)則波，根據(jù)某跨海大橋橋墩基礎(chǔ)的波流局部沖刷試驗[22],順流疊加波浪橋墩基礎(chǔ)沖深最大，試驗采用順流疊加波浪(表1)。洪潮組合水流動力采用數(shù)學(xué)模型計算得到的300年一遇潮位和流速過程，物理模型試驗結(jié)果與數(shù)模值誤差小于10%(圖3)，模型試驗漲急流速1.73 m/s，落急流速1.52 m/s，一個潮過程漲潮歷時10 h，落潮歷時15 h。

表1 試驗動力條件

圖3 300年一遇潮位和流速過程

2.4 模型試驗方法

沖刷試驗開始前，首先在試驗段周圍布置水位儀、流速儀和波高儀，調(diào)試水槽內(nèi)無結(jié)構(gòu)物時試驗要求的水流和波要素，將調(diào)試好的參數(shù)存入計算機控制系統(tǒng)，供沖刷試驗時采用。

進(jìn)行局部沖刷試驗時，放水前先在試驗段的中部放置橋墩，模型沙平鋪于試驗段。在試驗過程中，監(jiān)測橋墩周圍局部地形，當(dāng)?shù)匦螠y量值基本不變時，認(rèn)為沖刷達(dá)到平衡，然后測量橋墩沖刷地形。

恒定流橋墩沖刷試驗從沖刷坑的形成到基本沖刷平衡有一個大致固定沖刷時間，而潮流橋墩沖刷試驗時，需進(jìn)行特征水文年長系列潮周期的循環(huán)試驗才能使橋墩周圍的沖刷達(dá)到基本平衡。當(dāng)潮流最大流速與恒定流流速一致時，恒定流比潮流橋墩沖刷平衡所需時間少。通過300年一遇水文年潮型循環(huán)試驗，模型試驗時間14小時，不間斷監(jiān)測橋墩周圍沖刷地形，圖4為模型幾何比尺100沖刷坑發(fā)展過程線，從圖中可以看出，潮流沖刷歷經(jīng)8小時(模型時間)后橋墩最大沖深基本不變，為了保障各試驗達(dá)到?jīng)_刷平衡和節(jié)約試驗時間，初步確定橋墩沖刷試驗時間進(jìn)行10小時。在系列模型試驗過程中監(jiān)測發(fā)現(xiàn)10小時試驗橋墩沖刷均達(dá)到了平衡。

圖4 300年一遇潮型動力橋墩沖刷過程

3 試驗結(jié)果與分析

3.1 水流作用下橋墩基礎(chǔ)沖刷

在300年一遇和20年一遇設(shè)計水流作用下，水流受到橋墩基礎(chǔ)阻礙后，墩前壅水，橋墩兩側(cè)形成繞流向構(gòu)成表面漩輥(圖5)；中部以下水流遇樁群邊壁轉(zhuǎn)而向下，與下層水平方向行進(jìn)水流構(gòu)成底部的向下漩輥，底部的向下漩輥是產(chǎn)生局部沖刷的主要動力。局部沖刷達(dá)到相對平衡后，索塔基礎(chǔ)周圍的地形呈現(xiàn)機翼狀的沖刷形態(tài)，兩側(cè)沖刷地形基本對稱，沖刷區(qū)主要位于橋基承臺迎水面底部樁群間，最大沖深點處于迎水面的前三排樁群之間。從索塔基礎(chǔ)中部、緊靠承臺邊壁處的水流受到多重樁群的阻水消能，沖刷能力和攜沙能力同時遞減，形成類似蟹螯狀的淤積帶；橋基的背水面則呈現(xiàn)與橋基迎水寬度基本等寬、淤厚逐漸遞減的長條形淤積帶(圖6)。300年一遇和20年一遇設(shè)計水流條件下的橋墩局部沖刷形態(tài)基本相似，差別主要表現(xiàn)在沖深和范圍的大小。

圖5 索塔基礎(chǔ)水流流態(tài)(300年一遇設(shè)計水流)

圖6 索塔基礎(chǔ)周邊沖淤形態(tài)(300年一遇設(shè)計水流)

3.2 潮流作用下橋墩基礎(chǔ)沖刷

300年一遇水文年潮型潮流試驗中，潮流受漲落潮的影響呈往復(fù)流，橋基承臺及底部樁群的兩端既是迎水面，又是背水面。在漲潮和落潮期橋墩基礎(chǔ)兩端分別會形成單向流的水流特征，潮流過程有憩流、轉(zhuǎn)流期，在若干時段內(nèi)潮流流速小于300年一遇的設(shè)計水流流速，此時潮流受到橋墩基礎(chǔ)阻礙產(chǎn)生的繞流和漩輥范圍較小。青州航道橋橋墩縱軸線與潮流主流向約14°夾角，潮流作用下橋墩基礎(chǔ)沖刷分別進(jìn)行了無夾角和14°夾角試驗。橋墩縱軸線與潮流向無夾角時，橋墩周邊繞流與單向水流情況相似；橋墩縱軸線與潮流向有夾角時，橋墩周邊繞流不對稱與無夾角時略有區(qū)別(圖7)，試驗潮型的漲潮流速大于落潮流速，落潮流繞流范圍小于漲潮流，繞流形態(tài)與漲潮流基本相似。

圖7 漲急流速流態(tài)

潮流作用下橋墩基礎(chǔ)達(dá)到?jīng)_刷平衡所需要的時間明顯長于單向流情況。與水流作用下橋墩局部沖刷形態(tài)相同的是橋基附近形成局部沖坑，最大沖深點發(fā)生在流速大的迎水面樁群內(nèi)，不同特點是局部沖坑由一個發(fā)展為兩端各有一個，沖刷形態(tài)由勺形變?yōu)闃蚨諆啥松钪虚g淺的馬鞍形，橋墩兩端外緣形成與承臺基本等寬、堆高隨距離逐漸遞減的淤積帶(圖8),橋墩受潮流弱側(cè)淤積體大。

圖8 潮流作用下索塔基礎(chǔ)周邊沖淤形態(tài)

3.3 波流作用下橋墩基礎(chǔ)沖刷

在100年一遇和20年一遇設(shè)計水流與波浪共同作用下，橋墩水流結(jié)構(gòu)與單向流基本相同，由于波浪受橋墩的影響，在墩柱前、墩側(cè)面和墩后面分別形成波浪反射區(qū)、波浪散射區(qū)和波浪掩護(hù)區(qū)(圖9)。

圖9 波流共同作用下水流流態(tài)(100年一遇)

波流共同作用下的橋墩沖深范圍呈機翼狀向墩基兩側(cè)伸展，其伸展的角度與橋墩的迎水寬度、流速大小和與水流的交角有關(guān)。橋墩縱軸線與水流有交角時，墩基的沖刷形態(tài)兩側(cè)不對稱，沖刷范圍和最大沖深點偏向墩基與水流產(chǎn)生夾角的方位，同時產(chǎn)生了次沖深區(qū)(圖10)，沖出的泥沙落淤在橋墩掩護(hù)區(qū)域。在沖刷過程中，波浪主要起掀沙的作用，水流不僅對床面進(jìn)行沖刷，并將波浪掀起的泥沙向后輸移。當(dāng)水流導(dǎo)致沖坑深度較大時，波浪對沖刷坑的作用明顯減弱，在試驗條件下橋區(qū)的波流屬強流弱波，水流在沖刷中起主導(dǎo)作用。100年一遇設(shè)計水流疊加波浪條件下的橋墩周圍沖刷深度和范圍大于20年一遇設(shè)計水流疊加波浪條件下的情況。

圖10 波流共同作用下沖刷形態(tài)(100年一遇)

3.4 橋墩沖刷深度和范圍

由于水流受橋墩阻水的影響，在橋墩周圍水流結(jié)構(gòu)發(fā)生了變化，引起橋墩基礎(chǔ)周圍床面沖刷，形成局部沖刷坑，隨著沖刷加劇，水流動力對床底的作用逐漸減弱，而且河床較易沖刷的細(xì)顆粒泥沙沖走后，沖刷區(qū)泥沙也相對變?yōu)檩^粗顆粒的泥沙，抗沖能力增加，最后橋墩基礎(chǔ)周圍來流挾帶的落淤泥沙和沖刷泥沙達(dá)到平衡，形成橋墩基礎(chǔ)最大沖刷深度。采用系列模型延伸法可求出原型橋墩基礎(chǔ)最大沖刷深度，圖11和圖12分別為水流和波流作用下的橋墩沖刷系列延伸圖。

圖11 水流條件橋墩沖刷系列延伸圖

圖12 波流條件下橋墩沖刷系列延伸圖(夾角14°)

采用系列模型延伸法得到青州航道橋橋墩局部最大沖刷深度和范圍見表2，從表中數(shù)據(jù)可以看出，當(dāng)橋墩縱軸線與水流無夾角時，隨著水流流速的增大局部沖刷最大深度和范圍也增大；潮流作用下，最大流速與單向水流流速相同時，橋墩局部最大沖深和兩側(cè)沖刷范圍基本相近，橋墩兩端均出現(xiàn)淤積帶；波流共同作用下，由于受到波浪掀沙的影響，橋墩局部最大沖深和沖刷范圍大于單向水流情況。當(dāng)橋墩縱軸線與水流有夾角時，由于橋墩阻水面積的增大和夾角的影響，潮流、波流共同作用下，橋墩局部最大沖深和沖刷范圍大于無夾角情況，而且流速大和波浪強橋墩局部最大沖深和沖刷范圍均大。

表2 青州航道橋橋墩沖刷深度和范圍

4 結(jié) 語

橋墩局部沖刷特性和沖刷深度一直是研究的重點問題，是橋梁基礎(chǔ)設(shè)計和防護(hù)的重要參數(shù)。針對青州航道橋橋墩基礎(chǔ)周圍水動力條件復(fù)雜和橋墩基礎(chǔ)較大特點，采用系列物理模型試驗方法，分別進(jìn)行了水流、潮流和波流共同作用下的橋墩局部沖刷研究。試驗結(jié)果表明：1)潮流作用下橋墩沖刷過程由于有效流速作用時間較少，沖刷達(dá)到平衡需要的時間比恒定水流時間要長，橋墩基礎(chǔ)兩端均會發(fā)生迎流，沖刷形態(tài)有別于單向流沖刷形態(tài)，漲落潮形成兩個沖刷坑中的最大沖深和沖刷形態(tài)由流速較大的一側(cè)控制，橋墩局部沖刷達(dá)到?jīng)_刷平衡后，其最大沖深與相同水流流速的情況基本一致；2)波流共同作用下橋墩局部最大沖刷深度比單純水流沖刷深度增加10%左右，兩者的沖刷形態(tài)相近；3)相同動力條件下橋墩與水流有夾角時，橋墩局部沖刷深度和沖刷范圍會有所增加。

所建橋梁海床表層為淤泥質(zhì)亞黏土，其黏性強抗沖性能較好，目前物理模型試驗較難模擬，有待更深入的研究。從安全角度考慮，采用散粒體進(jìn)行試驗研究能夠滿足工程設(shè)計需求。港珠澳大橋2018年10月開通運營，青州航道橋索塔基礎(chǔ)由原來長八邊形調(diào)整為啞鈴型，直徑為36.5 m，根據(jù)2014年2月至2019年9月實測水下地形資料，橋墩周圍沖刷深度約1.5 m。