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

    南京北郊冬季PM2.5中重金屬的形態(tài)、來(lái)源及健康風(fēng)險(xiǎn)

    2023-10-30 07:15:28孫明洋蓋鑫磊
    關(guān)鍵詞:成人軌跡重金屬

    趙 騰,孫明洋,蓋鑫磊

    (南京信息工程大學(xué)環(huán)境科學(xué)與工程學(xué)院/大氣環(huán)境與裝備技術(shù)協(xié)同創(chuàng)新中心/江蘇省大氣環(huán)境監(jiān)測(cè)與污染控制高技術(shù)研究重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210044)

    重金屬是PM2.5中的有毒成分,可通過(guò)呼吸吸入、皮膚接觸、手口攝入等途徑進(jìn)入人體,引發(fā)多種疾病[1]。研究表明,Cr6+可能引起基因突變,損害肝臟和腎臟;車(chē)輛排放產(chǎn)生的顆粒Pb可能造成DNA損傷,V會(huì)影響神經(jīng)系統(tǒng),Cd會(huì)對(duì)肺部及支氣管造成損傷[2-4]。因此,PM2.5中重金屬的環(huán)境風(fēng)險(xiǎn)受到廣泛關(guān)注。

    重金屬對(duì)環(huán)境的危害不僅與其總量有關(guān),還取決于其化學(xué)形態(tài),不同形態(tài)的重金屬具有不同的環(huán)境行為、環(huán)境毒性及生物可利用性[5]。重金屬的弱酸溶解態(tài)與PM2.5結(jié)合較弱,最容易被生物吸收,毒性較強(qiáng),被認(rèn)為是生物直接可利用形態(tài);可還原態(tài)與可氧化態(tài)與PM2.5結(jié)合較強(qiáng),在一定條件下可被生物吸收,毒性較弱,被認(rèn)為是生物潛在可利用形態(tài);這三者比例越高,其生物可利用性越大,危害越高[6]。而殘?jiān)鼞B(tài)與PM2.5結(jié)合穩(wěn)定,基本不會(huì)被生物吸收,毒性最弱,被認(rèn)為是生物不可利用形態(tài)。研究PM2.5重金屬的化學(xué)形態(tài)對(duì)于科學(xué)評(píng)估其環(huán)境風(fēng)險(xiǎn)具有重要意義。

    來(lái)源及氣象條件是影響PM2.5中重金屬化學(xué)形態(tài)的重要因素,弄清PM2.5及其化學(xué)組分的排放源對(duì)于制定有效的大氣污染控制對(duì)策具有重要意義。正定矩陣因子分解模型(PMF)[7-8]、主成分分析(PCA)[9-12]、化學(xué)質(zhì)量平衡模型(CMB)[13]是目前PM2.5中重金屬來(lái)源解析的常用方法。其中PCA法得到的源解析結(jié)果往往不可避免地得到負(fù)值,從而導(dǎo)致物理意義上的損失,CMB法需要提供排放源成分譜信息,限制了其應(yīng)用。而PMF是美國(guó)國(guó)家環(huán)境保護(hù)局(EPA)推薦的一種定量源解析模型,它不需要詳細(xì)的源成分譜信息,且對(duì)因子分解矩陣添加了非負(fù)約束,使得到的源解析結(jié)果更加合理,已被廣泛應(yīng)用于大氣污染物的源解析[14]。此外,后向軌跡模型(HYSPLIT)及濃度權(quán)重軌跡模型(CWT)也被應(yīng)用于源解析領(lǐng)域,PMF、PCA和CMB等僅能用于解析本地源,HYSPLIT是美國(guó)國(guó)家海洋和大氣管理局的空氣資源實(shí)驗(yàn)室(NOAA)研發(fā)的一種研究空氣污染物區(qū)域傳輸和擴(kuò)散軌跡的源解析方法,但HYSPLIT只能定性解析污染來(lái)源,CWT則可以半定量地給出研究區(qū)域的貢獻(xiàn)大小[15]。

    PM2.5中重金屬健康風(fēng)險(xiǎn)評(píng)價(jià)多采用美國(guó)EPA推薦的健康風(fēng)險(xiǎn)評(píng)價(jià)模型,如LIU等[16]發(fā)現(xiàn)西安PM2.5中As、Cr有潛在致癌風(fēng)險(xiǎn),任萬(wàn)輝等[17]發(fā)現(xiàn)沈陽(yáng)PM2.5中Cr對(duì)兒童有明顯致癌風(fēng)險(xiǎn),KERMANI等[18]發(fā)現(xiàn)伊朗Karaj市PM2.5中Cd有明顯致癌風(fēng)險(xiǎn)。以上研究大多是基于重金屬總量進(jìn)行評(píng)估的,并未考慮化學(xué)形態(tài)的影響。如前文所述,只有弱酸溶解態(tài)、可還原態(tài)和可氧化態(tài)3種形態(tài)的重金屬才會(huì)被生物體直接吸收利用,而殘?jiān)鼞B(tài)重金屬在自然條件下基本不會(huì)被生物體吸收利用[19-21];此外,對(duì)特定污染源及區(qū)域傳輸?shù)慕】碉L(fēng)險(xiǎn)研究較少,但同一種元素往往有多個(gè)來(lái)源,且PM2.5中重金屬來(lái)源不僅有本地貢獻(xiàn),還包括周邊區(qū)域的傳輸,這就意味著周邊地區(qū)對(duì)研究區(qū)域的健康風(fēng)險(xiǎn)有一定影響。因此,筆者在南京北郊采集PM2.5樣品,利用逐級(jí)化學(xué)提取法研究其重金屬的化學(xué)形態(tài),采用基于生物可利用形態(tài)濃度修正的健康風(fēng)險(xiǎn)評(píng)價(jià)模型評(píng)估其健康風(fēng)險(xiǎn),并結(jié)合PMF、HYSPLIT及CWT分析特定污染源及區(qū)域傳輸對(duì)當(dāng)?shù)亟】碉L(fēng)險(xiǎn)的貢獻(xiàn)。研究結(jié)果可為南京北郊PM2.5中重金屬的綜合防治及居民健康指南的制定提供一定依據(jù)。

    1 材料與方法

    1.1 樣品采集

    采樣點(diǎn)位于南京信息工程大學(xué)學(xué)科一號(hào)樓樓頂(32.20° N,118.71° E)(距離地面約20 m),周?chē)鸁o(wú)高大建筑物遮擋。該點(diǎn)位于南京江北新區(qū),南側(cè)為海拔高度約100 m的龍王山風(fēng)景區(qū),北方為盤(pán)城住宅區(qū)和農(nóng)田,東側(cè)為寧六路主干道及中國(guó)石化揚(yáng)子石油化工有限公司、南京鋼鐵集團(tuán)有限公司等大型企業(yè)。該地區(qū)處于西風(fēng)帶,屬于季風(fēng)性氣候區(qū),春夏季多為東南風(fēng),秋冬季盛行西北風(fēng)。采用PM2.5大流量采樣器(嶗應(yīng)2031型,青島嶗山應(yīng)用技術(shù)研究所)采集PM2.5樣品,采樣日期為2021年1月8日—2月28日,每次采樣時(shí)間為早上10:00至次日早上08:00,約22 h,共采集樣品26個(gè)。采樣濾膜為whatman石英濾膜,實(shí)驗(yàn)前需將濾膜放入450 ℃馬弗爐中烘烤4 h,以排除有機(jī)物和水分的干擾。濾膜采樣前后放入干燥器平衡24 h,采用1/100 000天平稱(chēng)重。采樣前后石英濾膜均用錫箔紙包裹并放于自封袋中,采集后的樣品放入干燥器內(nèi)保存直至化學(xué)分析。

    1.2 重金屬的化學(xué)形態(tài)提取

    采用逐級(jí)連續(xù)提取法(BCR)[21]將PM2.5中重金屬化學(xué)形態(tài)分為弱酸溶解態(tài)(F1)、可還原態(tài)(F2)、可氧化態(tài)(F3)、殘?jiān)鼞B(tài)(F4)。使用電感耦合等離子體原子發(fā)射光譜儀(ICP-OES,美國(guó)PE,Optima8000)測(cè)定重金屬濃度。為保證測(cè)量結(jié)果的準(zhǔn)確性,每組樣品設(shè)置3個(gè)空白樣品,得到重金屬元素的相對(duì)標(biāo)準(zhǔn)偏差小于10%。

    1.3 富集因子

    富集因子(FE)是評(píng)估不同環(huán)境中某元素富集程度的指標(biāo),其計(jì)算公式[15]為

    FE=(Ci,s/Cn,s)/(Ci,b/Cn,b)。

    (1)

    式(1)中,Ci,s和Cn,s分別為樣品中i元素與參比元素含量,mg·m-3;Ci,b和Cn,b分別為背景中i元素與參比元素含量,mg·m-3。選擇Al作為參比元素,各元素土壤背景值取自中國(guó)土壤背景值[22]。當(dāng)FE<1時(shí),認(rèn)為該元素主要來(lái)自自然源;當(dāng)1≤FE<10時(shí),認(rèn)為該元素來(lái)自自然源與人為源;當(dāng)FE>10時(shí),認(rèn)為該元素主要來(lái)自人為源,且FE越大,受到人為源的影響越強(qiáng)烈。

    1.4 正定矩陣因子模型(PMF)

    PMF是一種多變量因子分析方法,該方法將樣本數(shù)據(jù)矩陣分解為源成分譜矩陣、源貢獻(xiàn)矩陣,通過(guò)尋求最小目標(biāo)函數(shù)Q的值,確定污染源成分譜矩陣和污染源貢獻(xiàn)矩陣,獲得解析結(jié)果。計(jì)算公式[8]為

    (2)

    (3)

    U=L×5/6 (c≤L),

    (4)

    (5)

    式(2)~(5)中,Xij為樣本數(shù)據(jù)矩陣;Gik為源成分譜矩陣;Fkj為源貢獻(xiàn)矩陣;Eij為觀(guān)測(cè)值與樣本值的殘差矩陣;Uij為不確定度;n為樣本數(shù);m為物種數(shù);P為污染源數(shù);L為儀器檢測(cè)限;f表示測(cè)定結(jié)果不確定度的百分比,一般設(shè)為10%;c為元素濃度,μg·m-3。

    1.5 后向軌跡模型(HYSPLIT)

    以南京信息工程大學(xué)學(xué)科一號(hào)樓(32.20° N,118.71° E)為起點(diǎn),模擬采樣期間的48 h后向氣團(tuán)軌跡,并使用TrajStat軟件提供的Euclidean distance算法對(duì)氣團(tuán)軌跡進(jìn)行聚類(lèi)分析。該算法的原理是:根據(jù)氣流的移動(dòng)速度和方向?qū)Σ煌壽E以組內(nèi)差異和組間差異極化的原則進(jìn)行聚類(lèi),聚類(lèi)的數(shù)目采用總空間方差法(TSV)確定,即一開(kāi)始隨著聚類(lèi)數(shù)量的減少,TSV曲線(xiàn)僅呈小幅度的增加趨勢(shì),并且?guī)缀醣3植蛔?但到達(dá)某一時(shí)間點(diǎn)之后,TSV曲線(xiàn)會(huì)突然迅速上升,表明聚類(lèi)軌跡之間相差較大,因此把到達(dá)這一點(diǎn)所生成的聚類(lèi)軌跡數(shù)目作為最終的聚類(lèi)數(shù)量[15]。

    1.6 濃度權(quán)重軌跡模型(CWT)

    以HYSPLIT中得到的后向軌跡為基礎(chǔ),使用TrajStat軟件將研究區(qū)域(25°~53° N、97°~127° E)劃分為0.5°×0.5°的網(wǎng)格,利用濃度權(quán)重軌跡分析法(CWT)計(jì)算軌跡的平均權(quán)重濃度(TCW,ij),并引入經(jīng)驗(yàn)權(quán)重函數(shù)(Wij)來(lái)減小不確定性誤差,得到源區(qū)網(wǎng)格的平均權(quán)重濃度(TWCW,ij),確定其對(duì)南京北郊PM2.5中重金屬健康風(fēng)險(xiǎn)的貢獻(xiàn)。計(jì)算公式[15]為

    (6)

    (7)

    式(6)~(7)中,TWCW,ij為網(wǎng)格(i,j)的平均權(quán)重濃度;TCW,ij為軌跡的平均權(quán)重濃度;l為網(wǎng)格(i,j)內(nèi)第l條軌跡;M為網(wǎng)格(i,j)總氣流軌跡數(shù)目;Cl為軌跡l經(jīng)過(guò)網(wǎng)格(i,j)時(shí)對(duì)應(yīng)的健康風(fēng)險(xiǎn);τijl為軌跡l在網(wǎng)格(i,j)停留時(shí)間;nij為氣流軌跡通過(guò)網(wǎng)格(i,j)的端點(diǎn)數(shù)。

    1.7 健康風(fēng)險(xiǎn)評(píng)價(jià)

    采用EPA健康風(fēng)險(xiǎn)評(píng)價(jià)模型評(píng)估通過(guò)呼吸吸入途徑進(jìn)入人體的重金屬的健康風(fēng)險(xiǎn)。重金屬的健康風(fēng)險(xiǎn)分為非致癌風(fēng)險(xiǎn)與致癌風(fēng)險(xiǎn)。根據(jù)美國(guó)環(huán)境保護(hù)署綜合危險(xiǎn)度信息數(shù)據(jù)庫(kù)(IRIS)和國(guó)際癌癥研究機(jī)構(gòu)(IARC)的研究成果,Cr、Co、As、Cd為致癌元素。其計(jì)算公式[23]為

    Qi=Di/DR,i,

    (8)

    (9)

    Ri=Di×FS,i,

    (10)

    Di=Ci×(F1i+F2i+F3i)×R×F×D/(W×T)。

    (11)

    式(8)~(11)中,I為重金屬的綜合非致癌風(fēng)險(xiǎn);Qi為i元素的非致癌風(fēng)險(xiǎn);Ri為i元素的致癌風(fēng)險(xiǎn);Di為i元素的日均暴露劑量(終身日均暴露劑量),mg·kg-1·d-1;DR,i為i元素的參考劑量,mg·kg-1·d-1;FS,i為i元素的致癌斜率因子,kg·d·mg-1;Ci為i元素的濃度,mg·m-3,Cr有3價(jià)和6價(jià)2種價(jià)態(tài)存在,因USEPA及我國(guó)環(huán)境大氣標(biāo)準(zhǔn)均以大氣Cr6+作為標(biāo)準(zhǔn)限值依據(jù),故式(11)中Cr元素的濃度(CCr)以Cr總量的1/6參與計(jì)算[24];F1i、F2i、F3i分別為i元素弱酸溶解態(tài)、可還原態(tài)、可氧化態(tài)含量的百分?jǐn)?shù),%;R為呼吸速率,m3·d-1;F為暴露頻率,d·y-1;D為暴露年限,y;W為體重,kg;T為平均暴露時(shí)間,d。參數(shù)取值見(jiàn)表1~2[23]。

    表1 健康風(fēng)險(xiǎn)暴露參數(shù)數(shù)值

    表2 重金屬的參考劑量(DR)與致癌斜率因子(FS)

    當(dāng)Qi大于1時(shí),認(rèn)為i元素存在明顯的非致癌風(fēng)險(xiǎn),當(dāng)Qi小于1,認(rèn)為i元素非致癌風(fēng)險(xiǎn)可忽略或不存在非致癌風(fēng)險(xiǎn);當(dāng)I大于1時(shí),認(rèn)為重金屬存在明顯的非致癌風(fēng)險(xiǎn),當(dāng)I小于1,認(rèn)為重金屬非致癌風(fēng)險(xiǎn)可忽略或不存在非致癌風(fēng)險(xiǎn);當(dāng)Ri小于10-6時(shí),認(rèn)為i元素?zé)o致癌風(fēng)險(xiǎn),當(dāng)Ri在10-6~10-4之間,認(rèn)為i元素具有潛在致癌風(fēng)險(xiǎn),當(dāng)Ri大于10-4時(shí),認(rèn)為i元素存在明顯致癌風(fēng)險(xiǎn)[23]。

    1.8 氣象數(shù)據(jù)來(lái)源

    溫度、相對(duì)濕度及風(fēng)速數(shù)據(jù)來(lái)自于南京信息工程大學(xué)觀(guān)測(cè)基地自動(dòng)氣象站數(shù)據(jù),后向軌跡模型中的全球資料同化氣象數(shù)據(jù)來(lái)自于ftp:∥arlftp.arlhq.noaa.gov。

    2 結(jié)果與分析

    2.1 PM2.5及重金屬的污染水平

    由圖1可知,采樣期間溫度、相對(duì)濕度及風(fēng)速變化范圍分別為-1.0~13.0 ℃、34.0%~92.9%、0~3.6 m·s-1,變化范圍均較大,其平均值為5.9 ℃、62.0%、1.7 m·s-1,溫度、相對(duì)濕度及風(fēng)速均較低,表明采樣期間大氣擴(kuò)散度和稀釋度均較低,不利于污染物擴(kuò)散。采樣期間PM2.5日平均質(zhì)量濃度范圍為51.9~143.0 μg·m-3,平均值為(99.0±18.3)μg·m-3,是GB 3095—2012二級(jí)標(biāo)準(zhǔn)(75 μg·m-3)的1.3倍,采樣期間有92.3%的天數(shù)超過(guò)二級(jí)標(biāo)準(zhǔn)。

    圖1 采樣期間PM2.5和氣象參數(shù)時(shí)間序列

    由表3可知,采樣期間11種元素的總?cè)站|(zhì)量濃度為(1 999.0±307.2)ng·m-3,占PM2.5日均質(zhì)量濃度的2.0%。其中Al最高,為(786.4±243.2)ng·m-3,其他依次為Pb>V>Mn>Cu>Ba>Cr>Co>Sr>Cd>As,As和Cd日均濃度分別為(7.0±3.5)和(22.0±14.0)ng·m-3,是標(biāo)準(zhǔn)值的1.2和4.4倍。各元素富集因子順序?yàn)镃d>Pb>Cu>Co>V>Cr>As>Sr>Mn>Ba>Al,根據(jù)富集因子分類(lèi)標(biāo)準(zhǔn)[15],除參比元素Al外,其余元素的富集因子均大于10,表明這些元素均受到不同程度人為活動(dòng)的影響,其中Mn、Sr、Ba富集因子較小,表明其受到的人為活動(dòng)的影響較小,而As、Cr、V、Co、Cu、Pb、Cd富集因子均大于50,尤其是Cd和Pb富集因子遠(yuǎn)高于50,說(shuō)明采樣期間大部分重金屬元素受人類(lèi)活動(dòng)影響較大。變異系數(shù)(CV)可反映元素含量的離散程度,CV值越大,說(shuō)明元素分布越不均勻,越可能受人為來(lái)源的控制和影響。采樣期間11種重金屬總量的變異系數(shù)為15.5%,除Cr和Pb外,其余元素的變異系數(shù)均大于15.0%,根據(jù)變異系數(shù)分類(lèi)標(biāo)準(zhǔn)[25],重金屬總量及大部分元素均為中等變異及以上,表明采樣期間大部分重金屬元素受人類(lèi)活動(dòng)影響較大。

    表3 采樣期間PM2.5中金屬元素質(zhì)量濃度和富集因子

    2.2 重金屬化學(xué)形態(tài)和生物可利用性

    采用改進(jìn)的BCR逐級(jí)連續(xù)提取法分析PM2.5中重金屬的形態(tài)(圖2)。V、Cr、Mn、As、Sr、Cd、Ba主要存在形態(tài)為弱酸溶解態(tài)(F1),其F1質(zhì)量分?jǐn)?shù)分別為51.8%、36.4%、44.6%、72.4%、64.5%、60.2%、63.9%,除Cr和Mn外,其余元素的F1質(zhì)量分?jǐn)?shù)高于50%,說(shuō)明這7種元素遷移活性較強(qiáng),環(huán)境風(fēng)險(xiǎn)較高[27]。Co、Cu、Pb主要存在形態(tài)為可還原態(tài)(F2),其F2質(zhì)量分?jǐn)?shù)分別為58.1%、48.7%、44.0%,表明Co、Cu、Pb具有較強(qiáng)的潛在可遷移性,環(huán)境風(fēng)險(xiǎn)相對(duì)較高[28]。Al主要存在形態(tài)為殘?jiān)鼞B(tài)(F4),其F4質(zhì)量分?jǐn)?shù)為46.3%,接近50%,說(shuō)明Al遷移活性較差,環(huán)境風(fēng)險(xiǎn)相對(duì)較低[27]。

    圖2 PM2.5中重金屬的化學(xué)形態(tài)和生物可利用性

    重金屬在生物體內(nèi)的溶解和吸收很大程度上取決于其生物可利用性。通常情況下,弱酸溶解態(tài)被認(rèn)為是生物直接可利用形態(tài),可還原態(tài)、可氧化態(tài)被認(rèn)為是生物潛在可利用形態(tài),利用三者之和與重金屬總量的比值來(lái)表示其生物可利用性系數(shù)[6]。結(jié)果表明,采樣期間Cd(98.2%)的生物可利用形態(tài)含量最高,其次為Ba(97.8%)>Co(97.1%)>Mn(96.6)>As(95.3%)>Pb(93.8%)>Sr(91.0%)>Cr(90.9%)>Cu(90.4%)>V(87.2%)>Al(53.7%)。除Al外,其余元素的生物可利用態(tài)含量均高于85%,表明采樣期間大部分重金屬元素生物可利用性較高。

    2.3 來(lái)源解析

    將重金屬元素濃度數(shù)據(jù)輸入PMF模型,共解析得到3類(lèi)源,結(jié)果如圖3所示。因子1中Sr、Ba和Al占比最高,分別為85.0%、69.1%和34.7%。Sr主要來(lái)自于土壤揚(yáng)塵,Ba主要來(lái)源于道路揚(yáng)塵和建筑揚(yáng)塵[4],Al通常也被認(rèn)為是揚(yáng)塵源的標(biāo)志性元素,但是Al與Sr和Ba的相關(guān)性較差,且其CV>30%,考慮Al還有一部分來(lái)自人為排放源,而B(niǎo)a和Sr的EF值較低,受人為源影響相對(duì)較小,故因子1可歸結(jié)為揚(yáng)塵源,貢獻(xiàn)率為31.4%(圖3)。因子2中Al、V、Cr、Mn、Co、Cu、Pb為特征組分,其占比分別為56.3%、52.3%、87.4%、69.7%、56.6%、46.1%、51.0%。V是燃油的標(biāo)識(shí)性組分,Cr主要來(lái)源于工業(yè)排放和固體廢棄物的燃燒,Mn主要來(lái)自于機(jī)動(dòng)車(chē)汽油添加劑和鋼鐵冶煉,Co主要來(lái)源于金屬冶煉和土壤塵,Cu與機(jī)動(dòng)車(chē)剎車(chē)片的磨損有關(guān),Pb主要來(lái)自機(jī)動(dòng)車(chē)尾氣和工業(yè)排放[4,26]。Al與土壤揚(yáng)塵的標(biāo)志性元素相關(guān)性較差,推測(cè)其除了來(lái)自土壤揚(yáng)塵之外還有一部分來(lái)自于鋁合金冶煉等工業(yè)排放。Mn是汽車(chē)燃油中的重要添加劑,GB 17930—2013《車(chē)用汽油國(guó)家標(biāo)準(zhǔn)》規(guī)定,國(guó)IV汽車(chē)升級(jí)到國(guó)V,Mn的含量標(biāo)準(zhǔn)也下降至0.002 g·L-1,具體到Mn的暴露頻率為15.8,受人為影響較小,說(shuō)明標(biāo)準(zhǔn)的提高顯著降低了機(jī)動(dòng)車(chē)Mn的排放量。故因子2可定義為交通和工業(yè)混合源,貢獻(xiàn)率為50.9%(圖3)。因子3中As、Cd占比最高,為67.3%和88.2%。As是燃煤的標(biāo)識(shí)性元素,Cd來(lái)源于煤炭燃燒和冶金工業(yè)[4],因此認(rèn)為因子3代表燃煤源,貢獻(xiàn)率為17.7%(圖3)。

    2.4 健康風(fēng)險(xiǎn)評(píng)價(jià)

    2.4.1重金屬對(duì)人群的健康風(fēng)險(xiǎn)

    為探究PM2.5中重金屬對(duì)兒童和成人的健康影響,采用美國(guó)EPA健康風(fēng)險(xiǎn)評(píng)價(jià)模型計(jì)算其非致癌風(fēng)險(xiǎn)與致癌風(fēng)險(xiǎn),結(jié)果見(jiàn)表4。Cr、Co、Cu、As、Cd、Pb的非致癌風(fēng)險(xiǎn)均為兒童>成人,6種元素對(duì)兒童和成人的綜合非致癌風(fēng)險(xiǎn)分別為7.24×10-1和3.06×10-1,均小于1,表明采樣期間Cr、Co、Cu、As、Cd、Pb對(duì)兒童和成人的非致癌風(fēng)險(xiǎn)可忽略;重金屬的非致癌風(fēng)險(xiǎn)順序?yàn)镃o>Cr>Pb>As>Cd>Cu,其中Co(55.3%)的貢獻(xiàn)最高,表明采樣期間Co是非致癌風(fēng)險(xiǎn)的主要貢獻(xiàn)元素。Cr、Co、As、Cd的致癌風(fēng)險(xiǎn)均為成人>兒童,4種元素對(duì)兒童和成人的致癌風(fēng)險(xiǎn)均為10-6~10-4,表明采樣期間Cr、Co、As、Cd對(duì)兒童和成人有潛在的致癌風(fēng)險(xiǎn);不同金屬的致癌風(fēng)險(xiǎn)順序?yàn)镃r>Co>Cd>As,其中Cr(43.9%)的貢獻(xiàn)最高,表明采樣期間Cr是致癌風(fēng)險(xiǎn)的主要貢獻(xiàn)元素。

    表4 PM2.5中重金屬元素的非致癌風(fēng)險(xiǎn)與致癌風(fēng)險(xiǎn)

    2.4.2基于PMF的健康風(fēng)險(xiǎn)

    將PMF結(jié)果代入EPA健康風(fēng)險(xiǎn)評(píng)價(jià)模型,計(jì)算各類(lèi)排放源的健康風(fēng)險(xiǎn),結(jié)果見(jiàn)表5。3類(lèi)排放源對(duì)兒童和成人的非致癌風(fēng)險(xiǎn)均小于1,表明采樣期間3類(lèi)排放源對(duì)兒童和成人的非致癌風(fēng)險(xiǎn)可忽略;各排放源對(duì)兒童和成人的非致癌風(fēng)險(xiǎn)排序?yàn)榻煌ê凸I(yè)混合源(53.1%)>揚(yáng)塵源(25.8%)>燃煤源(21.0%),表明采樣期間交通和工業(yè)混合源是非致癌風(fēng)險(xiǎn)的主要貢獻(xiàn)源。3類(lèi)排放源對(duì)兒童和成人的致癌風(fēng)險(xiǎn)均為10-6~10-4,表明采樣期間3類(lèi)排放源對(duì)兒童和成人均有潛在的致癌風(fēng)險(xiǎn);3類(lèi)排放源對(duì)兒童和成人的致癌風(fēng)險(xiǎn)排序?yàn)榻煌ê凸I(yè)混合源(50.1%)>燃煤源(26.9%)>揚(yáng)塵源(23.0%),表明采樣期間交通和工業(yè)混合源是致癌風(fēng)險(xiǎn)的主要貢獻(xiàn)源。

    表5 不同排放源的非致癌風(fēng)險(xiǎn)與致癌風(fēng)險(xiǎn)

    2.4.3基于聚類(lèi)軌跡的健康風(fēng)險(xiǎn)

    為評(píng)估采樣期間區(qū)域傳輸對(duì)健康風(fēng)險(xiǎn)的影響,結(jié)合后向軌跡模型(HYSPLIT)計(jì)算聚類(lèi)軌跡及其健康風(fēng)險(xiǎn),結(jié)果見(jiàn)表6。由表6可知,4條軌跡的非致癌風(fēng)險(xiǎn)為兒童>成人,4條軌跡對(duì)兒童和成人的非致癌風(fēng)險(xiǎn)(7.2×10-1、7.2×10-1、8.1×10-1、7.0×10-1)均小于1,表明采樣期間4條軌跡對(duì)兒童和成人的非致癌風(fēng)險(xiǎn)可忽略;4條軌跡對(duì)兒童的非致癌風(fēng)險(xiǎn)排序?yàn)檐壽E3(27.3%)>軌跡2(24.6%)>軌跡1(24.2%)>軌跡4(23.9%),對(duì)成人的非致癌風(fēng)險(xiǎn)排序?yàn)檐壽E3(27.0%)>軌跡2=軌跡1(24.6%)>軌跡4(23.8%),表明采樣期間第3條軌跡是南京北郊非致癌風(fēng)險(xiǎn)的主要貢獻(xiàn)軌跡。采樣期間4條軌跡的致癌風(fēng)險(xiǎn)為成人>兒童,4條軌跡對(duì)兒童的致癌風(fēng)險(xiǎn)為10-6~10-4,對(duì)成人的致癌風(fēng)險(xiǎn)均略高于10-4,表明采樣期間4條軌跡對(duì)兒童有潛在致癌風(fēng)險(xiǎn),對(duì)成人有明顯的致癌風(fēng)險(xiǎn);4條軌跡對(duì)兒童和成人的致癌風(fēng)險(xiǎn)排序?yàn)檐壽E3>軌跡2>軌跡1>軌跡4,表明采樣期間第3條軌跡是南京北郊致癌風(fēng)險(xiǎn)的主要貢獻(xiàn)軌跡。

    表6 采樣期間聚類(lèi)軌跡的非致癌風(fēng)險(xiǎn)與致癌風(fēng)險(xiǎn)

    2.4.4基于CWT的健康風(fēng)險(xiǎn)

    為探究采樣期間潛在源區(qū)對(duì)南京北郊健康風(fēng)險(xiǎn)的具體貢獻(xiàn)大小,將健康風(fēng)險(xiǎn)數(shù)據(jù)代入CWT模型,結(jié)果見(jiàn)圖4。各潛在源區(qū)對(duì)兒童(<0.88)和成人(<0.36)的非致癌風(fēng)險(xiǎn)均小于1,表明采樣期間各潛在源區(qū)的非致癌風(fēng)險(xiǎn)可忽略;各潛在源區(qū)的非致癌風(fēng)險(xiǎn)排序?yàn)榘不帐|部(0.82~0.88、0.34~0.36)>河南省信陽(yáng)市(0.76~0.82、0.32~0.34)>江蘇省西部(0.70~0.76、0.30~0.32)>其他區(qū)域(<0.70、<0.30),表明采樣期間安徽省東部是南京北郊非致癌風(fēng)險(xiǎn)的主要貢獻(xiàn)源區(qū)。各潛在源區(qū)對(duì)兒童的致癌風(fēng)險(xiǎn)均為10-6~10-4,河南省信陽(yáng)市、安徽省東部、江蘇省西部對(duì)成人的致癌風(fēng)險(xiǎn)高于10-4,其他區(qū)域?qū)Τ扇说闹掳╋L(fēng)險(xiǎn)為10-6~10-4,表明采樣期間各潛在源區(qū)對(duì)兒童有潛在的致癌風(fēng)險(xiǎn),河南省信陽(yáng)市、安徽省東部、江蘇省西部對(duì)成人有明顯的致癌風(fēng)險(xiǎn)。潛在源區(qū)對(duì)兒童和成人的致癌風(fēng)險(xiǎn)排序?yàn)楹幽鲜⌒抨?yáng)市(6.40×10-5~7.20×10-5、1.36×10-4~1.54×10-4)>安徽省東部(5.60×10-5~6.40×10-5、1.18×10-4~1.36×10-4)>江蘇省西部(4.80×10-5~5.60×10-5、1.00×10-4~1.18×10-4)>其他區(qū)域(<4.80×10-5、<1.00×10-4),表明采樣期間河南省信陽(yáng)市是南京北郊致癌風(fēng)險(xiǎn)的主要貢獻(xiàn)源區(qū)。

    3 討論

    3.1 PM2.5中重金屬濃度、形態(tài)及來(lái)源

    從重金屬的濃度水平看,As、Cd濃度均高于GB 3095—2012標(biāo)準(zhǔn)限值,與北京、武漢、唐山、鄭州、新鄉(xiāng)[29-32]等國(guó)內(nèi)城市相比,V、Cr、Mn、Co、Cu、Cd、Pb濃度均偏高,結(jié)合其富集特征、形態(tài)、生物可利用性及來(lái)源解析結(jié)果,發(fā)現(xiàn)這些元素大部分是交通和工業(yè)混合源的特征組分,且其富集性顯著,其主要存在形態(tài)均為弱酸溶解態(tài)(F1)或可還原態(tài)(F2),生物可利用性較高,可能對(duì)當(dāng)?shù)匚廴矩暙I(xiàn)較為明顯,故在大氣重金屬風(fēng)險(xiǎn)防控中需重點(diǎn)關(guān)注。

    3.2 PM2.5中重金屬的健康風(fēng)險(xiǎn)

    從健康風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果看,PM2.5中Cr和Co、交通和工業(yè)混合源及信陽(yáng)市和安徽省東部的區(qū)域傳輸對(duì)南京北郊健康風(fēng)險(xiǎn)的影響較大,其原因一方面可能是采樣點(diǎn)地理位置及氣象條件的影響。采樣點(diǎn)緊鄰安徽省,且周?chē)写笮突@區(qū),南京是典型的季風(fēng)氣候區(qū),冬季受到西伯利亞寒流影響較多,盛行西北風(fēng),而信陽(yáng)市、安徽省東部均位于南京偏西方向,故其產(chǎn)生的污染較易擴(kuò)散到南京。另一方面可能與其城市的工業(yè)化水平有關(guān)。信陽(yáng)市明港鎮(zhèn)產(chǎn)業(yè)集聚區(qū),有鋼鐵廠(chǎng)信鋼公司、聯(lián)強(qiáng)實(shí)業(yè)、中聯(lián)制管廠(chǎng)、微粉廠(chǎng)、河南杭氧公司等一系列工業(yè)、重工業(yè)廠(chǎng)區(qū),安徽省的銅陵市、馬鞍山市是安徽工業(yè)化水平最高的2個(gè)城市,其有色金屬冶煉業(yè)及鋼鐵工業(yè)所占比重均超過(guò)60%,合肥市、蕪湖市等工業(yè)化水平也較高,導(dǎo)致其輸送的工業(yè)源較多。從以往的研究看,大氣重金屬污染治理中多強(qiáng)調(diào)對(duì)本地源的管控,忽略了關(guān)注區(qū)域傳輸?shù)慕】碉L(fēng)險(xiǎn),但從研究結(jié)果看,區(qū)域傳輸?shù)慕】碉L(fēng)險(xiǎn)要高于本地源的貢獻(xiàn),說(shuō)明有時(shí)區(qū)域傳輸?shù)慕】碉L(fēng)險(xiǎn)不容忽視。因此,為保障南京北郊居民安全和減少健康風(fēng)險(xiǎn),不僅要加強(qiáng)對(duì)本地工業(yè)排放重金屬含量的評(píng)估及控制,尤其是Co和Cr,也要密切關(guān)注周邊區(qū)域?qū)Ρ镜氐挠绊憽?/p>

    4 結(jié)論

    (1)采樣期間PM2.5日均質(zhì)量濃度為51.9~143.0 μg·m-3,平均值為(99.0±18.3)μg·m-3。11種元素總?cè)站|(zhì)量濃度為(1 999.0±307.2)ng·m-3,As、Cd日均質(zhì)量濃度分別為(7.0±3.5)和(22.0±14.0)ng·m-3,是標(biāo)準(zhǔn)值的1.2、4.4倍。

    (2)BCR提取結(jié)果表明,大部分重金屬主要存在于弱酸溶解態(tài)(F1),Co、Cu、Pb主要存在于可還原態(tài)(F2),Al主要存在于殘?jiān)鼞B(tài)(F4)。人為排放的Co、Cd、Ba、Mn、As、Pb、Cr、Cu和V等元素富集性顯著,且生物可利用性較高,應(yīng)引起關(guān)注。

    (3)PMF解析結(jié)果表明,南京北郊PM2.5中重金屬主要有三大來(lái)源,分別為揚(yáng)塵源(31.4%)、燃煤源(17.7%)、交通和工業(yè)混合源(50.9%)。

    (4)健康風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果表明,重金屬的非致癌風(fēng)險(xiǎn)可忽略,Cr、Co、As、Cd對(duì)兒童和成人有潛在的致癌風(fēng)險(xiǎn),Cr和Co是主要貢獻(xiàn)元素。結(jié)合健康風(fēng)險(xiǎn)評(píng)價(jià)與正定矩陣因子分解模型、后向軌跡模型、濃度權(quán)重軌跡模型發(fā)現(xiàn),交通和工業(yè)混合源及周邊區(qū)域傳輸對(duì)南京北郊健康風(fēng)險(xiǎn)影響不容忽視。

    猜你喜歡
    成人軌跡重金屬
    軌跡
    軌跡
    重金屬對(duì)膨潤(rùn)土膨脹性的影響
    成人不自在
    測(cè)定不同產(chǎn)地寬筋藤中5種重金屬
    中成藥(2018年8期)2018-08-29 01:28:16
    軌跡
    進(jìn)化的軌跡(一)——進(jìn)化,無(wú)盡的適應(yīng)
    ICP-AES、ICP-MS測(cè)定水中重金屬的對(duì)比研究
    再生水回灌中DOM對(duì)重金屬遷移與保留問(wèn)題研究
    Un rite de passage
    中文字幕精品免费在线观看视频 | 国产精品国产三级专区第一集| 久久久国产欧美日韩av| 久久久久精品性色| 亚洲精品av麻豆狂野| 国产国拍精品亚洲av在线观看| 成人亚洲欧美一区二区av| 国产高清不卡午夜福利| 午夜福利在线观看免费完整高清在| 国产亚洲av片在线观看秒播厂| 内地一区二区视频在线| 欧美日韩国产mv在线观看视频| 啦啦啦视频在线资源免费观看| 久久99热这里只频精品6学生| 少妇的丰满在线观看| 免费不卡的大黄色大毛片视频在线观看| 一区二区av电影网| 一本色道久久久久久精品综合| 五月玫瑰六月丁香| 日韩中文字幕视频在线看片| 欧美少妇被猛烈插入视频| 欧美性感艳星| 97在线视频观看| 精品久久蜜臀av无| 日韩成人av中文字幕在线观看| 又黄又粗又硬又大视频| 久久人人爽人人片av| 亚洲图色成人| 久久亚洲国产成人精品v| 久久久久久久久久久免费av| 国产高清不卡午夜福利| 国产有黄有色有爽视频| 国产毛片在线视频| 亚洲精品一二三| 亚洲,欧美精品.| 夫妻性生交免费视频一级片| 人人妻人人澡人人爽人人夜夜| 肉色欧美久久久久久久蜜桃| 寂寞人妻少妇视频99o| 亚洲一区二区三区欧美精品| 91午夜精品亚洲一区二区三区| 欧美人与性动交α欧美精品济南到 | 亚洲精品久久午夜乱码| 国产乱来视频区| 考比视频在线观看| 久久99蜜桃精品久久| 国产高清国产精品国产三级| 99热这里只有是精品在线观看| 少妇被粗大的猛进出69影院 | 亚洲国产最新在线播放| av黄色大香蕉| 男女下面插进去视频免费观看 | 人妻一区二区av| 国产精品国产三级专区第一集| 一级a做视频免费观看| 咕卡用的链子| 少妇人妻精品综合一区二区| 国产精品久久久久久精品古装| 亚洲图色成人| 成人午夜精彩视频在线观看| 亚洲国产精品999| 黄色配什么色好看| 国产精品久久久久久久电影| 欧美精品人与动牲交sv欧美| 99热网站在线观看| 人人澡人人妻人| 色婷婷久久久亚洲欧美| 熟妇人妻不卡中文字幕| 亚洲图色成人| 激情五月婷婷亚洲| 国产成人av激情在线播放| 成人漫画全彩无遮挡| 在线观看免费高清a一片| 又黄又粗又硬又大视频| 国产精品久久久久久精品古装| 亚洲成av片中文字幕在线观看 | freevideosex欧美| 高清不卡的av网站| 国产欧美日韩一区二区三区在线| 女人久久www免费人成看片| 亚洲美女视频黄频| 国产在视频线精品| 男女无遮挡免费网站观看| 精品福利永久在线观看| 成人二区视频| 三级国产精品片| 亚洲av综合色区一区| 韩国精品一区二区三区 | 日本午夜av视频| 国产一区二区在线观看av| 九色成人免费人妻av| 丰满迷人的少妇在线观看| 亚洲婷婷狠狠爱综合网| 亚洲欧美中文字幕日韩二区| 午夜精品国产一区二区电影| tube8黄色片| 赤兔流量卡办理| 久久鲁丝午夜福利片| 国产精品久久久久久久电影| 国产熟女午夜一区二区三区| 有码 亚洲区| 在线 av 中文字幕| 欧美激情 高清一区二区三区| 最新的欧美精品一区二区| 精品久久蜜臀av无| 搡老乐熟女国产| 久久国产亚洲av麻豆专区| 2021少妇久久久久久久久久久| 97精品久久久久久久久久精品| 女的被弄到高潮叫床怎么办| 又黄又爽又刺激的免费视频.| 18禁国产床啪视频网站| 亚洲成人av在线免费| 国产黄色免费在线视频| 极品人妻少妇av视频| 青春草亚洲视频在线观看| av卡一久久| 亚洲精品av麻豆狂野| 在线免费观看不下载黄p国产| 亚洲伊人色综图| 大香蕉久久网| 七月丁香在线播放| 国产乱来视频区| 在线天堂中文资源库| 国产成人91sexporn| 我要看黄色一级片免费的| 精品视频人人做人人爽| 久久99精品国语久久久| 久久午夜福利片| 日韩视频在线欧美| av免费观看日本| 亚洲色图综合在线观看| 久久久精品区二区三区| 亚洲精华国产精华液的使用体验| 国产有黄有色有爽视频| av片东京热男人的天堂| videossex国产| 欧美成人午夜免费资源| 久久久国产欧美日韩av| 亚洲人成77777在线视频| 亚洲国产av新网站| 天美传媒精品一区二区| 国产精品一二三区在线看| 精品人妻在线不人妻| 少妇高潮的动态图| 26uuu在线亚洲综合色| 男男h啪啪无遮挡| 黑人巨大精品欧美一区二区蜜桃 | 国产亚洲午夜精品一区二区久久| 欧美 亚洲 国产 日韩一| 国产在线一区二区三区精| 国产伦理片在线播放av一区| 国产熟女欧美一区二区| 只有这里有精品99| 777米奇影视久久| 国语对白做爰xxxⅹ性视频网站| 亚洲人与动物交配视频| 国产精品麻豆人妻色哟哟久久| 欧美日本中文国产一区发布| 男女啪啪激烈高潮av片| 国产成人免费观看mmmm| 日韩成人av中文字幕在线观看| 狂野欧美激情性bbbbbb| 亚洲国产毛片av蜜桃av| 人人妻人人澡人人爽人人夜夜| 国产一区二区三区av在线| 成人影院久久| 日韩一区二区三区影片| 秋霞伦理黄片| 亚洲伊人久久精品综合| 亚洲欧美中文字幕日韩二区| 亚洲精品美女久久av网站| 成人国产av品久久久| 黑人欧美特级aaaaaa片| 国产日韩欧美在线精品| 日韩一本色道免费dvd| 日韩av免费高清视频| 在线 av 中文字幕| 国产免费又黄又爽又色| 永久免费av网站大全| 久久人人爽人人爽人人片va| 精品少妇内射三级| 欧美日韩亚洲高清精品| 色5月婷婷丁香| 男女下面插进去视频免费观看 | 在线亚洲精品国产二区图片欧美| 日韩一区二区视频免费看| 狠狠婷婷综合久久久久久88av| 五月伊人婷婷丁香| 国国产精品蜜臀av免费| 2021少妇久久久久久久久久久| 纯流量卡能插随身wifi吗| 亚洲人成网站在线观看播放| 亚洲人成77777在线视频| 少妇被粗大猛烈的视频| 国产精品欧美亚洲77777| 高清欧美精品videossex| 丁香六月天网| 婷婷色综合大香蕉| 一二三四在线观看免费中文在 | 中文字幕最新亚洲高清| av在线app专区| 天天操日日干夜夜撸| 黄色怎么调成土黄色| 欧美精品高潮呻吟av久久| 亚洲精品,欧美精品| 尾随美女入室| 自拍欧美九色日韩亚洲蝌蚪91| 蜜桃国产av成人99| 日韩av不卡免费在线播放| 亚洲欧洲精品一区二区精品久久久 | 久久精品国产鲁丝片午夜精品| kizo精华| 国产一区二区三区av在线| 久久久久人妻精品一区果冻| 亚洲精品第二区| 欧美精品高潮呻吟av久久| 制服丝袜香蕉在线| 黄片无遮挡物在线观看| 国产综合精华液| 一二三四中文在线观看免费高清| 国产无遮挡羞羞视频在线观看| 少妇熟女欧美另类| 午夜视频国产福利| 亚洲精品美女久久av网站| 中文字幕人妻熟女乱码| xxxhd国产人妻xxx| 亚洲国产精品成人久久小说| 尾随美女入室| 最黄视频免费看| 亚洲婷婷狠狠爱综合网| 伊人久久国产一区二区| 国产精品一区二区在线观看99| 亚洲av.av天堂| 午夜福利视频精品| 亚洲国产欧美日韩在线播放| 90打野战视频偷拍视频| 日韩精品有码人妻一区| 黄片播放在线免费| 免费黄网站久久成人精品| 在线看a的网站| 亚洲成国产人片在线观看| 少妇人妻精品综合一区二区| 亚洲久久久国产精品| 亚洲精品自拍成人| 黄片播放在线免费| 亚洲av综合色区一区| 女性生殖器流出的白浆| 日韩制服骚丝袜av| 咕卡用的链子| 久久久久久久大尺度免费视频| 午夜福利,免费看| 亚洲综合色惰| 色婷婷久久久亚洲欧美| 亚洲精品久久成人aⅴ小说| 亚洲av电影在线进入| 国产亚洲欧美精品永久| 春色校园在线视频观看| 国产亚洲一区二区精品| 久久精品熟女亚洲av麻豆精品| 五月伊人婷婷丁香| 大片免费播放器 马上看| 日本91视频免费播放| 香蕉丝袜av| 欧美国产精品一级二级三级| 亚洲国产成人一精品久久久| 国产伦理片在线播放av一区| 亚洲天堂av无毛| 99久久精品国产国产毛片| 成人午夜精彩视频在线观看| 亚洲av中文av极速乱| 日日爽夜夜爽网站| 久久久国产欧美日韩av| 欧美日韩成人在线一区二区| 性色avwww在线观看| 国产精品不卡视频一区二区| 香蕉国产在线看| 在线观看免费高清a一片| 国产色婷婷99| 亚洲美女黄色视频免费看| 两个人免费观看高清视频| 人成视频在线观看免费观看| 国产精品国产三级国产av玫瑰| 色婷婷久久久亚洲欧美| 日韩欧美一区视频在线观看| 日产精品乱码卡一卡2卡三| 久久人人爽av亚洲精品天堂| 久久人人97超碰香蕉20202| 香蕉丝袜av| av天堂久久9| 欧美日韩亚洲高清精品| 免费女性裸体啪啪无遮挡网站| 多毛熟女@视频| 久久久久久人人人人人| 亚洲第一区二区三区不卡| 亚洲精品国产色婷婷电影| 久热久热在线精品观看| 你懂的网址亚洲精品在线观看| 色吧在线观看| 老司机亚洲免费影院| 中文字幕免费在线视频6| 国产毛片在线视频| 一区二区av电影网| www日本在线高清视频| 国产探花极品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 国产av国产精品国产| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 91精品伊人久久大香线蕉| 日韩欧美一区视频在线观看| 26uuu在线亚洲综合色| 王馨瑶露胸无遮挡在线观看| 国产色爽女视频免费观看| 色94色欧美一区二区| 一级片'在线观看视频| 在线观看美女被高潮喷水网站| 久久 成人 亚洲| 亚洲第一av免费看| 精品熟女少妇av免费看| 多毛熟女@视频| 欧美成人精品欧美一级黄| 精品少妇黑人巨大在线播放| 日本午夜av视频| 国产亚洲精品久久久com| 亚洲国产精品999| 亚洲欧美成人综合另类久久久| 99久久精品国产国产毛片| 亚洲三级黄色毛片| 在线观看免费日韩欧美大片| 国精品久久久久久国模美| 中文字幕免费在线视频6| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区| 亚洲图色成人| 久久久亚洲精品成人影院| 黄色毛片三级朝国网站| 国产日韩欧美在线精品| xxx大片免费视频| 老女人水多毛片| 五月伊人婷婷丁香| 国产欧美亚洲国产| 亚洲精品久久午夜乱码| 在线天堂中文资源库| 午夜免费鲁丝| 黄色一级大片看看| 精品少妇内射三级| 人体艺术视频欧美日本| 我要看黄色一级片免费的| 99国产综合亚洲精品| 好男人视频免费观看在线| 伦理电影大哥的女人| 国产av国产精品国产| 欧美日本中文国产一区发布| 久久国产亚洲av麻豆专区| 九草在线视频观看| 丝袜美足系列| 最近中文字幕2019免费版| videossex国产| 老司机影院毛片| 王馨瑶露胸无遮挡在线观看| 香蕉丝袜av| 亚洲综合精品二区| 亚洲久久久国产精品| 一级片'在线观看视频| 中文字幕免费在线视频6| 亚洲精品视频女| 春色校园在线视频观看| 日本黄色日本黄色录像| 18禁裸乳无遮挡动漫免费视频| 人妻系列 视频| 最近中文字幕2019免费版| 久久热在线av| 亚洲欧美色中文字幕在线| 免费观看a级毛片全部| 性色avwww在线观看| 18禁国产床啪视频网站| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品古装| 成人手机av| 国产精品久久久久久精品古装| 制服诱惑二区| 内地一区二区视频在线| 十分钟在线观看高清视频www| 日本91视频免费播放| 18禁裸乳无遮挡动漫免费视频| 亚洲国产毛片av蜜桃av| 一级a做视频免费观看| 国产av一区二区精品久久| 夫妻午夜视频| 亚洲,欧美精品.| 欧美日韩视频高清一区二区三区二| 欧美+日韩+精品| 日韩成人av中文字幕在线观看| 不卡视频在线观看欧美| 一二三四中文在线观看免费高清| 国产av精品麻豆| 亚洲国产最新在线播放| 国产福利在线免费观看视频| 亚洲五月色婷婷综合| 国产黄色免费在线视频| 国产一级毛片在线| 国产一区二区在线观看av| 精品亚洲成a人片在线观看| 少妇的丰满在线观看| 久久狼人影院| 成人18禁高潮啪啪吃奶动态图| 亚洲,欧美,日韩| 亚洲 欧美一区二区三区| 黄色一级大片看看| 国语对白做爰xxxⅹ性视频网站| 免费观看在线日韩| 一二三四在线观看免费中文在 | 欧美国产精品一级二级三级| 国产毛片在线视频| 午夜福利视频精品| 免费观看在线日韩| 欧美xxxx性猛交bbbb| 亚洲精品久久久久久婷婷小说| 国产69精品久久久久777片| 婷婷色综合大香蕉| 极品少妇高潮喷水抽搐| 九色成人免费人妻av| 国产成人午夜福利电影在线观看| 观看美女的网站| 亚洲中文av在线| 人妻人人澡人人爽人人| 精品人妻偷拍中文字幕| 国产精品久久久久久精品古装| 日韩伦理黄色片| 精品人妻一区二区三区麻豆| 国产男女内射视频| 亚洲,欧美,日韩| 卡戴珊不雅视频在线播放| 熟女电影av网| 国产精品99久久99久久久不卡 | 满18在线观看网站| a级毛色黄片| 最黄视频免费看| av在线app专区| 午夜激情久久久久久久| 蜜桃在线观看..| 啦啦啦中文免费视频观看日本| 国产精品一区二区在线观看99| 高清欧美精品videossex| 久久久久久久久久久久大奶| 亚洲精品中文字幕在线视频| 国产一区二区在线观看av| 免费少妇av软件| 国产亚洲最大av| 免费看光身美女| 嫩草影院入口| 中文天堂在线官网| 美女视频免费永久观看网站| av免费在线看不卡| 少妇的逼水好多| 久久99热这里只频精品6学生| 最近中文字幕2019免费版| 一级片免费观看大全| 妹子高潮喷水视频| 免费不卡黄色视频| 两个人免费观看高清视频| 午夜免费鲁丝| 国内毛片毛片毛片毛片毛片| 亚洲专区国产一区二区| 一级作爱视频免费观看| 国产在线观看jvid| 国产日韩一区二区三区精品不卡| 国产一区二区激情短视频| 国产在视频线精品| 18在线观看网站| 国产高清视频在线播放一区| 老司机午夜福利在线观看视频| av视频免费观看在线观看| 脱女人内裤的视频| 精品免费久久久久久久清纯 | 精品国产一区二区久久| 精品视频人人做人人爽| 制服人妻中文乱码| 成人免费观看视频高清| 在线十欧美十亚洲十日本专区| 亚洲精品自拍成人| 欧美+亚洲+日韩+国产| 老司机福利观看| 亚洲av美国av| 午夜福利影视在线免费观看| 这个男人来自地球电影免费观看| 三上悠亚av全集在线观看| 久久草成人影院| 欧美老熟妇乱子伦牲交| 午夜91福利影院| 中亚洲国语对白在线视频| 国产高清视频在线播放一区| 777米奇影视久久| 国产91精品成人一区二区三区| 一级毛片高清免费大全| 电影成人av| 美女视频免费永久观看网站| 久久久久视频综合| 精品乱码久久久久久99久播| 咕卡用的链子| 最新的欧美精品一区二区| 日韩制服丝袜自拍偷拍| 夫妻午夜视频| 亚洲av日韩精品久久久久久密| 女性被躁到高潮视频| tocl精华| 天天躁狠狠躁夜夜躁狠狠躁| 国产在线精品亚洲第一网站| 国产有黄有色有爽视频| 亚洲熟女毛片儿| 很黄的视频免费| av视频免费观看在线观看| 黄色毛片三级朝国网站| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利免费观看在线| 老司机在亚洲福利影院| 久99久视频精品免费| 精品一区二区三卡| 人妻丰满熟妇av一区二区三区 | 一区二区三区国产精品乱码| 欧美精品人与动牲交sv欧美| 999久久久国产精品视频| av网站在线播放免费| 捣出白浆h1v1| 热99国产精品久久久久久7| 成年人黄色毛片网站| 在线视频色国产色| 午夜老司机福利片| 久久精品aⅴ一区二区三区四区| 国产精品电影一区二区三区 | 欧美日韩黄片免| 在线av久久热| 久久久精品免费免费高清| 免费观看精品视频网站| 亚洲av第一区精品v没综合| 国产精品久久久久久精品古装| 久久久久视频综合| 天堂√8在线中文| 美女福利国产在线| 免费观看人在逋| 下体分泌物呈黄色| 国产成人一区二区三区免费视频网站| 少妇被粗大的猛进出69影院| 91在线观看av| 大香蕉久久成人网| 色在线成人网| 精品国产乱子伦一区二区三区| 国产精品自产拍在线观看55亚洲 | 一进一出抽搐动态| 另类亚洲欧美激情| 91九色精品人成在线观看| 这个男人来自地球电影免费观看| 午夜影院日韩av| 黄色a级毛片大全视频| 久久中文字幕人妻熟女| 亚洲三区欧美一区| 老司机午夜福利在线观看视频| 捣出白浆h1v1| 午夜激情av网站| 久久这里只有精品19| 久久精品国产清高在天天线| 亚洲专区字幕在线| 午夜精品久久久久久毛片777| 亚洲片人在线观看| 精品国产国语对白av| av超薄肉色丝袜交足视频| 99国产极品粉嫩在线观看| 制服人妻中文乱码| 男女之事视频高清在线观看| 一进一出抽搐gif免费好疼 | 建设人人有责人人尽责人人享有的| 夜夜夜夜夜久久久久| 一二三四社区在线视频社区8| 老熟妇乱子伦视频在线观看| 国产麻豆69| 91麻豆精品激情在线观看国产 | 午夜福利欧美成人| 亚洲中文日韩欧美视频| 大香蕉久久成人网| 国产又色又爽无遮挡免费看| 精品高清国产在线一区| 桃红色精品国产亚洲av| 精品一区二区三卡| 亚洲中文日韩欧美视频| 国产精华一区二区三区| 久久狼人影院| netflix在线观看网站| 午夜精品在线福利| 欧美黑人精品巨大| 在线观看舔阴道视频| 后天国语完整版免费观看| 国产又色又爽无遮挡免费看| 王馨瑶露胸无遮挡在线观看| 香蕉国产在线看| 国产主播在线观看一区二区| 天天躁日日躁夜夜躁夜夜| 成人黄色视频免费在线看| 一区二区日韩欧美中文字幕| 国产亚洲精品一区二区www | 国产一区在线观看成人免费| 国产成人精品久久二区二区免费| 一区二区三区精品91| 啦啦啦免费观看视频1| 在线国产一区二区在线| 香蕉国产在线看| 国产av一区二区精品久久| 中文字幕精品免费在线观看视频| 99热国产这里只有精品6| 91麻豆精品激情在线观看国产 | 欧美黑人精品巨大| 国产成+人综合+亚洲专区| 日本一区二区免费在线视频| 欧美老熟妇乱子伦牲交| 亚洲午夜精品一区,二区,三区|