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

    基于LSTM預估補償?shù)幕痣姍C組主蒸汽溫度控制系統(tǒng)

    2023-02-03 13:15:02張國斌郭瑞君杜榮華辛曉鋼牛玉廣
    發(fā)電設備 2023年1期
    關鍵詞:溫水設定值溫度控制

    張國斌, 郭瑞君, 杜榮華, 辛曉鋼, 張 謙, 牛玉廣

    (1. 內(nèi)蒙古電力(集團)有限責任公司 內(nèi)蒙古電力科學研究院分公司,呼和浩特 010020;2. 華北電力大學 控制與計算機工程學院,北京 102206)

    主蒸汽溫度是電站機組運行過程中需要控制的重要參數(shù)之一。主蒸汽溫度過高,會使汽輪機中汽缸、汽閥、葉片等部件的壽命縮短,進而影響汽輪機的安全性[1]。對于300 MW火電機組,主蒸汽溫度降低10 K,汽輪機熱耗率升高約0.33%(約25 kJ/(kW·h)),并且機組容量越大,主蒸汽溫度對汽輪機熱耗率的影響越大[2],同時汽輪機末級處蒸汽含水量增加,會造成汽輪機葉片壽命縮短。此外,溫度的波動會導致汽輪機內(nèi)汽缸和轉(zhuǎn)子脹差的變化,影響機組的安全運行。因此,快速、準確地控制主蒸汽溫度并將其穩(wěn)定在設定值十分有必要。

    隨著智能控制的發(fā)展,出現(xiàn)了模糊控制、智能比例積分微分(PID)控制、Smith預估控制和預測控制等先進控制方法[3-10]。黃宇等[3]針對600 MW直流爐主蒸汽溫度控制系統(tǒng),提出了模糊自適應內(nèi)模控制(FAIMC)策略,通過4種典型工況下的仿真表明,該策略的動態(tài)性能優(yōu)于串級控制策略的動態(tài)性能。LI H等[4]設計了一種基于免疫遺傳算法和反向傳播(BP)神經(jīng)網(wǎng)絡的智能PID控制系統(tǒng),利用免疫遺傳算法全局尋優(yōu)能力和良好的收斂性,優(yōu)化神經(jīng)網(wǎng)絡的權值,同時利用BP神經(jīng)網(wǎng)絡調(diào)整PID參數(shù),仿真結(jié)果表明,該系統(tǒng)在控制品質(zhì)和魯棒性方面均優(yōu)于常規(guī)PID控制系統(tǒng)。WANG G L等[5]介紹了一種用于超臨界機組主蒸汽溫度控制的多變量約束預測控制方法,將主蒸汽溫度階躍響應作為多變量約束預測控制算法的預測模型,仿真結(jié)果表明在負荷變化過程中,過熱蒸汽溫度可以控制在設定值附近。然而,這些方法主要是基于反饋調(diào)節(jié)的改進,反饋調(diào)節(jié)的特點是出現(xiàn)偏差才會調(diào)節(jié),對于主蒸汽溫度控制系統(tǒng)這類具有大慣性、大遲延特點的系統(tǒng)的應用效果有限。如果能夠提前預測主蒸汽溫度的變化方向,使得減溫水系統(tǒng)能夠提前動作,將能夠有效克服各種擾動對蒸汽溫度的影響。張麗香等[6]提出一種自適應Smith預估補償?shù)闹髡羝麥囟瓤刂葡到y(tǒng)。仿真研究結(jié)果表明,該系統(tǒng)設定值跟蹤性能和抗干擾能力都明顯優(yōu)于經(jīng)典串級汽主蒸汽壓力及溫度控制系統(tǒng)。但是,Smith預估器過于依賴過程精確模型,主蒸汽溫度受到機組負荷、各層磨煤機出力、給水流量和減溫水流量等眾多因素的影響,通過傳統(tǒng)的數(shù)學建模方法很難精確地構(gòu)建各個變量對主蒸汽溫度影響的綜合模型,而神經(jīng)網(wǎng)絡算法在理論上可以逼近對任意非線性函數(shù)的擬合,因此考慮采用機器學習模型對溫度進行預測,建立主蒸汽溫度控制系統(tǒng)預測模型。

    隨著機器學習算法的發(fā)展,支持向量機(SVM)、BP神經(jīng)網(wǎng)絡、循環(huán)神經(jīng)網(wǎng)絡(RNN)、長短期記憶網(wǎng)絡(LSTM)及其變體均被應用到時間序列預測中。梅倩[7]提出了一種有選擇性的增量學習與快速剪枝算法結(jié)合的最小二乘法支持向量機(LSSVM)預測方法。該方法既能保證預測的精度,又能滿足在線預測的實時性要求。LSSVM具有理論體系完善、適合小樣本訓練的特點,然而對于大規(guī)模樣本的訓練效果不足。神經(jīng)網(wǎng)絡具備擬合對象非線性特性的能力,在大規(guī)模樣本的時間序列預測上的效果更好。王德明等[11]采用優(yōu)化后的BP神經(jīng)網(wǎng)絡分別建立短期風速預測模型。然而,傳統(tǒng)的BP神經(jīng)網(wǎng)絡所有的輸入彼此獨立,沒有跨時間傳遞信息的能力,不能很好地利用其中的信息。

    RNN由Jordan M. I.于1986年提出[12],它不僅在層與層之間建立了權值連接,而且在隱含層之間的神經(jīng)元中也進行了權值傳遞,這樣就具有跨時間傳遞信息的能力,因此其更擅長對時間序列數(shù)據(jù)進行預測。李潔等[13]建立RNN預測模型,對未來時段的日客流量進行預測。然而,RNN存在梯度消失和長期依賴的問題,即不能學習到長期信息,不能選擇性地遺忘不重要的信息。LSTM采用門控機制解決了RNN中梯度消失和長期依賴的問題,能更有效地被應用到時間序列預測中。楊國田等[14]建立基于LSTM的火電廠氮氧化物(NOx)排放預測模型,得出LSTM模型數(shù)據(jù)預測效果明顯優(yōu)于RNN和LSSVM模型,并且預測結(jié)果的穩(wěn)定性和準確率較高。

    針對上述問題,筆者提出一種基于LSTM預估補償?shù)闹髡羝麥囟瓤刂葡到y(tǒng)。一方面,考慮各個變量對主蒸汽溫度的影響,建立基于LSTM的時間序列預測模型,將大部分可能引起主蒸汽溫度變化的變量作為預測模型的輸入,預測主蒸汽溫度的變化;另一方面,以主蒸汽溫度預測模型為基礎,結(jié)合階梯式預測控制算法,構(gòu)建基于預估補償?shù)闹髡羝麥囟瓤刂葡到y(tǒng),實現(xiàn)對主蒸汽溫度的精準超前調(diào)節(jié),提升主蒸汽溫度控制系統(tǒng)的效果。

    1 傳統(tǒng)主蒸汽溫度控制系統(tǒng)

    以某330 MW火電機組汽包爐兩級噴水減溫為研究對象,其結(jié)構(gòu)見圖1,其中:PV為反饋值,SP為設定值。

    圖1 某330 MW火電機組汽包爐主蒸汽溫度控制系統(tǒng)結(jié)構(gòu)

    汽包出口的飽和蒸汽分別經(jīng)過低溫過熱器加熱、一級減溫器降溫、屏式過熱器加熱、二級減溫器降溫和末級過熱器加熱后,進入高壓缸做功。主蒸汽溫度控制系統(tǒng)由兩段相對獨立的控制系統(tǒng)組成,分別為一級減溫水和二級減溫水系統(tǒng),并且每級減溫水系統(tǒng)可分為A/B兩側(cè),以二級減溫水系統(tǒng)A側(cè)控制為例說明傳統(tǒng)主蒸汽溫度控制系統(tǒng)。

    傳統(tǒng)主蒸汽溫度控制一般采用選用串級-前饋控制結(jié)構(gòu)。在二級減溫控制系統(tǒng)中,選用主蒸汽溫度作為主控制器的反饋值,將主控制器輸出作為副控制器設定值,將二級減溫后蒸汽溫度作為副控制器反饋值,副控制器的輸出直接控制減溫閥,二級減溫前蒸汽溫度作為二級減溫控制系統(tǒng)前饋。二級減溫前蒸汽溫度更能反映鍋爐燃燒狀況的變化,能更快地反映機組燃燒狀態(tài)對主蒸汽溫度的影響,補償燃燒狀態(tài)對主蒸汽溫度造成的干擾。

    然而,主蒸汽溫度受到多種因素的影響,難以采用單一信號表征主蒸汽溫度的變化,并且二級減溫前蒸汽溫度與被控蒸汽溫度在不同的受熱區(qū),不能只采用二級減溫前蒸汽溫度判斷主蒸汽溫度變化的趨勢。采用LSTM算法,訓練相關變量對主蒸汽溫度的時序關系模型,將主蒸汽溫度的預測值與當前值的差作為控制器的前饋信號,使控制器實現(xiàn)超前調(diào)節(jié)。

    2 LSTM算法原理

    LSTM算法由Hochreiter于1997年提出,并于2012年被Alex Graves進行了改良和推廣,是一種特殊的RNN算法[15]。LSTM算法采用門控機制解決了RNN算法中梯度消失和長期依賴的問題,能更有效地預測時間序列。LSTM網(wǎng)絡結(jié)構(gòu)見圖2,其中:t為當前時刻;σ為Sigmoid激活函數(shù);h(t)為當前時刻隱含層;C(t)為當前時刻存儲單元,也稱為細胞狀態(tài);f(t)、i(t)、o(t)分別為遺忘門、輸入門、輸出門;a(t)為輸入層和隱含層向存儲單元提供的更新值;X(t)為輸入值;y(t)為輸出值。

    圖2 LSTM網(wǎng)絡結(jié)構(gòu)圖

    在不同時刻的模塊中,LSTM通過隱含層和存儲單元傳遞權值信息,并通過門控單元選擇性地使信息通過,從而控制信息的遺忘和更新。門控單元不提供額外信息,只是起到信息過濾的作用。細胞狀態(tài)記錄了長時記憶傳播的信息,各個門控單元在細胞狀態(tài)的基礎上丟棄和增添信息,隱含層狀態(tài)記錄了短期記憶信息。

    每個LSTM模塊可以分為遺忘層、輸入層、更新層和輸出層四個部分。

    (1) 第一個部分為遺忘層,將上一時刻隱含層的輸出h(t-1)和當前時刻的輸入x(t)作為輸入,通過Sigmoid激活函數(shù)得到細胞狀態(tài)被選擇的概率f(t),其取值在[0,1](1表示完全保留細胞狀態(tài),0表示完全遺忘細胞狀態(tài))。

    (2) 第二個部分為輸入層,來確定將哪些信息存放到細胞狀態(tài)中,將輸入通過tanh激活函數(shù)后產(chǎn)生更新值a(t)。采用tanh激活函數(shù)是因為輸出在[-1,1],可以在更新細胞狀態(tài)的過程中加強或減弱記憶信息。通過a(t)與遺忘門輸出相乘的結(jié)果,可以表明細胞狀態(tài)需要更新的程度。

    (3) 第三個部分為更新層,表示細胞狀態(tài)的更新過程。將上一時刻的細胞狀態(tài)經(jīng)過遺忘層丟失一部分信息,再經(jīng)過輸入更新一部分信息,生成新的細胞狀態(tài)。

    (4) 第四個部分為輸出層,將細胞狀態(tài)經(jīng)過tanh激活函數(shù)得到輸出值,并經(jīng)過新的遺忘門篩選后送到輸出層預備值中,輸出層預備值也是下一時刻的隱含層輸入,將預備值送入Sigmoid激活函數(shù)得到預測輸出,從而實現(xiàn)時間序列的預測。

    LSTM算法相比于RNN算法主要增加了3個門控單元控制信息的流通以更新細胞狀態(tài),從而在網(wǎng)絡結(jié)構(gòu)上實現(xiàn)對長短期歷史數(shù)據(jù)信息的利用,既能學到長期信息,又能對不重要的近期信息選擇遺忘。

    圖3 LSTM算法的整體結(jié)構(gòu)模型

    在帶前饋的預測控制算法中,需要得到優(yōu)化時域內(nèi)所有的預測值,因此LSTM算法的整體結(jié)構(gòu)考慮到反向傳播時域內(nèi)所有時刻的誤差,而不是只考慮反向傳播最后時刻的誤差,從而保證同時得到優(yōu)化時域內(nèi)所有的預測值。

    采用這種結(jié)構(gòu),相比于BP算法,除了較高的準確度外,還具備另一個優(yōu)點,即能同時輸出優(yōu)化時域內(nèi)的所有預測值,并且每秒都在更新優(yōu)化時域內(nèi)的所有值,只需要在訓練時將整體結(jié)構(gòu)模型的長度設置成優(yōu)化時域,而BP算法每次只能輸出單時間的預測值,不能更新過去的預測值。

    采用Adam算法更新學習率,由Diederik Kingma在2015年提出的自適應學習率的算法,是一種可以替代傳統(tǒng)隨機梯度下降的優(yōu)化算法,根據(jù)損失函數(shù)梯度的一階矩估計和二階矩估計調(diào)整參數(shù)的學習率,具體公式[16]為:

    (1)

    3 基于LSTM預估補償?shù)闹髡羝麥囟瓤刂?/h2>

    3.1 主蒸汽溫度控制系統(tǒng)

    廣義預測控制由Clarke等于1987年在最小方差自校正的基礎上提出,采用多步預測優(yōu)化,最小化優(yōu)化時域內(nèi)的性能指標函數(shù),既不需要精確知道系統(tǒng)的滯后,也不需要對象的精準建模,適用于難以精確建模的滯后控制系統(tǒng)。

    圖4 帶預估補償前饋的主蒸汽溫度控制系統(tǒng)

    性能指標J的公式為:

    (2)

    式中:yr(t+j)為柔化后的設定值;y(t+j)為模型預測輸出;[N0,N1]為優(yōu)化時域;Nu為控制時域;λ為控制加權系數(shù),使控制增量的變化不要過于劇烈。控制算子由控制系統(tǒng)的優(yōu)化性能指標最小值計算得到,令dJ/d(Δu)=0,得:

    dT=(1,0,…,0)×(GTG+λ)-1GT

    (3)

    Δu=dT(Yr-Y)

    (4)

    Yr=[yr(t+N0),yr(t+N0+1),…,yr(t+N1)]T

    (5)

    Y=[y(t+N0),y(t+N0+1),…,y(t+N1)]T

    (6)

    式中:G為預測模型階躍響應系數(shù)組成的矩陣;Yr和Y分別為柔化后設定值和模型預測輸出組成的矩陣。

    基于LSTM預估補償?shù)念A測控制器從兩個方面對主蒸汽溫度控制進行改進。一方面,廣義預測控制器能根據(jù)對象模型預測溫度,從而提前控制減溫閥穩(wěn)定溫度。另一方面,對于可測不可控干擾,選用LSTM算法預測被控溫度,并通過前饋預先加以補償,使調(diào)節(jié)器能提前動作,克服控制遲延,進一步提升控制效果。

    3.2 主蒸汽溫度預測

    3.2.1 主蒸汽溫度對象特性

    為了使神經(jīng)網(wǎng)絡建立的預測器模型替代干擾對被控溫度的傳遞函數(shù),需要確定減溫水流量對主蒸汽溫度控制的遲延時間。在實際運行時主蒸汽溫度會受到各種因素的影響,因此選取穩(wěn)定工況的歷史曲線(見圖5)。從圖5中可以看出:一級減溫水流量對二級減溫前蒸汽溫度的控制遲延時間為70 s,二級減溫水流量對主蒸汽溫度的控制遲延時間為90 s。

    圖5 一級減溫和二級減溫對象的特性

    3.2.2 主蒸汽溫度預測模型

    一級減溫系統(tǒng)的控制目標為二級減溫前蒸汽溫度,二級減溫系統(tǒng)的控制目標為主蒸汽溫度,由于二級減溫前蒸汽溫度是否穩(wěn)定對主蒸汽溫度的影響很大,因此要將主蒸汽溫度控制好,對二級減溫前蒸汽溫度和主蒸汽溫度的預測和控制都很重要??刂葡到y(tǒng)的控制效果取決于前饋預測的準確性,需要對比BP算法和LSTM算法的預測效果,一級減溫和二級減溫控制系統(tǒng)參數(shù)見表1,其中給煤量和入口風量為5臺磨煤機(A、B、C、D、E磨煤機)的數(shù)據(jù),因此每組的輸入變量為16組,預測70 s后的二級減溫前蒸汽溫度和90 s后的主蒸汽溫度。

    表1 一級和二級減溫控制系統(tǒng)參數(shù)

    選取10 d的運行數(shù)據(jù),每隔5 s取1個點,共有172 800組數(shù)據(jù),將數(shù)據(jù)歸一化后輸入網(wǎng)絡進行訓練。進行預測訓練時,選取75%的數(shù)據(jù)作為訓練數(shù)據(jù),取25%的數(shù)據(jù)作為測試數(shù)據(jù),對比BP算法和LSTM算法的訓練效果。BP算法訓練采用Sigmoid激活函數(shù),輸入時域為100 s,即將16組輸入變量的100 s(20個點)的數(shù)據(jù)堆疊起來作為輸入,2層隱含層個數(shù)分別為64和16,訓練代數(shù)為100,初始學習率為0.5,學習輸入歷史數(shù)據(jù)的特征。LSTM算法訓練將16組數(shù)據(jù)直接輸入到當前時刻的模塊中,不需要將歷史數(shù)據(jù)堆疊,隱含層個數(shù)為64,訓練代數(shù)為100,初始學習率為0.5。2種方法設置成相同的訓練代數(shù)和初始學習率,預測70 s后的二級減溫前蒸汽溫度和90 s后的主蒸汽溫度。

    為了評價和對比預測效果,采用平均絕對值誤差EMA、均方誤差EMS、平均絕對百分誤差EMAP3個指標說明預測效果,一般認為EMA小于10%時的預測效果可以被接受。3個指標的公式分別為:

    (7)

    (8)

    (9)

    將BP算法和LSTM算法的神經(jīng)網(wǎng)絡訓練結(jié)果反歸一化后,比較兩者的預測效果(見表2)。

    表2 BP算法與LSTM算法的效果對比

    由表2可以看出:主蒸汽溫度的預測誤差比二級減溫前蒸汽溫度的小,并且LSTM算法的效果明顯優(yōu)于BP算法,適合離線訓練網(wǎng)絡參數(shù),深入學習歷史數(shù)據(jù)中的隱含信息。選取3 000 個樣本的訓練集和2 000 個樣本的測試集作為LSTM算法的訓練效果(見圖6)。

    圖6 LSTM算法的預測效果

    由圖6可以看出:二級減溫前蒸汽溫度和主蒸汽溫度的測試集中實際值與預測值變化趨勢的吻合程度較高,誤差較小,有很好的魯棒性。

    圖7為電廠中采用LSTM算法的實際預測效果。圖7(a)中將預測值推遲70 s后,其曲線與二級減溫前蒸汽溫度實際值曲線相吻合;圖7(b)中將預測值推遲90 s后,其曲線與主蒸汽溫度實際值曲線相吻合。這說明電廠實際運行時,預測值能提前反映機組參數(shù)狀態(tài),可將其用于控制系統(tǒng)的預估補償前饋中,進而起到提前控制的作用。

    圖7 電廠實際預測效果

    3.3 主蒸汽溫度控制效果

    改進前電廠采用常規(guī)的串級PID控制策略對主蒸汽溫度進行控制,選取其中1 h的控制效果圖(見圖8),其中:機組負荷(MW)由[100,300]歸一化到[0,100],一級減溫水質(zhì)量流量(t/h)由[20,135]歸一化到[0,100],二級減溫水質(zhì)量流量(t/h)由[0,50]歸一化到[0,100]。

    圖8 改進前蒸汽溫度控制效果

    由圖8可以看出:在10:54:00至11:06:00,機組負荷在211~228 MW波動,主蒸汽溫度在設定值附近振蕩,二級減溫前蒸汽溫度波動最大偏差為7.3 K,由于二級減溫前蒸汽溫度存在偏差,即使二級減溫水的閥門完全關閉,主蒸汽溫度也一直低于設定值。在11:06:00后,機組負荷由220 MW上升到250 MW,運行人員手動降低二級減溫前蒸汽溫度的設定值防止主蒸汽溫度超溫??偠灾捎迷械目刂撇呗詴r溫度偏差大,甚至需要運行人員的手動干預。將基于LSTM預估補償?shù)闹髡羝麥囟瓤刂撇呗詰糜陔姀S控制中,分析電廠的實際控制效果,其歸一化方式與常規(guī)的串級PID控制策略相同。

    改進后穩(wěn)態(tài)時主蒸汽溫度的控制效果見圖9。

    圖9 改進后穩(wěn)態(tài)時蒸汽溫度控制效果

    由圖9可以看出:在改進后的控制策略中,LSTM預估器能夠預測蒸汽溫度的變化方向,從而利用控制器對減溫水進行超前調(diào)節(jié),避免了主蒸汽溫度的大幅度波動。該電廠中一級減溫控制系統(tǒng)響應相對較快,并且一級減溫水流量的調(diào)節(jié)范圍廣,因此主要對二級減溫前蒸汽溫度超前調(diào)節(jié),而二級減溫水流量起到穩(wěn)定主蒸汽溫度的作用。以圖9中的二級減溫前蒸汽溫度為例進行分析。在06:32:00時,二級減溫前蒸汽溫度高于設定值且在溫度的高點時,一級減溫水流量提前降低,防止后續(xù)二級減溫前蒸汽溫度欠溫。同理,在06:12:00時,主蒸汽溫度高于設定值且在高點時,二級減溫水流量提前降低,防止06:14:00時主蒸汽溫度欠溫,最終使主蒸汽溫度更快地收斂到設定值上。當機組負荷在210~232 MW波動時,二級減溫前蒸汽溫度波動最大偏差為3.48 K,主蒸汽溫度波動最大偏差為3.06 K,說明了基于LSTM預估補償?shù)闹髡羝麥囟瓤刂撇呗缘挠行浴?/p>

    改進后動態(tài)時主蒸汽溫度的控制效果見圖10。

    圖10 改進后動態(tài)時蒸汽溫度控制效果

    由圖10可以看出:在09:19:00時,二級減溫前蒸汽溫度的設定值由497 ℃調(diào)整到490 ℃,主蒸汽溫度設定值由538 ℃調(diào)整到533 ℃,此時減溫水流量快速提升,然而機組負荷下降,二級減溫前蒸汽溫度和主蒸汽溫度會隨著機組負荷的下降而下降。因此,減溫水的閥門需要提前關小一部分,防止蒸汽溫度超調(diào),同時在調(diào)節(jié)過程中抑制燃燒狀態(tài)造成的干擾,說明了動態(tài)過程中基于LSTM預估補償?shù)闹髡羝麥囟瓤刂撇呗缘挠行浴?/p>

    4 結(jié)語

    針對具有大慣性、大延遲特點的主蒸汽溫度控制系統(tǒng),基于LSTM預估補償?shù)闹髡羝麥囟瓤刂葡到y(tǒng)能夠預測蒸汽溫度變化的方向,對減溫水進行超前調(diào)節(jié),避免了主蒸汽溫度大幅度波動,提升控制系統(tǒng)的動態(tài)特性和控制效果。電廠的實際應用表明基于LSTM預估補償?shù)闹髡羝麥囟瓤刂葡到y(tǒng)的應用效果優(yōu)于常規(guī)控制系統(tǒng)的應用效果,能使主蒸汽溫度更快穩(wěn)定在設定值。

    猜你喜歡
    溫水設定值溫度控制
    為什么要用溫水洗腳
    船舶主機溫度控制方法優(yōu)化
    冷軋鍍鋅光整機延伸率控制模式的解析與優(yōu)化
    為什么南極會有溫水湖泊?
    鍛造溫度控制
    等溫水移熱變換裝置生產(chǎn)運行總結(jié)
    孫紅雷:請把冷水換成溫水
    中學生(2016年13期)2016-12-01 07:03:52
    基于模糊控制的PLC在溫度控制中的應用
    電子制作(2016年21期)2016-05-17 03:53:01
    西氣東輸二線人機交互界面的防呆系統(tǒng)
    自動化儀表(2015年5期)2015-06-15 19:01:34
    目標設定值傳感器的原理與檢修
    亚洲欧美成人精品一区二区| 亚洲激情五月婷婷啪啪| 精品一区二区免费观看| 王馨瑶露胸无遮挡在线观看| 亚洲伊人色综图| 亚洲久久久国产精品| 免费av不卡在线播放| 国产av精品麻豆| 免费大片黄手机在线观看| 观看美女的网站| 免费不卡的大黄色大毛片视频在线观看| 亚洲性久久影院| 精品亚洲成a人片在线观看| 欧美成人精品欧美一级黄| 欧美最新免费一区二区三区| 久久毛片免费看一区二区三区| 乱码一卡2卡4卡精品| 亚洲国产看品久久| 亚洲国产日韩一区二区| 色5月婷婷丁香| 巨乳人妻的诱惑在线观看| 亚洲av在线观看美女高潮| 久久精品久久久久久久性| 嫩草影院入口| 丰满乱子伦码专区| 综合色丁香网| 国产免费一区二区三区四区乱码| 亚洲欧美日韩卡通动漫| 五月玫瑰六月丁香| 国产精品国产三级专区第一集| 五月伊人婷婷丁香| 国产av一区二区精品久久| 久久亚洲国产成人精品v| 久久久久久久国产电影| 韩国精品一区二区三区 | 日韩视频在线欧美| 秋霞在线观看毛片| 另类精品久久| 在线观看www视频免费| 国产精品一区www在线观看| 日日摸夜夜添夜夜爱| 日韩欧美精品免费久久| 视频区图区小说| 欧美精品一区二区免费开放| 亚洲伊人久久精品综合| 久久99蜜桃精品久久| 最近最新中文字幕大全免费视频 | 免费高清在线观看视频在线观看| 高清黄色对白视频在线免费看| 99re6热这里在线精品视频| 少妇精品久久久久久久| 乱人伦中国视频| 国产精品嫩草影院av在线观看| 丝袜脚勾引网站| av女优亚洲男人天堂| 久久综合国产亚洲精品| 日韩免费高清中文字幕av| 精品一区在线观看国产| 国国产精品蜜臀av免费| 国产成人精品婷婷| 又粗又硬又长又爽又黄的视频| 成人18禁高潮啪啪吃奶动态图| 亚洲高清免费不卡视频| 久久久久网色| 如何舔出高潮| 美女中出高潮动态图| 精品人妻在线不人妻| 秋霞在线观看毛片| 国产色爽女视频免费观看| 黑人巨大精品欧美一区二区蜜桃 | 日本与韩国留学比较| 成人国产av品久久久| 午夜免费鲁丝| av一本久久久久| 久久综合国产亚洲精品| 亚洲一级一片aⅴ在线观看| 午夜av观看不卡| 日本黄大片高清| 亚洲精品视频女| 欧美精品人与动牲交sv欧美| 亚洲欧美一区二区三区黑人 | 欧美精品一区二区大全| 亚洲欧美一区二区三区黑人 | 桃花免费在线播放| 国产精品三级大全| 精品国产国语对白av| 在线观看免费高清a一片| 老司机亚洲免费影院| 激情视频va一区二区三区| 亚洲国产精品国产精品| 亚洲欧美成人精品一区二区| 亚洲性久久影院| 熟女电影av网| 最近中文字幕2019免费版| 国产精品偷伦视频观看了| 少妇人妻久久综合中文| 久久97久久精品| 亚洲国产精品成人久久小说| 午夜精品国产一区二区电影| 久久毛片免费看一区二区三区| 男女免费视频国产| 视频区图区小说| 色哟哟·www| 一二三四中文在线观看免费高清| 另类精品久久| 国产探花极品一区二区| 国产永久视频网站| 一级毛片黄色毛片免费观看视频| 亚洲av日韩在线播放| 欧美精品一区二区免费开放| 久久这里只有精品19| 桃花免费在线播放| 中国国产av一级| 美国免费a级毛片| 国产亚洲一区二区精品| 国产色婷婷99| 黄色配什么色好看| 1024视频免费在线观看| 桃花免费在线播放| 国产av国产精品国产| 欧美日韩av久久| 亚洲国产精品国产精品| 亚洲久久久国产精品| 国产精品 国内视频| 亚洲一级一片aⅴ在线观看| 国产成人a∨麻豆精品| 中文字幕精品免费在线观看视频 | 波多野结衣一区麻豆| 亚洲情色 制服丝袜| 天天躁夜夜躁狠狠躁躁| 人成视频在线观看免费观看| 在线观看美女被高潮喷水网站| 香蕉精品网在线| 国产av精品麻豆| 高清av免费在线| 麻豆精品久久久久久蜜桃| 免费大片黄手机在线观看| 一级毛片 在线播放| 亚洲精品日本国产第一区| 亚洲欧美成人精品一区二区| 亚洲综合色惰| 丝袜脚勾引网站| 伦精品一区二区三区| 中文字幕人妻熟女乱码| av女优亚洲男人天堂| 人妻一区二区av| 国产极品粉嫩免费观看在线| 五月天丁香电影| 亚洲欧美成人精品一区二区| 成人国产av品久久久| 久久综合国产亚洲精品| 国产一区有黄有色的免费视频| 男女边摸边吃奶| 人人澡人人妻人| 黑人巨大精品欧美一区二区蜜桃 | 婷婷成人精品国产| 最近最新中文字幕大全免费视频 | 老司机影院毛片| 亚洲av福利一区| 亚洲五月色婷婷综合| 久久久久精品人妻al黑| 亚洲国产精品专区欧美| 人妻系列 视频| 90打野战视频偷拍视频| 91久久精品国产一区二区三区| 国产无遮挡羞羞视频在线观看| videosex国产| 激情五月婷婷亚洲| 高清毛片免费看| 最近最新中文字幕大全免费视频 | 全区人妻精品视频| 午夜免费鲁丝| 精品亚洲成国产av| 男女边摸边吃奶| 一二三四中文在线观看免费高清| 精品国产一区二区三区四区第35| 欧美亚洲 丝袜 人妻 在线| 中文字幕av电影在线播放| 丝袜喷水一区| 欧美精品一区二区免费开放| 欧美日韩av久久| 国产在线视频一区二区| 亚洲内射少妇av| 女的被弄到高潮叫床怎么办| 欧美精品人与动牲交sv欧美| 最后的刺客免费高清国语| 久久精品aⅴ一区二区三区四区 | 男女边吃奶边做爰视频| 乱码一卡2卡4卡精品| 这个男人来自地球电影免费观看 | 国产成人av激情在线播放| 精品亚洲乱码少妇综合久久| 一区二区三区精品91| 亚洲精品av麻豆狂野| 男人添女人高潮全过程视频| 视频区图区小说| 国产精品一区www在线观看| 亚洲精品美女久久久久99蜜臀 | 午夜福利,免费看| 亚洲国产av新网站| 99热这里只有是精品在线观看| 夫妻性生交免费视频一级片| 亚洲情色 制服丝袜| 久久久a久久爽久久v久久| 18在线观看网站| 精品少妇内射三级| 99国产精品免费福利视频| 成人国产av品久久久| 午夜福利影视在线免费观看| 国产欧美亚洲国产| 2018国产大陆天天弄谢| 老女人水多毛片| 亚洲av中文av极速乱| 国产精品久久久久久av不卡| 亚洲性久久影院| 日韩三级伦理在线观看| 国产爽快片一区二区三区| 伊人亚洲综合成人网| 亚洲 欧美一区二区三区| 国产成人精品一,二区| 黄片无遮挡物在线观看| 国产精品久久久久成人av| 中国美白少妇内射xxxbb| 久久久久视频综合| 中文字幕制服av| 永久免费av网站大全| av国产精品久久久久影院| 午夜影院在线不卡| 菩萨蛮人人尽说江南好唐韦庄| 51国产日韩欧美| 亚洲经典国产精华液单| 久久久久久久精品精品| 午夜免费观看性视频| 国产极品粉嫩免费观看在线| 亚洲av综合色区一区| 欧美xxⅹ黑人| 毛片一级片免费看久久久久| 欧美激情 高清一区二区三区| 国产成人免费观看mmmm| 久久久久久久大尺度免费视频| 宅男免费午夜| 日本色播在线视频| 欧美日韩av久久| 亚洲av电影在线进入| 激情五月婷婷亚洲| 成年女人在线观看亚洲视频| 七月丁香在线播放| 中文字幕免费在线视频6| 一级,二级,三级黄色视频| 精品少妇内射三级| 亚洲国产精品999| 日韩三级伦理在线观看| 男女午夜视频在线观看 | 久久精品国产鲁丝片午夜精品| 亚洲av国产av综合av卡| 中文字幕av电影在线播放| 欧美变态另类bdsm刘玥| 国产精品国产三级国产av玫瑰| 亚洲成人一二三区av| 国产免费一级a男人的天堂| 午夜福利视频精品| 99国产精品免费福利视频| 亚洲欧洲日产国产| 亚洲欧美色中文字幕在线| 久久精品国产自在天天线| 午夜激情av网站| 成人亚洲欧美一区二区av| 欧美性感艳星| 秋霞伦理黄片| 99热网站在线观看| 高清毛片免费看| 黄片无遮挡物在线观看| 国产精品国产三级国产专区5o| 制服丝袜香蕉在线| 我要看黄色一级片免费的| 伦理电影大哥的女人| 亚洲av电影在线观看一区二区三区| 国产色爽女视频免费观看| 国产欧美日韩综合在线一区二区| 国产无遮挡羞羞视频在线观看| 伊人亚洲综合成人网| 精品久久久精品久久久| 国产高清不卡午夜福利| 一二三四中文在线观看免费高清| 日韩,欧美,国产一区二区三区| 91久久精品国产一区二区三区| 少妇精品久久久久久久| 日日撸夜夜添| 亚洲图色成人| 久久精品久久久久久噜噜老黄| 国产精品一二三区在线看| 午夜福利网站1000一区二区三区| www日本在线高清视频| 免费高清在线观看日韩| 免费观看性生交大片5| 赤兔流量卡办理| 中文精品一卡2卡3卡4更新| 亚洲美女搞黄在线观看| 中文字幕人妻熟女乱码| 中文字幕制服av| 欧美精品一区二区免费开放| 侵犯人妻中文字幕一二三四区| 香蕉丝袜av| 亚洲成人手机| 国产精品一区二区在线观看99| 国产成人av激情在线播放| 国产精品国产av在线观看| 亚洲综合色网址| 日日啪夜夜爽| 国产成人免费观看mmmm| 丁香六月天网| 亚洲一级一片aⅴ在线观看| 亚洲图色成人| 国产一区二区在线观看av| 999精品在线视频| 1024视频免费在线观看| 秋霞在线观看毛片| av在线观看视频网站免费| 亚洲国产精品成人久久小说| 亚洲av欧美aⅴ国产| 欧美 亚洲 国产 日韩一| 大陆偷拍与自拍| 亚洲成人手机| 精品酒店卫生间| xxxhd国产人妻xxx| 久久久久国产网址| 亚洲精品av麻豆狂野| 国产探花极品一区二区| av国产精品久久久久影院| 欧美人与善性xxx| 伦理电影免费视频| 亚洲国产看品久久| 国产国拍精品亚洲av在线观看| 亚洲欧美色中文字幕在线| 女人精品久久久久毛片| 涩涩av久久男人的天堂| 美女国产视频在线观看| 伊人久久国产一区二区| 成人国产av品久久久| 亚洲国产看品久久| 考比视频在线观看| 欧美日韩视频精品一区| 咕卡用的链子| 欧美激情 高清一区二区三区| 精品一区二区三区四区五区乱码 | 啦啦啦啦在线视频资源| 女人久久www免费人成看片| av一本久久久久| 丝袜喷水一区| 国产 一区精品| 9热在线视频观看99| 五月开心婷婷网| 精品午夜福利在线看| 一二三四在线观看免费中文在 | 日韩av不卡免费在线播放| 校园人妻丝袜中文字幕| 新久久久久国产一级毛片| 一区二区av电影网| 午夜精品国产一区二区电影| 亚洲美女视频黄频| 韩国精品一区二区三区 | 精品视频人人做人人爽| 欧美另类一区| 亚洲精品一二三| 制服人妻中文乱码| 免费在线观看完整版高清| 成人影院久久| 色5月婷婷丁香| 免费大片18禁| 久久鲁丝午夜福利片| 国精品久久久久久国模美| 99re6热这里在线精品视频| 视频中文字幕在线观看| 深夜精品福利| 久久 成人 亚洲| 欧美激情极品国产一区二区三区 | 成年美女黄网站色视频大全免费| 老司机影院成人| 高清黄色对白视频在线免费看| 在线 av 中文字幕| 国产精品久久久久久久电影| 国产成人av激情在线播放| 国产色婷婷99| 亚洲精品日韩在线中文字幕| 一区二区三区四区激情视频| 午夜影院在线不卡| 免费观看a级毛片全部| 久久国内精品自在自线图片| 亚洲精品一二三| 成人国产av品久久久| 午夜福利,免费看| 亚洲成国产人片在线观看| 国产精品久久久久久精品古装| 国产有黄有色有爽视频| 国产精品女同一区二区软件| 午夜福利视频在线观看免费| 国产一区二区三区av在线| 国内精品宾馆在线| 欧美亚洲日本最大视频资源| 欧美日本中文国产一区发布| av又黄又爽大尺度在线免费看| a级毛片黄视频| 99re6热这里在线精品视频| 免费久久久久久久精品成人欧美视频 | 在线观看人妻少妇| 在线观看三级黄色| 午夜免费男女啪啪视频观看| 黄色怎么调成土黄色| 精品少妇黑人巨大在线播放| 国产免费视频播放在线视频| 中文字幕人妻丝袜制服| 久久97久久精品| 在线观看免费视频网站a站| 又大又黄又爽视频免费| 啦啦啦视频在线资源免费观看| 日本欧美视频一区| 婷婷色综合大香蕉| 国产 精品1| 国产av精品麻豆| 看免费av毛片| 成年女人在线观看亚洲视频| 国产一区二区在线观看av| 99热这里只有是精品在线观看| 欧美xxⅹ黑人| 嫩草影院入口| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧美精品自产自拍| 精品少妇久久久久久888优播| 成年女人在线观看亚洲视频| 久久这里只有精品19| 最近2019中文字幕mv第一页| 免费av不卡在线播放| 最近中文字幕高清免费大全6| 精品视频人人做人人爽| 久久这里只有精品19| 亚洲精品成人av观看孕妇| 精品熟女少妇av免费看| 国产成人a∨麻豆精品| 黑丝袜美女国产一区| 日本欧美国产在线视频| 一级,二级,三级黄色视频| 热re99久久精品国产66热6| 99热国产这里只有精品6| 久久精品国产a三级三级三级| 日韩中文字幕视频在线看片| 80岁老熟妇乱子伦牲交| 欧美日韩视频高清一区二区三区二| 国产精品.久久久| 校园人妻丝袜中文字幕| av视频免费观看在线观看| 99香蕉大伊视频| 日韩 亚洲 欧美在线| 亚洲成人av在线免费| 91久久精品国产一区二区三区| 又黄又爽又刺激的免费视频.| 观看av在线不卡| 少妇熟女欧美另类| 成年美女黄网站色视频大全免费| 卡戴珊不雅视频在线播放| 日韩成人伦理影院| √禁漫天堂资源中文www| 久久久久久久亚洲中文字幕| 制服诱惑二区| 久久久久国产精品人妻一区二区| 国产av码专区亚洲av| 黄色毛片三级朝国网站| 国产成人一区二区在线| 多毛熟女@视频| 欧美97在线视频| av免费在线看不卡| 高清毛片免费看| 国产片特级美女逼逼视频| 久久婷婷青草| 日韩不卡一区二区三区视频在线| 麻豆乱淫一区二区| 国产有黄有色有爽视频| 国产成人免费无遮挡视频| 欧美激情 高清一区二区三区| 国产成人91sexporn| 午夜精品国产一区二区电影| 久久免费观看电影| 免费日韩欧美在线观看| 人体艺术视频欧美日本| 午夜福利视频精品| 欧美激情极品国产一区二区三区 | 国产淫语在线视频| 一本色道久久久久久精品综合| 亚洲av在线观看美女高潮| 一级毛片 在线播放| 少妇被粗大的猛进出69影院 | 国产精品不卡视频一区二区| 少妇猛男粗大的猛烈进出视频| 日本黄大片高清| 欧美丝袜亚洲另类| 中国三级夫妇交换| 中国美白少妇内射xxxbb| 美女内射精品一级片tv| 亚洲av日韩在线播放| 亚洲精品第二区| 蜜臀久久99精品久久宅男| 99热全是精品| 亚洲激情五月婷婷啪啪| 色婷婷久久久亚洲欧美| 美女脱内裤让男人舔精品视频| 99九九在线精品视频| 99视频精品全部免费 在线| 少妇 在线观看| 国产一区二区在线观看日韩| 男人添女人高潮全过程视频| 看免费av毛片| 亚洲第一区二区三区不卡| 天堂俺去俺来也www色官网| 色婷婷av一区二区三区视频| 国产精品国产av在线观看| av卡一久久| 国产永久视频网站| 内地一区二区视频在线| av又黄又爽大尺度在线免费看| 精品一区二区三卡| 精品午夜福利在线看| 少妇精品久久久久久久| 亚洲av男天堂| 日韩制服丝袜自拍偷拍| 亚洲美女视频黄频| 日本欧美视频一区| 最近中文字幕高清免费大全6| 成年美女黄网站色视频大全免费| 亚洲av综合色区一区| 日韩三级伦理在线观看| 亚洲在久久综合| av不卡在线播放| 国产精品偷伦视频观看了| 91成人精品电影| 两个人看的免费小视频| 亚洲第一区二区三区不卡| 国产无遮挡羞羞视频在线观看| 一级毛片 在线播放| 午夜福利乱码中文字幕| 国产高清国产精品国产三级| 色5月婷婷丁香| 免费大片18禁| 久久久久久久久久成人| 国产视频首页在线观看| 午夜免费男女啪啪视频观看| 搡老乐熟女国产| 中国美白少妇内射xxxbb| 青青草视频在线视频观看| 不卡视频在线观看欧美| 久久国产亚洲av麻豆专区| 色网站视频免费| 亚洲欧美清纯卡通| 你懂的网址亚洲精品在线观看| 水蜜桃什么品种好| 精品久久久久久电影网| 尾随美女入室| 啦啦啦中文免费视频观看日本| 成人免费观看视频高清| 午夜老司机福利剧场| 99视频精品全部免费 在线| 咕卡用的链子| 免费女性裸体啪啪无遮挡网站| 亚洲四区av| 91午夜精品亚洲一区二区三区| 亚洲欧美中文字幕日韩二区| 精品人妻偷拍中文字幕| 日韩一区二区三区影片| 国产熟女午夜一区二区三区| 狠狠精品人妻久久久久久综合| 极品人妻少妇av视频| 免费少妇av软件| 伦理电影免费视频| 九草在线视频观看| 日韩 亚洲 欧美在线| 久久久久久伊人网av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 两个人看的免费小视频| 日本黄色日本黄色录像| 91在线精品国自产拍蜜月| 国产精品.久久久| 99久久精品国产国产毛片| 精品人妻偷拍中文字幕| 热re99久久精品国产66热6| 免费播放大片免费观看视频在线观看| 乱码一卡2卡4卡精品| 国产精品免费大片| 高清视频免费观看一区二区| 少妇猛男粗大的猛烈进出视频| 日韩人妻精品一区2区三区| 黑人巨大精品欧美一区二区蜜桃 | 日本欧美视频一区| 91久久精品国产一区二区三区| 午夜福利,免费看| 一区二区三区四区激情视频| 日韩熟女老妇一区二区性免费视频| 国语对白做爰xxxⅹ性视频网站| 欧美精品一区二区大全| 99久久人妻综合| 亚洲国产精品一区三区| 亚洲一码二码三码区别大吗| 老熟女久久久| 国产成人av激情在线播放| 美女大奶头黄色视频| 乱码一卡2卡4卡精品| av国产精品久久久久影院| 亚洲精品自拍成人| 午夜影院在线不卡| av卡一久久| 香蕉国产在线看| 精品视频人人做人人爽| 一级毛片我不卡| 最黄视频免费看| 国产一区有黄有色的免费视频| 国产视频首页在线观看| 国产成人精品福利久久| 国产在视频线精品|