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

基于改進快速排序算法的MMC均壓控制策略

2017-07-10 10:26:20焦曉鵬劉青
電力建設 2017年7期
關鍵詞:排序方法

焦曉鵬,劉青

(華北電力大學電氣與電子工程學院,河北省保定市 071003)

基于改進快速排序算法的MMC均壓控制策略

焦曉鵬,劉青

(華北電力大學電氣與電子工程學院,河北省保定市 071003)

模塊化多電平換流器(modular multilevel converter,MMC)由于其自身具有輸出電平數高、開關頻率低、波形質量好等優勢而被廣泛研究和使用。子模塊電容的電壓均衡問題是MMC的重點研究方向之一。傳統均壓方法隨著子模塊數目增加,將極大增加開關元件頻率損耗和控制器運算量。該文提出了一種基于改進快速排序算法的均壓策略,通過實時監測子模塊電容電壓,設置子模塊電壓間的離散度指標,繼而控制排序模塊的開通和保持。同時,通過排序算法的優化,使控制器在多模塊時計算效率大大增加,降低硬件要求。在PSCAD/EMTDC仿真平臺搭建MMC-HVDC模型進行仿真驗證。仿真結果表明,改進的均壓控制方法能夠在維持系統特性相對均衡的同時,有效提高運行速度,降低子模塊開關頻率。

模塊化多電平換流器(MMC);高壓直流輸電;電容電壓平衡;開關頻率;快速排序

0 引 言

全控型電力電子元件的快速發展,使得對電壓源型換流器(voltage-sourced converter,VSC)的研究越來越深入。模塊化多電平換流器(modular multilevel converter,MMC)是一種新型的電壓變化電路,其通過多個子模塊級聯的方式,疊加輸出更高的電壓,具有開關頻率低、電壓諧波小、易于模塊化設計等優勢。因此,MMC技術在電力系統中具有廣泛的應用前景[1-2]。基于MMC的諸多優點,國內外學者對其進行了廣泛的研究。

目前,針對MMC的研究多集中在拓撲結構和基本原理、元件參數選擇、脈寬調制策略、環流和電容電壓平衡等方面。MMC不同于傳統三電平拓撲結構,MMC沒有集中的直流側電容,而用每個子模塊中分散的直流電容進行替代。因此,電容的參數差異、損耗、充放電時間等會造成電容電壓差異,使模塊間電壓不平衡,危害換流器的正常運行。所以,橋臂內子模塊電容電壓平衡是各子模塊電壓處于相同水平,各功率器件承受相同應力的重要保證。目前,MMC常采取的均壓方法可以分為2類:(1)通過實時采集各子模塊的電容電壓,應用算法進行排序,然后結合上下橋臂電流參數選擇投切相應子模塊;(2)采用載波移相調制技術,添加輔助電容電壓平衡的控制量,實現電容電壓平衡的效果[3]。

針對MMC子模塊電容電壓問題,文獻[4]提出以橋臂電流為標志,軟件監測各模塊電容電壓,在電流為正時投入電壓較低的若干模塊,電流為負時,投入電壓相對較高的若干模塊,原理簡單,是普遍采用的方法,稱之為傳統方法。但是,當子模塊數量提高時,大量的子模塊電容電壓數據會使得控制器計算量過大,并且開關頻率過高,對硬件要求極高。文獻[5]提出了設定保持因子的電容電壓平衡控制方法,以保持因子乘以電壓未越限的子模塊的電壓值,然后將經過處理的子模塊電壓進行排序,選擇性投切相應模塊。但是,這種方法對于保持因子的選取原則并未研究,僅進行了簡單定量選擇。文獻[6]提出變頻子模塊電容電壓排序方法,通過引入標準差來反映各子模塊電壓均衡程度,根據電容電壓差異度自動調整排序頻率。但是排序過程依然采用傳統排序方法,當模塊數量龐大時,控制器運算量過大,效率較低。文獻[7]提出將子模塊進行分組,并利用質因子分解方法優化分組,以減小控制器運算量,提高計算效率。文獻[8]對文獻[7]的方法進行改進,將質因子分解與希爾排序算法相結合,從而提高排序效率并降低開關頻率,但是算法較為復雜并且隨著分組層數的增加將使混合算法優化效率大大降低。

基于上述問題,為提高系統動態響應速度,減少模塊過多排序對控制器造成的負擔。本文首先提出一種均衡排序控制策略,通過設定一個離散度反映模塊間電容電壓的差異情況,以離散度指標的變化來決定排序模塊啟停,從而實現變頻率排序調節,降低開關頻率損耗。同時,引入效率較高的改進快速排序算法替代傳統冒泡排序,提高計算效率,從而降低控制器運算量,保障了系統高效可靠運行,并能在較低開關頻率下實現電壓平衡控制。最后在PSCAD/EMTDC上搭建MMC模型,驗證本方法的有效性。

1 MMC電容電壓平衡原理

模塊化多電平換流器拓撲結構如圖1所示。

圖1 MMC拓撲結構Fig.1 MMC topology

MMC有6個橋臂,每個橋臂都是由n個完全相同的子模塊相互連接,并且與一個換流電抗器L串聯而成,每個子模塊(sub module,SM)包含2個絕緣柵雙極型晶閘管(insulated gate bipolar transistor,IGBT)和續流二極管組以及1個儲能電容[5]。每個相單元由其所在的整個橋臂構成,總的直流端電壓則由3 個構成相同的相單元并聯提供。當給子模塊上部IGBT一個觸發開通脈沖,下部IGBT一個關斷信號時,子模塊電容將接入橋臂,稱投入該模塊;反之,子模塊電容將從橋臂切除,輸出電壓為零(忽略開關器件的開通壓降),稱切除該模塊[9]。

所以,每個相單元中的投入子模塊數量必須是n個,而通過調整這n個子模塊在上下橋臂的數量比例可以實現不同電平的輸出,這樣才能保證直流電壓的穩定性。因此,子模塊電容電壓的均衡是維持直流電壓穩定的關鍵。

理想情況下,以a相為例,上下橋臂電流為

(1)

式中:I1為上橋臂電流;I2為下橋臂電流;Id為直流端電流;Ia為a相換流器交流側電流。所以,傳統子模塊電容電壓平衡方法是根據某時刻橋臂電流方向及應開通模塊數,確定投入子模塊。當橋臂電流為正方向時,子模塊電容將進行充電,此刻應投入電容電壓相對較小的子模塊;當橋臂電流方向為負方向時,子模塊電容處于放電狀態,此時應投入電容電壓相對較大的子模塊。每一次方向改變時,都需要對子模塊電容電壓進行排序。在電容電壓算法的作用下,使得子模塊電容電壓之間差異很小,從而實現直流電壓的穩定。

MMC的調制方法有很多種,大體分為計算子模塊個數脈寬調制(submodules unified PWM,SUPWM)和載波相移脈寬調制(carrier phase shifted SPWM,CPS-SPWM)2種[10]。本文采取的是最近電平調制(nearest level modulation,NLM)。

最近電平調制即以調制波um與子模塊電容電壓平均值Uc比值最接近的整數來確定電平數,設某時刻a相上橋臂與下橋臂投入的子模塊個數分別為naT和naB。

(2)

式中:N為單個橋臂中上(下)橋臂的子模塊數;round( )為取整函數,取括號內數字最接近的整數。由此可知,任一時刻,NLM方法可以保證naT+naB=N[11]。換流器橋臂控制層接受到由最近電平控制得到的投入子模塊數目后,采用電容電壓均衡控制策略,選取合適數目的子模塊投入。

2 改進的子模塊電容電壓均衡方法

2.1 子模塊電容電壓離散度

傳統子模塊電容電壓均衡方法一般是實時監測子模塊電容電壓并排序,根據橋臂電流方向判斷充放電狀態進行投切子模塊。因此,通過改變子模塊電容充放電和旁路時長即可達到子模塊電容電壓調整效果,此原理即電容電壓排序均衡算法,由于排序頻率恒定,有文獻定義為恒頻排序均衡算法(constant frequency sorting,CFS)[8]。在這種情況下,子模塊的投切只考慮排序結果進行改變,并未考慮開關頻率的損耗,當有較多的模塊數時,減少開關頻率損耗就具有一定意義。

MMC電容電壓均衡的目的是使每個子模塊的電壓波動在工程允許范圍之內,保持近似一致并非要求每個模塊電壓完全一樣。因此,可以通過設置閾值,當子模塊電壓在一定允許范圍內波動時,閉鎖排序模塊,從而保證輸出符合要求的波形的情況下降低開關頻率損耗。有學者利用子模塊電壓的均值作為判斷閾值,可以一定程度降低開關頻率損耗[12-13],但是平均值受極值影響較大,并不能顯著降低開關頻率。

為表征模塊電容電壓均衡程度,引入方差概念。方差是用來度量隨機變量和其數學期望(即均值)之間的偏離程度。統計中的方差(樣本方差)是各個數據分別與其平均數之差的平方的和的平均數。方差可以用來衡量數據間的離散程度。

首先定義某相橋臂子模塊電容電壓為un(n=1,2,…,N),則方差

(3)

式中:ue為電容電壓平均值;uk為第k個子模塊的電容電壓。

D(un)即可反應模塊間電容電壓的分散程度。通過設定D(un)的給定值Dset來進行是否排序的控制。定義離散度:

(4)

當某橋臂子模塊電容電壓值的離散度在設定的允許范圍內時,閉鎖排序模塊,子模塊投入方式不發生變化;離散度超出設定允許范圍時,解鎖排序控制器模塊,對子模塊電容電壓進行重新排序。流程如圖2所示。

圖2 離散度控制流程Fig.2 Flow chart of discrete degree control

通過設定方差閾值,可以在電容電壓變化離散較小的情況下維持原有的觸發脈沖,不進行排序,從而減小開關頻率,降低損耗。

2.2 子模塊電容電壓排序算法優化

子模塊電容電壓的排序算法,關系著控制器運算的時間和效率。當子模塊數量增加時,如果排序算法不夠好,將增加控制器運算量,降低系統效率[14]。

快速排序算法的基本思想是:通過一趟排序將待排序的數列劃分為2個部分,其中一部分的數據都要比另外一部分的所有數據大(或者不小于),然后以該方法對這2個部分數據分別進行排序,通過不斷遞歸進行排序,使整個數據變成有序序列。快速排序算法的時間復雜度在一般情形下可以認為近似等于O(nlog2n)[15]。

快速排序的核心在于能在任何情形下做出最佳的數列分塊,保證分出的 2 個子序列具有大體相似的長度,這樣才可以保證算法的執行時間穩定在O(nlog2n)[15]。

所以,對于快速排序算法的改進,關鍵在于分割數列數據的選取。當使用合適數據作為劃分依據時,可以保證2個部分子數列相對均勻。

假設有n個點x1,x2,…,xn,這n個點與某一分割點x分布的標準差為

(5)

當σ取最小值時,此時的x可以滿足對于點集x1,x2,…,xn,所產生的2個子集A={xi|1≤i≤n}和B={xj|1≤j≤n}滿足A、B無交集,即切分后的點集總體均勻分布在x兩側,x即可作為最佳分割點。

g=(x1-x)2+(x2-x)2+…+ (xn-x)2

(6)

則整理可得

(7)

改進快速排序算法示意圖如圖3所示,假設目標排序序列長度為L,則取該數列數據的平均值X作為比較值temp,設定2個指針i,j對應首末端數據。

圖3 改進快速排序算法示意圖Fig.3 Schematic chart of improved quick sort algorithm

依次比較j指針所指數據與temp的大小,如果大于temp則指針向前移動1次,直到小于temp時,交換i,j所指數據的位置。

依次比較i指針所指數據與temp值的大小,當小于temp時指針向后移動1次,直到大于temp時,交換i,j所指數據的位置。

當i與j相遇時,則完成1次排序過程。然后,對temp之前和之后的數列元素再次遞歸排序,直到分組內元素個數為1時,結束快速排序。

由于改進快速排序算法中需要采用平均值劃分,設對n個排序數求均值的時間為tavg(n),對n個排序數進行1次均值劃分所需時間為f(n),所以改進排序算法的總體時間代價為

T(n)=tavg(n)+f(n)=tavg(n)+n+2f(n/2)=

tavg(n)+n+2[tavg(n/2)+n/2+2f(n/4)]=…=

ktavg(n)+kn+2kf(n/2k)=lbn[tavg(n)+n]+nf(1)=

lbn[tavg(n)+n]=O(lbn·n)=O(nlog2n)

(8)

因此,對于改進快速排序算法,每次對數據求均值并未影響算法效率,整體算法時間性能為O(nlog2n)。 同時,由于劃分均勻,孤點數據和有序數據對算法影響都很小。

3 均衡控制模型

由上述分析可知,改進的子模塊電容電壓排序模型如圖4所示。

圖4 子模塊電容電壓控制原理Fig.4 Schematic diagram of sub-module capacitor voltage balancing control system

每個周期內,對子模塊電容電壓進行離散度計算,如果未超出設定范圍,則保持原有脈沖觸發信號不變。如果超出設定范圍,觸發排序模塊,對子模塊電容電壓進行排序,并將結果輸送給脈沖觸發模塊。脈沖觸發模塊根據NLC調制確定開通模塊數,根據橋臂電流方向確定排序前后順序,進而生成觸發脈沖。通過這樣的方式,可以降低不必要的器件開關,降低開關損耗,并且提高控制器運算效率。

4 仿真分析

通過在 PSCAD/EMTDC上搭建兩端 21電平MMC-HVDC 系統模型,對所提的改進快速排序算法控制的MMC子模塊電容電壓均衡模型的正確性進行驗證。在PSCAD/EMTDC中設置仿真時間為5 s,仿真步長為50 μs。設定單個橋臂上的子模塊數為20個,離散度設置為2%。表 1 給出了模型中具體參數設置。

表1 MMC-HVDC系統參數
Table 1 MMC-HVDC system parameters

4.1 改進排序算法效率仿真

為驗證改進快速排序算法的電壓均衡策略的有效性和準確性,在Viscual C++ 6.0環境下對快速排序算法和改進快速排序算法進行仿真對比。針對系統生成的1 000個隨機數進行排序,結果如表2所示。

表2 不同算法下的排序效率
Table 2 Sorting efficiency under different algorithms

由表2可知,改進快速排序方法與傳統快速排序方法相比,由于均值劃分使數據前后長度大體一致,使得比較次數和交換次數都大大減小,運算效率提高了近17%。

4.2 子模塊開關頻率仿真

圖5表示不同均壓控制方法下a相橋臂的上橋臂子模塊工作狀態,1表示投入,0表示切出。圖5(a)是采用傳統控制方法,圖5(b)采用文獻[13]提出的優化偏差控制方法,圖5(c)使用本文提出的離散度控制方法。由圖5可以看出,采用本文所提出的改進快速排序算法的子模塊均衡策略能夠顯著降低子模塊的開關頻率,效果優于前2種方法,從而可以減小過高開關頻率帶來的損耗。

為了更直觀看出離散度與平衡效果的關系,本文從開關頻率和子模塊電容電壓最大波動幅度(與子模塊電壓額定值相比)2個方面對不同離散度取值進行比較,引入平均開關頻率favg[16]:

(9)

式中:nj為橋臂中子模塊開關次數;n為單個上(下)橋臂模塊數。

圖5 子模塊工作狀態Fig.5 Working condition of sub module

表3給出了不同離散度設置時,子模塊平均開關頻率和電壓最大偏離值百分比的變化情況。其中,子模塊電壓基準值上下限取66 kV和64 kV。由表3可知,當離散度為0時,即不設置控制閾值,開關頻率為570 Hz,隨著離散度的設置,開關頻率開始顯著降低,而子模塊電壓最大波動幅值有所增加。當離散度為0.04時,開關頻率降到92 Hz,而電壓波動幅值達到6.1%。由于實際工程中,MMC正常工作要求模塊電壓偏差不超過5%[17],因此離散度并非越大越好,根據表3的結果,本文取離散度為0.02,較為適宜。

表3 不同離散度下子模塊相關特性變化
Table 3 Performance change of sub-module with different dispersion degrees

圖6為改進均壓方法下(離散度取0.02)的子模塊電容電壓情況,可以看出,改進后電容電壓波動并未受到太大影響,仍在工程允許的范圍內。

圖6 子模塊電容電壓仿真Fig.6 Simulated waveform of sub-module capacitor voltage

4.3 均衡控制下的系統特性仿真

圖7為傳統均壓和改進均壓措施下系統特性的對比,由圖7可知,采用改進均壓措施后系統的閥側交流電壓、直流側電壓與傳統方法相比并未發生較大變化,系統特性較為穩定。

圖7 改進電容均壓控制下的系統特性Fig.7 System performance with improved voltage balancing control

5 結 論

(1)本文深入研究了MMC-HVDC子模塊電容電壓均衡方法,提出了一種適用于多模塊的基于改進快速排序算法的均衡控制策略。該方法通過設置離散度閾值,有效降低了子模塊開關頻率,從而降低系統損耗,提高控制器運算速度。

(2)通過將改進快速排序算法應用到子模塊的電壓排序中,可以將排序時間效率控制在nlog2n數量級,從而提高控制器運算效率,減少控制器運算時間。

(3)仿真結果表明,采用改進算法的MMC子模塊電容均壓模型與傳統均壓模型相比,MMC系統特性效果基本相同,對HVDC系統的運行特性影響很小。

[1]丁冠軍,湯廣福,丁明,等.新型多電平電壓源換流器模塊的拓撲機制與調制策略[J].中國電機工程學報,2009,29(36):1-8. 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-8.

[2]劉欣和,吳金龍,韓坤,等.基于電壓向量合成的模塊化多電平換流器控制策略[J].電力建設,2014,35(3):36-42. LIU Xinhe, WU Jinlong, HAN Kun, et al. Research of modular multi-level converter control strategy based on voltage vector composition[J]. Electric Power Construction, 2014,35(3):36-42.

[3]辛業春,王朝斌,李國慶,等.模塊化多電平換流器子模塊電容電壓平衡改進控制方法[J].電網技術,2014,38(5):1291-1299. XIN Yechun, WANG Chaobin, LI Guoqing, et al. An improved balance control for sub-module capacitor voltage of modular multilevel converter[J]. Power System Technology, 2014,38(5):1291-1299.

[4]丁冠軍,丁明,湯廣福,等.新型多電平 VSC子模塊電容參數與均壓策略[J].中國電機工程學報,2009,29(30):1-6. DING Guanjun, DING Ming, TANG Guangfu, et al. Submodule capacitance parameter and voltage balancing scheme of a new multilevel VSC modular[J]. Proceedings of the CSEE,2009,29(30):1-6.

[5]管敏淵,徐政.MMC 型 VSC-HVDC 系統電容電壓的優化平衡控制[J].中國電機工程學報,2011,31(12):9-14. GUAN Minyuan, XU Zheng. Optimized capacitor voltage balancing control for modular multilevel converter based VSC-HVDC system[J]. Proceedings of the CSEE,2011,31(12):9-14.

[6]張建坡,趙成勇,孫海峰.基于改進拓撲的 MMC-HVDC 控制策略仿真[J].中國電機工程學報,2015,35(5):1032-1038. ZHANG Jianpo, ZHAO Chengyong, SUN Haifeng. Simulations of control strategies in MMC-HVDC with an improved topology[J]. Proceedings of the CSEE,2015,35(5):1032-1038.

[7]彭茂蘭,趙成勇,劉興華,等.采用質因子分解法的模塊化多電平換流器電容電壓平衡優化算法[J].中國電機工程學報,2014,34(33):5847-5853. PENG Maolan, ZHAO Chengyong, LIU Xinghua, et al. An optimized capacitor voltage balancing control algorithm for modular multilevel converter employing prime factorization method[J]. Proceedings of the CSEE,2014,34(33):5847-5853.

[8]何智鵬,許建中,苑賓,等.采用質因子分解法與希爾排序算法的MMC電容均壓策略[J].中國電機工程學報,2015,35(12):5842-5849. HE Zhipeng, XU Jianzhong, YUAN Bin, et al. A capacitor voltage balancing strategy adopting prime factorization method and shell sorting algorithm for modular multilevel converter[J]. Proceedings of the CSEE,2015,35(12):5842-5849.

[9]徐政.柔性直流輸電系統[M]. 北京: 機械工業出版社,2013:75-80.

[10]KONSTANTINOU G, CIOBOTARU M, AGELIDIS V.Selective harmonic elimination pulse-width modulation of modular multilevel converters[J].IET Power Electronics,2013,6(1):96-107.

[11]ANTONIOS A,LENNART A,HANS-PETER N.On dynamics and voltage control of the modular multilevel converter[C]//Power Electronics and Applications, 2009. EPE′09. 13th European Conference on.Barcelona,Spain:IEEE,2009:1-10.

[12]張興輝,武宏波,王朝亮.基于自抗擾控制技術的MMC環流抑制器[J].電力建設,2014,35(11):97-101. ZHANG Xinghui, WU Hongbo, WANG Chaoliang. Circulating current suppressing controller for MMC based on active disturbance rejection control technology[J]. Electric Power Costruction,2014,35(11):97-101.

[13]付華,于翔.希爾排序算法優化偏差控制的MMC電容均壓模型[J].電源學報,2016,12(11):21-29. FU Hua, YU Xiang. The capacitor voltage balancing strategy adopting bias control and shell sorting algorithm for modular multilevel converter[J].Journal of Power Supply, 2016,12(11):21-29.

[14]許建中.模塊化多電平換流器電磁暫態高效建模方法研究[D]. 北京:華北電力大學,2014. XU Jianzhong. Research on the electromagnetic transient efficcient modelling method of modular multilevel converter[D]. Beijing: North China Electric Power University, 2014.

[15]湯亞玲,秦鋒.高效快速排序算法研究[J]. 計算機工程,2011,37(6):77-78. TANG Yaling, QIN Feng. Study of high efficiency quick sort algorithm[J]. Computer Engineering, 2011, 37(6):77-78.

[16]ILVES K,ANTONOPOULOS A,NORRGA S,et al.A new modulation method for the modular multilevel converter allowing fundamental switching frequency[J]. IEEE Transactions on Power Electronics, 2012, 27(8):3482-3494.

[17]湯廣福.基于電壓源換流器的高壓直流輸電術[M].北京:中國電力出版社,2010.

(編輯 張小飛)

Voltage Balancing Control Strategy for Modular Multilevel Converter Based on Improved Quick Sort Algorithm

JIAO Xiaopeng,LIU Qing

(School of Electrical and Electronic Engineering, North China Electric Power University, Baoding 071003, Hebei Province, China)

Modular multilevel converter (MMC) has a widely application and research because of its advantages of high output level, low switching frequency and good waveform quality, etc. The voltage balance of sub-module capacitor is one of the key research directions of MMC. While traditional voltage balancing method will greatly increase the frequency loss of switching element and controller operation with the increasing of the number of sub-modules. This paper puts forward an optimized capacitor voltage balancing strategy for MMC based on improved quick sort algorithm, which can trigger or maintain control scheduling module by real-time monitoring the voltage of sub-module capacitor, setting the discrete degree index of sub-modules voltage. At the same time, through the optimization of the sorting algorithm, the calculation efficiency of controller in multi-modules is greatly increased and the hardware requirements is reduced. Finally, this paper builds the MMC-HVDC model in PSCAD/EMTDC simulation platform to simulate and verify the algorithm. The simulation results show that the improved voltage balancing control method can significantly improve the running speed and reduce the switching frequency of sub-module without noticeablly influencing the system performance.

modular multilevel converter (MMC); high-voltage direct current; capacitor voltage balancing; switching frequency; quick sort

TM 46

A

1000-7229(2017)07-0146-07

10.3969/j.issn.1000-7229.2017.07.018

2017-03-15

焦曉鵬 (1991),男,碩士研究生,研究方向為柔性交流輸電和智能電網;

劉青(1974),女,博士,副教授,主要研究方向為電力系統繼電保護、電力系統安全防御與恢復控制。

猜你喜歡
排序方法
排排序
排序不等式
恐怖排序
學習方法
節日排序
刻舟求劍
兒童繪本(2018年5期)2018-04-12 16:45:32
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
賺錢方法
捕魚
主站蜘蛛池模板: 在线播放国产99re| 亚洲精品天堂自在久久77| 香蕉eeww99国产精选播放| 国产亚洲高清在线精品99| 一区二区三区成人| 欧美中文字幕无线码视频| 欧美亚洲第一页| 国产精品视频a| 高h视频在线| 国产精品自在拍首页视频8| 国产成人久久综合777777麻豆| 中文字幕日韩丝袜一区| 青青青国产视频| 成人综合网址| 国产精品不卡片视频免费观看| 激情無極限的亚洲一区免费| 青草视频免费在线观看| a毛片在线| 久久伊人色| 日韩麻豆小视频| 亚洲香蕉久久| 亚洲天堂网在线播放| 亚洲国产91人成在线| 毛片网站观看| 婷婷亚洲天堂| 美女被操黄色视频网站| 香港一级毛片免费看| 国产视频入口| 性欧美在线| 91久久国产热精品免费| 亚洲福利网址| 国产第八页| 亚洲国产av无码综合原创国产| 91九色视频网| 亚国产欧美在线人成| 欧美精品H在线播放| 久久成人18免费| 欧美日本在线一区二区三区| 美女视频黄频a免费高清不卡| 久草国产在线观看| 久久动漫精品| 国产成人凹凸视频在线| 日本影院一区| 91蜜芽尤物福利在线观看| 国产一区亚洲一区| 在线另类稀缺国产呦| 久久免费看片| 久久夜色精品国产嚕嚕亚洲av| 国产成人精品午夜视频'| 九九久久精品免费观看| 久久精品波多野结衣| 色婷婷成人网| av尤物免费在线观看| 欧美日韩国产成人高清视频| 欧美色亚洲| 婷婷综合色| 一级黄色欧美| 国产尤物在线播放| 欧洲亚洲一区| 国产爽歪歪免费视频在线观看| 亚洲国产日韩一区| 日韩欧美成人高清在线观看| 国产亚洲精品97AA片在线播放| 亚洲αv毛片| 亚洲国内精品自在自线官| 一级看片免费视频| 欧美国产日韩在线播放| 精品无码视频在线观看| 国产精品亚洲一区二区在线观看| 亚洲国产成人在线| 国产男人的天堂| 国产精品成人第一区| 少妇被粗大的猛烈进出免费视频| 91破解版在线亚洲| 国产一区二区精品福利| 特级毛片免费视频| 999精品在线视频| 欧美精品高清| 亚洲精品无码AⅤ片青青在线观看| 欧美日韩亚洲国产主播第一区| 国内丰满少妇猛烈精品播| 人妻免费无码不卡视频|