基于逸度方法的多介質(zhì)環(huán)境模型的應用與展望

張子慕, 楊崇銘, 全紅梅, 陶 然, 劉雪萌, 金光洙,3*

(1.延邊大學 農(nóng)學院;2.延邊大學 理學院;3.延邊大學 地理與海洋科學學院：吉林 延吉 133002)

多介質(zhì)環(huán)境模型可用于解析化學物質(zhì)在環(huán)境中的遷移轉(zhuǎn)化和歸趨研究。1979年,Mackay教授在多介質(zhì)環(huán)境模型框架基礎上,應用質(zhì)量平衡原理,將傳統(tǒng)濃度模型中的濃度用逸度替換,首次提出了一種基于逸度方法的多介質(zhì)環(huán)境模型[1-3]。基于逸度方法的多介質(zhì)環(huán)境模型將復雜的模型結(jié)構(gòu)與計算變得簡便,并能更準確地對污染物環(huán)境行為進行描述,已被廣泛應用到各種尺度的環(huán)境中,自其問世以來就受到人們的重視,多年來被廣泛應用并得到迅速的發(fā)展,目前已經(jīng)被北美、歐洲等國家廣泛應用于全球、地區(qū)與局部的環(huán)境歸趨[4-6]。近年來, 污染物多介質(zhì)環(huán)境歸趨的研究多側(cè)重于有機污染物的多介質(zhì)環(huán)境歸趨模型,由于污染物的環(huán)境行為受到研究地區(qū)的氣候條件、地理位置和工業(yè)活動等多因素的影響,應用基于逸度方法的多介質(zhì)模型可對不同環(huán)境介質(zhì)中污染物的濃度水平和遷移歸趨進行有效模擬,尤其是對持久性有機污染物(POPs)這種具有揮發(fā)性、持久性及長距離遷移能力的污染物,探究其在環(huán)境中的遷移歸趨行為及其生態(tài)風險具有重要意義[7]。該研究對逸度模型進行了簡要介紹,著重對基于逸度方法的多介質(zhì)環(huán)境模型的基本原理,以及在環(huán)境污染物多介質(zhì)歸趨研究中的應用和前景進行闡述。

1 逸度模型的原理及應用

逸度(Fugacity)由Lewis教授于1901年提出,用于描述物質(zhì)從某相中逃逸到相鄰相的趨勢[8]。逸度f與濃度C的關(guān)系表達式為：

C=Z·f,

式中,C為濃度,mol/m3;f為逸度,Pa;Z為逸度容量,mol/(m3·Pa)。

當物質(zhì)在各相間達到平衡時,相鄰相間的逸度相等。如果逸度不等,污染物則從高逸度相向低逸度相移動。Mackay教授基于逸度理論及質(zhì)量平衡原理,開發(fā)了4個等級的逸度模型,分別是Ⅰ級、Ⅱ級、Ⅲ級和Ⅳ級模型,用來模擬化學品在多介質(zhì)環(huán)境中的行為變化,為人們估算化學品在環(huán)境中的暴露量提供了更簡便的方法[9,10]。各模型的分類及特征見表1。Ⅰ級模型為平衡、穩(wěn)態(tài)、非流動的環(huán)境系統(tǒng),系統(tǒng)中無化學品的輸入輸出,沒有化學反應的發(fā)生,各相間分配平衡,是最簡單、理想的逸度模型。Ⅱ級模型為平衡、穩(wěn)態(tài)、流動的環(huán)境系統(tǒng),系統(tǒng)中既有化學品的輸入輸出,也存在一系列的化學反應,與Ⅰ級模型的共同點是各環(huán)境相中化學品處于平衡狀態(tài)。Ⅲ級模型為非平衡、穩(wěn)態(tài)、流動的環(huán)境系統(tǒng),系統(tǒng)中存在穩(wěn)態(tài)輸入輸出,化學品從高逸度介質(zhì)向低逸度介質(zhì)流動,不同介質(zhì)間的逸度不同,但各介質(zhì)內(nèi)的逸度不隨時間的變化而變化。Ⅳ級模型為非平衡、非穩(wěn)態(tài)、流動的環(huán)境系統(tǒng),系統(tǒng)中化學品的輸入速率、濃度、逸度隨時間變化,在Ⅲ級模型的基礎上,更為全面地描述了化學品在環(huán)境介質(zhì)中的遷移轉(zhuǎn)化及降解等過程。相較于Ⅰ級、Ⅱ級模型,Ⅲ級和Ⅳ級模型能更好地表示化學品在真實環(huán)境中的行為變化,目前應用較多的是Ⅲ級和Ⅳ級模型。

在平衡方程中,Vi為介質(zhì)i的體積,m3;Ci為介質(zhì)i中化合物的濃度,mol/m3;Zi為化合物在介質(zhì)i的逸度容量,mol/m3·Pa;I為化合物輸入量,mol/h;Gi為介質(zhì)i的流速,m3/h;Ki為化合物在介質(zhì)i中的降解速率系數(shù),1/h;Ei為介質(zhì)i的排放速率,mol/h;CAi為輸入介質(zhì)i中的化合物濃度,mol/m3;Dij為各環(huán)境相的遷移速率,mol/(h·Pa)。

谷曉悅等[11]應用Ⅲ級逸度模型探究了大連區(qū)域環(huán)境多介質(zhì)中全氟辛烷磺酸(PFOS)的分布及遷移規(guī)律,結(jié)果顯示大氣—土壤、土壤—水是研究區(qū)域環(huán)境相間遷移的主要途徑,大氣和水的平流輸入及輸出是污染的主要來源和損失途徑;Su等[12]利用Ⅲ級逸度模型研究了船舶尾管PAHs的運輸和歸宿,并對年排放量進行了估算,建立了環(huán)境排放清單。陳春麗等[13]基于Ⅲ級逸度模型預測和模擬了鄱陽湖區(qū)環(huán)境多介質(zhì)中PAHs的分布和歸趨行為,結(jié)果顯示土壤相和沉積物相是鄱陽湖區(qū)PAHs主要的匯,由土壤和沉積物的內(nèi)源釋放可能會引起二次污染。廖婷等[14]采用Ⅲ級逸度模型對神農(nóng)架大九湖地區(qū)PAHs的多介質(zhì)歸趨進行了研究,將溫度、降水等氣候因素對環(huán)境介質(zhì)中PAHs轉(zhuǎn)化的影響整合進了模型中,探討了氣候變化對PAHs歸趨的影響,并提出應重視因氣候變化而引起的土壤、水體、沉積物的內(nèi)源釋放所造成的二次污染。隨著研究的逐漸深入,越來越多的專家學者利用Ⅳ級逸度模型對污染物在一個時間段內(nèi)的連續(xù)變化進行模擬,逐漸從穩(wěn)態(tài)模擬轉(zhuǎn)變?yōu)閯討B(tài)模擬。Gao等[15]構(gòu)建了Ⅳ級多介質(zhì)逸度模型,模擬了p,p'-DDT和 γ-HCH在1952—2030年珠江三角洲地區(qū)高溫高濕等特殊氣候條件下隨時間和溫度變化的規(guī)律,結(jié)果表明,有機氯農(nóng)藥濃度分布存在季節(jié)性變化,溫度變化對污染物在環(huán)境中的分配有一定的影響。Ki等[16]應用Ⅳ級逸度模型模擬了韓國8個行政區(qū)域內(nèi)的16種優(yōu)控PAHs的動態(tài)歸趨,并預測到2050年,PAHs造成的潛在致癌風險將增加50%以上,西南部地區(qū)將超過風險閾值,表明了增加PAHs監(jiān)測站和控制化石燃料使用的必要性。Diamond等[17]用Ⅳ級逸度模型模擬了As、PCP和PCBs在某海灣中的環(huán)境歸趨行為,并分析了影響污染源行為的關(guān)鍵因素,為污染物的控制提供了有效措施,同時建立了一個基于Ⅲ級穩(wěn)態(tài)模型的食物鏈模型,用來描述有機和無機污染物在食物鏈中的命運。Foster等[18]將Ⅲ級穩(wěn)態(tài)逸度模型和Ⅳ級非穩(wěn)態(tài)逸度模型用在評估PAHs在城市環(huán)境歸趨研究中,穩(wěn)態(tài)模型發(fā)現(xiàn)較易揮發(fā)的PAHs主要分布在氣態(tài)空氣中,容易發(fā)生水平對流和光化學反應,不宜揮發(fā)的PAHs則在土壤和水體中的比例明顯增加,而且很少一部分會揮發(fā)到空氣中;非穩(wěn)態(tài)模型可以用來研究空氣中羥基自由基的變化以及水平對流速率對于PAHs濃度水平的影響,并提出環(huán)境多介質(zhì)分配以及傳輸參數(shù)的準確量化將有利于提高逸度模型對于PAHs存在狀態(tài)的評估。

2 基于逸度方法的多介質(zhì)環(huán)境模型的應用進展

近年來,研究學者在傳統(tǒng)逸度模型的基礎上對其進行了不斷改進,將逸度模型與地理信息系統(tǒng)、大氣傳輸模型、食物鏈模型及人體暴露模型等相結(jié)合,開發(fā)了多個基于逸度理論的多介質(zhì)環(huán)境模型,使其更適用于不同地區(qū)的化學品歸趨模擬,并將其應用于不同的環(huán)境系統(tǒng)中,不斷對研究區(qū)域進行細化與整合,實現(xiàn)了對化學品空間遷移轉(zhuǎn)化分析[19-21]。

2.1 QWASI模型

1983年,Mackay等[22]基于逸度方法建立了QWASI(Quantitative Water Air Sediment Interaction)模型,該模型是一種穩(wěn)態(tài)、非平衡的多介質(zhì)逸度模型,最初被用于模擬河流及湖泊中污染物的歸趨行為,系統(tǒng)由水、沉積物和空氣3個環(huán)境介質(zhì)組成,污染物的遷移過程包括平流、相間傳輸和降解過程。模型需要輸入?yún)?shù)主要包括環(huán)境參數(shù)及污染物的物化參數(shù),以計算化學物質(zhì)在湖泊環(huán)境體系各介質(zhì)中濃度水平、停留時間及遷移速率等。Woodfine等[23]對QWASI模型進行改善,動態(tài)模擬了某流域內(nèi)重金屬的歸趨與命運,并確定了關(guān)鍵的輸入和損失過程,模型準確性良好,對區(qū)域重金屬污染的治理與污染流域的修復提供了理論依據(jù)支持。Yiu等[24]為了量化某海域冰相生長期汞在冰、水、沉積物之間的遷移,在傳統(tǒng)QWASI模型的基礎上增加了一個冰相,建立了QWASI+ice模型,并建立了冰相質(zhì)量平衡方程,以模擬汞在冰相中的遷移轉(zhuǎn)化規(guī)律。Wang等[25]在原有QWASI模型的基礎上增加了光解層,對廢水潟湖系統(tǒng)中6種新興污染物(Emerging Contaminants)的產(chǎn)生遷移過程進行了探究,并對ECs在系統(tǒng)中的命運進行了預測,研究結(jié)果表明,經(jīng)改進的QWASI模型能夠精準地預測潟湖系統(tǒng)中ECs的潛在濃度,實現(xiàn)了污水處理過程的優(yōu)化與對廢水的風險評估,該模型同樣適用于模擬預測其他同樣暴露在陽光下的水體中其他有機污染物的歸趨命運。

2.2 CoZMo-POP模型

Wania等[26]于2000年開發(fā)了CoZMo-POP模型(Coastal Zone Model for Persistent Organic Pollutants),該模型是一種基于逸度方法的非穩(wěn)態(tài)多介質(zhì)環(huán)境模型,能準確描述持久性有機污染物的跨介質(zhì)行為及季節(jié)變化規(guī)律,可用于研究持久性有機污染物在大氣、水、土壤、植物等多種環(huán)境介質(zhì)中的歸趨命運。隨后,于2006年對模型進行了改進,提出了更為靈活的CoZMo-POP 2模型[27],該模型可模擬有機污染物在大氣、土壤、植物、淡水、沉積物等多達19個介質(zhì)組成的環(huán)境中的行為。Parajulee等[28]利用CoZMO-POP模型評估了加拿大阿爾伯塔頁巖油開發(fā)地區(qū)PAHs的環(huán)境排放量和歸趨,為當?shù)丨h(huán)境監(jiān)測及PAHs污染治理提供了數(shù)據(jù)支持。Lee等[29]利用CoZMo-POP 2模型評估了氣候變化對韓國蔚山市持久性有機污染物的歸趨影響,對不同溫度下PCBs在不同環(huán)境介質(zhì)中的含量進行了模擬,結(jié)果表明,溫度是直接影響持久性有機污染物在各環(huán)境介質(zhì)中分布及行為的重要因素,研究結(jié)果可用于了解包括溫度升高在內(nèi)的各種氣候變化效應對持久性有機污染物環(huán)境行為的影響,對制定持久性有機污染物的長期監(jiān)測計劃和監(jiān)管政策有一定的貢獻。Dayyani等[30]在CoZMo-POP模型的基礎上進行了改進,用來評估某流域內(nèi)多環(huán)芳烴及重金屬在積雪與融雪淡水中的濃度及潛在影響,模型預測結(jié)果與研究區(qū)域的實測濃度相近,驗證了該模型的可行性。

2.3 ACC-HUMAN模型

Czub等[31]基于逸度方法建立了非穩(wěn)態(tài)的ACC-HUMAN模型,該模型可用來評估持久性有機污染物從空氣、水和土壤到人體的生物累積潛勢,根據(jù)生物體的不同,考慮的吸收機制包括食物的攝入、土壤的攝入,吸入、呼吸攝入、飲水、大氣沉降和根部吸收,而消除途徑是消化、排尿、呼氣、泌乳和新陳代謝。同時考慮了生理和環(huán)境參數(shù)以及生長的時間變化。ACC-HUMAN模型成功預測了魚、牛肉和牛奶中由水、空氣和土壤污染引起的化學殘留物,將這些食品、水和空氣中的濃度輸入到人類的質(zhì)量平衡模型中來計算人體組織和乳汁中的化學物質(zhì)濃度。經(jīng)驗證,ACC-HUMAN是預測人類接觸有機化合物生物累積的有效工具,能很好地與其他現(xiàn)有的多介質(zhì)模型聯(lián)系起來模擬預測不同污染物產(chǎn)生的人體暴露。ACC-HUMAN模型曾經(jīng)用于模擬美國中西部地區(qū)[32]和加拿大北極地區(qū)的多氯聯(lián)苯人體暴露狀況。McLachlan等[33]利用ACC-HUMAN模型及BETR Global模型預測了全球范圍內(nèi)人體對空氣、水和土壤中PCB153的暴露量,并探究了母乳中PCB153濃度隨時間的變化,經(jīng)過驗證模型的預測值與實際值偏差較小,準確性較高,并提出污染物歸趨、生物累積和人類暴露的綜合建模是在全球范圍內(nèi)研究持久性有機污染物的有效工具。Whelan等[34]基于非穩(wěn)態(tài)的Oslofjord POP模型和ACC-HUMAN模型對挪威某海灣中的cVMS的濃度及遷移規(guī)律進行了探究,并對魚體內(nèi)中污染物的歸趨命運進行了模擬,以此來評估cVMS在區(qū)域海洋內(nèi)食物鏈中的健康風險。Norstr?m等[35]根據(jù)不同暴露途徑,結(jié)合ACC-HUMAN模型分別從飲食攝入(魚、肉和蛋、乳制品)及呼吸吸入等途徑預測了多氯聯(lián)苯對美國中西部成年男性的暴露量。

2.4 其他模型

Breivik等[36]將逸度模型、CoZMo-POP 2模型、BETR-Global模型及ACC-Human模型進行了整合,用來模擬某一特定環(huán)境區(qū)域的污染物行為與人體暴露,并對其在空間和時間上的變化進行預測,這種集成嵌套的暴露模型(NEM)實用性較高,模擬輸出結(jié)果與實際監(jiān)測結(jié)果擬合度良好,將有助于污染物區(qū)域生態(tài)風險評估和環(huán)境管理。Li等[37]整合了排放源模型(CiP-CAFé)[38]、室內(nèi)環(huán)境模型(ICECRM)[39]、農(nóng)村地區(qū)環(huán)境多介質(zhì)模型(CoZMo-POP 2)[26]、城市環(huán)境多介質(zhì)模型(MUM)[40]和食物鏈暴露模型(ACC-HUMAN[31]RAIDAR[41]),新開發(fā)了一種集成嵌套的環(huán)境多介質(zhì)化學物質(zhì)歸趨和暴露模型(Supper-Expo),首次描繪了排放源到室內(nèi)、城市和農(nóng)村環(huán)境中的多氯聯(lián)苯的分布特征和主要途徑的遷移動態(tài)(圖1)。該研究表明不同的環(huán)境在區(qū)域尺度上對多氯聯(lián)苯的總體歸趨和人體暴露具有不同的控制作用。

圖1 基于逸度方法的排放源、室內(nèi)、城市和農(nóng)村環(huán)境的多介質(zhì)歸趨暴露模型[37]

目前我國在污染物環(huán)境行為歸趨研究中也有了很大的發(fā)展,進行了一些大尺度和區(qū)域性的模型研究。其中以北京大學陶澍課題組在這方面的研究居多,該課題組建立了具有較高空間分辨率的區(qū)域歸趨模擬模型,大大降低了傳統(tǒng)模型的不確定性,并探討了多介質(zhì)模型的參數(shù)靈敏度和不確定性,提出了輸入?yún)?shù)變異系數(shù)標化的靈敏度系數(shù)概念,提出了相對于穩(wěn)態(tài)靈敏度的動態(tài)靈敏度的概念,通過模型的建立獲得了持久性有機污染物在各種環(huán)境介質(zhì)中的殘留水平,確定了主要源匯關(guān)系,闡明了污染物在區(qū)域內(nèi)遷移轉(zhuǎn)化的主要途徑和界面遷移速率[42]。曹紅英等[43]以天津地區(qū)為研究區(qū)域,以PAHs的代表物質(zhì)菲為研究對象,建立了穩(wěn)態(tài)多介質(zhì)逸度模型,主要考慮了水、土壤、沉積物和大氣4個主相,結(jié)果表明,土壤和沉積物是天津地區(qū)菲主要的匯,其中,沉積物中積累了70%的菲,大氣中菲的空間分布主要受人為排放的影響,土壤有機質(zhì)的含量對菲降解速率的影響決定了其在土壤中的分布。Lang等[44]在珠江三角洲地區(qū)建立了16種PAHs的Ⅳ級逸度模型模擬得到的季節(jié)變化趨勢與實測結(jié)果相似,氣態(tài)中的PAHs夏天濃度要高于冬天,而顆粒態(tài)中的濃度小于冬天,其中,溫度和沉降是影響PAHs季節(jié)變化的主要原因,其他因素像干沉降速率、水體徑流、水體和大氣中顆粒分配等也會對季節(jié)差異產(chǎn)生影響,該模型沒有擴展到國家尺度。哈爾濱工業(yè)大學曾開發(fā)了一個國家尺度的多介質(zhì)環(huán)境模型,已成功應用于我國太湖地區(qū)和我國環(huán)境介質(zhì)中α-HCH的長期模擬[45]。

3 研究展望

相較于基于濃度方法的多介質(zhì)模型,基于逸度方法的多介質(zhì)環(huán)境模型在化合物的多介質(zhì)歸趨模擬中有著明顯的優(yōu)勢,根據(jù)其目前發(fā)展及應用情況,該文對基于逸度方法的多介質(zhì)環(huán)境模型提出以下展望：

1) 隨著現(xiàn)代化工業(yè)的急劇發(fā)展,有機污染物的排放已經(jīng)在全球范圍內(nèi)被檢測到,了解和預測其在各種環(huán)境介質(zhì)中的來源、分布、遷移、擴散等過程十分必要。但是對于這類在各種環(huán)境介質(zhì)中都廣泛存在的污染物來說,由于設備、技術(shù)和資金等的限制想要進行全面的監(jiān)測是很難實現(xiàn)的,對于這類污染物,應用多介質(zhì)環(huán)境模型對其在不同環(huán)境介質(zhì)中的濃度水平和遷移歸趨進行模擬是一種有效的手段。

2) 用逸度的概念代替?zhèn)鹘y(tǒng)的濃度,在模型結(jié)構(gòu)上進行了簡化,模型只需要輸入研究區(qū)域環(huán)境參數(shù)和化學品物化參數(shù)即可,極大地簡便了運算,某些參數(shù)可通過熱力學計算獲得,減輕了實測及試驗工作。但隨著環(huán)境模型結(jié)構(gòu)越來越趨于復雜化,模型結(jié)構(gòu)、參數(shù)的優(yōu)化及優(yōu)化后參數(shù)的可識別等問題將會成為未來研究中的重要方面。

3) 基于逸度方法的多介質(zhì)環(huán)境模型可以適用于任意個數(shù)介質(zhì)組成的環(huán)境系統(tǒng),具有廣泛的通用性,可以通過對參數(shù)的不斷優(yōu)化對不同研究地區(qū)、不同研究區(qū)域面積、不同環(huán)境介質(zhì)以及各類污染物進行很好的模擬預測,因此,未來的研究工作應著重在調(diào)查和收集整理全球范圍內(nèi)各種污染物的環(huán)境背景值數(shù)據(jù),并建立數(shù)據(jù)庫,以便研究查閱。

4) 基于逸度方法的多介質(zhì)環(huán)境模型將會與越來越多的模型相結(jié)合,如地理信息系統(tǒng)(GIS)、大氣傳輸模型、食物鏈模型、人體暴露模型等,開發(fā)具有多功能性的多介質(zhì)模型,結(jié)合研究區(qū)域的地理信息分析化學品的多介質(zhì)歸趨,將有助于掌握多介質(zhì)環(huán)境系統(tǒng)內(nèi)污染物污染的時空分布特征,識別主要污染源,進行有效的污染防治,并對其進行生態(tài)與健康風險評價,為化學品的排放及污染控制等提供科學依據(jù)與技術(shù)支持。