鄭志國　張秀敏　岳　超　趙　鵬　李　斌

?

模糊數(shù)學法在水質評價中的應用

鄭志國1張秀敏2岳超1趙鵬1李斌1

對水功能區(qū)水質進行客觀、真實和合理的評價是實行最嚴格水資源管理的重要考核依據(jù)，由于水功能區(qū)水質污染是由排入水體的各種污染物共同作用的結果，因此對水功能區(qū)水質的評價應是多因子的綜合評價。本文根據(jù)地表水環(huán)境質量標準，給出了用模糊數(shù)學綜合評判法進行水功能區(qū)水質評價，其結果令人滿意。

一、模糊數(shù)學綜合評判模型的建立

設立兩個有限論域

其中U是有影響水質的m個因素組成的論域（因素集）；V是由n個評價級別鎖組成的論域（評價集）。

令rij表示因素ui子抉擇級別vi上的可能性（即ui對vi的隸屬度），則論域V上的模糊子集：

表示第i個因素單因子決策矩陣。而m×n模糊矩陣：

則表示m個單因素總決策矩陣。

表示了評判對象的模糊綜合評價結果。其中“·”表示某種合成運算。將做歸一化處理得=（b1'，b2'，…bn）'

由最大隸屬度原則，評判對象定為第j類水質。在式（4）中，運算“·”要根據(jù)具體的問題做一定的選擇。設置為加權平均型模型：

這種合成運算對所有因素依權重大小均衡兼顧，在考慮總體因素時效果明顯，比較適合用于要求整體指標的情形。

二、評價算法

設用《地表水環(huán)境質量標準》（GB3838-2002）等級標準對T（T表示被評價水質的某個參數(shù)）項目進行評價，標準中等級數(shù)為Gk，k=1，2，3，…，s，即有s個等級，假設某水質有m個評價因素（參數(shù)）uj，j=1，2，3，…，m。每個評價參數(shù)有n個定性的評價等級Vi=1，2，3，…，n。這些等級按評價要求具體劃分，可以定為I、II、III、IV、V等級別。

對照標準，可以確定某水質的每個評價參數(shù)uj所在的評價等級標準，記為Rmn，得到的評價表格如表1。

表1反映了各單項參數(shù)與等級之間的關系，這種關系用隸屬度表示稱作模糊關系。表中Rij表示被評價水質的第j個因素（參數(shù)）u可能為等級Vi的概率（即隸屬度）。用模糊矩陣R表示：

設被評價水質的參數(shù)評價矩陣為B，則B=A·R即：

A與B是兩個模糊矩陣，所以以上的矩陣的運算遵循模糊矩陣的復合運算法。得：

表1　評價表

B矩陣表示水質中的某評價中屬于V1等級的程度（比例）是b1，屬于V2等級的程度是b2，…，依次類推。根據(jù)矩陣B可以綜合評價水質所屬的等級。

三、隸屬度

根據(jù)《地表水環(huán)境質量標準》（GB3838-2002），水質類別共劃分為五類。利用線性函數(shù)來建立各類別的隸屬rij。設：

xi是第i種污染因素的實測濃度值（mg/L）；

Sij是第i種污染因素的第j類水質的標準濃度值（mg/L）；

Rij是第i種污染因素對第j類水質的隸屬函數(shù)值。

i=1，2…，m；j=1，2，…n。

當j=1時：

四、權重及歸一化運算

根據(jù)各參數(shù)超標情況進行加權，超標越多，加權越大。權重值為：

式中：Wi—第i種污染物以平均標準為基準的超標指數(shù)，即為權重；

Ci—第i種污染物實測濃度；

Si—第i種污染物各級標準值的算術平均值。

為進行模糊運算，將各單項權重再進行歸一化運算：

式中：Vi—第i種污染物的歸一化權重；Ci—同上；Si—同上。

五、計算評價結果

取2014年3月大汶河泰安工業(yè)用水區(qū)大汶口斷面三個參數(shù)來評價，水質分為五個等級。

1.用隸屬度刻畫水質分級界限

以氨氮（NH3-N）為例，監(jiān)測值為：X=0.83，其相鄰地兩類水質標準值為：X0=0.5（II類），X1=1.0（Ⅲ類）。

對II類水的隸屬度為：

即有34%可能劃為II類水，66%可能劃為Ⅲ類水。很明顯，就氨氮（NH3-N）單項污染物而言，該水質不可能劃為Ⅰ、IV、V類，故它們地隸屬度為零。同樣可得化學需氧量（CODcr）和生化需氧量（BOD5）對各等級水得隸屬度并構成m×n＝3×5的模糊矩陣R。

2.計算權重

計算權重并賦予各參數(shù)：對上述集合中U中m項參數(shù)給予權重，組成一個1×m矩陣A。

對氨氮（NH3-N）污染物各級標準值的算術平均值為：

此值介于III、IV級之間，同樣可得化學需氧量（CODcr）和生化需氧量（BOD5）的各級標準值的算術平均值SCODcr=24、SBOD5=5.2。

對氨氮（NH3-N）污染物的超標指數(shù)，即權重：

同樣可得化學需氧量（CODcr）和生化需氧量（BOD5）的權重分別為WCODcr=1. 76、WBOD5=1.06。

對氨氮（NH3-N）污染物的歸一化權重VNH3-N：

同樣可得化學需氧量（CODcr）和生化需氧量（BOD5）的歸一化權重分別為VCODcr=0.485、VBOD5=0.292。

各污染物歸一化權重計算結果見表2，表中：Ci—同上；Si—同上，該例為五類標準的平均值。

由此算出的各污染物的權重，構成一個1×3的行矩陣：A=（0.223，0.485，0.292）。

3.模糊矩陣復合運算及評價結果

為了進行綜合評價，將上述A和R矩陣進行復合運算。通過經(jīng)復合運算后得出的結果進行綜合評價。

由以上的計算結果可得綜合評價結果是該水質屬Ⅴ類■

表2　各污染物歸一化權重計算表

（作者單位：1.山東省泰安市水文局2710002.山東省濟寧市水文局272000）