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

    基于PCA和ARIMA的食品安全風(fēng)險(xiǎn)因素辨識(shí)及預(yù)測(cè)方法

    2021-11-04 12:35:58沈霞芬蔡強(qiáng)吳利芳肖革新
    食品工業(yè) 2021年10期
    關(guān)鍵詞:起數(shù)預(yù)測(cè)因素

    沈霞芬,蔡強(qiáng)*,吳利芳,肖革新

    1.浙江清華長(zhǎng)三角研究院(嘉興 314006);2.國(guó)家食品安全風(fēng)險(xiǎn)評(píng)估中心(北京 100022)

    食品安全是個(gè)綜合概念,涉及衛(wèi)生、質(zhì)量、營(yíng)養(yǎng)等多方面內(nèi)容,以及種植、養(yǎng)殖、生產(chǎn)、包裝、儲(chǔ)存、運(yùn)輸、流通、消費(fèi)等多個(gè)環(huán)節(jié),單一的食品衛(wèi)生、食品質(zhì)量、食品營(yíng)養(yǎng)等概念無(wú)法涵蓋以上全部?jī)?nèi)容[1]。因此,關(guān)于食品安全,除了食品生產(chǎn)過(guò)程中所有指標(biāo)需要滿足特定的標(biāo)準(zhǔn)外,許多食品安全風(fēng)險(xiǎn)的存在和發(fā)展由食品供應(yīng)鏈其他因素驅(qū)動(dòng),如空氣、土壤、水等農(nóng)業(yè)資源環(huán)境的影響,種養(yǎng)殖過(guò)程中化肥、農(nóng)獸藥等施用帶來(lái)的影響,人體對(duì)食物的攝入量等[2]。這些因素與供應(yīng)鏈之間的相互作用復(fù)雜,需要使用一種系統(tǒng)或整體的分析方法揭示其中的關(guān)系。目前,國(guó)內(nèi)在食品安全風(fēng)險(xiǎn)領(lǐng)域的分析,主要以監(jiān)管部門監(jiān)督抽檢和日常檢查數(shù)據(jù)的風(fēng)險(xiǎn)分析為主,食品加工流通環(huán)節(jié)外其他因素的分析尚未形成深入挖掘。黃湘鷺等[3]對(duì)2016—2017年國(guó)家食品安全監(jiān)督抽檢的不合格項(xiàng)目進(jìn)行歸類分析,結(jié)果發(fā)現(xiàn)抽檢合格率呈現(xiàn)穩(wěn)中有升的趨勢(shì)。李宗亮[4]以食品的污染物和營(yíng)養(yǎng)物質(zhì)為評(píng)價(jià)對(duì)象,以不合格率和不合格程度為評(píng)價(jià)指標(biāo)構(gòu)建食品風(fēng)險(xiǎn)系數(shù)評(píng)價(jià)體系,并結(jié)合BP神經(jīng)網(wǎng)絡(luò)模型,將食品風(fēng)險(xiǎn)系數(shù)應(yīng)用于各類食品及抽檢時(shí)待測(cè)參數(shù)的安全趨勢(shì)預(yù)測(cè)預(yù)警等領(lǐng)域,為監(jiān)管部門提供決策依據(jù)。相比之下,歐美發(fā)達(dá)國(guó)家的政府信息公開(kāi)工作具有起步早、數(shù)據(jù)多、深度開(kāi)放等特點(diǎn),在食品安全分析方面也有著更進(jìn)一步的研究[5]。Yamine等[6]使用貝葉斯網(wǎng)絡(luò)模型來(lái)識(shí)別和量化RASFF中關(guān)于源自印度、土耳其和荷蘭的水果和蔬菜的所有食品安全通告數(shù)據(jù)與產(chǎn)品原產(chǎn)國(guó)氣候因素、農(nóng)業(yè)因素和經(jīng)濟(jì)因素之間的關(guān)系,研究結(jié)果表明不能確定氣候因素對(duì)蔬果化學(xué)食品安全危害產(chǎn)生影響。Alberto等[7]根據(jù)聯(lián)合國(guó)糧食及農(nóng)業(yè)組織提供的關(guān)于進(jìn)出口和生產(chǎn)的公開(kāi)數(shù)據(jù),構(gòu)建數(shù)學(xué)模型模擬糧食產(chǎn)品的分配情況,首次從食品安全風(fēng)險(xiǎn)的角度論證網(wǎng)格結(jié)構(gòu)對(duì)食品分布影響。

    從食品安全多個(gè)影響因素的綜合分析與關(guān)鍵風(fēng)險(xiǎn)因子確定,以及單一指標(biāo)數(shù)據(jù)的預(yù)測(cè)2個(gè)角度探討食品安全風(fēng)險(xiǎn)分析的數(shù)據(jù)模型,并給出相應(yīng)的應(yīng)用實(shí)例。

    1 材料與方法

    1.1 材料

    1.1.1 數(shù)據(jù)來(lái)源

    食品安全涉及從農(nóng)田到餐桌的整個(gè)過(guò)程,理論上需要收集覆蓋食品供應(yīng)鏈上的所有數(shù)據(jù),同時(shí)融合時(shí)空、地理、環(huán)境、氣象及醫(yī)療健康等信息。但國(guó)內(nèi)目前尚未建立統(tǒng)一的數(shù)據(jù)標(biāo)準(zhǔn),也尚未解決跨層級(jí)、跨部門、跨業(yè)務(wù)、全鏈條多源異構(gòu)數(shù)據(jù)的融合,所以數(shù)據(jù)收集十分困難。為既兼顧各指標(biāo)的合理性、時(shí)間上的統(tǒng)一又不失數(shù)據(jù)的準(zhǔn)確性,從國(guó)家統(tǒng)計(jì)局官網(wǎng)上收集1999—2018年關(guān)于中國(guó)環(huán)境、農(nóng)業(yè)因素及消費(fèi)量的6個(gè)食品安全評(píng)價(jià)指標(biāo)數(shù)據(jù),數(shù)據(jù)維度為6×20。同時(shí),以“食物中毒事件”為關(guān)鍵詞,通過(guò)國(guó)家衛(wèi)生健康委員會(huì)官網(wǎng)(http://www.nhc.gov.cn)檢索的“衛(wèi)生部辦公廳關(guān)于XXXX年全國(guó)食物中毒事件情況的通報(bào)”文件中,收集2008—2015年公開(kāi)發(fā)布的中國(guó)食物中毒事件起數(shù)的月度數(shù)據(jù),數(shù)據(jù)長(zhǎng)度為96條。

    1.1.2 分析軟件

    采用R對(duì)食品安全6個(gè)影響因素進(jìn)行主成分分析,并對(duì)食物中毒事件起數(shù)季度數(shù)據(jù)進(jìn)行時(shí)間序列分析。以α=0.05為檢驗(yàn)水準(zhǔn),p<0.05說(shuō)明具有顯著性差異。

    1.2 方法

    1.2.1 多影響因素的綜合分析

    采用主成分分析法從多個(gè)影響因素入手,提取主要成分,并基于成分的方差貢獻(xiàn)率計(jì)算綜合評(píng)分,以此作為各影響因素對(duì)年度食品安全的影響程度。采用灰色關(guān)聯(lián)度分析,計(jì)算食品安全影響因素與綜合評(píng)分之間的關(guān)聯(lián)度,確定最為關(guān)鍵的要素。

    1) 構(gòu)建一個(gè)m×n的矩陣X,行表示年份,列表示食品安全的影響因素。

    2) 采用Z-score對(duì)以上n個(gè)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化使其量綱統(tǒng)一,得到1組新變量(z1,z2,…,zn)[8],并計(jì)算這組新變量的相關(guān)系數(shù)矩陣R。

    3) 由相關(guān)系數(shù)矩陣R計(jì)算得到其特征值λ,并從大到小排序,λ1≥λ2≥…≥λn≥0,且計(jì)算各特征值對(duì)應(yīng)的特征向量εi=(εi1,εi2,…,εin),i=1,2,…,n,確定X的第i個(gè)主成分。

    4) 采用灰色關(guān)聯(lián)度計(jì)算各評(píng)價(jià)指標(biāo)Xi與綜合評(píng)分Y之間關(guān)聯(lián)度[9],關(guān)聯(lián)系數(shù)定義為:

    式中:i為第i個(gè)評(píng)價(jià)指標(biāo);k為各指標(biāo)的第k個(gè)數(shù)值,故Y(k)為綜合評(píng)分Y的第k個(gè)數(shù)值,Xi(k)為第i個(gè)評(píng)價(jià)指標(biāo)的第k個(gè)數(shù)值;ρ為分辨系數(shù),取值為(0,1),ρ越小,區(qū)分度越低,一般取值為0.5。

    由于關(guān)聯(lián)系數(shù)是比較數(shù)列與參考數(shù)列在各個(gè)時(shí)刻的關(guān)聯(lián)程度值,所以它是一個(gè)序列,信息過(guò)于分散不便于進(jìn)行整體比較。因此需要將各時(shí)刻的關(guān)聯(lián)系數(shù)求其平均值,作為比較數(shù)列xi與參考數(shù)列Y之間關(guān)聯(lián)度的數(shù)量表示,關(guān)聯(lián)度ri定義為:

    ri越大,表明指標(biāo)i對(duì)綜合評(píng)分的影響越大。

    1.2.2 單一指標(biāo)數(shù)據(jù)的預(yù)測(cè)分析

    對(duì)單位時(shí)間內(nèi)某個(gè)指標(biāo)數(shù)據(jù)進(jìn)行統(tǒng)計(jì),如在農(nóng)業(yè)生產(chǎn)中,農(nóng)獸藥和化肥的使用量可按季度、按月記錄;加工運(yùn)輸環(huán)節(jié)中,食品生產(chǎn)件數(shù)、檢驗(yàn)檢測(cè)信息等可按周、按天記錄,這些數(shù)據(jù)便為時(shí)間序列。對(duì)于數(shù)據(jù)都為正值的時(shí)間序列,可通過(guò)5個(gè)步驟完成模型構(gòu)建及分析。

    1) 對(duì)于非線性的變量進(jìn)行對(duì)數(shù)變換,消除數(shù)據(jù)的非線性,且不改變數(shù)據(jù)之間的統(tǒng)計(jì)性質(zhì)。

    2) 對(duì)于非平穩(wěn)的變量進(jìn)行差分變換,弱化隨機(jī)性使其平穩(wěn)化。利用ADF單位根檢驗(yàn),檢驗(yàn)變換后的序列是否平穩(wěn)。具體方法:原假設(shè)H0,序列非平穩(wěn);備擇假設(shè)H1,序列平穩(wěn)。檢驗(yàn)統(tǒng)計(jì)量小于5%臨界值時(shí)拒絕原假設(shè),即認(rèn)為該序列平穩(wěn)。

    3) 對(duì)d階差分后的平穩(wěn)時(shí)間序列構(gòu)建ARIMA模型,若模型殘差的Ljung-Box統(tǒng)計(jì)量的p值大于0.05,則模型通過(guò)檢驗(yàn)[10],再利用所構(gòu)建的模型進(jìn)行預(yù)測(cè)。

    2 結(jié)果與分析

    2.1 食品安全綜合評(píng)價(jià)及關(guān)鍵影響因素確定

    2.1.1 食品安全風(fēng)險(xiǎn)的影響因素

    選取的食品安全風(fēng)險(xiǎn)指標(biāo)包含1999—2018年環(huán)境、農(nóng)業(yè)因素和消費(fèi)量等數(shù)據(jù)。食品安全風(fēng)險(xiǎn)與各影響因素值的大小呈正比。

    2.1.2 食品安全風(fēng)險(xiǎn)的綜合評(píng)分

    從標(biāo)準(zhǔn)化后各影響因素的散點(diǎn)圖矩陣來(lái)看,除了廢氣排放量與其他幾個(gè)指標(biāo)的相關(guān)性較弱之外(各點(diǎn)越接近矩形的對(duì)角線表明相關(guān)性越強(qiáng)),其余指標(biāo)之間都具有一定的相關(guān)性,存在一定的信息重疊(圖1)。主成分分析選取的前2個(gè)主成分的特征值分別為3.34和2.27,前2個(gè)主成分的累計(jì)貢獻(xiàn)率達(dá)93%,基本包括原來(lái)那些變量所包含信息。

    圖1 食品安全影響因素的散點(diǎn)圖矩陣

    由2個(gè)主成分加權(quán)所得的綜合得分,代表著環(huán)境、農(nóng)業(yè)和消費(fèi)因素對(duì)食品安全的影響程度(表2)。

    表2 1999—2018年食品安全影響因素的綜合評(píng)分

    結(jié)果表明,環(huán)境、農(nóng)業(yè)及消費(fèi)量對(duì)中國(guó)食品安全風(fēng)險(xiǎn)的影響具有階段性:1999—2002年中國(guó)食品安全風(fēng)險(xiǎn)逐年降低;2002—2012年食品安全風(fēng)險(xiǎn)基本上呈現(xiàn)逐年升高趨勢(shì);直到2012年之后又基本上恢復(fù)到逐年降低趨勢(shì)。產(chǎn)生該結(jié)果的主要原因在于工業(yè)與經(jīng)濟(jì)的發(fā)展加速了環(huán)境污染及農(nóng)業(yè)生產(chǎn)過(guò)程中農(nóng)獸藥與化肥的使用量,導(dǎo)致2002—2012年食品安全的風(fēng)險(xiǎn)逐漸增高,但以犧牲環(huán)境與國(guó)民健康的工業(yè)和經(jīng)濟(jì)發(fā)展打破中國(guó)可持續(xù)發(fā)展的基本理念?!笆舜蟆敝笳_(kāi)始出臺(tái)各種環(huán)保與民生健康的政策[11-12],使得工業(yè)廢棄物及農(nóng)業(yè)生產(chǎn)過(guò)程中化學(xué)品的投入量有所控制,因此食品安全風(fēng)險(xiǎn)在2012年后呈現(xiàn)明顯降低趨勢(shì)。

    表1 食品安全風(fēng)險(xiǎn)因素

    2.1.3 食品安全關(guān)鍵影響要素確定

    灰色關(guān)聯(lián)度分析得到各影響因素與食品安全風(fēng)險(xiǎn)綜合評(píng)分之間的關(guān)聯(lián)度分別為0.637,0.663,0.682,0.741,0.714和0.598。因此1999—2018年影響中國(guó)食品安全的因素排名依次為單位面積農(nóng)獸藥使用量、單位面積農(nóng)用化肥施用折純量、廢氣(二氧化硫+煙粉塵)排放量、廢水排放量、工業(yè)固體廢物傾倒丟棄量、居民主要食品人均消費(fèi)量。該結(jié)果表明,單從環(huán)境、農(nóng)業(yè)和消費(fèi)層面來(lái)看,農(nóng)業(yè)生產(chǎn)過(guò)程中投入的化學(xué)品量對(duì)于中國(guó)食品安全來(lái)說(shuō)是最為顯著的影響因素,其次為環(huán)境因素。為減輕病蟲(chóng)草害對(duì)產(chǎn)量的威脅并獲得高產(chǎn),農(nóng)民對(duì)于農(nóng)獸藥和化肥不合理使用的行為比較普遍。陳曉明等[13]表明1991—2013年中國(guó)農(nóng)藥使用量增長(zhǎng)135.5%,不合理的使用行為導(dǎo)致其殘留物在食物鏈中不斷傳遞和遷移[13]。史常亮等[14]指出中國(guó)每公頃耕地化肥施用量由1980年的127.8 kg增加到2013年的437.4 kg,是國(guó)際公認(rèn)化肥施用安全上限225 kg/hm2的1.94倍。過(guò)度使用化肥產(chǎn)生的污染土壤的重金屬和放射性物質(zhì),通過(guò)食物鏈不斷在植物和生物體內(nèi)富集,最終在人體內(nèi)積累危害人體健康。

    2.2 時(shí)間序列分析預(yù)測(cè)食品安全風(fēng)險(xiǎn)

    2.2.1 季節(jié)性分布結(jié)果

    2008—2015年全國(guó)食物中毒事件起數(shù)總體上呈現(xiàn)先降低后增高趨勢(shì)(圖2),并且存在明顯季度模式。第3季度的食物中毒事件呈現(xiàn)年度最高,第2季度次之,第1季度最低。因此,以季度作為時(shí)間單位,對(duì)食物中毒事件做時(shí)間序列預(yù)測(cè)將更具準(zhǔn)確性。

    圖2 2008—2015年全國(guó)食物中毒事件起數(shù)的月度數(shù)據(jù)時(shí)序圖

    2.2.2 平穩(wěn)性檢驗(yàn)結(jié)果

    消除非線性后的季度食物中毒事件序列,經(jīng)過(guò)1階差分后還是存在季節(jié)與非季節(jié)相關(guān)性(圖3a2);經(jīng)過(guò)1階季節(jié)差分后的序列消除相關(guān)性(圖3a3);經(jīng)過(guò)1階差分及1階季節(jié)差分后的序列在滯后1階和5階上還是存在一定程度相關(guān)性(圖3a4)。

    圖3 食物中毒事件對(duì)數(shù)序列及其差分后序列的自相關(guān)(ACF,a)和偏自相關(guān)(PACF,b)函數(shù)圖

    平穩(wěn)性檢驗(yàn)結(jié)果表明,除了食物中毒事件起數(shù)的對(duì)數(shù)序列外,其他經(jīng)過(guò)差分處理后序列都平穩(wěn)(p<0.05),符合ARIMA模型建立的條件(表3)。

    表3 所有序列的平穩(wěn)性檢驗(yàn)結(jié)果

    2.2.3 ARIMA模型及預(yù)測(cè)結(jié)果

    根據(jù)差分后序列的自相關(guān)(ACF)和偏自相關(guān)(PACF)函數(shù)圖(圖3)建立3個(gè)ARIMA模型(表4),模型殘差的Ljung-Box統(tǒng)計(jì)量的p值及殘差檢驗(yàn)圖均表明3個(gè)模型通過(guò)檢驗(yàn)。

    表4 根據(jù)自相關(guān)(ACF)和偏自相關(guān)(PACF)圖確定ARIMA模型的擬合優(yōu)度指標(biāo)

    表5為由R自動(dòng)選取的ARIMA模型。其中,ARIMA (3,1,0)等同于表4中第一個(gè)模型,但由于調(diào)用的是R中的不同函數(shù),所以在結(jié)果上略顯差異;考慮季節(jié)性的模型,R自動(dòng)選擇的是ARIMA(0,1,1)×(0,1,0)4,雖然模型殘差的Ljung-Box統(tǒng)計(jì)量的p值大于0.05,但相比于其他幾個(gè)模型的p值小很多,且殘差檢驗(yàn)圖顯示,殘差在滯后5階及其后面的階數(shù)上存在較小相關(guān)性。

    表5 R自動(dòng)選擇ARIMA模型的擬合優(yōu)度指標(biāo)

    以AIC最小原則,選擇最好解釋數(shù)據(jù)但含最少自由參數(shù)的模型ARIMA(5,0,0)×(0,1,0)4,即:

    作為全國(guó)食物中毒事件起數(shù)的時(shí)間序列模型。模型的擬合值與真實(shí)值的變動(dòng)具有較高的一致性(圖4a)。該模型對(duì)2015年度中毒事件起數(shù)的預(yù)測(cè)顯示出較強(qiáng)的季節(jié)性,預(yù)測(cè)2015年第3季度的中毒事件數(shù)仍然是全年最高,達(dá)到66件;第1季度和第4季度的中毒事件起數(shù)相對(duì)較低,分別為16和27件;第2季度的中毒事件起數(shù)全年第二,為35件。產(chǎn)生該結(jié)果的主要原因是春夏秋季節(jié)天氣溫?zé)?、氣候濕?rùn),適宜細(xì)菌和有毒植物的生長(zhǎng)繁殖,成為食物中毒事件的高發(fā)期。從預(yù)測(cè)效果上來(lái)看,該模型對(duì)2015年第2季度的預(yù)測(cè)完全吻合,除了第1季度稍有偏差外,其余3個(gè)季度預(yù)測(cè)效果較好,真實(shí)值都位于95%的預(yù)測(cè)區(qū)間內(nèi)[7](圖4b)。

    圖4 模型ARIMA(5,0,0)×(0,1,0)4的擬合結(jié)果(a)與預(yù)測(cè)結(jié)果(b)

    3 結(jié)論

    食品安全風(fēng)險(xiǎn)分析是一個(gè)針對(duì)復(fù)雜系統(tǒng)的綜合性評(píng)價(jià)問(wèn)題。主成分分析與灰色關(guān)聯(lián)度結(jié)合能夠客觀、全面地反映各項(xiàng)影響因素的信息,同時(shí)通過(guò)估計(jì)綜合評(píng)分與各影響因素之間的一致性程度判斷出最為關(guān)鍵的影響因素。在食品領(lǐng)域,由于氣候變化、飲食習(xí)慣等因素的影響,很多食品安全相關(guān)的數(shù)據(jù)具有周期性。時(shí)間序列ARIMA模型能夠根據(jù)某個(gè)指標(biāo)的歷史數(shù)據(jù),通過(guò)計(jì)量的手段反映時(shí)序中動(dòng)態(tài)依存關(guān)系,不僅從數(shù)量上揭示該指標(biāo)的發(fā)展變化規(guī)律,而且能夠準(zhǔn)確預(yù)測(cè)未來(lái)走勢(shì)。

    然而,在食品領(lǐng)域,食品安全的影響因素不僅是環(huán)境、農(nóng)業(yè)和消費(fèi)因素,還關(guān)系到食品生產(chǎn)和流通過(guò)程中食品質(zhì)量等問(wèn)題,理論上需要實(shí)現(xiàn)對(duì)食品從農(nóng)田到餐桌的全面評(píng)價(jià)。但由于目前國(guó)內(nèi)食品鏈上各環(huán)節(jié)信息化水平不一,一些環(huán)節(jié)數(shù)據(jù)采集缺失,數(shù)據(jù)存儲(chǔ)沒(méi)有統(tǒng)一的標(biāo)準(zhǔn),使得食品安全的綜合評(píng)價(jià)缺乏全面性。但從分析方法上來(lái)說(shuō),主成分分析與時(shí)間序列分析對(duì)于食品安全領(lǐng)域多因素的綜合評(píng)價(jià)和單一指標(biāo)的預(yù)測(cè)分析都具有較好的指示作用。

    猜你喜歡
    起數(shù)預(yù)測(cè)因素
    國(guó)家礦山安全監(jiān)察局:2022年煤礦百萬(wàn)噸死亡率比2012年下降86%
    國(guó)家森防辦:一季度全國(guó)森林草原火災(zāi)起數(shù)較往年顯著下降
    綠葉(2022年4期)2022-11-26 04:20:43
    無(wú)可預(yù)測(cè)
    黃河之聲(2022年10期)2022-09-27 13:59:46
    選修2-2期中考試預(yù)測(cè)卷(B卷)
    選修2-2期中考試預(yù)測(cè)卷(A卷)
    去年全國(guó)生產(chǎn)安全事故起數(shù)下降15.5%
    解石三大因素
    不必預(yù)測(cè)未來(lái),只需把握現(xiàn)在
    短道速滑運(yùn)動(dòng)員非智力因素的培養(yǎng)
    怎樣理解人是戰(zhàn)爭(zhēng)的決定因素?
    軍事歷史(1981年2期)1981-08-14 08:27:58
    一区二区三区免费毛片| 亚洲成av人片免费观看| 免费黄网站久久成人精品 | 国产欧美日韩一区二区三| 国产精品爽爽va在线观看网站| 国产成人影院久久av| 男人狂女人下面高潮的视频| 欧美乱色亚洲激情| 精品人妻1区二区| 国产亚洲欧美在线一区二区| 国产麻豆成人av免费视频| 国内精品久久久久精免费| 女同久久另类99精品国产91| 国产精品,欧美在线| 嫩草影院入口| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品乱码久久久久久99久播| 18禁裸乳无遮挡免费网站照片| 亚洲第一欧美日韩一区二区三区| 国产精品国产高清国产av| 午夜福利18| 午夜福利免费观看在线| 国产精品亚洲av一区麻豆| 国产精品野战在线观看| 97碰自拍视频| 国产精品精品国产色婷婷| 日韩 亚洲 欧美在线| 亚洲在线观看片| 欧美最新免费一区二区三区 | 亚洲成a人片在线一区二区| 精品国内亚洲2022精品成人| 亚洲 国产 在线| 床上黄色一级片| 99在线视频只有这里精品首页| 级片在线观看| 国产精品98久久久久久宅男小说| 亚洲av.av天堂| 日本 欧美在线| 亚洲人成网站在线播| 久久久久国产精品人妻aⅴ院| 国产真实伦视频高清在线观看 | 午夜激情福利司机影院| 91在线观看av| 欧美bdsm另类| 一区二区三区免费毛片| 精品久久久久久久久久免费视频| 亚洲一区二区三区色噜噜| 久久久久久久久久黄片| 一进一出抽搐动态| 欧美+日韩+精品| 国产午夜精品久久久久久一区二区三区 | 午夜精品久久久久久毛片777| 久久久久久久亚洲中文字幕 | 能在线免费观看的黄片| 午夜日韩欧美国产| 亚洲欧美日韩无卡精品| 日本 欧美在线| 亚洲人成网站在线播| 在线十欧美十亚洲十日本专区| 精品人妻视频免费看| 偷拍熟女少妇极品色| x7x7x7水蜜桃| 亚洲最大成人中文| 国产蜜桃级精品一区二区三区| 91在线精品国自产拍蜜月| 日本黄色片子视频| 日韩国内少妇激情av| 色在线成人网| 中亚洲国语对白在线视频| 欧美国产日韩亚洲一区| 精品不卡国产一区二区三区| 色哟哟·www| 日韩有码中文字幕| 午夜日韩欧美国产| 少妇被粗大猛烈的视频| 久久久久九九精品影院| av在线观看视频网站免费| a级毛片a级免费在线| 国产麻豆成人av免费视频| 午夜日韩欧美国产| 久久精品久久久久久噜噜老黄 | 精品日产1卡2卡| 欧美激情久久久久久爽电影| 日韩欧美一区二区三区在线观看| 一级作爱视频免费观看| 精品久久久久久久人妻蜜臀av| 精品人妻偷拍中文字幕| 精品免费久久久久久久清纯| av黄色大香蕉| 国内少妇人妻偷人精品xxx网站| 九色成人免费人妻av| 1024手机看黄色片| 在线观看舔阴道视频| 99国产极品粉嫩在线观看| 国产精品自产拍在线观看55亚洲| 日本黄色片子视频| 99久久九九国产精品国产免费| 男女做爰动态图高潮gif福利片| 国模一区二区三区四区视频| 国产精品电影一区二区三区| 91av网一区二区| 一个人观看的视频www高清免费观看| 人妻久久中文字幕网| 成年版毛片免费区| 亚洲国产欧洲综合997久久,| 日本撒尿小便嘘嘘汇集6| 国产伦在线观看视频一区| 亚洲乱码一区二区免费版| 一级毛片久久久久久久久女| 内地一区二区视频在线| 日本黄色视频三级网站网址| 久久久久久久久大av| 黄片小视频在线播放| 国产探花在线观看一区二区| 精品久久久久久久人妻蜜臀av| 国产欧美日韩精品一区二区| 国产亚洲精品综合一区在线观看| 国产欧美日韩精品一区二区| 三级国产精品欧美在线观看| 美女黄网站色视频| 99国产精品一区二区蜜桃av| 亚洲18禁久久av| 97碰自拍视频| 一a级毛片在线观看| 99国产极品粉嫩在线观看| 99热只有精品国产| 哪里可以看免费的av片| 色哟哟·www| 性色av乱码一区二区三区2| 亚洲成av人片免费观看| 亚洲欧美激情综合另类| 色尼玛亚洲综合影院| 美女被艹到高潮喷水动态| 白带黄色成豆腐渣| 久久热精品热| 亚洲国产精品999在线| 免费高清视频大片| 色精品久久人妻99蜜桃| 欧美中文日本在线观看视频| 日韩欧美国产在线观看| 韩国av一区二区三区四区| 国内精品美女久久久久久| 久久婷婷人人爽人人干人人爱| 亚洲欧美日韩卡通动漫| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 成人午夜高清在线视频| 日韩欧美国产在线观看| 悠悠久久av| 丰满乱子伦码专区| 中文字幕免费在线视频6| 小说图片视频综合网站| 麻豆国产97在线/欧美| 1024手机看黄色片| 欧美绝顶高潮抽搐喷水| 99热只有精品国产| 成人鲁丝片一二三区免费| 欧美最新免费一区二区三区 | 最近在线观看免费完整版| 麻豆一二三区av精品| 无遮挡黄片免费观看| 国产一区二区三区在线臀色熟女| 中文字幕免费在线视频6| 国产黄片美女视频| www.999成人在线观看| 中文字幕精品亚洲无线码一区| 国产精品99久久久久久久久| 欧美黄色片欧美黄色片| 久久精品国产99精品国产亚洲性色| 看十八女毛片水多多多| 久久久久久久午夜电影| 丰满人妻熟妇乱又伦精品不卡| 日韩中文字幕欧美一区二区| 国产精品日韩av在线免费观看| 男人狂女人下面高潮的视频| 欧美日韩黄片免| 久久久精品欧美日韩精品| 极品教师在线视频| 99久久久亚洲精品蜜臀av| 久久久久久久午夜电影| 亚洲久久久久久中文字幕| 动漫黄色视频在线观看| 88av欧美| 蜜桃久久精品国产亚洲av| 身体一侧抽搐| 嫩草影院入口| 亚洲中文日韩欧美视频| 免费av毛片视频| 一进一出抽搐动态| 国产乱人伦免费视频| 免费高清视频大片| 丰满人妻熟妇乱又伦精品不卡| 一区二区三区高清视频在线| 99久久无色码亚洲精品果冻| 小蜜桃在线观看免费完整版高清| 宅男免费午夜| 中文字幕免费在线视频6| 国产白丝娇喘喷水9色精品| 久久精品国产亚洲av香蕉五月| 亚洲最大成人手机在线| 婷婷亚洲欧美| 亚洲激情在线av| x7x7x7水蜜桃| 国产免费一级a男人的天堂| 免费看美女性在线毛片视频| 免费人成在线观看视频色| 最好的美女福利视频网| 麻豆国产av国片精品| 国产精品国产高清国产av| 国产黄a三级三级三级人| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久久久亚洲 | 成人一区二区视频在线观看| 亚洲国产精品sss在线观看| 国产亚洲av嫩草精品影院| 亚洲,欧美,日韩| 久久精品国产清高在天天线| 国产麻豆成人av免费视频| 在线国产一区二区在线| 国产主播在线观看一区二区| 欧美乱妇无乱码| 久久久国产成人精品二区| 欧美成人性av电影在线观看| 日韩亚洲欧美综合| 国产久久久一区二区三区| 好看av亚洲va欧美ⅴa在| 午夜免费成人在线视频| 美女xxoo啪啪120秒动态图 | 嫁个100分男人电影在线观看| 久久久久九九精品影院| av黄色大香蕉| 欧美日韩黄片免| 香蕉av资源在线| 亚洲av熟女| 亚洲美女搞黄在线观看 | a级毛片免费高清观看在线播放| 国产成人aa在线观看| 在线国产一区二区在线| 在线a可以看的网站| 嫩草影视91久久| 最近视频中文字幕2019在线8| 我的老师免费观看完整版| 国产野战对白在线观看| 日韩欧美免费精品| 老司机午夜福利在线观看视频| 在线播放国产精品三级| 小说图片视频综合网站| 婷婷六月久久综合丁香| 国产亚洲av嫩草精品影院| 欧美潮喷喷水| 中文字幕av成人在线电影| 久久久久久大精品| 国产69精品久久久久777片| 一个人看视频在线观看www免费| 少妇丰满av| 一本久久中文字幕| 国产激情偷乱视频一区二区| 亚洲美女黄片视频| 狠狠狠狠99中文字幕| 欧美色欧美亚洲另类二区| 免费观看的影片在线观看| 亚洲三级黄色毛片| 自拍偷自拍亚洲精品老妇| 亚洲五月天丁香| 免费人成在线观看视频色| 日韩成人在线观看一区二区三区| 成人午夜高清在线视频| 色视频www国产| 国内精品久久久久久久电影| 日日摸夜夜添夜夜添av毛片 | 动漫黄色视频在线观看| 国产三级黄色录像| 欧美丝袜亚洲另类 | 国产探花极品一区二区| 男女做爰动态图高潮gif福利片| 国产亚洲精品久久久com| 亚洲国产日韩欧美精品在线观看| 3wmmmm亚洲av在线观看| 男插女下体视频免费在线播放| 91狼人影院| 我的女老师完整版在线观看| 国内精品久久久久久久电影| 国产乱人伦免费视频| 久久精品影院6| 久久久国产成人精品二区| 久久精品国产亚洲av涩爱 | 男女之事视频高清在线观看| 日韩欧美在线二视频| 精品国产亚洲在线| 午夜福利免费观看在线| 中出人妻视频一区二区| 女同久久另类99精品国产91| 长腿黑丝高跟| 老司机午夜十八禁免费视频| 国产精品一区二区三区四区免费观看 | 90打野战视频偷拍视频| 久99久视频精品免费| 伊人久久精品亚洲午夜| 天美传媒精品一区二区| 特级一级黄色大片| 国产精品女同一区二区软件 | 日韩精品青青久久久久久| 给我免费播放毛片高清在线观看| 亚洲 欧美 日韩 在线 免费| 一边摸一边抽搐一进一小说| 麻豆成人av在线观看| 久久6这里有精品| 国产伦在线观看视频一区| 久久久久久久午夜电影| 成人精品一区二区免费| 久久久久久大精品| 免费av观看视频| 亚洲熟妇熟女久久| 国产成人a区在线观看| 国产蜜桃级精品一区二区三区| 婷婷丁香在线五月| 精华霜和精华液先用哪个| 免费在线观看亚洲国产| 日韩欧美国产一区二区入口| 亚洲人成伊人成综合网2020| 日韩欧美国产一区二区入口| 老司机午夜十八禁免费视频| 人人妻,人人澡人人爽秒播| 麻豆成人午夜福利视频| 国产精品1区2区在线观看.| 日日干狠狠操夜夜爽| 亚洲,欧美,日韩| 久久久精品欧美日韩精品| 啦啦啦韩国在线观看视频| 欧美黄色片欧美黄色片| 欧美国产日韩亚洲一区| 日韩免费av在线播放| 国产精品久久视频播放| 国产精品乱码一区二三区的特点| 波多野结衣高清作品| 99久国产av精品| 亚洲欧美日韩无卡精品| 婷婷丁香在线五月| 最新中文字幕久久久久| 色尼玛亚洲综合影院| 一区福利在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久精品久久久久久噜噜老黄 | 黄色日韩在线| 中国美女看黄片| 一个人看的www免费观看视频| 一区二区三区四区激情视频 | 一本一本综合久久| 99久久久亚洲精品蜜臀av| 真人一进一出gif抽搐免费| 亚洲国产高清在线一区二区三| 日韩欧美精品免费久久 | 欧美日韩福利视频一区二区| 日韩欧美 国产精品| 真人一进一出gif抽搐免费| 国产极品精品免费视频能看的| 桃红色精品国产亚洲av| 身体一侧抽搐| 国产美女午夜福利| 午夜免费成人在线视频| 男人狂女人下面高潮的视频| 久久精品久久久久久噜噜老黄 | 在线观看66精品国产| 中亚洲国语对白在线视频| 亚洲av二区三区四区| 欧美成人a在线观看| 99久久九九国产精品国产免费| 欧美成人性av电影在线观看| 亚洲av二区三区四区| 中亚洲国语对白在线视频| 免费av不卡在线播放| 亚洲精品影视一区二区三区av| 成年版毛片免费区| 99riav亚洲国产免费| 久久香蕉精品热| 俺也久久电影网| 亚洲国产高清在线一区二区三| 国产亚洲av嫩草精品影院| 国产极品精品免费视频能看的| 成年人黄色毛片网站| 精品国产亚洲在线| 久99久视频精品免费| 中文字幕av成人在线电影| 国产成人av教育| 久久欧美精品欧美久久欧美| 桃色一区二区三区在线观看| 精华霜和精华液先用哪个| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产自在天天线| 亚洲美女搞黄在线观看 | 午夜福利欧美成人| 亚洲精品乱码久久久v下载方式| 国产精品av视频在线免费观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲激情在线av| 国产高清三级在线| 国产一区二区亚洲精品在线观看| 欧美日本视频| 成人午夜高清在线视频| 国产老妇女一区| 久久久久九九精品影院| 日韩欧美 国产精品| 狂野欧美白嫩少妇大欣赏| 欧美又色又爽又黄视频| 亚洲最大成人手机在线| 看十八女毛片水多多多| 看十八女毛片水多多多| 亚洲av免费在线观看| 最新在线观看一区二区三区| 久久国产精品影院| 最近视频中文字幕2019在线8| 男人狂女人下面高潮的视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久中文看片网| 成年免费大片在线观看| 亚洲第一区二区三区不卡| av视频在线观看入口| 欧美高清成人免费视频www| 国产成年人精品一区二区| 国产久久久一区二区三区| 午夜福利18| 欧美xxxx性猛交bbbb| 动漫黄色视频在线观看| 老鸭窝网址在线观看| 熟妇人妻久久中文字幕3abv| 亚洲熟妇熟女久久| 国产在视频线在精品| av国产免费在线观看| x7x7x7水蜜桃| 国产午夜福利久久久久久| 99久久无色码亚洲精品果冻| 99久久99久久久精品蜜桃| 90打野战视频偷拍视频| 人妻制服诱惑在线中文字幕| 欧美黑人巨大hd| 亚洲av五月六月丁香网| 男人的好看免费观看在线视频| 国产精品亚洲一级av第二区| 国产淫片久久久久久久久 | 国产三级黄色录像| netflix在线观看网站| 一个人免费在线观看电影| 99久久九九国产精品国产免费| 18禁裸乳无遮挡免费网站照片| 国产免费一级a男人的天堂| 精品无人区乱码1区二区| 丁香欧美五月| 日韩欧美国产在线观看| 精品久久国产蜜桃| 特级一级黄色大片| 一边摸一边抽搐一进一小说| 久久热精品热| 99久国产av精品| 久久精品国产清高在天天线| 麻豆一二三区av精品| 日本免费a在线| av福利片在线观看| 好男人在线观看高清免费视频| 亚洲精品亚洲一区二区| 精品日产1卡2卡| 丰满人妻一区二区三区视频av| 日日摸夜夜添夜夜添小说| 18禁在线播放成人免费| 亚洲精品一卡2卡三卡4卡5卡| 在现免费观看毛片| 欧美日本亚洲视频在线播放| 日本在线视频免费播放| 国产一区二区在线av高清观看| 国产精品不卡视频一区二区 | 免费在线观看影片大全网站| 日本黄色片子视频| 日韩欧美一区二区三区在线观看| 国产伦精品一区二区三区四那| 天美传媒精品一区二区| 日韩中文字幕欧美一区二区| 免费在线观看亚洲国产| 男女做爰动态图高潮gif福利片| 午夜精品一区二区三区免费看| 国产精品精品国产色婷婷| 欧美激情久久久久久爽电影| 久久久色成人| 亚洲国产精品sss在线观看| 天堂影院成人在线观看| 淫秽高清视频在线观看| 国语自产精品视频在线第100页| 午夜精品在线福利| 人人妻人人澡欧美一区二区| 国产在视频线在精品| 欧美一区二区国产精品久久精品| 国产精品一区二区三区四区免费观看 | 男女做爰动态图高潮gif福利片| 啦啦啦观看免费观看视频高清| 舔av片在线| 国产精品1区2区在线观看.| 欧美黑人巨大hd| h日本视频在线播放| 俺也久久电影网| 欧美又色又爽又黄视频| 又爽又黄无遮挡网站| 自拍偷自拍亚洲精品老妇| 国产成人欧美在线观看| 88av欧美| 欧美激情国产日韩精品一区| 俺也久久电影网| 99国产极品粉嫩在线观看| 久久伊人香网站| 国产一区二区亚洲精品在线观看| 欧美不卡视频在线免费观看| 午夜精品久久久久久毛片777| 丁香欧美五月| 啦啦啦观看免费观看视频高清| 观看美女的网站| 亚洲国产高清在线一区二区三| xxxwww97欧美| 亚洲av.av天堂| 小蜜桃在线观看免费完整版高清| 男女视频在线观看网站免费| 午夜福利高清视频| 特大巨黑吊av在线直播| 久久人人爽人人爽人人片va | 中文字幕熟女人妻在线| 在线观看66精品国产| 午夜激情欧美在线| 黄色女人牲交| 成人精品一区二区免费| 一区二区三区免费毛片| 精品人妻1区二区| 国产在视频线在精品| 我要看日韩黄色一级片| 中文亚洲av片在线观看爽| 久久久久性生活片| 色吧在线观看| 亚洲,欧美,日韩| 午夜福利视频1000在线观看| 日韩成人在线观看一区二区三区| 亚洲在线观看片| 69人妻影院| 亚洲成人久久性| 亚洲av电影不卡..在线观看| 久久精品影院6| 黄色一级大片看看| 国产麻豆成人av免费视频| 热99在线观看视频| 亚洲中文日韩欧美视频| 午夜老司机福利剧场| 国产三级在线视频| 精品一区二区三区av网在线观看| 亚洲七黄色美女视频| 亚洲成人久久性| 99热这里只有精品一区| 成人特级av手机在线观看| 精品福利观看| 老司机午夜十八禁免费视频| 美女黄网站色视频| 99国产精品一区二区三区| 一本一本综合久久| 少妇丰满av| 小蜜桃在线观看免费完整版高清| 97热精品久久久久久| 亚洲真实伦在线观看| 久久国产精品人妻蜜桃| 日本免费a在线| 99久久久亚洲精品蜜臀av| 内地一区二区视频在线| 亚洲av日韩精品久久久久久密| 一级黄色大片毛片| 首页视频小说图片口味搜索| 亚洲18禁久久av| 国产精品乱码一区二三区的特点| 成人性生交大片免费视频hd| 国内揄拍国产精品人妻在线| 18+在线观看网站| 美女cb高潮喷水在线观看| 丁香欧美五月| 免费一级毛片在线播放高清视频| 色精品久久人妻99蜜桃| 午夜影院日韩av| 国产探花在线观看一区二区| 亚洲精品影视一区二区三区av| 亚洲真实伦在线观看| 天堂网av新在线| 99久久无色码亚洲精品果冻| 日本五十路高清| 久久久久久久亚洲中文字幕 | 日本黄大片高清| 亚洲欧美精品综合久久99| 久久精品综合一区二区三区| 一级黄色大片毛片| 亚洲精品456在线播放app | 精品无人区乱码1区二区| 丰满人妻一区二区三区视频av| 亚洲欧美日韩东京热| 一进一出抽搐动态| 国产精品,欧美在线| 欧美日韩瑟瑟在线播放| 51国产日韩欧美| 丰满人妻一区二区三区视频av| 国内精品一区二区在线观看| 国产中年淑女户外野战色| 亚洲精品色激情综合| 婷婷亚洲欧美| 99热这里只有是精品在线观看 | 内射极品少妇av片p| 亚洲成人久久爱视频| 国产精品免费一区二区三区在线| 757午夜福利合集在线观看| 无人区码免费观看不卡| 88av欧美| 亚洲av二区三区四区| 女人被狂操c到高潮| 成年女人永久免费观看视频| 全区人妻精品视频| 亚洲av不卡在线观看| 噜噜噜噜噜久久久久久91| 麻豆久久精品国产亚洲av|