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

    基于改進(jìn)蜉蝣算法的光伏多峰值最大功率跟蹤特性研究

    2024-02-13 00:00:00趙智勇李衛(wèi)軍鄭秒朱運(yùn)奇張翼
    太陽(yáng)能學(xué)報(bào) 2024年12期
    關(guān)鍵詞:光伏發(fā)電

    摘 要:使用一種改進(jìn)蜉蝣算法(EMA)解決在局部陰影條件下,太陽(yáng)電池P-U圖呈現(xiàn)多峰值,傳統(tǒng)最大功率點(diǎn)跟蹤算法易陷入局部最優(yōu)解,導(dǎo)致無(wú)法快速準(zhǔn)確找到光伏最大功率跟蹤點(diǎn)的問(wèn)題。首先通過(guò)雄性蜉蝣的全局搜索,跳出局部最優(yōu)解;然后利用雌性蜉蝣的隨機(jī)局部搜索和交配,減少系統(tǒng)震蕩;最后通過(guò)設(shè)置兩組不同的算例,驗(yàn)證所使用算法的跟蹤精度和跟蹤速度。研究結(jié)果表明:與粒子群算法(PSO)和灰狼算法(GWO)相比,改進(jìn)蜉蝣算法在搜索過(guò)程中可減少光伏輸出功率震蕩。相較于PSO和GWO,在EWA在跟蹤速度上得到顯著提升,使得太陽(yáng)電池在局部陰影下仍能保持高效的輸出功率。

    關(guān)鍵詞:光伏發(fā)電;最大功率點(diǎn)跟蹤;改進(jìn)蜉蝣算法;局部陰影;多峰值;光伏輸出特性

    中圖分類號(hào):TM615 " " " " " " " " " " " "文獻(xiàn)標(biāo)志碼:A

    0 引 言

    據(jù)國(guó)際能源署(International Energy Agency,IEA)預(yù)測(cè),到2050年,風(fēng)能和太陽(yáng)能發(fā)電量將占據(jù)全球總發(fā)電量的近70%[1]。目前,隨著太陽(yáng)能發(fā)電技術(shù)的不斷發(fā)展及成本的下降,研究人員已將注意力轉(zhuǎn)移到如何提升太陽(yáng)電池的轉(zhuǎn)換效率方面[2]。太陽(yáng)電池須根據(jù)環(huán)境變化實(shí)時(shí)跟蹤其最大功率的電壓值,傳統(tǒng)最大功率點(diǎn)跟蹤(maximum power point tracking,MPPT)控制方法有擾動(dòng)觀察法(perturbation and observation,Pamp;O)、電導(dǎo)增量法(incremental conductance,INC)等[3]。同時(shí),現(xiàn)代智能算法如粒子群算法(particle swarm optimization,PSO)、遺傳算法(genetic algorithm,GA)等也得到較為廣泛的應(yīng)用。文獻(xiàn)[4]針對(duì)預(yù)先給定步長(zhǎng)因子的光伏單峰MPPT方法無(wú)法保證在多種運(yùn)行條件下的適應(yīng)性問(wèn)題,提出一種基于序列二次規(guī)劃(sequential quadratic programming,SQP)信賴域算法的光伏單峰MPPT方法,實(shí)現(xiàn)了光伏單峰MPPT控制;文獻(xiàn)[5]針對(duì)最大功率點(diǎn)搜索精度不夠高的問(wèn)題,提出一種基于最優(yōu)梯度法的MPPT全數(shù)字控制方法,實(shí)現(xiàn)了對(duì)最大功率點(diǎn)跟蹤的高精度控制;文獻(xiàn)[6]針對(duì)太陽(yáng)電池的輸出性能較低的問(wèn)題,提出一種基于擬合曲線法和最優(yōu)梯度法相結(jié)合的MPPT算法,有效提高了太陽(yáng)電池的光電轉(zhuǎn)換效率。

    上述研究工作部分解決了現(xiàn)代智能算法應(yīng)用于單峰值尋優(yōu)時(shí)的響應(yīng)速度和響應(yīng)精度不夠的問(wèn)題。工程實(shí)際中,往往會(huì)遇到如落葉、灰塵、鳥糞、云朵等局部遮擋的工況。導(dǎo)致太陽(yáng)電池輸出功率與電壓P-U曲線出現(xiàn)多峰值衰退的現(xiàn)象,傳統(tǒng)算法會(huì)陷入局部最優(yōu)解從而降低了光伏發(fā)電效率。現(xiàn)代智能算法應(yīng)用于多峰值的問(wèn)題,如文獻(xiàn)[7]針對(duì)傳統(tǒng)MPPT控制方法易陷入局部最優(yōu)解的問(wèn)題,提出一種基于卡爾曼濾波和改進(jìn)粒子群的MPPT算法,有效提高了光伏系統(tǒng)的發(fā)電效率;文獻(xiàn)[8]針對(duì)光電轉(zhuǎn)換效率較低的問(wèn)題,提出一種基于擾動(dòng)觀察法的光伏MPPT技術(shù),實(shí)現(xiàn)了最大功率點(diǎn)的追蹤;文獻(xiàn)[9]針對(duì)光伏系統(tǒng)發(fā)電穩(wěn)定度和發(fā)電效率較低的問(wèn)題,提出一種基于回歸算法與粒子群算法相結(jié)合的最大功率跟蹤方法,提高了對(duì)最大功率點(diǎn)的跟蹤精度。

    針對(duì)局部遮陰太陽(yáng)電池多峰值問(wèn)題,且為克服傳統(tǒng)智能優(yōu)化算法在跟蹤光伏發(fā)電系統(tǒng)最大功率點(diǎn)時(shí)易陷入局部最優(yōu)的難題,本文使用一種改進(jìn)的蜉蝣算法(enhanced mayfly algorithm,EMA),其通過(guò)雌性的尋偶特性,跳出局部最優(yōu),并在蜉蝣算法中加入帕累托(Pareto)算法增強(qiáng)其占空比的搜索能力,使其在運(yùn)行過(guò)程中,具有更好的穩(wěn)定性和響應(yīng)速度。

    1 局部遮陰太陽(yáng)電池輸出特性

    當(dāng)光伏組件受到局部陰影遮擋時(shí),其輸出特性曲線將從單峰值變?yōu)槎喾逯担ぷ髟诰植空陉幑r下太陽(yáng)電池電壓[10]:

    [Io=Ip-IsexpαUo-IoReqgkT-1-Uo+IoReqRb] (1)

    式中:[Io]——輸出電流,A;[Ip]——光生電流,A;[Is]——反向飽和電流,A;α——二極管理想因子;[Uo]——輸出電壓,V;[Req]——太陽(yáng)電池等效串聯(lián)電阻,Ω;[q]——電子電荷,[l=1.602×10-19 C];[k]——玻爾茲曼常數(shù),1.381×[10-23] J/K;[T]——電池溫度,℃;[Rb]——太陽(yáng)電池等效并聯(lián)電阻,Ω。

    太陽(yáng)電池輸出特性主要受輻照強(qiáng)度和溫度的影響,具有非線性特征[15]。在局部遮陰條件下光伏陣列的U-I輸出特性曲線呈多膝現(xiàn)象,其P-U曲線呈多峰現(xiàn)象。由[N]塊太陽(yáng)電池所構(gòu)成并、串聯(lián)的光伏陣列輸出電流[I]:[11]。

    [I=Isc-N11-C1expUmC2Uoc-1Isc-N21-C1expUmC2Uoc-1?Isc-NM1-C1expUmC2Uoc-1] (2)

    式中:[Isc-NM]——對(duì)應(yīng)輻照度的短路電流,A;[Um]——最大功率點(diǎn)電壓,V;[Uoc]——開(kāi)路電壓,V。

    輻照度和環(huán)境溫度是影響太陽(yáng)電池輸出特性的主要因素。圖1所示為太陽(yáng)電池[P-U]和[U-I]輸出特性曲線。其中圖1a及圖1b的測(cè)試條件為:溫度為25 ℃、輻照度分別為900—900—900 W/m2、900—500—900 W/m2和900—500—400 W/m2;圖1c及圖1d的測(cè)試條件為:輻照度為500 W/m2,溫度分別為25、35、45 ℃。

    如圖1a所示,隨著輻照度的減小,太陽(yáng)電池伏安輸出特性曲線向下偏移,短路電流與輻照度成正比,表明短路電流隨輻照度的增加而增大;而開(kāi)路電壓在橫軸方向上的交點(diǎn)向左有較小偏移,表明開(kāi)路電壓隨輻照度的減小而減小。由圖1b可知,在跟蹤至最大功率點(diǎn)工作范圍內(nèi)時(shí),輸出電流基本保持不變,但如果超出此范圍則會(huì)快速降低。輻照度越大,最大功率點(diǎn)越高,最大功率點(diǎn)的工作范圍會(huì)越大越穩(wěn)定。由圖1c可知,隨溫度的升高,電池的伏安特性曲線向左側(cè)偏移,此現(xiàn)象表明開(kāi)路電壓與溫度成反比關(guān)系,而在縱軸方向上,曲線隨著溫度的升高向上有較小程度的偏移,說(shuō)明短路電流與溫度存在一定正比關(guān)系,可知溫度對(duì)短路電流影響較小,對(duì)開(kāi)路電壓影響較大。在橫軸方向上,曲線隨著溫度的升高向左移動(dòng),即溫度越高開(kāi)路電壓越小。由圖1d可知,在跟蹤至最大功率點(diǎn)工作范圍時(shí),溫度對(duì)最大功率點(diǎn)的影響幾乎可忽略不計(jì),但隨著溫度的升高,最大功率點(diǎn)的工作范圍越小,超出工作范圍,輸出電流會(huì)快速降低,可知光伏陣列受不同遮擋條件的影響,對(duì)應(yīng)的峰值也不同,但太陽(yáng)電池或整個(gè)陣列只有唯一的最大功率點(diǎn)。因此,跳出局部最優(yōu)解,并在變化中尋優(yōu)成為MPPT控制算法需解決的關(guān)鍵問(wèn)題。

    2 增強(qiáng)蜉蝣算法

    2.1 Mayfly模型

    蜉蝣算法(mayfly algorithm,MA)模擬蜉蝣的飛行和交配行為,并結(jié)合群體智能和迭代算法特點(diǎn),大大提高了在全局搜索中的尋優(yōu)能力[12]。MA工作原理與其生物特性一樣,通過(guò)蜉蝣在交配產(chǎn)生后代的過(guò)程,以及迭代后每個(gè)蜉蝣在搜索空間中的位置,提供了如何跳出局部最優(yōu)的解決思路。其基本原理是隨機(jī)產(chǎn)生兩組蜉蝣,分別代表雄性蜉蝣和雌性蜉蝣,并提前確定好目標(biāo)函數(shù)以便于對(duì)其性能進(jìn)行準(zhǔn)確預(yù)判,蜉蝣的新位置是由其移動(dòng)速度和當(dāng)前位置來(lái)決定。每個(gè)單獨(dú)的蜉蝣都可看作局部和全局的動(dòng)態(tài)交互,使其朝向個(gè)體最佳([pb])以及全局最佳([gb])迭代。

    2.1.1 雄性蜉蝣

    當(dāng)在一定有界區(qū)域內(nèi)投放一群雄性蜉蝣時(shí),雄性蜉蝣會(huì)成群聚集在一起,且它們的速度和位置會(huì)按一定的社會(huì)作用進(jìn)行變化,且各自的位置會(huì)根據(jù)自身以及相鄰蜉蝣個(gè)體做出調(diào)整。在時(shí)間步長(zhǎng)為[t]時(shí),其空間位置可由式(3)確定。雄性蜉蝣的速度計(jì)算[13]為:

    [xt+1i=xti+vt+1i] (3)

    [yt+1ni=vtni+A1pbni-xtni+A2gbni-xtni] (4)

    [A1=ae-βr2pA2=be-βr2g] (5)

    式中:[xt+1i]——蜉蝣在第[i]維度[t+1]時(shí)刻的位置;[xti]——蜉蝣在第[i]維度[t]時(shí)刻的位置;[vt+1ni]——蜉蝣在第[i]維度[t+1]時(shí)刻的速度;[yt+1ni]——蜉蝣[n]在第[i]維度[t+1]時(shí)刻的速度;[vtni]——蜉蝣[n]在第[i]維度[t]時(shí)刻的速度;[yni]——蜉蝣在[i]維度[t]時(shí)刻的速度;[A1]——與歷史最優(yōu)位置[pbni]相關(guān)的吸引系數(shù);[pbni]——蜉蝣[n]在第[i]維度歷史上的最優(yōu)位置;[A2]——與歷史最優(yōu)位置[gbni]相關(guān)的吸引系數(shù);[gbni]——蜉蝣[n]在第[i]維度歷史上的最優(yōu)位置;[xtni]——蜉蝣[n]在第[i]維度[t]時(shí)刻的位置;[a、b]——社會(huì)作用正吸引系數(shù);[β]——蜉蝣的能見(jiàn)度系數(shù),用來(lái)控制蜉蝣的能見(jiàn)范圍;[rp]——當(dāng)前位置與[pbni]的距離;[rg]——當(dāng)前位置與[gbni]的距離。

    2.1.2 雌性蜉蝣

    雌性蜉蝣與雄性蜉蝣相互吸引,但同性蜉蝣之間不會(huì)主動(dòng)靠近。針對(duì)這一特性,當(dāng)[ygt;x]時(shí),雌性蜉蝣的速度更新如式(6)所示;[ylt;x]時(shí),雌性蜉蝣的速度更新如式(7)所示。雌性蜉蝣的位置更新如式(8)所示[13]。

    [vt+1ni=vtni+A2e-βr2mxtni-ytni] (6)

    [vt+1ni=vtni+fl?r] (7)

    [yt+1n=ytn+vtn] (8)

    式中:[vt+1ni]——蜉蝣[n]在第[i]維度[t+1]時(shí)刻的速度;[vtni]——蜉蝣[n]在第[i]維度[t]時(shí)刻的速度;[e-βr2m]——指數(shù)衰減因子,其中[β]為影響衰減速度的參數(shù),[rm]為雄性和雌性蜉蝣之間的歐氏距離;[xtni]——第[n]個(gè)雄性蜉蝣在第[i]維度、第[t]個(gè)時(shí)間步的位置;[ytni]——第[n]個(gè)雌性蜉蝣在第[i]維度、第[t]個(gè)時(shí)間步的位置;[vt+1ni]——蜉蝣[n]在第[t]維度[t+1]時(shí)刻的速度;[fl]——隨機(jī)游走系數(shù),當(dāng)雌性蜉蝣沒(méi)有被雄性蜉蝣吸引時(shí)起作用;[r]——[r∈[-1,1]];[yt+1n]——蜉蝣[n]在[t+1]時(shí)刻的參考位置或當(dāng)前追蹤的位置,[ytn]——蜉蝣[n]在[t]時(shí)刻的參考位置或當(dāng)前追蹤的位置。

    2.1.3 蜉蝣交配

    雌雄性蜉蝣交配過(guò)程可用交叉算子表征,依據(jù)雌雄性蜉蝣相互吸引的原則,從雌雄群體中按隨機(jī)選擇或適應(yīng)度函數(shù)各選擇一個(gè)親本。在依據(jù)適應(yīng)度函數(shù)選擇方式中,最優(yōu)雌雄性個(gè)體交配,次好雌雄性個(gè)體交配。交配所產(chǎn)生的后代如式(9)所示[13]。

    [ch1=L·ma+1-Lfe, ch2=L·fe+1-L·ma] (9)

    式中:[ch1]、[ch2]——產(chǎn)生的兩個(gè)后代;[ma]——雄性親本;[fe]——雌性親本;[L]——[∈[-1,1]]。

    2.2 帕累托模型

    在實(shí)際追蹤情況下,為實(shí)現(xiàn)對(duì)占波比的快速響應(yīng)輸出,需充分考慮最大功率點(diǎn)電壓與對(duì)應(yīng)占波比的關(guān)系。在MA迭代過(guò)程中不斷更新Pareto解,在每次迭代過(guò)程中產(chǎn)生新解集,然后逐一與Pareto最優(yōu)解集進(jìn)行比較并判斷新解集是否支配Pareto最優(yōu)解集合中元素,進(jìn)而更新Pareto最優(yōu)解集合。為避免出現(xiàn)Pareto最優(yōu)解分布密集的情況,可對(duì)最優(yōu)解進(jìn)行部分剔除[14]:

    [FZxn-FZxmlt;DzDz=FZmax-FZminR] (10)

    式中:[Fzxn]——第[n]個(gè)蜉蝣個(gè)體在第[Z]個(gè)目標(biāo)函數(shù)下適應(yīng)值;[Dz]——第[Z]個(gè)目標(biāo)函數(shù)值的Pareto前沿距離閾值;[FZmax]和[FZmin]——當(dāng)前Pareto最優(yōu)解集中第[Z]個(gè)目標(biāo)函數(shù)的最大值與最小值;[R]——最優(yōu)解集集合中Pareto最優(yōu)解個(gè)數(shù)上限;[Z∈{1,2,3}]。

    2.3 改進(jìn)MA

    在一定環(huán)境變化條件下,部分遮陰的MPPT優(yōu)化模型可表示為:

    [fVoutmax=PoutVout=Vout?IoutVout] (11)

    式中:[Pout]——光伏輸出功率,kW;[Iout]——光伏輸出電流,A;[Vout]——光伏輸出電壓,V。

    蜉蝣種群初始化可表示為:

    [Xni=Vmin+φ(Vmax-Vmin)] (12)

    式中:[Vmax]、[Vmin]——光伏輸出電壓的最大值與最小值,V,其中[Vout]介于極值之間;[Xni]——第[n]個(gè)蜉蝣群中的第[i]個(gè)蜉蝣個(gè)體的位置;[φ]——位于區(qū)間[0,1]中的隨機(jī)數(shù)。

    改進(jìn)蜉蝣算法的MPPT流程圖如圖2所示。

    2.4 算法性能測(cè)試

    為驗(yàn)證EMA的尋優(yōu)性能,采用4個(gè)標(biāo)準(zhǔn)測(cè)試函數(shù)[15],如表1所列,對(duì)比分析PSO與GWO性能。在算法性能測(cè)試中,統(tǒng)一將種群數(shù)設(shè)定為30,迭代次數(shù)為200,且忽略計(jì)算機(jī)性能對(duì)算法性能測(cè)試的影響。圖3為各算法性能測(cè)試結(jié)果。表2為PSO、GWO和EMA迭代10次再取平均值的尋優(yōu)結(jié)果。結(jié)合圖3各算法性能測(cè)試結(jié)果觀察可知,EMA在單峰值測(cè)試函數(shù)尋優(yōu)中可更快搜索逼近理論最優(yōu)值。同時(shí),在多峰值測(cè)試函數(shù)尋優(yōu)中,能更快跳出局部最優(yōu)且具有更高的尋優(yōu)精度。

    3 仿真環(huán)境與模型

    本文基于Matlab/Simulink軟件平臺(tái),搭建3組件串聯(lián)光伏陣列,其中太陽(yáng)電池由北京合田方泰公司生產(chǎn),型號(hào)為HTMO-36-175[16-17],其具體參數(shù)如表3所示。仿真模型如圖4所示。設(shè)置兩個(gè)不同場(chǎng)景,場(chǎng)景1為太陽(yáng)電池部分遮擋,場(chǎng)景2為輻照度不均勻。在上述兩個(gè)場(chǎng)景中,環(huán)境溫度T均保持不變。分別利用PSO、GWO和EMA跟蹤光伏功率最大功率點(diǎn),并進(jìn)行跟蹤特性對(duì)比分析。

    4 討論與分析

    4.1 場(chǎng)景1:太陽(yáng)電池部分遮擋

    在此場(chǎng)景下,太陽(yáng)電池處于部分遮擋輻照強(qiáng)度下,其環(huán)境溫度為25 ℃,輻照度為350、450、250 W/m2。根據(jù)太陽(yáng)電池P-I曲線可知,此場(chǎng)景下,光伏發(fā)電系統(tǒng)有兩個(gè)功率峰值,其峰值均為112.2 W。圖5為通過(guò)PSO算法、GWO算法和EMA得到的太陽(yáng)電池部分遮擋條件下最大功率跟蹤曲線。由圖5可知,PSO算法跟蹤至第一峰值時(shí),陷入局部最優(yōu)解,在第一峰值附近不停波動(dòng),未能跟蹤到最大峰值的輸出功率。GWO算法最開(kāi)始上升至第一峰值后,產(chǎn)生局部擾動(dòng),搜索方向不確定,導(dǎo)致其輸出功率下降,當(dāng)輸出功率下降到接近零時(shí),重新觸發(fā)跟蹤。又因GWO算法的搜索方向錯(cuò)誤導(dǎo)致輸出功率的大幅下降,然后在反復(fù)上下擾動(dòng)的過(guò)程中定位至第二峰,但不能穩(wěn)定輸出功率。

    從搜索性能上,PSO算法能跳出局部最優(yōu)解,但因全局搜索時(shí),其輸出功率會(huì)上下劇烈波動(dòng),導(dǎo)致其輸出效率下降。EMA能在0.04 s內(nèi)完成全局的搜索并穩(wěn)定功率的輸出,且在搜索過(guò)程中未出現(xiàn)劇烈振蕩。對(duì)比傳統(tǒng)PSO和GWO算法,EMA能快速跳出局部最優(yōu)解,且能快速搜索至全局最優(yōu)處。同時(shí),由于EMA通過(guò)雌性蜉蝣代替雄性蜉蝣搜索,能減少無(wú)效搜索空域,從而穩(wěn)定輸出功率。由圖5可知,PSO、GWO、EMA收斂時(shí)間分別為0.15、0.58、0.04 s。比較可知,EMA能在更短時(shí)間內(nèi)達(dá)到最大功率點(diǎn),輸出功率更接近光伏陣列的最大功率,且在跟蹤過(guò)程中功率振蕩更小。因此,部分遮擋光照條件下,EMA在跟蹤光伏發(fā)電系統(tǒng)最大功率點(diǎn)時(shí)的性能優(yōu)于PSO和GWO。

    4.2 場(chǎng)景2:輻照強(qiáng)度不均勻

    此場(chǎng)景中,太陽(yáng)電池處于不均勻輻照度下,其中環(huán)境溫度為25 ℃,輻照度為900、500、400 W/m2。根據(jù)太陽(yáng)電池P-I曲線可知,此場(chǎng)景下有3個(gè)功率峰值,但只有一個(gè)最大功率點(diǎn),為224.4 W。圖6為通過(guò)PSO算法、GWO算法和EMA得到的太陽(yáng)電池不均勻光照下最大功率跟蹤曲線。由圖6可知,PSO算法和GWO算法當(dāng)識(shí)別到功率的異常變動(dòng)后即作為搜索的開(kāi)始,并在搜索區(qū)間不停進(jìn)行全局搜索。因?yàn)榇藞?chǎng)景下有3個(gè)功率峰值,導(dǎo)致其輸出功率頻繁波動(dòng),不能穩(wěn)定輸出。太陽(yáng)電池在不均勻光照下導(dǎo)致輸出功率波動(dòng)時(shí),EMA上下搜索的幅度較小且會(huì)及時(shí)收斂至最大輸出功率點(diǎn),從而驗(yàn)證了EMA在搜索速度上的優(yōu)越性。由圖6可知,PSO、GWO和EMA算法收斂時(shí)間分別為0.14、0.21、0.08 s,比較可知EMA可在更短時(shí)間內(nèi)跟蹤到最大功率點(diǎn),輸出功率更接近光伏陣列的最大功率且在跟蹤過(guò)程中功率振蕩較小。因此,驗(yàn)證了在不均勻輻照強(qiáng)度下,EMA對(duì)光伏發(fā)電系統(tǒng)MPPT的跟蹤性能優(yōu)于PSO算法和GWO算法。

    4.3 性能對(duì)比分析

    通過(guò)跟蹤效率進(jìn)一步評(píng)價(jià)各優(yōu)化算法MPPT模型的優(yōu)越性,跟蹤效率[η][18-19]可表示為:

    [η=1-Pa-PMPM×100%] (13)

    式中:[Pa]——標(biāo)準(zhǔn)測(cè)試條件功率,kW;[PM]——跟蹤功率,kW。

    MPPT技術(shù)的魯棒性由跟蹤時(shí)間和穩(wěn)定時(shí)間兩個(gè)參數(shù)定義?;赟I的MPPT技術(shù)有探索階段和開(kāi)發(fā)階段兩個(gè)階段,需在這兩個(gè)階段之間保持平衡。探索階段允許粒子探索搜索空間,這是技術(shù)的跟蹤時(shí)間;在開(kāi)發(fā)階段,粒子在由沉降時(shí)間定義的全局解處沉降。探索階段的波動(dòng)可在圖7中觀察到。從表4可知,任何輻照度下,改進(jìn)蜉蝣算法優(yōu)化的效率均高于PSO算法和GWO算法。因此,改進(jìn)蜉蝣算法在均勻或非均勻輻照強(qiáng)度條件下對(duì)于光伏發(fā)電系統(tǒng)最大功率點(diǎn)的跟蹤性能均優(yōu)于PSO算法和GWO算法。

    為驗(yàn)證光照遮擋條件下,EMA的光伏MPPT跟蹤效率,將光伏陣列的輻照度變化設(shè)置為:在0~0.5 s期間分別為800、800、600 W/m2;在0.5~1.0 s期間分別為900、800和900 W/m2。溫度均為25 ℃。EMA最大功率跟蹤有效性驗(yàn)證結(jié)果如圖8所示。在0.5 s時(shí),輻照度發(fā)生變化,在很短時(shí)間后,EMA再次跟蹤到最大功率點(diǎn)。因此,進(jìn)一步驗(yàn)證了EMA在輻照度發(fā)生變化的條件下,跟蹤光伏最大功率點(diǎn)的有效性。

    采用相對(duì)誤差(RE)、平均絕對(duì)誤差(MAE)和均方根誤差(RMSE)來(lái)評(píng)價(jià)各方法的性能:

    [ER=i=1n(Pa-PM)PM×100%] (14)

    [EM=i=1n(Pa-PM)n] (15)

    [ERM=i=1n(Pa-PM)2n] (16)

    式中:[ER]——相對(duì)誤差;[EM]——平均絕對(duì)誤差;[ERM]——均方根誤差;[n]——樣本數(shù)量。

    從圖9可知,EMA跟蹤精度優(yōu)于PSO算法和GWO算法,從而進(jìn)一步檢驗(yàn)了EMA在光伏發(fā)電系統(tǒng)最大功率點(diǎn)跟蹤過(guò)程中的有效性和優(yōu)越性。

    5 結(jié) 論

    本文針對(duì)串聯(lián)光伏陣列因輻照強(qiáng)度分布不均導(dǎo)致輸出P-U曲線呈現(xiàn)多峰值的現(xiàn)象,使用一種增強(qiáng)蜉蝣尋優(yōu)算法,在追蹤光伏發(fā)電系統(tǒng)最大功率點(diǎn)時(shí)可快速跳出局部最優(yōu)并穩(wěn)定在全局最優(yōu)解。通過(guò)設(shè)置不同的場(chǎng)景,對(duì)比不同算法的功率跟蹤結(jié)果,得到如下主要結(jié)論:

    1)本文在MA中加入Pareto算法增強(qiáng)其占空比搜索能力,利用同性蜉蝣位置互換原理并將突變機(jī)制引入蜉蝣交配階段,在增強(qiáng)MA中雌性蜉蝣局部搜索能力的同時(shí)也提升了MA的全局搜索性能。通過(guò)Pareto最優(yōu)解集本身的全局尋優(yōu)特性,使其在復(fù)雜環(huán)境中具有更快的收斂速度和較小的功率振蕩幅度,從而比PSO和GWO擁有更高的功率跟蹤效率。

    2)不同場(chǎng)景下,與PSO相比,EMA可有效減少無(wú)效搜索空域,從而減少振蕩。與GWO相比,EMA在跟蹤精度上限相同的情況下,其跟蹤速度得到顯著提升,證明了本文使用的EMA在光伏發(fā)電系統(tǒng)最大功率點(diǎn)跟蹤全局搜索能力上較好的跟蹤精度和跟蹤速度。

    [參考文獻(xiàn)]

    [1] 伍浩松, 戴定, 趙暢. 國(guó)際能源署發(fā)布2050年凈零排放路線圖[J]. 國(guó)外核新聞, 2021(6): 17-22.

    WU H S, DAI D, ZHAO C. The International Energy Agency released the road map of net zero emission in 2050[J]. Foreign nuclear news, 2021(6): 17-22.

    [2] KIM J, HO-BAILLIE A, HUANG S J. Review of novel passivation techniques for efficient and stable perovskite solar cells[J]. Solar RRL, 2019, 3(4): 1800302.

    [3] 劉文鋒, 李昂. 基于變步長(zhǎng)擾動(dòng)觀察法結(jié)合改進(jìn)天牛群優(yōu)化算法的三步復(fù)合光伏最大功率跟蹤算法[J]. 熱力發(fā)電, 2022, 51(10): 138-144.

    LIU W F, LI A. Three-step composite photovoltaic MPPT algorithm "based "on "IP amp; "O-IBSO[J]. "Thermal "power generation, 2022, 51(10): 138-144.

    [4] 屠卿瑞, 李一泉, 曾耿暉, 等. 一種基于SQP信賴域算法的光伏單峰MPPT方法: CN107272815B[P]. 2018-09-11.

    TU Q R, LI Y Q, ZENG G H, et al. A kind of SQP-based trust domain computing photovoltaic unimodal MPPT method: CN107272815B[P]. 2018-09-11.

    [5] 林文立, 劉治鋼, 馬亮. 基于最優(yōu)梯度法的MPPT全數(shù)字控制仿真[J]. 航天器工程, 2013, 22(4): 82-86.

    LIN W L, LIU Z G, MA L. Simulation of maximum power point tracking digital control based on optimized gradient method[J]. Spacecraft engineering, 2013, 22(4): 82-86.

    [6] 羅嘉, 郭興眾, 吳飛. 基于擬合曲線法和最優(yōu)梯度法相結(jié)合的MPPT算法[J]. 南陽(yáng)理工學(xué)院學(xué)報(bào), 2016, 8(2): 27-32.

    LUO J, GUO X Z, WU F. MPPT algorithm based on fitting curve method and optimal gradient method[J]. Journal of Nanyang Institute of Technology, 2016, 8(2): 27-32.

    [7] 潘海鵬, 李明華, 雷建峰, 等. 基于卡爾曼濾波和改進(jìn)粒子群的MPPT算法[J]. 輕工機(jī)械, 2022, 40(6): 44-51.

    PAN H P, LI M H, LEI J F, et al. MPPT algorithm based on Kalman filter and improved particle swarm optimization algorithm[J]. Light industry machinery, 2022, 40(6): 44-51.

    [8] 趙夢(mèng)玨. 基于擾動(dòng)觀察法的光伏MPPT技術(shù)的仿真分析[J]. 電工技術(shù), 2022(16): 54-56, 60.

    ZHAO M J. Analysis and simulation of photovoltaic MPPT based "on "disturbance " observation "method[J]. "Electric engineering, 2022(16): 54-56, 60.

    [9] 葉國(guó)敏, 肖文波, 吳華明. 回歸算法與粒子群算法融合算法的光伏系統(tǒng)多峰值MPPT研究[J]. 現(xiàn)代電子技術(shù), 2022, 45(15): 146-150.

    YE G M, XIAO W B, WU H M. Research on PV system multi-peak MPPT based on fusion algorithm of regression algorithm and particle swarm optimization algorithm[J]. Modern electronics technique, 2022, 45(15): 146-150.

    [10] 付文龍, 孟嘉鑫, 張赟寧, 等. 復(fù)雜遮蔭下基于改進(jìn)GWO的光伏多峰MPPT控制[J]. 太陽(yáng)能學(xué)報(bào), 2023, 44(3): 435-442.

    FU W L, MENG J X, ZHANG Y N, et al. Photovoltaic multi-peak MPPT control based on improved GWO under complex shade[J]. Acta energiae solaris sinica, 2023, 44(3): 435-442.

    [11] 韓思鵬, 蔣曉艷, 羅意, 等. 遮陰條件下光伏MPPT自適應(yīng)粒子群算法優(yōu)化[J]. 太陽(yáng)能學(xué)報(bào), 2022, 43(6): 99-105.

    HAN S P, JIANG X Y, LUO Y, et al. Photovoltaic MPPT adaptive particle swarm optimization optimization under shading " conditions[J]. "Acta "energiae "solaris "sinica, 2022, 43(6): 99-105.

    [12] 王加健, 帕孜來(lái)·馬合木提, 孔博龍. 基于改進(jìn)樽海鞘群算法的光伏系統(tǒng)MPPT研究[J]. 太陽(yáng)能學(xué)報(bào), 2022, 43(4): 191-197.

    WANG J J, PAZLAI ·MAHEMUTI, KONG B L. Research on MPPT of photovoltaic system based on improved salp swarm algorithm[J]. Acta energiae solaris sinica, 2022, 43(4): 191-197.

    [13] ZERVOUDAKIS K, TSAFARAKIS S. A mayfly optimization " "algorithm[J]. " Computers " "amp; " "industrial engineering, 2020, 145: 106559.

    [14] SUHARDI D, SYAFAAH L, IRFAN M, et al. Improvement of maximum power point tracking (MPPT) efficiency using grey wolf optimization (GWO) algorithm in photovoltaic (PV) system[J]. IOP conference series: materials science and engineering, 2019, 674(1): 012038.

    [15] 趙泰祥, 廖華, 馬遜, 等. 局部陰影下光伏組件的Matlab/Simulink仿真模擬與特性分析[J]. 太陽(yáng)能學(xué)報(bào), 2019, 40(11): 3110-3118.

    ZHAO T X, LIAO H, MA X, et al. Matlab/Simulink simulation model and characteristic analysis of PV module under partial shadow[J]. Acta energiae solaris sinica, 2019, 40(11): 3110-3118.

    [16] PERVEZ I, SHAMS I, MEKHILEF S, et al. Most valuable player algorithm based maximum power point tracking for a partially shaded PV generation system[J]. IEEE transactions on sustainable energy, 2021, 12(4): 1876-1890.

    [17] 顏景斌, 佟堯, 曹雷, 等. 局部陰影條件下光伏電池建模與輸出特性仿真[J]. 電源技術(shù), 2018, 42(5): 685-688, 692.

    YAN J B, TONG Y, CAO L, et al. Modeling and output characteristics simulation of photovoltaic cells under partial "shading "condition[J]. "Chinese " journal "of "power sources, 2018, 42(5): 685-688, 692.

    [18] 洪艷, 潘東方, 姚海峰, 等. 基于黃金分割搜索的光伏陣列最大功率點(diǎn)跟蹤方法[J]. 合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 37(4): 416-419.

    HONG Y, PAN D F, YAO H F, et al. Method of photovoltaic array maximum power point tracking based on golden section search[J]. Journal of Hefei University of Technology (natural science), 2014, 37(4): 416-419.

    [19] MO S X, YE Q T, JIANG K P, et al. An improved MPPT method for photovoltaic systems based on mayfly optimization "algorithm[J]. "Energy "reports, "2022, "8: 141-150.

    RESEARCH ON PHOTOVOLTAIC MULTI-PEAK MAXIMUM POWER TRACKING CHARACTERISTICS BASED ON IMPROVED

    MAYFLY ALGORITHM

    Zhao Zhiyong,Li Weijun,Zheng Miao,Zhu Yunqi,Zhang Yi

    (State Grid Henan Electric Power Company Zhumadian Power Supply Company, Zhumadian 463000, China)

    Abstract:An improved mayfly algorithm (EMA) is used to solve the problem of multiple peaks in solar ceu P-U graphs under partial shading conditions, where traditional maximum power point tracking (MPPT) algorithms tend to fall into local optimal solution, making it difficult to quickly and accurately find the global maximum power point of the PV system. Firstly, the male mayflies perform a global search to escape local optimal solution. Then, the female mayflies reduce system oscillations through random local search and mating. Finally, two different scenarios are set to verify the tracking accuracy and speed of the proposed algorithm. The results show that compared with particle swarm optimization (PSO) and grey wolf optimization (GWO), the EMA reduces PV output power oscillations during the search process. Compared to PSO and GWO, EMA has significantly improved the tracking speed, allowing the solar cells to maintain high output power even under partial shading.

    Keywords:PV generation; maximum power point tracking; improved Mayfly algorithm; partial shading; multi-peak; PV output characteristics

    猜你喜歡
    光伏發(fā)電
    光伏發(fā)電面臨的難題及建議分析
    科技傳播(2016年16期)2017-01-03 22:31:26
    電動(dòng)汽車換電站與光伏發(fā)電集成系統(tǒng)的容量配置研究
    我國(guó)光伏發(fā)電控制技術(shù)的發(fā)展現(xiàn)狀研究
    我國(guó)分布式光伏發(fā)電的現(xiàn)狀與展望
    淺談太陽(yáng)能光伏發(fā)電系統(tǒng)
    太陽(yáng)能光伏發(fā)電與廠房的一體化設(shè)計(jì)
    基于單片機(jī)的太陽(yáng)能路燈控制器
    光伏發(fā)電電氣系統(tǒng)設(shè)計(jì)優(yōu)化管理分析
    基于ARIMA模型的光伏客戶需求服務(wù)預(yù)測(cè)研究
    商情(2016年40期)2016-11-28 11:20:51
    大型并網(wǎng)光伏發(fā)電站選址分析
    欧美三级亚洲精品| 国产成人精品婷婷| 毛片一级片免费看久久久久| 嘟嘟电影网在线观看| 少妇精品久久久久久久| 高清午夜精品一区二区三区| 久久久成人免费电影| 国产黄片视频在线免费观看| 亚洲人与动物交配视频| 午夜免费男女啪啪视频观看| av在线观看视频网站免费| 九色成人免费人妻av| 99久久精品国产国产毛片| 久久久久久久大尺度免费视频| 亚洲欧美成人精品一区二区| 天天躁夜夜躁狠狠久久av| 在线观看三级黄色| 麻豆成人午夜福利视频| 国产淫语在线视频| 亚洲第一区二区三区不卡| 99久久精品热视频| 欧美日韩国产mv在线观看视频 | 国产精品福利在线免费观看| 国产精品一区二区在线不卡| 性色av一级| 免费看不卡的av| 国产高潮美女av| 久久久久视频综合| 久久久久网色| 精品一区二区三卡| 久久久久久九九精品二区国产| 久久久久久伊人网av| 久久午夜福利片| 如何舔出高潮| 国产成人精品福利久久| 少妇人妻久久综合中文| 天堂8中文在线网| 最近的中文字幕免费完整| 成人综合一区亚洲| 亚洲精品乱码久久久久久按摩| 亚洲成人手机| 国产精品久久久久久精品古装| 中国美白少妇内射xxxbb| 这个男人来自地球电影免费观看 | 久久99热这里只有精品18| 久久6这里有精品| 国产精品久久久久久久久免| 久久精品国产亚洲av涩爱| 国产精品偷伦视频观看了| 国产精品一区二区在线不卡| 免费看av在线观看网站| 性高湖久久久久久久久免费观看| 国产成人免费观看mmmm| 日本一二三区视频观看| 一个人免费看片子| 欧美精品人与动牲交sv欧美| 成人18禁高潮啪啪吃奶动态图 | 亚洲内射少妇av| 在线 av 中文字幕| 青春草视频在线免费观看| 亚洲av日韩在线播放| 在线观看三级黄色| 免费在线观看成人毛片| 中文资源天堂在线| 99热国产这里只有精品6| 啦啦啦视频在线资源免费观看| 美女视频免费永久观看网站| 大香蕉久久网| 亚洲精品乱久久久久久| 久久精品夜色国产| 国产又色又爽无遮挡免| 婷婷色综合大香蕉| 欧美日韩国产mv在线观看视频 | 久久ye,这里只有精品| 狂野欧美白嫩少妇大欣赏| 亚洲av中文字字幕乱码综合| 欧美日韩在线观看h| 欧美日韩视频高清一区二区三区二| 免费看av在线观看网站| 大陆偷拍与自拍| 亚洲综合色惰| 精品人妻熟女av久视频| 亚洲久久久国产精品| 亚洲精品日韩av片在线观看| 插阴视频在线观看视频| 午夜免费观看性视频| 天天躁夜夜躁狠狠久久av| 国产精品爽爽va在线观看网站| 秋霞在线观看毛片| 看十八女毛片水多多多| 自拍欧美九色日韩亚洲蝌蚪91 | 成人午夜精彩视频在线观看| 欧美日韩在线观看h| 久久青草综合色| 午夜精品国产一区二区电影| 精品国产露脸久久av麻豆| 丰满人妻一区二区三区视频av| 国产精品久久久久成人av| 99热这里只有是精品50| 国产亚洲5aaaaa淫片| 欧美日韩视频高清一区二区三区二| 如何舔出高潮| 欧美xxⅹ黑人| 18禁在线播放成人免费| 亚洲图色成人| 亚洲电影在线观看av| 亚洲自偷自拍三级| 日韩一区二区三区影片| 黄色配什么色好看| 日本vs欧美在线观看视频 | av在线老鸭窝| 午夜日本视频在线| 免费观看无遮挡的男女| 亚洲人成网站高清观看| 免费高清在线观看视频在线观看| 毛片一级片免费看久久久久| 国产一区有黄有色的免费视频| 在线看a的网站| 欧美少妇被猛烈插入视频| 天堂俺去俺来也www色官网| 一个人看视频在线观看www免费| 亚洲欧美精品专区久久| 王馨瑶露胸无遮挡在线观看| 人妻 亚洲 视频| 蜜臀久久99精品久久宅男| 黄色怎么调成土黄色| 在线亚洲精品国产二区图片欧美 | 人妻制服诱惑在线中文字幕| 夜夜爽夜夜爽视频| 一区二区三区精品91| 久久久久久久精品精品| 国产乱人偷精品视频| 永久网站在线| 视频中文字幕在线观看| 亚洲av在线观看美女高潮| 免费观看av网站的网址| 国产精品精品国产色婷婷| 美女视频免费永久观看网站| 亚洲欧美精品自产自拍| 男女国产视频网站| 国产深夜福利视频在线观看| 男人舔奶头视频| 综合色丁香网| 国产成人a区在线观看| 在线观看av片永久免费下载| 偷拍熟女少妇极品色| 亚洲av不卡在线观看| 老司机影院毛片| 日韩大片免费观看网站| 下体分泌物呈黄色| 国产精品一区二区三区四区免费观看| 又粗又硬又长又爽又黄的视频| 一本一本综合久久| 欧美激情国产日韩精品一区| 亚洲伊人久久精品综合| 久久精品夜色国产| 欧美日韩综合久久久久久| 啦啦啦在线观看免费高清www| 51国产日韩欧美| 亚洲久久久国产精品| 久久久久久伊人网av| 国产精品一区www在线观看| 国产毛片在线视频| 国产一区亚洲一区在线观看| 国产69精品久久久久777片| 十八禁网站网址无遮挡 | 最近手机中文字幕大全| 91狼人影院| 久久99热这里只有精品18| 国产精品久久久久久精品电影小说 | 国产黄色免费在线视频| 能在线免费看毛片的网站| 制服丝袜香蕉在线| 久久6这里有精品| 国产日韩欧美在线精品| 哪个播放器可以免费观看大片| 丰满少妇做爰视频| 97超碰精品成人国产| 日韩精品有码人妻一区| 97精品久久久久久久久久精品| 99国产精品免费福利视频| 3wmmmm亚洲av在线观看| 热99国产精品久久久久久7| 国产一区有黄有色的免费视频| 国产爱豆传媒在线观看| 久久97久久精品| 精品人妻一区二区三区麻豆| 午夜福利在线在线| 国产免费一区二区三区四区乱码| 国产午夜精品一二区理论片| 久久精品久久久久久噜噜老黄| 丰满少妇做爰视频| 欧美日本视频| 成人国产av品久久久| 91精品国产九色| 国产精品.久久久| 亚洲一级一片aⅴ在线观看| 熟女人妻精品中文字幕| 日本黄色日本黄色录像| 精品午夜福利在线看| 亚洲人与动物交配视频| 久久国产精品大桥未久av | 十八禁网站网址无遮挡 | 少妇的逼好多水| 久久精品国产亚洲网站| 99精国产麻豆久久婷婷| 亚洲国产欧美人成| 国产无遮挡羞羞视频在线观看| 国产黄片美女视频| 狂野欧美激情性xxxx在线观看| 亚洲av不卡在线观看| 国产精品久久久久久久电影| 欧美国产精品一级二级三级 | 黄色欧美视频在线观看| 日韩伦理黄色片| 亚洲国产色片| 欧美日韩综合久久久久久| 激情 狠狠 欧美| videos熟女内射| 蜜臀久久99精品久久宅男| 最后的刺客免费高清国语| 日韩在线高清观看一区二区三区| 大码成人一级视频| 欧美xxxx性猛交bbbb| 国产免费又黄又爽又色| 亚洲av成人精品一二三区| 18禁在线无遮挡免费观看视频| 丝瓜视频免费看黄片| 精品人妻偷拍中文字幕| 夜夜骑夜夜射夜夜干| 色5月婷婷丁香| 免费不卡的大黄色大毛片视频在线观看| 色视频www国产| 国产有黄有色有爽视频| 人人妻人人添人人爽欧美一区卜 | 亚洲怡红院男人天堂| 亚洲精品成人av观看孕妇| 夫妻性生交免费视频一级片| 色视频在线一区二区三区| 夜夜爽夜夜爽视频| 草草在线视频免费看| 少妇人妻精品综合一区二区| 制服丝袜香蕉在线| 97热精品久久久久久| 亚洲欧美成人精品一区二区| 日韩三级伦理在线观看| 久久久久久久久久久丰满| av在线观看视频网站免费| 亚洲四区av| 免费看光身美女| 国产亚洲欧美精品永久| 精品亚洲成国产av| 国产真实伦视频高清在线观看| 日本欧美国产在线视频| 国产男女超爽视频在线观看| 91久久精品国产一区二区成人| 成人午夜精彩视频在线观看| 午夜激情久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲5aaaaa淫片| 老司机影院成人| 成人免费观看视频高清| 日本爱情动作片www.在线观看| 久久韩国三级中文字幕| 国产伦理片在线播放av一区| 一级毛片电影观看| 一级毛片aaaaaa免费看小| 一区二区av电影网| 极品教师在线视频| 欧美亚洲 丝袜 人妻 在线| 中文字幕精品免费在线观看视频 | 中文乱码字字幕精品一区二区三区| 久久久久久久久久久免费av| 成人18禁高潮啪啪吃奶动态图 | 国产精品99久久99久久久不卡 | 成人免费观看视频高清| 久久精品久久精品一区二区三区| 色5月婷婷丁香| 精品久久久久久久久av| 涩涩av久久男人的天堂| 搡女人真爽免费视频火全软件| 麻豆乱淫一区二区| 亚洲精品自拍成人| 亚洲欧美一区二区三区国产| 两个人的视频大全免费| 王馨瑶露胸无遮挡在线观看| 黄色配什么色好看| 啦啦啦中文免费视频观看日本| 国产精品国产av在线观看| 中文精品一卡2卡3卡4更新| 亚洲精品乱码久久久v下载方式| 国产黄频视频在线观看| 午夜激情福利司机影院| 中文天堂在线官网| 亚洲国产精品成人久久小说| 国产免费一级a男人的天堂| 99久久综合免费| 精品少妇黑人巨大在线播放| 欧美日韩视频高清一区二区三区二| 久久久精品免费免费高清| 一级毛片 在线播放| 亚洲精品国产av蜜桃| 免费高清在线观看视频在线观看| 少妇熟女欧美另类| 中文字幕av成人在线电影| 久久国内精品自在自线图片| 免费观看av网站的网址| 99久久精品国产国产毛片| 国产精品一区二区在线观看99| 欧美日韩视频高清一区二区三区二| 少妇的逼水好多| 国产精品人妻久久久影院| 丝瓜视频免费看黄片| 日本午夜av视频| 三级国产精品片| 中文天堂在线官网| 久久久久精品性色| 成人毛片60女人毛片免费| 一级a做视频免费观看| 最近最新中文字幕大全电影3| 久久精品国产自在天天线| 国产无遮挡羞羞视频在线观看| 国产永久视频网站| 99久国产av精品国产电影| 国产精品一区二区在线观看99| av国产精品久久久久影院| 在线播放无遮挡| 97在线视频观看| 男人添女人高潮全过程视频| 国产成人精品婷婷| 夜夜爽夜夜爽视频| 99久久中文字幕三级久久日本| 亚洲国产成人一精品久久久| 在线观看一区二区三区| 免费看av在线观看网站| 日本色播在线视频| 特大巨黑吊av在线直播| 国产精品无大码| 久久99热6这里只有精品| 亚洲精品自拍成人| 麻豆精品久久久久久蜜桃| 亚洲精品乱码久久久久久按摩| 国产人妻一区二区三区在| 永久网站在线| 亚洲国产色片| 日韩av不卡免费在线播放| 亚洲欧美清纯卡通| av.在线天堂| 精品久久久精品久久久| 美女脱内裤让男人舔精品视频| 亚洲婷婷狠狠爱综合网| 国产熟女欧美一区二区| 国产在视频线精品| 91精品一卡2卡3卡4卡| 国产一区亚洲一区在线观看| 成人综合一区亚洲| 国产亚洲午夜精品一区二区久久| 国产伦精品一区二区三区四那| 一边亲一边摸免费视频| 在线观看免费视频网站a站| 国产欧美另类精品又又久久亚洲欧美| 国产一区有黄有色的免费视频| 中国三级夫妇交换| 日本欧美视频一区| 九色成人免费人妻av| 日韩电影二区| 黑人猛操日本美女一级片| 久久久精品免费免费高清| 国产精品欧美亚洲77777| 国产淫语在线视频| 国产亚洲最大av| 国产av一区二区精品久久 | 国产乱人偷精品视频| 久久久成人免费电影| 亚洲国产精品国产精品| 国产精品国产三级国产av玫瑰| 精品亚洲成a人片在线观看 | 亚洲,欧美,日韩| 国产综合精华液| 在线精品无人区一区二区三 | 人体艺术视频欧美日本| 美女cb高潮喷水在线观看| 99视频精品全部免费 在线| 午夜福利影视在线免费观看| 高清黄色对白视频在线免费看 | 久久精品熟女亚洲av麻豆精品| 2018国产大陆天天弄谢| 午夜免费观看性视频| 久久久精品94久久精品| 成人午夜精彩视频在线观看| 一级二级三级毛片免费看| 一级毛片aaaaaa免费看小| 夫妻午夜视频| 亚洲丝袜综合中文字幕| 亚洲国产精品专区欧美| 亚洲国产成人一精品久久久| 亚洲国产精品专区欧美| 亚洲不卡免费看| 国产极品天堂在线| 少妇 在线观看| 久久99精品国语久久久| 亚洲av不卡在线观看| 精品国产三级普通话版| 极品少妇高潮喷水抽搐| 日韩三级伦理在线观看| 国产av一区二区精品久久 | 精品一区二区三卡| 亚洲精品国产成人久久av| 国产伦理片在线播放av一区| 高清不卡的av网站| 人体艺术视频欧美日本| 成年人午夜在线观看视频| 男的添女的下面高潮视频| 久久综合国产亚洲精品| 80岁老熟妇乱子伦牲交| 欧美一区二区亚洲| h视频一区二区三区| 色吧在线观看| 黄色视频在线播放观看不卡| 男人和女人高潮做爰伦理| 久久久久久久久久成人| 少妇的逼好多水| 中文字幕人妻熟人妻熟丝袜美| 插逼视频在线观看| 七月丁香在线播放| 下体分泌物呈黄色| 成人毛片60女人毛片免费| 久久97久久精品| 国产精品99久久久久久久久| 久久影院123| av专区在线播放| 亚洲精品乱久久久久久| 视频区图区小说| 大香蕉97超碰在线| 一个人看视频在线观看www免费| 欧美另类一区| 日韩视频在线欧美| 18禁裸乳无遮挡动漫免费视频| 丝袜喷水一区| 97超视频在线观看视频| 亚洲人成网站在线播| av播播在线观看一区| 美女高潮的动态| 欧美日韩在线观看h| 久久人人爽人人片av| 国产亚洲5aaaaa淫片| 美女高潮的动态| 亚洲熟女精品中文字幕| 亚洲精品456在线播放app| 一区在线观看完整版| 国产一级毛片在线| 亚洲精品亚洲一区二区| 国内精品宾馆在线| 国产黄片美女视频| 久久国产乱子免费精品| 亚洲欧美清纯卡通| 你懂的网址亚洲精品在线观看| 日韩免费高清中文字幕av| 免费观看在线日韩| 婷婷色综合www| 五月开心婷婷网| 午夜免费鲁丝| 伦精品一区二区三区| 久久99热这里只频精品6学生| 免费看光身美女| 91精品国产国语对白视频| 伦精品一区二区三区| 97超视频在线观看视频| 久久久久久久亚洲中文字幕| 日韩不卡一区二区三区视频在线| 欧美 日韩 精品 国产| 男女边摸边吃奶| 简卡轻食公司| 成人亚洲精品一区在线观看 | 国产美女午夜福利| 亚洲精品亚洲一区二区| 亚洲精品成人av观看孕妇| 亚洲精品,欧美精品| 亚洲精品第二区| 在线观看av片永久免费下载| 99久久精品国产国产毛片| 久久久久久久久久久免费av| av在线播放精品| 国产毛片在线视频| 中国国产av一级| 18+在线观看网站| 麻豆国产97在线/欧美| 免费黄频网站在线观看国产| 又大又黄又爽视频免费| 免费看日本二区| 国产极品天堂在线| 少妇人妻精品综合一区二区| 天天躁夜夜躁狠狠久久av| 丝袜脚勾引网站| 国产免费一级a男人的天堂| 99久国产av精品国产电影| 肉色欧美久久久久久久蜜桃| 中文字幕精品免费在线观看视频 | 国产人妻一区二区三区在| 91午夜精品亚洲一区二区三区| 最近2019中文字幕mv第一页| 国产高清不卡午夜福利| 国产av精品麻豆| 国精品久久久久久国模美| 久久av网站| 视频中文字幕在线观看| 亚洲久久久国产精品| 在线观看美女被高潮喷水网站| 久久av网站| 观看美女的网站| 人妻制服诱惑在线中文字幕| 国精品久久久久久国模美| 国产精品.久久久| 18禁在线无遮挡免费观看视频| av天堂中文字幕网| 一本久久精品| 大话2 男鬼变身卡| 丝袜脚勾引网站| 亚洲av中文av极速乱| 99热6这里只有精品| av在线观看视频网站免费| 美女高潮的动态| av国产免费在线观看| 极品教师在线视频| 最新中文字幕久久久久| 亚洲不卡免费看| 国产精品一及| 国产成人a区在线观看| 少妇的逼好多水| 午夜视频国产福利| 成人18禁高潮啪啪吃奶动态图 | 亚洲国产欧美在线一区| 久久精品熟女亚洲av麻豆精品| 欧美国产精品一级二级三级 | 国产男女超爽视频在线观看| 亚洲va在线va天堂va国产| 三级国产精品欧美在线观看| 男人爽女人下面视频在线观看| 国产乱人视频| 亚洲精品成人av观看孕妇| 国产亚洲5aaaaa淫片| 国产精品伦人一区二区| 观看av在线不卡| 丝瓜视频免费看黄片| 亚洲国产av新网站| 日日啪夜夜爽| 五月玫瑰六月丁香| av专区在线播放| 狂野欧美激情性xxxx在线观看| 中国美白少妇内射xxxbb| 观看美女的网站| www.av在线官网国产| 大片免费播放器 马上看| 在线天堂最新版资源| 五月天丁香电影| 亚洲精品亚洲一区二区| 妹子高潮喷水视频| 久久精品久久久久久噜噜老黄| 高清在线视频一区二区三区| 黑人高潮一二区| h视频一区二区三区| 国产精品一区二区在线不卡| 国产一区二区三区综合在线观看 | 亚洲美女搞黄在线观看| 高清日韩中文字幕在线| 伦精品一区二区三区| 中文乱码字字幕精品一区二区三区| 少妇丰满av| 久久毛片免费看一区二区三区| 日韩欧美精品免费久久| 国产真实伦视频高清在线观看| 精品一区二区三卡| 毛片一级片免费看久久久久| 国产精品国产三级专区第一集| videos熟女内射| 精品少妇久久久久久888优播| 亚洲精品aⅴ在线观看| 亚洲精品乱码久久久v下载方式| 搡老乐熟女国产| 成年美女黄网站色视频大全免费 | 91久久精品国产一区二区成人| 最近中文字幕2019免费版| 中文字幕精品免费在线观看视频 | 免费人成在线观看视频色| 精品一区在线观看国产| 国产一区二区三区综合在线观看 | 久久婷婷青草| 人人妻人人添人人爽欧美一区卜 | 国产亚洲欧美精品永久| 中文资源天堂在线| 嫩草影院入口| 欧美精品人与动牲交sv欧美| 妹子高潮喷水视频| 精品一区二区免费观看| 久久这里有精品视频免费| 免费大片18禁| 99精国产麻豆久久婷婷| 日韩 亚洲 欧美在线| 亚洲欧美精品专区久久| 亚洲图色成人| 夫妻午夜视频| 国产精品一区二区三区四区免费观看| 国内精品宾馆在线| 国产爱豆传媒在线观看| 热re99久久精品国产66热6| 国产成人91sexporn| 精品一品国产午夜福利视频| 国产综合精华液| 小蜜桃在线观看免费完整版高清| 亚洲精品久久午夜乱码| 大陆偷拍与自拍| 国产人妻一区二区三区在| 亚洲精品日韩av片在线观看| 涩涩av久久男人的天堂| 激情五月婷婷亚洲|