張 瑩，謝仕義＊，楊 鋒

（1.廣東海洋大學(xué) 海洋遙感與信息技術(shù)實驗室，廣東 湛江524088；2.湛江市海洋與漁業(yè)環(huán)境監(jiān)測站，廣東 湛江524039）

粵西海域水質(zhì)評價方法的研究及應(yīng)用＊

張 瑩1，謝仕義1＊，楊 鋒2

（1.廣東海洋大學(xué) 海洋遙感與信息技術(shù)實驗室，廣東 湛江524088；2.湛江市海洋與漁業(yè)環(huán)境監(jiān)測站，廣東 湛江524039）

針對海洋水質(zhì)評價問題，建立了一套先利用主成分分析法明確各水質(zhì)指標(biāo)權(quán)重關(guān)系，繼而參考權(quán)重關(guān)系依照海水水質(zhì)分類標(biāo)準(zhǔn)建立分類樣本，最終通過判別分析方法建立判別函數(shù)實現(xiàn)水體自動分類的方案。該方案在粵西海域，水體分類回判符合率達(dá)到91.4%。依據(jù)自動分類結(jié)果繪制了粵西海域2008—2010年水體分布專題圖，為粵西海域環(huán)境污染防治提供了決策依據(jù)。

粵西海域；水體分類；主成分分析；判別分析

（陳 靖 編輯）

海洋水質(zhì)評價是海洋環(huán)境管理和決策的重要組成部分，其評價內(nèi)容是根據(jù)海水水體的主要物質(zhì)成分及其含量，按照《國家海水水質(zhì)標(biāo)準(zhǔn)》方案，分析水質(zhì)的時空分布狀況，從而為海洋養(yǎng)殖業(yè)、臨海工業(yè)及地方經(jīng)濟(jì)的發(fā)展規(guī)劃及管理決策提供依據(jù)。

從1965年Jacobs提出水質(zhì)指數(shù)概念以及相應(yīng)公式［1］以來，有關(guān)各類水體評價方法的研究成果不斷涌現(xiàn)。早期主要有單因子指數(shù)法［2－4］、綜合指數(shù)法［5－6］、分級評價法［7］等，近期主要有模糊理論、灰色系統(tǒng)理論、投影尋蹤及神經(jīng)網(wǎng)絡(luò)等方法［8－11］。但由于水質(zhì)系統(tǒng)是一個隨時空變化的多維因子組成的復(fù)雜系統(tǒng)，依據(jù)不同區(qū)域，不同時期水體情況建立起來的評價模型很難通用，因而至今國內(nèi)外還沒有統(tǒng)一的水質(zhì)評價模型與方法。

我們針對南?；浳骱Ｓ蛩|(zhì)實測資料，為完成海水水質(zhì)評價工作，設(shè)計下述方案：首先，利用主成分分析法明確海水水質(zhì)各指標(biāo)的權(quán)重關(guān)系；其次，依據(jù)權(quán)重關(guān)系按照《國家海水水質(zhì)標(biāo)準(zhǔn)》建立水體分類樣本；最后，利用判別分析方法分析樣本數(shù)據(jù)，建立水質(zhì)自動分類模型。利用上述方案對粵西海域2010年35個站位水質(zhì)實測資料進(jìn)行了處理，其數(shù)點回判正確率達(dá)到91.4%。同時利用建立的自動分類模型對另外10個站位的水質(zhì)實測資料進(jìn)行了自動分類，分類結(jié)果全部正確。又將上述方案應(yīng)用到2008和2009年粵西海域水體分類中，通過繪制水體分類專題圖，對比分析了粵西海域2008—2010年水體類型分布的變化情況，為粵西海域的生態(tài)治理與資源規(guī)劃提供了依據(jù)。

圖1 粵西海域站位分布圖Fig.1 Distribution of monitoring stations in seawaters off western Guangdong

1 研究區(qū)及水質(zhì)資料概況

粵西海域重點研究范圍為20°20′～21°40′N，110°10′～111°00′E。沿海地區(qū)包括茂名和湛江兩市，是廣東海洋經(jīng)濟(jì)圈的重要組成部分。該海域水質(zhì)采樣站位共計45個（圖1），其中國定站位8個，省定站位37個。

水質(zhì)資料包括全氮（TN）、全磷（TP）、溶解氧（DO）、化學(xué)需氧量（COD）、酸堿度（pH）、石油類（Oils）共計6個水質(zhì)指標(biāo)，時間為2008—2010年每年8月。我們選用2010年粵西海域水質(zhì)資料為研究對象，選取其中平面分布均勻的35個站位水質(zhì)資料作為樣本，其余10個站位水質(zhì)資料作為預(yù)測（圖1）。在2010年水體分類方案應(yīng)用的基礎(chǔ)上，又利用該方案得出了2008和2009年粵西海域水體類型的分布情況。

2 水質(zhì)指標(biāo)關(guān)系

海洋水質(zhì)評價是1種多指標(biāo)的復(fù)雜評價問題，因此明確影響水體分類的主導(dǎo)指標(biāo)，對于后續(xù)的水質(zhì)評價研究尤為重要。主成分分析是研究如何通過少數(shù)幾個主成分來解釋多個變量的方差—協(xié)方差結(jié)構(gòu)的1種數(shù)學(xué)變換方法，它能夠?qū)⒍鄠€實測變量轉(zhuǎn)換為少數(shù)幾個不相關(guān)的綜合指標(biāo)，并描述原始實測變量與后期綜合指標(biāo)間的關(guān)系。由此我們可以通過主成分分析方法來描述海洋水質(zhì)各指標(biāo)間的權(quán)重關(guān)系，為后續(xù)的水體分類研究奠定基礎(chǔ)。

將粵西海域2010年45個采樣站位的6個水質(zhì)指標(biāo)（TN，TP，DO，COD，pH和Oils）利用主成分分析方法進(jìn)行處理。第1個主成分的特征根為3.05，解釋6個原始變量總方差的50.831%；第2個主成分的特征根為1.078，解釋6個原始變量總方差的17.968%；第3個主成分的特征根為0.874，解釋6個原始變量總方差的14.571%。截止到第3個主成分，累計方差貢獻(xiàn)率達(dá)到83.37%，說明前3個主成分已經(jīng)能夠很好的解釋原來6個水質(zhì)指標(biāo)的信息了。從主成分?jǐn)?shù)與特征值所繪制的碎石圖（圖2）也可以看出，第1，2主成分的特征值連線較之第2，3主成分特征值連線的斜率有較大差別，而第2，3主成分特征值連線與其后的各主成分特征值連線斜率基本一致。

利用回歸法計算得出各水質(zhì)指標(biāo)與各主成分的相關(guān)系數(shù)（表1）。其中第一主成分中TN的權(quán)重最高，繼而為TP。說明粵西海域2010年水體分類著重考慮的水質(zhì)指標(biāo)應(yīng)為TN和TP。

圖2 2010年主成分分析碎石圖Fig.2 Scree plot of principal component analysis in 2010

表1 2010年粵西海域水質(zhì)指標(biāo)與主成分關(guān)系Table 1 Relationship between the principal components and the water quality indexes in 2010in the sea area off western Guangdong

3 水體分類樣本的建立

擬采用判別分析方法建立粵西海域水質(zhì)自動分類模型，判別分析是對已知分類的數(shù)據(jù)建立判別函數(shù)，然后把這樣的規(guī)則應(yīng)用到未知分類的數(shù)據(jù)中去實現(xiàn)其分類。因此首先應(yīng)建立粵西海域已知分類的水質(zhì)資料樣本。建立水體分類樣本的步驟如下：

1）建立粵西海域水質(zhì)指標(biāo)標(biāo)準(zhǔn)

按照GB3097－1997《國家海水水質(zhì)標(biāo)準(zhǔn)》，將粵西海域水體分為4類（表2）。

表2 粵西海域水質(zhì)指標(biāo)標(biāo)準(zhǔn)Table 2 National standard for the sea water quality

2）依據(jù)水質(zhì)指標(biāo)權(quán)重關(guān)系建立樣本

針對2010年粵西海域35個站位水質(zhì)資料，按照粵西海域水質(zhì)指標(biāo)分類標(biāo)準(zhǔn)，重點考慮TN和TP質(zhì)量濃度的分類界限，建立粵西海域2010年由6個水質(zhì)指標(biāo)（TN，TP，DO，COD，pH和Oils）組成的對應(yīng)4種水體類型的水體分類樣本。其中第一類水體16個樣本點，第二類水體2個樣本點，第三類水體5個樣本點，第四類水體12個樣本點。

4 水質(zhì)自動分類方法

4.1 方法原理

判別分析在解決多元分析問題中應(yīng)用較廣，其基本思想是在已知若干個體的幾組屬性（或特征）觀測數(shù)據(jù)及個體所屬類別的前提下，根據(jù)某些準(zhǔn)則建立判別公式，然后對未知類型的個體進(jìn)行判別分析。判別分析法按判別準(zhǔn)則的不同可分為距離判別法、Fisher判別法、Bayes判別法和逐步判別法等。

設(shè)有n項水質(zhì)指標(biāo)，建立判別函數(shù)：

式中，Yi為第i個判別函數(shù)；a1，a2，…，an為各水質(zhì)指標(biāo)的系數(shù)；x1，x2，…，xn為樣本的各水質(zhì)指標(biāo)。本次研究的海水水質(zhì)指標(biāo)只有6項，則應(yīng)將所有指標(biāo)都進(jìn)行判別分析，故選取了應(yīng)用較為廣泛的Fisher判別法。

4.2 實際應(yīng)用

利用Fisher判別法對2010年粵西海域水質(zhì)樣本數(shù)據(jù)進(jìn)行了分析，從中共可獲得3個典型判別函數(shù)（表3），其特征值為該函數(shù)組間平方和與組內(nèi)平方和的比值。其中第1個判別函數(shù)對于分類的貢獻(xiàn)率最大，達(dá)到90.1%，也就是說利用第1個判別函數(shù)就已經(jīng)能夠較好的區(qū)分4類水體了。

表3 判別函數(shù)說明表Table 3 Necessary explanation for discriminant functions

用x1，x2，…，x6表示TN，TP，DO，COD，pH，Oils這6個水質(zhì)指標(biāo)，得到3個判別函數(shù)表達(dá)式：

由于前2個判別函數(shù)的貢獻(xiàn)率已經(jīng)達(dá)到98.6%，故選取前2個判別函數(shù)進(jìn)行分類。將全部35個水質(zhì)數(shù)據(jù)的第一判別函數(shù)值與第二判別函數(shù)值在平面上投點，再根據(jù)樣本中各點與類中心位置距離遠(yuǎn)近劃分各點歸屬哪類（圖3）。由Fisher判別法分析得出的2010年粵西海域水質(zhì)樣本分類問題的回判正確率情況：第一類水體共計16個樣本點，有1個樣本點分類錯誤，對于一類水體判斷符合率為93.8%；第二類水體共計2個樣本點，全部分類正確，對于二類水體判斷符合率為100%；第三類水體共計5個樣本點，全部分類正確，對于三類水體判斷符合率為100%；第四類水體共計12個樣本點，有2個樣本點分類錯誤，對于四類水體判斷符合率為83.3%（表4）。全部35個樣本點判斷符合率為91.4%。說明利用Fisher判別法所建立的判別函數(shù)，能夠有效的解決粵西海域水質(zhì)自動分類問題。為了進(jìn)一步驗證上述方案的有效性，將粵西海域其余10個站位的水質(zhì)數(shù)據(jù)代入式（2）和（3）中，計算F1和F2的值，投點判斷其歸屬哪類（圖4），其中5，6，7號站位為第一類水體，判斷正確；8號站位屬第二類水體，判斷正確；1，2號站位屬第三類水體，判斷正確；3，4，9，10號站位屬第四類水體，判斷正確，則全部10個站位的水質(zhì)預(yù)測分類正確率為100%。

表4 逐步判別分析正確率Table 5 Accuracy of stepwise discriminant analysis

在粵西海域2010年水質(zhì)資料的分析中該套自動分類方案同樣適用于研究不同年限同一采樣站位群所獲得的水質(zhì)資料，因此，按照先明確水體分類主導(dǎo)指標(biāo)，再依照海水水質(zhì)標(biāo)準(zhǔn)建立樣本，繼而利用判別分析方法建立判別函數(shù)進(jìn)行水質(zhì)自動分類的思路，對粵西海域2008—2010年45個站位的水質(zhì)資料進(jìn)行了自動分類，并繪制了水體分類專題圖（圖5）。

圖5 2008－2010年水體分類圖Fig.5 Water body classification from 2008to 2010

2008—2010年粵西海域水體分類圖表明：2008年時粵西海域大部分水體屬第一類水體，僅在特呈島附近海域出現(xiàn)了第二類水體和及小范圍的第三類水體；2009年時，粵西海域水質(zhì)狀況急劇變差，近岸海域出現(xiàn)了大面積的第三類及第四類水體；2010年時，粵西海域的水質(zhì)狀況出現(xiàn)了好轉(zhuǎn)，第二、三、四類水體的面積均變小，但是近岸仍然存在大面積的第三、四類水體，嚴(yán)重影響了海洋水產(chǎn)養(yǎng)殖業(yè)的發(fā)展。

5 結(jié) 論

依據(jù)粵西海域2008—2010年水體實測情況，建立了適用于該區(qū)域的水體分類方案，利用該分類方案得出了該海域水體分布的具體情況：

1）粵西海域的水體分類方案：先通過主成分分析確定水體分類主導(dǎo)指標(biāo)，再通過判別分析方法建立關(guān)于主導(dǎo)指標(biāo)的2個判別函數(shù)，通過各類水體兩函數(shù)交會區(qū)域范圍的限定實現(xiàn)水體自動分類。該模型對粵西海域2010年35個站位的回判符合率達(dá)到91.4%，10個預(yù)測站位的正確率達(dá)到100%。

2）通過上述水質(zhì)分類方案得出了粵西海域這3年中各類水體的分布狀況，可概括為2008年整體水質(zhì)狀況較好，2009年水質(zhì)狀況急劇變差，特別是近岸海域，到2010年時水質(zhì)狀況稍有好轉(zhuǎn)，但近岸的海水水質(zhì)狀況仍不利于海洋水產(chǎn)養(yǎng)殖業(yè)的發(fā)展。

（References）：

［1］JACOBS H L.Water quality criteria［J］.Joumal of Water Pollution Control Federation，1965，37（5）：292－300.

［2］LIU C M，WEN W Y，HE Q X.Actuality survey and assessment of marine environment in Qiongzhou Strait［J］.Marine Environmental Science，2002，21（4）：24－28.劉翠梅，溫偉英，何清溪.瓊州海峽海洋環(huán)境現(xiàn)狀調(diào)查與評價［J］.海洋環(huán)境科學(xué)，2002，21（4）：24－28.

［3］ZHANG Q F，ZHANG J L，DING J K.Assessment of marine environmental quality in Yanghe River Estuary［J］.Marine Environmental Science，1997，16（4）：38－42.張秋豐，張建樂，丁家坤.洋河口海域海洋環(huán)境質(zhì)量現(xiàn)狀評價［J］.海洋環(huán)境科學(xué)，1997，16（4）：38－42.

［4］SONG L，ZHOU Z C，WANG N B，et al.An appraisal of the environmental quality status and trend in Liaoning Coastal Area［J］.Fisheries Science，2007，26（11）：613－618.宋倫，周遵春，王年斌，等.遼寧省近岸海洋環(huán)境質(zhì)量狀況與趨勢評價［J］.水產(chǎn)科學(xué)，2007，26（11）：613－618.

［5］WALSKI T M，PARKER F L.Consumers water quality index［J］.Journal of Environmental Engineering Division，1974，（EE3）：593－601.

［6］SHANNON E E，BROEZONIK P L.Eutrophication analysis：Amultivariate approach［J］.Journal of the Sanitary Engineering Division，1972，98（SA1）：37－57.

［7］FU H，SUN Y L，SUN L，et al.The application of grey relational analysis in assessment of marine environmental quality［J］.Transactions of Oceanology and Limnology，2007，5（3）：127－131.付會，孫英蘭，孫磊，等.灰色關(guān)聯(lián)分析法在海洋環(huán)境質(zhì)量評價中的應(yīng)用［J］.海洋湖沼通報，2007，5（3）：127－131.

［8］QIN C B，ZHENG B H，QIN Y W，et al.Gray correlation assessment of water environment quality for Tianjin coastal zone in Bohai Bay［J］.Research of Environmental Sciences，2006，19（6）：94－99.秦昌波，鄭丙輝，秦延文，等.渤海灣天津段海岸帶水環(huán)境質(zhì)量灰色關(guān)聯(lián)度評價［J］.環(huán)境科學(xué)研究，2006，19（6）：94－99.

［9］LIU J Y，YANG T X，LI M，et al.A weight absolute gray correlation degree and it’s application in evaluation of water quality of Miyun Reservoir［J］.Journal of Jilin University：Earth Science Edition，2005，35（1）：54－58.劉金英，楊天行，李明，等.一種加權(quán)絕對灰色關(guān)聯(lián)度及其在密云水庫水質(zhì)評價中的應(yīng)用［J］.吉林大學(xué)學(xué)報：地球科學(xué)版，2005，35（1）：54－58.

［10］WANG B X，YANG L F.Bioassessment of Qianhuai River using a river biological index［J］.Acta Ecologica Sinica，2003，23（10）：2082－2091.王備新，楊蓮芳.用河流生物指數(shù)評價秦淮河上游水質(zhì)的研究［J］.生態(tài)學(xué)報，2003，23（10）：2082－2091.

［11］CAI R Y.Artificial neutral network predicting model of noctiluca scientillans density and physical and chemical parameters［J］.Environmental Monitoring in China，2001，17（3）：52－55.蔡如鈺.基于人工神經(jīng)網(wǎng)絡(luò)的夜光藻密度預(yù)測模型［J］.中國環(huán)境監(jiān)測，2001，17（3）：52－55.

Research and Application of the Method for Water Quality Evaluation in the Sea Area off Western Guangdong

ZHANG Ying1，XIE Shi－yi1，YANG Feng2
（1.LabofOceanRemoteSensingandInformationTechnology，GuangDongOceanUniversity，Zhanjiang 524088，China；2.ZhanjiangOceanicandFisheryEnvironmentalMonitoringStation，Zhanjiang 524039，China）

Aimed at the existing problems of seawater quality evaluation，aplan has been established，in which first is to determine the weight relationships among various water quality indexes by using principal component analysis，then to classify the samples according to the weight relationships and the water quality classification criteria，and last to establish a discriminant function by discriminant analysis and to realize the water body classification automatically.Good results have been obtained from the application of the plan for the seawaters off western Guangdong.The returned discriminant result of the water body classification reaches to a high compliance rate of 91.4%.Based on the automatic classification，thematic maps of water body distributions during the period from 2008to 2010have also been plotted for the seawaters off western Guangdong，which could provide a decision basis for the pollution control in this sea area.

sea area off western Guangdong；water quality classification；principal component analysis；discriminant analysis

February 14，2011

X824

A

1671－6647（2012）02－0198－07

2011－02－14

廣東省科技項目——粵西海域衛(wèi)星遙感水質(zhì)監(jiān)測速報與預(yù)警系統(tǒng)（2009B030600007）

張 瑩（1982－），女，吉林四平人，講師，博士，主要從事海洋探測信息處理與解釋方面研究.

E－mail：zhangying198206＠126.com

＊通訊作者，男，四川巴中人，教授，碩士.E－mail：shiyxie＠126.com