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

    基于XGBoost建模和Change-Point殘差處理的風(fēng)電機組齒輪箱溫度預(yù)警

    2018-10-13 02:38:18王桂松劉瑞華李韶武
    電力科學(xué)與工程 2018年9期
    關(guān)鍵詞:軸溫變點齒輪箱

    王桂松, 郭 鵬, 胥 佳, 劉瑞華, 李韶武

    (1.華北電力大學(xué) 控制與計算機工程學(xué)院,北京 102206;2.龍源(北京)風(fēng)電工程技術(shù)有限公司,北京 100034)

    0 引言

    風(fēng)電機組長期運行在室外惡劣的自然環(huán)境中,故障率高。齒輪箱、發(fā)電機軸承等大部件的故障檢修難度大、費用高,高額的運營維護成本降低了風(fēng)電的經(jīng)濟效益。通過有效方法,實時監(jiān)測風(fēng)電機組運行狀態(tài),定位狀態(tài)劣化點,回溯狀態(tài)劣化原因,做出預(yù)防性維修,對于減少風(fēng)電場運營成本、降低風(fēng)電機組的運行風(fēng)險具有重要意義[1]。

    齒輪箱是風(fēng)電機組的重要部件,其結(jié)構(gòu)復(fù)雜,包含齒輪、軸承、軸、箱體等多種構(gòu)件,某一結(jié)構(gòu)的失效都有可能引發(fā)齒輪箱的整體故障,進而造成嚴(yán)重的后果[2]。風(fēng)電機組SCADA數(shù)據(jù)中蘊藏著風(fēng)電機組運行狀態(tài)信息,應(yīng)用齒輪箱溫度數(shù)據(jù)分析成為近年來國內(nèi)外風(fēng)電機組齒輪箱狀態(tài)監(jiān)測的主流。BRANDAO.R.F.M等人利用神經(jīng)網(wǎng)絡(luò)的方法預(yù)測風(fēng)電機組齒輪箱油液溫度以實現(xiàn)對齒輪箱的故障預(yù)測[3]。Wang L等使用改進的深度神經(jīng)網(wǎng)絡(luò)進行齒輪箱故障預(yù)警,預(yù)測準(zhǔn)確率已高于支持向量機 (Support Vector Machines,SVM)[4]。但依靠神經(jīng)網(wǎng)絡(luò)等機器學(xué)習(xí)算法建立的監(jiān)測模型存在模型學(xué)習(xí)時間過長,學(xué)習(xí)效率過低問題,不適應(yīng)于在線工程應(yīng)用。郭鵬等人采用非線性狀態(tài)估計建立齒輪箱正常工作狀態(tài)下的溫度模型并用其進行溫度預(yù)測,已具有較高的建模效率和預(yù)測精度[5]。劉帥等人將高斯混合模型、動態(tài)時間規(guī)整及熵權(quán)值算法三者緊密結(jié)合,提出了一種基于群體多維特征相似性的故障預(yù)警策略[6]。祝文穎等人將改進的經(jīng)驗小波變換應(yīng)用于行星齒輪箱振動信號的解調(diào)分析,提出了一種單分量個數(shù)的估算方法,解決了經(jīng)驗小波變換中的Fourier頻譜劃分問題,提高了分析的針對性、準(zhǔn)確率[7]。以上學(xué)者分別利用深度學(xué)習(xí)、機理模型、信號分析的方法進行齒輪箱故障分析,在齒輪箱預(yù)警領(lǐng)域取得了相應(yīng)的成果。但以上方法或存在計算效率問題,或存在故障監(jiān)測時間不確定度問題,即無法給出準(zhǔn)確的狀態(tài)劣化時間點。

    本文使用XGBoost算法建立齒輪箱正常工作狀態(tài)的溫度模型。通過XGBoost模型與其它4類模型對齒輪箱溫度數(shù)據(jù)進行回歸預(yù)測模型實驗,結(jié)果表明XGBoost模型在齒輪箱溫度預(yù)測中綜合性能要優(yōu)于其它4類模型。對于XGBoost模型預(yù)測殘差序列,本文開展了Change-Point算法構(gòu)建,利用Change-Point技術(shù)進行XGBoost模型預(yù)測殘差序列分析并給出變點,可發(fā)現(xiàn)齒輪箱工作異常狀態(tài),定位狀態(tài)劣化時間點,回溯狀態(tài)劣化原因。

    1 XGBoost算法原理

    XGBoost方法為有監(jiān)督學(xué)習(xí)算法,是一種新的梯度提升實現(xiàn)方法,主要解決分類和回歸問題。XGBoost方法自提出以來,廣泛應(yīng)用于金融、醫(yī)療等領(lǐng)域[8-9]。該算法基于“提升”思想,通過將一系列“弱”學(xué)習(xí)機的所有預(yù)測結(jié)果結(jié)合起來,得到“強”的學(xué)習(xí)模型。與傳統(tǒng)的梯度增強回歸樹(Gradient Boosting Decision Tree,GBDT)算法相比,其優(yōu)勢在于XGBoost對于損失函數(shù)進行了二階泰勒展開,將預(yù)測項和正則化項結(jié)合起來,通過簡化目標(biāo)函數(shù)實現(xiàn)計算資源優(yōu)化,以權(quán)衡損失函數(shù)的下降和模型的復(fù)雜程度,并得到最優(yōu)解[10]。同時,XGBoost自動執(zhí)行并行計算,保證了最優(yōu)的計算速度。

    假設(shè)模型有k個決策樹,則模型的輸出表示為

    (1)

    式中:xi為數(shù)據(jù)集;F為決策樹函數(shù)空間;fk對應(yīng)獨立的樹結(jié)構(gòu)和葉子權(quán)重。其損失函數(shù)為:

    (2)

    (3)

    (4)

    式中:gi和hi分別為損失函數(shù)的第一階和第二階梯度統(tǒng)計。

    在移除常數(shù)項后,得到目標(biāo)函數(shù)的簡化形式為:

    (5)

    定義Ij={i|q(xi)=j}為葉子j的實例集,可以將式(5)擴展正則化項Ω為:

    (6)

    對上式進行求導(dǎo)并令求導(dǎo)結(jié)果等于0,可得到ωj的最優(yōu)解

    (7)

    (8)

    正常情況下,枚舉所有可能的樹結(jié)構(gòu)是無法實現(xiàn)的。XGBoost采用貪心算法,每一次對已有的葉子節(jié)點加入一個分割,假設(shè)IL和IR是左右子樹分?jǐn)?shù)分割后的節(jié)點,則分割后的損失函數(shù)為:

    Lsplit=

    (9)

    XGBoost算法為了限值樹的生長,加入了正則項系數(shù)γ,當(dāng)增益大于γ時,允許節(jié)點進行分割,在優(yōu)化目標(biāo)函數(shù)時相當(dāng)于做了預(yù)剪枝。

    2 基于XGBoost算法的風(fēng)電機組齒輪箱溫度建模

    2.1 齒輪箱XGBoost建模參數(shù)選取

    本文采用某風(fēng)電場1.5 MW機組的運行數(shù)據(jù),采樣周期為1分鐘。其切入風(fēng)速為3 m/s,額定風(fēng)速為12 m/s。運用相關(guān)性分析方法和風(fēng)電機組控制機理提取若干與齒輪箱軸溫密切相關(guān)聯(lián)的變量,并增加部分分類變量共同作為模型的輸入,模型的輸出為變量齒輪箱軸溫。相關(guān)系數(shù)的計算公式為:

    (10)

    式中:x為齒輪箱軸溫變量;y為其它變量,包括有功功率、風(fēng)速等與齒輪箱軸溫關(guān)聯(lián)密切的變量;n為數(shù)據(jù)總量。計算該機組其它75個變量與齒輪箱軸溫的相關(guān)性系數(shù),部分關(guān)聯(lián)變量的相關(guān)系數(shù)見表1。

    根據(jù)表1,齒輪箱池油溫、齒輪箱驅(qū)動端油溫和齒輪箱非驅(qū)動端溫升變量為強相關(guān)變量,不能作為模型輸入。選取以下變量作為模型輸入?yún)?shù)。

    (1)功率、風(fēng)速、槳葉槳角、高速端轉(zhuǎn)速和齒輪箱入口油溫為齒輪箱軸溫密切相關(guān)聯(lián)的變量。

    (2)環(huán)境溫度、機艙溫度的不同也會導(dǎo)致齒輪箱軸溫有較大差異。在不同時刻,即使機組的風(fēng)速和功率相同,齒輪箱軸承溫度同樣會受到晝夜溫差和季節(jié)性溫度變化影響。

    表1 相關(guān)性系數(shù)統(tǒng)計表

    (3)風(fēng)電機組運行過程中,當(dāng)齒輪箱軸溫升高,與之相關(guān)聯(lián)的油溫升高,油液粘度下降,導(dǎo)致齒輪箱入口油壓和齒輪箱濾網(wǎng)入口油壓等壓力下降。齒輪箱入口油壓、齒輪箱濾網(wǎng)入口油壓與齒輪箱軸溫為負(fù)相關(guān)。

    (4)風(fēng)機狀態(tài)、可利用率狀態(tài)為SCADA數(shù)據(jù)狀態(tài)編碼,選用風(fēng)機狀態(tài)、可利用率狀態(tài)兩個變量作為模型輸入?yún)?shù)。

    (5)增加是否達(dá)到額定功率(達(dá)到為1,未達(dá)到為0)、是否達(dá)到額定風(fēng)速(達(dá)到為1,未達(dá)到為0)、是否為停機狀態(tài)(停機為1,未停機為0)分類變量。

    2.2 齒輪箱XGBoost建模參數(shù)調(diào)優(yōu)與模型實驗

    本文采用的XGBoost模型調(diào)參選用該風(fēng)場某機組2014年1月份至4月份正常運行數(shù)據(jù)。進行數(shù)據(jù)預(yù)處理后,采用建模數(shù)據(jù)共91 608條,采用測試數(shù)據(jù)共30 067條。模型涉及線程、樹的最大深度、迭代次數(shù)參數(shù)調(diào)優(yōu)。本文利用柵格搜索與單一參數(shù)搜索結(jié)合方法進行參數(shù)尋優(yōu)。Bergstra等[11]研究表明,柵格搜索的效果弱于隨機搜索,采用柵格搜索進行初步組合參數(shù)選取后,固定其它參數(shù)不變,先后對主參數(shù)進行單一參數(shù)搜索以搜尋最優(yōu)參數(shù),進而提高模型預(yù)測精度。在XGBoost建模中,線程參數(shù)決定了算法的CPU使用率,CPU使用率進一步?jīng)Q定模型的計算效率,CPU使用率、計算耗時與線程參數(shù)關(guān)系如圖1所示。在線程參數(shù)達(dá)到8后,CPU使用率達(dá)到100%,且計算耗時基本不變,因此本文將線程參數(shù)設(shè)置為8。

    圖1 CPU使用率、計算耗時與線程關(guān)系圖

    樹的最大深度參數(shù)用來避免過擬合,樹的最大深度參數(shù)越大,模型會學(xué)到更具體更局部的樣本,該參數(shù)會影響模型的預(yù)測精度和計算效率。對于XGBoost模型的預(yù)測精度引入平均絕對誤差(Mean Absolute Error, MAE)、均值偏移誤差(Mean Bias Error, MBE)、根均值平方誤差(Root Mean Squared Error, RMSE)[12]3種殘差分析方法。MAE、MBE和RMSE計算公式如式(11)所示。

    (11)

    通過對齒輪箱軸溫預(yù)測值與實際值之間的殘差統(tǒng)計分析,給出預(yù)測精度。殘差統(tǒng)計、計算耗時與樹的最大深度關(guān)系如圖2所示。

    圖2 殘差統(tǒng)計、計算耗時與樹的最大深度關(guān)系圖

    由圖2可知,在樹的最大深度為8時殘差最小,模型計算耗時并不會隨著樹的深度的增加而呈某種特定趨勢,而是有一個擬合優(yōu)值。本文選擇的樹的最大深度為8,在保證最優(yōu)精度的同時有較高的計算效率。

    弱學(xué)習(xí)器的迭代次數(shù)大小控制算法的擬合優(yōu)度,直接影響XGBboost模型的計算效率。在模型實驗中發(fā)現(xiàn)梯度增強回歸樹(Gradient Boosting Decision Tree, GBDT)、隨機森林(Random Forest, RF)兩類模型同樣含有迭代次數(shù)參數(shù),且RF模型物理內(nèi)存占用受該參數(shù)影響較大。圖3給出了3種模型計算耗時、物理內(nèi)存占用率與迭代次數(shù)參數(shù)選取的關(guān)系。

    圖3 模型計算耗時、物理內(nèi)存占用與迭代次數(shù)關(guān)系圖

    由圖3可知,XGBoost模型計算耗時最少,計算效率最高。RF、XGBoost算法均能實現(xiàn)并行處理,實驗過程中,選用合適的線程參數(shù),以保證RF、XGBoost模型使用處理器的所有CPU。隨著迭代次數(shù)增大,XGBoost、GBDT模型內(nèi)存占用率維持在44%左右;而RF模型內(nèi)存占用率隨著迭代次數(shù)的增加不斷升高,且計算機出現(xiàn)卡頓現(xiàn)象。XGBoost算法本身考慮了當(dāng)數(shù)據(jù)量比較大、內(nèi)存不夠時怎么有效地使用磁盤。如此,在通過大數(shù)據(jù)量樣本建模并監(jiān)控風(fēng)電機組狀態(tài)時,XGBoost可以充分發(fā)揮其優(yōu)勢,提高計算效率。

    為了驗證XGBoost模型對于齒輪箱溫度預(yù)測的適應(yīng)性和有效性。在選取合適的建模參數(shù)后,采用GBDT、人工神經(jīng)網(wǎng)絡(luò)(Artificial Neural Network,ANN)、K-最近鄰(K-Nearest Neighbor,KNN)和RF作為類比模型進行齒輪箱溫度回歸預(yù)測。選取兩臺機組正常樣本數(shù)據(jù)進行模型實驗,實驗樣本分布如表2所示。其中,樣本1和樣本2選用該風(fēng)場不同機組1~4月份1分鐘平均數(shù)據(jù),采用1~3月份數(shù)據(jù)建立模型,應(yīng)用4月份數(shù)據(jù)進行測試。

    5種模型在兩個實驗樣本的最終計算結(jié)果對比如表3所示。從表3中可以發(fā)現(xiàn),在參數(shù)調(diào)優(yōu)后,對比5種模型計算結(jié)果可知,XGBoost模型預(yù)測精度要明顯優(yōu)于KNN、ANN模型;在該數(shù)據(jù)量條件下,與GBDT、RF模型預(yù)測精度相當(dāng),但計算效率更高。結(jié)合XGBoost算法原理,數(shù)據(jù)量越大,XGBoost分類效果會越好。實驗結(jié)果表明,XGBoost模型綜合性能要優(yōu)于其它4類模型,具有預(yù)測精度和計算效率高,資源占用低的特點,證明了XGBoost算法應(yīng)用于齒輪箱溫度預(yù)測的可行性和穩(wěn)定性。

    表2 實驗樣本分布

    表3 5種模型在兩個實驗樣本的計算結(jié)果對比

    2.3 基于XGBoost建模的齒輪箱故障預(yù)警

    利用XGBoost算法建立齒輪箱溫度模型并進行預(yù)測,當(dāng)齒輪箱運行在正常工作狀態(tài)下,預(yù)測殘差較小,預(yù)測精度較高;當(dāng)齒輪箱發(fā)生故障時,數(shù)據(jù)特性將偏離正常工作狀態(tài),預(yù)測殘差明顯增大。將上述建模方法應(yīng)用于某機組7月份齒輪箱溫度傳感器接線故障實例,其中采用該機組1月到5月上旬共194 102條正常數(shù)據(jù)建模,采用包含故障時段的5月下旬到7月上旬共87 020條數(shù)據(jù)進行測試。變量與上述模型實驗選取變量一致,按照上述參數(shù)調(diào)優(yōu)方法,XGBoost預(yù)測結(jié)果如圖4所示。

    圖4 XGBoost預(yù)測結(jié)果

    在6月10號至7月6號出現(xiàn)XGBoost模型預(yù)測殘差過大問題。核對風(fēng)電機組運行日志,該機組在6月27號至7月7號頻繁出現(xiàn)齒輪箱軸溫變化異常并發(fā)出報警。風(fēng)電場運維人員在7月7號進行故障檢修并發(fā)現(xiàn)傳感器接線松動問題,在緊固傳感器接線后齒輪箱軸溫恢復(fù)正常。根據(jù)風(fēng)電機組控制原理,XGBoost模型較傳統(tǒng)閾值限超溫報警可提前發(fā)現(xiàn)齒輪箱軸溫異常。

    3 Change-Point算法與殘差序列分析

    在區(qū)域數(shù)據(jù)平臺通過預(yù)警模型監(jiān)測數(shù)十個風(fēng)場上千臺機組狀態(tài)時,如何實現(xiàn)準(zhǔn)確的故障自動預(yù)警以提前采取相關(guān)聯(lián)的控制動作,是降低風(fēng)電場運維成本的關(guān)鍵。針對殘差序列進行分析處理,采用單一閾值或自適應(yīng)閾值超限報警方法往往存在時間不確定度問題。為了提高齒輪箱故障預(yù)警的可靠性和準(zhǔn)確度,本文通過CUSUM控制圖構(gòu)造Change-Point方法,通過Change-Point方法進行殘差序列分析。

    3.1 Change-Point算法

    Change-Point算法是利用一定的統(tǒng)計指標(biāo)或統(tǒng)計方法,對時間序列的狀態(tài)進行觀測,以便準(zhǔn)確有效的估計出變點的位置。變點問題分連續(xù)形式和離散形式兩種[13],對于風(fēng)電機組運行數(shù)據(jù)來說是時間上的連續(xù),如若對樣本所有觀察點進行變點提取,應(yīng)采用連續(xù)形式。

    自20世紀(jì)70年代以來,許多統(tǒng)計學(xué)家投入到變點問題研究領(lǐng)域,估計和檢測變點問題的方法也不斷發(fā)展完善[14-15],如最小二乘法、極大似然法、累積和法(Cumulative Sum, CUSUM)。CUSUM法是通過對觀測值與目標(biāo)值之差的累積和來描點,因其方法簡潔有效,在經(jīng)濟學(xué)等領(lǐng)域得到了廣泛的應(yīng)用。假設(shè)在一段時間內(nèi)統(tǒng)計值高于全體均值,高于均值的數(shù)據(jù)量會持續(xù)累積,CUSUM控制圖呈現(xiàn)穩(wěn)步增長趨勢;相反,假設(shè)在一段時間內(nèi)統(tǒng)計值低于全體均值,低于均值的數(shù)據(jù)量同樣會持續(xù)累積,CUSUM控制圖呈現(xiàn)穩(wěn)步下降趨勢。某些數(shù)據(jù)特性在不能直觀地發(fā)現(xiàn)顯著變化,通過原始數(shù)據(jù)無法準(zhǔn)確定位到變點時間域,而通過CUSUM控制圖累加數(shù)據(jù)特性后可總體定位到變點位置,本文選用CUSUM法實現(xiàn)Change-Point算法。

    CUSUM的主要原理為假設(shè)數(shù)據(jù)模型:

    xt=u(t)+εt,t=1,2,3,…,n

    (12)

    式中:u(t)是一個非隨機函數(shù);εt為線性過程??紤]比較簡單的情形,即u(t)只取兩個值:

    (13)

    式中:u1,u2和k*為未知量。變點k*的CUSUM估計定義為:

    (14)

    其中:

    0≤γ≤1

    (15)

    圖5 Change-Point變點搜尋流程圖

    使用CUSUM方法確定變點的總體位置,在搜尋變點位置時存在時間不確定度問題。為了解決不確定度問題,算法在實際應(yīng)用中使用了改進的CUSUM方法[16]。本文使用二分法進行分割搜索,通過CUSUM方法找到變點總體位置后進行二分法分割處理,通過計算置信度的大小定位變點發(fā)生區(qū)間,繼續(xù)使用二分法分割,直到準(zhǔn)確搜尋到變點時刻。以此方法保證定位時刻精度,算法具體流程如圖5所示。本文使用Change-Point方法對XGBoost模型預(yù)測殘差序列進行變點分析,對變點進行原因回溯,以此分析狀態(tài)變化點深層次誘因。

    3.2 基于Change-Point的殘差序列分析

    風(fēng)電機組運行過程中,某些原因?qū)е慢X輪箱故障,故障發(fā)生在一個持續(xù)的時間序列中。通過上述方法構(gòu)建的Change-Point方法對上述模型預(yù)測殘差序列進行變點分析。Change-Point每一個數(shù)據(jù)窗口長度定義為10分鐘,將數(shù)據(jù)窗口精度設(shè)置為12,即找到變點的誤差為兩小時。引導(dǎo)分析限值設(shè)置為1 000,引導(dǎo)分析限值為每一次分析的間隔尺度,如第一次分析第1個數(shù)據(jù)點,則下一次分析第1 001個數(shù)據(jù)點。置信度限值設(shè)置為0.99,通過計算找到變點4個,其殘差序列變點分析圖如圖6所示。

    圖6 殘差序列變點分析圖

    利用CUSUM控制圖統(tǒng)計殘差序列的變化趨勢,殘差CUSUM控制圖如圖7所示,并對變點進行編號。

    圖7 殘差CUSUM控制圖

    根據(jù)風(fēng)電機組運行日志,圖7中Change-Point1為對齒輪箱散熱片進行定期清理的時間點。Change-Point2為6月10號22:35:00,Change-Point4為7月6號1:10:00,分別對應(yīng)圖5中傳感器接線松動和緊固傳感器接線時間點。Change-Point3是齒輪箱軸溫異常超限報警時間點,對應(yīng)6月27號23:40:00。該方法可準(zhǔn)確定位到狀態(tài)變化時間點,且其CUSUM累積特性可規(guī)避因個別孤立的殘差較大的點導(dǎo)致的超限誤報現(xiàn)象,避免隨機因素的影響。核對該風(fēng)電機組運行數(shù)據(jù),在Change-Point2與Change-Point3時間段內(nèi)已出現(xiàn)數(shù)據(jù)特性變化異常,只是未觸發(fā)閾值限而引發(fā)報警。相對于機組控制系統(tǒng)閾值超溫報警,該方法可提前17天發(fā)現(xiàn)齒輪箱軸溫異常并發(fā)出報警,證明了該方法的有效性。

    4 結(jié)論

    (1)本文基于SCADA運行數(shù)據(jù),采用相關(guān)性分析法和風(fēng)電機組控制機理提取與齒輪箱軸溫密切相關(guān)聯(lián)的變量,并增加部分分類變量共同作為模型輸入。利用XGBoost算法可建立齒輪箱正常工作狀態(tài)溫度模型,進行齒輪箱溫度預(yù)警。

    (2)采用5種模型進行齒輪箱溫度預(yù)測模型實驗,最終得出將XGBoost算法應(yīng)用于齒輪箱溫度預(yù)測要優(yōu)于其它4類模型,具有預(yù)測精度和計算效率高,資源占用低的特點。該模型可提前發(fā)現(xiàn)齒輪箱軸溫異常。

    (3)構(gòu)建了新的Change-Point算法,通過CUSUM控制圖實現(xiàn)Change-Point算法,利用Change-Point技術(shù)進行XGBoost模型預(yù)測殘差序列分析并給出變點。在齒輪箱軸溫故障實例中,該方法可準(zhǔn)確定位狀態(tài)劣化時間點,回溯狀態(tài)劣化誘因。

    猜你喜歡
    軸溫變點齒輪箱
    某型車轉(zhuǎn)向架軸溫監(jiān)測裝置研究
    風(fēng)電齒輪箱軸承用鋼100CrMnSi6-4的開發(fā)
    山東冶金(2022年3期)2022-07-19 03:24:36
    回歸模型參數(shù)的變點檢測方法研究
    基于城際動車組軸溫實時檢測系統(tǒng)的優(yōu)化設(shè)計
    廣東科技(2021年2期)2021-03-06 09:46:56
    正態(tài)分布序列均值變點檢測的貝葉斯方法
    CR400BF 型動車組達(dá)速交路軸溫變化規(guī)律研究
    基于二元分割的多變點估計
    獨立二項分布序列變點的識別方法
    基于光纖光柵溫度傳感的車輛軸溫監(jiān)測系統(tǒng)
    提高齒輪箱式換檔機構(gòu)可靠性的改進設(shè)計
    春色校园在线视频观看| 尤物成人国产欧美一区二区三区| 亚洲18禁久久av| 亚洲精品一区av在线观看| 国产视频一区二区在线看| 99久久精品一区二区三区| 精品久久久久久成人av| 能在线免费观看的黄片| 日韩高清综合在线| 亚洲va在线va天堂va国产| 日韩av在线大香蕉| 久久久久久久久久黄片| 亚洲成人中文字幕在线播放| 国内少妇人妻偷人精品xxx网站| 国产精品一区二区三区四区免费观看 | 人妻夜夜爽99麻豆av| 亚洲人成网站在线播放欧美日韩| 一级黄色大片毛片| 人妻制服诱惑在线中文字幕| 波多野结衣巨乳人妻| 亚洲欧美日韩高清在线视频| 久久久国产成人精品二区| 三级毛片av免费| 真人一进一出gif抽搐免费| 少妇猛男粗大的猛烈进出视频 | 亚洲国产色片| 欧美在线一区亚洲| 亚洲性夜色夜夜综合| 免费高清视频大片| 一进一出好大好爽视频| 亚洲精品成人久久久久久| 婷婷精品国产亚洲av在线| 欧美三级亚洲精品| 一个人看的www免费观看视频| 波野结衣二区三区在线| 欧美精品啪啪一区二区三区| 国产伦在线观看视频一区| 国产精品野战在线观看| 人妻丰满熟妇av一区二区三区| 亚洲国产精品久久男人天堂| 欧美高清性xxxxhd video| 夜夜夜夜夜久久久久| 国产中年淑女户外野战色| 免费观看的影片在线观看| 欧美成人a在线观看| 亚洲av第一区精品v没综合| 女人十人毛片免费观看3o分钟| 免费大片18禁| 国产精品日韩av在线免费观看| 成人三级黄色视频| 精品一区二区三区视频在线观看免费| 中文亚洲av片在线观看爽| 精品一区二区三区人妻视频| 成年女人看的毛片在线观看| 熟女电影av网| 欧美日韩黄片免| 国产精品一区二区性色av| 免费在线观看影片大全网站| 国产高潮美女av| 精品人妻1区二区| 欧美又色又爽又黄视频| 亚洲精品一区av在线观看| 最好的美女福利视频网| 亚洲国产精品sss在线观看| 自拍偷自拍亚洲精品老妇| 精品人妻视频免费看| 午夜福利18| 日本与韩国留学比较| 久久久午夜欧美精品| 精品人妻偷拍中文字幕| 村上凉子中文字幕在线| 午夜福利18| 国产三级中文精品| 亚洲精品亚洲一区二区| 亚洲成人精品中文字幕电影| 日本色播在线视频| 国内揄拍国产精品人妻在线| 国产三级在线视频| 国产精品福利在线免费观看| 国产精品福利在线免费观看| 亚洲欧美日韩无卡精品| 美女大奶头视频| 精品国产三级普通话版| 亚洲精品色激情综合| 欧美黑人欧美精品刺激| 欧美精品国产亚洲| 免费高清视频大片| 男女边吃奶边做爰视频| 嫩草影院精品99| 99国产精品一区二区蜜桃av| 欧美在线一区亚洲| 18禁裸乳无遮挡免费网站照片| 精品久久久久久久末码| 久久人妻av系列| 欧美日韩亚洲国产一区二区在线观看| 91麻豆av在线| 天堂影院成人在线观看| 人妻丰满熟妇av一区二区三区| 午夜精品一区二区三区免费看| 99热这里只有是精品50| 欧美最黄视频在线播放免费| 老司机深夜福利视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 韩国av一区二区三区四区| 看黄色毛片网站| 伦精品一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 日韩精品中文字幕看吧| 亚洲色图av天堂| 日本成人三级电影网站| av天堂在线播放| 久久人人精品亚洲av| 亚洲av中文av极速乱 | 一个人免费在线观看电影| 深爱激情五月婷婷| 欧美性猛交黑人性爽| 五月玫瑰六月丁香| 在线免费观看的www视频| 黄色女人牲交| 久久精品国产亚洲av香蕉五月| 少妇丰满av| 精华霜和精华液先用哪个| 国产午夜精品久久久久久一区二区三区 | 69人妻影院| 欧美日韩乱码在线| 欧美日韩乱码在线| 国产黄色小视频在线观看| 日韩av在线大香蕉| 俺也久久电影网| 成年人黄色毛片网站| xxxwww97欧美| 日本在线视频免费播放| 啦啦啦韩国在线观看视频| 亚洲成av人片在线播放无| 看片在线看免费视频| 天堂网av新在线| 真人做人爱边吃奶动态| 欧美激情国产日韩精品一区| 99久久九九国产精品国产免费| 国产白丝娇喘喷水9色精品| 成人高潮视频无遮挡免费网站| netflix在线观看网站| 国产精品一区二区三区四区免费观看 | 可以在线观看的亚洲视频| 国产亚洲欧美98| 黄色丝袜av网址大全| 国产av麻豆久久久久久久| 搡老熟女国产l中国老女人| 女同久久另类99精品国产91| 九九在线视频观看精品| 国产极品精品免费视频能看的| 亚洲成人精品中文字幕电影| 草草在线视频免费看| 波多野结衣高清作品| 我要搜黄色片| 久久国内精品自在自线图片| 久久久久性生活片| 色av中文字幕| 日日夜夜操网爽| 中文资源天堂在线| 熟妇人妻久久中文字幕3abv| 看免费成人av毛片| 亚洲国产高清在线一区二区三| 亚洲一区二区三区色噜噜| 久久热精品热| 亚洲av中文av极速乱 | 亚洲av免费高清在线观看| 99久久精品国产国产毛片| 午夜精品久久久久久毛片777| 热99re8久久精品国产| 99热这里只有精品一区| 日韩欧美国产在线观看| 国产色婷婷99| 欧美高清成人免费视频www| 国产成人一区二区在线| 少妇熟女aⅴ在线视频| 波多野结衣高清无吗| 久久这里只有精品中国| 久久这里只有精品中国| 他把我摸到了高潮在线观看| 亚洲av成人av| 91久久精品国产一区二区三区| 99九九线精品视频在线观看视频| 99久久无色码亚洲精品果冻| 男插女下体视频免费在线播放| av女优亚洲男人天堂| 成人鲁丝片一二三区免费| 直男gayav资源| a级一级毛片免费在线观看| 国产午夜精品论理片| 看免费成人av毛片| 亚洲国产高清在线一区二区三| 91午夜精品亚洲一区二区三区 | 亚洲第一电影网av| 人妻少妇偷人精品九色| 精品欧美国产一区二区三| 在线看三级毛片| 九九爱精品视频在线观看| 黄色欧美视频在线观看| 久久久久久久久久成人| 国产精品福利在线免费观看| 有码 亚洲区| 亚洲国产精品成人综合色| 午夜爱爱视频在线播放| 99热这里只有精品一区| 亚洲av中文字字幕乱码综合| 婷婷亚洲欧美| 亚洲成人久久性| 少妇丰满av| 欧美一区二区国产精品久久精品| 精品免费久久久久久久清纯| 禁无遮挡网站| 91精品国产九色| 十八禁国产超污无遮挡网站| 久久久久久久久中文| 国产毛片a区久久久久| 男女啪啪激烈高潮av片| 国产免费av片在线观看野外av| 国产精品爽爽va在线观看网站| 在线观看午夜福利视频| 一级av片app| 日日摸夜夜添夜夜添小说| 中文亚洲av片在线观看爽| 日本色播在线视频| 国产精品日韩av在线免费观看| 亚洲五月天丁香| 国产人妻一区二区三区在| 特大巨黑吊av在线直播| 在线看三级毛片| 国产精品一区二区免费欧美| 亚洲人成网站在线播| 精品久久久久久久久久免费视频| 国产伦在线观看视频一区| 日本黄色视频三级网站网址| av国产免费在线观看| 成人特级av手机在线观看| 亚洲精华国产精华精| 日本在线视频免费播放| 欧美人与善性xxx| 国产老妇女一区| 日日啪夜夜撸| 国产高清视频在线观看网站| 亚洲欧美日韩卡通动漫| 国内精品久久久久久久电影| 国产av一区在线观看免费| 中国美白少妇内射xxxbb| 一级黄色大片毛片| 午夜福利在线在线| 91麻豆精品激情在线观看国产| 国产精品日韩av在线免费观看| 免费搜索国产男女视频| 成人性生交大片免费视频hd| 免费黄网站久久成人精品| 99久久精品一区二区三区| 亚洲色图av天堂| 国产乱人视频| 香蕉av资源在线| 18禁在线播放成人免费| 亚洲欧美清纯卡通| 99久国产av精品| 美女被艹到高潮喷水动态| 婷婷精品国产亚洲av| 男女啪啪激烈高潮av片| 黄色女人牲交| 18禁在线播放成人免费| 99精品久久久久人妻精品| 九九爱精品视频在线观看| 国产伦精品一区二区三区四那| 人妻制服诱惑在线中文字幕| 亚洲精品日韩av片在线观看| 亚洲无线在线观看| 国产精品一区www在线观看 | 99久久久亚洲精品蜜臀av| 97超视频在线观看视频| 免费搜索国产男女视频| 国产精品伦人一区二区| 韩国av一区二区三区四区| 国产三级在线视频| 亚洲av一区综合| 久9热在线精品视频| 国产高清三级在线| 国产乱人伦免费视频| 2021天堂中文幕一二区在线观| 天堂√8在线中文| 成人无遮挡网站| 亚洲av不卡在线观看| 亚洲,欧美,日韩| 国产精品久久久久久亚洲av鲁大| 色吧在线观看| 亚洲国产精品成人综合色| 麻豆av噜噜一区二区三区| 麻豆一二三区av精品| 国产精品久久久久久亚洲av鲁大| 亚洲熟妇中文字幕五十中出| 成人av在线播放网站| 老司机午夜福利在线观看视频| 少妇的逼水好多| 欧美日韩黄片免| 丰满的人妻完整版| 日日干狠狠操夜夜爽| 亚洲av不卡在线观看| 老熟妇乱子伦视频在线观看| 九九久久精品国产亚洲av麻豆| 狠狠狠狠99中文字幕| 欧美中文日本在线观看视频| 尾随美女入室| 成人三级黄色视频| 国产伦人伦偷精品视频| 麻豆av噜噜一区二区三区| 亚洲中文日韩欧美视频| 12—13女人毛片做爰片一| a级毛片免费高清观看在线播放| 亚洲欧美日韩高清专用| 欧美又色又爽又黄视频| 亚洲自偷自拍三级| 亚洲四区av| 亚洲欧美日韩高清在线视频| 亚洲中文字幕一区二区三区有码在线看| 有码 亚洲区| xxxwww97欧美| 成人永久免费在线观看视频| av中文乱码字幕在线| 国产黄a三级三级三级人| 日本欧美国产在线视频| 国产又黄又爽又无遮挡在线| h日本视频在线播放| www日本黄色视频网| 一级毛片久久久久久久久女| 无人区码免费观看不卡| 免费一级毛片在线播放高清视频| av.在线天堂| 黄片wwwwww| 最后的刺客免费高清国语| 国产成人a区在线观看| 永久网站在线| 国产国拍精品亚洲av在线观看| 成人av在线播放网站| 亚洲国产日韩欧美精品在线观看| 亚洲乱码一区二区免费版| 偷拍熟女少妇极品色| 久久草成人影院| 在线观看66精品国产| 国产成人福利小说| 亚洲国产色片| 国内精品久久久久精免费| 悠悠久久av| 国产一区二区亚洲精品在线观看| 午夜激情欧美在线| 色5月婷婷丁香| 午夜福利18| 少妇裸体淫交视频免费看高清| 在线免费观看不下载黄p国产 | 国产三级在线视频| 免费看美女性在线毛片视频| 两人在一起打扑克的视频| 国产伦精品一区二区三区四那| ponron亚洲| 亚洲av日韩精品久久久久久密| 精品无人区乱码1区二区| 日韩在线高清观看一区二区三区 | 亚洲综合色惰| 男插女下体视频免费在线播放| 黄色欧美视频在线观看| 午夜福利高清视频| 亚洲性夜色夜夜综合| 国产精品久久电影中文字幕| av天堂中文字幕网| 波多野结衣高清作品| 麻豆精品久久久久久蜜桃| 日韩,欧美,国产一区二区三区 | 精品一区二区免费观看| 性色avwww在线观看| 99热精品在线国产| 日本免费一区二区三区高清不卡| or卡值多少钱| 99久国产av精品| .国产精品久久| 国产精品,欧美在线| 人妻少妇偷人精品九色| 搡老熟女国产l中国老女人| 精品国产三级普通话版| 日韩亚洲欧美综合| 欧美日本视频| 久久这里只有精品中国| 久久久色成人| 真人做人爱边吃奶动态| 少妇被粗大猛烈的视频| 俺也久久电影网| 小蜜桃在线观看免费完整版高清| 国产精品福利在线免费观看| 欧美黑人巨大hd| 亚洲人成网站高清观看| 人人妻,人人澡人人爽秒播| 一级毛片久久久久久久久女| 欧美极品一区二区三区四区| 国产视频一区二区在线看| 长腿黑丝高跟| 欧美日韩亚洲国产一区二区在线观看| 又紧又爽又黄一区二区| 女生性感内裤真人,穿戴方法视频| 美女高潮的动态| 国产毛片a区久久久久| 99热只有精品国产| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区三区四区免费观看 | 天堂网av新在线| 色在线成人网| 俺也久久电影网| 欧美高清成人免费视频www| 久久精品夜夜夜夜夜久久蜜豆| 一本一本综合久久| 国内精品一区二区在线观看| 国产91精品成人一区二区三区| 男人和女人高潮做爰伦理| 成人精品一区二区免费| 搡老岳熟女国产| 午夜福利视频1000在线观看| 日日摸夜夜添夜夜添小说| 日本色播在线视频| 亚洲精品久久国产高清桃花| 欧美最黄视频在线播放免费| 99国产极品粉嫩在线观看| 日日干狠狠操夜夜爽| 日韩大尺度精品在线看网址| 99久久精品热视频| 日本黄色片子视频| 日韩中文字幕欧美一区二区| 亚洲无线在线观看| 人妻久久中文字幕网| 在线天堂最新版资源| 琪琪午夜伦伦电影理论片6080| 女人被狂操c到高潮| 日本成人三级电影网站| 国产一区二区激情短视频| 一级av片app| 国产熟女欧美一区二区| 哪里可以看免费的av片| 美女高潮的动态| 亚洲专区国产一区二区| 国产成人福利小说| 极品教师在线视频| av.在线天堂| 免费看美女性在线毛片视频| 色尼玛亚洲综合影院| 免费在线观看影片大全网站| 丝袜美腿在线中文| 色综合婷婷激情| 欧美bdsm另类| 亚洲四区av| 九九热线精品视视频播放| 级片在线观看| 国产高清视频在线播放一区| 欧美激情在线99| 99久久精品一区二区三区| 久久亚洲精品不卡| 午夜激情福利司机影院| 亚洲一区二区三区色噜噜| 国产高清有码在线观看视频| 可以在线观看毛片的网站| 久久精品夜夜夜夜夜久久蜜豆| 99热6这里只有精品| 久久久久精品国产欧美久久久| 色综合色国产| 亚洲精品粉嫩美女一区| 蜜桃亚洲精品一区二区三区| 日本黄色视频三级网站网址| 露出奶头的视频| 一本精品99久久精品77| 亚洲精品一区av在线观看| 一本久久中文字幕| 国产免费男女视频| 久久人人爽人人爽人人片va| 精品99又大又爽又粗少妇毛片 | 超碰av人人做人人爽久久| 国产又黄又爽又无遮挡在线| 久久人人精品亚洲av| 1024手机看黄色片| 又紧又爽又黄一区二区| 美女xxoo啪啪120秒动态图| 久久久久久久亚洲中文字幕| 国产精品av视频在线免费观看| 欧美xxxx性猛交bbbb| .国产精品久久| 亚洲av一区综合| 国产探花在线观看一区二区| 一区福利在线观看| 中文字幕精品亚洲无线码一区| 日韩中文字幕欧美一区二区| 天美传媒精品一区二区| 可以在线观看毛片的网站| 他把我摸到了高潮在线观看| 成年免费大片在线观看| 国内精品美女久久久久久| 久久亚洲真实| 国产精品亚洲美女久久久| 精品午夜福利在线看| 99热这里只有精品一区| av中文乱码字幕在线| 麻豆成人午夜福利视频| 变态另类成人亚洲欧美熟女| 亚洲精品粉嫩美女一区| 午夜精品一区二区三区免费看| 国产免费av片在线观看野外av| 亚洲成a人片在线一区二区| 国产一级毛片七仙女欲春2| 91久久精品国产一区二区成人| 18+在线观看网站| 波多野结衣高清作品| 国内久久婷婷六月综合欲色啪| 精品久久久久久,| 18禁黄网站禁片午夜丰满| 在线天堂最新版资源| 亚洲aⅴ乱码一区二区在线播放| 欧美日本视频| 淫秽高清视频在线观看| 午夜精品一区二区三区免费看| 男女做爰动态图高潮gif福利片| 国产69精品久久久久777片| 久久久久久九九精品二区国产| 亚洲精品久久国产高清桃花| 免费看光身美女| 嫩草影院精品99| 联通29元200g的流量卡| 高清日韩中文字幕在线| 高清在线国产一区| 亚洲综合色惰| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看| x7x7x7水蜜桃| 老熟妇仑乱视频hdxx| 午夜福利视频1000在线观看| 变态另类成人亚洲欧美熟女| 日日干狠狠操夜夜爽| av在线亚洲专区| 综合色av麻豆| 啦啦啦观看免费观看视频高清| 亚洲一级一片aⅴ在线观看| 欧洲精品卡2卡3卡4卡5卡区| 美女大奶头视频| 夜夜夜夜夜久久久久| 在线观看av片永久免费下载| 久久午夜亚洲精品久久| 欧美日韩综合久久久久久 | 亚洲成a人片在线一区二区| 国产爱豆传媒在线观看| 真人一进一出gif抽搐免费| 成年女人看的毛片在线观看| 国模一区二区三区四区视频| 免费无遮挡裸体视频| 亚洲av熟女| h日本视频在线播放| 免费人成视频x8x8入口观看| 欧美不卡视频在线免费观看| 18禁在线播放成人免费| 国产探花极品一区二区| 亚洲七黄色美女视频| 不卡视频在线观看欧美| 天美传媒精品一区二区| 成人高潮视频无遮挡免费网站| 日韩av在线大香蕉| 日韩精品中文字幕看吧| 女人十人毛片免费观看3o分钟| 老熟妇乱子伦视频在线观看| 色精品久久人妻99蜜桃| 变态另类成人亚洲欧美熟女| 最后的刺客免费高清国语| 熟女电影av网| 久久国内精品自在自线图片| 中文亚洲av片在线观看爽| 午夜免费激情av| 全区人妻精品视频| 嫁个100分男人电影在线观看| 欧美日韩综合久久久久久 | 欧美极品一区二区三区四区| www日本黄色视频网| 欧美日韩亚洲国产一区二区在线观看| 男女做爰动态图高潮gif福利片| 99热这里只有精品一区| 男女做爰动态图高潮gif福利片| 高清在线国产一区| 最新中文字幕久久久久| 内地一区二区视频在线| 99久久精品热视频| 色尼玛亚洲综合影院| 精品一区二区免费观看| 老师上课跳d突然被开到最大视频| 亚洲av.av天堂| 精品久久久久久久久久免费视频| 免费高清视频大片| 免费人成在线观看视频色| 国国产精品蜜臀av免费| 国产精品不卡视频一区二区| 精品免费久久久久久久清纯| 国产毛片a区久久久久| 色播亚洲综合网| 国产高清不卡午夜福利| 美女cb高潮喷水在线观看| 99久久久亚洲精品蜜臀av| 亚洲人成伊人成综合网2020| av中文乱码字幕在线| 成年女人毛片免费观看观看9| 色噜噜av男人的天堂激情| 亚洲图色成人| 国产精品一及| 日韩大尺度精品在线看网址| 色哟哟·www| 在线播放无遮挡| 久久6这里有精品| 女同久久另类99精品国产91| 99精品久久久久人妻精品| 色在线成人网| 伦理电影大哥的女人| 亚洲va在线va天堂va国产| 亚洲第一区二区三区不卡| 国产成人aa在线观看|