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

海上風電傳輸中MMC的模型預測控制

2015-10-11 08:59:16周云夏向陽孔祥霽范家銘趙相濤
中南大學學報(自然科學版) 2015年10期
關鍵詞:交流模型

周云,夏向陽,孔祥霽,范家銘,趙相濤

?

海上風電傳輸中MMC的模型預測控制

周云1,夏向陽1,孔祥霽1,范家銘1,趙相濤2

(1. 長沙理工大學電氣與信息工程學院,湖南長沙,410077;2. 國網(wǎng)山東省電力公司濟寧供電公司,山東濟寧,272000)

針對海上風電柔性直流輸電系統(tǒng)(HVDC)中模塊化多電平換流器(MMC)的循環(huán)電流抑制及MMC中子模塊電容電壓平衡問題,提出一種模型預測控制(MPC)策略。該方法根據(jù)系統(tǒng)的離散時間數(shù)學模型,開發(fā)對應的預測控制方式,將每個MMC單元中最佳的開關狀態(tài)問題轉(zhuǎn)化為求目標函數(shù)的最優(yōu)化問題,以抑制環(huán)流,并實現(xiàn)MMC單元的電容電壓平衡。最后利用MATLAB/SIMULINK進行仿真驗證,研究結(jié)果表明:該模型預測控制策略是有效的和可行的,且實現(xiàn)原理簡單,適用于海上風電傳輸?shù)腗MC-HVDC系統(tǒng)。

模塊化多電平換流器;模型預測控制;電容電壓平衡;循環(huán)電流控制;海上風電

在可再生能源開發(fā)中,風力發(fā)電技術(shù)是最成熟、最具商業(yè)化發(fā)展前景的發(fā)電形式之一[1],近年來,由于土地資源短缺,海上風電憑借其風速高、噪音小、風力穩(wěn)定及發(fā)電量大等諸多優(yōu)勢,在世界范圍內(nèi)得到迅速發(fā)展[2?3]。截至2011年底,全球海上風電總裝機容量達到4.954 GW,其中,我國2011年海上風電累計為241.3 MW。另外,我國提出了到2020年海上風電總裝機容量為30 GW的目標,發(fā)展海上風電已成為未來風電發(fā)展的主要方向。隨著海上風電機組容量和風電場的規(guī)模不斷增大,基于模塊化多電平換流器的柔性直流輸電(modular multilevel converter-high voltage direct current, MMC-HVDC)系統(tǒng)得到越來越多的關注和研究[4?6],與VSC-HVDC相比,MMC-HVDC系統(tǒng)有效地減少了輸出電壓的諧波含量,降低了功率器件開關頻率和損耗,實現(xiàn)高電壓大功率的運行[7?8],非常適合于長距離、大容量的海上風電傳輸。然而,MMC拓撲自身也存在一些不足:由于分布式布置的儲能電容,其功率單元電容電壓脈動較大以及各相之間能量分配不平衡,這不僅會使本來正弦的橋臂電流發(fā)生畸變,而且還提高了器件的額定電流容量,進一步增加了系統(tǒng)成本等[9]。MMC子模塊間電容電壓的均衡分配與橋臂間存在的循環(huán)電流問題成為MMC亟需解決的難點,有必要采用合適的控制策略對電容電壓平衡和環(huán)流進行協(xié)同控制。丁冠軍等[10]介紹了MMC子模塊的拓撲機制與調(diào)制策略,采用基頻開關調(diào)制時子模塊導通角度的計算方法,但沒有給出與其配合的電容均壓策略。Hagiwara等[11]提出對調(diào)制波進行微調(diào)的形式來實現(xiàn)電容電壓的平衡,但此方法需要與脈沖寬度調(diào)制相配合,這使得換流器每個子模塊所用的調(diào)制波不一樣,子模塊越多就越難以實現(xiàn)。屠卿瑞等[9]推導了MMC在運行中相間環(huán)流的解析表達式,并針對環(huán)流抑制給出了MMC橋臂電抗器的計算方法。梁帥奇等[12?14]給出了一種基于上橋臂和下橋臂的平均直流電容的電壓差值進行環(huán)流抑制的方法,并通過實驗驗證,但該方法需用到所有功率單元的直流電容電壓,對系統(tǒng)信號采集環(huán)節(jié)提出了較高要求。為此,本文針對應用于海上風電柔性直流輸電領域的MMC子模塊電容電壓平衡及換流器內(nèi)部存在的循環(huán)電流抑制問題,提出一種模型預測控制(MPC)策略,根據(jù)柔性直流輸電系統(tǒng)的離散時間數(shù)學模型,開發(fā)對應于離散時間模型的預測模型。利用它的目標函數(shù)最優(yōu)化技術(shù),并使用預測模型來選擇每個MMC單元中最佳的開關狀態(tài),來抑制循環(huán)電流,并通過冗余開關狀態(tài)達到電容電壓平衡。最后,通過MATLAB/SIMULINK仿真驗證該策略的有效性和可行性。

1 MMC-HVDC的基本結(jié)構(gòu)

圖1所示為MMC-HVDC系統(tǒng)的結(jié)構(gòu)圖。該系統(tǒng)由2個換流站及直流線路組成,送端、受端換流器均采用MMC拓撲,則2個換流器具有相同的結(jié)構(gòu),分別為圖1中MMC1和MMC2。MMC的交流側(cè)通過串聯(lián)電阻器和電感器以及三相變壓器連接到公用電網(wǎng),交流系統(tǒng)由等值電源s和等值電感AC表示,直流側(cè)電壓用dc表示。每個MMC單元的同步控制信號從對應的變壓器低壓側(cè)提取。

圖1 MMC-HVDC系統(tǒng)結(jié)構(gòu)

圖2所示為MMC主電路結(jié)構(gòu)圖。由圖2可知:每個MMC單元各相均有2個橋臂(即上橋臂和下橋臂),其中每個橋臂由個相互連接且結(jié)構(gòu)相同的半橋式子模塊SM(sub-module)與一系列電抗器串聯(lián)構(gòu)成,相臂內(nèi)的電抗器是用來抑制橋臂間因總直流電壓差異引起的相間環(huán)流。系統(tǒng)結(jié)構(gòu)具有高度模塊化特點,子模塊之間可以互換,具有很強的擴展性和靈活性。

圖2 MMC主電路結(jié)構(gòu)圖

2 MMC的模型預測控制策略

2.1 MMC基本運行原理

構(gòu)成模塊化多電平的基本功率單元SM由1個直流儲能電容和2個絕緣柵雙極晶體管(insulated gate bipolar transistor,IGBT)半橋組成,其拓撲如圖3所示。圖3中:SM為子模塊的輸出電壓;SM為子模塊的充放電電流;c為子模塊的瞬時電容電壓。

圖3 單個子模塊結(jié)構(gòu)

從圖3可知:每個功率單元都是兩端元件,通過2個開關器件T1和T2的作用;SM可以實現(xiàn)在電容電壓c與0之間的切換,從而調(diào)整橋臂輸出電平的個數(shù)。子模塊4種常用的工作狀態(tài)如表1所示。

表1 MMC子模塊的常用工作狀態(tài)

理論上,各子模塊構(gòu)成的橋臂電壓均可等效成1個可控電壓源Pjk和Njk(為a,b和c,表示三相;為1和2,分別對應MMC1和MMC2系統(tǒng))[15];下標P表示上橋臂;下標N表示下橋臂;Pjk和Njk為相應橋臂電流。

為了維持直流電壓穩(wěn)定,一般要求MMC同相上、下2個橋臂的子模塊互補對稱投入,即滿足

在任意時刻,每一相有且僅有個子模塊導通個子模塊關斷,則開關管的導通個數(shù)必須滿足

假設每個SM電容電壓理想調(diào)節(jié)為dc/,在MMC中的各相端子,1個理想的(+1)的電平波形,相對于1個虛構(gòu)的直流側(cè)中點,交流側(cè)的輸出電壓滿足如下關系[9]:

式中:v為MMC中交流側(cè)的相電壓。由式(3)可知:v在dc/2和?dc/2范圍內(nèi),以步長為dc/呈階梯狀變化。在實際運行過程中,MMC中各個SM電容電壓應被實時監(jiān)測且維持平衡,即通過使用有源電壓平衡技術(shù)使各個子模塊電容電壓維持其額定值dc/。

根據(jù)MMC的三相對稱性,以MMC1為例,環(huán)流等值電路圖如圖4所示。其中:i為相電流。相應的臂電流描述如下:

其中:Cir jk為流過MMC中相的循環(huán)電流。環(huán)流分量疊加在橋臂電流中,增加了器件的額定容量,增大了換流器的損耗,導致子模塊電容電壓的波動,嚴重時還會造成系統(tǒng)不穩(wěn)定。因此,為了最大限度地減少其負面影響,應該設計相應的控制策略對環(huán)流進行有效抑制。

圖4 MMC1環(huán)流等值電路圖

2.2 MMC的數(shù)學模型

據(jù)基爾霍夫電壓定理,MMC中相動態(tài)數(shù)學方程如下:

式中:v為變壓器低壓側(cè)的電網(wǎng)電壓。根據(jù)式(4)~(6),推導交流側(cè)相電流i與環(huán)流Cir jk分別為

由式(7)和(8)可知:交流側(cè)相電流受上、下橋臂電壓之差的影響,而上、下橋臂電壓之和與直流電壓不相等是內(nèi)部環(huán)流產(chǎn)生的根本原因。

圖2中MMC單元的SM電容電壓的動態(tài)開關描述為

其中:=1,2,…,;=,,;=1,2。若在上橋臂中的子模塊,則i=Pjk;若在下橋臂中的子模塊,i=Njk。

2.3 模型預測控制策略

模型預測控制(MPC)是一類特殊的控制,具有控制效果好、魯棒性強等優(yōu)點,其本質(zhì)是在每一個采樣瞬間通過求解1個有限時域開環(huán)最優(yōu)控制問題。在前面MMC的數(shù)學模型基礎上,本文提出一種改進的MPC控制策略,控制交流側(cè)電流,同時調(diào)節(jié)SM的電容電壓平衡,并抑制循環(huán)電流。其實現(xiàn)過程如下:首先,為MMC變量正向預測其離散時間模型;然后,定義與控制目標相關的目標函數(shù);最后,選擇換流器所有可能的最好狀態(tài)來評估定義的目標函數(shù),得到所定義的目標函數(shù)的最小值[16]。

2.3.1 交流側(cè)電流控制

電流控制的目的是調(diào)節(jié)MMC交流側(cè)電流與參考值一致(即參考電流跟蹤)。在式(7)的基礎上,假設1個采樣周期s,用歐拉近似推導出的電流導數(shù)作為MMC交流側(cè)電流的預測模型如下:

式中:′=+/2,′=+′/s;i(+s)為交流側(cè)電流的預測值;i()為實際測定值;v(+s)為變壓器低壓側(cè)電網(wǎng)電壓的預測值,假設s足夠小時,可近似為測定值v();Pjk(+s)和Njk(+s)分別為上、下橋臂電壓的預測值,它的計算可通過增加一步向前預測上橋臂和下橋臂投入子模塊的電容電壓來獲得。

為了減小預測電流和參考電流之間的誤差,定義與電流誤差有關的目標函數(shù)為

式中:iref為參考電流,它根據(jù)預先設定傳遞到交流系統(tǒng)的有功和無功功率求得。在理想情況下,若交流側(cè)電流跟蹤參考值,目標函數(shù)J達到其最小值0,則這種情況被作為交流側(cè)電流控制的理想狀態(tài)。在每個采樣期間,J為MMC所有可能的開關狀態(tài)計算并比較,從而得到式(11)的最小值,作為下一個開關周期最佳的開關狀態(tài)。

2.3.2 電容電壓平衡

(13)

2.3.3 循環(huán)電流控制

基于式(8),推導出循環(huán)電流的離散時間預測模型如下:

改進的MPC策略目的是通過添加式(15)中第3項來抑制循環(huán)循環(huán)電流,與循環(huán)電流相關的目標函數(shù)如下:

圖5 MMC單元的MPC策略框圖

3 仿真分析

3.1 仿真模型及參數(shù)

為了驗證海上風電柔性直流輸電系統(tǒng)中模型預測控制策略的有效性和可行性,利用MATLAB/SIMULINK軟件對三相MMC模型進行仿真驗證。采用MMC單元精確開關模型,主電路結(jié)構(gòu)如圖1所示。選定每相橋臂級聯(lián)的功率單元數(shù)為4,系統(tǒng)參數(shù)如下。

1) MMC系統(tǒng)參數(shù):N=50 MV?A,額定頻率1=2=50 Hz,子模塊電容=3 300 μF,直流側(cè)電壓dc/2=30 kV,橋臂串聯(lián)電抗3 mH。

2) 交流系統(tǒng)參數(shù):交流系統(tǒng)電壓幅值N= 138 kV,線路阻抗參數(shù)=0.03 ΩmH,采用換流變壓器Yn/Δ接法,變壓器額定功率為55 MV?A,漏抗標么值T=0.05,兩側(cè)交流系統(tǒng)短路比為5;采樣周期s=100 μs

3.2 仿真結(jié)果

最初,海上風電傳輸MMC-HVDC系統(tǒng)處于穩(wěn)定狀態(tài),40 MW的電能從交流系統(tǒng)1流向交流系統(tǒng)2。Vref=60 kV,所有的MMC單元在單位功率因數(shù)條件下運行,模型預測控制的環(huán)流抑制功能在初始狀態(tài)是禁用的(即式(16)中的系數(shù))。在為0.1 s和 0.4 s時,開啟MMC單元的環(huán)流抑制功能(系數(shù)),仿真結(jié)果如圖6和圖7所示。

(a) 環(huán)流iCir;(b) 橋臂電流i;(c) 電容電壓Vcaps; (d) 交流輸出相電流iabc;(e) 交流輸出線電壓Vtabe

(a) 環(huán)流iCir;(b) 橋臂電流i;(c) 電容電壓Vcaps; (d) 交流輸出相電流iabc;(e) 交流輸出線電壓Vtabe

MMC1中的環(huán)流和A相上、下橋臂的電流圖分別如圖6(a)和6(b)所示。從圖6(a)和6(b)可見:在=0.1 s開啟模型預測控制的環(huán)流抑制功能之后,環(huán)流和上下橋臂的電流顯著減小,抑制效果明顯。

圖6(c)顯示了MMC1設置在參考值(dc/4時其子模塊的電容電壓,驗證了所提出的模型預測控制策略對維持子模塊電容電壓在其額定值的有效性。從電容電壓波形可以看出,子模塊的電容電壓紋波在=0.1 s后顯著減少,這是由于抑制了環(huán)流。

MMC1的交流輸出相電流波形和交流輸出線電壓波形圖分別如圖6(d)和6(e)所示。從圖6(d)和6(e)可見:當=0.1 s時,由于環(huán)流的影響,子模塊的電容電壓紋波增加,導致交流側(cè)電壓和電流波形失真;在=0.4 s開啟模型預測控制功能后,改善了交流側(cè)電壓和電流的波形質(zhì)量,得到了理想的波形。

在=0.4 s開啟環(huán)流抑制功能時,MMC2中相應的波形如圖7所示。圖7(a)~(e)所示圖形分別與圖6(a)~(e)所示波形有密切聯(lián)系。

仿真圖驗證了所提出的模型預測控制策略抑制MMC單元環(huán)流的有效性,并實現(xiàn)MMC單元的電容電壓均壓,同時也反映了環(huán)流和電容電壓不平衡對MMC單元的運行有不利影響,突出了環(huán)流抑制和電容電壓平衡的重要性。

4 結(jié)論

1) 提出一種模型預測控制(MPC)策略用于海上風電的MMC-HVDC系統(tǒng)。在模型預測控制理論的基礎上,根據(jù)MMC的離散時間數(shù)學模型,開發(fā)對應于離散時間模型的預測模型,并定義相對應的目標函數(shù),通過求解這個目標函數(shù)的最優(yōu)化問題,得到每個MMC單元中開關器件的最佳開關狀態(tài)來抑制循環(huán)電流,并實現(xiàn)MMC單元的電容電壓平衡。

2) 利用MATLAB/SIMULINK軟件對五電平MMC-HVDC系統(tǒng)的進行仿真,結(jié)果證明了該模型預測控制策略的有效性和可行性。該策略可大大降低輸出波形的畸變程度,使其更逼近正弦波;實現(xiàn)原理簡單,非常適用于海上風電傳輸?shù)腗MC-HVDC系統(tǒng)。

[1] Noro Y, Arai J, Takagi K, et al. System study of direct current power transmission system connected to wind farm[C]// Transmission and Distribution Conference and Exhibition. Yokohama, Japan: IEEE, 2002: 944?949.

[2] 宋礎, 劉漢中. 海上風力發(fā)電廠開發(fā)現(xiàn)狀及發(fā)展趨勢[J]. 上海電力, 2008, 2(1): 59?63. SONG Chu, LIU Hanzhong. Development situation and trend of offshore wind farm[J]. Shanghai Electric Power, 2008, 2(1): 59?63.

[3] 王志新. 海上風力發(fā)電技術(shù)[M]. 北京: 機械工業(yè)出版社, 2012: 1?16. WANG Zhixin. Offshore wind energy generation technology[M]. Beijing: China Machine Press, 2012: 1?16.

[4] Marquard T R, Lesnicar A. A new modular voltage source inverter topology[C/CD]//European Conference on Power Electronics and Applications. Toulouse, France: EPE, 2003: 72?75.

[5] 張振華, 江道灼. 基于模塊化多電平變流器的STATCOM研究[J]. 電力自動化設備, 2012, 32(2): 62?66. ZHANG Zhenhua, JIANG Daozhuo. STATCOM based on modularized multilevel converters[J]. Electric Power Automation Equipment, 2012, 32(2): 62?66.

[6] 何大清, 蔡旭. 模塊化多電平變流器的限幅控制和混合調(diào)制[J]. 電力自動化設備, 2012, 32(4): 63?66. HE Daqing, CAI Xu. Limitation control and hybrid modulation of modularized multilevel converter[J]. Electric Power Automation Equipment, 2012, 32(4): 63?66.

[7] Glinka M, Marquardt R. A new AC/AC multilevel converter family[J]. IEEE Transactions on Industrial Electronic, 2005, 52(3): 662?669.

[8] Allebrod S, HamerskI R, Marquardt R. New transformerless,scalable modular multilevel converters for HVDC-transmission[C]//Power Electronics Specialists Conference. Chicago, USA, 2008: 174?179.

[9] 屠卿瑞, 徐政, 鄭翔, 等. 模塊化多電平換流器型直流輸電內(nèi)部環(huán)流機理分析[J]. 高電壓技術(shù), 2010, 36(2): 547?552. TU Qingrui, XU Zheng, ZHENG Xiang, et al. Mechanism analysis on the circulating current in modular multilevel converter based HVDC[J]. High Voltage Engineering, 2010, 36(2): 547?552.

[10] 丁冠軍, 湯廣福, 丁明, 等. 新型多電平電壓源換流器模塊的拓撲機制與調(diào)制策略[J]. 中國電機工程學報, 2009, 29(36): 1?6. DING Guanjun, TANG Guangfu, DING Ming, et al. Topology mechanism and modulation scheme of a new multilevel voltage source converter modular[J]. Proceedings of the CSEE, 2009, 29(36): 1?6.

[11] Hagiwara M, Akagi H. Control and experiment of pulse width-modulated modular multilever converters[J]. IEEE Trans on Power Electronics, 2009, 24(7): 1737?1746.

[12] 梁帥奇, 牟曉春, 趙雪, 等. 含有儲能單元的微電網(wǎng)運行控制技術(shù)[J]. 電力科學與技術(shù)學報, 2011, 26(4): 74?79. LIANG Shuaiqi, MOU Xiaochun, ZHAO Xue, et al. Operation and control techniques of micro-grids with energy storage[J]. Journal of Electric Power Science and Technology, 2011, 26(4): 74?79.

[13] Hagiwara M, Maeda R, Akagi H. Control and analysis of the modular multilevel cascade converter based on double-star chopper-cells(MMCC-DSCC)[J]. IEEE Transactions on Power Electronics, 2011, 26(6): 1649?1658.

[14] Hagiwara M, Maeda R, Akagi H. Negative-sequence reactive-power control by the modular multilevel cascade converter based on double-star chopper-cells (MMCC- DSCC)[C]//IEEE Energy Conversion Congress and Exposition (ECCE). Atlanta, USA: IEEE, 2010: 3949?3954.

[15] Marquardt R, Lesnicar A, Hildinger J.Modulares stromrichterkonzept for netzkupplung san-wendung bei hohen spannungen[C]//Proceeding of ETG-Fachtagung. Bad Nauheim, Germany: IEEE, 2002: 1?7.

[16] Qin J, Saeedifard M. Predictive control of a modular multilevel converter for a back-to-back HVDC system[J]. IEEE Transactions on Power Delivery, 2012, 27(3): 1538?1547.

[17] Cortes P, Kouro R, Rocca B L, et al. Guidelines for weighting factors design in model predictive control of power converters and drives[C]//IEEE International Conference on Industrial Technology. Gippsland, Australia: IEEE, 2009: 1?7.

Model predictive control of modular multilevel converter for offshore wind power transmission

ZHOU Yun1, XIA Xiangyang1, KONG Xiangji1, FAN Jiaming1, ZHAO Xiangtao2

(1. College of Electrical and Information Engineering,Changsha University of Science and Technology, Changsha 410077, China; 2. Jining Power Supply Company, State Grid Shandong Electric Power Company, Jining 272000, China)

For the circulating currents suppression of modular multilevel converter (MMC) and the equalizing voltages of sub-module capacitor of flexible high voltage direct current(HVDC) system for offshore wind, a model predictive control (MPC) strategy was proposed, the corresponding prediction model was developed according to discrete-time mathematical model of the system, and the best switching state in each MMC unit was obtained by solving an objective function optimization problem to eliminate the circulating currents and achieve the voltage balancing task of DC capacitor. Finally, MATLAB/SIMULINK was used to simulate the control strategy. The results show that the model predictive control strategy is effective and feasible, its implementation principle iseasy and simple, and it is suitable for MMC-HVDC system of offshore wind power transmission.

modular multilevel converter (MMC); model predictive control(MPC);capacitor voltage balancing; circulating current control; offshore wind power

10.11817/j.issn.1672-7207.2015.10.021

TM76

A

1672?7207(2015)10?3703?07

2014?10?10;

2014?12?23

國家自然科學基金資助項目(51207048);湖南省科技計劃項目(2014GK3006);長沙市科技計劃項目(K1501012-11)(Project (51207048) supported by the National Natural Science Foundation of China; Project (2014GK3006) supported by Hunan Science and Technology Program; Project (K1501012-11) supported by Changsha Science and Technology Program)

夏向陽,博士,教授,從事電力電子技術(shù)在電力系統(tǒng)的應用研究;E-mail:xia_xy@126.com

(編輯 陳燦華)

猜你喜歡
交流模型
一半模型
如此交流,太暖!
科教新報(2022年12期)2022-05-23 06:34:16
加強交流溝通 相互學習借鑒
重要模型『一線三等角』
周櫻 兩岸交流需要更多“對畫”
海峽姐妹(2020年8期)2020-08-25 09:30:18
重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
交流&イベント
交流&イベント
主站蜘蛛池模板: 国产乱子伦一区二区=| 成·人免费午夜无码视频在线观看| 久久精品国产亚洲AV忘忧草18| 亚洲色婷婷一区二区| 欧美日韩理论| 亚洲欧美日韩中文字幕在线| 精品国产99久久| 欧美国产精品不卡在线观看| 一级毛片免费观看不卡视频| 色亚洲成人| 97se亚洲| 香蕉网久久| 99久久无色码中文字幕| 成年看免费观看视频拍拍| 2021亚洲精品不卡a| 97超碰精品成人国产| 一级黄色网站在线免费看| 国产成人1024精品| 国产欧美精品专区一区二区| 国产精品亚洲综合久久小说| 久久人搡人人玩人妻精品一| 久久精品丝袜高跟鞋| 毛片网站在线播放| 尤物特级无码毛片免费| 欧美日韩福利| 免费无遮挡AV| 欧美有码在线观看| 欧美一级专区免费大片| 91po国产在线精品免费观看| 四虎永久免费地址| 国产免费精彩视频| 手机成人午夜在线视频| 在线视频亚洲欧美| 欧美一级99在线观看国产| 亚洲日韩图片专区第1页| 激情成人综合网| 国内精自线i品一区202| 欧美激情视频二区| 色综合a怡红院怡红院首页| 亚洲综合天堂网| 欧美另类视频一区二区三区| 国产导航在线| 日本高清免费一本在线观看 | 亚洲娇小与黑人巨大交| 亚洲码一区二区三区| 激情综合激情| 中文字幕在线观看日本| 四虎影视8848永久精品| 欧美福利在线观看| 欧美无遮挡国产欧美另类| 久久精品66| 伊人久久大香线蕉综合影视| 亚洲最大看欧美片网站地址| 五月婷婷伊人网| 欧美三级视频在线播放| 亚洲精品桃花岛av在线| 亚洲成在人线av品善网好看| 国产十八禁在线观看免费| 欧美一级夜夜爽www| 四虎永久在线视频| 亚洲中字无码AV电影在线观看| 国产资源站| 福利视频一区| 国产毛片网站| 波多野结衣视频网站| 色婷婷成人网| 日韩欧美国产综合| 亚洲婷婷在线视频| 无码国内精品人妻少妇蜜桃视频| 人人澡人人爽欧美一区| 免费看a级毛片| 极品av一区二区| 亚洲女人在线| 99热在线只有精品| 男女猛烈无遮挡午夜视频| 日本91在线| 国产香蕉国产精品偷在线观看| www精品久久| 免费A∨中文乱码专区| 一级毛片免费的| 亚洲天堂网在线播放| 97青草最新免费精品视频|