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

    以設(shè)計(jì)精準(zhǔn)修復(fù)方案為目標(biāo)的土壤重金屬形態(tài)分布研究

    2018-02-06 02:47:58王佳佳劉永兵盧一富苗向前呂利光
    關(guān)鍵詞:弱酸全量重金屬

    王佳佳, 李 翔, 羅 楠, 何 躍, 劉永兵②, 盧一富, 苗向前, 呂利光

    (1.輕工業(yè)環(huán)境保護(hù)研究所工業(yè)場(chǎng)地污染與修復(fù)北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100089; 2.環(huán)境保護(hù)部南京環(huán)境科學(xué)研究所, 江蘇 南京 210042; 3.濟(jì)源市環(huán)境科學(xué)研究所河南省土壤重金屬污染監(jiān)測(cè)與修復(fù)重點(diǎn)實(shí)驗(yàn)室, 河南 濟(jì)源 459001)

    重金屬在土壤中不斷累積、難以降解,不僅會(huì)影響農(nóng)作物的產(chǎn)量、品質(zhì)和食品安全,而且會(huì)隨著食物鏈的積累進(jìn)入人體,引發(fā)人體健康風(fēng)險(xiǎn)[1],因此,重金屬污染農(nóng)田土壤修復(fù)迫在眉睫。原位鈍化/穩(wěn)定化是針對(duì)土壤重金屬污染且應(yīng)用較為普遍的一種化學(xué)修復(fù)技術(shù),它通過(guò)向重金屬污染農(nóng)田土壤加入穩(wěn)定劑使污染元素的化學(xué)賦存形態(tài)發(fā)生轉(zhuǎn)變,即通過(guò)多種物理-化學(xué)作用使土壤中重金屬固定在化學(xué)物質(zhì)的晶格中或使其化學(xué)形態(tài)向不活潑的方向轉(zhuǎn)變[2-3],阻止重金屬元素在環(huán)境中遷移、擴(kuò)散,從而降低毒害作用[4]。該技術(shù)僅改變了重金屬的存在形態(tài),而不能徹底去除土壤中重金屬污染物,土壤經(jīng)修復(fù)后仍存在一定的潛在生態(tài)危害,因此,土壤穩(wěn)定化修復(fù)的方案設(shè)計(jì)至關(guān)重要,既要防止因試劑過(guò)量而破壞土壤生態(tài)功能,又要避免因用量不足而導(dǎo)致修復(fù)效果不佳。

    精確的修復(fù)方案設(shè)計(jì)有賴于對(duì)待修復(fù)區(qū)污染特征與污染程度的準(zhǔn)確識(shí)別,因此,對(duì)研究區(qū)域重金屬元素分布特征與污染程度的準(zhǔn)確認(rèn)識(shí)是設(shè)計(jì)穩(wěn)定化修復(fù)方案的重要基礎(chǔ),能夠?yàn)楹笃谛迯?fù)方案的設(shè)計(jì)提供詳實(shí)的數(shù)據(jù)與理論依據(jù),是制定修復(fù)方案的必要前提。目前,國(guó)內(nèi)外對(duì)于重金屬污染評(píng)價(jià)方法主要包括單因子污染指數(shù)法[5-6]、內(nèi)梅羅綜合污染指數(shù)法[7-8]、地積累指數(shù)法[9-10]和Hakanson潛在生態(tài)風(fēng)險(xiǎn)指數(shù)法[11-12]等,這些方法分別從不同角度評(píng)價(jià)了重金屬的污染程度。然而,重金屬的生物毒性不僅與其總量有關(guān),更大程度上由其形態(tài)分布所決定[13],也就是說(shuō)重金屬元素的生物有效量越高,毒性就越大。土壤元素生物有效態(tài)通常指土壤中生物可吸收的元素形態(tài)[14]。目前,土壤生物有效量的測(cè)定方法中以選擇性提取最常見(jiàn),選擇性提取一般分為單步提取和連續(xù)提取。單步提取按試劑種類可分為去離子水、鹽溶液(如CaCl2、MgCl2、NH4OAc等)、稀酸(如硝酸、鹽酸、醋酸、檸檬酸等)、絡(luò)合劑(EDTA、DTPA等)等提取方法。連續(xù)提取方法包括Tessier法和BCR提取法,Tessier法將土壤重金屬元素分為可交換態(tài)、碳酸鹽結(jié)合態(tài)、鐵錳氧化物結(jié)合態(tài)、有機(jī)物結(jié)合態(tài)和殘?jiān)鼞B(tài)[15],BCR法將重金屬分為弱酸提取態(tài)、可還原態(tài)、可氧化態(tài)和殘?jiān)鼞B(tài)。近年來(lái),經(jīng)RAURET等[16]改進(jìn)的BCR法多用于測(cè)定重金屬的遷移性,能夠更具體地呈現(xiàn)重金屬對(duì)環(huán)境的影響,而且穩(wěn)定性好,精度高,可比性強(qiáng),得到了較廣泛的應(yīng)用[17]。

    以河南省某重金屬污染農(nóng)田土壤為研究對(duì)象,在測(cè)定土壤重金屬全量與其形態(tài)分布的基礎(chǔ)上,利用統(tǒng)計(jì)學(xué)方法分析不同重金屬元素的弱酸提取態(tài)含量隨全量的變化方式,并應(yīng)用Kriging插值法獲得重金屬全量與弱酸提取態(tài)含量的空間分布特征,分別基于重金屬全量與弱酸提取態(tài)含量評(píng)價(jià)研究區(qū)域重金屬污染程度與生態(tài)風(fēng)險(xiǎn)。最后,結(jié)合特征污染物的空間分布與評(píng)價(jià)結(jié)果對(duì)研究區(qū)域劃分小區(qū),以期為穩(wěn)定化修復(fù)方案的準(zhǔn)確設(shè)計(jì)提供理論依據(jù)。

    1 材料與方法

    1.1 樣點(diǎn)布設(shè)與樣品采集

    研究區(qū)域位于河南省某重金屬污染農(nóng)田修復(fù)試點(diǎn)區(qū),土壤為褐土,碳酸鈣含量較高,pH呈弱堿性,研究區(qū)域土壤理化性質(zhì)如下:pH值為8.3,陽(yáng)離子交換量為13.1 cmol·kg-1,w(有機(jī)質(zhì))為19.1 g·kg-1,w(全氮)為1.0 g·kg-1,w(全磷)為0.7 g·kg-1,w(全鉀)為10.0 g·kg-1,w(堿解氮)為81.7 mg·kg-1,w(速效磷)為24.8 mg·kg-1,w(速效鉀)為107.7 mg·kg-1。農(nóng)作物分兩季種植,每年6月至9月種植玉米,10月至次年5月種植冬小麥,屬華北平原的邊緣地帶,氣候冷熱分明,干旱或半干旱季節(jié)明顯,年降水量為640~680 mm,夏季炎熱,降雨集中,冬季寒冷,雨雪稀少。

    研究區(qū)域地形平整,形狀規(guī)則,面積約為1 hm2,周邊3~7 km范圍內(nèi)分布有3個(gè)鉛鋅冶煉廠,是重金屬污染防治重點(diǎn)區(qū)域,常年受到鉛鋅冶煉廠的大氣沉降污染。于2015年10月底按圖1所示布點(diǎn)方法對(duì)研究區(qū)域劃定32個(gè)采樣單元(采樣單元面積為18.5 m×16.25 m),每個(gè)采樣單元內(nèi)采集5個(gè)樣品混勻后作為1個(gè)樣品,共采集32個(gè)0~20 cm土壤樣品。

    圖1 研究區(qū)域位置與采樣布點(diǎn)Fig.1 Location of the study area and distribution of soil sampling sites

    1.2 樣品測(cè)試

    將新鮮的土壤樣品均勻鋪在干凈的聚乙烯塑料袋上,置于陰涼通風(fēng)處自然晾干,剔除其中的礫石、植物殘根等大塊異物,按四分對(duì)角線取樣法分取一半樣品并用瑪瑙研缽碾磨過(guò)0.15 mm孔徑篩,用于土壤重金屬元素全量及形態(tài)分析。

    土壤Cd、Pb、Cu、Zn、Ni和Cr含量采用普通酸分解法消解、電感耦合等離子體質(zhì)譜儀(ICP-MS)測(cè)定[18],As含量采用王水消解、原子熒光法測(cè)定[19],土壤重金屬形態(tài)分析采用改進(jìn)的BCR法,具體操作方法參見(jiàn)文獻(xiàn)[2],實(shí)驗(yàn)過(guò)程設(shè)置空白實(shí)驗(yàn),每隔5個(gè)樣品設(shè)置1個(gè)平行實(shí)驗(yàn)并保證相對(duì)標(biāo)準(zhǔn)偏差(RSD) < 10%。

    1.3 數(shù)據(jù)處理

    采用Excel 2016軟件進(jìn)行數(shù)據(jù)處理,采用Kriging插值、Surfer 9.0軟件制作重金屬污染分布圖,采用SPSS 19.0軟件進(jìn)行統(tǒng)計(jì)分析。

    1.4 評(píng)價(jià)方法

    1.4.1單因子污染指數(shù)與內(nèi)梅羅綜合污染指數(shù)方法

    單因子污染指數(shù)計(jì)算簡(jiǎn)便,反映了土壤中單一污染物的污染程度。內(nèi)梅羅綜合污染指數(shù)法兼顧單因子污染指數(shù)的平均值和最高值,能夠較全面地反映研究區(qū)域的污染程度,計(jì)算公式如下:

    (1)

    (2)

    式(1)~(2)中,pi為污染物i的單因子污染指數(shù);P為污染物的綜合污染指數(shù);ci為污染物i的實(shí)測(cè)含量,mg·kg-1;Si為HJ/T 332—2006《食用農(nóng)產(chǎn)品產(chǎn)地環(huán)境質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)》[20]中污染物i的限值,mg·kg-1,由于研究區(qū)域?yàn)闈?jì)源市主要的農(nóng)作物種植區(qū)域,且該標(biāo)準(zhǔn)的要求更嚴(yán)格,因此,以此作為重金屬污染程度評(píng)價(jià)的依據(jù)。pi,max和pi,ave分別為土壤中不同污染物i單因子污染指數(shù)的最大值和平均值,評(píng)價(jià)方法分級(jí)標(biāo)準(zhǔn)[21]見(jiàn)表1。

    表1單因子污染指數(shù)和內(nèi)梅羅污染指數(shù)分級(jí)標(biāo)準(zhǔn)[21]

    Table1StandardforgradingofSingleFactorIndexandNemeroPollutionIndex

    等級(jí)單因子污染指數(shù)(pi)內(nèi)梅羅綜合污染指數(shù)(P)污染等級(jí)污染水平1pi≤07P≤07安全清潔2073P>3重污染重污染

    1.4.2潛在生態(tài)危害指數(shù)評(píng)價(jià)法

    潛在生態(tài)危害指數(shù)(RI,IR)法是Hakanson根據(jù)重金屬的元素豐度與釋放性能,結(jié)合不同元素的生物毒性和環(huán)境響應(yīng)特點(diǎn),從沉積學(xué)角度提出的用于評(píng)價(jià)水生生態(tài)系統(tǒng)與沉積物重金屬污染程度的常用方法[22]。盡管該方法依據(jù)水生生態(tài)系統(tǒng)重金屬的沉積原理與遷移途徑提出,但目前也被用來(lái)評(píng)價(jià)土壤污染程度,且具有一定的參考價(jià)值[11,23],該方法的計(jì)算公式為

    (3)

    式(3)中,IR為整體潛在生態(tài)風(fēng)險(xiǎn)指數(shù);Ei,r為重金屬i的生態(tài)風(fēng)險(xiǎn)系數(shù);Ti,r為重金屬i的毒性響應(yīng)系數(shù);Ci,f為重金屬i的污染系數(shù);Ci和Ci,n分別為重金屬i的實(shí)測(cè)含量和所對(duì)應(yīng)的參比值,mg·kg-1,以河南省土壤元素背景值為參比值。各元素的毒性響應(yīng)系數(shù)分別為Zn,1;Cr,2;Cu,Ni,Pb,5;As,10;Cd,30[8,22]。根據(jù)Ei,r和IR值劃分不同的潛在生態(tài)危害水平,分級(jí)標(biāo)準(zhǔn)見(jiàn)表2。

    表2潛在生態(tài)危害指數(shù)分級(jí)標(biāo)準(zhǔn)[22]

    Table2StandardforgradingofPotentialEcologicalRiskIndex

    生態(tài)風(fēng)險(xiǎn)系數(shù)(Ei,r)潛在生態(tài)風(fēng)險(xiǎn)指數(shù)(IR)生態(tài)危害程度Ei,r<40IR<150輕微(Ⅰ)40≤Ei,r<80150≤IR<300中等(Ⅱ)80≤Ei,r<160300≤IR<600強(qiáng)(Ⅲ)160≤Ei,r<320很強(qiáng)(Ⅳ)Ei,r≥320IR≥600極強(qiáng)(Ⅴ)

    1.4.3風(fēng)險(xiǎn)評(píng)價(jià)編碼法(RAC)

    重金屬元素在土壤中以多種形態(tài)存在,其中,弱酸提取態(tài)鍵和作用較弱,易于遷移,更容易被動(dòng)植物所利用,因此,常被看作可交換態(tài)與碳酸鹽結(jié)合態(tài)之和以表征土壤重金屬的有效部分[13]。RAC基于土壤重金屬的不同存在形態(tài)對(duì)其不同的結(jié)合力而提出,通過(guò)計(jì)算重金屬有效部分含量占總量的比例高低來(lái)評(píng)價(jià)土壤重金屬的風(fēng)險(xiǎn),土壤重金屬有效性越高,對(duì)環(huán)境構(gòu)成的風(fēng)險(xiǎn)越大,反之亦然。

    RAC考慮了重金屬存在形態(tài)對(duì)其環(huán)境風(fēng)險(xiǎn)的影響,應(yīng)用范圍較廣[24-25],該方法的風(fēng)險(xiǎn)等級(jí)標(biāo)準(zhǔn)[24]:若弱酸提取態(tài)含量占總量的比例(%)<1,風(fēng)險(xiǎn)等級(jí)為無(wú)風(fēng)險(xiǎn);若為1~<10,屬低風(fēng)險(xiǎn);若為10~<30,屬中等風(fēng)險(xiǎn);若為30~<50,屬高風(fēng)險(xiǎn);若為≥50,屬極高風(fēng)險(xiǎn)。

    2 結(jié)果與討論

    2.1 土壤重金屬元素含量

    由表3可知,除Cr以外,研究區(qū)域土壤其他重金屬元素含量的變異系數(shù)均小于0.5,分布較均勻[26],這不僅與周邊的污染源類型有關(guān),當(dāng)?shù)貧夂驐l件對(duì)其也有一定影響。研究區(qū)域土壤As、Cu、Zn、Ni和Cr含量符合GB 15618—1995《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》二級(jí)標(biāo)準(zhǔn)[27]和HJ/T 332—2006,Cd、Pb含量超標(biāo)情況較嚴(yán)重。

    土壤Cd含量為0.40~2.32 mg·kg-1,96.88%的樣品超過(guò)GB 15618—1995二級(jí)標(biāo)準(zhǔn)和HJ/T 332—2006,超標(biāo)1.39~2.87倍;Pb含量為19.23~170.30 mg·kg-1,符合GB 15618—1995二級(jí)標(biāo)準(zhǔn),96.88%的樣品超過(guò)HJ/T 332—2006,超標(biāo)0.35~1.13倍。研究區(qū)域Cd、Pb含量不僅超標(biāo)嚴(yán)重,而且遠(yuǎn)超出河南省土壤元素背景值[28],這與周邊鉛鋅冶煉廠多年的廢氣排放有關(guān)。

    表3土壤重金屬含量及其統(tǒng)計(jì)分析

    Table3Heavymetalconcentrationinsoilandstatisticalanalysis

    重金屬w/(mg·kg-1)平均值最大值最小值變異系數(shù)w/(mg·kg-1)河南省土壤元素背景值[28]《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》二級(jí)標(biāo)準(zhǔn)[27]《食用農(nóng)產(chǎn)品產(chǎn)地環(huán)境質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)》[20]Cd 190 232 04001700650606Pb1449117030192301822335080As733814366010982525Cu322436702399009200100100Zn9140104205200013625300300Ni2880513520770182746060Cr19767935777064632250250

    n=32。

    2.2 重金屬形態(tài)分析

    重金屬元素在土壤中以不同形態(tài)存在,存在形態(tài)的不同對(duì)重金屬的遷移規(guī)律、變化形式和生物有效性都有一定影響[29]。原位穩(wěn)定化主要通過(guò)降低重金屬活性、減少植物對(duì)重金屬的吸收作用以降低重金屬對(duì)農(nóng)作物的毒害作用,因此,土壤重金屬的形態(tài)分布與活性含量研究對(duì)制定穩(wěn)定化方案意義較大。由于研究區(qū)土壤重金屬分布較均勻,故隨機(jī)選取12個(gè)樣品分析不同元素的形態(tài)分布,結(jié)果見(jiàn)圖2。

    BCR形態(tài)分析方法通過(guò)逐級(jí)提取將土壤重金屬分為弱酸提取態(tài)、可還原態(tài)、可氧化態(tài)和殘?jiān)鼞B(tài)4種形態(tài)。弱酸提取態(tài)一般指水或酸可溶態(tài)、可交換態(tài),活性最高,可還原態(tài)以鐵錳氧化物結(jié)合態(tài)為主,可氧化態(tài)以有機(jī)物結(jié)合態(tài)為主,殘?jiān)鼞B(tài)最穩(wěn)定,不能被植物吸收利用[14]。

    由圖2可知,研究區(qū)域土壤Cd超標(biāo)情況最嚴(yán)重,且弱酸提取態(tài)Cd含量所占比例較高(5.27%~28.55%),平均值為24.06%,此形態(tài)活性最高,容易被植物吸收利用,對(duì)植物有效態(tài)的貢獻(xiàn)率較大[30],表明研究區(qū)域Cd活性強(qiáng),危害較大。

    圖2 重金屬元素BCR形態(tài)分析Fig.2 Fractionation of heavy metals by BCR

    土壤Pb主要以可還原態(tài)存在(57.76%~83.19%),平均值為80.35%,其活性稍差于弱酸提取態(tài),還原條件下容易釋放出來(lái),對(duì)植物具有較強(qiáng)的潛在毒害作用[31]。Pb的弱酸提取態(tài)含量較低,占0.47%~3.17%,由此可見(jiàn),土壤Pb在還原條件下對(duì)植物的危害性較大。

    土壤中As、Cu、Zn、Ni、Cr含量符合相應(yīng)的國(guó)家標(biāo)準(zhǔn),且主要以殘?jiān)鼞B(tài)存在,各元素殘?jiān)鼞B(tài)含量分別為As,87.09%~94.61%;Cu,68.11%~76.06%;Zn,69.42%~82.98%;Ni,79.21%~ 84.35%;Cr,69.74% ~84.80%。由于殘?jiān)鼞B(tài)能夠在土壤中穩(wěn)定存在,表明這5種重金屬元素的環(huán)境危害性較小。

    研究區(qū)域土壤中所有元素的可氧化態(tài)含量均較低,也就是說(shuō)以有機(jī)物結(jié)合態(tài)存在的重金屬含量較低,Cr的可氧化態(tài)含量較高,為12.33%~24.67%,其余元素含量平均值小于10%,這與該區(qū)域土壤有機(jī)質(zhì)含量較低(19.1 g·kg-1)有關(guān)。

    2.3 土壤重金屬含量與其弱酸提取態(tài)含量的回歸分析

    通過(guò)對(duì)12個(gè)點(diǎn)位重金屬全量和弱酸提取態(tài)含量作回歸分析以預(yù)測(cè)土壤中弱酸提取態(tài)重金屬含量隨其全量的變化方式,結(jié)果見(jiàn)表4。

    表4重金屬弱酸提取態(tài)含量與全量回歸分析預(yù)測(cè)模型

    Table4Regressionanalysisbasedpredictionmodelforweakacidextractableandtotalofheavymetals

    弱酸提取態(tài)含量/(mg·kg-1)全量/(mg·kg-1)回歸方程相關(guān)系數(shù)rYCdXCdYCd=0317XCd-01110982YPbXPbYPb=0032XPb-06910926YAsXAsYAs=0052XAs-01740760YCuXCuYCu=0030XCu-01400741YZnXZnYZn=0099XZn-38350933YNiXNiYNi=0009XNi+01670626YCrXCrYCr=0003XCr+00030752

    顯著性水平為α=0.05,0.05水平下相關(guān)系數(shù)的臨界值為0.576。

    7種重金屬元素的相關(guān)系數(shù)均大于0.05水平下的臨界值,表明弱酸提取態(tài)含量與重金屬全量之間存在顯著正相關(guān)關(guān)系。Cd、Pb和Zn的相關(guān)系數(shù)均大于0.9,表明重金屬全量的變化能夠很好地解釋弱酸提取態(tài)含量的變化;其次為As、Cu、Cr 3種元素,相關(guān)系數(shù)為0.741~0.760;Ni的相關(guān)系數(shù)稍小,為0.626,表明土壤弱酸提取態(tài)Ni含量變化的62.6%可由土壤Ni全量解釋。利用表4所示回歸方程可以預(yù)測(cè)研究區(qū)域土壤弱酸提取態(tài)含量隨重金屬全量的變化方式。

    2.4 重金屬全量與弱酸提取態(tài)含量的空間分布

    土壤重金屬的弱酸提取態(tài)活性最強(qiáng),最容易被植物吸收利用,結(jié)合重金屬元素的超標(biāo)情況對(duì)比分析研究區(qū)域Cd、Pb全量與弱酸提取態(tài)含量的空間分布。采用Kriging插值法對(duì)研究區(qū)域土壤12個(gè)點(diǎn)位Cd、Pb全量和弱酸提取態(tài)含量數(shù)據(jù)進(jìn)行插值得到空間分布圖,該方法是土壤科學(xué)研究中常用的插值方法,插值效果較好[32-33]。

    全量Cd與弱酸提取態(tài)Cd含量的空間分布見(jiàn)圖3,根據(jù)全量Cd的空間分布情況可將研究區(qū)域分為2大部分,北側(cè)Cd全量為2~2.5 mg·kg-1,南側(cè)為1.5~2 mg·kg-1。圖3還表明根據(jù)弱酸提取態(tài)Cd含量的空間分布也可將研究區(qū)域大致分為2個(gè)部分,含量為0.2~0.5 mg·kg-1的區(qū)域占主要部分,位于南側(cè),其次為0.5~1 mg·kg-1的區(qū)域,位于北側(cè)。

    全量Pb和弱酸提取態(tài)Pb含量的空間分布見(jiàn)圖3,兩者的空間分布差異較大。全量Pb整體分布較均勻,無(wú)需根據(jù)Pb含量的空間分布開展分區(qū)修復(fù)。而弱酸提取態(tài)Pb含量的空間異質(zhì)性則較大,含量為2~4 mg·kg-1的區(qū)域?yàn)橹饕糠?位于研究區(qū)域南側(cè),其次為4~6 mg·kg-1區(qū)域,位置相對(duì)靠北。

    根據(jù)重金屬全量與弱酸提取態(tài)含量的空間分布有必要對(duì)研究區(qū)域進(jìn)行分區(qū)修復(fù),如圖3實(shí)線所示。基于全量Cd和弱酸提取態(tài)Cd含量的空間分布對(duì)研究區(qū)域的分區(qū)結(jié)果存在一定差異,而基于弱酸提取態(tài)Pb和弱酸提取態(tài)Cd含量空間分布差異所得的分區(qū)結(jié)果則大致相同。由于活性重金屬含量對(duì)土壤和農(nóng)作物的影響較大,故后期開展原位穩(wěn)定化修復(fù)時(shí),可參照土壤弱酸提取態(tài)Cd含量的分區(qū)結(jié)果設(shè)計(jì)修復(fù)方案,根據(jù)不同分區(qū)的具體情況設(shè)計(jì)相應(yīng)的穩(wěn)定化方案,以實(shí)現(xiàn)精準(zhǔn)修復(fù)。

    Ⅰ—Cd-全量; Ⅱ—Cd-弱酸提取態(tài); Ⅲ—Pb-全量; Ⅳ—Pb-弱酸提取態(tài)。 重金屬含量單位為mg·kg-1。

    2.5 重金屬污染評(píng)價(jià)與風(fēng)險(xiǎn)評(píng)估

    表5顯示,土壤Cd含量超標(biāo)情況最嚴(yán)重,單因子污染指數(shù)最大,為重污染;Pb含量超標(biāo)情況次之,屬輕污染;其余重金屬元素符合HJ/T 332—2006,未受到污染。由此可見(jiàn),研究區(qū)域的主要污染物為Cd,其次為Pb。

    內(nèi)梅羅綜合污染指數(shù)評(píng)價(jià)結(jié)果表明研究區(qū)域整體受到中度污染,綜合污染指數(shù)為2.33。

    表5土壤重金屬污染的單因子污染指數(shù)評(píng)價(jià)結(jié)果

    Table5SingleFactorIndexevaluationofheavymetalspollutionofthesoil

    重金屬元素平均值/(mg·kg-1)單因子污染指數(shù)污染程度Cd190316重污染Pb14491181輕污染As733029未受污染Cu3224032未受污染Zn9140031未受污染Ni2880048未受污染Cr1976008未受污染

    單因子污染指數(shù)與內(nèi)梅羅綜合污染指數(shù)按照相應(yīng)的計(jì)算方法得到了研究區(qū)域的污染程度,但并未考慮土壤重金屬的累積效應(yīng)和不同元素的生理生態(tài)毒性,因此,兩者的評(píng)價(jià)結(jié)果僅能表示研究區(qū)域當(dāng)前的污染程度,而未能明確所存在的生態(tài)風(fēng)險(xiǎn)。Hakanson潛在生態(tài)風(fēng)險(xiǎn)指數(shù)法綜合考慮了重金屬元素的毒性效應(yīng)、環(huán)境響應(yīng)和沉積特點(diǎn),為研究區(qū)域土壤重金屬的潛在生態(tài)風(fēng)險(xiǎn)提供參考,評(píng)價(jià)結(jié)果見(jiàn)表6。

    表6土壤重金屬生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果

    Table6Ecologicalriskassessmentofheavymetalsinthesoil

    重金屬元素平均值/(mg·kg-1)單元素生態(tài)危害風(fēng)險(xiǎn)污染程度Cd19087734極強(qiáng)生態(tài)危害Pb144913266輕微As733748輕微Cu3224806輕微Zn9140146輕微Ni2880525輕微Cr1976063輕微

    研究區(qū)域土壤Cd含量不僅遠(yuǎn)遠(yuǎn)高于河南省土壤元素背景值,且其毒性響應(yīng)系數(shù)最大,因此,Cd的生態(tài)風(fēng)險(xiǎn)也較高,生態(tài)危害指數(shù)為877.34,存在極強(qiáng)生態(tài)危害。Pb含量與河南省土壤元素背景值相差較大,由于它毒性效應(yīng)較低,因而存在輕微生態(tài)風(fēng)險(xiǎn)。雖然As的毒性效應(yīng)較強(qiáng),研究區(qū)域土壤As含量低于河南省土壤元素背景值,因此,土壤As對(duì)環(huán)境的潛在生態(tài)危害也較輕。其余重金屬元素不僅含量較低,毒性也較小,因此,所產(chǎn)生的生態(tài)風(fēng)險(xiǎn)均較低。

    研究區(qū)域Cd的生態(tài)危害極強(qiáng),導(dǎo)致土壤整體也存在極強(qiáng)生態(tài)危害,整體風(fēng)險(xiǎn)為932.88。

    單因子污染指數(shù)和內(nèi)梅羅綜合污染指數(shù)以及Hakanson風(fēng)險(xiǎn)評(píng)價(jià)方法基于土壤重金屬含量,從不同角度詮釋了研究區(qū)域土壤重金屬的污染程度與存在風(fēng)險(xiǎn),然而上述評(píng)價(jià)方法均未考慮土壤重金屬存在形態(tài)對(duì)其風(fēng)險(xiǎn)的影響,風(fēng)險(xiǎn)評(píng)價(jià)編碼法基于重金屬弱酸提取態(tài)含量占總量的比例來(lái)評(píng)價(jià)重金屬風(fēng)險(xiǎn)的高低,充分考慮了活性重金屬對(duì)重金屬風(fēng)險(xiǎn)的貢獻(xiàn)。

    表7顯示,弱酸提取態(tài)Cd含量占總量比例較高,土壤Cd風(fēng)險(xiǎn)較高,為中等風(fēng)險(xiǎn);Cr的弱酸提取態(tài)含量占比小于1%,無(wú)風(fēng)險(xiǎn);其余5種元素的弱酸提取態(tài)含量占比介于1%~10%之間,對(duì)環(huán)境具有低風(fēng)險(xiǎn)。

    表7采用風(fēng)險(xiǎn)評(píng)價(jià)編碼法對(duì)土壤重金屬的評(píng)價(jià)結(jié)果

    Table7RiskassessmentofheavymetalsinthesoilwithRAC

    重金屬元素弱酸提取態(tài)含量占總量比例/%風(fēng)險(xiǎn)等級(jí)Cd2406中等Pb249低As266低Cu253低Zn531低Ni151低Cr027無(wú)

    Hakanson潛在生態(tài)危害指數(shù)法和風(fēng)險(xiǎn)評(píng)價(jià)編碼法分別基于土壤重金屬全量和有效態(tài)含量對(duì)研究區(qū)域的生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià),兩者的評(píng)價(jià)角度不同,結(jié)果也有一定差異。前者顯示研究區(qū)域Cd存在極強(qiáng)生態(tài)風(fēng)險(xiǎn),而后者則表明Cd為中等風(fēng)險(xiǎn),2種評(píng)價(jià)結(jié)果間的差異提示土壤重金屬修復(fù)需綜合考慮不同評(píng)價(jià)方法的評(píng)價(jià)結(jié)果。

    研究區(qū)域重金屬污染農(nóng)田土壤擬采用原位穩(wěn)定化技術(shù)進(jìn)行修復(fù),旨在減少土壤重金屬有效態(tài)含量,進(jìn)而抑制重金屬向農(nóng)作物的遷移作用,因此,以風(fēng)險(xiǎn)評(píng)價(jià)編碼法的評(píng)價(jià)結(jié)果為主。此外,由于Hakanson潛在生態(tài)危害指數(shù)主要依據(jù)水生系統(tǒng)與湖泊中污染物質(zhì)的沉積特性、暴露與遷移途徑而提出的用于沉積物污染程度的評(píng)價(jià),這與土壤中重金屬污染的累積效應(yīng)、暴露途徑和遷移途徑存在一定差異,因此,采用潛在生態(tài)危害指數(shù)法對(duì)研究區(qū)域進(jìn)行評(píng)價(jià)得到的結(jié)果僅可作為參考。

    3 結(jié)論

    (1)研究區(qū)域土壤Cd、Pb、As、Cu、Zn和Ni含量分布均勻,Cr含量變異系數(shù)稍大。土壤Cd含量平均值為1.90 mg·kg-1,超出GB 15618—1995二級(jí)標(biāo)準(zhǔn)和HJ/T 332—2006,弱酸提取態(tài)Cd含量占總量比例較高,平均值為24.06%;Pb含量平均值為144.91 mg·kg-1,符合GB 15618—1995二級(jí)標(biāo)準(zhǔn),超出HJ/T 332—2006,可還原態(tài)Pb含量占比較高,平均值為80.35%;土壤As、Cu、Zn、Ni和Cr主要以殘?jiān)鼞B(tài)存在,含量符合GB 15618—1995二級(jí)標(biāo)準(zhǔn)和HJ/T 332—2006。

    (2)回歸分析結(jié)果表明研究區(qū)域土壤重金屬含量與弱酸提取態(tài)含量存在顯著正相關(guān)性,重金屬弱酸提取態(tài)含量隨該元素全量的變化呈線性變化。

    (3)以土壤全量Cd與弱酸提取態(tài)Cd含量的空間分布為依據(jù)的分區(qū)結(jié)果存在一定差異,而根據(jù)弱酸提取態(tài)Cd和弱酸提取態(tài)Pb含量的分區(qū)結(jié)果大致相同,可參考弱酸提取態(tài)Cd含量的分區(qū)結(jié)果制定穩(wěn)定化修復(fù)方案。

    (4)單因子污染指數(shù)與內(nèi)梅羅綜合污染指數(shù)評(píng)價(jià)結(jié)果表明研究區(qū)域土壤整體受到中度污染,其中,Cd為重污染,Pb為輕污染,其余元素為未受污染;潛在生態(tài)危害指數(shù)評(píng)價(jià)結(jié)果顯示研究區(qū)域土壤整體存在極強(qiáng)生態(tài)危害,其中,Cd對(duì)環(huán)境的生態(tài)危害極強(qiáng),除Cd以外的其他元素為輕微生態(tài)危害;風(fēng)險(xiǎn)評(píng)價(jià)編碼法評(píng)價(jià)結(jié)果顯示土壤Cd為中等風(fēng)險(xiǎn),Cr為無(wú)風(fēng)險(xiǎn),其余元素為低風(fēng)險(xiǎn)。不同評(píng)價(jià)方法所導(dǎo)致的結(jié)果差異表明土壤重金屬修復(fù)需綜合考慮不同評(píng)價(jià)方法的評(píng)價(jià)結(jié)果。

    [1] 劉永兵,賈斌,李翔,等.海南省南渡江新坡河塘底泥養(yǎng)分狀況及重金屬污染評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(3):213-224.[LIU Yong-bing,JIA Bin,LI Xiang,etal.Characteristic of Nutrients and Evaluation of Heavy Metal Contamination on Sediments Among Xinpo Pond,Nandu River in Hainan Province[J].Transactions of the Chinese Society of Agricultural Engineering,2013,29(3):213-224.]

    [2] 李翔,劉永兵,宋云,等.石灰干化污泥對(duì)土壤重金屬穩(wěn)定化處理的效果[J].環(huán)境工程學(xué)報(bào),2014,8(8):3461-3470.[LI Xiang,LIU Yong-bing,SONG Yun,etal.Effect of Lime-Stabilized Sludge on Stabilization of Heavy Metals in Soil[J].Chinese Journal of Environmental Engineering,2014,8(8):3461-3470.]

    [3] 李翔,劉永兵,程言君,等.穩(wěn)定化處理對(duì)底泥利用后土壤重金屬形態(tài)及蔬菜重金屬含量的影響[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2016,35(7):1278-1285.[LI Xiang,LIU Yong-bing,CHENG Yan-jun,etal.Effects of Stabilization on Soil Heavy Metal Fractions and Vegetable Heavy Metal Content Under Amendment With Dredged River Sediments[J].Journal of Agro-Environment Science,2016,35(7):1278-1285.]

    [4] 郝漢舟,陳同斌,靳孟貴,等.重金屬污染土壤穩(wěn)定/固化修復(fù)技術(shù)研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報(bào),2011,22(3):816-824.[HAO Han-zhou,CHEN Tong-bin,JIN Meng-gui,etal.Recent Advance in Solidification/ Stabilization Technology for the Remediation of Heavy Metals-Contaminated Soil[J].Chinese Journal of Applied Ecology,2011,22(3):816-824.]

    [5] 韓君,徐應(yīng)明,溫兆飛,等.重慶某廢棄電鍍工業(yè)園農(nóng)田土壤重金屬污染調(diào)查與生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J].環(huán)境化學(xué),2014,33(3):432-439.[HAN Jun,XU Ying-ming,WEN Zhao-fei,etal.Investigation and Ecological Risk Assessment of Heavy Metal Pollution in Agriculture Soil of an Abandoned Electroplating Industrial Park of Chongqing,China[J].Environmental Chemistry,2014,33(3):432-439.]

    [6] 柳云龍,章立佳,韓曉非,等.上海城市樣帶土壤重金屬空間變異特征及污染評(píng)價(jià)[J].環(huán)境科學(xué),2012,33(2):599-605.[LIU Yun-long,ZHANG Li-jia,HAN Xiao-fei,etal.Spatial Variability and Evaluation of Soil Heavy Metal Contamination in the Urban-Transect of Shanghai[J].Environmental Science,2012,33(2):599-605.]

    [7] 吳瓊,趙同科,鄒國(guó)元,等.北京東南郊農(nóng)田土壤重金屬含量與環(huán)境質(zhì)量評(píng)價(jià)[J].中國(guó)土壤與肥料,2016(1):7-12.[WU Qiong,ZHAO Tong-ke,ZOU Guo-yuan,etal.Heavy Metals Contents and Environment Quality Evaluation of Farmland Soil in East-South Area of Beijing[J].Soil and Fertilizer Sciences in China,2016(1):7-12.]

    [8] 劉亞納,朱書法,魏學(xué)鋒,等.河南洛陽(yáng)市不同功能區(qū)土壤重金屬污染特征及評(píng)價(jià)[J].環(huán)境科學(xué),2016,37(6):2322-2328.[LIU Ya-na,ZHU Shu-fa,WEI Xue-feng,etal.Assessment and Pollution Characteristics of Heavy Metals in Soil of Different Functional Areas in Luoyang[J].Environmental Science,2016,37(6):2322-2328.]

    [9] 劉衍君,張保華,劉子亭,等.基于ArcGIS的濰坊市土壤重金屬污染現(xiàn)狀與空間分布研究[J].水土保持通報(bào),2016,36(6):150-154.[LIU Yan-jun,ZHANG Bao-hua,LIU Zi-ting,etal.A Study on Distribution of Soil Heavy Metal Pollution in Weifang City Based on Geography Information System[J].Bulletin of Soil and Water Conservation,2016,36(6):150-154.]

    [10] 王佳佳,弓曉峰,向洪銳,等.城市內(nèi)湖疏浚前重金屬的沉積與污染特征研究[J].環(huán)境科學(xué)與技術(shù),2015,38(10):161-167.[WANG Jia-jia,GONG Xiao-feng,XIANG Hong-rui,etal.Study of Heavy Metals Deposition and Pollution Characteristics in Urban Lake Before Dredging[J].Environmental Science and Technology,2015,38(10):161-167.]

    [11] SHEN F,LIAO R M,ALI A,etal.Spatial Distribution and Risk Assessment of Heavy Metals in Soil Near a Pb/Zn Smelter in Feng County,China[J].Ecotoxicology and Environmental Safety,2017,139:254-262.

    [12] 趙倩,馬琳,劉翼飛,等.北京東南郊典型地層重金屬分布特征與潛在生態(tài)風(fēng)險(xiǎn)[J].環(huán)境科學(xué),2016,37(5):1931-1937.[ZHAO Qian,MA Lin,LIU Yi-fei,etal.Distribution Characteristics and Potential Ecological Hazards Assessment of Soil Heavy Metals in Typical Soil Profiles in Southeast Suburb of Beijing[J].Environmental Science,2016,37(5):1931-1937.]

    [13] 韓春梅,王林山,鞏宗強(qiáng),等.土壤中重金屬形態(tài)分析及其環(huán)境學(xué)意義[J].生態(tài)學(xué)雜志,2005,25(12):1499-1502.[HAN Chun-mei,WANG Lin-shan,GONG Zong-qiang,etal.Chemical Forms of Soil Heavy Metals and Their Environmental Significance[J].Chinese Journal of Ecology,2005,24(12):1499-1502.]

    [14] 周國(guó)華.土壤重金屬生物有效性研究進(jìn)展[J].物探與化探,2014,38(6):1097-1106.[ZHOU Guo-hua.Recent Progress in the Study of Heavy Metal Bioavailability in Soil[J].Geophysical & Geochemical Exploration,2014,38(6):1097-1106.]

    [15] TESSIER A,CAMPBELL P G C,BISSON M.Sequential Extraction Procedure for the Speciation of Particulate Trace Metal[J].Analytical Chemistry,1979,51(7):844-851.

    [16] RAURET G,LPEZ-SNCHEZ J F,SAHUQUILLO A,etal.Improvement of the BCR Three Step Sequential Extraction Procedure Prior to the Certification of New Sediment and Soil Reference Materials[J].Journal of Environmental Monitoring,1999,1(1):57-61.

    [17] 陳璐,文方,程艷,等.鉛鋅尾礦中重金屬形態(tài)分布與毒性浸出特征研究[J].干旱區(qū)資源與環(huán)境,2017,31(3):89-94.[CHEN Lu,WEN Fang,CHENG Yan,etal.Characteristics of Speciation Distribution and Toxicity Leaching of Heavy Metals in Pb-Zn Tailings[J].Journal of Arid Land Resources and Environment,2017,31(3):89-94.]

    [18] HJ/T 166—2004,土壤環(huán)境監(jiān)測(cè)技術(shù)規(guī)范[S].[HJ/T 166-2004,Technical Specification for Soil Environmental Monitoring[S].]

    [19] GB/T 22105.2—2008,土壤質(zhì)量 總汞、總砷、總鉛的測(cè)定 原子熒光法 第2部分:土壤中總砷的測(cè)定[S].[GB/T 22105.2—2008,Soil Quality:Analysis of Total Mercury,Arsenic and Lead Contents:Atomic Fluorescence Spectrometry:Part 2:Analysis of Total Arsenic Contents in Soils[S].]

    [20] HJ/T 332—2006,食用農(nóng)產(chǎn)品產(chǎn)地環(huán)境質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)[S].[HJ/T 332—2006,Environmental Quality Evaluation Standards for Farmland of Edible Agricultural Products[S].]

    [21] 周廣柱,楊鋒杰,程建光,等.土壤環(huán)境質(zhì)量綜合評(píng)價(jià)方法探討[J].山東科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,24(4):113-115.[ZHOU Guang-zhu,YANG Feng-jie,CHENG Jian-guang,etal.Research on Synthetic Assessment Method for Soil Environmental Quality[J].Journal of Shandong University of Science and Technology (Natural Science),2005,24(4):113-115.]

    [22] HAKANSON L.An Ecological Risk Index for Aquatic Pollution Control:A Sedimentological Approach[J].Water Research,1980,14(8):975-1001.

    [23] SOLTANI N,KESHAVARZI B,MOORE F,etal.Ecological and Human Health Hazards of Heavy Metals and Polycyclic Aromatic Hydrocarbons (PAHs) in Road Dust of Isfahan Metropolis,Iran[J].Science of the Total Environment,2015,505:712-723.

    [24] 范明毅,楊皓,黃先飛,等.典型山區(qū)燃煤型電廠周邊土壤重金屬形態(tài)特征及污染評(píng)價(jià)[J].中國(guó)環(huán)境科學(xué),2016,36(8):2425-2436.[FAN Ming-yi,YANG Hao,HUANG Xian-fei,etal.Chemical Forms and Risk Assessment of Heavy Metals in Soils Around a Typical Coal-Fired Power Plant Located in the Mountainous Area[J].China Environmental Science,2016,36(8):2425-2436.]

    [25] 李如忠,姜艷敏,潘成榮,等.典型有色金屬礦山城市小河流沉積物重金屬形態(tài)分布及風(fēng)險(xiǎn)評(píng)估[J].環(huán)境科學(xué),2013,34(3):1067-1075.[LI Ru-zhong,JIANG Yan-min,PAN Cheng-rong,etal.Fraction Distribution and Risk Assessment of Heavy Metals in Stream Sediments From a Typical Nonferrous Metals Mining City[J].Environmental Science,2013,34(3):1067-1075.]

    [26] 雷國(guó)建,陳志良,劉千鈞,等.廣州郊區(qū)土壤重金屬污染程度及潛在生態(tài)危害評(píng)價(jià)[J].中國(guó)環(huán)境科學(xué),2013,33(增刊1):49-53.[LEI Guo-jian,CHEN Zhi-liang,LIU Qian-jun,etal.The Assessments of Polluted Degree and Potential Ecological Hazards of Heavy Metals in Suburban Soil of Guangzhou City[J].China Environmental Science,2013,33(Suppl.1):49-53.]

    [27] GB 15618—1995,土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)[S].[GB 15618—1995,Environmental Quality Standard for Soils[S].]

    [28] 邵豐收,周皓韻.河南省主要元素的土壤環(huán)境背景值[J].河南農(nóng)業(yè),1998(10):29.[SHAO Feng-shou,ZHOU Hao-yun.The Main Elements of the Soil Environmental Background Value in Henan Province[J].Henan Agriculture,1998(10):29.]

    [29] 張金蓮,丁疆峰,盧桂寧,等.廣東清遠(yuǎn)電子垃圾拆解區(qū)農(nóng)田土壤重金屬污染評(píng)價(jià)[J].環(huán)境科學(xué),2015,36(7):2633-2640.[ZHANG Jin-lian,DING Jiang-feng,LU Gui-ning,etal.Heavy Metal Contamination in Farmland Soils at an E-Waste Disassembling Site in Qingyuan,Guangdong,South China[J].Environmental Science,2015,36(7):2633-2640.]

    [30] 李瑛,張桂銀,李洪軍,等.根際土壤中Cd、Pb形態(tài)轉(zhuǎn)化及其植物效應(yīng)[J].生態(tài)環(huán)境,2004,13(3):316-319.[LI Ying,ZHANG Gui-yin,LI Hong-jun,etal.Speciation Transformation of Cd,Pb in Rhizosphere and Their Effects on Plant[J].Ecology and Environment,2004,13(3):316-319.]

    [31] MERTENS J,VERVAEKE P,SCHRIJVER A D,etal.Metal Uptake by Young Trees From Dredged Brackish Sediment:Limitations and Possibilities for Phytoextraction and Phytostabilisation[J].Science of the Total Environment,2004,326(1/2/3):209-215.

    [32] 李?yuàn)檴?曹廣超,石平超,等.青島城區(qū)土壤重金屬元素空間分布及其現(xiàn)狀評(píng)價(jià)[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2015,31(1):112-117.[LI Shan-shan,CAO Guang-chao,SHI Ping-chao,etal.Status Quo and Evaluation of the Spatial Distribution of Heavy Metals in Urban Soil of Qingdao City[J].Journal of Ecology and Rural Environment,2015,31(1):112-117.]

    [33] 姜俊彥,黃星,李秀珍,等.潮灘濕地土壤有機(jī)碳儲(chǔ)量及其與土壤理化因子的關(guān)系:以崇明東灘為例[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2015,31(4):540-547.[JIANG Jun-yan,HUANG Xing,LI Xiu-zhen,etal.Soil Organic Carbon Storage in Tidal Wetland and Its Relationships With Soil Physico-Chemical Factors:A Case Study of Dongtan of Chongming,Shanghai[J].Journal of Ecology and Rural Environment,2015,31(4):540-547.]

    猜你喜歡
    弱酸全量重金屬
    運(yùn)用Ka與Ksp計(jì)算分析難溶弱酸鹽的溶解與沉淀
    成都市溫江區(qū)全力推進(jìn)醫(yī)保全量數(shù)據(jù)采集試點(diǎn)工作
    重金屬對(duì)膨潤(rùn)土膨脹性的影響
    測(cè)定不同產(chǎn)地寬筋藤中5種重金屬
    中成藥(2018年8期)2018-08-29 01:28:16
    多元弱酸的準(zhǔn)確滴定
    土壤重金屬污染調(diào)查與評(píng)估的誤區(qū)及其改進(jìn)方法
    ICP-AES、ICP-MS測(cè)定水中重金屬的對(duì)比研究
    麥秸全量還田下氮肥運(yùn)籌對(duì)水稻產(chǎn)量及其產(chǎn)量構(gòu)成的影響
    再生水回灌中DOM對(duì)重金屬遷移與保留問(wèn)題研究
    水稻冠層光譜變化特征的土壤重金屬全量反演研究
    遙感信息(2015年3期)2015-12-13 07:26:54
    国产 一区精品| 99在线视频只有这里精品首页| 久久精品国产清高在天天线| 国产蜜桃级精品一区二区三区| 97碰自拍视频| 黄色一级大片看看| 中文资源天堂在线| 欧美色视频一区免费| 欧美精品国产亚洲| 中文字幕av成人在线电影| 国产精品亚洲一级av第二区| 一个人看的www免费观看视频| 联通29元200g的流量卡| 国产大屁股一区二区在线视频| 色视频www国产| www.www免费av| 少妇裸体淫交视频免费看高清| 亚洲狠狠婷婷综合久久图片| 美女高潮喷水抽搐中文字幕| 色综合亚洲欧美另类图片| 久久99热这里只有精品18| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久亚洲av鲁大| 亚洲经典国产精华液单| 国语自产精品视频在线第100页| 久久久久久久久中文| 亚洲欧美日韩高清在线视频| 成人一区二区视频在线观看| 大又大粗又爽又黄少妇毛片口| 成人欧美大片| 99精品久久久久人妻精品| 国产成年人精品一区二区| 人人妻人人澡欧美一区二区| 精品人妻1区二区| 日本五十路高清| 久久久久久久久久成人| 99久久成人亚洲精品观看| 三级毛片av免费| 美女cb高潮喷水在线观看| 嫩草影院新地址| 有码 亚洲区| 丰满人妻一区二区三区视频av| 少妇丰满av| 最新中文字幕久久久久| 中文亚洲av片在线观看爽| 真人做人爱边吃奶动态| 99热这里只有精品一区| 欧美色视频一区免费| 成人午夜高清在线视频| 亚洲av成人av| 少妇丰满av| 精品人妻偷拍中文字幕| h日本视频在线播放| 赤兔流量卡办理| 国产久久久一区二区三区| 熟女人妻精品中文字幕| 免费人成视频x8x8入口观看| 国产一区二区亚洲精品在线观看| 久久精品夜夜夜夜夜久久蜜豆| 99久久久亚洲精品蜜臀av| 免费观看在线日韩| 麻豆成人av在线观看| 国产精品乱码一区二三区的特点| 亚洲午夜理论影院| 亚洲一区高清亚洲精品| 午夜爱爱视频在线播放| 午夜福利在线观看吧| 欧美xxxx黑人xx丫x性爽| 亚洲欧美清纯卡通| 亚洲最大成人手机在线| 91久久精品国产一区二区成人| 嫩草影院精品99| 最近视频中文字幕2019在线8| 国产精品自产拍在线观看55亚洲| 麻豆国产97在线/欧美| av国产免费在线观看| 一进一出抽搐gif免费好疼| 小说图片视频综合网站| 久久久久久九九精品二区国产| 亚洲 国产 在线| 99国产极品粉嫩在线观看| 在线播放无遮挡| 久久久久久久亚洲中文字幕| 日本色播在线视频| 亚洲国产精品合色在线| 99久久九九国产精品国产免费| 麻豆成人av在线观看| 两个人的视频大全免费| 又紧又爽又黄一区二区| 永久网站在线| av在线蜜桃| 免费观看在线日韩| av福利片在线观看| 欧美一区二区精品小视频在线| 国产三级在线视频| 国产精品久久久久久精品电影| 欧美精品国产亚洲| 亚洲自偷自拍三级| 日日撸夜夜添| 五月玫瑰六月丁香| 变态另类丝袜制服| 免费电影在线观看免费观看| 88av欧美| 亚洲av免费在线观看| 有码 亚洲区| 国产精品自产拍在线观看55亚洲| 日本精品一区二区三区蜜桃| 久久99热这里只有精品18| 亚洲美女视频黄频| 在线观看一区二区三区| 黄片wwwwww| 97人妻精品一区二区三区麻豆| 2021天堂中文幕一二区在线观| 午夜福利高清视频| 国产av在哪里看| 国产午夜精品久久久久久一区二区三区 | 欧洲精品卡2卡3卡4卡5卡区| 日本黄色片子视频| 国产成人av教育| 男人和女人高潮做爰伦理| 在线观看午夜福利视频| 免费搜索国产男女视频| 超碰av人人做人人爽久久| 欧美性猛交黑人性爽| 国产av不卡久久| 亚洲自偷自拍三级| 久久久久久久久久黄片| 男女做爰动态图高潮gif福利片| 97人妻精品一区二区三区麻豆| 国产男靠女视频免费网站| ponron亚洲| 国产精品一区二区性色av| 身体一侧抽搐| 麻豆成人av在线观看| 亚洲中文字幕日韩| 欧美又色又爽又黄视频| 五月伊人婷婷丁香| 亚洲精品色激情综合| 亚洲综合色惰| 3wmmmm亚洲av在线观看| 日韩精品青青久久久久久| 婷婷色综合大香蕉| 亚洲av成人av| 夜夜爽天天搞| 精品一区二区三区视频在线| 日韩欧美精品免费久久| 欧美色视频一区免费| 啦啦啦啦在线视频资源| 精品久久久久久久久久免费视频| 成人高潮视频无遮挡免费网站| 欧美不卡视频在线免费观看| 久久久久免费精品人妻一区二区| 午夜日韩欧美国产| 熟妇人妻久久中文字幕3abv| 日本成人三级电影网站| 在线观看av片永久免费下载| 少妇被粗大猛烈的视频| 亚洲综合色惰| 美女高潮喷水抽搐中文字幕| 99在线人妻在线中文字幕| 国产av不卡久久| 亚洲精品粉嫩美女一区| 亚洲内射少妇av| 国产黄色小视频在线观看| 最近在线观看免费完整版| 老司机福利观看| 国产高清三级在线| 日韩欧美在线乱码| 国产老妇女一区| 亚洲成人久久性| 午夜福利欧美成人| 亚洲国产欧美人成| 性色avwww在线观看| 成年免费大片在线观看| 国产aⅴ精品一区二区三区波| 真人做人爱边吃奶动态| 国产欧美日韩精品一区二区| av专区在线播放| 特级一级黄色大片| 国产一区二区三区视频了| 18禁黄网站禁片免费观看直播| 亚洲国产精品sss在线观看| eeuss影院久久| 国产午夜精品久久久久久一区二区三区 | 99在线人妻在线中文字幕| 亚洲精品国产成人久久av| 性色avwww在线观看| 少妇人妻精品综合一区二区 | 日日摸夜夜添夜夜添小说| 在线国产一区二区在线| 特级一级黄色大片| 国产白丝娇喘喷水9色精品| 日本成人三级电影网站| 色哟哟哟哟哟哟| 午夜福利在线在线| 婷婷色综合大香蕉| 男女下面进入的视频免费午夜| 日韩欧美三级三区| 欧美在线一区亚洲| 免费人成视频x8x8入口观看| 日韩欧美精品免费久久| 亚洲精品影视一区二区三区av| 狂野欧美白嫩少妇大欣赏| 高清在线国产一区| 免费av观看视频| 国产午夜精品论理片| 免费无遮挡裸体视频| 成人av一区二区三区在线看| 国产91精品成人一区二区三区| 超碰av人人做人人爽久久| 可以在线观看的亚洲视频| 免费黄网站久久成人精品| 国产成人a区在线观看| 国产伦在线观看视频一区| 动漫黄色视频在线观看| 午夜福利欧美成人| 欧美一区二区国产精品久久精品| 欧美日韩乱码在线| 亚洲aⅴ乱码一区二区在线播放| 少妇猛男粗大的猛烈进出视频 | 熟女人妻精品中文字幕| 蜜桃亚洲精品一区二区三区| 中文在线观看免费www的网站| 亚洲av中文字字幕乱码综合| 国产免费男女视频| 精品人妻视频免费看| 亚洲18禁久久av| 99久久无色码亚洲精品果冻| 两人在一起打扑克的视频| 精品国内亚洲2022精品成人| 中文亚洲av片在线观看爽| 精品久久久久久,| 十八禁网站免费在线| 最近最新免费中文字幕在线| 成人二区视频| 淫妇啪啪啪对白视频| 一个人看的www免费观看视频| 午夜免费成人在线视频| 免费高清视频大片| 99riav亚洲国产免费| 特大巨黑吊av在线直播| 人妻丰满熟妇av一区二区三区| 国产激情偷乱视频一区二区| 欧美三级亚洲精品| 在线观看一区二区三区| 一个人免费在线观看电影| 亚洲成人久久爱视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美一级a爱片免费观看看| 日韩一区二区视频免费看| 国语自产精品视频在线第100页| 亚洲精品亚洲一区二区| 久久久成人免费电影| 国产精品自产拍在线观看55亚洲| 亚洲五月天丁香| 三级毛片av免费| 老师上课跳d突然被开到最大视频| 在线观看66精品国产| 亚洲美女黄片视频| 亚洲av免费在线观看| 少妇裸体淫交视频免费看高清| 高清毛片免费观看视频网站| 国产视频内射| 最近最新中文字幕大全电影3| 人人妻人人看人人澡| 看十八女毛片水多多多| 在线观看av片永久免费下载| 久久精品夜夜夜夜夜久久蜜豆| 成年女人永久免费观看视频| 免费搜索国产男女视频| 51国产日韩欧美| 亚洲不卡免费看| 国产激情偷乱视频一区二区| 变态另类成人亚洲欧美熟女| 日韩一本色道免费dvd| 高清日韩中文字幕在线| 久久精品综合一区二区三区| 午夜精品久久久久久毛片777| 99热精品在线国产| 亚洲一区二区三区色噜噜| 最好的美女福利视频网| 精品午夜福利在线看| 99riav亚洲国产免费| 国产久久久一区二区三区| 亚洲色图av天堂| av国产免费在线观看| 啪啪无遮挡十八禁网站| 老司机福利观看| 久久久久久久亚洲中文字幕| av福利片在线观看| 丰满人妻一区二区三区视频av| 中亚洲国语对白在线视频| av中文乱码字幕在线| 久久中文看片网| 午夜老司机福利剧场| 99久国产av精品| 色综合婷婷激情| 日韩欧美 国产精品| 国产欧美日韩精品亚洲av| 国产成人影院久久av| 淫妇啪啪啪对白视频| 12—13女人毛片做爰片一| 国产一区二区在线av高清观看| 国产成人aa在线观看| 九色国产91popny在线| 成年版毛片免费区| 亚洲成人中文字幕在线播放| 精品久久国产蜜桃| 成人性生交大片免费视频hd| 国产高清三级在线| 波野结衣二区三区在线| 亚州av有码| 国产淫片久久久久久久久| 国产欧美日韩精品一区二区| 最好的美女福利视频网| 嫩草影院新地址| 欧美色欧美亚洲另类二区| 国产精品电影一区二区三区| 在线观看免费视频日本深夜| 国产精品免费一区二区三区在线| 高清毛片免费观看视频网站| 九色国产91popny在线| 欧美区成人在线视频| 免费av毛片视频| 最近在线观看免费完整版| 中文字幕av在线有码专区| 无遮挡黄片免费观看| 一区二区三区四区激情视频 | 日韩欧美在线二视频| 免费av观看视频| 久久人人精品亚洲av| 亚洲aⅴ乱码一区二区在线播放| 免费黄网站久久成人精品| 国语自产精品视频在线第100页| 国产高清视频在线观看网站| 久久九九热精品免费| 亚洲天堂国产精品一区在线| 亚洲国产精品合色在线| 亚洲自拍偷在线| 俄罗斯特黄特色一大片| 成年免费大片在线观看| 欧美日韩国产亚洲二区| 婷婷精品国产亚洲av在线| 91久久精品国产一区二区三区| 久久久国产成人免费| 九色国产91popny在线| 永久网站在线| 综合色av麻豆| 欧美另类亚洲清纯唯美| 白带黄色成豆腐渣| 欧美日本亚洲视频在线播放| 久久久久久大精品| 成人性生交大片免费视频hd| 国产色婷婷99| 欧美三级亚洲精品| 国产91精品成人一区二区三区| 中亚洲国语对白在线视频| 免费黄网站久久成人精品| 国产毛片a区久久久久| 村上凉子中文字幕在线| 噜噜噜噜噜久久久久久91| 51国产日韩欧美| 麻豆精品久久久久久蜜桃| 亚洲精品乱码久久久v下载方式| 在线播放无遮挡| 男女下面进入的视频免费午夜| 成人特级av手机在线观看| 国产中年淑女户外野战色| 免费一级毛片在线播放高清视频| 美女大奶头视频| 欧美成人a在线观看| 级片在线观看| 日本黄色视频三级网站网址| 床上黄色一级片| 又爽又黄无遮挡网站| 又紧又爽又黄一区二区| 亚洲av成人av| 一级av片app| 网址你懂的国产日韩在线| 成人国产麻豆网| 欧美激情在线99| 别揉我奶头 嗯啊视频| 我的老师免费观看完整版| 欧美一区二区亚洲| 此物有八面人人有两片| 国产精品国产三级国产av玫瑰| 精品人妻熟女av久视频| 精品久久久久久久人妻蜜臀av| 婷婷精品国产亚洲av在线| 99九九线精品视频在线观看视频| 欧美中文日本在线观看视频| 欧美不卡视频在线免费观看| 久久久久久伊人网av| 欧美色视频一区免费| 看黄色毛片网站| 亚洲欧美日韩无卡精品| 欧美zozozo另类| 久久久精品大字幕| 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩卡通动漫| 国产成人a区在线观看| 午夜福利18| 69av精品久久久久久| 尤物成人国产欧美一区二区三区| 国产老妇女一区| 亚洲自偷自拍三级| 在现免费观看毛片| 国产伦精品一区二区三区四那| 他把我摸到了高潮在线观看| 最后的刺客免费高清国语| 久久精品国产亚洲av天美| 黄色日韩在线| 免费电影在线观看免费观看| 国产淫片久久久久久久久| 午夜免费男女啪啪视频观看 | 国产伦一二天堂av在线观看| 亚洲精品日韩av片在线观看| 中亚洲国语对白在线视频| 精品一区二区三区av网在线观看| 老熟妇仑乱视频hdxx| 淫秽高清视频在线观看| 欧美丝袜亚洲另类 | 欧美日韩中文字幕国产精品一区二区三区| 在线免费观看不下载黄p国产 | 久久久久久久午夜电影| 女人被狂操c到高潮| 精品福利观看| 国产亚洲91精品色在线| 国产亚洲欧美98| 黄色日韩在线| 亚洲精品粉嫩美女一区| 亚洲久久久久久中文字幕| 我的老师免费观看完整版| 一个人观看的视频www高清免费观看| 波多野结衣巨乳人妻| 国产精品国产三级国产av玫瑰| 女人被狂操c到高潮| 精品不卡国产一区二区三区| 亚洲精品色激情综合| 在线看三级毛片| 女同久久另类99精品国产91| 嫩草影院精品99| 九色国产91popny在线| 搡老岳熟女国产| 热99re8久久精品国产| 日本-黄色视频高清免费观看| 久久久久久久精品吃奶| 日韩欧美在线乱码| 国产私拍福利视频在线观看| 能在线免费观看的黄片| 少妇高潮的动态图| 欧美日韩中文字幕国产精品一区二区三区| 天堂影院成人在线观看| 欧美极品一区二区三区四区| 91在线观看av| 欧美黑人欧美精品刺激| 男人狂女人下面高潮的视频| 精品一区二区三区视频在线观看免费| 欧美成人性av电影在线观看| 真人做人爱边吃奶动态| 色综合亚洲欧美另类图片| 久久精品国产亚洲网站| 国产精品久久久久久久电影| 69人妻影院| 中文字幕av在线有码专区| 亚洲一级一片aⅴ在线观看| 99热这里只有是精品在线观看| 在线观看66精品国产| 国产精品女同一区二区软件 | 国产一区二区在线观看日韩| 午夜激情欧美在线| 国产极品精品免费视频能看的| av视频在线观看入口| 欧美潮喷喷水| 久久精品夜夜夜夜夜久久蜜豆| 在线免费观看不下载黄p国产 | 制服丝袜大香蕉在线| 天堂av国产一区二区熟女人妻| 国产精品一及| 久久精品国产亚洲av天美| 久久人人精品亚洲av| 一级黄片播放器| 免费在线观看影片大全网站| 色精品久久人妻99蜜桃| 小蜜桃在线观看免费完整版高清| 国产精品av视频在线免费观看| 亚洲人成伊人成综合网2020| 成人特级黄色片久久久久久久| 嫩草影院新地址| 小蜜桃在线观看免费完整版高清| 91麻豆精品激情在线观看国产| 又爽又黄无遮挡网站| 国产精品综合久久久久久久免费| 免费搜索国产男女视频| 国产亚洲精品综合一区在线观看| 99热网站在线观看| 成人二区视频| 日本免费一区二区三区高清不卡| 日韩欧美精品免费久久| 白带黄色成豆腐渣| 91在线精品国自产拍蜜月| 中国美白少妇内射xxxbb| 免费高清视频大片| 亚洲欧美精品综合久久99| 久久久久久九九精品二区国产| 久久久久久久精品吃奶| 看片在线看免费视频| 亚洲综合色惰| 成年女人毛片免费观看观看9| 免费在线观看影片大全网站| 亚洲国产精品成人综合色| 少妇的逼水好多| 久久精品国产清高在天天线| 干丝袜人妻中文字幕| 22中文网久久字幕| 日韩一本色道免费dvd| 久久6这里有精品| 日韩在线高清观看一区二区三区 | 岛国在线免费视频观看| 成人欧美大片| 91麻豆精品激情在线观看国产| 男人狂女人下面高潮的视频| 一进一出抽搐gif免费好疼| 久久午夜福利片| 免费看光身美女| 日韩强制内射视频| 亚洲在线自拍视频| 国产伦精品一区二区三区视频9| 国产亚洲欧美98| 22中文网久久字幕| 国产精华一区二区三区| 亚洲色图av天堂| 天美传媒精品一区二区| 在线国产一区二区在线| 88av欧美| 久久精品国产鲁丝片午夜精品 | 在线看三级毛片| 人妻夜夜爽99麻豆av| 日韩高清综合在线| 欧美日韩瑟瑟在线播放| 制服丝袜大香蕉在线| av在线亚洲专区| 天美传媒精品一区二区| 成人亚洲精品av一区二区| 丰满的人妻完整版| 真实男女啪啪啪动态图| 美女黄网站色视频| 亚洲第一电影网av| 免费av毛片视频| 别揉我奶头 嗯啊视频| 亚洲美女视频黄频| 亚洲精品乱码久久久v下载方式| 亚洲乱码一区二区免费版| 午夜激情欧美在线| 91久久精品电影网| 最近在线观看免费完整版| 日韩欧美 国产精品| xxxwww97欧美| 3wmmmm亚洲av在线观看| 久久久久久久久久久丰满 | 免费看美女性在线毛片视频| 日韩在线高清观看一区二区三区 | 两人在一起打扑克的视频| 欧美+亚洲+日韩+国产| 午夜激情欧美在线| 午夜免费成人在线视频| 午夜激情福利司机影院| 久久久久久九九精品二区国产| 久9热在线精品视频| 麻豆久久精品国产亚洲av| 亚洲av免费高清在线观看| 三级国产精品欧美在线观看| 成人性生交大片免费视频hd| 国内少妇人妻偷人精品xxx网站| 国产麻豆成人av免费视频| 此物有八面人人有两片| 国产中年淑女户外野战色| 又粗又爽又猛毛片免费看| 国产亚洲精品久久久com| 国语自产精品视频在线第100页| 亚洲av熟女| 欧美绝顶高潮抽搐喷水| 国产一区二区三区在线臀色熟女| 人人妻人人看人人澡| 欧美性猛交黑人性爽| 亚洲精华国产精华液的使用体验 | 老司机深夜福利视频在线观看| 一区二区三区四区激情视频 | 两个人视频免费观看高清| 久久久久国产精品人妻aⅴ院| 日韩av在线大香蕉| 亚洲av二区三区四区| 狠狠狠狠99中文字幕| 亚洲国产精品合色在线| 免费人成视频x8x8入口观看| 日韩精品中文字幕看吧| 丝袜美腿在线中文| 日本成人三级电影网站| 国产亚洲av嫩草精品影院| 老熟妇仑乱视频hdxx| x7x7x7水蜜桃| 国产一区二区三区av在线 | 美女免费视频网站| 校园春色视频在线观看| 国产蜜桃级精品一区二区三区| 高清在线国产一区| 丰满人妻一区二区三区视频av| 午夜久久久久精精品| 我的女老师完整版在线观看| 九色国产91popny在线| 哪里可以看免费的av片| 人人妻人人看人人澡| 超碰av人人做人人爽久久| 欧美不卡视频在线免费观看| 在线观看午夜福利视频|