沱江水質(zhì)模糊綜合評(píng)價(jià)及主要污染物的預(yù)測(cè)研究

符東，吳雪菲，易珍言，陳永燦，3*

（1.西南科技大學(xué)環(huán)境與資源學(xué)院，四川 綿陽(yáng) 621010；2.四川文理學(xué)院化學(xué)化工學(xué)院，四川 達(dá)州 635000；3.清華大學(xué)水沙科學(xué)與水利水電工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100084）

沱江位于四川省中部，是長(zhǎng)江上游的一條支流。沱江自北向南流經(jīng)金堂縣趙鎮(zhèn)，穿過(guò)金堂峽進(jìn)入簡(jiǎn)陽(yáng)市，經(jīng)過(guò)資陽(yáng)市、內(nèi)江市、自貢市，最后在瀘州市匯入長(zhǎng)江。作為四川盆地重要的自然資源，在過(guò)去幾十年，由于受到工農(nóng)業(yè)活動(dòng)的影響，沱江水質(zhì)呈現(xiàn)出了不同程度的惡化[1]。雖然沱江的水資源對(duì)整個(gè)四川盆地的經(jīng)濟(jì)社會(huì)發(fā)展起著重要作用，但是對(duì)整個(gè)沱江的水質(zhì)評(píng)價(jià)、污染源解析及水質(zhì)模擬的綜合研究卻較少。因此，對(duì)沱江水體進(jìn)行水質(zhì)綜合評(píng)價(jià)和預(yù)測(cè)顯得尤為重要。水質(zhì)評(píng)價(jià)及預(yù)測(cè)可以為沱江水環(huán)境綜合治理提供科學(xué)合理的分析，有助于制定水環(huán)境綜合治理針對(duì)性方案，同時(shí)可以模擬水質(zhì)的變化趨勢(shì)，為水環(huán)境規(guī)劃與管理提供科學(xué)支撐。

目前，水質(zhì)評(píng)價(jià)方法主要有模糊綜合評(píng)價(jià)（Fuzzy comprehensive assessment，F(xiàn)CA）法[2]、人工神經(jīng)網(wǎng)絡(luò)（Artificial neural network，ANN）評(píng)價(jià)法[3]、地理信息系統(tǒng)方法[4]以及多種方法的相互耦合[5]等。由于水環(huán)境的連續(xù)性及不確定性，國(guó)內(nèi)外學(xué)者常選擇模糊綜合評(píng)價(jià)法應(yīng)用于地下水和湖庫(kù)河流、海水等地表水環(huán)境的評(píng)價(jià)中[6-10]。例如，徐曉云等[11]利用模糊綜合評(píng)價(jià)法對(duì)京杭大運(yùn)河揚(yáng)州段的水質(zhì)進(jìn)行了評(píng)價(jià)，并分析了水質(zhì)的主要污染來(lái)源。樊慶鋅等[12]利用灰色關(guān)聯(lián)度和主成分分析法優(yōu)化評(píng)價(jià)指標(biāo)，并結(jié)合模糊綜合評(píng)價(jià)對(duì)松花江哈爾濱段的水質(zhì)進(jìn)行了評(píng)價(jià)。向文英等[13]通過(guò)改進(jìn)權(quán)重賦值方法，利用模糊綜合評(píng)價(jià)對(duì)某水庫(kù)水質(zhì)進(jìn)行了評(píng)價(jià)。在水質(zhì)評(píng)價(jià)的基礎(chǔ)上進(jìn)行污染源解析和水質(zhì)預(yù)測(cè)，可以對(duì)主要污染因子的時(shí)空變化趨勢(shì)進(jìn)行分析，掌握水體水質(zhì)的變化狀況。人工神經(jīng)網(wǎng)絡(luò)能夠精確地計(jì)算復(fù)雜的非線性輸入輸出關(guān)系，因此被廣泛用于水質(zhì)的預(yù)測(cè)[14-16]。當(dāng)前水質(zhì)預(yù)測(cè)主要包括以當(dāng)前數(shù)據(jù)預(yù)測(cè)未來(lái)水質(zhì)變化趨勢(shì)[17-18]和以當(dāng)前采樣點(diǎn)數(shù)據(jù)預(yù)測(cè)其他采樣點(diǎn)水質(zhì)變化趨勢(shì)[19-20]兩個(gè)方面。雖然模糊綜合評(píng)價(jià)和人工神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)在水環(huán)境中應(yīng)用廣泛，但大多數(shù)研究對(duì)模糊綜合評(píng)價(jià)中評(píng)價(jià)因子的篩選沒(méi)有詳細(xì)的描述和論證，同時(shí)對(duì)人工神經(jīng)網(wǎng)絡(luò)輸入指標(biāo)的選擇大多都基于主觀性，使得計(jì)算結(jié)果不能真實(shí)反映水體的狀況。鑒于沱江的重要性，本文在前人研究的基礎(chǔ)上，擬采用沱江流域31個(gè)監(jiān)測(cè)斷面2018年1月—2019年10月的逐月水質(zhì)數(shù)據(jù)，通過(guò)篩選評(píng)價(jià)因子對(duì)各斷面使用模糊綜合評(píng)價(jià)方法進(jìn)行水質(zhì)評(píng)價(jià)。在水質(zhì)評(píng)價(jià)的基礎(chǔ)上，利用主成分分析（Principal component analysis，PCA）識(shí)別沱江的主要污染源和污染因子，然后構(gòu)建BP（Back Propagation，BP）神經(jīng)網(wǎng)絡(luò)，并利用內(nèi)江二水廠斷面數(shù)據(jù)預(yù)測(cè)申家溝斷面主要污染因子濃度。以期為沱江的綜合管理提供一定的參考。

1 研究區(qū)域及指標(biāo)

圖1 沱江及各監(jiān)測(cè)斷面位置分布圖Figure 1 Tuo River and location distribution of the monitoring sections

本文選擇了沱江上游至下游的31個(gè)（SS1～SS31）監(jiān)測(cè)斷面，沱江流域及各監(jiān)測(cè)斷面分布如圖1所示。水質(zhì)數(shù)據(jù)來(lái)源于中國(guó)環(huán)境監(jiān)測(cè)總站在各監(jiān)測(cè)斷面2018年1月—2019年10月的逐月監(jiān)測(cè)濃度。水質(zhì)指標(biāo)包括化學(xué)需氧量（CODCr）、溶解氧（DO）、五日生化需氧量（BOD5）、高錳酸鹽指數(shù)（CODMn）、總氮（TN）、總磷（TP）、氨氮（NH3-N）、硝酸鹽（NO-3）、硫酸鹽（SO24-）、氟化物（F-）、氯化物（Cl-）、pH、氰化物（CN-）、硫化物（S）、陰離子表面活性劑（LAS）、石油類(lèi)（Oils）、揮發(fā)酚（V-phen）、電導(dǎo)率（EC）、糞大腸菌群（E.coli）、汞（Hg）、鉛（Pb）、銅（Cu）、鋅（Zn）、鐵（Fe）、錳（Mn）、六價(jià)鉻（Cr6+）、鎘（Cd）和溫度（T）等。為保證所選指標(biāo)能真實(shí)反映水體狀況，需對(duì)模糊綜合評(píng)價(jià)的評(píng)價(jià)因子集進(jìn)行篩選和優(yōu)化。GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中[21]未對(duì)pH、T、、、Cl-、Fe、Mn和EC的濃度做出等級(jí)劃分和明確要求，故將這些指標(biāo)排除。CN-、S、LAS、Oils、V-phen、Hg、Pb、Cu、Zn、Fe、Mn、Cd、Cr6+濃度均符合Ⅰ類(lèi)水質(zhì)標(biāo)準(zhǔn)且濃度長(zhǎng)期穩(wěn)定無(wú)變化，因此將這些指標(biāo)排除。E.coli由于監(jiān)測(cè)次數(shù)較少，數(shù)據(jù)不連續(xù)，所以排除該指標(biāo)。經(jīng)過(guò)篩選后的評(píng)價(jià)因子集包含 CODCr、DO、BOD5、CODMn、TN、NH3-N、TP和F-，各指標(biāo)的統(tǒng)計(jì)分析如表1所示。在模糊綜合評(píng)價(jià)和主成分分析的基礎(chǔ)上，選取主要污染物，利用BP神經(jīng)網(wǎng)絡(luò)對(duì)其進(jìn)行預(yù)測(cè)。構(gòu)建BP網(wǎng)絡(luò)前，利用各指標(biāo)的相關(guān)性，確定BP網(wǎng)絡(luò)的輸入變量。

2 方法及原理

2.1 模糊綜合評(píng)價(jià)模型

2.1.1 評(píng)價(jià)因子集與評(píng)價(jià)標(biāo)準(zhǔn)

沱江流域工業(yè)集中、農(nóng)業(yè)發(fā)達(dá)、食品業(yè)和養(yǎng)殖業(yè)密集。針對(duì)各行業(yè)廢水特點(diǎn)以及檢測(cè)數(shù)據(jù)，兼顧相關(guān)標(biāo)準(zhǔn)的水質(zhì)要求，經(jīng)過(guò)對(duì)28個(gè)物理化學(xué)水質(zhì)參數(shù)篩選，最后使用選定的8個(gè)指標(biāo)構(gòu)建因子集，因子集可以表示為CTP}。評(píng)價(jià)標(biāo)準(zhǔn)依據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》，具體如表1所示。

2.1.2 評(píng)價(jià)因子權(quán)重

模糊綜合評(píng)價(jià)過(guò)程中，需要對(duì)每一個(gè)評(píng)價(jià)因子賦予相應(yīng)的權(quán)重。本文選擇污染因子貢獻(xiàn)率法計(jì)算各評(píng)價(jià)因子的權(quán)重：

式中：xi為污染物i的實(shí)測(cè)濃度；Si為各評(píng)價(jià)因子的第Ⅲ類(lèi)水質(zhì)標(biāo)準(zhǔn)；wi為歸一化后的i因子權(quán)重。評(píng)價(jià)因子的權(quán)重向量可表示為 W={w1，w2，...，wn}，n為評(píng)價(jià)因子個(gè)數(shù)。由于DO屬于數(shù)值越大，水質(zhì)越好，所以對(duì)于 DO，ai=Si∕xi。

2.1.3 隸屬函數(shù)與模糊矩陣

建立各評(píng)價(jià)因子的隸屬函數(shù)，得到其對(duì)每一類(lèi)水質(zhì)的隸屬度，進(jìn)而得到模糊關(guān)系矩陣R。目前，隸屬度一般采用“降半梯形分布法”計(jì)算[22]。對(duì)于數(shù)值越大污染越重的因子，可根據(jù)公式（2）～公式（4）計(jì)算其對(duì)應(yīng)評(píng)價(jià)標(biāo)準(zhǔn)各等級(jí)的隸屬度。

Ⅰ類(lèi)水的隸屬函數(shù)，即j=1時(shí)：

Ⅱ～Ⅳ類(lèi)水的隸屬函數(shù)，即j=2～4時(shí)：

Ⅴ類(lèi)水的隸屬函數(shù)，即j=5時(shí)：

式中：xi為評(píng)價(jià)因子i的實(shí)際監(jiān)測(cè)濃度（i=1，2，…，8）；Sij為第i個(gè)因子的第j類(lèi)水質(zhì)標(biāo)準(zhǔn)（j=1，2，…，5）。由于DO屬于數(shù)值越小，污染越嚴(yán)重，其隸屬函數(shù)如公式（5）～公式（7）所示。

表1 評(píng)價(jià)因子統(tǒng)計(jì)分析及水質(zhì)標(biāo)準(zhǔn)（mg·L-1）Table 1 Statistical analysis of evaluation factors and water quality standards（mg·L-1）

Ⅰ類(lèi)水的隸屬函數(shù)，j=1時(shí)：

Ⅱ～Ⅳ類(lèi)水的隸屬函數(shù)，即j=2～4時(shí)：

Ⅴ類(lèi)水的隸屬函數(shù)，即j=5時(shí)：

將各監(jiān)測(cè)斷面的評(píng)價(jià)因子濃度帶入相應(yīng)的隸屬函數(shù)中，計(jì)算得到不同評(píng)價(jià)標(biāo)準(zhǔn)的隸屬度。然后建立每個(gè)斷面的單因素評(píng)價(jià)矩陣R：

式中：rij為第i個(gè)評(píng)價(jià)因子對(duì)第j類(lèi)水的隸屬度。

2.1.4 模糊綜合評(píng)價(jià)

根據(jù)上述求得的權(quán)重向量和單因素評(píng)價(jià)矩陣，選取相乘相加算子[23]計(jì)算評(píng)價(jià)結(jié)果：B=W·R={b1，b2，b3，b4，b5}，選取bmax對(duì)應(yīng)的水質(zhì)類(lèi)別作為評(píng)價(jià)結(jié)果。

2.1.5 主成分分析

主成分分析是一種數(shù)學(xué)方法，它可以減少水質(zhì)數(shù)據(jù)集的維數(shù)，同時(shí)使數(shù)據(jù)集中有用信息的損失程度較小，進(jìn)而使數(shù)據(jù)集更容易理解。本文對(duì)31個(gè)監(jiān)測(cè)斷面在2018—2019年監(jiān)測(cè)的12種變量（增加Cl-、和EC 4個(gè)指標(biāo)）進(jìn)行了主成分分析，以確定沱江的潛在污染源。

2.2 BP神經(jīng)網(wǎng)絡(luò)模型

2.2.1 網(wǎng)絡(luò)參數(shù)選擇

在構(gòu)建BP神經(jīng)網(wǎng)絡(luò)時(shí)，需對(duì)選用的輸入輸出數(shù)據(jù)進(jìn)行訓(xùn)練，以沱江內(nèi)江二水廠斷面（斷面編號(hào)SS23）的水質(zhì)數(shù)據(jù)對(duì)沱江申家溝斷面（斷面編號(hào)SS24）主要污染指標(biāo)進(jìn)行預(yù)測(cè)。BP網(wǎng)絡(luò)相關(guān)的參數(shù)選擇如下：輸入層與隱含層之間的傳遞函數(shù)為tansig函數(shù)，隱含層與輸出層之間的傳遞函數(shù)為purelin函數(shù)，數(shù)據(jù)歸一化函數(shù)為mapminmax，訓(xùn)練函數(shù)為trainlm，最大學(xué)習(xí)次數(shù)1 000，目標(biāo)誤差為0.000 5，學(xué)習(xí)速率為0.01，其余參數(shù)為默認(rèn)值。

2.2.2 網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)

本研究輸入層神經(jīng)元個(gè)數(shù)的確定方法為，根據(jù)模糊綜合評(píng)價(jià)結(jié)果確定沱江的主要污染因子，然后利用各水質(zhì)參數(shù)與主要污染的相關(guān)性分析確定輸入層神經(jīng)元個(gè)數(shù)。輸出層的神經(jīng)元為沱江流域下游某斷面的主要污染物濃度，所以輸出層神經(jīng)元個(gè)數(shù)為1。輸入層和輸出層數(shù)均為1。對(duì)于含有1個(gè)隱含層的BP神經(jīng)網(wǎng)絡(luò)，其可以逼近任意一個(gè)非線性函數(shù)，所以本研究采用1個(gè)隱含層。隱含層神經(jīng)元個(gè)數(shù)的確定目前沒(méi)有明確的方法，因此本文采取試錯(cuò)法確定隱含層神經(jīng)元的個(gè)數(shù)[24]。

3 結(jié)果與討論

3.1 沱江水質(zhì)模糊綜合評(píng)價(jià)

本文經(jīng)過(guò)指標(biāo)篩選并構(gòu)建評(píng)價(jià)因子，利用公式（1）計(jì)算各監(jiān)測(cè)斷面評(píng)價(jià)因子的歸一化權(quán)重。通過(guò)復(fù)合運(yùn)算B=W·R，以最大隸屬度原則可得出當(dāng)前斷面的水質(zhì)類(lèi)別。以斷面SS2為例，單因素評(píng)價(jià)矩陣為：

斷面SS2的權(quán)重向量為W={0.057，0.076，0.075，0.130，0.060，0.075，0.356，0.171}，則斷面SS2的評(píng)價(jià)結(jié)果為 B={0.365，0.125，0.154，0，0.356}，根據(jù)最大隸屬度原則，評(píng)價(jià)結(jié)果為Ⅰ類(lèi)。沱江各斷面的評(píng)價(jià)結(jié)果如表2所示。

從表2可以看出，沱江31個(gè)監(jiān)測(cè)斷面中有9個(gè)斷面評(píng)價(jià)為Ⅰ類(lèi)水，占比29.03%，其余22個(gè)斷面均為Ⅴ類(lèi)水，占比70.97%。按照現(xiàn)有沱江水質(zhì)功能區(qū)劃分，普遍水質(zhì)標(biāo)準(zhǔn)為Ⅲ類(lèi)水[25]，因此目前沱江已被嚴(yán)重污染。從各監(jiān)測(cè)斷面評(píng)價(jià)結(jié)果可知，斷面SS2、SS3、SS18～SS22、SS30和SS31為Ⅰ類(lèi)水。綜合圖1可以發(fā)現(xiàn)，這部分?jǐn)嗝婢挥谶h(yuǎn)離城郊的農(nóng)村地區(qū)，森林覆蓋率較高，同時(shí)周?chē)鷽](méi)有工業(yè)和服務(wù)業(yè)，受人為活動(dòng)的影響較小。其余22個(gè)斷面為Ⅴ類(lèi)水，上游、中游和下游均有分布。此部分?jǐn)嗝婢哂锌拷袇^(qū)和城鎮(zhèn)、人口密度大、種植業(yè)發(fā)達(dá)等特點(diǎn)，因此受人為活動(dòng)的影響較大。從權(quán)重賦值可以得出，在挑選的8個(gè)水質(zhì)評(píng)價(jià)因子中，TN的濃度很高，幾乎所有斷面的TN濃度都超過(guò)了Ⅴ類(lèi)水的水質(zhì)標(biāo)準(zhǔn)，所以導(dǎo)致沱江屬于Ⅴ類(lèi)水的監(jiān)測(cè)斷面比例較高。但是在沱江9個(gè)Ⅰ類(lèi)水?dāng)嗝嬷?，?個(gè)斷面（SS2、SS3）的TN濃度超過(guò)Ⅴ類(lèi)水水質(zhì)標(biāo)準(zhǔn)，有4個(gè)斷面（SS19、SS20，SS30、SS31）的TN濃度超過(guò)了Ⅳ類(lèi)水水質(zhì)標(biāo)準(zhǔn)。因此，即使該監(jiān)測(cè)斷面評(píng)價(jià)結(jié)果為Ⅰ類(lèi)水，也需對(duì)其中的TN濃度加以控制。從各斷面監(jiān)測(cè)數(shù)據(jù)可以看出，除了TN以外，其他水質(zhì)指標(biāo)均能滿足Ⅲ類(lèi)水水質(zhì)標(biāo)準(zhǔn)。

表2 沱江水質(zhì)模糊綜合評(píng)價(jià)結(jié)果Table 2 Results of fuzzy comprehensive assessment of water quality in Tuo River

3.2 沱江潛在污染源

利用沱江12個(gè)水質(zhì)參數(shù)的主成分分析法對(duì)不同污染物的來(lái)源進(jìn)行識(shí)別。Bartlett球形檢驗(yàn)的值為785，KMO值為0.709（P＜0.01），說(shuō)明數(shù)據(jù)集適合主成分分析。以特征值大于1為基礎(chǔ)，從沱江水質(zhì)數(shù)據(jù)集中提取4個(gè)主成分[26]，如表3所示。

第一個(gè)主成分（F1）對(duì)NO-3、SO2-4、EC和TN具有較強(qiáng)的正載荷，對(duì)TP具有中等的正載荷，對(duì)NH3-N具有較弱的正載荷，F(xiàn)1占總方差的28.263%。根據(jù)模糊綜合評(píng)價(jià)結(jié)果可知，整個(gè)沱江的TN濃度很高。同時(shí)，沱江流域農(nóng)業(yè)相對(duì)發(fā)達(dá)，氮和磷可以通過(guò)氮肥和磷肥的徑流進(jìn)入河流[27-30]?？赡軄?lái)自礦物巖石[31]，但是在所有斷面的濃度低且較為穩(wěn)定，因此綜合氮和磷的來(lái)源分析，F(xiàn)1可以看作是農(nóng)業(yè)面源污染。第二個(gè)主成分（F2）對(duì)CODMn有強(qiáng)的正載荷，對(duì)Cl-和TP有中等的正載荷，對(duì)DO有中等的負(fù)載荷，F(xiàn)2解釋了總方差的19.644%。根據(jù)CODMn和DO的載荷，F(xiàn)2可能是一種耗氧污染源[32]。同時(shí)，沱江流域的化工廠和制造業(yè)較多，Cl-可能來(lái)自工業(yè)廢水。因此F2可能代表工業(yè)廢水污染[29]。第三個(gè)主成分（F3）占總變異量的16.987%，對(duì)CODCr和BOD5的正載荷較大，對(duì)NH3-N的負(fù)載荷較小。F3包括營(yíng)養(yǎng)物質(zhì)和有機(jī)污染物，可歸因于生活污水和工業(yè)廢水[33]。由于大多數(shù)采樣點(diǎn)的大腸桿菌濃度很高，基本超過(guò)地表水Ⅴ類(lèi)標(biāo)準(zhǔn)（由于數(shù)據(jù)的不連續(xù)性，未對(duì)大腸桿菌濃度做其他分析），因此F3更有可能是生活污水污染源。第四個(gè)主成分（F4）對(duì)F-有較強(qiáng)的正載荷，對(duì)CODMn有較弱的負(fù)載荷，解釋了總方差的9.195%。F-通常來(lái)自氟化工廠、水泥廠和冶煉廠，但沱江各采樣點(diǎn)F-的平均濃度均未超過(guò)地表水Ⅰ類(lèi)水質(zhì)標(biāo)準(zhǔn)值，污染程度幾乎為零或極低，因此沱江中的F-可能是當(dāng)?shù)赝寥离S徑流進(jìn)入的結(jié)果[34]。所以F4可以看作是土壤風(fēng)化。綜合以上分析，可以確認(rèn)農(nóng)業(yè)面源和工業(yè)廢水是沱江的主要污染源[35]。

表3 12種指標(biāo)的旋轉(zhuǎn)因子載荷矩陣Table 3 Loadings of 12 measured variables on VARIMAX rotated factors

3.3 沱江水質(zhì)預(yù)測(cè)

根據(jù)模糊綜合評(píng)價(jià)和主成分分析結(jié)果，選取TN為預(yù)測(cè)指標(biāo)，即輸出層為1個(gè)神經(jīng)元。將上述用于主成分分析的12個(gè)指標(biāo)作為BP神經(jīng)網(wǎng)絡(luò)的輸入備選集。數(shù)據(jù)統(tǒng)計(jì)技術(shù)可以幫助確定BP神經(jīng)網(wǎng)絡(luò)的輸入?yún)?shù)[36-37]。本文通過(guò)相關(guān)性分析，選取與TN濃度顯著相關(guān)且相關(guān)系數(shù)大于0.3的指標(biāo)作為最終的BP輸入?yún)?shù)，具體如表4所示。最終選擇CODCr、NH3-N、NO-3、SO2-4、EC和TN為輸入變量，即利用內(nèi)江二水廠斷面的CODCr、NH3-N、NO-3、SO2-4、EC和TN濃度預(yù)測(cè)申家溝斷面的TN濃度。通過(guò)試錯(cuò)法確定了隱含層的神經(jīng)元個(gè)數(shù)為4，所以BP網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)為6-4-1。根據(jù)相關(guān)系數(shù)r、決定系數(shù)R2和相對(duì)誤差評(píng)價(jià)模型預(yù)測(cè)精度[38-40]，模型預(yù)測(cè)結(jié)果及相對(duì)誤差如圖2和圖3所示。

表4 各指標(biāo)的相關(guān)性Table 4 The correlation of each index

圖2 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值與實(shí)際值的散點(diǎn)圖Figure 2 The scatter plot of predicted value and measured value of BP neural network

圖3 TN線性回歸方程Figure 3 Linear regression equation of TN

BP神經(jīng)網(wǎng)絡(luò)達(dá)到收斂時(shí)的MSE為0.000 49。另外從圖2和圖3可以看出，本文構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)在預(yù)測(cè)TN濃度時(shí)的決定系數(shù)R2為0.999，模型在訓(xùn)練、測(cè)試和驗(yàn)證階段的相對(duì)誤差均小于1.5%。雖然在預(yù)測(cè)階段的相對(duì)誤差有明顯增大，但均小于3%，這有可能與用于BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練的水質(zhì)樣本數(shù)量較少有關(guān)。同時(shí)利用SPSS計(jì)算實(shí)測(cè)值與預(yù)測(cè)值的Pearson相關(guān)系數(shù)，其結(jié)果表明TN的實(shí)測(cè)值與預(yù)測(cè)值極顯著相關(guān)（r=0.99，P＜0.01）。以上評(píng)價(jià)參數(shù)都表明構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)符合水質(zhì)預(yù)測(cè)要求，在樣本較少的情況下能夠?qū)︺|(zhì)進(jìn)行有效預(yù)測(cè)。

4 結(jié)論

（1）沱江31個(gè)評(píng)價(jià)斷面中有9個(gè)斷面水質(zhì)符合Ⅰ類(lèi)水，占29.03%，其余22個(gè)斷面均為Ⅴ類(lèi)水，占比70.97%。Ⅴ類(lèi)水水質(zhì)斷面在沱江上游、中游和下游均有分布，表明沱江整體污染較為嚴(yán)重。

（2）沱江各監(jiān)測(cè)斷面TN濃度均超過(guò)Ⅳ類(lèi)水質(zhì)標(biāo)準(zhǔn)，其中27個(gè)監(jiān)測(cè)斷面超過(guò)Ⅴ類(lèi)水質(zhì)標(biāo)準(zhǔn)，說(shuō)明沱江的主要污染物為T(mén)N。通過(guò)主成分分析，確定沱江TN的主要來(lái)源為農(nóng)業(yè)面源和工業(yè)廢水。

（3）通過(guò)對(duì)輸入數(shù)據(jù)篩選，構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)性能較好，預(yù)測(cè)精度高。在對(duì)沱江申家溝斷面TN濃度的預(yù)測(cè)時(shí)，平均相對(duì)誤差為2.041%?；诒疚牡膬?nèi)容，可構(gòu)建沱江其他斷面的BP神經(jīng)網(wǎng)絡(luò)模型，并用于以沱江上游水質(zhì)數(shù)據(jù)預(yù)測(cè)下游水質(zhì)斷面的TN濃度。