黃廷杰， 吳 忠， 邵 勇， 陸 彥， 吳 攀

(1.南京水利科學(xué)研究院， 江蘇 南京 210029； 2.江蘇省太湖治理工程建設(shè)管理局， 江蘇 常州 213022)

新孟河屬太湖流域湖西區(qū)沿江水系，新孟河延伸拓浚工程是《太湖流域水環(huán)境綜合治理總體方案》安排的提高水環(huán)境容量(納污能力)引排工程關(guān)鍵實(shí)施項(xiàng)目之一，河道全長116.69 km，平水年引長江水入湖水量達(dá)到25.2億m3，枯水年流域水資源配置引江入湖水量21.4億m3，100年一遇流域防洪排江水量7.9億m3，對改善太湖西北部水域和上游地區(qū)河湖水環(huán)境具有重要的意義。

1 區(qū)域概況

新孟河延伸拓浚工程設(shè)計(jì)中長江—北干河段拓寬河道底寬至50～80 m，底高程為-3.0 m(圖1)。該段工程中長江至大運(yùn)河段沿線兩岸支河口門實(shí)施控制，而新孟河穿京杭運(yùn)河后平地開河至北干河之間的鶴溪河、夏溪河、湟里河等主要支河沒有實(shí)施有效控制。支河現(xiàn)狀底寬8～12 m，底高程0～1.8 m。無閘控支河將影響新孟河引水工程運(yùn)行的穩(wěn)定與工程效益的發(fā)揮，同時引水也將直接影響支河的水動力泥沙過程。為保障新孟河引水后工程有效運(yùn)行，同時減少引水對支河生境的負(fù)面影響，需要探明新孟河引水沿線交無閘控支河交匯河口的水動力泥沙變化過程及其特征，制定有效的修復(fù)對策。

圖1 新孟河引水沿線無閘控支河概化圖

2 二維水動力泥沙數(shù)學(xué)模型建立與驗(yàn)證

2.1 控制方程

(1)二維水動力模型

連續(xù)性方程為

(1)

x方向動量方程為

(2)

y方向動量方程為

(3)

式中：h為水深；u為x方向的流速；v為y方向的流速；t為時間；g為重力加速度。

初始條件為二維模型給定各計(jì)算網(wǎng)格點(diǎn)上水位、流速初值，即

(4)

(2)二維泥沙輸運(yùn)模型

在水流模型的基礎(chǔ)上進(jìn)行泥沙運(yùn)動數(shù)值模擬，懸移質(zhì)不平衡輸沙控制方程為

(5)

式中：s為泥沙的含沙量；E為泥沙的沖刷通量；D為泥沙的沉降通量。

泥沙輸運(yùn)模型采用隱式求解，含沙量定義在單元的中心。在單元i上積分方程(5)有

(6)

2.2 模型建立

在二維河網(wǎng)計(jì)算中，河道斷面比較順直，采用矩形網(wǎng)格計(jì)算效率較高。而在河口區(qū)域計(jì)算中，由于滆湖北側(cè)區(qū)域島嶼較多，岸線曲折復(fù)雜，采用非結(jié)構(gòu)網(wǎng)格對計(jì)算域進(jìn)行剖分是非常合適的。采用非結(jié)構(gòu)網(wǎng)格剖分計(jì)算域，既可以克服矩形網(wǎng)格鋸齒形邊界所造成的流動失真，也可以避免生成有結(jié)構(gòu)貼體曲線網(wǎng)格的復(fù)雜計(jì)算和其他困難。因此，為了更好地擬合研究區(qū)域的形狀以及島嶼等，采用的網(wǎng)格為三角形與四邊形混合的非結(jié)構(gòu)網(wǎng)格單元組成。

為了更好地反映各支河與新孟河交匯情況、支河入滆湖水沙情況，模型計(jì)算區(qū)域進(jìn)口邊界為界牌樞紐支河口門和扁擔(dān)河在京杭運(yùn)河口門，出口邊界為鶴溪河交新孟河西側(cè)匯流、夏溪河交新孟河西側(cè)出流、湟里河交新孟河西側(cè)出流、北干河交新孟河西側(cè)出流、扁擔(dān)河、孟津河的匯流，滆湖以及武南河、太滆運(yùn)河、漕橋河、殷村港、燒香港等出湖河道的出流。此處采用的網(wǎng)格為三角形與四邊形混合的非結(jié)構(gòu)網(wǎng)格單元組成。模型網(wǎng)格如圖2所示。

圖2 模型網(wǎng)格

2.3 模型驗(yàn)證

根據(jù)《新孟河初步設(shè)計(jì)總報告》的設(shè)計(jì)資料對模型進(jìn)行驗(yàn)證(表1～2)，在設(shè)計(jì)流量下各斷面水位流量為瞬時最大值，只驗(yàn)證日均引水流量的水位。結(jié)果表明，在日均引水流量(Q=320 m3/s)條件下，運(yùn)河以北段(長江—運(yùn)河)的入江口斷面(界牌樞紐)、過京杭運(yùn)河斷面(立交地涵)、南延段(運(yùn)河—北干河)的3條支河以及北干河斷面計(jì)算水位與設(shè)計(jì)值相差0.10 m以內(nèi)。本研究建立的數(shù)學(xué)模型計(jì)算水位、流量結(jié)果與初設(shè)報告結(jié)果比較吻合。

表1 水位模型驗(yàn)證結(jié)果

表2 流量模型驗(yàn)證結(jié)果

3 引水對無閘控支河水動力的影響

水動力計(jì)算結(jié)果表明，新孟河引水后與鶴溪河、夏溪河和湟里河的交匯河口水位平均抬高10～20 cm，但交匯區(qū)域流速明顯增加(表3)。

日均引水流量(Q=320 m3/s)條件下，新孟河沿線與各支河交匯河口西側(cè)平均流速為0.14～0.26 m/s，東側(cè)平均流速為0.14～0.19 m/s。其中，新孟河沿線與鶴溪河、夏溪河交匯斷面的西側(cè)平均流速高于東側(cè)，而與湟里河交匯斷面東側(cè)平均流速較西側(cè)略大(增幅約0.05 m/s)。以新孟河為界，從3條支河與新孟河的交匯地貌和流場分布來看(圖3)。鶴溪河上游流場分布均勻，且水流仍然向下游流動，表明新孟河與鶴溪河上游入?yún)R水流較為平順，且受到新孟河引水影響，上游水流受到頂托作用，造成流速減小。而鶴溪河下游水流受到新孟河截流及引水影響，口門區(qū)域流速分布較為不均，平均流速也較上游減半，這也可能受到下游急彎段水流頂托的影響；夏溪河河型較為平直，與新孟河呈垂直交匯，因此交匯口門兩側(cè)流場均有不同程度的擾動，而夏溪河上游水流溯回，表明新孟河分流水流功率大于原河道，下游流場分布較為均勻，上游流速普遍高于下游(小于0.1 m/s)，這可能是受到交匯口門高程影響；湟里河與新孟河交匯區(qū)域恰處于2個連續(xù)彎道段，兩側(cè)入?yún)R水流受一定程度彎道阻力的影響，從而使流速減小，口門區(qū)域流場分布也不均勻。但上游水流依然溯回，且上下游流速大小相當(dāng)，下游受到彎道形態(tài)影響，流場分布較為不均。交匯區(qū)域東側(cè)平均流速大于西側(cè)，則是由于下游順?biāo)鞣较蛉雲(yún)R，且水流運(yùn)動軌跡與彎道凹岸深槽相一致。

表3 不同引水條件下新孟河沿線與支河交匯河口斷面平均流速 單位：m/s

圖3 不同引水流量條件下新孟河引水沿線與無閘控支河交匯河口的流速分布

設(shè)計(jì)引水流量(Q=745 m3/s)條件下，相比日均引水流量，新孟河沿線與各支河交匯河口及上下游河道流速普遍增加，西側(cè)平均流速為0.16～0.42 m/s，東側(cè)平均流速為0.19～0.25 m/s，各支河交匯河口的流速特征與日均引水流量條件下的流速分布變化又有不同，各交匯河口流速分布不均勻范圍繼續(xù)向兩側(cè)擴(kuò)大，這當(dāng)然是受到各入?yún)R河口處水流動力增強(qiáng)的緣故。鶴溪河上游段水流溯回，而各支河下游段水流流速分布則基本以河道深槽為界，下游水流順河道入?yún)R與原河道地貌形態(tài)變化相一致。湟里河上游河道流速變化不大，這也可能是由于水位抬升后河道斷面大幅展寬的影響。

另外，新孟河引水后在支河交匯河口的北岸區(qū)域存在低速區(qū)，流速0.1 m/s左右，而交匯處南岸區(qū)域受到來流頂沖，流速較大，流速最高可達(dá)到0.5 m/s。水流進(jìn)入支河后，沿程逐漸恢復(fù)為均勻水流。流量不同時，河口段橫斷面的流速分配不均勻程度也不同，日均流量時低速區(qū)域較大，而設(shè)計(jì)流量下支河口口門區(qū)域斷面流速差異性較小。

4 引水對無閘控支河泥沙輸移的影響

長江水質(zhì)較好，但含沙量相比各支河河道較大，從新孟河引水干流沿線河道沖淤變化來看(圖4)，河道泥沙總體呈現(xiàn)由快速沉降到緩慢沉降的趨勢，且隨著含沙量增加，沿線泥沙呈現(xiàn)加速沉降趨勢。

圖4 新孟河引水沿線與無閘控支河交匯河口的沖淤分布

新孟河干線從上游向下游主要表現(xiàn)為先淤后沖的趨勢，大部分泥沙淤積在新孟河口至鶴溪河河段。鶴溪河、夏溪河與新孟河引水沿線交匯河口兩側(cè)地形受到入?yún)R水流動力較強(qiáng)，且河勢較為平順，因此呈沖刷趨勢。而湟里河與新孟河沿線交匯河口受到兩側(cè)彎道形態(tài)影響，口門呈現(xiàn)微淤趨勢，年沖淤幅度約3～9 cm。鶴溪河整體呈現(xiàn)微淤趨勢(約0～3 cm)。各支河上游段以夏溪河淤積程度最高(約6～12 cm)，考慮其水動力條件受到引水后分流影響更強(qiáng)，口門的拓寬也使得上游分沙增大，造成上游普遍淤積，而下游交匯河口高程阻止了更多的分沙進(jìn)入，河道整體呈現(xiàn)微沖趨勢。隨著泥沙向下游運(yùn)動并不斷淤積，湟里河與新孟河交匯區(qū)域水體含沙量已經(jīng)大幅下降，彎道地貌形態(tài)阻止了更多的含沙量進(jìn)入上下游河道，但水流慣性作用繼續(xù)增強(qiáng)，造成上下游河段普遍沖刷(約0～9 cm)。下游新孟河干線沖刷強(qiáng)度則繼續(xù)增大(最大約12 cm)。

由此看來，鶴溪河與夏溪河支河河口底質(zhì)將從淤泥質(zhì)基質(zhì)向粉砂質(zhì)基底發(fā)展，而湟里河口門的不斷淤積將有可能影響上下游引水，從而影響該河段沿河生態(tài)條件。夏溪河下游及湟里河河道沖刷也會增加兩岸岸坡穩(wěn)固及生態(tài)風(fēng)險。

5 引水工程交匯河口的修復(fù)對策

5.1 岸坡防護(hù)

新孟河引水對無閘控支河交匯河口的沖刷直接影響了河岸地質(zhì)結(jié)構(gòu)的穩(wěn)定性，在今后應(yīng)特別注意加強(qiáng)對各支河與新孟河沿線交匯河口兩岸及各支河下游右岸岸坡的防護(hù)，保障新孟河引水干流通道的暢通。同時，引水后水位上升導(dǎo)致河水對支河岸坡的滲透作用增強(qiáng)，覆蓋植被對岸坡的“錨固”作用被削弱，岸坡穩(wěn)定性下降。

交匯河口的岸坡可采用防護(hù)性強(qiáng)、生態(tài)友好型的護(hù)坡形式，如生態(tài)護(hù)坡磚、鋼絲網(wǎng)格護(hù)坡等。生態(tài)護(hù)坡磚具有較強(qiáng)的整體性，可采用連鎖式設(shè)計(jì)，塊體之間拼裝緊密，能抵抗較大流速的沖刷；鋼絲網(wǎng)格采用耐腐蝕、高強(qiáng)度、柔性好的鋼絲編織成網(wǎng)籠，其內(nèi)充填石料而成，具有整體性和透水性好、可適應(yīng)變形、防水流沖刷等優(yōu)點(diǎn)。兩種護(hù)坡形式均可配合覆土和種植耐淹植物[2-3]，營造一定的生態(tài)景觀。

5.2 生態(tài)修復(fù)

2018年與2019年分別對交匯河口及支河的生態(tài)環(huán)境進(jìn)行了全面的調(diào)查，支河沿線河濱帶分布有大片蘆葦、菖蒲、菹草等水生植物。引水后適宜生長環(huán)境上移將使得現(xiàn)有挺水植物部分消亡，新擴(kuò)展的河岸水深1.0 m區(qū)域內(nèi)將逐漸成為挺水植物的主要分布區(qū)[4]。同時，在支河交匯河口流速增加(可高達(dá)0.5 m/s)，對柔性沉水植物可能造成割裂[5]?？煽紤]在支河交匯河口附近營造深潭地形，防止口門泥沙淤積；在交匯口下游的自然岸帶，采用如植物扦插生態(tài)護(hù)坡、土工織物草皮護(hù)坡、杉木樁護(hù)岸等形式，重建河岸帶生境，創(chuàng)造出適宜于水生生物生長的棲息環(huán)境。

6 結(jié) 語

新孟河延伸拓浚工程建成后引水將對沿線無閘控連通支河主要帶來水動力泥沙、生態(tài)等方面的影響，同時反過來水動力泥沙和生態(tài)也影響引水工程的穩(wěn)定運(yùn)行，為此需要提前制定修復(fù)對策。水動力泥沙模型研究表明，新孟河引水后將使得無閘控支河水位升高、流速增大，對岸坡穩(wěn)定性和生態(tài)產(chǎn)生影響，需引起足夠重視。建議在無閘控支河與新孟河引水沿線交匯河口及其下游河段進(jìn)行岸坡防護(hù)和生態(tài)修復(fù)，保證提高新孟河引水工程效益。