• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于序阻抗的虛擬同步機(jī)同步頻率諧振現(xiàn)象

    2022-05-20 07:44:32于彥雪關(guān)萬琳陳曉光李浩昱
    電工技術(shù)學(xué)報(bào) 2022年10期
    關(guān)鍵詞:環(huán)控制諧振阻尼

    于彥雪 關(guān)萬琳 陳曉光 李浩昱

    基于序阻抗的虛擬同步機(jī)同步頻率諧振現(xiàn)象

    于彥雪1關(guān)萬琳2陳曉光2李浩昱1

    (1. 哈爾濱工業(yè)大學(xué)電氣工程及自動(dòng)化學(xué)院 哈爾濱 150001 2. 國(guó)網(wǎng)黑龍江省電力有限公司電力科學(xué)研究院 哈爾濱 150030)

    針對(duì)虛擬同步機(jī)(VSG)的同步頻率諧振(SFR)現(xiàn)象,鑒于多輸入多輸出功率耦合模型的不足,通過建立可描述VSG寬頻域動(dòng)態(tài)特性的單輸入單輸出序阻抗模型,提出一種基于序?qū)Ъ{和序阻抗的SFR分析方法,并對(duì)單環(huán)和多環(huán)控制VSG的SFR現(xiàn)象進(jìn)行對(duì)比研究。基于序?qū)Ъ{的極點(diǎn)分布,有效研究SFR與VSG功率環(huán)控制參數(shù)、虛擬阻抗、內(nèi)環(huán)控制及電網(wǎng)阻抗之間的關(guān)系;基于序阻抗的頻率特性曲線,對(duì)比內(nèi)環(huán)控制引入前后SFR的阻尼特性,從“負(fù)阻尼”的角度揭示單環(huán)和多環(huán)控制VSG的SFR諧振機(jī)理;搭建基于RT-Box的實(shí)驗(yàn)平臺(tái),驗(yàn)證VSG在不同控制方案及不同電網(wǎng)阻抗下的SFR現(xiàn)象,指明VSG控制的優(yōu)化設(shè)計(jì)方向。

    虛擬同步機(jī) 同步頻率諧振 序阻抗 負(fù)阻尼 內(nèi)環(huán)控制

    0 引言

    作為一種模擬同步發(fā)電機(jī)外特性的并網(wǎng)逆變器控制技術(shù),虛擬同步機(jī)(Virtual Synchronous Gener- ator, VSG)不僅可為系統(tǒng)提供電壓和頻率支撐[1-2],還可改善高滲透率新能源發(fā)電并網(wǎng)系統(tǒng)的穩(wěn)定性[3]。此外,虛擬慣性和阻尼的加入,給系統(tǒng)動(dòng)態(tài)特性和穩(wěn)態(tài)性能的提高增加了更多的自由度[4]。目前有關(guān)VSG功率振蕩現(xiàn)象的研究大多忽略了電感磁鏈的動(dòng)態(tài)過程[5],其慣性和阻尼參數(shù)的選取依靠功率環(huán)簡(jiǎn)化的二階模型。這種忽略多建立在VSG的整體響應(yīng)帶寬遠(yuǎn)小于系統(tǒng)同步頻率的前提下[6],此時(shí)虛擬慣性系數(shù)較大,電感磁鏈的動(dòng)態(tài)變化量較小。然而,在新能源滲透率較大的新型電網(wǎng)中,要求VSG的響應(yīng)速度能夠滿足源載快速變化的要求,若VSG功率控制參數(shù)設(shè)計(jì)不當(dāng),則會(huì)激發(fā)其電磁動(dòng)態(tài)過程,引起同步頻率諧振(Synchronous Frequency Resonance, SFR)問題。

    根據(jù)內(nèi)部電壓電流雙閉環(huán)控制存在與否,可將VSG分為只含功率環(huán)的單環(huán)控制VSG以及含功率環(huán)和內(nèi)部雙閉環(huán)的多環(huán)控制VSG兩種。目前有關(guān)SFR問題的研究?jī)H局限于單環(huán)控制VSG,研究?jī)?nèi)容涉及SFR建模、諧振機(jī)理分析與諧振抑制[7-10]。其中,文獻(xiàn)[7]通過分析電磁動(dòng)態(tài),建立了考慮功率耦合時(shí)VSG功率閉環(huán)控制的傳遞函數(shù),揭示了SFR的諧振機(jī)理及其與功率控制參數(shù)之間的關(guān)系。文獻(xiàn)[8]通過VSG功率閉環(huán)控制傳遞函數(shù),指出線路阻抗在dq坐標(biāo)系下的頻移是導(dǎo)致功率產(chǎn)生SFR的原因,并提出了一種相位-幅值交互調(diào)節(jié)方法來抑制SFR。文獻(xiàn)[9]針對(duì)阻感比較小時(shí),VSG控制在同步頻率處相位裕度不足的特點(diǎn),提出了一種相位補(bǔ)償加前饋補(bǔ)償?shù)目刂品椒▉砀纳芐FR。文獻(xiàn)[10]針對(duì)傳統(tǒng)虛擬電阻所引起的功率耦合加強(qiáng)的副作用,提出了一種暫態(tài)虛擬電阻方法來抑制SFR。對(duì)于多環(huán)控制VSG,文獻(xiàn)[6]指出功率環(huán)、虛擬阻抗、內(nèi)環(huán)控制和電磁動(dòng)態(tài)之間均存在相互作用,鑒于多環(huán)控制VSG技術(shù)的普遍應(yīng)用,尚需對(duì)多環(huán)控制VSG的SFR現(xiàn)象進(jìn)行研究,揭示虛擬阻抗和內(nèi)環(huán)控制對(duì)SFR的影響。

    此外,已有研究表明,功率耦合會(huì)惡化SFR[7-8],因此,SFR的建模分析需考慮VSG的功率耦合現(xiàn)象。然而,當(dāng)考慮功率耦合時(shí),現(xiàn)有的SFR分析模型均為基于功率閉環(huán)傳遞函數(shù)的多輸入多輸出或等效的單輸入單輸出模型,模型較為復(fù)雜,且難以應(yīng)用于VSG的寬頻域穩(wěn)定性分析。對(duì)VSG控制并網(wǎng)逆變器來說,其穩(wěn)定性問題可能產(chǎn)生于多個(gè)頻率段[7, 11-13],因此,有必要將對(duì)SFR的研究統(tǒng)一于VSG其他穩(wěn)定性問題的研究模型中,建立可描述VSG控制并網(wǎng)逆變器寬頻域動(dòng)態(tài)特性的穩(wěn)定性分析模型。

    序阻抗作為一種描述系統(tǒng)端口特性的模型,被廣泛應(yīng)用在逆變器并網(wǎng)系統(tǒng)的小信號(hào)交互穩(wěn)定性研究中[14]。在VSG控制逆變器并網(wǎng)系統(tǒng)的穩(wěn)定性分析中,序阻抗也得到了有效應(yīng)用[15]。但現(xiàn)有的VSG序阻抗模型應(yīng)用局限于分析弱電網(wǎng)下VSG控制逆變器并網(wǎng)系統(tǒng)的交互穩(wěn)定性。鑒于SFR與線路阻尼密切相關(guān),而序阻抗能反映VSG的寬頻域阻尼特性。因此,本文旨在將SFR的建模分析統(tǒng)一于序阻抗模型中,通過序阻抗模型揭示SFR的諧振機(jī)理。

    具體來說,本文以基于VSG的電壓控制型逆變器(VSG-based Voltage Controlled Inverter, VSG- VCI)為研究對(duì)象,綜合考慮有功-頻率、無功-電壓控制、虛擬阻抗、內(nèi)環(huán)控制和電磁動(dòng)態(tài),采用諧波線性化方法,詳細(xì)建立單環(huán)和多環(huán)控制型VSG-VCI的序阻抗模型;提出通過序?qū)Ъ{的極點(diǎn)分布對(duì)SFR現(xiàn)象進(jìn)行分析的方法,并站在序阻抗“負(fù)阻尼”的角度,闡明SFR的產(chǎn)生及抑制原理;對(duì)比單環(huán)和多環(huán)控制型VSG-VCI功率諧振的異同,全面研究SFR與VSG功率控制參數(shù)、虛擬阻抗及電網(wǎng)阻抗之間的關(guān)系,并進(jìn)行實(shí)驗(yàn)驗(yàn)證。

    1 VSG控制系統(tǒng)及其SFR現(xiàn)象介紹

    圖1所示為本文所研究的VSG-VCI并網(wǎng)系統(tǒng)結(jié)構(gòu),包含功率環(huán)、虛擬阻抗和內(nèi)部電壓電流雙閉環(huán),是一個(gè)多環(huán)控制系統(tǒng)。主電路中,dc為恒定的直流母線電壓;1、1、2組成LCL濾波器,1為其阻尼電阻;vi分別為系統(tǒng)采樣電壓和電流,v()和i()分別為電壓和電流采樣函數(shù);uu分別為逆變橋輸出電壓和電容支路電壓;g和g分別為電網(wǎng)電壓及其線路阻抗,g=g+jLg。VSG功率環(huán)中,N、N為有功、無功功率給定;e、e為經(jīng)式(1)計(jì)算得到的有功、無功功率,p()為功率濾波器;N、N為電壓、角頻率給定值;為虛擬慣量;p為有功阻尼系數(shù);q為無功下垂系數(shù);v為有功環(huán)輸出的相位;v為無功環(huán)輸出的電壓幅值,它們通過式(2)為內(nèi)環(huán)提供參考電壓v。v和v分別為虛擬電感和虛擬電阻。VSG內(nèi)環(huán)控制中,rv()和ri()分別為ab坐標(biāo)系下的電壓和電流控制器。多環(huán)控制型VSG-VCI的PWM信號(hào)取自內(nèi)環(huán)控制的輸出電壓[iaib]T;而單環(huán)控制型VSG-VCI不含內(nèi)環(huán)控制,其PWM信號(hào)取自內(nèi)環(huán)輸入電壓[rarb]T。為了全面認(rèn)識(shí)VSG的SFR現(xiàn)象,本文將分別對(duì)單環(huán)和多環(huán)控制型VSG-VCI進(jìn)行SFR分析。此外,=a, b, c,代表a、b、c三相。

    (2)

    圖2 VSG-VCI功率閉環(huán)控制框圖

    式中,和+v為VSG輸出電動(dòng)勢(shì)v與電網(wǎng)電壓g之間的等效線路電阻和電感,對(duì)單環(huán)控制型VSG- VCI來說,=v+g,=1+2+g;對(duì)多環(huán)控制型VSG-VCI來說,=v+g,=g。因此,單環(huán)控制型VSG-VCI的電磁動(dòng)態(tài)不僅與電網(wǎng)阻抗g有關(guān),還關(guān)系到濾波總電感1+2;而多環(huán)控制型VSG- VCI的電磁動(dòng)態(tài)只與電網(wǎng)阻抗g有關(guān)。式(3)中的小信號(hào)模型具有共同的共軛極點(diǎn),即

    根據(jù)圖2所示的功率閉環(huán)控制框圖,功率環(huán)、內(nèi)環(huán)控制、電磁動(dòng)態(tài)之間存在相互作用,影響式(4)所描述的功率諧振狀態(tài)。現(xiàn)有的基于功率閉環(huán)傳遞函數(shù)的SFR分析需分別推導(dǎo)圖2中各環(huán)節(jié)的表達(dá)式,較為復(fù)雜。此外,功率耦合作用使得SFR分析模型需以多輸入多輸出或等效為單輸入單輸出的形式表達(dá),這為SFR的分析帶來不便[7]。因此,本文提出通過單輸入單輸出序阻抗模型,從VSG-VCI并網(wǎng)端口特性的角度對(duì)SFR問題進(jìn)行研究。

    2 VSG-VCI的序阻抗建模及驗(yàn)證

    2.1 序阻抗建模

    首先,由圖1可得逆變器的主電路平均模型為

    其中

    根據(jù)諧波線性化建模過程[14],在v中注入小信號(hào)電壓擾動(dòng),則A相輸出電壓和電流可表示為

    式中,1、1分別為基波電壓的幅值、頻率;p、p、vp分別為正序電壓擾動(dòng)的幅值、頻率、相位;n、n、vn分別為負(fù)序電壓擾動(dòng)的幅值、頻率、相位;1、i1分別為逆變器輸出基波電流的幅值、相位;p、ip分別為輸出正序電流諧波的幅值、相位;n、in分別為輸出負(fù)序電流諧波的幅值、相位。

    式(6)和式(7)對(duì)應(yīng)的頻域信號(hào)經(jīng)Clarke變換后可表示為

    式中,“*”為其對(duì)應(yīng)復(fù)數(shù)的共軛。

    將式(9)和式(10)代入圖1中的虛擬阻抗控制,可得內(nèi)環(huán)輸入a軸電壓ra的頻域表達(dá)式為

    式中,v=v±jNv。

    將式(8)、式(9)和式(11)代入內(nèi)環(huán)控制,可得內(nèi)環(huán)輸出a軸電壓ia的頻域表達(dá)式為

    對(duì)單環(huán)控制型VSG-VCI來說,頻域下調(diào)制信號(hào)滿足a[]=ra[],將其代入到主電路表達(dá)式(5),計(jì)算可得正序阻抗sp()和負(fù)序阻抗sn()的表達(dá)式如式(13)所示。對(duì)多環(huán)控制型VSG-VCI來說,頻域下滿足a[]=ia[],將其代入到主電路表達(dá)式(5),計(jì)算可得正序阻抗mp()和負(fù)序阻抗mn()的表達(dá)式如式(14)所示。

    2.2 序阻抗模型驗(yàn)證

    VSG-VCI的系統(tǒng)參數(shù)見表1。本文借助Plecs軟件通過阻抗掃描的方式對(duì)所建立的序阻抗模型進(jìn)行分析驗(yàn)證。圖3和圖4分別為單環(huán)和多環(huán)控制型VSG- VCI的序阻抗模型及其仿真測(cè)量結(jié)果,可見本文建立的序阻抗模型與仿真得到的結(jié)果相吻合,證明了該建模的準(zhǔn)確性。

    表1 VSG-VCI的系統(tǒng)參數(shù)

    Tab.1 Parameters of the studied VSG-VCI

    此外,從整體模型來看,單環(huán)控制型VSG-VCI的輸出阻抗在低頻段呈感性,而多環(huán)控制型VSG- VCI的輸出阻抗在低頻段呈阻感特性,表明內(nèi)環(huán)控制的引入改變了VSG-VCI的阻抗特性。

    3 基于序?qū)Ъ{和序阻抗的SFR分析

    3.1 基于序?qū)Ъ{的SFR諧振特性分析

    3.1.1 SFR與序?qū)Ъ{之間的關(guān)系

    根據(jù)序阻抗的建模過程,可以看到序阻抗中包含了VSG-VCI的控制環(huán)路及主電路信息。在并網(wǎng)模式下,注入固定電壓擾動(dòng),引起的電流響應(yīng)實(shí)際也代表了相應(yīng)的功率響應(yīng),因此,電流與電壓之比,即序?qū)Ъ{(序阻抗的倒數(shù))可以代表整個(gè)系統(tǒng)的閉環(huán)行為,其零極點(diǎn)特性與功率閉環(huán)特性密切相關(guān)。

    圖3 單環(huán)控制型VSG-VCI序阻抗模型驗(yàn)證

    圖4 多環(huán)控制型VSG-VCI序阻抗模型驗(yàn)證

    圖5和圖6分別為理想電網(wǎng)下(g=0時(shí)),單環(huán)和多環(huán)控制型VSG-VCI序?qū)Ъ{的零極點(diǎn)。其中,單環(huán)控制型VSG-VCI的v=0mHv=0.3W,多環(huán)控制型VSG-VCI的v=5mHv=0.3W,其他參數(shù)與表1保持一致。如圖5和圖6所示,VSG-VCI輸出序?qū)Ъ{的零極點(diǎn)主要分布在三個(gè)頻率段:與LCL濾波和數(shù)字控制延時(shí)相關(guān)的高頻段(kHz)、與控制環(huán)參數(shù)相關(guān)的中頻段(50~100Hz)和基頻段(50Hz)。

    圖5 單環(huán)控制型VSG-VCI序?qū)Ъ{的零極點(diǎn)

    圖6 多環(huán)控制型VSG-VCI序?qū)Ъ{的零極點(diǎn)

    此外,圖5所示的序?qū)Ъ{中還存在一對(duì)特殊的共軛極點(diǎn),它們的頻率遠(yuǎn)小于基頻。考慮到并網(wǎng)功率與電壓、電流、序阻抗之間存在50Hz頻移,所以,這對(duì)共軛極點(diǎn)可反映VSG并網(wǎng)功率在(10.3+50)Hz處的穩(wěn)定性。然而,由式(4)可知,當(dāng)v=0mH時(shí),由電磁動(dòng)態(tài)引發(fā)的功率諧振頻率為50Hz,其與(10.3+50)Hz存在一定的偏差,這是因?yàn)榍罢呶磳⒖刂骗h(huán)的交互影響考慮在內(nèi)。因此,通過虛線所指出的極點(diǎn)可以分析SFR現(xiàn)象,圖中,該SFR極點(diǎn)實(shí)部小于0,表明0.3W的虛擬電阻可完全抑制SFR;反之,當(dāng)SFR極點(diǎn)移動(dòng)到右半平面時(shí),表明SFR問題產(chǎn)生。

    以上是對(duì)理想電網(wǎng)下單環(huán)和多環(huán)控制型VSG- VCI輸出序?qū)Ъ{的零極點(diǎn)分析,只涵蓋了VSG等效內(nèi)電動(dòng)勢(shì)與并網(wǎng)點(diǎn)之間的等效線路動(dòng)態(tài)。當(dāng)考慮電網(wǎng)阻抗g的動(dòng)態(tài)時(shí),需根據(jù)式(15)所示的等效序阻抗evp()和evn()來進(jìn)行SFR分析。其中,vp()和vn()分別為單環(huán)或多環(huán)控制型VSG-VCI的正序和負(fù)序阻抗。

    由圖5可得,SFR極點(diǎn)在正負(fù)序?qū)Ъ{中互為共軛極點(diǎn)。因此,只需根據(jù)正序(或負(fù)序)導(dǎo)納即可分析SFR與系統(tǒng)運(yùn)行參數(shù)之間的關(guān)系。

    3.1.2 單環(huán)控制型VSG-VCI的SFR諧振特性分析

    單環(huán)控制型VSG-VCI的電磁動(dòng)態(tài)同時(shí)與濾波電路和電網(wǎng)阻抗動(dòng)態(tài)有關(guān),因此,通過正序?qū)Ъ{sp()= 1/sp()和正序等效導(dǎo)納evp()=1/[sp()+gp()]分別對(duì)理想電網(wǎng)和實(shí)際電網(wǎng)下VSG的SFR進(jìn)行分析,研究虛擬阻抗、功率控制和電網(wǎng)阻抗對(duì)單環(huán)控制VSG功率諧振現(xiàn)象的影響。

    圖7所示為理想電網(wǎng)下,當(dāng)=0.03kg·m2、p=5、q=0.001時(shí),sp的SFR極點(diǎn)隨虛擬電阻v(通過v的變化來改變式(4)中的阻感比/)、虛擬電感v變化的分布情況。由圖7a可知,隨v減小,極點(diǎn)向右半平面移動(dòng),說明VSG-VCI并網(wǎng)系統(tǒng)穩(wěn)定性降低,在所選的v取值范圍內(nèi),當(dāng)v<0.1W時(shí),SFR極點(diǎn)位于右半平面,表明系統(tǒng)失穩(wěn),根據(jù)極點(diǎn)頻率進(jìn)一步可知,失穩(wěn)系統(tǒng)的并網(wǎng)功率存在55.7Hz左右的諧振。由圖7b可知,因電磁動(dòng)態(tài)引發(fā)的功率諧振頻率會(huì)隨v增大而增大,這與式(4)所示的功率諧振結(jié)論一致。比如,當(dāng)v=4mH時(shí),由圖可知,功率諧振頻率約為100Hz,而根據(jù)式(4)計(jì)算可得功率諧振頻率為50(1+2+v)/(1+2)≈100Hz,兩者相符,表明基于序?qū)Ъ{模型進(jìn)行SFR現(xiàn)象分析的準(zhǔn)確性。

    圖7 Ysp的SFR極點(diǎn)隨虛擬阻抗Zv變化的分布

    圖7顯示,理想電網(wǎng)下,當(dāng)=0.03kg·m2、p=5、q=0.001且v=0mH、v=0.1W時(shí),單環(huán)控制型VSG- VCI的SFR可被有效抑制。若改變功率環(huán)的控制參數(shù),基于序?qū)Ъ{SFR極點(diǎn)分布的分析結(jié)果表明(見附圖1):即使v=0.1W的虛擬電阻可提供抑制SFR的阻尼,太?。ㄈ纾?.02kg·m2)或q太大(如q>0.001)均會(huì)導(dǎo)致極點(diǎn)移到右半平面,產(chǎn)生SFR問題,證明了功率環(huán)與電磁動(dòng)態(tài)會(huì)發(fā)生相互作用,從而影響SFR阻尼;而有功阻尼系數(shù)p對(duì)SFR的影響較小,當(dāng)p發(fā)生大范圍變化時(shí),VSG-VCI并網(wǎng)系統(tǒng)的穩(wěn)定性保持不變。

    進(jìn)一步地,考慮電網(wǎng)阻抗g的動(dòng)態(tài),當(dāng)= 0.03kg·m2、p=5、q=0.001時(shí),基于序?qū)Ъ{SFR極點(diǎn)分布的分析結(jié)果表明(見附圖2):只要單環(huán)控制型VSG-VCI本身的控制參數(shù)設(shè)計(jì)得當(dāng)(即理想電網(wǎng)下VSG-VCI不存在SFR問題),那么電網(wǎng)電感g(shù)變化不會(huì)影響SFR,即SFR問題與電網(wǎng)強(qiáng)度無關(guān);而電網(wǎng)電阻g可有效阻尼SFR,且其阻尼效果近似于v,因此,當(dāng)并網(wǎng)線路阻尼不足時(shí),可通過虛擬電阻的方式來削弱電磁動(dòng)態(tài)引發(fā)的功率諧振。

    3.1.3 多環(huán)控制型VSG-VCI的SFR諧振特性分析

    多環(huán)控制型VSG-VCI的電磁動(dòng)態(tài)與濾波電路無關(guān)(即與1+2無關(guān)),僅與電網(wǎng)阻抗g有關(guān)。因此,通過正序等效導(dǎo)納evp()=1/[mp()+gp()]對(duì)實(shí)際電網(wǎng)下VSG的SFR進(jìn)行分析,研究電網(wǎng)阻抗、虛擬阻抗和功率控制對(duì)多環(huán)控制VSG功率諧振現(xiàn)象的影響。

    圖8所示為當(dāng)=0.03kg·m2、p=5、q=0.001時(shí),SFR相關(guān)極點(diǎn)隨電網(wǎng)電感g(shù)和電網(wǎng)電阻g變化的分布,此外,圖中虛擬電阻v和虛擬電感v分別為0W和4mH,因v>0mH,由式(4)可知,與電磁動(dòng)態(tài)有關(guān)的極點(diǎn)頻率大于50Hz,為超SFR極點(diǎn)。由圖8a可見,超SFR極點(diǎn)與g之間呈非線性關(guān)系。在圖8所示的g變化范圍內(nèi),雖然系統(tǒng)總的等效電阻=g+v=0W,但只有當(dāng)g取值為3~6mH時(shí),evp()才存在右半平面極點(diǎn),且極點(diǎn)頻率約為50Hz,說明VSG-VCI并網(wǎng)系統(tǒng)存在100Hz左右的功率諧振。當(dāng)g=4mH時(shí),由圖8b可見,將g增加到0.2W,超SFR極點(diǎn)移動(dòng)到左半平面,且g越大,超SFR極點(diǎn)離虛軸越遠(yuǎn),說明增加電網(wǎng)線路的阻性成分有助于抑制該功率諧振。

    圖8 Yevp的超SFR極點(diǎn)隨電網(wǎng)阻抗Zg變化的分布

    圖8顯示,當(dāng)=0.03kg·m2、p=5、q=0.001時(shí),若保持=0W且v=4mH不變,只有當(dāng)g取值為3~6mH時(shí),電磁動(dòng)態(tài)才會(huì)引起100Hz附近的功率諧振。然而,若保持=0W且g=6mH不變,基于序?qū)Ъ{超SFR極點(diǎn)分布的分析結(jié)果表明(見附圖3a),只有當(dāng)v取值為4~7mH時(shí),超SFR極點(diǎn)位于右半平面,系統(tǒng)產(chǎn)生100Hz左右的功率諧振。由此可得,對(duì)多環(huán)控制型VSG-VCI來說,當(dāng)g≈v時(shí),VSG-VCI并網(wǎng)功率易產(chǎn)生100Hz附近的超SFR現(xiàn)象。此外,若保持g=0W且g=v=4mH不變,改變v,基于序?qū)Ъ{超SFR極點(diǎn)分布的分析結(jié)果表明(見附圖3b),當(dāng)v≥0.2W時(shí),超SFR極點(diǎn)位于左半平面,將其與圖8b對(duì)比可知,v與g的阻尼效果近似相同,表明當(dāng)多環(huán)控制型VSG-VCI在近似純感性的電網(wǎng)中運(yùn)行時(shí),可通過虛擬電阻來有效抑制超SFR。

    進(jìn)一步地,當(dāng)v=4mH、g=3mH且=0W時(shí),若改變功率環(huán)的控制參數(shù),基于序?qū)Ъ{超SFR極點(diǎn)分布的分析結(jié)果表明(見附圖4),減小q有助于抑制超SFR,而和p變化對(duì)超SFR的影響較小,這與單環(huán)控制型VSG-VCI的功率諧振特性有所不同(見附圖1)。究其原因,單環(huán)控制VSG功率諧振頻率約為50Hz,而圖8顯示多環(huán)控制VSG易產(chǎn)生100Hz左右的功率諧振,后者的功率諧振頻率遠(yuǎn)超出有功環(huán)控制參數(shù)的作用范圍。

    3.2 基于序阻抗的SFR阻尼特性分析

    根據(jù)以上對(duì)SFR現(xiàn)象的分析結(jié)果可知,等效線路阻尼與SFR有直接的關(guān)系,而VSG功率環(huán)控制會(huì)影響阻尼效果,進(jìn)而產(chǎn)生惡化或抑制SFR的作用。本節(jié)將通過對(duì)序阻抗“負(fù)阻尼”特性的研究,更清晰直觀地揭示VSG-VCI的同步頻率諧振機(jī)理。由于單環(huán)和多環(huán)控制型VSG-VCI的SFR現(xiàn)象均與電網(wǎng)阻抗g的動(dòng)態(tài)有關(guān),下面將通過分析不同電網(wǎng)阻抗下等效序阻抗的頻率特性,對(duì)比單環(huán)和多環(huán)控制VSG之間SFR“負(fù)阻尼”特性的差異。

    圖9所示為單環(huán)控制型VSG-VCI并網(wǎng)等效序阻抗evp()=sp()+gp()隨g變化的頻率特性曲線。圖中,灰色陰影區(qū)相位大于90°,代表負(fù)阻尼區(qū)。圖9a所示為當(dāng)g=0W、v=0.1W時(shí),evp()隨g變化的頻率特性曲線,可見在低頻范圍內(nèi),阻抗均具有正阻尼特性,雖然其相位隨g增加而增加,但仍小于90°。圖9b所示為當(dāng)g=0mH、v=0W時(shí),evp()隨g變化的頻率特性曲線,由圖可知,當(dāng)g=0.05W時(shí),序阻抗存在低頻負(fù)阻尼區(qū),表明其低頻電阻分量為負(fù),VSG-VCI并網(wǎng)存在SFR問題;隨g增加,序阻抗在低頻段由負(fù)阻尼區(qū)向正阻尼區(qū)轉(zhuǎn)移,表明低頻電阻分量由負(fù)變正,SFR得以抑制。因此,單環(huán)控制型VSG-VCI的SFR“負(fù)阻尼”以低頻負(fù)電阻的形式存在,會(huì)導(dǎo)致VSG-VCI并網(wǎng)功率產(chǎn)生50Hz左右的諧振(v=0mH),該圖對(duì)SFR的分析結(jié)果與附圖2一致。

    圖10所示為當(dāng)v=0W、v=4mH時(shí),多環(huán)控制型VSG-VCI的并網(wǎng)等效序阻抗evp()=mp()+gp()隨電網(wǎng)阻抗g變化的頻率特性曲線。圖中,相位大于90°或小于-90°的灰色陰影區(qū),代表負(fù)阻尼區(qū)。如圖10所示,evp()在50Hz附近存在諧振峰,若在諧振頻率處相位變化超過180°,則代表系統(tǒng)阻尼

    圖9 單環(huán)控制型VSG-VCI的Zevp隨Zg變化的頻率特性

    不足,即具有“負(fù)阻尼”特性,從而引發(fā)系統(tǒng)產(chǎn)生100Hz左右的功率諧振。由圖10a可得,當(dāng)g處于4mH左右(包括3mH, 4mH, 5mH, 6mH)時(shí),相位發(fā)生超出180°或近似180°的跳變,表明系統(tǒng)存在諧振問題,且諧振頻率隨g增大而減小。如圖10b所示,電網(wǎng)電阻g可為諧振峰提供有效的阻尼,且g越大,阻尼效果越好。因此,多環(huán)控制型VSG-VCI的SFR“負(fù)阻尼”以諧振峰的形式存在,會(huì)導(dǎo)致VSG- VCI并網(wǎng)功率產(chǎn)生100Hz左右的諧振,當(dāng)電網(wǎng)電感與虛擬電感相差較大時(shí),功率諧振頻率N(g+v)/g離100Hz較遠(yuǎn),對(duì)應(yīng)序阻抗離50Hz較遠(yuǎn),序阻抗呈正電阻特性,表明系統(tǒng)阻尼充分。該圖對(duì)SFR的分析結(jié)果與圖8一致。

    圖10 多環(huán)控制型VSG-VCI的Zevp隨Zg變化的頻率特性

    4 實(shí)驗(yàn)驗(yàn)證

    本文通過一臺(tái)5kW的并網(wǎng)逆變器對(duì)SFR現(xiàn)象進(jìn)行實(shí)驗(yàn)驗(yàn)證,采用基于RT-Box的硬件在環(huán)實(shí)驗(yàn)平臺(tái)如圖11所示,即通過圖中的RT-Box運(yùn)行并網(wǎng)逆變器的主電路部分,控制部分采用芯片TMS320F28335。其中,RT-Box中運(yùn)行的主電路來源于Plecs軟件中搭建的三相并網(wǎng)逆變器功率電路模型(見圖1),主電路的離散時(shí)間為10ms,開關(guān)頻率和控制電路的采樣頻率均為10kHz,其他實(shí)驗(yàn)系統(tǒng)參數(shù)與表1保持一致。

    圖11 基于RT-Box的硬件在環(huán)實(shí)驗(yàn)平臺(tái)

    實(shí)驗(yàn)通過采集逆變器并網(wǎng)有功功率e的波形來監(jiān)測(cè)功率諧振現(xiàn)象,將在RT-Box中實(shí)時(shí)計(jì)算得到的e傳送至模擬輸出口,進(jìn)而通過轉(zhuǎn)接板和刺刀螺母連接器(BNC)接口的傳輸線將功率數(shù)據(jù)傳輸?shù)绞静ㄆ髦羞M(jìn)行觀測(cè)分析。

    4.1 單環(huán)控制型VSG-VCI的有功輸出波形

    圖12所示為理想電網(wǎng)下,單環(huán)控制型VSG-VCI有功功率e隨控制環(huán)參數(shù)變化的實(shí)驗(yàn)波形。由圖12a可見,當(dāng)虛擬電阻v=0.05W時(shí),e發(fā)生諧振,且由放大波形可知諧振周期為20ms,表明e的諧振頻率為50Hz,系統(tǒng)存在SFR問題(圖12所示波形是因?qū)嶒?yàn)平臺(tái)模擬輸出限幅所致,不影響SFR分析結(jié)果);當(dāng)v增加到0.1W后,SFR現(xiàn)象消失,系統(tǒng)輸出穩(wěn)定的有功功率波形。當(dāng)v=0.05W時(shí),若將虛擬電感v的值由0mH改變?yōu)?mH,如圖12b所示,e的諧振周期減小到約10ms,表明e的諧振頻率約為100Hz,這與圖7闡明的現(xiàn)象相符。當(dāng)v=0.1W、v=0mH時(shí),由圖12c和圖12e可知,減小慣性系數(shù)或增加無功下垂系數(shù)q會(huì)影響VSG-VCI并網(wǎng)功率的穩(wěn)定性,即使0.1W的虛擬電阻可提供抑制SFR的阻尼(見圖12a),但當(dāng)減小到0.01或q增加到0.002時(shí),e波形再次出現(xiàn)SFR現(xiàn)象。由圖12d可見,有功阻尼系數(shù)p的大范圍變化不會(huì)影響系統(tǒng)的穩(wěn)定性。實(shí)驗(yàn)結(jié)果驗(yàn)證了理論分析的準(zhǔn)確性。

    實(shí)際電網(wǎng)中,同理可以通過實(shí)驗(yàn)驗(yàn)證電網(wǎng)阻抗g變化對(duì)SFR現(xiàn)象的影響,這里不再贅述。

    圖12 理想電網(wǎng)下,單環(huán)控制型VSG-VCI隨控制參數(shù)變化的有功功率

    4.2 多環(huán)控制型VSG-VCI的有功輸出波形

    保持功率環(huán)控制參數(shù)=0.03、p=5、q=0.001不變,圖13所示為當(dāng)v=0W、v=4mH時(shí),多環(huán)控制型VSG-VCI的有功功率e隨電網(wǎng)阻抗g變化的實(shí)驗(yàn)波形;圖14所示為當(dāng)g=0W、g=6mH時(shí),多環(huán)控制型VSG-VCI的有功功率e隨虛擬阻抗v變化的實(shí)驗(yàn)波形。

    當(dāng)g=0W時(shí),由圖13可知,當(dāng)g為2mH和7mH時(shí),即使系統(tǒng)阻尼=v+g=0W,e波形仍可保持穩(wěn)定;當(dāng)g為3~6mH時(shí),由放大的功率波形可見,e存在100Hz的諧振。此外,由圖13d可知,加入0.3W的電網(wǎng)電阻g可抑制功率諧振。實(shí)驗(yàn)結(jié)果驗(yàn)證了圖8所示的理論分析。

    圖13 多環(huán)控制型VSG-VCI隨Zg變化的有功波形

    圖14 多環(huán)控制型VSG-VCI隨Zv變化的有功波形

    由圖14a可見,當(dāng)v=2mH和v=3mH時(shí),e波形穩(wěn)定在5kW;當(dāng)v=4mH和v=5mH時(shí),功率失穩(wěn),加入0.2W的虛擬電阻v可使功率波形重新恢復(fù)穩(wěn)定。實(shí)驗(yàn)結(jié)果驗(yàn)證了附圖3所示的理論分析。此外,結(jié)合圖13的實(shí)驗(yàn)結(jié)果可證明,對(duì)多環(huán)控制型VSG-VCI并網(wǎng)系統(tǒng)來說,只有當(dāng)g和v的值接近時(shí),系統(tǒng)才容易產(chǎn)生因電磁動(dòng)態(tài)引發(fā)的功率諧振現(xiàn)象,且諧振頻率約為100Hz。

    同理可以通過實(shí)驗(yàn)驗(yàn)證功率環(huán)控制參數(shù)變化對(duì)超SFR現(xiàn)象的影響,這里不再贅述。

    5 結(jié)論

    針對(duì)VSG電磁動(dòng)態(tài)引發(fā)的SFR問題,為了對(duì)其進(jìn)行有效的建模分析且將其統(tǒng)一于VSG寬頻域穩(wěn)定性的研究中,本文提出了基于單輸入單輸出序阻抗模型的SFR分析方法,并全面對(duì)比研究了單環(huán)和多環(huán)控制VSG的SFR特性與功率環(huán)控制參數(shù)、虛擬阻抗、電網(wǎng)阻抗之間的關(guān)系。研究表明:

    1)虛擬阻抗和內(nèi)環(huán)控制的引入會(huì)極大地改變SFR表現(xiàn)形式。單環(huán)控制VSG的功率諧振頻率會(huì)隨v的增大而增大,而多環(huán)控制VSG的功率諧振頻率均在100Hz附近。這是因?yàn)閮?nèi)環(huán)控制的引入改變了VSG-VCI的輸出阻抗特性,引發(fā)了50Hz附近的諧振峰,從而改變了系統(tǒng)的負(fù)阻尼特性。

    2)單環(huán)控制VSG的SFR與濾波線路動(dòng)態(tài)密切相關(guān),在VSG-VCI自身設(shè)計(jì)參數(shù)合理、阻尼充分的情況下,電網(wǎng)阻抗變化不會(huì)影響SFR;而多環(huán)控制VSG的SFR與濾波線路動(dòng)態(tài)無關(guān)、與電網(wǎng)阻抗動(dòng)態(tài)密切相關(guān),但是只有當(dāng)虛擬電感與電網(wǎng)電感值接近時(shí),系統(tǒng)才會(huì)出現(xiàn)100Hz附近的功率諧振。

    3)虛擬電阻v與電網(wǎng)電阻g都可有效阻尼SFR,且同樣阻值下它們的阻尼效果近似相同。因此,在近似純感性的電網(wǎng)中,可通過引入虛擬電阻v的方式有效抑制SFR。

    附圖1sp的SFR極點(diǎn)隨功率環(huán)控制參數(shù)變化的分布

    App.Fig.1 Map of the SFR poles ofspwith varied power control parameters

    附圖2evp的SFR極點(diǎn)隨電網(wǎng)阻抗g變化的分布

    App.Fig.2 Map of the SFR poles ofevpwith variedg

    附圖3evp的超SFR極點(diǎn)隨虛擬阻抗v變化的分布

    App.Fig.3 Map of the super-SFR poles ofevpwith variedv

    附圖4evp的超SFR極點(diǎn)隨功率環(huán)控制參數(shù)變化的分布

    App.Fig.4 Map of the super-SFR poles ofevpwith varied power control parameters

    [1] 姜靜雅, 王瑋, 吳學(xué)智, 等. 基于自適應(yīng)無功功率補(bǔ)償?shù)奶摂M同步機(jī)功率解耦策略[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(13): 2747-2756.

    Jiang Jingya, Wang Wei, Wu Xuezhi, et al. Power decoupling strategy in virtual synchronous generator based on adaptive reactive power compensation[J]. Transactions of China Electrotechnical Society, 2020, 35(13): 2747-2756.

    [2] 章艷, 高晗, 張萌. 不同虛擬同步機(jī)控制下雙饋風(fēng)機(jī)系統(tǒng)頻率響應(yīng)差異研究[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(13): 2889-2900.

    Zhang Yan, Gao Han, Zhang Meng. Research on frequency response difference of DFIG system controlled by different virtual synchronous generator controls[J]. Transactions of China Electrotechnical Society, 2020, 35(13): 2889-2900.

    [3] 張琛, 蔡旭, 李征. 具有自主電網(wǎng)同步與弱網(wǎng)穩(wěn)定運(yùn)行能力的雙饋風(fēng)電機(jī)組控制方法[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2017, 37(2): 476-485.

    Zhang Chen, Cai Xu, Li Zheng. Control of DFIG- based wind turbines with the capability of automatic grid-synchronization and stable operation under weak grid condition[J]. Proceedings of the CSEE, 2017, 37(2): 476-485.

    [4] Liu Jia, Miura Y, Ise T. Comparison of dynamic characteristics between virtual synchronous generator and droop control in inverter-based distributed gener- ators[J]. IEEE Transactions on Power Electronics, 2016, 31(5): 3600-3611.

    [5] 顏湘武, 賈焦心, 王德勝, 等. 基于/“導(dǎo)納”的并聯(lián)虛擬同步機(jī)功頻響應(yīng)建模與分析[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(15): 3191-3202.

    Yan Xiangwu, Jia Jiaoxin, Wang Desheng, et al. Modeling and analysis of active power-frequency response of parallel VSGs using a/“Admittance”[J]. Transactions of China Electrotechnical Society, 2020, 35(15): 3191-3202.

    [6] Chen Junru, O'Donnell T. Parameter constraints for virtual synchronous generator considering stability[J]. IEEE Transactions on Power Systems, 2019, 34(3): 2479-2481.

    [7] 李武華, 王金華, 楊賀雅, 等. 虛擬同步發(fā)電機(jī)的功率動(dòng)態(tài)耦合機(jī)理及同步頻率諧振抑制策略[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2017, 37(2): 381-390.

    Li Wuhua, Wang Jinhua, Yang Heya, et al. Power dynamic coupling mechanism and resonance suppression of synchronous frequency for virtual synchronous generators[J]. Proceedings of the CSEE, 2017, 37(2): 381-390.

    [8] Yang Dongsheng, Wu Heng, Wang Xiongfei, et al. Suppression of synchronous resonance for VSGs[J]. The Journal of Engineering, 2017, 13(1): 2574- 2579.

    [9] Li Xin, Hu Yaowei, Shao Yuting, et al. Mechanism analysis and suppression strategies of power oscilla- tion for virtual synchronous generator[C]//IECON 2017-43rd Annual Conference of the IEEE Industrial Electronics Society, Beijing, China, 2017: 4955-4960.

    [10] 張行, 張興, 李明, 等. 基于暫態(tài)虛擬電阻VSG同步頻率諧振抑制策略[J]. 電力電子技術(shù), 2019, 53(3): 68-70, 81.

    Zhang Hang, Zhang Xing, Li Ming, et al. Synchronous frequency resonance suppression strategy for virtual synchronous generators based on transient virtual resistance[J]. Power Electronics, 2019, 53(3): 68-70, 81.

    [11] 涂春鳴, 謝偉杰, 肖凡, 等. 控制延時(shí)對(duì)虛擬同步機(jī)全局穩(wěn)定性的影響分析[J]. 電網(wǎng)技術(shù), 2020, 44(8): 3021-3028.

    Tu Chunming, Xie Weijie, Xiao Fan, et al. Analysis of influence of control delay on stability of virtual synchronous machine[J]. Power System Technology, 2020, 44(8): 3021-3028.

    [12] 李輝, 王坤, 胡玉, 等. 雙饋風(fēng)電系統(tǒng)虛擬同步控制的阻抗建模及穩(wěn)定性分析[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2019, 39(12): 3434-3442.

    Li Hui, Wang Kun, Hu Yu, et al. Impedance modeling and stability analysis of virtual synchronous control based on doubly-fed wind generation systems[J]. Proceedings of the CSEE, 2019, 39(12): 3434-3442.

    [13] 涂春鳴, 謝偉杰, 肖凡, 等. 多虛擬同步發(fā)電機(jī)并聯(lián)系統(tǒng)控制參數(shù)對(duì)穩(wěn)定性的影響分析[J]. 電力系統(tǒng)自動(dòng)化, 2020, 44(15): 77-86.

    Tu Chunming, Xie Weijie, Xiao Fan, et al. Influence analysis of control parameters of parallel system with multiple virtual synchronous generators on stability[J]. Automation of Electric Power Systems, 2020, 44(15): 77-86.

    [14] Cespedes M, Sun Jian. Adaptive control of grid- connected inverters based on online grid impedance measurements[J]. IEEE Transactions on Sustainable Energy, 2014, 5(2): 516-523.

    [15] Wu Wenhua, Chen Yandong, Zhou Leming, et al. Sequence impedance modeling and stability com- parative analysis of voltage-controlled VSGs and current-controlled VSGs[J]. IEEE Transactions on Industrial Electronics, 2019, 66(8): 6460-6472.

    Synchronous Frequency Resonance in Virtual Synchronous Generator Based on Sequence-Impedance

    1221

    (1. School of Electrical Engineering & Automation Harbin Institute of Technology Harbin 150001 China 2. Electric Power Research Institute of State Grid Heilongjiang Electric Power Co. Ltd Harbin 150030 China)

    As for the synchronous frequency resonance (SFR) phenomenon existing in virtual synchronous generator (VSG), considering the shortcomings of multi-input multi-output power coupling SFR analysis models, a sequence-admittance and sequence-impedance based SFR analysis method is proposed by modeling the single-input single-output sequence impedance model that can describe the wide-frequency dynamics of the VSG, and the SFR phenomena of single-loop and multi-loop VSG are comparatively studied. Based on the pole distribution map of sequence admittance, the relationships between SFR and the VSG power control loop parameters, virtual impedance, inner-loop control and grid impedance are effectively studied. Based on the Bode diagrams of sequence-impedance, the damping characteristics of SFR with and without the inner-loop control are compared, and then the resonance mechanism of SFR in single-loop and multi-loop VSG is revealed from the perspective of negative damping. By building the experimental platform based on the RT-Box, the SFR phenomena of VSG under different control schemes and grid impedances are verified, indicating the optimal design direction of VSG control.

    Virtual synchronous generator, synchronous frequency resonance, sequence-impedance, negative damping, inner-loop control

    10.19595/j.cnki.1000-6753.tces.210697

    TM464

    于彥雪 女,1991年生,博士研究生,研究方向?yàn)槿蹼娋W(wǎng)下并網(wǎng)逆變器建模及穩(wěn)定性研究。E-mail: yuyanxueppkz@163.com

    李浩昱 男,1974年生,教授,博士生導(dǎo)師,研究方向?yàn)槲⒕W(wǎng)環(huán)境變換器協(xié)調(diào)控制技術(shù)、極端環(huán)境電能變換與控制。E-mail: lihy@hit.edu.cn(通信作者)

    2021-05-18

    2021-12-10

    黑龍江電力公司科技資助項(xiàng)目(522437200038)。

    (編輯 陳 誠(chéng))

    猜你喜歡
    環(huán)控制諧振阻尼
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項(xiàng)的擴(kuò)散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    基于諧振開關(guān)技術(shù)的低相噪LC VCO的設(shè)計(jì)
    自適應(yīng)模糊多環(huán)控制在慣導(dǎo)平臺(tái)穩(wěn)定回路中的應(yīng)用
    諧振式單開關(guān)多路輸出Boost LED驅(qū)動(dòng)電源
    具阻尼項(xiàng)的Boussinesq型方程的長(zhǎng)時(shí)間行為
    滯環(huán)控制的變步長(zhǎng)MPPT算法實(shí)驗(yàn)研究
    基于CM6901 的LLC半橋諧振開關(guān)電源設(shè)計(jì)
    一種雙饋風(fēng)電機(jī)組網(wǎng)側(cè)變流器抑制低次諧波的多環(huán)控制
    亚洲性夜色夜夜综合| 中文字幕av电影在线播放| 色老头精品视频在线观看| 纯流量卡能插随身wifi吗| 国产淫语在线视频| 叶爱在线成人免费视频播放| 极品人妻少妇av视频| 精品一区二区三区av网在线观看 | 人人妻,人人澡人人爽秒播| 亚洲,欧美精品.| 99国产精品一区二区蜜桃av | av电影中文网址| 亚洲专区中文字幕在线| 国产在线免费精品| 国产国语露脸激情在线看| 男女高潮啪啪啪动态图| 老熟妇乱子伦视频在线观看| 久久亚洲精品不卡| 麻豆av在线久日| 欧美日韩亚洲国产一区二区在线观看 | 日本vs欧美在线观看视频| 99精品在免费线老司机午夜| 天天躁夜夜躁狠狠躁躁| 国产男女超爽视频在线观看| 两人在一起打扑克的视频| 真人做人爱边吃奶动态| 97在线人人人人妻| 黄频高清免费视频| 精品欧美一区二区三区在线| 欧美精品啪啪一区二区三区| 一区福利在线观看| 一边摸一边抽搐一进一小说 | 黄色成人免费大全| 国产91精品成人一区二区三区 | 国产高清视频在线播放一区| 色94色欧美一区二区| 欧美人与性动交α欧美精品济南到| 久久国产精品人妻蜜桃| 国产日韩欧美视频二区| 欧美黑人欧美精品刺激| 99精品久久久久人妻精品| 亚洲精品久久午夜乱码| 每晚都被弄得嗷嗷叫到高潮| 国产精品.久久久| 免费观看a级毛片全部| 欧美久久黑人一区二区| 18禁观看日本| 俄罗斯特黄特色一大片| 色综合欧美亚洲国产小说| 热99国产精品久久久久久7| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人| 国产成人啪精品午夜网站| 成年人午夜在线观看视频| 国产精品麻豆人妻色哟哟久久| 99精国产麻豆久久婷婷| 777米奇影视久久| 99国产精品一区二区三区| 精品亚洲成国产av| av网站免费在线观看视频| 老司机福利观看| 女性被躁到高潮视频| 巨乳人妻的诱惑在线观看| 亚洲精品中文字幕一二三四区 | 欧美激情 高清一区二区三区| 久久精品91无色码中文字幕| 黄色丝袜av网址大全| 亚洲精品成人av观看孕妇| 国产av一区二区精品久久| 亚洲av日韩在线播放| 一级片'在线观看视频| 国产精品熟女久久久久浪| 国产一区有黄有色的免费视频| 少妇裸体淫交视频免费看高清 | 久久人妻熟女aⅴ| 亚洲精品乱久久久久久| 久久久久网色| 欧美老熟妇乱子伦牲交| 人人妻人人澡人人看| 国产免费视频播放在线视频| 国产成人系列免费观看| 久久人妻熟女aⅴ| 色综合欧美亚洲国产小说| 下体分泌物呈黄色| 国产欧美日韩精品亚洲av| 欧美日韩视频精品一区| 精品少妇内射三级| 免费黄频网站在线观看国产| a级毛片在线看网站| 正在播放国产对白刺激| 黄色 视频免费看| 在线天堂中文资源库| 免费在线观看黄色视频的| 国产1区2区3区精品| 久久精品国产a三级三级三级| 久久国产精品大桥未久av| 精品久久久久久电影网| 国产老妇伦熟女老妇高清| 午夜福利乱码中文字幕| tocl精华| 国产成+人综合+亚洲专区| 亚洲欧美色中文字幕在线| 亚洲av第一区精品v没综合| 欧美一级毛片孕妇| 亚洲欧美激情在线| 90打野战视频偷拍视频| 亚洲免费av在线视频| av不卡在线播放| 国产淫语在线视频| 久久久国产成人免费| 久久 成人 亚洲| 亚洲全国av大片| 露出奶头的视频| 中文字幕色久视频| 一区二区三区激情视频| 免费高清在线观看日韩| 国产成人啪精品午夜网站| 在线观看免费视频日本深夜| 欧美激情极品国产一区二区三区| aaaaa片日本免费| 欧美人与性动交α欧美软件| 1024香蕉在线观看| 69精品国产乱码久久久| 大片免费播放器 马上看| 一个人免费看片子| 亚洲av国产av综合av卡| 免费一级毛片在线播放高清视频 | 亚洲精品国产一区二区精华液| 变态另类成人亚洲欧美熟女 | 美国免费a级毛片| 岛国毛片在线播放| 免费不卡黄色视频| 亚洲欧美日韩另类电影网站| 久久久国产精品麻豆| 一本色道久久久久久精品综合| 成年人午夜在线观看视频| kizo精华| 性少妇av在线| 亚洲成人免费电影在线观看| 91成年电影在线观看| 国产成人欧美在线观看 | 两个人免费观看高清视频| 一区二区三区精品91| 久久99热这里只频精品6学生| 久久久久久久国产电影| 狠狠狠狠99中文字幕| www.熟女人妻精品国产| 悠悠久久av| 俄罗斯特黄特色一大片| 人妻久久中文字幕网| 脱女人内裤的视频| av线在线观看网站| 美女高潮到喷水免费观看| 黄片大片在线免费观看| 美女福利国产在线| 男女无遮挡免费网站观看| 国产黄色免费在线视频| 亚洲国产欧美网| 精品一区二区三卡| 别揉我奶头~嗯~啊~动态视频| 99久久99久久久精品蜜桃| 日韩 欧美 亚洲 中文字幕| 这个男人来自地球电影免费观看| 色播在线永久视频| 蜜桃在线观看..| 女警被强在线播放| 国产亚洲精品一区二区www | 国产精品成人在线| 国产精品久久久久成人av| 亚洲视频免费观看视频| e午夜精品久久久久久久| 丁香六月天网| 两性夫妻黄色片| 午夜91福利影院| 久久亚洲精品不卡| 日韩一卡2卡3卡4卡2021年| 午夜激情av网站| 捣出白浆h1v1| 人妻久久中文字幕网| 久久久久久免费高清国产稀缺| 露出奶头的视频| 欧美中文综合在线视频| 国产深夜福利视频在线观看| 欧美日韩成人在线一区二区| 又黄又粗又硬又大视频| 精品乱码久久久久久99久播| 国产成人欧美在线观看 | 国产欧美日韩精品亚洲av| 国产区一区二久久| 国产成人欧美| 久久久久久久久久久久大奶| 午夜日韩欧美国产| 亚洲成人国产一区在线观看| 成年女人毛片免费观看观看9 | 黑人巨大精品欧美一区二区蜜桃| 汤姆久久久久久久影院中文字幕| 又黄又粗又硬又大视频| 国产精品九九99| 大码成人一级视频| 18禁观看日本| 欧美日韩成人在线一区二区| a级毛片黄视频| 成年人免费黄色播放视频| 18禁裸乳无遮挡动漫免费视频| 侵犯人妻中文字幕一二三四区| 性高湖久久久久久久久免费观看| 在线播放国产精品三级| 涩涩av久久男人的天堂| 美女主播在线视频| 麻豆成人av在线观看| 国产成人精品在线电影| 9色porny在线观看| 久久精品熟女亚洲av麻豆精品| 欧美人与性动交α欧美精品济南到| 国产精品熟女久久久久浪| 天堂动漫精品| 国产精品免费大片| 丁香六月天网| 国产精品美女特级片免费视频播放器 | 女性生殖器流出的白浆| 久久影院123| 国产日韩欧美亚洲二区| 成人黄色视频免费在线看| 色婷婷久久久亚洲欧美| 日韩熟女老妇一区二区性免费视频| 国产亚洲欧美精品永久| av线在线观看网站| 丝袜人妻中文字幕| 不卡av一区二区三区| 91精品三级在线观看| 久久久久视频综合| 自线自在国产av| 久久久久国内视频| 国产精品熟女久久久久浪| 亚洲精品美女久久av网站| 美女扒开内裤让男人捅视频| 十八禁高潮呻吟视频| 18在线观看网站| 亚洲av成人一区二区三| www.999成人在线观看| 婷婷丁香在线五月| 亚洲伊人色综图| 18在线观看网站| 又黄又粗又硬又大视频| 国产精品一区二区在线观看99| 在线播放国产精品三级| 日韩 欧美 亚洲 中文字幕| 国产精品美女特级片免费视频播放器 | 黄色视频,在线免费观看| 欧美大码av| 久久九九热精品免费| 欧美国产精品va在线观看不卡| 国产一区二区三区综合在线观看| 法律面前人人平等表现在哪些方面| www.自偷自拍.com| 一区在线观看完整版| 一边摸一边抽搐一进一出视频| 国产亚洲欧美在线一区二区| 午夜免费鲁丝| 宅男免费午夜| 亚洲国产欧美在线一区| 丁香六月天网| 菩萨蛮人人尽说江南好唐韦庄| 午夜两性在线视频| 色精品久久人妻99蜜桃| 国产精品熟女久久久久浪| 国产亚洲欧美在线一区二区| 久久99热这里只频精品6学生| 久久亚洲精品不卡| 十分钟在线观看高清视频www| 美女视频免费永久观看网站| 亚洲黑人精品在线| www日本在线高清视频| 岛国毛片在线播放| 国产人伦9x9x在线观看| 国产精品香港三级国产av潘金莲| 精品久久蜜臀av无| 免费看a级黄色片| 香蕉国产在线看| 下体分泌物呈黄色| 大型av网站在线播放| 美女福利国产在线| 美女高潮到喷水免费观看| 最新美女视频免费是黄的| 成年动漫av网址| 亚洲国产av新网站| 成人特级黄色片久久久久久久 | 一区二区三区乱码不卡18| 久久久久精品人妻al黑| 最新美女视频免费是黄的| 精品视频人人做人人爽| 国产97色在线日韩免费| www.自偷自拍.com| bbb黄色大片| 国产一区二区三区视频了| 国产精品 国内视频| 超碰成人久久| 免费观看av网站的网址| 又黄又粗又硬又大视频| 久久人妻福利社区极品人妻图片| 9191精品国产免费久久| 亚洲av日韩精品久久久久久密| 国产日韩一区二区三区精品不卡| 欧美午夜高清在线| av一本久久久久| 国产精品久久久av美女十八| 夜夜夜夜夜久久久久| 少妇粗大呻吟视频| 国产熟女午夜一区二区三区| 999久久久精品免费观看国产| 久久精品91无色码中文字幕| 国产精品免费视频内射| 国产成人精品久久二区二区91| 咕卡用的链子| 亚洲九九香蕉| 一本久久精品| 久久中文字幕人妻熟女| 日日摸夜夜添夜夜添小说| 国产精品免费一区二区三区在线 | 国产精品电影一区二区三区 | 亚洲熟女精品中文字幕| 欧美日韩黄片免| 日韩大码丰满熟妇| 精品少妇一区二区三区视频日本电影| 欧美在线一区亚洲| 成年人午夜在线观看视频| 中文字幕人妻丝袜制服| 一本久久精品| 黄频高清免费视频| 久久精品亚洲精品国产色婷小说| 中文字幕人妻丝袜制服| av不卡在线播放| 免费观看人在逋| 亚洲欧美一区二区三区久久| 黄色视频,在线免费观看| 精品高清国产在线一区| 90打野战视频偷拍视频| 99久久99久久久精品蜜桃| 国产亚洲欧美在线一区二区| 亚洲中文日韩欧美视频| 黑丝袜美女国产一区| 午夜福利视频精品| 日韩熟女老妇一区二区性免费视频| 啦啦啦中文免费视频观看日本| 高清视频免费观看一区二区| 国产成人av激情在线播放| 黄片小视频在线播放| 国产片内射在线| 国产精品香港三级国产av潘金莲| 多毛熟女@视频| av欧美777| av又黄又爽大尺度在线免费看| 极品少妇高潮喷水抽搐| 国产成人免费观看mmmm| 国产精品久久久av美女十八| 国产亚洲欧美精品永久| 午夜福利一区二区在线看| 人妻久久中文字幕网| 久久国产精品男人的天堂亚洲| 少妇的丰满在线观看| 午夜福利视频在线观看免费| 777米奇影视久久| 久久国产精品男人的天堂亚洲| 精品久久蜜臀av无| h视频一区二区三区| 亚洲综合色网址| 在线观看一区二区三区激情| 午夜激情久久久久久久| av欧美777| 国产亚洲午夜精品一区二区久久| 欧美变态另类bdsm刘玥| 国产高清国产精品国产三级| 1024香蕉在线观看| 日本av手机在线免费观看| 深夜精品福利| 久久 成人 亚洲| 国产精品.久久久| 在线观看免费高清a一片| 国产一区二区激情短视频| 国产1区2区3区精品| 日韩熟女老妇一区二区性免费视频| 少妇猛男粗大的猛烈进出视频| 国产欧美亚洲国产| 亚洲成人手机| 中文字幕另类日韩欧美亚洲嫩草| 久久久久网色| 18在线观看网站| 日韩中文字幕欧美一区二区| 亚洲色图 男人天堂 中文字幕| 婷婷成人精品国产| 午夜精品久久久久久毛片777| 女人被躁到高潮嗷嗷叫费观| 国产麻豆69| 十八禁高潮呻吟视频| 在线av久久热| 精品一品国产午夜福利视频| 在线播放国产精品三级| kizo精华| 国产99久久九九免费精品| 99在线人妻在线中文字幕 | 亚洲专区中文字幕在线| 桃花免费在线播放| 欧美av亚洲av综合av国产av| 亚洲avbb在线观看| 亚洲成人国产一区在线观看| 啦啦啦中文免费视频观看日本| 大香蕉久久网| 在线观看舔阴道视频| 国产亚洲精品久久久久5区| 免费在线观看视频国产中文字幕亚洲| 欧美精品一区二区大全| 黄网站色视频无遮挡免费观看| 中文欧美无线码| 老司机靠b影院| 亚洲人成电影免费在线| 精品亚洲乱码少妇综合久久| 成人三级做爰电影| 久久久欧美国产精品| 国产日韩一区二区三区精品不卡| 亚洲精品国产精品久久久不卡| 热re99久久国产66热| 黑丝袜美女国产一区| 我要看黄色一级片免费的| 男女免费视频国产| 一级黄色大片毛片| 波多野结衣一区麻豆| 亚洲专区中文字幕在线| 久久精品国产99精品国产亚洲性色 | 中文字幕人妻丝袜制服| 国产精品99久久99久久久不卡| 亚洲人成电影观看| 999久久久国产精品视频| 国产麻豆69| 丁香欧美五月| 在线播放国产精品三级| 99在线人妻在线中文字幕 | 国产亚洲精品第一综合不卡| 午夜激情久久久久久久| 亚洲,欧美精品.| 久久青草综合色| 19禁男女啪啪无遮挡网站| 精品卡一卡二卡四卡免费| 久久人妻福利社区极品人妻图片| 老熟妇乱子伦视频在线观看| 国产色视频综合| 亚洲 国产 在线| 亚洲免费av在线视频| 男女无遮挡免费网站观看| av天堂久久9| www.999成人在线观看| 亚洲中文日韩欧美视频| 一夜夜www| 在线 av 中文字幕| 大香蕉久久成人网| 丝袜在线中文字幕| 一本—道久久a久久精品蜜桃钙片| 99re在线观看精品视频| av天堂久久9| 岛国毛片在线播放| 久久精品亚洲熟妇少妇任你| 亚洲一区二区三区欧美精品| 成人av一区二区三区在线看| 最新在线观看一区二区三区| 日韩欧美一区二区三区在线观看 | 精品国产国语对白av| 亚洲一码二码三码区别大吗| 中文亚洲av片在线观看爽 | 国产成人精品久久二区二区免费| 亚洲精品中文字幕一二三四区 | 国产麻豆69| 99re在线观看精品视频| 国产成人av激情在线播放| 成人亚洲精品一区在线观看| 国产精品国产高清国产av | 精品国产乱子伦一区二区三区| 午夜福利免费观看在线| 欧美日韩精品网址| 午夜精品久久久久久毛片777| 精品亚洲成a人片在线观看| 欧美另类亚洲清纯唯美| 亚洲久久久国产精品| 亚洲中文日韩欧美视频| 人成视频在线观看免费观看| 国产成人欧美| 日韩中文字幕视频在线看片| 亚洲第一青青草原| 午夜91福利影院| 母亲3免费完整高清在线观看| 妹子高潮喷水视频| 免费不卡黄色视频| 另类亚洲欧美激情| 大片免费播放器 马上看| 久久久久久久精品吃奶| 飞空精品影院首页| 777米奇影视久久| 欧美黄色片欧美黄色片| 一区二区三区乱码不卡18| 一级,二级,三级黄色视频| 纯流量卡能插随身wifi吗| 日韩中文字幕视频在线看片| 人人妻人人爽人人添夜夜欢视频| 亚洲国产欧美网| 一级片免费观看大全| 悠悠久久av| 老司机午夜福利在线观看视频 | 亚洲欧洲日产国产| 大香蕉久久网| 下体分泌物呈黄色| 日韩有码中文字幕| h视频一区二区三区| 亚洲 欧美一区二区三区| 男女之事视频高清在线观看| 妹子高潮喷水视频| 狂野欧美激情性xxxx| 国产不卡av网站在线观看| 在线观看免费高清a一片| 大型av网站在线播放| 老司机靠b影院| 操美女的视频在线观看| 国产91精品成人一区二区三区 | 色精品久久人妻99蜜桃| 日韩制服丝袜自拍偷拍| 国产亚洲午夜精品一区二区久久| av有码第一页| 成年版毛片免费区| 国产国语露脸激情在线看| av欧美777| 黄色毛片三级朝国网站| 亚洲一码二码三码区别大吗| 欧美日韩亚洲高清精品| 老司机福利观看| 日本黄色日本黄色录像| 热99国产精品久久久久久7| 如日韩欧美国产精品一区二区三区| videos熟女内射| 成人18禁在线播放| 日韩视频一区二区在线观看| 麻豆乱淫一区二区| 在线av久久热| 91精品三级在线观看| 亚洲va日本ⅴa欧美va伊人久久| 欧美乱妇无乱码| 欧美精品高潮呻吟av久久| 国产精品一区二区免费欧美| 黄频高清免费视频| 亚洲国产欧美日韩在线播放| 2018国产大陆天天弄谢| 亚洲中文av在线| av免费在线观看网站| 国产极品粉嫩免费观看在线| 午夜福利一区二区在线看| 精品卡一卡二卡四卡免费| 这个男人来自地球电影免费观看| 极品少妇高潮喷水抽搐| 国产在线精品亚洲第一网站| 18禁美女被吸乳视频| 妹子高潮喷水视频| 老司机福利观看| 亚洲精品国产精品久久久不卡| 亚洲色图av天堂| 日本一区二区免费在线视频| 国产高清激情床上av| 久久国产亚洲av麻豆专区| 久久天躁狠狠躁夜夜2o2o| 十八禁高潮呻吟视频| 国产日韩一区二区三区精品不卡| 国产伦人伦偷精品视频| 老汉色∧v一级毛片| 欧美成人免费av一区二区三区 | 大型av网站在线播放| 纯流量卡能插随身wifi吗| 色尼玛亚洲综合影院| 人人妻人人澡人人爽人人夜夜| 亚洲专区国产一区二区| 极品人妻少妇av视频| 日韩欧美三级三区| 欧美人与性动交α欧美软件| 美女主播在线视频| av视频免费观看在线观看| 亚洲午夜理论影院| 国产精品成人在线| 国产真人三级小视频在线观看| 免费一级毛片在线播放高清视频 | 日韩制服丝袜自拍偷拍| 久久午夜综合久久蜜桃| 亚洲欧美一区二区三区黑人| 狠狠精品人妻久久久久久综合| 国产片内射在线| 国产亚洲精品久久久久5区| 女人被躁到高潮嗷嗷叫费观| 色尼玛亚洲综合影院| 亚洲九九香蕉| 欧美 日韩 精品 国产| www.熟女人妻精品国产| 肉色欧美久久久久久久蜜桃| 不卡一级毛片| 午夜激情久久久久久久| 国产成人精品在线电影| 999久久久精品免费观看国产| 高清av免费在线| 露出奶头的视频| 精品一区二区三区视频在线观看免费 | 另类亚洲欧美激情| 无限看片的www在线观看| 成人黄色视频免费在线看| 一区二区日韩欧美中文字幕| 黑人操中国人逼视频| 欧美乱妇无乱码| 欧美黑人精品巨大| av电影中文网址| 欧美日韩福利视频一区二区| av视频免费观看在线观看| 欧美人与性动交α欧美软件| 女警被强在线播放| 午夜福利,免费看|