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

    基于調(diào)頻層面的風(fēng)電棄風(fēng)分析

    2015-04-06 06:33:58田書欣程浩忠曾平良馬洲俊
    電工技術(shù)學(xué)報(bào) 2015年7期
    關(guān)鍵詞:調(diào)頻電量風(fēng)電場(chǎng)

    田書欣 程浩忠 曾平良 柳 璐 王 侃 馬洲俊

    (1.電力傳輸與功率變換控制教育部重點(diǎn)實(shí)驗(yàn)室(上海交通大學(xué)) 上海 200240 2.中國電力科學(xué)研究院 北京 100192)

    ?

    基于調(diào)頻層面的風(fēng)電棄風(fēng)分析

    田書欣1程浩忠1曾平良2柳 璐1王 侃1馬洲俊1

    (1.電力傳輸與功率變換控制教育部重點(diǎn)實(shí)驗(yàn)室(上海交通大學(xué)) 上海 200240 2.中國電力科學(xué)研究院 北京 100192)

    棄風(fēng)已經(jīng)成為制約風(fēng)電產(chǎn)業(yè)發(fā)展的關(guān)鍵問題。棄風(fēng)造成的原因主要有調(diào)峰調(diào)頻能力和輸電能力不足等。從調(diào)頻層面來看,當(dāng)電網(wǎng)頻率越上限時(shí),可依據(jù)電網(wǎng)調(diào)度部門指令控制風(fēng)電場(chǎng)有功功率輸出,甚至分步切除風(fēng)電機(jī)組。為了分析調(diào)頻層面分鐘級(jí)的應(yīng)急性棄風(fēng),利用有功功率盈余時(shí)風(fēng)電場(chǎng)切機(jī)的調(diào)頻靈敏度指標(biāo),選擇最合適的風(fēng)電切機(jī)策略,并基于頻率越限過渡時(shí)間與負(fù)荷功率變化量的關(guān)系選取棄風(fēng)的調(diào)頻困難典型時(shí)段,再結(jié)合所選的風(fēng)電切機(jī)策略,估算調(diào)頻層面全年規(guī)劃方案棄風(fēng)電量。最后,基于DIgSILENT/PowerFactory分析軟件,利用我國某區(qū)域2020年規(guī)劃電網(wǎng)算例驗(yàn)證了所提方法的可行性和有效性。

    風(fēng)電切機(jī)策略 調(diào)頻 靈敏度分析 棄風(fēng)

    0 引言

    中國正在規(guī)劃建設(shè)9個(gè)千萬千瓦級(jí)風(fēng)電基地,但這些大型風(fēng)電基地一般都遠(yuǎn)離負(fù)荷中心,處于電網(wǎng)末端,且考慮到風(fēng)電的間歇性、隨機(jī)性、波動(dòng)性和不可控性,風(fēng)電并網(wǎng)必然會(huì)給維持系統(tǒng)正常頻率帶來一定困難,進(jìn)而影響電力系統(tǒng)的安全穩(wěn)定運(yùn)行?;诖耍覈L(fēng)電并網(wǎng)技術(shù)標(biāo)準(zhǔn)對(duì)風(fēng)電場(chǎng)有功功率要求為:①風(fēng)電場(chǎng)應(yīng)具有有功功率調(diào)節(jié)能力,配置有功功率控制系統(tǒng),接收并自動(dòng)執(zhí)行調(diào)度部門遠(yuǎn)方發(fā)送的有功功率控制信號(hào);②在風(fēng)電場(chǎng)并網(wǎng)以及風(fēng)速增長過程中,宜控制風(fēng)電場(chǎng)每分鐘有功功率變化率不超過2%~5%;在風(fēng)速降低過程中或超出切機(jī)風(fēng)速情況下,允許有功功率變化率超過該范圍(在超出切機(jī)風(fēng)速情況下,宜分步切除風(fēng)電機(jī)組);③當(dāng)電網(wǎng)頻率高于50.2 Hz時(shí),依據(jù)電網(wǎng)調(diào)度部門指令降低風(fēng)電場(chǎng)有功功率[1]。依據(jù)上述要求,在電網(wǎng)頻率越上限時(shí),切除風(fēng)電機(jī)組,快速降低風(fēng)電場(chǎng)有功功率,確保系統(tǒng)頻率的穩(wěn)定性,從而產(chǎn)生調(diào)頻層面的棄風(fēng)電量。文獻(xiàn)[2,3]從機(jī)組組合或機(jī)組運(yùn)行狀態(tài)分類角度分別提出了風(fēng)電場(chǎng)有功功率控制策略,保證風(fēng)電場(chǎng)有功功率控制的平穩(wěn)性。文獻(xiàn)[4,5]在分析風(fēng)電棄風(fēng)機(jī)理的基礎(chǔ)上,從電網(wǎng)的角度研究了輸電能力不足對(duì)系統(tǒng)棄風(fēng)的影響。文獻(xiàn)[6-8]基于電網(wǎng)調(diào)峰能力提出了分析電網(wǎng)棄風(fēng)情況評(píng)估方法。文獻(xiàn)[9,10]研究了風(fēng)電接入對(duì)系統(tǒng)頻率的影響,并提出相應(yīng)的調(diào)頻策略。但以上文獻(xiàn)均未提及調(diào)頻層面棄風(fēng)電量的估算。

    因此,本文所提的調(diào)頻層面的棄風(fēng)電量估算方法,分析了風(fēng)電場(chǎng)棄風(fēng)與系統(tǒng)一次調(diào)頻運(yùn)行之間的相互作用機(jī)理,實(shí)現(xiàn)了規(guī)劃和運(yùn)行這兩方面問題的耦合,是對(duì)現(xiàn)有考慮風(fēng)電接入下系統(tǒng)運(yùn)行理論和規(guī)劃理論的有效補(bǔ)充和完善。首先給出了計(jì)算各節(jié)點(diǎn)風(fēng)電機(jī)組有功功率調(diào)頻靈敏度的數(shù)學(xué)模型;其次推導(dǎo)了頻率越限過渡時(shí)間與負(fù)荷功率變化量的數(shù)學(xué)關(guān)系式,并設(shè)置調(diào)頻事件場(chǎng)景,選取出現(xiàn)棄風(fēng)的調(diào)頻困難典型時(shí)段;然后,利用蒙特卡洛法計(jì)算每種功率偏差事件在不同典型時(shí)段出現(xiàn)的概率,結(jié)合系統(tǒng)可調(diào)容量,估算調(diào)頻層面的全年規(guī)劃方案棄風(fēng)電量。最后,結(jié)合DIgSILENT/PowerFactory軟件,構(gòu)建含風(fēng)電場(chǎng)的某區(qū)域2020年規(guī)劃電網(wǎng)仿真算例,依據(jù)提出的風(fēng)電機(jī)組切機(jī)策略,實(shí)現(xiàn)對(duì)含風(fēng)電場(chǎng)的電力系統(tǒng)頻率的有效控制,并給出調(diào)頻棄風(fēng)電量的計(jì)算結(jié)果。

    1 切除風(fēng)機(jī)的調(diào)頻靈敏度分析

    風(fēng)電場(chǎng)有功功率控制的準(zhǔn)確性和平穩(wěn)性是高風(fēng)電滲透率電網(wǎng)穩(wěn)定運(yùn)行的基礎(chǔ)。降低不同風(fēng)電場(chǎng)有功功率對(duì)系統(tǒng)頻率[11,12]的影響不同。因而當(dāng)電網(wǎng)頻率越上限時(shí),可按照預(yù)先算出的各風(fēng)電機(jī)組切機(jī)靈敏度[13,14],及時(shí)選擇最合適的切除風(fēng)機(jī)方法來降低風(fēng)電有功功率。下面提出的切風(fēng)機(jī)優(yōu)化方法就是應(yīng)用這一策略。

    靈敏度[15,16]是電力網(wǎng)絡(luò)中某一參量變化量與擾動(dòng)量之比。因而可定義所求取的靈敏度為當(dāng)系統(tǒng)出現(xiàn)較大功率盈余時(shí),有功功率盈余對(duì)所切除風(fēng)機(jī)節(jié)點(diǎn)功率的靈敏系數(shù)。

    設(shè)某時(shí)刻系統(tǒng)出現(xiàn)較大功率盈余,用Pa表示有功功率盈余,假設(shè)系統(tǒng)中有m個(gè)常規(guī)發(fā)電機(jī)節(jié)點(diǎn),n個(gè)風(fēng)力發(fā)電機(jī)節(jié)點(diǎn),k個(gè)負(fù)荷節(jié)點(diǎn),則

    (1)

    由式(1)可看出,求解有功功率盈余對(duì)所切除風(fēng)機(jī)節(jié)點(diǎn)功率的靈敏度必須先求出有功網(wǎng)損靈敏度。有功網(wǎng)損靈敏度的物理意義為:節(jié)點(diǎn)注入功率(包括有功和無功)變化所引起的系統(tǒng)有功損耗改變量。由潮流計(jì)算可知,各節(jié)點(diǎn)注入功率的代數(shù)和就等于網(wǎng)絡(luò)的總有功損耗[17]

    (2)

    式中,i、j為節(jié)點(diǎn)號(hào);Ui、Uj為節(jié)點(diǎn)i、j的電壓幅值,kV;Gij、Bij分別為節(jié)點(diǎn)導(dǎo)納矩陣元素的實(shí)部與虛部,S;θij為節(jié)點(diǎn)i,j的電壓相角差值,θij=θi-θj。

    假設(shè)系統(tǒng)中只有平衡節(jié)點(diǎn)和實(shí)施切風(fēng)機(jī)措施的節(jié)點(diǎn)發(fā)生功率變化,若風(fēng)機(jī)所在節(jié)點(diǎn)i線路實(shí)施切除后,該節(jié)點(diǎn)有功功率變化為ΔPi,無功功率變化為ΔQi,由于平衡節(jié)點(diǎn)的調(diào)節(jié)作用,使網(wǎng)絡(luò)總有功損耗也發(fā)生相應(yīng)變化。i節(jié)點(diǎn)有功功率的減少將使平衡節(jié)點(diǎn)注入有功功率增加,并使網(wǎng)絡(luò)損耗發(fā)生變化。

    可推導(dǎo)風(fēng)機(jī)所在節(jié)點(diǎn)有功功率變化對(duì)網(wǎng)損的靈敏度為

    (3)

    可推導(dǎo)風(fēng)機(jī)所在節(jié)點(diǎn)無功功率變化對(duì)網(wǎng)損的靈敏度為

    (4)

    式中,s為平衡節(jié)點(diǎn)號(hào);i為風(fēng)機(jī)所在節(jié)點(diǎn),i=1,2,…,n。

    聯(lián)立式(3)和式(4),并進(jìn)行相應(yīng)轉(zhuǎn)換,可得到

    (5)

    在所需切除點(diǎn)附近將系統(tǒng)潮流方程線性化,可得

    (6)

    可得到

    (7)

    可以看出,式(7)中系數(shù)矩陣就是極坐標(biāo)牛頓-拉夫遜法潮流計(jì)算中的雅可比矩陣。結(jié)合式(1),則有功功率盈余對(duì)風(fēng)機(jī)所接節(jié)點(diǎn)功率的靈敏度為

    (8)

    可得到

    (9)

    (10)

    基于式(9)和式(10)得到切除風(fēng)機(jī)靈敏度,其數(shù)值越大的風(fēng)機(jī)機(jī)組節(jié)點(diǎn),對(duì)快速恢復(fù)系統(tǒng)頻率的影響越大。因此可根據(jù)靈敏度大小,確定調(diào)頻過程中所需切除風(fēng)機(jī)的順序。

    2 調(diào)頻困難時(shí)段棄風(fēng)電量的計(jì)算

    2.1 頻率越限過渡時(shí)間的計(jì)算

    電力系統(tǒng)中發(fā)電機(jī)組的投切或負(fù)荷的增減,使某時(shí)段內(nèi)系統(tǒng)中的能量供求關(guān)系失去平衡,系統(tǒng)頻率將從正常的穩(wěn)態(tài)值過渡到另一個(gè)使能量供求達(dá)到新的平衡的穩(wěn)態(tài)值,或失去穩(wěn)定即發(fā)生頻率崩潰事故,該時(shí)段即調(diào)頻困難時(shí)段。在含風(fēng)電接入的電力系統(tǒng)中,風(fēng)速上升或有功負(fù)荷下降所造成的頻率波動(dòng)越限的時(shí)間,即頻率越限過渡時(shí)間,對(duì)計(jì)算風(fēng)電棄風(fēng)電量具有重要意義。風(fēng)速上升或負(fù)荷下降造成的頻率變化與系統(tǒng)中發(fā)電機(jī)組備用容量大小、負(fù)荷調(diào)節(jié)效應(yīng)、發(fā)電機(jī)組的機(jī)械慣性和負(fù)荷的機(jī)械慣性等多個(gè)因素有關(guān),暫不考慮負(fù)荷隨電壓的變化、AGC控制策略等,則可基于圖1給出的負(fù)荷下降情況下電力系統(tǒng)功率-頻率靜態(tài)特性曲線分析其調(diào)頻特性。

    圖1 電力系統(tǒng)功率-頻率靜態(tài)特性Fig.1 Power-frequency static characteristic curve

    由于風(fēng)力發(fā)電出力完全由天氣狀況決定,具有隨機(jī)性、間歇性的特點(diǎn),在風(fēng)力接入電網(wǎng)研究中,通常將風(fēng)力電站出力視為負(fù)的負(fù)荷。因此計(jì)算頻率越限時(shí)間所用的負(fù)荷功率PL可表示為

    PL=PGL-PWL

    (11)

    式中,PGL為電力系統(tǒng)中的普通負(fù)荷,MW;PWL為電力系統(tǒng)中風(fēng)電機(jī)組有功出力,MW。

    Δf∞=f1-f0

    (12)

    式中,f1為c點(diǎn)對(duì)應(yīng)的新的穩(wěn)態(tài)運(yùn)行頻率,Hz;f0為起始時(shí)刻電力系統(tǒng)頻率的正常運(yùn)行值,Hz。

    在風(fēng)速上升或負(fù)荷下降后的短時(shí)間內(nèi)

    ω=ω0+Δω

    (13)

    PL=PL0+ΔPL

    (14)

    式中,ω為系統(tǒng)發(fā)電機(jī)轉(zhuǎn)速,rad/s;ω0為系統(tǒng)發(fā)電機(jī)轉(zhuǎn)速的正常運(yùn)行值,rad/s;Δω為系統(tǒng)發(fā)電機(jī)轉(zhuǎn)速的變化量,rad/s;PL為負(fù)荷功率值,MW;PL0為不平衡初始的負(fù)荷功率,MW;ΔPL為不平衡出現(xiàn)后的負(fù)荷功率變化量,MW。又因頻率隨時(shí)間的變化率可記為

    (15)

    式中,ω*為系統(tǒng)發(fā)電機(jī)轉(zhuǎn)速標(biāo)幺值;Δω*為系統(tǒng)發(fā)電機(jī)轉(zhuǎn)速的變化量標(biāo)幺值;Δf*為頻率變化量的標(biāo)幺值。

    因而,電力系統(tǒng)頻率變化時(shí)等值機(jī)組的運(yùn)動(dòng)方程表達(dá)式可寫為

    (16)

    式中,TJ為系統(tǒng)等效機(jī)組慣性時(shí)間常數(shù);PT*為以系統(tǒng)發(fā)電機(jī)總功率額定值PGN為基準(zhǔn)的原動(dòng)機(jī)功率標(biāo)幺值;PL*為以系統(tǒng)發(fā)電機(jī)總功率額定值PGN為基準(zhǔn)的負(fù)荷功率標(biāo)幺值。設(shè)風(fēng)速上升或負(fù)荷下降時(shí)系統(tǒng)出現(xiàn)不平衡功率ΔPa=ΔPL,且出現(xiàn)ΔPL后任意時(shí)刻t的系統(tǒng)頻率值為ft,負(fù)荷功率PL=PLt,可參見圖1得

    PT*-PL*=ΔPL*-KL*Δf*

    (17)

    式中,KL*為考慮風(fēng)電接入的負(fù)載調(diào)節(jié)效應(yīng)系數(shù)標(biāo)幺值;Δf*為t時(shí)刻系統(tǒng)頻率與起始頻率之差的標(biāo)幺值,即Δf*=(ft-f0)/f0。 公式推導(dǎo)過程中ω和f均是以事件開始發(fā)生時(shí)的ω0和f0為基準(zhǔn)。將式(17)代入式(16)可得

    (18)

    式中,Tf為頻率上升過程的時(shí)間常數(shù),Tf與整個(gè)系統(tǒng)的轉(zhuǎn)動(dòng)部分機(jī)械慣性時(shí)間常數(shù)及負(fù)載調(diào)節(jié)效應(yīng)系數(shù)相關(guān),即Tf=TJ/KL*;f*為出現(xiàn)不平衡功率后任意時(shí)刻t的系統(tǒng)頻率標(biāo)幺值;ΔPL*為不平衡功率ΔPL的標(biāo)幺值,且ΔPL<0。

    式(18)表明,簡化的系統(tǒng)動(dòng)態(tài)頻率特性是按指數(shù)規(guī)律衰減的曲線。令Δf∞*=ΔPL*/KL*, 其表達(dá)式為

    f*=f0*-Δf∞*(1-e-t/Tf)

    (19)

    式中,f0*為不平衡初始時(shí)刻(t=0)的系統(tǒng)頻率標(biāo)幺值;Δf∞*為系統(tǒng)新穩(wěn)態(tài)頻率與起始頻率之差的標(biāo)幺值。

    由式(19)易知

    (20)

    將Δf∞*=ΔPL*/KL*代入式(20),可得

    (21)

    由以上公式可計(jì)算風(fēng)速上升和負(fù)荷下降所引起頻率越限過渡時(shí)間tv。

    2.2 調(diào)頻層面棄風(fēng)電量的估算

    從電網(wǎng)運(yùn)行角度,風(fēng)電的棄風(fēng)是由調(diào)峰調(diào)頻能力不足造成的。調(diào)峰層面[18]的棄風(fēng)電量計(jì)算通常是在負(fù)荷預(yù)測(cè)的基礎(chǔ)上,經(jīng)過修改機(jī)組組合而得到;而調(diào)頻層面的棄風(fēng)是在考慮負(fù)荷突變的基礎(chǔ)上得到的[19],兩者之間沒有交集,所以全年的總體棄風(fēng)電量可分別從調(diào)峰層面和調(diào)頻層面進(jìn)行計(jì)算。下面著重對(duì)調(diào)頻層面全年棄風(fēng)電量的估算方法進(jìn)行研究。

    電力系統(tǒng)的調(diào)頻場(chǎng)景有很多種,其中可能需要棄風(fēng)的場(chǎng)景有3類:風(fēng)速上升導(dǎo)致的風(fēng)電場(chǎng)輸出功率的增加、負(fù)荷下降、風(fēng)速上升與負(fù)荷下降同時(shí)發(fā)生。在調(diào)頻層面棄風(fēng)分析中,可將風(fēng)電場(chǎng)輸出功率作為負(fù)的負(fù)荷處理,所以以上3個(gè)場(chǎng)景均可認(rèn)為是負(fù)荷下降,但考慮到調(diào)頻事件中存在更多的不確定性,風(fēng)電和負(fù)荷波動(dòng)性和預(yù)測(cè)誤差更大,因而可采用組合的方法將各種不確定性因素的可能取值組成一個(gè)未來可能場(chǎng)景,比如負(fù)荷波動(dòng)和風(fēng)電功率波動(dòng)組成了一個(gè)調(diào)頻場(chǎng)景,在場(chǎng)景中的各種不確定性信息具有確定的數(shù)值,如負(fù)荷波動(dòng)1%和風(fēng)電功率波動(dòng)1%等。每個(gè)場(chǎng)景代表了某種概率下可能的未來環(huán)境,通過合理選擇一系列的場(chǎng)景可近似表示未來可能出現(xiàn)的所有不確定性環(huán)境。

    2.2.1 調(diào)頻事件場(chǎng)景設(shè)置

    在調(diào)頻事件的場(chǎng)景設(shè)置中,考慮了以下幾個(gè)因素:

    1)頻率越限過渡時(shí)間只與負(fù)荷以及負(fù)荷的變化率有關(guān)。

    2)負(fù)荷變化率超過10%的概率較低,所以只選擇負(fù)荷下降10%以內(nèi)作為計(jì)算棄風(fēng)電量的調(diào)頻事件。

    3)由于負(fù)荷變化率為連續(xù)變量,不易計(jì)算,所以將連續(xù)的負(fù)荷變化率離散化處理。

    4)各時(shí)段負(fù)荷值雖為離散值,但由于負(fù)荷值較多,所以將小時(shí)級(jí)負(fù)荷曲線分區(qū)間處理。

    基于以上幾點(diǎn),調(diào)頻事件設(shè)置方式為:根據(jù)負(fù)荷水平對(duì)全年8 760 h劃分為10個(gè)區(qū)間,每個(gè)區(qū)間由該區(qū)間聚類數(shù)值代表,每個(gè)區(qū)間的負(fù)荷變化率離散化為10個(gè)調(diào)頻事件,因此共有100個(gè)調(diào)頻事件。

    2.2.2 計(jì)算系統(tǒng)可調(diào)容量

    負(fù)荷下降時(shí),系統(tǒng)的頻率不會(huì)突然上升,而是會(huì)以一定的規(guī)律上升,最終達(dá)到一個(gè)穩(wěn)態(tài)值。因此可通過仿真找到系統(tǒng)頻率越限過渡時(shí)間tv。在這個(gè)越限過渡時(shí)間內(nèi),所有機(jī)組的可調(diào)容量之和大于或等于ΔP時(shí),認(rèn)為該時(shí)刻不需要棄風(fēng);反之,則需要按照靈敏度優(yōu)先順序逐步棄風(fēng)。所有調(diào)頻機(jī)組可調(diào)容量之和的計(jì)算公式為

    (22)

    式中,N為火電機(jī)組和水電機(jī)組的數(shù)量;M為抽水蓄能機(jī)組的數(shù)量;Ramp為各類型機(jī)組的下坡速率,%/min;tv為頻率越限時(shí)間,s;pi、pj為各機(jī)組該時(shí)刻出力,MW;pi,min為機(jī)組最小出力,MW;PN為抽水蓄能機(jī)組的額定出力,MW。

    2.2.3 計(jì)算棄風(fēng)電量

    1)單一事件發(fā)生在某一時(shí)刻下棄風(fēng)電量的計(jì)算

    在調(diào)頻時(shí)間達(dá)到tv時(shí),系統(tǒng)各機(jī)組出力無法與負(fù)荷相平衡,所以此時(shí)需要棄風(fēng)。風(fēng)力發(fā)電的燃料費(fèi)用為零,所以應(yīng)多利用風(fēng)電,減小棄風(fēng)電量。在tv時(shí)刻之后,系統(tǒng)各機(jī)組仍需以最大的下坡速率[20]進(jìn)行降速運(yùn)行,同時(shí)被切除的風(fēng)機(jī)以相應(yīng)速率并網(wǎng),直到所有機(jī)組減小出力之和與ΔP相等,此時(shí)所有被切除的風(fēng)機(jī)均可并網(wǎng)運(yùn)行。

    圖2中的曲線y表示各時(shí)刻系統(tǒng)各機(jī)組出力之和。在時(shí)刻tv之后,曲線y上的點(diǎn)與(P0-ΔP)之差為此時(shí)的棄風(fēng)功率,因此圖2中陰影部分面積為棄風(fēng)電量。

    圖2 單一調(diào)頻事件棄風(fēng)電量計(jì)算圖示Fig.2 Wind power curtailment in a single frequency adjustment event

    圖2中曲線y的斜率即為機(jī)組的總爬坡率。在tv之后,如果所有可調(diào)機(jī)組的出力均為達(dá)到其最小出力,那么曲線y即為一直線;如果某些可調(diào)機(jī)組的出力達(dá)到最小,則不能參與調(diào)節(jié),那么曲線y為分段函數(shù)。下面給出曲線y為一直線時(shí)的陰影部分的計(jì)算公式

    Cw=0.5×(ΔP-Pad)2/Ramptotal

    (23)

    式中,ΔP為風(fēng)電或負(fù)荷突變引起的凈負(fù)荷的變化量,MW;Pad為調(diào)頻機(jī)組的可調(diào)容量,MW;Ramptotal為各類型發(fā)電機(jī)組的總下坡率,%/min。

    2)全年棄風(fēng)電量的計(jì)算

    棄風(fēng)電量的計(jì)算過程如下:

    1)計(jì)算出每個(gè)事件在某一調(diào)頻困難典型時(shí)段發(fā)生時(shí)的棄風(fēng)電量。

    2)計(jì)算每個(gè)事件發(fā)生的概率。

    3)計(jì)算所有事件在某一調(diào)頻困難典型時(shí)段發(fā)生時(shí)產(chǎn)生的棄風(fēng)電量。

    4)將全年8 760h中所有調(diào)頻困難典型時(shí)段的數(shù)據(jù)相加即得全年棄風(fēng)電量。

    全年調(diào)頻層面棄風(fēng)電量為

    (24)

    3 算例分析

    3.1 切除風(fēng)機(jī)靈敏度的計(jì)算

    以我國某區(qū)域2020年規(guī)劃電網(wǎng)為算例,應(yīng)用所提的調(diào)頻層面棄風(fēng)電量的計(jì)算方法對(duì)其進(jìn)行評(píng)估。該區(qū)域2020年規(guī)劃電網(wǎng)主網(wǎng)架共有86個(gè)節(jié)點(diǎn),其系統(tǒng)示意圖如圖3所示。圖3中電源基本情況為火電147 000MW,水電21 500MW,核電2 000MW,抽水蓄能1 600MW,風(fēng)電23 600MW。其中,風(fēng)電分別從節(jié)點(diǎn)58、60、65的匯集站接入主網(wǎng)。

    圖3 某區(qū)域2020年規(guī)劃電網(wǎng)86節(jié)點(diǎn)系統(tǒng)圖Fig.3 China’s a certain regional power network planning of 2020 (86-bus system)

    基于圖2給出的規(guī)劃電網(wǎng)86節(jié)點(diǎn)系統(tǒng)圖,結(jié)合DIgSILENT/PowerFactory電力系統(tǒng)分析軟件,建立了該系統(tǒng)的調(diào)頻仿真模型,如圖4所示。根據(jù)接入系統(tǒng)的風(fēng)電基地的功率匯集方案,將各風(fēng)電基地的每個(gè)風(fēng)電場(chǎng)均等效為一臺(tái)風(fēng)機(jī),風(fēng)機(jī)選用2MW的雙饋感應(yīng)電機(jī)模型,根據(jù)每個(gè)風(fēng)電場(chǎng)實(shí)際裝機(jī)容量設(shè)置相關(guān)參數(shù)。

    圖4 基于DIgSILENT的86節(jié)點(diǎn)系統(tǒng)仿真模型Fig.4 Simulation model of 86-bus system based on DIgSILENT/PowerFactory

    根據(jù)式(1)~式(10)計(jì)算風(fēng)電匯入節(jié)點(diǎn)58、60、65的切機(jī)調(diào)頻靈敏度,結(jié)果如表1所示。

    表1 風(fēng)電匯入節(jié)點(diǎn)的切機(jī)調(diào)頻靈敏度Tab.1 The frequency adjustment sensitivities of the bus integrated wind power

    3.2 頻率越限過渡時(shí)間的計(jì)算

    由式(21)可知,在已知Tf、KL*等參數(shù)條件下,頻率越限過渡時(shí)間tv主要取決于頻率變化量Δf*和不平衡功率ΔPL*。假設(shè)從不平衡功率出現(xiàn)時(shí)刻開始經(jīng)tv時(shí)間后ft=50.2 Hz,即Δf*=0.004, 則頻率越限過渡時(shí)間tv僅與負(fù)荷功率變化量ΔPL有關(guān)。同時(shí)負(fù)荷功率變化量ΔPL又可表示為

    ΔPL=x%PLN

    (25)

    式中,PLN為負(fù)荷功率參考值,MW;x%為負(fù)荷功率百分比,通過在圖4中的SimulationEvents設(shè)置不同x%的LoadEvent,即x=-1、x=-2、…、x=-10等。

    因而可根據(jù)圖5給出的電網(wǎng)規(guī)劃年內(nèi)各小時(shí)段的負(fù)荷預(yù)測(cè)曲線計(jì)算x%PLN,并基于Kmeans方法對(duì)調(diào)頻時(shí)段區(qū)間負(fù)荷水平進(jìn)行等效聚類,得到每類區(qū)間的典型時(shí)段,結(jié)果參見表2。

    根據(jù)表2給出的典型時(shí)段,結(jié)合圖4中的主網(wǎng)仿真模型,輸入典型時(shí)段的機(jī)組安排及負(fù)荷水平以及該時(shí)段下的初始風(fēng)速,分別對(duì)x=-1、x=-2、…、x=-10等Load Event進(jìn)行仿真,記錄所需要的頻率越限過渡時(shí)間,參見表3。

    圖5 某區(qū)域2020年規(guī)劃電網(wǎng)年負(fù)荷預(yù)測(cè)曲線Fig.5 The year load curve of a certain regional power network planning of 2020

    表2 典型時(shí)段表格Tab.2 Typical difficult time interval

    表3 頻率越限時(shí)間tv

    由表3可知,該區(qū)域規(guī)劃電網(wǎng)僅在系統(tǒng)功率偏差變化較大時(shí),即x=-9、x=-10時(shí)發(fā)生頻率越限,其余功率偏差事件,頻率仍可保持在穩(wěn)定范圍內(nèi)。

    3.3 調(diào)頻層面棄風(fēng)電量的計(jì)算

    假設(shè)風(fēng)速服從威布爾分布而負(fù)荷增長服從正態(tài)分布,然后利用已有的全年的風(fēng)電和負(fù)荷的預(yù)測(cè)數(shù)據(jù)進(jìn)行分布的參數(shù)估計(jì),并采用蒙特卡洛模擬法計(jì)算每種功率偏差事件在不同時(shí)段出現(xiàn)的概率,可參見表4。

    對(duì)于調(diào)頻困難典型時(shí)段內(nèi)的每種功率偏差事件,按照切機(jī)靈敏度順序依次逐步棄風(fēng),并基于式(22)~式(24),計(jì)算調(diào)頻層面的全年棄風(fēng)電量。結(jié)合文獻(xiàn)[21]的調(diào)峰層面的風(fēng)電棄風(fēng)機(jī)理,計(jì)算了調(diào)峰層面的全年棄風(fēng)電量,并將二者進(jìn)行了比較。具體結(jié)果參見表5。

    表4 調(diào)頻事件概率Tab.4 Probability of power deviation event

    表5 全年棄風(fēng)電量匯總Tab.5 Annual total wind energy curtailment

    基于表5的全年棄風(fēng)電量的估算結(jié)果可知,該規(guī)劃年全年調(diào)頻棄風(fēng)電量占全年棄風(fēng)電量的比例為28.9%,因而對(duì)調(diào)頻層面的棄風(fēng)電量進(jìn)行評(píng)估,將有助于促進(jìn)風(fēng)電消納和系統(tǒng)頻率的恢復(fù)。

    為了確定調(diào)頻棄風(fēng)中按照調(diào)頻靈敏度選擇不同匯入節(jié)點(diǎn)風(fēng)電場(chǎng)線路切除策略對(duì)頻率的影響,基于表3和表4,針對(duì)出現(xiàn)頻率越限的x=-9、x=-10功率偏差事件,并分別選取概率最高的第2 019小時(shí)和第1 236小時(shí)典型時(shí)段,對(duì)不同風(fēng)電匯入節(jié)點(diǎn)進(jìn)行棄風(fēng)調(diào)頻仿真。

    1)第2 019小時(shí)棄風(fēng)調(diào)頻仿真

    具體仿真參數(shù)設(shè)置為:風(fēng)速:11m/s,隨機(jī);仿真時(shí)間:50s;事件:5s時(shí)發(fā)生x=-10功率偏差;條件:槳距角不動(dòng)作,切1 000MW風(fēng)機(jī)(當(dāng)頻率偏差超過0.2Hz)。仿真結(jié)果如圖6~圖8所示。

    2)第1 236小時(shí)棄風(fēng)調(diào)頻仿真

    具體仿真參數(shù)設(shè)置為:風(fēng)速:11m/s,隨機(jī);仿真時(shí)間:50s;事件:5s時(shí)發(fā)生x=-9功率偏差;條件:槳距角不動(dòng)作,切900MW風(fēng)機(jī)(當(dāng)頻率偏差超過0.2Hz)。仿真結(jié)果如圖9~圖11所示。

    可看出,切除靈敏度數(shù)值越大的58節(jié)點(diǎn)風(fēng)電場(chǎng)所在線路,更有利于系統(tǒng)頻率的平穩(wěn)恢復(fù)。

    圖6 切除58節(jié)點(diǎn)1 000 MW風(fēng)機(jī)的調(diào)頻仿真結(jié)果Fig.6 The frequency adjustment simulation result for curtailing 1 000 MW wind generators at bus 58

    圖7 切除60節(jié)點(diǎn)1 000 MW風(fēng)機(jī)的調(diào)頻仿真結(jié)果Fig.7 The frequency adjustment simulation result for curtailing 1 000 MW wind generators at bus 60

    圖8 切除65節(jié)點(diǎn)1 000 MW風(fēng)機(jī)的調(diào)頻仿真結(jié)果Fig.8 The frequency adjustment simulation result for curtailing 1 000 MW wind generators at bus 65

    圖9 切除58節(jié)點(diǎn)900 MW風(fēng)機(jī)的調(diào)頻仿真結(jié)果Fig.9 The frequency adjustment simulation result for curtailing 900 MW wind generators at bus 58

    圖10 切除60節(jié)點(diǎn)900 MW風(fēng)機(jī)的調(diào)頻仿真結(jié)果Fig.10 The frequency adjustment simulation result for curtailing 900 MW wind generators at bus 60

    圖11 切除65節(jié)點(diǎn)900 MW風(fēng)機(jī)的調(diào)頻仿真結(jié)果Fig.11 The frequency adjustment simulation result for curtailing 900 MW wind generators at bus 65

    4 結(jié)論

    調(diào)頻層面棄風(fēng)考慮的是風(fēng)電的分鐘級(jí)波動(dòng),相比較調(diào)峰層面有更多的不確定性,波動(dòng)性和預(yù)測(cè)誤差更大。通過聯(lián)合DIgSILENT/PowerFactory軟件,仿真模擬基于切機(jī)靈敏度的棄風(fēng)調(diào)頻過程,給出頻率越限過渡時(shí)間,進(jìn)而提出調(diào)頻層面全年棄風(fēng)電量的估算方法,并以我國某區(qū)域2020年規(guī)劃電網(wǎng)為算例,計(jì)算了全年調(diào)頻層面棄風(fēng)電量及其占全年棄風(fēng)電量的比例。算例結(jié)果表明全年調(diào)頻層面棄風(fēng)電量對(duì)全年總體棄風(fēng)電量影響較大。雖然2020年該區(qū)域主網(wǎng)規(guī)劃了多條大容量長距離的外送線路,其中包含了多條特高壓交直流線路,從一定程度上很好地解決了2020年大規(guī)模風(fēng)電的消納問題,但全年調(diào)頻層面棄風(fēng)電量評(píng)估方法和基于調(diào)頻靈敏度切除風(fēng)機(jī)策略仍有助于合理配置系統(tǒng)的調(diào)頻容量,從而將整體的風(fēng)電消納提高到一個(gè)較合理的規(guī)模,可為解決風(fēng)電接入后電力系統(tǒng)規(guī)劃和運(yùn)行的耦合問題提供一定的借鑒和參考。另外,所提的調(diào)頻層面棄風(fēng)電量估算方法對(duì)風(fēng)電調(diào)頻棄風(fēng)問題只是做了些探索性的、初步的研究工作,尚有很多問題有待進(jìn)一步深入研究。

    [1] 國家電網(wǎng)公司企業(yè)標(biāo)準(zhǔn):Q/GDW392-2009,風(fēng)電場(chǎng)接入電網(wǎng)技術(shù)規(guī)定[S].

    [2]MoyanoCF,LopesJAP.Unitcommitmentanddispatchstrategiesforawindpark[C].InternationalConferenceonPowerEngineering,EnergyandElectricalDrives,Setubal,Portugal,2007:12-17.

    [3] 鄒見效,李丹,鄭剛,等.基于機(jī)組狀態(tài)分類的風(fēng)電場(chǎng)有功功率控制策略[J].電力系統(tǒng)自動(dòng)化,2011,35(24):28-32.ZouJianxiao,LiDan,ZhengGang,etal.Anactivepowercontrolschemeforwindfarmsbasedonstateclassificationalgorithm[J].AutomationofElectricPowerSystems,2011,35(24):28-32.

    [4]BurkeDJ,O’MalleyMJ.Transmissionconnectedwindcurtailmentwithincreasingwindcapacityconnection[C].IEEEPower&EnergySocietyGeneralMeeting,Calgary,AB,2009:1-5.

    [5]BurkeDanielJ,O’MalleyMJ.FactorsInfluencingWindEnergyCurtailment[J].IEEETransactionsonSustainableEnergy,2011,2(2):185-193.

    [6] 劉新東,方科,陳煥遠(yuǎn),等.利用合理?xiàng)夛L(fēng)提高大規(guī)模風(fēng)電消納能力的理論研究[J].電力系統(tǒng)保護(hù)與控制,2012,40(6):35-39.LiuXindong,F(xiàn)angKe,ChenHuanyuan,etal.Researchonrationalwindpowercastingtheoryforlarge-scalewindpowerintegrationimprovement[J].PowerSystemProtectionandControl,2012,40(6):35-39.

    [7] 侯婷婷,婁素華,吳耀武,等.含大型風(fēng)電場(chǎng)的電力系統(tǒng)調(diào)峰運(yùn)行特性分析[J].電工技術(shù)學(xué)報(bào),2013,28(5):105-111.HouTingting,LouSuhua,WuYaowu,etal.Analysisonpeakloadregulationoperationcharacteristicsofpowersystemintegratedwithlarge-scalewindpower[J].TransactionsofChinaElectrotechnicalSociety,2013,28(5):105-111.

    [8]ZhangNing,KangChongqing,KirschenDS,etal.Planningpumpedstoragecapacityforwindpowerintegration[J].IEEETransactionsonSustainableEnergy,2013,4(2):393-401.

    [9] 吳子雙,于繼來,彭喜云.高風(fēng)速段次優(yōu)功率追蹤方式的風(fēng)電調(diào)頻方法[J].電工技術(shù)學(xué)報(bào),2013,28(5):112-119.WuZishuang,YuJilai,PengXiyun.DFIG’sfrequencyregulationmethodonlyforhighwindspeedwithsuboptimalpowertracking[J].TransactionsofChinaElectrotechnicalSociety,2013,28(5):112-119.

    [10]HaqueME,PaulS,SheikhMRI.Gridfrequencyanalysiswiththeissueofhighwindpowerpenetration[C].InternationalConferenceonElectricalInformationandCommunicationTechnology(EICT),Khulna,2014:1-7.

    [11]VaidyaA,ChowdhuryBH,ChamanaM.Enablingwindpowerplantswithfrequencyandvoltageregulationcapabilitybyforcedcurtailment[C].NorthAmericanPowerSymposium(NAPS),Champaign,IL,USA,2012:1-6.

    [12]KleinSA.Accomodatingwindpowervariabilitybyreal-timeadjustmentofcontrollableloads[C].IEEEPowerandEnergySocietyGeneralMeeting,Pittsburgh,PA,USA,2008:1-5.

    [13]GuYingzhong,XieLe.Fastsensitivityanalysisapproachtoassessingcongestioninducedwindcurtailment[J].IEEETransactionsonPowerSystems,2014,29(1):101-110.

    [14]BenidrisMohammed,MitraJoydeep.Reliabilityandsensitivityanalysisofcompositepowersystemsunderemissionconstraints[J].IEEETransactionsonPowerSystems,2014,29(1):404-412.

    [15]洪蘆誠,石立寶,姚良忠,等.計(jì)及風(fēng)電場(chǎng)發(fā)電功率不確定性的電力系統(tǒng)模糊潮流[J].電工技術(shù)學(xué)報(bào),2010,25(8):116-122.HongLucheng,ShiLibao,YaoLiangzhong,etal.Fuzzymodellingandsolutionofloadflowincorporatinguncertaintiesofwindfarmgeneration[J].TransactionsofChinaElectrotechnicalSociety,2010,25(8):116-122.

    [16]荊朝霞,胡仁.電能質(zhì)量評(píng)估靈敏度分析及在改進(jìn)密切值法中的應(yīng)用[J].電力系統(tǒng)保護(hù)與控制,2013,41(9):60-66.JingZhaoxia,HuRen.Sensitivityanalysisforpowerqualityevaluationanditsapplicationinmodifiedclosevaluemethod[J].PowerSystemProtectionandControl,2013,41(9):60-66.

    [17]HuWei,DongShimeng,LuQiuyu,etal.CoordinatedoptimizingcontrolstrategybetweenpowergeneratorsandpowergridinInnerMongolia[C].UniversitiesPowerEngineeringConference,LondonUK,2012:1-7.

    [18]劉暢,吳浩,高長征,等.風(fēng)電消納能力分析方法的研究[J].電力系統(tǒng)保護(hù)與控制,2014,42(4):61-66.LiuChang,WuHao,GaoChangzheng,etal.Studyonanalysismethodofaccommodatedcapacityforwindpower[J].PowerSystemProtectionandControl,2014,42(4):61-66.

    [19]Wang-HansenM,JosefssonR,MehmedovicH.Frequencycontrollingwindpowermodelingofcontrolstrategies[J].IEEETransactionsonSustainableEnergy,2013,4(4):954-959.

    [20]KalantariA,GalianaFD.Theimpactofwindpowervariabilityandcurtailmentonrampingrequirements[C].IEEE/PESTransmissionandDistributionConferenceandExposition,SaoPaulo,USA,2010:133-138.

    [21]呂泉,王偉,韓水,等.基于調(diào)峰能力分析的電網(wǎng)棄風(fēng)情況評(píng)估方法[J].電網(wǎng)技術(shù),2013,37(7):1887-1894.LvQuan,WangWei,HanShui,etal.Anewevaluationmethodforwindpowercurtailmentbasedonanalysisofsystemregulationcapability[J].PowerSystemTechnology,2013,37(7):1887-1894.

    程浩忠 男,1962生,博士,教授,博士生導(dǎo)師,研究方向?yàn)殡娏ο到y(tǒng)規(guī)劃、電壓穩(wěn)定、電能質(zhì)量等。

    Analysis on Wind Power Curtailment at Frequency Adjustment Level

    TianShuxin1ChengHaozhong1ZengPingliang2LiuLu1WangKan1MaZhoujun1

    (1.Key Laboratory of Control of Power Transmission and Conversion Ministry of Education(Shanghai Jiao Tong University) Shanghai 200240 China 2.China Electric Power Research Institute Beijing 100192 China)

    Wind power curtailment has been one of the key restricting factors behind the development of the wind power industry.The main reasons for wind power curtailment are inadequate peak and frequency adjustment capability as well as inadequate transmission capacity.As seen from the aspect of the frequency adjustment,the active power output of the wind farm should be reduced step by step according to the command from the dispatch department of the power grid when the grid frequency exceeds the upper limit.To study urgent wind power curtailment at the minute-level,the most appropriate shedding method for wind generators is selected based on the frequency adjustment sensitivity index of wind generators’ tripping as lopsided power exists.Typical difficult time intervals of the frequency adjustment are chosen according to the relationship between frequency overshoot time and load variation.Combined with the above wind generator tripping strategies,annual wind power curtailment in the planning scheme at the frequency adjustment level is estimated.Based on DIgSILENT/PowerFactory software,a certain regional power network planning in China of 2020 is used to verify the correctness and validity of the proposed method.

    Wind generator tripping strategy,power grid frequency,sensitivity analysis,wind power curtailment

    國家自然科學(xué)基金(51337005),國家電網(wǎng)公司科技項(xiàng)目(SGCC-MPLG018-2012)和國家電網(wǎng)公司“千人計(jì)劃”專項(xiàng)(XT71-12-028)資助項(xiàng)目。

    2014-11-15 改稿日期2015-01-26

    TM73

    田書欣 男,1985年生,博士研究生,研究方向?yàn)榭紤]風(fēng)電的輸電網(wǎng)規(guī)劃。(通信作者)

    猜你喜歡
    調(diào)頻電量風(fēng)電場(chǎng)
    電量越低越透明的手機(jī)
    考慮頻率二次跌落抑制的風(fēng)火聯(lián)合一次調(diào)頻控制
    能源工程(2021年5期)2021-11-20 05:50:42
    基于PSS/E的風(fēng)電場(chǎng)建模與動(dòng)態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    四川2018年7月轉(zhuǎn)讓交易結(jié)果:申報(bào)轉(zhuǎn)讓電量11.515 63億千瓦時(shí)
    電量隔離傳感器測(cè)試儀的研制
    含風(fēng)電場(chǎng)電力系統(tǒng)的潮流計(jì)算
    調(diào)頻發(fā)射機(jī)技術(shù)改造
    調(diào)頻激勵(lì)器干擾的排除方法
    調(diào)頻引信中噪聲調(diào)幅干擾的自適應(yīng)抑制
    探求風(fēng)電場(chǎng)的遠(yuǎn)景
    風(fēng)能(2015年9期)2015-02-27 10:15:25
    18禁在线无遮挡免费观看视频| 少妇丰满av| 三级国产精品欧美在线观看| 久久久欧美国产精品| 性色av一级| 91精品国产国语对白视频| 亚洲国产av新网站| 91精品一卡2卡3卡4卡| 精品久久久噜噜| 晚上一个人看的免费电影| 99九九在线精品视频 | 18禁在线播放成人免费| 最近最新中文字幕免费大全7| 亚洲av中文av极速乱| 久久热精品热| 国产成人精品无人区| 一级毛片黄色毛片免费观看视频| 看免费成人av毛片| 人妻一区二区av| 国产精品不卡视频一区二区| 日韩视频在线欧美| 午夜影院在线不卡| videos熟女内射| 亚洲av不卡在线观看| 青春草国产在线视频| 亚洲欧美一区二区三区国产| 97精品久久久久久久久久精品| 乱系列少妇在线播放| 乱系列少妇在线播放| 美女cb高潮喷水在线观看| 热99国产精品久久久久久7| 欧美日韩国产mv在线观看视频| 美女福利国产在线| av一本久久久久| 国产精品女同一区二区软件| 亚洲国产精品一区三区| 国产精品人妻久久久影院| 各种免费的搞黄视频| 大码成人一级视频| 啦啦啦在线观看免费高清www| 精品午夜福利在线看| 精品午夜福利在线看| 久久久久久久久大av| 七月丁香在线播放| 成人国产麻豆网| 一级毛片 在线播放| 七月丁香在线播放| 九九久久精品国产亚洲av麻豆| 亚洲美女视频黄频| 欧美日韩视频精品一区| 哪个播放器可以免费观看大片| 精品人妻熟女毛片av久久网站| 建设人人有责人人尽责人人享有的| 日韩免费高清中文字幕av| 精品久久久精品久久久| 在线 av 中文字幕| 高清欧美精品videossex| 久久国内精品自在自线图片| 成年人午夜在线观看视频| 黄色视频在线播放观看不卡| 水蜜桃什么品种好| 久久人人爽人人片av| 欧美丝袜亚洲另类| 麻豆精品久久久久久蜜桃| 免费看日本二区| 欧美3d第一页| 中文字幕人妻丝袜制服| 久久精品国产亚洲av天美| 亚洲精品乱码久久久v下载方式| 国产精品麻豆人妻色哟哟久久| 日韩精品免费视频一区二区三区 | 国产有黄有色有爽视频| 不卡视频在线观看欧美| 黑丝袜美女国产一区| 亚洲欧美一区二区三区黑人 | 永久网站在线| 交换朋友夫妻互换小说| 亚洲精品日韩在线中文字幕| av在线观看视频网站免费| 久久99精品国语久久久| 国产极品天堂在线| 蜜臀久久99精品久久宅男| 国产黄频视频在线观看| 日韩制服骚丝袜av| 各种免费的搞黄视频| 国产男人的电影天堂91| 国产精品一区二区在线不卡| 人人妻人人添人人爽欧美一区卜| 91精品国产国语对白视频| 国产精品熟女久久久久浪| 亚洲欧美一区二区三区国产| 嘟嘟电影网在线观看| 在线观看人妻少妇| 日韩精品免费视频一区二区三区 | 在线看a的网站| 永久网站在线| 男人舔奶头视频| 少妇丰满av| 亚洲综合色惰| 日韩av免费高清视频| 成人影院久久| 久久6这里有精品| 高清欧美精品videossex| 久久人人爽人人片av| 久久精品国产亚洲网站| 亚洲人成网站在线播| 久久国产精品男人的天堂亚洲 | 精品久久久久久久久av| 丰满迷人的少妇在线观看| 热re99久久国产66热| 曰老女人黄片| 国产精品久久久久久精品古装| 国产一级毛片在线| 国产91av在线免费观看| 国产日韩欧美亚洲二区| 色婷婷久久久亚洲欧美| 久久韩国三级中文字幕| 久久人人爽av亚洲精品天堂| 午夜福利影视在线免费观看| 人妻系列 视频| 久久精品国产亚洲av天美| 亚洲av成人精品一区久久| 少妇被粗大的猛进出69影院 | 精品久久久久久久久亚洲| 亚洲av男天堂| 天美传媒精品一区二区| 夜夜骑夜夜射夜夜干| 免费看av在线观看网站| 丁香六月天网| 国产伦理片在线播放av一区| 国产有黄有色有爽视频| 亚洲激情五月婷婷啪啪| 亚洲伊人久久精品综合| 国产日韩一区二区三区精品不卡 | 亚洲国产av新网站| 精品一区在线观看国产| 好男人视频免费观看在线| 秋霞在线观看毛片| 自线自在国产av| 国产午夜精品一二区理论片| 欧美激情极品国产一区二区三区 | 韩国av在线不卡| 亚洲精品乱码久久久v下载方式| 男女无遮挡免费网站观看| a级一级毛片免费在线观看| 亚洲国产日韩一区二区| 精品熟女少妇av免费看| 日韩中字成人| 久久久久视频综合| 国产日韩一区二区三区精品不卡 | 色吧在线观看| 国产亚洲最大av| 嫩草影院入口| 久久久久网色| 亚洲精品乱码久久久v下载方式| 国产欧美日韩精品一区二区| 免费大片黄手机在线观看| 日本免费在线观看一区| 久久久久久久久久久久大奶| 美女主播在线视频| 亚洲av成人精品一二三区| 2018国产大陆天天弄谢| 精品国产国语对白av| 嘟嘟电影网在线观看| 亚洲图色成人| 亚洲人成网站在线观看播放| 亚洲经典国产精华液单| 亚洲国产精品一区二区三区在线| 99热网站在线观看| 国产高清有码在线观看视频| 国产精品国产三级国产专区5o| 少妇被粗大的猛进出69影院 | 午夜福利影视在线免费观看| 又爽又黄a免费视频| 国产白丝娇喘喷水9色精品| 69精品国产乱码久久久| 色视频在线一区二区三区| 免费看av在线观看网站| 国产在线视频一区二区| 国产午夜精品久久久久久一区二区三区| 七月丁香在线播放| 国产成人免费观看mmmm| 色94色欧美一区二区| 熟女电影av网| 精品久久久久久久久av| 少妇熟女欧美另类| 国产伦在线观看视频一区| 亚洲av成人精品一二三区| 久久国内精品自在自线图片| 精品久久久噜噜| 成人毛片a级毛片在线播放| 午夜激情久久久久久久| 男人爽女人下面视频在线观看| 精品久久久噜噜| 国产极品天堂在线| 国产黄片视频在线免费观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲成人手机| 蜜桃久久精品国产亚洲av| 欧美日本中文国产一区发布| 3wmmmm亚洲av在线观看| 国产免费一区二区三区四区乱码| 嘟嘟电影网在线观看| 黄色日韩在线| 国产亚洲91精品色在线| 人妻人人澡人人爽人人| 国产日韩欧美在线精品| 国产免费福利视频在线观看| 有码 亚洲区| 在线观看www视频免费| 日韩成人伦理影院| 亚洲精品久久久久久婷婷小说| 国产精品欧美亚洲77777| 欧美xxxx性猛交bbbb| 夜夜爽夜夜爽视频| 精品熟女少妇av免费看| 在线观看三级黄色| 高清欧美精品videossex| 男人添女人高潮全过程视频| 国产成人精品无人区| 午夜激情久久久久久久| 精品一区二区三区视频在线| 日韩电影二区| 久久热精品热| 中国国产av一级| 韩国高清视频一区二区三区| 日日摸夜夜添夜夜添av毛片| 欧美亚洲 丝袜 人妻 在线| 久久久久久久精品精品| 久久精品国产a三级三级三级| 99精国产麻豆久久婷婷| 一边亲一边摸免费视频| 久久精品国产鲁丝片午夜精品| 日日啪夜夜撸| 最近手机中文字幕大全| 亚洲中文av在线| 91久久精品国产一区二区三区| 又爽又黄a免费视频| 久久99蜜桃精品久久| 中文字幕精品免费在线观看视频 | 欧美一级a爱片免费观看看| 亚洲av不卡在线观看| 国产成人精品福利久久| 中国国产av一级| 成人美女网站在线观看视频| 国产精品秋霞免费鲁丝片| 51国产日韩欧美| 国产精品一区二区在线观看99| 久久ye,这里只有精品| 免费黄频网站在线观看国产| 亚洲av成人精品一二三区| 人妻系列 视频| 亚洲欧美一区二区三区黑人 | 永久免费av网站大全| 美女福利国产在线| 成年人免费黄色播放视频 | av在线老鸭窝| 国内揄拍国产精品人妻在线| 国产无遮挡羞羞视频在线观看| 大话2 男鬼变身卡| 最近手机中文字幕大全| 国产高清不卡午夜福利| 国精品久久久久久国模美| a级一级毛片免费在线观看| 亚洲国产精品专区欧美| 国产精品久久久久久精品古装| 中文天堂在线官网| 高清欧美精品videossex| 欧美激情极品国产一区二区三区 | 一级片'在线观看视频| 久久精品国产鲁丝片午夜精品| 亚洲美女黄色视频免费看| 国产 精品1| 亚洲成人av在线免费| 偷拍熟女少妇极品色| 久久国产精品男人的天堂亚洲 | 一区在线观看完整版| 51国产日韩欧美| 精品熟女少妇av免费看| 国产在视频线精品| 久久久久久久久久人人人人人人| 啦啦啦视频在线资源免费观看| 亚洲欧美成人精品一区二区| 一级a做视频免费观看| 国产爽快片一区二区三区| 亚洲精品久久久久久婷婷小说| 五月天丁香电影| 亚洲经典国产精华液单| 美女福利国产在线| 女人久久www免费人成看片| 偷拍熟女少妇极品色| 91在线精品国自产拍蜜月| 自线自在国产av| 三级国产精品欧美在线观看| 高清在线视频一区二区三区| kizo精华| 亚洲av中文av极速乱| 人人澡人人妻人| 欧美日韩亚洲高清精品| 九九在线视频观看精品| 夫妻午夜视频| 美女视频免费永久观看网站| 国产有黄有色有爽视频| 成人影院久久| 日日啪夜夜爽| 一级,二级,三级黄色视频| av国产精品久久久久影院| 国产高清三级在线| 久久免费观看电影| 久久久久久久精品精品| 涩涩av久久男人的天堂| 精品人妻偷拍中文字幕| 国产高清有码在线观看视频| 免费少妇av软件| 日本与韩国留学比较| 嫩草影院入口| 少妇的逼水好多| 日韩不卡一区二区三区视频在线| av不卡在线播放| 亚洲性久久影院| 各种免费的搞黄视频| 日本黄色片子视频| 亚洲一区二区三区欧美精品| 在线免费观看不下载黄p国产| 色网站视频免费| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲色图综合在线观看| 亚洲精品第二区| 大又大粗又爽又黄少妇毛片口| 国产免费一级a男人的天堂| 97在线视频观看| 国产无遮挡羞羞视频在线观看| 成人免费观看视频高清| 亚洲精品乱久久久久久| 国产中年淑女户外野战色| www.av在线官网国产| 午夜激情久久久久久久| 三上悠亚av全集在线观看 | 2021少妇久久久久久久久久久| 久久婷婷青草| 少妇人妻一区二区三区视频| 老司机亚洲免费影院| 亚洲性久久影院| 最近最新中文字幕免费大全7| 韩国av在线不卡| 最近中文字幕2019免费版| a级毛片免费高清观看在线播放| 内地一区二区视频在线| 亚洲精品乱码久久久久久按摩| 国产免费视频播放在线视频| 国产精品麻豆人妻色哟哟久久| 国产精品无大码| 亚洲欧美中文字幕日韩二区| 三级经典国产精品| 欧美另类一区| 国产91av在线免费观看| 欧美日韩精品成人综合77777| 亚洲精品456在线播放app| 自拍偷自拍亚洲精品老妇| 少妇 在线观看| 在线观看av片永久免费下载| 国产成人aa在线观看| 欧美精品人与动牲交sv欧美| 黄色视频在线播放观看不卡| 午夜日本视频在线| 欧美3d第一页| 3wmmmm亚洲av在线观看| 亚洲av在线观看美女高潮| 草草在线视频免费看| 国产在线一区二区三区精| 亚洲伊人久久精品综合| 亚洲av成人精品一区久久| 国产一区亚洲一区在线观看| h视频一区二区三区| √禁漫天堂资源中文www| 国产av一区二区精品久久| 激情五月婷婷亚洲| 色5月婷婷丁香| 久久综合国产亚洲精品| 久久久久久久久久成人| 久久国产亚洲av麻豆专区| www.av在线官网国产| 91精品国产国语对白视频| 色94色欧美一区二区| 久久午夜福利片| 国产av码专区亚洲av| 美女cb高潮喷水在线观看| 午夜福利在线观看免费完整高清在| 99久久精品国产国产毛片| 亚洲欧美一区二区三区黑人 | 五月伊人婷婷丁香| 亚洲图色成人| 国产一区亚洲一区在线观看| 80岁老熟妇乱子伦牲交| 欧美日韩视频高清一区二区三区二| 欧美+日韩+精品| 免费观看性生交大片5| 黄色配什么色好看| 18禁在线无遮挡免费观看视频| 精品国产一区二区三区久久久樱花| 婷婷色综合www| 亚洲国产精品一区二区三区在线| 97精品久久久久久久久久精品| 国产精品国产三级专区第一集| 国产视频内射| 精品久久久久久久久亚洲| 亚洲不卡免费看| 国产成人免费观看mmmm| 精品视频人人做人人爽| 久久午夜综合久久蜜桃| 九草在线视频观看| 亚洲综合精品二区| 伊人亚洲综合成人网| 欧美xxxx性猛交bbbb| 自拍偷自拍亚洲精品老妇| 亚洲国产欧美日韩在线播放 | 国产色婷婷99| av天堂久久9| 一个人看视频在线观看www免费| 国产高清有码在线观看视频| 国产av精品麻豆| 欧美 亚洲 国产 日韩一| 久久99一区二区三区| 精品人妻熟女av久视频| 亚洲欧美精品专区久久| 99热这里只有精品一区| 成人国产av品久久久| 国产在线视频一区二区| av国产久精品久网站免费入址| av.在线天堂| 亚洲国产精品成人久久小说| 国产精品国产三级专区第一集| 免费看日本二区| 亚洲国产av新网站| 欧美97在线视频| av天堂久久9| 色视频www国产| 日本午夜av视频| 精品99又大又爽又粗少妇毛片| 亚洲精品日韩在线中文字幕| 欧美 日韩 精品 国产| 国产精品无大码| 亚洲人成网站在线播| 十分钟在线观看高清视频www | 少妇人妻 视频| 精品亚洲成a人片在线观看| 国产精品人妻久久久影院| 久久婷婷青草| 午夜免费鲁丝| 精品国产露脸久久av麻豆| 在现免费观看毛片| 免费人妻精品一区二区三区视频| 自线自在国产av| 国产一区亚洲一区在线观看| 在线播放无遮挡| 亚洲,一卡二卡三卡| 国产av码专区亚洲av| 国产精品一区www在线观看| 自拍偷自拍亚洲精品老妇| 搡老乐熟女国产| 最近2019中文字幕mv第一页| 男女免费视频国产| 精品少妇久久久久久888优播| 国产精品国产av在线观看| 国产精品人妻久久久影院| 久久久国产欧美日韩av| 亚洲国产精品专区欧美| 在线观看一区二区三区激情| 国产视频首页在线观看| 黄色配什么色好看| 热99国产精品久久久久久7| 亚洲国产精品一区二区三区在线| 一本大道久久a久久精品| 你懂的网址亚洲精品在线观看| 一区二区三区乱码不卡18| 精品国产露脸久久av麻豆| 中文字幕精品免费在线观看视频 | 狂野欧美激情性bbbbbb| 纵有疾风起免费观看全集完整版| 亚洲婷婷狠狠爱综合网| 九草在线视频观看| 免费观看性生交大片5| 高清视频免费观看一区二区| 国产日韩一区二区三区精品不卡 | 伊人久久国产一区二区| 亚洲欧洲国产日韩| 中文字幕制服av| 日本午夜av视频| 国产精品久久久久久精品古装| 2022亚洲国产成人精品| 51国产日韩欧美| 久久免费观看电影| 人妻少妇偷人精品九色| 在线看a的网站| 亚洲四区av| 人人妻人人澡人人爽人人夜夜| 久久6这里有精品| 黄色配什么色好看| 黄色毛片三级朝国网站 | 一二三四中文在线观看免费高清| 伦精品一区二区三区| 大片电影免费在线观看免费| 国产精品一区二区在线不卡| 精品一品国产午夜福利视频| 亚洲国产欧美在线一区| 国产高清三级在线| 亚洲国产欧美在线一区| 大又大粗又爽又黄少妇毛片口| 色吧在线观看| 亚洲精品乱码久久久v下载方式| 男女啪啪激烈高潮av片| 日韩亚洲欧美综合| 熟女电影av网| 午夜老司机福利剧场| 黑丝袜美女国产一区| videossex国产| av有码第一页| 男男h啪啪无遮挡| 成人免费观看视频高清| 在线观看免费高清a一片| 久久久久国产网址| 日韩成人av中文字幕在线观看| 内射极品少妇av片p| 免费大片18禁| 国产午夜精品一二区理论片| 91久久精品国产一区二区三区| 蜜桃在线观看..| 国产高清国产精品国产三级| 91精品一卡2卡3卡4卡| 我要看日韩黄色一级片| 国产午夜精品一二区理论片| 成人特级av手机在线观看| 伊人久久国产一区二区| 免费人妻精品一区二区三区视频| 国产日韩欧美视频二区| 中文天堂在线官网| 婷婷色综合www| 国产高清不卡午夜福利| 2021少妇久久久久久久久久久| 中文欧美无线码| 在线观看人妻少妇| 青青草视频在线视频观看| 高清午夜精品一区二区三区| 国产亚洲一区二区精品| 丰满人妻一区二区三区视频av| 日韩成人伦理影院| 免费观看的影片在线观看| 国产精品福利在线免费观看| 亚洲av成人精品一区久久| 久久99热这里只频精品6学生| 日本黄大片高清| 国产精品久久久久久久久免| 国产男女内射视频| h日本视频在线播放| 久久久久国产精品人妻一区二区| 少妇人妻久久综合中文| 午夜福利,免费看| 国产精品熟女久久久久浪| 婷婷色综合大香蕉| 日韩欧美一区视频在线观看 | 久久精品国产亚洲av涩爱| 国产老妇伦熟女老妇高清| 国产欧美日韩综合在线一区二区 | 一本色道久久久久久精品综合| 人人妻人人爽人人添夜夜欢视频 | 成人亚洲欧美一区二区av| a级毛色黄片| 一本色道久久久久久精品综合| 十八禁网站网址无遮挡 | 五月伊人婷婷丁香| 日韩欧美一区视频在线观看 | 男女边摸边吃奶| 又爽又黄a免费视频| 亚洲情色 制服丝袜| 三级国产精品欧美在线观看| 熟妇人妻不卡中文字幕| 我要看日韩黄色一级片| 永久网站在线| 久久久久久人妻| av一本久久久久| 夜夜看夜夜爽夜夜摸| 亚洲精品456在线播放app| 午夜久久久在线观看| 国产高清不卡午夜福利| 婷婷色av中文字幕| 一边亲一边摸免费视频| 久久国产精品大桥未久av | 色吧在线观看| 久久鲁丝午夜福利片| 久久精品夜色国产| 99热网站在线观看| 国产探花极品一区二区| 下体分泌物呈黄色| 国产成人91sexporn| 国产毛片在线视频| 亚洲欧美中文字幕日韩二区| 久久久久久伊人网av| 亚洲欧美日韩东京热| 精品人妻熟女av久视频| 女人精品久久久久毛片| 看免费成人av毛片| 久久人人爽人人片av| 亚洲国产欧美在线一区| 人人妻人人澡人人爽人人夜夜| 性色av一级| 国产精品三级大全| 全区人妻精品视频| 亚洲精品国产av蜜桃| 亚洲精品国产色婷婷电影| 天堂中文最新版在线下载| 99热这里只有精品一区| 成人漫画全彩无遮挡| 交换朋友夫妻互换小说| 乱码一卡2卡4卡精品| 最近中文字幕高清免费大全6|