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

    融合粗糙集與灰色模型的道岔故障預(yù)測

    2019-03-23 05:59:26張友鵬江雪瑩趙斌
    關(guān)鍵詞:粗糙集道岔故障診斷

    張友鵬,江雪瑩,趙斌

    融合粗糙集與灰色模型的道岔故障預(yù)測

    張友鵬,江雪瑩,趙斌

    (蘭州交通大學(xué) 自動化與電氣工程學(xué)院,甘肅 蘭州 730070)

    以高鐵常用S700K型轉(zhuǎn)轍機(jī)為例,通過預(yù)測電流特征實現(xiàn)對未發(fā)生故障的預(yù)先檢修。結(jié)合粗糙集與灰色理論提出一種新方法用于道岔故障預(yù)測,該方法通過粗糙集的知識獲取與規(guī)則約簡,獲得最小診斷規(guī)則;通過離散灰色預(yù)測模型的建模方法,實時建立不同故障類型的預(yù)測模型。隨機(jī)抽取30組故障進(jìn)行診斷,其中96.67%與實際情況相符,可滿足診斷準(zhǔn)確率的要求;隨機(jī)抽取1組預(yù)測情況,預(yù)測值與實際值之間的殘差較小,可滿足預(yù)測準(zhǔn)確率的要求。所提出方法能夠有效減少由故障帶來的安全問題,可操作性高,更具實用性。

    道岔;故障預(yù)測;粗糙集;離散灰色預(yù)測模型

    道岔是鐵路信號設(shè)備中最重要同時也是極易發(fā)生故障的設(shè)備[1],其安全性和可靠性直接關(guān)系鐵路的正常運營,因此,對于道岔的故障診斷是十分關(guān)鍵且必要的。經(jīng)過現(xiàn)場調(diào)研發(fā)現(xiàn),道岔故障常常是造成重大事故的主要原因之一,一旦發(fā)生故障,即可能造成脫軌、追尾等險性事故[2]。因此,研究道岔故障預(yù)測具有重要意義,可在故障發(fā)生前,根據(jù)預(yù)測結(jié)果采取針對性措施,確保道岔安全可靠運行。研究表明[3?4],道岔動作電流可以清晰反映道岔當(dāng)前所處狀態(tài),通過預(yù)測,分析道岔動作電流的變化,能得出道岔的變化趨勢。近些年,基于Adaboost算法和神經(jīng)網(wǎng)絡(luò)等智能技術(shù)的故障預(yù)測,均取得了一定的成果[5?6],但由于Adaboost算法迭代次數(shù)難以確定且訓(xùn)練時間長;神經(jīng)網(wǎng)絡(luò)的收斂速度慢有時還會發(fā)散,且需要大量的歷史數(shù)據(jù),因此,需要尋找一種更為優(yōu)越的方法。離散灰色模型所需數(shù)據(jù)少,計算量小、預(yù)測精度高,是一種較為理想的預(yù)測模型[7]。然而,該模型并不能反映與故障類型的關(guān)系,因此,不能單獨用于故障預(yù)測。研究發(fā)現(xiàn)[8?9],道岔故障類型與其動作電流之間存在很強(qiáng)的對應(yīng)關(guān)系,但通過動作電流獲取的數(shù)據(jù)并不能夠很好地滿足現(xiàn)有道岔故障診斷方法所需數(shù)據(jù)量的要求,使得診斷結(jié)果的準(zhǔn)確率較低,而通過粗糙集知識獲取方法,能夠在數(shù)據(jù)量較少的情況下,更加準(zhǔn)確地對故障進(jìn)行分類。因此,提出基于粗糙集與離散灰色模型的道岔故障預(yù)測方法,利用粗糙集建立決策表,并對電流特征與故障類型間的規(guī)則進(jìn)行簡化;再建立模型對電流特征進(jìn)行預(yù)測,獲得預(yù)測特征后,對照規(guī)則得出最終預(yù)測結(jié)果。該方法結(jié)合預(yù)測的動作電流特征和診斷規(guī)則,盡早發(fā)現(xiàn)還未發(fā)生的故障,能預(yù)先有針對性地對故障處進(jìn)行檢修。

    1 基于粗糙集理論的知識提取與規(guī)則約簡

    粗糙集理論由Pawlak提出[10],該理論能夠有效地對信息進(jìn)行分析與處理,并從中發(fā)現(xiàn)隱含的知識,揭示潛在的規(guī)律[11],尤其適用于數(shù)據(jù)量較少的故障診斷。

    粗糙集理論能有效處理不精確、不一致的信息數(shù)據(jù),發(fā)現(xiàn)內(nèi)部潛在關(guān)系,得到有用知識[12],并把研究對象抽象為一個信息系統(tǒng),在這個信息系統(tǒng)中存在屬性與值[13],而這種屬性與值的關(guān)系構(gòu)成一張二維表,稱之為決策表。

    決策表的生成、化簡方法如下:

    1) 構(gòu)建。將選出的電流特征作為條件屬性,故障作為決策屬性,列出初步?jīng)Q策表。

    2) “替換”。將每個數(shù)據(jù)與其正常狀態(tài)時的數(shù)據(jù)進(jìn)行比較:大于記為2,等于記為1(將屬于正常值的±0.1范圍內(nèi)記作等于),小于記為0。

    3) 屬性約簡。刪除多余的條件屬性,得到簡化決策表。

    4) 規(guī)則生成。簡化決策表的條件屬性為規(guī)則前件,決策屬性為規(guī)則后件,生成規(guī)則。

    2 離散灰色模型的預(yù)測原理

    離散灰色模型(Discrete Grey Model, DGM)是以微分方程為核心的建模方法,其根據(jù)已有信息,通過對初始數(shù)據(jù)的累加、累減等計算建立模型,利用提供的已知或不確定的數(shù)據(jù)來計算各參量,并對數(shù)據(jù)進(jìn)行預(yù)測。

    DGM(1,1)的定義形式[14]為:

    DGM(1,1)建模及預(yù)測具體步驟[15]如下:

    2) 建立矩陣和,

    本文利用DGM(1,1)對道岔動作電流特征進(jìn)行預(yù)測,根據(jù)預(yù)測結(jié)果,結(jié)合診斷規(guī)則,確定故障類型,為現(xiàn)場道岔故障的預(yù)先檢修提供依據(jù)。

    3 基于粗糙集與DGM(1,1)的道岔故障預(yù)測方法

    利用粗糙集獲得最小診斷規(guī)則,利用DGM(1,1)獲取高精度的道岔動作電流特征的預(yù)測值。具體預(yù)測過程如圖1所示。

    圖1 具體預(yù)測過程

    采用直接多步預(yù)測,設(shè)定某一時刻動作電流特征預(yù)測點數(shù)目為,通過預(yù)測點建模得到預(yù)測值,但不宜過大,因為DGM(1,1) 長期預(yù)測會出現(xiàn)較大偏差,而短期預(yù)測效果好[16]。每一組預(yù)測點都反映了道岔動作電流在某一時刻的特征,將特征“替換”后,對照診斷規(guī)則,若符合則可知預(yù)測的故障類型。

    其中:b為類間離散距離,如式(12)所示;w為類內(nèi)離散距離,如式(13)所示。

    根據(jù)式(11)~(15),即可選出最優(yōu)特征。

    選出最優(yōu)特征之后,進(jìn)行“替換”。即將特征值用“0,1和2”來替換,“替換”規(guī)則為:大于基準(zhǔn)值替為2,小于替為0,在基準(zhǔn)值正負(fù)0.1范圍內(nèi)替為1(基準(zhǔn)值選取為道岔正常動作電流曲線的特征值)。

    “替換”后,建立決策表,約簡得到最小決策表和最小診斷規(guī)則。

    根據(jù)系統(tǒng)已監(jiān)測到的動作電流數(shù)據(jù),計算其特征值,建立DGM(1,1)預(yù)測模型(模型中每5個實際值為1組進(jìn)行預(yù)測,預(yù)測10次),并對相應(yīng)故障類型特征進(jìn)行預(yù)測。截取其中1組預(yù)測值,結(jié)合“替換”方法,與最小診斷規(guī)則的條件屬性進(jìn)行比對,判斷符合哪一種,即可得出預(yù)測的道岔故障類型。

    4 應(yīng)用實例

    本文主要研究高速鐵路道岔故障,因此提到轉(zhuǎn)轍機(jī)相關(guān)故障時,均為S700K型。由于S700K型為交流轉(zhuǎn)轍機(jī),因此對其三相動作電流進(jìn)行故障預(yù)測。根據(jù)現(xiàn)場調(diào)研可知幾種常見的故障如表1所示。

    表1 常見故障

    以正常動作電流為例,對其進(jìn)行分區(qū):1DQJ吸起后,曲線開始記錄:0?1為啟動區(qū);1?2為解鎖區(qū),2?3為轉(zhuǎn)換區(qū),3?4為鎖閉區(qū);4?5為緩放區(qū),但在其緩放時間內(nèi),啟動電路中仍會有兩相小電流存在;最后1DQJ落下,曲線終止,如圖2 所示。

    按照圖2所示,對其余8種故障類型(2?9)按相應(yīng)時間區(qū)段對電流曲線進(jìn)行區(qū)域劃分,提取時域特征參數(shù),利用Fisher準(zhǔn)則選出最優(yōu)特征分別如下。A相:1,7,2,9,3,9,4,9和5,6;B相:1,8,2,9,3,4,4,9和5,2;C相:1,2,2,9,3,9,4,9和5,6。這里做一個簡單的解釋:A相的1,7就代表A相動作電流曲線在第1個區(qū)段內(nèi)的第9個時域特征參數(shù)所對應(yīng)的值,以此類推,然后通過前文所述的“替換”方法建立決策表,如表2所示。其中,為序號集,為決策屬性集,1?15分別為三相動作電流的最優(yōu)特征所對應(yīng)的時域特征參數(shù)經(jīng)過“替換”之后得到的狀態(tài)特征,這些狀態(tài)特征的集合稱為條件屬性集。例如:故障8的3,其對應(yīng)的時域特征參數(shù)數(shù)值為1.01,經(jīng)過與正常狀態(tài)的對比,“替換”為1,即3所對應(yīng)的值。

    分別去掉1?15,對比發(fā)現(xiàn),去掉其中某些條件屬性會對診斷結(jié)果有影響,根據(jù)粗糙集規(guī)則簡化方法,可得最小診斷規(guī)則如表3所示。

    從表3可以看出,原來需要15個條件屬性得到的決策屬性,在經(jīng)過規(guī)則簡化之后僅需4個條件屬性即可得到,如第8行獲取的規(guī)則:3(1)∩5(2)∩10(0)∩11(2)T(f8)。

    由于條件屬性的減少,減少了后續(xù)故障預(yù)測的復(fù)雜性。從原有故障數(shù)據(jù)集中抽取30組進(jìn)行診斷準(zhǔn)確率的驗證,結(jié)果表明,其中96.67%是符合實際情況的。

    圖2 正常(f1)動作電流分區(qū)示意

    表2 決策表

    表3 最小診斷規(guī)則

    武漢鐵路局電務(wù)段管轄內(nèi),設(shè)備名稱為1J1的道岔于2015?02?06,12:23:03發(fā)生道岔表示整流匣開路故障,取其故障前15次動作的電流數(shù)據(jù),將這15次中的第1~5次作為一組初始序列,分別對3,5,10和11建立模型,建立模型分別如式(16)~(19)所示。

    3的DGM(1,1):

    5的DGM(1,1):

    10的DGM(1,1):

    11的DGM(1,1):

    利用式(16)~(19)所建立的模型以及第1~5次的實際動作電流數(shù)據(jù)所提取的時域特征參數(shù),對后續(xù)電流特征參數(shù)進(jìn)行10步預(yù)測,得到10次預(yù)測的電流特征參數(shù)后,對其進(jìn)行“替換”,然后根據(jù)“替換”后的狀態(tài)特征結(jié)合最小診斷規(guī)則,得到最后預(yù)測的道岔故障類型。

    將第1次至第5次的實際值與其相應(yīng)預(yù)測值(這里的預(yù)測值其實相當(dāng)于是模擬值,即根據(jù)遞推公式以及還原值公式推出來的值)進(jìn)行對比,并得出殘差(模擬誤差),如表4所示。

    從表4可以看出,殘差(%)均較小,說明預(yù)測值與實際值之間的差別較小,滿足預(yù)測的準(zhǔn)確率要求。

    根據(jù)所建立模型計算預(yù)測故障類型的時域特征參數(shù),對其進(jìn)行“替換”,并最終得出故障預(yù)測結(jié)果如表5所示,其中預(yù)測次數(shù)共為10次,即故障發(fā)生前的10次道岔動作電流特征的預(yù)測,但需要說明的一點是:預(yù)測出的并不是準(zhǔn)確值,而是接近于準(zhǔn)確值的一種趨勢。

    表4 實際值與預(yù)測值對比

    表5 故障預(yù)測數(shù)據(jù)結(jié)果

    從表5可以看出,第1~3次預(yù)測時,故障代碼分別為“1211”和“1212”,根據(jù)最小診斷規(guī)則(表3)可知無此種代碼對應(yīng)的故障類型;當(dāng)?shù)?~10次動作時,故障代碼為“1202”,由表3可知故障類型為8,因而表示電路故障的概率為70%,無對應(yīng)故障代碼概率為30%。

    由于無對應(yīng)故障代碼概率低于表示電路故障概率,因此,最終此次預(yù)測結(jié)果為表示電路故障(未構(gòu)通)。預(yù)測結(jié)果與實際情況相符。

    若通過早期動作電流特征直接進(jìn)行故障診斷,則診斷結(jié)果為60%正常,40%無對應(yīng)故障代碼,不能正確及時的診斷出故障,如表6所示。這是由于早期電流特征不足以直接準(zhǔn)確診斷出道岔的故障,而通過預(yù)測則可以獲得更多后期可能出現(xiàn)的特征,幫助故障的早期識別。

    為進(jìn)一步驗證本文方法的有效性,試驗分析了117個實例,其中110例符合實際情況,預(yù)測所得故障中,概率最高的特征所對應(yīng)的故障即為實際發(fā)生的故障,預(yù)測準(zhǔn)確率達(dá)到94.02%,而在預(yù)測前同時刻基于粗糙集故障診斷方法的準(zhǔn)確率僅為17.09%。

    據(jù)此對比,充分說明,在還未發(fā)生故障之時直接利用故障診斷方法是無法有效準(zhǔn)確地得出診斷結(jié)果,而經(jīng)過預(yù)測后續(xù)特征的方法得出可能出現(xiàn)的結(jié)果,再對預(yù)測結(jié)果進(jìn)行故障診斷能夠更為高效、準(zhǔn)確的針對還未發(fā)生故障的設(shè)備進(jìn)行檢修。

    部分較為典型的實例如表7所示。對于表7中提出的任一故障,其中實際值特征代碼均來自于用于預(yù)測的5次實際動作中的最后一次。預(yù)測次數(shù)為10次。所得故障類型用表1中各故障類型的對應(yīng)代碼表示。

    表6 早期故障診斷

    表7 實例分析表

    從表7可以看出,在早期故障未發(fā)生時,利用粗糙集的診斷方法直接對故障進(jìn)行診斷,準(zhǔn)確率是相當(dāng)?shù)偷?,原因就在于早期故障特征并不明顯,不能準(zhǔn)確作出診斷。

    5 結(jié)論

    1) 基于粗糙集的知識獲取與規(guī)則化簡方法,得到用于診斷故障的最小決策表及最小診斷規(guī)則,能夠更好地反映電流與故障之間的關(guān)系,提高了診斷的準(zhǔn)確率。

    2) 選用的樣本是從眾多實例中選取的典型案例,因而獲取的規(guī)則可信度較高。

    3) 本文所提方法結(jié)合預(yù)測電流特征與最小診斷規(guī)則,能盡早且準(zhǔn)確地發(fā)現(xiàn)還未發(fā)生的故障,根據(jù)預(yù)測結(jié)果,有針對地對其進(jìn)行檢修,可有效減少故障的發(fā)生,避免故障真實發(fā)生所帶來的安全性問題。

    4) 由于電流數(shù)據(jù)可隨道岔動作獲得,因此規(guī)則庫的建立是容易的;又因為預(yù)測模型的建立是隨數(shù)據(jù)而變換的,因此是實時建立的,更具實用性。

    [1] 林海香, 李陽慶. 面向高速鐵路道岔的隱患分析與故障定位改進(jìn)[J]. 鐵道科學(xué)與工程學(xué)報, 2018, 15(9): 2217?2223. LIN Haixiang, LI Yangqing. Hidden-security-risk analysis and fault location improvement for high-speed railway switch[J]. Journal of Railway Science and Engineering, 2018, 15(9): 2217?2223.

    [2] 劉明明, 王良順, 趙輝, 等. 基于群決策的道岔控制電路故障診斷方法研究[J]. 自動化學(xué)報, 2017, 43(5): 1? 10. LIU Mingming, WANG Liangshun, ZHAO Hui, et al. Fault diagnosis method of turnout control circuit based on group decision making[J]. Acta Automatica Sinica, 2017, 43(5): 1?10.

    [3] 張星. 基于相似度的道岔故障電流曲線識別方法研究與實現(xiàn)[D]. 蘭州: 蘭州交通大學(xué), 2016: 10?15. ZHANG Xing. Research and implementation of switch fault current curve recognition method based on similarity[D]. Lanzhou: Lanzhou Jiaotong University, 2016: 10?15.

    [4] 王鐵軍, 董昱, 馬彩霞, 等. 基于BP神經(jīng)網(wǎng)絡(luò)的道岔智能故障診斷技術(shù)[J]. 鐵道運營技術(shù), 2011, 17(2): 4?7. WANG Tiejun, DONG Yu, MA Caixia, et al. Fault diagnosis of turnout based on BP neural network[J]. Railway Operation Technology, 2011, 17(2): 4?7.

    [5] Yilboga H, Eker O F, Gü?lü A, et al. Failure prediction on railway turnouts using time delay neural networks[C]// Computational Intelligence for Measurement Systems and Applications (CIMSA), Taranto: IEEE Xplore, 2010: 1521?1524.

    [6] WANG Guang, XU Tianhua, WANG Haifeng, et al. AdaBoost and least square based failure prediction of railway turnouts[C]// 2016 9th International Symposium on Computational Intelligence and Design (ISCID), Hangzhou: IEEE Xplore, 2016: 2488?2491.

    [7] 丁松, 黨耀國, 徐寧, 等. 基于時滯效應(yīng)的多變量離散灰色預(yù)測模型[J]. 控制與決策, 2017, 32(11): 1997? 2004. DING Song, DANG Yaoguo, XU Ning, et al. Multi- variable time-delayed discrete grey model[J]. Control and Decision, 2017, 32(11): 1997?2004.

    [8] 趙林海, 陸橋. 基于灰關(guān)聯(lián)的道岔故障診斷方法[J]. 鐵道學(xué)報, 2014, 36(2): 69?74. ZHAO Linhai, LU Qiao. Fault diagnosis method of turnout based on grey relation[J]. Journal of the China Railway Society, 2014, 36(2): 69?74.

    [9] 程宇佳. 基于核方法的高速鐵路道岔故障診斷方法研究[D]. 北京: 北京交通大學(xué), 2016: 3?9. CHENG Yujia. Research on fault diagnosis method of high-speed railway turnout based on kernel method[D]. Beijing: Beijing Jiaotong University, 2016: 3?9.

    [10] 黃文濤, 趙學(xué)增, 王偉杰, 等. 基于粗糙集理論的故障診斷決策規(guī)則提取方法[J]. 中國電機(jī)工程學(xué)報, 2003, 23(11): 150?154. HUANG Wentao, ZHAO Xuezeng, WANG Weijie, et al. Fault diagnosis decision rules extraction method based on rough set theory[J]. Journal of Chinese Electrical Engineering, 2003, 23(11): 150?154.

    [11] 費勝巍, 孫宇. 融合粗糙集與灰色理論的電力變壓器故障預(yù)測[J]. 中國電機(jī)工程學(xué)報, 2008, 28(16): 154? 160. FEI Shengwei, SUN Yu. Fault prediction of power transformer by combination of rough sets and grey theory[J]. Journal of Chinese Electrical Engineering, 2008, 28(16): 154?160.

    [12] 張明. 粗糙集理論中的知識獲取與約簡方法的研究[D].南京: 南京理工大學(xué), 2012: 7?13. ZHANG Ming. Research on knowledge acquisition and reduction in rough set theory[D]. Nanjing: Nanjing University of Science & Technology, 2012: 7?13.

    [13] 朱敏玲, 吳海艋, 石磊. 粗糙集規(guī)則匹配算法及其在文本分類中的應(yīng)用[J]. 計算機(jī)系統(tǒng)應(yīng)用, 2018, 27(4): 131?137. ZHU Minling, WU Haimeng, SHI Lei. Rough set rule matching algorithm and its application in text categorization[J]. Computer Systems & Applications, 2018, 27(4): 131?137.

    [14] 蘇先娜, 謝富紀(jì). DGM(1,1)模型的特性及其在科技創(chuàng)新領(lǐng)域中的應(yīng)用[J]. 系統(tǒng)工程理論與實踐, 2016, 36(3): 635?641. SU Xianna, XIE Fuji. The properties of model DGM(1,1) and its application in technology innovation[J]. Systems Engineering-Theory & Practice, 2016, 36(3): 635?641.

    [15] 謝乃明,張可.離散灰色預(yù)測模型及其應(yīng)用[M]. 北京: 科學(xué)出版社, 2016: 51?57. XIE Naiming, ZHANG Ke. Discrete grey forecasting model and its application[M]. Beijing: Science Press, 2016: 51?57.

    [16] WU Lifeng, WANG Yinao. Modelling DGM(1,1) under the criterion of the minimization of mean absolute percentage error[C]// 2009 Second International Symposium on Knowledge Acquisition and Modeling. Wuhan: IEEE Xplore, 2009: 123?126.

    [17] 武漢鐵路局電務(wù)處. 信號集中監(jiān)測信息分析指南[M].北京: 中國鐵道出版社, 2016: 24?26. Wuhan Railway Bureau Electric Service. Analysis of signal centralized monitoring information[M]. Beijing: China Railway Press, 2016: 24?26.

    Turnout fault prediction based on Rough Set and Grey Model

    ZHANG Youpeng, JIANG Xueying, ZHAO Bin

    (School of Automatic & Electric Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China)

    Taking S700K switch machine which is often used in high speed railway as an example, by predicting current characteristics, the fault of turnout which has not happened can be examined and repaired in advance. A new method for fault prediction of turnout based on Rough Set and Grey Theory was presented. And the minimum diagnosis rules were obtained by knowledge acquisition and rule reduction methods which are based on Rough Set. Then through modeling method, discrete grey prediction models of different fault types were established in real time. Thirty groups of faults were randomly selected for diagnosis, and 96.67% of them were consistent with the actual situation, which could meet the requirements of diagnostic accuracy. A group of cases were randomly selected for prediction, and it could be known that the residual between the predictive value and the actual value is small, which could meet the requirement of predictive accuracy. The proposed method can effectively reduce accident rate caused by faults, and also have high operability and practicability.

    turnout; fault prediction; Rough Set; Discrete Grey Prediction Models

    TP206+.3;U213.6

    A

    1672 ? 7029(2019)09? 2331 ? 08

    10.19713/j.cnki.43?1423/u.2019.09.027

    2018?11?22

    中國鐵路總公司科技研究開發(fā)計劃項目(2016J006-A,2015X007-H);蘭州交通大學(xué)青年科學(xué)基金資助項目(2017052)

    張友鵬(1965?),男,甘肅慶陽人,教授,從事故障診斷及預(yù)測研究;E?mail:Jiangxy_105@163.com

    (編輯 陽麗霞)

    猜你喜歡
    粗糙集道岔故障診斷
    基于Pawlak粗糙集模型的集合運算關(guān)系
    中低速磁浮道岔與輪軌道岔的差異
    場間銜接道岔的應(yīng)用探討
    既有線站改插鋪臨時道岔電路修改
    多?;植诩再|(zhì)的幾個充分條件
    雙論域粗糙集在故障診斷中的應(yīng)用
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    兩個域上的覆蓋變精度粗糙集模型
    KJH101-127型氣動司控道岔的改造
    基于LCD和排列熵的滾動軸承故障診斷
    国产 精品1| 天天躁夜夜躁狠狠躁躁| 午夜免费观看性视频| 叶爱在线成人免费视频播放| 精品一区二区三区四区五区乱码 | 久久久精品免费免费高清| 成人亚洲精品一区在线观看| 18在线观看网站| 高清av免费在线| 国语对白做爰xxxⅹ性视频网站| 晚上一个人看的免费电影| 亚洲三级黄色毛片| 午夜福利在线免费观看网站| 精品国产一区二区久久| 大片免费播放器 马上看| 青春草国产在线视频| 免费少妇av软件| 五月开心婷婷网| 大香蕉久久成人网| 日本爱情动作片www.在线观看| 韩国高清视频一区二区三区| kizo精华| 两性夫妻黄色片| 丰满乱子伦码专区| 建设人人有责人人尽责人人享有的| 成年动漫av网址| 免费观看在线日韩| 一边亲一边摸免费视频| 国产成人aa在线观看| 欧美精品高潮呻吟av久久| 色网站视频免费| 搡女人真爽免费视频火全软件| 一区二区三区四区激情视频| 狠狠精品人妻久久久久久综合| 国产成人一区二区在线| 91国产中文字幕| av天堂久久9| 久久精品夜色国产| 伊人久久国产一区二区| 欧美黄色片欧美黄色片| 老鸭窝网址在线观看| 街头女战士在线观看网站| 激情五月婷婷亚洲| 秋霞在线观看毛片| 欧美老熟妇乱子伦牲交| 美女主播在线视频| 久久 成人 亚洲| 欧美成人精品欧美一级黄| 欧美人与善性xxx| 大香蕉久久成人网| 在线观看免费高清a一片| 纯流量卡能插随身wifi吗| 精品人妻熟女毛片av久久网站| 亚洲人成电影观看| 看免费av毛片| 亚洲av在线观看美女高潮| 90打野战视频偷拍视频| 成人手机av| 成人漫画全彩无遮挡| 一二三四在线观看免费中文在| 免费日韩欧美在线观看| 亚洲一码二码三码区别大吗| 日本免费在线观看一区| 午夜免费观看性视频| 精品人妻在线不人妻| 国产伦理片在线播放av一区| av有码第一页| 老鸭窝网址在线观看| 色94色欧美一区二区| 交换朋友夫妻互换小说| av免费观看日本| 熟女少妇亚洲综合色aaa.| 夫妻午夜视频| 中文字幕精品免费在线观看视频| 日韩中文字幕欧美一区二区 | 亚洲第一区二区三区不卡| 亚洲av欧美aⅴ国产| 香蕉国产在线看| 日韩三级伦理在线观看| 在线观看美女被高潮喷水网站| 9191精品国产免费久久| 亚洲五月色婷婷综合| 黑人巨大精品欧美一区二区蜜桃| 亚洲三级黄色毛片| 免费大片黄手机在线观看| 午夜福利网站1000一区二区三区| 亚洲天堂av无毛| 国产成人91sexporn| 国产亚洲最大av| 卡戴珊不雅视频在线播放| 亚洲一区中文字幕在线| 老汉色∧v一级毛片| 国产在线视频一区二区| 亚洲国产av影院在线观看| 咕卡用的链子| 一区福利在线观看| 一区二区三区精品91| 免费观看性生交大片5| 18禁观看日本| 亚洲av日韩在线播放| 女人精品久久久久毛片| 一级a爱视频在线免费观看| 亚洲精品一二三| 成人手机av| 一区二区三区四区激情视频| 免费少妇av软件| 精品人妻偷拍中文字幕| 黄色怎么调成土黄色| 美女视频免费永久观看网站| 欧美日韩一区二区视频在线观看视频在线| 欧美老熟妇乱子伦牲交| 99热国产这里只有精品6| 国产成人a∨麻豆精品| 久久久久久久大尺度免费视频| 91精品三级在线观看| 成人黄色视频免费在线看| 久久精品久久精品一区二区三区| 久久久a久久爽久久v久久| 午夜福利网站1000一区二区三区| 精品少妇久久久久久888优播| 国产精品蜜桃在线观看| 黄频高清免费视频| www.熟女人妻精品国产| 三上悠亚av全集在线观看| 久久精品人人爽人人爽视色| 丰满乱子伦码专区| 亚洲av在线观看美女高潮| 永久网站在线| 国产一区二区在线观看av| 一本大道久久a久久精品| 亚洲精品久久成人aⅴ小说| 久久精品亚洲av国产电影网| 99九九在线精品视频| 国产精品久久久av美女十八| 超色免费av| 看免费av毛片| 秋霞伦理黄片| 视频区图区小说| 日韩欧美精品免费久久| 亚洲国产看品久久| av福利片在线| 国产片特级美女逼逼视频| 成人午夜精彩视频在线观看| 国产国语露脸激情在线看| 欧美 亚洲 国产 日韩一| 亚洲精品久久成人aⅴ小说| 久久99蜜桃精品久久| 日日摸夜夜添夜夜爱| videos熟女内射| 成人毛片60女人毛片免费| 青春草亚洲视频在线观看| 国产精品久久久久久久久免| 熟女电影av网| www.精华液| 狂野欧美激情性bbbbbb| 不卡视频在线观看欧美| 国产精品免费大片| 在线亚洲精品国产二区图片欧美| 成人毛片a级毛片在线播放| 搡老乐熟女国产| 纯流量卡能插随身wifi吗| 久久人妻熟女aⅴ| 久久ye,这里只有精品| 国产成人欧美| 熟女电影av网| 久久久久久久精品精品| 久久国产精品男人的天堂亚洲| 这个男人来自地球电影免费观看 | 少妇熟女欧美另类| 看十八女毛片水多多多| √禁漫天堂资源中文www| 国产精品久久久久久精品古装| 久久久久久久久久久免费av| 午夜久久久在线观看| 激情视频va一区二区三区| 久久久久国产网址| 色哟哟·www| 亚洲av综合色区一区| 黑丝袜美女国产一区| 国产男女内射视频| 亚洲图色成人| 18在线观看网站| 少妇猛男粗大的猛烈进出视频| 美女午夜性视频免费| 免费高清在线观看视频在线观看| av不卡在线播放| 人妻一区二区av| 精品一区二区免费观看| 青春草国产在线视频| 色网站视频免费| 91久久精品国产一区二区三区| 熟妇人妻不卡中文字幕| 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 亚洲av综合色区一区| 熟妇人妻不卡中文字幕| 男女免费视频国产| 亚洲成人av在线免费| 美女大奶头黄色视频| 国产精品三级大全| 黄色视频在线播放观看不卡| 日韩制服丝袜自拍偷拍| 欧美人与善性xxx| 2022亚洲国产成人精品| √禁漫天堂资源中文www| 春色校园在线视频观看| 伊人久久大香线蕉亚洲五| 尾随美女入室| 人成视频在线观看免费观看| 亚洲一级一片aⅴ在线观看| 成人18禁高潮啪啪吃奶动态图| av不卡在线播放| 亚洲欧美色中文字幕在线| 欧美激情高清一区二区三区 | 国产免费福利视频在线观看| 亚洲av男天堂| 一二三四中文在线观看免费高清| 亚洲精品国产av蜜桃| 亚洲精品美女久久久久99蜜臀 | 午夜av观看不卡| 人妻人人澡人人爽人人| 亚洲国产日韩一区二区| 精品第一国产精品| 午夜免费观看性视频| 久久久精品国产亚洲av高清涩受| 国产一级毛片在线| 亚洲av中文av极速乱| av国产久精品久网站免费入址| 亚洲欧美成人精品一区二区| 免费久久久久久久精品成人欧美视频| videos熟女内射| 久久精品亚洲av国产电影网| 最近2019中文字幕mv第一页| 波野结衣二区三区在线| av在线播放精品| 久热这里只有精品99| 久久精品久久久久久噜噜老黄| www日本在线高清视频| 不卡视频在线观看欧美| 丰满迷人的少妇在线观看| 在线观看人妻少妇| 波多野结衣一区麻豆| 下体分泌物呈黄色| 中文字幕另类日韩欧美亚洲嫩草| xxxhd国产人妻xxx| 蜜桃在线观看..| 色婷婷av一区二区三区视频| 中文字幕av电影在线播放| 美女国产高潮福利片在线看| 中文乱码字字幕精品一区二区三区| 春色校园在线视频观看| 80岁老熟妇乱子伦牲交| 看非洲黑人一级黄片| 国语对白做爰xxxⅹ性视频网站| 久久精品熟女亚洲av麻豆精品| 亚洲精品aⅴ在线观看| 国产片内射在线| 香蕉精品网在线| 国产有黄有色有爽视频| 青春草亚洲视频在线观看| 999久久久国产精品视频| 黄频高清免费视频| 亚洲国产av影院在线观看| 水蜜桃什么品种好| 国产在线视频一区二区| 日日撸夜夜添| 亚洲欧美中文字幕日韩二区| 秋霞伦理黄片| 亚洲国产精品一区二区三区在线| 国产精品偷伦视频观看了| 国产精品欧美亚洲77777| 欧美成人午夜精品| 午夜免费男女啪啪视频观看| 成年美女黄网站色视频大全免费| 又粗又硬又长又爽又黄的视频| 亚洲精品日本国产第一区| 久久久国产欧美日韩av| 成人毛片60女人毛片免费| 免费久久久久久久精品成人欧美视频| 一二三四在线观看免费中文在| 免费在线观看完整版高清| 欧美精品一区二区免费开放| 综合色丁香网| 丝瓜视频免费看黄片| 久久人妻熟女aⅴ| 有码 亚洲区| 亚洲av综合色区一区| 日韩一本色道免费dvd| 日韩免费高清中文字幕av| 不卡视频在线观看欧美| 搡老乐熟女国产| av在线老鸭窝| 日韩伦理黄色片| 嫩草影院入口| 亚洲在久久综合| 大片电影免费在线观看免费| 天天躁狠狠躁夜夜躁狠狠躁| 一区二区av电影网| 18禁动态无遮挡网站| 亚洲精品日韩在线中文字幕| 欧美日韩av久久| 亚洲成av片中文字幕在线观看 | 亚洲一区中文字幕在线| 天堂中文最新版在线下载| 国产免费视频播放在线视频| 在线天堂中文资源库| 亚洲国产精品一区三区| 欧美变态另类bdsm刘玥| 多毛熟女@视频| 巨乳人妻的诱惑在线观看| 性少妇av在线| 亚洲图色成人| 日本vs欧美在线观看视频| 三上悠亚av全集在线观看| 男人操女人黄网站| 亚洲精品中文字幕在线视频| 街头女战士在线观看网站| 亚洲欧美日韩另类电影网站| 久久免费观看电影| 欧美日韩精品网址| 一级,二级,三级黄色视频| 欧美 亚洲 国产 日韩一| 久久国产精品大桥未久av| 中文字幕av电影在线播放| 在线精品无人区一区二区三| 在线观看美女被高潮喷水网站| kizo精华| 亚洲情色 制服丝袜| 考比视频在线观看| 久久久精品94久久精品| 亚洲视频免费观看视频| 国产一区二区 视频在线| 国产又爽黄色视频| 爱豆传媒免费全集在线观看| 免费少妇av软件| 亚洲精品乱久久久久久| 亚洲国产精品国产精品| 精品国产一区二区久久| 精品人妻一区二区三区麻豆| 成人国产av品久久久| 成年女人在线观看亚洲视频| 大片电影免费在线观看免费| 亚洲精品久久成人aⅴ小说| xxx大片免费视频| 久久久a久久爽久久v久久| 美女国产视频在线观看| 18禁动态无遮挡网站| 一级毛片黄色毛片免费观看视频| 不卡视频在线观看欧美| 青草久久国产| 久久国产亚洲av麻豆专区| 日韩制服丝袜自拍偷拍| 久久久久久久久久久久大奶| 啦啦啦在线免费观看视频4| 国产乱人偷精品视频| 亚洲国产精品一区二区三区在线| 夫妻午夜视频| 国产精品久久久久久久久免| 永久网站在线| 中国国产av一级| 久热久热在线精品观看| 久久精品人人爽人人爽视色| 又粗又硬又长又爽又黄的视频| 欧美亚洲日本最大视频资源| 国产免费现黄频在线看| 亚洲精品乱久久久久久| 999精品在线视频| 最近中文字幕2019免费版| av网站在线播放免费| 久久精品久久久久久久性| 国产av码专区亚洲av| 丰满迷人的少妇在线观看| 精品亚洲成国产av| 国产精品无大码| 一本久久精品| 欧美黄色片欧美黄色片| av片东京热男人的天堂| 老女人水多毛片| 99热网站在线观看| 天天影视国产精品| 一边亲一边摸免费视频| 免费看不卡的av| 黑人猛操日本美女一级片| 精品一品国产午夜福利视频| 久久久久人妻精品一区果冻| 不卡视频在线观看欧美| 亚洲成人手机| 久久人妻熟女aⅴ| 老司机影院成人| 亚洲一码二码三码区别大吗| 亚洲精品第二区| 在线观看免费日韩欧美大片| 18禁动态无遮挡网站| 久久久久精品久久久久真实原创| 国产亚洲精品第一综合不卡| 久久久国产一区二区| 久久精品aⅴ一区二区三区四区 | 久久久精品区二区三区| 青草久久国产| 久久人人97超碰香蕉20202| 精品久久久久久电影网| 欧美av亚洲av综合av国产av | 热99国产精品久久久久久7| 久久精品国产亚洲av高清一级| 亚洲一区中文字幕在线| 亚洲人成电影观看| 纯流量卡能插随身wifi吗| 久久久精品国产亚洲av高清涩受| 久久久久久久亚洲中文字幕| 欧美少妇被猛烈插入视频| 少妇的逼水好多| 天天躁夜夜躁狠狠躁躁| 日本wwww免费看| 在线观看人妻少妇| 亚洲欧洲日产国产| a级片在线免费高清观看视频| 两个人看的免费小视频| 亚洲国产日韩一区二区| 一级毛片我不卡| 亚洲欧美精品自产自拍| 久久精品夜色国产| 国产熟女午夜一区二区三区| 老熟女久久久| av网站在线播放免费| 国产午夜精品一二区理论片| 国产综合精华液| 国产一区二区激情短视频 | 色94色欧美一区二区| 一区二区三区乱码不卡18| 日本-黄色视频高清免费观看| 在线 av 中文字幕| 免费人妻精品一区二区三区视频| 亚洲国产成人一精品久久久| 18禁裸乳无遮挡动漫免费视频| 国产精品香港三级国产av潘金莲 | 日韩三级伦理在线观看| 少妇被粗大的猛进出69影院| 在线观看国产h片| 亚洲国产日韩一区二区| 最近手机中文字幕大全| 1024视频免费在线观看| 香蕉精品网在线| 精品一区二区免费观看| 成年动漫av网址| 两个人看的免费小视频| 波野结衣二区三区在线| av女优亚洲男人天堂| 欧美日韩成人在线一区二区| www.自偷自拍.com| 在线观看三级黄色| 欧美激情 高清一区二区三区| 一区二区三区四区激情视频| 欧美日韩一区二区视频在线观看视频在线| 最近2019中文字幕mv第一页| 欧美人与性动交α欧美精品济南到 | 国产黄色视频一区二区在线观看| 五月开心婷婷网| videossex国产| 国产成人精品无人区| 三级国产精品片| 亚洲熟女精品中文字幕| 纯流量卡能插随身wifi吗| 91aial.com中文字幕在线观看| 秋霞在线观看毛片| 大话2 男鬼变身卡| 9191精品国产免费久久| 人妻 亚洲 视频| 99香蕉大伊视频| 国产无遮挡羞羞视频在线观看| 在线观看免费日韩欧美大片| 99久国产av精品国产电影| 国产精品av久久久久免费| 91精品伊人久久大香线蕉| 久久人妻熟女aⅴ| 亚洲五月色婷婷综合| av福利片在线| 日韩视频在线欧美| 久久热在线av| 老司机影院毛片| 国产老妇伦熟女老妇高清| 国产黄色免费在线视频| 亚洲精品国产av蜜桃| 久久久久久伊人网av| 国产在线一区二区三区精| 免费在线观看完整版高清| 中文天堂在线官网| 99久久精品国产国产毛片| 成年美女黄网站色视频大全免费| 欧美人与性动交α欧美软件| 成人亚洲精品一区在线观看| 亚洲 欧美一区二区三区| 成年人免费黄色播放视频| 最近中文字幕高清免费大全6| 热99久久久久精品小说推荐| av一本久久久久| a级片在线免费高清观看视频| 2018国产大陆天天弄谢| 欧美av亚洲av综合av国产av | 日本av免费视频播放| 亚洲国产欧美日韩在线播放| 街头女战士在线观看网站| 国产亚洲一区二区精品| 亚洲精品自拍成人| 少妇熟女欧美另类| 欧美日韩成人在线一区二区| 99热全是精品| 久久免费观看电影| 人妻少妇偷人精品九色| 中文天堂在线官网| 秋霞在线观看毛片| 色婷婷av一区二区三区视频| 欧美变态另类bdsm刘玥| 亚洲婷婷狠狠爱综合网| 亚洲av国产av综合av卡| 97在线人人人人妻| 成人18禁高潮啪啪吃奶动态图| 性高湖久久久久久久久免费观看| 久久久久网色| 国产一区二区三区av在线| 精品亚洲成国产av| 久久久久久久久免费视频了| 丝袜美腿诱惑在线| 麻豆乱淫一区二区| 日韩一卡2卡3卡4卡2021年| av视频免费观看在线观看| 毛片一级片免费看久久久久| 亚洲激情五月婷婷啪啪| 伊人久久大香线蕉亚洲五| 国产精品三级大全| 亚洲成色77777| 午夜激情av网站| 亚洲欧美成人综合另类久久久| 精品视频人人做人人爽| 亚洲色图综合在线观看| 寂寞人妻少妇视频99o| 国产无遮挡羞羞视频在线观看| 各种免费的搞黄视频| a级片在线免费高清观看视频| 一级毛片我不卡| 免费人妻精品一区二区三区视频| 波多野结衣av一区二区av| 国产综合精华液| 欧美xxⅹ黑人| 久久久久久人人人人人| 97在线视频观看| 夫妻性生交免费视频一级片| 国产乱来视频区| 国产精品亚洲av一区麻豆 | 久久精品亚洲av国产电影网| 日本91视频免费播放| 久久久久人妻精品一区果冻| 欧美激情高清一区二区三区 | 亚洲国产日韩一区二区| a 毛片基地| 久久精品久久久久久噜噜老黄| 久热这里只有精品99| 王馨瑶露胸无遮挡在线观看| 国产黄色免费在线视频| 女人高潮潮喷娇喘18禁视频| 成人漫画全彩无遮挡| 亚洲第一区二区三区不卡| 亚洲国产欧美日韩在线播放| 有码 亚洲区| 波野结衣二区三区在线| 日韩中文字幕欧美一区二区 | 亚洲国产av影院在线观看| 久久久欧美国产精品| 咕卡用的链子| 男人操女人黄网站| 亚洲精品久久成人aⅴ小说| 大片免费播放器 马上看| 欧美亚洲 丝袜 人妻 在线| 涩涩av久久男人的天堂| 黄色怎么调成土黄色| 欧美精品一区二区大全| 乱人伦中国视频| 国产精品一区二区在线不卡| 欧美激情高清一区二区三区 | 两个人免费观看高清视频| 少妇猛男粗大的猛烈进出视频| 亚洲一码二码三码区别大吗| 亚洲美女搞黄在线观看| 国产亚洲最大av| 免费观看在线日韩| 久久久久久人人人人人| 91久久精品国产一区二区三区| 精品国产露脸久久av麻豆| 性色av一级| 熟女少妇亚洲综合色aaa.| 亚洲人成77777在线视频| 青春草国产在线视频| 在线观看三级黄色| 亚洲成色77777| 亚洲精品在线美女| 国产激情久久老熟女| 国产av精品麻豆| 亚洲av综合色区一区| 国产男女内射视频| 欧美日韩综合久久久久久| 欧美亚洲日本最大视频资源| 国产精品欧美亚洲77777| 在线天堂最新版资源| 久久久久久久大尺度免费视频| av一本久久久久| 午夜福利在线观看免费完整高清在| 飞空精品影院首页| 精品一区在线观看国产| 天堂中文最新版在线下载| 自线自在国产av| 成人18禁高潮啪啪吃奶动态图| 国产福利在线免费观看视频| 亚洲精品一区蜜桃| 日本爱情动作片www.在线观看| 男女边摸边吃奶| 天天影视国产精品|