毛 志，陳亞麗

（ 銅仁學(xué)院 大數(shù)據(jù)學(xué)院，貴州 銅仁 554300 ）

基于貝葉斯公式的河流水質(zhì)綜合評價(jià)
——以銅仁市錦江河為例

毛 志，陳亞麗

（ 銅仁學(xué)院 大數(shù)據(jù)學(xué)院，貴州 銅仁 554300 ）

以《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838-2002）為評價(jià)標(biāo)準(zhǔn)，選取10個水質(zhì)評價(jià)指標(biāo)，根據(jù)銅仁市錦江河13個監(jiān)測斷面多年的水質(zhì)監(jiān)測數(shù)據(jù)平均值，利用貝葉斯公式對銅仁市錦江河進(jìn)行水質(zhì)評價(jià)。通過評價(jià)結(jié)果可以看出在多年平均情況下除錦江東門橋下斷面符合Ⅲ類標(biāo)準(zhǔn)外，其余12個斷面水質(zhì)均優(yōu)于規(guī)定的Ⅲ類標(biāo)準(zhǔn)，符合Ⅱ類水質(zhì)標(biāo)準(zhǔn)。

貝葉斯公式； 河流水質(zhì)評價(jià)； 錦江河

水是生命之源，人們?nèi)粘Ｉ钣盟芏嗍峭ㄟ^河流中的水凈化而來。錦江河是銅仁市碧江區(qū)最大河流，全流域面積為4157平方公里，境內(nèi)4069平方公里，境外湖南省88平方公里，東至辰溪縣城郊入沅江，經(jīng)洞庭湖后匯入長江，是長江上游河流的上游支流[1]。因此，其水質(zhì)是非常重要的，于是尋找出一種簡便、有效的評價(jià)錦江河水質(zhì)的方法成為一項(xiàng)十分迫切而且具有重要意義的工作[2]。由于在水質(zhì)評價(jià)中的不確定因素較多，所以在水質(zhì)管理中，水環(huán)境的不確定性是一個難題。近些年，貝葉斯統(tǒng)計(jì)方法越來越多地應(yīng)用到水環(huán)境不確定性問題研究中，其優(yōu)勢在于將水環(huán)境的不確定性問題轉(zhuǎn)化成對模型參數(shù)的估計(jì)量，建立不確定性和估計(jì)概率之間的關(guān)聯(lián)，以聯(lián)合后驗(yàn)分布的數(shù)值代表不確定事件發(fā)生的概率大小，根據(jù)估計(jì)概率的大小，決策者可以對不確定性問題進(jìn)行評估和判斷，從而有針對性地開展水質(zhì)管理工作，提高水質(zhì)管理的準(zhǔn)確度和有效性。本論文旨在通過貝葉斯公式對銅仁市錦江河的水質(zhì)進(jìn)行綜合評價(jià)，同時給水務(wù)和環(huán)保等部門的管理和決策提供科學(xué)參考。

1.數(shù)據(jù)采集

根據(jù)不同的用水要求，我國已制定頒布了不同的水資源質(zhì)量評價(jià)標(biāo)準(zhǔn)，本研究采用《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838-2002）[3]。結(jié)合水質(zhì)質(zhì)量評價(jià)選擇內(nèi)容和實(shí)際監(jiān)測數(shù)據(jù)情況，我們選用總氮、化學(xué)需氧量、高錳酸鹽指數(shù)、五日生化需氧量、氨氮、總磷、揮發(fā)酚、汞、砷、鉛共10個水質(zhì)評價(jià)指標(biāo)[4]。

監(jiān)測站點(diǎn)的布設(shè)在總體和宏觀上必須能反映出水系或所在區(qū)域的水環(huán)境質(zhì)量狀況[5]。各監(jiān)測站點(diǎn)的具體位置要能反映所在區(qū)域環(huán)境的污染特征，要盡可能以最少的監(jiān)測站點(diǎn)獲取最有足夠代表性的環(huán)境信息，而且還要考慮實(shí)際采樣時的可行性和方便性。因此監(jiān)測站點(diǎn)選取的位置應(yīng)避開死水區(qū)、回水區(qū)、排污口處，選擇順直河段、河床穩(wěn)定、水流平穩(wěn)或水面寬闊、無急流、無淺灘處[6]。綜合以上的條件，我們把本文所研究的錦江河分為大江段、小江段和錦江段，選取了13個監(jiān)測站點(diǎn)（如圖1所示），各監(jiān)測站點(diǎn)的地名及其編號見表1。通過銅仁市環(huán)保局，獲得各監(jiān)測站點(diǎn)水質(zhì)監(jiān)測指標(biāo)的多年平均值，

其中化學(xué)需氧量、氨氮、高錳酸鹽指數(shù)的測量值是2008年至2014年內(nèi)的平均值[7]，其他的水質(zhì)監(jiān)測指標(biāo)測量值是2010年至2014年內(nèi)的平均值。

圖1 錦江河水質(zhì)監(jiān)測斷面示意圖

表1 各監(jiān)測站點(diǎn)的地名及其編號

2.地表水水質(zhì)評價(jià)的貝葉斯模型及方法

2.1.地表水水質(zhì)評價(jià)的貝葉斯模型

根據(jù)水質(zhì)評價(jià)的情況，貝葉斯公式可改寫為：

式中yji表示水質(zhì)類型Bi的值，即水質(zhì)級別為i的j指標(biāo)的標(biāo)準(zhǔn)值，xjk表示站點(diǎn)k的j指標(biāo)的測量值，i為水質(zhì)級別，j為指標(biāo)，k為代表站點(diǎn)，在本論文中上式即為銅仁市碧江區(qū)錦江河水質(zhì)評價(jià)依據(jù)的貝葉斯模型。

2.2.地表水水質(zhì)評價(jià)的貝葉斯方法

（1）計(jì)算P( yji)。在沒有任何水質(zhì)信息的條件下，某站點(diǎn)的水質(zhì)屬于某水質(zhì)級別最為人們接受的

原則是概率相等，即先驗(yàn)概率取

（2）P( xjk|yji)采用幾何概率中的距離法計(jì)算。以代表站點(diǎn)指標(biāo)與水質(zhì)類型標(biāo)準(zhǔn)之間的距離絕對值的倒數(shù)進(jìn)行計(jì)算[8]，即

式中j=1,2,···,10，i=1,2,···,5，k=1,2,···,13。如果代表站點(diǎn)的xji距yji越遠(yuǎn)，則代表站點(diǎn)k的水質(zhì)級別屬于i的可能性就越小。

（3）多指標(biāo)下綜合水質(zhì)后驗(yàn)概率的求解如下：

式中wj為不同水質(zhì)指標(biāo)的權(quán)重，因本文是綜合研究，所以就采用等權(quán)重[9]，即w1=w2=···=w10=1/10，Pi是綜合考慮了10個水質(zhì)指標(biāo)信息后代表站點(diǎn)的水質(zhì)級別為i的概率。

（4）最終的級別h是以最大概率原則決定，即

3.基于貝葉斯模型的錦江河河流水質(zhì)評價(jià)

根據(jù)水質(zhì)監(jiān)測站點(diǎn)的區(qū)域分布和監(jiān)測資料，本文選取10個水質(zhì)指標(biāo)，利用銅仁市錦江河13個監(jiān)測站點(diǎn)多年的水質(zhì)監(jiān)測數(shù)據(jù)平均值[10]，以《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838-2002）作為標(biāo)準(zhǔn)，通過上述步驟進(jìn)行計(jì)算，得到P( xjk|yji)及Pi值，Pi的值見表2（Pi保留三位小數(shù)）。由以上的各個站點(diǎn)Pi的值根據(jù)步驟(4)以概率最大原則取五個Pi的最大值（Ph）可以得到13個監(jiān)測站點(diǎn)的水質(zhì)級別（見表3）[11]。

4.結(jié)論與建議

根據(jù)《地表水環(huán)境質(zhì)量評價(jià)辦法》（GB3838-2002）要求，通過貝葉斯水質(zhì)評價(jià)模型，計(jì)算結(jié)果表明銅仁市錦江河的13個代表站點(diǎn)多年來的平均水質(zhì)評價(jià)結(jié)果為Ⅱ類、Ⅲ類，沒有Ⅰ類、Ⅳ類和Ⅴ類水質(zhì)。在錦江河的13個監(jiān)測站點(diǎn)中，除錦江東門橋下屬于Ⅲ類水質(zhì)外，其余各監(jiān)測斷面都屬于Ⅱ類水質(zhì)標(biāo)準(zhǔn)。整體來說，錦江河的水質(zhì)達(dá)標(biāo)，且通過錦江東門橋旁邊的污水處理廠的處理后，錦江東關(guān)碼頭和大明邊城下游的水質(zhì)有所改善。但通過計(jì)算發(fā)現(xiàn)大江的木杉河楚溪村上游、龍洞河匯入木杉河前200米、木杉河匯入大江前200米和大江坪橋下的水質(zhì)指標(biāo)中總氮含量接近Ⅳ、Ⅴ類的總氮標(biāo)準(zhǔn)值的概率較大，而且錦江的三個監(jiān)測站點(diǎn)的總氮接近Ⅳ類的總氮標(biāo)準(zhǔn)值的概率非常大。

表2 各監(jiān)測站點(diǎn)的不同水質(zhì)級別的概率

表3 各監(jiān)測站點(diǎn)水質(zhì)級別的最大概率hP的值及相應(yīng)的水質(zhì)級別

通過貝葉斯公式對錦江河的水質(zhì)得出的評價(jià)結(jié)果表明錦江河的河水符合地表水水質(zhì)標(biāo)準(zhǔn)，但還應(yīng)加強(qiáng)水質(zhì)的檢測，特別注意保護(hù)東門橋下的河水。此外，大江的木杉河楚溪村上游、龍洞河匯入木杉河前200米、木杉河匯入大江前200米和大江坪橋下和錦江監(jiān)測斷面的總氮含量偏高，建議拒絕污染重、耗能高、耗水大的項(xiàng)目進(jìn)入該區(qū)域。通過錦江東門橋旁邊的污水處理廠的處理后，錦江東關(guān)碼頭和大明邊城下游的水質(zhì)有所改善?？梢?，污水處理廠的建設(shè)是非常有必要的，應(yīng)加快污水處理廠的建設(shè)及垃圾處置設(shè)施建設(shè)。

[1] 夏鵬章.長江志·水系篇編撰記述[J].長江志季刊，2003，（2）：41-45.

[2] MO Angelidis, G Kamizoulis. A Rapid Decision-Making Method for the Evaluationof Pollution-Sensitive Coastal Areas in the Mediterranean Sea[J]. Environmental Management, 2005, 35 (6): 811-820.

[3] 齊培培.基于貝葉斯網(wǎng)絡(luò)的水質(zhì)污染評價(jià)及預(yù)測[D].武漢：武漢大學(xué)，2009.

[4] 李貴寶，郝紅，張燕.我國水環(huán)境質(zhì)量標(biāo)準(zhǔn)的發(fā)展[J].水利技術(shù)監(jiān)督，2003，（3）.

[5] 劉紅亮.淺談水環(huán)境監(jiān)測的質(zhì)量保證與質(zhì)量控制[J].水文，2012，（10）：320-320.

[6] 李德榮，劉鶴.關(guān)于全概率公式和貝葉斯公式的一種新講解[J].內(nèi)蒙古統(tǒng)計(jì)，2009，18（6）：514-517.

[7] RM Singh, B Datta. Artificial neural network modeling for identification of unknown pollution sources in groundwater with partially missing concentration observation data[J]. Water Resources Management, 2007, 21(3): 557-572.

[8] GF Cooper, E Herskovits. A Bayesian Method for the Induction of Probabilistic Networks from Data[J]. Machine Learning, 1992, 9 (4): 309-347.

[9] 何鎮(zhèn)邦，李桂榮.概率論與數(shù)理統(tǒng)計(jì)[M].北京：北京理工大學(xué)出版社，1998.

[10] 郭敬輝.川西滇北地區(qū)水文地理[M].北京：科學(xué)出版社，1985.

[11] 趙英，郭亮.貝葉斯方法的污染源季節(jié)性排放量控制和管理[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào)，2014，46（12）：26-32.

The Comprehensive Assessment of Water Quality Based on Bayesian Formula—Taking Jinjiang River as an Example

Mao Zhi, Chen Yali
( School of Data Sciences, Tongren University, Tongren, Guizhou 411105, China )

In this paper, with Standard GB3838-2002 as the evaluation criteria, we select 10 water quality evaluation index，and use the Bayesian formula to assess the water quality of Jinjiang river according to the mean of monitoring data in 13 monitoring sections. The evaluation result shows that the East Gate Bridge section meets Ⅲ class water quality standard and the remaining 12 section water quality get toⅡclass standard.

Bayesian formula, river water quality assessment, Jinjiang river

O211；X824

A

1673-9639 (2016) 04-0150-04

（責(zé)任編輯 佘彥超）（責(zé)任校對 印有家）

2016-05-10

貴州省普通高等學(xué)校創(chuàng)新人才團(tuán)隊(duì)項(xiàng)目（黔教合人才團(tuán)隊(duì)字[2015]64）。

毛 志（1982-），男，湖南汨羅人，副教授，博士，研究方向：微分方程數(shù)值計(jì)算，生物信息學(xué)。陳亞麗（1991-），女，貴州德江人，銅仁學(xué)院數(shù)學(xué)科學(xué)學(xué)院2012級數(shù)學(xué)與應(yīng)用數(shù)學(xué)（2）班學(xué)生。