• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2015-06-05 09:51:45劉丹張圣虎劉濟(jì)寧姚成石利利
    生態(tài)毒理學(xué)報(bào) 2015年2期
    關(guān)鍵詞:多氯聯(lián)苯平流沉積物

    劉丹,張圣虎,劉濟(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)

    猜你喜歡
    多氯聯(lián)苯平流沉積物
    2007—2020年威海機(jī)場(chǎng)4—9月平流低云的變化特征及影響因子分析
    晚更新世以來(lái)南黃海陸架沉積物源分析
    水環(huán)境中多氯聯(lián)苯來(lái)源與風(fēng)險(xiǎn)評(píng)價(jià)
    云南化工(2021年7期)2021-12-21 07:27:26
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    水體表層沉積物對(duì)磷的吸收及釋放研究進(jìn)展
    氣相色譜-串聯(lián)質(zhì)譜法測(cè)定PM2.5中7種指示性多氯聯(lián)苯和16種多環(huán)芳烴
    基于衛(wèi)星和Argo觀測(cè)的阿拉伯海中北部海表鹽度季節(jié)和年際變化
    荊州市一次局地濃霧天氣特征分析
    浦東機(jī)場(chǎng)一次低云低能見(jiàn)度天氣氣象服務(wù)總結(jié)
    氣相色譜-串聯(lián)質(zhì)譜法測(cè)定牛奶中多氯聯(lián)苯及多環(huán)芳烴
    好看av亚洲va欧美ⅴa在| 免费电影在线观看免费观看| 午夜亚洲福利在线播放| 成年女人毛片免费观看观看9| 国产激情欧美一区二区| 18禁观看日本| 日韩欧美三级三区| 久久天躁狠狠躁夜夜2o2o| 女人被狂操c到高潮| 精品久久久久久久人妻蜜臀av| 亚洲熟女毛片儿| 国内揄拍国产精品人妻在线| 亚洲国产精品成人综合色| 一进一出抽搐动态| 亚洲 欧美一区二区三区| 日韩中文字幕欧美一区二区| 中文字幕人妻丝袜一区二区| 一本一本综合久久| 欧美一级a爱片免费观看看| 在线a可以看的网站| www日本黄色视频网| 欧美一区二区国产精品久久精品| 一级作爱视频免费观看| 日日夜夜操网爽| 夜夜看夜夜爽夜夜摸| 国产蜜桃级精品一区二区三区| 欧美zozozo另类| 又爽又黄无遮挡网站| 人人妻人人看人人澡| 日本免费一区二区三区高清不卡| 九九热线精品视视频播放| 亚洲国产色片| 在线国产一区二区在线| 精品久久久久久,| 久久中文字幕人妻熟女| 精品午夜福利视频在线观看一区| 亚洲av五月六月丁香网| 18禁国产床啪视频网站| 中国美女看黄片| 国产精品 欧美亚洲| 国产精品国产高清国产av| 亚洲成人中文字幕在线播放| 中亚洲国语对白在线视频| 国产激情久久老熟女| 又黄又爽又免费观看的视频| 日本 av在线| 日韩av在线大香蕉| 久久久久精品国产欧美久久久| 欧美日韩精品网址| 18禁国产床啪视频网站| 一本一本综合久久| 成人欧美大片| 男女午夜视频在线观看| 欧美在线一区亚洲| 国产又黄又爽又无遮挡在线| 久久香蕉精品热| 亚洲激情在线av| 性色avwww在线观看| 高清在线国产一区| 日本黄大片高清| 深夜精品福利| av女优亚洲男人天堂 | 在线观看一区二区三区| 一边摸一边抽搐一进一小说| xxx96com| 国产极品精品免费视频能看的| 久久久久久久久久黄片| 可以在线观看的亚洲视频| 国内精品久久久久精免费| 亚洲精品在线美女| 88av欧美| 美女 人体艺术 gogo| 波多野结衣巨乳人妻| 男人舔奶头视频| a级毛片a级免费在线| 精品一区二区三区四区五区乱码| 国产伦在线观看视频一区| 天天躁狠狠躁夜夜躁狠狠躁| а√天堂www在线а√下载| 成人鲁丝片一二三区免费| 国产伦在线观看视频一区| 搡老岳熟女国产| 国产成人精品无人区| 国产激情久久老熟女| 变态另类成人亚洲欧美熟女| 日韩三级视频一区二区三区| 两人在一起打扑克的视频| 亚洲国产欧美网| 男女午夜视频在线观看| 中国美女看黄片| 男女那种视频在线观看| 国内精品美女久久久久久| 日韩三级视频一区二区三区| 欧美中文日本在线观看视频| 老熟妇乱子伦视频在线观看| 久久精品国产亚洲av香蕉五月| 麻豆av在线久日| 久久国产精品人妻蜜桃| 神马国产精品三级电影在线观看| 一夜夜www| 成人亚洲精品av一区二区| 网址你懂的国产日韩在线| 国产免费av片在线观看野外av| 俺也久久电影网| 天堂网av新在线| 一进一出抽搐gif免费好疼| 19禁男女啪啪无遮挡网站| 国内精品久久久久久久电影| h日本视频在线播放| 他把我摸到了高潮在线观看| 热99re8久久精品国产| 中文资源天堂在线| 九九久久精品国产亚洲av麻豆 | 欧美日韩瑟瑟在线播放| 草草在线视频免费看| 亚洲精品久久国产高清桃花| 熟妇人妻久久中文字幕3abv| 久久天堂一区二区三区四区| 精品熟女少妇八av免费久了| 成人高潮视频无遮挡免费网站| 在线观看一区二区三区| 久久久色成人| 国产伦一二天堂av在线观看| 国产91精品成人一区二区三区| 欧美乱妇无乱码| 国产伦精品一区二区三区四那| 99国产极品粉嫩在线观看| 我要搜黄色片| 亚洲天堂国产精品一区在线| 人人妻,人人澡人人爽秒播| 免费大片18禁| xxx96com| 97人妻精品一区二区三区麻豆| 欧美黑人巨大hd| 久久国产精品人妻蜜桃| 无遮挡黄片免费观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品永久免费网站| 老熟妇乱子伦视频在线观看| 99精品欧美一区二区三区四区| 亚洲成人久久爱视频| 一本综合久久免费| 1024香蕉在线观看| 桃色一区二区三区在线观看| 又紧又爽又黄一区二区| 欧美性猛交黑人性爽| 欧美3d第一页| 欧美日韩黄片免| 成人特级av手机在线观看| 精品久久蜜臀av无| 1000部很黄的大片| 99久久无色码亚洲精品果冻| 成人精品一区二区免费| 国产免费av片在线观看野外av| 精品欧美国产一区二区三| 亚洲中文日韩欧美视频| 久久精品国产清高在天天线| 国产亚洲欧美98| 国产黄片美女视频| 亚洲av日韩精品久久久久久密| 99久国产av精品| 99视频精品全部免费 在线 | 国语自产精品视频在线第100页| 99热只有精品国产| 久久久久国产一级毛片高清牌| 一级黄色大片毛片| 久久久国产精品麻豆| 欧美极品一区二区三区四区| av在线蜜桃| 免费看美女性在线毛片视频| 制服人妻中文乱码| 女人高潮潮喷娇喘18禁视频| 亚洲精品久久国产高清桃花| 欧美日韩国产亚洲二区| 97超视频在线观看视频| 亚洲人成伊人成综合网2020| 国产视频一区二区在线看| 亚洲av日韩精品久久久久久密| 成熟少妇高潮喷水视频| 午夜福利高清视频| 日本 欧美在线| 一级毛片精品| 久久99热这里只有精品18| 国产淫片久久久久久久久 | 亚洲精品456在线播放app | 老司机午夜十八禁免费视频| 亚洲av五月六月丁香网| 国产一区二区三区在线臀色熟女| 99久久成人亚洲精品观看| 国产高清视频在线播放一区| 免费人成视频x8x8入口观看| 中文字幕高清在线视频| 日韩欧美一区二区三区在线观看| 夜夜夜夜夜久久久久| 变态另类成人亚洲欧美熟女| 日本免费a在线| 99久久成人亚洲精品观看| 悠悠久久av| 欧美成人性av电影在线观看| 99国产精品99久久久久| 91麻豆精品激情在线观看国产| 午夜激情福利司机影院| 午夜两性在线视频| 51午夜福利影视在线观看| 真实男女啪啪啪动态图| www日本黄色视频网| 好男人在线观看高清免费视频| 又黄又粗又硬又大视频| 国产一区二区在线av高清观看| 精品久久蜜臀av无| 神马国产精品三级电影在线观看| 欧美成狂野欧美在线观看| 国产精品亚洲美女久久久| 日韩欧美三级三区| 国产成人啪精品午夜网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲黑人精品在线| 久久精品亚洲精品国产色婷小说| 久久香蕉国产精品| 国产高清videossex| 国产男靠女视频免费网站| 国模一区二区三区四区视频 | 国产精品久久久久久人妻精品电影| 欧美一级a爱片免费观看看| 美女大奶头视频| 一本一本综合久久| av福利片在线观看| 日本a在线网址| 淫妇啪啪啪对白视频| 在线视频色国产色| 99久久国产精品久久久| 身体一侧抽搐| 欧美成人一区二区免费高清观看 | 免费搜索国产男女视频| 一本综合久久免费| 99在线人妻在线中文字幕| 91麻豆av在线| 免费看a级黄色片| 国产亚洲av嫩草精品影院| 天天一区二区日本电影三级| a级毛片a级免费在线| 制服人妻中文乱码| 搡老岳熟女国产| 操出白浆在线播放| 身体一侧抽搐| 日本一二三区视频观看| 欧美日韩精品网址| 欧美成人性av电影在线观看| 国产伦一二天堂av在线观看| 波多野结衣高清作品| 午夜福利成人在线免费观看| 国产野战对白在线观看| 欧美成狂野欧美在线观看| 国产精品国产高清国产av| 欧美日韩中文字幕国产精品一区二区三区| 男女那种视频在线观看| 免费av毛片视频| 国产精品亚洲av一区麻豆| 亚洲中文字幕一区二区三区有码在线看 | 国产久久久一区二区三区| 欧美+亚洲+日韩+国产| 嫁个100分男人电影在线观看| 无限看片的www在线观看| 老汉色av国产亚洲站长工具| 国产成人啪精品午夜网站| 免费看a级黄色片| 欧美色视频一区免费| 国产乱人伦免费视频| 极品教师在线免费播放| 色综合婷婷激情| 99热只有精品国产| 精品免费久久久久久久清纯| 国产av麻豆久久久久久久| 他把我摸到了高潮在线观看| 免费看日本二区| 欧美一区二区精品小视频在线| 两个人视频免费观看高清| 国产精品电影一区二区三区| 久久久成人免费电影| 国产97色在线日韩免费| 国产三级中文精品| 给我免费播放毛片高清在线观看| 成人高潮视频无遮挡免费网站| 午夜福利免费观看在线| 啦啦啦韩国在线观看视频| 午夜免费观看网址| 欧美成人免费av一区二区三区| 国产乱人视频| 麻豆国产av国片精品| 一级黄色大片毛片| 中亚洲国语对白在线视频| 日韩国内少妇激情av| 欧美性猛交黑人性爽| 夜夜躁狠狠躁天天躁| 一区福利在线观看| 成熟少妇高潮喷水视频| 女人高潮潮喷娇喘18禁视频| 精品熟女少妇八av免费久了| 成年女人看的毛片在线观看| 91久久精品国产一区二区成人 | 国产精品98久久久久久宅男小说| 国内毛片毛片毛片毛片毛片| 啦啦啦韩国在线观看视频| 久久中文看片网| 校园春色视频在线观看| 18禁黄网站禁片免费观看直播| 久久这里只有精品19| 精品国产超薄肉色丝袜足j| 国产精品久久久久久人妻精品电影| 一本综合久久免费| 日韩欧美免费精品| 国产伦在线观看视频一区| 欧美日韩中文字幕国产精品一区二区三区| 国产伦精品一区二区三区四那| 亚洲精品色激情综合| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩无卡精品| 成年免费大片在线观看| 亚洲乱码一区二区免费版| 午夜激情欧美在线| 伊人久久大香线蕉亚洲五| 麻豆国产av国片精品| 精品国产亚洲在线| 两性夫妻黄色片| 亚洲人成网站在线播放欧美日韩| 天天躁日日操中文字幕| 成人国产综合亚洲| 男女视频在线观看网站免费| 亚洲国产精品合色在线| 国产成人精品久久二区二区91| 香蕉久久夜色| 俄罗斯特黄特色一大片| 国产亚洲精品一区二区www| 综合色av麻豆| 国产亚洲精品av在线| 国产一区二区在线av高清观看| 国产成人av激情在线播放| 免费av毛片视频| 日本在线视频免费播放| 国产真人三级小视频在线观看| 亚洲成a人片在线一区二区| 欧美+亚洲+日韩+国产| 精品熟女少妇八av免费久了| 老汉色∧v一级毛片| 香蕉国产在线看| 一区二区三区激情视频| 毛片女人毛片| 国产99白浆流出| 美女cb高潮喷水在线观看 | 美女 人体艺术 gogo| 国内少妇人妻偷人精品xxx网站 | 久久久国产欧美日韩av| 精品久久蜜臀av无| 国产又色又爽无遮挡免费看| 又爽又黄无遮挡网站| 色在线成人网| 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| 亚洲欧美激情综合另类| 午夜福利18| 在线看三级毛片| 欧美+亚洲+日韩+国产| 国产麻豆成人av免费视频| 一进一出抽搐gif免费好疼| 最近在线观看免费完整版| 日本黄大片高清| 真人一进一出gif抽搐免费| 精品一区二区三区视频在线观看免费| 精品欧美国产一区二区三| 999精品在线视频| 好男人在线观看高清免费视频| 亚洲一区二区三区色噜噜| 国产亚洲av嫩草精品影院| 日本黄色视频三级网站网址| 嫁个100分男人电影在线观看| 国产亚洲精品久久久久久毛片| 麻豆av在线久日| 小说图片视频综合网站| 男人的好看免费观看在线视频| 精品不卡国产一区二区三区| 亚洲精品久久国产高清桃花| 国产亚洲精品久久久久久毛片| 午夜免费激情av| 午夜两性在线视频| 国产欧美日韩一区二区精品| 麻豆久久精品国产亚洲av| 精品欧美国产一区二区三| 两个人的视频大全免费| 一进一出抽搐动态| 久久人人精品亚洲av| 国产亚洲欧美98| 成人三级黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美一区二区三区在线观看| 日韩欧美国产在线观看| 五月玫瑰六月丁香| 天堂av国产一区二区熟女人妻| 九九热线精品视视频播放| 亚洲av片天天在线观看| 一夜夜www| 手机成人av网站| 久久亚洲真实| 亚洲中文字幕日韩| 精品一区二区三区视频在线 | 久久久国产精品麻豆| 欧美高清成人免费视频www| av在线天堂中文字幕| 熟妇人妻久久中文字幕3abv| av视频在线观看入口| 国产乱人伦免费视频| 又大又爽又粗| 两个人的视频大全免费| 亚洲欧美激情综合另类| 五月伊人婷婷丁香| 精品久久久久久久久久久久久| 国产91精品成人一区二区三区| 亚洲在线自拍视频| 午夜福利免费观看在线| 日韩国内少妇激情av| 一个人观看的视频www高清免费观看 | 色综合亚洲欧美另类图片| 国产 一区 欧美 日韩| 在线永久观看黄色视频| 亚洲在线自拍视频| a级毛片a级免费在线| 亚洲专区字幕在线| 欧美黑人巨大hd| 亚洲中文字幕日韩| 18禁黄网站禁片午夜丰满| 精品人妻1区二区| 精品久久蜜臀av无| 亚洲激情在线av| 亚洲欧美日韩高清在线视频| 床上黄色一级片| 亚洲av成人一区二区三| 国产成人精品久久二区二区91| 变态另类成人亚洲欧美熟女| 久久午夜亚洲精品久久| 神马国产精品三级电影在线观看| 精品久久蜜臀av无| 1024香蕉在线观看| 日日夜夜操网爽| 国产精品精品国产色婷婷| 熟妇人妻久久中文字幕3abv| 男女视频在线观看网站免费| 成人特级黄色片久久久久久久| 别揉我奶头~嗯~啊~动态视频| 欧美+亚洲+日韩+国产| 又爽又黄无遮挡网站| 日本 av在线| 久久久久久久久免费视频了| 女生性感内裤真人,穿戴方法视频| 国产亚洲精品av在线| 欧美zozozo另类| 亚洲中文字幕日韩| 99热6这里只有精品| 亚洲成人中文字幕在线播放| 国产一区二区在线av高清观看| 18禁黄网站禁片午夜丰满| av中文乱码字幕在线| 成人亚洲精品av一区二区| 国产精品免费一区二区三区在线| 麻豆一二三区av精品| 免费观看人在逋| 国产精品1区2区在线观看.| 国产亚洲精品久久久com| ponron亚洲| 久久久国产成人精品二区| 国产伦精品一区二区三区视频9 | 国产精品99久久久久久久久| 男女之事视频高清在线观看| 中文亚洲av片在线观看爽| 成人性生交大片免费视频hd| 亚洲国产精品999在线| 亚洲国产精品久久男人天堂| 久久久水蜜桃国产精品网| 99久久99久久久精品蜜桃| 欧美zozozo另类| 国产极品精品免费视频能看的| x7x7x7水蜜桃| 久久婷婷人人爽人人干人人爱| 最新美女视频免费是黄的| 国产成人福利小说| 99国产精品一区二区三区| 国模一区二区三区四区视频 | 一卡2卡三卡四卡精品乱码亚洲| 欧美午夜高清在线| 麻豆一二三区av精品| av视频在线观看入口| 精品国产超薄肉色丝袜足j| 免费在线观看成人毛片| 久久精品国产亚洲av香蕉五月| 啦啦啦韩国在线观看视频| 观看美女的网站| 99视频精品全部免费 在线 | 99久久精品热视频| 国产1区2区3区精品| 最近最新中文字幕大全免费视频| 欧美三级亚洲精品| 国产97色在线日韩免费| 色吧在线观看| 欧美高清成人免费视频www| 一个人观看的视频www高清免费观看 | 国产成年人精品一区二区| 久久精品国产综合久久久| 亚洲精华国产精华精| 国产高清视频在线播放一区| 99久久综合精品五月天人人| 一二三四社区在线视频社区8| 国产亚洲精品av在线| 欧美成人一区二区免费高清观看 | 校园春色视频在线观看| 99在线人妻在线中文字幕| 日本熟妇午夜| 亚洲色图 男人天堂 中文字幕| av天堂在线播放| 国内久久婷婷六月综合欲色啪| 变态另类丝袜制服| 又黄又爽又免费观看的视频| av女优亚洲男人天堂 | 国产成人aa在线观看| 91av网一区二区| 亚洲在线自拍视频| 性色av乱码一区二区三区2| 国产亚洲精品久久久久久毛片| aaaaa片日本免费| 国产精品久久久久久精品电影| 国产av麻豆久久久久久久| cao死你这个sao货| 精品电影一区二区在线| 熟女少妇亚洲综合色aaa.| 十八禁网站免费在线| 欧美乱色亚洲激情| 美女cb高潮喷水在线观看 | 最新美女视频免费是黄的| 亚洲性夜色夜夜综合| 亚洲专区国产一区二区| 99热6这里只有精品| 色综合亚洲欧美另类图片| 国产亚洲精品综合一区在线观看| 这个男人来自地球电影免费观看| 午夜免费成人在线视频| 中文在线观看免费www的网站| 中文字幕熟女人妻在线| 国产精品亚洲一级av第二区| 最新在线观看一区二区三区| 久久中文看片网| 国产精品99久久久久久久久| 好男人在线观看高清免费视频| av国产免费在线观看| 脱女人内裤的视频| 成人午夜高清在线视频| 日本免费一区二区三区高清不卡| 亚洲国产精品久久男人天堂| 国产一区二区三区在线臀色熟女| 亚洲熟女毛片儿| 欧美中文日本在线观看视频| 长腿黑丝高跟| 婷婷亚洲欧美| 非洲黑人性xxxx精品又粗又长| 草草在线视频免费看| 亚洲午夜理论影院| 美女cb高潮喷水在线观看 | 亚洲av成人精品一区久久| 日本免费a在线| av在线蜜桃| 精品久久久久久成人av| 露出奶头的视频| 国产精品亚洲一级av第二区| 久久久国产成人精品二区| 长腿黑丝高跟| h日本视频在线播放| 在线看三级毛片| 亚洲精品色激情综合| 欧美激情久久久久久爽电影| bbb黄色大片| 色综合亚洲欧美另类图片| 国产伦在线观看视频一区| 好男人在线观看高清免费视频| 男女下面进入的视频免费午夜| 看片在线看免费视频| 99riav亚洲国产免费| 亚洲国产精品sss在线观看| 窝窝影院91人妻| 亚洲精华国产精华精| 亚洲乱码一区二区免费版| 精品国产乱子伦一区二区三区| 夜夜爽天天搞| 99国产极品粉嫩在线观看| 亚洲在线自拍视频| 国产精品久久久久久人妻精品电影| 亚洲av成人不卡在线观看播放网| 国产成人啪精品午夜网站| 黄色视频,在线免费观看| 色哟哟哟哟哟哟| 国产久久久一区二区三区| 亚洲九九香蕉| 国产高潮美女av| 级片在线观看| 亚洲欧美一区二区三区黑人| 老熟妇仑乱视频hdxx| 99久久成人亚洲精品观看| 国产黄色小视频在线观看| 国内少妇人妻偷人精品xxx网站 | 国产一区二区三区在线臀色熟女| 成年女人看的毛片在线观看| 精品99又大又爽又粗少妇毛片 | 国语自产精品视频在线第100页| av片东京热男人的天堂| 天堂√8在线中文| 在线观看66精品国产| 亚洲第一欧美日韩一区二区三区| 久久香蕉国产精品| 亚洲中文字幕日韩| 美女扒开内裤让男人捅视频|