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

虛擬同步控制的港口岸電電源阻抗建模及穩定性分析

2022-04-29 00:00:00陳輝潘雄曾燦林宋先勇蒲貞洪郭健陳燕東
湖南大學學報·自然科學版 2022年4期

摘要:為增強系統慣量和阻尼,虛擬同步控制被廣泛應用于港口岸電電源中,但虛擬同步控制的港口岸電電源與船舶PWM整流器負荷之間可能存在交互失穩問題.因此,本文首先根據其多時間尺度控制特性,提出了虛擬同步控制的港口岸電電源的分頻段dq阻抗模型.其次,基于所建dq阻抗和廣義奈奎斯特穩定判據的穩定性分析表明,港口岸電電源的交流電壓環與船舶PWM整流器負荷的直流電壓環之間存在控制交互作用,進而會誘發系統振蕩.增加港口岸電電源的交流電壓比例和諧振系數,或減小船舶PWM整流器負荷的直流側電壓比例系數可增強港口岸電供電系統的穩定性.最后基于硬件在環實驗平臺,驗證了阻抗模型和穩定性分析結果的有效性.

關鍵詞:港口岸電電源;虛擬同步控制;阻抗建模;系統穩定性

中圖分類號:TM910文獻標志碼:A

Impedance Modeling and Stability Analysis of Port Shore Power Supply Based on Virtual Synchronous Control

CHEN Hui PAN Xiongwen ZENG Canlin SONG Xianyong PU Zhenhong GUO Jian CHEN Yandong2

(1.Yueyang Power Supply Company,State Grid,Yueyang 414000,China;

2. National Electric Energy Conversion and Control Engineering Technology Research Center,Hunan University,Changsha 410000,China)

Abstract:Virtual synchronous control is widely studied in port shore power sources to enhance system inertia and damping,but there may be control interaction between the port shore power source and the ship PWM rectifier load. Thus,firstly,according to the multi-time scale control characteristics,a frequency-division dq-frame impedance model of the port shore power source with virtual synchronous control is proposed. Secondly,the stability analyses based on the established dq-frame impedance and generalized Nyquist stability criterion show that there is a control interaction between the AC voltage loop of port shore power source and the DC voltage loop of ship PWM rectifier load,inducing the system oscillation. Increasing the AC voltage proportional and resonance gain of port shore power source or reducing the DC voltage proportional gain of ship PWM rectifier load can enhance the stability of the portshore power supply system. Finally,the effectiveness of the impedance model and stability analysis results are verified by the experimental results based on the hardware in the loop experimental platform.

Key words :port shore power supply;virtual synchronous control;impedance modeling;stability analysis

水路運輸作為傳統的交通運輸方式,對生態環境的影響極大.國際海事組織的數據顯示,全世界以柴油為動力的各類艦船每年向大氣排放1 000萬噸氮氧化物和850萬噸硫氧化物.由于常年采用柴油發電機作為浮吊的動力,油污、油煙、噪音給環境和大氣造成嚴重影響.并且船舶在待閘錨泊期間,也需要通過輔機發電滿足船舶值班、生活、照明設備等用電需求,持續排放出有毒有害物質,如硫化物、碳氧化物、PM2.5,這些嚴重影響了內陸湖泊的生態環境和浮吊產業的經濟發展.我國內陸湖泊急需從現在低效、粗放、污染的柴油機供電體系逐步轉變成潔凈、高效、節約、多元、安全的現代化岸電供電體系[1-2].

船舶岸電供電系統一般由10 kV架空線路作為輸入電壓,經過變壓器和變頻器等裝置輸出400 V/ 50 Hz、440 V/60 Hz和6.6 kV/60 Hz三種電壓等級的船舶供電電壓.如圖1所示,本文主要針440 V/60 Hz 的情況進行研究.為了增強慣性和阻尼,虛擬同步控制已逐漸應用于港口岸電供電電源[3],以模擬同步機的外部特性,提供慣量和阻尼[4-7]。然而,隨著大量電力電子化船舶負載接入岸電系統,恒功率負荷的負阻尼會削弱岸電電源的穩定裕度[8],并可能導致系統振蕩.因此,虛擬同步控制的港口岸電電源和船舶PWM整流器負荷之間的控制相互作用有待進一步研究.

基于阻抗的穩定性分析是處理此類問題的有效方法,其基本原理是將廣義奈奎斯特穩定判據(GNC)應用于系統的源荷阻抗比[9-15].根據坐標系的不同,阻抗形式可分為如下幾類:αβ軸阻抗,極坐標阻抗,序阻抗[15-16]和dq阻抗.其中,dq阻抗和序阻抗的研究較為廣泛,考慮頻率耦合后,序阻抗和dq阻抗都為二維矩陣,且基于序阻抗和dq阻抗的穩定性分析方法本質上是相同的[17].然而,在dq坐標系下,三相平衡的變流器是線性時不變系統.因此可對各環節直接線性化,對所得的線性時不變模型進行拉普拉斯變換即可推導dq阻抗.因此,dq阻抗建模更簡單,阻抗模型表達也更簡潔.

目前,考慮鎖相環(PLL)、直流電壓等環節[18],已有大量文獻建立了變流器的dq阻抗模型.并且,文獻[19]分析了逆變電源、PWM整流器負荷和電網之間的控制作用,并從PWM整流器的角度提出了抑制振蕩的方法.文獻[20]利用測量的dq阻抗探討了逆變電源和PWM整流器負荷之間的控制相互作用.此外,文獻[21]建立了含dq軸控制器的虛擬同步機的dq阻抗模型.一方面,上述文獻都忽略了控制延時和采樣濾波器在dq坐標系和靜止坐標系下的模型差異性.另一方面,由于港口供電電源存在多時間尺度控制特性,在各個頻段的阻抗特性呈現顯著的差異性.因此,在保證模型精度的前提下,可對基于虛擬同步控制的港口岸電電源進行分頻段阻抗建模,進一步探索虛擬同步控制的港口岸電供電系統的小信號穩定性.

為解決上述問題,本文基于虛擬同步控制的港口岸電電源的多時間尺度控制特性,建立了其分頻段dq阻抗模型.此外,通過阻抗測量[22-23]直接驗證了所提阻抗模型的準確性和有效性.最后,本文基于所建dq阻抗模型和廣義奈奎斯特穩定判據(GNC)分析了虛擬同步機控制的港口岸電系統的小信號穩定性.

1虛擬同步控制的港口岸電系統

1.1系統描述

圖2展示了基于虛擬同步控制的港口岸電供電系統的簡化電路,其中,Lf、Rf和Cf分別為港口岸電電源逆變側的交流側濾波電感、電阻和電容;ea、eb和ec為岸電供電電源的輸出電壓;ia、ib和ic是岸電供電電源的電感電流;Uab和Ubc是PCC電壓.

如圖2所示,Udc1是船舶PWM整流器的直流側電壓;ial、ibl和icl是PWM整流器的電感電流;Lfl、Rfl和Cd為PWM整流器濾波電感、寄生電阻和直流側電容.本文研究岸電電源與船舶負荷的控制交互作用,因此對實際系統進行簡化.考慮到港口岸電電源的逆變器直流側電壓由強電網經二極管整流及電容濾波所得,視為常數;船舶負荷類型多樣,含大量PWM 整流器負荷(電動機負荷),在后文的分析中,考慮最惡劣的情況,即船舶負荷選取為PWM整流器,忽略其他類型負荷(Rload,Lload).港口岸電電源逆變側的主電路可表示如下:

式中,esd,esq是港口岸電電源的端口電壓;iLd,iLq為輸出電流的dq分量;ud,uq是PCC電壓的dq分量.

1.2港口岸電電源的虛擬同步控制

港口岸電電源的逆變器采用虛擬同步控制,其有功環模擬了同步發電機的慣性和一次調頻特性.

式中:J為虛擬轉動慣量;ω和ωn分別為虛擬同步控制的輸出角頻率和額定角頻率;Te是電磁轉矩的給定值;Dp為有功阻尼系數;θ是VSG的相位.

無功環模擬同步電機的一次調壓特性.

Es=(Q*+Dq(U*-Un)-Q)/Kqs(3)

式中:Es為內電勢有效值;Q*為瞬時輸出無功功率Q的指令值;Dq是無功阻尼系數;Kq是無功環慣性系數;Un是電壓幅值U的額定值.

其中,瞬時有功和無功功率的計算方式如下:

由無功環輸出的電壓幅值和有功環輸出的相角可得αβ軸電壓的給定值.

電壓環采用準比例諧振控制器.

式中:kpv,krv和ωr分別為比例系數、諧振系數和低通截止頻率.

電流環采用比例控制.

式中:kpi是電流控制器的比例系數.

2港口岸電電源的阻抗建模

2.1控制延時與采樣濾波器的dq建模

控制延時在αβ軸下的小信號模型表示如下

式中:dαc(s)和dβc(s)控制延時前αβ軸上的占空比;dαs(s)和dβs(s)控制延時后αβ軸上的占空比;Ts=1.5/fs,fs指開關頻率.

根據靜止坐標系到dq坐標系下傳遞函數的轉換關系,可推導港口岸電電源控制延時的dq模型如式(9)所示.

電壓或電流采樣低通濾波器的表達式如下:

式中:x表示電流(i)或電壓(v);Tx=1/ωxc,ωxc是電壓或電流信號低通濾波器的截止頻率.

近似地,電壓或電流的低通濾波器的dq模型可推導如式(11)所示.

2.2虛擬同步控制的港口岸電電源的分頻段dq阻抗建模

岸電電源的虛擬同步控制具有多時間尺度特性,不同頻段下對應不同的環節動態.在1~20 Hz的低頻段主要的動態環節為功率控制器與電壓控制器;20~100 Hz的中低頻段主要動態環節為電壓環、功率環和電流環;100 Hz至幾百赫茲的中頻段主要動態環節為電壓環和電流環;幾百赫茲至2 kHz的高頻段主要的動態環節為電流環、采樣濾波器、控制延時和電壓環.因此,在dq阻抗建模過程中,可在不同頻段內考慮主要動態環節進行建模.

圖3展示了基于虛擬同步控制的港口岸電電源的寬頻帶小信號模型,港口岸電電源的輸出電壓與PCC點電壓之間存在穩態相位差δ0,因此港口岸電電源的輸出電壓轉換為PCC點電壓時存在如下轉換關系:

式中:δ00=P*/(3U0E0ωnLf);U0是PCC點的額定電壓;E0是港口岸電電源輸出端口的額定電壓.

當港口岸電電源的端口注入dq軸的電壓擾動時,其輸出電壓可表示為式(13).

通過抵消式(13)中的穩態分量,并消除二次擾動分量可得dq軸向的電壓擾動變化量.

由式(14)可得有功無功功率的小信號模型如式(15)所示.

式中:Ud0和Uq0為PCC點電壓的dq分量;Id0和Iq0是輸出電流的dq分量.

根據式(2)和式(3),有功和無功功率控制器的小信號模型表示如下:

根據控制器從αβ軸到dq軸的轉變方式[24],電流比例控制器及電壓準比例諧振控制器的dq模型可推導如式(17)和式(18)所示.

在低頻區域,控制延時和電壓電流信號的低通濾波器等環節可忽略.此時,電流環可等效為“1”.港口岸電電源的dq模型可簡化為圖4所示.

根據圖 在低頻區域,港口岸電電源的dq阻抗模型推導如下:

式中:Zc和Zl分別是濾波電容和電感的dq阻抗.

在中低頻區域,控制延時和采樣濾波器的dq模型可以忽略.此時,港口岸電電源的dq阻抗模型如式(20)所示.

在中頻區域,考慮電壓電流控制環,忽略功率環、采樣濾波器和控制延時,港口岸電電源的δθ阻抗可推導如下:

在高頻區域,虛擬同步控制的港口岸電電源的功率控制環可以忽略,此時,其dq阻抗模型可推導如下:

為驗證本文所建分頻段dq阻抗模型的準確性,基于圖5所示的硬件在環實驗平臺搭建了虛擬同步控制的港口岸電電源給PWM整流器負荷供電的系統,其中虛擬同步控制和PWM整流器的控制分別在兩個TI DSP28335控制器中執行,控制器通過輸入輸出接口與RT_LAB連接,系統主電路及擾動電流源由RT_LAB仿真實現.表1和表2分別展示了港口岸電電源和PWM整流器的控制及系統參數,采用并聯注入理想電流源擾動的方式測量虛擬同步控制的港口岸電電源的dq阻抗.圖6展示了港口岸電電源的分頻段阻抗模型和基于RT_LAB實驗平臺的測量結果.由圖6可知,盡管低頻段模型與測量值有細微差異,港口岸電電源的分頻段阻抗模型與測量值基本吻合,這證實了所建分頻段dq阻抗模型的準確性.

3港口岸電系統的穩定性分析

3.1基于分頻段dq阻抗模型的系統穩定性分析

受篇幅限制,本文不展示PWM整流器的阻抗模型Zvsr,讀者可參考文獻[24].由圖7可知,PWM整流器的dq阻抗模型與仿真測量結果是吻合的.

港口岸電電源和PWM整流器負荷的dq軸耦合阻抗不可忽略.因此,GNC需應用于其阻抗比:

Q(s)=Zvsi/Zvsr(23)

當且僅當Q(s)的特征根逆時針包圍(- 0)的次數等于Zvsi和Zvsr-1的右半平面極點個數時,系統才穩定.考慮到Zvsi和Zvsr-1不存在右半平面極點,當Q(s)特征根的奈奎斯特曲線不逆時針包含(- j0)時,系統是穩定的.而Q(s)的特征根的解如下:

det(λI-Q(s))=0(24)

3.2港口岸電電源與船舶PWM整流器負荷的控制相互作用分析

將式(19)~式(23)代入式(24),圖8展示了參數變化時,Zvsi/Zvsr的特征根.同時,為了驗證系統的穩定性分析,在RT_LAB中進行了實驗驗證.參數基于表1~表 圖9顯示了與圖8對應的系統交流接口電壓的仿真波形.圖8(a)顯示了僅改變港口岸電電源的交流電壓比例系數時,Zvsi/Zvsr的特征根的奈奎斯特曲線.隨著港口岸電電源的kp_vsg的降低,λ1和λ2逐漸包圍(- j0),這意味著減小港口岸電電源的電壓比例系數,系統變得不穩定.圖9(a)、(c)證明圖8(a)分析的正確性,kp_vsg等于0.0 系統振蕩;當僅增大kp_vsg為0.03時,系統由不穩定變為穩定.圖8(b)顯示了僅改變港口岸電電源的電壓諧振系數時,Zvsi/Zvsr的特征根的奈奎斯特曲線.隨著kr_vsg的降低,λ1和λ2逐漸包圍(- j0),這意味著減小港口岸電電源的電壓諧振系數,系統同樣變得不穩定.圖9(a)、(d)證實了圖8(b)分析的正確性.圖8(c)顯示僅改變船舶PWM整流器負荷的直流側電壓比例系數時,Zvsi/Zvsr的特征根的奈奎斯特曲線.隨著kr_vst的增大,λ1和λ2逐漸包圍(- j0),這意味著增大船舶PWM整流器負荷的直流側電壓比例系數,系統變得不穩定.圖9(a)、(b)證實了圖8(c)分析的正確性.

4結論

本文采用dq阻抗分析法對虛擬同步控制的港口岸電系統進行了穩定性分析,得出以下結論:

1)根據虛擬同步控制的港口岸電電源的多時間尺度控制特性,本文提出了虛擬同步控制的港口岸電電源的分頻段dq阻抗模型,并通過阻抗測量驗證了模型的準確性.

2)船舶PWM整流器負荷呈現負阻尼特性,虛擬同步控制的港口岸電電源的交流電壓環與船舶PWM整流器負荷的直流電壓環之間的相互作用可能導致系統不穩定.增加港口岸電電源電壓比例和諧振系數,或減小船舶PWM整流器負荷的直流側電壓比例系數可增強系統的穩定性.

參考文獻

[1]劉杜,孫佳星,喬康恒,等.港口船舶岸電電源系統研究綜述[J].船電技術,202 41(6):29-34.

LIU D,SUN J X,QIAO K H,et al. Review of development of the shore power system for port ship[J].Marine Electric amp; Electronic Engineering,202 41(6):29-34.(In Chinese)

[2]陳澤西,孫玉樹,張妍,等.考慮風電波動率的儲能系統優化配置策略[J].湖南大學學報(自然科學版),2020,47(8):60-68.

CHEN Z X,SUN Y S,ZHANG Y,et al. Optimal configuration strategy of energy storage system considering wind power fluctuation rate [J]. Journal of Hunan University (Natural Sciences Edition),2020,47(8):60-68.(In Chinese)

[3]邱長青,朱強,羅文斌.基于虛擬同步發電機的岸電并網控制策略研究[J].船電技術,2020,40(2):22-27.

QIU C Q,ZHU Q,LUO W B.A control strategy study on shore and ship power grid-connection based on virtual synchronous generator [J]. Marine Electric amp; Electronic Engineering,2020,40(2):22- 27.(In Chinese)

[4] ZHONG Q C,WEISS G. Synchronverters :inverters that mimic synchronous generators [J].IEEE Transactions on Industrial Electronics ,201 58(4):1259-1267.

[5] TORRES L M A,LOPES LA C,MORAN TLA,et al. Self-tuning virtual synchronous machine:a control strategy for energy storage systems to support dynamic frequency control [J ]. IEEE Transactions on Energy Conversion,201 29(4):833-840.

[6] ZHOU L M,SHUAI Z K,CHEN Y D,et al. Impedance-based harmonic current suppression method for VSG connected to distorted grid[J]. IEEE Transactions on Industrial Electronics,2020,67 (7):5490-5502.

[7] WU H,RUAN X B,YANG D S,et al. Small-signal modeling and parameters design for virtual synchronous generators[J]. IEEE Transactions on Industrial Electronics,2016,63(7):4292-4303.

[8]熊義勇,趙鏡紅,劉洋,等.典型中壓直流岸電系統小干擾穩定性分析[J].海軍工程大學學報,202 33(3):75-80.

XIONG Y Y,ZHAO J H,LIU Y,et al. Small signal stability analysis of typical medium voltage DC shore power system[J] .Journal of Naval University of Engineering,202 33(3):75-80.(In Chinese)

[9] YUE X L,WANG X F,BLAABJERG F. Review of small-signal modeling methods including frequency-coupling dynamics of power converters [J] . IEEE Transactions on Power Electronics,2019,34(4):3313-3328.

[10] SUN J. Small-signal methods for AC distributed power systems—a review [J] . IEEE Transactions on Power Electronics,2009,24 (11):2545-2554.

[11] BELKHAYAT M. Stability criteria for AC power systems with regulated loads [D] . West Lafayette,IN,USA:Purdue University,1997.

[12] SUN J. Impedance-based stability criterion for grid-connected inverters[J] . IEEE Transactions on Power Electronics,201 26(11):3075-3078.

[13] BURGOS R,BOROYEVICH D ,WANG F,et al. On the Ac stability of high power factor three-phase rectifiers [C]//2010 IEEE Energy Conversion Congress and Exposition. Atlanta,GA,USA:IEEE,2010:2047-2054 .

[ 14] WANG X F,BLAABJERG F,WU W M. Modeling and analysis of harmonic stability in an AC power-electronics-based power system[J] . IEEE Transactions on Power Electronics,201 29(12):6421-6432.

[15] CESPEDES M,SUN J. Impedance modeling and analysis of grid- connected voltage-source converters[J] . IEEE Transactions on Power Electronics,201 29(3):1254-1261.

[16]伍文華,陳燕東,周樂明,等.虛擬同步發電機接入弱電網的序阻抗建模與穩定性分析[J].中國電機工程學報,2019,39(6):1560-1571.

WU W H,CHEN Y D,ZHOU L M,et al. Sequence impedance modeling and stability analysis for virtual synchronous generator connected to the weak grid [J].Proceedings of the CSEE ,2019,39 (6):1560-1571(. In Chinese)

[17] RYGG A,MOLINAS M,ZHANG C,et al. On the equivalence and impact on stability of impedance modeling of power electronic converters in different domains[J]. IEEE Journal of Emerging and Selected Topics in Power Electronics,2017,5(4):1444-1454.

[18] RADWAN A AA,MOHAMED Y A R I. Analysis and activeimpedance-based stabilization of voltage-source-rectifier loads in grid-connected and isolated microgrid applications [J] . IEEE Transactions on Sustainable Energy,201 4(3):563-576.

[19] LU D P,WANG X F,BLAABJERG F. Impedance-based analysis of DC-link voltage dynamics in voltage-source converters[J]. IEEE Transactions on Power Electronics,2019,34(4):3973- 3985.

[20] WEN B,BOROYEVICH D,BURGOS R,et al. Small-signal stability analysis of three-phase AC systems in the presence of constant power loads based on measured d-q frame impedances[J] . IEEE Transactions on Power Electronics,2015,30(10):5952-5963.

[21] WANG S K,LIU Z,LIU J J. Modeling of D-Q small-signal impe dance of virtual synchronous generator [C]〃2018 IEEE International Power Electronics and Application Conference and Exposition. Shenzhen,China:IEEE,2018:1-6.

[22] GONG H,YANG D S,WANG X F. Impact analysis and mitigation of synchronization dynamics for DQ impedance measurement[J] . IEEE Transactions on Power Electronics,2019,34(9):8797- 8807.

[23]黃瑞,余敏琪,劉奕玹,等.基于Karrenbauer相模變換的低壓供電線路阻抗量測方法[J].湖南大學學報(自然科學版),202 48(8):90-95.

HUANG R,YU M Q,LIU Y X,et al. Impedance measurement method of low-voltage power supply line based on karrenbauer phase-mode transformation[J] .Journal of Hunan University(Natural Sciences),202 48(8):90-95.(In Chinese)

[24]郭健,陳燕東,王翔宇,等.負荷虛擬同步機的寬頻帶dq阻抗建模及弱電網下與傳統PWM整流器的穩定性對比分析[J].中國電機工程學報,2020,40(15):4758-4770.

GUO J,CHEN Y D,WANG X Y,et al. Wideband dq-frame impedance modeling of load-side virtual synchronous machine and its stability analysis in comparison with conventional PWM rectifier in weak grid[J]. Proceedings of the CSEE,2020,40(15):4758-4770.(In Chinese)

主站蜘蛛池模板: 国内精品视频| 国产经典免费播放视频| 在线观看无码av五月花| 精品午夜国产福利观看| 玖玖免费视频在线观看| 免费一级全黄少妇性色生活片| 欧洲日本亚洲中文字幕| 亚洲成人黄色网址| 国产日韩精品欧美一区灰| 国产精品久久久久鬼色| 国产第一色| 香蕉视频在线观看www| 国产精品主播| av在线人妻熟妇| 91精品国产情侣高潮露脸| 曰AV在线无码| 亚洲午夜天堂| 无码人中文字幕| 国产人成在线观看| 91黄色在线观看| 亚洲综合色区在线播放2019| 狂欢视频在线观看不卡| 欧美精品啪啪一区二区三区| 无码国产偷倩在线播放老年人| 免费一级毛片不卡在线播放| 色亚洲成人| 国产精品青青| 首页亚洲国产丝袜长腿综合| 特级精品毛片免费观看| 日韩欧美一区在线观看| 国产福利微拍精品一区二区| 国产成人1024精品| 亚洲第一中文字幕| 国产精品亚洲五月天高清| 永久毛片在线播| 国产在线观看一区二区三区| 久草视频中文| 99热亚洲精品6码| 亚洲成A人V欧美综合天堂| 99精品高清在线播放| 国产成人精品一区二区三在线观看| 任我操在线视频| 国产微拍一区二区三区四区| 中文字幕有乳无码| AV不卡在线永久免费观看| 欧美日韩在线国产| 日本三级黄在线观看| 日韩在线成年视频人网站观看| 日韩国产另类| 萌白酱国产一区二区| 网久久综合| 日韩AV无码免费一二三区| 色偷偷一区二区三区| 亚洲日韩高清在线亚洲专区| 欧美亚洲日韩不卡在线在线观看| 欧美日韩国产系列在线观看| 亚洲一级毛片免费观看| 精品人妻无码中字系列| 欧美国产在线一区| 久久免费精品琪琪| 五月婷婷丁香综合| 国产成人AV男人的天堂| 亚洲综合九九| 激情无码字幕综合| 久久精品丝袜高跟鞋| 免费国产无遮挡又黄又爽| 国产高清在线观看| 久久频这里精品99香蕉久网址| 伊人天堂网| 在线观看亚洲国产| 一区二区日韩国产精久久| 欧美一区日韩一区中文字幕页| 婷婷六月综合网| 欧美一区日韩一区中文字幕页| 国产一区二区免费播放| 永久免费av网站可以直接看的| 国产资源站| 精品国产网站| 色综合国产| 婷婷综合在线观看丁香| 伊人色综合久久天天| 欧美第一页在线|