南京地區(qū)PCB52多介質(zhì)遷移歸趨行為模擬及環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)

劉丹，張圣虎，劉濟(jì)寧,*，姚成，石利利

1. 環(huán)境保護(hù)部南京環(huán)境科學(xué)研究所，南京 210042 2. 南京工業(yè)大學(xué)理學(xué)院，南京 210009

南京地區(qū)PCB52多介質(zhì)遷移歸趨行為模擬及環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)

劉丹1,2，張圣虎1，劉濟(jì)寧1,*，姚成2，石利利1

1. 環(huán)境保護(hù)部南京環(huán)境科學(xué)研究所，南京 210042 2. 南京工業(yè)大學(xué)理學(xué)院，南京 210009

運(yùn)用多介質(zhì)逸度模型對(duì)PCB52在南京地區(qū)空氣、水體、土壤和沉積物中的濃度分布與多介質(zhì)間的遷移、歸趨行為進(jìn)行模擬研究并分析PCB52在研究區(qū)域環(huán)境多介質(zhì)間的遷移通量，確定其在環(huán)境中的主要遷移過(guò)程；結(jié)合文獻(xiàn)中的監(jiān)測(cè)數(shù)據(jù)，對(duì)模型的可靠性進(jìn)行驗(yàn)證；對(duì)模型的輸入?yún)?shù)進(jìn)行靈敏度分析；通過(guò)模擬濃度與實(shí)際濃度的對(duì)比，表明模型在該地區(qū)具有很好的適用性。結(jié)果表明，環(huán)境系統(tǒng)達(dá)到平衡時(shí)，PCB52在沉積物中的含量占其在環(huán)境系統(tǒng)總含量的96.12%；PCB52從環(huán)境系統(tǒng)中的消失途徑主要為空氣平流輸出和空氣降解；環(huán)境溫度和空氣平流輸入是影響化合物在環(huán)境相中濃度分布的最主要因素。此外，PCB52初步的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)表明其對(duì)生態(tài)環(huán)境未構(gòu)成嚴(yán)重威脅，但潛在危害仍不容忽視。

PCB52；遷移歸趨；多介質(zhì)模型；風(fēng)險(xiǎn)評(píng)價(jià)

多氯聯(lián)苯(polychlorinated biphenyls, PCBs)是目前污染嚴(yán)重的持久性有機(jī)污染物(POPs)之一[1-4]，由于具有優(yōu)良的化學(xué)惰性、不可燃性和抗熱性，20世紀(jì)70年代前曾被廣泛應(yīng)用于化工、電力工業(yè)和塑料工業(yè)等領(lǐng)域。目前，全球絕大多數(shù)國(guó)家停止生產(chǎn)和禁止使用PCBs，雖然我國(guó)也已經(jīng)禁止生產(chǎn)和使用PCBs，但是由于我國(guó)PCBs貯存量較大，相應(yīng)的處置技術(shù)與保管措施還不完善，導(dǎo)致環(huán)境污染問(wèn)題比較嚴(yán)重，給人體健康和生態(tài)系統(tǒng)造成了潛在的威脅[5-6]。

長(zhǎng)江三角洲地區(qū)是我國(guó)工農(nóng)業(yè)發(fā)達(dá)、人口密集的區(qū)域之一，同時(shí)也是我國(guó)典型的PCBs污染區(qū)[7]，而南京作為長(zhǎng)江下游承東啟西的樞紐城市，南臨浙江，東濱上海，不僅容易受到周邊工業(yè)化程度較高城市的污染影響，而且區(qū)域內(nèi)有多家電機(jī)工廠、造紙廠和塑料廠，這些潛在的污染源也會(huì)造成PCBs的環(huán)境污染，因此，對(duì)南京地區(qū)PCBs歸趨行為的研究，對(duì)于了解南京地區(qū)PCBs歸趨行為規(guī)律及生態(tài)風(fēng)險(xiǎn)具有重要影響。國(guó)內(nèi)外已有較多PCBs環(huán)境歸趨行為的相關(guān)研究，Devi等[8]對(duì)印度某地區(qū)的PCBs含量進(jìn)行研究，結(jié)果發(fā)現(xiàn)在檢測(cè)的PCBs同系物中，PCB52的檢出濃度最高，其在空氣中的濃度為(0.96～1.6)×10-10mg·m-3，而在土壤中的濃度是空氣中的2～3倍。Pan等[9]研究結(jié)果表明PCB52在滇池水體和沉積物中的濃度分別為(1.3～7.2)×10-5mg·m-3，(0.6～2.4)×10-3mg·g-1，其含量?jī)H次于PCB28。Donald Mackay的逸度模型也被廣泛的應(yīng)用于污染物歸趨行為的研究[10]。Mackay等[11]應(yīng)用多介質(zhì)逸度模型模擬PCB180在丕平湖和安大略湖的環(huán)境歸趨行為，結(jié)果發(fā)現(xiàn)PCB180在不同的湖泊中表現(xiàn)出不同的趨勢(shì)，主要原因在于降解半衰期和空氣-水分配系數(shù)的不同。國(guó)內(nèi)學(xué)者李英明等[12]研究東北某鋼鐵廠及其附近地區(qū)PCBs的遷移規(guī)律和交換通量，構(gòu)建的模型較好地揭示該地區(qū)PCBs的遷移規(guī)律。蓋楠等[13]通過(guò)測(cè)定若爾蓋地區(qū)不同季節(jié)和海拔的大氣和土壤樣品，結(jié)果發(fā)現(xiàn)無(wú)論是冬季還是夏季，PCB28和PCB52的檢出率都是最高，同時(shí)構(gòu)建的模型模擬該地區(qū)環(huán)境介質(zhì)中PCBs的遷移過(guò)程。目前，尚未見(jiàn)利用逸度模型對(duì)南京地區(qū)PCBs歸趨規(guī)律的研究報(bào)道。因此，本文擬運(yùn)用Mackay的三級(jí)多介質(zhì)逸度模型方法，針對(duì)環(huán)境中檢出率和檢出濃度較高的PCB52在南京地區(qū)各環(huán)境介質(zhì)中的主要遷移過(guò)程進(jìn)行模擬，分析PCB52在南京地區(qū)的遷移規(guī)律和濃度分布，并在此基礎(chǔ)上對(duì)南京地區(qū)環(huán)境介質(zhì)中的PCB52進(jìn)行初步的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)，為該地區(qū)的持久性有機(jī)污染物控制和風(fēng)險(xiǎn)評(píng)價(jià)提供科學(xué)依據(jù)[14-15]。

1 研究方法(Research method)

1.1 研究區(qū)域簡(jiǎn)介

南京地處長(zhǎng)江下游中部地區(qū)，位于北緯31°14″至32°37″，東經(jīng)118°22″至119°14″，屬亞熱帶季風(fēng)氣候，雨量充沛，年降水1 200 mm，四季分明，夏季盛行東南風(fēng)，冬季盛行西北風(fēng)，年平均溫度15.3 °C。全市土地面積6 597 km2，長(zhǎng)江南京段長(zhǎng)度約97 km，江南有秦淮河，江北有滁河，為南京市境內(nèi)2條主要長(zhǎng)江支流，其河谷平原為重要農(nóng)業(yè)區(qū)[16]。

1.2 模型框架

三級(jí)多介質(zhì)環(huán)境逸度模型模擬的是穩(wěn)態(tài)、非平衡、流動(dòng)的環(huán)境系統(tǒng)下，假定物質(zhì)在各環(huán)境相間處于非平衡狀態(tài)，模型中考慮污染物的穩(wěn)態(tài)輸入、輸出和在各個(gè)環(huán)境相內(nèi)發(fā)生的反應(yīng)，以及物質(zhì)在相鄰環(huán)境相間的擴(kuò)散與非擴(kuò)散過(guò)程[11]。本文構(gòu)建的PCBs多介質(zhì)逸度模型的基本框架如圖1所示[17]。模型將整個(gè)研究區(qū)域分為空氣、水、土壤和沉積物四個(gè)主環(huán)境相，每個(gè)主相又由若干個(gè)子相組成，其中空氣相包括氣相和氣溶膠，水體相包括水相和懸浮顆粒物，土壤相包括固相、水相和氣相，沉積物相包括固相和水相。

根據(jù)穩(wěn)態(tài)假設(shè)和質(zhì)量平衡關(guān)系，其質(zhì)量平衡方程組可用如下形式表示：

(1)空氣

df1/dt=(E1+GA1CB1+f2D21+f3D31-f1(D12+D13+DR1+DA1))/V1Z1

(2)水

df2/dt=(E2+GA2CB2+f1D12+f3D32+f4D42-f2(D21+D24+DR2+DA2))/V2Z2

(3)土壤

df3/dt=(E3+f1D13-f3(D31+D32+DR3))/V3Z3

(4)沉積物

df4/dt=(E4+f2D24+f3D31-f4(D42+DR4+DA4))/V4Z4式中：Ei為化合物向環(huán)境的排放速率，mol·h-1；GA為平流流入速率，m·h-1；CBi為平流流入濃度，mol·m-3；DRi為反應(yīng)速率D值，mol·(pa·h)-1；DAi為平流流速D值，mol·(pa·h)-1；V為環(huán)境相體積，m3；Z為逸度容量，mol·m-3·pa；f為逸度，pa；i=1～4，分別表示空氣相、水相、土壤相和沉積物相。穩(wěn)態(tài)非平衡模型可以表述為df/dt=0的形式。具體的符號(hào)意義參照文獻(xiàn)[10]。

圖1 多介質(zhì)逸度模型概念圖注：遷移通量的數(shù)值均為-lg T(mol·h-1)。Fig. 1 Mechanisms incorporated into the multimedia fugacity mode1Note: All values of transfer fluxes are defined as -1g T(mol·h-1).

1.3 過(guò)程和參數(shù)識(shí)別

PCB52在環(huán)境系統(tǒng)中所涉及的環(huán)境相和基本的環(huán)境過(guò)程為：氣相所涉及的環(huán)境過(guò)程主要有污染源排放、鄰近區(qū)域的平流流入和流出、土壤和水界面的擴(kuò)散、降解過(guò)程、懸浮顆粒的干濕沉降、向土壤和水的擴(kuò)散；水相所涉及的環(huán)境過(guò)程主要有污染源排放、鄰近水相的平流流入和流出、氣-水界面的擴(kuò)散、氣相中顆粒物的干濕沉降、土壤徑流、水-沉積物界面的擴(kuò)散、沉降和再懸浮以及自身的降解過(guò)程；土壤相所涉及的環(huán)境過(guò)程主要有污染源排放、地表徑流的流入和流出，氣-土壤界面的擴(kuò)散，降雨氣相中顆粒物的干濕沉降土壤水相水的滲漏和自身的降解過(guò)程；沉積物相所涉及的環(huán)境過(guò)程主要有水-沉積物界面的擴(kuò)散、沉降和再懸浮。

收集的數(shù)據(jù)包括用于模擬計(jì)算和模擬驗(yàn)證的2套數(shù)據(jù)，模型計(jì)算的參數(shù)分為PCB52的物理化學(xué)參數(shù)(表1)、南京地區(qū)的主要環(huán)境參數(shù)(表2)以及PCB52排放參數(shù)(表3)。用于模型驗(yàn)證的數(shù)據(jù)來(lái)自于文獻(xiàn)。我國(guó)于1965年開(kāi)始生產(chǎn)多氯聯(lián)苯，到20世紀(jì)80年代初國(guó)內(nèi)基本已停止生產(chǎn)PCBs，因此本文所研究的環(huán)境系統(tǒng)中假設(shè)沒(méi)有直接污染源[18]。大氣相中的污染源全部來(lái)自于上風(fēng)向大氣的輸入，采用南京夏季上風(fēng)向地區(qū)空氣中PCB52的濃度作為模型中大氣平流輸入污染源濃度，用GA1CB1表示；水體平流輸入濃度參考長(zhǎng)江口水體中PCB52的濃度，用GA2CB2表示；土壤相和沉積物相不直接輸入[19]。

表1 PCB52主要物理化學(xué)性質(zhì)(25 ℃)[20]Table 1 Main physicochemical properties of PCB52 (25 ℃)[20]

注：“-”表示無(wú)量綱。

Note: “-”represents dimensionless parameter.

表2 南京地區(qū)主要環(huán)境參數(shù)Table 2 The environmental parameters of Nanjing

注：*表示實(shí)驗(yàn)測(cè)得數(shù)據(jù)。

Note:*obtained from experiments.

表3 南京地區(qū)PCB52的排放參數(shù)Table 3 Emission parameters of PCB52 in Nanjing

1.4 模型求解與驗(yàn)證

用Matlab編程對(duì)上述方程進(jìn)行求解，計(jì)算出PCB52在各環(huán)境相中的濃度，將其與文獻(xiàn)中的實(shí)測(cè)濃度進(jìn)行比較，驗(yàn)證模型結(jié)果，分析模型的可靠性和準(zhǔn)確性。多介質(zhì)模型是由方程和參數(shù)組成的，由于環(huán)境的復(fù)雜性和可變性，隨著時(shí)間的推移，模型參數(shù)不可避免地會(huì)發(fā)生變化。因此，在研究過(guò)程中為了評(píng)價(jià)輸入?yún)?shù)對(duì)模型輸出結(jié)果的影響，對(duì)模型進(jìn)行必要的靈敏度分析。

2 結(jié)果與討論(Results and discussions)

2.1 模型驗(yàn)證及濃度分布

為驗(yàn)證模型的適用性，從文獻(xiàn)中搜集了PCB52在南京地區(qū)大氣、水、土壤和沉積物的實(shí)測(cè)濃度，對(duì)模型的實(shí)測(cè)值和模擬值進(jìn)行比較。由于缺乏研究區(qū)域內(nèi)的測(cè)試數(shù)據(jù)，故參考了研究區(qū)域周?chē)貐^(qū)的測(cè)試數(shù)據(jù)[23-26]，大氣中PCB52的濃度為7.4×10-7mg·m-3，水體中PCB52的濃度為6.0×10-5mg·m-3，土壤中PCB52的濃度為3.0×10-9mg·g-1，沉積物中PCB52的濃度為7.9×10-6mg·g-1。經(jīng)過(guò)模擬計(jì)算，PCB52在大氣、水、土壤和沉積物相間的濃度分別為7.3×10-7mg·m-3、3.2×10-6mg·m-3、1.4×10-6mg·g-1和6.2×10-6mg·g-1。圖2為該區(qū)域PCB52的實(shí)測(cè)濃度和模擬濃度的對(duì)比。由于區(qū)域環(huán)境的復(fù)雜性使模型受到不確定因素的影響，模型預(yù)測(cè)值與實(shí)測(cè)值在同一數(shù)量級(jí)內(nèi)均認(rèn)為合理[27-28]，本次模擬結(jié)果只有土壤相濃度殘差略高于一個(gè)對(duì)數(shù)，其余各相濃度殘差均控制在一個(gè)對(duì)數(shù)范圍內(nèi)，說(shuō)明模擬值與實(shí)測(cè)值具有較好的一致性，表明所建立的三級(jí)多介質(zhì)逸度模型適用于PCB52在南京地區(qū)各環(huán)境相濃度分布的模擬計(jì)算。

模擬結(jié)果還表明，當(dāng)系統(tǒng)達(dá)到穩(wěn)態(tài)時(shí)，PCB52在環(huán)境各相中的濃度分布為沉積物>土壤>水體>大氣，總量分布為沉積物>土壤>大氣>水體。在大氣、水體、土壤和沉積物相的量占總量的比例分別為：0.67%、0.01%、3.20%和96.12%。其中沉積物是PCB52最主要的儲(chǔ)庫(kù)，幾乎占據(jù)了整個(gè)環(huán)境系統(tǒng)中總量。由于PCB52是一種疏水性化合物，使其更容易吸附在沉積物和土壤中而不利于遷移，且PCB52在水體中有直接的輸入污染源，隨著水-沉積物界面的擴(kuò)散和沉降過(guò)程而進(jìn)入沉積物，使得PCB52主要分布于沉積物中[29-30]。

2.2 相間遷移通量和主要遷移過(guò)程

PCB52在空氣、水體、土壤和沉積物相間的遷移過(guò)程主要包括：在氣相中，輸入過(guò)程以空氣的平流輸入為主，輸出過(guò)程以空氣的平流輸出為主；在水體中，輸入過(guò)程以空氣-水界面的干濕沉降為主要途徑，輸出則以水-空氣界面的擴(kuò)散為主；土壤相的大氣中顆粒物的干濕沉降和雨水溶解為主要輸入途徑，輸出途徑以土壤-大氣界面的擴(kuò)散為主；在沉積物相中，污染物來(lái)源為水體-沉積物間的擴(kuò)散、沉降和再懸浮，消失途徑為沉積物相降解。

PCB52在環(huán)境中的遷移通量如圖3所示。由圖可知，輸入途徑主要是空氣的平流輸入，輸出途徑主要是空氣的平流輸出，其界面遷移過(guò)程主要空氣-土壤沉降、土壤-空氣揮發(fā)和水-空氣揮發(fā)[31]。

2.3 靈敏度分析

靈敏度分析是一種分析模型的輸出結(jié)果對(duì)參數(shù)變化的敏感程度的方法，根據(jù)靈敏度系數(shù)的絕對(duì)值大小可以判斷輸入?yún)?shù)對(duì)輸出結(jié)果的影響程度，符號(hào)的正負(fù)則表明參數(shù)對(duì)輸出結(jié)果的影響方向，即正值表示輸出結(jié)果隨輸入?yún)?shù)的增加而增加，負(fù)值則反之。本文通過(guò)對(duì) 58個(gè)主要模型輸入?yún)?shù)(除常數(shù)外)的靈敏度分析表明：影響大氣相中PCB52濃度的參數(shù)主要有環(huán)境溫度、降雨量、亨利常數(shù)和液態(tài)蒸發(fā)焓，靈敏度系數(shù)分別為-19.5%、-9.0%、8.3%和6.7%，如圖4所示；影響水體相中PCB52大氣平流速率、水體中降解速率常數(shù)、環(huán)境溫度和沉積物空隙水分子有效擴(kuò)散系數(shù)，靈敏度系數(shù)分別為18.4%、16.4%、16.0%和-1.4%，如圖5所示；影響土壤相中PCB52濃度的因素主要有土壤中降解速率常數(shù)、大氣平流速率、大氣平流濃度和環(huán)境溫度，靈敏度系數(shù)分別為18.2%、16.5%、14.1%和1.2%，如圖6所示；影響沉積相中PCB52濃度的因素主要有沉積物中降解速率常數(shù)、大氣平流速率、水體中降解速率常數(shù)和分子在沉積物中擴(kuò)散路徑長(zhǎng)，靈敏度系數(shù)分別為13.7%、10.0%、10.0%和-6.9%，如圖7所示。由此可知，對(duì)關(guān)鍵性模型參數(shù)可靠性的保證也是對(duì)模擬結(jié)果可靠性的保證[32]。

2.4 環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)

目前，美國(guó)環(huán)境保護(hù)局(United States Environmental Protection Agency, USEPA)大多采用風(fēng)險(xiǎn)熵值法進(jìn)行風(fēng)險(xiǎn)評(píng)估，即利用多介質(zhì)逸度模型估算出的預(yù)測(cè)暴露濃度(predicted environmental concentration, PEC)與預(yù)測(cè)無(wú)效應(yīng)濃度值(predicted no effect concentration, PNEC)相比，獲得目標(biāo)污染物的風(fēng)險(xiǎn)熵值(risk quotient, RQ)，如果計(jì)算得到的RQ值<1，那么就說(shuō)明這種化學(xué)品生態(tài)風(fēng)險(xiǎn)低，不需要作進(jìn)一步的評(píng)估；反之，如果任何一種生物級(jí)別計(jì)算得到的RQ值均>1，那么需要進(jìn)行更高層次的風(fēng)險(xiǎn)評(píng)估，以確定其風(fēng)險(xiǎn)是否可接受或不可接受，并根據(jù)生態(tài)風(fēng)險(xiǎn)評(píng)估結(jié)果指導(dǎo)風(fēng)險(xiǎn)減少措施的制定[32-33]。

圖3 PCB52在環(huán)境中的遷移通量Fig. 3 Transfer fluxes of PCB52 in the environment

圖4 氣相中PCB52濃度的靈敏度分析Fig. 4 Sensitivities of PCB52 concentration in air

圖5 水相中PCB52濃度的靈敏度分析Fig. 5 Sensitivities of PCB52 concentration in water

圖6 土壤相中PCB52濃度的靈敏度分析Fig. 6 Sensitivities of PCB52 concentration in soil

圖7 沉積物相中PCB52濃度的靈敏度分析Fig. 7 Sensitivities of PCB52 concentration in sediment

表4 南京地區(qū)PCB52的生態(tài)風(fēng)險(xiǎn)熵Table 4 The ecological risk quotient of PCB52 in Nanjing

注：AF為評(píng)估因子，PNEC為預(yù)測(cè)無(wú)效應(yīng)濃度值，PEC為預(yù)測(cè)暴露濃度，RQ為風(fēng)險(xiǎn)熵值。

Note: AF stands for assessment factor; PNEC stands for predicted no effect concentration; PEC stands for predicted environmental concentration; RQ stands for risk quotient.

為正確對(duì)南京地區(qū)PCB52進(jìn)行初步的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)，選擇蚯蚓作為陸生生物受體，大型溞類(lèi)作為水生生物受體，根據(jù)US EPA的ECOTOX數(shù)據(jù)庫(kù)[34]中PCB52對(duì)生物的急性毒性數(shù)據(jù)，用評(píng)估因子(assessment factor, AF) 法[35]，得到其對(duì)不同生物的PNEC值(表4)。

由此可知水生和陸生生態(tài)風(fēng)險(xiǎn)熵值(RQ為9×10-5和9×10-3)都在可接受范圍內(nèi)，雖然PCB52對(duì)生態(tài)環(huán)境未構(gòu)成嚴(yán)重威脅，但由于該物質(zhì)的生態(tài)累積效應(yīng)，其潛在危害仍不容忽視。

綜上，多介質(zhì)逸度模型較好地模擬了南京地區(qū)PCB52的遷移和歸趨行為，結(jié)果表明：沉積物是PCB52的主要儲(chǔ)庫(kù)，占整個(gè)環(huán)境系統(tǒng)總含量的96.12%；PCB52的空氣平流輸入是其主要輸入途徑，而空氣平流輸出和空氣降解是其主要消失途徑，其界面遷移過(guò)程主要是空氣-土壤沉降、土壤-空氣揮發(fā)和水-空氣揮發(fā)；模型的靈敏度分析顯示環(huán)境溫度和空氣平流輸入是影響PCB52在環(huán)境相中濃度分布的最主要參數(shù)；PCB52的生態(tài)風(fēng)險(xiǎn)熵值(RQ<1=表明其對(duì)生態(tài)環(huán)境未構(gòu)成嚴(yán)重影響。

[1] Wu S S, Chan A K Y, Fang J K H, et al. Measuring and monitoring persistent organic pollutants in the context of risk assessment [J]. Marine Pollution Bulletin, 2008, 57(6): 236-244

[2] 楊彥, 王宗慶, 王瓊, 等. 電子垃圾拆解場(chǎng)多環(huán)境介質(zhì)多氯聯(lián)苯(PCBs)污染特征及風(fēng)險(xiǎn)評(píng)估[J]. 生態(tài)毒理學(xué)報(bào), 2014, 9(1): 133-144

Yang Y, Wang Z Q, Wang Q, et al. Characteristics and environment risk assessment of PCBs in multi-media environment of the electronic waste dismantling venues [J]. Asian Journal of Ecotoxicology, 2014, 9(1): 133-144 (in Chinese)

[3] 連子如, 王江濤, 譚麗菊, 等. 青島近海生物體內(nèi)多環(huán)芳烴、多氯聯(lián)苯和有機(jī)氯農(nóng)藥的含量和分布特征[J]. 生態(tài)毒理學(xué)報(bào), 2010, 5(5): 746-751

Lian Z R, Wang J T, Tan L J, et al. Concentrations and distribution characteristic of PAHs, PCBs and OCPs in the marine organisms of Qingdao coastal area [J]. Asian Journal of Ecotoxicology, 2010, 5(5): 746-751 (in Chinese)

[4] 楊中藝, 鄭晶, 陳社軍, 等. 廣東電子廢物處理處置地區(qū)環(huán)境介質(zhì)污染研究進(jìn)展[J]. 生態(tài)毒理學(xué)報(bào), 2008, 3(6): 533-544

Yang Z Y, Zheng J, Chen S J, et al. Advantage of researches on environmental pollutions from E-Waste treating activities in Guangdong, China [J]. Asian Journal of Ecotoxicology, 2008, 3(6): 533-544 (in Chinese)

[5] 李志勇. 大氣中多氯聯(lián)苯的研究進(jìn)展[J]. 安徽農(nóng)業(yè)科學(xué), 2011, 39(18): 11043-11045, 11050

Li Z Y. Research progress of PCBs in the atmosphere [J]. Journal of Anhui Agriculture and Science, 2011, 39(18): 11043-11045, 11050 (in Chinese)

[6] Tasdemir Y, Vardar N, Odabasi M, et al. Concentrations and gas-particle partitioning of PCBs in Chicago [J]. Environment Pollution, 2004, 131: 35-44

[7] 高軍. 長(zhǎng)江三角洲典型污染農(nóng)田土壤多氯聯(lián)苯分布、微生物效應(yīng)和生物修復(fù)研究[D]. 杭州: 浙江大學(xué), 2005: 1-16

Gao J. Distribution, microbial effects and bioremediation of PCBs in typical polluted agricultural soils in the Yangtze River Delta [D]. Hangzhou: Zhejiang University, 2005: 1-16 (in Chinese)

[8] Ningombam L D, Ishwar C Y, Qi S H, et al. Distribution and risk assessment of polychlorinated biphenyls(PCBs) in the remote air and soil of Manipur, India [J]. Environment Earth Science, 2014, 72: 3955-3967

[9] Wan X, Pan X J, Wang B, et al. Distributions, historical trends, and source investigation of polychlorinated biphenyls in Dianchi Lake China [J]. Chemosphere, 2011, 85: 361-367

[10] Mackay D. Multimedia Environment Models: The Fugacity Approach [M]. Second edition. Boca Raton: Lewis Publishers, 2001: 1-141

[11] Mackay D, Hughes L, Powell D E, et al. An updated Quantitative Water Air Sediment Interaction (QWASI) model for evaluating chemical fate and input parameter sensitivities in aquatic systems: Application to D5 (decamethylcyclopentasiloxane) and PCB-180 in two lakes [J]. Chemosphere, 2014, 111: 359-365

[12] 李英明, 張慶華, 楊洪彪, 等. 鋼鐵廠及其周?chē)h(huán)境中二惡英和多氯聯(lián)苯土-氣交換研究[J]. 中國(guó)環(huán)境監(jiān)測(cè), 2010, 26(4): 1-4

Li Y M, Zhang Q H, Yang H B, et al. Estimation on the soil-air exchange of PCDD/Fs and PCBs around an iron and steel making plant [J]. Environmental Monitoring in China, 2010, 26(4): 1-4 (in Chinese)

[13] 蓋楠. 若爾蓋濕地有機(jī)氯農(nóng)藥和多氯聯(lián)苯的生態(tài)地球化學(xué)研究[D]. 北京: 中國(guó)地質(zhì)科學(xué)院, 2012: 33-115

Gai N. Ecological geochemistry of organochorine pesticides and polychlorinated biphenyls in Zoige Alpine Wetland [D]. Beijing: China Academy of Geochemistry, 2012: 33-115 (in Chinese)

[14] 燕文明, 劉凌. 長(zhǎng)江流域生態(tài)環(huán)境問(wèn)題及其成因[J]. 河海大學(xué)學(xué)報(bào): 自然科學(xué)版, 2006, 34(6): 610-613

Yan W M, Liu L. Eco-environmental problems in Yangtze River Basin and their causes of formation [J]. Journal of Hohai University: Natural Science, 2006, 34(6): 610-613 (in Chinese)

[15] 逄勇, 趙棣華, 姚琪, 等. 長(zhǎng)江江蘇段區(qū)域供水水源地水質(zhì)可達(dá)性研究[J]. 水科學(xué)進(jìn)展, 2003, 14(2): 184-188

Pang X, Zhao L H, Yao Q, et al. Reaching function analysis of water quality in large water supply planning in the Yangtze River in Jiangsu Province [J]. Advances in Water Science, 2003, 14(2): 184-188 (in Chinese)

[16] 方利江, 柳敏, 李向陽(yáng), 等. 南京地區(qū)DDTs跨界面遷移與歸趨的研究[J]. 中國(guó)環(huán)境科學(xué), 2011, 31(2): 307-315

Fang L J, Liu M, Li X Y, et al. Simulating the transfer and fate of DDTs in Nanjing [J]. China Environment Science, 2011, 31(2): 307-315 (in Chinese)

[17] 劉丹, 張圣虎, 劉濟(jì)寧, 等. 化學(xué)品多介質(zhì)逸度模型軟件研究進(jìn)展[J]. 環(huán)境化學(xué), 2014, 33(6): 891-900

Liu D, Zhang S H, Liu J N, et al. Studies on software for multimedia fugacity models of chemicals [J]. Environmental Chemistry, 2014, 33(6): 891-900 (in Chinese)

[18] 張志. 中國(guó)大氣和土壤中多氯聯(lián)苯空間分布特征及規(guī)律研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2010: 16-20

Zhang Z. Polychlorinated biphenyls in Chinese air and surface soil: Spatial distribution characteristics and their inherent causes [D]. Harbin: Harbin Institute of Technology, 2010: 16-20 (in Chinese)

[19] 王敏杰. 遼河口濕地PAHs和PCBs的多介質(zhì)歸趨行為研究[D]. 上海: 中國(guó)海洋大學(xué), 2012: 30-32

Wang M J. Study on the environmental behavior of PAHs and PCBs in Liaohe estuarine wetland [D]. Shanghai: Ocean University of China, 2012: 30-32 (in Chinese)

[20] 劉振宇. 遼河流域殘留多氯有機(jī)物的多介質(zhì)環(huán)境模擬[D]. 大連: 大連理工大學(xué), 2005: 39-45

Liu Z Y. Multimedia environment modeling of poylehlorinated organic compounds (PCOCs) in Liaohe River watershed [D]. Dalian: Dalian University of Technology, 2005: 39-45 (in Chinese)

[21] 林承沖, 吳小根. 長(zhǎng)江徑流量特性及其重要意義的研究[J]. 自然雜志, 1999, 21(4): 200-2005

Lin C C, Wu X G. A study on characteristics of runoff amount and its significance at Yangtze [J]. Ziran Zazhi, 1999, 21(4): 200-2005 (in Chinese)

[22] 李昆, 趙高峰, 周懷東, 等. 枯、平、豐水期長(zhǎng)江3條支流表層水中多氯聯(lián)苯的分布特征及風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境科學(xué), 2012, 33(8): 2574-2579

Li K, Zhao G F, Zhou H D, et al. Distribution characteristics and potential risk of PCBs in surface water from three tributaries of Yangtze River in different periods [J]. Environmental Science, 2012, 33(8): 2574-2579 (in Chinese)

[23] 丁倩. 遼東灣地區(qū)大氣和土壤中多氯聯(lián)苯組分特征分析[D]. 大連: 大連海事大學(xué), 2008: 9-47

Ding Q. Component characteristics of PCBs in the atmosphere and soil of Liaodong Bay Area [D]. Dalian: Dalian Maritime University, 2008: 9-47

[24] 王英, 何松潔, 金軍, 等. 無(wú)錫市大氣中類(lèi)二惡英多氯聯(lián)苯污染水平[J]. 環(huán)境科學(xué)與技術(shù), 2013, 6(4): 44-48

Wang Y, He S J, Jin J, et al. Level of dioxin-like polychlorinated biphenyls in atmosphere of Wuxi [J]. 2013, 6(4): 44-48 (in Chinese)

[25] 張利飛, 史雙昕, 董亮, 等. 江蘇省南通市地表水和大氣中多氯聯(lián)苯的污染特征[C]. 持久性有機(jī)污染物全國(guó)學(xué)術(shù)研討會(huì)論文集, 2010: 72-73

[26] 黃宏, 尹方, 吳瑩. 長(zhǎng)江口表層沉積物中多氯聯(lián)苯殘留和風(fēng)險(xiǎn)評(píng)價(jià)[J]. 同濟(jì)大學(xué)學(xué)報(bào): 自然科學(xué)版, 2011, 39(10): 1500-1505

Huang H, Yin F, Wu Y. Residual characteristics and ecological risk assessment of polychlorinated biphenyls in surface sediments from Yangtze estuary and near shore of East China Sea [J]. Journal of Tongji University: Natural Science, 2011, 39(10): 1500-1505 (in Chinese)

[27] 田慧, 郭強(qiáng), 方利江, 等. 半干旱地區(qū)典型POPs硫丹環(huán)境多介質(zhì)遷移歸趨模擬研究-以蘭州河谷盆地為例[J]. 環(huán)境科學(xué)學(xué)報(bào), 2013, 33(4): 1110-1117

Tian H, Guo Q, Fang L J, et al. Simulation of multimedia transfer and fate of endosulfan in semi-arid area: A case study in Lanzhou, a valley-basin city in Western China [J]. Journal of Environmental Sciences, 2013, 33(4): 1110-1117 (in Chinese)

[28] 董繼元, 王式功, 高宏, 等. 蘭州地區(qū)DDT的環(huán)境多介質(zhì)遷移和歸趨模擬[J]. 生態(tài)環(huán)境學(xué)報(bào), 2009, 18(2): 519-522

Dong J Y, Wang S G, Gao H, et al. Simulation of multi-media transfer and fate of DDT in the Lanzhou Region [J]. Ecology and Environmental Sciences, 2009, 18(2): 519-522 (in Chinese)

[29] 畢新慧, 徐曉白. 多氯聯(lián)苯的環(huán)境行為[J]. 化學(xué)進(jìn)展, 2000, 12(2): 152-160

Bi X H, Xu X B. Behaviors of PCBs in environment [J]. Progress in Chemistry, 2000, 12(2): 152-160 (in Chinese)

[30] 聶海峰, 趙傳冬, 劉應(yīng)漢, 等. 松花江流域河流沉積物中多氯聯(lián)苯的分布、來(lái)源及風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境科學(xué), 2012, 33(10): 3434-3442

Nie H F, Zhao C D, Liu Y H, et al. Distribution, sources and ecological risk assessment of polychlorinated biphenyl in sediments from Songhua River Basin [J]. Environmental Science, 2012, 33(10): 3434-3442 (in Chinese)

[31] 田慧, 郭強(qiáng), 高宏, 等. 廣州地區(qū)典型多溴聯(lián)苯醚遷移和歸趨行為模擬[J]. 中國(guó)環(huán)境科學(xué), 2014, 34(3): 758-765

Tian H, Guo Q, Gao H, et al. Simulating the transfer and fate of typical PBDEs in Guangzhou [J]. China Environmental Science, 2014, 34(3): 758-765 (in Chinese)

[32] 董繼元, 王式功, 高宏, 等. 蘭州地區(qū)苯并(а)芘的環(huán)境多介質(zhì)遷移和歸趨模擬[J]. 生態(tài)環(huán)境, 2008, 17(6): 2150-2153

Dong J Y, Wang S G, Gao H, et al. Simulation of multi-media transfer and fate of benzo[а]pyrene in Lanzhou Region [J]. Ecology and Environment, 2008, 17(6): 2150-2153 (in Chinese)

[33] Bilitewski B. Global risk-based management of chemical additives II: Risk-based assessment and management strategies [J]. Emerging Organic Contaminants and Human Health, 2013, 23: 47-72

[34] Leeuwen C J, van Vermeire T G. Risk Assessment of Chemicals: An Introduction [M]. Second edition. The Netherlands: Springer, 2007: 1-13

[35] US EPA. ECOTOX Database. [EB/OL]. [2008-07-01]. http://cfpub.epa.gov/ecotox/

[36] 湯薪瑤, 左劍惡, 余忻, 等. 制藥廢水中頭孢類(lèi)抗生素殘留檢測(cè)方法及環(huán)境風(fēng)險(xiǎn)評(píng)估[J]. 中國(guó)環(huán)境科學(xué), 2014, 34(9): 2273-2278

Tang X Y, Zuo J E, Yu X, et al. Residue analysis method and environmental risk assessment of cephalosporin antibiotics in pharmaceutical wastewater [J]. China Environmental Science, 2014, 34(9): 2273-2278 (in Chinese)

◆

Simulating the Transfer and Fate and Environmental Risk Assessment of PCB52 in Nanjing

Liu Dan1,2, Zhang Shenghu1, Liu Jining1,*, Yao Cheng2, Shi Lili1

1. Nanjing Institute of Environmental Sciences, MEP, Nanjing 210042, China 2. College of Sciences, Nanjing Tech University, Nanjing 210009, China

17 November 2014 accepted 19 January 2015

A multimedia fugacity model was applied to simulate the concentrations distribution, transfer and fate of PCB52 in air, water, soil and sediment of Nanjing Area. Meanwhile, the transfer fluxes between different compartments were analyzed in order to infer the main transfer process. Moreover, the simulated concentrations in air, soil and sediment were compared to monitoring data for validation purpose. Input parameters of the model were tested and the key parameters were identified using sensitivity analysis method. The reliability of the model was verified by the agreement between calculated and measured concentrations. When the system reaches equilibrium, sediment compartments were the main reservoirs of PCB52 in Nanjing area because of the mass fraction of PCB52 was 96.12%. The air advection outflow and air degradation were the major elimination routes for PCB52 in Nanjing. The results of sensitivity analysis in this study also indicated that the temperature and air advection inflow had significant influence on concentrations of PCB52 in various media. The preliminary risk assessment suggested that the PCB52 in the study area might pose little harm to the environment, but the latent damage should not be disregarded.

PCB52; transfer and fate; multimedia model; risk assessment

國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)(2013AA06A308)；中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專(zhuān)項(xiàng)(2014)；科研院所技術(shù)開(kāi)發(fā)研究專(zhuān)項(xiàng)(2012)

劉丹(1989-)，女，博士，研究方向?yàn)槎嘟橘|(zhì)暴露預(yù)測(cè)模型與風(fēng)險(xiǎn)評(píng)價(jià)，E-mail: liudan0809@njut.edu.cn；

*通訊作者(Corresponding author)， E-mail: ljn@nies.org

10.7524/AJE.1673-5897.20141117004

2014-11-17 錄用日期：2015-01-19

1673-5897(2015)2-353-09

X171.5

A

劉濟(jì)寧(1977-)，男，博士，副研究員，主要從事化學(xué)品風(fēng)險(xiǎn)評(píng)估研究工作。

劉丹, 張圣虎, 劉濟(jì)寧, 等. 南京地區(qū)PCB52多介質(zhì)遷移歸趨行為模擬及環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)[J]. 生態(tài)毒理學(xué)報(bào), 2015, 10(2): 353-361

Liu D, Zhang S H, Liu J N, et al. Simulating the transfer and fate and environmental risk assessment of PCB52 in Nanjing [J]. Asian Journal of Ecotoxicology, 2015, 10(2): 353-361(in Chinese)