張澤宇,梁存峰
(1.太原理工大學水利科學與工程學院,山西太原 030024;2.晉中水文水資源勘測分局, 山西晉中 030600)
汾河一直是山西省經(jīng)濟發(fā)展的重要資源保障之一。隨著經(jīng)濟的發(fā)展,迅速增長的工農(nóng)業(yè)生產(chǎn)和城市生活取水量使得大部分汾河水量被引用消耗,河道實際流量迅速減少。據(jù)統(tǒng)計,汾河流域1980—2000 年平均地表水利用消耗率為64.8%,同期平均匯入黃河的水量為58 471 萬m3。根據(jù)水文站歷年觀測資料分析,汾河水庫以上各監(jiān)測斷面未有河干和斷流現(xiàn)象出現(xiàn),汾河水庫以下各站除趙城、河津兩站外均不同程度出現(xiàn)河干和斷流現(xiàn)象,其中斷流最嚴重的汾河二壩站平均斷流138 d/a,義棠站平均斷流100 d/a[1]。大量未經(jīng)處理達標的污水直接排入河道,使汾河干流水質遠低于山西省水功能區(qū)劃規(guī)定的水質目標。2008年山西省政府提出再現(xiàn)“汾河流水嘩啦啦”和汾河水生態(tài)明顯改善的目標,在汾河上中游和下游分別建設了萬家寨引黃入晉工程南干線和臨汾馬房溝引沁入汾工程兩個大型調水工程。引黃入晉工程南干線向汾河引水3.2 億m3/a,臨汾馬房溝引沁入汾工程設計引水量5 902 萬m3,除滿足生產(chǎn)、生活用水外,用于維持汾河生態(tài)需水,保證汾河不斷流。污染治理是清水復流工程的一個重要方面。筆者重點針對功能區(qū)水質目標,根據(jù)1956—2006 年系列不同水文條件分析計算汾河復流前后的納污能力,旨在為汾河清水復流后水環(huán)境的綜合治理提供參考依據(jù)。
根據(jù)《山西省水功能區(qū)劃》, 汾河干流總長693.8 km,共劃分5 個一級水功能區(qū):靜樂源頭保護區(qū)、靜樂婁煩開發(fā)利用區(qū)、古交保留區(qū)、太原運城開發(fā)利用區(qū)和河津緩沖區(qū)。在開發(fā)利用區(qū)內又劃分了20 個二級水功能區(qū),具體見表1。
表1 汾河復流前水功能區(qū)納污能力
汾河干流有入河排污口110 處(包括入汾支流口15 處),其中入河廢污水量小于300 萬m3/a 的排污口81 處,300 萬~500 萬m3/a 的排污口8 處,大于500 萬m3/a 的排污口21 處,廢污水入河量達33 654萬m3/a,COD 入河量達49 843 t/a,NH3-N 入河量達9 015 t/a。在不同水功能區(qū)上,污廢水排放主要集中在排污控制區(qū), 其污廢水年排放量達17 413萬m3/a,COD 入河量達31 326 t/a,NH3-N 入河量達6 858 t/a。
根據(jù)2006 年汾河水功能區(qū)水質監(jiān)測資料,汾河干流23 個功能區(qū)中,水質達到Ⅱ類的功能區(qū)為1 個;達到Ⅲ類的為1 個;達到Ⅳ類的2 個;達到劣Ⅴ類的19個,占汾河水功能區(qū)總數(shù)的82.6%。對照水功能區(qū)水質目標,僅3 個功能區(qū)水質達標,達標率僅為13%。
納污能力計算的重點區(qū)是水功能區(qū)的開發(fā)利用區(qū)(飲用水源區(qū)除外),對于需要改善水質的保留區(qū)和緩沖區(qū),也同時計算其納污能力。根據(jù)各二級水功能區(qū)的設計條件和水質目標,以各水功能區(qū)為計算單元,自上而下逐段進行水域納污能力計算。
根據(jù)污染物在各河段橫斷面上均勻混合的特點,采用河流一維模型[2-4]:
式中:W 為水功能區(qū)納污能力,t/a;Q 為水功能區(qū)設計流量,m3/s;q 為水功能區(qū)入河污水量,m3/s;ρs 為水功能區(qū)水質目標,mg/L;ρ0為水功能區(qū)上斷面污染物質量濃度,mg/L;k 為污染物綜合降解系數(shù),d-1;x 為水功能區(qū)上斷面到下斷面的距離,km;x1為簡化后排污口到下計算斷面的距離,km;u 為水功能區(qū)設計流量下的河段平均流速,m/s。
水功能區(qū)水質目標按照《山西省水功能區(qū)劃》規(guī)定的水質目標對照相應標準予以確定。以COD、NH3-N 為主要污染控制指標。
①設計流量。汾河復流前,有水文站的水功能區(qū),設計流量采用1956—2006 年系列實測資料進行分析計算。無水文資料的河段,設計流量按以下原則確定:距水文站較近,區(qū)間無較大支流加入或大的取水口,直接借用鄰近水文站的資料推求設計流量;距水文站較遠,區(qū)間有較大支流加入或大的取水口,通過水量平衡計算,確定設計流量。汾河復流后,各控制站設計流量采用“汾河清水復流工程水量調度方案”中的相應資料確定。 ②設計流速。有水文資料的根據(jù)水文站斷面水位流速關系確定;無水文資料的,采用曼寧公式[5]計算。
綜合降解系數(shù)根據(jù)汾河實測資料確定,COD 的k 值取值范圍在0.20 ~0.48d-1,平均為0.34d-1;NH3-N 的k 值取值范圍在0.12 ~0.32d-1,平均為0.22d-1[6]。
根據(jù)SL 348—2006《水域納污能力計算規(guī)程》,采用90%保證率最枯月平均流量對汾河復流前即汾河未實施引黃、引沁輸水前的納污能力進行分析計算。經(jīng)計算,復流前COD 和NH3-N 的納污能力分別為21 687 t/a 和975 t/a,其中開發(fā)利用區(qū)納污能力分別為20 875 t/a 和941 t/a,占各功能區(qū)納污能力總和的96.3%和96%(表1)。
為反映汾河在不同水量條件下的納污能力,結合北方季節(jié)性河流的特點,選取不同保證率(20%、50%、75%及95%)枯季(10 月至翌年5 月)平均流量進行納污能力的分析計算(表2)。結果顯示:在不同保證率來水條件下,COD 的平均納污能力由95%保證率來水條件下的29 495 t/a,增大到20%保證率來水條件下的85 375 t/a,增加了1.89 倍;NH3-N 平均納污能力由95%保證率來水條件下的1 343 t/a,增大到20%保證率來水條件下的3 054 t/a,增加了1.27 倍。
偏枯水年、枯水年枯季(75%、95%保證率來水條件)全河COD、NH3-N 總納污能力小于現(xiàn)狀污染物入河量,也即在偏枯水年、枯水年現(xiàn)狀污染物入河量大于汾河本身的水環(huán)境承載能力。平水年(保證率50%來水條件)全河COD 納污能力總體上大于現(xiàn)狀污染物入河量,但太原、介休、霍州、稷山等區(qū)縣污染物現(xiàn)狀入河量仍大于其納污能力,表明上述城市附近河段污染物入河量已超過其水環(huán)境承載能力。偏豐水年(20%保證率來水條件)全河總納污能力遠大于現(xiàn)狀污染物入河量,但就局部河段如太原市區(qū)、介休等區(qū)縣所在河段而言,其納污能力仍小于現(xiàn)狀入河量,尤其是太原市現(xiàn)狀入河量為其納污能力的1.6 倍。
汾河復流后各控制站設計流量采用山西省水利廳編制的《汾河清水復流工程水量調度方案》中相應資料。復流后各控制斷面設計流量平均較復流前增大0.5 m3/s 左右。汾河入黃口設計流量由復流前的0.12 m3/s 增加為復流后的1.34 m3/s。經(jīng)計算,復流后COD 和NH3-N 的納污能力分別為30 275 t/a和1 283 t/a(表3)。
表2 復流前不同水量條件下水功能區(qū)納污能力 t/a
表3 汾河復流后水功能區(qū)納污能力 t/a
考慮汾河中下游水質已遭受嚴重污染的實際,本次限制排污總量的原則是流域內各行政區(qū)排放的主要污染物總量不再增加,具體計算方法為:若功能區(qū)現(xiàn)狀入河量小于或等于該功能區(qū)允許納污量,則將該目標控制總量視為功能區(qū)入河污染物控制總量;反之,將該功能區(qū)允許納污量作為其入河污染物控制總量。分析計算得汾河干流復流后COD、NH3-N 限制排污總量分別為25 139 t/a 和1 158 t/a。
以各功能區(qū)納污能力以及現(xiàn)狀污染排放狀況為基礎,統(tǒng)計得出汾河干流復流后水力條件下,COD、NH3-N 應削減量分別為24 704 t/a 和7 857 t/a(表4),平均削減率分別為49.6%和87.2%?,F(xiàn)狀排放各行政分區(qū)主要污染物削減率見圖1。
表4 汾河復流后現(xiàn)狀排放行政分區(qū)污染物削減量 t/a
圖1 現(xiàn)狀排放各行政分區(qū)主要污染物削減率
達標排放削減率是指假設現(xiàn)有排污口都達到污水綜合Ⅱ類排放標準,按照水功能區(qū)納污能力和限制排污總量的要求計算的污染物削減率[7]。以汾河干流主要行政區(qū)為單元,達標排放后COD、NH3-N 入河量分別為38 247 t/a 和5510 t/a,仍不能滿足各水功能區(qū)的水質目標, 仍需削減15 936t/a和4 356 t/a。達標排放后COD、NH3-N 平均削減率41.7%和79.0%,各行政分區(qū)主要污染物削減率見圖2。
圖2 達標排放后各行政分區(qū)主要污染物削減率
經(jīng)計算,汾河復流前COD 和NH3-N 的納污能力分別為21 687 t/a 和975 t/a,其中開發(fā)利用區(qū)納污能力分別為20 875 t/a 和941 t/a,占各功能區(qū)納污能力總和的96.3%和96.6%。
復流后各功能區(qū)COD 和NH3-N 的納污能力分別為30 275 t/a 和1 283 t/a,總納污能力分別較復流前增加39.6%和31.6%。限制排污總量分別為25 139 t/a 和1 158 t/a,污染物削減量分別為24 704 t/a和7 857 t/a,平均削減率分別為49.6%和87.2%。若各入河排污口達標排放后,COD、NH3-N 仍要削減15 936 t/a 和4 356 t/a。
可見要達到“汾河清水復流”的水質目標要求,污染源達標排放是遠遠不夠的,必須依據(jù)限制排污總量進行控制:①要以水功能區(qū)管理為載體,依法向環(huán)保部門提出限制排污的意見;②要加強重要斷面的水質監(jiān)測以及入河排污總量的監(jiān)控,將確定的水功能區(qū)納污能力作為污染物減排依據(jù)和河道管理依據(jù);③強化入河排污口審批和監(jiān)督檢查,嚴格污染物總量控制和污染物濃度控制,確保水質目標的最終實現(xiàn)[8]。
[1] 張曉斌, 張建國,白繼中.汾河水環(huán)境容量計算及特征分析[J] .水資源與水工程學報,2010(2):135-138.
[2] 韓龍喜, 朱黨生,蔣莉華.中小型河流納污能力計算方法研究[J] .河海大學學報:自然科學版,2002(1):35-38.
[3] 張永良, 劉陪哲.水環(huán)境容量綜合手冊[M] .北京:清華大學出版社,1991.
[4] 宋宏杰, 馬軍霞,左其亭.鄭州市水環(huán)境承載能力計算及調控對策[J] .鄭州大學學報:工學版,2005(1):103-107.
[5] 盧偉.浙江省主要飲用水源地水環(huán)境承載能力與對策[J] .水文,2003(5):38-41.
[6] 陳永燦, 劉昭偉.三峽水庫水環(huán)境承載能力的評價和分析[J] .水科學進展,2005(5):715-719.
[7] 王秀蘭, 李紅亮.河北省水環(huán)境承載能力及污染物總量控制方案研究[J] .海河水利,2004(4):31-33.
[8] 汪恕誠.水環(huán)境承載能力分析與調控[J] .水利發(fā)展研究,2002(1):2-6.