周波

摘要： 運用MIKE11構建水量水質耦合模型，模擬現狀河網和增加航塘港閘后的河網。首先模擬引清調水中水質代表點的污染物濃度變化過程，研究航塘港水閘對周邊鎮(zhèn)區(qū)河道水質的改善作用，然后模擬莫拉克臺風以及50年一遇暴雨期的水位變化過程，分析航塘港水閘的防汛排澇作用。研究成果能為航塘港水閘建設提供參考依據。

關鍵詞： 水文模型 水閘 引清調水防汛排澇

1 前言

奉賢區(qū)位于上海市南部，地處黃浦江以南、杭州灣以北，兩面瀕水，東與浦東相鄰，西與松江、金山接壤，全區(qū)面積約720Km2。區(qū)域內市、區(qū)級骨干河道19條，鎮(zhèn)級河道260多條，河道面積為55.39 km2。

奉賢區(qū)水務局十三五規(guī)劃中，航塘港將延伸約3公里至杭州灣，并在該河出?？谠O置水閘。航塘港水閘的建設將對奉賢區(qū)的防汛以及調水產生作用，本文通過水文模型的模擬試對該作用進行分析，為航塘港水閘的建設提供參考依據。

2水利概況

奉賢區(qū)境內共有12座水閘，其中沿江7座，全部分布在金匯港以西。沿海4座，金匯港以西2座（金匯港南閘、南竹港南閘），以東2座（南門港閘、中港閘）。西部邊界1座。規(guī)劃在金匯港以東地區(qū)新建1座沿海閘門。水閘分布見圖1。

3水文模型簡介

MIKE11是丹麥水力研究所研制的一維河網模擬軟件，軟件包括水動力模塊、降雨徑流模塊、對流擴散模塊、可控水工建筑物模塊、水質實驗室模塊等，是進行河網水動力和環(huán)境狀態(tài)模擬的有力工具。

MIKE11軟件中的水動力模塊（HD）模擬河網水動力狀態(tài)。該模塊基于描述簡單明渠非恒定漸變流運動規(guī)律的圣維南方程組建立，包括反映動量守恒定律的運動方程和反映質量守恒定律的連續(xù)性方程：

（1）

（2）

式（1）～（2）中，A為河道內任意斷面的面積；Q為斷面流量；t、x分別為時間與位置的坐標；h為水位；g為重力加速度；q為旁側入流；R為水力半徑；C為謝才系數；為動量校正系數。模型采用6 點Abbott-Ionescu有限差分格式求解，在給定的邊界條件下推求近似結果。

降雨徑流模塊（RR）包含幾種模型，其中主要是NAM模型。該模型描述水文循環(huán)的地表過程，是概念性集總式模型，其將整個（子）流域視作一個整體單位，并考慮水文循環(huán)的各種物理過程，模型參數代表平均情況，率定目標是總水量平衡和徑流過程線形狀相似。NAM模型的結構圖見圖2。

對流擴散模塊（AD）可模擬溶解態(tài)或懸浮物質的對流擴散和一級降解過程。污染物質在水中的濃度與自身降解、水流運動以及擴散有關，其控制方程采用物質輸移的對流擴散方程：

（3）

式（3）中：C為某污染物的水質濃度；A為斷面面積；Q為斷面流量；Ex為縱向擴散系數；K為生物降解系數；C2為外部匯入同種污染物的濃度；q為帶有該污染物的旁側入流量。該方程為偏微分方程，通過與水動力方程聯合進行求解。

模塊（SO）為可控制建筑物模型，可設置參數為水閘類型、閘寬、閘底高程、閘門開啟高度、河網中所在位置、調度原則等。新建水閘可在該模塊中進行設置。

4 模型構建

利用上述模塊可構建水量水質耦合模型，根據實測資料對模型參數進行率定，使水位流量過程基本與實際吻合。假定某不可降解污染物質在河網中的初始狀態(tài)和邊界狀態(tài)，擴散系數取經驗值建立水質模型。通過水量水質模型耦合計算，能在一定程度上真實反映污染物質在河網中的傳輸過程。

4.1水量模型

模型對浦東水利控制片區(qū)的河網進行了概化，其中河道斷面采用2011年水利普查數據。

模型邊界：西部邊界河道交換水量很少，故將西部邊界處理為閉邊界。沿江沿海邊界以實測或插補潮位為邊界。初始狀態(tài)：水位為2.5m，流量為0m3/s。

模型計算時段為2013年6月1日0：00至12月31日0：00，模型中HD模塊計算步長為5min，RR模塊計算步長為1hr。

4.2水質模型

水質邊界初始狀態(tài)：假定不可降解污染物A在黃浦江和杭州灣中的濃度為0mg/L，在河網中的濃度為2mg/L。

模型參數：市、區(qū)級河道擴散系數為20，其他河道擴散系數為10；降解系數為0。

耦合計算時段：2013年9月16日0：00至9月23日0：00。模型計算步長為5min。

因缺乏足夠的連續(xù)的水質數據，故水質模型沒有進行率定。

5航塘港閘調水作用分析

5.1污染物對流擴散模擬

模型分別在現狀河網和規(guī)劃河網中進行水質模擬。規(guī)劃河網是在現狀河網的基礎上延長航塘港并設置出海水閘。航塘港延長段總長約3.5Km，設定河寬30m，底高程0m，邊坡比1：1.5，航塘港水閘底高程為0m，閘門寬分別假定10m、12m、15m，水頭損失系數0.5。

模型分別在常規(guī)和非常規(guī)調水狀態(tài)下進行模擬。常規(guī)調水時，水閘每次運行時間約3小時，非常規(guī)6小時。本文中常規(guī)調水狀態(tài)為：沿江沿海水閘常規(guī)運行，每日引排1潮次。非常規(guī)調水狀態(tài)為：沿江水閘常規(guī)引水每日2潮次，金南閘常規(guī)排水每日1潮次，其余沿海水閘非常規(guī)排水，每日2潮次。

模型對水質代表點進行模擬，代表點選取鎮(zhèn)內河道，為錢橋社區(qū)的錢塘港、塘外社區(qū)的航塘港和奉城鎮(zhèn)的城河。

圖4 水質代表點所在位置

5.2航塘港水閘的作用分析

假設片區(qū)內某污染物濃度初始值為2mg/L，模擬引清調水7天后得到的污染物濃度變化數據如下表1。

對代表點的污染物濃度值進行分析，可研究航塘港閘對周邊城鎮(zhèn)水質的改善作用。比較現狀河網與規(guī)劃河網的模擬數據，得到以下結果：

（1）常規(guī)調水時，航塘港閘寬10m時只對塘外水質有影響，但影響很小。當閘寬12m和15m時，對3個代表點水質都有影響，影響程度基本相同，分別為：0.7%（錢橋）、1.5%（塘外）、1.5%（奉城）。可見新建航塘港水閘后，在常規(guī)調水時對附近鎮(zhèn)區(qū)河道水質的改善作用較小。endprint

（2）非常規(guī)調水時，航塘港水閘對3個代表點水質都有影響。閘寬10m時，影響程度為：0.9%（錢橋）、12.3%（塘外）、3.7%（奉城）。閘寬12m和15m的影響基本相同，分別為：0.9%（錢橋）、12.8%（塘外）、4.0%（奉城）。可見新建航塘港水閘后，非常規(guī)調水對塘外鎮(zhèn)河道水質改善作用明顯，對奉城鎮(zhèn)水質有一定改善，對錢橋鎮(zhèn)水質改善較小。

（3）航塘港閘對河道水質的改善作用隨距離增加而減少。經過比較，無論常規(guī)還是非常規(guī)調水，航塘港閘寬12m的經濟效用最佳。

6 防汛排澇作用分析

6.1 菲特臺風水情模擬

2013年10月6日，臺風“菲特”登陸，奉賢地區(qū)受此影響開始降雨。6日至8日，三天累積雨量達到194.6mm。為了防澇，南橋水位先緊急降到2.29m，之后水位持續(xù)上漲至降雨結束，8日水位達到最高3.55m。主降雨期，沿江水閘由于黃浦江潮位較高無法排水，只能依靠沿海水閘排水。

模型分別對現狀河網和規(guī)劃河網的內河水位進行模擬，規(guī)劃河網中航塘港水閘閘寬設定為12m。模擬期：2013年10月1日至10月10日，包括了菲特臺風。各水位代表點模擬數據如下表2，水位代表點位置同水質代表點。

菲特臺風降雨峰值出現在8日4：00～10：00，期間沿海水閘排水1潮次，6小時內降雨119mm，查圖6，頻率為11%（9年一遇）。

6.2 50年一遇暴雨水情模擬

查圖6，50年一遇6小時最大降雨178mm。將菲特臺風主降雨期各時段的雨量同比例放大1.5倍，即從119mm放大到178mm。在其它因素條件不變的情況下進行模擬，得到各代表點的模擬水位值，見表3。

表350年一遇暴雨期間各代表點模擬的最高水位表單位：m

6.3 航塘港水閘的作用分析

奉賢區(qū)降雨量無論接近十年一遇還是五十年一遇，新建航塘港水閘后能對塘外鎮(zhèn)的最高水位消減4～5公分，防汛排澇作用明顯。對距離稍遠的錢橋和奉城鎮(zhèn)的最高水位只能消減2公分左右，排澇作用不是很大，對更遠的南橋鎮(zhèn)的最高水位消減作用幾乎為零。

7 結論

通過水量水質耦合模型的模擬，規(guī)劃的航塘港水閘對周邊鎮(zhèn)區(qū)的水質改善和防汛排澇都有一定作用。在引清調水時，水閘運行時間越長，頻次越高，對水質改善的作用越明顯。在常規(guī)調水時，航塘港閘對河道水質的改善作用并不明顯，非常規(guī)調水時對距離最近的塘外鎮(zhèn)河道水質改善作用明顯，次之是奉城鎮(zhèn)，對最遠的錢橋鎮(zhèn)水質改善作用很小。經比較后，無論常規(guī)還是非常規(guī)調水，航塘港閘寬定為12m比較合適。

在確定航塘港閘寬后，經模型模擬，降雨量無論接近十年一遇還是五十年一遇，航塘港閘對塘外鎮(zhèn)的最高水位能消減4到5公分，防汛作用較大，對距離稍遠的錢橋和奉城鎮(zhèn)的最高水位只能消減2公分左右，防汛作用并很大，對更遠的南橋鎮(zhèn)防汛作用幾乎為零。

參考文獻

[1]顧忠華.水文模型在引清調水方案優(yōu)化中的應用[J].上海水務，2016

[2]蔡金傍，等.常州市平原河網水動力模型構建[J].生態(tài)與農村環(huán)境學報，2010

[3]常旭，等.MIKE11模型在洋河流域水質預測中的應用[J].水電能源科學，2013

[4]高程程，等.上海市青松水利片引清調水方案研究[J].水電能源科學，2012

[5]杜鵑.MIKE11在引清調水效果分析中應用[J].東北水利水電，2014

[6]賈瑞華.利用水利工程調度改善水環(huán)境質量的探討[J].上海水利，1989

[7]黃琳煜，等.基于MIKE11的白蓮涇區(qū)域水量水質模型研究[J].水電能源科學，2011

[8]徐祖信，盧士強.平原感潮河網水質模型研究[J].水動力學研究與進展，2013

[9]陳凱，孫鴻萍等.上海市第一次全國水利普查暨第二次水資源普查奉賢區(qū)普查技術報告[M].奉賢，2013.8endprint