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

基于低場核磁共振技術和傅里葉變換紅外光譜的植物油種類判別可行性分析

2025-08-07 00:00:00向燕陸佳怡何強鄧莎董怡
中國調味品 2025年7期
關鍵詞:玉米油葵花籽橄欖油

中圖分類號:TS225.1 文獻標志碼:A 文章編號:1000-9973(2025)07-0191-07

DOI:10.3969/j.issn.1000-9973.2025.07.028

Feasibility Analysis of Discrimination of Plant Oil Types Based on Low-Field Nuclear Magnetic Resonance Technology and Fourier Transform Infrared Spectroscopy

XIANG Yan,LU Jia-yi,HE Qiang,DENG Sha,DONG Yi* (College of Biomass Science and Engineering,Sichuan University,Chengdu 6lOo65,China)

Abstract: The adulteration phenomenon of plant oils often occurs, so the use of non-destructive and eficient methods to identifythe types of plant oils is an important direction for the technological development in the oil industry. In this paper,with five types of plant oils (olive oil,peanut oil,corn oil,sunflower seed oil and mixed oil)as the research objects,the feasibility of using low-field nuclear magnetic resonance (LF-NMR) technology and Fourier transform infrared spectroscopy(FTIR) to distinguish the types of plant oils is analyzed.Hierarchical cluster analysis (HCA) is conducted on faty acid components. The results show that five types of plant oils can be divided into three groups: corn oil and sunflower seed oil, peanut oil and mixed oil,and olive oil. LF-NMR combined with principal component analysis (PCA) could identify the linoleic acid-type plant oils (sunflower seed oil and corn oil)and theoleic acid-type plant oils (olive oil and peanut oil). FTIR combined with PCA could further identify the types of linoleic acid-type plant oils. In conclusion,it is feasible to use LF-NMR and FTIR to identify the types of plant oils, and the combination of the two methods has complementary effects. This study can lay a foundation for the rapid identification of edible oil adulteration.

Key words: plant oil; without destructiveness;identification of types;low-field nuclear magnetic resonance; Fourier transform infrared spectroscopy;principal component analysis

食用油脂是人類日常膳食中重要的組成部分,植物油是食用油脂的主要來源,全球食用植物油的產量和消費量均在逐年增長[1]。然而,隨著植物油種類的增多,由于原料、工藝等原因,不同植物油的產量、成分和價格各不相同,使得在利益驅使下的植物油摻偽現象普遍存在,植物油的品質和安全面臨巨大挑戰[2]。傳統的油脂摻偽檢測主要采用理化檢驗方法,例如高效液相色譜法、氣相色譜法、薄層層析法等[3]。Gong等[4]采用液相色譜-質譜技術(LC-MS)分離食用油中的三酰甘油(TAG),鑒別出在植物油中摻偽的豬油。韓玉澤等[5]利用氣相色譜-質譜聯用技術(GC-MS)測定亞麻籽油的揮發性組成并構建指紋圖譜,建立的模型適用于亞麻籽油的摻偽鑒定。這些方法主要以油脂的內部特征組成為依據進行鑒別,但僅依靠某一特征指標難以獲取整個體系的準確信息[],且這些方法通常需要破壞油樣,操作較復雜,時間和試劑耗費較多,不能達到現代油脂行業的無損、實時、高效的檢測要求。尤其是目前缺少簡單有效的油脂鑒偽方法,不利于遏制市場上的植物油摻偽現象。因此,使用快捷、高效的檢測技術和儀器設備,實現對植物油種類的快速精準鑒別,是油脂行業技術發展的重要方向[7]

低場核磁共振(low-field nuclearmagneticresonance,LF-NMR)技術是一種新興的快速無損檢測技術,廣泛應用于材料、化學、食品等領域[8]。林曉浪等9基于LF-NMR開發的模型可用于油茶籽油的摻偽鑒別。Jin等[10]利用LF-NMR結合化學計量學分析,可快速篩選地溝油。成實等[1]和Paula等[12]研究發現,隨著脂肪酸碳鏈長度的增加,弛豫時間會縮短,而隨著脂肪酸不飽和度的增加,弛豫時間延長。傅里葉變換紅外光譜(Fourier transform inforared spectroscopy,FTIR)是基于對干涉后的紅外光做傅里葉變換而得到的光譜。Usman 等[13]研究表明FTIR對黃油中摻人植物油具有良好的預測效率。Mahboubifar等[14]發現不同種類油脂的紅外光譜差異與脂肪酸鏈長、不飽和度和酰基形式有關。以上研究結果均為LF-NMR和FTIR鑒別植物油種類提供了理論基礎。

現有的研究主要是通過單一檢測方法鑒別植物油的摻偽情況,但植物油體系復雜,不同種類植物油的組成接近,僅用一種方法進行規律分析難以對植物油的摻偽鑒別起到全面完整的指導作用,故本實驗結合橄欖油、花生油、玉米油、葵花籽油及其混合油的內部脂肪酸組成,采用LF-NMR和FTIR兩種方法進行相互驗證,對植物油的種類進行鑒別研究,旨在快速鑒別植物油的種類,進而更好地為食用油的質量安全提供理論指導與技術支持。

1材料與方法

1.1材料

特級橄欖油:西班牙伯爵Borges公司;一級花生油:邯鄲中好糧油加工有限公司;一級葵花籽油:南海油脂工業(赤灣)有限公司;一級玉米油:長壽花食品股份有限公司;混合油:4種植物油按等質量比混合均勻制成。

1.2試劑

37種混合脂肪酸甲酯標準品:上海安譜璀世標準技術服務有限公司;異辛烷(色譜純)、甲醇(色譜純)、氫氧化鉀(分析純)、硫氫化鈉(分析純):市科隆化學品有限公司。

1.3 儀器與設備

QP-2010SE氣相色譜-質譜聯用儀、IRTracer-100傅里葉變換紅外光譜儀日本島津公司;NMI20-040H-I核磁共振成像分析儀蘇州紐邁分析儀器股份有限公司。

1.4方法

1.4.1 脂肪酸組成分析

油樣甲酯化:取 50mg 油樣溶解在 4mL 異辛烷中,充分溶解后加入 200μL 氫氧化鉀-甲醇溶液 2mol/L) ,充分振蕩30s使溶液完全溶解,最后加入約 1g 硫酸氫鈉。取 1~1.5mL 上清液至樣品瓶中,采用GC-MS分析油樣的脂肪酸組成[15]。脂肪酸組成按GB5009.168—2016進行測定,脂肪酸組成的絕對含量按照外標法獲得,相對含量按照面積歸一化法獲得。

GC-MS條件:FFAP石英毛細管色譜柱( 30m× 0.32mm,0.5μm) ;升溫程序:初始溫度 130°C ,保持5min ,以 4°C/min 升至 240°C ,保持 20min ;進樣口溫度: 250°C ;質量掃描范圍 (m/z) : 25~550u ;進樣量1μL 。

1.4.2低場核磁回波衰減曲線測定及反演

量取 7mL 樣品放入低場核磁專用試管內,經37°C 恒溫水浴 15min 后置于LF-NMR樣品槽中,利用Carr-Purcell-Meiboom-Gill脈沖序列測定樣品的橫向弛豫時間 T2

儀器參數設置:采樣頻率 200kHz ;循環延時2s;半回波時間 0.5ms ;增益 50dB ;掃描次數8次;回波個數5 000個[16]。每個樣品進行多次采樣。

樣品的核磁共振衰減信號數據通過T-invfit軟件進行SIRT反演擬合,得到多組分分析 T2 弛豫圖譜。將反演擬合所得 T20 起始時間、 T20 頂點時間、 T20 結束時間 .S20 峰面積、 T20 峰比例、 T21 起始時間、 T21 頂點時間、T21 結束時間、 S21 峰面積、 T21 峰比例、 T22 起始時間、T22 頂點時間、 T22 結束時間、 S22 峰面積、 T22 峰比例、S2 峰面積總和預設為16個變量,進行主成分分析(PCA)。

1.4.3 FTIR采集

用移液槍將 2μL 油樣滴到KBr窗片上,形成均勻的液膜進行測定。采集空氣光譜作為背景光譜,以消除 CO2 和濕度的影響。數據的輸出格式為吸光度。采集范圍: 4000~400cm-1 ,每個光譜掃描32次,分辨率為 4cm-1[17] 。每個樣品的光譜收集多次并取平均值,對光譜數據進行PCA。

1.5 數據處理

所有實驗設置3個平行,結果表示為平均值士標準差。運用SPSS26進行系統聚類分析(HCA)和PCA,并通過單因素ANOVA檢驗對數據進行差異顯著性分析, Plt;0. 05 表示差異顯著。利用Excel處理數據,使用Origin2022繪圖。

2 結果與分析

2.1不同種類植物油的脂肪酸組成分析

脂肪酸是油脂中重要的組成成分,脂肪酸的不飽和度、種類、碳鏈長度等均會影響油脂的性質,是衡量油脂品質的重要指標[18]。4種植物油及其混合油的主要脂肪酸組成見表1。5種植物油的主要脂肪酸組成基本一致,共檢測出14種主要脂肪酸,其中棕櫚酸( (Cl6:0) 、棕櫚油酸( Cl6:1) 、油酸 (Cl8:1) 、亞油酸( Cl8:2) 和亞麻酸( Cl8:3) 的總相對含量均達到 70% 以上,其中橄欖油的該5種脂肪酸占比最高,達到 95.35% ,而花生油的占比最低,為 73.73% 。

表1不同種類植物油的脂肪酸組成Table1 Fatty acid composition ofdifferent types of plant oils
注:單不飽和脂肪酸包括棕櫚油酸、十七碳一烯酸和油酸;多不飽和脂肪酸包括亞油酸和亞麻酸;同行相同指標的不同小寫字母表示差異顯著 (Plt;0.05) ,下表同。

由于品種、種植環境、壓榨方式等因素不同,不同植物油的脂肪酸含量具有差異性。葵花籽油和玉米油的亞油酸 (Cl8:2) 絕對含量較高,兩者相對含量均大于50% ,且兩者的多不飽和脂肪酸(polyunsaturatedfattyacid,PUFA)的相對含量與絕對含量均高于飽和脂肪酸(saturatedfattyacid,SFA和單不飽和脂肪酸(monounsaturatedfattyacid,MUFA),屬于亞油酸型植物油。橄欖油和花生油的油酸( Cl8:1) 相對含量與絕對含量較高,且MUFA的相對含量與絕對含量高于SFA和PUFA,屬于油酸型植物油。而混合油的MUFA和PUFA相對含量與絕對含量差異不顯著。相較于其他植物油,橄欖油的棕櫚酸( Cl6:0. 、棕櫚油酸( C16:1) 和油酸1 Cl8:1) 相對含量更高 (82.21% ,且其MUFA相對含量高達 66.26% 。花生油的長鏈脂肪酸的相對含量顯著高于其他植物油 (Plt;0.05) ,如花生酸 (C20:0) )、山崙酸 (C22:0). 和木焦油酸 (24:0) 的相對含量均超過4% 。此外,葵花籽油和玉米油的PUFA含量顯著高于花生油和橄欖油 ?Plt;0.05) ,而其MUFA含量顯著低于花生油和橄欖油 (Plt;0.05) 。混合油的各脂肪酸含量基本介于亞油酸型植物油(葵花籽油和玉米油)與油酸型植物油(橄欖油和花生油)之間。

HCA可更直觀地表現出5種植物油在脂肪酸組成上的區別[19],結果見圖1。

圖1基于脂肪酸組成的系統聚類分析圖 Fig.1 Hierarchical clusteranalysis diagram based on fatty acid composition

根據脂肪酸組成的相似程度,可以將這5種植物油分為3組,由圖1可知,葵花籽油和玉米油的脂肪酸組成相似度較高,混合油的脂肪酸組成與花生油相似,橄欖油與其他植物油的脂肪酸組成差異顯著。

2.2 LF-NMR結果分析

2.2.1不同種類植物油的弛豫圖譜分析

反演擬合核磁共振衰減信號得到橫向弛豫圖譜,以分析不同種類植物油的LF-NMR弛豫特性,結果見圖2和表2。

圖2不同種類植物油的橫向弛豫圖譜Fig.2 Transverse relaxation spectra ofdifferent types of plant oils

由圖2可知,5種植物油的多組分弛豫圖譜相似,主要存在 T21 響應峰,這是由于植物油的LF-NMR弛豫特性與其主要脂肪酸種類和脂肪酸含量密切相關[20]。

表2不同種類植物油的LF-NMR弛豫特性指標結果Table2 Results of indexes of LF-NMR relaxation propertiesof different types of plant oils

由表2可知,以峰面積 S21 為評價標準時,按峰面積 S21 從小到大排序為葵花籽油 lt; 玉米油 lt; 混合油 lt; 花生油 lt; 橄欖油。當以 T21 峰頂點時間為評價標準時,按時間從早到晚排序為橄欖油、花生油、混合油、玉米油、葵花籽油,與峰面積 S21 排列順序剛好相反。研究表明,內部結構和外部環境均會影響脂質的弛豫特性,如油酸甲酯和油酸間形成的不同液晶狀態,可改變峰面積和弛豫時間[12]。不同種類植物油的脂肪酸組成和含量存在差異,從而引起多組分弛豫譜圖中弛豫時間和峰面積發生變化[21]。

2.2.2 PCA結果

PCA是一種常用的數據降維方法,可最大限度地保留并發掘原始數據的特征[22]。對5種植物油的低場核磁回波衰減數據進行PCA,結果見表3和圖3。

表3LF-NMR弛豫特性在主成分中的載荷和累計方差貢獻率 Table3Loadings and cumulative variance contribution rates ofLF-NMR relaxation properties in principal components

由表3可知,主成分1(PC1)的貢獻率為 46.54% 結合主成分2(PC2)、主成分3(PC3)和主成分4(PC4)的累計貢獻率可達到 91.77% ,保留了原始數據的大部分信息。以得分值繪制5種植物油的得分圖,結果見圖3。

圖3低場核磁回波衰減數據的PCA得分圖Fig.3 PCA score plot oflow-field nuclearmagnetic echo decay data

由圖3可知,橄欖油和花生油在PCA得分圖上的分布較接近,玉米油和葵花籽油的分布有部分重疊,橄欖油和花生油與玉米油和葵花籽油存在明顯的邊界。混合油介于中間,與其他4種油均有部分重疊。由此說明,LF-NMR的PCA能較好地鑒別出亞油酸型植物油(葵花籽油和玉米油)和油酸型植物油(橄欖油和花生油),這與彭丹等[6]的研究結果一致。

2.3 FTIR結果分析

2.3.1不同種類植物油的紅外吸收光譜分析

不同種類植物油在 4000~400cm-1 之間的紅外吸收光譜見圖4。

圖4不同種類植物油的紅外吸收光譜 Fig.4 Infrared absorption spectra of different types of plant oils

由圖4可知,5種植物油的紅外吸收光譜基本一致,在 3 008cm-1 處有不飽和碳的C—H伸縮振動峰,在 2 925cm-1 和 2854cm-1 處分別為飽和碳的C—H反對稱伸縮振動峰和對稱伸縮振動峰[23],在1 747cm-1 處有 C=O 鍵的伸縮振動峰,在 1465cm-1 處是亞甲基的變形振動,在 1163cm-1 處有酯類的C—O伸縮振動峰,在 723cm-1 處有碳鏈骨架振動峰[24]。5種植物油的紅外吸收光譜相似,即特征峰的數目、形狀和位置基本相同,說明其主要成分的種類一致。

紅外吸收光譜中的特征吸收峰強度以吸光度表示,以各植物油的 C=O 伸縮振動峰 (1747cm-1 處)強度為基準,計算其他特征吸收峰的相對強度,結果見表4。

表4紅外吸收光譜的特征吸收峰強度 Table 4 Intensity of characteristic absorption peaks in infrared absorption spectra
由表4可知,5種植物油的特征吸收峰強度與相對強

度均存在差異。 3008cm-1 處是不飽和碳的C—H伸縮振動峰,此處的峰相對強度由小到大排序為花生油 lt; 橄欖油 lt; 混合油 lt; 葵花籽油 lt; 玉米油,說明油樣中不飽和脂肪酸的含量依次增多[24]。 2925cm-1 處是飽和碳的C—H伸縮振動峰,此處的峰相對強度從小到大排序為玉米油 lt; 葵花籽油 lt; 混合油 lt; 橄欖油 lt; 花生油,說明油樣中的飽和脂肪酸含量依次增多,與不飽和脂肪酸含量結果相反。 1465cm-1 處是亞甲基的彎曲振動峰,723cm-1 處是碳鏈骨架振動峰,相較于其他植物油,花生油在這兩處的強度最強,說明花生油中亞甲基含量較高,即脂肪酸的碳鏈較長,這與前文中脂肪酸分析結果一致。 1163cm-1 處是甘油三酯中的C—〇伸縮振動峰,玉米油和葵花籽油在此處的峰相對強度更高,說明這兩種植物油的甘油三酯含量更高[25]。混合油各特征峰的相對強度始終介于葵花籽油和橄欖油之間,各成分含量也介于兩者之間,這與GC-MS結果和LF-NMR結果基本符合。

2.3.2 PCA結果

基于FTIR的差異,同時為了消除首位噪聲、避免信息冗余[26],截取 1800~720cm-1 范圍內的光譜數據,采用PCA對5種植物油的FTIR數據進行分析,結果見圖5。

圖5紅外吸收光譜數據的PCA得分圖Fig.5 PCA score plot of infrared absorption spectra data

由圖5可知,主成分(PC1)的貢獻率為 91.43% 主成分(PC2)的貢獻率為 8.05% ,兩個主成分的累計貢獻率達 99.48% ,能夠保留大部分的原始紅外吸收光譜信息。玉米油、葵花籽油與其他植物油均有明顯界限,能較好地通過FTIR的PCA鑒別出來。橄欖油和花生油有部分重疊,可能是由于脂肪酸組成的種類和含量較接近。混合油與其他植物油有較明顯的邊界,說明混合后油的紅外光譜在一定程度上發生改變。此結果與LF-NMR結果基本符合,并能進一步鑒別出亞油酸型植物油的種類,與LF-NMR聯用能更好地對油脂的種類進行無損鑒別。

3結論

本研究基于植物油中脂肪酸組成,分析通過LF-NMR和FTIR兩種方法鑒別植物油種類的可行性。利用GC-MS對5種植物油進行脂肪酸組成分析,并根據脂肪酸含量進行了系統聚類分析。結果表明,葵花籽油和玉米油的相似度較高,混合油的脂肪酸組成與花生油相似,橄欖油與其他植物油差異顯著。采用LF-NMR回波衰減曲線信息結合PCA,結果表明,由于5種植物油的主要脂肪酸種類一致,其橫向弛豫圖譜相似,主要存在 T21 響應峰;PCA回波衰減曲線數據能較好地鑒別出油酸型油脂(橄欖油和花生油)和亞油酸型油脂(葵花籽油和玉米油)。采用FTIR結合PCA,結果表明,5種植物油的紅外吸收光譜相似,其特征峰強度有差異;對1800~720cm-1 范圍內的紅外光譜數據進行PCA能進一步鑒別出亞油酸型植物油的種類。綜上,采用LF-NMR與FTIR技術鑒別植物油種類可行,且兩種方法聯用具有互補效應,能更好地對植物油的種類進行無損鑒別,為市場快速鑒別植物油種類、篩查摻偽植物油的快速鑒別提供了理論與方法基礎。

參考文獻:

[1朱銳,劉永輝,梁荊璞,等.基于智能手機的食用油快速檢測試紙分析平臺的構建[J].食品工業,2024,45(2):267-270.

[2]王建化,宋瑜,王彩惠,等.芝麻油的摻偽檢驗J」中國調味品,2018,43(11):137-140.

[3]閆軍劍,李潔艷,高健,等.我國食用油安全狀況調查[J].糧油與飼料科技,2020(6):18-20.

[4]GONGC,LUXT,ZHANGSD,etal.Detection oflardadulterationin3kinds of vegetable oilsby liquid chromatography-mass spectrometry with porous graphite carbon column[J].Analytical Sciences,2024,40(7):1289-1299.

[5]韓玉澤,王興瑞,李應霞,等.基于SPME-GC-MS的青海亞麻籽油揮發性組分指紋圖譜構建及摻偽識別[J].現代食品科技,2021,37(11):364-371.

[6彭丹,史翠熠,陳名揚,等.基于低場核磁共振技術的植物油種類快速識別[J].食品科學,2023,44(16):340-346.

[7]俞玥,宋嘉慧,魯玉杰,等.食用植物油摻偽鑒別技術研究進展[J].食品安全質量檢測學報,2021,12(13):5153-5161.

[8]王明洋.糜子自然發酵品質及其淀粉特性研究[D].哈爾濱:東北農業大學,2024.

[9]林曉浪,傅利斌,王欣.基于低場核磁弛豫特性的油茶籽油支持向量機摻偽鑒別模型的建立與評價[J].食品科學,2024,45(10):19-27.

[10]JINHQ,TULY,WANGYX,etal.Rapid detection ofwaste cooking oil using low-field nuclear magnetic resonanceJ].Food Control,2023,145(10) :109448.

[11成實,王欣,劉寶林.碳鏈長度及不飽合度對脂肪酸低場核磁弛豫特性的影響[J].分析化學,2018,46(2):281-287.

[12]PAULAB,NITZANM,ALBERTOLC,etal.Studyofliquid-phase molecular packing interactions and morphology offattyacid methyl esters(biodiesel)[J].BiotechnologyforBiofuels,2015,8(1):12.

[13]USMANA,AMNA S,AYSHA S,et al. Use of Fouriertransform infrared spectroscopy and multi-variant analysisfordetection of butteradulteration withvegetable oil[J].InternationalJournal of Food Properties,2023,26(1):167-178.

[14]MAHBOUBIFARM,YOUSEFINEJADS,ALIZADEHM,etal.Predictionoftheacid value,peroxidevalueand thepercentage of some fatty acids in edible oils during longheating time by chemometrics analysis of FTIR-ATR spectra[J].Journal of theIranian Chemical Society,2016,13(12):2291-2299.

[15]馬鋼,馬浩翔,童俊杰,等.適合中老年人群食用的芝麻調和油的工藝研究[J].中國調味品,2023,48(7):91-96.

[16]史翠熠.基于LF-NMR技術快速鑒別植物油種類研究及應用[D].鄭州:河南工業大學,2024.

[17]ZHOUYY,CUIYW,WANGC,etal.Rapidandaccuratemonitoring and modeling analysis of eight kinds of nut oilsduring oil oxidationprocessbased on Fourier transform infraredspectroscopy[J].Food Control,2021,130(13):108294.

[18]徐飛,陳金艷,姚美芹,等.三種不同產地香莢蘭油理化特性及抗氧化特性研究[J].中國調味品,2022,47(1):90-97.

[19]黃連琴,王興進,盧宗桂,等.基于脂肪酸和生育酚組成的油茶籽油摻假判別可行性分析[J].中國油脂,2019,44(1):105-109.

[20]ZHANG Y,WANG M Z, ZHANG X P,et al. Mechanism,indexes,methods,challenges,and perspectives of edible oiloxidation analysis[J].Critical Reviews in Food Science andNutrition,2021,63(21) :11-15.

[21]杜蘅,胡毓元,蓋爭艷,等.低場核磁共振技術在油脂摻偽鑒別中的應用[J].中國糧油學報,2019,34(3):105-110,116.

[22]韓艷紅,楊海棠,胡延嶺,等.基于主成分和聚類分析的高油花生品種遺傳多樣性評估及綜合評價[J].四川農業大學學報,2024,42(1):28-38.

[23]AHMMED F, GORDON KC,KILLEEN D P,et al. Detectionand quantification of adulteration in krill oil with Raman andinfrared spectroscopicmethods[J].Molecules,2023,28(9):3695.

[24]何優選,梁奇峰.食用油品質的紅外吸收光譜分析[J].化工技術與開發,2011,40(6):32-34.

[25]朱銳,李祥慧,易陽.菜籽油品質多樣性與紅外光譜特征的研究[J].中國油脂,2022,47(11):126-135.

[26]韓建勛,孫瑞雪,陳穎,等.傅里葉變換紅外光譜結合化學計量學用于山茶油中摻雜大豆油的鑒別[J].食品與發酵工業,2019,45(18):222-227.

猜你喜歡
玉米油葵花籽橄欖油
地中海飲食法:一人食的四季餐桌(三)
食品與生活(2025年8期)2025-08-13 00:00:00
長壽花 十年翼路同行,共繪健康“食光”
新航空(2025年7期)2025-07-31 00:00:00
惟有葵花向日傾
主站蜘蛛池模板: 99久久精品国产麻豆婷婷| 久久福利片| 久久久久久久久18禁秘| 日韩黄色大片免费看| 99手机在线视频| 狠狠干综合| 国产三级精品三级在线观看| 亚洲欧美成人综合| 国产亚洲精品资源在线26u| 国产拍揄自揄精品视频网站| 国产黑丝一区| 国产99视频在线| 中国成人在线视频| 91精品国产自产91精品资源| 欧美成人一级| 亚洲中字无码AV电影在线观看| 最新国语自产精品视频在| 久久久久久尹人网香蕉| 国产一区二区丝袜高跟鞋| 欧美日韩一区二区在线播放| 67194成是人免费无码| 天天婬欲婬香婬色婬视频播放| 欧美精品啪啪一区二区三区| 久久国产香蕉| 91在线免费公开视频| 亚洲女人在线| 精品小视频在线观看| 色久综合在线| 97在线观看视频免费| 亚洲一级色| 国产成人久久777777| 国产免费怡红院视频| 亚洲成av人无码综合在线观看| 暴力调教一区二区三区| 国产精品第页| 成人中文字幕在线| 久久亚洲国产一区二区| 性欧美在线| 中文字幕无码制服中字| 四虎影视永久在线精品| 亚洲黄色激情网站| 国产嫩草在线观看| 国产成人精品男人的天堂下载 | 国产成人一级| 色噜噜中文网| 欧美中出一区二区| 92午夜福利影院一区二区三区| 亚洲精品国产乱码不卡| 久久国产精品夜色| 国内毛片视频| 欧美精品啪啪| 91丝袜美腿高跟国产极品老师| 亚洲伦理一区二区| 欧美性色综合网| 亚洲最新在线| 久久婷婷五月综合97色| 精品国产成人国产在线| 成人一级黄色毛片| 国产精品成人久久| 亚洲视频一区| 天堂网国产| 91精品日韩人妻无码久久| 国产午夜人做人免费视频| 一级黄色欧美| 无码高潮喷水专区久久| 在线欧美日韩| 国产美女丝袜高潮| 午夜久久影院| 国产免费久久精品99re不卡 | 亚洲91在线精品| 国产真实二区一区在线亚洲| 欧美午夜性视频| 国产全黄a一级毛片| 91精品伊人久久大香线蕉| 欧美在线导航| 欧美一级高清免费a| 国产麻豆永久视频| 亚洲美女一区| 91成人在线免费视频| 波多野结衣无码AV在线| 国产自产视频一区二区三区| 国内精品九九久久久精品|