• <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
    亚洲精华国产精华精| 欧洲精品卡2卡3卡4卡5卡区| 欧美高清成人免费视频www| 欧美色欧美亚洲另类二区| 午夜免费激情av| 久久久久性生活片| 日本一二三区视频观看| 淫秽高清视频在线观看| 久久精品国产清高在天天线| 国产蜜桃级精品一区二区三区| av福利片在线观看| 日韩精品中文字幕看吧| 人妻制服诱惑在线中文字幕| 波多野结衣巨乳人妻| 最近最新中文字幕大全电影3| 久久精品国产99精品国产亚洲性色| 99在线视频只有这里精品首页| 99riav亚洲国产免费| 99久久无色码亚洲精品果冻| 搡老妇女老女人老熟妇| 国产一区二区三区在线臀色熟女| 丰满人妻一区二区三区视频av| 婷婷精品国产亚洲av在线| 观看免费一级毛片| 别揉我奶头 嗯啊视频| 天天躁日日操中文字幕| 超碰av人人做人人爽久久| 国产一区二区亚洲精品在线观看| 国产蜜桃级精品一区二区三区| 亚洲精品日韩av片在线观看| 精品久久久久久,| 婷婷精品国产亚洲av在线| 人人妻,人人澡人人爽秒播| 日韩免费av在线播放| 亚洲成人免费电影在线观看| 黄色丝袜av网址大全| 久久久久国内视频| 欧美成人a在线观看| 一区二区三区四区激情视频 | 欧美在线黄色| ponron亚洲| 日本黄色视频三级网站网址| 少妇高潮的动态图| 99国产精品一区二区三区| 久久久久国产精品人妻aⅴ院| 国产激情偷乱视频一区二区| 国产av麻豆久久久久久久| 国产av麻豆久久久久久久| 国产在线精品亚洲第一网站| 成人av在线播放网站| 久久九九热精品免费| 嫩草影院入口| 久久天躁狠狠躁夜夜2o2o| 一区二区三区免费毛片| 亚洲五月婷婷丁香| 午夜精品久久久久久毛片777| 国产精华一区二区三区| 免费观看的影片在线观看| 久久久色成人| 国产麻豆成人av免费视频| 亚洲人与动物交配视频| 99精品久久久久人妻精品| 国产精品美女特级片免费视频播放器| 久久九九热精品免费| 18禁黄网站禁片午夜丰满| 成人av在线播放网站| 国产日本99.免费观看| 特大巨黑吊av在线直播| 欧美日韩亚洲国产一区二区在线观看| 热99re8久久精品国产| 国内精品久久久久久久电影| 日本黄色片子视频| 欧洲精品卡2卡3卡4卡5卡区| 国产av一区在线观看免费| 欧美高清性xxxxhd video| 精品日产1卡2卡| 丰满人妻熟妇乱又伦精品不卡| www.熟女人妻精品国产| 最近最新中文字幕大全电影3| 成人一区二区视频在线观看| 欧美一级a爱片免费观看看| av专区在线播放| 9191精品国产免费久久| 两性午夜刺激爽爽歪歪视频在线观看| 婷婷亚洲欧美| 欧美不卡视频在线免费观看| 身体一侧抽搐| 日韩欧美 国产精品| 亚洲国产精品sss在线观看| 日韩成人在线观看一区二区三区| 午夜福利18| 99久久精品热视频| 久久精品综合一区二区三区| 美女xxoo啪啪120秒动态图 | 国产伦精品一区二区三区四那| 怎么达到女性高潮| 韩国av一区二区三区四区| 激情在线观看视频在线高清| 久久香蕉精品热| 国产亚洲精品久久久久久毛片| 听说在线观看完整版免费高清| 日韩欧美三级三区| 变态另类丝袜制服| 乱人视频在线观看| 亚洲av第一区精品v没综合| 欧美xxxx黑人xx丫x性爽| av女优亚洲男人天堂| 成年女人毛片免费观看观看9| 我要搜黄色片| 国产精品一区二区三区四区久久| 深夜精品福利| 亚洲性夜色夜夜综合| 久久亚洲真实| 最近在线观看免费完整版| 午夜激情福利司机影院| 免费黄网站久久成人精品 | 欧美又色又爽又黄视频| 美女免费视频网站| 少妇丰满av| 日韩欧美三级三区| 国内毛片毛片毛片毛片毛片| 18美女黄网站色大片免费观看| 久久久国产成人免费| 久久99热6这里只有精品| 日本三级黄在线观看| 国产欧美日韩精品一区二区| 欧美乱色亚洲激情| 少妇高潮的动态图| 老熟妇仑乱视频hdxx| 国内揄拍国产精品人妻在线| 国产一区二区亚洲精品在线观看| 制服丝袜大香蕉在线| 性插视频无遮挡在线免费观看| 国产精华一区二区三区| 两人在一起打扑克的视频| 中文字幕人妻熟人妻熟丝袜美| 美女被艹到高潮喷水动态| 成人高潮视频无遮挡免费网站| 亚洲 欧美 日韩 在线 免费| 亚洲自拍偷在线| 国产精品免费一区二区三区在线| 少妇的逼好多水| 黄色视频,在线免费观看| 国产精品久久视频播放| 久99久视频精品免费| 欧美成人免费av一区二区三区| 国产在视频线在精品| 亚洲三级黄色毛片| 午夜日韩欧美国产| 国产精品99久久久久久久久| 精品不卡国产一区二区三区| 一级a爱片免费观看的视频| 最新中文字幕久久久久| 亚洲av日韩精品久久久久久密| 日本a在线网址| 午夜精品久久久久久毛片777| 欧美3d第一页| 久久久精品大字幕| 在线十欧美十亚洲十日本专区| 日本成人三级电影网站| 身体一侧抽搐| 深夜a级毛片| 精品午夜福利在线看| 男女之事视频高清在线观看| 禁无遮挡网站| 国产成人欧美在线观看| 亚洲男人的天堂狠狠| 中文字幕久久专区| 特级一级黄色大片| 亚洲一区高清亚洲精品| 婷婷丁香在线五月| www.熟女人妻精品国产| 天堂网av新在线| 老女人水多毛片| 赤兔流量卡办理| 欧美激情国产日韩精品一区| 久久精品91蜜桃| 精品人妻视频免费看| 欧美精品国产亚洲| 一区二区三区四区激情视频 | 久久亚洲真实| 国产成人av教育| 国内久久婷婷六月综合欲色啪| 欧美潮喷喷水| 国内精品一区二区在线观看| 国内精品美女久久久久久| 久久人人精品亚洲av| 免费电影在线观看免费观看| 性色av乱码一区二区三区2| 禁无遮挡网站| 欧美成人免费av一区二区三区| 日韩大尺度精品在线看网址| 亚洲一区二区三区不卡视频| 黄色一级大片看看| 国内精品美女久久久久久| 在线观看一区二区三区| 国产精品三级大全| 一个人免费在线观看的高清视频| 男女之事视频高清在线观看| 国产中年淑女户外野战色| 尤物成人国产欧美一区二区三区| 欧美高清成人免费视频www| 99视频精品全部免费 在线| 身体一侧抽搐| 亚洲美女黄片视频| 精品久久国产蜜桃| 青草久久国产| 色综合欧美亚洲国产小说| 色5月婷婷丁香| 两个人的视频大全免费| 亚州av有码| 日韩高清综合在线| 免费看日本二区| 搡老熟女国产l中国老女人| av在线蜜桃| 狂野欧美白嫩少妇大欣赏| 国产高清三级在线| 网址你懂的国产日韩在线| 久久国产乱子伦精品免费另类| 国产91精品成人一区二区三区| 男女之事视频高清在线观看| 99久久精品国产亚洲精品| 成人精品一区二区免费| 精品国产亚洲在线| 男女做爰动态图高潮gif福利片| 一级毛片久久久久久久久女| 日韩有码中文字幕| 国产在线精品亚洲第一网站| 国产精品一区二区三区四区免费观看 | 99国产精品一区二区三区| 在线十欧美十亚洲十日本专区| 嫩草影视91久久| 国产精品三级大全| 国产黄色小视频在线观看| 桃色一区二区三区在线观看| 免费观看精品视频网站| 欧美最新免费一区二区三区 | 精品久久久久久久久av| 啪啪无遮挡十八禁网站| 成熟少妇高潮喷水视频| 麻豆成人午夜福利视频| 国产视频内射| 国产探花在线观看一区二区| 国产精品1区2区在线观看.| 日韩欧美精品v在线| 女生性感内裤真人,穿戴方法视频| 欧美成人一区二区免费高清观看| 亚洲第一欧美日韩一区二区三区| 亚洲第一电影网av| 亚洲久久久久久中文字幕| 国产精品永久免费网站| 别揉我奶头 嗯啊视频| 中文在线观看免费www的网站| 国产综合懂色| 成人av在线播放网站| 99riav亚洲国产免费| 一边摸一边抽搐一进一小说| 美女cb高潮喷水在线观看| 欧美乱妇无乱码| 成人性生交大片免费视频hd| 桃色一区二区三区在线观看| 亚洲av不卡在线观看| 久久伊人香网站| 国产大屁股一区二区在线视频| 又爽又黄a免费视频| 久久久精品大字幕| 可以在线观看毛片的网站| 18+在线观看网站| 在线国产一区二区在线| netflix在线观看网站| 熟妇人妻久久中文字幕3abv| 亚洲va日本ⅴa欧美va伊人久久| 欧美绝顶高潮抽搐喷水| 成人三级黄色视频| 中文字幕av成人在线电影| 最新中文字幕久久久久| 免费大片18禁| 非洲黑人性xxxx精品又粗又长| 18+在线观看网站| 久久人人精品亚洲av| 成人特级av手机在线观看| 精品人妻熟女av久视频| 啦啦啦韩国在线观看视频| 久久精品人妻少妇| 国产国拍精品亚洲av在线观看| 亚洲国产精品成人综合色| 久久草成人影院| 内射极品少妇av片p| 深夜a级毛片| 色播亚洲综合网| 高清在线国产一区| 亚洲人成网站在线播放欧美日韩| 男女做爰动态图高潮gif福利片| 男女视频在线观看网站免费| 日本三级黄在线观看| 九九在线视频观看精品| 国产成年人精品一区二区| 久久久国产成人免费| 99热精品在线国产| 久久热精品热| 黄色视频,在线免费观看| 女人被狂操c到高潮| 日韩国内少妇激情av| 亚洲美女黄片视频| 最近中文字幕高清免费大全6 | 深夜精品福利| 成人av在线播放网站| 欧美最黄视频在线播放免费| 午夜福利高清视频| 亚洲无线在线观看| 色视频www国产| 成人永久免费在线观看视频| 亚洲国产精品sss在线观看| 亚洲av熟女| 日韩欧美在线二视频| 日韩中文字幕欧美一区二区| 少妇丰满av| 97热精品久久久久久| 十八禁国产超污无遮挡网站| 我的老师免费观看完整版| 哪里可以看免费的av片| 无人区码免费观看不卡| 18+在线观看网站| 97人妻精品一区二区三区麻豆| 日本一本二区三区精品| 成年版毛片免费区| 国产91精品成人一区二区三区| 国产欧美日韩一区二区精品| 国产精品野战在线观看| 亚洲精品日韩av片在线观看| 97热精品久久久久久| 99久久精品一区二区三区| 91字幕亚洲| 在线观看午夜福利视频| 久久国产精品影院| 成人午夜高清在线视频| 757午夜福利合集在线观看| 舔av片在线| 国产激情偷乱视频一区二区| 国产精品日韩av在线免费观看| 国产aⅴ精品一区二区三区波| 一级a爱片免费观看的视频| 夜夜夜夜夜久久久久| 亚洲中文字幕一区二区三区有码在线看| 久久午夜福利片| 成人高潮视频无遮挡免费网站| 国产高清视频在线观看网站| 免费在线观看成人毛片| 伦理电影大哥的女人| 国产高清三级在线| 色av中文字幕| 久久精品人妻少妇| 欧美成人一区二区免费高清观看| 久久人妻av系列| 又粗又爽又猛毛片免费看| 色精品久久人妻99蜜桃| 国产在线男女| 在线免费观看的www视频| 国产欧美日韩一区二区精品| 亚洲,欧美,日韩| 五月伊人婷婷丁香| 欧美成人一区二区免费高清观看| 亚洲国产精品久久男人天堂| 麻豆久久精品国产亚洲av| 亚洲aⅴ乱码一区二区在线播放| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲电影在线观看av| 香蕉av资源在线| 国产人妻一区二区三区在| 一本综合久久免费| 丰满的人妻完整版| 国产欧美日韩一区二区精品| 亚洲成av人片在线播放无| 国产蜜桃级精品一区二区三区| av女优亚洲男人天堂| 伦理电影大哥的女人| 亚洲精品粉嫩美女一区| 99久国产av精品| 日本黄大片高清| 一级毛片久久久久久久久女| 九色国产91popny在线| 久久人人精品亚洲av| 欧美乱妇无乱码| 日本熟妇午夜| 亚洲av成人av| 国产成人欧美在线观看| 日本黄大片高清| 国产高清有码在线观看视频| 啦啦啦观看免费观看视频高清| 亚洲国产精品999在线| 一区二区三区免费毛片| 脱女人内裤的视频| 亚洲av日韩精品久久久久久密| 美女被艹到高潮喷水动态| 九色成人免费人妻av| 日日摸夜夜添夜夜添小说| 精品久久国产蜜桃| 日日摸夜夜添夜夜添av毛片 | 午夜免费激情av| 亚洲国产高清在线一区二区三| 亚洲欧美清纯卡通| 两个人的视频大全免费| 非洲黑人性xxxx精品又粗又长| 男人舔奶头视频| 国内精品久久久久久久电影| 91av网一区二区| 免费看日本二区| 欧美bdsm另类| 精品人妻一区二区三区麻豆 | 午夜两性在线视频| 国产欧美日韩一区二区精品| 亚洲av成人不卡在线观看播放网| 久久99热6这里只有精品| 村上凉子中文字幕在线| 亚洲欧美精品综合久久99| 午夜a级毛片| 日本三级黄在线观看| 国产私拍福利视频在线观看| 国产精品电影一区二区三区| 乱人视频在线观看| 黄色一级大片看看| 久久久久九九精品影院| 国产精品一区二区性色av| 一级作爱视频免费观看| 婷婷精品国产亚洲av在线| 两人在一起打扑克的视频| 亚洲成人精品中文字幕电影| 国产精品久久电影中文字幕| av国产免费在线观看| 午夜激情福利司机影院| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美日韩卡通动漫| 久久久久久久久久成人| 淫秽高清视频在线观看| 亚洲三级黄色毛片| 亚洲av免费在线观看| 欧美另类亚洲清纯唯美| 久久国产乱子免费精品| 免费在线观看亚洲国产| 日韩欧美 国产精品| 18美女黄网站色大片免费观看| 日本黄色片子视频| 每晚都被弄得嗷嗷叫到高潮| 禁无遮挡网站| 国产老妇女一区| а√天堂www在线а√下载| 中文字幕高清在线视频| 女生性感内裤真人,穿戴方法视频| 乱码一卡2卡4卡精品| 久久精品国产自在天天线| aaaaa片日本免费| 久久久久久九九精品二区国产| 国模一区二区三区四区视频| 99久久成人亚洲精品观看| 国产91精品成人一区二区三区| 啦啦啦观看免费观看视频高清| 九九热线精品视视频播放| 好看av亚洲va欧美ⅴa在| 女人十人毛片免费观看3o分钟| 在线a可以看的网站| 久9热在线精品视频| 日韩欧美国产一区二区入口| 婷婷丁香在线五月| 在线免费观看的www视频| 看片在线看免费视频| 国内毛片毛片毛片毛片毛片| 一级黄色大片毛片| 国产真实乱freesex| 色在线成人网| 中文字幕人妻熟人妻熟丝袜美| 亚洲av成人av| 亚洲,欧美,日韩| 国产精品久久视频播放| 午夜福利高清视频| 亚洲欧美精品综合久久99| 国产私拍福利视频在线观看| 午夜免费男女啪啪视频观看 | 亚洲精品粉嫩美女一区| 久久精品国产亚洲av涩爱 | 日本免费一区二区三区高清不卡| 国产视频内射| 精品国产三级普通话版| 免费人成在线观看视频色| 少妇裸体淫交视频免费看高清| 成年版毛片免费区| 成人av在线播放网站| 国产真实伦视频高清在线观看 | 在线观看66精品国产| 99精品久久久久人妻精品| 日韩成人在线观看一区二区三区| 又爽又黄a免费视频| 国产一区二区亚洲精品在线观看| 欧美精品啪啪一区二区三区| 国产精华一区二区三区| 亚洲真实伦在线观看| 久久精品国产亚洲av天美| 变态另类丝袜制服| 两性午夜刺激爽爽歪歪视频在线观看| 一边摸一边抽搐一进一小说| 日韩精品青青久久久久久| 精品99又大又爽又粗少妇毛片 | 真人一进一出gif抽搐免费| 人人妻,人人澡人人爽秒播| 国产精品久久电影中文字幕| 精品久久国产蜜桃| 老熟妇乱子伦视频在线观看| 国产人妻一区二区三区在| 美女免费视频网站| 午夜精品一区二区三区免费看| 俄罗斯特黄特色一大片| 又紧又爽又黄一区二区| 亚洲人成伊人成综合网2020| 丰满的人妻完整版| 99久久99久久久精品蜜桃| 亚洲五月天丁香| 老熟妇仑乱视频hdxx| 老司机福利观看| 天堂√8在线中文| 国产在线男女| 黄色视频,在线免费观看| 久久精品国产亚洲av涩爱 | 老司机福利观看| 国产不卡一卡二| 中文字幕免费在线视频6| 亚洲精品成人久久久久久| 色视频www国产| 亚洲美女黄片视频| 国产人妻一区二区三区在| 国产爱豆传媒在线观看| 简卡轻食公司| 91久久精品电影网| 国产美女午夜福利| 日韩人妻高清精品专区| 一级作爱视频免费观看| 人妻夜夜爽99麻豆av| 简卡轻食公司| 91字幕亚洲| 久久久色成人| 欧美不卡视频在线免费观看| 丰满乱子伦码专区| 一本精品99久久精品77| 久久国产乱子免费精品| 成人av在线播放网站| 国产精品1区2区在线观看.| 午夜福利免费观看在线| 日本一本二区三区精品| 成人三级黄色视频| 国产亚洲精品久久久com| 夜夜躁狠狠躁天天躁| 国产一区二区在线观看日韩| 国产麻豆成人av免费视频| 国产精品1区2区在线观看.| 最近视频中文字幕2019在线8| 日韩欧美国产在线观看| 亚洲成人精品中文字幕电影| 国产欧美日韩一区二区精品| 精品福利观看| 深爱激情五月婷婷| av专区在线播放| 一进一出好大好爽视频| 久久草成人影院| 此物有八面人人有两片| 免费观看的影片在线观看| 成人国产综合亚洲| 乱人视频在线观看| 99精品在免费线老司机午夜| 国产蜜桃级精品一区二区三区| 老司机福利观看| 亚洲综合色惰| 国内精品久久久久精免费| 亚洲av日韩精品久久久久久密| 在线观看舔阴道视频| 国产亚洲av嫩草精品影院| 无遮挡黄片免费观看| 十八禁国产超污无遮挡网站| 51午夜福利影视在线观看| 麻豆久久精品国产亚洲av| 欧美另类亚洲清纯唯美| 九九热线精品视视频播放| 九九在线视频观看精品| 又爽又黄无遮挡网站| 99在线人妻在线中文字幕| 精品福利观看| 欧美色欧美亚洲另类二区| 国产精品,欧美在线| 久久久久久久久久成人| 麻豆成人午夜福利视频| 国产精品av视频在线免费观看| 欧美高清成人免费视频www| 日本一二三区视频观看| 中文字幕人成人乱码亚洲影| 亚洲中文日韩欧美视频| 午夜老司机福利剧场| 亚洲精品日韩av片在线观看| 精品国内亚洲2022精品成人| 毛片一级片免费看久久久久 | 国产一区二区在线观看日韩| 国产激情偷乱视频一区二区| www.www免费av| 搡老岳熟女国产| 真人做人爱边吃奶动态| 18美女黄网站色大片免费观看| 亚洲乱码一区二区免费版| 亚洲在线观看片| 亚洲av第一区精品v没综合| 大型黄色视频在线免费观看| 99久久成人亚洲精品观看| 日本免费a在线| 精品人妻偷拍中文字幕| 国产又黄又爽又无遮挡在线| 一区二区三区四区激情视频 | 99久国产av精品| 丰满乱子伦码专区| 亚洲国产高清在线一区二区三| 人妻久久中文字幕网|