柯斌，勞國(guó)民

（浙江省水文局，浙江 杭州 310009）

1 納污能力計(jì)算模型

根據(jù)浙江省水功能區(qū)的實(shí)際情況，選用的計(jì)算模型包括：山區(qū)性河流一維均勻混合模型、湖泊和水庫(kù)零維模型、平原河網(wǎng)零維模型、感潮河段組合模型。

1.1 山區(qū)性河流—維均勻混合模型

該模型適用于寬深比不大的河流，在較短的時(shí)間內(nèi)，污染物能在河段橫斷面上均勻混合。

由于污染物一般是沿河岸分多處排放的，每一河段內(nèi)可能存在多個(gè)污染源，同時(shí)，規(guī)劃各排污口的位置具有不確定性，在此采用概化的方法，即認(rèn)為污染物排放口在同一功能區(qū)內(nèi)沿河均勻分布，此概化體現(xiàn)了污染源分布的一種平均狀況，混合過(guò)程也假定在排污口斷面瞬時(shí)完成均勻混合，圖1為某一功能區(qū)污染源均勻排放的示意圖。

圖1 功能區(qū)中污染源均勻排放示意圖

功能區(qū)終止斷面的污染物濃度按下式計(jì)算：

式中：Cs為終止斷面的污染物濃度（mg／L）；C0為起始斷面的污染物濃度（mg／L）；k為污染物綜合衰減系數(shù)（L／s）；L為功能區(qū)長(zhǎng)度（m）；u為設(shè)計(jì)流量下功能區(qū)的斷面平均流速（m／s）；W 為進(jìn)入功能區(qū)的污染物量（t／a）；Q 為功能區(qū)設(shè)計(jì)流量（m3／s）。

當(dāng)Cs剛好滿(mǎn)足該功能區(qū)水質(zhì)目標(biāo)時(shí)，進(jìn)入該功能區(qū)的污染物量即為其納污能力，計(jì)算公式如下：

1.2 湖泊和水庫(kù)零維模型

以圖2水庫(kù)為代表，設(shè)有n條支流匯入湖體，各匯入支流的流 量 及 污 染 物 濃 度 分 別 為q1，q2，q3，……，qn；C01，C02，……，C0n，水庫(kù)水質(zhì)目標(biāo)為Cs，在這種情況下，計(jì)算水庫(kù)的納污能力。

圖2 水庫(kù)污染物遷移示意圖

假定水量為穩(wěn)態(tài)，則出流流量q出＝Σqi，水庫(kù)納污能力計(jì)算公式為：

式中：qi為i條支流入流流量（m3／s）；C0i為i條支流入流濃度（mg／L）；B為水庫(kù)不均勻系數(shù)（無(wú)量綱）；V為設(shè)計(jì)水文條件下的水庫(kù)庫(kù)容（m3）；其余指標(biāo)同上。

由于水庫(kù)多位于河流上游，水質(zhì)較好，出流水質(zhì)與入流水質(zhì)接近，因此上式可簡(jiǎn)化為：

W＝31.536×BkVCs

1.3 平原河網(wǎng)零維模型

平原河網(wǎng)的功能區(qū)采用零維模型，根據(jù)不同水文條件下的河網(wǎng)容積，計(jì)算整個(gè)河網(wǎng)相應(yīng)的納污能力，再近似按河長(zhǎng)或河段蓄量分配到功能區(qū)，公式如下：

式中：W河網(wǎng)為整個(gè)河網(wǎng)的納污能力（t／a）；V河網(wǎng)為設(shè)計(jì)水文條件下整個(gè)河網(wǎng)槽蓄量（m3）；W功能區(qū)為功能區(qū)納污能力（t／a）；l為功能區(qū)河長(zhǎng)（m）；L河網(wǎng)為整個(gè)河網(wǎng)功能區(qū)總河長(zhǎng)（m）；v為功能區(qū)河段槽蓄量（m3）；其余符號(hào)意義同前。

1.4 感潮河段組合模型

感潮河段水量由上游下泄水量和下游潮水量組成 （見(jiàn)圖3）。

圖3 感潮河段水量組成示意圖

河段納污能力可由上游來(lái)水納污能力和下游潮水納污能力2部分組成，計(jì)算公式如下：

式中：b為感潮河段不均勻系數(shù)；V潮為潮水量（m3）；其余符號(hào)意義同前。

上游來(lái)水納污能力按照山區(qū)性河流一維均勻混合模型計(jì)算，設(shè)計(jì)流量和設(shè)計(jì)流速采用相應(yīng)保證率的數(shù)值；感潮河段納污能力采用零維模型計(jì)算，潮水量采用相應(yīng)保證率潮差進(jìn)行計(jì)算。

2 模型參數(shù)選擇

上述納污能力計(jì)算模型涉及的參數(shù)主要有綜合衰減系數(shù)k、湖庫(kù)不均勻系數(shù)B、感潮河段不均勻系數(shù)b。

2.1 綜合衰減系數(shù)k

綜合衰減系數(shù)的影響因素眾多，相當(dāng)復(fù)雜，一般采用下述方法進(jìn)行分析。

2.1.1 分析借用

對(duì)于以前在環(huán)評(píng)、環(huán)保規(guī)劃、環(huán)保科研等工作中可供利用的有關(guān)資料經(jīng)過(guò)分析檢驗(yàn)后采用。

無(wú)資料時(shí)，可借用水力特性、污染狀況、及地理、氣象條件相似的鄰近河流的資料。

2.1.2 實(shí)測(cè)法

選取一個(gè)河道順直、水流穩(wěn)定、中間無(wú)支流匯入、無(wú)排污口的河段，分別在河段上游 （A點(diǎn)）和下游 （B點(diǎn)）布設(shè)采樣點(diǎn)，監(jiān)測(cè)污染物濃度值，并同時(shí)測(cè)驗(yàn)水文參數(shù)以確定斷面平均流速，綜合衰減系數(shù)可由實(shí)測(cè)數(shù)據(jù)推算。

可以根據(jù)各種經(jīng)驗(yàn)公式推求綜合衰減系數(shù)。

本次綜合衰減系數(shù)主要采用分析借用和實(shí)測(cè)法，全省不同區(qū)域、不同水體類(lèi)型功能區(qū)的綜合衰減系數(shù)k取值范圍見(jiàn)表1。

表1 不同類(lèi)型水功能區(qū)的綜合衰減系數(shù)k的確定值表

續(xù)表1

2.2 湖庫(kù)不均勻系數(shù)B

湖泊、水庫(kù)的不均勻系數(shù)按表2進(jìn)行選用。

表2 湖泊、水庫(kù)不均勻系數(shù)表

2.3 感潮河段不均勻系數(shù)b

感潮河段的不均勻系數(shù)按表3進(jìn)行選用。

表3 感潮河段不均勻系數(shù)表

3 設(shè)計(jì)水文條件確定

根據(jù)SL 348—2006《水域納污能力計(jì)算規(guī)程》，納污能力計(jì)算應(yīng)采用90%保證率最枯月或近10a最枯月平均流量（水位、蓄量）作為設(shè)計(jì)水文條件。該條件偏嚴(yán)格，為反映不同設(shè)計(jì)水文條件下功能區(qū)的納污能力，從而為下一步污染物總量控制提供較為全面的科學(xué)數(shù)據(jù)，采用90%，75%和50%三種保證率最枯月平均流量 （水位、蓄量）作為設(shè)計(jì)水文條件，分別核定功能區(qū)的納污能力。

3.1 設(shè)計(jì)流量

山區(qū)性河流水功能區(qū)設(shè)計(jì)流量確定方法：

（1）選用全省代表性好、資料系列長(zhǎng)的流量站，統(tǒng)計(jì)其逐年各月平均流量，以歷年最枯月平均流量系列排頻，得出各站50%、75%及90%保證率的最枯月平均流量。根據(jù)河流、地形等下墊面情況，將全省劃分成32個(gè)區(qū)塊，將所在分區(qū)或臨近分區(qū)內(nèi)的水文站的設(shè)計(jì)流量換算為不同保證率的單位面積流量。

PBL教學(xué)模式已成為國(guó)際上較為流行的一種教學(xué)方法，也逐漸推廣成為我國(guó)醫(yī)學(xué)教育改革的一種主流模式[10-11]，其優(yōu)勢(shì)也被國(guó)內(nèi)許多的醫(yī)學(xué)院校教學(xué)實(shí)踐所證實(shí)。研究表明，運(yùn)用PBL 教學(xué)模式的學(xué)生在進(jìn)入臨床工作后其臨床思辨能力、解決實(shí)際問(wèn)題的能力以及與患者的溝通能力均有顯著提升 [10]。

山區(qū)性河流流量大小取決于上游降水量、下墊面條件、河流形態(tài)、集水面積等多種因素，降水量的大小對(duì)單次流量影響較大，但從長(zhǎng)系列資料來(lái)看，影響流量的主要因素是集水面積。繪制各站50%、75%及90%保證率的最枯月平均流量與集水面積的關(guān)系見(jiàn)圖4。由圖4可見(jiàn)，各站不同保證率流量與集水面積基本成正相關(guān)，隨著集水面積的增加，流量相應(yīng)增加。

（2）在電子地圖上根據(jù)等高線(xiàn)求取各功能區(qū)河段上斷面或下斷面以上匯水面積，并用水文站控制面積進(jìn)行復(fù)核。

（3）各功能區(qū)按所在區(qū)塊的不同保證率單位面積流量和功能區(qū)上斷面或下斷面流域面積或區(qū)間面積計(jì)算設(shè)計(jì)來(lái)水流量，并用水文站控制斷面設(shè)計(jì)流量進(jìn)行協(xié)調(diào)修正后作為納污能力計(jì)算的設(shè)計(jì)流量條件。

圖4 各站不同保證率流量與集水面積關(guān)系圖

3.2 設(shè)計(jì)流速

河流控制斷面的實(shí)測(cè)流與斷面平均流速有如下指數(shù)關(guān)系：

式中：u為斷面平均流速（m／s）；Q為斷面流量（m3／s）；a、b為系數(shù)。

以諸暨站為例，該站中低水實(shí)測(cè)流量～流速關(guān)系及擬合流速見(jiàn)圖5，由圖5可見(jiàn)，流速與流量指數(shù)關(guān)系擬合較好，能基本反映不同流量下的實(shí)際流速情況。

圖5 諸暨站中低水流量～流速關(guān)系擬合圖

根據(jù)全省68個(gè)水文站實(shí)測(cè)中、低水流量資料，逐站分析系數(shù)a、b，選擇相關(guān)關(guān)系好的代表站，分水系綜合后作為各功能區(qū)由流量推算流速的依據(jù)。

3.3 湖庫(kù)及平原河網(wǎng)設(shè)計(jì)水位

湖庫(kù)和平原河網(wǎng)采用50%、75%、90%最枯月平均水位下的蓄水量作為納污能力計(jì)算的水文條件。

全省平原河網(wǎng)分為杭嘉湖平原、蕭紹平原、姚慈平原、鄞東、鄞西、鎮(zhèn)海、金清水系7大河網(wǎng)，分別選取各河網(wǎng)水位代表站，統(tǒng)計(jì)各站歷年逐月平均水位，以歷年最枯月平均水位系列排頻，得出各站50%、75%及90%保證率的最枯月平均水位。根據(jù)各河網(wǎng)的水位～槽蓄量關(guān)系，得出各平原河網(wǎng)不同保證率的蓄水量。

選用的水位代表站見(jiàn)表4。

表4 各平原河網(wǎng)區(qū)選用水位代表站表

3.4 感潮河段設(shè)計(jì)潮量

浙江省水系多為獨(dú)流入海河流，河口段受潮汐影響較大，設(shè)計(jì)流量除上游來(lái)水外，還必須包括潮水量。其設(shè)計(jì)流量的確定方法如下：

（1）劃分潮區(qū)界；

（2）根據(jù)潮位站歷年潮水位資料，計(jì)算各站50%、75%及90%保證率潮差；

（3）根據(jù)河道實(shí)測(cè)斷面、潮差和河段長(zhǎng)度，計(jì)算感潮河段的潮水量。

浙江省感潮河段主要有錢(qián)塘江感潮河段、浦陽(yáng)江感潮河段、甬江感潮河段、椒江感潮河段、甌江感潮河段、飛云江感潮河段、鰲江感潮河段。各感潮河段潮區(qū)界及選用潮位站見(jiàn)表5。

表5 全省感潮河段潮區(qū)界及選用潮位站表

根據(jù)感潮河段選用潮位站50%，75%，90%潮差，結(jié)合功能區(qū)斷面、長(zhǎng)度資料計(jì)算各功能區(qū)的潮水量，作為納污能力計(jì)算的設(shè)計(jì)潮量。

4 計(jì)算分區(qū)及模型選用

根據(jù)全省水系狀況結(jié)合行政分區(qū)和水資源分區(qū)，納污能力分析將全省水功能區(qū)分為10個(gè)分區(qū)，分別為錢(qián)塘江區(qū)、杭嘉湖區(qū)、苕溪區(qū)、曹娥江區(qū)、甬江區(qū)、舟山區(qū)、椒江區(qū)、甌江區(qū)、飛云江區(qū)、鰲江區(qū)，獨(dú)流入海及出省小河流按地理屬性劃入上述10個(gè)分區(qū)內(nèi)。根據(jù)各計(jì)算分區(qū)內(nèi)的水域特性，納污能力計(jì)算選用4種模型，分別為山區(qū)性河流一維模型、湖泊水庫(kù)零維模型、平原河網(wǎng)零維模型、感潮河段混合模型。各計(jì)算分區(qū)不同水功能區(qū)選用模型見(jiàn)表6。

表6 各計(jì)算分區(qū)不同水功能區(qū)選用模型表

續(xù)表6

5 納污能力分析成果

根據(jù)浙江省水功能區(qū)選用的納污能力計(jì)算模型及相關(guān)參數(shù)、設(shè)計(jì)水文條件、功能區(qū)水質(zhì)目標(biāo)等，分析計(jì)算功能區(qū)納污能力 （見(jiàn)表7）。由表7可見(jiàn)，50%設(shè)計(jì)水文條件下，浙江省CODcr納污能力為150.25萬(wàn)t／a，其中感潮河段、平原河網(wǎng)、山區(qū)性河流分別為51.46萬(wàn)，17.14萬(wàn)，81.64萬(wàn)t／a，分別占全省的34.3%，11.4%，54.3%；NH3－N納污能力為6.13萬(wàn)t／a，其中感潮河段、平原河網(wǎng)、山區(qū)性河流分別為2.45萬(wàn)，0.90萬(wàn)，2.78萬(wàn)t／a，分別占全省的40.0%，14.6%，45.4%。

75%設(shè)計(jì)水文條件下，全省CODcr納污能力為123.44萬(wàn)t／a，其中感潮河段、平原河網(wǎng)、山區(qū)性河流分別為43.95萬(wàn)，15.50萬(wàn)，64.00萬(wàn)t／a，分別占全省的35.6%，12.6%，51.8%；NH3－N納污能力為4.95萬(wàn)t／a，其中感潮河段、平原河網(wǎng)、山區(qū)性河流分別為2.01萬(wàn)，0.81萬(wàn)，2.12萬(wàn)t／a，分別占全省的40.7%，16.5%，42.8%。

90%設(shè)計(jì)水文條件下，全省CODcr納污能力為102.86萬(wàn)t／a，其中感潮河段、平原河網(wǎng)、山區(qū)性河流分別為36.77萬(wàn)，14.24萬(wàn)，51.84萬(wàn)t／a，分別占全省的35.8%，13.8%，50.4%；NH3－N納污能力為4.19萬(wàn)t／a，其中感潮河段、平原河網(wǎng)、山區(qū)性河流分別為1.73萬(wàn)，0.75萬(wàn)，1.71萬(wàn)t／a，分別占全省的41.3%，18.0%，40.7%。

表7 各水系水功能區(qū)納污能力計(jì)算成果表

6 結(jié) 語(yǔ)

（1）根據(jù)浙江省水功能區(qū)涉及的水域類(lèi)型，納污能力計(jì)算模型主要可選用山區(qū)性河流一維均勻混合模型、湖泊和水庫(kù)零維模型、平原河網(wǎng)零維模型、感潮河段組合模型4類(lèi)，模型參數(shù)主要有綜合衰減系數(shù)k、湖庫(kù)不均勻系數(shù)B、感潮河段不均勻系數(shù)b，可根據(jù)不同功能區(qū)特性來(lái)確定。

（2）為反映不同設(shè)計(jì)水文條件下功能區(qū)的納污能力，從而為下一步污染物總量控制提供較為全面的科學(xué)數(shù)據(jù)，確定90%、75%、50%三種保證率的設(shè)計(jì)水文條件，分別計(jì)算全省水功能區(qū)的納污能力。

（3）根據(jù)全省納污能力分析計(jì)算成果，感潮河段、平原河網(wǎng)、山區(qū)性河流CODcr納污能力的分配情況大致為0.35∶0.13∶0.52；NH3－N納污能力的分配情況大致為0.41∶0.16∶0.43。

（4）功能區(qū)納污能力與該功能區(qū)的水質(zhì)目標(biāo)、設(shè)計(jì)水文條件等密切相關(guān)。隨著浙江省水利工程的變化，如河道疏浚、水庫(kù)建設(shè)、引調(diào)水工程的實(shí)施等，水文條件將相應(yīng)變化；同時(shí)由于經(jīng)濟(jì)社會(huì)發(fā)展對(duì)水功能區(qū)要求的變化，功能區(qū)水質(zhì)目標(biāo)也會(huì)有所調(diào)整。納污能力應(yīng)根據(jù)相應(yīng)條件的變化重新核定，以利于污染物總量控制方案的科學(xué)實(shí)施。