基于模糊聚類的PM2.5擬合組分選擇模型的研究

徐恒鵬,李 岳,史國(guó)良,王 瑋*,軒淑艷(1.南開(kāi)大學(xué)計(jì)算機(jī)與控制工程學(xué)院,天津 00071；2.南開(kāi)大學(xué)軟件學(xué)院,天津 00071；.南開(kāi)大學(xué)環(huán)境科學(xué)與工程學(xué)院,國(guó)家環(huán)境保護(hù)城市空氣顆粒物污染防治重點(diǎn)實(shí)驗(yàn)室,天津 00071；.河北省唐山市玉田縣環(huán)境保護(hù)局,河北 唐山 06199)

?

基于模糊聚類的PM2.5擬合組分選擇模型的研究

徐恒鵬1,2,李 岳1,2,史國(guó)良3,王 瑋1,2*,軒淑艷4(1.南開(kāi)大學(xué)計(jì)算機(jī)與控制工程學(xué)院,天津 300071；2.南開(kāi)大學(xué)軟件學(xué)院,天津 300071；3.南開(kāi)大學(xué)環(huán)境科學(xué)與工程學(xué)院,國(guó)家環(huán)境保護(hù)城市空氣顆粒物污染防治重點(diǎn)實(shí)驗(yàn)室,天津 300071；4.河北省唐山市玉田縣環(huán)境保護(hù)局,河北 唐山 064199)

摘要：提出了一種新的P M2.5源成分譜擬合組分選擇模型,在充分考慮擬合過(guò)程的物理意義的基礎(chǔ)上,采用聚類正確率作為組分選擇的依據(jù).實(shí)驗(yàn)驗(yàn)證,該模型能夠準(zhǔn)確獲取較好的擬合主組分, 相比與經(jīng)驗(yàn)選或者手動(dòng)盲選所得擬合結(jié)果,我們提出的模型將成功擬合(誤差范圍在0～0.05之間)的比例由40%提升到83%.

關(guān)鍵詞：PM2.5源成分譜；組分選擇；CMB受體模型；源解析；模糊聚類

* 責(zé)任作者, 副教授, kevinwangwei@nankai.edu.cn

近年來(lái),多地頻發(fā)的霧霾天氣,讓PM2.5成為時(shí)下最為熱門的詞匯之一,國(guó)家出臺(tái)了《國(guó)十條》以指導(dǎo)當(dāng)前大氣污染防治[1].在進(jìn)行大氣污染治理之前,必須明確各種大氣污染源類的分擔(dān)率.大氣顆粒物源解析的方法主要分為兩大類:擴(kuò)散模型和受體模型,其中發(fā)展最為迅速和成熟的是受體模型.自20世紀(jì)70年代以來(lái),提出了化學(xué)質(zhì)量平衡法(CMB)、因子分析法(FA)等受體模型,其中化學(xué)質(zhì)量平衡(CMB)模型由于其物理意義明確且能定量計(jì)算各個(gè)源類貢獻(xiàn)度而成為發(fā)展最迅速、應(yīng)用最廣泛的源解析技術(shù)[2],被我國(guó)環(huán)保部和美國(guó)EPA列為推薦模型,在中國(guó)得到了廣泛的發(fā)展和應(yīng)用[3].胡珊等[4]利用CMB模型對(duì)珠江三角洲PM2.5進(jìn)行致癌風(fēng)險(xiǎn)及源解析、馮銀廠等[5]采用CMB模型對(duì)烏魯木齊環(huán)境空氣中TSP和 PM10源解析、鄒長(zhǎng)武等[6]采用CMB模型提出了一種混合塵溯解析新方法.

雖然CMB模型在我國(guó)的顆粒物防治工作中發(fā)揮了重要作用,但是在利用CMB進(jìn)行擬合時(shí),如何選取擬合組分仍沒(méi)有得到很好的解決.靠人工或經(jīng)驗(yàn)挑選擬合組分,不僅工作量大且可能漏掉更優(yōu)結(jié)果,很多研究人員提出了不同的解決思路,如2001馮銀廠等[7]提出“窮舉法”方法,即通過(guò)增加一些診斷方法以期獲取“最優(yōu)”的擬合結(jié)果[7].

本文針對(duì)PM2.5源成分譜在CMB模型擬合過(guò)程中存在的人工或經(jīng)驗(yàn)選取組分時(shí)工作量大且較難獲取“最優(yōu)”擬合結(jié)果的問(wèn)題,提出了一種基于模糊聚類的PM2.5源成分譜擬合組分選擇模型,旨在為PM2.5源成分譜在CMB模型擬合過(guò)程中的組分選擇問(wèn)題提供一種新的解決思路,并對(duì)本文提出的組分選擇方法進(jìn)行了可行性研究.

1 基于模糊聚類的PM2.5源成分譜參與擬合組分的選擇模型描述

PM2.5源成分譜數(shù)據(jù)是典型的小樣本、高維度數(shù)據(jù),且不同維度含量的數(shù)量級(jí)有時(shí)相差極大,出現(xiàn)維數(shù)災(zāi)難[8],即隨著維數(shù)的增加,計(jì)算量呈指數(shù)倍增長(zhǎng).有鑒于此,可將PM2.5源成分譜擬合組分問(wèn)題轉(zhuǎn)換為機(jī)器學(xué)習(xí)中的特征選擇問(wèn)題.通過(guò)特征選擇方法,對(duì)源成分譜數(shù)據(jù)進(jìn)行處理,將得到的相應(yīng) “主組分”當(dāng)作參與擬合的組分.本文采用模糊聚類方法對(duì)源成分譜數(shù)據(jù)進(jìn)行聚類,將聚類正確率作為PM2.5源成分譜擬合組分選擇依據(jù)或標(biāo)準(zhǔn).

圖1 算法整體流程Fig.1 The flow chart of algorithm

本文提出的基于模糊聚類的PM2.5源成分譜的參與擬合組分選擇模型的實(shí)驗(yàn)流程如圖1所示.輸入為PM2.5源成分譜數(shù)據(jù),首先設(shè)定組分變異、組分含量橫向?qū)Ρ纫约敖M分含量縱向?qū)Ρ冗@3組參數(shù)的取值,進(jìn)而生成候選組分集集合,對(duì)所有組分候選集進(jìn)行聚類中心初始化,再用FSC聚類算法進(jìn)行模糊聚類,最后計(jì)算聚類正確率,將聚類正確率最高的參數(shù)對(duì)應(yīng)的候選組分集作為參與擬合的最佳組分集.在進(jìn)行聚類中心初始化時(shí)采用協(xié)方差矩陣進(jìn)行樣本篩選及聚類中心初始化操作,在后面章節(jié)2.1、2.2、2.3、2.4將詳細(xì)介紹.

1.1 組分候選集生成方法

PM2.5源成分譜中,樣本數(shù)目少,組分?jǐn)?shù)目較多,根據(jù)組分的物理意義,設(shè)定3組參數(shù)來(lái)生成擬合組分候選集.

組分變異參數(shù):對(duì)于每種源類,其中某組分變異系數(shù)越小,說(shuō)明該組分離散度越小.將該組分納入CMB等主流源解析模型中,擬合結(jié)果更準(zhǔn)確.組分變異系數(shù)其中,cvi,l表示對(duì)于源類i中組分l的變異系數(shù),σi,l表示源類i中組分l的標(biāo)準(zhǔn)差,ui,l表示源類i中組分l的均值.

組分含量橫向?qū)Ρ葏?shù):即組分含量大小.對(duì)于每種源類,組分含量較大的若干個(gè)組分,一般均為該源類的必選組分,需納入CMB等主流源解析模型,也可稱為組分含量橫向?qū)Ρ认禂?shù).

組分含量縱向?qū)Ρ葏?shù):對(duì)于某種源類,有一些標(biāo)識(shí)性組分是其他源類所不具有的,組分含量縱向?qū)Ρ认禂?shù)用來(lái)找出該種源類的標(biāo)識(shí)性組分.對(duì)于源類i,假設(shè)組分l是其標(biāo)識(shí)性組分,那么para_veri,l值就為1(測(cè)量精確的情況).對(duì)比系數(shù)如下:

式中,para_veri,l表示源類i中組分l的組分縱向?qū)Ρ认禂?shù);ui,l表示源類i中組分l的均值;n表示源類個(gè)數(shù).

設(shè)定上述3個(gè)參數(shù)取值范圍,組合生成不同的候選組分集η.例如,可以將求出的變異系數(shù)劃分為10個(gè)取值空間,即10個(gè)候選集.對(duì)于組分含量橫向?qū)Ρ燃翱v向?qū)Ρ葏?shù)也可以通過(guò)設(shè)定范圍來(lái)生成各自的候選組分集.最終,將上述生成的三類組分候選集進(jìn)行簡(jiǎn)單組合去重后,生成最終的候選組分集η.

1.2 FSC模糊加權(quán)軟子空間聚類算法

過(guò)去幾十年里,針對(duì)已有的傳統(tǒng)聚類方法,出現(xiàn)了很多基于特征加權(quán)或特征選擇的數(shù)據(jù)挖掘方法,其中子空間聚類是目前高維數(shù)據(jù)聚類分析中一個(gè)非常重要的研究領(lǐng)域[9-10].根據(jù)聚類方法的不同,主要分為兩種聚類形式[12-13]:硬子空間聚類、軟子空間聚類.

軟子空間聚類算法具有更好的適應(yīng)性與靈活性,可以分為模糊加權(quán)子空間聚類方法及熵加權(quán)子空間聚類算法.模糊加權(quán)子空間聚類算法包括FWKM、FSC[11]、AWA,該類方法對(duì)每一類的各組分進(jìn)行模糊加權(quán),對(duì)于每一個(gè)組分,樣本分布越緊湊,則組分獲得權(quán)值越大,與CMB擬合算法思想一致.此外,組分方差越小,所占權(quán)重越大,同時(shí)為權(quán)重賦予一個(gè)模糊指數(shù),使該聚類算法更具有適應(yīng)性.

在FSC算法中,一個(gè)源類在某個(gè)組分或特征上的密度越大,該組分或特征對(duì)于該源類的貢獻(xiàn)度或權(quán)值越大.對(duì)于給定的源成分譜樣本數(shù)據(jù)集,本文利用軟子空間聚類算法得到k個(gè)源類中心值,即每個(gè)源類的信息:.定義uji表示第i個(gè)樣本xi屬于第j個(gè)聚類中心zj的模糊隸屬度,定義U表示整個(gè)源類的模糊隸屬度矩陣,.

此外,以期更好發(fā)現(xiàn)各個(gè)數(shù)據(jù)簇相應(yīng)的子空間結(jié)構(gòu),軟子空間聚類算法在聚類過(guò)程中對(duì)每個(gè)數(shù)據(jù)簇的全部特征賦予一個(gè)特征加權(quán)系數(shù).因此,定義wjh表示第h個(gè)特征對(duì)于第j個(gè)數(shù)據(jù)簇的重要性或貢獻(xiàn)度,則W表示整個(gè)源分譜樣本數(shù)據(jù)集的特征加權(quán)系數(shù)矩陣,.

FSC模糊聚類的目標(biāo)函數(shù)如下[11]:

文獻(xiàn)[11]給出了上述公式的詳細(xì)算法步驟及W、Z和U的迭代公式.

在初始化聚類中心時(shí)Z時(shí),采用了基于協(xié)方差矩陣的樣本篩選方法[14].

1.3 基于模糊聚類的PM2.5源成分譜組分選擇模型描述

本文采用聚類正確率作為組分選擇衡量標(biāo)準(zhǔn).聚類正確率是指被正確劃分樣本占總體樣本的比率,定義為:

其中,numl表示源類l被正確劃分的樣本數(shù),N表示所有源成分譜樣本數(shù)目.

對(duì)于組分候選集中每一個(gè)組分集合ηt,利用文獻(xiàn)[11]樣本篩選方法進(jìn)行聚類中心初始化,然后根據(jù)算法FSC進(jìn)行模糊聚類,最后計(jì)算該ηt下的聚類正確率CAt.最后選取CA數(shù)組中最大的值對(duì)應(yīng)的組分候選集作為輸出.

測(cè)繪新技術(shù)在建筑工程測(cè)量中的地位越來(lái)越重要，尤其是隨著建筑工程行業(yè)的規(guī)模不斷壯大其重要性越來(lái)越突出。但與此同時(shí)，傳統(tǒng)的測(cè)繪技術(shù)因?yàn)樾省?zhǔn)確性和可靠性的問(wèn)題很難再滿足建筑工程發(fā)展的需要。在這種狀況下，越來(lái)越多的新的測(cè)繪技術(shù)不斷涌現(xiàn)對(duì)推動(dòng)測(cè)繪事業(yè)的發(fā)展以及提高建筑工程測(cè)量的準(zhǔn)確度和可靠性起到了重要作用。

算法描述如下:

算法.Dim_Sel algorithm

輸入:PM2.5源成分譜數(shù)據(jù)集,源類個(gè)數(shù)k,組分候選集η,模糊加權(quán)系數(shù)α,一個(gè)無(wú)窮小參數(shù)ε.

輸出:輸出最佳參與擬合組分

第1步:從候選集η選取一個(gè)組分集ηt;

第2步:按照文獻(xiàn)[14]樣本篩選算法進(jìn)行聚類中心初始化;

第3步:按照文獻(xiàn)[11]模糊聚類算法聚類; 第4步:依據(jù)公式(3)計(jì)算聚類正確率CA ;

第5步:計(jì)算出每個(gè)閾值對(duì)應(yīng)的CA,找到CA最大值對(duì)應(yīng)的組分候選集輸出.

2 方法驗(yàn)證與評(píng)估

2.1 受體生成方法

為了保證實(shí)驗(yàn)結(jié)果的準(zhǔn)確性、有效性,實(shí)驗(yàn)所用的PM2.5源成分譜數(shù)據(jù)是來(lái)自于實(shí)際監(jiān)測(cè)數(shù)據(jù),該數(shù)據(jù)集中,每個(gè)樣本包括101個(gè)屬性,含有3個(gè)源類:SOIL類、COAL類、VEHICLE類,即在進(jìn)行模糊聚類時(shí)源類個(gè)數(shù)k設(shè)定為3.實(shí)驗(yàn)所用源解析工具是美國(guó)環(huán)保EPA-CMB8.2[15]受體模型軟件.受體生成方法參考相關(guān)文獻(xiàn),從PM2.5源成分譜3種源類中隨機(jī)選取3個(gè)樣本,再按照隨機(jī)設(shè)定的貢獻(xiàn)度生成一條受體.實(shí)驗(yàn)時(shí),隨機(jī)生成30條受體數(shù)據(jù),對(duì)所有組分候選集進(jìn)行驗(yàn)證,進(jìn)而證明本文提出的組分選擇的方法可行性.

2.2 擬合結(jié)果衡量方法

CMB軟件擬合得出的各個(gè)源類的貢獻(xiàn)度需與真實(shí)貢獻(xiàn)度進(jìn)行比較,本文采用相似度系數(shù)及平均絕對(duì)誤差A(yù)AE[16]進(jìn)行衡量.相似度系數(shù)包括夾角余弦和相關(guān)系數(shù).

夾角余弦函數(shù)忽略兩數(shù)據(jù)點(diǎn)(向量)之間的絕對(duì)長(zhǎng)度而考慮其在方向上的相互關(guān)系,擬合貢獻(xiàn)度與真實(shí)貢獻(xiàn)度越相近,其值越大.

相關(guān)系數(shù)是關(guān)于向量標(biāo)準(zhǔn)差的夾角余弦,它表示兩個(gè)向量線性相關(guān)的程度,若兩個(gè)向量越相近,其值越大.

通過(guò)計(jì)算擬合貢獻(xiàn)值和真實(shí)貢獻(xiàn)度的平均絕對(duì)誤差(AAE)來(lái)表示擬合值和真實(shí)值的平均差異,若AAE的值較低,表明擬合值和真實(shí)值較接近.本文通過(guò)上述3種衡量標(biāo)準(zhǔn),對(duì)算法Dim_sel驗(yàn)證.

2.3 聚類正確率與組分選擇的關(guān)系驗(yàn)證

圖2為閾值集η與聚類正確率CA關(guān)系,其中,橫坐標(biāo)表示閾值集η,縱坐標(biāo)表示聚類正確率.從圖2可知第11組組分候選集相對(duì)應(yīng)的聚類正確率最高,約為84%.

圖2 聚類正確率變化Fig.2 Variation of CA

圖3、圖4、圖5中分別為采用相似系數(shù)、夾角余弦及平均絕對(duì)誤差對(duì)擬合結(jié)果進(jìn)行評(píng)估的折線圖,橫坐標(biāo)代表組分候選集η,縱坐標(biāo)分別為30條受體擬合結(jié)果與真實(shí)貢獻(xiàn)度的相似系數(shù)、夾角余弦及平均絕對(duì)誤差的平均值大小.圖3與圖4的走勢(shì)和圖2的走勢(shì)近乎一致,圖5的走勢(shì)和圖2的走勢(shì)幾乎相反,表明選擇聚類正確率較高的點(diǎn)對(duì)應(yīng)的組分候選集進(jìn)行CMB擬合,其擬合正確率較高,說(shuō)明了將聚類正確率作為擬合組分選擇標(biāo)準(zhǔn)的有效性.

圖3 相似系數(shù)衡量擬合結(jié)果變化Fig.3 Evaluation of fitting results by similarity factor

圖4 夾角余弦衡量擬合擬合結(jié)果變化Fig.4 Evaluation of fitting results by angle cosine

圖5 平均絕對(duì)誤差衡量擬合結(jié)果變化Fig.5 Evaluation of fitting results by AAE

2.4 參與擬合最佳組分與盲選法對(duì)比試驗(yàn)

圖6為候選集組分個(gè)數(shù)與聚類正確率對(duì)應(yīng)關(guān)系圖,從圖中可知聚類正確率最高點(diǎn)對(duì)應(yīng)的參與擬合組分是16個(gè),其中包含了Al、Si、OC、EC等3種源類的標(biāo)識(shí)性組分,也是擬合過(guò)程中人工經(jīng)驗(yàn)篩選的必選組分,也包含了Ca、K、Cr、Fe等主要組分[2,5,7].

圖6 候選集組分個(gè)數(shù)與聚類正確率Fig.6 The number of components corresponding to CA

本文設(shè)計(jì)了兩組實(shí)驗(yàn)進(jìn)行對(duì)比:對(duì)于實(shí)驗(yàn)隨機(jī)生成的30條受體,一組選取聚類正確率最高點(diǎn)對(duì)應(yīng)的16個(gè)組分進(jìn)行擬合;另一組采用k折交叉驗(yàn)證思想,對(duì)20次手動(dòng)盲選選取組分進(jìn)行擬合,其中選取組分中均含三種源類常用的擬合組分Al,Si,OC,EC等.手動(dòng)盲選實(shí)驗(yàn)中即包含隨機(jī)選取組分,又包含三種源類的常用組分,使對(duì)比實(shí)驗(yàn)結(jié)果更加客觀與穩(wěn)定.

圖7為利用AAE對(duì)20次手動(dòng)盲選實(shí)驗(yàn)的擬合結(jié)果進(jìn)行衡量,橫坐標(biāo)表示受體ID,縱坐標(biāo)表示AAE.三條線分別對(duì)應(yīng)20次手動(dòng)盲選實(shí)驗(yàn)中每個(gè)受體擬合結(jié)果的AAE的平均值,最大值及最小值.從圖7可知:最大值線和最小值線之間的浮動(dòng)較大,說(shuō)明采用手動(dòng)盲選組分選擇方法得到的擬合結(jié)果較為隨機(jī)且常出現(xiàn)無(wú)法得到擬合結(jié)果的情況;平均值線與最小值線浮動(dòng)較小,說(shuō)明平均值線可較好體現(xiàn)手動(dòng)盲選實(shí)驗(yàn)的擬合結(jié)果.

圖8為利用AAE對(duì)盲選及算法選組分?jǐn)M合所得到的結(jié)果進(jìn)行衡量對(duì)比折線圖.橫坐標(biāo)代表受體ID,縱坐標(biāo)是AAE,即30條受體擬合結(jié)果與真實(shí)貢獻(xiàn)度之間的平均絕對(duì)誤差,其值越小表示擬合結(jié)果越精確.采用聚類正確率最高的16個(gè)組分所得擬合結(jié)果的AAE均遠(yuǎn)低于20次手動(dòng)盲選擬合結(jié)果AAE的平均值,再次表明采用本文算法模型選取的組分所得擬合結(jié)果較大程度上優(yōu)于隨機(jī)選取組分所得擬合結(jié)果,并且一定程度上減少了擬合實(shí)驗(yàn)的次數(shù),同時(shí)減少了PM2.5源解析過(guò)程中的工作量.

圖7 AEE衡量盲選擬合結(jié)果Fig.7 Evaluation of Random selection fitting results by AEE

圖8 AEE衡量擬合結(jié)果對(duì)比Fig.8 Fitting results contrast based on the AEE

圖9 COAL源擬合結(jié)果對(duì)比Fig.9 The comparison diagram of COAL fitting results

圖9、圖10、圖11分別表示源類COAL、SOIL、VEHICLE的算法選及隨機(jī)手選的擬合結(jié)果與真實(shí)貢獻(xiàn)度的比較.橫坐標(biāo)表示受體ID,縱坐標(biāo)表示擬合結(jié)果減去真實(shí)貢獻(xiàn)度的差值的絕對(duì)值.觀察得知,手動(dòng)盲選對(duì)應(yīng)的86%以上的受體的擬合結(jié)果與真實(shí)貢獻(xiàn)度的誤差絕對(duì)值均遠(yuǎn)高于算法選所對(duì)應(yīng)的誤差絕對(duì)值,即采用隨機(jī)手選的組分較難得到精確度較高的CMB擬合結(jié)果.

圖10 SOIL源擬合誤差對(duì)比圖Fig.10 The comparison diagram of SOIL fitting results

圖11 VEHICLE源擬合結(jié)果Fig.11 The comparison diagram of VEHICLE fitting results

采用算法選組分進(jìn)行CMB擬合,所得3種源類的擬合結(jié)果絕對(duì)誤差值范圍在(0～0.05)之內(nèi)的比例達(dá)到83%,但采用手動(dòng)盲選組分進(jìn)行擬合的比例僅為40%,即采用本文提出的組分選擇模型篩選的組分參與的擬合效果與真實(shí)貢獻(xiàn)度的誤差較小,而手動(dòng)盲選的擬合結(jié)果則和真實(shí)貢獻(xiàn)度曲線誤差較大,進(jìn)一步表明了算法模型的有效性.

3 結(jié)論

本文提出了一種基于模糊聚類的PM2.5擬合組分選擇模型,通過(guò)選取聚類正確率較高的值所對(duì)應(yīng)的組分來(lái)進(jìn)行CMB擬合.實(shí)驗(yàn)結(jié)果表明,采用本文提出的PM2.5擬合組分選擇方法,有86%的受體所得擬合結(jié)果優(yōu)于手動(dòng)盲選所得擬合結(jié)果,其誤差范圍在0～0.05之間的比例由40%提升到83%,在一定程度上提高了PM2.5源成分譜CMB擬合結(jié)果的精確度.

參考文獻(xiàn)：

[1] 大氣污染防治行動(dòng)計(jì)劃 [J].中國(guó)環(huán)保產(chǎn)業(yè), 2013,(10):4-9.

[2] Shi G L, Tian Y Z, Zhang Y F, et al.Estimation of the concentrations of primary and secondary organic carbon in ambient particulate matter: Application of the CMB-Iteration method [J].Atmospheric Environment, 2011,45(32):5692-5698.

[3] 朱 坦,吳 琳,畢曉輝,等.大氣顆粒物源解析受體模型優(yōu)化技術(shù)研究 [J].中國(guó)環(huán)境科學(xué), 2010,30(7):865-870.

[4] 胡 珊,張遠(yuǎn)航,魏永杰.珠江三角洲大氣細(xì)顆粒物的致癌風(fēng)險(xiǎn)及源解析 [J].中國(guó)環(huán)境科學(xué), 2010,30(11):1202-1208.

[5] 馮銀廠,彭 林,吳建會(huì),等.烏魯木齊市環(huán)境空氣中TSP和P M10來(lái)源解析 [J].中國(guó)環(huán)境科學(xué), 2005,25(S1):30-33.

[6] 鄒長(zhǎng)武,印紅玲,劉盛余,等.大氣顆粒物混合塵溯源解析新方法[J].中國(guó)環(huán)境科學(xué), 2011,31(6):881-885.

[7] 馮銀廠.關(guān)于化學(xué)質(zhì)量平衡(CMB)受體模型應(yīng)用中若干技術(shù)問(wèn)題的研究 [D].天津:南開(kāi)大學(xué), 2002.

[8] Scott D W.Multivariate density estimation: theory, practice, and visualization [M].Wiley.com, 2009.

[9] Müller E, Günnemann S, Assent I, et al.Evaluating clustering in subspace projections of high dimensional data [J].Proceedings of the VLDB Endowment, 2009,2(1):1270-1281.

[10] Parsons L, Haque E, Liu H.Subspace clustering for high dimensional data: a review [J].ACM SIGKDD Explorations Newsletter, 2004,6(1):90-105.

[11] Gan G, Wu J.A convergence theorem for the fuzzy subspace clustering (FSC) algorithm [J].Pattern Recognition, 2008,41(6): 1939–1947.

[12] Deng Z, Choi K S, Chung F L, et al.Enhanced soft subspace clustering integrating within-cluster and between-cluster information [J].Pattern Recognition, 2010,43(3):767–781.

[13] Jing L, Ng M K, Huang J Z.An Entropy Weighting k-Means Algorithm for Subspace Clustering of High-Dimensional Sparse Data [J].IEEE Transactions on Knowledge & Data Engineering, 2007,19(8):1026-1041.

[14] 徐恒鵬.基于智能信息處理的PM2.5源解析問(wèn)題的研究 [D].天津:南開(kāi)大學(xué), 2014.

[15] Habre R, Coull B, Koutrakis P.Impact of source collinearity in simulated PM2.5data on the PMF receptor model solution [J].Atmospheric Environment, 2011,45(38):6938-6946.

[16] Christensen W F, Gunst R F.Measurement error models in chemical mass balance analysis of air quality data [J].Atmospheric Environment, 2004,38(5):733–744.

The fitting component selection model of PM2.5based on fuzzy clustering.

XU Heng-peng1,2, LI Yue1,2, SHI Guo-liang3, WANG Wei1,2*, XUAN Shu-yan4(1.College of Computer and Control Engineering, NanKai University, Tianjin 300071, China；2.College of Software, NanKai University, Tianjin 300071, China；3.State Environmental Protection Key Laboratory of Urban Ambient Air Particulate Matter Pollution Prevention and Control, College of Environmental Science and Engineering, NanKai University, Tianjin 300071, China；4.Yutian Environmental Protection Agency, Tangshan 064199, China).China Environmental Science, 2016,36(1)：12～17

Abstract：In current research, there is a lack of uniform standards for components selection in PM2.5source profile apportionment.Researchers tend to choose the component manually and empirically, leading to a subsequent poor fitting result, or even failures.Concerning on this problem, this paper has proposed an innovative component selection model of PM2.5source profiles apportionment.On the basis of the physical representative of each component, the proposed model calculates the accuracy of fuzzy clustering as the standard score for selection.The experiments prove that our model outperforms the traditional empirical models.The successful rate for fitting, measured by the fitting errors in 0 to 0.05, grows to 83% by implementing our model, in contrast to rate of 40% from the traditional selection model.

Key words：PM2.5source profile；components selection；CMB receptor model；source apportionment；fuzzy clustering

中圖分類號(hào)：X513

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-6923(2016)01-0012-06

收稿日期：2015-06-01

作者簡(jiǎn)介：徐恒鵬(1988-),男,山東臨沂人,在讀博士,主要從事PM2.5源成分譜擬合模型研究.發(fā)表論文1篇.