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

考慮低頻振蕩的MMC有源阻尼環流抑制方法

2015-10-25 05:52:50徐千鳴馬伏軍熊橋坡王凌志
電工技術學報 2015年24期
關鍵詞:系統

徐千鳴 羅 安 馬伏軍 熊橋坡 王凌志

(1. 湖南大學國家電能變換與控制工程技術研究中心 長沙 410082

2. 國網湖南省電力公司長沙供電公司 長沙 410000)

考慮低頻振蕩的MMC有源阻尼環流抑制方法

徐千鳴1羅安1馬伏軍1熊橋坡1王凌志2

(1. 湖南大學國家電能變換與控制工程技術研究中心長沙410082

2. 國網湖南省電力公司長沙供電公司長沙410000)

以直接調制方式的開關函數為出發點,建立MMC每相橋臂電容電壓之和與橋臂環流的暫態數學模型,揭示了系統在異常工況下會產生低頻振蕩環流,推導低頻振蕩環流的頻率,指出無環流控制時,二倍頻環流存在諧振點。為實現對低頻振蕩與二倍頻環流的同時抑制,提出一種基于有源阻尼控制的環流抑制方法,簡單實用,無需相間解耦與坐標變換。仿真與實驗結果表明該方法能同時抑制對暫態時的低頻振蕩與穩態時的二倍頻及以上頻次環流,可增強系統穩定性。

模塊化多電平低頻振蕩環流抑制有源阻尼

0 引言

由于具有模塊化易擴展、輸出電壓畸變率低、公共直流母線和無需工頻變壓器等優點,模塊化多電平換流器(Modular Multilevel Converter,MMC)自21世紀初被Rainer Marquardt等提出以來,在高壓直流輸電、電能質量治理以及中高壓變頻等領域得到了廣泛的應用[1-6]。但MMC在工程應用中也存在子模塊電容電壓不平衡[7]、環流[8]等問題。

MMC的環流是指在直流母線正負極之間或者不同相之間流動的電流。MMC的分布式儲能電容與橋臂電感相串聯的特殊結構,使其三相之間以及與外部電網之間容易發生諧振,進而產生低頻振蕩環流;實際運行中,MMC的子模塊電壓的波動隨調制反映在直流母線電壓與輸出電壓上,導致橋臂電流中存在偶次諧波環流,其中高次諧波環流相比于二倍頻環流很小,可以忽略;同時輸出電流中含有3次頻為主的奇次諧波電流分量,3次及其倍數次諧波電流會被三角形聯結的變壓器阻斷。

現有文獻對MMC環流的研究主要集中在二倍頻環流產生機理及其抑制方法,對低頻振蕩環流的機理分析及抑制方法少有報道。文獻[9]通過分析二倍頻環流的影響因素,指出增大橋臂電抗可以減小二倍頻環流,但這不僅增大了換流器成本與損耗,還會降低系統的響應速度;文獻[10]通過外加LC濾波電路構成LCL濾波支路,以抑制橋臂的二倍頻環流,但引入的LC支路會導致橋臂阻抗存在串聯諧振問題;文獻[11]建立MMC動態模型,分析了MMC存在的固有低頻振蕩頻率及抑制方法,但沒有考慮對二倍頻環流的影響;文獻[12]詳細研究了MMC內部環流和輸出電流諧波成分,并給出了二倍頻環流的諧振頻率;文獻[13,14]提出基于坐標變換的二倍頻環流抑制方法,但在電網不平衡時需用到多個同步旋轉坐標;文獻[15,16]提出了無需坐標變換的PI控制和PR控制的環流抑制策略,但增加了控制系統的難度。

針對上述問題,本文根據直接調制方式的開關函數,建立系統橋臂電容電壓之和與橋臂環流的數學模型,指出系統低頻振蕩環流的原因及其頻率和二倍頻環流諧振問題,提出基于有源阻尼控制的環流抑制方法(Circulating Current Suppressing Method,CCSM),能同時實現低頻振蕩與二倍頻環流的抑制,仿真和實驗結果驗證了本文提出的控制方法的有效性。

1 MMC環流分析

1.1MMC基本原理

單相MMC結構如圖1所示,分為上、下兩個橋臂,上、下橋臂相對于直流母線為串聯結構,而相對于交流電網是并聯結構。每個橋臂均由數目相同或者相近的子模塊與一個橋臂電抗串聯組成,子模塊結構為半H橋。

根據圖1,建立電路KVL與KCL方程

圖1 MMC單相結構示意圖Fig.1 Schematic configuration of single-phase MMC

式中,ixz與ixo分別為x(x=a,b,c)相環流與輸出電流;Ud為公共直流母線電壓;uxp和uxn分別為上、下橋臂的實時投入電壓;R、L和C分別為橋臂等效電阻、橋臂電感和子模塊電容。

分別將上、下橋臂電壓的差模部分與共模部分定義為

聯立式(1)~式(3)可得

式(4)包含MMC系統的交流接口微分方程與直流接口微分方程,描述了MMC系統的外部特性與內部特性,可看出MMC輸出電流與內部環流之間無耦合。

MMC子模塊直流側并不是理想的直流電源,而是有容值限制的電容。交、直流側進行功率傳輸或者相間功率互換時,必然造成電容電壓的波動。假設子模塊電容電壓波動一致,根據橋臂電容電壓變化與流過電流的關系,上、下橋臂N個子模塊電壓之和的波動可以表示為

將式(2)、式(3)、式(5)和式(6)代入式(7),可得描述MMC系統上、下橋臂電容電壓之和與橋臂環流的三元一階微分方程組為

為簡化分析,對微分方程組進行整理得

將式(9)代入式(8),微分方程組變為

式(10)描述了MMC橋臂電容電壓之和及電壓之差與環流的關系,設輸出參考電壓與輸出電流為

則換流器單相平均功率P0、脈沖功率P~與瞬時功率P表示為[17]

式中,Uom為輸出參考電壓幅值;Iom為輸出電流幅值;S為換流器單相視在功率;ω0為基波角頻率;φ為輸出參考電壓與輸出電流的相位差。

1.2MMC低頻振蕩環流

求得系統穩態工作點為

由式(13)的狀態矩陣可以解得系統特征方程根為

式中,α為振蕩的衰減速度;ωlf為振蕩的頻率,一般低于基波頻率,故稱之為低頻,其數值跟模塊個數N、橋臂電抗L、等效電阻R及子模塊電容C相關。低頻振蕩環流的本質是子模塊電容經由橋臂電抗及電阻充放電產生的沖擊電流。直流側電壓或傳輸功率的突變都可能導致低頻環流的發生,文獻[18]針對直流母線短路的工況進行了具體分析,本文不再贅述。

值得注意的是,電感L、電容C及模塊數N選擇不當時,ωlf可能處于MMC系統典型工作頻率如二倍頻和四倍頻等,使系統發生諧振。進行參數L和C設計時,一般要求ωlf低于基頻,忽略橋臂電阻R時,表示為

搭建基于PSIM9.0的三相MMC逆變模型驗證推導的振蕩環流頻率,圖2為ωlf與L、C關系三維圖,可知ωlf隨L和C的減小而增大,當L和C小至一定程度后,ωlf處于工頻及以上頻次的區域。圖3為三電平與五電平模型下,子模塊電容保持不變時,ωlf隨L變化的曲線,仿真結果和理論分析基本一致,說明理論推導的正確。

圖2 低頻振蕩頻率與橋臂電感、子模塊電容關系Fig.2 Relationship between oscillation frequency and arm inductor and sub module capacitance

圖3 低頻振蕩頻率與橋臂電感關系Fig.3 Relationship between oscillation frequency and arm inductor

1.3MMC二倍頻環流

考慮功率傳輸導致的二倍頻電壓與電流波動,將式(11)、式(12)、式(14)和式(15)代入式(10),構造系統小信號模型為

聯立式(20)和式(21),二元一階微分方程組化為一元二階微分方程

二倍頻環流即為微分方程式(22)的特解,解得二倍頻環流幅值為[19]

式中

由式(25)可知,當N、L、C滿足式(26)時,二倍頻環流幅值存在諧振點,其取最大值即式(27),此時二倍頻環流幅值完全由橋臂等效電阻R限制。

為避免二倍頻環流幅值的諧振放大,參數N、L和C選取需滿足約束[20]

2 MMC環流抑制

由式(16)與式(25)可知,增大系統阻尼比能有效抑制系統在非正常工作狀態下的低頻振蕩環流和二倍頻環流。增大系統阻尼比有兩種方法[21]:①增大橋臂的等效電阻R,但會導致換流器損耗的增加和效率的降低;②采用有源阻尼控制,通過控制算法增加系統虛擬電阻,在不增加損耗的情況下,增大系統阻尼比,即

式中,環流中交流分量的獲得需要先提取環流的直流分量,如何獲得橋臂環流中的直流分量成為有源阻尼環流抑制方法控制實現的關鍵。通常直流分量使用低通濾波器來獲取,但低通濾波器存在頻帶窄的缺點,濾除交流分量的同時,會影響系統動態性能,根據式(21),采用陷波器,將陷波角頻率設計為基波頻率的2倍,即可獲得環流中的直流分量,避免了對陷波頻率以外頻次的影響,所設計控制策略如圖4所示。

圖4 有源阻尼環流抑制方法框圖Fig.4 Schematic of the proposed CCSM based on active damping control

2.1低頻振蕩抑制

根據上節所述有源阻尼環流抑制方法,聯立式(10)、式(14)和式(15),忽略引入的基頻及以上頻次波動與高階部分,重構系統小信號模型可得

由狀態矩陣可得系統特征方程根為

為保證系統的穩定性,根據勞斯判據

由式(33)可知,Rv有很大的取值裕度,使系統為過阻尼狀態,即Rv?R,以抑制低頻振蕩,或在振蕩發生時加快其衰減速度,增強系統的穩定性。

圖5為加入CCSM前后,系統的伯德圖,可以看出加入環流抑制后,諧振峰值明顯下降,由L、C引起的低頻振蕩得到有效抑制。

圖5 CCSM加入前后系統伯德圖Fig.5 System Bode diagram of MMC with and without CCSM

2.2二倍頻環流抑制

由2.1節可知,加入CCSM后,能有效抑制MMC在暫態時的低頻振蕩。下面進一步研究加入CCSM對穩態時MMC二倍頻環流的影響。

加入CCSM后,將式(29)代入式(19),忽略高階分量,系統可簡化為

聯立式(20)和式(34),二元一階微分方程組化為一元二階微分方程

二倍頻環流即為式(35)微分方程的特解,解得二倍頻環流幅值為

同樣,二倍頻環流幅值在滿足式(37)時,存在諧振點,取得最大值見式(38)。

結合式(25)和式(36),可以看出加入CCSM后,二倍頻環流雖沒有完全消除,但幅值大幅度減小。

以上分析說明所提出的有源阻尼環流抑制方法不僅可以有效抑制系統暫態時的低頻振蕩環流,還能有效抑制系統穩態時的二倍頻環流。

3 仿真分析

為驗證本文提出的有源阻尼環流抑制方法的有效性,在PSIM9.0仿真環境下,搭建了基于圖1所示的單相四電平MMC無源逆變模型,Rv取為25Ω,仿真參數見下表。

表 仿真參數Tab. Simulation parameters

圖6 CCSM加入前、后橋臂電流、環流仿真波形Fig.6 Simulation results of arm currents and circulating current

圖7 橋臂電流、環流頻譜Fig.7 FFT analysis of arm currents and circulating current

t0時刻帶載起動MMC,子模塊電壓開始上升。由圖6與圖7可知,橋臂環流除含有直流分量外,還存在二倍頻分量與29Hz的低頻振蕩分量。

t1時刻加入有源阻尼環流抑制,低頻振蕩環流與二倍頻環流同時得到抑制,低頻振蕩環流幅值由18.8A降至0.08A,二倍頻環流幅值由9.3A減小為0.25A,而直流分量與基頻分量基本不變。表明所提出的有源阻尼環流控制方法能同時有效抑制系統暫態過程中的低頻振蕩環流與穩態時的二倍頻環流,提升系統暫態與穩態性能。

4 實驗

為驗證本文理論分析及控制策略的效果,在仿真結果正確的基礎上,以圖1所示的主電路結構為基礎,搭建單相四電平MMC樣機進行研究,系統設定為滯后無功發生工況,輸出參考電壓由直流電壓外環與無差拍電流內環控制獲得,環流抑制策略如圖4所示,子模塊電容電壓平衡采用比例控制,調制方式為載波移相調制。MMC交流側通過交流調壓器與電網AB相連接,調壓器一次側與二次側電壓有效值分別410V與190V,指令電流幅值為30A,模塊電容電壓設定為200V,載波頻率為2kHz,虛擬電阻Rv取20Ω,電壓、電流方向與圖1一致,其余參數與表1仿真參數相同。

單相四電平MMC實驗裝置見圖8a,圖8b為穩態時電網電壓eAB(調壓器一次側)與MMC輸出無功電流io波形,可以看出無功電流滯后于電網電壓約90°,有效值為20.47A。圖9a為下橋臂電壓波形,峰-峰值為600V左右,圖9b為下橋臂子模塊6電容電壓波形,在202V上下波動。

圖9 下橋臂電壓與模塊電容電壓波形Fig.9 Waveforms of lower arm voltage and sub module capacitance voltage

在系統已起動后,加入有源阻尼環流抑制,將DSP存儲的數據導入到Origin75工作空間,得到環流在CCSM加入前后的波形如圖10所示,抑制前存在低頻振蕩環流與二倍頻環流且幅值很大,抑制后不含低頻環流,二倍頻環流幅值大幅度減小。圖11為電能質量分析儀Fluke434記錄的環流不控實驗與環流控制實驗電流波形與頻譜的對比。環流不控組:橋臂電流含有較大二倍頻分量,電流畸變率(Total Harmonic Distortion,THD)約為15%,輸出電流THD為3.4%;環流控制組:二倍頻環流大幅度降低,橋臂電流THD降為5%,輸出電流THD降至2.6%,說明所述環流抑制方法的有效性。

圖10 CCSM加入前后環流波形Fig.10 Waveforms of circulating current

圖11 輸出電流、橋臂電流波形與頻譜Fig.11 Output current,arm currents and spectrum

5 結論

本文建立MMC橋臂電容電壓之和與環流的數學模型,得出以下結論:

(1)系統阻尼較小時,暫態過程容易產生低頻振蕩環流,其頻率與主電路參數相關。

(2)無環流控制時,二倍頻環流存在諧振點。

(3)基于有源阻尼控制的環流抑制方法能同時對MMC暫態時的低頻振蕩與穩態時的二倍頻及以上頻次環流進行抑制。

附錄

直流分量在后續計算中,不影響二倍頻分量大小,可以省略,其中調制度m=2UomUd,則

對式(A2)中積分項作以下變換

代入式(A2)后,可得

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

[2] Akagi H. Classification,terminology,and application of the modular multilevel cascade converter (MMCC)[J]. IEEE Transactions on Power Electronics,2011,26(11):3119-3130.

[3] 鄧旭,王東舉,沈揚,等. 舟山多端柔性直流輸電工程換流站內部暫態過電壓[J]. 電力系統保護與控制,2013,41(18): 111-119.

Deng Xu,Wang Dongju,Shen Yang,et al. Research on transient overvoltage for converter station of Zhoushan multi-terminal VSC-HVDC project[J]. Power System Protection and Control,2013,41(18):111-119.

[4] 楊曉峰,范文寶,王曉鵬,等. 基于模塊組合多電平變換器的STATCOM及其控制[J]. 電工技術學報,2011,26(8): 7-13.

Yang Xiaofeng,Fan Wenbao,Wang Xiaopeng,et al. Study on modular multilevel converter based STATCOM and its control[J]. Transactions of China Electrotechnical Society,2011,26(8): 7-13.

[5] 翟曉萌,趙成勇,李路遙,等. 模塊化多電平動態電壓恢復器的研究[J]. 電力系統保護與控制,2013,41(12): 86-91.

Zhai Xiaomeng,Zhao Chengyong,Li Luyao,et al. Study of modular multilevel dynamic voltage restorer[J]. Power System Protection and Control,2013,41(12):86-91.

[6] Baruschka L,Mertens A. A new 3-phase AC/AC modular multilevel converter with six branches in hexagonal configuration[C]. IEEE Energy Conversion Congress and Exposition (ECCE),Phoenix,2011:4005-4012.

[7] 王奎,鄭澤東,李永東,等. 新型模塊化多電平變換器電容電壓波動規律及抑制方法[J]. 電工技術學報,2011,26(5): 8-14.

Wang Kui,Zheng Zedong,Li Yongdong,et al. Voltage ripple principle and restrain method of floating capacitors in a new modular multilevel converter[J]. Transactions of China Electrotechnical Society,2011,26(5): 8-14.

[8] 楊曉峰,王曉鵬,范文寶,等. 模塊組合多電平變換器的環流模型[J]. 電工技術學報,2011,26(5):21-27.

Yang Xiaofeng,Wang Xiaopeng,Fan Wenbao,et al. Circulating current model of modular multilevel converters[J]. Transactions of China Electrotechnical Society,2011,26(5): 21-27.

[9] 屠卿瑞,徐政,鄭翔,等. 模塊化多電平換流器型直流輸電內部環流機理分析[J]. 高電壓技術,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] Ilves K,Norrga S,Harnefors L,et a1. Analysis of arm current harmonics in modular multilevel converters with main-circuit filters[C]. 9th International Multi-Conference on Systems,Signals and Devices(SSD),Chemnitz,2012: 1-6.

[11] Lennart H,Norrga S,Antonopoulos A,et a1. Dynamic modeling of modular multilevel converters[C]. 14th European Conference on Power Electronic and Applications (EPE),2011: 1-10.

[12] Ilves K,Antonopoulos A,Norrga S,et a1. Steadystate analysis of interaction between harmonic components of arm and line quantities of modular multilevel converters[J]. IEEE Transactions on Industrial Electronics,2012,27(1): 57-68.

[13] 卓谷穎,江道灼,連霄壤. 模塊化多電平換流器不平衡環流抑制研究[J]. 電力系統保護與控制,2012,40(24) : 118-124.

Zhuo Guying,Jiang Daozhuo,Lian Xiaorang. Study of unbalanced circular current suppressing for modular multilevel converter[J]. Power System Protection and Control,2012,40(24): 118-124.

[14] Tu Qingrui,Xu Zheng,Zhang Jing. Circulating current suppressing controller in modular multilevel converter[C]. 36th Annual Conference on IEEE Industrial Electronics Society (IECON),2010: 3198-3202.

[15] 楊曉峰,鄭瓊林. 基于MMC環流模型的通用環流抑制策略[J]. 中國電機工程學報,2012,32(18): 59-65.

Yang Xiaofeng,Zheng Qionglin. A novel universal circulating current suppressing strategy based on the MMC circulating current model[J]. Proceedings of the CSEE,2012,32(18): 59-65.

[16] Li Zixin,Wang Ping,Chu Zunfang,et a1. An inner current suppressing method for modular multilevel converters[J]. IEEE Transactions on Industrial Electronics,2013,28(11): 4873-4879.

[17] Brandon J,Pierquet,David J,et al. A single- phase photovoltaic inverter topology with a seriesconnected energy buffer[J]. IEEE Transactions on Industrial Electronics,2013,28(10): 4603-4611.

[18] 王姍姍,周孝信,湯廣福,等. 模塊化多電平換流器HVDC 直流雙極短路子模塊過電流分析[J]. 中國電機工程學報,2011,31(1): 1-7.

Wang Shanshan,Zhou Xiaoxin,Tang Guangfu,et al. Analysis of submodule overcurrent caused by DC pole-to-pole fault in modular multilevel converter HVDC system[J]. Proceedings of the CSEE,2011,31(1): 1-7.

[19] 徐政. 柔性直流輸電系統[M]. 北京: 機械工業出版社,2012.

[20] Zhou Yuebin,Jiang Daozhuo,Guo Jie,et al. Analysis and control of modular multilevel converters under unbalanced conditions[J]. IEEE Transactions on Industrial Electronics,2013,28(4): 1986-1995.

[21] 張興,張崇巍. PWM整流器及其控制[M]. 北京:機械工業出版社,2012.

Circulating Current Suppressing Method Based on Active Damping Control of MMC Considering Low-Frequency Oscillation

Xu Qianming1Luo An1Ma Fujun1Xiong Qiaopo1Wang Lingzhi2
(1. National Electric Power Conversion and Control Engineering Technology Research Center Hunan UniversityChangsha410082China
2. State Grid Hunan Electric Power CompanyChangsha410000China)

In order to analyze and suppress the low-frequency circulating currents of modular multilevel converter (MMC),a transient mathematical model for the sum of capacitor voltages and the circulating current of arm is established. The model reveals the principle of low-frequency oscillation under abnormal conditions. And then,the analytical expressions of oscillation frequency and the resonance point of 2nd circulating current without circulating current control are derived. In order to inhibit both low-frequency oscillation and 2nd circulating current,a circulating current suppressing method (CCSM) based on active damping control is proposed,which is simple without decoupling and coordinate transformation. Finally,the simulation and experimental results show that this method can suppress the low-frequency oscillation in transient state as well as 2nd and 2nd above circulating currents in steady state simultaneously.

Modular multilevel converter,low-frequency oscillation,circulating current suppression,active damping

TM712

徐千鳴男,1989年生,博士研究生,研究方向為輕型直流輸電與電能質量控制。

羅安男,1957年生,教授,博士生導師,研究方向為電能變換與控制、電能質量控制技術與裝備研制。

國家自然科學基金重點項目(51237003)和國家自然科學基金面上項目(51477045)資助。

2014-03-05改稿日期 2014-05-09

猜你喜歡
系統
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
WJ-700無人機系統
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
基于PowerPC+FPGA顯示系統
基于UG的發射箱自動化虛擬裝配系統開發
半沸制皂系統(下)
FAO系統特有功能分析及互聯互通探討
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
一德系統 德行天下
PLC在多段調速系統中的應用
主站蜘蛛池模板: 亚洲天堂福利视频| www.99在线观看| 欧美成人A视频| 国产精品九九视频| 色综合久久无码网| 精品福利一区二区免费视频| 一级毛片在线免费视频| a在线亚洲男人的天堂试看| 国产精品浪潮Av| 亚洲最大情网站在线观看| 中文字幕人妻av一区二区| 日韩精品一区二区三区中文无码| 操美女免费网站| 色综合中文| 国产在线观看91精品| 青青草原国产一区二区| 天天色综网| 国产精品自在在线午夜区app| 亚洲伊人天堂| 不卡国产视频第一页| 亚洲欧美人成电影在线观看| 国模私拍一区二区| 久久精品无码国产一区二区三区| 制服丝袜一区二区三区在线| 伊伊人成亚洲综合人网7777 | 在线免费不卡视频| 狼友视频国产精品首页| 久青草网站| 久久国产成人精品国产成人亚洲| 狂欢视频在线观看不卡| 热久久综合这里只有精品电影| 萌白酱国产一区二区| 3p叠罗汉国产精品久久| 精品国产三级在线观看| 日本一本正道综合久久dvd| 国产精品成人啪精品视频| 911亚洲精品| 第一页亚洲| 97视频精品全国在线观看| 亚洲男人天堂2018| 中文字幕中文字字幕码一二区| 国模视频一区二区| 久久久久国产精品熟女影院| 无码免费试看| 国产精品免费电影| 亚洲精品另类| 欧美精品黑人粗大| 成年人国产视频| 亚洲国产精品久久久久秋霞影院| 国产小视频在线高清播放| 伊人精品成人久久综合| 97av视频在线观看| 狠狠做深爱婷婷久久一区| 色综合天天娱乐综合网| 啊嗯不日本网站| 亚洲国产精品一区二区第一页免 | 欧美日韩国产精品综合| 91视频国产高清| 成年看免费观看视频拍拍| 2018日日摸夜夜添狠狠躁| 2021国产乱人伦在线播放 | 国产黄在线免费观看| 91小视频在线观看免费版高清| 1级黄色毛片| 国产福利一区视频| 精品视频一区在线观看| www.99在线观看| 亚洲性日韩精品一区二区| 亚洲欧美日韩中文字幕在线一区| 中国国产高清免费AV片| 成人在线观看不卡| 国产18页| 亚洲v日韩v欧美在线观看| 青草视频久久| 高清免费毛片| 午夜视频日本| 福利一区三区| 国产麻豆精品在线观看| 国产精品免费露脸视频| 永久免费精品视频| AV无码一区二区三区四区| 中文字幕久久精品波多野结|