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

    北京東南郊灌區(qū)多環(huán)芳烴污染風險與人體健康風險評估

    2017-10-11 01:34:39黃冠華黃權(quán)中劉洪祿
    農(nóng)業(yè)機械學(xué)報 2017年9期
    關(guān)鍵詞:表層總量籽粒

    李 艷 顧 華 黃冠華 黃權(quán)中 劉洪祿

    (1.北京市水科學(xué)技術(shù)研究院, 北京 100048; 2.北京市非常規(guī)水資源開發(fā)利用與節(jié)水工程技術(shù)研究中心, 北京 100048;3.中國農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院, 北京 100083)

    北京東南郊灌區(qū)多環(huán)芳烴污染風險與人體健康風險評估

    李 艷1,2顧 華1,2黃冠華3黃權(quán)中3劉洪祿1,2

    (1.北京市水科學(xué)技術(shù)研究院, 北京 100048; 2.北京市非常規(guī)水資源開發(fā)利用與節(jié)水工程技術(shù)研究中心, 北京 100048;3.中國農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院, 北京 100083)

    為明確北京東南郊灌區(qū)表層土壤和農(nóng)作物多環(huán)芳烴(PAHs)含量及污染水平,2015年在該灌區(qū)采集了31個土壤樣品和38個作物樣品,采用氣象色譜-質(zhì)譜儀分析了樣品中16種優(yōu)先控制PAHs含量。研究結(jié)果表明灌區(qū)表層土壤16種PAHs總量為113.5~449.8 μg/kg(均值為258.5 μg/kg),與荷蘭土壤環(huán)境質(zhì)量標準相比,本研究表層土壤中相應(yīng)10種PAHs的總量為83~307 μg/kg,土壤10種PAHs的苯并(a)芘毒性當量濃度均值為19.0 μg/kg,灌區(qū)土壤總體上受到PAHs輕度污染,但土壤PAHs含量遠低于荷蘭土壤環(huán)境質(zhì)量標準中的干預(yù)值,土壤PAHs潛在生態(tài)風險不顯著。小麥籽粒、玉米籽粒和蔬菜PAHs含量分別為368.7~389.9、87.7~113.5、51.8~291.8 μg/kg。糧食作物籽粒PAHs主要以苯并(g,h,i)苝為主,蔬菜PAHs以菲含量最高。糧食作物籽粒苯并(a)芘含量均低于食品安全國家標準限值。燃燒源是土壤和糧食作物PAHs主要污染源,而蔬菜PAHs主要污染源為燃燒源和石油源混合來源。成人和兒童PAHs致癌風險分別為4.02×10-5和1.76×10-5,非致癌指數(shù)分別為2.72×10-2和4.78×10-2,均低于US EPA相應(yīng)標準參考值。人體健康風險主要暴露途徑和介質(zhì)是口-作物,對致癌風險貢獻最大的組分為苯并(a)芘和二苯并(a,h)蒽,對非致癌風險貢獻最大的組分是苯并(g,h,i)苝和菲。

    北京東南郊灌區(qū); 多環(huán)芳烴; 土壤; 冬小麥; 夏玉米; 蔬菜; 風險評估

    引言

    多環(huán)芳烴(PAHs)作為有機污染物的一類,是較早被發(fā)現(xiàn)的致癌物質(zhì)[1],被美國環(huán)保局和我國列為環(huán)境中優(yōu)先控制的污染物[2]。PAHs主要來源于石油精煉、化石和煤炭的不完全燃燒,土壤是其主要儲庫,土壤中的PAHs主要來源于干/濕沉降、污泥農(nóng)用和污水灌溉等途徑[3-4]。土壤或大氣中的PAHs能被植物吸收積累,進而通過食物鏈危害人體健康[5]。因此,PAHs在土壤-作物系統(tǒng)中的污染效應(yīng)已逐漸成為環(huán)境科學(xué)等領(lǐng)域的研究熱點。

    國內(nèi)外一些學(xué)者調(diào)查研究了工業(yè)區(qū)(煤礦、化工、金屬冶煉等)、城區(qū)(公園、綠化等)、交通干線區(qū)以及郊區(qū)(城市周邊、遠郊農(nóng)業(yè)用地)等不同環(huán)境功能區(qū)表層土壤和作物PAHs污染情況[6-15],得出工業(yè)區(qū)、城區(qū)、交通區(qū)表層土壤PAHs含量通常高于郊區(qū),這3個區(qū)域大部分表層土壤受到了PAHs的污染,且部分區(qū)域土壤存在潛在生態(tài)風險;城市周邊農(nóng)田部分土壤受到PAHs污染,遠郊農(nóng)田土壤總體上未受到PAHs污染,農(nóng)產(chǎn)品PAHs含量低于歐盟與國內(nèi)相應(yīng)標準;城區(qū)、交通區(qū)及工業(yè)區(qū)表層土壤PAHs主要來源于化石燃料燃燒,郊區(qū)農(nóng)田土壤PAHs主要來源為燃燒源和石油源混合來源。另一些學(xué)者對國內(nèi)外污灌區(qū)土壤PAHs污染水平及健康風險進行了調(diào)查,得出大部分污灌區(qū)表層土壤受到PAHs的污染,部分土壤存在潛在生態(tài)風險,PAHs人體致癌風險尚在可接受范圍內(nèi)[16-19]。

    北京市東南郊灌區(qū)從20世紀50年代開始利用污水灌溉,主要污水來源為北京城區(qū)污水。該灌區(qū)從2003年逐步開始利用再生水灌溉農(nóng)作物,近年灌溉面積已發(fā)展到3.857萬hm2。目前已有學(xué)者對該灌區(qū)土壤鹽堿性、肥力水平以及重金屬遷移規(guī)律進行了研究[20-22]。也有學(xué)者對該灌區(qū)部分土壤PAHs和有機氯農(nóng)藥以及地下水壬基酚空間含量進行了研究[23-27],但這些關(guān)于土壤PAHs的研究僅集中在灌區(qū)某渠道周圍,整個灌區(qū)土壤-植物系統(tǒng)PAHs含量水平還未見報道,同時灌區(qū)PAHs對人體健康風險水平也不明確。北京市東南郊灌區(qū)是北京地區(qū)主要農(nóng)業(yè)生產(chǎn)基地之一,該灌區(qū)土壤和農(nóng)產(chǎn)品PAHs含量水平關(guān)系著人體健康和環(huán)境安全。本文對該灌區(qū)PAHs污染水平和人體健康風險進行調(diào)查研究,以期為該灌區(qū)以及其他相似區(qū)域土壤、糧食安全問題分析提供理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況及樣品采集

    研究區(qū)位于北京東南郊,其范圍為北緯39°26′~40°02′,東經(jīng)116°32′~116°43′(圖1)。該研究區(qū)屬暖溫帶半濕潤大陸性季風氣候,多年平均氣溫11~12℃,多年平均日照時數(shù)2 459 h,多年平均降水量565 mm,主要集中在6~8月。研究區(qū)位于永定河和潮白河沖積扇下游部分,從西北向東南地勢緩慢遞減。研究區(qū)表層土壤粘粒、粉粒及砂礫質(zhì)量百分數(shù)為10.5%~27.5%、46.0%~78.5%和1.5%~43.5%,土壤質(zhì)地為粘壤土和壤土。該研究區(qū)部分區(qū)域從20世紀50年代開始引用城市污水進行農(nóng)業(yè)灌溉,從2003年開始,該區(qū)域陸續(xù)引用黃村高碑店、小紅門和污水處理廠的出水(再生水)進行農(nóng)業(yè)灌溉,污灌歷史為20~40年(圖1中區(qū)域1、區(qū)域2和區(qū)域3污灌歷史分別為40 a、30 a和20 a)。

    圖1 研究區(qū)示意圖Fig.1 Schematic diagram of study area

    為探討研究區(qū)不同污灌歷史年限區(qū)域土壤和作物PAHs含量及污染水平,在各區(qū)域內(nèi)沿主要灌溉河道和渠道均勻布置采樣點,選擇種植面積相對較大的田塊為具體采樣監(jiān)測點位,研究區(qū)內(nèi)共設(shè)置31個土壤監(jiān)測點,分屬20年、30年和40年3種不同污灌歷史年限,如圖1所示。2015年6月和9月進行土壤和植物采樣,在各監(jiān)測點位內(nèi)采集土樣時,從10 m×10 m正方形的4個頂點處取表層0~20 cm土壤樣品組成混合樣品,土樣取回后于-4℃冷藏,在-20℃進行冷凍干燥。干燥后的樣品用研缽研磨,過50目篩。過篩后的樣品置于玻璃瓶中,保存在-20℃的環(huán)境下待測。

    作物收獲時(6月中旬和9月下旬)在各土壤監(jiān)測點位采集植物樣品,包括冬小麥籽粒、夏玉米籽粒和夏季蔬菜可食用部位,由于采樣時間問題導(dǎo)致某些監(jiān)測點位僅6月或9月采集到了作物樣品,不同區(qū)域作物采樣情況見表1。采集的作物樣品放入保溫箱中運輸至實驗室,依次用自來水、蒸餾水沖洗表面土壤和雜質(zhì),用濾紙沾干表面水分在4℃下保存。用于檢測有機污染物的土壤樣品和作物采用鋁箔紙包裹,以避免二次污染。

    表1 農(nóng)作物樣品數(shù)Tab.1 Number of crop samples

    1.2 土壤和植物中PAHs的索氏提取與測定方法

    取10 g干燥過的土樣,加入0.1 mg/kg的PAHs替代物(屈-d12和二萘嵌苯-d12,屈-d12作為2~4環(huán)PAHs的替代物,二萘嵌苯-d12作為5~6環(huán)PAHs的替代物),攪拌后密閉過夜,用濾紙包好后放入索氏提取器進行處理,以1∶1(體積比例)丙酮和甲醇混合液220 mL作為提取液。索氏提取12 h后,把提取液用50 g無水硫酸鈉過濾脫水,用約15 mL丙酮和甲醇混合液進行潤洗。脫水后的提取液用旋轉(zhuǎn)蒸發(fā)儀(50℃)和氮吹儀(50℃)濃縮至0.8~1.5 mL后經(jīng)過0.22 μm濾膜過濾,然后轉(zhuǎn)移至1.5 mL的樣品瓶(最后液體在0.5~1.0 mL之間),放入冰箱保存至待測。

    植物在65℃的真空干燥箱中干燥12 h,取一定量(小麥和玉米籽粒取5 g,蔬菜取2 g)干燥過的樣品,加入一定量的PAHs替代物(0.1 mg/kg),攪拌后密閉過夜,用濾紙包好后放入索氏提取器進行處理,220 mL正己烷作為提取液。索氏提取12 h后,將提取液用50 g無水硫酸鈉過濾脫水,用約15 mL正己烷溶液進行潤洗。脫水后的提取液用旋轉(zhuǎn)蒸發(fā)儀(50℃)和氮吹儀(50℃)濃縮至0.8~1.5 mL后經(jīng)0.22 μm濾膜過濾(采用柱凈化會增加目標物的損失率,采用膜過濾可以提高目標物的回收率),然后轉(zhuǎn)移至1.5 mL的樣品瓶(最后在0.5~1.0 mL之間),放入冰箱保存至待測。

    使用氣相色譜(7890A)/質(zhì)譜聯(lián)用儀(5975C)進行樣品PAHs分析,包括萘(Nap)、苊稀(Acy)、苊(Ace)、芴(Flu)、菲(Phe)、蒽(Ant)、熒蒽(Fth)、芘(Pye)、苯并(a)蒽(BaA)、屈(Chr)、苯并(b)熒蒽(BbF)、苯并(k)熒蒽(BkF)、苯并(a)芘(BaP)、茚苯(1,2,3-cd)芘(IcdP)、二苯并(a,h)蒽(DahA)、苯并(g,h,i)苝(BghiP)共16種化合物。氣相色譜(7890A)-質(zhì)譜(5975C)聯(lián)用儀采用DB-5MS型毛細管柱:30 m×0.25 mm×0.25 μm(美國安捷倫公司生產(chǎn))。進樣口溫度280℃,無分流進樣,進樣1 μL,離子源溫度230℃,四極桿溫度150℃,氣質(zhì)接口溫度280℃。GC爐溫采用程序升溫,40℃保持2 min,5℃/min升溫至290℃,保持4 min。樣品分析時采用SIM掃描模式,根據(jù)特征峰和保留時間進行定性分析,根據(jù)基峰面積進行定量分析,各物質(zhì)基峰和特征峰見表2。

    質(zhì)量保證和質(zhì)量控制:每批次試劑均分析試劑空白,每批次樣品至少做3個空白,空白測試結(jié)果均低于檢出限??瞻准訕嘶厥章试?0%~120%之間,樣品替代物屈-d12和二萘嵌苯-d12的絕對回收率為50%~70%(濃縮液中屈-d12和二萘嵌苯-d12的量與加入樣品中的替代物量的比值);儀器每12 h做一次溶劑空白,檢查儀器污染狀況;本研究方法中萘、苊稀、苊、芴、菲、蒽、熒蒽、芘、苯并(a)蒽、屈、苯并(b)熒蒽、苯并(k)熒蒽、苯并(a)芘、茚苯(1,2,3-cd)芘、二苯并(a,h)蒽、苯并(g,h,i)苝檢出限分別為0.238、0.036、0.101、0.097、0.034、0.065、0.064、0.110、0.044、0.078、0.048、0.090、0.025、0.067、0.047、0.049 μg/kg。

    表2 多環(huán)芳烴及替代物的基峰和特征峰Tab.2 Base peak and characteristic peak of PAHs and substitutes

    1.3 健康風險評價方法

    采用美國能源部風險評估信息系統(tǒng)中化學(xué)物質(zhì)風險模型用戶指南中的暴露模型計算人群經(jīng)口攝入、呼吸吸入和皮膚接觸等途徑暴露污染物的日平均暴露劑量[28],計算公式為

    ADD=

    式中ADD——化學(xué)物質(zhì)日均暴露劑量,mg/(kg·d)

    其余參數(shù)定義和參考值見表3,參考值取值參考美國能源部風險評估信息系統(tǒng)中化學(xué)物質(zhì)風險模型用戶指南[28]和美國環(huán)境保護署超級基金人體健康風險評估指導(dǎo)和土壤篩選水平補充指導(dǎo)[29-31]。

    表3 健康風險評價暴露參數(shù)Tab.3 Exposure parameters of health risk assessment

    注:VF計算方法參考美國環(huán)境保護署土壤篩選導(dǎo)則用戶指南[31]。

    根據(jù)PAHs各組分致癌性,進行非致癌風險評價和致癌風險評價,計算公式和相應(yīng)參數(shù)取值參考美國能源部風險評估信息系統(tǒng)中化學(xué)物質(zhì)風險模型用戶指南、美國環(huán)境保護署超級基金人體健康風險評估指導(dǎo)和土壤篩選水平補充指導(dǎo)和我國污染場地風險評估技術(shù)導(dǎo)則[28-30,32-33]。

    HI=ADD/RfD

    (5)

    (6)

    RfDs=RfDmABSgi

    (7)

    SFs=SFm/ABSgi

    (8)

    RfDb=RfCbIR3a/BWa

    (9)

    SFb=IURbBWa/IR3a

    (10)

    式中HI——非致癌PAHs風險指數(shù),比值超過1時,認為會對人體健康產(chǎn)生危害

    RfD——非致癌污染物參考劑量,mg/(kg·d)

    RfDs——經(jīng)皮膚參考劑量

    RfDm——經(jīng)口參考劑量

    RfDb——經(jīng)呼吸參考劑量

    RfCb——經(jīng)呼吸參考濃度

    Risk——致癌PAHs暴露風險,若低劑量暴露計算致癌風險值大于0.01時則用高劑量暴露計算

    SF——污染物的致癌斜率因子,(kg·d)/mg

    SFs——經(jīng)皮膚的致癌斜率因子

    SFm——經(jīng)口的致癌斜率因子

    SFb——經(jīng)呼吸的致癌斜率因子

    ABSgi——經(jīng)腸胃吸收的污染物分數(shù),取0.5

    IR3a——成人日均空氣攝入量

    BWa——成人體重

    IURb——呼吸吸入單位致癌因子

    PAHs的揮發(fā)因子、非致癌參考劑量和致癌斜率因子參考值見表4。

    用Excel 2010和Origin軟件處理數(shù)據(jù)和作圖,利用SPSS 20.0對數(shù)據(jù)進行統(tǒng)計分析,包括利用方差分析中LSD法分析3個區(qū)域玉米PAHs含量差異顯著性,利用獨立樣本t檢驗法分析2個域小麥PAHs含量差異顯著性,顯著性水平選取0.05。

    表4 PAHs揮發(fā)因子(VF)、非致癌參考劑量(RfD)和致癌斜率因子(SF)Tab.4 Volatilization factor (VF), reference dose(RfD) and slope factor (SF) of PAHs

    注:D表示沒有分類。

    2 結(jié)果與分析

    2.1 灌區(qū)表層土壤PAHs空間分布

    圖2顯示灌區(qū)表層土壤PAHs含量空間分布,其含量為113.5~449.8 μg/kg,均值為258.5 μg/kg??傮w上灌區(qū)中偏北部地區(qū)土壤PAHs含量最高,其次為灌區(qū)西部、西南及南部地區(qū),東南區(qū)域表層土壤PAHs含量最低。表層土壤PAHs含量空間分布特征與區(qū)域污灌歷史空間分布(圖1)相似:灌區(qū)中偏北以及東北部分區(qū)域污灌歷史達40年,灌區(qū)西部、西南及南部區(qū)域污灌歷史達30年,灌區(qū)東南部分區(qū)域污灌歷史達20年。這說明灌區(qū)表層土壤PAHs含量受污灌年限影響,其含量隨污灌歷史年限增加而增加。

    依據(jù)MALISZEWSKA-KORDYBACH[34](1996年)土壤污染程度劃分方法(PAHs小于200 μg/kg,表明土壤未受到污染;200~600 μg/kg,受輕度污染;600~1 000 μg/kg,受到中度污染;大于1 000 μg/kg,受到嚴重污染),本研究灌區(qū)土壤采樣點中80%受到了PAHs輕度污染。

    圖2 灌區(qū)表層土壤PAHs空間分布Fig.2 Spatial distribution of PAHs in topsoil in irrigation district

    表5顯示了國內(nèi)外部分區(qū)域土壤PAHs含量。與國內(nèi)研究成果相比,本研究灌區(qū)表層土壤PAHs

    含量高于上海市崇明島農(nóng)田土壤PAHs含量[35];與大興和通州區(qū)農(nóng)田果園、安徽省合肥和蕪湖周邊蔬菜基地土壤PAHs含量接近[7-8];低于山東壽光市周邊農(nóng)田、上海市郊區(qū)蔬菜地、遼寧污灌區(qū)土壤、河北邯鄲工業(yè)區(qū)麥田土壤PAHs含量[13-14,36-37];遠低于唐山煤礦區(qū)農(nóng)田和山西小店污灌區(qū)土壤PAHs含量[12,18]。

    與國外研究成果相比,研究區(qū)表層土壤PAHs含量與韓國農(nóng)業(yè)土壤PAHs含量相近[38],低于波蘭農(nóng)田和巴基斯坦污灌區(qū)土壤PAHs含量[6,19]。與國內(nèi)外農(nóng)田土壤相比,本研究表層土壤PAHs含量處于相對較低的水平。

    2.2 灌區(qū)表層土壤PAHs各組分含量及占總量比例

    表6顯示了灌區(qū)表層土壤PAHs各檢測組分含量以及占總量的比例。PAHs中16種檢測組分在所有采集的土壤樣品中均能檢測出。土樣PAHs各檢測組分中,Phe含量最高,為26.8~103.6 μg/kg,均值占PAHs總量的21.2%;其次為BbF、BghiP和Fth,含量分別為15.5~100.0、8.2~67.9、16.7~37.5 μg/kg,均值分別占PAHs總量的15.2%、10.9%、10.5%;剩余檢測組分占PAHs比例總體上低于5%。郜紅建等[8](2012年)和吳迪等[14](2016年)研究也得出土壤PAHs各檢測組分中以Phe含量最高,分別占總量的24.8%和14%;而周海軍等[9](2013年)研究得出重點污染源企業(yè)周邊農(nóng)田土壤PAHs以BkF和BbF為主。主要是不同區(qū)域土壤PAHs各組分比例與土壤自身功能區(qū)域或所處環(huán)境有關(guān),如朱媛媛等[10](2014年)調(diào)查天津市不同功能區(qū)表層土壤PAHs各組分含量,得出遠郊區(qū)土壤PAHs各組分以Phe、Ant和Fth含量為主,南開城區(qū)土壤以IcdP和BghiP為主,這主要是受土壤周圍工礦企業(yè)和降塵等影響。

    表5 不同地區(qū)土壤PAHs含量的比較Tab.5 Comparative analysis of PAHs in soil from different regions

    表6 灌區(qū)表層土壤PAHs各檢測組分含量、占總量比例以及BaP毒性當量濃度Tab.6 Concentrations, percentage and TEQBaP of PAHs in topsoil in irrigation district

    注:a: BaP毒性當量因子TEF選值參考NISBET和LAGOY[40](1992年)研究結(jié)果。

    我國目前還沒發(fā)布土壤PAHs含量質(zhì)量標準,國際上也沒有統(tǒng)一評價標準,通常以荷蘭土壤質(zhì)量標準進行評價(《荷蘭土壤和地下水環(huán)境質(zhì)量(DIV,2009)》)[39],該標準中10種PAHs(Nap、Ant、Fth、BaA、Chr、Phe、BkF、BaP、IcdP、BghiP)的干預(yù)值為40 mg/kg。本研究表層土壤中相應(yīng)10種PAHs的總量為0.083~0.307 mg/kg,均遠低于此標準。

    PAHs在土壤中長期累積可能會給生態(tài)系統(tǒng)帶來較大風險,因此需要對土壤PAHs生態(tài)風險進行評估。研究者常用BaP毒性當量濃度(Toxic Equivalence Quantity,TEQBap,為PAHs實測值與其BaP毒性當量因子TEF乘積)進行土壤PAHs潛在生態(tài)風險評估[10,40-41]。本文計算了PAHs各組分的毒性當量濃度,見表5。結(jié)果顯示灌區(qū)表層土壤16種PAHs總的毒性當量濃度均值為30.15 μg/kg。表層土壤PAHs各組分中毒性濃度較大的依次為BaP、DahA、BbF、IcdP、BkF,其毒性當量濃度均大于1 μg/kg,這5種組分對總毒性當量濃度貢獻率為96%,可見,主要是這5種組分其BaP毒性當量因子相對較高所致。周海軍等[9](2013年)在內(nèi)蒙古某化工企業(yè)周邊研究得出農(nóng)田土壤PAHs總量平均毒性當量濃度為26.8 μg/kg,且對毒性當量濃度貢獻最大為BaP、BkF、DahA和BbF,本研究結(jié)果與此相似。荷蘭土壤標準目標值規(guī)定的10種PAHs毒性當量濃度限值為33 μg/kg(由10種PAHs目標值與相應(yīng)毒性當量因子計算而得),本研究中相應(yīng)的10種PAHs總的毒性當量濃度均值為19.0 μg/kg,說明該灌區(qū)表層土壤總體上PAHs潛在生態(tài)風險不顯著。

    2.3 灌區(qū)各植物PAHs含量

    圖3顯示了灌區(qū)不同區(qū)域作物PAHs含量。其中,圖3a顯示了灌區(qū)區(qū)域1、區(qū)域2冬小麥籽粒PAHs總量均值分別為368.7 μg/kg和389.9 μg/kg,區(qū)域1、區(qū)域2、區(qū)域3夏玉米籽粒PAHs總量均值分別為93.9、87.7、113.5 μg/kg,顯著性分析顯示了小麥籽?;蛴衩鬃蚜2煌瑓^(qū)域PAHs含量無顯著差異(P>0.05),說明不同污灌歷史并未顯著影響糧食作物籽粒PAHs含量。本研究中小麥籽粒PAHs總量高于薛海全[42](2011年)在濟南商河得的小麥麥穗PAHs含量(240 μg/kg),本研究中玉米籽粒PAHs含量低于其研究中玉米種子PAHs含量(249 μg/kg),主要可能與研究區(qū)域土壤、灌溉水和大氣降塵等有關(guān)。

    圖3 灌區(qū)各作物PAHs總量Fig.3 Concentrations of PAHs in crops in irrigation district

    本研究采集的蔬菜樣品包含葉菜、莖菜、果菜和根菜,圖3b顯示蔬菜可食用部分PAHs總量為51.8~291.8 μg/kg(均值168 μg/kg),稍低于郭雪[13](2015年)研究得出的蔬菜PAHs含量(78.0~320.0 μg/kg,均值205.0 μg/kg),而高于郜紅建等[8](2012年)在安徽蔬菜基地研究得出的蔬菜PAHs含量(23.4~209.1 μg/kg,均值121.0 μg/kg)。本研究中的蔬菜PAHs含量由大到小為葉菜類、莖菜類、果菜類、根菜類,這與郭雪[13](2015年)得出蔬菜PAHs含量規(guī)律相似,主要是不同蔬菜品種對PAHs吸收除受土壤含量影響外,可能還受大氣沉降影響,而大氣沉降量與植物生長結(jié)構(gòu)密切相關(guān)。

    2.4 灌區(qū)作物PAHs各組分量及占總量比例

    表7顯示了作物PAHs各組分含量以及占總量的比例。冬小麥和夏玉米籽粒PAHs以BghiP為主,含量分別為157.7 μg/kg和35.7 μg/kg,分別占PAHs總量的41.6%和37.2%;其次為Phe,分別占總量的24.5%和20.2%。薛海全[42](2011年)在濟南商河研究得出小麥麥穗PAHs以BghiP和Phe含量最高,分別占PAHs總量的46.3%和12.9%;玉米種子PAHs含量以BghiP為主,占總量的80%;本研究結(jié)果與此類似。蔬菜PAHs各組分以Phe含量最高,均值為44.39 μg/kg,占總量比例26.4%;其次為Fth,占總量的13.5%。其他學(xué)者研究也得出蔬菜PAHs各組分中以Phe含量最高,其次為Fth、Pyr和Flu[8,13,43]。

    表7 灌區(qū)作物PAHs各檢測組分含量以及占總量比例Tab.7 Concentrations and percentage of PAHs in crops in irrigation district

    糧食作物主要吸收BghiP和Phe,蔬菜可食用部位主要吸收Phe和Fth。作物能吸收較多的Phe和Fth主要是土壤中這兩類組分含量相對較高(表6)且均屬于低環(huán)PAHs,水溶性較高,易于被作物吸收。糧食作物能吸收較多BghiP,除受土壤較高含量影響外,還受大氣降塵中PAHs各組分含量影響。有研究表明北京市大氣PM2.5和PM10中高分子PAHs含量高于低分子PAHs含量,且BghiP是PAHs含量最高的組分,PAHs含量由高到低排序為冬季、秋季或春季、夏季[44-45]。本研究中冬小麥生長期長且生長期內(nèi)降塵較多,玉米植株較高大能接觸較多的大氣沉降,所以糧食作物能吸收較多的BghiP。

    本研究中冬小麥籽粒、夏玉米籽粒BaP含量分別為0.4~4.4 μg/kg和0~1.9 μg/kg,均低于食品中糧食作物污染物限量(GB2762—2012)[46]規(guī)定限值(BaP濃度限值為5.0 μg/kg)。根據(jù)PAHs各組分BaP毒性當量因子計算得出冬小麥籽粒和夏玉米籽粒16種PAHs的BaP毒性當量濃度均值分別為6.12 μg/kg和2.45 μg/kg。冬小麥籽粒PAHs各組分中毒性濃度最大的為BaP和BghiP,其毒性當量濃度均值均大于1.5 μg/kg,主要是這2種組分的毒性當量因子或含量較高所致。夏玉米籽粒PAHs各組分中毒性濃度較大的為DahA和BkF,但均值均小于0.70 μg/kg,主要是各組分含量均較低。

    2.5 表層土壤和作物PAHs來源分析

    人為活動引起多環(huán)芳烴污染的兩類主要來源為石化燃料的不完全燃燒和未燃燒石油產(chǎn)品的泄漏,即燃燒源和石油源[47]。燃燒源產(chǎn)生的PAHs中,2、3環(huán)主要是在低溫或中溫下大量產(chǎn)生,而4環(huán)及以上PAHs主要在高溫下產(chǎn)生;石油源殘留化合物較多為2、3環(huán)PAHs;因此一般情況下4環(huán)及以上PAHs含量較低環(huán)PAHs含量呈相對優(yōu)勢分布時,表明高溫不完全燃燒為主要污染源,反之則以石油源輸入為主[48]。

    利用Origin軟件繪出土壤、作物不同環(huán)數(shù)PAHs組成三元圖(圖4)。從圖4可以看出灌區(qū)表層土壤、小麥籽粒、玉米籽粒各自PAHs組分分別相近,說明各自有相似來源;而蔬菜PAHs組分差別較大,污染來源為混合源。

    圖4 灌區(qū)表層土壤、糧食作物籽粒、蔬菜PAHs組成三元圖Fig.4 Ternary diagram of PAHs in topsoil, cereal grain and vegetables

    圖5顯示表層土壤、作物各環(huán)PAHs占總量的比例。表層土壤2、3、4、5、6環(huán)PAHs各自占總量的比例分別為1.8%、26.3%、23.9%、28.2%、19.8%,高環(huán)(4、5、6環(huán))PAHs含量占總量的71.9%,說明本研究中化石燃料的高溫燃燒是土壤PAHs主要污染源。陳素暖等[17](2010年)研究北京東南郊區(qū)再生水和清水灌溉得出土壤高環(huán)PAHs占總量比重為71.80%和58.94%,熱源為PAHs主要污染源,石鈺婷等[24](2011年)研究也得出PAHs污染主要來源為煤的燃燒和汽車尾氣排放等燃燒源。

    圖5 表層土壤、作物各環(huán)PAHs占總量比例Fig.5 Percentages of PAHs with different rings in topsoil and crops

    冬小麥和夏玉米籽粒中6環(huán)PAHs占總量比例最高,分別為43.1%和38.2%;其次為3環(huán),分別為37.9%和29.9%;4環(huán)占總量比例分別為12.1%和13.9%;2環(huán)和5環(huán)占總量比例均低于10%。糧食作物中高環(huán)(4、5、6環(huán))PAHs含量占總量58.6%~60.0%,說明本研究中高溫燃燒是糧食作物PAHs主要污染源。這與薛海全[42](2011年)在濟南研究得出結(jié)論類似,成熟期小麥籽粒和玉米籽粒均以6環(huán)PAHs含量最高,且PAHs來源均為燃燒源。

    蔬菜作物各環(huán)PAHs占總量比例規(guī)律相似,但與糧食作物不同,總體上蔬菜3環(huán)和4環(huán)PAHs占總量比例較高,分別為43.3%~67.0%和12.1%~39.6%;5、6環(huán)各自占總量比例較低,為0~16%。本研究結(jié)果與萬開等[43](2009年)在珠江三角洲和郜紅建等[8](2012年)在安徽研究結(jié)果相似,即3、4環(huán)PAHs在蔬菜體內(nèi)含量較高,而5、6環(huán)含量較低,主要是蔬菜體內(nèi)含水量較高,分子量相對較低的PAHs水溶性大,生物有效性高,容易被蔬菜吸收累積。蔬菜作物中高環(huán)PAHs含量占總量的28.0%~53.7%,說明本研究中蔬菜PAHs主要污染源為燃燒源和石油源(交通燃料泄漏)混合來源。郭雪[13](2015年)在上海市郊研究也得出蔬菜PAHs來源為燃燒源和石油類泄漏。

    2.6 人體健康風險評估

    表8顯示了灌區(qū)PAHs各組分人體健康致癌風險(土壤、作物PAHs各組分濃度含量取灌區(qū)均值)。US EPA提出一般可接受的致癌風險水平為10-6,可接受的致癌風險水平上限為10-4[49]。本研究中成人和兒童PAHs致癌風險分別為4.02×10-5和1.76×10-5,均低于致癌風險水平上限,說明癌癥風險水平在可接受范圍內(nèi);但從安全角度考慮,應(yīng)長期注意這些區(qū)域PAHs濃度水平。從暴露途徑和暴露介質(zhì)分析,PAHs致癌危害貢獻以口-作物(小麥、玉米和蔬菜)為主,其對人體健康所造成的致癌危害占總個人年致癌風險的比例分別為99.64%(成人)和98.39%(兒童),其中以口-小麥的貢獻率最高,分別為58.72%(成人)和46.22%(兒童)。本文通過土壤引起的綜合致癌風險為1.43×10-7(成人)和2.82×10-7(兒童),與馮煥銀等[50](2011年)在寧波得出土壤綜合致癌風險相近(3.17×10-7)。從PAHs各組分分析,對人體致癌風險貢獻最大的分別為BaP和DahA,各自對成人致癌風險貢獻率分別為46.72%和22.15%,各自對兒童致癌風險貢獻率分別為42.89%和26.36%,主要是這2個組分致癌斜率因子相對較高所致(表4)。

    表9顯示了灌區(qū)PAHs各組分人體健康非致癌風險(土壤、作物PAHs各組分濃度含量取灌區(qū)均值)。US EPA指出當非致癌危害指數(shù)大于1時,認為對人體健康產(chǎn)生危害。本研究中成人和兒童非致癌指數(shù)分別為2.72×10-2和4.78×10-2,均低于1,說明PAHs未對人群產(chǎn)生明顯的非致癌健康危害。從暴露途徑和暴露介質(zhì)分析,PAHs非致癌危害貢獻以口-作物(小麥、玉米和蔬菜)為主,其對人體健康所造成的非致癌危害占總個人年風險的比例分別為99.9%(成人)和99.8%(兒童),其中通過口-小麥的貢獻率分別為80.95%(成人)和71.78%(兒童)。從PAHs各組分分析,對人體非致癌風險貢獻最大的分別為BghiP和Phe,各自對成人非致癌風險貢獻率分別為46.28%和29.14%,各自對兒童非致癌風險貢獻率分別為44.28%和29.60%,主要是作物中這兩類組分含量相對較高所致(表7)。

    3 結(jié)論

    (1)整個灌區(qū)表層土壤PAHs含量為113.5~449.8 μg/kg,均值為258.5 μg/kg,依據(jù)MALISZEWSKA-KORDYBACH (1996年)土壤污染程度劃分方法,灌區(qū)土壤采樣點中80%受到了PAHs輕度污染。與荷蘭土壤環(huán)境質(zhì)量標準(10種PAHs的干預(yù)值為40 mg/kg)相比,本研究表層土壤中相應(yīng)10種PAHs的總量為0.083~0.307 mg/kg,遠低于荷蘭土壤PAHs修復(fù)干預(yù)值,說明灌區(qū)土壤PAHs不需干預(yù)修復(fù)。本研究中10種PAHs苯并(a)芘毒性當量濃度(TEQBap)均值為19.0 μg/kg,低于荷蘭土壤環(huán)境生態(tài)風險相應(yīng)標準值(33.0 μg/kg),說明灌區(qū)土壤潛在生態(tài)風險不顯著。

    表8 PAHs致癌健康風險計算結(jié)果Tab.8 Calculation results of carcinogenic health risk of PAHs

    表9 PAHs非致癌健康風險計算結(jié)果Tab. 9 Calculation results of non-carcinogenic health hazard index of PAHs

    (2)冬小麥籽粒、夏玉米籽粒和蔬菜可食用部分PAHs總量分別為368.7~389.9 μg/kg、87.7~113.5 μg/kg和51.8~291.8 μg/kg,總體上蔬菜PAHs含量由大到小為葉菜類、莖菜類、果菜類、根菜類。糧食作物籽粒PAHs主要以BghiP為主,其含量占總量的36.6%~41.6%。蔬菜PAHs各組分中以Phe含量最高,占總量比例24.0%。糧食作物籽粒BaP含量均低于食品中污染物限量(GB2762—2012)規(guī)定的限量值(5 μg/kg)。

    (3)表層土壤、糧食作物和蔬菜高環(huán)PAHs含量占總量比例分別為71.9%、58.6%~60.0%和28.0%~53.7%,說明高溫燃燒是土壤和糧食作物PAHs主要污染源,而蔬菜PAHs主要污染源為燃燒源和石油源混合來源。

    (4)成人和兒童PAHs致癌風險分別為4.02×10-5和1.76×10-5,非致癌指數(shù)分別為2.72×10-2和4.78×10-2,均低于UE EPA規(guī)定相應(yīng)標準參考值???作物是人體健康風險主要暴露途徑和介質(zhì),對致癌風險和非致癌指數(shù)貢獻率分別為98.39%~99.64%和99.8%~99.9%。對致癌風險貢獻最大的組分為BaP和DahA,貢獻率分別為42.89%~46.72%和22.15%~26.36%。對非致癌風險貢獻最大的組分是BghiP和Phe,貢獻率分別為44.28%~46.28%和29.14%~29.60%。

    1 EDWARDS N T. Polycyclic aromatic hydrocarbons (PAHs) in the terrestrial environment: A review[J]. Journal of Environmental Quality, 1983, 12(4): 427-441.

    2 JONES K C, DE VOOGT P. Persistent organic pollutants(POPs): state of the science[J]. Environmental Pollution, 1999, 100(1-3):209-221.

    3 JONES K C, STRAFORD J A, WATERHOUSE K S, et al. Increases in the polynuclear aromatic hydrocarbon content of an agricultural soil over last century[J]. Environmental Science & Technology, 1989,23(1):95-101.

    4 WILD S R, JONES K C. Polynuclear aromatic hydrocarbons in the United Kingdom environment: a preliminary source in inventory and budget[J]. Environment Pollution, 1995,101:91-108.

    5 JONES K C, GRIMMER G, JACOB J, et al. Changes in the polynuclear aromatic hydrocarbon content of wheat grain and pasture grassland and over the last century from one site in the UK[J]. Science of the Total Environment, 1989, 78:117-130.

    6 MALISZEWSKA-KORDYBACH B, SMRECZAK B, KLIMKOWICZ-PAWLAS A. Concentrations, sources, and spatial distribution of individual polycyclic aromatic hydrocarbons (PAHs) in agricultural soils in the Eastern part of the EU: Poland as a case study[J]. Science of the Total Environment, 2009, 407(12): 3746-3753.

    7 張枝煥,盧另,賀光秀,等. 北京地區(qū)表層土壤中多環(huán)芳烴的分布特征及污染源分析[J]. 生態(tài)環(huán)境學(xué)報,2011,20(4):668-675. ZHAND Z H, LU L, HE G X, et al. Distribution and sources of polycyclic aromatic hydrocarbon compounds in topsoil of Beijing, China[J]. Ecology and Environmental Sciences, 2011, 20(4): 668-675. (in Chinese)

    8 郜紅建,魏俊嶺,馬靜靜,等. 安徽省典型城市周邊土壤-蔬菜中PAHs的污染特征[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2012,31(10):1913-1919. GAO H J, WEI J L, MA J J, et al. Characteristics of PAHs contamination in soil and vegetable in typical city outskirt of Anhui Province, China[J]. Journal of Agro-Environment Science, 2012, 31(10): 1913-1919. (in Chinese)

    9 周海軍,孫文靜,團良,等.重點污染源企業(yè)周邊農(nóng)田土壤中多環(huán)芳烴污染水平與風險評價[J].環(huán)境化學(xué),2013,32(10):1976-1982. ZHOU H J, SUN W J, TUAN L, et al. Pollution level and risk assessment of polycyclic aromatic hydrocarbons in agricultural soil around pollution sources[J].Environmental Chemistry, 2013,32(10):1976-1982. (in Chinese)

    10 朱媛媛,田靖,魏恩琪,等.天津市土壤多環(huán)芳烴污染特征、源解析和生態(tài)風險評價[J].環(huán)境化學(xué),2014,33(2):248-255. ZHU Y Y, TIAN J, WEI E Q, et al. Characteristics, sources apportionment and ecological risks assessment of polycyclic aromatic hydrocarbons in soils of Tianjin, China[J]. Environmental Chemistry, 2014,33(2): 248-255. (in Chinese)

    11 KWON H O, CHOI S D. Polycyclic aromatic hydrocarbons (PAHs) in soils from a multi-industrial city, South Korea[J]. Science of the Total Environment, 2014, 470 -471: 1494-1501.

    12 趙歐亞,馮圣東,石維,等.煤礦區(qū)農(nóng)田土壤多環(huán)芳烴生態(tài)風險評估方法比較[J].安全與環(huán)境學(xué)報,2015,15(2):352-358. ZHAO O Y, FENG S D, SHI W, et al. A method for assessing the ecological risks due to the polycyclic aromatic hydrocarbons in the farming soil near coal mine area[J]. Journal of Safety and Environment, 2015, 15(2):352-358. (in Chinese)

    13 郭雪. 上海市郊區(qū)土壤-蔬菜系統(tǒng)中多環(huán)芳烴污染效應(yīng)研究[D]. 上海:華東師范大學(xué),2015. GUO X. Pollution effect of polycyclic aromatic hydrocarbons in soil-vegetable system in Suburban Shanghai[D]. Shanghai: East China Normal University, 2015. (in Chinese)

    14 吳迪,汪宜龍,劉偉健,等.河北邯鄲鋼鐵冶煉區(qū)周邊麥田土和小麥籽粒的多環(huán)芳烴含量及其組分譜特征[J].環(huán)境科學(xué),2016,37(2):740-749. WU D, WANG Y L, LIU W J, et al. Concentrations and component profiles PAHs in surface soils and wheat grains from the cornfields close to the steel smelting industry in Handan' Hebei Province[J]. Environmental Science, 2016, 37(2):740-749.(in Chinese)

    15 BI X, LUO W, GAO J J, et al. Polycyclic aromatic hydrocarbons in soils from the Central-Himalaya region: distribution, sources, and risks to humans and wildlife[J]. Science of the Total Environment,2016,556:12-22.

    16 曹云者,施烈焰,李麗和,等.渾蒲污灌區(qū)表層土壤中多環(huán)芳烴的健康風險評價[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2008,27(2):542-548. CAO Y Z, SHI L Y, LI L H, et al. Health risk assessment of polycyclic aromatic hydrocarbons in surface soils from Hunpu Wastewater Irrigated area[J].Journal of Agro-Environment Science, 2008,27(2):542-548. (in Chinese)

    17 陳素暖,何江濤,金愛芳,等.多環(huán)芳烴在不同灌區(qū)土壤剖面的分布特征研究[J].環(huán)境科學(xué)與技術(shù),2010,33(10):10-14. CHEN S N, HE J T, JIN A F, et al. Distribution Characteristics of PAHs in soil profiles of different irrigation areas[J]. Environmental Science & Technology, 2010,33(10):10-14. (in Chinese)

    18 黨晉華,劉利軍,趙穎,等.太原市小店污灌區(qū)土壤重金屬和多環(huán)芳烴復(fù)合污染現(xiàn)狀評價[J].山西農(nóng)業(yè)科學(xué),2013,41(10):1085-1088. DANG J H, LIU L J, ZHAO Y, et al. Assessment of Farmland with Heavy Metal and PAHs Pollution in Xiaodian Sewage Irrigation Area, Taiyuan[J]. Journal of Shanxi Agricultural Sciences, 2013,41(10):1085-1088. (in Chinese)

    19 MUHAMMAD W, SARDAR K, CAI C, et al. Quantification of PAHs and health risk via ingestion of vegetable in Khyber Pakhtunkhwa Province, Pakistan[J]. Science of the Total Environment, 2014, 497(1): 448-458.

    20 楊軍,陳同斌,雷梅,等. 北京市再生水灌溉對土壤、農(nóng)作物的重金屬污染風險[J]. 自然資源學(xué)報,2011,26(2):209-217. YANG J, CHEN T B, LEI M, et al. Assessing the effect of irrigation with reclaimed water: the soil and crop pollution risk of heavy metals[J]. Journal of Natural Resources, 2011,26(2):209-217. (in Chinese)

    21 高軍,王會肖,劉海軍,等. 北京市再生水灌溉對土壤質(zhì)量的影響研究[J]. 北京師范大學(xué)學(xué)報:自然科學(xué)版,2012,48(5):572-576. GAO J, WANG H X, LIU H J, et al. Effects of reclaimed water irrigation upon soil quality in Beijing[J]. Journal of Beijing Normal University:Natural Science, 2012, 48(5):572-576. (in Chinese)

    22 馬闖,楊軍,雷梅,等. 北京市再生水灌溉對地下水的重金屬污染風險[J]. 地理研究,2012,31(12):2250-2258. MA C, YANG J, LEI M, et al. Assessing the effect of reclaimed water irrigation on groundwater pollution of heavy metals in Beijing[J]. Geographical Research, 2012, 31(12): 2250-2258. (in Chinese)

    23 何江濤,金愛芳,陳素暖,等. 北京東南郊再生水灌區(qū)土壤PAHs污染特征[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2010,29(4):666-673. HE J T, JIN A F, CHEN S N, et al. Distribution characteristics of soil PAHs in reclaimed water irrigation area in the Southeastern Suburb of Beijing[J]. Journal of Agro-Environment Science, 2010, 29(4): 666-673. (in Chinese)

    24 石鈺婷,何江濤,金愛芳. 北京市東南郊不同灌區(qū)表層土壤中PAHs來源解析[J]. 現(xiàn)代地質(zhì),2011,25(2):393-400. SHI Y T, HE J T, JIN A F. Sources apportionment of PAHs in the surface soil of different irrigation areas in Southeast Suburb of Beijing[J]. Geoscience, 2011, 25(2): 393-400. (in Chinese)

    25 金愛芳,何江濤,陳素暖,等. 北京東南郊污灌區(qū)土壤有機氯農(nóng)藥的垂向分布特征[J]. 環(huán)境科學(xué)學(xué)報,2009,29(9):1970-1977. JIN A F, HE J T, CHEN S N, et al. Distribution of organochlorine pesticides in soil profiles from wastewater irrigated areas in the southeast suburbs of Beijing[J]. Acta Scientiae Circumstantiae, 2009, 29 (9): 1970-1977. (in Chinese)

    26 馬文潔,何江濤,金愛芳,等. 北京市郊再生水灌區(qū)土壤有機氯農(nóng)藥垂向分布特征[J]. 生態(tài)環(huán)境學(xué)報,2010,19(7):1675-1681. MA W J, HE J T, JIN A F, et al. Verticle distribution of organochlorine pesticides in reclaimed water irrigation area in Beijing suburbs[J]. Ecology and Environmental Sciences, 2010, 19(7): 1675-1681. (in Chinese)

    27 WANG S Y, WU W Y, LIU F, et al. Spatial distribution and migration of nonylphenol in groundwater following long-term wastewater irrigation[J]. Journal of Contaminant Hydrology, 2015,177-178: 85-92.

    28 RAIS(Risk Assessment Information System). Risk exposure models for chemicals user’s guide [EB/OL]. [2016-08-20]. http:∥rais.ornl.gov/tools/rais-chemical-risk-guide.html.

    29 U S. EPA. Risk assessment guidance for superfund Volume Ⅰ:human health evaluation manual (Part E, Supplemental guidance for dermal risk assessment)[EB/OL]. [2016-08-20]. http:∥www.epa.gov/oswer/riskassessment/ragse/index.htm.

    30 U S. EPA. Supplement guidance for developing soil screening levels for superfound sites[EB/OL].(2010-10-08)[2016-08-20]. http:∥www.epa.gov/superfund/health/conmedia/soil/pdfs/ssg-main.pdf

    31 U S. EPA. Soil screening guidance:user's guide[EB/OL]. [2016-08-20]. http:∥www.epa.gov/superfund/resources/soil/ssg496.pdf.

    32 U S. EPA. Risk assessment guidance for superfund Volume Ⅰ:human health evaluation manual(Part A) interim Final[EB/OL]. (2010-10-08)[2016-08-20].http:∥www.epa.gov/oswer/riskassessment/ragsa/index.htm.

    33 HJ25.3—2014 污染場地風險評估技術(shù)導(dǎo)則[S]. 2014.

    34 MALISZEWSKA-KORDYBACH B. Polycyclic aromatic hydrocarbons in agricultural soils in Poland: preliminary proposals for criteria to evaluate the level of soil contamination[J]. Applied Geochemistry, 1996, 11(12): 121-127.

    35 呂金剛,畢春娟,陳振樓,等. 上海市崇明島農(nóng)田土壤中多環(huán)芳烴分布和生態(tài)風險評價[J]. 環(huán)境科學(xué),2012,33(12):4270-4275. LV J G, BI C J, CHEN Z L, et al. Distribution and ecological risk assessment of polycyclic aromatic hydrocarbons in agricultural soil of the Chongming Island in Shanghai[J]. Environmental Science, 2012,33(12):4270-4275. (in Chinese)

    36 馮岸紅,朱智成,陳社軍,等. 壽光土壤中多環(huán)芳烴的污染特征及風險評估[J]. 中國環(huán)境科學(xué),2013,33(9):1607-1614. FENG A H, ZHU Z C, CHEN S J, et al. Spatial distribution and risk assessment of polycyclic aromatic hydrocarbons in soils of Shouguang City, North China[J]. China Environmental Science, 2013,33(9):1607-1614. (in Chinese)

    37 王洪濤,張麗華.遼寧省污灌區(qū)土壤中多環(huán)芳烴的測定[J].環(huán)境保護與循環(huán)經(jīng)濟,2009,29(5):30-31.

    38 NAM J J, SONG B H, EOM K C, et al. Distribution of polycyclic aromatic hydrocarbons in agricultural soil in South Korea[J]. Chemosphere, 2003, 50 (10): 1281-1289.

    39 MINISTERIE VAN VOLKSHURSVESTING. Rumtelyke orderning en Milieubeheer (VROM)[EB/OL].[2011-04-03]. http:∥www.rijksoverheid.nl/ministeries/vrom.

    40 NISBET J C T, LAGOY P K. Toxic Equivalence Factors(TEFs) for polycyclic aromatic hydrocarbons (PAHs)[J]. Regulatory Toxicology and Pharmacology,1992,16(3): 290-300.

    41 TSAI P J, SHIH T S, CHEN H L, et al. Assessing and predicting the exposures of polycyclic aromatic hydrocarbons (PAHs) and their carcinogenic potencies from vehicle engine exhausts to highway toll station workers[J]. Atmospheric Environment, 2004, 38(2): 333-343.

    42 薛海全.典型農(nóng)作物中多環(huán)芳烴和多氯聯(lián)苯的分布、累積規(guī)律[D]. 濟南:山東大學(xué),2011. XUE H Q. Accumulations and distributions of PAHs, PCBs in typical crops[D]. Ji’nan: Shandong University, 2011.(in Chinese)

    43 萬開,江明,楊國義,等.珠江三角洲典型城市蔬菜中多環(huán)芳烴分布特征[J].土壤,2009,41(4):583-587. WAN K, JIANG M, YANG G Y, et al. Distribution characteristics of PAHs in vegetables of typical city in Pearl River Delta[J].Soils, 2009,41(4):583-587.(in Chinese)

    44 于國光,王鐵冠,吳大鵬,等.北京市大氣氣溶膠中多環(huán)芳烴的研究[J].中國礦業(yè)大學(xué)學(xué)報,2008,37(1):72-78. YU G G, WANG T G, WU D P, et al. A study of polycyclic aromatic hydrocarbons present in Beijing aerosols[J].Journal of China University of Mining & Technology, 2008,37(1):72-78.(in Chinese)

    45 王超,張霖琳,刀谞,等.京津冀地區(qū)城市空氣顆粒物中多環(huán)芳烴的污染特征及來源[J].中國環(huán)境科學(xué),2015,35(1):1-6. WANG C, ZHANG L L, DAO X, et al. Pollution characteristics and source identification of polycyclic aromatic hydrocarbons in airborne particulates of Beijing-Tianjin-Hebei Region, China[J].China Environmental Science, 2015,35(1):1-6. (in Chinese)

    46 GB2672—2012 食品中污染物限量[S].2012.

    47 FABBRI D, VASSURA I, SUN C G, et al. Source apportionment of polycyclic aromatic hydrocarbons in a coastal lagoon by molecular and isotopic characterization[J]. Marine Chemistry, 2003, 84: 123-135.

    48 MCCREADY S, SLEE D J, BIRCH G F, et al. The distribution of polycyclic aromatic hydrocarbons in surficial sediments of Sydney Harbor, Australia[J]. Marine Pollution Bulletin, 2000, 40(11):999-1006.

    49 U S. EPA. Risk assessment guidance for superfund (RAGS), Volume Ⅰ:human health evaluation manual(HHEM), supplemental guidance, standard default exposure factors, interim guidance, Directive. Washington, DC, OSWER: Office of Emergency and Remedial Response, 9285.6-03.1,1991.

    50 馮煥銀,傅曉欽,趙倩,等.寧波土壤中多環(huán)芳烴的健康風險評價[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2011,30(10):1998-2004. FENG H Y, FU X Q, ZHAO Q, et al. Health risk assessment of polycyclic aromatic hydrocarbons in soils of Ningbo area, China[J]. Journal of Agro-Environment Science, 2011, 30(10): 1998-2004.(in Chinese)

    ContaminationandHealthRiskAssessmentofPAHsinIrrigationDistrictinSoutheasternSuburbofBeijing

    LI Yan1,2GU Hua1,2HUANG Guanhua3HUANG Quanzhong3LIU Honglu1,2

    (1.BeijingWaterScienceandTechnologyInstitute,Beijing100048,China2.BeijingEngineeringTechniqueResearchCenterforExplorationandUtilizationofNon-conventionalWaterResourcesandWaterUseEfficiency,Beijing100048,China3.CollegeofWaterResourcesandCivilEngineering,ChinaAgriculturalUniversity,Beijing100083,China)

    Wastewater irrigation was adopted in irrigation district in southeastern suburb of Beijing from 1950s to 2003, the main source of wastewater was industrial and domestic wastewater from Beijing. Reclaimed water was used for irrigation since 2003 in the irrigation district. To investigate the residual levels of polycyclic aromatic hydrocarbons (PAHs) in topsoil and crops in the irrigation district, totally 31 samples of soil and 38 samples of crops were collected, and the concentrations of 16 US EPA (United States Environmental Protection Agency) priority PAHs were determined by gas chromatography equipped with a mass spectrometry detector (GC-MS) in 2015. Results showed that the PAHs in topsoil in the district were ranged from 113.5 μg/kg to 449.8 μg/kg, with a mean value of 258.5 μg/kg. And 80% of the collected soil samples were slightly polluted by PAHs. The sums of 10 PAHs were 83~307 μg/kg, which were lower than the standard of Netherlands (the intervention value was 40mg/kg). The average value of bap toxic equivalence quantity (TEQBap) for 10 PAHs in topsoil was 19.0 μg/kg, which was lower than the standard of Netherland (33.0 μg/kg), indicating that there was no remarkable potential ecological risk. The concentrations of winter wheat grain, summer maize grain and vegetables were 368.7~389.9 μg/kg, 87.7~113.5 μg/kg and 51.8~291.8 μg/kg, respectively. Benzo(g,h,i)perylene was the dominant component in winter wheat grain and summer maize grain, which accounted for 36.6%~41.6% of the 16 PAHs concentrations. For vegetables, the highest individual PAH was phenanthrene, which represented 24.0% of the 16 PAHs. The concentrations of benzo(a)pyrene in cereal grains were lower than the Chinese standard limits of 5 μg/kg in food. The high temperature combustion was the major pollution source for PAHs in topsoil and cereal grains, while the mixed sources of fossil fuel combustion and oil source were the pollution source for PAHs in vegetables. The carcinogenic risk of adult and child caused by PAHs were 4.02×10-5and 1.76×10-5, respectively; the non-carcinogenic hazard index of PAHs for adult and child were 2.72×10-2and 4.78×10-2, respectively; all of them were lower than the threshold values. Dietary intake was the major route of human exposure, which accounted for 98.39%~99.64% of carcinogenic risk and 99.8%~99.9% of non-carcinogenic hazard index, respectively. The carcinogenic risk caused by benzo(a)pyrene and dibenzo(a,h)anthracene were relatively higher, which accounted for 42.89%~46.72% and 22.15%~26.36% of the total carcinogenic risk, respectively. The non-carcinogenic hazard index of benzo(g,h,i)perylene and phenanthrene were relatively higher, which accounted for 44.28%~46.28% and 29.14%~29.60% of the total non-carcinogenic hazard index, respectively.

    irrigation district in southeastern suburb of Beijing; polycyclic aromatic hydrocarbons; soil; winter wheat; summer maize; vegetables; risk assessment

    X171.5

    A

    1000-1298(2017)09-0237-13

    10.6041/j.issn.1000-1298.2017.09.030

    2017-06-15

    2017-07-27

    國家自然科學(xué)基金項目(51339007)

    李艷(1986—),女,博士生,主要從事農(nóng)業(yè)水肥高效利用機理和技術(shù)研究,E-mail: liyan7986@126.com

    劉洪祿(1963—),男,教授級高級工程師,博士,主要從事再生水灌溉和農(nóng)業(yè)節(jié)水研究,E-mail: liuhonglu@yeah.net

    猜你喜歡
    表層總量籽粒
    籽粒莧的飼用價值和高產(chǎn)栽培技術(shù)
    “十三五”期間山西省與10省簽約糧食總量2230萬噸
    2020年全國農(nóng)民工總量比上年減少517萬人
    籽粒莧的特性和種植技術(shù)
    半潛式平臺表層卡套管處理與認識
    海洋石油(2021年3期)2021-11-05 07:43:10
    水體表層沉積物對磷的吸收及釋放研究進展
    為何化肥淡儲總量再度增加
    玉米機械脫粒籽粒含水量與破碎率的相關(guān)研究
    商麥1619 籽粒灌漿的特性
    氬弧熔覆原位合成Ti(C,N)-WC增強鎳基表層復(fù)合材料的研究
    焊接(2015年6期)2015-07-18 11:02:25
    在线观看免费视频日本深夜| 亚洲av成人不卡在线观看播放网| 亚洲国产精品999在线| АⅤ资源中文在线天堂| www.www免费av| www国产在线视频色| 国内精品久久久久精免费| 国产视频一区二区在线看| 亚洲国产看品久久| 日本免费a在线| 午夜精品久久久久久毛片777| 最近最新中文字幕大全电影3 | 亚洲电影在线观看av| 午夜a级毛片| 免费女性裸体啪啪无遮挡网站| x7x7x7水蜜桃| 免费一级毛片在线播放高清视频| 在线观看免费视频日本深夜| 久久精品国产亚洲av香蕉五月| 一夜夜www| 久99久视频精品免费| 国产又黄又爽又无遮挡在线| av免费在线观看网站| 免费在线观看影片大全网站| 色综合亚洲欧美另类图片| 亚洲国产精品成人综合色| 成年免费大片在线观看| 久久 成人 亚洲| 国产99白浆流出| 亚洲成人国产一区在线观看| 18美女黄网站色大片免费观看| 精品久久久久久成人av| 久久久久免费精品人妻一区二区 | 午夜免费激情av| 法律面前人人平等表现在哪些方面| 欧美精品啪啪一区二区三区| 亚洲国产精品sss在线观看| 极品教师在线免费播放| 精品无人区乱码1区二区| 日日爽夜夜爽网站| 老司机靠b影院| 亚洲男人天堂网一区| 丝袜人妻中文字幕| 欧美大码av| 变态另类成人亚洲欧美熟女| 一区二区日韩欧美中文字幕| 亚洲国产精品成人综合色| 国产亚洲欧美精品永久| 女性被躁到高潮视频| 天天躁狠狠躁夜夜躁狠狠躁| 免费看美女性在线毛片视频| 欧美色欧美亚洲另类二区| 看黄色毛片网站| 伊人久久大香线蕉亚洲五| 美女高潮到喷水免费观看| 在线观看免费午夜福利视频| 一级毛片精品| 午夜视频精品福利| 亚洲精品美女久久av网站| 好男人在线观看高清免费视频 | 久久人妻福利社区极品人妻图片| 97超级碰碰碰精品色视频在线观看| 一级毛片精品| 日韩精品青青久久久久久| 成人免费观看视频高清| 午夜老司机福利片| 久久 成人 亚洲| 久久精品人妻少妇| 精品久久久久久久毛片微露脸| 国产麻豆成人av免费视频| 免费观看精品视频网站| 一区二区三区高清视频在线| 亚洲熟妇熟女久久| 免费在线观看黄色视频的| 国产1区2区3区精品| 国内少妇人妻偷人精品xxx网站 | 亚洲久久久国产精品| 国产av又大| 免费观看精品视频网站| 久久精品人妻少妇| 国产精品精品国产色婷婷| 老司机午夜十八禁免费视频| 国产精品亚洲美女久久久| 精华霜和精华液先用哪个| 亚洲av熟女| 夜夜夜夜夜久久久久| 无限看片的www在线观看| 午夜福利在线在线| 久久中文字幕一级| 中出人妻视频一区二区| 国产精品98久久久久久宅男小说| 国产精品久久久久久精品电影 | 国产精品久久久久久亚洲av鲁大| 欧美日韩亚洲综合一区二区三区_| 国产精品亚洲av一区麻豆| 在线十欧美十亚洲十日本专区| 美女午夜性视频免费| 夜夜爽天天搞| 久久久久精品国产欧美久久久| 国产精品亚洲一级av第二区| 久久久久久久久久黄片| 欧美色视频一区免费| 亚洲 欧美一区二区三区| ponron亚洲| 91成年电影在线观看| 免费在线观看成人毛片| 熟妇人妻久久中文字幕3abv| 久久狼人影院| 老鸭窝网址在线观看| 欧美大码av| 悠悠久久av| 色综合站精品国产| 首页视频小说图片口味搜索| 国产精品国产高清国产av| 婷婷六月久久综合丁香| 久久久久九九精品影院| 最近最新中文字幕大全电影3 | 国产成人影院久久av| 国产在线精品亚洲第一网站| 18禁观看日本| 久久人妻福利社区极品人妻图片| 国产激情久久老熟女| 日本成人三级电影网站| 国产精品一区二区三区四区久久 | 亚洲成人免费电影在线观看| 久久精品影院6| 无人区码免费观看不卡| 一本大道久久a久久精品| 国产一区在线观看成人免费| 这个男人来自地球电影免费观看| 久久青草综合色| 久久欧美精品欧美久久欧美| 亚洲七黄色美女视频| 日本在线视频免费播放| 亚洲精品国产区一区二| 天天添夜夜摸| 老司机福利观看| 国产精品1区2区在线观看.| 国产私拍福利视频在线观看| 国产精品美女特级片免费视频播放器 | 亚洲七黄色美女视频| 高清在线国产一区| 国产高清视频在线播放一区| 女性生殖器流出的白浆| 1024视频免费在线观看| 九色国产91popny在线| 在线十欧美十亚洲十日本专区| 日本五十路高清| 最近最新中文字幕大全免费视频| 中文在线观看免费www的网站 | 国产精品日韩av在线免费观看| 国产97色在线日韩免费| 国产精品99久久99久久久不卡| 免费在线观看完整版高清| 麻豆成人午夜福利视频| 亚洲一区二区三区不卡视频| 亚洲片人在线观看| 给我免费播放毛片高清在线观看| 国产又黄又爽又无遮挡在线| www.自偷自拍.com| 69av精品久久久久久| 久久久久精品国产欧美久久久| 少妇粗大呻吟视频| 中文字幕另类日韩欧美亚洲嫩草| svipshipincom国产片| 色哟哟哟哟哟哟| 国产一区二区三区在线臀色熟女| 久久久国产成人精品二区| 国产精品免费视频内射| 久久精品影院6| 成人国语在线视频| 女同久久另类99精品国产91| 91成年电影在线观看| 美女扒开内裤让男人捅视频| 在线观看www视频免费| 久久久久久大精品| 熟女少妇亚洲综合色aaa.| 精品国产亚洲在线| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久| 免费在线观看影片大全网站| 精品免费久久久久久久清纯| 久久热在线av| 99国产极品粉嫩在线观看| 男女视频在线观看网站免费 | 12—13女人毛片做爰片一| 欧美性猛交╳xxx乱大交人| 美女扒开内裤让男人捅视频| 老熟妇仑乱视频hdxx| 少妇裸体淫交视频免费看高清 | 色播在线永久视频| 国内揄拍国产精品人妻在线 | 亚洲专区国产一区二区| 91成年电影在线观看| 日本三级黄在线观看| 日韩欧美 国产精品| а√天堂www在线а√下载| 一本精品99久久精品77| 午夜福利在线在线| 99riav亚洲国产免费| 人人妻人人看人人澡| 久久婷婷人人爽人人干人人爱| 可以在线观看的亚洲视频| 日本熟妇午夜| 美女高潮喷水抽搐中文字幕| 欧美丝袜亚洲另类 | www日本在线高清视频| 99热只有精品国产| 亚洲成国产人片在线观看| 成人特级黄色片久久久久久久| 日韩视频一区二区在线观看| 色尼玛亚洲综合影院| 久久精品人妻少妇| 一边摸一边抽搐一进一小说| 色播在线永久视频| 免费av毛片视频| 国产日本99.免费观看| 精品不卡国产一区二区三区| 俄罗斯特黄特色一大片| xxx96com| 国产一区在线观看成人免费| 久99久视频精品免费| 国语自产精品视频在线第100页| 在线天堂中文资源库| 99久久综合精品五月天人人| 久久婷婷成人综合色麻豆| 国产伦人伦偷精品视频| 免费观看精品视频网站| 国产高清激情床上av| 日韩有码中文字幕| 亚洲欧美精品综合一区二区三区| 视频区欧美日本亚洲| 日本三级黄在线观看| 亚洲av成人不卡在线观看播放网| 十八禁人妻一区二区| 中文资源天堂在线| 天堂动漫精品| 日韩欧美在线二视频| 欧美又色又爽又黄视频| 嫁个100分男人电影在线观看| 精品国产国语对白av| 亚洲成人国产一区在线观看| 老司机午夜福利在线观看视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲成人精品中文字幕电影| 香蕉国产在线看| 国产精品国产高清国产av| 成人免费观看视频高清| 成人18禁高潮啪啪吃奶动态图| av中文乱码字幕在线| 嫁个100分男人电影在线观看| 超碰成人久久| 一二三四社区在线视频社区8| 欧美日韩亚洲国产一区二区在线观看| 亚洲av成人不卡在线观看播放网| 国产亚洲av嫩草精品影院| 久久精品亚洲精品国产色婷小说| 亚洲熟妇中文字幕五十中出| 亚洲avbb在线观看| 神马国产精品三级电影在线观看 | 人妻丰满熟妇av一区二区三区| 国产不卡一卡二| 中文字幕精品亚洲无线码一区 | 国产高清有码在线观看视频 | 国语自产精品视频在线第100页| 国产一区二区三区在线臀色熟女| 亚洲在线自拍视频| avwww免费| xxxwww97欧美| 亚洲国产欧美日韩在线播放| 亚洲五月天丁香| 国产亚洲精品av在线| 99精品欧美一区二区三区四区| 亚洲欧美精品综合久久99| 日本免费一区二区三区高清不卡| 日韩大码丰满熟妇| 少妇熟女aⅴ在线视频| 久久久国产精品麻豆| 成人18禁高潮啪啪吃奶动态图| 99re在线观看精品视频| 日韩一卡2卡3卡4卡2021年| 在线观看www视频免费| 国产精品久久久久久人妻精品电影| 热re99久久国产66热| 久久久久久久久久黄片| 国内久久婷婷六月综合欲色啪| 午夜免费成人在线视频| 中文字幕人妻熟女乱码| 夜夜看夜夜爽夜夜摸| 亚洲一码二码三码区别大吗| 亚洲中文日韩欧美视频| www.999成人在线观看| 亚洲avbb在线观看| 国产成人系列免费观看| 无人区码免费观看不卡| 国产成人欧美| 午夜福利视频1000在线观看| 亚洲精品一卡2卡三卡4卡5卡| 日韩视频一区二区在线观看| 国产精品野战在线观看| 一二三四社区在线视频社区8| 在线播放国产精品三级| 日韩精品免费视频一区二区三区| 日本撒尿小便嘘嘘汇集6| 身体一侧抽搐| www.www免费av| 麻豆国产av国片精品| 精品不卡国产一区二区三区| 国产精品久久久久久人妻精品电影| 黑人操中国人逼视频| 久久精品影院6| 最新美女视频免费是黄的| 99re在线观看精品视频| www国产在线视频色| 午夜福利在线在线| 三级毛片av免费| 亚洲人成网站在线播放欧美日韩| 99久久无色码亚洲精品果冻| av电影中文网址| 亚洲av日韩精品久久久久久密| 高潮久久久久久久久久久不卡| 国产精品久久久久久精品电影 | 老汉色av国产亚洲站长工具| 中文在线观看免费www的网站 | 又大又爽又粗| 国产精品99久久99久久久不卡| 日日干狠狠操夜夜爽| 精品国产国语对白av| 久久香蕉国产精品| 国产99久久九九免费精品| 美女高潮喷水抽搐中文字幕| 制服人妻中文乱码| 香蕉av资源在线| 日韩欧美在线二视频| 国产真人三级小视频在线观看| 国产97色在线日韩免费| 男女午夜视频在线观看| 午夜激情福利司机影院| 久久久久久人人人人人| 在线观看日韩欧美| 少妇粗大呻吟视频| 一级作爱视频免费观看| 国产一区二区三区视频了| 一区二区三区精品91| 在线观看免费视频日本深夜| 久久久久国内视频| 熟女电影av网| 久久精品成人免费网站| 久久国产精品影院| 亚洲精品久久国产高清桃花| 少妇 在线观看| 免费人成视频x8x8入口观看| 在线十欧美十亚洲十日本专区| 欧美日韩乱码在线| 无限看片的www在线观看| 国产精品av久久久久免费| 午夜日韩欧美国产| 亚洲无线在线观看| 99久久99久久久精品蜜桃| 色综合站精品国产| 色播在线永久视频| 极品教师在线免费播放| 国产精品亚洲美女久久久| 俄罗斯特黄特色一大片| 成人永久免费在线观看视频| 丝袜美腿诱惑在线| 老司机深夜福利视频在线观看| 欧美激情高清一区二区三区| 国产伦一二天堂av在线观看| 香蕉丝袜av| 亚洲av中文字字幕乱码综合 | 搞女人的毛片| 国产高清有码在线观看视频 | 亚洲一区二区三区色噜噜| 麻豆av在线久日| 久热爱精品视频在线9| 麻豆av在线久日| 少妇 在线观看| 免费看十八禁软件| 国产精品99久久99久久久不卡| av超薄肉色丝袜交足视频| 亚洲av五月六月丁香网| 久久精品国产亚洲av高清一级| 欧美精品啪啪一区二区三区| 国产精品久久久久久人妻精品电影| 激情在线观看视频在线高清| 老鸭窝网址在线观看| 国产精品免费视频内射| 大香蕉久久成人网| 亚洲成人久久性| 久9热在线精品视频| 12—13女人毛片做爰片一| 在线观看日韩欧美| 久久久久亚洲av毛片大全| 制服诱惑二区| 国产精品1区2区在线观看.| 免费看美女性在线毛片视频| 午夜精品在线福利| 男女视频在线观看网站免费 | 欧美最黄视频在线播放免费| 99久久久亚洲精品蜜臀av| a级毛片在线看网站| 国产成+人综合+亚洲专区| 变态另类丝袜制服| 国产久久久一区二区三区| 欧美日韩黄片免| 国产黄色小视频在线观看| 黄色毛片三级朝国网站| 国产一卡二卡三卡精品| 欧美中文日本在线观看视频| 精品电影一区二区在线| 丝袜人妻中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久午夜电影| 中文字幕最新亚洲高清| 搡老岳熟女国产| 最近最新中文字幕大全免费视频| 又黄又粗又硬又大视频| 88av欧美| 国产成人欧美在线观看| 黄色 视频免费看| 国产成人精品无人区| 精品高清国产在线一区| 淫秽高清视频在线观看| 国产成人影院久久av| 午夜激情av网站| 制服丝袜大香蕉在线| 九色国产91popny在线| 91麻豆精品激情在线观看国产| 国产激情久久老熟女| 欧美又色又爽又黄视频| 国产黄色小视频在线观看| 亚洲精品av麻豆狂野| 男女午夜视频在线观看| xxxwww97欧美| 好男人电影高清在线观看| 午夜免费激情av| 成人国产综合亚洲| 久久久精品欧美日韩精品| 国产又黄又爽又无遮挡在线| 99国产精品一区二区蜜桃av| 国内精品久久久久久久电影| 日韩中文字幕欧美一区二区| 一夜夜www| 黄色毛片三级朝国网站| 精品人妻1区二区| 欧美黑人巨大hd| 亚洲五月色婷婷综合| 波多野结衣高清作品| 国产激情久久老熟女| 亚洲欧美日韩无卡精品| 99久久国产精品久久久| 国产99白浆流出| 丝袜人妻中文字幕| 亚洲一区中文字幕在线| 丝袜美腿诱惑在线| 99国产精品99久久久久| 精品乱码久久久久久99久播| 国产精品 欧美亚洲| 欧美在线黄色| 性欧美人与动物交配| av天堂在线播放| 中文字幕人妻熟女乱码| 国产一级毛片七仙女欲春2 | 久久久久国产精品人妻aⅴ院| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3 | 1024视频免费在线观看| 曰老女人黄片| cao死你这个sao货| 亚洲一区中文字幕在线| 日韩欧美免费精品| 搞女人的毛片| 嫁个100分男人电影在线观看| 999久久久国产精品视频| 少妇的丰满在线观看| 欧美午夜高清在线| 中文字幕高清在线视频| 免费搜索国产男女视频| 精品不卡国产一区二区三区| 国产国语露脸激情在线看| 美女大奶头视频| 国产片内射在线| 国产三级黄色录像| 欧美日韩瑟瑟在线播放| 国产精华一区二区三区| 人人妻人人澡欧美一区二区| 久久久久精品国产欧美久久久| 黄色视频不卡| 日本熟妇午夜| 禁无遮挡网站| 99国产极品粉嫩在线观看| 法律面前人人平等表现在哪些方面| 国产精品一区二区三区四区久久 | 黄色视频不卡| www日本在线高清视频| 久久性视频一级片| 国产野战对白在线观看| 丰满人妻熟妇乱又伦精品不卡| 黄色视频不卡| 婷婷亚洲欧美| 午夜精品久久久久久毛片777| 欧美av亚洲av综合av国产av| 最新美女视频免费是黄的| 精品久久久久久久毛片微露脸| 精品福利观看| 一进一出抽搐动态| 久久精品aⅴ一区二区三区四区| 久久青草综合色| 亚洲av第一区精品v没综合| 亚洲中文字幕日韩| 国产97色在线日韩免费| 黄片小视频在线播放| 亚洲国产日韩欧美精品在线观看 | 亚洲 国产 在线| 亚洲欧美精品综合久久99| 欧美黑人欧美精品刺激| 老司机深夜福利视频在线观看| 男人的好看免费观看在线视频 | 极品教师在线免费播放| 熟女电影av网| 夜夜夜夜夜久久久久| 国产一卡二卡三卡精品| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品久久国产高清桃花| 日本五十路高清| 日韩欧美国产在线观看| 极品教师在线免费播放| 一区二区三区激情视频| 在线av久久热| 深夜精品福利| 亚洲最大成人中文| 中文字幕人妻熟女乱码| 久久亚洲精品不卡| 在线国产一区二区在线| 日韩大码丰满熟妇| 村上凉子中文字幕在线| 成人18禁高潮啪啪吃奶动态图| 亚洲性夜色夜夜综合| 亚洲人成网站在线播放欧美日韩| 夜夜夜夜夜久久久久| 在线免费观看的www视频| av电影中文网址| 成人手机av| 村上凉子中文字幕在线| 我的亚洲天堂| 两个人视频免费观看高清| 身体一侧抽搐| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品国产精品久久久不卡| 精品福利观看| 免费在线观看成人毛片| 亚洲第一青青草原| 不卡av一区二区三区| 久久午夜亚洲精品久久| 亚洲国产毛片av蜜桃av| 日韩大尺度精品在线看网址| 亚洲国产精品成人综合色| 好男人电影高清在线观看| 国产av一区在线观看免费| 黄色视频,在线免费观看| 亚洲 欧美 日韩 在线 免费| 久久精品国产99精品国产亚洲性色| 波多野结衣高清作品| 国产精品国产高清国产av| 人人妻人人澡人人看| 久久精品国产亚洲av香蕉五月| 亚洲av成人一区二区三| 一本久久中文字幕| 亚洲第一av免费看| 国产1区2区3区精品| 视频在线观看一区二区三区| 国产野战对白在线观看| 我的亚洲天堂| 一个人观看的视频www高清免费观看 | 亚洲成人久久爱视频| 久久99热这里只有精品18| 中文资源天堂在线| 香蕉国产在线看| 国产精品香港三级国产av潘金莲| 亚洲中文日韩欧美视频| 一区二区三区精品91| 巨乳人妻的诱惑在线观看| 真人做人爱边吃奶动态| a在线观看视频网站| 午夜福利18| 99精品在免费线老司机午夜| 日韩av在线大香蕉| 亚洲精品在线观看二区| 亚洲在线自拍视频| 99久久无色码亚洲精品果冻| 日本五十路高清| 黄色a级毛片大全视频| 国产精品久久久久久人妻精品电影| 国产精品久久久久久亚洲av鲁大| 日本一本二区三区精品| 色播亚洲综合网| 日本五十路高清| 高潮久久久久久久久久久不卡| 男女视频在线观看网站免费 | 国产伦人伦偷精品视频| 午夜福利高清视频| 国产又黄又爽又无遮挡在线| 成人午夜高清在线视频 | 特大巨黑吊av在线直播 | 国产精品99久久99久久久不卡| av片东京热男人的天堂| 男女床上黄色一级片免费看| 精品福利观看| 性欧美人与动物交配| 国产三级黄色录像| 夜夜看夜夜爽夜夜摸| 校园春色视频在线观看| 俄罗斯特黄特色一大片| 少妇粗大呻吟视频|