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

    基于不對稱正負(fù)反饋效應(yīng)的PQ功率控制并網(wǎng)逆變器穩(wěn)定性分析

    2023-02-08 08:15:08涂春鳴鄒凱星高家元葛平娟
    電工技術(shù)學(xué)報 2023年2期
    關(guān)鍵詞:負(fù)反饋環(huán)路穩(wěn)定性

    涂春鳴 鄒凱星 高家元 肖 凡 葛平娟

    基于不對稱正負(fù)反饋效應(yīng)的PQ功率控制并網(wǎng)逆變器穩(wěn)定性分析

    涂春鳴 鄒凱星 高家元 肖 凡 葛平娟

    (湖南大學(xué)國家電能變換與控制工程技術(shù)研究中心 長沙 410082)

    針對弱電網(wǎng)下PQ功率控制并網(wǎng)逆變器的穩(wěn)定性問題,該文通過構(gòu)建有功/無功環(huán)的小信號控制框圖,首次揭示弱電網(wǎng)下PQ功率控制是一種不對稱控制的本質(zhì),即無功環(huán)控制器會在有功和無功控制環(huán)中形成不對稱的正負(fù)反饋環(huán)路。同時,基于PQ功率環(huán)小信號控制框圖的不對稱現(xiàn)象,揭示有功/無功環(huán)控制器參數(shù)變化對并網(wǎng)逆變器穩(wěn)定性影響存在差異性的機理,指出無功環(huán)控制器參數(shù)變化對并網(wǎng)系統(tǒng)穩(wěn)定性的影響程度要明顯強于有功環(huán)控制器參數(shù)變化。在此基礎(chǔ)上,借助有功與無功環(huán)小信號控制框圖中的正負(fù)反饋效應(yīng),分析并網(wǎng)逆變器穩(wěn)定性隨無功環(huán)控制器參數(shù)增加而出現(xiàn)先增強后減弱的“拐點”現(xiàn)象。最后,通過仿真與實驗驗證了該文理論分析的正確性。

    弱電網(wǎng) 并網(wǎng)逆變器 PQ功率控制 小信號模型 不對稱正負(fù)反饋效應(yīng) 穩(wěn)定性

    0 引言

    隨著傳統(tǒng)化石能源的大量消耗及其對環(huán)境帶來的不利影響,分布式發(fā)電技術(shù)受到廣泛的關(guān)注和應(yīng)用[1]。而并網(wǎng)逆變器(Grid-Connected Inverter, GCI)作為分布式發(fā)電與電網(wǎng)之間的接口設(shè)備,將在可再生能源輸送到電網(wǎng)中發(fā)揮著重要作用[2-3]。但可再生能源分布不均導(dǎo)致其發(fā)電單元通常部署在偏遠地區(qū),長距離輸電線路的阻抗較大,使電網(wǎng)表現(xiàn)出弱電網(wǎng)特性[4-5]。因此,逆變器接入電網(wǎng)后極易與電網(wǎng)阻抗發(fā)生交互作用而導(dǎo)致并網(wǎng)系統(tǒng)出現(xiàn)寬頻振蕩等穩(wěn)定性問題,嚴(yán)重威脅新型電力系統(tǒng)的安全、可靠運行[6-8]。

    目前,大量研究成果表明弱電網(wǎng)下并網(wǎng)逆變器的外環(huán)控制對并網(wǎng)系統(tǒng)穩(wěn)定性的影響不容忽視[9-11]。文獻[9]指出在弱電網(wǎng)條件下,逆變器功率外環(huán)、電壓外環(huán)與鎖相環(huán)之間的交互作用會更加明顯,不利于GCI的穩(wěn)定性。文獻[10]研究了有功功率和無功功率參考值對GCI阻抗特性的影響,表明在PQ功率控制模式下,給定有功功率超過閾值時,可能誘發(fā)并網(wǎng)系統(tǒng)出現(xiàn)低頻振蕩現(xiàn)象。文獻[11]利用傳遞函數(shù)模型分析了PQ功率控制下GCI的穩(wěn)定性,指出增加有功功率輸出對GCI穩(wěn)定性的負(fù)面影響要強于增加無功功率輸出。此外,文獻[12]比較了不同外環(huán)控制方式下GCI的穩(wěn)定裕度,發(fā)現(xiàn)在定無功功率控制時GCI的穩(wěn)定性最差。文獻[13]重構(gòu)了有功和無功控制環(huán)路,并用d軸電流來表征功率耦合,同時提出一種基于自適應(yīng)無功補償?shù)慕怦畈呗砸蕴嵘到y(tǒng)穩(wěn)定性。但以上研究并未分析有功/無功環(huán)控制器參數(shù)變化對GCI穩(wěn)定性的影響。為此,文獻[14-16]通過構(gòu)建考慮無功功率外環(huán)影響的GCI阻抗模型,分析了無功環(huán)控制器參數(shù)對并網(wǎng)系統(tǒng)穩(wěn)定性的影響規(guī)律,指出并網(wǎng)系統(tǒng)的穩(wěn)定性會隨著無功環(huán)參數(shù)增加而逐漸減弱。然而,上述研究在分析PQ功率控制對GCI運行特性的影響時,并未同時考慮有功/無功環(huán)控制器參數(shù)變化對GCI穩(wěn)定性的影響,難以揭示PQ功率外環(huán)控制器參數(shù)對并網(wǎng)系統(tǒng)穩(wěn)定性影響的真實特性。

    為了充分研究PQ功率外環(huán)控制器參數(shù)對GCI穩(wěn)定性的影響,文獻[17]通過構(gòu)建PQ功率控制下GCI的直流側(cè)阻抗模型,分析了電流內(nèi)環(huán)與PQ外環(huán)控制器參數(shù)對直流側(cè)阻抗特性的影響,并指出電流內(nèi)環(huán)控制帶寬減小會減弱PQ外環(huán)控制器參數(shù)對直流側(cè)阻抗的影響。文獻[18]分析表明,采用PQ功率控制的GCI在單位功率因數(shù)運行下,功率外環(huán)控制器參數(shù)變化對并網(wǎng)系統(tǒng)穩(wěn)定性影響較小,系統(tǒng)是否穩(wěn)定主要取決于有功功率對d軸電流的靈敏度。文獻[19]構(gòu)建了模塊化多電平變換器(Modular Multilevel Converters, MMC)的阻抗模型,并基于該模型指出MMC的阻抗特性會隨著功率外環(huán)控制器比例系數(shù)增加而向較低頻段移動,從而惡化MMC的穩(wěn)定性。文獻[20]基于風(fēng)電變流器分析了功率外環(huán)控制器參數(shù)變化對額定風(fēng)速下并網(wǎng)系統(tǒng)穩(wěn)定性的影響,并提出一種提高風(fēng)電變流器穩(wěn)定性的控制器參數(shù)設(shè)計方法。文獻[21-22]構(gòu)建了考慮PQ功率外環(huán)的風(fēng)電變流器機組簡化阻抗模型,詳細(xì)分析了功率外環(huán)控制器參數(shù)變化對變流器阻抗特性的影響規(guī)律,并指出所建簡化阻抗模型和詳細(xì)阻抗模型的偏差會隨著功率外環(huán)控制帶寬與阻尼比的增加而增大。

    通過以上分析不難發(fā)現(xiàn),現(xiàn)有針對PQ功率控制下GCI穩(wěn)定性的研究依舊沒有深入分析有功/無功環(huán)控制器參數(shù)變化對GCI穩(wěn)定性影響存在的差異性,更沒有從不對稱正負(fù)反饋效應(yīng)的角度揭示導(dǎo)致這種差異性的內(nèi)在機理。

    因此,針對以上研究存在的不足,本文首先推導(dǎo)了PQ功率控制下并網(wǎng)逆變器的小信號數(shù)學(xué)模型,并構(gòu)建了計及電網(wǎng)阻抗影響的并網(wǎng)逆變器的小信號控制框圖;基于控制框圖首次揭示了弱電網(wǎng)下PQ功率控制是一種不對稱控制的本質(zhì),即無功環(huán)控制器不僅會出現(xiàn)在自身控制環(huán)中,而且還會以正負(fù)反饋環(huán)路的形式出現(xiàn)在有功控制環(huán)中。其次,通過阻抗模型分析了有功/無功環(huán)參數(shù)變化對并網(wǎng)逆變器穩(wěn)定性的影響特性,指出無功環(huán)控制器參數(shù)變化對并網(wǎng)系統(tǒng)穩(wěn)定性的影響程度要明顯強于有功環(huán)控制器參數(shù)變化,進一步驗證了弱電網(wǎng)下有功/無功環(huán)小信號控制框圖的不對稱性。在此基礎(chǔ)上,借助小信號控制框圖中形成的正負(fù)反饋環(huán)效應(yīng),分析了并網(wǎng)逆變器穩(wěn)定性隨無功環(huán)控制器參數(shù)增加而出現(xiàn)先增強后減弱的“拐點”現(xiàn)象。最后,通過仿真與實驗驗證了本文理論分析的正確性。

    1 PQ功率控制GCI的小信號數(shù)學(xué)模型構(gòu)建

    為了分析PQ功率控制的不對稱機理及其對GCI穩(wěn)定性的影響特性,本文以三相LCL型并網(wǎng)逆變器為例,構(gòu)建了GCI系統(tǒng)的小信號數(shù)學(xué)模型。圖1給出了PQ功率控制下GCI的系統(tǒng)結(jié)構(gòu),圖中,功率外環(huán)的輸出作為電流內(nèi)環(huán)的參考值。

    圖1 PQ功率控制下LCL型GCI的系統(tǒng)結(jié)構(gòu)

    圖1中,dc為直流側(cè)電壓;逆變器側(cè)濾波電感1、濾波電容f和網(wǎng)側(cè)濾波電感2構(gòu)成LCL濾波器;d為抑制LCL濾波器諧振尖峰的阻尼電阻;1為橋臂輸出電流,為電容支路電壓,2為并網(wǎng)電流,PCC為并網(wǎng)點電壓,g為電網(wǎng)電壓,g為電網(wǎng)阻抗;/、ref/ref分別為GCI有功/無功功率輸出值和參考值;2refd、2refq分別為GCI輸出電流在d軸和q軸的參考值;Md、Mq分別為調(diào)制信號的d、q軸分量;為鎖相環(huán)輸出相位;PWM為并網(wǎng)逆變器等效調(diào)制增益。

    1.1 GCI主電路的小信號數(shù)學(xué)模型構(gòu)建

    在系統(tǒng)dq坐標(biāo)系下,LCL型GCI主電路的小信號數(shù)學(xué)模型[23]為

    GCI系統(tǒng)輸出的有功/無功功率與電壓電流之間的關(guān)系為

    對式(3)進行小信號線性化處理后可得

    其中

    式中,1為無窮大電網(wǎng)的角頻率。

    1.2 GCI控制電路的小信號數(shù)學(xué)模型構(gòu)建

    1)鎖相環(huán)

    為了實現(xiàn)GCI與電網(wǎng)同步,本文采用同步參考坐標(biāo)系鎖相環(huán)(Synchronous Reference Frame Phase Locked Loop, SRF-PLL)結(jié)構(gòu)來實現(xiàn)相位跟蹤。通過對鎖相環(huán)進行小信號建模[24],可以得到表征其輸出特性的小信號數(shù)學(xué)模型為

    式中,p,PLL、i,PLL分別為鎖相環(huán)內(nèi)部PI控制器的比例和積分系數(shù);m0為PCC電壓幅值。

    2)PQ功率外環(huán)

    根據(jù)圖1可以得到外環(huán)PQ功率控制的輸出與內(nèi)環(huán)電流參考值之間的小信號模型為

    式中,上標(biāo)“c”表示控制器dq坐標(biāo)系;P()=pp+ip/為有功環(huán)PI控制器;Q()=qp+qi/為無功環(huán)PI控制器。

    3)電流內(nèi)環(huán)

    由于PLL小信號擾動影響,導(dǎo)致并網(wǎng)逆變器的狀態(tài)變量在系統(tǒng)dq坐標(biāo)系與控制器dq坐標(biāo)系之間存在關(guān)系[24]為

    比等狀態(tài)變量在系統(tǒng)dq坐標(biāo)系和控制器dq坐標(biāo)系的d、q軸分量。

    基于以上關(guān)系,可以得到系統(tǒng)dq坐標(biāo)系下調(diào)制信號Mdq的小信號數(shù)學(xué)模型為

    式中,i()=p,i+i,i/為電流內(nèi)環(huán)PI控制器。

    4)延時環(huán)節(jié)

    考慮PWM與計算延時,通常將其線性化[5]為

    1.3 PQ功率控制的不對稱正負(fù)反饋環(huán)路模型構(gòu)建

    利用控制電路輸出調(diào)制信號與逆變器橋臂輸出信號之間的關(guān)系,可以得到主電路與控制電路之間的數(shù)學(xué)表達式為

    結(jié)合式(4)~式(7)、式(9)~式(11),忽略有功/無功環(huán)功率參考值的擾動,可以得到逆變器橋臂輸出電壓與并網(wǎng)電流和PCC電壓之間的關(guān)系為

    其中

    考慮到在穩(wěn)態(tài)值附近,LCL濾波電感和電容引入的耦合量遠小于系統(tǒng)電壓和電流的穩(wěn)態(tài)值,為了簡化計算分析過程,暫且忽略這些耦合量的影響。結(jié)合式(1)與式(12)可以得到q軸并網(wǎng)電流與PCC電壓d軸、q軸分量之間的表達式為

    其中

    此外,本文忽略等效電網(wǎng)電壓g的小信號擾動。在弱電網(wǎng)下,由于電網(wǎng)阻抗會通過PLL結(jié)構(gòu)在控制環(huán)中形成并網(wǎng)電流正反饋環(huán),嚴(yán)重影響系統(tǒng)的穩(wěn)定性[24]。因此,在構(gòu)建PQ功率控制的小信號模型時,需要將電網(wǎng)阻抗引入到控制環(huán)的等效變換過程中,并以此來反映弱電網(wǎng)對GCI穩(wěn)定性的影響。基于此,聯(lián)立式(12)、式(13)可得d軸電壓與dq軸電流以及q軸電壓與dq軸電流之間的表達式分別為

    根據(jù)式(15),并結(jié)合式(4)、式(5),得到有功功率與d軸電流、無功功率與q軸電流之間的關(guān)系為

    其中

    基于以上分析,可得PQ功率控制下GCI有功環(huán)與無功環(huán)的小信號控制框圖,如圖2所示。由于電網(wǎng)阻抗的存在,PLL會在電流參考值和調(diào)制信號處誘發(fā)產(chǎn)生相應(yīng)的正負(fù)反饋環(huán)路[24],等效到圖2中即形成負(fù)反饋環(huán)路1中的asy,q()與正反饋環(huán)路4中的asy,d()。此外,當(dāng)并網(wǎng)逆變器以單位功率因數(shù)運行時,并網(wǎng)電流的q軸分量為0,式(15)可進一步簡化為

    圖2 PQ功率控制中有功環(huán)路和無功環(huán)路的等效小信號控制框圖

    與此同時,根據(jù)圖2可以看出,無功環(huán)控制器參數(shù)Q()不僅存在于自身控制環(huán)路中,而且還會以正負(fù)反饋環(huán)的形式存在于有功環(huán)路中,這意味著無功環(huán)控制器參數(shù)將同時影響有功環(huán)和無功環(huán)的穩(wěn)定性。然而,有功環(huán)控制器參數(shù)卻只出現(xiàn)在有功環(huán)中,僅對自身環(huán)路穩(wěn)定性產(chǎn)生影響。因此,在考慮電網(wǎng)阻抗影響時,PQ功率控制并網(wǎng)逆變器的小信號控制框圖表現(xiàn)出一種不對稱控制現(xiàn)象,揭示了PQ功率控制是一種不對稱控制的數(shù)學(xué)本質(zhì)。

    值得說明的是,上述推導(dǎo)的PQ功率控制下并網(wǎng)逆變器的小信號控制框圖是在考慮電網(wǎng)阻抗影響與PLL小信號擾動的情形下構(gòu)建的,圖2中所形成的正負(fù)反饋環(huán)路都與電網(wǎng)阻抗或PLL傳遞函數(shù)直接相關(guān)。若不考慮電網(wǎng)阻抗影響,則圖2中的正負(fù)反饋環(huán)路就不會被形成。因此,本文所推導(dǎo)的分析模型僅適用于新能源裝機容量不斷增加下考慮電網(wǎng)阻抗影響的并網(wǎng)逆變器穩(wěn)定性分析場景。

    2 不對稱正負(fù)反饋效應(yīng)對GCI穩(wěn)定性的影響分析

    2.1 有功/無功環(huán)控制器參數(shù)變化對GCI穩(wěn)定性影響的差異性分析

    通過第1節(jié)的分析不難發(fā)現(xiàn),PQ功率控制GCI小信號控制框圖中的正負(fù)反饋環(huán)路都只與無功環(huán)控制器有關(guān)。這意味著相較于有功環(huán)控制器參數(shù)而言,GCI的穩(wěn)定性對無功環(huán)參數(shù)變化會表現(xiàn)的更加敏感。為了驗證上述猜想,本文通過構(gòu)建PQ功率控制下GCI的阻抗模型,并利用基于阻抗的穩(wěn)定性判據(jù)分析有功/無功環(huán)參數(shù)變化對GCI穩(wěn)定性影響程度的差異性。圖3給出了矩陣形式下PQ功率控制GCI的小信號控制框圖模型。

    基于圖3可以求得GCI的等效輸出阻抗模型為

    其中

    需要說明的是,式(18)中LCL濾波環(huán)節(jié)、延時環(huán)節(jié)、內(nèi)環(huán)控制器及asy1()、asy2()的具體表達式詳見文獻[24]。同時,電網(wǎng)阻抗的矩陣表示為

    根據(jù)基于阻抗的穩(wěn)定性判據(jù)可知,弱電網(wǎng)下GCI的穩(wěn)定性可以通過回率矩陣()的特征值是否包圍點(-1, j0)來進行判定[24-25]。其中,GCI系統(tǒng)回率矩陣()的表達式為

    本文所采用LCL型GCI系統(tǒng)的主要參數(shù)見表1,并分四種情形來分析有功/無功環(huán)參數(shù)變化對并網(wǎng)系統(tǒng)穩(wěn)定性的影響。其中,PQ功率環(huán)控制器參數(shù)見表2。

    表1 系統(tǒng)主要參數(shù)

    表2 有功與無功環(huán)PI控制器參數(shù)

    在表1與表2參數(shù)下,圖4a給出了有功控制環(huán)比例系數(shù)pp變化時,即情形1的穩(wěn)定性分析結(jié)果。觀察圖4a可以看出:當(dāng)pp從0.000 6增大到0.005 2時,回率矩陣()的特征值1()的Nyquist曲線均不包圍點(-1, j0),此時GCI系統(tǒng)處于穩(wěn)定運行狀態(tài)。圖4b給出了有功控制環(huán)積分系數(shù)pi變化時,即情形2的穩(wěn)定性分析結(jié)果。觀察圖4b可以看出:當(dāng)pi從0.5增大到5時,回率矩陣()的特征值1()的Nyquist曲線都不包圍點(-1, j0),此時并網(wǎng)逆變器系統(tǒng)同樣處于穩(wěn)定運行狀態(tài)。值得說明的是:當(dāng)pi從0.5增大到1,再增大到1.6時,特征值1()與實軸的交點先遠離再靠近點(-1, j0);當(dāng)pi從1.6增大到5時,特征值1()與實軸的交點明顯靠近點(-1, j0)。這意味著,隨著pi的增大,GCI系統(tǒng)的穩(wěn)定裕度先小幅度提升,而后再降低。

    圖4 情形1與情形2下GCI的穩(wěn)定性分析結(jié)果

    為了驗證有功環(huán)參數(shù)變化對GCI系統(tǒng)穩(wěn)定性影響分析的正確性,圖5a與圖5b分別給出了pp與pi變化時GCI的輸出電流仿真結(jié)果。

    觀察圖5a可以看出:當(dāng)有功環(huán)比例系數(shù)pp從0.000 6增大到0.005 2時,GCI輸出電流波形都能夠保持穩(wěn)定運行。而且,通過對三組參數(shù)下的并網(wǎng)電流進行快速傅里葉變換(Fast Fourier Transform, FFT)分析,得到相應(yīng)的總諧波畸變率(Total Harmonic Distortion, THD)分別為2.21%、2.37%和3.80%。FFT分析結(jié)果表明,隨著pp的增大,并網(wǎng)電流的諧波含量也隨之增加,即GCI系統(tǒng)的穩(wěn)定裕度逐漸降低。以上分析驗證了圖4a理論分析的正確性。

    圖5 有功環(huán)控制器參數(shù)變化時的并網(wǎng)電流仿真結(jié)果

    觀察圖5b可以看出:當(dāng)有功環(huán)積分系數(shù)pi從0.5增大到5時,GCI輸出電流同樣能夠保持穩(wěn)定運行。同時,通過對pi=0.5、1.0、1.6、5.0時的并網(wǎng)電流進行FFT分析,可得相應(yīng)的THD分別為2.58%、1.85%、2.66%和4.42%。FFT分析結(jié)果表明:隨著pi增大,并網(wǎng)電流的諧波含量表現(xiàn)出先減小后增大的特性,這意味著GCI系統(tǒng)的穩(wěn)定裕度先增加后減小,與圖4b的理論分析結(jié)果一致。

    同理,圖6給出了無功環(huán)參數(shù)變化時GCI的穩(wěn)定性分析結(jié)果。

    觀察圖6a可以看出:當(dāng)比例系數(shù)qp=0.001 5時,回率矩陣()的特征值曲線都不包圍點(-1, j0),此時并網(wǎng)系統(tǒng)處于穩(wěn)定運行狀態(tài);而當(dāng)qp繼續(xù)增加到0.005 2時,特征值1()的Nyquist曲線包圍了點(-1, j0),此時系統(tǒng)失穩(wěn)。

    觀察圖6b可以看出,當(dāng)積分系數(shù)qi從0.5增大到1時,回率矩陣()的特征值曲線1()都不包圍點(-1, j0),此時系統(tǒng)能夠保持穩(wěn)定運行;而當(dāng)qi增大到1.6時,特征值1()的Nyquist曲線包圍點(-1, j0),此時GCI處于失穩(wěn)狀態(tài)。

    圖6 情形3與情形4下GCI的穩(wěn)定性分析結(jié)果

    為了驗證圖6理論分析的正確性,圖7給出了相應(yīng)情形下的并網(wǎng)電流仿真結(jié)果。

    圖7 無功環(huán)控制器參數(shù)變化時的并網(wǎng)電流仿真結(jié)果

    觀察圖7a可以看出:當(dāng)比例系數(shù)qp=0.001 5時,GCI輸出電流波形穩(wěn)定;而當(dāng)qp增加到0.005 2時,輸出電流波形發(fā)生明顯畸變,表現(xiàn)出諧波失穩(wěn)現(xiàn)象。同時,觀察圖7b可以看出:當(dāng)qi增大到1.6時,GCI輸出電流已經(jīng)明顯失穩(wěn)。上述仿真結(jié)果驗證了圖6理論分析的正確性。

    對比有功環(huán)和無功環(huán)參數(shù)變化對GCI系統(tǒng)穩(wěn)定性的影響分析可以發(fā)現(xiàn):當(dāng)有功環(huán)比例系數(shù)pp從0.000 6增大到0.005 2時,GCI依舊能夠保持穩(wěn)定運行;然而無功環(huán)比例系數(shù)qp增加到0.005 2,系統(tǒng)已經(jīng)失穩(wěn)。除此之外,當(dāng)有功環(huán)積分系數(shù)pi從0.5增大到5時,GCI同樣能夠保持穩(wěn)定運行;然而無功環(huán)積分系數(shù)qi增大到1.6時,系統(tǒng)就已經(jīng)表現(xiàn)出失穩(wěn)。以上對比分析表明,GCI系統(tǒng)的穩(wěn)定性對無功環(huán)參數(shù)變化更加敏感,換句話說就是無功環(huán)參數(shù)變化對系統(tǒng)穩(wěn)定性的影響程度要明顯強于有功環(huán)參數(shù),這也證明了弱電網(wǎng)下PQ功率控制在本質(zhì)上是一種不對稱控制。

    2.2 不對稱正負(fù)反饋效應(yīng)對GCI穩(wěn)定性影響的特性分析

    由以上分析可以看出,無功環(huán)控制器Q()對GCI的穩(wěn)定性將產(chǎn)生“雙重影響”。當(dāng)調(diào)整無功環(huán)控制器參數(shù)使“環(huán)路1”與“環(huán)路4”的負(fù)反饋效應(yīng)強于“環(huán)路2”與“環(huán)路3”的正反饋效應(yīng)時,GCI的穩(wěn)定性會增強;當(dāng)調(diào)整無功環(huán)控制器參數(shù)使“環(huán)路2”與“環(huán)路3”的正反饋效應(yīng)強于“環(huán)路1”與“環(huán)路4”的負(fù)反饋效應(yīng)時,GCI的穩(wěn)定性會減弱。

    為驗證由無功環(huán)控制器參數(shù)引入的不對稱正負(fù)反饋環(huán)路對GCI穩(wěn)定性影響特性分析的正確性,圖8給出了無功環(huán)控制器參數(shù)qp變化時的穩(wěn)定性分析結(jié)果。

    觀察圖8a可以看出:當(dāng)無功環(huán)參數(shù)qp從0.000 6增大到0.002 25時,回率矩陣()的特征值1()與實軸的交點向坐標(biāo)原點移動,即逐漸由包圍點(-1, j0)向不包圍點(-1, j0)移動,并網(wǎng)系統(tǒng)的穩(wěn)定裕度增加。這表明隨著qp的增加,由無功環(huán)控制器參數(shù)所形成的負(fù)反饋環(huán)效應(yīng)要強于正反饋環(huán)效應(yīng),GCI系統(tǒng)的穩(wěn)定性增強。當(dāng)無功環(huán)參數(shù)qp從0.002 25繼續(xù)增加到0.005 2時,回率矩陣()的特征值1()與實軸的交點向遠離坐標(biāo)原點的方向移動,即逐漸由不包圍點(-1, j0)向包圍點(-1, j0)移動,如圖8b所示。這意味著此時由無功環(huán)控制器參數(shù)形成的正反饋環(huán)效應(yīng)要強于負(fù)反饋環(huán)效應(yīng),GCI系統(tǒng)的穩(wěn)定性減弱。

    圖8 無功環(huán)控制器參數(shù)變化時的穩(wěn)定性分析結(jié)果

    為了更加直觀地展示無功環(huán)控制器參數(shù)對并網(wǎng)逆變器穩(wěn)定性的影響規(guī)律,圖9給出了用GCI輸出阻抗invqq()與電網(wǎng)阻抗gqq()表征的GCI穩(wěn)定裕度隨參數(shù)qp變化的Bode圖。

    觀察圖9可以看出,隨著無功環(huán)控制器參數(shù)qp的增大,GCI表現(xiàn)出由不穩(wěn)定向穩(wěn)定,再向不穩(wěn)定變化的規(guī)律,即GCI的穩(wěn)定性表現(xiàn)出一種“拐點”行為,而非現(xiàn)有研究中表明的系統(tǒng)穩(wěn)定性隨無功環(huán)參數(shù)增加而逐漸減弱[14-16]。

    圖9 無功環(huán)控制器參數(shù)變化時的Bode圖

    為了驗證上述理論分析的正確性,圖10給出了無功環(huán)控制器參數(shù)qp變化時的并網(wǎng)電流仿真結(jié)果。

    觀察圖10可以看出:當(dāng)qp=0.000 6時,GCI輸出電流發(fā)生明顯畸變,對輸出電流進行FFT分析可得,此時并網(wǎng)電流的THD=67.84%,系統(tǒng)失穩(wěn)。而當(dāng)qp增加到0.002 25時,GCI輸出電流波形穩(wěn)定。當(dāng)qp進一步增加到0.005 2時,與qp=0.000 6類似,GCI系統(tǒng)同樣出現(xiàn)明顯的諧波失穩(wěn)現(xiàn)象。以上仿真結(jié)果驗證了本文理論分析的正確性。

    圖10 無功環(huán)控制器參數(shù)變化時的并網(wǎng)電流仿真結(jié)果

    3 實驗驗證

    為驗證上述理論分析的正確性,在實驗室搭建了一套三相LCL型并網(wǎng)逆變器實驗平臺,其中平臺實驗參數(shù)與表1一致。圖11和圖12分別給出了有功/無功環(huán)參數(shù)變化時的并網(wǎng)實驗波形。

    觀察圖11a可以看出,在初始參數(shù)條件下,并網(wǎng)逆變器能夠輸出穩(wěn)定的并網(wǎng)電流波形,即系統(tǒng)穩(wěn)定。在圖11a初始參數(shù)基礎(chǔ)上,當(dāng)調(diào)整有功環(huán)比例系數(shù)pp從0.001 5變化到0.005 2時,并網(wǎng)電流依舊穩(wěn)定,如圖11b所示;而當(dāng)無功環(huán)比例系數(shù)qp從0.001 5變化到0.005 2時,并網(wǎng)電流發(fā)生明顯畸變,如圖11c所示。同樣觀察圖12可以看出,在初始參數(shù)條件下,當(dāng)調(diào)整有功環(huán)積分系數(shù)pi從1變化到1.6時,并網(wǎng)電流依舊穩(wěn)定,如圖12a所示;而當(dāng)無功環(huán)積分系數(shù)qi從1變化到1.6時,并網(wǎng)電流發(fā)生明顯畸變,如圖12b所示。通過以上實驗結(jié)果對比分析表明,無功環(huán)參數(shù)變化對系統(tǒng)穩(wěn)定性的影響程度要明顯強于有功環(huán)參數(shù),這與圖4、圖6理論分析結(jié)果相一致。

    圖11 有功/無功環(huán)比例系數(shù)變化時的實驗波形

    圖12 有功/無功環(huán)積分系數(shù)變化時的實驗波形

    為驗證無功環(huán)參數(shù)形成的正負(fù)反饋環(huán)效應(yīng)對GCI穩(wěn)定性影響特性分析的正確性,圖13給出了無功環(huán)參數(shù)分別為0.000 6與0.002 25時的實驗結(jié)果。

    圖13 無功環(huán)比例系數(shù)變化時的實驗波形

    結(jié)合圖11c和圖13不難發(fā)現(xiàn):當(dāng)無功環(huán)參數(shù)qp從0.000 6增加到0.005 2時,GCI會出現(xiàn)由最初的不穩(wěn)定運行過渡到穩(wěn)定運行,而后隨著參數(shù)增加又再次失穩(wěn)的一個變化過程,即GCI的穩(wěn)定性隨無功環(huán)參數(shù)增加會出現(xiàn)“拐點”現(xiàn)象。以上實驗結(jié)果進一步驗證了圖8與圖9理論分析的正確性。

    4 結(jié)論

    本文通過構(gòu)建PQ功率控制GCI的小信號數(shù)學(xué)模型,揭示了弱電網(wǎng)下PQ功率控制是一種不對稱

    控制的數(shù)學(xué)本質(zhì);同時利用有功與無功環(huán)路中的正負(fù)反饋環(huán)效應(yīng),揭示了有功/無功環(huán)控制器參數(shù)變化對GCI穩(wěn)定性影響存在差異性的機理,并得出以下重要結(jié)論:

    1)揭示了弱電網(wǎng)下PQ功率控制是一種不對稱控制的本質(zhì),即無功環(huán)控制器會在有功和無功控制環(huán)路中形成不對稱的正負(fù)反饋環(huán)路。

    2)無功環(huán)控制器參數(shù)變化對GCI穩(wěn)定性的影響程度要強于有功環(huán)控制器參數(shù)變化,即GCI的穩(wěn)定性對無功環(huán)參數(shù)變化更加敏感。

    3)調(diào)整無功環(huán)控制器參數(shù)時,不對稱正負(fù)反饋環(huán)效應(yīng)會出現(xiàn)彼此抵消的現(xiàn)象,導(dǎo)致系統(tǒng)穩(wěn)定性出現(xiàn)不穩(wěn)定到穩(wěn)定再到不穩(wěn)定的一個變化過程,即“拐點”現(xiàn)象,并非是并網(wǎng)系統(tǒng)穩(wěn)定性隨無功環(huán)控制器參數(shù)增加而逐漸減弱。

    值得說明的是,如何通過構(gòu)建解析表達式或者量化指標(biāo)的方式尋求正、負(fù)反饋環(huán)效應(yīng)相互抵消的臨界點,將是未來需要進一步探索的工作。

    [1] 張學(xué)廣, 夏丹妮, 陳文佳, 等. 三相并網(wǎng)變流器弱電網(wǎng)下頻率耦合抑制控制方法[J]. 電工技術(shù)學(xué)報, 2019, 34(21): 4559-4571.

    Zhang Xueguang, Xia Danni, Chen Wenjia, et al. Frequency coupling suppression control method for three-phase grid-connected converter under weak grid[J]. Transactions of China Electrotechnical Society, 2019, 34(21): 4559-4571.

    [2] 陳燕東, 王伊, 周樂明, 等. 弱電網(wǎng)下LCL逆變器阻尼諧振抑制與功率快速調(diào)節(jié)方法[J]. 電工技術(shù)學(xué)報, 2018, 33(11): 2564-2574.

    Chen Yandong, Wang Yi, Zhou Leming, et al. Damping resonance suppression and fast power regulation method for LCL-type inverter under weak grid[J]. Transactions of China Electrotechnical Society, 2018, 33(11): 2564-2574.

    [3] Zhu Donghai, Zhou Shiying, Zou Xudong, et al. Improved design of PLL controller for LCL-type grid-connected converter in weak grid[J]. IEEE Transactions on Power Electronics, 2020, 35(5): 4715-4727.

    [4] 張旸, 陳新, 王昀, 等. 弱電網(wǎng)下并網(wǎng)逆變器的阻抗相角動態(tài)控制方法[J]. 電工技術(shù)學(xué)報, 2017, 32(1): 97-106.

    Zhang Yang, Chen Xin, Wang Yun, et al. Impedance phased dynamic control method of grid-connected inverters under weak grid condition[J]. Transactions of China Electrotechnical Society, 2017, 32(1): 97- 106.

    [5] Zhang Xueguang, Xia Danni, Fu Zhichao, et al. An improved feedforward control method considering PLL dynamics to improve weak grid stability of grid-connected inverters[J]. IEEE Transactions on Industry Applications, 2018, 54(5): 5143-5151.

    [6] 涂春鳴, 高家元, 趙晉斌, 等. 弱電網(wǎng)下具有定穩(wěn)定裕度的并網(wǎng)逆變器阻抗重塑分析與設(shè)計[J]. 電工技術(shù)學(xué)報, 2020, 35(6): 1327-1335.

    Tu Chunming, Gao Jiayuan, Zhao Jinbin, et al. Analysis and design of grid-connected inverter impedance remodeling with fixed stability margin in weak grid[J]. Transactions of China Electrotechnical Society, 2020, 35(6): 1327-1335.

    [7] 李龍源, 付瑞清, 呂曉琴, 等. 接入弱電網(wǎng)的同型機直驅(qū)風(fēng)電場單機等值建模[J/OL]. 電工技術(shù)學(xué)報, 2022, DOI: 10.19595/j.cnki.1000-6753.tces.211555.

    Li Longyuan, Fu Ruiqing, Lü Xiaoqin, et al. Single machine equivalent modeling of weak grid connected wind farm with same type PMSGs[J/OL]. Transa- ctions of China Electrotechnical Society: 2022, DOI: 10.19595/j.cnki.1000-6753.tces.211555.

    [8] 張占俊, 李建文, 董耀, 等. 弱電網(wǎng)下多逆變器并網(wǎng)諧振失穩(wěn)分析方法[J]. 電氣技術(shù), 2020, 21(10): 21-28.

    Zhang Zhanjun, Li Jianwen, Dong Yao, et al. Method of resonance instability analysis of multiple grid- connected inverters in weak grid[J]. Electrical Engineering, 2020, 21(10): 21-28.

    [9] Arani M F M, Mohamed Y A I. Analysis and per- formance enhancement of vector-controlled VSC in HVDC links connected to very weak grids[J]. IEEE Transactions on Power Systems, 2017, 32(1): 684- 693.

    [10] Zhou Weihua, Wang Yanbo, Torres-Olguin R E, et al. Effect of reactive power characteristic of offshore wind power plant on low-frequency stability[J]. IEEE Transactions on Energy Conversion, 2020, 35(2): 837-853.

    [11] Wu Guanglu, Sun Huadong, Zhang Xi, et al. Parameter design oriented analysis of the current control stability of the weak-grid-tied VSC[J]. IEEE Transactions on Power Delivery, 2021, 36(3): 1458- 1470.

    [12] Huang Linbin, Xin Huanhai, Li Zhiyi, et al. Grid- synchronization stability analysis and loop shaping for PLL-based power converters with different reactive power control[J]. IEEE Transactions on Smart Grid, 2020, 11(1): 501-516.

    [13] 姜靜雅, 王瑋, 吳學(xué)智, 等. 基于自適應(yīng)無功功率補償?shù)奶摂M同步機功率解耦策略[J]. 電工技術(shù)學(xué)報, 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.

    [14] 徐路遙, 蘭洲, 陳飛, 等. 弱電網(wǎng)下無功控制對并網(wǎng)變流器穩(wěn)定性影響分析[J]. 電力系統(tǒng)自動化, 2019, 43(15): 56-63.

    Xu Luyao, Lan Zhou, Chen Fei, et al. Analysis on impact of reactive power control on stability of grid-connected converter in weak grid[J]. Automation of Electric Power Systems, 2019, 43(15): 56-63.

    [15] Xu Luyao, Yuan Hui, Lan Zhou, et al. Impedance- based analysis of constant reactive power control for voltage source converter in weak AC grids[C]//2019 IEEE Power & Energy Society General Meeting (PESGM), Atlanta, GA, USA, 2019: 1-5.

    [16] 任必興, 孫蓉, 李強, 等. 弱連接下并網(wǎng)VSC外環(huán)控制穩(wěn)定性分析[J]. 高電壓技術(shù), 2022, 48(2): 762- 771.

    Ren Bixing, Sun Rong, Li Qiang, et al. Stability analysis of power control outer loop of grid- connected VSC under weak connection conditions[J]. High Voltage Engineering, 2022, 48(2): 762-771.

    [17] Xue Danhong, Liu Jinjun, Liu Zeng, et al. Modeling and analysis of DC terminal impedance of voltage- source converters with different control modes[J]. IEEE Transactions on Power Electronics, 2020, 35(6): 5883-5896.

    [18] 邵冰冰, 趙書強, 高本鋒, 等. 連接弱交流電網(wǎng)的VSC-HVDC失穩(wěn)機理及判據(jù)研究[J]. 電工技術(shù)學(xué)報, 2019, 34(18): 3884-3889.

    Shao Bingbing, Zhao Shuqiang, Gao Benfeng, et al. Instability mechanism and criterion analysis of VSC- HVDC connected to the weak AC power grid[J]. Transactions of China Electrotechnical Society, 2019, 34(18): 3884-3889.

    [19] Yin Taiyuan, Wang Yue, Wang Xiaolei, et al. Mode- ling and analysis of high-frequency MMC impedance considering different control modes and voltage feed- forward[J]. IEEE Access, 2020, 8: 218575- 218584.

    [20] Jia Yanbing, Huang Tao, Li Yubo, et al. Parameter setting strategy for the controller of the DFIG wind turbine considering the small-signal stability of power grids[J]. IEEE Access, 2020, 8: 31287-31294.

    [21] 李光輝, 王偉勝, 張興, 等. 雙饋風(fēng)電場并網(wǎng)次/超同步振蕩建模與機理分析(一): 考慮功率外環(huán)的阻抗建模[J]. 中國電機工程學(xué)報, 2022, 42(7): 2438- 2448.

    Li Guanghui, Wang Weisheng, Zhang Xing, et al. Modeling and suppression method of sub/super- synchronous oscillation of grid-connected DFIG wind farms (part Ⅰ): sequence impedance modeling and analysis of DFIG wind turbines considering power loop[J]. Proceedings of the CSEE, 2022, 42(7): 2438-2448.

    [22] 李光輝, 王偉勝, 張興, 等. 雙饋風(fēng)電場并網(wǎng)次/超同步振蕩建模與機理分析(二): 阻抗特性與振蕩機理分析[J]. 中國電機工程學(xué)報, 2022, 42(10): 3614- 3627.

    Li Guanghui, Wang Weisheng, Zhang Xing, et al. Modeling and mechanism analysis of sub/super- synchronous oscillation of grid-connected DFIG wind farms (part Ⅱ): analysis of impedance characteristic and oscillation mechanism[J]. Proceedings of the CSEE, 2022, 42(10): 3614-3627.

    [23] 張琦, 李銳, 張永平, 等. 具有主動阻尼特征的LCL并網(wǎng)逆變器狀態(tài)反饋控制方法[J]. 電工技術(shù)學(xué)報, 2019, 34(8): 1698-1707.

    Zhang Qi, Li Rui, Zhang Yongping, et al. Active damping method based on state-feedback for grid- connected inverter with LCL filter[J]. Transactions of China Electrotechnical Society, 2019, 34(8): 1698- 1707.

    [24] 高家元, 涂春鳴, 肖凡, 等. 弱電網(wǎng)下PLL結(jié)構(gòu)引入的不對稱正反饋環(huán)路對并網(wǎng)逆變器穩(wěn)定性的影響分析[J]. 中國電機工程學(xué)報, 2022, 42(3): 1113- 1124.

    Gao Jiayuan, Tu Chunming, Xiao Fan, et al. Analysis of the influence of the asymmetric positive-feedback loop introduced by the PLL structure on the stability of grid-connected inverter in weak grids[J]. Pro- ceedings of the CSEE, 2022, 42(3): 1113-1124.

    [25] 劉一鋒, 周小平, 洪樂榮, 等. 虛擬慣性控制的負(fù)荷變換器接入弱電網(wǎng)的序阻抗建模與穩(wěn)定性分析[J]. 電工技術(shù)學(xué)報, 2021, 36(4): 843-856.

    Liu Yifeng, Zhou Xiaoping, Hong Lerong, et al. Sequence impedance modeling and stability analysis of load converter with virtual inertia control con- nected to weak grid[J]. Transactions of China Elec- trotechnical Society, 2021, 36(4): 843-856.

    Stability Analysis of Grid-Connected Inverter under PQ Power Control Based on Asymmetric Positive-Negative-Feedback Effects

    (National Electric Power Conversion and Control Engineering Technology Research Center Hunan University Changsha 410082 China)

    With the large consumption of traditional fossil energy and its adverse impact on the environment, distributed power generation technology has been widely concerned and applied. The stability problem brought by PQ power control to grid-connected inverter (GCI) system under a weak grid has been widely studied.

    Firstly, the small signal mathematical model of grid connected inverter under PQ power control is derived, and the small signal control block diagram of grid connected inverter considering the influence of grid impedance is constructed. In the small signal control block diagram, there are two positive feedback loops and two negative feedback loops at the same time. At the same time, from the composition of positive feedback and negative feedback loops, it can be seen that the parameters of the reactive power loop controller not only exist in the self control loop, but also exist in the active power loop in the form of positive and negative feedback loops, which means that the parameters of the reactive power loop controller will affect the stability of the active power loop and the reactive power loop at the same time. However, the parameters of the active power loop controller only appear in the active power loop, which only affects the stability of its own loop. Therefore, when considering the influence of grid impedance, the small signal control block diagram of PQ power control grid connected inverter shows an asymmetric control phenomenon, revealing the mathematical essence of PQ power control as an asymmetric control. Secondly, in order to verify the influence of active/reactive loop parameters on the stability of grid connected inverter, the impedance model of PQ power controlled grid connected inverter under weak current grid is built. It can be seen from the simulation and experimental results that, within the same parameter variation range, the influence of the parameter variation of the reactive power loop controller on the stability of the grid connected system is obviously stronger than the parameter variation of the active power loop controller, which further verifies the asymmetry of the small signal control block diagram of the active/reactive power loop under the weak current network. On this basis, with the help of the positive and negative feedback loop effect formed in the small signal control block diagram, the law of the stability of grid connected inverter changing with the parameters of the reactive power loop controller is analyzed. When the parameters of the reactive power loop controller are adjusted, the asymmetric positive and negative feedback loop effects will cancel each other, leading to a change process of system stability from instability to stability and then to instability. This is not a gradual weakening of the stability of the grid connected system as the parameters of the reactive power loop controller increase.

    The following conclusions can be drawn from the simulation analysis: ① It is revealed that PQ power control in weak current network is the essence of asymmetric control, that is, the reactive power loop controller will form an asymmetric positive and negative feedback loop in the active and reactive power control loops. ② The change of reactive power loop controller parameters has a stronger influence on the stability of GCI than the change of active power loop controller parameters, that is, the stability of GCI is more sensitive to the change of reactive power loop parameters. ③ When the parameters of the reactive power loop controller are adjusted, the asymmetric positive and negative feedback loop effects will cancel each other, leading to a change process of system stability from instability to stability and then to instability, that is, the “inflection point” phenomenon. It is not that the stability of the grid connected system gradually weakens with the increase of the parameters of the reactive power loop controller. It is worth noting that how to find the critical point of mutual cancellation of positive and negative feedback loop effects by constructing analytical expressions or quantitative indicators will be the work that needs further exploration in the future.

    Weak grid, grid-connected inverter, PQ power control, small-signal model, asymmetric positive- negative-feedback effects, stability

    TM46

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

    國家自然科學(xué)基金項目(52077072, 51907057)、湖南省科技創(chuàng)新領(lǐng)軍人才項目(2019RS3014)資助。

    2022-04-24

    2022-05-30

    涂春鳴 男,1976年生,教授,博士生導(dǎo)師,研究方向為電力電子技術(shù)及其在電力系統(tǒng)中的應(yīng)用。

    E-mail: chunming_tu@263.net

    高家元 男,1991年生,博士研究生,研究方向為電力電子技術(shù),分布式發(fā)電和逆變器控制技術(shù)。

    E-mail: Jiayuan_gao@163.com(通信作者)

    (編輯 陳 誠)

    猜你喜歡
    負(fù)反饋環(huán)路穩(wěn)定性
    全新的虛短虛斷概念與兩類集成運放之導(dǎo)出
    負(fù)反饋放大電路設(shè)計
    電子制作(2019年23期)2019-02-23 13:21:36
    非線性中立型變延遲微分方程的長時間穩(wěn)定性
    上海市中環(huán)路標(biāo)線調(diào)整研究
    上海公路(2018年4期)2018-03-21 05:57:46
    半動力系統(tǒng)中閉集的穩(wěn)定性和極限集映射的連續(xù)性
    基于Multisim的負(fù)反饋放大電路仿真分析
    Buck-Boost變換器的環(huán)路補償及仿真
    電測與儀表(2014年8期)2014-04-04 09:19:36
    單脈沖雷達導(dǎo)引頭角度跟蹤環(huán)路半實物仿真
    莫斯科地鐵計劃于2019—2020年推出第三換乘環(huán)路
    模糊微分方程的一致穩(wěn)定性
    精品国产乱码久久久久久小说| 人人妻人人爽人人添夜夜欢视频 | 午夜精品国产一区二区电影| 国产精品爽爽va在线观看网站| 一区在线观看完整版| 交换朋友夫妻互换小说| 黄色配什么色好看| 在线观看免费视频网站a站| 国产乱人偷精品视频| 免费看日本二区| 日韩一区二区三区影片| 国产精品久久久久久精品电影小说 | 免费观看性生交大片5| 欧美性感艳星| 爱豆传媒免费全集在线观看| 亚洲国产精品专区欧美| 一区在线观看完整版| 国产一区二区三区av在线| 22中文网久久字幕| 国产成人精品久久久久久| 色综合色国产| 国产欧美亚洲国产| 熟女av电影| 国产免费又黄又爽又色| av黄色大香蕉| 久久久久久久国产电影| 蜜桃亚洲精品一区二区三区| 精品人妻熟女av久视频| 男人和女人高潮做爰伦理| 18禁裸乳无遮挡免费网站照片| 亚洲婷婷狠狠爱综合网| av卡一久久| 卡戴珊不雅视频在线播放| 天堂8中文在线网| 国产黄片视频在线免费观看| 日本av手机在线免费观看| av一本久久久久| 免费av中文字幕在线| 国产国拍精品亚洲av在线观看| 色婷婷av一区二区三区视频| 精品久久久久久电影网| 在线观看一区二区三区| 成人黄色视频免费在线看| 亚洲精品日韩av片在线观看| 日产精品乱码卡一卡2卡三| 亚洲精品视频女| 男人舔奶头视频| 国产极品天堂在线| 亚洲精品自拍成人| 22中文网久久字幕| 交换朋友夫妻互换小说| 天天躁夜夜躁狠狠久久av| 国产高潮美女av| 身体一侧抽搐| 国产成人精品久久久久久| 国产精品一区二区性色av| av.在线天堂| 久久人妻熟女aⅴ| 国产成人精品久久久久久| 成年美女黄网站色视频大全免费 | 亚洲婷婷狠狠爱综合网| 欧美成人a在线观看| 亚洲成人手机| 王馨瑶露胸无遮挡在线观看| av国产久精品久网站免费入址| 插阴视频在线观看视频| 插阴视频在线观看视频| 亚洲精品日本国产第一区| 日韩 亚洲 欧美在线| 水蜜桃什么品种好| 少妇人妻精品综合一区二区| av不卡在线播放| 大香蕉97超碰在线| 成年人午夜在线观看视频| 成年人午夜在线观看视频| 人妻少妇偷人精品九色| 人妻少妇偷人精品九色| 午夜激情福利司机影院| 51国产日韩欧美| 欧美xxxx性猛交bbbb| 国产黄色视频一区二区在线观看| 少妇 在线观看| 免费黄色在线免费观看| 免费在线观看成人毛片| 性色av一级| 亚洲国产成人一精品久久久| 女的被弄到高潮叫床怎么办| 国产爽快片一区二区三区| 99热这里只有精品一区| 丝瓜视频免费看黄片| 国产在线一区二区三区精| 内地一区二区视频在线| 亚洲国产成人一精品久久久| 精品亚洲乱码少妇综合久久| 一个人免费看片子| 国产精品麻豆人妻色哟哟久久| 国产成人精品一,二区| 日韩电影二区| 国产伦理片在线播放av一区| 欧美成人精品欧美一级黄| 国产精品99久久99久久久不卡 | 51国产日韩欧美| a级毛色黄片| 国产精品久久久久久精品电影小说 | 久久精品国产鲁丝片午夜精品| 乱码一卡2卡4卡精品| 观看免费一级毛片| 大又大粗又爽又黄少妇毛片口| a级一级毛片免费在线观看| 亚洲欧美成人精品一区二区| 精品少妇黑人巨大在线播放| 国产欧美亚洲国产| 免费观看在线日韩| 91午夜精品亚洲一区二区三区| 亚洲真实伦在线观看| 久久精品国产自在天天线| 精品久久久久久电影网| 国产国拍精品亚洲av在线观看| 中国国产av一级| 卡戴珊不雅视频在线播放| 亚洲国产精品成人久久小说| 一级黄片播放器| 热re99久久精品国产66热6| 熟妇人妻不卡中文字幕| 国产日韩欧美亚洲二区| av专区在线播放| 国产伦在线观看视频一区| 国产在线视频一区二区| 亚洲经典国产精华液单| 国产精品久久久久久久电影| 22中文网久久字幕| 国产真实伦视频高清在线观看| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽人人爽人人片va| 国产亚洲最大av| 黄色怎么调成土黄色| 欧美成人一区二区免费高清观看| 成年人午夜在线观看视频| 麻豆乱淫一区二区| 如何舔出高潮| 亚洲国产精品国产精品| 久久国产精品男人的天堂亚洲 | 高清视频免费观看一区二区| 国产精品.久久久| 国产成人精品久久久久久| 欧美亚洲 丝袜 人妻 在线| 久久精品久久久久久噜噜老黄| 国产精品国产三级国产av玫瑰| 男女边吃奶边做爰视频| 国产精品精品国产色婷婷| 欧美成人一区二区免费高清观看| 国产人妻一区二区三区在| 青春草国产在线视频| 久久久久久人妻| 精品一区二区三卡| 91狼人影院| 亚洲经典国产精华液单| 1000部很黄的大片| 久久99蜜桃精品久久| 菩萨蛮人人尽说江南好唐韦庄| 久久久精品94久久精品| 一区二区三区四区激情视频| 国产亚洲精品久久久com| 国产美女午夜福利| 亚洲av综合色区一区| videos熟女内射| 午夜视频国产福利| av福利片在线观看| 精品人妻偷拍中文字幕| 一区在线观看完整版| 欧美另类一区| 久久 成人 亚洲| 久久久久网色| 欧美日韩亚洲高清精品| 中文字幕久久专区| 人妻制服诱惑在线中文字幕| 交换朋友夫妻互换小说| 亚洲欧美成人综合另类久久久| 日本猛色少妇xxxxx猛交久久| 一区二区三区免费毛片| 国产成人aa在线观看| 香蕉精品网在线| 欧美一级a爱片免费观看看| 91狼人影院| 国产成人精品久久久久久| 欧美精品一区二区免费开放| 777米奇影视久久| 啦啦啦视频在线资源免费观看| 2021少妇久久久久久久久久久| 亚洲精品中文字幕在线视频 | 国产亚洲精品久久久com| 一级av片app| 夜夜爽夜夜爽视频| 成人免费观看视频高清| 蜜臀久久99精品久久宅男| 国产精品一及| 舔av片在线| 插阴视频在线观看视频| 成人二区视频| 男人添女人高潮全过程视频| 黄色日韩在线| 这个男人来自地球电影免费观看 | 国产高清三级在线| 人人妻人人爽人人添夜夜欢视频 | 国产精品麻豆人妻色哟哟久久| 日韩成人伦理影院| 精品国产露脸久久av麻豆| 91aial.com中文字幕在线观看| 爱豆传媒免费全集在线观看| 校园人妻丝袜中文字幕| 国产精品嫩草影院av在线观看| 中文字幕免费在线视频6| av专区在线播放| 亚洲电影在线观看av| 少妇丰满av| 在线观看国产h片| 大香蕉97超碰在线| 欧美成人午夜免费资源| 成年av动漫网址| 国产淫语在线视频| 黄色欧美视频在线观看| 免费大片黄手机在线观看| 婷婷色麻豆天堂久久| 自拍欧美九色日韩亚洲蝌蚪91 | 国产av国产精品国产| av女优亚洲男人天堂| 99久久人妻综合| 欧美日本视频| 中文天堂在线官网| av免费在线看不卡| 啦啦啦视频在线资源免费观看| 人妻系列 视频| 亚洲第一av免费看| 成人亚洲精品一区在线观看 | 91久久精品国产一区二区三区| 久久韩国三级中文字幕| 亚洲av综合色区一区| 免费看不卡的av| 在线看a的网站| 亚洲综合色惰| 观看美女的网站| 狂野欧美白嫩少妇大欣赏| 国产 一区 欧美 日韩| 亚洲最大成人中文| 自拍偷自拍亚洲精品老妇| av免费在线看不卡| 美女内射精品一级片tv| 韩国高清视频一区二区三区| 大又大粗又爽又黄少妇毛片口| 大陆偷拍与自拍| 少妇丰满av| 欧美一区二区亚洲| 少妇 在线观看| 成年女人在线观看亚洲视频| 亚洲成人av在线免费| 国产一区二区在线观看日韩| 一级二级三级毛片免费看| a 毛片基地| 黄色欧美视频在线观看| 国产男女超爽视频在线观看| 免费观看a级毛片全部| 噜噜噜噜噜久久久久久91| 成年人午夜在线观看视频| 亚洲人成网站高清观看| 人妻制服诱惑在线中文字幕| 深夜a级毛片| 亚洲国产精品一区三区| 黄色视频在线播放观看不卡| 日韩视频在线欧美| 日本免费在线观看一区| 视频区图区小说| 亚洲欧美日韩卡通动漫| 久久精品夜色国产| 午夜视频国产福利| 永久免费av网站大全| 七月丁香在线播放| 一区二区三区乱码不卡18| 国产久久久一区二区三区| 大香蕉久久网| 国产熟女欧美一区二区| 日产精品乱码卡一卡2卡三| 七月丁香在线播放| 亚洲欧洲日产国产| 中文乱码字字幕精品一区二区三区| 99热这里只有是精品50| 夜夜看夜夜爽夜夜摸| 亚洲欧洲国产日韩| 高清不卡的av网站| 国产久久久一区二区三区| 日本黄色日本黄色录像| 欧美bdsm另类| 啦啦啦中文免费视频观看日本| 免费在线观看成人毛片| 蜜桃在线观看..| 亚洲精品日韩在线中文字幕| 国产亚洲午夜精品一区二区久久| 国产淫片久久久久久久久| 黄色欧美视频在线观看| 久久人妻熟女aⅴ| 亚洲av福利一区| 一区二区三区乱码不卡18| 久久精品国产a三级三级三级| av天堂中文字幕网| 99热国产这里只有精品6| 国产高清有码在线观看视频| 国产精品久久久久久精品古装| 亚洲久久久国产精品| 婷婷色麻豆天堂久久| 欧美少妇被猛烈插入视频| 最近的中文字幕免费完整| 亚洲丝袜综合中文字幕| 边亲边吃奶的免费视频| 亚洲av福利一区| 国产淫片久久久久久久久| 国产伦精品一区二区三区四那| 永久免费av网站大全| 成人漫画全彩无遮挡| 最近2019中文字幕mv第一页| 亚洲三级黄色毛片| 日韩欧美 国产精品| 欧美国产精品一级二级三级 | 成人免费观看视频高清| 搡女人真爽免费视频火全软件| 亚洲婷婷狠狠爱综合网| videos熟女内射| 永久网站在线| 各种免费的搞黄视频| 久久精品熟女亚洲av麻豆精品| 亚洲电影在线观看av| 又粗又硬又长又爽又黄的视频| 在线 av 中文字幕| 日韩一区二区视频免费看| 天天躁夜夜躁狠狠久久av| 黑人猛操日本美女一级片| 韩国av在线不卡| 国产伦在线观看视频一区| 精华霜和精华液先用哪个| 午夜福利视频精品| av在线app专区| 精品熟女少妇av免费看| 亚洲内射少妇av| 欧美xxxx性猛交bbbb| 久久久欧美国产精品| 只有这里有精品99| 久久精品久久精品一区二区三区| 九九爱精品视频在线观看| h日本视频在线播放| 美女国产视频在线观看| 下体分泌物呈黄色| 久久毛片免费看一区二区三区| 欧美日韩视频精品一区| 一二三四中文在线观看免费高清| 欧美 日韩 精品 国产| 国产黄片美女视频| 国产欧美另类精品又又久久亚洲欧美| 国产老妇伦熟女老妇高清| 成年女人在线观看亚洲视频| 日韩av不卡免费在线播放| 丰满迷人的少妇在线观看| 嘟嘟电影网在线观看| 欧美激情国产日韩精品一区| 日韩不卡一区二区三区视频在线| 夜夜骑夜夜射夜夜干| 丝袜喷水一区| 我的老师免费观看完整版| 偷拍熟女少妇极品色| 欧美成人精品欧美一级黄| 国产高清不卡午夜福利| 国产男人的电影天堂91| 国产视频内射| 九草在线视频观看| 中国三级夫妇交换| 少妇人妻精品综合一区二区| 99re6热这里在线精品视频| 免费播放大片免费观看视频在线观看| 性色avwww在线观看| 国产午夜精品一二区理论片| 在线观看免费日韩欧美大片 | 亚洲久久久国产精品| 成人国产麻豆网| 国产精品久久久久成人av| a级一级毛片免费在线观看| 国产黄频视频在线观看| 国产黄片视频在线免费观看| 蜜桃久久精品国产亚洲av| 日本免费在线观看一区| 欧美另类一区| 18禁动态无遮挡网站| 人体艺术视频欧美日本| 免费久久久久久久精品成人欧美视频 | 欧美xxⅹ黑人| 国产高清国产精品国产三级 | 国产亚洲5aaaaa淫片| 啦啦啦视频在线资源免费观看| 国产乱人偷精品视频| 高清黄色对白视频在线免费看 | 亚洲va在线va天堂va国产| 久久久精品免费免费高清| 欧美精品一区二区免费开放| 好男人视频免费观看在线| 成人美女网站在线观看视频| 国产精品国产三级专区第一集| 久久 成人 亚洲| 男人添女人高潮全过程视频| 91aial.com中文字幕在线观看| 日韩欧美精品免费久久| 成人一区二区视频在线观看| 街头女战士在线观看网站| 一二三四中文在线观看免费高清| 免费黄频网站在线观看国产| 午夜日本视频在线| 亚洲aⅴ乱码一区二区在线播放| tube8黄色片| 久久久午夜欧美精品| videossex国产| 欧美高清成人免费视频www| 日韩精品有码人妻一区| 国产片特级美女逼逼视频| 高清视频免费观看一区二区| av免费在线看不卡| a级毛色黄片| .国产精品久久| 欧美另类一区| 直男gayav资源| 亚洲一区二区三区欧美精品| 卡戴珊不雅视频在线播放| 一级毛片aaaaaa免费看小| 深夜a级毛片| 国产亚洲精品久久久com| 亚洲国产精品999| 亚洲av.av天堂| 免费在线观看成人毛片| 九九在线视频观看精品| a级毛片免费高清观看在线播放| 亚洲精品乱码久久久v下载方式| 国内少妇人妻偷人精品xxx网站| 亚洲,一卡二卡三卡| 男人添女人高潮全过程视频| 欧美日韩精品成人综合77777| 国产片特级美女逼逼视频| 一个人看视频在线观看www免费| 最近最新中文字幕免费大全7| 国产黄色视频一区二区在线观看| 国产欧美亚洲国产| 国产亚洲91精品色在线| 美女视频免费永久观看网站| 欧美xxⅹ黑人| 一区二区三区四区激情视频| 成人一区二区视频在线观看| 在线观看一区二区三区| 亚洲人与动物交配视频| 亚洲国产日韩一区二区| 亚洲精品亚洲一区二区| 亚洲aⅴ乱码一区二区在线播放| 97热精品久久久久久| 久久久久久久国产电影| 黄色日韩在线| 噜噜噜噜噜久久久久久91| 成年av动漫网址| 午夜福利在线观看免费完整高清在| 亚洲国产欧美在线一区| 国产精品久久久久久久久免| 成人毛片60女人毛片免费| 我要看日韩黄色一级片| 国产成人a区在线观看| 日本午夜av视频| 一本久久精品| 看免费成人av毛片| 一级片'在线观看视频| 国产v大片淫在线免费观看| 天天躁夜夜躁狠狠久久av| 亚洲高清免费不卡视频| 日韩电影二区| 国内揄拍国产精品人妻在线| a 毛片基地| 欧美日韩综合久久久久久| 99九九线精品视频在线观看视频| 国产精品麻豆人妻色哟哟久久| 哪个播放器可以免费观看大片| 精品人妻熟女av久视频| 亚洲第一区二区三区不卡| 噜噜噜噜噜久久久久久91| 又黄又爽又刺激的免费视频.| 黄片无遮挡物在线观看| 18禁裸乳无遮挡免费网站照片| 不卡视频在线观看欧美| 国产精品免费大片| 久久久a久久爽久久v久久| av播播在线观看一区| 国产 精品1| av天堂中文字幕网| 国产黄片视频在线免费观看| 秋霞伦理黄片| 亚州av有码| 亚洲欧美中文字幕日韩二区| 插逼视频在线观看| 国产黄片视频在线免费观看| 中国美白少妇内射xxxbb| 在线观看三级黄色| 亚洲精品一区蜜桃| 精品人妻视频免费看| av女优亚洲男人天堂| 高清午夜精品一区二区三区| 精品少妇久久久久久888优播| 午夜免费男女啪啪视频观看| 九色成人免费人妻av| 免费黄频网站在线观看国产| 最新中文字幕久久久久| 热re99久久精品国产66热6| 国内少妇人妻偷人精品xxx网站| 一本久久精品| 国产白丝娇喘喷水9色精品| 91在线精品国自产拍蜜月| 亚洲精品一二三| 97在线视频观看| 日韩电影二区| av又黄又爽大尺度在线免费看| 久久精品国产鲁丝片午夜精品| 亚洲精品日韩av片在线观看| a级毛色黄片| 日韩av在线免费看完整版不卡| 午夜福利影视在线免费观看| 亚洲成色77777| 久久97久久精品| 91狼人影院| 中国国产av一级| 久久这里有精品视频免费| 久久久久视频综合| 乱码一卡2卡4卡精品| 国产精品秋霞免费鲁丝片| 极品教师在线视频| 国产精品福利在线免费观看| 少妇人妻精品综合一区二区| 寂寞人妻少妇视频99o| 久久国内精品自在自线图片| 久久国产精品大桥未久av | 国产精品精品国产色婷婷| 精品视频人人做人人爽| 五月玫瑰六月丁香| 亚洲人成网站高清观看| av免费在线看不卡| 亚洲不卡免费看| 黄片wwwwww| 久久久久久久大尺度免费视频| 在线免费十八禁| 插逼视频在线观看| 亚洲成人手机| 大香蕉久久网| 欧美zozozo另类| 久久久久国产网址| 久久青草综合色| 欧美+日韩+精品| 久久ye,这里只有精品| 日韩欧美一区视频在线观看 | 亚洲伊人久久精品综合| 99热这里只有精品一区| 免费观看a级毛片全部| 亚洲精品456在线播放app| 午夜日本视频在线| 欧美xxxx性猛交bbbb| 久久久午夜欧美精品| 少妇人妻精品综合一区二区| 成人毛片a级毛片在线播放| 嫩草影院入口| 亚洲成人av在线免费| 能在线免费看毛片的网站| 欧美日韩综合久久久久久| 欧美日韩国产mv在线观看视频 | 91久久精品国产一区二区成人| 日韩国内少妇激情av| 人人妻人人添人人爽欧美一区卜 | freevideosex欧美| 少妇人妻精品综合一区二区| 久久99精品国语久久久| 少妇人妻精品综合一区二区| 亚洲精品视频女| 国产成人aa在线观看| 22中文网久久字幕| av女优亚洲男人天堂| 午夜福利高清视频| 少妇人妻精品综合一区二区| 午夜福利在线在线| 美女高潮的动态| 日本午夜av视频| 久久久久久人妻| 少妇人妻久久综合中文| 精品亚洲成国产av| 777米奇影视久久| 国产欧美日韩精品一区二区| 国产精品三级大全| 秋霞伦理黄片| 男女国产视频网站| 国产精品不卡视频一区二区| 国产高清国产精品国产三级 | 国产欧美另类精品又又久久亚洲欧美| videos熟女内射| 亚州av有码| 国产伦精品一区二区三区四那| 国产视频首页在线观看| 国产毛片在线视频| 久久精品久久精品一区二区三区| 黄色一级大片看看| 国模一区二区三区四区视频| 熟女av电影| 一级毛片久久久久久久久女| 十八禁网站网址无遮挡 | 丰满迷人的少妇在线观看| 久久久久久久久久成人| 亚洲精品久久午夜乱码| 色综合色国产| 欧美日韩国产mv在线观看视频 | 卡戴珊不雅视频在线播放| 网址你懂的国产日韩在线| 国模一区二区三区四区视频| 高清午夜精品一区二区三区|