數(shù)值模擬在嘉定蘊南中部地區(qū)引調(diào)水試運行準備工作中的應(yīng)用

包小燕，李念斌

（1.上海市嘉定區(qū)水文站，上海 201800；2.上海嘉定水務(wù)工程設(shè)計有限公司，上海 201800）

0 引言

上海市嘉定區(qū)蘊南中部地區(qū)涉及南翔、江橋、黃渡三鎮(zhèn)，是嘉定區(qū)經(jīng)濟最發(fā)達的地區(qū)之一，但也是嘉定區(qū)水環(huán)境質(zhì)量最差的地區(qū)。為此嘉定區(qū)水務(wù)局于2007年組織編制了《嘉定蘊南中部地區(qū)水環(huán)境治理——引調(diào)水專項規(guī)劃》，經(jīng)過2008年至2012年的建設(shè)，規(guī)劃中的水利控制工程已基本建設(shè)完成，嘉定區(qū)水務(wù)局擬定于2012年10月進行調(diào)水試運行。針對規(guī)劃中設(shè)置的調(diào)水試運行方案，使用數(shù)值模擬的方法進行運算對比，選取合理的水質(zhì)分析比對點，為調(diào)水試運行提供參考依據(jù)。

1 河網(wǎng)水動力數(shù)值模擬發(fā)展現(xiàn)狀

數(shù)值模擬也叫計算機模擬。它以電子計算機為手段，通過數(shù)值計算和圖像顯示的方法，達到對工程問題和物理問題乃至自然界各類問題研究的目的。20世紀50年代以前，數(shù)值模擬的基本理論已經(jīng)建立，運用這些理論也曾解決過一些簡單的工程問題。但是，這些基本理論真正應(yīng)用于解決工程問題卻是在電子計算機發(fā)明以后。1952年—1954 年，E.lsaacson、J.J.stoker和 B.A.Troesch 首次建立了俄亥俄河和密西西比河部分河段的數(shù)學模型，并進行了實際洪水過程的模擬。到上世紀60年代中期，為了解決各種各樣的設(shè)計和規(guī)劃問題，水動力數(shù)學模型再次得到重視。隨著計算機技術(shù)的發(fā)展，水動力數(shù)學模型的模擬功能也大大增強，可以對整個流域、洪泛區(qū)、已建或規(guī)劃中的水利工程進行系統(tǒng)模擬。水動力模型是其它模型(如水質(zhì)模型、水環(huán)境容量模型等)的基礎(chǔ)和前提。水動力模型的建立涉及控制方程組的簡化、方程組的離散和求解、初邊條件的確定、模型的率定和驗證等一系列問題。平原河網(wǎng)模型不同于單一河流的特點，在于河網(wǎng)錯綜復(fù)雜性，以及由此帶來的方程組離散和求解上的困難，這是多年來人們研究河網(wǎng)問題的一大難點。平原河網(wǎng)水動力數(shù)學模型目前大體可以分為節(jié)點—河道模型、單元劃分模型、混合模型、人工神經(jīng)網(wǎng)絡(luò)模型，以及蒙特卡羅隨機游動模型，其中最常用的是節(jié)點—河道模型。

2 河網(wǎng)非恒定流數(shù)值解法

2.1 水動力方程

2.1.1 基本方程組

水流在平底、棱柱形明渠中一維非恒定流動的基本方程組—圣維南方程組,包括反映動量守恒定律的運動方程和反映質(zhì)量守恒定律的連續(xù)性方程：

式中，A為河道內(nèi)任意斷面的面積，m2；Q為斷面水流量，m3；t、x分別為時間與距離的坐標；Z為水位，m；g為重力加速度，m·s-2；q 為旁側(cè)入流，m3；Rh為水力半徑，m；C k 為謝才系數(shù)，m1/2·s-1；a 為動量校正系數(shù)。

2.1.2 方程組求解方法

如圖1所示,利用Abbott六點隱式差分格式求解圣維南方程,在每個網(wǎng)格節(jié)點按順序交替計算水位和流量。該格式為無條件穩(wěn)定,可以在相當大的Courant數(shù)下保持計算穩(wěn)定,因此,可以取較長的時間步長以節(jié)省模型運算時間。

圖1 Abott格式水位點、流量點布置示意圖

2.1.3 定解條件

定解條件為水流的初值與邊界值,具體如下:水流初始條件：t=0，Z(x,t)=Z(x,0)，Q(x,t)=Q(x,0)。邊界條件：當 x=0 時，Z(x,t)=Z(0,t)；當 x=L 時，Z(x,t)=Z(L,t)。

河網(wǎng)的內(nèi)邊界條件比較復(fù)雜，其中河網(wǎng)支汊節(jié)點的處理最為關(guān)鍵，支汊點應(yīng)滿足Z1=Z2=……

2.2 水質(zhì)擴散方程

根據(jù)物質(zhì)在水體中對流擴散的過程，考慮污染物的對流擴散與線性降解，其基本方程為[1-2]：

式中:C為污染物濃度，mg/L；D為污染物彌散系數(shù)；A 為斷面過水面積，m2；Q 為流量，m3/s；K 為降解系數(shù)，s-1；C2為污染物的點源濃度，mg/L；q 為污染物的點源流量，m3/s；x為空間步長，m；t為時間步長，s。

3 區(qū)域概況及方案設(shè)置

3.1 區(qū)域情況

嘉定區(qū)位于上海市的西北部，市級河道蘊藻浜橫貫其南部，在穿越滬寧鐵路后折向西南與蘇州河交匯。由此，嘉定區(qū)水系形成了界線分明的南北兩片。其中，蘊藻浜以南地區(qū)的面積約為106.4 km2，其范圍是：東以嘉定行政區(qū)域為界；西以蘊藻浜為界；南北兩側(cè)各以蘇州河和蘊藻浜為界，形成了一個相對獨立的區(qū)域。這一區(qū)域俗稱嘉定的蘊南地區(qū)。

3.2 水環(huán)境現(xiàn)狀

該區(qū)屬平原感潮河網(wǎng)地區(qū)，內(nèi)部水系絕大部分天然形成，且彎曲流長；河道水流流向不定，主要受南部感潮河流蘇州河的水位影響，水流方向由北向南、由西向東居多。但由于區(qū)內(nèi)河道向蘇州河排水受到限制，阻斷了區(qū)域內(nèi)外水體自然交換的道路，改變了地區(qū)原先“由北向南流”的水力特性。此外，由于鐵路、公路的建設(shè)，交叉建筑物也阻礙了水流的暢通，使該地區(qū)基本成為死水區(qū)，水質(zhì)情況常年處于劣Ⅴ類水。

3.3 方案設(shè)置

根據(jù)《嘉定蘊南中部地區(qū)水環(huán)境治理——引調(diào)水專項規(guī)劃》的方案規(guī)劃，將通過沿蘊藻浜南側(cè)沿線設(shè)站，提抽內(nèi)部河道水進入蘊藻浜，待內(nèi)部河水位降低后，強行引入蘇州河清水實行水體交換與流動。并利用東側(cè)現(xiàn)有的連浦閘站向內(nèi)引水，走馬塘水閘開啟。這樣形成“北排、南引”的有序水流大格局。

該方案沿蘊藻浜在淺江、封浜、橫瀝設(shè)泵閘，淺江泵閘流量為5 m3/s、封浜泵閘流量為10 m3/s、橫瀝泵閘流量為10 m3/s。為形成良好的換水條件，與該方案配套的工程有：西線沿鹽鐵塘在狗頸骨河、北線沿蘊藻浜在花園江處設(shè)閘控制。沿蘊藻浜控制工程如圖2所示。

圖2 蘊南中部地區(qū)引調(diào)水控制工程分布圖

4 數(shù)值模擬方法

4.1 河網(wǎng)概化

在充分掌握詳細的天然河網(wǎng)、湖泊的水動力、水文資料的基礎(chǔ)上，以主干河道為基礎(chǔ)，按照河網(wǎng)概化的基本原則，對蘊南中部地區(qū)河網(wǎng)進行合理地概化。所涉及的主要河流有:吳淞江、蘊藻浜、鹽鐵塘、淺江、封浜、橫瀝、西虬江、中槎浦等。

4.2 邊界條件

4.2.1 外邊界條件

平原河網(wǎng)地區(qū)在洪水期往往同時遭遇平原區(qū)產(chǎn)流和潮汐等多重影響，其邊界條件大致分為3類：水位（潮位）邊界、流量邊界（水位流量關(guān)系邊界可以轉(zhuǎn)化為流量邊界），以及計算區(qū)域的凈雨輸入。該項研究采用吳淞江黃渡站水位及蘊藻浜閘內(nèi)站水位為外邊界。

4.2.2 內(nèi)邊界條件

在內(nèi)部邊界條件處，圣維南方程組不再適用，必須依據(jù)其水力特性作特殊處理，以保證河道追趕系數(shù)的連續(xù)傳遞。

4.3 時間空間步長的選取

考慮到模型計算穩(wěn)定性要求及模型運行時間要求，水動力模型的時間步長可以取5～15 min，空間步長取50 m～1500 m不等。

4.4 水質(zhì)分析點選取

該項計算主要選取封浜、橫瀝、西虬江的汊點斷面處，由于封浜、西虬江被鐵路阻隔，在封浜被鐵路阻隔的南北、西虬江被鐵路阻隔的東西分別選取水質(zhì)分析點，共計20個，詳見表1所列。

表1 水質(zhì)分析點一覽表

4.5 模擬結(jié)果

由于溶解氧、氨氮等水質(zhì)指標會沿程和隨時間變化，所以采用綜合的CODcr作為計算指標，用于換水效果的比較。限于篇幅，以橫瀝河為例，72 hCOD相對減少量的平均值（%）見表2所列，隨時間變化過程曲線如圖3所示。

表2 72 hCOD相對減少量的平均值一覽表（單位：%）

圖3 橫瀝選取斷面72 hCOD相對減少量隨時間變化過程曲線圖

5 建議

在蘊南中部地區(qū)引調(diào)水試運行前，針對一維河網(wǎng)非恒定流復(fù)雜的水流特性，建立了嘉定蘊南中部地區(qū)的水量水質(zhì)模型，選取了水質(zhì)分析點，并進行了數(shù)值計算，為試運行提供了參考依據(jù)。現(xiàn)針對將進行的試運行，提出以下建議：

（1）做好水質(zhì)取樣分析工作，掌握蘊南中部地區(qū)的水質(zhì)在引調(diào)水前、后的情況，尤其是數(shù)值計算涉及的這些分析點，為模型的校準提供數(shù)據(jù)。

（2）及時采集引調(diào)水的數(shù)據(jù)，分析試運行中存在的問題，為蘊南中部地區(qū)引調(diào)水常態(tài)化做好準備。

（3）在已建立的數(shù)值模型基礎(chǔ)上，通過試運行的數(shù)據(jù)，對此區(qū)域的水利控制工程調(diào)度問題開展探索。

[1]歐 劍,馬進榮,張行南,張曉艷.大通至長江口整體水動力模型[J].河海大學學報:自然科學版，2009,(5):258-262.

[2]匡翠萍,李行偉,劉曙光.大規(guī)模圍墾對香港維多利亞港水動力環(huán)流的影響 [J].同濟大學學報 (自然科學版)，2009，（2）：176-181.

[3]郜會彩,李義天,何用,陳建.改善漢陽湖群水環(huán)境的調(diào)水方案研究[J].水資源保護，2006,(9):41-44.

[4]黃偉,何用,李榮,陳鋼鋒.感潮地區(qū)引清調(diào)水方案研究[J].水電能源科學，2004,(6):19-22.