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

    典型高危害日化品的篩選及其水生生物預(yù)測無效應(yīng)濃度推導(dǎo)

    2022-09-28 08:07:50張璐瑩于洋黃怡王力明鄭玉婷張麗麗李晨男高子竣林軍于彩虹
    生態(tài)毒理學(xué)報 2022年3期
    關(guān)鍵詞:物種信息

    張璐瑩,于洋,黃怡,王力明,鄭玉婷,張麗麗,李晨男,高子竣,林軍,于彩虹,#

    1.中國礦業(yè)大學(xué)(北京)化學(xué)與環(huán)境工程學(xué)院,北京 100083

    2.生態(tài)環(huán)境部固體廢物與化學(xué)品管理技術(shù)中心,北京 100029

    3.北京市污染源管理事務(wù)中心,北京 100089

    日化品給人們生活帶來極大便利的同時,也對生態(tài)環(huán)境與人體健康造成了潛在的危害。例如,洗滌劑中含有表面活性劑和助劑,其中三聚磷酸鹽(STPP)是最理想的助劑,其在水體富營養(yǎng)化進(jìn)程中,起到了一定的促進(jìn)作用。1970年日本琵琶湖等封閉水域由于水中磷酸鹽含量超過正常值,水體富營養(yǎng)化,出現(xiàn)了水藻瘋長、魚類死亡的現(xiàn)象。20世紀(jì)80年代以后我國湖泊、水庫等水域也普遍受到磷、氮等物質(zhì)的污染,水域生態(tài)系統(tǒng)遭到破壞,魚類和貝類大量死亡,從而影響了人類生活[1]。目前世界上約有7 000多種化學(xué)物質(zhì)作為化妝品的原料,其中最常見的化學(xué)毒物為汞、鉛、砷、甲醇及雌激素等,能引起皮膚功能障礙的原料中最危險的是香料,其次為色素和防腐殺菌劑[2]。這些以原料和添加劑作為化妝品組分的化學(xué)物質(zhì),大多為人工合成,可能對人體產(chǎn)生一定的急、慢性危害,如可使人體過敏、刺激皮膚黏膜等,甚至不乏誘發(fā)人體疾病的物質(zhì),以至潛在的致癌性、致突變性和生殖毒性的物質(zhì)(carcinogenic,mutagenic or substances toxic to reproduction,CMR)類物質(zhì)。據(jù)調(diào)查顯示,染發(fā)劑中均含有20余種化學(xué)成分,其中有10種會引起人體細(xì)胞突變而致癌,這也是導(dǎo)致近年來國內(nèi)中老年急性白血病患者顯著增多的原因之一。

    歐盟化學(xué)品管理法規(guī)(Registration,Evaluation,Authorization and Restriction of Chemicals,REACH),將CMR類物質(zhì)劃分為有毒有害物質(zhì)[3]。部分日化品,例如洗滌劑、清潔用產(chǎn)品、油漆去除劑、香皂和沐浴露等制劑產(chǎn)品中均可能存在含有CMR類物質(zhì)的風(fēng)險,這類物質(zhì)的長期大量使用會對環(huán)境和健康帶來直接或間接的危害。因此,如何對CMR類日化品進(jìn)行篩選評估并加以管控十分重要。

    本研究基于文獻(xiàn)調(diào)研法及數(shù)據(jù)庫檢索法,收集了8類日化品成分清單,包括國際香料組織香水成分清單、國際化妝品原料標(biāo)準(zhǔn)中文名稱目錄清單和食品用洗滌劑原料(成分)清單等。通過與歐洲化學(xué)品管理局(ECHA)最新發(fā)布的高關(guān)注物質(zhì)清單、國際癌癥研究機(jī)構(gòu)(IARC)的致癌物清單、法國國家科學(xué)研究中心(CNRS)的CMR類物質(zhì)清單和內(nèi)分泌干擾物清單進(jìn)行比對篩重,識別出一種典型高危害CMR類物質(zhì)。通過數(shù)據(jù)質(zhì)量評估開展高質(zhì)量數(shù)據(jù)篩選,采用物種敏感度分布法(SSD)構(gòu)建該物質(zhì)對各物種的敏感度曲線,計算出5%物種危害濃度(HC5)值,進(jìn)而應(yīng)用評估系數(shù)法推導(dǎo)出各物種的預(yù)測無效應(yīng)濃度(PNEC)值,評估該污染物對水生生物的危害,以期為我國典型高危害日化品的環(huán)境風(fēng)險評估和水質(zhì)基準(zhǔn)的建立提供參考。

    1 研究方法(Research methodology)

    1.1 高危害日化品清單的建立及篩選

    1.1.1 高危害日化品清單的建立

    本研究采用文獻(xiàn)調(diào)研法以及數(shù)據(jù)庫檢索法,收集各類日化品成分清單為后續(xù)篩選工作提供數(shù)據(jù)支撐。

    1.1.2 高危害日化品清單的篩選

    本研究篩選的是同時具有CMR和內(nèi)分泌干擾物(EDC)特性的,可能會對人體健康帶來潛在危害的高危害日化品。我國任幸等[4]利用基于風(fēng)險的Football組合法對15種農(nóng)用地酞酸酯類內(nèi)源污染物進(jìn)行篩選,發(fā)現(xiàn)鄰苯二甲酸二丁酯(DBP)、鄰苯二甲酸二(2-乙基己基)酯(DEHP)等5種酞酸酯類污染物可能會對生態(tài)環(huán)境和人體健康造成風(fēng)險。朱曉晶等[5]等運用Copeland計分排序法和證據(jù)權(quán)重法對50種有毒有害大氣污染物開展危害特性優(yōu)先性排序。張麗麗等[6]采用兩步篩選法,篩選出高分值的污染物名單,進(jìn)而形成包括苯、環(huán)氧乙烷等在內(nèi)的12種污染物優(yōu)控名錄。以上學(xué)者在篩選高危害污染物時的方法雖然不完全相同,但關(guān)注的篩選指標(biāo)均為致癌性、致突變、生殖毒性和內(nèi)分泌干擾性等,且都是通過名錄對比法得到篩選結(jié)果。

    本研究利用EXCEL篩選重復(fù)值工具進(jìn)行名錄對比,首先將8類日化品成分清單中帶有CAS號的有效物質(zhì)信息與ECHA最新發(fā)布的高關(guān)注物質(zhì)(SVHC)清單和IARC的1類致癌物清單中的有效物質(zhì)信息進(jìn)行比對篩重,得到一個CMR類物質(zhì)的初篩清單。為保證篩選結(jié)果更為準(zhǔn)確,將上述初篩清單與法國國家科學(xué)研究中心化學(xué)研究院的化學(xué)品防范(PRC)管理部門發(fā)布的CMR類物質(zhì)清單對照,重復(fù)以上篩選操作,得到進(jìn)一步明確的物質(zhì)清單。最后,再與內(nèi)分泌干擾物物質(zhì)清單進(jìn)行篩重,識別出一個典型高危害CMR類物質(zhì)。

    1.2 水生生物毒性數(shù)據(jù)的獲取和篩選

    1.2.1 水生生物毒性數(shù)據(jù)的獲取

    本研究從國內(nèi)外公開發(fā)表的文獻(xiàn)和報告中收集典型高危害日化品對水生生物的急、慢性毒性數(shù)據(jù)。當(dāng)同一物種存在多組數(shù)據(jù)時,選擇最敏感試驗終點的數(shù)據(jù)。

    毒性數(shù)據(jù)指標(biāo)包括物種種類、拉丁文名稱、測試標(biāo)準(zhǔn)、介質(zhì)、測試時間、毒性終點、毒性值、數(shù)據(jù)來源和數(shù)據(jù)年份等內(nèi)容。其中,水生生物種類優(yōu)先采用藻類、溞類和魚類3個營養(yǎng)級的毒性數(shù)據(jù)。急性毒性數(shù)據(jù)選擇暴露時間96 h以內(nèi)的LC50或EC50數(shù)據(jù),其中魚類選擇96 h的數(shù)據(jù),溞類選擇48 h的數(shù)據(jù),藻類選擇72 h或96 h的數(shù)據(jù);慢性毒性數(shù)據(jù)選擇暴露時間最長的無可見效應(yīng)濃度(NOEC)以及最低有影響濃度(LOEC)。

    1.2.2 水生生物毒性數(shù)據(jù)的篩選

    本研究應(yīng)用Klimisch等[7]于1997年提出的生態(tài)毒理試驗可靠性評估一級及二級數(shù)據(jù)信息表,對所收集到的毒性數(shù)據(jù)進(jìn)行數(shù)據(jù)質(zhì)量評估。根據(jù)表中21個評估因素對文獻(xiàn)進(jìn)行打分,通過總分?jǐn)?shù)將文獻(xiàn)劃分為非??煽俊⒖煽亢筒豢煽?個等級,總分?jǐn)?shù)在18~21分為非??煽浚?3~17分為可靠;<13分為不可靠,進(jìn)而篩選出可供后續(xù)使用的高質(zhì)量數(shù)據(jù),具體評估信息如表1所示。

    表1 生態(tài)毒理試驗可靠性評估表Table 1 Factors for reliability assessment of ecotoxicological tests

    續(xù)表1數(shù)據(jù)質(zhì)量評估因素Data quality assessment factors組Ⅳ:試驗結(jié)果GroupⅣ: Test results測試終點和其分析方法是否闡述清晰?Are the end points and analysis methods clearly stated?所有測試終點相關(guān)的試驗結(jié)果描述是否透明和完整?Are the descriptions of test results related to all test endpoints transparent and complete?是否給出數(shù)據(jù)分析的統(tǒng)計方法(如不必要或不適用請直接給分)?Are statistical methods for data analysis provided (please give marks directly if unnecessary or not applicable)?組V:試驗設(shè)計和結(jié)果的合理性Group V: Experimental design and rationality of results*實驗設(shè)計對于旨在獲得的數(shù)據(jù)是否合適?*Is the experimental design appropriate for the data to be obtained?定量測試方法是否可靠Whether the quantitative test method is reliable備注Remarks1. 每條標(biāo)準(zhǔn)賦分分值均為1分Each criterion is assigned a score of 1 point2. *代表達(dá)到1或2級評分的最大分?jǐn)?shù)來源* Represents the maximum source of points to achieve a score of 1 or 23. 18~21分為非??煽?13~17分為可靠;<13分或若不滿足所有的*標(biāo)準(zhǔn)為不可靠18~21 score is very reliable; 13~17 is reliable; <13 points or unreliable if all * criteria are not met

    對于收集到的所有生態(tài)毒理學(xué)終點數(shù)據(jù),首先剔除存在支撐信息不充分、測試明顯不規(guī)范等情形的無效數(shù)據(jù)。對于有效數(shù)據(jù)中的所有生態(tài)毒理學(xué)數(shù)據(jù)進(jìn)行質(zhì)量評估,篩選滿足質(zhì)量要求的數(shù)據(jù)。選用最敏感物種數(shù)據(jù)作為關(guān)鍵效應(yīng)數(shù)據(jù),明確評估高危害典型日化品對水環(huán)境生物的危害性。

    1.3 HC5值的計算

    本研究通過數(shù)據(jù)庫檢索法和文獻(xiàn)調(diào)研法收集了典型高危害日化品的毒性數(shù)據(jù),利用數(shù)據(jù)質(zhì)量評估開展高質(zhì)量數(shù)據(jù)篩選,采用SSD法構(gòu)建該物質(zhì)對全部物種、魚類、藻類和浮游動物的敏感度曲線,計算出3個類別水生生物的HC5值。

    SSD曲線的擬合采用國家生態(tài)環(huán)境基準(zhǔn)委員會于2021年研制開發(fā)的“國家生態(tài)環(huán)境基準(zhǔn)計算軟件 物種敏感度分布法”(EEC-SSD)來完成,以此計算出HC5。EEC-SSD主要用于毒性數(shù)據(jù)的物種敏感度分布擬合,軟件涵蓋的擬合模型包括正態(tài)分布、對數(shù)正態(tài)分布、邏輯斯蒂分布和對數(shù)邏輯斯蒂分布4個模型。由于HC5更具有統(tǒng)計學(xué)意義,故一般選用HC5作為最大有害環(huán)境濃度[8]。在數(shù)據(jù)量上的選擇,一般認(rèn)為控制數(shù)據(jù)量在10~15個隨機(jī)選擇量就能符合統(tǒng)計分析的要求,此時數(shù)據(jù)變異性較小。這些數(shù)據(jù)常要求覆蓋3個營養(yǎng)級8個分類群,并選用標(biāo)準(zhǔn)毒性測試物種,數(shù)據(jù)過少(少于5個)或是數(shù)據(jù)質(zhì)量差(多個物種的毒性數(shù)據(jù)相同),都可能導(dǎo)致擬合度的降低。另外,數(shù)據(jù)并非越多越好,數(shù)據(jù)越多HC5就會越小,這是不合理的[9]。

    1.4 預(yù)測無效應(yīng)濃度的推導(dǎo)

    水環(huán)境危害表征通常采用統(tǒng)計外推法或評估系數(shù)法,推導(dǎo)水環(huán)境中生物的PNEC。對于不同方法獲得的不同PNEC值,應(yīng)予以比較分析,確定用于適合水環(huán)境危害表征的PNEC值。

    本研究獲得的水生生物毒性數(shù)據(jù)充分且滿足評估系數(shù)法對數(shù)據(jù)的基本要求,故采用評估系數(shù)法推導(dǎo)PNEC值。該方法由危害識別確定的水環(huán)境生物的生態(tài)毒理學(xué)關(guān)鍵效應(yīng)數(shù)據(jù)(HC5)除以評估系數(shù)(AFwater)推算得出PNEC值,由于數(shù)據(jù)來自實驗室,實驗物種單一,具有較好的保守性。具體見公式(1)。

    PNECwater=HC5/AFwater

    (1)

    式中:PNECwater為水環(huán)境生物的預(yù)測無效應(yīng)濃度(mg·L-1);HC5為水環(huán)境生物的生態(tài)毒理學(xué)關(guān)鍵效應(yīng)值,通常采用最敏感物種的半數(shù)致死濃度(LC50)、效應(yīng)濃度(EC50或EC10)或無觀察效應(yīng)濃度(NOEC)及最低有影響濃度(LOEC)等(mg·L-1);AFwater為評估系數(shù),取值范圍為1~5。AFwater的取值根據(jù)可獲得毒性數(shù)據(jù)的質(zhì)量、涉及生物物種的代表性和多樣性等情況選擇確定,無量綱。根據(jù)生態(tài)風(fēng)險評價中“最壞情況”(The Worstcase)原則,該研究中AFwater值取5。

    2 結(jié)果與分析(Results and analysis)

    2.1 初篩清單的建立及篩選

    通過數(shù)據(jù)庫檢索法收集了8類日化品成分清單,共28 515條物質(zhì)信息,僅選用具有CAS號的15 201條有效物質(zhì)信息(以下簡稱有效信息)供后續(xù)名錄對比使用,如表2所示。其中包括國際香料組織(IFRA)香水成分清單共3 963條物質(zhì)信息;國際化妝品原料標(biāo)準(zhǔn)中文名稱目錄2018(草)清單共24 039條物質(zhì)信息,其中10 741條有效信息;食品用洗滌劑原料(成分)名單(第一批)(草)共122條信息,其中120條有效信息;抗(抑)菌劑有效成分清單2019(草)96條物質(zhì)信息,其中95條有效信息;消毒劑原料清單及禁限用物質(zhì)名單共203條物質(zhì)信息,其中190條有效信息;通過查詢口腔衛(wèi)生用品成分表、洗發(fā)用品成分表、洗衣用品成分表,形成了3類日化品的成分清單,共92條物質(zhì)信息,均為有效信息,如表2所示。

    表2 日化品清單初篩表Table 2 Preliminary screening list of daily chemical products

    ECHA最新的高關(guān)注物質(zhì)(SVHC)清單共440條物質(zhì)信息,其中390條有效信息;IARC致癌物清單共1 112條物質(zhì)信息,其中包括1類、2A類、2B類以及3類致癌物有效信息共1 029條,僅選取1類致癌物的121條物質(zhì)信息,其中60條有效信息;法國國家科學(xué)研究中心化學(xué)研究院的化學(xué)品防范(PRC)管理部門CMR類物質(zhì)清單共1 238條物質(zhì)信息;內(nèi)分泌干擾物物質(zhì)清單共688條物質(zhì)信息,其中620條有效信息,如表3所示。

    表3 4類篩選條件物質(zhì)信息表Table 3 Material information sheet of 4 categories of screening conditions

    將8類日化品成分清單與SVHC清單和IARC的1類致癌物清單中的有效物質(zhì)信息進(jìn)行比對篩重,得到一個CMR類物質(zhì)的初篩清單。為保證篩選結(jié)果更為準(zhǔn)確,將上述初篩清單與法國CNRS的CMR類物質(zhì)清單對照重復(fù)以上篩選操作,得到進(jìn)一步明確的物質(zhì)清單。最后,再與內(nèi)分泌干擾物物質(zhì)清單進(jìn)行篩重,最終得到鄰苯二甲酸丁芐酯(BBP)、1,2-二氯乙烷、N-甲基吡咯烷酮、鄰苯二甲酸二(2-乙基己)酯(DEHP)、乙二醇甲醚、鄰苯二甲酸二丁酯(DBP)和酚酞7種1類CMR類物質(zhì),如圖1所示。

    圖1 典型高危害日化品篩選流程圖Fig. 1 Screening flow chart of typical high-hazard daily chemical products

    對以上7種CMR類物質(zhì)近10年的文獻(xiàn)發(fā)表數(shù)量進(jìn)行統(tǒng)計,結(jié)果表明有關(guān)酚酞的文獻(xiàn)發(fā)表量最高,但其環(huán)境風(fēng)險評估方面的文獻(xiàn)量幾乎為零,或因為其環(huán)境濃度較低或不造成生態(tài)危害。其次是鄰苯二甲酸酯類(PAEs),包括BBP、DEHP和DBP這3種物質(zhì),相對于BBP和DEHP,DBP在總文獻(xiàn)發(fā)表量與環(huán)境風(fēng)險評估方面的文獻(xiàn)發(fā)表量都較高,如圖2所示。又因其是一種新污染物,將是今后的研究熱點及重點治理管控對象之一,故本文選擇DBP作為研究對象。

    圖2 7種致癌、致突變和致生殖毒性(CMR)類物質(zhì)近10年的文獻(xiàn)發(fā)表數(shù)量Fig. 2 Number of literature publications for seven carcinogenic,mutagenic or toxic to reproduction (CMR)-like substances in the last decade

    2.2 數(shù)據(jù)質(zhì)量評估

    通過文獻(xiàn)調(diào)研法收集到23篇相關(guān)文獻(xiàn),其中,中文文獻(xiàn)10篇,英文文獻(xiàn)13篇,共計151條數(shù)據(jù),根據(jù)所需毒性數(shù)據(jù)信息進(jìn)一步在文獻(xiàn)中提取有效信息。通過數(shù)據(jù)庫檢索法分別從ECHA數(shù)據(jù)庫[10]、EPA-ECOTOX數(shù)據(jù)庫[11]中收集到共計24條數(shù)據(jù),均視為有效數(shù)據(jù)。通過對比表1中21項評估因素,對該23篇文獻(xiàn)進(jìn)行數(shù)據(jù)質(zhì)量評估,得到非??煽课墨I(xiàn)10篇(18~21分),可靠文獻(xiàn)5篇(13~17分),不可靠文獻(xiàn)8篇(<13分)。上述文獻(xiàn)共得到52條相關(guān)數(shù)據(jù),通過數(shù)據(jù)質(zhì)量評估最終得到17條有效數(shù)據(jù),而滿足數(shù)據(jù)質(zhì)量要求的數(shù)據(jù)庫數(shù)據(jù)24條,合計41條有效數(shù)據(jù)。

    根據(jù)PNEC值對數(shù)據(jù)的要求,對41條評估數(shù)據(jù)進(jìn)一步篩選,最終篩選出滿足生物毒性數(shù)據(jù)要求的17個物種的急性數(shù)據(jù),如表4所示。17個物種分為三大類,即魚類、藻類和浮游動物。其中魚類包含4種即三角魴(Megalobramaterminalis)、斑馬魚(Daniorerio)、紅鰭笛鯛(Lutjanuserythropterus)和大黃魚(Larimichthyscrocea);藻類包含9種即斜生柵藻(Scenedesmusobliquus)、埃氏小球藻(C.emersonii)、蛋白核小球藻(C.pyrenoidosa)、近具刺柵藻(S.subspicatus)、羊角月牙藻(Selenastrumcapricornutum)、短裸甲藻(Gymnodiniumbreve)、巴夫杜氏藻(D.parva)、蒙古裸腹藻(MoinamongolicaDaday)和三角褐指藻(PhaeodactylumtricornutumBohlin);浮游動物包含2種即多刺裸腹溞(Moinamacrocopa)、大型溞(Daphniamagna)。全部物種LC50/EC50的范圍為0.28~13 mg·L-1,魚類LC50/EC50的范圍為2.08~8.51 mg·L-1,藻類LC50/EC50的范圍為0.28~13 mg·L-1,浮游動物L(fēng)C50/EC50的范圍為9.68~10.35 mg·L-1。

    2.3 HC5值的計算

    本研究分別通過正態(tài)分布模型和邏輯斯蒂分布模型,對17個水生生物物種的毒性數(shù)據(jù)進(jìn)行擬合及分析,得到基于急性毒性的物種敏感度分析曲線,如圖3所示。由于SSD模型要求每組樣本數(shù)據(jù)量在5~500個范圍內(nèi),而浮游動物的數(shù)據(jù)量僅有2個,不滿足SSD曲線的要求,因此未對浮游動物擬合SSD曲線。

    生物分類Biological classification物種中文名/物種拉丁名Chinese name of the species/Latin name of the species毒性值Toxicity valueLC50/EC50/(mg·L-1)試驗時間/hTime duration/h數(shù)據(jù)來源Data sources魚類Fish三角魴 Megalobrama terminalis2.0896[12]紅鰭笛鯛 Lutjanus erythropterus6.6696[13] 虹鱒 Salmo gairdneri6.4796[14] 斑點叉尾 Ictalurus punctatus2.9196[15] 斑馬魚 Brachydanio rerio8.5196[15] 大黃魚 Larimichthys crocea5.2396[16] 藻類Alga斜生柵藻 Scenedesmus obliquus2.2196[17] 埃氏小球藻 C. emersonii2.7896[18] 蛋白核小球藻 C. pyrenoidosa1396[19] 近具刺柵藻 S. subspicatus3.572[20] 羊角月牙藻 Selenastrum capricornutum0.472[18] 短裸甲藻 Gymnodinium breve0.696[21] 巴夫杜氏藻 D. parva0.2872[22] 蒙古裸腹藻 Moina mongolica Daday6.3472[23] 三角褐指藻 Phaeodactylum tricornutum Bohlin7.596[24] 浮游動物Zooplankton多刺裸腹溞 Moina macrocopa9.6848[25] 大型溞 Daphnia magna10.3548[26]

    圖3 鄰苯二甲酸二丁酯(DBP)對不同物種的物種敏感度分布(SSD)曲線Fig. 3 Species sensitivity distribution (SSD) curves of dibutyl phthalate (DBP) for different species

    毒性值的高低對DBP在不同物種體內(nèi)的累積頻率具有一定影響。如圖3所示,當(dāng)毒性值較低時(<0.79 mg·L-1),不同物種對DBP的累積頻率由高到低的順序依次為:藻類>全部物種>魚類,表明藻類對DBP的敏感度高于其他物種;當(dāng)毒性值在0.79~0.90 mg·L-1之間時,不同物種對DBP的累積概率由高到低的順序依次為:藻類>魚類>全部物種,藻類的敏感度仍處于最高,但魚類的敏感度有所變化,高于了全部物種;當(dāng)毒性值>0.90 mg·L-1時,不同物種對DBP的累積頻率由高到低的順序依次為:魚類>全部物種>藻類,表明魚類對DBP的敏感度高于藻類,全部物種則上升到第二位??梢姰?dāng)DBP的毒性值濃度較低時,在藻類體內(nèi)的累積頻率較高,但當(dāng)其毒性值逐漸升高時,其在魚類體內(nèi)的累積頻率逐漸升高。綜上結(jié)果表明,當(dāng)毒性值較高時,相對藻類來說,DBP在魚類體內(nèi)更容易累積。

    采用正態(tài)分布模型和邏輯斯蒂分布模型擬合出SSD曲線,可見以上2種方法獲得了較為一致的結(jié)果。其中,正態(tài)分布模型中全部物種的HC5值為0.59 mg·L-1,邏輯斯蒂模型中的HC5值為0.63 mg·L-1;正態(tài)分布模型中魚類的HC5值為1.97 mg·L-1,邏輯斯蒂模型中的HC5值為2.08 mg·L-1;正態(tài)分布模型中藻類的HC5值為0.42 mg·L-1,邏輯斯蒂模型中的HC5值0.23 mg·L-1。SSD曲線擬合參數(shù)如表5和表6所示。全部物種的擬合度(R2)始終保持在0.9以上;魚類的R2均在0.8以上,但未超過0.9;藻類R2處于0.85~0.92之間。

    表5 正態(tài)分布模型模擬SSD曲線關(guān)鍵擬合參數(shù)Table 5 Key fitting parameters of SSD curve simulated by normal distribution model

    2.4 預(yù)測無效應(yīng)濃度

    根據(jù)公式(1)通過HC5及AFwater推導(dǎo)PNEC。由于正態(tài)分布模型和邏輯斯蒂分布模型擬合SSD曲線,可得出2組包含全部物種、魚類和藻類的HC5值,基于保守原則,本研究HC5選擇2組數(shù)據(jù)中較低的值。PNEC結(jié)果如表7所示。

    表7 不同種類水生生物的預(yù)測無效應(yīng)濃度(PNEC)推算結(jié)果Table 7 Predicted no effect concentration (PNEC) results of different aquatic organisms

    如表6所示,不同物種的PNEC由低到高順序依次為:藻類<全部物種<魚類。在2類不同物種之間,藻類的PNEC最小,數(shù)值為0.046 mg·L-1,魚類的PNEC最大,數(shù)值為0.394 mg·L-1。PNEC結(jié)果表明,DBP對不同物種的PNEC存在差異,魚類與全部物種的PNEC較接近,與藻類的PNEC相差較大。

    表6 邏輯斯蒂分布模型模擬SSD曲線關(guān)鍵擬合參數(shù)Table 6 Key fitting parameters of SSD curve simulated by Logistic distribution model

    3 結(jié)論與討論(Conclusion and discussion)

    3.1 結(jié)論

    本研究通過收集國際香料組織(IFRA)香水成分清單在內(nèi)的8類日化品成分清單,共得到15 201條有效物質(zhì)信息。將以上物質(zhì)信息與高關(guān)注物清單、致癌物清單以及內(nèi)分泌干擾物清單等進(jìn)行比對篩重,得到鄰苯二甲酸丁芐酯(BBP)、1,2-二氯乙烷、N-甲基吡咯烷酮和鄰苯二甲酸二丁酯(DBP)在內(nèi)的7種1類CMR類物質(zhì)。通過文獻(xiàn)檢索,統(tǒng)計結(jié)果表明DBP在總發(fā)文量與環(huán)境風(fēng)險評估方面的發(fā)文量都相對較高,又因其是一種新污染物,故選擇DBP作為本文研究對象。

    通過數(shù)據(jù)質(zhì)量評估,篩除了不可靠文獻(xiàn),得到15篇高質(zhì)量文獻(xiàn),最終獲得包括數(shù)據(jù)庫有效數(shù)據(jù)在內(nèi)的17個物種的急性毒性數(shù)據(jù),共計41條物質(zhì)信息??梢姅?shù)據(jù)質(zhì)量評估對于推導(dǎo)典型有毒有害化學(xué)物質(zhì)的PNEC具有十分重要的意義。

    正態(tài)分布模型和邏輯斯諦分布模型擬合結(jié)果表明,2種模型擬合出的SSD曲線變化趨勢相似,即DBP對水生生物的HC5值較為接近。全部物種的PNEC值為0.118 mg·L-1、魚類的PNEC值為0.394 mg·L-1、藻類的PNEC值為0.046 mg·L-1。全部物種與魚類的PNEC值相近,與藻類存在較大差異,得出DBP對藻類的敏感度最高,對魚類的敏感度相對較低。

    3.2 討論

    本研究使用物種敏感度分布法基于EEC-SSD模型推導(dǎo)出DBP的PNEC值,所用數(shù)據(jù)涉及了魚類、藻類和浮游動物3個營養(yǎng)層級生物在內(nèi)的17個物種,減少了單獨使用某一類生物數(shù)據(jù)推導(dǎo)的不確定性。且相比于沒有開展數(shù)據(jù)質(zhì)量評估所獲得毒性數(shù)據(jù)的研究,本文篩除了不可靠數(shù)據(jù),避免了可能出現(xiàn)的偏差。

    不同的SSD分析軟件對典型有毒有害化學(xué)物質(zhì)的PNEC擬合結(jié)果可能存在差異。我國有學(xué)者利用荷蘭公共健康與環(huán)境研究所(RIVM)開發(fā)的ETX 2.1風(fēng)險評估軟件進(jìn)行SSD分析,進(jìn)而推導(dǎo)出代森錳鋅、百菌清和敵敵畏3種農(nóng)藥水生態(tài)系統(tǒng)的PNEC[27]。該模型主要為對數(shù)正態(tài)分布,具有A-D初級檢驗、K-S中級檢驗和C-M高級檢驗3種階梯式檢驗方法。而本研究采用的EEC-SSD模型,檢驗方法較為單一,僅有K-S檢驗。因此,ETX2.1軟件可能具有更高的擬合度,結(jié)果可能更為精確,但其擬合曲線表現(xiàn)形式較為單一,無法更直觀形象地體現(xiàn)出不同物種之間的種平均毒性與累積頻率的關(guān)系。

    建立水質(zhì)基準(zhǔn)的目的在于保護(hù)水生生態(tài)系統(tǒng)中的生物免受污染。通過對比本研究不同水生生物物種的PNEC值:PNEC藻類

    本研究通過文獻(xiàn)檢索發(fā)現(xiàn),近10年國內(nèi)外對于DBP的毒性研究實驗相對較少,研究主要集中在20世紀(jì)80年代,且試驗物種類型有限,尚未有如林氏細(xì)鯽等我國特有物種的毒性數(shù)據(jù)。另外,本研究獲得的數(shù)據(jù)僅基于文獻(xiàn)調(diào)研與數(shù)據(jù)庫檢索,后續(xù)隨著更多毒性效應(yīng)研究的開展,能夠獲得更多的水生生物毒性數(shù)據(jù)或中國特有物種試驗數(shù)據(jù)時,本研究推算的PNEC值應(yīng)進(jìn)一步更新與完善。

    4 展望(Prospect)

    當(dāng)前,我國新污染物治理任務(wù)艱巨。國家發(fā)布了《中共中央關(guān)于制定國民經(jīng)濟(jì)和社會發(fā)展第十四個五年規(guī)劃和二〇三五年遠(yuǎn)景目標(biāo)的建議》等一系列文件,其中反復(fù)強(qiáng)調(diào)重視新污染物治理,并對此提出明確要求,可見新污染物治理將是未來“十四五”期間的熱點話題。

    新污染物是指由人類活動造成的、目前已明確存在、但尚無法律法規(guī)和標(biāo)準(zhǔn)予以規(guī)定或規(guī)定不完善、危害生活和生態(tài)環(huán)境的所有在生產(chǎn)建設(shè)或者其他活動中產(chǎn)生的污染物[28]。本研究中的鄰苯二甲酸二丁酯(DBP)作為一種典型新污染物,在環(huán)境中廣泛存在。最主要的危害在于環(huán)境激素作用,在極低的濃度下就可干擾人和動物的內(nèi)分泌系統(tǒng)[29]。美國環(huán)境保護(hù)局(United States Environmental Protection Agency,US EPA)于1977年將DBP類列為優(yōu)先控制的有毒污染物,同年世界野生動物基金會(World Wide Fund for Nature,WWF)列出的68種環(huán)境激素類物質(zhì)中也包括DBP[30]。目前,中國已將DBP列入優(yōu)先控制污染物名單[31],《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)也將DBP作為鄰苯二甲酸酯類污染物的檢測指標(biāo),并規(guī)定其含量不得超過3 μg·L-1。

    在環(huán)境危害方面,近年來我國農(nóng)膜的使用量不斷增加,由于塑料農(nóng)膜的穩(wěn)定性差,極易導(dǎo)致DBP從塑料薄膜中釋放。因此,越來越多的DBP進(jìn)入到周邊土壤中,導(dǎo)致農(nóng)業(yè)用地污染。另外,由于化肥、污灌等的不當(dāng)使用,也使土壤中DBP的含量有所增加。受DBP污染土壤的影響,導(dǎo)致種植的作物對DBP的積累,進(jìn)而通過食物鏈進(jìn)入到人體中,對人體內(nèi)分泌、神經(jīng)系統(tǒng)等產(chǎn)生危害[32]。不僅如此,DBP還會以吸附或者氣溶膠的形式存在于氣體環(huán)境中,借助灰塵等進(jìn)入大氣[33]。由于其化學(xué)性質(zhì)穩(wěn)定,不易光解,因此會造成污染的積累[34],并且距離地面越高的空氣中,含量越高[35]。

    在環(huán)境暴露方面,我國水環(huán)境中DBP含量嚴(yán)重超標(biāo),污染嚴(yán)重。隨著塑料的生產(chǎn)、應(yīng)用與處理,DBP不斷地向周圍環(huán)境釋放,水環(huán)境中殘余量持續(xù)增長,已經(jīng)在長江、黃河等重要水域及多種水生生物(海洋哺乳動物、魚類和水生無脊椎動物)中普遍檢出,且對水生生物產(chǎn)生影響[36],還可通過食物鏈的富集和傳遞,潛在威脅人體健康[37-39]??钻猾h和劉紅玲[40]在華東地區(qū)檢出DBP平均濃度為54.389 μg·L-1;卓麗等[41]在黃河中下游檢出DBP濃度為ND~26.0 μg·L-1,胡衛(wèi)星等[42]在渭河流域西安段檢出DBP平均濃度為13.73 μg·L-1;朱鵬利等[43]在曹娥江紹興段檢出DBP濃度為1 974.4 ng·L-1;吳自清等[44]在長江南京段檢出DBP濃度為10~500 μg·L-1;韓文輝等[45]在汾河流域檢出DBP濃度為1.01~16.53 μg·L-1。

    本文推導(dǎo)了新污染物DBP的PNEC值,接下來還需進(jìn)一步開展其暴露評估,識別其潛在環(huán)境風(fēng)險,從而為我國新污染物治理提供參考。同時,現(xiàn)階段我國應(yīng)加強(qiáng)對CMR類日化品風(fēng)險評估及管控等領(lǐng)域的研究。由于日化品種類多樣,在日常生活中暴露水平也相對較高,若其中含有CMR類物質(zhì),其毒性更應(yīng)引起重視,需要盡量降低其最大的暴露源,尤其是化妝品、洗滌劑以及增塑劑中CMR類物質(zhì)的含量。因此,建議各生產(chǎn)企業(yè)和監(jiān)管部門針對日化品的生產(chǎn)線、包裝材料等產(chǎn)業(yè)鏈各環(huán)節(jié)進(jìn)行分級排查,有效管控日化品中CMR物質(zhì)的污染風(fēng)險。

    猜你喜歡
    物種信息
    物種大偵探
    物種大偵探
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    回首2018,這些新物種值得關(guān)注
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    訂閱信息
    中華手工(2017年2期)2017-06-06 23:00:31
    瘋狂的外來入侵物種
    展會信息
    中外會展(2014年4期)2014-11-27 07:46:46
    信息
    健康信息
    祝您健康(1987年3期)1987-12-30 09:52:32
    99热这里只有是精品50| 麻豆成人av视频| 亚洲无线观看免费| 午夜精品一区二区三区免费看| 久久精品国产亚洲网站| 久久人人爽av亚洲精品天堂 | 男女那种视频在线观看| 日日摸夜夜添夜夜爱| 精品久久久久久电影网| 日韩欧美精品v在线| 街头女战士在线观看网站| 欧美日韩在线观看h| 亚洲一区二区三区欧美精品 | 国产精品人妻久久久影院| 国内精品美女久久久久久| 一级黄片播放器| 99热网站在线观看| 亚洲成人中文字幕在线播放| 七月丁香在线播放| 欧美日本视频| 中文字幕亚洲精品专区| 男男h啪啪无遮挡| 亚洲国产av新网站| 九九爱精品视频在线观看| 亚洲av电影在线观看一区二区三区 | 亚洲经典国产精华液单| 欧美极品一区二区三区四区| 丝袜脚勾引网站| 亚洲精品国产av成人精品| 激情 狠狠 欧美| 免费观看a级毛片全部| 中文天堂在线官网| 久久精品人妻少妇| 夜夜看夜夜爽夜夜摸| 特级一级黄色大片| 日日摸夜夜添夜夜添av毛片| 亚洲四区av| 成人黄色视频免费在线看| 成人毛片a级毛片在线播放| 免费观看在线日韩| 久热久热在线精品观看| 听说在线观看完整版免费高清| 街头女战士在线观看网站| 两个人的视频大全免费| 国产精品国产av在线观看| av国产久精品久网站免费入址| 日日撸夜夜添| 亚洲成人av在线免费| 国产精品国产av在线观看| 一级毛片aaaaaa免费看小| 又爽又黄a免费视频| 欧美极品一区二区三区四区| 九色成人免费人妻av| 免费看a级黄色片| 日韩人妻高清精品专区| 黄片无遮挡物在线观看| 国产老妇女一区| 国产人妻一区二区三区在| 欧美日韩在线观看h| 欧美xxxx性猛交bbbb| 在现免费观看毛片| 一本久久精品| 一级毛片久久久久久久久女| 亚洲av.av天堂| 22中文网久久字幕| 伊人久久精品亚洲午夜| 97热精品久久久久久| 在线看a的网站| 夜夜看夜夜爽夜夜摸| 精品亚洲乱码少妇综合久久| 在线亚洲精品国产二区图片欧美 | 在线精品无人区一区二区三 | 午夜福利视频1000在线观看| 美女脱内裤让男人舔精品视频| 干丝袜人妻中文字幕| 只有这里有精品99| 国产欧美亚洲国产| 毛片女人毛片| 欧美亚洲 丝袜 人妻 在线| 午夜福利在线观看免费完整高清在| 99精国产麻豆久久婷婷| 赤兔流量卡办理| 免费黄网站久久成人精品| 春色校园在线视频观看| 亚洲av一区综合| 亚洲欧美精品自产自拍| 中文字幕久久专区| 一本久久精品| 国产欧美另类精品又又久久亚洲欧美| 日韩欧美精品v在线| 2022亚洲国产成人精品| 国产精品一区二区三区四区免费观看| 美女被艹到高潮喷水动态| 国产黄片视频在线免费观看| 最近最新中文字幕免费大全7| 天天躁夜夜躁狠狠久久av| 欧美人与善性xxx| 一级av片app| 国产精品久久久久久av不卡| 王馨瑶露胸无遮挡在线观看| 国国产精品蜜臀av免费| av免费在线看不卡| 国产老妇女一区| 三级经典国产精品| 毛片女人毛片| 99九九线精品视频在线观看视频| 黑人高潮一二区| 久久99蜜桃精品久久| 两个人的视频大全免费| 午夜福利在线观看免费完整高清在| 亚洲色图av天堂| 可以在线观看毛片的网站| 亚洲精品影视一区二区三区av| 3wmmmm亚洲av在线观看| 99热国产这里只有精品6| 日本黄大片高清| 一个人看视频在线观看www免费| 秋霞伦理黄片| 成人综合一区亚洲| 国产精品久久久久久av不卡| 日日摸夜夜添夜夜爱| 又粗又硬又长又爽又黄的视频| 亚洲国产欧美人成| 国产国拍精品亚洲av在线观看| 亚洲天堂av无毛| 一级毛片aaaaaa免费看小| 欧美性猛交╳xxx乱大交人| 亚洲av免费在线观看| 中文天堂在线官网| 日韩人妻高清精品专区| av在线蜜桃| 亚洲丝袜综合中文字幕| 性插视频无遮挡在线免费观看| 身体一侧抽搐| 亚洲国产最新在线播放| 丰满少妇做爰视频| 欧美激情在线99| 欧美精品国产亚洲| 一级毛片电影观看| 亚洲av.av天堂| 交换朋友夫妻互换小说| 欧美+日韩+精品| 观看免费一级毛片| 免费少妇av软件| 伦精品一区二区三区| 欧美区成人在线视频| 日本午夜av视频| 亚洲国产欧美人成| 国产美女午夜福利| 亚洲天堂国产精品一区在线| 日韩在线高清观看一区二区三区| 国产精品一区二区三区四区免费观看| 国产真实伦视频高清在线观看| 久久久精品94久久精品| 久久国内精品自在自线图片| 久久精品国产a三级三级三级| 秋霞伦理黄片| 国产亚洲5aaaaa淫片| 亚洲av免费在线观看| 男男h啪啪无遮挡| 久久99热6这里只有精品| 成人欧美大片| 两个人的视频大全免费| 亚洲欧洲日产国产| 伦精品一区二区三区| 国产男人的电影天堂91| 精品视频人人做人人爽| 精品一区二区三卡| 国产精品三级大全| 日韩中字成人| 久久久久久九九精品二区国产| 人妻一区二区av| 国产在线男女| 国产精品99久久99久久久不卡 | 国产精品99久久99久久久不卡 | 欧美精品人与动牲交sv欧美| 亚洲精品乱久久久久久| 好男人在线观看高清免费视频| 99热这里只有是精品50| 国产乱人视频| 在线观看三级黄色| 亚洲四区av| 国产精品人妻久久久久久| 久久精品人妻少妇| 男女边吃奶边做爰视频| 免费观看的影片在线观看| 女的被弄到高潮叫床怎么办| 男女那种视频在线观看| 亚洲精品aⅴ在线观看| 97超视频在线观看视频| 日韩亚洲欧美综合| 中文字幕人妻熟人妻熟丝袜美| 大话2 男鬼变身卡| 中文精品一卡2卡3卡4更新| 女人十人毛片免费观看3o分钟| 美女高潮的动态| 99热这里只有精品一区| 精品熟女少妇av免费看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久热久热在线精品观看| 永久免费av网站大全| 爱豆传媒免费全集在线观看| 国产精品精品国产色婷婷| 成人二区视频| 三级经典国产精品| 午夜福利网站1000一区二区三区| 91在线精品国自产拍蜜月| 免费看日本二区| 搞女人的毛片| 国产精品一及| 精品久久久精品久久久| 免费黄网站久久成人精品| 岛国毛片在线播放| 视频区图区小说| 国产高清国产精品国产三级 | 免费观看的影片在线观看| 涩涩av久久男人的天堂| 美女脱内裤让男人舔精品视频| 九九久久精品国产亚洲av麻豆| 欧美成人a在线观看| 97精品久久久久久久久久精品| 亚洲av中文字字幕乱码综合| av播播在线观看一区| 69人妻影院| 成年av动漫网址| 国产黄片美女视频| 精品久久久久久电影网| 卡戴珊不雅视频在线播放| 中文字幕制服av| 日本熟妇午夜| 久久久久精品久久久久真实原创| 91精品伊人久久大香线蕉| 丝瓜视频免费看黄片| 内地一区二区视频在线| 最近中文字幕2019免费版| www.色视频.com| 色婷婷久久久亚洲欧美| 成人美女网站在线观看视频| 精品国产三级普通话版| 在线免费十八禁| 天堂网av新在线| 看十八女毛片水多多多| 精品久久久久久电影网| 婷婷色综合www| 国产高清三级在线| 亚洲精品日韩在线中文字幕| 嫩草影院新地址| 五月天丁香电影| 别揉我奶头 嗯啊视频| 日韩欧美一区视频在线观看 | 成年人午夜在线观看视频| 波多野结衣巨乳人妻| 国产乱人视频| 视频中文字幕在线观看| 久久久久久久久久久丰满| 你懂的网址亚洲精品在线观看| 97超视频在线观看视频| 国产黄色视频一区二区在线观看| 亚洲人成网站在线观看播放| 免费高清在线观看视频在线观看| av在线播放精品| av线在线观看网站| 精品国产一区二区三区久久久樱花 | 国产探花在线观看一区二区| 久久热精品热| 偷拍熟女少妇极品色| 免费看av在线观看网站| 成人免费观看视频高清| 久久久久久久久大av| 久久ye,这里只有精品| 黄片wwwwww| 在线播放无遮挡| 久久99热6这里只有精品| av女优亚洲男人天堂| 久久国内精品自在自线图片| 人妻少妇偷人精品九色| 成年免费大片在线观看| 新久久久久国产一级毛片| 久久久久国产网址| 18+在线观看网站| 久热这里只有精品99| 插逼视频在线观看| 亚洲精品自拍成人| 久久久久精品久久久久真实原创| 成人特级av手机在线观看| 免费观看性生交大片5| 人妻少妇偷人精品九色| 欧美丝袜亚洲另类| 欧美极品一区二区三区四区| 午夜免费男女啪啪视频观看| 免费av不卡在线播放| a级一级毛片免费在线观看| 亚洲欧美一区二区三区黑人 | 久久精品国产亚洲av涩爱| 亚洲天堂av无毛| 亚洲欧美日韩另类电影网站 | a级一级毛片免费在线观看| 只有这里有精品99| 国产免费一区二区三区四区乱码| 欧美激情久久久久久爽电影| 亚洲国产欧美在线一区| av女优亚洲男人天堂| 亚洲丝袜综合中文字幕| 亚洲av成人精品一二三区| 麻豆成人午夜福利视频| 欧美人与善性xxx| 亚洲欧美日韩无卡精品| 国产黄色视频一区二区在线观看| 三级经典国产精品| a级毛片免费高清观看在线播放| 亚洲成人一二三区av| av福利片在线观看| 91久久精品国产一区二区三区| 性色av一级| 亚洲av二区三区四区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 少妇人妻精品综合一区二区| 三级经典国产精品| 国产男女内射视频| av线在线观看网站| 国产爱豆传媒在线观看| 国产黄片视频在线免费观看| av线在线观看网站| 亚洲av中文av极速乱| 老司机影院成人| 国产伦理片在线播放av一区| 少妇人妻 视频| 午夜爱爱视频在线播放| .国产精品久久| 九色成人免费人妻av| 国产精品久久久久久久久免| 午夜福利在线观看免费完整高清在| www.色视频.com| 男人和女人高潮做爰伦理| 国产伦精品一区二区三区四那| 亚洲av男天堂| 国产大屁股一区二区在线视频| 亚洲怡红院男人天堂| a级毛色黄片| 永久网站在线| av国产精品久久久久影院| 亚洲图色成人| 日产精品乱码卡一卡2卡三| 亚洲欧美一区二区三区国产| 一区二区av电影网| 久热这里只有精品99| 91aial.com中文字幕在线观看| 啦啦啦中文免费视频观看日本| tube8黄色片| 99re6热这里在线精品视频| 中文资源天堂在线| 亚洲精品日韩在线中文字幕| 18禁裸乳无遮挡动漫免费视频 | 卡戴珊不雅视频在线播放| 亚洲欧美清纯卡通| 欧美xxxx性猛交bbbb| 网址你懂的国产日韩在线| 最近最新中文字幕免费大全7| 成人黄色视频免费在线看| 在线观看一区二区三区| 亚洲av免费在线观看| 少妇熟女欧美另类| kizo精华| 亚洲美女视频黄频| 亚洲精品色激情综合| 成人鲁丝片一二三区免费| 国产熟女欧美一区二区| 国产精品秋霞免费鲁丝片| 亚洲欧美清纯卡通| 伦精品一区二区三区| 人人妻人人看人人澡| 精华霜和精华液先用哪个| 成人无遮挡网站| 欧美97在线视频| 亚洲精品456在线播放app| 成人无遮挡网站| 久久久久网色| 精品99又大又爽又粗少妇毛片| 日韩一区二区视频免费看| 秋霞伦理黄片| 天堂网av新在线| 中文字幕免费在线视频6| 男女那种视频在线观看| av黄色大香蕉| 伦理电影大哥的女人| 精品人妻偷拍中文字幕| 欧美精品国产亚洲| 国产亚洲av嫩草精品影院| 人妻系列 视频| 成人亚洲精品av一区二区| 中文字幕亚洲精品专区| 男女那种视频在线观看| 国产伦精品一区二区三区四那| 高清午夜精品一区二区三区| 亚洲综合色惰| 亚洲av电影在线观看一区二区三区 | 男插女下体视频免费在线播放| 在线亚洲精品国产二区图片欧美 | 成人高潮视频无遮挡免费网站| 国语对白做爰xxxⅹ性视频网站| 丰满少妇做爰视频| 亚洲综合精品二区| 日日啪夜夜爽| 99热这里只有是精品在线观看| 美女xxoo啪啪120秒动态图| 丝袜脚勾引网站| 色播亚洲综合网| 毛片女人毛片| 亚洲精品一区蜜桃| 插阴视频在线观看视频| 中文字幕亚洲精品专区| 午夜视频国产福利| 日韩一本色道免费dvd| 国产伦精品一区二区三区视频9| 久久综合国产亚洲精品| 亚洲国产色片| 如何舔出高潮| 蜜桃亚洲精品一区二区三区| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品电影小说 | 啦啦啦在线观看免费高清www| 免费观看性生交大片5| 国产91av在线免费观看| 免费看光身美女| 老司机影院毛片| 欧美日韩视频高清一区二区三区二| 亚洲人成网站在线播| 在线亚洲精品国产二区图片欧美 | 成人高潮视频无遮挡免费网站| 777米奇影视久久| 国产精品嫩草影院av在线观看| 日韩成人伦理影院| 99热国产这里只有精品6| 成人鲁丝片一二三区免费| 亚洲国产精品国产精品| 国产午夜精品久久久久久一区二区三区| 国产高清国产精品国产三级 | av女优亚洲男人天堂| 男人添女人高潮全过程视频| 九九久久精品国产亚洲av麻豆| 国产免费一级a男人的天堂| 国产视频首页在线观看| 夜夜爽夜夜爽视频| 美女国产视频在线观看| 国产精品久久久久久久久免| 午夜日本视频在线| 99久久精品国产国产毛片| 大香蕉久久网| 国国产精品蜜臀av免费| 久久久久久久大尺度免费视频| 成人亚洲精品一区在线观看 | 欧美日韩视频精品一区| 日本欧美国产在线视频| 欧美变态另类bdsm刘玥| 精品国产三级普通话版| 在线观看美女被高潮喷水网站| 人妻 亚洲 视频| 日韩大片免费观看网站| 久久精品国产亚洲av天美| 少妇人妻一区二区三区视频| 99久久精品国产国产毛片| 国产白丝娇喘喷水9色精品| 国产伦精品一区二区三区视频9| 99久久中文字幕三级久久日本| 丝瓜视频免费看黄片| 国产黄频视频在线观看| 亚洲成色77777| 一级a做视频免费观看| 日本黄色片子视频| 欧美人与善性xxx| eeuss影院久久| 五月天丁香电影| 国产午夜精品一二区理论片| 婷婷色综合大香蕉| 啦啦啦在线观看免费高清www| 大又大粗又爽又黄少妇毛片口| 国产免费一级a男人的天堂| 777米奇影视久久| 欧美bdsm另类| 国产伦精品一区二区三区视频9| 国产一级毛片在线| 国产午夜精品一二区理论片| 亚洲精品成人久久久久久| 又粗又硬又长又爽又黄的视频| 欧美精品国产亚洲| 国产高潮美女av| 搡女人真爽免费视频火全软件| 国产精品一区二区性色av| 两个人的视频大全免费| 激情五月婷婷亚洲| 日本欧美国产在线视频| 成年女人在线观看亚洲视频 | 特级一级黄色大片| 国产精品熟女久久久久浪| 人妻少妇偷人精品九色| 伦精品一区二区三区| 亚洲av男天堂| 亚洲,一卡二卡三卡| 国产日韩欧美亚洲二区| 大陆偷拍与自拍| 久久99热这里只有精品18| 3wmmmm亚洲av在线观看| 高清欧美精品videossex| 午夜免费男女啪啪视频观看| 日韩欧美精品v在线| 欧美另类一区| 亚洲人成网站在线播| 免费看av在线观看网站| 人妻系列 视频| 国产精品99久久久久久久久| 久久综合国产亚洲精品| 午夜福利视频1000在线观看| 亚洲人与动物交配视频| 亚洲电影在线观看av| 一级毛片黄色毛片免费观看视频| 久久久久久久精品精品| a级一级毛片免费在线观看| 亚洲精品亚洲一区二区| 乱系列少妇在线播放| 汤姆久久久久久久影院中文字幕| 亚洲欧美中文字幕日韩二区| av天堂中文字幕网| 一区二区三区乱码不卡18| 三级男女做爰猛烈吃奶摸视频| 大码成人一级视频| 亚洲无线观看免费| 亚洲欧洲日产国产| 国产探花在线观看一区二区| av网站免费在线观看视频| 好男人在线观看高清免费视频| 日韩一区二区视频免费看| 国产真实伦视频高清在线观看| 777米奇影视久久| 久久久久精品性色| 日韩人妻高清精品专区| 99热国产这里只有精品6| 亚洲人成网站在线播| 欧美潮喷喷水| 亚洲精品国产成人久久av| 免费观看的影片在线观看| 国产黄片美女视频| 国产亚洲最大av| 尾随美女入室| 日日撸夜夜添| 久久久久久久久久成人| 久久影院123| 九草在线视频观看| 国产一区亚洲一区在线观看| 国产淫语在线视频| 国产女主播在线喷水免费视频网站| 亚洲国产色片| 亚洲成人av在线免费| 美女cb高潮喷水在线观看| 久久久久久久久久久免费av| 丰满少妇做爰视频| 伦精品一区二区三区| 爱豆传媒免费全集在线观看| 国产在视频线精品| 亚洲国产高清在线一区二区三| 国产黄频视频在线观看| 女人被狂操c到高潮| 色网站视频免费| 日韩一区二区三区影片| 国产精品一区二区在线观看99| 又粗又硬又长又爽又黄的视频| a级一级毛片免费在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 中文字幕久久专区| 美女脱内裤让男人舔精品视频| 在线观看国产h片| 少妇熟女欧美另类| h日本视频在线播放| 有码 亚洲区| 亚洲精品中文字幕在线视频 | 久久久久精品性色| 久久精品国产亚洲av涩爱| 国产女主播在线喷水免费视频网站| www.色视频.com| 国产片特级美女逼逼视频| 欧美激情国产日韩精品一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 97超视频在线观看视频| 性色av一级| 纵有疾风起免费观看全集完整版| 欧美另类一区| 国产精品一二三区在线看| 欧美极品一区二区三区四区| 美女高潮的动态| 国产伦在线观看视频一区| 日韩不卡一区二区三区视频在线| 亚洲精品影视一区二区三区av| 黄色欧美视频在线观看| 亚洲天堂av无毛| 日本wwww免费看| 日本色播在线视频| 久久综合国产亚洲精品| 极品教师在线视频| 久久久久久久久久久免费av| 亚洲怡红院男人天堂| 免费黄色在线免费观看| 大话2 男鬼变身卡| 国产免费一区二区三区四区乱码| 午夜免费男女啪啪视频观看| 大话2 男鬼变身卡| 久久久成人免费电影| 精品久久久久久久久av| 中文字幕久久专区| 久久久精品94久久精品| kizo精华| 日韩成人av中文字幕在线观看| 美女被艹到高潮喷水动态| av免费在线看不卡| 男人狂女人下面高潮的视频| 三级国产精品欧美在线观看| 少妇丰满av|