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

    基于同步輻射技術(shù)研究鋅在廢棄稀土礦區(qū)土壤中的吸附及垂向遷移特征

    2025-03-02 00:00:00李洋歐麗黃琨劉鵬劉燕紅李曼曼唐金平
    關(guān)鍵詞:吸附遷移土壤

    摘要: 為了研究鋅(Zn)在贛南稀土礦區(qū)土壤的吸附及垂向遷移特征,以及為礦山開(kāi)采地區(qū)土壤及水體Zn污染防治、修復(fù)及風(fēng)險(xiǎn)評(píng)估提供理論依據(jù),首先采用靜態(tài)批實(shí)驗(yàn),研究某廢棄稀土礦區(qū)垂直剖面不同深度土壤對(duì)Zn的吸附特征;然后采用X射線衍射(XRD)及X射線吸收譜(XAFS)研究土壤主要物相組成及土壤中Zn的賦存形態(tài)變化,探討Zn在土壤中的吸附機(jī)制及可遷移性。結(jié)果表明:剖面3層土壤對(duì)Zn的吸附平衡過(guò)程較快,土壤對(duì)Zn的吸附能力呈現(xiàn)出從表層土到深層土再到中層土逐漸減弱的趨勢(shì),準(zhǔn)二級(jí)動(dòng)力學(xué)模型及Langmuir方程能較好地描述土壤吸附Zn的過(guò)程;元素Fe、Mn在土壤中主要以Fe的氫氧化物及混合價(jià)態(tài)錳氧化物形式存在,剖面各層土壤理化性質(zhì)差異影響土壤對(duì)Zn的吸附行為,相比于黏土礦物,F(xiàn)e、Mn氧化物對(duì)土壤吸附Zn的影響更加顯著;吸附反應(yīng)前后土壤中Zn的局部配位模式由Zn-O四面體轉(zhuǎn)變?yōu)閆n-O四面體和八面體的混合配位模式,反應(yīng)后土壤表面可能形成了類似硫酸鹽的外層絡(luò)合物,反應(yīng)前后土壤中的Zn由堿式碳酸鹽轉(zhuǎn)變?yōu)橐粤蛩猁}與堿式碳酸鹽的混合形式存在;受土壤pH影響,土壤滯留Zn能力有限,土壤中的有效態(tài)Zn易伴隨降雨沖刷向下層土壤遷移,土壤中Zn質(zhì)量分?jǐn)?shù)隨深度增加而增加。

    關(guān)鍵詞:土壤;鋅;吸附;遷移;X射線吸收譜;稀土礦區(qū);同步輻射

    doi:10.13278/j.cnki.jjuese.20230110

    中圖分類號(hào):X53 文獻(xiàn)標(biāo)志碼:A

    收稿日期: 2023-04-29

    作者簡(jiǎn)介: 李洋(1999-),女,碩士研究生,主要從事地下水重金屬污染防治方面的研究,E-mail:melilonli@outlook.com

    通信作者: 劉鵬(1981-),男,教授,博士生導(dǎo)師,主要從事地下水污染防治方面的研究,E-mail: pengliu@cug.edu.cn

    基金項(xiàng)目: 江西省工業(yè)園區(qū)地下水監(jiān)測(cè)管控網(wǎng)絡(luò)體系建設(shè)示范項(xiàng)目(JXMD2022-Z002);中央環(huán)保投資儲(chǔ)備庫(kù)項(xiàng)目(2022B000563260)

    Supported by the Demonstration Project of Groundwater Control Network System Construction in Industrial Park of Jiangxi Province (JXMD2022-Z002) and the Central Investment Reserve Project for Environmental Protection(2022B000563260)

    Adsorption and Vertical Migration Characteristics" of Zinc in Soil from" Abandoned Rare Earth Mine Using Synchrotron Radiation-Based Technique

    Li Yang1, Ou Li2, Huang Kun1, Liu Peng1, Liu Yanhong2, Li Manman1, Tang Jinping1

    1. College of Environment, China University of Geosciences, Wuhan 430074, China

    2. Ecological Environmental Monitoring Centre of Jiangxi Province, Nanchang 330039, China

    Abstract: This study focuses on the adsorption characteristics of zinc in soil at different depths of a certain abandoned rare earth mine in the south of Ganzhou, Jiangxi Province. Static batch experiments were conducted. X-ray diffraction (XRD) and X-ray absorption fine structure spectroscopy (XAFS) were used to study the major phases and morphological changes of zinc in soil before and after adsorption reaction, in order to investigate the adsorption mechanism and mobility of zinc in soil. The results show that the adsorption equilibrium process of zinc in soil is fast. The adsorption capacity of soil for zinc gradually decreased from the surface soil to the deep soil and then to the middle soil. The pseudo-second-order kinetic model and Langmuir equation can describe the adsorption process of zinc in soil well. XRD and XAFS results show that iron and manganese exist mainly in the form of iron hydroxide and mixed valence state manganese oxide in soil. Differences in the physical and chemical properties of soil layers affect the adsorption behavior of soil towards zinc. Compared with clay minerals, iron and manganese oxides have a more significant effect on the adsorption of zinc in soil. XAFS indicates that the local coordination mode of zinc changes from Zn-O tetrahedron to a mixed coordination mode of Zn-O tetrahedron and octahedron after adsorption reaction. A complex similar to sulfate may form on the soil surface after the reaction. The form of zinc in soil changes from basic carbonate to a mixed form of sulfate and basic carbonate after adsorption reaction. Due to the influence of soil pH, the capacity of zinc retention in soil is limited. The readily available zinc in soil is prone to migrate downward with rainfall. The zinc content in soil increases with depth. The study provides a scientific basis for the prevention and control of zinc pollution in soil and groundwater in mining areas.

    Key words: soil; zinc; adsorption; migration; X-ray absorption fine structure spectroscopy; rare earth mine; synchrotron radiation

    0 引言

    鋅(Zn)是動(dòng)植物生長(zhǎng)發(fā)育所必須的微量元素,然而過(guò)量的Zn也可能對(duì)生物體具有毒性[1]。贛南地區(qū)是我國(guó)離子型稀土礦開(kāi)采的重點(diǎn)地區(qū)。開(kāi)采過(guò)程中產(chǎn)生的酸性廢水易攜帶Zn2+等離子進(jìn)入土壤[2],造成土壤酸化,或?qū)е峦寥乐兄亟饘匐x子活化[3],而且在降雨沖刷和淋濾作用下Zn2+等離子會(huì)可能向下層土壤遷移[4],對(duì)周邊地區(qū)土壤及水生生態(tài)系統(tǒng)產(chǎn)生不利影響[5]。土壤是由固相及流體相共同組成的復(fù)雜系統(tǒng)[6],進(jìn)入土壤系統(tǒng)的重金屬可與各相介質(zhì)之間發(fā)生交互作用[7],該過(guò)程受不同反應(yīng)類型控制,如吸附-解吸作用、溶解-沉淀作用及氧化-還原作用等[8-9]。其中,吸附作為重金屬在土壤中循環(huán)的主要過(guò)程[10],能夠引起土壤中重金屬賦存形態(tài)變化,是影響重金屬在環(huán)境中遷移、生物有效性及毒理性的關(guān)鍵因素[11]。開(kāi)展重金屬在土壤組分上宏觀吸附熱力學(xué)、動(dòng)力學(xué)行為和分子機(jī)制的研究是評(píng)價(jià)重金屬環(huán)境行為的基礎(chǔ),對(duì)于預(yù)測(cè)土壤中Zn的遷移性、開(kāi)展土壤及水相Zn污染防治具有重要意義。

    長(zhǎng)期以來(lái),由于分析儀器和實(shí)驗(yàn)手段的局限性,對(duì)Zn在土壤及其典型固相組分中的化學(xué)行為研究多集中在宏觀動(dòng)力學(xué)和熱力學(xué)層面。宏觀吸附結(jié)果表明土壤對(duì)Zn的吸附行為與土壤pH、鐵錳氧化物、黏土礦物等密切相關(guān)[12-13]。近年來(lái),隨著同步輻射技術(shù)的快速發(fā)展,X射線吸收譜(XAFS)技術(shù)在研究固-液界面吸附態(tài)分子微觀結(jié)構(gòu)方面得到了實(shí)際應(yīng)用,研究人員利用XAFS技術(shù)從分子水平上獲取重金屬在固-液界面的化學(xué)形態(tài)和結(jié)合方式等信息[14],研究了重金屬在土壤組分的界面化學(xué)行為[15]。有學(xué)者[16]通過(guò)開(kāi)展靜態(tài)批實(shí)驗(yàn)結(jié)合XAFS研究發(fā)現(xiàn),不同理化性質(zhì)土壤吸附Zn2+的能力與土壤pH呈顯著正相關(guān),吸附反應(yīng)后土壤中Zn的構(gòu)型為Zn-O八面體(六配位)和Zn-O四面體(四配位)的混合物。Fan等[17]利用XAFS結(jié)合通徑分析手段研究了Zn2+在江蘇省土壤中的吸附機(jī)理,發(fā)現(xiàn)土壤pH、黏土礦物是影響土壤吸附Zn的主要因素,Zn以Zn-Al層狀氫氧化物(LDH)、ZnCO3、Zn2(OH)2CO3沉淀及其與黏土礦物或鐵氧化物相結(jié)合的形式滯留于堿性土壤中。目前有關(guān)XAFS技術(shù)應(yīng)用于土壤重金屬吸附行為的研究多集中在單一礦物或單層土壤對(duì)Zn吸附行為的探討[10, 17-19]方面,但對(duì)Zn在同一剖面不同深度土壤中的吸附及垂向遷移特征的研究較少。礦區(qū)土壤大多經(jīng)過(guò)回填,填土孔隙率較高,有利于重金屬垂向下滲遷移[20],研究Zn在剖面不同深度土壤中的吸附及遷移特征更能真實(shí)評(píng)價(jià)礦區(qū)土壤中重金屬Zn的環(huán)境安全性。

    本研究以贛南地區(qū)某廢棄稀土礦區(qū)域垂直剖面土壤為對(duì)象,首先通過(guò)開(kāi)展靜態(tài)批實(shí)驗(yàn),研究不同反應(yīng)時(shí)間及初始Zn濃度條件下剖面土壤對(duì)Zn的吸附特征;然后結(jié)合X射線衍射(XRD)及X射線吸收精細(xì)結(jié)構(gòu)譜(XAFS)等手段,分析土壤主要物相組成及吸附反應(yīng)前后土壤中Zn的賦存形態(tài),揭示Zn在土壤中的吸附固持機(jī)制及垂向遷移特征,以期為礦山開(kāi)采地區(qū)土壤及水體Zn污染防治、修復(fù)及風(fēng)險(xiǎn)評(píng)估提供理論依據(jù)。

    1 材料與方法

    1.1 實(shí)驗(yàn)材料

    土壤樣品采自江西省贛南某廢棄稀土礦區(qū)剖面土層,從上至下依次為灰色雜填土層S1(0~20 cm)、灰黃色風(fēng)化巖層S2(350~450 cm)及灰色風(fēng)化巖層S3(700~750 cm)。土壤樣品采集后經(jīng)混合、自然風(fēng)干及研磨后,過(guò)100目篩保存?zhèn)溆?。土壤樣品的基本理化性質(zhì)見(jiàn)表1。

    1.2 實(shí)驗(yàn)設(shè)計(jì)

    1.2.1 Zn的吸附動(dòng)力學(xué)實(shí)驗(yàn)

    分別稱取剖面土壤樣品1.0 g于離心管中,以0.01 mol/L的NaCl2溶液作為背景電解質(zhì),水土質(zhì)量比為20∶1,加入初始Zn質(zhì)量濃度為10 mg/L的ZnCl2溶液20 mL,在25℃下避光振蕩后分別于0、2、5、10、30、60、120、240、480、720、1 440 min后,取上清液過(guò)濾(孔徑0.22 μm水系濾膜)后使用火焰原子吸收分光光度計(jì)(GGX-830)測(cè)試溶液Zn濃度,得到吸附平衡時(shí)間作為等溫吸附實(shí)驗(yàn)的反應(yīng)時(shí)間。在不同的實(shí)驗(yàn)條件下設(shè)置空白實(shí)驗(yàn)以排除土壤中固有吸附態(tài)Zn對(duì)實(shí)驗(yàn)結(jié)果的影響。

    1.2.2 Zn的吸附熱力學(xué)實(shí)驗(yàn)

    每組實(shí)驗(yàn)設(shè)置6個(gè)Zn初始質(zhì)量濃度(5、10、20、50、100、200 mg/L),在動(dòng)力學(xué)實(shí)驗(yàn)確定的平衡時(shí)間(8 h)取樣,其他實(shí)驗(yàn)條件及步驟同吸附動(dòng)力學(xué)實(shí)驗(yàn)一致。

    1.3 樣品分析

    土壤pH值采用電位法測(cè)定,陽(yáng)離子交換量(CEC)采用滴定法測(cè)定,有機(jī)質(zhì)質(zhì)量濃度采用重鉻酸鉀氧化-容量法測(cè)定。土壤重金屬元素總量采用電感耦合等離子體質(zhì)譜儀(ICP-MS)測(cè)定。

    土壤樣品冷凍干燥、研磨并過(guò)200目篩后,利用X射線衍射儀(XRD,D8 Advance,德國(guó))分析土壤物相組成,工作電壓為40 kV,電流為40 mA,步長(zhǎng)范圍為10 °~80 °。

    吸附動(dòng)力學(xué)實(shí)驗(yàn)中反應(yīng)24 h后離心收集土樣放入真空干燥機(jī)中烘干保存,即為反應(yīng)后的土壤樣品(#S1、#S2、#S3)。采用同步輻射X射線吸收精細(xì)結(jié)構(gòu)(XAFS)譜分析原始土壤中Fe、Mn形態(tài)及反應(yīng)前后土壤中的Zn形態(tài)。將土壤樣品冷凍干燥后,在手套箱(COY-7000220A,美國(guó))中過(guò)2 mm篩并除去其中的生物碎屑,研磨成粉末(200目篩)后壓實(shí)填入聚四氟乙烯板上長(zhǎng)方形的孔中,孔長(zhǎng)15 mm,高3 mm,厚2 mm。樣品填入后,用聚酰亞胺膠帶將其兩面密封,在美國(guó)阿貢國(guó)家實(shí)驗(yàn)室先進(jìn)光子源(APS)的20BM線站進(jìn)行分析。室溫下使用熒光模式獲取土壤的XAFS譜,各能量點(diǎn)的數(shù)據(jù)采集時(shí)間為1 s。同時(shí)使用透射模式測(cè)量Fe、Mn及Zn的金屬箔標(biāo)準(zhǔn)樣品并用于能量校準(zhǔn),校準(zhǔn)后Fe、Mn及Zn的X射線吸收近邊結(jié)構(gòu)(XANES)譜的一

    階導(dǎo)數(shù)最大值分別為7 112、6 539及9 659 eV。對(duì)各樣品進(jìn)行3~6次掃描,利用Demeter軟件包對(duì)樣品

    XAFS譜進(jìn)行處理及分析,利用Athena從XAFS譜的邊前區(qū)域扣除由線性函數(shù)定義的背景,對(duì)吸收邊的高度進(jìn)行歸一化得到XANES譜數(shù)據(jù)。將原始XAFS譜扣除一個(gè)平滑的、近似于吸收背景的三次多項(xiàng)式后獲得能量空間的XAFS數(shù)據(jù)(χ(E)),將數(shù)據(jù)從X射線的能量空間轉(zhuǎn)換為(電子)k空間(即波矢空間)得到擴(kuò)展X射線精細(xì)吸收結(jié)構(gòu)(EXAFS)譜χ(k)(即k空間的EXAFS譜數(shù)據(jù))[16]。為放大高k區(qū)域處的信號(hào)振蕩,將χ(k)乘k2獲得加權(quán)擴(kuò)展EXAFS譜。使用Athena對(duì)樣品的XANES譜進(jìn)行線性擬合(linear combination fitting, LCF)。Fe、Mn及Zn的K邊對(duì)應(yīng)范圍分別為7 107~7 157、6 527~6 577和9 644~9 694 eV[11,16]。采集水鐵礦、針鐵礦、黑云母、MnO、MnSO4·4H2O、MnCO3、Mn2O3、γ-MnOOH、Mn3O4、MnO2、Zn2(OH)2CO3、Zn(OH)2、ZnO、ZnSO4·7H2O和ZnCO3的標(biāo)準(zhǔn)譜圖數(shù)據(jù)用于線性擬合。

    利用Artemis對(duì)土壤Fe、Mn及Zn的EXAFS譜進(jìn)行建模,在2×10-10~10×10-10 m-1的范圍內(nèi)應(yīng)用Hanning窗函數(shù)后進(jìn)行傅里葉變換,得到EXAFS譜的徑向分布函數(shù),即R空間數(shù)據(jù)(FT|χ(R)|)(即加權(quán)后的k空間數(shù)據(jù)轉(zhuǎn)換得到徑向分布函數(shù))[16]。使用最小二乘法對(duì)R空間數(shù)據(jù)進(jìn)行基于EXAFS方程的非線性擬合。EXAF譜R空間中的配位數(shù)(CNs)、原子間距(ΔR)、Debye-Waller因子(σ2)為變量。擬合中保持每條路徑的能量位移(ΔE0)相同。通過(guò)擬合優(yōu)度因子評(píng)估模型與實(shí)驗(yàn)數(shù)據(jù)的相似性。

    1.4 數(shù)據(jù)處理及分析方法

    Zn的吸附量表達(dá)式為

    式中:Ce為反應(yīng)平衡時(shí)溶液中溶質(zhì)質(zhì)量濃度(mg/L);qm為吸附介質(zhì)(單分子層)對(duì)溶質(zhì)的飽和吸附量(mg/g);KL為表征吸附結(jié)合能的參數(shù)(L/mg);KF和1/n為Freundlich等溫吸附模型中與吸附能力相關(guān)的參數(shù),一般1/n<1,當(dāng)1/n>2時(shí),吸附反應(yīng)難以進(jìn)行。

    2 結(jié)果與討論

    2.1 Zn在土壤中的吸附動(dòng)力學(xué)

    根據(jù)吸附動(dòng)力學(xué)的實(shí)驗(yàn)結(jié)果(圖1),在Zn初始質(zhì)量濃度為10 mg/L條件下,3層土壤(S1、S2、S3)對(duì)Zn的吸附反應(yīng)動(dòng)力學(xué)過(guò)程存在相當(dāng)多的共性,吸附過(guò)程主要分為快速吸附、慢速吸附及平衡吸附

    階段。反應(yīng)初期(10 min內(nèi)),吸附量快速增加。隨著反應(yīng)時(shí)間的增加,3組樣品對(duì)Zn吸附速率明顯減緩。在8 h左右,樣品對(duì)Zn的吸附趨于平衡,平衡吸附量從大到小依次為S1、S3、S2。

    采用準(zhǔn)一級(jí)和準(zhǔn)二級(jí)動(dòng)力學(xué)模型對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合,擬合結(jié)果顯示,2個(gè)方程均能較好地?cái)M合Zn的吸附動(dòng)力學(xué)數(shù)據(jù)(圖1、表2),其中準(zhǔn)二級(jí)動(dòng)力學(xué)方程擬合的理論吸附量與實(shí)驗(yàn)結(jié)果更接近(擬合優(yōu)度>0.99),說(shuō)明化學(xué)吸附可能是剖面土壤吸附Zn的速率控制步驟[21]。

    2.2 Zn在土壤中的吸附熱力學(xué)

    3層土壤吸附Zn的等溫吸附曲線如圖2所示。隨著背景液中Zn質(zhì)量濃度增加,各剖面土壤對(duì)Zn的吸附量均呈現(xiàn)出先增加后趨于平緩的趨勢(shì),這與吉林地區(qū)剖面土壤對(duì)Zn的吸附熱力學(xué)特征一致[22]。這可能歸因于在吸附初始階段,低濃度Zn離子優(yōu)先與土壤表面高吸附位點(diǎn)結(jié)合;隨著溶液質(zhì)量濃度的增加,高吸附點(diǎn)位逐漸飽和,此時(shí)土壤中對(duì)Zn吸附的區(qū)域由專性吸附轉(zhuǎn)變?yōu)榉菍P晕綖橹鳎剿俾首兟敝吝_(dá)到動(dòng)態(tài)平衡[20,22]。

    采用Langmuir和Freundlich模型對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合,Langmuir方程常用于描述在相同吸附位點(diǎn)上的有限數(shù)量表面的單層吸附,而Freundlich方程通常描述異質(zhì)固體表面上的吸附[23]。由表3可知,兩個(gè)模型均能較好地模擬3層土壤對(duì)Zn的吸附行為(擬合優(yōu)度>0.9),且Langmuir擬合效果更好,這與Bauer等的研究結(jié)果一致[24],模型的非線性擬合結(jié)果也印證了上述實(shí)驗(yàn)現(xiàn)象(圖2)。Freundlich等溫吸附系數(shù)KF越大,吸附介質(zhì)對(duì)溶質(zhì)的吸附能力越強(qiáng)。根據(jù)Freundlich模型擬合結(jié)果中KF及Langmuir模型擬合的飽和吸附量從大到小排序均為S1、S3、S2,實(shí)驗(yàn)觀測(cè)和模型模擬結(jié)果一致,表明3層土壤對(duì)Zn的吸附能力差異較小,從大到小排序?yàn)镾1、S3、S2。

    研究[25]表明土壤pH、CEC及有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)越高,其表面負(fù)電荷越多,吸附重金屬Zn2+的能力也越強(qiáng)。中層土S2對(duì)Zn吸附量較低可能與其較低的pH及有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)有關(guān),對(duì)Zn的吸附容量小于其他樣品。表層雜填土S1的pH、CEC及有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)均高于S3(表1),這可能是土壤S1對(duì)Zn的吸附能力相對(duì)較強(qiáng)的原因。

    2.3 土壤XRD及Fe、Mn的同步輻射

    2.3.1 土壤XRD

    由圖3可知,同一剖面土壤(S1、S2、S3)的礦物相組成較為相似,主要包括石英、長(zhǎng)石和黏土礦物(高嶺石)等,每層土壤礦物種類及質(zhì)量分?jǐn)?shù)略有差異,較大顆粒物為石英,高嶺石等黏土礦物為較細(xì)顆粒,在較深的土壤中質(zhì)量分?jǐn)?shù)較高。

    重金屬在土壤中的溶解度、遷移性和毒性取決于土壤中與其結(jié)合的固相組分,如黏土礦物、鐵錳氧化物及有機(jī)物等[26]。黏土礦物可為重金屬提供活性吸附位點(diǎn)[24],土壤中相當(dāng)一部分Zn賦存于黏土礦物中,特別是高嶺石族礦物[27]。一方面,剖面中層(S2)及深層(S3)土壤Zn質(zhì)量分?jǐn)?shù)明顯高于表層(S1)土壤(表1),這可能歸因于下層黏土礦物對(duì)Zn的吸附;另一方面,靜態(tài)吸附實(shí)驗(yàn)結(jié)果表明,黏土礦物質(zhì)量分?jǐn)?shù)最少的表層土壤對(duì)Zn的吸附能力最強(qiáng),進(jìn)一步說(shuō)明黏土礦物對(duì)于土壤吸附Zn能力的影響程度有限,并不是首要影響因素。

    2.3.2 土壤Fe、Mn的同步輻射

    Fe XANES譜的最大峰值位于7 133 eV,表明3層土壤樣品(S1、S2、S3)中的Fe主要以Fe3+形式存在(圖4a)。Fe XANES譜的LCF結(jié)果(表4)表明,土壤中的Fe主要以Fe3+的氫氧化物(水鐵礦和針鐵礦)(體積分?jǐn)?shù)為82%~93%)形式存在,此外還有少量黑云母(7%~13%)。

    傅里葉變換的EXAFS譜(圖4b、c)表明,土壤具有與水鐵礦相似的原子近鄰環(huán)境,即共邊Fe-O八面體配位結(jié)構(gòu)(Fe-Fe間距為3×10-10 m)和共角Fe-O八面體配位結(jié)構(gòu)(Fe-Fe間距為3.4×10-10 m),上述特征與水鐵礦及其可能的相變礦物針鐵礦特征對(duì)應(yīng)。從圖5可以看出,表層土S1和深層土S3的Mn XANES譜的最大峰值位于6 558 eV,而S2的Mn XANES譜則出現(xiàn)雙峰結(jié)構(gòu),峰值位置分別位于6 552和6 556 eV;這表明Mn在土壤可能以混合價(jià)態(tài)錳氧化物形式存在,土壤S1和S3中的含Mn礦物主要為Mn2O3和Mn3O4(體積分?jǐn)?shù)73.1%~74.2%),還有少量的MnO2(18.5%~22.2%),此外,中層土S2中的Mn可能以Mn2+和Mn3+混合形式存在。上述結(jié)果表明,元素Fe及Mn在剖面不同深度土壤中的存在形式基本一致,主要以鐵的氫氧化物(水鐵礦和針鐵礦)及混合價(jià)態(tài)錳氧化物形式存在。結(jié)合前文,剖面土壤對(duì)Zn吸附能力與土壤中Fe、Mn質(zhì)量分?jǐn)?shù)呈正相關(guān)(表1、表3),兩者隨深度變化趨勢(shì)一致(土壤吸附Zn能力及Fe、Mn質(zhì)量分?jǐn)?shù)從大到小為S1、S3、S2)。此外,三層土壤中S2吸附Zn的能力最低,可能與其中Mn XANES譜中獨(dú)特的雙峰結(jié)構(gòu)及Mn的賦存形態(tài)特征有關(guān)。以上分析表明,土壤pH較低時(shí),相比于黏土礦物,以水鐵礦、針鐵礦及混合價(jià)態(tài)錳氧化物為代表的鐵錳氧化物對(duì)土壤吸附重金屬Zn的影響更加顯著。其中,相對(duì)于鐵氧化物而言,錳氧化物質(zhì)量分?jǐn)?shù)較低,但由于其比表面積大、表面電荷低、對(duì)重金屬離子親和力強(qiáng)等特性,能夠在一定程度上深刻影響土壤對(duì)于重金屬Zn的吸附能力,這與以往研究結(jié)論一致[29]。

    Zn2+可以通過(guò)多種作用力吸附于鐵錳氧化物表面,依靠物理吸附、靜電引力以及離子交換吸附所形成的產(chǎn)物通常被稱為外層吸附產(chǎn)物,依靠化學(xué)吸附所形成的產(chǎn)物一般被稱為內(nèi)層吸附產(chǎn)物。吸附類型直接影響金屬離子在礦物表面的微觀構(gòu)型、活性、吸附速率及可逆性[30],與重金屬在環(huán)境中的可遷移性及毒性密切相關(guān)。研究[31]表明,Zn2+吸附于鐵和錳的非晶態(tài)水合氧化物上時(shí)保留了其初級(jí)水合層,形成外層絡(luò)合物,與礦物表面的結(jié)合力較弱,在宏觀上表現(xiàn)出較高的吸附可逆性。一方面,水鐵礦具有極大的比表面積(120~850 m2/g)和高表面活性[32],Zn2+吸附于水鐵礦等亞穩(wěn)態(tài)鐵氧化物與針鐵礦等晶態(tài)鐵氧化物的過(guò)程中可完全或部分失去其水合層,形成內(nèi)層絡(luò)合物,與礦物表面的結(jié)合力較強(qiáng)。另一方面,水鐵礦由于其熱力學(xué)不穩(wěn)定性成為晶質(zhì)鐵礦物的前驅(qū)體,在向針鐵礦等晶質(zhì)氧化物轉(zhuǎn)變的過(guò)程中,其表面積和對(duì)重金屬的螯合能力下降,其中保留的重金屬組分將被再度活化或重新分配[19],從而增加土壤中重金屬向水相遷移的風(fēng)險(xiǎn)。

    2.4 吸附反應(yīng)前后Zn在土壤中的賦存形態(tài)及遷移性

    吸附反應(yīng)前土壤(S1、S2、S3)及反應(yīng)后土壤([XC5號(hào).TIF]S1、[XC5號(hào).TIF]S2及[XC5號(hào).TIF]S3)Zn XANES譜的最大峰值均位于9 668 eV,表明反應(yīng)前土壤中Zn以Zn2+形態(tài)存在(圖6a)。土壤XANES譜的定性分析結(jié)果表明,反應(yīng)前后土壤中Zn局部配位結(jié)構(gòu)及形態(tài)組成發(fā)生了一定變化(圖6b、c,圖7)。土壤與標(biāo)準(zhǔn)物質(zhì)XANES譜圖特征的相似性表明反應(yīng)前后土壤中Zn的第一殼層配位模式分別與Zn2(OH)2CO3及ZnSO4·7H2O相近。XANES譜的LCF結(jié)果(表5)進(jìn)一步對(duì)反應(yīng)前后土壤與標(biāo)準(zhǔn)物質(zhì)譜圖之間的相似性進(jìn)行定量描述,結(jié)果表明,反應(yīng)前土壤譜與Zn2(OH)2CO3譜的失配小于2%(擬合因子為0.012~0.015),反應(yīng)后土壤中Zn表現(xiàn)出類似于ZnSO4·7H2O(體積分?jǐn)?shù)47.4%~59.6%)和Zn2(OH)2CO3(36.4%~46.8%)的局部配位環(huán)境。以上分析結(jié)果表明,反應(yīng)前后土壤中Zn可能由四配位的Zn-O四面體配位轉(zhuǎn)變?yōu)橐粤湮坏腪n-O八面體和四配位的Zn-O四面體混合形式存在,反應(yīng)后土壤中Zn的配位模式與以往研究中Zn2+吸附于全國(guó)不同地區(qū)土壤所觀察到的結(jié)果相似[16]。

    傅里葉變換的EXAFS譜(圖6c、圖7c)表明,反應(yīng)前后土壤的第一殼層在約2.1×10-10 m處有對(duì)應(yīng)于近鄰O原子的散射信號(hào),而在3×10-10~4×10-10 m處存在雜亂的信號(hào),考慮到土壤中Zn的質(zhì)量分?jǐn)?shù)較低,該處信號(hào)可能是低信噪比信號(hào)帶來(lái)的“假峰”。Zn EXAFS譜建模結(jié)果表明,反應(yīng)前后Zn-O原子對(duì)間距分別為2.05×10-10~2.07×10-10 m及2.06×10-10~2.08×10-10 m,配位數(shù)分別為3.4~4.3及4.3~5.6,分別對(duì)應(yīng)于前文提及的Zn-O四面體配位及Zn-O四面體和八面體混合配位模式。LCF分析結(jié)果表明,原始土壤中的Zn主要以堿式碳酸鋅類礦物(如水鋅礦)的形式存在,而反應(yīng)后土壤中Zn可能以硫酸鹽及堿式碳酸鹽的混合形式存在(表5)。LCF結(jié)果用于物相分析時(shí)應(yīng)當(dāng)結(jié)合標(biāo)準(zhǔn)物質(zhì)的物理化學(xué)特性,ZnSO4·7H2O在水中溶解極高(960 g/L),結(jié)晶ZnSO4·7H2O的EXAFS信號(hào)與水溶態(tài)Zn的光譜相似度高[33],酸性環(huán)境下水溶態(tài)Zn2+以六配位的水合離子八面體結(jié)構(gòu)存在,Zn-O原子間距為(2.08+0.02)×10-10 m[34]。反應(yīng)后土壤中Zn與水溶態(tài)Zn2+局部配位結(jié)構(gòu)的相似性為形成外層絡(luò)合物提供了證據(jù)[1],這與土壤吸附Zn熱力學(xué)模型擬合結(jié)果略有不同,這可能歸因于土壤作為一個(gè)多介質(zhì)多界面的復(fù)合體系,具有復(fù)雜的吸附活性中心[35],其對(duì)Zn2+的吸附與物理吸附和化學(xué)吸附有關(guān),是二者共同作用的結(jié)果。LCF結(jié)果中ZnCO3的存在表明反應(yīng)過(guò)程中可能生成了沉淀。此外,不能排除Zn與土壤中的有機(jī)質(zhì)結(jié)合的可能性。

    土壤中Zn的遷移、轉(zhuǎn)化及生物有效性與其在土壤孔隙水中的化學(xué)形態(tài)密切相關(guān)。一般而言,重金屬的可交換態(tài)、碳酸鹽結(jié)合態(tài)和鐵錳氧化物結(jié)合態(tài)為有效態(tài),其可遷移性及生物有效性較高[36]。有效態(tài)組分對(duì)pH值變化敏感,pH降低會(huì)促進(jìn)有效態(tài)Zn向水體中釋放,酸性土壤中的Zn具有較高的遷移率[37]。受歷史礦山開(kāi)采等活動(dòng)影響,剖面土壤呈酸性,pH值為4.46~4.78(表1),場(chǎng)地土壤吸附重金屬能力下降,土壤中Zn的遷移能力提高,從而增加地下水污染風(fēng)險(xiǎn)。此外,強(qiáng)降雨條件會(huì)加速土壤侵蝕,提高土壤中Zn的遷移能力。

    贛南地區(qū)位于亞熱帶季風(fēng)區(qū),雨量充沛,酸雨出現(xiàn)頻率較高,研究地塊表層土壤Zn質(zhì)量分?jǐn)?shù)較低可能與土壤pH值和Zn的賦存形態(tài)及降雨淋濾有關(guān)。歷史礦山開(kāi)采等活動(dòng)導(dǎo)致外源Zn進(jìn)入土壤時(shí),土壤對(duì)Zn表現(xiàn)出一定的滯留能力。相比于原始土壤,吸附反應(yīng)后Zn以類似于硫酸鹽的外層絡(luò)合物或碳酸鋅沉淀形式存在時(shí),其在土壤中的可遷移性降低。反應(yīng)前后Zn部分以堿式碳酸鹽形態(tài)存在時(shí),在酸性環(huán)境下易于溶解,隨降雨沖刷向深層土壤遷移,導(dǎo)致深層土壤Zn質(zhì)量分?jǐn)?shù)增加,增強(qiáng)地下水Zn污染風(fēng)險(xiǎn)。

    3 結(jié)論

    1)贛南廢棄礦區(qū)同一剖面表層雜填土、中層風(fēng)化巖及深層風(fēng)化巖3層土壤對(duì)Zn的吸附曲線存在相當(dāng)多的共性,土壤吸附Zn的平衡時(shí)間約為8 h,剖面土壤對(duì)Zn的吸附能力由

    大到小依次排序?yàn)楸韺与s填土、深層風(fēng)化巖、中層風(fēng)化巖,準(zhǔn)二級(jí)動(dòng)力學(xué)模型和Langmuir吸附模型較好地描述了土壤吸附Zn的過(guò)程。

    2)剖面不同深度土壤主要物相組成包含石英、長(zhǎng)石及黏土礦物等,元素Fe、Mn在土壤中主要以鐵的氫氧化物(水鐵礦及針鐵礦)及混合價(jià)態(tài)錳氧化物形式存在。剖面各層土壤理化性質(zhì)差異影響土壤對(duì)Zn的吸附行為,Zn在不同深度土壤中吸附能力大小與土壤中Fe、Mn質(zhì)量分?jǐn)?shù)呈正相關(guān),相比于黏土礦物,鐵錳氧化物對(duì)土壤吸附Zn能力影響更加顯著。

    3)吸附反應(yīng)前后土壤中Zn的局部配位結(jié)構(gòu)及形態(tài)組成發(fā)生了變化,局部配位結(jié)構(gòu)由Zn-O四面體配位轉(zhuǎn)變?yōu)閆n-O四面體和八面體的混合配位模式,反應(yīng)后土壤表面可能形成了類似硫酸鹽的外層絡(luò)合物,物理吸附及化學(xué)吸附共同參與了剖面土壤吸附Zn的過(guò)程。受土壤pH影響,土壤吸附滯留Zn能力有限,土壤中的有效態(tài)Zn易伴隨降雨沖刷向下層土壤遷移,土壤中Zn質(zhì)量分?jǐn)?shù)隨深度增加而增加,增強(qiáng)地下水污染風(fēng)險(xiǎn)。

    參考文獻(xiàn)(References):

    [1] Wang Z, Kwon K D, Peacock C, et al. Zn Stable Isotope Fractionation During Adsorption onto Todorokite: A Molecular Perspective from X-Ray Absorption Spectroscopy and Density Functional Theory[J]. Geochimica et Cosmochimica Acta, 2022, 327: 116-136.

    [2] 郭鐘群, 趙奎, 金解放, 等. 離子型稀土礦環(huán)境風(fēng)險(xiǎn)評(píng)估及污染治理研究進(jìn)展[J]. 稀土, 2019, 40(3): 115-126.

    Guo Zhongqun, Zhao Kui, Jin Jiefang, et al. Reviews on Environmental Assessment and Pollution Prevention of Ion Adsorption Type Rare Earth Ores[J]. Chinese Rare Earths, 2019, 40(3): 115-126.

    [3] Qiao J, Tang J, Xue Q. Study on Pb Release by Several New Lixiviants in Weathered Crust Elution-Deposited Rare Earth Ore Leaching Process: Behavior and Mechanism[J]. Ecotoxicology and Environmental Safety, 2020, 190: 110138.

    [4] Zhang Q, Ren F, Li F, et al. Ammonia Nitrogen Sources and Pollution Along Soil Profiles in an in-Situ Leaching Rare Earth Ore[J]. Environmental Pollution, 2020, 267(12): 115449.

    [5] 李曉艷, 張青偉, 洪松濤, 等. 不同氧化還原條件下鉛鋅礦尾砂中重金屬元素活化遷移規(guī)律[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2022, 52(2): 434-441.

    Li Xiaoyan, Zhang Qingwei, Hong Songtao, et al.Activation and Migration of Heavy Metal Elements in Lead Zinc Ore Tailings Under Different Redox Conditions[J]. Journal of Jilin University (Earth Science Edition), 2022, 52(2): 434-441.

    [6] 李英鈞, 宋澤峰, 楊悅鎖, 等. 石棉尾礦場(chǎng)地土壤污染地球化學(xué)特征及其對(duì)微生物多樣性的影響[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2024,54(3):980-992.

    Li Yingjun, Song Zefeng, Yang Yuesuo, et al. Effects of Asbestos Tailings Site on Geochemistry and Microbial Diversity of Contaminated Soils[J]. Journal of Jilin University (Earth Science Edition),2024,54(3):980-992.

    [7] Sparks D L, Singh B, Siebecker M G. Environmental Soil Chemistry[M]. Elsevier:[s.n.], 2022.

    [8] Zeng J, Tabelin C B, Gao W, et al. Heterogeneous Distributions of Heavy Metals in the Soil-Groundwater System Empowers the Knowledge of the Pollution Migration at a Smelting Site[J]. Chemical Engineering Journal, 2023, 454(2): 140307.

    [9] Dewey C, Juillot F, Fendorf S, et al. Seasonal Oxygenation of Contaminated Floodplain Soil Releases Zn to Porewater[J]. Environmental Science amp; Technology, 2023, 57(12): 4841-4851.

    [10] Huang B, Yuan Z, Li D, et al. Effects of Soil Particle Size on the Adsorption, Distribution, and Migration Behaviors of Heavy Metal(loid)s in Soil: A Review[J]. Environmental Science: Processes amp; Impacts, 2020, 22(8): 1596-1615.

    [11] Lopes G, Li W, Siebecker M G, et al. Combining Zinc Desorption with EXAFS Speciation Analysis to Understand Zn Mobility in Mining and Smelting Affected Soils in Minas Gerais, Brazil[J]. Science of the Total Environment, 2021, 754: 142450.

    [12] Roberts D R, Scheinost A C, Sparks D L. Zinc Speciation in a Smelter-Contaminated Soil Profile Using Bulk and Microspectroscopic Techniques[J]. Environmental Science amp; Technology, 2002, 36(8): 1742-1750.

    [13] Yu H, Li C, Yan J, et al. A Review on Adsorption Characteristics and Influencing Mechanism of Heavy Metals in Farmland Soil[J]. RSC Advances, 2023, 13(6): 3505-3519.

    [14] Li Q, Wang Y, Li Y, et al. Speciation of Heavy Metals in Soils and Their Immobilization at Micro-Scale Interfaces Among Diverse Soil Components[J]. Science of the Total Environment, 2022, 825: 153862.

    [15] Mo X, Siebecker M G, Gou W, et al. A Review of Cadmium Sorption Mechanisms on Soil Mineral Surfaces Revealed from Synchrotron-Based X-Ray Absorption Fine Structure Spectroscopy: Implications for Soil Remediation[J]. Pedosphere, 2021, 31(1): 11-27.

    [16] Gu J, Cui P, Zeng Z, et al. Soil Sorption as Influenced by Soil Properties and the Zn-O Bond Distance: Insights from EXAFS and MC-ICP-MS[J]. Applied Geochemistry, 2022, 145: 105427.

    [17] Fan T T, Cui P X, Sun Q, et al. Combining Path Analysis and X-Ray Absorption Spectroscopy to Unravel the Zn Sorption Mechanism on Soils[J]. Soil Science Society of America Journal, 2018, 82(4): 796-802.

    [18] Sipos P, Németh T, Kis V K, et al. Sorption of Copper, Zinc and Lead on Soil Mineral Phases[J]. Chemosphere, 2008, 73(4): 461-469.

    [19] Shi M, Min X, Ke Y, et al. Recent Progress in Understanding the Mechanism of Heavy Metals Retention by Iron (Oxyhydr) Oxides[J]. Science of the Total Environment, 2021, 752: 141930.

    [20] 劉凌青, 肖細(xì)元, 郭朝暉, 等. 鋅冶煉地塊剖面土壤對(duì)鎘、鉛的吸附特征及機(jī)制[J]. 環(huán)境科學(xué), 2021, 42(8): 4015-4023.

    Liu Lingqing, Xiao Xiyuan, Guo Zhaohui, et al. Adsorption Characteristics and Mechanism of Cd and Pb in Tiered Soil Profiles from a Zinc Smelting Site[J]. Environmental Science, 2021, 42(8): 4015-4023.

    [21] Zhang J, Wang X, Zhu Y, et al. The Influence of Heavy Metals in Road Dust on the Surface Runoff Quality: Kinetic, Isotherm, and Sequential Extraction Investigations[J]. Ecotoxicology and Environmental Safety, 2019, 176: 270-278.

    [22] 姜延, 馬秀蘭, 郭亞馨, 等. 不同分層土壤對(duì)鋅的吸附解吸特性[J]. 東北林業(yè)大學(xué)學(xué)報(bào), 2021, 49(3): 99-107.

    Jiang Yan, Ma Xiulan, Guo Yaxin, et al.Adsorption and Desorption Characteristics of Zinc in Different Layered Soils[J]. Journal of Northeast Forestry University, 2021, 49( 3) : 99-107.

    [23] Chen X, Hossain M F, Duan C, et al. Isotherm Models for Adsorption of Heavy Metals from Water: A Review[J]. Chemosphere, 2022: 135545.

    [24] Bauer T V, Pinskii D L, Minkina T M, et al. Application of XAFS and XRD Methods for Describing the Copper and Zinc Adsorption Characteristics in Hydromorphic Soils[J]. Environmental Geochemistry and Health, 2022, 44(2): 335-347.

    [25] Qin L, Wang M, Zhao S, et al. Effect of Soil Leaching on the Toxicity Thresholds (ECx) of Zn in Soils with Different Properties[J]. Ecotoxicology and Environmental Safety, 2021, 228: 112999.

    [26] Pan Y, Fu Y, Liu S, et al. Spatial and Temporal Variations of Metal Fractions in Paddy Soil Flooding with Acid Mine Drainage[J]. Environmental Research, 2022, 212: 113241.

    [27] Tian L, Fu K bin, Chen S, et al. Comparison of Microscopic Adsorption Characteristics of Zn(Ⅱ), Pb(Ⅱ), and Cu(Ⅱ) on Kaolinite[J]. Scientific Reports, 2022, 12(1): 15936.

    [28] Minkina T, Nevidomskaya D, Bauer T, et al. Determining the Speciation of Zn in Soils Around the Sediment Ponds of Chemical Plants by XRD and XAFS Spectroscopy and Sequential Extraction[J]. Science of the Total Environment, 2018, 634: 1165-1173.

    [29] Covelo E F, lvarez N, Andrade Couce M L, et al. Zn Adsorption by Different Fractions of Galician Soils[J]. Journal of Colloid and Interface Science, 2004, 280(2): 343-349.

    [30] Pan G, Qin Y, Li X, et al. EXAFS Studies on Adsorption-Desorption Reversibility at Manganese Oxides-Water Interfaces: I: Irreversible Adsorption of Zinc onto Manganite (γ-MnOOH)[J]. Journal of Colloid and Interface Science, 2004, 271(1): 28-34.

    [31] Trivedi P, Dyer J A, Sparks D L, et al. Mechanistic and Thermodynamic Interpretations of Zinc Sorption onto Ferrihydrite[J]. Journal of Colloid and Interface Science, 2004, 270(1): 77-85.

    [32] Li B, Wei D, Zhou Y, et al. Mechanisms of Arsenate and Cadmium Co-Immobilized on Ferrihydrite Inferred from Ternary Surface Configuration[J]. Chemical Engineering Journal, 2021, 424: 130410.

    [33] Kuzmin A, Obst S, Purans J. X-Ray Absorption Spectroscopy and Molecular Dynamics Studies of Hydration in Aqueous Solutions[J]. Journal of Physics: Condensed Matter, 1997, 9(46): 10065.

    [34] Ha J, Trainor T P, Farges F, et al. Interaction of

    Aqueous Zn(Ⅱ) with Hematite Nanoparticles and Microparticles. Part 1: EXAFS Study of Zn(Ⅱ) Adsorption and Precipitation[J]. Langmuir, 2009, 25(10): 5574-5585.

    [35] 胡世文, 劉同旭, 李芳柏, 等. 土壤鐵礦物的生物-非生物轉(zhuǎn)化過(guò)程及其界面重金屬反應(yīng)機(jī)制的研究進(jìn)展[J].土壤學(xué)報(bào), 2022, 59(1): 54-65.

    Hu Shiwen, Liu Tongxu, Li Fangbai, et al.The Abiotic and Biotic Transformation Processes of Soil Iron-Bearing Minerals and Its Interfacial Reaction Mechanisms of Heavy Metals: A Review [J]. Acta Pedologica Sinica,2022, 59(1):54-65.

    [36] Liu Y, Xiao T, Zhu Z, et al. Geogenic Pollution, Fractionation and Potential Risks of Cd and Zn in Soils from a Mountainous Region Underlain by Black Shale[J]. Science of the Total Environment, 2021, 760: 143426.

    [37] Fernández-Calvio D, Pateiro-Moure M, Nóvoa-Muoz J C, et al. Zinc Distribution and Acid–Base Mobilisation in Vineyard Soils and Sediments[J]. Science of the Total Environment, 2012, 414: 470-479.

    猜你喜歡
    吸附遷移土壤
    土壤
    靈感的土壤
    識(shí)破那些優(yōu)美“擺拍”——鏟除“四風(fēng)”的土壤
    靈感的土壤
    Fe(Ⅲ)負(fù)載改性核桃殼對(duì)Cu2+吸附研究
    炭氣凝膠的制備與吸附氫氣性能
    核桃殼基吸附材料的制備及性能測(cè)試
    科技視界(2016年25期)2016-11-25 12:30:06
    淺析遷移規(guī)律在足球教學(xué)中的影響
    資治文摘(2016年7期)2016-11-23 01:50:57
    現(xiàn)代人類遷移與信息、文化流動(dòng)的雙向不對(duì)稱現(xiàn)象
    活化粉煤灰在煤焦油加工廢水預(yù)處理中的應(yīng)用
    一级作爱视频免费观看| 日日夜夜操网爽| 欧美 亚洲 国产 日韩一| 午夜两性在线视频| 777久久人妻少妇嫩草av网站| 久久精品成人免费网站| 国产主播在线观看一区二区| 国产99久久九九免费精品| 精品久久久久久久久久免费视频 | 正在播放国产对白刺激| 啦啦啦在线免费观看视频4| 超碰成人久久| 亚洲自偷自拍图片 自拍| 国产精品1区2区在线观看.| 精品电影一区二区在线| 久久久久久人人人人人| 三上悠亚av全集在线观看| 亚洲欧美日韩无卡精品| 一进一出好大好爽视频| 亚洲狠狠婷婷综合久久图片| 人成视频在线观看免费观看| 视频在线观看一区二区三区| 国产精品秋霞免费鲁丝片| 中文字幕色久视频| 97超级碰碰碰精品色视频在线观看| 午夜老司机福利片| 欧美成狂野欧美在线观看| 香蕉国产在线看| 国产高清国产精品国产三级| 亚洲国产欧美网| 精品卡一卡二卡四卡免费| 看免费av毛片| 男男h啪啪无遮挡| 亚洲中文日韩欧美视频| 老汉色∧v一级毛片| 黄片小视频在线播放| 免费搜索国产男女视频| av天堂久久9| 国产精品1区2区在线观看.| 一a级毛片在线观看| 51午夜福利影视在线观看| 国产区一区二久久| 久久精品国产综合久久久| 18禁观看日本| 免费观看人在逋| 色婷婷av一区二区三区视频| 欧美国产精品va在线观看不卡| 亚洲欧美日韩无卡精品| av超薄肉色丝袜交足视频| 国产亚洲欧美精品永久| 亚洲中文字幕日韩| 国产三级黄色录像| 国产97色在线日韩免费| 又黄又爽又免费观看的视频| 国产不卡一卡二| 精品免费久久久久久久清纯| 日韩欧美在线二视频| 亚洲专区字幕在线| 国产一区在线观看成人免费| a级毛片在线看网站| 午夜久久久在线观看| 欧美日韩瑟瑟在线播放| 欧美午夜高清在线| 久久精品亚洲av国产电影网| 国产在线精品亚洲第一网站| 天天躁夜夜躁狠狠躁躁| 午夜亚洲福利在线播放| 中文字幕色久视频| 亚洲熟女毛片儿| 免费在线观看完整版高清| 亚洲人成电影观看| 欧美最黄视频在线播放免费 | 久久中文字幕一级| 夫妻午夜视频| 久久国产精品影院| 成年人免费黄色播放视频| 国产精品美女特级片免费视频播放器 | 午夜精品久久久久久毛片777| 我的亚洲天堂| 成人av一区二区三区在线看| 成人亚洲精品一区在线观看| 婷婷丁香在线五月| 午夜激情av网站| 80岁老熟妇乱子伦牲交| 久久香蕉激情| 国产伦一二天堂av在线观看| 日韩欧美一区二区三区在线观看| 99精品久久久久人妻精品| 久久国产亚洲av麻豆专区| 日韩欧美国产一区二区入口| 最近最新中文字幕大全免费视频| 在线播放国产精品三级| 99久久久亚洲精品蜜臀av| 欧美激情极品国产一区二区三区| 亚洲欧美精品综合久久99| 手机成人av网站| 另类亚洲欧美激情| 国产欧美日韩一区二区精品| 在线观看免费高清a一片| av中文乱码字幕在线| 久久人妻熟女aⅴ| 精品日产1卡2卡| 亚洲色图av天堂| 国产av一区二区精品久久| 日韩大码丰满熟妇| 精品卡一卡二卡四卡免费| 成人三级黄色视频| 亚洲五月婷婷丁香| 9191精品国产免费久久| 日本 av在线| 久久精品成人免费网站| 热re99久久国产66热| 久久久国产欧美日韩av| 一边摸一边抽搐一进一小说| 国产成人免费无遮挡视频| 乱人伦中国视频| 国产有黄有色有爽视频| 亚洲男人的天堂狠狠| 国产视频一区二区在线看| 十八禁网站免费在线| 操出白浆在线播放| 成人永久免费在线观看视频| 国产一卡二卡三卡精品| 免费在线观看亚洲国产| 变态另类成人亚洲欧美熟女 | 香蕉国产在线看| 99国产精品99久久久久| 男女做爰动态图高潮gif福利片 | 在线av久久热| 亚洲午夜精品一区,二区,三区| 国产免费av片在线观看野外av| 两性夫妻黄色片| 在线观看免费高清a一片| 日韩成人在线观看一区二区三区| 这个男人来自地球电影免费观看| 中文字幕av电影在线播放| av天堂在线播放| 日本a在线网址| 男女高潮啪啪啪动态图| 人人妻,人人澡人人爽秒播| 黄频高清免费视频| 精品一区二区三区av网在线观看| 久久久国产欧美日韩av| 制服诱惑二区| 欧美 亚洲 国产 日韩一| 国产精品久久久人人做人人爽| 国产av又大| 精品国产国语对白av| 国产不卡一卡二| 成人手机av| 久久久久久久久中文| 精品国产一区二区三区四区第35| 久久久国产成人免费| 亚洲av片天天在线观看| 国产欧美日韩一区二区三| 久99久视频精品免费| 亚洲va日本ⅴa欧美va伊人久久| 在线永久观看黄色视频| 两人在一起打扑克的视频| 久久久久国产精品人妻aⅴ院| 国产91精品成人一区二区三区| 精品熟女少妇八av免费久了| 亚洲国产精品合色在线| 亚洲中文av在线| 久久精品国产99精品国产亚洲性色 | 热re99久久国产66热| 一级毛片精品| 神马国产精品三级电影在线观看 | 亚洲欧美激情在线| 亚洲精品美女久久久久99蜜臀| aaaaa片日本免费| 午夜激情av网站| 99在线人妻在线中文字幕| 欧美中文日本在线观看视频| 国产一区二区三区视频了| 99国产综合亚洲精品| 午夜成年电影在线免费观看| 久久人人97超碰香蕉20202| 18禁美女被吸乳视频| 丰满迷人的少妇在线观看| 亚洲av成人不卡在线观看播放网| 大香蕉久久成人网| 国产1区2区3区精品| 天天影视国产精品| 波多野结衣一区麻豆| 午夜福利影视在线免费观看| 99国产精品免费福利视频| 精品久久久久久电影网| 丰满饥渴人妻一区二区三| 国产精品爽爽va在线观看网站 | 亚洲国产欧美一区二区综合| 黑人操中国人逼视频| 久久精品aⅴ一区二区三区四区| 国产熟女午夜一区二区三区| 亚洲专区字幕在线| 欧美激情极品国产一区二区三区| 成人三级做爰电影| 嫩草影视91久久| 欧美在线一区亚洲| 亚洲av片天天在线观看| 亚洲成av片中文字幕在线观看| 国产亚洲欧美98| 少妇 在线观看| 色综合站精品国产| 人人澡人人妻人| 村上凉子中文字幕在线| 久久香蕉激情| 中文亚洲av片在线观看爽| 久久精品人人爽人人爽视色| 国产高清videossex| 精品电影一区二区在线| 这个男人来自地球电影免费观看| 色婷婷久久久亚洲欧美| 久久久久久大精品| 级片在线观看| 男女午夜视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 一二三四在线观看免费中文在| 色综合站精品国产| 我的亚洲天堂| 亚洲视频免费观看视频| 操美女的视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 免费看十八禁软件| 久久精品亚洲熟妇少妇任你| 亚洲av日韩精品久久久久久密| 亚洲avbb在线观看| 亚洲成人免费av在线播放| 久久性视频一级片| 亚洲专区国产一区二区| 亚洲av电影在线进入| 在线播放国产精品三级| 国产真人三级小视频在线观看| 99精品久久久久人妻精品| 97人妻天天添夜夜摸| 亚洲av电影在线进入| 制服人妻中文乱码| 久久亚洲精品不卡| 女人被狂操c到高潮| 999久久久精品免费观看国产| 少妇裸体淫交视频免费看高清 | 一进一出好大好爽视频| 校园春色视频在线观看| 亚洲欧美日韩无卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 成人特级黄色片久久久久久久| 男女之事视频高清在线观看| 在线观看舔阴道视频| 日韩欧美三级三区| 高清在线国产一区| 午夜福利欧美成人| 桃色一区二区三区在线观看| 午夜精品国产一区二区电影| 黄片大片在线免费观看| 激情在线观看视频在线高清| 天堂中文最新版在线下载| 免费在线观看日本一区| 亚洲成人国产一区在线观看| 色尼玛亚洲综合影院| 婷婷精品国产亚洲av在线| 国产精华一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 大陆偷拍与自拍| 国产成人av激情在线播放| 亚洲片人在线观看| www国产在线视频色| 少妇粗大呻吟视频| 中文欧美无线码| 曰老女人黄片| 高清黄色对白视频在线免费看| 久久国产亚洲av麻豆专区| 亚洲精品一二三| 精品日产1卡2卡| 久久久久亚洲av毛片大全| 国产精华一区二区三区| 国产亚洲精品久久久久久毛片| 亚洲国产精品一区二区三区在线| 高清欧美精品videossex| 日韩三级视频一区二区三区| 精品一区二区三区av网在线观看| 18禁国产床啪视频网站| 男女高潮啪啪啪动态图| 久热这里只有精品99| 首页视频小说图片口味搜索| 涩涩av久久男人的天堂| 精品国产一区二区久久| 12—13女人毛片做爰片一| 亚洲精华国产精华精| 一本大道久久a久久精品| 成人免费观看视频高清| 亚洲 欧美一区二区三区| 91麻豆精品激情在线观看国产 | 人人妻人人爽人人添夜夜欢视频| 国产aⅴ精品一区二区三区波| 一a级毛片在线观看| 国产97色在线日韩免费| 999久久久国产精品视频| 在线观看免费日韩欧美大片| 亚洲熟女毛片儿| a级毛片在线看网站| 不卡av一区二区三区| 国产成人精品在线电影| 热re99久久国产66热| 天堂√8在线中文| 国产人伦9x9x在线观看| 亚洲视频免费观看视频| 法律面前人人平等表现在哪些方面| 国产男靠女视频免费网站| 99在线视频只有这里精品首页| 变态另类成人亚洲欧美熟女 | 自拍欧美九色日韩亚洲蝌蚪91| 乱人伦中国视频| 亚洲一卡2卡3卡4卡5卡精品中文| 每晚都被弄得嗷嗷叫到高潮| 亚洲一卡2卡3卡4卡5卡精品中文| 女人被躁到高潮嗷嗷叫费观| 欧美日韩福利视频一区二区| 悠悠久久av| 久久久久九九精品影院| 亚洲 欧美 日韩 在线 免费| 亚洲人成电影观看| 黄网站色视频无遮挡免费观看| a级毛片在线看网站| 国产精品乱码一区二三区的特点 | 久久精品国产99精品国产亚洲性色 | 丰满饥渴人妻一区二区三| 最近最新免费中文字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 高清黄色对白视频在线免费看| 校园春色视频在线观看| 日本欧美视频一区| www国产在线视频色| 亚洲少妇的诱惑av| 免费不卡黄色视频| 国产亚洲精品第一综合不卡| 三级毛片av免费| 999久久久精品免费观看国产| 电影成人av| 老司机深夜福利视频在线观看| 亚洲人成77777在线视频| 怎么达到女性高潮| 丰满的人妻完整版| www.精华液| 国产三级在线视频| 麻豆国产av国片精品| 国产精品乱码一区二三区的特点 | 亚洲国产欧美一区二区综合| 亚洲一区中文字幕在线| 国产伦人伦偷精品视频| 午夜91福利影院| 欧美日本中文国产一区发布| 夜夜躁狠狠躁天天躁| 国产高清激情床上av| 性少妇av在线| 午夜免费成人在线视频| 亚洲欧美激情综合另类| 日韩欧美国产一区二区入口| 热99re8久久精品国产| 精品国产超薄肉色丝袜足j| 免费搜索国产男女视频| 激情在线观看视频在线高清| 男女下面插进去视频免费观看| 精品国产国语对白av| 国产一区二区在线av高清观看| 亚洲情色 制服丝袜| 久热这里只有精品99| 久久久久国产一级毛片高清牌| 日韩精品青青久久久久久| 久久九九热精品免费| 黄色成人免费大全| 一个人观看的视频www高清免费观看 | 啦啦啦 在线观看视频| 精品一区二区三区av网在线观看| 午夜福利,免费看| 国内久久婷婷六月综合欲色啪| 精品一区二区三卡| 亚洲第一欧美日韩一区二区三区| 亚洲成人免费av在线播放| 亚洲人成伊人成综合网2020| www.自偷自拍.com| 国产精品美女特级片免费视频播放器 | 成人国产一区最新在线观看| 国产精品秋霞免费鲁丝片| 波多野结衣高清无吗| 午夜老司机福利片| 亚洲成人精品中文字幕电影 | www.自偷自拍.com| 欧美黄色淫秽网站| 国产xxxxx性猛交| 婷婷丁香在线五月| 日韩欧美一区视频在线观看| 男女高潮啪啪啪动态图| 99精品欧美一区二区三区四区| 国产精品一区二区免费欧美| 国产一区二区三区综合在线观看| 色播在线永久视频| 18禁美女被吸乳视频| 又黄又爽又免费观看的视频| 在线观看午夜福利视频| 麻豆久久精品国产亚洲av | 无遮挡黄片免费观看| 女人被狂操c到高潮| 亚洲第一av免费看| 在线播放国产精品三级| 亚洲熟妇中文字幕五十中出 | 亚洲精品中文字幕一二三四区| 又黄又爽又免费观看的视频| 这个男人来自地球电影免费观看| 成人影院久久| 日韩有码中文字幕| 久久久久久免费高清国产稀缺| 久久久久九九精品影院| 欧美日韩亚洲高清精品| 国产麻豆69| 国产野战对白在线观看| 亚洲avbb在线观看| 欧美老熟妇乱子伦牲交| 日日爽夜夜爽网站| 精品国产乱码久久久久久男人| 亚洲va日本ⅴa欧美va伊人久久| 人成视频在线观看免费观看| 亚洲中文字幕日韩| 国产极品粉嫩免费观看在线| 涩涩av久久男人的天堂| 美女扒开内裤让男人捅视频| 午夜福利,免费看| 国产精品电影一区二区三区| 免费高清在线观看日韩| 不卡av一区二区三区| 电影成人av| 夜夜看夜夜爽夜夜摸 | 无限看片的www在线观看| 又黄又粗又硬又大视频| 老鸭窝网址在线观看| 国产视频一区二区在线看| 男人操女人黄网站| 午夜免费成人在线视频| 成人亚洲精品av一区二区 | 看片在线看免费视频| 精品久久久久久久毛片微露脸| 中文字幕精品免费在线观看视频| 夜夜爽天天搞| 超碰97精品在线观看| 在线观看免费午夜福利视频| 天天躁狠狠躁夜夜躁狠狠躁| 免费av中文字幕在线| 国产精品免费一区二区三区在线| 日日夜夜操网爽| 99精国产麻豆久久婷婷| 手机成人av网站| 国产精品 国内视频| 一区福利在线观看| 精品欧美一区二区三区在线| 亚洲精品一二三| 亚洲 欧美一区二区三区| 亚洲狠狠婷婷综合久久图片| 9热在线视频观看99| www.www免费av| 午夜视频精品福利| 国产精品亚洲av一区麻豆| 男人舔女人下体高潮全视频| 久久国产精品人妻蜜桃| 久久久久久大精品| 色老头精品视频在线观看| 一边摸一边抽搐一进一小说| 精品电影一区二区在线| 91麻豆av在线| 深夜精品福利| 99久久99久久久精品蜜桃| 国产成人精品无人区| 成人特级黄色片久久久久久久| 国产国语露脸激情在线看| 久久久国产欧美日韩av| 成人影院久久| 淫秽高清视频在线观看| 亚洲国产精品合色在线| 两性夫妻黄色片| av天堂在线播放| 美女扒开内裤让男人捅视频| 老汉色av国产亚洲站长工具| 欧美乱色亚洲激情| 视频在线观看一区二区三区| 好男人电影高清在线观看| 91麻豆av在线| 天堂动漫精品| 一本大道久久a久久精品| 中国美女看黄片| 精品电影一区二区在线| 欧美日韩亚洲国产一区二区在线观看| 天堂动漫精品| 国产免费男女视频| 日韩大码丰满熟妇| 亚洲伊人色综图| 国产激情久久老熟女| 日本精品一区二区三区蜜桃| 午夜91福利影院| 一级片免费观看大全| 欧美日韩亚洲高清精品| 精品日产1卡2卡| 男女床上黄色一级片免费看| 亚洲专区字幕在线| 亚洲五月婷婷丁香| 欧美日本亚洲视频在线播放| 国产精品一区二区在线不卡| 久久久久久大精品| 一二三四在线观看免费中文在| 精品熟女少妇八av免费久了| 国产又色又爽无遮挡免费看| 国产99久久九九免费精品| 免费女性裸体啪啪无遮挡网站| 麻豆成人av在线观看| 黄色片一级片一级黄色片| 男女下面进入的视频免费午夜 | 欧美色视频一区免费| 久久人人精品亚洲av| 亚洲精品国产色婷婷电影| videosex国产| av免费在线观看网站| 午夜精品久久久久久毛片777| 搡老岳熟女国产| 亚洲在线自拍视频| 欧美日韩乱码在线| 黄网站色视频无遮挡免费观看| 日本vs欧美在线观看视频| 在线永久观看黄色视频| 在线国产一区二区在线| 中亚洲国语对白在线视频| 日本免费一区二区三区高清不卡 | 水蜜桃什么品种好| 在线观看免费午夜福利视频| av欧美777| 亚洲色图综合在线观看| 高清黄色对白视频在线免费看| 亚洲国产欧美日韩在线播放| 老司机午夜福利在线观看视频| 亚洲,欧美精品.| 午夜91福利影院| 又大又爽又粗| 亚洲美女黄片视频| 伦理电影免费视频| 99精品欧美一区二区三区四区| 精品国产乱子伦一区二区三区| 国产黄a三级三级三级人| 国产1区2区3区精品| 亚洲av电影在线进入| 午夜91福利影院| 亚洲精品中文字幕一二三四区| 777久久人妻少妇嫩草av网站| 麻豆av在线久日| 中文字幕精品免费在线观看视频| 欧美乱色亚洲激情| 国产精品98久久久久久宅男小说| 男人舔女人的私密视频| 日韩大码丰满熟妇| 老司机在亚洲福利影院| 欧美成人性av电影在线观看| 一本大道久久a久久精品| 黑丝袜美女国产一区| 亚洲午夜理论影院| 国产三级黄色录像| 美女福利国产在线| 国产一区二区激情短视频| 亚洲 欧美一区二区三区| 99在线人妻在线中文字幕| 免费一级毛片在线播放高清视频 | 自线自在国产av| 人人妻人人澡人人看| 嫩草影视91久久| 色婷婷av一区二区三区视频| 亚洲欧美日韩无卡精品| 欧美另类亚洲清纯唯美| 国产一卡二卡三卡精品| xxxhd国产人妻xxx| 亚洲av电影在线进入| 女性被躁到高潮视频| 国产精品美女特级片免费视频播放器 | 日韩欧美一区二区三区在线观看| 韩国av一区二区三区四区| 午夜免费激情av| 搡老岳熟女国产| 日日摸夜夜添夜夜添小说| 老鸭窝网址在线观看| 精品人妻1区二区| 淫妇啪啪啪对白视频| 亚洲精品中文字幕一二三四区| 夜夜看夜夜爽夜夜摸 | 成人影院久久| 欧美日韩av久久| 1024视频免费在线观看| 中文亚洲av片在线观看爽| 一区在线观看完整版| 国产三级在线视频| 多毛熟女@视频| 国产午夜精品久久久久久| 亚洲欧美一区二区三区久久| 黄频高清免费视频| 国产乱人伦免费视频| 午夜老司机福利片| 亚洲男人天堂网一区| 欧美日韩瑟瑟在线播放| 国产麻豆69| 亚洲久久久国产精品| 国产精品国产高清国产av| 亚洲人成伊人成综合网2020| 国产成年人精品一区二区 | 男女之事视频高清在线观看| 日本免费一区二区三区高清不卡 | 精品国产超薄肉色丝袜足j| 成人亚洲精品av一区二区 | 久久国产精品影院| www国产在线视频色| videosex国产| 欧美日韩黄片免| 黑人猛操日本美女一级片| 叶爱在线成人免费视频播放|