999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

二階LADRC在風電并網(wǎng)逆變器網(wǎng)側直流母線電壓控制中的運用

2024-03-21 03:17:22周雪松葛建鵬馬幼捷
電測與儀表 2024年3期
關鍵詞:系統(tǒng)

周雪松,葛建鵬,馬幼捷

(1. 天津理工大學 天津市復雜系統(tǒng)控制理論與應用重點實驗室,天津 300384;2.天津理工大學 電氣電子工程學院, 天津 300384)

0 引 言

在風力發(fā)電系統(tǒng)中,并網(wǎng)變流器通常運用電壓電流雙閉環(huán)控制來穩(wěn)定并網(wǎng)網(wǎng)側直流母線電壓[1]。隨著現(xiàn)代控制技術的發(fā)展和新形式下更加嚴格的并網(wǎng)要求,雖然電壓電流的雙閉環(huán)PI控制方法在實際工程中得到了廣泛的應用并取得了一定的效果,但這種基于誤差消除誤差的被動控制方法與擾動的影響相比具有一定的滯后性,也可能導致系統(tǒng)產生過大的超調和振蕩次數(shù)的增加,進而影響控制性能,使得控制器輸出的能量不能滿足控制要求[2-3]。

為了改進傳統(tǒng)控制方法的不足,提高直流側母線電壓的控制性能,在文獻[4]中將滑動濾波算法與傳統(tǒng)PI控制器相結合,雖然該方法提高了并網(wǎng)電流的正弦飽和和直流側母線電壓的抗干擾性能,但控制器的參數(shù)整定比較復雜,不利于工程應用。在文獻[5]中,將滑膜控制法運用到電壓外環(huán)上,而電流內環(huán)在PI的控制的基礎上,增加了電流預測的方法,雖然可以有效減弱網(wǎng)側諧波影響,讓母線電壓更加穩(wěn)定,但由于控制器結構復雜,硬件實現(xiàn)較困難。

為了提高風力發(fā)電并網(wǎng)逆變器直流母線電壓為控制的控制性能,本文設計了一種線性二階自抗擾控制系統(tǒng)[6-7]。一方面,系統(tǒng)還可以利用線性微分跟蹤器(linear differential tracker,LTD)提高系統(tǒng)的跟蹤性能;另一方面,通過引入線性擴張狀態(tài)觀測器(linear extended state observer, LESO),增加一個狀態(tài)變量,實現(xiàn)了對擾動進行觀測與干擾補償,避免了積分環(huán)節(jié)的使用所帶來的響應慢、易振蕩、控制量積分飽和等負面影響。該本文以風力發(fā)電并網(wǎng)逆變器直流母線電壓為控制對象。首先介紹了風力發(fā)電逆變器的建模過程和傳統(tǒng)LADRC控制器的基本結構。其次,加入TD增強對反饋信號的跟蹤效果,使得進入LESO的信號保真性大,干擾小,精度高;接著在LESO中加入全擾動差分信號,增加其觀測帶寬。通過合理配置極點來參數(shù)整定。再次,利用頻域分析法對比本文所設計LADRC與傳統(tǒng)PI控制下系統(tǒng)的跟蹤、穩(wěn)定、抗擾性能。最后,通過3.6 MW風力發(fā)電機組全真模擬平臺驗證了該LADRC的控制性能。

1 風電并網(wǎng)逆變器的數(shù)學模型及其傳統(tǒng)PI控制策略

1.1 風電并網(wǎng)逆變器的數(shù)學模型

直驅式發(fā)電機(permanent magnet synchronous generator, PMSG)通過全功率整流器并網(wǎng),其機側變流器通過對風機轉速控制,實現(xiàn)最大風能利用;網(wǎng)側變流器用來穩(wěn)定網(wǎng)側直流母線電壓及并網(wǎng)的電能質量[8-10]。風能轉換系統(tǒng)如圖1所示。

圖1 直驅式永磁同步發(fā)電機原理圖

典型的三相電壓型PWM風力發(fā)電并網(wǎng)逆變器的電路拓撲如圖2所示。

圖2 并網(wǎng)逆變器拓撲圖

由圖2可以得到abc坐標系下網(wǎng)側變流器的數(shù)學模型。

(1)

由于在三相靜止坐標系下電流存在耦合,且在都是時變量,不利于控制器的設計,因此進行坐標變換,得到[11]:

(2)

式中ω是鎖相環(huán)角頻率;egd和egq是三相電壓在d、q軸的投影分量;id和iq是網(wǎng)側電流在d、q軸的投影分量;Sd和Sq是開關函數(shù)的d,q軸分量。

并網(wǎng)變流器的輸出電壓為:

(3)

則在d-q坐標系下的并網(wǎng)變流器的數(shù)學模型為:

(4)

其中,ωLid、ωLiq為解耦項。

由式(4)知,abc坐標系下的交流量經(jīng)坐標變換后轉換成為d-q坐標系中的直流量,從而簡化了數(shù)學模型,雖然存在耦合項,但經(jīng)過變換后使系統(tǒng)控制器的設計更加便捷。

1.2 雙閉環(huán)PI控制分析

傳統(tǒng)的并網(wǎng)變流器的控制方法是電壓外環(huán)和電流內環(huán)均為PI控制。當系統(tǒng)為三相對稱時,電網(wǎng)電壓投影到d軸時,若在功率因數(shù)為1時并網(wǎng),則其無功為零,egd=E,egq=0,在穩(wěn)態(tài)時,idiq均為直流,則導數(shù)項均為零[10]。

根據(jù)式(4),得到穩(wěn)態(tài)控制方程,即:

(5)

式中E為三相電壓的相電壓峰值。

在同步旋轉坐標系中,變流器輸出到電網(wǎng)的有功功率和無功功率可以表示為:

(6)

三相電壓型PWM并網(wǎng)逆變器在傳統(tǒng)PI控制方式下的控制框圖如圖3所示。

圖3 雙閉環(huán)PI控制框圖

設ud、uq的控制方程為:

(7)

式中KiP和KiI為內環(huán)的比例系數(shù)和積分系數(shù),id_ref和iq_ref分別是內環(huán)d軸和q軸的給定值。將方程(7)代入方程(4),得到方程(8)為:

(8)

由式(8)可知,電網(wǎng)側d軸輸出電流控制模型和q軸輸出電流控制模型均只包含各自的分量,從而實現(xiàn)了系統(tǒng)的解耦控制。這樣,在設計系統(tǒng)的PI控制器時,單獨控制d軸電流就可控制輸出有功率,控制q軸電流就可控制無功功率。同時利用零d軸電流控制策略,當定子電流的d軸分量為零時,此時電樞反應與轉子磁鏈垂直,從而避免了永磁體的退磁現(xiàn)象,進而定子電流全部用來產生電磁轉矩,就可實現(xiàn)發(fā)電機的最大轉矩控制,減少機側PWM變流器各部分的功率損耗,為此設計出來的控制器將更加精確和簡單,從而使系統(tǒng)的控制變得更加精確和簡單。之后設計控制器,根據(jù)“一階最優(yōu)”原則選擇電流內環(huán)控制器參數(shù),得到最優(yōu)階躍響應;根據(jù)“模式最優(yōu)”原則設置電壓外環(huán)控制器參數(shù),得到最優(yōu)調節(jié)性能,保證系統(tǒng)穩(wěn)定。從圖3可以看出,ud和uq作為控制變量輸出。經(jīng)過坐標變換后,將PWM策略接口連接起來,得到最終的開關函數(shù)來控制網(wǎng)側逆變器的開關量,從而實現(xiàn)對直流側母線電壓的控制。

2 二階LADRC控制原理

所設計的二階LADRC控制器如圖4所示[12]。

圖4 二階LADRC示意圖

2.1 二階線性跟蹤微分器的設計

將給定信號v(t)經(jīng)過LTD輸出兩個信號:r1(t)和r2(t)。其中,r1(t)跟蹤v(t),又r2(t)=r1(t),因此r2(t)可近似為v(t)的導數(shù)。LTD的作用:(1)為安排過渡過程;(2)消除超調與快速性之間的矛盾;(3)消除或減弱噪聲的放大效應,將D控制運用至控制器中[12,15]。

給定信號v(t)經(jīng)過二階LTD在頻域下得式(9)。

(9)

若時間常數(shù)τ1、τ2很接近常數(shù)τ,那么傳遞函數(shù)為:

(10)

(11)

其狀態(tài)變量實現(xiàn)為:

(12)

當參數(shù)r適當大時具有很好的高階微分功能。

2.2 三階線性擴展狀態(tài)觀測器的設計

一般二階系統(tǒng)為:

(13)

式中y,u分別為控制器的輸出和輸入;w為擾動。而a1、a2以及w均未知,b部分已知(已知部分記為b0),則式(13)可寫為:

(14)

將系統(tǒng)內外的總擾動擴展為系統(tǒng)的新狀態(tài)變量x3=f。系統(tǒng)的狀態(tài)方程為:

(15)

設計LESO為:

(16)

2.3 干擾補償方法

傳統(tǒng)的PI控制器通過積分消除誤差[13]。然而,增加積分環(huán)節(jié)導致其相位滯后,從而降低系統(tǒng)的穩(wěn)定性,增長系統(tǒng)的響應時間。二階LADRC可以利用三階LESO對廣義擾動進行實時估計,并對估計值進行補償,避免了積分環(huán)節(jié)帶來的副作用。系統(tǒng)控制律取式(17)[13]。

(17)

式中u0是線性誤差反饋控制率的輸出。將上述公式代入式可得出:

(18)

由式(18)可以看出,當f≈z3時,系統(tǒng)可以變成一個積分串聯(lián)結構,從而簡化被控對象,提高了控制性能。這就是干擾補償?shù)姆椒ā?/p>

2.4 線性誤差反饋率的設計

線性誤差反饋控制率的輸出u0為:

u0=kp(r1-z1)+kD(r2-z2)

(19)

式中r1為給定v經(jīng)TD的跟蹤信號;r2為v的近似導數(shù);kp和kD是控制器參數(shù)。根據(jù)式(18)和式(19),系統(tǒng)的閉環(huán)傳遞函數(shù)為:

G(s)=

(20)

二階LADRC的結構可以通過上述分析得到,如圖5所示。

圖5 二階LADRC的結構圖

2.5 二階LADRC控制器的參數(shù)整定

二階LADRC中二階LTD、三階LESO和LSEF的參數(shù)設置是否合理,將直接影響整個控制器的性能。每一部分需要設定的參數(shù)有:二階LTD的參數(shù)r、三階LESO的觀測器增益β1、β2、β3和LSEF的控制器參數(shù)kP、kD。

2.5.1 TD與LSEF的參數(shù)整定

根據(jù)式(20),將系統(tǒng)特征方程的極點放在同一個位置,即:

s4+(kD+2r)s3+(kP+2rkD+r2)s2+(2rkP+kDr2)s+kPr2=(s+r)2(s2+kDs+kP)=(s+ωc)4

(21)

因此,TD與LSEF中的控制器參數(shù)分別表示為:

(22)

式中ωc是控制器的帶寬。

由此可見,TD與PD控制器中唯一需要設置的參數(shù)是ωc,ωc越大,系統(tǒng)的輸出響應越快,動態(tài)過程時間越短。然而,在實際工程參數(shù)設置過程中,雖然ωc值越大,TD的跟蹤效果越好,但考慮到PD控制器的負載越大,導致系統(tǒng)對噪聲的抑制能力減弱,甚至會導致失穩(wěn)。因此,實際工程中的參數(shù)整定需要平衡系統(tǒng)的快速性和穩(wěn)定性[14-16]。

2.5.2 三階LESO的參數(shù)整定

經(jīng)過參數(shù)化,將所有的觀測器極點都配置在同一個位置上,即取觀測器的增益矩陣為:

(23)

式中ω0為觀測器帶寬。因此ω0是LESO中唯一要設置的參數(shù)。ω0越大,LESO的帶寬越寬,跟蹤狀態(tài)變量的精度越高,因此控制器的控制質量越好。但在實際參數(shù)設置中,ω0太大也會導致測量噪聲的放大,不利于系統(tǒng)的控制。因此,ω0在實際工程中不宜過大。觀測噪聲對系統(tǒng)的影響應綜合考慮,適當調整參數(shù)大小。

在文中的參數(shù)設置過程中,ωc保持不變,ω0逐漸增大,直到噪聲的影響滿足系統(tǒng)的要求。然后逐漸增大ωc,當噪聲影響不可承受時減小ω0,然后增大ωc,以達到預期的控制效果。最后,選擇外電壓環(huán)的二階LADRC參數(shù)為ωc=3 000和ω0=500。

3 LADRC的穩(wěn)定性、跟蹤性及抗擾性分析

LADRC具有很強的穩(wěn)定性、跟蹤性及抗擾能力。LADRC的設計是設計ωc和ω0的參數(shù),然后分析ωc和ω0對控制性能的影響[17]。

3.1 LADRC的穩(wěn)定性分析

結合式(20)與式(22),可知系統(tǒng)電壓外環(huán)的閉環(huán)傳遞函數(shù)為:

(24)

則系統(tǒng)的特征方程為:

(25)

根據(jù)赫爾維茨穩(wěn)定性判據(jù),系統(tǒng)穩(wěn)定的充要條件是在a0>0的前提下,各階主子式Δi均大于零,即:

由以上公式可知所設計的控制器具有良好的穩(wěn)定性。

3.2 LADRC的跟蹤特性分析

3.2.1 LTD的跟蹤特性分析

二階跟蹤微分器具有比經(jīng)典微分器更低的噪聲放大增益。雖然在相同濾波器時間常數(shù)下,二階微分器的相位滯后比經(jīng)典微分器略大一些,但其增加采樣頻率來抑制噪聲放大的特性,使得其具有比經(jīng)典微分器更好的微分預測性能。根據(jù)式(11),利用頻率Bode圖分析了二階跟蹤微分器的跟蹤特性,如圖6所示。

圖6 對給定信號的抗擾能力

將信號v(t)作為系統(tǒng)的輸入,相當于信號經(jīng)過一個二階慣性環(huán)節(jié),其阻尼比ξ≥1,近似一個低通濾波器,r越大,即T越小,其帶寬就越大,跟蹤效果就越好,對高頻噪聲的敏感性也低。

由于系統(tǒng)是臨界阻尼或過阻尼的,因此二階TD系統(tǒng)的帶寬可以理解為轉折頻率,其中參數(shù)r越大,即帶寬越大,系統(tǒng)更快到達設定。因此參數(shù)r可以當作決定系統(tǒng)跟蹤速度的速度因子。

由參數(shù)整定可知r=ωc,即調整ωc大小即可,雖然ωc值越大,TD的跟蹤效果越好,但考慮到PD控制器的負載越大,導致系統(tǒng)抗擾能力減弱,進而導致系統(tǒng)不穩(wěn)定。

3.2.2 LESO的跟蹤特性分析

三階LESO是影響二階LADRC控制器性能的核心部件,因此采用經(jīng)典的頻域分析方法分析了三階LESO的收斂性以及總擾動對三階LESO性能的影響。將式(16)轉換成傳遞函數(shù)的形式并代入式(23),則z1、z2和z3的傳遞函數(shù)為:

(26)

(27)

(28)

由上述方程可以看出,LESO能夠實現(xiàn)對系統(tǒng)狀態(tài)變量的跟蹤和對擾動的無偏估計。

3.3 LADRC的抗擾特性分析

LADRC具有很強的抗干擾能力。LADRC的設計是設計ωc和ω0的參數(shù),然后分析ωc和ω0的對控制性能的影響。

3.3.1 總擾動對二階LTD性能的影響分析

根據(jù)式(11),利用頻率Bode圖,分析了二階跟蹤微分器對擾動微分信號的抗擾能力分析,如圖7所示。

圖7 對微分信號的抗擾能力

3.3.2 總擾動對三階LESO性能的影響分析

根據(jù)第3.1節(jié)中定義的系統(tǒng)各狀態(tài)變量的三階LESO估計誤差為:

(29)

三階LESO估計誤差的狀態(tài)方程可由方程(16)減去方程(15)得到,得出方程為:

(30)

在對上述公式進行拉氏變換后,得到各狀態(tài)變量估計誤差對總擾動的傳遞函數(shù)為:

(31)

從方程(31)可以看出,總擾動是影響三階LESO估計誤差的唯一因素。利用頻率Bode圖,詳細分析了總擾動對三階LESO各估計誤差的影響領域通過改變ω0比較每個狀態(tài)變量的估計誤差的頻域曲線如圖8所示。

圖8 各狀態(tài)變量估計誤差的頻域曲線

從圖8所示的頻域特性曲線可以看出,隨著觀測器帶寬ω0的增大,總擾動對LESO狀態(tài)估計誤差的影響逐漸減小。其中,中低頻段總干擾引起的估計誤差值比較明顯。但對高頻段的估計誤差沒有影響;隨著觀測器帶寬ω0的增加,各誤差傳遞函數(shù)的相位裕度沒有明顯變化。頻率在ω0附近的總擾動信號對三階LESO各態(tài)的估計誤差影響最大。對于頻率大于觀測器帶寬頻率的總干擾,三階LESO不能準確估計,導致LESO的估計性能受到帶寬的影響。

4 基于二階LADRC的電壓外環(huán)控制系統(tǒng)設計

根據(jù)以上分析,建立基于二階LADRC的風電并網(wǎng)逆變器電壓外環(huán)控制器,建立了網(wǎng)側變流器直流側數(shù)學模型,如式(32)所示[17-18]。

(32)

通過對上述方程的兩邊進行求導取其微分,得到:

(33)

式(33)的狀態(tài)空間描述形式,如式(34)所示。

(34)

根據(jù)式(16)和式(23),電壓外環(huán)的三階LESO可由式(35)得到:

(35)

LTD與LSEF的設計如式(36)所示。

(36)

其中,udc_ref是直流母線的給定,id_ref是電流內環(huán)的d軸給定。本文設計的控制框圖如圖9所示。

圖9 雙閉環(huán)控制圖

5 仿真分析

為驗證控制器的有效性,先將控制器模型離散化,并在3.6 MW級全真風場模擬實驗平臺中進行對比實驗。系統(tǒng)的采樣時間設置為6×10-6s。控制器參數(shù)見表1,系統(tǒng)主要參數(shù)如表2所示。

表1 控制器參數(shù)

表2 系統(tǒng)主要參數(shù)

為了驗證設計的基于二階LADRC控制器的電壓外環(huán)比傳統(tǒng)PI控制器具有更好的抗干擾特性,將兩種控制方式下并網(wǎng)逆變器直流側母線電壓的動態(tài)響應速度和穩(wěn)態(tài)性能與系統(tǒng)穩(wěn)態(tài)運行時低穿、機側加載和機側減載的擾動所造成實驗進行了比較。

5.1 電網(wǎng)電壓低穿時的比較

5.1.1 電網(wǎng)電壓對稱跌落10%

將故障引起的電網(wǎng)電壓對稱跌落量設為10%,故障開始時間為1.5 s,故障結束時間為2.0 s,系統(tǒng)實驗時間為0.5 s,其它條件相同,比較兩種控制方式下直流側母線電壓波形。電網(wǎng)電壓波形如圖10(a)所示,兩種控制方式下的直流側母線電壓波形如圖10(b)所示。

圖10 電網(wǎng)平衡跌落10%并網(wǎng)點與直流母線電壓波形圖

5.1.2 電網(wǎng)電壓平衡跌落25%

將故障引起的電網(wǎng)電壓對稱跌落量設為25%。電網(wǎng)電壓波形如圖11(a)所示,兩種控制方式下的直流側母線電壓波形如圖11(b)所示。

圖11 電網(wǎng)平衡跌落25%并網(wǎng)點與直流母線電壓波形圖

5.2 風機機側加載50%

將故障引起的風機機側加載量設為50%,設置加載時間為1.0 s,其它條件相同,比較兩種控制方式下直流側母線電壓波形。兩種控制方式下的直流側母線電壓波形如圖12所示。

圖12 風機機側加載50%直流母線電壓波形圖

在電機突然加載50%(即機側有功增加50%)時,并網(wǎng)點母線的電壓會因為機側功率增大導致電壓升高,PI控制下的并網(wǎng)點的母線電壓超調量為6%,調整時間為50 ms; 在LADRC控制下,母線電壓僅超調量為2%,調整時間為25 ms,其控制效果是優(yōu)于PI控制的。

5.3 風機機側甩載100%

將故障引起的風機機側甩載量設為100%,設置加載時間為1.0 s,其它條件相同,比較兩種控制方式下直流側母線電壓波形。兩種控制方式下的直流側母線電壓波形如圖13所示。

圖13 風機機側甩載100%直流母線電壓波形圖

在電機突然甩載100%(即機側有功減少100%)時,并網(wǎng)點母線的電壓會因為機側功率降低導致電壓降低,PI控制下的并網(wǎng)點的母線電壓最大跌落至0.85 p.u.,其調整時間約為30 ms; 在LADRC控制下,母線電壓最大跌落至0.95 p.u,調整時間為20 ms,其控制效果仍然是優(yōu)于PI控制的。

圖10、11為并網(wǎng)點的電壓降至0.9 p.u.及0.75p.u.時,二階LADRC與傳統(tǒng)PI控制母線電壓的仿真對比。在1.5 s時發(fā)生電壓跌落故障,電壓跌落深度為額定電壓的10%,持續(xù)0.5 s后恢復正常,如圖9所示。若在此期間風力發(fā)電機繼續(xù)保持正常運行,變頻器仍保持控制不變。以低穿25%時進入誤差為±0.05%為穩(wěn)態(tài)舉例,如表3所示。

表3 低穿25%下兩個控制方式的性能指標對比

由表3可知,電網(wǎng)電壓降至0.75 p.u.時,風機側的有功率將繼續(xù)向直流母線傳遞,此時母線電壓突增,經(jīng)過一段時間后將重新進入穩(wěn)定狀態(tài),采用LADRC控制時其性能指標均優(yōu)于在傳統(tǒng)PI控制下的母線電壓。

由此可知,在相同的運行工況中,結合LADRC控制策略會大幅降低直流測電壓波動的峰值、減少系統(tǒng)過渡時間,增強系統(tǒng)的跟蹤性和抗擾性,從而可以更好地提升電壓故障穿越的能力。

6 結束語

本文以風電逆變器母線電壓為控制對象,針對電網(wǎng)電壓波動與機側擾動對母線電壓存在影響這一情況,為了提高風電并網(wǎng)逆變器直流側母線電壓的動態(tài)響應速度和穩(wěn)定性,設計了基于二階LADRC的電壓外環(huán)控制器。這項工作的主要貢獻是:

1)利用跟蹤微分器(LTD)對輸入信號安排過渡過程并提取其微分信號;

2)采用三階LESO實時估計和跟蹤系統(tǒng)狀態(tài)的變化,不僅得到各個狀態(tài)變量的估計,還能對擾動進行觀測和估計,提高了對直流側母線電壓的精確控制;

3)采用頻域分析方法分析了LADRC的穩(wěn)定性、跟蹤性和抗擾性能;

4)在全真風場模擬實驗平臺上建立了3.6 MW級直驅永磁風力發(fā)電系統(tǒng)。與傳統(tǒng)的PI控制效果相比,二階LADRC比PI控制器更能處理系統(tǒng)中的各種干擾,具有較強的適應性與魯棒性。

實驗結果表明,基于二階LADRC電壓外環(huán)控制器的風電并網(wǎng)逆變器控制方法可以有效地提高直流側母線電壓的動態(tài)響應速度和穩(wěn)態(tài)性能,為復雜情況下提高風力發(fā)電并網(wǎng)質量提供了有效的參考和幫助。

猜你喜歡
系統(tǒng)
Smartflower POP 一體式光伏系統(tǒng)
WJ-700無人機系統(tǒng)
ZC系列無人機遙感系統(tǒng)
北京測繪(2020年12期)2020-12-29 01:33:58
基于PowerPC+FPGA顯示系統(tǒng)
基于UG的發(fā)射箱自動化虛擬裝配系統(tǒng)開發(fā)
半沸制皂系統(tǒng)(下)
FAO系統(tǒng)特有功能分析及互聯(lián)互通探討
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
一德系統(tǒng) 德行天下
PLC在多段調速系統(tǒng)中的應用
主站蜘蛛池模板: 成人国产小视频| 亚洲成人77777| 国产99精品久久| 久久久久亚洲Av片无码观看| 92午夜福利影院一区二区三区| 日韩精品无码一级毛片免费| 日韩人妻精品一区| 亚洲色欲色欲www在线观看| 国产人免费人成免费视频| 国产黄色免费看| 91色老久久精品偷偷蜜臀| 蜜桃视频一区| 精品伊人久久久香线蕉| 91欧美在线| 亚洲色中色| 亚洲久悠悠色悠在线播放| 99re精彩视频| 精品国产一区91在线| 久久精品人妻中文视频| 国产亚洲欧美日本一二三本道| 色综合五月婷婷| 日韩欧美中文| 国产第一页免费浮力影院| 国产AV毛片| 国产成人a在线观看视频| 小说区 亚洲 自拍 另类| 国产在线自乱拍播放| 欧美精品一区二区三区中文字幕| 中文字幕中文字字幕码一二区| 狠狠干综合| 无码高潮喷水在线观看| 欧美成在线视频| 无码有码中文字幕| 午夜国产小视频| 久久这里只有精品免费| 国产免费一级精品视频| 欧美日韩中文国产| 国内精品视频在线| 亚洲天堂免费| 国产网友愉拍精品视频| 中美日韩在线网免费毛片视频 | 国内精品久久久久鸭| 国产微拍一区二区三区四区| 99无码中文字幕视频| 国产www网站| 91午夜福利在线观看| 亚洲天堂网2014| 全部毛片免费看| 亚洲视频影院| 囯产av无码片毛片一级| 国产精品亚洲日韩AⅤ在线观看| 天天色综合4| 四虎影视国产精品| 国产sm重味一区二区三区| 麻豆国产在线观看一区二区| 久热精品免费| 国产精品手机视频一区二区| av手机版在线播放| 国产美女精品人人做人人爽| 欧美成人免费午夜全| 欧美午夜小视频| 美女被躁出白浆视频播放| 国产精品hd在线播放| 成人亚洲天堂| 国产91丝袜在线播放动漫| 精品无码一区二区三区在线视频| 片在线无码观看| 久久精品66| 无码有码中文字幕| 最新亚洲人成无码网站欣赏网| 欧美伦理一区| 色AV色 综合网站| 亚洲最大看欧美片网站地址| 无码一区中文字幕| 嫩草影院在线观看精品视频| 奇米影视狠狠精品7777| 成人亚洲国产| 日韩精品毛片| 97超爽成人免费视频在线播放| 欧美一级色视频| 亚洲清纯自偷自拍另类专区| 中文字幕av无码不卡免费|