• <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.中國電力科學(xué)研究院,北京 100192;3.國網(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.中國電力科學(xué)研究院,北京 100192;3.國網(wǎng)福建省電力有限公司經(jīng)濟(jì)技術(shù)研究院,福州 350012)

    摘 要:針對現(xiàn)有優(yōu)化調(diào)度方法在考慮風(fēng)險(xiǎn)因素時(shí)過于單一、片面的現(xiàn)狀,該文考慮風(fēng)電波動、機(jī)組停運(yùn)、線路故障、負(fù)荷波動等多種不確定因素的影響,以電壓越限和潮流越限衡量系統(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)度中,并針對不同的天氣狀況,對線路故障概率進(jìn)行實(shí)時(shí)修正以期更加接近實(shí)際運(yùn)行環(huán)境。結(jié)合多目標(biāo)粒子群算法及模糊集理論對所建立的模型進(jìn)行求解,算例結(jié)果表明所提模型和方法能夠處理不確定因素對系統(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ī)潮流

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

    文獻(xiàn)[3-5]以節(jié)能、經(jīng)濟(jì)、環(huán)保三大原則為目標(biāo),建立了多目標(biāo)優(yōu)化調(diào)度模型。但這些方法均未對系統(tǒng)的安全可靠性給予太多的重視。目前有關(guān)電網(wǎng)風(fēng)險(xiǎn)的研究主要圍繞靜態(tài)安全評估[6-7]、輸電斷面輸電能力分析[8]、暫態(tài)安全評估[9-10]等方面,取得了大量成果。而近年來在優(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ī)出力對系統(tǒng)風(fēng)險(xiǎn)的影響,而忽略了影響系統(tǒng)風(fēng)險(xiǎn)的其他因素。文獻(xiàn)[12]考慮發(fā)電資源消耗最少、能源環(huán)境效益最好、系統(tǒng)安全穩(wěn)定程度最高等因素,提出了含風(fē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),針對線路故障,構(gòu)建了基于風(fēng)險(xiǎn)的多目標(biāo)經(jīng)濟(jì)調(diào)度模型,但未考慮風(fēng)電出力的波動性對于系統(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)害對輸電線路的影響,建立了考慮天氣因素的機(jī)組優(yōu)化方法。綜上,現(xiàn)有的優(yōu)化調(diào)度方法在考慮風(fēng)險(xiǎn)因素時(shí)大多只考慮了個(gè)別不確定因素的影響。

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

    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í)段動態(tài)優(yōu)化調(diào)度。多時(shí)段動態(tài)優(yōu)化調(diào)度是在單時(shí)段靜態(tài)優(yōu)化調(diào)度的基礎(chǔ)上考慮了不同時(shí)間斷面的耦合性。為簡化計(jì)算,驗(yàn)證本文所提方法的可行性,文中僅針對單時(shí)段建立了優(yōu)化調(diào)度模型。

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

    電力系統(tǒng)中節(jié)點(diǎn)電壓和支路潮流是否處于合理的區(qū)間內(nèi),對電力系統(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ī)有功出力和無功出力;Pw、Qw分別為風(fēng)機(jī)的有功出力和無功出力;Pdi、Qdi分別為節(jié)點(diǎn)i的有功負(fù)荷和無功負(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有功和無功出力的上下限。

    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ù)測誤差對正旋轉(zhuǎn)備用的需求;Pwmax為風(fēng)機(jī)的最大出力;Pw為風(fēng)機(jī)的有功出力;wu%、wd%分別為風(fēng)機(jī)出力預(yù)測誤差對正、負(fù)旋轉(zhuǎn)備用的需求。

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

    隨機(jī)潮流可以考慮多種不確定因素的影響,以其數(shù)學(xué)模型為基礎(chǔ),通過概率理論得到系統(tǒng)運(yùn)行狀況的宏觀統(tǒng)計(jì)信息。相比于傳統(tǒng)的確定性潮流分析,能夠以較小的計(jì)算量更加全面、深刻地展現(xiàn)系統(tǒng)的運(yùn)行特性。因此本節(jié)采用隨機(jī)潮流處理風(fēng)電波動、負(fù)荷波動、機(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)行,因此其無功出力的概率密度函數(shù)可表示為

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

    2)負(fù)荷的波動性

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

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

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

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

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

    4)線路故障

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

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

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

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

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

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

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

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

    2.2隨機(jī)潮流

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

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

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

    式中:ΔW、ΔX、ΔZ表示注入功率的隨機(jī)擾動量以及對應(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級數(shù)相結(jié)合,對隨機(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í),可能會造成相應(yīng)的保護(hù)裝置動作,導(dǎo)致在保護(hù)裝置動作期間系統(tǒng)出現(xiàn)一定程度的負(fù)荷損失。在電力市場環(huán)境下,可以采用貨幣的形式來衡量系統(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為因供電合同中斷而帶來的懲罰成本。

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

    式中:Uim

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

    pmlia

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

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

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

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

    將當(dāng)前找到的非支配解按以下規(guī)則保存在外部檔案中:當(dāng)外部檔案個(gè)體數(shù)小于設(shè)定值時(shí),非支配解直接加入外部檔案;當(dāng)外部檔案個(gè)體數(shù)達(dá)到設(shè)定值時(shí),如果新解支配了外部檔案的部分個(gè)體,則新解取代受支配的那些個(gè)體,否則根據(jù)自適應(yīng)網(wǎng)格法[22]移除一個(gè)個(gè)體,將新解加入檔案中。隨著迭代過程不斷維護(hù)外部檔案,最后,根據(jù)一定的策略從外部檔案中選取一個(gè)折中解作為優(yōu)化問題的最優(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]。具體做法如下:

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

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

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

    4 算例分析

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

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

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

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

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

    在優(yōu)化前,分別針對未接入風(fēng)電、接入風(fēng)電、接入風(fēng)電并考慮線路故障3種情況對系統(tǒng)進(jìn)行隨機(jī)潮流分析,以節(jié)點(diǎn)4的電壓為例,其概率密度函數(shù)(PDF)和累積分布函數(shù)(CDF)如圖2和圖3所示??梢钥闯?,風(fē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)粒子群算法對系統(tǒng)進(jìn)行優(yōu)化,得到的Pareto最優(yōu)解集如圖4所示。從Pareto最優(yōu)解的分布來看,前沿分布相對均勻,具有良好的多樣性,證明了本文所設(shè)計(jì)程序的有效性。

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

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

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

    表5 不同場景下各發(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)卻很高,在不考慮線路故障的情況下,對系統(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)。而對比不同天氣狀況下的優(yōu)化情況,可以看出系統(tǒng)在惡劣天氣下的風(fēng)險(xiǎn)要高于正常天氣下的,與實(shí)際相符,并且系統(tǒng)各個(gè)機(jī)組的出力狀況也有所不同,體現(xiàn)了不同天氣狀況下調(diào)度方案的差異性。

    5 結(jié)語

    電力系統(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)因素對于系統(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ì)特性來處理風(fēng)電出力、負(fù)荷波動、發(fā)電機(jī)強(qiáng)迫停運(yùn)、線路故障等一系列不確定因素的影響,并針對不同的天氣狀況,對線路的故障概率進(jìn)行實(shí)時(shí)修正從而更加接近實(shí)際運(yùn)行環(huán)境。求解模型時(shí)為體現(xiàn)不同目標(biāo)間的競爭關(guān)系,采用多目標(biāo)粒子群算法及模糊決策理論進(jìn)行求解。實(shí)例證明所提模型及方法可很好地處理各種不確定因素的影響,并體現(xiàn)不同天氣狀況下調(diào)度方案的差異性,適合于實(shí)際的短期優(yōu)化調(diào)度。

    本文僅針對單個(gè)時(shí)段建立了優(yōu)化調(diào)度模型,且在評估系統(tǒng)越限后果時(shí),采用按一定百分比切除越限節(jié)點(diǎn)負(fù)荷和過載支路兩端節(jié)點(diǎn)負(fù)荷的方法。如何考慮不同時(shí)間斷面的耦合性,將其應(yīng)用到動態(tài)優(yōu)化調(diào)度之中,并且更加準(zhǔn)確地評估系統(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)自動化(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)及其自動化學(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)及其自動化學(xué)報(bào)(Proceedings of the CSU-EPSA),2015,27(7):24-29.

    [6]張理(Zhang Li).含風(fēng)電場的發(fā)輸電系統(tǒng)運(yùn)行風(fēng)險(xiǎn)評估(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)線路過載風(fēng)險(xiǎn)控制(Risk control of transmission line overload for wind power inter?connected grid with correlated load)[J].電力系統(tǒ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)自動化(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)評估(Transient security risk assess?ment of power system based on risk theory and fuzzy rea?soning)[J].電力自動化設(shè)備(Electric Power Automation Equipment),2009,29(2):15-20.

    [10]宋曉喆,汪震,甘德強(qiáng),等(Song Xiaozhe,Wang Zhen,Gan Deqiang,et al).臺風(fēng)天氣條件下的電網(wǎng)暫態(tài)穩(wěn)定風(fēng)險(xiǎn)評估(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)自動化(Automation of Electric Power Systems),2014,38(21):114-120.

    [12]陳道君,龔慶武,張茂林,等(Chen Daojun,Gong Qing?wu,Zhang Maolin,et al).考慮能源環(huán)境效益的含風(fēng)電場多目標(biāo)優(yōu)化調(diào)度(Multi-objective optimal dispatch in wind power integrated system incorporating energy-envi?ronmental efficiency)[J].中國電機(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].中國電機(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]姚瑤,于繼來(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)自動化(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].長沙:長沙理工大學(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)量化評估技術(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].中國電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE),2011,31(1):59-64.

    [20]高立志(Gao Lizhi).考慮風(fēng)電場間歇性的電力系統(tǒng)風(fēng)險(xiǎn)評估(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)定量分析的電壓安全評估方法(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].中國電機(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)評估。Email:chaoqin@tju.edu.cn

    中圖分類號:TM73

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

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

    DOI:10.3969/j.issn.1003-8930.2016.06.013

    作者簡介:

    收稿日期: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

    69av精品久久久久久| 欧美变态另类bdsm刘玥| 美女脱内裤让男人舔精品视频| 蜜臀久久99精品久久宅男| av在线蜜桃| 国产极品天堂在线| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美高清成人免费视频www| 成年版毛片免费区| 日韩人妻高清精品专区| 最新中文字幕久久久久| 校园人妻丝袜中文字幕| 熟女人妻精品中文字幕| 国产大屁股一区二区在线视频| 国产淫语在线视频| 亚洲av日韩在线播放| 国产一区二区在线观看日韩| 一级a做视频免费观看| av卡一久久| 麻豆成人午夜福利视频| 国产综合懂色| 女人久久www免费人成看片| 国内少妇人妻偷人精品xxx网站| 男女下面进入的视频免费午夜| 久久久精品94久久精品| 男人和女人高潮做爰伦理| 成人亚洲精品av一区二区| 午夜视频国产福利| 交换朋友夫妻互换小说| 欧美变态另类bdsm刘玥| 亚洲精品自拍成人| 免费电影在线观看免费观看| 性插视频无遮挡在线免费观看| 免费在线观看成人毛片| 又爽又黄a免费视频| 伦精品一区二区三区| 1000部很黄的大片| 欧美日韩亚洲高清精品| 观看美女的网站| 免费观看性生交大片5| 最后的刺客免费高清国语| 国产精品一及| 亚洲精品久久午夜乱码| 色婷婷久久久亚洲欧美| 特大巨黑吊av在线直播| 欧美激情国产日韩精品一区| 亚洲av中文av极速乱| 97人妻精品一区二区三区麻豆| 精品久久久久久久久av| 久久久色成人| 久久精品熟女亚洲av麻豆精品| 三级经典国产精品| 3wmmmm亚洲av在线观看| 免费观看性生交大片5| 成年免费大片在线观看| 国产伦理片在线播放av一区| 国产黄频视频在线观看| 中国国产av一级| 久久国内精品自在自线图片| 国产一区有黄有色的免费视频| 国产黄片美女视频| 亚洲怡红院男人天堂| 最近中文字幕高清免费大全6| 水蜜桃什么品种好| 欧美变态另类bdsm刘玥| 国产高潮美女av| 国产精品国产三级国产av玫瑰| 男男h啪啪无遮挡| 免费看av在线观看网站| 国产伦在线观看视频一区| 国产成人精品久久久久久| 99热网站在线观看| 婷婷色综合www| 亚洲国产精品国产精品| 欧美最新免费一区二区三区| 国产精品av视频在线免费观看| 校园人妻丝袜中文字幕| 人妻系列 视频| 人妻 亚洲 视频| 纵有疾风起免费观看全集完整版| 久久人人爽人人爽人人片va| 国产精品99久久99久久久不卡 | av卡一久久| 欧美日韩亚洲高清精品| 久久久久久久久久人人人人人人| 国产成人免费观看mmmm| 亚洲va在线va天堂va国产| 蜜桃久久精品国产亚洲av| 久久鲁丝午夜福利片| 91aial.com中文字幕在线观看| 精品久久国产蜜桃| 99热全是精品| 久久久久国产网址| 亚洲欧美成人综合另类久久久| 亚洲,一卡二卡三卡| 一级av片app| 国产成人a区在线观看| 内地一区二区视频在线| 国产黄a三级三级三级人| 日日摸夜夜添夜夜添av毛片| 18+在线观看网站| 国产精品一区二区三区四区免费观看| 免费看不卡的av| 人人妻人人看人人澡| 久久热精品热| 中国国产av一级| 日本黄色片子视频| 国产成人freesex在线| 亚洲真实伦在线观看| 国产高清国产精品国产三级 | 亚洲国产成人一精品久久久| 三级男女做爰猛烈吃奶摸视频| 久久99热这里只有精品18| av国产精品久久久久影院| 97精品久久久久久久久久精品| 美女被艹到高潮喷水动态| 成人毛片a级毛片在线播放| 日韩三级伦理在线观看| 亚洲欧美中文字幕日韩二区| tube8黄色片| 国产探花极品一区二区| 看黄色毛片网站| 久久久成人免费电影| 纵有疾风起免费观看全集完整版| kizo精华| 人人妻人人爽人人添夜夜欢视频 | 少妇的逼水好多| 国产成人精品福利久久| 久久99热这里只有精品18| 性色avwww在线观看| 又爽又黄a免费视频| 美女脱内裤让男人舔精品视频| 99热网站在线观看| 在线免费十八禁| 免费少妇av软件| 午夜福利在线观看免费完整高清在| 成人国产av品久久久| 日韩国内少妇激情av| 中文欧美无线码| 少妇裸体淫交视频免费看高清| 色视频www国产| 国产伦在线观看视频一区| 少妇被粗大猛烈的视频| 国产老妇伦熟女老妇高清| 人妻 亚洲 视频| 国产真实伦视频高清在线观看| 亚洲精品久久午夜乱码| 亚洲在线观看片| 午夜福利高清视频| 亚洲欧美日韩卡通动漫| 国产日韩欧美在线精品| 久久99蜜桃精品久久| 国产高清有码在线观看视频| 国产成人精品一,二区| 国内少妇人妻偷人精品xxx网站| 狂野欧美激情性xxxx在线观看| 国产视频内射| 亚洲成色77777| 国产成人免费观看mmmm| 日本黄色片子视频| 一级毛片黄色毛片免费观看视频| 嫩草影院精品99| 超碰97精品在线观看| 伦理电影大哥的女人| 国产精品久久久久久精品电影| 成人无遮挡网站| 男女那种视频在线观看| 2022亚洲国产成人精品| 成人无遮挡网站| 成人欧美大片| 99久久精品热视频| 午夜激情福利司机影院| 九九在线视频观看精品| 成人漫画全彩无遮挡| 亚洲久久久久久中文字幕| 寂寞人妻少妇视频99o| 男女国产视频网站| 亚洲国产精品国产精品| 自拍偷自拍亚洲精品老妇| 亚洲国产精品成人综合色| 国产成人精品一,二区| 汤姆久久久久久久影院中文字幕| 婷婷色综合www| 国产精品人妻久久久久久| 哪个播放器可以免费观看大片| 久久午夜福利片| 青春草视频在线免费观看| 亚洲真实伦在线观看| 欧美成人a在线观看| 免费人成在线观看视频色| 国产一区二区三区av在线| 色播亚洲综合网| 国产色爽女视频免费观看| 国产高潮美女av| 我要看日韩黄色一级片| 精品一区二区三区视频在线| 中文乱码字字幕精品一区二区三区| 日本欧美国产在线视频| 精品久久久噜噜| 亚洲真实伦在线观看| 久久综合国产亚洲精品| 午夜精品一区二区三区免费看| 日本-黄色视频高清免费观看| 国国产精品蜜臀av免费| 午夜精品国产一区二区电影 | 午夜福利在线在线| 一边亲一边摸免费视频| 香蕉精品网在线| 伊人久久精品亚洲午夜| 精品久久久久久电影网| 国产亚洲一区二区精品| 亚洲av二区三区四区| 午夜激情久久久久久久| 一级黄片播放器| 久久亚洲国产成人精品v| 国产精品福利在线免费观看| 亚洲最大成人手机在线| 亚洲,欧美,日韩| 色吧在线观看| 简卡轻食公司| 最近中文字幕2019免费版| 免费不卡的大黄色大毛片视频在线观看| 午夜激情福利司机影院| 国产毛片a区久久久久| 久久人人爽人人爽人人片va| 国产精品久久久久久精品古装| 国产精品久久久久久久电影| 亚洲欧洲国产日韩| 三级国产精品片| 人妻 亚洲 视频| 日本色播在线视频| 亚洲精品第二区| 80岁老熟妇乱子伦牲交| 男女边摸边吃奶| 午夜亚洲福利在线播放| 亚洲人成网站高清观看| 国产黄频视频在线观看| 亚洲国产色片| 国产成人一区二区在线| 精品久久国产蜜桃| 在线亚洲精品国产二区图片欧美 | 搡老乐熟女国产| 99热这里只有是精品在线观看| 人体艺术视频欧美日本| 久久久久久国产a免费观看| 欧美高清性xxxxhd video| 你懂的网址亚洲精品在线观看| 国产精品爽爽va在线观看网站| 久久精品久久久久久久性| 少妇的逼水好多| 成人国产av品久久久| 大香蕉久久网| 亚洲精品日韩av片在线观看| 乱系列少妇在线播放| 国产日韩欧美在线精品| 五月伊人婷婷丁香| 欧美3d第一页| 高清午夜精品一区二区三区| 亚洲国产精品成人综合色| 又黄又爽又刺激的免费视频.| 日韩欧美精品v在线| 国产成人a∨麻豆精品| 午夜免费男女啪啪视频观看| 男女下面进入的视频免费午夜| 亚洲国产精品成人久久小说| 久久精品国产亚洲av涩爱| 国产男人的电影天堂91| 欧美3d第一页| av在线观看视频网站免费| 久久精品久久精品一区二区三区| 日韩成人伦理影院| 亚洲av男天堂| 中文字幕制服av| 成人免费观看视频高清| 欧美性猛交╳xxx乱大交人| 在线亚洲精品国产二区图片欧美 | av女优亚洲男人天堂| 成人欧美大片| 国产欧美另类精品又又久久亚洲欧美| 成年免费大片在线观看| 久久久久久九九精品二区国产| 丰满人妻一区二区三区视频av| 国产美女午夜福利| 人人妻人人澡人人爽人人夜夜| 美女国产视频在线观看| 在线免费十八禁| 小蜜桃在线观看免费完整版高清| 亚洲一级一片aⅴ在线观看| 人妻系列 视频| 黄色日韩在线| 日韩成人av中文字幕在线观看| 久久精品国产鲁丝片午夜精品| 亚洲欧美精品专区久久| 一区二区三区乱码不卡18| 十八禁网站网址无遮挡 | 大话2 男鬼变身卡| 精品久久久精品久久久| 97超碰精品成人国产| 九色成人免费人妻av| 久久精品夜色国产| 亚洲精品视频女| 舔av片在线| 国产大屁股一区二区在线视频| 26uuu在线亚洲综合色| 亚洲av免费高清在线观看| 久久久久久久久久久免费av| 搡老乐熟女国产| 日韩一区二区视频免费看| 亚洲美女视频黄频| 青春草国产在线视频| 在线观看av片永久免费下载| 性插视频无遮挡在线免费观看| 91狼人影院| 国产精品一及| 最近最新中文字幕大全电影3| 自拍偷自拍亚洲精品老妇| 亚洲人与动物交配视频| av福利片在线观看| 日本爱情动作片www.在线观看| 久久久欧美国产精品| 99久久精品热视频| 人人妻人人看人人澡| 亚洲一区二区三区欧美精品 | 精品熟女少妇av免费看| 美女xxoo啪啪120秒动态图| 看黄色毛片网站| 亚洲伊人久久精品综合| 亚洲精品一区蜜桃| 久久久亚洲精品成人影院| 日本欧美国产在线视频| 午夜福利在线观看免费完整高清在| 国产高清有码在线观看视频| 久久精品国产亚洲av涩爱| 又粗又硬又长又爽又黄的视频| av在线app专区| 神马国产精品三级电影在线观看| 国产精品一区www在线观看| 在线精品无人区一区二区三 | 男女边吃奶边做爰视频| 成人综合一区亚洲| 欧美精品人与动牲交sv欧美| 人妻少妇偷人精品九色| 亚洲欧美精品自产自拍| 精品久久久久久久末码| 久久久久久久久久人人人人人人| 久久影院123| 日本一本二区三区精品| 国产精品av视频在线免费观看| 在线 av 中文字幕| 99精国产麻豆久久婷婷| 精品久久久精品久久久| 国产毛片a区久久久久| 久久97久久精品| 成年女人在线观看亚洲视频 | 亚洲最大成人手机在线| 国产成人a区在线观看| 99热网站在线观看| 欧美 日韩 精品 国产| 亚洲欧美精品专区久久| 97精品久久久久久久久久精品| 日韩人妻高清精品专区| av国产久精品久网站免费入址| 女人久久www免费人成看片| 各种免费的搞黄视频| av在线app专区| 国产精品福利在线免费观看| 欧美 日韩 精品 国产| 欧美日韩精品成人综合77777| 午夜爱爱视频在线播放| 亚洲精品自拍成人| 日本欧美国产在线视频| 国内少妇人妻偷人精品xxx网站| 97热精品久久久久久| 网址你懂的国产日韩在线| 在线观看三级黄色| 亚洲精品一区蜜桃| 久久国产乱子免费精品| 久久久久久久大尺度免费视频| 欧美日韩视频精品一区| 国产精品国产三级国产专区5o| 国产黄色免费在线视频| 在现免费观看毛片| 十八禁网站网址无遮挡 | 成年版毛片免费区| 久久久精品94久久精品| 91精品一卡2卡3卡4卡| 国内揄拍国产精品人妻在线| 乱码一卡2卡4卡精品| 美女主播在线视频| 国产91av在线免费观看| 国产精品成人在线| 少妇人妻精品综合一区二区| 日韩欧美 国产精品| 久久久久精品性色| 全区人妻精品视频| 久久国内精品自在自线图片| 麻豆久久精品国产亚洲av| 男男h啪啪无遮挡| 在线a可以看的网站| 色视频在线一区二区三区| 成人高潮视频无遮挡免费网站| 狂野欧美激情性xxxx在线观看| 尾随美女入室| 少妇的逼好多水| 80岁老熟妇乱子伦牲交| 欧美日韩视频高清一区二区三区二| 久久久a久久爽久久v久久| 亚洲一区二区三区欧美精品 | 少妇人妻一区二区三区视频| 免费观看的影片在线观看| 伊人久久精品亚洲午夜| 午夜亚洲福利在线播放| 久久久欧美国产精品| 激情 狠狠 欧美| av线在线观看网站| 亚洲一级一片aⅴ在线观看| 熟妇人妻不卡中文字幕| www.色视频.com| 嘟嘟电影网在线观看| 精品久久久久久久久亚洲| 麻豆成人午夜福利视频| 一区二区三区精品91| 又爽又黄无遮挡网站| 久久精品人妻少妇| 极品教师在线视频| 亚洲激情五月婷婷啪啪| 欧美日韩亚洲高清精品| 日日撸夜夜添| av线在线观看网站| 一级片'在线观看视频| 另类亚洲欧美激情| a级毛片免费高清观看在线播放| 人妻一区二区av| 欧美性感艳星| 日韩成人伦理影院| 91精品伊人久久大香线蕉| a级一级毛片免费在线观看| 亚洲国产欧美在线一区| 另类亚洲欧美激情| 精品一区在线观看国产| 不卡视频在线观看欧美| 一区二区三区免费毛片| 国产白丝娇喘喷水9色精品| 婷婷色综合www| av网站免费在线观看视频| 午夜福利视频1000在线观看| 久久久久久九九精品二区国产| 成人无遮挡网站| 日韩欧美 国产精品| 国产精品一二三区在线看| 天天躁夜夜躁狠狠久久av| 人妻系列 视频| 久久久久久久久久久免费av| 国产亚洲午夜精品一区二区久久 | 人妻制服诱惑在线中文字幕| 青青草视频在线视频观看| 免费黄频网站在线观看国产| 免费电影在线观看免费观看| 日韩一区二区视频免费看| 卡戴珊不雅视频在线播放| 亚洲天堂国产精品一区在线| 一级毛片 在线播放| 国产久久久一区二区三区| 国产精品久久久久久精品电影| 九九爱精品视频在线观看| 99热国产这里只有精品6| 特大巨黑吊av在线直播| 日韩中字成人| 国产精品伦人一区二区| 日产精品乱码卡一卡2卡三| 国产人妻一区二区三区在| 久久久久国产网址| 国产视频首页在线观看| 久久人人爽av亚洲精品天堂 | 国产精品一区www在线观看| 欧美日韩亚洲高清精品| 国产有黄有色有爽视频| 亚洲精品久久午夜乱码| 亚洲国产精品999| 最新中文字幕久久久久| 欧美xxxx性猛交bbbb| 国产精品国产av在线观看| 国内揄拍国产精品人妻在线| 国产视频内射| 三级国产精品欧美在线观看| 成人漫画全彩无遮挡| 国产精品麻豆人妻色哟哟久久| 久久99蜜桃精品久久| 热re99久久精品国产66热6| 欧美xxⅹ黑人| 国产精品99久久99久久久不卡 | 美女主播在线视频| 精品久久久久久电影网| 欧美3d第一页| 中文字幕av成人在线电影| 国产精品久久久久久av不卡| 亚洲性久久影院| 各种免费的搞黄视频| 国产亚洲av片在线观看秒播厂| 国产亚洲5aaaaa淫片| 最近中文字幕高清免费大全6| 精品久久久久久久人妻蜜臀av| 亚洲精品日韩在线中文字幕| 欧美激情在线99| 亚洲av不卡在线观看| 欧美精品一区二区大全| 亚洲国产精品成人综合色| 久热这里只有精品99| 日韩大片免费观看网站| 欧美xxxx性猛交bbbb| 国产午夜福利久久久久久| 亚洲伊人久久精品综合| 好男人在线观看高清免费视频| 免费看不卡的av| 九九久久精品国产亚洲av麻豆| 中文乱码字字幕精品一区二区三区| 成人综合一区亚洲| 草草在线视频免费看| 久久久色成人| 国产精品一区二区三区四区免费观看| 久久影院123| 亚洲三级黄色毛片| 一边亲一边摸免费视频| 国产欧美亚洲国产| 男男h啪啪无遮挡| 日韩人妻高清精品专区| 最近的中文字幕免费完整| 久久精品国产自在天天线| 在线亚洲精品国产二区图片欧美 | 内地一区二区视频在线| 成人亚洲精品一区在线观看 | 狠狠精品人妻久久久久久综合| 亚洲精品乱码久久久久久按摩| 一二三四中文在线观看免费高清| 黄色欧美视频在线观看| 黑人高潮一二区| 狠狠精品人妻久久久久久综合| 国产精品无大码| 久久精品国产亚洲av涩爱| 嫩草影院新地址| 女人久久www免费人成看片| 久久韩国三级中文字幕| 亚洲欧美成人精品一区二区| 综合色av麻豆| 国产 一区 欧美 日韩| 3wmmmm亚洲av在线观看| 成人黄色视频免费在线看| 一级毛片我不卡| 国产欧美亚洲国产| 六月丁香七月| 人人妻人人看人人澡| 又粗又硬又长又爽又黄的视频| 秋霞伦理黄片| 亚洲一级一片aⅴ在线观看| 亚洲四区av| 国产老妇女一区| 中国国产av一级| 欧美潮喷喷水| 国产欧美另类精品又又久久亚洲欧美| 蜜桃亚洲精品一区二区三区| 九九在线视频观看精品| 大香蕉97超碰在线| 国产精品一区二区在线观看99| 亚洲av电影在线观看一区二区三区 | 成人亚洲精品av一区二区| 2021天堂中文幕一二区在线观| 亚洲成人一二三区av| 午夜日本视频在线| 建设人人有责人人尽责人人享有的 | 嫩草影院精品99| 国产成人精品一,二区| 精品久久国产蜜桃| 中文字幕久久专区| 丰满乱子伦码专区| 青春草视频在线免费观看| 夫妻午夜视频| 成人综合一区亚洲| 久久99热这里只频精品6学生| 人人妻人人澡人人爽人人夜夜| 青春草国产在线视频| 26uuu在线亚洲综合色| 午夜免费男女啪啪视频观看| 欧美一区二区亚洲| 亚洲精品aⅴ在线观看| av线在线观看网站| 久久久国产一区二区| 女人久久www免费人成看片| 欧美成人一区二区免费高清观看| 欧美三级亚洲精品| 亚洲人成网站在线播| 国产成人精品久久久久久| 亚洲人成网站在线播| 综合色av麻豆| 三级经典国产精品| 久久精品久久久久久久性| 高清av免费在线| 亚洲丝袜综合中文字幕| 天天躁日日操中文字幕| 国产精品成人在线| 亚洲精品乱码久久久v下载方式| 亚洲精品乱久久久久久| 99久久精品一区二区三区| 久久久久久国产a免费观看| 三级男女做爰猛烈吃奶摸视频| 又爽又黄无遮挡网站| av免费观看日本| 秋霞伦理黄片| 春色校园在线视频观看| 欧美最新免费一区二区三区| av女优亚洲男人天堂| 国产成人免费观看mmmm| 国产在线男女| 亚洲综合精品二区| 三级男女做爰猛烈吃奶摸视频| 国产精品久久久久久久久免|