環(huán)境空氣顆粒物來(lái)源解析正定矩陣因子分解模型計(jì)算技術(shù)指南

1環(huán)境空氣顆粒物來(lái)源解析正定矩陣因子分解模型計(jì)算技術(shù)指南1適用范圍本標(biāo)準(zhǔn)規(guī)定了正定矩陣因子分解（PMF）模型計(jì)算的數(shù)據(jù)準(zhǔn)備、基礎(chǔ)計(jì)算、旋轉(zhuǎn)計(jì)算、顆粒物源貢獻(xiàn)計(jì)算、結(jié)果合理性判斷等技術(shù)要求。本標(biāo)準(zhǔn)適用于應(yīng)用PMF模型對(duì)環(huán)境空氣顆粒物進(jìn)行來(lái)源解析相關(guān)的計(jì)算工作。2規(guī)范性引用文件本標(biāo)準(zhǔn)引用了下列文件或其中的條款。凡是注明日期的引用文件，僅注日期的版本適用于本標(biāo)準(zhǔn)。凡是未注日期的引用文件，其最新版本（包括所有的修改單）適用于本標(biāo)準(zhǔn)。HJ630環(huán)境監(jiān)測(cè)質(zhì)量管理技術(shù)導(dǎo)則HJ1352環(huán)境空氣顆粒物來(lái)源解析基于手工監(jiān)測(cè)的受體模型法監(jiān)測(cè)數(shù)據(jù)處理與檢驗(yàn)技術(shù)規(guī)范《環(huán)境空氣顆粒物來(lái)源解析監(jiān)測(cè)技術(shù)方法指南》（監(jiān)測(cè)函〔2020〕8號(hào)）3術(shù)語(yǔ)和定義下列術(shù)語(yǔ)和定義適用于本標(biāo)準(zhǔn)。3.1環(huán)境受體ambientreceptor受到污染源排放影響的環(huán)境空氣，簡(jiǎn)稱受體。3.2正定矩陣因子分解模型positivematrixfactorizationmodel（PMF）將樣本數(shù)據(jù)的X矩陣分解為因子貢獻(xiàn)G矩陣和因子譜F矩陣，對(duì)因子譜F矩陣進(jìn)行識(shí)別，并定量計(jì)算樣本的因子貢獻(xiàn)的一種多元因子分析類模型。3.3污染源化學(xué)成分譜sourcechemicalprofile特定污染源類排放顆粒物的相對(duì)穩(wěn)定的化學(xué)組分信息（g/g），簡(jiǎn)稱源成分譜。源成分譜應(yīng)包含該類污染源的標(biāo)識(shí)組分。3.4污染源貢獻(xiàn)sourcecontribution通過(guò)受體模型計(jì)算得到的污染源對(duì)環(huán)境空氣顆粒物的貢獻(xiàn)，用貢獻(xiàn)濃度(μg/m3）或貢獻(xiàn)占比表示。3.5污染源標(biāo)識(shí)組分sourcechemicaltracer污染源化學(xué)成分譜中對(duì)該源類有指示作用的一種或多種顆粒物化學(xué)組分，又稱示蹤組分。源標(biāo)識(shí)組分是區(qū)別該源類與其他源類的重要標(biāo)識(shí)物，每種污染源有各自的標(biāo)識(shí)組分。23.6因子factorPMF模型計(jì)算提取的G矩陣和F矩陣的每列信息，其中F矩陣的每列可識(shí)別為特定污染源類，G矩陣的每列可計(jì)算得到特定污染源類的貢獻(xiàn)。3.7Q值函數(shù)functionQvalue反映PMF模型計(jì)算結(jié)果質(zhì)量的一種目標(biāo)函數(shù)。3.8殘差residual顆粒物化學(xué)組分的實(shí)際測(cè)試值與PMF模型計(jì)算值之間的差值。3.9基礎(chǔ)計(jì)算baserun在計(jì)算過(guò)程中不加入任何數(shù)學(xué)、物理意義的約束條件，只根據(jù)模型Q值函數(shù)的收斂條件（Q值最?。┑玫揭蜃迂暙I(xiàn)G矩陣和因子譜F矩陣的一種PMF模型計(jì)算方法。3.10旋轉(zhuǎn)計(jì)算rotationalrun對(duì)于基礎(chǔ)計(jì)算得到的因子貢獻(xiàn)G矩陣和因子譜F矩陣，加入一定的約束條件再次進(jìn)行計(jì)算，轉(zhuǎn)換成另一對(duì)矩陣（G*和F*使各因子譜的標(biāo)識(shí)組分更加突出，從而解決因子譜難以識(shí)別問(wèn)題的一種PMF模型計(jì)算方法。4模型原理和計(jì)算流程4.1模型原理在環(huán)境空氣顆粒物來(lái)源解析中，PMF模型計(jì)算是將受體顆粒物化學(xué)組分濃度矩陣Xn×m因子化，分解為兩個(gè)因子矩陣，G(n×p)和F(p×m)，以及一個(gè)殘差矩陣E(n×m)（見(jiàn)公式1在G矩陣和F矩陣中的數(shù)值都是正值（即非負(fù)限制）的限定下，通過(guò)定義一個(gè)目標(biāo)Q值函數(shù)（見(jiàn)公式2），選擇Q值函數(shù)值最小時(shí)的G矩陣（因子貢獻(xiàn)矩陣）和F矩陣（因子譜矩陣）計(jì)算結(jié)果，再利用污染源標(biāo)識(shí)組分將F矩陣識(shí)別為不同的源類，對(duì)G矩陣進(jìn)行多元線性回歸計(jì)算，從而得出不同源類對(duì)環(huán)境樣品的貢獻(xiàn)。式中：X(n×m)——受體顆粒物化學(xué)組分濃度X矩陣；G(n×p)——因子貢獻(xiàn)G矩陣；F(p×m)——因子譜F矩陣；E(n×m)——?dú)埐頔矩陣；n——樣品個(gè)數(shù)；m——化學(xué)成分種數(shù)；p——解析的因子（污染源）數(shù)目。2式中：Q——目標(biāo)Q值函數(shù)；Eij——第i個(gè)樣品中第j個(gè)化學(xué)成分的殘差；σij——第i個(gè)樣品中第j個(gè)化學(xué)成分的不確定度；3n——樣品個(gè)數(shù)；m——化學(xué)成分種數(shù)。4.2計(jì)算流程應(yīng)用PMF模型進(jìn)行顆粒物來(lái)源解析的計(jì)算流程包括數(shù)據(jù)處理和檢驗(yàn)、基礎(chǔ)計(jì)算及結(jié)果分析、旋轉(zhuǎn)計(jì)算及結(jié)果分析、模型計(jì)算結(jié)果合理性判斷等環(huán)節(jié)，見(jiàn)圖1。PMF模型計(jì)算需要進(jìn)行反復(fù)優(yōu)化計(jì)算，調(diào)整計(jì)算數(shù)據(jù)和參數(shù)，以得到相對(duì)較優(yōu)的計(jì)算結(jié)果。圖1PMF模型計(jì)算流程5數(shù)據(jù)準(zhǔn)備5.1數(shù)據(jù)種類和數(shù)量要求。PMF模型需要輸入環(huán)境空氣顆粒物質(zhì)量濃度以及無(wú)機(jī)元素、碳組分、水溶性無(wú)機(jī)離子等化學(xué)組分濃度數(shù)據(jù)，應(yīng)不少于100個(gè)有效樣品的數(shù)據(jù)量。顆粒物化學(xué)組分需涵蓋主要污染4源的標(biāo)識(shí)組分，常用的標(biāo)識(shí)組分包括有機(jī)碳（OC）、元素碳（EC）、SO42-、NO3-、Cl-、NH4+、Si、Al、Na、Mg、K、Ca、Mn、Fe、Cu、Zn等。為解析更多的源類，可增加多環(huán)芳烴、左旋葡聚糖、正構(gòu)烷烴、脂肪酸、甾醇類、有機(jī)酸等一次和二次有機(jī)標(biāo)識(shí)組分。5.2數(shù)據(jù)質(zhì)量要求。用于顆粒物來(lái)源解析工作的監(jiān)測(cè)數(shù)據(jù)，應(yīng)按照采用國(guó)家標(biāo)準(zhǔn)、行業(yè)標(biāo)準(zhǔn)、地方標(biāo)準(zhǔn)的優(yōu)先次序獲得；無(wú)國(guó)家標(biāo)準(zhǔn)、行業(yè)標(biāo)準(zhǔn)、地方標(biāo)準(zhǔn)的，可采用國(guó)際標(biāo)準(zhǔn)、國(guó)外標(biāo)準(zhǔn)或研究建立的方法，經(jīng)評(píng)估確認(rèn)滿足方法要求后使用，同時(shí)監(jiān)測(cè)工作需按照HJ630的要求納入質(zhì)量管理體系，開(kāi)展質(zhì)量保證和質(zhì)量控制工作并符合相關(guān)要求。5.3數(shù)據(jù)格式要求。PMF模型需要輸入兩個(gè)數(shù)據(jù)文件1）顆粒物及化學(xué)組分質(zhì)量濃度的數(shù)據(jù)文件，（2）顆粒物及化學(xué)組分質(zhì)量濃度不確定度的數(shù)據(jù)文件。兩個(gè)數(shù)據(jù)文件內(nèi)容需要對(duì)應(yīng)，格式符合模型軟件的要求。5.4數(shù)據(jù)檢驗(yàn)。對(duì)于每個(gè)顆粒物樣品的顆粒物及化學(xué)組分質(zhì)量濃度數(shù)據(jù)，按照HJ1352，進(jìn)行顆粒物化學(xué)組分質(zhì)量平衡、陰陽(yáng)離子電荷平衡等數(shù)據(jù)檢驗(yàn)。5.5不確定度計(jì)算。不確定度是顆粒物來(lái)源解析工作中采樣、實(shí)驗(yàn)分析等環(huán)節(jié)引起的誤差，直接影響輸入PMF模型中樣品顆粒物及化學(xué)組分質(zhì)量濃度數(shù)據(jù)的計(jì)算權(quán)重，根據(jù)以下公式進(jìn)行計(jì)算：Unc=MDL式中：Unc——顆粒物及化學(xué)組分的不確定度，μg/m3；MDL——顆粒物及化學(xué)組分的檢出限，μg/m3?！?EFConc)2+(0.5MDL)2式中：Unc——顆粒物及化學(xué)組分的不確定度，μg/m3；MDL——顆粒物及化學(xué)組分的檢出限，μg/m3；Conc——顆粒物及化學(xué)組分的質(zhì)量濃度，μg/m3；EF——顆粒物及化學(xué)組分質(zhì)量濃度的誤差系數(shù)（Errorfraction）。當(dāng)Conc小于等于MDL時(shí)，采用公式（3）計(jì)算；當(dāng)Conc大于MDL時(shí)，采用公式（4）計(jì)算。EF大小范圍一般為0.1～0.6，如顆粒物及化學(xué)組分質(zhì)量濃度接近檢出限或數(shù)據(jù)缺失較多時(shí)，則相應(yīng)采用較大的數(shù)值。6基礎(chǔ)計(jì)算6.1基礎(chǔ)計(jì)算運(yùn)行6.1.1數(shù)據(jù)的模型檢驗(yàn)。通過(guò)PMF模型的濃度散點(diǎn)圖、時(shí)間序列圖等，檢驗(yàn)不同化學(xué)組分之間關(guān)系、化學(xué)組分隨時(shí)間變化情況。對(duì)于異常高或者低的化學(xué)組分?jǐn)?shù)據(jù)，經(jīng)核實(shí)沒(méi)有特定污染過(guò)程（如沙塵天氣等），可予以剔除。6.1.2數(shù)據(jù)計(jì)算權(quán)重的選擇。顆粒物的關(guān)鍵標(biāo)識(shí)組分（如OC、EC、SO42-、NO3-、NH4+、Al、Si、Ca、K、Fe等化學(xué)組分）的計(jì)算權(quán)重需設(shè)為“強(qiáng)”。其他非關(guān)鍵標(biāo)識(shí)組分（如Cl-、Na、Mg、Cu、Cd、Cr等化學(xué)組分）的計(jì)算權(quán)重通常設(shè)為“強(qiáng)”，如測(cè)量不確定性較高，可設(shè)為“弱”計(jì)算權(quán)重。6.1.3因子個(gè)數(shù)的確定。因子個(gè)數(shù)需要基于對(duì)顆粒物可能來(lái)源、樣品數(shù)量、采集時(shí)間等，通過(guò)模型反復(fù)計(jì)算確定。首先根據(jù)工作區(qū)域的基本情況（如當(dāng)?shù)嘏欧旁辞鍐?、污染源普查?shù)據(jù)、工業(yè)布局、能源結(jié)構(gòu)、機(jī)動(dòng)車保有量等信息初步確定因子個(gè)數(shù)的范圍。然后通過(guò)設(shè)置不同的因子個(gè)數(shù)嘗試計(jì)算，對(duì)各5計(jì)算情況進(jìn)行分析，結(jié)合工作區(qū)域基本情況，驗(yàn)證和確定適合的因子個(gè)數(shù)，常用方法如下：1）根據(jù)Q值隨因子個(gè)數(shù)變化關(guān)系判斷。在設(shè)定合適的不確定度后，可逐漸增加因子數(shù)，觀察Q值的變化。當(dāng)因子數(shù)逐漸增加，而對(duì)應(yīng)結(jié)果的Q值沒(méi)有明顯變化，則此時(shí)因子數(shù)增加為合理。當(dāng)因子數(shù)增加到某一特定值（例如pQ值開(kāi)始有明顯變化，則p-1則可能為合適的因子數(shù)。2）根據(jù)化學(xué)組分的加權(quán)殘差大小判斷?；瘜W(xué)組分的加權(quán)殘差（尤其是標(biāo)識(shí)組分）應(yīng)在正負(fù)3以內(nèi)，如果加權(quán)殘差有較大的分布，則應(yīng)重新檢查因子個(gè)數(shù)。3）根據(jù)多次計(jì)算結(jié)果的差異情況判斷。多次運(yùn)行不應(yīng)該得到多個(gè)Q值或因子譜差異性比較大的結(jié)果（如20次運(yùn)行的結(jié)果中，有多個(gè)結(jié)果之間差異較大），如有，則說(shuō)明因子數(shù)的設(shè)定可能不合理。4）根據(jù)回歸分析結(jié)果判斷。將G矩陣與顆粒物濃度數(shù)據(jù)進(jìn)行回歸分析。如果回歸系數(shù)出現(xiàn)負(fù)值，說(shuō)明因子個(gè)數(shù)設(shè)置可能過(guò)多或存在共線性源；將F矩陣的每一個(gè)因子除以相應(yīng)的回歸系數(shù)后分別各自加和，如果大于1.2，則說(shuō)明選擇的因子個(gè)數(shù)不合理。6.2基礎(chǔ)計(jì)算的結(jié)果判斷6.2.1對(duì)基礎(chǔ)計(jì)算得到的多個(gè)計(jì)算結(jié)果，需進(jìn)一步判斷分析。首先，通過(guò)Q理論值（QTheo）和Q計(jì)算值（QTrue）的大小、是否收斂等，初步選擇查看特定的計(jì)算結(jié)果。然后，從Q計(jì)算值與Q理論值差異分析、觀測(cè)值/預(yù)測(cè)值差異分析、殘差分析（Residualanalysis）等對(duì)計(jì)算結(jié)果進(jìn)一步判斷。6.2.2Q計(jì)算值與Q理論值越接近（85%～115%之間），表明計(jì)算結(jié)果更合理。6.2.3通過(guò)“觀測(cè)值/預(yù)測(cè)值差異分析”，分析輸入的觀測(cè)值和計(jì)算得到的預(yù)測(cè)值的相關(guān)性。如果某化學(xué)組分的觀測(cè)值和預(yù)測(cè)值之間有強(qiáng)相關(guān)性，則表明該化學(xué)組分計(jì)算較好；反之，則考慮是否應(yīng)該降低該化學(xué)組分計(jì)算權(quán)重或從計(jì)算中排除。此外，還可以在時(shí)間序列上分析化學(xué)組分的觀測(cè)值和預(yù)測(cè)值的差異。對(duì)于預(yù)測(cè)值明顯高于觀測(cè)值，需要進(jìn)一步分析確定是否排除該數(shù)據(jù)。6.2.4通過(guò)“殘差分析”，查看每一個(gè)化學(xué)組分的加權(quán)殘差（通過(guò)不確定度加權(quán)）。比如，通過(guò)殘差直方圖中不同化學(xué)組分加權(quán)殘差的百分比，判斷各化學(xué)組分在基礎(chǔ)計(jì)算的情況。如果殘差直方圖顯示殘差范圍在（-3，+3）之間，則該化學(xué)組分的模型計(jì)算在數(shù)學(xué)意義上較好；如果化學(xué)組分具有許多較大的殘差或顯示非正常高斯分布曲線，則表示該化學(xué)組分的模型計(jì)算在數(shù)學(xué)意義上較差。6.3基礎(chǔ)計(jì)算的結(jié)果分析6.3.1基礎(chǔ)計(jì)算的結(jié)果主要包括因子譜F矩陣和因子貢獻(xiàn)G矩陣。因子譜識(shí)別是PMF模型計(jì)算的關(guān)鍵，可通過(guò)分析不同因子譜中標(biāo)識(shí)組分，將因子譜識(shí)別為具體的污染源類。顆粒物源類識(shí)別方法及主要源類的標(biāo)識(shí)組分見(jiàn)附錄A。6.3.2當(dāng)通過(guò)因子譜標(biāo)識(shí)組分難以識(shí)別因子時(shí)，可分析因子貢獻(xiàn)G矩陣的時(shí)間序列貢獻(xiàn)情況，從而識(shí)別出在時(shí)間序列上有特征貢獻(xiàn)變化的污染源類。6.4基礎(chǔ)計(jì)算的誤差評(píng)估6.4.1通過(guò)誤差評(píng)估，分析計(jì)算結(jié)果的不確定性。結(jié)果誤差評(píng)估方法主要有Bootstrap（BS）誤差評(píng)估、Displacement（DISP）誤差評(píng)估、Bootstrap-Displacement（BS-DISP）誤差評(píng)估等。優(yōu)先采用BS誤差評(píng)估方法，評(píng)估模型計(jì)算結(jié)果的可靠性。6.4.2通過(guò)分析BS誤差評(píng)估得到的因子和基礎(chǔ)計(jì)算得到的因子匹配程度，評(píng)估因子中主要化學(xué)組分分配的合理性。如果兩者匹配不到80%，則這個(gè)因子的主要化學(xué)組分分配可能不當(dāng)，需要通過(guò)調(diào)整因子數(shù)量或調(diào)整輸入數(shù)據(jù)等提高匹配程度。6.4.3通過(guò)分析BS誤差評(píng)估得到的Q修正值分布情況，評(píng)估基礎(chǔ)計(jì)算中Q值的合理性。BS計(jì)算可得到Q修正值的最小、最大、中位數(shù)，以及第25和第75百分位數(shù)值?；A(chǔ)計(jì)算得到的Q值一般要求6在第25和第75百分位數(shù)值之間。6.4.4通過(guò)分析BS誤差評(píng)估得到的每個(gè)化學(xué)組分計(jì)算結(jié)果分布情況，評(píng)估化學(xué)組分計(jì)算結(jié)果的不確定性。比如一個(gè)化學(xué)組分的箱型圖上下距（第25分位與第75百分位數(shù)）較寬（超過(guò)20%），則表明該化學(xué)組分計(jì)算結(jié)果的不確定性可能較高；反之，則表明該計(jì)算結(jié)果的不確定性可能較低。7旋轉(zhuǎn)計(jì)算7.1當(dāng)模型基礎(chǔ)計(jì)算得到的因子譜難以識(shí)別為實(shí)際源類的時(shí)候，可使用旋轉(zhuǎn)計(jì)算（Rotationalrun），使因子譜中的標(biāo)識(shí)組分更加突出，促進(jìn)因子譜識(shí)別為具體污染源類。旋轉(zhuǎn)計(jì)算包括F矩陣峰值模型計(jì)算（Fpeakmodelrun）、約束模型計(jì)算（Constrainedmodelrun）。7.2F矩陣峰值模型計(jì)算是常用的旋轉(zhuǎn)計(jì)算方法，適用于標(biāo)識(shí)組分的占比在各因子譜中的占比比較平均，不利于源類識(shí)別的情況。通過(guò)F矩陣峰值模型旋轉(zhuǎn)計(jì)算，調(diào)整F矩陣中的各因子譜中化學(xué)組分占比，使其中的標(biāo)識(shí)組分的占比更加突出。7.3F矩陣峰值模型計(jì)算中，需要嘗試不同的Fpeak參數(shù)。當(dāng)Fpeak參數(shù)設(shè)為正值或負(fù)值時(shí)，F(xiàn)矩陣中的化學(xué)組分在某些因子中占比會(huì)分別更加突出或趨向平均。查看不同F(xiàn)peak參數(shù)下Q值的變化情況，選擇Q值變化在5%以內(nèi)，并且Q值拐點(diǎn)之前Fpeak值對(duì)應(yīng)的計(jì)算結(jié)果，參照“6.2、6.3、6.4”進(jìn)行相應(yīng)的結(jié)果判斷、分析和誤差評(píng)估。7.4PMF共線性問(wèn)題及解決方法。經(jīng)過(guò)旋轉(zhuǎn)計(jì)算后，如果提取的單個(gè)因子明顯包含不同污染源類或子源類（如道路揚(yáng)塵、施工揚(yáng)塵等）的信息，則PMF模型計(jì)算出現(xiàn)了共線性問(wèn)題，可進(jìn)行二次解析。另外，可采用源成分譜約束的PMF模型進(jìn)一步解析。8顆粒物源貢獻(xiàn)計(jì)算8.1通過(guò)因子貢獻(xiàn)G矩陣與顆粒物質(zhì)量濃度進(jìn)行多元線性擬合分析，將得到的回歸系數(shù)乘以G矩陣各列，即為F矩陣識(shí)別的各污染源貢獻(xiàn)。擬合公式如下：PM2.5(j)=skgjk式中：PM2.5(j)——第j個(gè)樣品的顆粒物質(zhì)量濃度；sk——回歸系數(shù)；gjk——PMF模型計(jì)算的G矩陣中第k源對(duì)第j個(gè)樣品的歸一化貢獻(xiàn)；p——污染源個(gè)數(shù)。8.2多個(gè)環(huán)境受體監(jiān)測(cè)點(diǎn)位的模型計(jì)算。將所有受體監(jiān)測(cè)點(diǎn)位數(shù)據(jù)一起輸入PMF模型，計(jì)算得到該工作區(qū)域的總體源貢獻(xiàn)結(jié)果；通過(guò)對(duì)其中特定點(diǎn)位全部樣品的源貢獻(xiàn)結(jié)果進(jìn)行平均，得到該點(diǎn)位的源貢獻(xiàn)結(jié)果。8.3填寫(xiě)PMF模型計(jì)算記錄表，并附上參與計(jì)算的全部數(shù)據(jù)。PMF模型計(jì)算記錄表格式參見(jiàn)附錄B。9結(jié)果合理性判斷9.1結(jié)果需符合模型計(jì)算要求，即各化學(xué)組分的模型計(jì)算結(jié)果與實(shí)測(cè)結(jié)果接近?？赏ㄟ^(guò)直接比較顆粒物及化學(xué)組分質(zhì)量濃度計(jì)算值與實(shí)測(cè)值、分析主要化學(xué)組分或標(biāo)識(shí)組分的殘差、對(duì)比Q計(jì)算值與Q理論值等方法，判斷模型計(jì)算結(jié)果的好壞。如果計(jì)算值與實(shí)測(cè)值越接近、化學(xué)組分殘差越小、Q計(jì)算值與Q理論值越接近，則表明計(jì)算結(jié)果更合理。79.2結(jié)果需符合當(dāng)?shù)貙?shí)際情況，即經(jīng)過(guò)識(shí)別的源成分譜特征及其源貢獻(xiàn)要相對(duì)合理。可通過(guò)PMF模型計(jì)算得到的源成分譜與實(shí)測(cè)源成分譜相關(guān)性比較、PMF模型計(jì)算的平均源貢獻(xiàn)排序與CMB模型（Chemicalmassbalancemodel）計(jì)算結(jié)果一致性比較、源貢獻(xiàn)變化趨勢(shì)分析等方法來(lái)評(píng)估模型計(jì)算結(jié)果的合理性。此外，還可通過(guò)其他監(jiān)測(cè)結(jié)果（比如PM2.5和PM10的濃度比值，自動(dòng)監(jiān)測(cè)的SO2和NOx濃度比值，自動(dòng)監(jiān)測(cè)的O3、CO等污染物濃度等）或污染源信息等來(lái)輔助判斷結(jié)果的合理性。8（資料性附錄）顆粒物源類識(shí)別方法及主要源類的標(biāo)識(shí)組分可通過(guò)以下方法，將PMF模型計(jì)算得到的因子譜識(shí)別為具體的顆粒物源類：A.1根據(jù)各種源類的主要標(biāo)識(shí)組分來(lái)識(shí)別。某個(gè)因子譜的標(biāo)識(shí)組分與其他因子譜相比占比較高，則可識(shí)別為相應(yīng)的源類。例如，某因子譜中NH4+、SO42-的占比分別大于30%、50%，高于其他因子譜中NH4+、SO42-的占比，則該因子譜可識(shí)別為二次硫酸鹽；某因子譜中Ca、Si、Al的占比最高，可識(shí)別為揚(yáng)塵。常見(jiàn)源類的標(biāo)識(shí)組分見(jiàn)表A.1。A.2根據(jù)源類排放特有的時(shí)間變化特征來(lái)識(shí)別。某些源類排放具有顯著時(shí)間變化特征，比如機(jī)動(dòng)車源在早晚高峰時(shí)段貢獻(xiàn)大，季節(jié)變化不明顯；北方民用