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

    土壤重金屬的健康風(fēng)險(xiǎn)評價(jià)及其參數(shù)不確定性的量化研究

    2015-07-20 19:08:51李飛王曉鈺李雪
    關(guān)鍵詞:實(shí)例不確定性重金屬

    李飛 王曉鈺 李雪

    摘要:基于美國環(huán)保署現(xiàn)行健康風(fēng)險(xiǎn)評價(jià)體系的理論、方法和實(shí)例污染場址情景構(gòu)建了土壤重金屬污染的多途徑健康風(fēng)險(xiǎn)評價(jià)模型,而后利用MonteCarlo模擬法和敏感度分析技術(shù)初步量化研究了評價(jià)中參數(shù)不確定性的影響程度,并篩選得出相關(guān)敏感參數(shù).結(jié)果表明:確定性評價(jià)下實(shí)例區(qū)域土壤中Cd和Ni的非致癌風(fēng)險(xiǎn)(HI)均小于1(HICd>HINi),對目標(biāo)受體暫時(shí)不會造成非致癌健康風(fēng)險(xiǎn);而Cd的致癌風(fēng)險(xiǎn)(RCd)大于1.0×10-6,已可能對目標(biāo)受體造成致癌風(fēng)險(xiǎn);在參數(shù)MonteCarlo模擬的輔助下,參數(shù)不確定性對于非致癌風(fēng)險(xiǎn)值的影響波動在0.1~0.2以內(nèi),而其引起致癌風(fēng)險(xiǎn)值的波動在1個(gè)數(shù)量級以內(nèi),實(shí)例中RCd的模擬值區(qū)間橫跨1.0×10-6,可能誤導(dǎo)決策;敏感度分析結(jié)果表明應(yīng)重點(diǎn)對實(shí)例區(qū)域蔬菜和土壤中的Cd濃度及區(qū)域目標(biāo)受體體重這3個(gè)敏感參數(shù)進(jìn)行資料搜集,如此可以在有效提高風(fēng)險(xiǎn)評價(jià)可信度的同時(shí)降低不必要的預(yù)算開支.

    關(guān)鍵詞:土壤;重金屬;風(fēng)險(xiǎn)評價(jià);不確定性分析;MonteCarlo模擬

    中圖分類號:X820.2 文獻(xiàn)標(biāo)識碼:A

    Abstract:Based on USEPA framework of health risk assessment, the multi-pathway health risk assessment model was developed for heavy metals in the case soils. MonteCarlo simulation was adopted to quantitatively study the effect degree of parameter uncertainty on assessment result, and sensitive parameters were further identified by sensitivity analysis. The results without MonteCarlo simulation showed that there was no noncarcinogenic risk because of 1>HICd>HINi while it was under carcinogenic risk due to RCd>1.0×10-6. However, with the MonteCarlo simulation of selected variables, it indicated that the parameter uncertainty made the non-carcinogenic risk assessment results varying within 0.1~0.2, and the carcinogenic risk assessment results varying in the range within 1 order of magnitudes. Therefore, the parameter uncertainty made the value range of RCd span 1.0×10-6 and it was probable to mislead the corresponding decisionmaking. Sensitivity analysis indicated that the contents of Cd in case soils and vegetables, and the weight for the target receptor should be treated as sensitive variables, which were the main source of parameter uncertainty. To further improve the credibility of assessment with cost-benefit consideration, data of the sensitive variables should be preferentially gathered.

    Key words:soils;heavy metals;risk assessment;uncertainty analysis;MonteCarlo simulation

    重金屬是一類具有富集性,并很難在環(huán)境中降解的有毒污染物[1].由于土壤與水體、大氣等環(huán)境介質(zhì)有著密切關(guān)系[2],土壤重金屬污染已經(jīng)引起了公眾的廣泛關(guān)注,并且近年來多起公眾污染事件與土壤重金屬污染有關(guān)[3-4].土壤重金屬污染健康風(fēng)險(xiǎn)評價(jià)作為國際廣泛認(rèn)可的污染預(yù)防和修復(fù)決策的輔助工具,在當(dāng)下中國有著重要的現(xiàn)實(shí)意義.近年來,有一些學(xué)者基于健康風(fēng)險(xiǎn)評價(jià)的理論框架[5-7],分別對一些污染區(qū)域進(jìn)行了健康風(fēng)險(xiǎn)評價(jià)的研究,得出了一些具有指導(dǎo)性的結(jié)論,但不確定性卻在健康風(fēng)險(xiǎn)評價(jià)中貫穿始終,威脅著風(fēng)險(xiǎn)評價(jià)結(jié)論的可信度[8].鑒于模型不確定性和情景不確定性的固有性和隨機(jī)性特點(diǎn)[8],多數(shù)學(xué)者在進(jìn)行風(fēng)險(xiǎn)評價(jià)中重點(diǎn)考慮到了參數(shù)不確定性對評價(jià)結(jié)果的影響,并嘗試采用模糊數(shù)學(xué)、盲數(shù)理論、神經(jīng)網(wǎng)絡(luò)理論和隨機(jī)模擬理論等來降低評價(jià)過程中的參數(shù)不確定性,其中MonteCarlo隨機(jī)參數(shù)模擬法對于參數(shù)不確定性的良好控制具有廣泛的認(rèn)可度[9-11],并被美國環(huán)保署(US EPA)、中國環(huán)保部等權(quán)威機(jī)構(gòu)所推薦使用.但大多研究中僅定性或是半定量地描述了評價(jià)結(jié)果可能存在的參數(shù)不確定性大小及其影響,并提醒決策者加以注意[5,7,12],文獻(xiàn)中鮮有報(bào)道關(guān)于風(fēng)險(xiǎn)評價(jià)中參數(shù)不確定性影響程度的定量研究.

    鑒于參數(shù)不確定性的大小對決策的可信度有著不同程度的影響,本研究基于一個(gè)農(nóng)用土壤重金屬污染場址實(shí)例和US EPA現(xiàn)行的健康風(fēng)險(xiǎn)評價(jià)體系,根據(jù)實(shí)例情景構(gòu)建了土壤重金屬污染的多途徑健康風(fēng)險(xiǎn)評價(jià)模型;而后,借助參數(shù)MonteCarlo模擬下的風(fēng)險(xiǎn)評價(jià)結(jié)果和確定性風(fēng)險(xiǎn)評價(jià)結(jié)果的定量比對,初步量化研究了參數(shù)不確定性可能的影響程度,并通過敏感度分析進(jìn)一步開展了敏感參數(shù)排選,以期為不確定環(huán)境下我國土壤健康風(fēng)險(xiǎn)評價(jià)與相關(guān)環(huán)境管理決策的制定提供理論基礎(chǔ)和實(shí)踐經(jīng)驗(yàn).

    1實(shí)例區(qū)域概況與相關(guān)采樣分析

    1.1實(shí)例區(qū)域概況

    實(shí)例區(qū)域所在城市位于華北大平原的南端,地勢西高東低,南鄰黃河.該城市季風(fēng)特征明顯,冬季盛行東北風(fēng),夏季盛行西南風(fēng).在地理環(huán)境、大氣環(huán)流、地形,地勢等因子的綜合作用下,形成了暖溫大陸性季風(fēng)型氣候.該市域面積8 269平方公里、總?cè)丝诩s570萬,其中市區(qū)建成區(qū)面積140平方公里、人口120萬,并且該市是全國重要的商品糧基地和優(yōu)質(zhì)小麥生產(chǎn)基地,優(yōu)質(zhì)糧比重達(dá)85%.

    1.2實(shí)例區(qū)域的采樣分析

    實(shí)例區(qū)域農(nóng)用土壤的采樣監(jiān)測(土壤環(huán)境監(jiān)測技術(shù)規(guī)范HJ/166-2004)采用網(wǎng)格系統(tǒng)布點(diǎn)法,以100 m×100 m網(wǎng)格為樣品單元,設(shè)置20個(gè)采樣點(diǎn),每個(gè)樣品設(shè)3個(gè)平行樣,混勻后用四分法分別留取1 kg樣品,于室溫下自然風(fēng)干,細(xì)磨后過200目篩備用.采用四酸法進(jìn)行消解,用原子吸收分光光度法(日立Z5000原子吸收光譜測定儀)測定土壤中Cd和Ni的總量,結(jié)果見表1[13].本文數(shù)據(jù)均使用SPSS 16軟件進(jìn)行處理.將土壤重金屬含量實(shí)測數(shù)據(jù)進(jìn)行ShapiroWilk檢驗(yàn),Cd和Ni的sig.值均大于0.05,表明這2種重金屬含量的實(shí)測數(shù)據(jù)都呈正態(tài)分布.根據(jù)美國環(huán)保署的綜合風(fēng)險(xiǎn)信息系統(tǒng)(Integrated Risk Information System,IRIS)可知,研究中2種重金屬均具有非致癌健康風(fēng)險(xiǎn)和致癌健康風(fēng)險(xiǎn)[12-15].

    2.2人體暴露評估

    暴露評估的定義為“測量或估計(jì)人體暴露到目前存在于環(huán)境中物質(zhì)的程度、頻率和持續(xù)期間,或估計(jì)新化學(xué)物進(jìn)入環(huán)境中所可能引起的假設(shè)性暴露之過程”.沒有人體的暴露就沒有風(fēng)險(xiǎn),所以準(zhǔn)確的暴露途徑分析是構(gòu)建特征風(fēng)險(xiǎn)評價(jià)模型的關(guān)鍵前提,暴露途徑主要包括經(jīng)口食入,吸入(空氣),經(jīng)皮膚和黏膜吸收等.經(jīng)過對實(shí)例污染場址的調(diào)查分析,成年職業(yè)人群(當(dāng)?shù)剞r(nóng)民)其主要的暴露途徑為蔬菜經(jīng)口攝入暴露途徑、皮膚直接接觸暴露途徑、誤食土壤暴露途徑和土壤顆粒呼吸攝入暴露途徑.不同的暴露部位可能會有不同的吸收和代謝.全部吸收劑量(total absorbed dose)是每種暴露途徑所吸收劑量的總和.暴露評估需要本土的人體暴露參數(shù)及影響污染物在環(huán)境介質(zhì)中傳輸?shù)膮?shù),但由于我國尚未建立起各級別的暴露參數(shù)的數(shù)據(jù)庫,本文在盡量搜集我國現(xiàn)有相關(guān)文獻(xiàn)資料的基礎(chǔ)上,也參考US EPA的部分推薦參數(shù).

    參考USEPA的健康風(fēng)險(xiǎn)評估體系,下面是各個(gè)暴露途徑下的單位時(shí)間單位體重的暴露量(mg/(kg·d))估算式[15-19]:

    式中:wf為蔬菜的污染物濃度,mg/kg;IR為攝入率,kg/d;FI為被攝入污染源的比例,范圍0.0~1.0;w為土壤中化學(xué)物質(zhì)濃度,mg/kg;CF為轉(zhuǎn)換因子;SA為皮膚接觸面積,cm2/d;AF為皮膚(對土壤)黏附因子,mg/cm2;ABS為皮膚對化學(xué)物質(zhì)的吸收因子;IR′為攝取速率,mg/d;EF為暴露頻率,d·a-1;ED為暴露持續(xù)時(shí)間,a;BW為目標(biāo)受體體重,kg;AT為平均接觸時(shí)間,d;PEF為土壤塵擴(kuò)散因子,m3/kg;IRb為空氣的吸入量,m3/d.

    2.3風(fēng)險(xiǎn)特性描述

    定量的風(fēng)險(xiǎn)評價(jià)可更直觀、有效的表述風(fēng)險(xiǎn)的大小,同時(shí)也便于對污染因子的風(fēng)險(xiǎn)進(jìn)行篩選排序,為決策者提供科學(xué)的參考.本研究中2種重金屬都具有非致癌風(fēng)險(xiǎn)和致癌風(fēng)險(xiǎn)效應(yīng),故將風(fēng)險(xiǎn)定量表征公式列出如下.

    2.3.1非致癌風(fēng)險(xiǎn)的定量表征

    非致癌效應(yīng)假定具有閾值的機(jī)制,從而推導(dǎo)出參考劑量(reference dose, RfD).參考劑量的定義是“估計(jì)人類族群每天的暴露,此暴露在一生之中可能不會造成可察覺到有害健康效應(yīng)的風(fēng)險(xiǎn)”.非致癌風(fēng)險(xiǎn)的特征算式如下[15-19]:

    HI=HQf+HQa+HQs+HQb=

    (ADDfRfDf)+(ADDaRfDa)+(ADDsRfDs)+(ADDbRfDb).(5)

    式中:HI為某種污染物的非致癌污染指數(shù);HQf,HQa,HQs和HQb分別為蔬菜經(jīng)口攝入、皮膚直接接觸攝入、誤食土壤攝入和土壤顆粒呼吸攝入暴露途徑的非致癌風(fēng)險(xiǎn)商數(shù);ADD為某一非致癌物在某暴露途徑下目標(biāo)受體一生中平均每天每人每千克的暴露劑量,mg/(kg·d);RfD為某一非致癌物在某暴露途徑下的參考劑量,mg/(kg·d).

    關(guān)于非致癌風(fēng)險(xiǎn)的評判標(biāo)準(zhǔn),一般當(dāng)HI和HQ小于1時(shí),認(rèn)為風(fēng)險(xiǎn)較小或可以忽略;HI和HQ大于1時(shí),認(rèn)為存在風(fēng)險(xiǎn).

    2.3.2致癌風(fēng)險(xiǎn)的定量表征

    在人體或動物會造成癌癥的物質(zhì)被認(rèn)為具有非閾值效應(yīng),亦即并沒有安全暴露的水平.致癌風(fēng)險(xiǎn)的定量算法如式(6)[15-19]:

    R=∑kADIk×CSFk. (6)

    式中:ADIk為經(jīng)由暴露途徑k的每日平均暴露量,mg/(kg·d);CSFk為暴露途徑k的致癌斜率因子,(kg·d)/mg;其中ADIk與式(5)中的ADD在本研究中意義相同.

    關(guān)于致癌風(fēng)險(xiǎn)的評判標(biāo)準(zhǔn),目前采用的主要有美國環(huán)保局的健康風(fēng)險(xiǎn)評價(jià)標(biāo)準(zhǔn)、國際防輻射委員會推薦標(biāo)準(zhǔn),以及瑞典環(huán)保局、荷蘭建設(shè)和環(huán)境部推薦的評價(jià)標(biāo)準(zhǔn).其中美國環(huán)保局的致癌風(fēng)險(xiǎn)評價(jià)指南認(rèn)為風(fēng)險(xiǎn)水平處于10-4~10-6時(shí)的風(fēng)險(xiǎn)是可以接受的,其認(rèn)為的最大可接受風(fēng)險(xiǎn)為1.0×10-4;瑞典環(huán)保局、荷蘭建設(shè)和環(huán)境部推薦的健康危害風(fēng)險(xiǎn)度最大可接受限值為1.0×10-6.本研究選取最嚴(yán)格的1.0×10-6作為判別標(biāo)準(zhǔn).

    2.4相關(guān)確定性評價(jià)參數(shù)的選取

    表2列出了本研究暴露評估中選取的模型參數(shù).研究中2種重金屬不同暴露途徑的參考劑量(RfD)值(mg/(kg·d))、致癌重金屬的致癌強(qiáng)度系數(shù)值((kg·d)/mg)參考IRIS和風(fēng)險(xiǎn)評價(jià)信息系統(tǒng)(RAIS)的相關(guān)取值[15,19].

    2.5MonteCarlo模擬

    MonteCarlo模擬又稱隨機(jī)抽樣或統(tǒng)計(jì)試驗(yàn)方法,屬于計(jì)算數(shù)學(xué)的一個(gè)分支,是由Nicholas Metropolis在二次世界大戰(zhàn)期間提出的,而Von Neumann是Monte Carlo方法的正式奠基者,他與Stanislaw Ulam合作建立了概率密度函數(shù)、反累積分布函數(shù)的數(shù)學(xué)基礎(chǔ),以及偽隨機(jī)數(shù)產(chǎn)生器,現(xiàn)此方法在金融工程學(xué),宏觀經(jīng)濟(jì)學(xué),生物醫(yī)學(xué)和計(jì)算物理學(xué)(如粒子輸運(yùn)計(jì)算、量子熱力學(xué)計(jì)算、空氣動力學(xué)計(jì)算)等領(lǐng)域已得到應(yīng)用廣泛,效果良好[10-11].

    故本研究將MonteCarlo模擬引入構(gòu)建的土壤健康風(fēng)險(xiǎn)評價(jià)模型,并基于評價(jià)中參數(shù)的MonteCarlo模擬可以良好的控制參數(shù)不確定性的前提下,利用確定性風(fēng)險(xiǎn)評價(jià)結(jié)果與基于MonteCarlo模擬的風(fēng)險(xiǎn)評價(jià)結(jié)果的對比分析,量化研究實(shí)例風(fēng)險(xiǎn)評價(jià)中參數(shù)不確定性的影響大小,并進(jìn)一步進(jìn)行參數(shù)的敏感度分析,其主要模擬步驟為[10-11]:1)確定模型隨機(jī)變量,即確定影響評價(jià)結(jié)果的隨機(jī)因素,本研究中選取表3中的8個(gè)特征參數(shù)為隨機(jī)變量(包括CSCd,CSNi,BW,AT,IRb,IR′,CFfCd和CFfNi);2)構(gòu)建隨機(jī)因素的概率分布模型,主要通過實(shí)地采樣檢測和歷史經(jīng)驗(yàn)判定等方法,在本研究中采用歷史經(jīng)驗(yàn)和實(shí)地采樣檢測相結(jié)合的方法;3)將所得到的隨機(jī)數(shù)轉(zhuǎn)化為輸入?yún)?shù)的抽樣值,主要方法為MonteCarlo抽樣和拉丁超立方抽樣(Latin Hypercube Sampling,LHS),其中MonteCarlo抽樣一般從樣本分布較少的低概率區(qū)進(jìn)行抽樣,即為偏尾端抽樣;LHS抽樣則是由樣本整體分布考慮[10,12],這說明LHS方法更適合構(gòu)建小樣本的概率分布,故本文采用LHS法;4)整理分析所得模擬評價(jià)結(jié)果,其中主要包括評價(jià)指標(biāo)的期望值、概率分布和累積概率分布等;5)參數(shù)的敏感度分析,主要目的是篩選關(guān)鍵影響因子,并以此提高相關(guān)數(shù)據(jù)收集預(yù)算的效用率.

    3結(jié)果與討論

    3.1確定性參數(shù)下實(shí)例區(qū)域土壤污染健康風(fēng)險(xiǎn)評價(jià)

    根據(jù)前文中所建的特征土壤重金屬污染的多途徑健康風(fēng)險(xiǎn)評價(jià)模型,基于式(1)~(4),代入表2中的各參數(shù)值,本節(jié)參數(shù)均取確定性值或均值,計(jì)算出實(shí)例區(qū)域農(nóng)用土壤重金屬經(jīng)各暴露途徑可能引起成人受體的重金屬攝入量,計(jì)算結(jié)果見表3.

    由表3,對于Cd和Ni來說,其各途徑的日均暴露量的高低排序均為:蔬菜經(jīng)口暴露途徑>誤食土壤暴露途徑>皮膚接觸途徑的重金屬暴露劑量>土壤顆粒呼吸攝入暴露途徑.其中,Cd的各途徑的日均暴露量差異相對較大,而對于Ni來說各途徑的日均暴露量差異在3個(gè)數(shù)量級以內(nèi),并且其蔬菜經(jīng)口暴露途徑和誤食土壤暴露途徑的日均暴露量處于同一數(shù)量級.

    為進(jìn)一步研究該地區(qū)健康風(fēng)險(xiǎn)的現(xiàn)狀,利用式(5)和表2中的數(shù)據(jù),經(jīng)計(jì)算得出土壤重金屬在4種暴露途徑下的危害商數(shù)值(HQ),見表3.由表3,2種重金屬不同途徑的非致癌風(fēng)險(xiǎn)HI均<1.Cd各暴露途徑HQ的高低排序?yàn)椋篐Qf>HQt>HQa>HQb;Ni各暴露途徑HQ的高低排序則為:HQb≈HQf>HQt>HQa.根據(jù)式(5)計(jì)算得到HICd和HINi分別為0.138和0.058,故可知該污染場址中的Cd和Ni的非致癌風(fēng)險(xiǎn)暫時(shí)較小或可以忽略,相比之下,Cd的非致癌風(fēng)險(xiǎn)相對嚴(yán)重,并且高風(fēng)險(xiǎn)貢獻(xiàn)途徑均主要為蔬菜經(jīng)口暴露途徑和誤食土壤暴露途徑.

    根據(jù)式(6)計(jì)算得出RCd和RNi分別為5.26×10-5和9.93×10-7.參比本研究選定的致癌風(fēng)險(xiǎn)標(biāo)準(zhǔn)值1.0×10-6可知該場址土壤中的Cd已存在較高的致癌風(fēng)險(xiǎn),需要有關(guān)部門立即采取相應(yīng)的修復(fù)治理措施.相比之下,土壤中Ni的引起致癌風(fēng)險(xiǎn)暫時(shí)在可接受范圍內(nèi).

    3.2參數(shù)不確定性對評價(jià)結(jié)果的定量影響分析

    本研究嘗試?yán)肕onteCarlo隨機(jī)參數(shù)模擬下的風(fēng)險(xiǎn)評價(jià)結(jié)果和確定性風(fēng)險(xiǎn)評價(jià)結(jié)果的定量比對,初步量化研究參數(shù)不確定性對風(fēng)險(xiǎn)評價(jià)結(jié)果的影響程度.

    故在其他計(jì)算過程保持不變的前提下,將表3中的帶概率分布的參數(shù)利用水晶球(Crystal Ball) 2000軟件設(shè)置為對應(yīng)的概率分布模型,而后又基于非致癌風(fēng)險(xiǎn)(HQ和HI)和致癌風(fēng)險(xiǎn)(R)設(shè)置了定義預(yù)測單元.而后設(shè)置MonteCarlo模擬的最大實(shí)驗(yàn)量為1 000,置信區(qū)間為95%,抽樣方法為拉丁超立方(Latin hypercube sampling),其它參數(shù)取軟件的默認(rèn)值.運(yùn)行模擬得出對研究區(qū)域土壤中各種金屬的評價(jià)模擬預(yù)測圖,如圖2~5所示.圖中Probability代表概率可信度,F(xiàn)requency代表頻數(shù),并且圖4,圖5中橫坐標(biāo)值分別需要乘以10-5和10-6.由圖2,圖3可知,參比于確定性評價(jià)結(jié)果HICd(0.138)和HINi(0.058),在參數(shù)的MonteCarlo模擬下的HICd和HINi的值區(qū)間分別為[0.01, 0.38]和[0, 0.15].在參數(shù)不確定性的影響下,Cd和Ni的非致癌風(fēng)險(xiǎn)商數(shù)值在0.1~0.2波動,但在此實(shí)例下參數(shù)不確定性的影響程度暫時(shí)不會改變最終非致癌風(fēng)險(xiǎn)評價(jià)的結(jié)論.

    由圖4,5可知,Cd的致癌風(fēng)險(xiǎn)評價(jià)的區(qū)間值[0.6×10-6, 1.5×10-5]可能低于或高于所選定的風(fēng)險(xiǎn)標(biāo)準(zhǔn)值1.0×10-6,此時(shí)參數(shù)不確定性可能會誤導(dǎo)Cd的污染防控決策.而對于Ni來說,其致癌風(fēng)險(xiǎn)的RNi值區(qū)間小于標(biāo)準(zhǔn)值1.0×10-6,故參數(shù)不確定性暫不會影響Ni的污染防控決策.由圖2~5,研究中參數(shù)不確定性對于致癌風(fēng)險(xiǎn)值造成的波動均在1個(gè)數(shù)量級以內(nèi),所以如果致癌風(fēng)險(xiǎn)評價(jià)結(jié)果R的數(shù)量級比風(fēng)險(xiǎn)標(biāo)準(zhǔn)值的數(shù)量級高或低2個(gè)數(shù)量級以上,則參數(shù)不確定性將可能不會影響到最終的污染防控決策.但是如果致癌風(fēng)險(xiǎn)評價(jià)結(jié)果R的數(shù)量級與風(fēng)險(xiǎn)標(biāo)準(zhǔn)值的數(shù)量級的差異在1個(gè)數(shù)量級以內(nèi)的話,建議進(jìn)行進(jìn)一步的資料搜集以提高評價(jià)結(jié)果的可信度.

    3.3參數(shù)的敏感度分析

    敏感性分析是指從眾多不確定性因素中找出對投資項(xiàng)目經(jīng)濟(jì)效益指標(biāo)有重要影響的敏感性因素,并分析、測算其對項(xiàng)目經(jīng)濟(jì)效益指標(biāo)的影響程度和敏感性程度,進(jìn)而判斷項(xiàng)目承受風(fēng)險(xiǎn)能力的一種不確定性分析方法[25].參數(shù)的小幅度變化能導(dǎo)致經(jīng)濟(jì)效果指標(biāo)的較大變化,則稱此參數(shù)為敏感性因素,反之則稱其為非敏感性因素.

    由于在參數(shù)不確定性的影響下,本研究中總風(fēng)險(xiǎn)值(R)可能會誤導(dǎo)決策;并且實(shí)例情景場址下的資料調(diào)查與搜集花費(fèi)往往占整個(gè)風(fēng)險(xiǎn)評價(jià)項(xiàng)目的總預(yù)算的50%以上,故本研究在成本控制和評價(jià)可靠性的綜合考慮下,進(jìn)一步對評價(jià)過程進(jìn)行了參數(shù)敏感度分析(結(jié)果見表4),以期篩選出對于R來說的敏感變量,從而盡可能提高預(yù)算的有效利用率.

    由表4可知,在本研究中選取的8個(gè)隨機(jī)變量關(guān)于R的敏感度的高低排序?yàn)椋菏卟酥蠧d的濃度(78.7%)>目標(biāo)受體的體重(-13.3%)>土壤中Cd的濃度(7.5%)>土壤經(jīng)口攝入量(0.1%)≈平均接觸時(shí)間(-0.1%)≈蔬菜中的Ni的濃度(-0.1%)≈空氣吸入量(-0.1%)>蔬菜中Ni的濃度(0.0%).上述數(shù)據(jù)說明,對于評價(jià)結(jié)果R,蔬菜和土壤中Cd的濃度和目標(biāo)受體的體重這3個(gè)變量具有高敏感度,對評價(jià)結(jié)果起決定性作用;而相比之下,土壤經(jīng)口攝入量、平均接觸時(shí)間和空氣吸入量這些參數(shù)的感敏度相對較低,對R影響較小.故在R的值可能會誤導(dǎo)決策的前提下,需要進(jìn)一步通過對實(shí)例區(qū)域蔬菜和土壤中的Cd濃度和區(qū)域目標(biāo)受體體重這3個(gè)敏感參數(shù)的資料再搜集(包括歷史參數(shù)整理和實(shí)地檢測分析等)來降低評價(jià)中的參數(shù)不確定性,進(jìn)而有的放矢地提高評價(jià)結(jié)論的可信度和預(yù)算的有效使用率.

    4結(jié)論

    1)確定性健康風(fēng)險(xiǎn)評價(jià)結(jié)果顯示實(shí)例區(qū)域Cd和Ni的HI均小于1(HICd>HINi),對人體暫時(shí)不會造成非致癌健康風(fēng)險(xiǎn);而RCd>1.0×10-6,已經(jīng)對該區(qū)域造成致癌風(fēng)險(xiǎn),需要引起有關(guān)部門的注意.

    2)在參數(shù)MonteCarlo模擬的輔助下,定量研究表明參數(shù)不確定性對于總非致癌風(fēng)險(xiǎn)商數(shù)值的波動在0.1~0.2以內(nèi),而其對致癌風(fēng)險(xiǎn)值造成的波動在1個(gè)數(shù)量級以內(nèi),實(shí)例中RCd的模擬值區(qū)間橫跨1.0×10-6,可能誤導(dǎo)決策.

    3)敏感度分析的結(jié)果顯示對于總致癌風(fēng)險(xiǎn)值來說,實(shí)例區(qū)域蔬菜和土壤中的Cd濃度及區(qū)域目標(biāo)受體體重應(yīng)作為敏感參數(shù),重點(diǎn)搜集這3個(gè)參數(shù)的信息可有效地提高結(jié)論可信度和預(yù)算的效用率.

    4)本研究中方法在技術(shù)參數(shù)、特征模型架構(gòu)等方面仍需要進(jìn)一步完善,并且需要更系統(tǒng)的區(qū)域流行病學(xué)調(diào)查研究.

    參考文獻(xiàn)

    [1]LI F, HUANG J H, ZENG G M, et al. Spatial risk assessment and sources identification of heavy metals in surface sediments from the Dongting Lake, Middle China[J]. Journal of Geochemical Exploration, 2013, 132: 75-83.

    [2]陳懷滿. 環(huán)境土壤學(xué)[M].2版. 北京: 科學(xué)出版社, 2010: 2-235.

    CHEN Huaiman. Environmental soil science[M].2nd ed. Peking: Science Press, 2010: 2-235. (In Chinese)

    [3]王東宇, 張勇. 2006年中國城市飲用水源突發(fā)污染事件統(tǒng)計(jì)及分析[J]. 安全與環(huán)境學(xué)報(bào), 2007, 7(6): 150-155.

    WANG Dongyu, ZHANG Yong. Statistical analysis on drinking water source and supply system contamination threats and incidents for urban areas in China during 2006[J]. Journal of Safety and Environment, 2007, 7(6): 150-155. (In Chinese)

    [4]李秋虹, 方欣, 曹兆進(jìn). 突發(fā)環(huán)境污染事件應(yīng)急管理信息系統(tǒng)研究進(jìn)展[J]. 環(huán)境與健康雜志, 2008, 25(2):177-179.

    LI Qiuhong, FANG Xin, CAO Zhaojin. The development of managing information system in initial response to environmental incidents of public health[J]. Journal of Environment and Health, 2008, 25(2): 177-179. (In Chinese)

    [5]ZHENG N, LIU J S, WANG Q C, et al. Health risk assessment of heavy metal exposure to street dust in the zinc smelting district, Northeast of China[J]. Science of the Total Environment, 2010, 408(4): 726-733.

    [6]CHEN J Q, WANG Z X, WU X, et al. Source and hazard identification of heavy metals in soils of Changsha based on TIN model and direct exposure method[J]. Transactions of Nonferrous Metals Society of China, 2011, 21(3): 642-651.

    [7]?;勰?, 袁興中, 曾光明, 等. 基于動態(tài)聚類分析的水環(huán)境健康風(fēng)險(xiǎn)綜合評價(jià)[J]. 湖南大學(xué)學(xué)報(bào):自然科學(xué)版, 2010, 37(9): 73-78.

    ZHU Huina, YUAN Xingzhong, ZENG Guangming, et al. Integrated assessment of the environmental health risks of water sources based on dynamic cluster analysis[J]. Journal of Hunan University: Natural Science, 2010, 37(9): 73-78. (In Chinese)

    [8]LI F, HUANG J H, ZENG G M, et al. Multimedia health impact assessment: A study of the scenariouncertainty[J]. Journal of Central South University, 2012, 19(10): 2901-2909.

    [9]張應(yīng)華, 劉志全, 李廣賀, 等. 基于不確定性分析的健康環(huán)境風(fēng)險(xiǎn)評價(jià)[J]. 環(huán)境科學(xué), 2007, 28(7): 1409-1414.

    ZHANG Yinghua, LIU Zhiquan, LI Guanghe, et al. Uncertainty analysis of health risk assessment caused by benzene contamination in a contaminated site[J]. Environmental Science, 2007, 28(7): 1409-1414. (In Chinese)

    [10]EVANS J R, OLSON D L. Statistics, data analysis, and decision modeling[M]4th ed.New Jersey: Prentice Hall, 2002: 93-427.

    [11]李飛, 黃瑾輝, 曾光明, 等. 基于MonteCarlo模擬的土壤環(huán)境重金屬污染評價(jià)法與實(shí)例研究[J]. 湖南大學(xué)學(xué)報(bào):自然科學(xué)版, 2013, 40(9): 103-108.

    LI Fei, HUANG Jinhui, ZENG Guangming, et al. Probabilistic environmental quality assessment method and case study of soil heavy metals based on MonteCarlo simulation[J]. Journal of Hunan University: Natural Science, 2013, 40(9): 103-108. (In Chinese)

    [12]車飛, 于云江, 胡成, 等. 沈撫灌區(qū)土壤重金屬污染健康風(fēng)險(xiǎn)初步評價(jià)[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2009, 28(7): 1439-1443.

    CHE Fei, YU Yunjiang, HU Cheng, et al. Preliminary health assessment of heavy metals in soil in Shenfu irrigation area[J]. Journal of Agroenvironment Science, 2009, 28(7): 1439-1443. (In Chinese)

    [13]王曉鈺, 李飛. 農(nóng)用土壤重金屬多受體健康風(fēng)險(xiǎn)評價(jià)模型及實(shí)例應(yīng)用[J]. 環(huán)境工程, 2014, 32(1): 120-125.

    WANG Xiaoyu, LI Fei. Environmental health risk assessment model of soil heavy metals for multireceptor and its practical application[J]. Environmental Engineering, 2014, 32(1): 120-125.(In Chinese)

    [14]曾光明, 卓利, 鐘政林, 等. 水環(huán)境健康風(fēng)險(xiǎn)評價(jià)模型[J]. 水科學(xué)進(jìn)展, 1998, 9(3): 212-217.

    ZENG Guangming, ZHUO Li, ZHONG Zhenglin, et al. Assessment models for water environmental health risk analysis[J]. Advances in Water Science, 1998, 9(3): 212-217. (In Chinese)

    [15]CHEN H, LU X W, LI L Y. Spatial distribution and risk assessment of metals in dust based on samples from nursery and primary schools of Xian, China[J]. Atmospheric Environment, 2014, 88: 172-182.

    [16]陳鴻漢, 諶宏偉, 何江濤, 等. 污染場地健康風(fēng)險(xiǎn)評價(jià)的理論和方法[J]. 地學(xué)前緣, 2006, 13(1): 216-223.

    CHEN Honghan, CHEN Hongwei, HE Jiangtao, et al. Healthbased risk assessment of contaminated sites: principles and methods[J]. Earth Science Frontiers, 2006, 13(1): 216-223. (In Chinese)

    [17]US EPA. Risk assessment guidance forsuperfund. volume 1, human health evaluation manual. Part A[R]. Washington DC: Office of Emergency and Remedial Response, U.S. Environmental Protection Agency, 1990.

    [18]MCKONE T E. CalTOX: A multimedia totalexposure model for hazardous waste sites: Part Ⅰ~Ⅲ. Prepared for the state of California[R]. Livermore: Department of Toxic Substances Control, Lawrence Livermore National Laboratory, 1993.

    [19]于云江. 環(huán)境污染的健康風(fēng)險(xiǎn)評估與管理技術(shù)[M]. 北京: 中國環(huán)境科學(xué)出版社, 2011: 1-220.

    YU Yunjiang. Health risk assessment and management techniques of environmental pollution[M]. Peking: China Environmental Science Press, 2011: 1-220. (In Chinese)

    [20]段小麗. 暴露參數(shù)的研究方法及其在環(huán)境健康風(fēng)險(xiǎn)評價(jià)中的應(yīng)用[M]. 北京: 科學(xué)出版社, 2012: 1-179.

    DUAN Xiaoli. The research method of the exposure parameters and its application in environmental health risk assessment[M]. Peking: Science Press, 2012: 1-179. (In Chinese)

    [21]趙秀閣, 黃楠, 段小麗, 等. 環(huán)境健康風(fēng)險(xiǎn)評價(jià)中的皮膚暴露參數(shù)[J]. 環(huán)境與健康雜志, 2012, 29(2): 124-126.

    ZHAO Xiuge, HUANG Nan, DUAN Xiaoli, et al.Dermal exposure factors in environmental health risk assessment[J]. Journal of Environment and Health, 2012, 29(2): 124-126. (In Chinese)

    [22]王宗爽, 段小麗, 王貝貝, 等. 土壤/塵健康風(fēng)險(xiǎn)評價(jià)中的暴露參數(shù)[J]. 環(huán)境與健康雜志, 2012, 29(2): 114-117.

    WANG Zongshuang, DUAN Xiaoli, WANG Beibei, et al. Exposure factors in health risk assessment for soil and dust exposure[J]. Journal of Environment and Health, 2012, 29(2): 114-117. (In Chinese)

    [23]王學(xué)鋒, 馮穎俊, 林海, 等. 新鄉(xiāng)市部分市售蔬菜中重金屬污染狀況與質(zhì)量評價(jià)[J]. 河南師范大學(xué)學(xué)報(bào):自然科學(xué)版, 2006, 34(3): 120-123.

    WANG Xuefeng, FENG Yingjun, LIN Hai, et al. Evaluation and research of heavy metal on marketed vegetables in Xinxiang[J]. Journal of Henan Normal University: Natural Science, 2006, 34(3): 120-123. (In Chinese)

    [24]師榮光, 趙玉杰, 高懷友, 等. 天津市郊蔬菜重金屬污染評價(jià)與特征分析[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2005, 24(S1): 169-173.

    SHI Rongguang, ZHAO Yujie, GAO Huaiyou, et al. Metals influencing the level of heavy metal contamination of vegetable s in Tianjin city as derived by principal components analysis[J]. Journal of Agroenvironment Science, 2005, 24(S1): 169-173. (In Chinese)

    [25]劉慧, 韓兵, 戴峰. 決策結(jié)論的敏感度分析[J]. 中國管理科學(xué), 2008, 16(S1): 104-107.

    LIU Hui, HAN Bing, DAI Feng. Susceptibility analysis of decisionmaking result[J]. Chinese Journal of Management Science, 2008, 16(S1): 104-107. (In Chinese)

    猜你喜歡
    實(shí)例不確定性重金屬
    法律的兩種不確定性
    法律方法(2022年2期)2022-10-20 06:41:56
    重金屬對膨潤土膨脹性的影響
    英鎊或繼續(xù)面臨不確定性風(fēng)險(xiǎn)
    中國外匯(2019年7期)2019-07-13 05:45:04
    測定不同產(chǎn)地寬筋藤中5種重金屬
    中成藥(2018年8期)2018-08-29 01:28:16
    具有不可測動態(tài)不確定性非線性系統(tǒng)的控制
    ICP-AES、ICP-MS測定水中重金屬的對比研究
    再生水回灌中DOM對重金屬遷移與保留問題研究
    完形填空Ⅱ
    完形填空Ⅰ
    從翻譯的不確定性看譯者主體性
    成人手机av| 国产精品综合久久久久久久免费 | 一进一出好大好爽视频| 中文字幕制服av| 成人av一区二区三区在线看| 建设人人有责人人尽责人人享有的| 国产成+人综合+亚洲专区| 国产成人精品久久二区二区免费| 一本综合久久免费| 丝袜美足系列| 精品卡一卡二卡四卡免费| 在线观看免费高清a一片| 国产激情欧美一区二区| 亚洲精品av麻豆狂野| √禁漫天堂资源中文www| 成年女人毛片免费观看观看9 | 午夜免费鲁丝| svipshipincom国产片| 精品一区二区三区四区五区乱码| 中文字幕色久视频| 十分钟在线观看高清视频www| 色婷婷av一区二区三区视频| 亚洲一区高清亚洲精品| 久久精品国产99精品国产亚洲性色 | 亚洲av片天天在线观看| av国产精品久久久久影院| 啦啦啦视频在线资源免费观看| 亚洲精品av麻豆狂野| 美女福利国产在线| 天天躁日日躁夜夜躁夜夜| 久久人人爽av亚洲精品天堂| 99国产精品99久久久久| 亚洲午夜理论影院| 午夜91福利影院| 亚洲成人免费电影在线观看| 午夜亚洲福利在线播放| 老鸭窝网址在线观看| 极品教师在线免费播放| 狠狠狠狠99中文字幕| 欧美精品一区二区免费开放| 操出白浆在线播放| 免费不卡黄色视频| 最新在线观看一区二区三区| 在线播放国产精品三级| 国产在线精品亚洲第一网站| 久久久久久久久免费视频了| av国产精品久久久久影院| 看黄色毛片网站| 韩国av一区二区三区四区| 超碰成人久久| 久久久久国产精品人妻aⅴ院 | 国产亚洲欧美在线一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲熟妇中文字幕五十中出 | 午夜福利视频在线观看免费| 在线观看一区二区三区激情| 国产野战对白在线观看| 国产野战对白在线观看| 亚洲精品久久成人aⅴ小说| 在线免费观看的www视频| 久久天堂一区二区三区四区| 精品人妻1区二区| 一本大道久久a久久精品| 亚洲欧美一区二区三区久久| 又黄又爽又免费观看的视频| 国产99白浆流出| 亚洲国产欧美日韩在线播放| 久久久国产欧美日韩av| 国产高清激情床上av| 亚洲久久久国产精品| 老汉色av国产亚洲站长工具| 中文字幕人妻丝袜制服| cao死你这个sao货| 久久精品国产亚洲av高清一级| 色精品久久人妻99蜜桃| 一区二区三区精品91| 亚洲熟女精品中文字幕| 精品久久蜜臀av无| 亚洲综合色网址| 国内久久婷婷六月综合欲色啪| 99久久综合精品五月天人人| e午夜精品久久久久久久| 美女扒开内裤让男人捅视频| 最新美女视频免费是黄的| 国产高清国产精品国产三级| 精品久久久久久久久久免费视频 | 99香蕉大伊视频| 国产极品粉嫩免费观看在线| 精品一区二区三区四区五区乱码| 色婷婷av一区二区三区视频| 国产亚洲欧美精品永久| av网站在线播放免费| www.精华液| 欧美久久黑人一区二区| 一级毛片高清免费大全| 欧美乱码精品一区二区三区| 成年版毛片免费区| 女同久久另类99精品国产91| 久久久久久人人人人人| 999久久久国产精品视频| 亚洲国产毛片av蜜桃av| 国产精品香港三级国产av潘金莲| 在线免费观看的www视频| 日韩熟女老妇一区二区性免费视频| 国产精品av久久久久免费| 51午夜福利影视在线观看| 在线观看午夜福利视频| 亚洲色图 男人天堂 中文字幕| 久9热在线精品视频| 国产国语露脸激情在线看| 亚洲欧美激情综合另类| 国产单亲对白刺激| 制服人妻中文乱码| 看免费av毛片| 国产在视频线精品| 国产精品98久久久久久宅男小说| 天堂俺去俺来也www色官网| 少妇粗大呻吟视频| 嫩草影视91久久| 国产一区有黄有色的免费视频| avwww免费| 人人妻人人爽人人添夜夜欢视频| 久久香蕉国产精品| 欧美激情高清一区二区三区| 亚洲精华国产精华精| 免费观看a级毛片全部| 狂野欧美激情性xxxx| 最新美女视频免费是黄的| 久久精品国产综合久久久| 亚洲国产精品合色在线| 亚洲专区国产一区二区| 在线观看www视频免费| 日韩 欧美 亚洲 中文字幕| 久久天堂一区二区三区四区| 少妇裸体淫交视频免费看高清 | 人成视频在线观看免费观看| 国产av精品麻豆| 老鸭窝网址在线观看| 超色免费av| 日本撒尿小便嘘嘘汇集6| 国产精品永久免费网站| 亚洲精品av麻豆狂野| 国产精品乱码一区二三区的特点 | 99国产综合亚洲精品| 男女免费视频国产| 欧美成狂野欧美在线观看| 最新在线观看一区二区三区| 国产成人欧美在线观看 | 欧美 日韩 精品 国产| 十八禁高潮呻吟视频| 99久久人妻综合| 亚洲第一av免费看| 欧美午夜高清在线| 女性生殖器流出的白浆| 国产伦人伦偷精品视频| 久久九九热精品免费| 免费不卡黄色视频| 国内久久婷婷六月综合欲色啪| 黄色怎么调成土黄色| 国产精品乱码一区二三区的特点 | 视频区欧美日本亚洲| 91av网站免费观看| 高清毛片免费观看视频网站 | 欧美日本中文国产一区发布| 日韩精品免费视频一区二区三区| 午夜福利一区二区在线看| 一级a爱片免费观看的视频| 国产91精品成人一区二区三区| 欧美日韩乱码在线| 欧美日韩成人在线一区二区| 首页视频小说图片口味搜索| 黄色a级毛片大全视频| 手机成人av网站| 男女午夜视频在线观看| 国产91精品成人一区二区三区| 男女之事视频高清在线观看| а√天堂www在线а√下载 | 国产精品99久久99久久久不卡| 久久久久久免费高清国产稀缺| 女人高潮潮喷娇喘18禁视频| 国产视频一区二区在线看| 咕卡用的链子| 下体分泌物呈黄色| 三上悠亚av全集在线观看| 国产淫语在线视频| 亚洲,欧美精品.| 欧美激情高清一区二区三区| 午夜福利在线观看吧| 女警被强在线播放| 国产精品欧美亚洲77777| 中文字幕av电影在线播放| 99国产精品一区二区三区| 日本黄色视频三级网站网址 | 99精品久久久久人妻精品| 久久久久国产一级毛片高清牌| 国产成人啪精品午夜网站| 精品国内亚洲2022精品成人 | 欧美国产精品va在线观看不卡| 亚洲午夜精品一区,二区,三区| 热re99久久精品国产66热6| 久久久国产欧美日韩av| 在线视频色国产色| 国产欧美亚洲国产| 久久精品成人免费网站| av福利片在线| 在线播放国产精品三级| 国产片内射在线| 国产亚洲精品一区二区www | 午夜福利一区二区在线看| 90打野战视频偷拍视频| 亚洲aⅴ乱码一区二区在线播放 | 精品久久久久久久毛片微露脸| 久久久久久久精品吃奶| 一区二区三区精品91| 国产亚洲欧美98| 国产片内射在线| 成人永久免费在线观看视频| 欧美日韩av久久| 精品国产亚洲在线| 天天躁狠狠躁夜夜躁狠狠躁| 在线永久观看黄色视频| 国产成人精品在线电影| 亚洲一卡2卡3卡4卡5卡精品中文| 免费av中文字幕在线| 一边摸一边抽搐一进一出视频| 777米奇影视久久| 欧美人与性动交α欧美精品济南到| 婷婷丁香在线五月| 久久热在线av| 一边摸一边抽搐一进一出视频| 国产又爽黄色视频| 久久精品国产亚洲av香蕉五月 | 在线播放国产精品三级| 欧美 亚洲 国产 日韩一| 老熟妇仑乱视频hdxx| tube8黄色片| 成年人免费黄色播放视频| 性色av乱码一区二区三区2| 母亲3免费完整高清在线观看| 亚洲熟女精品中文字幕| 国产成人啪精品午夜网站| 国产亚洲精品一区二区www | 男女之事视频高清在线观看| av福利片在线| 午夜免费鲁丝| 一级黄色大片毛片| 日韩欧美三级三区| 国产成人影院久久av| 欧美国产精品一级二级三级| 国产在视频线精品| 国产男靠女视频免费网站| 欧美日本中文国产一区发布| 国产精品偷伦视频观看了| av一本久久久久| 麻豆成人av在线观看| 宅男免费午夜| 亚洲av欧美aⅴ国产| 黑人猛操日本美女一级片| 国产免费男女视频| 在线观看一区二区三区激情| 亚洲性夜色夜夜综合| 午夜精品久久久久久毛片777| av国产精品久久久久影院| 妹子高潮喷水视频| 另类亚洲欧美激情| 在线看a的网站| 99riav亚洲国产免费| 亚洲欧美激情综合另类| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美日韩高清在线视频| 亚洲欧美激情在线| 国产一区在线观看成人免费| 男女免费视频国产| 熟女少妇亚洲综合色aaa.| 午夜激情av网站| 久久久精品免费免费高清| 亚洲aⅴ乱码一区二区在线播放 | 五月开心婷婷网| 免费在线观看亚洲国产| 国产亚洲精品一区二区www | 欧美丝袜亚洲另类 | 久久这里只有精品19| 国产激情欧美一区二区| 乱人伦中国视频| 国产精品av久久久久免费| 久久精品国产99精品国产亚洲性色 | 老司机亚洲免费影院| 欧美在线黄色| 交换朋友夫妻互换小说| 成人18禁在线播放| 欧美成人午夜精品| 捣出白浆h1v1| 99久久国产精品久久久| 国产1区2区3区精品| 成人av一区二区三区在线看| 久久久水蜜桃国产精品网| 黄色片一级片一级黄色片| 黑人操中国人逼视频| 午夜免费成人在线视频| 亚洲av日韩精品久久久久久密| 操出白浆在线播放| 韩国av一区二区三区四区| 亚洲国产毛片av蜜桃av| 精品免费久久久久久久清纯 | 波多野结衣av一区二区av| 51午夜福利影视在线观看| 人妻久久中文字幕网| 亚洲精品成人av观看孕妇| 日韩欧美一区二区三区在线观看 | 成人影院久久| 妹子高潮喷水视频| 美女高潮到喷水免费观看| 精品无人区乱码1区二区| 咕卡用的链子| 国产一区在线观看成人免费| 亚洲五月天丁香| 国产aⅴ精品一区二区三区波| 999久久久精品免费观看国产| 国产不卡av网站在线观看| a在线观看视频网站| 一级毛片女人18水好多| 成人国产一区最新在线观看| 国产av精品麻豆| 可以免费在线观看a视频的电影网站| 老司机亚洲免费影院| 国产欧美日韩一区二区三| 青草久久国产| 国产精品电影一区二区三区 | x7x7x7水蜜桃| 免费日韩欧美在线观看| 人人妻人人爽人人添夜夜欢视频| 69精品国产乱码久久久| 精品国产国语对白av| videos熟女内射| 啪啪无遮挡十八禁网站| 91老司机精品| 看黄色毛片网站| 国产深夜福利视频在线观看| 丝袜美腿诱惑在线| 免费av中文字幕在线| 每晚都被弄得嗷嗷叫到高潮| 精品国产一区二区三区久久久樱花| 大码成人一级视频| 国产亚洲精品久久久久久毛片 | 我的亚洲天堂| 亚洲 国产 在线| 美女高潮到喷水免费观看| 亚洲av熟女| 在线永久观看黄色视频| 精品第一国产精品| 免费高清在线观看日韩| 动漫黄色视频在线观看| 午夜成年电影在线免费观看| 男男h啪啪无遮挡| 久久久久久人人人人人| 1024香蕉在线观看| 久久精品成人免费网站| 免费观看精品视频网站| 不卡一级毛片| 在线天堂中文资源库| 国产区一区二久久| 亚洲精品一卡2卡三卡4卡5卡| 捣出白浆h1v1| 久久久久视频综合| 亚洲国产精品合色在线| 18禁国产床啪视频网站| 三级毛片av免费| 老司机深夜福利视频在线观看| 国产不卡av网站在线观看| 中文字幕人妻熟女乱码| 热99国产精品久久久久久7| 国产91精品成人一区二区三区| 亚洲人成电影免费在线| 欧美 亚洲 国产 日韩一| 国产蜜桃级精品一区二区三区 | 国产视频一区二区在线看| 视频区欧美日本亚洲| 久久午夜综合久久蜜桃| 亚洲精品在线美女| 欧美激情极品国产一区二区三区| 国产成人av教育| 国产精品乱码一区二三区的特点 | 精品国产乱子伦一区二区三区| av线在线观看网站| 99国产综合亚洲精品| 欧美另类亚洲清纯唯美| 亚洲精品国产区一区二| 女人久久www免费人成看片| 狂野欧美激情性xxxx| 妹子高潮喷水视频| 国产精品一区二区在线观看99| 亚洲一区高清亚洲精品| 777久久人妻少妇嫩草av网站| 亚洲av美国av| 日本一区二区免费在线视频| av不卡在线播放| 久久香蕉国产精品| 国产精品免费大片| 欧美日韩黄片免| 日韩欧美免费精品| 可以免费在线观看a视频的电影网站| 亚洲欧美日韩另类电影网站| videosex国产| 伦理电影免费视频| 亚洲全国av大片| 亚洲欧美一区二区三区黑人| 黄色怎么调成土黄色| 日韩 欧美 亚洲 中文字幕| 欧美精品一区二区免费开放| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久人人做人人爽| 午夜精品在线福利| 老熟妇乱子伦视频在线观看| 19禁男女啪啪无遮挡网站| 久久久国产成人精品二区 | 老司机在亚洲福利影院| 国产成人精品在线电影| 在线播放国产精品三级| 999久久久国产精品视频| 99国产精品一区二区蜜桃av | 国产av一区二区精品久久| 别揉我奶头~嗯~啊~动态视频| 老司机亚洲免费影院| 国产欧美日韩精品亚洲av| 少妇的丰满在线观看| 婷婷成人精品国产| 亚洲欧美一区二区三区黑人| 亚洲色图综合在线观看| 亚洲国产看品久久| 老司机午夜十八禁免费视频| 757午夜福利合集在线观看| 99精品欧美一区二区三区四区| 成人国语在线视频| 国产深夜福利视频在线观看| 99国产极品粉嫩在线观看| 日韩人妻精品一区2区三区| 黄片播放在线免费| 一区福利在线观看| 久久久精品国产亚洲av高清涩受| 在线观看66精品国产| 欧美日韩亚洲综合一区二区三区_| 91大片在线观看| 久久国产精品男人的天堂亚洲| 精品无人区乱码1区二区| 亚洲精品中文字幕在线视频| 色精品久久人妻99蜜桃| 在线观看舔阴道视频| 人妻一区二区av| 黑人猛操日本美女一级片| 人人妻,人人澡人人爽秒播| tocl精华| 国产av一区二区精品久久| 亚洲人成伊人成综合网2020| 久久久精品区二区三区| 免费人成视频x8x8入口观看| 啦啦啦视频在线资源免费观看| 国产蜜桃级精品一区二区三区 | 久久人人97超碰香蕉20202| 黄片大片在线免费观看| 热99国产精品久久久久久7| 欧美日韩中文字幕国产精品一区二区三区 | 人人澡人人妻人| 一a级毛片在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲精品中文字幕在线视频| av超薄肉色丝袜交足视频| 欧美乱码精品一区二区三区| 捣出白浆h1v1| 国产人伦9x9x在线观看| 免费少妇av软件| 精品无人区乱码1区二区| 12—13女人毛片做爰片一| 啪啪无遮挡十八禁网站| 老司机深夜福利视频在线观看| 国产免费现黄频在线看| 久久人人爽av亚洲精品天堂| 另类亚洲欧美激情| 老司机福利观看| 新久久久久国产一级毛片| 国产精品美女特级片免费视频播放器 | 欧美乱色亚洲激情| 露出奶头的视频| 亚洲精品久久午夜乱码| 久久天堂一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看 | 纯流量卡能插随身wifi吗| 一区二区三区精品91| 国产精品99久久99久久久不卡| 久久99一区二区三区| 欧美人与性动交α欧美软件| 国产免费现黄频在线看| 看免费av毛片| 一本一本久久a久久精品综合妖精| 夫妻午夜视频| 欧美黑人精品巨大| 亚洲avbb在线观看| 精品亚洲成国产av| 人人妻人人澡人人看| 国产精品综合久久久久久久免费 | 搡老熟女国产l中国老女人| 女性生殖器流出的白浆| 怎么达到女性高潮| 亚洲视频免费观看视频| 亚洲精品乱久久久久久| av中文乱码字幕在线| 免费在线观看黄色视频的| 不卡一级毛片| 亚洲五月婷婷丁香| 美女国产高潮福利片在线看| 精品免费久久久久久久清纯 | 一个人免费在线观看的高清视频| 亚洲精品久久午夜乱码| 欧美黑人欧美精品刺激| 91老司机精品| 国产一卡二卡三卡精品| 亚洲国产中文字幕在线视频| 一a级毛片在线观看| 国产精品免费大片| 在线观看免费视频网站a站| 欧美成人免费av一区二区三区 | 黄频高清免费视频| 无限看片的www在线观看| 中文字幕另类日韩欧美亚洲嫩草| 精品卡一卡二卡四卡免费| 大香蕉久久网| 精品国产乱子伦一区二区三区| 亚洲精品一二三| 正在播放国产对白刺激| 在线永久观看黄色视频| 色尼玛亚洲综合影院| 国产高清激情床上av| 欧美中文综合在线视频| 欧美在线黄色| 人妻久久中文字幕网| 99re在线观看精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产在线观看jvid| 在线av久久热| 日韩成人在线观看一区二区三区| e午夜精品久久久久久久| 91字幕亚洲| 欧美日韩一级在线毛片| 欧美日韩亚洲高清精品| 欧美日韩视频精品一区| 亚洲三区欧美一区| 久久精品成人免费网站| 亚洲精品一卡2卡三卡4卡5卡| av网站免费在线观看视频| 激情在线观看视频在线高清 | 757午夜福利合集在线观看| 久久人妻熟女aⅴ| 在线十欧美十亚洲十日本专区| 久久中文字幕一级| 高清欧美精品videossex| 国产精品自产拍在线观看55亚洲 | 亚洲av电影在线进入| 成人精品一区二区免费| 露出奶头的视频| 亚洲欧美一区二区三区久久| 亚洲av电影在线进入| 亚洲av成人一区二区三| 一个人免费在线观看的高清视频| 亚洲午夜理论影院| 丝袜美腿诱惑在线| 不卡av一区二区三区| 在线看a的网站| 黄色成人免费大全| 高清毛片免费观看视频网站 | 看片在线看免费视频| 老司机影院毛片| 无限看片的www在线观看| 在线看a的网站| 精品福利观看| 丝袜人妻中文字幕| 亚洲精品国产区一区二| 99香蕉大伊视频| 丝瓜视频免费看黄片| av片东京热男人的天堂| 免费黄频网站在线观看国产| 中文字幕色久视频| 日韩精品免费视频一区二区三区| 搡老乐熟女国产| 一本一本久久a久久精品综合妖精| 怎么达到女性高潮| 久久久久久久午夜电影 | 黑人操中国人逼视频| 狠狠婷婷综合久久久久久88av| 岛国在线观看网站| 国产精品久久视频播放| 日韩中文字幕欧美一区二区| 欧美人与性动交α欧美精品济南到| 国产精品综合久久久久久久免费 | 成年人午夜在线观看视频| 亚洲av欧美aⅴ国产| 国产精品亚洲一级av第二区| 91国产中文字幕| 老司机靠b影院| 久久久精品国产亚洲av高清涩受| 久久久国产欧美日韩av| 午夜久久久在线观看| 精品国产美女av久久久久小说| 9色porny在线观看| 午夜福利免费观看在线| 久久精品亚洲熟妇少妇任你| 18在线观看网站| 每晚都被弄得嗷嗷叫到高潮| 日本wwww免费看| 欧美在线黄色| 国产高清国产精品国产三级| 中文字幕人妻丝袜一区二区| 啦啦啦视频在线资源免费观看| 久久人妻福利社区极品人妻图片| 久久久久国产精品人妻aⅴ院 | av不卡在线播放| 亚洲精华国产精华精|