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

    考慮不確定因素的關(guān)鍵輸電斷面分析

    2016-12-15 03:22:24靳冰潔梁晨
    廣東電力 2016年11期
    關(guān)鍵詞:介數(shù)支路分區(qū)

    靳冰潔,梁晨

    (1.廣東電網(wǎng)發(fā)展研究院有限責(zé)任公司,廣東 廣州 510080;2.中國南方電網(wǎng)超高壓輸電公司檢修試驗(yàn)中心,廣東 廣州 510663)

    ?

    考慮不確定因素的關(guān)鍵輸電斷面分析

    靳冰潔1,梁晨2

    (1.廣東電網(wǎng)發(fā)展研究院有限責(zé)任公司,廣東 廣州 510080;2.中國南方電網(wǎng)超高壓輸電公司檢修試驗(yàn)中心,廣東 廣州 510663)

    計(jì)及負(fù)荷與間歇式電源波動在網(wǎng)絡(luò)中的傳播特性,提出基于綜合介數(shù)指標(biāo)的系統(tǒng)關(guān)鍵支路辨識方法,用以全面反映特定運(yùn)行環(huán)境、網(wǎng)絡(luò)結(jié)構(gòu)及運(yùn)行方式下系統(tǒng)的薄弱環(huán)節(jié)。然后提出基于節(jié)點(diǎn)電氣距離聚類的電網(wǎng)分區(qū)算法及其修正措施,避免漏選區(qū)內(nèi)重要線路。最后,在電網(wǎng)分區(qū)的基礎(chǔ)上根據(jù)輸電斷面特征選擇關(guān)鍵支路主導(dǎo)的輸電斷面作為系統(tǒng)關(guān)鍵輸電斷面。IEEE 39節(jié)點(diǎn)系統(tǒng)算例驗(yàn)證了該方法的有效性。實(shí)際應(yīng)用時(shí),可根據(jù)不同時(shí)段的負(fù)荷或間歇式能源波動特性以及系統(tǒng)運(yùn)行方式的變化隨時(shí)搜索系統(tǒng)關(guān)鍵輸電斷面,為有針對性地采取系統(tǒng)安全控制措施提供依據(jù)。

    綜合介數(shù);關(guān)鍵支路;電網(wǎng)分區(qū);輸電斷面

    現(xiàn)代電力系統(tǒng)具有規(guī)模龐大、元件復(fù)雜、隨機(jī)因素眾多等特性,給運(yùn)行研究帶來了困難和挑戰(zhàn)[1]。輸電斷面作為大電網(wǎng)的重要安全特征而受到運(yùn)行人員的廣泛關(guān)注[2-5]。我國地域遼闊,能源分布不均勻,地區(qū)發(fā)展不平衡,輸電線路承擔(dān)著大規(guī)模、遠(yuǎn)距離送電的重要任務(wù),容易成為電網(wǎng)運(yùn)行的薄弱環(huán)節(jié)。關(guān)鍵輸電斷面的識別與分析,有利于運(yùn)行人員掌握電網(wǎng)安全運(yùn)行的重要信息,有針對性地實(shí)施安全保護(hù)措施。

    負(fù)荷性質(zhì)的多樣化和間歇式能源的接入帶來許多不確定因素,時(shí)刻威脅著系統(tǒng)的安全穩(wěn)定運(yùn)行。系統(tǒng)中各線路受負(fù)荷及間歇式電源波動特性的影響程度不同,將導(dǎo)致關(guān)鍵輸電斷面的識別結(jié)果與確定性環(huán)境下有所差異。因此,關(guān)鍵輸電斷面分析應(yīng)計(jì)及不確定因素帶來的影響?,F(xiàn)有輸電斷面識別方法大多僅考慮電網(wǎng)結(jié)構(gòu)特征及運(yùn)行方式,對含不確定因素環(huán)境下的輸電斷面特征分析鮮有論及。

    針對傳統(tǒng)基于運(yùn)行方式專家經(jīng)驗(yàn)的輸電斷面選取方法難以適應(yīng)運(yùn)行方式變化的問題,已有很多文獻(xiàn)提出輸電斷面自動識別方法。電網(wǎng)分區(qū)是發(fā)現(xiàn)關(guān)鍵輸電斷面的有效手段[6-7],能夠簡化網(wǎng)絡(luò),降低全網(wǎng)搜索的工作量。然而現(xiàn)有分區(qū)方法大多基于某一特定原則,缺乏針對誤差的修正措施,容易存在漏選區(qū)內(nèi)關(guān)鍵線路的情況?;谶^載或故障支路選擇相關(guān)輸電斷面[8-9]是從電網(wǎng)安全控制角度分析問題的可行思路,此類方法主要依賴關(guān)鍵支路的選取,在通過關(guān)鍵支路的選取來甄別系統(tǒng)薄弱環(huán)節(jié)方面仍有待進(jìn)一步深入研究。

    本文在以上研究內(nèi)容的基礎(chǔ)上提出考慮不確定因素的關(guān)鍵輸電斷面識別方法。計(jì)及負(fù)荷與間歇式電源波動在網(wǎng)絡(luò)中的傳播特性,提出基于綜合介數(shù)指標(biāo)的系統(tǒng)關(guān)鍵支路辨識方法;在此基礎(chǔ)上提出基于節(jié)點(diǎn)電氣距離聚類的電網(wǎng)分區(qū)算法及其修正措施,避免漏選區(qū)內(nèi)重要線路;最后,在電網(wǎng)分區(qū)的基礎(chǔ)上根據(jù)輸電斷面特征選擇關(guān)鍵支路主導(dǎo)的輸電斷面作為系統(tǒng)關(guān)鍵輸電斷面。值得指出的是,目前文獻(xiàn)對輸電斷面的細(xì)節(jié)定義仍有差異,如一些文獻(xiàn)認(rèn)為輸電斷面是電網(wǎng)中的割集[6-7],而另一些文獻(xiàn)不這樣認(rèn)為[8-9],本文采取前者說法。

    1 考慮不確定因素的關(guān)鍵支路辨識

    輸電斷面由電網(wǎng)中的部分支路組成,能反映電網(wǎng)中的薄弱環(huán)節(jié)[8],對電力系統(tǒng)進(jìn)行關(guān)鍵支路辨識能夠發(fā)現(xiàn)系統(tǒng)中的薄弱支路,進(jìn)而引導(dǎo)輸電斷面的選擇?;陉P(guān)鍵支路選擇的輸電斷面能夠集中反映電力系統(tǒng)中易受攻擊或干擾而故障,或者故障后嚴(yán)重影響電能傳輸并直接威脅電網(wǎng)安全的薄弱環(huán)節(jié)。

    介數(shù)指標(biāo)已被廣泛應(yīng)用于電力系統(tǒng)關(guān)鍵線路辨識中[7,10-13]。傳統(tǒng)介數(shù)方法通常基于最短路徑傳輸?shù)募僭O(shè),無法準(zhǔn)確反映電力系統(tǒng)實(shí)際運(yùn)行狀態(tài);而加權(quán)介數(shù)[11]加雖考慮了電網(wǎng)自身拓?fù)浣Y(jié)構(gòu)的特點(diǎn),但未反應(yīng)系統(tǒng)當(dāng)前運(yùn)行狀態(tài);一些文獻(xiàn)將網(wǎng)絡(luò)結(jié)構(gòu)與運(yùn)行狀態(tài)相結(jié)合,提出改進(jìn)方法,如電氣介數(shù)[12]、輸電介數(shù)[7]、潮流介數(shù)[13]等。然而,以上方法均未考慮不確定性因素對系統(tǒng)安全運(yùn)行的影響。

    負(fù)荷或間歇式能源的波動性會隨其在電網(wǎng)中的傳播而威脅系統(tǒng)安全運(yùn)行,而系統(tǒng)中各支路受其影響的程度不同。因此,關(guān)鍵支路的辨識應(yīng)在網(wǎng)架結(jié)構(gòu)及運(yùn)行方式的基礎(chǔ)上,計(jì)及不確定因素在網(wǎng)絡(luò)中傳播帶來的影響?;诖?,本文提出線路的綜合介數(shù)指標(biāo)Bc,其計(jì)算公式為

    (1)

    式中:i為發(fā)電節(jié)點(diǎn),G為發(fā)電節(jié)點(diǎn)集合;j為負(fù)荷節(jié)點(diǎn),D為負(fù)荷節(jié)點(diǎn)集合;eδi、eδj為反應(yīng)發(fā)電節(jié)點(diǎn)i和負(fù)荷節(jié)點(diǎn)j功率誤差傳播影響的誤差系數(shù),δi、δj分別為換算到基準(zhǔn)容量下的發(fā)電節(jié)點(diǎn)i和負(fù)荷節(jié)點(diǎn)j功率誤差百分比;Gk.ij為發(fā)電-負(fù)荷節(jié)點(diǎn)對ij對支路k的功率傳輸轉(zhuǎn)移分布因子(power transmission distribution factors,PTDF)[14],反映節(jié)點(diǎn)對之間功率傳輸對線路的占用情況;Cij為發(fā)電-負(fù)荷節(jié)點(diǎn)對ij的最大傳輸容量。

    (2)

    (3)

    式中:m、n為支路k兩端節(jié)點(diǎn)編號;Xmi、Xmj、Xni和Xnj均為系統(tǒng)電納矩陣的逆矩陣X中的元素;xk為支路k的電抗。

    (4)

    發(fā)電-負(fù)荷節(jié)點(diǎn)對之間的最大傳輸容量受限于發(fā)電節(jié)點(diǎn)功率和負(fù)荷節(jié)點(diǎn)功率中較小的值,反應(yīng)了系統(tǒng)的發(fā)電和負(fù)荷水平對節(jié)點(diǎn)對ij傳輸容量的限制,因此將Cij作為發(fā)電-負(fù)荷節(jié)點(diǎn)對ij功率傳輸轉(zhuǎn)移分布因子的傳輸容量權(quán)重。

    綜合介數(shù)指標(biāo)綜合考慮了電網(wǎng)的結(jié)構(gòu)特征、運(yùn)行方式以及系統(tǒng)中不確定因素的影響,其中,通過功率傳輸轉(zhuǎn)移分布因子反映支路在系統(tǒng)各發(fā)電-負(fù)荷節(jié)點(diǎn)對之間功率傳輸?shù)碾娋W(wǎng)結(jié)構(gòu)特征重要性,通過傳輸容量權(quán)重考慮當(dāng)前運(yùn)行方式下發(fā)電和負(fù)荷水平對節(jié)點(diǎn)對傳輸容量的限制,通過誤差系數(shù)計(jì)及發(fā)電-負(fù)荷節(jié)點(diǎn)對功率傳輸過程中分布式電源和負(fù)荷波動在網(wǎng)絡(luò)中傳播給支路帶來的影響。

    基于綜合介數(shù)的關(guān)鍵支路辨識應(yīng)用于電力系統(tǒng)關(guān)鍵輸電斷面選取主要體現(xiàn)在兩方面:一方面,將支路綜合介數(shù)指標(biāo)按照從大到小的順序排列,選擇綜合介數(shù)最大的若干條線路作為主導(dǎo)線路,并選取其主導(dǎo)的輸電斷面作為關(guān)鍵輸電斷面;另一方面,選擇Bc>ε1(ε1為預(yù)先設(shè)定的關(guān)鍵線路綜合介數(shù)閾值)的線路形成關(guān)鍵線路集Lkey,用以對電網(wǎng)分區(qū)算法進(jìn)行修正,使得關(guān)鍵支路被包含于分區(qū)間聯(lián)絡(luò)線中,而不至于在選擇輸電斷面時(shí)被忽略。

    2 基于節(jié)點(diǎn)聚類算法的電網(wǎng)分區(qū)

    輸電斷面的選取方法很多,其本質(zhì)仍是對電網(wǎng)結(jié)構(gòu)特征的辨識過程[7]。對電網(wǎng)進(jìn)行分區(qū),在分區(qū)聯(lián)絡(luò)線中選擇滿足輸電斷面特征的線路集合,是對電力系統(tǒng)的一種有效降維,與在全網(wǎng)所有線路中搜索輸電斷面相比,工作量大為減少,利于實(shí)際應(yīng)用。

    已有的基于系統(tǒng)分區(qū)的輸電斷面選取算法中,文獻(xiàn)[6]以地理分區(qū)作為電網(wǎng)初始分區(qū),對系統(tǒng)運(yùn)行方式的適應(yīng)性稍差;文獻(xiàn)[7]基于分裂類算法對電網(wǎng)進(jìn)行劃分,雖然能夠保證關(guān)鍵線路在分區(qū)之間,但無法保證分區(qū)內(nèi)的節(jié)點(diǎn)之間相互聯(lián)系緊密。電網(wǎng)分區(qū)既要保證分區(qū)內(nèi)節(jié)點(diǎn)之間相互聯(lián)系緊密,又要避免將關(guān)鍵線路和長程線路[7]劃分區(qū)內(nèi)而在選取輸電斷面時(shí)將其忽略。為此,有必要對基于聚類算法的分區(qū)結(jié)果進(jìn)行修正,使其同時(shí)滿足以上兩點(diǎn)要求。

    綜上,本文提出基于節(jié)點(diǎn)聚類算法及其修正方法的電網(wǎng)分區(qū)方法,分為初步分區(qū)和最終分區(qū)兩個(gè)過程。其中,初步分區(qū)根據(jù)節(jié)點(diǎn)之間的電氣距離進(jìn)行聚類,最終分區(qū)根據(jù)系統(tǒng)關(guān)鍵線路集和長程線路集對初步分區(qū)結(jié)果進(jìn)行修正。

    2.1 電網(wǎng)初步分區(qū)

    基于電氣距離的系統(tǒng)節(jié)點(diǎn)聚類算法既能體現(xiàn)網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)特征,又能考慮到網(wǎng)絡(luò)的電氣特點(diǎn)。具體表現(xiàn)在:電氣距離較小的節(jié)點(diǎn)在拓?fù)浣Y(jié)構(gòu)中距離較近,通常在同一片區(qū)域內(nèi)。

    節(jié)點(diǎn)電氣距離矩陣D中元素計(jì)算如下:

    (5)

    式中:Dst為節(jié)點(diǎn)s和節(jié)點(diǎn)t之間的電氣距離;Zss、Ztt和Zst分別為阻抗矩陣Z中對應(yīng)元素。

    按照下述步驟對系統(tǒng)節(jié)點(diǎn)進(jìn)行聚類:

    a) 將第一個(gè)節(jié)點(diǎn)歸為第一類。

    b) 對于除第一個(gè)節(jié)點(diǎn)外的任一節(jié)點(diǎn)s,根據(jù)節(jié)點(diǎn)電氣距離矩陣D得到已有各個(gè)分類中的節(jié)點(diǎn)與節(jié)點(diǎn)s之間的最大電氣距離Ds.rmax,其中r=1,…,u為已存在的類別編號,u為已存在的類別數(shù)目。若對于任何一個(gè)已有類別r,均有Ds.rmax>ε2(ε2為預(yù)先設(shè)定的電氣距離閾值),則節(jié)點(diǎn)s不屬于任何一個(gè)已有類,將其自身作為一個(gè)新的類;否則,節(jié)點(diǎn)s屬于一個(gè)已有的類別。

    c) 若存在一個(gè)已有類別r滿足Ds.rmax<ε1,則節(jié)點(diǎn)s屬于類別r;若有兩個(gè)或兩個(gè)以上已有類別滿足該條件,則分別計(jì)算節(jié)點(diǎn)s與滿足條件的各個(gè)類別的平均電氣距離Ds.ravg(即類別內(nèi)各個(gè)節(jié)點(diǎn)與節(jié)點(diǎn)s之間電氣距離的平均值),取其最小值所對應(yīng)的類別作為節(jié)點(diǎn)s的所屬類別。

    d) 遍歷系統(tǒng)所有節(jié)點(diǎn),同一類內(nèi)節(jié)點(diǎn)為一個(gè)分區(qū),不同類內(nèi)節(jié)點(diǎn)為不同分區(qū),形成電網(wǎng)初步分區(qū)。

    2.2 電網(wǎng)最終分區(qū)

    聚類方法能從整體上反應(yīng)所研究問題的關(guān)鍵特點(diǎn),但通常存在誤差,如系統(tǒng)中部分長程線路或關(guān)鍵線路本應(yīng)是重點(diǎn)分析對象,若將其歸為某一分區(qū)內(nèi),則會被忽略。本步驟旨在通過劃分子區(qū)的方法辨識區(qū)域內(nèi)長程線路或關(guān)鍵線路,并加以分析。

    首先形成長程線路集Llong,Llong中的線路應(yīng)滿足xk>ε3(ε3為預(yù)先設(shè)定的線路長度閾值);關(guān)鍵線路集Lkey按第一節(jié)所述方法形成。

    對任一分區(qū),若其中包含長程線路集Llong或關(guān)鍵線路集Lkey中的線路,則將線路兩端節(jié)點(diǎn)分別劃到兩個(gè)子類中,該分區(qū)中其他節(jié)點(diǎn)按照與該兩端節(jié)點(diǎn)之間電氣距離的大小劃分到兩個(gè)子類中;若子類中仍包含長程線路集或關(guān)鍵線路集中的線路,則繼續(xù)進(jìn)行子類劃分,直到每個(gè)類別中均不含有長程線路或關(guān)鍵線路。

    考慮線路潮流方向,形成電網(wǎng)有向分區(qū)拓?fù)鋱D。由于輸電斷面有潮流方向一致性的約束,可能存在兩個(gè)或以上分區(qū)組成的分區(qū)組,其中任何兩個(gè)分區(qū)由于潮流方向的限值均不可能被分在輸電斷面的兩側(cè),則可將這樣的分區(qū)組中的幾個(gè)分區(qū)進(jìn)行合并,以減少后續(xù)工作中不必要的工作量。

    至此,電網(wǎng)最終分區(qū)形成,后續(xù)只需在分區(qū)間聯(lián)絡(luò)線中選擇滿足輸電斷面特征的線路集合構(gòu)成輸電斷面,分析效率大大提高。

    3 基于關(guān)鍵支路和電網(wǎng)分區(qū)的關(guān)鍵斷面分析

    前文闡述了關(guān)鍵線路辨識及電網(wǎng)分區(qū)方法,基于此,選擇系統(tǒng)中最關(guān)鍵的若干條輸電線路作為主導(dǎo)線路,并在電網(wǎng)最終分區(qū)的區(qū)域間聯(lián)絡(luò)線中選擇主導(dǎo)線路所在的輸電斷面為關(guān)鍵輸電斷面。本文中輸電斷面所包含的線路集合應(yīng)滿足3個(gè)條件:輸電斷面包含的線路集合構(gòu)成電網(wǎng)的割集;輸電斷面包含的各條線路有功潮流基本一致;主導(dǎo)線路開斷對斷面內(nèi)其他線路影響較大。

    以下為基于關(guān)鍵支路辨識和電網(wǎng)分區(qū)的關(guān)鍵輸電斷面選取步驟。

    表1 不同情況下IEEE 39節(jié)點(diǎn)系統(tǒng)關(guān)鍵線路辨識結(jié)果(本文方法)

    關(guān)鍵程度排序情況1支路Bc情況2支路Bc情況3支路Bc情況4支路Bc116-171664619-201738015-162387429-3896922216-191307116-191717216-172386516-1764311315-161210516-171314514-152255726-2760555414-151108715-16101923-42168117-2759295517-271059914-159687.94-52065619-205860162-25102823-48040.216-191907116-1956060726-27987317-187477.46-111584728-2947180819-2096863-187415.417-181555326-2846412917-1896012-257385.13-181550826-2946301103-4945817-277326.35-8154572-2539365113-18932026-277030.610-111481025-2638182122-391544-56848.413-141477715-1631101

    注:“a-b”表示兩端節(jié)點(diǎn)編號為a、b的支路,下同。

    a) 選取主導(dǎo)線路。計(jì)算所有線路的綜合介數(shù),選擇綜合介數(shù)最高的若干條關(guān)鍵線路作為主導(dǎo)線路。

    b) 形成系統(tǒng)分區(qū)。按照基于節(jié)點(diǎn)聚類算法的電網(wǎng)分區(qū)方法形成電網(wǎng)最終分區(qū),長程線路及關(guān)鍵線路均包含在最終分區(qū)的區(qū)域間聯(lián)絡(luò)線路中。

    c) 確定受主導(dǎo)線路開斷影響較大的線路集合。對于每一條主導(dǎo)線路l,計(jì)算其與其他任一區(qū)域間聯(lián)絡(luò)線h的支路開斷分布因子

    (6)

    其中

    (7)

    式中:Mh和Ml分別為聯(lián)絡(luò)線h和支路l的節(jié)點(diǎn)-支路關(guān)聯(lián)矢量;ηl=XMl;xh和xl分別為聯(lián)絡(luò)線h和支路l的電抗。

    若Dh-l>ε4(ε4為預(yù)先設(shè)定的支路開斷分布因子閾值),則將聯(lián)絡(luò)線h列入支路l的強(qiáng)相關(guān)支路集,即聯(lián)絡(luò)線h與支路l之間的開斷靈敏度較大。

    d) 根據(jù)潮流方向基本一致性確定輸電斷面。在上一步選出的線路集合中搜索割集斷面,并選擇支路潮流方向基本一致的斷面作為輸電斷面,所謂基本一致,就是在一定程度上允許個(gè)別支路潮流反向,即只要滿足Pl.h<ε5Pl.max(Pl.h為反向支路的有功潮流,Pl.max為該斷面各支路有功潮流中的最大值,ε5為預(yù)先設(shè)定的反向潮流閾值系數(shù)),則仍然認(rèn)為其構(gòu)成輸電斷面。

    e) 遍歷所有主導(dǎo)線路完成關(guān)鍵輸電斷面的選取。

    以上方法共涉及5個(gè)閾值,即關(guān)鍵線路綜合介數(shù)閾值ε1、電氣距離閾值ε2、線路長度閾值ε3、支路開斷分布因子閾值ε4和反向潮流閾值系數(shù)ε5,這5個(gè)閾值可根據(jù)實(shí)際情況和具體分析需要設(shè)置,用以對關(guān)鍵線路、系統(tǒng)分區(qū)、長程連接線及關(guān)鍵斷面的數(shù)目進(jìn)行控制。

    4 算例分析

    為驗(yàn)證本文方法的有效性,在IEEE 39節(jié)點(diǎn)測試系統(tǒng)上進(jìn)行仿真,該系統(tǒng)包含10個(gè)發(fā)電機(jī)節(jié)點(diǎn)和46條輸電線路,如圖1所示。

    圖1 IEEE 39節(jié)點(diǎn)系統(tǒng)

    4.1 考慮不確定因素的電力系統(tǒng)關(guān)鍵線路辨識

    本算例考慮以下4種情況:

    a) 情況1,不考慮負(fù)荷和風(fēng)電的功率波動;

    b) 情況2,各個(gè)負(fù)荷節(jié)點(diǎn)功率的誤差取5%;

    c) 情況3,部分負(fù)荷節(jié)點(diǎn)功率的誤差取8%,另一部分負(fù)荷節(jié)點(diǎn)功率的誤差取3%;

    d) 情況4,各個(gè)負(fù)荷節(jié)點(diǎn)功率的誤差取5%,同時(shí)在節(jié)點(diǎn)38添加300 MW的風(fēng)電功率,其誤差取10%,該節(jié)點(diǎn)的常規(guī)機(jī)組出力相應(yīng)減少300 MW。

    4種情況下關(guān)鍵線路的辨識結(jié)果見表1。表2列舉了幾種現(xiàn)有方法的計(jì)算結(jié)果,這幾種方法均未計(jì)及不確定因素的影響,故將其作為本例中情況1辨識結(jié)果的對比。

    表2 不同方法下IEEE 39節(jié)點(diǎn)系統(tǒng)關(guān)鍵線路辨識結(jié)果

    從表1情況1的結(jié)果和表2可以看出,在忽略不確定因素影響時(shí),本文方法與文獻(xiàn)[11,15-16]方法所得結(jié)果具有一定的相似性,其中支路15-16、16-17、16-19、14-15在3種對比方法中同時(shí)出現(xiàn),說明其在系統(tǒng)中具有關(guān)鍵地位,本文方法不僅得出了同樣的結(jié)論,而且在本文的排序中,這4條支路位列前4。從網(wǎng)絡(luò)結(jié)構(gòu)來看,支路26-27是發(fā)電機(jī)38功率外送的關(guān)鍵通道,也是發(fā)電機(jī)37功率外送的通道之一;從運(yùn)行方式來看,發(fā)電機(jī)38的凈輸出功率是所有發(fā)電機(jī)中最大的,而其所在分區(qū)對系統(tǒng)主干網(wǎng)絡(luò)的供電功率幾乎全部通過支路26-27傳送,綜合考慮網(wǎng)絡(luò)結(jié)構(gòu)及運(yùn)行方式,支路26-27應(yīng)當(dāng)成為關(guān)鍵支路,而該支路在表2所示3種對比方法中均被忽略??梢?,本文所選關(guān)鍵線路在系統(tǒng)發(fā)電-負(fù)荷節(jié)點(diǎn)對之間的功率傳輸過程中,要么在網(wǎng)絡(luò)結(jié)構(gòu)中占據(jù)關(guān)鍵位置,要么在當(dāng)前運(yùn)行方式下承擔(dān)了較重的傳輸任務(wù),因而其重要程度較高。

    從表1可看出,負(fù)荷節(jié)點(diǎn)和發(fā)電節(jié)點(diǎn)的功率波動會對關(guān)鍵支路的辨識結(jié)果產(chǎn)生影響。如支路 19-20 在情況2中關(guān)鍵程度排第1位,而在情況1中排在第8位,在情況3中沒有排進(jìn)前12位。從蒙特卡洛模擬得到的支路功率方差情況來看,情況2中該支路的方差排在第3位,表明該情況下其受負(fù)荷波動的影響較大,因此與情況1相比排名靠前,而情況3中該支路的方差排在第15位,表明該情況下其受負(fù)荷波動的影響不大,因此該情況下其關(guān)鍵程度未排在前面。支路16-19與支路19-20類似。支路29-38在情況1到情況3中關(guān)鍵程度排名均靠后,而在情況4中排在第1位,原因是在情況4中節(jié)點(diǎn)38接入了300 MW風(fēng)電,其波動給該支路帶來很大影響,從蒙特卡洛模擬結(jié)果來看,在情況2和情況3負(fù)荷波動情況下,該支路功率基本未受影響,兩種情況下其功率波動方差均基本為零,而在情況4中其方差躍居第6位,說明其受風(fēng)電功率波動影響很大。

    4.2 電力系統(tǒng)分區(qū)

    本算例中,取電氣距離閾值ε1=0.018,關(guān)鍵線路綜合介數(shù)閾值ε2=9 000,線路長度閾值ε3=0.017 7。電網(wǎng)分區(qū)結(jié)果如圖2所示,關(guān)鍵線路集Lkey見表1情況1,長程線路集Llong見表3。

    圖2 IEEE 39節(jié)點(diǎn)系統(tǒng)分區(qū)圖

    表3 IEEE 39節(jié)點(diǎn)系統(tǒng)長程線路集

    支路電抗標(biāo)幺值支路電抗標(biāo)幺值支路電抗標(biāo)幺值1-20.04119-390.025025-260.03231-390.025014-150.021726-280.04743-40.021316-190.019526-290.06258-90.036323-240.0350

    本算例將系統(tǒng)分為15個(gè)區(qū),區(qū)域間聯(lián)絡(luò)線路包含了除線路19-20和線路23-24以外所有關(guān)鍵線路和長程線路(其中19-20為變壓器支路,而23-24由于和16-21、21-22及22-23的潮流方向不一致而無法構(gòu)成輸電斷面)。輸電斷面從區(qū)域間聯(lián)絡(luò)線中選取,既不會遺漏關(guān)鍵線路和長程線路,又有效降低了從全網(wǎng)中選取輸電斷面的工作量。本文所提聚類方法及其修正措施有效避免了關(guān)鍵線路的遺漏,同時(shí),電氣距離閾值ε1、關(guān)鍵線路綜合介數(shù)閾值ε2和線路長度閾值ε3可根據(jù)實(shí)際情況和具體分析需要靈活設(shè)置,具有較強(qiáng)的適應(yīng)性。

    4.3 考慮不確定因素的電力系統(tǒng)關(guān)鍵輸電斷面選取

    前面兩個(gè)算例已驗(yàn)證了本文關(guān)鍵線路辨識和系統(tǒng)分區(qū)方法的合理性,在此基礎(chǔ)上,本算例驗(yàn)證關(guān)鍵輸電斷面選取的合理性。支路開斷分布因子閾值ε4和反向潮流閾值系數(shù)ε5可以根據(jù)實(shí)際情況和具體分析需要進(jìn)行設(shè)置,本算例中ε4=0.2、ε5=0.18。4種情況下的關(guān)鍵輸電斷面選取結(jié)果見表4,其中粗體代表輸電斷面的主導(dǎo)線路。

    表4 不同情況下IEEE 39節(jié)點(diǎn)系統(tǒng)關(guān)鍵輸電斷面選取結(jié)果

    從表4可以看出,負(fù)荷節(jié)點(diǎn)和發(fā)電節(jié)點(diǎn)的功率波動會對關(guān)鍵輸電斷面的選取結(jié)果產(chǎn)生影響。通過蒙特卡洛模擬得到各斷面功率在不同情況下的方差及排序,用以說明不同波動下系統(tǒng)斷面受不確定因素的影響程度。在所有關(guān)鍵輸電斷面中,情況2、情況3中斷面11在該情況下方差排在第一位,表明其功率不確定程度很高,受負(fù)荷波動影響很大,因此該情況下將其列入關(guān)鍵斷面;情況4中,斷面11和斷面12在該情況下方差顯著增大,表明其受風(fēng)電接入的影響很嚴(yán)重,因此該情況下將其列入關(guān)鍵斷面。而以往方法中4種情況選取的關(guān)鍵斷面結(jié)果相同,無法體現(xiàn)出系統(tǒng)中不確定因素帶來的影響。

    由表4還可看出,16-17和16-19主導(dǎo)的斷面以及斷面{2-25,15-16,17-18}和{2-25,15-16,3-18}在4種情況下均被選為關(guān)鍵輸電斷面,表明由網(wǎng)架結(jié)構(gòu)及運(yùn)行方式?jīng)Q定了其在不同運(yùn)行環(huán)境下均對系統(tǒng)安全運(yùn)行起到重要作用。

    當(dāng)一個(gè)斷面功率的不確定性很大時(shí),其安全程度的可信性將降低,其對系統(tǒng)安全帶來的潛在危害不容忽視。算例說明本文方法所選斷面不僅體現(xiàn)了其在常規(guī)潮流分析中的重要性,還反映了其受不確定因素影響的程度。

    5 結(jié)論

    本文計(jì)及負(fù)荷與間歇式電源波動在網(wǎng)絡(luò)中的傳播特性,提出考慮不確定因素的關(guān)鍵輸電斷面識別方法,該方法具有以下特點(diǎn):

    a) 綜合介數(shù)指標(biāo)綜合考慮了網(wǎng)絡(luò)結(jié)構(gòu)、運(yùn)行方式及環(huán)境中的不確定因素,將其應(yīng)用于引導(dǎo)關(guān)鍵輸電斷面的選擇,能夠體現(xiàn)出不同運(yùn)行環(huán)境下系統(tǒng)關(guān)鍵斷面的差異。根據(jù)不同時(shí)段的負(fù)荷或間歇式能源波動特性以及系統(tǒng)運(yùn)行方式的變化搜索系統(tǒng)關(guān)鍵輸電斷面,調(diào)整電網(wǎng)運(yùn)行中的關(guān)注點(diǎn),為有針對性地采取系統(tǒng)安全控制措施提供依據(jù)。

    b) 所提電網(wǎng)分區(qū)算法,既保證了分區(qū)內(nèi)節(jié)點(diǎn)之間相互聯(lián)系緊密,又避免了將關(guān)鍵線路和長程線路劃分區(qū)內(nèi)而在選取輸電斷面時(shí)將其忽略。通過分區(qū)對電網(wǎng)進(jìn)行簡化,使搜索輸電斷面的工作量大為減少。

    c) 關(guān)鍵線路綜合介數(shù)閾值、電氣距離閾值、線路長度閾值、支路開斷分布因子閾值和反向潮流閾值系數(shù)可以根據(jù)實(shí)際情況和具體分析需要進(jìn)行設(shè)置,增加了本文方法的靈活性和適用性。

    [1] MEI Shengwei,ZHANG Xuemin,CAO Ming.Power Grid Complexity[M].Beijing:Tsinghua University Press,2011.

    [2] LUO G,SHI D,CHEN J,et al.Automatic Identification of Transmission Sections Based on Complex Network Theory[J].Iet Generation Transmission & Distribution,2014,8(8):1203-1210.

    [3] 甘德樹,歐思源,彭顯剛,等.考慮風(fēng)險(xiǎn)和收益平衡的輸電斷面供電方案優(yōu)化[J].廣東電力,2014(3):63-67.

    GAN Deshu,OU Siyuan,PENG Xiangang,et al.Optimization for Power Transmission Section Power Supply Scheme Considering Risk and Profit Balance[J].Guangdong Electric Power,2014(3):63-67.

    [4] 張文朝,何玉龍,顧雪平,等.單輸電通道中輸電斷面靜穩(wěn)極限的快速估算[J].電網(wǎng)技術(shù),2012,36(5):92-95.

    ZHANG Wenchao,HE Yulong,GU Xueping,et al.A Fast Method to Estimate Steady-state Stability Limit of Tie Line in Single Power Transmission Channel[J].Power System Technology,2012,36(5):92-95.

    [5] 崔曉丹,方勇杰,周霞,等.基于風(fēng)險(xiǎn)的輸電斷面經(jīng)濟(jì)輸電功率計(jì)算[J].電力系統(tǒng)自動化,2014,38(8):30-34.

    CUI Xiaodan,F(xiàn)ANG Yongjie,ZHOU Xia,et al.Risk-based Economical Transmission Power Calculation of Transmission Sections[J].Automation of Electric Power Systems,2014,38(8):30-34.

    [6] 趙峰,孫宏斌,張伯明.基于電氣分區(qū)的輸電斷面及其自動發(fā)現(xiàn)[J].電力系統(tǒng)自動化,2011,35(5):42-47.

    ZHAO Feng,SUN Hongbin,ZHANG Boming.Electrical Zone Division Based Automatic Discovery of Flowgates[J].Automation of Electric Power Systems,2011,35(5):42-47.

    [7] 羅鋼,陳金富,石東源,等.基于復(fù)雜網(wǎng)絡(luò)理論的關(guān)鍵輸電斷面分析[J].中國電機(jī)工程學(xué)報(bào),2013,33(25):147-155.

    LUO Gang,CHEN Jinfu,SHI Dongyuan,et al.Key Transmission Sections Analysis Based on Complex Network Theory[J].Proceedings of the CSEE,2013,33(25):147-155.

    [8] 周德才,張保會,姚峰,等.基于圖論的輸電斷面快速搜索[J].中國電機(jī)工程學(xué)報(bào),2006,26(12):32-38.

    ZHOU Decai,ZHANG Baohui,YAO Feng,et al.Fast Search for Transmission Section Based on Graph Theory[J].Proceedings of the CSEE,2006,26(12):32-38.

    [9] 苗世洪,馬帥,尚亞男,等.基于割點(diǎn)和路徑搜索的輸電斷面快速識別方法[J].電力系統(tǒng)自動化,2014,38(2):39-45.

    MIAO Shihong,MA Shuai,SHANG Yanan,et al.A Fast Recognition Method of Transmission Section Based on Cut-vertex and Path Search[J].Automation of Electric Power Systems,2014,38(2):39-45.

    [10] KINNEY R,CRUCITTI P,ALBERT R,et al.Modeling Cascading Failures in the North American Power Grid[J].The European Physical Journal B:Condensed Matter Physics,2005,46(1):101-107.

    [11] 曹一家,陳曉鋼,孫可.基于復(fù)雜網(wǎng)絡(luò)理論的大型電力系統(tǒng)脆弱線路識別[J].電力自動化設(shè)備,2006,26(12):1-5,31.

    CAO Yijia,CHEN Xiaogang,SUN Ke.Identification of Vulnerable Lines in Power Grid Based on Complex Network Theory[J].Electric Power Automation Equipment,2006,26(12):1-5,31.

    [12] 程臨燕,張保會,李光輝,等.采用有向電氣介數(shù)的脆弱線路選取[J].西安交通大學(xué)學(xué)報(bào),2011,45(6):91-96.

    CHENG Linyan,ZHANG Baohui,LI Guanghui,et al.Search for Vulnerable Line Based on Directed Electric Betweenness[J].Journal of Xi’an Jiaotong University,2011,45(6):91-96.

    [13] 劉文穎,梁才,徐鵬,等.基于潮流介數(shù)的電力系統(tǒng)關(guān)鍵線路辨識[J].中國電機(jī)工程學(xué)報(bào),2013,33(31):90-99.

    LIU Wenying,LIANG Cai,XU Peng,et al.Identification of Critical Line in Power Systems Based on Flow Betweenness[J].Proceedings of the CSEE,2013,33(31):90-99.

    [14] 張伯明,陳壽孫,嚴(yán)正.高等電力網(wǎng)絡(luò)分析[M].北京:清華大學(xué)出版社,2011:212-214.

    [15] 鞠文云,李銀紅.基于最大流傳輸貢獻(xiàn)度的電力網(wǎng)關(guān)鍵線路和節(jié)點(diǎn)辨識[J].電力系統(tǒng)自動化,2012,36(9):6-12.

    JU Wenyun,LI Yinhong.Identification of Critical Lines and Nodes in Power Grid Based on Maximum Flow Transmission Contribution Degree[J].Automation of Electric Power Systems,2012,36(9):6-12.

    [16] 丁明,過羿,張晶晶.基于效用風(fēng)險(xiǎn)熵的復(fù)雜電網(wǎng)連鎖故障脆弱性辨識[J].電力系統(tǒng)自動化,2013,37(17):52-57.

    DING Ming,GUO Yi,ZHANG Jingjing.Vulnerability Identification for Cascading Failures of Complex Power Grid Based on Effect Risk Entropy[J].Automation of Electric Power Systems,2013,37(17):52-57.

    (編輯 彭艷)

    Analysis on Key Power Transmission Section Considering Uncertain Factors

    JIN Bingjie1, LIANG Chen2

    (1.Guangdong Power Grid Development Research Institute Co., Ltd., Guangzhou, Guangdong 510080, China; 2.Maintenance & Test Center, UHV Power Transmission Company, CSG, Guangzhou, Guangdong 510663, China)

    Considering propagation characteristics of load and batch power source fluctuation in network, this paper presents a method for key branch identification based on comprehensive betweenness which is used to fully reflect weak links of the system in particular operational environment or with specific network structure and running model.Then it proposes a power grid partitioning method based on node electrical distance clustering and its modification measures so as to avoid missing important lines.Finally, on the basis of grid partitioning and according to features of power transmission section, it selects key branch leading power transmission section as the key section of system.IEEE 39 node system example is used to verify effectiveness of this method.It is able to search for the key power transmission section according to load or batch energy fluctuation characteristic and system running mode at any time in actual application, so as to provide basis for adopting security control measures for the system.

    comprehensive betweenness; key branch; power grid partitioning; power transmission section

    2016-04-17

    10.3969/j.issn.1007-290X.2016.11.023

    TM 73

    A

    1007-290X(2016)11-0120-08

    靳冰潔(1989),女,河北蠡縣人。工學(xué)碩士,主要研究方向?yàn)殡娏ο到y(tǒng)分析、規(guī)劃。

    梁晨(1989),男,廣東高州人。助理工程師,工學(xué)學(xué)士,主要研究方向?yàn)樽儔浩黝愒O(shè)備試驗(yàn)、技術(shù)監(jiān)督和技術(shù)分析。

    猜你喜歡
    介數(shù)支路分區(qū)
    上海實(shí)施“分區(qū)封控”
    基于限流可行方案邊界集的最優(yōu)支路投切
    能源工程(2020年6期)2021-01-26 00:55:22
    浪莎 分區(qū)而治
    多支路兩跳PF協(xié)作系統(tǒng)的誤碼性能
    利用支路參數(shù)的狀態(tài)估計(jì)法辨識拓?fù)溴e誤
    基于電氣介數(shù)的電力系統(tǒng)脆弱線路辨識
    基于SAGA聚類分析的無功電壓控制分區(qū)
    電測與儀表(2015年8期)2015-04-09 11:50:16
    基于多種群遺傳改進(jìn)FCM的無功/電壓控制分區(qū)
    電測與儀表(2015年7期)2015-04-09 11:40:16
    多并聯(lián)支路型可控電抗器短路電抗對支路電抗和電流的影響
    樹形網(wǎng)絡(luò)的平均介數(shù)*
    久久亚洲国产成人精品v| 久久久久性生活片| 国产精品伦人一区二区| 免费人成在线观看视频色| 99视频精品全部免费 在线| 91av网一区二区| 又粗又硬又长又爽又黄的视频 | 日本一二三区视频观看| 国内久久婷婷六月综合欲色啪| 亚洲精华国产精华液的使用体验 | 国产片特级美女逼逼视频| 国产又黄又爽又无遮挡在线| 日本熟妇午夜| 男人舔女人下体高潮全视频| 日本av手机在线免费观看| 99热全是精品| 亚洲性久久影院| 久久精品人妻少妇| 中文字幕熟女人妻在线| 国产亚洲欧美98| 人人妻人人看人人澡| 日韩欧美 国产精品| 久久久久久久久久黄片| 1024手机看黄色片| 天堂av国产一区二区熟女人妻| 欧美最新免费一区二区三区| eeuss影院久久| 天美传媒精品一区二区| 欧美色欧美亚洲另类二区| 一级毛片aaaaaa免费看小| 国产黄片美女视频| 99久久中文字幕三级久久日本| 一个人观看的视频www高清免费观看| 国产一级毛片在线| 看十八女毛片水多多多| 国产精品一及| 1024手机看黄色片| 人人妻人人看人人澡| 大型黄色视频在线免费观看| 狂野欧美激情性xxxx在线观看| 精品人妻熟女av久视频| 国产伦理片在线播放av一区 | 最后的刺客免费高清国语| 精品久久久久久成人av| 99视频精品全部免费 在线| 色综合亚洲欧美另类图片| 偷拍熟女少妇极品色| 插阴视频在线观看视频| 亚洲在线观看片| 亚洲中文字幕一区二区三区有码在线看| 国产精品免费一区二区三区在线| 亚洲欧美清纯卡通| 欧美潮喷喷水| av在线观看视频网站免费| 舔av片在线| 久久精品国产亚洲网站| 亚洲国产日韩欧美精品在线观看| 日韩欧美精品v在线| 一级毛片久久久久久久久女| 女的被弄到高潮叫床怎么办| 亚洲最大成人中文| 久久久午夜欧美精品| 色哟哟·www| 精品人妻熟女av久视频| 日韩成人伦理影院| 亚洲成人久久性| 国产亚洲91精品色在线| 午夜老司机福利剧场| 又爽又黄a免费视频| 在线a可以看的网站| av免费观看日本| 91在线精品国自产拍蜜月| 国产一级毛片在线| 亚洲经典国产精华液单| 国产成人91sexporn| 日韩欧美精品免费久久| 久久久久久久亚洲中文字幕| 三级国产精品欧美在线观看| a级毛片免费高清观看在线播放| 亚洲国产精品成人久久小说 | 国产精品女同一区二区软件| 成人二区视频| 久久久久久久久久成人| 最新中文字幕久久久久| 亚洲四区av| 啦啦啦观看免费观看视频高清| 亚洲中文字幕日韩| 九九热线精品视视频播放| 美女xxoo啪啪120秒动态图| 国产视频内射| 黄色日韩在线| 亚洲av成人av| 少妇人妻精品综合一区二区 | 亚洲精品色激情综合| 久久久色成人| 精品国内亚洲2022精品成人| 日韩精品青青久久久久久| 乱人视频在线观看| 亚洲国产精品国产精品| ponron亚洲| 91久久精品国产一区二区三区| 非洲黑人性xxxx精品又粗又长| 亚洲成av人片在线播放无| 村上凉子中文字幕在线| 一级毛片久久久久久久久女| 免费看av在线观看网站| 不卡视频在线观看欧美| 亚洲av免费在线观看| 好男人在线观看高清免费视频| 欧美潮喷喷水| 少妇人妻一区二区三区视频| 卡戴珊不雅视频在线播放| 欧洲精品卡2卡3卡4卡5卡区| 国产成人精品久久久久久| 九色成人免费人妻av| 日韩 亚洲 欧美在线| 免费av观看视频| 深夜a级毛片| 亚洲欧美日韩卡通动漫| 成人综合一区亚洲| 亚洲欧美精品专区久久| 欧美成人一区二区免费高清观看| 99精品在免费线老司机午夜| 国产精品无大码| 日韩av在线大香蕉| 99热这里只有是精品50| 国产爱豆传媒在线观看| 99国产精品一区二区蜜桃av| 国产日韩欧美在线精品| 国产精品无大码| 夜夜爽天天搞| 国产日韩欧美在线精品| 日本色播在线视频| 男人和女人高潮做爰伦理| 看片在线看免费视频| 免费看a级黄色片| 一区二区三区四区激情视频 | 狂野欧美激情性xxxx在线观看| 性欧美人与动物交配| 网址你懂的国产日韩在线| 久久久久久久久久黄片| 亚洲精品456在线播放app| 一个人看的www免费观看视频| 丰满的人妻完整版| 毛片女人毛片| 国产免费男女视频| 国产高清视频在线观看网站| 又爽又黄无遮挡网站| 看十八女毛片水多多多| 亚洲精品国产成人久久av| 欧美成人免费av一区二区三区| av福利片在线观看| 激情 狠狠 欧美| 老司机福利观看| 日本-黄色视频高清免费观看| 男女边吃奶边做爰视频| 国产精品久久电影中文字幕| 一进一出抽搐gif免费好疼| 我要看日韩黄色一级片| 国产大屁股一区二区在线视频| 在线观看免费视频日本深夜| 最近视频中文字幕2019在线8| 特大巨黑吊av在线直播| 日韩欧美 国产精品| 一个人免费在线观看电影| 国内久久婷婷六月综合欲色啪| 边亲边吃奶的免费视频| 夜夜爽天天搞| 尾随美女入室| 亚洲精品色激情综合| 久久久久性生活片| 亚洲中文字幕一区二区三区有码在线看| 我要看日韩黄色一级片| 老女人水多毛片| 午夜福利高清视频| 亚洲五月天丁香| 欧美一级a爱片免费观看看| 亚洲精品亚洲一区二区| 午夜福利在线观看吧| 可以在线观看毛片的网站| 亚洲自拍偷在线| 99riav亚洲国产免费| 午夜视频国产福利| 国产精品一及| 国产日本99.免费观看| 一个人看视频在线观看www免费| 国产精品,欧美在线| 久久精品夜色国产| 精品熟女少妇av免费看| 日韩av在线大香蕉| 中国美白少妇内射xxxbb| 亚洲综合色惰| 亚洲欧美中文字幕日韩二区| 看十八女毛片水多多多| 偷拍熟女少妇极品色| 乱码一卡2卡4卡精品| 亚洲精品乱码久久久久久按摩| 欧美zozozo另类| 黄色欧美视频在线观看| 国产蜜桃级精品一区二区三区| 亚洲成人久久爱视频| 国产精品一区二区三区四区免费观看| 国产一区二区激情短视频| 亚洲av电影不卡..在线观看| 日韩av不卡免费在线播放| 国产成人a∨麻豆精品| 看片在线看免费视频| 欧美日韩一区二区视频在线观看视频在线 | 日产精品乱码卡一卡2卡三| 卡戴珊不雅视频在线播放| 国产欧美日韩精品一区二区| 色综合亚洲欧美另类图片| 非洲黑人性xxxx精品又粗又长| av天堂中文字幕网| 欧美性猛交╳xxx乱大交人| 欧美成人免费av一区二区三区| 人妻夜夜爽99麻豆av| а√天堂www在线а√下载| 国产色爽女视频免费观看| 国产av在哪里看| 久久午夜福利片| 看黄色毛片网站| 国产真实伦视频高清在线观看| 一本一本综合久久| 三级经典国产精品| 免费观看人在逋| 日韩视频在线欧美| 日日撸夜夜添| 免费黄网站久久成人精品| ponron亚洲| 好男人视频免费观看在线| kizo精华| 国产成人福利小说| 真实男女啪啪啪动态图| 少妇被粗大猛烈的视频| 午夜激情欧美在线| 99热全是精品| 不卡一级毛片| 天天躁日日操中文字幕| 国产伦理片在线播放av一区 | 日本三级黄在线观看| 91午夜精品亚洲一区二区三区| 18禁在线无遮挡免费观看视频| 看片在线看免费视频| 91麻豆精品激情在线观看国产| 免费看av在线观看网站| 国产精品综合久久久久久久免费| 18禁裸乳无遮挡免费网站照片| 欧美+日韩+精品| 国产男人的电影天堂91| 亚洲美女搞黄在线观看| 午夜福利在线在线| 午夜爱爱视频在线播放| 三级国产精品欧美在线观看| 超碰av人人做人人爽久久| 亚洲三级黄色毛片| 99久久中文字幕三级久久日本| 久久精品国产99精品国产亚洲性色| 不卡一级毛片| 亚洲欧美清纯卡通| 一级黄色大片毛片| 亚洲五月天丁香| 中文资源天堂在线| 亚洲国产欧美人成| 亚洲av成人精品一区久久| avwww免费| 午夜激情福利司机影院| 亚洲中文字幕日韩| 一边摸一边抽搐一进一小说| 三级经典国产精品| 丝袜美腿在线中文| 日本熟妇午夜| 一个人免费在线观看电影| 国内精品久久久久精免费| 国内揄拍国产精品人妻在线| 国产v大片淫在线免费观看| 精品人妻偷拍中文字幕| 国产精华一区二区三区| 九草在线视频观看| 日韩成人av中文字幕在线观看| 一级毛片aaaaaa免费看小| 久久久久久伊人网av| 女的被弄到高潮叫床怎么办| 精品人妻视频免费看| 成人av在线播放网站| 亚洲经典国产精华液单| 一区二区三区四区激情视频 | 少妇猛男粗大的猛烈进出视频 | 国模一区二区三区四区视频| 亚洲经典国产精华液单| 免费搜索国产男女视频| 亚洲国产精品合色在线| av在线播放精品| av专区在线播放| 99久久中文字幕三级久久日本| 欧美成人精品欧美一级黄| 亚洲精品乱码久久久久久按摩| 久久九九热精品免费| 欧美不卡视频在线免费观看| 日日摸夜夜添夜夜爱| 国产午夜精品论理片| 国国产精品蜜臀av免费| 成人欧美大片| 色哟哟·www| 亚洲在线自拍视频| 亚洲av免费在线观看| 日本在线视频免费播放| 国内揄拍国产精品人妻在线| 亚洲欧美精品综合久久99| 亚洲激情五月婷婷啪啪| 午夜激情福利司机影院| 亚洲精品亚洲一区二区| 一区二区三区免费毛片| 国产精品1区2区在线观看.| 嫩草影院新地址| kizo精华| 听说在线观看完整版免费高清| 最近最新中文字幕大全电影3| 日本色播在线视频| 特大巨黑吊av在线直播| 欧美bdsm另类| 成人亚洲精品av一区二区| 国产一区二区三区av在线 | 91aial.com中文字幕在线观看| 国产麻豆成人av免费视频| 两个人的视频大全免费| 国产人妻一区二区三区在| 成年女人看的毛片在线观看| 日本欧美国产在线视频| 亚洲欧美精品自产自拍| 午夜精品在线福利| 身体一侧抽搐| 男人和女人高潮做爰伦理| 中文亚洲av片在线观看爽| 免费搜索国产男女视频| 国产精品一区二区在线观看99 | 亚洲欧美成人综合另类久久久 | 少妇高潮的动态图| 69av精品久久久久久| 国产精品99久久久久久久久| 精品免费久久久久久久清纯| 免费看光身美女| 深爱激情五月婷婷| 91狼人影院| 久久久久久大精品| 在线国产一区二区在线| 久久久精品欧美日韩精品| 成人午夜精彩视频在线观看| 国产精品日韩av在线免费观看| 青春草视频在线免费观看| 亚洲国产欧洲综合997久久,| 啦啦啦韩国在线观看视频| 精品一区二区免费观看| 尤物成人国产欧美一区二区三区| 国产精品,欧美在线| 天堂影院成人在线观看| 免费不卡的大黄色大毛片视频在线观看 | 一夜夜www| 色5月婷婷丁香| 中国国产av一级| 亚洲欧美成人精品一区二区| 深爱激情五月婷婷| 久久鲁丝午夜福利片| 综合色丁香网| 久久久久国产网址| 国产成人a∨麻豆精品| 国产白丝娇喘喷水9色精品| av视频在线观看入口| 成人高潮视频无遮挡免费网站| 亚洲av一区综合| 国产精品99久久久久久久久| 女人十人毛片免费观看3o分钟| 免费看光身美女| 日日干狠狠操夜夜爽| 黄色配什么色好看| 亚洲婷婷狠狠爱综合网| 国产熟女欧美一区二区| 99久国产av精品国产电影| 国产一级毛片七仙女欲春2| 亚洲最大成人手机在线| 国产免费一级a男人的天堂| 老司机影院成人| 国内精品美女久久久久久| 亚洲av免费高清在线观看| 三级国产精品欧美在线观看| 午夜爱爱视频在线播放| 国产av在哪里看| 日本黄色片子视频| 亚洲精品成人久久久久久| 国产精品久久久久久久久免| 黄片wwwwww| 久久精品国产99精品国产亚洲性色| 亚洲欧美精品专区久久| 精品99又大又爽又粗少妇毛片| 少妇丰满av| 亚洲精品456在线播放app| www日本黄色视频网| 成人鲁丝片一二三区免费| 少妇人妻一区二区三区视频| 亚洲美女搞黄在线观看| 免费不卡的大黄色大毛片视频在线观看 | 插阴视频在线观看视频| 国产免费一级a男人的天堂| 精品人妻熟女av久视频| 国产精品福利在线免费观看| 校园春色视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美成人精品一区二区| 搡老妇女老女人老熟妇| h日本视频在线播放| videossex国产| 欧美三级亚洲精品| 最近的中文字幕免费完整| 久久久午夜欧美精品| 欧美在线一区亚洲| 小说图片视频综合网站| 麻豆成人午夜福利视频| 日本五十路高清| 美女高潮的动态| eeuss影院久久| 99久久人妻综合| 久久精品国产清高在天天线| 国产日韩欧美在线精品| 欧美性猛交黑人性爽| 麻豆成人av视频| 日日干狠狠操夜夜爽| 黄色日韩在线| 99在线视频只有这里精品首页| 日本免费一区二区三区高清不卡| 免费在线观看成人毛片| 亚洲国产精品sss在线观看| 亚洲欧美精品综合久久99| 国产伦精品一区二区三区视频9| 久久人妻av系列| АⅤ资源中文在线天堂| 亚洲无线观看免费| 亚洲aⅴ乱码一区二区在线播放| 波野结衣二区三区在线| 久久久久九九精品影院| 日韩av不卡免费在线播放| 亚洲一区高清亚洲精品| 久久99热这里只有精品18| 观看免费一级毛片| 97在线视频观看| 国产亚洲精品av在线| 精品免费久久久久久久清纯| 女同久久另类99精品国产91| 亚洲国产精品成人久久小说 | 久久久久久久久久久丰满| 亚洲一级一片aⅴ在线观看| 啦啦啦观看免费观看视频高清| 日韩强制内射视频| 亚洲国产精品合色在线| 国产亚洲精品久久久com| 日韩人妻高清精品专区| 嘟嘟电影网在线观看| 18禁黄网站禁片免费观看直播| 久久精品夜色国产| 日韩精品青青久久久久久| 中出人妻视频一区二区| 3wmmmm亚洲av在线观看| 欧美成人精品欧美一级黄| 91aial.com中文字幕在线观看| 久久精品国产亚洲av天美| 亚州av有码| 国产午夜精品一二区理论片| 国产黄片视频在线免费观看| 午夜免费男女啪啪视频观看| 长腿黑丝高跟| 91狼人影院| 亚洲五月天丁香| 悠悠久久av| 欧美3d第一页| 亚洲av.av天堂| 国产激情偷乱视频一区二区| 国产精品精品国产色婷婷| 亚洲最大成人中文| 国内精品宾馆在线| 人人妻人人澡欧美一区二区| 精品99又大又爽又粗少妇毛片| 免费不卡的大黄色大毛片视频在线观看 | 国产精品美女特级片免费视频播放器| 久久99蜜桃精品久久| 亚洲av不卡在线观看| 日日撸夜夜添| 亚洲精华国产精华液的使用体验 | 99久久久亚洲精品蜜臀av| 日本黄色片子视频| 日日撸夜夜添| 97超碰精品成人国产| 91久久精品国产一区二区成人| 伦理电影大哥的女人| 99久国产av精品国产电影| av在线观看视频网站免费| 国产久久久一区二区三区| 床上黄色一级片| 日韩国内少妇激情av| 男插女下体视频免费在线播放| 亚洲精品亚洲一区二区| 亚洲不卡免费看| 国产av不卡久久| 欧美激情国产日韩精品一区| 插阴视频在线观看视频| 日产精品乱码卡一卡2卡三| 日本三级黄在线观看| 有码 亚洲区| 国内精品一区二区在线观看| 麻豆国产av国片精品| 成人三级黄色视频| 日本一本二区三区精品| 亚洲av一区综合| 亚洲中文字幕日韩| 啦啦啦啦在线视频资源| 国产精品免费一区二区三区在线| 中国美白少妇内射xxxbb| 成人无遮挡网站| 久久欧美精品欧美久久欧美| 国产av一区在线观看免费| 色综合站精品国产| 白带黄色成豆腐渣| 日韩 亚洲 欧美在线| 最新中文字幕久久久久| 亚洲丝袜综合中文字幕| 卡戴珊不雅视频在线播放| 一进一出抽搐gif免费好疼| 成熟少妇高潮喷水视频| 哪里可以看免费的av片| 国产毛片a区久久久久| 亚洲国产精品成人久久小说 | 国内精品久久久久精免费| 一个人免费在线观看电影| 精品熟女少妇av免费看| 婷婷色av中文字幕| 精品午夜福利在线看| 久久久久久久久久成人| 亚洲最大成人手机在线| 免费观看a级毛片全部| 国内精品久久久久精免费| 亚洲av第一区精品v没综合| 不卡一级毛片| 人妻久久中文字幕网| 国产成人福利小说| 婷婷精品国产亚洲av| 少妇人妻精品综合一区二区 | 精品日产1卡2卡| 国产精品人妻久久久影院| 一个人看视频在线观看www免费| 久久午夜亚洲精品久久| www日本黄色视频网| 欧美最黄视频在线播放免费| av视频在线观看入口| 久久久久久九九精品二区国产| 内地一区二区视频在线| 99在线人妻在线中文字幕| 婷婷亚洲欧美| 少妇人妻精品综合一区二区 | 亚洲av电影不卡..在线观看| 亚洲国产欧美人成| 一级黄色大片毛片| 成年免费大片在线观看| 爱豆传媒免费全集在线观看| 国产成人freesex在线| 好男人在线观看高清免费视频| 久久精品国产自在天天线| 美女被艹到高潮喷水动态| 免费人成视频x8x8入口观看| 国产精品久久久久久亚洲av鲁大| 中文在线观看免费www的网站| 免费看美女性在线毛片视频| 身体一侧抽搐| 国产精品久久电影中文字幕| 激情 狠狠 欧美| 天堂中文最新版在线下载 | 自拍偷自拍亚洲精品老妇| 99热这里只有是精品50| 亚洲欧美日韩卡通动漫| 国产三级中文精品| 乱码一卡2卡4卡精品| 搡老妇女老女人老熟妇| 一区二区三区高清视频在线| 99热这里只有精品一区| 国产一区二区三区在线臀色熟女| 中文亚洲av片在线观看爽| 国产一级毛片在线| 久久久久久国产a免费观看| 女同久久另类99精品国产91| 男人狂女人下面高潮的视频| 国产一区二区在线av高清观看| 国产精品1区2区在线观看.| 成年免费大片在线观看| 99在线人妻在线中文字幕| 变态另类成人亚洲欧美熟女| 性欧美人与动物交配| 亚洲精品粉嫩美女一区| 国产老妇女一区| 中出人妻视频一区二区| 久久久精品94久久精品| 日本黄色视频三级网站网址| 蜜桃亚洲精品一区二区三区| 有码 亚洲区| 99视频精品全部免费 在线| 极品教师在线视频| 91久久精品国产一区二区三区| 日本黄色视频三级网站网址| 午夜精品国产一区二区电影 | 精品国内亚洲2022精品成人| 好男人视频免费观看在线| 亚洲自拍偷在线| 国产伦精品一区二区三区四那| 插阴视频在线观看视频| 3wmmmm亚洲av在线观看| 亚洲精华国产精华液的使用体验 | 亚洲欧美日韩无卡精品| 天堂网av新在线| 夜夜爽天天搞|