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

    雙饋風機附加頻率控制對系統(tǒng)調頻動態(tài)的影響

    2021-10-20 02:12:16張政彭曉濤李少林梁愷許饒琪張銳
    電力建設 2021年10期
    關鍵詞:雙饋調頻特征值

    張政,彭曉濤,李少林,梁愷,許饒琪,張銳

    (1.武漢大學電氣與自動化學院,武漢市430072;2.中國電力科學研究院有限公司,北京市 100192;3.合肥陽光智維科技有限公司,合肥市 230088)

    0 引 言

    雙饋風電機組(doubly-fed induction generator based wind turbine,DFIG-WT)由于采用最大功率跟蹤控制(maximum power point tracking,MPPT)模式,其轉子轉速與電網(wǎng)頻率解耦,不能為系統(tǒng)提供類似于同步機慣性響應特性的動態(tài)頻率支撐,因此大規(guī)模風電并網(wǎng)會削弱系統(tǒng)慣性,使系統(tǒng)在負荷擾動后產生更大的頻率波動[1-2]。為改善雙饋風機有功頻率響應特性,使風機參與到系統(tǒng)頻率控制中,雙饋風機的附加頻率控制策略得到了廣泛的研究[3-5],其中,利用組合頻率微分和偏差信號實現(xiàn)風機慣性支撐和下垂控制的附加頻率控制,成為改善風電系統(tǒng)頻率響應特性的重要手段。但是這些附加控制在使風機具有頻率支撐能力的同時,也對系統(tǒng)動態(tài)特性產生了一定的影響[6-7]。風機與系統(tǒng)通過有功調頻產生相互作用的通道,從而可能使系統(tǒng)調頻動態(tài)環(huán)節(jié)與雙饋風電機組的多物理控制環(huán)節(jié)發(fā)生耦合作用,進而對系統(tǒng)負荷-頻率控制的動態(tài)特性產生影響[8-9]。因此,有必要對雙饋風機附加頻率控制給系統(tǒng)調頻動態(tài)帶來的影響進行進一步研究。

    近年來,已有很多文獻圍繞利用小干擾分析方法針對風機附加頻率控制對系統(tǒng)小干擾穩(wěn)定性的影響進行了研究。文獻[10]采用阻尼轉矩分析法,分析了雙饋風機采用基于頻率偏差的下垂控制對多機風電并網(wǎng)系統(tǒng)小干擾功角振蕩的影響,以及雙饋風機接入位置和同步機慣性時間常數(shù)與下垂控制對系統(tǒng)小擾動功角穩(wěn)定影響特性的關系。文獻[11-14]則研究了風電虛擬慣量控制對風電并網(wǎng)多機系統(tǒng)小干擾穩(wěn)定性的影響。其中,文獻[11]通過理論推導和基于DigSILENT仿真的特征值模態(tài)分析,研究了工作在MPPT模式下的雙饋風電機組采用附加虛擬慣性控制對四機兩區(qū)系統(tǒng)機電振蕩模式的影響;文獻[12]通過時域仿真,借助模式分析法研究了雙饋風電機組采用基于頻率微分反饋和基于變流器控制環(huán)節(jié)時間常數(shù)調整的附加虛擬慣量控制,以及控制策略不同控制參數(shù)對風電并網(wǎng)區(qū)域系統(tǒng)的發(fā)電機主導振蕩模態(tài)小干擾穩(wěn)定的影響;文獻[13]研究了雙饋風電機組鎖相環(huán)動態(tài)與虛擬慣量控制對風機并網(wǎng)四機兩區(qū)系統(tǒng)區(qū)間振蕩模態(tài)的影響;文獻[14]通過建立考慮鎖相環(huán)動態(tài)的直驅風電機組虛擬慣性控制模型,針對鎖相環(huán)振蕩模態(tài)與系統(tǒng)機電模式在復平面分布位置的不同情況,利用特征值分析研究了虛擬慣性控制對系統(tǒng)機電振蕩模式動態(tài)特性的影響。但是,用于這些影響研究的系統(tǒng)對象都沒有考慮互聯(lián)系統(tǒng)的調頻動態(tài)特性,并且在反映風電機組多物理控制環(huán)節(jié)對系統(tǒng)小干擾穩(wěn)定動態(tài)特性的耦合關系方面,由于考慮側重點不同,并未對風電機組動態(tài)進行詳細建模。此外文獻[15]通過建立一次調頻的直驅永磁風電系統(tǒng)小信號模型,基于特征值分析研究了不同風速下風電參與系統(tǒng)一次調頻的小擾動穩(wěn)定性,但也未研究對多機多區(qū)系統(tǒng)調頻動態(tài)的影響。

    基于上述分析,從分析雙饋風機附加調頻控制對互聯(lián)系統(tǒng)負荷-頻率控制動態(tài)特性的影響出發(fā),同時考慮風電并網(wǎng)電力系統(tǒng)的運動模態(tài)通常由多動態(tài)環(huán)節(jié)耦合作用而成的特點[16-17],本文通過建立詳細描述雙饋風機多物理控制環(huán)節(jié)耦合特性的狀態(tài)空間模型,以及描述互聯(lián)系統(tǒng)二次調頻動態(tài)特性的數(shù)學模型,利用特征值模態(tài)和參與因子分析方法,研究風機附加調頻控制對風機參與系統(tǒng)調頻的小干擾穩(wěn)定動態(tài)特性影響,并利用時域仿真驗證理論分析結果的合理性。

    1 考慮附加頻率控制的雙饋風機數(shù)學模型

    1.1 雙饋風電機組數(shù)學模型

    DFIG-WT是利用雙饋感應發(fā)電機(doubly-fed induction generator,DFIG)將風力機捕獲的風能轉換為電能的機組,與傳統(tǒng)同步發(fā)電機不同的是,DFIG除定子與電網(wǎng)相連外,轉子也經(jīng)背靠背變流器與電網(wǎng)相連,其結構如圖1所示。

    圖1 雙饋風電機組結構圖Fig.1 Structure diagram of a DFIG-WT

    DFIG-WT由DFIG、背靠背變流器及其直流環(huán)節(jié)、濾波電路、軸系等多動態(tài)物理環(huán)節(jié),以及轉子側變流器(rotor side converter,RSC)控制、網(wǎng)側變流器(grid side converter,GSC)控制、直流環(huán)節(jié)控制、鎖相環(huán)控制、轉速控制等多時間尺度控制環(huán)節(jié)構成[18]。各物理控制環(huán)節(jié)的數(shù)學模型如圖1所示。圖中,Tm、Te和Teref分別表示風電機組傳動鏈軸系的機械轉矩、電磁轉矩及其參考值;ωrref和ωr分別表示發(fā)電機轉子轉速參考值和實際值;Ps和Psref分別表示定子輸出有功及其指令值;Qsref表示定子輸出無功指令值;ir和is分別表示發(fā)電機轉子和定子繞組電流;ig表示GSC交流側輸入電流;Udc和Cdc分別表示變流器直流側端電壓和濾波電容;Lg和Rg分別表示GSC交流側濾波電感和電阻;ur表示發(fā)電機轉子繞組端電壓;usd和usq分別表示風電機組并網(wǎng)點電壓d、q軸分量;ωPLL和θPLL分別表示并網(wǎng)點電壓的角頻率和相位。

    1)DFIG數(shù)學模型。

    DFIG數(shù)學模型由同步旋轉dq坐標系的電壓方程和磁鏈方程組成??紤]感應電機的三相平衡,并不考慮磁鏈飽和效應和零序分量影響,將磁鏈方程代入電壓方程,經(jīng)整理可建立發(fā)電機四階微分方程[19],如式(1)所示。

    (1)

    式中:urd、urq和ird、irq分別表示RSC交流側輸出電壓和電流的d、q軸分量;usd、usq和isd、isq分別表示發(fā)電機定子繞組電壓和電流的d、q軸分量;Rs、Rr分別表示發(fā)電機定子和轉子的電阻;Xss=Xs+Xm;Xrr=Xr+Xm;Xs表示定子繞組感抗;Xr表示轉子繞組感抗;Xm表示定子和轉子繞組間的互感抗;ωb和ωs分別表示基準角頻率和系統(tǒng)同步電氣角頻率,且ωs=ω/ωb,當參考坐標系的旋轉速度ω設置為ωb時,同步電氣角頻率與基準角頻率相等;s為ωr與ωs之間的轉差率。

    2)變流器控制數(shù)學模型。

    為實現(xiàn)雙饋風電機組的變速恒頻并網(wǎng)運行,DFIG轉子需采用交-直-交變流器實現(xiàn)交流勵磁控制。RSC根據(jù)DFIG輸出有功和無功指令參考值Psref、Qsref,經(jīng)比例-積分(PI)控制確定轉子d、q軸電流分量參考值irdref、irqref,利用此電流參考指令經(jīng)PI控制產生用于PWM控制的調制電壓指令urdref和urqref,在為轉子勵磁繞組提供變頻三相交流電壓時,也分別利用q、d軸電壓分量實現(xiàn)定子輸出有功和無功的解耦控制。GSC利用直流環(huán)節(jié)端電壓控制偏差,經(jīng)PI控制產生交流側輸入電流d軸分量參考值igdref,利用定子端電壓變化量產生輸入電流q軸分量參考值igqref,再利用此2個電流參考值,經(jīng)PI控制產生用于PWM控制的調制電壓信號ugdref和ugqref,實現(xiàn)直流環(huán)節(jié)電容端電壓恒定控制的同時,并可以利用無功調節(jié)阻尼定子端電壓變化。上述RSC和GSC控制過程的數(shù)學模型分別如式(2)和式(3)所示。

    (2)

    (3)

    3)軸系數(shù)學模型。

    軸系數(shù)學模型用于描述連接風力機和感應發(fā)電機的轉子角速度動態(tài)變化過程,通??刹捎脙少|量塊進行建模,其數(shù)學模型如式(4)所示。

    (4)

    式中:ωt為風力機轉子機械角速度;ωr為發(fā)電機的轉子轉速;θ為軸系扭轉角度;Ht為風力機的慣性質量常數(shù);Hg為發(fā)電機的慣性質量常數(shù);風力機和發(fā)電機的轉子之間通過剛性系數(shù)為K、阻尼為D的彈簧軸連接。

    4)轉速控制數(shù)學模型。

    轉速控制用于產生定子輸出功率指令,實現(xiàn)風電機組功率跟蹤運行,其動態(tài)數(shù)學方程如式(5)所示。

    (5)

    5)直流環(huán)節(jié)數(shù)學模型。

    直流環(huán)節(jié)數(shù)學模型主要用于描述電容端電壓的動態(tài)變化過程,其數(shù)學模型如式(6)所示。

    式中:ugd、ugq和igd、igq分別為GSC交流側輸入電壓和電流的d、q軸分量;urd、urq和ird、irq分別為RSC交流側輸出電壓和電流的d、q軸分量。

    6)濾波電路數(shù)學模型。

    濾波電路數(shù)學模型用于描述GSC動態(tài)調節(jié)輸入功率時,GSC交流輸入電流的動態(tài)變化過程,其數(shù)學模型如式(7)所示。

    (7)

    7)鎖相環(huán)數(shù)學模型。

    鎖相環(huán)主要用于通過跟蹤端電壓相位來實現(xiàn)變流器與交流電網(wǎng)的同步,其同步dq坐標下的數(shù)學模型如式(8)所示[6]。

    (8)

    式中:kp7和ki7是鎖相環(huán)的PI控制參數(shù)。根據(jù)由鎖相環(huán)得到的相角θPLL,可將雙饋風機的端電壓由xy坐標系變換到dq坐標系,坐標變換的公式如式(9)所示。

    (9)

    式中:usx和usy是以αβ靜止坐標系表示的風機并網(wǎng)點電壓。

    綜合上述微分方程可知,DFIG-WT模型包括20個狀態(tài)變量,分別為[isd,isq,ird,irq,x1,x2,x3,x4,x5,x6,x7,θ,ωt,ωr,x8,UDC,igd,igq,ωPLL,θPLL]。

    現(xiàn)代大型風電場包含數(shù)量眾多的風電機組,在研究中若采用包含每臺風機及其控制模塊的詳細模型,不僅會増加運算的復雜程度,而且在利用特征值法進行小干擾穩(wěn)定分析時可能會出現(xiàn)“維數(shù)災”的問題?,F(xiàn)有研究對這一問題的解決方法是先將風電場聚合等效為1臺或少數(shù)幾臺等值風電機組,同時根據(jù)研究需要對模型的某些部分進行簡化處理。本文同樣利用聚合法將風電場中多臺風機用1臺等值風電機組代替,同時由于模型總體階數(shù)不算太高,出于對模型準確性和計算精確度的考慮,未對風機進行進一步簡化處理。

    1.2 風電機組附加調頻控制數(shù)學模型

    由于雙饋風機利用電力電子裝置并網(wǎng)運行時,其功率跟蹤運行方式使轉子動能與電網(wǎng)頻率變化解耦,無法為頻率擾動提供轉動慣量支撐和頻率調節(jié)響應[3]。隨著同步發(fā)電機被風電并網(wǎng)替代而退出運行,這種解耦也使電力系統(tǒng)的同步旋轉慣量減少,對頻率擾動時的慣量支撐能力減弱,因此風電機組的調頻控制研究也日益受到關注。利用頻率微分和頻率偏差實現(xiàn)雙饋風機慣性支撐和下垂控制的組合頻率控制是一種典型附加頻率控制方案[5],如圖2所示。該控制策略可使雙饋風機在2種頻率支撐時間尺度上發(fā)揮調頻作用,從而改善風電并網(wǎng)系統(tǒng)的頻率響應特性。

    圖2 雙饋風電機組附加頻率控制原理圖Fig.2 Additional frequency control scheme of DFIG-WT

    如圖2所示,附加頻率控制策略是根據(jù)電網(wǎng)頻率f及其參考值fref的偏差在風電機組功率跟蹤控制指令PsPPT上疊加慣量響應有功指令Pinertia和調頻響應有功指令Pdroop,使風機在頻率變化時能利用變槳減載、超速減載、轉子動能等有功儲備進行頻率支撐的有功調節(jié)[20]。其中,Pinertia由頻率偏差微分信號經(jīng)比例環(huán)節(jié)kd確定,即:

    (10)

    Pdroop由頻率偏差信號經(jīng)下垂控制比例環(huán)節(jié)kp確定,即:

    Pdroop=kp(f-fref)

    (11)

    數(shù)學模型表明附加頻率控制給系統(tǒng)動態(tài)特性帶來的影響主要體現(xiàn)在慣性響應的微分環(huán)節(jié)增加了系統(tǒng)的階數(shù),并通過引入新狀態(tài)變量Δf而產生一個新的模態(tài)。同時,由于系統(tǒng)狀態(tài)空間描述的輸入系數(shù)矩陣元素值與狀態(tài)變量穩(wěn)態(tài)值有關,采用附加頻率控制后,風機的有功指令值在頻率擾動時與不采用附加頻率控制時不同,因此狀態(tài)變量穩(wěn)態(tài)值的不同也會使采用附加頻率控制前后的狀態(tài)系數(shù)矩陣元素值不同,從而使求解的特征值也不同。

    2 無附加調頻的系統(tǒng)調頻動態(tài)特性

    為進一步分析頻率擾動下的風電并網(wǎng)系統(tǒng)的調頻動態(tài)特性,以及風機附加頻率控制對其動態(tài)特性的影響,基于1.1節(jié)討論的雙饋風機數(shù)學模型和等值風電機組建模方法[21],采用MATLAB/Simulink建立含雙饋風機參與互聯(lián)電力系統(tǒng)負荷頻率控制調頻動態(tài)數(shù)學模型,如圖3所示。

    圖3 含風電互聯(lián)電力系統(tǒng)調頻控制原理圖Fig.3 Frequency-regulation scheme of interconnected power system with wind power generator

    模型中TG和TCH分別表示調速器和原動機一階慣性響應時間常數(shù);M表示同步發(fā)電機轉子角動量,等于轉子轉動慣量J乘以同步電氣角速度ω0;KD表示負荷的頻率調節(jié)效應系數(shù);B1、B2和R1、R2分別表示系統(tǒng)1、2的頻率調節(jié)效應系數(shù)和等值發(fā)電機的頻率調差系數(shù);K和T都表示比例放大系數(shù)。由圖3可見,依據(jù)前述風電機組狀態(tài)空間模型建立的等值雙饋風電機組接入互聯(lián)系統(tǒng)1時,由于互聯(lián)電力系統(tǒng)包含x21—x30所示的10個狀態(tài)變量,因此圖3所示含風電互聯(lián)電力系統(tǒng)共包含30個狀態(tài)變量,如表1所示。

    圖3中等值同步發(fā)電機和等值風電機組出力分別設為500 MW和300 MW,等值負荷設為1 300 MW,互聯(lián)系統(tǒng)調頻動態(tài)數(shù)學模型的相關參數(shù)和雙饋風機參數(shù)參見附錄A。采用MATLAB/Simulink的Control design toolbox對仿真模型在穩(wěn)態(tài)運行時進行線性化,基于所得系統(tǒng)狀態(tài)系數(shù)矩陣求解特征值,并利用復數(shù)特征值計算相應振蕩頻率和阻尼比,所得系統(tǒng)調頻動態(tài)的運動模態(tài)分析結果如表2所示。結果表明,系統(tǒng)共有7對共軛復特征根和16個實特征根,分別對應7個振蕩模態(tài)和16個非振蕩模態(tài)。所有特征值均位于虛軸左側,表明系統(tǒng)調頻動態(tài)具有較好小干擾穩(wěn)定性。

    表2 無附加頻率控制時系統(tǒng)調頻動態(tài)的模態(tài)分析結果Table 2 Modal analysis results of frequency-regulation dynamic of power system without additional frequency control

    進一步對各振蕩模態(tài)的參與因子進行分析,所得結果如圖4所示。

    圖4 振蕩模態(tài)參與因子分析Fig.4 Participation factor analysis of oscillation modes

    由圖4中振蕩模態(tài)的參與因子分析可見,λ1,2所確定模態(tài)的振蕩頻率是51.768 Hz,雙饋感應發(fā)電機的狀態(tài)變量[isdisqirdirq]是最大參與因子,表明該振蕩模態(tài)由感應發(fā)電機轉子和定子電流控制的相互耦合作用產生。λ3,4所確定模態(tài)的振蕩頻率是9.639 Hz,狀態(tài)變量[ωPLLδPLL]的參與因子分析表明,該模態(tài)由鎖相環(huán)動態(tài)主導。λ5,6和λ7,8的振蕩頻率分別是0.42 Hz和1.145 Hz,根據(jù)參與因子分析可知,2個模態(tài)分別是由RSC的無功和有功外環(huán)控制與感應發(fā)電機相互作用主導的低頻振蕩。λ9,10所確定模態(tài)的振蕩頻率是3.182 Hz,僅有軸系部分的3個狀態(tài)變量[θωtωr]參與,屬于軸系振蕩。λ11,12所確定模態(tài)的振蕩頻率是0.324 Hz,狀態(tài)變量[x5x6igdUdc] 的參與因子分析表明,該模態(tài)由濾波環(huán)節(jié)、直流環(huán)節(jié)和GSC相互作用的直流電壓恒定控制產生。λ13,14所確定模態(tài)的振蕩頻率是0.016 Hz,是由互聯(lián)系統(tǒng)調頻動態(tài)相互耦合作用產生。并且風機與系統(tǒng)振蕩模態(tài)的耦合作用環(huán)節(jié)相互獨立,λ1—λ12所確定的6種振蕩模態(tài)僅與雙饋風機的狀態(tài)變量相關,互聯(lián)系統(tǒng)負荷-頻率控制的振蕩模態(tài)也僅與調頻動態(tài)環(huán)節(jié)的狀態(tài)變量相關,這也表明無附加頻率控制的雙饋風機接入系統(tǒng),不會對互聯(lián)系統(tǒng)負荷-頻率控制的振蕩模態(tài)產生影響,并且有功輸出控制與系統(tǒng)頻率變化解耦。對比各振蕩模態(tài)的阻尼比可知,系統(tǒng)調頻控制的小干擾穩(wěn)定動態(tài)特性優(yōu)于雙饋風機的小干擾穩(wěn)定動態(tài)特性。

    系統(tǒng)16個非振蕩模態(tài)的參與因子分析結果如表3所示。

    表3 非振蕩模態(tài)的參與因子分析結果Table 3 Participation factor analysis of non-oscillation modes

    表3中,λ17、λ18、λ22、λ23、λ26和λ29所確定非振蕩模態(tài)是僅由互聯(lián)系統(tǒng)調頻動態(tài)環(huán)節(jié)主導的模態(tài),而其他非振蕩模態(tài)都只與雙饋風電機組的狀態(tài)變量相關。分析結果同樣表明,無附加頻率控制的雙饋風機沒有對系統(tǒng)負荷-頻率控制的非振蕩模態(tài)產生影響。

    根據(jù)各模態(tài)的參與因子分析的結果,總結了DFIG各組成物理控制環(huán)節(jié)耦合作用的調頻運動模態(tài),如表4所示。結果表明,感應電機與其他組成部分相比,在振蕩模態(tài)和非振蕩模態(tài)中的參與程度均較高,分別與51.768、0.42、1.145 Hz的振蕩模態(tài)和2個非振蕩模態(tài)存在耦合作用。因此在DFIG各組成物理控制環(huán)節(jié)中,感應電機對系統(tǒng)調頻動態(tài)的影響作用最大。

    表4 DFIG各組成環(huán)節(jié)對系統(tǒng)調頻動態(tài)耦合作用的情況對比Table 4 Comparison of the dynamic coupling effect of each DFIG component on the system frequency regulation

    3 有風電調頻的系統(tǒng)調頻動態(tài)特性

    3.1 附加調頻控制對調頻動態(tài)的影響

    將圖3所示仿真模型的DFIG-WT采用圖2所示附加頻率控制,并設5 s時刻等值發(fā)電機1的有功負荷增加100 MW,雙饋風機采用附加頻率控制前后系統(tǒng)頻率和風機輸出有功的動態(tài)響應如圖5和圖6所示。

    圖5 DFIG-WT采用附加頻率控制前后的系統(tǒng)頻率響應Fig.5 System frequency response when DFIG-WT with or without additional frequency control

    圖6 DFIG-WT采用附加頻率控前后的有功輸出Fig.6 Output power of DFIG-WT with or without additional frequency control

    對比結果驗證了附加頻率控制的有效性。圖5的對比表明,采用附加頻率控制后,系統(tǒng)的頻率最低點從48.44 Hz提高到49.62 Hz,系統(tǒng)頻率恢復穩(wěn)態(tài)的動態(tài)過程對比說明附加控制在提高系統(tǒng)頻率支撐能力的同時,也能夠優(yōu)化系統(tǒng)的頻率調節(jié)特性。圖6的對比則說明,附加頻率控制能使雙饋風機在頻率擾動期間利用附加調頻有功指令Pinertia和Pdroop增加有功輸出,以提供慣量支撐和調頻響應。

    上述系統(tǒng)中,由于附加頻率慣量支撐使系統(tǒng)引入了頻率f作為新增狀態(tài)變量x31,因此基于特征值分析得到的運動模態(tài)也增加為31個,如表5所示。

    表5 有附加頻率控制時系統(tǒng)調頻動態(tài)的模態(tài)分析Table 5 Mode analysis results of frequency-regulation dynamic of power system with additional frequency control

    與表2相比,風機附加頻率控制的引入使系統(tǒng)負荷-頻率控制動態(tài)增加一個特征值λ31=-94.104,并且此特征值在負實軸上離虛軸較遠,其所對應非振蕩模態(tài)將不會對系統(tǒng)負荷-頻率控制的小干擾穩(wěn)定動態(tài)特性產生影響。對比其他特征值,采用附加頻率控制后所計算的特征值都在未采用附加頻率控制時的特征值附近,說明附加頻率控制雖然對系統(tǒng)各模態(tài)的小干擾穩(wěn)定動態(tài)特性產生影響,但影響程度不明顯。

    對各模態(tài)進行參與因子分析可確定采用附加頻率控制前后,各模態(tài)的耦合作用環(huán)節(jié)的變化。通過對比2.2節(jié)中各模態(tài)參與因子的變化可知,λ31、λ7,8和λ17所確定模態(tài)的參與因子發(fā)生變化。采用附加頻率控制后,上述3種模態(tài)參與因子的分析結果如圖7所示。其中,λ31是采用附加頻率控制后新增的模態(tài),由模態(tài)的耦合狀態(tài)變量可知,該模態(tài)由慣量支撐控制、系統(tǒng)1的二次調頻、雙饋感應發(fā)電機定子轉子有功電流分量控制相互耦合作用產生。

    圖7 λ31、λ7,8和λ17的參與因子分析Fig.7 Participation factor analysis of λ31,λ7,8 and λ17

    對比圖4和圖7中λ7,8所確定模態(tài)的參與因子變化可見,附加頻率控制使系統(tǒng)調頻動態(tài)環(huán)節(jié)與風電機組動態(tài)環(huán)節(jié)發(fā)生耦合作用,對比阻尼比變化可知,該耦合作用使雙饋風機的該振蕩模態(tài)的小干擾穩(wěn)定動態(tài)特性變差。而λ17所確定非振蕩模態(tài)的參與因子表明,x31也具有一定的參與作用,并且該狀態(tài)變量使雙饋風機的感應發(fā)電機和RSC控制環(huán)節(jié)參與耦合作用,結合特征值在負實軸的位置變化可知,附加頻率控制提高了該模態(tài)的動態(tài)穩(wěn)定性。在其他模態(tài)中,因x31的參與因子接近于0,使附加頻率控制在其他模態(tài)中沒有使雙饋風機的物理控制環(huán)節(jié)與系統(tǒng)調頻動態(tài)發(fā)生耦合作用。圖7仿真結果也同時表明,附加頻率控制使風機對與其電氣聯(lián)系較強區(qū)域系統(tǒng)的調頻動態(tài)產生影響,對電氣聯(lián)系較弱系統(tǒng)的調頻動態(tài)則影響不大。

    3.2 附加調頻控制參數(shù)對動態(tài)特性的影響

    在附加頻率控制中,kd和kp這2個參數(shù)不僅會影響慣性支撐和下垂控制的作用效果,而且對風機動態(tài)特性也會產生影響,為此,基于參數(shù)攝動法進一步研究2個參數(shù)變化對調頻動態(tài)特性的影響。

    分別改變kd和kp研究特征值變化情況,小信號穩(wěn)定性分析表明,系統(tǒng)特征值在參數(shù)大范圍變化時變化情況較為復雜,不僅不同模態(tài)對參數(shù)變化的敏感度不同,而且不同模態(tài)的穩(wěn)定性變化也不同。圖8分別對比了模態(tài)7、8和模態(tài)1、2在保持kp=10、kd=1,振蕩頻率和阻尼比分別對kd和kp變化的響應。由圖8可見,模態(tài)7、8的阻尼比在kd增加到12.04時減小為0,模態(tài)1、2的阻尼比在kp增加到1 327時減小至0,系統(tǒng)在此2種參數(shù)作用下達到臨界穩(wěn)定狀態(tài)。分析結果也表明,附加調頻控制參數(shù)的變化范圍存在使系統(tǒng)保持穩(wěn)定的限制。

    圖8 臨界穩(wěn)定附近的模態(tài)阻尼比和振蕩頻率變化Fig.8 Modal damping ratio and oscillation frequency variation near critical stability

    在參數(shù)限定范圍內,進一步研究各模態(tài)振蕩頻率和阻尼比對kd和kp小范圍變化的響應靈敏度,以明確系統(tǒng)調頻動態(tài)特性在參數(shù)設定值附近的穩(wěn)定性。

    保持kp=10,令kd以0.1的步長從0.1增加到1,系統(tǒng)7個振蕩模態(tài)的振蕩頻率和阻尼比變化情況如圖9所示。

    圖9 λ1-14振蕩頻率和阻尼比隨kd的變化Fig.9 Affection of oscillation frequency and damping ratio of λ1-14 with kd change

    圖9中,λ1,2和λ13,14所確定振蕩模態(tài)的振蕩頻率和阻尼比隨kd增大而分別增大和減小。λ5—λ8所確定振蕩模態(tài)的振蕩頻率和阻尼比隨kd增大而分別減小和增大。λ3,4所確定振蕩模態(tài)的動態(tài)特性不隨kd變化而變化,λ9—λ12所確定振蕩模態(tài)的動態(tài)特性也對kd變化不敏感。

    對于非振蕩模態(tài),選取僅與互聯(lián)系統(tǒng)調頻動態(tài)環(huán)節(jié)相關的6個模態(tài),研究kd變化時附加頻率控制對系統(tǒng)調頻動態(tài)特性的影響。相關模態(tài)特征值變化情況如圖10所示,圖中紅色和藍色叉號分別表示加入附加頻率控制前后的特征值。其中λ17、λ18、λ22、λ23、λ26所確定的模態(tài),其特征值隨kd的增大而向左移動,λ29所確定模態(tài)的特征值隨kd增大而向右移動。

    圖10 部分非振蕩模態(tài)特征值隨kd的變化Fig.10 Affection of eigenvalues of part of non-oscillation modes with kd change

    保持kd=1,kp以0.1的步長從1增加到10,系統(tǒng)7個振蕩模態(tài)的振蕩頻率和阻尼比變化情況如圖11所示。圖11中,λ1,2和λ5,6所確定振蕩模態(tài)的振蕩頻率和阻尼比隨kp增大而分別增大和減小。λ7,8所對應振蕩模態(tài)的振蕩頻率和阻尼比都隨kp增大而增大。λ13,14所確定振蕩模態(tài)的振蕩頻率和阻尼比隨kp增大而分別減小和增大。λ3,4所對應振蕩模態(tài)的振蕩頻率和阻尼比不隨kp變化而變化。λ3,4所確定振蕩模態(tài)的動態(tài)特性也不隨kp變化而變化,λ9—λ12所確定振蕩模態(tài)的動態(tài)特性也對kp變化不敏感。

    圖11 λ1-14振蕩頻率和阻尼比隨kp的變化Fig.11 Affection of oscillation frequency and damping ratio of λ1-14with kp change

    對于非振蕩模態(tài),選取僅與互聯(lián)系統(tǒng)調頻動態(tài)環(huán)節(jié)相關的6個模態(tài),研究kp變化時雙饋風機附加頻率控制對互聯(lián)系統(tǒng)調頻動態(tài)特性的影響。相關模態(tài)特征值的變化情況如圖12所示,圖中紅色和藍色叉號分別表示加入附加頻率控制前后的特征值。其中λ17、λ18、λ22、λ23、λ26所確定的模態(tài),其特征值隨kp的增大而向右移動,λ29所確定模態(tài)的特征值隨kp增大而向左移動。

    圖12 部分非振蕩模態(tài)特征值隨kp的變化Fig.12 Affection of eigenvalues of part of non-oscillation modes with kp change

    綜上分析可見,kd和kp的參數(shù)變化對λ1,2和λ5,6所確定的雙饋風機振蕩模態(tài)、對λ13,14所確定的系統(tǒng)調頻振蕩模態(tài)、對λ7,8所確定的雙饋風機和系統(tǒng)調頻動態(tài)相互耦合作用振蕩模態(tài),以及部分非振蕩模態(tài)的小干擾穩(wěn)定動態(tài)特性都會產生影響。

    3.3 時域仿真驗證

    分別在kd恒定、kp變化和kp恒定、kd變化2種情況下對圖3所示系統(tǒng)進行時域仿真,并將風機有功輸出與無附加頻率控制時的情況進行對比,仿真結果如圖13所示。在31個模態(tài)中,由于λ1,2所確定模態(tài)的耦合狀態(tài)變量包括定子電流,因此風機輸出有功的動態(tài)響應將反映該模態(tài)小干擾穩(wěn)定動態(tài)特性對控制參數(shù)變化響應。由3.1和3.2節(jié)的分析可知,加入附加頻率控制后,該振蕩模態(tài)的阻尼比會隨控制參數(shù)kd和kp增大而減小,圖13的仿真結果驗證了此動態(tài)特性的變化特點。

    圖13 風機輸出功率動態(tài)響應Fig.13 Dynamic response of the output of DFIG

    4 結 論

    本文主要討論了等值風電機組附加調頻控制對區(qū)域系統(tǒng)調頻動態(tài)特性的影響。在建立雙饋風機詳細狀態(tài)方程基礎上,研究了采用頻率支撐控制風機對互聯(lián)電力系統(tǒng)調頻動態(tài)特性的影響?;谔卣髦的B(tài)分析和仿真驗證,得到以下結論:

    1)含雙饋風機互聯(lián)電力系統(tǒng)的負荷-頻率控制包含分別由20個風機內部多物理控制環(huán)節(jié)和10個系統(tǒng)調頻動態(tài)環(huán)節(jié)耦合作用形成的多個模態(tài),其中包括7個振蕩模態(tài)和16個非振蕩模態(tài),并且這些模態(tài)都具有小干擾穩(wěn)定性。

    2)由慣性支撐和下垂控制組成的雙饋風機附加頻率支撐控制能有效提高系統(tǒng)的調頻動態(tài)特性,將系統(tǒng)受擾后的頻率最低點提高0.18 Hz。雖然頻率支撐控制對系統(tǒng)各模態(tài)小干擾穩(wěn)定動態(tài)特性的影響不明顯,但使雙饋風機多物理控制環(huán)節(jié)和系統(tǒng)調頻動態(tài)環(huán)節(jié)在特征值為-4.022±j7.184的振蕩模態(tài)和特征值為-40.183的非振蕩模態(tài)上發(fā)生耦合作用,并使系統(tǒng)調頻動態(tài)增加一個特征值為-94.104的非振蕩模態(tài)。

    3)在所有振蕩模態(tài)中,4個振蕩模態(tài)的振蕩頻率和阻尼比對頻率支撐控制的kd和kp參數(shù)變化具有較強敏感性,其余3個振蕩模態(tài)對參數(shù)變化敏感性較弱。同時,部分非振蕩模態(tài)的動態(tài)特性對頻率支撐控制的參數(shù)變化也具有敏感性。

    當研究多風場接入系統(tǒng)不同位置時,可先根據(jù)電氣聯(lián)系強弱對系統(tǒng)進行分區(qū),再分別研究不同區(qū)域中風機附加控制對系統(tǒng)調頻動態(tài)的影響。同時需要指出,模態(tài)分析結論僅針對本文提出的系統(tǒng)模型和運行工況,如何從理論角度將分析結論推廣到其他系統(tǒng)和不同運行工況,是后續(xù)研究值得深入開展的工作。

    猜你喜歡
    雙饋調頻特征值
    雙饋式可變速抽水蓄能機組運行控制
    大電機技術(2022年3期)2022-08-06 07:48:18
    基于Motor CAD的雙饋發(fā)電機溫升速算方法研究
    防爆電機(2021年6期)2022-01-17 02:40:18
    考慮頻率二次跌落抑制的風火聯(lián)合一次調頻控制
    能源工程(2021年5期)2021-11-20 05:50:42
    一類帶強制位勢的p-Laplace特征值問題
    單圈圖關聯(lián)矩陣的特征值
    調頻發(fā)射機技術改造
    調頻激勵器干擾的排除方法
    基于商奇異值分解的一類二次特征值反問題
    受平均斜率控制的Crowbar雙饋異步電機低電壓穿越
    調頻引信中噪聲調幅干擾的自適應抑制
    青草久久国产| 欧美午夜高清在线| videosex国产| 久久天躁狠狠躁夜夜2o2o| 女警被强在线播放| 欧美一级毛片孕妇| 亚洲 欧美一区二区三区| 国产亚洲精品一区二区www| 色婷婷久久久亚洲欧美| 黄色a级毛片大全视频| 日韩欧美在线二视频| 一进一出抽搐动态| 亚洲人成电影免费在线| 青草久久国产| 国产精品成人在线| 在线观看免费日韩欧美大片| 国产精品偷伦视频观看了| 两性夫妻黄色片| 欧美精品一区二区免费开放| 亚洲国产精品一区二区三区在线| 国产一区二区三区综合在线观看| 在线天堂中文资源库| 亚洲中文av在线| 极品教师在线免费播放| 五月开心婷婷网| 亚洲av日韩精品久久久久久密| 91成年电影在线观看| 国产一区二区在线av高清观看| 一二三四社区在线视频社区8| 精品一区二区三卡| 亚洲国产精品999在线| 久久欧美精品欧美久久欧美| √禁漫天堂资源中文www| 老司机福利观看| 成人三级做爰电影| 欧美日韩亚洲国产一区二区在线观看| 两人在一起打扑克的视频| 亚洲avbb在线观看| √禁漫天堂资源中文www| 久久热在线av| 久久精品国产综合久久久| 热re99久久国产66热| 怎么达到女性高潮| 久久狼人影院| 男女下面进入的视频免费午夜 | 亚洲一区二区三区欧美精品| 成熟少妇高潮喷水视频| 在线十欧美十亚洲十日本专区| 亚洲三区欧美一区| 国产欧美日韩一区二区三| 欧美在线一区亚洲| 欧美日韩亚洲国产一区二区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 女人高潮潮喷娇喘18禁视频| 88av欧美| 香蕉国产在线看| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩另类电影网站| 午夜免费观看网址| 99riav亚洲国产免费| 欧美中文综合在线视频| 欧美日韩精品网址| 一区福利在线观看| 一个人观看的视频www高清免费观看 | 亚洲精品美女久久久久99蜜臀| 国产黄色免费在线视频| www.熟女人妻精品国产| 侵犯人妻中文字幕一二三四区| 日韩欧美一区视频在线观看| 日本a在线网址| 首页视频小说图片口味搜索| 成人永久免费在线观看视频| 亚洲欧美日韩无卡精品| 国产一区二区三区在线臀色熟女 | 色在线成人网| 日本a在线网址| 精品久久蜜臀av无| 欧美激情久久久久久爽电影 | 69av精品久久久久久| 久久久久国产一级毛片高清牌| 精品卡一卡二卡四卡免费| 亚洲全国av大片| 久久香蕉国产精品| 黄色片一级片一级黄色片| 制服人妻中文乱码| 国产无遮挡羞羞视频在线观看| 欧美 亚洲 国产 日韩一| 欧美av亚洲av综合av国产av| 美女高潮到喷水免费观看| 后天国语完整版免费观看| 高潮久久久久久久久久久不卡| 12—13女人毛片做爰片一| 99在线人妻在线中文字幕| 欧美日韩一级在线毛片| 丁香六月欧美| 乱人伦中国视频| 午夜福利在线免费观看网站| 女人精品久久久久毛片| 不卡av一区二区三区| 免费少妇av软件| 在线看a的网站| 校园春色视频在线观看| 欧美性长视频在线观看| 欧美av亚洲av综合av国产av| 在线看a的网站| 高清av免费在线| 亚洲av成人一区二区三| 丁香欧美五月| 99国产精品一区二区蜜桃av| 国内久久婷婷六月综合欲色啪| 在线观看66精品国产| av视频免费观看在线观看| 电影成人av| 成熟少妇高潮喷水视频| 国产视频一区二区在线看| 久久精品国产亚洲av香蕉五月| cao死你这个sao货| 亚洲国产中文字幕在线视频| av免费在线观看网站| 国产精品国产av在线观看| 欧美在线一区亚洲| 级片在线观看| 精品一区二区三区视频在线观看免费 | 黄频高清免费视频| 色综合欧美亚洲国产小说| 99久久99久久久精品蜜桃| 欧美日韩黄片免| 亚洲国产看品久久| 亚洲熟妇中文字幕五十中出 | 国产成人免费无遮挡视频| 手机成人av网站| 在线观看免费视频日本深夜| 亚洲国产中文字幕在线视频| 亚洲在线自拍视频| xxx96com| 热re99久久国产66热| 一二三四在线观看免费中文在| 91成人精品电影| 成人亚洲精品一区在线观看| 在线观看免费日韩欧美大片| 国产精品亚洲一级av第二区| 欧美+亚洲+日韩+国产| 美女高潮到喷水免费观看| 18禁裸乳无遮挡免费网站照片 | 日韩免费av在线播放| 国产精品99久久99久久久不卡| 欧美成人午夜精品| 亚洲av五月六月丁香网| 国产亚洲av高清不卡| 亚洲国产欧美网| 中文字幕精品免费在线观看视频| 亚洲欧美日韩另类电影网站| 黄片播放在线免费| 久久久精品国产亚洲av高清涩受| 一边摸一边抽搐一进一小说| 大码成人一级视频| 午夜免费激情av| 亚洲久久久国产精品| 丰满的人妻完整版| 天天影视国产精品| 亚洲欧美日韩无卡精品| 黄色成人免费大全| 欧美 亚洲 国产 日韩一| 99精品在免费线老司机午夜| 天堂影院成人在线观看| 国产精品永久免费网站| 人妻久久中文字幕网| www.精华液| x7x7x7水蜜桃| 亚洲免费av在线视频| 亚洲成a人片在线一区二区| 男人的好看免费观看在线视频 | 国产亚洲欧美在线一区二区| 美女高潮到喷水免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 大陆偷拍与自拍| 久久99一区二区三区| 十八禁人妻一区二区| 99香蕉大伊视频| 国产区一区二久久| svipshipincom国产片| 欧美久久黑人一区二区| 国产99久久九九免费精品| 不卡av一区二区三区| 亚洲精品在线美女| 国产男靠女视频免费网站| 国产亚洲精品一区二区www| 亚洲精品一卡2卡三卡4卡5卡| 老司机靠b影院| 国产xxxxx性猛交| 美女午夜性视频免费| 国产高清激情床上av| 黄片小视频在线播放| 午夜免费观看网址| a级片在线免费高清观看视频| 美女午夜性视频免费| 欧美一区二区精品小视频在线| 国产激情欧美一区二区| 日本三级黄在线观看| 很黄的视频免费| 久久久国产精品麻豆| 国产成人欧美在线观看| 99久久国产精品久久久| 天堂动漫精品| 啦啦啦免费观看视频1| 国产欧美日韩一区二区三区在线| 夫妻午夜视频| 午夜免费激情av| 欧美成人性av电影在线观看| 日韩三级视频一区二区三区| 国产一卡二卡三卡精品| 男人舔女人下体高潮全视频| 色老头精品视频在线观看| 免费久久久久久久精品成人欧美视频| 一夜夜www| 国产一区二区激情短视频| 亚洲av美国av| 国产男靠女视频免费网站| 午夜视频精品福利| 夜夜看夜夜爽夜夜摸 | 1024视频免费在线观看| 一进一出抽搐动态| 精品久久久久久电影网| 狂野欧美激情性xxxx| √禁漫天堂资源中文www| 两个人看的免费小视频| 国产片内射在线| svipshipincom国产片| 午夜免费鲁丝| 亚洲中文日韩欧美视频| 久久人人爽av亚洲精品天堂| 一边摸一边抽搐一进一小说| 极品教师在线免费播放| 丝袜美腿诱惑在线| 俄罗斯特黄特色一大片| 亚洲av熟女| 啪啪无遮挡十八禁网站| 波多野结衣一区麻豆| 久久久国产精品麻豆| 精品久久久久久久毛片微露脸| 最近最新免费中文字幕在线| 久久精品亚洲av国产电影网| 久久这里只有精品19| 国产单亲对白刺激| 精品午夜福利视频在线观看一区| 伦理电影免费视频| 久久精品亚洲熟妇少妇任你| 一边摸一边抽搐一进一出视频| 在线十欧美十亚洲十日本专区| 一级黄色大片毛片| 91大片在线观看| 国产1区2区3区精品| 妹子高潮喷水视频| 波多野结衣高清无吗| 无遮挡黄片免费观看| 久久久国产欧美日韩av| 999久久久精品免费观看国产| 男人舔女人的私密视频| 一边摸一边做爽爽视频免费| 午夜精品国产一区二区电影| 男人舔女人下体高潮全视频| 黄色女人牲交| 国产一区二区激情短视频| 在线看a的网站| 亚洲一区二区三区欧美精品| 母亲3免费完整高清在线观看| 麻豆国产av国片精品| 久久久国产成人免费| 大码成人一级视频| 两个人看的免费小视频| 亚洲自偷自拍图片 自拍| 中文亚洲av片在线观看爽| 国产欧美日韩一区二区三区在线| 亚洲成人精品中文字幕电影 | 国产精品野战在线观看 | 50天的宝宝边吃奶边哭怎么回事| 精品国内亚洲2022精品成人| 国产伦一二天堂av在线观看| 国产亚洲av高清不卡| 国产深夜福利视频在线观看| 波多野结衣av一区二区av| 亚洲成av片中文字幕在线观看| 亚洲欧美日韩无卡精品| 老汉色av国产亚洲站长工具| 国产成年人精品一区二区 | 两性夫妻黄色片| 色婷婷av一区二区三区视频| 大香蕉久久成人网| 亚洲欧美激情综合另类| 黄色丝袜av网址大全| 黄色a级毛片大全视频| 一边摸一边抽搐一进一出视频| 免费在线观看影片大全网站| 日韩免费高清中文字幕av| 亚洲欧美一区二区三区久久| 亚洲精品国产区一区二| 国产深夜福利视频在线观看| 午夜a级毛片| 日本 av在线| 19禁男女啪啪无遮挡网站| 精品久久久久久久毛片微露脸| 久久99一区二区三区| 黄色片一级片一级黄色片| 伦理电影免费视频| 最近最新中文字幕大全电影3 | 欧美在线一区亚洲| 欧美成狂野欧美在线观看| 欧美黄色淫秽网站| e午夜精品久久久久久久| 91麻豆av在线| www.精华液| 亚洲在线自拍视频| 国内久久婷婷六月综合欲色啪| 午夜免费鲁丝| 日韩精品中文字幕看吧| 不卡一级毛片| av天堂在线播放| 宅男免费午夜| 在线天堂中文资源库| 狠狠狠狠99中文字幕| 最新在线观看一区二区三区| 国产日韩一区二区三区精品不卡| 欧美性长视频在线观看| 午夜精品国产一区二区电影| 精品一区二区三卡| 亚洲国产精品999在线| 黄色视频不卡| 色精品久久人妻99蜜桃| 狂野欧美激情性xxxx| 成人精品一区二区免费| 在线观看www视频免费| 精品一品国产午夜福利视频| 日本黄色日本黄色录像| 亚洲av熟女| 女警被强在线播放| 国产熟女午夜一区二区三区| 日本a在线网址| 大型av网站在线播放| 如日韩欧美国产精品一区二区三区| 99久久精品国产亚洲精品| 午夜久久久在线观看| 亚洲精品av麻豆狂野| 国产免费男女视频| 精品电影一区二区在线| 变态另类成人亚洲欧美熟女 | 亚洲成人精品中文字幕电影 | 激情视频va一区二区三区| 日韩欧美国产一区二区入口| 成人三级黄色视频| 欧美在线黄色| 亚洲激情在线av| 狠狠狠狠99中文字幕| 亚洲一码二码三码区别大吗| 大香蕉久久成人网| 久久久久国产一级毛片高清牌| 亚洲av片天天在线观看| 成年版毛片免费区| 最好的美女福利视频网| 女性生殖器流出的白浆| 国产精品国产av在线观看| 国产伦人伦偷精品视频| 男人的好看免费观看在线视频 | 一区二区日韩欧美中文字幕| 国产亚洲欧美在线一区二区| 亚洲五月婷婷丁香| 多毛熟女@视频| 国产精品一区二区精品视频观看| 国内久久婷婷六月综合欲色啪| 高清av免费在线| 亚洲国产欧美网| 水蜜桃什么品种好| 国产亚洲精品久久久久久毛片| 精品久久久久久,| 亚洲熟妇熟女久久| 成人国语在线视频| 桃色一区二区三区在线观看| 黄片播放在线免费| 夜夜看夜夜爽夜夜摸 | cao死你这个sao货| 日本一区二区免费在线视频| 9色porny在线观看| 国产乱人伦免费视频| 亚洲国产精品sss在线观看 | 久久香蕉国产精品| 99热只有精品国产| 国产高清激情床上av| 交换朋友夫妻互换小说| 真人一进一出gif抽搐免费| 国产精品一区二区免费欧美| 国产亚洲精品一区二区www| 一区福利在线观看| 久久久久久久精品吃奶| 久久精品成人免费网站| 国产精品av久久久久免费| 99riav亚洲国产免费| 久久久久久久精品吃奶| 久久这里只有精品19| 少妇 在线观看| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区四区五区乱码| 中文字幕高清在线视频| av天堂在线播放| 久久中文字幕一级| 高清av免费在线| 国产激情欧美一区二区| 人妻丰满熟妇av一区二区三区| 777久久人妻少妇嫩草av网站| 成人av一区二区三区在线看| 亚洲精品国产精品久久久不卡| 久久久国产成人免费| 人妻久久中文字幕网| 久久国产亚洲av麻豆专区| 老司机靠b影院| 免费高清视频大片| videosex国产| 欧美日韩视频精品一区| 成人18禁高潮啪啪吃奶动态图| 国产亚洲av高清不卡| 久久精品国产亚洲av香蕉五月| 国产精品久久视频播放| 亚洲欧洲精品一区二区精品久久久| 在线播放国产精品三级| 国产精品国产高清国产av| 久久亚洲精品不卡| 亚洲欧美一区二区三区黑人| 日本撒尿小便嘘嘘汇集6| 欧美中文日本在线观看视频| 激情在线观看视频在线高清| 在线观看一区二区三区| 性色av乱码一区二区三区2| 日韩欧美在线二视频| 亚洲三区欧美一区| a级片在线免费高清观看视频| 国产精华一区二区三区| 日本wwww免费看| 国产精品永久免费网站| 久久精品亚洲熟妇少妇任你| av在线播放免费不卡| 叶爱在线成人免费视频播放| 97超级碰碰碰精品色视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 在线十欧美十亚洲十日本专区| 国产xxxxx性猛交| 一个人观看的视频www高清免费观看 | 亚洲国产精品合色在线| 亚洲自拍偷在线| 叶爱在线成人免费视频播放| 久久久久国产一级毛片高清牌| 热99国产精品久久久久久7| 91九色精品人成在线观看| 精品福利观看| 一区二区日韩欧美中文字幕| 最近最新中文字幕大全免费视频| 亚洲一区二区三区色噜噜 | 国产亚洲精品一区二区www| 午夜视频精品福利| 村上凉子中文字幕在线| 亚洲精品国产精品久久久不卡| 两性夫妻黄色片| 欧美大码av| 久久国产精品影院| 宅男免费午夜| 两性夫妻黄色片| 97超级碰碰碰精品色视频在线观看| 成人精品一区二区免费| 国产高清视频在线播放一区| 成人免费观看视频高清| 夜夜爽天天搞| 精品久久久久久成人av| 一级黄色大片毛片| 波多野结衣高清无吗| 搡老岳熟女国产| 国产精品自产拍在线观看55亚洲| 成年女人毛片免费观看观看9| 久久精品亚洲熟妇少妇任你| 精品国产一区二区三区四区第35| 免费人成视频x8x8入口观看| 久久久国产一区二区| 激情在线观看视频在线高清| 亚洲精品一二三| 嫁个100分男人电影在线观看| 欧美最黄视频在线播放免费 | 欧美黄色淫秽网站| 麻豆一二三区av精品| 亚洲国产精品999在线| 狠狠狠狠99中文字幕| √禁漫天堂资源中文www| 国产亚洲av高清不卡| 欧美在线黄色| 国产精品久久久久久人妻精品电影| 国产激情久久老熟女| 午夜激情av网站| 88av欧美| 精品熟女少妇八av免费久了| 韩国av一区二区三区四区| 国产三级在线视频| 免费观看人在逋| 精品久久久久久成人av| 中文字幕av电影在线播放| 九色亚洲精品在线播放| 亚洲激情在线av| 久99久视频精品免费| 女人爽到高潮嗷嗷叫在线视频| 又紧又爽又黄一区二区| 亚洲成av片中文字幕在线观看| aaaaa片日本免费| 好看av亚洲va欧美ⅴa在| 久热这里只有精品99| 欧美乱妇无乱码| 亚洲欧美一区二区三区久久| 在线观看一区二区三区激情| 最近最新免费中文字幕在线| 午夜免费鲁丝| 亚洲精品在线观看二区| 欧美在线一区亚洲| 91大片在线观看| 亚洲伊人色综图| 久久人妻熟女aⅴ| 极品人妻少妇av视频| 国产精品乱码一区二三区的特点 | 伦理电影免费视频| 999久久久国产精品视频| 性色av乱码一区二区三区2| 在线观看免费午夜福利视频| 精品国产国语对白av| 男人舔女人的私密视频| 欧美日韩亚洲综合一区二区三区_| 欧美日韩亚洲国产一区二区在线观看| 一级a爱视频在线免费观看| 亚洲精品美女久久久久99蜜臀| 亚洲一区高清亚洲精品| 欧美久久黑人一区二区| 人人妻,人人澡人人爽秒播| 亚洲国产看品久久| 久久天躁狠狠躁夜夜2o2o| 国产乱人伦免费视频| 久久精品国产亚洲av高清一级| 亚洲专区字幕在线| 日韩视频一区二区在线观看| 一区二区三区国产精品乱码| av有码第一页| 精品卡一卡二卡四卡免费| 嫩草影院精品99| 精品卡一卡二卡四卡免费| 午夜福利免费观看在线| 国产成人免费无遮挡视频| 亚洲精品美女久久av网站| 日韩大尺度精品在线看网址 | 91字幕亚洲| 男人操女人黄网站| 午夜精品国产一区二区电影| 午夜影院日韩av| 热99re8久久精品国产| 亚洲av五月六月丁香网| 性色av乱码一区二区三区2| 超碰成人久久| 男女下面进入的视频免费午夜 | 国产片内射在线| 在线观看免费视频日本深夜| 19禁男女啪啪无遮挡网站| 99riav亚洲国产免费| 看片在线看免费视频| 男人舔女人下体高潮全视频| 久久午夜亚洲精品久久| 女性被躁到高潮视频| 久久香蕉国产精品| 熟女少妇亚洲综合色aaa.| 亚洲av成人不卡在线观看播放网| 大型av网站在线播放| 熟女少妇亚洲综合色aaa.| 欧美日本亚洲视频在线播放| 女人精品久久久久毛片| 男女午夜视频在线观看| 操美女的视频在线观看| 国产精品影院久久| 在线看a的网站| 色婷婷久久久亚洲欧美| 欧美成人性av电影在线观看| 国产99久久九九免费精品| 美女午夜性视频免费| 国产精品 国内视频| 91在线观看av| 亚洲九九香蕉| 日韩大尺度精品在线看网址 | 国产av一区二区精品久久| 亚洲欧美激情综合另类| 精品国产超薄肉色丝袜足j| 满18在线观看网站| √禁漫天堂资源中文www| av中文乱码字幕在线| 国产熟女午夜一区二区三区| 精品国产乱子伦一区二区三区| 十八禁网站免费在线| 国产亚洲欧美98| 老司机亚洲免费影院| 欧美日本中文国产一区发布| av片东京热男人的天堂| 成人18禁在线播放| 老司机靠b影院| 日韩大码丰满熟妇| 亚洲国产精品合色在线| 成人免费观看视频高清| 正在播放国产对白刺激| 老熟妇仑乱视频hdxx| 中文字幕人妻丝袜制服| 国产激情欧美一区二区| 国产一卡二卡三卡精品| 久久久久久亚洲精品国产蜜桃av| 中文字幕av电影在线播放| 午夜福利一区二区在线看| 国产精品亚洲av一区麻豆| 老司机靠b影院| 久久久久九九精品影院| 热re99久久精品国产66热6| 国产黄色免费在线视频|