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

    計(jì)及風(fēng)電—負(fù)荷耦合關(guān)系的含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化

    2017-12-22 10:13:45吳耀武汪昌霜婁素華劉寶林王曦冉
    電力系統(tǒng)自動(dòng)化 2017年21期
    關(guān)鍵詞:調(diào)峰出力風(fēng)電

    吳耀武, 汪昌霜, 婁素華, 劉寶林, 王曦冉

    (1. 強(qiáng)電磁工程與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(華中科技大學(xué)), 湖北省武漢市 430074;2. 華中科技大學(xué)電氣與電子工程學(xué)院, 湖北省武漢市 430074; 3. 云南電網(wǎng)有限責(zé)任公司, 云南省昆明市 650011)

    計(jì)及風(fēng)電—負(fù)荷耦合關(guān)系的含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化

    吳耀武1,2, 汪昌霜1,2, 婁素華1,2, 劉寶林3, 王曦冉1,2

    (1. 強(qiáng)電磁工程與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(華中科技大學(xué)), 湖北省武漢市 430074;2. 華中科技大學(xué)電氣與電子工程學(xué)院, 湖北省武漢市 430074; 3. 云南電網(wǎng)有限責(zé)任公司, 云南省昆明市 650011)

    隨機(jī)性風(fēng)電出力的準(zhǔn)確建模是含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化的關(guān)鍵環(huán)節(jié)之一。針對(duì)中國(guó)系統(tǒng)運(yùn)行的實(shí)際特點(diǎn),提出計(jì)及風(fēng)電出力與日負(fù)荷時(shí)序耦合關(guān)系的交互特性評(píng)價(jià)指標(biāo),并據(jù)此建立基于調(diào)峰問題驅(qū)動(dòng)的風(fēng)電出力模型。該模型以關(guān)鍵場(chǎng)景、聚類場(chǎng)景及其概率分布表征風(fēng)電接入對(duì)系統(tǒng)運(yùn)行和調(diào)峰需求的綜合影響。基于該風(fēng)電出力模型,建立了含大規(guī)模風(fēng)電系統(tǒng)的調(diào)峰運(yùn)行優(yōu)化模型。通過多場(chǎng)景下的運(yùn)行優(yōu)化,獲得各項(xiàng)運(yùn)行指標(biāo)的期望值,實(shí)現(xiàn)對(duì)系統(tǒng)調(diào)峰運(yùn)行特性的綜合分析。對(duì)中國(guó)西部某省進(jìn)行了算例分析,驗(yàn)證了所提方法的有效性和實(shí)用性。

    耦合關(guān)系; 關(guān)鍵場(chǎng)景; 聚類場(chǎng)景; 調(diào)峰; 運(yùn)行優(yōu)化

    0 引言

    中國(guó)作為一個(gè)貧油、少氣、相對(duì)富煤的國(guó)家[1],燃煤發(fā)電在未來相當(dāng)長(zhǎng)的一段時(shí)期仍將是多種能源發(fā)電系統(tǒng)的主體電源之一。煤電機(jī)組由于動(dòng)力設(shè)備繁多,控制操作復(fù)雜,其發(fā)電出力調(diào)節(jié)范圍比較有限,調(diào)節(jié)速度相對(duì)較慢,且頻繁啟停影響運(yùn)行經(jīng)濟(jì)性,因此系統(tǒng)的調(diào)峰能力受到極大限制。隨著風(fēng)電并網(wǎng)規(guī)模不斷增大,其固有的出力不確定性、反調(diào)峰特性給中國(guó)電力系統(tǒng)調(diào)峰帶來了更大挑戰(zhàn)[2-4]。國(guó)家能源局統(tǒng)計(jì)數(shù)據(jù)顯示,截至2016年底,中國(guó)風(fēng)電累計(jì)并網(wǎng)容量為149 GW,全年累計(jì)棄風(fēng)電量約為4.97×104GW·h,全國(guó)平均棄風(fēng)率達(dá)17%[5]??梢?受限于現(xiàn)有電力系統(tǒng)的調(diào)峰能力[6],中國(guó)多數(shù)以火電為主體電源的省級(jí)系統(tǒng)對(duì)大容量風(fēng)電的有效接納仍是個(gè)難題。因此,研究含大規(guī)模風(fēng)電系統(tǒng)的調(diào)峰運(yùn)行特性,客觀評(píng)價(jià)風(fēng)電并網(wǎng)系統(tǒng)的風(fēng)電消納能力,對(duì)于促進(jìn)可再生能源發(fā)電產(chǎn)業(yè)持續(xù)健康發(fā)展具有重要意義。

    隨機(jī)性風(fēng)電出力的準(zhǔn)確建模是電力系統(tǒng)調(diào)峰運(yùn)行特性研究的關(guān)鍵環(huán)節(jié)。文獻(xiàn)[7-8]基于日前風(fēng)功率預(yù)測(cè),以預(yù)測(cè)誤差反映風(fēng)電出力的不確定性,研究風(fēng)電對(duì)日前機(jī)組組合的影響,分析系統(tǒng)次日的調(diào)峰能力。文獻(xiàn)[9-10]通過風(fēng)速或風(fēng)功率數(shù)據(jù)假設(shè)其服從的概率分布或隨機(jī)過程,模擬風(fēng)電中長(zhǎng)期時(shí)序出力,并據(jù)此提出風(fēng)電并網(wǎng)系統(tǒng)的調(diào)峰裕度評(píng)估方法,但現(xiàn)有模擬方法難以準(zhǔn)確表征風(fēng)電出力的時(shí)序相關(guān)性。文獻(xiàn)[11]考慮對(duì)調(diào)峰影響較大的風(fēng)電出力極端場(chǎng)景,提出了一種負(fù)荷低谷時(shí)段調(diào)峰能力計(jì)算方法,計(jì)算結(jié)果較為保守。文獻(xiàn)[12]基于風(fēng)電加入后的凈負(fù)荷曲線的調(diào)峰需求特征,應(yīng)用場(chǎng)景削減技術(shù)獲得三類風(fēng)電出力典型場(chǎng)景對(duì)風(fēng)電隨機(jī)性建模,但由于場(chǎng)景數(shù)量較少,還不能全面反映風(fēng)電對(duì)系統(tǒng)調(diào)峰運(yùn)行優(yōu)化的綜合影響。電力系統(tǒng)調(diào)峰問題本身就是負(fù)荷與電源之間的動(dòng)態(tài)平衡導(dǎo)致的[13],因而風(fēng)電和日負(fù)荷間的時(shí)序耦合關(guān)系是影響系統(tǒng)調(diào)峰平衡能力的主要因素,目前的研究文獻(xiàn)在風(fēng)電出力建模中均未能深入地考慮這種耦合關(guān)系。

    針對(duì)現(xiàn)有研究的不足,考慮煤電比例較高系統(tǒng)的發(fā)電出力及調(diào)節(jié)特性,在充分計(jì)及風(fēng)電與日負(fù)荷間的時(shí)序耦合關(guān)系基礎(chǔ)上,本文提出了一種基于調(diào)峰問題驅(qū)動(dòng)的風(fēng)電出力模型。以此為基礎(chǔ),提出以系統(tǒng)總計(jì)算費(fèi)用最小為目標(biāo)的含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化模型,模型中不僅考慮了風(fēng)電帶來的能源節(jié)約效益和環(huán)境效益,也考慮了其出力變化特征對(duì)系統(tǒng)供電可靠性的影響。對(duì)中國(guó)西部某省進(jìn)行了算例分析,驗(yàn)證了本文所提方法的有效性和實(shí)用性。

    1 基于調(diào)峰問題驅(qū)動(dòng)的風(fēng)電出力模型

    為適應(yīng)含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化的需要,本文以計(jì)及風(fēng)電—負(fù)荷耦合關(guān)系的交互特性評(píng)價(jià)指標(biāo)為基礎(chǔ),建立基于調(diào)峰問題驅(qū)動(dòng)的風(fēng)電出力模型。

    1.1 計(jì)及風(fēng)電—負(fù)荷耦合關(guān)系的交互特性評(píng)價(jià)指標(biāo)

    在運(yùn)行優(yōu)化中,為系統(tǒng)、定量地表征風(fēng)電和負(fù)荷間的時(shí)序耦合關(guān)系,須根據(jù)負(fù)荷的變化特征將日負(fù)荷曲線劃分為高峰時(shí)段TH、低谷時(shí)段TL、腰荷時(shí)段TM,即

    TH={t|L(t)≥ρpLmax}

    (1)

    TL={t|L(t)≤ρbLmin}

    (2)

    TM={t|ρbLmin

    (3)

    式中:L(t)表示日負(fù)荷曲線t時(shí)刻負(fù)荷;Lmax為日最大負(fù)荷;Lmin為日最小負(fù)荷;ρp<1,ρb>1。

    根據(jù)風(fēng)電在日負(fù)荷不同時(shí)段的出力特性,考慮大規(guī)模風(fēng)電接入對(duì)系統(tǒng)調(diào)峰運(yùn)行優(yōu)化的影響,本文從三個(gè)方面建立計(jì)及風(fēng)電—負(fù)荷耦合關(guān)系的交互特性評(píng)價(jià)指標(biāo)。

    (4)

    (5)

    (6)

    式中:max(·)為最大值函數(shù);min(·)為最小值函數(shù);PW(dt)為研究周期內(nèi)第d日風(fēng)電出力曲線PW(d)在t時(shí)刻的出力;N為研究周期內(nèi)風(fēng)電出力曲線數(shù)量;Nsize(TH)為高峰時(shí)段TH包含的時(shí)刻數(shù)。

    (7)

    (8)

    (9)

    (10)

    (11)

    這三類指標(biāo)反映了風(fēng)電出力與負(fù)荷變化的耦合關(guān)系,據(jù)此可衍生出:風(fēng)電出力最大調(diào)峰需求CWmax(d)、最小調(diào)峰需求CWmin(d),二者共同表征了風(fēng)電出力的調(diào)峰特性,CWmax(d)值越大且CWmin(d)不小于0,表明風(fēng)電的反調(diào)峰特性越強(qiáng)。

    (12)

    (13)

    結(jié)合風(fēng)電自身的出力特征和上述交互特性評(píng)價(jià)指標(biāo),可建立基于調(diào)峰問題驅(qū)動(dòng)的風(fēng)電出力模型。該模型以關(guān)鍵場(chǎng)景、聚類場(chǎng)景及其概率分布表征隨機(jī)性風(fēng)電接入對(duì)系統(tǒng)運(yùn)行和調(diào)峰需求的綜合影響。

    1.2 基于調(diào)峰問題驅(qū)動(dòng)的風(fēng)電出力關(guān)鍵場(chǎng)景

    風(fēng)電出力關(guān)鍵場(chǎng)景是指出現(xiàn)概率較小,但對(duì)系統(tǒng)調(diào)峰運(yùn)行優(yōu)化有重大影響的場(chǎng)景,此類場(chǎng)景需予以單獨(dú)考慮。考慮風(fēng)電出力對(duì)常規(guī)電源機(jī)組組合和系統(tǒng)調(diào)峰容量需求的顯著影響,本文將優(yōu)選出兩個(gè)關(guān)鍵場(chǎng)景:①電力平衡場(chǎng)景;②調(diào)峰平衡場(chǎng)景。

    1.2.1風(fēng)電出力電力平衡場(chǎng)景

    由于風(fēng)電出力的不確定性,為保證系統(tǒng)發(fā)電充裕度,制定合理的開停機(jī)計(jì)劃,須按照負(fù)荷高峰時(shí)段風(fēng)電出力較小的情況確定電力平衡場(chǎng)景,方法如下。

    α∈[0,1)

    (14)

    (15)

    PWΩ1′={PW(D1′)|p(PWave≤PWave(D1′))>

    β&PW(D1′)∈PWΩ1}β∈[0,1)

    (16)

    3)確定研究周期內(nèi)置信水平α下的風(fēng)電出力電力平衡場(chǎng)景PW1,計(jì)算其概率pW1。

    (17)

    (18)

    1.2.2風(fēng)電出力調(diào)峰平衡場(chǎng)景

    為反映風(fēng)電出力對(duì)系統(tǒng)調(diào)峰平衡的顯著影響,須按照風(fēng)電反調(diào)峰特性較強(qiáng),即風(fēng)電出力調(diào)峰需求較大的情況確定調(diào)峰平衡場(chǎng)景,方法如下。

    1)給定置信水平γ,篩選出風(fēng)電曲線PW(D2),使以超過γ的概率確信研究周期內(nèi)風(fēng)電出力最大調(diào)峰需求不大于CWmax(D2),PW(D2)構(gòu)成的集合記為PWΩ2,其中,最小調(diào)峰需求特性指標(biāo)最大的風(fēng)電曲線記為Pμ2。

    PWΩ2={PW(D2)|p(CWmax≤CWmax(D2))>γ}

    γ∈[0,1)

    (19)

    2)為充分保證調(diào)峰平衡場(chǎng)景的反調(diào)峰特性,從PWΩ2中篩選出最小調(diào)峰需求特性指標(biāo)不小于0的風(fēng)電曲線集合PWΩ2′。PWΩ2′中最大調(diào)峰需求CWmax最小的風(fēng)電曲線記為Pν2,該曲線為CWmax指標(biāo)最接近置信水平γ的場(chǎng)景。

    PWΩ2′={PW(D2′)|CWmin(D2′)≥

    0&PW(D2′)∈PWΩ2}

    (20)

    3)確定研究周期內(nèi)置信水平γ下的風(fēng)電出力調(diào)峰平衡場(chǎng)景PW2,計(jì)算其概率pW2。

    (21)

    (22)

    在電力系統(tǒng)調(diào)峰運(yùn)行優(yōu)化中,采用關(guān)鍵場(chǎng)景中的電力平衡場(chǎng)景確定常規(guī)電源開機(jī)組合,能夠充分保證供電可靠性;采用調(diào)峰平衡場(chǎng)景進(jìn)行系統(tǒng)的運(yùn)行模擬可客觀反映系統(tǒng)的調(diào)峰容量需求。

    1.3 基于調(diào)峰問題驅(qū)動(dòng)及改進(jìn)K-means聚類算法的風(fēng)電出力聚類場(chǎng)景

    K-means聚類算法的基本思想是:對(duì)于給定的包含n個(gè)m維數(shù)據(jù)點(diǎn)的數(shù)據(jù)集X={x1,x2,…,xd,…,xn},其中xd為m維數(shù)組,將其劃分為K個(gè)數(shù)據(jù)類Ck,每個(gè)數(shù)據(jù)類Ck對(duì)應(yīng)一個(gè)聚類中心Uk,計(jì)算類中各數(shù)據(jù)樣本到Uk的距離,通過反復(fù)迭代使得所有數(shù)據(jù)類的總距離平方和最小[14-15]。

    本文針對(duì)調(diào)峰問題對(duì)K-means聚類算法適當(dāng)改進(jìn),生成風(fēng)電出力聚類場(chǎng)景,以較高的精度和較低的計(jì)算耗費(fèi)表征原始風(fēng)電曲線對(duì)系統(tǒng)調(diào)峰運(yùn)行的影響。

    1.3.1距離準(zhǔn)則函數(shù)建立方法

    常規(guī)聚類一般僅根據(jù)風(fēng)電電源自身出力特征建立距離準(zhǔn)則函數(shù),應(yīng)用于調(diào)峰運(yùn)行優(yōu)化不夠準(zhǔn)確。本文根據(jù)調(diào)峰運(yùn)行優(yōu)化的特點(diǎn),采用計(jì)及風(fēng)電—負(fù)荷時(shí)序耦合關(guān)系的交互特性評(píng)價(jià)指標(biāo),建立距離準(zhǔn)則函數(shù),這樣可以確保將各負(fù)荷時(shí)段出力水平最相似的風(fēng)電曲線聚為一類,從而保證聚類場(chǎng)景應(yīng)用于調(diào)峰運(yùn)行優(yōu)化的適應(yīng)性,具體方法流程如下。

    (23)

    式中:d=1,2,…,N-2;k=1,2,…,K;λθ為權(quán)重,根據(jù)不同負(fù)荷時(shí)段風(fēng)電出力對(duì)運(yùn)行優(yōu)化的重要性賦予不同的權(quán)重,使聚類結(jié)果具有一定傾向性,更好地適應(yīng)調(diào)峰問題研究,一般取λTH=λTL>λTM。

    3)將風(fēng)電曲線PW(d)分配給使得S(xd,Uk)值最小的類Ck,更新Uk,計(jì)算距離準(zhǔn)則函數(shù)J:

    (24)

    (25)

    不斷迭代至J收斂,即可確定適應(yīng)調(diào)峰運(yùn)行優(yōu)化的風(fēng)電出力聚類結(jié)果。

    1.3.2確定風(fēng)電出力聚類場(chǎng)景

    (26)

    c(b,b′)=

    (27)

    式中:Ck′為削減后的目標(biāo)場(chǎng)景集合;LKD(Ck,Ck′)為Ck和Ck′的LKD;b和b′分別表示Ck和Ck′中的風(fēng)電場(chǎng)景;pb為風(fēng)電場(chǎng)景b在Ck中的概率;λσ為權(quán)重,一般取λσ≥0.6。

    (28)

    各聚類場(chǎng)景的概率為:

    (29)

    以關(guān)鍵場(chǎng)景和聚類場(chǎng)景構(gòu)成的風(fēng)電出力典型場(chǎng)景集S對(duì)風(fēng)電出力隨機(jī)性建模,可以表征風(fēng)電接入對(duì)系統(tǒng)運(yùn)行和調(diào)峰需求的綜合影響。

    2 含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化模型

    2.1 目標(biāo)函數(shù)

    為確保電力系統(tǒng)安全和可靠供電,本文首先采用風(fēng)電出力電力平衡場(chǎng)景確定常規(guī)電源開機(jī)組合,然后針對(duì)同一開機(jī)方式,優(yōu)化各風(fēng)電場(chǎng)景下各類型機(jī)組的工作位置和發(fā)電出力,使電力系統(tǒng)總計(jì)算費(fèi)用最小。含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化的目標(biāo)函數(shù)如下:

    minF=

    (30)

    (31)

    (32)

    (33)

    (34)

    2.2 約束條件

    含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化模型的約束條件包括:系統(tǒng)電力平衡約束、系統(tǒng)旋轉(zhuǎn)/事故備用約束、機(jī)組運(yùn)行狀態(tài)邏輯約束、機(jī)組出力約束、火電最小開停機(jī)時(shí)間約束、火電保安開機(jī)約束、水電電量平衡約束等,其表達(dá)式限于篇幅不再贅述。

    2.3 調(diào)峰棄電協(xié)調(diào)策略

    在庫(kù)容式水電站及常規(guī)火電的調(diào)峰能力已全部利用后,若系統(tǒng)出現(xiàn)調(diào)峰不足,考慮到大規(guī)模風(fēng)力發(fā)電的可控性較差,而水電站來水在短期內(nèi)比較穩(wěn)定,為保證系統(tǒng)的調(diào)峰平衡及供電質(zhì)量,優(yōu)先限制風(fēng)電出力,其次限制水電站出力。

    3 算例分析

    為驗(yàn)證本文所提模型,以中國(guó)西部某省2020年電源規(guī)劃為基礎(chǔ),不考慮光伏和生物質(zhì)發(fā)電(比例很小),進(jìn)行系統(tǒng)調(diào)峰能力評(píng)估。該系統(tǒng)2020年火電裝機(jī)容量為23 100 MW,水電裝機(jī)容量為12 780 MW,風(fēng)電裝機(jī)容量為11 000 MW,其中,水電站共18個(gè),含徑流式水電裝機(jī)容量為8 270 MW,庫(kù)容式水電裝機(jī)容量為4 510 MW;火電機(jī)組技術(shù)特性見表1;風(fēng)電特性采用該省某地區(qū)2014年未棄風(fēng)的標(biāo)幺出力歷史數(shù)據(jù),可利用小時(shí)數(shù)約為1 875 h;系統(tǒng)年最大負(fù)荷為28 150 MW。

    表1 火電機(jī)組技術(shù)特性Table 1 Technical characteristics of thermal power unit

    其他參數(shù)取值:標(biāo)煤?jiǎn)蝺r(jià)取600元/t,排污費(fèi)標(biāo)準(zhǔn)取0.6元/污染當(dāng)量,失負(fù)荷懲罰系數(shù)yB取8元/(kW·h),棄風(fēng)懲罰系數(shù)yW取0.5元/(kW·h),碳排放懲罰系數(shù)yCO2取120元/t,系統(tǒng)允許棄風(fēng)率ξ取10%,碳排放配額MΔ取6×107t。

    3.1 風(fēng)電出力典型場(chǎng)景集有效性分析

    由于該地區(qū)風(fēng)電具有較顯著的季節(jié)特性,本文對(duì)風(fēng)電場(chǎng)景的研究周期取為季度,即分季節(jié)生成風(fēng)電出力典型場(chǎng)景集。以春季為例,原始場(chǎng)景含92條風(fēng)電出力曲線,采用本文方法可以獲得2個(gè)風(fēng)電出力關(guān)鍵場(chǎng)景和9個(gè)風(fēng)電出力聚類場(chǎng)景,分別如圖1和圖2所示。圖1中2個(gè)關(guān)鍵場(chǎng)景的概率均為0.011,圖2中各聚類場(chǎng)景的概率標(biāo)于圖例中。

    圖1 風(fēng)電出力關(guān)鍵場(chǎng)景Fig.1 Key scenarios of wind power output

    圖2 風(fēng)電出力聚類場(chǎng)景Fig.2 Clustering scenarios of wind power output

    考慮到原始風(fēng)電曲線數(shù)量龐大,參與含多類型常規(guī)電源的實(shí)際大系統(tǒng)運(yùn)行優(yōu)化非常復(fù)雜,故本節(jié)暫將水電替換為等容量火電機(jī)組,驗(yàn)證通過本文方法提取的風(fēng)電出力典型場(chǎng)景集應(yīng)用于調(diào)峰計(jì)算的有效性,典型場(chǎng)景集和原始場(chǎng)景的計(jì)算結(jié)果對(duì)比見表2,關(guān)鍵場(chǎng)景對(duì)運(yùn)行優(yōu)化的影響分析見表3和表4。

    表2 風(fēng)電典型場(chǎng)景集與原始場(chǎng)景計(jì)算結(jié)果對(duì)比Table 2 Results comparison of wind power typical scenarios set and original scenarios

    表3 電力平衡場(chǎng)景對(duì)調(diào)峰運(yùn)行優(yōu)化的影響Table 3 Utility of power balance scenario on operation optimization for peak load regulation

    表4 調(diào)峰平衡場(chǎng)景對(duì)調(diào)峰運(yùn)行優(yōu)化的影響Table 4 Utility of peak shaving balance scenario on operation optimization for peak load regulation

    從圖1可以看出,兩個(gè)關(guān)鍵場(chǎng)景與負(fù)荷的高峰、低谷出力特性變化密切相關(guān),且均為明顯的反調(diào)峰特性。從表3可以看出,運(yùn)行優(yōu)化中采用電力平衡場(chǎng)景確定開機(jī)組合,保證了系統(tǒng)的供電可靠性,同時(shí),當(dāng)風(fēng)電滲透率提高到40%時(shí),由于失負(fù)荷懲罰費(fèi)用大幅下降,系統(tǒng)總計(jì)算費(fèi)用比按風(fēng)電平均場(chǎng)景開機(jī)更低;由表4可以看出,采用含調(diào)峰平衡場(chǎng)景的調(diào)峰運(yùn)行模擬,可客觀反映風(fēng)電接入后系統(tǒng)的最大調(diào)峰容量需求。由此可見,基于關(guān)鍵場(chǎng)景的調(diào)峰運(yùn)行模擬更加符合中國(guó)電力調(diào)度的特點(diǎn),模擬計(jì)算得到的數(shù)據(jù)更接近實(shí)際運(yùn)行情況。

    從圖2可知,對(duì)原始場(chǎng)景聚類后得到的9個(gè)聚類場(chǎng)景涵蓋了原始風(fēng)電場(chǎng)景的正調(diào)峰、反調(diào)峰、平調(diào)峰等全部特征,其中反調(diào)峰場(chǎng)景概率和達(dá)0.63。由表2可知,采用本文風(fēng)電出力典型場(chǎng)景集進(jìn)行系統(tǒng)調(diào)峰分析計(jì)算具有較高的精度,其中棄電率誤差不超過0.15%,調(diào)峰深度誤差不超過0.01%。由此可見,本文提出的風(fēng)電場(chǎng)景提取方法在含大規(guī)模風(fēng)電系統(tǒng)的調(diào)峰運(yùn)行特性分析研究中是十分有效的。

    3.2 含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行特性分析

    在風(fēng)電出力典型場(chǎng)景集的基礎(chǔ)上,對(duì)該省實(shí)際系統(tǒng)2020年全年調(diào)峰運(yùn)行進(jìn)行了模擬計(jì)算,其中,水電站按各自的水文特性和調(diào)節(jié)性能參與系統(tǒng)運(yùn)行。獲得各項(xiàng)運(yùn)行特性指標(biāo)的年期望值見表5,各月風(fēng)電棄電率分布如圖3所示。

    表5 系統(tǒng)調(diào)峰運(yùn)行特性指標(biāo)年期望值Table 5 Annual expected value of operation characteristic indices for peak load regulation

    圖3 各月風(fēng)電棄電率期望值Fig.3 Expected value of wind power abandon rate in each month

    由表5可以看出,風(fēng)電接入后火電利用小時(shí)數(shù)減少783 h,系統(tǒng)總能耗、排污量、碳排放亦隨之降低,總計(jì)算費(fèi)用減少了16.7%。由于風(fēng)電出力的反調(diào)峰特性,系統(tǒng)調(diào)峰壓力增大,火電機(jī)組調(diào)峰深度由14.5%增大至32.8%,這亦導(dǎo)致火電實(shí)際單耗增大。由于系統(tǒng)調(diào)峰能力有限,風(fēng)電的大規(guī)模接入導(dǎo)致系統(tǒng)調(diào)峰不足概率增大,風(fēng)電棄電率達(dá)12.1%。

    由圖3可見,7,8月份風(fēng)電棄電率明顯高于其他月份,主要原因在于7,8月份是豐水期,來水量最大,水電發(fā)電量增大,調(diào)峰能力大幅下降,在負(fù)荷低谷時(shí)段風(fēng)電大發(fā)時(shí),不得不大量棄風(fēng)。這也從側(cè)面說明了分季節(jié)建立風(fēng)電出力典型場(chǎng)景集的必要性。

    4 結(jié)語(yǔ)

    風(fēng)電和負(fù)荷間的時(shí)序耦合關(guān)系是影響系統(tǒng)調(diào)峰平衡和風(fēng)電接納的主要因素。本文從多個(gè)角度建立了計(jì)及風(fēng)電—負(fù)荷時(shí)序耦合關(guān)系的交互特性評(píng)價(jià)指標(biāo),并據(jù)此提出了一種以關(guān)鍵場(chǎng)景、聚類場(chǎng)景及其概率分布對(duì)風(fēng)電出力建模的方法。在此基礎(chǔ)上,建立了一種考慮失負(fù)荷懲罰、棄風(fēng)懲罰、碳排放懲罰費(fèi)用的含大規(guī)模風(fēng)電系統(tǒng)調(diào)峰運(yùn)行優(yōu)化模型。算例分析表明,本文方法應(yīng)用于含風(fēng)電系統(tǒng)的調(diào)峰分析計(jì)算具有較高的精度,模擬計(jì)算結(jié)果更接近系統(tǒng)實(shí)際運(yùn)行情況。當(dāng)電力系統(tǒng)含多類型高比例可再生能源時(shí),如何對(duì)本文模型進(jìn)行改進(jìn)和擴(kuò)展,使其具有更廣泛的適應(yīng)性,有待下一步深入研究。

    [1] 劉金朋.基于資源與環(huán)境約束的中國(guó)能源供需格局發(fā)展研究[D].北京:華北電力大學(xué),2013.

    [2] XU Yanping, SUN Yong, ZHENG Taiyi, et al. Time-series characteristics of wind power and its impact on Jilin power grid[J]. Journal of Power and Energy Engineering, 2014, 2(4): 203-212.

    [3] BILLINTON R, KARKI B. Effect of hourly wind trends on the peak load-carrying capability of wind-integrated power systems[J]. Proceedings of the Institution of Mechanical Engineers, 2009, 223(4): 279-287.

    [4] 曹陽(yáng),李鵬,袁越,等.基于時(shí)序仿真的新能源消納能力分析及其低碳效益評(píng)估[J].電力系統(tǒng)自動(dòng)化,2014,38(17):60-66.DOI:10.7500/AEPS20140329002.

    CAO Yang, LI Peng, YUAN Yue, et al. Analysis on peak load regulation capability of power grid integrated with wind farms in valley load period[J]. Automation of Electric Power Systems, 2014, 38(17): 60-66. DOI: 10.7500/AEPS 20140329002.

    [5] 國(guó)家能源局.2016年風(fēng)電并網(wǎng)運(yùn)行情況[EB/OL].[2017-01-26].http://www.nea.gov.cn/2017-01/26/c_13601 4615.htm.

    [6] 魯宗相,李海波,喬穎.含高比例可再生能源電力系統(tǒng)靈活性規(guī)劃及挑戰(zhàn)[J].電力系統(tǒng)自動(dòng)化,2016,40(13):147-158.DOI:10.7500/AEPS20151215008.

    LU Zongxiang, LI Haibo, QIAO Ying. Power system flexibility planning and challenges considering high proportion of renewable energy[J]. Automation of Electric Power Systems, 2016, 40(13): 147-158. DOI: 10.7500/AEPS 20151215008.

    [7] 楊宏,苑津莎,吳立增.基于風(fēng)電功率預(yù)測(cè)的電網(wǎng)靜態(tài)調(diào)峰能力極限研究[J].電工技術(shù)學(xué)報(bào),2014,29(4):266-273.

    YANG Hong, YUAN Jinsha, WU Lizeng. Research of static regulation capacity limit based on wind power prediction[J]. Transactions of China Electrotechnical Society, 2014, 29(4): 266-273.

    [8] 賈文昭,康重慶,李丹,等.基于日前風(fēng)功率預(yù)測(cè)的風(fēng)電消納能力評(píng)估方法[J].電網(wǎng)技術(shù),2012,36(8):69-75.

    JIA Wenzhao, KANG Chongqing, LI Dan, et al. Evaluation on capability of wind power accommodation based on its day-ahead forecasting[J]. Power System Technology, 2012, 36(8): 69-75.

    [9] 張宏宇,印永華,申洪,等.基于序貫蒙特卡洛方法的風(fēng)電并網(wǎng)系統(tǒng)調(diào)峰裕度評(píng)估[J].電力系統(tǒng)自動(dòng)化,2012,36(1):32-37.

    ZHANG Hongyu, YIN Yonghua, SHEN Hong, et al. Peaking-shaving margin evaluation associated with wind power integrated system based on sequential Monte-Carlo method[J]. Automation of Electric Power Systems, 2012, 36(1): 32-37.

    [10] 張文韜,黃鑌,王秀麗,等.基于多狀態(tài)新能源發(fā)電模擬的日運(yùn)行方式及調(diào)峰安排研究[J].電網(wǎng)技術(shù),2015,39(8):2100-2106.

    ZHANG Wentao, HUANG Bin, WANG Xiuli, et al. Research on daily operation mode and peak load regulation scheduling based on generation simulation of multi-state renewable energy resources[J]. Power System Technology, 2015, 39(8): 2100-2106.

    [11] 楊冬鋒,周蘇荃,鮑鋒.風(fēng)電并網(wǎng)系統(tǒng)低谷時(shí)段的調(diào)峰能力分析[J].電網(wǎng)技術(shù),2014,38(6):1446-1451.

    YANG Dongfeng, ZHOU Suquan, BAO Feng. Analysis on peak load regulation capability of power grid integrated with wind farms in valley load period[J]. Power System Technology, 2014, 38(6): 1446-1451.

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

    HOU Tingting, LOU Suhua, WU Yaowu, et al. Analysis on peak load regulation operation characteristics of power system integrated with large-scale wind power[J]. Transactions of China Electrotechnical Society, 2013, 28(5): 105-111.

    [13] BASIT A, HANSEN A D, ALTIN M, et al. Compensating active power imbalances in power system with large-scale wind power penetration[J]. Journal of Modern Power Systems and Clean Energy, 2016, 4(2): 229-237.

    [14] GOYAL M, KUMAR S. Improving the initial centroids ofk-means clustering algorithm to generalize its applicability[J]. Journal of The Institution of Engineers (India): Series B, 2014, 95(4): 345-350.

    [15] 王千,王成,馮振元,等.K-means聚類算法研究綜述[J].電子設(shè)計(jì)工程,2012,20(7):21-24.

    WANG Qian, WANG Cheng, FENG Zhenyuan, et al. Review ofK-means clustering algorithm[J]. Electronic Design Engineering, 2012, 20(7): 21-24.

    [16] 周涓,熊忠陽(yáng),張玉芳,等.基于最大最小距離法的多中心聚類算法[J].計(jì)算機(jī)應(yīng)用,2006,26(6):1425-1427.

    ZHOU Juan, XIONG Zhongyang, ZHANG Yufang, et al. Multiseed clustering algorithm based on max-min distance means[J]. Computer Applications, 2006, 26(6): 1425-1427.

    [17] MORALES J M, PINEDA S, CONEJO A J, et al. A scenario reduction for futures market trading in electricity markets[J]. IEEE Trans on Power Systems, 2009, 24(2): 878-888.

    PeakLoadRegulatingOperationandOptimizationinPowerSystemswithLarge-scaleWindPowerandConsideringCouplingRelationbetweenWindPowerandLoad

    WUYaowu1,2,WANGChangshuang1,2,LOUSuhua1,2,LIUBaolin3,WANGXiran1,2

    (1. State Key Laboratory of Advanced Electromagnetic Engineering and Technology (Huazhong University of Science and Technology), Wuhan 430074, China; 2. School of Electrical and Electronics Engineering, Huazhong University of Science and Technology, Wuhan 430074, China; 3. Yunnan Electric Power Co. Ltd., Kunming 650011, China)

    Accurate modeling of stochastic wind power is one of the key links for peak load regulating operation and optimization in systems with large-scale wind power. According to the actual operation characteristics in China’s power systems, the interactive characteristic indices considering the coupling between wind power and load are proposed, and a wind power output model based on peak load regulation problems is put forward. This model uses key and clustering scenarios to characterize the influence of wind power on the power system operation and peak demand. Afterwards, an operation optimization model for peak load regulation with large-scale wind power based on the wind power output model is established. The expected value of various operation indices through operation optimization in multiple scenarios is obtained, and the analysis on peak load regulating operation characteristics is comprehensively accomplished. A case analysis is made on a province in West China that verified the validity and practicality of the proposed method.

    This work is supported by National Key Research and Development Program of China (No. 2016YFB0900102) and National Natural Science Foundation of China (No. 51677076).

    coupling relation; key scenarios; clustering scenarios; peak load regulation; operation and optimization

    2016-12-29;

    2017-04-23。

    上網(wǎng)日期: 2017-06-30。

    國(guó)家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(2016YFB0900102);國(guó)家自然科學(xué)基金資助項(xiàng)目(51677076)。

    吳耀武(1963—),男,博士,副教授,主要研究方向:電力系統(tǒng)規(guī)劃與優(yōu)化運(yùn)行。E-mail: ywwu@mail.hust.edu.cn

    汪昌霜(1992—),男,碩士研究生,主要研究方向:電力系統(tǒng)規(guī)劃與可靠性,新能源發(fā)電。E-mail: CS_Wang@hust.edu.cn

    婁素華(1974—),女,通信作者,博士、副教授,主要研究方向:電力系統(tǒng)規(guī)劃與優(yōu)化運(yùn)行、新能源發(fā)電與電力技術(shù)經(jīng)濟(jì)。E-mail: shlou@mail.hust.edu.cn

    (編輯孔麗蓓)

    猜你喜歡
    調(diào)峰出力風(fēng)電
    新常態(tài)下電站鍋爐深度調(diào)峰改造與調(diào)試實(shí)踐
    調(diào)峰保供型和普通型LNG接收站罐容計(jì)算
    煤氣與熱力(2021年5期)2021-07-22 09:02:14
    重慶市天然氣調(diào)峰儲(chǔ)氣建設(shè)的分析
    煤氣與熱力(2021年2期)2021-03-19 08:56:04
    海上風(fēng)電躍進(jìn)隱憂
    能源(2018年6期)2018-08-01 03:42:00
    分散式風(fēng)電破“局”
    能源(2018年6期)2018-08-01 03:41:56
    風(fēng)電:棄風(fēng)限電明顯改善 海上風(fēng)電如火如荼
    能源(2018年8期)2018-01-15 19:18:24
    關(guān)于寶雞市天然氣調(diào)峰問題的分析
    上海煤氣(2016年1期)2016-05-09 07:12:37
    風(fēng)電場(chǎng)有功出力的EEMD特性分析
    重齒風(fēng)電
    風(fēng)能(2016年12期)2016-02-25 08:46:38
    要爭(zhēng)做出力出彩的黨員干部
    河南電力(2016年5期)2016-02-06 02:11:35
    精品午夜福利在线看| 日韩一区二区视频免费看| 老司机影院毛片| 国产亚洲精品第一综合不卡| 伦精品一区二区三区| 最近最新中文字幕大全免费视频 | 欧美日本中文国产一区发布| 男的添女的下面高潮视频| 少妇人妻久久综合中文| 高清不卡的av网站| 亚洲成人一二三区av| 精品人妻在线不人妻| 精品国产一区二区三区久久久樱花| 26uuu在线亚洲综合色| 国产精品女同一区二区软件| 亚洲精品第二区| 成年女人在线观看亚洲视频| 日本免费在线观看一区| 久久99精品国语久久久| 激情五月婷婷亚洲| 男人爽女人下面视频在线观看| videos熟女内射| 女性生殖器流出的白浆| 亚洲美女黄色视频免费看| 国产视频首页在线观看| 九色亚洲精品在线播放| 精品一区二区三卡| 亚洲av电影在线进入| 欧美日韩一级在线毛片| 国产精品一区二区在线不卡| 狂野欧美激情性bbbbbb| 在线观看免费视频网站a站| 满18在线观看网站| 各种免费的搞黄视频| 国产免费一区二区三区四区乱码| 999久久久国产精品视频| 国产一区二区三区综合在线观看| 日本av免费视频播放| 在线观看国产h片| 少妇 在线观看| 久久这里有精品视频免费| 最近最新中文字幕大全免费视频 | 国产国语露脸激情在线看| 久久精品久久精品一区二区三区| 一区二区三区激情视频| 蜜桃在线观看..| 一级,二级,三级黄色视频| 大片电影免费在线观看免费| 五月天丁香电影| 丁香六月天网| 国产有黄有色有爽视频| 亚洲av男天堂| 亚洲欧美一区二区三区久久| 成年av动漫网址| 黄色配什么色好看| 观看av在线不卡| 最近中文字幕高清免费大全6| 国产亚洲一区二区精品| 国产伦理片在线播放av一区| 超色免费av| 国产亚洲最大av| 最近2019中文字幕mv第一页| 国产精品国产av在线观看| 1024视频免费在线观看| 久久ye,这里只有精品| 一级黄片播放器| 亚洲精华国产精华液的使用体验| 视频在线观看一区二区三区| www.熟女人妻精品国产| 在线免费观看不下载黄p国产| 欧美变态另类bdsm刘玥| 亚洲精品视频女| 91久久精品国产一区二区三区| 一二三四在线观看免费中文在| 少妇被粗大的猛进出69影院| 18在线观看网站| av一本久久久久| 久久久久国产精品人妻一区二区| 免费看不卡的av| 丝袜人妻中文字幕| 久久久久人妻精品一区果冻| 热re99久久精品国产66热6| 亚洲男人天堂网一区| 国产精品偷伦视频观看了| 久久综合国产亚洲精品| 美女视频免费永久观看网站| 九草在线视频观看| 观看av在线不卡| 久久久欧美国产精品| 人妻一区二区av| 久久99蜜桃精品久久| 人人妻人人添人人爽欧美一区卜| 999精品在线视频| 国产在线视频一区二区| 日韩不卡一区二区三区视频在线| 欧美日韩国产mv在线观看视频| 涩涩av久久男人的天堂| 亚洲国产精品一区二区三区在线| 高清av免费在线| 九色亚洲精品在线播放| 免费在线观看完整版高清| 欧美精品一区二区免费开放| 免费观看无遮挡的男女| 涩涩av久久男人的天堂| 国产女主播在线喷水免费视频网站| 色哟哟·www| 狂野欧美激情性bbbbbb| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜福利视频精品| 亚洲精品一区蜜桃| 中文字幕人妻丝袜制服| 午夜激情久久久久久久| 秋霞在线观看毛片| 精品酒店卫生间| 看非洲黑人一级黄片| 男男h啪啪无遮挡| tube8黄色片| 边亲边吃奶的免费视频| 伦理电影免费视频| 久久久精品区二区三区| 国产成人精品福利久久| 久久久久久伊人网av| 日本欧美国产在线视频| 成人午夜精彩视频在线观看| 国产精品国产三级专区第一集| 人人澡人人妻人| 中文字幕亚洲精品专区| 精品国产乱码久久久久久男人| 国产亚洲欧美精品永久| 日韩一区二区三区影片| 黄色一级大片看看| 日本vs欧美在线观看视频| 国产精品嫩草影院av在线观看| 80岁老熟妇乱子伦牲交| 一区在线观看完整版| 卡戴珊不雅视频在线播放| 91成人精品电影| 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区| 丝袜在线中文字幕| 99精国产麻豆久久婷婷| 亚洲精品国产一区二区精华液| 秋霞在线观看毛片| 大香蕉久久网| 中国三级夫妇交换| 97人妻天天添夜夜摸| 国产毛片在线视频| 成人国产麻豆网| 99久国产av精品国产电影| 欧美变态另类bdsm刘玥| 婷婷色综合大香蕉| 高清欧美精品videossex| 日韩一区二区视频免费看| 国产精品成人在线| 精品少妇内射三级| 久久狼人影院| 亚洲精品,欧美精品| 国产免费视频播放在线视频| 在线观看一区二区三区激情| 1024视频免费在线观看| 久久国产亚洲av麻豆专区| 波多野结衣av一区二区av| 日韩一本色道免费dvd| 在线观看美女被高潮喷水网站| 自线自在国产av| 成年美女黄网站色视频大全免费| 国产人伦9x9x在线观看 | 日韩中文字幕视频在线看片| 国产乱来视频区| 成人免费观看视频高清| 欧美精品国产亚洲| 妹子高潮喷水视频| 色网站视频免费| √禁漫天堂资源中文www| 蜜桃国产av成人99| 亚洲第一区二区三区不卡| 亚洲精品国产色婷婷电影| 国产无遮挡羞羞视频在线观看| 国产精品一区二区在线观看99| 免费久久久久久久精品成人欧美视频| 久久精品久久久久久久性| 亚洲人成77777在线视频| 亚洲内射少妇av| 婷婷色av中文字幕| 国产熟女欧美一区二区| 制服人妻中文乱码| 97在线人人人人妻| 亚洲精品久久久久久婷婷小说| 久久韩国三级中文字幕| 精品人妻在线不人妻| 少妇被粗大猛烈的视频| 亚洲一区中文字幕在线| 亚洲精品美女久久久久99蜜臀 | 中文字幕人妻熟女乱码| 日本欧美国产在线视频| 国产在线视频一区二区| 高清av免费在线| 久久久久久久大尺度免费视频| 黄频高清免费视频| 色视频在线一区二区三区| 丁香六月天网| 国产精品成人在线| 99久久精品国产国产毛片| 黄片小视频在线播放| 日韩一卡2卡3卡4卡2021年| 欧美日韩av久久| 大香蕉久久网| 欧美少妇被猛烈插入视频| 精品99又大又爽又粗少妇毛片| 亚洲一码二码三码区别大吗| 一二三四中文在线观看免费高清| 午夜福利乱码中文字幕| 少妇人妻 视频| 一级a爱视频在线免费观看| 成人免费观看视频高清| 91久久精品国产一区二区三区| 亚洲国产精品国产精品| 丝袜人妻中文字幕| 夜夜骑夜夜射夜夜干| 人妻一区二区av| 各种免费的搞黄视频| 亚洲av综合色区一区| 中文字幕人妻丝袜制服| 人人妻人人添人人爽欧美一区卜| 美女中出高潮动态图| 欧美人与性动交α欧美精品济南到 | 美女xxoo啪啪120秒动态图| 精品亚洲成国产av| 成人手机av| 2021少妇久久久久久久久久久| 欧美人与善性xxx| 99久久中文字幕三级久久日本| 国产亚洲一区二区精品| 亚洲国产精品一区二区三区在线| 啦啦啦在线免费观看视频4| 秋霞伦理黄片| 日日撸夜夜添| 亚洲av综合色区一区| 肉色欧美久久久久久久蜜桃| 在线观看国产h片| 一边摸一边做爽爽视频免费| 久久久久久久久久久免费av| 18禁裸乳无遮挡动漫免费视频| 日韩免费高清中文字幕av| 尾随美女入室| 自拍欧美九色日韩亚洲蝌蚪91| 女的被弄到高潮叫床怎么办| 母亲3免费完整高清在线观看 | 91精品三级在线观看| 亚洲精品国产一区二区精华液| 九九爱精品视频在线观看| 亚洲激情五月婷婷啪啪| 欧美 亚洲 国产 日韩一| 亚洲国产日韩一区二区| 欧美精品av麻豆av| 国产精品久久久av美女十八| 欧美日韩一级在线毛片| 精品少妇一区二区三区视频日本电影 | 亚洲欧美成人精品一区二区| av又黄又爽大尺度在线免费看| 亚洲精品国产av蜜桃| 亚洲精品自拍成人| 国产在视频线精品| 免费大片黄手机在线观看| av电影中文网址| 亚洲精品美女久久av网站| 多毛熟女@视频| 老女人水多毛片| 午夜av观看不卡| 欧美日韩精品网址| 美女脱内裤让男人舔精品视频| 最近中文字幕高清免费大全6| 80岁老熟妇乱子伦牲交| 国产精品人妻久久久影院| 狠狠精品人妻久久久久久综合| 欧美日韩成人在线一区二区| 9热在线视频观看99| 美女脱内裤让男人舔精品视频| 国产精品一区二区在线不卡| 亚洲精品国产av蜜桃| 成人亚洲欧美一区二区av| 精品国产一区二区久久| 亚洲,一卡二卡三卡| 国产精品一二三区在线看| 久久久久久久久久久免费av| 久久久久国产一级毛片高清牌| 久久久久精品性色| 国产精品不卡视频一区二区| 亚洲人成77777在线视频| 男人操女人黄网站| 中文字幕av电影在线播放| 久久精品国产鲁丝片午夜精品| 成年美女黄网站色视频大全免费| 日韩三级伦理在线观看| 九色亚洲精品在线播放| 日韩在线高清观看一区二区三区| 午夜福利在线观看免费完整高清在| 国产精品秋霞免费鲁丝片| 2021少妇久久久久久久久久久| 国产亚洲av片在线观看秒播厂| 久久99一区二区三区| 91久久精品国产一区二区三区| 国产一级毛片在线| 久久女婷五月综合色啪小说| 久久久久国产一级毛片高清牌| 少妇人妻精品综合一区二区| 久久精品国产亚洲av涩爱| 26uuu在线亚洲综合色| 99久久精品国产国产毛片| 80岁老熟妇乱子伦牲交| 日本色播在线视频| 国产精品久久久久久av不卡| 亚洲第一av免费看| 国产精品偷伦视频观看了| 丰满少妇做爰视频| 日韩电影二区| 欧美精品一区二区免费开放| 亚洲国产最新在线播放| 午夜福利乱码中文字幕| 26uuu在线亚洲综合色| 老司机亚洲免费影院| 99精国产麻豆久久婷婷| 91精品伊人久久大香线蕉| 国产在线视频一区二区| 黄片播放在线免费| 成人午夜精彩视频在线观看| 国产精品三级大全| 亚洲国产看品久久| 最近最新中文字幕免费大全7| 精品99又大又爽又粗少妇毛片| 国产精品一国产av| 麻豆精品久久久久久蜜桃| 男女啪啪激烈高潮av片| 免费日韩欧美在线观看| 日本免费在线观看一区| 天天躁日日躁夜夜躁夜夜| 人人妻人人澡人人爽人人夜夜| 国精品久久久久久国模美| 欧美人与善性xxx| 亚洲国产精品一区三区| 热99久久久久精品小说推荐| 国产男人的电影天堂91| 一个人免费看片子| 久久久久久伊人网av| 久久精品国产a三级三级三级| 国产欧美日韩综合在线一区二区| 黄色 视频免费看| 香蕉国产在线看| 性色av一级| 视频在线观看一区二区三区| 秋霞在线观看毛片| 亚洲精品视频女| 最近的中文字幕免费完整| 一级片免费观看大全| 亚洲欧美中文字幕日韩二区| 亚洲av日韩在线播放| 亚洲精品国产av成人精品| 欧美日韩综合久久久久久| 99九九在线精品视频| 人人妻人人澡人人爽人人夜夜| 国产免费视频播放在线视频| 国产成人精品在线电影| 丝袜美足系列| 国产成人精品在线电影| 99久久综合免费| 亚洲人成77777在线视频| 99久久综合免费| 人人妻人人澡人人爽人人夜夜| av卡一久久| 国产免费现黄频在线看| 国产成人免费无遮挡视频| 在线天堂中文资源库| 久久精品久久精品一区二区三区| 国产一区二区激情短视频 | 五月开心婷婷网| 黄色毛片三级朝国网站| 性高湖久久久久久久久免费观看| 在线观看免费高清a一片| 国产成人精品福利久久| 在线观看免费高清a一片| 久久精品人人爽人人爽视色| 午夜91福利影院| 人体艺术视频欧美日本| 午夜老司机福利剧场| 大陆偷拍与自拍| 日韩欧美精品免费久久| 中文字幕另类日韩欧美亚洲嫩草| 99久久精品国产国产毛片| 国产免费福利视频在线观看| 人妻少妇偷人精品九色| av网站免费在线观看视频| 黄色视频在线播放观看不卡| 亚洲激情五月婷婷啪啪| 韩国av在线不卡| 91国产中文字幕| 妹子高潮喷水视频| 观看av在线不卡| 熟妇人妻不卡中文字幕| 少妇人妻 视频| 亚洲少妇的诱惑av| 亚洲精品乱久久久久久| 美女中出高潮动态图| 色哟哟·www| 蜜桃在线观看..| 国产综合精华液| 捣出白浆h1v1| 大片免费播放器 马上看| 人人妻人人澡人人爽人人夜夜| 午夜91福利影院| 免费观看在线日韩| 国产老妇伦熟女老妇高清| 伊人亚洲综合成人网| 少妇 在线观看| 久久99蜜桃精品久久| 欧美成人精品欧美一级黄| 欧美bdsm另类| 十八禁高潮呻吟视频| 国产97色在线日韩免费| 国产色婷婷99| 国产精品一区二区在线不卡| 精品人妻熟女毛片av久久网站| 亚洲,一卡二卡三卡| 亚洲av电影在线进入| 又黄又粗又硬又大视频| 国产成人精品无人区| 国产黄色视频一区二区在线观看| 欧美精品一区二区免费开放| www.自偷自拍.com| 一级毛片 在线播放| 国产免费又黄又爽又色| av一本久久久久| 91精品三级在线观看| 午夜福利视频精品| 欧美激情高清一区二区三区 | 久久精品aⅴ一区二区三区四区 | 国产av一区二区精品久久| 亚洲国产精品一区三区| 黄色配什么色好看| 亚洲欧美一区二区三区久久| 边亲边吃奶的免费视频| 国产精品国产三级国产专区5o| 哪个播放器可以免费观看大片| 狠狠精品人妻久久久久久综合| 亚洲av男天堂| 久热这里只有精品99| 亚洲精品乱久久久久久| 午夜影院在线不卡| 亚洲精品国产一区二区精华液| 国产成人午夜福利电影在线观看| 校园人妻丝袜中文字幕| 精品国产一区二区三区久久久樱花| 老司机影院成人| 91午夜精品亚洲一区二区三区| 免费看av在线观看网站| 一区二区av电影网| 欧美在线黄色| 9色porny在线观看| 国产精品一区二区在线观看99| 不卡av一区二区三区| 亚洲精品视频女| 欧美 日韩 精品 国产| av片东京热男人的天堂| 日韩精品免费视频一区二区三区| 人人澡人人妻人| 90打野战视频偷拍视频| 老司机影院成人| 欧美日韩成人在线一区二区| 韩国精品一区二区三区| 熟女少妇亚洲综合色aaa.| 久久av网站| 99久久综合免费| 成人毛片60女人毛片免费| 黄网站色视频无遮挡免费观看| 久久精品久久久久久噜噜老黄| videos熟女内射| 欧美另类一区| 伦理电影大哥的女人| 超碰成人久久| 亚洲精品久久久久久婷婷小说| 在线亚洲精品国产二区图片欧美| 日韩三级伦理在线观看| 免费在线观看完整版高清| 嫩草影院入口| 国产女主播在线喷水免费视频网站| 久久久精品区二区三区| 日日撸夜夜添| 午夜免费男女啪啪视频观看| 啦啦啦视频在线资源免费观看| av.在线天堂| 99热全是精品| 日韩一卡2卡3卡4卡2021年| 日韩精品免费视频一区二区三区| 日本av手机在线免费观看| 老司机影院毛片| 欧美精品亚洲一区二区| 婷婷色综合大香蕉| 国产xxxxx性猛交| 国产日韩一区二区三区精品不卡| 欧美最新免费一区二区三区| 最近2019中文字幕mv第一页| 久久久久精品久久久久真实原创| 国产淫语在线视频| 校园人妻丝袜中文字幕| 午夜福利在线观看免费完整高清在| 国产精品国产三级国产专区5o| 国产精品秋霞免费鲁丝片| 成年女人在线观看亚洲视频| av在线播放精品| 寂寞人妻少妇视频99o| 在线观看www视频免费| 七月丁香在线播放| 宅男免费午夜| 一级毛片我不卡| 国产黄色免费在线视频| 伊人久久国产一区二区| 国产人伦9x9x在线观看 | 国产精品久久久久久精品电影小说| 国产精品一国产av| 国产精品久久久av美女十八| 综合色丁香网| 寂寞人妻少妇视频99o| 超碰97精品在线观看| 亚洲一级一片aⅴ在线观看| a级片在线免费高清观看视频| 国产精品久久久久久av不卡| 男女啪啪激烈高潮av片| 乱人伦中国视频| 国产精品三级大全| 久久久久久久久久久免费av| 狂野欧美激情性bbbbbb| 精品一区二区三区四区五区乱码 | 麻豆乱淫一区二区| av片东京热男人的天堂| 久久亚洲国产成人精品v| 欧美日韩国产mv在线观看视频| 亚洲av国产av综合av卡| 久久精品国产亚洲av天美| √禁漫天堂资源中文www| 伦理电影大哥的女人| 水蜜桃什么品种好| 精品一区二区三区四区五区乱码 | 一区二区三区激情视频| 中国三级夫妇交换| 中文字幕人妻丝袜一区二区 | 男人添女人高潮全过程视频| xxxhd国产人妻xxx| 成人毛片60女人毛片免费| 久久久久久人人人人人| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产最新在线播放| 只有这里有精品99| 丝袜喷水一区| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品一区二区三区在线| 国产精品蜜桃在线观看| 亚洲精品国产一区二区精华液| 美女高潮到喷水免费观看| 亚洲欧美一区二区三区国产| 精品一区二区免费观看| 亚洲精品一二三| 一本久久精品| 男人操女人黄网站| 久久精品熟女亚洲av麻豆精品| 啦啦啦在线免费观看视频4| 美女大奶头黄色视频| 自线自在国产av| 欧美精品av麻豆av| 欧美激情高清一区二区三区 | 国产成人精品福利久久| 亚洲欧美精品综合一区二区三区 | 国产一级毛片在线| 日韩制服骚丝袜av| 波多野结衣av一区二区av| 国产麻豆69| 色哟哟·www| 午夜免费男女啪啪视频观看| 丰满迷人的少妇在线观看| 伦精品一区二区三区| 国精品久久久久久国模美| 最近中文字幕2019免费版| 欧美人与性动交α欧美软件| 边亲边吃奶的免费视频| 亚洲精品成人av观看孕妇| 国产一区二区三区av在线| 欧美成人午夜免费资源| 亚洲一级一片aⅴ在线观看| 国产精品香港三级国产av潘金莲 | 亚洲第一青青草原| 国产1区2区3区精品| 日本-黄色视频高清免费观看| 9热在线视频观看99| 中文欧美无线码| 亚洲美女搞黄在线观看| 亚洲av福利一区| 欧美激情 高清一区二区三区| 一二三四中文在线观看免费高清| 狠狠婷婷综合久久久久久88av| 黄色视频在线播放观看不卡| 亚洲色图 男人天堂 中文字幕| 亚洲伊人色综图| www.精华液| 亚洲精品美女久久av网站| av卡一久久| 伊人亚洲综合成人网| av在线app专区| 午夜免费鲁丝| 天天躁夜夜躁狠狠久久av| 免费黄频网站在线观看国产| 高清不卡的av网站| 超色免费av| 成年女人在线观看亚洲视频| 咕卡用的链子| 欧美日韩一级在线毛片| 日韩在线高清观看一区二区三区| 国产精品嫩草影院av在线观看|