• <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)用
    国产一区二区在线观看日韩| 亚洲成人久久性| 黄色日韩在线| 亚洲中文字幕一区二区三区有码在线看| 亚洲国产高清在线一区二区三| 悠悠久久av| 亚洲精品日韩av片在线观看| 久久精品国产亚洲av香蕉五月| 一区福利在线观看| 少妇高潮的动态图| 亚洲在线观看片| 国产伦精品一区二区三区四那| 亚洲午夜理论影院| 不卡一级毛片| 亚洲国产精品sss在线观看| 人人妻人人澡欧美一区二区| 人人妻人人澡欧美一区二区| 男女下面进入的视频免费午夜| 免费av毛片视频| 中文字幕人成人乱码亚洲影| 国产蜜桃级精品一区二区三区| 午夜免费成人在线视频| 人人妻人人澡欧美一区二区| 国产精品av视频在线免费观看| 日韩中字成人| 国产欧美日韩一区二区三| 欧美精品国产亚洲| 国产欧美日韩精品亚洲av| 在线免费观看的www视频| 床上黄色一级片| 最近中文字幕高清免费大全6 | 亚洲成人久久性| 精品一区二区三区人妻视频| 精品一区二区三区人妻视频| 国产精品一区二区性色av| 久久热精品热| 村上凉子中文字幕在线| 3wmmmm亚洲av在线观看| 亚洲成人中文字幕在线播放| 神马国产精品三级电影在线观看| 国产探花在线观看一区二区| 欧美激情久久久久久爽电影| 亚洲经典国产精华液单 | 老鸭窝网址在线观看| 欧美xxxx性猛交bbbb| 蜜桃久久精品国产亚洲av| 日日摸夜夜添夜夜添av毛片 | 91字幕亚洲| 国产高清激情床上av| 成年女人永久免费观看视频| а√天堂www在线а√下载| 18美女黄网站色大片免费观看| 好男人在线观看高清免费视频| 亚洲av电影不卡..在线观看| 欧美日本视频| 成人av一区二区三区在线看| 国产欧美日韩一区二区三| 九九热线精品视视频播放| 久久6这里有精品| 成人性生交大片免费视频hd| 色视频www国产| 国产成人啪精品午夜网站| 久久久久久久久中文| 757午夜福利合集在线观看| 欧美区成人在线视频| 精品久久久久久,| 久久久成人免费电影| 国产一区二区亚洲精品在线观看| 久久6这里有精品| 成年人黄色毛片网站| 免费在线观看亚洲国产| 黄色丝袜av网址大全| 在线播放无遮挡| 伦理电影大哥的女人| 91字幕亚洲| 国产精品久久久久久精品电影| 免费大片18禁| 久久久久国内视频| 午夜精品久久久久久毛片777| 欧洲精品卡2卡3卡4卡5卡区| 免费观看精品视频网站| 波多野结衣巨乳人妻| 免费黄网站久久成人精品 | 91在线观看av| 波野结衣二区三区在线| 好男人电影高清在线观看| 3wmmmm亚洲av在线观看| 久久精品国产亚洲av香蕉五月| 成年女人永久免费观看视频| 亚洲最大成人av| 亚洲av成人av| 亚洲七黄色美女视频| 日韩高清综合在线| 国产视频内射| 老女人水多毛片| netflix在线观看网站| 日韩欧美在线乱码| 亚洲电影在线观看av| 欧美又色又爽又黄视频| 草草在线视频免费看| 亚洲片人在线观看| 精品久久久久久久久av| 亚洲七黄色美女视频| 国产人妻一区二区三区在| 亚洲av免费高清在线观看| 成人性生交大片免费视频hd| 亚洲人成网站在线播| 久久人妻av系列| 亚洲在线观看片| 搞女人的毛片| 国产极品精品免费视频能看的| 此物有八面人人有两片| 一区福利在线观看| 亚洲精品粉嫩美女一区| 黄色一级大片看看| 成人特级黄色片久久久久久久| 成人三级黄色视频| 久久久久久大精品| 午夜福利在线观看吧| 国产成人欧美在线观看| 成人欧美大片| 久久久久久久久久成人| 欧美日韩综合久久久久久 | 99热这里只有是精品50| 国产亚洲精品av在线| 成年女人看的毛片在线观看| 精品久久久久久久久亚洲 | 亚洲18禁久久av| 综合色av麻豆| 最近最新免费中文字幕在线| 日本撒尿小便嘘嘘汇集6| 夜夜躁狠狠躁天天躁| 搞女人的毛片| 精品久久久久久久人妻蜜臀av| 欧美国产日韩亚洲一区| 成人鲁丝片一二三区免费| 桃色一区二区三区在线观看| 欧美成狂野欧美在线观看| 久久久国产成人免费| 成人亚洲精品av一区二区| 亚洲av电影不卡..在线观看| 久久精品国产亚洲av天美| 级片在线观看| 欧美+日韩+精品| 人人妻人人看人人澡| 国产探花在线观看一区二区| 波多野结衣巨乳人妻| 白带黄色成豆腐渣| ponron亚洲| 亚洲精品在线观看二区| 国产真实伦视频高清在线观看 | 少妇熟女aⅴ在线视频| 日日摸夜夜添夜夜添小说| 啦啦啦韩国在线观看视频| 亚洲av成人精品一区久久| 亚洲三级黄色毛片| 丁香六月欧美| 欧美不卡视频在线免费观看| 国产一级毛片七仙女欲春2| 97碰自拍视频| 日韩中字成人| 日本精品一区二区三区蜜桃| 欧美一区二区亚洲| 亚洲国产精品成人综合色| 制服丝袜大香蕉在线| 99国产精品一区二区蜜桃av| 给我免费播放毛片高清在线观看| 国产日本99.免费观看| 又黄又爽又刺激的免费视频.| 成人国产综合亚洲| 神马国产精品三级电影在线观看| 欧美午夜高清在线| 在线观看美女被高潮喷水网站 | 直男gayav资源| 一个人看视频在线观看www免费| 99久久久亚洲精品蜜臀av| 国产高清激情床上av| 一级毛片久久久久久久久女| 久9热在线精品视频| www.熟女人妻精品国产| 韩国av一区二区三区四区| 999久久久精品免费观看国产| 国产不卡一卡二| 欧美性猛交╳xxx乱大交人| 麻豆成人午夜福利视频| 淫妇啪啪啪对白视频| h日本视频在线播放| 99视频精品全部免费 在线| 日韩人妻高清精品专区| 精品久久久久久,| 亚洲av免费在线观看| 亚洲国产欧洲综合997久久,| a级毛片a级免费在线| 久久久色成人| 婷婷丁香在线五月| 国产精品一区二区免费欧美| 欧美激情久久久久久爽电影| 又紧又爽又黄一区二区| av在线观看视频网站免费| 亚洲片人在线观看| 99在线人妻在线中文字幕| 国产精品久久视频播放| 欧美色欧美亚洲另类二区| 国产精品一区二区免费欧美| 成人精品一区二区免费| 国内精品一区二区在线观看| 91久久精品电影网| 少妇被粗大猛烈的视频| 免费高清视频大片| 色精品久久人妻99蜜桃| 亚洲精品久久国产高清桃花| 麻豆国产av国片精品| 中文字幕精品亚洲无线码一区| 成人av一区二区三区在线看| 91麻豆精品激情在线观看国产| 九九热线精品视视频播放| 免费人成视频x8x8入口观看| 欧美日韩福利视频一区二区| 亚洲精品乱码久久久v下载方式| 在线十欧美十亚洲十日本专区| 国产精品1区2区在线观看.| 久久久久国产精品人妻aⅴ院| 欧美日韩综合久久久久久 | 中文字幕人妻熟人妻熟丝袜美| 精品福利观看| 丁香欧美五月| 亚洲无线在线观看| 久久性视频一级片| 两人在一起打扑克的视频| 国产精品电影一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 人人妻,人人澡人人爽秒播| 床上黄色一级片| 嫩草影院入口| 99久久成人亚洲精品观看| av天堂中文字幕网| 久久精品久久久久久噜噜老黄 | 国内精品一区二区在线观看| 99久久久亚洲精品蜜臀av| 亚洲精品色激情综合| 桃红色精品国产亚洲av| 桃色一区二区三区在线观看| 国产av麻豆久久久久久久| 欧美日韩黄片免| 香蕉av资源在线| 十八禁人妻一区二区| 国产成人啪精品午夜网站| 757午夜福利合集在线观看| 在线观看免费视频日本深夜| 亚洲成人精品中文字幕电影| 91麻豆精品激情在线观看国产| 成人特级av手机在线观看| 十八禁人妻一区二区| 一级av片app| av黄色大香蕉| 少妇熟女aⅴ在线视频| 久久欧美精品欧美久久欧美| 日本精品一区二区三区蜜桃| 嫩草影院入口| 国产精品亚洲av一区麻豆| 国产在线男女| 亚洲精品亚洲一区二区| 哪里可以看免费的av片| 最近中文字幕高清免费大全6 | 亚洲国产精品sss在线观看| 精品久久国产蜜桃| 亚洲欧美日韩高清专用| 真实男女啪啪啪动态图| 亚洲人成网站在线播放欧美日韩| 午夜a级毛片| 天美传媒精品一区二区| 国产成人av教育| 欧美一级a爱片免费观看看| 国产精品一区二区三区四区免费观看 | 波多野结衣巨乳人妻| 国产精品一区二区三区四区免费观看 | 2021天堂中文幕一二区在线观| 亚洲激情在线av| 最近最新免费中文字幕在线| 欧美中文日本在线观看视频| 中文资源天堂在线| 日韩人妻高清精品专区| 五月玫瑰六月丁香| 日韩欧美精品免费久久 | 美女大奶头视频| 国产高清视频在线观看网站| 免费搜索国产男女视频| 国产精品亚洲av一区麻豆| 窝窝影院91人妻| 1000部很黄的大片| 琪琪午夜伦伦电影理论片6080| 午夜老司机福利剧场| 精品午夜福利视频在线观看一区| 在线观看舔阴道视频| 91久久精品电影网| 综合色av麻豆| 两人在一起打扑克的视频| 人人妻人人看人人澡| 悠悠久久av| 天堂动漫精品| 亚洲人成电影免费在线| av专区在线播放| 国产高清视频在线播放一区| 亚洲国产欧洲综合997久久,| 亚洲成av人片免费观看| 九色国产91popny在线| 亚洲国产高清在线一区二区三| 在线观看一区二区三区| 天天躁日日操中文字幕| 婷婷丁香在线五月| 亚洲中文日韩欧美视频| 国内精品美女久久久久久| 日韩欧美一区二区三区在线观看| 噜噜噜噜噜久久久久久91| 午夜精品在线福利| 亚洲在线自拍视频| 三级男女做爰猛烈吃奶摸视频| 久9热在线精品视频| 欧美一级a爱片免费观看看| 熟女电影av网| 深夜精品福利| 午夜福利视频1000在线观看| 99久国产av精品| 久久人妻av系列| 午夜福利在线观看免费完整高清在 | 国产高清激情床上av| 国产精品久久久久久精品电影| 成人无遮挡网站| 一卡2卡三卡四卡精品乱码亚洲| 黄色配什么色好看| 丰满的人妻完整版| 精品免费久久久久久久清纯| 夜夜夜夜夜久久久久| 国语自产精品视频在线第100页| 国产精品久久久久久精品电影| 午夜福利视频1000在线观看| 怎么达到女性高潮| 久久精品人妻少妇| 日本a在线网址| xxxwww97欧美| 国产精品一区二区免费欧美| 国产精品久久久久久久电影| 久久精品夜夜夜夜夜久久蜜豆| 99riav亚洲国产免费| av在线老鸭窝| 国产精品久久久久久久电影| 国产高清视频在线观看网站| 国产淫片久久久久久久久 | 三级国产精品欧美在线观看| 色综合站精品国产| 欧美成人免费av一区二区三区| 亚洲成人免费电影在线观看| 亚洲国产欧美人成| 性欧美人与动物交配| 一个人看的www免费观看视频| 久久久久久大精品| 99热这里只有是精品50| 亚洲av成人不卡在线观看播放网| 我要搜黄色片| 麻豆av噜噜一区二区三区| 欧美日本视频| 深夜a级毛片| 中亚洲国语对白在线视频| 赤兔流量卡办理| 日本与韩国留学比较| 午夜福利在线观看吧| 国产探花在线观看一区二区| 日韩精品青青久久久久久| 亚洲欧美日韩卡通动漫| 最近最新中文字幕大全电影3| 好男人电影高清在线观看| 午夜亚洲福利在线播放| 国产aⅴ精品一区二区三区波| 色哟哟哟哟哟哟| 免费观看人在逋| 中文在线观看免费www的网站| 嫩草影视91久久| 99久久无色码亚洲精品果冻| 免费人成在线观看视频色| 亚洲成人中文字幕在线播放| 99在线视频只有这里精品首页| 亚洲va日本ⅴa欧美va伊人久久| 精品乱码久久久久久99久播| 国产中年淑女户外野战色| 亚洲经典国产精华液单 | 麻豆成人午夜福利视频| 国产免费男女视频| 五月玫瑰六月丁香| 亚洲精品456在线播放app | 亚洲欧美日韩东京热| 亚洲av美国av| 成人鲁丝片一二三区免费| 中文亚洲av片在线观看爽| 深爱激情五月婷婷| 中文字幕免费在线视频6| 男人舔奶头视频| 嫁个100分男人电影在线观看| 内地一区二区视频在线| 老女人水多毛片| 国产aⅴ精品一区二区三区波| 日本五十路高清| 国产精品久久视频播放| 亚洲欧美日韩东京热| 我的女老师完整版在线观看| 国产成人欧美在线观看| 精品人妻偷拍中文字幕| 99riav亚洲国产免费| 天天一区二区日本电影三级| 国产精品爽爽va在线观看网站| 精品一区二区三区av网在线观看| 九九久久精品国产亚洲av麻豆| 非洲黑人性xxxx精品又粗又长| 国产精品一区二区三区四区免费观看 | 69人妻影院| 国产欧美日韩精品一区二区| 亚洲激情在线av| 成人鲁丝片一二三区免费| 日本免费一区二区三区高清不卡| 欧美在线一区亚洲| 免费一级毛片在线播放高清视频| eeuss影院久久| 一进一出抽搐动态| 黄色日韩在线| 制服丝袜大香蕉在线| av在线老鸭窝| 欧美又色又爽又黄视频| 国产乱人视频| 欧美成人性av电影在线观看| 午夜久久久久精精品| av天堂在线播放| 国产成人影院久久av| 少妇高潮的动态图| 国产黄a三级三级三级人| 婷婷丁香在线五月| 中文字幕av成人在线电影| 久久午夜福利片| 麻豆一二三区av精品| 最新中文字幕久久久久| 51国产日韩欧美| 婷婷精品国产亚洲av| 亚洲av第一区精品v没综合| 可以在线观看毛片的网站| 五月伊人婷婷丁香| 别揉我奶头~嗯~啊~动态视频| 国产精品亚洲美女久久久| 一a级毛片在线观看| 嫁个100分男人电影在线观看| 99精品在免费线老司机午夜| 免费看a级黄色片| a在线观看视频网站| 国产熟女xx| 国产精品,欧美在线| 日本黄色片子视频| 欧美中文日本在线观看视频| 在线看三级毛片| 黄色日韩在线| 偷拍熟女少妇极品色| 精品人妻一区二区三区麻豆 | 精品一区二区免费观看| 欧美中文日本在线观看视频| 国产一区二区三区视频了| 深爱激情五月婷婷| 十八禁国产超污无遮挡网站| 午夜福利在线在线| 俺也久久电影网| 久久精品国产亚洲av涩爱 | 麻豆国产97在线/欧美| 蜜桃久久精品国产亚洲av| 国产精品久久久久久亚洲av鲁大| 久久午夜福利片| 欧美乱色亚洲激情| 久久精品综合一区二区三区| 天堂影院成人在线观看| 亚洲成人精品中文字幕电影| 国产美女午夜福利| 一个人免费在线观看的高清视频| 少妇高潮的动态图| 精品午夜福利在线看| 亚洲av不卡在线观看| 每晚都被弄得嗷嗷叫到高潮| 成人午夜高清在线视频| 日日夜夜操网爽| 日韩欧美一区二区三区在线观看| xxxwww97欧美| 12—13女人毛片做爰片一| 欧美成人a在线观看| 好看av亚洲va欧美ⅴa在| 天堂√8在线中文| 校园春色视频在线观看| av天堂在线播放| 亚洲国产色片| 亚洲中文字幕日韩| 日日干狠狠操夜夜爽| 国产精品不卡视频一区二区 | 99riav亚洲国产免费| 欧美激情久久久久久爽电影| 国产一区二区激情短视频| 国产高清三级在线| 国产蜜桃级精品一区二区三区| aaaaa片日本免费| 成年女人看的毛片在线观看| 69人妻影院| 91九色精品人成在线观看| 国产国拍精品亚洲av在线观看| 中文字幕免费在线视频6| 成人三级黄色视频| 99热这里只有是精品在线观看 | 一个人观看的视频www高清免费观看| a级一级毛片免费在线观看| 国产欧美日韩一区二区三| 国产精品久久久久久久久免 | 欧美精品国产亚洲| 欧美高清性xxxxhd video| av中文乱码字幕在线| 亚洲人与动物交配视频| 欧美日韩黄片免| 啪啪无遮挡十八禁网站| 日韩国内少妇激情av| a级毛片a级免费在线| 欧美最黄视频在线播放免费| 国产私拍福利视频在线观看| 色尼玛亚洲综合影院| 亚洲美女视频黄频| 色视频www国产| 午夜福利高清视频| 国产真实乱freesex| 亚洲 欧美 日韩 在线 免费| 啪啪无遮挡十八禁网站| 久久久久久久久久成人| 一个人免费在线观看的高清视频| 国产探花在线观看一区二区| 91麻豆精品激情在线观看国产| 亚洲第一电影网av| 少妇被粗大猛烈的视频| 亚洲成a人片在线一区二区| 身体一侧抽搐| 色精品久久人妻99蜜桃| 少妇高潮的动态图| 91久久精品国产一区二区成人| 亚洲,欧美,日韩| 听说在线观看完整版免费高清| 看片在线看免费视频| 亚洲熟妇熟女久久| 成人欧美大片| 久久精品综合一区二区三区| 精品国内亚洲2022精品成人| 亚洲国产欧洲综合997久久,| 真实男女啪啪啪动态图| 亚洲av.av天堂| 国产大屁股一区二区在线视频| 搡女人真爽免费视频火全软件 | 熟妇人妻久久中文字幕3abv| 亚洲七黄色美女视频| 国产黄片美女视频| 性色avwww在线观看| 身体一侧抽搐| 亚洲人与动物交配视频| 国产男靠女视频免费网站| 国产精品久久久久久久久免 | 亚洲成人免费电影在线观看| 国产亚洲精品久久久com| 一区二区三区高清视频在线| 中文字幕精品亚洲无线码一区| av中文乱码字幕在线| 午夜福利成人在线免费观看| 国产免费av片在线观看野外av| 国产 一区 欧美 日韩| 国产av不卡久久| 在线国产一区二区在线| 18禁裸乳无遮挡免费网站照片| 国产男靠女视频免费网站| 国产在线精品亚洲第一网站| 综合色av麻豆| 一区二区三区激情视频| 亚洲精品456在线播放app | 成人av在线播放网站| 一a级毛片在线观看| 亚洲中文日韩欧美视频| 嫁个100分男人电影在线观看| 久久国产乱子伦精品免费另类| 怎么达到女性高潮| 搡老妇女老女人老熟妇| 午夜两性在线视频| 婷婷丁香在线五月| 国产av不卡久久| 少妇高潮的动态图| 国产人妻一区二区三区在| 日本 av在线| 色综合婷婷激情| 青草久久国产| 午夜免费成人在线视频| 国产成人啪精品午夜网站| 天堂动漫精品| 亚洲午夜理论影院| 亚洲欧美日韩卡通动漫| 免费av不卡在线播放| a在线观看视频网站| 午夜福利在线观看免费完整高清在 | 黄片小视频在线播放| 亚洲av成人精品一区久久| 亚洲国产欧洲综合997久久,| 窝窝影院91人妻| 久久久精品欧美日韩精品| 91狼人影院| 亚洲av熟女| 在线天堂最新版资源| 简卡轻食公司| 亚洲av成人不卡在线观看播放网| 在线看三级毛片| 成人国产一区最新在线观看| 99热6这里只有精品| 精品人妻视频免费看| 亚洲男人的天堂狠狠| 少妇被粗大猛烈的视频| 日韩欧美国产在线观看|