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

    基于氣象數(shù)據(jù)外推法和顯著性分析的光伏自適應功率預測模型

    2025-03-14 00:00:00王麗婕張青山郝穎周穎邱敏孫沖
    太陽能學報 2025年2期
    關鍵詞:外推法光伏發(fā)電

    摘 要:分布式光伏電站裝機容量較小,一般不進行實時功率統(tǒng)計,難以直接建立功率預測模型。針對分布式光伏電站間光伏組件安裝型號和安裝方式多樣性的問題,基于氣象數(shù)據(jù)外推法和顯著性分析提出一種自適應功率預測模型。首先,利用顯著性分析篩選光伏組件標稱參數(shù)和氣象數(shù)據(jù)外推法特征參數(shù)之間的相關性,確定與特征參數(shù)顯著相關的標稱參數(shù)集合;然后,建立最小二乘支持向量機模型,擬合標稱參數(shù)集合和特征參數(shù)之間的自適應函數(shù)關系;最后,建立自適應功率預測模型,基于光伏組件安裝方式將環(huán)境溫度和輻照度轉(zhuǎn)換為板面溫度和輻照度,基于自適應函數(shù)選擇適合當前光伏組件型號的特征參數(shù),根據(jù)氣象數(shù)據(jù)外推法得到預測功率。利用光伏電站實際數(shù)據(jù)進行驗證,結(jié)果顯示,自適應功率預測模型能夠為不同型號的光伏組件選擇合適的特征參數(shù),通用性更強,相比于目前流行的簡化功率預測模型,多云天氣下的預測精度提升約2.34%。

    關鍵詞:光伏發(fā)電;功率預測;外推法;顯著性分析;自適應函數(shù)

    中圖分類號:TM615 " " " " 文獻標志碼:A

    0 引 言

    太陽能作為一種易于獲得的可再生清潔能源,具有巨大的資源潛力,是當下發(fā)展最迅猛,各國扶持力度最大的能源形式[1-2]。截至2022年,全球光伏裝機容量共計1073.1 GW,新增裝機容量230 GW,同比增長35.3%。據(jù)英國石油公司《2022年世界能源統(tǒng)計回顧》顯示,從2017年開始,每年新增光伏裝機容量近100 GW,光伏發(fā)電正逐漸成為電力供應的重要方式。

    光伏發(fā)電的功率存在顯著的波動和難以確定的特點,一旦大規(guī)模接入電網(wǎng),將給電力系統(tǒng)的安全穩(wěn)定運行帶來巨大挑戰(zhàn)。為此,建立高準確度的光伏功率超短期預測系統(tǒng),并根據(jù)預測結(jié)果科學制定調(diào)度計劃、預留足夠的消納空間,是提升光伏發(fā)電消納能力、確保電力系統(tǒng)安全運行的重要舉措。隨著光伏發(fā)電技術和功率預測技術的不斷發(fā)展,許多專家學者將研究重點轉(zhuǎn)向太陽輻照度預測,通過提高太陽輻照度的預測精度進一步提高光伏發(fā)電的預測精度[3-4]。分布式光伏電站由于裝機容量較小,一般不進行實時功率統(tǒng)計,難以直接建立功率預測模型。因此,首先建立太陽輻照度預測模型,再在預測太陽輻照度的基礎上進行功率預測是分布式光伏功率預測必要的環(huán)節(jié)。

    基于太陽輻照度的功率預測模型可分為統(tǒng)計模型和物理模型兩大類。統(tǒng)計模型是指采用神經(jīng)網(wǎng)絡等人工智能算法直接擬合太陽輻照度和光伏功率之間的非線性函數(shù)關系;物理模型是通過剖析光伏發(fā)電系統(tǒng)內(nèi)部各物理模塊,基于太陽電池的伏安特性曲線來建立預測模型。文獻[5-6]使用回歸神經(jīng)網(wǎng)絡建立太陽輻照度和光伏功率之間的預測模型,雖然取得不錯的精度,但預測效果并不穩(wěn)定,偶爾會出現(xiàn)偏差較大的預測點。物理模型基于光伏發(fā)電原理建立,對發(fā)電功率的追蹤效果較穩(wěn)定,因而被廣泛應用。太陽電池工程用數(shù)學模型提出光伏組件實際運行時的電壓和電流,可由參考狀態(tài)下光伏組件的電壓和電流修正得到,其中修正系數(shù)(也稱補償系數(shù))經(jīng)實驗數(shù)據(jù)擬合得到[7-9]。文獻[10-11]采用國家標準GB 50797—2012《光伏發(fā)電站設計規(guī)范》中的修正系數(shù)法對光伏電站的理論發(fā)電量進行了預測,修正系數(shù)法是通過對比光伏組件實際接收輻照度與標況輻照度(1000 W/m2)的比值修正裝機容量,同時依據(jù)經(jīng)驗參數(shù)選取合適的綜合效率系數(shù),對光伏電站的理論發(fā)電量進行預測,該模型參數(shù)簡單易求,是目前國內(nèi)采用范圍最廣的光伏理論出力預測模型之一。文獻[12-13]在使用修正系數(shù)法預測理論發(fā)電量時,除考慮實際輻照度和標況輻照度的比值外,還進一步考慮了由于輻照度和光伏組件溫度變化帶來的轉(zhuǎn)換效率的變化,從而進一步提升了功率預測精度。氣象數(shù)據(jù)外推法同時考慮了實際輻照度和實際溫度相對于標況輻照度和標況溫度(25 ℃)的變化,采用多個補償系數(shù)對光伏組件實際運行電壓、電流都進行修正,是當下一種較為完善的功率預測模型[14]。

    以上這些方法雖具有較高的預測精度,但模型一旦訓練完成,只能適用于當前目標光伏電站,一旦電力系統(tǒng)屬性發(fā)生改變,甚至還需重新擬合。尤其分布式光伏電站裝機容量小,相鄰電站的實際裝機模式和光伏組件型號都存在差異,應用此類模型存在工作量巨大的困難。為進一步提升效率和模型的通用性,簡化功率預測模型應運而生,該模型僅采用一個補償系數(shù)考慮實際發(fā)電量與實際輻照度和裝機容量之間的效率,通過大量的實際運行數(shù)據(jù)對該轉(zhuǎn)換效率進行擬合。簡化功率預測模型在缺失實驗數(shù)據(jù)的情況下可采用固定損失,即系統(tǒng)效率85%、逆變器效率95%、并網(wǎng)效率98%[15-16]。簡化功率預測模型參數(shù)簡單、通用性強,但無法具體量化多種因素對光伏組件運行的不同影響,因此損失了一部分精度。文獻[17]在簡化功率預測模型的基礎上,進一步考慮了溫度變化對發(fā)電功率的影響。文獻[18]在簡化功率預測模型的基礎上,具體分析了不同光伏組件材料和轉(zhuǎn)換效率的關系,進一步提出了優(yōu)化模型。

    針對分布式光伏電站間光伏組件安裝型號和安裝方式多樣性的問題,本文在現(xiàn)有功率預測模型的基礎上,對氣象數(shù)據(jù)外推法中需要擬合的特征參數(shù)進行研究,通過分析光伏組件自身特性和功率預測模型特征參數(shù)之間的函數(shù)關系,在保證功率預測精度的基礎上,提出通用性更高的自適應功率預測模型,以滿足分布式光伏電站功率預測系統(tǒng)在線建模需求。首先,對氣象數(shù)據(jù)外推法中的特征參數(shù)進行研究,通過實驗確定典型型號光伏組件的特征參數(shù);然后,對光伏組件的標稱參數(shù)和氣象數(shù)據(jù)外推法特征參數(shù)進行顯著性分析,選擇對特征參數(shù)有顯著影響的光伏組件標稱參數(shù)組成自變量集合;隨后,建立最小二乘支持向量機模型,確定自變量集合和特征參數(shù)之間的自適應函數(shù)關系,從而得到自適應功率預測模型;最后,將該模型應用到吉林省實際運行的分布式光伏場站中,驗證該模型的有效性。

    1 氣象數(shù)據(jù)外推法

    氣象數(shù)據(jù)外推法是一種基于綜合考慮環(huán)境條件、光伏電站位置、不同光伏組件特性和安裝方式等因素的方法,用于建立光伏發(fā)電單元的光電轉(zhuǎn)換模型[14],具體計算步驟如下:

    1)首先根據(jù)氣象監(jiān)測設備的實測水平輻照度和環(huán)境溫度,將水平輻照度轉(zhuǎn)換為光伏組件斜面的有效輻照度,將環(huán)境溫度轉(zhuǎn)換為光伏組件的有效溫度。

    根據(jù)太陽光的入射角度,光伏陣列的架設方式可劃分為固定式和跟蹤式兩種模式。顧名思義,固定式就是將光伏組件以一個固定的傾角進行安裝,不可移動。跟蹤式是通過安裝在可旋轉(zhuǎn)支架上的光伏組件,實現(xiàn)追蹤太陽高度角或方位角的功能。這樣一來,光伏組件能夠在太陽輻射時刻始終處于最佳角度,以獲得更高的太陽輻射吸收效率。

    以固定式安裝的光伏組件接收到的板面輻照度為:

    [G'=G·sin(α+β)·cosθ] (1)

    式中:[G']——板面接收輻照度,W/m2;[G]——太陽直接輻射,W/m2;[α]——太陽高度角,( °);[β]——光伏組件安裝傾角,( °);[θ]——太陽天頂角,( °)。

    跟蹤式分為單軸跟蹤和雙軸跟蹤,單軸跟蹤又分為跟蹤太陽高度角和跟蹤太陽天頂角。在考慮跟蹤式能完全跟蹤時,單軸跟蹤太陽高度角的板面接收輻照度的大小如式(2)所示,單軸跟蹤太陽天頂角的板面接收輻照度的大小如式(3)所示,雙軸跟蹤的板面接收輻照度的大小等同于實際直接輻射。

    [G'=G·cosθ] (2)

    [G'=G·sin(α+β)] (3)

    光伏組件的板面溫度由環(huán)境溫度和光伏組件的受熱特性確定:

    [T'=T+KG'] (4)

    式中:[T']——板面溫度,℃;[T]——環(huán)境溫度,℃;[K]——溫度修正系數(shù),根據(jù)光伏組件在標況下的板面溫度計算所得。

    2)根據(jù)光伏組件標準工況(輻照度為1000 W/m2,溫度為25 ℃)下的設備參數(shù),計算當前氣象條件下光伏組件輸出的直流功率。

    [Pm=UmIm] (5)

    [Im=IrGsGr(1+aΔT)] (6)

    [Um=Urln(e+bΔG)(1-cΔT)] (7)

    式中:[Pm]——最大理論出力功率,W;[Um]——最佳工作電壓,V;[Im]——最佳工作電流,A;[Ir]——標準工況下的最佳輸出電流,A;[Gs]——實際板面輻照度,W/m2;[Gr]——標準板面輻照度,取值為1000 W/m2;[a、b、c]——補償系數(shù),即需要根據(jù)實驗數(shù)據(jù)擬合得到的特征參數(shù);[ΔT]——板面實際溫度與標準溫度之差,℃;[Ur]——標準工況下的最佳輸出電壓,V;[ΔG]——實際輻照度與標準輻照度之差,W/m2。

    [a]由標況輻照度(1000 W/m2)下不同溫度的最佳工作電流擬合得到:

    [a=Ir-ImIrΔT] (8)

    [b]由標況溫度(25 ℃)下不同輻照度的最佳工作電壓擬合得到:

    [b=eUmUr-eΔG] (9)

    [c]由標況輻照度(1000 W/m2)下不同溫度的最佳工作電壓擬合得到:

    [c=Ur-UmUrΔT] (10)

    3)綜合考慮光伏組件的有效數(shù)量、光伏組件的老化、光伏組件的失配損失、光伏組件表面的塵埃遮擋、光伏組件至并網(wǎng)點的線路傳輸及電站用電損失、逆變器效率等因素,得到光伏電站并網(wǎng)點的交流功率。

    [Pth=N·Pm·K1·K2·K3·K4·η] (11)

    式中:[N]——并網(wǎng)運行的光伏組件的數(shù)量;[K1]——光伏組件老化損失系數(shù);[K2]——光伏組件失配損失系數(shù);[K3]——塵埃遮擋損失系數(shù);[K4]——線路傳輸及電站用電損失系數(shù);[η]——并網(wǎng)逆變器效率。一般情況下,綜合效率可取為95%。

    氣象數(shù)據(jù)外推法進行功率預測時,對于不同型號的光伏組件需重新擬合特征參數(shù)[a、b、c]以保證預測精度,因此通用性較差。分布式光伏電站裝機容量小,電站間光伏組件安裝型號及安裝方式差異較大,若使用氣象數(shù)據(jù)外推法進行功率預測,特征參數(shù)擬合工作量較大,預測成本較高。

    2 顯著性分析

    顯著性檢驗是回歸分析中用于判斷自變量集合[xi(i=1,2,…,m)]中每個自變量對因變量[y]的影響是否顯著的一種方式[19],本文通過顯著性分析光伏組件標稱參數(shù)和氣象數(shù)據(jù)外推法特征參數(shù)之間的相關關系,從而確定最優(yōu)的自變量集合。

    一般回歸分析的模型為:

    [y=β0+β1x1+…+βmxm+ζζ~N(0,σ2)] (12)

    式中:[β0, β1,…, βm]——回歸系數(shù)。

    因變量[y]與自變量[x1,x2,…, xm]之間是否存在如模型(12)所示的線性關系是需要檢驗的。顯然,如果所有的[β]都很小,則說明因變量[y]和自變量[x]之間的函數(shù)關系并不明顯。

    對于[n]個獨立觀測數(shù)據(jù)[(yi,xi1,…,xim),i=1,…,n,ngt;m],采用最小二乘法估計求解參數(shù)[β0, β1,…, βm]的值,即選擇估計值[βm=βm]時殘差平方和[Q]最小,此時y的估計值為:

    [y=β0+β0x1+…+βmxm] (13)

    [Q=i=1nζ2=i=1n(yi-β0-β1xi1-…-βmxim)2] (14)

    [β]服從正態(tài)分布:

    [β~N(β,σ2(XTX)-1)] (15)

    式中:[X]——自變量[x]構成的矩陣。

    對于殘差平方和[Q],其期望[E(Q)]滿足:

    [E(Q)=(n-m-1)σ2] (16)

    [Qσ2~χ2(n-m-1)] (17)

    若定義總平方和為[S=i=1n(yi-y)2],對[S]進行分解:

    [S=Q+U] (18)

    [U=i=1n(yi-y)2] (19)

    式中:[Q]——由式(14)定義的殘差平方和,反映隨機誤差對[y]的影響;[U]——回歸平方和,反映自變量對[y]的影響;[y]——基于[β]得到的擬合值。

    假設[H0:βj=0(j=1,…,m)],當該假設成立時,分解式(18)定義的Q、U有:

    [F=U/mQ/(n-m-1)~F(m,n-m-1)] (20)

    若[Flt;Fα(m,n-m-1)],則接受[H0];否則,拒絕[H0]。

    接受[H0]說明因變量[y]與自變量[x1,x2,…,xm]之間的線性關系不明顯,但可能存在非線性關系,如平方關系。

    因此定義用于衡量因變量[y]與自變量[x1,x2,…,xm]相關程度的指標,如用回歸平方和在總平方和中的比值定義復判定系數(shù),如式(21)所示。

    [R2=URSST] (21)

    [R=R2]稱為復相關系數(shù),[R]越大,表示[y]與[x1,x2,…,xm]的相關關系越密切,通常[R]大于0.8才認為相關關系成立。

    3 最小二乘支持向量機模型

    最小二乘支持向量機(least squares support vector machine,LS-SVM)是支持向量機(support vector machine,SVM)的一種擴展。最小二乘支持向量機在非線性函數(shù)擬合及數(shù)據(jù)處理方面具有顯著優(yōu)勢。首先,LS-SVM能夠?qū)⒃继卣饔成涞礁呔S空間,因此面對復雜的非線性問題具有高效的擬合能力,且LS-SVM避免了復雜的二次規(guī)劃問題,能更快速找到全局最優(yōu)解。其次,LS-SVM在有限的樣本數(shù)量下也能保持高精度的預測和泛化能力,且能在存在噪聲或異常值的情況下保持穩(wěn)定的預測性能。

    支持向量機采用非線性映射將自變量映射到一個高維特征空間,建立自變量與因變量的回歸模型為:

    [y=f(x)=ωTx+d] (22)

    式中:[ω]——特征空間中的權向量;d——偏置量。

    然后定義一個損失函數(shù):

    [Lf(xi)-yi=0,f(xi)-yilt;εf(xi)-yi-ε," 其他] (23)

    式中:[ε]——不敏感函數(shù)。

    根據(jù)結(jié)構風險最小化原則,得到線性回歸的凸二次規(guī)劃問題。

    [min12ω2+Ci=1Nl(f(xi)+yi)s.t. -ε≤yi-f(xi)≤ε, i=1,…,N] (24)

    式中:[C]——懲罰因子。

    求解該問題時引入拉格朗日函數(shù),可得: [L(ω,b,ξi,ξ*i)=12ω2+Ci=1N(ξi+ξ*i)-i=1Nαi(ε+ξi-yi+ωxi+b)- i=1Nα*i(ε+ξ*i-yi+ωxi-b)-i=1N(ηiξi+η*iξ*i)] (25)

    式中:[ξi]、[ξ*i]——引入的松弛變量;[αi]、[α*i]、[ηi]、[η*i]——拉格朗日乘子。

    根據(jù)KKT(Ksrush-Kuhn-Tucker)條件求解公式可得回歸擬合函數(shù)[f(x)=i=1N(-αi-αi*)K(x,xi)+b],[K(x,xi)]為核函數(shù)。

    SVM求解需借助二次規(guī)劃,當樣本數(shù)據(jù)過大時必將增加規(guī)劃問題的復雜程度,耗費大量的計算機內(nèi)存和時間。LS-SVM最大的改變就是引入誤差平方項,不僅極大簡化了運算過程,對異常數(shù)據(jù)擁有更好的包容性,并進一步提升了運算精度。LS-SVM簡化運算后即變?yōu)樽钚《饲蠼獾男问?,LS-SVM的目標函數(shù)為:

    [minJ(ω,e)=12ω2+r2i=1Ne2is.t. yi=ωT?(xi)+b+ei," i=1,2,…,N] (26)

    式中:[r]——懲罰因子;[e2]——誤差平方和。

    引入拉格朗日函數(shù)求解得:

    [L(ω,b,e,a)=J(ω,e)-i=1mαiω?(xi)+b+ei-yi] (27)

    最后根據(jù)KKT條件求解得回歸擬合函數(shù)為:

    [f(x)=i=1mαiK(x,xi)+b] (28)

    4 自適應功率預測模型

    自適應功率預測模型基于氣象數(shù)據(jù)外推法建立,通過尋找光伏組件標稱參數(shù)和氣象數(shù)據(jù)外推法特征參數(shù)之間的自適應函數(shù)關系,進一步構建自適應功率預測模型。自適應功率預測模型的流程如圖1所示。

    首先根據(jù)光伏組件在不同條件下的工作狀態(tài),求得其氣象數(shù)據(jù)外推法特征參數(shù);然后對光伏組件標稱參數(shù)和氣象數(shù)據(jù)外推法特征參數(shù)進行顯著性分析,確定最優(yōu)的標稱參數(shù)自變量集合;之后建立最小二乘支持向量機模型,擬合標稱參數(shù)自變量集合和特征參數(shù)之間的自適應函數(shù)關系;最后將自適應函數(shù)關系式代入氣象數(shù)據(jù)外推法,完成自適應功率預測模型的建立,并對功率進行預測。

    5 應用研究

    本文實驗分為兩個步驟進行:1)自適應功率預測模型的訓練;2)基于輻照度的功率預測。為保證自適應功率預測模型的通用性,本文使用品牌A和B共11個型號的光伏組件進行訓練,輸入選擇光伏組件的標稱參數(shù)及其非線性升維組合,共9個特征,輸出為在實驗條件下計算得到的光伏組件特征參數(shù),共3個。該11個型號光伏組件標稱參數(shù)及其氣象數(shù)據(jù)外推法特征參數(shù)如表1所示,其中短路電流和開路電壓主要受光伏組件材料影響,是反映光伏組件轉(zhuǎn)換效率的重要參數(shù)。

    在基于輻照度進行功率預測時,輸入特征為輻照度(實際輻照度或預測輻照度)、光伏組件標稱參數(shù)(峰值功率、開路電壓、短路電流)、安裝參數(shù)(安裝傾角、安裝數(shù)量、安裝方式),共7維,其中光伏組件的安裝參數(shù)用于確定板面能夠接收的輻照度大小。本文選用2022年12月吉林省兩個光伏電站的實際運行數(shù)據(jù)進行驗證,每個光伏電站安裝的光伏組件型號和方式是統(tǒng)一的,且都不在訓練選用的光伏組件之列。另外,使用預測輻照度進行功率預測時需先訓練太陽輻照度預測模型。本文選用的太陽輻照度預測模型來自之前的研究成果[20],訓練數(shù)據(jù)集為2020和2021兩年的光伏電站實際運行數(shù)據(jù),輸入特征使用晴空輻照度、數(shù)值天氣預報(numerical weather prediction,NWP)特征(溫度、濕度、降水量、短波輻射、熱感通量)和衛(wèi)星云圖紋理特征(像素均值、熵、相關性、對比度),共10維,時間尺度為4 h,時間分辨率為15 min。

    結(jié)果采用歸一化均方根誤差(normalized root mean square error,NRMSE)和歸一化平均絕對誤差(normalized mean absolute error, NMAE)進行評判。

    [eNRMSE=i=1n(Yi-Yi)2n×100%] (29)

    [eNMAE=i=1nYi-Yin×100%] (30)

    式中:[n]——樣本點的數(shù)量;[Yi]——實際值;[Yi]——預測值。這里[Yi]和[Yi]都是除以實際功率的最大值后歸一化的值。

    5.1 顯著性分析結(jié)果

    從氣象數(shù)據(jù)外推法模型求取特征參數(shù)的過程中可看出特征參數(shù)與電壓電流并非呈線性關系。結(jié)合求取特征參數(shù)的公式,對標稱參數(shù)進行非線性處理,組成新的自變量集合[{Voc, Isc, Pm, V2oc, I2sc, P2m, ln(Voc), ln(Isc), ln(Pm)}]。然后將該集合與特征參數(shù)[a、b、c]進行顯著性分析,得其復相關系數(shù)分別為0.98、0.97、0.94,說明該自變量集合與[a、b、c]相關關系成立,該自變量集合與[a、b、c]的殘差分布如圖2所示,可看出殘差的置信區(qū)間均包含零點,且不存在異常值。

    5.2 特征參數(shù)擬合結(jié)果

    在5.1節(jié)顯著性分析的基礎上,得到光伏組件標稱參數(shù)自變量集合[{Voc, Isc, Pm, V2oc, I2sc, P2m, ln(Voc), ln(Isc), ln(Pm)}]和氣象數(shù)據(jù)外推法特征參數(shù)([a、b、c])具有顯著相關性。之后建立LS-SVM模型,對標稱參數(shù)自變量集合[{Voc, Isc, Pm, V2oc, I2sc, P2m, ln(Voc), ln(Isc), ln(Pm)}]和特征參數(shù)([a、b、c])進行擬合,求其特征參數(shù)自適應函數(shù)關系式,即:

    [a=La(Voc,Isc,Pm,V2oc,I2sc,P2m,ln(Voc),ln(Isc),ln(Pm))] (31)

    [b=Lb(Voc,Isc,Pm,V2oc,I2sc,P2m,ln(Voc),ln(Isc),ln(Pm))] (32)

    [c=Lc(Voc,Isc,Pm,V2oc,I2sc,P2m,ln(Voc),ln(Isc),ln(Pm))] (33)

    式中:[L(·)]——LS-SVM擬合函數(shù)。

    采用表1中的數(shù)據(jù)運用LS-SVM模型進行擬合,擬合結(jié)果如圖3所示(按參數(shù)升序進行排列),可看出基于LS-SVM的擬合曲線和實際曲線具有良好的對應關系,進一步論證了顯著性分析篩選的標稱參數(shù)自變量集合和特征參數(shù)之間具有良好的相關性,且LS-SVM模型在函數(shù)擬合方面也具有較好的性能。

    5.3 功率預測結(jié)果

    上文僅對比了訓練樣本的函數(shù)擬合效果,為進一步驗證自適應功率預測模型的通用性,本節(jié)以訓練樣本之外的吉林省兩個分布式光伏電站S1和S2的實際運行數(shù)據(jù)對本文所提自適應功率預測模型進行驗證,并與目前領域內(nèi)常使用的另一種通用模型即簡化功率預測模型進行對比。

    簡化功率預測模型為:

    [P=ηSI[1-0.005(t+25)]] (34)

    式中:[P]——理論發(fā)電量,W;[η]——轉(zhuǎn)換效率;[S]——光伏系統(tǒng)有效面積,m2;[I]——輻射量,W/m2;[t]——環(huán)境溫度,℃。

    根據(jù)上文對氣象數(shù)據(jù)外推法的介紹,氣象數(shù)據(jù)外推法中的變量包括太陽輻照度、環(huán)境溫度、溫度修正系數(shù)、特征參數(shù)[a、b、c],因此可將氣象數(shù)據(jù)外推法公式簡化為:

    [P=0.95N·f(a,b,c,G,T,K)] (35)

    式中:0.95——綜合損失;[N]——太陽電池數(shù)量;[f(·)]——氣象數(shù)據(jù)外推法計算方法。

    將式(31)~式(33)代入式(35),通過光伏組件標稱的開路電壓、短路電流、峰值功率為不同型號的光伏組件自適應的選擇特征參數(shù)[a、b、c],進一步簡化模型步驟,得到自適應功率預測模型為:

    [P=0.95N·f(La,Lb,Lc,G,T,K)] (36)

    為排除太陽輻照度預測精度對功率預測結(jié)果造成干擾,首先使用實際輻照度輸入自適應功率預測模型進行驗證,將S1和S2電站的實際輻照度和實際溫度代入式(34)和式(36),結(jié)果如表2和圖4所示。

    從圖4可看出,兩種功率預測模型計算得到的功率趨勢大致相同,而自適應功率預測模型的結(jié)果更接近實際值。由于簡化功率預測模型采用相同的功率損失,所以該預測模型對不同型號光伏組件的適用性不同,而自適應功率預測模型可根據(jù)光伏組件的型號自適應的調(diào)節(jié)特征參數(shù),從而能更好適應不同的光伏組件。

    下面利用預測的輻照度來預測光伏功率,首先使用文獻[20]中的太陽輻照度預測模型對S1電站的太陽輻照度進行預測,然后將S1電站的預測輻照度分別代入簡化功率預測模型和自適應功率預測模型進行功率預測,結(jié)果如圖5和表3所示。

    通過對比晴天和陰天天氣下的預測效果,發(fā)現(xiàn)無論在何種天氣下,自適應功率預測模型的預測效果都具有更明顯優(yōu)勢。自適應功率預測模型在晴天天氣下的均方根誤差為2.90%,多云天氣下的均方根誤差為5.87%。與簡化功率預測模型相比,自適應功率預測模型在正午時預測效果更好,尤其是晴天天氣下,自適應功率預測模型的預測功率與實際功率曲線明顯更接近。再一次說明了自適應功率預測模型能夠通過自適應的選擇特征參數(shù),匹配不同光伏組件的不同特點。

    6 結(jié) 論

    針對基于輻照度的光伏功率預測環(huán)節(jié),本文提出一種自適應功率預測模型。自適應功率預測模型通過顯著性分析光伏組件標稱參數(shù)和氣象數(shù)據(jù)外推法特征參數(shù)之間的相關性,利用支持向量機模型為不同型號的光伏組件選取合適的特征參數(shù)。通過實例驗證,得出以下主要結(jié)論:

    1)顯著性分析發(fā)現(xiàn)光伏組件標稱參數(shù)的組合與氣象數(shù)據(jù)外推法特征參數(shù)之間具有相關性。

    2)LS-SVM模型可很好地擬合上述相關性,為不同光伏組件自適應選取合適的特征參數(shù)。

    3)與目前常用的簡化功率預測模型相比,自適應功率預測模型通用性和適用性更強,預測精度更高。

    4)該自適應功率預測模型不僅適用于已有光伏電站根據(jù)預測輻照度進行功率預測,也適用于光伏電站的選址中根據(jù)實際輻照度進行功率計算。

    [參考文獻]

    [1] 賴昌偉, 黎靜華, 陳博, 等. 光伏發(fā)電出力預測技術研究綜述[J]. 電工技術學報, 2019, 34(6): 1201-1217.

    LAI C W, LI J H, CHEN B, et al. Review of photovoltaic power output prediction technology[J]. Transactions of China Electrotechnical Society, 2019, 34(6): 1201-1217.

    [2] 陳元峰, 馬溪原, 程凱, 等. 基于氣象特征量選取與SVM模型參數(shù)優(yōu)化的新能源超短期功率預測[J]. 太陽能學報, 2023, 44(12): 568-576.

    CHEN Y F, MA X Y, CHENG K, et al. Ultra-short-term power forecast of new energy based on meteorological feature selection and SVM model parameter optimization[J]. Acta energiae solaris sinica, 2023, 44(12): 568-576.

    [3] 倪超, 王聰, 朱婷婷, 等. 基于CNN-Bi-LSTM的太陽輻照度超短期預測[J]. 太陽能學報, 2022, 43(3): 197-202.

    NI C, WANG C, ZHU T T, et al. Super-short-term forecast of solar irradiance based on CNN-Bi-LSTM[J]. Acta energiae solaris sinica, 2022, 43(3): 197-202.

    [4] 吳小濤, 袁曉輝, 袁艷斌, 等. 基于二階數(shù)據(jù)分解算法和蝗蟲優(yōu)化混合核LSSVM的太陽輻照度預測模型研究[J]. 可再生能源, 2021, 39(7): 899-907.

    WU X T, YUAN X H, YUAN Y B, et al. Solar irradiance prediction model based on two-stage data decomposition algorithm and hybrid kernel LSSVM optimized by grasshopper optimization algorithm[J]. Renewable energy resources, 2021, 39(7): 899-907.

    [5] RILEY D M, VENAYAGAMOORTHY G K. Characterization and modeling of a grid-connected photovoltaic system using a recurrent neural network[C]//The 2011 International Joint Conference on Neural Networks. San Jose, CA, USA, 2011: 1761-1766.

    [6] NGUYEN D D, LEHMAN B, KAMARTHI S. Performance evaluation of solar photovoltaic arrays including shadow effects" "using" "neural" "network[C]//2009 IEEE Energy Conversion Congress and Exposition. San Jose, CA, 2009: 3357-3362.

    [7] 蘇建徽, 余世杰, 趙為, 等. 硅太陽電池工程用數(shù)學模型[J]. 太陽能學報, 2001, 22(4): 409-412.

    SU J H, YU S J, ZHAO W, et al. Investigation on engineering analytical model of silicon solar cells[J]. Acta energiae solaris sinica, 2001, 22(4): 409-412.

    [8] 傅望, 周林, 郭珂, 等. 光伏電池工程用數(shù)學模型研究[J]. 電工技術學報, 2011, 26(10): 211-216.

    FU W, ZHOU L, GUO K, et al. Research on engineering analytical model of solar cells[J]. Transactions of China Electrotechnical Society, 2011, 26(10): 211-216.

    [9] 彭湃, 程漢湘, 陳杏燦, 等. 光伏電池工程用數(shù)學模型及其模型的應用研究[J]. 電源技術, 2017, 41(5): 780-782, 789.

    PENG P, CHENG H X, CHEN X C, et al. Engineering analytical model of photovoltaic cell and its application research[J]. Chinese journal of power sources, 2017, 41(5): 780-782, 789.

    [10] 楊娜. 國電內(nèi)蒙古察右前旗光伏電站的出力預測研究[D]. 北京: 華北電力大學, 2016.

    YANG N. Prediction of the output of Guodian Inner Mongolia" "Chayouqianqi" "photovoltaic" "power" "plant[D]. Beijing: North China Electric Power University, 2016.

    [11] 李恩茂. 10 MW光伏發(fā)電項目設計[D]. 大連: 大連海事大學, 2019.

    LI E M. Design of 10 MW photovoltaic power generation project[D]. Dalian: Dalian Maritime University, 2019.

    [12] DREWS A, DE KEIZER A C, BEYER H G, et al. Monitoring and remote failure detection of grid-connected PV" systems" based" on" satellite" observations[J]." Solar energy, 2007, 81(4): 548-564.

    [13] LORENZ E, SCHEIDSTEGER T, HURKA J, et al. Regional PV power prediction for improved grid integration[J]. Progress in photovoltaics: research and applications, 2011, 19(7): 757-771.

    [14] 國家能源局. 光伏電站消納監(jiān)測統(tǒng)計管理辦法[EB/OL].(2021-12-3). http://zfxxgk.nea.gov.cn/2021-12/03/c_1310383862. htm.

    National Energy Administration. Management measures for consumption monitoring and statistics of photovoltaic power stations[EB/OL]." (2021-12-3)." http://zfxxgk.nea.gov.cn/2021-12/03/c_1310383862.htm.

    [15] 陳祥. 基于機理模型的并網(wǎng)光伏電站實時效率分析[J]. 太陽能, 2011(17): 12-14, 49.

    CHEN X. Real-time efficiency analysis of grid-connected photovoltaic power station based on mechanism model[J]. Solar energy, 2011(17): 12-14, 49.

    [16] 張峻嶺, 殷建英, 王文軍. 鄂爾多斯新能源產(chǎn)業(yè)示范區(qū)的風光互補最優(yōu)容量匹配[J]. 電力與能源, 2011, 32(3): 224-227.

    ZHANG J L, YIN J Y, WANG W J. Optimal capacity for the complementation of wind-solar power located in Ordos new" energy" demonstration" zone[J]." Power" amp;" energy, 2011, 32(3): 224-227.

    [17] 張雪莉, 劉其輝, 馬會萌, 等. 光伏電站輸出功率影響因素分析[J]. 電網(wǎng)與清潔能源, 2012, 28(5): 75-81.

    ZHANG X L, LIU Q H, MA H M, et al. Analysis of influencing factors of output power of photovoltaic power plant[J]. Power system and clean energy, 2012, 28(5): 75-81.

    [18] 王利珍, 譚洪衛(wèi), 莊智, 等. 基于GIS平臺的我國太陽能光伏發(fā)電潛力研究[J]. 上海理工大學學報, 2014, 36(5): 491-496.

    WANG L Z, TAN H W, ZHUANG Z, et al. Evaluation of the photovoltaic solar energy potential in China based on GIS platform[J]. Journal of University of Shanghai for Science and Technology, 2014, 36(5): 491-496.

    [19] 田佳垚, 馮自平, 胡亞飛, 等. 基于多元回歸模型的燃氣熱泵系統(tǒng)制熱性能分析[J]. 熱能動力工程, 2023, 38(6): 129-136.

    TIAN J Y, FENG Z P, HU Y F, et al. Analysis of heating performance of gas engine-driven heat pump system based on multiple regression model[J]. Journal of engineering for thermal energy and power, 2023, 38(6): 129-136.

    [20] 張青山, 王麗婕, 郝穎, 等. 基于衛(wèi)星云圖和晴空模型的分布式光伏電站太陽輻照度超短期預測[J]. 高電壓技術, 2022, 48(8): 3271-3281.

    ZHANG Q S, WANG L J, HAO Y, et al. Ultra-short-term solar irradiance prediction of distributed photovoltaic power stations based on satellite cloud images and clear sky model[J]. High voltage engineering, 2022, 48(8): 3271-3281.

    PHOTOVOLTAIC ADAPTIVE POWER PREDICTION MODEL BASED ON METEOROLOGICAL DATA EXTRAPOLATION AND

    SIGNIFICANCE ANALYSIS

    Wang Lijie1,Zhang Qingshan1,Hao Ying2,Zhou Ying3,Qiu Min3,Sun Chong4

    (1. School of Automation, Beijing Information Science amp; Technology University, Beijing 100192, China;

    2. Tangshan Research Institute, Beijing Institute of Technology, Tangshan 063000, China;

    3. Beijing Key Laboratory of Demand Side Multi?Energy Carriers Optimization and Interaction Technique(China Electric Power Research Institute), Beijing 100192, China;

    4. State Grid Hebei Marketing Service Center, Shijiazhuang 050081, China)

    Abstract:Distributed photovoltaic power stations typically possess a relatively small installed capacity and generally do not conduct real-time power statistics, posing challenges to directly establish power prediction models. This paper proposes an adaptive power prediction model based on meteorological data extrapolation and significance analysis aiming at the diversity of installation models and methods for photovoltaic panels in distributed photovoltaic power stations. Firstly, significance analysis is conducted on the nominal parameters of the photovoltaic panel and the characteristic parameters in the meteorological data extrapolation method to determine the set of nominal parameters. Then, establish a least squares support vector machine model to fit the adaptive function relationship between the nominal parameter set and the feature parameters. Finally, an adaptive power prediction model is established, which converts the ambient temperature and irradiance into the temperature and irradiance of the panel based on the installation mode of photovoltaic panels, selects the characteristic parameters suitable for the current photovoltaic cell model based on the adaptive function, and inputs the meteorological data extrapolation method to obtain the predicted power. The research results indicate that the adaptive power prediction model can adaptively select appropriate feature parameters for different types of photovoltaic panels, with stronger universality. Compared with the current popular simplified power prediction model, the prediction accuracy is improved by approximately 2.34%.

    Keywords:photovoltaic power generation; power forecasting; extrapolation; significance analysis; adaptive function

    猜你喜歡
    外推法光伏發(fā)電
    一種基于機械臂的外推天線增益標定系統(tǒng)
    中國新通信(2023年3期)2023-06-24 03:00:06
    I 型裂紋應力強度因子的有限元計算方法
    濱海相河口大型沉井基礎設計荷載效應估計
    人民黃河(2018年8期)2018-09-10 16:05:55
    基于單片機的太陽能路燈控制器
    光伏發(fā)電電氣系統(tǒng)設計優(yōu)化管理分析
    基于ARIMA模型的光伏客戶需求服務預測研究
    商情(2016年40期)2016-11-28 11:20:51
    大型并網(wǎng)光伏發(fā)電站選址分析
    中國市場(2016年41期)2016-11-28 05:37:35
    中國西北地區(qū)光伏發(fā)電的環(huán)境影響評估
    時代金融(2016年27期)2016-11-25 16:14:21
    應用于教學的太陽能電源裝置的設計
    用熱線法測定氣體熱導率的儀器改進
    国产欧美日韩精品亚洲av| 伊人久久精品亚洲午夜| 亚洲国产欧美网| 成人一区二区视频在线观看| 中文资源天堂在线| 人人妻,人人澡人人爽秒播| 97超视频在线观看视频| 亚洲成av人片免费观看| 亚洲在线观看片| 日本成人三级电影网站| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看电影| 亚洲人成电影免费在线| aaaaa片日本免费| 麻豆国产av国片精品| 一进一出抽搐gif免费好疼| 午夜福利高清视频| 91麻豆精品激情在线观看国产| 国产一区二区在线av高清观看| 19禁男女啪啪无遮挡网站| 夜夜爽天天搞| 夜夜夜夜夜久久久久| xxx96com| 欧美乱色亚洲激情| 日本成人三级电影网站| 亚洲 欧美 日韩 在线 免费| 少妇熟女aⅴ在线视频| 岛国在线观看网站| 又紧又爽又黄一区二区| 国产高潮美女av| 亚洲自拍偷在线| 嫩草影院精品99| 最近在线观看免费完整版| 久久人妻av系列| 中文亚洲av片在线观看爽| 久久婷婷人人爽人人干人人爱| 成人午夜高清在线视频| 成人18禁在线播放| 欧美三级亚洲精品| 国产精品一区二区三区四区久久| 欧美乱妇无乱码| 亚洲av成人精品一区久久| 九九在线视频观看精品| 久久中文看片网| 日韩高清综合在线| 亚洲人成网站高清观看| 国产精品综合久久久久久久免费| 欧美日韩中文字幕国产精品一区二区三区| 欧美色欧美亚洲另类二区| 久久久久久久久中文| 中文字幕av成人在线电影| 床上黄色一级片| 麻豆成人午夜福利视频| 国产精品野战在线观看| 2021天堂中文幕一二区在线观| 90打野战视频偷拍视频| 欧美中文日本在线观看视频| 美女被艹到高潮喷水动态| 美女被艹到高潮喷水动态| 观看美女的网站| 久久这里只有精品中国| 久久99热这里只有精品18| 久久精品影院6| 一个人免费在线观看电影| 美女cb高潮喷水在线观看| 久9热在线精品视频| 国产97色在线日韩免费| 欧美乱色亚洲激情| 男人舔女人下体高潮全视频| 国产精品精品国产色婷婷| 桃色一区二区三区在线观看| 特级一级黄色大片| 婷婷六月久久综合丁香| 欧美色视频一区免费| 国产成人av激情在线播放| 成年女人看的毛片在线观看| 长腿黑丝高跟| av片东京热男人的天堂| 国产黄片美女视频| 精品乱码久久久久久99久播| 国产亚洲av嫩草精品影院| svipshipincom国产片| 岛国在线观看网站| 男人舔女人下体高潮全视频| 久久久久国内视频| 国产精品永久免费网站| 岛国在线观看网站| 久久中文看片网| 色精品久久人妻99蜜桃| 天堂影院成人在线观看| 欧美日韩黄片免| 色吧在线观看| 免费观看的影片在线观看| 中出人妻视频一区二区| 国产高清videossex| 久久久久性生活片| 亚洲精品成人久久久久久| 亚洲成人免费电影在线观看| 亚洲成av人片在线播放无| 日韩国内少妇激情av| 一级黄色大片毛片| 亚洲专区中文字幕在线| 91字幕亚洲| 午夜免费成人在线视频| 最后的刺客免费高清国语| 成人一区二区视频在线观看| 欧美另类亚洲清纯唯美| 亚洲人成伊人成综合网2020| 久久亚洲真实| 99在线视频只有这里精品首页| 精品日产1卡2卡| 日韩大尺度精品在线看网址| 三级国产精品欧美在线观看| 日韩欧美三级三区| 色播亚洲综合网| 久久久国产成人精品二区| 免费一级毛片在线播放高清视频| 色哟哟哟哟哟哟| 欧美bdsm另类| 狠狠狠狠99中文字幕| 黄色片一级片一级黄色片| 国产成人啪精品午夜网站| 欧美日韩精品网址| 岛国视频午夜一区免费看| 国产精品久久视频播放| 国内精品久久久久精免费| 欧美最新免费一区二区三区 | 色老头精品视频在线观看| eeuss影院久久| 亚洲av第一区精品v没综合| 中出人妻视频一区二区| 国产老妇女一区| 久久久久久久久久黄片| 少妇人妻精品综合一区二区 | 国产精品99久久99久久久不卡| 草草在线视频免费看| 中国美女看黄片| 真实男女啪啪啪动态图| 97碰自拍视频| 久久久久久久精品吃奶| 国产在线精品亚洲第一网站| 国产精品久久久久久亚洲av鲁大| 99国产精品一区二区蜜桃av| 欧美大码av| 三级国产精品欧美在线观看| 欧美日韩乱码在线| 天堂影院成人在线观看| 毛片女人毛片| 国内精品久久久久久久电影| 成人国产综合亚洲| 国产精品久久久久久精品电影| 国产精品女同一区二区软件 | 两个人看的免费小视频| 久久草成人影院| 少妇熟女aⅴ在线视频| 国内少妇人妻偷人精品xxx网站| 久久久久久九九精品二区国产| 免费在线观看影片大全网站| 欧美成人a在线观看| 亚洲国产欧美人成| 哪里可以看免费的av片| 国产伦人伦偷精品视频| 18禁裸乳无遮挡免费网站照片| 51国产日韩欧美| 久久草成人影院| 久久国产精品人妻蜜桃| 中文在线观看免费www的网站| 亚洲一区高清亚洲精品| 真人一进一出gif抽搐免费| 国产亚洲欧美在线一区二区| 有码 亚洲区| 亚洲精品成人久久久久久| 日本与韩国留学比较| 99久久久亚洲精品蜜臀av| 亚洲精品亚洲一区二区| www日本黄色视频网| 国产精品一区二区三区四区久久| 九九久久精品国产亚洲av麻豆| 少妇的逼好多水| 中文字幕熟女人妻在线| 香蕉丝袜av| 99国产精品一区二区蜜桃av| 男人和女人高潮做爰伦理| 国产精品久久久久久久久免 | 免费高清视频大片| 日本三级黄在线观看| 午夜影院日韩av| 欧美高清成人免费视频www| 精品国产亚洲在线| 精品国产三级普通话版| 精品久久久久久,| 亚洲一区二区三区色噜噜| 欧美区成人在线视频| 亚洲久久久久久中文字幕| 国产精品一区二区免费欧美| 国产精品99久久久久久久久| 亚洲av五月六月丁香网| 波多野结衣高清无吗| 精品国产美女av久久久久小说| 国产精品自产拍在线观看55亚洲| 18禁美女被吸乳视频| av视频在线观看入口| 国产精品久久久久久人妻精品电影| 国产成人影院久久av| 18禁黄网站禁片午夜丰满| 亚洲成a人片在线一区二区| 国产精品亚洲av一区麻豆| 中文字幕人妻熟人妻熟丝袜美 | 最新中文字幕久久久久| 国产精品三级大全| 国产一区二区激情短视频| 久久久精品欧美日韩精品| 亚洲狠狠婷婷综合久久图片| 国模一区二区三区四区视频| 成人性生交大片免费视频hd| 欧美一级毛片孕妇| 桃色一区二区三区在线观看| 精品日产1卡2卡| 脱女人内裤的视频| 国产不卡一卡二| 好男人在线观看高清免费视频| 18+在线观看网站| 啦啦啦免费观看视频1| 男人的好看免费观看在线视频| a级一级毛片免费在线观看| 五月伊人婷婷丁香| 久久久色成人| 欧美最新免费一区二区三区 | 叶爱在线成人免费视频播放| 看免费av毛片| 黄色女人牲交| 女人被狂操c到高潮| 老司机深夜福利视频在线观看| 国产91精品成人一区二区三区| 少妇人妻精品综合一区二区 | 9191精品国产免费久久| 国产私拍福利视频在线观看| 两人在一起打扑克的视频| 好看av亚洲va欧美ⅴa在| 亚洲国产精品sss在线观看| 国产精华一区二区三区| 免费看光身美女| 黄片小视频在线播放| 亚洲成a人片在线一区二区| 国产视频一区二区在线看| 国产三级黄色录像| 最近最新中文字幕大全电影3| 亚洲av成人不卡在线观看播放网| 午夜福利视频1000在线观看| www日本在线高清视频| 中文字幕人妻丝袜一区二区| 日本黄色片子视频| 在线观看免费视频日本深夜| 51午夜福利影视在线观看| 精品久久久久久久久久免费视频| eeuss影院久久| 熟女人妻精品中文字幕| 在线观看66精品国产| 黄色视频,在线免费观看| xxx96com| 久久这里只有精品中国| 久久久国产成人精品二区| 12—13女人毛片做爰片一| av视频在线观看入口| 波多野结衣巨乳人妻| 久久精品国产亚洲av香蕉五月| 少妇人妻精品综合一区二区 | 亚洲精品亚洲一区二区| 麻豆国产av国片精品| 久久久久国产精品人妻aⅴ院| 岛国视频午夜一区免费看| 亚洲第一欧美日韩一区二区三区| 亚洲av成人不卡在线观看播放网| 麻豆成人午夜福利视频| 动漫黄色视频在线观看| 一个人观看的视频www高清免费观看| 日日摸夜夜添夜夜添小说| 热99在线观看视频| 成年女人毛片免费观看观看9| 国产99白浆流出| av视频在线观看入口| 亚洲中文日韩欧美视频| 亚洲18禁久久av| 色av中文字幕| 亚洲成人精品中文字幕电影| 热99re8久久精品国产| 国产视频一区二区在线看| 狠狠狠狠99中文字幕| 亚洲av中文字字幕乱码综合| 天堂动漫精品| 午夜久久久久精精品| 久久国产乱子伦精品免费另类| 麻豆国产av国片精品| 三级男女做爰猛烈吃奶摸视频| 无遮挡黄片免费观看| 少妇丰满av| 三级国产精品欧美在线观看| 国产精品久久久久久亚洲av鲁大| 精华霜和精华液先用哪个| 欧美bdsm另类| 制服丝袜大香蕉在线| 无人区码免费观看不卡| 精品一区二区三区av网在线观看| 国产精品亚洲一级av第二区| 亚洲精品美女久久久久99蜜臀| 狂野欧美白嫩少妇大欣赏| 国产精品99久久久久久久久| 搡女人真爽免费视频火全软件 | 极品教师在线免费播放| 成年女人毛片免费观看观看9| 精品国产美女av久久久久小说| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播| 日本成人三级电影网站| 久久亚洲精品不卡| 男女之事视频高清在线观看| 国内久久婷婷六月综合欲色啪| 日韩免费av在线播放| 午夜亚洲福利在线播放| 日韩欧美国产在线观看| 国产三级中文精品| 两个人视频免费观看高清| 美女被艹到高潮喷水动态| 18禁国产床啪视频网站| 精品人妻一区二区三区麻豆 | 亚洲成人久久爱视频| 国产精品久久久久久久电影 | 欧美一区二区精品小视频在线| 美女cb高潮喷水在线观看| 一个人看的www免费观看视频| 我要搜黄色片| 亚洲va日本ⅴa欧美va伊人久久| a级毛片a级免费在线| xxx96com| 老汉色av国产亚洲站长工具| 成人欧美大片| 亚洲精品国产精品久久久不卡| 日本一本二区三区精品| 在线免费观看的www视频| 亚洲精品乱码久久久v下载方式 | 真人做人爱边吃奶动态| 久久香蕉国产精品| 国产伦精品一区二区三区四那| 亚洲精品美女久久久久99蜜臀| 日本一本二区三区精品| 国产主播在线观看一区二区| 国产精品永久免费网站| 男女做爰动态图高潮gif福利片| 亚洲avbb在线观看| 日本 av在线| 人妻久久中文字幕网| 神马国产精品三级电影在线观看| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼| 午夜福利在线观看吧| 久久久久免费精品人妻一区二区| 琪琪午夜伦伦电影理论片6080| 日本一本二区三区精品| 免费高清视频大片| 真人做人爱边吃奶动态| 久久久久久国产a免费观看| 日日夜夜操网爽| 少妇的逼好多水| 最近最新中文字幕大全电影3| 香蕉av资源在线| 精品午夜福利视频在线观看一区| 色在线成人网| 色视频www国产| 午夜福利在线在线| 亚洲成a人片在线一区二区| 国产av在哪里看| x7x7x7水蜜桃| 国产成人系列免费观看| 最新在线观看一区二区三区| 91久久精品电影网| 免费在线观看亚洲国产| 精品久久久久久,| 国产中年淑女户外野战色| 中文字幕av成人在线电影| 国产亚洲精品久久久com| 又紧又爽又黄一区二区| 欧美乱妇无乱码| 精品福利观看| 淫秽高清视频在线观看| 美女 人体艺术 gogo| 国产精品一区二区三区四区久久| 国产免费一级a男人的天堂| 一区福利在线观看| 女生性感内裤真人,穿戴方法视频| www国产在线视频色| h日本视频在线播放| 国产精品日韩av在线免费观看| 日本撒尿小便嘘嘘汇集6| 在线十欧美十亚洲十日本专区| 国产一区二区在线av高清观看| 亚洲第一电影网av| 成年女人看的毛片在线观看| 亚洲性夜色夜夜综合| 中文字幕av在线有码专区| 欧美黑人欧美精品刺激| 精品午夜福利视频在线观看一区| 亚洲一区二区三区色噜噜| 深爱激情五月婷婷| 成人欧美大片| www.色视频.com| www国产在线视频色| 国产av在哪里看| 51国产日韩欧美| 午夜日韩欧美国产| 成人三级黄色视频| 无人区码免费观看不卡| 国产极品精品免费视频能看的| 午夜亚洲福利在线播放| 最近最新中文字幕大全电影3| 国产男靠女视频免费网站| 中国美女看黄片| 18+在线观看网站| 精品久久久久久久毛片微露脸| 亚洲av中文字字幕乱码综合| 亚洲内射少妇av| 日韩亚洲欧美综合| 99久久99久久久精品蜜桃| 天堂√8在线中文| 校园春色视频在线观看| 国产真实伦视频高清在线观看 | 动漫黄色视频在线观看| 国产精品精品国产色婷婷| 丰满人妻熟妇乱又伦精品不卡| 午夜福利18| 在线观看免费午夜福利视频| 国产爱豆传媒在线观看| 1024手机看黄色片| 国产精品一区二区三区四区免费观看 | 精品福利观看| 国产精品久久久久久人妻精品电影| 老司机福利观看| 欧美激情久久久久久爽电影| 日韩欧美一区二区三区在线观看| 无人区码免费观看不卡| 麻豆久久精品国产亚洲av| 一区二区三区免费毛片| 午夜福利高清视频| 日本 av在线| 首页视频小说图片口味搜索| 老汉色∧v一级毛片| 欧美zozozo另类| 中文字幕av在线有码专区| 亚洲,欧美精品.| 精品人妻一区二区三区麻豆 | 国产真实伦视频高清在线观看 | 综合色av麻豆| 国产成人av教育| 两个人的视频大全免费| 日本一二三区视频观看| av福利片在线观看| 免费大片18禁| 欧美大码av| 久久久久久久久大av| 亚洲精品乱码久久久v下载方式 | 51国产日韩欧美| 欧美日韩黄片免| 最近最新免费中文字幕在线| 亚洲国产欧美网| 长腿黑丝高跟| 啦啦啦韩国在线观看视频| 在线免费观看的www视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲自拍偷在线| 国产美女午夜福利| 色哟哟哟哟哟哟| 禁无遮挡网站| 啦啦啦韩国在线观看视频| 精品国产三级普通话版| 可以在线观看的亚洲视频| 国产探花极品一区二区| 国产一区二区激情短视频| 免费看日本二区| 国产视频一区二区在线看| 床上黄色一级片| 国产精品一区二区免费欧美| 国产真人三级小视频在线观看| 无人区码免费观看不卡| 成年版毛片免费区| 18禁在线播放成人免费| 国产v大片淫在线免费观看| 欧美一区二区精品小视频在线| 久9热在线精品视频| 嫩草影视91久久| 99久国产av精品| 国产91精品成人一区二区三区| 国产探花在线观看一区二区| 搡老妇女老女人老熟妇| 久久香蕉精品热| 女警被强在线播放| 成人三级黄色视频| 亚洲av成人精品一区久久| 九九热线精品视视频播放| 成人国产综合亚洲| 久久久国产成人免费| 亚洲成人精品中文字幕电影| 国产 一区 欧美 日韩| 久久久久久久午夜电影| 亚洲国产精品合色在线| 久9热在线精品视频| 欧美av亚洲av综合av国产av| 久久久久久久亚洲中文字幕 | 熟女人妻精品中文字幕| 国产真人三级小视频在线观看| 亚洲天堂国产精品一区在线| 久久精品综合一区二区三区| 欧美+日韩+精品| 一本久久中文字幕| 欧美又色又爽又黄视频| 内地一区二区视频在线| 少妇的逼好多水| 在线十欧美十亚洲十日本专区| 久久九九热精品免费| 亚洲成人精品中文字幕电影| 两个人看的免费小视频| 最近视频中文字幕2019在线8| 亚洲色图av天堂| 麻豆国产av国片精品| 国产男靠女视频免费网站| 精品久久久久久久久久久久久| 欧美日韩中文字幕国产精品一区二区三区| www日本在线高清视频| 真人一进一出gif抽搐免费| 亚洲国产精品sss在线观看| 国产色爽女视频免费观看| 久久精品综合一区二区三区| 国产精品98久久久久久宅男小说| 黄色片一级片一级黄色片| 色噜噜av男人的天堂激情| 日日夜夜操网爽| 免费人成视频x8x8入口观看| 大型黄色视频在线免费观看| 99热这里只有精品一区| 亚洲精品粉嫩美女一区| 日韩精品中文字幕看吧| 久久久国产成人免费| 欧美bdsm另类| 亚洲av中文字字幕乱码综合| 精品久久久久久久末码| 日韩中文字幕欧美一区二区| 国产中年淑女户外野战色| 他把我摸到了高潮在线观看| 91九色精品人成在线观看| 日本一二三区视频观看| 美女大奶头视频| 啦啦啦观看免费观看视频高清| 欧美日韩乱码在线| 久久香蕉精品热| 国产精品亚洲av一区麻豆| 日韩欧美三级三区| 天天一区二区日本电影三级| 好男人在线观看高清免费视频| eeuss影院久久| 看片在线看免费视频| 一区二区三区免费毛片| 日本与韩国留学比较| 日本一本二区三区精品| 狠狠狠狠99中文字幕| 国产高潮美女av| 国产久久久一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| av天堂中文字幕网| 国产97色在线日韩免费| 日韩欧美三级三区| 757午夜福利合集在线观看| 欧美日韩福利视频一区二区| 免费在线观看影片大全网站| 法律面前人人平等表现在哪些方面| 淫秽高清视频在线观看| 亚洲国产高清在线一区二区三| 老司机午夜福利在线观看视频| 波多野结衣高清作品| 99久久九九国产精品国产免费| 亚洲人成电影免费在线| 国产高清激情床上av| 精品久久久久久久久久免费视频| 日本免费a在线| 老汉色∧v一级毛片| 国产成人啪精品午夜网站| 天天添夜夜摸| 中文亚洲av片在线观看爽| 国产成人欧美在线观看| 国产欧美日韩精品亚洲av| 国语自产精品视频在线第100页| 在线看三级毛片| 精品乱码久久久久久99久播| 两性午夜刺激爽爽歪歪视频在线观看| 真人一进一出gif抽搐免费| 国产高清激情床上av| 好看av亚洲va欧美ⅴa在| 欧美乱妇无乱码| 精品国产三级普通话版| 91九色精品人成在线观看| 1024手机看黄色片| 九九久久精品国产亚洲av麻豆| 国产熟女xx| 亚洲成av人片免费观看| 久久亚洲精品不卡| 欧美日韩国产亚洲二区| 久久午夜亚洲精品久久| 国产男靠女视频免费网站| 欧美国产日韩亚洲一区| 亚洲精品美女久久久久99蜜臀| 亚洲无线观看免费| 天美传媒精品一区二区| 三级男女做爰猛烈吃奶摸视频| 亚洲国产色片| 亚洲精品亚洲一区二区| 亚洲第一欧美日韩一区二区三区| 久久精品国产亚洲av香蕉五月| 真实男女啪啪啪动态图|