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

雙饋風(fēng)電場(chǎng)抑制電網(wǎng)低頻振蕩的自適應(yīng)附加控制策略

2023-09-26 11:58:44劉新宇原紹恒郝正航逯芯妍
關(guān)鍵詞:控制策略發(fā)電機(jī)系統(tǒng)

劉新宇, 王 森, 曾 龍, 原紹恒, 郝正航, 逯芯妍

(1. 華北水利水電大學(xué) 電力學(xué)院,鄭州 450011; 2. 廣東工業(yè)大學(xué) 自動(dòng)化學(xué)院,廣州 510006;3. 貴州大學(xué) 電氣工程學(xué)院,貴陽(yáng) 550025)

現(xiàn)代工業(yè)的迅猛發(fā)展已經(jīng)造成對(duì)煤和石油等化石能源的過度依賴和嚴(yán)重的環(huán)境污染等一系列問題.由于風(fēng)力發(fā)電(簡(jiǎn)稱風(fēng)電)不會(huì)污染環(huán)境且風(fēng)能屬于可再生能源,風(fēng)能的充分利用對(duì)能源結(jié)構(gòu)的調(diào)整和實(shí)現(xiàn)“雙碳”目標(biāo)具有重要意義.近些年,伴隨著大規(guī)模風(fēng)電并網(wǎng)的快速發(fā)展,建設(shè)更大規(guī)模的風(fēng)電場(chǎng)以及遠(yuǎn)距離輸電問題應(yīng)運(yùn)而生[1].由于風(fēng)能具有隨機(jī)、不穩(wěn)定的特點(diǎn),而我國(guó)的用電負(fù)荷中心往往距離風(fēng)電場(chǎng)比較遠(yuǎn),輸電線路負(fù)荷較重,加之風(fēng)電場(chǎng)又處于相對(duì)薄弱的電網(wǎng)末端,所以給風(fēng)電并網(wǎng)系統(tǒng)的穩(wěn)定性造成較大影響[2].

大規(guī)模的風(fēng)電并網(wǎng)、長(zhǎng)距離的輸送線路、風(fēng)能不穩(wěn)定這些問題大大降低了系統(tǒng)的阻尼,因此當(dāng)系統(tǒng)發(fā)生跳閘、三相接地短路或是小擾動(dòng)干擾時(shí),易引起傳輸線路的功率低頻振蕩.迄今為止,已經(jīng)發(fā)現(xiàn)多次電力系統(tǒng)低頻振蕩現(xiàn)象,給電力系統(tǒng)帶來很大影響.在國(guó)內(nèi),隨著中國(guó)新疆維吾爾自治區(qū)北部地區(qū)風(fēng)電裝機(jī)容量的迅速增加,自2014年6月以來,廣域測(cè)量系統(tǒng)已反復(fù)捕獲次同步頻率下的持續(xù)功率振蕩[3].在國(guó)外,德克薩斯電力可靠性委員會(huì)電網(wǎng)中已經(jīng)發(fā)現(xiàn)多次4 Hz的低頻振蕩[4].在此類嚴(yán)重事故中,低頻振蕩甚至?xí)て鸶浇鼫u輪發(fā)電機(jī)的強(qiáng)烈扭轉(zhuǎn)振動(dòng),從而導(dǎo)致發(fā)電廠的所有發(fā)電機(jī)跳閘.因此,有必要增強(qiáng)包含風(fēng)電的混合電力系統(tǒng)的阻尼,以提高系統(tǒng)頻率穩(wěn)定性.

為保證風(fēng)電并網(wǎng)系統(tǒng)穩(wěn)定安全運(yùn)行,很多國(guó)家已經(jīng)出臺(tái)相關(guān)的政策規(guī)定在風(fēng)電機(jī)組中必須帶有可以增強(qiáng)系統(tǒng)阻尼的附加環(huán)節(jié).文獻(xiàn)[5]中提出基于區(qū)域極點(diǎn)配置的風(fēng)電系統(tǒng)弱阻尼低頻振蕩抑制控制策略,采用線性矩陣不等式的區(qū)域極點(diǎn)配置方法設(shè)計(jì)反饋控制器,將系統(tǒng)中的弱阻尼振蕩模式配置到所控制區(qū)域,提高系統(tǒng)的阻尼比.文獻(xiàn)[6]中以系統(tǒng)的頻率偏差值作為輸入,根據(jù)不同的風(fēng)速設(shè)置槳距角,設(shè)計(jì)基于變系數(shù)的下垂控制調(diào)節(jié)風(fēng)力發(fā)電機(jī)的有功出力,為系統(tǒng)注入正阻尼.文獻(xiàn)[7]中根據(jù)兩區(qū)域電力系統(tǒng)模型,選擇合適的輸入輸出變量,對(duì)輸入函數(shù)和狀態(tài)函數(shù)求導(dǎo)得到反饋控制律,提出基于非線性控制算法的附加阻尼控制策略,可以有效地處理阻尼控制器中的非線性環(huán)節(jié).文獻(xiàn)[8]中將雙饋風(fēng)力發(fā)電機(jī)(DFIG)和電力系統(tǒng)靜態(tài)穩(wěn)定器(PSS)放在一起,建立系統(tǒng)的狀態(tài)方程,定義系統(tǒng)振蕩模式對(duì)傳遞函數(shù)的靈敏度,然后運(yùn)用定義的靈敏度優(yōu)化控制參數(shù),該控制策略能有效抑制系統(tǒng)低頻振蕩.文獻(xiàn)[9]中分析含DFIG的風(fēng)電系統(tǒng)可能會(huì)引起低頻振蕩的原因,對(duì)比基于復(fù)轉(zhuǎn)矩系數(shù)法的4種控制策略,為開發(fā)最優(yōu)控制策略以抑制DFIG系統(tǒng)在不同條件下的軸系振動(dòng)提供了參考.文獻(xiàn)[10]中提出一種有源-無源協(xié)調(diào)雙通道功率調(diào)制阻尼控制器,該控制器可以較好地發(fā)揮DFIG轉(zhuǎn)子側(cè)控制器的阻尼控制能力.文獻(xiàn)[11]中提出基于電力系統(tǒng)穩(wěn)定器和功率振蕩阻尼器的魯棒協(xié)調(diào)控制策略改善低頻振蕩阻尼,該方法使用改進(jìn)的基于特征值的目標(biāo)函數(shù)實(shí)現(xiàn),并用灰狼優(yōu)化器進(jìn)行優(yōu)化.文獻(xiàn)[12]中提出基于阻耗系數(shù)的附加阻尼模糊控制策略,通過控制轉(zhuǎn)子側(cè)換流器調(diào)節(jié)DFIG無功功率來增強(qiáng)系統(tǒng)阻尼,進(jìn)而抑制系統(tǒng)低頻振蕩.文獻(xiàn)[13]中提出一種基于主動(dòng)抗擾控制的有功無功附加阻尼控制器,能起到很好的抑制負(fù)阻尼的作用.文獻(xiàn)[14]中提出一種基于自抗擾控制器的廣域阻尼控制器協(xié)調(diào)優(yōu)化策略,采用人工蜂群算法對(duì)自抗擾控制器和廣域阻尼控制器進(jìn)行協(xié)調(diào)優(yōu)化,提高系統(tǒng)阻尼水平.文獻(xiàn)[15]中設(shè)計(jì)基于頻率反饋的大規(guī)模柔性負(fù)荷控制策略和統(tǒng)一潮流控制器串聯(lián)側(cè)功率反饋附加控制策略,最后運(yùn)用遺傳算法優(yōu)化協(xié)同控制參數(shù),得到最優(yōu)控制策略.文獻(xiàn)[16]中分析有功功率環(huán)節(jié)和無功功率環(huán)節(jié)附加阻尼控制的差異,提出一種基于bang-bang調(diào)制的混合快速阻尼控制策略和基于二階滑模的雙饋風(fēng)電機(jī)組魯棒阻尼控制策略,對(duì)電力系統(tǒng)區(qū)間振蕩有很好的抑制效果.文獻(xiàn)[17]中利用狀態(tài)反饋的方法將系統(tǒng)輸入、輸出作線性化處理,將非線性系統(tǒng)線性化,最后提出基于超螺旋算法的阻尼控制器,該控制方法不依賴于系統(tǒng)的詳細(xì)建模參數(shù),具有良好的工程應(yīng)用價(jià)值.

綜上所述,國(guó)內(nèi)外學(xué)者針對(duì)電力系統(tǒng)區(qū)域間功率低頻振蕩已經(jīng)進(jìn)行了深入的研究,也提出大量的附加阻尼控制器設(shè)計(jì)方法.然而,大多數(shù)設(shè)計(jì)方案將系統(tǒng)中聯(lián)絡(luò)線功率或發(fā)電機(jī)功角偏差等參數(shù)作為輸入,忽略系統(tǒng)平衡點(diǎn)變化后的動(dòng)態(tài)特性對(duì)低頻振蕩的影響,在系統(tǒng)平衡點(diǎn)附近線性化的基礎(chǔ)上設(shè)計(jì)類似于PSS結(jié)構(gòu)的附加阻尼控制器,而很少有研究將DFIG自身轉(zhuǎn)子電壓和磁鏈的關(guān)系應(yīng)用到阻尼控制器設(shè)計(jì)上.本文對(duì)DFIG的動(dòng)態(tài)頻率響應(yīng)和轉(zhuǎn)子勵(lì)磁電壓與磁鏈之間的關(guān)系進(jìn)行研究,結(jié)合終端滑模變結(jié)構(gòu)控制的思想,設(shè)計(jì)基于轉(zhuǎn)子磁鏈偏差的快速終端滑模附加阻尼控制器,并證明該系統(tǒng)的收斂性和穩(wěn)定性.該控制器將磁鏈的偏差作為輸入,輸出一個(gè)附加信號(hào)到轉(zhuǎn)子側(cè)功率控制環(huán)節(jié),進(jìn)而控制轉(zhuǎn)子側(cè)換流器,調(diào)節(jié)DFIG出力,抑制系統(tǒng)低頻振蕩現(xiàn)象.最后,通過仿真驗(yàn)證本文所提控制方法對(duì)系統(tǒng)低頻功率振蕩的影響.

1 DFIG響應(yīng)阻尼調(diào)節(jié)能力分析

DFIG系統(tǒng)在三相靜止坐標(biāo)系中的模型階是一個(gè)多變量、高階次、強(qiáng)耦合、非線性的復(fù)雜時(shí)變系統(tǒng).為了便于計(jì)算和控制,運(yùn)用坐標(biāo)變換將其轉(zhuǎn)換到兩相同步旋轉(zhuǎn)d-q坐標(biāo)系下進(jìn)行分析[18].忽略定子繞組壓降和電磁暫態(tài)變化過程,將定子電壓定向于d軸上,DFIG定子輸出有功功率(Ps)、無功功率(Qs)和轉(zhuǎn)子磁鏈方程可表示為

(1)

式中:Ls、Lr、Lm分別代表定子自感、轉(zhuǎn)子自感和定轉(zhuǎn)子互感;ω1為同步轉(zhuǎn)速;ψqs為d-q坐標(biāo)系下定子q軸磁鏈分量;ψdr、ψqr分別為d-q坐標(biāo)系下轉(zhuǎn)子磁鏈分量.

(2)

式中:Rr為轉(zhuǎn)子電阻;ωs=ω1-ωr為發(fā)電機(jī)轉(zhuǎn)差角速度,其中ωr為發(fā)電機(jī)轉(zhuǎn)子轉(zhuǎn)速;udr、uqr分別為轉(zhuǎn)子d-q軸電壓分量;Ug為電網(wǎng)側(cè)電壓.

圖1 兩區(qū)域電力系統(tǒng)模型Fig.1 Two-region power system model

2 終端滑模附加控制器設(shè)計(jì)

由上述分析可知,調(diào)節(jié)DFIG輸出有功功率可以提高系統(tǒng)的阻尼水平,因此可以通過調(diào)節(jié)DFIG的有功功率改善系統(tǒng)阻尼水平,進(jìn)而抑制系統(tǒng)區(qū)域間的低頻振蕩.

2.1 系統(tǒng)滑模面與控制律設(shè)計(jì)

由式(2)可得,DFIG運(yùn)行過程中,轉(zhuǎn)子外加電壓和轉(zhuǎn)子磁鏈之間的關(guān)系如圖2所示.圖中:s為DFIG轉(zhuǎn)差率.DFIG在發(fā)電狀態(tài)時(shí),將轉(zhuǎn)子磁鏈分量ψdr、ψqr作為系統(tǒng)的受控輸出,則轉(zhuǎn)子勵(lì)磁電壓udr、uqr為系統(tǒng)的控制輸入.為達(dá)到系統(tǒng)快速控制的目的,將終端滑模變結(jié)構(gòu)控制方法運(yùn)用到DFIG轉(zhuǎn)子磁鏈控制方案之中,設(shè)計(jì)轉(zhuǎn)子磁鏈控制器,使系統(tǒng)磁鏈跟蹤期望磁鏈,確保系統(tǒng)在指定時(shí)間內(nèi)迅速達(dá)到穩(wěn)定狀態(tài).

圖2 DFIG轉(zhuǎn)子外加電壓與磁鏈之間的關(guān)系圖Fig.2 Relationship between applied voltage and flux linkage of DFIG rotor

(3)

基于以上分析,設(shè)計(jì)如下滑模面:

(4)

式中:參數(shù)ρ1>0;p1、q1為正奇數(shù),且0

式(4)對(duì)時(shí)間求導(dǎo)得:

(5)

將式(3)代入式(5)得:

(6)

則轉(zhuǎn)子d軸磁鏈采用的控制律可設(shè)計(jì)為

udr=udreq+udrn

(7)

式中:udrep、udrn分別為等效控制律和切換控制律,且滿足

(8)

2.2 系統(tǒng)收斂性與穩(wěn)定性分析

(1) 系統(tǒng)收斂性證明.

證明當(dāng)Sd=0時(shí),系統(tǒng)位于滑模面上,此時(shí)有

(9)

求解可得:

(10)

則DFIG的d軸磁鏈狀態(tài)可在有限時(shí)間Δtd內(nèi)快速收斂到0,證畢.

(2) 系統(tǒng)穩(wěn)定性證明.

證明取李雅普諾夫函數(shù):

V(t)=Sd(t)/2

(11)

求導(dǎo)可得:

(12)

將式(6)代入式(12)可得:

(13)

所以V(t)≤0,系統(tǒng)在李雅普諾夫穩(wěn)定性理論上全局漸進(jìn)穩(wěn)定,證畢.

同理可得轉(zhuǎn)子q軸磁鏈控制律為

uqr=uqreq+uqrn

(14)

(15)

式中:參數(shù)ρ2>0;p2、q2為正奇數(shù),且0

圖3 機(jī)側(cè)控制器模型Fig.3 Control model of machine side

系統(tǒng)中使能環(huán)節(jié)表達(dá)式如下:

(16)

3 系統(tǒng)離線仿真與分析

根據(jù)前文推導(dǎo)的快速終端滑模附加控制的原理設(shè)計(jì)附加阻尼控制器,在MATLAB/Simulink中建立如圖4所示的四機(jī)兩區(qū)域模型并進(jìn)行仿真分析.圖中:G2、G3和G4為同步發(fā)電機(jī),L1、L2分別代表區(qū)域1和區(qū)域2負(fù)荷.

仿真環(huán)境設(shè)置:設(shè)t=15 s時(shí),兩區(qū)域聯(lián)絡(luò)線發(fā)生0.1 s的三相短路,此時(shí)系統(tǒng)由于阻尼不足,在聯(lián)絡(luò)線上發(fā)生功率振蕩;t=15.1 s時(shí)故障消除,系統(tǒng)恢復(fù)正常運(yùn)行狀態(tài).觀察無附加阻尼控制下和本文提出的附加阻尼控制策略下系統(tǒng)中各項(xiàng)參數(shù)的響應(yīng)曲線.

圖5為同步發(fā)電機(jī)轉(zhuǎn)速振蕩曲線.圖中:ω為轉(zhuǎn)速.對(duì)比有、無附加控制的響應(yīng)曲線可知,在考慮附加阻尼控制情況下的同步發(fā)電機(jī)轉(zhuǎn)子轉(zhuǎn)速響應(yīng)曲線的振幅更小,更快趨于穩(wěn)定,具有更好的阻尼特性.由圖6所示的同步機(jī)G1和G3轉(zhuǎn)子角速度的差值(Δω13)振蕩曲線可知,在考慮附加阻尼的情況下,G1和G3轉(zhuǎn)子角速度偏差振蕩的振幅和時(shí)間均減小.說明所提出的附加控制策略對(duì)系統(tǒng)同步發(fā)電機(jī)的振蕩有一定抑制作用.

圖5 同步發(fā)電機(jī)轉(zhuǎn)速變化曲線Fig.5 Speed curves of synchronous generator

圖6 同步發(fā)電機(jī)G1和G3轉(zhuǎn)子角速度偏差振蕩曲線Fig.6 Angular velocity deviation oscillation mode of synchronous generator G1 and G3

圖7為DFIG轉(zhuǎn)子勵(lì)磁電流(Ir)的變化曲線.對(duì)比圖7(a)和7(b)可以看出,當(dāng)系統(tǒng)發(fā)生低頻振蕩時(shí),在有附加阻尼控制情況下,故障期間DFIG轉(zhuǎn)子電流增大并有一段時(shí)間的振蕩響應(yīng),使DFIG及時(shí)發(fā)出有功功率支撐電網(wǎng);在沒有附加阻尼控制情況下,故障期間DFIG轉(zhuǎn)子電流變化較小,不能很好地響應(yīng)系統(tǒng)振蕩.圖8為兩區(qū)域聯(lián)絡(luò)線功率(P12)變化曲線,從圖中可以看出在考慮了附加控制的情況下,聯(lián)絡(luò)線功率振蕩的振幅明顯減小,并且更快地趨于穩(wěn)定.

圖7 雙饋發(fā)電機(jī)轉(zhuǎn)子電流變化曲線Fig.7 Rotor current curves of DFIG

圖8 聯(lián)絡(luò)線功率振蕩曲線Fig.8 Power oscillation curves of transmission line

圖9為DFIG定子電流(Is)變化曲線.由圖可知,在無附加控制下,故障期間DFIG定子電流幾乎沒有變化.而在采用附加控制時(shí),故障期間定子會(huì)短時(shí)振蕩現(xiàn)象,但其電流值仍在安全范圍之內(nèi),說明所提的附加阻尼控制策略在及時(shí)調(diào)節(jié)有功功率輸出的同時(shí),并不會(huì)對(duì)雙饋風(fēng)電系統(tǒng)帶來負(fù)面影響.圖10為DFIG有功功率(Pw)響應(yīng)曲線.從圖中可以看出,在沒有考慮附加控制情況下,系統(tǒng)發(fā)生低頻振蕩之后,DFIG有功功率并沒有響應(yīng)系統(tǒng)振蕩,而在考慮本文所設(shè)計(jì)的附加阻尼控制情況下,DFIG有功功率輸出能快速做出動(dòng)態(tài)響應(yīng)以提高系統(tǒng)的阻尼水平來抑制低頻振蕩.

圖9 雙饋發(fā)電機(jī)定子電流響應(yīng)變化Fig.9 Stator current curve of DFIG

圖10 DFIG有功輸出振蕩曲線Fig.10 Response curves of DFIG active power

4 系統(tǒng)實(shí)時(shí)仿真實(shí)驗(yàn)驗(yàn)證與分析

為了充分驗(yàn)證所提附加控制策略在工程應(yīng)用上的可行性,開發(fā)基于RTDS的風(fēng)電并網(wǎng)四機(jī)兩區(qū)域?qū)崟r(shí)仿真系統(tǒng),主要仿真參數(shù)和離線仿真參數(shù)相同.

在RTDS平臺(tái)實(shí)驗(yàn)中,系統(tǒng)運(yùn)行至t=1.2 s時(shí),在兩區(qū)域聯(lián)絡(luò)線上設(shè)置三相瞬時(shí)短路故障,傳輸線路出現(xiàn)功率低頻振蕩現(xiàn)象;t=1.4 s時(shí),去除故障;監(jiān)測(cè)DFIG定子電流、轉(zhuǎn)子電流和有功輸出響應(yīng)情況以及系統(tǒng)振蕩頻率在是否考慮附加阻尼兩種情況下變化曲線,結(jié)果如圖11~14所示.

圖11 系統(tǒng)頻率變化曲線Fig.11 Frequency change curves of system

圖11為系統(tǒng)頻率(f)響應(yīng)曲線.從圖中可以看出,在系統(tǒng)發(fā)生低頻振蕩故障期間,無附加阻尼控制情況下的系統(tǒng)頻率衰減較慢,振蕩持續(xù)時(shí)間相對(duì)較長(zhǎng).而采用所提附加控制之后,系統(tǒng)頻率振蕩曲線衰減較快,且振蕩的幅值更低,系統(tǒng)整體表現(xiàn)出更好的阻尼特性.圖12為DFIG有功出力響應(yīng)曲線.由圖可見,當(dāng)聯(lián)絡(luò)線發(fā)生功率振蕩時(shí),如果考慮附加控制,DFIG就可以及時(shí)輸出更多有功功率響應(yīng)其振蕩,為系統(tǒng)提供更多阻尼.

圖12 DFIG有功出力變化曲線Fig.12 Active output oscillation curves of DFIG

圖13為轉(zhuǎn)子勵(lì)磁電流實(shí)時(shí)仿真曲線.圖中:I為電流;Ia、Ib、Ic分別為abc三相電流.對(duì)比實(shí)驗(yàn)結(jié)果可知:在采用附加阻尼控制策略的情況下,在故障期間,轉(zhuǎn)子勵(lì)磁電流增大,以便使雙饋發(fā)電機(jī)及時(shí)發(fā)出更多有功功率;而無附加阻尼控制下,轉(zhuǎn)子勵(lì)磁電流在故障期間幾乎沒有變化.而且從圖14可以看出,定子電流實(shí)時(shí)仿真結(jié)果與離線仿真時(shí)的動(dòng)態(tài)響應(yīng)曲線一致,均能夠保證DFIG的安全運(yùn)行.

圖13 DFIG轉(zhuǎn)子電流變化曲線Fig.13 Rotor current curves of DFIG

圖14 DFIG定子電流變化曲線Fig.14 Stator current variation curves of DFIG

5 結(jié)語(yǔ)

研究風(fēng)電功率調(diào)節(jié)對(duì)系統(tǒng)低頻振蕩抑制問題,在DFIG的機(jī)側(cè)電流控制環(huán)上引入以增加電力系統(tǒng)阻尼為目標(biāo)、基于李雅普諾夫穩(wěn)定性理論的快速終端滑模附加阻尼控制器.利用MATLAB仿真測(cè)試與實(shí)驗(yàn)驗(yàn)證,在電力系統(tǒng)出現(xiàn)低頻振蕩故障時(shí),所提出的快速終端滑模控制不僅可以快速調(diào)節(jié)DFIG有功功率,提高系統(tǒng)的阻尼水平,有效抑制低頻振蕩;而且在系統(tǒng)響應(yīng)過程中,風(fēng)力發(fā)電機(jī)的有功功率響應(yīng)速度快、變化平滑,定子電流振蕩持續(xù)時(shí)間很短且振蕩幅值在安全范圍內(nèi),不會(huì)影響風(fēng)電場(chǎng)安全運(yùn)行,具有一定的應(yīng)用參考價(jià)值.

猜你喜歡
控制策略發(fā)電機(jī)系統(tǒng)
Smartflower POP 一體式光伏系統(tǒng)
WJ-700無人機(jī)系統(tǒng)
考慮虛擬慣性的VSC-MTDC改進(jìn)下垂控制策略
能源工程(2020年6期)2021-01-26 00:55:22
ZC系列無人機(jī)遙感系統(tǒng)
工程造價(jià)控制策略
山東冶金(2019年3期)2019-07-10 00:54:04
現(xiàn)代企業(yè)會(huì)計(jì)的內(nèi)部控制策略探討
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
大型發(fā)電機(jī)勵(lì)磁用旋轉(zhuǎn)變換器的開發(fā)和應(yīng)用
隨身攜帶的小發(fā)電機(jī)
軍事文摘(2016年16期)2016-09-13 06:15:49
容錯(cuò)逆變器直接轉(zhuǎn)矩控制策略
主站蜘蛛池模板: 无码区日韩专区免费系列 | 国产精品lululu在线观看| 在线观看免费人成视频色快速| 四虎永久在线视频| a毛片免费观看| 丁香婷婷在线视频| 在线视频一区二区三区不卡| 午夜综合网| 亚洲Va中文字幕久久一区| 国产一在线观看| 九月婷婷亚洲综合在线| 国产凹凸一区在线观看视频| 在线日本国产成人免费的| 色偷偷一区二区三区| 精品国产一二三区| 亚洲无码视频喷水| 色久综合在线| 亚洲区欧美区| 亚洲色图欧美一区| 久久毛片网| 国产成人一区免费观看| 国产农村精品一级毛片视频| 亚洲开心婷婷中文字幕| 91 九色视频丝袜| 人妻一区二区三区无码精品一区| 国产99欧美精品久久精品久久| 国产成人亚洲无吗淙合青草| 国产又爽又黄无遮挡免费观看| av大片在线无码免费| 91色在线视频| 日本欧美午夜| 日本手机在线视频| 在线观看亚洲国产| 亚洲AV人人澡人人双人| 国产成人精品无码一区二 | 国产精品yjizz视频网一二区| 欧美精品亚洲精品日韩专| 午夜精品久久久久久久2023| 久久a毛片| 午夜精品一区二区蜜桃| 成人免费午间影院在线观看| 国产男女免费完整版视频| 国产人成在线视频| 永久免费av网站可以直接看的| 精品伊人久久久久7777人| 亚洲一区二区三区国产精华液| 91精品啪在线观看国产91九色| 久久99国产乱子伦精品免| 日韩区欧美区| 美女无遮挡拍拍拍免费视频| 国产精品一线天| 波多野结衣一区二区三区四区视频 | 波多野结衣爽到高潮漏水大喷| 国产精品美女自慰喷水| 欧美中文字幕一区| 五月激情婷婷综合| 亚洲自偷自拍另类小说| 美女高潮全身流白浆福利区| 国产三级精品三级在线观看| 欧美精品不卡| 欧美在线视频a| 婷婷六月在线| 欧美精品成人一区二区视频一| 国产精品网址你懂的| 国产99热| 午夜少妇精品视频小电影| 91www在线观看| 久久国产毛片| 亚洲欧美日韩成人高清在线一区| 国产成人精品视频一区视频二区| 国产在线日本| 四虎影视库国产精品一区| 色哟哟国产精品一区二区| 亚洲天堂成人在线观看| 国产三级成人| 九色在线观看视频| 亚洲精品图区| 国产门事件在线| 无码中文字幕精品推荐| 午夜视频免费一区二区在线看| 日韩欧美国产综合| 中国丰满人妻无码束缚啪啪|