董彭旭 杜慧玲 于曉英

摘要? ? 針對(duì)松花江流域水環(huán)境容量季節(jié)性變化情況，在動(dòng)態(tài)的水文設(shè)計(jì)條件下，基于二維岸邊排放模型，采用段尾控制法對(duì)松花江佳木斯段COD和NH3-N的實(shí)際水環(huán)境容量進(jìn)行了數(shù)值計(jì)算與分析。結(jié)果表明，松花江佳木斯段COD與NH3-N實(shí)際水環(huán)境容量變化趨勢(shì)相似，NH3-N波動(dòng)幅度大于COD的波動(dòng)幅度;實(shí)際水環(huán)境容量季節(jié)性變化較大，豐水期最大，枯水期最小，應(yīng)實(shí)施動(dòng)態(tài)的管理。研究結(jié)果可為松花江水環(huán)境管理和保護(hù)提供基礎(chǔ)理論依據(jù)。

關(guān)鍵詞? ? 動(dòng)態(tài)水環(huán)境容量;COD;NH3-N;松花江佳木斯段

中圖分類號(hào)? ? X832? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼? ? A? ? ? ? 文章編號(hào)? ?1007-5739（2019）11-0166-03

Abstract? ? In view of the seasonal variations of water environmental capacity in Songhua River，under dynamic hydrological design conditions，based on two-dimensional side-discharge model，the COD and NH3-N actual water environmental capacity of Jiamusi section in Songhua River was calculated and analyzed by using section-ending control method.The results showed that the variation trend of COD and NH3-N actual water environmental capacity was similar in Songhua River Jiamusi section，and the fluctuation range of NH3-N was greater than that of COD.Seasonal difference of the capacity was obvious，the maximum and minimum value of the COD and NH3-N actual water environmental capacity appeared in flood season and dry season，respectively.Dynamic management should be implemented in Songhua River.This research can provide theoretical accordance for water environmental management and protection of Songhua River.

Key words? ? dynamic water environmental capacity;chemical oxygen demand（COD）;ammonia nitrogen（NH3-N）;Jiamusi section in Songhua River

水環(huán)境容量是指在一定的環(huán)境目標(biāo)條件下，水體能夠容納污染物的最大量或自身調(diào)節(jié)凈化并保持生態(tài)平衡的能力，反映了水體納污能力的大小[1]。隨著我國水環(huán)境管理中對(duì)以水質(zhì)改善為目標(biāo)的容量總量控制認(rèn)識(shí)的不斷加深，科學(xué)合理地確定河湖水庫水環(huán)境容量的技術(shù)方法成為諸多學(xué)者們不斷研究的方向[2]。影響河流水環(huán)境容量的因素很多，主要包括水域水質(zhì)環(huán)境功能區(qū)要求、河道的水文條件、污染物的特性等。我國對(duì)于水環(huán)境容量的研究已進(jìn)行了40多年[3-4]，在理論及總量控制方法2個(gè)方面取得了豐碩成果，并提出了可分配水環(huán)境容量的概念及辦法[5]。水環(huán)境容量已在我國許多重要水體中開展了研究和應(yīng)用，為水污染治理和保護(hù)提供了科學(xué)依據(jù)，豐富了水環(huán)境容量的理論和研究方法。

松花江佳木斯江段位于松花江下游，在同江附近匯入中俄界河黑龍江，該江段的水質(zhì)直接影響松花江的出境水質(zhì)[6]。因此，對(duì)其水環(huán)境容量進(jìn)行研究，可以為流域環(huán)境規(guī)劃、污染防治與環(huán)境管理提供重要依據(jù)，對(duì)加強(qiáng)松花江流域水資源保護(hù)、防治松花江水質(zhì)污染具有十分重要的意義。松花江流域位于我國北方高緯度寒冷地區(qū)，存在冰封期與非冰封期的差異，河流徑流量的季節(jié)性變化較大且冬季水溫低，因而松花江水環(huán)境容量存在季節(jié)性差異變化，需進(jìn)行動(dòng)態(tài)水環(huán)境容量分析并實(shí)施動(dòng)態(tài)管理。近年來，國內(nèi)外關(guān)于動(dòng)態(tài)水環(huán)境容量的研究較多，但對(duì)松花江流域動(dòng)態(tài)水環(huán)境容量的研究相對(duì)較少且大部分研究集中在2013年以前[7]。

1? ? 研究方法

1.1? ? 環(huán)境容量計(jì)算模型的選取

松花江佳木斯江段多年平均流量約2 149 m3/s，河槽寬度約500～2 000 m，平均水深約4 m，平均流速約0.4 m/s，屬于寬淺大型河流，其主要支流從上游至下游主要包括湯旺河、梧桐河和安邦河。

根據(jù)《水域納污能力計(jì)算規(guī)程》（GB/T 25173—2010），河流水環(huán)境容量計(jì)算模型由簡(jiǎn)單到復(fù)雜依次為河流零維模型、河流一維模型和河流二維模型。其中，河流二維模型適用于大型河段（流量Q≥150 m3/s），污染物在河段非均勻混合，故本文選取河流二維模型。研究江段沿岸的污染源為岸邊排放，因而采用二維岸邊排放模型對(duì)污染物在水體中的遷移擴(kuò)散規(guī)律進(jìn)行模擬。敏感點(diǎn)位置取岸邊，并采用段尾控制法，推導(dǎo)出研究江段的水環(huán)境容量計(jì)算公式，見式（1）。

式（1）中，cs—污染物在研究河段的水質(zhì)標(biāo)準(zhǔn)（mg/L）;c0—污染物在研究河段的本底值（mg/L）;xi—敏感點(diǎn)到排污口的距離（m）;My—河流縱向、橫向的彌散系數(shù)（m2/s）;ux—河流的平均流速（m/s）;K1—污染物的衰減速度常數(shù)（d-1）;B—河流水面寬度（m）;H—河流平均水深（m）。

1.2? ? 研究江段計(jì)算單元的劃分

根據(jù)松花江佳木斯江段的國家水質(zhì)考核斷面及其功能區(qū)目標(biāo)，將研究江段劃分為3個(gè)計(jì)算單元，分別為佳木斯上—佳木斯下、佳木斯下—江南屯、江南屯—同江，計(jì)算單元屬性見表1。松花江佳木斯江段監(jiān)測(cè)斷面位置見圖1。

1.3? ? 水文參數(shù)的選取

水文參數(shù)根據(jù)2009—2013年松花江佳木斯站的水文監(jiān)測(cè)數(shù)據(jù)，取每個(gè)月的最小實(shí)測(cè)流量作為當(dāng)月的設(shè)計(jì)流量，流速、平均水深、水面寬度等水文參數(shù)取相應(yīng)月份的最小流量對(duì)應(yīng)的值。由此計(jì)算出的水環(huán)境容量隨月份變化呈現(xiàn)動(dòng)態(tài)變化。

1.4? ? 橫向彌散系數(shù)My及污染物綜合衰減系數(shù)

橫向彌散系數(shù)My采用Taylor公式進(jìn)行求解，見式（2）。

式（2）中，My—河流橫向彌散系數(shù)（m2/s）;B—河流水面平均寬度（m）;h—河流平均水深（m）;g—重力加速度，取9.8 m/s2;i—平均縱向坡降，松花江佳木斯段取0.1‰。

松花江佳木斯段的橫向彌散系數(shù)My見表2。污染物綜合衰減系數(shù)K1依據(jù)已有文獻(xiàn)以及經(jīng)驗(yàn)系數(shù)法初步確定，再根據(jù)已有斷面COD和NH3-N的監(jiān)測(cè)數(shù)據(jù)進(jìn)行擬合分析和修正，經(jīng)過率定后的污染物綜合衰減系數(shù)K1值見表2。

1.5? ? 主要支流與排污口概化

松花江佳木斯段主要支流從上游至下游依次為湯旺河、格金河、阿凌達(dá)河、梧桐河、安邦河5個(gè)支流，其各計(jì)算單元內(nèi)主要排污口排放量及其與終點(diǎn)斷面的距離見表3。該江段主要支流、排污口與斷面位置關(guān)系的概化示意圖見圖2。

2? ? 結(jié)果與分析

結(jié)合松花江佳木斯段的實(shí)際情況，將各個(gè)計(jì)算單元視為獨(dú)立的河段，由式（1）計(jì)算佳木斯上—佳木斯下、佳木斯下—江南屯及江南屯—同江3個(gè)計(jì)算單元的COD和NH3-N的實(shí)際水環(huán)境容量。實(shí)際水環(huán)境容量需考慮上游來水水質(zhì)的影響，根據(jù)實(shí)際水質(zhì)情況以及支流、排污口對(duì)水環(huán)境容量的影響進(jìn)行綜合計(jì)算。根據(jù)松花江佳木斯段2016年各監(jiān)測(cè)斷面的實(shí)際水質(zhì)監(jiān)測(cè)結(jié)果，佳木斯下和江南屯2個(gè)斷面的COD和NH3-N的監(jiān)測(cè)值與其功能區(qū)劃要求不一致，均為Ⅲ類水，優(yōu)于水環(huán)境功能區(qū)要求，對(duì)COD和NH3-N的實(shí)際水環(huán)境容量產(chǎn)生影響。另外，支流和排污口處的匯水將改變水體的流量、流速、水深、水面寬度等要素從而影響水環(huán)境容量。本文受水文資料的限制，假設(shè)支流和排污口匯水并不會(huì)改變水體的寬度和深度，只對(duì)流速產(chǎn)生影響。根據(jù)流域水文資料，佳木斯站位于佳木斯市西區(qū)污水處理廠排污口與泉林生態(tài)農(nóng)業(yè)有限公司排污口之間，將各支流及排污口匯入的位置作為流速變化的分界斷面，重新核算不同斷面沿江流速變化，并據(jù)此計(jì)算實(shí)際COD和NH3-N水環(huán)境容量。松花江佳木斯段考慮上游來水水質(zhì)和支流、排污口匯入影響的COD實(shí)際水環(huán)境容量見表4和圖3，NH3-N實(shí)際水環(huán)境容量見表5和圖4。

根據(jù)圖3和圖4可知，松花江佳木斯段COD和NH3-N的實(shí)際水環(huán)境容量呈現(xiàn)出動(dòng)態(tài)變化趨勢(shì)，變化趨勢(shì)相似，在夏季7—9月明顯高于其他月份，最大值均出現(xiàn)在8月。佳木斯上—佳木斯下段COD實(shí)際水環(huán)境容量最大值與最小值的比值為3.38，NH3-N實(shí)際水環(huán)境容量最大值與最小值的比值為4.14;佳木斯下—江南屯段COD實(shí)際水環(huán)境容量最大值與最小值的比值為4.11，NH3-N實(shí)際水環(huán)境容量最大值與最小值的比值為4.49;江南屯—同江段COD實(shí)際水環(huán)境容量最大值與最小值的比值為3.36，NH3-N實(shí)際水環(huán)境容量最大值與最小值的比值為4.17。由此說明，松花江佳木斯段的實(shí)際環(huán)境容量NH3-N波動(dòng)幅度大于COD的波動(dòng)幅度。

松花江佳木斯段各水期COD和NH3-N的實(shí)際水環(huán)境容量平均值分別見圖5和圖6?？芍苫ń涯舅苟蜟OD和NH3-N的實(shí)際水環(huán)境容量在豐水期最大，其次為平水期，枯水期最小。這受季節(jié)影響，不同水期水體溫度不同，導(dǎo)致污染物綜合衰減系數(shù)不同：豐水期的水溫高，生化反應(yīng)速率快，衰減系數(shù)大;而枯水期正值冰封期，水溫低，衰減系數(shù)小。水體的實(shí)際水環(huán)境容量與污染物綜合衰減系數(shù)成正相關(guān)，受溫度影響。另外，根據(jù)2009—2013年松花江佳木斯站的水文監(jiān)測(cè)數(shù)據(jù)，松花江佳木斯段3個(gè)水期的平均流速呈現(xiàn)豐水期（1.07 m/s）>平水期（0.71 m/s）>枯水期（0.27 m/s）的特征，水環(huán)境容量與水體流速隨水期的變化趨勢(shì)相似。因此，松花江佳木斯段各水期的實(shí)際水環(huán)境容量與水體流速成正相關(guān)。

3? ? 結(jié)論

（1）松花江佳木斯段COD與NH3-N實(shí)際水環(huán)境容量的變化趨勢(shì)相似，夏季明顯高于其他月份，最大值均出現(xiàn)在?8月，NH3-N波動(dòng)幅度大于COD的波動(dòng)幅度。

（2）松花江佳木斯段COD和NH3-N的實(shí)際水環(huán)境容量在豐水期最大，枯水期最小。除受溫度影響外，也與水體流速有關(guān)，實(shí)際水環(huán)境容量與水體流速成正相關(guān)。實(shí)際水環(huán)境容量的季節(jié)性變化較大，應(yīng)進(jìn)行動(dòng)態(tài)管理，以保障水質(zhì)安全。

4? ? 參考文獻(xiàn)

[1] 程功.大慶市安肇新河流域水環(huán)境容量核算[J].黑龍江水利科技，2017，45（8）：75-77.

[2] 陶亞，陳宇軒，趙喜亮，等.基于EFDC模型的阿什河水環(huán)境容量季節(jié)性分析[J].環(huán)境工程，2017，35（7）：65-69.

[3] 馮民權(quán)，胡芳，吳波，等.漢江安康段排污口近區(qū)污染擴(kuò)散規(guī)律[J].西北大學(xué)學(xué)報(bào)（自然科學(xué)版），2009，39（5）：887-890.

[4] 梁兢波，張國成.漢江安康段水質(zhì)衛(wèi)生評(píng)價(jià)[J].環(huán)境與健康雜志，2001，18（5）：283-285.

[5] 張濤，任昭，王鳳金，等.漢江漢中段水質(zhì)調(diào)查與分析[J].監(jiān)測(cè)分析，2011（2）：86-89.

[6] 李永亮，牟學(xué)軍，李健，等.松花江佳木斯江段休養(yǎng)生息前后水環(huán)境質(zhì)量變化研究[J].中國水利，2014（11）：50-52.

[7] 杜慧玲，于曉英，曲茉莉.松花江干流哈爾濱段COD和氨氮?jiǎng)討B(tài)水環(huán)境容量研究[J].水資源與水工程學(xué)報(bào)，2018，29（2）：69-75.

基金項(xiàng)目? ?國家水體污染控制與治理科技重大專項(xiàng)（2014ZX07201012）。

作者簡(jiǎn)介? ?董彭旭（1988-），女，黑龍江哈爾濱人，工程師，碩士，從事水環(huán)境保護(hù)研究工作。

收稿日期? ?2019-03-04