藍(lán)霄峰，張金明，鄭國(guó)賢

（珠江水利科學(xué)研究院，廣州 510611）

0 前言

過江隧道工程具有占地少、不影響河道正常通航等特點(diǎn)，近年來廣受重視，在進(jìn)行越江隧道的設(shè)計(jì)和修建中，隧道的埋深往往是影響隧道建設(shè)的關(guān)鍵因素，不同埋深對(duì)合理選取施工組織方式、掘進(jìn)速度等極為重要[1]。埋深小，可一定程度上減少投資，但一旦河床沖刷較大，出現(xiàn)隧道上方覆土層厚度小于設(shè)計(jì)標(biāo)準(zhǔn)或使隧道頂部出露，則工程運(yùn)營(yíng)安全將受到嚴(yán)重威脅，同時(shí)河床、堤岸穩(wěn)定也受到威脅；埋深過大，工程安全可得以保障，工程投資及施工難度會(huì)增加[2]。

目前國(guó)內(nèi)外河道極限沖刷這類研究較少，近年來，史英標(biāo)等[3-4]人對(duì)錢塘江河口過江隧道河段極端洪水條件下的沖刷深度開展了深入研究，取得了較好的成果，但研究主要針對(duì)河道受洪潮共同作用，底、懸泥沙粒徑差異不大，河底表層覆蓋層為粉質(zhì)黏土層情況；張為等[5]對(duì)長(zhǎng)江下游南京河段過江隧道最大沖刷深度采用數(shù)學(xué)模型進(jìn)行了研究，主要考慮三峽水庫(kù)的影響，基于河床覆蓋層為床沙質(zhì)和沖瀉質(zhì)。楊芳等[6]對(duì)珠江三角洲東平水道河段最大沖刷進(jìn)行相關(guān)研究，基于感潮河道特性及全風(fēng)化泥層。吳門伍等[7]對(duì)內(nèi)河邕江最大沖刷進(jìn)行相關(guān)研究，基于典型內(nèi)河特性及河床覆蓋層為圓礫層。南渡江是我國(guó)唯一的一條大型熱帶河流，其水文等自然條件與我國(guó)其他大型河流存在較大差異，河口更是經(jīng)常受到風(fēng)暴潮的影響，河道坡降平緩，沙洲、小丘及淺灘較多，河床覆蓋層以粗砂為主，目前尚未見報(bào)道相關(guān)研究工作，因此有必要對(duì)南渡江河口處穿越南渡江隧道河段進(jìn)行極限沖刷深度預(yù)測(cè)，為過江隧道的合理埋設(shè)提供科學(xué)依據(jù)。

1 工程概況

1.1 河道概況

?？谑续P翔東路越江通道項(xiàng)目橫跨南渡江，連接江西主城和江東新區(qū)，主城區(qū)接線為現(xiàn)狀主干路鳳翔東路，江東接線為規(guī)劃四號(hào)路，距離上游海瑞大橋1400 m，距離下游瓊州大橋2200 m，為東西向城市主干道。

南渡江是海南省第一大河流。發(fā)源于白沙縣南峰山，流經(jīng)白沙縣、儋州市、瓊中市、屯昌縣、澄邁縣、定安縣和海口市，于海口市三聯(lián)村向北流入瓊州海峽。干流全長(zhǎng)333.8 km，落差703 m，平均坡降0.72‰。多年平均流量219 m3/s[8]。

松濤水庫(kù)壩址以上為南渡江上游，河長(zhǎng)137 km，為中低山地區(qū)，河谷狹窄，坡降大，急灘多，兩岸地形陡峻，高程都在500 m 以上，最高點(diǎn)鸚歌嶺1812 m。松濤水庫(kù)壩址以下至九龍灘為南渡江中游，河長(zhǎng)83 km，屬低山丘陵，南高北低，一般山頂高程200～500 m，最高點(diǎn)黎母嶺1411 m，山間谷溝發(fā)育，河道迂回彎曲，兩岸坡陡。九龍灘以下為南渡江下游，河長(zhǎng)114 km，屬丘陵臺(tái)地及濱海平原三角洲，地勢(shì)南高北低，河道寬闊，坡降平緩，沙洲、小丘及淺灘較多，兩岸是平坦的臺(tái)地，大部分為農(nóng)田，其中潭口以下為河口段，最下游的梯級(jí)龍?zhí)链髩坞x河口26 km。

1.2 水文特征

南渡江為山區(qū)性河流，流域內(nèi)徑流由降雨形成，由于流域上中下游地理位置和氣候特性不同，造成徑流在時(shí)間和空間分布上差異較大。汛期（5～10月）徑流量集中，占全年總徑流量的70%以上，夏秋臺(tái)風(fēng)暴雨時(shí)，常出現(xiàn)潮水倒灌現(xiàn)象，時(shí)有洪潮災(zāi)害。據(jù)龍?zhí)琳舅馁Y料統(tǒng)計(jì)：實(shí)測(cè)最大年平均流量296 m3/s，徑流量93.3億m3（1973年）；實(shí)測(cè)最小年平均流量74.5 m3/s，徑流量23.5 億m3（2004 年）；實(shí)測(cè)多年平均流量182 m3/s，徑流量57.4億m3（1955—2014年系列）；實(shí)測(cè)最小流量1.40 m3/s（1976年）；實(shí)測(cè)最大流量9300 m3/s（2000年）。

1.3 泥沙特征

南渡江具有豐水少沙特點(diǎn)，據(jù)龍?zhí)了恼窘y(tǒng)計(jì)，多年平均徑流量為66.81×10 m3，平均輸沙量為53.12×104t。流域水沙主要來源于南渡江干流和龍州河[9]。南渡江徑流和輸沙主要集中在汛期，季節(jié)分配不均，其中汛期徑流量占到全年的80%左右，輸沙量占全年的90%左右，汛期河流攜帶大量泥沙入海[10]。根據(jù)龍?zhí)了恼緫乙瀑|(zhì)泥沙觀測(cè)資料統(tǒng)計(jì)（見表1），1959年松濤水庫(kù)建成后，平均含沙量呈減小趨勢(shì)，懸移輸沙量明顯減少。南渡江河口不具備“夏儲(chǔ)冬輸”的泥沙特征，而是具有臺(tái)風(fēng)季節(jié)特征，河口水體含沙量近岸大于離岸，底層大于表層，洪季含沙量明顯高于枯季[11]。南渡江泥沙變化特點(diǎn)使得隧址河段極限沖刷深度研究變得非常必要。

表1 龍?zhí)了恼緫乙瀑|(zhì)含沙量和輸沙量變化

1.4 河床邊界條件

南渡江下游流域地貌特征屬河流沖積階地與玄武巖剝蝕臺(tái)地地貌，地形總體趨勢(shì)南高北低。南渡江下游段向北匯入瓊州海峽，在入海口處形成現(xiàn)代河口。河道大致自南渡江大橋后加入開闊的一級(jí)階地三角洲平原區(qū)，其上游為玄武巖噴出巖剝蝕二級(jí)階地，高程5～10 m，二級(jí)階地表面普遍履蓋了幾米厚的紅色全風(fēng)化黏土，地形略有起伏。河道河漫灘及江心洲發(fā)育完整。河床及河漫灘表面為現(xiàn)代河砂沉積，東岸土堤大部分座落于河漫灘松軟沉積層上。

工程地勘資料顯示，河床表層粗砂，厚度為4.2～9.5 m，平均厚度為5.8 m，泥沙粒徑為1.0～0.5 mm；第二層為粉質(zhì)黏土，厚度為15.0～17.0 m，平均厚度為16.4 m，泥沙粒徑為0.03～0.05 mm；第三層為礫砂，厚度為4.5～5.6 m，平均厚度為5.0 m，泥沙粒徑為1～2 mm；第四層為粉質(zhì)黏土，泥沙粒徑為0.04～0.05 mm。隧道穿越的土層為第二、三、四層。從河道河床覆蓋層來看，其泥沙粒徑及構(gòu)成相對(duì)于其他大型河流河口河床的有明顯不同。

2 河床演變

工程區(qū)處于南渡江河口段，屬低階地和平原區(qū)，河床充填砂礫質(zhì)沉積物，江心洲和淺灘發(fā)育。根據(jù)歷史遙感圖影響顯示，河口段左、右岸堤防未建之前，南渡江出龍?zhí)翇{谷后，洪水一經(jīng)出槽便似放射狀漫流入海，自然狀況下南渡江河口段具游蕩性河流的特點(diǎn)，主槽常有偏移，歷史上河流改道時(shí)有發(fā)生。

（1）河床沖淤變化。隧址河段沖淤變化如圖1所示。2008—2011年，整體上隧址河段呈現(xiàn)槽沖灘淤的態(tài)勢(shì)，工程上下游有幾處沖刷深度達(dá)數(shù)米，其中工程下游左岸局部沖刷達(dá)4 m，工程下游司馬坡島左汊總體呈淤積態(tài)勢(shì)，右汊呈沖刷趨勢(shì)；2011—2018 年，工程下游司馬坡島左汊繼續(xù)呈淤積態(tài)勢(shì)，右汊輕微沖刷，“左支右干”態(tài)勢(shì)更加明顯；2018—2021 年，隧址河段總體呈現(xiàn)緩慢回淤的態(tài)勢(shì)，沖刷坑基本消失，工程下游司馬坡島左汊邊灘淤積明顯，右汊也有輕微淤積?？傮w來看，隧址河段呈“槽沖灘淤”的態(tài)勢(shì)。

圖1 隧址河段2008—2021年沖淤變化圖

（2）深泓線高程變化（見圖2）。由圖2可見，隧址河段近13年來深泓線有抬升亦有下切，整體以下切為主。2008—2011 年，深泓線整體明顯下切，最大下切幅度約1.3 m，平均下切幅度約0.9 m；2011—2018 年，深泓線有抬升也有下切，最大下切幅度約為2.6 m，最大抬升幅度約為0.5 m，整體下切約為0.6 m；2018—2021年，深泓線有抬升也有下切，最大下切幅度約為0.1 m，最大抬升幅度約為0.6 m，整體抬升約為0.3 m。

圖2 隧址河段沿程深泓線高程變化圖

（3）隧址斷面高程變化（見圖3）。由圖3可見，2008—2011 年，左支主槽明顯出現(xiàn)沖刷，最大沖刷深度約為4.1 m，右汊主槽整體也呈略沖刷形態(tài)，平均沖刷深度約0.4 m；2011—2018年，左支主槽略微淤積，最大淤積幅度約為0.3 m，右汊主槽出現(xiàn)很大沖刷，最大沖刷深度約為3.4 m；2018—2021 年，斷面略微回淤。

圖3 隧址斷面多年地形變化圖

由于隧址河段近年來呈“槽沖灘淤”的態(tài)勢(shì)，深泓下切幅度很大，隧址斷面也出現(xiàn)沖刷深度很大現(xiàn)象，因此進(jìn)行隧址河段的極限沖刷試驗(yàn)研究是非常必要的。

3 物理模型設(shè)計(jì)及驗(yàn)證

3.1 模型設(shè)計(jì)

過江隧道工程位于河流分汊段，為確保水流相似，物理模型上、下游邊界須選取河道順直段。根據(jù)研究?jī)?nèi)容和相關(guān)規(guī)范，考慮到工程河段的水流特性、河道地形等因素，選定物理模型（見圖4）模擬范圍為距離工程位置上、下游各約3 km 原型；動(dòng)床模擬范圍為隧址軸線上下游各0.8 km 的河道。根據(jù)研究問題的特點(diǎn)確定模型為變態(tài)模型，模型平面比尺λl為180，垂直比尺λh為100，相應(yīng)的變率η為1.8，這樣可以減輕因變率偏大引起的泥沙運(yùn)動(dòng)相似性偏離之程度，同時(shí)經(jīng)計(jì)算亦可滿足避免表面張力影響的要求。根據(jù)相似原理可得相應(yīng)的模型流量比尺λ0為180 000，流速比尺λv為18，糙率比尺λn為1.61。

圖4 物理模型示意圖

根據(jù)試驗(yàn)任務(wù)及要求，試驗(yàn)按最不利因素考慮河床的極限沖刷，因此，模型試驗(yàn)不考慮懸移質(zhì)泥沙運(yùn)動(dòng)，而只考慮泥沙推移質(zhì)運(yùn)動(dòng)。通過對(duì)多種不同模型沙起動(dòng)流速進(jìn)行分析比較，確定選用中值粒徑為0.3 mm，容重為1.15 t/m3的塑料沙作為模型沙。通過水槽起動(dòng)試驗(yàn)測(cè)得，該模型沙在水深5～10 cm時(shí)（相當(dāng)原型5～10 m）的起動(dòng)流速為5.0～7.0 cm/s，與試驗(yàn)所需起動(dòng)流速基本一致，基本滿足研究的起動(dòng)相似要求。按文獻(xiàn)[12]設(shè)計(jì)模型的有關(guān)相似條件。

3.2 清水定床模型驗(yàn)證及動(dòng)床相似性分析

物理模型對(duì)清水定床和渾水動(dòng)床兩方面進(jìn)行驗(yàn)證。清水定床試驗(yàn)采用數(shù)值模擬計(jì)算結(jié)果[13]，100年一遇洪峰流量對(duì)模型進(jìn)行水動(dòng)力相似驗(yàn)證，試驗(yàn)布共設(shè)了11個(gè)流速流向測(cè)點(diǎn)和5個(gè)水位測(cè)點(diǎn)，測(cè)點(diǎn)位置見圖5，驗(yàn)證結(jié)果見表2 和表3，驗(yàn)證結(jié)果符合規(guī)范要求[14]。動(dòng)床模型試驗(yàn)主要采用規(guī)范公式計(jì)算來確定各類比尺，試驗(yàn)主要從以下兩方面考慮：①在模型控制上，動(dòng)床模型與定床驗(yàn)證試驗(yàn)水流控制條件相同；②合理選擇模型沙，本次試驗(yàn)進(jìn)行了模型沙的起動(dòng)流速水槽試驗(yàn)，模型選用模型沙的特性與原型沙在起動(dòng)相似較為一致。因此，河道沖淤變形的精度滿足河工模型試驗(yàn)規(guī)程的要求，為過江隧道河段的動(dòng)床沖刷試驗(yàn)奠定了可靠基礎(chǔ)。

圖5 驗(yàn)證點(diǎn)位置示意圖

表2 模型流速驗(yàn)證成果

表3 模型水位驗(yàn)證成果m

4 隧址河段極限沖刷深度預(yù)測(cè)

4.1 邊界條件

根據(jù)項(xiàng)目水文測(cè)驗(yàn)及分析計(jì)算成果[13]可知，河口感潮河段水位受洪潮共同影響，但河口分汊段以下河段百年一遇洪水位普遍高于百年一遇風(fēng)暴潮位，即使洪水與天文潮遭遇，河口段仍由洪水控制。風(fēng)暴潮發(fā)生在前，洪峰出現(xiàn)在后，河口段洪峰總滯后于風(fēng)暴潮1～2 d，當(dāng)洪峰到達(dá)?？跁r(shí)，此時(shí)風(fēng)暴潮亦消失，因此，龍?zhí)磷罡吆榉逅慌c風(fēng)暴潮水位遭遇的機(jī)率極小。南渡江區(qū)域洪水主要由暴雨形成，洪水發(fā)生時(shí)間與暴雨一致，而天文大潮持續(xù)時(shí)間長(zhǎng)，洪水遭遇外海低潮位的幾率較大。結(jié)合隧道工程防洪標(biāo)準(zhǔn)按100年一遇設(shè)計(jì)，300年一遇校核，本次物理模型試驗(yàn)上游水文邊界選取50 年（10 500 m3/s）、100年（11 800 m3/s）和300年一遇典型頻率洪水（13 500 m3/s）遭遇河口多年平均低水位進(jìn)行河床沖刷試驗(yàn)分析。

4.2 極限沖刷深度預(yù)測(cè)

隧址河段水流流態(tài)見圖6～8，隧址河段河床沖刷形態(tài)見圖9～11。

圖6 隧址河段水流流態(tài)（P=2%）

圖7 隧址河段水流流態(tài)（P=1%）

圖8 隧址河段水流流態(tài)（P=0.33%）

圖9 隧址河段河床沖刷形態(tài)（P=2%）

圖10 隧址河段河床沖刷形態(tài)（P=1%）

圖11 隧址河段河床沖刷形態(tài)（P=0.33%）

（1）隧址上游斷面沖刷。隧址上游河床沖刷主要發(fā)生于主槽內(nèi)，兩側(cè)高邊灘高程基本沒有變化，主槽平均沖刷深度在2.52～3.01 m之間，最大沖刷深度在3.93～4.43 m 之間，最深點(diǎn)高程在-6.6～-7.2 m之間，斷面最深點(diǎn)向下游有增大趨勢(shì)。

（2）隧址斷面沖刷。兩側(cè)高邊灘地形變化較小，主槽河床整體下切，平均沖刷深度2.58 m；河床最大沖刷深度位置出現(xiàn)在主槽左側(cè)，最大沖刷深度約3.91 m，沖刷深度后河床高程為-6.80 m。斷面最深點(diǎn)河床沖刷深度約3.06 m，沖刷后最深點(diǎn)高程至-7.90 m。

（3）隧址下游斷面沖刷。隧址下游750 m范圍內(nèi)，河床沖刷主要發(fā)生在司馬坡島的兩側(cè)，這主要受司馬坡島分流、堤防硬邊界以及紊亂的水流綜合影響，主槽平均沖刷深度在1.98～3.38 m 之間，最大沖刷深度在2.97～5.25 m之間，最深點(diǎn)高程在-8.4 m左右，斷面最深點(diǎn)向下游有增大趨勢(shì)。

5 結(jié)語(yǔ)

洪潮及臺(tái)風(fēng)作用下的南渡江河口河床沖淤幅度較大，在此該河段埋設(shè)隧道，極端洪水情況下極限沖刷深度是設(shè)計(jì)關(guān)鍵的參數(shù)之一。本文通過河演分析和物理模型試驗(yàn)多學(xué)科研究手段對(duì)過江隧道河段極限沖刷深度進(jìn)行預(yù)測(cè)。結(jié)果表明：在不同水文年條件下工程引起的河床沖淤變化，隧址段河床沖刷主要發(fā)生在現(xiàn)狀主槽內(nèi)及司馬坡島的兩岸附近，河道等高線總體趨勢(shì)表現(xiàn)為向兩側(cè)推進(jìn)，主槽寬度加大；受司馬坡島分流、堤防硬邊界以及紊亂水流的綜合影響，斷面最大沖刷深度向下游有增大趨勢(shì)；隧址河段主槽最深點(diǎn)沖刷深度3.06 m，對(duì)應(yīng)的最小底高程為-7.90 m。預(yù)測(cè)成果可作為過江隧道的合理埋設(shè)提供科學(xué)依據(jù)。