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

    直流微電網(wǎng)集群的大信號穩(wěn)定性分析

    2022-06-25 07:02:42劉宿城秦強(qiáng)棟夏夢宇劉曉東
    電工技術(shù)學(xué)報(bào) 2022年12期
    關(guān)鍵詞:信號系統(tǒng)

    劉宿城 李 響 秦強(qiáng)棟 夏夢宇 劉曉東

    直流微電網(wǎng)集群的大信號穩(wěn)定性分析

    劉宿城1,2李 響1,2秦強(qiáng)棟1,2夏夢宇1,2劉曉東1,2

    (1. 安徽工業(yè)大學(xué)電氣與信息工程學(xué)院 馬鞍山 243000 2. 安徽工業(yè)大學(xué)電力電子與運(yùn)動控制安徽省重點(diǎn)實(shí)驗(yàn)室 馬鞍山 243000)

    直流微電網(wǎng)集群通常由多個(gè)直流微電網(wǎng)互聯(lián)而成,通過靈活的功率流動控制實(shí)現(xiàn)區(qū)域能源共享和優(yōu)化利用,以充分發(fā)揮直流分布式發(fā)電系統(tǒng)的優(yōu)勢。然而,小規(guī)模直流微電網(wǎng)具有低慣性和高阻抗的“弱電網(wǎng)”特性,“弱-弱”互聯(lián)會降低集群系統(tǒng)的阻尼,甚至出現(xiàn)振蕩或系統(tǒng)崩潰等嚴(yán)重后果。同時(shí),直流微電網(wǎng)集群的高階、強(qiáng)耦合以及非線性動態(tài)特性也對其穩(wěn)定性分析帶來了巨大挑戰(zhàn)。為此,基于Brayton-Moser混合勢理論,提出針對直流微電網(wǎng)集群的大信號穩(wěn)定性分析方法。建立集群的大信號降階模型,詳細(xì)推導(dǎo)用于系統(tǒng)大信號穩(wěn)定性判據(jù)的混合勢函數(shù),并分析關(guān)鍵參數(shù)對穩(wěn)定區(qū)間的影響,最后通過實(shí)時(shí)仿真結(jié)果驗(yàn)證了該分析的正確性。

    直流微電網(wǎng)集群 大信號穩(wěn)定性 恒功率負(fù)載 Brayton-Moser混合勢理論

    0 引言

    隨著直流功率變換技術(shù)的發(fā)展,基于直流技術(shù)的電力應(yīng)用得到了復(fù)蘇,尤其是在分布式發(fā)電領(lǐng) 域[1]。在此背景下,直流微電網(wǎng)(DC Microgrids, DCMGs)作為整合可再生能源的有效解決方案應(yīng)運(yùn)而生。與交流微電網(wǎng)相比,直流微電網(wǎng)因具有高效率、低成本、控制設(shè)計(jì)簡單等諸多優(yōu)點(diǎn)而受到越來越多的關(guān)注[2-3]。為進(jìn)一步發(fā)揮分布式發(fā)電系統(tǒng)的優(yōu)勢,地理上毗鄰的多個(gè)直流微電網(wǎng)可互聯(lián)構(gòu)成直流微電網(wǎng)集群(DCMG Clusters, DCMGC),通過網(wǎng)間靈活的功率流動控制實(shí)現(xiàn)電源系統(tǒng)高彈性、高可靠性、經(jīng)濟(jì)性以及最優(yōu)負(fù)載響應(yīng)等目標(biāo)[4-6]。

    然而,由于直流微電網(wǎng)通常由分布式可再生能源(如光伏和風(fēng)電)、儲能系統(tǒng)(如蓄電池和超級電容)以及各類負(fù)載等異構(gòu)單元通過電力電子接口連接至母線,其功率等級為kW級的小規(guī)模電源系統(tǒng),導(dǎo)致系統(tǒng)等效慣量較小、網(wǎng)絡(luò)阻抗較大,易受各類暫態(tài)擾動(如新能源的間歇性、負(fù)載投切、工作模式的切換及短路故障等)的影響,為典型的弱電網(wǎng)特性[7-11]。從這個(gè)意義上講,直流微電網(wǎng)集群可看作為多個(gè)直流微電網(wǎng)的“弱-弱”互聯(lián),系統(tǒng)阻尼將因此進(jìn)一步降低,可能導(dǎo)致振蕩或系統(tǒng)崩潰,嚴(yán)重威脅系統(tǒng)的安全穩(wěn)定運(yùn)行。同時(shí),直流微電網(wǎng)集群的動態(tài)分析涉及眾多狀態(tài)變量及多種非線性因素,如開關(guān)行為和恒功率負(fù)載(Constant Power Load, CPL)等,其高階動態(tài)非線性特性也對集群系統(tǒng)的穩(wěn)定性分析帶來了巨大挑戰(zhàn)。

    一方面,目前對于直流微電網(wǎng)集群的關(guān)注不夠,僅有少量文獻(xiàn)研究其穩(wěn)定性,且重點(diǎn)針對小信號擾動,將系統(tǒng)的非線性問題簡化為線性時(shí)不變系統(tǒng)。文獻(xiàn)[12]中,基于直流微電網(wǎng)集群的多輸入多輸出線性狀態(tài)空間模型,推導(dǎo)出了局部Lyapunov函數(shù),保證了系統(tǒng)的漸近穩(wěn)定性。文獻(xiàn)[13]中,通過簡化各直流微電網(wǎng)的傳遞函數(shù)建立了直流微電網(wǎng)集群的小信號模型,通過根軌跡可以觀察到聯(lián)絡(luò)線的電感和CPL等參數(shù)對系統(tǒng)小信號穩(wěn)定性的影響。相同的方法也被應(yīng)用于交流微電網(wǎng)集群的小信號穩(wěn)定性評估與參數(shù)優(yōu)化設(shè)計(jì)[14]。然而,小信號分析只能預(yù)測平衡點(diǎn)附近一個(gè)小區(qū)域內(nèi)的系統(tǒng)動態(tài)行為,不適于大擾動分析。另一方面,與直流微電網(wǎng)及其系統(tǒng)相關(guān)的大信號穩(wěn)定性研究尚集中于單個(gè)直流微電網(wǎng),對于集群層面的大信號研究不足[15]。文獻(xiàn)[16]針對含有CPL的分布式電力系統(tǒng)的大信號穩(wěn)定性進(jìn)行了分析,并在Lyapunov穩(wěn)定性背景下,討論了混合勢理論穩(wěn)定性判據(jù)。此后,CPL引起的大信號穩(wěn)定性問題在復(fù)雜電力電子系統(tǒng)中得到了廣泛的關(guān)注和研究[17]。文獻(xiàn)[18]運(yùn)用混合勢函數(shù)、Takagi-Sugeno(T-S)模糊模型與逆軌跡等方法對比分析了含有CPL的直流微電網(wǎng)的大信號穩(wěn)定性。在文獻(xiàn)[7, 19-20]中,對直流微電網(wǎng)的大信號穩(wěn)定性分析中重點(diǎn)考慮了下垂控制的影響。文獻(xiàn)[20]同時(shí)運(yùn)用分岔圖輔助分析了直流微電網(wǎng)下垂控制在平衡點(diǎn)處的穩(wěn)定性??紤]到CPL的不確定性,文獻(xiàn)[21]研究了直流微電網(wǎng)在公共耦合點(diǎn)的大信號穩(wěn)定問題,并對電網(wǎng)阻抗進(jìn)行估計(jì)。文獻(xiàn)[22]中,以并網(wǎng)電壓源換流器為主電源的直流微電網(wǎng)為研究對象,建立了混合勢函數(shù)模型,并分別分析了恒阻抗、恒電流及恒功率負(fù)載所對應(yīng)的穩(wěn)定性判據(jù)。

    這些研究工作促進(jìn)了單個(gè)直流微電網(wǎng)大信號穩(wěn)定性的研究,但目前尚欠缺對于由多個(gè)直流微電網(wǎng)互聯(lián)形成的集群系統(tǒng)的建模和分析。僅有文獻(xiàn)[23]提出了基于T-S多模型法分析直流微電網(wǎng)集群的漸進(jìn)穩(wěn)定性區(qū)域(Region of Asymptotic Stability, RAS),但由于集群網(wǎng)絡(luò)的復(fù)雜性,對于系統(tǒng)大信號穩(wěn)定判定條件不易得到其全解析解,而采取部分解析、部分?jǐn)?shù)值的混合分析方式。

    為進(jìn)一步簡化直流微電網(wǎng)集群的大信號穩(wěn)定性解析,本文提出基于Brayton-Moser混合勢理論的分析方法。首先,針對直流微電網(wǎng)集群系統(tǒng),將其等效成降階模型,借助混合勢理論,直接構(gòu)造直流微電網(wǎng)集群系統(tǒng)等效降階模型的能量泛函,在Lyapunov理論框架下進(jìn)行系統(tǒng)的大信號穩(wěn)定性評估,為直流微電網(wǎng)集群的大信號穩(wěn)定性解析提供簡便且有效的理論工具。本文將詳細(xì)推導(dǎo)用于系統(tǒng)大信號穩(wěn)定性判據(jù)的混合勢函數(shù),分析關(guān)鍵參數(shù)對穩(wěn)定區(qū)間的影響,并給出實(shí)時(shí)仿真驗(yàn)證結(jié)果。

    1 直流微電網(wǎng)集群拓?fù)浼捌浞謱涌刂?/h2>

    1.1 直流微電網(wǎng)集群拓?fù)?/h3>

    根據(jù)具體需求,直流微電網(wǎng)集群的組網(wǎng)可存在多種形式。例如,在直流配電網(wǎng)環(huán)境下,可采用串聯(lián)、并聯(lián)以及混聯(lián)等結(jié)構(gòu)[24]。文獻(xiàn)[25]進(jìn)一步指出,并聯(lián)結(jié)構(gòu)也分為直接并聯(lián)、并聯(lián)接入母線以及分段母線并聯(lián)等方式。本文考慮集群中各直流微電網(wǎng)為同一電壓等級,即采用直接互聯(lián)的方式組網(wǎng)。

    圖1給出了研究對象的拓?fù)浣Y(jié)構(gòu)。圖1a所示為單個(gè)低壓直流微電網(wǎng)的電路拓?fù)?,由可再生能源發(fā)電系統(tǒng)、蓄電池儲能系統(tǒng)和恒功率負(fù)載構(gòu)成,各異構(gòu)單元均需通過DC-DC或DC-AC電力電子接口連接至直流母線。圖1b為采用多個(gè)直流微電網(wǎng)通過聯(lián)絡(luò)線互聯(lián)形成且具有環(huán)形拓?fù)浣Y(jié)構(gòu)的直流微電網(wǎng)集群,而環(huán)形拓?fù)浼耗軌驅(qū)崿F(xiàn)高可靠性和高靈活性的供電目標(biāo)。圖中,rLC分別為蓄電池儲能單元線路電阻、蓄電池儲能單元濾波器電感和濾波器電容;rt和Lt分別為第條網(wǎng)間聯(lián)絡(luò)線的等效電阻和等效電感,it為第條網(wǎng)間聯(lián)絡(luò)線上的電流。

    圖1 直流微電網(wǎng)集群拓?fù)浣Y(jié)構(gòu)

    1.2 直流微電網(wǎng)集群的分層控制策略

    分層控制策略廣泛應(yīng)用于公共大電網(wǎng)和交流微電網(wǎng)[26],隨著直流微電網(wǎng)技術(shù)的發(fā)展[27-28],直流微電網(wǎng)集群系統(tǒng)也逐漸引入分層控制策略以實(shí)現(xiàn)電壓電流的調(diào)節(jié)及網(wǎng)間的能量管理[29-32]。使得直流微電網(wǎng)集群系統(tǒng)能夠在不同的時(shí)間尺度上實(shí)現(xiàn)設(shè)備級控制與系統(tǒng)級解耦控制,保證了靈活的電壓電流調(diào)節(jié)和功率均衡,實(shí)現(xiàn)集群內(nèi)各直流微電網(wǎng)間區(qū)域能源共享和優(yōu)化利用。下面針對本文所研究的直流微電網(wǎng)集群系統(tǒng),以光伏-蓄電池儲能型(Photovoltaic- Battery Energy Storage System, PV-BESS)直流微電網(wǎng)為研究對象,首先說明單個(gè)微電網(wǎng)內(nèi)各單元的主電路拓?fù)浼捌淇刂平Y(jié)構(gòu)。

    蓄電池儲能系統(tǒng)(BESS)及其接口為組網(wǎng)(Grid- Forming, GF)單元,主要起到控制母線電壓的目的,其主電路拓?fù)浼捌淇刂撇呗匀鐖D2所示。圖2a給出了儲能單元的主電路為雙向Buck-Boost變換器和LC濾波器串聯(lián)的結(jié)構(gòu)。圖中,batik、oik、BLik和oik分別為第個(gè)DCMG內(nèi)第個(gè)蓄電池儲能單元變換器的輸入電壓、輸出電壓、電感電流和輸出電流;Bik、Bik、Bik和busi分別為蓄電池儲能單元變換器的電感寄生電阻、電感、輸出電容和第個(gè)DCMG的直流母線電壓。圖2b為儲能單元的分層控制結(jié)構(gòu)框圖。一次控制采用下垂控制,通過虛擬電阻來改變變換器的輸出特性,以實(shí)現(xiàn)各變換器間的功率分配;二次控制主要目的是補(bǔ)償由下垂控制所帶來的電壓跌落,一般通過一致性算法來實(shí)現(xiàn);而各直流微電網(wǎng)間的功率流動通常由三次控制所決定,以實(shí)現(xiàn)功率均衡和應(yīng)力分擔(dān)等目標(biāo)。圖中,oi和oj分別為本網(wǎng)內(nèi)各儲能單元輸出總電流和除本網(wǎng)外其他網(wǎng)內(nèi)各儲能單元輸出總電流;avg-i和avg-j分別為由一致性算法求出的本單元和相鄰單元的平均輸出電壓;dik、oref和vm分別為下垂電阻、參考電壓和參考電壓調(diào)整量。

    (a)蓄電池儲能單元主電路拓?fù)?/p>

    (b)分層控制策略

    圖2 蓄電池儲能系統(tǒng)主電路拓?fù)浼捌浞謱涌刂撇呗?/p>

    Fig.2 Topology and hierarchical control strategy of the BESS

    光伏發(fā)電單元的主電路拓?fù)浼捌淇刂瓶驁D如圖3所示,主電路為Boost變換器。為了提高發(fā)電效率,光伏發(fā)電單元通常工作在最大功率點(diǎn)追蹤(Maximum Power Point Tracking, MPPT)模式。當(dāng)光伏發(fā)電單元工作在MPPT模式時(shí),從直流母線側(cè)看,整個(gè)光伏發(fā)電單元可以等效成一個(gè)恒功率源(Constant Power Source, CPS),其功率值對應(yīng)光伏陣列最大功率點(diǎn)功率。圖中,pvi和pvi分別為第個(gè)DCMG內(nèi)光伏陣列的輸出電壓和輸出電流;Pi、Pi和Pi分別為光伏單元變換器的電感寄生電阻、電感和輸出電容。

    圖3 光伏單元主電路拓?fù)浼捌淇刂瓶驁D

    在實(shí)際應(yīng)用中,直流微電網(wǎng)的負(fù)載通常由一個(gè)嚴(yán)格控制的DC-DC或DC-AC變換器加上直流或者交流負(fù)載[15]。本文所考慮的恒功率負(fù)載主電路及其控制框圖如圖4所示,主電路為單向Buck變換器加電阻性負(fù)載,控制環(huán)路為傳統(tǒng)電壓電流雙閉環(huán)比例積分(Proportional Integral, PI)控制。圖中,CPLi、iLi和Li分別為第個(gè)DCMG內(nèi)CPL的輸入電流、電感電流和輸出電壓;Li、Li、Li和Li分別為CPL變換器的電感寄生電阻、電感、輸出電容和負(fù)載電阻。在一定的電壓、電流范圍內(nèi),輸入端口可以被看成是一個(gè)恒功率負(fù)載(Constant Power Load, CPL)。換言之,從直流母線側(cè)看,CPL又可被等效成一個(gè)受直流母線電壓控制的電流源,CPL輸入電流CPL與直流母線電壓bus呈非線性關(guān)系,可表示為

    式中,CPL為恒功率負(fù)載的功率;min、max和min、max分別為CPL輸入電流和電壓的最小、最大值。

    圖4 恒功率負(fù)載主電路拓?fù)浼捌淇刂瓶驁D

    圖5 恒功率負(fù)載特性曲線

    1.3 直流微電網(wǎng)集群的大信號降階模型

    通過圖1可知,由于單個(gè)直流微電網(wǎng)內(nèi)存在多個(gè)電力電子接口電路,且考慮了連接線路阻抗、恒功率負(fù)載等多種因素,因此整體集群系統(tǒng)為典型的高階、非線性復(fù)雜電路網(wǎng)絡(luò),對系統(tǒng)的大信號建模與分析帶來了挑戰(zhàn)。

    然而,通過研究發(fā)現(xiàn),在直流微電網(wǎng)的分層控制中,直流微電網(wǎng)集群的功率流動通常由系統(tǒng)的三次控制層所決定。相對于一次和二次控制,三次控制環(huán)路具有更低截止頻率的低通濾波特性,因此系統(tǒng)的高頻動態(tài)特性可以被忽略。這為通過低頻降階模型研究直流微電網(wǎng)集群的動態(tài)特性奠定了基礎(chǔ),能夠從一定程度上簡化分析。

    考慮直流微電網(wǎng)內(nèi)BESS由并聯(lián)的多個(gè)蓄電池儲能單元構(gòu)成,一般通過一次和二次控制來維持母線電壓的穩(wěn)定;而三次控制動態(tài)調(diào)整電壓參考值以實(shí)現(xiàn)網(wǎng)間功率流動,且電壓偏差被控制在一定的范圍內(nèi),其輸出特性保持一致。因此,可將BESS等效為一個(gè)受控電壓源。

    另一方面,由光伏發(fā)電單元的輸出特性及恒功率負(fù)載的輸入特性可知,導(dǎo)出其等效大信號模型如圖6所示。光伏單元作為CPS與CPL可以看成一個(gè)新的等效恒功率負(fù)載(Equivalent Constant Power Load, ECPL),其等效功率ECPLi可表示為

    式中,PCPLi為恒功率負(fù)載功率之和;PCPSi為光伏單元在MPPT模式下作為CPS所提供的功率之和。

    對于單個(gè)直流微電網(wǎng),目前已有多篇文獻(xiàn)采用二階簡化電路模型對系統(tǒng)穩(wěn)定性進(jìn)行了分析[7, 21, 23]?;谏鲜龅刃Х治觯瑘D7給出了本文所考慮的直流微電網(wǎng)等效降階電路模型,單個(gè)直流微電網(wǎng)被簡化為RLC電路網(wǎng)絡(luò)。如第個(gè)DCMG由等效受控電壓源V、線路電阻r、線路電感L、母線電容C和等效恒功率負(fù)載ECPLi組成。

    圖7 直流微電網(wǎng)等效降階電路模型

    根據(jù)電路參數(shù)取值區(qū)間,考慮設(shè)計(jì)條件

    式中,等效線路電阻r可推導(dǎo)為

    等效線路電感L

    等效母線電容C

    根據(jù)圖7所示直流微電網(wǎng)的等效電路模型,可導(dǎo)出集群網(wǎng)絡(luò)的動態(tài)模型,其狀態(tài)方程可描述為

    式中,i為DCMG內(nèi)輸出總電流。

    觀察式(7)發(fā)現(xiàn):①等效電路模型為3×階動態(tài)系統(tǒng);②由式(7)中直流母線電壓的微分方程可知,系統(tǒng)的非線性因素重點(diǎn)表現(xiàn)為各直流微電網(wǎng)內(nèi)等效恒功率負(fù)載的特性。運(yùn)用通用Lyapunov穩(wěn)定性理論分析非線性系統(tǒng)穩(wěn)定性,關(guān)鍵在于構(gòu)造一個(gè)合適的Lyapunov函數(shù)。由于非線性系統(tǒng)的復(fù)雜性,目前尚無構(gòu)造Lyapunov函數(shù)的統(tǒng)一方法,通常采用克拉索夫斯基法(Krasovsky Method, KM)和變量梯度法(Variable Gradient Method, VGM)構(gòu)造Lyapunov函數(shù)。對于直流微電網(wǎng)集群這樣一個(gè)高階非線性系統(tǒng),通過現(xiàn)有方法構(gòu)造Lyapunov函數(shù)對其進(jìn)行大信號穩(wěn)定性分析會帶來龐大的計(jì)算量;隨著集群微電網(wǎng)個(gè)數(shù)的增加,計(jì)算量將呈指數(shù)增長,維數(shù)災(zāi)難將最終導(dǎo)致函數(shù)無法求解的問題。

    2 直流微電網(wǎng)集群的混合勢函數(shù)構(gòu)造

    在基于Lyapunov函數(shù)的非線性系統(tǒng)穩(wěn)定性分析中,一般有兩種途徑可達(dá)到簡化分析的目的: ①得到有效的簡化模型以便于Lyapunov函數(shù)的構(gòu)造;②提出簡化的Lyapunov函數(shù)構(gòu)造形式。本文提出基于Brayton-Moser混合勢理論的直流微電網(wǎng)集群大信號穩(wěn)定性分析方法屬于兩者結(jié)合,首先對直流微電網(wǎng)集群系統(tǒng)進(jìn)行模型降階,而后基于降階模型利用混合勢理論對系統(tǒng)進(jìn)行大信號穩(wěn)定性分析。混合勢理論就是一種簡化的Lyapunov類型函數(shù)構(gòu)造方法,通過建立一個(gè)特殊形式的能量函數(shù)(勢函數(shù)),為非線性RLC電路網(wǎng)絡(luò)的動態(tài)特性分析提供了有效的解決途徑,近年來也被應(yīng)用于各類電力電子系統(tǒng)的穩(wěn)定性分析[17-19, 22, 34]。

    2.1 混合勢理論及其穩(wěn)定性判據(jù)

    根據(jù)Brayton-Moser混合勢理論的定義,混合勢函數(shù)可以在特定的條件下構(gòu)造Lyapunov型函數(shù)以分析和證明非線性電路網(wǎng)絡(luò)的穩(wěn)定性[35]。其中,“勢函數(shù)”代表函數(shù)項(xiàng)包含電壓和電流的乘積,“混合勢函數(shù)”進(jìn)一步表明其構(gòu)造直接來源于電感、電容和電阻等多個(gè)相關(guān)混合變量勢函數(shù)的疊加。同時(shí),混合勢理論也特別適于包含負(fù)阻抗元件的非線性電路網(wǎng)絡(luò)的穩(wěn)定性分析。

    假設(shè)對一個(gè)非線性電路網(wǎng)絡(luò),其動態(tài)特性可表述為

    式中,為電路中電感元件;為電路中電容元件。

    那么對于這個(gè)非線性電路網(wǎng)絡(luò),其混合勢函數(shù)可構(gòu)造為

    式中,=1,…,i為流經(jīng)電感的電流;=v+1,…,v+s為電容支路的電壓,其用于非線性電路網(wǎng)絡(luò)穩(wěn)定性分析標(biāo)量函數(shù)的一般形式可改寫為

    式中,()為電路中非儲能元件的電流勢函數(shù);()為電路中非儲能元件的電壓勢函數(shù);(,)為電路中電容的能量以及部分非儲能元件的能量,它由電路的拓?fù)浣Y(jié)構(gòu)決定,為常數(shù)矩陣。

    混合勢函數(shù)穩(wěn)定性定理[35]:對于穩(wěn)定性分析,可以構(gòu)造函數(shù)為

    其中

    如果對于電路中所有、,滿足條件

    2.2 直流微電網(wǎng)集群的混合勢函數(shù)推導(dǎo)

    根據(jù)1.3節(jié)建立的直流微電網(wǎng)集群等效降階電路模型,對其構(gòu)造混合勢函數(shù),可依據(jù)電路元件分為以下幾部分:

    (1)等效電壓源V的電流勢為

    (2)網(wǎng)內(nèi)線路電阻r的電流勢為

    (3)網(wǎng)間聯(lián)絡(luò)線電阻rt的電流勢為

    (4)等效恒功率負(fù)載ECPLi的電流勢為

    (5)網(wǎng)內(nèi)母線電容C的能量為

    綜合直流微電網(wǎng)集群等效降階模型的各部分電流勢和能量函數(shù),整體集群系統(tǒng)的混合勢函數(shù)可表示為

    根據(jù)式(19)對系統(tǒng)各狀態(tài)變量求偏導(dǎo)可得

    通過對比可知,式(20)滿足式(8),說明系統(tǒng)混合勢函數(shù)模型建立正確。

    2.3 基于混合勢理論的穩(wěn)定性判據(jù)

    由式(19),根據(jù)2.1節(jié)所介紹的Brayton-Moser混合勢理論可導(dǎo)出

    進(jìn)而,由穩(wěn)定性定理導(dǎo)出確保系統(tǒng)大信號穩(wěn)定的充分條件為

    3 系統(tǒng)大信號穩(wěn)定性分析

    3.1 大信號穩(wěn)定性判據(jù)

    本文以4個(gè)直流微電網(wǎng)通過聯(lián)絡(luò)線互聯(lián)形成具有環(huán)形拓?fù)浣Y(jié)構(gòu)的集群為具體研究對象,其結(jié)構(gòu)如圖1所示,其中每個(gè)直流微電網(wǎng)由一個(gè)光伏發(fā)電單元、一個(gè)蓄電池儲能系統(tǒng)(包含2個(gè)儲能單元)以及一個(gè)恒功率負(fù)載構(gòu)成,則=4,=2。

    因此,基于混合勢理論的穩(wěn)定性定理,得到確保系統(tǒng)大信號穩(wěn)定的充分條件為

    為便于分析,假設(shè)集群內(nèi)各微電網(wǎng)的參數(shù)保持一致,本文中主要變量參數(shù)的單位分別為電壓單位為V、電阻單位為W、電感單位為mH、電容單位為mF、功率單位為kW;光伏發(fā)電單元、蓄電池儲能單元和恒功率負(fù)載的主電路及控制環(huán)路參數(shù)分別見表1~表3,其中光伏發(fā)電單元采用增量電導(dǎo)法(Incre- mental Conductance Method, ICM)實(shí)現(xiàn)MPPT控制,在標(biāo)準(zhǔn)工況下(輻照強(qiáng)度r=1kW/m2,溫度=25℃),輸出最大功率約為1kW(即CPSi=1kW)。

    根據(jù)1.3節(jié)建立直流微電網(wǎng)集群系統(tǒng)等效降階電路模型的方法,針對上述具體研究對象,可以得出整個(gè)直流微電網(wǎng)集群系統(tǒng)等效降階模型和聯(lián)絡(luò)線阻抗的參數(shù)見表4。

    表1 光伏發(fā)電單元主電路及控制環(huán)路參數(shù)

    Tab.1 Circuit and control loop parameters of the photovoltaic power generation unit

    表2 蓄電池儲能系統(tǒng)單元主電路及控制環(huán)路參數(shù)

    Tab.2 Parameters of the power circuit and the control loop of the BESS unit

    表3 恒功率負(fù)載主電路及控制環(huán)路參數(shù)

    Tab.3 Circuit and control loop parameters of the CPL

    表4 直流微電網(wǎng)集群等效降階模型參數(shù)

    Tab.4 Equivalent reduced-order model parameters of the DCMG cluster

    3.2 拓?fù)渥兓瘜Ψ€(wěn)定性的影響

    由圖1給出的環(huán)形集群拓?fù)淇芍粢蚓W(wǎng)絡(luò)故障其中一條網(wǎng)間聯(lián)絡(luò)線發(fā)生斷路的情況,仍能夠保持各直流微電網(wǎng)互聯(lián),集群拓?fù)鋭t由環(huán)形演變?yōu)殒溞?,但此種拓?fù)涞淖兓瘜ο到y(tǒng)大信號穩(wěn)定性的影響尚未可知。

    針對本文具體研究對象,假設(shè)4號聯(lián)絡(luò)線斷路,系統(tǒng)動態(tài)從階數(shù)而言由12階變?yōu)?1階系統(tǒng)。此時(shí),根據(jù)混合勢理論導(dǎo)出滿足系統(tǒng)大信號穩(wěn)定性的充分條件為

    假設(shè)在集群內(nèi)各直流微電網(wǎng)參數(shù)及各網(wǎng)間聯(lián)絡(luò)線保持一致的情況下,對比式(23)與式(24)發(fā)現(xiàn),環(huán)形拓?fù)浜玩溞瓮負(fù)涞拇笮盘柗€(wěn)定性條件保持一致,均為

    式中,當(dāng)rt/Lt>r/L時(shí),則k=r/L,此時(shí)r/L為最小特征值;當(dāng)rt/Lt<r/L時(shí),則k=rt/Lt,此時(shí)rt/Lt為最小特征值,可分別考慮=或≠的情況。

    將各參數(shù)代入式(25)中的穩(wěn)定性判據(jù)可得

    需要指出,盡管基于混合勢函數(shù)得到的穩(wěn)定性條件為充分非必要條件,但仍舊為直流微電網(wǎng)集群的大信號穩(wěn)定性判定與分析提供了可供參考的依據(jù),同時(shí)有效降低了分析計(jì)算量。

    3.3 系統(tǒng)參數(shù)對穩(wěn)定性的影響

    根據(jù)穩(wěn)定性判據(jù)式(23)、式(24)可知,只需取穩(wěn)定性判據(jù)中兩部分最小值的和即可得到最終判據(jù),所以分析不同電路參數(shù)對系統(tǒng)穩(wěn)定性的影響時(shí),需區(qū)分兩種情況:①rt/Lt>r/L;②rt/Lt<r/L

    圖8 當(dāng)rit/Lit>ri/Li時(shí)的參數(shù)穩(wěn)定性邊界

    圖9 當(dāng)rit/Lit>ri/Li時(shí)的三維穩(wěn)定空間

    圖10 當(dāng)rit/Lit<ri/Li時(shí)的參數(shù)穩(wěn)定性邊界

    圖11 當(dāng)rit/Lit<ri/Li時(shí)的三維穩(wěn)定空間

    4 實(shí)時(shí)仿真證明

    為了驗(yàn)證上述直流微電網(wǎng)集群大信號穩(wěn)定性分析的正確性,搭建了實(shí)時(shí)仿真平臺,仿真模型參數(shù)與表1~表4所列一致。圖12給出了實(shí)時(shí)仿真測試平臺及實(shí)時(shí)仿真過程。如圖12所示,實(shí)時(shí)仿真過程主要包括三個(gè)步驟:①在Matlab/Simulink環(huán)境中搭建所研究直流微電網(wǎng)集群系統(tǒng)的仿真模型并初步進(jìn)行數(shù)值仿真證明;②將仿真模型編譯成代碼文件,然后通過上位機(jī)(dSPACE Control Desk)下載到dSPACE(DS1202)MicroLabBox硬件平臺上;③進(jìn)行實(shí)時(shí)仿真證明。

    4.1 分層控制策略有效性驗(yàn)證

    本節(jié)重點(diǎn)對直流微電網(wǎng)集群分層控制策略的有效性進(jìn)行驗(yàn)證。根據(jù)3.2節(jié)的穩(wěn)定性判據(jù)求出系統(tǒng)大信號穩(wěn)定的充分條件為ECPLi<0.921 6kW,由于光伏發(fā)電單元在MPPT模式下輸出最大功率為CPSi=1kW,因此,假設(shè)在光伏單元輸出功率不變的情況下,可得出恒功率負(fù)載的穩(wěn)定范圍為CPLi<1.921 6kW。圖13給出了集群系統(tǒng)各網(wǎng)內(nèi)恒功率負(fù)載超過功率界限后系統(tǒng)不穩(wěn)定運(yùn)行又跳變回穩(wěn)定界限內(nèi)的系統(tǒng)響應(yīng)波形。由圖13可知,系統(tǒng)首先運(yùn)行在初始工況:各網(wǎng)內(nèi)恒功率負(fù)載均為CPLi=1.2kW、光伏在MPPT模式下輸出最大功率CPSi=1kW;而后各網(wǎng)內(nèi)恒功率負(fù)載跳變至CPLi=2.3kW,超出穩(wěn)定界限,系統(tǒng)處于不穩(wěn)定運(yùn)行狀態(tài);最后,各網(wǎng)內(nèi)恒功率負(fù)載分別跳變至CPL1=1.2kW、CPL2=1.4kW、CPL3=1.6kW和CPL4=1.8kW,此時(shí)各網(wǎng)內(nèi)負(fù)荷均不相同,而系統(tǒng)在分層控制的作用下實(shí)現(xiàn)了網(wǎng)間的功率均衡,由圖13b可知,各直流微電網(wǎng)內(nèi)蓄電池儲能系統(tǒng)的出力保持一致。由此可知,系統(tǒng)在超過穩(wěn)定界限后再次跳變回穩(wěn)定界限內(nèi)依然能夠穩(wěn)定地實(shí)現(xiàn)分層控制策略。因此,系統(tǒng)在滿足大信號穩(wěn)定性判據(jù)的情況,且直流微電網(wǎng)集群在分層控制的作用下,有效地保證了靈活的電壓電流調(diào)節(jié)和功率流動控制,實(shí)現(xiàn)了區(qū)域能源共享和優(yōu)化利用。

    圖13 系統(tǒng)在分層控制策略下功率均衡響應(yīng)

    4.2 恒功率負(fù)載跳變驗(yàn)證

    由4.1節(jié)分析結(jié)果可知,恒功率負(fù)載的穩(wěn)定范圍為CPLi<1.921 6kW。圖14所示為系統(tǒng)在環(huán)形拓?fù)湎轮绷髂妇€電壓隨CPL跳變的響應(yīng)波形,圖14分別為4個(gè)CPL功率(CPLi,=1, 2, 3, 4)同時(shí)從1.2kW跳變到2.1kW、2.2kW和2.3kW直流母線電壓響應(yīng)波形。圖15所示為系統(tǒng)在鏈形拓?fù)湎轮绷髂妇€電壓隨CPL跳變的響應(yīng)波形,圖15同樣分別為4個(gè)CPL功率(CPLi,=1, 2, 3, 4)同時(shí)從1.2kW跳變到2.1kW、2.2kW和2.3kW直流母線電壓響應(yīng)波形。由圖14a、圖15a可知,當(dāng)CPL功率從1.2kW跳變到2.1kW時(shí),雖然已超出穩(wěn)定邊界,但系統(tǒng)依舊能夠維持母線電壓穩(wěn)定運(yùn)行。由圖14b和圖15b可知,當(dāng)CPL功率從1.2kW跳變到2.2kW時(shí),直流母線電壓開始小幅振蕩,系統(tǒng)運(yùn)行開始趨于不穩(wěn)定。如圖14c和圖15c所示,當(dāng)CPL功率從1.2kW跳變到2.3kW時(shí),直流母線電壓開始大幅振蕩,系統(tǒng)不穩(wěn)定運(yùn)行。

    另外,綜合比較圖14和圖15可知,直流微電網(wǎng)集群在環(huán)形拓?fù)浜玩溞瓮負(fù)湎碌拇笮盘柗€(wěn)定性基本一致,這與前文大信號穩(wěn)定性理論分析結(jié)果相吻合。結(jié)合大信號穩(wěn)定性判據(jù)及實(shí)時(shí)仿真結(jié)果可知,基于混合勢函數(shù)得到的大信號穩(wěn)定性判據(jù)雖為充分非必要條件,預(yù)測具有一定的保守性,但分析偏差仍在可接受范圍。

    圖14 環(huán)形拓?fù)湎轮绷髂妇€電壓隨CPL跳變的響應(yīng)

    圖15 鏈形拓?fù)湎轮绷髂妇€電壓隨CPL跳變的響應(yīng)

    4.3 光伏輸出功率變化驗(yàn)證

    根據(jù)1.3節(jié)等效恒功率負(fù)載建??芍?,當(dāng)光伏單元工作于MPPT模式,由于光照、溫度等因素而導(dǎo)致輸出功率發(fā)生變化時(shí),等效恒功率負(fù)載的功率也隨之變化。為此,本節(jié)主要驗(yàn)證系統(tǒng)隨光伏發(fā)電單元輸出功率變化的響應(yīng)情況。圖16所示為系統(tǒng)隨光伏發(fā)電單元輸出功率變化的響應(yīng),集群系統(tǒng)首先運(yùn)行在初始工況:各網(wǎng)內(nèi)恒功率負(fù)載均為CPLi= 1.2kW、光伏在MPPT模式下輸出最大功率CPSi= 1kW;然后,各網(wǎng)內(nèi)恒功率負(fù)載跳變?yōu)镃PLi=1.9kW,由于滿足大信號穩(wěn)定性條件,系統(tǒng)還能夠穩(wěn)定運(yùn)行;假設(shè)輻照強(qiáng)度從r=1kW/m2減弱到原來的一半r= 0.5kW/m2,此時(shí)各光伏單元出力約為0.46kW,各直流母線電壓開始劇烈振蕩,系統(tǒng)不穩(wěn)定運(yùn)行。根據(jù)3.2節(jié)的穩(wěn)定性判據(jù)求出系統(tǒng)大信號穩(wěn)定的充分條件ECPLi<0.9216kW,在網(wǎng)內(nèi)各恒功率負(fù)載功率為1.9kW的情況下可得光伏出力的穩(wěn)定性邊界為CPSi>0.978 4kW。由上述條件可知,在輻照強(qiáng)度減半后,光伏單元出力減弱,不滿足系統(tǒng)大信號穩(wěn)定性條件,因此無法保證系統(tǒng)穩(wěn)定。

    圖16 系統(tǒng)隨光伏發(fā)電單元輸出功率變化的響應(yīng)

    4.4 拓?fù)浣Y(jié)構(gòu)變化驗(yàn)證

    基于第3.2節(jié)的理論分析可知,直流微電網(wǎng)集群在環(huán)形拓?fù)浜玩溞瓮負(fù)浣Y(jié)構(gòu)下穩(wěn)定性判據(jù)保持一致。因此,可推論得到,環(huán)形拓?fù)浜玩溞瓮負(fù)浣Y(jié)構(gòu)間的變化對系統(tǒng)穩(wěn)定性無本質(zhì)影響。另外,根據(jù)4.1節(jié)驗(yàn)證結(jié)果,當(dāng)各直流微電網(wǎng)內(nèi)負(fù)荷相同時(shí),則各網(wǎng)間無功率流動,即網(wǎng)間聯(lián)絡(luò)線上電流為零,此時(shí)網(wǎng)間聯(lián)絡(luò)線是否連接則對系統(tǒng)運(yùn)行無影響。

    為了得到有效的驗(yàn)證分析結(jié)果,本節(jié)采用了各網(wǎng)內(nèi)負(fù)荷均不相同的情況(CPL1=1.2kW、CPL2= 1.4kW、CPL3=1.6kW和CPL4=1.8kW)對拓?fù)渥兓M(jìn)行了驗(yàn)證。圖17給出了直流微電網(wǎng)集群在環(huán)形拓?fù)淝袚Q到鏈形拓?fù)浜笤偾袚Q回環(huán)形拓?fù)涞南到y(tǒng)響應(yīng)波形。由圖可知,當(dāng)4號聯(lián)絡(luò)線斷路時(shí),系統(tǒng)從環(huán)形拓?fù)渥兓癁殒溞瓮負(fù)?,此時(shí)網(wǎng)間聯(lián)絡(luò)線上的電流進(jìn)行了重新分配。而當(dāng)4號聯(lián)絡(luò)線恢復(fù)時(shí),系統(tǒng)又從鏈形拓?fù)淝袚Q回環(huán)形拓?fù)?,?jīng)過短暫調(diào)整,系統(tǒng)又穩(wěn)定運(yùn)行在環(huán)形拓?fù)湎?。同時(shí),各微電網(wǎng)的直流母線電壓和蓄電池儲能系統(tǒng)的出力做了微調(diào)以適應(yīng)拓?fù)浣Y(jié)構(gòu)的變化,從而達(dá)到新的功率均衡狀態(tài)。因此可知,只要系統(tǒng)滿足大信號穩(wěn)定判據(jù),環(huán)形拓?fù)渑c鏈形拓?fù)涞那袚Q能夠保證直流微電網(wǎng)集群的穩(wěn)定運(yùn)行。

    圖17 拓?fù)浣Y(jié)構(gòu)變化下的系統(tǒng)響應(yīng)

    5 結(jié)論

    針對直流微電網(wǎng)集群,本文提出了基于混合勢理論的大信號穩(wěn)定性分析方法,通過分析與驗(yàn)證可得如下結(jié)論:

    1)混合勢理論適用于直流微電網(wǎng)集群系統(tǒng),為通過構(gòu)造能量函數(shù)分析此類復(fù)雜網(wǎng)絡(luò)的大信號穩(wěn)定性提供了簡便且有效的依據(jù)。

    2)導(dǎo)出的穩(wěn)定性判據(jù)能夠有效預(yù)測系統(tǒng)的大信號穩(wěn)定區(qū)間,其保守性在可接受范圍之內(nèi)。

    3)恒功率負(fù)載為影響集群系統(tǒng)穩(wěn)定性的主要因素;系統(tǒng)大信號穩(wěn)定性對聯(lián)絡(luò)線電感的敏感度相對較低,穩(wěn)定區(qū)域?qū)β?lián)絡(luò)線電感的相對變化量并不顯著。

    4)環(huán)形與鏈形拓?fù)淝袚Q并未對系統(tǒng)的大信號穩(wěn)定性產(chǎn)生本質(zhì)影響,鏈形拓?fù)淠軌虮WC系統(tǒng)大信號穩(wěn)定運(yùn)行;但從組網(wǎng)結(jié)構(gòu)而言,環(huán)形拓?fù)浔孺溞瓮負(fù)涠嘁粭l潮流通道,因此環(huán)形拓?fù)渚哂懈玫娜哂嘈院透叩目煽啃浴?/p>

    后續(xù)研究將重點(diǎn)針對直流微電網(wǎng)集群系統(tǒng)的實(shí)證,同時(shí)考慮分布式通信網(wǎng)絡(luò)對系統(tǒng)穩(wěn)定性的影響,進(jìn)一步深入評估與驗(yàn)證大信號穩(wěn)定性的分析結(jié)果,形成有效的系統(tǒng)控制與設(shè)計(jì)方法。

    [1] Dragi?evi? T, Lu Xiaonan, Vasquez J C, et al. DC microgrids-part I: a review of control strategies and stabilization techniques[J]. IEEE Transactions on Power Electronics, 2016, 31(7): 4876-4891.

    [2] Lasseter R H. Microgrids[C]//2002 IEEE Power Engineering Society Winter Meeting, New York, NY, USA, 2002: 305-308.

    [3] Ito Y, Yang Zhongqing, Akagi H. DC microgrid based distribution power generation system[C]//the 4th International Power Electronics and Motion Control Conference, Xi’an, China, 2004: 1740-1745.

    [4] Li Xialin, Guo Li, Li Yunwei, et al. Flexible interlinking and coordinated power control of multiple DC microgrids clusters[J]. IEEE Transa- ctions on Sustainable Energy, 2018, 9(2): 904-915.

    [5] Meng Lexuan, Shafiee Q, Trecate G F, et al. Review on control of DC microgrids and multiple microgrid clusters[J]. IEEE Journal of Emerging and Selected Topics in Power Electronics, 2017, 5(3): 928-948.

    [6] 張偉亮, 張輝, 支娜, 等. 考慮網(wǎng)絡(luò)損耗的基于模型預(yù)測直流微電網(wǎng)群能量優(yōu)化策略[J]. 電力系統(tǒng)自動化, 2021, 45(13): 49-56.

    Zhang Weiliang, Zhang Hui, Zhi Na, et al. Model prediction based energy optimization strategy for DC microgrid groups considering network loss[J]. Auto- mation of Electric Power Systems, 2021, 45(13): 49-56.

    [7] Xie Wenqiang, Han Minxiao, Cao Wenyuan, et al. System-level large-signal stability analysis of droop- controlled DC microgrids[J]. IEEE Transactions on Power Electronics, 2021, 36(4): 4224-4236.

    [8] Liu Sucheng, Zheng Jiazhu, Li Run, et al. Multiple Lyapunov function-based large signal stability analysis of DC microgrid with coordinated control[C]// 2019 22nd International Conference on Electrical Machines and Systems, Harbin, China, 2019: 1-6.

    [9] 朱曉榮, 李錚, 孟凡奇. 基于不同網(wǎng)架結(jié)構(gòu)的直流微電網(wǎng)穩(wěn)定性分析[J]. 電工技術(shù)學(xué)報(bào), 2021, 36(1): 166-178.

    Zhu Xiaorong, Li Zheng, Meng Fanqi. Stability analysis of DC microgrid based on different grid structures[J]. Transactions of China Electrotechnical Society, 2021, 36(1): 166-178.

    [10] 施靜容, 李勇, 賀悝, 等. 一種提升交直流混合微電網(wǎng)動態(tài)特性的綜合慣量控制方法[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(2): 337-345.

    Shi Jingrong, Li Yong, He Li, et al. A comprehensive inertia control method for improving the dynamic characteristics of hybrid AC-DC microgrid[J]. Transa- ctions of China Electrotechnical Society, 2020, 35(2): 337-345.

    [11] 張偉亮, 張輝, 支娜, 等. 環(huán)形直流微電網(wǎng)故障分析與保護(hù)[J]. 電力系統(tǒng)自動化, 2020, 44(24): 105- 110.

    Zhang Weiliang, Zhang Hui, Zhi Na, et al. Fault analysis and protection of ring DC microgrid[J]. Automation of Electric Power Systems, 2020, 44(24): 105-110.

    [12] Han Renke, Tucci M, Martinelli A, et al. Stability analysis of primary plug-and-play and secondary leader-based controllers for DC microgrid clusters[J]. IEEE Transactions on Power Systems, 2019, 34(3): 1780-1800.

    [13] Shafiee Q, Dragicevic T, Vasquez J C, et al. Modeling, stability analysis and active stabilization of multiple DC-microgrid clusters[C]//2014 IEEE international energy conference, Cavtat, Croatia, 2014: 1284-1290.

    [14] He Jinghan, Wu Xiaoyu, Wu Xiangyu, et al. Small- signal stability analysis and optimal parameters design of microgrid clusters[J]. IEEE Access, 2019, 7: 36896-36909.

    [15] Marx D, Magne P, Nahid-Mobarakeh B, et al. Large signal stability analysis tools in DC power systems with constant power loads and variable power loads-a review[J]. IEEE Transactions on Power Electronics, 2012, 27(4): 1773-1787.

    [16] Belkhayat M, Cooley R, Witulski A. Large signal stability criteria for distributed systems with constant power loads[C]//Proceedings of PESC'95-Power Elec- tronics Specialist Conference, Atlanta, GA, USA, 1995: 1333-1338.

    [17] Jiang Jianbo, Liu Fei, Pan Shangzhi, et al. A conservatism-free large signal stability analysis method for DC microgrid based on mixed potential theory[J]. IEEE Transactions on Power Electronics, 2019, 34(11): 11342-11351.

    [18] 滕昌鵬, 王玉斌, 周博愷, 等. 含恒功率負(fù)載的直流微網(wǎng)大信號穩(wěn)定性分析[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(5): 973-982.

    Teng Changpeng, Wang Yubin, Zhou Bokai, et al. Large-signal stability analysis of DC microgrid with constant power loads[J]. Transactions of China Elec- trotechnical Society, 2019, 34(5): 973-982.

    [19] Li Zekun, Pei Wei, Ye Hua, et al. Large signal stability analysis for DC microgrid under droop control based on mixed potential theory[J]. The Journal of Engineering, 2019, 2019(16): 1189-1193.

    [20] Liu Haiyuan, Guo Wenzhong, Cheng Dong, et al. Stability and bifurcation analysis of DC microgrid with multiple droop control sources and loads[J]. IEEE Transactions on Power Electronics, 2021, 36(2): 2361-2372.

    [21] Sanchez S, Molinas M. Large signal stability analysis at the common coupling point of a DC microgrid: a grid impedance estimation approach based on a recursive method[J]. IEEE Transactions on energy conversion, 2015, 30(1): 122-131.

    [22] 厲澤坤, 孔力, 裴瑋. 直流微電網(wǎng)大擾動穩(wěn)定判據(jù)及關(guān)鍵因素分析[J]. 高電壓技術(shù), 2019, 45(12): 3993-4002.

    Li Zekun, Kong Li, Pei Wei. Analyses of stability criterion and key factors of DC microgrid under large disturbance[J]. High Voltage Engineering, 2019, 45(12): 3993-4002.

    [23] Liu Sucheng, Li Xiang, Xia Mengyu, et al. Takagi- Sugeno multimodeling-based large signal stability analysis of DC microgrid clusters[J]. IEEE Transa- ctions on Power Electronics, 2021, 36(11): 12670- 12684.

    [24] 劉海濤, 熊雄, 季宇, 等. 直流配電下多微網(wǎng)系統(tǒng)集群控制研究[J]. 中國電機(jī)工程學(xué)報(bào), 2019, 39(24): 7159-7167, 7489.

    Liu Haitao, Xiong Xiong, Ji Yu, et al. Cluster control research of multi-microgrids system under DC distribution system[J]. Proceedings of the CSEE, 2019, 39(24): 7159-7167, 7489.

    [25] 顧后生, 李霞林, 郭力, 等. 多直流微電網(wǎng)群柔性互聯(lián)與控制[J]. 電力系統(tǒng)及其自動化學(xué)報(bào), 2020, 32(4): 1-8.

    Gu Housheng, Li Xialin, Guo Li, et al. Flexible interconnection and control of multiple DC microgrid clusters[J]. Proceedings of the CSU-EPSA, 2020, 32(4): 1-8.

    [26] Tan D. Structured microgrids (SμGs) and flexible electronic large power transformers (FeLPTs)[J]. CES Transactions on Electrical Machines and Systems, 2020, 4(4): 255-263.

    [27] 胡長斌, 王慧圣, 羅珊娜, 等. 計(jì)及直流微電網(wǎng)擾動抑制的殘差動態(tài)分散補(bǔ)償控制策略[J]. 電工技術(shù)學(xué)報(bào), 2021, 36(21): 4493-4507, 4543.

    Hu Changbin, Wang Huisheng, Luo Shanna, et al. Sesidual dynamic decentralized compensation control strategy considering disturbance suppression[J]. Transactions of China Electrotechnical Society, 2021, 36(21): 4493-4507, 4543.

    [28] 喻禮禮, 張兆云, 劉藝濤. 基于改進(jìn)自適應(yīng)下垂的直流微電網(wǎng)穩(wěn)定分析與研究[J]. 電氣技術(shù), 2020, 21(5): 28-32.

    Yu Lili, Zhang Zhaoyun, Liu Yitao. Stability analysis and research of DC microgrid based on improved adaptive drooping[J]. Electrical Engineering, 2020, 21(5): 28-32.

    [29] Shafiee Q, Dragi?evi? T, Vasquez J C, et al. Hierarchical control for multiple DC-microgrids clusters[J]. IEEE Transactions on Energy Conversion, 2014, 29(4): 922-933.

    [30] Liu Sucheng, Liu Rui, Zheng Jiazhu, et al. Predictive function control in tertiary level for power flow management of DC microgrid clusters[J]. Electronics Letters, 2020, 56(13): 675-676.

    [31] 李忠文, 程志平, 張書源, 等. 考慮經(jīng)濟(jì)調(diào)度及電壓恢復(fù)的直流微電網(wǎng)分布式二次控制[J]. 電工技術(shù)學(xué)報(bào), 2021, 36(21): 4482-4492.

    Li Zhongwen, Cheng Zhiping, Zhang Shuyuan, et al. Distributed secondary control for economic dispatch and voltage restoration of DC microgrid[J]. Transa- ctions of China Electrotechnical Society, 2021, 36(21): 4482-4492.

    [32] 楊美輝, 周念成, 王強(qiáng)鋼, 等. 基于分布式協(xié)同的雙極直流微電網(wǎng)不平衡電壓控制策略[J]. 電工技術(shù)學(xué)報(bào), 2021, 36(3): 634-645.

    Yang Meihui, Zhou Niancheng, Wang Qianggang, et al. Unbalanced voltage control strategy of bipolar DC microgrid based on distributed cooperation[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 634-645.

    [33] Emadi A, Khaligh A, Rivetta C H, et al. Constant power loads and negative impedance instability in automotive systems: definition, modeling, stability, and control of power electronic converters and motor drives[J]. IEEE Transactions on Vehicular Tech- nology, 2006, 55(4): 1112-1125.

    [34] Peng Dongdong, Huang Meng, Li Jinhua, et al. Large-signal stability criterion for parallel-connected DC-DC converters with current source equivalence[J]. IEEE Transactions on Circuits and Systems II: Express Briefs, 2019, 66(12): 2037-2041.

    [35] Brayton R K, Moser J K. A theory of nonlinear networks. I[J]. Quarterly of Applied Mathematics, 1964, 22(1): 1-33.

    Large Signal Stability Analysis for DC Microgrid Clusters

    1,21,21,21,21,2

    (1. School of Electrical and Information Engineering Anhui University of Technology Maanshan 243000 China 2. Key Lab of Power Electronics & Motion Control Anhui University of Technology Maanshan 243000 China)

    DC microgrid (DCMG) clusters are in general formed by interconnecting multiple DCMGs to achieve zonal energy sharing and optimized utilization through flexible power flow control, and thus the advantages of DC-based distributed generation systems can be fully exploited. However, small-scale DCMGs are weak grids of low inertia and high impedance, and hence the weak-weak interconnection will reduce the damping of DCMG clusters, and even lead to severe consequences like oscillation and system collapse. In the meantime, the dynamic characteristics of higher-order, strong coupling, and nonlinearity bring great challenges to the stability analysis of DCMG clusters. To address this issue, a method for large signal stability analysis of DCMG clusters was proposed based on Brayton-Moser mixed potential theory. The large signal reduced-order model of the DCMG cluster was built, and the mixed potential function that facilitates large signal stability criterion was derived in detail, and the influence of the critical parameters on the stability region was analyzed. The correctness of the analysis is verified by real-time simulation results.

    DC microgrid (DCMG) cluster, large signal stability, constant power load (CPL), Brayton-Moser mixed potential theory

    10.19595/j.cnki.1000-6753.tces.210977

    TM46; TM712

    國家自然科學(xué)基金資助項(xiàng)目(51407003)。

    2021-07-05

    2021-08-04

    劉宿城 男,1981年生,博士,副教授,研究方向?yàn)殡娏﹄娮酉到y(tǒng)建模與控制、直流微電網(wǎng)。E-mail: liusucheng@126.com(通信作者)

    李 響 男,1996年生,碩士研究生,研究方向?yàn)橹绷魑㈦娋W(wǎng)。E-mail: lixiang960120@gmail.com

    (編輯 陳 誠)

    猜你喜歡
    信號系統(tǒng)
    Smartflower POP 一體式光伏系統(tǒng)
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    WJ-700無人機(jī)系統(tǒng)
    ZC系列無人機(jī)遙感系統(tǒng)
    北京測繪(2020年12期)2020-12-29 01:33:58
    完形填空二則
    基于PowerPC+FPGA顯示系統(tǒng)
    半沸制皂系統(tǒng)(下)
    孩子停止長個(gè)的信號
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    基于LabVIEW的力加載信號采集與PID控制
    久久久午夜欧美精品| 97在线人人人人妻| 日韩大片免费观看网站| 成年美女黄网站色视频大全免费 | 99国产精品免费福利视频| 啦啦啦啦在线视频资源| 国产成人a∨麻豆精品| 男女下面进入的视频免费午夜| 国产av码专区亚洲av| 国产免费又黄又爽又色| 美女中出高潮动态图| 国产精品99久久99久久久不卡 | 又粗又硬又长又爽又黄的视频| 日韩欧美 国产精品| 人人妻人人看人人澡| 欧美日韩亚洲高清精品| 精品少妇久久久久久888优播| 国产精品成人在线| 国产精品国产av在线观看| 中文字幕久久专区| 免费黄网站久久成人精品| 午夜激情久久久久久久| 波野结衣二区三区在线| 午夜视频国产福利| 在线观看一区二区三区激情| 亚洲国产精品国产精品| 国产人妻一区二区三区在| 久久99热这里只频精品6学生| 国产高清国产精品国产三级 | 欧美成人精品欧美一级黄| 国产精品久久久久成人av| 亚洲av成人精品一二三区| 久久人妻熟女aⅴ| 少妇高潮的动态图| 亚洲最大成人中文| 大香蕉久久网| 午夜视频国产福利| 国产在线一区二区三区精| 日韩 亚洲 欧美在线| 搡老乐熟女国产| 97超视频在线观看视频| 亚洲熟女精品中文字幕| av在线观看视频网站免费| 只有这里有精品99| 精品酒店卫生间| 久久久久久久久久久免费av| 在线播放无遮挡| 老女人水多毛片| 91精品国产九色| 嫩草影院入口| 美女xxoo啪啪120秒动态图| 丰满迷人的少妇在线观看| 激情五月婷婷亚洲| xxx大片免费视频| 少妇被粗大猛烈的视频| 国产精品久久久久久久电影| 看免费成人av毛片| 免费大片黄手机在线观看| 大香蕉久久网| 免费人成在线观看视频色| 性色av一级| 18禁裸乳无遮挡免费网站照片| 最近中文字幕高清免费大全6| 麻豆国产97在线/欧美| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品,欧美精品| 精品熟女少妇av免费看| 观看免费一级毛片| 精品人妻一区二区三区麻豆| 天堂俺去俺来也www色官网| 亚洲婷婷狠狠爱综合网| 青春草亚洲视频在线观看| 丰满人妻一区二区三区视频av| 自拍偷自拍亚洲精品老妇| 日本-黄色视频高清免费观看| 国产探花极品一区二区| 免费人成在线观看视频色| 丝瓜视频免费看黄片| 国产 精品1| 婷婷色综合www| 高清不卡的av网站| 久久精品久久精品一区二区三区| 青青草视频在线视频观看| 天堂8中文在线网| 高清av免费在线| 又大又黄又爽视频免费| av网站免费在线观看视频| 日韩av不卡免费在线播放| 精品久久久久久久久亚洲| 亚洲真实伦在线观看| 中文天堂在线官网| 国产熟女欧美一区二区| 成年美女黄网站色视频大全免费 | 国产男女内射视频| 日本vs欧美在线观看视频 | 国产精品女同一区二区软件| 美女xxoo啪啪120秒动态图| 免费观看无遮挡的男女| 亚洲精品国产av成人精品| 日本欧美国产在线视频| 嫩草影院新地址| 一本—道久久a久久精品蜜桃钙片| av在线蜜桃| 亚洲久久久国产精品| 成人一区二区视频在线观看| 人妻 亚洲 视频| 十八禁网站网址无遮挡 | 最近最新中文字幕大全电影3| 亚洲av在线观看美女高潮| 亚洲成人一二三区av| 久久久久国产精品人妻一区二区| 国产精品一区www在线观看| 国产精品一二三区在线看| 女的被弄到高潮叫床怎么办| 国产无遮挡羞羞视频在线观看| 免费观看的影片在线观看| 女人十人毛片免费观看3o分钟| 免费观看无遮挡的男女| 18禁裸乳无遮挡免费网站照片| 麻豆精品久久久久久蜜桃| 免费在线观看成人毛片| 99久国产av精品国产电影| 午夜福利视频精品| 国产成人免费无遮挡视频| 国精品久久久久久国模美| 校园人妻丝袜中文字幕| 超碰av人人做人人爽久久| 深爱激情五月婷婷| 久久99热这里只频精品6学生| 成人特级av手机在线观看| 激情五月婷婷亚洲| 欧美日韩视频高清一区二区三区二| 一本—道久久a久久精品蜜桃钙片| 菩萨蛮人人尽说江南好唐韦庄| 欧美bdsm另类| 亚洲精品日韩av片在线观看| av在线观看视频网站免费| 这个男人来自地球电影免费观看 | 国产精品久久久久久久电影| 成人高潮视频无遮挡免费网站| 免费看不卡的av| 亚洲性久久影院| 国产白丝娇喘喷水9色精品| 男女无遮挡免费网站观看| 丝瓜视频免费看黄片| 女的被弄到高潮叫床怎么办| 亚洲精品一区蜜桃| 亚洲三级黄色毛片| 色5月婷婷丁香| 少妇高潮的动态图| 国产精品一区二区三区四区免费观看| 黄色一级大片看看| 最近中文字幕2019免费版| 三级国产精品欧美在线观看| 国产男女超爽视频在线观看| 十八禁网站网址无遮挡 | 亚洲丝袜综合中文字幕| 黄色日韩在线| 日韩av在线免费看完整版不卡| 男女下面进入的视频免费午夜| 国产成人精品婷婷| 亚洲内射少妇av| 一级毛片电影观看| 激情 狠狠 欧美| 午夜老司机福利剧场| 欧美日韩国产mv在线观看视频 | 国产精品一区二区三区四区免费观看| 三级国产精品欧美在线观看| 亚洲成人中文字幕在线播放| 精品人妻视频免费看| 亚洲精品乱久久久久久| 狂野欧美激情性xxxx在线观看| 久久ye,这里只有精品| 综合色丁香网| 一本一本综合久久| 特大巨黑吊av在线直播| 国产伦精品一区二区三区四那| 亚洲性久久影院| 欧美激情极品国产一区二区三区 | 内射极品少妇av片p| 最近中文字幕2019免费版| h日本视频在线播放| 免费大片黄手机在线观看| 免费av不卡在线播放| 欧美xxⅹ黑人| freevideosex欧美| 精品国产乱码久久久久久小说| 一级毛片 在线播放| 亚洲成人一二三区av| 国产91av在线免费观看| 国产日韩欧美在线精品| 99久久精品一区二区三区| av福利片在线观看| 亚洲成色77777| 精品国产一区二区三区久久久樱花 | 久久午夜福利片| 久久午夜福利片| 在现免费观看毛片| 最近2019中文字幕mv第一页| 韩国高清视频一区二区三区| 男人狂女人下面高潮的视频| 熟妇人妻不卡中文字幕| 极品少妇高潮喷水抽搐| 国产成人a区在线观看| 小蜜桃在线观看免费完整版高清| 性高湖久久久久久久久免费观看| 黄色日韩在线| 男女边摸边吃奶| 日韩大片免费观看网站| 日日摸夜夜添夜夜爱| 性色avwww在线观看| 亚洲国产最新在线播放| 国产老妇伦熟女老妇高清| 日本黄色片子视频| 人人妻人人看人人澡| 边亲边吃奶的免费视频| 久久久久久伊人网av| 久久99热这里只频精品6学生| 欧美xxxx黑人xx丫x性爽| 大香蕉97超碰在线| 99热这里只有是精品在线观看| 久久久久久久亚洲中文字幕| 黄色一级大片看看| 亚洲电影在线观看av| 午夜福利网站1000一区二区三区| 亚洲欧美一区二区三区黑人 | 亚洲国产最新在线播放| 最近2019中文字幕mv第一页| 人妻系列 视频| 国产乱人偷精品视频| 亚洲欧美成人精品一区二区| 肉色欧美久久久久久久蜜桃| 我的老师免费观看完整版| 国产精品无大码| 97超视频在线观看视频| 国产黄片视频在线免费观看| 亚洲精品日本国产第一区| 熟妇人妻不卡中文字幕| 国产欧美亚洲国产| 免费观看av网站的网址| 免费高清在线观看视频在线观看| 国产一区亚洲一区在线观看| 亚洲色图综合在线观看| av女优亚洲男人天堂| 91精品国产九色| 色5月婷婷丁香| 亚洲一区二区三区欧美精品| 日本猛色少妇xxxxx猛交久久| 丰满少妇做爰视频| 人妻系列 视频| 久久精品国产亚洲av天美| 蜜桃久久精品国产亚洲av| 国产成人一区二区在线| 欧美亚洲 丝袜 人妻 在线| 欧美国产精品一级二级三级 | 最新中文字幕久久久久| 久久久国产一区二区| 最近中文字幕高清免费大全6| 女性被躁到高潮视频| 晚上一个人看的免费电影| 久热这里只有精品99| 久久国产精品大桥未久av | 夫妻性生交免费视频一级片| 中文天堂在线官网| 久久99热这里只有精品18| 亚洲欧美日韩卡通动漫| 成年人午夜在线观看视频| 久久久久国产精品人妻一区二区| 亚洲美女搞黄在线观看| 精品一品国产午夜福利视频| 亚洲中文av在线| 永久网站在线| 日本午夜av视频| 色哟哟·www| 亚洲高清免费不卡视频| 狂野欧美白嫩少妇大欣赏| 久久久久久人妻| 国产精品人妻久久久久久| 九九在线视频观看精品| 国产亚洲5aaaaa淫片| 免费不卡的大黄色大毛片视频在线观看| 成人漫画全彩无遮挡| 热99国产精品久久久久久7| 亚洲精品中文字幕在线视频 | 高清不卡的av网站| 99热全是精品| 精品亚洲乱码少妇综合久久| 国产伦在线观看视频一区| 日韩欧美 国产精品| 中国三级夫妇交换| 久久国内精品自在自线图片| 丝袜脚勾引网站| 性色avwww在线观看| 狂野欧美激情性bbbbbb| 亚洲第一av免费看| 国产免费福利视频在线观看| 97超视频在线观看视频| 一本一本综合久久| 国产一区有黄有色的免费视频| 人人妻人人爽人人添夜夜欢视频 | 免费少妇av软件| 你懂的网址亚洲精品在线观看| 国产成人免费观看mmmm| 色视频在线一区二区三区| 下体分泌物呈黄色| 在线亚洲精品国产二区图片欧美 | 五月玫瑰六月丁香| 99久国产av精品国产电影| 亚洲第一av免费看| 这个男人来自地球电影免费观看 | 在线观看三级黄色| 一本久久精品| 亚洲美女搞黄在线观看| 少妇高潮的动态图| 亚洲av.av天堂| 亚洲精品国产av成人精品| 国产久久久一区二区三区| 国产日韩欧美亚洲二区| 欧美97在线视频| 国产精品久久久久久久电影| 只有这里有精品99| 视频区图区小说| 2018国产大陆天天弄谢| 视频中文字幕在线观看| 亚洲av不卡在线观看| 女性生殖器流出的白浆| 国产成人免费观看mmmm| 国内精品宾馆在线| 尤物成人国产欧美一区二区三区| 日本黄色日本黄色录像| 中文字幕av成人在线电影| 国产成人精品婷婷| 欧美精品亚洲一区二区| 国产精品国产av在线观看| 深夜a级毛片| 我要看日韩黄色一级片| 国产老妇伦熟女老妇高清| 丝袜喷水一区| 在线天堂最新版资源| 国产精品av视频在线免费观看| 亚洲欧美日韩卡通动漫| 国产v大片淫在线免费观看| 久久av网站| 精品久久久久久电影网| 久久久色成人| 小蜜桃在线观看免费完整版高清| 高清黄色对白视频在线免费看 | 99视频精品全部免费 在线| 精品久久久久久久久亚洲| 天天躁日日操中文字幕| 乱码一卡2卡4卡精品| 九九在线视频观看精品| 搡女人真爽免费视频火全软件| 黄色怎么调成土黄色| 国产乱来视频区| 欧美高清成人免费视频www| 午夜福利高清视频| 国产国拍精品亚洲av在线观看| 国模一区二区三区四区视频| 九九爱精品视频在线观看| www.av在线官网国产| 男女下面进入的视频免费午夜| 好男人视频免费观看在线| 少妇精品久久久久久久| 亚洲激情五月婷婷啪啪| 校园人妻丝袜中文字幕| 国产精品久久久久久精品电影小说 | 国产 一区 欧美 日韩| 永久网站在线| 国产精品精品国产色婷婷| 国产v大片淫在线免费观看| 全区人妻精品视频| 国产亚洲欧美精品永久| 黄色视频在线播放观看不卡| 少妇的逼水好多| 免费高清在线观看视频在线观看| 一区二区三区四区激情视频| 涩涩av久久男人的天堂| av网站免费在线观看视频| 国产高清国产精品国产三级 | 网址你懂的国产日韩在线| 街头女战士在线观看网站| 国产高清国产精品国产三级 | 亚洲精品乱码久久久久久按摩| 人人妻人人添人人爽欧美一区卜 | 国产成人精品福利久久| 精品国产三级普通话版| 国产免费一区二区三区四区乱码| 国产日韩欧美亚洲二区| 成人毛片60女人毛片免费| 免费av中文字幕在线| 国产精品一二三区在线看| 成人高潮视频无遮挡免费网站| 亚洲国产av新网站| 国产人妻一区二区三区在| 十八禁网站网址无遮挡 | 久久久久久久久大av| 亚洲综合精品二区| 精品一品国产午夜福利视频| 免费观看a级毛片全部| 大香蕉97超碰在线| 久久久久久久久久成人| 五月伊人婷婷丁香| 亚洲高清免费不卡视频| 成年女人在线观看亚洲视频| 精品久久久久久久久亚洲| 毛片女人毛片| 久久鲁丝午夜福利片| 大码成人一级视频| 国产精品伦人一区二区| 人妻制服诱惑在线中文字幕| 街头女战士在线观看网站| 国产伦在线观看视频一区| 亚洲欧美精品专区久久| 能在线免费看毛片的网站| 欧美一级a爱片免费观看看| 亚洲精品,欧美精品| 一本久久精品| 亚洲无线观看免费| 五月伊人婷婷丁香| 一级毛片我不卡| 久久这里有精品视频免费| 国产精品无大码| 欧美成人a在线观看| 一级爰片在线观看| 亚洲真实伦在线观看| 国产精品久久久久久av不卡| 亚洲国产毛片av蜜桃av| 亚洲精品日韩在线中文字幕| 亚洲av电影在线观看一区二区三区| 建设人人有责人人尽责人人享有的 | 在线免费观看不下载黄p国产| 久久毛片免费看一区二区三区| 成人国产av品久久久| 在线观看免费视频网站a站| 日产精品乱码卡一卡2卡三| 久久6这里有精品| 又黄又爽又刺激的免费视频.| 中文字幕av成人在线电影| 天美传媒精品一区二区| 免费黄网站久久成人精品| 亚洲欧美精品自产自拍| 搡老乐熟女国产| 欧美高清性xxxxhd video| 日本黄大片高清| 青春草亚洲视频在线观看| 成年av动漫网址| freevideosex欧美| a级毛片免费高清观看在线播放| 国产精品嫩草影院av在线观看| 久久久欧美国产精品| 日本wwww免费看| 麻豆精品久久久久久蜜桃| 日本-黄色视频高清免费观看| 久久久久久久精品精品| 97超碰精品成人国产| 一级片'在线观看视频| 丝瓜视频免费看黄片| 国产一区二区三区av在线| 最近的中文字幕免费完整| 黄色日韩在线| 免费观看性生交大片5| 国产探花极品一区二区| 男女国产视频网站| 精品视频人人做人人爽| 亚洲av福利一区| 国产深夜福利视频在线观看| 全区人妻精品视频| 国产黄频视频在线观看| 国产黄色视频一区二区在线观看| 身体一侧抽搐| 欧美日韩国产mv在线观看视频 | 五月伊人婷婷丁香| 在线 av 中文字幕| 人妻制服诱惑在线中文字幕| 亚洲国产精品一区三区| 亚洲丝袜综合中文字幕| 搡女人真爽免费视频火全软件| 99热全是精品| 亚洲av日韩在线播放| 王馨瑶露胸无遮挡在线观看| 黄色配什么色好看| 18禁在线播放成人免费| 80岁老熟妇乱子伦牲交| 国产亚洲欧美精品永久| 97在线人人人人妻| av免费观看日本| 成年免费大片在线观看| 久久精品人妻少妇| 亚洲av免费高清在线观看| 91精品一卡2卡3卡4卡| 午夜激情久久久久久久| 日日啪夜夜爽| 日本wwww免费看| 蜜桃久久精品国产亚洲av| 免费观看a级毛片全部| 国产69精品久久久久777片| 一区二区三区精品91| 国产乱来视频区| 亚洲精品,欧美精品| 韩国av在线不卡| 亚洲欧美一区二区三区黑人 | 国产 一区精品| 精品人妻熟女av久视频| 亚洲成人av在线免费| 黄色怎么调成土黄色| 精品视频人人做人人爽| 亚洲国产av新网站| 国产无遮挡羞羞视频在线观看| 日韩欧美 国产精品| 日韩欧美精品免费久久| 国产探花极品一区二区| 欧美少妇被猛烈插入视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲熟女精品中文字幕| 极品教师在线视频| 最新中文字幕久久久久| 欧美成人午夜免费资源| 久久99精品国语久久久| 久久久久网色| 高清黄色对白视频在线免费看 | 99久久精品国产国产毛片| 黄色一级大片看看| 天天躁日日操中文字幕| 最近中文字幕2019免费版| 毛片一级片免费看久久久久| 黄色配什么色好看| av在线播放精品| 国产精品不卡视频一区二区| 欧美另类一区| 在线观看一区二区三区激情| 日本wwww免费看| 成人二区视频| 久久99蜜桃精品久久| 狂野欧美白嫩少妇大欣赏| 简卡轻食公司| 美女内射精品一级片tv| 久久久精品免费免费高清| 亚洲一区二区三区欧美精品| 能在线免费看毛片的网站| 九草在线视频观看| 久久人人爽av亚洲精品天堂 | 全区人妻精品视频| 色婷婷久久久亚洲欧美| 身体一侧抽搐| 国产老妇伦熟女老妇高清| 人妻夜夜爽99麻豆av| 国产精品熟女久久久久浪| 精品熟女少妇av免费看| 日韩免费高清中文字幕av| 最近中文字幕2019免费版| 毛片女人毛片| 国内少妇人妻偷人精品xxx网站| 美女cb高潮喷水在线观看| 免费av不卡在线播放| 欧美成人一区二区免费高清观看| 国产成人精品福利久久| 大片免费播放器 马上看| 最近中文字幕2019免费版| 亚洲欧美清纯卡通| 国产av码专区亚洲av| 亚洲国产最新在线播放| 99视频精品全部免费 在线| 插阴视频在线观看视频| 日本-黄色视频高清免费观看| 丝袜脚勾引网站| 人妻夜夜爽99麻豆av| 亚洲国产色片| 亚洲精品亚洲一区二区| a级毛色黄片| 久久99蜜桃精品久久| 18禁动态无遮挡网站| 免费播放大片免费观看视频在线观看| 在线亚洲精品国产二区图片欧美 | 极品教师在线视频| 看非洲黑人一级黄片| 久久久久久久国产电影| 国产伦精品一区二区三区视频9| 高清在线视频一区二区三区| 日韩中文字幕视频在线看片 | 麻豆国产97在线/欧美| 成人午夜精彩视频在线观看| av一本久久久久| 尾随美女入室| 一级二级三级毛片免费看| 18+在线观看网站| 青春草亚洲视频在线观看| 亚洲精品一二三| 亚洲欧洲日产国产| 久久久欧美国产精品| 人人妻人人添人人爽欧美一区卜 | 久久热精品热| 最近的中文字幕免费完整| av在线观看视频网站免费| 嘟嘟电影网在线观看| 国产精品秋霞免费鲁丝片| 在线观看人妻少妇| 国产亚洲欧美精品永久| 亚洲欧美成人精品一区二区| av在线播放精品| av福利片在线观看| 成年人午夜在线观看视频| 亚洲av国产av综合av卡| 干丝袜人妻中文字幕| 黄片wwwwww| 国产日韩欧美亚洲二区| 亚洲国产欧美在线一区| 在线观看一区二区三区激情| 日韩av免费高清视频| 亚洲精品一二三| 午夜免费鲁丝| 国产男女超爽视频在线观看| 亚洲婷婷狠狠爱综合网| 日本-黄色视频高清免费观看| 国产精品久久久久久精品电影小说 |