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

    基于交替方向乘子法的分布式負(fù)荷頻率控制策略

    2021-11-17 11:50:06李玲芳陳義宣許巖文福拴
    電力建設(shè) 2021年11期
    關(guān)鍵詞:區(qū)域間增益控制策略

    李玲芳,陳義宣,許巖,文福拴

    (1. 云南電網(wǎng)有限責(zé)任公司電網(wǎng)規(guī)劃建設(shè)研究中心,昆明市 650011;2. 浙江大學(xué)電氣工程學(xué)院,杭州市 310027)

    0 引 言

    負(fù)荷頻率控制 (load frequency control, LFC) 指根據(jù)系統(tǒng)頻率變化和聯(lián)絡(luò)線傳輸功率變化調(diào)節(jié)發(fā)電機(jī)的有功出力,從而使電力系統(tǒng)頻率和聯(lián)絡(luò)線交換功率維持在設(shè)定的允許范圍內(nèi)[1]。由于很難精確預(yù)測負(fù)荷需求,需要實(shí)時(shí)調(diào)整發(fā)電機(jī)有功出力以跟蹤負(fù)荷需求的動(dòng)態(tài)變化,從而維持系統(tǒng)頻率在給定的允許范圍之內(nèi)。導(dǎo)致系統(tǒng)實(shí)際頻率偏離額定頻率。LFC對(duì)抑制系統(tǒng)頻率偏差和區(qū)域間聯(lián)絡(luò)線上的功率偏差從而維持電力系統(tǒng)安全穩(wěn)定運(yùn)行具有重要意義。

    電力系統(tǒng)頻率調(diào)節(jié)包括一次調(diào)頻和二次調(diào)頻[2]。在二次調(diào)頻中,通常采用比例-積分(proportional-integral,PI)控制。不過,PI控制的動(dòng)態(tài)性能高度依賴其增益,而高增益一方面可以提高系統(tǒng)的響應(yīng)速度,另一方面會(huì)惡化系統(tǒng)動(dòng)態(tài)行為從而可能導(dǎo)致系統(tǒng)頻率振蕩和失穩(wěn)。在此背景下,本文研究如何設(shè)計(jì)分布式最優(yōu)LFC策略來進(jìn)行二次調(diào)頻以實(shí)現(xiàn)頻率偏差為0且系統(tǒng)動(dòng)態(tài)性能最優(yōu)。

    為提升LFC的性能,現(xiàn)有文獻(xiàn)中采用了一些先進(jìn)控制方法設(shè)計(jì)負(fù)荷頻率控制器。文獻(xiàn)[3-5]采用魯棒控制策略抑制負(fù)荷變化導(dǎo)致的系統(tǒng)頻率偏差。文獻(xiàn)[6-7]研究了計(jì)及通信延遲情況下負(fù)荷頻率控制器的設(shè)計(jì)問題。文獻(xiàn)[8-9]采用自適應(yīng)控制策略提升系統(tǒng)運(yùn)行點(diǎn)發(fā)生變化時(shí)的系統(tǒng)控制性能。文獻(xiàn)[10]采用最優(yōu)控制策略抑制系統(tǒng)頻率偏差和聯(lián)絡(luò)線功率偏差。此外,模型預(yù)測控制[11]、自抗擾控制[12]、滑??刂芠13-14]和事件驅(qū)動(dòng)控制[15]也被用于設(shè)計(jì)負(fù)荷頻率控制器以改善控制性能。不過,這些方法應(yīng)用于實(shí)際多區(qū)域電力系統(tǒng)負(fù)荷頻率控制器設(shè)計(jì)時(shí),參數(shù)整定比較困難?;诰€性矩陣不等式的魯棒控制策略的控制器參數(shù)整定相對(duì)容易,但這種方法通常需要對(duì)系統(tǒng)模型進(jìn)行降階預(yù)處理,從而可能導(dǎo)致系統(tǒng)部分動(dòng)態(tài)特性的丟失。

    也有一些文獻(xiàn)采用啟發(fā)式優(yōu)化算法[16-21]設(shè)計(jì)負(fù)荷頻率控制器。文獻(xiàn)[22]采用粒子群優(yōu)化算法整定控制器參數(shù)。文獻(xiàn)[23]采用灰狼優(yōu)化算法確定分?jǐn)?shù)階控制器的參數(shù)。不過從理論上無法保證這些啟發(fā)式優(yōu)化算法能夠獲得控制器參數(shù)的最優(yōu)解,尤其是對(duì)于實(shí)際多區(qū)域電力系統(tǒng),需要優(yōu)化的控制器參數(shù)較多,啟發(fā)式優(yōu)化算法就未必能求得最優(yōu)解。

    在上述背景下,本文提出一種基于分布式優(yōu)化策略的負(fù)荷頻率控制器設(shè)計(jì)方法,以優(yōu)化系統(tǒng)頻率控制性能,并提高控制器參數(shù)整定的計(jì)算效率。在所構(gòu)建的優(yōu)化模型中,目標(biāo)函數(shù)包含兩項(xiàng),分別表征控制性能和控制器增益矩陣稀疏度。采用交替方向乘子法(alternating direction method of multiplier, ADMM)求解該優(yōu)化問題。所提方法具有下述優(yōu)點(diǎn):1) 可以獲得與集中控制策略相當(dāng)?shù)念l率控制性能,且通訊復(fù)雜度較低;2) 利用ADMM的特點(diǎn)可把優(yōu)化目標(biāo)分解為可用解析方式求解的2個(gè)子優(yōu)化問題,進(jìn)而提高參數(shù)整定效率,且適用于實(shí)際多區(qū)域電力系統(tǒng)。最后,以三區(qū)域電力系統(tǒng)模型為例對(duì)所提方法的可行性和有效性進(jìn)行驗(yàn)證。

    1 負(fù)荷頻率控制模型

    多區(qū)域互聯(lián)電力系統(tǒng)中的每個(gè)區(qū)域都設(shè)有一個(gè)調(diào)度控制中心,系統(tǒng)頻率和區(qū)域間聯(lián)絡(luò)線上的交換功率均在其監(jiān)控之中。目前大多數(shù)汽輪發(fā)電機(jī)和水輪發(fā)電機(jī)都裝有調(diào)速裝置,可以對(duì)系統(tǒng)頻率變化做出響應(yīng)。

    圖1為多區(qū)域互聯(lián)電力系統(tǒng)中第i個(gè)區(qū)域的負(fù)荷頻率控制系統(tǒng)框圖。盡管電力系統(tǒng)的穩(wěn)態(tài)和動(dòng)態(tài)模型都是非線性的,但負(fù)荷頻率控制是針對(duì)小擾動(dòng)的,可以近似采用線性化模型[24]。圖1中的KBi和KEi為傳統(tǒng)的PI控制器參數(shù)。

    圖1所示的第i個(gè)區(qū)域的模型可用下述微分方程組描述[24]:

    圖1 互聯(lián)電力系統(tǒng)中第i個(gè)區(qū)域的模型Fig.1 The model of the i-th area in an interconnected power system

    (1)

    (2)

    (3)

    (4)

    式中:Δfi(t)、ΔPgi(t)、ΔXgi(t)和Δδi(t)分別表示t

    時(shí)刻的頻率偏差、有功輸出調(diào)整量、調(diào)速器閥門位置調(diào)整量和轉(zhuǎn)子角偏差;N表示所研究的互聯(lián)電力系統(tǒng)中包含的區(qū)域數(shù)目;ΔPdi(t)表示t時(shí)刻的負(fù)荷擾動(dòng);TGi、TTi和TPi分別表示調(diào)速器、汽輪機(jī)和電力系統(tǒng)的時(shí)間常數(shù);KPi和Ri分別表示電力系統(tǒng)的增益和速度調(diào)節(jié)系數(shù);Ksij表示區(qū)域i和區(qū)域j之間的連接增益,如果這2個(gè)區(qū)域間沒有功率交換則Ksij為0;ui(t)表示區(qū)域i的輔助控制輸入。

    為便于表述,可將微分方程組式(1)—(4)改寫成矩陣形式:

    (5)

    式中:Xi(t)和Xj(t)分別表示區(qū)域i和區(qū)域j在t時(shí)刻的狀態(tài)向量,Xi(t)=[Δfi(t) ΔPgi(t) ΔXgi(t) Δδi(t)]T;Ui(t)表示控制器在t時(shí)刻的輸入向量;

    2 分布式最優(yōu)負(fù)荷頻率控制器設(shè)計(jì)

    根據(jù)式(5),包含n個(gè)區(qū)域的互聯(lián)電力系統(tǒng)的負(fù)荷頻率控制問題可描述為:

    (6)

    在互聯(lián)電力系統(tǒng)中,一個(gè)區(qū)域內(nèi)的負(fù)荷擾動(dòng)可能會(huì)對(duì)另一個(gè)區(qū)域造成明顯影響,嚴(yán)重時(shí)甚至?xí)a(chǎn)生級(jí)聯(lián)故障(cascade failure)。為盡可能降低這種區(qū)域間的耦合影響,可以采用最優(yōu)控制策略實(shí)現(xiàn)負(fù)荷頻率控制[10]。雖然采用集中最優(yōu)控制策略設(shè)計(jì)控制器增益K可以實(shí)現(xiàn)良好的負(fù)荷頻率控制性能,但其對(duì)區(qū)域間的信息實(shí)時(shí)交互要求較高。另一方面,描述多區(qū)域互聯(lián)電力系統(tǒng)動(dòng)態(tài)過程的微分方程組相當(dāng)復(fù)雜,這對(duì)采用集中優(yōu)化策略實(shí)現(xiàn)負(fù)荷頻率控制帶來了很大困難。

    本文所研究的負(fù)荷頻率控制器設(shè)計(jì)方法可以用下述優(yōu)化問題描述:

    minJ(K)+γg(K)

    (7)

    s.t.(A-B2K)TP+P(A-B2K)=
    -(Q-KTRK)

    (8)

    式(8)所描述的等式約束為系統(tǒng)穩(wěn)定條件,可通過動(dòng)態(tài)系統(tǒng)穩(wěn)定性理論分析得到[26]。

    式(7)中的優(yōu)化目標(biāo)包含兩項(xiàng):J和g。J表示傳統(tǒng)的控制性能指標(biāo)。僅以指標(biāo)J作為優(yōu)化目標(biāo)時(shí),可實(shí)現(xiàn)對(duì)頻率控制性能的優(yōu)化。但這種集中控制策略需要反饋的信號(hào)數(shù)量較多,因此增加了區(qū)域間信息實(shí)時(shí)交互的復(fù)雜度。g為控制器增益矩陣稀疏度指標(biāo),g的值越小,控制器增益矩陣的非零元素越少,需要反饋的信號(hào)數(shù)量也就越少。本文采用J和g之和作為優(yōu)化目標(biāo),以兼顧最優(yōu)控制性能和降低區(qū)域間信息交互的復(fù)雜程度這兩方面的期望。

    由于式(7)中的優(yōu)化目標(biāo)同時(shí)包含了指標(biāo)J和指標(biāo)g,在求解由式(7)和(8)所描述的優(yōu)化模型時(shí),采用傳統(tǒng)基于梯度的求解方法一般只能求得局部最優(yōu)解。采用ADMM算法可以有效求解這個(gè)問題。ADMM算法將要求解的優(yōu)化問題分解為2個(gè)子優(yōu)化問題,然后對(duì)這2個(gè)子優(yōu)化問題交替迭代求解,最終求得原始優(yōu)化問題的解。理論上已證明采用ADMM算法可求得控制器參數(shù)的最優(yōu)解[27]。

    采用ADMM算法求解式(7)和(8)所描述的優(yōu)化模型的步驟如下:

    步驟1:將式(7)和(8)所描述優(yōu)化模型分解為2個(gè)子優(yōu)化問題。

    采用g(G)表示式(7)中指標(biāo)g(K),則由式(7)和(8)所描述的優(yōu)化模型可分解為2個(gè)子問題:針對(duì)K優(yōu)化系統(tǒng)的頻率偏差指標(biāo)J(·)

    (9)

    和針對(duì)G優(yōu)化控制器結(jié)構(gòu)指標(biāo)g(·)

    (10)

    Lp(K,G,Λ)=J(K)+γg(G)+

    tr[ΛT(K-G)]+ (ρ/2)‖K-G‖2

    (11)

    式中:Λ表示拉格朗日乘子;ρ表示正的標(biāo)量。

    步驟2:交替求解子優(yōu)化問題式(9)和式(10)。

    基于ADMM的方法架構(gòu),通過對(duì)子優(yōu)化問題式(9)和式(10)交替求解,最終可求得分布式最優(yōu)負(fù)荷頻率控制策略;這種方法的計(jì)算效率較高,適用于決策變量較多的優(yōu)化問題。采用ADMM算法求解這2個(gè)子優(yōu)化問題的過程如圖2所示,其具體實(shí)現(xiàn)流程可參看文獻(xiàn)[27],因篇幅所限這里不再贅述。

    圖2 基于ADMM的分布式最優(yōu)負(fù)荷頻率控制器設(shè)計(jì)過程Fig.2 Design procedure of the ADMM-based distributed optimal load frequency controller

    3 算例與結(jié)果

    3.1 算例參數(shù)

    采用一個(gè)三區(qū)域電力系統(tǒng)算例來說明所提分布式負(fù)荷頻率控制策略的有效性。IEEE 10機(jī)39節(jié)點(diǎn)電力系統(tǒng)即新英格蘭測試系統(tǒng)[28]可劃分為3個(gè)區(qū)域,連接母線3和4、14和15、16和17為3個(gè)區(qū)域之間的聯(lián)絡(luò)線。根據(jù)文獻(xiàn)[29],經(jīng)過抽象建模之后的三區(qū)域控制系統(tǒng)連接框圖如圖3所示。u1、u2和u3為輸入的控制信號(hào),由所提出的控制策略優(yōu)化確定。區(qū)域1和區(qū)域3采用圖1所示的汽輪機(jī)調(diào)速系統(tǒng)模型。區(qū)域2采用水輪機(jī)調(diào)速器系統(tǒng),其數(shù)學(xué)模型如下詳述[28]。

    圖3 三區(qū)域電力系統(tǒng)模型框圖Fig.3 Block diagram of a three-area power system model

    (12)

    (13)

    (14)

    (15)

    (16)

    (17)

    式中:ΔPg2(t)表示水輪機(jī)有功輸出調(diào)整量;ΔXz2(t)表示閥門位置調(diào)整量;ΔXm2(t)表示引導(dǎo)閥位置調(diào)整量;ΔXy2(t)表示閥門動(dòng)作速度補(bǔ)償調(diào)整量;Tw2表示水流慣性時(shí)間常數(shù);TG2表示調(diào)速器時(shí)間常數(shù);Tr2表示軟反饋時(shí)間常數(shù);rt2表示軟反饋增益;rp2表示下垂系數(shù);Ks2為增益。

    本文在MATLAB2016b編程環(huán)境下進(jìn)行仿真研究。參數(shù)設(shè)置為:KP1=120 Hz/pu,KP2=112.5 Hz/pu,KP3=115 Hz/pu;TP1=20 s,TP2=25 s,TP3=20 s;TG1=0.08 s,TG2=0.20 s,TG3=0.07 s;TT1=0.3 s,Tw2=1 s,TT3=0.35 s;R1=2.4 Hz/pu,rp2=0.04 pu,R3=2.5 Hz/pu;Tr2=6 s;rt2=0.5 pu;Ks2=0.2 pu;Ks12=Ks13=Ks21=Ks23=Ks31=Ks32=0.5 pu/Hz。

    3.2 控制器增益矩陣

    通過圖2所示算法,可得到控制器增益矩陣K,其結(jié)構(gòu)如圖4所示。在圖4中,藍(lán)色的點(diǎn)表示該處矩陣元素不為0。從圖4中可以看出,當(dāng)γ=0.000 10時(shí),K的所有元素都不為0,即控制結(jié)構(gòu)為集中模式;這意味著對(duì)每個(gè)區(qū)域的控制器設(shè)計(jì)而言都需要用其相鄰區(qū)域的系統(tǒng)狀態(tài)信息來計(jì)算式(6)中的控制信號(hào)U。隨著γ值的增加,K中的零元素越來越多;當(dāng)γ=0.390 690時(shí),K中只有9個(gè)非零元素,此時(shí)控制器結(jié)構(gòu)呈現(xiàn)分布式特征,且區(qū)域間需要交互的系統(tǒng)狀態(tài)信息最少。

    圖4 控制器增益矩陣K的結(jié)構(gòu)Fig.4 Structure of the gain matrix K

    3.3 控制性能及與現(xiàn)有方法比較

    本節(jié)通過算例對(duì)所提分布式負(fù)荷頻率控制策略和傳統(tǒng)基于PI的負(fù)荷頻率控制策略進(jìn)行對(duì)比分析,以說明所提方法的有效性。給定3個(gè)區(qū)域所受到的負(fù)荷擾動(dòng)分別為ΔPd1=0.010 pu、ΔPd2=0.020 pu和ΔPd3=0.015 pu。

    當(dāng)采用基于PI的負(fù)荷頻率控制策略時(shí),獲得的系統(tǒng)頻率響應(yīng)偏差和區(qū)域間聯(lián)絡(luò)線功率響應(yīng)偏差如圖5所示??梢钥闯觯@2個(gè)偏差在15 s時(shí)均仍存在微小振蕩。

    采用集中式負(fù)荷頻率控制策略時(shí),控制器增益矩陣結(jié)構(gòu)如圖4(a)所示。系統(tǒng)的頻率響應(yīng)偏差和聯(lián)絡(luò)線功率響應(yīng)偏差如圖6所示。從圖中可以看出,3個(gè)區(qū)域的頻率偏差在15 s時(shí)幾乎為0。與圖5對(duì)比可知,采用集中式負(fù)荷頻率控制策略時(shí)的頻率響應(yīng)超調(diào)量更小。比較圖5和圖6中的功率響應(yīng)曲線可得類似結(jié)論。因此,采用集中式負(fù)荷頻率控制能有效抑制負(fù)荷擾動(dòng)造成的系統(tǒng)頻率偏差和聯(lián)絡(luò)線功率偏差,并改善頻率控制性能。

    圖5 基于PI的負(fù)荷頻率控制響應(yīng)曲線Fig.5 Response curves of the three-area interconnected power system with the PI-based LFC strategy under load disturbances

    圖6 集中式負(fù)荷頻率控制響應(yīng)曲線Fig.6 Response curves of the three-area interconnected power system with the centralized LFC strategy under load disturbances

    采用所提分布式負(fù)荷頻率控制策略時(shí),控制器增益矩陣結(jié)構(gòu)如圖4(d)所示。系統(tǒng)的頻率響應(yīng)偏差和聯(lián)絡(luò)線功率響應(yīng)偏差如圖7所示??梢钥闯?個(gè)區(qū)域的頻率偏差在15 s時(shí)幾乎為0。與圖6對(duì)比可知,采用所提分布式負(fù)荷頻率控制策略時(shí)的頻率控制性能與集中式負(fù)荷頻率控制一致,只是區(qū)域1和區(qū)域3的頻率偏差響應(yīng)的超調(diào)量在2.5 s附近存在微小差異。比較圖7和圖6中的功率響應(yīng)曲線可得類似結(jié)論。因此,由3.2節(jié)的分析可知,采用圖4(d)所示控制器結(jié)構(gòu)的信息交換復(fù)雜程度低于采用圖4(a)所示控制器結(jié)構(gòu)。因此,所提分布式負(fù)荷頻率控制策略具有較低的通信復(fù)雜度。

    圖7 分布式負(fù)荷頻率控制響應(yīng)曲線Fig.7 Response curves of the three-area interconnected power system with the distributed LFC strategy under load disturbances

    4 結(jié) 語

    本文提出基于交替方向乘子法的分布式負(fù)荷頻率控制策略,以優(yōu)化頻率控制性能,且降低對(duì)區(qū)域間信息的交互需要從而提高計(jì)算效率。在所提優(yōu)化方法中,目標(biāo)函數(shù)包含二次最優(yōu)控制性能指標(biāo)和控制器增益矩陣稀疏化指標(biāo);該優(yōu)化問題被轉(zhuǎn)化為2個(gè)子優(yōu)化問題并采用ADMM有效求取控制器增益矩陣,從而降低了控制器實(shí)現(xiàn)的復(fù)雜程度。采用含三區(qū)域的互聯(lián)電力系統(tǒng)進(jìn)行仿真驗(yàn)證,計(jì)算結(jié)果表明針對(duì)負(fù)荷擾動(dòng)和時(shí)變參數(shù),所提方法能把系統(tǒng)頻率偏差和區(qū)域聯(lián)絡(luò)線功率偏差控制到0。

    后續(xù)研究工作將針對(duì)區(qū)域間實(shí)時(shí)交換的信息數(shù)據(jù)不完整或有錯(cuò)誤時(shí),對(duì)所提方法進(jìn)行測試并基于相關(guān)結(jié)果對(duì)所提分布式負(fù)荷頻率控制設(shè)計(jì)方法進(jìn)行改進(jìn)和發(fā)展。

    猜你喜歡
    區(qū)域間增益控制策略
    基于增益調(diào)度與光滑切換的傾轉(zhuǎn)旋翼機(jī)最優(yōu)控制
    考慮虛擬慣性的VSC-MTDC改進(jìn)下垂控制策略
    能源工程(2020年6期)2021-01-26 00:55:22
    基于單片機(jī)的程控增益放大器設(shè)計(jì)
    電子制作(2019年19期)2019-11-23 08:41:36
    工程造價(jià)控制策略
    山東冶金(2019年3期)2019-07-10 00:54:04
    基于Multisim10和AD603的程控增益放大器仿真研究
    電子制作(2018年19期)2018-11-14 02:37:02
    現(xiàn)代企業(yè)會(huì)計(jì)的內(nèi)部控制策略探討
    中國電子信息產(chǎn)業(yè)轉(zhuǎn)移特征及驅(qū)動(dòng)因素——基于區(qū)域間投入產(chǎn)出表
    容錯(cuò)逆變器直接轉(zhuǎn)矩控制策略
    結(jié)合區(qū)域間差異性的水平集演化模型
    自增益電路在激光測距中的應(yīng)用
    精品国内亚洲2022精品成人| 少妇人妻一区二区三区视频| 久久久精品欧美日韩精品| 最近手机中文字幕大全| 久久久久久久国产电影| 大香蕉久久网| 日本色播在线视频| 国产精品女同一区二区软件| 久久精品久久久久久久性| 日韩欧美国产在线观看| 插逼视频在线观看| 亚洲国产精品专区欧美| 精品久久久久久久末码| 极品教师在线视频| 欧美日韩国产亚洲二区| 亚洲精品,欧美精品| 国产精品三级大全| 国产真实乱freesex| 亚洲天堂国产精品一区在线| 插逼视频在线观看| 亚洲18禁久久av| 一区二区三区免费毛片| 老师上课跳d突然被开到最大视频| 黄色配什么色好看| 久久精品国产自在天天线| 美女xxoo啪啪120秒动态图| 乱系列少妇在线播放| 最近最新中文字幕大全电影3| 国产精品乱码一区二三区的特点| 国产三级在线视频| 青青草视频在线视频观看| 最近2019中文字幕mv第一页| 国产精品蜜桃在线观看| 国产伦在线观看视频一区| 日韩精品青青久久久久久| 亚洲国产欧洲综合997久久,| 五月玫瑰六月丁香| 激情 狠狠 欧美| 精品欧美国产一区二区三| videos熟女内射| 夫妻性生交免费视频一级片| 亚洲伊人久久精品综合 | 亚洲国产精品成人久久小说| 久久久久久久久久久丰满| 黄色配什么色好看| 日日撸夜夜添| 亚洲婷婷狠狠爱综合网| 国产又黄又爽又无遮挡在线| 国产免费福利视频在线观看| 亚洲内射少妇av| 在现免费观看毛片| a级一级毛片免费在线观看| 寂寞人妻少妇视频99o| 我的女老师完整版在线观看| 日本欧美国产在线视频| 亚洲人成网站高清观看| 插逼视频在线观看| 长腿黑丝高跟| 99热6这里只有精品| 久久这里有精品视频免费| 国产精品一区二区性色av| 亚洲一级一片aⅴ在线观看| av视频在线观看入口| 蜜桃亚洲精品一区二区三区| h日本视频在线播放| 哪个播放器可以免费观看大片| 亚洲av电影在线观看一区二区三区 | 99久久精品热视频| 亚洲四区av| 久久久久久久亚洲中文字幕| 日本免费一区二区三区高清不卡| 国产精品三级大全| 一级av片app| 一个人看视频在线观看www免费| 欧美日本亚洲视频在线播放| 欧美激情在线99| 中文字幕亚洲精品专区| 两性午夜刺激爽爽歪歪视频在线观看| 精品国产三级普通话版| 99久久精品国产国产毛片| 老师上课跳d突然被开到最大视频| 熟女电影av网| 国产视频首页在线观看| 欧美成人精品欧美一级黄| 亚洲欧美日韩东京热| 美女大奶头视频| 最近中文字幕高清免费大全6| 久久草成人影院| 天天一区二区日本电影三级| 久久久久久久久大av| 性色avwww在线观看| 乱码一卡2卡4卡精品| 熟女电影av网| 日韩强制内射视频| 老司机福利观看| 久久久久国产网址| 美女黄网站色视频| 欧美xxxx性猛交bbbb| 日韩欧美在线乱码| 少妇丰满av| 国产免费视频播放在线视频 | 国语自产精品视频在线第100页| 久久精品久久久久久噜噜老黄 | 人妻夜夜爽99麻豆av| 长腿黑丝高跟| 蜜桃亚洲精品一区二区三区| 欧美三级亚洲精品| 国产国拍精品亚洲av在线观看| 亚洲国产色片| kizo精华| 婷婷六月久久综合丁香| 18禁裸乳无遮挡免费网站照片| 老师上课跳d突然被开到最大视频| 精品一区二区三区人妻视频| 变态另类丝袜制服| 青春草国产在线视频| 国产精品乱码一区二三区的特点| www.色视频.com| 午夜免费激情av| 午夜精品一区二区三区免费看| 人妻夜夜爽99麻豆av| 色综合站精品国产| 超碰97精品在线观看| 国内精品一区二区在线观看| 国产一区亚洲一区在线观看| 久久精品国产99精品国产亚洲性色| 亚洲最大成人中文| 国产成人免费观看mmmm| 少妇人妻精品综合一区二区| 少妇高潮的动态图| 国产精品一区二区三区四区久久| 桃色一区二区三区在线观看| 美女脱内裤让男人舔精品视频| 日韩三级伦理在线观看| 乱系列少妇在线播放| 亚洲无线观看免费| 亚洲国产精品专区欧美| 三级国产精品片| 国产成人freesex在线| 亚洲国产色片| 欧美人与善性xxx| 男女下面进入的视频免费午夜| 一区二区三区免费毛片| 日韩在线高清观看一区二区三区| 国产成人免费观看mmmm| 免费在线观看成人毛片| 精品少妇黑人巨大在线播放 | 大又大粗又爽又黄少妇毛片口| 日韩av在线免费看完整版不卡| 18禁在线播放成人免费| 亚洲国产精品久久男人天堂| 日韩av在线免费看完整版不卡| 七月丁香在线播放| 国产精品久久久久久精品电影小说 | 国产综合懂色| 亚洲第一区二区三区不卡| 国产黄a三级三级三级人| 青春草国产在线视频| 九草在线视频观看| 国产亚洲91精品色在线| 成人鲁丝片一二三区免费| 国产精品国产高清国产av| 久久久久久久久久黄片| 美女高潮的动态| 免费无遮挡裸体视频| 久久这里只有精品中国| 热99re8久久精品国产| 精品少妇黑人巨大在线播放 | a级一级毛片免费在线观看| 亚洲自拍偷在线| 一边摸一边抽搐一进一小说| 亚洲成av人片在线播放无| 黄片wwwwww| 春色校园在线视频观看| 欧美激情国产日韩精品一区| 神马国产精品三级电影在线观看| 五月玫瑰六月丁香| 欧美一区二区亚洲| 精品人妻视频免费看| 少妇裸体淫交视频免费看高清| or卡值多少钱| 麻豆成人午夜福利视频| 国产成人精品一,二区| 亚洲经典国产精华液单| 日本午夜av视频| 一级毛片久久久久久久久女| 日本黄色片子视频| 青春草视频在线免费观看| 精品酒店卫生间| 午夜日本视频在线| 亚洲伊人久久精品综合 | 国产精华一区二区三区| 一级av片app| 97超碰精品成人国产| 日本黄大片高清| www日本黄色视频网| 亚洲欧洲日产国产| 亚洲最大成人手机在线| 午夜精品国产一区二区电影 | 九色成人免费人妻av| 一区二区三区免费毛片| 国产91av在线免费观看| 国产精品一区二区性色av| 国产高清国产精品国产三级 | 九九热线精品视视频播放| 成人亚洲精品av一区二区| 日本-黄色视频高清免费观看| 中文字幕熟女人妻在线| 久久久久久久亚洲中文字幕| 搡女人真爽免费视频火全软件| 成人av在线播放网站| 国产精品乱码一区二三区的特点| 国产精品爽爽va在线观看网站| 99久久成人亚洲精品观看| 91午夜精品亚洲一区二区三区| 69av精品久久久久久| 欧美97在线视频| 国产午夜精品久久久久久一区二区三区| 亚洲精品国产av成人精品| 亚洲欧美成人精品一区二区| 少妇的逼好多水| 特大巨黑吊av在线直播| 亚洲av日韩在线播放| 91在线精品国自产拍蜜月| 大又大粗又爽又黄少妇毛片口| 久久国产乱子免费精品| 国产中年淑女户外野战色| 网址你懂的国产日韩在线| 亚洲国产欧美人成| 黄色欧美视频在线观看| 一级毛片aaaaaa免费看小| 22中文网久久字幕| 日本欧美国产在线视频| 黄色日韩在线| 偷拍熟女少妇极品色| 久久久精品大字幕| 国产精品av视频在线免费观看| 国产成人aa在线观看| 三级男女做爰猛烈吃奶摸视频| 成人漫画全彩无遮挡| 老女人水多毛片| 亚洲真实伦在线观看| 亚洲综合精品二区| 伦精品一区二区三区| 国产成人精品婷婷| 天堂影院成人在线观看| 只有这里有精品99| 久99久视频精品免费| h日本视频在线播放| 听说在线观看完整版免费高清| 国产又黄又爽又无遮挡在线| 国产伦精品一区二区三区视频9| 51国产日韩欧美| 老司机影院毛片| 综合色丁香网| 午夜免费男女啪啪视频观看| 亚洲av福利一区| 久久精品国产99精品国产亚洲性色| 午夜福利网站1000一区二区三区| 18禁裸乳无遮挡免费网站照片| 两性午夜刺激爽爽歪歪视频在线观看| 一级毛片电影观看 | 91aial.com中文字幕在线观看| 99久国产av精品国产电影| 亚洲最大成人中文| 综合色丁香网| 欧美日韩国产亚洲二区| 狂野欧美激情性xxxx在线观看| 乱码一卡2卡4卡精品| 黄片wwwwww| 久久久久性生活片| 亚洲av中文av极速乱| 一本久久精品| 久久鲁丝午夜福利片| 91久久精品国产一区二区三区| 久久精品国产鲁丝片午夜精品| 亚洲精品自拍成人| 亚洲欧美一区二区三区国产| 男的添女的下面高潮视频| 18禁裸乳无遮挡免费网站照片| 久久鲁丝午夜福利片| 中文字幕熟女人妻在线| 欧美97在线视频| 美女国产视频在线观看| 亚洲av电影不卡..在线观看| 中文欧美无线码| 51国产日韩欧美| 午夜福利在线观看免费完整高清在| 成人午夜高清在线视频| 久久久久久久久久成人| 国产亚洲最大av| 免费观看精品视频网站| 亚洲精品一区蜜桃| 久久精品国产自在天天线| 久久99热6这里只有精品| 国产极品精品免费视频能看的| 亚洲精品影视一区二区三区av| 午夜免费男女啪啪视频观看| 淫秽高清视频在线观看| 成人美女网站在线观看视频| 亚洲欧美成人综合另类久久久 | 亚洲av免费在线观看| 久久精品国产亚洲av天美| 日本午夜av视频| 色视频www国产| 国产精品综合久久久久久久免费| 国内精品宾馆在线| 小说图片视频综合网站| 日韩高清综合在线| 一级黄色大片毛片| 一个人观看的视频www高清免费观看| 亚洲国产日韩欧美精品在线观看| 亚州av有码| 尤物成人国产欧美一区二区三区| 亚洲精品,欧美精品| 在线免费十八禁| 高清午夜精品一区二区三区| 一边摸一边抽搐一进一小说| 国产黄片美女视频| 最近中文字幕高清免费大全6| 欧美日本视频| 国产高清国产精品国产三级 | 变态另类丝袜制服| 真实男女啪啪啪动态图| 日产精品乱码卡一卡2卡三| 在现免费观看毛片| 精品一区二区三区人妻视频| 99在线视频只有这里精品首页| 亚洲成人av在线免费| 91在线精品国自产拍蜜月| 国产午夜精品论理片| 97在线视频观看| 黄色欧美视频在线观看| 一边亲一边摸免费视频| 天堂av国产一区二区熟女人妻| 亚洲av福利一区| 国产一区有黄有色的免费视频 | 一级毛片电影观看 | av又黄又爽大尺度在线免费看 | 欧美bdsm另类| 女人被狂操c到高潮| 日韩av在线免费看完整版不卡| 成人av在线播放网站| 国产精品一区二区性色av| 毛片女人毛片| 我的老师免费观看完整版| 国产在视频线精品| 亚洲中文字幕日韩| 好男人在线观看高清免费视频| 国产一级毛片七仙女欲春2| 美女cb高潮喷水在线观看| 99视频精品全部免费 在线| 国模一区二区三区四区视频| 网址你懂的国产日韩在线| 91久久精品国产一区二区成人| 九草在线视频观看| 亚洲婷婷狠狠爱综合网| 高清午夜精品一区二区三区| 一本一本综合久久| 22中文网久久字幕| 国产乱来视频区| 亚洲精品自拍成人| 国产69精品久久久久777片| 2021天堂中文幕一二区在线观| 久久99热6这里只有精品| 免费看a级黄色片| 非洲黑人性xxxx精品又粗又长| 22中文网久久字幕| 最近2019中文字幕mv第一页| 久久99精品国语久久久| 成年免费大片在线观看| 日韩欧美精品v在线| 国产精品一二三区在线看| 日韩精品青青久久久久久| 精品无人区乱码1区二区| 久久精品人妻少妇| 最近的中文字幕免费完整| 亚洲精品亚洲一区二区| 少妇人妻精品综合一区二区| 黄片无遮挡物在线观看| 蜜桃亚洲精品一区二区三区| 搡女人真爽免费视频火全软件| 国产久久久一区二区三区| 免费播放大片免费观看视频在线观看 | 麻豆成人av视频| 欧美性猛交╳xxx乱大交人| 人人妻人人看人人澡| 男女国产视频网站| 91狼人影院| 日韩成人av中文字幕在线观看| 国产亚洲精品av在线| 激情 狠狠 欧美| 国产伦在线观看视频一区| 精品免费久久久久久久清纯| 亚洲成色77777| 美女内射精品一级片tv| 高清午夜精品一区二区三区| 成人av在线播放网站| 又粗又爽又猛毛片免费看| 成人性生交大片免费视频hd| 亚洲精品日韩在线中文字幕| 嘟嘟电影网在线观看| 国产不卡一卡二| 99久久人妻综合| 国产伦精品一区二区三区视频9| 天堂影院成人在线观看| 99在线人妻在线中文字幕| 99久久精品一区二区三区| 97在线视频观看| 日韩成人伦理影院| 日韩av在线大香蕉| 久久精品91蜜桃| 午夜福利视频1000在线观看| 国产女主播在线喷水免费视频网站 | 韩国高清视频一区二区三区| 日本色播在线视频| 国产在视频线在精品| 五月玫瑰六月丁香| 看黄色毛片网站| 我要搜黄色片| 亚洲色图av天堂| 插阴视频在线观看视频| 1000部很黄的大片| 99热全是精品| av国产免费在线观看| 久久精品国产亚洲av天美| 日韩大片免费观看网站 | 九九久久精品国产亚洲av麻豆| 直男gayav资源| 国产精品,欧美在线| 啦啦啦韩国在线观看视频| 国产国拍精品亚洲av在线观看| 久久久久久久久中文| 午夜亚洲福利在线播放| 在线a可以看的网站| 国产精品嫩草影院av在线观看| 中文字幕av成人在线电影| 日韩一区二区三区影片| 日韩精品有码人妻一区| 又爽又黄a免费视频| 亚洲av成人av| 国产午夜精品论理片| 一二三四中文在线观看免费高清| 欧美区成人在线视频| 亚洲最大成人手机在线| 午夜福利在线在线| 国产精品久久电影中文字幕| 97超视频在线观看视频| 成人午夜精彩视频在线观看| 秋霞伦理黄片| 丰满乱子伦码专区| 伦精品一区二区三区| 天天躁日日操中文字幕| 国产女主播在线喷水免费视频网站 | 18禁动态无遮挡网站| 国产伦精品一区二区三区视频9| 国产日韩欧美在线精品| 国产视频首页在线观看| 女人久久www免费人成看片 | av在线观看视频网站免费| 欧美成人免费av一区二区三区| 婷婷六月久久综合丁香| 一级毛片aaaaaa免费看小| 舔av片在线| 白带黄色成豆腐渣| 99久国产av精品国产电影| 久久久午夜欧美精品| 久热久热在线精品观看| 黑人高潮一二区| 99久久精品一区二区三区| 日韩成人av中文字幕在线观看| 一级黄色大片毛片| 国产麻豆成人av免费视频| 国产一区二区三区av在线| 精品久久久久久久久av| 日韩视频在线欧美| 听说在线观看完整版免费高清| 欧美日韩国产亚洲二区| 一级毛片久久久久久久久女| 亚洲无线观看免费| 午夜老司机福利剧场| 国产女主播在线喷水免费视频网站 | 亚洲国产精品久久男人天堂| 床上黄色一级片| 亚洲无线观看免费| 一区二区三区乱码不卡18| 中文天堂在线官网| 欧美成人精品欧美一级黄| 丰满人妻一区二区三区视频av| 国产午夜精品一二区理论片| 一个人观看的视频www高清免费观看| 国产成人精品久久久久久| 精品无人区乱码1区二区| 欧美日韩精品成人综合77777| 又黄又爽又刺激的免费视频.| 狂野欧美激情性xxxx在线观看| 亚洲色图av天堂| 亚洲国产色片| 少妇被粗大猛烈的视频| 大又大粗又爽又黄少妇毛片口| 舔av片在线| 夫妻性生交免费视频一级片| 三级毛片av免费| 免费播放大片免费观看视频在线观看 | 永久网站在线| 桃色一区二区三区在线观看| 高清日韩中文字幕在线| 少妇的逼水好多| 中文欧美无线码| 中文字幕免费在线视频6| 久久99热这里只有精品18| 国产又黄又爽又无遮挡在线| 成人漫画全彩无遮挡| 色网站视频免费| 在线免费十八禁| www日本黄色视频网| 午夜精品一区二区三区免费看| 亚洲av一区综合| 在线a可以看的网站| 热99re8久久精品国产| 亚洲美女视频黄频| 亚洲三级黄色毛片| 亚洲欧美日韩卡通动漫| 欧美一区二区精品小视频在线| 一级av片app| 国产高清视频在线观看网站| 久久精品国产亚洲av涩爱| 嫩草影院入口| 亚洲av中文av极速乱| 久久久久久久久久久丰满| 欧美性猛交黑人性爽| 国产精品1区2区在线观看.| 51国产日韩欧美| 国产伦精品一区二区三区视频9| a级一级毛片免费在线观看| 99久久无色码亚洲精品果冻| 男女啪啪激烈高潮av片| 内地一区二区视频在线| 国产精品一区二区三区四区久久| 天天一区二区日本电影三级| 国产精品久久久久久精品电影小说 | 亚洲欧美中文字幕日韩二区| 在线观看一区二区三区| 欧美一级a爱片免费观看看| 国产黄片视频在线免费观看| 边亲边吃奶的免费视频| 欧美高清性xxxxhd video| 亚州av有码| 不卡视频在线观看欧美| 男人和女人高潮做爰伦理| 欧美日韩精品成人综合77777| 一级毛片aaaaaa免费看小| 99久久无色码亚洲精品果冻| 久久欧美精品欧美久久欧美| 久久精品国产鲁丝片午夜精品| 嫩草影院精品99| 欧美日韩一区二区视频在线观看视频在线 | 久久精品夜夜夜夜夜久久蜜豆| 麻豆一二三区av精品| 91久久精品国产一区二区成人| 日韩人妻高清精品专区| 亚洲不卡免费看| 亚洲精品国产成人久久av| 一级黄色大片毛片| 成人鲁丝片一二三区免费| 97超视频在线观看视频| 少妇熟女aⅴ在线视频| 熟女人妻精品中文字幕| 久久综合国产亚洲精品| 综合色丁香网| 欧美成人一区二区免费高清观看| 少妇人妻精品综合一区二区| 波多野结衣巨乳人妻| 国产极品天堂在线| 淫秽高清视频在线观看| 又黄又爽又刺激的免费视频.| 国产一区二区在线观看日韩| 美女cb高潮喷水在线观看| 国产精品乱码一区二三区的特点| 噜噜噜噜噜久久久久久91| 欧美潮喷喷水| 婷婷色综合大香蕉| 国产久久久一区二区三区| ponron亚洲| 99久国产av精品国产电影| 最近中文字幕2019免费版| 青春草亚洲视频在线观看| 国产探花在线观看一区二区| 国产成人免费观看mmmm| 秋霞在线观看毛片| 欧美又色又爽又黄视频| 精品人妻熟女av久视频| 欧美日韩一区二区视频在线观看视频在线 | 色吧在线观看| 久久精品国产亚洲av天美| 国产一区二区在线观看日韩| 中文字幕制服av| 国产精品无大码| 国产精品一区www在线观看| 身体一侧抽搐| 国产亚洲精品久久久com| 插逼视频在线观看| 丰满少妇做爰视频| 精品不卡国产一区二区三区| av国产久精品久网站免费入址| 国产黄片美女视频| 亚洲高清免费不卡视频| 日本熟妇午夜| 久久精品久久精品一区二区三区| 尾随美女入室| 综合色丁香网| 伦理电影大哥的女人| 日韩av在线大香蕉| 成年版毛片免费区| 亚洲精品影视一区二区三区av|