范 巍,李興春

(中國(guó)石油集團(tuán)安全環(huán)保技術(shù)研究院,北京100049)

1 引言

2 受體模型的定義

源解析所用數(shù)學(xué)模型可以分為2種。一種是以污染源為研究對(duì)象的擴(kuò)散模型(diffusion m odel);另一種是以污染區(qū)域?yàn)檠芯繉?duì)象的受體模型(receptor mode1)。對(duì)前者的研究及其應(yīng)用較多,而對(duì)后者的研究及應(yīng)用主要集中在大氣污染物。

(1)擴(kuò)散模型根據(jù)各個(gè)污染源的排放量、研究區(qū)域與排放源的水平與垂直距離、污染物的理化性質(zhì),以及風(fēng)速風(fēng)向、湍流等環(huán)境因素來(lái)計(jì)算各個(gè)源對(duì)研究區(qū)域的影響程度。模型必需的輸入數(shù)據(jù)是源排放清單,如果輸入數(shù)據(jù)不準(zhǔn)確性,會(huì)造成預(yù)測(cè)誤差,而由于排放源清單存在很大的不準(zhǔn)確性,模型的預(yù)測(cè)效果不能使人信服。

(2)受體模型就是通過(guò)測(cè)量排放源和受體樣品的物理、化學(xué)性質(zhì),定性識(shí)別對(duì)受體有貢獻(xiàn)的污染源并定量計(jì)算各污染源的分擔(dān)率,著眼于研究排放源對(duì)受到排放源影響的局部環(huán)境介質(zhì)的貢獻(xiàn)。受體模型不依賴于排放源的排放條件、氣象、地形等數(shù)據(jù),不追蹤污染物的遷移過(guò)程,避開(kāi)了擴(kuò)散模型遇到的困難,因此得到了更廣泛的應(yīng)用。

3 受體模型技術(shù)的研究進(jìn)展

受體模型的基礎(chǔ)是質(zhì)量守恒原理。如果各源排放物之間沒(méi)有相互反應(yīng),傳輸過(guò)程中沒(méi)有引起某種物質(zhì)的消除和形成,則受體中元素(或化合物)含量是各源貢獻(xiàn)的線性加和。

3.1 化學(xué)質(zhì)量平衡法(CMB)

CMB是是美國(guó)環(huán)?？偸?USEPA)推薦使用的受體模型,并開(kāi)發(fā)出了相應(yīng)的應(yīng)用軟件包[2]。

3.1.1 CMB的原理

設(shè)采樣分析測(cè)得受體區(qū)域中污染物質(zhì)i的濃度為di(μg/m3),該區(qū)域排放物質(zhì)i的源有P種,若已知某排放源 j排放污染物中物質(zhì)i的含量為xij(μg/mg),則源j對(duì)受體的貢獻(xiàn)值g j(mg/m3)應(yīng)滿足：

源j的貢獻(xiàn)率η為：

研究組患者的住院時(shí)長(zhǎng)、臥床時(shí)長(zhǎng)、醫(yī)療費(fèi)用、健康知識(shí)評(píng)分等一般情況均優(yōu)于常規(guī)組患者的(P<0.05)。具體見(jiàn)表1:

測(cè)定n種污染物就可建立n個(gè)方程,只要測(cè)定項(xiàng)目數(shù)量大于或等于排放源數(shù)目,就可解出一組g j,即各排放源的貢獻(xiàn)值。

3.1.2 CMB的模擬效果檢驗(yàn)

檢驗(yàn)?zāi)Ｐ偷哪康氖球?yàn)證排放源貢獻(xiàn)值的有效性以及方程擬合的好壞程度。方法主要包括源貢獻(xiàn)值擬合優(yōu)度檢驗(yàn)、不定性/相似性組檢驗(yàn)、元素濃度計(jì)算值擬合優(yōu)度檢驗(yàn)等。通常用所得的gj和排放源各組分濃度計(jì)算出預(yù)測(cè)濃度,從而得到預(yù)測(cè)濃度與觀測(cè)濃度的比值,評(píng)價(jià)模型的模擬效果。結(jié)果越接近1,說(shuō)明計(jì)算值與測(cè)量值擬合得越好[3～4]。如果某非揮發(fā)性物質(zhì)的模擬效果較差,可能是由于含有這種物質(zhì)的重要排放源未被考慮。

3.1.3 CM B的特點(diǎn)

CMB從一個(gè)受體樣品的分析項(xiàng)目出發(fā)就可以得到結(jié)果,可以避免大量的樣品采集所帶來(lái)的資金等方面的壓力,且能夠檢測(cè)出是否遺漏了某重要源,可以檢驗(yàn)其他方法的適用性。但同一類排放源排放的成分是有差別的,同一排放源在不同的時(shí)間排放物質(zhì)也不同,而CMB模型對(duì)此沒(méi)有加以區(qū)別;CMB模型假設(shè)從排放源到受體之間,排放的物質(zhì)組成沒(méi)有發(fā)生變化,而實(shí)際上某些物質(zhì)并不滿足該條件;排放源的選擇上存在主觀性和經(jīng)驗(yàn)性。

3.1.4 CMB的應(yīng)用

CMB模型已應(yīng)用于水體[5～7]、沉積物[8～9]等介質(zhì)的研究。污染物種類主要集中在多環(huán)芳烴(PAH s)等有機(jī)污染物。Li等[9]通過(guò)文獻(xiàn)搜集了28種源成分譜并對(duì)其進(jìn)行了修正,運(yùn)用USEPA開(kāi)發(fā)的軟件包CMB8.2對(duì)9種模型進(jìn)行了422次運(yùn)行,對(duì)Lake Calum et沉積物中PAH s的來(lái)源進(jìn)行了解析。這一研究為CMB模型在土壤污染物源解析上的應(yīng)用提供了寶貴的經(jīng)驗(yàn)。

3.2 多元統(tǒng)計(jì)模型(multivariate statistics)

多元方法的基本思路是利用觀測(cè)信息中物質(zhì)間的相互關(guān)系來(lái)產(chǎn)生源成分譜或產(chǎn)生暗示重要排放源類型的因子,主要包括因子分析法、多元線性回歸法等。

3.2.1 因子分析法(factor analysis,FA)及主成分分析(principal component analysis,PCA)

因子分析能把具有復(fù)雜關(guān)系的變量歸結(jié)為數(shù)量較少的幾個(gè)綜合因子。在污染物來(lái)源研究中,通常采集大量(設(shè)為N個(gè))樣品,從每一樣品中分析出若干種(設(shè)為M)化學(xué)成分的濃度,這樣就構(gòu)成了一個(gè)包含N×M個(gè)數(shù)據(jù)的集合。由于同一環(huán)境樣品的組成成分并不相互獨(dú)立,來(lái)自于同一類源的那些成分間存在較強(qiáng)的相關(guān)性,因此可以用P個(gè)因子(P＜M)來(lái)描述原來(lái)的樣品集合[10]。

3.2.2 多元線性回歸法(mu ltivariate linearregression,M LR)

不同行業(yè)污染源排放出的污染物質(zhì)含量差別很大。這些主要由某一類排放源決定的元素(物質(zhì))就稱為示蹤元素(物質(zhì))。測(cè)量受體的某物質(zhì)濃度,并對(duì)各示蹤元素的濃度進(jìn)行多元回歸,回歸系數(shù)可以用于計(jì)算各示蹤元素對(duì)應(yīng)排放源對(duì)受體中該物質(zhì)的貢獻(xiàn)。該方法需要知道各類排放源的最佳示蹤元素。示蹤元素之間必須相互獨(dú)立,否則會(huì)存在共線性問(wèn)題,導(dǎo)致多元線性回歸技術(shù)不穩(wěn)定,但可以通過(guò)因子分析(varimax旋轉(zhuǎn))來(lái)實(shí)現(xiàn)。

3.2.3 多元統(tǒng)計(jì)方法的特點(diǎn)

與CMB相比,多元統(tǒng)計(jì)方法不用事先假設(shè)排放源的數(shù)目和類型,排放源的判定相對(duì)比較客觀;能夠解決次生或易變化物質(zhì)的來(lái)源,能利用除濃度以外的一些參數(shù);研究者只需對(duì)排放源組成有大致的了解,并不需要準(zhǔn)確的源成分譜數(shù)據(jù)。但是這種方法也具有一定的局限性：方法不是對(duì)具體數(shù)值進(jìn)行分析而是對(duì)偏差進(jìn)行處理,如果某重要排放源比較恒定,而其他非重要源具有較大的排放強(qiáng)度變異,可能會(huì)忽略排放強(qiáng)度較大的排放源;氣象學(xué)因素的變化影響較大,為了得到準(zhǔn)確的結(jié)果,需要采集足夠多數(shù)量的樣品;在實(shí)際中一般只能鑒別出5～8個(gè)因子,如果重要排放源類型大于10,這種方法不能提供較好的結(jié)果[11]。

3.2.4 多元統(tǒng)計(jì)方法的應(yīng)用

王學(xué)松等[12]對(duì)徐州城市表層土壤中重金屬含量進(jìn)行了因子分析,鑒別出自然、交通及燃煤等3個(gè)因子,從而把該地區(qū)重金屬元素分為“自然因子”類別元素、“交通因子”類別元素、“燃煤因子”類別元素及混合源類別元素等4種。多元線性回歸的研究報(bào)道并不多見(jiàn)[13～14],還未見(jiàn)到大氣環(huán)境以外的研究實(shí)例。

3.3 分子標(biāo)志物及同位素

分子標(biāo)志物(molecu lar markers,MM),是指可用于來(lái)源相關(guān)關(guān)系研究的指標(biāo)化合物,它們具有確定的化學(xué)結(jié)構(gòu),直接或間接地與母源的變化作用有關(guān)[15～16]。介質(zhì)中特定組分的比值可以作為污染源識(shí)別的線索,即為分子標(biāo)志物參數(shù)(MMI)。MM已經(jīng)成為了追蹤有機(jī)污染物污染源的重要手段[17～18],劉瑞民等[19]利用6種分子標(biāo)志物鑒別出天津地區(qū)土壤中主要的PAH s污染源是煤燃燒源,且在天津絕大部分地區(qū)都有分布,其次是石油源和交通污染源,主要分布在中部和濱海地區(qū),在遠(yuǎn)郊縣的邊界地區(qū)有一定的木柴污染源分布。

4 結(jié)語(yǔ)

我國(guó)在源解析方面的研究工作起步較晚,直到20世紀(jì)80年代后期才開(kāi)始進(jìn)行空氣中懸浮顆粒物來(lái)源解析的研究,取得了一定的成績(jī),特別是在CMB法的應(yīng)用方面。但目前研究?jī)?nèi)容主要體現(xiàn)在大氣環(huán)境的研究,研究介質(zhì)較單一,對(duì)土壤的研究還十分不足,各種方法相結(jié)合的趨勢(shì)也并不明顯。土壤環(huán)境介質(zhì)污染物源解析的研究起步較晚,很多方面需要借鑒大氣顆粒物研究的經(jīng)驗(yàn)。受體模型只能反推出對(duì)某受體貢獻(xiàn)較大的幾類排放源,而不能求出某一具體排放源對(duì)受體的貢獻(xiàn)大小,這對(duì)土壤污染防治工作是不利的。而擴(kuò)散模型可能會(huì)成為解決這一困難的途徑,因此如何找到擴(kuò)散模型與受體模型的結(jié)合點(diǎn)是未來(lái)研究的重要課題之一。

[1]Simeonov Einax J,Tsakovski S,K raft J.Multivariate statistical assessment of polluted soils[J].Cen tral Eu ropean Science Jou rnals,2005,3(1)：1～ 9.

[2]Desert Research Institute(DRI).CM B8 User's Manual：EPA's O ffice of A ir Quality Planning&Standards fOAQPs)[M].W ashington DC：United States Environmen tal Protection Agency,Research T riange Park,2001.

[3]王帥杰,朱 坦.大氣顆粒物源解析技術(shù)研究進(jìn)展[J].環(huán)境污染治理技術(shù)與設(shè)備,2002,3(8)：8～12.

[4]朱 坦,白志鵬,陳 威,等.化學(xué)質(zhì)量平衡受體模型新技術(shù)的應(yīng)用——TEDA大氣顆粒物來(lái)源解析實(shí)例[J].城市環(huán)境與城市生態(tài),1996,9(1)：9～ 14.

[5]黃國(guó)蘭,蕭 航,陳春江,等.化學(xué)質(zhì)量平衡法在水體污染物源解析中的應(yīng)用[J].環(huán)境科學(xué),1999,20(6)：14～17.

[6]Booty WG,Resler O,M cCrimmon C.Mass balancemodelling of priority toxic chemicals within the g reat lakes toxic chemicaldecision support system：RateCon m odel resu lts for Lake Ontario and Lake Erie[J].Environm en talModelling and Softw are,2005,20(6)：671～ 688.

[7]Xie ZQ,Sun LG,Zhang PF,et al.Prelim inary geochem ical eviden ce of groundw ater contam ination in coral islands of Xisha[J].Sou th China Sea：Applied Geochem istry,2005,20(10)：1 848～ 1 856.

[8]Kelley D,Nater EA.Sou rce apportionmen t of lake bed sediments to watersheds in an Upper M ississippi basin using a chem ical mass balancemethod[J].Catena,2000,41(4)：277～ 292.

[9]Li A,Jang JK,Schef PA.Ap plication of EPA CMB8.2m odel for sou rce apportionment of sedim en t PAH s in Lake Calumet[J].Environm en tal Science and Technology,2003,37(13)：2 958～2 965.

[10]柯昌華,金文剛,鐘 秦,等.玉溪市中心城區(qū)TSP源解析研究[J].安全與環(huán)境學(xué)報(bào),2003,3(5)：40～43.

[11]Gordon GE.Recep tor models[J].Environmental Science and Technology,1980,14(7)：792～ 799.

[12]王學(xué)松,秦 勇.徐州城市表層土壤中重金屬環(huán)境風(fēng)險(xiǎn)測(cè)度與源解析[J].地球化學(xué),2006,35(1)：88～94.

[13]胡 勇.空氣中總懸浮顆粒物來(lái)源的研究[J].云南環(huán)境科學(xué),1999,18(4)：37～ 38.

[14]Christensen ER,Zhang X.Sources of polycyclic aromatichyd rocarbons to Lake M ichigan determined from sedimentary records[J].Environm ental Science and Technology,1993,27(1)：139～146.

[15]Sim oneit BRT,Sheng G,Chen X.M olecu lar marker study of extractab le organicmatter in aerosols from urban areas of China[J].A tmospheric Environmen t.General Topics,1991,25(10)：2111～ 2129.

[16]祁士華,傅家謨,盛國(guó)英,等.大氣氣溶膠物質(zhì)來(lái)源研究進(jìn)展[J].環(huán)境污染治理技術(shù)與設(shè)備,1999,7(6)：26～31.

[17]Fabb ri D,Vassura I,Sun CG.Sou rce apportionment of polycyclic aromatic hyd rocarbons in a coastal lagoon bym olecu lar and isotopic characterization[J].Marine Chem istry,2003,84(2)：123～135.

[18]Walker SE,Dickhut RM,Chisholm-BrauseC.M olecu lar and isotopic identification of PAH sources in a highly industrialized urban estuary[J].Organic Geochem istry,2005,36(4)：619～632.

[19]劉瑞民,王學(xué)軍,陶 澍,等.天津表層土壤 PAH s分子標(biāo)志物參數(shù)的空間特征[J].中國(guó)環(huán)境科學(xué),2004,24(6)：684～687.