• <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雙饋異步電機低電壓穿越
    調頻引信中噪聲調幅干擾的自適應抑制
    国产亚洲精品久久久久久毛片| 精品久久久久久电影网| 国产精品野战在线观看 | 视频区欧美日本亚洲| 满18在线观看网站| 国产精品亚洲av一区麻豆| 成人精品一区二区免费| 神马国产精品三级电影在线观看 | 99国产精品99久久久久| 老汉色∧v一级毛片| 在线永久观看黄色视频| 一边摸一边抽搐一进一小说| 国产aⅴ精品一区二区三区波| 久久久久国产精品人妻aⅴ院| 大型黄色视频在线免费观看| 亚洲人成77777在线视频| 国产熟女午夜一区二区三区| 一级毛片女人18水好多| 欧美一区二区精品小视频在线| 欧美黄色片欧美黄色片| 在线观看免费日韩欧美大片| 亚洲在线自拍视频| 成人av一区二区三区在线看| 久久久久久久午夜电影 | 亚洲成人精品中文字幕电影 | 精品国产乱子伦一区二区三区| 久热爱精品视频在线9| 日韩 欧美 亚洲 中文字幕| 少妇裸体淫交视频免费看高清 | 女人被躁到高潮嗷嗷叫费观| 久久精品人人爽人人爽视色| 成人国语在线视频| 看黄色毛片网站| 精品电影一区二区在线| 免费高清在线观看日韩| 91麻豆av在线| 超碰成人久久| 国产高清激情床上av| 亚洲av熟女| 十分钟在线观看高清视频www| 精品卡一卡二卡四卡免费| 男人舔女人的私密视频| 大陆偷拍与自拍| 精品一区二区三区四区五区乱码| 国产成+人综合+亚洲专区| 久久欧美精品欧美久久欧美| 免费高清在线观看日韩| 国产单亲对白刺激| 国产熟女xx| 久热爱精品视频在线9| 国产精品成人在线| 久久久国产精品麻豆| 亚洲激情在线av| 中文字幕人妻丝袜一区二区| 大香蕉久久成人网| 国产99久久九九免费精品| 国产一区二区三区视频了| 欧美黄色片欧美黄色片| 99久久国产精品久久久| 亚洲一区二区三区不卡视频| 国产野战对白在线观看| 高清黄色对白视频在线免费看| 可以免费在线观看a视频的电影网站| 国产免费男女视频| 国产精品香港三级国产av潘金莲| 久久亚洲真实| 久久热在线av| 18禁观看日本| 国产伦人伦偷精品视频| √禁漫天堂资源中文www| 97人妻天天添夜夜摸| 男女做爰动态图高潮gif福利片 | 18禁裸乳无遮挡免费网站照片 | 后天国语完整版免费观看| 91老司机精品| 久久欧美精品欧美久久欧美| 中文字幕另类日韩欧美亚洲嫩草| 欧美中文综合在线视频| 亚洲精品在线美女| 亚洲中文av在线| 一进一出抽搐动态| 午夜福利影视在线免费观看| 丁香欧美五月| 超色免费av| 精品第一国产精品| 男人舔女人的私密视频| 久久中文字幕人妻熟女| 欧洲精品卡2卡3卡4卡5卡区| 一级黄色大片毛片| 国产91精品成人一区二区三区| 十八禁人妻一区二区| 久久国产精品人妻蜜桃| 成人国产一区最新在线观看| 久久久精品欧美日韩精品| 99久久精品国产亚洲精品| 国产伦一二天堂av在线观看| 他把我摸到了高潮在线观看| 日韩有码中文字幕| 久久精品亚洲熟妇少妇任你| 久久中文看片网| 99在线视频只有这里精品首页| 可以在线观看毛片的网站| 天天躁夜夜躁狠狠躁躁| 国产伦人伦偷精品视频| 国产日韩一区二区三区精品不卡| 999久久久精品免费观看国产| 老汉色av国产亚洲站长工具| 欧美丝袜亚洲另类 | 免费少妇av软件| 首页视频小说图片口味搜索| 黄色毛片三级朝国网站| 美女大奶头视频| 香蕉丝袜av| 国产亚洲欧美98| 99精品在免费线老司机午夜| 久久久久亚洲av毛片大全| 午夜福利一区二区在线看| 亚洲av第一区精品v没综合| 欧美乱妇无乱码| 搡老乐熟女国产| 水蜜桃什么品种好| 伊人久久大香线蕉亚洲五| 两个人免费观看高清视频| av网站免费在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 精品国内亚洲2022精品成人| 看片在线看免费视频| 免费日韩欧美在线观看| 色哟哟哟哟哟哟| 亚洲国产精品sss在线观看 | 日韩av在线大香蕉| 无人区码免费观看不卡| 国产片内射在线| 精品高清国产在线一区| 欧美精品亚洲一区二区| 神马国产精品三级电影在线观看 | 电影成人av| 热99re8久久精品国产| 天堂影院成人在线观看| 国产1区2区3区精品| 亚洲午夜精品一区,二区,三区| 12—13女人毛片做爰片一| 99re在线观看精品视频| 琪琪午夜伦伦电影理论片6080| 性少妇av在线| 一边摸一边做爽爽视频免费| 国产精品自产拍在线观看55亚洲| 亚洲精品国产色婷婷电影| 国产成人欧美| а√天堂www在线а√下载| 中文字幕最新亚洲高清| 国产野战对白在线观看| 久久精品亚洲av国产电影网| 欧美 亚洲 国产 日韩一| 一区二区日韩欧美中文字幕| 国产精品美女特级片免费视频播放器 | 国产精品永久免费网站| 国产深夜福利视频在线观看| 精品福利永久在线观看| 黄色怎么调成土黄色| 亚洲精品在线观看二区| 交换朋友夫妻互换小说| 91精品国产国语对白视频| 美女大奶头视频| 午夜免费激情av| 日日摸夜夜添夜夜添小说| 国产野战对白在线观看| 午夜精品在线福利| 久久精品亚洲熟妇少妇任你| 19禁男女啪啪无遮挡网站| 久久久国产成人精品二区 | 99在线视频只有这里精品首页| 国产熟女xx| 少妇被粗大的猛进出69影院| 日本vs欧美在线观看视频| 手机成人av网站| 99热只有精品国产| 精品午夜福利视频在线观看一区| 夜夜看夜夜爽夜夜摸 | 夜夜躁狠狠躁天天躁| 亚洲av成人不卡在线观看播放网| 1024视频免费在线观看| 欧美乱码精品一区二区三区| 午夜福利,免费看| 久99久视频精品免费| 久久人人爽av亚洲精品天堂| 免费在线观看黄色视频的| 久久精品人人爽人人爽视色| 中文欧美无线码| 亚洲免费av在线视频| 久久国产精品人妻蜜桃| www日本在线高清视频| 国产精品一区二区精品视频观看| 亚洲国产精品999在线| 国产亚洲精品久久久久5区| 操美女的视频在线观看| 亚洲av片天天在线观看| 国产精品影院久久| 久久久久久大精品| 热99re8久久精品国产| 日日摸夜夜添夜夜添小说| 亚洲专区中文字幕在线| 久久精品影院6| 国产av又大| 91老司机精品| 免费在线观看影片大全网站| 天堂√8在线中文| 久久久久国产一级毛片高清牌| 成人特级黄色片久久久久久久| 国内毛片毛片毛片毛片毛片| 免费一级毛片在线播放高清视频 | 久久精品影院6| 精品日产1卡2卡| 老司机靠b影院| 性色av乱码一区二区三区2| 精品久久久久久久久久免费视频 | 精品熟女少妇八av免费久了| 777久久人妻少妇嫩草av网站| 亚洲av日韩精品久久久久久密| 免费久久久久久久精品成人欧美视频| 国产三级在线视频| 别揉我奶头~嗯~啊~动态视频| 超碰成人久久| 啦啦啦 在线观看视频| 免费久久久久久久精品成人欧美视频| 亚洲美女黄片视频| 国产国语露脸激情在线看| 亚洲av五月六月丁香网| 亚洲性夜色夜夜综合| 欧美大码av| 国产男靠女视频免费网站| 一区二区三区精品91| av电影中文网址| 天堂动漫精品| 久久久久精品国产欧美久久久| 黑人巨大精品欧美一区二区蜜桃| 黄片大片在线免费观看| 好看av亚洲va欧美ⅴa在| 深夜精品福利| 久久伊人香网站| 亚洲国产欧美日韩在线播放| 视频在线观看一区二区三区| 日韩国内少妇激情av| 免费高清视频大片| 国产人伦9x9x在线观看| 国产三级在线视频| 三级毛片av免费| 在线看a的网站| 我的亚洲天堂| 久久精品国产清高在天天线| 桃色一区二区三区在线观看| 国产深夜福利视频在线观看| 99香蕉大伊视频| 少妇 在线观看| 国产亚洲精品久久久久5区| 人人妻人人添人人爽欧美一区卜| 婷婷六月久久综合丁香| 精品欧美一区二区三区在线| 亚洲 国产 在线| 日韩大码丰满熟妇| 777久久人妻少妇嫩草av网站| 国产一区二区三区在线臀色熟女 | 久久久精品欧美日韩精品| 最近最新中文字幕大全电影3 | 中文字幕精品免费在线观看视频| 久久久国产精品麻豆| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品999在线| 99香蕉大伊视频| 黑人巨大精品欧美一区二区mp4| 国产亚洲精品久久久久5区| 可以在线观看毛片的网站| 久久天堂一区二区三区四区| 精品日产1卡2卡| 婷婷精品国产亚洲av在线| 久久午夜亚洲精品久久| 中文亚洲av片在线观看爽| 高清黄色对白视频在线免费看| 久久久久久久久中文| 男女床上黄色一级片免费看| 级片在线观看| 丁香欧美五月| 日本vs欧美在线观看视频| 老司机靠b影院| 午夜福利在线观看吧| 国产成+人综合+亚洲专区| av有码第一页| 日韩大尺度精品在线看网址 | 久久精品亚洲av国产电影网| 国产亚洲av高清不卡| 狠狠狠狠99中文字幕| 99国产精品99久久久久| 女性生殖器流出的白浆| 亚洲精品久久午夜乱码| 国产蜜桃级精品一区二区三区| 性欧美人与动物交配| 日韩免费av在线播放| 男人的好看免费观看在线视频 | 欧美人与性动交α欧美精品济南到| 九色亚洲精品在线播放| 99国产精品一区二区蜜桃av| 老司机深夜福利视频在线观看| 久久久久久大精品| 看黄色毛片网站| 亚洲精品美女久久久久99蜜臀| 免费av中文字幕在线| 巨乳人妻的诱惑在线观看| 国产精品秋霞免费鲁丝片| 亚洲精品中文字幕在线视频| 国产av又大| 国产亚洲精品久久久久久毛片| 免费搜索国产男女视频| 亚洲av五月六月丁香网| 交换朋友夫妻互换小说| 欧美日韩精品网址| 男女下面插进去视频免费观看| 欧美日韩国产mv在线观看视频| 99在线人妻在线中文字幕| 最近最新中文字幕大全电影3 | 亚洲成人免费电影在线观看| 桃红色精品国产亚洲av| 久久精品91蜜桃| 最近最新免费中文字幕在线| 一级黄色大片毛片| 电影成人av| 欧美激情久久久久久爽电影 | 欧美一级毛片孕妇| 亚洲国产精品999在线| 夜夜夜夜夜久久久久| 好男人电影高清在线观看| 成人黄色视频免费在线看| 夜夜看夜夜爽夜夜摸 | www日本在线高清视频| 久久人妻福利社区极品人妻图片| 亚洲国产欧美网| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品亚洲av一区麻豆| 久久久久亚洲av毛片大全| 18禁黄网站禁片午夜丰满| 超碰97精品在线观看| 日韩高清综合在线| 久久香蕉精品热| 又大又爽又粗| 国产精品一区二区精品视频观看| 麻豆av在线久日| 又黄又爽又免费观看的视频| 最近最新中文字幕大全免费视频| 高清欧美精品videossex| 亚洲性夜色夜夜综合| 99久久99久久久精品蜜桃| 美女高潮到喷水免费观看| 黄色怎么调成土黄色| 精品国产美女av久久久久小说| 国产精品 欧美亚洲| 在线观看一区二区三区激情| 精品欧美一区二区三区在线| 精品一品国产午夜福利视频| 国产片内射在线| 另类亚洲欧美激情| 制服人妻中文乱码| 人人妻人人爽人人添夜夜欢视频| 免费在线观看日本一区| 美女高潮到喷水免费观看| 国产成人免费无遮挡视频| 久久久水蜜桃国产精品网| 人人妻人人添人人爽欧美一区卜| 欧美激情极品国产一区二区三区| 看黄色毛片网站| 午夜福利免费观看在线| 国产伦一二天堂av在线观看| 亚洲性夜色夜夜综合| 国产熟女xx| 久久人妻福利社区极品人妻图片| 女性被躁到高潮视频| 久久精品国产亚洲av高清一级| 久久精品亚洲熟妇少妇任你| 免费搜索国产男女视频| 在线观看免费视频日本深夜| 久久精品亚洲av国产电影网| 欧美日韩亚洲高清精品| 老司机午夜十八禁免费视频| 男女午夜视频在线观看| 99久久综合精品五月天人人| 日韩精品免费视频一区二区三区| 成人影院久久| 人人妻人人爽人人添夜夜欢视频| 啦啦啦在线免费观看视频4| 长腿黑丝高跟| 色综合婷婷激情| 国产精品自产拍在线观看55亚洲| 欧美性长视频在线观看| 亚洲精华国产精华精| 9热在线视频观看99| 亚洲aⅴ乱码一区二区在线播放 | 黄色片一级片一级黄色片| 丝袜在线中文字幕| 女警被强在线播放| 国产亚洲精品第一综合不卡| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久久人人做人人爽| 91大片在线观看| 日韩欧美三级三区| 久久 成人 亚洲| 亚洲av日韩精品久久久久久密| 黑人操中国人逼视频| 天天添夜夜摸| 亚洲国产精品合色在线| 日韩高清综合在线| 欧美日韩亚洲高清精品| 97超级碰碰碰精品色视频在线观看| 日本三级黄在线观看| 男人舔女人下体高潮全视频| 国产精品永久免费网站| 91精品三级在线观看| 免费观看精品视频网站| 午夜福利,免费看| 欧美最黄视频在线播放免费 | 人妻丰满熟妇av一区二区三区| 亚洲熟妇中文字幕五十中出 | 久久久久久久精品吃奶| 午夜两性在线视频| 99热国产这里只有精品6| 人人妻,人人澡人人爽秒播| 欧美黄色淫秽网站| av在线播放免费不卡| 成人三级做爰电影| 日韩欧美一区二区三区在线观看| 黄色a级毛片大全视频| 亚洲,欧美精品.| av国产精品久久久久影院| 美女高潮到喷水免费观看| 韩国av一区二区三区四区| 人人澡人人妻人| 男女床上黄色一级片免费看| 热re99久久国产66热| 中文字幕人妻丝袜制服| 水蜜桃什么品种好| 成人黄色视频免费在线看| 国产深夜福利视频在线观看| 欧美成人免费av一区二区三区| 80岁老熟妇乱子伦牲交| 我的亚洲天堂| 99精国产麻豆久久婷婷| 精品福利观看| 后天国语完整版免费观看| www.自偷自拍.com| 老司机午夜福利在线观看视频| 日韩欧美一区视频在线观看| 怎么达到女性高潮| 亚洲人成77777在线视频| 免费日韩欧美在线观看| 黄片大片在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 成人亚洲精品av一区二区 | 看黄色毛片网站| 一夜夜www| 777久久人妻少妇嫩草av网站| 国产精华一区二区三区| 人人妻,人人澡人人爽秒播| 色综合欧美亚洲国产小说| 精品国产美女av久久久久小说| 免费在线观看日本一区| 免费观看人在逋| 一二三四在线观看免费中文在| 一区二区三区精品91| 伊人久久大香线蕉亚洲五| 大码成人一级视频| 免费在线观看完整版高清| 免费一级毛片在线播放高清视频 | 成熟少妇高潮喷水视频| 欧美黑人欧美精品刺激| 91九色精品人成在线观看| 在线看a的网站| 亚洲五月色婷婷综合| 欧美日韩福利视频一区二区| 村上凉子中文字幕在线| 国产人伦9x9x在线观看| 国产成年人精品一区二区 | 成人免费观看视频高清| 亚洲国产精品合色在线| 欧美成狂野欧美在线观看| 欧美+亚洲+日韩+国产| av片东京热男人的天堂| 久久亚洲精品不卡| 精品久久久久久久久久免费视频 | 日本vs欧美在线观看视频| 男女下面进入的视频免费午夜 | 后天国语完整版免费观看| 亚洲少妇的诱惑av| 国产av一区在线观看免费| 久久精品影院6| 久久精品亚洲精品国产色婷小说| 十八禁人妻一区二区| 韩国精品一区二区三区| 国产精品 欧美亚洲| 国产午夜精品久久久久久| 不卡av一区二区三区| 黄片大片在线免费观看| 涩涩av久久男人的天堂| 中文字幕最新亚洲高清| 成在线人永久免费视频| 亚洲五月色婷婷综合| 国产国语露脸激情在线看| 一a级毛片在线观看| 黄色怎么调成土黄色| 免费观看人在逋| 午夜两性在线视频| 亚洲成a人片在线一区二区| 国产精品久久久人人做人人爽| 国产1区2区3区精品| 日韩av在线大香蕉| 纯流量卡能插随身wifi吗| 婷婷丁香在线五月| 老汉色av国产亚洲站长工具| 免费日韩欧美在线观看| 麻豆国产av国片精品| 丁香欧美五月| xxxhd国产人妻xxx| 欧美成人性av电影在线观看| 悠悠久久av| 老司机福利观看| 91精品三级在线观看| 欧美人与性动交α欧美软件| 黄片大片在线免费观看| 亚洲成人免费av在线播放| 久久影院123| 久久中文字幕人妻熟女| 欧美成人午夜精品| 久久久久久亚洲精品国产蜜桃av| 18禁国产床啪视频网站| 午夜福利影视在线免费观看| 99riav亚洲国产免费| 新久久久久国产一级毛片| 久久午夜综合久久蜜桃| 成在线人永久免费视频| 日本黄色视频三级网站网址| 久久精品国产99精品国产亚洲性色 | 国产男靠女视频免费网站| 天堂影院成人在线观看| 久久中文字幕一级| 亚洲精品在线美女| e午夜精品久久久久久久| 91麻豆精品激情在线观看国产 | 视频在线观看一区二区三区| 国产av在哪里看| 亚洲精品国产色婷婷电影| 国产精品自产拍在线观看55亚洲| 中文亚洲av片在线观看爽| 久久久久久亚洲精品国产蜜桃av| 人人妻人人爽人人添夜夜欢视频| 日韩国内少妇激情av| 波多野结衣一区麻豆| 国产在线观看jvid| 国产精品电影一区二区三区| 一个人观看的视频www高清免费观看 | 一二三四在线观看免费中文在| 国产精品秋霞免费鲁丝片| 久久精品国产清高在天天线| 亚洲专区字幕在线| 国产欧美日韩一区二区精品| 国产不卡一卡二| 中文字幕av电影在线播放| 亚洲精华国产精华精| 婷婷精品国产亚洲av在线| 最新美女视频免费是黄的| 欧美日韩乱码在线| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美在线二视频| 免费高清在线观看日韩| 日本黄色视频三级网站网址| 婷婷丁香在线五月| 欧美日韩中文字幕国产精品一区二区三区 | 国产熟女xx| 欧美成狂野欧美在线观看| 国产又色又爽无遮挡免费看| 精品高清国产在线一区| 午夜免费激情av| 国产又爽黄色视频| 真人做人爱边吃奶动态| 又大又爽又粗| 亚洲国产精品sss在线观看 | 亚洲男人天堂网一区| 女性被躁到高潮视频| 三级毛片av免费| 精品乱码久久久久久99久播| 日韩免费高清中文字幕av| 69精品国产乱码久久久| 男女下面进入的视频免费午夜 | 在线av久久热| 国产精品一区二区三区四区久久 | 级片在线观看| 日本免费a在线| 12—13女人毛片做爰片一| 成人18禁在线播放| 香蕉丝袜av| 久久久久久人人人人人| 亚洲欧洲精品一区二区精品久久久| 免费观看人在逋| 国产欧美日韩一区二区三区在线| 亚洲欧洲精品一区二区精品久久久| 日本免费a在线| 亚洲专区中文字幕在线| 成人18禁高潮啪啪吃奶动态图| 高清黄色对白视频在线免费看| 亚洲av成人不卡在线观看播放网| 亚洲情色 制服丝袜| 国产成人精品久久二区二区91| 神马国产精品三级电影在线观看 | 国产精品二区激情视频| 女性被躁到高潮视频| 免费在线观看亚洲国产| 美女高潮喷水抽搐中文字幕| 757午夜福利合集在线观看| 操出白浆在线播放| 欧美激情 高清一区二区三区|