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

    基于廣義向量模長的食品安全風險評價體系

    2022-11-22 07:03:46王建新李騰旭李中天王曄茹
    中國食品學報 2022年10期
    關鍵詞:指標值廣義分值

    王建新,李騰旭,李中天,王曄茹

    (1 北京林業(yè)大學信息學院 北京 100083 2 國家林草局智能信息處理中心 北京 100083 3 國家食品安全風險評估中心 北京 100022)

    諸多發(fā)達國家已經(jīng)把食品安全問題提升到國家安全的戰(zhàn)略高度[1]。我國也十分重視食品安全問題,已把食品安全上升至與國防安全同等重要的重大民生工程之一[2]。而食品安全評價是實現(xiàn)食品安全管理的前提和首要條件,它對食品的安全性給出判斷性描述,最終是一個分類和等級評價的問題,這是目前研究的熱點和難點。由于食品安全風險評價是其它食品安全管理活動的前提和基礎,能為食品安全管理部門提供決策輔助[1],也能為社會大眾提供食品安全的知情權,因此是食品安全管理的關鍵環(huán)節(jié)。公平公正、科學合理的評價體系使得后續(xù)的風險管理、預警和響應機制有章可循、有據(jù)可依。

    食品安全風險評價體系的主要議題有4 個,這都是當前研究的熱點。第一是建立指標體系,即評價體系中應該包含什么指標,它們的組織結(jié)構是什么。第二是權重設置,即設置各個指標在評價體系中的數(shù)量化的重要程度。第三是指標值到分值的映射關系。由于一般的評價結(jié)果最終是分值或等級這樣的數(shù)量值,因此,需要把指標中類別化的或數(shù)量化的取值映射為分值或等級值。第四是分項分值到總項分值之間的匯總合成關系。由于指標體系中總項會由若干個分項構成,而每個分項又可能會有若干個更細的分項構成,因此需要建立分項分值到總項之間的匯總合成映射關系。有時這種合成匯總關系比較簡單,包括簡單的加權相加、求均值等,有時按照實際需要會非常復雜。

    要進行食品安全評價,首先需要建立合理的指標體系。指標體系分為單層指標體系和多層指標體系。在單層指標體系中,決定食品安全風險綜合程度或指數(shù)的是若干個變量,這些變量之間的關系是平行的,它們或者是等權重的,或者是非等權重的,所有變量的分值最終決定綜合的食品安全風險程度。為建立乳制品安全風險評估指數(shù),有研究[3]選取了14 個指標,用3 種綜合評價法施加于這14 個指標,結(jié)果具有高度一致性,驗證了指標選取的合理性。呂學東等[4]構建了出口食品生產(chǎn)企業(yè)“三同”成熟度評價體系,也是典型的單層評價指標體系,文中通過百分制給每個指標賦予了一個權重,通過加權得分評價一個企業(yè)的“三同”成熟度。

    然而,在食品安全評價中,單層評價指標體系相對較少,這是因為食品安全評價是一個復雜的系統(tǒng)工程,因此多層指標體系在文獻中更為常見,也更為頻繁地在實踐中使用,常見的層次分析法(AHP 法)也是典型的多層結(jié)構。AHP 法是由美國匹茲堡大學的運籌學教授薩蒂提出的一種定性與定量結(jié)合的多層權重分析方法,該方法把一個復雜問題先看作一個整體,然后逐層分解,一個高層的目標可以分解為若干個低一層級的目標[5]。

    建立指標體系之后,需要為各個指標設置權重,方能通過指標變量的值計算風險得分。設置指標權重通常分為主觀設定和客觀計算兩種方式,有時也會把二者結(jié)合起來操作。有研究[4]由專家直接給出各個指標的權重值。然而,AHP 法中的權重一般由若干個專家分別給出,收集專家給出的權重分值后,需要進行一致性驗證[5]。如果驗證通過,就說明專家們對權重的分配達成基本一致的意見,權重保留;否則,需要專家重新對權重進行分配。除了采用層次分析法[5-6]外,多個應用中用層次分析法建立指標體系和設置權重,包括對農(nóng)副產(chǎn)品供應商的評價體系研究[7],對基于供應鏈的區(qū)域性食品安全狀況的評價[8],對食品生產(chǎn)企業(yè)質(zhì)量安全的評價體系研究[9],以及對我國食品企業(yè)信用評價體系的建設[10]等。AHP 法中的權重設置是典型的主客觀相結(jié)合的方法。熵權法通過指標數(shù)據(jù)確定指標的權重[11-13],是純粹客觀的權重設置法,然而,前提條件是各個待評價的對象在每個指標變量上的取值已經(jīng)全部獲取,因此該方法的適應性較弱,對一些指標值超出已知數(shù)據(jù)范圍的新對象無法評價。主成分分析方法也是一種客觀設置權重的方法,其利用歷史數(shù)據(jù),通過各指標的線性貢獻判定指標權重。例如,研究[3]利用主成分分析方法實現(xiàn)了乳制品安全風險評估指數(shù)中的權重設置。

    在建立指標體系、設置各指標權重后,需要建立指標值到分值的映射關系,映射關系建立的方法有很多種,如等級離散劃分法[2],把評價結(jié)果值區(qū)間[1,5]劃分為5 個安全等級,包括“很差”“很好”等,例如后者的評價值區(qū)間可以是(4.2,5]等。而為了度量哪些變量會影響消費者舉報網(wǎng)絡外賣餐飲店的食品安全違法違規(guī)行為,王浩等[14]自定義了若干個原因變量,這些變量大部分都是5 等級取值的。更多使用的是指標值到得分值之間的線性映射,如采用線性映射[15]計算低碳生產(chǎn)、經(jīng)濟增收、安全供給等方面的綠色發(fā)展指數(shù),最后綜合計算總體的綠色發(fā)展指數(shù)。

    食品安全風險評價體系中另外一個重要的議題就是建立分指標項(分項)分值到總指標項(總項)之間的匯總合成映射關系。不管是單層指標體系[3],還是多層指標體系[4-10],總項分值是n 個分項分值的線性加權和,如公式(1)所示。

    式中,T——總項分值;vi——各分項分值;αi——各分項權重,i=1,2,…,n,且要滿足αi=1。

    1 分項到總項的映射關系

    除了線性映射外,分項到總項分值的映射方式可以有多種,可以更復雜,而不僅僅限于公式(1)這種映射方式。更一般的形式如公式(2)所示。

    式中,f——任意的線性或非線性映射。

    雖然線性映射方式在分項到總項的映射中比較常用,但是它存在一個重要缺陷。在食品安全風險評價中,如果某個分指標的分值較高,即風險較高,則會認為總指標中存在這樣一個“短板”:即使其它分指標的數(shù)量很多,風險都比較小,也不能把總指標判定為低風險。也就是說,在風險評估和評價中,分值較高的分項的影響較大,分指標的分值不適合按照權重線性相加。

    食品安全風險的大小,以及它對監(jiān)管部門和輿情的影響,并不簡單是各個因素嚴重程度的加權和,而是更偏重于突出事件、出現(xiàn)嚴重問題的因素或分指標。實現(xiàn)這一映射關系的傳統(tǒng)做法是使用模糊隸屬度函數(shù)(簡稱為模糊函數(shù))。例如,公式(3)是一個簡單的模糊函數(shù),也稱為最大模糊算子:

    根據(jù)公式(3),總指標分值(代表總體風險)是分指標分值的最大值。也就是說,總指標風險值由分指標風險值中最大的那個決定。這個映射關系比較符合人們對風險的認知:系統(tǒng)風險由系統(tǒng)中的短板決定。有研究應用模糊綜合[16]的方式對高校食堂的食品衛(wèi)生做了評價和分析,易于發(fā)現(xiàn)食堂食品衛(wèi)生的短板問題,以便督促改進,有研究應用更一般的三角模糊法[12]評價冷鏈貨架周期,然而,通過施加三角模糊的運算法則,最終的映射結(jié)果仍呈現(xiàn)為非此即彼的形式。還有研究[11]首先用熵權法確定分指標的風險分值,然后也通過模糊綜合的方式獲得總指標的風險分值。

    另一方面,與風險評估評價相反,若要對安全性進行評價,計算安全指數(shù)等,則分指標的安全值中最小的那個對總指標具有決定作用,這也是一個模糊函數(shù),稱為最小模糊算子:

    公式(4)的意義是:總指標的安全性由最不安全的分指標決定。

    然而,上述模糊評價的方式存在另一個重要的缺陷,即總指標的分值必須是某一個分指標的分值,非此即彼,具有排他性,并不滿足實踐應用中的需求,也就是嚴格單調(diào)性范式要求。所謂嚴格單調(diào)性范式要求是:對于兩個評價對象,若除第k個分指標外,其它分指標的值都相同,則兩個對象的總指標值相對于第k 個分指標的值要滿足嚴格單調(diào)性。嚴格單調(diào)性包括嚴格單調(diào)增和嚴格單調(diào)減。這里以嚴格單調(diào)增為例,形式化地描述嚴格單調(diào)性范式要求。

    假設兩個待評價對象O'和O″在n-1 個分指標上的分值都相同,而在第k 個分指標上的分值不同,也就是二者的分指標值分別是:S1,S2,…Sk-1,Sk',Sk+1…Sn和S1,S2,…Sk-1,Sk″,Sk+1…Sn,且假設Sk'>Sk″,那么嚴格單調(diào)增范式要求需要兩個總指標的分值滿足Tk'>Tk″,然而,最大模糊算子和最小模糊算子顯然都不滿足嚴格單調(diào)范式要求,會出現(xiàn)Tk'=Tk″的情況。例如,假設兩個待評價對象O'和O″只有兩個分指標,它們的分值向量分別是(0.3,0.8)和(0.4,0.8)。根據(jù)最大模糊算子,這兩個待評價對象的總指標分值都是0.8,顯然不滿足嚴格單調(diào)性范式要求。事實上,就風險而言,O″比O'風險更高,二者風險值不應該相等。同理,對最小模糊算子也存在大量不滿足嚴格單調(diào)性范式的例子。

    模糊評價方式雖然突出了重點分指標對總指標分值的決定作用,但是完全忽略了其它分指標對總指標分值的影響,在邏輯上并不合理,在實踐應用中是有缺陷的。

    為了融合和兼顧重點分指標對總指標的決定性影響和嚴格單調(diào)性范式兩個要求,需要新的映射方式,為此,作者提出廣義向量模長的映射方法,并從數(shù)學證明和實踐應用兩個角度說明其優(yōu)勢。

    2 數(shù)據(jù)與方法

    2.1 歐式向量模長映射方法

    有研究[12]采用歐式距離比較兩個商品的貨架期差異,并取得良好的效果。兩個商品的歐式距離由分項計算,如公式(5)所示:

    借鑒這個距離定義,可以建立基于歐式模長的分指標值到總指標值的映射方式,即通過計算分指標值構成的向量與0 向量的距離獲取總指標的分值,如公式(6)所示:

    顯然,歐式向量模長的映射方式是符合嚴格單調(diào)范式的(證明略)。另外,該評價方式也滿足重點分指標對總指標的決定性影響的要求。例如,假設食品安全風險由4 個因素構成,其分值分別是v1,v2,v3,v4,各因素的權重均為0.25,對于某食品,前3 個分指標的分值都是0.1,第4 個分指標的值是0.9,則按照線性加權和映射方式,綜合風險的分值是0.25×0.1+0.25×0.1+0.25×0.1+0.25×0.9=0.3。用歐式向量模長的方式可以做如下計算:

    歐式向量模長的映射方式也有其不足之處。從公式(7)的例子可以看出,突出的分值對總指標的影響程度是固定不變的,因此,歐式模長不具有普適性,不能隨需要而調(diào)整,究其關鍵原因是由于該評價方式?jīng)]有可調(diào)整的超參數(shù)。另外,公式(7)所示歐式模長評價方式?jīng)]有考慮各個分指標的權重,需要在歐式模長的基礎上提出更普適的映射方式。

    2.2 廣義向量模長映射方法

    定義從分指標值到總指標值的廣義向量模長映射方式如公式(8)所示:

    式中,β——超參數(shù),取值范圍可以從負無窮(包括)到正無窮(包括);αi——各個分指標的權重,i=1,2,…,n,且要滿足廣義向量模長的廣義性體現(xiàn)在兩個方面:第一,該模型具有超參數(shù)β,有更廣泛的適應性;第二,與歐式模長映射方式不同,它兼顧了各個分項的權重值。

    首先,證明如下有關漸近性的命題:

    命題1.當超參數(shù)β 趨于正無窮時,廣義向量模長映射退化為最大模糊算子映射。

    證明:為了證明過程的簡潔,先假設共有2 個分指標(有更多分指標時,證明過程類似),也就是n=2,兩個分指標的值分別是v1和v2,且v1>v2。

    從上面的證明結(jié)果可以看出,經(jīng)過廣義向量模長映射,總指標值就是分指標的最大值。這種映射實際上就是最大模糊算子映射。

    命題2.當超參數(shù)β 趨于正0 時,廣義向量模長評價方式退化為幾何加權平均。

    根據(jù)洛必達法則,上式分子分母同時求導:

    因此,

    公式(9)其實是幾何平均的一般形式。特殊地,當α1=α2=0.5 時,

    上式恰好是各個分指標值傳統(tǒng)意義上的幾何平均。下面通過一個例子比較算術平均映射和幾何平均映射。仍然假設有2 個分指標,每個分指標所占權重都是0.5,兩個分指標的值分別是0.3 和0.9,那么它們的算術平均值是0.6,然而,幾何平均值是0.5196,二者雖接近,但又有不同,幾何平均映射的結(jié)果更接近于分指標中較小的值。

    命題3.當超參數(shù)β 趨于負無窮時,廣義向量模長評價方式退化為最小模糊算子。

    于是有,

    因此,

    所得結(jié)果即為分指標中的最小值。這說明,當超參β 趨于負無窮時,廣義向量模長映射即最小模糊算子映射。

    2.3 廣義向量模長映射方法的超參取值區(qū)間

    公式(8)表示廣義向量模長映射方法,其超參β 在正無窮、0、負無窮處均沒有定義,除此之外,在實數(shù)軸上的任何其它地方都有定義,可計算。為此,前文通過求極限的方式,重新在這3 處給出了公式(8)的擴展定義,使得公式(8)的超參β 在包括正、負無窮的整個實數(shù)軸上都有了定義。進一步地,通過前文的證明,可以得出:廣義向量模長映射方法能夠涵蓋前述的幾乎所有的映射方法,包括最大模糊算子法、最小模糊算子法、歐式模長法、算術平均法、幾何平均法,也包括這些方法之間的過渡狀態(tài)方法。廣義向量模長映射法隨著超參β 的變化,與其它映射方法的關系如圖1所示。

    圖1 廣義向量模長映射法與其他映射方法的關系Fig.1 The relationship between generalized vector modulus length mapping method and other methods

    首先,超參β 的值越大,映射結(jié)果就越接近于分指標項中的最大值。極端情況,當β 為正無窮時,映射結(jié)果即分指標項的最大值。相反地,超參β 的值越小,映射結(jié)果就越接近于分指標項中的最小值。極端情況,當β 為負無窮時,映射結(jié)果即分指標項的最小值。特殊地,當β 的值為0,1,2時,廣義向量模長映射法分別退化為幾何平均映射法、算術平均映射法和歐式模長映射法。

    除了超參β 在正負無窮的極端情況,廣義向量模長映射方法均滿足嚴格單調(diào)范式。

    除了上述特殊情況,超參β 在過渡狀態(tài)情況下的映射方法更有意義。當β 位于區(qū)間(1,+∞)中時,映射結(jié)果會更偏向于較大的分項值,并且隨著β 的增加,這種偏向的程度更深。當考慮風險指數(shù)等綜合指標時,可以讓超參β 選擇這個區(qū)間的值。相反地,若考慮安全指數(shù)等指標,則可讓超參取值在區(qū)間(-∞,0)中。超參β 在這個區(qū)間中時,映射獲得的總指標的值更傾向于各個分指標的最小值,并且隨著β 的減小,這種傾向性就更明顯。當超參β 在區(qū)間[0,1]中時,既可以映射風險指數(shù),也可以映射安全指數(shù),然而,由于不能很好地體現(xiàn)關鍵子指標對總指標的決定作用,因此建議盡量不采用。不過,當要把分指標信息盡量完整地反映到總指標中,而又不考慮分指標值的大小對映射方式的影響時,可以考慮用這個區(qū)間中的超參值。與眾多文獻中的映射方式類似,可以讓超參β 取值為1,這時廣義向量模長映射法就退化為線性加權平均法。

    由于廣義向量模長映射法完全由超參數(shù)β 決定,因此,用Gβ表示一個特定的映射方法。例如,G10表示超參為10 的廣義向量模長映射法,G∞就是最大模糊算子映射法。

    3 結(jié)果與分析

    3.1 超參數(shù)對映射關系的影響實例

    根據(jù)廣義向量模長的公式(8),表1列出部分超參數(shù)對總指標分值的影響程度。這里仍然假設只有兩個分指標,其權重α1=α2=0.5,且分項值v1=0.9,v2=0.2。根據(jù)廣義向量模長的公式(8),以及它在一些特殊點上的擴展定義,得到對應關系,見表1。

    從表1可以看出,隨著超參數(shù)β 從正無窮變化到負無窮,總指標項的值從分指標的最大值逐漸變化為分指標的最小值。在實際應用中,可以參考表1中的數(shù)據(jù),按照實際需求,選取合適的映射關系,即選取重要值的分指標對總指標的合適的影響程度。也可以根據(jù)實際需求,選取表中兩個相鄰超參之間的某個中間值。

    表1 廣義向量模長方法中超參數(shù)對總指標分值的影響關系Table 1 Impact of the hyper-parameter to total metric score in the method of generalized vector modulus length

    3.2 廣義向量模長映射法的應用

    為了評估蘭州市柑橘類水果的農(nóng)藥殘留風險,文獻[17]建立了指標體系,其中包括甲拌磷等6 個分指標。本文把廣義向量模長映射法應用于這個評價案例,驗證該方法的有效性。比較結(jié)果如表2所示。

    表2 廣義向量模長映射法與其它方法的應用比較Table 2 Application comparison between the method of generalized vector modulus length and other methods

    觀察表2中線性加權映射法(G1),它的最終結(jié)果與分指標最大風險值0.8 相差甚遠,因此不能反映具有重要值的分指標對總指標的決定作用。相反地,最大模糊映射法(G∞)和廣義向量模長法(G10)的結(jié)果都是或比較接近分指標的最大值,符合實際情況,更有應用價值。

    然而,最小模糊映射法(G-∞)在所有情境中的結(jié)果都是0,因此在考慮風險指數(shù)而不是安全指數(shù)時,不具有應用價值。相反地,在考慮安全指數(shù)計算時,該方法會更為實用。

    情景1 和情景2 的6 個分指標風險值除了“苯醚甲環(huán)唑”外,全部相同,廣義向量模長映射法(G10)捕捉到這個差異,總指標的風險值也有所不同,且滿足嚴格單調(diào)性范式。然而,最大模糊映射法(G∞)不具備捕捉這個差異的能力,其最終的風險值均為0.8,與分指標的最大值完全相同。

    兩個分指標“敵敵畏”和“苯醚甲環(huán)唑”的權重不同,在情景2 和情景3 中交換了風險值,這在線性加權映射法(G1)中得以體現(xiàn),從情景2 的0.32變化為情景3 的0.23,在廣義向量模長映射法(G10)也得以體現(xiàn),從情景2 的0.70947 變化為0.63661,而最大模糊映射法(G∞)不具備捕捉這個變化的能力,它對權重的變化并不敏感。

    綜上分析,廣義向量模長映射法從各個角度看,具有比較大的綜合優(yōu)勢。

    4 結(jié)論

    食品安全風險評價體系的科學合理性和公平公正性對食品安全管理具有重要作用,是其不可或缺的組成部分。好的評價體系能夠使食品安全監(jiān)督檢查有據(jù)可依,可以更好地區(qū)分食品企業(yè)和食品產(chǎn)品,結(jié)合監(jiān)督抽檢和檢驗檢查等起到引導和督促的作用。評價體系建設的第4 個議題(由分到總的映射)中,存在不能體現(xiàn)重要分指標值對總指標值的決定作用、不滿足嚴格單調(diào)性范式等問題。本文提出的廣義向量模長映射法不但能夠很好地體現(xiàn)重要分指標值對總指標值的決定作用,而且能夠滿足嚴格單調(diào)性范式。更重要的是,通過嚴格的數(shù)學推導,證明廣義向量模長映射法是傳統(tǒng)映射方法的統(tǒng)一形式,可通過超參的選取,確定適合應用場景的廣義向量模長映射實例。通過實際的應用場景,驗證了這些結(jié)論。

    在多層的指標體系中,存在多次的從分指標值到總指標值的映射過程,需要有多個超參β。一般來說,同一個指標體系下的超參不要求統(tǒng)一,應根據(jù)實際需求,決定不同映射關系中的不同超參的值。

    雖然提出廣義向量模長的分指標值到總指標值的映射方法,但是指標體系建設中仍然存在諸多問題和難點,后續(xù)將把映射方法與另外3 個議題結(jié)合,形成有機統(tǒng)一的指標體系方法論和建設指南。

    猜你喜歡
    指標值廣義分值
    一起來看看交通違法記分分值有什么變化
    工會博覽(2022年8期)2022-06-30 12:19:30
    Rn中的廣義逆Bonnesen型不等式
    從廣義心腎不交論治慢性心力衰竭
    淺談食品中大腸菌群檢測方法以及指標值的對應關系
    維修性定性要求評價指標融合模型研究
    有限群的廣義交換度
    1995年—2013年地方預算內(nèi)財力、中央返還及上解情況
    宿遷城鎮(zhèn)居民醫(yī)保按病種分值結(jié)算初探
    病種分值結(jié)算模式下的醫(yī)療監(jiān)管之實踐與啟示
    一種新的武器裝備體系能力指標變權算法
    中文字幕熟女人妻在线| 好男人在线观看高清免费视频| 久久天躁狠狠躁夜夜2o2o| 国产精品久久视频播放| 国产亚洲精品久久久com| 亚洲av免费高清在线观看| 99在线视频只有这里精品首页| 毛片一级片免费看久久久久| 国产一区二区三区av在线 | 日本-黄色视频高清免费观看| 在线a可以看的网站| 亚洲中文字幕一区二区三区有码在线看| 1000部很黄的大片| 国产视频一区二区在线看| 亚洲第一区二区三区不卡| 中文亚洲av片在线观看爽| 免费观看在线日韩| 天天躁夜夜躁狠狠久久av| 国产三级中文精品| 日韩中字成人| 亚洲国产精品sss在线观看| 亚洲真实伦在线观看| 欧美日本亚洲视频在线播放| 一a级毛片在线观看| 色综合亚洲欧美另类图片| 村上凉子中文字幕在线| 99久国产av精品| 蜜桃亚洲精品一区二区三区| 国产91av在线免费观看| 日本在线视频免费播放| 国产v大片淫在线免费观看| 91精品国产九色| 中文字幕人妻熟人妻熟丝袜美| 欧美最黄视频在线播放免费| 日韩欧美在线乱码| 国产成人freesex在线 | 精品午夜福利在线看| 国产一区二区在线av高清观看| 欧美一区二区亚洲| 欧美成人a在线观看| 亚洲婷婷狠狠爱综合网| 国产欧美日韩一区二区精品| 一进一出抽搐动态| 在线免费观看不下载黄p国产| 久久草成人影院| 午夜久久久久精精品| 一级毛片aaaaaa免费看小| 伦精品一区二区三区| 狂野欧美白嫩少妇大欣赏| 日韩亚洲欧美综合| 久久精品国产亚洲网站| 免费在线观看成人毛片| 非洲黑人性xxxx精品又粗又长| 综合色丁香网| 一级毛片久久久久久久久女| 黑人高潮一二区| 国产免费男女视频| 97超碰精品成人国产| 免费人成视频x8x8入口观看| 精品免费久久久久久久清纯| 高清毛片免费观看视频网站| 久久精品国产亚洲网站| 亚洲久久久久久中文字幕| 国产高清三级在线| 美女被艹到高潮喷水动态| 免费人成在线观看视频色| 国产亚洲精品av在线| 一夜夜www| 91av网一区二区| 在线免费观看的www视频| 99久久精品热视频| 久久中文看片网| 国产精品人妻久久久影院| 深夜精品福利| 中国国产av一级| av卡一久久| 日本欧美国产在线视频| 热99re8久久精品国产| 欧美成人精品欧美一级黄| 久久精品国产亚洲av天美| 日日摸夜夜添夜夜添小说| 欧美激情久久久久久爽电影| 日本与韩国留学比较| 久久久久久大精品| 亚洲人与动物交配视频| 91麻豆精品激情在线观看国产| 两个人视频免费观看高清| 91精品国产九色| 日韩国内少妇激情av| 亚洲精品一区av在线观看| 九九热线精品视视频播放| 亚洲久久久久久中文字幕| 日本熟妇午夜| 精品人妻视频免费看| 十八禁网站免费在线| 亚洲激情五月婷婷啪啪| 国产av一区在线观看免费| 久久久精品大字幕| 国产成人freesex在线 | 综合色av麻豆| 午夜老司机福利剧场| 一级黄片播放器| 精品久久久久久久久久久久久| 午夜日韩欧美国产| 中文亚洲av片在线观看爽| 国产精品福利在线免费观看| 日韩三级伦理在线观看| 在线观看美女被高潮喷水网站| 18禁裸乳无遮挡免费网站照片| aaaaa片日本免费| 亚洲天堂国产精品一区在线| 国产在线男女| 中文亚洲av片在线观看爽| 无遮挡黄片免费观看| 哪里可以看免费的av片| 毛片一级片免费看久久久久| 欧美中文日本在线观看视频| 亚洲第一电影网av| 国产亚洲精品综合一区在线观看| 国产精品一区二区三区四区免费观看 | 国产精品一区二区三区四区免费观看 | 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久大av| 男女边吃奶边做爰视频| 18禁在线播放成人免费| or卡值多少钱| 老熟妇仑乱视频hdxx| 婷婷精品国产亚洲av在线| 3wmmmm亚洲av在线观看| 欧美性猛交黑人性爽| 国产亚洲av嫩草精品影院| 亚洲欧美精品综合久久99| 国产av麻豆久久久久久久| 亚洲专区国产一区二区| 国产精品福利在线免费观看| 长腿黑丝高跟| 18禁在线无遮挡免费观看视频 | 国产精品亚洲一级av第二区| 此物有八面人人有两片| 国产精品一及| 亚洲av熟女| 免费在线观看成人毛片| 中文字幕久久专区| 黄色日韩在线| 美女大奶头视频| 亚洲精品国产成人久久av| 免费观看的影片在线观看| 亚洲美女搞黄在线观看 | 亚洲av第一区精品v没综合| av女优亚洲男人天堂| 99久久成人亚洲精品观看| 天天一区二区日本电影三级| 久久久久久九九精品二区国产| 少妇被粗大猛烈的视频| 秋霞在线观看毛片| 国产亚洲av嫩草精品影院| 成熟少妇高潮喷水视频| 国产精品不卡视频一区二区| 亚洲一区二区三区色噜噜| 国产蜜桃级精品一区二区三区| 在线观看美女被高潮喷水网站| 如何舔出高潮| 亚洲成人久久性| 免费看美女性在线毛片视频| 日本一本二区三区精品| 22中文网久久字幕| av视频在线观看入口| 国产aⅴ精品一区二区三区波| 国产高清视频在线播放一区| 18禁在线无遮挡免费观看视频 | 成人漫画全彩无遮挡| 亚洲精品乱码久久久v下载方式| 在线观看66精品国产| 精品久久久久久久久av| 别揉我奶头 嗯啊视频| 亚洲一区二区三区色噜噜| 乱人视频在线观看| 日本一本二区三区精品| 精品人妻视频免费看| 一本久久中文字幕| 人妻制服诱惑在线中文字幕| 国国产精品蜜臀av免费| 在线播放国产精品三级| 亚洲中文字幕一区二区三区有码在线看| 丰满人妻一区二区三区视频av| 精品免费久久久久久久清纯| 成人av在线播放网站| 精品人妻一区二区三区麻豆 | 丰满的人妻完整版| 卡戴珊不雅视频在线播放| 晚上一个人看的免费电影| 国产熟女欧美一区二区| or卡值多少钱| 人妻丰满熟妇av一区二区三区| 亚洲欧美日韩无卡精品| 精品一区二区三区视频在线| 久久久久国产精品人妻aⅴ院| 久久人人爽人人爽人人片va| 在线播放国产精品三级| 日韩一本色道免费dvd| 欧美zozozo另类| 99热精品在线国产| 草草在线视频免费看| 亚洲欧美中文字幕日韩二区| 国产成人freesex在线 | h日本视频在线播放| 精品久久国产蜜桃| 尤物成人国产欧美一区二区三区| 日日啪夜夜撸| 精品欧美国产一区二区三| 国产精品99久久久久久久久| 久久人人爽人人爽人人片va| 校园春色视频在线观看| 国产爱豆传媒在线观看| 久久久久国内视频| 亚洲欧美日韩东京热| 久久精品影院6| 国产一区二区三区av在线 | 九色成人免费人妻av| 91久久精品电影网| 精品久久久久久久久久免费视频| 一级黄色大片毛片| 九九久久精品国产亚洲av麻豆| 日韩成人av中文字幕在线观看 | av专区在线播放| 别揉我奶头 嗯啊视频| 黄色一级大片看看| www日本黄色视频网| 亚洲七黄色美女视频| 老司机福利观看| 在线播放无遮挡| 日韩精品青青久久久久久| 亚洲内射少妇av| 国产成人aa在线观看| 精品久久久久久久人妻蜜臀av| 日韩一本色道免费dvd| 嫩草影院精品99| 国产淫片久久久久久久久| 日韩大尺度精品在线看网址| 精品一区二区三区人妻视频| 久久久久久久午夜电影| 久久精品夜色国产| 日韩成人伦理影院| 久久国产乱子免费精品| 欧美日韩综合久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 精品一区二区免费观看| 高清毛片免费观看视频网站| 亚洲一级一片aⅴ在线观看| 亚洲av第一区精品v没综合| 亚洲不卡免费看| 不卡视频在线观看欧美| 2021天堂中文幕一二区在线观| 国产一区二区三区在线臀色熟女| 日韩 亚洲 欧美在线| 久久久欧美国产精品| av在线观看视频网站免费| 亚洲性夜色夜夜综合| 日本黄色片子视频| 搞女人的毛片| 国内揄拍国产精品人妻在线| 性色avwww在线观看| 亚洲高清免费不卡视频| 在线播放国产精品三级| 亚洲中文字幕一区二区三区有码在线看| 床上黄色一级片| 51国产日韩欧美| 国产熟女欧美一区二区| 国产高清视频在线观看网站| 国产黄色视频一区二区在线观看 | 午夜福利在线观看吧| 中文在线观看免费www的网站| 国产精品久久久久久久电影| 亚洲熟妇中文字幕五十中出| 亚洲欧美中文字幕日韩二区| 欧美高清成人免费视频www| 久久久午夜欧美精品| 超碰av人人做人人爽久久| 一区二区三区四区激情视频 | 亚洲人成网站在线播| 亚洲精品一卡2卡三卡4卡5卡| 99久久中文字幕三级久久日本| 日本一二三区视频观看| a级一级毛片免费在线观看| 人人妻人人澡欧美一区二区| 最近中文字幕高清免费大全6| 国内久久婷婷六月综合欲色啪| 亚洲国产高清在线一区二区三| 国产成人福利小说| 国产日本99.免费观看| 国产精品国产高清国产av| 久久精品久久久久久噜噜老黄 | 国语自产精品视频在线第100页| 亚洲一级一片aⅴ在线观看| 国产黄色视频一区二区在线观看 | 十八禁网站免费在线| 一级黄片播放器| 99久久九九国产精品国产免费| 国产精品一及| 久久精品夜色国产| 国产探花极品一区二区| 国产伦精品一区二区三区视频9| 日本一二三区视频观看| 日韩欧美在线乱码| 变态另类丝袜制服| 99热6这里只有精品| 精品少妇黑人巨大在线播放 | 搞女人的毛片| 亚洲综合色惰| 久久久精品94久久精品| 一级av片app| 女同久久另类99精品国产91| 成年女人永久免费观看视频| 国产69精品久久久久777片| 亚洲av五月六月丁香网| 深夜精品福利| 一本久久中文字幕| 国产成人a区在线观看| 男人舔奶头视频| 欧美激情在线99| 99热精品在线国产| 少妇人妻一区二区三区视频| 嫩草影院入口| 亚洲丝袜综合中文字幕| 久久人人精品亚洲av| 免费人成视频x8x8入口观看| 久久人妻av系列| 天天一区二区日本电影三级| 村上凉子中文字幕在线| 日日摸夜夜添夜夜爱| 韩国av在线不卡| 老熟妇乱子伦视频在线观看| 亚洲av成人av| 久久99热6这里只有精品| 男人舔奶头视频| 一个人观看的视频www高清免费观看| 久久久久久久亚洲中文字幕| 久久欧美精品欧美久久欧美| 色播亚洲综合网| 99精品在免费线老司机午夜| 无遮挡黄片免费观看| 俺也久久电影网| 亚洲欧美成人精品一区二区| 最新中文字幕久久久久| 久久久国产成人免费| 欧美日韩国产亚洲二区| 日韩欧美三级三区| 国产日本99.免费观看| 国产高清不卡午夜福利| 久久久久久久午夜电影| 国产片特级美女逼逼视频| 中文字幕av在线有码专区| 久久热精品热| 国产成人福利小说| 在线看三级毛片| 尾随美女入室| 香蕉av资源在线| 国产真实乱freesex| 欧美激情国产日韩精品一区| 99热只有精品国产| 国产亚洲欧美98| 在线天堂最新版资源| 中国国产av一级| 国产乱人偷精品视频| 嫩草影视91久久| 亚洲精品影视一区二区三区av| 免费人成在线观看视频色| 看非洲黑人一级黄片| 最近在线观看免费完整版| 伊人久久精品亚洲午夜| 99久久中文字幕三级久久日本| 欧美中文日本在线观看视频| 永久网站在线| 偷拍熟女少妇极品色| 中文在线观看免费www的网站| av天堂在线播放| 69av精品久久久久久| 99热精品在线国产| 国内揄拍国产精品人妻在线| 国产麻豆成人av免费视频| .国产精品久久| 午夜影院日韩av| 中国美女看黄片| 亚洲,欧美,日韩| 免费高清视频大片| 成人亚洲精品av一区二区| 一进一出抽搐gif免费好疼| 听说在线观看完整版免费高清| 乱码一卡2卡4卡精品| 国产一区二区三区av在线 | 可以在线观看的亚洲视频| 国产白丝娇喘喷水9色精品| 亚洲av美国av| 亚洲精品456在线播放app| 亚洲第一区二区三区不卡| 亚洲av一区综合| 十八禁国产超污无遮挡网站| 老师上课跳d突然被开到最大视频| 中文字幕av成人在线电影| 国模一区二区三区四区视频| 亚洲人成网站在线观看播放| 国产伦精品一区二区三区视频9| 亚洲国产精品成人久久小说 | 亚洲第一区二区三区不卡| 日本爱情动作片www.在线观看 | 国产精品三级大全| 久久久久精品国产欧美久久久| 免费看a级黄色片| 日韩亚洲欧美综合| 亚洲经典国产精华液单| 亚洲最大成人av| 精品久久久久久久久久免费视频| 欧美日韩精品成人综合77777| 在线免费十八禁| 女同久久另类99精品国产91| 国产黄色小视频在线观看| 亚洲三级黄色毛片| 男人舔女人下体高潮全视频| 欧美xxxx性猛交bbbb| 国产熟女欧美一区二区| 嫩草影院新地址| 中文字幕av在线有码专区| 久久久久久久久久黄片| 成人亚洲精品av一区二区| 99国产极品粉嫩在线观看| 久久久成人免费电影| 精品一区二区三区视频在线| 欧美激情久久久久久爽电影| 久久人人爽人人爽人人片va| 国产91av在线免费观看| 此物有八面人人有两片| 亚洲av免费高清在线观看| 网址你懂的国产日韩在线| 18禁在线播放成人免费| 欧美潮喷喷水| 日日撸夜夜添| 亚洲人成网站在线播| 村上凉子中文字幕在线| 亚洲欧美中文字幕日韩二区| 欧美区成人在线视频| 国产老妇女一区| 女生性感内裤真人,穿戴方法视频| av黄色大香蕉| 色在线成人网| 国产黄色视频一区二区在线观看 | 欧美丝袜亚洲另类| 99视频精品全部免费 在线| 成人午夜高清在线视频| 欧美日本视频| 床上黄色一级片| 午夜久久久久精精品| 日韩精品有码人妻一区| 欧美中文日本在线观看视频| 嫩草影视91久久| 简卡轻食公司| 身体一侧抽搐| 美女cb高潮喷水在线观看| 一级黄色大片毛片| 卡戴珊不雅视频在线播放| 搞女人的毛片| 欧美性感艳星| 成人精品一区二区免费| 亚洲精品色激情综合| 色尼玛亚洲综合影院| 在线观看一区二区三区| 性欧美人与动物交配| 婷婷精品国产亚洲av在线| 可以在线观看的亚洲视频| av在线蜜桃| 一本一本综合久久| 欧美日韩一区二区视频在线观看视频在线 | 日本一本二区三区精品| 久久亚洲精品不卡| 女的被弄到高潮叫床怎么办| 无遮挡黄片免费观看| 97超级碰碰碰精品色视频在线观看| 欧美一区二区精品小视频在线| 亚洲四区av| 日本黄色视频三级网站网址| 欧美xxxx黑人xx丫x性爽| 成人高潮视频无遮挡免费网站| 国产成人一区二区在线| 亚洲av免费在线观看| 直男gayav资源| 欧美bdsm另类| 91精品国产九色| 欧美日韩在线观看h| 国产精品一二三区在线看| 91午夜精品亚洲一区二区三区| 成人性生交大片免费视频hd| 色哟哟哟哟哟哟| 日本欧美国产在线视频| 夜夜看夜夜爽夜夜摸| 国产成人影院久久av| 在线观看66精品国产| 亚洲精品影视一区二区三区av| 婷婷亚洲欧美| 亚洲国产精品成人综合色| 亚洲最大成人中文| 亚洲七黄色美女视频| 免费看光身美女| 69人妻影院| 不卡视频在线观看欧美| 日本免费a在线| 伦精品一区二区三区| 在线观看美女被高潮喷水网站| 搡老熟女国产l中国老女人| 最近视频中文字幕2019在线8| 97超碰精品成人国产| 嫩草影院精品99| 欧美日韩国产亚洲二区| 欧美性感艳星| 国产精品99久久久久久久久| 在线观看美女被高潮喷水网站| 成人永久免费在线观看视频| 超碰av人人做人人爽久久| 亚洲无线观看免费| 欧美一区二区亚洲| 亚洲,欧美,日韩| 搡老妇女老女人老熟妇| 一个人看视频在线观看www免费| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 欧美激情在线99| 网址你懂的国产日韩在线| 亚洲精品影视一区二区三区av| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久视频播放| 欧美区成人在线视频| 亚洲精品国产av成人精品 | 蜜桃亚洲精品一区二区三区| 精品午夜福利在线看| 99热这里只有是精品50| 国产男人的电影天堂91| 黄色欧美视频在线观看| 欧美精品国产亚洲| 国产激情偷乱视频一区二区| 久久久久久国产a免费观看| 国产激情偷乱视频一区二区| 亚洲精品国产成人久久av| 亚洲国产日韩欧美精品在线观看| 偷拍熟女少妇极品色| 日本熟妇午夜| av卡一久久| 国产熟女欧美一区二区| 久久午夜亚洲精品久久| 国产免费一级a男人的天堂| 精品人妻熟女av久视频| 神马国产精品三级电影在线观看| 久久99热6这里只有精品| 久久国产乱子免费精品| 精品久久久久久成人av| 久久久久久伊人网av| 亚洲在线自拍视频| 精品久久国产蜜桃| 亚洲欧美成人精品一区二区| 一个人免费在线观看电影| 久久鲁丝午夜福利片| 少妇猛男粗大的猛烈进出视频 | 亚洲av中文字字幕乱码综合| 我要看日韩黄色一级片| 国产精品99久久久久久久久| 欧美在线一区亚洲| 99久久九九国产精品国产免费| 欧美成人一区二区免费高清观看| 亚洲人成网站在线播放欧美日韩| 禁无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄 | 成人国产麻豆网| 国产综合懂色| 亚洲熟妇熟女久久| 久久久久久久久中文| 免费大片18禁| 看十八女毛片水多多多| 日韩欧美精品免费久久| 欧美日韩综合久久久久久| 99久久精品国产国产毛片| 国产女主播在线喷水免费视频网站 | 久久精品夜色国产| 美女免费视频网站| 色尼玛亚洲综合影院| 免费高清视频大片| 尾随美女入室| 男女视频在线观看网站免费| 欧美日韩一区二区视频在线观看视频在线 | 人妻制服诱惑在线中文字幕| 美女内射精品一级片tv| 国产精品一区二区性色av| 精品无人区乱码1区二区| 搡老岳熟女国产| 国内久久婷婷六月综合欲色啪| 欧美精品国产亚洲| 日韩成人伦理影院| 岛国在线免费视频观看| 中文亚洲av片在线观看爽| 人妻少妇偷人精品九色| 精品人妻视频免费看| 亚洲国产精品久久男人天堂| 亚洲熟妇中文字幕五十中出| 亚洲欧美成人精品一区二区| 91久久精品国产一区二区成人| 九九在线视频观看精品| 久99久视频精品免费| 一区福利在线观看| 日日摸夜夜添夜夜添av毛片| 毛片女人毛片| 热99在线观看视频| 国产午夜精品论理片| 国产人妻一区二区三区在| 久久99热这里只有精品18| 国产免费男女视频| 国产黄色视频一区二区在线观看 | 久久精品久久久久久噜噜老黄 | 午夜视频国产福利| 午夜福利视频1000在线观看|