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

    計(jì)及運(yùn)行風(fēng)險(xiǎn)的電力系統(tǒng)優(yōu)化調(diào)度方法

    2016-08-11 09:29:50宋云亭

    李 燦,曾 沅,秦 超,宋云亭,吉 平,吳 威

    (1.天津大學(xué)智能電網(wǎng)教育部重點(diǎn)實(shí)驗(yàn)室,天津 300072;2.中國(guó)電力科學(xué)研究院,北京 100192;3.國(guó)網(wǎng)福建省電力有限公司經(jīng)濟(jì)技術(shù)研究院,福州 350012)

    ?

    計(jì)及運(yùn)行風(fēng)險(xiǎn)的電力系統(tǒng)優(yōu)化調(diào)度方法

    李燦1,曾沅1,秦超1,宋云亭2,吉平2,吳威3

    (1.天津大學(xué)智能電網(wǎng)教育部重點(diǎn)實(shí)驗(yàn)室,天津 300072;2.中國(guó)電力科學(xué)研究院,北京 100192;3.國(guó)網(wǎng)福建省電力有限公司經(jīng)濟(jì)技術(shù)研究院,福州 350012)

    摘 要:針對(duì)現(xiàn)有優(yōu)化調(diào)度方法在考慮風(fēng)險(xiǎn)因素時(shí)過(guò)于單一、片面的現(xiàn)狀,該文考慮風(fēng)電波動(dòng)、機(jī)組停運(yùn)、線路故障、負(fù)荷波動(dòng)等多種不確定因素的影響,以電壓越限和潮流越限衡量系統(tǒng)整體的運(yùn)行風(fēng)險(xiǎn),構(gòu)建了能夠計(jì)及系統(tǒng)面臨風(fēng)險(xiǎn)的多目標(biāo)優(yōu)化調(diào)度模型。該模型以運(yùn)行風(fēng)險(xiǎn)和發(fā)電成本為目標(biāo)函數(shù),將隨機(jī)潮流理論引入到發(fā)電優(yōu)化調(diào)度中,并針對(duì)不同的天氣狀況,對(duì)線路故障概率進(jìn)行實(shí)時(shí)修正以期更加接近實(shí)際運(yùn)行環(huán)境。結(jié)合多目標(biāo)粒子群算法及模糊集理論對(duì)所建立的模型進(jìn)行求解,算例結(jié)果表明所提模型和方法能夠處理不確定因素對(duì)系統(tǒng)運(yùn)行風(fēng)險(xiǎn)的影響,并反映不同天氣狀況下調(diào)度方案的差異性,為實(shí)施短期優(yōu)化調(diào)度提供有益參考。

    關(guān)鍵詞:優(yōu)化調(diào)度;運(yùn)行風(fēng)險(xiǎn);不確定因素;隨機(jī)潮流

    隨著我國(guó)區(qū)域電網(wǎng)的互聯(lián),風(fēng)電并網(wǎng)容量的逐年增加以及電力市場(chǎng)的逐步建設(shè),電力系統(tǒng)優(yōu)化調(diào)度問(wèn)題面臨的復(fù)雜性日益增加。傳統(tǒng)的電力系統(tǒng)優(yōu)化調(diào)度問(wèn)題通常是在考慮系統(tǒng)功率平衡、機(jī)組出力不越限等一系列約束條件的前提下,以系統(tǒng)總發(fā)電費(fèi)用最小為目標(biāo),對(duì)發(fā)電廠各機(jī)組出力進(jìn)行優(yōu)化分配[1-2]。但往往過(guò)于關(guān)注經(jīng)濟(jì)性而忽略了調(diào)度運(yùn)行中的其他因素,無(wú)法適應(yīng)電力系統(tǒng)的發(fā)展需要。

    文獻(xiàn)[3-5]以節(jié)能、經(jīng)濟(jì)、環(huán)保三大原則為目標(biāo),建立了多目標(biāo)優(yōu)化調(diào)度模型。但這些方法均未對(duì)系統(tǒng)的安全可靠性給予太多的重視。目前有關(guān)電網(wǎng)風(fēng)險(xiǎn)的研究主要圍繞靜態(tài)安全評(píng)估[6-7]、輸電斷面輸電能力分析[8]、暫態(tài)安全評(píng)估[9-10]等方面,取得了大量成果。而近年來(lái)在優(yōu)化調(diào)度領(lǐng)域引入風(fēng)險(xiǎn)理論已成為研究人員的研究熱點(diǎn)。文獻(xiàn)[11]引入棄風(fēng)懲罰因子,建立了考慮系統(tǒng)安全風(fēng)險(xiǎn)的優(yōu)化調(diào)度模型,但該模型只關(guān)注了風(fēng)機(jī)出力對(duì)系統(tǒng)風(fēng)險(xiǎn)的影響,而忽略了影響系統(tǒng)風(fēng)險(xiǎn)的其他因素。文獻(xiàn)[12]考慮發(fā)電資源消耗最少、能源環(huán)境效益最好、系統(tǒng)安全穩(wěn)定程度最高等因素,提出了含風(fēng)電場(chǎng)的多目

    5)系統(tǒng)旋轉(zhuǎn)備用約束標(biāo)優(yōu)化調(diào)度模型,但該模型僅以電壓穩(wěn)定指標(biāo)衡量系統(tǒng)的安全性。文獻(xiàn)[13]定義了運(yùn)行風(fēng)險(xiǎn),針對(duì)線路故障,構(gòu)建了基于風(fēng)險(xiǎn)的多目標(biāo)經(jīng)濟(jì)調(diào)度模型,但未考慮風(fēng)電出力的波動(dòng)性對(duì)于系統(tǒng)運(yùn)行風(fēng)險(xiǎn)的影響。文獻(xiàn)[14-15]建立了計(jì)及風(fēng)險(xiǎn)的風(fēng)電并網(wǎng)系統(tǒng)的優(yōu)化調(diào)度模型,但涉及的風(fēng)險(xiǎn)均為考慮風(fēng)電并網(wǎng)影響的系統(tǒng)備用風(fēng)險(xiǎn),并未考慮系統(tǒng)整體的運(yùn)行風(fēng)險(xiǎn)。文獻(xiàn)[16-17]計(jì)及冰凍災(zāi)害對(duì)輸電線路的影響,建立了考慮天氣因素的機(jī)組優(yōu)化方法。綜上,現(xiàn)有的優(yōu)化調(diào)度方法在考慮風(fēng)險(xiǎn)因素時(shí)大多只考慮了個(gè)別不確定因素的影響。

    本文考慮風(fēng)電波動(dòng)、機(jī)組停運(yùn)、線路故障、負(fù)荷波動(dòng)等多種不確定因素的影響,以電壓越限和潮流越限衡量系統(tǒng)的運(yùn)行風(fēng)險(xiǎn),構(gòu)建了計(jì)及運(yùn)行風(fēng)險(xiǎn)的多目標(biāo)優(yōu)化調(diào)度模型。該模型以系統(tǒng)運(yùn)行風(fēng)險(xiǎn)和發(fā)電成本為目標(biāo)函數(shù),通過(guò)將隨機(jī)潮流理論引入優(yōu)化調(diào)度問(wèn)題,以其宏觀統(tǒng)計(jì)特性處理系統(tǒng)中不確定因素對(duì)運(yùn)行風(fēng)險(xiǎn)的影響,并在計(jì)算中針對(duì)不同天氣狀況對(duì)線路故障概率進(jìn)行修正,從而更加接近實(shí)際運(yùn)行環(huán)境。為求解該模型,采用多目標(biāo)粒子群算法與模糊集理論來(lái)確定最優(yōu)方案,最后以IEEE-14節(jié)點(diǎn)測(cè)試系統(tǒng)為例來(lái)證明本文所述方法的有效性。

    1 計(jì)及運(yùn)行風(fēng)險(xiǎn)的優(yōu)化調(diào)度模型

    電力系統(tǒng)優(yōu)化調(diào)度分為單時(shí)段靜態(tài)優(yōu)化調(diào)度和多時(shí)段動(dòng)態(tài)優(yōu)化調(diào)度。多時(shí)段動(dòng)態(tài)優(yōu)化調(diào)度是在單時(shí)段靜態(tài)優(yōu)化調(diào)度的基礎(chǔ)上考慮了不同時(shí)間斷面的耦合性。為簡(jiǎn)化計(jì)算,驗(yàn)證本文所提方法的可行性,文中僅針對(duì)單時(shí)段建立了優(yōu)化調(diào)度模型。

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

    電力系統(tǒng)中節(jié)點(diǎn)電壓和支路潮流是否處于合理的區(qū)間內(nèi),對(duì)電力系統(tǒng)的運(yùn)行安全性至關(guān)重要,因此節(jié)點(diǎn)電壓和支路潮流的越限風(fēng)險(xiǎn)可以很好地衡量系統(tǒng)整體的運(yùn)行風(fēng)險(xiǎn)。本文以系統(tǒng)的運(yùn)行風(fēng)險(xiǎn)和各火電廠的發(fā)電成本為目標(biāo)函數(shù),構(gòu)建了計(jì)及運(yùn)行風(fēng)險(xiǎn)的多目標(biāo)優(yōu)化調(diào)度模型,其數(shù)學(xué)表達(dá)式為

    式中:Ng為系統(tǒng)中常規(guī)發(fā)電機(jī)組的數(shù)量;pi為機(jī)組i的有功出力;fi(pi)為機(jī)組i的耗量特性,一般用二次函數(shù)近似表示,即 fi(pi)=aip2i+bipi+ci,其中ai、bi、ci為機(jī)組i的耗量特性參數(shù);Nl為系統(tǒng)中支路的數(shù)量;Rli為第i條支路的越限風(fēng)險(xiǎn);Nn為系統(tǒng)中PQ節(jié)點(diǎn)的數(shù)量;Rui為第i個(gè)PQ節(jié)點(diǎn)的電壓越限風(fēng)險(xiǎn)。

    1.2約束條件

    1)系統(tǒng)功率平衡約束

    式中:Pi、Qi分別為節(jié)點(diǎn)i的發(fā)電機(jī)有功出力和無(wú)功出力;Pw、Qw分別為風(fēng)機(jī)的有功出力和無(wú)功出力;Pdi、Qdi分別為節(jié)點(diǎn)i的有功負(fù)荷和無(wú)功負(fù)荷;φij為節(jié)點(diǎn)i和 j的電壓相角差;Gij、Bij分別為支路i、j的電導(dǎo)和電納。

    2)常規(guī)發(fā)電機(jī)出力約束

    式中,Pim

    ax、Pimin和Qmiax、Qmiin分別為發(fā)電機(jī)i有功和無(wú)功出力的上下限。

    3)節(jié)點(diǎn)電壓約束

    式中,Uim

    ax、Uimin分別為節(jié)點(diǎn)i的電壓上下限。

    4)線路傳輸功率約束

    式中,Smaxli為線路i的傳輸極限。

    式中:Sui、Sd

    i分別為常規(guī)機(jī)組i提供的正、負(fù)旋轉(zhuǎn)備用容量;PL為系統(tǒng)的總負(fù)荷;L%分別為系統(tǒng)總負(fù)荷預(yù)測(cè)誤差對(duì)正旋轉(zhuǎn)備用的需求;Pwmax為風(fēng)機(jī)的最大出力;Pw為風(fēng)機(jī)的有功出力;wu%、wd%分別為風(fēng)機(jī)出力預(yù)測(cè)誤差對(duì)正、負(fù)旋轉(zhuǎn)備用的需求。

    2 基于隨機(jī)潮流的運(yùn)行風(fēng)險(xiǎn)

    隨機(jī)潮流可以考慮多種不確定因素的影響,以其數(shù)學(xué)模型為基礎(chǔ),通過(guò)概率理論得到系統(tǒng)運(yùn)行狀況的宏觀統(tǒng)計(jì)信息。相比于傳統(tǒng)的確定性潮流分析,能夠以較小的計(jì)算量更加全面、深刻地展現(xiàn)系統(tǒng)的運(yùn)行特性。因此本節(jié)采用隨機(jī)潮流處理風(fēng)電波動(dòng)、負(fù)荷波動(dòng)、機(jī)組停運(yùn)、線路故障等不確定因素,并以此為基礎(chǔ)計(jì)算優(yōu)化調(diào)度模型中的運(yùn)行風(fēng)險(xiǎn)F2。

    2.1不確定因素的數(shù)學(xué)模型

    1)風(fēng)電出力的不確定性

    風(fēng)機(jī)的種類多種多樣,其中雙饋風(fēng)電機(jī)組是目前主流的風(fēng)電機(jī)組,其有功出力和風(fēng)速的關(guān)系為

    式中:Pr為風(fēng)電機(jī)組額定有功出力;vci為切入風(fēng)速;vr為額定風(fēng)速;vco為切出風(fēng)速;k1、k2為風(fēng)速介于切入風(fēng)速和額定風(fēng)速之間時(shí)風(fēng)機(jī)有功出力函數(shù)的相關(guān)系數(shù),具體表達(dá)式為

    根據(jù)實(shí)際的統(tǒng)計(jì),一般情況下風(fēng)速都是維持在切入風(fēng)速和額定風(fēng)速的范圍內(nèi),根據(jù)在此區(qū)間內(nèi)風(fēng)機(jī)有功出力與風(fēng)速之間的近似線性關(guān)系,并結(jié)合風(fēng)速的威布爾分布模型,可以得到風(fēng)機(jī)有功出力的概率密度函數(shù)為

    式中,k、c分別為威布爾分布的形狀參數(shù)和尺度參數(shù)。

    雙饋風(fēng)機(jī)一般以恒功率因數(shù)方式運(yùn)行,因此其無(wú)功出力的概率密度函數(shù)可表示為

    式中,α為功率因數(shù)角。

    2)負(fù)荷的波動(dòng)性

    大量相關(guān)文獻(xiàn)以及長(zhǎng)期的實(shí)踐表明,負(fù)荷的隨機(jī)波動(dòng)服從正態(tài)分布。其概率密度函數(shù)可表示為

    式中:μP、σP為有功功率的期望值和標(biāo)準(zhǔn)差;μQ、σQ為無(wú)功功率的期望值和標(biāo)準(zhǔn)差。

    3)常規(guī)發(fā)電機(jī)的強(qiáng)迫停運(yùn)

    常規(guī)發(fā)電機(jī)模型只有正常運(yùn)行和故障停運(yùn)兩種狀態(tài),因此其隨機(jī)特性可以用二項(xiàng)分布來(lái)描述:

    式中:p為發(fā)電機(jī)組的可用概率;SN為發(fā)電機(jī)組的額定容量。

    4)線路故障

    在中長(zhǎng)期風(fēng)險(xiǎn)評(píng)估中通常是根據(jù)多年的統(tǒng)計(jì)數(shù)據(jù)得到線路的故障概率,并未考慮到不同時(shí)間不同區(qū)域線路故障概率的差別;在實(shí)時(shí)運(yùn)行風(fēng)險(xiǎn)評(píng)估中為更加接近系統(tǒng)的實(shí)際運(yùn)行環(huán)境,有必要對(duì)線路的故障概率進(jìn)行實(shí)時(shí)修正。

    針對(duì)不同的天氣狀況,對(duì)線路故障概率按如下方法進(jìn)行修正[18],即

    式中:p(line)為修正后的線路故障概率;p0(line)為歷史統(tǒng)計(jì)線路故障概率;w為天氣影響因子。

    天氣影響因子的取值根據(jù)天氣的不同而有所不同,具體如表1所示。

    表1 不同天氣下的天氣影響因子Tab.1 Weather factors under different weather

    通過(guò)修正后的線路故障概率,可以得到第i種線路故障情況發(fā)生的概率p(li)為

    式中,m、n分別為非故障線路和故障線路。

    對(duì)于由指定的故障線路組成的預(yù)想事故集,各種線路故障情況發(fā)生的概率之和必然為1,因此第i種線路故障情況發(fā)生的概率可以定義為

    2.2隨機(jī)潮流

    電力系統(tǒng)的節(jié)點(diǎn)功率和支路潮流可概括為

    式中:W為節(jié)點(diǎn)注入功率;Z為支路潮流;X為節(jié)點(diǎn)電壓;f為功率方程;g為支路潮流方程。

    利用泰勒級(jí)數(shù)將式(16)在系統(tǒng)基準(zhǔn)運(yùn)行點(diǎn)處展開(kāi),忽略二次及以上高次項(xiàng)可得

    式中:ΔW、ΔX、ΔZ表示注入功率的隨機(jī)擾動(dòng)量以及對(duì)應(yīng)的節(jié)點(diǎn)電壓和支路潮流的隨機(jī)變化量;J0表示潮流計(jì)算最后一次迭代所用的雅可比矩陣;,X0表示基準(zhǔn)運(yùn)行點(diǎn)處電壓的期望值。

    式(17)即為隨機(jī)潮流的數(shù)學(xué)模型,而為了處理線路故障,本文借鑒文獻(xiàn)[19]的全概率理論。

    全概率理論的表達(dá)式為

    將其應(yīng)用到電力系統(tǒng)隨機(jī)潮流中,表達(dá)式各變量被賦予電力系統(tǒng)新的含義:P(A)為考慮線路故障的系統(tǒng)隨機(jī)潮流;Bi為第i種線路故障情況;P(Bi)為第i種線路故障情況發(fā)生的概率p′(li);P(A|Bi)為第i種線路故障情況下的系統(tǒng)隨機(jī)潮流。

    本文為減少計(jì)算量,采用半不變量法與Gram-Charlier級(jí)數(shù)相結(jié)合,對(duì)隨機(jī)潮流模型進(jìn)行求解,近似擬合出系統(tǒng)各節(jié)點(diǎn)電壓、支路潮流的概率密度函數(shù)和累積分布函數(shù)。

    2.3系統(tǒng)運(yùn)行風(fēng)險(xiǎn)

    在實(shí)際電網(wǎng)運(yùn)行中,當(dāng)某一節(jié)點(diǎn)的電壓或者某一支路的潮流發(fā)生越限時(shí),可能會(huì)造成相應(yīng)的保護(hù)裝置動(dòng)作,導(dǎo)致在保護(hù)裝置動(dòng)作期間系統(tǒng)出現(xiàn)一定程度的負(fù)荷損失。在電力市場(chǎng)環(huán)境下,可以采用貨幣的形式來(lái)衡量系統(tǒng)的負(fù)荷損失,因此,當(dāng)系統(tǒng)出現(xiàn)越限情況時(shí),其后果的嚴(yán)重程度可表示[20]為

    式中:當(dāng)越限情況i為節(jié)點(diǎn)i電壓越限時(shí),Li為節(jié)點(diǎn)i處的切負(fù)荷量,當(dāng)越限情況i為支路i潮流越限時(shí),Li為支路i越限導(dǎo)致的切負(fù)荷量;Ti、Ci分別為越限情況i導(dǎo)致的停電事故的恢復(fù)時(shí)間和單位平均成本;Di為因供電合同中斷而帶來(lái)的懲罰成本。

    節(jié)點(diǎn)電壓和支路潮流越限的概率[20-21]分別為

    式中:Uim

    ax、Uimin分別為節(jié)點(diǎn)i的電壓上下限;F(x)為累積分布函數(shù);1-F(Uimax)為節(jié)點(diǎn)i電壓超過(guò)上限的概率;F(Uimin)為節(jié)點(diǎn)i電壓低于下限的概率;

    pmlia

    x為線路i輸送有功功率的上限;1-F(pmliax)為線路i潮流超過(guò)上限的概率。

    根據(jù)上述節(jié)點(diǎn)電壓和支路潮流的越限概率及相應(yīng)的后果嚴(yán)重度,可以得到節(jié)點(diǎn)電壓越限和支路潮流越限的風(fēng)險(xiǎn)分別為

    通過(guò)式(22)即可求得所有PQ節(jié)點(diǎn)的電壓越限風(fēng)險(xiǎn)和所有支路的潮流越限風(fēng)險(xiǎn),從而得到系統(tǒng)整體的運(yùn)行風(fēng)險(xiǎn)F2。

    3 優(yōu)化調(diào)度模型的求解

    許多文獻(xiàn)在求解多目標(biāo)優(yōu)化問(wèn)題時(shí)大多沒(méi)有考慮各優(yōu)化目標(biāo)之間的沖突性,只是簡(jiǎn)單地采用加權(quán)的方法轉(zhuǎn)化為單目標(biāo)問(wèn)題進(jìn)行求解。因此并不能準(zhǔn)確地求解多目標(biāo)優(yōu)化調(diào)度問(wèn)題。多目標(biāo)粒子群算法(MPSO)是一種基于群體智能的演化算法,針對(duì)多個(gè)目標(biāo)函數(shù),采用Pareto支配的概念構(gòu)造非支配解集,即

    若決策向量P1、P2都屬于可行解,那么當(dāng)且僅當(dāng)式(23)成立時(shí),稱P1支配P2。

    將當(dāng)前找到的非支配解按以下規(guī)則保存在外部檔案中:當(dāng)外部檔案?jìng)€(gè)體數(shù)小于設(shè)定值時(shí),非支配解直接加入外部檔案;當(dāng)外部檔案?jìng)€(gè)體數(shù)達(dá)到設(shè)定值時(shí),如果新解支配了外部檔案的部分個(gè)體,則新解取代受支配的那些個(gè)體,否則根據(jù)自適應(yīng)網(wǎng)格法[22]移除一個(gè)個(gè)體,將新解加入檔案中。隨著迭代過(guò)程不斷維護(hù)外部檔案,最后,根據(jù)一定的策略從外部檔案中選取一個(gè)折中解作為優(yōu)化問(wèn)題的最優(yōu)解,將其應(yīng)用到本文的電力系統(tǒng)優(yōu)化調(diào)度模型中,具體流程如圖1所示。最終得到的外部檔案即為優(yōu)化調(diào)度方案的Pareto最優(yōu)解集。

    本文采用模糊決策方法從Pareto最優(yōu)解集中選取一個(gè)最優(yōu)折中解作為最終調(diào)度方案[23]。具體做法如下:

    對(duì)每個(gè)目標(biāo)函數(shù)進(jìn)行模糊化處理,其隸屬度函數(shù)為

    式中,fimax、fimin分別為第i個(gè)目標(biāo)函數(shù)的最大值和最小值,其中i∈No,No為目標(biāo)函數(shù)的個(gè)數(shù)。

    每個(gè)Pareto最優(yōu)解的每個(gè)目標(biāo)函數(shù)的隸屬度對(duì)應(yīng)于調(diào)度運(yùn)行人員對(duì)該目標(biāo)函數(shù)的滿意度。而調(diào)度運(yùn)行人員對(duì)于每個(gè)Pareto最優(yōu)解的滿意度可表示為

    通過(guò)計(jì)算Pareto最優(yōu)解集中每個(gè)解的M值,并進(jìn)行大小排序,選出M值最大的Pareto最優(yōu)解作為該優(yōu)化調(diào)度問(wèn)題的最優(yōu)折中解。

    圖1 多目標(biāo)粒子群算法流程Fig.1 Flow chart of MPSO

    4 算例分析

    采用IEEE-14節(jié)點(diǎn)系統(tǒng),對(duì)本文所提出的方法進(jìn)行驗(yàn)證。該系統(tǒng)的接線圖和具體參數(shù)見(jiàn)文獻(xiàn)[23],所有功率數(shù)據(jù)都是以100 MVA為功率基值的標(biāo)幺值。對(duì)于發(fā)電機(jī)節(jié)點(diǎn),除2號(hào)節(jié)點(diǎn)的發(fā)電機(jī)的強(qiáng)迫停運(yùn)率為0.09,其余的發(fā)電機(jī)節(jié)點(diǎn)的發(fā)電機(jī)的強(qiáng)迫停運(yùn)率均為0.08。負(fù)荷均以IEEE-14節(jié)點(diǎn)系統(tǒng)的負(fù)荷值為均值,取其10%為標(biāo)準(zhǔn)差,且各負(fù)荷相互獨(dú)立。假定負(fù)荷預(yù)測(cè)最大誤差和風(fēng)電功率最大誤差對(duì)旋轉(zhuǎn)備用的需求分別為10%和30%,各節(jié)點(diǎn)電壓的閾值設(shè)為[0.95,1.05],各支路潮流的傳輸極限為優(yōu)化前其傳輸功率的1.2倍。對(duì)于節(jié)點(diǎn)電壓越限,假定切除該節(jié)點(diǎn)30%的負(fù)荷后,系統(tǒng)恢復(fù)正常運(yùn)行[21];對(duì)于支路潮流越限,則假定就近切除線路兩端30%的負(fù)荷后,系統(tǒng)恢復(fù)正常運(yùn)行[20]。

    為計(jì)及風(fēng)電場(chǎng)對(duì)于系統(tǒng)的影響,在系統(tǒng)的14號(hào)節(jié)點(diǎn)處通過(guò)變壓器和110 kV輸電線路引入風(fēng)電場(chǎng)。風(fēng)電場(chǎng)采用雙饋異步風(fēng)機(jī),切入風(fēng)速是3 m/s,切出風(fēng)速是25 m/s,額定風(fēng)速是14 m/s,其接入部分的參數(shù)見(jiàn)表2。

    表2 風(fēng)電場(chǎng)的相關(guān)參數(shù)Tab.2 Related parameters of wind farm

    對(duì)于線路故障,假設(shè)該系統(tǒng)中只有線路2~4和線路7~9易發(fā)生故障,且前者易遭遇雷雨天氣,后者易遭遇大霧天氣,表3給出了兩條線路在正常天氣下以及在惡劣天氣下經(jīng)天氣影響因子修正后的線路故障概率。

    表3 不同天氣下線路故障概率Tab.3 Line fault probabilities under different weather conditions

    在優(yōu)化前,分別針對(duì)未接入風(fēng)電、接入風(fēng)電、接入風(fēng)電并考慮線路故障3種情況對(duì)系統(tǒng)進(jìn)行隨機(jī)潮流分析,以節(jié)點(diǎn)4的電壓為例,其概率密度函數(shù)(PDF)和累積分布函數(shù)(CDF)如圖2和圖3所示??梢钥闯?,風(fēng)電的接入以及線路故障增大了電壓的期望值及波動(dòng)范圍,從而增大了系統(tǒng)的越限風(fēng)險(xiǎn)。

    圖2 節(jié)點(diǎn)4電壓概率密度函數(shù)Fig.2 Voltage PDF of node 4

    圖3 節(jié)點(diǎn)4電壓累積分布函數(shù)Fig.3 Voltage CDF of node 4

    以系統(tǒng)接入風(fēng)電未考慮線路故障的情況為例,采用多目標(biāo)粒子群算法對(duì)系統(tǒng)進(jìn)行優(yōu)化,得到的Pareto最優(yōu)解集如圖4所示。從Pareto最優(yōu)解的分布來(lái)看,前沿分布相對(duì)均勻,具有良好的多樣性,證明了本文所設(shè)計(jì)程序的有效性。

    圖4 未考慮線路故障的Pareto最優(yōu)解集Fig.4 Pareto optimal solution set without considering line fault

    系統(tǒng)接入風(fēng)電后,針對(duì)不同的場(chǎng)景,得到優(yōu)化后各運(yùn)行風(fēng)險(xiǎn)和發(fā)電成本的結(jié)果如表4所示,對(duì)應(yīng)的各個(gè)發(fā)電機(jī)出力情況如表5所示。其中,優(yōu)化前不考慮線路故障,場(chǎng)景1為不考慮線路故障,場(chǎng)景2為考慮正常天氣下線路故障,場(chǎng)景3為考慮惡劣天氣下線路故障。

    表4 不同場(chǎng)景下系統(tǒng)的運(yùn)行風(fēng)險(xiǎn)和發(fā)電成本Tab.4 Operation risk and generation cost of the system under different scenarios

    表5 不同場(chǎng)景下各發(fā)電機(jī)出力情況Tab.5 Output of the generators under different scenarios

    從上述優(yōu)化結(jié)果中可以看出,優(yōu)化之前雖然系統(tǒng)的發(fā)電成本較低,但系統(tǒng)的運(yùn)行風(fēng)險(xiǎn)卻很高,在不考慮線路故障的情況下,對(duì)系統(tǒng)進(jìn)行發(fā)電優(yōu)化,雖然發(fā)電成本有所增加,但系統(tǒng)的運(yùn)行風(fēng)險(xiǎn)卻大大降低,優(yōu)化方案取得了很好的效果。而在考慮線路故障的情況下,相比不考慮線路故障,運(yùn)行風(fēng)險(xiǎn)的增幅明顯高于發(fā)電成本的增幅,主要是因?yàn)楫?dāng)這些線路發(fā)生故障時(shí),系統(tǒng)結(jié)構(gòu)發(fā)生改變,導(dǎo)致潮流轉(zhuǎn)移,與這些線路相鄰的線路和節(jié)點(diǎn)發(fā)生潮流越限和電壓越限的概率隨之增加,從而增大了系統(tǒng)的運(yùn)行風(fēng)險(xiǎn)。而對(duì)比不同天氣狀況下的優(yōu)化情況,可以看出系統(tǒng)在惡劣天氣下的風(fēng)險(xiǎn)要高于正常天氣下的,與實(shí)際相符,并且系統(tǒng)各個(gè)機(jī)組的出力狀況也有所不同,體現(xiàn)了不同天氣狀況下調(diào)度方案的差異性。

    5 結(jié)語(yǔ)

    電力系統(tǒng)的運(yùn)行環(huán)境日益嚴(yán)峻,忽略系統(tǒng)風(fēng)險(xiǎn)的電力系統(tǒng)優(yōu)化調(diào)度愈發(fā)顯得保守。而現(xiàn)有研究成果往往只注重個(gè)別風(fēng)險(xiǎn)因素對(duì)于系統(tǒng)運(yùn)行風(fēng)險(xiǎn)的影響,不能很好地反映系統(tǒng)整體的運(yùn)行風(fēng)險(xiǎn)。本文提出了一種計(jì)及運(yùn)行風(fēng)險(xiǎn)的優(yōu)化調(diào)度模型,以系統(tǒng)運(yùn)行風(fēng)險(xiǎn)和發(fā)電成本作為目標(biāo)函數(shù)。該模型將隨機(jī)潮流理論引入優(yōu)化調(diào)度,以其宏觀統(tǒng)計(jì)特性來(lái)處理風(fēng)電出力、負(fù)荷波動(dòng)、發(fā)電機(jī)強(qiáng)迫停運(yùn)、線路故障等一系列不確定因素的影響,并針對(duì)不同的天氣狀況,對(duì)線路的故障概率進(jìn)行實(shí)時(shí)修正從而更加接近實(shí)際運(yùn)行環(huán)境。求解模型時(shí)為體現(xiàn)不同目標(biāo)間的競(jìng)爭(zhēng)關(guān)系,采用多目標(biāo)粒子群算法及模糊決策理論進(jìn)行求解。實(shí)例證明所提模型及方法可很好地處理各種不確定因素的影響,并體現(xiàn)不同天氣狀況下調(diào)度方案的差異性,適合于實(shí)際的短期優(yōu)化調(diào)度。

    本文僅針對(duì)單個(gè)時(shí)段建立了優(yōu)化調(diào)度模型,且在評(píng)估系統(tǒng)越限后果時(shí),采用按一定百分比切除越限節(jié)點(diǎn)負(fù)荷和過(guò)載支路兩端節(jié)點(diǎn)負(fù)荷的方法。如何考慮不同時(shí)間斷面的耦合性,將其應(yīng)用到動(dòng)態(tài)優(yōu)化調(diào)度之中,并且更加準(zhǔn)確地評(píng)估系統(tǒng)越限后果,將是本文下一步研究工作的重點(diǎn)。

    參考文獻(xiàn):

    [1]胡澤春,丁華杰,孔濤(Hu Zechun,Ding Huajie,Kong Tao).風(fēng)電-抽水蓄能聯(lián)合日運(yùn)行優(yōu)化調(diào)度模型(A joint daily operational optimization model for wind power and pumped-storage plant)[J].電力系統(tǒng)自動(dòng)化(Automation ofElectricPowerSystems),2012,36(2):36-41,57.

    [2]Vlachogiannis J G,Lee K Y.Quantum-inspired evolution?ary algorithm for real and reactive power dispatch[J]. IEEE Trans on Power Systems,2008,23(4):1627-1636.

    [3]李學(xué)斌(Li Xuebin).機(jī)組負(fù)荷分配的多目標(biāo)優(yōu)化和多屬性決策(Multi-objective optimization and multi-attri?bute decision making analysis for economic load dispatch)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings of the CSUEPSA),2010,22(1):142-145.

    [4]Niknam T,Azizipanah-Abarghooee R,Zare M,et al.Re?serve constrained dynamic environmental/economic dis?patch:a new multiobjective self-adaptive learning bat al?gorithm[J].IEEE Systems Journal,2013,7(4):763-776.

    [5]陳碧云,韋杏秋,陳紹南,等(Chen Biyun,Wei Xingqiu,Chen Shaonan,et al).基于多種群遺傳算法的電力系統(tǒng)多目標(biāo)優(yōu)化(Power system multi-objective optimization based on multi-population genetic algorithm)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings of the CSU-EPSA),2015,27(7):24-29.

    [6]張理(Zhang Li).含風(fēng)電場(chǎng)的發(fā)輸電系統(tǒng)運(yùn)行風(fēng)險(xiǎn)評(píng)估(Operation Risk Assessment of Generation and Transmis?sion System Containing Wind Farms)[D].杭州:浙江大學(xué)電氣工程學(xué)院(Hangzhou:College of Electrical Engi?neering,Zhejiang University),2013.

    [7]路攀,李雪,王春亮(Lu Pan,Li Xue,Wang Chunliang).負(fù)荷相關(guān)的風(fēng)電并網(wǎng)系統(tǒng)線路過(guò)載風(fēng)險(xiǎn)控制(Risk control of transmission line overload for wind power inter?connected grid with correlated load)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings of the CSU-EPSA),2015,27 (11):16-20.

    [8]崔曉丹,方勇杰,周霞,等(Cui Xiaodan,F(xiàn)ang Yongjie,Zhou Xia,et al).基于風(fēng)險(xiǎn)的輸電斷面經(jīng)濟(jì)輸電功率計(jì)算(Risk-based economical transmission power calcula?tion of transmission sections)[J].電力系統(tǒng)自動(dòng)化(Auto?mation of Electric Power Systems),2014,38(8):30-34.

    [9]劉新東,江全元,曹一家,等(Liu Xindong,Jiang Quanyu?an,Cao Yijia,et al).基于風(fēng)險(xiǎn)理論和模糊推理的電力系統(tǒng)暫態(tài)安全風(fēng)險(xiǎn)評(píng)估(Transient security risk assess?ment of power system based on risk theory and fuzzy rea?soning)[J].電力自動(dòng)化設(shè)備(Electric Power Automation Equipment),2009,29(2):15-20.

    [10]宋曉喆,汪震,甘德強(qiáng),等(Song Xiaozhe,Wang Zhen,Gan Deqiang,et al).臺(tái)風(fēng)天氣條件下的電網(wǎng)暫態(tài)穩(wěn)定風(fēng)險(xiǎn)評(píng)估(Transient stability risk assessment of power grid under typhoon weather)[J].電力系統(tǒng)保護(hù)與控制(Power System Protection and Control),2012,40(24):1-8.

    [11]徐帆,王穎,楊建平,等(Xu Fan,Wang Ying,Yang Jian?ping,et al).考慮電網(wǎng)安全的風(fēng)電火電協(xié)調(diào)優(yōu)化調(diào)度模型及其求解(Generation scheduling model and applica?tion for wind-thermal power system considering security constraints)[J].電力系統(tǒng)自動(dòng)化(Automation of Electric Power Systems),2014,38(21):114-120.

    [12]陳道君,龔慶武,張茂林,等(Chen Daojun,Gong Qing?wu,Zhang Maolin,et al).考慮能源環(huán)境效益的含風(fēng)電場(chǎng)多目標(biāo)優(yōu)化調(diào)度(Multi-objective optimal dispatch in wind power integrated system incorporating energy-envi?ronmental efficiency)[J].中國(guó)電機(jī)工程學(xué)報(bào)(Proceed?ings of the CSEE),2011,31(13):10-17.

    [13]邱威,張建華,劉念,等(Qiu Wei,Zhang Jianhua,Liu Nian,et al).計(jì)及運(yùn)行風(fēng)險(xiǎn)的多目標(biāo)發(fā)電優(yōu)化調(diào)度(Multi-objective optimal generation dispatch with consid?eration of operation risk)[J].中國(guó)電機(jī)工程學(xué)報(bào)(Pro?ceedings of the CSEE),2012,32(22):64-72.

    [14]欒士巖(Luan Shiyan).計(jì)及風(fēng)力發(fā)電風(fēng)險(xiǎn)的電力系統(tǒng)優(yōu)化調(diào)度研究(Research on Optimal Scheduling of Pow?er System Considering Wind Power Risk)[D].上海:上海交通大學(xué)電子信息與電氣工程學(xué)院(Shanghai:School of Electronic Information and Electrical Engineering,Shanghai Jiao Tong University),2010.

    [15]姚瑤,于繼來(lái)(Yao Yao,Yu Jilai).計(jì)及風(fēng)電備用風(fēng)險(xiǎn)的電力系統(tǒng)多目標(biāo)混合優(yōu)化調(diào)度(Multi-objective hybrid optimal dispatch of power systems considering reserve risk due to wind power)[J].電力系統(tǒng)自動(dòng)化(Automation of Electric Power Systems),2011,35(22):118-124.

    [16]陳皓勇,胡斌,葉榮(Chen Haoyong,Hu Bin,Ye Rong).計(jì)及氣象災(zāi)害因素的電網(wǎng)規(guī)劃方法(Transmission net?work planning considering weather disaster)[J].電網(wǎng)技術(shù)(Power System Technology),2009,33(20):82-87.

    [17]張永熙(Zhang Yongxi).極端氣象災(zāi)害下考慮不確定斷線故障的電力系統(tǒng)隨機(jī)優(yōu)化調(diào)度(Stochastic Optimal Dispatch of Power System Considering Uncertain Line Failures under Extreme Weather Disaster)[D].長(zhǎng)沙:長(zhǎng)沙理工大學(xué)電氣與信息工程學(xué)院(Changsha:School of Electrical and Information Engineering,Changsha Uni?versity of Science&Technology),2011.

    [18]Q/CSG 11104002—2012,南方電網(wǎng)運(yùn)行安全風(fēng)險(xiǎn)量化評(píng)估技術(shù)規(guī)范[S].

    [19]劉怡芳,張步涵,李俊芳,等(Liu Yifang,Zhang Buhan,Li Junfang,et al).考慮電網(wǎng)靜態(tài)安全風(fēng)險(xiǎn)的隨機(jī)潮流計(jì)算(Probabilistic load flow algorithm considering static security risk of the power system)[J].中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE),2011,31(1):59-64.

    [20]高立志(Gao Lizhi).考慮風(fēng)電場(chǎng)間歇性的電力系統(tǒng)風(fēng)險(xiǎn)評(píng)估(Power System Risk Assessment Considering Wind Farm Intermittency)[D].北京:北京交通大學(xué)電氣工程學(xué)院(Beijing:School of Electrical Engineering,Beijing Jiaotong University),2012.

    [21]王紅印,張明亮,孫素琴,等(Wang Hongyin,Zhang Min?gliang,Sun Suqin,et al).基于風(fēng)險(xiǎn)定量分析的電壓安全評(píng)估方法(Research on voltage security evaluation meth?od based on risk quantitative analysis)[J].電網(wǎng)技術(shù)(Pow?er System Technology),2007,31(S2):98-102.

    [22]Coello C A C,Pulido G T,Lechuga M S.Handling multi?ple objectives with particle swarm optimization[J].IEEE Trans on Evolutionary Computation,2004,8(3):256-279.

    [23]彭春華,孫惠娟(Peng Chunhua,Sun Huijuan).基于非劣排序微分進(jìn)化的多目標(biāo)優(yōu)化發(fā)電調(diào)度(Multi-objective optimization power dispatch based on non-dominated sort?ing differential evolution)[J].中國(guó)電機(jī)工程學(xué)報(bào)(Pro?ceedings of the CSEE),2009,29(34):71-76.

    李燦(1990—),男,碩士研究生,研究方向?yàn)殡娏ο到y(tǒng)優(yōu)化運(yùn)行。Email:215638621@qq.com

    曾沅(1975—),男,博士,副教授,研究方向?yàn)殡娏ο到y(tǒng)安全性與穩(wěn)定性、電力系統(tǒng)規(guī)劃、電力系統(tǒng)優(yōu)化運(yùn)行。Email:zengyuan@tju.edu.cn

    秦超(1986—),男,博士,講師,研究方向?yàn)殡娏ο到y(tǒng)風(fēng)險(xiǎn)評(píng)估。Email:chaoqin@tju.edu.cn

    中圖分類號(hào):TM73

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1003-8930(2016)06-0073-07

    DOI:10.3969/j.issn.1003-8930.2016.06.013

    作者簡(jiǎn)介:

    收稿日期:2015-11-13;修回日期:2015-12-28

    Optimization Method for Power System Dispatching Considering Operation Risk

    LI Can1,ZENG Yuan1,QIN Chao1,SONG Yunting2,JI Ping2,WU Wei3
    (1.Key Laboratory of Smart Grid of Ministry of Education,Tianjin University,Tianjin 300072,China;2.China Elctric Power Research Institute,Beijing 100192,China;3.State Grid Fujian Electric Power Co.,Ltd.,Economic and Technology Institute,F(xiàn)uzhou 350012,China)

    Abstract:In view of the situation that the consideration of risk factors in the existing optimal dispatch methods is too sin?gle and one-sided,the entire operation risk of the system is measured by branch overload and node voltage violation of limit,and a model based on operation risk is formulated with operation risk and generation cost as objective functions,allowing for the influence of wind power fluctuation,unit outage,line fault,load fluctuation and other uncertain factors. In this model,the stochastic load flow theory is introduced into the optimal dispatch to deal with the effect of many un?certainties on the operation risk with its Macro statistical properties,and for different weather conditions,the line fault probability is under real-time correction to get closer to the actual operating environment.This model is solved by the multi-objective particle swarm algorithm and the fuzzy set theory to seek a balance between economy and security.The IEEE-14 system is used to test the model and the results show that it is able to deal with the impact of uncertainties on the operation risk and reflect the difference of the dispatch scheme under different weather conditions which can provide a reference for the practical short-term optimal dispatch.

    Key words:optimal dispatch;operation risk;uncertainties;stochastic load flow

    视频区图区小说| 精品国产超薄肉色丝袜足j| 久久久久久久久久人人人人人人| 男女边吃奶边做爰视频| 秋霞伦理黄片| 啦啦啦在线免费观看视频4| 午夜福利在线观看免费完整高清在| 亚洲精品视频女| 国产一区有黄有色的免费视频| av国产久精品久网站免费入址| 国产视频首页在线观看| 欧美精品一区二区大全| 人人妻人人爽人人添夜夜欢视频| 熟妇人妻不卡中文字幕| 伊人久久大香线蕉亚洲五| 韩国精品一区二区三区| 亚洲人成网站在线观看播放| 美女高潮到喷水免费观看| 免费观看在线日韩| 久久久久久久久久久久大奶| 午夜福利网站1000一区二区三区| 街头女战士在线观看网站| av福利片在线| 国产精品亚洲av一区麻豆 | 捣出白浆h1v1| 国产成人一区二区在线| 男人添女人高潮全过程视频| 国产免费又黄又爽又色| 成人国产麻豆网| 亚洲国产精品国产精品| 男人舔女人的私密视频| 97人妻天天添夜夜摸| 丰满迷人的少妇在线观看| 黄片无遮挡物在线观看| 交换朋友夫妻互换小说| h视频一区二区三区| 男人舔女人的私密视频| 青春草国产在线视频| av一本久久久久| 一本—道久久a久久精品蜜桃钙片| 五月伊人婷婷丁香| 日韩不卡一区二区三区视频在线| av网站免费在线观看视频| 90打野战视频偷拍视频| 黄色配什么色好看| kizo精华| 人人妻人人爽人人添夜夜欢视频| 成人免费观看视频高清| 亚洲,欧美,日韩| 久久国产精品大桥未久av| 久久精品人人爽人人爽视色| 这个男人来自地球电影免费观看 | 婷婷色麻豆天堂久久| 精品久久久久久电影网| 亚洲综合色网址| 制服人妻中文乱码| 青青草视频在线视频观看| www.熟女人妻精品国产| 91久久精品国产一区二区三区| 精品一区在线观看国产| 午夜久久久在线观看| 国产亚洲午夜精品一区二区久久| 黑人猛操日本美女一级片| 男女高潮啪啪啪动态图| 精品一区在线观看国产| 久久综合国产亚洲精品| 成人影院久久| 黄片播放在线免费| 欧美 亚洲 国产 日韩一| 国产精品久久久久久久久免| 亚洲欧美色中文字幕在线| 亚洲在久久综合| 欧美人与善性xxx| 少妇人妻精品综合一区二区| 国产在线一区二区三区精| 国产伦理片在线播放av一区| 免费女性裸体啪啪无遮挡网站| 国产精品一国产av| 亚洲国产av新网站| 多毛熟女@视频| 久久 成人 亚洲| 一级a爱视频在线免费观看| 一边亲一边摸免费视频| 色哟哟·www| av免费在线看不卡| 啦啦啦在线观看免费高清www| 午夜免费鲁丝| 国产 精品1| 在线亚洲精品国产二区图片欧美| 久久国产精品男人的天堂亚洲| 国产白丝娇喘喷水9色精品| 汤姆久久久久久久影院中文字幕| 中国三级夫妇交换| 国产免费视频播放在线视频| 国产精品国产三级国产专区5o| 一个人免费看片子| www.av在线官网国产| 日韩av在线免费看完整版不卡| 久久99蜜桃精品久久| 亚洲欧洲日产国产| 亚洲国产欧美网| 中文字幕亚洲精品专区| 国产又色又爽无遮挡免| 在线天堂最新版资源| 永久网站在线| 久久毛片免费看一区二区三区| 久久狼人影院| 色网站视频免费| 久久免费观看电影| 色婷婷av一区二区三区视频| 天堂8中文在线网| 欧美 亚洲 国产 日韩一| 一区二区三区激情视频| 亚洲欧美色中文字幕在线| 久久久国产精品麻豆| 日韩中字成人| 国产精品人妻久久久影院| 男女国产视频网站| 午夜福利视频精品| 一级a爱视频在线免费观看| 国产深夜福利视频在线观看| 美女中出高潮动态图| a级片在线免费高清观看视频| 最近中文字幕2019免费版| 丝袜人妻中文字幕| 久久久久久伊人网av| 亚洲 欧美一区二区三区| 亚洲国产色片| 久久人人97超碰香蕉20202| 精品一品国产午夜福利视频| 免费大片黄手机在线观看| 久久久久久久久久久免费av| 成人亚洲精品一区在线观看| 人人妻人人爽人人添夜夜欢视频| 精品人妻在线不人妻| 国产深夜福利视频在线观看| 秋霞在线观看毛片| 亚洲一区中文字幕在线| 亚洲色图综合在线观看| 欧美国产精品一级二级三级| 国产成人免费观看mmmm| 97在线视频观看| 卡戴珊不雅视频在线播放| 久久久久网色| 侵犯人妻中文字幕一二三四区| 久久精品国产鲁丝片午夜精品| 亚洲人成网站在线观看播放| 亚洲欧美中文字幕日韩二区| videosex国产| 免费看不卡的av| 日本欧美视频一区| 赤兔流量卡办理| 国产男女超爽视频在线观看| 视频在线观看一区二区三区| 久久精品亚洲av国产电影网| 欧美日韩精品成人综合77777| 日韩大片免费观看网站| 久久综合国产亚洲精品| 一本大道久久a久久精品| 9191精品国产免费久久| 免费看av在线观看网站| 久久人人97超碰香蕉20202| 天天躁狠狠躁夜夜躁狠狠躁| 男人添女人高潮全过程视频| h视频一区二区三区| 波多野结衣av一区二区av| 日本91视频免费播放| 日韩制服丝袜自拍偷拍| 女人被躁到高潮嗷嗷叫费观| 日韩一区二区三区影片| 一级a爱视频在线免费观看| 一级a爱视频在线免费观看| 亚洲国产色片| 9191精品国产免费久久| 国精品久久久久久国模美| 久久精品久久精品一区二区三区| 黄片播放在线免费| 国产精品一区二区在线不卡| 成人影院久久| 国产精品 欧美亚洲| 日韩大片免费观看网站| 精品一区二区三区四区五区乱码 | 国产爽快片一区二区三区| 丰满乱子伦码专区| 亚洲av综合色区一区| 少妇的逼水好多| 亚洲精品日本国产第一区| 欧美精品国产亚洲| 永久网站在线| 一区二区日韩欧美中文字幕| 国产精品久久久久成人av| 午夜免费鲁丝| 一二三四中文在线观看免费高清| 啦啦啦中文免费视频观看日本| 日韩成人av中文字幕在线观看| 日韩成人av中文字幕在线观看| 亚洲av电影在线观看一区二区三区| 日韩av免费高清视频| 亚洲国产av新网站| 国产一区二区在线观看av| 久久久精品国产亚洲av高清涩受| 精品一区二区三区四区五区乱码 | 中文字幕人妻丝袜制服| 国产男女超爽视频在线观看| 亚洲人成77777在线视频| 欧美日韩综合久久久久久| 大陆偷拍与自拍| 成人免费观看视频高清| 欧美激情 高清一区二区三区| 亚洲第一青青草原| 亚洲第一青青草原| 国产毛片在线视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产免费又黄又爽又色| 日韩欧美一区视频在线观看| 人妻一区二区av| 亚洲久久久国产精品| 日韩,欧美,国产一区二区三区| 亚洲,欧美,日韩| 国产一区二区激情短视频 | 超碰97精品在线观看| 十八禁高潮呻吟视频| 成人黄色视频免费在线看| 婷婷色综合大香蕉| 日韩中字成人| 亚洲伊人久久精品综合| 高清av免费在线| 一级毛片电影观看| 国产精品一区二区在线不卡| 成人影院久久| 啦啦啦在线免费观看视频4| 毛片一级片免费看久久久久| 中文字幕av电影在线播放| 丰满迷人的少妇在线观看| 午夜免费鲁丝| 精品少妇内射三级| 春色校园在线视频观看| 老司机影院成人| 丝袜在线中文字幕| 国产乱人偷精品视频| 午夜免费观看性视频| 亚洲在久久综合| 欧美日韩综合久久久久久| 国产精品 国内视频| 波多野结衣一区麻豆| 一本色道久久久久久精品综合| 亚洲综合色网址| 9色porny在线观看| 免费观看无遮挡的男女| 超碰成人久久| 啦啦啦啦在线视频资源| 精品国产超薄肉色丝袜足j| 国产乱人偷精品视频| 欧美变态另类bdsm刘玥| 欧美+日韩+精品| 国产免费现黄频在线看| 看十八女毛片水多多多| 高清黄色对白视频在线免费看| 国产精品国产三级专区第一集| 26uuu在线亚洲综合色| 最近最新中文字幕免费大全7| av一本久久久久| 热re99久久国产66热| 丰满少妇做爰视频| 夜夜骑夜夜射夜夜干| 欧美成人午夜精品| 久久久久久人妻| 国产免费又黄又爽又色| 亚洲男人天堂网一区| 日本免费在线观看一区| 丝袜人妻中文字幕| 亚洲伊人久久精品综合| 免费看av在线观看网站| 妹子高潮喷水视频| www.自偷自拍.com| 久热这里只有精品99| 欧美激情极品国产一区二区三区| 国产免费福利视频在线观看| 免费播放大片免费观看视频在线观看| 一本—道久久a久久精品蜜桃钙片| 最新中文字幕久久久久| 久久精品国产综合久久久| 精品一品国产午夜福利视频| 日本av免费视频播放| 色视频在线一区二区三区| 99久久综合免费| 天美传媒精品一区二区| 亚洲精品国产av蜜桃| 成人免费观看视频高清| 中文精品一卡2卡3卡4更新| 人妻少妇偷人精品九色| 99国产精品免费福利视频| 人成视频在线观看免费观看| 午夜福利,免费看| 天天躁夜夜躁狠狠久久av| 国产亚洲欧美精品永久| 好男人视频免费观看在线| 久久人人爽人人片av| 精品人妻一区二区三区麻豆| 国产午夜精品一二区理论片| 色94色欧美一区二区| 亚洲精品视频女| 亚洲一级一片aⅴ在线观看| 久久热在线av| 中文字幕精品免费在线观看视频| 激情视频va一区二区三区| 国产亚洲精品第一综合不卡| 久久亚洲国产成人精品v| 制服丝袜香蕉在线| 成人黄色视频免费在线看| 国产男人的电影天堂91| 精品少妇久久久久久888优播| 国产在视频线精品| 激情视频va一区二区三区| 久久精品国产鲁丝片午夜精品| 国产日韩欧美亚洲二区| 国产成人一区二区在线| 丝袜美腿诱惑在线| 国产一区二区三区av在线| 十分钟在线观看高清视频www| 国产精品 欧美亚洲| 精品少妇一区二区三区视频日本电影 | 精品国产一区二区三区四区第35| 中文字幕人妻熟女乱码| 在线 av 中文字幕| 尾随美女入室| xxx大片免费视频| 黄色 视频免费看| 女人久久www免费人成看片| 日韩精品有码人妻一区| 久久久久久久精品精品| 色婷婷久久久亚洲欧美| 如何舔出高潮| 亚洲男人天堂网一区| 男的添女的下面高潮视频| 少妇被粗大猛烈的视频| 男女国产视频网站| 99久国产av精品国产电影| 国产亚洲精品第一综合不卡| 色94色欧美一区二区| 亚洲国产av新网站| 精品久久久精品久久久| 亚洲精品国产av蜜桃| 美女中出高潮动态图| 久久精品人人爽人人爽视色| 国产免费福利视频在线观看| 最新中文字幕久久久久| 成人18禁高潮啪啪吃奶动态图| 国产深夜福利视频在线观看| 国产亚洲午夜精品一区二区久久| 天天躁日日躁夜夜躁夜夜| 老汉色∧v一级毛片| 人人妻人人爽人人添夜夜欢视频| 亚洲激情五月婷婷啪啪| 精品久久蜜臀av无| 午夜福利影视在线免费观看| av线在线观看网站| 久久鲁丝午夜福利片| 久久99蜜桃精品久久| www.熟女人妻精品国产| 国语对白做爰xxxⅹ性视频网站| 日日啪夜夜爽| 久久久久精品人妻al黑| 国产 一区精品| 亚洲欧美中文字幕日韩二区| 巨乳人妻的诱惑在线观看| 狠狠婷婷综合久久久久久88av| 如日韩欧美国产精品一区二区三区| 国产人伦9x9x在线观看 | h视频一区二区三区| 欧美xxⅹ黑人| 18禁裸乳无遮挡动漫免费视频| 亚洲成人手机| 久久精品国产a三级三级三级| 久久久久久久久久久免费av| 中文字幕人妻熟女乱码| av国产精品久久久久影院| 黄色视频在线播放观看不卡| 亚洲精品久久成人aⅴ小说| 日韩av不卡免费在线播放| 多毛熟女@视频| 超碰成人久久| h视频一区二区三区| 中文天堂在线官网| 欧美日韩亚洲国产一区二区在线观看 | 免费观看无遮挡的男女| 亚洲国产精品成人久久小说| 在线亚洲精品国产二区图片欧美| 2021少妇久久久久久久久久久| 高清在线视频一区二区三区| 99re6热这里在线精品视频| 久久99蜜桃精品久久| 最近中文字幕高清免费大全6| 欧美成人精品欧美一级黄| 午夜福利在线观看免费完整高清在| 日韩伦理黄色片| 午夜日韩欧美国产| 国产精品久久久av美女十八| 高清黄色对白视频在线免费看| 一本色道久久久久久精品综合| 国产精品久久久久久精品电影小说| 欧美黄色片欧美黄色片| 国产一区二区在线观看av| 又黄又粗又硬又大视频| 成年女人在线观看亚洲视频| 欧美日韩视频精品一区| 国语对白做爰xxxⅹ性视频网站| 精品视频人人做人人爽| 午夜av观看不卡| 一本久久精品| a级片在线免费高清观看视频| 精品国产乱码久久久久久小说| 性色avwww在线观看| av视频免费观看在线观看| 免费在线观看黄色视频的| 在线 av 中文字幕| 成年人免费黄色播放视频| 搡女人真爽免费视频火全软件| 中文乱码字字幕精品一区二区三区| 考比视频在线观看| 国产熟女午夜一区二区三区| 咕卡用的链子| 宅男免费午夜| 日韩大片免费观看网站| 国产乱来视频区| 欧美激情极品国产一区二区三区| 久久久国产一区二区| 国产一级毛片在线| 午夜福利在线观看免费完整高清在| 综合色丁香网| 国产福利在线免费观看视频| 久久久精品94久久精品| 成年av动漫网址| 婷婷色麻豆天堂久久| 久久午夜福利片| 午夜福利一区二区在线看| 观看av在线不卡| 久久久久网色| 在线天堂中文资源库| 熟女少妇亚洲综合色aaa.| 伊人久久大香线蕉亚洲五| 男女啪啪激烈高潮av片| 色播在线永久视频| 18+在线观看网站| 精品一区二区免费观看| 最近的中文字幕免费完整| 人体艺术视频欧美日本| 精品第一国产精品| 国产精品麻豆人妻色哟哟久久| 69精品国产乱码久久久| 欧美日本中文国产一区发布| 一区二区日韩欧美中文字幕| 视频区图区小说| 亚洲一区中文字幕在线| 久久女婷五月综合色啪小说| av一本久久久久| 热re99久久国产66热| 日韩欧美一区视频在线观看| 精品一区二区三区四区五区乱码 | 老汉色av国产亚洲站长工具| 亚洲欧洲日产国产| 精品亚洲成a人片在线观看| 一区在线观看完整版| 人妻系列 视频| 美女大奶头黄色视频| 国产在视频线精品| 岛国毛片在线播放| 天美传媒精品一区二区| 少妇被粗大的猛进出69影院| 啦啦啦在线观看免费高清www| 久久精品国产鲁丝片午夜精品| 日本av免费视频播放| 侵犯人妻中文字幕一二三四区| 最近最新中文字幕大全免费视频 | 日韩中字成人| 日韩中文字幕欧美一区二区 | 老司机影院毛片| 中文字幕精品免费在线观看视频| 久久免费观看电影| 美女视频免费永久观看网站| 亚洲综合色惰| 久久这里有精品视频免费| 亚洲国产精品一区三区| 老司机影院毛片| a 毛片基地| 91国产中文字幕| 久久久久久人妻| 亚洲国产欧美在线一区| 熟女av电影| 久久 成人 亚洲| 肉色欧美久久久久久久蜜桃| 狂野欧美激情性bbbbbb| 婷婷色综合大香蕉| 丝袜喷水一区| 精品一区在线观看国产| 亚洲在久久综合| 9191精品国产免费久久| 成人毛片a级毛片在线播放| 欧美人与性动交α欧美软件| 欧美激情高清一区二区三区 | 在线观看免费视频网站a站| 我的亚洲天堂| 国产熟女欧美一区二区| 最近的中文字幕免费完整| 精品国产乱码久久久久久小说| 天天影视国产精品| 精品第一国产精品| 国产熟女欧美一区二区| 建设人人有责人人尽责人人享有的| 免费播放大片免费观看视频在线观看| 巨乳人妻的诱惑在线观看| 人人澡人人妻人| 国产av码专区亚洲av| 国产精品免费视频内射| 久久精品久久久久久噜噜老黄| 国产精品麻豆人妻色哟哟久久| 亚洲国产精品一区三区| 人体艺术视频欧美日本| 麻豆精品久久久久久蜜桃| 亚洲美女视频黄频| 考比视频在线观看| 99久久综合免费| 精品福利永久在线观看| 亚洲国产欧美日韩在线播放| 中文字幕色久视频| 国产成人一区二区在线| 亚洲综合精品二区| 青草久久国产| 新久久久久国产一级毛片| 激情五月婷婷亚洲| 日韩熟女老妇一区二区性免费视频| 国产精品久久久久久精品电影小说| 欧美人与性动交α欧美精品济南到 | 精品少妇黑人巨大在线播放| 国产在视频线精品| 日本av免费视频播放| 亚洲欧美色中文字幕在线| 成人漫画全彩无遮挡| 久久女婷五月综合色啪小说| 韩国av在线不卡| 亚洲国产av影院在线观看| 一区二区日韩欧美中文字幕| 精品少妇内射三级| 2022亚洲国产成人精品| 亚洲国产精品国产精品| 91久久精品国产一区二区三区| 久久久国产精品麻豆| a 毛片基地| 一区二区av电影网| 在线观看免费视频网站a站| 久久韩国三级中文字幕| 午夜日本视频在线| 乱人伦中国视频| 一级毛片电影观看| 女人久久www免费人成看片| 国产精品久久久久久精品电影小说| 国产极品粉嫩免费观看在线| av在线app专区| 国产亚洲欧美精品永久| 欧美日本中文国产一区发布| 下体分泌物呈黄色| 又大又黄又爽视频免费| 一区二区日韩欧美中文字幕| 王馨瑶露胸无遮挡在线观看| 男女边摸边吃奶| 国产深夜福利视频在线观看| 久久影院123| 国产熟女欧美一区二区| 999精品在线视频| 在线观看三级黄色| 五月开心婷婷网| 日韩 亚洲 欧美在线| 日韩不卡一区二区三区视频在线| 亚洲视频免费观看视频| 丰满少妇做爰视频| 大香蕉久久成人网| 欧美精品高潮呻吟av久久| 综合色丁香网| 天天操日日干夜夜撸| 狠狠精品人妻久久久久久综合| 欧美中文综合在线视频| 色94色欧美一区二区| 国产精品一二三区在线看| 一本—道久久a久久精品蜜桃钙片| 老汉色av国产亚洲站长工具| 国产成人aa在线观看| 国产成人免费观看mmmm| 国产在线视频一区二区| 久久久久精品性色| 亚洲av成人精品一二三区| 国产精品国产av在线观看| 亚洲国产最新在线播放| 久久99蜜桃精品久久| 久久狼人影院| 午夜91福利影院| 青草久久国产| 一区福利在线观看| 色播在线永久视频| 香蕉国产在线看| 巨乳人妻的诱惑在线观看| av卡一久久| 超色免费av| 日韩电影二区| 男男h啪啪无遮挡| 欧美精品人与动牲交sv欧美| 国产欧美日韩一区二区三区在线| 国产人伦9x9x在线观看 | 天堂俺去俺来也www色官网| 欧美精品国产亚洲| 久久毛片免费看一区二区三区| 免费在线观看视频国产中文字幕亚洲 | av国产久精品久网站免费入址| 免费人妻精品一区二区三区视频| 午夜影院在线不卡| 免费观看无遮挡的男女|