南京市高淳區(qū)城市水網(wǎng)水環(huán)境改善模擬研究

黃金鳳，宋云浩，董慶華

(1. 長(zhǎng)江勘測(cè)規(guī)劃設(shè)計(jì)研究院，武漢 430010；2. 武漢大學(xué)水利水電學(xué)院，武漢430072)

0 引 言

近年來(lái)，隨著城鎮(zhèn)化進(jìn)程的加快、經(jīng)濟(jì)的快速發(fā)展，導(dǎo)致城市水系現(xiàn)狀十分嚴(yán)峻，如水質(zhì)惡化、生態(tài)功能退化等[1]。于此同時(shí)，也導(dǎo)致河湖水系連通格局與狀況發(fā)生了重大變化，部分地區(qū)產(chǎn)生了連通性減弱、水資源與水環(huán)境承載能力不足、洪水宣泄不暢、水安全風(fēng)險(xiǎn)加大等問題，已經(jīng)成為人水關(guān)系不和諧的重要表現(xiàn)，成為影響經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展、水生態(tài)系統(tǒng)健康的關(guān)鍵制約因素[2-4]。河湖水網(wǎng)作為水資源調(diào)配，水生態(tài)修復(fù)和改善、水災(zāi)害防御的重要手段[5,6]在國(guó)內(nèi)外得到了廣泛應(yīng)用?？茖W(xué)利用水網(wǎng)中的閘門、泵站等控制建筑物調(diào)節(jié)雨洪資源，實(shí)現(xiàn)雨洪資源優(yōu)化利用，是發(fā)揮水網(wǎng)綜合效益的關(guān)鍵[7]。

南京市高淳區(qū)西部水網(wǎng)密布、水系發(fā)育，是首批國(guó)家級(jí)生態(tài)示范區(qū)，然而隨著經(jīng)濟(jì)的發(fā)展，區(qū)域流域水系水質(zhì)惡化，無(wú)法達(dá)到相應(yīng)水功能區(qū)水質(zhì)目標(biāo)，水質(zhì)惡化嚴(yán)重制約了高淳區(qū)生態(tài)文明建設(shè)與經(jīng)濟(jì)發(fā)展。本文章選取南京市高淳區(qū)西部圩區(qū)主要河流為研究對(duì)象，研究區(qū)內(nèi)主要水系連通之后，在枯水年(95%保證率)保障區(qū)域防洪安全、生態(tài)基流的基礎(chǔ)上，雨洪資源最大化利用對(duì)區(qū)域水環(huán)境的改善效果分析。

1 研究區(qū)概況

研究區(qū)域位于高淳區(qū)西部，涉及陽(yáng)江鎮(zhèn)、磚墻鎮(zhèn)、淳溪鎮(zhèn)、古柏鎮(zhèn)、漆橋鎮(zhèn)、椏溪鎮(zhèn)6個(gè)城鎮(zhèn)，水域面積約24.3 km2；高淳區(qū)河湖水系連通涉及固城湖、水陽(yáng)江、水碧橋河、石固河、官溪河、胥河，其中固城湖是高淳區(qū)飲用水水源地(見圖1)。水質(zhì)改善研究?jī)H涉及固城湖、水碧橋河、石固河、官溪河，其中水陽(yáng)江、胥河僅作為邊界調(diào)控閘壩運(yùn)行。

近些年來(lái)，隨著城鎮(zhèn)化進(jìn)程的加速及工業(yè)發(fā)展，研究區(qū)域內(nèi)點(diǎn)源排放量大；研究區(qū)域內(nèi)水資源豐富，養(yǎng)殖業(yè)及水稻種植業(yè)繁榮，面源污染也較嚴(yán)重；大量污染物排入研究水體，導(dǎo)致污染物承接量遠(yuǎn)大于其納污能力。通過分析當(dāng)?shù)貙?shí)測(cè)水質(zhì)數(shù)據(jù)，固城湖水質(zhì)相對(duì)較好，水質(zhì)類別以Ⅱ～Ⅳ類為主；固城湖西圩區(qū)主干河道(水碧橋河、石固河)次之，水質(zhì)類別以Ⅲ～Ⅴ類為主；官溪河水質(zhì)較差，水質(zhì)類別以Ⅳ～Ⅴ類為主；污染因子主要為COD、NH3-N。

圖1 研究范圍示意圖Fig.1 Research scope

2 研究方法

2.1 模型概況

MIKE系列模型結(jié)構(gòu)清晰、界面友好、考慮涉水要素全面，在河口、河流、河湖水網(wǎng)模型構(gòu)建方面應(yīng)用廣泛[8,9]，本文章采用由丹麥水利研究所(簡(jiǎn)稱DHI，下同)獨(dú)立開發(fā)的MIKE11/21/FLOOD來(lái)搭建高淳區(qū)水碧橋河、官溪河、石固河、水陽(yáng)江、胥河、固城湖水動(dòng)力、水質(zhì)(對(duì)流擴(kuò)散)模型。MIKE11/21/FLOOD包含水動(dòng)力、對(duì)流擴(kuò)散、水環(huán)境等多個(gè)模塊，具有計(jì)算穩(wěn)定、精度高、可靠性強(qiáng)等特點(diǎn)，能方便靈活地模擬復(fù)雜河網(wǎng)的水流特性，模擬閘門、水泵等水工建筑物的運(yùn)營(yíng)調(diào)度，是目前世界上應(yīng)用最為廣泛的軟件之一，其中水質(zhì)(對(duì)流擴(kuò)散)模塊基于水動(dòng)力模塊建立。

研究區(qū)域內(nèi)河流均采用一維水動(dòng)力水質(zhì)模型；固城湖湖面面積較大，平均水深約在3 m左右，屬于寬淺型湖泊，采用平面二維水動(dòng)力水質(zhì)模型；最終通過MIKE FLOOD進(jìn)行耦合計(jì)算。

2.2 模型原理

(1)水動(dòng)力模塊(HD)。一維水動(dòng)力模型方程組[10]為：

(1)

式中：x為笛卡爾坐標(biāo)系坐標(biāo)，m；g為重力加速度，m/s2；t為時(shí)間，s；A為過水?dāng)嗝婷娣e，m2；R水力半徑，m；C為謝才系數(shù)，m1/2/s；Q為流量，m3/s；q為旁側(cè)入流量，m3/(s/m)；h為水位，m;α為動(dòng)量修正系數(shù)。

上式中第一個(gè)方程為質(zhì)量守恒方程，第二個(gè)方程為動(dòng)量守恒方程。

二維非恒定淺水運(yùn)動(dòng)方程組[11]為：

(2)

式中：t為時(shí)間,s；x、y為笛卡爾坐標(biāo)系坐標(biāo)，m；h水深為，m；η為水位，m；d=η-h，m；u、ν分別為x、y方向上的速度分量，m/s；g為重力加速度，m/s2；ρ為水的密度，kg/m3;C為謝才系數(shù)，m1/2/s；τxx、τxy、τyy分別有效剪應(yīng)力分量，N/m2；S為源項(xiàng)，m3/(s/m2)；Ω為科氏力系數(shù)，s-1；f(V)是風(fēng)摩擦力系數(shù)；V、Vx、Vy分別為風(fēng)速及風(fēng)速分量，m/s；Pa為大氣壓，kg/(m/s2)；p=hu，q=hν。

(2)對(duì)流擴(kuò)散模塊(AD)。

一維對(duì)流擴(kuò)散方程為[12]：

(3)

二維對(duì)流擴(kuò)散方程為[13]:

(4)

式中：C為污染物濃度，mg/L;Q為河道流量，m3/s;A為過水?dāng)嗝婷娣e，m2;t為時(shí)間，s;x、y為笛卡爾坐標(biāo)系坐標(biāo)，m；Dx、Dy為河道x、y方向上的污染物擴(kuò)散系數(shù)，m2/s，K為污染物的衰減系數(shù)，s-1；C2為旁側(cè)入流的污染物的濃度，mg/L；q為旁側(cè)入流，m3/(s·m);u、ν分別為x、y方向的流速，m/s;h為水深，m;QS為源項(xiàng)入流，m3/(s·m2)；CS為源項(xiàng)濃度，mg/L。

3 模型構(gòu)建

(1)水動(dòng)力模塊(HD)。模型概化：收集研究區(qū)域的湖泊地形、河道斷面等數(shù)據(jù)，利用MIKE生成河網(wǎng)文件(河網(wǎng)連接見圖1，固城湖網(wǎng)格劃分見圖2)。

邊界條件：根據(jù)實(shí)測(cè)數(shù)據(jù)，主要參數(shù)設(shè)定如下：固城湖初始水位9.5 m；水碧橋河、官溪河、石固河與固城湖在模型中為標(biāo)準(zhǔn)連接；糙率取0.03；河床阻力取0.32 m1/3/s；河流上游采用95%保證率枯水年下的日流量過程，下游采用95%保證率下的日水位過程。

控制建筑物及調(diào)度規(guī)則：根據(jù)研究需要，本模型主要設(shè)定的控制建筑物包括水碧橋閘、楊家灣閘、楊家灣泵站、蛇山泵站；蛇山閘、茅東閘在模型中一直處于關(guān)閉狀態(tài)。本次研究主要包括兩種工況(見表1)，雨洪資源最大化利用工況是指在枯水期保證生態(tài)基流、豐水期保證區(qū)域防洪安全的前提下，科學(xué)利用水網(wǎng)中的閘門、泵站等控制建筑物，調(diào)節(jié)雨洪資源盡可能多地補(bǔ)給區(qū)域生產(chǎn)、生活用水，發(fā)揮水網(wǎng)綜合效益。

圖2 固城湖網(wǎng)格劃分示意圖Fig.2 mesh generation of Gucheng Lake(right)

圖3 水碧橋河、石固河、官溪河代表斷面示意圖Fig.3 Typical cross-sections of Shuibiqiao River、Shigu River、Guanxi River

表1 控制建筑物調(diào)度規(guī)則設(shè)定Tab.1 Regulation principles of control structures

河流生態(tài)基流根據(jù)現(xiàn)狀污染物入河情況，依據(jù)《水域納污能力計(jì)算規(guī)程 SL348-2006》中的河流一維模型計(jì)算90%保證率下，為保證研究區(qū)域內(nèi)河流COD、NH3-N達(dá)到《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB 3838-2002)》規(guī)定的地表水Ⅲ類水質(zhì)標(biāo)準(zhǔn)所需要的基礎(chǔ)流量，計(jì)算結(jié)果見表2。蘭林等[14]研究提出固城湖湖泊需控制水位7.0 m以上，以防止出現(xiàn)生態(tài)災(zāi)難；固城湖是高淳區(qū)的飲用水水源地，取水能力10 萬(wàn)m3/d，取水口位于固城湖大湖區(qū)紅砂嘴處，設(shè)計(jì)最低取水水位7.03 m；綜合分析，固城湖最低水位應(yīng)為取水口最低取水水位7.03 m。

表2 河流生態(tài)基流表Tab.2 Ecological basic flow of rivers in study area

(2)對(duì)流擴(kuò)散模塊(AD)。對(duì)流擴(kuò)散模型在模擬過程中需要在水動(dòng)力模型的基礎(chǔ)上準(zhǔn)備以下文件:邊界文件、模型參數(shù)文件。研究區(qū)域內(nèi)水體主要超標(biāo)因子為COD、NH3-N，因此本次水質(zhì)分析主要對(duì)COD、NH3-N進(jìn)行模擬分析；其污染主要為外源污染。

邊界文件：主要包括源匯項(xiàng)和開邊界的水質(zhì)濃度，模型的計(jì)算以參考《高淳區(qū)水生態(tài)文明建設(shè)實(shí)施方案》，考慮社會(huì)經(jīng)濟(jì)發(fā)展、截污減排之后,在保證水體基礎(chǔ)流量、防洪安全、雨洪資源最大化利用的基礎(chǔ)上，分析河流湖泊水環(huán)境的改善效果分析。根據(jù)分析，污染物入河量如下表3，排污口設(shè)置見圖1。

表3 計(jì)算斷面污染物入河量Tab.3 Inflow of pollutants

參數(shù)文件：主要包括模擬污染物的基本信息，初始濃度，對(duì)流系數(shù)，擴(kuò)散系數(shù)，附加輸出結(jié)果等。其中綜合衰減系數(shù)及污染物擴(kuò)散系數(shù)參考河海大學(xué)逄勇教授關(guān)于江蘇省水環(huán)境容量計(jì)算的實(shí)驗(yàn)值[15]，結(jié)合《高淳區(qū)水生態(tài)文明建設(shè)實(shí)施方案》中的值確定，污染物綜合衰減系數(shù)COD為0.15 d-1，氨氮為0.05 d-1；擴(kuò)散系數(shù)為20 m2/s。上下游邊界處水質(zhì)濃度主要根據(jù)水質(zhì)監(jiān)測(cè)數(shù)據(jù)輸入。

4 結(jié)果與分析

通過《高淳區(qū)水生態(tài)文明建設(shè)實(shí)施方案》，水質(zhì)改善涉及水域水功能區(qū)水質(zhì)目標(biāo)如表4所示。

(1) 現(xiàn)狀實(shí)際調(diào)度工況。河流：通過模擬分析，三條河流全河段COD全年滿足Ⅲ類水水質(zhì)目標(biāo)要求，三條河流入湖口COD全年可達(dá)到Ⅱ類水濃度，具體如圖4所示。石固河排污口以下3 km范圍內(nèi)，有33%的時(shí)間NH3-N僅達(dá)到Ⅳ類水濃度；官溪河排污口以下1 km范圍內(nèi)，有24%的時(shí)間NH3-N僅達(dá)到Ⅳ類水濃度；三條河流其他河段、其他時(shí)間NH3-N均可滿足Ⅲ類水水質(zhì)目標(biāo)要求；官溪河、水碧橋河入固城湖湖口NH3-N全年滿足Ⅲ類水水質(zhì)目標(biāo)要求，石固河入湖口NH3-N達(dá)標(biāo)率為98%；具體如圖5所示。

表4 研究區(qū)域水功能區(qū)水質(zhì)目標(biāo)表Tab.4 Water quality objectives of water functional areas in study area

圖4 現(xiàn)狀實(shí)際調(diào)度工況河流入湖口斷面 COD濃度分布(1986年)Fig.4 Concentration of COD at the entrance cross-section of Gucheng Lake under current regulation principle

圖5 現(xiàn)狀實(shí)際調(diào)度工況河流入湖口斷面 NH3-N濃度分布(1986年)Fig.5 Concentration of NH3-N at the entrance cross-section of Gucheng Lake under current regulation principle

固城湖：通過模擬分析，固城湖COD全年幾乎全湖泊可滿足Ⅱ類水水質(zhì)目標(biāo)要求，NH3-N全湖大面積水域可滿足Ⅱ類水水質(zhì)目標(biāo)要求，最不利分布時(shí)刻分別見圖6。

圖6 現(xiàn)狀實(shí)際調(diào)度工況下湖泊COD、NH3-N濃度最不利時(shí)刻分布Fig.6 Concentration of COD and NH3-N in Gucheng Lake at the worst moment under current regulation principle

(2) 雨洪資源最大化利用工況。河流：通過模擬分析，三條河流全河段COD全年滿足Ⅲ類水水質(zhì)目標(biāo)要求，三條河流入湖口COD全年可達(dá)到Ⅱ類水濃度，具體如圖7所示。石固河排污口以下2 km范圍內(nèi)，有21%的時(shí)間NH3-N僅達(dá)到Ⅳ類水濃度；官溪河排污口以下1 km范圍內(nèi)，有15%的時(shí)間NH3-N僅達(dá)到Ⅳ類水濃度；其他河段全年水質(zhì)均可滿足Ⅲ類水水質(zhì)目標(biāo)要求；官溪河排污口以下1 km范圍內(nèi)，有15%的時(shí)間NH3-N僅滿足Ⅳ類水水質(zhì)目標(biāo)要求；三條河流其他河段、其他時(shí)間NH3-N均可滿足Ⅲ類水水質(zhì)目標(biāo)要求；官溪河、水碧橋河入固城湖湖口NH3-N全年滿足Ⅲ類水水質(zhì)目標(biāo)要求，石固河入湖口NH3-N達(dá)標(biāo)率為98%；具體如圖8所示。

圖7 雨洪資源最大化利用工況下河流入湖口斷面 COD濃度分布(1986年)Fig.7 Concentration of COD at the entrance cross-section of Gucheng Lake under optimized regulation principle

圖8 雨洪資源最大化利用工況下河流入湖口斷面 NH3-N濃度分布(1986年)Fig.8 Concentration of NH3-N at the entrance cross-section of Gucheng Lake under optimized regulation principle

固城湖：固城湖COD全年幾乎全湖泊可滿足Ⅱ類水水質(zhì)目標(biāo)要求，NH3-N全湖大面積水域可滿足Ⅱ類水水質(zhì)目標(biāo)要求，最不利分布時(shí)刻分別見圖9。

圖9 雨洪資源最大化利用工況下湖泊 COD、NH3-N濃度最不利時(shí)刻分布Fig.9 Concentration of COD(left) and NH3-N(right) in Gucheng Lake at the worst moment under optimized regulation principle

(3) 結(jié)論分析。通過模擬分析，兩種工況下，三條河流全河段COD、水碧橋河全河段NH3-N均能達(dá)到Ⅲ類水水質(zhì)目標(biāo)。但是官溪河、石固河排污口下游3 km河段內(nèi)，雨洪資源最大化利用工況下，河段水質(zhì)達(dá)到地表水Ⅲ類水水質(zhì)要求的比例均大于現(xiàn)狀實(shí)際調(diào)度工況；河段水質(zhì)不達(dá)標(biāo)比例均小于現(xiàn)狀實(shí)際調(diào)度工況，具體分析見圖10、圖11。兩種工況下，固城湖水質(zhì)均能滿足地表水Ⅱ類水水質(zhì)要求，通過對(duì)比，雨洪資源最大化利用工況相比于現(xiàn)狀實(shí)際調(diào)度工況而言，水質(zhì)相對(duì)更好的湖域面積更多。

圖10 不同工況官溪河不同河段水質(zhì)不達(dá)標(biāo)比例Fig.10 The ratio of substandard water quality along Guanxi River under current and optimized regulation principle

圖11 不同工況石固河不同河段水質(zhì)不達(dá)標(biāo)比例Fig.11 The ratio of substandard water quality along Shigu River under current and optimized regulation principle

5 結(jié) 論

通過建立水碧橋河、官溪河、石固河、固城湖對(duì)流擴(kuò)散模型分析，相比現(xiàn)狀閘壩調(diào)度工況，雨洪資源最大化利用工況不僅可以優(yōu)化當(dāng)?shù)厮Y源調(diào)配，更好利用雨洪資源滿足區(qū)域用水需求，還可通過提高河湖枯期水量，提高河湖枯期納污能力，從而提高區(qū)域水質(zhì)。綜合來(lái)看，在本研究區(qū)域內(nèi)的，河湖水網(wǎng)可通過水量調(diào)配，改善高淳區(qū)水體水質(zhì)。

□