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

基于2 種原理近紅外光譜儀對輻照花生的快速鑒別比較

2016-11-14 11:20:09靳欣欣馬智宏王紀(jì)華潘立剛
食品科學(xué) 2016年8期
關(guān)鍵詞:劑量檢測模型

王 冬,李 安,靳欣欣,馬智宏,王紀(jì)華,潘立剛

(北京農(nóng)業(yè)質(zhì)量標(biāo)準(zhǔn)與檢測技術(shù)研究中心,農(nóng)業(yè)部農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評估實(shí)驗(yàn)室,北京 100097)

基于2 種原理近紅外光譜儀對輻照花生的快速鑒別比較

王 冬,李 安,靳欣欣,馬智宏,王紀(jì)華,潘立剛*

(北京農(nóng)業(yè)質(zhì)量標(biāo)準(zhǔn)與檢測技術(shù)研究中心,農(nóng)業(yè)部農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評估實(shí)驗(yàn)室,北京 100097)

為比較不同近紅外光譜儀在對輻照花生鑒別方面的差異,采用傅里葉變換、線性漸變分光2 種原理的近紅外光譜儀針對不同輻照劑量的花生采集了近紅外光譜數(shù)據(jù),結(jié)合一階導(dǎo)數(shù)、標(biāo)準(zhǔn)正態(tài)變量變換預(yù)處理,對0、1、2 kGy輻照劑量的花生建立偏最小二乘-判別分析模型,采用外部驗(yàn)證集量本對模型正確率進(jìn)行驗(yàn)證。結(jié)果顯示,線性漸變分光-近紅外光譜儀的模型校正集正確率皆大于92%,交互驗(yàn)證集正確率皆大于85%,外部驗(yàn)證集正確率皆大于83%,普遍高于傅里葉變換-近紅外光譜儀模型的正確率。因此,可得出結(jié)論:預(yù)處理對模型正確率的影響相對較小,不同儀器對模型正確率影響相對較大。

線性漸變分光;傅里葉變換;近紅外光譜;輻照;花生

王冬, 李安, 靳欣欣, 等. 基于2 種原理近紅外光譜儀對輻照花生的快速鑒別比較[J]. 食品科學(xué), 2016, 37(8): 212-216.

WANG Dong, LI An, JIN Xinxin, et al. Comparison of rapid identification of irradiated peanuts based on two near infrared(NIR) spectrometers[J]. Food Science, 2016, 37(8): 212-216. (in Chinese with English abstract) DOI:10.7506/spkx1002-6630-201608038. http://www.spkx.net.cn

放射源所產(chǎn)生的電離輻射含有巨大的能量和強(qiáng)大的殺傷力,因此可用于一些食品的防腐等加工處理[1],并且具有延緩谷物老化、延長保鮮期的作用[2]。然而,花生等富含脂質(zhì)的食品一旦經(jīng)過電離輻射的照射,其中所含的脂肪酸,特別是含有雙鍵的不飽和脂肪酸極易發(fā)生斷裂,造成脂肪酸分解,其結(jié)果類似于強(qiáng)烈氧化,導(dǎo)致食品品質(zhì)下降,甚至發(fā)生酸敗,嚴(yán)重影響飲食者的身體健康。

對食品是否經(jīng)過電離輻射照射的檢測,傳統(tǒng)方法是采用電子自旋共振光譜法[3-4]、揮發(fā)性碳?xì)浠衔锓ǎ?]、熱釋光/光釋光分析法[6]、標(biāo)志物檢測法[7]、提取礦物質(zhì)法及自由基檢測法[8]等,不僅分析時(shí)間長,而且在其分析過程中使用大量化學(xué)試劑,對環(huán)境造成污染。

近年來,近紅外光譜(near-infrared spectroscopy,NIR)技術(shù)隨著計(jì)算機(jī)技術(shù)的進(jìn)步和化學(xué)計(jì)量學(xué)的發(fā)展在工農(nóng)業(yè)生產(chǎn)等諸多領(lǐng)域得到成功應(yīng)用[9-11]。由于具有操作簡單、幾乎不需要量品前處理以及可實(shí)現(xiàn)多指標(biāo)同時(shí)測定等特點(diǎn),NIR技術(shù)在產(chǎn)品質(zhì)量控制[12]、過程分析[13-14]、品質(zhì)分級[15]、真?zhèn)巫R別[16]等諸多方面具有成功的應(yīng)用;尤其在食品加工[17]、食品質(zhì)量快速測定[18-19]、農(nóng)作物病害快速檢測[20]等方面,NIR技術(shù)以其快速[21]、無損[22-23]、高效[24-25]的特點(diǎn)得到了廣泛的應(yīng)用。

至今,NIR儀主要?dú)v經(jīng)了濾光片型、棱鏡/光柵分光型、傅里葉變換(Fourier transformation,F(xiàn)T)型、線性漸變分光(linear variable filter,LVF)型等幾種不同類型。由于分光原理不同,各種儀器所獲取的數(shù)據(jù)會有不同程度的差異,進(jìn)而會影響NIR分析模型的效果。目前在工農(nóng)業(yè)生產(chǎn)中應(yīng)用最廣泛的屬FT型。該類型NIR儀通過邁克爾遜干涉儀中動鏡位置的改變產(chǎn)生光程差,從而產(chǎn)生干涉光,進(jìn)而實(shí)現(xiàn)分光的目的;通過FT將時(shí)域信號轉(zhuǎn)換為頻域信號,即得NIR圖。一般地,F(xiàn)T-NIR儀的光譜分辨率、光譜范圍用波數(shù)表示。FT-NIR儀由于使用邁克爾遜干涉儀,一般不可隨意移動,對環(huán)境要求較高。近年來,微制造技術(shù)的興起為NIR儀增加了新的類型。LVF-NIR儀采用精確加工的楔形透光材料作為分光器件,特定波長的光通過楔形透明材料時(shí),只有與指定位置相長干涉的光方可通過,并被微陣列型檢測器接收。一般情況下LVF-NIR儀的光譜分辨率、光譜范圍用波長表示。LVF-NIR儀無移動部件,可隨意移動;該類型光譜儀容易實(shí)現(xiàn)設(shè)備便攜化、微型化,可以實(shí)現(xiàn)對量品的現(xiàn)場、原位的無損快速分析。

本研究以不同輻照劑量的花生為研究對象,針對FT、LVF兩種原理的NIR儀所獲取的數(shù)據(jù),用2 種不同原理的NIR儀所獲取的數(shù)據(jù)分別建立不同輻照劑量花生的鑒別模型;比較了使用不同原理的NIR儀鑒別不同輻照劑量花生的正確率。本研究可為NIR儀的選型提供一定的參考,并可為輻照食品的快速鑒別提供一定的技術(shù)支持。

1 材料與方法

1.1材料與試劑

收集產(chǎn)自山東省的白花生量本共50 份,收集產(chǎn)自吉林省的紅花生量本50 份,共收集花生量本100 份。本研究中,花生的貯藏溫度(15±2) ℃,花生含水率6%~8%。

1.2儀器與設(shè)備

Spectrum 400 FT-NIR儀 美國Perkin Elmer公司;MicroNIR 1700 LVF-NIR儀 美國JDSU公司。

1.3方法

1.3.1分組處理方法

將每份花生量本均分為A、B、C三組,A組不經(jīng)過輻照處理,B組、C組分別送入60Co放射源輻照室進(jìn)行輻照處理,輻照劑量分別為1、2 kGy。

1.3.2儀器條件

FT-NIR儀:波數(shù)范圍12 800~4 000 cm-1,波數(shù)分辨率8 cm-1,累加平均次數(shù)64 次。

LVF-NIR儀:波長范圍908~1 676 nm,中心波長分辨率10 nm,累加平均次數(shù)50 次。

1.4數(shù)據(jù)處理

針對2 種NIR儀所獲取的數(shù)據(jù),分別采用一階導(dǎo)數(shù)、標(biāo)準(zhǔn)正態(tài)變量變換(standard normal variable,SNV)算法對數(shù)據(jù)進(jìn)行預(yù)處理,以消除數(shù)據(jù)平移、突出數(shù)據(jù)特征。針對2 種不同原理的NIR儀所獲取的數(shù)據(jù),以白花生、紅花生及上述2 種花生的輻照值為化學(xué)值,采用偏最小二乘-判別分析(partial least square-discrimination analysis,PLS-DA)算法結(jié)合全交互驗(yàn)證算法用2 種原理的NIR儀所獲取數(shù)據(jù)分別建立不同輻照劑量花生的鑒別模型,并用外部驗(yàn)證集量本對所建模型進(jìn)行盲量驗(yàn)證。

2 結(jié)果與分析

圖1 基于FT-NIR儀的白花生、紅花生的原光譜和一階導(dǎo)數(shù)譜Fig.1 FT-NIR original and first-derivative spectra of white-skinned peanut and red-skinned peanut

由圖1可見,紅花生原光譜在短波處離散程度較大,一階導(dǎo)數(shù)在一定程度上減弱了上述光譜數(shù)據(jù)的離散程度,SNV處理后,紅花生在全域離散程度大于白花生。

圖2 基于LVF-NIR儀的白花生、紅花生的原光譜和一階導(dǎo)數(shù)譜Fig.2 LVF-NIR original and first-derivative spectra of white-skinned peanut and red-skinned peanut

由圖2可見,白花生、紅花生原光譜離散程度差異不大,一階導(dǎo)數(shù)、SNV處理后,白花生、紅花生的光譜數(shù)據(jù)離散差異未見明顯差異。

通過上述光譜數(shù)據(jù)可見,采用FT-NIR儀所采集數(shù)據(jù),白花生的數(shù)據(jù)變異程度較小;采用LVF-NIR儀所采集數(shù)據(jù),白花生、紅花生光譜數(shù)據(jù)的變異程度相當(dāng)。FTNIR儀譜區(qū)范圍較大,而LVF-NIR儀譜區(qū)范圍較小。2 種儀器所采集的NIR數(shù)據(jù)在譜形上具有一定的差異,僅從NIR的外觀上無法對各種輻照劑量的花生進(jìn)行區(qū)分,需要進(jìn)一步對數(shù)據(jù)信息進(jìn)行特征提取,建立校正模型。

分別采用一階導(dǎo)數(shù)、SNV兩種方法做數(shù)據(jù)預(yù)處理,以消除光譜平移、突出數(shù)據(jù)特征,針對白花生、紅花生和上述2 種花生,采用PLS-DA算法結(jié)合全交互驗(yàn)證算法建立0、1、2 kGy輻照劑量花生的鑒別模型;采用外部驗(yàn)證集對模型進(jìn)行盲量驗(yàn)證;分別針對校正集、交互驗(yàn)證集、外部驗(yàn)證集統(tǒng)計(jì)3 種輻照劑量及總體的正確率,結(jié)果如表1所示。

表1 2 種NIR儀模型正確率匯總表Table 1 Summary of the accuracies of the two NIR models

表1顯示,對于FT-NIR儀,白花生、紅花生以及上述2 種花生3 個(gè)模型的校正集正確率皆不低于81%,而交互驗(yàn)證集、外部驗(yàn)證集正確率皆較低,除個(gè)別情況,普遍不超過80%;對于LVF-NIR儀,3 個(gè)模型的校正集正確率皆高于92%,而交互驗(yàn)證集正確率皆高于85%,外部驗(yàn)證集正確率皆高于83%。其中,LVF-NIR儀數(shù)據(jù)結(jié)合一階導(dǎo)數(shù)預(yù)處理的模型,其校正集、交互驗(yàn)證集、外部驗(yàn)證集的正確率皆不低于90%。上述結(jié)果說明,針對3 種輻照劑量花生的分辨,LVF-NIR儀可以得到正確率更高的鑒別結(jié)果。比較不同的數(shù)據(jù)預(yù)處理方法可見,F(xiàn)T-NIR儀所測數(shù)據(jù)經(jīng)過一階導(dǎo)數(shù)、SNV預(yù)處理后建模,模型正確率基本持平;LVF-NIR儀所測數(shù)據(jù)經(jīng)一階導(dǎo)數(shù)預(yù)處理后建模,模型正確率普遍優(yōu)于SNV預(yù)處理的結(jié)果,但數(shù)據(jù)預(yù)處理的影響較儀器原理的影響較小。

造成這一結(jié)果的原因是,基于微制造技術(shù)的LVFNIR儀整體無移動部件,所獲取的光譜數(shù)據(jù)更穩(wěn)定;而FT-NIR儀中采用了邁克爾遜干涉儀作為分光器件,其構(gòu)造以及動鏡的移動增加了所獲取的光譜數(shù)據(jù)的不確定性。此外,F(xiàn)T-NIR儀直接獲取時(shí)域譜,經(jīng)過FT后轉(zhuǎn)換為頻域譜;而LVF-NIR儀沒有上述數(shù)據(jù)轉(zhuǎn)換步驟,從而進(jìn)一步減小了數(shù)據(jù)的不確定性。

3 結(jié) 論

本研究采用FT、LVF 2 種原理的NIR儀針對不同輻照劑量的花生采集了NIR,針對白花生、紅花生及上述2 種花生建立了0、1、2 kGy三種輻照劑量的PLS-DA模型,采用外部驗(yàn)證集量本對模型正確率進(jìn)行了盲量驗(yàn)證。結(jié)果表明,采用NIR技術(shù)對輻照花生進(jìn)行快速鑒別具有一定的可行性;不同儀器對模型正確率影響相對較大,而數(shù)據(jù)預(yù)處理對模型正確率的影響相對較小;針對本研究中的花生量本,LVF-NIR儀的模型正確率普遍較高。

從本研究的結(jié)果看,針對本研究中所涉及的0、1、2 kGy輻照花生的鑒別,基于LVF-NIR儀數(shù)據(jù)的模型具有更好的鑒別結(jié)果。但是目前的制作工藝尚無法實(shí)現(xiàn)在一臺LVF-NIR儀上更改光譜分辨率和波長范圍。關(guān)于光譜分辨率和波長范圍對LVF-NIR儀的模型性能的影響仍需要進(jìn)一步深入研究。

[1] 徐麗婧, 高麗樸, 王清, 等. 輻照保鮮技術(shù)及其在雙孢蘑菇保鮮中的應(yīng)用[J]. 食品工業(yè)科技, 2014, 35(9): 392-395. DOI:10.13386/ j.issn1002-0306.2014.09.078.

[2] 陳銀基, 陳霞, 蔣偉鑫, 等. γ-輻照處理對谷物儲藏及品質(zhì)特性的影響研究進(jìn)展[J]. 食品工業(yè)科技, 2014, 35(4): 358-362; 368. DOI:10.13386/j.issn1002-0306.2014.04.038.

[3] 王鋒, 王志東, 周洪杰, 等. 輻照食品檢測技術(shù)方法與標(biāo)準(zhǔn)[J]. 食品與機(jī)械, 2009, 25(5): 142-144. DOI:10.13652/ j.issn.1003-5788.2009.05.018.

[4] 周鑫瑤, 徐芙蓉, 趙小俊, 等. ESR波譜技術(shù)檢測輻照農(nóng)產(chǎn)品[J]. 中國食品學(xué)報(bào), 2013, 13(3): 202-207. DOI:1009-7848(2013)03-0202-06.

[5] 朱珍, 王寧, 李振興. 輻照食品檢測技術(shù)的研究進(jìn)展[J]. 中國漁業(yè)質(zhì)量與標(biāo)準(zhǔn), 2015, 5(2): 56-64. DOI:2095-1833(2015)02-0056-09.

[6] 趙良娟, 張海濱, 曲鵬, 等. 輻照食品檢測標(biāo)準(zhǔn)及檢測方法研究進(jìn)展[J].食品研究與開發(fā), 2012, 33(9): 208-213.

[7] 胡娟, 汪美鳳, 趙國華. 輻照食品標(biāo)志物及其檢測技術(shù)研究進(jìn)展[J]. 食品與發(fā)酵工業(yè), 2010, 36(7): 127-130. DOI:10.13995/ j.cnki.11-1802/ts.2010.07.022.

[8] 哈益明, 周洪杰. 輻照食品分析檢測技術(shù)的研究進(jìn)展[J]. 食品科學(xué),2005, 26(6): 260-265.

[9] ALEXANDRINO G L, KHORASANI M R, AMIGO J M, et al. Monitoring of multiple solid-state transformations at tablet surfaces using multi-series near-infrared hyperspectral imaging and multivariate curve resolution[J]. European Journal of Pharmaceutics and Biopharmaceutics, 2015, 93: 224-230. DOI:10.1016/ j.ejpb.2015.03.034.

[10] SULUB Y, LOBRUTTO R, VIVILECCHIA R, et al. Content unifirmity determination of pharmaceutical tablets using five nearinfrared reflectance spectrometers: a process analytical technology(PAT) approach using robust multivariate calibration transfer algorithms[J]. Analytica Chimica Acta, 2008, 611(2): 143-150. DOI:10.1016/j.aca.2008.02.016.

[11] VIEIRA F S, PASQUINI C. Near infrared emission photometer for measuring the oxidative stability of edible oils[J]. Analytica Chimica Acta, 2013, 796: 101-107. DOI:10.1016/j.aca.2013.08.028.

[12] FERREIRA D S, GAL?O O F, PALLONE J A L, et al. Comparison and application of near-infrared (NIR) and mid-infrared (MIR)spectroscopy for determination of quality parameters in soybean samples[J]. Food Control, 2014, 35(1): 227-232. DOI:10.1016/ j.foodcont.2013.07.010.

[13] SCHAEFER C, LECONTE C, CLICQ D, et al. On-line near infrared spectroscopy as a process analytical technology (PAT) tool to control an industrial seeded API crystallization[J]. Journal of Pharmaceutical and Biomedical Analysis, 2013, 83: 194-201. DOI:10.1016/ j.jpba.2013.05.015.

[14] 劉逸, 張兆斌, 張永剛, 等. 近紅外技術(shù)快速測定裂解焦油芳烴含量[J]. 光學(xué)精密工程, 2014, 22(7): 1723-1731. DOI:10.3788/ OPE.20142207.1723.

[15] VASQUES G M, DEMATTê J A M, ROSSEL R A V, et al. Soil classification using visible/near-infrared diffuse reflectance spectra from multiple depths[J]. Geoderma, 2014, 223/225: 73-78. DOI:10.1016/j.geoderma.2014.01.019.

[16] 王冬, 潘立剛, 王紀(jì)華, 等. 基于線性漸變分光微型近紅外儀的西湖龍井真?zhèn)文P筒贿m應(yīng)性析因及修正研究[J]. 光譜學(xué)與光譜分析, 2014, 34(11): 2938-2943. DOI:10.3964/j.is sn.1000-0593(2014)11-2938-06.

[17] 王麗, 韓曉祥, 徐歡, 等. 近紅外光譜法測定金華火腿生產(chǎn)過程中的品質(zhì)參數(shù)[J]. 中國食品學(xué)報(bào), 2011, 11(7): 181-184. DOI:1009-7848(2011)07-0181-04.

[18] 陳天華, 雷春寧, 李月. 基于近紅外光譜特性分析的花生含水率檢測[J]. 食品科學(xué)技術(shù)學(xué)報(bào), 2013, 31(5): 50-54.

[19] 董小玲, 孫旭東. 基于小波壓縮的馬鈴薯全粉還原糖近紅外光譜檢測研究[J]. 光譜學(xué)與光譜分析, 2013, 33(12): 3216-3220. DOI:10.3964/j.issn.1000-0593(2013)12-3216-05.

[20] 周竹, 李小昱, 高海龍, 等. 漫反射和透射光譜檢測馬鈴薯黑心病的比較[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2012, 28(11): 237-242. DOI:10.3969/ j.issn.1002-6819.2012.11.038.

[21] 賴立群, 麻望瓊, 陳華才. 近紅外漫反射光譜結(jié)合神經(jīng)網(wǎng)絡(luò)檢測小麥蛋白質(zhì)含量[J]. 中國計(jì)量學(xué)院學(xué)報(bào), 2015, 26(1): 55-59. DOI:10.3969/j.issn.1004-1540.2015.01.010.

[22] 夏立婭, 申世剛, 劉崢顥, 等. 基于近紅外光譜和模式識別技術(shù)鑒別大米產(chǎn)地的研究[J]. 光譜學(xué)與光譜分析, 2013, 33(1): 102-105. DOI:10.3964/j.issn.1000-0593(2013)01-0102-04.

[23] 王明, 于峰, 劉新, 等. 采用近紅外漫反射技術(shù)對牛奶中蛋白質(zhì)、脂肪檢測[J]. 激光雜志, 2015, 36(1): 70-73. DOI:10.14016/j.cnki. jgzz.2015.01.070.

[24] 黃偉, 楊秀娟, 曹志勇, 等. 近紅外反射光譜快速檢測滇南小耳豬肉中水分、粗脂肪及粗蛋白含量的研究[J]. 中國畜牧雜志, 2015,51(7): 73-77. DOI:0258-7033(2015)07-0073-05.

[25] 李琳琳, 金華麗, 崔彬彬, 等. 基于近紅外透射光譜的大豆蛋白質(zhì)和粗脂肪含量快速檢測[J]. 糧食與油脂, 2014, 27(12): 57-60. DOI:1008-9578(2014)1-0057-04.

Comparison of Rapid Identification of Irradiated Peanuts Based on Two Near Infrared (NIR) Spectrometers

WANG Dong, LI An, JIN Xinxin, MA Zhihong, WANG Jihua, PAN Ligang*
(Risk Assessment Laboratory for Agro-products (Beijing), Ministry of Agriculture,Beijing Research Center for Agricultural Standards and Testing, Beijing 100097, China)

In order to compare the difference in the identification of60Co γ-ray irradiated peanuts by different near-infrared(NIR) spectrometers, NIR spectra of peanuts exposed to different dosages of irradiation were collected by Fourier transform(FT) spectrometer and linear variable filter (LVF) spectrometer, respectively. The spectra were pretreated by the first derivative and the standard normal variable (SNV), respectively. The partial least square-discriminat analysis (PLS-DA) calibration models were developed and validated using external validation sets. The results indicated that the accuracy of the LVF-NIR model for calibration, cross validation and external validation sets were larger than 92%, larger than 85% and larger than 83% respectively,which were all larger than those of the FT-NIR model. It could be concluded that different data preprocessing methods have less influence on the model accuracy, while the influence of different NIR spectrometers is significant.

linear variable filter; Fourier transform; near-infrared spectroscopy; irradiation; peanut

10.7506/spkx1002-6630-201608038

O657.33

A

1002-6630(2016)08-0212-05

10.7506/spkx1002-6630-201608038. http://www.spkx.net.cn

2015-05-22

北京市農(nóng)林科學(xué)院青年基金項(xiàng)目(QNJJ201431);“十二五”國家科技支撐計(jì)劃項(xiàng)目(2014BAD04B05)

王冬(1982—),男,副研究員,博士,研究方向?yàn)檎駝庸庾V分析及化學(xué)計(jì)量學(xué)。E-mail:nirphd@163.com

潘立剛(1964—),男,研究員,博士,研究方向?yàn)檗r(nóng)藥毒理學(xué)。E-mail:panlig@126.com

引文格式:

猜你喜歡
劑量檢測模型
一半模型
結(jié)合劑量,談輻射
·更正·
“不等式”檢測題
“一元一次不等式”檢測題
“一元一次不等式組”檢測題
90Sr-90Y敷貼治療的EBT3膠片劑量驗(yàn)證方法
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
3D打印中的模型分割與打包
主站蜘蛛池模板: 亚欧美国产综合| 丰满人妻被猛烈进入无码| 国产成人综合在线观看| 波多野结衣在线一区二区| 亚洲熟女偷拍| 色欲色欲久久综合网| 欧美精品H在线播放| 久久中文无码精品| 日韩经典精品无码一区二区| 在线观看国产一区二区三区99| 综合亚洲色图| 一级毛片无毒不卡直接观看 | 亚洲va欧美va国产综合下载| 国产乱子伦一区二区=| 91小视频在线| 色哟哟色院91精品网站| 国产门事件在线| 麻豆精品在线视频| 中文字幕 91| 五月天久久综合| 亚洲一区二区三区国产精华液| 欧美日韩一区二区三区在线视频| 亚洲久悠悠色悠在线播放| 91伊人国产| 亚洲精品777| 米奇精品一区二区三区| 丝袜亚洲综合| 亚洲高清无码精品| 精品无码一区二区在线观看| 国产精品浪潮Av| 国产在线一二三区| 婷婷色在线视频| 伊人久热这里只有精品视频99| 国产一级妓女av网站| 免费看一级毛片波多结衣| 亚洲精品国产自在现线最新| 69精品在线观看| 精品91自产拍在线| 亚洲精品成人片在线观看| 国产原创第一页在线观看| 99re经典视频在线| 欧美无遮挡国产欧美另类| 日韩高清在线观看不卡一区二区 | 亚洲国产精品人久久电影| 这里只有精品国产| 一本大道香蕉久中文在线播放| 中文无码精品a∨在线观看| 色综合天天视频在线观看| 国产日韩欧美视频| 亚洲精品第一在线观看视频| 91黄色在线观看| 国产成人亚洲无码淙合青草| 国产亚洲视频免费播放| 精品国产www| 精品一区二区三区水蜜桃| 色网站在线免费观看| 国产精品免费露脸视频| 亚洲精选无码久久久| 伊人成人在线视频| 免费观看精品视频999| 97免费在线观看视频| 波多野结衣一二三| 人妻精品全国免费视频| 久久天天躁狠狠躁夜夜2020一| 国产精品手机在线播放| 日韩高清一区 | 久久综合丝袜长腿丝袜| 日本成人精品视频| 天天色天天操综合网| 中文字幕日韩视频欧美一区| 欧美一道本| 激情无码字幕综合| 日韩a在线观看免费观看| 亚洲区一区| 亚洲天堂网在线视频| 国产美女91呻吟求| 天堂成人av| 亚洲黄色视频在线观看一区| 日韩成人在线网站| 五月综合色婷婷| 欧美在线三级| 欧美国产三级|