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

    基于不同學習算法的RBF神經(jīng)網(wǎng)絡(luò)在故障診斷中的應(yīng)用

    2016-02-05 07:30:58胡克磊
    上海電力大學學報 2016年6期
    關(guān)鍵詞:征兆凝汽器訓練方法

    李 玲, 胡克磊

    (1.上海電力學院 自動化工程學院, 上海 200090; 2.安徽華電六安電廠有限公司 熱控分廠, 安徽 六安 237000)

    基于不同學習算法的RBF神經(jīng)網(wǎng)絡(luò)在故障診斷中的應(yīng)用

    李 玲1, 胡克磊2

    (1.上海電力學院 自動化工程學院, 上海 200090; 2.安徽華電六安電廠有限公司 熱控分廠, 安徽 六安 237000)

    基于對RBF神經(jīng)網(wǎng)絡(luò)常用的3種學習算法的研究,通過對凝汽器典型故障類型與故障征兆分析,提出了基于不同學習算法的RBF神經(jīng)網(wǎng)絡(luò)凝汽器故障診斷,并對診斷結(jié)果進行比較.診斷結(jié)果表明,基于3種常見學習算法的RBF神經(jīng)網(wǎng)絡(luò)都可以準確診斷出凝汽器的各種故障,但聚類方法和OLS算法學習速度要快得多,梯度訓練方法速度較慢.研究還表明,RBF神經(jīng)網(wǎng)絡(luò)在故障診斷領(lǐng)域具有很好的實用性.

    RBF神經(jīng)網(wǎng)絡(luò); 學習算法; 凝汽器; 故障診斷

    作為現(xiàn)代大型電站汽輪機組的主要輔助設(shè)備,凝汽器設(shè)備在汽輪機組的熱力循環(huán)中起著冷源的作用,其工作性能直接影響到整個汽輪機組的安全性和經(jīng)濟性,因此,凝汽器設(shè)備狀態(tài)的監(jiān)測和故障診斷一直是研究熱點[1-2].凝汽器設(shè)備發(fā)生故障的原因與其故障征兆之間是非線性關(guān)系,具有隨機性、模糊性、復雜性等,精確的數(shù)學公式很難將其表達出來.而人工神經(jīng)網(wǎng)絡(luò)的特點是自學習、自組織能力突出,具有很強非線性模式分類能力,非常適合對設(shè)備進行故障診斷.

    目前,反向傳播(Back Propagation,BP)神經(jīng)網(wǎng)絡(luò)在故障診斷領(lǐng)域中應(yīng)用最為廣泛,但其在故障診斷中具有容易陷入局部最小、訓練收斂速度慢等缺點[3].而徑向基(Radical Basis Function,RBF)神經(jīng)網(wǎng)絡(luò)具有收斂速度快、不存在局部極小值、能夠以任意精度逼近任意非線性函數(shù)、泛化能力強等特點,目前它是最引起研究者關(guān)注的神經(jīng)網(wǎng)絡(luò)模型之一[4].

    本文將RBF神經(jīng)網(wǎng)絡(luò)應(yīng)用于凝汽器故障診斷,同時基于對RBF神經(jīng)網(wǎng)絡(luò)常見學習算法的學習和研究,將不同學習算法的RBF神經(jīng)網(wǎng)絡(luò)在凝汽器故障診斷中的應(yīng)用進行了比較.

    1 凝汽器故障類型和征兆集

    基于神經(jīng)網(wǎng)絡(luò)故障診斷的基本思想就是:第一,利用一組故障樣本訓練網(wǎng)絡(luò),確定網(wǎng)絡(luò)結(jié)構(gòu)和參數(shù);第二,進行故障診斷,即向網(wǎng)絡(luò)輸入一組征兆,實現(xiàn)征兆集到故障集之間的非線性映射[5].凝汽器系統(tǒng)結(jié)構(gòu)非常復雜且運行環(huán)境特殊,發(fā)生故障時具有多種故障原因和故障征兆.本文確定了凝汽器的14種典型故障類型,通過分析這些故障發(fā)生的原因,選取了16種故障征兆,分別如表1和表2所示.對這些故障征兆進行規(guī)一化處理并根據(jù)系統(tǒng)運行參數(shù)和結(jié)構(gòu)參數(shù),最終使征兆取值范圍為[0,1][6-7].結(jié)合文獻[6]得到凝汽器系統(tǒng)的故障樣本征兆數(shù)據(jù)表,如表3所示.

    表1 凝汽器故障類型集合

    表2 凝汽器故障征兆集合

    在RBF神經(jīng)網(wǎng)絡(luò)訓練中,訓練樣本集包含征兆集和故障集.征兆集作為樣本輸入,即表3中數(shù)據(jù);故障集作為樣本輸出,如表4所示,每個輸出分別取兩個狀態(tài),0 代表無故障狀態(tài),1代表故障狀態(tài).RBF神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)可以確定為:輸入層n=16,每一個故障征兆參數(shù)對應(yīng)一個神經(jīng)元的輸入;隱含層h由RBF神經(jīng)網(wǎng)絡(luò)的學習算法確定;輸出層m=14,每一個故障類型對應(yīng)一個神經(jīng)元輸出.

    表3 凝汽器系統(tǒng)故障征兆數(shù)據(jù)

    表4 故障類型對應(yīng)的網(wǎng)絡(luò)輸出

    2 不同學習算法及故障診斷應(yīng)用

    RBF神經(jīng)網(wǎng)絡(luò)學習算法的選擇十分重要,網(wǎng)絡(luò)的隱層節(jié)點數(shù)h,數(shù)據(jù)中心ci,擴展常數(shù)σi以及輸出權(quán)值修正都由學習算法在網(wǎng)絡(luò)訓練時確定.RBF神經(jīng)網(wǎng)絡(luò)常用的學習算法是聚類方法(最常用的是k-means聚類),梯度訓練方法,正交最小二乘(orthogonalleastsquares,OLS)算法.這3種學習算法各有優(yōu)缺點,聚類方法可以根據(jù)各聚類中心之間的距離得出隱節(jié)點的寬度,但是確定數(shù)據(jù)中心時沒有用到樣本輸出信息,只用到了樣本輸入信息;梯度訓練方法具有更新隱含層數(shù)據(jù)中心及寬度的能力,顯著地提高了RBF網(wǎng)絡(luò)的性能,但其增加了訓練算法的復雜度和網(wǎng)絡(luò)訓練時間;正交最小二乘算法能夠保證學習誤差不大于給定值,且能在權(quán)值學習的同時確定隱節(jié)點的數(shù)目,較容易實現(xiàn),但該算法并不一定能設(shè)計出具有最小結(jié)構(gòu)的RBF神經(jīng)網(wǎng)絡(luò)[8-10].

    由于輸出層的神經(jīng)元之間沒有相互作用,基于這個原因,不失一般性,可以考慮單輸出的RBF神經(jīng)網(wǎng)絡(luò).有多于一個輸出的網(wǎng)絡(luò)可以看作幾個單輸出網(wǎng)絡(luò)共享一個隱含層的疊加.在以下RBF神經(jīng)網(wǎng)絡(luò)學習算法中,x1,x2,x3,…,xN為樣本輸入,相應(yīng)的樣本輸出為y1,y2,y3,…,yN,其中N為樣本個數(shù).

    2.1 基于聚類方法的RBF神經(jīng)網(wǎng)絡(luò)凝汽器故障診斷

    2.1.1 聚類方法

    在RBF神經(jīng)網(wǎng)絡(luò)學習算法中,聚類方法是最經(jīng)典的,該算法的實現(xiàn)過程如下[11-13].

    (1) 設(shè)k為迭代次數(shù),此時的聚類中心為c1(k),c2(k),c3(k),…,ch(k),相應(yīng)的聚類域為θ1(k),θ2(k),θ3(k),…,θh(k).初始化算法,令k=1,選擇h個不同的初始聚類中心.

    (2) 計算輸入樣本與數(shù)據(jù)中心的距離,即‖xj-ci(k)‖,i=1,2,3,…,h;j=1,2,3,…,N.

    (3) 對輸入xj按最小距離原則進行分類.當i(xj)=mini‖xj-ci(k)‖時,xj被歸為第i類,即xj∈θi(k).

    (4) 重新計算聚類中心,即:

    (1)

    式中:x——輸入樣本;Ni——聚類域;θi(k)——所含有的樣本數(shù).

    (5) 如果ci(k+1)≠ci(k),則轉(zhuǎn)到步驟(2),否則進行下一步驟.

    (6) 計算各隱節(jié)點的寬度.隱節(jié)點的寬度為:

    (2)

    式中:α——重疊系數(shù);di=mini‖cj-ci(k)‖.

    (7) 數(shù)據(jù)中心和寬度確定后,輸出權(quán)值矢量ω=[ω1,ω2,ω3,…,ωh]T用最小二乘方法計算得到.當輸入為xj時,第i個隱節(jié)點的輸出為:

    (3)

    隱層輸出矩陣為:

    (4)

    (5)

    (6)

    2.1.2 凝汽器故障診斷

    采用表3中的故障類型和相關(guān)的征兆作為凝汽器故障診斷的訓練樣本,根據(jù)上述聚類算法的實現(xiàn)步驟,利用MATLAB編程設(shè)計訓練出基于聚類算法的RBF神經(jīng)網(wǎng)絡(luò),實現(xiàn)對這些故障的識別能力.其中,隱節(jié)點重疊系數(shù)α取1;該RBF神經(jīng)網(wǎng)絡(luò)的隱含層由14個神經(jīng)元組成,這是經(jīng)過多次訓練得出的.本文進行凝汽器故障診斷所用的故障樣本如表5所示,即凝汽器故障數(shù)據(jù)集.其中,數(shù)據(jù)1故障征兆對應(yīng)的故障是Y13,即循環(huán)水量不足;數(shù)據(jù)2故障征兆對應(yīng)的故障是Y9,即冷卻管堵塞;數(shù)據(jù)3故障征兆對應(yīng)的故障是Y4,即真空系統(tǒng)管路破裂.

    利用訓練好的基于聚類算法的RBF神經(jīng)網(wǎng)絡(luò)對表5數(shù)據(jù)進行凝汽器故障診斷,診斷結(jié)果如表6所示.

    表5 故障數(shù)據(jù)集

    其中,網(wǎng)絡(luò)輸出值第g維的取值為相應(yīng)第g種故障發(fā)生的概率(g=1,2,3,…,14).由表6可知:數(shù)據(jù)1最大故障類型值為1.000 0,遠超過其他故障類型值,因此可判斷該故障類型為Y13,即循環(huán)水量不足;數(shù)據(jù)2最大故障類型值為0.950 4,遠大于其他故障類型值,因此可判斷該故障類型為Y9,即冷卻管堵塞;數(shù)據(jù)3最大故障類型值為0.958 1,比其他故障類型值都大許多,因此可判斷該故障類型為Y4,真空系統(tǒng)管路破裂.診斷結(jié)果與實際故障情況完全相符.

    表6 基于聚類方法的RBF故障診斷結(jié)果

    由聚類算法的描述中看到,基于聚類方法的RBF網(wǎng)絡(luò)在訓練確定數(shù)據(jù)中心時沒有用到樣本輸出信息,而是只用到了樣本輸入信息;網(wǎng)絡(luò)訓練中可調(diào)參數(shù)僅為輸出層的權(quán)值,該方法導致訓練算法過于簡單.而梯度訓練方法允許調(diào)整所有的3組網(wǎng)絡(luò)參數(shù),即隱節(jié)點數(shù)據(jù)中心、徑向基函數(shù)擴展常數(shù)、輸出層權(quán)值,極大地提高了RBF神經(jīng)網(wǎng)絡(luò)的性能.

    2.2 基于梯度訓練方法的RBF神經(jīng)網(wǎng)絡(luò)凝汽器故障診斷

    2.2.1 梯度訓練方法

    RBF神經(jīng)網(wǎng)絡(luò)的梯度訓練方法通過最小化目標函數(shù)調(diào)節(jié)各隱節(jié)點聚類中心、擴展常數(shù)和輸出權(quán)值,類似于BP算法訓練多層感知器的原理.本文給出一種單輸出RBF網(wǎng)絡(luò)學習方法[14-15],其網(wǎng)絡(luò)學習的目標函數(shù)為:

    (7)

    誤差信號為:

    (8)

    梯度訓練方法具體描述如下.

    (1) 選擇隱節(jié)點數(shù)據(jù)中心,從輸入樣本集選擇數(shù)據(jù)中心.

    (2) 設(shè)定目標誤差E0,最大訓練次數(shù),適當?shù)膶W習率參數(shù)ηc,ησ,ηω,初始化RBF函數(shù)擴展常數(shù),網(wǎng)絡(luò)輸出層權(quán)值為一些小的隨機值.

    (3) 計算網(wǎng)絡(luò)輸出:

    (9)

    (4) 更新RBF神經(jīng)網(wǎng)絡(luò)參數(shù):

    (10)

    (11)

    (12)

    式中:n——訓練次數(shù).

    (5) 如果網(wǎng)絡(luò)學習的目標函數(shù)滿足E

    2.2.2 凝汽器故障診斷

    與基于聚類方法的RBF神經(jīng)網(wǎng)絡(luò)訓練時一樣,訓練時采用相同的樣本輸入,并使每種故障對應(yīng)的網(wǎng)絡(luò)輸出如表4所示.根據(jù)梯度訓練方法的原理,利用Matlab編程訓練構(gòu)建基于梯度訓練方法的RBF神經(jīng)網(wǎng)絡(luò),經(jīng)過多次調(diào)試確定該RBF神經(jīng)網(wǎng)絡(luò)隱含層由13個神經(jīng)元組成.利用基于梯度訓練方法的RBF神經(jīng)網(wǎng)絡(luò)對表4中的凝汽器故障數(shù)據(jù)進行故障診斷,結(jié)果如表7所示.

    表7 基于梯度訓練方法的RBF故障診斷結(jié)果

    由表7可知:數(shù)據(jù)1最大故障類型值為0.951 2,其余很小,最大的為0.089 6,因此可判斷該故障類型為Y13,即循環(huán)水量不足;數(shù)據(jù)2最大故障類型值為1.000 0,其余遠小于該值,因此可判斷該故障類型為Y9,即冷卻管堵塞;數(shù)據(jù)3最大故障類型值為0.820 5,其余故障類型值很小,最大的為0.174 2,因此可判斷該故障類型為Y4,即真空系統(tǒng)管路破裂.該網(wǎng)絡(luò)的診斷結(jié)果與基于聚類方法的RBF神經(jīng)網(wǎng)絡(luò)的診斷結(jié)果一致,與實際故障情況完全相符.

    在基于該算法的RBF神經(jīng)網(wǎng)絡(luò)訓練中,該算法不僅在確定數(shù)據(jù)中心時用到了樣本輸入輸出信息,同時具有更新隱含層數(shù)據(jù)中心及擴展常數(shù)的能力,這極大地提高了RBF神經(jīng)網(wǎng)絡(luò)的性能,但該學習算法增加了訓練算法的復雜度,因此也增加了訓練網(wǎng)絡(luò)需要的時間.于是本文又提出了基于正交最小二乘算法的RBF神經(jīng)網(wǎng)絡(luò)凝汽器故障診斷.

    2.3 基于正交最小二乘算法的RBF神經(jīng)網(wǎng)絡(luò)凝汽器故障診斷

    2.3.1 OLS算法

    (13)

    式中:ε——指定誤差; ω0——最優(yōu)權(quán)矢量ω的值.

    該算法選擇數(shù)據(jù)中心的步驟如下:

    (4) 采取步驟(2)的方法,選擇第二個數(shù)據(jù)中心,使用步驟(3)方法判斷算法是否終止.

    (5) 重復上面步驟,直至找到M個數(shù)據(jù)中心,使網(wǎng)絡(luò)滿足式(13).

    2.3.2 凝汽器故障診斷

    采用表3和表4中的數(shù)據(jù)作為訓練基于OLS算法的RBF神經(jīng)網(wǎng)絡(luò)的樣本和網(wǎng)絡(luò)輸出,根據(jù)OLS算法步驟,利用Matlab訓練構(gòu)建基于OLS算法的RBF神經(jīng)網(wǎng)絡(luò),經(jīng)多次訓練調(diào)試確認此時神經(jīng)網(wǎng)絡(luò)隱含層由13個神經(jīng)元構(gòu)成.利用基于OLS算法RBF神經(jīng)網(wǎng)絡(luò)對表5中凝汽器故障數(shù)據(jù)進行故障診斷,輸出如表8所示.

    表8 基于OLS算法的RBF故障診斷結(jié)果

    由表8可知:數(shù)據(jù)1最大故障類型值為0.959 0,大于0.6,另外一個較大故障類型值為0.459 7,且兩個故障類型值相差大于0.3,其余的都比后者小,故可判斷該故障類型為Y13,即循環(huán)水量不足;數(shù)據(jù)2最大故障類型值為0.952 8,其余很小,最大的為0.199 8,因此可判斷該故障類型為Y9,即冷卻管堵塞;數(shù)據(jù)3最大故障類型值為1.000 0,其余很小,最大的絕對值為0.231 9,因此可判斷該故障類型為Y4,真空系統(tǒng)管路破裂.診斷結(jié)果與前面兩種RBF神經(jīng)網(wǎng)絡(luò)的診斷結(jié)果一致,與實際故障情況完全相符.

    基于OLS算法的RBF神經(jīng)網(wǎng)絡(luò)不僅在RBF神經(jīng)網(wǎng)絡(luò)訓練確定數(shù)據(jù)中心時用到樣本輸入輸出信息,且訓練方法的復雜度比基于梯度訓練方法的RBF神經(jīng)網(wǎng)絡(luò)簡單,訓練時間短.該算法能夠自動設(shè)計滿足精度要求的網(wǎng)絡(luò)結(jié)構(gòu),有時其結(jié)構(gòu)不一定是最優(yōu),但網(wǎng)絡(luò)規(guī)模確實是相對較小的.

    3 結(jié) 語

    本文針對電廠中的重要輔助設(shè)備凝汽器故障診斷進行了研究,在分析凝汽器故障征兆和故障類型的基礎(chǔ)上,重點討論了基于3種不同算法的RBF神經(jīng)網(wǎng)絡(luò)在凝汽器故障診斷中的應(yīng)用.本文利用表4中訓練樣本集訓練得到基于3種不同算法的RBF神經(jīng)網(wǎng)絡(luò)所用的時間:基于聚類算法的RBF神經(jīng)網(wǎng)絡(luò)用時為0.005 880 s;基于梯度訓練方法的RBF神經(jīng)網(wǎng)絡(luò)用時3.755 148 s;基于OLS算法的RBF神經(jīng)網(wǎng)絡(luò)用時0.010 684 s.由于梯度訓練方法算法較其他兩種算法復雜,所以3種學習算法中,梯度訓練方法速度最慢,聚類方法和OLS方法的學習速度比較快,因此在凝汽器故障診斷中可以采用以這兩種算法.本文為相關(guān)算法在后續(xù)實際工程中的應(yīng)用提供了研究依據(jù).

    [1] 單旭昇.神經(jīng)網(wǎng)絡(luò)專家系統(tǒng)在凝汽器故障診斷中的應(yīng)用研究[D].北京:華北電力大學,2009.

    [2] 侯國蓮,孫曉剛,張建華,等.基于非線性主元分析和概率神經(jīng)網(wǎng)絡(luò)的凝汽器故障診斷方法研究[J].中國電機工程學報,2005,25(18):104-108.

    [3] 王光研,許寶杰.RBF 神經(jīng)網(wǎng)絡(luò)在旋轉(zhuǎn)機械故障診斷中的應(yīng)用[J].機械設(shè)計與制造,2008(9):57-58.

    [4] 吳偉.基于改進RBF神經(jīng)網(wǎng)絡(luò)的故障診斷技術(shù)研究[D].重慶:重慶大學,2012.

    [5] 張敬齋.RBF神經(jīng)網(wǎng)絡(luò)在水輪發(fā)電機故障診斷中的應(yīng)用[J].計算機仿真,2011,28(12):314-317;

    [6] 王晉權(quán),李磊,唐國瑞.基于RBF神經(jīng)網(wǎng)絡(luò)的凝汽器故障診斷研究[J].電力科學與工程,2007,23(4):27-31.

    [7] 熊晉魁,謝春玲,施小成,等.RBF人工神經(jīng)網(wǎng)絡(luò)在核電廠故障診斷中的應(yīng)用[J].核動力工程,2006,27(3):57-60.

    [8] 田雨波.混合神經(jīng)網(wǎng)絡(luò)[M].北京:科學出版社,2009:60-73.

    [9] 葉世偉,王海娟.神經(jīng)計算原理[M].北京:機械工業(yè)出版社,2007:101-109.

    [10] 姜鵬飛,蔡之華.基于遺傳算法和梯度下降的RBF神經(jīng)網(wǎng)絡(luò)組合訓練方法[J].計算機應(yīng)用,2007,27(2),366-372.

    [11] 張慶新,崔展博,馬睿,等.基于k-means 聚類與徑向基神經(jīng)(RBF)網(wǎng)絡(luò)的電力系統(tǒng)日負荷預(yù)測[J].科學技術(shù)與工程,2013,13(34):10 177-10 181.

    [12] 吳秋明.基于RBF神經(jīng)網(wǎng)絡(luò)的故障診斷[D].鎮(zhèn)江:江蘇大學,2006.

    [13] 張彬.基于聚類算法的RBF神經(jīng)網(wǎng)絡(luò)設(shè)計綜述[J].微型機與應(yīng)用,2012,31(12):1-3.

    [14] 魏海坤,李奇,宋文忠.梯度算法下RBF網(wǎng)的參數(shù)變化動態(tài)[J].控制理論與應(yīng)用,2007,24(3):357-360.

    [15] 吳瀟瀟,陳機林,王力.基于聚類與梯度混合學習算法RBF神經(jīng)網(wǎng)絡(luò)的電液伺服系統(tǒng)建模及仿真[J].電氣自動化,2011,33(1):20-22.

    [16] MICCHELLI C.Interpolation of scattered data:distance matrices and conditionally positive definite function [J].Constructive Approximation,1986,2(1):11-22.

    (編輯 桂金星)

    Application of Fault Diagnosis Based on RBF Neural Network withDifferent Learning Algorithms

    LI Ling1, HU Kelei2

    (1.SchoolofAutomationEngineering,ShanghaiUniversityofElectricPower,Shanghai200090,China; 2.AnhuiHuadianLuanPowerGenrationCo.,Ltd.C&I,Lu’an237000,China)

    Through the research of three kinds of learning algorithms for RBF neural network and analysis of typical faults and symptoms of the condenser,RBF neural network based on different learning algorithms for the condenser fault diagnosis is presented.Finally,the result of the diagnosis is compared.Diagnosis results show RBF neural network based on three kinds of learning algorithms can accurately diagnose various fault diagnosis of the condenser.But the study speed of the clustering method and the OLS algorithm is faster;on the contrary,the gradient training method is slower.The research also shows that the RBF neural network has good practicability in the field of fault diagnosis.

    RBF neural network; learning algorithm; condenser; fault diagnosis

    10.3969/j.issn.1006-4729.2016.06.017

    2015-10-20

    簡介:李玲(1977-),女,碩士,講師,江蘇射陽人.主要研究方向為優(yōu)化算法的研究與應(yīng)用,火電站與核電站安全及故障診斷.E-mail:liling@shiep.edu.cn.

    TP183;TK264.11

    A

    1006-4729(2016)06-0583-06

    猜你喜歡
    征兆凝汽器訓練方法
    凝汽器汽側(cè)真空泵冷卻系統(tǒng)提效改造
    凝汽器配套液環(huán)真空泵的應(yīng)用及相關(guān)設(shè)計計算
    談高中數(shù)學習題訓練方法與答題技巧
    甘肅教育(2020年18期)2020-10-28 09:07:12
    壁球反手擊球技術(shù)及其訓練方法
    跳遠運動員專項力量訓練方法
    簡論1min跳繩訓練方法
    運動(2016年7期)2016-12-01 06:34:36
    軸向排汽凝汽器膨脹節(jié)的分析研究
    電站輔機(2016年3期)2016-05-17 03:52:29
    地震征兆
    地震征兆
    熱電廠凝汽器鈦管板焊接工藝評定與應(yīng)用
    性色av乱码一区二区三区2| 美女 人体艺术 gogo| 天堂影院成人在线观看| 婷婷精品国产亚洲av| 伊人久久大香线蕉亚洲五| 国产精品久久久人人做人人爽| videosex国产| 一级作爱视频免费观看| 日日摸夜夜添夜夜添小说| 国产精品1区2区在线观看.| 国产成人系列免费观看| 19禁男女啪啪无遮挡网站| 99riav亚洲国产免费| 欧美一级毛片孕妇| 美女国产高潮福利片在线看| 亚洲成人精品中文字幕电影| 国产精品九九99| 亚洲成av人片免费观看| 免费在线观看日本一区| 久久久久久九九精品二区国产 | 亚洲一码二码三码区别大吗| 久久草成人影院| 国产亚洲精品av在线| 久久久久久人人人人人| 91成年电影在线观看| 久9热在线精品视频| 久久热在线av| 12—13女人毛片做爰片一| 免费在线观看完整版高清| 精品久久久久久久末码| 丰满人妻熟妇乱又伦精品不卡| 天天一区二区日本电影三级| 一区福利在线观看| 免费搜索国产男女视频| 欧美另类亚洲清纯唯美| 女同久久另类99精品国产91| 美女高潮喷水抽搐中文字幕| 久久性视频一级片| 欧美最黄视频在线播放免费| 大型黄色视频在线免费观看| 亚洲九九香蕉| 校园春色视频在线观看| 黑人巨大精品欧美一区二区mp4| 一a级毛片在线观看| 亚洲精品中文字幕一二三四区| 两人在一起打扑克的视频| 热re99久久国产66热| 中文字幕人妻熟女乱码| 看片在线看免费视频| ponron亚洲| 国产精品一区二区三区四区久久 | 好看av亚洲va欧美ⅴa在| 人人妻,人人澡人人爽秒播| 国产三级黄色录像| 亚洲成a人片在线一区二区| 亚洲国产欧美网| 久久精品人妻少妇| 日韩三级视频一区二区三区| 成年人黄色毛片网站| 观看免费一级毛片| 国产精品综合久久久久久久免费| 免费高清在线观看日韩| 国产日本99.免费观看| 亚洲午夜理论影院| 国产又爽黄色视频| 国产成人影院久久av| 正在播放国产对白刺激| 亚洲人成网站在线播放欧美日韩| 精品乱码久久久久久99久播| 制服诱惑二区| 9191精品国产免费久久| 啦啦啦免费观看视频1| 亚洲七黄色美女视频| 久久草成人影院| 亚洲国产精品合色在线| 黄片小视频在线播放| 嫁个100分男人电影在线观看| 欧美久久黑人一区二区| 91字幕亚洲| 久久久国产精品麻豆| 99久久国产精品久久久| 超碰成人久久| 亚洲激情在线av| 久久性视频一级片| 韩国av一区二区三区四区| 亚洲美女黄片视频| 欧美黑人巨大hd| 国产片内射在线| 亚洲欧美精品综合一区二区三区| 国产精品一区二区三区四区久久 | 国产亚洲精品久久久久久毛片| 91麻豆精品激情在线观看国产| 99riav亚洲国产免费| xxx96com| 波多野结衣高清无吗| 欧美性猛交╳xxx乱大交人| xxxwww97欧美| 白带黄色成豆腐渣| 啦啦啦 在线观看视频| 18禁黄网站禁片午夜丰满| 韩国av一区二区三区四区| 国产成人系列免费观看| 在线观看免费日韩欧美大片| 国产精品电影一区二区三区| 久久精品人妻少妇| 啪啪无遮挡十八禁网站| 看免费av毛片| 欧美av亚洲av综合av国产av| 又大又爽又粗| 波多野结衣巨乳人妻| 日韩欧美国产一区二区入口| 成人午夜高清在线视频 | 国产午夜精品久久久久久| 免费在线观看成人毛片| 欧美日韩精品网址| 老汉色av国产亚洲站长工具| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品久久男人天堂| 亚洲色图 男人天堂 中文字幕| 午夜免费激情av| 国产区一区二久久| 黄色丝袜av网址大全| 中文字幕av电影在线播放| 亚洲午夜理论影院| 在线观看66精品国产| 日韩欧美国产一区二区入口| 久久99热这里只有精品18| 久久婷婷人人爽人人干人人爱| 1024香蕉在线观看| 成在线人永久免费视频| 中国美女看黄片| 99在线视频只有这里精品首页| 国产区一区二久久| 男女那种视频在线观看| 少妇熟女aⅴ在线视频| 一区福利在线观看| 可以在线观看的亚洲视频| 国产亚洲精品综合一区在线观看 | 亚洲av美国av| 成年免费大片在线观看| 国产视频一区二区在线看| 波多野结衣高清作品| 日本一区二区免费在线视频| 国产精品一区二区精品视频观看| 19禁男女啪啪无遮挡网站| 欧美丝袜亚洲另类 | 色综合婷婷激情| 成人三级黄色视频| 久久性视频一级片| 日日夜夜操网爽| 最好的美女福利视频网| 高清毛片免费观看视频网站| 看黄色毛片网站| 天天躁狠狠躁夜夜躁狠狠躁| 免费看十八禁软件| 韩国av一区二区三区四区| 日本 av在线| 日韩高清综合在线| 观看免费一级毛片| 精品国产亚洲在线| 国产亚洲精品久久久久5区| 久久香蕉激情| 99久久99久久久精品蜜桃| 亚洲精品久久成人aⅴ小说| 午夜福利视频1000在线观看| 波多野结衣巨乳人妻| 美女午夜性视频免费| 黄片大片在线免费观看| 国产亚洲精品av在线| 日韩欧美一区视频在线观看| 好男人在线观看高清免费视频 | 亚洲欧美一区二区三区黑人| 男女那种视频在线观看| 亚洲欧美日韩无卡精品| 欧美乱妇无乱码| 99久久99久久久精品蜜桃| 亚洲片人在线观看| 精品一区二区三区四区五区乱码| 精品不卡国产一区二区三区| 叶爱在线成人免费视频播放| 精品国内亚洲2022精品成人| 国产在线观看jvid| 在线观看日韩欧美| 精品不卡国产一区二区三区| 久久久精品欧美日韩精品| 精品久久蜜臀av无| 老熟妇仑乱视频hdxx| 制服丝袜大香蕉在线| 日本在线视频免费播放| 好男人电影高清在线观看| 色播亚洲综合网| 怎么达到女性高潮| 少妇被粗大的猛进出69影院| 国产精品亚洲av一区麻豆| 又紧又爽又黄一区二区| 香蕉丝袜av| 亚洲男人的天堂狠狠| 老司机靠b影院| 精品免费久久久久久久清纯| 国产1区2区3区精品| 99re在线观看精品视频| 国内精品久久久久久久电影| 夜夜看夜夜爽夜夜摸| 亚洲狠狠婷婷综合久久图片| 男女那种视频在线观看| 午夜a级毛片| 可以在线观看的亚洲视频| 欧美不卡视频在线免费观看 | 精品熟女少妇八av免费久了| 侵犯人妻中文字幕一二三四区| 视频区欧美日本亚洲| 99热这里只有精品一区 | 国产精品野战在线观看| 一a级毛片在线观看| 99精品久久久久人妻精品| 岛国在线观看网站| 久久久国产精品麻豆| 成人18禁在线播放| 国产精品影院久久| 在线观看午夜福利视频| 久久久久国产精品人妻aⅴ院| 老汉色av国产亚洲站长工具| 校园春色视频在线观看| 91成人精品电影| 午夜福利视频1000在线观看| 亚洲国产精品合色在线| 国产一区二区激情短视频| 中出人妻视频一区二区| 久久国产精品影院| 亚洲精品粉嫩美女一区| 91麻豆av在线| 欧美成人午夜精品| 精品国产乱码久久久久久男人| 日本一本二区三区精品| 国产免费av片在线观看野外av| 女警被强在线播放| 亚洲三区欧美一区| 亚洲aⅴ乱码一区二区在线播放 | 2021天堂中文幕一二区在线观 | av天堂在线播放| 中亚洲国语对白在线视频| 亚洲av片天天在线观看| 亚洲熟妇熟女久久| 91大片在线观看| 男人舔奶头视频| 欧美成人一区二区免费高清观看 | xxx96com| 久久国产精品男人的天堂亚洲| 精品欧美国产一区二区三| 色老头精品视频在线观看| 久久香蕉国产精品| 日本一本二区三区精品| 久久精品aⅴ一区二区三区四区| 国产真实乱freesex| 757午夜福利合集在线观看| 啦啦啦 在线观看视频| www.熟女人妻精品国产| 久久精品夜夜夜夜夜久久蜜豆 | 午夜久久久在线观看| 自线自在国产av| 丝袜人妻中文字幕| www日本在线高清视频| 久久久久久久午夜电影| 精品国产乱子伦一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 欧美大码av| 人人妻,人人澡人人爽秒播| 成年人黄色毛片网站| 老司机在亚洲福利影院| 男女之事视频高清在线观看| 一本综合久久免费| av在线天堂中文字幕| 中文字幕久久专区| 深夜精品福利| 亚洲成国产人片在线观看| 少妇 在线观看| 两个人免费观看高清视频| 黄片大片在线免费观看| 九色国产91popny在线| 久久香蕉精品热| 亚洲男人天堂网一区| 国产精华一区二区三区| 女人被狂操c到高潮| 窝窝影院91人妻| 动漫黄色视频在线观看| 日本 av在线| 欧美zozozo另类| 国产真实乱freesex| 黄色a级毛片大全视频| 亚洲熟女毛片儿| 一区福利在线观看| 亚洲av日韩精品久久久久久密| 国产亚洲av嫩草精品影院| 亚洲在线自拍视频| 亚洲中文字幕一区二区三区有码在线看 | 久久热在线av| 国产高清激情床上av| 亚洲免费av在线视频| 悠悠久久av| 久久精品成人免费网站| 老熟妇乱子伦视频在线观看| 欧美日韩亚洲综合一区二区三区_| 成人欧美大片| 国产精品九九99| 亚洲在线自拍视频| 他把我摸到了高潮在线观看| 国产精品永久免费网站| 亚洲第一av免费看| 精品日产1卡2卡| 免费看a级黄色片| 一级作爱视频免费观看| 一进一出抽搐gif免费好疼| 日韩大码丰满熟妇| 麻豆成人av在线观看| 亚洲av熟女| 女人爽到高潮嗷嗷叫在线视频| videosex国产| 亚洲av成人一区二区三| 午夜福利免费观看在线| 一区二区三区国产精品乱码| 女生性感内裤真人,穿戴方法视频| 制服人妻中文乱码| 99热只有精品国产| 两个人看的免费小视频| 亚洲第一青青草原| 日本在线视频免费播放| 精品国产美女av久久久久小说| 成熟少妇高潮喷水视频| av欧美777| 欧美一级毛片孕妇| 在线观看一区二区三区| 亚洲成国产人片在线观看| 日本五十路高清| 国产视频一区二区在线看| 日本 欧美在线| 又黄又爽又免费观看的视频| 免费人成视频x8x8入口观看| 亚洲色图av天堂| 亚洲欧美一区二区三区黑人| 日韩精品中文字幕看吧| 久久久久久久精品吃奶| 天天一区二区日本电影三级| 久久久久久久精品吃奶| 少妇熟女aⅴ在线视频| 国产精品自产拍在线观看55亚洲| 中文资源天堂在线| 淫秽高清视频在线观看| 一级毛片精品| 女性生殖器流出的白浆| 欧美日本视频| 国产黄色小视频在线观看| 久久精品亚洲精品国产色婷小说| 俄罗斯特黄特色一大片| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩乱码在线| 欧美日本视频| 精品免费久久久久久久清纯| 色老头精品视频在线观看| 欧美日韩乱码在线| 757午夜福利合集在线观看| а√天堂www在线а√下载| 美国免费a级毛片| 俄罗斯特黄特色一大片| 欧美乱码精品一区二区三区| 精品免费久久久久久久清纯| 免费无遮挡裸体视频| 欧美在线一区亚洲| 在线观看舔阴道视频| 嫩草影院精品99| 亚洲国产日韩欧美精品在线观看 | 性欧美人与动物交配| 亚洲电影在线观看av| 免费搜索国产男女视频| 女同久久另类99精品国产91| 这个男人来自地球电影免费观看| 夜夜看夜夜爽夜夜摸| 禁无遮挡网站| 一区二区三区高清视频在线| 三级毛片av免费| 好男人电影高清在线观看| 国产成人一区二区三区免费视频网站| 国产三级在线视频| 性欧美人与动物交配| 国产精品久久久av美女十八| 亚洲男人的天堂狠狠| 久久婷婷人人爽人人干人人爱| 高清在线国产一区| 国产精品香港三级国产av潘金莲| 国产精品 欧美亚洲| 亚洲国产看品久久| √禁漫天堂资源中文www| 国产高清激情床上av| 久久久久精品国产欧美久久久| 国产私拍福利视频在线观看| 免费在线观看日本一区| 亚洲中文av在线| 99国产精品99久久久久| 亚洲人成77777在线视频| 久久久久免费精品人妻一区二区 | 免费一级毛片在线播放高清视频| 欧美性猛交╳xxx乱大交人| 国产成+人综合+亚洲专区| 免费观看人在逋| 免费在线观看视频国产中文字幕亚洲| 免费人成视频x8x8入口观看| 黄网站色视频无遮挡免费观看| 国产色视频综合| 老司机靠b影院| 国产麻豆成人av免费视频| 日本五十路高清| 久久精品国产99精品国产亚洲性色| 亚洲国产毛片av蜜桃av| 熟女少妇亚洲综合色aaa.| 久久婷婷人人爽人人干人人爱| 日韩欧美国产在线观看| 男女床上黄色一级片免费看| 国产精品野战在线观看| 亚洲狠狠婷婷综合久久图片| www日本在线高清视频| 日本 欧美在线| 亚洲三区欧美一区| 麻豆成人午夜福利视频| 欧美又色又爽又黄视频| 少妇的丰满在线观看| 精品国产一区二区三区四区第35| 午夜久久久久精精品| 亚洲国产欧洲综合997久久, | 一夜夜www| 午夜激情av网站| 一区二区三区精品91| 国产精品综合久久久久久久免费| 国产97色在线日韩免费| 搡老岳熟女国产| 国产精品久久电影中文字幕| 亚洲av美国av| 精品卡一卡二卡四卡免费| 黄频高清免费视频| 国产精品99久久99久久久不卡| 男人的好看免费观看在线视频 | 亚洲精品美女久久av网站| 男男h啪啪无遮挡| 少妇裸体淫交视频免费看高清 | 日韩免费av在线播放| aaaaa片日本免费| 国产又爽黄色视频| 久久久久精品国产欧美久久久| 操出白浆在线播放| 亚洲专区中文字幕在线| 欧美中文日本在线观看视频| 午夜福利一区二区在线看| 久久人妻av系列| 久久久久久国产a免费观看| 欧美乱色亚洲激情| 一边摸一边抽搐一进一小说| 日韩欧美 国产精品| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品sss在线观看| 国内久久婷婷六月综合欲色啪| 日韩欧美三级三区| 亚洲午夜理论影院| 国产精品 国内视频| 不卡一级毛片| 十分钟在线观看高清视频www| 国内揄拍国产精品人妻在线 | 久久婷婷人人爽人人干人人爱| 国产成人av激情在线播放| 亚洲人成77777在线视频| 亚洲国产精品合色在线| 亚洲成av人片免费观看| 亚洲精品一区av在线观看| 成人三级做爰电影| 亚洲精品色激情综合| 亚洲av中文字字幕乱码综合 | 在线观看www视频免费| 久久久精品国产亚洲av高清涩受| 男人舔女人下体高潮全视频| 精品国产乱子伦一区二区三区| 亚洲性夜色夜夜综合| 两性午夜刺激爽爽歪歪视频在线观看 | 久久九九热精品免费| 成人国产一区最新在线观看| 91麻豆av在线| 国产精品一区二区三区四区久久 | 19禁男女啪啪无遮挡网站| 久久精品亚洲精品国产色婷小说| 日韩大码丰满熟妇| 久久久久亚洲av毛片大全| 老司机福利观看| 亚洲精品色激情综合| 91在线观看av| 久久精品成人免费网站| 日韩精品青青久久久久久| 亚洲av电影不卡..在线观看| 亚洲男人天堂网一区| 久久精品国产综合久久久| 亚洲国产欧美一区二区综合| 精品电影一区二区在线| 亚洲avbb在线观看| 岛国视频午夜一区免费看| 黑人操中国人逼视频| 日韩av在线大香蕉| 悠悠久久av| 亚洲第一欧美日韩一区二区三区| 免费电影在线观看免费观看| 91九色精品人成在线观看| 欧美成人一区二区免费高清观看 | 亚洲一区二区三区色噜噜| 精品一区二区三区四区五区乱码| 黄色女人牲交| 国产1区2区3区精品| 亚洲aⅴ乱码一区二区在线播放 | 国产精品九九99| 亚洲国产精品合色在线| 麻豆成人午夜福利视频| 中文字幕人成人乱码亚洲影| 夜夜躁狠狠躁天天躁| 黄频高清免费视频| 亚洲精品中文字幕在线视频| 欧美一级毛片孕妇| АⅤ资源中文在线天堂| 岛国在线观看网站| 一级作爱视频免费观看| av在线播放免费不卡| 亚洲 欧美一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲av高清不卡| av有码第一页| 搞女人的毛片| 好男人在线观看高清免费视频 | 91老司机精品| 免费女性裸体啪啪无遮挡网站| 国产亚洲精品久久久久久毛片| 亚洲欧美激情综合另类| 国产私拍福利视频在线观看| 香蕉久久夜色| 高清毛片免费观看视频网站| 2021天堂中文幕一二区在线观 | 亚洲狠狠婷婷综合久久图片| 国产亚洲av高清不卡| 在线看三级毛片| 国产欧美日韩一区二区三| 亚洲国产欧美网| 又大又爽又粗| 久久人人精品亚洲av| 最新美女视频免费是黄的| 免费看美女性在线毛片视频| 999久久久精品免费观看国产| 欧美丝袜亚洲另类 | 国产高清激情床上av| 精品久久久久久久人妻蜜臀av| av超薄肉色丝袜交足视频| 91字幕亚洲| 国产片内射在线| 极品教师在线免费播放| 精品一区二区三区视频在线观看免费| 人人妻人人澡人人看| 他把我摸到了高潮在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲五月婷婷丁香| 亚洲自拍偷在线| 高潮久久久久久久久久久不卡| 性色av乱码一区二区三区2| 麻豆久久精品国产亚洲av| 熟女少妇亚洲综合色aaa.| 久久久久久久午夜电影| 日本成人三级电影网站| 女人被狂操c到高潮| 国产精品,欧美在线| 精品卡一卡二卡四卡免费| 欧美成人性av电影在线观看| 亚洲欧美精品综合一区二区三区| tocl精华| 一a级毛片在线观看| 黄色毛片三级朝国网站| 少妇裸体淫交视频免费看高清 | 亚洲 欧美 日韩 在线 免费| 久久久国产成人精品二区| 日本精品一区二区三区蜜桃| 亚洲午夜精品一区,二区,三区| 三级毛片av免费| 亚洲第一青青草原| 少妇熟女aⅴ在线视频| 黄色片一级片一级黄色片| 少妇被粗大的猛进出69影院| 国产伦一二天堂av在线观看| 亚洲成av人片免费观看| 国产精品香港三级国产av潘金莲| 亚洲国产欧美网| 欧美成狂野欧美在线观看| av欧美777| 韩国av一区二区三区四区| 精品欧美国产一区二区三| 欧美黑人精品巨大| 久久久久免费精品人妻一区二区 | 女同久久另类99精品国产91| 欧美中文日本在线观看视频| 久久香蕉国产精品| 国产亚洲精品久久久久久毛片| 免费看日本二区| 国产成年人精品一区二区| 亚洲av五月六月丁香网| 免费在线观看成人毛片| 亚洲精品美女久久久久99蜜臀| 欧美日韩精品网址| 欧美黄色淫秽网站| 国产一区二区在线av高清观看| 亚洲欧美精品综合一区二区三区| 在线av久久热| 亚洲精品粉嫩美女一区| 熟女电影av网| 国产亚洲欧美精品永久| 露出奶头的视频| 老汉色∧v一级毛片| www.熟女人妻精品国产| 亚洲av第一区精品v没综合|