阜陽市潁西片區活水工程方案研究

董勝男， 張龍陽

(1.上海市政工程設計研究總院(集團)第六設計院有限公司，安徽 合肥 230031；2.安徽省水文局，安徽 合肥 230031)

1 項目概況

1.1 項目背景

阜陽市位于安徽省西北部，淮河以北，阜陽市區位于市域中部，淮河主要支流潁河自西北向東南流經市區，泉河自西向東注入潁河，將城市自然分割成潁西、潁東、泉北三大分區，城水相依。

近年來，隨著阜陽經濟的飛速發展，水環境問題也日趨惡化。根據檢測資料，阜陽市潁西片區七里長溝、七漁河、西城河、東城河、南城河等河道水質均為劣Ⅴ類，被列為黑臭水體。河道污染嚴重，一方面由于城區河道的雨污合流，導致大量污水未經處理直排入河，另一方面由于地處平原河網區，河道水體流動性差，河道水源補充不足導致河道納污能力下降[1]。隨著雨污分流改造的同時，通過水系溝通，建設活水工程，加強水體流動，可以有效提高河道的納污能力[2，3]，增強河流生態修復功能，從而改善潁西片區河道水質和人居水環境。

1.2 項目范圍

根據潁西片區所在區域及水系情況，本次潁西片區活水工程研究范圍北到泉河，東到潁河，南到蘆橋溝、躍進溝，西到七漁河和人民河，研究面積為201 km2，范圍如圖1所示。

根據潁西片區蓄水工程規劃，工程范圍內劃分了兩個蓄水分區，其中潁西Ⅰ區面積為81 km2，正常蓄水位27.50 m；潁西Ⅱ區面積為120 km2，正常蓄水位27.00 m。

圖1 項目范圍圖

1.3 活水現狀及存在問題

阜陽市泉河、潁河穿城而過，阜陽閘樞紐位于阜陽市三里灣泉河入潁河口處，閘上設計正常蓄水位28.00～28.50 m，閘下最低蓄水位22.00 m，但由于上游農田片區排澇問題，阜陽閘上蓄水位一直未達到正常蓄水位，近年來隨著上游片區河道治理，阜陽閘上正常蓄水位在逐漸抬高。阜陽閘上下游較明顯的水位差為潁西區提供較為有利的自流活水條件。

目前潁西區主要利用西城河涵、七漁河涵自泉河引水，經內部河道，排入阜陽閘下。但由于現狀引水補水點較少，補水水體主要沿流程較短的河道流出，部分河道形成死水區，水體流動性不足，活水效果較差。

2 生態需水量分析

目前，生態需水量計算還處于探索階段，其計算方法眾多，不同方法間計算成果差別較大[4,5]。根據《河湖生態環境需水計算規范》(SL/Z 712-2014)附錄A[6]，生態環境需水量計算通常有Tennant法、Qp法(又稱不同頻率最枯月平均值法)、水量平衡法、7Q10法(又稱最小流量法)等。本次河道生態水量的確定采用Tennant法、保證率分析法和水量平衡法[7]，并按照最基本的生態要求和生態效果非常好的兩種工況進行分析，選取每種工況下不同方法的大值作為本次河道生態需水量(表1)。

表1 河道生態用水量成果分析(單位：萬m3/d)

三種方法在兩種工況下的生態條件需水量的計算成果：為保證河道能夠維持最基本的生態需求，日生態需水量為26.94萬m3，其中潁西Ⅰ區10.27萬m3，潁西Ⅱ區16.67萬m3；為保證生態效果非常好時目標生態環境需水量要求，最大日補水量為100.46萬m3，其中潁西Ⅰ區44萬m3，潁西Ⅱ區56.46萬m3。

3 活水條件分析

3.1 外水水源條件分析

根據區域水系分析，潁西片區活水工程的外水水源為泉河水。根據阜陽閘水文站(1959～2017年)歷年日平均水位資料分析，阜陽閘上歷年日平均水位在27.50m以上天數占總天數的百分比為55%，其中汛期日平均水位在27.50m以上所占百分比為55%，非汛期日平均水位在27.50m以上所占百分比為54%。

但資料顯示，2011年以來，阜陽閘上蓄水位明顯抬高，蓄水位27.50m以上天數年均在300 d以上，其中2013年、2015年全年都在27.50m以上。根據阜陽閘上蓄水容量曲線，阜陽閘上蓄水位在27.50～28.00m間蓄水容積約為1 400萬m3。在上游無來水情況下，可保障潁西區在基本生態環境需求下補水約50 d，在目標生態環境需求下補水約14 d。

3.2 引排水工程條件分析

(1)引水涵閘條件分析。根據調查，潁西片區可利用的現有引水涵閘4座，分別為永豐河閘、七漁河引水閘、西城河引水閘、東城河引水閘，根據現狀涵閘規模計算不同內外水位條件下過流能力，當外河水位高于27.80m時，同時打開各引水閘是能夠滿足不同活水方案的引水規模要求。

(2)排水涵閘條件分析。潁西片區現狀主要出口涵閘主要有6座，包括七里長溝出口閘、華橋溝出口閘、蘆橋溝出口閘、東清河出口閘、中清河出口閘和七漁河出口閘。根據各排水閘現狀涵閘規模計算不同內外水位條件下過流能力，現狀出口閘的總規模能夠滿足活水方案的排水要求。

4 活水工程方案

工程區除引水、排水涵閘外，還有15座水位、流量控制閘，由于涉及河道及閘壩工程較多，調度復雜，本次活水工程方案論證采用MIKE11模型進行分析[8]。

4.1 MIKE11模型簡介

MIKE11水動力計算模型是基于垂向積分的物質和動量守恒方程，即一維非恒定流圣維南Saint-Venant方程組來模擬河流或河口的水流狀態。

式中：x、t分別為計算點空間和時間的坐標，A為過水斷面面積，Q為過流流量，h為水位，q為旁側入流流量，C為謝才系數，R為水力半徑，α為動量校正系數，g為重力加速度。

方程組利用Abbott-Ionescu六點隱式有限差分格式求解。該格式在每一個網格點不同時計算水位和流量，而是按順序交替計算水位或流量，分別稱為h點和Q點。Abbott-Ionescu格式具有穩定性好、計算精度高的特點。離散后的線形方程組用追趕法求解。

4.2 模型參數的選取

MIKE11模型參數文件主要是定義模擬的初始條件和河床糙率。初始條件設定的一個很重要目的是讓模型平穩啟動，所以原則上初始水位和流量的設定應盡可能與模擬開始時刻的河網水動力條件一致，本次設置初始水位同各分區正常蓄水位。

河床糙率根據地質資料及護坡材料，查糙率表選定。根據河道不同斷面及護坡型式，本次各河道糙率選擇范圍為0.017～0.028。

4.3 模型邊界條件

入流邊界：由于本次為活水工程，主要考慮通過泉河水引水，因此本次入流邊界為設置固定流量。

出流邊界：本次出流邊界為對外水系，且為非汛期的自由出流，因此出流邊界均為按照蓄水水位控制要求的自流溢流。

4.4 河網模型概化

MIKE模型的建立關鍵是河網模型的搭建，河網模型搭建合理則模型穩定性好，運算速度快。本次MIKE模型河網搭建過程中，充分考慮模型的穩定性及運算速度，適當對河道水系進行概化，在不影響總體活水方案的前提下，有效改善模型運行條件，概化后河網模型及各閘壩位置如圖2所示。

圖2 活水閘壩位置圖

4.5 活水工程方案

由于區內涉及河道、閘壩較多，不同閘壩之間的調度組合方案多達上千種，本次活水工程提出優化條件下的活水方案，該方案主要利用已建成的閘壩，通過不同閘壩的控制調度，實現河道活水的有效控制。通過分區、各閘壩分功能組合后，將論證方案主要劃分為11大類，51個方案依次進行分析，經過對河道在不同方案中的活水效果分析，確定河道活水的最佳調度方案，并通過調度方案的有效歸集，形成可靠易行的調度方案。

本次活水工程分析河道33條，由于同一河道在不同河段受與其他河道交叉后形成的分流與匯流影響，同一河道又可分為多個河段，此次共分析河段102段。對所有河段分別在不同活水調度方案中的活水流量進行分析，由此確定在引水流量一定的前提下，各河段的最大活水流量。

根據MIKE11模型模擬結果(圖3)，除斷頭河外，其他河段均能進行不同程度的活水，因此該優化條件下的河道活水方案基本是滿足活水要求的，對于不能活水的斷頭河新建斷頭河水體自循環活水系統來改善該河段的活水效果。

圖3 河道活水流量大小示意圖

5 結論及建議

潁西片區水系活水工程通過利用阜陽閘上、閘下的水位差，形成較好的自流活水條件，利用現狀已建涵閘、配套河道連通工程并結合片區水系內部的控制閘調度，實現不同河段的科學合理調度，有效改善河段活水效果，基本能滿足潁西片區水系的活水要求。由于本次活水工程分析中涉及的活水工程調度方案較多，且在一種調度方案中有效活水河段是比較固定的，因此為確保潁西片區各條河道的活水效果，建議按照調度不同規則，進行定期調整。