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

    一種適合實(shí)際風(fēng)電接入電力系統(tǒng)的小干擾概率穩(wěn)定分析方法

    2016-11-16 08:39:13楊升峰
    現(xiàn)代電力 2016年5期
    關(guān)鍵詞:分析程序阻尼比出力

    蘇 展,徐 謙,周 明,楊升峰,郭 尊

    (1.新能源電力系統(tǒng)國家重點(diǎn)實(shí)驗(yàn)室(華北電力大學(xué)),北京 102206;2.國網(wǎng)浙江省電力公司經(jīng)濟(jì)技術(shù)研究院,浙江杭州 310000)

    ?

    一種適合實(shí)際風(fēng)電接入電力系統(tǒng)的小干擾概率穩(wěn)定分析方法

    蘇展1,徐謙2,周明1,楊升峰2,郭尊1

    (1.新能源電力系統(tǒng)國家重點(diǎn)實(shí)驗(yàn)室(華北電力大學(xué)),北京102206;2.國網(wǎng)浙江省電力公司經(jīng)濟(jì)技術(shù)研究院,浙江杭州310000)

    0 引 言

    在我國以風(fēng)電為代表的可再生能源發(fā)電正大規(guī)模接入電力系統(tǒng)。風(fēng)電出力的隨機(jī)性對(duì)大規(guī)模風(fēng)電接入電網(wǎng)的小干擾穩(wěn)定的影響凸顯出來[1]。風(fēng)電出力的隨機(jī)變化導(dǎo)致系統(tǒng)運(yùn)行點(diǎn)也具有隨機(jī)性,傳統(tǒng)的確定性分析方法[2]已難以全面分析含大規(guī)模風(fēng)電電力系統(tǒng)的小干擾穩(wěn)定問題,對(duì)電力系統(tǒng)小干擾穩(wěn)定進(jìn)行概率性分析是非常有必要的。

    目前,概率方法在電力系統(tǒng)小干擾穩(wěn)定分析中的應(yīng)用已有報(bào)道,這些研究主要采用蒙特卡羅(Monte Carlo)法[3-5]、以累積量法[2,6,7]為代表的解析法以及以點(diǎn)估計(jì)法[8-9]為代表的近似法分析電力系統(tǒng)的小干擾穩(wěn)定。但由于缺乏實(shí)際數(shù)據(jù)的合適形式以及分析方法受到系統(tǒng)規(guī)模限制等原因,上述研究都只是對(duì)一些小規(guī)模的測試系統(tǒng)進(jìn)行概率分析,對(duì)實(shí)際大系統(tǒng)的適用性不足。使用概率分析來研究風(fēng)電接入的實(shí)際大型電力系統(tǒng)的小干擾穩(wěn)定問題的手段和方法尚未實(shí)現(xiàn)。

    本文針對(duì)大規(guī)模風(fēng)電接入實(shí)際電力系統(tǒng)的穩(wěn)定分析需要,對(duì)于隨機(jī)波動(dòng)的風(fēng)電出力,采用Monte Carlo方法,利用Matlab開發(fā)了以BPA穩(wěn)定性分析程序?yàn)閮?nèi)核的風(fēng)電接入電力系統(tǒng)小干擾概率穩(wěn)定分析系統(tǒng),并用3機(jī)9節(jié)點(diǎn)系統(tǒng)進(jìn)行了風(fēng)電接入容量、接入點(diǎn)及不同的消納方式對(duì)系統(tǒng)小干擾穩(wěn)定性影響的概率分析,然后對(duì)某省電力系統(tǒng)2020年風(fēng)電接入情況的穩(wěn)定性進(jìn)行概率分析。

    1 基于BPA的小干擾概率穩(wěn)定分析

    1.1基于BPA的小干擾穩(wěn)定分析系統(tǒng)構(gòu)成

    所提出的基于BPA的小干擾穩(wěn)定分析方法是通過一個(gè)利用Matlab開發(fā)的風(fēng)電出力Monte Carlo采樣程序、調(diào)用BPA接口程序和結(jié)果統(tǒng)計(jì)分析程序模塊組成的閉環(huán)系統(tǒng)實(shí)現(xiàn)的,其中BPA內(nèi)部的小干擾穩(wěn)定分析程序采用的是特征值分析法[10]。

    BPA中的小干擾穩(wěn)定分析程序只能做確定性的分析,為達(dá)到利用BPA小干擾穩(wěn)定分析程序進(jìn)行概率分析的目的,解決思路是采用循環(huán)調(diào)用執(zhí)行該程序。本文提出用Matlab作為開發(fā)平臺(tái),開發(fā)Matlab自動(dòng)調(diào)用BPA小干擾穩(wěn)定計(jì)算的分析系統(tǒng),實(shí)現(xiàn)風(fēng)電接入電力系統(tǒng)的小干擾概率穩(wěn)定分析。

    1.2BPA與Matlab的數(shù)據(jù)交互

    實(shí)現(xiàn)BPA小干擾穩(wěn)定分析程序與Matlab間的數(shù)據(jù)交互是關(guān)鍵。BPA的小干擾穩(wěn)定分析程序是一個(gè)封裝好的程序,可以將它看作一個(gè)黑盒子,其輸入輸出如圖1所示。

    圖1 BPA小干擾穩(wěn)定分析程序輸入輸出圖

    從圖1可知,該分析程序共有3個(gè)輸入文件,其中潮流文件給出了系統(tǒng)發(fā)電機(jī)出力、網(wǎng)絡(luò)參數(shù)等數(shù)據(jù);穩(wěn)定文件給出了動(dòng)態(tài)元件選用的數(shù)學(xué)模型和詳細(xì)參數(shù);控制文件則決定了分析采用的算法和結(jié)果輸出范圍。而小干擾穩(wěn)定分析結(jié)果則以圖形、文本和表格3種形式給出。由以上分析來看,只要能夠?qū)崿F(xiàn)Matlab與BPA輸入輸出文件的數(shù)據(jù)交互,就可以在Matlab中調(diào)用BPA進(jìn)行電力系統(tǒng)小干擾穩(wěn)定的概率分析。

    在輸入交互上,對(duì)于只考慮風(fēng)電隨機(jī)出力情況,每次計(jì)算的樣本只有包括風(fēng)電機(jī)和同步發(fā)電機(jī)在內(nèi)的各個(gè)機(jī)組的實(shí)際出力是變化的。因此只需每次改變潮流數(shù)據(jù)文件即可,穩(wěn)定文件和控制文件可以保持不變。在輸出交互上,BPA每次的分析結(jié)果可以利用Matlab讀取BPA文本結(jié)果來獲得。

    1.3潮流文件的修改方法

    每次仿真時(shí)都需要對(duì)潮流文件中的風(fēng)電出力和調(diào)節(jié)電源出力進(jìn)行修改。實(shí)際系統(tǒng)中風(fēng)電出力的變化通常由該區(qū)域內(nèi)的火電機(jī)組來平衡。BPA潮流文件中通常都將電源按其類型做歸類處理,這為每次自動(dòng)改變火電出力保證系統(tǒng)潮流平衡帶來了方便。

    BPA中的CHANGES控制語句和PZ卡(發(fā)電出力負(fù)荷百分?jǐn)?shù)修改卡)可以實(shí)現(xiàn)對(duì)同一類型的電源按比例調(diào)節(jié)其有功和無功出力。PZ卡的填寫方式和各數(shù)據(jù)意義如圖2所示。

    圖2 PZ卡形式

    利用上述原則,采用如下方式對(duì)潮流文件進(jìn)行修改。以字符串的形式將整個(gè)初始潮流文件讀入Matlab中記為OrigFile,獲取初始潮流文件中所有風(fēng)電場總出力的原始值Pw_orig,并統(tǒng)計(jì)其中常規(guī)機(jī)組總出力PG。以后每次仿真的潮流文件都在OrigFile的基礎(chǔ)上進(jìn)行修改。對(duì)每一次仿真,首先根據(jù)事先生成的風(fēng)電出力樣本將該次各個(gè)風(fēng)電場出力寫入OrigFile,計(jì)算風(fēng)電場總出力Pw∑,利用式(1)得到發(fā)電有功修改因子k3,最后將k3填入OrigFile中的常規(guī)機(jī)組PZ卡里得到修改后的潮流文件。

    (1)

    1.4小擾動(dòng)概率穩(wěn)定分析算法

    Matlab是一款強(qiáng)大的數(shù)學(xué)分析軟件,其中的system函數(shù)可以對(duì)Matlab以外的可執(zhí)行程序進(jìn)行調(diào)用,并在所調(diào)用的可執(zhí)行程序運(yùn)行完后返回到Matlab中繼續(xù)運(yùn)行Matlab其他代碼。這里就是采用system函數(shù)對(duì)BPA小干擾穩(wěn)定分析程序多次重復(fù)調(diào)用,在Matlab中編寫具體算法如下:

    ① 確定樣本個(gè)數(shù)Ns,利用Matlab工具箱生成各個(gè)風(fēng)電場隨機(jī)出力樣本序列;

    ② 讀取BPA原始潮流數(shù)據(jù)文件并命名為OrigFile,確定原始風(fēng)電場總出力,并令計(jì)數(shù)變量n=1;

    ③ 若n>Ns,轉(zhuǎn)至⑦;若n≤Ns,根據(jù)第n個(gè)風(fēng)電隨機(jī)出力樣本修改OrigFile中的風(fēng)電機(jī)組出力;計(jì)算風(fēng)電機(jī)組出力與原始出力的變化量,進(jìn)一步修改OrigFile中的參與調(diào)節(jié)的其他電源出力。

    ④ 利用修改后OrigFile生成新的BPA潮流數(shù)據(jù)文件。

    ⑤ 利用system函數(shù)調(diào)用BPA小干擾穩(wěn)定分析程序,對(duì)當(dāng)前的BPA潮流數(shù)據(jù)文件做確定性小干擾穩(wěn)定性分析;

    ⑥ 讀取BPA小干擾穩(wěn)定分析程序的文本結(jié)果文件,獲得當(dāng)前樣本的小干擾穩(wěn)定計(jì)算結(jié)果,計(jì)數(shù)n加1,轉(zhuǎn)至③;

    ⑦ 統(tǒng)計(jì)計(jì)算小干擾概率穩(wěn)定指標(biāo)。

    綜上,可將整個(gè)算法流程用圖3表示。

    圖3 小干擾概率穩(wěn)定分析算法流程圖

    2 風(fēng)電隨機(jī)出力樣本的生成方法

    風(fēng)電出力的隨機(jī)性主要來自風(fēng)速的不確定,描述風(fēng)速不確定性常用的概率模型有Weibull、Rayleigh 和Beta分布等[3]。其中,兩參數(shù)的Weibull 分布應(yīng)用最為廣泛。

    本文采用兩參數(shù)的Weibull 分布描述風(fēng)速v的概率分布特性,其概率密度函數(shù)為

    (2)

    式中:c為尺度參數(shù),它表征了平均風(fēng)速的大小;k為形狀參數(shù),它能夠反映風(fēng)速分布的特點(diǎn),對(duì)應(yīng)著Weibull分布的形狀。利用(2)可以得到風(fēng)速的隨機(jī)樣本。在本文中是通過Matlab中的Weibull 分布樣本生成工具箱生成的。

    為了得到風(fēng)電的隨機(jī)出力樣本,還需利用風(fēng)電出力與風(fēng)速之間的相互關(guān)系。工程實(shí)際中,風(fēng)機(jī)并網(wǎng)的最小出力限制和對(duì)其施加的控制措施使得風(fēng)電出力P與風(fēng)速之間存在如圖4所示的近似關(guān)系。圖中vci、vr和vco分別為風(fēng)機(jī)的切入風(fēng)速、額定風(fēng)速和切出風(fēng)速,Pr為風(fēng)電機(jī)組的額定有功出力。a1為風(fēng)機(jī)并網(wǎng)的最小出力占額定出力的比例系數(shù),理論上0≤a1<1。

    圖4 風(fēng)電出力與風(fēng)速關(guān)系圖

    由圖4可得風(fēng)電出力與風(fēng)速的數(shù)學(xué)關(guān)系近似為

    (3)

    式中:a2=1-a1。

    利用式(3)即可將風(fēng)速樣本轉(zhuǎn)換成風(fēng)電的隨機(jī)出力樣本。

    3 小干擾概率穩(wěn)定評(píng)估指標(biāo)

    本文主要考慮以下兩類概率指標(biāo):

    ①系統(tǒng)重要模式的特征值實(shí)部、虛部和阻尼比的均值以及它們的概率分布情況。這些模式主要是那些反映同步發(fā)電機(jī)轉(zhuǎn)子間振蕩頻率在 0.1~2.0Hz 之間的低頻振蕩模式,故稱為機(jī)電振蕩模式。對(duì)于實(shí)際大型電力系統(tǒng),這類振蕩模式的阻尼最小,也最危險(xiǎn),是需要著重關(guān)注的模式[10]。本文主要關(guān)注的也是風(fēng)電對(duì)這類模式的影響。

    ②小干擾穩(wěn)定概率(probability of small signal stability,PSSS)。根據(jù)Lyapunov穩(wěn)定判據(jù),當(dāng)系統(tǒng)狀態(tài)矩陣A的特征值實(shí)部全小于0時(shí),系統(tǒng)在該運(yùn)行點(diǎn)是小干擾穩(wěn)定的,否則系統(tǒng)就是小干擾失穩(wěn)的。因此,以A的特征值實(shí)部最大值λrmax小于0為穩(wěn)定依據(jù),則系統(tǒng)小干擾穩(wěn)定概率可定義為

    (4)

    式中:Nls為λrmax小于0的樣本數(shù);Ns為總樣本數(shù)。

    4 算例分析

    4.1測試系統(tǒng)算例驗(yàn)證

    本文將由雙饋風(fēng)機(jī)構(gòu)成的風(fēng)電場等值后接入BPA自帶的3機(jī)9節(jié)點(diǎn)系統(tǒng)中,對(duì)所開發(fā)的系統(tǒng)進(jìn)行測試,如圖5所示。

    圖5 含風(fēng)電的3機(jī)9節(jié)點(diǎn)系統(tǒng)圖

    表1給出了風(fēng)電接入前系統(tǒng)的機(jī)電振蕩模式(它是根據(jù)振蕩模式的機(jī)電回路相關(guān)比ρ選擇出來的,ρ?1為機(jī)電模式,ρ<1為非機(jī)電模式),可以看出原系統(tǒng)有2個(gè)機(jī)電模式,即表中模式1和模式2。在BPA中對(duì)這兩個(gè)模式進(jìn)一步做模態(tài)分析,發(fā)現(xiàn)模式1表現(xiàn)為機(jī)組G1、G2的功角與機(jī)組G3功角的相對(duì)振蕩,模式2表現(xiàn)為機(jī)組G2、G3的功角與機(jī)組G1功角的相對(duì)振蕩。

    表1 風(fēng)電接入前系統(tǒng)的機(jī)電振蕩模式

    4.1.1風(fēng)電滲透率對(duì)系統(tǒng)小干擾穩(wěn)定性影響

    選用單機(jī)容量為3.6MW的雙饋風(fēng)電機(jī)組接入該系統(tǒng)中的母線6,研究接入25.2MW(對(duì)應(yīng)表2中8%滲透率。當(dāng)前滲透率的定義有很多,本文中滲透率定義為風(fēng)電裝機(jī)容量與系統(tǒng)總負(fù)荷的比值)、50.4MW和100.8MW 3種容量雙饋風(fēng)場對(duì)系統(tǒng)的小干擾穩(wěn)定性的影響。實(shí)際系統(tǒng)中,消納風(fēng)電的方式有3種[1]:①增加旋轉(zhuǎn)備用,減少常規(guī)機(jī)組出力;②用風(fēng)電來滿足新增的負(fù)荷需求;③常規(guī)機(jī)組退出運(yùn)行,即為停機(jī)接納。后兩種接納方式是等價(jià)的,都是通過降低常規(guī)機(jī)組裝機(jī)容量在系統(tǒng)總負(fù)荷中占比實(shí)現(xiàn)對(duì)風(fēng)電的接納。為此,本文在研究滲透率對(duì)系統(tǒng)小干擾穩(wěn)定性的影響中采用如下兩種場景:

    場景1:每次仿真中通過減小發(fā)電機(jī)2和發(fā)電機(jī)3的出力來平衡系統(tǒng)功率。

    場景2:在原有系統(tǒng)基礎(chǔ)上,增加與風(fēng)電裝機(jī)容量相同大小的負(fù)荷,每次仿真中通過依次調(diào)節(jié)發(fā)電機(jī)3和發(fā)電機(jī)1出力來平衡系統(tǒng)功率。

    在場景1下對(duì)上述系統(tǒng)做小干擾概率穩(wěn)定分析,調(diào)用BPA做1 000次計(jì)算模擬風(fēng)電出力的隨機(jī)波動(dòng),得到不同滲透率下系統(tǒng)機(jī)電模式和靠近虛軸的新增振蕩模式的參數(shù)均值,如表2所示,它反映了接入不同滲透率的風(fēng)電時(shí)系統(tǒng)小干擾穩(wěn)定特性的平均水平。對(duì)比表1,從表2中可以看出:①雙饋型風(fēng)場接入后系統(tǒng)的機(jī)電振蕩模式仍為2個(gè)(模式1、2),并未增加與同步發(fā)電機(jī)功角強(qiáng)相關(guān)的振蕩模式。②隨著滲透率的增加,原有機(jī)電模式的阻尼比均值變大。這說明隨滲透率的增大,系統(tǒng)的小干擾穩(wěn)定性總體上是增強(qiáng)的。③雙饋型風(fēng)場接入系統(tǒng)會(huì)出現(xiàn)新的振蕩模式,新出現(xiàn)的模式中存在很靠近虛軸的模式(模式3),意味著增加了系統(tǒng)不穩(wěn)定的風(fēng)險(xiǎn)。

    表2 不同滲透率下系統(tǒng)振蕩模式均值(場景1)

    利用前面提出概率指標(biāo)對(duì)小干擾穩(wěn)定做全局分析,圖6給出了不同滲透率下的1 000次仿真中全部特征值最大實(shí)部和最小阻尼比的概率分布情況,圖6中每一縱列的子圖的橫坐標(biāo)含義和取值完全相同,下文中的其他概率分布圖也如此。可以看出由于風(fēng)電的隨機(jī)出力,即使在同一滲透率下系統(tǒng)小干擾穩(wěn)定特性也有不同且呈現(xiàn)特定的分布,這是確定性分析難以得到的,也恰恰是概率分析的優(yōu)點(diǎn)。從圖中可看出:①在本文場景1中所研究的3種滲透率下,系統(tǒng)特征值的最大實(shí)部都分布在原點(diǎn)左側(cè),說明在1 000次仿真中沒有出現(xiàn)小干擾失穩(wěn)的情況,系統(tǒng)小干擾穩(wěn)定概率指標(biāo)PSSS都為1。②從特征值最小阻尼比的概率分布來看,隨著滲透率的增加,最小阻尼比的分布范圍逐漸增大,這說明,滲透率越大對(duì)系統(tǒng)的小干擾穩(wěn)定影響越大。而分布范圍擴(kuò)大的方向也表明滲透率增大,系統(tǒng)最小阻尼比整體上有向好的趨勢(shì),對(duì)系統(tǒng)小干擾穩(wěn)定性總體上更有利。

    圖6 特征值最大實(shí)部和最小阻尼比的概率分布(場景1)

    在場景2下對(duì)上述系統(tǒng)做小干擾概率穩(wěn)定分析,調(diào)用BPA做1 000次計(jì)算,模擬風(fēng)電出力的隨機(jī)波動(dòng),得到不同滲透率下系統(tǒng)振蕩模式的參數(shù)均值,如表3所示。表3表明,在場景2下隨著滲透率的增加,系統(tǒng)機(jī)電模式的阻尼比均值不斷下降,系統(tǒng)的小干擾穩(wěn)定性總體上是逐漸變差的,這與場景1中的結(jié)果截然相反。

    表3 不同滲透率下系統(tǒng)振蕩模式均值(場景2)

    圖7 特征值最大實(shí)部和最小阻尼比的概率分布(場景2)

    圖7給出的是在場景2中不同滲透率下的1 000次仿真中全部特征值最大實(shí)部和最小阻尼比的概率分布情況。從圖7可以看出:3種滲透率下系統(tǒng)特征值的最大實(shí)部都小于0,說明在1 000次仿真中沒有出現(xiàn)失穩(wěn)的情況,系統(tǒng)小干擾穩(wěn)定概率指標(biāo)PSSS都為1。而根據(jù)系統(tǒng)最小阻尼比的概率分布結(jié)果可知:隨著滲透率增加,最小阻尼比分布逐漸向左移,表明系統(tǒng)小干擾失穩(wěn)的總體風(fēng)險(xiǎn)在加大。

    結(jié)合場景1和場景2的結(jié)果可知,風(fēng)電滲透率對(duì)系統(tǒng)小干擾穩(wěn)定的影響與系統(tǒng)消納風(fēng)電的方式密切相關(guān),當(dāng)通過常規(guī)機(jī)組減出力方式時(shí),滲透率的增加對(duì)系統(tǒng)的小干擾穩(wěn)定性起到改善作用,而采用減少常規(guī)機(jī)組裝機(jī)容量的方式,即風(fēng)電裝機(jī)容量替代常規(guī)機(jī)組時(shí),風(fēng)電接入越多對(duì)小擾動(dòng)穩(wěn)定帶來負(fù)面影響越大。

    4.1.2風(fēng)電接入位置對(duì)系統(tǒng)小干擾穩(wěn)定性影響

    前面的場景1和2是兩種接納風(fēng)電的極端情況,實(shí)際系統(tǒng)中通常是這兩種方式的結(jié)合,因此在分析接入位置時(shí)采用接入50.4MW并新增40MW負(fù)荷的方式,此時(shí)滲透率為14.2%。在此場景下分析了風(fēng)電接入5、6、8三個(gè)負(fù)荷節(jié)點(diǎn)時(shí)系統(tǒng)小干擾穩(wěn)定性,其中節(jié)點(diǎn)5負(fù)荷最大,節(jié)點(diǎn)6負(fù)荷最小。表4給出了各接入點(diǎn)系統(tǒng)振蕩模式的參數(shù)均值,從結(jié)果上看,不同接入位置對(duì)系統(tǒng)振蕩模式的影響是不同的,其中母線6接入雙饋風(fēng)機(jī)對(duì)機(jī)電振蕩模式阻尼比的改善最大。

    圖8給出了不同接入位置下的1 000次仿真中全部特征值最大實(shí)部和最小阻尼比的概率分布情況。從特征值的分布情況可以看出3個(gè)接入位置下系統(tǒng)都沒有出現(xiàn)小干擾失穩(wěn)的情況。而從系統(tǒng)最小阻尼比的概率分布來看,雙饋風(fēng)場接入負(fù)荷最輕的節(jié)點(diǎn)母線6時(shí)系統(tǒng)最小阻尼比分布整體偏右,而接入負(fù)荷最重的節(jié)點(diǎn)母線5的結(jié)果則最靠左,這說明風(fēng)電接入負(fù)荷較輕的節(jié)點(diǎn)對(duì)系統(tǒng)的小干擾穩(wěn)定性更有利。

    表4 不同接入位置下系統(tǒng)部分振蕩模式均值

    圖8 全部特征值最大實(shí)部和最小阻尼比的概率分布

    4.2某省電網(wǎng)實(shí)際算例

    某省的風(fēng)電規(guī)劃主要為海上風(fēng)電,規(guī)劃給出NB地區(qū)的風(fēng)電接入量最大,達(dá)到2 800MW。本文利用所提概率分析方法分析NB地區(qū)接入風(fēng)電后對(duì)整個(gè)省電網(wǎng)的小干擾穩(wěn)定的影響。NB地區(qū)風(fēng)電規(guī)劃如圖9所示。風(fēng)電場1~5經(jīng)過500kV交流匯集點(diǎn)接入句章500kV節(jié)點(diǎn),其接入量分別為700MW、800MW、400MW、300MW和200MW。風(fēng)場6經(jīng)升壓變接入春曉500kV節(jié)點(diǎn),接入量為250MW。風(fēng)場7經(jīng)升壓變接入寧海500kV節(jié)點(diǎn),接入量為150MW。圖9只給出了NB地區(qū)部分節(jié)點(diǎn)。

    圖9 NB地區(qū)風(fēng)電規(guī)劃

    利用BPA對(duì)未接入風(fēng)電的該省電網(wǎng)做小干擾穩(wěn)定性分析得到整個(gè)省電網(wǎng)振蕩模式如圖10所示,其中共有54個(gè)靠近虛軸的機(jī)電振蕩模式,且弱阻尼(阻尼比小于0.05)的只有2個(gè)。

    圖10 風(fēng)電接入前系統(tǒng)振蕩模式分布

    圖11 風(fēng)電接入后1 000次仿真的系統(tǒng)振蕩模式分布

    利用所提方法對(duì)風(fēng)電接入該省電網(wǎng)的小干擾穩(wěn)定進(jìn)行概率分析,圖11為接入風(fēng)電后,做1 000次小干擾穩(wěn)定分析后系統(tǒng)的振蕩模式分布。對(duì)比圖10可以看出,風(fēng)電接入對(duì)系統(tǒng)的振蕩模式并沒有太大影響,但是風(fēng)電接入后該省網(wǎng)系統(tǒng)出現(xiàn)了一些新的靠近原點(diǎn)的振蕩模式,如圖11中的右下角所圈部分所示。分析這些特征值發(fā)現(xiàn)每個(gè)特征值的機(jī)電回路相關(guān)比都小于0.484,表明這些特征值都為非機(jī)電振蕩模式。這說明NB地區(qū)接入風(fēng)電沒有給該系統(tǒng)增加新的功角振蕩問題,但會(huì)增加一些靠近原點(diǎn)的振蕩模式,可能會(huì)對(duì)局部的小干擾穩(wěn)定性造成負(fù)面影響。

    表5給出了風(fēng)電接入前后系統(tǒng)弱阻尼振蕩模式參數(shù)的均值。對(duì)比發(fā)現(xiàn)規(guī)劃的風(fēng)電接入對(duì)系統(tǒng)弱阻尼功角振蕩的影響非常小。

    表5 系統(tǒng)弱阻尼振蕩模式

    圖12是1 000次仿真下,系統(tǒng)全部特征值最大實(shí)部的概率分布情況。它也表明在NB地區(qū)接入風(fēng)電后該系統(tǒng)并未出現(xiàn)失穩(wěn)情況,說明該省電網(wǎng)對(duì)風(fēng)電有很好的接納能力。

    圖12 全部特征值最大實(shí)部的概率分布

    5 結(jié)束語

    本文開發(fā)的基于BPA的風(fēng)電接入電力系統(tǒng)小干擾穩(wěn)定概率型分析方法能夠準(zhǔn)確刻畫隨機(jī)性風(fēng)電功率接入電力系統(tǒng)的小擾動(dòng)穩(wěn)定特征。所開發(fā)的系統(tǒng)可以作為大型實(shí)際電力系統(tǒng)的小干擾概率穩(wěn)定分析工具,對(duì)波動(dòng)性新能源電力系統(tǒng)的規(guī)劃和穩(wěn)定分析也有借鑒作用。

    [1]譚謹(jǐn),王曉茹,李龍?jiān)?含大規(guī)模風(fēng)電的電力系統(tǒng)小擾動(dòng)穩(wěn)定研究綜述[J].電力系統(tǒng)保護(hù)與控制,2014,42(3):15-23.

    [2]邊曉燕,李廣躍,王克文,等.多運(yùn)行方式下含風(fēng)電場電力系統(tǒng)的小干擾概率穩(wěn)定性研究[J].電網(wǎng)技術(shù),2013,37(11):3046-3054.

    [3]李國慶,馬建功,裴建楠,等.基于拉丁超立方采樣的含風(fēng)電小干擾穩(wěn)定概率分析[J].黑龍江電力,2015,37(1):5-10.

    [4]Xu Z, Dong Z Y, Zhang P. Probabilistic small signal analysis using Monte Carlo simulation[C]∥Power Engineering Society General Meeting, 2005. IEEE. IEEE, 2005: 1658-1664.

    [5]Rueda J L, Erlich I. Impacts of large scale integration of wind power on power system small-signal stability[C]∥ Electric Utility Deregulation and Restructuring and Power Technologies (DRPT), 2011 4th International Conference on. IEEE, 2011: 673-681.

    [6]杜文娟,卜思齊,王海風(fēng).考慮并網(wǎng)風(fēng)電隨機(jī)波動(dòng)的電力系統(tǒng)小干擾概率穩(wěn)定性分析[J].中國電機(jī)工程學(xué)報(bào),2011,31(z1):7-11.

    [7]Bu S Q, Du W, Wang H F, et al. Small-signal probabilistic stability of power systems considering the stocastic uncertainty of grid-connected wind farm[C]∥ Renewable Power Generation (RPG 2011), IET Conference on. IET, 2011: 1-6.

    [8]岳昊.考慮并網(wǎng)風(fēng)電隨機(jī)波動(dòng)的電力系統(tǒng)小干擾概率穩(wěn)定研究[D].北京:華北電力大學(xué),2014.

    [9]徐遐齡,查曉明,林濤.微電網(wǎng)小干擾穩(wěn)定概率分析[J].高電壓技術(shù),2009,35 (12):3117-3122.

    [10]仲悟之.大型電力系統(tǒng)小干擾穩(wěn)定性分析方法研究和軟件開發(fā)[D].北京:中國電力科學(xué)研究院,2005.

    (責(zé)任編輯:楊秋霞)

    A Small-signal Probabilistic Stability Analysis Approach for Practical Power System Integrated with Wind Power

    SU Zhan1, XU Qian2, ZHOU Ming1, YANG Shengfeng2, GUO Zun1

    (1.State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources (North China Electric Power University), Beijing 102206, China;2.State Grid Zhejiang Electric Power Company Power Economic Research Institute,Hangzhou 310000, China)

    風(fēng)電出力的隨機(jī)性導(dǎo)致風(fēng)電接入電力系統(tǒng)運(yùn)行點(diǎn)具有明顯的不確定性,也使得確定性小干擾穩(wěn)定分析方法難以全面分析含風(fēng)電的電力系統(tǒng)小干擾穩(wěn)定問題。為此,本文提出了一種實(shí)用化的風(fēng)電接入電力系統(tǒng)小干擾概率穩(wěn)定分析方法。以BPA為小擾動(dòng)穩(wěn)定分析核心程序,通過利用Matlab開發(fā)風(fēng)電出力Monte Carlo采樣程序、調(diào)用BPA接口程序和結(jié)果概率統(tǒng)計(jì)分析程序組成的閉環(huán)系統(tǒng)實(shí)現(xiàn)。用3機(jī)9節(jié)點(diǎn)系統(tǒng)進(jìn)行了風(fēng)電接入容量、接入點(diǎn)及不同的風(fēng)電消納方式對(duì)系統(tǒng)小干擾概率穩(wěn)定的影響分析;然后對(duì)某省實(shí)際電力系統(tǒng)2020年風(fēng)電規(guī)劃情況進(jìn)行了小干擾概率穩(wěn)定分析,驗(yàn)證了所提方法的有效性和實(shí)用性。

    蒙特卡羅;BPA;小干擾穩(wěn)定;概率分析;風(fēng)力發(fā)電

    The randomness of wind power output leads to the uncertainty of power system operating point, which makes it difficult to accurately analyze small signal stability (SSS) of power systems integrated with wind farm. Thus a practical small-signal probabilistic stability analysis approach for power systems integrated with large-scale wind power is proposed. By taking BPA as key procedure for analyzing small-signal stability, the proposed analysis system is designed as a closed loop system, which consist of Matlab based Monte Carlo sampling procedure of wind power output, BPA interface procedure and results probability analysis procedure. The influence of wind power penetration amount, access point and the system accommodating mode on SSS is studied on IEEE 3 machine 9-bus system by using the proposed approach. Then the SSS analysis is tested on certain provincial grid in 2020 year to illustrate the validity and practicability of the proposed approach.

    Monte Carlo; BPA; small signal stability; probabilistic analysis; wind power

    1007-2322(2016)05-0045-07

    A

    TM712

    國家自然科學(xué)基金重大項(xiàng)目(51190103)

    2015-11-21

    蘇展(1990—),男,碩士研究生,主要研究方向?yàn)殡娏ο到y(tǒng)穩(wěn)定性分析,E-mail:1090787474@qq.com;

    徐謙(1963—),男,碩士,高級(jí)工程師,主要研究方向?yàn)殡娏ο到y(tǒng)規(guī)劃和運(yùn)行;

    周明(1967—),女,教授,博士生導(dǎo)師,主要研究方向?yàn)殡娏ο到y(tǒng)優(yōu)化運(yùn)行、電力系統(tǒng)規(guī)劃與可靠性及電力市場等。

    猜你喜歡
    分析程序阻尼比出力
    管控經(jīng)營風(fēng)險(xiǎn),以分析程序提升企業(yè)財(cái)務(wù)報(bào)表審計(jì)效能
    中國商人(2022年9期)2022-09-02 08:30:42
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    黏滯阻尼器在時(shí)程分析下的附加有效阻尼比研究
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    結(jié)構(gòu)構(gòu)件阻尼比對(duì)大跨度懸索橋地震響應(yīng)的影響
    風(fēng)電場有功出力的EEMD特性分析
    要爭做出力出彩的黨員干部
    河南電力(2016年5期)2016-02-06 02:11:35
    基于小波包變換的樂音時(shí)—頻綜合分析程序的開發(fā)
    音樂探索(2015年1期)2015-04-27 01:01:52
    風(fēng)電場群出力的匯聚效應(yīng)分析
    電測與儀表(2015年5期)2015-04-09 11:31:12
    試論分析程序的應(yīng)用
    黑丝袜美女国产一区| 久久久久国产一级毛片高清牌| 久久久久久久久免费视频了| 久久久久国产精品人妻aⅴ院 | 久久久国产成人精品二区 | 国产高清视频在线播放一区| 国产欧美亚洲国产| 99精品久久久久人妻精品| 国产亚洲欧美精品永久| 咕卡用的链子| 91字幕亚洲| 99国产精品免费福利视频| 香蕉丝袜av| 国产精品偷伦视频观看了| 黑人猛操日本美女一级片| 久久精品亚洲熟妇少妇任你| 久久久国产欧美日韩av| 狠狠婷婷综合久久久久久88av| 久久精品亚洲av国产电影网| 国产高清videossex| 不卡av一区二区三区| av天堂在线播放| 国产午夜精品久久久久久| 老鸭窝网址在线观看| 亚洲成a人片在线一区二区| 看片在线看免费视频| 国产成人欧美| 男女床上黄色一级片免费看| 69精品国产乱码久久久| 美女福利国产在线| 少妇粗大呻吟视频| 亚洲成av片中文字幕在线观看| 麻豆av在线久日| 国产精品久久电影中文字幕 | 美女高潮到喷水免费观看| 男人的好看免费观看在线视频 | 麻豆成人av在线观看| 精品国产一区二区三区四区第35| 高清黄色对白视频在线免费看| 免费不卡黄色视频| 啦啦啦视频在线资源免费观看| 日韩中文字幕欧美一区二区| 91麻豆精品激情在线观看国产 | 黑人猛操日本美女一级片| 欧美激情久久久久久爽电影 | 999久久久国产精品视频| 欧美黄色片欧美黄色片| 在线观看午夜福利视频| 久久热在线av| 午夜激情av网站| 涩涩av久久男人的天堂| 国产伦人伦偷精品视频| 国产精品免费大片| 一级黄色大片毛片| 国产蜜桃级精品一区二区三区 | 一进一出抽搐gif免费好疼 | 免费不卡黄色视频| 午夜91福利影院| 久久久久久久国产电影| 欧美激情极品国产一区二区三区| 精品无人区乱码1区二区| 国精品久久久久久国模美| 动漫黄色视频在线观看| 亚洲久久久国产精品| 超色免费av| 黄色丝袜av网址大全| 欧美激情久久久久久爽电影 | 日韩欧美国产一区二区入口| 久久中文字幕一级| 狂野欧美激情性xxxx| 好看av亚洲va欧美ⅴa在| 天天躁夜夜躁狠狠躁躁| 91av网站免费观看| 欧美一级毛片孕妇| 手机成人av网站| 啦啦啦免费观看视频1| 性色av乱码一区二区三区2| 99热网站在线观看| 超碰成人久久| 久热这里只有精品99| 久久精品成人免费网站| 国产精品一区二区免费欧美| 三上悠亚av全集在线观看| www.999成人在线观看| 国产熟女午夜一区二区三区| 国产成人av激情在线播放| 丰满饥渴人妻一区二区三| 热99久久久久精品小说推荐| tocl精华| 咕卡用的链子| 久久精品熟女亚洲av麻豆精品| 亚洲av熟女| svipshipincom国产片| 国产乱人伦免费视频| 久久狼人影院| 久久久久国内视频| 99国产精品99久久久久| 极品少妇高潮喷水抽搐| 国产成人精品久久二区二区91| 成熟少妇高潮喷水视频| 国产亚洲av高清不卡| 少妇被粗大的猛进出69影院| 一边摸一边抽搐一进一小说 | 国产男女超爽视频在线观看| 欧美一级毛片孕妇| 一级黄色大片毛片| 欧美黑人欧美精品刺激| 热99re8久久精品国产| 91九色精品人成在线观看| 新久久久久国产一级毛片| 国产av精品麻豆| 亚洲免费av在线视频| 高清欧美精品videossex| 在线视频色国产色| 一级,二级,三级黄色视频| 成人18禁在线播放| 亚洲国产欧美日韩在线播放| 国产亚洲精品久久久久5区| 少妇的丰满在线观看| 不卡av一区二区三区| 久久亚洲精品不卡| 国产精品久久视频播放| 免费在线观看日本一区| 午夜视频精品福利| 国产精品av久久久久免费| 成人国语在线视频| 麻豆成人av在线观看| 美女高潮喷水抽搐中文字幕| 日本一区二区免费在线视频| 免费在线观看视频国产中文字幕亚洲| 一进一出好大好爽视频| 极品人妻少妇av视频| 国产精品久久久久成人av| 精品国产一区二区三区四区第35| 一区二区三区精品91| 中文字幕av电影在线播放| 人人妻,人人澡人人爽秒播| 亚洲九九香蕉| 黑人欧美特级aaaaaa片| 婷婷成人精品国产| 香蕉丝袜av| 校园春色视频在线观看| 午夜精品国产一区二区电影| 国产有黄有色有爽视频| 伦理电影免费视频| 国产在线一区二区三区精| 国内毛片毛片毛片毛片毛片| 国产精品永久免费网站| 三上悠亚av全集在线观看| 91av网站免费观看| 淫妇啪啪啪对白视频| 国产欧美日韩一区二区三区在线| 18禁美女被吸乳视频| 多毛熟女@视频| 午夜两性在线视频| 香蕉丝袜av| 丰满人妻熟妇乱又伦精品不卡| 建设人人有责人人尽责人人享有的| 亚洲熟妇中文字幕五十中出 | 精品国产国语对白av| 国产在线观看jvid| 久久九九热精品免费| 日日爽夜夜爽网站| 久久精品aⅴ一区二区三区四区| 国产成人影院久久av| 高潮久久久久久久久久久不卡| 亚洲全国av大片| 性色av乱码一区二区三区2| 国内久久婷婷六月综合欲色啪| 老汉色∧v一级毛片| 久久精品人人爽人人爽视色| www.熟女人妻精品国产| 亚洲少妇的诱惑av| 亚洲精品自拍成人| 一级毛片高清免费大全| 国产又色又爽无遮挡免费看| 国产极品粉嫩免费观看在线| av免费在线观看网站| 男女床上黄色一级片免费看| 精品欧美一区二区三区在线| 天堂俺去俺来也www色官网| 欧美激情极品国产一区二区三区| 首页视频小说图片口味搜索| 国产亚洲av高清不卡| 午夜福利免费观看在线| 男女高潮啪啪啪动态图| 国产在线一区二区三区精| cao死你这个sao货| 精品国产国语对白av| 国产主播在线观看一区二区| bbb黄色大片| 色婷婷久久久亚洲欧美| 一二三四社区在线视频社区8| 老司机午夜十八禁免费视频| 国产精品98久久久久久宅男小说| 亚洲av日韩精品久久久久久密| 国产男女内射视频| 91字幕亚洲| 欧美在线一区亚洲| 黑人操中国人逼视频| 国产淫语在线视频| 操出白浆在线播放| 老司机深夜福利视频在线观看| 亚洲男人天堂网一区| 操美女的视频在线观看| 在线观看www视频免费| 天天添夜夜摸| 中文字幕色久视频| 欧美日韩亚洲高清精品| 满18在线观看网站| 热re99久久国产66热| 欧美在线黄色| bbb黄色大片| 欧美精品亚洲一区二区| 亚洲熟妇熟女久久| 日韩欧美在线二视频 | 久久久久视频综合| 精品国产亚洲在线| 啪啪无遮挡十八禁网站| 美国免费a级毛片| 黄色女人牲交| 国产单亲对白刺激| 乱人伦中国视频| 久久亚洲真实| 亚洲欧美一区二区三区久久| 国产欧美亚洲国产| 中文字幕最新亚洲高清| tube8黄色片| 大香蕉久久成人网| 热99国产精品久久久久久7| 黄片小视频在线播放| 热99re8久久精品国产| 99精品在免费线老司机午夜| 亚洲片人在线观看| 久久午夜亚洲精品久久| 一级毛片精品| 人妻 亚洲 视频| 一进一出抽搐动态| 9热在线视频观看99| 亚洲精品国产精品久久久不卡| 一级片'在线观看视频| 欧美日韩黄片免| 91av网站免费观看| 亚洲午夜理论影院| 亚洲av美国av| 国产激情久久老熟女| 大型黄色视频在线免费观看| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看 | 免费日韩欧美在线观看| 五月开心婷婷网| 黑人猛操日本美女一级片| 乱人伦中国视频| 丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av日韩精品久久久久久密| 夫妻午夜视频| 午夜福利视频在线观看免费| 亚洲色图 男人天堂 中文字幕| 老司机深夜福利视频在线观看| 亚洲av成人av| 国产成人啪精品午夜网站| 日韩欧美国产一区二区入口| 麻豆国产av国片精品| 手机成人av网站| 午夜福利免费观看在线| 看免费av毛片| 精品人妻在线不人妻| 操美女的视频在线观看| 18禁观看日本| 久久精品亚洲av国产电影网| 亚洲精品成人av观看孕妇| 天天添夜夜摸| 91精品国产国语对白视频| 国产欧美日韩一区二区三| 99热网站在线观看| 国产一区二区三区视频了| av超薄肉色丝袜交足视频| 亚洲国产欧美网| 宅男免费午夜| 丰满人妻熟妇乱又伦精品不卡| 国产男靠女视频免费网站| 中国美女看黄片| 在线十欧美十亚洲十日本专区| 亚洲精品美女久久av网站| 成熟少妇高潮喷水视频| 91成人精品电影| 99香蕉大伊视频| 一本综合久久免费| 黄片小视频在线播放| e午夜精品久久久久久久| 亚洲熟女毛片儿| 满18在线观看网站| 757午夜福利合集在线观看| 亚洲av成人av| 国产在线精品亚洲第一网站| 他把我摸到了高潮在线观看| 老司机福利观看| 国产日韩欧美亚洲二区| 99热国产这里只有精品6| 欧美成人午夜精品| 免费人成视频x8x8入口观看| 国产亚洲精品一区二区www | 九色亚洲精品在线播放| 亚洲欧美一区二区三区黑人| 搡老熟女国产l中国老女人| 久久久久精品人妻al黑| 人人澡人人妻人| 国产精品香港三级国产av潘金莲| 老司机福利观看| 日本五十路高清| 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看影片大全网站| 欧美乱妇无乱码| 黄色成人免费大全| 免费看十八禁软件| av免费在线观看网站| 男女免费视频国产| 王馨瑶露胸无遮挡在线观看| 亚洲熟妇中文字幕五十中出 | 欧美日韩黄片免| 色精品久久人妻99蜜桃| 97人妻天天添夜夜摸| 亚洲av熟女| 国产亚洲精品一区二区www | www日本在线高清视频| 1024香蕉在线观看| 久久ye,这里只有精品| 国产不卡一卡二| 亚洲成a人片在线一区二区| 18禁观看日本| 午夜老司机福利片| 伊人久久大香线蕉亚洲五| 久久性视频一级片| 日韩中文字幕欧美一区二区| 高清毛片免费观看视频网站 | 午夜视频精品福利| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲三区欧美一区| 久久精品91无色码中文字幕| 一级a爱视频在线免费观看| а√天堂www在线а√下载 | 高清毛片免费观看视频网站 | 飞空精品影院首页| 久久久国产成人免费| 日本黄色视频三级网站网址 | 国产精品综合久久久久久久免费 | 日本精品一区二区三区蜜桃| 午夜久久久在线观看| 精品视频人人做人人爽| 国产成人啪精品午夜网站| 精品视频人人做人人爽| 色尼玛亚洲综合影院| 一a级毛片在线观看| 很黄的视频免费| 国产精品影院久久| av片东京热男人的天堂| 涩涩av久久男人的天堂| 老司机午夜福利在线观看视频| 十八禁高潮呻吟视频| 欧美日韩福利视频一区二区| 亚洲av成人不卡在线观看播放网| bbb黄色大片| 一进一出抽搐gif免费好疼 | netflix在线观看网站| 国产亚洲欧美在线一区二区| 免费一级毛片在线播放高清视频 | 午夜福利欧美成人| 人妻久久中文字幕网| 国产成人系列免费观看| 女人精品久久久久毛片| 亚洲欧美色中文字幕在线| 最新在线观看一区二区三区| 亚洲国产精品sss在线观看 | 中文字幕最新亚洲高清| 1024香蕉在线观看| 一级毛片女人18水好多| 99在线人妻在线中文字幕 | 狠狠狠狠99中文字幕| 国产精品久久久久成人av| 一夜夜www| 一本综合久久免费| 亚洲精品中文字幕一二三四区| 欧美国产精品va在线观看不卡| 极品教师在线免费播放| 国产免费男女视频| 18禁美女被吸乳视频| 免费在线观看视频国产中文字幕亚洲| 人成视频在线观看免费观看| 午夜激情av网站| 国产精品 国内视频| 亚洲五月色婷婷综合| 国产精品久久视频播放| 亚洲中文字幕日韩| 欧美日韩福利视频一区二区| 乱人伦中国视频| 久久久久久亚洲精品国产蜜桃av| 国产在线观看jvid| 亚洲国产中文字幕在线视频| 热99久久久久精品小说推荐| 人人澡人人妻人| 日韩三级视频一区二区三区| 热re99久久国产66热| 露出奶头的视频| netflix在线观看网站| 99精国产麻豆久久婷婷| 欧美日韩瑟瑟在线播放| 麻豆av在线久日| 黄网站色视频无遮挡免费观看| 国产精华一区二区三区| 女人被躁到高潮嗷嗷叫费观| 精品无人区乱码1区二区| 精品人妻1区二区| 久久久久久亚洲精品国产蜜桃av| 国产一区在线观看成人免费| www.精华液| 久久天堂一区二区三区四区| 高潮久久久久久久久久久不卡| 日本黄色日本黄色录像| 韩国av一区二区三区四区| 国产一区二区三区在线臀色熟女 | 无遮挡黄片免费观看| 十分钟在线观看高清视频www| 女人高潮潮喷娇喘18禁视频| 777米奇影视久久| av欧美777| 成人18禁在线播放| 久久久久精品人妻al黑| 色播在线永久视频| 国产成人精品久久二区二区91| 精品福利观看| 18禁裸乳无遮挡免费网站照片 | a级毛片在线看网站| 久久国产精品影院| 香蕉国产在线看| 中文字幕制服av| 精品一区二区三区四区五区乱码| 欧美+亚洲+日韩+国产| 亚洲成a人片在线一区二区| 亚洲欧美色中文字幕在线| 日韩欧美在线二视频 | 亚洲av成人一区二区三| 久久人妻福利社区极品人妻图片| 色老头精品视频在线观看| 日韩大码丰满熟妇| 欧美亚洲日本最大视频资源| 日韩有码中文字幕| 丁香欧美五月| 黑人猛操日本美女一级片| 国产精品二区激情视频| 在线国产一区二区在线| cao死你这个sao货| 国产亚洲一区二区精品| av福利片在线| 亚洲欧美色中文字幕在线| 国产成人啪精品午夜网站| 黄片小视频在线播放| 丁香六月欧美| 极品教师在线免费播放| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩成人在线一区二区| 国产成人欧美| 久久久国产成人精品二区 | 精品福利观看| 欧美亚洲 丝袜 人妻 在线| 妹子高潮喷水视频| 很黄的视频免费| 亚洲欧美色中文字幕在线| 多毛熟女@视频| 一二三四社区在线视频社区8| 男男h啪啪无遮挡| 中文字幕另类日韩欧美亚洲嫩草| 麻豆成人av在线观看| 欧美色视频一区免费| 新久久久久国产一级毛片| 久久久国产一区二区| 亚洲成人国产一区在线观看| 亚洲中文日韩欧美视频| 男女之事视频高清在线观看| 日韩中文字幕欧美一区二区| 午夜老司机福利片| 久久亚洲精品不卡| 亚洲精品久久午夜乱码| 岛国毛片在线播放| 亚洲va日本ⅴa欧美va伊人久久| 在线看a的网站| 久热爱精品视频在线9| 亚洲熟女精品中文字幕| 99国产极品粉嫩在线观看| 国产精品久久电影中文字幕 | 国内久久婷婷六月综合欲色啪| 一本一本久久a久久精品综合妖精| 少妇猛男粗大的猛烈进出视频| 少妇裸体淫交视频免费看高清 | 久久精品国产a三级三级三级| 中出人妻视频一区二区| 免费在线观看完整版高清| 亚洲精品av麻豆狂野| 久久精品国产亚洲av高清一级| 一边摸一边做爽爽视频免费| 视频区图区小说| av在线播放免费不卡| 国产高清国产精品国产三级| 亚洲人成77777在线视频| 色老头精品视频在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲美女黄片视频| 国产日韩一区二区三区精品不卡| 黄色视频,在线免费观看| 水蜜桃什么品种好| 99香蕉大伊视频| 美女午夜性视频免费| 欧美大码av| 久久精品人人爽人人爽视色| 美女国产高潮福利片在线看| 99精品久久久久人妻精品| 欧美av亚洲av综合av国产av| 亚洲在线自拍视频| 久久国产精品大桥未久av| 亚洲专区国产一区二区| 女警被强在线播放| 日韩精品免费视频一区二区三区| 亚洲欧美日韩高清在线视频| 国产伦人伦偷精品视频| 亚洲精品国产区一区二| 女人精品久久久久毛片| 日韩欧美三级三区| 国产亚洲精品一区二区www | 亚洲美女黄片视频| 国产aⅴ精品一区二区三区波| 午夜福利免费观看在线| 麻豆成人av在线观看| 视频区图区小说| 日韩欧美一区视频在线观看| 啦啦啦在线免费观看视频4| 日本a在线网址| 久久久精品区二区三区| 99久久精品国产亚洲精品| 国产精品久久久人人做人人爽| 欧美最黄视频在线播放免费 | 一区在线观看完整版| 久久狼人影院| 国产在线精品亚洲第一网站| 国产伦人伦偷精品视频| 精品一区二区三区视频在线观看免费 | 欧美精品人与动牲交sv欧美| 一进一出抽搐gif免费好疼 | 欧美成人午夜精品| aaaaa片日本免费| 久久中文看片网| 王馨瑶露胸无遮挡在线观看| 99riav亚洲国产免费| 午夜免费成人在线视频| 中文字幕最新亚洲高清| 国产av精品麻豆| 激情视频va一区二区三区| aaaaa片日本免费| av欧美777| 国产成人免费无遮挡视频| 亚洲精华国产精华精| 纯流量卡能插随身wifi吗| 欧美另类亚洲清纯唯美| 国产三级黄色录像| 欧美精品亚洲一区二区| 热99re8久久精品国产| 久久天堂一区二区三区四区| 久久精品国产综合久久久| 18禁裸乳无遮挡免费网站照片 | 18禁黄网站禁片午夜丰满| 国产高清激情床上av| 99久久99久久久精品蜜桃| 亚洲一区二区三区欧美精品| av中文乱码字幕在线| 最近最新中文字幕大全电影3 | 亚洲中文日韩欧美视频| 久久精品91无色码中文字幕| av天堂在线播放| 亚洲在线自拍视频| 亚洲国产精品一区二区三区在线| 亚洲专区中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 极品教师在线免费播放| 在线观看www视频免费| 亚洲综合色网址| 国产亚洲欧美在线一区二区| 亚洲精品成人av观看孕妇| 女人精品久久久久毛片| 日本黄色日本黄色录像| 成年版毛片免费区| 日韩一卡2卡3卡4卡2021年| 丰满的人妻完整版| 咕卡用的链子| 如日韩欧美国产精品一区二区三区| avwww免费| 看片在线看免费视频| 高清毛片免费观看视频网站 | 人妻一区二区av| 久久精品亚洲熟妇少妇任你| 国产成人av教育| 18禁国产床啪视频网站| 欧美日韩瑟瑟在线播放| 国产国语露脸激情在线看| 一进一出好大好爽视频| 又黄又粗又硬又大视频| av福利片在线| 在线永久观看黄色视频| 老司机深夜福利视频在线观看| 午夜精品久久久久久毛片777| 国产精品国产av在线观看| 久久精品亚洲精品国产色婷小说| 国产熟女午夜一区二区三区| 一区二区日韩欧美中文字幕| 国产在视频线精品| 日韩大码丰满熟妇| 丁香欧美五月|