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

    基于ARMA模型的在線油液監(jiān)測故障預(yù)警研究 *

    2019-12-26 07:07:38李美威1謝小鵬1賀石中
    潤滑與密封 2019年12期
    關(guān)鍵詞:界限油液殘差

    李美威1 謝小鵬1 馮 偉 賀石中

    (1.華南理工大學(xué)機(jī)械與汽車工程學(xué)院 廣東廣州 510640; 2.廣州機(jī)械科學(xué)研究院有限公司設(shè)備狀態(tài)檢測研究所 廣東廣州 510700)

    在線油液監(jiān)測是指在設(shè)備不停機(jī)的情況下,通過各種在線傳感器對(duì)油液理化指標(biāo)和鐵磁磨粒等進(jìn)行實(shí)時(shí)監(jiān)測,根據(jù)監(jiān)測參數(shù)來判定設(shè)備的運(yùn)行情況,并進(jìn)行相應(yīng)的故障診斷[1]。早期的異常檢測能夠避免更多的嚴(yán)重后果,保證系統(tǒng)性能和效率并減少維修費(fèi)用。在整個(gè)在線監(jiān)測過程中,當(dāng)系統(tǒng)嚴(yán)重偏離其正常情形時(shí)就會(huì)導(dǎo)致異常發(fā)生[2],異常檢測就是識(shí)別出監(jiān)測數(shù)據(jù)偏離正常情況的程度。對(duì)在線油液監(jiān)測數(shù)據(jù)進(jìn)行趨勢分析和預(yù)測,以提前發(fā)現(xiàn)異常并報(bào)警,可有效避免故障發(fā)生。目前,常用的趨勢分析方法有回歸擬合、統(tǒng)計(jì)分析、支持向量機(jī)、時(shí)間序列分析、灰色理論和人工神經(jīng)網(wǎng)絡(luò)等。趨勢預(yù)測的方法有很多,但是適用于油液監(jiān)測系統(tǒng)的較少[3]。

    針對(duì)油液檢測特征的趨勢分析,國內(nèi)外學(xué)者都展開了大量研究。林麗等人[4]基于油液在線監(jiān)測的磨粒信息對(duì)齒輪箱磨損狀態(tài)進(jìn)行了診斷和預(yù)測。張紅和龔玉[5]以光譜分析為例,運(yùn)用灰色理論對(duì)磨損趨勢進(jìn)行了預(yù)測。高經(jīng)緯等[6]對(duì)內(nèi)燃機(jī)潤滑油中的各元素運(yùn)用時(shí)序模型進(jìn)行了趨勢分析,以判斷內(nèi)燃機(jī)的磨損狀態(tài)。TOBON-MEJIA等[7]使用小波分解技術(shù)和基于高斯混合隱馬爾科夫來估計(jì)軸承的壽命。SOUALHI等[8]運(yùn)用基于向量回歸機(jī)的方法對(duì)軸承進(jìn)行了監(jiān)測分析。上述研究主要是通過離線油液分析或是對(duì)一個(gè)周期內(nèi)某個(gè)特征進(jìn)行實(shí)驗(yàn)擬合分析再運(yùn)用于在線,這往往需要大量的時(shí)間去收集足夠全面的數(shù)據(jù),而且用某個(gè)周期內(nèi)的特征變化代表所有情況是不準(zhǔn)確的。

    在線監(jiān)測數(shù)據(jù)往往是不平衡的,能監(jiān)測到的故障數(shù)據(jù)比正常數(shù)據(jù)通常要少得多,這使得許多傳統(tǒng)的診斷方法并不適用于這種問題[9]。而各種設(shè)備故障的發(fā)生往往會(huì)有一個(gè)潛在過程。針對(duì)在線數(shù)據(jù)的這些特點(diǎn),本文作者對(duì)設(shè)備運(yùn)行平穩(wěn)期的監(jiān)測數(shù)據(jù)構(gòu)建ARMA模型,然后對(duì)整個(gè)時(shí)期的模型殘差進(jìn)行分析,將監(jiān)測數(shù)據(jù)特征量在故障發(fā)生之前的一段時(shí)間內(nèi)的殘差劃分為平穩(wěn)期和故障潛伏期,設(shè)定故障潛伏期殘差界限值,一旦殘差越界即可提前報(bào)警。

    1 ARMA模型

    1.1 模型定義

    自回歸滑動(dòng)平均(Auto Regression Moving Average,ARMA)模型是研究平穩(wěn)隨機(jī)過程的典型方法。ARMA模型認(rèn)為一個(gè)時(shí)間序列的相互依存關(guān)系表現(xiàn)在原始數(shù)據(jù)的延續(xù)性上,在某時(shí)刻的值受到歷史值和噪聲的影響[10-11]。對(duì)于一個(gè)離散的時(shí)間序列{x1,x2,......,xn,......},ARMA數(shù)學(xué)模型表示為

    xt=μ+φ1xt-1+φ2xt-2+......+φpxt-p+εt-

    θ1εt-1-......-θqεt-q

    (1)

    式中:xt是當(dāng)前值;μ是常數(shù)項(xiàng);p和q是模型階數(shù);φi和θi是模型參數(shù);{εt}是白噪聲序列。

    當(dāng)q=0時(shí),ARMA(p,q)模型就退化成自回歸AR(p)模型,即

    xt=μ+φ1xt-1+φ2xt-2+......+φpxt-p+εt

    (2)

    當(dāng)p=0時(shí),ARMA(p,q)模型就退化成移動(dòng)平均MA(q)模型,即

    xt=μ+εt-θ1εt-1-......-θqεt-q

    (3)

    1.2 模型定階和參數(shù)估計(jì)

    在使用ARMA(p,q)模型對(duì)時(shí)間序列進(jìn)行擬合時(shí),首先需要對(duì)時(shí)間序列進(jìn)行平穩(wěn)非白噪聲檢驗(yàn)。可使用單位根檢驗(yàn)和Ljung-Box檢驗(yàn)進(jìn)行序列的平穩(wěn)非白噪聲檢驗(yàn)。ADF檢驗(yàn)通過檢驗(yàn)序列中是否存在單位根來判斷序列的平穩(wěn)性[12],如果序列存在單位根則為非平穩(wěn)序列,否則就是平穩(wěn)序列。Ljung-Box是對(duì)時(shí)間序列是否存在滯后相關(guān)的一種檢驗(yàn)方法[13],如果序列存在滯后相關(guān)則為非白噪聲序列,否則為白噪聲序列。

    檢驗(yàn)完序列為平穩(wěn)非白噪聲后,就可使用ARMA(p,q)模型進(jìn)行分析。首先需要確定模型中的p和q值。ARMA(p,q)模型的階數(shù)主要是根據(jù)序列的自相關(guān)系數(shù)和偏自相關(guān)系數(shù)的拖尾性和截尾性來判斷。其對(duì)應(yīng)的模型如表1所示。

    表1 ARMA模型定階基本原則

    確定模型階數(shù)后,估計(jì)模型中的參數(shù)。使用最小二乘估計(jì)法對(duì)模型參數(shù)進(jìn)行估計(jì)[14],最小二乘法能充實(shí)使用序列值,精度較高。此時(shí),記

    a=(μ,φ1,φ2,......,φp,θ1,θ2,......,θq)T

    (4)

    Ft(a)=μ+φ1xt-1+φ2xt-2+......+φpxt-p-

    θ1εt-1-......-θqεt-q

    (5)

    計(jì)算殘差:

    εt=xt-Ft(a)

    (6)

    計(jì)算殘差平方和:

    φpxt-p+θ1εt-1+......+θqεt-q)2

    (7)

    使式(7)達(dá)到最小的參數(shù)值a即為參數(shù)的估計(jì)值。由于白噪聲序列的值無法確定,所以需要運(yùn)用迭代的方法進(jìn)行計(jì)算。

    1.3 模型檢驗(yàn)與優(yōu)化

    在模型構(gòu)建完成以后,就可以得到擬合殘差序列。如果殘差序列為白噪聲,則認(rèn)為擬合模型充分提取了原始序列中的所有信息,這樣模型才顯著有效。所以需要使用Ljung-Box方法對(duì)殘差序列進(jìn)行白噪聲檢驗(yàn)以判斷模型的有效性。

    在進(jìn)行模型定階時(shí)只是主觀判斷相關(guān)系數(shù)的拖尾性和截尾性來確定p和q的值,但該模型不一定是最優(yōu)模型。選用最小信息量準(zhǔn)則(Akaike Information Criterion,AIC)來進(jìn)行模型優(yōu)化,AIC準(zhǔn)則通過模型的似然函數(shù)值和參數(shù)個(gè)數(shù)來衡量模型的擬合效果[15],即

    (8)

    通過比較多個(gè)模型的AIC值,即可選取AIC最小的模型作為最優(yōu)模型。

    1.4 故障潛伏期界限值設(shè)定

    K均值是一種聚類方法。對(duì)于一維K均值,其算法如下:

    (1)隨機(jī)選取序列中的2個(gè)數(shù)作為初始中心點(diǎn)k1、k2;

    (2)分別計(jì)算所有點(diǎn)到k1和k2的距離d1和d2,若d1≤d2,則記為1類,否則為2類;

    (3)更新中心點(diǎn),將所有1類點(diǎn)的均值賦給k1,所有2類點(diǎn)的均值賦給k2;

    (4)重復(fù)步驟(2)、(3),直至k1和k2不再變化。

    由于數(shù)據(jù)是一維時(shí)間序列,即可使用兩中心點(diǎn)的均值作為界限值,即

    thresholdε=(k1+k2)/2

    (9)

    2 實(shí)驗(yàn)驗(yàn)證與分析

    2.1 實(shí)驗(yàn)數(shù)據(jù)

    以某水電站水輪機(jī)組為研究對(duì)象,用實(shí)時(shí)獲取的推力油槽油液含水量進(jìn)行模型驗(yàn)證和分析。該機(jī)組安裝了廣研檢測的在線油液監(jiān)測儀,可實(shí)時(shí)采集機(jī)組油液的各項(xiàng)特征,采樣周期為1 h。該設(shè)備在運(yùn)行過程中遭受了水污染,油液含水量大幅上升。水分會(huì)使油液乳化,降低油液黏度和油膜厚度,還會(huì)促使油品氧化,加速腐蝕,惡化油質(zhì),所以應(yīng)保證油液中的含水量盡可能地低。文中選取一次含水量超標(biāo)的數(shù)據(jù)作為樣本數(shù)據(jù)進(jìn)行分析,總共454條數(shù)據(jù),根據(jù)油液標(biāo)準(zhǔn)設(shè)定的含水量界限值為1×10-4(質(zhì)量分?jǐn)?shù),下同),含水量實(shí)時(shí)趨勢如圖1所示。

    圖1 原始數(shù)據(jù)時(shí)序Fig 1 Original time series

    2.2 平穩(wěn)序列建模

    使用ADF檢驗(yàn)序列平穩(wěn)性,其檢驗(yàn)結(jié)果如表2所示。

    表2 序列檢驗(yàn)結(jié)果

    表3 序列檢驗(yàn)結(jié)果

    2.2.2 模型定階與參數(shù)估計(jì)

    圖2 序列自相關(guān)函數(shù)Fig 2 Time autocorrelation function (ACF)

    圖3 序列偏相關(guān)函數(shù)Fig 3 Time partial autocorrelation function (ACF)

    從圖2和圖3中可看出,其自相關(guān)表現(xiàn)出明顯的拖尾性,偏相關(guān)在9階之后取值均在置信區(qū)間內(nèi),可認(rèn)為偏相關(guān)9階截尾,根據(jù)表1可判定為AR(9)模型。再分別對(duì)AR(7)、AR(8)、AR(9)、AR(10)、AR(11)進(jìn)行分析,比較它們的AIC,選取最優(yōu)模型。5種模型的AIC值如表4所示。

    表4 各模型的AIC值對(duì)比

    從表4可知,AR(9)模型的AIC值最小,則選取AR(9)模型進(jìn)行序列擬合,并運(yùn)用最小二乘法對(duì)AR(9)模型進(jìn)行參數(shù)估計(jì),其參數(shù)估計(jì)結(jié)果如表5所示

    表5 參數(shù)估計(jì)結(jié)果

    則模型的表達(dá)式為

    xt=26.288 9+0.639 8(xt-1-26.288 9)+

    0.123 6(xt-2-26.288 9)-0.004 1(xt-3-26.288 9)-

    0.020 6(xt-4-26.288 9)+0.003 8(xt-5-26.288 9)+

    0.027 3(xt-6-26.288 9)+0.378 8(xt-7-26.288 9)-

    0.243 5(xt-8-26.288 9)-0.160 6(xt-9-26.288 9)+

    εt=6.716 8+0.639 8xt-1+0.123 6xt-2-0.004 1xt-3-

    0.020 6xt-4+0.003 8xt-5+0.027 3xt-6+0.378 8xt-7-

    0.243 5xt-8-0.160 6xt-9+εt

    圖4 AR(9)模型擬合結(jié)果和原始數(shù)據(jù)比較Fig 4 Comparison of fitting results of AR(9) model and original data

    2.2.3 模型檢驗(yàn)

    驗(yàn)證擬合模型的殘差序列是否為白噪聲序列。同理,使用Ljung-Box檢驗(yàn)方法對(duì)殘差序列的3、6、9期和12期滯后進(jìn)行自相關(guān)分析,其結(jié)果如表6所示。

    表6 殘差Ljung-Box檢驗(yàn)結(jié)果

    從表6中可知?dú)埐钚蛄械?、6、9和12期滯后的檢驗(yàn)概率遠(yuǎn)大于0.05,并且均接近于1,則可認(rèn)為殘差序列為白噪聲序列,模型有效。

    3 殘差序列分析

    圖5 原始時(shí)間序列及其擬合結(jié)果Fig 5 Original time series and its fitting results

    3.1 對(duì)殘差序列進(jìn)行統(tǒng)計(jì)分析

    圖6 不同時(shí)期殘差統(tǒng)計(jì)直方圖Fig 6 Residual statistical histogram of different period

    從圖6中可以看出,平穩(wěn)期殘差是一個(gè)白噪聲序列,大致符合正態(tài)分布。而故障潛伏期的殘差分布發(fā)生了較大變化,對(duì)故障潛伏期的殘差序列進(jìn)一步分析劃分界限值。

    3.2 對(duì)殘差進(jìn)行K均值聚類

    圖7 原始序列報(bào)警示意圖Fig 7 Alarming schematic diagram of original time sequence

    從圖7中可以看出,油液含水量的界限值為1×10-4, 2019年2月14日4∶22時(shí)的實(shí)測值為9.685×10-5,尚未達(dá)到報(bào)警界限值,而且此時(shí)下一步的預(yù)測值為7.83×10-5,也未達(dá)到報(bào)警界限值,即不會(huì)提前報(bào)警。只有在2019年2月14日5∶22時(shí)的實(shí)測值為1.157 2×10-4,超過報(bào)警界限值,系統(tǒng)及時(shí)報(bào)警。即在這種情況下,系統(tǒng)在測量值超過報(bào)警界限值時(shí)會(huì)及時(shí)報(bào)警,但無法提前報(bào)警。

    圖8 殘差報(bào)警示意圖Fig 8 Alarming schematic diagram of residual

    從圖8中可以看出,根據(jù)3sigma原則確定殘差界限值2.643×10-5,在2019年2月14日4∶22時(shí)的殘差值為2.987×10-5,超過殘差界限值,系統(tǒng)發(fā)出警報(bào),即系統(tǒng)提前1 h發(fā)出警報(bào)。而根據(jù)K均值得到的界限值1.271×10-5,在2019年2月14日2∶22時(shí)的殘差值為1.877×10-5,超過殘差界限值,系統(tǒng)發(fā)出警報(bào)。相對(duì)于3sigma原則,使用K均值可在相同的采樣周期下大大增加提前報(bào)警的時(shí)間,而且報(bào)警時(shí)的含水量相對(duì)更小,可有效地防止設(shè)備故障發(fā)生和進(jìn)一步惡化。

    4 結(jié)論

    (1)利用ARMA模型對(duì)設(shè)備平穩(wěn)運(yùn)行時(shí)期的在線油液監(jiān)測數(shù)據(jù)進(jìn)行建模,利用該模型可及時(shí)掌握監(jiān)測特征量的變化趨勢,運(yùn)用殘差來表示設(shè)備偏離平穩(wěn)狀態(tài)的程度。

    (2)針對(duì)ARMA模型對(duì)非平穩(wěn)期的預(yù)測精度較低且只能進(jìn)行一步預(yù)測的特點(diǎn),利用模型擬合殘差判定設(shè)備處于正常期還是故障潛伏期,并通過K均值聚類將設(shè)備發(fā)生故障前的殘差進(jìn)行二分類,將兩類中心點(diǎn)的均值作為正常期與故障潛伏期間的界限值。

    (3)通過實(shí)際的監(jiān)測數(shù)據(jù)進(jìn)行模型驗(yàn)證和分析,結(jié)果表明,通過此方法可提前系統(tǒng)警報(bào)時(shí)間,有效控制油液監(jiān)測特征量的值并且防止設(shè)備故障發(fā)生。

    猜你喜歡
    界限油液殘差
    基于在線監(jiān)測的油液污染智能控制系統(tǒng)
    界限
    十幾歲(2022年21期)2022-11-19 11:14:42
    基于雙向GRU與殘差擬合的車輛跟馳建模
    間隙
    基于殘差學(xué)習(xí)的自適應(yīng)無人機(jī)目標(biāo)跟蹤算法
    通過油液分析查找發(fā)動(dòng)機(jī)抱軸故障原因
    破次元
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    承諾是跨越時(shí)間界限的恒久
    中國寶玉石(2016年6期)2017-01-03 09:37:07
    合理潤滑與油液監(jiān)測技術(shù)的應(yīng)用
    中文字幕最新亚洲高清| 三上悠亚av全集在线观看| 中文字幕制服av| 女人被躁到高潮嗷嗷叫费观| 一级黄片播放器| 亚洲精品第二区| 又黄又粗又硬又大视频| 精品卡一卡二卡四卡免费| 久久午夜综合久久蜜桃| 青春草亚洲视频在线观看| 成人国语在线视频| 亚洲欧美成人综合另类久久久| 久久人妻熟女aⅴ| 成人午夜精彩视频在线观看| 建设人人有责人人尽责人人享有的| av电影中文网址| 亚洲综合色惰| 欧美日韩精品成人综合77777| 一区二区三区精品91| 国产精品秋霞免费鲁丝片| 黑人猛操日本美女一级片| 制服人妻中文乱码| 2018国产大陆天天弄谢| av福利片在线| 最新中文字幕久久久久| 91午夜精品亚洲一区二区三区| 国产一区二区激情短视频 | 亚洲精品在线美女| 99热全是精品| videossex国产| 精品少妇久久久久久888优播| 777米奇影视久久| av福利片在线| 日韩一本色道免费dvd| 免费观看在线日韩| 欧美最新免费一区二区三区| 亚洲精品美女久久久久99蜜臀 | 亚洲av男天堂| 国产一区二区三区av在线| 人体艺术视频欧美日本| a级片在线免费高清观看视频| 午夜福利一区二区在线看| 国产精品麻豆人妻色哟哟久久| 下体分泌物呈黄色| 国产免费一区二区三区四区乱码| av天堂久久9| 免费久久久久久久精品成人欧美视频| 精品一区二区免费观看| 免费在线观看黄色视频的| 1024视频免费在线观看| 一级毛片黄色毛片免费观看视频| 在线精品无人区一区二区三| 久热这里只有精品99| 青春草视频在线免费观看| 久久久a久久爽久久v久久| 国产精品国产三级国产专区5o| 久久99热这里只频精品6学生| av又黄又爽大尺度在线免费看| 日本黄色日本黄色录像| 另类精品久久| av不卡在线播放| av一本久久久久| 十分钟在线观看高清视频www| 欧美人与性动交α欧美精品济南到 | 9191精品国产免费久久| 精品少妇久久久久久888优播| www.熟女人妻精品国产| 久久午夜综合久久蜜桃| 亚洲精品国产一区二区精华液| 可以免费在线观看a视频的电影网站 | 国产亚洲av片在线观看秒播厂| 999久久久国产精品视频| 在线观看www视频免费| 精品国产超薄肉色丝袜足j| 午夜福利影视在线免费观看| 啦啦啦中文免费视频观看日本| 男女国产视频网站| 男女国产视频网站| 激情五月婷婷亚洲| 成人影院久久| 亚洲欧美一区二区三区久久| 一二三四在线观看免费中文在| 精品国产超薄肉色丝袜足j| 午夜福利视频在线观看免费| 欧美人与善性xxx| 中文天堂在线官网| 一级毛片 在线播放| 日韩免费高清中文字幕av| 天天躁狠狠躁夜夜躁狠狠躁| 边亲边吃奶的免费视频| 久久久久久久久免费视频了| 久久精品国产亚洲av天美| 叶爱在线成人免费视频播放| 97在线视频观看| 韩国av在线不卡| 人人妻人人爽人人添夜夜欢视频| 99热网站在线观看| 欧美日韩精品网址| 欧美日本中文国产一区发布| 日韩一卡2卡3卡4卡2021年| 一级毛片我不卡| 熟妇人妻不卡中文字幕| 男人添女人高潮全过程视频| 久久精品aⅴ一区二区三区四区 | 99国产精品免费福利视频| 久久久久久久久久久免费av| 亚洲成人av在线免费| 性高湖久久久久久久久免费观看| 成人亚洲精品一区在线观看| 一级片免费观看大全| 国产精品一国产av| 中文精品一卡2卡3卡4更新| 99久久中文字幕三级久久日本| 啦啦啦在线观看免费高清www| 久久综合国产亚洲精品| 丰满乱子伦码专区| 久久久久久久久免费视频了| 国产免费视频播放在线视频| 亚洲情色 制服丝袜| 激情五月婷婷亚洲| av国产精品久久久久影院| 国产av国产精品国产| 久久韩国三级中文字幕| 老司机亚洲免费影院| 一本久久精品| 中文字幕精品免费在线观看视频| 久久午夜福利片| 春色校园在线视频观看| 一区二区三区四区激情视频| 久久99热这里只频精品6学生| 91久久精品国产一区二区三区| 免费久久久久久久精品成人欧美视频| 欧美精品国产亚洲| 亚洲综合色网址| 精品午夜福利在线看| 极品人妻少妇av视频| 国产成人91sexporn| 人人妻人人澡人人看| 日本vs欧美在线观看视频| 美女国产视频在线观看| 亚洲av.av天堂| 搡老乐熟女国产| 制服人妻中文乱码| 亚洲精品久久午夜乱码| 久久99热这里只频精品6学生| 校园人妻丝袜中文字幕| 日日爽夜夜爽网站| 亚洲国产欧美网| 嫩草影院入口| 亚洲成人手机| 欧美日韩国产mv在线观看视频| 国产成人a∨麻豆精品| 卡戴珊不雅视频在线播放| 啦啦啦中文免费视频观看日本| 大片免费播放器 马上看| 9色porny在线观看| 热re99久久国产66热| 毛片一级片免费看久久久久| 91精品国产国语对白视频| 伊人亚洲综合成人网| 亚洲精品美女久久av网站| 久久精品国产a三级三级三级| 麻豆精品久久久久久蜜桃| 丝袜美足系列| 午夜日韩欧美国产| 国语对白做爰xxxⅹ性视频网站| 老司机影院成人| 麻豆精品久久久久久蜜桃| 日韩制服丝袜自拍偷拍| 青春草视频在线免费观看| 亚洲,欧美精品.| 日韩熟女老妇一区二区性免费视频| 久久久久久久久久人人人人人人| a级毛片在线看网站| 欧美亚洲 丝袜 人妻 在线| 久久99精品国语久久久| 少妇人妻 视频| 午夜久久久在线观看| 在线观看美女被高潮喷水网站| 精品卡一卡二卡四卡免费| 日韩 亚洲 欧美在线| 中文字幕制服av| 免费少妇av软件| 亚洲中文av在线| 永久免费av网站大全| 制服诱惑二区| 一本—道久久a久久精品蜜桃钙片| 中文字幕色久视频| 搡老乐熟女国产| 校园人妻丝袜中文字幕| 女性被躁到高潮视频| 超碰97精品在线观看| videos熟女内射| 十八禁网站网址无遮挡| 91午夜精品亚洲一区二区三区| 观看av在线不卡| 又粗又硬又长又爽又黄的视频| 91久久精品国产一区二区三区| 亚洲av电影在线观看一区二区三区| 午夜老司机福利剧场| 99久国产av精品国产电影| 亚洲精品第二区| 亚洲美女搞黄在线观看| 国产亚洲av片在线观看秒播厂| 精品一品国产午夜福利视频| 黄色配什么色好看| 在线天堂最新版资源| 久久久久久久久久久久大奶| 青草久久国产| 亚洲av男天堂| 久久久精品94久久精品| 日本爱情动作片www.在线观看| 亚洲精品av麻豆狂野| 欧美日韩国产mv在线观看视频| 最近2019中文字幕mv第一页| 亚洲美女搞黄在线观看| av在线老鸭窝| 少妇精品久久久久久久| 热99久久久久精品小说推荐| 波野结衣二区三区在线| 交换朋友夫妻互换小说| 美女国产视频在线观看| 日本猛色少妇xxxxx猛交久久| 99九九在线精品视频| 美女国产高潮福利片在线看| 久热久热在线精品观看| 久久精品国产亚洲av天美| 国产xxxxx性猛交| 免费在线观看完整版高清| 你懂的网址亚洲精品在线观看| 日韩不卡一区二区三区视频在线| 大片电影免费在线观看免费| 午夜免费鲁丝| 亚洲欧美精品综合一区二区三区 | 欧美在线黄色| 90打野战视频偷拍视频| 亚洲av欧美aⅴ国产| 国产 一区精品| 熟女少妇亚洲综合色aaa.| 嫩草影院入口| 久久精品久久久久久久性| 老熟女久久久| 久久久精品免费免费高清| 亚洲中文av在线| 视频在线观看一区二区三区| 一级毛片 在线播放| 日日摸夜夜添夜夜爱| 一级a爱视频在线免费观看| 秋霞伦理黄片| 交换朋友夫妻互换小说| 国产精品 欧美亚洲| 秋霞在线观看毛片| 亚洲男人天堂网一区| 哪个播放器可以免费观看大片| 高清不卡的av网站| 十八禁高潮呻吟视频| 久久久久久久久久久免费av| 伊人久久大香线蕉亚洲五| 国产毛片在线视频| 免费人妻精品一区二区三区视频| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品成人久久小说| 久热这里只有精品99| av在线观看视频网站免费| 欧美另类一区| 国产精品久久久久久精品电影小说| 欧美激情 高清一区二区三区| 免费人妻精品一区二区三区视频| 久久人人爽av亚洲精品天堂| 久久久久久免费高清国产稀缺| 人人妻人人澡人人看| 国产成人a∨麻豆精品| 涩涩av久久男人的天堂| 在线看a的网站| 亚洲国产精品一区三区| 成年av动漫网址| 香蕉丝袜av| 日日摸夜夜添夜夜爱| a级片在线免费高清观看视频| 亚洲成国产人片在线观看| 亚洲国产成人一精品久久久| 一区二区三区精品91| 人妻少妇偷人精品九色| 久久精品国产亚洲av高清一级| 国产成人a∨麻豆精品| 欧美+日韩+精品| 免费av中文字幕在线| 在线观看人妻少妇| 美女国产视频在线观看| 亚洲av在线观看美女高潮| 香蕉国产在线看| 狂野欧美激情性bbbbbb| 日韩一本色道免费dvd| 不卡av一区二区三区| 夫妻性生交免费视频一级片| 丝袜美足系列| 亚洲av电影在线进入| 免费观看性生交大片5| 国产成人午夜福利电影在线观看| 啦啦啦在线观看免费高清www| 伊人久久国产一区二区| 一区在线观看完整版| 麻豆精品久久久久久蜜桃| 亚洲av欧美aⅴ国产| 免费在线观看视频国产中文字幕亚洲 | 肉色欧美久久久久久久蜜桃| 大片免费播放器 马上看| 久久久久精品人妻al黑| 亚洲国产精品一区三区| 久热这里只有精品99| 欧美精品国产亚洲| 青春草视频在线免费观看| 18禁裸乳无遮挡动漫免费视频| www.自偷自拍.com| 久久这里只有精品19| 老汉色∧v一级毛片| 国产黄色免费在线视频| 黄色视频在线播放观看不卡| 伊人亚洲综合成人网| 在现免费观看毛片| 久久国产精品大桥未久av| 一二三四中文在线观看免费高清| 亚洲中文av在线| 日本爱情动作片www.在线观看| 亚洲第一青青草原| 青青草视频在线视频观看| 1024香蕉在线观看| 黑人欧美特级aaaaaa片| 91精品三级在线观看| 搡老乐熟女国产| 赤兔流量卡办理| 国产精品免费大片| 欧美97在线视频| 久久久久视频综合| 男人添女人高潮全过程视频| 国产精品免费大片| 少妇被粗大的猛进出69影院| 搡老乐熟女国产| 国产精品国产av在线观看| 色哟哟·www| 亚洲国产欧美日韩在线播放| 永久网站在线| 青春草亚洲视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 天堂中文最新版在线下载| 成人漫画全彩无遮挡| av在线老鸭窝| 超色免费av| 嫩草影院入口| kizo精华| 亚洲综合色惰| 建设人人有责人人尽责人人享有的| 美女高潮到喷水免费观看| 亚洲精品乱久久久久久| 久久精品久久久久久久性| 国产av精品麻豆| 精品第一国产精品| 波多野结衣av一区二区av| 亚洲美女黄色视频免费看| 国语对白做爰xxxⅹ性视频网站| 在线精品无人区一区二区三| 夫妻午夜视频| 十八禁高潮呻吟视频| 91午夜精品亚洲一区二区三区| 91精品伊人久久大香线蕉| 午夜福利一区二区在线看| 国产精品秋霞免费鲁丝片| 侵犯人妻中文字幕一二三四区| videosex国产| 国产精品 欧美亚洲| 亚洲伊人色综图| 亚洲人成电影观看| 丰满少妇做爰视频| 亚洲男人天堂网一区| 各种免费的搞黄视频| 性色av一级| 永久网站在线| 国产1区2区3区精品| 久久久精品免费免费高清| 极品人妻少妇av视频| 人妻系列 视频| 黄频高清免费视频| 亚洲精品久久午夜乱码| 国产精品香港三级国产av潘金莲 | 午夜福利一区二区在线看| 桃花免费在线播放| 精品少妇一区二区三区视频日本电影 | 色婷婷久久久亚洲欧美| 国产一区二区三区综合在线观看| 中文精品一卡2卡3卡4更新| 狠狠精品人妻久久久久久综合| 国产精品 国内视频| 一本久久精品| 国产一区二区三区综合在线观看| 2018国产大陆天天弄谢| av线在线观看网站| 国产精品成人在线| a级毛片黄视频| 夫妻性生交免费视频一级片| 久久热在线av| 成年动漫av网址| av片东京热男人的天堂| 天堂俺去俺来也www色官网| 在线天堂最新版资源| 久久青草综合色| 午夜免费观看性视频| 国产成人91sexporn| 亚洲精品中文字幕在线视频| 边亲边吃奶的免费视频| 亚洲国产最新在线播放| 精品少妇黑人巨大在线播放| 国产精品三级大全| 老女人水多毛片| 哪个播放器可以免费观看大片| 久久精品夜色国产| 国产高清不卡午夜福利| 夫妻午夜视频| 国产野战对白在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲国产色片| 精品国产露脸久久av麻豆| 国产综合精华液| 亚洲国产看品久久| 少妇人妻久久综合中文| 免费看av在线观看网站| 日韩人妻精品一区2区三区| 亚洲欧美清纯卡通| 黄色配什么色好看| 青春草亚洲视频在线观看| 美女中出高潮动态图| 色婷婷av一区二区三区视频| 黄色 视频免费看| 亚洲激情五月婷婷啪啪| 99久国产av精品国产电影| 亚洲av电影在线进入| 国产精品熟女久久久久浪| 精品人妻熟女毛片av久久网站| 丰满乱子伦码专区| 日日啪夜夜爽| 久久久精品国产亚洲av高清涩受| 建设人人有责人人尽责人人享有的| 寂寞人妻少妇视频99o| av卡一久久| 大香蕉久久网| 丰满迷人的少妇在线观看| 男女无遮挡免费网站观看| av片东京热男人的天堂| 久久久久久久国产电影| 两性夫妻黄色片| 亚洲精品成人av观看孕妇| 亚洲精品aⅴ在线观看| 又黄又粗又硬又大视频| 亚洲av电影在线进入| 最近中文字幕高清免费大全6| 亚洲精华国产精华液的使用体验| 欧美精品一区二区免费开放| 人妻人人澡人人爽人人| 久久久久久久久久久免费av| 久久久久久久精品精品| 宅男免费午夜| 国产片内射在线| 成年动漫av网址| 国产精品人妻久久久影院| www.自偷自拍.com| 五月伊人婷婷丁香| 性色av一级| 欧美日韩成人在线一区二区| 亚洲欧美精品综合一区二区三区 | 中国三级夫妇交换| av线在线观看网站| 久久精品久久久久久噜噜老黄| 亚洲av综合色区一区| 最新的欧美精品一区二区| 9热在线视频观看99| 久久精品国产亚洲av高清一级| 午夜老司机福利剧场| 久久女婷五月综合色啪小说| 一区二区三区四区激情视频| 99re6热这里在线精品视频| 成人18禁高潮啪啪吃奶动态图| av视频免费观看在线观看| 国产亚洲欧美精品永久| 在线观看一区二区三区激情| 久久久久久伊人网av| 999久久久国产精品视频| 精品国产乱码久久久久久小说| 中文字幕制服av| 国产 一区精品| 妹子高潮喷水视频| 一区二区三区精品91| 亚洲精品成人av观看孕妇| a级毛片在线看网站| 人妻少妇偷人精品九色| 交换朋友夫妻互换小说| a级毛片黄视频| 高清黄色对白视频在线免费看| 最近手机中文字幕大全| 国产精品av久久久久免费| 嫩草影院入口| 午夜福利在线免费观看网站| 少妇猛男粗大的猛烈进出视频| 黄色一级大片看看| 日本爱情动作片www.在线观看| 精品少妇黑人巨大在线播放| 天堂8中文在线网| 精品少妇黑人巨大在线播放| 国产成人av激情在线播放| 欧美激情 高清一区二区三区| 久久久久人妻精品一区果冻| 国产成人午夜福利电影在线观看| 国产视频首页在线观看| 亚洲精品一二三| 丰满迷人的少妇在线观看| 欧美激情 高清一区二区三区| 成人漫画全彩无遮挡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一级a爱视频在线免费观看| 一级,二级,三级黄色视频| 亚洲五月色婷婷综合| 这个男人来自地球电影免费观看 | 国产爽快片一区二区三区| 成人国产麻豆网| 99久久人妻综合| 美女国产视频在线观看| 日本91视频免费播放| 国产一区亚洲一区在线观看| 天天影视国产精品| xxx大片免费视频| 久久狼人影院| 亚洲欧洲日产国产| 久久精品aⅴ一区二区三区四区 | 丰满饥渴人妻一区二区三| 亚洲一区二区三区欧美精品| 国产黄色视频一区二区在线观看| 久久久亚洲精品成人影院| 亚洲欧美日韩另类电影网站| 免费人妻精品一区二区三区视频| 免费看av在线观看网站| 成人黄色视频免费在线看| 亚洲国产成人一精品久久久| 免费播放大片免费观看视频在线观看| 亚洲美女搞黄在线观看| 日本vs欧美在线观看视频| 巨乳人妻的诱惑在线观看| 国产男女内射视频| 亚洲人成网站在线观看播放| 国产高清国产精品国产三级| 亚洲男人天堂网一区| 亚洲 欧美一区二区三区| 看十八女毛片水多多多| 午夜老司机福利剧场| 高清黄色对白视频在线免费看| 丝袜脚勾引网站| 2022亚洲国产成人精品| 亚洲av综合色区一区| 欧美日韩一级在线毛片| 亚洲图色成人| 免费不卡的大黄色大毛片视频在线观看| 免费少妇av软件| 丰满乱子伦码专区| 高清视频免费观看一区二区| 啦啦啦中文免费视频观看日本| 亚洲成色77777| 日本欧美国产在线视频| 1024视频免费在线观看| 成人二区视频| 久久这里只有精品19| 天天影视国产精品| 男女国产视频网站| 欧美黄色片欧美黄色片| 日日爽夜夜爽网站| 国产精品av久久久久免费| 女性被躁到高潮视频| 免费高清在线观看视频在线观看| 老汉色∧v一级毛片| 欧美bdsm另类| 免费少妇av软件| 国产一区二区 视频在线| 国产一区二区激情短视频 | 巨乳人妻的诱惑在线观看| 欧美少妇被猛烈插入视频| 国产成人一区二区在线| 日韩成人av中文字幕在线观看| 伦理电影大哥的女人| 中文字幕精品免费在线观看视频| 男的添女的下面高潮视频| 在线观看www视频免费| 亚洲欧美成人综合另类久久久| 欧美av亚洲av综合av国产av | 91精品国产国语对白视频| 免费女性裸体啪啪无遮挡网站| 国产精品99久久99久久久不卡 | 久久久久久伊人网av| 精品一区二区三卡| 精品第一国产精品| videosex国产| 中文乱码字字幕精品一区二区三区| 欧美日韩成人在线一区二区| 亚洲图色成人| av网站在线播放免费| 中文天堂在线官网| 亚洲精品在线美女| 极品人妻少妇av视频| 国产精品香港三级国产av潘金莲 | 亚洲在久久综合| videossex国产| 成人毛片a级毛片在线播放| 三级国产精品片| 一区二区av电影网| 男女无遮挡免费网站观看| 黄片播放在线免费| 成人亚洲欧美一区二区av| 欧美日韩av久久| 五月伊人婷婷丁香| 久久久a久久爽久久v久久|