• <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葉方慧熊國(guó)江楊曌昕徐波
    現(xiàn)代電子技術(shù) 2024年18期

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

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

    中圖分類號(hào): TN322.8+?34; TM73 " " " " " " " " " 文獻(xiàn)標(biāo)識(shí)碼: A " " " " " " " " " 文章編號(hào): 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 "引 "言

    近年來(lái),隨著新能源和儲(chǔ)能技術(shù)[1]的不斷發(fā)展,用戶對(duì)直流負(fù)荷的需求日益增長(zhǎng),傳統(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)通過(guò)電壓源型換流器(Voltage Source Converter, VSC)[4]將交流側(cè)與直流側(cè)分隔,并將多個(gè)交直流轉(zhuǎn)換環(huán)節(jié)組合起來(lái)。交直流互聯(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算法對(duì)安全域邊界點(diǎn)進(jìn)行求解。在交直流配電臺(tái)區(qū),文獻(xiàn)[8]中設(shè)計(jì)一種考慮源荷不確定性的日前調(diào)度方案,平抑系統(tǒng)的電壓波動(dòng),減少了網(wǎng)絡(luò)損耗。文獻(xiàn)[9]中提出一個(gè)基于動(dòng)態(tài)網(wǎng)絡(luò)重構(gòu)以及多能源協(xié)同的交直流配電網(wǎng)優(yōu)化調(diào)度方案,但沒(méi)有考慮換流站的損耗成本。在考慮電價(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)方法,但未考慮儲(chǔ)能裝置的接入及其損耗成本。文獻(xiàn)[13]中提出了考慮多種分布式電源的柔性互聯(lián)配電系統(tǒng)多區(qū)域、多目標(biāo)分布式優(yōu)化模型,但未考慮柔性互聯(lián)開關(guān)的傳輸功率損耗。文獻(xiàn)[14]中提出了一個(gè)基于源荷不確定性的交直流混合主動(dòng)配電網(wǎng)的分層?分布式優(yōu)化調(diào)度方案,有效提高了配電網(wǎng)的系統(tǒng)可控性以及對(duì)分布式能源的消納能力,但沒(méi)有考慮交直流雙向換流器的功率損耗。在配電網(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)化算法。與其他智能算法相比,該算法具有原理簡(jiǎn)單、設(shè)置參數(shù)少、尋優(yōu)能力強(qiáng)的優(yōu)點(diǎn)。

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

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

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

    在實(shí)現(xiàn)系統(tǒng)功率平衡的過(guò)程中,儲(chǔ)能和換流器有著重要的作用,所以在建立交直流互聯(lián)配電網(wǎng)經(jīng)濟(jì)調(diào)度模型時(shí),應(yīng)考慮儲(chǔ)能系統(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)出力大小與太陽(yáng)輻射強(qiáng)度、環(huán)境溫度、電池表面溫度以及相對(duì)光學(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]為太陽(yáng)輻射強(qiáng)度;[TC]為蓄電池發(fā)電時(shí)其表面的實(shí)際溫度,取值為[TC]=(48±2)℃;[ε]為功率溫度系數(shù),取值為[ε=-0.004 2 ℃]。

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

    交直流混合配電網(wǎng)中鉛酸蓄電池作為儲(chǔ)能裝置被廣泛使用,其具有“低儲(chǔ)高發(fā)”的作用,可提高新能源的消納率,抑制電網(wǎng)波動(dòng)。儲(chǔ)能裝置運(yùn)行過(guò)程中的荷電狀態(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í)段儲(chǔ)能蓄電池充電、放電功率;[Δt]為一個(gè)調(diào)度周期的時(shí)間步長(zhǎng),其中自放電功率一般很小,忽略不計(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]表示儲(chǔ)能系統(tǒng)運(yùn)行成本;[f3]表示交直流雙向換流器損耗成本;[f4]表示配電網(wǎng)和大電網(wǎng)之間購(gòu)售電成本。

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

    交直流互聯(lián)配電網(wǎng)運(yùn)行在不同狀態(tài)時(shí),由于設(shè)備的損耗或年限增加,會(huì)出現(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) 儲(chǔ)能系統(tǒng)運(yùn)行成本

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

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

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

    VSC在運(yùn)行過(guò)程中,傳輸功率在交流側(cè)和直流側(cè)之間雙向流動(dòng),會(huì)存在一定功率損耗。在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í)換流功率取絕對(duì)值。

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

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

    式中[Pgrid(t)]表示t時(shí)段配電網(wǎng)向大電網(wǎng)購(gòu)售電的功率。當(dāng)配電網(wǎng)向大電網(wǎng)購(gòu)電時(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)購(gòu)電時(shí),選擇分時(shí)購(gòu)電價(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在功率交互的過(guò)程中必須滿足其最大交互功率約束。

    [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) 儲(chǔ)能裝置約束

    [μ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í)段儲(chǔ)能裝置的充放電狀態(tài);SOC(t)為t時(shí)段的荷電狀態(tài);[SOCmin]表示儲(chǔ)能蓄電池最小荷電狀態(tài);[SOCmax]表示儲(chǔ)能蓄電池最大荷電狀態(tài);儲(chǔ)能電池在一次充放電的起始時(shí)刻和終止時(shí)刻的荷電狀態(tài)相等,即[SOC0=SOCend];[Pmaxc]、[Pmaxd]分別為儲(chǔ)能電池充放電的最大功率。

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

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

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

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

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

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

    假設(shè)鯨魚捕捉的獵物為最優(yōu)解,鯨魚所在位置為潛在解,在每次迭代的過(guò)程中,每頭鯨魚的位置更新策略通過(guò)隨機(jī)數(shù)p的取值和系數(shù)向量A的模來(lái)決定,隨著迭代次數(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為限定對(duì)數(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ì)算而來(lái)。

    [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算法對(duì)WOA算法進(jìn)行改進(jìn),即從限定對(duì)數(shù)螺旋形狀的常數(shù)考慮,將b值設(shè)成動(dòng)態(tài)值,并隨著迭代次數(shù)更新,隨機(jī)地改變螺旋路徑。在迭代開始時(shí)期,螺旋形狀偏大,可以增大鯨魚的搜索空間,使遍歷種群中的個(gè)體增多,增強(qiáng)算法的全局尋優(yōu)能力;在迭代后期,隨著螺線形狀逐漸減小,鯨魚會(huì)在最優(yōu)個(gè)體的附近進(jìn)行搜索,從而加快算法的速度,同時(shí)提高算法的尋優(yōu)精度。對(duì)限定對(duì)數(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)作為研究對(duì)象進(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)的儲(chǔ)能電池裝置容量為300 kW·h,充電功率最大值為30 kW,放電功率最大值為45 kW;VSC換流器傳輸功率的最大值為100 kW,換流效率為0.95;交直流互聯(lián)配電網(wǎng)與大電網(wǎng)交互功率的最大值為100 kW。儲(chǔ)能裝置的單位損耗成本為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é)果分析

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

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

    通過(guò)EWOA算法仿真得出各單元出力的調(diào)度策略,對(duì)調(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í)儲(chǔ)能系統(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)購(gòu)電,儲(chǔ)能系統(tǒng)放電供用戶使用。13:00—16:00負(fù)荷量減少,輻射強(qiáng)度較大,此時(shí)段燃?xì)廨啓C(jī)發(fā)電量、光伏和風(fēng)電發(fā)電量大于負(fù)荷量,配電網(wǎng)向大電網(wǎng)售電,儲(chǔ)能系統(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)購(gòu)電,儲(chǔ)能系統(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)售電,儲(chǔ)能裝置進(jìn)行充電。由此可得:儲(chǔ)能系統(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ā)電及儲(chǔ)能放電功率與直流負(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ā)電和儲(chǔ)能放電功率之和大于直流負(fù)荷,反之則小于直流負(fù)荷。交直流側(cè)功率流向關(guān)系如圖6所示。

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

    5.2.2 "算法對(duì)比分析

    EWOA算法與WOA算法、TLBO算法的優(yōu)化調(diào)度結(jié)果對(duì)比如表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ù)成本、儲(chǔ)能系統(tǒng)運(yùn)行成本、交直流雙向換流器損耗成本以及配電網(wǎng)與大電網(wǎng)之間購(gòu)售電成本之和最小的經(jīng)濟(jì)調(diào)度模型,在滿足配電網(wǎng)安全穩(wěn)定運(yùn)行約束的條件下進(jìn)行優(yōu)化調(diào)度。文中提出一種有效求解模型的EWOA算法,將對(duì)數(shù)螺旋形狀的常數(shù)b進(jìn)行動(dòng)態(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)化性能。

    注:本文通訊作者為熊國(guó)江。

    參考文獻(xiàn)

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

    [2] 魏斌,韓肖清,李雯,等.融合多場(chǎng)景分析的交直流混合微電網(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與光?儲(chǔ)?充協(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] 孫碩愷,張妍婷,戚星宇.考慮源荷不確定性的交直流配電臺(tái)區(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)自動(dòng)化,2021,45(8):95?103.

    [14] 梁海平,王巖,劉英培,等.計(jì)及源荷不確定性的混合交直流主動(dòng)配電網(wǎng)分層?分布式優(yōu)化調(diào)度[J].電力自動(dòng)化設(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)無(wú)功調(diào)度和儲(chǔ)能優(yōu)化方法[J].電力科學(xué)與技術(shù)學(xué)報(bào),2022,37(4):134?142.

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

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

    [19] 鄧長(zhǎng)征,馮朕,邱立,等.基于混沌灰狼算法的交直流混合微網(wǎng)經(jīng)濟(jì)調(diào)度[J].電測(cè)與儀表,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ì)及儲(chǔ)能損耗和換流成本的交直流混合微網(wǎng)區(qū)域協(xié)調(diào)調(diào)度[J].電力系統(tǒng)自動(dòng)化,2017,41(23):29?37.

    [22] 張?chǎng)赫鸩?,郭毅,?含混合儲(chǔ)能的交直流配電網(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.

    亚洲经典国产精华液单| 亚洲va在线va天堂va国产| 日韩欧美一区二区三区在线观看| 久久久欧美国产精品| 免费看av在线观看网站| 狠狠狠狠99中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 免费在线观看影片大全网站| 色尼玛亚洲综合影院| 国产亚洲av嫩草精品影院| 日日撸夜夜添| 亚洲一区二区三区色噜噜| 色5月婷婷丁香| 大又大粗又爽又黄少妇毛片口| 干丝袜人妻中文字幕| 男女视频在线观看网站免费| 亚洲综合色惰| 我的老师免费观看完整版| 少妇熟女aⅴ在线视频| 欧美另类亚洲清纯唯美| 精品久久久久久久末码| 长腿黑丝高跟| 麻豆久久精品国产亚洲av| 国产免费男女视频| 欧美成人精品欧美一级黄| 国模一区二区三区四区视频| 一进一出好大好爽视频| av.在线天堂| 欧美人与善性xxx| 青春草视频在线免费观看| 国产精品久久久久久久久免| 18禁裸乳无遮挡免费网站照片| 最近视频中文字幕2019在线8| 国产伦精品一区二区三区视频9| 色哟哟哟哟哟哟| 97超级碰碰碰精品色视频在线观看| 亚洲最大成人av| 日本色播在线视频| 国产精品亚洲一级av第二区| 国内精品久久久久精免费| 非洲黑人性xxxx精品又粗又长| 一个人免费在线观看电影| 免费av毛片视频| 三级国产精品欧美在线观看| 日韩,欧美,国产一区二区三区 | 日本一本二区三区精品| 久久久久性生活片| 卡戴珊不雅视频在线播放| 天堂动漫精品| 美女大奶头视频| 精品熟女少妇av免费看| 国产探花极品一区二区| 国产高清三级在线| 色av中文字幕| 久久国内精品自在自线图片| 91在线精品国自产拍蜜月| 有码 亚洲区| 黑人高潮一二区| 亚洲国产欧美人成| 丰满的人妻完整版| 日韩一本色道免费dvd| 成年版毛片免费区| 欧美潮喷喷水| 免费黄网站久久成人精品| 色视频www国产| 一边摸一边抽搐一进一小说| 亚洲性夜色夜夜综合| 神马国产精品三级电影在线观看| 久久6这里有精品| 狂野欧美白嫩少妇大欣赏| 麻豆av噜噜一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 如何舔出高潮| 免费看av在线观看网站| 狂野欧美白嫩少妇大欣赏| 免费观看的影片在线观看| 久久午夜福利片| 成人午夜高清在线视频| 91在线精品国自产拍蜜月| 校园人妻丝袜中文字幕| 有码 亚洲区| 精品久久久久久久久av| 成年女人永久免费观看视频| 亚洲精品色激情综合| 日韩三级伦理在线观看| 少妇丰满av| 亚洲成人中文字幕在线播放| 99国产极品粉嫩在线观看| 精品人妻偷拍中文字幕| 国产熟女欧美一区二区| 少妇高潮的动态图| 尾随美女入室| 久久亚洲国产成人精品v| 国产在线男女| 少妇的逼水好多| 国产大屁股一区二区在线视频| 欧洲精品卡2卡3卡4卡5卡区| 精品日产1卡2卡| 国产亚洲欧美98| 成人鲁丝片一二三区免费| 在现免费观看毛片| 毛片女人毛片| 亚洲婷婷狠狠爱综合网| 日韩av在线大香蕉| 91在线精品国自产拍蜜月| 日韩欧美免费精品| 亚洲熟妇中文字幕五十中出| 久久这里只有精品中国| 久久99热这里只有精品18| 男女做爰动态图高潮gif福利片| 国产精品福利在线免费观看| 永久网站在线| 搞女人的毛片| 丝袜喷水一区| 国产久久久一区二区三区| 国产精品一区www在线观看| 日本一二三区视频观看| 精品免费久久久久久久清纯| 国产一区二区在线av高清观看| 色吧在线观看| 日韩欧美精品v在线| 国产极品精品免费视频能看的| av福利片在线观看| 在线天堂最新版资源| 99热这里只有精品一区| 国产日本99.免费观看| 精品欧美国产一区二区三| 亚洲五月天丁香| 亚洲aⅴ乱码一区二区在线播放| 欧美一区二区国产精品久久精品| 久久人人爽人人爽人人片va| 久久久久久伊人网av| 欧美日韩在线观看h| 99国产极品粉嫩在线观看| 欧美三级亚洲精品| 国产高清不卡午夜福利| 免费无遮挡裸体视频| 搞女人的毛片| 一级毛片aaaaaa免费看小| 久久久精品大字幕| 免费黄网站久久成人精品| 老司机福利观看| 久久久国产成人精品二区| 久久久色成人| 国内精品宾馆在线| 毛片女人毛片| 国产真实乱freesex| 成人av在线播放网站| 久99久视频精品免费| 久99久视频精品免费| aaaaa片日本免费| 丝袜喷水一区| 中文字幕av在线有码专区| 国产精品99久久久久久久久| 看十八女毛片水多多多| 国产高潮美女av| 看十八女毛片水多多多| 少妇高潮的动态图| 成人亚洲欧美一区二区av| 一级黄色大片毛片| 日本黄色视频三级网站网址| 麻豆久久精品国产亚洲av| 日本免费一区二区三区高清不卡| 12—13女人毛片做爰片一| 国产免费一级a男人的天堂| 久久久久久九九精品二区国产| 国产三级中文精品| 日韩一本色道免费dvd| 一级a爱片免费观看的视频| 久久国产乱子免费精品| 看非洲黑人一级黄片| 免费av毛片视频| 久久精品影院6| a级毛片a级免费在线| 噜噜噜噜噜久久久久久91| 免费观看人在逋| av在线播放精品| 午夜老司机福利剧场| 麻豆成人午夜福利视频| 最近2019中文字幕mv第一页| 久久精品国产99精品国产亚洲性色| 色尼玛亚洲综合影院| 级片在线观看| 国产探花在线观看一区二区| a级一级毛片免费在线观看| 国产av麻豆久久久久久久| 成人永久免费在线观看视频| 人人妻人人澡人人爽人人夜夜 | 欧美性猛交╳xxx乱大交人| 久久韩国三级中文字幕| 国内精品宾馆在线| 狂野欧美白嫩少妇大欣赏| 午夜激情福利司机影院| 亚洲人成网站在线播| 欧美最黄视频在线播放免费| 欧美日韩国产亚洲二区| 欧美日韩综合久久久久久| 国产精品一区二区性色av| ponron亚洲| 自拍偷自拍亚洲精品老妇| avwww免费| 丰满人妻一区二区三区视频av| 亚洲精品国产成人久久av| 亚洲七黄色美女视频| 嫩草影院入口| 欧美xxxx性猛交bbbb| 午夜福利成人在线免费观看| av黄色大香蕉| 久久亚洲国产成人精品v| av福利片在线观看| 欧美一区二区精品小视频在线| 日韩中字成人| 久久亚洲精品不卡| 精品久久久久久久久av| 我要搜黄色片| 一本久久中文字幕| 久久久国产成人精品二区| 国产又黄又爽又无遮挡在线| 日韩在线高清观看一区二区三区| 久久久久免费精品人妻一区二区| 高清午夜精品一区二区三区 | 国产伦精品一区二区三区四那| 欧美日韩国产亚洲二区| 老熟妇乱子伦视频在线观看| 一级毛片久久久久久久久女| videossex国产| 一进一出抽搐动态| 晚上一个人看的免费电影| 国产伦精品一区二区三区四那| 老熟妇乱子伦视频在线观看| 全区人妻精品视频| 色哟哟哟哟哟哟| ponron亚洲| 亚洲精品成人久久久久久| 在线免费观看的www视频| 丰满乱子伦码专区| 精品免费久久久久久久清纯| 性插视频无遮挡在线免费观看| 午夜爱爱视频在线播放| 成人亚洲精品av一区二区| 青春草视频在线免费观看| 亚洲国产精品合色在线| 日本成人三级电影网站| 国产av在哪里看| 亚洲欧美成人综合另类久久久 | 女人被狂操c到高潮| 老女人水多毛片| 成人永久免费在线观看视频| av天堂中文字幕网| 国产三级中文精品| 国产精品av视频在线免费观看| 露出奶头的视频| 精品人妻熟女av久视频| 国产色爽女视频免费观看| 特级一级黄色大片| 免费观看在线日韩| 久久人人精品亚洲av| 免费不卡的大黄色大毛片视频在线观看 | 色av中文字幕| 又黄又爽又刺激的免费视频.| 国产伦一二天堂av在线观看| 22中文网久久字幕| 我的女老师完整版在线观看| 日韩成人av中文字幕在线观看 | 午夜日韩欧美国产| 国产成人精品久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 久久精品综合一区二区三区| 国内精品宾馆在线| 欧美国产日韩亚洲一区| 久久精品国产亚洲av涩爱 | 国产白丝娇喘喷水9色精品| 神马国产精品三级电影在线观看| 国产成人91sexporn| 中国国产av一级| av在线天堂中文字幕| 午夜免费激情av| 亚洲va在线va天堂va国产| 精品日产1卡2卡| 亚洲欧美成人精品一区二区| 精品不卡国产一区二区三区| 成人国产麻豆网| 精品免费久久久久久久清纯| 久久精品国产亚洲网站| 精品久久久噜噜| 婷婷精品国产亚洲av| 少妇人妻一区二区三区视频| 日韩欧美免费精品| 久久人人爽人人片av| 真人做人爱边吃奶动态| 三级国产精品欧美在线观看| 久久久久久伊人网av| 国产成人a∨麻豆精品| 精品午夜福利视频在线观看一区| 国产亚洲精品av在线| 男女啪啪激烈高潮av片| 亚洲av美国av| 亚洲精品国产av成人精品 | 男女啪啪激烈高潮av片| 国产黄a三级三级三级人| 精品欧美国产一区二区三| 真实男女啪啪啪动态图| 日本成人三级电影网站| 国产高清激情床上av| 丝袜美腿在线中文| 女人被狂操c到高潮| 在线播放无遮挡| 男人舔奶头视频| 偷拍熟女少妇极品色| 狂野欧美白嫩少妇大欣赏| 人妻丰满熟妇av一区二区三区| 黄色配什么色好看| 日韩中字成人| 久久久精品94久久精品| 人妻少妇偷人精品九色| 天堂网av新在线| 色哟哟哟哟哟哟| 亚洲自偷自拍三级| 久久久色成人| 黄色欧美视频在线观看| 欧美三级亚洲精品| 久久中文看片网| 在线观看午夜福利视频| 欧美另类亚洲清纯唯美| 不卡视频在线观看欧美| 国产精品福利在线免费观看| a级毛片a级免费在线| 在现免费观看毛片| 欧美高清性xxxxhd video| 国语自产精品视频在线第100页| 国产三级在线视频| 欧美性猛交黑人性爽| 麻豆成人午夜福利视频| 亚洲成av人片在线播放无| 婷婷六月久久综合丁香| 日韩 亚洲 欧美在线| 白带黄色成豆腐渣| 日韩,欧美,国产一区二区三区 | 女人十人毛片免费观看3o分钟| 国产精品伦人一区二区| 内射极品少妇av片p| 亚洲成人精品中文字幕电影| 69人妻影院| 欧美成人精品欧美一级黄| 国产白丝娇喘喷水9色精品| 久久这里只有精品中国| 美女内射精品一级片tv| 真人做人爱边吃奶动态| 精品一区二区免费观看| 日本撒尿小便嘘嘘汇集6| 一级黄色大片毛片| 精品久久久久久久久av| 国产精品伦人一区二区| 精品一区二区三区视频在线观看免费| 波多野结衣高清作品| 女人被狂操c到高潮| 亚洲专区国产一区二区| 床上黄色一级片| 日本在线视频免费播放| 蜜臀久久99精品久久宅男| 日韩 亚洲 欧美在线| 22中文网久久字幕| 最近2019中文字幕mv第一页| 日本免费一区二区三区高清不卡| 国产黄色视频一区二区在线观看 | 99久久九九国产精品国产免费| 最近中文字幕高清免费大全6| 熟妇人妻久久中文字幕3abv| 男女之事视频高清在线观看| 国产成人freesex在线 | 在线天堂最新版资源| 99riav亚洲国产免费| 精品不卡国产一区二区三区| 毛片女人毛片| av在线蜜桃| 国产精品嫩草影院av在线观看| 在线a可以看的网站| 成人午夜高清在线视频| 国产亚洲精品av在线| 中文字幕熟女人妻在线| 内射极品少妇av片p| 日本a在线网址| 特级一级黄色大片| 久久6这里有精品| 国产精品一区二区三区四区免费观看 | 直男gayav资源| 男女做爰动态图高潮gif福利片| 欧美国产日韩亚洲一区| 丰满的人妻完整版| 国产成人a区在线观看| 精品午夜福利在线看| 99热这里只有是精品50| 女生性感内裤真人,穿戴方法视频| 亚洲av中文字字幕乱码综合| 中国国产av一级| 国产成人一区二区在线| av在线播放精品| 国产精品一区二区三区四区免费观看 | 联通29元200g的流量卡| 成年女人永久免费观看视频| 久久久久久久久中文| 最近视频中文字幕2019在线8| 欧美日本视频| 少妇裸体淫交视频免费看高清| 日韩国内少妇激情av| 天美传媒精品一区二区| 亚洲av二区三区四区| 中国美女看黄片| 国产亚洲精品久久久久久毛片| 99久久精品国产国产毛片| 寂寞人妻少妇视频99o| 老司机影院成人| 午夜福利高清视频| 99热全是精品| 3wmmmm亚洲av在线观看| 久久久久久久久久久丰满| 香蕉av资源在线| 久久久久久久亚洲中文字幕| 国产精品三级大全| 色5月婷婷丁香| 男女之事视频高清在线观看| 亚洲av二区三区四区| 久久久久九九精品影院| 成人av一区二区三区在线看| 欧美色欧美亚洲另类二区| 亚洲成人精品中文字幕电影| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美三级三区| 亚洲av中文字字幕乱码综合| 一级毛片久久久久久久久女| 一个人看的www免费观看视频| 最后的刺客免费高清国语| 亚洲中文字幕日韩| 亚洲人与动物交配视频| 亚洲精品在线观看二区| 亚洲欧美精品综合久久99| 成人国产麻豆网| 又爽又黄a免费视频| 国产高清视频在线观看网站| 人人妻人人看人人澡| 日本免费a在线| 国产白丝娇喘喷水9色精品| 午夜老司机福利剧场| 人妻夜夜爽99麻豆av| 国产美女午夜福利| 日本色播在线视频| 尤物成人国产欧美一区二区三区| 九色成人免费人妻av| 免费av观看视频| 亚洲七黄色美女视频| 亚洲av二区三区四区| 久久九九热精品免费| 日本免费一区二区三区高清不卡| 干丝袜人妻中文字幕| 欧美日韩精品成人综合77777| 成人亚洲精品av一区二区| 久久亚洲精品不卡| 久久久久久久亚洲中文字幕| 精品久久国产蜜桃| 色5月婷婷丁香| 日本a在线网址| 大香蕉久久网| 18禁裸乳无遮挡免费网站照片| 在线天堂最新版资源| 国内精品宾馆在线| 在线观看av片永久免费下载| 国产亚洲欧美98| 美女免费视频网站| 欧美国产日韩亚洲一区| 美女黄网站色视频| 欧美不卡视频在线免费观看| 亚洲av成人av| 国产精品不卡视频一区二区| 国产高清有码在线观看视频| 欧美区成人在线视频| 99久国产av精品国产电影| 91午夜精品亚洲一区二区三区| 97超级碰碰碰精品色视频在线观看| 搡老岳熟女国产| 国产精品久久电影中文字幕| 在线免费观看的www视频| 免费人成在线观看视频色| 国产一区二区三区在线臀色熟女| 免费看光身美女| 一边摸一边抽搐一进一小说| 我的老师免费观看完整版| 精品无人区乱码1区二区| 村上凉子中文字幕在线| 欧美高清性xxxxhd video| 99热只有精品国产| 欧美潮喷喷水| 日韩成人av中文字幕在线观看 | 性插视频无遮挡在线免费观看| 国产老妇女一区| 国产欧美日韩一区二区精品| 亚洲av中文av极速乱| 干丝袜人妻中文字幕| 男女下面进入的视频免费午夜| av黄色大香蕉| 床上黄色一级片| 97超碰精品成人国产| 少妇猛男粗大的猛烈进出视频 | 天天一区二区日本电影三级| 一个人观看的视频www高清免费观看| 免费观看在线日韩| 在线免费观看不下载黄p国产| 国产爱豆传媒在线观看| 好男人在线观看高清免费视频| 91久久精品电影网| 欧美中文日本在线观看视频| 日韩欧美一区二区三区在线观看| 亚洲av成人精品一区久久| 亚洲国产欧美人成| 国产精品乱码一区二三区的特点| 一边摸一边抽搐一进一小说| 欧美一区二区精品小视频在线| 免费看av在线观看网站| 丰满的人妻完整版| 99久久精品一区二区三区| 亚洲av一区综合| av天堂中文字幕网| 麻豆成人午夜福利视频| 亚洲成人精品中文字幕电影| av黄色大香蕉| 欧美一区二区国产精品久久精品| 黄色配什么色好看| 国产欧美日韩一区二区精品| 色在线成人网| 级片在线观看| 久久久久久久午夜电影| 日韩国内少妇激情av| 一边摸一边抽搐一进一小说| 日韩人妻高清精品专区| 极品教师在线视频| 日韩三级伦理在线观看| 国产精品野战在线观看| 国产黄色小视频在线观看| 变态另类成人亚洲欧美熟女| 欧美xxxx性猛交bbbb| 国产精品av视频在线免费观看| 99久久精品一区二区三区| 亚洲国产精品成人综合色| 久久久久久九九精品二区国产| 老司机福利观看| 精品久久久久久久久av| 久久鲁丝午夜福利片| 麻豆成人午夜福利视频| 色噜噜av男人的天堂激情| 亚洲丝袜综合中文字幕| 欧美激情在线99| 国产高清视频在线播放一区| 午夜精品国产一区二区电影 | 久久精品综合一区二区三区| 看十八女毛片水多多多| 熟女人妻精品中文字幕| 搡女人真爽免费视频火全软件 | av免费在线看不卡| 色5月婷婷丁香| 久久久久国内视频| 久久人人精品亚洲av| 国产aⅴ精品一区二区三区波| 深爱激情五月婷婷| 日本黄色视频三级网站网址| 色视频www国产| 精品午夜福利在线看| 精品少妇黑人巨大在线播放 | 欧美zozozo另类| 国产高清视频在线播放一区| 精品99又大又爽又粗少妇毛片| 精品久久国产蜜桃| 人人妻人人澡人人爽人人夜夜 | 欧美日韩一区二区视频在线观看视频在线 | 在线观看午夜福利视频| 免费看a级黄色片| 国产在线精品亚洲第一网站| 亚洲七黄色美女视频| 国内少妇人妻偷人精品xxx网站| 九九爱精品视频在线观看| 一级毛片电影观看 | 亚洲国产精品久久男人天堂| a级毛片a级免费在线| 国产熟女欧美一区二区| 成人二区视频| 99在线人妻在线中文字幕| 国产黄片美女视频| 久久精品国产亚洲网站| 深夜精品福利| 欧美激情久久久久久爽电影| 狂野欧美白嫩少妇大欣赏| 香蕉av资源在线| 搞女人的毛片| 午夜爱爱视频在线播放| av天堂在线播放| 亚洲精品乱码久久久v下载方式| 日本精品一区二区三区蜜桃| 少妇高潮的动态图| 久久人人爽人人片av| 日韩欧美三级三区| 黄色配什么色好看| 秋霞在线观看毛片| 久久久久久久久久黄片| 国内精品一区二区在线观看| 国产麻豆成人av免费视频| 免费电影在线观看免费观看| 伦精品一区二区三区| 久久久国产成人免费| 国产精品av视频在线免费观看| 久久亚洲国产成人精品v| 男人和女人高潮做爰伦理| 精华霜和精华液先用哪个| 亚洲人成网站在线观看播放| 久久久久九九精品影院| 3wmmmm亚洲av在线观看| 看黄色毛片网站| 一本一本综合久久|