徐 飛， 周 飛， 沈 樂

(江蘇省水文水資源勘測局南京分局，江蘇 南京 210008)

平原感潮河網(wǎng)地區(qū)因其特殊的地理位置,經(jīng)濟發(fā)達,污染物排放量較大;但與山丘河道相較,又存在水動力條件不足,河道自凈能力較弱,水流往復(fù),污染物易于擴散等特點。因此,針對其水動力和水環(huán)境特征,進行合理的水環(huán)境容量計算,進而為污染物總量控制提供科學(xué)的依據(jù),對緩解當?shù)亟?jīng)濟發(fā)展與水環(huán)境保護之間的矛盾意義重大。張家港市是典型的平原感潮河網(wǎng)地區(qū),經(jīng)濟發(fā)達、人口密集、河道密布,高頻的人類生產(chǎn)生活也帶來了較為嚴重的水環(huán)境污染問題。且由于河網(wǎng)縱橫、河床比降小、水流往復(fù),水體往往自凈能力較差。因此,對張家港市進行水環(huán)境容量計算研究十分必要,可為污染物總量控制提供科學(xué)依據(jù)。

1 研究方法

1.1 構(gòu)建水環(huán)境數(shù)學(xué)模型

基于當?shù)睾拥罃嗝?、水位、流量等水文資料及污染物調(diào)查情況,概化河網(wǎng),建立研究區(qū)域水量水質(zhì)模型。本次水動力模型采用圣維南方程組(式1)進行計算,水質(zhì)模型采用污染物對流擴散方程(式2、式3)進行計算。

(1)

其中,Q為河流流量,m3/s;x為沿水流方向的空間坐標,m;BW為調(diào)蓄寬度,m;Z為河道水位,m;t為時間坐標,s;q為支流入流流量(流入為正,流出為負),m3/s;u為斷面平均流速,m/s;g為重力加速度,m/s2;A為主槽過水斷面面積,m2;B為主流斷面寬度,m;n為河道糙率;R為水力半徑,m。

(2)

(3)

其中,Q為河道流量,m3/s;Z為河道水位,m;A為河道面積,m2;Ex為縱向分散系數(shù);C為水流輸送的物質(zhì)濃度,mg/L;Ω為河道叉點節(jié)點的水面面積,m2;j為節(jié)點編號;I為與節(jié)點j相關(guān)聯(lián)的河道編號;Sc為與輸送物質(zhì)濃度相關(guān)的衰減項;Kd為衰減因子;S為外部的源匯項。

1.2 參數(shù)率定與驗證

通過野外調(diào)水試驗,利用先進的在線監(jiān)測設(shè)備進行水量水質(zhì)同步監(jiān)測,從而獲得一系列的水文、水質(zhì)實測數(shù)據(jù),在此基礎(chǔ)上對張家港市河道的糙率、水質(zhì)綜合降解系數(shù)等參數(shù)進行率定和驗證。

1.3 水環(huán)境容量計算

本文采用總體達標法[12]對張家港市水環(huán)境容量進行計算。選取90%保證率枯水年最枯月平均流量作為設(shè)計水文條件,根據(jù)研究區(qū)水環(huán)境功能區(qū)劃水質(zhì)達標要求,計算概化主干河網(wǎng)的水環(huán)境容量,水環(huán)境容量計算公式如式4、式5所示:

(4)

Wij環(huán)境容量=Q0ij(Csij-Coij)+KVijCsij

(5)

其中,αij為不均勻系數(shù);Q為進口斷面入流流量,m3/s;Csij為該水體的目標水質(zhì),mg/L;Coij為進口斷面水質(zhì)濃度,mg/L;K為水質(zhì)降解系數(shù),d-1;V為水體體積,m3。

由于研究區(qū)河網(wǎng)受潮汐影響,水流往復(fù)方向不固定,故取設(shè)計水文條件下的最枯月正向流、逆向流和流量為0時候的平均值作為水環(huán)境容量,詳細計算如公式6所示:

W=(W1×d1+W2×d2+W3×d3)/d

(6)

其中,W為往復(fù)流河道環(huán)境容量,t/d;d為最枯月天數(shù),d;W1、W2、W3為分別為正向流、逆向流和流量為0時候的水環(huán)境容量,t/d;d1、d2、d3為分別為正向流、逆向流和流量為0時候的天數(shù),d。

2 水環(huán)境容量計算

2.1 模型構(gòu)建

本文采用節(jié)點、河網(wǎng)、斷面三種參數(shù)來進行河網(wǎng)概化,為提高河網(wǎng)計算的精度,本文在確保河道長度、坡降、流量、水位、水利工程基本與實際相符的前提下,將張家港復(fù)雜水系概化成367條河道。

2.2 模型率定和驗證

本文設(shè)定模型外部邊界為朝東圩港、一干河等,采用野外同步水量水質(zhì)監(jiān)測值進行計算,長江邊界來水按Ⅲ類水標準值輸入模型。通過模型率定,部分代表監(jiān)測點位水量、水位實測值和模擬值結(jié)果對比詳見圖1、圖2。

圖1 水動力模型水量實測值與模擬值對比圖

圖2 水動力模型水位實測值與模擬值對比圖

由圖1、圖2可知，代表監(jiān)測點位流量和水位模型計算值與實測值均較為吻合，由圖1、圖2可知,各監(jiān)測點位流量和水位模型計算值與實測值均較為吻合,其中流量的相對誤差在30%以內(nèi),故本文判定該水動力模型可基本模擬研究區(qū)的水動力過程,且根據(jù)率定結(jié)果得出現(xiàn)狀河道糙率取值范圍為0.01～0.04。

本文采用原型引調(diào)水實驗的同步監(jiān)測值對水質(zhì)模型進行率定。為調(diào)高計算的精確性和穩(wěn)定度,時間步長t設(shè)置為15s。根據(jù)實地污染源調(diào)查,將全市范圍內(nèi)共計127家企業(yè)、22家生活及工業(yè)污水處理廠排污口概化成102個集中排污口,將約514家畜禽養(yǎng)殖場排污口概化成173個排污口,加入水質(zhì)模型中,調(diào)整各參數(shù),使代表監(jiān)測點的模擬值與實測值趨于吻合。最終通過水質(zhì)模型率定,COD、氨氮的模擬值和實測值對比詳見圖3、圖4。

圖3 COD濃度實測值與模擬值對比圖

圖4 COD濃度實測值與模擬值對比圖

根據(jù)圖3、圖4可知，率定后的水質(zhì)模型COD、氨氮模擬值與野外實測結(jié)果吻合度良好,相對誤差均小于40%,結(jié)果表明該水質(zhì)模型基本可模擬研究區(qū)污染物遷移過程。根據(jù)水質(zhì)模型率定結(jié)果得出COD和氨氮降解系數(shù)取值范圍分別為0.04～0.08/d和0.03～0.08/d,縱向擴散系數(shù)為2.5m2/s。

3.3 水環(huán)境容量計算

張家港市共有37個水環(huán)境功能區(qū),其中長江張家港飲用、工業(yè)用水區(qū)和一干河張家港飲用水源區(qū)水質(zhì)要求為Ⅱ類,不參與本次水環(huán)境容量計算,故本文只有35個功能區(qū)的水環(huán)境容量計入最終測算結(jié)果,水質(zhì)目標詳見表1。根據(jù)張家港市水污染調(diào)查結(jié)果和水環(huán)境管理目標,本文選取COD和氨氮作為特征污染物計算水環(huán)境容量。

表1 張家港市各水環(huán)境功能區(qū)環(huán)境容量

本文采用P-III型曲線對獲取的張家港市1967-2013年的多年月平均降水量長序列資料進行水文頻率分析,計算得到90%枯水年典型年為1971年,再依據(jù)太湖流域相關(guān)測站水位資料選取1971年最枯月2月的平均水位作為計算水位。由于一些河道受到潮汐影響易于形成往復(fù)流,故對于流向不定的河道采用流量為0時的設(shè)計水位對應(yīng)的水量作為設(shè)計流量。

最終,利用總體達標法計算得出張家港市各主干河道的水環(huán)境容量,進而得出各水環(huán)境功能區(qū)環(huán)境容量,詳細結(jié)果見表1。此外功能區(qū)外主要河道COD和氨氮的環(huán)境容量分別為9132.5t/a和564.3t/a,故張家港全市COD和氨氮水環(huán)境容量為分別為43742.9t/a和3527.6t/a。

4 結(jié) 論

平原河網(wǎng)地區(qū)由于污染物排放量較大且水體自凈能力不足,往往水污染嚴重,對其進行污染物總量控制十分必要。本文通過建立張家港市水量水質(zhì)模型,并利用野外高頻在線監(jiān)測數(shù)據(jù)對水環(huán)境數(shù)學(xué)模型進行了率定驗證,結(jié)果顯示模型擬合度較高,基本可用于描述研究區(qū)水量水質(zhì)過程,從而確定了現(xiàn)狀河道糙率和水質(zhì)綜合降解系數(shù),其中COD和氨氮降解系數(shù)取值范圍分別為0.04～0.08/d和0.03-0.08/d。在此基礎(chǔ)上,通過對長序列資料的水文頻率分析得出90%保證率典型年最枯月為1971年2月,并選取其作為計算水文條件,根據(jù)境容量,最終得出張家港市現(xiàn)狀水環(huán)境容量COD為43742.9t/a,氨氮為3527.6t/a。建議以此為約束條件,通過加快城鎮(zhèn)污水處理配套設(shè)施建設(shè)、促進工業(yè)企業(yè)優(yōu)化工藝達標排放、加強農(nóng)業(yè)面源污染整治力度等有效手段,一方面從源頭減少污染物排放,另一方面提高污水處理率,從而有效削減污染物排放總量,切實改善張家港市水環(huán)境質(zhì)量。