段 震 張新明 楊豐云

一、前言

大汶河發(fā)源于山東省萊蕪市鋼城區(qū)黃莊鎮(zhèn)臺(tái)子村，迂回西流，途經(jīng)山東萊蕪市、泰安市，在泰安市東平縣馬口村注入東平湖，經(jīng)陳山口出湖閘出東平湖，進(jìn)入黃河。大汶河全長(zhǎng)208km，流域面積9069km2。大汶河汶口壩以上為大汶河上游，是大汶河的主要集水區(qū)，汶口壩至戴村壩為大汶河中游，戴村壩以下至東平湖為大汶河下游，其中上游主要支流有柴汶河，中下游主要支流有漕濁河和匯河。

東平湖是全省第二大淡水湖，位于大汶河下游東平縣境內(nèi)，上承大汶河來(lái)水，南與運(yùn)河相連，北有小清河與黃河相通，是黃河下游最大的滯洪區(qū)，多年平均水面面積119km2。位于東平縣中部的稻屯洼，三面環(huán)山，既是東平北部山丘區(qū)河流的匯集區(qū)，又擔(dān)負(fù)著大汶河滯洪蓄任務(wù)。

大汶河水系布設(shè)國(guó)家水質(zhì)監(jiān)測(cè)站多處，自上而下設(shè)有萊蕪站、北望站、大汶口站、戴村壩站，大汶河南支柴汶河設(shè)有樓德站，大汶河支流匯河設(shè)有白樓站。

二、評(píng)價(jià)數(shù)學(xué)模型

設(shè)待評(píng)價(jià)分類(lèi)的有n個(gè)污染水體樣本(斷面)，每個(gè)樣本中有m個(gè)污染指標(biāo)，其實(shí)測(cè)值為y'。作為評(píng)價(jià)標(biāo)準(zhǔn)有k個(gè)水質(zhì)功能分類(lèi)，對(duì)應(yīng)于m項(xiàng)污染指標(biāo)，均有國(guó)家規(guī)定的標(biāo)準(zhǔn)值為x'，并設(shè)：

k類(lèi)水質(zhì)標(biāo)準(zhǔn)污染物的特征值：

n個(gè)污染水體樣本的污染指數(shù)：

評(píng)價(jià)樣本對(duì)k類(lèi)標(biāo)準(zhǔn)之間的符合程度，即貼近度定義為：

式中：dij—某評(píng)價(jià)樣本j對(duì)類(lèi)別i的貼近度；

Σxim—水質(zhì)標(biāo)準(zhǔn)污染物的特征值之和；

Σyjm—樣本j的污染指數(shù)之和。

若滿(mǎn)足：

樣本j的綜合污染指數(shù)與水質(zhì)標(biāo)準(zhǔn)類(lèi)別污染物的特征值相差越小，說(shuō)明越貼近k類(lèi)。從而可判斷評(píng)價(jià)樣本j屬于k類(lèi)。

三、評(píng)價(jià)步驟

1.評(píng)價(jià)標(biāo)準(zhǔn)

依據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838-2002），將河流斷面的水質(zhì)綜合污染指數(shù)與標(biāo)準(zhǔn)類(lèi)別的綜合特征值比較，對(duì)斷面水質(zhì)進(jìn)行評(píng)價(jià)分類(lèi)。

2.評(píng)價(jià)參數(shù)

評(píng)價(jià)參數(shù)對(duì)評(píng)價(jià)結(jié)果影響較大，選用有代表性的參數(shù)極為重要。為使評(píng)價(jià)結(jié)果更為合理地符合水質(zhì)標(biāo)準(zhǔn)類(lèi)別，選用水質(zhì)標(biāo)準(zhǔn)建議選用的參數(shù)，且在水質(zhì)監(jiān)測(cè)中檢出頻率高、超標(biāo)嚴(yán)重的10項(xiàng)指標(biāo)進(jìn)行評(píng)價(jià)。

3.確定水質(zhì)標(biāo)準(zhǔn)的特征值和樣本污染指數(shù)

用標(biāo)準(zhǔn)類(lèi)別同一參數(shù)的平均值作為標(biāo)準(zhǔn)，即：

由各污染物的平均值按下面兩種情況求出水質(zhì)標(biāo)準(zhǔn)的特征值和樣本的污染指數(shù)。

第一種情況：污染物的危害程度隨其濃度的增加而降低的評(píng)價(jià)參數(shù)（如溶解氧），計(jì)算式為：

第二種情況：污染物的危害程度隨其濃度的增加而增加的評(píng)價(jià)參數(shù)，計(jì)算式為：

表1 水質(zhì)監(jiān)測(cè)實(shí)測(cè)值表（mg/L）

表3 各斷面的貼近度與評(píng)價(jià)結(jié)果表

4.求貼近度確定樣本評(píng)價(jià)類(lèi)別

根據(jù)（3）式求出評(píng)價(jià)樣本綜合污染指數(shù)對(duì)k類(lèi)標(biāo)準(zhǔn)特征值的貼近度。在評(píng)價(jià)樣本對(duì)各類(lèi)標(biāo)準(zhǔn)的貼近度中，將最小的貼近度挑選出來(lái)，該貼近度所對(duì)應(yīng)的類(lèi)別即為樣本應(yīng)屬的水質(zhì)類(lèi)別。

四、大汶河水系水質(zhì)監(jiān)測(cè)站水質(zhì)評(píng)價(jià)分析

大汶河流域6個(gè)河道水質(zhì)監(jiān)測(cè)斷面，現(xiàn)舉例分析評(píng)價(jià)水質(zhì)監(jiān)測(cè)成果見(jiàn)表1。

1.確定標(biāo)準(zhǔn)類(lèi)別的特征值

根據(jù)（5）式求出水質(zhì)標(biāo)準(zhǔn)中同一參數(shù)的平均值，計(jì)算結(jié)果見(jiàn)表2最后一行。以平均值為基準(zhǔn)，按（6）（8）式計(jì)算水質(zhì)標(biāo)準(zhǔn)類(lèi)別的特征值。

第一種情況特征值的計(jì)算如DO：

第二種情況特征值的計(jì)算如CODcr：

仿此可求得其他類(lèi)別參數(shù)的特征值。

2.確定斷面實(shí)測(cè)值的污染指數(shù)

以水質(zhì)標(biāo)準(zhǔn)參數(shù)的平均值為基準(zhǔn)，按（7）（9）式確定相應(yīng)實(shí)測(cè)值的污染指數(shù)。以斷面l為例，如DO、CODcr的污染指數(shù)分別為：

仿此可求得其他斷面各測(cè)定值的污染指數(shù)。

3.求貼近度，確定評(píng)價(jià)結(jié)果

水質(zhì)標(biāo)準(zhǔn)各類(lèi)特征值之和如表3右側(cè)一欄所示，各斷面實(shí)測(cè)值的污染指數(shù)之和如表4右側(cè)一欄所示。依據(jù)（3）式分別求出各斷面綜合污染指數(shù)對(duì)k類(lèi)標(biāo)準(zhǔn)特征值的貼近度，依據(jù)（4）式確定各斷面的評(píng)價(jià)結(jié)果。

如斷面1的5個(gè)貼近度分別為：

d11=|Σx1m-Σy1m|=|2.075-5.942|=3.867

d21=|Σx2m-Σy1m|=|3.494-5.942|=2.448

……

仿此可求得其他各斷面的貼近度，結(jié)果見(jiàn)表3。

由評(píng)價(jià)結(jié)果可以看出，萊蕪站控制斷面水質(zhì)較好，其次為戴村壩站、大汶口站、北望站、樓德站；污染最嚴(yán)重的斷面為白樓站，主要超標(biāo)物質(zhì)為氨氮、總磷。各斷面主要超標(biāo)為揮發(fā)酚、氨氮、總磷等，其余7項(xiàng)評(píng)價(jià)指標(biāo)也存在不同程度的污染。

五、結(jié)論

綜合污染指數(shù)法得出的結(jié)論較為合理，與實(shí)際污染狀況較為吻合，不存在人為主觀分類(lèi)現(xiàn)象，比較符合地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)分類(lèi)要求，且概念清晰、直觀，信息利用較完全，這是本法的優(yōu)點(diǎn)，值得推廣■