• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    糧油食品中微生物和真菌毒素污染預(yù)測(cè)模型研究進(jìn)展

    2020-05-24 07:05:36裴鵬鋼王小藝趙峙堯柳佳蕓
    中國(guó)糧油學(xué)報(bào) 2020年2期
    關(guān)鍵詞:糧油食品黃曲霉霉菌

    裴鵬鋼 熊 科 葉 宏 王小藝 趙峙堯 柳佳蕓 鄧 蕾

    (北京工商大學(xué)北京市食品添加劑工程技術(shù)研究中心1,北京 100048)(北京工商大學(xué)北京市風(fēng)味化學(xué)重點(diǎn)實(shí)驗(yàn)室2,北京 100048)(北京工商大學(xué)北京市質(zhì)量與安全重點(diǎn)實(shí)驗(yàn)室3,北京 100048)(北京工商大學(xué)北京食品營(yíng)養(yǎng)與人類(lèi)健康高精尖創(chuàng)新中心4,北京 100048)

    1 糧油食品預(yù)測(cè)微生物學(xué)

    稻谷、小麥、大豆等大宗糧油食品是我國(guó)糧油食品加工的主要種類(lèi)。而在這些糧油食品生產(chǎn)加工過(guò)程中,供應(yīng)鏈環(huán)節(jié)鏈條長(zhǎng),各環(huán)節(jié)污染物種類(lèi)多、來(lái)源復(fù)雜,致使糧油食品從原料種植到儲(chǔ)藏運(yùn)輸?shù)炔煌A段均易遭受微生物尤以霉菌及其產(chǎn)生毒素污染。迄今為止在糧油食品中已發(fā)現(xiàn)300多種真菌毒素,主要為黃曲霉毒素、脫氧雪腐鐮刀菌烯醇、赭曲霉毒素、展青霉素、單端孢霉烯族毒素、玉米赤霉烯酮、伏馬毒素等[1]。不同類(lèi)型糧油食品中的產(chǎn)毒霉菌并非完全相同,即使同一種糧油食品中加工或儲(chǔ)藏不同階段或不同的環(huán)境條件下的有害菌也不盡相同[2]。

    以小麥供應(yīng)鏈主要環(huán)節(jié)為例,在其種植期間容易遭受赤霉病、葉銹病、紋枯病等真菌病害的侵染,造成小麥減產(chǎn),并且產(chǎn)生的真菌毒素在小麥中形成積累;其次,在小麥原糧倉(cāng)儲(chǔ)階段由于糧堆各部分貯藏溫度不同,再加上外界條件變化導(dǎo)致濕度變化使得小麥也極易遭受鐮刀菌、青霉菌、交聯(lián)孢菌及其次級(jí)代謝產(chǎn)物如脫氧雪腐鐮刀菌烯醇(DON)、玉米赤霉烯酮(ZEN)、黃曲霉毒素(AFB1)的污染[3]。同時(shí)小麥制粉階段潤(rùn)麥環(huán)節(jié)兩次潤(rùn)麥處理,導(dǎo)致糧堆水分含量增加,潤(rùn)麥倉(cāng)倉(cāng)儲(chǔ)環(huán)節(jié)為微生物的生長(zhǎng)繁殖創(chuàng)造良好的條件。同時(shí)小麥粉在貯藏階段衛(wèi)生狀況也會(huì)影響小麥粉中的微生物數(shù)量;在產(chǎn)品運(yùn)輸環(huán)節(jié),環(huán)境的差異同樣也會(huì)導(dǎo)致小麥粉霉菌等腐敗菌的數(shù)量增加。因此小麥供應(yīng)鏈中存在真菌毒素污染物影響產(chǎn)品安全的較大風(fēng)險(xiǎn),而對(duì)于供應(yīng)鏈的污染物難以逐一進(jìn)行分析研究,只能聚焦于關(guān)鍵階段對(duì)其進(jìn)行控制分析。

    預(yù)測(cè)食品微生物學(xué)(Predictive Food Microbiology)是結(jié)合微生物學(xué)、數(shù)學(xué)、統(tǒng)計(jì)學(xué)和計(jì)算機(jī)等學(xué)科,采用數(shù)學(xué)的方法描述微生物數(shù)量變化和外部環(huán)境因素之間的關(guān)系,并對(duì)微生物的生長(zhǎng)動(dòng)態(tài)情況做出預(yù)測(cè)。因此為了確保食品安全性,科學(xué)家們通過(guò)引入預(yù)測(cè)食品微生物學(xué)理論對(duì)食品中的微生物進(jìn)行量化和行為預(yù)測(cè)。通過(guò)計(jì)算機(jī)軟件進(jìn)行病原菌和腐敗菌生長(zhǎng)及其產(chǎn)生的污染物在特定條件下的概率模型等的擬合[4],快速判斷食品內(nèi)主要病原菌和腐敗菌生長(zhǎng)或殘存的動(dòng)態(tài)變化以及有害物含量的高低,從而依據(jù)標(biāo)準(zhǔn)對(duì)特定階段食品的安全性作出快速評(píng)價(jià)。表1為近年來(lái)國(guó)內(nèi)外對(duì)食品如玉米、水稻、麥芽等食品中的霉菌生長(zhǎng)動(dòng)力學(xué)模擬研究情況。

    因此,針對(duì)糧油食品供應(yīng)鏈污染物特點(diǎn),采用預(yù)測(cè)微生物學(xué)模型的方式對(duì)其供應(yīng)鏈中關(guān)鍵環(huán)節(jié)的毒素污染產(chǎn)生菌構(gòu)建動(dòng)力學(xué)模型以及通過(guò)不同條件下產(chǎn)毒概率模型,探究關(guān)鍵環(huán)節(jié)毒素變化及降解遷移規(guī)律,從而實(shí)現(xiàn)對(duì)糧油產(chǎn)品關(guān)鍵環(huán)節(jié)品質(zhì)的安全評(píng)估與預(yù)測(cè)。通過(guò)模型的反饋來(lái)優(yōu)化糧油產(chǎn)品貯藏及加工的方法和手段,實(shí)現(xiàn)“綠色儲(chǔ)糧”和“生態(tài)儲(chǔ)糧”,為確保糧油食品安全提供重要的理論依據(jù)。這也為未來(lái)開(kāi)展糧油食品供應(yīng)鏈真菌毒素污染的安全風(fēng)險(xiǎn)評(píng)估以及對(duì)其供應(yīng)鏈中真菌毒素危害物預(yù)警追溯等奠定模型基礎(chǔ),滿(mǎn)足糧油食品供應(yīng)鏈安全與社會(huì)可持續(xù)性發(fā)展的需求。

    表1 近年國(guó)內(nèi)外關(guān)于食品霉菌生長(zhǎng)動(dòng)力學(xué)模擬情況

    2 糧油食品中真菌毒素污染情況

    國(guó)內(nèi)目前大宗糧油食品真菌毒素污染超標(biāo)問(wèn)題嚴(yán)重,對(duì)消費(fèi)者的身體健康帶來(lái)極大的安全隱患。稻谷、小麥、大豆、花生等大宗糧油食品中均發(fā)現(xiàn)真菌毒素不同程度污染的情況。

    表2 近年國(guó)內(nèi)糧油食品真菌毒素污染調(diào)查

    注:ND指的是未檢出相應(yīng)毒素。

    3 微生物預(yù)測(cè)模型

    微生物預(yù)測(cè)模型最早應(yīng)用于食品貯藏階段細(xì)菌生長(zhǎng)模型,同時(shí)期微生物學(xué)預(yù)測(cè)模型對(duì)霉菌并未有同等程度的關(guān)注研究,這與霉菌生長(zhǎng)的量化復(fù)雜性有關(guān)。細(xì)菌裂變繁殖通常發(fā)生在食物表面或者液體基質(zhì)而霉菌菌絲則可穿透食品的物理三維空間,所以對(duì)其定量分析更加復(fù)雜。并且產(chǎn)生有用且可靠的模型的先決條件必須是包含大量相關(guān)數(shù)據(jù)的數(shù)據(jù)庫(kù),這對(duì)于研究真菌毒素來(lái)說(shuō)無(wú)疑是巨大的挑戰(zhàn)。一般來(lái)說(shuō),微生物模型分為三類(lèi)。初級(jí)模型是用來(lái)描述微生物在一定生長(zhǎng)條件下其數(shù)量隨時(shí)間的變化關(guān)系,如孢子的萌發(fā)、菌株生長(zhǎng)、菌體的失活等[22]。Zwietering等[23]比較了幾種S形函數(shù)(logistic,Gompertz,Richards,Schnute和Stannard)來(lái)描述細(xì)菌生長(zhǎng)曲線(xiàn),通過(guò)t檢驗(yàn)、F檢驗(yàn)發(fā)現(xiàn)Gompertz曲線(xiàn)10 ℃條件下植物乳桿菌數(shù)量隨時(shí)間變化的擬合度最高(見(jiàn)圖1);而二級(jí)模型在一級(jí)模型的基礎(chǔ)上增加了環(huán)境因素,如:溫度、水分含量、pH等因素,表征環(huán)境條件的改變?nèi)绾斡绊懢甑拿劝l(fā)和生長(zhǎng)。圖2為評(píng)估水分活度為變量菌株徑向生長(zhǎng)速率二級(jí)模型擬合圖;三級(jí)模型是一種軟件化模型,軟件包括了1個(gè)或者多個(gè)一級(jí)模型和二級(jí)模型,作為完整的預(yù)測(cè)工具為用戶(hù)提供服務(wù)[24]。

    圖1 10 ℃生長(zhǎng)的植物乳桿菌相對(duì)種群大小隨時(shí)間變化Gompertz曲線(xiàn)

    圖2 22 ℃條件下模擬UG AF06菌株在全黑胡椒粒不同水分活度下菌株徑向最大生長(zhǎng)率實(shí)測(cè)模型擬合圖

    概率模型用于描述某一特定時(shí)間段某些事件發(fā)生的可能性,例如,特定條件下菌體是否生長(zhǎng)、是否產(chǎn)生毒素。通過(guò)概率模型可以對(duì)某種食品基質(zhì)在特定時(shí)間點(diǎn)的安全性做出快速評(píng)價(jià)[25]。

    3.1 微生物生長(zhǎng)初級(jí)模型

    常用的一級(jí)模型主要有Logistic方程[25]、Gompertz模型[25]、Baranyi & Roberts模型[26]、Boltzmann模型[27]、致死模型[28]。

    3.1.1 Logistic方程:

    3.1.2 Gompertz模型:

    y=a·exp[-exp(b-cx)]

    式中兩個(gè)模型中的參數(shù)只是數(shù)學(xué)意義上的參數(shù)(a,b,c,…),并不具有生物學(xué)意義,不能直接用于微生物生長(zhǎng)曲線(xiàn)的擬合。需要對(duì)其中的參數(shù)進(jìn)行修改,使用具有生物學(xué)意義的參數(shù)如(λ,μm,A)來(lái)代替純數(shù)學(xué)參數(shù)(a,b,c,…)。

    修正后的Logistic模型:

    修正后的Gompertz模型:

    式中:μm為最大比生長(zhǎng)速率,定義為生長(zhǎng)曲線(xiàn)拐點(diǎn)處的斜率;λ為遲滯期時(shí)間,滯后期被定義為定義最大特定生長(zhǎng)速率的線(xiàn)與時(shí)間軸的交點(diǎn);t為培養(yǎng)時(shí)間;A指微生物數(shù)量達(dá)最大值;N0為初始孢子數(shù);N為t時(shí)間后存活的孢子數(shù)量。Dantigny等[29]通過(guò)實(shí)驗(yàn)比較Logistic模型和Gompertz模型,通過(guò)比較RMSE值(均方根誤差)發(fā)現(xiàn)兩者擬合優(yōu)度類(lèi)似,無(wú)法確定哪種模型在特定實(shí)驗(yàn)條件下為首選,繼而通過(guò)評(píng)估其他標(biāo)準(zhǔn)如參數(shù)估計(jì)值的準(zhǔn)確性來(lái)進(jìn)一步判斷如:評(píng)估孢子萌發(fā)時(shí)間(50%孢子萌發(fā)所需時(shí)間)時(shí),Logistic模型的精確度更高;而Gompertz模型可以更準(zhǔn)確地估計(jì)遲滯期時(shí)間。

    3.1.3 Boltzmann模型

    3.1.4 Baranyi & Roberts模型

    式中:μm為最大比生長(zhǎng)速率(定義為拐點(diǎn)處生長(zhǎng)曲線(xiàn)的斜率)(mm/d,μm/h);λ為遲滯期時(shí)間;t為培養(yǎng)時(shí)間;A為微生物數(shù)量最大值。

    Lo′pez等[30]評(píng)估微生物生長(zhǎng)模型(光密度方法定量)確定包括真菌和細(xì)菌在內(nèi)物種,采用包括Gompertz和Baranyi & Roberts模型在內(nèi)的各種模型進(jìn)行擬合,發(fā)現(xiàn)Baranyi & Roberts模型擬合度最高;Marin等[31]研究lineal模型、Gompertz模型及Baranyi & Roberts模型對(duì)14種常見(jiàn)食物腐敗真菌生長(zhǎng)一級(jí)模型的適用性,采用菌落直徑進(jìn)行生長(zhǎng)速率表征發(fā)現(xiàn):當(dāng)所獲數(shù)據(jù)沒(méi)有發(fā)現(xiàn)漸進(jìn)趨勢(shì)時(shí),使用liner模型和Baranyi & Roberts模型可以更好地估計(jì)最大生長(zhǎng)速率和遲滯期;當(dāng)觀察到生長(zhǎng)速率隨時(shí)間降低時(shí),Baranyi & Roberts模型的擬合度更高。

    3.1.5 致死模型

    模擬溫度高低或致死劑濃度對(duì)微生物群體致死處理模型時(shí),使用經(jīng)典的一級(jí)方程對(duì)孢子滅活建立模型:

    式中:N為孢子數(shù);N0為初始孢子數(shù);N(t)為經(jīng)過(guò)處理后孢子數(shù);ki為致死率/min-1。

    現(xiàn)今大多數(shù)孢子滅活模型研究采用線(xiàn)性滅活方法,但是目前對(duì)于孢子滅活的研究還較少且只有少部分應(yīng)用于真菌[32]。

    3.2 微生物生長(zhǎng)二級(jí)模型

    二級(jí)模型主要有Ratkowsky平方根模型[32]、Arrhenius-Davey模型[33]、Rosso cardinal模型[34]、Polynomial模型[35]、Gamma Concept模型[36]。

    3.2.1 Ratkowsky平方根模型

    該模型最初是為細(xì)菌開(kāi)發(fā)的。

    式中:μm為最大比生長(zhǎng)速率;b2、C2均為假定系數(shù);Tmin2為最低生長(zhǎng)溫度;Tmax2為最高生長(zhǎng)溫度。通常該模型僅有溫度變量,Tassou等[37]在建立溫度和水分活度對(duì)合成葡萄汁培養(yǎng)基上希臘釀酒葡萄中2種產(chǎn)生赭曲霉A的炭曲霉菌株模型時(shí),將Ratkowsky平方根模型變量拓展到其他變量因素上如水分活度。

    3.2.2 修正的線(xiàn)性Arrhenius-Davey方程

    最初用于描述溫度對(duì)細(xì)菌生長(zhǎng)的影響,Panagou等[38]模擬溫度、pH和水分活度對(duì)從綠茶橄欖提取的紅曲霉生長(zhǎng)速度影響時(shí),將該模型變量因素延伸到pH和水分活度。

    lnμmax=C0+C1/T+C2/T2

    式中:μ為菌落生長(zhǎng)速率;T為開(kāi)氏溫度;C0、C1、C2為3個(gè)參數(shù)。

    3.2.3 Rosso cardinal模型

    Rosso等[34]在1993年提出了這個(gè)溫度模型,使用3個(gè)基本溫度(Tmin、Tmax、Topt),和最佳溫度下的生長(zhǎng)速率μopt。該模型是在經(jīng)驗(yàn)基礎(chǔ)上建立的作為描述數(shù)據(jù)的工具,模型的優(yōu)點(diǎn)是所有參數(shù)都具有生理學(xué)意義,有利于對(duì)初始參數(shù)的估計(jì)。

    μR(T,aw)=μopt·τ(T)·ρ(aw)

    式中:μR為菌落徑向生長(zhǎng)速率;μopt菌落直徑最大值;Tmin、Tmax、awmin為菌株可以生長(zhǎng)的最低和最高溫度以及最小水分活度;Topt、awopt為菌落直徑達(dá)到最大值時(shí)對(duì)應(yīng)的溫度和水分活度;awmax為1,方便計(jì)算無(wú)特殊意義。

    3.2.4 Gibson模型

    Gibson模型是第一個(gè)專(zhuān)門(mén)為霉菌開(kāi)發(fā)的模型,Gibson等[35]發(fā)現(xiàn)真菌生長(zhǎng)速率的對(duì)數(shù)與(1-aw)的平方根呈拋物線(xiàn)關(guān)系,所以將多項(xiàng)式模型應(yīng)用于黃曲霉的生長(zhǎng)。

    式中:μ為菌落生長(zhǎng)速率;a0、a1、a2為3個(gè)需要估計(jì)的參數(shù)。

    Tassou等[39]研究環(huán)境溫度(25 ℃)和選定溫度(45、55 ℃)下的高靜水壓力對(duì)火腿中金黃色葡萄球菌的滅活動(dòng)力學(xué)建模,在Gibson模型中增加溫度變量。

    式中:a0、a1、a2、a3、a4、a5為3個(gè)需要估計(jì)的參數(shù)。

    3.2.5 Gamma Concept模型

    Gamma Concept模型由Zwietering等[36]提出,認(rèn)為影響微生物生長(zhǎng)速率的各種因素可以通過(guò)單獨(dú)效應(yīng)的乘積來(lái)組合。使用這種方法的優(yōu)點(diǎn)是,對(duì)于確定影響生長(zhǎng)速率每個(gè)變量,可以通過(guò)分離各種變量的影響來(lái)計(jì)算相對(duì)效應(yīng)。這個(gè)模型有一般形式:

    其中給定變量的相對(duì)影響可以通過(guò)該變量的gamma因子來(lái)描述:

    3.3 概率模型

    當(dāng)實(shí)驗(yàn)?zāi)繕?biāo)是觀察在特定條件下微生物是否生長(zhǎng)時(shí),即可采用概率模型來(lái)進(jìn)行擬合。概率模型在很多情況下可確保食物不受微生物污染至關(guān)重要,尤其適用于涉及病原微生物或者產(chǎn)毒微生物污染的情況。概率模型允許在各種條件下預(yù)測(cè)特定事件(例如微生物的生長(zhǎng)和毒素的產(chǎn)生)是否可能發(fā)生,其中增長(zhǎng)概率模型是根據(jù)特定時(shí)間點(diǎn)在整個(gè)實(shí)驗(yàn)設(shè)計(jì)中“增長(zhǎng)/無(wú)增長(zhǎng)”響應(yīng)的比例構(gòu)建的[40]。

    Logistic回歸模型是判斷菌株增長(zhǎng)/無(wú)增長(zhǎng)邊界的有用工具,將事件發(fā)生概率Y與條件變量進(jìn)行關(guān)聯(lián)。

    式中:xi為因子變量(T,aw,pH,…);bi為要估計(jì)的常數(shù)。

    4 糧油食品中的微生物預(yù)測(cè)模型的應(yīng)用

    現(xiàn)階段對(duì)于模型的應(yīng)用大多集中在產(chǎn)品的儲(chǔ)藏階段:真菌的動(dòng)力學(xué)模型、外界條件影響的二級(jí)模型、毒素積累模型以及產(chǎn)毒概率模型。目前人工模擬的儲(chǔ)藏實(shí)驗(yàn)一般通過(guò)控制糧油食品儲(chǔ)藏環(huán)境中的溫度和濕度來(lái)控制微生物的生長(zhǎng)以及毒素的產(chǎn)生。通過(guò)人工模擬儲(chǔ)藏實(shí)驗(yàn)?zāi)軐?duì)實(shí)際微生物的生長(zhǎng)和產(chǎn)毒的情況進(jìn)行準(zhǔn)確預(yù)測(cè),以達(dá)到預(yù)警的目的。對(duì)于產(chǎn)品的加工階段由于沒(méi)有合適的可控外界條件,因此目前加工階段關(guān)于模型的研究還較少。

    4.1 儲(chǔ)藏階段動(dòng)力學(xué)模型

    小麥在儲(chǔ)藏中更易滋長(zhǎng)霉菌,最容易感染嘔吐毒素(DON)和玉米赤霉烯酮(ZEN)等真菌毒素[41]。彭堅(jiān)等[42]以小麥為基質(zhì)建立了黃曲霉在小麥貯藏階段的初級(jí)和二級(jí)模型,分別利用Gompertz模型、Logistic模型對(duì)黃曲霉在小麥中生長(zhǎng)曲線(xiàn)進(jìn)行擬合,通過(guò)t檢驗(yàn)和F檢驗(yàn)發(fā)現(xiàn)Gompertz模型擬合效果更好,符合黃曲霉生長(zhǎng)規(guī)律。利用Ratkowsky方程擬合不同溫度下黃曲霉生長(zhǎng)二級(jí)模型發(fā)現(xiàn)先R2相關(guān)性比較好, 達(dá)到0.99以上。

    玉米是我國(guó)三大糧油食品品種之一,其在儲(chǔ)藏期間極易遭受灰綠曲霉和黃曲霉的污染[43],短期內(nèi)可能發(fā)生霉腐變質(zhì)現(xiàn)象。岳曉禹等[44]進(jìn)行溫度對(duì)儲(chǔ)藏玉米中霉菌生長(zhǎng)動(dòng)力學(xué)模型構(gòu)建,以?xún)?chǔ)藏玉米霉菌為研究對(duì)象,在不同儲(chǔ)藏溫度、不同時(shí)間條件下用修正的Gompertz方程擬合腐敗霉菌的生長(zhǎng)動(dòng)力學(xué)模型,采用修正的線(xiàn)性Arrhenius-Davey方程和二階多項(xiàng)式方程擬合溫度對(duì)玉米中霉菌生長(zhǎng)的比生長(zhǎng)速率和遲滯期的影響,通過(guò)模型驗(yàn)證表明修正的線(xiàn)性Arrhenius-Davey方程可以較好地描述、預(yù)測(cè)不同溫度對(duì)玉米中霉菌的比生長(zhǎng)速率和遲滯期的影響。

    房保海等[45]在不同溫度和相對(duì)濕度條件下對(duì)花生中的黃曲霉菌株生長(zhǎng)特性進(jìn)行測(cè)定,并采用Boltzmann 和Logistic模型對(duì)其生長(zhǎng)曲線(xiàn)進(jìn)行擬合,結(jié)果表明在適宜溫濕度條件花生黃曲霉生長(zhǎng)曲線(xiàn)更適合用Boltzmann 模型擬合。

    Yue等[46]在大米干燥過(guò)程對(duì)大米黃曲霉增長(zhǎng)率和遲滯期影響的建模研究中,利用Baranyi & Roberts模型較好地?cái)M合了黃曲霉在大米中的生長(zhǎng)曲線(xiàn),利用線(xiàn)性Arrhenius-Davey模型描述不同水分活度下黃曲霉生長(zhǎng)速率和遲滯期,通過(guò)模型驗(yàn)證具有較高的R2和較低的MSE值,模型具有良好的擬合度。利用二階多項(xiàng)式模型來(lái)擬合溫度和水分活度對(duì)生長(zhǎng)率和遲滯期的共同作用,通過(guò)模型成功預(yù)測(cè)出該類(lèi)大米中黃曲霉生長(zhǎng)的最適溫度,對(duì)水稻的貯藏保鮮至關(guān)重要。

    Mousa等[47]以?huà)伖獯竺缀筒诿诪榛|(zhì),通過(guò)控制溫度和水分活度來(lái)對(duì)2種不同基質(zhì)條件黃曲霉徑向生長(zhǎng)速率和黃曲霉毒素的產(chǎn)生分別進(jìn)行建模分析。通過(guò)Baranyi模型較好地?cái)M合菌落半徑隨時(shí)間變化模型,利用線(xiàn)性Arrhenius-Davey模型和多項(xiàng)式模型擬合溫度和aw對(duì)黃曲霉生長(zhǎng)速率的影響,并通過(guò)獨(dú)立的實(shí)驗(yàn)數(shù)據(jù)對(duì)模型加以驗(yàn)證,結(jié)果表明線(xiàn)性Arrhenius-Davey模型是拋光大米和糙米上黃曲霉生長(zhǎng)速率的最佳預(yù)測(cè)模型,預(yù)測(cè)黃曲霉的最適生長(zhǎng)溫度為30 ℃。

    Marin等[48]研究溫度、濕度等外界條件對(duì)開(kāi)心果貯藏的影響,人工模擬不同溫度、水分活度對(duì)開(kāi)心果中的黃曲霉生長(zhǎng)影響進(jìn)行建模分析。利用Davey模型、Gibson模型、一般多項(xiàng)式模型、多因子cardinal模型和改進(jìn)的cardinal模型五種方法評(píng)價(jià)黃曲霉徑向生長(zhǎng)速率μR對(duì)環(huán)境的響應(yīng)值。通過(guò)均方根誤差(RMSE值)判斷發(fā)現(xiàn)cardinal模型對(duì)黃曲霉的徑向生長(zhǎng)速率擬合度較高。同時(shí)利用改進(jìn)的cardinal模型(減少觀測(cè)實(shí)驗(yàn)值數(shù)量)擬合的同時(shí)發(fā)現(xiàn)預(yù)測(cè)值更接近于觀測(cè)值,但因?yàn)榇嬖谄罟烙?jì),r2只有0.67而多因素cardinal模型r2為0.78。改進(jìn)的cardinal模型是基于減少的實(shí)驗(yàn)觀測(cè)點(diǎn)開(kāi)發(fā)的,從而使驗(yàn)證步驟的預(yù)測(cè)煩瑣程度得到改善[49]。

    4.2 菌株生長(zhǎng)邊界模型及毒素積累預(yù)測(cè)模型

    Battey等[50]利用Logistic回歸模型預(yù)測(cè)腐敗霉菌生長(zhǎng)概率,以0、1進(jìn)行賦值。0代表無(wú)霉菌生長(zhǎng),1代表有霉菌生長(zhǎng)。發(fā)現(xiàn)pH、可滴定酸度、糖含量、苯甲酸鈉和山梨酸鉀水平都是影響霉菌隨時(shí)間存在生長(zhǎng)可能性的重要因素。該Logistic回歸模型針對(duì)14種新條件進(jìn)行了驗(yàn)證,并預(yù)測(cè)了8周后霉菌的生長(zhǎng)狀況,準(zhǔn)確度超過(guò)96%。產(chǎn)品開(kāi)發(fā)人員可以使用該模型來(lái)預(yù)測(cè)即飲飲料中的霉菌生長(zhǎng)狀況。

    Marín等[51]將概率模型首次應(yīng)用于霉菌毒素的積累,在不同水分含量和貯存溫度條件下預(yù)測(cè)從開(kāi)心果中分離出的A.carbonarius菌株的生長(zhǎng)和OTA生產(chǎn)邊界。在評(píng)估水分含量對(duì)辣椒中黃曲霉生長(zhǎng)的影響利用logistic回歸模型對(duì)黃曲霉生長(zhǎng)邊界進(jìn)行預(yù)測(cè),結(jié)果表明,將辣椒或者辣椒粉的水分含量保持在30%以下,黃曲霉生長(zhǎng)概率低于50%,辣椒含水量低于1%時(shí),黃曲霉不生長(zhǎng)。

    Lindblad等[52]在溫度10~25 ℃,水分活度0.77~0.95 條件下對(duì)谷物接種疣孢青霉孢子展開(kāi)谷物貯藏實(shí)驗(yàn)。通過(guò)檢測(cè)青霉菌代謝產(chǎn)物OTA的水平是否超過(guò)法定限制5 mg/kg同時(shí)根據(jù)檢測(cè)值將疣孢青霉生長(zhǎng)觀察值定值0或1。即0代表不生長(zhǎng),1代表生長(zhǎng)。在不同的溫度和水分活度條件下,谷物疣孢青霉菌落達(dá)到一定數(shù)量時(shí),利用Logistic 回歸模型計(jì)算超過(guò)法定限度的概率,通過(guò)模型有效地將谷物霉菌生長(zhǎng)與產(chǎn)毒結(jié)合起來(lái)。

    5 問(wèn)題及應(yīng)用展望

    5.1 存在的問(wèn)題

    糧油食品供應(yīng)鏈鏈條長(zhǎng)環(huán)節(jié)多,危害物種類(lèi)多,產(chǎn)生機(jī)理復(fù)雜。霉菌毒素可以在田間以及倉(cāng)庫(kù)貯藏期間產(chǎn)生,并且不同的氣象條件、種植環(huán)境以及農(nóng)藝技術(shù)等都會(huì)影響霉菌的生長(zhǎng)及其次級(jí)代謝產(chǎn)物的產(chǎn)生。由于其危害物產(chǎn)生機(jī)理的復(fù)雜性,源頭控制危害物難度較大,因此宜采用預(yù)測(cè)模型的方式對(duì)關(guān)鍵環(huán)節(jié)如貯藏階段、加工環(huán)節(jié)等進(jìn)行優(yōu)化。

    而目前對(duì)于預(yù)測(cè)模型研究存在的主要問(wèn)題是模型數(shù)據(jù)源不完整,覆蓋不全面。以小麥真菌毒素模型構(gòu)建為例,模型的構(gòu)建不可能針對(duì)于所有的小麥品種,只能針對(duì)于特定區(qū)域小麥品種進(jìn)行分析。同時(shí)預(yù)測(cè)模型的可信度和精準(zhǔn)度也存在風(fēng)險(xiǎn)?,F(xiàn)今大多數(shù)微生物模型的擬合都集中在實(shí)驗(yàn)室基礎(chǔ)上進(jìn)行仿真模擬,針對(duì)相應(yīng)食品基質(zhì)或者對(duì)應(yīng)基質(zhì)培養(yǎng)基來(lái)進(jìn)行,但是經(jīng)研究發(fā)現(xiàn)兩者擬合出的數(shù)據(jù)存在較大差別,后者擬合模型普遍存在可信度低的問(wèn)題。且實(shí)驗(yàn)室基礎(chǔ)的模型擬合只能是在有限的條件下,基于所獲得的有限數(shù)據(jù)基礎(chǔ)上的擬合,數(shù)據(jù)的不足是制約模型擬合精準(zhǔn)度的最大難題。目前只能通過(guò)盡可能多的獲取數(shù)據(jù)來(lái)進(jìn)行模型構(gòu)建,但大量數(shù)據(jù)的獲取和數(shù)據(jù)庫(kù)的構(gòu)建本身就是難度很大的挑戰(zhàn)。

    5.2 展望

    目前糧油食品中危害物的預(yù)測(cè)模型還主要是集中在對(duì)危害霉菌的生長(zhǎng)預(yù)測(cè)模型研究,霉菌產(chǎn)毒素積累預(yù)測(cè)模型、產(chǎn)毒概率模型研究還比較少。今后應(yīng)重點(diǎn)關(guān)注毒素預(yù)測(cè)模型,通過(guò)毒素含量及特定條件下的產(chǎn)毒概率直觀地反映食品的安全狀況。現(xiàn)階段對(duì)糧油食品加工過(guò)程中真菌毒素消減規(guī)律及加工階段毒素模型研究還較少。未來(lái)需要對(duì)加工階段的某些關(guān)鍵點(diǎn)進(jìn)行大量模型參數(shù)擬合研究,以通過(guò)模型對(duì)加工階段的關(guān)鍵參數(shù)進(jìn)行控制改進(jìn),嘗試進(jìn)行加工環(huán)節(jié)的預(yù)警溯源分析。同時(shí)模型本身的構(gòu)建基于大量實(shí)驗(yàn)數(shù)據(jù)的獲取,這在實(shí)驗(yàn)室研究基礎(chǔ)上是比較棘手的問(wèn)題。因此只有通過(guò)在糧食食品供應(yīng)鏈各環(huán)節(jié)中大量獲取真實(shí)有效的數(shù)據(jù)才能建立起較完整的糧油產(chǎn)業(yè)供應(yīng)鏈微生物預(yù)測(cè)模型,這對(duì)建立適合我國(guó)國(guó)情的微生物預(yù)報(bào)體系具有極其重要的意義。

    猜你喜歡
    糧油食品黃曲霉霉菌
    揭示水霉菌繁殖和侵染過(guò)程
    IAC-HPLC-ESI-MS/MS法測(cè)定不同產(chǎn)地柏子仁中4種黃曲霉毒素
    中成藥(2021年5期)2021-07-21 08:38:40
    雞黃曲霉毒素中毒的臨床表現(xiàn)、實(shí)驗(yàn)室診斷與防治
    《糧油食品科技》辦刊理念視覺(jué)表達(dá)
    中糧集團(tuán)有限公司:努力打造糧油食品行業(yè)人才培養(yǎng)高地
    歡迎訂閱2019年《糧油食品科技》
    2017 年《糧油食品科技》總目錄
    霉菌的新朋友—地衣
    地衣和霉菌
    青霉菌柚苷酶的分離純化及酶學(xué)性質(zhì)
    中文字幕人妻丝袜一区二区| 国产精品乱码一区二三区的特点 | 熟女少妇亚洲综合色aaa.| 多毛熟女@视频| xxx96com| 水蜜桃什么品种好| 久久人人97超碰香蕉20202| 男女午夜视频在线观看| 国产高清激情床上av| 亚洲男人天堂网一区| 亚洲精品美女久久av网站| 免费观看人在逋| 丝袜人妻中文字幕| 久久久精品国产亚洲av高清涩受| 最新在线观看一区二区三区| 无人区码免费观看不卡| bbb黄色大片| av片东京热男人的天堂| 久久精品国产综合久久久| 欧美日韩亚洲高清精品| 久久久久久大精品| 性少妇av在线| 欧美 亚洲 国产 日韩一| 亚洲伊人色综图| 女生性感内裤真人,穿戴方法视频| 亚洲五月婷婷丁香| 天堂影院成人在线观看| 动漫黄色视频在线观看| 美女高潮喷水抽搐中文字幕| 成人黄色视频免费在线看| 精品第一国产精品| 亚洲精品中文字幕在线视频| 成年人免费黄色播放视频| 制服诱惑二区| 久久青草综合色| 精品国产超薄肉色丝袜足j| 精品国产美女av久久久久小说| 人成视频在线观看免费观看| 国内久久婷婷六月综合欲色啪| 女性被躁到高潮视频| 国产在线观看jvid| 夜夜爽天天搞| 久久久久国产一级毛片高清牌| 精品久久久精品久久久| 99国产极品粉嫩在线观看| videosex国产| 女人被狂操c到高潮| 亚洲午夜理论影院| 无遮挡黄片免费观看| 天天添夜夜摸| 首页视频小说图片口味搜索| 香蕉久久夜色| 一a级毛片在线观看| 久久中文字幕人妻熟女| 美女扒开内裤让男人捅视频| 在线观看66精品国产| 男女做爰动态图高潮gif福利片 | 午夜精品久久久久久毛片777| 9色porny在线观看| 丰满的人妻完整版| 亚洲人成电影观看| 国产高清videossex| 亚洲av片天天在线观看| 丰满人妻熟妇乱又伦精品不卡| 88av欧美| 午夜免费观看网址| 久久亚洲精品不卡| 国产欧美日韩精品亚洲av| 国产精品自产拍在线观看55亚洲| 欧美在线一区亚洲| 欧美乱妇无乱码| 久久草成人影院| 91大片在线观看| 久久久久久久久中文| 国产乱人伦免费视频| 亚洲精品在线观看二区| 久久中文看片网| 亚洲九九香蕉| 多毛熟女@视频| 在线观看免费视频日本深夜| 亚洲熟妇熟女久久| 国产成人欧美在线观看| av免费在线观看网站| av欧美777| a级毛片在线看网站| 国产成人啪精品午夜网站| 欧美日韩精品网址| 99精品欧美一区二区三区四区| 欧美成狂野欧美在线观看| 亚洲欧美日韩高清在线视频| 91字幕亚洲| 国产一区二区三区在线臀色熟女 | 日本免费一区二区三区高清不卡 | 国产国语露脸激情在线看| 国产精品野战在线观看 | 成熟少妇高潮喷水视频| 亚洲精品成人av观看孕妇| 国产一卡二卡三卡精品| 99热国产这里只有精品6| 成人精品一区二区免费| av欧美777| 日韩高清综合在线| 看黄色毛片网站| 亚洲自拍偷在线| а√天堂www在线а√下载| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩精品网址| 国产一区二区激情短视频| 国产成人一区二区三区免费视频网站| 久久精品亚洲熟妇少妇任你| 欧美黄色片欧美黄色片| e午夜精品久久久久久久| 国产免费男女视频| 国产免费av片在线观看野外av| 大香蕉久久成人网| 精品久久久精品久久久| 久久精品国产清高在天天线| 精品国产一区二区久久| 精品福利观看| 视频区欧美日本亚洲| 首页视频小说图片口味搜索| 动漫黄色视频在线观看| 90打野战视频偷拍视频| 亚洲男人的天堂狠狠| 国产精品亚洲一级av第二区| 成人免费观看视频高清| 成年女人毛片免费观看观看9| 老熟妇乱子伦视频在线观看| 亚洲熟妇熟女久久| 国产97色在线日韩免费| 国产视频一区二区在线看| 在线免费观看的www视频| 中文欧美无线码| 黄色毛片三级朝国网站| 男女下面进入的视频免费午夜 | 免费在线观看完整版高清| 国产熟女午夜一区二区三区| 国产精品乱码一区二三区的特点 | 嫁个100分男人电影在线观看| 欧美最黄视频在线播放免费 | 午夜精品久久久久久毛片777| 9191精品国产免费久久| 成人三级黄色视频| 一区在线观看完整版| 亚洲午夜理论影院| 两个人免费观看高清视频| 午夜精品在线福利| 超碰成人久久| 欧美午夜高清在线| 两人在一起打扑克的视频| 亚洲国产精品合色在线| 亚洲成人国产一区在线观看| 性欧美人与动物交配| 国产乱人伦免费视频| 欧美黑人欧美精品刺激| 婷婷丁香在线五月| 久久久久久人人人人人| 夜夜爽天天搞| 91在线观看av| 婷婷六月久久综合丁香| 久久香蕉激情| 欧美日韩精品网址| 亚洲中文av在线| 亚洲第一欧美日韩一区二区三区| 老司机靠b影院| 丝袜美腿诱惑在线| 50天的宝宝边吃奶边哭怎么回事| 无遮挡黄片免费观看| 1024视频免费在线观看| 91在线观看av| 国产精品二区激情视频| 免费久久久久久久精品成人欧美视频| 男人舔女人的私密视频| 99国产精品一区二区蜜桃av| 色综合婷婷激情| 国产主播在线观看一区二区| 国产欧美日韩综合在线一区二区| 精品少妇一区二区三区视频日本电影| 丝袜美腿诱惑在线| 国产99久久九九免费精品| 欧美日本亚洲视频在线播放| 亚洲色图 男人天堂 中文字幕| 精品久久久精品久久久| 亚洲精品美女久久av网站| 老司机午夜福利在线观看视频| 久久精品91蜜桃| 亚洲人成77777在线视频| www.精华液| 久久婷婷成人综合色麻豆| 亚洲狠狠婷婷综合久久图片| 色尼玛亚洲综合影院| 超碰97精品在线观看| 国产免费现黄频在线看| 久久久久精品国产欧美久久久| 日本免费一区二区三区高清不卡 | 高清av免费在线| 国产亚洲精品久久久久5区| 免费久久久久久久精品成人欧美视频| 神马国产精品三级电影在线观看 | www.熟女人妻精品国产| 国产精品一区二区免费欧美| 麻豆一二三区av精品| 激情在线观看视频在线高清| 国产在线观看jvid| 午夜亚洲福利在线播放| 国产成人欧美在线观看| 高清欧美精品videossex| 亚洲自拍偷在线| 欧美日韩亚洲综合一区二区三区_| 极品人妻少妇av视频| 亚洲欧美一区二区三区久久| 韩国精品一区二区三区| 久久香蕉国产精品| 精品一区二区三区四区五区乱码| 美国免费a级毛片| 久久人人精品亚洲av| 性色av乱码一区二区三区2| xxx96com| 两个人免费观看高清视频| 欧美日韩一级在线毛片| 欧美日韩精品网址| 亚洲人成77777在线视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人不卡在线观看播放网| a级毛片在线看网站| 精品卡一卡二卡四卡免费| 婷婷六月久久综合丁香| 精品久久蜜臀av无| 神马国产精品三级电影在线观看 | 韩国精品一区二区三区| 国产精品偷伦视频观看了| 一本综合久久免费| 亚洲精品国产精品久久久不卡| 亚洲avbb在线观看| 亚洲欧美日韩无卡精品| x7x7x7水蜜桃| 国产无遮挡羞羞视频在线观看| 少妇裸体淫交视频免费看高清 | 手机成人av网站| 久久精品91蜜桃| av视频免费观看在线观看| 夜夜爽天天搞| svipshipincom国产片| 国产成人精品久久二区二区免费| 少妇粗大呻吟视频| 88av欧美| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩福利视频一区二区| 亚洲精品久久午夜乱码| 欧美乱色亚洲激情| 又紧又爽又黄一区二区| 国产野战对白在线观看| 在线天堂中文资源库| 男人的好看免费观看在线视频 | 日韩欧美在线二视频| 精品高清国产在线一区| 亚洲精品成人av观看孕妇| 国产在线观看jvid| 久久精品亚洲av国产电影网| 9热在线视频观看99| 亚洲狠狠婷婷综合久久图片| 国产99久久九九免费精品| avwww免费| 别揉我奶头~嗯~啊~动态视频| 老熟妇乱子伦视频在线观看| 久久人妻熟女aⅴ| 精品国内亚洲2022精品成人| 91麻豆精品激情在线观看国产 | 精品国产美女av久久久久小说| 欧美丝袜亚洲另类 | 国产黄a三级三级三级人| svipshipincom国产片| 亚洲美女黄片视频| 久久久精品国产亚洲av高清涩受| 纯流量卡能插随身wifi吗| 18美女黄网站色大片免费观看| netflix在线观看网站| 日韩欧美三级三区| 黑人猛操日本美女一级片| 新久久久久国产一级毛片| 午夜精品久久久久久毛片777| 亚洲精品一卡2卡三卡4卡5卡| 精品国产乱码久久久久久男人| 看免费av毛片| 亚洲男人天堂网一区| 免费久久久久久久精品成人欧美视频| 免费少妇av软件| 一级作爱视频免费观看| 在线观看www视频免费| 亚洲国产欧美一区二区综合| 亚洲av五月六月丁香网| 日本a在线网址| 丝袜美足系列| 女性生殖器流出的白浆| 9色porny在线观看| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 国产又色又爽无遮挡免费看| 可以免费在线观看a视频的电影网站| 久久久久久人人人人人| 老熟妇仑乱视频hdxx| 国产91精品成人一区二区三区| 国产单亲对白刺激| 无遮挡黄片免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 激情视频va一区二区三区| cao死你这个sao货| 又黄又粗又硬又大视频| 久久天堂一区二区三区四区| 免费在线观看黄色视频的| 久久精品亚洲精品国产色婷小说| 国产精品影院久久| 满18在线观看网站| 国产精品久久电影中文字幕| 国产视频一区二区在线看| av有码第一页| 欧美成人免费av一区二区三区| 国产亚洲精品久久久久久毛片| 亚洲精品国产精品久久久不卡| 亚洲一区高清亚洲精品| 老司机在亚洲福利影院| 国产免费av片在线观看野外av| 在线观看日韩欧美| 免费在线观看视频国产中文字幕亚洲| 久久精品aⅴ一区二区三区四区| 欧美成狂野欧美在线观看| 一本综合久久免费| 另类亚洲欧美激情| 久久这里只有精品19| 变态另类成人亚洲欧美熟女 | 国产精品av久久久久免费| xxxhd国产人妻xxx| 啦啦啦 在线观看视频| 我的亚洲天堂| 国产精品一区二区三区四区久久 | 午夜精品在线福利| 老汉色∧v一级毛片| 午夜免费激情av| 国产精品1区2区在线观看.| 久久99一区二区三区| 国产免费男女视频| 久久精品国产亚洲av高清一级| 亚洲一区二区三区不卡视频| 国产亚洲精品久久久久5区| 欧美日韩国产mv在线观看视频| 女生性感内裤真人,穿戴方法视频| 我的亚洲天堂| 日韩三级视频一区二区三区| 黄色毛片三级朝国网站| videosex国产| 91在线观看av| 欧美激情高清一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 高清av免费在线| 日本欧美视频一区| 人妻久久中文字幕网| 国产又爽黄色视频| 亚洲成国产人片在线观看| 丝袜在线中文字幕| 久久天躁狠狠躁夜夜2o2o| av电影中文网址| 一级毛片高清免费大全| 国产一区二区三区视频了| 夜夜爽天天搞| 热99国产精品久久久久久7| 十八禁网站免费在线| 国产麻豆69| 黄网站色视频无遮挡免费观看| 亚洲av熟女| 黄色毛片三级朝国网站| 国产精品香港三级国产av潘金莲| 午夜老司机福利片| 亚洲精品一卡2卡三卡4卡5卡| 成在线人永久免费视频| 夜夜躁狠狠躁天天躁| 国产蜜桃级精品一区二区三区| 精品国产一区二区久久| 久久香蕉精品热| 日日干狠狠操夜夜爽| 看免费av毛片| av在线播放免费不卡| 国产97色在线日韩免费| 亚洲一区二区三区色噜噜 | 夜夜爽天天搞| 精品少妇一区二区三区视频日本电影| 中文欧美无线码| 色哟哟哟哟哟哟| 久久午夜亚洲精品久久| 男人舔女人的私密视频| 夜夜夜夜夜久久久久| 久久久久久人人人人人| 久久精品91无色码中文字幕| 看黄色毛片网站| 亚洲精品久久成人aⅴ小说| 99久久国产精品久久久| 9热在线视频观看99| 19禁男女啪啪无遮挡网站| 人人澡人人妻人| 十八禁网站免费在线| 国产精品电影一区二区三区| 成人亚洲精品一区在线观看| 久久精品亚洲精品国产色婷小说| 亚洲aⅴ乱码一区二区在线播放 | 女警被强在线播放| www.自偷自拍.com| 少妇的丰满在线观看| 国产欧美日韩一区二区三区在线| 高清黄色对白视频在线免费看| 国产在线精品亚洲第一网站| 99re在线观看精品视频| 少妇的丰满在线观看| 亚洲一区二区三区色噜噜 | 国产高清国产精品国产三级| 国产真人三级小视频在线观看| 国产麻豆69| 麻豆久久精品国产亚洲av | 又黄又爽又免费观看的视频| 中文字幕人妻丝袜一区二区| 国产精品一区二区在线不卡| 日本黄色日本黄色录像| 丝袜美足系列| 最近最新中文字幕大全免费视频| 亚洲精品一卡2卡三卡4卡5卡| 午夜免费激情av| 亚洲成人久久性| 很黄的视频免费| 亚洲美女黄片视频| √禁漫天堂资源中文www| 国产精品国产av在线观看| 黑人欧美特级aaaaaa片| 91精品三级在线观看| 国产高清videossex| 涩涩av久久男人的天堂| 午夜福利在线观看吧| 在线国产一区二区在线| 午夜免费观看网址| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美一区视频在线观看| 黄色a级毛片大全视频| 99在线人妻在线中文字幕| 看黄色毛片网站| 精品国产美女av久久久久小说| 人人妻人人澡人人看| 又紧又爽又黄一区二区| 中文字幕人妻丝袜制服| 亚洲av片天天在线观看| 亚洲精品中文字幕一二三四区| av免费在线观看网站| 久久久久九九精品影院| 纯流量卡能插随身wifi吗| 桃色一区二区三区在线观看| 国产高清激情床上av| 久久人人爽av亚洲精品天堂| 久久久久久大精品| 久久久久国产精品人妻aⅴ院| 这个男人来自地球电影免费观看| 国产欧美日韩一区二区三区在线| 少妇的丰满在线观看| 黄色女人牲交| 亚洲精品av麻豆狂野| 国产av精品麻豆| 亚洲午夜精品一区,二区,三区| 久久久久九九精品影院| 97碰自拍视频| 亚洲国产精品一区二区三区在线| 亚洲精品一卡2卡三卡4卡5卡| 欧美成狂野欧美在线观看| 999久久久国产精品视频| 亚洲成国产人片在线观看| 男女做爰动态图高潮gif福利片 | 免费人成视频x8x8入口观看| 欧美乱妇无乱码| 岛国在线观看网站| 久久久国产精品麻豆| 我的亚洲天堂| 亚洲激情在线av| 国产亚洲欧美98| www.自偷自拍.com| 成人黄色视频免费在线看| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av高清一级| 搡老熟女国产l中国老女人| 长腿黑丝高跟| 黄色丝袜av网址大全| 午夜a级毛片| 亚洲第一欧美日韩一区二区三区| 婷婷精品国产亚洲av在线| 男人操女人黄网站| 亚洲三区欧美一区| 久久亚洲真实| 手机成人av网站| 中文字幕人妻丝袜制服| 国产精品野战在线观看 | 美女扒开内裤让男人捅视频| 国产精品一区二区三区四区久久 | 老司机午夜福利在线观看视频| 1024香蕉在线观看| 色播在线永久视频| 色精品久久人妻99蜜桃| 久久久久久人人人人人| 成熟少妇高潮喷水视频| 亚洲精品美女久久av网站| 老司机深夜福利视频在线观看| www.www免费av| 国产深夜福利视频在线观看| 美女扒开内裤让男人捅视频| 亚洲精品粉嫩美女一区| 岛国在线观看网站| 欧美人与性动交α欧美精品济南到| 波多野结衣av一区二区av| 国产av一区在线观看免费| 色综合站精品国产| 欧美不卡视频在线免费观看 | 成人三级黄色视频| 久久久久久久久免费视频了| 大陆偷拍与自拍| 一级a爱视频在线免费观看| 女警被强在线播放| 久久亚洲精品不卡| 真人做人爱边吃奶动态| 老鸭窝网址在线观看| 午夜精品国产一区二区电影| 校园春色视频在线观看| 久久香蕉国产精品| 亚洲色图av天堂| 午夜精品久久久久久毛片777| 免费少妇av软件| 人人妻人人添人人爽欧美一区卜| 少妇裸体淫交视频免费看高清 | 国产精品久久视频播放| 亚洲专区中文字幕在线| 夫妻午夜视频| 淫秽高清视频在线观看| 女同久久另类99精品国产91| 极品人妻少妇av视频| 精品国产美女av久久久久小说| 亚洲av片天天在线观看| 18禁裸乳无遮挡免费网站照片 | 不卡av一区二区三区| 如日韩欧美国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费鲁丝| 黄频高清免费视频| 免费在线观看影片大全网站| 熟女少妇亚洲综合色aaa.| 久久香蕉国产精品| 水蜜桃什么品种好| 精品一区二区三区四区五区乱码| 亚洲视频免费观看视频| 欧美日韩国产mv在线观看视频| 亚洲人成电影观看| 精品久久久久久,| 夜夜躁狠狠躁天天躁| 欧美成人性av电影在线观看| 国产欧美日韩一区二区三区在线| 夫妻午夜视频| 99精国产麻豆久久婷婷| 国产99白浆流出| 欧美成人免费av一区二区三区| 亚洲精品粉嫩美女一区| 亚洲精品一卡2卡三卡4卡5卡| 精品无人区乱码1区二区| 国产一卡二卡三卡精品| 无限看片的www在线观看| 真人一进一出gif抽搐免费| 精品国产一区二区三区四区第35| 国产1区2区3区精品| 亚洲一区二区三区欧美精品| 午夜影院日韩av| 高潮久久久久久久久久久不卡| 99国产精品一区二区蜜桃av| 97碰自拍视频| 国产日韩一区二区三区精品不卡| 18禁美女被吸乳视频| 亚洲国产欧美日韩在线播放| 午夜成年电影在线免费观看| 露出奶头的视频| 日本黄色视频三级网站网址| 精品日产1卡2卡| av有码第一页| 欧美成人午夜精品| 国产成人欧美在线观看| 国产黄a三级三级三级人| 黄色成人免费大全| 免费不卡黄色视频| 91精品三级在线观看| 欧美日韩精品网址| 精品一区二区三区四区五区乱码| 免费在线观看完整版高清| 女性被躁到高潮视频| 久久久精品国产亚洲av高清涩受| 黄片小视频在线播放| 日韩国内少妇激情av| 国产欧美日韩综合在线一区二区| 欧美午夜高清在线| 日韩国内少妇激情av| 18禁黄网站禁片午夜丰满| 国产91精品成人一区二区三区| 欧美性长视频在线观看| 日韩有码中文字幕| 亚洲成人国产一区在线观看| 免费日韩欧美在线观看| 美女 人体艺术 gogo| 国产成年人精品一区二区 | 亚洲人成电影免费在线| 亚洲激情在线av| 老司机亚洲免费影院| 精品无人区乱码1区二区| 欧美黑人欧美精品刺激| 亚洲色图av天堂| 在线av久久热| av有码第一页| 亚洲 欧美 日韩 在线 免费| 免费久久久久久久精品成人欧美视频| 极品教师在线免费播放|