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

    風(fēng)電并網(wǎng)系統(tǒng)區(qū)域間概率可用輸電能力計(jì)算

    2014-08-16 12:51:48陳厚合李國(guó)慶張芳晶
    電力系統(tǒng)保護(hù)與控制 2014年21期
    關(guān)鍵詞:蒙特卡羅出力風(fēng)電場(chǎng)

    陳厚合,李國(guó)慶,張芳晶

    ?

    風(fēng)電并網(wǎng)系統(tǒng)區(qū)域間概率可用輸電能力計(jì)算

    陳厚合,李國(guó)慶,張芳晶

    (東北電力大學(xué),吉林 吉林 132012)

    為了準(zhǔn)確評(píng)估風(fēng)電場(chǎng)接入電網(wǎng)對(duì)系統(tǒng)可用輸電能力(Available Transfer Capability, ATC)的影響,針對(duì)風(fēng)電并網(wǎng)系統(tǒng)的概率可用輸電能力計(jì)算展開(kāi)研究,詳細(xì)分析不同風(fēng)電并網(wǎng)情況下ATC的變化規(guī)律。首先基于風(fēng)速Weibull分布,建立了大型風(fēng)電場(chǎng)輸出功率數(shù)學(xué)模型;進(jìn)而采用原-對(duì)偶內(nèi)點(diǎn)法完成風(fēng)電并網(wǎng)系統(tǒng)可用輸電能力單一樣板值的求解。在此基礎(chǔ)上,采用蒙特卡羅仿真法從廣義角度對(duì)風(fēng)電并網(wǎng)系統(tǒng)的ATC進(jìn)行概率評(píng)估。仿真結(jié)果表明,所提算法能夠有效評(píng)估風(fēng)電這種波動(dòng)性電源對(duì)ATC的影響。研究成果可為風(fēng)電并網(wǎng)系統(tǒng)安全經(jīng)濟(jì)性能評(píng)估提供有效參考信息,對(duì)未來(lái)電網(wǎng)規(guī)劃擴(kuò)建具有指導(dǎo)意義。

    電力系統(tǒng);風(fēng)電場(chǎng);可用輸電能力;蒙特卡羅仿真;原對(duì)偶內(nèi)點(diǎn)法

    0 引言

    風(fēng)能是當(dāng)前世界上最具大規(guī)模商業(yè)化開(kāi)發(fā)潛力的可再生能源,風(fēng)電并網(wǎng)運(yùn)行是大規(guī)模利用風(fēng)能的最重要途徑。然而,由于風(fēng)能所特有的隨機(jī)性、間歇性和波動(dòng)性特點(diǎn),受自然風(fēng)力驅(qū)動(dòng)的風(fēng)電機(jī)組輸出功率難以日前調(diào)度,這給互聯(lián)系統(tǒng)區(qū)域間輸電能力(Available Transfer Capability,ATC)計(jì)算和分析帶來(lái)了新的困難。因此,在大規(guī)模風(fēng)電集中并網(wǎng)條件下,如何準(zhǔn)確地把握風(fēng)電接入對(duì)系統(tǒng)區(qū)域間可用輸電能力的影響,在系統(tǒng)安全可靠運(yùn)行極限的范圍內(nèi),盡可能地滿(mǎn)足用電負(fù)荷的電力需求,對(duì)維持電力市場(chǎng)的正常運(yùn)營(yíng)及電力系統(tǒng)的安全穩(wěn)定運(yùn)行有著重要的意義。

    可用輸電能力是反映電網(wǎng)安全穩(wěn)定裕度并確保電力市場(chǎng)交易順利進(jìn)行的重要指標(biāo)。其現(xiàn)有的研究方法可分為兩類(lèi):基于概率的求解方法和確定性的求解方法?;诟怕市缘那蠼夥椒軌蛉娴孛枋霾淮_定性因素影響下的電網(wǎng)輸電能力,從而得到各種反映ATC充裕度統(tǒng)計(jì)指標(biāo)及其概率分布的曲線,是系統(tǒng)長(zhǎng)期規(guī)劃研究的有效工具;確定性的求解方法是概率方法中單一樣板的求解方法,即采用優(yōu)化技術(shù)或其他方法直接獲得所描述問(wèn)題的解。文獻(xiàn)[8]以準(zhǔn)確地計(jì)及系統(tǒng)中不確定性因素對(duì)ATC影響為目標(biāo),提出了蒙特卡羅仿真結(jié)合線性規(guī)劃的方法對(duì)ATC 進(jìn)行概率評(píng)估的思想,準(zhǔn)確地模擬了系統(tǒng)各種可能出現(xiàn)的運(yùn)行狀態(tài),詳細(xì)地分析了系統(tǒng)中不確定因素(系統(tǒng)元件故障、負(fù)荷波動(dòng)等)對(duì)電網(wǎng)可用輸電能力計(jì)算的影響。

    隨著我國(guó)風(fēng)電并網(wǎng)規(guī)模的不斷擴(kuò)增,電網(wǎng)風(fēng)電穿透功率的不斷增大,風(fēng)電出力的不確定性也成為影響ATC的一個(gè)重要因素。一些學(xué)者已經(jīng)展開(kāi)關(guān)于風(fēng)電并網(wǎng)條件下的電力系統(tǒng)區(qū)域間可用輸電能力的探索研究工作,文獻(xiàn)[9-10]建立了加入異步風(fēng)力發(fā)電機(jī)模型的含參潮流模型,推導(dǎo)出含風(fēng)電場(chǎng)注入功率項(xiàng)的全注入空間靜態(tài)電壓穩(wěn)定域邊界局部切平面解析式,著重分析了風(fēng)速概率分布參數(shù)對(duì)TTC的影響;文獻(xiàn)[11]采用 ARMA 時(shí)間序列模型模擬風(fēng)速,采用序貫蒙特卡羅仿真的方法對(duì)包含風(fēng)電場(chǎng)的 ATC 進(jìn)行概率評(píng)估,對(duì)抽樣狀態(tài)的 ATC 采用關(guān)鍵約束下的連續(xù)潮流來(lái)計(jì)算;文獻(xiàn)[12]建立了大型風(fēng)電場(chǎng)經(jīng)VSC-HVDC并網(wǎng)的ATC計(jì)算模型,分析了VSC控制策略,采用weibull分布模擬風(fēng)速,通過(guò)連續(xù)潮流法完成風(fēng)電并網(wǎng)系統(tǒng)ATC的確定性計(jì)算。

    本文首先對(duì)風(fēng)電場(chǎng)輸出功率波動(dòng)特性進(jìn)行分析,建立了大型風(fēng)電場(chǎng)輸出功率數(shù)學(xué)模型;進(jìn)而通過(guò)原-對(duì)偶內(nèi)點(diǎn)法完成風(fēng)電并網(wǎng)系統(tǒng)可用輸電能力單一樣板值的求解;在此基礎(chǔ)上,綜合考慮風(fēng)電出力隨機(jī)性、發(fā)輸電元件隨機(jī)故障、負(fù)荷波動(dòng)性和發(fā)電調(diào)度的不確定性等因素,采用非序貫蒙特卡羅仿真法對(duì)風(fēng)電并網(wǎng)系統(tǒng)的ATC進(jìn)行概率評(píng)估,建立了包括ATC期望、方差、最大值、最小值、風(fēng)險(xiǎn)概率和風(fēng)電影響因子等在內(nèi)的一系列評(píng)估指標(biāo),詳細(xì)分析了不同風(fēng)電并網(wǎng)情況下系統(tǒng)可用輸電能力的變化規(guī)律。

    1 風(fēng)電場(chǎng)出力模型

    1.1 風(fēng)速模型

    威布爾分布模型是國(guó)際上使用最多的風(fēng)速分布模型,大量實(shí)測(cè)數(shù)據(jù)表明,很多地區(qū)的風(fēng)速特征都可以用兩參數(shù)威布爾分布體現(xiàn),其概率密度函數(shù)和概率分布函數(shù)可以描述為

    (2)

    式中:為風(fēng)速;為形狀系數(shù),取值范圍為 1.5 ~ 2.5;為尺度系數(shù),表示該地區(qū)的平均風(fēng)速,取值一般在5~10 m/s。

    設(shè)是[0,1]上服從均勻分布的隨機(jī)變量,令=(),對(duì)風(fēng)速概率分布函數(shù)()進(jìn)行反變換,推得公式如式(3)。

    由于和1-都是在 [0,1]上均勻分布的隨機(jī)變量,用代替1-,即可得到風(fēng)速發(fā)生器模型,進(jìn)而產(chǎn)生風(fēng)速樣本。

    (4)

    1.2 風(fēng)電機(jī)組出力模型

    風(fēng)電機(jī)組的輸出功率主要取決于風(fēng)速的大小,二者之間的關(guān)系曲線即為風(fēng)機(jī)的功率-風(fēng)速特性曲線,該曲線可以由式(5)所示的分段函數(shù)表示。

    式中:為風(fēng)機(jī)輪轂高度處的風(fēng)速;為切入風(fēng)速;為切出風(fēng)速;為額定風(fēng)速;為額定輸出功率。當(dāng)風(fēng)速高于時(shí),風(fēng)機(jī)啟動(dòng)并網(wǎng)運(yùn)行;當(dāng)風(fēng)速大于時(shí),風(fēng)機(jī)保持額定功率輸出;當(dāng)風(fēng)速低于或高于時(shí),風(fēng)電機(jī)組停機(jī)并從電網(wǎng)解列出來(lái)。式(5)中、、為模型參數(shù),表達(dá)式如式(6)。

    (6)

    當(dāng)風(fēng)電機(jī)組額定功率為1.5 MW,切入風(fēng)速取4 m/s,切出風(fēng)速取25 m/s,額定風(fēng)速取13 m/s時(shí),風(fēng)力發(fā)電機(jī)組的功率-風(fēng)速關(guān)系曲線如圖1所示。

    1.3 風(fēng)電場(chǎng)出力模型

    本文以單臺(tái)機(jī)組輸出功率特性的線性放大來(lái)等效整個(gè)風(fēng)電場(chǎng)出力特性的基礎(chǔ)上,綜合考慮尾流效應(yīng)與風(fēng)電機(jī)組故障檢修對(duì)風(fēng)電場(chǎng)輸出功率的影響,得到風(fēng)電場(chǎng)的整體輸出功率如式(7)所示。

    式中:為尾流效應(yīng)系數(shù),此處取典型值10%;x為風(fēng)機(jī)狀態(tài),正常工作時(shí)為1,故障或檢修時(shí)為0;p為單臺(tái)風(fēng)電機(jī)組的輸出功率;為風(fēng)電場(chǎng)內(nèi)風(fēng)電機(jī)組的個(gè)數(shù)。

    目前我國(guó)風(fēng)電場(chǎng)應(yīng)用比較廣泛的風(fēng)力發(fā)電機(jī)有普通異步電機(jī)、雙饋異步電機(jī)和直驅(qū)永磁同步電機(jī)。普通異步電機(jī)是我國(guó)風(fēng)電場(chǎng)最早采用的機(jī)型,在現(xiàn)有的風(fēng)電場(chǎng)中也占有一定的比例,其結(jié)構(gòu)簡(jiǎn)單,不具備無(wú)功調(diào)節(jié)能力,但可以通過(guò)風(fēng)電場(chǎng)配置的無(wú)功補(bǔ)償設(shè)備來(lái)保證無(wú)功功率的平衡。雙饋感應(yīng)風(fēng)力發(fā)電機(jī)組和直驅(qū)永磁同步電機(jī)占據(jù)新增機(jī)組的很大比例,具有靈活的無(wú)功調(diào)節(jié)能力,可保持并網(wǎng)點(diǎn)的功率因數(shù)恒定。因此,在進(jìn)行穩(wěn)態(tài)分析時(shí)風(fēng)電場(chǎng)能夠按著系統(tǒng)操作人員給定的功率因數(shù)運(yùn)行,因此本文將風(fēng)電場(chǎng)視為恒功率因數(shù)方式下的PQ節(jié)點(diǎn)。

    2 基于蒙特卡羅法的ATC概率計(jì)算

    2.1 電力系統(tǒng)狀態(tài)確定

    電力系統(tǒng)的運(yùn)行狀態(tài)是系統(tǒng)內(nèi)所有設(shè)備狀態(tài)的組合,且每一個(gè)設(shè)備的狀態(tài)可由設(shè)備出現(xiàn)在該狀態(tài)的概率進(jìn)行抽樣來(lái)確定。因此,本文主要分析以下幾種不確定因素的概率分布規(guī)律,進(jìn)而確定系統(tǒng)的運(yùn)行狀態(tài)。

    (1) 系統(tǒng)中的常規(guī)發(fā)電機(jī)組、風(fēng)電機(jī)組和輸電線路有故障和運(yùn)行兩種運(yùn)行狀態(tài),其概率分布函數(shù)服從兩點(diǎn)分布,如式(8)所示。

    其中:x為發(fā)電機(jī)或線路的運(yùn)行狀態(tài),1代表該元件正常運(yùn)行,0代表退出運(yùn)行;A是在[0, 1]間抽取的隨機(jī)變量,若A小于指定的故障率,則此時(shí)該元件故障,若A大于故障率,則元件正常運(yùn)行。

    (2) 對(duì)于負(fù)荷及常規(guī)發(fā)電機(jī)出力的波動(dòng),認(rèn)為其服從正態(tài)分布N(,) 的規(guī)律。其中,為節(jié)點(diǎn)負(fù)荷或發(fā)電機(jī)出力的預(yù)測(cè)值;參數(shù)一般根據(jù)經(jīng)驗(yàn)值給定,它代表了系統(tǒng)實(shí)際負(fù)荷值或發(fā)電機(jī)出力值偏離給定值的大小。

    (3) 對(duì)于風(fēng)電場(chǎng)出力的不確定性,假定風(fēng)速服從weibull分布,根據(jù)式(5)、式(7)確定風(fēng)電場(chǎng)隨機(jī)出力。

    2.2 基于原-對(duì)偶內(nèi)點(diǎn)法的ATC計(jì)算模型

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

    本文ATC采用最優(yōu)潮流進(jìn)行求解,目標(biāo)函數(shù)為區(qū)域間A、B間聯(lián)絡(luò)線在基態(tài)潮流的基礎(chǔ)上,可進(jìn)一步傳輸?shù)淖畲箅娔?,如?9)所示。

    式中:為送電區(qū)域;為受電區(qū)域;P為從節(jié)點(diǎn)到節(jié)點(diǎn)的線路傳輸?shù)挠泄β?;上?biāo)“0”代表基態(tài)值;表示容量效益裕度CBM大小,取值5%。

    (2) 等式約束條件

    等式約束為系統(tǒng)常規(guī)潮流方程約束。

    式中:L、L為節(jié)點(diǎn)的負(fù)荷;g、g為風(fēng)電機(jī)組有功/無(wú)功發(fā)電容量;G、G為常規(guī)發(fā)電機(jī)組有功/無(wú)功發(fā)電容量;UU分別為節(jié)點(diǎn)和的電壓幅值;為節(jié)點(diǎn)和電壓的相位差;GB為節(jié)點(diǎn)導(dǎo)納陣中的元素;α為風(fēng)電場(chǎng)輸出的功率因數(shù)。

    (3) 不等式約束條件

    不等式約束條件主要包含發(fā)電機(jī)出力約束、熱穩(wěn)定約束、節(jié)點(diǎn)電壓約束和負(fù)荷出力約束等,相應(yīng)的數(shù)學(xué)表達(dá)形式如式(11)。

    式中:是受電區(qū)的所有負(fù)荷節(jié)點(diǎn)集合;和分別是送電區(qū)的所有常規(guī)機(jī)組節(jié)點(diǎn)集合和風(fēng)電場(chǎng)節(jié)點(diǎn)集合;為線路集合;為總節(jié)點(diǎn)集合;P對(duì)應(yīng)線路有功潮流;變量上角標(biāo)“0”、“min”、“max”分別表示基態(tài)潮流中的值、變量的下限和上限值。

    2.3 ATC概率評(píng)估指標(biāo)

    為研究風(fēng)電并網(wǎng)對(duì)系統(tǒng)ATC的影響,利用統(tǒng)計(jì)分析的方法,求取下列ATC概率統(tǒng)計(jì)指標(biāo)和風(fēng)險(xiǎn)指標(biāo),以便定量分析電網(wǎng)的可用輸電能力。

    1) ATC最大值

    2) ATC最小值

    (13)

    3) ATC期望值

    4) ATC方差

    (15)

    5) ATC變異系數(shù)

    6) ATC 等于零的概率

    (17)

    7) ATC概率密度

    輸電能力概率密度是指,在蒙特卡羅模擬中,某一ATC值出現(xiàn)的次數(shù)與抽樣總次數(shù)的比值。

    8) ATC不足概率

    輸電能力不足概率是對(duì)概率密度函數(shù)的積分,表示ATC值小于等于某一設(shè)定ATC值的累積概率,也稱(chēng)為風(fēng)險(xiǎn)概率。

    9)風(fēng)電影響因子

    風(fēng)電場(chǎng)影響因子表示風(fēng)電接入對(duì)ATC風(fēng)險(xiǎn)概率的影響,[C<] 表示無(wú)風(fēng)場(chǎng)時(shí)的ATC的值小于特定值的仿真次數(shù),[C<] 表示加入風(fēng)場(chǎng)后ATC的值小于特定值的抽樣次數(shù)。

    2.4 基于蒙特卡羅法的ATC計(jì)算流程

    基于非序貫蒙特卡羅法的ATC評(píng)估流程如下所示。

    1) 輸入系統(tǒng)中發(fā)電機(jī)、線路、負(fù)荷以及風(fēng)電場(chǎng)的相關(guān)參數(shù),形成系統(tǒng)的基本信息。

    2) 確定最大抽樣次數(shù)。

    3) 采用蒙特卡羅法對(duì)系統(tǒng)各種不確定因素進(jìn)行隨機(jī)抽樣,包括風(fēng)電場(chǎng)出力、常規(guī)發(fā)電機(jī)出力波動(dòng)、負(fù)荷出力波動(dòng)及發(fā)輸電元件隨機(jī)故障等。

    4) 針對(duì)每次抽樣得到的系統(tǒng)狀態(tài),若系統(tǒng)出現(xiàn)了系統(tǒng)解列(孤島運(yùn)行)、或發(fā)電量與負(fù)荷量不匹配的情況,則認(rèn)為此種情況下的ATC值為零,返回步驟3),重新抽樣;如果系統(tǒng)能夠正常運(yùn)行,則繼續(xù)下一步。

    5) 對(duì)于可行的結(jié)構(gòu),采用原-對(duì)偶內(nèi)點(diǎn)法對(duì)單一樣板情況下的可用輸電能力進(jìn)行求解;若抽樣次數(shù)小于指定抽樣次數(shù),則返回到步驟3),為下一次抽樣做準(zhǔn)備。

    6) 若達(dá)到了指定抽樣次數(shù),則記錄ATC相關(guān)數(shù)據(jù),并計(jì)算其概率評(píng)估指標(biāo)。

    3 算例分析

    為驗(yàn)證所建立ATC計(jì)算模型及采用優(yōu)化算法的正確性和有效性,本文采用Matlab語(yǔ)言實(shí)現(xiàn)本課題算法,對(duì)IEEE-30節(jié)點(diǎn)系統(tǒng)進(jìn)行仿真分析,網(wǎng)絡(luò)相關(guān)參數(shù)詳見(jiàn)文獻(xiàn)[17]。IEEE-30節(jié)點(diǎn)系統(tǒng)是電力系統(tǒng)研究領(lǐng)域中廣泛應(yīng)用的標(biāo)準(zhǔn)測(cè)試系統(tǒng)之一,共有6臺(tái)發(fā)電機(jī),41條線路,劃分為3個(gè)區(qū),系統(tǒng)網(wǎng)絡(luò)接線如圖2所示。其中,風(fēng)電場(chǎng)的相關(guān)參數(shù)見(jiàn)表1,此時(shí)假定風(fēng)速服從尺度參數(shù)為10、形狀參數(shù)為2的雙參數(shù)Weibull分布;蒙特卡羅抽樣過(guò)程中的不確定因素參數(shù)設(shè)置如表2所示。設(shè)定抽樣次數(shù)為5 000次,采用標(biāo)幺值進(jìn)行計(jì)算,基準(zhǔn)容量取100 MW。

    圖2 IEEE30節(jié)點(diǎn)系統(tǒng)結(jié)構(gòu)圖

    表1 風(fēng)場(chǎng)相關(guān)參數(shù)

    表2 蒙特卡羅抽樣過(guò)程中系統(tǒng)參數(shù)

    根據(jù)風(fēng)電場(chǎng)并網(wǎng)條件的不同,大型風(fēng)電場(chǎng)對(duì)ATC 的影響呈現(xiàn)出不同的特點(diǎn)。本文著重評(píng)估風(fēng)電場(chǎng)的出力的隨機(jī)性和波動(dòng)性對(duì)ATC的影響,分以下五種情況對(duì)比研究。

    1)案例 1是原系統(tǒng)總裝機(jī)容量不變,無(wú)風(fēng)電場(chǎng)的情況;

    2)案例 2 是在送電區(qū)域節(jié)點(diǎn)5并入一個(gè)容量為20×1.5 MW的風(fēng)場(chǎng);

    3)案例 3是在送電區(qū)域節(jié)點(diǎn)5并入兩個(gè)容量為20×1.5 MW的風(fēng)場(chǎng);

    4)案例 4是在受電區(qū)域節(jié)點(diǎn)8并入一個(gè)容量為20×1.5 MW的風(fēng)場(chǎng);

    5)案例 5 是在送電區(qū)域節(jié)點(diǎn)5直接并入一個(gè)容量為30 MW火電場(chǎng)。

    分別對(duì)以上五種情況下的ATC進(jìn)行統(tǒng)計(jì)分析,得到ATC的各項(xiàng)概率指標(biāo)如表3所示。此時(shí)在計(jì)算ATC不足概率時(shí),設(shè)定ATC指定值的大小為基態(tài)潮流下的可用輸電能力值79.73 MW。

    表3 五種風(fēng)電并網(wǎng)條件下ATC統(tǒng)計(jì)指標(biāo)對(duì)比

    (1) 方案1、方案2和方案3的對(duì)比分析體現(xiàn)了風(fēng)電場(chǎng)裝機(jī)容量的變化對(duì)ATC的影響,其概率密度分布曲線和累積分布曲線,如圖3、圖4所示。

    圖3 不同并網(wǎng)風(fēng)電場(chǎng)容量下的ATC概率密度分布曲線

    圖4 不同并網(wǎng)風(fēng)電場(chǎng)容量下的ATC累積分布曲線

    通過(guò)對(duì)以上ATC各項(xiàng)指標(biāo)的對(duì)比分析,可以得出如下結(jié)論:

    a) 風(fēng)電場(chǎng)接入系統(tǒng)后,ATC的期望值增大了,而且隨著風(fēng)電場(chǎng)容量的增加,ATC期望值有增大的趨勢(shì),這說(shuō)明風(fēng)電場(chǎng)并網(wǎng)能夠提高系統(tǒng)的可用輸電能力值。但由于受到輸電線路熱容量的限制,系統(tǒng)ATC增大數(shù)值并不與風(fēng)電場(chǎng)容量成正比關(guān)系,且當(dāng)風(fēng)場(chǎng)裝機(jī)容量超過(guò)一定值后,ATC值將不會(huì)再隨之增大。

    b) 風(fēng)電場(chǎng)并網(wǎng)將導(dǎo)致ATC的方差增大,系統(tǒng)的波動(dòng)性增強(qiáng),這是由于風(fēng)電功率的隨機(jī)性和間歇性使系統(tǒng)的不確定因素增多的原因,風(fēng)電并網(wǎng)容量越大,ATC值的波動(dòng)幅度就越強(qiáng)。此時(shí)變異系數(shù)的變化則說(shuō)明了風(fēng)電并網(wǎng)影響了ATC計(jì)算精度。

    c) 由表3或圖4都可以看出,隨著風(fēng)電并網(wǎng)容量的增大,系統(tǒng)ATC不足概率逐漸減少,說(shuō)明風(fēng)電并網(wǎng)后雖然引起ATC的波動(dòng)程度,但也一定量地提高了系統(tǒng)ATC值。同時(shí),隨著風(fēng)電并網(wǎng)容量的提升,風(fēng)電影響因子對(duì)電網(wǎng)的影響也不斷增大,這也說(shuō)明了風(fēng)電并網(wǎng)在一定程度上減小了ATC風(fēng)險(xiǎn)概率,提高了系統(tǒng)的可靠性。

    (2) 方案2和方案4的對(duì)比分析體現(xiàn)了風(fēng)電場(chǎng)并網(wǎng)位置的變化對(duì)ATC的影響。其概率密度分布曲線和累積分布曲線,如圖5、圖6所示。

    圖5 不同風(fēng)電場(chǎng)并網(wǎng)位置時(shí)的ATC概率密度分布曲線

    圖6 不同風(fēng)電場(chǎng)并網(wǎng)位置時(shí)的ATC累積分布曲線

    結(jié)果分析:風(fēng)電場(chǎng)并入系統(tǒng)位置的不同對(duì)ATC各項(xiàng)概率指標(biāo)有著不同的影響。與送電側(cè)(節(jié)點(diǎn)5)相比,在受電側(cè)(節(jié)點(diǎn)8)并入風(fēng)電場(chǎng),會(huì)令A(yù)TC的期望值增加效果明顯,但是ATC的方差和風(fēng)電影響因子的變化量也較大,這說(shuō)明風(fēng)電在受電側(cè)并網(wǎng)時(shí)對(duì)系統(tǒng)可用輸電能力的影響較大。因此,在系統(tǒng)規(guī)劃時(shí)應(yīng)該考慮在系統(tǒng)的送電側(cè)進(jìn)行風(fēng)電并網(wǎng),減少風(fēng)電并網(wǎng)對(duì)可用輸電能力的影響。

    (3) 方案2和方案5表示相同容量的風(fēng)電機(jī)組和火電機(jī)組并入系統(tǒng)的情況。其概率密度分布曲線和累積分布曲線,如圖7、圖8所示。

    圖7 相同容量風(fēng)電場(chǎng)和火電機(jī)組并網(wǎng)時(shí)ATC概率密度分布曲線

    圖8 相同容量風(fēng)電場(chǎng)和火電機(jī)組并網(wǎng)時(shí)ATC累積分布曲線

    結(jié)果分析:風(fēng)電場(chǎng)或與其額定容量相等的火電機(jī)組接入節(jié)點(diǎn)5以后,系統(tǒng)的ATC期望值均增大,而且火電機(jī)組并網(wǎng)后ATC的增量要大于風(fēng)電機(jī)組并網(wǎng)情況;同時(shí),火電機(jī)組并網(wǎng)后系統(tǒng)ATC的方差、變異系數(shù)等均小于風(fēng)電接入的情況,這體現(xiàn)了風(fēng)電這種波動(dòng)性的電源特性,其可靠性低于常規(guī)發(fā)電機(jī)組。

    4 結(jié)論

    本文采用蒙特卡羅仿真法從廣義的角度對(duì)風(fēng)電并網(wǎng)后系統(tǒng)區(qū)域間的輸電能力進(jìn)行求解分析。研究了風(fēng)速模型、尾流效應(yīng)和風(fēng)電機(jī)組有功出力三者的對(duì)應(yīng)關(guān)系,建立了大型風(fēng)電場(chǎng)輸出功率數(shù)學(xué)模型;并通過(guò)原-對(duì)偶內(nèi)點(diǎn)法完成風(fēng)電并網(wǎng)系統(tǒng)可用輸電能力單一樣板值的求解;最后通過(guò)所建立的一系列可用輸電能力評(píng)估指標(biāo)詳細(xì)分析了不同風(fēng)電并網(wǎng)情況下系統(tǒng)可用輸電能力的變化規(guī)律,對(duì)風(fēng)電并網(wǎng)系統(tǒng)的運(yùn)行規(guī)劃具有一定的指導(dǎo)意義。

    [1] 國(guó)家中長(zhǎng)期科學(xué)和技術(shù)發(fā)展規(guī)劃綱要(2006~2020年). [EB/OL]. http://ghs.ndrc.gov. cn/ghwb/default.htm.

    National medium and long-term plan for science and technology development (2006~2020)[EB/OL]. http://ghs. ndrc. gov. cn/ghwb/default.htm.

    [2] 李俊峰,蔡豐波,喬黎明,等. 中國(guó)風(fēng)電發(fā)展報(bào)告2013 [M]. 北京: 中國(guó)環(huán)境科學(xué)出版社, 2013.

    LI Jun-feng, CAI Feng-bo, QIAO Li-ming, et al. 2013 annual review and outlook on China wind power[M]. Beijing: China Environmental Science Press, 2013.

    [3] North American Electric Reliability Council (NERC). Available transfer capability definitions and determination[R]. 1996.

    [4] 陳厚合,李國(guó)慶,姜濤. 計(jì)及靜態(tài)電壓穩(wěn)定約束的交直流系統(tǒng)可用輸電能力[J]. 電網(wǎng)技術(shù), 2012, 36(2): 75-81.

    CHEN Hou-he, LI Guo-qing, JIANG Tao. Available transfer capability of hybrid AC/DC power system considering constraint of static voltage stability[J]. Power System Technology, 2012, 36(2): 75-81.

    [5] 梁海峰,郭然,高亞靜. 計(jì)及天氣因素的風(fēng)電并網(wǎng)系統(tǒng)可用輸電能力評(píng)估[J]. 電力系統(tǒng)保護(hù)與控制, 2013, 41(1): 164-168.

    LIANG Hai-feng, GUO Ran, GAO Ya-jing. Assessment on ATC of wind farm incorporated system considering weather[J]. Power System Protection and Control, 2013, 41(1): 164-168.

    [6] 崔雅莉,別朝紅,王錫凡. 輸電系統(tǒng)可用輸電能力的概率模型及計(jì)算[J]. 電力系統(tǒng)自動(dòng)化, 2003, 27(14): 36-40.

    CUI Ya-li, BIE Zhao-hong, WANG Xi-fan. Study on probabilistic models and algorithms of available transfer capability[J]. Automation of Electric Power Systems, 2003, 27(14): 36-40.

    [7] 王艷玲,韓學(xué)山,周曉峰. 基于擴(kuò)展潮流的輸電斷面最大傳輸能力[J]. 電力系統(tǒng)保護(hù)與控制, 2011, 39(13): 20-24.

    WANG Yan-ling, HAN Xue-shan, ZHOU Xiao-feng. The total transfer capability of cross-section of transmission line based on expanded power flow[J]. Power System Protection and Control, 2011, 39(13): 20-24.

    [8] 高雅靜. 考慮不確定因素的電網(wǎng)可用輸電能力的研究[D]. 北京: 華北電力大學(xué), 2007.

    GAO Ya-jing. Research on available transfer capability considering uncertainties[D]. Beijing: North China Electric Power University, 2007.

    [9] 王成山,孫瑋,王興剛. 含大型風(fēng)電場(chǎng)的電力系統(tǒng)最大輸電能力計(jì)算[J]. 電力系統(tǒng)自動(dòng)化, 2007, 31(2): 17-31.

    WANG Cheng-shan, SUN Wei, WANG Xing-gang. Total transfer capability calculation of power system including large scale wind farm[J]. Automation of Electric Power Systems, 2007, 31(2): 17-31.

    [10] 王成山,王興剛,孫瑋. 含大型風(fēng)電場(chǎng)的電力系統(tǒng)概率最大輸電能力快速計(jì)算[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2008, 28(10): 56-62.

    WANG Cheng-shan, WANG Xing-gang, SUN Wei. Fast calculation and analysis of probabilistic total transfer capability in power system including large-scale wind farms[J]. Proceedings of the CSEE, 2008, 28(10): 56-62.

    [11] 周明,冉瑞江,李庚銀. 風(fēng)電并網(wǎng)系統(tǒng)可用輸電能力的評(píng)估[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2010, 30(22): 14-21.

    ZHOU Ming, RAN Rui-jiang, LI Geng-yin. Assessment on available transfer capability of wind farm incorporated system[J]. Proceedings of the CSEE, 2010, 30(22): 14-21.

    [12] GAO Ya-jing, WANG Zhi, LIANG Hai-feng. Available transfer capability calculation with large offshore wind farms connected by VSC-HVDC[C] // Innovative Smart Grid Technologies-Asia (ISGT Asia) Conference, Tianjin, China, 1-6 May, 2012.

    [13] 鐘浩,唐民富. 風(fēng)電場(chǎng)發(fā)電可靠性及容量可信度評(píng)估[J]. 電力系統(tǒng)保護(hù)與控制, 2012, 40(18): 75-80.

    ZHONG Hao, TANG Min-fu. Assessment on the reliability and capacity credibility of wind farms[J]. Power System Protection and Control, 2012, 40(18): 75-80.

    [14] 彭鵠,謝開(kāi)貴,王進(jìn),等. 含風(fēng)力發(fā)電的發(fā)輸配電系統(tǒng)可靠性綜合評(píng)估[J]. 電氣應(yīng)用, 2009, 28(17): 44-49.

    PENG Hu, XIE Kai-gui, WANG Jin, et al. Reliability assessment of generation and transmission systems considering wind farm effects[J]. Electrotechnical Application, 2009, 28(17): 44-49.

    [15] 王俊,蔡興國(guó). 考慮分布式發(fā)電不確定性的可用輸電能力評(píng)估[J]. 電力系統(tǒng)自動(dòng)化, 2011, 35(21): 21-25.

    WANG Jun, CAI Xing-guo. Assessment of available transfer capability considering uncertainties distributed generator[J]. Automation of Electric Power Systems, 2011, 35(21): 21-25.

    [16] CHUNG C Y, YAN Wei, LIU Fang. Decomposed predictor-corrector interior point method for dynamic optimal power flow[J]. IEEE Trans on Power Systems, 2011, 26(3): 1030-1039.

    [17] 吳際舜, 侯志儉. 電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)算法[M]. 上海: 上海交通大學(xué)出版社, 2000.

    WU Ji-shun, HOU Zhi-jian. Calculation algorithms of power system flow[M]. Shanghai: Shanghai Jiao Tong University Press, 2000.

    Calculation of probabilistic available transfer capability of wind farm incorporated system

    CHEN Hou-he, LI Guo-qing, ZHANG Fang-jing

    (Northeast Dianli University, Jinlin 132012, China)

    In order to accurately assess the effect of wind farms on the available transfer capability (ATC), probabilistic available transfer capability of wind farm incorporated system is researched, and the change law of ATC with different wind farm incorporated system is analyzed in detail. First, the mathematical model of large scale wind farm is established based on the Weibull model of wind speed. Then, primal-dual interior point method is used to calculate the value of ATC in certain system operation mode. Finally, ATC probability assessment of wind power integrated system is achieved by Monte-carlo simulation method in a broad perspective. The simulation results show that the proposed algorithm can effectively assess the influence of wind farms to ATC, which has guiding functions for future expansion of the power grid and can provide effective reference information for the safety and economics evaluation of the system with wind farms.

    This work is supported by National Natural Science Foundation of China (No. 51377016).

    power system; wind farm; available transfer capability; Monte-carlo simulation; interior-point algorithm

    TM76

    A

    1674-3415(2014)21-0059-07

    2014-01-27

    陳厚合(1978-),男,博士研究生,講師,主要從事電力系統(tǒng)安全性和穩(wěn)定性分析的研究和教學(xué)工作;E-mail: chenhouhe_tju@163.com

    李國(guó)慶(1963-),男,教授,博士生導(dǎo)師,主要從事電力系統(tǒng)安全性和穩(wěn)定性分析、控制與決策、配電系統(tǒng)自動(dòng)化的研究工作;

    張芳晶(1988-),女,通訊作者,碩士研究生,主要研究方向?yàn)殡娏ο到y(tǒng)安全性和穩(wěn)定性分析。E-mail: zhangfangjing074@126.com

    國(guó)家自然科學(xué)基金項(xiàng)目(51377016)

    猜你喜歡
    蒙特卡羅出力風(fēng)電場(chǎng)
    利用蒙特卡羅方法求解二重積分
    基于PSS/E的風(fēng)電場(chǎng)建模與動(dòng)態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    風(fēng)電場(chǎng)有功出力的EEMD特性分析
    含風(fēng)電場(chǎng)電力系統(tǒng)的潮流計(jì)算
    要爭(zhēng)做出力出彩的黨員干部
    河南電力(2016年5期)2016-02-06 02:11:35
    風(fēng)電場(chǎng)群出力的匯聚效應(yīng)分析
    探討蒙特卡羅方法在解微分方程邊值問(wèn)題中的應(yīng)用
    探求風(fēng)電場(chǎng)的遠(yuǎn)景
    風(fēng)能(2015年9期)2015-02-27 10:15:25
    代力吉風(fēng)電場(chǎng)的我們
    風(fēng)能(2015年7期)2015-02-27 10:15:02
    復(fù)合型種子源125I-103Pd劑量場(chǎng)分布的蒙特卡羅模擬與實(shí)驗(yàn)測(cè)定
    同位素(2014年2期)2014-04-16 04:57:20
    国产真实伦视频高清在线观看| 国模一区二区三区四区视频| 国产极品精品免费视频能看的| 久久午夜亚洲精品久久| 神马国产精品三级电影在线观看| 中文字幕免费在线视频6| 深爱激情五月婷婷| 22中文网久久字幕| 真实男女啪啪啪动态图| 久久久久性生活片| 久久久久久久久久黄片| 赤兔流量卡办理| 成人av在线播放网站| 黑人高潮一二区| 日本-黄色视频高清免费观看| 最近最新中文字幕大全电影3| 91久久精品国产一区二区成人| 精品久久久久久久久亚洲| 国产精品不卡视频一区二区| 久久久久久久久久久免费av| 日本撒尿小便嘘嘘汇集6| 亚州av有码| 欧美日韩一区二区视频在线观看视频在线 | 99国产极品粉嫩在线观看| 精品人妻一区二区三区麻豆| 偷拍熟女少妇极品色| 国产精品无大码| av在线亚洲专区| 免费看av在线观看网站| 最近中文字幕高清免费大全6| 深夜a级毛片| 亚洲av二区三区四区| 91精品国产九色| 国产精品美女特级片免费视频播放器| 久久久久久久久久黄片| 你懂的网址亚洲精品在线观看 | 深爱激情五月婷婷| 人妻制服诱惑在线中文字幕| 亚洲国产色片| 日韩制服骚丝袜av| 国产av一区在线观看免费| 日日撸夜夜添| 看黄色毛片网站| 一个人免费在线观看电影| 国产亚洲5aaaaa淫片| 99精品在免费线老司机午夜| 好男人视频免费观看在线| 亚洲乱码一区二区免费版| 久久精品国产鲁丝片午夜精品| 成人午夜高清在线视频| 亚洲真实伦在线观看| 在线国产一区二区在线| 国产精品美女特级片免费视频播放器| а√天堂www在线а√下载| 亚洲欧洲国产日韩| 久久精品91蜜桃| 99久久精品一区二区三区| 亚洲精品久久国产高清桃花| 天堂√8在线中文| 久久99蜜桃精品久久| kizo精华| 啦啦啦啦在线视频资源| 亚洲av一区综合| 国产女主播在线喷水免费视频网站 | 国产麻豆成人av免费视频| 99在线人妻在线中文字幕| 神马国产精品三级电影在线观看| 久久久久久久久久久丰满| 六月丁香七月| 又粗又硬又长又爽又黄的视频 | 草草在线视频免费看| 女同久久另类99精品国产91| 久久久久久久久中文| 大又大粗又爽又黄少妇毛片口| 久久精品夜色国产| 岛国在线免费视频观看| 99久久中文字幕三级久久日本| 欧美精品一区二区大全| 尤物成人国产欧美一区二区三区| 久久九九热精品免费| 我要看日韩黄色一级片| 欧美日本亚洲视频在线播放| 婷婷色综合大香蕉| 2022亚洲国产成人精品| 午夜视频国产福利| 麻豆精品久久久久久蜜桃| 超碰av人人做人人爽久久| 26uuu在线亚洲综合色| 人体艺术视频欧美日本| 赤兔流量卡办理| 国内久久婷婷六月综合欲色啪| 欧美色视频一区免费| 久久久精品欧美日韩精品| 久久人妻av系列| 成年免费大片在线观看| 免费人成视频x8x8入口观看| 久久精品影院6| 亚洲成人精品中文字幕电影| 在线观看66精品国产| 国产在线精品亚洲第一网站| 可以在线观看的亚洲视频| 免费看美女性在线毛片视频| 免费av毛片视频| 国产一区二区三区av在线 | 麻豆国产97在线/欧美| 国产精品av视频在线免费观看| 国产精品一区二区性色av| 中国美白少妇内射xxxbb| 久久久a久久爽久久v久久| 成年av动漫网址| 麻豆乱淫一区二区| 久久久国产成人免费| 亚洲欧美日韩卡通动漫| 国产成人精品婷婷| 国产精品无大码| 免费观看在线日韩| 成人高潮视频无遮挡免费网站| 可以在线观看毛片的网站| 免费黄网站久久成人精品| 人妻久久中文字幕网| 国产亚洲欧美98| 久久精品夜色国产| 99久国产av精品国产电影| 卡戴珊不雅视频在线播放| 黄色视频,在线免费观看| 国产成人福利小说| 中文精品一卡2卡3卡4更新| 亚洲精品成人久久久久久| 黄色配什么色好看| 日日摸夜夜添夜夜爱| 老司机影院成人| 九九久久精品国产亚洲av麻豆| 国产精品麻豆人妻色哟哟久久 | 黄色欧美视频在线观看| 国产中年淑女户外野战色| 在线观看美女被高潮喷水网站| 久久久久久久久久黄片| 18禁黄网站禁片免费观看直播| 22中文网久久字幕| 夫妻性生交免费视频一级片| 亚洲成av人片在线播放无| 最近最新中文字幕大全电影3| 久久久久性生活片| 午夜免费激情av| 99国产精品一区二区蜜桃av| 国内揄拍国产精品人妻在线| 婷婷亚洲欧美| 国内精品久久久久精免费| 国产成人aa在线观看| 听说在线观看完整版免费高清| 亚洲欧美精品专区久久| 一个人看视频在线观看www免费| 又爽又黄无遮挡网站| 毛片一级片免费看久久久久| 精品国产三级普通话版| 亚洲精品乱码久久久久久按摩| 久久鲁丝午夜福利片| 赤兔流量卡办理| 免费看日本二区| 边亲边吃奶的免费视频| 神马国产精品三级电影在线观看| 国产精品一区二区三区四区免费观看| 亚洲性久久影院| 国产精品人妻久久久久久| 日韩大尺度精品在线看网址| 免费黄网站久久成人精品| 99久国产av精品| 干丝袜人妻中文字幕| 亚洲国产精品久久男人天堂| 免费黄网站久久成人精品| 蜜桃久久精品国产亚洲av| 内地一区二区视频在线| 亚洲成a人片在线一区二区| 在线免费观看不下载黄p国产| 99久国产av精品国产电影| 男插女下体视频免费在线播放| 99久久无色码亚洲精品果冻| 久久久久久大精品| 亚洲自偷自拍三级| 丰满的人妻完整版| 精品一区二区免费观看| 免费观看a级毛片全部| 亚洲欧美精品专区久久| 欧美日韩在线观看h| 亚洲精品自拍成人| 两个人的视频大全免费| 亚洲电影在线观看av| 亚洲国产高清在线一区二区三| 人妻制服诱惑在线中文字幕| 中文字幕av在线有码专区| 18+在线观看网站| 小蜜桃在线观看免费完整版高清| 麻豆一二三区av精品| 国产av麻豆久久久久久久| 亚洲国产日韩欧美精品在线观看| 午夜精品一区二区三区免费看| 精品久久久久久久久久久久久| 直男gayav资源| 免费看av在线观看网站| 亚洲精品456在线播放app| 三级毛片av免费| 国产成人aa在线观看| 三级毛片av免费| 午夜激情欧美在线| 国产真实乱freesex| 久久国内精品自在自线图片| 一区二区三区免费毛片| 男人和女人高潮做爰伦理| 亚洲av免费高清在线观看| 99在线视频只有这里精品首页| 午夜激情福利司机影院| 午夜激情福利司机影院| 久久久国产成人精品二区| 美女高潮的动态| 干丝袜人妻中文字幕| 九九爱精品视频在线观看| 99久久人妻综合| 亚洲av电影不卡..在线观看| 99久久人妻综合| 高清在线视频一区二区三区 | 91aial.com中文字幕在线观看| 蜜臀久久99精品久久宅男| 我的女老师完整版在线观看| 少妇的逼好多水| 青春草亚洲视频在线观看| 有码 亚洲区| 午夜亚洲福利在线播放| 亚洲欧美日韩卡通动漫| 丝袜美腿在线中文| 高清毛片免费观看视频网站| 久久久色成人| 韩国av在线不卡| 老司机福利观看| 在线播放国产精品三级| 99热精品在线国产| 中文资源天堂在线| 欧美不卡视频在线免费观看| 日本黄大片高清| 蜜臀久久99精品久久宅男| 国产一区二区在线观看日韩| 五月玫瑰六月丁香| 九九热线精品视视频播放| а√天堂www在线а√下载| 国产成人a∨麻豆精品| 久久久久性生活片| 日韩一区二区视频免费看| 日日摸夜夜添夜夜添av毛片| 国产伦理片在线播放av一区 | 日本黄大片高清| 一区二区三区高清视频在线| 午夜福利在线在线| 婷婷六月久久综合丁香| 久久精品国产清高在天天线| 久久精品综合一区二区三区| 欧美日韩乱码在线| 午夜爱爱视频在线播放| 老司机福利观看| 国产麻豆成人av免费视频| 麻豆成人午夜福利视频| 狂野欧美白嫩少妇大欣赏| 国产精品蜜桃在线观看 | 啦啦啦韩国在线观看视频| 亚洲,欧美,日韩| 99视频精品全部免费 在线| 亚洲一区二区三区色噜噜| 精品人妻视频免费看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲人成网站高清观看| 少妇的逼水好多| 国产精品.久久久| 成人特级黄色片久久久久久久| 欧美日本亚洲视频在线播放| 久久久久久久久大av| 我要看日韩黄色一级片| 99riav亚洲国产免费| 日韩一区二区三区影片| 在线播放国产精品三级| 欧美日韩国产亚洲二区| 在线观看免费视频日本深夜| 亚洲av中文av极速乱| 国产极品天堂在线| 国产精品av视频在线免费观看| 狂野欧美激情性xxxx在线观看| 一边摸一边抽搐一进一小说| 一级黄片播放器| 国内精品久久久久精免费| 国产精品国产高清国产av| 午夜福利高清视频| 日韩欧美精品v在线| 99热网站在线观看| 噜噜噜噜噜久久久久久91| 中国美女看黄片| 最近2019中文字幕mv第一页| 中文字幕免费在线视频6| 国产av麻豆久久久久久久| 国产探花极品一区二区| 久久精品91蜜桃| 精品人妻视频免费看| 变态另类丝袜制服| 夜夜爽天天搞| 国产精品一区www在线观看| 内地一区二区视频在线| 精品久久久久久久久久免费视频| 赤兔流量卡办理| 99久久九九国产精品国产免费| 在线观看av片永久免费下载| 最新中文字幕久久久久| 国产一区亚洲一区在线观看| 国产黄片视频在线免费观看| 能在线免费观看的黄片| 国产在视频线在精品| 亚洲欧美成人精品一区二区| 变态另类成人亚洲欧美熟女| 人人妻人人澡欧美一区二区| 久久久成人免费电影| 亚洲欧美日韩卡通动漫| 精品无人区乱码1区二区| 久久久久久久亚洲中文字幕| 老熟妇乱子伦视频在线观看| 日韩精品青青久久久久久| av免费在线看不卡| 国产黄片视频在线免费观看| 久久精品夜色国产| 国产色婷婷99| 成人亚洲欧美一区二区av| 欧美bdsm另类| 亚洲av一区综合| 99在线视频只有这里精品首页| 亚洲不卡免费看| 久久久久久久亚洲中文字幕| 国产老妇女一区| 亚洲欧洲国产日韩| 久久久精品欧美日韩精品| 一本一本综合久久| 国产乱人偷精品视频| 亚洲,欧美,日韩| av天堂在线播放| 国产精品三级大全| 亚洲av免费高清在线观看| 中文在线观看免费www的网站| kizo精华| 国产乱人视频| 亚洲国产欧洲综合997久久,| 一级毛片我不卡| 亚洲av不卡在线观看| 亚洲国产日韩欧美精品在线观看| 人体艺术视频欧美日本| 亚洲成人中文字幕在线播放| 深爱激情五月婷婷| 观看美女的网站| 自拍偷自拍亚洲精品老妇| 久久久欧美国产精品| 国产精品女同一区二区软件| 国产精品久久久久久av不卡| 女的被弄到高潮叫床怎么办| 美女被艹到高潮喷水动态| 亚洲人成网站高清观看| 免费人成视频x8x8入口观看| 日本免费a在线| 国产老妇伦熟女老妇高清| 亚洲人成网站在线播| 欧美激情国产日韩精品一区| 日韩一本色道免费dvd| 国产成年人精品一区二区| 亚洲精品影视一区二区三区av| 国产69精品久久久久777片| 国产亚洲精品av在线| 久久久午夜欧美精品| 精品午夜福利在线看| 成人无遮挡网站| 国产亚洲欧美98| 国产精品不卡视频一区二区| 久久99精品国语久久久| 99九九线精品视频在线观看视频| 亚洲欧美日韩高清专用| 中文字幕人妻熟人妻熟丝袜美| 色综合色国产| 国产午夜精品久久久久久一区二区三区| 久久久精品94久久精品| 国产成人aa在线观看| 成人午夜精彩视频在线观看| 不卡一级毛片| 国产精品人妻久久久影院| 男人的好看免费观看在线视频| 色尼玛亚洲综合影院| 三级经典国产精品| 日韩制服骚丝袜av| 性插视频无遮挡在线免费观看| 成人亚洲精品av一区二区| 内射极品少妇av片p| 日韩欧美精品v在线| 69人妻影院| 亚洲一区高清亚洲精品| 黄色配什么色好看| 最新中文字幕久久久久| 少妇猛男粗大的猛烈进出视频 | 久久久国产成人免费| 成人综合一区亚洲| 国产爱豆传媒在线观看| 国产一区二区亚洲精品在线观看| 99热6这里只有精品| 不卡视频在线观看欧美| 哪个播放器可以免费观看大片| 丝袜美腿在线中文| 高清在线视频一区二区三区 | 性欧美人与动物交配| 欧美日韩国产亚洲二区| 深夜a级毛片| 69av精品久久久久久| 精品人妻一区二区三区麻豆| 亚洲欧美精品自产自拍| 久久精品影院6| 欧美xxxx性猛交bbbb| 搞女人的毛片| 人人妻人人澡欧美一区二区| av在线播放精品| 中文字幕人妻熟人妻熟丝袜美| av.在线天堂| 一级毛片电影观看 | 欧美精品国产亚洲| 国产麻豆成人av免费视频| av在线亚洲专区| 亚洲欧美成人综合另类久久久 | 日日摸夜夜添夜夜爱| 国产亚洲av片在线观看秒播厂 | 日日撸夜夜添| 亚洲成人久久爱视频| 日本熟妇午夜| 久久亚洲精品不卡| 午夜老司机福利剧场| 国产日韩欧美在线精品| 国产乱人视频| 乱人视频在线观看| 麻豆乱淫一区二区| 国内精品一区二区在线观看| 国产精品无大码| 69av精品久久久久久| 99久久人妻综合| 熟女电影av网| av又黄又爽大尺度在线免费看 | 欧美+亚洲+日韩+国产| 亚洲精品国产成人久久av| 夜夜夜夜夜久久久久| 看免费成人av毛片| 国产私拍福利视频在线观看| 日本一二三区视频观看| 国产激情偷乱视频一区二区| 久久99精品国语久久久| а√天堂www在线а√下载| 99久国产av精品| 亚洲av二区三区四区| 99热网站在线观看| 麻豆国产av国片精品| 中国美白少妇内射xxxbb| 哪个播放器可以免费观看大片| 免费人成视频x8x8入口观看| 国产在视频线在精品| 午夜福利高清视频| 一个人免费在线观看电影| 色尼玛亚洲综合影院| 国产黄片美女视频| 亚洲成人av在线免费| 国产成人91sexporn| 桃色一区二区三区在线观看| 一个人看的www免费观看视频| 色哟哟·www| 蜜桃久久精品国产亚洲av| 欧美一区二区精品小视频在线| 嫩草影院新地址| 丝袜喷水一区| 亚洲成人精品中文字幕电影| 久久久久九九精品影院| 亚洲av中文字字幕乱码综合| 午夜a级毛片| 能在线免费观看的黄片| 色综合色国产| 男女边吃奶边做爰视频| 亚洲美女搞黄在线观看| ponron亚洲| 91久久精品国产一区二区成人| 欧美最黄视频在线播放免费| 日本五十路高清| 亚洲欧美清纯卡通| 亚洲成av人片在线播放无| 国产精品一及| 国产69精品久久久久777片| 99视频精品全部免费 在线| 久久久精品欧美日韩精品| 亚洲精品乱码久久久v下载方式| 精品人妻偷拍中文字幕| 99久久无色码亚洲精品果冻| 午夜免费男女啪啪视频观看| 国产蜜桃级精品一区二区三区| 最近的中文字幕免费完整| 欧美日韩综合久久久久久| 乱人视频在线观看| 乱码一卡2卡4卡精品| 亚洲乱码一区二区免费版| 高清毛片免费观看视频网站| 欧美最新免费一区二区三区| 天天躁日日操中文字幕| 国产精品国产高清国产av| 嫩草影院新地址| 欧美人与善性xxx| 午夜亚洲福利在线播放| 熟妇人妻久久中文字幕3abv| 18禁在线播放成人免费| 伦精品一区二区三区| 欧美性感艳星| 免费人成在线观看视频色| 免费不卡的大黄色大毛片视频在线观看 | 男人和女人高潮做爰伦理| 色5月婷婷丁香| 99在线人妻在线中文字幕| 99热这里只有精品一区| 五月伊人婷婷丁香| 男人狂女人下面高潮的视频| 2021天堂中文幕一二区在线观| 桃色一区二区三区在线观看| 欧美日本亚洲视频在线播放| 日本撒尿小便嘘嘘汇集6| 国产精品福利在线免费观看| 色视频www国产| 美女大奶头视频| 中文字幕制服av| 国产伦精品一区二区三区视频9| 在线播放国产精品三级| 日韩欧美一区二区三区在线观看| 亚洲图色成人| 久久这里只有精品中国| 色视频www国产| 成人毛片60女人毛片免费| 夜夜爽天天搞| av在线亚洲专区| 欧美精品国产亚洲| 蜜桃久久精品国产亚洲av| 久久欧美精品欧美久久欧美| 久久久国产成人免费| 悠悠久久av| 亚洲人与动物交配视频| 亚洲av中文av极速乱| 久久99热6这里只有精品| 日韩欧美精品v在线| 国产精品蜜桃在线观看 | av视频在线观看入口| 黄色一级大片看看| 亚洲精品影视一区二区三区av| 变态另类成人亚洲欧美熟女| 国产成人影院久久av| 成人特级黄色片久久久久久久| 精品一区二区三区视频在线| 两个人的视频大全免费| 一级黄片播放器| 少妇人妻一区二区三区视频| 亚洲久久久久久中文字幕| 久久久成人免费电影| 欧洲精品卡2卡3卡4卡5卡区| 91麻豆精品激情在线观看国产| 国产日韩欧美在线精品| 99久久精品热视频| 十八禁国产超污无遮挡网站| 国产精品一区二区三区四区久久| 成人美女网站在线观看视频| 国产亚洲91精品色在线| 成人欧美大片| 亚洲av免费在线观看| 免费在线观看成人毛片| 22中文网久久字幕| 18禁裸乳无遮挡免费网站照片| 欧美+日韩+精品| 91久久精品国产一区二区三区| 国内精品久久久久精免费| 久久久久久久久大av| 亚洲不卡免费看| 亚洲欧美成人综合另类久久久 | 毛片女人毛片| 九九在线视频观看精品| 午夜老司机福利剧场| h日本视频在线播放| 国产黄片视频在线免费观看| av在线亚洲专区| 91久久精品国产一区二区三区| 一本精品99久久精品77| 国产精品av视频在线免费观看| 国产精品伦人一区二区| 亚洲欧美精品专区久久| 日本黄色片子视频| 欧美日本亚洲视频在线播放| 麻豆av噜噜一区二区三区| 午夜福利高清视频| 午夜爱爱视频在线播放| 国产高潮美女av| 亚洲欧美日韩高清专用| 身体一侧抽搐| 麻豆国产av国片精品| 国产成人精品久久久久久| 精品欧美国产一区二区三| 日本黄色视频三级网站网址| 久久国内精品自在自线图片| 日本-黄色视频高清免费观看| 亚洲精品国产av成人精品| 亚洲人成网站在线播放欧美日韩| 亚洲精品国产av成人精品| 欧美日韩综合久久久久久| 一边摸一边抽搐一进一小说| 国产高清不卡午夜福利| 一级毛片我不卡| 国内精品美女久久久久久| 久久午夜亚洲精品久久| 国产一级毛片七仙女欲春2| 99久国产av精品国产电影| 国产精品一区二区性色av| 欧美色视频一区免费| 欧美精品一区二区大全| 老司机影院成人|