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

    電網(wǎng)故障下全功率風(fēng)電場(chǎng)內(nèi)部機(jī)組的暫態(tài)穩(wěn)定性分析

    2022-03-16 00:30:08柳飛揚(yáng)
    關(guān)鍵詞:變流器暫態(tài)風(fēng)電場(chǎng)

    柳飛揚(yáng),曾 平,李 征

    電網(wǎng)故障下全功率風(fēng)電場(chǎng)內(nèi)部機(jī)組的暫態(tài)穩(wěn)定性分析

    柳飛揚(yáng),曾 平,李 征

    (東華大學(xué)信息科學(xué)與技術(shù)學(xué)院,上海 201620)

    全功率風(fēng)電場(chǎng)中機(jī)組并網(wǎng)運(yùn)行的穩(wěn)定性往往與網(wǎng)側(cè)變流器的控制和并網(wǎng)阻抗有關(guān),因而風(fēng)電場(chǎng)中機(jī)組的失穩(wěn)并不是同時(shí)發(fā)生的。為揭示風(fēng)電場(chǎng)中機(jī)組的位置與其暫態(tài)穩(wěn)定性的關(guān)系以及機(jī)組間穩(wěn)定性的相互影響,建立了基于鎖相環(huán)的多風(fēng)電機(jī)組非線性降階模型?;诘让娣e法的思想,分析了不同運(yùn)行條件下風(fēng)電場(chǎng)內(nèi)部不同機(jī)組的暫態(tài)穩(wěn)定裕度以及場(chǎng)內(nèi)線路阻抗對(duì)穩(wěn)定性的影響。通過(guò)Matlab/Simulink時(shí)域仿真,進(jìn)行了定量分析和驗(yàn)證。結(jié)果表明,風(fēng)電場(chǎng)內(nèi)饋線阻抗將減少輸出電流流經(jīng)該阻抗所聯(lián)接機(jī)組的暫態(tài)穩(wěn)定域,而對(duì)輸出電流不流經(jīng)該阻抗的其他機(jī)組的暫態(tài)穩(wěn)定性影響很小。阻抗上流過(guò)的電流越大,其對(duì)穩(wěn)定裕度的縮減也越明顯,因而,隨著鏈接機(jī)組臺(tái)數(shù)增多,不能簡(jiǎn)單忽略場(chǎng)內(nèi)阻抗對(duì)穩(wěn)定性的影響。此外,在相同條件下,有功出力大的機(jī)組將更易失穩(wěn)。

    電網(wǎng)故障;全功率風(fēng)電場(chǎng);暫態(tài)穩(wěn)定性;鎖相環(huán);等面積法;穩(wěn)定裕度;場(chǎng)內(nèi)阻抗

    0 引言

    風(fēng)能作為一種清潔高效且儲(chǔ)量豐富的可再生能源,在能源危機(jī)和環(huán)境保護(hù)問(wèn)題日益嚴(yán)峻的當(dāng)下,受到了世界各國(guó)的廣泛關(guān)注[1-4]。大量分布式能源通過(guò)并網(wǎng)變流器接入電網(wǎng),降低了系統(tǒng)的慣量,由此導(dǎo)致的并網(wǎng)穩(wěn)定性問(wèn)題時(shí)有發(fā)生[5-8]。其中,諸如文獻(xiàn)[9]中發(fā)生的風(fēng)電場(chǎng)并網(wǎng)的暫態(tài)穩(wěn)定性問(wèn)題引起了相關(guān)從業(yè)者的廣泛重視,在風(fēng)電滲透率不斷增加的背景下,風(fēng)電并網(wǎng)的暫態(tài)穩(wěn)定性研究具有極其重要的意義。

    全功率變速風(fēng)電機(jī)組變流器多采用矢量控制,鎖相環(huán)(PLL)能夠估計(jì)電網(wǎng)端電壓的頻率和角度,為矢量控制提供基準(zhǔn)。然而,基于矢量控制下的網(wǎng)側(cè)電壓源型變流器(VSC)在連接到較弱電網(wǎng)或電網(wǎng)嚴(yán)重故障時(shí)很可能會(huì)產(chǎn)生同步穩(wěn)定性問(wèn)題[10-14]。相關(guān)研究表明,PLL對(duì)于網(wǎng)側(cè)VSC的小信號(hào)穩(wěn)定性極其重要。文獻(xiàn)[11-12]采用特征值法研究了PLL增益和帶寬對(duì)VSC穩(wěn)定性的影響。文獻(xiàn)[15]采用復(fù)轉(zhuǎn)矩系數(shù)法研究了PLL和交流控制系統(tǒng)的阻尼特性對(duì)VSC穩(wěn)定性的影響。文獻(xiàn)[16]采用阻抗分析法研究了弱電網(wǎng)下PLL和電流環(huán)控制參數(shù)對(duì)VSC穩(wěn)定性的影響。不僅如此,在實(shí)際電力系統(tǒng)中,電網(wǎng)故障會(huì)對(duì)PLL產(chǎn)生不同程度的影響,這取決于電網(wǎng)故障的嚴(yán)重程度、故障位置、線路配置以及系統(tǒng)的控制和電氣參數(shù)等[17-19]。文獻(xiàn)[20]指出,電網(wǎng)故障下PLL和電流控制環(huán)的非線性相互作用會(huì)對(duì)系統(tǒng)的暫態(tài)穩(wěn)定性產(chǎn)生很大影響。越來(lái)越多的研究表明,鎖相環(huán)的動(dòng)態(tài)特性會(huì)與機(jī)組有功電流及機(jī)組到電網(wǎng)之間的阻抗特性耦合,在電網(wǎng)故障期間,由鎖相環(huán)主導(dǎo)的非線性動(dòng)力學(xué)行為很有可能導(dǎo)致系統(tǒng)暫態(tài)失穩(wěn),這一點(diǎn)在較弱電網(wǎng)條件下尤為明顯。因此,若要研究風(fēng)電場(chǎng)的暫態(tài)穩(wěn)定性,則必須考慮電網(wǎng)故障期間PLL與網(wǎng)側(cè)VSC的相互作用,這意味著必須對(duì)PLL進(jìn)行合理建模。

    針對(duì)電網(wǎng)故障下PLL的動(dòng)態(tài)行為已有一些研究。文獻(xiàn)[21]研究了電網(wǎng)故障期間風(fēng)電場(chǎng)注入電流對(duì)風(fēng)電場(chǎng)暫態(tài)同步穩(wěn)定性的影響,得出了系統(tǒng)穩(wěn)定運(yùn)行的電流邊界條件,但該判據(jù)是基于穩(wěn)態(tài)下建立的,沒(méi)有考慮鎖相環(huán)的動(dòng)態(tài)行為對(duì)系統(tǒng)暫態(tài)穩(wěn)定性的影響,因此該穩(wěn)定性判據(jù)是不充分的。文獻(xiàn)[22]指出了風(fēng)電場(chǎng)和電力系統(tǒng)在弱連接情況下,PLL會(huì)與風(fēng)電場(chǎng)注入電流及線路阻抗等因素相互作用從而導(dǎo)致暫態(tài)同步問(wèn)題,并由此提出了一種電網(wǎng)故障下的電流注入策略。但該研究是基于風(fēng)場(chǎng)聚合模型進(jìn)行的,無(wú)法反映電網(wǎng)故障期間場(chǎng)內(nèi)各機(jī)組的暫態(tài)穩(wěn)定性。文獻(xiàn)[23]采用相圖法對(duì)網(wǎng)側(cè)VSC的暫態(tài)穩(wěn)定性進(jìn)行分析,提出了一種自適應(yīng)PLL以保證暫態(tài)穩(wěn)定性和矢量角的追蹤精度。文獻(xiàn)[24]利用修正的等面積法研究了電網(wǎng)故障期間PLL的動(dòng)態(tài)特性,結(jié)果表明PLL的積分調(diào)節(jié)器不利于系統(tǒng)的再同步,并由此提出一種變結(jié)構(gòu)PLL以保證系統(tǒng)在故障期間的暫態(tài)同步穩(wěn)定性。文獻(xiàn)[25]的研究表明增大PLL阻尼比可以提高系統(tǒng)的暫態(tài)穩(wěn)定性,但尚未對(duì)其進(jìn)行深入的定量分析,阻尼比的設(shè)計(jì)準(zhǔn)則依舊是缺失的。文獻(xiàn)[26-28]研究了PLL非線性動(dòng)力學(xué)對(duì)全功率機(jī)組的暫態(tài)同步穩(wěn)定性的影響,建立了嚴(yán)格的PLL非線性模型,類比同步發(fā)電機(jī)的轉(zhuǎn)動(dòng)方程,基于等面積法揭示了電網(wǎng)故障下機(jī)組暫態(tài)失穩(wěn)的機(jī)理,結(jié)果表明PLL帶寬、短路比、故障持續(xù)時(shí)間對(duì)系統(tǒng)的暫態(tài)穩(wěn)定性有影響。其中文獻(xiàn)[27]基于PLL在電網(wǎng)故障下的動(dòng)態(tài)行為提出一種有功電流、無(wú)功電流時(shí)序控制策略,提高了全功率機(jī)組的頻率穩(wěn)定性。文獻(xiàn)[28]進(jìn)一步研究了阻抗比、有功無(wú)功電流比、短路點(diǎn)相移對(duì)全功率機(jī)組暫態(tài)同步穩(wěn)定性的影響。

    總的來(lái)說(shuō),針對(duì)電網(wǎng)故障期間PLL的非線性動(dòng)力學(xué)行為的機(jī)理研究已初步取得了一些成果,對(duì)注入電流、并網(wǎng)阻抗這些與PLL耦合的因素以及PLL的控制參數(shù)已有了一些定量分析,但風(fēng)電場(chǎng)內(nèi)部機(jī)組的暫態(tài)穩(wěn)定性及場(chǎng)內(nèi)各機(jī)組之間的阻抗對(duì)機(jī)組暫態(tài)穩(wěn)定性的影響研究還比較匱乏。

    本文將從故障分析的角度,建立全功率風(fēng)場(chǎng)的并網(wǎng)等效電路以及場(chǎng)內(nèi)各機(jī)組的PLL非線性降階模型,在考慮場(chǎng)內(nèi)機(jī)組之間的相互作用的前提下著重研究電網(wǎng)嚴(yán)重故障下場(chǎng)內(nèi)任一單臺(tái)機(jī)組的暫態(tài)穩(wěn)定性,基于等面積法思想定量分析場(chǎng)內(nèi)各機(jī)組的暫態(tài)穩(wěn)定裕度,然后定性分析風(fēng)電場(chǎng)中機(jī)組之間阻抗大小對(duì)機(jī)組暫態(tài)穩(wěn)定性的影響。建立了基于Matlab/Simulink平臺(tái)的多機(jī)時(shí)域仿真模型,驗(yàn)證了理論分析的正確性。

    1 全功率風(fēng)電場(chǎng)的暫態(tài)穩(wěn)定性建模

    1.1 全功率風(fēng)電場(chǎng)的并網(wǎng)與控制系統(tǒng)

    圖1 全功率風(fēng)電場(chǎng)的接線示意圖

    圖2 全功率風(fēng)電機(jī)組單機(jī)控制結(jié)構(gòu)圖

    1.2 全功率機(jī)組的簡(jiǎn)化模型及分析

    電網(wǎng)故障下風(fēng)電機(jī)組的完整數(shù)學(xué)模型較為復(fù)雜,往往需要根據(jù)研究的問(wèn)題進(jìn)行簡(jiǎn)化。時(shí)間尺度分析方法能夠?yàn)槟P秃?jiǎn)化提供依據(jù)。在電網(wǎng)故障期間,由于全功率機(jī)組發(fā)電機(jī)的機(jī)械動(dòng)態(tài)特性在時(shí)間尺度上往往遠(yuǎn)大于其變流器控制環(huán)節(jié)[29],機(jī)組的機(jī)械動(dòng)態(tài)特性可以不予考慮;此外,故障期間由于卸荷電路的投入,變流器直流母線電壓能夠保持在工作范圍內(nèi),不會(huì)影響到機(jī)側(cè)變流器的運(yùn)行,因而全功率變流器對(duì)電網(wǎng)故障起到了隔離的作用,在暫態(tài)穩(wěn)定分析過(guò)程中機(jī)側(cè)變流器可以簡(jiǎn)化掉;在網(wǎng)側(cè)變流器中,PLL和電壓外環(huán)響應(yīng)速度往往遠(yuǎn)小于電流內(nèi)環(huán),因此在討論基于PLL動(dòng)態(tài)行為的時(shí)間尺度下的穩(wěn)定性問(wèn)題時(shí),可認(rèn)為電流內(nèi)環(huán)保持恒定,網(wǎng)側(cè)變流器在電網(wǎng)故障期間可視為電流源。此時(shí)PLL在

    電網(wǎng)故障下的動(dòng)態(tài)行為需要較為精確的描述。

    基于上述分析,可以得出電網(wǎng)故障期間全功率風(fēng)電機(jī)組暫態(tài)穩(wěn)定性分析模型的簡(jiǎn)化假設(shè)條件:

    1) 機(jī)組的機(jī)械動(dòng)態(tài)特性時(shí)間尺度較大,在電網(wǎng)故障期間不予考慮。

    2) 電網(wǎng)故障期間由于卸荷電路的投入,假設(shè)直流側(cè)電壓保持穩(wěn)定,機(jī)側(cè)變流器被簡(jiǎn)化。

    3) 電網(wǎng)故障期間網(wǎng)側(cè)變流器電流內(nèi)環(huán)保持恒定,網(wǎng)側(cè)變流器等效為電流源,網(wǎng)側(cè)變流器控制可采用恒定電流基準(zhǔn)。

    4) PLL的動(dòng)態(tài)行為對(duì)機(jī)組的暫態(tài)穩(wěn)定性至關(guān)重要,應(yīng)用非線性模型對(duì)PLL詳細(xì)建模。

    圖3 全功率機(jī)組暫態(tài)穩(wěn)定性分析簡(jiǎn)化模型

    1.3 全功率風(fēng)電場(chǎng)的并網(wǎng)等效電路

    對(duì)于圖1中的全功率風(fēng)場(chǎng)接線,可用圖4所示的三臺(tái)機(jī)組接線表示其機(jī)組之間的連接關(guān)系。圖4中1、2為風(fēng)場(chǎng)內(nèi)機(jī)組之間的等效線路阻抗,line為機(jī)組A一次升壓變壓器T1與二次升壓變壓器T2之間的等效線路阻抗。若研究35 kV低壓母線上不同饋線上的第一臺(tái)機(jī)組在電網(wǎng)故障下的暫態(tài)穩(wěn)定性關(guān)系,考慮到這些機(jī)組之間的電氣距離很小,認(rèn)為機(jī)組之間的阻抗為0,即設(shè)1、2為0;若需研究同一饋線上鏈接的多臺(tái)機(jī)組,此時(shí)這些機(jī)組之間的電氣距離較大,機(jī)組之間的阻抗不能忽略,1、2可取機(jī)組之間線路的阻抗,1、2的不同取值可以體現(xiàn)出不同機(jī)組之間的電氣距離的大小。

    研究并聯(lián)在低壓母線上的風(fēng)電機(jī)組穩(wěn)定性時(shí),考慮到風(fēng)場(chǎng)中機(jī)組之間及變壓器T1與T2之間的電氣距離一般不會(huì)太大,1、2將遠(yuǎn)小于變壓器阻抗,可在機(jī)理分析中,忽略它們?cè)陔娋W(wǎng)故障下對(duì)機(jī)組的暫態(tài)穩(wěn)定性影響,本文的原理分析著重討論這種情況,結(jié)論可推廣到任意情況。機(jī)組之間的阻抗對(duì)各機(jī)組的暫態(tài)穩(wěn)定性影響將在仿真分析中作更全面的展現(xiàn)。

    圖4 全功率風(fēng)電場(chǎng)并網(wǎng)結(jié)構(gòu)圖

    當(dāng)電網(wǎng)在點(diǎn)PCC處發(fā)生故障時(shí),網(wǎng)側(cè)變流器采用恒定電流基準(zhǔn)控制,風(fēng)場(chǎng)內(nèi)機(jī)組可以等效為電流源,依據(jù)圖3和圖4,忽略1、2對(duì)機(jī)組的暫態(tài)穩(wěn)定性影響,并將線路阻抗line折算到變壓器T2上,可以得到如圖5所示的全功率風(fēng)電場(chǎng)并網(wǎng)暫態(tài)穩(wěn)定性分析等效電路。

    由于變壓器繞組中的電阻與電感相比很小,為簡(jiǎn)化分析,忽略T1和T2中的電阻,認(rèn)為兩者是近似純感性的。由圖5可得電網(wǎng)故障后PCC點(diǎn)的電壓矢量為

    圖5 全功率風(fēng)電場(chǎng)并網(wǎng)暫態(tài)穩(wěn)定分析等效電路

    Fig. 5 Equivalent circuit for transient stability analysis of grid-connected full-scale wind farm

    1.4 基于PLL的風(fēng)電場(chǎng)暫態(tài)穩(wěn)定性建模

    為了聚焦電網(wǎng)故障期間PLL的動(dòng)態(tài)行為,可假定故障期間電流內(nèi)環(huán)保持穩(wěn)定,則網(wǎng)側(cè)VSC輸出電流值等于電流基準(zhǔn)值,而通常網(wǎng)側(cè)VSC采用恒定的電流基準(zhǔn)值。PLL的控制結(jié)構(gòu)圖如圖6所示,下文的分析將在PLL參考坐標(biāo)系下進(jìn)行,相關(guān)電壓電流變量用“pll”上標(biāo)表示。

    圖6 PLL的控制結(jié)構(gòu)圖

    由圖6可以得到

    當(dāng)全功率風(fēng)電場(chǎng)中機(jī)組A、B、C出力不相同,機(jī)組都以單位功率因數(shù)運(yùn)行,故障前后機(jī)組A、B、C的電壓相量圖如圖7、圖8所示,其中T1為機(jī)組A、B、C經(jīng)變壓器T1后共同連接點(diǎn)處電壓,LT1為變壓器T1的電感。

    圖7 故障前機(jī)組A、B、C電壓相量圖

    圖8 故障后機(jī)組A、B、C電壓相量圖

    假定在風(fēng)電場(chǎng)機(jī)組間物理距離范圍內(nèi)風(fēng)速的變化較為平緩,此時(shí)機(jī)組有功出力相差不是太大,若故障后PLL不會(huì)失穩(wěn),則故障前后作近似處理有:

    以機(jī)組A為例,在機(jī)組A的PLL參考坐標(biāo)系下,有以下關(guān)系:

    綜合以上分析,以機(jī)組A為例,推導(dǎo)全功率風(fēng)電場(chǎng)中單臺(tái)機(jī)組的PLL暫態(tài)穩(wěn)定性非線性模型。對(duì)于機(jī)組A,有

    式中,pcc為PCC電壓矢量的模值。

    將式(5)代入式(2)中并標(biāo)幺化后可得到

    式(6)中各變量具體含義如式(7)所示。

    由式(7)可以看出,機(jī)組A的非線性動(dòng)力學(xué)行為不僅與機(jī)組A的輸出電流有關(guān),還與場(chǎng)內(nèi)其他機(jī)組(B、C)的輸出電流有關(guān),也即故障期間各機(jī)組之間存在著交互作用,而這一點(diǎn)是機(jī)組單機(jī)模型和風(fēng)場(chǎng)聚合模型所無(wú)法體現(xiàn)的。風(fēng)場(chǎng)中其他機(jī)組的PLL非線性暫態(tài)模型也可類似得到,在此不再贅述。

    2 風(fēng)場(chǎng)內(nèi)單機(jī)的暫態(tài)穩(wěn)定性機(jī)理分析

    在1.4節(jié)中得到的PLL非線性暫態(tài)模型與同步發(fā)電機(jī)的轉(zhuǎn)動(dòng)方程類似,因此,可以考慮用同步發(fā)電機(jī)中功角穩(wěn)定的相關(guān)概念來(lái)解釋全功率風(fēng)電場(chǎng)在電網(wǎng)故障下出現(xiàn)的暫態(tài)失穩(wěn)問(wèn)題。

    暫態(tài)穩(wěn)定性問(wèn)題是指當(dāng)系統(tǒng)遭遇大擾動(dòng)后,能否回到原來(lái)的平衡點(diǎn)或者運(yùn)行于新平衡點(diǎn)的問(wèn)題。下文將基于同步發(fā)電機(jī)中“等面積法”的思想,對(duì)全功率風(fēng)電場(chǎng)中各機(jī)組的暫態(tài)穩(wěn)定性進(jìn)行分析。

    2.1 暫態(tài)穩(wěn)定性機(jī)理分析

    圖9 暫態(tài)穩(wěn)定性分析示意圖

    需要指出的是,盡管此類情況不會(huì)直接導(dǎo)致系統(tǒng)出現(xiàn)暫態(tài)同步穩(wěn)定性問(wèn)題,但故障后系統(tǒng)從原來(lái)的平衡點(diǎn)過(guò)渡到新平衡點(diǎn)往往會(huì)伴隨著劇烈的電壓和電流變化,這將會(huì)損壞系統(tǒng)中的電力電子器件或觸發(fā)器件保護(hù),這是實(shí)際情況中不允許發(fā)生的。因此,此種情況雖沒(méi)有導(dǎo)致系統(tǒng)失步,但也屬于暫態(tài)失穩(wěn)的范疇,這一點(diǎn)將在下文的仿真中得到體現(xiàn)。

    2.2 暫態(tài)穩(wěn)定裕度分析

    2.3 機(jī)組之間阻抗對(duì)機(jī)組暫態(tài)穩(wěn)定性的影響

    如圖4所示,全功率風(fēng)電場(chǎng)內(nèi)機(jī)組之間用電纜相連,其等值電路如圖10所示。

    圖10 電纜的等值電路

    由圖10可知,與忽略場(chǎng)內(nèi)機(jī)組之間的阻抗相比,考慮電纜等值阻抗將會(huì)使得機(jī)組到PCC點(diǎn)的總阻抗增大,同時(shí)會(huì)有電流損耗。

    對(duì)于如圖4所示的風(fēng)電場(chǎng),若考慮場(chǎng)內(nèi)阻抗,則距離場(chǎng)站升壓變壓器T2電氣距離最近的機(jī)組A 由于電流損耗,其暫態(tài)穩(wěn)定裕度將比不考慮場(chǎng)內(nèi)阻抗的情況要大,然而該影響實(shí)際上很小,大多數(shù)時(shí)候無(wú)需著重考慮;而對(duì)于場(chǎng)內(nèi)其他機(jī)組B和C,考慮到電流損耗較小,此時(shí)機(jī)組與到PCC點(diǎn)的阻抗成為主導(dǎo)因素,其暫態(tài)穩(wěn)定裕度將會(huì)減小。

    3 暫態(tài)穩(wěn)定裕度的仿真分析

    本節(jié)以Matlab/Simulink為仿真平臺(tái),建立了如圖4所示包含3臺(tái)6.25 MW永磁同步發(fā)電機(jī)(PMSG)的風(fēng)電場(chǎng)模型,模型中未加入場(chǎng)內(nèi)機(jī)組之間的阻抗,三臺(tái)機(jī)組控制參數(shù)一致,主要仿真參數(shù)由表1給出。

    3.1 不考慮場(chǎng)內(nèi)阻抗的暫態(tài)穩(wěn)定性仿真

    3.1.1三臺(tái)機(jī)組有功電流相同的情況

    A、B、C三臺(tái)機(jī)組均設(shè)置有功電流為0.9 p.u.,無(wú)功電流為0 p.u.的恒定電流基準(zhǔn),PCC處用受控電壓源模擬電網(wǎng)故障,在1 s時(shí)設(shè)置90%無(wú)相移的電壓跌落,故障持續(xù)時(shí)間分別為125.8 ms(機(jī)組A的極限切除時(shí)間)、126 ms、240 ms,機(jī)組A的暫態(tài)穩(wěn)定性仿真圖如圖11—圖13所示。

    表1 主要仿真參數(shù)

    圖11 機(jī)組A的暫態(tài)穩(wěn)定性仿真圖(125.8 ms)

    圖12 機(jī)組A的暫態(tài)穩(wěn)定性仿真圖(126 ms)

    圖13 機(jī)組A的暫態(tài)穩(wěn)定性仿真圖(240 ms)

    由圖13可以看出,當(dāng)故障持續(xù)時(shí)間太長(zhǎng)(240 ms)時(shí),故障切除后PLL將失穩(wěn),系統(tǒng)與電網(wǎng)失去同步,機(jī)組的有功/無(wú)功電流將發(fā)生振蕩,機(jī)組發(fā)生了暫態(tài)同步穩(wěn)定性問(wèn)題。

    3.1.2 三臺(tái)機(jī)組有功電流不同的情況

    A、B、C三臺(tái)機(jī)組分別設(shè)置有功電流為0.9 p.u.、0.8 p.u.、0.7 p.u.,無(wú)功電流均為0 p.u.的恒定電流基準(zhǔn),PCC處用受控電壓源模擬電網(wǎng)故障,在1 s時(shí)設(shè)置90%無(wú)相移的電壓跌落,故障持續(xù)時(shí)間分別為134 ms(機(jī)組A極限切除時(shí)間)、142.1 ms(機(jī)組B極限切除時(shí)間)、152.5 ms(機(jī)組C極限切除時(shí)間),三種情況下機(jī)組A、B、C的暫態(tài)穩(wěn)定性仿真圖(暫態(tài)虛擬功角)如圖14—圖16所示。

    圖14 機(jī)組A、B、C的暫態(tài)穩(wěn)定性仿真圖(134 ms)

    圖15 機(jī)組A、B、C的暫態(tài)穩(wěn)定性仿真圖(142.1 ms)

    圖16 機(jī)組A、B、C的暫態(tài)穩(wěn)定性仿真圖(152.5 ms)

    3.1.3 仿真結(jié)果分析

    可以看出,本文提出的暫態(tài)穩(wěn)定裕度計(jì)算方法較為保守,所得結(jié)果比實(shí)際值要小,這可能是由以下因素引起的:

    1) 系統(tǒng)本身存在阻尼,而理論分析中沒(méi)有且很難考慮系統(tǒng)自身阻尼的影響,這將導(dǎo)致實(shí)際的暫態(tài)穩(wěn)定裕度比理論計(jì)算值要大。

    2) 理論分析忽略了系統(tǒng)中電阻的影響,所得 結(jié)果偏保守。

    3) 風(fēng)場(chǎng)中機(jī)組出力不同時(shí)采用了近似處理。

    4) 等值電路參數(shù)的計(jì)算采用了近似計(jì)算,因而結(jié)果略有差異。

    3.2 機(jī)組之間阻抗對(duì)暫態(tài)穩(wěn)定性影響仿真驗(yàn)證

    仿真設(shè)置參考3.1.2節(jié),故障持續(xù)時(shí)間分別設(shè)置為134.1 ms、142.1 ms、152.5 ms,分別在忽略與考慮機(jī)組之間阻抗1、2(見(jiàn)圖4)的情況下進(jìn)行仿真。假設(shè)三臺(tái)機(jī)組是在同一條饋線上的不同位置,其中1、2的參數(shù)依據(jù)實(shí)際風(fēng)電場(chǎng)內(nèi)的機(jī)組分布計(jì)算得出,如表2所示,機(jī)組A、B、C的暫態(tài)穩(wěn)定性仿真圖(暫態(tài)虛擬功角)如圖17—圖19所示。

    圖17 機(jī)組A的暫態(tài)穩(wěn)定性仿真圖(134.1 ms)

    表2 阻抗仿真參數(shù)

    由圖17可以看出,在134.1 ms(略大于A機(jī)組的極限切除時(shí)間)的故障持續(xù)時(shí)間下,若忽略機(jī)組之間阻抗,機(jī)組A暫態(tài)失穩(wěn),若考慮機(jī)組之間的阻抗,機(jī)組A暫態(tài)穩(wěn)定。此時(shí),B、C機(jī)組均是穩(wěn)定的??梢?jiàn)考慮饋線阻抗增加了前端機(jī)組的穩(wěn)定性。

    圖18、圖19分別為142.1 ms、152.5 ms的故障持續(xù)時(shí)間下的仿真波形,可見(jiàn),若忽略機(jī)組之間阻抗,則機(jī)組B、C暫態(tài)穩(wěn)定,若考慮機(jī)組之間阻抗,則機(jī)組B、C暫態(tài)失穩(wěn)。

    上述仿真結(jié)果驗(yàn)證了理論分析的正確性,即考慮機(jī)組之間的阻抗與不考慮機(jī)組之間的阻抗相比,機(jī)組A的暫態(tài)穩(wěn)定性提高,機(jī)組B、C的暫態(tài)穩(wěn)定性會(huì)降低。

    依據(jù)表2中阻抗參數(shù),在3.1.2節(jié)的仿真設(shè)置條件下,機(jī)組A、B、C 的仿真結(jié)果如表3所示,其中A、B、C分別為機(jī)組A、B、C的故障極限切除時(shí)間。

    圖18 機(jī)組B的暫態(tài)穩(wěn)定性仿真圖(142.1 ms)

    圖19 機(jī)組C的暫態(tài)穩(wěn)定性仿真圖(152.5 ms)

    表3 考慮阻抗的仿真結(jié)果

    其他條件相同的情況下,增大機(jī)組A、B之間的電氣距離,即增大1,機(jī)組A、B、C的仿真結(jié)果如表4所示。

    表4 增大Z1的仿真結(jié)果

    由表3和表4的仿真結(jié)果可知,增大機(jī)組A、B之間的電氣距離,機(jī)組B、C穩(wěn)定的極限切除時(shí)間減少,因而暫態(tài)穩(wěn)定域顯著降低,而機(jī)組A的暫態(tài)穩(wěn)定性基本不變。

    其他條件相同的情況下,增大機(jī)組B、C之間的電氣距離,即增大2,機(jī)組A、B、C的仿真結(jié)果如表5所示。

    表5 增大Z2的仿真結(jié)果

    對(duì)比表3和表5的仿真結(jié)果可知,增大機(jī)組B、C之間的電氣距離會(huì)降低機(jī)組C的暫態(tài)穩(wěn)定域,而機(jī)組A、B的暫態(tài)穩(wěn)定性基本不變。

    此外,對(duì)比表4和表5的仿真結(jié)果可知,盡管2的增加量大于1的增加量,但機(jī)組C的穩(wěn)定極限切除時(shí)間卻較長(zhǎng)。這表明,機(jī)組A、B之間電氣距離的增加,比機(jī)組B、C之間電氣距離增加對(duì)機(jī)組C的暫態(tài)穩(wěn)定域的減小作用更加顯著。

    根據(jù)上述仿真結(jié)果可以推論:

    1) 機(jī)組之間的阻抗對(duì)輸出電流不流經(jīng)該等效阻抗的其他機(jī)組的暫態(tài)穩(wěn)定性影響很小,但對(duì)輸出電流流經(jīng)該等效阻抗的其他機(jī)組的暫態(tài)穩(wěn)定性影響較大。

    2) 當(dāng)風(fēng)場(chǎng)中與機(jī)組之間的等效阻抗上有多臺(tái)機(jī)組的輸出電流流過(guò)時(shí),則該等效阻抗對(duì)場(chǎng)內(nèi)各機(jī)組的暫態(tài)穩(wěn)定性影響會(huì)很大,隨著場(chǎng)內(nèi)機(jī)組數(shù)的增多,此時(shí)將不能簡(jiǎn)單忽略場(chǎng)內(nèi)阻抗的影響。

    4 結(jié)論

    本文主要研究了在考慮機(jī)組有功出力不同且各機(jī)組之間相互作用的情況下,風(fēng)電場(chǎng)及其內(nèi)部各機(jī)組的暫態(tài)穩(wěn)定性及其穩(wěn)定裕度,同時(shí)對(duì)場(chǎng)內(nèi)機(jī)組之間的阻抗對(duì)機(jī)組暫態(tài)穩(wěn)定性的影響作了定性分析與仿真驗(yàn)證。主要結(jié)論有:

    1) 風(fēng)電場(chǎng)暫態(tài)失穩(wěn)是由PLL的非線性動(dòng)力學(xué)行為引起的,PLL平衡點(diǎn)發(fā)生改變及失穩(wěn)都將引起風(fēng)電場(chǎng)的暫態(tài)穩(wěn)定性問(wèn)題。

    2) 在不考慮場(chǎng)內(nèi)阻抗的前提下,對(duì)于同一低壓饋線上的機(jī)組來(lái)說(shuō),電網(wǎng)故障后若有功出力最大的機(jī)組暫態(tài)失穩(wěn),則其他機(jī)組早已暫態(tài)失穩(wěn);若有功出力最小的機(jī)組暫態(tài)穩(wěn)定,則其他機(jī)組也暫態(tài)穩(wěn)定。

    3) 機(jī)組之間的阻抗對(duì)輸出電流不流經(jīng)該等效阻抗的其他機(jī)組的暫態(tài)穩(wěn)定性影響很小,但對(duì)輸出電流流經(jīng)該等效阻抗的其他機(jī)組的暫態(tài)穩(wěn)定性影響較大。對(duì)于不同低壓饋線上的機(jī)組而言,某一低壓饋線上機(jī)組之間的阻抗對(duì)另一低壓饋線上機(jī)組的暫態(tài)穩(wěn)定性影響很小。

    4) 若機(jī)組之間的等效阻抗上有多臺(tái)機(jī)組的輸出電流流過(guò),則該等效阻抗對(duì)這些機(jī)組的暫態(tài)穩(wěn)定性影響會(huì)很大。在場(chǎng)內(nèi)同一饋線上機(jī)組數(shù)量較多的情況下,此類阻抗對(duì)于機(jī)組的暫態(tài)穩(wěn)定性影響將不能被簡(jiǎn)單忽略。

    本文得到的暫態(tài)穩(wěn)定性相關(guān)結(jié)論能夠?qū)﹄娋W(wǎng)故障期間場(chǎng)內(nèi)單臺(tái)機(jī)組的暫態(tài)穩(wěn)定性的判斷以及場(chǎng)內(nèi)阻抗對(duì)機(jī)組的暫態(tài)穩(wěn)定性影響提供一定的參考依據(jù)。需要指出的是,本文針對(duì)全功率機(jī)組模型做了簡(jiǎn)化處理,并未考慮機(jī)組的非理想電流源特性,即電壓外環(huán)在電網(wǎng)故障期間的影響。同時(shí),基于各機(jī)組有功電流的近似處理存在著一些誤差,本文未考慮系統(tǒng)自身阻尼對(duì)暫態(tài)穩(wěn)定性的影響,導(dǎo)致所得暫態(tài)穩(wěn)定裕度較為保守。此外,本文定性分析了場(chǎng)內(nèi)機(jī)組之間的阻抗對(duì)各機(jī)組的暫態(tài)穩(wěn)定性影響并進(jìn)行了仿真驗(yàn)證,但未進(jìn)行深入討論。上述問(wèn)題是本文的不足之處,也是將來(lái)此類暫態(tài)穩(wěn)定性問(wèn)題的研究重點(diǎn)。

    [1] 張雪娟, 束洪春, 孫士云, 等. 雙饋風(fēng)機(jī)參與系統(tǒng)調(diào)頻對(duì)系統(tǒng)暫態(tài)功角穩(wěn)定性的影響分析[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(2): 18-29.

    ZHANG Xuejuan, SHU Hongchun, SUN Shiyun, at al. Effect analysis of a doubly-fed induction generator participating in system frequency modulation on system transient power angle stability[J]. Power System Protection and Control, 2021, 49(2): 18-29.

    [2] 盧一菲, 陳沖, 金成日, 等. 直驅(qū)永磁風(fēng)電機(jī)組高電壓穿越協(xié)調(diào)控制策略[J]. 電力系統(tǒng)保護(hù)與控制, 2020, 48(15): 50-60.

    LU Yifei, CHEN Chong, JIN Chengri, et al. HVRT coordinated control strategy of a direct-driven PMSG[J]. Power System Protection and Control, 2020, 48(15): 50-60.

    [3] NADOUR M, ESSADKI A, NASSER T. Improving low-voltage ride-through capability of a multimegawatt DFIG based wind turbine under grid faults[J]. Protection and Control of Modern Power Systems, 2020, 5(4): 370-382.

    [4] 李東東, 孫雅茹, 徐波, 等. 考慮頻率穩(wěn)定的新能源高滲透率電力系統(tǒng)最小慣量與一次調(diào)頻容量評(píng)估方法[J].電力系統(tǒng)保護(hù)與控制, 2021, 49(23): 54-61.

    LI Dongdong, SUN Yaru, XU Bo, et al. Minimum inertia and primary frequency capacity assessment for a new energy highpermeability power system considering frequency stability[J]. Power System Protection and Control, 2021, 49(23): 54-61.

    [5] 晁凱云, 苗世洪, 劉子文, 等. 基于虛擬同步電機(jī)控制的微電網(wǎng)穩(wěn)定性分析[J]. 電力系統(tǒng)保護(hù)與控制, 2019, 47(3): 9-16.

    CHAO Kaiyun, MIAO Shihong, LIU Ziwen, et al. Stability analysis of microgrid based on virtual synchronous generator control[J]. Power System Protection and Control, 2019, 47(3): 9-16.

    [6] 孫亮, 楊曉飛, 孫立國(guó), 等. 基于改進(jìn)虛擬同步發(fā)電機(jī)的多逆變器頻率無(wú)差控制策略[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(11): 18-27.

    SUN Liang, YANG Xiaofei, SUN Liguo, et al. Frequent deviation-free control for microgrid multi-inverters based on improving a virtual synchronous generator[J]. Power System Protection and Control, 2021, 49(11): 18-27.

    [7] 張成, 趙濤, 朱愛(ài)華, 等. 提高弱電網(wǎng)下并網(wǎng)逆變器穩(wěn)定性的復(fù)合補(bǔ)償策略[J]. 智慧電力, 2021, 49(5): 42-47.

    ZHANG Cheng, ZHAO Tao, ZHU Aihua, et al. Composite compensation strategy for improving stability of grid connected inverter in weak current network[J]. Smart Power, 2021, 49(5): 42-47.

    [8] 歐陽(yáng)逸風(fēng), 鄒宇. 弱電網(wǎng)條件下并網(wǎng)逆變器的鎖相環(huán)靜態(tài)穩(wěn)定分析[J]. 電力系統(tǒng)保護(hù)與控制, 2018, 46(18): 74-79.

    OUYANG Yifeng, ZOU Yu. Static stability analysis of phase-locked loop in grid-connected inverters under weak grid condition[J]. Power System Protection and Control, 2018, 46(18): 74-79.

    [9] 孫華東, 張振宇, 林偉芳, 等. 2011年西北電網(wǎng)風(fēng)機(jī)脫網(wǎng)事故分析及啟示[J]. 電網(wǎng)技術(shù), 2012, 36(10): 76-80.

    SUN Huadong, ZHANG Zhenyu, LIN Weifang, et al. Analysis on serious wind turbine generators tripping accident in northwest China power grid in 2011 and its lessons[J]. Power System Technology, 2012, 36(10): 76-80.

    [10] ZHOU J Z, GOLE A M. VSC transmission limitations imposed by AC system strength and AC impedance characteristics[C] // 10th IET International Conference on AC and DC Power Transmission, December 4-5, 2012, Birmingham, Britain: 1-6.

    [11] ZHOU J Z, DING H, FAN S, et al. Impact of short-circuit ratio and phase-locked-loop parameters on the small-signal behavior of a VSC-HVDC converter[J]. IEEE Transactions on Power Delivery, 2014, 29(5): 2287-2296.

    [12] HUANG Y, YUAN X, HU J, et al. Modeling of VSC connected to weak grid for stability analysis of DC-link voltage control[J]. IEEE Journal of Emerging & Selected Topics in Power Electronics, 2017, 3(4): 1193-1204.

    [13] GOKSU O, TEODORESCU R, BAK C L, et al. Instability of wind turbine converters during current injection to low voltage grid faults and PLL frequency based stability solution[J]. IEEE Transactions on Power Systems, 2014, 29(4): 1683-1691.

    [14] HU J, WANG B, WANG W, et al. Small signal dynamics of DFIG-based wind turbines during riding through symmetrical faults in weak AC grid[J]. IEEE Transactions on Energy Conversion, 2017, 32(2): 720-730.

    [15] HU J, HU Q, WANG B, et al. Small signal instability of PLL-synchronized type-4 wind turbines connected to high-impedance AC grid during LVRT[J]. IEEE Transactions on Energy Conversion, 2016, 31(4): 1676-1687.

    [16]張學(xué)廣, 馬彥, 王天一, 等. 弱電網(wǎng)下雙饋發(fā)電機(jī)輸入導(dǎo)納建模及穩(wěn)定性分析[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2017, 37(5): 1507-1515.

    ZHANG Xueguang, MA Yan, WANG Tianyi, et al. Input admittance modeling and stability analysis of DFIG under weak grid condition[J]. Proceedings of the CSEE, 2017, 37(5): 1507-1515.

    [17] ZAMANI M A, YAZDANI A, SIDHU T S. A control strategy for enhanced operation of inverter-based microgrids under transient disturbances and network faults[J]. IEEE Transactions on Power Delivery, 2012, 27(4): 1737-1747.

    [18] GUO X, ZHANG X, WANG B, et al. Asymmetrical grid fault ride-through strategy of three-phase grid-connected inverter considering network impedance impact in low-voltage grid[J]. IEEE Transactions on Power Electronics, 2013, 29(3): 1064-1068.

    [19] HUANG M, PENG Y, TSE C K, et al. Bifurcation and large-signal stability analysis of three-phase voltage source converter under grid voltage dips[J]. IEEE Transactions on Power Electronics, 2017, 32(11): 8868-8879.

    [20] MIDTSUND T, SUUL J A, UNDELAND T. Evaluation of current controller performance and stability for voltage source converters connected to a weak grid[C] // The 2nd International Symposium on Power Electronics for Distributed Generation Systems (PEDG), August 12-15, 2010, Hefei, China: 382-388.

    [21] ERLICH I, SHEWAREGA F, ENGELHARDT S, et al. Effect of wind turbine output current during faults on grid voltage and the transient stability of wind parks[C] // IEEE Power & Energy Society General Meeting, October 2-5, 2009, Calgary, AB, Canada: 1-8.

    [22] MA S K, GENG H, LIU L, et al. Grid-synchronization stability improvement of large scale wind farm during severe grid fault[J]. IEEE Transactions on Power Systems, 2018, 33(1): 216-226.

    [23] WU H, WANG X. Design-oriented transient stability analysis of PLL-synchronized voltage-source converters[J]. IEEE Transactions on Power Electronics, 2020, 35(4): 3573-3599.

    [24] HE X, GENG H, XI J, et al. Resynchronization analysis and improvement of grid-connected VSCs during grid faults[J]. IEEE Journal of Emerging and Selected Topics in Power Electronics, 2021, 9(1): 438-450.

    [25] WU H, WANG X. Transient stability impact of the phase-locked loop on grid-connected voltage source converters[C] // International Power Electronics Conference (IPEC-ECCE Asia), October 25-28, 2018, Niigata, Japan: 2673-2680.

    [26] 張琛, 蔡旭, 李征. 全功率變換風(fēng)電機(jī)組的暫態(tài)穩(wěn)定性分析[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2017, 37(14): 4018-4026.

    ZHANG Chen, CAI Xu, LI Zheng. Transient stability analysis of wind turbines with full-scale voltage source converter[J]. Proceedings of the CSEE, 2017, 37(14): 4018-4026.

    [27] 韓剛, 張琛, 蔡旭. 電網(wǎng)短路故障引發(fā)的全功率風(fēng)電機(jī)組頻率失穩(wěn)機(jī)理與控制方法[J]. 電工技術(shù)學(xué)報(bào), 2018, 33(10): 2167-2175.

    HAN Gang, ZHANG Chen, CAI Xu. Mechanism of frequency instability of full-scale wind turbines caused by grid short circuit fault and its control method[J]. Transactions of China Electrotechnical Society, 2018, 33(10): 2167-2175.

    [28] ZHANG C, CAI X, RYGG A, et al. Modeling and analysis of grid-synchronizing stability of a type-IV wind turbine under grid faults[J]. International Journal of Electrical Power & Energy Systems, 2019, 117.

    [29] 張琛, 李征, 蔡旭, 等. 面向電力系統(tǒng)暫態(tài)穩(wěn)定分析的雙饋風(fēng)電機(jī)組動(dòng)態(tài)模型[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2016, 36(20): 5449-5460.

    ZHANG Chen, LI Zheng, CAI Xu, et al. Dynamic model of DFIG wind turbines for power system transient stability analysis[J]. Proceedings of the CSEE, 2016, 36(20): 5449-5460.

    Transient stability analysis of wind turbines with a full-scale converter under grid fault

    LIU Feiyang, ZENG Ping, LI Zheng

    (College of Information Science and Technology, Donghua University, Shanghai 201620, China)

    The stability of grid-connected operation of full-scale wind turbines in wind farms is often related to the control of the grid-side converter of the turbine and the impedance between grid and the turbine. Therefore, the instability of wind turbines in a wind farm does not always occur at the same time. In order to reveal the relationship between the location of wind turbines and their transient stability, as well as the interaction on stability between wind turbines, a nonlinear reduced order model of multi wind turbines based on phase-locked loop is established. Based on the theory of equal-area-principle, the transient stability margin of different turbines in the wind farm under different operating conditions and the influence of line impedance on the stability are then analyzed. Through Matlab/Simulink time domain simulation, quantitative analysis and verification are carried out. The results show that the feeder impedance in the wind farm will reduce the transient stability region of the connected turbines whose output current flow through the impedance, and has little effect on the transient stability of other turbines whose output current does not flow through the impedance. The greater the current flowing through the impedance is, the smaller the stability margin. Therefore, with the increase of the number of turbines linked on a feeder line, the impact of the internal impedance of a wind farm on stability cannot be simply ignored. In addition, the larger output of a wind turbine, the easier it is for it to become unstable under the same conditions.

    grid fault; full-scale wind farm; transient stability; phase-locked loop; equal-area-principle; stability margin;internal impedance of wind farm

    10.19783/j.cnki.pspc.210362

    國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目資助(51837007)

    This work is supported by the Key Project of the Natural Science Foundation of China (No. 51837007).

    2021-04-06;

    2021-12-18

    柳飛揚(yáng)(1990—),男,碩士,研究方向?yàn)轱L(fēng)電機(jī)組的建模及穩(wěn)定性分析;E-mail: 2191365@mail.dhu.edu.cn

    曾 平(1998—),男,碩士,研究方向?yàn)轱L(fēng)電并網(wǎng)的暫態(tài)穩(wěn)定性分析;E-mail: 150910218@mail.dhu.edu.cn

    李 征(1961—),女,通信作者,博士,教授,研究方向?yàn)樾履茉聪到y(tǒng)建模與并網(wǎng)控制、儲(chǔ)能系統(tǒng)應(yīng)用、智能微電網(wǎng)。E-mail: lizheng@dhu.edu.cn

    (編輯 周金梅)

    猜你喜歡
    變流器暫態(tài)風(fēng)電場(chǎng)
    300Mvar空冷隱極同步調(diào)相機(jī)暫態(tài)特性仿真分析
    基于PSS/E的風(fēng)電場(chǎng)建模與動(dòng)態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    電力系統(tǒng)全網(wǎng)一體化暫態(tài)仿真接口技術(shù)
    電子制作(2018年14期)2018-08-21 01:38:28
    除氧器暫態(tài)計(jì)算研究
    中壓一體化儲(chǔ)能變流器的設(shè)計(jì)
    含風(fēng)電場(chǎng)電力系統(tǒng)的潮流計(jì)算
    基于背靠背變流器的并網(wǎng)控制研究
    改進(jìn)PR控制在直驅(qū)風(fēng)機(jī)變流器中的應(yīng)用
    探求風(fēng)電場(chǎng)的遠(yuǎn)景
    風(fēng)能(2015年9期)2015-02-27 10:15:25
    代力吉風(fēng)電場(chǎng)的我們
    風(fēng)能(2015年7期)2015-02-27 10:15:02
    亚洲午夜精品一区,二区,三区| 久久久久久久久免费视频了| 交换朋友夫妻互换小说| 琪琪午夜伦伦电影理论片6080| 国产av精品麻豆| 成年人黄色毛片网站| 黑人猛操日本美女一级片| 国产亚洲精品综合一区在线观看 | 国产精华一区二区三区| 国产高清国产精品国产三级| 亚洲 欧美 日韩 在线 免费| 99热只有精品国产| 麻豆av在线久日| 国产亚洲欧美精品永久| 欧美激情久久久久久爽电影 | 91av网站免费观看| 久久精品国产99精品国产亚洲性色 | 国产精品美女特级片免费视频播放器 | 精品国产国语对白av| 亚洲国产欧美日韩在线播放| 黄色视频不卡| 男人操女人黄网站| 日日摸夜夜添夜夜添小说| 久久久久久免费高清国产稀缺| 国产精品久久久人人做人人爽| 欧美黄色片欧美黄色片| 日韩人妻精品一区2区三区| 欧美成狂野欧美在线观看| 97超级碰碰碰精品色视频在线观看| 日日摸夜夜添夜夜添小说| 免费在线观看影片大全网站| 日韩欧美一区视频在线观看| 在线观看午夜福利视频| 人人妻人人添人人爽欧美一区卜| 欧美一区二区精品小视频在线| 十八禁网站免费在线| 欧美成狂野欧美在线观看| 久久久久久久久免费视频了| 999久久久国产精品视频| 一二三四社区在线视频社区8| 国产黄色免费在线视频| 50天的宝宝边吃奶边哭怎么回事| 国产无遮挡羞羞视频在线观看| 日韩欧美三级三区| 欧美日韩福利视频一区二区| 在线国产一区二区在线| 91字幕亚洲| 国产一区二区在线av高清观看| 一区在线观看完整版| 夜夜躁狠狠躁天天躁| 一区二区三区激情视频| 亚洲精品国产精品久久久不卡| 69精品国产乱码久久久| 丝袜美足系列| 国产在线观看jvid| 人成视频在线观看免费观看| 成熟少妇高潮喷水视频| 在线观看免费视频日本深夜| 亚洲精品美女久久久久99蜜臀| 免费看十八禁软件| 最新在线观看一区二区三区| 少妇的丰满在线观看| 久久久国产一区二区| 欧美最黄视频在线播放免费 | 18美女黄网站色大片免费观看| 午夜久久久在线观看| 777久久人妻少妇嫩草av网站| 亚洲精品久久成人aⅴ小说| 人成视频在线观看免费观看| 黄片播放在线免费| 好男人电影高清在线观看| 中文字幕人妻丝袜一区二区| 欧美日韩视频精品一区| 看免费av毛片| 女人爽到高潮嗷嗷叫在线视频| 精品国产超薄肉色丝袜足j| 日韩欧美三级三区| 真人一进一出gif抽搐免费| 淫妇啪啪啪对白视频| 久久婷婷成人综合色麻豆| 十八禁网站免费在线| 99香蕉大伊视频| 免费在线观看日本一区| 波多野结衣高清无吗| 俄罗斯特黄特色一大片| 身体一侧抽搐| 国产成人精品在线电影| 久久久国产一区二区| 九色亚洲精品在线播放| 免费高清视频大片| 久久久精品国产亚洲av高清涩受| 国产精品久久久久成人av| 成人三级黄色视频| 69av精品久久久久久| 国产成人啪精品午夜网站| 久久国产精品影院| 欧美激情极品国产一区二区三区| 久久精品国产99精品国产亚洲性色 | 国产1区2区3区精品| 人人妻,人人澡人人爽秒播| 丝袜美腿诱惑在线| 国产激情久久老熟女| 水蜜桃什么品种好| 国产黄a三级三级三级人| 国产极品粉嫩免费观看在线| 国产精品国产高清国产av| 80岁老熟妇乱子伦牲交| 99香蕉大伊视频| 亚洲av片天天在线观看| 亚洲国产精品sss在线观看 | 777久久人妻少妇嫩草av网站| 午夜免费观看网址| 日韩一卡2卡3卡4卡2021年| 日韩免费高清中文字幕av| 国产又爽黄色视频| 男人操女人黄网站| 热99国产精品久久久久久7| 欧美黑人精品巨大| 久9热在线精品视频| 男男h啪啪无遮挡| 欧美日本亚洲视频在线播放| 十八禁人妻一区二区| 一级毛片女人18水好多| 在线观看一区二区三区激情| 久久久久九九精品影院| 亚洲午夜理论影院| 一进一出好大好爽视频| 亚洲欧美激情综合另类| 不卡av一区二区三区| 日本 av在线| www.www免费av| 亚洲 欧美一区二区三区| 国产av在哪里看| 99国产精品一区二区蜜桃av| 精品一区二区三区av网在线观看| 色在线成人网| 精品福利永久在线观看| 日日爽夜夜爽网站| 国产av一区二区精品久久| 欧美精品一区二区免费开放| 国产精品久久久久成人av| 久热这里只有精品99| 在线观看免费日韩欧美大片| 在线观看舔阴道视频| 午夜福利欧美成人| 国产欧美日韩综合在线一区二区| 美女国产高潮福利片在线看| 中文字幕人妻丝袜一区二区| 欧美成人免费av一区二区三区| 欧美成人午夜精品| 久久香蕉国产精品| 高清在线国产一区| 新久久久久国产一级毛片| 大型av网站在线播放| 久久人妻av系列| 国产深夜福利视频在线观看| 十八禁人妻一区二区| 国产精品野战在线观看 | 9色porny在线观看| 日韩成人在线观看一区二区三区| 电影成人av| 人人澡人人妻人| www.自偷自拍.com| av网站免费在线观看视频| 久久热在线av| av天堂在线播放| 老鸭窝网址在线观看| 午夜免费鲁丝| av在线播放免费不卡| 青草久久国产| 欧美乱码精品一区二区三区| 少妇被粗大的猛进出69影院| 18禁黄网站禁片午夜丰满| 99国产精品一区二区三区| 久久中文看片网| 美女 人体艺术 gogo| 久久精品亚洲精品国产色婷小说| 岛国视频午夜一区免费看| 日本撒尿小便嘘嘘汇集6| 亚洲av五月六月丁香网| 国产一区二区三区在线臀色熟女 | 亚洲精品粉嫩美女一区| 久久久国产成人精品二区 | 国产精品秋霞免费鲁丝片| www日本在线高清视频| 亚洲第一青青草原| 精品久久久久久电影网| 999久久久精品免费观看国产| 淫秽高清视频在线观看| 高清在线国产一区| 久久欧美精品欧美久久欧美| 麻豆av在线久日| 国产成人免费无遮挡视频| 国产又色又爽无遮挡免费看| 很黄的视频免费| 国产精品一区二区精品视频观看| 大陆偷拍与自拍| 免费av中文字幕在线| 久久人人97超碰香蕉20202| 亚洲欧美一区二区三区黑人| 热re99久久国产66热| 日韩精品免费视频一区二区三区| 精品一区二区三区四区五区乱码| 一区二区三区激情视频| 免费搜索国产男女视频| 麻豆一二三区av精品| 国产主播在线观看一区二区| 国产熟女午夜一区二区三区| 午夜激情av网站| 国产三级黄色录像| 欧美大码av| 乱人伦中国视频| 国产精品久久视频播放| 99国产精品一区二区三区| 一二三四社区在线视频社区8| 亚洲第一欧美日韩一区二区三区| 老司机靠b影院| 天堂中文最新版在线下载| 精品久久久久久久久久免费视频 | 人人妻人人澡人人看| 欧美不卡视频在线免费观看 | 女人爽到高潮嗷嗷叫在线视频| 在线永久观看黄色视频| 女人精品久久久久毛片| 久久婷婷成人综合色麻豆| 一夜夜www| 亚洲欧美日韩高清在线视频| 一级a爱片免费观看的视频| 69精品国产乱码久久久| 亚洲av五月六月丁香网| 丝袜美腿诱惑在线| 最新在线观看一区二区三区| 妹子高潮喷水视频| a级毛片在线看网站| 色播在线永久视频| 免费高清视频大片| 精品午夜福利视频在线观看一区| 嫩草影视91久久| 这个男人来自地球电影免费观看| 亚洲午夜理论影院| 亚洲精品国产精品久久久不卡| 多毛熟女@视频| 嫩草影视91久久| 老鸭窝网址在线观看| 精品卡一卡二卡四卡免费| 国产一卡二卡三卡精品| 国产精品久久电影中文字幕| 美女大奶头视频| 99精品久久久久人妻精品| 午夜激情av网站| av在线天堂中文字幕 | 久久精品国产亚洲av香蕉五月| 视频区图区小说| 一级黄色大片毛片| 免费搜索国产男女视频| 国产免费男女视频| 亚洲av电影在线进入| 制服人妻中文乱码| 久久精品91蜜桃| 日本精品一区二区三区蜜桃| 久99久视频精品免费| 欧美 亚洲 国产 日韩一| 91麻豆精品激情在线观看国产 | 宅男免费午夜| 国产成人啪精品午夜网站| 亚洲精品一区av在线观看| 在线观看免费日韩欧美大片| 国产精品野战在线观看 | 狠狠狠狠99中文字幕| 国产亚洲欧美在线一区二区| www.www免费av| 国产精品电影一区二区三区| 日韩欧美国产一区二区入口| 丝袜美腿诱惑在线| 成人三级做爰电影| 精品国产超薄肉色丝袜足j| 男女之事视频高清在线观看| 亚洲情色 制服丝袜| 久久久久久人人人人人| 久久久久久久精品吃奶| 亚洲在线自拍视频| 精品久久久久久,| 在线观看免费日韩欧美大片| 亚洲欧美一区二区三区久久| 一区二区日韩欧美中文字幕| 99久久国产精品久久久| 亚洲av五月六月丁香网| 黄片小视频在线播放| 日本撒尿小便嘘嘘汇集6| 又黄又爽又免费观看的视频| 亚洲成人国产一区在线观看| 精品国产一区二区久久| 一本大道久久a久久精品| 成人三级黄色视频| 69av精品久久久久久| 美女扒开内裤让男人捅视频| av免费在线观看网站| 欧美日韩乱码在线| 丁香六月欧美| 国产在线观看jvid| 日本免费a在线| 亚洲久久久国产精品| 亚洲欧美激情综合另类| 久久精品亚洲熟妇少妇任你| 国产精品综合久久久久久久免费 | 在线观看一区二区三区| 久久久精品国产亚洲av高清涩受| 国产视频一区二区在线看| 久久久国产精品麻豆| 人妻久久中文字幕网| 色老头精品视频在线观看| 成人亚洲精品一区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 少妇 在线观看| 男女床上黄色一级片免费看| 亚洲人成伊人成综合网2020| 美女国产高潮福利片在线看| 男女床上黄色一级片免费看| x7x7x7水蜜桃| 亚洲一区中文字幕在线| 岛国视频午夜一区免费看| 国产精品秋霞免费鲁丝片| 亚洲国产精品合色在线| 亚洲,欧美精品.| 中文字幕人妻熟女乱码| 亚洲欧美一区二区三区久久| 亚洲精品久久成人aⅴ小说| 欧美成狂野欧美在线观看| 99久久99久久久精品蜜桃| 国产又色又爽无遮挡免费看| 久久久久久久久久久久大奶| 午夜免费鲁丝| 免费少妇av软件| 亚洲国产精品一区二区三区在线| av网站在线播放免费| 亚洲一区中文字幕在线| 精品国产一区二区三区四区第35| a级毛片在线看网站| 在线观看免费高清a一片| 在线观看免费视频网站a站| 91国产中文字幕| 免费一级毛片在线播放高清视频 | avwww免费| 日韩大码丰满熟妇| 一级毛片女人18水好多| 国产男靠女视频免费网站| 18禁美女被吸乳视频| 伦理电影免费视频| ponron亚洲| 最近最新中文字幕大全免费视频| 中文字幕最新亚洲高清| 女人被躁到高潮嗷嗷叫费观| 天天添夜夜摸| 中文字幕最新亚洲高清| 久久久久亚洲av毛片大全| 大型av网站在线播放| 制服人妻中文乱码| 男人舔女人的私密视频| 天堂中文最新版在线下载| 正在播放国产对白刺激| 少妇粗大呻吟视频| а√天堂www在线а√下载| 成熟少妇高潮喷水视频| 久久久久久免费高清国产稀缺| 免费高清在线观看日韩| www国产在线视频色| 久久精品91蜜桃| 精品久久久久久久毛片微露脸| 国产单亲对白刺激| 亚洲欧美日韩另类电影网站| 长腿黑丝高跟| 国产97色在线日韩免费| 在线观看66精品国产| 日本黄色日本黄色录像| 免费在线观看黄色视频的| 亚洲国产欧美网| 一a级毛片在线观看| 制服诱惑二区| 久久青草综合色| 久久国产精品人妻蜜桃| 精品电影一区二区在线| 天堂中文最新版在线下载| 级片在线观看| 大型黄色视频在线免费观看| 男女之事视频高清在线观看| svipshipincom国产片| 老司机靠b影院| 亚洲精品中文字幕一二三四区| 国产激情久久老熟女| 久久久久久免费高清国产稀缺| 天天添夜夜摸| 中文字幕最新亚洲高清| www.自偷自拍.com| 欧美日韩亚洲综合一区二区三区_| 夫妻午夜视频| 嫩草影视91久久| 一夜夜www| 50天的宝宝边吃奶边哭怎么回事| 亚洲全国av大片| 免费搜索国产男女视频| 黄色a级毛片大全视频| 午夜免费成人在线视频| 在线观看午夜福利视频| 91老司机精品| 久久欧美精品欧美久久欧美| 夜夜爽天天搞| 免费在线观看日本一区| 一区二区日韩欧美中文字幕| 黄色片一级片一级黄色片| 手机成人av网站| 长腿黑丝高跟| 午夜成年电影在线免费观看| 999久久久国产精品视频| 中文字幕人妻丝袜制服| 丰满迷人的少妇在线观看| 99国产极品粉嫩在线观看| 亚洲精华国产精华精| 亚洲中文日韩欧美视频| 少妇被粗大的猛进出69影院| 久久久精品国产亚洲av高清涩受| 国产精品99久久99久久久不卡| 美女高潮到喷水免费观看| 老司机午夜福利在线观看视频| 丰满饥渴人妻一区二区三| 大型av网站在线播放| 一级毛片精品| 国产精品久久久人人做人人爽| 国产精品秋霞免费鲁丝片| 国产精品98久久久久久宅男小说| av有码第一页| 成年女人毛片免费观看观看9| 精品一区二区三区视频在线观看免费 | 久久精品国产清高在天天线| 国产极品粉嫩免费观看在线| 欧美 亚洲 国产 日韩一| 欧美日本中文国产一区发布| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区视频在线观看免费 | 亚洲国产欧美一区二区综合| 无限看片的www在线观看| 亚洲国产看品久久| 中文字幕人妻熟女乱码| 精品卡一卡二卡四卡免费| 日韩欧美一区视频在线观看| 国产精品秋霞免费鲁丝片| 又紧又爽又黄一区二区| 大陆偷拍与自拍| 最近最新中文字幕大全电影3 | 亚洲一码二码三码区别大吗| 国产99久久九九免费精品| 男女下面进入的视频免费午夜 | 成人永久免费在线观看视频| 午夜免费观看网址| 欧美一级毛片孕妇| 午夜精品国产一区二区电影| 夜夜爽天天搞| 在线观看舔阴道视频| 男女床上黄色一级片免费看| 欧美丝袜亚洲另类 | 久久精品国产清高在天天线| 欧美丝袜亚洲另类 | 国产欧美日韩精品亚洲av| 无人区码免费观看不卡| 亚洲熟女毛片儿| 亚洲av电影在线进入| 制服人妻中文乱码| 99re在线观看精品视频| 男人操女人黄网站| 精品福利观看| 免费观看人在逋| 亚洲av片天天在线观看| 99精品欧美一区二区三区四区| 精品日产1卡2卡| 纯流量卡能插随身wifi吗| 免费人成视频x8x8入口观看| 一级毛片高清免费大全| 在线免费观看的www视频| www国产在线视频色| 国产精品综合久久久久久久免费 | 一区二区日韩欧美中文字幕| 99久久国产精品久久久| 欧美亚洲日本最大视频资源| 久久99一区二区三区| 看免费av毛片| 精品人妻在线不人妻| 亚洲欧美精品综合久久99| 欧美+亚洲+日韩+国产| 如日韩欧美国产精品一区二区三区| 热99国产精品久久久久久7| 欧美激情极品国产一区二区三区| 男女高潮啪啪啪动态图| 精品一区二区三卡| 国产精品秋霞免费鲁丝片| 国产99白浆流出| av超薄肉色丝袜交足视频| 高清欧美精品videossex| 国产亚洲欧美在线一区二区| 黑人猛操日本美女一级片| 天天影视国产精品| 色婷婷av一区二区三区视频| 老熟妇仑乱视频hdxx| 国产黄色免费在线视频| 黄色成人免费大全| 99国产综合亚洲精品| 亚洲 欧美一区二区三区| 国产av一区二区精品久久| 欧美一级毛片孕妇| 午夜影院日韩av| 国产精品秋霞免费鲁丝片| 欧美+亚洲+日韩+国产| 亚洲人成电影观看| 黄色成人免费大全| 欧美成人性av电影在线观看| 色综合婷婷激情| 久久久久久久午夜电影 | 一级a爱片免费观看的视频| 日本撒尿小便嘘嘘汇集6| ponron亚洲| 超色免费av| 9热在线视频观看99| 亚洲九九香蕉| 国产精品野战在线观看 | xxx96com| 人人妻人人爽人人添夜夜欢视频| 黑人猛操日本美女一级片| 国产一区在线观看成人免费| 成年女人毛片免费观看观看9| 黑人巨大精品欧美一区二区蜜桃| 桃红色精品国产亚洲av| 侵犯人妻中文字幕一二三四区| 99国产精品一区二区蜜桃av| 欧美乱色亚洲激情| 国产一区二区激情短视频| 久久国产亚洲av麻豆专区| 久久久水蜜桃国产精品网| 久久人人爽av亚洲精品天堂| 国产一卡二卡三卡精品| 国内久久婷婷六月综合欲色啪| 日本欧美视频一区| 欧美中文日本在线观看视频| 91成人精品电影| 国产精品亚洲av一区麻豆| 国产精品影院久久| 国产国语露脸激情在线看| 亚洲国产欧美日韩在线播放| 国产av又大| 妹子高潮喷水视频| 超色免费av| 一级a爱视频在线免费观看| 免费人成视频x8x8入口观看| 欧美日韩瑟瑟在线播放| 亚洲欧美一区二区三区黑人| 免费在线观看视频国产中文字幕亚洲| 俄罗斯特黄特色一大片| 好看av亚洲va欧美ⅴa在| 亚洲av电影在线进入| 色播在线永久视频| 久久久久久大精品| 日韩大尺度精品在线看网址 | 少妇被粗大的猛进出69影院| 国产97色在线日韩免费| 一进一出抽搐gif免费好疼 | 亚洲av成人一区二区三| 国产蜜桃级精品一区二区三区| 欧美日韩福利视频一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 精品无人区乱码1区二区| 亚洲国产精品sss在线观看 | 女人爽到高潮嗷嗷叫在线视频| 成人亚洲精品一区在线观看| 亚洲中文日韩欧美视频| 91成年电影在线观看| 色综合婷婷激情| 极品教师在线免费播放| 91在线观看av| 亚洲性夜色夜夜综合| 香蕉久久夜色| 国产有黄有色有爽视频| 国产免费男女视频| 日日干狠狠操夜夜爽| 亚洲国产精品合色在线| 亚洲aⅴ乱码一区二区在线播放 | 国产精品99久久99久久久不卡| 一个人观看的视频www高清免费观看 | 日韩免费高清中文字幕av| 黑丝袜美女国产一区| 久久精品国产综合久久久| 五月开心婷婷网| 国产成人精品久久二区二区免费| 免费高清在线观看日韩| 国产伦一二天堂av在线观看| 18禁美女被吸乳视频| 99国产精品一区二区三区| 久久久久久久精品吃奶| 国产精品98久久久久久宅男小说| 视频在线观看一区二区三区| 亚洲五月色婷婷综合| 亚洲五月婷婷丁香| 在线观看www视频免费| 久9热在线精品视频| 18禁观看日本| 国产麻豆69| 亚洲成人免费av在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 久久人人精品亚洲av| 国产区一区二久久| 亚洲国产精品sss在线观看 | 一级毛片高清免费大全| 国内毛片毛片毛片毛片毛片| 免费在线观看视频国产中文字幕亚洲| 国产亚洲精品久久久久5区| 欧美性长视频在线观看| 一级片免费观看大全| 国产亚洲精品久久久久久毛片| 国产精品偷伦视频观看了|