楊朝軍，應澤林，陳 飛

(1.紹興市越城區(qū)河道水域管理中心，浙江 紹興 312000；2.紹興市水利工程管理中心，浙江 紹興 312000；3.紹興市水利水電勘測設計院有限公司，浙江 紹興 312000)

0 引 言

紹興城區(qū)位于紹興平原，包括越城區(qū)、柯橋區(qū)及上虞區(qū)的一部分區(qū)域。紹興平原河網(wǎng)河流眾多，水系密布，素有“東方威尼斯”之稱。紹興平原河網(wǎng)地區(qū)人口密集，社會經(jīng)濟發(fā)達，工業(yè)企業(yè)和服務業(yè)聚集，污染源較多，并且由于地勢平坦，平原河網(wǎng)水動力條件較差，存在著往復流和滯流等現(xiàn)象，使得水體抗污染能力較弱，水環(huán)境整體質(zhì)量較差[1]。引水工程是改善平原河網(wǎng)水質(zhì)較為經(jīng)濟、可行性較高的技術方法之一[2]，主要通過稀釋水體，提升水體流動性以增強水體自凈能力的方式來提高水環(huán)境質(zhì)量，已在我國海河、長江等流域中得到了廣泛應用。

由于平原河網(wǎng)復雜的水系、水環(huán)境狀況，引調(diào)水最終產(chǎn)生的實際效果往往受到多重因素影響，引水水量是其中一個重要因素。如何科學有效地根據(jù)特定目標情況，按照經(jīng)濟合理的原則，選擇能夠達到目標水質(zhì)所需要的最小引水流量，即“最優(yōu)引水流量”，至關重要。

本研究以紹興平原河網(wǎng)區(qū)為研究對象，以一維水動力模型為基礎，利用水環(huán)境數(shù)值模擬技術，研究不同引水流量對區(qū)域水環(huán)境的改善效果，根據(jù)不同目標組合對引調(diào)水方案進行精準評估，以期為平原河網(wǎng)地區(qū)水環(huán)境綜合治理提供有效支持。

1 模型基本方程

本次水動力計算采用mike11中的HD模塊[3-5]，本模型演進計算的對象是平原河網(wǎng)，河道交錯，水流流向復雜。計算采用一維非恒定流方法，圣維南偏微分方程組為：

(1)

(2)

式中：x—計算點的空間坐標；t—計算點的時間坐標；A—過水斷面面積,m2;Q—過流流量,m3/s;h—水位,m;q—旁側(cè)入流流量,m3/s;C—謝才系數(shù),m1/2/s;R—水力半徑,m;α—動量校正系數(shù)，g—重力加速度,m/s2。

方程組利用Abbott-ionescu六點隱式有限差分格式求解圣維南方程(見圖1)。該格式在每一個網(wǎng)格點按順序交替計算水位或流量，分別稱為h點和Q點。Abbott-ionescu格式具有穩(wěn)定性好、計算精度高的特點。離散后的線形方程組用追趕法求解。

圖1 Abbott格式水位點、流量點布置圖

2 水文和水動力模型構建

2.1 水文和水動力模型

紹興平原河網(wǎng)模型概化了“六縱五橫”骨干河道及其間連接的重要河道共計200余條，共布設約3 900多個計算斷面，包含有重要湖泊8個，排澇閘、節(jié)制閘12個，鐵路公路涵洞50個，規(guī)劃排澇泵站3個。模型基本上涵蓋紹興平原區(qū)主要河道、全部的排水擋潮閘以及曹娥江上浦閘以下河段，綜合了現(xiàn)有河網(wǎng)、水閘和規(guī)劃河道、泵站、排水擋潮閘的調(diào)蓄與排澇的作用。

紹興平原河網(wǎng)水利計算概化圖見圖2。

圖2 紹興平原河網(wǎng)水利計算概化圖

2.2 水質(zhì)模型

Mike11 ecolab模型是一個開放性的通用工具，可用來制定水生態(tài)模型，模擬水質(zhì)、富營養(yǎng)化和生態(tài)狀況等。水質(zhì)模型所需要的輸入數(shù)據(jù)包括流域近年來水質(zhì)監(jiān)測數(shù)據(jù)，點源和非點源污染負荷數(shù)據(jù)等。本次研究的污染物主要包括COD、NH4-N等，影響河流水體中生物需氧量的主要過程為生物降解，其平衡方程為：

式中:CB—水中BOD5的濃度;Kd—BOD5的降解系數(shù);Θd—BOD5降解過程溫度系數(shù);T—水溫;Ks—半飽和常數(shù);KD—溶解氧的降解系數(shù)。

3 紹興平原河網(wǎng)的引水工程

3.1 曹娥江引水工程

紹興城區(qū)曹娥江引水工程是一項將曹娥江水引至紹興市區(qū)，改善紹興市城區(qū)水環(huán)境的綜合性工程。整個工程東西向橫穿上虞區(qū)、柯橋區(qū)和越城區(qū)，至紹興市城區(qū)，再向北匯流入曹娥江。該工程于2010年底正式通水運行。

3.2 浙東引水工程

浙東引水工程的主要任務是引富春江水向蕭紹寧平原及舟山市補充工業(yè)和農(nóng)灌等一般用水，并兼顧改善水環(huán)境，設計引水流量為50 m3/s。

4 河道水質(zhì)情況分析

4.1 引水工程對水質(zhì)改善的影響分析

引水工程對城區(qū)水質(zhì)改善的效果非常明顯。受水區(qū)范圍內(nèi)的河道監(jiān)測水質(zhì)指標明顯低于外部區(qū)域。受水區(qū)水質(zhì)指標主要受曹娥江引水工程、浙東引水工程、以及南部山區(qū)徑流影響。以實測高錳酸鹽指數(shù)為例，從圖3中可以看出，曹娥江引水對區(qū)域內(nèi)水質(zhì)高錳酸鹽指數(shù)有明顯的改善效果，浙東引水對水體的改善主要集中在柯橋西北部，對環(huán)城河區(qū)域影響較弱，曹娥江引水則對南環(huán)河等越城區(qū)周邊河道產(chǎn)生明顯影響。

圖3 曹娥江引水工程影響范圍示意圖

4.2 受益區(qū)河道水質(zhì)的月際變化規(guī)律

城區(qū)河道監(jiān)測站點的水質(zhì)監(jiān)測數(shù)據(jù)存在明顯的月際變化規(guī)律[6]?？傮w上，各項水質(zhì)指標自每年在1—6月份逐漸升高，并自9月份開始至年底逐漸降低，表現(xiàn)為夏季濃度高，冬季濃度低。圖4為2018—2020年間高錳酸鹽指數(shù)和總磷指標數(shù)值隨著月度實測數(shù)據(jù)變化的規(guī)律。

4.3 河道水質(zhì)變化與降雨關系分析

承前所述，城區(qū)河道水質(zhì)監(jiān)測數(shù)據(jù)的月際變化規(guī)律與紹興地區(qū)降雨的月際變化規(guī)律存在很高的關聯(lián)性，即隨著降雨量的增加，河道水質(zhì)指標數(shù)數(shù)值會相應增加；反之，河道水質(zhì)指標值相應減少。圖5所示為高錳酸鹽指數(shù)和總磷指標數(shù)值與降雨量的關系，其中，水質(zhì)指標同上文，降雨量為紹興市平均月降雨量。

圖5 高錳酸鹽指數(shù)和總磷指標數(shù)值與降雨量的關系

4.4 不同引水流量對受水區(qū)的影響分析

引水流量與水質(zhì)指標的相應關系采用數(shù)學模型法進行研究。對河道情況和邊界條件進行合理概化，建立紹興平原河網(wǎng)水質(zhì)模型，分析曹娥江引水對受水區(qū)水質(zhì)改善范圍的影響。采用年平均徑流量和負荷量作為邊界，以環(huán)水質(zhì)指標數(shù)值高錳酸鹽指數(shù)4.0 mg/L，氨氮0.5 mg/L，總磷0.12 mg/L為依據(jù)，引水工程的影響范圍結(jié)果見圖6。從模型計算結(jié)果可以分析，當引水量在10～15 m3/s之間，受益范圍主要包括古城內(nèi)河、環(huán)城河及南部的南環(huán)河、南池江、坡塘江河道，受益河道水面面積1.05 m2，受益水量約3 600 000 m3；當引水量為5 m3/s時，受益范圍主要包括南環(huán)河(坡塘江以東)、南池江、禹陵江、平水西江，而南環(huán)河以北至馬山閘西江其余河道則處于浙東引水工程與曹娥江引水工程的影響范圍之間，水體流動緩慢，置換周期較長。

圖6 曹娥江引水工程影響范圍示意圖

4.5 階段性停引對城區(qū)水質(zhì)的影響分析

引水工程階段性停引對城區(qū)河道水質(zhì)影響明顯。以圖7主要代表站點南渡橋、渡東橋、西郭鐵路橋、葛山頭、偏門公路橋的實測數(shù)據(jù)為分析對象，當工程停引之后，河道水質(zhì)指標會明顯提高；停引一段時間后再啟動引水工程，河道水質(zhì)指標會明顯下降。

當曹娥江引水工程持續(xù)30 d引水后，高錳酸鹽指數(shù)、氨氮和總磷三項水質(zhì)指標的平均數(shù)值較停引期間分別低32%、70%、65%，指標值可達到地表水Ⅱ類水標準。當引水工程停止引水超過30 d時，高錳酸鹽指數(shù)、氨氮、總磷三項水質(zhì)指標的平均數(shù)值會逐步上升?；謴鸵?，根據(jù)引水流量大小，水質(zhì)指標數(shù)值有不同程度的降低。圖8所示為氨氮指標數(shù)值變化與引水量關系。

圖8 氨氮指標數(shù)值與引水天數(shù)關系

5 實現(xiàn)紹興平原河網(wǎng)精確引水方案研究

5.1 精準引水方案總體目標

受益區(qū)主要代表站點的高錳酸鹽指數(shù)、氨氮、總磷三項水質(zhì)指標數(shù)值分別低于4.0 mg/L、0.5 mg/L、0.12 mg/L。減少水質(zhì)指標的月際波動，維持水質(zhì)相對穩(wěn)定。以維持受益水體最低生態(tài)流量5.0 m3/s為引水下限,以設計最大引水流量15 m3/s為上限，綜合考慮梅汛期、臺汛期調(diào)度情況，采用模型試算法確定引水量[7-8]。

(1)引水水量調(diào)度。每年11月至次年2月，日平均引水流量為5 m3/s，可以達到預期目標(Ⅱ類水)；3—10月(6月份除外)，日平均引水流量為10～15 m3/s，可以達到預期目標(Ⅱ類水)；6月，日平均引水流量為15 m3/s(理論推算建議20 m3/s)，可以達到預期目標(Ⅱ類水)。

(2)引水啟停調(diào)度。根據(jù)實際水質(zhì)情況，進行調(diào)度優(yōu)化，適當減少引水天數(shù)，每年11月至次年2月，按照10 m3/s引7 d停3 d間歇引水，可停引9 d左右；3月、7月、8月和10月，按照15 m3/s引5 d停2 d間歇引水，可停引8 d左右；4月、5月、6月和9月，除了特殊條件下停止引水外，均保持每天引水。

(3)根據(jù)上述建議，在保持持續(xù)引水的情況下，年總引水量約為2.81 億m3，引水天數(shù)為332 d；考慮間歇性引水的情況下，年總引水量為3.26億m3，引水天數(shù)為280 d。

逐月的引水天數(shù)(汛期5—9月引水天數(shù)的依據(jù)是2018年實際引水天數(shù))和具體引水規(guī)模見表1。

表1 曹娥江引水工程2021年度引水方案表

5.2 降雨的應急引水方案

降雨徑流攜帶的地表污染物質(zhì)進入水體是導致水質(zhì)污染和指標數(shù)值提高的重要因素之一，因此需要在產(chǎn)生明顯的降雨徑流事件之后，在滿足區(qū)域防汛、防臺任務的前提下，加強應急引水量，以改善河道水質(zhì)狀況。

5.3 枯水期引水方案

枯水期為了維持水體的生態(tài)流量[9]，需要保證5 m3/s以上的生態(tài)基流作為基本的引水流量，并且每個月保證引水不少于20 d。

圖9以南渡橋為代表站，為5 m3/s持續(xù)引水、以10 m3/s引7 d，停7 d間歇引水、以15 m3/s引5 d，停10 d間歇引水的氨氮濃度的計算數(shù)值與引水流量的變化過程示意圖。從圖9可以看出5 m3/s持續(xù)引水情形下，氨氮濃度數(shù)值基本保持穩(wěn)定，10 m3/s、15 m3/s間歇引水一旦停止引水后，約3～5 d氨氮濃度數(shù)值達到5 m3/s持續(xù)引水的濃度數(shù)值。

圖9 不同引水方式下氨氮濃度計算數(shù)值變化示意圖

對間歇引水的調(diào)度方式進行了改進后，從圖10可知，經(jīng)計算得出：以10 m3/s流量引7 d，停3 d間歇引水，或以15 m3/s流量引5 d，停3 d間歇引水時,可使得氨氮濃度數(shù)值低于目標濃度數(shù)值。

圖10 不同引水方式下氨氮濃度計算數(shù)值變化示意圖

6 結(jié) 論

本次研究構建紹興平原河網(wǎng)區(qū)一維水動力水質(zhì)模型，探討不同引水流量對區(qū)域水環(huán)境的改善效果，并從決策目標水質(zhì)、引水流量等多個層面綜合考慮，構建以環(huán)境效益與經(jīng)濟效益相結(jié)合的評價方式，對不同引水方案進行評估優(yōu)選，確定出不同目標組合下的最優(yōu)化調(diào)度方案。

(1)較小引水流量下對研究區(qū)域河網(wǎng)水環(huán)境改善效果十分有限。引水流量較大時，能夠有效改善區(qū)域水環(huán)境狀況，適當?shù)牧髁考扔欣诟纳颇繕怂h(huán)境，又能將費用降到最低。

(2)在平原河網(wǎng)引水方案制定過程中，引水水量和調(diào)度模式都會對河網(wǎng)地區(qū)水質(zhì)產(chǎn)生差異性影響，采用數(shù)值模擬方法進行前置規(guī)劃十分必要。在模擬時還應結(jié)合控源截污工程實施進度，將上游污染源分布精準納入到模型邊界條件中。

(3)本研究建立引水多目標與評價方式，能夠從多個層面對不同引調(diào)水方案進行評估與優(yōu)選，得到綜合考慮環(huán)境效益與經(jīng)濟效益的最優(yōu)方案，為平原河網(wǎng)地區(qū)水環(huán)境長效管理與科學決策提供理論參考。