• <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国产亚洲站长工具| 夜夜骑夜夜射夜夜干| 中亚洲国语对白在线视频| 黑人巨大精品欧美一区二区蜜桃| 人人妻人人添人人爽欧美一区卜| 国产亚洲精品久久久久5区| 国产高清videossex| 视频区欧美日本亚洲| 日韩欧美国产一区二区入口| 美国免费a级毛片| 99久久99久久久精品蜜桃| 国精品久久久久久国模美| 首页视频小说图片口味搜索| 久久香蕉激情| 国产av又大| 制服诱惑二区| 十分钟在线观看高清视频www| 久久青草综合色| 天堂8中文在线网| 三级毛片av免费| 亚洲中文日韩欧美视频| 国产在视频线精品| a级毛片黄视频| 国产精品免费视频内射| 午夜日韩欧美国产| 久久精品国产综合久久久| 首页视频小说图片口味搜索| 国产精品秋霞免费鲁丝片| 又大又爽又粗| 成年人黄色毛片网站| 黄片大片在线免费观看| 亚洲国产欧美一区二区综合| 老司机影院毛片| 90打野战视频偷拍视频| 精品久久久久久久毛片微露脸| 视频区欧美日本亚洲| 成人免费观看视频高清| 啦啦啦中文免费视频观看日本| 中文字幕最新亚洲高清| 午夜激情久久久久久久| 精品福利观看| 久久av网站| 久久久久国内视频| 人人妻人人添人人爽欧美一区卜| 老司机深夜福利视频在线观看| 精品少妇内射三级| 精品欧美一区二区三区在线| 国产一区二区三区综合在线观看| av线在线观看网站| 亚洲精品国产区一区二| 久久性视频一级片| av电影中文网址| 高清视频免费观看一区二区| 久久国产精品男人的天堂亚洲| 在线观看www视频免费| 飞空精品影院首页| 新久久久久国产一级毛片| 免费少妇av软件| 久久中文看片网| 99精品欧美一区二区三区四区| 国产精品久久久久久人妻精品电影 | 日韩成人在线观看一区二区三区| 99久久人妻综合| 久久久久久久久久久久大奶| 久久天堂一区二区三区四区| 99re6热这里在线精品视频| 男女床上黄色一级片免费看| 亚洲国产毛片av蜜桃av| 91成人精品电影| 丰满饥渴人妻一区二区三| 亚洲av成人不卡在线观看播放网| 亚洲伊人色综图| 99国产极品粉嫩在线观看| 亚洲第一青青草原| 国产av又大| 啦啦啦在线免费观看视频4| 亚洲国产成人一精品久久久| 国产成人精品在线电影| 亚洲精品国产一区二区精华液| 美女国产高潮福利片在线看| 亚洲男人天堂网一区| 欧美成狂野欧美在线观看| 久久人人爽av亚洲精品天堂| 成人国产一区最新在线观看| 亚洲精品中文字幕一二三四区 | 十八禁高潮呻吟视频| 丝袜美腿诱惑在线| 色尼玛亚洲综合影院| 又紧又爽又黄一区二区| 免费在线观看完整版高清| 久久九九热精品免费| 亚洲欧美日韩高清在线视频 | 日韩大片免费观看网站| 国产老妇伦熟女老妇高清| 啦啦啦免费观看视频1| 啪啪无遮挡十八禁网站| 50天的宝宝边吃奶边哭怎么回事| 国产男靠女视频免费网站| 午夜福利在线免费观看网站| 高清av免费在线| 国产精品免费视频内射| 99国产精品免费福利视频| 亚洲人成电影免费在线| 女人爽到高潮嗷嗷叫在线视频| 亚洲第一av免费看| 超色免费av| 亚洲精品久久午夜乱码| 嫩草影视91久久| 俄罗斯特黄特色一大片| 一级毛片电影观看| 亚洲专区中文字幕在线| 高清黄色对白视频在线免费看| 午夜激情久久久久久久| 超碰97精品在线观看| 国产欧美亚洲国产| aaaaa片日本免费| 不卡av一区二区三区| 国产精品香港三级国产av潘金莲| 亚洲av国产av综合av卡| 国产免费视频播放在线视频| 极品少妇高潮喷水抽搐| 欧美日韩成人在线一区二区| 老司机福利观看| 国产高清videossex| 黄网站色视频无遮挡免费观看| 亚洲专区字幕在线| 老司机靠b影院| 日韩免费av在线播放| 久久狼人影院| 国产真人三级小视频在线观看| 日韩 欧美 亚洲 中文字幕| 欧美日韩中文字幕国产精品一区二区三区 | 免费观看人在逋| 老司机影院毛片| 黄色视频,在线免费观看| 久久亚洲真实| 欧美一级毛片孕妇| 99国产精品一区二区三区| 老熟妇仑乱视频hdxx| 99久久人妻综合| 午夜精品国产一区二区电影| 91字幕亚洲| 精品一品国产午夜福利视频| 欧美老熟妇乱子伦牲交| 中亚洲国语对白在线视频| 久久久精品国产亚洲av高清涩受| 电影成人av| 国产有黄有色有爽视频| 亚洲欧美色中文字幕在线| 新久久久久国产一级毛片| 99国产精品免费福利视频| 人成视频在线观看免费观看| 麻豆国产av国片精品| 午夜日韩欧美国产| 黄片小视频在线播放| 午夜精品久久久久久毛片777| 在线观看舔阴道视频| 国产淫语在线视频| 精品国产乱码久久久久久小说| 天天躁日日躁夜夜躁夜夜| h视频一区二区三区| 高清毛片免费观看视频网站 | 国产成人影院久久av| 99精品在免费线老司机午夜| av在线播放免费不卡| 少妇裸体淫交视频免费看高清 | 亚洲黑人精品在线| 中文字幕高清在线视频| 国产免费av片在线观看野外av| 欧美人与性动交α欧美精品济南到| 啦啦啦视频在线资源免费观看| 免费黄频网站在线观看国产| 国产欧美日韩一区二区精品| 欧美日韩精品网址| 国产成人精品在线电影| 久久国产精品人妻蜜桃| 91字幕亚洲| 精品国产乱码久久久久久男人| www.精华液| 每晚都被弄得嗷嗷叫到高潮| 男女下面插进去视频免费观看| 国产一区二区 视频在线| 亚洲精品在线观看二区| 视频在线观看一区二区三区| 夜夜骑夜夜射夜夜干| 亚洲成人国产一区在线观看| 精品欧美一区二区三区在线| 999久久久国产精品视频| 满18在线观看网站| 丝袜美足系列| 国产欧美日韩一区二区三| 欧美日韩福利视频一区二区| 少妇 在线观看| 国产在线免费精品| 国产精品欧美亚洲77777| 美女视频免费永久观看网站| 日本a在线网址| 国产精品一区二区免费欧美| 人人妻,人人澡人人爽秒播| 电影成人av| 99国产精品99久久久久| 欧美黄色片欧美黄色片| av网站在线播放免费| 国产无遮挡羞羞视频在线观看| 国内毛片毛片毛片毛片毛片| 精品一区二区三卡| 大香蕉久久成人网| 淫妇啪啪啪对白视频| 亚洲中文字幕日韩| 精品一区二区三卡| 91av网站免费观看| 免费人妻精品一区二区三区视频| 美女扒开内裤让男人捅视频| 中文字幕高清在线视频| 亚洲成人免费电影在线观看| 丝瓜视频免费看黄片| 亚洲精品美女久久久久99蜜臀| 黄片大片在线免费观看| 中文字幕色久视频| 国产成人啪精品午夜网站| 窝窝影院91人妻| 在线观看www视频免费| 久久精品熟女亚洲av麻豆精品| 叶爱在线成人免费视频播放| 在线播放国产精品三级| 欧美av亚洲av综合av国产av| 中文欧美无线码| 国产在视频线精品| 19禁男女啪啪无遮挡网站| av网站在线播放免费| 老熟妇仑乱视频hdxx| 啪啪无遮挡十八禁网站| 色综合婷婷激情| 国产不卡一卡二| 成人黄色视频免费在线看| 一本综合久久免费| 亚洲精品一二三| 桃花免费在线播放| 国产成人欧美| 这个男人来自地球电影免费观看| 我要看黄色一级片免费的| 91麻豆精品激情在线观看国产 | 国产伦人伦偷精品视频| 大片免费播放器 马上看| 欧美变态另类bdsm刘玥| 操美女的视频在线观看| 中文字幕制服av| 久久久久久久久久久久大奶| 夜夜骑夜夜射夜夜干| 老汉色av国产亚洲站长工具| 正在播放国产对白刺激| 亚洲精品在线观看二区| 别揉我奶头~嗯~啊~动态视频| 成年版毛片免费区| 久久人妻福利社区极品人妻图片| 久久久精品94久久精品| 国产在线观看jvid| 性少妇av在线| 亚洲精品国产区一区二| 欧美日韩成人在线一区二区| 亚洲欧美精品综合一区二区三区| 成在线人永久免费视频| 亚洲精品乱久久久久久| 男男h啪啪无遮挡| 老汉色∧v一级毛片| 欧美黄色片欧美黄色片| 亚洲精华国产精华精| 免费在线观看黄色视频的| 啦啦啦中文免费视频观看日本| 日韩中文字幕欧美一区二区| 搡老岳熟女国产| 国内毛片毛片毛片毛片毛片| 国产精品美女特级片免费视频播放器 | 欧美精品亚洲一区二区| 精品国产一区二区三区四区第35| 纵有疾风起免费观看全集完整版| 熟女少妇亚洲综合色aaa.| 一级黄色大片毛片| 丁香欧美五月| 18禁美女被吸乳视频| 黄色丝袜av网址大全| 老司机靠b影院| 国产伦理片在线播放av一区| 亚洲av第一区精品v没综合| 国产成人免费观看mmmm| 窝窝影院91人妻| 少妇粗大呻吟视频| 日韩视频在线欧美| 久久中文字幕一级| 91麻豆av在线| 另类精品久久| 国产黄频视频在线观看| 亚洲五月色婷婷综合| 亚洲av日韩在线播放| 天天躁日日躁夜夜躁夜夜| 19禁男女啪啪无遮挡网站| 制服人妻中文乱码| 国产一卡二卡三卡精品| 国产片内射在线| 午夜福利欧美成人| 一级片'在线观看视频| 亚洲伊人色综图| 中文字幕人妻熟女乱码| 久久狼人影院| 黄色片一级片一级黄色片| 亚洲精品在线美女| 国产黄频视频在线观看| 下体分泌物呈黄色| av网站免费在线观看视频| 亚洲午夜理论影院| 黑丝袜美女国产一区| 啦啦啦中文免费视频观看日本| 久久中文看片网| 免费av中文字幕在线| 日本欧美视频一区| 成人三级做爰电影| 精品视频人人做人人爽| 亚洲国产中文字幕在线视频| 精品一品国产午夜福利视频| 午夜成年电影在线免费观看| 国产主播在线观看一区二区| 国产男女内射视频| 操美女的视频在线观看| 亚洲专区中文字幕在线| 99久久国产精品久久久| 男人舔女人的私密视频| 国产极品粉嫩免费观看在线| 欧美亚洲日本最大视频资源| 午夜精品久久久久久毛片777| 日本黄色视频三级网站网址 | 国产成人啪精品午夜网站| 搡老岳熟女国产| 在线av久久热| 欧美精品高潮呻吟av久久| 亚洲伊人久久精品综合| 新久久久久国产一级毛片| 久久久精品国产亚洲av高清涩受| 国产在线一区二区三区精| 国产国语露脸激情在线看| 成人av一区二区三区在线看| 久久天躁狠狠躁夜夜2o2o| 热re99久久国产66热| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久久久精品古装| 亚洲av电影在线进入| 成人黄色视频免费在线看| 黄色a级毛片大全视频| 在线亚洲精品国产二区图片欧美| 日韩视频一区二区在线观看| 久久人人爽av亚洲精品天堂| 欧美精品高潮呻吟av久久| 成人手机av| 窝窝影院91人妻| 久久精品熟女亚洲av麻豆精品| 后天国语完整版免费观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产看品久久| 国产熟女午夜一区二区三区| 国产一区二区 视频在线| 男女午夜视频在线观看| 免费看十八禁软件| 黄色成人免费大全| 黄片大片在线免费观看| 美女视频免费永久观看网站| 国产精品美女特级片免费视频播放器 | 1024香蕉在线观看| √禁漫天堂资源中文www| 在线观看免费视频网站a站| 欧美日韩亚洲国产一区二区在线观看 | 久久久久久人人人人人| 午夜老司机福利片| 亚洲少妇的诱惑av| 欧美一级毛片孕妇| 国产av一区二区精品久久| 男女高潮啪啪啪动态图| 美女高潮喷水抽搐中文字幕| 日日爽夜夜爽网站| 日韩熟女老妇一区二区性免费视频| 日日摸夜夜添夜夜添小说| 亚洲中文av在线| 热re99久久国产66热| 老司机影院毛片| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩亚洲综合一区二区三区_| 美国免费a级毛片| 亚洲九九香蕉| 国产老妇伦熟女老妇高清| 欧美精品人与动牲交sv欧美| 国产精品一区二区精品视频观看| 亚洲人成电影观看| 欧美 日韩 精品 国产| 99热国产这里只有精品6| 18禁观看日本| 色视频在线一区二区三区| 国产精品香港三级国产av潘金莲| 亚洲av美国av| 丰满人妻熟妇乱又伦精品不卡| 女人高潮潮喷娇喘18禁视频| 69精品国产乱码久久久| 动漫黄色视频在线观看| 青草久久国产| 黑人巨大精品欧美一区二区mp4| 满18在线观看网站| 亚洲精品一卡2卡三卡4卡5卡| 大型av网站在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久精品国产欧美久久久| 涩涩av久久男人的天堂| 久久精品aⅴ一区二区三区四区| 性少妇av在线| 老熟女久久久| 色婷婷av一区二区三区视频| 男女午夜视频在线观看| 日韩大片免费观看网站| 日韩 欧美 亚洲 中文字幕| 一级片'在线观看视频| 男女无遮挡免费网站观看| 国产男女内射视频| 考比视频在线观看| 9色porny在线观看| 欧美激情极品国产一区二区三区| 亚洲精品av麻豆狂野| 亚洲成av片中文字幕在线观看| 91老司机精品| 女人久久www免费人成看片| 国产精品电影一区二区三区 | 免费在线观看黄色视频的| 亚洲全国av大片| 女警被强在线播放| 91成人精品电影| 国产单亲对白刺激| h视频一区二区三区| 欧美国产精品一级二级三级| 蜜桃在线观看..| 久久热在线av| 欧美激情极品国产一区二区三区| 国产亚洲精品一区二区www | √禁漫天堂资源中文www| 久久精品熟女亚洲av麻豆精品| 多毛熟女@视频| 少妇被粗大的猛进出69影院| 人人澡人人妻人| 久久久水蜜桃国产精品网| 少妇猛男粗大的猛烈进出视频| 99久久国产精品久久久| 免费在线观看影片大全网站| 亚洲人成伊人成综合网2020| 国产精品久久久av美女十八| 国产亚洲午夜精品一区二区久久| 精品乱码久久久久久99久播| 啦啦啦在线免费观看视频4| 国产成人精品在线电影| 久久精品亚洲熟妇少妇任你| 久久精品91无色码中文字幕| 另类精品久久| 久久这里只有精品19| 丁香六月天网| 欧美中文综合在线视频| 国产伦理片在线播放av一区| 久久国产亚洲av麻豆专区| 天天添夜夜摸| 一边摸一边抽搐一进一出视频| 久久免费观看电影| 蜜桃国产av成人99| 极品人妻少妇av视频| 国产精品久久久av美女十八| 国产aⅴ精品一区二区三区波| 精品少妇黑人巨大在线播放| 9191精品国产免费久久| 国产成人精品久久二区二区91| 亚洲欧美激情在线| 女警被强在线播放| 男女无遮挡免费网站观看| 精品第一国产精品| 欧美日韩中文字幕国产精品一区二区三区 | 老熟妇仑乱视频hdxx| aaaaa片日本免费| 国产日韩欧美在线精品| 久久ye,这里只有精品| 十八禁网站网址无遮挡| 成人永久免费在线观看视频 | 久久久久久久大尺度免费视频| 亚洲国产av新网站| 国产又色又爽无遮挡免费看| 国产一区二区 视频在线| 国产精品九九99| 午夜福利影视在线免费观看| 一级,二级,三级黄色视频| 欧美日韩黄片免| 午夜老司机福利片| 国产又色又爽无遮挡免费看| 亚洲精品久久午夜乱码| 国产精品av久久久久免费| 欧美激情久久久久久爽电影 | 国产精品一区二区免费欧美| 国产免费现黄频在线看| 在线 av 中文字幕| 成人手机av| 老熟女久久久| 黑人猛操日本美女一级片| 精品少妇黑人巨大在线播放| 国产成+人综合+亚洲专区| 久久亚洲真实| 女人被躁到高潮嗷嗷叫费观| 亚洲精品乱久久久久久| 久久天堂一区二区三区四区| 91精品国产国语对白视频| 性色av乱码一区二区三区2| 侵犯人妻中文字幕一二三四区| 国产欧美日韩一区二区精品| 亚洲成人免费av在线播放| 免费av中文字幕在线| 捣出白浆h1v1| 免费在线观看视频国产中文字幕亚洲| 亚洲av国产av综合av卡| 久久精品国产a三级三级三级| 一级片'在线观看视频| 亚洲成av片中文字幕在线观看| 美国免费a级毛片| 中国美女看黄片| 午夜精品国产一区二区电影| 色综合欧美亚洲国产小说| 国产伦理片在线播放av一区| 国产在视频线精品| 亚洲人成电影观看| av天堂在线播放| 国产精品 欧美亚洲| 19禁男女啪啪无遮挡网站| 国产亚洲精品一区二区www | 怎么达到女性高潮| 久久人妻福利社区极品人妻图片| 国产亚洲av高清不卡| 日韩有码中文字幕| 亚洲国产看品久久| 久久青草综合色| 性色av乱码一区二区三区2| 妹子高潮喷水视频| 首页视频小说图片口味搜索| 男女床上黄色一级片免费看| 国产在视频线精品| 香蕉国产在线看| 亚洲第一av免费看| 色视频在线一区二区三区| 亚洲欧美一区二区三区黑人| 亚洲精品在线美女| 亚洲精品乱久久久久久| 十八禁网站网址无遮挡| 激情在线观看视频在线高清 | 久久 成人 亚洲| 男女下面插进去视频免费观看| 亚洲精品中文字幕在线视频| 狂野欧美激情性xxxx| 午夜福利,免费看| 不卡av一区二区三区| 成年人黄色毛片网站| 午夜福利视频在线观看免费| 国产人伦9x9x在线观看| 国产男女超爽视频在线观看| 在线观看66精品国产| 免费不卡黄色视频| 亚洲成人国产一区在线观看| 女人爽到高潮嗷嗷叫在线视频| 又大又爽又粗| 男女午夜视频在线观看| 我的亚洲天堂| 51午夜福利影视在线观看| 欧美 日韩 精品 国产| 老司机福利观看| 国产亚洲欧美在线一区二区| 麻豆成人av在线观看| 菩萨蛮人人尽说江南好唐韦庄| 在线永久观看黄色视频| 色播在线永久视频| 十八禁人妻一区二区| 丁香六月天网| 看免费av毛片| 日韩人妻精品一区2区三区| www.自偷自拍.com| 黄色丝袜av网址大全| 成在线人永久免费视频| 中文亚洲av片在线观看爽 | 不卡av一区二区三区| 欧美人与性动交α欧美精品济南到| 巨乳人妻的诱惑在线观看| 国产成人免费无遮挡视频| 国产黄色免费在线视频| videos熟女内射| 91老司机精品| 国产精品久久久久久人妻精品电影 | 中亚洲国语对白在线视频| 老熟女久久久| 青青草视频在线视频观看| 亚洲七黄色美女视频| 亚洲 国产 在线| 亚洲精品国产色婷婷电影| 亚洲九九香蕉| 女人被躁到高潮嗷嗷叫费观| 在线播放国产精品三级| 狠狠婷婷综合久久久久久88av| 超碰97精品在线观看| 国产激情久久老熟女| 亚洲人成77777在线视频| 建设人人有责人人尽责人人享有的| 成人国语在线视频| 捣出白浆h1v1| 又紧又爽又黄一区二区| 日韩大片免费观看网站| 亚洲精品国产精品久久久不卡| 亚洲欧美精品综合一区二区三区| 精品久久久久久电影网| 午夜福利,免费看| 老司机在亚洲福利影院| 菩萨蛮人人尽说江南好唐韦庄| 亚洲久久久国产精品| 老熟妇乱子伦视频在线观看|