999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

沒食子酸還原六價鉻反應動力學規律研究

2015-11-19 07:42:52陳子方焦維琦趙勇勝杜鵬程吉林大學地下水資源與環境教育部重點實驗室吉林長春130021
中國環境科學 2015年5期

陳子方,焦維琦,白 靜,趙勇勝,杜鵬程 (吉林大學,地下水資源與環境教育部重點實驗室,吉林 長春 130021)

沒食子酸還原六價鉻反應動力學規律研究

陳子方,焦維琦,白 靜,趙勇勝*,杜鵬程 (吉林大學,地下水資源與環境教育部重點實驗室,吉林 長春 130021)

實驗研究了不同pH值,溫度和六價鉻初始濃度等條件對沒食子酸還原六價鉻的反應動力學的影響,分別建立六價鉻還原反應速率常數關于氫離子濃度和溫度的函數方程.結果表明:pH值 為2.0~5.0時,沒食子酸與六價鉻發生的氧化還原反應符合準一級動力學反應;而pH值為6時,因[H+]不足,二者的反應不符合準一級動力學反應.在pH值為2.0~5.0時,六價鉻反應速率常數(103kobs)分別是是951.6、103.6、17.3和7.5h-1.pH值 為2.5時,溫度升高,反應速率常數快速增加;溫度為30℃時,沒食子酸還原六價鉻的反應速率常數分別是溫度為20℃和10℃時反應速率常數的1.61倍和3.68倍.不同初始濃度六價鉻的反應體系中,六價鉻反應速率常數隨著反應體系中六價鉻與沒食子酸初始濃度比值的增加而逐漸減小.利用最小二乘線性回歸法和阿倫尼烏斯公式分別建立了六價鉻還原反應速率常數關于氫離子濃度和溫度的函數方程,為預測和分析沒食子酸還原六價鉻反應過程中六價鉻的反應速率常數和濃度變化情況提供動力學模型參考.

沒食子酸;六價鉻還原;動力學;影響因素

可變價重金屬鉻造成的土壤和地下水污染是環境中最嚴重,最普遍的重金屬污染問題之一.在全球范圍內,中國和美國等許多國家將鉻列為優先控制污染物[1-2].環境中鉻的主要來源有2種:自然來源和人為污染;其中,鉻在工業生產被廣泛應用,如染料,化工,鍍鉻,木材防腐,工業造紙,制革加工等行業產生的廢水含有大量的鉻,這些工業廢水排放到環境中導致土壤,地表水和地下水污染[3-4].鉻在環境中普遍且穩定存在的形態主要有2種:三價鉻和六價鉻[5]; Nakayama等[6]報道了環境中的重金屬鉻的總量中有25%~40%以六價鉻或)形態存在,10%~20%以三價鉻(Cr3+)形態存在,另外有45%~60%是以三價鉻的有機絡合物形式存在.

鉻在環境中對動植物和人類的毒害性主要由其存在的價態及其濃度決定.六價鉻是環境中毒性最強的一種價態形式,可導致生物癌變,細胞突變以及人體器官損傷;且在水體中溶解性好,易遷移,污染易擴散[7-8].而三價鉻的毒性小,環境中溶解度低,并且三價鉻是動植物代謝的必需元素并且在糖和脂肪的代謝中占有非常重要的作用,但過量的三價鉻也會導致人體器官病變[9-10].由于六價鉻和三價鉻的毒害性和理化性質的差異,將毒害性大,遷移性強的六價鉻還原為毒性小,遷移性差的三價鉻是去除地下水中六價鉻污染的主要方法.

化學還原是將環境中六價鉻轉化為三價鉻的主要方法.目前,在實驗室研究和場地應用中常用的化學還原劑可以分為:基于鐵的還原劑(鐵屑、納米零價鐵和改性納米零價鐵等)[11-14],基于硫化物的還原劑(硫代硫酸鈉,硫化鈉等)[15-17]和基于亞鐵離子的還原劑(硫化亞鐵,亞鐵鹽等)[18-19].由于還原劑的物理化學性質和實際污染場地的水文地質條件的限制,很多化學還原劑在實際應用中修復效果與實驗室中研究結果相差很大,而且不少還原劑在修復中會產生經濟成本過高以及環境二次污染問題.

然而,植物派生性有機物因其在植物組織及腐殖質中廣泛分布,且具有環境友好性以及成本低廉等優點,而成為還原六價鉻的新型處理劑.其中, Arakawa等[20]以橡樹落葉的堿性提取液為處理劑,在酸性條件下將六價鉻還原為三價鉻. Okello等[21]研究了槲皮素(黃酮醇類化合物)還原六價鉻的反應機理,利用PbNP(鈀納米顆粒)作為催化劑,槲皮素分子中的羥基能夠還原六價鉻為三價鉻,降低反應體系pH值或升高反應溫度可以有效提高六價鉻反應速率常數. Xu等[22]利用維生素C將六價鉻還原為三價鉻,并對影響反應的濃度、溫度、pH值、光照、反應時間等因素進行了分析,認為維生素C能夠在一個較為廣泛的pH值范圍內將受污染的土壤和地下水中的六價鉻還原為三價鉻,且不受光照影響.沒食子酸是自然界存在的一種多酚類化合物,廣泛存在于掌葉大黃、大葉桉、山茱萸等植物中,在食品、生物、醫藥、化工等領域有廣泛的應用. Nakayasu等[23]確定了沒食子酸作為還原劑將六價鉻還原為三價鉻的可行性,并檢測到大分子聚合物的形成. Chen等[23]確定沒食子酸還原六價鉻的反應機理,論證了降低pH值以及升高溫度有利于促進六價鉻還原反應的進行.然而,沒食子酸還原六價鉻的反應速率與pH值和溫度等因素的定量關系還沒有文獻進行詳細的研究.

沒食子酸的分子結構中酚羥基在酸性反應體系中會失去氫原子而將六價鉻還原為三價鉻,自身形成苯醌基.本文通過實驗研究了不同pH值、溫度和六價鉻初始濃度對沒食子酸還原六價鉻的反應動力學規律影響,分別建立六價鉻還原反應速率常數與[H+]和溫度的函數關系,為分析和預測沒食子酸還原六價鉻的反應過程提供數學模型參考.

1 材料與方法

1.1 主要藥品

重鉻酸鉀、濃硫酸、濃磷酸等為優級純試劑,二苯碳酰二肼、無水碳酸鈉、沒食子酸為分析純試劑,以上藥品購自北京化工廠或國藥集團化學試劑有限公司;福林酚試劑購自北京鼎國生物科技有限公司等

1.2 分析方法及儀器

六價鉻測定方法為二苯碳酰二肼分光光度法,儀器型號:尤尼克7200可見光分光光度計;總鉻測定方法為原子吸收分光光度法,儀器型號:島津AA-6300;沒食子酸測定方法為:福林酚法,儀器型號為:WFZUV-2802H紫外分光光度計.

1.3 實驗方法

1.3.1 pH值 室溫條件下(20℃),設置5組實驗,分別稱量濃度為50mg/L溶液體積為100mL的沒食子酸溶液和濃度為20mg/L體積為100mL的重鉻酸鉀溶液先后置于棕色反應瓶中,然后利用0.5mol/L的硫酸溶液和1.0mol/L的氫氧化鈉溶液將五個反應瓶中反應體系的pH值分別調整為2.0、3.0、4.0、5.0和6.5(±0.1),混合均勻,根據不同反應時間依次取樣分析六價鉻和沒食子酸濃度.

1.3.2 溫度 實驗前將沒食子酸溶液的pH值調整為2.5(±0.1),設置3組實驗,分別量取100mL濃度為50mg/L的沒食子酸溶液和100mL濃度為20mg/L的重鉻酸鉀溶液置于棕色反應瓶中,混合均勻,將3組實驗分別置于10, 20和30℃(±1℃),根據不同的反應時間取樣分析反應體系中六價鉻和沒食子酸的濃度.

1.3.3 六價鉻初始濃度 室溫條件下(20℃),先將沒食子酸溶液的pH值調整為2.5(±0.1),設置3組實驗,分別量取100mL濃度為50mg/L的沒食子酸溶液和100mL濃度分別為10, 25和50mg/L的重鉻酸鉀溶液置于3個棕色反應瓶中,混合均勻,根據不同的反應時間取樣分析反應體系中六價鉻和沒食子酸的濃度.

1.4 動力學模型

很多學者[25-27]已經確定了六價鉻還原反應能夠用準一級動力學反應描述,因此假設沒食子酸還原六價鉻的反應符合準一級動力學規律,方程式如下:

式中: ν是六價鉻還原的反應速率,mg/(L·h); c是還原反應中六價鉻濃度, mg/L; t是反應時間, h;kobs是一級反應速率常數, h-1.式(1)兩邊求自然對數,整理可以得到式(2):

因此,六價鉻還原反應的半衰期能夠通過式(3)計算:

2 結果與討論

2.1 pH值對六價鉻還原反應的影響

如圖1所示, pH值為2.0~5.0時,沒食子酸還原六價鉻的反應基本符合準一級動力學反應.在酸性條件下特別是氫離子濃度足夠大時,沒食子酸與六價鉻的反應能夠很好的用準一級動力學反應描述.隨著反應體系中pH值逐漸減小,六價鉻還原反應的速率常數迅速增大; 其中pH值為2.0時六價鉻還原反應的速率常數分別是pH值為3.0, 4.0, 5.0時六價鉻還原反應速率常數的9.19, 55.01和126.88倍(表1).

當pH值為6.5時,實驗結果擬合得到的動力學曲線的相關系數R2為0.5005,因此沒食子酸還原六價鉻的反應在pH值為6.5時不能用準一級動力學反應表述.這說明pH值在沒食子酸還原六價鉻的反應過程中起到了非常重要的作用,足夠的氫離子濃度有利于沒食子酸還原六價鉻的反應進行,而pH值接近7.0時,沒食子酸自身的不穩定將影響六價鉻還原反應的進行.

圖1 不同pH值對沒食子酸還原六價鉻反應動力學影響Fig.1 Kinetics of Cr (VI) reduction at various initial pHvaluesa.2.0; b.3.0; c.4.0; d.5.0

表1 不同初始pH值下六價鉻反應速率常數Table 1 Rate constant and half-time at different initial pH

氫離子濃度與反應速率常數的關系無法用簡單的線性關系描述,并且該反應的反應速率常數隨著pH值的增大而以指數關系減小;因此,需要建立一個函數方程來表征六價鉻還原反應的速率常數與[H+]間的關系.通過加權最小二乘線性回歸方法得到一個合適的擬合方程能夠說明氫離子濃度與六價鉻的反應速率常數之間的定量關系[28],方程式如下:式(4)兩邊求對數,整理后得到式(5):

根據不同pH值條件下六價鉻的反應速率常數作圖(表1),可得到lg(k)和lg([H+])之間的線性關系,如圖2所示.通過直線的斜率和截距可以計算得到式(4)中a和b分別是18.09和0.7087.所以,氫離子濃度與六價鉻的反應速率常數間的函數方程式如下:

通過式(6)和動力學方程(2)能夠在不同pH值條件下(2.0~5.0)較為準確地預測和分析沒食子酸與六價鉻的反應速率常數和六價鉻濃度隨時間的變化規律.

圖2 lgk~lg[H+]關系Fig.2 Relationship of lgk and lg[H+]in the reaction

2.2 溫度對六價鉻還原反應的影響

溫度是影響沒食子酸與六價鉻氧化還原反應的主要因素之一,在10, 20和30℃研究沒食子酸與六價鉻在pH值為2.5時動力學反應規律.通過對實驗數據進行動力學擬合,該反應符合準一級動力學反應,其相關系數R2均大于98.5%(圖4和表2).在不同溫度的反應體系中六價鉻的去除率都接近100%,但是六價鉻去除所需要的時間隨著溫度的升高而逐漸減小.隨著溫度升高,反應速率常數逐漸增大;其中溫度為30℃時,沒食子酸還原六價鉻的反應速率常數分別是溫度為20℃和10℃時反應速率常數的1.61倍和3.68倍(表2).這是由于升高溫度可以增加反應體系中沒食子酸和六價鉻的活性分子數量,進而加快了沒食子酸還原六價鉻的反應速率.

圖3 溫度對六價鉻還原反應動力學影響Fig.3 Kinetics of Cr(VI) reduction at various temperatures

表2 溫度對六價鉻反應速率常數的影響Table 2 Rate constant and half-time at different temperatures

為進一步研究在沒食子酸還原六價鉻的反應過程中溫度對反應速率常數影響的定量關系,利用阿倫尼烏斯經驗公式建立了溫度與反應速率常數的函數關系式[29].阿倫尼烏斯方程式如下:

式中:Ea是反應活化能,J/mol;A是指前因子,s-1;方程式兩邊取對數,式(7)可以整理得到式(8):

根據表2中的不同溫度時六價鉻還原反應的速率常數,得到溫度與反應速率常數的lgk與1/T關系圖,如圖4所示.由直線的斜率和截距可以得到六價鉻還原反應的活化能Ea為47.16kJ/ mol,指前因子A為2.6×104s-1.因此,溫度對反應速率常數的函數關系如下:

利用式(9)和式(2),在不同溫度條件下(10~ 30℃)可以較為準確的預測和分析沒食子酸與六價鉻的反應速率常數和六價鉻濃度隨時間的變化情況,其中式(9)中活化能(Ea)越大,隨溫度升高六價鉻還原反應的速率常數增加得越快.

圖4 lgk-1/TFig.4 Relationships between lgk and 1/T in the reaction

2.3 六價鉻初始濃度對六價鉻還原反應的影響

在避光,20℃時,沒食子酸與不同初始濃度六價鉻的反應體系中六價鉻還原反應進入拖尾現象的時間幾乎是相同的,均為9h左右(圖5).反應體系中,不同初始濃度的六價鉻濃度分別減少了9.51, 17.47和17.59mg/L.這說明隨著六價鉻初始濃度增加單位時間內六價鉻還原速率增大;而當沒食子酸濃度不足時,六價鉻還原速率受沒食子酸濃度的限制.

pH值為2.5,不同初始濃度六價鉻的反應體系中,沒食子酸與六價鉻的反應基本符合一級動力學反應,該反應的相關系數R2均大于0.85.隨著六價鉻初始濃度與沒食子酸濃度比值增大,六價鉻還原反應速率常數逐漸減小;當六價鉻初始濃度與沒食子酸濃度比為1:2.5時,六價鉻還原反應的速率常數分別是濃度比為1:1和2:1時反應速率常數的9.82和24.81倍.這主要是因為在六價鉻初始濃度較大時沒食子酸濃度不足,并且隨著三價鉻的生成,沒食子酸與三價鉻在溶液中發生絡合反應形成有機絡合物,妨礙了沒食子酸還原六價鉻反應過程中電子傳遞速率[30].

沒食子酸在不同反應體系中消耗的物質的量分別是0.092, 0.168和0.168μmol,所以,在沒食子酸還原不同初始濃度的六價鉻的反應的過程中,沒食子酸的消耗的物質的量與六價鉻還原的物質的量比值分別為1:1.98, 1:2和1:2.01.這說明在強酸性條件下,不同初始濃度的六價鉻與沒食子酸反應原理基本是相同的,每摩爾沒食子酸能夠提供6個電子,將2mol的六價鉻還原為三價鉻.

表3 六價鉻初始濃度對反應速率常數的影響Table 3 Rate constant at different initial Cr (VI)concentration

3 結論

3.1 pH值為2.0~5.0時,沒食子酸還原六價鉻反應符合準一級動力學反應;隨體系pH值降低,六價鉻反應速率常數增大,通過最小二乘法得到六價鉻反應速率常數關于與[H+]的函數關系為k=18.09[H+]0.7087,該方程式的pH值適用范圍為2.0~5.0.

3.2 升高溫度,六價鉻還原反應速率常數增大;擬合阿倫尼烏斯公式可得到該反應活化能為47.16kJ/mol,指前因子A為2.6×104s-1,該方程式的溫度適用范圍為10~30℃.

3.3 pH值為2.0時,沒食子酸還原六價鉻反應的摩爾比約為1:2;反應體系中沒食子酸濃度不變,隨六價鉻與沒食子酸初始濃度比值的增大,六價鉻還原反應速率常數減小,而還原反應速率增大.

[1]聶 寧,丁遠昭,李喜青.制備斜發沸石和零價鐵復合材料處理水中的六價鉻污染 [J]. 中國環境科學, 2014,33(3):443-447.

[2]Dhal B, Thatoi H N, Das N N, et al., Chemical and microbial remediation of hexavalent chromium from contaminated soil and mining/metallurgical solid waste: a review [J]. Hazardous Materials, 2013,250-251:272-291.

[3]王玲玲,劉 琰,鄭丙輝,等.納污水體中鉻的遷移規律及各形態分配系數 [J]. 環境科學研究, 2010,23(4):426-429.

[4]Fu F, Wang Q. Removal of heavy metal ions from wastewaters: a review [J]. Environmental Management, 2011,92(3):407-18.

[5]Eary L E, Davis A. Geochemistry of an acidic chromium sulfate plume [J]. Applied Geochemistry, 2007,22(2):357-369.

[6]Nakayama E, Tokoro H, Kuwamoto T, et al. Dissolved state of chromium in seawater [J]. Nature, 1981,290:768-770.

[7]張義賢.三價鉻和六價鉻對大麥毒害效應的比較 [J]. 中國環境科學, 1997,17(6):565-568.

[8]陳子方,趙勇勝,孫家強,等.重金屬鉛和鉻在包氣帶中遷移轉化規律研究 [J]. 中國環境科學, 2014,34(9):2211-2216.

[9]裘波音.水稻鉻脅迫耐性的遺傳分析與還原型谷胱甘肽緩解鉻毒害的機理研究 [D]. 杭州:浙江大學, 2012.

[10]Shafaqat Ali.大麥鉻與鹽,鋁的互作和減輕鉻毒害的化學途徑研究 [D]. 杭州:浙江大學, 2010.

[11]秦澤敏,董黎明,劉 平,等.零價納米鐵吸附去除水中六價鉻的研究 [J]. 中國環境科學, 2014,34(12):3106-3111.

[12]Qiu X, Fang Z, Yan X, et al. Chemical stability and toxicity of nanoscale zero-valent iron in the remediation of chromiumcontaminated water shed [J]. Chemical Engineering Journal, 2013,22:61-66.

[13]朱文會,王興潤,董良飛,等.海藻酸鈉固定化Fe-Cu雙金屬去除Cr(Ⅵ)的作用機制 [J]. 中國環境科學, 2013,33(11):1965-1971.

[14]張 珍,沈珍妙,范天恩,等.腐植酸對nZVI去除水中Cr(VI)的抑制及抑制作用的消除 [J]. 中國環境科學, 2013,33(1):63-68.

[15]楊俊香,蘭葉青.硫化物還原Cr(VI)的反應動力學研究 [J]. 環境科學學報, 2005,25(3):356-360.

[16]謝翼飛,李旭東,李福德.生物硫鐵復合材料處理含鉻廢水及鉻資源化研究 [J]. 中國環境科學, 2009,29(12):1269-1275.

[17]Henderson A D, Demond A H. Permeability of iron sulfide(FeS)-based materials for groundwater remediation [J]. Water Reserch, 2013,47:1267-1276.

[18]Ludwig R D, Su C M, Lee T R, et al. In situ chemical reduction of Cr (VI) in groundwater using a combination of ferrous sulfate and sodium dithionite: A field investigation [J]. Environmental Science Technology, 2007,41:5299-5305.

[19]韓曉剛,黃廷林,陳秀珍.硫酸亞鐵還原法在突發性鉻(VI)污染應急處理中的應用研究 [J]. 工業水處理, 2013,33(6):40-46.

[20]Arakawa H, Tsushima M, Kishi M, et al. Reduction of chromium(VI) by water-extracts from withered oak leaves [J]. Chem Lett,(1993) 2113-2116.

[21]Okello V A, Mwilu S, Noah N, et al. Reduction of hexavalent chromium using naturally-derived flavonoids [J]. Environmental Science and Technology, 2012,46(19):10743-10751.

[22]Xu XR, Li HB, Li XY, et al. Reduction of hexavalent chromium by ascorbic acid in aqueous solutions [J]. Chemosphere, 2004,57:609-613.

[23]Nakayasu K, Fukushima M, Sasaki K, et al. Comparative studies of the reduction behavior of chromium (VI) by humic substances and their precursors [J]. Environmental Toxicology and Chemistry,1999,18(6):1085-1090.

[24]Chen Z F, Zhao Y S. Characteristics and kinetics of hexavalent chromium reduction by gallic acid in aqueous solutions [J]. Unpublished, 2015.

[25]Wittbrodt P R, Palmer C D. Reduction of Cr (VI) in the presence of excess soil fulvic acid [J]. Environmental Science and Technology, 1995,29:255-263.

[26]Nakayasu K, Fukushima M, Sasaki K, et al. Comparative studies of the reduction behavior of chromium (VI) by humic substances and their precursors [J]. Environmental Toxicology and Chemistry,1999,18:1085-1090.

[27]Kassem T S. Kinetics and thermodynamic treatments of the reduction of hexavalent to trivalent chromium in presence of organic sulphide compounds [J]. Desalination, 2010,258(1-3):206-218.

[28]Park D, Yun Y S, Ahn C K, et al. Kinetics of the reduction of hexavalent chromium with the brown seaweed Ecklonia biomass[J]. Chemosphere, 2007,66:939-946.

[29]Park D, Yun Y S, Jo J H, et al. Mechanism of hexavalent chromium removal by dead fungal biomass of Aspergillus niger[J]. Water Res, 2005,39:533-540.

[30]Xu G R, In M Y, Yuan Y, et al. In situ spectroelectrochemical study of quercetin oxidation and complexation with metal ions in acidic solutions [J]. Bull. Korean Chem. Soc., 28(2007) 889-892.

Kinetics of the reduction of hexavalent chromium with Gallic acid.

CHEN Zi-Fang, JIAO Wei-qi, BAI Jing, ZHAO Yong-sheng*, DU Peng-cheng (Key Lab of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China). China Environmental Science, 2015,35(5):1391~1396

The objective of this work was to determine the effect of pH, temperature and initial Cr (VI) concentration on kinetics of Cr (VI) reduction, and to develop two kinetics models for the rate constants upon pH and temperature,respectively. Results indicated that the redox reaction of gallic acid and hexavalent chromium was the pseudo-first-order kinetics reaction at pH 2.0~5.0; because H+is insufficient, the reaction was not pseudo-first-order at pH 6.5. The rate constants (103kobs) were 951.6, 103.6, 17.3and 7.5h-1respectively when pH was from 2.0to 5.0. At pH 2.5, increasing temperature could increase the kobs, and the kobsincreased by 0.61 and 2.68 times when temperature was adjusted from 10℃ to 20℃ and 30℃. The rate constants of Cr (VI) reduction were gradually decreased with increasing initial concentration ratio of Cr (VI) and gallic acid. The kinetics equation of the rate constants and [H+]was obtained through a weighted least-squares linear regression method, and Arrhenius formula was used to develop the kinetics equation of the rate constants and temperature. These provided a model reference for analysis and forecast the reaction rates and Cr (VI)concentration on reaction of Cr (VI) reduction with gallic acid.

gallic acid;Cr (VI) reduction;kinetics;influence factor

X523

A

1000-6923(2015)05-1391-06

陳子方(1987-),男,山東聊城人,吉林大學博士研究生,主要研究方向為污染場地污染控制與修復.

2014-09-15

國家環保公益項目(2013A073);北京市教育委員會市屬高校創新能力提升計劃項目(TJSH201310772028)

* 責任作者, 教授, zhaoyongsheng@jlu.edu.cn

主站蜘蛛池模板: 为你提供最新久久精品久久综合| 香蕉久人久人青草青草| 自慰网址在线观看| 色综合国产| 午夜日韩久久影院| 国产成人91精品| 国产丝袜无码精品| 国产区成人精品视频| 中文字幕欧美日韩高清| 久久6免费视频| 91免费在线看| 日本a∨在线观看| 免费 国产 无码久久久| 国产chinese男男gay视频网| 怡春院欧美一区二区三区免费 | 在线欧美日韩| 国产人成网线在线播放va| 高清视频一区| 亚洲毛片网站| 国产美女主播一级成人毛片| 日韩在线第三页| 日韩不卡高清视频| 无遮挡国产高潮视频免费观看| 国产一级小视频| 国产欧美日韩在线一区| 国产午夜不卡| 91亚洲精选| 国产十八禁在线观看免费| 久操中文在线| 国产1区2区在线观看| 一边摸一边做爽的视频17国产| 亚洲国产天堂久久综合| 欧美一区日韩一区中文字幕页| 国产高清在线观看91精品| 亚卅精品无码久久毛片乌克兰| 精品综合久久久久久97超人| 人妻91无码色偷偷色噜噜噜| 亚洲福利片无码最新在线播放| 国产精品免费p区| 国产精品自在线拍国产电影| 免费欧美一级| 亚洲中文字幕日产无码2021| 成人福利在线看| 毛片国产精品完整版| 国产性生大片免费观看性欧美| 91毛片网| 91国语视频| 国产素人在线| 欧美一级黄片一区2区| 国产在线观看成人91 | 影音先锋丝袜制服| 国产精品 欧美激情 在线播放| 亚洲欧美精品在线| 久久www视频| 亚洲黄色成人| 一级看片免费视频| 亚洲第一区在线| 日韩中文无码av超清| 制服丝袜 91视频| 国产欧美又粗又猛又爽老| 日本久久久久久免费网络| 少妇人妻无码首页| 久久精品aⅴ无码中文字幕 | 日韩中文字幕亚洲无线码| 任我操在线视频| 综合久久五月天| 亚洲AV无码久久精品色欲| 一级毛片免费不卡在线| 一级毛片不卡片免费观看| 日韩黄色精品| 国产激情无码一区二区APP| 色综合久久无码网| 亚洲福利片无码最新在线播放 | 亚洲第一视频免费在线| 国产精品手机在线观看你懂的 | 亚洲精品视频网| 国产精品太粉嫩高中在线观看| 精品在线免费播放| 天天色综网| 欧美不卡在线视频| 国产日韩欧美一区二区三区在线| 国产凹凸视频在线观看|