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

    基于間隙度量的多模型過熱汽溫預(yù)測(cè)控制

    2022-02-17 02:53:04齊衛(wèi)祎鄧拓宇
    電力科學(xué)與工程 2022年1期
    關(guān)鍵詞:熱汽被控階躍

    齊衛(wèi)祎,鄧拓宇

    (華北電力大學(xué) 控制與計(jì)算機(jī)工程學(xué)院,河北 保定 071003)

    0 引 言

    隨著2030年“碳達(dá)峰”目標(biāo)的提出,大力發(fā)展新能源發(fā)電已經(jīng)成為必然。新能源發(fā)電存在隨機(jī)性與不穩(wěn)定性,需要火電機(jī)組來承擔(dān)更多深度調(diào)峰的任務(wù)。傳統(tǒng)的火電機(jī)組主汽溫串級(jí)控制器的控制參數(shù)是在一定工況下整定得來的。當(dāng)機(jī)組處于深度調(diào)峰的變負(fù)荷工況下,機(jī)組的運(yùn)行參數(shù)變化幅度較大;同時(shí),汽溫對(duì)象本身具有大慣性、時(shí)滯性特性,再加上擾動(dòng)的影響,這些因素會(huì)導(dǎo)致系統(tǒng)控制品質(zhì)的惡化,使得過熱汽溫的動(dòng)態(tài)偏差易超過允許安全范圍,進(jìn)而影響整個(gè)火電機(jī)組的安全運(yùn)行[1]。

    由于傳統(tǒng)的PID串級(jí)控制難以適應(yīng)當(dāng)前機(jī)組深度變負(fù)荷的情況,故諸多先進(jìn)的控制策略被引入主汽溫控制中來,如模糊控制[2]、自抗擾控制[3]、神經(jīng)網(wǎng)絡(luò)控制[4]和預(yù)測(cè)控制[5-6]。預(yù)測(cè)控制[7]自被提出以來,已經(jīng)在復(fù)雜工業(yè)過程中取得廣泛的應(yīng)用,其在處理工業(yè)過程中的大遲延、大慣性問題時(shí)顯示出強(qiáng)大的優(yōu)勢(shì)[8]。

    目前,有關(guān)學(xué)者已經(jīng)對(duì)預(yù)測(cè)控制在主汽溫調(diào)節(jié)中的應(yīng)用做了大量的研究。文獻(xiàn)[9]提出了一種基于ARX-Laguerre預(yù)測(cè)模型的預(yù)測(cè)控制算法,把PID控制器的性能指標(biāo)引入到預(yù)測(cè)控制的優(yōu)化指標(biāo)中去,通過在線最小二乘辨識(shí)來提升系統(tǒng)的自適應(yīng)性。文獻(xiàn)[10]將迭代預(yù)測(cè)控制器引入過熱汽溫控制中,通過歷史誤差來修正控制器,實(shí)現(xiàn)了模型失配情況下的設(shè)定值跟蹤。文獻(xiàn)[11]提出了多模型切換的預(yù)測(cè)控制方法,即根據(jù)不同工況預(yù)先設(shè)計(jì)不同的控制器,再在變負(fù)荷時(shí)切換到最匹配的預(yù)測(cè)函數(shù)控制器。文獻(xiàn)[12]提出了一種基于RBFNN的PFC控制策略,即通過徑向基神經(jīng)網(wǎng)絡(luò)實(shí)時(shí)辨識(shí)被控對(duì)象來克服被控對(duì)象的時(shí)變;但文中方法神經(jīng)網(wǎng)絡(luò)在線計(jì)算量大,不易于實(shí)際應(yīng)用。文獻(xiàn)[13]提出了主汽溫多模型預(yù)測(cè)控制策略,各控制器的權(quán)值計(jì)算方法是:計(jì)算實(shí)際汽溫與預(yù)測(cè)汽溫的差值,然后用遞推貝葉斯公式計(jì)算得到各子模型與實(shí)際被控模型的匹配概率。文獻(xiàn)[14]將PI控制與隱式GPC算法相結(jié)合,增強(qiáng)了控制器的魯棒性。文獻(xiàn)[15]基于多模型切換的鍋爐汽溫預(yù)測(cè)控制,首先選取多個(gè)典型工況設(shè)計(jì)動(dòng)態(tài)矩陣控制器,通過各子模型輸出與實(shí)際輸出的誤差計(jì)算 ITAE性能指標(biāo),并將其作為子控制器的切換準(zhǔn)則;但該方法使用時(shí)需要選取的子模型足夠多,否則各控制器在切換時(shí)會(huì)引入控制輸出的振蕩。

    間隙度量理論能準(zhǔn)確地描述不同模型之間的動(dòng)態(tài)特性差異。本文提出基于間隙度量理論的多模型過熱汽溫預(yù)測(cè)控制策略:根據(jù)若干典型工況的固定模型,設(shè)計(jì)相應(yīng)的局部控制器,計(jì)算得到每一個(gè)時(shí)刻的控制增量;通過間隙度量理論加權(quán)各個(gè)局部控制器在該時(shí)刻的控制增量,得到該時(shí)刻總控制增量并作用到被控對(duì)象上。仿真結(jié)果表明,該方法可以滿足火電機(jī)組全工況運(yùn)行的要求。

    1 過熱汽溫噴水減溫控制系統(tǒng)

    過熱汽溫控制系統(tǒng)由汽包、一級(jí)過熱器、屏式過熱器、二級(jí)過熱器和2級(jí)噴水減溫器組成。

    采用2級(jí)噴水減溫結(jié)構(gòu)對(duì)末級(jí)出口主汽溫進(jìn)行控制:一級(jí)噴水減溫起粗調(diào)作用,二級(jí)噴水減溫起細(xì)調(diào)作用。2個(gè)噴水減溫系統(tǒng)的原理相同。圖1所示為一級(jí)噴水減溫示意圖。

    圖1 過熱汽溫控制結(jié)構(gòu)Fig.1 Structure of superheated steam temperature control

    通過控制噴水減溫閥的開度來控制減溫水流量,進(jìn)而調(diào)控過熱蒸汽溫度。圖1中,減溫水流量W和噴水減溫器出口溫度θ1為導(dǎo)前區(qū),θ1到θ2之間的過熱器管道為惰性區(qū)。

    導(dǎo)前區(qū)傳遞函數(shù):

    惰性區(qū)傳遞函數(shù):

    2 多模型預(yù)測(cè)控制策略

    考慮過熱汽溫對(duì)象具有時(shí)變性,本文采用多模型切換方法。通過間隙度量理論加權(quán)多個(gè)控制器的輸出。對(duì)于大遲延、大慣性的被控對(duì)象,采用動(dòng)態(tài)矩陣控制。圖2示出本文采用的多模型預(yù)測(cè)控制策略。

    圖2 多模型預(yù)測(cè)控制策略Fig.2 Strategy of multi-model predictive control

    首先選取可以覆蓋整個(gè)電廠負(fù)荷變化的4個(gè)典型工況。根據(jù)這4個(gè)典型工況分別建立局部動(dòng)態(tài)矩陣控制器。局部預(yù)測(cè)控制器輸出為當(dāng)前這個(gè)時(shí)刻的控制增量Δui(k)(i=1,2,3,4)。當(dāng)負(fù)荷偏離這4個(gè)典型工況點(diǎn)的時(shí)候,通過間隙度量理論計(jì)算當(dāng)前工況偏離4個(gè)典型工況的程度,求得4個(gè)局部預(yù)測(cè)控制器的控制量權(quán)重wj(j=1,2,3,4)。

    2.1 多個(gè)子模型的選取

    對(duì)于過熱汽溫被控對(duì)象,影響模型參數(shù)變化最大的擾動(dòng)是負(fù)荷變化。當(dāng)火電機(jī)組負(fù)荷變化的時(shí)候,機(jī)組的燃煤量和給水量都會(huì)相應(yīng)的發(fā)生變化。機(jī)組給水量的變化會(huì)影響過熱器中的換熱,導(dǎo)致過熱汽溫發(fā)生變化;機(jī)組的燃煤量發(fā)生改變以后,直接改變產(chǎn)生的熱量,使得煙氣溫度變化,導(dǎo)致過熱汽溫改變。所以,需要選取幾個(gè)典型工況下的過熱汽溫被控對(duì)象,建立模型集。

    文獻(xiàn)[16]給出了在4個(gè)典型工況下,過熱汽溫被控對(duì)象導(dǎo)前區(qū)和惰性區(qū)的傳遞函數(shù),如表1所示。

    表1 不同負(fù)荷過熱汽溫參數(shù)Tab.1 Parameters of superheated steam temperature under different loads

    2.2 多模型切換控制權(quán)重的確定

    本文采用間隙度量理論[17]計(jì)算不同負(fù)荷下被控對(duì)象之間的距離,從而確定各控制器的輸出權(quán)重。

    假設(shè)M,N為巴那赫空間Z中的2個(gè)線性子空間。若K(s)是傳遞函數(shù)矩陣,則可以得到其正規(guī)化右互質(zhì)分解:

    線性模型K1和K2之間的間隙度量值:

    2個(gè)傳遞函數(shù)越相似,則間隙度量值越接近于0。間隙度量值與工況的相似程度呈反比。通過間隙度量值來確定各個(gè)控制器輸出的權(quán)重:

    由式(9)可知,當(dāng)運(yùn)行工況接近某一工況的時(shí)候,則該工況的局部控制器輸出最大。k為可調(diào)參數(shù),k大于等于1,通過k以調(diào)整輸出的權(quán)值。

    2.3 動(dòng)態(tài)矩陣控制

    動(dòng)態(tài)矩陣控制是一種基于系統(tǒng)階躍響應(yīng)的預(yù)測(cè)控制算法。因其以被控對(duì)象的階躍響應(yīng)系數(shù)為模型,所以避免了通常的傳遞函數(shù)模型參數(shù)的辨識(shí)過程,從而減少了計(jì)算量。采用多步預(yù)估技術(shù)能有效解決被控對(duì)象的遲延問題;按照使預(yù)估輸出和給定值偏差最小的二次性能指標(biāo)來計(jì)算控制增量,從而實(shí)現(xiàn)最優(yōu)控制;通過計(jì)算預(yù)測(cè)值和實(shí)際的偏差,不斷校正每一步的預(yù)測(cè)值,構(gòu)成一個(gè)閉環(huán),增強(qiáng)系統(tǒng)的穩(wěn)定性。

    2.3.1 預(yù)測(cè)模型

    圖3為被控對(duì)象的階躍響應(yīng)曲線圖。對(duì)象的動(dòng)態(tài)特性通過單位階躍響應(yīng)在采樣時(shí)刻的值來描述,如a1,a2,···,ap。p為模型的時(shí)域長(zhǎng)度,ap為接近穩(wěn)態(tài)值的系數(shù)。

    圖3 被控對(duì)象階躍響應(yīng)Fig.3 Step response of the controlled object

    通過階躍響應(yīng)可以采樣得到p個(gè)時(shí)刻的采樣值,進(jìn)而可以構(gòu)建反映系統(tǒng)動(dòng)態(tài)特性的動(dòng)態(tài)矩陣A。根據(jù)線性系統(tǒng)疊加的原理,可以由給定輸入的增量來預(yù)測(cè)系統(tǒng)未來的輸出。

    式中:n為最大預(yù)測(cè)長(zhǎng)度;m為控制長(zhǎng)度;ΔU為控制增量;Y0為控制作用保持不變時(shí)的初始預(yù)測(cè)值;為n步預(yù)測(cè)值。

    2.3.2 滾動(dòng)優(yōu)化

    滾動(dòng)優(yōu)化的目的是為了使輸出預(yù)測(cè)值逼近設(shè)定值,并且盡可能減少控制量的波動(dòng)。如式(15),構(gòu)建二次性能指標(biāo)函數(shù)來求取使得誤差最小和控制量最小的控制增量。

    式中:q(j)為誤差的加權(quán);r(j)為控制增量的加權(quán)。

    將式(15)用向量形式表示:

    式中:W為設(shè)定值的向量形式;Q誤差加權(quán)矩陣;R為控制增量加權(quán)矩陣。

    通過偏導(dǎo)為零求取極值點(diǎn),獲得控制增量值。

    優(yōu)化過程不是一次完成的。每次計(jì)算出最優(yōu)控制增量以后,取當(dāng)前的控制增量施加控制;每一步控制量只取當(dāng)前的控制增量,滾動(dòng)進(jìn)行優(yōu)化。

    2.3.3 誤差校正

    由于實(shí)際運(yùn)行中系統(tǒng)會(huì)存在擾動(dòng),使得預(yù)測(cè)值和對(duì)象的實(shí)際值產(chǎn)生較大的偏差,所以引入誤差校正。

    式(22)為k+1時(shí)刻實(shí)際輸出和預(yù)測(cè)輸出構(gòu)成的預(yù)測(cè)誤差。

    式中:h為誤差校正向量。

    在k+1時(shí)刻,時(shí)間基點(diǎn)的變動(dòng)導(dǎo)致預(yù)測(cè)未來時(shí)間點(diǎn)也將移動(dòng)到k+2,···,k+1+N;因此k+1時(shí)刻的初始預(yù)測(cè)值可以通過校正后預(yù)測(cè)值移位得到。

    誤差反饋校正的引入,使系統(tǒng)成為一個(gè)反饋控制系統(tǒng),提升了系統(tǒng)的穩(wěn)定性。

    3 仿真實(shí)驗(yàn)

    為了驗(yàn)證本文提出的多模型預(yù)測(cè)控制策略的有效性,通過 MATLAB進(jìn)行仿真分析:分別與PID串級(jí)控制、單模型預(yù)測(cè)控制進(jìn)行對(duì)比,從穩(wěn)態(tài)性能、動(dòng)態(tài)性能角度進(jìn)行分析。

    過熱汽溫模型的內(nèi)回路采用比例控制,可以快速消除內(nèi)回路的擾動(dòng)。通過 ZN整定法求取比例系數(shù)。將整定好的內(nèi)回路與惰性區(qū)傳遞函數(shù)作為廣義被控對(duì)象,對(duì)其進(jìn)行控制器的設(shè)計(jì)。多模型預(yù)測(cè)控制選用表1的4個(gè)典型工況作為預(yù)測(cè)模型。構(gòu)建動(dòng)態(tài)矩陣,建立4個(gè)局部動(dòng)態(tài)矩陣控制器??刂破鲄?shù)為采樣時(shí)間T=5 s,預(yù)測(cè)步長(zhǎng)P=200,控制步長(zhǎng)M=10,誤差加權(quán)矩陣Q=IP×P,控制量加權(quán)矩陣R=IM×M,誤差校正向量hi=1(i=1,···,P)。

    3.1 定負(fù)荷情況

    當(dāng)機(jī)組運(yùn)行在典型工況負(fù)荷下,多模型預(yù)測(cè)控制與單模型預(yù)測(cè)控制效果一致。將多模型預(yù)測(cè)控制與PID串級(jí)控制進(jìn)行比較。

    在負(fù)荷37%、50%、75%、100% 4種典型工況下,分別進(jìn)行設(shè)定值階躍變化實(shí)驗(yàn)。設(shè)定值起始值為10,在2 500 s變?yōu)?,仿真時(shí)間為5 000 s。PID串級(jí)控制的主控制器的參數(shù)采用ZN整定法。各負(fù)荷下,PID串級(jí)控制的主控制器參數(shù)如下:在 37%負(fù)荷時(shí),Kp=0.97,Ti=541.1,Td=135.27;在 50%負(fù)荷時(shí),Kp=0.99,Ti=371.69,Td=92.92;在 75%負(fù)荷時(shí),Kp=0.95,Ti=251.03,Td=62.76;在100%負(fù)荷時(shí),Kp=1.05,Ti=161.33,Td=40.33。圖4所示為2種控制方法的仿真結(jié)果。

    圖4 各典型工況階躍變化仿真曲線Fig.4 Simulation curves of step changes under typical operating conditions

    參考圖4,對(duì)比不同負(fù)荷下多模預(yù)測(cè)控制相較于 PID串級(jí)控制的超調(diào)量與調(diào)節(jié)時(shí)間:在37%負(fù)荷下,多模型預(yù)測(cè)控制的超調(diào)量為4.48%,減少了1.17%,調(diào)節(jié)時(shí)間為780 s,減少了1 570 s;在50%負(fù)荷下,多模型預(yù)測(cè)控制的超調(diào)量為3.06%,減少了4.63%,調(diào)節(jié)時(shí)間為685 s,減少了920 s;在75%負(fù)荷下,多模型預(yù)測(cè)控制的超調(diào)量為3.76%,減少了4.91%,調(diào)節(jié)時(shí)間為360 s,減少了720 s;在100%負(fù)荷下,多模型預(yù)測(cè)控制的超調(diào)量為4.90%,減少了8.11%,調(diào)節(jié)時(shí)間為255 s,減少了420 s。

    由上可以看出,多模型預(yù)測(cè)控制的調(diào)節(jié)時(shí)間和超調(diào)量都明顯小于PID串級(jí)控制,顯示出更好的動(dòng)態(tài)性能和穩(wěn)態(tài)性能。

    3.2 變工況情況

    為了驗(yàn)證多模型預(yù)測(cè)控制在負(fù)荷變化時(shí)的控制特性,做機(jī)組負(fù)荷遞減時(shí)的設(shè)定值階躍仿真實(shí)驗(yàn)。假設(shè)機(jī)組剛開始運(yùn)行在100%負(fù)荷下,在時(shí)間2 500 s,5 000 s,7 500 s,機(jī)組負(fù)荷分別遞減為 75%,50%,37%。將獲取的系統(tǒng)響應(yīng)曲線與在 100%負(fù)荷下設(shè)計(jì)的傳統(tǒng)串級(jí)控制器和單模型預(yù)測(cè)控制器響應(yīng)曲線進(jìn)行對(duì)比,如圖5所示。

    圖5 工況遞減下不同控制策略階躍響應(yīng)Fig.5 Step response of different control strategies under decreasing operating conditions

    在機(jī)組負(fù)荷突變時(shí),被控對(duì)象參數(shù)的變化導(dǎo)致模型失配,使得被控量發(fā)生波動(dòng),偏離設(shè)定值。由圖5可見:多模型預(yù)測(cè)控制在機(jī)組負(fù)荷每次驟變之后,能夠快速使被控量穩(wěn)定到設(shè)定值;而單模型預(yù)測(cè)控制隨著模型失配程度的增加,在每次機(jī)組負(fù)荷驟變之后,系統(tǒng)的被控量波動(dòng)越來越大,系統(tǒng)的穩(wěn)定性降低;傳統(tǒng)的串級(jí)控制器由于其魯棒性,在負(fù)荷大范圍波動(dòng)時(shí),控制器可以保證被控值穩(wěn)定不發(fā)散,但隨著負(fù)荷變化越大,其控制性能開始變差。

    考慮到在實(shí)際生產(chǎn)過程中,機(jī)組多數(shù)時(shí)間處于小范圍負(fù)荷波動(dòng)運(yùn)行工況,所以做工況小范圍波動(dòng)實(shí)驗(yàn),以觀察多模型預(yù)測(cè)控制、單模型預(yù)測(cè)控制和PID串級(jí)控制對(duì)于小范圍變工況的適應(yīng)程度。假設(shè)機(jī)組在100%負(fù)荷穩(wěn)定,而此時(shí)被控對(duì)象參數(shù)如表1所示。在仿真時(shí)間0時(shí)刻,負(fù)荷發(fā)生改變,被控對(duì)象參數(shù)發(fā)生變化:導(dǎo)前區(qū)參數(shù)K增益變?yōu)?.2,時(shí)間常數(shù)T變?yōu)?9;惰性區(qū)參數(shù)保持不變。做設(shè)定值階躍響應(yīng)擾動(dòng)實(shí)驗(yàn),階躍5個(gè)單位,獲得系統(tǒng)響應(yīng)曲線如圖6所示。

    圖6 變工況下不同控制策略階躍響應(yīng)Fig.6 Step response of different control strategies under variable conditions

    由圖6可以看出,多模型預(yù)測(cè)控制的超調(diào)量為1.51%,調(diào)節(jié)時(shí)間為480 s;單模型預(yù)測(cè)控制的超調(diào)量為23.92%,調(diào)節(jié)時(shí)間為780 s;PID串級(jí)控制的超調(diào)量為25.42%,調(diào)節(jié)時(shí)間為805 s。多模型預(yù)測(cè)控制在變工況的情況下能夠更加快速平穩(wěn)地使過熱汽溫達(dá)到設(shè)定值。

    可以看出,本文所提出的多模型預(yù)測(cè)控制能夠適應(yīng)被控對(duì)象參數(shù)的時(shí)變,具有較好的魯棒性。

    3.3 擾動(dòng)干擾情況

    為了驗(yàn)證多模型預(yù)測(cè)控制的抗干擾性能,選用100%負(fù)荷下的過熱汽溫模型,在2 500 s給設(shè)定值加10倍(t)的階躍擾動(dòng),并與PID串級(jí)控制進(jìn)行對(duì)比。響應(yīng)曲線如圖7所示。

    圖7 過熱汽溫輸入端擾動(dòng)階躍響應(yīng)Fig.7 Step response of input disturbance of superheated steam temperature

    從系統(tǒng)響應(yīng)曲線可以看出,在系統(tǒng)受到外界干擾的情況下,多模型預(yù)測(cè)控制和PID串級(jí)控制均存在超調(diào)現(xiàn)象;但多模型預(yù)測(cè)控制的調(diào)節(jié)時(shí)間更短,可以使系統(tǒng)更快恢復(fù)設(shè)定值,擁有更好的抗干擾性能,魯棒性更強(qiáng)。

    4 結(jié)論

    為應(yīng)對(duì)火電機(jī)組深度調(diào)峰,提出了基于間隙度量的多模型過熱汽溫預(yù)測(cè)控制策略。

    (1)在負(fù)荷不變的情況,多模型預(yù)測(cè)控制比PID串級(jí)控制更能快速穩(wěn)定完成對(duì)設(shè)定值的響應(yīng)。

    (2)在負(fù)荷變化的情況下,多模型預(yù)測(cè)控制可以更好地適應(yīng)負(fù)荷快速變化的情況,使過熱汽溫穩(wěn)定在一定范圍內(nèi)。

    (3)過熱汽溫系統(tǒng)輸入擾動(dòng)實(shí)驗(yàn)結(jié)果表明,相較與PID串級(jí)控制,多模型預(yù)測(cè)控制能夠使過熱汽溫值更加快速回到設(shè)定值。

    猜你喜歡
    熱汽被控階躍
    基于階躍雙包層光纖的螺旋型光纖傳感器
    基于遺傳算法的模糊控制在過熱汽溫控制系統(tǒng)優(yōu)化中的應(yīng)用
    電子制作(2019年16期)2019-09-27 09:34:44
    大慣量系統(tǒng)位置控制策略研究
    探討單位階躍信號(hào)的教學(xué)
    600MW亞臨界汽包爐再熱汽溫低的治理及應(yīng)用
    600MW超臨界直流爐過熱汽溫動(dòng)態(tài)特性試驗(yàn)與控制系統(tǒng)優(yōu)化
    一種階躍函數(shù)在矩形時(shí)間窗口頻域特性的分析方法
    單位階躍響應(yīng)的時(shí)域求解方法
    600 MW火電機(jī)組再熱汽溫的調(diào)節(jié)
    浙江電力(2011年5期)2011-05-29 03:48:10
    女警被强在线播放| 精品国产乱子伦一区二区三区| 丝袜美足系列| 少妇粗大呻吟视频| 亚洲av美国av| 在线观看免费日韩欧美大片| 亚洲国产中文字幕在线视频| 亚洲一区中文字幕在线| 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 精品欧美一区二区三区在线| 久久久久国产精品人妻aⅴ院| 神马国产精品三级电影在线观看 | 在线观看免费视频日本深夜| 久热这里只有精品99| 午夜亚洲福利在线播放| 性欧美人与动物交配| 黄色a级毛片大全视频| 国产精品亚洲一级av第二区| 精品乱码久久久久久99久播| 后天国语完整版免费观看| 校园春色视频在线观看| 天天一区二区日本电影三级 | 黄色丝袜av网址大全| 少妇熟女aⅴ在线视频| 国产av在哪里看| av有码第一页| 欧美人与性动交α欧美精品济南到| 成人18禁高潮啪啪吃奶动态图| 九色亚洲精品在线播放| 一进一出好大好爽视频| 久久九九热精品免费| 午夜成年电影在线免费观看| av免费在线观看网站| 一a级毛片在线观看| 国产伦人伦偷精品视频| www.精华液| 亚洲午夜精品一区,二区,三区| 欧美av亚洲av综合av国产av| 亚洲av熟女| 日本vs欧美在线观看视频| 国产1区2区3区精品| 国产成人欧美| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日韩一级在线毛片| 99久久99久久久精品蜜桃| 韩国精品一区二区三区| 久久久久久国产a免费观看| 国产欧美日韩一区二区精品| 丝袜美足系列| 男女午夜视频在线观看| а√天堂www在线а√下载| 女人被狂操c到高潮| 一本综合久久免费| av在线播放免费不卡| 久久精品国产亚洲av高清一级| 性欧美人与动物交配| 欧美绝顶高潮抽搐喷水| 欧美黑人精品巨大| 成人手机av| 亚洲情色 制服丝袜| 成人特级黄色片久久久久久久| 宅男免费午夜| 精品熟女少妇八av免费久了| av电影中文网址| 村上凉子中文字幕在线| 两个人视频免费观看高清| 人妻久久中文字幕网| 看免费成人av毛片| 日日撸夜夜添| 能在线免费观看的黄片| 中文资源天堂在线| 午夜a级毛片| 少妇猛男粗大的猛烈进出视频 | 久久99热6这里只有精品| 中文资源天堂在线| 国产精品av视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 午夜福利成人在线免费观看| 在线观看66精品国产| 日日撸夜夜添| 三级毛片av免费| 国产亚洲精品综合一区在线观看| 亚洲自拍偷在线| 天美传媒精品一区二区| 成熟少妇高潮喷水视频| 全区人妻精品视频| 成人特级av手机在线观看| 三级毛片av免费| 国产精品不卡视频一区二区| 久久精品久久久久久噜噜老黄 | 中文字幕av成人在线电影| 国产一区二区激情短视频| 麻豆国产av国片精品| 欧美潮喷喷水| 色在线成人网| 欧美成人性av电影在线观看| 国内久久婷婷六月综合欲色啪| 亚洲经典国产精华液单| 亚洲专区中文字幕在线| 一本久久中文字幕| 国产伦精品一区二区三区四那| 成人国产一区最新在线观看| 久久久久久久精品吃奶| 免费观看精品视频网站| 久久国产乱子免费精品| 亚洲va在线va天堂va国产| 男女那种视频在线观看| 波多野结衣高清无吗| 国产精品女同一区二区软件 | 少妇丰满av| 一边摸一边抽搐一进一小说| 赤兔流量卡办理| 亚洲真实伦在线观看| 精品免费久久久久久久清纯| 国内久久婷婷六月综合欲色啪| 亚洲国产高清在线一区二区三| АⅤ资源中文在线天堂| 中出人妻视频一区二区| 999久久久精品免费观看国产| 美女高潮喷水抽搐中文字幕| 欧美色欧美亚洲另类二区| 内射极品少妇av片p| 国内揄拍国产精品人妻在线| 午夜精品在线福利| 中文字幕av在线有码专区| 欧美日本亚洲视频在线播放| 国产成人影院久久av| 国内精品久久久久精免费| 大型黄色视频在线免费观看| 69人妻影院| 久久精品国产鲁丝片午夜精品 | 色视频www国产| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品国产三级国产av玫瑰| 1024手机看黄色片| 国产午夜精品论理片| 国内揄拍国产精品人妻在线| 国产麻豆成人av免费视频| 国产伦在线观看视频一区| 男人舔女人下体高潮全视频| 亚洲国产色片| 特级一级黄色大片| 国产精品人妻久久久影院| 国产精品美女特级片免费视频播放器| 熟妇人妻久久中文字幕3abv| 天堂动漫精品| 亚洲欧美精品综合久久99| 国产视频一区二区在线看| 亚洲熟妇中文字幕五十中出| 精品无人区乱码1区二区| 一个人看视频在线观看www免费| 久久精品夜夜夜夜夜久久蜜豆| 99久久精品国产国产毛片| 人人妻人人看人人澡| 黄片wwwwww| 桃红色精品国产亚洲av| 91在线观看av| 精品午夜福利在线看| 国产成人av教育| 亚洲最大成人中文| 日本精品一区二区三区蜜桃| 夜夜夜夜夜久久久久| h日本视频在线播放| 深夜精品福利| 精品午夜福利视频在线观看一区| 国产一级毛片七仙女欲春2| a级毛片a级免费在线| 午夜福利成人在线免费观看| 91久久精品电影网| 国产精品久久久久久av不卡| 人人妻人人看人人澡| 久久精品久久久久久噜噜老黄 | 精华霜和精华液先用哪个| 99riav亚洲国产免费| 欧美潮喷喷水| www.色视频.com| 成人av在线播放网站| 内射极品少妇av片p| 亚洲一级一片aⅴ在线观看| 午夜激情福利司机影院| 国语自产精品视频在线第100页| 亚洲国产欧洲综合997久久,| 色av中文字幕| 又粗又爽又猛毛片免费看| 亚洲欧美日韩卡通动漫| 午夜日韩欧美国产| 精品乱码久久久久久99久播| 亚洲欧美精品综合久久99| 永久网站在线| 69人妻影院| 精品一区二区三区av网在线观看| 无遮挡黄片免费观看| 午夜老司机福利剧场| 天堂av国产一区二区熟女人妻| 一级a爱片免费观看的视频| 久久久久九九精品影院| 人妻制服诱惑在线中文字幕| 久久精品综合一区二区三区| 九色国产91popny在线| 久久久午夜欧美精品| 午夜福利在线观看吧| 小蜜桃在线观看免费完整版高清| 一进一出好大好爽视频| 国产v大片淫在线免费观看| 成年女人看的毛片在线观看| 深爱激情五月婷婷| 级片在线观看| 欧美日本亚洲视频在线播放| 99精品久久久久人妻精品| 国产毛片a区久久久久| 变态另类丝袜制服| 九九爱精品视频在线观看| 欧美一区二区精品小视频在线| 嫩草影院新地址| 人妻丰满熟妇av一区二区三区| 久久国产乱子免费精品| 国产中年淑女户外野战色| 国产黄片美女视频| 国产精品女同一区二区软件 | 18禁在线播放成人免费| 日韩大尺度精品在线看网址| 99久久精品热视频| 极品教师在线视频| 免费看a级黄色片| 免费观看的影片在线观看| 亚洲,欧美,日韩| 亚州av有码| 特大巨黑吊av在线直播| 欧美又色又爽又黄视频| 亚洲最大成人手机在线| 网址你懂的国产日韩在线| 欧美日韩中文字幕国产精品一区二区三区| 十八禁国产超污无遮挡网站| 日本色播在线视频| 热99re8久久精品国产| 在线观看66精品国产| 乱人视频在线观看| 日韩大尺度精品在线看网址| 99久久中文字幕三级久久日本| 美女高潮的动态| 三级国产精品欧美在线观看| 久久午夜福利片| 色5月婷婷丁香| 久久99热这里只有精品18| 国产中年淑女户外野战色| 九九在线视频观看精品| 精品免费久久久久久久清纯| 国产激情偷乱视频一区二区| 国产男人的电影天堂91| 亚洲第一电影网av| 欧美日韩精品成人综合77777| 成人性生交大片免费视频hd| 成年女人看的毛片在线观看| 国产主播在线观看一区二区| 久久99热6这里只有精品| 国产一区二区三区视频了| 狠狠狠狠99中文字幕| 美女被艹到高潮喷水动态| 九九久久精品国产亚洲av麻豆| 精品无人区乱码1区二区| 一区二区三区激情视频| 男女视频在线观看网站免费| 成人高潮视频无遮挡免费网站| 琪琪午夜伦伦电影理论片6080| 日韩欧美国产在线观看| 国产一区二区在线观看日韩| 欧美中文日本在线观看视频| 欧美又色又爽又黄视频| 久久国内精品自在自线图片| .国产精品久久| 免费观看人在逋| 啦啦啦韩国在线观看视频| 热99在线观看视频| 亚洲av电影不卡..在线观看| 国产精品福利在线免费观看| aaaaa片日本免费| 人妻夜夜爽99麻豆av| 久久天躁狠狠躁夜夜2o2o| 欧美+日韩+精品| 久久久久久国产a免费观看| 九色国产91popny在线| 国产不卡一卡二| 中文字幕精品亚洲无线码一区| 国产老妇女一区| 少妇裸体淫交视频免费看高清| 成人欧美大片| 亚洲内射少妇av| 九色国产91popny在线| 日韩精品有码人妻一区| 夜夜夜夜夜久久久久| 一进一出抽搐gif免费好疼| 午夜福利在线观看吧| 日韩国内少妇激情av| 露出奶头的视频| 精品久久久噜噜| 国产高清激情床上av| 国产单亲对白刺激| 日韩欧美精品v在线| 热99re8久久精品国产| 色哟哟哟哟哟哟| 18禁裸乳无遮挡免费网站照片| 看免费成人av毛片| 亚洲欧美日韩卡通动漫| 日本熟妇午夜| 亚洲性夜色夜夜综合| 国产探花在线观看一区二区| 久久人人爽人人爽人人片va| 啦啦啦观看免费观看视频高清| 久久久久久久久大av| 91麻豆av在线| 黄片wwwwww| 性欧美人与动物交配| 99国产极品粉嫩在线观看| 乱系列少妇在线播放| 夜夜看夜夜爽夜夜摸| 亚洲精品影视一区二区三区av| 欧美+日韩+精品| 久久久精品欧美日韩精品| 看片在线看免费视频| 亚洲精品色激情综合| 色尼玛亚洲综合影院| 一a级毛片在线观看| 男人狂女人下面高潮的视频| 国产午夜精品久久久久久一区二区三区 | 丝袜美腿在线中文| 亚洲av五月六月丁香网| 毛片一级片免费看久久久久 | 狂野欧美激情性xxxx在线观看| 日韩大尺度精品在线看网址| 哪里可以看免费的av片| 免费观看精品视频网站| 有码 亚洲区| 联通29元200g的流量卡| 精品久久国产蜜桃| 国产伦精品一区二区三区四那| 欧美精品国产亚洲| 人妻丰满熟妇av一区二区三区| 美女黄网站色视频| 淫妇啪啪啪对白视频| 亚洲一级一片aⅴ在线观看| 黄色丝袜av网址大全| 久久这里只有精品中国| 亚洲图色成人| 亚洲av.av天堂| 九九爱精品视频在线观看| 日韩中字成人| 蜜桃久久精品国产亚洲av| 黄片wwwwww| 欧美成人免费av一区二区三区| 午夜影院日韩av| 夜夜夜夜夜久久久久| av中文乱码字幕在线| 欧美色欧美亚洲另类二区| 波多野结衣高清无吗| 三级毛片av免费| 国产乱人视频| 国产视频内射| 亚洲乱码一区二区免费版| 亚洲人成伊人成综合网2020| 午夜影院日韩av| 色视频www国产| 午夜影院日韩av| 夜夜夜夜夜久久久久| 亚洲精品一区av在线观看| 欧美色欧美亚洲另类二区| av天堂在线播放| 午夜激情欧美在线| 国产精品福利在线免费观看| 国产v大片淫在线免费观看| 精品人妻视频免费看| 欧美中文日本在线观看视频| 少妇丰满av| 蜜桃久久精品国产亚洲av| 成人国产麻豆网| 九九爱精品视频在线观看| 久久亚洲精品不卡| 国产综合懂色| 国产av不卡久久| 中文字幕久久专区| 国产黄a三级三级三级人| 免费观看精品视频网站| АⅤ资源中文在线天堂| 中文字幕av在线有码专区| a级毛片a级免费在线| 亚洲成av人片在线播放无| 俄罗斯特黄特色一大片| 黄色配什么色好看| 嫩草影院入口| 国产熟女欧美一区二区| av女优亚洲男人天堂| 国产高清视频在线播放一区| 热99re8久久精品国产| 亚洲精品在线观看二区| 欧美中文日本在线观看视频| 男女下面进入的视频免费午夜| 校园春色视频在线观看| 高清日韩中文字幕在线| 成人综合一区亚洲| 九九热线精品视视频播放| 国产精品不卡视频一区二区| 欧美成人性av电影在线观看| 国产高清不卡午夜福利| 国内精品宾馆在线| 22中文网久久字幕| 欧美3d第一页| 亚洲美女视频黄频| 亚洲在线自拍视频| 精品久久久久久久久av| 真人做人爱边吃奶动态| 99九九线精品视频在线观看视频| 久久精品国产亚洲av涩爱 | 国产一区二区三区在线臀色熟女| av在线亚洲专区| 嫩草影院新地址| 国产白丝娇喘喷水9色精品| 中国美女看黄片| .国产精品久久| 久久国内精品自在自线图片| 日本免费一区二区三区高清不卡| 日本黄色片子视频| 国产不卡一卡二| 黄色丝袜av网址大全| 深夜精品福利| 在线国产一区二区在线| 两个人视频免费观看高清| 婷婷精品国产亚洲av| 亚洲欧美清纯卡通| 久久精品人妻少妇| 91精品国产九色| 99久久精品热视频| 免费在线观看日本一区| 国产精品爽爽va在线观看网站| 亚洲性久久影院| 69人妻影院| 免费av不卡在线播放| 美女cb高潮喷水在线观看| 中文字幕高清在线视频| 亚洲自偷自拍三级| 国产精品人妻久久久久久| 九九热线精品视视频播放| 成人永久免费在线观看视频| 国产精品久久久久久av不卡| 精品一区二区三区人妻视频| 非洲黑人性xxxx精品又粗又长| 欧美+日韩+精品| 色播亚洲综合网| 全区人妻精品视频| 国产成人av教育| 少妇人妻一区二区三区视频| 日本在线视频免费播放| 哪里可以看免费的av片| 欧美激情国产日韩精品一区| 一边摸一边抽搐一进一小说| 久久热精品热| 成年女人毛片免费观看观看9| 草草在线视频免费看| 真实男女啪啪啪动态图| 国产一区二区亚洲精品在线观看| 日韩国内少妇激情av| 免费人成视频x8x8入口观看| 国产一区二区三区视频了| 在线免费观看不下载黄p国产 | 久久人妻av系列| 亚洲熟妇熟女久久| 搡老熟女国产l中国老女人| 丝袜美腿在线中文| 欧美日韩黄片免| 韩国av在线不卡| 一a级毛片在线观看| 十八禁国产超污无遮挡网站| 中文字幕熟女人妻在线| 白带黄色成豆腐渣| 日韩亚洲欧美综合| 婷婷精品国产亚洲av在线| 如何舔出高潮| 亚州av有码| 亚洲午夜理论影院| 男女视频在线观看网站免费| 国产高潮美女av| 干丝袜人妻中文字幕| 久久99热这里只有精品18| 精品国内亚洲2022精品成人| 搡老妇女老女人老熟妇| 欧美三级亚洲精品| 国内精品久久久久久久电影| 久久这里只有精品中国| 99久久精品一区二区三区| 亚洲av第一区精品v没综合| 中文亚洲av片在线观看爽| 熟女人妻精品中文字幕| www日本黄色视频网| 日韩欧美国产在线观看| 亚洲最大成人av| 国产乱人伦免费视频| 欧美色视频一区免费| 国产午夜精品久久久久久一区二区三区 | 韩国av在线不卡| 久久亚洲精品不卡| 国产av麻豆久久久久久久| 亚洲七黄色美女视频| 免费观看在线日韩| 成人鲁丝片一二三区免费| www.色视频.com| 看片在线看免费视频| 99国产精品一区二区蜜桃av| 三级国产精品欧美在线观看| 国产精品伦人一区二区| 美女高潮喷水抽搐中文字幕| 国产伦在线观看视频一区| 午夜老司机福利剧场| 国产三级中文精品| 日本在线视频免费播放| 欧美性猛交黑人性爽| 欧美激情久久久久久爽电影| 欧美色视频一区免费| 国产 一区精品| 99热这里只有是精品50| 欧美性感艳星| 国产精品亚洲一级av第二区| 久久久久性生活片| 婷婷六月久久综合丁香| 亚洲中文日韩欧美视频| 精品久久久久久久末码| 神马国产精品三级电影在线观看| 高清在线国产一区| 久久久精品欧美日韩精品| 国产高清不卡午夜福利| 免费在线观看成人毛片| 高清日韩中文字幕在线| 色5月婷婷丁香| 男人的好看免费观看在线视频| 人人妻人人澡欧美一区二区| 一本精品99久久精品77| 欧美一区二区亚洲| 日本-黄色视频高清免费观看| 国产视频内射| 丰满人妻一区二区三区视频av| 日日摸夜夜添夜夜添av毛片 | 日日啪夜夜撸| 乱码一卡2卡4卡精品| av在线天堂中文字幕| 全区人妻精品视频| 久久久国产成人精品二区| 亚洲精品色激情综合| 久久精品国产亚洲av涩爱 | 一边摸一边抽搐一进一小说| 日本五十路高清| 尤物成人国产欧美一区二区三区| 我要看日韩黄色一级片| 日本成人三级电影网站| 两人在一起打扑克的视频| 国产精品国产三级国产av玫瑰| 国产精品亚洲美女久久久| 成年人黄色毛片网站| 精华霜和精华液先用哪个| 免费人成在线观看视频色| 免费人成视频x8x8入口观看| 亚洲精品亚洲一区二区| 亚洲男人的天堂狠狠| 亚洲无线在线观看| 九九在线视频观看精品| 性欧美人与动物交配| 亚洲av成人av| 国产白丝娇喘喷水9色精品| 美女 人体艺术 gogo| 久久久成人免费电影| 亚洲精华国产精华液的使用体验 | 免费搜索国产男女视频| 91久久精品电影网| 别揉我奶头 嗯啊视频| 乱人视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 老司机午夜福利在线观看视频| 九色成人免费人妻av| 老熟妇乱子伦视频在线观看| x7x7x7水蜜桃| 亚州av有码| 日韩精品中文字幕看吧| 国产麻豆成人av免费视频| 香蕉av资源在线| 国产精品99久久久久久久久| 亚洲欧美日韩高清专用| 黄片wwwwww| 中文字幕高清在线视频| 美女cb高潮喷水在线观看| av国产免费在线观看| 欧美潮喷喷水| 在线观看午夜福利视频| 精品人妻1区二区| 观看免费一级毛片| 校园人妻丝袜中文字幕| 一区二区三区高清视频在线| 高清在线国产一区| 观看美女的网站| 精品人妻偷拍中文字幕| 黄色丝袜av网址大全| 国产伦人伦偷精品视频| 国产爱豆传媒在线观看| 一本精品99久久精品77| 欧美成人性av电影在线观看| 男人舔奶头视频| 中文字幕人妻熟人妻熟丝袜美| 观看美女的网站| 日韩一本色道免费dvd| 日韩欧美国产在线观看| 亚洲七黄色美女视频| 日韩欧美国产一区二区入口| 亚洲av成人精品一区久久| 毛片一级片免费看久久久久 | 狂野欧美白嫩少妇大欣赏| 色吧在线观看| 婷婷亚洲欧美| 最新在线观看一区二区三区| 麻豆国产97在线/欧美|