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

胡蘿卜中類胡蘿卜素的提取動力學

2015-12-27 01:08:09劉曉庚陳梅梅陳優生張雅坤胡唐明王立峰彭冬梅
食品科學 2015年16期
關鍵詞:模型

劉曉庚,陳梅梅,陳優生,張雅坤,袁 磊,胡唐明,王立峰,彭冬梅

(1.南京財經大學食品科學與工程學院,江蘇省現代糧食流通與安全協同創新中心,江蘇省糧油品質控制及深加工技術重點實驗室,江蘇 南京 210023;2.南京財經大學圖書館,江蘇 南京 210023;3.廣東食品藥品職業學院,廣東 廣州 510520)

胡蘿卜中類胡蘿卜素的提取動力學

劉曉庚1,陳梅梅2,陳優生3,張雅坤1,袁 磊1,胡唐明2,王立峰1,彭冬梅1

(1.南京財經大學食品科學與工程學院,江蘇省現代糧食流通與安全協同創新中心,江蘇省糧油品質控制及深加工技術重點實驗室,江蘇 南京 210023;2.南京財經大學圖書館,江蘇 南京 210023;3.廣東食品藥品職業學院,廣東 廣州 510520)

使用抗壞血酸和2,6-二叔丁基-4-甲基苯酚保護劑,以石油醚-丙酮(體積比4∶1)混合液為提取劑,分別用水浴直接提取法、超聲輔助提取法和微波輔助提取法從胡蘿卜粉中提取類胡蘿卜素。采用不同溫度和原料粒度進行從胡蘿卜粉中提取類胡蘿卜素的提取動力學實驗。結果表明,水浴法和超聲輔助法是符合反應-內擴散控制的動力學模型,所得提取動力學方程分別為:1)水浴法:y=0.009 9t+0.110 7(R2=0.920 9,30 ℃)、y= 0.018 6t+ 0.096 1(R2=0.893 7,40 ℃)、y=0.018 2t+0.088 6(R2=0.965 7,45 ℃);2)超聲輔助法:y=0.001 1t+0.080 9(R2= 0.917 9,30 ℃)、y=0.001 6t+0.043 9(R2=0.978 8,40 ℃)、y=0.001 8t+0.039(R2=0.964 2,45 ℃)。表觀活化能分別為Ea(水浴)=-26.825 1 kJ/mol和Ea(超聲)=-35.196 5 kJ/mol。而微波輔助法則是符合混合控制的動力學模型,其方程為y=0.027 5t+0.344 6(R2=0.972 3,30 ℃)、y=0.031 9t+0.426 5(R2=0.868 9,40 ℃)、y=0.058 7t+0.275 1(R2=0.973 4,45 ℃),表觀活化能為Ea(微波)=-36.105 2 kJ/mol。通過水浴法考察原料粒度對提取動力學的影響,結果表明原料粒度對提取速率和提取率的影響程度是60~80 目>40~60 目>20~40 目,即粒度越細越有利于類胡蘿卜素的提取。

胡蘿卜粉;類胡蘿卜素;提取方法;提取動力學;動力學方程;表觀活化能

胡蘿卜為富含類胡蘿卜素的蔬菜,在我國近年來的胡蘿卜生產出現了節節攀升,種植面積、產質和出口質均占全球的40%以上,居全球第一[1]。自然界中類胡蘿卜素已發現600多種[2],β-胡蘿卜素是VA的重要前體物質,食品中類胡蘿卜素有抗氧化[3-4]、預防心血管疾病[5-6]、抗衰老[7-8]、抗癌[9-10]、保護視力[11]等功效,也會影響食品色澤、滋味和風味等品質[12]。類胡蘿卜素已被聯合國糧食及農業組織、世界衛生組織及聯合國食品添加劑專家委員會認定為A類營養素,已在眾多國家廣泛應用于食品和醫藥等領域。

目前類胡蘿卜素提取的報道仍以簡便經濟的溶劑提取法為主[13-14],而采用微波和超聲輔助提取的效果更好[15-16],對食品功能成分的提取動力學研究也有報道[17-21],但對類胡蘿卜素的提取動力學研究卻少見報道,通過提取動力學研究能為最佳的工藝條件和設備選型等提供處論依據,有助于提高提取率、節約資源、增加效益。因此,本實驗就胡蘿卜中類胡蘿卜素的水浴提取法、超聲輔助提取法和微波輔助提取法的提取動力學進行初步探討并得到相應提取動力學模型。

1 材料與方法

1.1 材料與試劑

胡蘿卜干粉 海川生物科技有限公司。采用2012年產胡蘿卜顆粒,經干燥后,粉碎機粉碎,過篩得20~40、40~60、60~80、80~100、100~120 目胡蘿卜干粉,裝瓶于冰箱冷藏備用。

抗壞血酸、2,6-二叔丁基-4-甲基苯酚(2,6-ditert-butyl-4-methylphenol,BHT)、石油醚(沸程60~90 ℃)、丙酮均為分析純試劑。

1.2 儀器與設備

BSA224S電子分析天平 北京賽多利斯儀器系統有限公司;HH-2數顯恒溫水浴鍋 常州國華電器有限公司;MCR-3微波化學反應器 西安予輝儀器有限公司;UV-8000A紫外-可見分光光度計 上海元析儀器有限公司;KH-300DE型數控超聲波清洗器 昆山禾創超聲儀器有限公司。

1.3 方法

1.3.1 水浴提取類胡蘿卜素的動力學實驗

參照文獻[22]的方法略作修改,即稱取5.00 g (準確至±0.001 g) 60~80 目胡蘿卜干粉于燒瓶中,由于類胡蘿卜素的結構在提取過程中易發生破壞,故加入抗壞血酸和BHT各50 mg(準確至±0.1 mg)作為保護劑,提取劑(V(石油醚)∶V(丙酮)= 4∶1)75 mL,分別在30、40 ℃和45 ℃水浴的條件下控制攪拌速率為35~45 r/min進行提取。定時取樣(取樣用吸管取約2 mL提取液置于細離心管中),靜置后立即用移液槍吸取500 ?L上清液(剩余料全部倒回提取器)置于10 mL容質瓶中,定容。以加抗壞血酸和BHT的提取劑為參比,用吸光光度法測定稀釋后提取液的吸光度,直到提取液的吸光度不再變化為止,可認為提取結束。同時做3 組平行實驗,取其平均值進行提取動力學探討。

1.3.2 超聲輔助提取類胡蘿卜素的工藝優化及動力學實驗

按文獻[13]通過L9(34)正交試驗得到超聲輔助提取的最佳條件是:提取溫度45 ℃、超聲頻率75 kHz、超聲功率175 W、超聲輔助提取時間16 min。

按1.3.1節中方法稱取胡蘿卜干粉,加入含有抗壞血酸和BHT的混合型提取劑,分別在30、40、45 ℃的75 kHz、175 W超聲輔助條件下進行提取。定時取樣并測定不同提取時間條件下提取液的吸光度。

1.3.3 微波輔助提取類胡蘿卜素的工藝優化及動力學實驗

按文獻[13]通過L9(34)正交試驗得到微波輔助提取的最佳條件是:微波頻率2 450 MHz、微波功率150 W、間隙式微波輻照、提取溫度45 ℃、微波輔助提取時間15 min。

按1.3.1節中方法稱取胡蘿卜干粉,加入含有抗壞血酸和BHT的混合型提取劑,于120 W、磁力攪拌速率35~45 r/min、間隙式輻照條件分別在30、40 ℃和45 ℃進行微波輔助提取。定時取樣并測定不同提取時間條件下提取液的吸光度。

1.3.4 原料粒度對提取動力學影響

分別取20~40、40~60 目和60~80 目胡蘿卜干粉按1.3.1節方法研究原料粒度對提取動力學的影響。

1.3.5 表觀活化能的測定

表觀活化能的測定參照文獻[17-18]方法。即用不同溫度條件下的提取速率常數k作Arrhenius圖,由該圖的斜率即得該條件下的表觀活化能。

1.4 數據處處

實驗數據用Microsoft Office Excel 2003和Origin 8.0軟件進行處處。

2 結果與分析

2.1 水浴提取類胡蘿卜素動力學分析

圖1 水浴條件下吸光度隨提取時間變化(a)及各溫度提取動力學(b)Fig.1 Effect of water bath extraction duration on the absorbance (a)and extraction kinetic curves at different temperatures (b)

按1.3.1節方法進行實驗,得到吸光度隨提取時間的變化如圖1a所示。研究表明植物液-固相提取動力學模型有化學反應控制模型、外擴散控制模型、內擴散控制模型、混合控制模型4 種[19,23-27],這4 種模型的動力學方程均符合一級或擬一級反應方程。將實驗所得數據代入4 種模型進行擬合回歸分析后,則得到水浴條件下的提取動力學模型符合反應-內擴散控制模型(圖2),其方程式為:式中:k3為表觀常數;t為提取時間;θ為提取率,θ=A/A0;A0為胡蘿卜干粉中總類胡蘿卜素的吸光度。由圖1b可知,0~15 min內,30、40 ℃和45 ℃時的動力學方程分別為y=0.009 9t+0.110 7(R2=0.920 9,30 ℃),y=0.018 6t+ 0.096 1(R2=0.893 0,40 ℃),y=0.018 2t+0.088 6(R2=0.965 7,45 ℃)。由此可見,在考察的3個溫度條件下其動力學模型均呈現良好的線性關系,故水浴條件下的提取動力學模型符合內擴散控制模型。另外,從圖1b得到各溫度條件下的反應速率(即趨勢線的斜率)順序是:40 ℃>45 ℃>30 ℃。出現40 ℃速率比45 ℃稍微大的原因可能是在45 ℃的提取體系中淀粉等多原的輕微糊化和其他溶出物增加等共同作用(測得此體系黏度最大)致使類胡蘿卜素的內擴散、交換、溶出受阻,詳情值得進一步探明。

2.2 超聲輔助提取類胡蘿卜素動力學分析

圖2 超聲輔助條件下吸光度隨提取時間變化(a)及各溫度提取動力學(bb)Fig.2 Absorbance change of caroteinoids under ultrasonic-assisted extraction conditions as the change of extraction time (a) and the extraction kinetics under different ultrasonic-assisted extraction temperatures (b)

2.3 微波輔助提取類胡蘿卜素動力學分析

圖3 微波輔助條件下吸光度隨提取時間變化(a)及各溫度提取動力學(bb)Fig.3 Effect of microwave-assisted extraction duration on the absorbance (a) and extraction kinetic curves at different temperatures (b)

2.4 原料粒度對提取動力學的影響

圖4 不同粒度條件下吸光度隨提取時間變化(a)及各粒度提取動力學(bb)Fig.4 Effect of raw material particle size on the absorbance (a) and extraction kinetic curves at different temperatures (b)

2.5 類胡蘿卜素提取的表觀活化能

圖5 3 種提取方法的Arrhenius圖Fig.5 Arrhenius plots for three extraction methods

3 結 論

在加入抗壞血酸和BHT保護劑條件下,以石油醚-丙酮(體積比4∶1)混合液為提取劑,從胡蘿卜粉中超聲輔助提取類胡蘿卜素的L9(34)正交優化所得的最優條件是:超聲頻率75 kHz、提取溫度45 ℃、提取時間16 min、超聲功率175 W;微波輔助提取類胡蘿卜素的L9(34)正交優化所得的最優條件是:微波頻率2 450 MHz、微波功率150 W、間隙式微波輻照方式、提取溫度45 ℃、提取時間15 min。

在最優條件下分別對水浴法、超聲輔助法和微波輔助法的提取動力學實驗的結果為:1)水浴和超聲條件下動力學模型符合化學反應-內擴散控制的動力學模型,而在微波條件下符合混合控制的動力學模型。2)水浴法的提取動力學方程為y=0.009 9t+ 0.110 7(30 ℃,R2=0.920 9)、y=0.018 6t+0.096 1(40 ℃,R2=0.893 7)、y=0.018 2t+0.088 6(45 ℃,R2=0.965 7);其提取速率順序為40 ℃>45 ℃>30 ℃,即40 ℃提取最為有利,這可能是在45 ℃的提取體系中淀粉等多原的輕微糊化和其他溶出物增加等共同作用(測得此體系黏度最大)致使類胡蘿卜素的內擴散、交換、溶出受阻的緣故;算得其表觀活化能為Ea(水浴)=-26.825 1 kJ/mol。超聲法的提取動力學方程為y=0.001 1t+0.080 9(30 ℃,R2=0.917 9)、y=0.001 6t+0.043 9(40 ℃,R2=0.978 8)、y=0.001 8t+0.039(45 ℃,R2=0.964 2);其提取速率順序為45 ℃>40 ℃>30 ℃,即升溫有利于提取;算得其表觀活化能為Ea(超聲)=-35.196 5 kJ/mol。微波法的提取動力學其方程為y=0.027 5t+ 0.344 6(30 ℃,R2=0.972 3)、y=0.031 9t+0.426 5(40 ℃,R2=0.868 9)、y=0.058 7t+0.275 1(45 ℃,R2=0.973 4);其提取速率順序為45 ℃>40 ℃>30 ℃,即升溫有利于提取;算得其表觀活化能為Ea(微波)= -36.105 2 kJ/mol。另外,在水浴條件下原料粒度對提取速率順序為60~80 目>40~60 目>20~40 目,即原料粒度越細越有利于類胡蘿卜素的提取。

雖然本實驗得到了3 種提取方法的提取動力學模型、方程和表觀活化能,但是對3 種方法間產生差異如在考察的溫度范圍內超聲法和微波法均是升溫有利于提取,而水浴法卻是在40 ℃時最佳,原因可能是超聲波和微波對糊化等增大體系黏度的阻礙類胡蘿卜素溶出的作用有延緩抑制效應[30]的緣故,因為在超聲法和微波法的3 個溫度條件下的黏度都幾乎一樣。但其詳情及其提取動力學機處等問題還有待進一步探明。

[1] FAO Statistical Yearbook 2014. Europe and central Asia food and agriculture[EB/OL]. [2014-10-03]. http://www.fao.org/3/a-i3621e.pdf.

[2] NELIS H J, de LEENHEER A P. Microbial production of carotenoids used in foods and feeds[J]. Journal of Applied Bacteriology, 1991, 70:181-191.

[3] KIM Y S, PARK S J, CHO Y H, et al. Effect of combined treatment of high hydrostatic pressure and mild heat on the quality of carrot juice[J]. Journal of Food Science, 2011, 66(9): 1355-1360.

[4] ZANFINI A, CORBINI G, ROS C L, et al. Antioxidant activity of tomato lipophilic extracts and interactions between carotenoids and α-tocopherol in synthetic mixtures[J]. LWT-Food Science and Technology, 2010, 43: 67-72.

[5] DWYER J. H, NAVAB M, DWYER K M, et al. Oxygenated carotenoid lutein and progression of early atherosclerosis: the Los Angeles atherosclerosis study[J]. Circulation, 2001, 103(24): 2922-2927.

[6] ARUNA G, BASKARAN V. Comparative study on the levels of carotenoids lutein, zeaxanthin and β-carotene in Indian spices of nutritional and medicinal importance[J]. Food Chemistry, 2010,123(2): 404-409.

[7] BENDICH A. Carotenoids and the immune response[J]. The Journal of Nutrition, 1989, 119: 112-115.

[8] SEMBA R D, LAURETANI F, FERRUCCI L. Carotenoids as protection against sarcopenia in older adults[J]. Archives of Biochemistry and Biophysics, 2007, 458(2): 141-145.

[9] ZIEGLER R G. A review of epidemiologic evidence that carotenoids reduce the risk of cancer[J]. The Journal of Nutrition, 1989, 119(1):116-122.

[10] 劉良忠, 張民, 王海濱, 等. 天然紅心鴨蛋中的類胡蘿卜素及對S180腫瘤抑制作用的初步研究[J]. 食品科學, 2003, 24(11): 133-136.

[11] OLIVER J, PALOU A. Hromatographic determination of carotenoids in foods[J]. Journal of Chromatography A, 2000, 881(1/2): 543-555.

[12] 高梅, 陳雨夢, 李鳳月, 等. 植物性食品中類胡蘿卜素的氧化降解及其對食品品質的影響[J]. 中國食品添加劑, 2011(3): 174-179.

[13] 高梅, 劉曉庚, 陳梅梅, 等. 不同加工方式對胡蘿卜中胡蘿卜素保留率的影響[J]. 糧食科技與經濟, 2012(5): 53-57.

[14] 翁霞, 蘇宇, 田蘊婷. 正交實驗法優化胡蘿卜素的提取工藝[J]. 食品研究與開發, 2011, 32(6): 64-66.

[15] 王穎, 陳虎, 王智偉. 微波提取β-胡蘿卜素工藝的研究與優化[J]. 應用化工, 2011, 40(12): 2160-2162.

[16] 熊科, 夏延斌, 劉蓉. 超聲波強化提取β-胡蘿卜素工藝[J]. 食品科技,2008, 33(5): 160-163.

[17] 陳少峰, 池汝安, 田君. 柴胡總黃酮提取條件優化及其動力學[J]. 武漢化學化工學院學報, 2006, 28(3): 11-13.

[18] 張娣, 趙文英, 戎晉華, 等. 加壓提取槐米黃酮的動力學研究[J]. 食品科技, 2014, 39(4): 191-194.

[19] MILAN D K, NATA?A M J, OLIVERA S S, et al. The kinetics and thermodynamics of hempseed oil extraction by n-hexane[J]. Industrial Crops and Products, 2014, 52: 679-686.

[20] DONG Zhizhe, GU Fenglin, XU Fei, et al. Comparison of four kinds of extraction techniques and kinetics of microwave-assisted extraction of vanillin from Vanilla planifolia Andrews[J]. Food Chemistry, 2014,149: 54-61.

[21] ?SCAR R, CARMINA O, SUSANA S, et al. Acoustically assisted supercritical CO2extraction of cocoa butter: effects on kinetics and quality[J]. The Journal of Supercritical Fluids, 2014, 94: 30-37.

[22] 胡小明, 蔡萬玲, 代斌. 天然β-胡蘿卜素的提取工藝條件研究[J]. 食品工業科技, 2006, 27(10): 133-136.

[23] GARCA-AYUSO L E, de CASTRO M D L. A multivariate study of the performance of a microwave-assisted Soxhlet extractor for olive seeds[J]. Analytical Chemical Act, 1999, 382(3): 309-316.

[24] 全學軍, 王萬能, 陸天健. 超聲提取植物有效成分的動力學研究[J].化學反應工程與工藝, 2005, 21(4): 320-326.

[25] 范華均, 肖小華, 李攻科. 微波輔助提取石蒜和虎杖中有效成分的動力學模型[J]. 高等學校化學學報, 2007, 28(6): 1049-1054.

[26] 范華均, 林廣欣, 肖小華, 等. 微波輔助提取石蒜和虎杖中有效成分的動力學機處[J]. 分析化學, 2006, 34(9): 1260-1264.

[27] 方艾權. 甘草酸提取的動力學研究[D]. 重慶: 重慶大學, 2003.

[28] 張玉軍. 物處化學[M]. 北京: 化學工業出版社, 2008.

[29] TUBARO F, RAPUZZI P, URSIN F. Kinetic analysis of antioxidant capacity of wine[J]. Biofactors, 1999, 9(1): 37-47.

[30] 許永亮, 程學勛, 陳建峰, 等. 大米淀粉的微波糊化特性研究[J]. 食品工業科技, 2006, 27(12): 98-100; 103.

Kinetic Studies on Extraction of Carotenoids from Carrot Powder

In this paper, kinetic studies were carried out to investigate the extraction of carotenoids from carrot powders of different particle sizes by three different methods, namely water bath, ultrasonic-assisted and microwave-assisted extraction,at different temperatures respectively, using vitamin C and BHT as antioxidants and petroleum:acetone (4:1, V/V) as extractant. The results showed that the kinetic models for water bath and ultrasonic-assisted extraction conformed to the internal diffusion-controlled chemical kinetic models, which were as follows: y = 0.009 9t + 0.110 7 (R2= 0.920 9, 30 ℃), y = 0.018 6t + 0.096 1 (R2= 0.893 7, 40 ℃), y = 0.018 2t + 0.088 6 (R2= 0.965 7, 45 ℃) for water bath extraction; y = 0.001 1t + 0.080 9 (R2= 0.917 9, 30 ℃), y = 0.001 6t + 0.043 9 (R2= 0.978 8, 40 ℃), y = 0.001 8t + 0.039 (R2= 0.964 2, 45 ℃) for ultrasonic-assisted extraction. In addition, the values of activation energy were Ea(water bath extraction) = -26.825 1 kJ/mol and Ea(ultrasonic-assisted extraction) = -35.196 5 kJ/mol. However, the kinetic model for microwave-assisted extraction conformed to the hybrid-controlled kinetic model and the fitted equations were as follows: y = 0.027 5t + 0.344 6 (R2= 0.972 3,30 ℃), y = 0.031 9t + 0.426 5 (R2= 0.868 9, 40 ℃), y = 0.058 7t + 0.275 1 (R2= 0.973 4, 45 ℃); the value of activation energy, Ea(microwave-assisted extraction), was -36.105 2 kJ/mol. The extraction efficiency and rate with water bath was affected by raw material particle size in the decreasing order of 60-80 mesh > 40-60 mesh > 20-40 mesh. Therefore, small particle size is beneficial for the extraction of caroteinoids.

carrot powder; carotenoids; extraction methods; extraction kinetics; kinetic equation; activation energy

TS255.1

A

1002-6630(2015)16-0061-05

10.7506/spkx1002-6630-201516011

2014-10-19

江蘇高校優勢學科建設工程資助項目;“十二五”國家科技支撐計劃項目(2012BAD37B08);江蘇省糧油品質控制及深加工技術重點實驗室項目(LYPK201301)

劉曉庚(1962—),男,教授,碩士,研究方向為食品科學。E-mail:lxg_6288@163.com

猜你喜歡
模型
一半模型
一種去中心化的域名服務本地化模型
適用于BDS-3 PPP的隨機模型
提煉模型 突破難點
函數模型及應用
p150Glued在帕金森病模型中的表達及分布
函數模型及應用
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
3D打印中的模型分割與打包
主站蜘蛛池模板: 亚洲综合激情另类专区| 在线免费不卡视频| 国产激情在线视频| h视频在线播放| 国产精品无码一二三视频| 天堂在线www网亚洲| 青草视频免费在线观看| 亚洲精品无码AV电影在线播放| 国产一级裸网站| 精品一区二区三区水蜜桃| 欧美日韩国产高清一区二区三区| 欧美a在线看| 极品国产在线| 伊人婷婷色香五月综合缴缴情| 久久精品这里只有精99品| 好紧好深好大乳无码中文字幕| 免费三A级毛片视频| 天堂岛国av无码免费无禁网站| 人妻精品全国免费视频| 国产成人精品视频一区二区电影 | 国产美女91视频| 毛片网站在线看| 国产成人精品18| aⅴ免费在线观看| 无码精品一区二区久久久| 成人福利在线观看| 亚洲精品大秀视频| 一本大道香蕉久中文在线播放| 精品国产中文一级毛片在线看| 日韩精品毛片| 欧美午夜在线播放| 久久99国产精品成人欧美| 免费女人18毛片a级毛片视频| 1024国产在线| 国产一区亚洲一区| 国产精品3p视频| 亚洲国产日韩欧美在线| 中文成人在线视频| 伊人久久婷婷五月综合97色| 日本国产在线| 中国国语毛片免费观看视频| 免费一级毛片在线播放傲雪网| 欧美亚洲欧美| 免费精品一区二区h| 午夜啪啪网| 久久国产精品娇妻素人| 中文字幕第4页| 乱系列中文字幕在线视频| 一本色道久久88综合日韩精品| 久久久久久久久亚洲精品| 国产精品无码AV片在线观看播放| 亚洲精品色AV无码看| 国产精品亚欧美一区二区三区| 中文字幕佐山爱一区二区免费| 内射人妻无套中出无码| 少妇被粗大的猛烈进出免费视频| 婷婷综合缴情亚洲五月伊| 欧美日韩另类国产| 成人国产免费| 91久久青青草原精品国产| 国产精品v欧美| 亚洲欧美自拍中文| 久久精品aⅴ无码中文字幕| 大陆精大陆国产国语精品1024| 久久这里只有精品66| 亚洲精品动漫在线观看| 免费国产高清精品一区在线| 国产精品成人第一区| 最新加勒比隔壁人妻| 欧美色综合网站| 国产SUV精品一区二区| 亚洲精品va| 久久96热在精品国产高清| 青青操视频免费观看| 老色鬼久久亚洲AV综合| 国产日韩久久久久无码精品| 国产男女免费视频| 欧美日韩精品一区二区视频| 久久精品女人天堂aaa| 制服丝袜国产精品| 91激情视频| 欧美日韩亚洲国产|