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

    基于增強(qiáng)型鯨魚優(yōu)化算法的交直流互聯(lián)配電網(wǎng)經(jīng)濟(jì)調(diào)度

    2024-09-21 00:00:00葉方慧熊國江楊曌昕徐波
    現(xiàn)代電子技術(shù) 2024年18期

    摘 "要: 隨著新能源的不斷發(fā)展和電力電子技術(shù)的廣泛應(yīng)用,交直流互聯(lián)配電網(wǎng)技術(shù)成為未來配電網(wǎng)的發(fā)展方向。合理有效的交直流互聯(lián)配電網(wǎng)調(diào)度策略可以提高新能源消納比例,降低運(yùn)行成本。文中建立了一種以設(shè)備運(yùn)行維護(hù)成本、儲能系統(tǒng)運(yùn)行成本、交直流雙向換流器損耗成本以及配電網(wǎng)與大電網(wǎng)之間購售電成本最小的交直流互聯(lián)配電網(wǎng)經(jīng)濟(jì)調(diào)度模型。在滿足配電網(wǎng)系統(tǒng)安全運(yùn)行的約束條件下,為了有效求解該模型,設(shè)計(jì)一種增強(qiáng)型鯨魚優(yōu)化算法,對限定對數(shù)螺旋形狀的常數(shù)b值進(jìn)行動態(tài)更新,提高了算法的精度并加快了其全局收斂。最后進(jìn)行仿真計(jì)算,驗(yàn)證所構(gòu)建模型的正確性以及算法的優(yōu)勢。結(jié)果表明,所提算法的收斂速度更快,全局收斂能力更強(qiáng),并且有效降低了系統(tǒng)的運(yùn)行經(jīng)濟(jì)成本。

    關(guān)鍵詞: 交直流互聯(lián)配電網(wǎng); 經(jīng)濟(jì)調(diào)度; 鯨魚優(yōu)化算法; 新能源消納; 儲能; 電壓源型換流器; 光伏發(fā)電系統(tǒng)

    中圖分類號: TN322.8+?34; TM73 " " " " " " " " " 文獻(xiàn)標(biāo)識碼: A " " " " " " " " " 文章編號: 1004?373X(2024)18?0015?07

    AC/DC interconnected distribution network economic dispatch based on EWOA

    YE Fanghui, XIONG Guojiang, YANG Zhaoxin, XU Bo

    (College of Electrical Engineering, Guizhou University, Guiyang 550025, China)

    Abstract: With the continuous development of new energy and the wide application of power electronic technology, AC/DC interconnected distribution network technology has become the development direction of future distribution network. The reasonable and effective dispatch strategy of AC/DC interconnected distribution network can improve the proportion of new energy consumption and reduce the operation cost. An economic dispatch model of AC/DC interconnected distribution network is established, which takes the cost of equipment operation and maintenance, the operating cost of energy storage system, the cost of AC/DC bi?directional converter loss, and the power purchasing and selling cost between the distribution network and the power grid as the minimum. In order to effectively solve the model under the constraint of ensuring the safe operation of the distribution network system, an enhanced whale optimization algorithm (EWOA) is designed to dynamically update the constant b value of the limited logarithmic spiral shape, which improves the accuracy of the algorithm and accelerates its global convergence. The simulation calculation is conducted to verify the correctness of the constructed model and the advantages of the algorithm. The results show that the proposed algorithm has a faster convergence speed, stronger global convergence ability, and effectively reduces the economic cost of system operation.

    Keywords: AC/DC interconnected distribution network; economic dispatch; whale optimization algorithm; new energy consumption; energy storage; voltage source converter; photovoltaic power generation system

    0 "引 "言

    近年來,隨著新能源和儲能技術(shù)[1]的不斷發(fā)展,用戶對直流負(fù)荷的需求日益增長,傳統(tǒng)交流配電網(wǎng)中接入直流電源及負(fù)荷的種類也逐漸增多,使得交流配電網(wǎng)中整流逆變環(huán)節(jié)增多,不僅影響了配電網(wǎng)的電能質(zhì)量,而且降低了配電網(wǎng)運(yùn)行的經(jīng)濟(jì)性[2?3]。交直流互聯(lián)配電網(wǎng)通過電壓源型換流器(Voltage Source Converter, VSC)[4]將交流側(cè)與直流側(cè)分隔,并將多個(gè)交直流轉(zhuǎn)換環(huán)節(jié)組合起來。交直流互聯(lián)配電網(wǎng)比傳統(tǒng)交流配電網(wǎng)傳輸容量更大,運(yùn)行更靈活,更有利于分布式光伏接入[5]。

    設(shè)計(jì)合理的經(jīng)濟(jì)調(diào)度方案[6]可以減小配電網(wǎng)的運(yùn)行成本,是合理高效利用可再生能源并提高系統(tǒng)運(yùn)行可靠性的基礎(chǔ)。文獻(xiàn)[7]中建立一個(gè)交直流互聯(lián)配電網(wǎng)的穩(wěn)態(tài)安全域模型,采用JAYA算法對安全域邊界點(diǎn)進(jìn)行求解。在交直流配電臺區(qū),文獻(xiàn)[8]中設(shè)計(jì)一種考慮源荷不確定性的日前調(diào)度方案,平抑系統(tǒng)的電壓波動,減少了網(wǎng)絡(luò)損耗。文獻(xiàn)[9]中提出一個(gè)基于動態(tài)網(wǎng)絡(luò)重構(gòu)以及多能源協(xié)同的交直流配電網(wǎng)優(yōu)化調(diào)度方案,但沒有考慮換流站的損耗成本。在考慮電價(jià)型需求響應(yīng)的情況下,文獻(xiàn)[10]中提出用求解器求解交直流混合配電網(wǎng)的優(yōu)化調(diào)度策略,在考慮電價(jià)型需求響應(yīng)的交直流混合配電系統(tǒng)運(yùn)行中,建立交直流混合配電網(wǎng)的優(yōu)化調(diào)度模型??紤]調(diào)節(jié)范圍、速度和時(shí)間的靈活性,文獻(xiàn)[11?12]中提出了一個(gè)交直流配電網(wǎng)分層分區(qū)優(yōu)化調(diào)度的實(shí)現(xiàn)方法,但未考慮儲能裝置的接入及其損耗成本。文獻(xiàn)[13]中提出了考慮多種分布式電源的柔性互聯(lián)配電系統(tǒng)多區(qū)域、多目標(biāo)分布式優(yōu)化模型,但未考慮柔性互聯(lián)開關(guān)的傳輸功率損耗。文獻(xiàn)[14]中提出了一個(gè)基于源荷不確定性的交直流混合主動配電網(wǎng)的分層?分布式優(yōu)化調(diào)度方案,有效提高了配電網(wǎng)的系統(tǒng)可控性以及對分布式能源的消納能力,但沒有考慮交直流雙向換流器的功率損耗。在配電網(wǎng)經(jīng)濟(jì)調(diào)度中較常用的算法有傳統(tǒng)優(yōu)化方法和智能優(yōu)化方法。鯨魚優(yōu)化算法(Whale Optimization Algorithm, WOA)[15]是2016 年由Seyedali Mirjalili教授提出的,它是模仿自然界中鯨魚捕食行為的新型群體智能優(yōu)化算法。與其他智能算法相比,該算法具有原理簡單、設(shè)置參數(shù)少、尋優(yōu)能力強(qiáng)的優(yōu)點(diǎn)。

    基于上述分析,本文建立一個(gè)以設(shè)備運(yùn)行維護(hù)成本、儲能系統(tǒng)運(yùn)行成本、交直流雙向換流器損耗成本以及配電網(wǎng)與大電網(wǎng)之間購售電成本之和最小的交直流互聯(lián)配電網(wǎng)經(jīng)濟(jì)調(diào)度模型,設(shè)計(jì)一種增強(qiáng)型鯨魚優(yōu)化算法(Enhanced Whale Optimization Algorithm, EWOA)對調(diào)度模型進(jìn)行求解。EWOA通過引入自適應(yīng)螺旋參數(shù),將對數(shù)螺旋形狀的常數(shù)b值設(shè)置成動態(tài)值,進(jìn)一步提升算法在種群內(nèi)的搜索能力,使算法的全局收斂速度加快,精度得到提高。以某交直流混合配電網(wǎng)示范工程為例,將所提方法和教與學(xué)優(yōu)化算法(Teaching?Learning?Based Optimization, TLBO)進(jìn)行對比分析,結(jié)果表明所提方案可以獲得更優(yōu)的調(diào)度策略。

    1 "交直流互聯(lián)配電網(wǎng)結(jié)構(gòu)

    交直流互聯(lián)配電網(wǎng)運(yùn)行方式有與大電網(wǎng)并網(wǎng)和離網(wǎng)兩種,本文采用并網(wǎng)運(yùn)行方式。交直流互聯(lián)配電網(wǎng)系統(tǒng)中,交流分布式電源(Distributed Generator, DG)和交流負(fù)荷接入交流母線,直流DG和直流負(fù)荷接入直流母線,交流子系統(tǒng)和直流子系統(tǒng)通過雙向AC/DC換流器連接實(shí)現(xiàn)功率平衡,其結(jié)構(gòu)圖如圖1所示。

    當(dāng)配網(wǎng)正常運(yùn)行時(shí),交流側(cè)和直流側(cè)的電源出力能滿足各側(cè)的負(fù)荷需求。主要運(yùn)行方式有以下4種。

    1) 交流側(cè)和直流側(cè)的電源出力相應(yīng)滿足各系統(tǒng)的負(fù)荷需求。

    2) 當(dāng)交流子系統(tǒng)或直流子系統(tǒng)存在功率不平衡時(shí),可通過VSC進(jìn)行功率交互,交直流兩側(cè)功率互為補(bǔ)充,滿足交直流互聯(lián)配電網(wǎng)內(nèi)交流負(fù)荷和直流負(fù)荷的總的供電需求。

    3) 利用儲能裝置進(jìn)行充放電以調(diào)整功率,當(dāng)配電網(wǎng)中電源的出力大于負(fù)荷所需功率時(shí),儲能進(jìn)行充電;當(dāng)電源出力之和不滿足負(fù)荷需求時(shí),儲能進(jìn)行放電。

    4) 當(dāng)配電網(wǎng)中電源和儲能裝置在出力不能滿足配電網(wǎng)負(fù)荷需求時(shí),配電網(wǎng)向交流大電網(wǎng)購電;反之則向大電網(wǎng)售電。

    在實(shí)現(xiàn)系統(tǒng)功率平衡的過程中,儲能和換流器有著重要的作用,所以在建立交直流互聯(lián)配電網(wǎng)經(jīng)濟(jì)調(diào)度模型時(shí),應(yīng)考慮儲能系統(tǒng)運(yùn)行成本和換流器損耗成本。

    2 "分布式電源數(shù)學(xué)模型

    2.1 "風(fēng)力發(fā)電機(jī)出力模型

    在交直流互聯(lián)配電網(wǎng)中,風(fēng)力發(fā)電機(jī)(Wind Turbine, WT)的出力大小與風(fēng)力發(fā)電機(jī)的瞬時(shí)風(fēng)速有關(guān)。風(fēng)速概率采用威布爾分布模擬,其概率密度函數(shù)如式(1)所示,風(fēng)電輸出功率與風(fēng)速關(guān)系式如式(2)所示[16?17]。

    [f(v)=βαvαβ-1e-vαβ] (1)

    [PWT(v)=Pr, vrlt;vlt;voPrv-vivr-vi, vilt;vlt;vr0, v≤vi 或 v≥vo] (2)

    式中:[α]和[β]分別為威布爾分布的尺度和形狀參數(shù);[vi]為風(fēng)機(jī)的切入風(fēng)速;[vo]為風(fēng)機(jī)的切出風(fēng)速;[vr]為額定風(fēng)速;[v]是風(fēng)機(jī)的輪轂高處的實(shí)際風(fēng)速;[Pr]表示風(fēng)機(jī)的額定功率。

    2.2 "光伏發(fā)電系統(tǒng)出力模型

    光伏發(fā)電系統(tǒng)(Photovoltaic Power Generation System, PV)出力大小與太陽輻射強(qiáng)度、環(huán)境溫度、電池表面溫度以及相對光學(xué)質(zhì)量有關(guān)。光伏發(fā)電系統(tǒng)的出力可由額定條件下的輻射強(qiáng)度、環(huán)境溫度、輸出功率得到,其出力表達(dá)式[18?20]如下:

    [PPV=PSTGCGST1+ε(TC-TST)] (3)

    式中:[PST]為光伏發(fā)電系統(tǒng)在額定條件下的額定輸出功率;[PPV]為光伏發(fā)電系統(tǒng)的實(shí)際輸出功率;標(biāo)準(zhǔn)條件[GST=1 kW/m2];[TST]為電池溫度,取值為[TST=25 ℃];[GC]為太陽輻射強(qiáng)度;[TC]為蓄電池發(fā)電時(shí)其表面的實(shí)際溫度,取值為[TC]=(48±2)℃;[ε]為功率溫度系數(shù),取值為[ε=-0.004 2 ℃]。

    2.3 "儲能系統(tǒng)模型

    交直流混合配電網(wǎng)中鉛酸蓄電池作為儲能裝置被廣泛使用,其具有“低儲高發(fā)”的作用,可提高新能源的消納率,抑制電網(wǎng)波動。儲能裝置運(yùn)行過程中的荷電狀態(tài)應(yīng)滿足如下關(guān)系[21]:

    [S(t+1)=S(t)(1-σ)+ηcPc(t)Δt , "充電S(t)(1-σ)+Pd(t)ηdΔt , " 放電] (4)

    式中:[S(t+1)]代表t+1時(shí)段內(nèi)蓄電池的荷電狀態(tài)(SOC)值;σ代表蓄電池自放電功率;[ηc]為蓄電池充電效率;[ηd]為蓄電池放電效率;[Pc(t)]、[Pd(t)]分別為t時(shí)段儲能蓄電池充電、放電功率;[Δt]為一個(gè)調(diào)度周期的時(shí)間步長,其中自放電功率一般很小,忽略不計(jì)。

    3 "配電網(wǎng)經(jīng)濟(jì)調(diào)度模型

    3.1 "目標(biāo)函數(shù)

    配電網(wǎng)優(yōu)化調(diào)度的目標(biāo)函數(shù)如下:

    [minf=f1+f2+f3+f4 ] (5)

    式中:[f1]表示設(shè)備運(yùn)行維護(hù)成本;[f2]表示儲能系統(tǒng)運(yùn)行成本;[f3]表示交直流雙向換流器損耗成本;[f4]表示配電網(wǎng)和大電網(wǎng)之間購售電成本。

    1) 設(shè)備運(yùn)行維護(hù)成本

    交直流互聯(lián)配電網(wǎng)運(yùn)行在不同狀態(tài)時(shí),由于設(shè)備的損耗或年限增加,會出現(xiàn)需要維修系統(tǒng)的情況。維修設(shè)備所花費(fèi)的維修費(fèi)用利用t時(shí)段的功率大小與設(shè)備運(yùn)行單位損耗成本系數(shù)表示。本文考慮光伏系統(tǒng)、風(fēng)力發(fā)電機(jī)組、微型燃?xì)廨啓C(jī)和AC/DC換流器的維護(hù)成本,其成本函數(shù)的表達(dá)式如下:

    [f1=t=124αPVPPV(t)+αWTPWT(t)+αMTPMT(t)+αTPT(t)] (6)

    式中:[αPV]、[αWT]、[αMT]和[αT]分別為光伏系統(tǒng)、風(fēng)力發(fā)電機(jī)組、微型燃?xì)廨啓C(jī)及VSC的設(shè)備單位運(yùn)行維護(hù)成本系數(shù);[PPV(t)]為t時(shí)段光伏發(fā)電系統(tǒng)出力;[PWT(t)]為t時(shí)段風(fēng)力發(fā)電機(jī)組出力;[PMT(t)]為t時(shí)段微型燃?xì)廨啓C(jī)出力;[PT(t)]為t時(shí)段VSC的注入功率。

    2) 儲能系統(tǒng)運(yùn)行成本

    [f2=t=124KBTPBT(t)] (7)

    式中:[PBT(t)]為t時(shí)段儲能的平均充放電功率,當(dāng)儲能充電時(shí),[PBT(t)]為負(fù);當(dāng)儲能放電時(shí),[PBT(t)]為正。[KBT]為單位儲能充放電的成本系數(shù)[22],取值為0.01。

    3) 交直流雙向換流器損耗成本

    VSC在運(yùn)行過程中,傳輸功率在交流側(cè)和直流側(cè)之間雙向流動,會存在一定功率損耗。在t時(shí)段內(nèi)其損耗成本與VSC傳輸?shù)墓β蔥PT(t)]之間的關(guān)系式為:

    [f3=mt·(1-ηt)·t=124PT(t)] (8)

    式中:[mt]為VSC損耗成本系數(shù);[ηt]表示換流效率。當(dāng)VSC注入功率為從直流側(cè)流入交流側(cè)時(shí),[PT(t)]為正;注入功率從交流側(cè)流入直流側(cè)時(shí),[PT(t)]為負(fù)。所以,計(jì)算換流器損耗成本時(shí)換流功率取絕對值。

    4) 配電網(wǎng)與大電網(wǎng)之間購售電成本

    [f4=t=1NCgrid(t)Pgrid(t)] (9)

    式中[Pgrid(t)]表示t時(shí)段配電網(wǎng)向大電網(wǎng)購售電的功率。當(dāng)配電網(wǎng)向大電網(wǎng)購電時(shí),[Pgrid(t)≥0];當(dāng)配電網(wǎng)向大電網(wǎng)售電時(shí),[Pgrid(t)lt;0]。

    [Cgrid(t)=Cbuy(t), " "Pgrid≥0Csell(t), " "Pgridlt;0] (10)

    式中[Cgrid(t)]表示t時(shí)段配電網(wǎng)與大電網(wǎng)之間的交易價(jià)格。當(dāng)配電網(wǎng)分向大電網(wǎng)購電時(shí),選擇分時(shí)購電價(jià)格[Cbuy(t)];當(dāng)配電網(wǎng)分向大電網(wǎng)售電時(shí),選擇分時(shí)售電價(jià)格[Csell(t)]。

    3.2 "約束條件

    1) 有功功率平衡約束

    系統(tǒng)的功率平衡約束如式(11)所示,交流側(cè)和直流側(cè)的有功功率平衡約束分別如式(12)和式(13)所示。

    [PWT(t)+PPV(t)+Pgrid(t)+Pd(t)-Pc(t)= " " " "PLAC(t)+PLDC(t)+(1-ηt)PT(t)] (11)

    [PWT(t)+PMT(t)+Pgrid(t)=PLAC(t)+PACT(t)-ηtPDCT(t)] (12)

    [PPV(t)+Pd(t)-Pc(t)=PLDC(t)+PDCT(t)-ηtPACT(t)] (13)

    式中:[PLAC(t)]表示t時(shí)段配電網(wǎng)交流負(fù)荷功率;[PLDC(t)]表示t時(shí)段配電網(wǎng)直流負(fù)荷功率;[PACT(t)]和[PDCT(t)]分別為t時(shí)段交流側(cè)和直流側(cè)注入VSC的功率。

    2) 光伏系統(tǒng)和風(fēng)力發(fā)電機(jī)的出力約束

    交直流互聯(lián)配電網(wǎng)中光伏系統(tǒng)和風(fēng)力發(fā)電機(jī)的出力應(yīng)小于最大輸出功率,其約束條件如下:

    [0≤PPV(t)≤PmaxPV(t)] (14)

    [0≤PWT(t)≤PmaxWT(t)] (15)

    式中:[PmaxPV(t)]、[PmaxWT(t)]分別為光伏系統(tǒng)和風(fēng)力發(fā)電機(jī)在t時(shí)段的最大輸出功率。

    3) 微型燃?xì)廨啓C(jī)的功率約束

    [PMT_min≤PMT(t)≤PMT_max] (16)

    式中:[PMT_min]、[PMT_max]分別為微型燃?xì)廨啓C(jī)出力的下限和上限。

    4) VSC功率約束

    交直流混合配電網(wǎng)中連接交直流兩側(cè)的VSC在功率交互的過程中必須滿足其最大交互功率約束。

    [PT(t)≤PT_max] (17)

    [UACBC,t+UDCBC,t≤1] (18)

    式中:[PT_max]為VSC的交換功率最大允許值;[UACBC,t]、[UDCBC,t]分別為t時(shí)段VSC傳輸功率的流向情況。兩者均為0?1變量,[UACBC,t=1]表示t時(shí)段VSC傳輸功率由交流側(cè)流向直流側(cè);[UDCBC,t=1]表示t時(shí)段VSC傳輸功率由直流側(cè)流向交流側(cè)。

    5) 儲能裝置約束

    [μc,t+μd,t≤1μc,tμd,t∈0,1] (19)

    [SOCmin≤SOC(t)≤SOCmaxSOC0=SOCend] (20)

    [Pc(t)≤μc,tPmaxcPd(t)≤μd,tPmaxd] (21)

    式中:[μc,t]、[μd,t]都為0?1變量,表示t時(shí)段儲能裝置的充放電狀態(tài);SOC(t)為t時(shí)段的荷電狀態(tài);[SOCmin]表示儲能蓄電池最小荷電狀態(tài);[SOCmax]表示儲能蓄電池最大荷電狀態(tài);儲能電池在一次充放電的起始時(shí)刻和終止時(shí)刻的荷電狀態(tài)相等,即[SOC0=SOCend];[Pmaxc]、[Pmaxd]分別為儲能電池充放電的最大功率。

    6) 配電網(wǎng)向大電網(wǎng)購售電功率約束

    [0≤Pbuy(t)≤λgPmaxgrid] (22)

    [0≤Psell(t)≤(1-λg)Pmaxgrid] (23)

    式中:[Pbuy(t)]和[Psell(t)]分別表示t時(shí)段配電網(wǎng)向大電網(wǎng)購電和售電的有功功率;[λg]為0?1變量,當(dāng)[λg]為1時(shí),表示配電網(wǎng)向大電網(wǎng)購電;[Pmaxgrid]為配電網(wǎng)向大電網(wǎng)購電或售電的最大功率。

    4 "增強(qiáng)型鯨魚優(yōu)化算法

    4.1 "鯨魚優(yōu)化算法

    假設(shè)鯨魚捕捉的獵物為最優(yōu)解,鯨魚所在位置為潛在解,在每次迭代的過程中,每頭鯨魚的位置更新策略通過隨機(jī)數(shù)p的取值和系數(shù)向量A的模來決定,隨著迭代次數(shù)更新,鯨魚種群不斷靠近最優(yōu)解。鯨魚位置的更新公式如下:

    [X(t+1)=X*(t)-A?D , plt;0.5XDp?ebl?cos(2πl(wèi))+X*(t), p≥0.5 ] (24)

    式中:t為迭代次數(shù);p為[0,1]之間的隨機(jī)數(shù),當(dāng)p≥0.5時(shí),采用螺旋更新的方式,當(dāng)plt;0.5時(shí),根據(jù)種群中其他鯨魚的位置進(jìn)行更新;l為[-1,1]之間的隨機(jī)數(shù);b為限定對數(shù)螺旋形狀的常數(shù),一般取值為1即可;[X*(t)]表示獵物位置;D和[Dp]分別由式(25)和式(26)得到。

    [D=C?X*(t)-X(t)] (25)

    [Dp=X*(t)-X(t)] (26)

    式(24)中的A和式(25)中的C為系數(shù)向量,可由式(27)和式(28)計(jì)算而來。

    [A=2a1r1-2a1] (27)

    [C=2r2] (28)

    式中:[r1]、[r2]為[0,1]之間的隨機(jī)向量;[a1]表示收斂因子,隨著迭代次數(shù)的增加,從2線性減小到0。[a1]表達(dá)式如下:

    [a1=2-2TTmax] (29)

    式中:T為當(dāng)前迭代次數(shù);[Tmax]為最大迭代次數(shù)。

    當(dāng)[plt;0.5、A≥1]時(shí),隨機(jī)選擇種群中一頭鯨魚的位置進(jìn)行更新,計(jì)算公式如下:

    [X(t+1)=Xrand(t)-A?D] (30)

    [D=C?Xrand(t)-X(t)] (31)

    式中[Xrand]表示種群中隨機(jī)選擇的位置。當(dāng)[plt;0.5、Alt;1]時(shí),用歷史最優(yōu)鯨魚的位置更新。

    4.2 "增強(qiáng)型鯨魚優(yōu)化算法

    WOA算法存在求解速度慢、求解精度較低、容易陷入局部最優(yōu)的缺點(diǎn)。因此,本文提出EWOA算法對WOA算法進(jìn)行改進(jìn),即從限定對數(shù)螺旋形狀的常數(shù)考慮,將b值設(shè)成動態(tài)值,并隨著迭代次數(shù)更新,隨機(jī)地改變螺旋路徑。在迭代開始時(shí)期,螺旋形狀偏大,可以增大鯨魚的搜索空間,使遍歷種群中的個(gè)體增多,增強(qiáng)算法的全局尋優(yōu)能力;在迭代后期,隨著螺線形狀逐漸減小,鯨魚會在最優(yōu)個(gè)體的附近進(jìn)行搜索,從而加快算法的速度,同時(shí)提高算法的尋優(yōu)精度。對限定對數(shù)螺旋形狀的系數(shù)b進(jìn)行更新的公式為:

    [b=0.5-0.5TTmax] (32)

    EWOA的執(zhí)行步驟如下。

    1) 初始化算法參數(shù),即鯨魚種群規(guī)模大小N以及最大迭代次數(shù)[Tmax]。

    2) 初始化鯨魚種群的位置。

    3) 計(jì)算每頭鯨魚的相應(yīng)適應(yīng)度值,按照適應(yīng)度值的大小排序。

    4) 判斷是否達(dá)到終止條件。若滿足終止條件,則輸出最優(yōu)個(gè)體,即算法仿真找到的最優(yōu)解;否則,根據(jù)式(27)~式(29)分別計(jì)算參數(shù)A、C和[a1],根據(jù)式(32)更新b值,以更新鯨魚位置。

    5) 計(jì)算出種群適應(yīng)度值的大小,并將適應(yīng)度值最小的個(gè)體位置作為種群的最優(yōu)位置。

    6) 更新迭代次數(shù),返回步驟4)。

    EWOA算法的流程如圖2所示。

    5 "算例仿真結(jié)果分析

    5.1 "算例相關(guān)參數(shù)

    選擇圖1所示的并網(wǎng)型交直流混合配電網(wǎng)作為研究對象進(jìn)行調(diào)度分析。風(fēng)力發(fā)電機(jī)、光伏發(fā)電系統(tǒng)以及微型燃?xì)廨啓C(jī)的參數(shù)如表1所示。其中,風(fēng)力發(fā)電機(jī)的額定風(fēng)速為13 m/s,切入風(fēng)速為3 m/s,切出風(fēng)速為25 m/s;光伏發(fā)電系統(tǒng)的發(fā)電效率為0.9。

    系統(tǒng)內(nèi)的儲能電池裝置容量為300 kW·h,充電功率最大值為30 kW,放電功率最大值為45 kW;VSC換流器傳輸功率的最大值為100 kW,換流效率為0.95;交直流互聯(lián)配電網(wǎng)與大電網(wǎng)交互功率的最大值為100 kW。儲能裝置的單位損耗成本為0.01元/(kW·h);VSC換流器的單位損耗成本為0.4元/(kW·h)。

    典型日風(fēng)速與輻射強(qiáng)度和風(fēng)電光伏日運(yùn)行曲線如圖3所示。

    交直流負(fù)荷日運(yùn)行功率曲線如圖4所示。

    5.2 "結(jié)果分析

    本文分別對EWOA算法、WOA算法、TLBO算法[23]進(jìn)行仿真計(jì)算,三種算法的種群規(guī)模設(shè)置為50,最大迭代次數(shù)為100。種群規(guī)模和最大迭代次數(shù)的值根據(jù)經(jīng)驗(yàn)設(shè)置,若設(shè)置的值較小,則無法達(dá)到較好的調(diào)度結(jié)果;若設(shè)置的值過大,則造成計(jì)算資源浪費(fèi)。通過算例仿真驗(yàn)證,并對三種算法的計(jì)算結(jié)果進(jìn)行對比分析。

    5.2.1 "優(yōu)化調(diào)度結(jié)果分析

    通過EWOA算法仿真得出各單元出力的調(diào)度策略,對調(diào)度結(jié)果進(jìn)行分析,其功率平衡關(guān)系調(diào)度結(jié)果如圖5所示。

    由圖5可知,0:00—5:00時(shí)段,用戶用電量處于低谷期,風(fēng)速較大,風(fēng)力發(fā)電充足,但在此時(shí)段輻射強(qiáng)度為0,光伏系統(tǒng)出力為0。此時(shí)風(fēng)力發(fā)電和微型燃?xì)廨啓C(jī)滿足用戶用電的量,且發(fā)電多出用電的量向電網(wǎng)進(jìn)行售電,同時(shí)儲能系統(tǒng)進(jìn)行充電。6:00—12:00時(shí)段,負(fù)荷量開始增多,風(fēng)速減小,光照增強(qiáng)。此時(shí)段風(fēng)光和燃?xì)廨啓C(jī)發(fā)電量小于負(fù)荷需求量,配電網(wǎng)向大電網(wǎng)購電,儲能系統(tǒng)放電供用戶使用。13:00—16:00負(fù)荷量減少,輻射強(qiáng)度較大,此時(shí)段燃?xì)廨啓C(jī)發(fā)電量、光伏和風(fēng)電發(fā)電量大于負(fù)荷量,配電網(wǎng)向大電網(wǎng)售電,儲能系統(tǒng)進(jìn)行充電。17:00—19:00時(shí)段,輻射強(qiáng)度減小,風(fēng)速增大,負(fù)荷量處于高峰時(shí)段,燃?xì)廨啓C(jī)發(fā)電量、光伏和風(fēng)電發(fā)電量小于用戶用電量,系統(tǒng)向大電網(wǎng)購電,儲能系統(tǒng)放電。20:00—22:00時(shí)段,輻射強(qiáng)度為0,風(fēng)速增大,負(fù)荷量逐漸減小,風(fēng)力發(fā)電和微型燃?xì)廨啓C(jī)出力之和等于用戶用電量。23:00—24:00時(shí)段,用電量減少,系統(tǒng)發(fā)電量大于負(fù)荷量,配電網(wǎng)向大電網(wǎng)售電,儲能裝置進(jìn)行充電。由此可得:儲能系統(tǒng)在負(fù)荷低谷時(shí)進(jìn)行充電,吸收多余電量;在負(fù)荷高峰時(shí)進(jìn)行放電,具有“削峰填谷”的作用,以維持系統(tǒng)的功率平衡。

    交流側(cè)功率差為風(fēng)力發(fā)電功率及微型燃?xì)廨啓C(jī)出力與交流負(fù)荷的差值,其表達(dá)式如式(33)所示;直流側(cè)功率差為光伏系統(tǒng)發(fā)電及儲能放電功率與直流負(fù)荷的差值,其表達(dá)式如式(34)所示。

    [ΔPAC(t)=PWT(t)+PMT(t)-PLAC(t)] (33)

    [ΔPDC(t)=PPV(t)+PBT(t)-PLDC(t)] (34)

    式中:[ΔPAC(t)]為交流側(cè)功率差,若[ΔPAC(t)]為正,則t時(shí)段風(fēng)力發(fā)電和微型燃?xì)廨啓C(jī)出力之和大于交流負(fù)荷,反之小于交流負(fù)荷;[ΔPDC(t)]為直流側(cè)功率差,若[ΔPDC(t)]為正,則t時(shí)段光伏發(fā)電和儲能放電功率之和大于直流負(fù)荷,反之則小于直流負(fù)荷。交直流側(cè)功率流向關(guān)系如圖6所示。

    由圖6交直流側(cè)功率流向關(guān)系可知:交直流側(cè)功率流向關(guān)系曲線和直流側(cè)功率差曲線趨勢一致,和交流側(cè)功率差曲線趨勢形成互補(bǔ);當(dāng)直流側(cè)功率差為正時(shí)交流側(cè)功率差為負(fù),則VSC的傳輸功率為直流側(cè)流向交流側(cè),反之則由交流側(cè)流向直流側(cè)。仿真結(jié)果符合實(shí)際。

    5.2.2 "算法對比分析

    EWOA算法與WOA算法、TLBO算法的優(yōu)化調(diào)度結(jié)果對比如表2所示。

    由表2可知,EWOA算法得到的綜合成本為697.097元,WOA算法得到綜合成本為848.882元,TLBO算法得到的綜合成本為805.749元。EWOA算法的綜合成本比WOA的計(jì)算結(jié)果減少了17.9%,比TLBO計(jì)算結(jié)果減少了13.5%,表明EWOA的優(yōu)化效果優(yōu)于WOA和TLBO。幾種算法的計(jì)算用時(shí)如表3所示。表中,EWOA算法仿真用時(shí)比WOA算法減少15.95 s,比TLBO減少7.979 s,計(jì)算速度最快。

    6 "結(jié) "論

    本文為減小交直流互聯(lián)配電網(wǎng)的運(yùn)行成本,構(gòu)建以設(shè)備運(yùn)行維護(hù)成本、儲能系統(tǒng)運(yùn)行成本、交直流雙向換流器損耗成本以及配電網(wǎng)與大電網(wǎng)之間購售電成本之和最小的經(jīng)濟(jì)調(diào)度模型,在滿足配電網(wǎng)安全穩(wěn)定運(yùn)行約束的條件下進(jìn)行優(yōu)化調(diào)度。文中提出一種有效求解模型的EWOA算法,將對數(shù)螺旋形狀的常數(shù)b進(jìn)行動態(tài)更新后擴(kuò)大了算法的搜索空間,增強(qiáng)了算法的全局尋優(yōu)能力,加快了收斂速度。算例仿真結(jié)果表明,與鯨魚優(yōu)化算法和教與學(xué)優(yōu)化算法相比,EWOA算法獲得的運(yùn)行成本分別減少了17.9%和13.5%,計(jì)算效率分別提升了44.7%和22.4%,表現(xiàn)出了很好的優(yōu)化性能。

    注:本文通訊作者為熊國江。

    參考文獻(xiàn)

    [1] 侯帥輝,張新燕,代玥,等.計(jì)及儲能配合與需求側(cè)響應(yīng)的主動配電網(wǎng)分布式電源規(guī)劃[J].現(xiàn)代電子技術(shù),2022,45(23):132?137.

    [2] 魏斌,韓肖清,李雯,等.融合多場景分析的交直流混合微電網(wǎng)多時(shí)間尺度隨機(jī)優(yōu)化調(diào)度策略[J].高電壓技術(shù),2020,46(7):2359?2369.

    [3] 鄭重,苗世洪,李超,等.面向微型能源互聯(lián)網(wǎng)接入的交直流配電網(wǎng)協(xié)同優(yōu)化調(diào)度策略[J].電工技術(shù)學(xué)報(bào),2022,37(1):192?207.

    [4] 朱建昆,高紅均,賀帥佳,等.考慮VSC與光?儲?充協(xié)同配置的交直流混合配電網(wǎng)規(guī)劃[J].智慧電力,2023,51(11):7?14.

    [5] 王揚(yáng),蔡永翔,付宇,等.適應(yīng)分布式光伏并網(wǎng)的低壓交直流混合配電工程示范應(yīng)用[J].供用電,2023,40(3):40?46.

    [6] "QI C, WANG K, FU Y, et al. A decentralized optimal operation of AC/DC hybrid distribution grids [J].IEEE transaction on smart grid, 2017, 9(6): 6095?6105.

    [7] 張盧珍,田易之.基于JAYA算法的交直流柔性配電網(wǎng)穩(wěn)態(tài)安全域研究[J].現(xiàn)代電子技術(shù),2023,46(18):113?118.

    [8] 孫碩愷,張妍婷,戚星宇.考慮源荷不確定性的交直流配電臺區(qū)日前優(yōu)化調(diào)度[J].供用電,2023,40(1):25?32.

    [9] 孟明,魏怡,劉晗,等.含三端柔性環(huán)網(wǎng)裝置的交直流混合配電網(wǎng)分層控制策略研究[J].華北電力大學(xué)學(xué)報(bào)(自然科學(xué)版),2019,46(3):9?16.

    [10] 王文超,龐丹,成龍,等.考慮電價(jià)型需求響應(yīng)的交直流混合配電網(wǎng)優(yōu)化調(diào)度[J].電網(wǎng)技術(shù),2019,43(5):1675?1682.

    [11] GAO S, LIU S, LIU Y, et al. Flexible and economic dispatching of AC/DC distribution networks considering uncertainty of wind power [J]. IEEE access, 2019, 7: 100051?100065.

    [12] ZHAO D, WANG H, TAO R. Multi?time scale dispatch approach for an AC/DC hybrid distribution system considering the response uncertainty of flexible loads [J]. Electric power systems research, 2021, 199: 107394.

    [13] 賈善翔,彭克,李喜東,等.基于智能軟開關(guān)的柔性互聯(lián)交直流配電系統(tǒng)分布式優(yōu)化[J].電力系統(tǒng)自動化,2021,45(8):95?103.

    [14] 梁海平,王巖,劉英培,等.計(jì)及源荷不確定性的混合交直流主動配電網(wǎng)分層?分布式優(yōu)化調(diào)度[J].電力自動化設(shè)備,2021,41(12):62?69.

    [15] MIRJALILI S, LEWIS A. The whale optimization algorithm [J]. Advances in engineering software, 2016, 95: 51?67.

    [16] 王秀茹,劉剛,黃華峰,等.考慮分布式電源的配電網(wǎng)無功調(diào)度和儲能優(yōu)化方法[J].電力科學(xué)與技術(shù)學(xué)報(bào),2022,37(4):134?142.

    [17] 劉一格,趙振宙,馬遠(yuǎn)卓,等.基于鯨魚優(yōu)化算法的串列風(fēng)力機(jī)主動尾流控制策略[J].中國電機(jī)工程學(xué)報(bào),2024,44(9):3702?3710.

    [18] 裴蕾,衛(wèi)志農(nóng),陳勝,等.計(jì)及光伏出力波動的交直流混合配電網(wǎng)安全域模型[J].電網(wǎng)技術(shù),2021,45(10):4084?4093.

    [19] 鄧長征,馮朕,邱立,等.基于混沌灰狼算法的交直流混合微網(wǎng)經(jīng)濟(jì)調(diào)度[J].電測與儀表,2020,57(4):99?107.

    [20] 楊佳霖,趙鵬翔,竇真蘭,等.計(jì)及經(jīng)濟(jì)性與可靠性的綜合能源系統(tǒng)優(yōu)化調(diào)度[J].供用電,2024,41(1):90?99.

    [21] 邱海峰,趙波,林達(dá),等.計(jì)及儲能損耗和換流成本的交直流混合微網(wǎng)區(qū)域協(xié)調(diào)調(diào)度[J].電力系統(tǒng)自動化,2017,41(23):29?37.

    [22] 張雯雯,魏震波,郭毅,等.含混合儲能的交直流配電網(wǎng)日經(jīng)濟(jì)優(yōu)化運(yùn)行[J].高電壓技術(shù),2022,48(2):565?576.

    [23] RAO R V, SAVSANI V J, VAKHARIA D P. Teaching?learning?based optimization: a novel method for constrained mechanical design optimization problems [J]. Computer?aided design, 2011, 43(3): 303?315.

    精品国产亚洲在线| 美女高潮到喷水免费观看| 在线亚洲精品国产二区图片欧美| 精品福利观看| 午夜日韩欧美国产| 多毛熟女@视频| 十分钟在线观看高清视频www| 91国产中文字幕| 老司机影院毛片| 精品乱码久久久久久99久播| 99re在线观看精品视频| 欧美亚洲日本最大视频资源| 一区二区三区乱码不卡18| 国产激情久久老熟女| 成人手机av| 日韩视频一区二区在线观看| kizo精华| 国产亚洲精品第一综合不卡| 日本黄色视频三级网站网址 | 黄色丝袜av网址大全| 国产欧美日韩综合在线一区二区| 国产福利在线免费观看视频| 国产精品香港三级国产av潘金莲| 少妇精品久久久久久久| 黄色 视频免费看| 亚洲五月色婷婷综合| 欧美人与性动交α欧美软件| 王馨瑶露胸无遮挡在线观看| 国产av国产精品国产| 大码成人一级视频| 在线观看免费午夜福利视频| 亚洲av国产av综合av卡| av有码第一页| 人人妻人人爽人人添夜夜欢视频| 丁香六月天网| 午夜精品国产一区二区电影| a级毛片黄视频| 中文字幕精品免费在线观看视频| 中文字幕av电影在线播放| 淫妇啪啪啪对白视频| 在线播放国产精品三级| 老司机午夜福利在线观看视频 | 久久影院123| 久久久久久亚洲精品国产蜜桃av| 麻豆国产av国片精品| 在线永久观看黄色视频| 久久久国产欧美日韩av| 三上悠亚av全集在线观看| 成人特级黄色片久久久久久久 | 麻豆乱淫一区二区| 五月天丁香电影| 午夜91福利影院| 精品高清国产在线一区| 日日爽夜夜爽网站| 怎么达到女性高潮| 人人澡人人妻人| www.熟女人妻精品国产| 国产亚洲欧美精品永久| 777米奇影视久久| 欧美另类亚洲清纯唯美| 99精品久久久久人妻精品| 亚洲国产av影院在线观看| 国产精品亚洲av一区麻豆| 亚洲第一欧美日韩一区二区三区 | 久久久久国产一级毛片高清牌| 别揉我奶头~嗯~啊~动态视频| 久久天躁狠狠躁夜夜2o2o| 日韩人妻精品一区2区三区| 又大又爽又粗| 一区二区三区精品91| 国产有黄有色有爽视频| 狠狠婷婷综合久久久久久88av| 欧美日韩中文字幕国产精品一区二区三区 | 成人亚洲精品一区在线观看| 人妻 亚洲 视频| 国产xxxxx性猛交| 日韩精品免费视频一区二区三区| 国产高清激情床上av| 精品视频人人做人人爽| 99re在线观看精品视频| 国产精品麻豆人妻色哟哟久久| 日韩欧美国产一区二区入口| 欧美日韩亚洲高清精品| 每晚都被弄得嗷嗷叫到高潮| 曰老女人黄片| 亚洲五月色婷婷综合| 妹子高潮喷水视频| 亚洲国产中文字幕在线视频| 国产单亲对白刺激| 最新在线观看一区二区三区| 免费不卡黄色视频| 国产97色在线日韩免费| 在线av久久热| e午夜精品久久久久久久| 精品久久久久久电影网| 日韩人妻精品一区2区三区| 亚洲七黄色美女视频| 黄频高清免费视频| 欧美日韩亚洲综合一区二区三区_| 露出奶头的视频| 美女高潮到喷水免费观看| 亚洲av日韩精品久久久久久密| 精品人妻熟女毛片av久久网站| 十八禁网站免费在线| 国产伦理片在线播放av一区| 在线天堂中文资源库| 久久久国产精品麻豆| 黄色视频,在线免费观看| 99国产综合亚洲精品| 在线看a的网站| 日韩人妻精品一区2区三区| 亚洲黑人精品在线| 满18在线观看网站| 久久亚洲真实| 嫁个100分男人电影在线观看| 麻豆国产av国片精品| 日本wwww免费看| 99久久精品国产亚洲精品| av片东京热男人的天堂| 人人妻人人澡人人看| 视频在线观看一区二区三区| 亚洲av国产av综合av卡| 大型av网站在线播放| 久久人妻熟女aⅴ| 大型av网站在线播放| 狠狠狠狠99中文字幕| 男女之事视频高清在线观看| 男女下面插进去视频免费观看| 69av精品久久久久久 | 精品少妇久久久久久888优播| 久久人人爽av亚洲精品天堂| 一级毛片电影观看| 成人18禁在线播放| 老司机福利观看| av又黄又爽大尺度在线免费看| 精品久久蜜臀av无| 狂野欧美激情性xxxx| 国产亚洲午夜精品一区二区久久| 高清黄色对白视频在线免费看| 国产在线免费精品| 人人妻人人爽人人添夜夜欢视频| tube8黄色片| 国产男靠女视频免费网站| 日韩制服丝袜自拍偷拍| 亚洲专区中文字幕在线| 国产亚洲一区二区精品| 国产精品亚洲一级av第二区| 久久国产精品影院| 在线观看www视频免费| 久久久久久久国产电影| 亚洲成人免费电影在线观看| 国产精品久久久久久精品古装| 激情视频va一区二区三区| 亚洲精品中文字幕在线视频| 午夜精品国产一区二区电影| 亚洲成人免费av在线播放| 亚洲美女黄片视频| 国产又色又爽无遮挡免费看| 深夜精品福利| 精品人妻熟女毛片av久久网站| 大香蕉久久网| 久久中文字幕一级| 人成视频在线观看免费观看| 性色av乱码一区二区三区2| 香蕉丝袜av| 精品福利永久在线观看| 午夜福利在线观看吧| 亚洲色图 男人天堂 中文字幕| 99久久国产精品久久久| 性少妇av在线| 亚洲avbb在线观看| 国产男女内射视频| 精品一区二区三区av网在线观看 | 老熟妇乱子伦视频在线观看| 一级a爱视频在线免费观看| 一边摸一边做爽爽视频免费| 中文字幕制服av| 老熟妇仑乱视频hdxx| 午夜福利视频在线观看免费| 啦啦啦免费观看视频1| 黄片小视频在线播放| 色94色欧美一区二区| 交换朋友夫妻互换小说| 这个男人来自地球电影免费观看| 成年人黄色毛片网站| 午夜福利欧美成人| 999久久久国产精品视频| 91成人精品电影| 国产不卡一卡二| 99在线人妻在线中文字幕 | 国产野战对白在线观看| tube8黄色片| 亚洲久久久国产精品| 男女之事视频高清在线观看| 高清毛片免费观看视频网站 | e午夜精品久久久久久久| 国产精品亚洲av一区麻豆| 亚洲伊人久久精品综合| 欧美中文综合在线视频| 男女高潮啪啪啪动态图| 国精品久久久久久国模美| 不卡av一区二区三区| 亚洲专区国产一区二区| 日韩一区二区三区影片| 成人永久免费在线观看视频 | 亚洲天堂av无毛| 欧美 亚洲 国产 日韩一| 一本一本久久a久久精品综合妖精| 丰满人妻熟妇乱又伦精品不卡| 大码成人一级视频| 人妻久久中文字幕网| 精品少妇久久久久久888优播| 大片免费播放器 马上看| 午夜福利影视在线免费观看| 香蕉久久夜色| 777久久人妻少妇嫩草av网站| 成人特级黄色片久久久久久久 | 一进一出好大好爽视频| 91精品三级在线观看| 精品久久久久久久毛片微露脸| 电影成人av| 操出白浆在线播放| 2018国产大陆天天弄谢| 中文字幕av电影在线播放| 在线十欧美十亚洲十日本专区| 777久久人妻少妇嫩草av网站| 日韩中文字幕欧美一区二区| av线在线观看网站| 色在线成人网| 国产精品久久久人人做人人爽| 欧美精品啪啪一区二区三区| av国产精品久久久久影院| 757午夜福利合集在线观看| 久久久久久久国产电影| 久久精品aⅴ一区二区三区四区| 久久av网站| 成人亚洲精品一区在线观看| 国产精品av久久久久免费| 欧美亚洲 丝袜 人妻 在线| 日本wwww免费看| 天堂8中文在线网| 久久久久网色| 久久久久久免费高清国产稀缺| 国产成人欧美| 免费少妇av软件| 成人免费观看视频高清| 久久亚洲真实| 日本一区二区免费在线视频| 婷婷丁香在线五月| 高清av免费在线| 久久精品国产综合久久久| 在线观看www视频免费| 久久国产精品男人的天堂亚洲| www.999成人在线观看| 男女之事视频高清在线观看| 欧美日韩一级在线毛片| 欧美另类亚洲清纯唯美| 91精品国产国语对白视频| 亚洲熟女毛片儿| 国产高清激情床上av| 肉色欧美久久久久久久蜜桃| 国产成人av教育| 看免费av毛片| 五月天丁香电影| 亚洲精品国产区一区二| av欧美777| 国内毛片毛片毛片毛片毛片| 电影成人av| tocl精华| 亚洲精品国产一区二区精华液| 午夜福利视频精品| 亚洲成人免费电影在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产欧美日韩精品亚洲av| cao死你这个sao货| 99久久精品国产亚洲精品| 一级毛片女人18水好多| 欧美人与性动交α欧美软件| 精品人妻1区二区| 波多野结衣av一区二区av| 91成人精品电影| 欧美激情 高清一区二区三区| 美女高潮喷水抽搐中文字幕| 天天操日日干夜夜撸| 久久久久视频综合| 青青草视频在线视频观看| av不卡在线播放| 999久久久国产精品视频| 欧美午夜高清在线| 极品教师在线免费播放| 大香蕉久久成人网| 国产精品免费大片| 大片免费播放器 马上看| 亚洲成人免费电影在线观看| 成人三级做爰电影| www.999成人在线观看| 人成视频在线观看免费观看| 少妇的丰满在线观看| av国产精品久久久久影院| 9热在线视频观看99| 波多野结衣av一区二区av| 叶爱在线成人免费视频播放| 国产欧美日韩一区二区精品| 正在播放国产对白刺激| 老司机亚洲免费影院| 国产成人欧美| 午夜福利,免费看| 久久九九热精品免费| 久久人妻av系列| 飞空精品影院首页| 亚洲男人天堂网一区| a在线观看视频网站| 欧美国产精品一级二级三级| 99国产精品一区二区三区| 国产真人三级小视频在线观看| 捣出白浆h1v1| 人人妻人人澡人人看| 啦啦啦在线免费观看视频4| 国产成人欧美| 99国产精品免费福利视频| 久久天躁狠狠躁夜夜2o2o| 精品欧美一区二区三区在线| 色综合婷婷激情| 999久久久国产精品视频| 啦啦啦 在线观看视频| 99香蕉大伊视频| 久久国产亚洲av麻豆专区| 精品一区二区三区视频在线观看免费 | 亚洲色图综合在线观看| 曰老女人黄片| 在线永久观看黄色视频| 中亚洲国语对白在线视频| 久久久久精品国产欧美久久久| 在线观看免费视频日本深夜| 久久精品国产亚洲av香蕉五月 | 国产福利在线免费观看视频| 婷婷丁香在线五月| 免费看a级黄色片| 免费在线观看完整版高清| 制服诱惑二区| 久久青草综合色| 国产成人欧美| 十八禁网站免费在线| 国产欧美日韩综合在线一区二区| 国产av国产精品国产| 午夜精品国产一区二区电影| 十八禁网站网址无遮挡| 69精品国产乱码久久久| 亚洲色图综合在线观看| 久久久久久久大尺度免费视频| 久久精品人人爽人人爽视色| 日本一区二区免费在线视频| 在线观看免费高清a一片| 999精品在线视频| 啦啦啦 在线观看视频| 久久久精品国产亚洲av高清涩受| 丰满少妇做爰视频| 亚洲av成人不卡在线观看播放网| 精品少妇黑人巨大在线播放| 久久午夜亚洲精品久久| 黄色a级毛片大全视频| 色老头精品视频在线观看| 日韩欧美国产一区二区入口| 久久人妻福利社区极品人妻图片| 中文字幕人妻熟女乱码| 久久青草综合色| 成年人午夜在线观看视频| 国产又爽黄色视频| 精品久久久久久久毛片微露脸| 别揉我奶头~嗯~啊~动态视频| 亚洲 欧美一区二区三区| 国产成人欧美| 90打野战视频偷拍视频| 精品久久蜜臀av无| 黄色视频在线播放观看不卡| 国产精品久久久人人做人人爽| 青草久久国产| 脱女人内裤的视频| 亚洲全国av大片| 亚洲久久久国产精品| 久久精品亚洲av国产电影网| 久久99热这里只频精品6学生| 亚洲专区中文字幕在线| 男女之事视频高清在线观看| www.熟女人妻精品国产| 动漫黄色视频在线观看| 国产在线视频一区二区| 黄色片一级片一级黄色片| 国产在视频线精品| 少妇被粗大的猛进出69影院| 老鸭窝网址在线观看| 黄色丝袜av网址大全| 欧美亚洲日本最大视频资源| 亚洲人成77777在线视频| 亚洲综合色网址| 老司机靠b影院| 国内毛片毛片毛片毛片毛片| svipshipincom国产片| 一边摸一边做爽爽视频免费| 精品人妻熟女毛片av久久网站| 久久久精品国产亚洲av高清涩受| 久久精品亚洲精品国产色婷小说| 午夜福利欧美成人| 女警被强在线播放| 欧美国产精品一级二级三级| 国产日韩欧美视频二区| 精品乱码久久久久久99久播| 操出白浆在线播放| 老司机亚洲免费影院| 欧美黄色片欧美黄色片| 成人黄色视频免费在线看| 我的亚洲天堂| 国产精品久久久av美女十八| 日韩人妻精品一区2区三区| 脱女人内裤的视频| 午夜福利视频在线观看免费| 午夜福利在线观看吧| 人人澡人人妻人| 十八禁高潮呻吟视频| 男女下面插进去视频免费观看| 国产成人av激情在线播放| 老熟女久久久| 天天操日日干夜夜撸| 精品一区二区三区av网在线观看 | 欧美激情 高清一区二区三区| 国产成人av教育| 国产精品免费视频内射| 性高湖久久久久久久久免费观看| 大码成人一级视频| 亚洲中文字幕日韩| 免费观看人在逋| 精品一区二区三区av网在线观看 | 老汉色av国产亚洲站长工具| 18禁裸乳无遮挡动漫免费视频| 女同久久另类99精品国产91| 人人妻人人澡人人爽人人夜夜| 欧美大码av| 午夜免费成人在线视频| 18在线观看网站| 日韩免费高清中文字幕av| 91麻豆精品激情在线观看国产 | 精品亚洲成国产av| 久久久久国内视频| 国产成人免费无遮挡视频| 欧美在线黄色| 欧美日韩一级在线毛片| 男人操女人黄网站| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕av电影在线播放| 欧美黑人精品巨大| 精品久久久久久电影网| 1024视频免费在线观看| 欧美人与性动交α欧美精品济南到| 99精品在免费线老司机午夜| 亚洲美女黄片视频| 午夜福利在线观看吧| 成人18禁高潮啪啪吃奶动态图| 亚洲 国产 在线| 69av精品久久久久久 | 精品国产一区二区三区四区第35| 两性夫妻黄色片| 9热在线视频观看99| 国产欧美亚洲国产| 午夜久久久在线观看| 伊人久久大香线蕉亚洲五| 中国美女看黄片| 一区二区三区国产精品乱码| 国产xxxxx性猛交| 中文字幕人妻熟女乱码| 亚洲专区国产一区二区| 国产av精品麻豆| 看免费av毛片| av欧美777| 操美女的视频在线观看| 久久国产精品影院| 极品教师在线免费播放| 成人18禁在线播放| 波多野结衣av一区二区av| 亚洲人成电影免费在线| 侵犯人妻中文字幕一二三四区| 一个人免费在线观看的高清视频| 国产精品欧美亚洲77777| 久久人人97超碰香蕉20202| 欧美日韩成人在线一区二区| 欧美性长视频在线观看| 国产成人精品在线电影| 久久国产精品人妻蜜桃| 一区二区日韩欧美中文字幕| 日韩大片免费观看网站| 纯流量卡能插随身wifi吗| 超色免费av| 欧美日韩精品网址| 免费女性裸体啪啪无遮挡网站| 亚洲 欧美一区二区三区| 欧美午夜高清在线| 最新在线观看一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 桃红色精品国产亚洲av| 九色亚洲精品在线播放| 亚洲国产精品一区二区三区在线| 丝袜美足系列| 中文字幕人妻丝袜一区二区| 最近最新中文字幕大全电影3 | 国产欧美日韩综合在线一区二区| 一二三四在线观看免费中文在| 黄频高清免费视频| 午夜日韩欧美国产| 久久久水蜜桃国产精品网| 久久影院123| 一区福利在线观看| 国产亚洲欧美精品永久| 国产精品影院久久| 操出白浆在线播放| 婷婷成人精品国产| 国产精品偷伦视频观看了| 午夜激情av网站| 久久人妻福利社区极品人妻图片| 男男h啪啪无遮挡| 精品少妇黑人巨大在线播放| 欧美精品av麻豆av| 日本欧美视频一区| 可以免费在线观看a视频的电影网站| 日韩免费av在线播放| 免费观看a级毛片全部| 建设人人有责人人尽责人人享有的| 亚洲精品中文字幕一二三四区 | 午夜福利在线免费观看网站| 啦啦啦中文免费视频观看日本| 99精品欧美一区二区三区四区| 国产高清videossex| 国产高清激情床上av| 777米奇影视久久| 丝袜喷水一区| 精品久久蜜臀av无| 日韩欧美国产一区二区入口| 又大又爽又粗| 一边摸一边抽搐一进一出视频| 丝袜美腿诱惑在线| 热re99久久国产66热| 99国产精品免费福利视频| 欧美精品av麻豆av| 久久久欧美国产精品| √禁漫天堂资源中文www| avwww免费| 亚洲欧洲日产国产| 日韩 欧美 亚洲 中文字幕| 2018国产大陆天天弄谢| 国产成人系列免费观看| 在线av久久热| 国产又爽黄色视频| 成人三级做爰电影| 国产精品香港三级国产av潘金莲| 欧美一级毛片孕妇| 日本欧美视频一区| 18禁裸乳无遮挡动漫免费视频| 日韩成人在线观看一区二区三区| 亚洲免费av在线视频| 免费看a级黄色片| 色婷婷av一区二区三区视频| 精品人妻1区二区| e午夜精品久久久久久久| 精品亚洲成a人片在线观看| 亚洲成a人片在线一区二区| 国产日韩欧美亚洲二区| 女警被强在线播放| 纯流量卡能插随身wifi吗| 一进一出抽搐动态| 丝瓜视频免费看黄片| 国产精品秋霞免费鲁丝片| 国产高清视频在线播放一区| av福利片在线| 黄色 视频免费看| 国产黄频视频在线观看| 动漫黄色视频在线观看| 亚洲熟女精品中文字幕| 亚洲国产成人一精品久久久| 日韩有码中文字幕| 欧美日韩亚洲高清精品| 别揉我奶头~嗯~啊~动态视频| 久久久久久人人人人人| 最黄视频免费看| 亚洲国产成人一精品久久久| 国产成人av教育| 91成人精品电影| av又黄又爽大尺度在线免费看| 久久国产精品男人的天堂亚洲| 777米奇影视久久| 在线永久观看黄色视频| 国产精品一区二区在线不卡| 久久天堂一区二区三区四区| 搡老乐熟女国产| 国产一卡二卡三卡精品| 国产高清videossex| 男女床上黄色一级片免费看| 国产aⅴ精品一区二区三区波| 天堂动漫精品| 亚洲中文日韩欧美视频| 免费观看人在逋| 欧美日韩精品网址| 亚洲人成电影观看| 亚洲精品中文字幕在线视频| 99精品久久久久人妻精品| 亚洲专区中文字幕在线| 操美女的视频在线观看| 亚洲国产看品久久| 中文字幕最新亚洲高清| 精品少妇黑人巨大在线播放| 免费在线观看日本一区| 日韩熟女老妇一区二区性免费视频| 在线观看免费高清a一片| 成在线人永久免费视频| 欧美性长视频在线观看| 男女边摸边吃奶| 免费在线观看影片大全网站| 日本vs欧美在线观看视频|