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

    主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)的時(shí)延穩(wěn)定性分析

    2022-11-22 06:39:38趙聯(lián)港廖小兵
    電力系統(tǒng)保護(hù)與控制 2022年21期
    關(guān)鍵詞:裕度時(shí)延分布式

    樂(lè) 健,綦 淦,趙聯(lián)港,廖小兵

    主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)的時(shí)延穩(wěn)定性分析

    樂(lè) 健1,綦 淦1,趙聯(lián)港1,廖小兵2

    (1.武漢大學(xué)電氣與自動(dòng)化學(xué)院,湖北 武漢 430072;2.武漢工程大學(xué),湖北 武漢 430205)

    分布式經(jīng)濟(jì)調(diào)度模式具有更好的魯棒性、更快的計(jì)算速度和更少的通信量等優(yōu)點(diǎn),適用于高滲透率分布式電源接入的主動(dòng)配電網(wǎng)。但分布式模式下各單元之間通信的時(shí)延將可能破壞系統(tǒng)的穩(wěn)定性。因此對(duì)采用分布式經(jīng)濟(jì)調(diào)度策略的主動(dòng)配電網(wǎng)的時(shí)延穩(wěn)定性進(jìn)行了研究。首先,引入虛擬分布式電源,設(shè)計(jì)了主動(dòng)配電網(wǎng)完全分布式的經(jīng)濟(jì)調(diào)度策略。然后,建立了采用該分布式經(jīng)濟(jì)調(diào)度策略下主動(dòng)配電網(wǎng)對(duì)稱/不對(duì)稱時(shí)延下的穩(wěn)定性分析模型。同時(shí)利用Lyapunov穩(wěn)定性定理和自由權(quán)矩陣方法得出具有更低保守性的系統(tǒng)時(shí)延穩(wěn)定判據(jù),并通過(guò)線性矩陣不等式計(jì)算時(shí)延系統(tǒng)的穩(wěn)定裕度。最后,通過(guò)IEEE-14節(jié)點(diǎn)配電網(wǎng)系統(tǒng)中的仿真算例驗(yàn)證了所提分布式經(jīng)濟(jì)調(diào)度控制系統(tǒng)時(shí)延穩(wěn)定性分析和時(shí)延穩(wěn)定裕度計(jì)算方法的正確性和有效性。

    分布式經(jīng)濟(jì)調(diào)度;通信時(shí)延;時(shí)延穩(wěn)定裕度;Lyapunov穩(wěn)定性判據(jù);線性矩陣不等式

    0 引言

    經(jīng)濟(jì)調(diào)度(economic dispatch, ED)是電力系統(tǒng)中最基本的問(wèn)題之一,它本質(zhì)上是一個(gè)優(yōu)化問(wèn)題,即在滿足電力系統(tǒng)各種約束的條件下,決定各個(gè)發(fā)電機(jī)組的出力以滿足總負(fù)荷需求,同時(shí)使得總發(fā)電成本最小。傳統(tǒng)經(jīng)濟(jì)調(diào)度方法是集中式的,即由調(diào)度中心接收各個(gè)發(fā)電單元上傳的信息,進(jìn)行優(yōu)化計(jì)算后下發(fā)出調(diào)度指令,各個(gè)單元接收并執(zhí)行相應(yīng)調(diào)度指令[1-2]。但是隨著數(shù)量眾多、地理位置分散且特性各異的分布式電源(distributed generation, DG)接入主動(dòng)配電網(wǎng)(active distribution network, ADN),集中式經(jīng)濟(jì)調(diào)度中心的通信和計(jì)算壓力巨大,且無(wú)法提供即插即用功能導(dǎo)致靈活性和可擴(kuò)展性較低[3-5]。

    為克服集中式經(jīng)濟(jì)調(diào)度系統(tǒng)的不足,分布式經(jīng)濟(jì)調(diào)度方法得到重視。分布式調(diào)度模式無(wú)需調(diào)度中心,各發(fā)電單元交換一定的信息,自主計(jì)算并做出決策。因此分布式經(jīng)濟(jì)調(diào)度系統(tǒng)具有更好的魯棒性、更快的計(jì)算速度和更少的通信量[6-11]。文獻(xiàn)[12]提出了Leader-Follower分布式經(jīng)濟(jì)調(diào)度方法,以發(fā)電成本微增率為一致性變量,在所有發(fā)電單元成本微增率相等時(shí)得到最優(yōu)解,但該方法需要獲得功率差額等全局信息以確保功率平衡。文獻(xiàn)[13]提出了一種兩級(jí)迭代算法,避免了對(duì)領(lǐng)導(dǎo)者和全局信息的依賴,但其需要設(shè)置發(fā)電單元出力初值總和等于系統(tǒng)總負(fù)荷,因此不適用于負(fù)荷變化的系統(tǒng)。對(duì)于主動(dòng)配電網(wǎng),其網(wǎng)絡(luò)運(yùn)行管理者可方便地實(shí)時(shí)測(cè)量公共耦合點(diǎn)(common coupling point, PCC)處的功率信息,由此解決了Leader-Follower模式下分布式經(jīng)濟(jì)調(diào)度系統(tǒng)全局信息獲取困難的問(wèn)題[14]。

    考慮到分布式經(jīng)濟(jì)調(diào)度系統(tǒng)中各單元需要經(jīng)由局部通信網(wǎng)絡(luò)交換信息,通信時(shí)延將不可避免。文獻(xiàn)[15]通過(guò)將時(shí)延應(yīng)用到分布式經(jīng)濟(jì)調(diào)度算法中,評(píng)估通信網(wǎng)絡(luò)存在時(shí)間延遲時(shí)算法的有效性和性能,結(jié)果表明時(shí)間延遲將影響算法的收斂性,導(dǎo)致算法要么收斂至一個(gè)錯(cuò)誤的值,要么完全不能收斂,但并沒(méi)有進(jìn)行相應(yīng)的理論分析。文獻(xiàn)[16]給出了一種帶通信固定的全分布式算法,并通過(guò)廣義Nyquist準(zhǔn)則得到了最大允許延遲邊界,但計(jì)算復(fù)雜且保守性較高。文獻(xiàn)[17]分析了具有時(shí)變時(shí)延分布式經(jīng)濟(jì)調(diào)度策略的收斂性問(wèn)題,并給出了算法收斂到全局最優(yōu)點(diǎn)的較低保守性的判據(jù)。上述研究均是從算法和數(shù)值仿真的層面分析分布式經(jīng)濟(jì)調(diào)度算法的收斂性和收斂速度等,而尚無(wú)從分布式控制的角度出發(fā)開(kāi)展系統(tǒng)穩(wěn)定性的研究。對(duì)于實(shí)際電力系統(tǒng)而言,其運(yùn)行方式發(fā)生變化或遭受虛假數(shù)據(jù)攻擊之后,僅僅以經(jīng)濟(jì)性最優(yōu)為目標(biāo)的調(diào)度策略能否滿足系統(tǒng)持續(xù)穩(wěn)定運(yùn)行的要求亟待研究[18-19]。

    文獻(xiàn)[20]提出了考慮聯(lián)絡(luò)線功率穩(wěn)定的微電網(wǎng)魯棒經(jīng)濟(jì)調(diào)度模型,通過(guò)協(xié)調(diào)優(yōu)化其他可調(diào)度分布式電源和儲(chǔ)能系統(tǒng)的出力,降低可再生能源不確定性帶來(lái)的風(fēng)險(xiǎn),從而保證系統(tǒng)運(yùn)行的穩(wěn)定性,但是只研究了聯(lián)絡(luò)線功率穩(wěn)定,且沒(méi)有考慮通信時(shí)延。針對(duì)時(shí)延系統(tǒng)的穩(wěn)定性分析,目前除了傳統(tǒng)的時(shí)域分析法[21]、頻域分析方法[22-23]外,基于Lyapunov函數(shù)的分析方法計(jì)算簡(jiǎn)單且可以得到更低保守性的穩(wěn)定性條件[24-25]。文獻(xiàn)[26]引入線性矩陣不等式(linear matrix inequality, LMI)直接沿系統(tǒng)軌跡列寫(xiě)Lyapunov泛函,利用自由權(quán)矩陣法提供列解Lyapunov函數(shù)和求解最優(yōu)牛頓-萊布尼茲公式權(quán)矩陣的系統(tǒng)性方法,大大降低了原有方法的保守性,將時(shí)延系統(tǒng)的穩(wěn)定性判據(jù)表示為一組線性矩陣不等式,由此得到求解時(shí)延穩(wěn)定性的系統(tǒng)性方法。文獻(xiàn)[27]同樣利用此系統(tǒng)性方法,分析了電力系統(tǒng)的時(shí)延穩(wěn)定性問(wèn)題,并針對(duì)時(shí)延系統(tǒng)的研究表明,采用自由權(quán)矩陣法所求得的系統(tǒng)狀態(tài)量的最大允許時(shí)延與其真值非常接近,所得結(jié)果的保守性低于其他方法。

    但是,目前尚無(wú)針對(duì)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)的時(shí)延穩(wěn)定性問(wèn)題的研究。ADN分布經(jīng)濟(jì)調(diào)度作為典型的復(fù)雜信息-物理耦合系統(tǒng),需要考慮分布式控制系統(tǒng)參數(shù)和通信等的影響,此外分布經(jīng)濟(jì)調(diào)度作為優(yōu)化問(wèn)題需要把經(jīng)濟(jì)成本納入系統(tǒng)方程,給時(shí)延穩(wěn)定性分析帶來(lái)了困難。本文以上述研究為基礎(chǔ),建立ADN分布式經(jīng)濟(jì)調(diào)度系統(tǒng)整體狀態(tài)方程,從控制角度對(duì)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)進(jìn)行穩(wěn)定性分析,使得系統(tǒng)在滿足經(jīng)濟(jì)性最優(yōu)目標(biāo)的同時(shí),也滿足系統(tǒng)持續(xù)穩(wěn)定運(yùn)行的要求;此外還將推導(dǎo)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)時(shí)延穩(wěn)定性判據(jù),探索改善分布控制模式下主動(dòng)配電網(wǎng)系統(tǒng)穩(wěn)定性,提高系統(tǒng)時(shí)延穩(wěn)定裕度的方法。

    本文首先以微增成本為一致性變量,設(shè)計(jì)了Leader-Follower模式下的主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度策略,考慮通信時(shí)延建立了主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度控制系統(tǒng)無(wú)時(shí)延、具有對(duì)稱時(shí)延和非對(duì)稱時(shí)延時(shí)的系統(tǒng)狀態(tài)模型。利用Lyapunov穩(wěn)定性定理和自由權(quán)矩陣方法得出具有更低保守性的系統(tǒng)時(shí)延穩(wěn)定判據(jù),同時(shí)利用線性矩陣不等式計(jì)算時(shí)延系統(tǒng)的穩(wěn)定裕度。最后通過(guò)IEEE-14節(jié)點(diǎn)配電網(wǎng)中的仿真算例驗(yàn)證本文分布式經(jīng)濟(jì)調(diào)度控制系統(tǒng)時(shí)延穩(wěn)定性分析和時(shí)延穩(wěn)定裕度計(jì)算方法的正確性和有效性。

    1 主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)設(shè)計(jì)

    1.1 分布式經(jīng)濟(jì)調(diào)度系統(tǒng)結(jié)構(gòu)

    主動(dòng)配電網(wǎng)經(jīng)濟(jì)調(diào)度系統(tǒng)需實(shí)現(xiàn)兩個(gè)目標(biāo):目標(biāo)1是配電網(wǎng)內(nèi)的所有DG總輸出有功功率加上從外部電網(wǎng)輸入的有功功率之和,需平衡配電網(wǎng)總負(fù)荷的有功功率需求;目標(biāo)2是在實(shí)現(xiàn)目標(biāo)1的基礎(chǔ)上使得配電網(wǎng)內(nèi)的DG總發(fā)電成本最小。本文采用如圖1所示分布式經(jīng)濟(jì)調(diào)度系統(tǒng)來(lái)實(shí)現(xiàn)上述目標(biāo)。

    圖1 主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)

    1.2 分布式經(jīng)濟(jì)調(diào)度模型

    本文采用的分布式經(jīng)濟(jì)調(diào)度模型的目標(biāo)函數(shù)為

    調(diào)度模型的等式約束條件為功率平衡,即

    當(dāng)PCC處實(shí)際功率與指令一致時(shí),表明系統(tǒng)達(dá)到功率平衡,等式約束(3)可改寫(xiě)為

    經(jīng)濟(jì)調(diào)度模型的線性不等式約束為

    在網(wǎng)損變化相對(duì)各DG功率變化較小時(shí),經(jīng)濟(jì)調(diào)度問(wèn)題式(1)、式(3)和式(5)為凸優(yōu)化問(wèn)題且具有唯一最優(yōu)解,且最優(yōu)解可表示為

    1.3 分布式經(jīng)濟(jì)調(diào)度系統(tǒng)控制模型

    本文采用分布式一致性協(xié)議,以分布式電源的微增成本為一致性變量實(shí)現(xiàn)分布式經(jīng)濟(jì)調(diào)度。

    本文分布式經(jīng)濟(jì)調(diào)度系統(tǒng)為多智能體系統(tǒng),其控制策略可表示為

    將式(7)以矩陣形式表示,得到該控制系統(tǒng)的狀態(tài)方程為

    考慮虛擬DG的分布式經(jīng)濟(jì)調(diào)度策略為

    根據(jù)式(10)得到分布式調(diào)度系統(tǒng)整體狀態(tài)方程為

    其中:

    若一個(gè)DG接收到其自身信息以及接收到其鄰居的信息均具有相同的時(shí)延,即對(duì)稱通信時(shí)延,此時(shí)ADN的系統(tǒng)空間狀態(tài)模型應(yīng)修改為

    實(shí)際中,DG接收自身信息通常不存在時(shí)延,而接收到其鄰居的信息具有時(shí)延,即不對(duì)稱通信時(shí)延,此時(shí)系統(tǒng)空間狀態(tài)模型為

    式(13)、式(14)和式(15)分別為分布式經(jīng)濟(jì)調(diào)度控制系統(tǒng)在無(wú)時(shí)延、對(duì)稱時(shí)延和不對(duì)稱時(shí)延情況下的狀態(tài)方程,相應(yīng)的控制框圖如圖2所示。

    圖2 分布式經(jīng)濟(jì)調(diào)度系統(tǒng)控制框圖

    2 系統(tǒng)穩(wěn)定性分析

    2.1 無(wú)時(shí)延系統(tǒng)穩(wěn)定性分析

    2.2 時(shí)延系統(tǒng)穩(wěn)定性分析

    本文基于Lyapunov穩(wěn)定性判斷方法,采用自由權(quán)矩陣,結(jié)合LMI數(shù)值分析法分析式(14)和式(15)所示時(shí)延系統(tǒng)的穩(wěn)定性,并求解時(shí)延穩(wěn)定裕度。以不對(duì)稱時(shí)延系統(tǒng)式(15)為例進(jìn)行說(shuō)明。

    首先利用二次型方法構(gòu)建式(15)所示時(shí)延系統(tǒng)的Lyapunov泛函為

    沿系統(tǒng)軌跡求Lyapunov泛函對(duì)時(shí)間的導(dǎo)數(shù)可得

    擴(kuò)展牛頓-萊布尼茲公式可得

    式中:

    其中:

    本文給出線性矩陣不等式表示的系統(tǒng)判穩(wěn)條件,進(jìn)一步采用LMI數(shù)值分析方法對(duì)判穩(wěn)條件進(jìn)行求解,以確定最優(yōu)的自由權(quán)矩陣系數(shù)和可接受的系統(tǒng)時(shí)延范圍。式(15)所示不對(duì)稱時(shí)延系統(tǒng)穩(wěn)定性判據(jù)為

    其中:

    3 仿真驗(yàn)證

    采用IEEE-14節(jié)點(diǎn)配電網(wǎng)系統(tǒng)進(jìn)行仿真計(jì)算,系統(tǒng)結(jié)構(gòu)和局部通信網(wǎng)絡(luò)如圖3所示。配電網(wǎng)共接入5個(gè)DG,分別位于節(jié)點(diǎn)2、3、7、10和12。整個(gè)網(wǎng)絡(luò)總有功負(fù)荷為28.7 MW,表1給出了各DG的具體信息。

    圖3 仿真系統(tǒng)拓?fù)浣Y(jié)構(gòu)及通信網(wǎng)絡(luò)

    表1 分布式電源參數(shù)

    3.1 系統(tǒng)穩(wěn)定性分析

    圖5 無(wú)時(shí)延通信故障和“即插即用”仿真結(jié)果

    圖6 無(wú)時(shí)延領(lǐng)導(dǎo)者通信故障仿真結(jié)果

    圖7 系統(tǒng)根軌跡變化圖

    不同反饋系數(shù)時(shí)的微增成本的對(duì)比如圖8所示。

    圖8 不同反饋系數(shù)時(shí)微增成本的比較

    3.2 系統(tǒng)的時(shí)延上界分析

    考慮系統(tǒng)存在通信時(shí)延,按照式(14)和式(15)分別得出系統(tǒng)不對(duì)稱以及對(duì)稱時(shí)延下的狀態(tài)方程,根據(jù)定理1,利用Matlab LMI Toolbox求解式(23),即可得出系統(tǒng)的時(shí)延穩(wěn)定裕度(反饋系數(shù)=0.5)。

    (1) 求解得出系統(tǒng)對(duì)稱時(shí)延下的穩(wěn)定裕度為0.2952 s,分別取時(shí)延為上界0.2952 s、超過(guò)時(shí)延上界的0.3 s和不超過(guò)時(shí)延上界的0.1 s和0.2 s,仿真得到分布式電源的微增成本的變化如圖9所示。

    圖9 不同對(duì)稱時(shí)延下微增成本

    (2) 求解此時(shí)系統(tǒng)不對(duì)稱時(shí)延下的穩(wěn)定裕度為0.5478 s,分別取時(shí)延為上界0.5478 s、超過(guò)時(shí)延上界的0.7 s和不超過(guò)時(shí)延上界的0.2 s和0.4 s,仿真得到DG微增成本的變化如圖10所示。

    (3) 當(dāng)采用文獻(xiàn)[20]使用的方法計(jì)算以上時(shí)延穩(wěn)定裕度時(shí),求解得出系統(tǒng)對(duì)稱和不對(duì)稱時(shí)延下的穩(wěn)定裕度分別為0.2877 s、0.5248 s,分別取時(shí)延為以上結(jié)果,仿真得到分布式電源的微增成本的變化如圖11所示。

    圖10 不同不對(duì)稱時(shí)延下微增成本變化

    圖11 文獻(xiàn)[20]方法時(shí)延裕度下微增成本變化

    相比本文方法,文獻(xiàn)[20]使用的方法所求的時(shí)延穩(wěn)定裕度較小,通過(guò)圖11對(duì)比圖9和圖10的仿真結(jié)果,本文所采用的方法所求得系統(tǒng)時(shí)延穩(wěn)定裕度更接近系統(tǒng)實(shí)際仿真的裕度,即本文的方法有效降低了結(jié)果的保守性。

    3.3 時(shí)延穩(wěn)定裕度的影響因素分析

    (1) 由圖5和圖6可以看出:當(dāng)系統(tǒng)通信網(wǎng)絡(luò)不連通時(shí),無(wú)法達(dá)成控制目標(biāo),也就會(huì)影響系統(tǒng)的穩(wěn)定性;而當(dāng)拓?fù)溥B通時(shí),通信故障等導(dǎo)致拓?fù)浒l(fā)生變化,節(jié)點(diǎn)的度數(shù)發(fā)生變化,雖然不會(huì)影響系統(tǒng)控制目標(biāo),但也對(duì)變化過(guò)程產(chǎn)生了影響。因此,改變系統(tǒng)的通信拓?fù)?,探究其?duì)系統(tǒng)時(shí)延穩(wěn)定裕度的影響。

    表2 不同通信拓?fù)湎碌膶?duì)稱時(shí)延穩(wěn)定裕度

    圖12 不同通信拓?fù)湎碌膶?duì)稱時(shí)延穩(wěn)定裕度

    分析表2和圖12可知,系統(tǒng)通信拓?fù)溥B通下,當(dāng)總度數(shù)逐步增大時(shí),系統(tǒng)的時(shí)延穩(wěn)定裕度均逐漸增大;相同度數(shù)下不同的連接方式導(dǎo)致的拓?fù)洳町悓?duì)系統(tǒng)時(shí)延穩(wěn)定裕度也有一定的影響。從提高系統(tǒng)時(shí)延穩(wěn)定裕度的角度出發(fā),應(yīng)該適當(dāng)增加系統(tǒng)的通信連接。

    圖6表明,盡管主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)的反饋系數(shù)不影響系統(tǒng)的穩(wěn)定性,但將影響系統(tǒng)收斂至穩(wěn)定狀態(tài)的過(guò)渡時(shí)間以及穩(wěn)定前的振蕩幅度;而3.2節(jié)中仿真結(jié)果表明系統(tǒng)具有不同的對(duì)稱和不對(duì)稱時(shí)延穩(wěn)定裕度。改變反饋系數(shù)得到系統(tǒng)對(duì)稱和不對(duì)稱時(shí)延穩(wěn)定裕度分別如表3和表4所示,圖13給出了兩種時(shí)延穩(wěn)定裕度的對(duì)比。

    分析表3、表4和圖13可知,當(dāng)反饋系數(shù)逐步增大時(shí),系統(tǒng)的對(duì)稱和不對(duì)稱時(shí)延穩(wěn)定裕度均逐漸減小;相同反饋系數(shù)下非對(duì)稱時(shí)延穩(wěn)定裕度要大于對(duì)稱時(shí)延下的穩(wěn)定裕度,但隨著反饋系數(shù)的增大,非對(duì)稱時(shí)延穩(wěn)定裕度減小的速度更快,當(dāng)>1時(shí),系統(tǒng)的非對(duì)稱與對(duì)稱時(shí)延穩(wěn)定裕度將基本相同。從提高系統(tǒng)時(shí)延穩(wěn)定裕度的角度出發(fā),分布式經(jīng)濟(jì)調(diào)度策略應(yīng)選擇較小的反饋系數(shù)。

    表3 不同反饋系數(shù)下的對(duì)稱時(shí)延穩(wěn)定裕度

    表4 不同反饋系數(shù)的不對(duì)稱時(shí)延穩(wěn)定裕度

    圖13 不同反饋系數(shù)下的系統(tǒng)時(shí)延穩(wěn)定裕度

    表4和圖13還表明,當(dāng)反饋系數(shù)較小時(shí),系統(tǒng)具有很大的非對(duì)稱時(shí)延裕度。圖14給出了較小時(shí)非對(duì)稱時(shí)延穩(wěn)定裕度的變化情況。

    圖14 反饋系數(shù)較小時(shí)的非對(duì)稱時(shí)延穩(wěn)定裕度

    圖15 很長(zhǎng)不對(duì)稱時(shí)延下的微增成本

    從圖15中可以看出,由于存在100 s時(shí)延,各DG在100 s內(nèi)將無(wú)法實(shí)現(xiàn)一致性,但各DG僅依靠自身的狀態(tài)信息也可滿足系統(tǒng)穩(wěn)定性的需求。在很長(zhǎng)的時(shí)間范圍內(nèi),各DG開(kāi)始接收到的其他DG的時(shí)延信息,經(jīng)過(guò)長(zhǎng)時(shí)間的調(diào)整后可以最終達(dá)成一致,實(shí)現(xiàn)最優(yōu)經(jīng)濟(jì)的目標(biāo)。

    綜上所述,從提高系統(tǒng)時(shí)延穩(wěn)定性的角度出發(fā),應(yīng)盡量為分布式經(jīng)濟(jì)調(diào)度控制策略式(10)選取較小的反饋系數(shù),以增大系統(tǒng)的時(shí)滯穩(wěn)定裕度;同時(shí)由于相同反饋系數(shù)下非對(duì)稱時(shí)延穩(wěn)定裕度要大于對(duì)稱時(shí)延下的穩(wěn)定裕度,因此應(yīng)盡可能減小各DG獲取自身狀態(tài)信息的時(shí)延,因?yàn)橄嗤答佅禂?shù)下非對(duì)稱時(shí)延穩(wěn)定裕度要大于對(duì)稱時(shí)延下的穩(wěn)定裕度。

    4 結(jié)論

    本文以微增成本為一致性變量,設(shè)計(jì)了Leader- Follower模式下的主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度策略,建立了考慮通信時(shí)延的主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度控制系統(tǒng)穩(wěn)定性分析模型。利用Lyapunov穩(wěn)定性定理和自由權(quán)矩陣方法得出系統(tǒng)考慮通信時(shí)延時(shí)的穩(wěn)定條件并求解時(shí)延穩(wěn)定裕度。仿真分析結(jié)果表明:

    (1) 本文提出的分布式經(jīng)濟(jì)調(diào)度控制策略可同時(shí)有效滿足主動(dòng)配電網(wǎng)功率平衡和發(fā)電成本最小的要求,且主動(dòng)配電網(wǎng)可保持穩(wěn)定運(yùn)行,可以應(yīng)對(duì)負(fù)荷變動(dòng)、DG“即插即用”、PCC功率指令變化以及通信拓?fù)溥B通下的通信故障等情況,提高了主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度策略的實(shí)用性與可靠性。

    (2) 本文給出的時(shí)延穩(wěn)定判據(jù)可有效分析主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)對(duì)稱和非對(duì)稱的時(shí)延穩(wěn)定性,以此為基礎(chǔ)利用LMI可求得系統(tǒng)保持穩(wěn)定的時(shí)延上界,降低了結(jié)果的保守性。當(dāng)實(shí)際時(shí)延小于穩(wěn)定裕度時(shí),隨著時(shí)延的增大各DG的微增成本進(jìn)入穩(wěn)定的時(shí)間將增加且波動(dòng)幅度增大,但系統(tǒng)可保持穩(wěn)定;而當(dāng)時(shí)延超過(guò)穩(wěn)定裕度時(shí)系統(tǒng)將失去穩(wěn)定性。

    (3) 分布式經(jīng)濟(jì)調(diào)度策略中通信拓?fù)浜头答佅禂?shù)將影響系統(tǒng)的時(shí)延穩(wěn)定裕度。當(dāng)拓?fù)溥B通的度數(shù)增大時(shí),系統(tǒng)時(shí)延穩(wěn)定裕度增大;而當(dāng)反饋系數(shù)增加時(shí),系統(tǒng)對(duì)稱和非對(duì)稱時(shí)延穩(wěn)定裕度均減小,反饋系數(shù)相同時(shí)系統(tǒng)不對(duì)稱時(shí)延穩(wěn)定裕度要大于對(duì)稱時(shí)延下的系統(tǒng)穩(wěn)定裕度。因此,從提高主動(dòng)配電網(wǎng)分布式經(jīng)濟(jì)調(diào)度系統(tǒng)穩(wěn)定性的角度而言,應(yīng)適當(dāng)增加系統(tǒng)的通信連接,選擇較小的反饋系數(shù),并盡可能減小DG檢測(cè)自身狀態(tài)信息的時(shí)延。

    本文研究中考慮的是具有二次項(xiàng)形式的分布式電源成本-功率函數(shù),而實(shí)際成本函數(shù)可能是非凸的,甚至是不連續(xù)的,給設(shè)計(jì)分布式經(jīng)濟(jì)調(diào)度策略帶來(lái)了困難。此外,分布式經(jīng)濟(jì)調(diào)度系統(tǒng)運(yùn)行時(shí)還需要考慮除了功率平衡約束和分布電源出力上下限約束之外的安全性約束,例如節(jié)點(diǎn)電壓約束和支路潮流約束等,如何在分布式框架下在經(jīng)濟(jì)調(diào)度模型中對(duì)這些全局性因素進(jìn)行處理也是需要進(jìn)一步深入研究的問(wèn)題。

    [1] 樂(lè)健, 周謙, 王曹, 等. 無(wú)功補(bǔ)償設(shè)備和分布式電源協(xié)同的配電網(wǎng)優(yōu)化控制策略研究[J]. 電力系統(tǒng)保護(hù)與控制, 2020, 48(18): 38-47.

    LE Jian, ZHOU Qian, WANG Cao, et al. Research on distributed optimal control strategy for a distribution network based on the cooperation of DGs and Var compensators[J]. Power System Protection and Control, 2020, 48(18): 38-47.

    [2] 張瑩, 孟潤(rùn)泉, 王子昂, 等. 一種基于一致性算法的改進(jìn)下垂控制策略[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(14): 104-111.

    ZHANG Ying, MENG Runquan, WANG Ziang, et al. An improved droop control strategy based on a consensus algorithm[J]. Power System Protection and Control, 2021, 49(14): 104-111.

    [3] 唐程輝, 張凡, 張寧, 等. 基于風(fēng)電場(chǎng)總功率條件分布的電力系統(tǒng)經(jīng)濟(jì)調(diào)度二次規(guī)劃方法[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(10): 2069-2078.

    TANG Chenghui, ZHANG Fan, ZHANG Ning, et al. Quadraticprogramming for power system economic dispatch based on the conditional probability distribution of wind farms sum power[J]. Transactions of China Electrotechnical Society, 2019, 34(10): 2069-2078.

    [4] 李曉明, 劉翔宇, 李安昌, 等. 配電網(wǎng)電壓控制的分布式光伏無(wú)功仿射可調(diào)魯棒優(yōu)化方法[J].電力系統(tǒng)保護(hù)與控制, 2021, 49(12): 124-131.

    LI Xiaoming, LIU Xiangyu, LI Anchang, et al. Distributed photovoltaic reactive power affine adjustable robust optimization method for voltage control of a distribution network[J]. Power System Protection and Control, 2021, 49(12): 124-131.

    [5] 李哲, 付明, 張曉燕, 等. 提高分布式儲(chǔ)能系統(tǒng)放電效率的分散協(xié)調(diào)優(yōu)化方法[J]. 電測(cè)與儀表, 2020, 57(22): 98-104, 119.

    LI Zhe, FU Ming, ZHANG Xiaoyan, et al. Decentralized cooperative optimization method enhancing discharging efficiency of distributed energy storage system[J]. Electrical Measurement & Instrumentation, 2020, 57(22): 98-104, 119.

    [6] 莊懷東, 吳紅斌, 劉海濤, 等. 含電動(dòng)汽車的微網(wǎng)系統(tǒng)多目標(biāo)經(jīng)濟(jì)調(diào)度[J]. 電工技術(shù)學(xué)報(bào), 2014, 29(增刊): 365-373.

    ZHUANG Huaidong, WU Hongbin, LIU Haitao, et al. Multi-objective economic dispatch of microgrid system considering electric vehicles[J]. Transactions of China Electrotechnical Society 2014, 29(S1): 365-373.

    [7] LI Q, GAO W, ZHANG H W, et al. Consensus based distributed economic dispatch control method in power systems[J]. IEEE Transactions on Smart Grid, 2019, 10(1): 941-954.

    [8] 李慧琴, 陳燕, 王崗紅, 等. 計(jì)及功率預(yù)測(cè)誤差的主動(dòng)配電網(wǎng)實(shí)時(shí)經(jīng)濟(jì)調(diào)度方法[J]. 電力系統(tǒng)保護(hù)與控制, 2020, 48(16): 100-107.

    LI Huiqin, CHEN Yan, WANG Ganghong, et al. Real-time economic dispatching method for active distribution networks considering power prediction errors[J]. Power System Protection and Control, 2020, 48(16): 100-107.

    [9] 張志榮, 邱曉燕, 孫旭, 等. 協(xié)調(diào)柔性負(fù)荷與儲(chǔ)能的交直流配電網(wǎng)經(jīng)濟(jì)優(yōu)化調(diào)度[J]. 電力建設(shè), 2020, 41(5): 116-123.

    ZHANG Zhirong, QIU Xiaoyan, SUN Xu, et al. Economic optimal scheduling of AC/DC distribution network coordinating flexible load and energy storage systems[J]. Electric Power Construction, 2020, 41(5): 116-123.

    [10] 張劉冬, 丁昊, 袁曉冬, 等. 考慮價(jià)格需求響應(yīng)的主動(dòng)配電網(wǎng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度[J]. 電力工程技術(shù), 2017, 36(4): 31-35, 42.

    ZHANG Liudong, DING Hao, YUAN Xiaodong, et al. Active and reactive power coordinated economic dispatch of active distribution networks with consideration of price-based demand response[J]. Electric Power Engineering Technology, 2017, 36(4): 31-35, 42.

    [11] 溫廣輝, 余星火, 劉智偉. 智能電網(wǎng)中分布式經(jīng)濟(jì)調(diào)度研究進(jìn)展: 綜述(英文)[J]. Frontiers of Information Technology & Electronic Engineering, 2021, 22(1): 25-40.

    WEN Guanghui, YU Xinghuo, LIU Zhiwei. Recent progress on the study of distributed economic dispatch in smart grid: an overview[J]. Frontiers of Information Technology & Electronic Engineering, 2021, 22(1): 25-40.

    [12] 樂(lè)健, 周謙, 趙聯(lián)港, 等. 考慮個(gè)體欺騙的有源配電網(wǎng)完全分布式經(jīng)濟(jì)調(diào)度策略[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2020, 40(17): 5445-5454.

    LE Jian, ZHOU Qian, ZHAO Liangang, et al. Fully distributed economic dispatch of active distribution network considering individual cheating[J]. Proceedings of the CSEE, 2020, 40(17): 5445-5454.

    [13] YANG S, TAN S, XU J X. Consensus based approach for economic dispatch problem in a smart grid[J]. IEEE Transactions on Power Systems, 2013, 28(4): 4416-4426.

    [14] LE J, ZHOU Q, ZHAO L, et al. Research on the real-time distributed economic dispatch strategy for microgrids[J]. Energies, 2019, 12(20): 1-16.

    [15] YANG T, WU D, SUN Y, et al. Impacts of time delays on distributed algorithms for economic dispatch[C] // 2015 IEEE Power & Energy Society General Meeting, July 26-28, 2015, Denver CO: 1-5.

    [16] CHEN G, ZHAO Z. Delay effects on consensus-based distributed economic dispatch algorithm in microgrid[J]. IEEE Transactions on Power Systems, 2018, 33(1): 602-612.

    [17] HUANG B, LIU L, ZHANG H, et al. Distributed optimal economic dispatch for microgrids considering communication delays[J]. IEEE Transactions on Systems, Man, and Cybernetics: Systems, 2019, 49(8): 1634-1642.

    [18] 聞秉科. 考慮穩(wěn)定性的微電網(wǎng)經(jīng)濟(jì)運(yùn)行優(yōu)化策略研究[D]. 武漢: 武漢大學(xué), 2016.

    WEN Bingke. Economic operation optimization strategy of microgrid considering stability[D]. Wuhan: Wuhan University, 2016.

    [19] XU Yan. A review of cyber security risks of power systems: from static to dynamic false data attacks[J]. Protection and Control of Modern Power Systems, 2020, 5(3): 190-201.

    [20] 康慨, 鄧少平, 張超, 等. 考慮聯(lián)絡(luò)線功率穩(wěn)定的微電網(wǎng)魯棒經(jīng)濟(jì)調(diào)度模型[J]. 中國(guó)電力, 2019, 52(6): 60-67.

    KANG Kai, DENG Shaoping, ZHANG Chao, et al. Robust economic dispatch model of interconnected microgrid considering tie line power stability[J]. Electric Power, 2019, 52(6): 60-67.

    [21] NI B, XIAO D, SHAH S L. Time delay estimation for MIMO dynamical systems with time-frequency domain analysis[J]. Journal of Process Control, 2010, 20(1): 83-94.

    [22] 陳萌, 肖湘寧, 羅超. 基于虛擬同步發(fā)電機(jī)的微電網(wǎng)延時(shí)補(bǔ)償二次頻率控制[J]. 電工技術(shù)學(xué)報(bào), 2018, 33(16): 3845-3854.

    CHEN Meng, XIAO Xiangning, LUO Chao. Secondary frequency control including delay compensation in microgrids based on virtual synchronous generator[J]. Transactions of China Electrotechnical Society, 2018, 33(16): 3845-3854.

    [23] 賈宏杰, 安海云, 余曉丹. 電力系統(tǒng)時(shí)滯依賴型魯棒穩(wěn)定判據(jù)及其應(yīng)用[J]. 電力系統(tǒng)自動(dòng)化, 2010, 34(3): 6-11.

    JIA Hongjie, AN Haiyun, YU Xiaodan. A delay-dependent robust stability criterion for power system and its application[J]. Automation of Electric Power Systems, 2010, 34(3): 6-11.

    [24] 錢偉, 王晨晨, 費(fèi)樹(shù)岷. 區(qū)間變時(shí)滯廣域電力系統(tǒng)穩(wěn)定性分析與控制器設(shè)計(jì)[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(17): 3640-3650.

    QIAN Wei, WANG Chenchen, FEI Shumin. Stability analysis and controller design of wide area power system with time varying delay[J]. Transactions of China Electrotechnical Society, 2019, 34(17): 3640-3650.

    [25] 丁貴立, 林濤, 陳汝斯, 等. 抗飽和廣域時(shí)滯阻尼控制器參量Lyapunov設(shè)計(jì)策略[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(20): 181-191.

    DING Guili, LIN Tao, CHEN Rusi, et al. Research on design method of anti-windup wide-area time-delay damping controller[J]. Transactions of China Electrotechnical Society, 2019, 34(20): 181-191.

    [26] HE Y, WU M, SHE J H, et al. Parameter-dependent Lyapunov functional for stability of time-delay systems with polytopic-type uncertainties[J]. IEEE Transactions on Automatic Control, 2004, 49(5): 828-832.

    [27] YAO W, JIANG L, WU Q H, et al. Delay-dependent stability analysis of the power system with a wide-area damping controller embedded[J]. IEEE Transactions on Power Systems, 2011, 26(1): 233-240.

    Time-delay stability analysis of an active distribution network adopting a distributed economic dispatch strategy

    LE Jian1, QI Gan1, ZHAO Liangang1, LIAO Xiaobing2

    (1. School of Electrical Engineering and Automation, Wuhan University, Wuhan 430072, China;2. Wuhan Institute of Technology, Wuhan 430205, China)

    The distributed economic dispatch model has advantages of better robustness, faster computation and less communication, and is suitable for active distribution networks with a high penetration of distributed power sources. However, the time delay in communication between units in the distributed model will likely undermine the stability of the system. So in this paper the time delay stability of active distribution networks adopting a distributed economic dispatching strategy is studied. First, a fully distributed economic dispatch strategy is designed by introducing virtual distributed power sources. Then a stability analysis model under symmetric and asymmetric time delay of active distribution networks using this strategy is established. The Lyapunov stability theorem and the free power matrix method are also used to derive a lower conservative criterion for the time delay stability of the system, and the stability margin of the time delay system is calculated by means of linear matrix inequalities. Finally, the correctness and validity of the delay stability analysis and delay stability margin calculation methods of the system are verified by simulation examples in the IEEE-14 node distribution network.

    distributed economic dispatch; communication delay; delay stability margin; Lyapunov stability criterion; LMI

    10.19783/j.cnki.pspc.220010

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

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

    2022-01-02;

    2022-03-17

    樂(lè) 健(1975—),男,博士,副教授,研究方向?yàn)殡娔苜|(zhì)量問(wèn)題及其控制技術(shù);E-mail: lej01@tsinghua.org.cn

    綦 淦(1999—),男,通信作者,碩士,研究方向?yàn)殡娏ο到y(tǒng)經(jīng)濟(jì)調(diào)度及優(yōu)化控制。E-mail: 1679264596@qq.com

    (編輯 魏小麗)

    猜你喜歡
    裕度時(shí)延分布式
    基于GCC-nearest時(shí)延估計(jì)的室內(nèi)聲源定位
    電子制作(2019年23期)2019-02-23 13:21:12
    基于改進(jìn)二次相關(guān)算法的TDOA時(shí)延估計(jì)
    分布式光伏熱錢洶涌
    能源(2017年10期)2017-12-20 05:54:07
    分布式光伏:爆發(fā)還是徘徊
    能源(2017年5期)2017-07-06 09:25:54
    FRFT在水聲信道時(shí)延頻移聯(lián)合估計(jì)中的應(yīng)用
    基于DFIG可用無(wú)功裕度的風(fēng)電場(chǎng)無(wú)功電壓控制方法
    基于分段CEEMD降噪的時(shí)延估計(jì)研究
    基于DDS的分布式三維協(xié)同仿真研究
    三環(huán)路核電廠的抗震裕度評(píng)價(jià)
    基于ANN模型的在線電壓穩(wěn)定裕度評(píng)估
    亚洲av美国av| 精品一区二区三区视频在线观看免费| 国产精品一区二区三区四区免费观看 | 亚洲精品国产av成人精品 | 欧美绝顶高潮抽搐喷水| 亚洲最大成人中文| 99久久久亚洲精品蜜臀av| 国产精品福利在线免费观看| 变态另类丝袜制服| 性插视频无遮挡在线免费观看| 亚洲欧美中文字幕日韩二区| 久久久久精品国产欧美久久久| 日韩欧美在线乱码| 亚洲电影在线观看av| 日韩一本色道免费dvd| 亚洲美女视频黄频| 搡老妇女老女人老熟妇| 最近的中文字幕免费完整| 精品久久久久久成人av| 免费av观看视频| 黄色日韩在线| 搡老岳熟女国产| 一a级毛片在线观看| 变态另类丝袜制服| 熟妇人妻久久中文字幕3abv| 国内久久婷婷六月综合欲色啪| 村上凉子中文字幕在线| 国产精品福利在线免费观看| av黄色大香蕉| 一级毛片aaaaaa免费看小| 小蜜桃在线观看免费完整版高清| 久久久精品大字幕| 亚洲无线观看免费| 悠悠久久av| 99久久久亚洲精品蜜臀av| 欧美日本亚洲视频在线播放| 国产国拍精品亚洲av在线观看| 国产精品女同一区二区软件| 国产片特级美女逼逼视频| 欧美丝袜亚洲另类| videossex国产| 国产亚洲av嫩草精品影院| 美女 人体艺术 gogo| 亚洲美女搞黄在线观看 | 最近视频中文字幕2019在线8| 日日干狠狠操夜夜爽| 亚洲四区av| 精品久久久久久久末码| 久久久国产成人精品二区| 国产 一区 欧美 日韩| 午夜日韩欧美国产| 在线天堂最新版资源| 蜜臀久久99精品久久宅男| 精品一区二区免费观看| 色综合站精品国产| 欧美zozozo另类| 能在线免费观看的黄片| 国产综合懂色| 99久久九九国产精品国产免费| 1000部很黄的大片| 国产蜜桃级精品一区二区三区| 国产国拍精品亚洲av在线观看| 日韩av在线大香蕉| 97超级碰碰碰精品色视频在线观看| 级片在线观看| 久久久久国产网址| 国产三级在线视频| 久久久久久久午夜电影| 亚洲精品国产av成人精品 | 亚洲精品日韩av片在线观看| 少妇的逼好多水| 女人十人毛片免费观看3o分钟| av在线播放精品| 精品一区二区三区视频在线| 国产乱人偷精品视频| 亚洲欧美成人精品一区二区| 人妻夜夜爽99麻豆av| 看非洲黑人一级黄片| 亚洲熟妇熟女久久| 一区二区三区四区激情视频 | 国产精品伦人一区二区| 久久久欧美国产精品| 国产视频一区二区在线看| 亚洲自拍偷在线| 草草在线视频免费看| 久久精品综合一区二区三区| 毛片女人毛片| 国产乱人偷精品视频| 亚洲精华国产精华液的使用体验 | av天堂中文字幕网| 亚洲欧美日韩无卡精品| 少妇的逼好多水| 久久天躁狠狠躁夜夜2o2o| 丰满乱子伦码专区| 欧美xxxx黑人xx丫x性爽| 一a级毛片在线观看| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图| 深夜a级毛片| 精品少妇黑人巨大在线播放 | 成人性生交大片免费视频hd| 在线观看美女被高潮喷水网站| 九九久久精品国产亚洲av麻豆| 老熟妇仑乱视频hdxx| 网址你懂的国产日韩在线| 欧美一区二区精品小视频在线| 少妇被粗大猛烈的视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲在线观看片| 好男人在线观看高清免费视频| 国产精品综合久久久久久久免费| 亚洲在线观看片| 夜夜夜夜夜久久久久| 国产成人a区在线观看| 国产免费男女视频| 亚洲精品在线观看二区| 女生性感内裤真人,穿戴方法视频| 免费人成视频x8x8入口观看| 老司机午夜福利在线观看视频| 人妻丰满熟妇av一区二区三区| 大香蕉久久网| 日本与韩国留学比较| 精品午夜福利在线看| 日日干狠狠操夜夜爽| 三级毛片av免费| 精品久久久久久久久久免费视频| 97碰自拍视频| 免费人成视频x8x8入口观看| 国内精品久久久久精免费| 国产女主播在线喷水免费视频网站 | 亚洲人成网站在线播放欧美日韩| 国产一区亚洲一区在线观看| 久久久久性生活片| 色综合色国产| 成人三级黄色视频| 亚洲va在线va天堂va国产| 十八禁国产超污无遮挡网站| 22中文网久久字幕| 日韩成人av中文字幕在线观看 | 一进一出好大好爽视频| 麻豆国产av国片精品| 欧美潮喷喷水| 亚洲av中文av极速乱| 亚洲av免费在线观看| 校园人妻丝袜中文字幕| 国产精华一区二区三区| 丝袜喷水一区| 成人亚洲欧美一区二区av| 亚洲自拍偷在线| 一级毛片电影观看 | 久久人妻av系列| 婷婷亚洲欧美| av国产免费在线观看| 久久久久久大精品| 韩国av在线不卡| 赤兔流量卡办理| 不卡一级毛片| 免费av毛片视频| 岛国在线免费视频观看| 晚上一个人看的免费电影| 欧美人与善性xxx| 男女啪啪激烈高潮av片| 日韩强制内射视频| 99久久无色码亚洲精品果冻| 日本三级黄在线观看| 熟妇人妻久久中文字幕3abv| 97人妻精品一区二区三区麻豆| 亚洲性夜色夜夜综合| 在线免费十八禁| 久久精品国产99精品国产亚洲性色| 国产探花在线观看一区二区| 精品熟女少妇av免费看| 亚洲av美国av| 麻豆成人午夜福利视频| 99热这里只有是精品50| 免费在线观看影片大全网站| 色哟哟·www| 男人的好看免费观看在线视频| 桃色一区二区三区在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区精品小视频在线| 99久久精品热视频| 深爱激情五月婷婷| 免费在线观看成人毛片| 亚洲第一电影网av| 大香蕉久久网| 日韩三级伦理在线观看| 在现免费观看毛片| 亚洲第一电影网av| 欧美三级亚洲精品| 精品乱码久久久久久99久播| 亚洲欧美清纯卡通| 亚洲三级黄色毛片| 国产极品精品免费视频能看的| 我的女老师完整版在线观看| 国产黄色视频一区二区在线观看 | 我要搜黄色片| 午夜福利成人在线免费观看| h日本视频在线播放| 成人特级黄色片久久久久久久| 国产伦精品一区二区三区视频9| 69av精品久久久久久| 亚洲人成网站在线播| 欧美绝顶高潮抽搐喷水| 最新中文字幕久久久久| 99久久无色码亚洲精品果冻| 亚洲国产精品国产精品| 寂寞人妻少妇视频99o| 十八禁网站免费在线| 久久午夜亚洲精品久久| 欧美最黄视频在线播放免费| 97碰自拍视频| 91在线精品国自产拍蜜月| 欧洲精品卡2卡3卡4卡5卡区| 黄色欧美视频在线观看| 国产极品精品免费视频能看的| 亚洲av美国av| 好男人在线观看高清免费视频| 亚洲丝袜综合中文字幕| 国产视频一区二区在线看| 免费不卡的大黄色大毛片视频在线观看 | 婷婷六月久久综合丁香| 97热精品久久久久久| 久久精品国产亚洲网站| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 一边摸一边抽搐一进一小说| 日本三级黄在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲美女视频黄频| 久久久色成人| 亚洲在线观看片| 老熟妇仑乱视频hdxx| 日韩欧美三级三区| 狂野欧美白嫩少妇大欣赏| 国产在线男女| 美女大奶头视频| 欧美激情在线99| 在线播放国产精品三级| 精品不卡国产一区二区三区| 亚洲精品粉嫩美女一区| 精品无人区乱码1区二区| 亚洲,欧美,日韩| 日本色播在线视频| 最近的中文字幕免费完整| 两个人视频免费观看高清| 亚洲精品在线观看二区| 日本与韩国留学比较| 亚洲七黄色美女视频| 久久久久久久亚洲中文字幕| 亚洲精品成人久久久久久| 1024手机看黄色片| 欧美最新免费一区二区三区| 一级毛片久久久久久久久女| 淫妇啪啪啪对白视频| 成人三级黄色视频| 秋霞在线观看毛片| 亚洲经典国产精华液单| 校园人妻丝袜中文字幕| 久久久欧美国产精品| www.色视频.com| 国产日本99.免费观看| 中文亚洲av片在线观看爽| 一级毛片久久久久久久久女| 午夜日韩欧美国产| 高清日韩中文字幕在线| 日本黄色片子视频| 久久精品人妻少妇| 99久国产av精品国产电影| 一个人看视频在线观看www免费| 精品乱码久久久久久99久播| 九九爱精品视频在线观看| 此物有八面人人有两片| 久久久色成人| 九九在线视频观看精品| 免费搜索国产男女视频| 国产精品乱码一区二三区的特点| 乱人视频在线观看| 日韩一区二区视频免费看| 色综合亚洲欧美另类图片| 在线免费十八禁| 欧美另类亚洲清纯唯美| 中文资源天堂在线| 一个人看的www免费观看视频| 国产精品爽爽va在线观看网站| 国产大屁股一区二区在线视频| 国产精品无大码| 国产人妻一区二区三区在| 亚洲av二区三区四区| 中国国产av一级| 又黄又爽又刺激的免费视频.| 18禁黄网站禁片免费观看直播| 成人av在线播放网站| 在线观看66精品国产| 少妇裸体淫交视频免费看高清| 青春草视频在线免费观看| 成人国产麻豆网| 久久99热6这里只有精品| 身体一侧抽搐| 美女cb高潮喷水在线观看| 国产精品国产高清国产av| 欧美性猛交黑人性爽| 久久久a久久爽久久v久久| 黄片wwwwww| 一进一出抽搐动态| 国产高潮美女av| 中国美女看黄片| 搞女人的毛片| 国产精品精品国产色婷婷| 大又大粗又爽又黄少妇毛片口| 色哟哟哟哟哟哟| 国产成人91sexporn| 3wmmmm亚洲av在线观看| 搡女人真爽免费视频火全软件 | 99热6这里只有精品| 午夜激情欧美在线| www日本黄色视频网| 精品一区二区三区人妻视频| 欧美人与善性xxx| 男女啪啪激烈高潮av片| 国产精品av视频在线免费观看| 欧美一区二区亚洲| 国产成人aa在线观看| 一级毛片aaaaaa免费看小| 国产午夜福利久久久久久| 国产精品av视频在线免费观看| 免费人成视频x8x8入口观看| 男人狂女人下面高潮的视频| 中文资源天堂在线| 12—13女人毛片做爰片一| 18禁在线播放成人免费| 国产精品人妻久久久久久| 国产亚洲91精品色在线| 精品久久久久久久久久免费视频| 久久中文看片网| 亚洲人成网站在线播| 国产午夜精品久久久久久一区二区三区 | 中文字幕熟女人妻在线| 看片在线看免费视频| 国产三级在线视频| 欧美一区二区国产精品久久精品| 亚洲精品国产av成人精品 | 色5月婷婷丁香| 中文亚洲av片在线观看爽| 日韩欧美一区二区三区在线观看| 天天躁日日操中文字幕| 欧美高清性xxxxhd video| 特级一级黄色大片| 亚洲欧美日韩无卡精品| 免费看美女性在线毛片视频| 欧美3d第一页| 日本熟妇午夜| 欧美bdsm另类| 1024手机看黄色片| 免费观看在线日韩| 欧美xxxx黑人xx丫x性爽| 日韩,欧美,国产一区二区三区 | 少妇的逼好多水| 真实男女啪啪啪动态图| 日韩欧美三级三区| 老女人水多毛片| 久久99热这里只有精品18| 一卡2卡三卡四卡精品乱码亚洲| 日本色播在线视频| 亚洲av不卡在线观看| 国产精品无大码| 免费看av在线观看网站| 国产一区二区三区av在线 | 内射极品少妇av片p| 乱人视频在线观看| 精品少妇黑人巨大在线播放 | 亚洲精品粉嫩美女一区| av专区在线播放| 变态另类丝袜制服| 日日摸夜夜添夜夜添小说| 国产亚洲91精品色在线| a级一级毛片免费在线观看| 国产一区二区激情短视频| 国产色婷婷99| 亚洲内射少妇av| 国产av麻豆久久久久久久| 变态另类丝袜制服| 男人和女人高潮做爰伦理| 精品不卡国产一区二区三区| 欧美xxxx黑人xx丫x性爽| 人人妻人人澡人人爽人人夜夜 | 在线国产一区二区在线| 嫩草影院精品99| 亚洲欧美成人精品一区二区| av免费在线看不卡| 精品久久久久久久末码| 又爽又黄无遮挡网站| 少妇人妻精品综合一区二区 | 麻豆国产av国片精品| 啦啦啦韩国在线观看视频| 两个人的视频大全免费| 国产av不卡久久| 俄罗斯特黄特色一大片| 国产成人影院久久av| 色av中文字幕| 国产单亲对白刺激| 老女人水多毛片| 看片在线看免费视频| 自拍偷自拍亚洲精品老妇| 两个人视频免费观看高清| 久久久久精品国产欧美久久久| 悠悠久久av| 欧美性猛交╳xxx乱大交人| 日本黄色片子视频| 日韩中字成人| 亚洲图色成人| 日本色播在线视频| 免费观看的影片在线观看| 成人一区二区视频在线观看| 国产伦精品一区二区三区视频9| 99久久无色码亚洲精品果冻| 天堂网av新在线| 国产亚洲精品久久久久久毛片| 在线观看av片永久免费下载| 女人十人毛片免费观看3o分钟| 午夜精品在线福利| 在线国产一区二区在线| 午夜福利高清视频| 综合色av麻豆| 亚洲欧美日韩高清专用| 久久精品国产清高在天天线| 国产精品久久久久久av不卡| 91在线观看av| 97超视频在线观看视频| 日韩高清综合在线| 波多野结衣巨乳人妻| 97人妻精品一区二区三区麻豆| 毛片一级片免费看久久久久| 亚洲图色成人| 国产精品国产高清国产av| 亚洲成人中文字幕在线播放| 丝袜美腿在线中文| 久久久久久久久久黄片| 日韩制服骚丝袜av| 国产亚洲精品久久久com| 美女高潮的动态| 欧美一级a爱片免费观看看| 亚洲人成网站在线播放欧美日韩| 国产极品精品免费视频能看的| 国产一区二区亚洲精品在线观看| 亚洲美女搞黄在线观看 | 嫩草影院入口| 国产成人freesex在线 | 国产精品国产高清国产av| 嫩草影院精品99| 91在线观看av| 天天躁日日操中文字幕| 永久网站在线| 三级毛片av免费| 亚洲成人精品中文字幕电影| 国内久久婷婷六月综合欲色啪| 亚州av有码| 午夜日韩欧美国产| 69av精品久久久久久| 国产熟女欧美一区二区| 熟女人妻精品中文字幕| 国产真实伦视频高清在线观看| 国产精品人妻久久久影院| 午夜免费激情av| ponron亚洲| 美女高潮的动态| 亚洲美女视频黄频| 99精品在免费线老司机午夜| 天美传媒精品一区二区| 国产熟女欧美一区二区| 99热6这里只有精品| 国产一区二区在线观看日韩| 淫妇啪啪啪对白视频| 国产精品日韩av在线免费观看| 日韩国内少妇激情av| 久久午夜福利片| 91久久精品电影网| 久久久国产成人精品二区| 日本 av在线| 久久久精品大字幕| 我要搜黄色片| 女同久久另类99精品国产91| 午夜激情欧美在线| 久久亚洲精品不卡| 美女被艹到高潮喷水动态| 综合色丁香网| 国产精品亚洲一级av第二区| 非洲黑人性xxxx精品又粗又长| 97人妻精品一区二区三区麻豆| 乱码一卡2卡4卡精品| 亚洲无线在线观看| 大香蕉久久网| 日韩亚洲欧美综合| 久久人人精品亚洲av| 国产av一区在线观看免费| 变态另类丝袜制服| 久久久久国产网址| 日韩亚洲欧美综合| 97人妻精品一区二区三区麻豆| 国产 一区 欧美 日韩| 成年女人永久免费观看视频| 国产精品国产高清国产av| 在现免费观看毛片| 亚洲自拍偷在线| 老师上课跳d突然被开到最大视频| 亚洲国产欧美人成| 又黄又爽又免费观看的视频| 精品一区二区三区视频在线| 亚洲精品成人久久久久久| 亚洲在线自拍视频| 校园人妻丝袜中文字幕| 亚洲一区二区三区色噜噜| 欧美丝袜亚洲另类| 亚洲综合色惰| 午夜日韩欧美国产| 一边摸一边抽搐一进一小说| 久久久午夜欧美精品| 亚洲人成网站在线播| 变态另类丝袜制服| 老熟妇乱子伦视频在线观看| 熟女人妻精品中文字幕| 日本 av在线| 99久久精品一区二区三区| av视频在线观看入口| 国产成人精品久久久久久| 桃色一区二区三区在线观看| 亚洲一级一片aⅴ在线观看| а√天堂www在线а√下载| 亚洲精品一卡2卡三卡4卡5卡| 日韩 亚洲 欧美在线| 男人舔女人下体高潮全视频| 国产精品电影一区二区三区| 成人午夜高清在线视频| 亚洲成人久久性| 老司机福利观看| 五月玫瑰六月丁香| 国产成人freesex在线 | 91av网一区二区| 精品一区二区三区av网在线观看| 校园人妻丝袜中文字幕| 免费在线观看影片大全网站| 欧美极品一区二区三区四区| 色哟哟·www| 日韩精品中文字幕看吧| 国产淫片久久久久久久久| 男人舔女人下体高潮全视频| 看片在线看免费视频| 天天躁日日操中文字幕| 久久欧美精品欧美久久欧美| 亚洲第一电影网av| 俄罗斯特黄特色一大片| 日本-黄色视频高清免费观看| 成年女人毛片免费观看观看9| 亚洲av美国av| 日本熟妇午夜| 精品人妻熟女av久视频| 久久精品国产自在天天线| 国产激情偷乱视频一区二区| 偷拍熟女少妇极品色| 人人妻人人澡人人爽人人夜夜 | 亚洲第一区二区三区不卡| 成人二区视频| 五月玫瑰六月丁香| 婷婷亚洲欧美| 国产美女午夜福利| 日韩欧美一区二区三区在线观看| 久久久久久久久中文| 亚洲va在线va天堂va国产| 日韩欧美精品v在线| 99热精品在线国产| 岛国在线免费视频观看| 国内少妇人妻偷人精品xxx网站| 直男gayav资源| or卡值多少钱| 亚洲电影在线观看av| 欧美成人精品欧美一级黄| 精品久久国产蜜桃| 精品不卡国产一区二区三区| 欧美中文日本在线观看视频| 十八禁网站免费在线| 人妻制服诱惑在线中文字幕| 国产精品一区www在线观看| 亚洲欧美成人综合另类久久久 | 国产伦精品一区二区三区四那| 精品久久久久久久末码| 村上凉子中文字幕在线| 联通29元200g的流量卡| 色播亚洲综合网| 日产精品乱码卡一卡2卡三| 精品久久久久久久久久免费视频| 女人十人毛片免费观看3o分钟| 嫩草影视91久久| 午夜久久久久精精品| 综合色丁香网| 国产女主播在线喷水免费视频网站 | 亚洲电影在线观看av| 老司机福利观看| a级毛片a级免费在线| 欧美日韩一区二区视频在线观看视频在线 | 人妻少妇偷人精品九色| 婷婷六月久久综合丁香| 亚洲国产日韩欧美精品在线观看| 好男人在线观看高清免费视频| 成人国产麻豆网| 欧美日韩国产亚洲二区| 在线免费十八禁| 香蕉av资源在线| 久久久久国产网址| 麻豆乱淫一区二区| 国产男人的电影天堂91| 亚洲色图av天堂| 欧美色欧美亚洲另类二区| 国产精品一区二区三区四区久久| 黄色日韩在线| 美女高潮的动态|