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

    基于數(shù)據(jù)挖掘與信息融合的制冷設(shè)備故障診斷?

    2021-04-28 16:25:58周旖鋆楊曉燕樓曉華
    振動、測試與診斷 2021年2期
    關(guān)鍵詞:壓縮機故障診斷樣本

    周旖鋆,武 凱,孫 宇,楊曉燕,樓曉華

    (1.南京理工大學機械工程學院 南京,210094) (2.南通四方科技集團股份有限公司 南通,226300)

    引言

    制冷設(shè)備幾種故障同時發(fā)生時,各單發(fā)故障間會相互影響、相互誘導、相互耦合,增加故障診斷系統(tǒng)的誤報率和漏報率。因此,開展制冷設(shè)備并發(fā)故障的診斷具有重要的應用價值。Macalister等[1]基于圖論建立了制冷設(shè)備的故障診斷模型,但這種模型面對較為復雜的系統(tǒng)時可能給出無效診斷結(jié)果。Salvatore等[2]將專家經(jīng)驗的知識庫與推理機制相結(jié)合對制冷設(shè)備故障進行診斷,但這類方法推理過程中存在匹配沖突等問題。Flore等[3]基于主元分析方法在制冷設(shè)備故障檢測中建立主元模型,通過對比測量數(shù)據(jù)與正常數(shù)據(jù)的統(tǒng)計量判斷制冷設(shè)備是否發(fā)生故障,然而這種診斷方法只能執(zhí)行單一故障診斷。Kumar[4]基于指定元分析(designated cell analysis,簡稱DCA)方法對制冷系統(tǒng)故障進行了診斷,但由于DCA方法僅在各特征向量相互正交的情況適用,因此當用于并發(fā)故障診斷時,診斷結(jié)果不夠準確。韓 華 等[5]把 支 持 向 量 機(support vector machines,簡稱SVM)模型用于制冷設(shè)備故障診斷研究并得到了較為準確的診斷結(jié)果。

    為準確檢測出制冷設(shè)備的并發(fā)故障,提出一種基于數(shù)據(jù)挖掘與信息融合的并發(fā)故障診斷方法。首先,采集設(shè)備正常工況數(shù)據(jù)和各類單發(fā)故障與并發(fā)故障下的數(shù)據(jù),并分別建立指定元分析模型和支持向量機模型,為了使指定元分析方法適用于非正交模式,對原來的算法進行改良,提出非完全正交指定元分析方法(non-fully orthogonal designated cell analysis,簡稱NFODCA);其次,采用上述兩種數(shù)據(jù)挖掘方法診斷制冷設(shè)備故障并進行實驗驗證;最后,采用加權(quán)證據(jù)理論對兩種方法的診斷結(jié)果進行信息融合以提高診斷結(jié)果的可信度。

    1 故障數(shù)據(jù)獲取

    制冷設(shè)備可能產(chǎn)生的故障類型較多,根據(jù)經(jīng)驗知識對制冷系統(tǒng)的工作原理進入深入分析后,選取了6類較為典型的單發(fā)故障,通過故障隨機組合形成并發(fā)故障進行研究,所選取的6類典型故障如下:①壓縮機的吸排氣閥片受損,當這類故障發(fā)生時,壓縮機的實際輸氣量下降,制冷效果降低;②制冷劑不足,制冷劑短缺的原因通常有兩個,一種是制冷劑在出廠前充注不足,另一種是由于閥門或焊點松動導致的制冷劑泄漏;③循環(huán)水泵不轉(zhuǎn),這種故障會造成壓縮機吸排氣溫度迅速上升,同時冷凝器的冷卻效果也會急劇下降;④冷卻水的流量太小,冷卻塔的水垢太大或冷卻塔內(nèi)部布置不當可能導致冷卻水流量太?。虎菖蛎涢y開度過小,當這種故障發(fā)生時會造成制冷設(shè)備循環(huán)的制冷劑不足;⑥過濾器堵塞,過濾網(wǎng)用于過濾灰塵和金屬碎屑等,使用時間過長易被堵住。

    實驗對象為一臺需配冷量為55 kW的制冷機組,在原有設(shè)備基礎(chǔ)上引入制冷設(shè)備故障模擬及參數(shù)測試組件,通過實驗采集故障數(shù)據(jù)。各故障的實現(xiàn)方法如下:①壓縮機吸排氣閥片損壞,在壓縮機吸、排氣管路之間設(shè)旁通通路,并安裝針閥或計量閥作為旁通量調(diào)節(jié)裝置;②制冷劑不足,先將系統(tǒng)抽真空,然后通過逐步加注實現(xiàn)不同程度制冷劑不足的故障模擬;③循環(huán)水泵不轉(zhuǎn),斷開水泵控制電路即可模擬循環(huán)水泵不轉(zhuǎn)的故障;④冷卻水流量過小,降低旁通流量或減少冷卻水閘閥開度均可改變冷卻水量;⑤膨脹閥開度過小,將手動調(diào)節(jié)閥置于電磁閥與蒸發(fā)器之間以模擬膨脹閥開度過?。虎捱^濾器堵塞,將一定比例的濾網(wǎng)面積用均布孔洞的厚紙擋住可以模擬過濾器堵塞引起的故障。故障模擬及參數(shù)測試組件安裝位置如圖1所示。

    圖1 故障模擬及參數(shù)測試組件安裝位置Fig.1 Fault simulation and parameter test component installation location

    2 非完全正交指定元分析

    2.1 非完全正交指定元分析模型建立

    由于NFODCA模型在建立時無需并發(fā)故障數(shù)據(jù)作為訓練集,僅在驗證模型階段需要并發(fā)故障的數(shù)據(jù)。因此,將2/3的單發(fā)故障數(shù)據(jù)和2/3正常工況數(shù)據(jù)作為訓練集用于模型建立,將并發(fā)故障數(shù)據(jù)、剩下的1/3類單發(fā)故障數(shù)據(jù)和1/3正常工況數(shù)據(jù)作為測試集用于模型驗證。

    制冷設(shè)備的故障現(xiàn)象、故障原因之間的關(guān)系可以用故障征兆集描述,征兆用測量值與標稱值的偏差度表示。根據(jù)實際運行經(jīng)驗和論證分析,在制冷設(shè)備系統(tǒng)運行過程中對壓縮機吸氣壓力上升、冷凝器進口溫度上升等12種征兆進行監(jiān)測,則每種征兆ur(r=1,2,…,12)可表示為

    將故障種類的數(shù)目記為s,則s種故障Ds可定義為論域U中的6種指定模式

    其中

    選取6種單發(fā)故障:壓縮機吸排氣閥片損壞d1、制冷劑不足d2、循環(huán)水泵不轉(zhuǎn)d3、冷卻水流量過小d4、膨脹閥開度過小d5、過濾器堵塞d6;2種單發(fā)故障組合成的7類并發(fā)故障:d7=d1+d3,d8=d1+d6,d9=d2+d5,d10=d3+d4,d11=d3+d6,d12=d4+d5,d13=d5+d6;3種 單 發(fā) 故 障 組 合 而成的3類并發(fā)故障d14=d1+d2+d4,d15=d3+d4+d6,d16=d4+d5+d6共 計16種 故 障 和 正 常工況進行實驗數(shù)據(jù)采集。以上故障類型根據(jù)經(jīng)驗知識結(jié)合實際工況選出,其他類型的單發(fā)故障及其他組合形成的并發(fā)故障也可采用文中所述的研究方法。

    把故障征兆定義為壓縮機排氣壓力上升u1、下降u2,壓縮機吸氣壓力上升u3、下降u4,壓縮機排氣溫度上升u5、下降u6,壓縮機吸氣溫度上升u7、下降u8,冷凝器出口溫度上升u9、下降u10,過濾器溫度上升u11、下降u12,分析實驗結(jié)果可得故障與征兆間關(guān)系如表1所示。

    表1 制冷裝置故障與征兆間關(guān)系Tab.1 The relationship between faults and symptoms of this quick freezing device

    由表1可知,按照式(2)定義的故障模式并非全部正交,例如:發(fā)生壓縮機吸氣閥片損壞時會導致壓縮機吸氣壓力下降,而發(fā)生循環(huán)水泵不轉(zhuǎn)的故障時會間接導致壓縮機吸氣壓力上升;當這兩種故障同時發(fā)生時,二者的故障特征會相互影響甚至抵消,而傳統(tǒng)指定元分析算法僅在正交模式下適用;這兩種不相互正交的故障同時發(fā)生時,將造成故障漏報。所以對傳統(tǒng)指定元分析方法進行改良,使之適用于非完全正交模式。

    當樣本數(shù)據(jù)Y處于完全正交模式時,Y具有指定元的分解式[6]

    其中:Y由正常子空間和故障子空間組成;di為正交指定模式;wi為相應指定元。

    當Y處于非完全正交模式時可以表示為

    其中:E為殘差空間。

    以EET最小為原則將指定模式集劃分為m個組內(nèi)正交的子集,樣本數(shù)據(jù)Y就可以表示為組內(nèi)正交子集與殘差陣的和

    各子集指定元顯著性可用式(7)計算得出

    在指定模式di所表示方向上用式(8)將樣本數(shù)據(jù)Y做投影,即可得到相應的指定元wi

    其中:n為樣本數(shù)量。

    例如:[y11,y21,…,ym1]T表示第1組樣本數(shù)據(jù)第m個指定模式下的特征向量。計算得到的win表示第i個指定元下第n個樣本所對應的映射值,組成如下的特征向量

    最后,根據(jù)正常工況下的數(shù)據(jù)計算出每個指定元對應特征值的控制上下限并得到Shewhart圖,若圖中特征值超過了控制限,則認為制冷設(shè)備發(fā)生了相應故障。

    2.2 實例驗證

    從測試集中隨機選取實驗數(shù)據(jù),用NFODCA模型進行故障診斷,并將診斷結(jié)果與故障模擬實驗所采集數(shù)據(jù)的故障記錄對比,驗證模型的可靠性。

    根據(jù)表1中的制冷裝置故障與征兆間關(guān)系和EET最小原則,將非完全正交的模式集D={d1,d2,…,d6}表示成3組組內(nèi)正交的模式子集,根據(jù)式(6),樣本數(shù)據(jù)Y可以表示為

    即當非完全正交模式集D劃分為3個組內(nèi)相互正交的模式子集時,Y可由各投影空間完全表示。

    隨機選取一組實驗數(shù)據(jù),分別關(guān)于D1,D2,D3做指定元分析,根據(jù)式(7)計算出該組數(shù)據(jù)對應的各個指定模式的顯著性如表2所示。

    表2 對應各指定模式的顯著性Tab.2 Corresponding to the saliency of each desig?nated mode

    從表2可以看出,該組數(shù)據(jù)對應的6種指定模式中,d1,d4,d6對系統(tǒng)的影響較大,據(jù)此可以初步判斷系統(tǒng)中可能發(fā)生了壓縮機吸排氣閥片損壞、冷卻水流量過小、過濾器堵塞三種故障。為了進一步驗證該方法的合理性,根據(jù)正常工況數(shù)據(jù)在置信度為95%的條件下計算出各個指定模式d1的控制下限L1及上限U1為

    將Y在所有指定元上進行投影,根據(jù)上文計算出的控制限在圖2中給出各指定元的Shewhart圖,圖中U,L分別為控制上、下限,橫坐標為樣本編號,縱坐標表示由式(8),(9)得出的第i個指定元下第n個樣本所對應映射值組成的特征向量的特征值,該量無單位。

    從 圖2可 以 看 出,第1,4,6個 指 定 元 的Shewhart圖在第270個采樣點以后超出了控制限,而第2,3,5個指定元的特征值仍穩(wěn)定在控制限以內(nèi),證明制冷設(shè)備從第270個點開始出現(xiàn)了d1,d4,d63種類型的故障。

    將故障診斷結(jié)果與故障模擬實驗采集數(shù)據(jù)的故障記錄對比,結(jié)果表明,診斷結(jié)果與采集數(shù)據(jù)時的記錄一致。同理可得到其他測試集數(shù)據(jù)經(jīng)NFODCA模型分析后的診斷結(jié)果,并與故障模擬實驗采集數(shù)據(jù)的故障記錄對比,最終得出NFODCA模型的故障診斷準確率為96.94%,虛警率為0.48%,證明文中提出的非完全正交指定元分析方法能準確判斷制冷設(shè)備的并發(fā)故障,是一種有效的多故障診斷方法。

    圖2 各指定元的Shewhart圖Fig.2 Shewhart diagram of each designated cell

    3 支持向量機

    3.1 SVM模型建立

    考慮到SVM模型的診斷效果需要與前面的NFODCA模型對比,仍采用前面的2/3單發(fā)故障數(shù)據(jù)和2/3正常工況數(shù)據(jù)作為訓練集,并發(fā)故障數(shù)據(jù)和剩下的1/3單發(fā)故障、1/3正常工況數(shù)據(jù)作為測試集。通過對訓練集進行訓練,優(yōu)化SVM參數(shù),建立模型,并用測試集數(shù)據(jù)驗證SVM模型的可靠性。

    SVM的主要參數(shù)包含懲罰參數(shù)c以及核函數(shù)參數(shù)g,通常對于這兩個參數(shù)的選擇是經(jīng)過大量的實驗分析得出的[7],耗時耗力。為了快速得到最佳的參數(shù)c和g,使SVM的分類效果更好,采用網(wǎng)格尋優(yōu)和K-折交叉驗證方法。網(wǎng)格尋優(yōu)即先在一定范圍上粗略選擇找出c和g,然后在c,g周圍進行參數(shù)的精細選擇[8]。K-折交叉驗證方法的主要參數(shù)為折疊次數(shù)K,文中取K=5,即5折交叉驗證方法。首先把制冷設(shè)備的數(shù)據(jù)樣本分類成5組子集,將前4組樣本作為訓練集,最后1組樣本用于測試,然后循環(huán)4次,每次都得出一個均方誤差(mean square error,簡稱MSE)。將數(shù)據(jù)進行5次訓練,且每次訓練后將得到的數(shù)據(jù)進行整理分析,然后將5次的MSE取平均值。最后以MSE最小原則選取c值和g值作為SVM建模的最優(yōu)參數(shù)[9]。

    根據(jù)前面得出的最優(yōu)參數(shù),按照分而治之的原則,總故障樣本集Y首先按d1故障分為正樣本(含有d1故障樣本)Fd1和負樣本(不含d1的故障樣本)FNd1,然后對這個二值分類問題建立SVM-d1分類器[10];同理,再將Y分別按故障d2,d3,d4,d5,d6建立二值分類器。如果Y由多個分類器診斷為包含故障,則根據(jù)表決原則,判定故障樣本Y包含多個故障[11]。

    假設(shè)樣本數(shù)據(jù)Y為含有d2,d3故障的數(shù)據(jù),首先用SVM-d1分類器對其進行診斷,Y通過診斷表現(xiàn)為不含d1故障,則其決策函數(shù)R1可以定義為0;然后將Y用SVM-d2分類器診斷,Y通過診斷表現(xiàn)含有d2故障,則其決策函數(shù)R2可以定義為1;同理得到?jīng)Q策函數(shù)R3,R4,R5,R6的值。這樣診斷下來綜合的評定決策函數(shù)R表現(xiàn)為0/1/1/0/0/0,由此可以得到此時的故障Y同時含有d2,d3故障,同理可以得到其他各類故障所對應的決策函數(shù)。

    3.2 SVM模型驗證

    首先找出SVM的主要參數(shù)即懲罰參數(shù)c和核函數(shù)參數(shù)g,SVM模型在數(shù)據(jù)處理過程中自動選擇了參數(shù)c=10,g=0.1,這是隨機得到的結(jié)果,具有主觀性。因此通過網(wǎng)格尋優(yōu)的方法優(yōu)化參數(shù),圖3為參數(shù)c和g的粗略優(yōu)化和精細優(yōu)化過程,橫坐標表示log2c的范圍,縱坐標表示log2g的范圍,均無量綱。

    圖3 參數(shù)粗略和精細優(yōu)化過程Fig.3 Rough and fine optimization of parameters

    圖3 (a)為參數(shù)c和g的粗略優(yōu)化過程:在2-8~28中 粗 略 得 出c=8.031 4,g=4.358 97,MSE=0.109 92,粗略確定了c和g的數(shù)據(jù)范圍;圖3(b)為精細優(yōu)化過程:縮小范圍在2-4~24中查找c和g,最后 得到最 優(yōu)的c和g為cb=6.062 9,gb=2.297 4,并將其作為SVM此次運算的參數(shù)。使用5折交叉驗證方法對并發(fā)故障進行診斷,并將故障診斷加過與故障模擬實驗采集數(shù)據(jù)的故障記錄對比得出故障命中率和虛警率。表3對比了5折交叉驗證下經(jīng)網(wǎng)格優(yōu)化后使用的參數(shù)優(yōu)化前后對并發(fā)故障的診斷結(jié)果。

    由表3可知,選擇默認設(shè)置(c=10,g=0.1)時,制冷設(shè)備故障診斷的平均正確率為46.75%,MSE平均值為0.634;而使用優(yōu)化得到的參數(shù)(c=6.062 9,g=2.297 4)后,制冷設(shè)備故障診斷的正確率有了明顯提高,而且每個折次下都有較高的正確率,平均正確率達到了96.80%,MSE平均值為0.027,證明該SVM模型對制冷設(shè)備故障診斷具有較高的準確率。

    表3 5折交叉驗證下參數(shù)優(yōu)化前后故障診斷結(jié)果Tab.3 Fault diagnosis results before and after pa?rameter optimization under 5-fold cross vali?dation

    4 信息融合

    4.1 SVM模 型 與NFODCA模 型 對 比

    將測試集數(shù)據(jù)分別用支持向量機(SVM)模型和非完全正交指定元分析(NFODCA)模型診斷,各類故障命中率和虛警率對比如表4所示。

    由表4可見,非完全正交指定元分析(NFODCA)和支持向量機(SVM)兩種模型對制冷設(shè)備的故障均有較高的準確率,且兩種方法對不同故障識別度不同,各自在不同類型故障的識別有一定優(yōu)勢。

    4.2 加權(quán)證據(jù)理論信息融合

    對于制冷裝置的并發(fā)故障診斷,如果只采用單一的診斷結(jié)論,極易造成誤診。前面分別從非完全正交指定元分析(NFODCA)和支持向量機(SVM)兩種模型研究了制冷設(shè)備的并發(fā)故障,如果采用多專家集成診斷的方式,將不同方法的診斷結(jié)論進行融合決策,那么將會大大提高制冷設(shè)備故障診斷結(jié)論的準確性和全面性[12]。

    文獻[13]證明基于加權(quán)證據(jù)理論融合具有較好的聚焦能力,因此采用該方法對非完全正交指定元分析(NFODCA)和支持向量機(SVM)兩種模型的診斷結(jié)果進行信息融合。將訓練集中的每組數(shù)據(jù)視為一個證據(jù)組,將16種不同類型的制冷設(shè)備故障定義為一個空間Θ{F1,F(xiàn)2,…,F(xiàn)16}作為辨識框架,Θ是由相互排斥的命題組成的有限完備集。2Θ為Θ的冪集,表示Θ所有子集的集合。在2Θ上定義基本概率分配函數(shù)m:2Θ?[0,1],m滿足下列條件

    表4 SVM和NFODCA模型對各類故障命中率和虛警率對比Tab.4 Comparison of SVM and NFODCA models for various failure hit ratios and false alarm rates %

    筆者采用了兩種故障診斷方法,診斷結(jié)果中有兩種待組合證據(jù)體E1,E2,對應基本概率分配函數(shù)分別為m1,m2,對應焦元分別為Ai,Bj,Ai,Bj的公共焦元為Cl,m1,m2之間的沖突程度系數(shù)為K。根據(jù)不同證據(jù)源的可靠程度分配權(quán)重系數(shù),對相應故障敏感度較高的證據(jù)源賦以較高的權(quán)值,另一個模型賦以較低的權(quán)值,具體權(quán)值由多次實驗后對比平均命中率和虛警率擇優(yōu)選出[14]。設(shè)E1,E2的權(quán)重系數(shù)分別為w1,w2,滿足w1+w2=1,定義N個證據(jù)源對焦元的平均支持程度為

    因此,加權(quán)證據(jù)理論融合公式為

    將NFODCA與SVM模型的證據(jù)體經(jīng)上述方法融合后得到的故障診斷結(jié)果與故障模擬實驗采集數(shù)據(jù)的故障記錄對比,得出故障診斷命中率和虛警率如表5所示。由表5可以看出,制冷設(shè)備的故障識別結(jié)果經(jīng)加權(quán)證據(jù)理論信息融合后更為精確,平均命中率達到了99.10%,平均虛警率降低至0.21%。

    表5 信息融合后各類故障命中率和虛警率Tab.5 Hit rate and false alarm rate of various faults after information fusion %

    5 結(jié)論

    1)非完全指定元分析(NFODCA)和支持向量機(SVM)兩種數(shù)據(jù)挖掘方法在訓練集無并發(fā)故障數(shù)據(jù),僅在測試集中包含并發(fā)故障數(shù)據(jù)的條件下對制冷設(shè)備的故障診斷均有較高的準確率,表明這兩種方法均具有對未知并發(fā)故障類型進行判斷的能力。其中非完全正交指定元分析(NFODCA)的平均命中率達到了96.94%,平均虛警率為0.48%;支持向量機(SVM)的平均命中率達到了96.80%,平均虛警率為0.50%。這兩種方法對不同故障識別度不同,各自在不同類型故障的識別有一定優(yōu)勢。

    2)非完全正交指定元分析(NFODCA)和支持向量機(SVM)兩種方法的故障識別結(jié)果經(jīng)加權(quán)證據(jù)理論信息融合后,平均命中率提高了2.23%,平均虛警率降低了0.28%。說明這兩種數(shù)據(jù)挖掘方法和信息融合技術(shù)結(jié)合應用后可以對制冷設(shè)備的并發(fā)故障得到更為準確的判斷。

    猜你喜歡
    壓縮機故障診斷樣本
    用樣本估計總體復習點撥
    BOG壓縮機在小型LNG船舶上的應用
    淺析特大型往復壓縮機組的安裝
    壓縮機組智能控制與節(jié)能增效
    推動醫(yī)改的“直銷樣本”
    隨機微分方程的樣本Lyapunov二次型估計
    村企共贏的樣本
    因果圖定性分析法及其在故障診斷中的應用
    基于LCD和排列熵的滾動軸承故障診斷
    4M20型壓縮機活塞改造
    一级毛片电影观看| 欧美三级亚洲精品| 男女下面进入的视频免费午夜| 在线观看三级黄色| 一级毛片我不卡| 91精品国产九色| 成人一区二区视频在线观看| 国产免费福利视频在线观看| 丰满人妻一区二区三区视频av| 69人妻影院| 99热全是精品| 80岁老熟妇乱子伦牲交| 日韩av不卡免费在线播放| 国产伦精品一区二区三区视频9| 狂野欧美激情性xxxx在线观看| 色婷婷久久久亚洲欧美| 青春草视频在线免费观看| 色网站视频免费| 看非洲黑人一级黄片| 18禁在线播放成人免费| 亚洲国产最新在线播放| 青春草国产在线视频| 国产91av在线免费观看| 国产乱人偷精品视频| 特大巨黑吊av在线直播| 欧美激情国产日韩精品一区| 亚洲熟女精品中文字幕| 99九九线精品视频在线观看视频| 国产精品一二三区在线看| 日韩视频在线欧美| 大话2 男鬼变身卡| 国产成人精品婷婷| 成年人午夜在线观看视频| videos熟女内射| 黄色日韩在线| 日韩大片免费观看网站| 好男人在线观看高清免费视频| 亚洲精品一区蜜桃| 你懂的网址亚洲精品在线观看| 哪个播放器可以免费观看大片| 丰满乱子伦码专区| 国产老妇伦熟女老妇高清| 欧美高清性xxxxhd video| 在线播放无遮挡| 寂寞人妻少妇视频99o| 亚洲伊人久久精品综合| 亚洲最大成人手机在线| 成年女人在线观看亚洲视频 | 在线观看免费高清a一片| 99热这里只有是精品在线观看| 大又大粗又爽又黄少妇毛片口| 精品国产露脸久久av麻豆| 精华霜和精华液先用哪个| 国产高清国产精品国产三级 | videos熟女内射| 夫妻性生交免费视频一级片| 精品久久久久久久久亚洲| 午夜亚洲福利在线播放| 99热全是精品| 在线精品无人区一区二区三 | 视频中文字幕在线观看| 久久久久久久国产电影| 国产熟女欧美一区二区| 各种免费的搞黄视频| 久久久久九九精品影院| 中文字幕人妻熟人妻熟丝袜美| 国产精品一区二区在线观看99| 男女下面进入的视频免费午夜| 免费观看性生交大片5| 美女视频免费永久观看网站| 国产精品一区二区三区四区免费观看| 国产成人91sexporn| 国产高清有码在线观看视频| 直男gayav资源| 久久久久久久久久人人人人人人| 国产精品人妻久久久影院| 大香蕉97超碰在线| 亚洲av中文字字幕乱码综合| 深夜a级毛片| 偷拍熟女少妇极品色| 一级黄片播放器| 赤兔流量卡办理| 日韩一本色道免费dvd| 国产男女内射视频| 色视频www国产| 久久亚洲国产成人精品v| 国产黄色免费在线视频| 成人漫画全彩无遮挡| 韩国av在线不卡| 国产成人午夜福利电影在线观看| 国产午夜精品久久久久久一区二区三区| 搡老乐熟女国产| 亚洲精品成人av观看孕妇| 赤兔流量卡办理| 男的添女的下面高潮视频| 午夜激情久久久久久久| 三级国产精品片| 国产精品不卡视频一区二区| 国产精品福利在线免费观看| 欧美zozozo另类| 亚洲国产成人一精品久久久| 熟女电影av网| 亚洲内射少妇av| 亚洲欧美中文字幕日韩二区| 精品视频人人做人人爽| 亚洲精品,欧美精品| 久久影院123| 三级国产精品片| 国产精品99久久久久久久久| 草草在线视频免费看| 嫩草影院入口| 成年女人在线观看亚洲视频 | 精品99又大又爽又粗少妇毛片| 成人免费观看视频高清| 日本色播在线视频| 一本久久精品| 真实男女啪啪啪动态图| 2022亚洲国产成人精品| 免费观看性生交大片5| 亚洲美女搞黄在线观看| 三级国产精品欧美在线观看| 高清av免费在线| 好男人在线观看高清免费视频| 久久精品熟女亚洲av麻豆精品| 1000部很黄的大片| 97在线视频观看| 国产美女午夜福利| 中文资源天堂在线| 日本午夜av视频| 色播亚洲综合网| 亚洲美女搞黄在线观看| 白带黄色成豆腐渣| 啦啦啦在线观看免费高清www| 亚洲欧美日韩卡通动漫| 丝袜喷水一区| 亚洲一区二区三区欧美精品 | 三级经典国产精品| 亚洲欧洲日产国产| 日韩 亚洲 欧美在线| 亚洲av中文字字幕乱码综合| 少妇丰满av| 亚洲高清免费不卡视频| 美女被艹到高潮喷水动态| 亚洲av国产av综合av卡| 性插视频无遮挡在线免费观看| 老女人水多毛片| 国产免费一级a男人的天堂| 舔av片在线| 女人被狂操c到高潮| av专区在线播放| 亚洲精品视频女| 久久综合国产亚洲精品| av黄色大香蕉| av在线观看视频网站免费| 久久久精品欧美日韩精品| 国产v大片淫在线免费观看| 在线观看人妻少妇| 久久ye,这里只有精品| 色视频www国产| 亚洲色图av天堂| 18禁裸乳无遮挡动漫免费视频 | 精品熟女少妇av免费看| 国产精品女同一区二区软件| 国产精品麻豆人妻色哟哟久久| 亚洲精品日韩在线中文字幕| 性色avwww在线观看| 秋霞伦理黄片| 久久久久久久久大av| 亚洲四区av| 丝袜喷水一区| 成人黄色视频免费在线看| 禁无遮挡网站| 天堂网av新在线| 亚洲精品,欧美精品| 亚洲av电影在线观看一区二区三区 | 人妻制服诱惑在线中文字幕| 国产一区二区亚洲精品在线观看| 又爽又黄无遮挡网站| 久久久久久久久久人人人人人人| 深爱激情五月婷婷| 日本一本二区三区精品| 美女主播在线视频| 乱系列少妇在线播放| 特级一级黄色大片| 亚洲一级一片aⅴ在线观看| 欧美xxxx性猛交bbbb| 国产精品.久久久| 2018国产大陆天天弄谢| 精品亚洲乱码少妇综合久久| 久久99精品国语久久久| 亚洲四区av| 伦精品一区二区三区| 少妇人妻久久综合中文| 成年免费大片在线观看| 三级国产精品欧美在线观看| 好男人视频免费观看在线| 美女被艹到高潮喷水动态| 3wmmmm亚洲av在线观看| 亚洲精品第二区| 色5月婷婷丁香| 亚洲国产精品专区欧美| 女人被狂操c到高潮| 一级毛片黄色毛片免费观看视频| 色视频在线一区二区三区| 国产大屁股一区二区在线视频| 欧美精品人与动牲交sv欧美| 国模一区二区三区四区视频| 亚洲欧美日韩东京热| 五月玫瑰六月丁香| 精品99又大又爽又粗少妇毛片| 91久久精品电影网| 国产免费一区二区三区四区乱码| 国产伦在线观看视频一区| 国产精品一及| 99热这里只有是精品在线观看| 99视频精品全部免费 在线| 青春草亚洲视频在线观看| 老师上课跳d突然被开到最大视频| 亚洲av欧美aⅴ国产| 国产av国产精品国产| 国产免费又黄又爽又色| 深夜a级毛片| 老司机影院成人| 国产欧美亚洲国产| 青春草亚洲视频在线观看| 老女人水多毛片| 欧美日韩在线观看h| 午夜福利视频1000在线观看| 色5月婷婷丁香| 国产欧美另类精品又又久久亚洲欧美| 欧美极品一区二区三区四区| 国内少妇人妻偷人精品xxx网站| 国精品久久久久久国模美| 精品视频人人做人人爽| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产女主播在线喷水免费视频网站| 天美传媒精品一区二区| 青春草视频在线免费观看| 久久久久精品久久久久真实原创| 女人被狂操c到高潮| 两个人的视频大全免费| 男插女下体视频免费在线播放| 色5月婷婷丁香| 精品午夜福利在线看| 国产精品久久久久久久久免| 日本免费在线观看一区| 国产成人a∨麻豆精品| 免费看不卡的av| 一级毛片 在线播放| 国产一级毛片在线| freevideosex欧美| 精品人妻熟女av久视频| 美女脱内裤让男人舔精品视频| 午夜亚洲福利在线播放| 熟女av电影| 国产精品国产三级国产av玫瑰| 色视频www国产| 最近中文字幕2019免费版| 天美传媒精品一区二区| 国产伦理片在线播放av一区| 日韩强制内射视频| 国产淫语在线视频| 97超视频在线观看视频| av线在线观看网站| 日本色播在线视频| 欧美3d第一页| 欧美精品人与动牲交sv欧美| 一级毛片我不卡| 天天一区二区日本电影三级| a级毛片免费高清观看在线播放| 新久久久久国产一级毛片| 哪个播放器可以免费观看大片| 国产精品一区www在线观看| .国产精品久久| 黄色怎么调成土黄色| 在线观看一区二区三区激情| 亚洲国产精品专区欧美| 狠狠精品人妻久久久久久综合| 男人舔奶头视频| 天美传媒精品一区二区| 精品亚洲乱码少妇综合久久| 亚洲,一卡二卡三卡| 99热6这里只有精品| 日韩欧美一区视频在线观看 | 免费观看性生交大片5| 久久精品国产亚洲av天美| 99re6热这里在线精品视频| 国产爱豆传媒在线观看| 国产精品精品国产色婷婷| 国产乱来视频区| 亚洲国产精品999| 国产精品秋霞免费鲁丝片| 亚洲最大成人手机在线| 美女被艹到高潮喷水动态| 亚洲av日韩在线播放| 亚洲欧美中文字幕日韩二区| 另类亚洲欧美激情| 久久热精品热| 美女xxoo啪啪120秒动态图| 久久久久国产精品人妻一区二区| 国产白丝娇喘喷水9色精品| 成人国产av品久久久| 亚洲av日韩在线播放| 国产精品秋霞免费鲁丝片| 春色校园在线视频观看| 欧美97在线视频| 插逼视频在线观看| 麻豆国产97在线/欧美| 一级黄片播放器| 成人高潮视频无遮挡免费网站| 午夜福利在线在线| 插阴视频在线观看视频| 亚洲av成人精品一区久久| freevideosex欧美| 国产成人免费无遮挡视频| 精品人妻一区二区三区麻豆| 午夜福利在线观看免费完整高清在| 国产亚洲午夜精品一区二区久久 | h日本视频在线播放| www.色视频.com| 男人爽女人下面视频在线观看| 欧美+日韩+精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品中文字幕在线视频 | 超碰av人人做人人爽久久| 日日啪夜夜爽| 国产黄片美女视频| 国产精品爽爽va在线观看网站| 午夜精品国产一区二区电影 | 免费高清在线观看视频在线观看| 一级毛片久久久久久久久女| 精品99又大又爽又粗少妇毛片| 夜夜看夜夜爽夜夜摸| 一级毛片久久久久久久久女| 亚洲av在线观看美女高潮| 亚洲欧美精品自产自拍| 在线观看国产h片| 亚洲精品一区蜜桃| 91午夜精品亚洲一区二区三区| 欧美激情在线99| 国产色婷婷99| 日日啪夜夜撸| 久久97久久精品| 中文在线观看免费www的网站| 亚洲精品日韩在线中文字幕| 亚洲欧美清纯卡通| 亚洲国产日韩一区二区| 色哟哟·www| 国产精品蜜桃在线观看| 欧美xxⅹ黑人| 国精品久久久久久国模美| av专区在线播放| www.色视频.com| 国产免费视频播放在线视频| 久久久久久久久久久丰满| 岛国毛片在线播放| 日韩不卡一区二区三区视频在线| 五月天丁香电影| 伊人久久国产一区二区| 麻豆久久精品国产亚洲av| 国产精品人妻久久久影院| 一级毛片我不卡| 亚洲人成网站高清观看| av.在线天堂| 亚洲怡红院男人天堂| 乱系列少妇在线播放| 国产欧美日韩精品一区二区| 视频区图区小说| 亚洲综合精品二区| 欧美潮喷喷水| 99热国产这里只有精品6| 毛片一级片免费看久久久久| 99久久精品一区二区三区| 国产真实伦视频高清在线观看| 亚洲婷婷狠狠爱综合网| 国产成人a区在线观看| 日韩伦理黄色片| 校园人妻丝袜中文字幕| 中文字幕久久专区| 精华霜和精华液先用哪个| 国产黄色视频一区二区在线观看| 国产探花在线观看一区二区| 久久99热这里只有精品18| 亚洲av免费高清在线观看| 日韩视频在线欧美| 国产精品无大码| 少妇猛男粗大的猛烈进出视频 | 国产精品三级大全| 高清在线视频一区二区三区| 身体一侧抽搐| 亚洲精品影视一区二区三区av| 午夜免费男女啪啪视频观看| 日韩中字成人| 久久久久九九精品影院| 亚洲丝袜综合中文字幕| 国产精品成人在线| 国产综合精华液| 香蕉精品网在线| 你懂的网址亚洲精品在线观看| 欧美一区二区亚洲| 中文字幕久久专区| 国产极品天堂在线| 国产老妇伦熟女老妇高清| 免费av不卡在线播放| 在线观看美女被高潮喷水网站| 午夜精品一区二区三区免费看| 久久久亚洲精品成人影院| 久久6这里有精品| 国产毛片在线视频| 亚洲国产av新网站| 亚洲精品456在线播放app| 成人毛片a级毛片在线播放| 青青草视频在线视频观看| 成人美女网站在线观看视频| 日韩不卡一区二区三区视频在线| 久久久久精品性色| 色5月婷婷丁香| 丰满乱子伦码专区| 精品99又大又爽又粗少妇毛片| 22中文网久久字幕| 久久精品久久精品一区二区三区| 国产亚洲av片在线观看秒播厂| 国产乱人视频| 国产免费一区二区三区四区乱码| 人妻一区二区av| 欧美国产精品一级二级三级 | 啦啦啦在线观看免费高清www| 久久久久久久久久久丰满| 男女啪啪激烈高潮av片| 我的老师免费观看完整版| 青青草视频在线视频观看| 久热久热在线精品观看| 天天躁夜夜躁狠狠久久av| 嫩草影院入口| 亚洲精品一区蜜桃| 天天躁夜夜躁狠狠久久av| 99九九线精品视频在线观看视频| 亚洲不卡免费看| 黄色日韩在线| 狂野欧美激情性xxxx在线观看| 欧美高清性xxxxhd video| 亚洲精品久久久久久婷婷小说| 日本爱情动作片www.在线观看| 观看免费一级毛片| 少妇人妻一区二区三区视频| 老司机影院成人| h日本视频在线播放| 欧美日韩视频高清一区二区三区二| 欧美变态另类bdsm刘玥| 美女被艹到高潮喷水动态| 国产伦精品一区二区三区视频9| 亚洲av不卡在线观看| 国产一区二区三区av在线| 欧美激情在线99| 搞女人的毛片| 精品国产乱码久久久久久小说| 欧美xxxx性猛交bbbb| 国产av码专区亚洲av| 午夜免费鲁丝| 99久久精品一区二区三区| 国产精品人妻久久久影院| 国产精品伦人一区二区| 国产老妇女一区| 免费观看无遮挡的男女| 在线播放无遮挡| 中文精品一卡2卡3卡4更新| 97在线人人人人妻| 少妇人妻 视频| 女人被狂操c到高潮| 免费高清在线观看视频在线观看| 亚洲欧美中文字幕日韩二区| 亚洲第一区二区三区不卡| 日日摸夜夜添夜夜爱| 搞女人的毛片| 听说在线观看完整版免费高清| 精品久久久噜噜| 国产成人免费观看mmmm| 午夜精品国产一区二区电影 | 黄色怎么调成土黄色| 日韩成人伦理影院| 在线观看国产h片| 看免费成人av毛片| 热99国产精品久久久久久7| 男女边摸边吃奶| 男人狂女人下面高潮的视频| 国产精品久久久久久精品电影| 欧美潮喷喷水| 午夜免费鲁丝| 99视频精品全部免费 在线| 日韩 亚洲 欧美在线| 深夜a级毛片| 乱码一卡2卡4卡精品| 少妇的逼好多水| 欧美激情国产日韩精品一区| 美女内射精品一级片tv| 欧美精品一区二区大全| a级毛色黄片| 九九久久精品国产亚洲av麻豆| 国产一区有黄有色的免费视频| 欧美国产精品一级二级三级 | 精品久久久久久久人妻蜜臀av| 超碰av人人做人人爽久久| 成年版毛片免费区| 久久久久精品性色| 91aial.com中文字幕在线观看| 久久精品国产a三级三级三级| 各种免费的搞黄视频| 中国美白少妇内射xxxbb| 国产成人a∨麻豆精品| 亚洲av欧美aⅴ国产| 国产69精品久久久久777片| 99久久精品热视频| 成人一区二区视频在线观看| 午夜免费鲁丝| 亚洲欧美精品自产自拍| 国产乱来视频区| 亚洲国产精品专区欧美| 国产在线男女| 午夜福利视频精品| 亚洲av欧美aⅴ国产| 亚洲av不卡在线观看| 国产欧美另类精品又又久久亚洲欧美| 91狼人影院| 丰满少妇做爰视频| 汤姆久久久久久久影院中文字幕| 99热这里只有精品一区| 国产黄a三级三级三级人| 人妻一区二区av| 国产在视频线精品| 亚洲激情五月婷婷啪啪| 国产91av在线免费观看| 亚洲av福利一区| 日韩,欧美,国产一区二区三区| 插阴视频在线观看视频| 中文欧美无线码| 伊人久久国产一区二区| 听说在线观看完整版免费高清| 成人免费观看视频高清| 在线免费观看不下载黄p国产| 国内揄拍国产精品人妻在线| 美女被艹到高潮喷水动态| 日本av手机在线免费观看| 欧美xxⅹ黑人| 国产69精品久久久久777片| 97超视频在线观看视频| 欧美+日韩+精品| 99久国产av精品国产电影| 亚洲真实伦在线观看| 伊人久久国产一区二区| 亚洲最大成人av| 熟女电影av网| 日本黄色片子视频| 777米奇影视久久| 日本熟妇午夜| 黄色视频在线播放观看不卡| 欧美日韩精品成人综合77777| 国产免费又黄又爽又色| 午夜亚洲福利在线播放| 亚洲精品aⅴ在线观看| 啦啦啦啦在线视频资源| 在线观看三级黄色| 欧美激情国产日韩精品一区| 国产老妇伦熟女老妇高清| 少妇人妻一区二区三区视频| 精品一区二区三卡| 观看美女的网站| 两个人的视频大全免费| 午夜精品国产一区二区电影 | 国产高清三级在线| 午夜免费鲁丝| 成人毛片a级毛片在线播放| 亚洲精品国产av蜜桃| av一本久久久久| 国产高清不卡午夜福利| 麻豆精品久久久久久蜜桃| 少妇人妻久久综合中文| 欧美人与善性xxx| 99久国产av精品国产电影| 精华霜和精华液先用哪个| 午夜免费观看性视频| 日日摸夜夜添夜夜爱| 日韩av在线免费看完整版不卡| 熟妇人妻不卡中文字幕| 久久影院123| 亚洲电影在线观看av| 18禁动态无遮挡网站| 久久久国产一区二区| 成人亚洲精品一区在线观看 | 欧美成人精品欧美一级黄| 在线观看三级黄色| h日本视频在线播放| 一级黄片播放器| 国产成年人精品一区二区| 日韩,欧美,国产一区二区三区| 丝袜脚勾引网站| www.色视频.com| 国产毛片在线视频| 久久久久精品性色| 免费观看a级毛片全部| 久久精品熟女亚洲av麻豆精品| 国产免费一区二区三区四区乱码| 亚洲成人精品中文字幕电影| 亚洲av成人精品一区久久| 少妇猛男粗大的猛烈进出视频 | 亚洲精品乱码久久久v下载方式| 久久久精品94久久精品| 免费观看av网站的网址| 久久久久久久大尺度免费视频| 国国产精品蜜臀av免费| 夫妻性生交免费视频一级片| 国产精品久久久久久精品古装| 啦啦啦在线观看免费高清www| .国产精品久久| 97精品久久久久久久久久精品| 成人国产av品久久久|