亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于熱重法的植物油氧化動力學研究

        2016-09-14 07:13:38張巧智齊寶坤隋曉楠馬文君江連洲東北農(nóng)業(yè)大學食品學院黑龍江哈爾濱150030
        食品工業(yè)科技 2016年2期
        關鍵詞:植物油油脂組分

        張巧智,李 楊,齊寶坤,隋曉楠,馬文君,江連洲(東北農(nóng)業(yè)大學食品學院,黑龍江哈爾濱150030)

        基于熱重法的植物油氧化動力學研究

        張巧智,李 楊*,齊寶坤,隋曉楠,馬文君,江連洲
        (東北農(nóng)業(yè)大學食品學院,黑龍江哈爾濱150030)

        采用熱重(TG)及微商熱重(DTG)法研究了6種市售植物油的氧化特性,探究了不同升溫速率對TG及DTG曲線的影響,通過DTG曲線確定了不同植物油在不同升溫速率下的起始氧化溫度,并利用Ozawa-Flynn-Wall方程及Arrhenius方程對油脂的氧化進行了動力學參數(shù)的推導,同時分析了油脂熱降解特性與其主要脂肪酸組成間的相關性。研究表明:隨著升溫速率的升高,氧化分解的溫度隨之升高,油脂的起始氧化溫度、活化能等參數(shù)均可用于表征油脂的氧化穩(wěn)定性,且植物油的熱降解特性與其脂肪酸組成間有著不同程度的相關性;其中油脂的起始氧化溫度與其飽和脂肪酸、單不飽和脂肪酸及代表脂肪酸含量成正比,與多不飽和脂肪酸含量成反比;油樣的活化能及前指數(shù)因子同樣也與飽和脂肪酸含量成正比且前者相關性顯著,而速率常數(shù)κ與樣品組分相關性與Ton-組分相關性趨勢相反。熱重分析法采樣量少,操作簡單快速,精密度高,可為油脂加工及儲藏條件的控制、新型抗氧化劑研發(fā)等提供實驗依據(jù)。

        植物油,氧化穩(wěn)定性,熱重分析,微商熱重分析,氧化動力學

        油脂在人體膳食結(jié)構中扮演著必不可少的角色,同時也是一種重要的食品工業(yè)原料,然而它在加工、運輸、儲藏的過程中易發(fā)生氧化酸敗,生成的過氧化物及其進一步降解/聚合生成的低分子量醛、酮及環(huán)氧化物、羥基化合物等其他有毒物質(zhì),破壞食品的營養(yǎng)價值和感官特性[1]。食用氧化的油脂會造成機體損傷,加速衰老,誘發(fā)多種疾病,嚴重危害人體健康。因此,油脂酸敗程度是評價油脂品質(zhì)的一項重要指標。近些年來,研究人員提出了一些測定油脂氧化程度的檢測方法,如硫代巴比妥酸法(TBA法)、AOM活性氧法、schall烘箱法、OSI法、Rancimat法、氣液色譜法(GLC)、紫外吸收光譜法(UV)等[2]。這些方法或所需樣品量多,或需使用有毒有害試劑、或操作繁瑣,且均忽略了在氧化過程中油樣重量的變化。熱重法是在程序控溫的情況下通過高度靈敏的電子天平記錄物質(zhì)的質(zhì)量隨溫度或時間變化的一種熱分析技術。利用該方法可研究油脂的氧化穩(wěn)定性,開發(fā)新型抗氧化劑,測定油脂氧化分解溫度、氧化誘導時間等[3-4],并與傳統(tǒng)Rancimat法所測結(jié)果有較好的相關性[2,4-5]。熱重法采樣少、操作簡單方便、測定時間短,在油脂生產(chǎn)加工及質(zhì)量監(jiān)控領域逐漸受到人們的重視[6]。

        油脂在酸敗過程中會發(fā)生一系列氧化、裂變、聚合等反應,重量發(fā)生變化,熱重儀通過對這些重量變化的精確記錄反映出油脂的氧化進程,結(jié)合動力學推導及與組分相關性的研究可進一步闡述油脂氧化特性[7-8],這一方法在油脂理論研究和實際生產(chǎn)中有著廣闊的應用前景[9],目前國外已有相關研究[7,10-12],但國內(nèi)鮮有報道。本文通過測定不同升溫速率下植物油的熱重(TG)及微商熱重(DTG)曲線,比較不同品種植物油的氧化特性,并通過動力學求解及與樣品脂肪酸組成的相關性研究,進一步探討油脂的氧化特性。

        1 材料與方法

        1.1 材料與儀器

        橄欖油 上海佳格食品有限公司蘇州分公司;葡萄籽油 上海勝夯貿(mào)易有限公司;葵花籽油 上海佳格食品有限公司蘇州分公司;玉米油 天津中糧佳悅有限公司;大豆油 九三集團哈爾濱惠康食品有限公司;芝麻油 天津嘉里糧油有限公司;油樣購于當?shù)爻?,生產(chǎn)日期相近,加工精度等級一致,商品標簽上均標明不含有人工添加的抗氧化劑,分析前將油樣置于棕色瓶內(nèi)于低溫(4℃)避光無頂空條件下保存1周;脂肪酸甲酯標準品 購于Sigma公司;乙醚、正己烷、甲醇、氫氧化鉀 均為色譜純。

        STA449 F3型STA同步熱分析儀 德國耐馳公司;Agilent GC7890A-5975C型氣相色譜-質(zhì)譜聯(lián)用儀 美國安捷倫公司。

        1.2 實驗方法

        1.2.1 熱重測定方法 有研究表明當樣品質(zhì)量在1~4 mg范圍內(nèi)時可保證凹盤內(nèi)樣品與吹掃氣體的充分接觸,避免溫度梯度的產(chǎn)生[9],因此,稱?。?±0.5)mg樣品置于密封中間有孔的鋁凹盤中,放入儀器加熱爐。溫度控制程序如下:30℃保溫5 min,分別以5個不同升溫速率(5、7.5、10、12.5、15℃/min)升至400℃。吹掃氣體為O2(99.99%,80 mL/min),天平保護氣體為N2(99.99%,20 mL/min)。

        1.2.2 起始氧化溫度的確定及動力學方程的推導

        在氧化初始階段,熱重曲線變化細微,對其取一階導數(shù)得到相應的微商熱重曲線,表征重量變化速率對溫度的關系[7,10],以微商熱重曲線偏離基線時的溫度為起始氧化溫度(Ton)[8],見圖1。

        在本文的非等溫條件下,O2始終處于過量狀態(tài),因此可忽略熱氧化過程中O2的消耗,即油脂氧化過程與O2的量無關,氧化屬于一級動力學[5,8,11]。將同一樣品在5個不同升溫速率下的Ton帶入Ozawa-Flynn-Wall[12-13]方程中:

        圖1 熱重曲線起始氧化溫度的確定Fig.1 Confirmation of the onset oxidation temperature of TG curve

        方程中β為升溫速率(℃/min),R為摩爾氣體常數(shù)(8.314 J/mol·K),T為Ton(K)。以logβ對,所得直線的斜率和截距可求得活化能Ea(J/mol)和前指數(shù)因子A(min-1),結(jié)合Arrhenius方程可進一步計算得到反應速率κ(min-1):

        1.2.3 脂肪酸組成測定方法 脂肪酸甲酯的制備方法參考GB/T 17377-2008的方法,色譜柱溫控程序為:進樣口溫度和色譜-質(zhì)譜接口溫度分別為260℃和280℃;初始溫度70℃保持1 min,以5℃/min升至100℃,保持2 min,再以10℃/min升至175℃,保持40 min,最后以15℃/min升至220℃,保持30 min。通過對照標準品保留時間對樣品中脂肪酸進行定性,面積歸一化定量,確定每種油樣中脂肪酸的種類和含量(%)。

        1.2.4 數(shù)據(jù)分析 實驗均平行測定三次,對實驗數(shù)據(jù)的平均值、標準差、差異顯著性分析、熱氧化特性與脂肪酸組成的皮爾森相關性分析均采用SPSS 22.0統(tǒng)計分析軟件。差異顯著性分析采用Duncan多重檢驗,95%置信水平(p<0.05)。動力學直線擬合采用OriginPro 8.5.1。

        2 結(jié)果與討論

        2.1 不同升溫速率對熱重及微商熱重曲線的影響

        按照1.2.1中敘述的方法,分別測定6個油脂樣品在5個升溫速率(5、7.5、10、12.5、15℃/min)下的熱重曲線見圖2,為了更清晰的對比曲線間的差異,各圖內(nèi)從上至下TG曲線按升溫速率遞增的順序排列。

        由圖2可知,不同升溫速率對6種植物油的TG曲線有顯著影響。在熱氧化初期油脂中脂肪酸鏈吸收氧氣,重量會有所增加,體現(xiàn)在熱重曲線上則為增重峰,這一“增重”現(xiàn)象在低加熱速率時較明顯,而在高加熱速率(如15℃/min)時變得平坦而幾乎消失,這是因為在較低升溫速率的情況下,溫度緩慢上升,樣品與氧氣充分接觸反應后才開始分解反應,增重明顯;而在較高升溫速率時,坩堝內(nèi)熱氧化反應劇烈,樣品未能充分進行氧化就開始分解[1-2]。同時,隨著升溫速率的升高,增重峰向較高溫度方向移動,由TG和DTG曲線確定的起始氧化溫度(Ton)也隨之增加(表1),如在5℃/min升溫速率下,大豆油的Ton為173.55℃,而當升溫速率為15℃/min時,其Ton升至194.58℃。這些現(xiàn)象都說明升溫速率越快,發(fā)生氧化分解的溫度就越高。其他樣品的熱重及微商熱重曲線也呈現(xiàn)類似趨勢。

        2.2 熱重曲線的動力學推導

        按照1.2.2的方法,由樣品在不同升溫速率下的熱重及微商熱重曲線確定的起始氧化溫度(Ton)見表1。根據(jù)已有文獻報道,Ton表征油脂對熱氧化分解的耐受程度[6],因此,在相同實驗條件下,起始氧化溫度的高低可以比較各種油脂的氧化穩(wěn)定性,從表1可知,不論在哪一升溫速率下,6種植物油的氧化穩(wěn)定性強弱順序均為:葡萄籽油<葵花籽油≈玉米油≈大豆油<橄欖油<芝麻油,其中,芝麻油雖然含有較多的不飽和脂肪酸(表3),但其具有最高的Ton,穩(wěn)定性最強,可能是由于芝麻油中含有活性較強的天然抗氧化劑—芝麻素和芝麻酚[14],而葡萄籽油氧化穩(wěn)定性最低的原因可能是其多不飽和脂肪酸含量較高(表3)。

        圖2 6種植物油在不同升溫速率下的熱重曲線Fig.2 TG curves of 6 vegetable oils under different heating rates

        將表1中的Ton數(shù)值帶入Ozawa-Flynn-Wall方程中,以升溫速率logβ對起始氧化溫度1/Ton作圖,圖3以芝麻油為例表示其logβ與1/Ton間的線性關系。由圖3中可知,logβ與1/Ton間具有很好的線性關系,可用線性方程y=-4631x+10.89表示,相關系數(shù)為0.9978。其他樣品的線性擬合相關系數(shù)也均大于0.89(表2),利用方程(1)~(4),求得的動力學參數(shù)見表2。通過動力學求解,6種油脂樣品的活化能分別為:橄欖油,74.13 kJ/mol;葡萄籽油,68.03 kJ/mol;葵花籽油,73.66 kJ/mol;玉米油,72.47 kJ/mol;大豆油,87.45 kJ/ mol;芝麻油,84.28 kJ/mol。所得的活化能為油脂樣品進行復雜熱氧化反應的表觀活化能[6],表征樣品進行氧化劣變的難易程度。以活化能的高低對油脂樣品的氧化穩(wěn)定性排序為:葡萄籽油<玉米油<葵花籽油<橄欖油<芝麻油<大豆油,與以Ton高低排序略有差異,說明單一指標表征并不完全,要結(jié)合多個指標從多個方面綜合評價油脂的氧化穩(wěn)定性。

        圖3 芝麻油起始氧化溫度與升溫速率間的線性關系Fig.3 The linear relationship between onset oxidation temperature and heating rate of sesame oil

        2.3 樣品熱氧化特性與其脂肪酸組成的相關性

        對油脂樣品的脂肪酸分析結(jié)果見表3,按不飽和程度可將脂肪酸分為飽和脂肪酸(saturated fatty acid, SFA),單不飽和脂肪酸(monounturated fatty acid,MUFA)和多不飽和脂肪酸(polyunturated fatty acid,PUFA)。由表3可知,6種油脂樣品中均主要含有四種脂肪酸,即軟脂酸(C16∶0),棕櫚酸(C18∶0),油酸(C18∶1)和亞油酸(C18∶2)。對油樣組分(上述四種脂肪酸含量,SFA,MUFA,PUFA含量)與其熱氧化過程中的參數(shù)(起始氧化溫度Ton,活化能Ea,前指數(shù)因子A及速率常數(shù)κ)進行皮爾森相關性分析,相關系數(shù)的大小和正負代表不同組分對油脂氧化穩(wěn)定性的影響,具體結(jié)果見表4。

        表1 不同升溫速率下植物油樣品的起始氧化溫度(℃)Table1 Onset oxidation temperatures of oil samples under different heating rates(℃)

        表2 植物油樣品的動力學參數(shù)Table2 Kinetic parameters of oil samples

        表3 植物油樣品的脂肪酸組成(%)Table3 Fatty acid compositions of oil samples(%)

        由表4可知,在油脂組分與其起始氧化溫度及動力學參數(shù)間有著不同程度的相關性。其中,油脂的起始氧化溫度Ton與飽和脂肪酸含量及軟脂酸(C16∶0)含量成正比,同時也與單不飽和脂肪酸及油酸含量(C18∶1)成正比,但與多不飽和脂肪酸及亞油酸(C18∶2)含量成反比。這也解釋了葡萄籽油具有較低Ton的原因是其含有較多的PUFA和亞油酸,證明了前面的假設。油樣的活化能Ea同樣也與其飽和脂肪酸含量及硬脂酸含量(C18∶0)成正比且與前者相關性顯著(0.05水平),因此可由油脂的SFA含量推斷其活化能的高低,但Ea與其他組分的相關性較弱,其中緣由還需進一步的探究。由Ea推得的前指數(shù)因子A同樣也與飽和脂肪酸含量成正比。200℃時的速率常數(shù)κ與樣品組分含量的相關性與Ton-組分相關性有相反趨勢。這些數(shù)據(jù)可進一步明確植物油的熱降解特性及影響油脂熱穩(wěn)定性的內(nèi)部因素,對“熱特性-組分”相關性的研究可為油脂儲藏、加工條件控制及質(zhì)量監(jiān)控提供一定的實驗依據(jù)。此外,本研究所用植物油樣中含有天然抗氧化劑,其中主要為酚類物質(zhì)(多酚、生育酚等),Vecchio等通過對單品種初榨橄欖油的組分-熱特性相關性的研究,發(fā)現(xiàn)油樣的熱降解參數(shù)與其總酚含量并無顯著相關,原因可能是這些物質(zhì)在200℃前已發(fā)生降解,對后續(xù)的熱曲線影響微弱[7,15]。從樣品脂肪酸組成出發(fā)研究其與氧化穩(wěn)定性的相關性是“組分-熱降解特性”關系的初步探索,接下來對組分的進一步細致化研究如對甘三酯組成、微量氧化產(chǎn)物等與熱特性的關系均有待于進一步探索。

        表4 植物油樣品熱氧化特性與其脂肪酸的相關性分析Table4 Correlation analysis between thermal oxidation properties and fatty acid compositions of oil samples

        3 結(jié)論

        熱重及微商熱重分析可用于研究油脂的氧化特性,本文通過探究升溫速率對熱重曲線的影響,發(fā)現(xiàn)隨著升溫速率的升高,起始氧化溫度和“增重峰”均向高溫方向移動,表明氧化分解的溫度隨之升高。利用OFW方程及阿倫尼烏斯方程對油脂樣品進行了動力學參數(shù)的求解,發(fā)現(xiàn)起始氧化溫度、表觀活化能、前指數(shù)因子及速率常數(shù)均可表征油脂的氧化穩(wěn)定性及熱降解反應的特性。同時,通過對油脂熱降解特性與其主要脂肪酸組成間的相關性分析,發(fā)現(xiàn)油脂的Ton,Ea,A,κ與SFA,MUFA,PUFA及代表脂肪酸含量間有著不同程度的相關性。其中,油脂的Ton與其飽和脂肪酸、單不飽和脂肪酸及代表脂肪酸含量成正比,與多不飽和脂肪酸及代表脂肪酸含量成反比;油樣的活化能及前指數(shù)因子同樣也與飽和脂肪酸含量成正比且前者相關性顯著,而速率常數(shù)κ與樣品組分相關性與Ton-組分相關性趨勢相反。關于熱特性-組分相關性的研究需進一步擴大油脂樣品范圍,通過主成分分析等統(tǒng)計方法進行進一步的探索。熱重分析法具有采樣量少,操作簡單,測定時間短等優(yōu)點,可為油脂加工、儲藏條件的控制、抗氧化劑研發(fā)提供實驗依據(jù),在油脂生產(chǎn)加工及質(zhì)量監(jiān)控等領域具有廣闊的應用前景。

        [1]E Choe,D B Min.Mechanisms and factors for edible oil oxidation[J].Comprehensive Reviews in Food Science and Food Safety,2006,5(4):169-186.

        [2]B Barriuso,I Astiasara'n,D Ansorena.A review of analytical methods measuring lipid oxidation status in foods:a challenging task[J].European Food Research and Technology,2013,236 (1):1-15.

        [3]張振山,張麗霞,李棟,等.熱重法對亞麻籽油的氧化動力學研究[J].食品工業(yè)科技,2013(11):62-64,68.

        [4]楊國燕.DSC和Rancimat法測定亞麻籽油氧化穩(wěn)定性研究[J].糧食與油脂,2014,27(8):29-32.

        [5]E O Ligeza,W Bekas,D Kowalska,et al.Kinetics of commercial olive oil oxidation:Dynamic differential scanning calorimetry and Rancimat studies[J].European Journal of Lipid Science and Technology,2010,112(2):268-274.

        [6]顧敏芬,王昉,金宜英.熱重分析法在油脂氧化穩(wěn)定性研究中的應用[J].中國油脂,2012,37(9):85-87.

        [7]S Vecchio,L Cerretani,A Bendinib,et al.Thermal decomposition study of monovarietal extra virgin olive oil by simultaneous thermogravimetry/differential scanning calorimetry:Relation with chemical composition[J].Journal of Agricultural and Food Chemistry,2009,57(11):4793-4800.

        [8]S I M Monteagudo,M D A Saldan,J J Kennelly.Kinetics of non-isothermal oxidation of anhydrous milk fat rich in conjugated linoleic acid using differential scanning calorimetry[J].Journal of Thermal Analysis and Calorimetry,2012,107(3):973-981.

        [9]W Tsang,J A Walker.Characterization of RDF properties through high-pressure differential scanning calorimetry[J].Resources and Conservation,1982,9(1-4):355-365.

        [10]S Vecchio,L Campanella,A Nuccilli,et al.Kinetic study of thermal breakdown of triglycerides contained in extra-virgin olive oil[J].Journal of Thermal Analysis and Calorimetry,2008,91(1):51-56.

        [11]A Adhvaryua,S Z Erhana,Z S Liu,et al.Oxidation kinetic studies of oils derived from unmodified and genetically modified vegetables using pressurized differential scanning calorimetry and nuclear magnetic resonance spectroscopy[J].Thermochimica Acta,2000,364(1-2):87-97.

        [12]J H Flynn,L A Wall.A quick direct method for determination of activation energy from thermogravimetric data[J].Journal of Polymer Science Part B-Polymer Letters,1966,4(5PB):323-328.

        [13]T Ozawa.A new method of analyzing thermogravimetric data [J].Bulletin of the Chemical Society of Japan,1965,38(11):1881-1886.

        [14]劉菁,侯利霞,劉玉蘭,等.芝麻及芝麻油中主要抗氧化物質(zhì)的研究進展[J].中國調(diào)味品,2013,38(10):16-19(45).

        [15]A C Pancorbo,L Cerretani,A Bendini,et al.Evaluation of the influence of thermal oxidation in the phenolic composition and in antioxidant activity of extra-virgin olive oils.Journal of Agricultural and Food Chemistry,2007,55:4771-4780.

        Study on the oxidation kinetics of vegetable oils by thermogravimetry

        ZHANG Qiao-zhi,LI Yang*,QI Bao-kun,SUI Xiao-nan,MA Wen-jun,JIANG Lian-zhou
        (College of Food Science,Northeast Agricultural University,Harbin 150030,China)

        In this paper,thermogravimetry(TG)and derivative thermogravimetry(DTG)methods were applied to study the oxidative properties of 6 commercial vegetable oils.The effect of heating rate on the TG and DTG curves was investigated and the onset oxidation temperatures(Ton)of oils under different heating rates were confirmed using DTG curves.Through applying Ozawa-Flynn-Wall and Arrhenius equation,the kinetic parameters of oil oxidation were calculated.A correlation analysis between oil degradation properties and their major fatty acids compositions was also studied.The results showed that the oxidation decomposition temperature increased as the heating rate went up.The onset oxidation temperature,activation energy etc.could be used to feature the oxidation stability of vegetable oils.Correlations at different levels were found between thermal-degradation properties of oils and their fatty acid compositions:the onset oxidation temperatures of oil were found to be positive correlated with their amount of saturated fatty acid(SFA),monounsaturated fatty acid(MUFA)as well as typical FAs,while negative correlation was found between Tonand polyunsaturated fatty acid(PUFA).As for Ea and A,both were positive associated with SFA content and the former was statistically significant.With regard to rate content,contrary trend was found compared with Toncomposition relationship.TG analysis was known as an easy-operating,time-saving method with less sample mass and high precision.It could be used to offer experimental bases for the condition control of oil processing and storage as well as the exploration of new antioxidants.

        vegetable oil;oxidation stability;thermogravimetry;derivative thermogravimetry;oxidative kinetics

        TS221

        A

        1002-0306(2016)02-0103-05

        10.13386/j.issn1002-0306.2016.02.012

        2015-07-17

        張巧智(1990-),女,碩士研究生,研究方向:糧食、油脂及植物蛋白工程,E-mail:miumiulovelife@163.com。

        *通訊作者:李楊(1981-),男,博士,副教授,研究方向:糧食、油脂及植物蛋白工程,E-mail:liyanghuangyu@163.com。

        國家自然基金重點項目(3143000560);農(nóng)業(yè)部崗位科學家(CARS-04-PS25)。

        猜你喜歡
        植物油油脂組分
        《中國油脂》征訂啟事
        組分分發(fā)管理系統(tǒng)在天然氣計量的應用
        QuEChERS-氣相色譜-質(zhì)譜法測定植物油中16種鄰苯二甲酸酯
        一種難溶難熔未知組分板材的定性分析
        A taste of Peking duck
        中國油脂
        黑順片不同組分對正常小鼠的急性毒性
        中成藥(2018年8期)2018-08-29 01:28:26
        金雀花中黃酮苷類組分鑒定及2種成分測定
        中成藥(2018年2期)2018-05-09 07:20:09
        歡迎訂閱2019年《中國油脂》
        中國食用植物油供需平衡表
        私人毛片免费高清影视院| 亚洲激情一区二区三区视频| 中文字幕一区二区综合| 五月丁香六月综合缴清无码| 国内老熟妇对白xxxxhd| 91超碰在线观看免费| 在线观看高清视频一区二区三区| 免费av片在线观看网址| 国产午夜精品理论片| 国产激情久久99久久| 日韩性感av一区二区三区| 天堂视频在线观看一二区| 欧美裸体xxxx极品少妇| 色综合另类小说图片区| 国内国外日产一区二区| 亚洲午夜久久久精品影院| 中文字幕av无码一区二区三区| 色噜噜狠狠色综合中文字幕| 久久开心婷婷综合中文| 国产精品亚洲片在线观看不卡| 亚洲依依成人亚洲社区| 欧美久久中文字幕| 国产一区二区黑丝美胸| 亚洲国产精品无码专区在线观看| 丰满多毛少妇做爰视频| 超级少妇一区二区三区| 亚洲一区二区三区地址| 狠狠综合久久av一区二区| 国产精品久久久久久2021| 中文字幕日韩精品亚洲精品| 人妻体内射精一区二区三区| 少妇寂寞难耐被黑人中出 | 亚洲国产av剧一区二区三区| 麻神在线观看免费观看| 国产裸体xxxx视频在线播放| 欧美亚洲国产人妖系列视| 一边摸一边做爽的视频17国产| 一色桃子中文字幕人妻熟女作品| 日韩精品无码一区二区三区免费| 日韩熟女一区二区三区| 日本一区二区三区视频在线观看|