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

    風(fēng)電采用MGP并網(wǎng)的小干擾建模和阻尼特性

    2017-12-22 08:49:44衛(wèi)思明黃永章
    電力系統(tǒng)自動(dòng)化 2017年22期
    關(guān)鍵詞:同步機(jī)功角狀態(tài)變量

    衛(wèi)思明,黃永章

    (新能源電力系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室(華北電力大學(xué)),北京市 102206)

    風(fēng)電采用MGP并網(wǎng)的小干擾建模和阻尼特性

    衛(wèi)思明,黃永章

    (新能源電力系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室(華北電力大學(xué)),北京市 102206)

    針對(duì)高比例新能源給電網(wǎng)帶來(lái)的低慣量和弱阻尼問(wèn)題,提出新能源同步機(jī)的并網(wǎng)方式,即新能源驅(qū)動(dòng)同步電動(dòng)機(jī)—同步發(fā)電機(jī)(MGP)后并網(wǎng)。詳細(xì)闡述了MGP的小干擾建模方法。首先,分析MGP的功角關(guān)系并考慮勵(lì)磁系統(tǒng)動(dòng)態(tài),建立包含自動(dòng)電壓調(diào)節(jié)器和電力系統(tǒng)穩(wěn)定器的小干擾模型,指出MGP存在阻尼轉(zhuǎn)矩疊加效應(yīng)。然后,考慮雙饋風(fēng)機(jī)狀態(tài)變量,建立其與MGP對(duì)接的完整的小干擾狀態(tài)方程。在此基礎(chǔ)上,推導(dǎo)了基于K系數(shù)的動(dòng)態(tài)特性和轉(zhuǎn)矩分量表示法,揭示MGP的阻尼轉(zhuǎn)矩疊加效應(yīng)的本質(zhì)是雙勵(lì)磁耦合特性。最后,基于模型建立單機(jī)無(wú)窮大算例,對(duì)比了MGP和傳統(tǒng)機(jī)組的小干擾穩(wěn)定特性;在雙饋風(fēng)機(jī)采用傳統(tǒng)和MGP兩種并網(wǎng)方式下,研究了雙勵(lì)磁耦合特性對(duì)風(fēng)機(jī)主導(dǎo)的振蕩模態(tài)的影響。結(jié)果表明,提出的小干擾建模方法可以準(zhǔn)確反映MGP系統(tǒng)的雙勵(lì)磁耦合特性,該特性能夠有效提升阻尼轉(zhuǎn)矩分量和改善系統(tǒng)振蕩模態(tài)。

    高比例新能源;同步電動(dòng)機(jī)—同步發(fā)電機(jī);小干擾建模;阻尼特性

    0 引言

    新能源以高比例接入電網(wǎng)給系統(tǒng)穩(wěn)定性帶來(lái)巨大挑戰(zhàn),其中最受關(guān)注的是慣性和阻尼問(wèn)題。眾所周知,同步發(fā)電機(jī)(SG)的轉(zhuǎn)動(dòng)慣量和勵(lì)磁系統(tǒng)可分別提供慣性響應(yīng)和阻尼控制,能夠有效維持頻率和小干擾穩(wěn)定[1]。而以變速風(fēng)機(jī)和光伏為代表的新能源由于自身結(jié)構(gòu)、波動(dòng)性和運(yùn)行控制策略等方面的特點(diǎn),導(dǎo)致慣性和阻尼比較欠缺,這在接入比例較高時(shí)尤為突出。

    國(guó)內(nèi)外圍繞該問(wèn)題開展了大量研究。文獻(xiàn)[2-4]提出虛擬同步機(jī)技術(shù),在換流器控制環(huán)節(jié)引入同步機(jī)的方程,使其能夠模擬同步機(jī)的慣性、阻尼和調(diào)壓等特性,對(duì)風(fēng)機(jī)來(lái)說(shuō)主要是利用轉(zhuǎn)子動(dòng)能,對(duì)光伏來(lái)說(shuō)主要是利用儲(chǔ)能環(huán)節(jié)中的能量。但虛擬同步機(jī)本質(zhì)上是電力電子裝置,無(wú)法真正像同步機(jī)一樣具有高短路電流支撐能力和強(qiáng)電壓控制能力[5-6],為提供充足的慣性能量而額外配置的儲(chǔ)能環(huán)節(jié)也提高了成本;文獻(xiàn)[7-10]基于相似的思路提出風(fēng)機(jī)的一次調(diào)頻和多端直流系統(tǒng)參與慣性響應(yīng)的控制策略;文獻(xiàn)[11-12]研究了在各類風(fēng)機(jī)中引入附加阻尼控制器,使風(fēng)機(jī)的功率輸出中包含與轉(zhuǎn)速成比例的分量來(lái)增加阻尼,本質(zhì)上仍是通過(guò)改變風(fēng)機(jī)轉(zhuǎn)速為系統(tǒng)提供阻尼轉(zhuǎn)矩。

    可以看出,不管是虛擬同步技術(shù)還是附加阻尼控制,本質(zhì)都是通過(guò)新能源與電網(wǎng)的能量交互實(shí)現(xiàn)的。但新能源換流器不僅受到自身運(yùn)行條件(如最大功率追蹤)和暫態(tài)過(guò)流能力、短路容量等方面的制約,其與系統(tǒng)潮流的動(dòng)態(tài)交互也具有很強(qiáng)的不確定性,如文獻(xiàn)[13-15]指出虛擬慣量控制在某些工況下對(duì)小干擾穩(wěn)定產(chǎn)生不利影響,文獻(xiàn)[16-18]指出風(fēng)機(jī)類型、控制策略、滲透率和運(yùn)行區(qū)域等諸多因素都會(huì)對(duì)穩(wěn)定性造成不同程度的影響。上述原因使新能源往往難以提供可靠的慣性響應(yīng)和阻尼轉(zhuǎn)矩以維持電網(wǎng)穩(wěn)定運(yùn)行。

    綜合來(lái)看,以上解決方案都在一定程度將同步機(jī)的運(yùn)動(dòng)方程、慣性響應(yīng)、下垂特性和阻尼控制等思想借鑒到換流器的控制中。可見,雖然電網(wǎng)的電力電子化程度越來(lái)越高,但傳統(tǒng)同步機(jī)的相關(guān)理論仍發(fā)揮著重要作用,同步機(jī)的一些優(yōu)勢(shì)也是換流器不具備的。未來(lái)的高比例新能源電網(wǎng)是否一定要拋棄同步機(jī)、如何利用同步機(jī)提升電網(wǎng)可靠性是值得探索的。

    因此,也有諸多方案嘗試采用真正的電機(jī)系統(tǒng)解決新能源電網(wǎng)穩(wěn)定問(wèn)題,如調(diào)相機(jī)和電動(dòng)機(jī)—發(fā)電機(jī)系統(tǒng)。調(diào)相機(jī)不僅能夠進(jìn)行無(wú)功補(bǔ)償,還能提供短路容量和轉(zhuǎn)動(dòng)慣量,文獻(xiàn)[19-20]考察了其對(duì)頻率響應(yīng)和短路比的提升作用;電動(dòng)機(jī)—發(fā)電機(jī)系統(tǒng)將兩臺(tái)電機(jī)同軸連接,在電網(wǎng)發(fā)展早期已有應(yīng)用[21-25],所采用的電機(jī)類型包括同步、異步和直流電機(jī),主要用于能量變換。近年來(lái)隨著新能源的發(fā)展,電動(dòng)機(jī)—發(fā)電機(jī)系統(tǒng)由于轉(zhuǎn)動(dòng)慣量、過(guò)載能力和電氣隔離等方面的優(yōu)勢(shì)在提升電網(wǎng)穩(wěn)定性的場(chǎng)合中得到應(yīng)用。文獻(xiàn)[21-22]研究了同步型電動(dòng)機(jī)—發(fā)電機(jī)系統(tǒng)的結(jié)構(gòu)特點(diǎn)、運(yùn)行方式和并聯(lián)運(yùn)行的勵(lì)磁控制等問(wèn)題;文獻(xiàn)[23-25]采用異步型電動(dòng)機(jī)—發(fā)電機(jī)系統(tǒng),研究了改善系統(tǒng)動(dòng)態(tài)過(guò)程的控制策略,為含高比例新能源的微網(wǎng)提供轉(zhuǎn)動(dòng)慣量和頻率支撐;文獻(xiàn)[26]將各類新能源匯聚到直流母線上,再通過(guò)直流電動(dòng)機(jī)和同步發(fā)電機(jī)構(gòu)成的系統(tǒng)并網(wǎng),研究了其對(duì)諧波、電壓和頻率的改善作用。

    鑒于同步機(jī)的諸多優(yōu)勢(shì),本文提出由同步電動(dòng)機(jī)(SM)和同步發(fā)電機(jī)構(gòu)成的新能源同步機(jī)系統(tǒng)實(shí)現(xiàn)風(fēng)電并網(wǎng)的方案。在前期工作的基礎(chǔ)上,詳細(xì)闡述了新能源同步機(jī)的小干擾建模方法,通過(guò)理論推導(dǎo)揭示能反映其阻尼特性的雙勵(lì)磁耦合機(jī)理,最后用算例驗(yàn)證了小干擾模型的有效性,對(duì)比研究了雙勵(lì)磁耦合特性對(duì)阻尼轉(zhuǎn)矩的提升作用和對(duì)風(fēng)機(jī)主導(dǎo)的振蕩模態(tài)的影響。

    1 新能源同步機(jī)

    新能源(以風(fēng)電為例)采用同步電動(dòng)機(jī)—同步發(fā)電機(jī)(MGP)實(shí)現(xiàn)并網(wǎng)的結(jié)構(gòu)如圖1(a)所示,多臺(tái)風(fēng)機(jī)作為一個(gè)單元,發(fā)出的電能經(jīng)匯聚和升壓等環(huán)節(jié)后驅(qū)動(dòng)同步電動(dòng)機(jī),同步電動(dòng)機(jī)作為原動(dòng)機(jī),帶動(dòng)同步發(fā)電機(jī)發(fā)電和并網(wǎng)。圖中:AVR表示自動(dòng)電壓調(diào)節(jié)器;PSS表示電力系統(tǒng)穩(wěn)定器。

    圖1(a)中,MGP環(huán)節(jié)的能量傳輸形式是從電能(電動(dòng)機(jī)電磁轉(zhuǎn)矩TeM)到機(jī)械能(軸系機(jī)械轉(zhuǎn)矩Tm)再到電能(發(fā)電機(jī)電磁轉(zhuǎn)矩TeG),而傳統(tǒng)火電機(jī)組是從化學(xué)能到機(jī)械能再到電能,這個(gè)主要的區(qū)別使MGP在運(yùn)行效率上具有優(yōu)勢(shì):由于鍋爐中煤燃燒的不穩(wěn)定性,傳統(tǒng)火力發(fā)電機(jī)組在低于50%額定功率運(yùn)行時(shí)的效率將大幅降低,MGP系統(tǒng)只包含兩臺(tái)同步機(jī),根據(jù)效率實(shí)測(cè)曲線,即使大幅度降額運(yùn)行也不會(huì)使效率過(guò)多地偏離額定運(yùn)行點(diǎn),這對(duì)應(yīng)于風(fēng)機(jī)單元出力較低時(shí)的工況;而當(dāng)MGP額定運(yùn)行,即對(duì)應(yīng)風(fēng)機(jī)單元出力較高時(shí),該環(huán)節(jié)的運(yùn)行效率可以達(dá)到96%甚至更高,以200 MW大型同步發(fā)電機(jī)為例,占比最高的機(jī)械損耗為1 332 kW,占0.67%,總損耗為3 404 kW,占1.7%[27]。因此,MGP并網(wǎng)方式相當(dāng)于在新能源和大電網(wǎng)之間加入兩臺(tái)同步機(jī)環(huán)節(jié),就這個(gè)環(huán)節(jié)本身的效率而言,MGP系統(tǒng)運(yùn)行時(shí)不會(huì)過(guò)多地增加總的電能損耗。

    圖1 MGP及其構(gòu)成的高比例新能源電網(wǎng)形態(tài)Fig.1 MGP system and a grid configuration with high penetration of renewable energy

    基于MGP并網(wǎng)方式構(gòu)建的高比例新能源電網(wǎng)的一種可能形態(tài)如圖1(b)所示,一部分新能源通過(guò)MGP系統(tǒng)并網(wǎng),另一部分仍采用傳統(tǒng)方式并網(wǎng),協(xié)調(diào)兩種并網(wǎng)方式的容量配比以實(shí)現(xiàn)最優(yōu)穩(wěn)定運(yùn)行特性。這時(shí),電網(wǎng)呈現(xiàn)出以同步機(jī)為核心的形態(tài),同步機(jī)的轉(zhuǎn)動(dòng)慣量和短路電流支撐能力彌補(bǔ)了新能源的不足,勵(lì)磁控制系統(tǒng)能夠改善電網(wǎng)動(dòng)態(tài)特性;經(jīng)濟(jì)性方面,大容量同步電機(jī)易于制造,同容量下,同步機(jī)的成本遠(yuǎn)低于儲(chǔ)能系統(tǒng),退役的火電機(jī)組可以用于MGP系統(tǒng),避免了再次投資。MGP系統(tǒng)的這些優(yōu)勢(shì)為提升高比例新能源電網(wǎng)穩(wěn)定性提供了一種可能的解決方案。

    2 考慮勵(lì)磁系統(tǒng)動(dòng)態(tài)的MGP系統(tǒng)小干擾建模方法

    2.1 MGP系統(tǒng)的功角關(guān)系

    功角關(guān)系是同步機(jī)建模、控制和運(yùn)行的基礎(chǔ)。MGP系統(tǒng)中的兩臺(tái)同步機(jī)在機(jī)械上通過(guò)轉(zhuǎn)軸實(shí)現(xiàn)連接,反映到電氣上表現(xiàn)為兩機(jī)空載電勢(shì)EM′和EG′的同時(shí)、同向和同速旋轉(zhuǎn),如圖2(a)所示。將電動(dòng)機(jī)和發(fā)電機(jī)兩側(cè)的母線電壓UBM和UBG作為測(cè)量功角的參考量,忽略轉(zhuǎn)子位置偏差,功角δM和δG滿足如下關(guān)系:

    δM+δG=δMG

    (1)

    式中:δMG為兩機(jī)功角之和。

    圖2 MGP功角關(guān)系及其單機(jī)無(wú)窮大系統(tǒng)Fig.2 Rotor angle relation of MGP and its single machine infinite bus system

    一方面,在圖2中規(guī)定的功角參考方式下,已開展的研究工作提出了基于源網(wǎng)相位差的有功功率控制策略[28],在這種控制模式下,源網(wǎng)相位差與δMG相等,δMG作為控制參數(shù)隨有功功率的變化而變化,如圖2(a)所示;另一方面,功角變化導(dǎo)致有功功率變化的本質(zhì)是引起了磁通的變化,從而造成了氣隙電磁轉(zhuǎn)矩的變化,反之亦然,式(1)中的功角耦合關(guān)系使MGP這一特殊的雙同步機(jī)系統(tǒng)有別于單同步機(jī),每個(gè)功角的變化都會(huì)通過(guò)轉(zhuǎn)軸對(duì)另一個(gè)功角產(chǎn)生影響,進(jìn)而影響兩臺(tái)電機(jī)的磁場(chǎng)特性。綜上,為了便于討論控制參數(shù)對(duì)小干擾穩(wěn)定的影響并體現(xiàn)功角耦合關(guān)系,本文在進(jìn)行小干擾建模時(shí)選取δMG作為狀態(tài)變量。

    2.2 考慮勵(lì)磁系統(tǒng)動(dòng)態(tài)的MGP小干擾模型

    新能源(以雙饋風(fēng)機(jī)為例)與MGP構(gòu)成的單機(jī)無(wú)窮大系統(tǒng)如圖2(b)所示。MGP穩(wěn)定運(yùn)行時(shí),兩臺(tái)同步機(jī)同軸連接,不考慮扭振等特殊工況時(shí),各自的機(jī)械轉(zhuǎn)矩是相等的,體現(xiàn)為轉(zhuǎn)速相同;兩臺(tái)同步機(jī)的容量和轉(zhuǎn)子結(jié)構(gòu)相同,體現(xiàn)為兩者的慣性常數(shù)相等,因此,電動(dòng)機(jī)和發(fā)電機(jī)的運(yùn)動(dòng)方程可分別表示為:

    (2)

    (3)

    式中:D為微分算子d/dt;Δωr為轉(zhuǎn)速偏差;ω0為基準(zhǔn)轉(zhuǎn)子電氣速度;H′為MGP中單臺(tái)同步機(jī)的慣性常數(shù);KDM和KDG分別為電動(dòng)機(jī)和發(fā)電機(jī)的阻尼系數(shù)。

    將式(2)和式(3)中轉(zhuǎn)速的微分方程相加并線性化,得到以Δωr和ΔδG為狀態(tài)變量的狀態(tài)方程為:

    (4)

    從式(4)可以直觀地看出,MGP的阻尼轉(zhuǎn)矩分量由兩機(jī)各自的阻尼轉(zhuǎn)矩分量疊加而成,疊加效應(yīng)的必然結(jié)果是整體的阻尼水平高于單臺(tái)同步機(jī)。

    為了將式(4)中的狀態(tài)變量ΔδG替換為ΔδMG,還需利用MGP系統(tǒng)的功率平衡關(guān)系,在采用標(biāo)幺值的穩(wěn)定性分析中,假設(shè)轉(zhuǎn)速變化對(duì)電壓沒(méi)有明顯的影響(ωr=ω0),則電磁轉(zhuǎn)矩等于電磁功率,即

    TeM-TeG=pmloss

    (5)

    式中:pmloss為MGP系統(tǒng)中間軸系環(huán)節(jié)的損耗。

    可認(rèn)為pmloss在受到小擾動(dòng)時(shí)是恒定的,因此,對(duì)式(1)和式(5)線性化可得到:

    (6)

    考慮勵(lì)磁系統(tǒng)動(dòng)態(tài)時(shí),式(6)中的兩個(gè)電磁轉(zhuǎn)矩可寫成如下形式:

    (7)

    式中:具有不同下標(biāo)的系數(shù)K表示與電機(jī)參數(shù)和磁鏈初始條件有關(guān)的常數(shù),下同。

    聯(lián)立式(6)和式(7)可求得:

    (8)

    可以看出,考慮勵(lì)磁后,δM和δG與兩機(jī)的勵(lì)磁磁鏈ψfdM和ψfdG都建立了聯(lián)系,這反映了MGP中每臺(tái)電機(jī)的磁場(chǎng)變化都會(huì)影響另一臺(tái)電機(jī)的功角,從而影響另一臺(tái)電機(jī)的磁場(chǎng)。

    將式(7)和式(8)代入式(4)中可將δMG的微分方程表示為:

    (9)

    可以看出,式(9)中包含DΔψfdM和DΔψfdG兩項(xiàng),而ΔψfdM和ΔψfdG也是狀態(tài)變量,表示兩機(jī)勵(lì)磁磁場(chǎng)動(dòng)態(tài)的影響,對(duì)應(yīng)的微分方程如式(10)所示,為簡(jiǎn)化推導(dǎo)過(guò)程,以下均采用線性化后的方程:

    (10)

    式中:EfdM和EfdG分別為兩機(jī)勵(lì)磁電壓;ifdM和ifdG分別為兩機(jī)勵(lì)磁電流;具有不同下標(biāo)的系數(shù)A表示與電機(jī)參數(shù)有關(guān)的常數(shù),下同。

    式(10)中的勵(lì)磁電流可表示為:

    (11)

    式中:具有不同下標(biāo)的系數(shù)F表示與電機(jī)參數(shù)有關(guān)的常數(shù)。

    利用式(8)將式(11)中的ΔδM和ΔδG替換為ΔδMG,再代入式(10)中可得到包含ΔEfdM和ΔEfdG的勵(lì)磁磁鏈方程??紤]AVR和PSS后,需將ΔEfdM和ΔEfdG用狀態(tài)變量表示,下面以附錄A所示的簡(jiǎn)單勵(lì)磁系統(tǒng)為例進(jìn)行推導(dǎo)。

    首先僅考慮AVR,ΔEfdM和ΔEfdG可表示為:

    (12)

    式中:KAM和KAG為各自的勵(lì)磁機(jī)增益。

    將式(12)代入式(10)中即可將勵(lì)磁磁鏈的微分方程用全部狀態(tài)變量來(lái)表示,把求得的DΔψfdM和DΔψfdG兩項(xiàng)代入式(9),可得到完整的功角微分方程。

    另外,對(duì)于附錄A中的勵(lì)磁系統(tǒng),還需要增加兩個(gè)分別屬于電動(dòng)機(jī)和發(fā)電機(jī)的微分方程:

    (13)

    式中:具有不同下標(biāo)的Δv表示與電機(jī)參數(shù)有關(guān)的電壓傳感器輸出變化量,下同;ΔEtM和ΔEtG為端電壓變化量;具有不同下標(biāo)的T表示與電機(jī)參數(shù)有關(guān)的時(shí)間常數(shù),下同。

    結(jié)合網(wǎng)絡(luò)方程,式(13)中的端電壓可表示為:

    (14)

    式中:具有不同下標(biāo)的系數(shù)B表示與電機(jī)參數(shù)有關(guān)的常數(shù),下同。

    用式(8)替換掉式(14)中的ΔδM和ΔδG,并代入式(13),即可得到用所有狀態(tài)變量表示的微分方程。

    綜上,考慮含AVR的勵(lì)磁系統(tǒng)后,狀態(tài)矩陣變?yōu)榱A,如式(15)所示,矩陣中系數(shù)的詳細(xì)表達(dá)式見附錄B。

    (15)

    從上述建模過(guò)程可以看出,描述功角和磁鏈關(guān)系的式(8)是推導(dǎo)所有微分方程的基礎(chǔ);ΔδMG的微分方程是MGP和傳統(tǒng)發(fā)電機(jī)建模時(shí)最主要的區(qū)別:對(duì)于傳統(tǒng)發(fā)電機(jī)的任意高階模型,狀態(tài)矩陣第2行除了a21以外的元素均為0,而式(15)中第2行元素都不為0,說(shuō)明ΔδMG與每個(gè)狀態(tài)變量都產(chǎn)生耦合。

    造成這個(gè)現(xiàn)象的原因可通過(guò)式(6)和式(7)來(lái)解釋。式(6)中的電磁轉(zhuǎn)矩平衡方程代表了兩臺(tái)同步機(jī)在機(jī)械上的耦合,而式(7)中每臺(tái)同步機(jī)的電磁轉(zhuǎn)矩都與各自磁場(chǎng)變化和功角變化相關(guān),式(6)中的δMG包含了兩個(gè)功角的信息,因此δMG必然與兩機(jī)的磁場(chǎng)都產(chǎn)生耦合。從廣義上講,凡是能引起磁場(chǎng)變化的狀態(tài)變量,都將與δMG耦合,這體現(xiàn)在狀態(tài)矩陣上就是第2行元素均不為0。

    因此,MGP作為一個(gè)整體,采用δMG作為狀態(tài)變量能夠描述兩臺(tái)同步機(jī)的電磁耦合關(guān)系,而用δM或δG作為狀態(tài)變量則會(huì)造成電磁耦合信息的缺失。本文將這種關(guān)系稱為“雙勵(lì)磁耦合特性”。

    此外從式(15)還可以看出,對(duì)于含勵(lì)磁系統(tǒng)的MGP小干擾模型,電動(dòng)機(jī)和發(fā)電機(jī)的狀態(tài)變量總是成對(duì)出現(xiàn)的,這導(dǎo)致在電機(jī)模型的階數(shù)相同時(shí),MGP的狀態(tài)矩陣總是比單臺(tái)同步機(jī)的狀態(tài)矩陣多出兩階。

    根據(jù)以上建模思路,為了加入雙機(jī)的PSS,需要增加4個(gè)微分方程,這里僅以同步電動(dòng)機(jī)為例闡述建模方法,微分方程如下:

    (16)

    (17)

    式中:KSTABM為PSS增益。

    將式(4)中的轉(zhuǎn)速方程代入式(16)中,再將式(16)代入式(17)中消去DΔv2M,即可將式(16)和式(17)用全部狀態(tài)變量表示。

    此外,考慮PSS后的勵(lì)磁電壓可表示為:

    ΔEfdM=KAM(ΔvsM-Δv1M)

    (18)

    將式(18)代入式(9)和式(10)中對(duì)功角和勵(lì)磁磁鏈的微分方程進(jìn)行修正,即可得到考慮PSS的狀態(tài)方程:

    (19)

    2.3 考慮雙饋風(fēng)機(jī)狀態(tài)變量的小干擾模型

    以上建立了以MGP的運(yùn)行控制參數(shù)為狀態(tài)變量、TeM為輸入變量的狀態(tài)方程。在圖2(b)所示的系統(tǒng)中,雙饋風(fēng)機(jī)連接MGP這一環(huán)節(jié)整體作為一臺(tái)單機(jī)接入無(wú)窮大母線,必須考慮風(fēng)機(jī)自身的運(yùn)行控制參數(shù)對(duì)小干擾穩(wěn)定的影響。因此,本節(jié)對(duì)式(19)進(jìn)行擴(kuò)充和修正以包含風(fēng)機(jī)所有狀態(tài)變量。

    風(fēng)機(jī)通過(guò)MGP并網(wǎng)后的輸出功率即為MGP系統(tǒng)電動(dòng)機(jī)側(cè)的輸入功率,滿足轉(zhuǎn)矩平衡關(guān)系:

    TeDFIG=TeM

    (20)

    式中:TeDFIG為風(fēng)機(jī)輸出電磁轉(zhuǎn)矩。

    類似的,雙饋風(fēng)機(jī)的定子電流與換流器網(wǎng)側(cè)電流之和即為MGP系統(tǒng)電動(dòng)機(jī)的定子電流,在統(tǒng)一的參考坐標(biāo)下滿足:

    (21)

    式中:ids和iqs分別為風(fēng)機(jī)定子dq軸電流;idg和iqg分別為換流器網(wǎng)側(cè)dq軸電流;idM和iqM分別為同步電動(dòng)機(jī)定子dq軸電流。

    通常將雙饋風(fēng)機(jī)小干擾模型用下式描述,矩陣元素詳見附錄D。

    DxDFIG=ADFIGxDFIG+BDFIGuDFIG

    (22)

    對(duì)于MGP中的同步電動(dòng)機(jī)來(lái)說(shuō),用擾動(dòng)值表示的定子電流具有如下形式[1]:

    (23)

    式中:系數(shù)m1,m2,n1,n2為與電動(dòng)機(jī)初始條件和電氣參數(shù)有關(guān)的系數(shù),具體表示見附錄D。

    仍然利用式(8)消去式(23)中的ΔδM,再聯(lián)立式(21)和式(23),可將風(fēng)機(jī)狀態(tài)方程中的輸入變量idg和iqg表示為:

    (24)

    式中:a為修正后的系數(shù)矩陣。

    從式(24)可以看出,通過(guò)聯(lián)立轉(zhuǎn)矩平衡方程,建立了網(wǎng)側(cè)電流與MGP的功角和磁場(chǎng)狀態(tài)變量的聯(lián)系,將式(24)代入式(22)中消去輸入變量idg和iqg,從而建立xDFIG與xMGP的關(guān)系。

    另外,從式(15)和式(19)中可以看出,在MGP的狀態(tài)方程中,只有Δωr和ΔTeM相關(guān)(b11≠0),因此僅需要修正轉(zhuǎn)速方程使其與xDFIG建立聯(lián)系。對(duì)雙饋風(fēng)機(jī)來(lái)說(shuō),線性化的輸出電磁轉(zhuǎn)矩可以表示為[29-30]:

    ΔTeDFIG=k1ΔEd′+k2ΔEq′+k3Δidr+k4Δiqr

    (25)

    式中:系數(shù)k1,k2,k3,k4為與雙饋風(fēng)機(jī)參數(shù)有關(guān)的系數(shù),具體表示見附錄D。

    將式(25)與轉(zhuǎn)矩平衡方程式(20)聯(lián)立,并代入式(19)中消去輸入變量ΔTeM,得到修正后的轉(zhuǎn)速微分方程:

    (26)

    式中:a′為修正后的系數(shù)矩陣。

    綜合上述推導(dǎo)可以得到雙饋風(fēng)機(jī)通過(guò)MGP并網(wǎng)的狀態(tài)方程為:

    (27)

    式中:AMGP為式(19)中的狀態(tài)矩陣;A12′和A21′為新增加的參數(shù)矩陣,詳見附錄D;AMGP為MGP原有的狀態(tài)矩陣。

    由式(27)可以看出,A12′和A21′的元素都不全為零,這表明了包括MGP的勵(lì)磁磁場(chǎng)在內(nèi)的部分狀態(tài)變量對(duì)DFIG產(chǎn)生的影響。

    3 MGP系統(tǒng)的動(dòng)態(tài)特性描述和雙勵(lì)磁耦合機(jī)理

    為揭示雙勵(lì)磁耦合特性的機(jī)理,重新定義描述MGP系統(tǒng)動(dòng)態(tài)特性的K系數(shù)。首先聯(lián)立式(7)和式(8)得到發(fā)電機(jī)輸出電磁轉(zhuǎn)矩:

    ΔTeG=K1ΔδMG+K2ΔψfdG+K3ΔψfdM

    (28)

    其中

    可以看出,考慮勵(lì)磁后,電磁轉(zhuǎn)矩同時(shí)包含了雙機(jī)磁場(chǎng)變化的信息,這正是雙勵(lì)磁耦合特性在輸出電磁轉(zhuǎn)矩上的最終體現(xiàn)。

    對(duì)相關(guān)方程進(jìn)行類似的變換,得到圖3所示的MGP系統(tǒng)動(dòng)態(tài)特性框圖,各K系數(shù)表達(dá)式見附錄E。

    圖3 MGP動(dòng)態(tài)特性的框圖Fig.3 Block diagram of dynamic characteristics for MGP

    圖3中,雙機(jī)勵(lì)磁磁通的變化不僅由ΔδMG的反饋分別通過(guò)系數(shù)K10和K11引起,每一臺(tái)同步機(jī)的勵(lì)磁磁通還由另一臺(tái)同步機(jī)分別通過(guò)系數(shù)K8和K9引起,這種雙機(jī)勵(lì)磁耦合特性最終體現(xiàn)在電磁轉(zhuǎn)矩ΔTeG上,與式(28)對(duì)應(yīng)。因此,對(duì)于MGP系統(tǒng),ΔTeG中的同步轉(zhuǎn)矩分量和阻尼轉(zhuǎn)矩分量是由雙機(jī)勵(lì)磁共同疊加引起的,這從本質(zhì)上解釋了MGP系統(tǒng)機(jī)械耦合背后的勵(lì)磁耦合關(guān)系,雙勵(lì)磁耦合特性的存在將顯著提升MGP系統(tǒng)的阻尼轉(zhuǎn)矩。

    根據(jù)上述K系數(shù)表示法,由δMG變化引起的雙機(jī)磁通變化而造成的氣隙轉(zhuǎn)矩變化可以表示為:

    (29)

    式中:系數(shù)L為各個(gè)K系數(shù)和微分因子s的函數(shù),具體表達(dá)式見附錄E。

    根據(jù)式(29)可求出在穩(wěn)態(tài)時(shí),由ΔψfdM和ΔψfdG引起的轉(zhuǎn)矩分量為:

    ΔTeG=ΔδMG·

    (30)

    對(duì)于某一振蕩模態(tài),轉(zhuǎn)矩分量可以表示為:

    ΔTeG|s=a+jb=

    (31)

    式中:a+jb為該振蕩模態(tài)對(duì)應(yīng)的特征值。

    在小干擾穩(wěn)定分析中,可用以上各式計(jì)算MGP系統(tǒng)輸出的電磁轉(zhuǎn)矩分量。

    4 算例分析

    為驗(yàn)證小干擾模型的有效性和MGP的雙勵(lì)磁耦合特性對(duì)風(fēng)機(jī)并網(wǎng)振蕩模態(tài)的影響,對(duì)圖2(b)中的單機(jī)無(wú)窮大系統(tǒng)進(jìn)行算例分析,系統(tǒng)參數(shù)如附錄A所示。

    4.1 MGP與傳統(tǒng)機(jī)組的小干擾穩(wěn)定對(duì)比

    首先,對(duì)式(4)中采用同步機(jī)經(jīng)典模型描述的MGP小干擾模型進(jìn)行驗(yàn)證。式(4)中的阻尼系數(shù)KD表示了所有阻尼轉(zhuǎn)矩分量的效果,為分析MGP的阻尼特性,設(shè)置每臺(tái)發(fā)電機(jī)和傳統(tǒng)機(jī)組有相同的阻尼系數(shù),不同工況下的小干擾計(jì)算結(jié)果如附錄A所示。對(duì)功角施加5°的擾動(dòng),ΔδMG的時(shí)間響應(yīng)如圖4(a)所示。

    當(dāng)KDM=KDM=10時(shí),表示雙機(jī)勵(lì)磁共同作用,從附錄A和圖4可以看出:此時(shí)的阻尼比(0.276)是傳統(tǒng)機(jī)組(0.112)的2倍以上,MGP阻尼振蕩的能力顯著增強(qiáng);當(dāng)KDM=0或KDM=0時(shí),表示只有一臺(tái)同步機(jī)的勵(lì)磁提供阻尼轉(zhuǎn)矩分量,兩種情況的結(jié)果是一致的,即阻尼比(0.138)略高于傳統(tǒng)機(jī)組,但明顯低于雙機(jī)勵(lì)磁,阻尼振蕩的效果介于兩者之間。以上結(jié)果說(shuō)明雙勵(lì)磁耦合特性帶來(lái)的轉(zhuǎn)矩疊加效應(yīng)使MGP具有更高的阻尼水平。

    圖4 MGP和傳統(tǒng)機(jī)組小干擾穩(wěn)定結(jié)果對(duì)比Fig.4 Comparison of MGP and traditional power unit in small signal stability

    此外,從圖4(a)還可以看出MGP的慣性響應(yīng)能力,在擾動(dòng)初期,機(jī)組依靠慣性限制功角變化的速率,MGP時(shí)間響應(yīng)曲線的斜率略高于傳統(tǒng)機(jī)組,說(shuō)明MGP的總轉(zhuǎn)動(dòng)慣量比同容量的單臺(tái)同步機(jī)組略低,但整體水平接近,這對(duì)高比例新能源電網(wǎng)的頻率穩(wěn)定是有利的。

    其次,對(duì)考慮雙勵(lì)磁系統(tǒng)的MGP小干擾模型進(jìn)行驗(yàn)證,對(duì)應(yīng)式(15)所示的模型。該模型下的阻尼效果不能通過(guò)直接設(shè)置KD得到,而是主要受到勵(lì)磁機(jī)增益KA的影響,阻尼轉(zhuǎn)矩分量需要通過(guò)將特征值代入式(31)中求得。

    選取與MGP的功角和轉(zhuǎn)速?gòu)?qiáng)相關(guān)的振蕩模態(tài),對(duì)電動(dòng)機(jī)和發(fā)電機(jī)分別設(shè)置勵(lì)磁機(jī)增益,計(jì)算輸出阻尼轉(zhuǎn)矩分量,并與傳統(tǒng)機(jī)組計(jì)算結(jié)果進(jìn)行對(duì)比,計(jì)算結(jié)果見附錄A和圖4(b)。

    從圖4(b)可以看出,對(duì)于傳統(tǒng)機(jī)組來(lái)說(shuō),隨著KA的增大,阻尼轉(zhuǎn)矩分量先降低,對(duì)于中等或快速反應(yīng)的勵(lì)磁系統(tǒng)(KA較高)則引入負(fù)的阻尼轉(zhuǎn)矩分量,當(dāng)KA高于200后,阻尼轉(zhuǎn)矩分量有所提升,但仍然是負(fù)阻尼;對(duì)于雙勵(lì)磁作用的MGP系統(tǒng)(KAM=KAG=KA),阻尼轉(zhuǎn)矩分量隨KA的整體變化趨勢(shì)與傳統(tǒng)機(jī)組是一致的,都是先降低再升高。但是,在相同的KA下,MGP的阻尼轉(zhuǎn)矩分量明顯高于傳統(tǒng)機(jī)組,下降的幅度也較低,當(dāng)KA高于200后呈現(xiàn)出正的阻尼轉(zhuǎn)矩分量,與傳統(tǒng)機(jī)組的特性相反。由此可見,MGP的雙勵(lì)磁耦合特性使系統(tǒng)阻尼顯著提升,兼顧勵(lì)磁機(jī)快速性能的同時(shí)引入正阻尼,這對(duì)于系統(tǒng)穩(wěn)定是有利的。

    為了進(jìn)一步驗(yàn)證雙勵(lì)磁耦合特性,撤去同步電動(dòng)機(jī)的AVR(KAM=0,KAG=KA),這時(shí)阻尼轉(zhuǎn)矩分量不僅在變化趨勢(shì)上仍然與傳統(tǒng)機(jī)組保持一致,KA較高時(shí)的負(fù)阻尼特性也呈現(xiàn)出來(lái),與傳統(tǒng)機(jī)組更加相似,這是因?yàn)榇藭r(shí)的MGP只有發(fā)電機(jī)包含AVR,與傳統(tǒng)機(jī)組的勵(lì)磁效果相近。但是,相同的KA下,MGP的阻尼轉(zhuǎn)矩分量仍然高于傳統(tǒng)機(jī)組而低于雙機(jī)勵(lì)磁,這是因?yàn)榇藭r(shí)MGP中的同步電動(dòng)機(jī)沒(méi)有AVR作用,但它的電樞反應(yīng)是始終存在而且是自發(fā)的,電樞反應(yīng)引起磁通變化的作用是引入一個(gè)正的阻尼轉(zhuǎn)矩分量,導(dǎo)致MGP阻尼轉(zhuǎn)矩分量總的效果仍然優(yōu)于傳統(tǒng)機(jī)組。

    從上述分析可以看出,MGP系統(tǒng)的雙勵(lì)磁耦合特性帶來(lái)的阻尼轉(zhuǎn)矩疊加效應(yīng),不僅由參數(shù)可控的勵(lì)磁系統(tǒng)提供,還由雙同步機(jī)各自的電樞反應(yīng)提供,其本質(zhì)都是通過(guò)影響雙機(jī)的勵(lì)磁磁場(chǎng)起作用,這與式(28)和圖3所反映的雙勵(lì)磁耦合機(jī)理是一致的。MGP的雙勵(lì)磁耦合特性再結(jié)合其自身的轉(zhuǎn)動(dòng)慣量,可以為高比例新能源電網(wǎng)提供有效的阻尼和慣性,從而維持穩(wěn)定性。

    4.2 風(fēng)機(jī)不同并網(wǎng)方式下的小干擾穩(wěn)定對(duì)比

    首先根據(jù)式(22)計(jì)算雙饋風(fēng)機(jī)按傳統(tǒng)方式接入無(wú)窮大母線系統(tǒng)的振蕩模態(tài),再根據(jù)式(27)計(jì)算雙饋風(fēng)機(jī)通過(guò)MGP并網(wǎng)的振蕩模態(tài)。選取傳統(tǒng)并網(wǎng)方式和KA=100時(shí)的振蕩模態(tài)作為比較基準(zhǔn),對(duì)應(yīng)的特征值見附錄A,風(fēng)機(jī)按傳統(tǒng)方式并網(wǎng)時(shí)有三種振蕩模態(tài),其中λ1,2與風(fēng)機(jī)軸系參數(shù)、轉(zhuǎn)差率和暫態(tài)電抗后電勢(shì)d軸分量強(qiáng)相關(guān),λ3,4和λ7,8與暫態(tài)電抗后電勢(shì)和換流器控制參數(shù)強(qiáng)相關(guān);風(fēng)機(jī)通過(guò)MGP并網(wǎng)后,系統(tǒng)增加了一個(gè)與MGP轉(zhuǎn)速?gòu)?qiáng)相關(guān)的特征值λMGP,同時(shí)由2.3節(jié)的推導(dǎo)可知,這時(shí)風(fēng)機(jī)原有的三種振蕩模態(tài)不僅和自身參數(shù)相關(guān),也和MGP狀態(tài)變量,特別是勵(lì)磁磁場(chǎng)相關(guān)聯(lián)。

    結(jié)合實(shí)際運(yùn)行參數(shù),設(shè)置雙勵(lì)磁增益KA在100到400區(qū)間內(nèi)調(diào)節(jié),上述振蕩模態(tài)變化見圖5。

    圖5 振蕩模態(tài)隨雙勵(lì)磁增益的變化Fig.5 Change of oscillation modes with different gains of dual excitation

    從圖5可以看出:λ1,2在數(shù)值上幾乎不變,阻尼比只在KA較大時(shí)有微小提升,這是由于在風(fēng)機(jī)模型(22)中,軸系參數(shù)和轉(zhuǎn)差率本來(lái)就和輸入變量idg和iqg無(wú)關(guān),即矩陣BDFIG中對(duì)應(yīng)元素為0,因此和MGP狀態(tài)方程聯(lián)立后,MGP的狀態(tài)變量幾乎不參與該振蕩模態(tài);而根據(jù)式(22)、式(24)和附錄D,ΔEd′,ΔEq′,Δx1,Δx3均通過(guò)idg和iqg與MGP的勵(lì)磁磁鏈建立聯(lián)系,因此在圖5中,隨著雙勵(lì)磁增益的升高,λ3,4和λ7,8逐漸向遠(yuǎn)離虛軸方向移動(dòng),對(duì)應(yīng)的阻尼比也逐漸提高,表明在雙勵(lì)磁特性作用下,這兩種風(fēng)機(jī)主導(dǎo)模態(tài)的阻尼比和穩(wěn)定域得到有效改善;對(duì)于MGP引入的振蕩模態(tài),阻尼比明顯提升,對(duì)應(yīng)的特征值λMGP向坐標(biāo)平面左側(cè)移動(dòng),這是MGP的雙勵(lì)磁特性直接作用在輸出電磁轉(zhuǎn)矩上的結(jié)果,與圖4(b)中的得到的規(guī)律相符。

    5 結(jié)語(yǔ)

    本文針對(duì)高比例新能源電網(wǎng)阻尼和慣性的缺失問(wèn)題,提出了新能源同步機(jī)用于風(fēng)電并網(wǎng)的方案并研究其小干擾穩(wěn)定性,詳細(xì)建立了能夠反映MGP系統(tǒng)雙勵(lì)磁耦合特性的小干擾模型,揭示雙勵(lì)磁耦合特性對(duì)阻尼的提升作用,并通過(guò)算例進(jìn)行驗(yàn)證。本文的研究有助于進(jìn)一步理解MGP系統(tǒng)在新能源并網(wǎng)中的小干擾建模方法與運(yùn)行特性,但仍然有諸多工作有待深入開展,如完善控制策略、應(yīng)對(duì)新能源的波動(dòng)性、暫態(tài)穩(wěn)定和開展動(dòng)態(tài)模擬實(shí)驗(yàn)等,是下一步的研究重點(diǎn)。

    附錄見本刊網(wǎng)絡(luò)版(http://www.aeps-info.com/aeps/ch/index.aspx)。

    [1] KUNDUR P.Power system stability and control[M].New York:McGraw-hill,1994.

    [2] BECK H P,HESSE R.Virtual synchronous machine[C]// 9th International Conference on Electrical Power Quality and Utilization,October 9-11,2007,Barcelona,Spain:1-6.

    [3] VISSCHER K,de HAAN S W H.Virtual synchronous machines (VSG’s) for frequency stabilisation in future grids with a significant share of decentralized generation[C]// CIRED Seminar SmartGrids for Distribution,June 23-24,2008,Frankfurt,Germany:1-4.

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

    [5] 鄭天文,陳來(lái)軍,陳天一,等.虛擬同步發(fā)電機(jī)技術(shù)及展望[J].電力系統(tǒng)自動(dòng)化,2015,39(21):165-175.DOI:10.7500/AEPS20150508006.

    ZHENG Tianwen,CHEN Laijun,CHEN Tianyi,et al.Review and prospect of virtual synchronous generator technologies[J].Automation of Electric Power Systems,2015,39(21):165-175.DOI:10.7500/AEPS20150508006.

    [6] 呂志鵬,盛萬(wàn)興,劉海濤,等.虛擬同步機(jī)技術(shù)在電力系統(tǒng)中的應(yīng)用與挑戰(zhàn)[J].中國(guó)電機(jī)工程學(xué)報(bào),2017,37(2):349-359.

    Lü Zhipeng,SHENG Wanxing,LIU Haitao,et al.Application and challenge of virtual synchronous machine technology in power system[J].Proceedings of the CSEE,2017,37(2):349-359.

    [7] 劉彬彬,楊健維,廖凱,等.基于轉(zhuǎn)子動(dòng)能控制的雙饋風(fēng)電機(jī)組頻率控制改進(jìn)方案[J].電力系統(tǒng)自動(dòng)化,2016,40(16):17-22.DOI:10.7500/AEPS20150930009.

    LIU Binbin,YANG Jianwei,LIAO Kai,et al.Improved frequency control strategy for DFIG based wind turbines based on rotor kinetic energy control[J].Automation of Electric Power Systems,2016,40(16):17-22.DOI:10.7500/AEPS20150930009.

    [8] ZHANG Y,GEVORGIAN V,ELA E,et al.Role of wind power in primary frequency response of an interconnection[R].Golden,CO,USA:National Renewable Energy Laboratory (NREL),2013.

    [9] 丁磊,尹善耀,王同曉,等.考慮慣性調(diào)頻的雙饋風(fēng)電機(jī)組主動(dòng)轉(zhuǎn)速保護(hù)控制策略[J].電力系統(tǒng)自動(dòng)化,2015,39(24):29-34.DOI:10.7500/AEPS20150529018.

    DING Lei,YIN Shanyao,WANG Tongxiao,et al.Active rotor speed protection strategy for DFIG based wind turbines with inertia control[J].Automation of Electric Power Systems,2015,39(24):29-34.DOI:10.7500/AEPS20150529018.

    [10] WIGET R,ANDERSSON G,ANDREASSON M,et al.Dynamic simulation of a combined AC and MTDC grid with decentralized controllers to share primary frequency control reserves[C]// IEEE PowerTech,June 29-July 2,2015,Eindhoven,Netherlands:1-6.

    [11] 董曉亮,李江,侯金鳴,等.基于雙饋風(fēng)機(jī)轉(zhuǎn)子側(cè)變流器的次同步諧振抑制方法[J].電力系統(tǒng)自動(dòng)化,2016,40(8):92-97.DOI:10.7500/AEPS20150129006.

    DONG Xiaoliang,LI Jiang,HOU Jinming,et al.Supression method of SSR based on rotor side converter of DFIG[J].Automation of Electric Power Systems,2016,40(8):92-97.DOI:10.7500/AEPS20150129006.

    [12] 王毅,張祥宇,李和明,等.永磁直驅(qū)風(fēng)電機(jī)組對(duì)系統(tǒng)功率振蕩的阻尼控制[J].電工技術(shù)學(xué)報(bào),2012,27(12):162-171.

    WANG Yi,ZHANG Xiangyu,LI Heming,et al.Damping control of PMSG-based wind turbines for power system oscillations[J].Transactions of China Electrotechnical Society,2012,27(12):162-171.

    [13] XI X,GENG H,YANG G.Small signal stability of weak power system integrated with inertia tuned large scale wind farm[C]// IEEE Innovative Smart Grid Technologies-Asia,May 20-23,2014,Kuala Lumpur,Malaysia:514-518.

    [14] 馬靜,李益楠,邱揚(yáng),等.雙饋風(fēng)電機(jī)組虛擬慣量控制對(duì)系統(tǒng)小干擾穩(wěn)定性的影響[J].電力系統(tǒng)自動(dòng)化,2016,40(16):1-7.DOI:10.7500/AEPS20150827007.

    MA Jing,LI Yinan,QIU Yang,et al.Impact of virtual inertia control of DFIG wind turbines on system small-signal stability[J].Automation of Electric Power Systems,2016,40(16):1-7.DOI:10.7500/AEPS20150827007.

    [15] 陳潤(rùn)澤,吳文傳,孫宏斌,等.雙饋風(fēng)電機(jī)組慣量控制對(duì)系統(tǒng)小干擾穩(wěn)定的影響[J].電力系統(tǒng)自動(dòng)化,2014,38(23):6-12.DOI:10.7500/AEPS20131213001.

    CHEN Runze,WU Wenchuan,SUN Hongbin,et al.Impact of inertia control of DFIG wind turbines on small-signal stability[J].Automation of Electric Power Systems,2014,38(23):6-12.DOI:10.7500/AEPS20131213001.

    [16] 謝小榮,劉華坤,賀靜波,等.直驅(qū)風(fēng)機(jī)風(fēng)電場(chǎng)與交流電網(wǎng)相互作用引發(fā)次同步振蕩的機(jī)理與特性分析[J].中國(guó)電機(jī)工程學(xué)報(bào),2016,36(9):2366-2372.

    XIE Xiaorong,LIU Huakun,HE Jingbo,et al.Mechanism and characteristics of subsynchronous oscillation caused by the interaction between full-converter wind turbines and AC systems[J].Proceedings of the CSEE,2016,36(9):2366-2372.

    [17] 王清,畢天姝,薛安成.基于相似性的雙饋風(fēng)機(jī)小干擾穩(wěn)定性分析及鎖相環(huán)的影響[J].電力系統(tǒng)自動(dòng)化,2016,40(6):30-34.DOI:10.7500/AEPS20150602006.

    WANG Qing,BI Tianshu,XUE Ancheng.Similarity-based small-signal stability analysis of DFIG-based wind turbines and influence of phase-locked loop[J].Automation of Electric Power Systems,2016,40(6):30-34.DOI:10.7500/AEPS20150602006.

    [18] 和萍,文福拴,薛禹勝,等.不同類型風(fēng)電機(jī)組對(duì)小干擾和暫態(tài)穩(wěn)定性的影響[J].電力系統(tǒng)自動(dòng)化,2013,37(17):23-29.

    HE Ping,WEN Fushuan,XUE Yusheng,et al.Impacts of different wind power generators on power system small signal and transient stability[J].Automation of Electric Power Systems,2013,37(17):23-29.

    [19] MASOOD N A,YAN R,SAHA T K,et al.Post-retirement utilisation of synchronous generators to enhance security performances in a wind dominated power system[J].IET Generation Transmission &Distribution,2016,10(13):3314-3321.

    [20] NAHID A M,YAN R,SAHA T K,et al.Frequency response and its enhancement using synchronous condensers in presence of high wind penetration[C]// IEEE Power &Energy Society General Meeting,July 26-30,2015,Denver,USA:1-5.

    [21] TAYLOR J B.Some features affecting the parallel operation of synchronous motor-generator sets[J].Proceedings of the American Institute of Electrical Engineers,1906,25(3):121-144.

    [22] WANG X,CHAI J,WANG Z.Research on parallel operation characteristics of variable frequency synchronous-motor-generator (SMG)[C]// Proceedings of the Eighth International Conference on Electrical Machines and Systems,September 27-29,2005,Nanjing,China:72-75.

    [23] SAE-KOK W,YOKOYAMA A,OGAWA S,et al.Nonlinear excitation control for rotary type frequency converter using two sets of adjustable speed generators/motors[C]// IEEE Transmission and Distribution Conference and Exhibition 2002:Asia Pacific,October 6-10,2002,Yokohama,Japan:195-200.

    [24] SAE-KOK W,YOKOYAMA A,VERMA S C,et al.Excitation control system design of rotary type frequency converter for performance improvement of power system dynamics[J].IEEE Trans on Energy Conversion,2006,21(1):210-220.

    [25] CRISPO R,BREKKEN T K A.A motor-generator and supercapacitor based system for microgrid frequency stabilization[C]// 1st Conference on IEEE Technologies for Sustainability,August 1-2,2013,Portland,Oregon:162-166.

    [26] REDDY Y J,KUMAR Y V P,RAJU K P,et al.Retrofitted hybrid power system design with renewable energy sources for buildings[J].IEEE Trans on Smart Grid,2012,3(4):2174-2187.

    [27] 湯蘊(yùn)璆,史乃,沈文豹.電機(jī)理論與運(yùn)行[M].北京:水利電力出版社,1983:49-52.

    [28] WEI Siming,ZHOU Yingkun,LI Song,et al.A possible configuration with motor-generator pair for renewable energy integration[J].CSEE Journal of Power and Energy Systems,2017,3(1):93-100.

    [29] 張寶群,張伯明,吳文傳.基于小干擾穩(wěn)定的雙饋感應(yīng)電機(jī)電磁模型降階分析[J].電力系統(tǒng)自動(dòng)化,2010,34(10):51-55.

    ZHANG Baoqun,ZHANG Boming,WU Wenchuan.Small signal based electromagnetic model analysis of double fed induction generator[J].Automation of Electric Power Systems,2010,34(10):51-55.

    [30] 魏巍,劉瑩,丁理杰,等.改進(jìn)的雙饋風(fēng)力發(fā)電系統(tǒng)小干擾模型[J].電網(wǎng)技術(shù),2013,37(10):2904-2911.

    WEI Wei,LIU Ying,DING Lijie,et al.An improved small perturbation model for doubly fed wind power system[J].Power System Technology,2013,37(10):2904-2911.

    Small Signal Modeling and Damping Characteristic of MGP for Wind Power Integration

    WEISiming,HUANGYongzhang

    (State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources (North China Electric Power University),Beijing 102206,China)

    Focusing on low inertia and damping problems brought about by high penetration of renewable energy,this paper proposes the motor-generator pair (MGP) system of renewable energy resource for wind power integration and a detailed modeling method for small-signal stability.First,according to the rotor angle relation of MGP and dynamic process of excitation systems,state equations are established considering the automatic voltage regulator (AVR) and power system stabilizer (PSS),while pointing out that MGP has superimposed effect for damping torque.Then,a comprehensive model consisting of all state variables of doubly-fed induction generator (DFIG) and MGP is developed.On this basis,representation of the torque component based onKcoefficients is derived to reveal that the coupling characteristic of dual excitation is the essence of superimposed effect.Finally,a single machine infinite bus system is tested based on established models.Characteristics of small-signal stability of MGP and traditional generation unit are compared.The oscillation modes dominated by wind turbine are also analyzed with dual excitation coupling taken into account when the wind farm is connected to the grid the traditional way using MGP.Results show that the proposed modeling method can accurately reflect the dual excitation characteristic of MGP to effectively enhance the output component of damping torque and improve oscillation modes.

    high penetration of renewable energy resource;motor-generator pair (MGP);small signal modeling;damping characteristic

    2016-12-16;

    2017-04-02。

    上網(wǎng)日期:2017-06-14。。

    衛(wèi)思明(1990—),男,通信作者,博士研究生,主要研究方向:新能源電力系統(tǒng)穩(wěn)定與控制。E-mail: siming_wei@ncepu.edu.cn

    黃永章(1962—),男,教授,博士生導(dǎo)師,主要研究方向:新能源電力系統(tǒng)、大功率電力電子器件及應(yīng)用、電動(dòng)汽車與電網(wǎng)大數(shù)據(jù)。E-mail:huang_y_z@ncepu.edu.cn

    (編輯孔麗蓓)

    猜你喜歡
    同步機(jī)功角狀態(tài)變量
    一階動(dòng)態(tài)電路零狀態(tài)響應(yīng)公式的通用拓展
    基于TwinCAT3控制系統(tǒng)的YB518型小盒透明紙包裝機(jī)運(yùn)行速度的控制分析
    基于虛擬同步機(jī)的新能源并網(wǎng)智能控制研究
    虛擬調(diào)速器對(duì)VSG暫態(tài)功角穩(wěn)定影響機(jī)理分析
    新能源同步機(jī)與調(diào)相機(jī)在高比例新能源電網(wǎng)中應(yīng)用的比較
    基于嵌套思路的飽和孔隙-裂隙介質(zhì)本構(gòu)理論
    同步機(jī)和異步機(jī)之磁勢(shì)和電勢(shì)的比較
    基于改進(jìn) shapelet 挖掘的風(fēng)電并網(wǎng)系統(tǒng)暫態(tài)功角穩(wěn)定評(píng)估
    能源工程(2019年6期)2019-12-02 01:58:20
    基于FPGA的同步機(jī)軸角編碼器
    基于功角測(cè)量和等面積法則的發(fā)電機(jī)暫態(tài)穩(wěn)定在線判別研究
    嫁个100分男人电影在线观看| 一进一出抽搐gif免费好疼| 国产又爽黄色视频| 中文亚洲av片在线观看爽| 亚洲av电影不卡..在线观看| 一夜夜www| 国产精品一区二区三区四区久久 | 在线av久久热| 国产精品,欧美在线| 搡老岳熟女国产| 亚洲一区二区三区色噜噜| 一级毛片女人18水好多| 99久久综合精品五月天人人| 中国美女看黄片| 亚洲久久久国产精品| 欧美大码av| 亚洲一码二码三码区别大吗| 日韩视频一区二区在线观看| 女同久久另类99精品国产91| 最近最新中文字幕大全电影3 | 成年免费大片在线观看| 人成视频在线观看免费观看| 国产精品日韩av在线免费观看| 久久久久久人人人人人| 中文在线观看免费www的网站 | 久久香蕉国产精品| 日韩成人在线观看一区二区三区| 色婷婷久久久亚洲欧美| 九色国产91popny在线| bbb黄色大片| 99国产综合亚洲精品| 少妇被粗大的猛进出69影院| 久久精品人妻少妇| 天天一区二区日本电影三级| 国产av一区在线观看免费| 午夜精品久久久久久毛片777| 精品国产乱码久久久久久男人| 欧美av亚洲av综合av国产av| 精品少妇一区二区三区视频日本电影| 欧美 亚洲 国产 日韩一| 成人三级做爰电影| 制服丝袜大香蕉在线| 国产不卡一卡二| 久久人人精品亚洲av| 亚洲真实伦在线观看| 啦啦啦韩国在线观看视频| 精品欧美国产一区二区三| 国产亚洲精品久久久久5区| 亚洲人成网站在线播放欧美日韩| 99久久99久久久精品蜜桃| 精品卡一卡二卡四卡免费| 国产亚洲精品综合一区在线观看 | 欧美av亚洲av综合av国产av| 久久久久国产精品人妻aⅴ院| 亚洲精品av麻豆狂野| 亚洲av片天天在线观看| 91九色精品人成在线观看| 日本一本二区三区精品| a在线观看视频网站| 男人操女人黄网站| 国产成人一区二区三区免费视频网站| aaaaa片日本免费| 免费电影在线观看免费观看| 欧美又色又爽又黄视频| 看免费av毛片| 性色av乱码一区二区三区2| 久久99热这里只有精品18| 亚洲精品一区av在线观看| 久久香蕉精品热| 日韩国内少妇激情av| 我的亚洲天堂| 亚洲一区二区三区不卡视频| 看免费av毛片| 最新美女视频免费是黄的| 97超级碰碰碰精品色视频在线观看| 亚洲精品久久国产高清桃花| 18禁美女被吸乳视频| 琪琪午夜伦伦电影理论片6080| 男女床上黄色一级片免费看| 人人妻人人澡欧美一区二区| 欧美一区二区精品小视频在线| 欧美三级亚洲精品| 日本黄色视频三级网站网址| 一边摸一边做爽爽视频免费| 一本精品99久久精品77| 欧美一级毛片孕妇| 久99久视频精品免费| 国产精品美女特级片免费视频播放器 | 热re99久久国产66热| 亚洲国产中文字幕在线视频| 亚洲人成电影免费在线| 久久久久九九精品影院| 国产精华一区二区三区| 男人操女人黄网站| 国产av又大| 亚洲国产欧美一区二区综合| 国产一区二区三区在线臀色熟女| 宅男免费午夜| 美女午夜性视频免费| 国产精品久久久久久亚洲av鲁大| 欧美日韩黄片免| 黄色成人免费大全| 熟妇人妻久久中文字幕3abv| 热re99久久国产66热| 国产在线精品亚洲第一网站| www.精华液| 18禁黄网站禁片午夜丰满| 91国产中文字幕| 男女视频在线观看网站免费 | 国产午夜精品久久久久久| 亚洲精品久久成人aⅴ小说| 波多野结衣av一区二区av| 亚洲国产毛片av蜜桃av| 亚洲片人在线观看| 久久久久久久久久黄片| 亚洲人成网站在线播放欧美日韩| 欧美日本视频| 午夜成年电影在线免费观看| 国产aⅴ精品一区二区三区波| 国产精品一区二区免费欧美| 国产成人av激情在线播放| 亚洲在线自拍视频| 亚洲在线自拍视频| 色综合婷婷激情| 色老头精品视频在线观看| 免费一级毛片在线播放高清视频| av福利片在线| 一级作爱视频免费观看| 精品无人区乱码1区二区| 亚洲成人久久爱视频| 制服丝袜大香蕉在线| 亚洲性夜色夜夜综合| 国产精品av久久久久免费| 亚洲人成77777在线视频| 午夜福利免费观看在线| 久久久久久久午夜电影| 亚洲精品国产精品久久久不卡| 亚洲精品美女久久av网站| 午夜福利免费观看在线| 欧美激情高清一区二区三区| 国产一区二区激情短视频| 亚洲久久久国产精品| 深夜精品福利| 久久精品国产清高在天天线| 一进一出好大好爽视频| 欧美黑人精品巨大| 国产免费av片在线观看野外av| 欧美在线黄色| 久久人人精品亚洲av| 天天躁狠狠躁夜夜躁狠狠躁| 一a级毛片在线观看| 男人舔女人的私密视频| 国产精品综合久久久久久久免费| 在线观看免费视频日本深夜| ponron亚洲| 精华霜和精华液先用哪个| 欧美性猛交╳xxx乱大交人| 中文字幕人妻丝袜一区二区| 国产午夜精品久久久久久| 黄色视频不卡| 中文字幕人妻丝袜一区二区| 叶爱在线成人免费视频播放| 俺也久久电影网| 最新在线观看一区二区三区| 制服人妻中文乱码| 可以在线观看毛片的网站| 国产激情久久老熟女| 首页视频小说图片口味搜索| 欧美另类亚洲清纯唯美| 热99re8久久精品国产| 亚洲精品国产一区二区精华液| 欧美最黄视频在线播放免费| 日韩精品中文字幕看吧| 久热爱精品视频在线9| 身体一侧抽搐| 国产成人系列免费观看| 女生性感内裤真人,穿戴方法视频| 最近最新中文字幕大全免费视频| av超薄肉色丝袜交足视频| 亚洲五月天丁香| 一本久久中文字幕| 欧美乱妇无乱码| 国产极品粉嫩免费观看在线| 国产麻豆成人av免费视频| tocl精华| 国产激情久久老熟女| 国产激情偷乱视频一区二区| 色老头精品视频在线观看| 亚洲欧美精品综合一区二区三区| 757午夜福利合集在线观看| 国产亚洲精品av在线| 亚洲自偷自拍图片 自拍| 一进一出抽搐gif免费好疼| 久热这里只有精品99| 欧美激情极品国产一区二区三区| 国产精品香港三级国产av潘金莲| 香蕉av资源在线| 美女高潮到喷水免费观看| 香蕉丝袜av| 满18在线观看网站| 久久久久九九精品影院| 久久久国产成人精品二区| 男女那种视频在线观看| 亚洲av中文字字幕乱码综合 | 黑人巨大精品欧美一区二区mp4| 妹子高潮喷水视频| 91字幕亚洲| 国产又爽黄色视频| 亚洲第一青青草原| 亚洲av日韩精品久久久久久密| 国产精品av久久久久免费| 亚洲aⅴ乱码一区二区在线播放 | 亚洲熟妇中文字幕五十中出| 国产精品日韩av在线免费观看| 日韩av在线大香蕉| 1024视频免费在线观看| 男女之事视频高清在线观看| 亚洲七黄色美女视频| 中文字幕人妻熟女乱码| 亚洲欧洲精品一区二区精品久久久| 精品卡一卡二卡四卡免费| 日本黄色视频三级网站网址| 亚洲午夜精品一区,二区,三区| 久久久久久人人人人人| 久久精品国产清高在天天线| 成在线人永久免费视频| 两性夫妻黄色片| 国产成人欧美在线观看| 在线av久久热| 国产精品久久电影中文字幕| 免费一级毛片在线播放高清视频| 日本成人三级电影网站| 亚洲国产看品久久| 男女做爰动态图高潮gif福利片| 国产三级在线视频| 久热这里只有精品99| 变态另类丝袜制服| 18禁美女被吸乳视频| 真人一进一出gif抽搐免费| 欧美黄色淫秽网站| 国内久久婷婷六月综合欲色啪| 国产片内射在线| 中文资源天堂在线| 一区二区三区高清视频在线| 日韩成人在线观看一区二区三区| 黄色女人牲交| 国产亚洲欧美在线一区二区| 51午夜福利影视在线观看| 国产亚洲av嫩草精品影院| 国产视频内射| 久久午夜综合久久蜜桃| 国产成人精品久久二区二区免费| 亚洲精品一卡2卡三卡4卡5卡| 成人av一区二区三区在线看| 巨乳人妻的诱惑在线观看| 国产欧美日韩精品亚洲av| 亚洲精品中文字幕一二三四区| 欧美乱妇无乱码| 成年女人毛片免费观看观看9| 久久人妻av系列| 亚洲专区字幕在线| 妹子高潮喷水视频| 成年人黄色毛片网站| 很黄的视频免费| 久久久久免费精品人妻一区二区 | 黄色丝袜av网址大全| av电影中文网址| 熟女电影av网| 巨乳人妻的诱惑在线观看| 99热只有精品国产| 亚洲人成77777在线视频| 中国美女看黄片| 久久久久国产一级毛片高清牌| 啦啦啦韩国在线观看视频| 动漫黄色视频在线观看| 婷婷六月久久综合丁香| 亚洲欧美一区二区三区黑人| 欧美黑人巨大hd| 日韩一卡2卡3卡4卡2021年| 国产aⅴ精品一区二区三区波| 国产精品野战在线观看| 久久久久久免费高清国产稀缺| 久久国产精品男人的天堂亚洲| 久久久久久大精品| 别揉我奶头~嗯~啊~动态视频| 国产av一区在线观看免费| 精品福利观看| 波多野结衣av一区二区av| 国产一区二区在线av高清观看| 久久性视频一级片| 亚洲天堂国产精品一区在线| 国产乱人伦免费视频| 18禁黄网站禁片免费观看直播| 成在线人永久免费视频| 男人操女人黄网站| xxx96com| 中文亚洲av片在线观看爽| 国产av一区二区精品久久| 满18在线观看网站| 成年女人毛片免费观看观看9| 看免费av毛片| 一级a爱视频在线免费观看| 人人妻人人澡欧美一区二区| 十八禁网站免费在线| 夜夜爽天天搞| 国产视频内射| 可以在线观看毛片的网站| 又紧又爽又黄一区二区| 在线视频色国产色| 美女高潮到喷水免费观看| 久久久久九九精品影院| 婷婷丁香在线五月| 色播在线永久视频| 欧美三级亚洲精品| 久久精品国产亚洲av香蕉五月| 亚洲九九香蕉| 国产国语露脸激情在线看| 国产在线精品亚洲第一网站| 欧美不卡视频在线免费观看 | 国产黄片美女视频| 嫩草影视91久久| 深夜精品福利| 久久午夜亚洲精品久久| 精品午夜福利视频在线观看一区| 色综合婷婷激情| 伊人久久大香线蕉亚洲五| 90打野战视频偷拍视频| 亚洲人成伊人成综合网2020| 啦啦啦韩国在线观看视频| 国产亚洲av高清不卡| 国产不卡一卡二| 亚洲狠狠婷婷综合久久图片| e午夜精品久久久久久久| 韩国av一区二区三区四区| 亚洲国产精品sss在线观看| а√天堂www在线а√下载| 国产极品粉嫩免费观看在线| 淫妇啪啪啪对白视频| 欧美日韩黄片免| 久久精品aⅴ一区二区三区四区| 最近最新免费中文字幕在线| 亚洲成国产人片在线观看| 欧美性长视频在线观看| 国内揄拍国产精品人妻在线 | 在线看三级毛片| 日日摸夜夜添夜夜添小说| 国产成人欧美| 天堂动漫精品| 成年人黄色毛片网站| 久久青草综合色| 成人av一区二区三区在线看| 男女床上黄色一级片免费看| 欧美色视频一区免费| 50天的宝宝边吃奶边哭怎么回事| 久久久久久九九精品二区国产 | 国产高清视频在线播放一区| 天堂√8在线中文| 1024手机看黄色片| 久久 成人 亚洲| 亚洲国产毛片av蜜桃av| 夜夜爽天天搞| 成人18禁在线播放| 曰老女人黄片| 1024香蕉在线观看| 制服诱惑二区| 国产激情偷乱视频一区二区| 村上凉子中文字幕在线| 精品免费久久久久久久清纯| 欧美日韩亚洲国产一区二区在线观看| 一级毛片精品| 99久久久亚洲精品蜜臀av| 亚洲av成人av| 18禁黄网站禁片午夜丰满| 搞女人的毛片| 久久狼人影院| 亚洲黑人精品在线| 久久伊人香网站| 亚洲成人免费电影在线观看| 国产精品二区激情视频| 欧美日韩一级在线毛片| 亚洲午夜理论影院| 国产在线精品亚洲第一网站| 午夜视频精品福利| 亚洲av电影在线进入| 很黄的视频免费| 成人永久免费在线观看视频| 国内揄拍国产精品人妻在线 | 女警被强在线播放| 日日爽夜夜爽网站| 好男人电影高清在线观看| 久久久久久大精品| 国产亚洲欧美98| 日本成人三级电影网站| 免费在线观看完整版高清| 99久久国产精品久久久| 看片在线看免费视频| 成人国产一区最新在线观看| or卡值多少钱| 首页视频小说图片口味搜索| 日本熟妇午夜| 久久婷婷人人爽人人干人人爱| 18禁裸乳无遮挡免费网站照片 | 又黄又爽又免费观看的视频| 国产精品免费一区二区三区在线| 国产精华一区二区三区| 国内精品久久久久精免费| 亚洲真实伦在线观看| 亚洲av熟女| 男女那种视频在线观看| av天堂在线播放| 久久中文字幕一级| 性欧美人与动物交配| 欧美绝顶高潮抽搐喷水| 亚洲国产精品999在线| a在线观看视频网站| 国产亚洲精品综合一区在线观看 | 亚洲自拍偷在线| 国产精品 国内视频| 国产高清videossex| 一区二区三区高清视频在线| 亚洲性夜色夜夜综合| 亚洲国产中文字幕在线视频| 91字幕亚洲| 视频区欧美日本亚洲| 国产亚洲欧美精品永久| 免费在线观看影片大全网站| 精品免费久久久久久久清纯| 国产在线精品亚洲第一网站| 亚洲真实伦在线观看| 欧美黑人巨大hd| 丁香六月欧美| 无遮挡黄片免费观看| 国产精品久久视频播放| 国产激情偷乱视频一区二区| 亚洲自拍偷在线| 午夜精品在线福利| 日韩国内少妇激情av| 久久久久久久久久黄片| 成年女人毛片免费观看观看9| 亚洲精品在线美女| 欧美日韩乱码在线| 国产又黄又爽又无遮挡在线| 国产高清视频在线播放一区| 国产黄色小视频在线观看| 一级毛片女人18水好多| 午夜福利欧美成人| 无人区码免费观看不卡| 日本成人三级电影网站| 国产一区二区三区在线臀色熟女| 亚洲中文日韩欧美视频| 国内少妇人妻偷人精品xxx网站 | 男男h啪啪无遮挡| 一区福利在线观看| 香蕉国产在线看| 少妇熟女aⅴ在线视频| 亚洲熟妇熟女久久| 亚洲第一电影网av| 午夜福利高清视频| av视频在线观看入口| 国产又黄又爽又无遮挡在线| www.自偷自拍.com| av视频在线观看入口| 国产精品电影一区二区三区| 日本 av在线| aaaaa片日本免费| 欧美成人午夜精品| 免费女性裸体啪啪无遮挡网站| 精品久久久久久久久久久久久 | 中文字幕人成人乱码亚洲影| 午夜精品久久久久久毛片777| 色老头精品视频在线观看| 亚洲中文av在线| 亚洲国产中文字幕在线视频| 日韩免费av在线播放| 丝袜人妻中文字幕| 欧美精品啪啪一区二区三区| 国产精品日韩av在线免费观看| 精品久久久久久久久久久久久 | 十八禁网站免费在线| 老司机午夜福利在线观看视频| 亚洲成人免费电影在线观看| 国产精品精品国产色婷婷| 久久午夜综合久久蜜桃| 99国产极品粉嫩在线观看| www.www免费av| 成人永久免费在线观看视频| 岛国视频午夜一区免费看| 久久精品国产99精品国产亚洲性色| 久久久久久免费高清国产稀缺| 中文字幕另类日韩欧美亚洲嫩草| 国产高清videossex| 色在线成人网| 99在线视频只有这里精品首页| av免费在线观看网站| 欧美性猛交╳xxx乱大交人| 中文字幕久久专区| 一区二区日韩欧美中文字幕| 久99久视频精品免费| 欧美不卡视频在线免费观看 | 一区二区日韩欧美中文字幕| 亚洲欧美激情综合另类| 悠悠久久av| 男人的好看免费观看在线视频 | 自线自在国产av| 老鸭窝网址在线观看| 亚洲国产日韩欧美精品在线观看 | 国产午夜精品久久久久久| 美女扒开内裤让男人捅视频| 91字幕亚洲| 一边摸一边抽搐一进一小说| 午夜日韩欧美国产| 男女床上黄色一级片免费看| 可以免费在线观看a视频的电影网站| 男人舔奶头视频| 99精品久久久久人妻精品| 好男人电影高清在线观看| 精品欧美国产一区二区三| 变态另类成人亚洲欧美熟女| 欧美日韩中文字幕国产精品一区二区三区| 国产aⅴ精品一区二区三区波| 少妇被粗大的猛进出69影院| 日韩一卡2卡3卡4卡2021年| 女同久久另类99精品国产91| 俄罗斯特黄特色一大片| 看免费av毛片| 久久中文看片网| √禁漫天堂资源中文www| 亚洲国产精品sss在线观看| 热99re8久久精品国产| 999精品在线视频| 亚洲精华国产精华精| 亚洲va日本ⅴa欧美va伊人久久| 露出奶头的视频| 人妻丰满熟妇av一区二区三区| 国产免费av片在线观看野外av| 国产精品亚洲av一区麻豆| 欧美亚洲日本最大视频资源| 午夜福利在线观看吧| 欧美 亚洲 国产 日韩一| 国产三级黄色录像| 正在播放国产对白刺激| 国产亚洲精品第一综合不卡| 久久久国产成人免费| 亚洲欧美精品综合一区二区三区| 免费观看人在逋| 亚洲国产欧美网| 国产精品九九99| 中文在线观看免费www的网站 | 一区二区三区高清视频在线| www.www免费av| 欧美日韩精品网址| 一夜夜www| 一本综合久久免费| 国产一区二区在线av高清观看| 国产aⅴ精品一区二区三区波| 男人舔女人下体高潮全视频| 最新在线观看一区二区三区| 18禁观看日本| 久久香蕉激情| 亚洲精品中文字幕在线视频| 国产成人影院久久av| 国产精品98久久久久久宅男小说| 精品免费久久久久久久清纯| 久久精品国产综合久久久| 琪琪午夜伦伦电影理论片6080| 精品国产乱子伦一区二区三区| 色综合亚洲欧美另类图片| 国产激情久久老熟女| www日本黄色视频网| 久久精品国产99精品国产亚洲性色| 日韩欧美国产一区二区入口| 免费在线观看黄色视频的| 成人av一区二区三区在线看| 精品欧美国产一区二区三| 亚洲中文字幕日韩| 又大又爽又粗| 婷婷精品国产亚洲av| 夜夜爽天天搞| 成年版毛片免费区| 最近最新中文字幕大全免费视频| 老司机深夜福利视频在线观看| 久久人妻福利社区极品人妻图片| 窝窝影院91人妻| 成人免费观看视频高清| 精品卡一卡二卡四卡免费| 久久久久久久午夜电影| 久久国产亚洲av麻豆专区| 亚洲无线在线观看| 精品高清国产在线一区| 成人免费观看视频高清| 国产精品电影一区二区三区| 国产精品九九99| 99国产精品99久久久久| 久久久久久大精品| 久久精品国产综合久久久| 极品教师在线免费播放| 香蕉久久夜色| 成人午夜高清在线视频 | 黑丝袜美女国产一区| 无限看片的www在线观看| 51午夜福利影视在线观看| 久久久久国产精品人妻aⅴ院| 满18在线观看网站| 黄色成人免费大全| 在线观看免费日韩欧美大片| 不卡一级毛片| 热99re8久久精品国产| 国产野战对白在线观看| 欧美中文日本在线观看视频| 在线观看一区二区三区| 精品久久久久久成人av| 色综合亚洲欧美另类图片| 精品第一国产精品| 在线天堂中文资源库| 久久久久久久久中文|