曹　榮　王仁雷

（華電電力科學(xué)研究院，浙江杭州　310030）

一種基于SOM神經(jīng)網(wǎng)絡(luò)的污水處理工藝方案比選方法

曹榮王仁雷

（華電電力科學(xué)研究院，浙江杭州310030）

目前越來越多的電廠開始應(yīng)用城鎮(zhèn)污水處理廠出水作為生產(chǎn)用水水源，而其出水水質(zhì)直接取決于所采用的污水處理工藝。本文嘗試將自組織神經(jīng)網(wǎng)絡(luò)模型應(yīng)用于城鎮(zhèn)污水處理工藝方案的決策，實(shí)例分析結(jié)果證明模型和方法是實(shí)用、有效的。且與常規(guī)比選方法相比，該模型編程簡(jiǎn)潔，實(shí)現(xiàn)容易，效率更高。

城鎮(zhèn)污水處理廠自組織神經(jīng)網(wǎng)絡(luò)方案比選

1　引言

圖1　自組織映射神經(jīng)網(wǎng)絡(luò)模型

表1　4種工藝方案和理想方案及淘汰方案的評(píng)價(jià)指標(biāo)方案

目前越來越多的電廠開始應(yīng)用城鎮(zhèn)污水處理廠出水作為生產(chǎn)用水水源，而其出水水質(zhì)直接取決于所采用的污水處理工藝。城鎮(zhèn)污水處理廠工藝方案比選是一個(gè)十分復(fù)雜的系統(tǒng)分析問題，涉及經(jīng)濟(jì)、技術(shù)、社會(huì)、生態(tài)環(huán)境等諸多因素，而這些因素各自的屬性不同、重要性差異較大，許多因素的評(píng)價(jià)結(jié)果難于定量化［1］。近年來國內(nèi)外許多學(xué)者對(duì)該問題進(jìn)行了研究，提出了灰色關(guān)聯(lián)度比選模型［2-3］、層次分析法［4-5］、模糊數(shù)學(xué)評(píng)價(jià)模型［6-8］、物元分析模型［9］等眾多方法，對(duì)城鎮(zhèn)污水廠工藝方案比選起到了積極的作用。但是鑒于城鎮(zhèn)污水處理工程關(guān)系重大，而現(xiàn)有方法各有其缺陷，因此在進(jìn)行工藝方案比選時(shí)，有必要采用多種方法進(jìn)行評(píng)估，綜合比較，從而得出較為客觀與合理的結(jié)論。為此，本文提出一種新的比選方法，即自組織映射（Self-Organizing Mapping，簡(jiǎn)稱SOM）神經(jīng)網(wǎng)絡(luò)模型進(jìn)行工藝方案的比選，以求豐富和完善城鎮(zhèn)污水處理工藝方案比選的模型和方法。

2　自組織映射神經(jīng)網(wǎng)絡(luò)原理

自組織映射模型（SOM）是由芬蘭學(xué)者Kohonen于1982年提出的，這種網(wǎng)絡(luò)能模擬大腦腦皮層自組織的功能，它是一種競(jìng)爭(zhēng)式學(xué)習(xí)網(wǎng)絡(luò)，能無監(jiān)督地進(jìn)行自組織學(xué)習(xí)。自組織映射神經(jīng)網(wǎng)絡(luò)模型如圖1所示，它由兩層神經(jīng)元構(gòu)成，即輸入層和競(jìng)爭(zhēng)層（輸出層），沒有隱含層。輸入層的神經(jīng)元排成一列，其個(gè)數(shù)由輸入網(wǎng)絡(luò)的向量個(gè)數(shù)而定，競(jìng)爭(zhēng)層（輸出層）的若干神經(jīng)元排成一個(gè)二維陣列。輸入層與競(jìng)爭(zhēng)層之間實(shí)行全互連接，競(jìng)爭(zhēng)層的神經(jīng)元之間還實(shí)行側(cè)抑制連接。

Kohonen認(rèn)為當(dāng)一個(gè)神經(jīng)網(wǎng)絡(luò)接收外界輸入模式時(shí)，將會(huì)分為不同的區(qū)域，各區(qū)域?qū)斎肽Ｊ骄哂胁煌捻憫?yīng)特征，同時(shí)這一過程是自動(dòng)完成的。各神經(jīng)元的連接權(quán)值具有一定的分布，最近的神經(jīng)元互相誘導(dǎo)，而較遠(yuǎn)的神經(jīng)元?jiǎng)t相互抑制，而更遠(yuǎn)一些的則具有較弱的誘導(dǎo)作用。這樣某個(gè)輸出結(jié)點(diǎn)能對(duì)某個(gè)模式作出特別的反應(yīng)以代表該模式，而且二維平面上相鄰的結(jié)點(diǎn)能對(duì)實(shí)際模式分布中相近的模式作出特別的反應(yīng)。當(dāng)某數(shù)據(jù)模式輸入時(shí)，對(duì)其某點(diǎn)給予最大的誘導(dǎo)。以指示該類模式所屬區(qū)域，而同時(shí)對(duì)該點(diǎn)周圍的點(diǎn)也給予較大的誘導(dǎo)。因此，自組織映射神經(jīng)網(wǎng)絡(luò)的輸出狀況，不但能判斷輸入模式所屬的類別并使輸出結(jié)點(diǎn)代表某一模式，還能夠得到整個(gè)數(shù)據(jù)區(qū)域的大體分布情況，即從樣本數(shù)據(jù)中抽取到所有數(shù)據(jù)分布的大體本質(zhì)特征?？傊?，自組織特征映射神經(jīng)網(wǎng)絡(luò)是一種無教師的聚類方法，它能將高維模式映射到一平面上，而保持其拓?fù)浣Y(jié)構(gòu)不變，亦即距離相近的模式點(diǎn)，其映射點(diǎn)的距離也相近。其學(xué)習(xí)過程可描述為：對(duì)于每一個(gè)網(wǎng)絡(luò)的輸入，只調(diào)整一部分權(quán)值，使權(quán)向量更接近或更偏離輸入矢量，這一調(diào)整過程，即為競(jìng)爭(zhēng)學(xué)習(xí)。隨著不斷學(xué)習(xí)，權(quán)矢量在輸入空間相互分離，形成了各自代表輸入空間的一類模式，即實(shí)現(xiàn)了聚類功能。SOM神經(jīng)網(wǎng)絡(luò)的詳細(xì)內(nèi)容見文獻(xiàn)［10-11］。下面給出這種模型的通用算法：

（1）初始化：將權(quán)值向量Wi用小的隨機(jī)值進(jìn)行初始化。設(shè)置初始學(xué)習(xí)速率η（0）。

（2）采樣：從輸入空間中選取訓(xùn)練樣本X。

（3）近似匹配：通過歐式距離最小的標(biāo)準(zhǔn)：

來選取獲勝神經(jīng)元c，從而實(shí)現(xiàn)了神經(jīng)元的競(jìng)爭(zhēng)過程。

（4）更新：對(duì)獲勝神經(jīng)元拓?fù)溧徲騈c（n）內(nèi)的興奮神經(jīng)元，以Hebb學(xué)習(xí)規(guī)則：

更新神經(jīng)元的權(quán)值向量，從而實(shí)現(xiàn)神經(jīng)元的合作和更新過程。

（5）更新學(xué)習(xí)速率η（n）及拓?fù)溧徲騈bi（n），其中一種更新方法為：

式中rb和ri分別是輸出單元b和i在映射平面的位置；σ（n）反映了鄰域的范圍。η（n）是可變學(xué)習(xí)速度，隨迭代次數(shù)的增加而減小。也就是說，隨著訓(xùn)練過程的進(jìn)行，權(quán)值的調(diào)整幅度越來越小。Nbi（n）也隨著迭代次數(shù)而收縮，當(dāng)n足夠大時(shí)，只訓(xùn)練獲勝結(jié)點(diǎn)本身。η（n）和Nbi（n）有多種不同的表達(dá)形式，在具體訓(xùn)練過程中可以根據(jù)不同的要求不同的數(shù)據(jù)分布進(jìn)行設(shè)計(jì)。

（6）判斷迭代次數(shù)n是否超過N，如果n≤N，轉(zhuǎn)到第二步，否則結(jié)束迭代過程。

3　實(shí)例研究

3.1實(shí)例概況

南方某城鎮(zhèn)污水處理廠的設(shè)計(jì)規(guī)模為5×104m3/d，工藝方案準(zhǔn)備在三溝式氧化溝、SBR法、A2/O法和AB法4種工藝中選擇，并依次記為方案1、方案2、方案3、方案4。選擇的評(píng)價(jià)指標(biāo)有工藝技術(shù)成熟度、費(fèi)用現(xiàn)值（主要為主體構(gòu)筑物和污泥處理部分的基建投資和運(yùn)行電費(fèi)現(xiàn)值之和）、占地面積、操作管理、出水穩(wěn)定性和除磷脫氮效果6項(xiàng)。除費(fèi)用現(xiàn)值和占地面積根據(jù)概算得到外，其余4項(xiàng)模糊指標(biāo)經(jīng)多位專家打分后取平均值得到。最后得到4種工藝方案和理想方案（記為方案5）及淘汰方案（記為方案6）的評(píng)價(jià)指標(biāo)值，如表1所示。數(shù)據(jù)來自文獻(xiàn)［3］。

3.2原始數(shù)據(jù)預(yù)處理

城鎮(zhèn)污水處理工藝方案決策的影響因素錯(cuò)綜復(fù)雜，包括多個(gè)評(píng)價(jià)指標(biāo)。由于各評(píng)價(jià)指標(biāo)量綱不一，范圍大小差別較大，因此有必要對(duì)評(píng)價(jià)指標(biāo)原始數(shù)據(jù)作預(yù)處理，以改善數(shù)據(jù)質(zhì)量，保證工藝方案的可比性。常見的數(shù)據(jù)預(yù)處理方法有均值中心化、對(duì)數(shù)變換、正規(guī)化變換和標(biāo)準(zhǔn)化處理等，其中標(biāo)準(zhǔn)化處理方法是應(yīng)用最廣泛和最有效的數(shù)據(jù)預(yù)處理技術(shù)。本文就選用該方法對(duì)評(píng)價(jià)指標(biāo)原始數(shù)據(jù)作預(yù)處理。變換公式為：

由層次分析法得到6項(xiàng)評(píng)價(jià)指標(biāo)的權(quán)重為［0.3269，0.3269，0.0639，0.1165，0.1019，0.0639］，數(shù)據(jù)來自文獻(xiàn)［3］。按照式（1）、式（2）對(duì)評(píng)價(jià)指標(biāo)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，再乘以相應(yīng)的權(quán)重得到樣本矩陣X（輸入SOM網(wǎng)絡(luò)時(shí)需要進(jìn)行轉(zhuǎn)置）。

3.3模型建立與方案比選

按照SOM原理，首先將上述的各個(gè)方案標(biāo)準(zhǔn)化和加權(quán)后的評(píng)價(jià)指標(biāo)矩陣（轉(zhuǎn)置后）作為輸入樣本，則輸出為各個(gè)方案的分類，然后根據(jù)分類結(jié)果進(jìn)行排序，從中選出最優(yōu)方案。

利用MATLAB6.5神經(jīng)網(wǎng)絡(luò)工具箱［12］編制了一個(gè)程序，主要過程如下：

（1）確定輸入樣本，此處輸入樣本是標(biāo)準(zhǔn)化和加權(quán)后方案的評(píng)價(jià)指標(biāo)集。

（2）利用神經(jīng)網(wǎng)絡(luò)工具箱函數(shù)newsom建立一個(gè)SOM網(wǎng)絡(luò)，先確定競(jìng)爭(zhēng)層（輸出層）的神經(jīng)元個(gè)數(shù)為兩個(gè)，即將方案分為2類（優(yōu)類與劣類）。

（3）利用train函數(shù)對(duì)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，最大訓(xùn)練次數(shù)定為500次。

（4）最后通過sim函數(shù)實(shí)現(xiàn)輸入樣本的分類。

運(yùn)行上述自編的程序，得到輸出結(jié)果為：

式（3）中括弧內(nèi)第1個(gè)數(shù)字為神經(jīng)元編號(hào)，第2個(gè)數(shù)字為方案編號(hào)。即方案1、方案3、方案6在第2個(gè)神經(jīng)元響應(yīng)，它們屬于同一類，方案2、方案4、方案5在第1個(gè)神經(jīng)元響應(yīng)，它們也屬于同一類。由于方案5是理想方案，方案6是淘汰方案，故方案2和4要優(yōu)于1和3。

由于還沒有選出最優(yōu)方案，需要利用SOM網(wǎng)絡(luò)對(duì)樣本進(jìn)行進(jìn)一步的細(xì)分。考慮將競(jìng)爭(zhēng)層（輸出層）的神經(jīng)元個(gè)數(shù)定為3個(gè)，即將方案分為3類（優(yōu)類、中類和劣類）。

運(yùn)行程序得到輸出為：

可見方案2和5為一類，方案1和4為一類，方案3和6為一類。綜合以上兩次分類所得結(jié)果可以知道方案2是最優(yōu)方案，方案3是最差方案，4個(gè)方案的排序是2＞4＞1＞3，即SBR法＞AB法＞三溝式氧化溝＞A2/O法。故本方法推薦的比選方案是SBR法，這與實(shí)際結(jié)果是一致的，同時(shí)與文獻(xiàn)［3］灰色關(guān)聯(lián)度法比選的結(jié)果也是一致的。

4　結(jié)論與建議

（1）利用自組織神經(jīng)網(wǎng)絡(luò)模型對(duì)城鎮(zhèn)污水處理工藝方案進(jìn)行比（2）建立自組織神經(jīng)網(wǎng)絡(luò)模型時(shí)，利用MATLAB軟件的神經(jīng)網(wǎng)絡(luò)工具箱進(jìn)行編程，簡(jiǎn)潔方便，且實(shí)現(xiàn)較為容易，比常規(guī)比選方法效率更高。

選，實(shí)質(zhì)是一個(gè)對(duì)候選方案與理想開發(fā)方案及淘汰方案進(jìn)行模式識(shí)別的過程，實(shí)例證明它是合理的、可行的，該方法也適用于其它類似問題的比選決策。

（3）由于城鎮(zhèn)污水處理工藝方案的比選關(guān)系重大，因此可同時(shí)考慮幾種決策方法，對(duì)它們的結(jié)果進(jìn)行綜合比較，從而得到更為客觀、合理的方案。

［1］周律，米金套，吳春旭，等.中小城鎮(zhèn)污水處理投資決策與工藝技術(shù)［M］.北京：化學(xué)工業(yè)出版社，2002.

［2］凌猛，杭世珺.城鎮(zhèn)污水處理廠工藝方案模糊決策方法的應(yīng)用［J］.給水排水，1998，24（3）：6－9.

［3］王浙明，史惠祥，蘇雨生，等.灰色關(guān)聯(lián)模型用于工程方案優(yōu)化［J］.中國給水排水，2002，18（1）：81－84.

［4］胡天覺，陳維平，曾光明，等.運(yùn)用層次分析法對(duì)株洲霞灣污水處理廠污水處理工藝方案擇優(yōu)［J］.環(huán)境工程，2000，18（1）：61－63.

［5］張建鋒，黃廷林.關(guān)中地區(qū)污水處理工藝選擇的系統(tǒng)分析［J］.環(huán)境工程，1999，17（3）：61.

［6］李如忠.多層次模糊綜合評(píng)判模型在城鎮(zhèn)污水處理中的應(yīng)用［J］.淮南工業(yè)學(xué)院學(xué)報(bào)，2000，20（4）：1－4.

［7］鄢恒珍，龔文琪，晏全香.小城鎮(zhèn)污水處理方案的模糊相似優(yōu)先比決策模型［J］.環(huán)境科學(xué)與管理，2009，08：62-65.

［8］劉慶玉，柴麗.污水處理工藝方案的熵權(quán)模糊灰色關(guān)聯(lián)優(yōu)選法［J］.沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào)，2008.02：125-127.

［9］慕金波，酒濟(jì)明.物元分析在環(huán)境工程方案評(píng)標(biāo)中的應(yīng)用［J］.重慶環(huán)境科學(xué)，1997，19（3）：28-32.

［10］Kohonen T. Self-Organizing Maps［M］.Spring.Berlin.1995.

［11］袁曾任.人工神經(jīng)元網(wǎng)絡(luò)及其應(yīng)用［M］.北京：清華大學(xué)出版社，1999.

［12］許東，吳錚.基于MATLAB6.x的系統(tǒng)分析與設(shè)計(jì)—神經(jīng)網(wǎng)絡(luò)（第二版）［M］.西安：西安電子科技大學(xué)出版社，2002.

Self-organizing neural network model is used to make decision about process selection of urban WTP. An example shows that the method is effective and valuable. Compared with conventional methods of optimization，the new model is easier to program and realize.

Urban WTP；Self-organizing neural network；Project；Decision-making

曹榮（1982—），男，碩士，畢業(yè)于浙江工業(yè)大學(xué)，主要從事電力化學(xué)及環(huán)保技術(shù)研究工作；王仁雷（1980—），男，碩士，畢業(yè)于重慶大學(xué)，主要從事電力化學(xué)及環(huán)保技術(shù)研究工作。