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

    基于斜面輻射算法的短期光伏功率預(yù)測方法研究

    2020-11-09 07:36:50李照榮王小勇閆曉敏趙文婧
    干旱氣象 2020年5期
    關(guān)鍵詞:實況斜面線性

    李 遙,李照榮,王小勇,閆曉敏,趙文婧

    (1.甘肅省氣象服務(wù)中心,甘肅 蘭州 730020;2.甘肅省氣象局,甘肅 蘭州 730020)

    引 言

    太陽能作為一種可再生清潔能源,其開發(fā)與利用已成為減緩全球氣候變暖、實現(xiàn)我國能源戰(zhàn)略轉(zhuǎn)型的重要舉措之一,光伏發(fā)電是太陽能利用最主要的方式之一,但由于其波動性強、隨機性大等特點,大規(guī)模并網(wǎng)勢必會對電網(wǎng)系統(tǒng)造成嚴重沖擊[1-2],因此,準確預(yù)測光伏發(fā)電功率是提高太陽能資源利用率的關(guān)鍵問題。

    光伏功率預(yù)測從預(yù)測時間上主要分為超短期(未來4 h)、短期(未來72 h)以及中長期(未來168 h)預(yù)測[3],短期預(yù)測是電力系統(tǒng)考核的重要指標之一,直接影響電站當(dāng)日的發(fā)電計劃,因此,短期光伏功率的精準預(yù)測尤為重要。從預(yù)測方法上主要分為原理法和統(tǒng)計法。原理法是基于輻射預(yù)測和光伏發(fā)電物理原理進行功率預(yù)測的一種間接預(yù)測方法[4],該方法計算原理簡單,不需要依靠歷史數(shù)據(jù)建立預(yù)測模型,但是對電站系統(tǒng)數(shù)據(jù)的精度要求較高,包括光伏逆變器安裝的地理位置以及光電轉(zhuǎn)換效率等信息[5]。統(tǒng)計法基于歷史數(shù)據(jù)資料,通過統(tǒng)計學(xué)方法建立預(yù)測模型[4],主要有人工神經(jīng)網(wǎng)絡(luò)[6]、回歸分析[7]、支持向量機[8]等方法,然后輸入數(shù)值預(yù)報數(shù)據(jù)得到預(yù)測結(jié)果,是一種直接預(yù)測方法。該方法不需要考慮光伏發(fā)電原理以及詳細的電站系統(tǒng)信息數(shù)據(jù),但是對電站的歷史數(shù)據(jù)質(zhì)量要求較高,增加了數(shù)據(jù)前期收集和預(yù)處理工作的難度[5]。

    近年來,隨著我國光伏發(fā)電產(chǎn)業(yè)的迅速發(fā)展和裝機容量的不斷增長,光伏功率預(yù)測技術(shù)取得了重要進展。華中科技大學(xué)及湖北省氣象服務(wù)中心先后建立了光伏發(fā)電實驗平臺和示范光伏電站,積累了可靠的長期實驗數(shù)據(jù)資料,獲得一系列相關(guān)研究成果,使得光伏發(fā)電的隨機化問題得到改善[9],并針對不同電站運行情況建立多種短期功率預(yù)測模型以滿足實際業(yè)務(wù)需求[10]。此外,多種數(shù)學(xué)模型及算法組合的預(yù)測方法相繼出現(xiàn),其預(yù)測精度往往較傳統(tǒng)單一模型有不同程度的提高。葛樂等[11]基于改進相似日和人工蜂群算法對原有支持向量機模型進行了優(yōu)化。其他類似的組合預(yù)測方法還有GA-BP算法[12]、IWD-Elman神經(jīng)網(wǎng)絡(luò)模型[13]和熵權(quán)法[14]等,均使得功率預(yù)測效果有所改進。

    甘肅省太陽能資源豐富,年太陽總輻射量為4700~6350 MJ·m-2[15],具有優(yōu)越的光伏產(chǎn)業(yè)發(fā)展優(yōu)勢,近年來,該省光伏市場規(guī)模增長迅速,根據(jù)國家能源局統(tǒng)計結(jié)果,截止2018年底光伏發(fā)電累計裝機容量達到8 280 000 kW。然而,對于該省太陽能光伏發(fā)電預(yù)測技術(shù)的研究相對較少,預(yù)測精度還有待進一步提高。基于此,本文選取甘肅北部典型光伏電站,在改進原有斜面輻射算法的基礎(chǔ)上,基于多元線性回歸法和經(jīng)驗公式法,建立適用于光伏電站實際業(yè)務(wù)化運行的短期光伏發(fā)電功率預(yù)測模型,通過對結(jié)果的檢驗和評估,為電站選取更優(yōu)的預(yù)測模型提供參考。

    1 資料來源

    選取2017年3月至2019年2月中電蘆陽扶貧光伏電站(ZDLYFP,103.58°E、37.23°N)15 min間隔的功率與環(huán)境觀測數(shù)據(jù),其中環(huán)境數(shù)據(jù)為光伏電站內(nèi)氣象站觀測資料,包括總輻射、風(fēng)速、風(fēng)向、溫度、相對濕度和氣壓。該電站裝機5000 kW,位于甘肅省白銀市景泰縣蘆陽鎮(zhèn),地處黃土高原與騰格里沙漠過渡地帶,地勢較平坦,年日照時數(shù)約2726 h,年太陽總輻射量約 5870 MJ·m-2[16],屬于太陽能資源較豐富地區(qū)。該站為扶貧光伏電站,不存在限電情況,且運行狀況良好,設(shè)備維護較為及時,數(shù)據(jù)質(zhì)量較高,適宜作為研究站點。

    另外,選取同期中國氣象局公共氣象服務(wù)中心提供的BJ-RUC數(shù)值預(yù)報數(shù)據(jù),包括總輻射、直接輻射、散射輻射、風(fēng)速、風(fēng)向、溫度、相對濕度和氣壓,預(yù)報時間長度為72 h,時間分辨率為15 min。

    經(jīng)過篩選,剔除停電、檢修時段以及儀器故障引起的異常值,有效數(shù)據(jù)占總數(shù)據(jù)量的95.6%,數(shù)據(jù)資料相對較為完整。其中,2017年3月至2018年2月的數(shù)據(jù)資料用于前期數(shù)據(jù)分析,2018年3月至2019年2月的數(shù)據(jù)資料用于短期光伏功率預(yù)測方法研究和結(jié)果分析與檢驗。

    2 方 法

    2.1 斜面輻射算法的改進

    為提高太陽能資源的利用效率,電站的光伏面板一般傾斜放置[17],傾斜面板接收到的太陽輻射由直接輻射、散射輻射和反射輻射三部分組成[18],斜面總輻射 IT的計算公式為[4,19]:

    式中:IT,b、IT,d、IT,g(W·m-2)分別表示斜面接收到的直接輻射、散射輻射和地面反射輻射。

    利用Liu-Jordan模型進行斜面輻射計算,該模型是典型的天空各向同性模型,假設(shè)天空中太陽散射輻射各向同性且均勻分布,基于該模型計算IT的表達式為[20]:

    式中:Ib、Id、I(W·m-2)分別表示水平面直接輻射、散射輻射和總輻射;ρ為地表反照率,一般取0.2;Rb為 IT,b與 Ib的比值;β為斜面傾角。

    Rb通過以下幾何關(guān)系計算得到[4,21]:

    式中:δ=23.45sin[360×(284+n)/365],δ為當(dāng)天赤緯,n為當(dāng)日在一年中的日序號;γT為光伏面板斜面方位角;φ為緯度;ω為時角。

    在進行光伏功率預(yù)測時,利用水平面直接輻射、散射輻射和總輻射的模式預(yù)報數(shù)據(jù),根據(jù)公式(2)可計算斜面輻射預(yù)報值進而實現(xiàn)功率預(yù)報。然而,在實際業(yè)務(wù)中由此計算得到的預(yù)測值往往與實際功率存在差距,預(yù)報效果不理想,因此,本文利用直散分離模型對原有的計算方法進行改進。直散分離模型將水平面總輻射分解為水平面直接輻射和散射輻射,其中較為典型的直散分離模型有Erbs模型[22]、DIRINT模型[23]、Reindl模型[24]等,本文選取 Erbs模型進行直散分離計算。

    Erbs模型因其準確率高且參數(shù)較易獲取而廣泛應(yīng)用,該模型是利用晴空指數(shù)KT反演散射分數(shù)DF的分段多項式模型[19],晴空指數(shù)KT表征大氣的通透程度,是水平面總輻射與大氣層外太陽輻射的比值,散射分數(shù)DF為水平面散射輻射與總輻射的比值,Erbs模型的表達式為[22]:

    大氣層外太陽輻射I0的計算公式為[25]:

    式中:Esc為太陽常數(shù),取1367 W·m-2;γ為日地距離訂正系數(shù),γ=1+0.033cos(360×n/365)。

    基于Erbs模型得到Ib和Id的預(yù)報值,再根據(jù)Liu-Jordan模型得到斜面輻射的預(yù)報結(jié)果。

    2.2 短期光伏功率預(yù)測方法

    基于上述改進后的斜面總輻射計算方法和BJ-RUC資料,采用多元線性回歸法[4]和經(jīng)驗公式法分別建立短期光伏功率預(yù)測模型。光伏電站每日需要上傳未來3 d的短期功率預(yù)測結(jié)果,但考核時只針對未來第一天的預(yù)測結(jié)果,因此僅選取第一天的模型輸出結(jié)果進行檢驗。

    2.2.1 多元線性回歸法

    為選取對發(fā)電功率影響最為顯著的預(yù)報因子,利用2017年3月至2018年2月實況功率與BJ-RUC數(shù)值預(yù)報數(shù)據(jù),先將總輻射換算成斜面總輻射,然后對實況功率與各預(yù)報變量分季節(jié)進行相關(guān)性分析。從表1相關(guān)分析結(jié)果看出,實況功率與斜面總輻射、風(fēng)速、溫度呈正相關(guān)關(guān)系,而與風(fēng)向、相對濕度、氣壓呈負相關(guān)關(guān)系,除氣壓外,其他氣象要素各季節(jié)的相關(guān)性均通過0.01的顯著性檢驗。其中,實況功率與斜面總輻射相關(guān)系數(shù)各季節(jié)均遠大于其他要素,冬季最小為0.862,春季最大為0.924;溫度次之,夏季溫度對發(fā)電功率的影響最大,相關(guān)系數(shù)為0.471。另外,實況功率與斜面總輻射和溫度總的相關(guān)系數(shù)分別為0.896和0.386,明顯大于其他要素,考慮到預(yù)報因子對光伏發(fā)電功率的物理意義,最終確定斜面總輻射和溫度作為構(gòu)建多元線性回歸模型的預(yù)報因子。

    表1 不同季節(jié)各因子與實況功率的相關(guān)系數(shù)Tab.1 Correlation coefficients between real power and each factor in different seasons

    采用滾動系數(shù)法建立預(yù)測模型,即利用過去一段時間的實況和預(yù)報資料進行回歸分析,建立二者之間的動態(tài)統(tǒng)計模型,采用滾動樣本的方式每天對模型參數(shù)進行更新,考慮到季節(jié)氣候變化帶來的輻射、溫度變化,滾動系數(shù)法建模一般使用20 d左右的樣本量較為合理[4,10]。為選取具有代表性的樣本量,首先對不同樣本量建模的預(yù)測效果進行評估,發(fā)現(xiàn)15 d的均方根誤差最小,20 d與15 d相差不大,而10 d和25 d以上均方根誤差較大且比較接近(表2),因此建模時選擇15 d的樣本容量。滾動過程:第i天進行第i+1天功率預(yù)測時,利用第i-15至第i-1天的BJ-RUC數(shù)值預(yù)報數(shù)據(jù)和實況功率數(shù)據(jù)建立模型,結(jié)合第i天生成的BJ-RUC預(yù)報數(shù)據(jù)即可獲得第i+1天逐15 min的預(yù)測功率,并利用第i+1天的實況功率檢驗預(yù)測效果。該方法可以減少季節(jié)氣候變化對光伏發(fā)電的影響以及光伏電站長期運營而產(chǎn)生的衰減效應(yīng)。

    為評估多元線性回歸方程的預(yù)測效果,對構(gòu)建的方程進行擬合優(yōu)度檢驗、方程整體及各項回歸系數(shù)的顯著性檢驗,發(fā)現(xiàn)每日建立的多元線性回歸方程整體及斜面總輻射項均通過顯著性檢驗,溫度項t檢驗的Sig.值小于0.05的天數(shù)占總數(shù)的83.4%(表略),大部分天數(shù)溫度對功率的影響顯著,可見斜面總輻射、溫度與功率的擬合效果比較理想,回歸方程較為可靠。以2019年3月16日的顯著性檢驗結(jié)果為例,發(fā)現(xiàn)方程的擬合優(yōu)度 R2為0.926,接近1,F(xiàn)檢驗的 Sig.值接近0,方程整體顯著,且各項回歸系數(shù)t檢驗的 Sig.值均小于0.05(表3),說明回歸方程中各項對功率的影響均顯著。

    表2 不同樣本量建模的預(yù)測值與實況的均方根誤差Tab.2 The root mean square error between real value and forecasted value by models based on different sample sizes

    表3 各項回歸系數(shù)顯著性檢驗Tab.3 Significance test of regression coefficients of each item

    2.2.2 經(jīng)驗公式法

    經(jīng)驗公式法是根據(jù)太陽能電池光伏發(fā)電的物理原理、光電轉(zhuǎn)換效率和逆變器轉(zhuǎn)換效率的定義,建立影響光電轉(zhuǎn)換效率及逆變器效率的經(jīng)驗公式和合理的經(jīng)驗系數(shù),輸入輻射預(yù)報,進行功率預(yù)測的方法[10]。基于經(jīng)驗公式法建立的光伏電站輸出功率的預(yù)測模型為[4]:

    式中:IT(W·m-2)為斜面總輻射;NOCT(℃)為額定太陽能電池工作溫度;Ta(℃)為氣溫;Tc(℃)為陣列板溫度;Pdc(kW)為15 min間隔的發(fā)電功率;ηs為標準測試條件的光電轉(zhuǎn)換效率;α(℃-1)為溫度系數(shù);S(m2)為光伏組件有效面積;K1為光伏陣列老化損失系數(shù);K2為光伏陣列失配損失系數(shù);K3為塵埃遮擋損失系數(shù);K4為直流回路線路損失系數(shù);k為并網(wǎng)光伏電站投入使用的年數(shù);ya為不同太陽能電池材料年衰減率。

    利用相關(guān)系數(shù)(CC)、均方根誤差(RMSE)和相對均方根誤差(RRMSE)等評價指標[4],對斜面輻射算法改進后多元線性回歸法、經(jīng)驗公式法的功率預(yù)測效果進行檢驗及對比,兩者的預(yù)測功率分別用P1、P2表示,實況功率用Pr表示。

    3 結(jié)果分析

    3.1 斜面輻射算法改進后的效果檢驗

    基于改進后的斜面輻射算法,利用多元線性回歸法和經(jīng)驗公式法分別計算2018年3月至2019年2月逐15 min的短期功率。對比算法改進前后兩種模型預(yù)測功率的RRMSE(圖1),發(fā)現(xiàn)各月RRMSE在算法改進后均有不同程度的降低,P1和P2平均RRMSE分別降低0.066和0.040,斜面輻射算法的改進使得功率預(yù)測效果有所提高。

    圖1 斜面輻射算法改進前后RRMSE對比(a)多元線性回歸法,(b)經(jīng)驗公式法Fig.1 Comparison of RRMSE before and after the improvement of inclined plane solar radiation method(a)multiple linear regression method,(b)empirical formula method

    3.2 兩種方法預(yù)測結(jié)果對比

    圖2是多元線性回歸法和經(jīng)驗公式法模擬功率與實況功率的日變化??梢钥闯觯骷竟?jié)Pr與P1、P2的日變化特征均呈單峰結(jié)構(gòu),變化趨勢較為一致。春、夏季節(jié)由于日照時數(shù)相對較長,曲線開口程度大于秋、冬季節(jié),光伏電站日平均工作時間較長,尤其是夏季最長,近14 h,冬季日照時數(shù)最短,曲線開口程度最小,約10 h;Pr春季峰值最大,秋季次之,夏季最低,夏季峰值低于其他季節(jié)是由于日間受高溫影響,光伏面板輸出功率的效率降低,導(dǎo)致實際發(fā)電功率低于其他季節(jié)。P1、P2與Pr的偏差表現(xiàn)為正午前后最大,向兩側(cè)呈現(xiàn)逐漸減小的趨勢,且各季節(jié)偏差程度存在差異,P1秋季偏差最小,日最大偏差為378.8 kW,春季最大,日最大偏差為888.7 kW,這是由于春季天氣過程較多,輻射變化較為劇烈,預(yù)報輻射未能準確預(yù)測實際輻射的變化特征,導(dǎo)致功率預(yù)測偏差較大;P2冬季偏差最小,日最大偏差為863.0 kW,春季最大,日最大偏差為1197.3 kW,且各季節(jié)P2的偏差程度較P1更為明顯。

    3.3 兩種方法預(yù)測誤差對比

    對比兩種方法的功率預(yù)測效果(表4)發(fā)現(xiàn),P1和P2與實況功率的相關(guān)系數(shù)各月均較為接近,且均通過0.01的顯著性檢驗,總的相關(guān)系數(shù)分別為0.823和 0.817,P1與實況功率的相關(guān)性略高于 P2;P1與實況功率的RMSE和RRMSE各月均小于P2,P1和P2總RMSE分別為940.917、1147.172 kW,總RRMSE分別為0.188和0.229,多元線性回歸法的預(yù)報效果優(yōu)于經(jīng)驗公式法。

    圖3是兩種方法RRMSE的月變化??梢钥闯?,二者全年均呈現(xiàn)波動變化趨勢,P1的RRMSE在4—7月以及10—12月較小且較為穩(wěn)定,而在3月、8月和1月出現(xiàn)明顯的峰值,這是由于在季節(jié)交替時,輻射及溫度變化較為劇烈,導(dǎo)致預(yù)測結(jié)果誤差增大;P2的RRMSE在9、11月以及2月出現(xiàn)明顯的峰值,P2與P1差值最小和最大的月份分別為3月和11月,差值分別為0.009和0.067,且 P2的變化幅度大于P1,P2各月之間的差異更為顯著,說明多元線性回歸法的預(yù)報效果較經(jīng)驗公式法更為穩(wěn)定。

    圖2 多元線性回歸法、經(jīng)驗公式法的模擬功率與實況功率的日變化(a)春季,(b)夏季,(c)秋季,(d)冬季Fig.2 Diurnal change of simulated power by MLR method and empirical formula and actual power(a)spring,(b)summer,(c)autumn,(d)winter

    表4 兩種方法預(yù)報結(jié)果檢驗Tab.4 Test of forecasted results based on two methods

    圖3 兩種方法預(yù)報結(jié)果與實況功率的RRMSE月變化Fig.3 Monthly change of RRMSE between actual power and predicted power by two methods

    3.4 不同天氣條件下兩種方法的預(yù)報效果對比

    將2018年3月至2019年2月按天氣條件進行分類,可分為晴天、多云、陰、雨、雪、沙塵(以及污染)等6種天氣類型,對比不同天氣條件下兩種方法的預(yù)報效果(圖4)。可以看出,兩種方法預(yù)測功率與實況功率的RRMSE值均呈現(xiàn)自晴天、多云天、陰天、雨天、沙塵天、雪天依次增大的趨勢。晴天時,P1與P2的預(yù)報效果最好,RRMSE分別為0.112和0.143,預(yù)報效果優(yōu)于其他天氣狀況;多云天氣次之,RRMSE分別為0.188和0.228,這兩種天氣狀況占所有天數(shù)的71.2%。陰天和雨天所占比例為21.9%,二者預(yù)報效果相當(dāng),RRMSE值均大于晴天及多云天氣;雪天和沙塵天氣預(yù)報效果最不理想,沙塵天氣下P1與 P2的 RRMSE值分別為 0.273和 0.337,雪天分別達0.277和0.341,遠高于晴天。

    為了進一步對比兩種方法的預(yù)報效果,選取晴天(2018年10月9日)、多云天(2018年5月27日)、陰天(2019年1月29日)、雨天(2018年8月7日)、雪天(2019年1月15日)、沙塵天(2018年4月4日)個例,細致分析不同天氣狀況下兩種方法預(yù)測功率與實況的差異(圖5)??梢钥闯觯缣鞎r,兩種方法均能準確預(yù)測實際發(fā)電功率的變化趨勢,但P1與實況曲線基本重合,擬合程度最好;在多云及陰、雨天氣狀況時,P1與P2能大致預(yù)測出實際功率的變化趨勢和變化幅度,但無法準確把握實際功率的波動特征;沙塵、雪天氣時,預(yù)測功率未能捕捉到實際功率的變化趨勢和變化幅度,預(yù)測效果不理想。

    在數(shù)值預(yù)報相同的情況下,P1在各種天氣條件下RRMSE值均小于P2,預(yù)報效果相對更佳。在有天氣過程時,由于數(shù)值預(yù)報與實況存在偏差,導(dǎo)致P1和P2均未能準確捕捉到發(fā)電功率的實際變化特征。下一步考慮結(jié)合衛(wèi)星云圖、沙塵天氣顆粒物濃度等數(shù)據(jù)資料,深入研究有天氣過程時光伏發(fā)電功率的日變化特征,同時考慮使用集合預(yù)報集成不同方法的優(yōu)勢提高功率預(yù)測準確率。

    圖4 不同天氣占比(a)及不同天氣狀況下兩種方法預(yù)測效果(b)Fig.4 The percentage of different weather conditions(a)and forecast effect based on two methods under different weather conditions(b)

    圖5 不同天氣條件下兩種方法預(yù)測功率與實況的日變化(a)晴天,(b)多云天,(c)陰天,(d)雨天,(e)雪天,(f)沙塵天Fig.5 Diurnal change of predicted power by two methods and actual power under different weather conditions(a)sunny weather,(b)cloudy weather,(c)overcast weather,(d)rainy weather,(e)snow weather,(f)dust weather

    4 結(jié) 論

    (1)各季節(jié)斜面總輻射和溫度與光伏發(fā)電功率相關(guān)系數(shù)明顯高于其他變量,總相關(guān)系數(shù)分別為0.896和0.386,可作為基于多元線性回歸方法的發(fā)電功率預(yù)測模型的預(yù)報因子。

    (2)斜面輻射計算方法改進后,多元線性回歸法和經(jīng)驗公式法的預(yù)測功率與實況各月RRMSE值均有所降低,總RRMSE值分別降低0.066和0.040,改進后兩種方法的預(yù)報效果均有所提高。

    (3)多元線性回歸法與經(jīng)驗公式法的預(yù)測功率日變化曲線均呈單峰結(jié)構(gòu),峰值與開口程度季節(jié)差異明顯,二者總體與實況功率變化趨勢較為一致,但存在不同程度的偏差,多元線性回歸法的偏差小于經(jīng)驗公式法,更接近實況功率。

    (4)多元線性回歸法各月RMSE、RRMSE均小于經(jīng)驗公式法,CC均略高于經(jīng)驗公式法,兩種方法總RRMSE值分別為0.188和0.229,總 CC值分別為0.823和0.817,且前者RRMSE月變化幅度小于后者,表明多元線性回歸法的預(yù)報精度更高,且預(yù)報結(jié)果更為穩(wěn)定。

    (5)不同天氣狀況下,多元線性回歸法和經(jīng)驗公式法的RRMSE值均自晴天、多云天、陰天、雨天、沙塵天、雪天依次增大,兩種方法對晴天的預(yù)測效果均最優(yōu),RRMSE分別為0.112和0.143,而對雪天和沙塵天氣的預(yù)測效果均不理想,雪天RRMSE值最大,分別為0.277和0.341。各種天氣條件下,多元線性回歸法的RRMSE值均小于經(jīng)驗公式法。

    猜你喜歡
    實況斜面線性
    漸近線性Klein-Gordon-Maxwell系統(tǒng)正解的存在性
    斜面之上探動能
    巧用“相對”求解光滑斜面體問題
    巧用“相對”求解光滑斜面體問題
    線性回歸方程的求解與應(yīng)用
    鄉(xiāng)村小學(xué)的愿望與現(xiàn)實——宜君鄉(xiāng)村教育實況
    二階線性微分方程的解法
    天舟一號貨運飛船發(fā)射實況掠影
    可愛潮咖們的獨門彩妝實況直播
    一題多變 搞定斜面上的運動
    在线看a的网站| 国产高清国产精品国产三级| 最后的刺客免费高清国语| 国产无遮挡羞羞视频在线观看| 一区二区日韩欧美中文字幕 | 中国三级夫妇交换| 午夜91福利影院| 超碰97精品在线观看| 80岁老熟妇乱子伦牲交| 人成视频在线观看免费观看| 免费不卡的大黄色大毛片视频在线观看| 毛片一级片免费看久久久久| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人精品一区二区| 免费人妻精品一区二区三区视频| 久久久久人妻精品一区果冻| 狠狠精品人妻久久久久久综合| 午夜免费观看性视频| 国产免费福利视频在线观看| 日韩av免费高清视频| 卡戴珊不雅视频在线播放| 99热这里只有是精品在线观看| 久久久精品免费免费高清| 九九在线视频观看精品| 十分钟在线观看高清视频www| 国产精品久久久久久av不卡| 久久精品国产a三级三级三级| 黄色一级大片看看| 少妇被粗大猛烈的视频| 9色porny在线观看| 亚洲成人一二三区av| 欧美激情 高清一区二区三区| 天天躁夜夜躁狠狠久久av| 大片电影免费在线观看免费| 欧美xxxx性猛交bbbb| 亚洲丝袜综合中文字幕| 丝袜人妻中文字幕| 久久人人爽av亚洲精品天堂| 韩国av在线不卡| 午夜视频国产福利| 女人久久www免费人成看片| 在线天堂最新版资源| 91aial.com中文字幕在线观看| 91久久精品国产一区二区三区| 欧美日韩av久久| 日韩电影二区| 寂寞人妻少妇视频99o| 久久久久视频综合| 亚洲av免费高清在线观看| 国产色爽女视频免费观看| 熟女电影av网| 我的女老师完整版在线观看| 在现免费观看毛片| 日韩精品免费视频一区二区三区 | 日日爽夜夜爽网站| 亚洲国产精品国产精品| 国产福利在线免费观看视频| 久久久久视频综合| a级毛片黄视频| 国产色婷婷99| 亚洲精品乱码久久久久久按摩| 男女边吃奶边做爰视频| 蜜桃国产av成人99| 亚洲国产日韩一区二区| 国产国语露脸激情在线看| 国产免费视频播放在线视频| av在线老鸭窝| 久久久亚洲精品成人影院| 有码 亚洲区| 美女国产高潮福利片在线看| 亚洲精华国产精华液的使用体验| 国产极品天堂在线| 久久青草综合色| 如何舔出高潮| 欧美日韩一区二区视频在线观看视频在线| 777米奇影视久久| 婷婷色综合www| 久久久a久久爽久久v久久| 天堂中文最新版在线下载| 亚洲成国产人片在线观看| 日韩免费高清中文字幕av| 91成人精品电影| 国产熟女欧美一区二区| 九色成人免费人妻av| 一二三四在线观看免费中文在 | 色吧在线观看| 熟女人妻精品中文字幕| 欧美老熟妇乱子伦牲交| 精品久久久久久电影网| 曰老女人黄片| 亚洲,欧美,日韩| 婷婷色综合大香蕉| 热re99久久国产66热| 亚洲少妇的诱惑av| 22中文网久久字幕| 制服人妻中文乱码| 美女视频免费永久观看网站| 香蕉丝袜av| 亚洲国产欧美在线一区| 国产黄色视频一区二区在线观看| 日本黄色日本黄色录像| 久久久久久伊人网av| 高清欧美精品videossex| 欧美97在线视频| 母亲3免费完整高清在线观看 | 9191精品国产免费久久| 天天影视国产精品| 男女边摸边吃奶| 亚洲国产精品一区三区| 999精品在线视频| 国产在线免费精品| 自拍欧美九色日韩亚洲蝌蚪91| 下体分泌物呈黄色| 久久久久精品人妻al黑| av又黄又爽大尺度在线免费看| 男人爽女人下面视频在线观看| 久久久精品免费免费高清| 成人手机av| 啦啦啦中文免费视频观看日本| 国产爽快片一区二区三区| 久久久久久久久久人人人人人人| 多毛熟女@视频| 在线观看www视频免费| 精品酒店卫生间| 我的女老师完整版在线观看| 99热全是精品| 成人亚洲欧美一区二区av| 欧美精品国产亚洲| 中文字幕亚洲精品专区| 成人毛片a级毛片在线播放| 国国产精品蜜臀av免费| 三上悠亚av全集在线观看| 成人亚洲欧美一区二区av| 久久久国产一区二区| 麻豆乱淫一区二区| 99视频精品全部免费 在线| 九色成人免费人妻av| 99热全是精品| 香蕉丝袜av| 国国产精品蜜臀av免费| 久久精品国产a三级三级三级| 天堂中文最新版在线下载| 男女国产视频网站| 十八禁高潮呻吟视频| 日本黄大片高清| 成人18禁高潮啪啪吃奶动态图| 如何舔出高潮| 伊人亚洲综合成人网| 久久综合国产亚洲精品| 国产免费又黄又爽又色| 99香蕉大伊视频| 亚洲综合精品二区| 成人午夜精彩视频在线观看| 国产国语露脸激情在线看| 日本黄大片高清| 日产精品乱码卡一卡2卡三| 又大又黄又爽视频免费| 一本—道久久a久久精品蜜桃钙片| 捣出白浆h1v1| 熟女电影av网| 制服诱惑二区| 久久ye,这里只有精品| 国产69精品久久久久777片| 日韩,欧美,国产一区二区三区| 国产精品一国产av| 黄色怎么调成土黄色| 免费看av在线观看网站| 国产不卡av网站在线观看| 18在线观看网站| 青春草亚洲视频在线观看| 99久久精品国产国产毛片| av黄色大香蕉| 成人亚洲精品一区在线观看| 超色免费av| 高清av免费在线| 亚洲欧美中文字幕日韩二区| 精品一品国产午夜福利视频| 久久久亚洲精品成人影院| 国产亚洲最大av| 9热在线视频观看99| 三上悠亚av全集在线观看| 久久av网站| 成人手机av| 黑人猛操日本美女一级片| 国产精品一区www在线观看| 成人午夜精彩视频在线观看| 亚洲精品视频女| 一级a做视频免费观看| 国产 精品1| 国产又爽黄色视频| 亚洲av日韩在线播放| 成人国语在线视频| videossex国产| 欧美成人午夜免费资源| 亚洲成av片中文字幕在线观看 | 国产深夜福利视频在线观看| 亚洲在久久综合| 大香蕉久久成人网| 国产男人的电影天堂91| 久久国内精品自在自线图片| 国产日韩欧美视频二区| 精品亚洲乱码少妇综合久久| 看十八女毛片水多多多| 一区二区av电影网| 9色porny在线观看| 欧美3d第一页| 中国美白少妇内射xxxbb| 大话2 男鬼变身卡| 久久人人爽人人爽人人片va| 一级毛片我不卡| 黄网站色视频无遮挡免费观看| 亚洲三级黄色毛片| 欧美xxⅹ黑人| 亚洲伊人久久精品综合| 午夜福利视频精品| 老熟女久久久| 欧美+日韩+精品| 国产永久视频网站| 午夜久久久在线观看| a级毛片黄视频| 亚洲av在线观看美女高潮| 成人亚洲精品一区在线观看| 街头女战士在线观看网站| 国产成人精品婷婷| 在线 av 中文字幕| 成年动漫av网址| 日产精品乱码卡一卡2卡三| 亚洲国产欧美在线一区| 欧美 亚洲 国产 日韩一| av在线app专区| 亚洲久久久国产精品| 母亲3免费完整高清在线观看 | 五月伊人婷婷丁香| 天堂8中文在线网| 在线观看www视频免费| 欧美+日韩+精品| 男人操女人黄网站| 亚洲情色 制服丝袜| 精品一区二区三区视频在线| 精品国产一区二区三区四区第35| 国产成人精品久久久久久| 五月玫瑰六月丁香| 春色校园在线视频观看| 国产xxxxx性猛交| 中文天堂在线官网| 七月丁香在线播放| 中文字幕av电影在线播放| 亚洲国产欧美日韩在线播放| 男女无遮挡免费网站观看| 91午夜精品亚洲一区二区三区| 亚洲成国产人片在线观看| 亚洲精品国产av成人精品| 国产深夜福利视频在线观看| av有码第一页| 69精品国产乱码久久久| 国产1区2区3区精品| 久久精品国产鲁丝片午夜精品| 99精国产麻豆久久婷婷| 一级毛片 在线播放| 成人无遮挡网站| 欧美成人精品欧美一级黄| 日韩,欧美,国产一区二区三区| 国产熟女午夜一区二区三区| 国产日韩一区二区三区精品不卡| 欧美日韩av久久| 宅男免费午夜| 欧美日韩成人在线一区二区| 综合色丁香网| 亚洲精品第二区| 国产探花极品一区二区| 99香蕉大伊视频| 亚洲国产最新在线播放| av线在线观看网站| 国产免费现黄频在线看| 我的女老师完整版在线观看| av国产久精品久网站免费入址| 亚洲av在线观看美女高潮| 大香蕉久久网| 七月丁香在线播放| 最近手机中文字幕大全| 男女午夜视频在线观看 | av播播在线观看一区| 男女边吃奶边做爰视频| 国产精品熟女久久久久浪| 免费少妇av软件| 人成视频在线观看免费观看| 男女边吃奶边做爰视频| 欧美xxxx性猛交bbbb| 秋霞伦理黄片| 亚洲欧美精品自产自拍| 免费久久久久久久精品成人欧美视频 | av免费在线看不卡| 熟女人妻精品中文字幕| 亚洲精品美女久久av网站| 久久久久精品久久久久真实原创| 精品少妇内射三级| 久久久a久久爽久久v久久| 色5月婷婷丁香| 看非洲黑人一级黄片| 精品久久久精品久久久| 免费看不卡的av| 国产精品一国产av| 精品少妇久久久久久888优播| 男女午夜视频在线观看 | 黄色怎么调成土黄色| 男人爽女人下面视频在线观看| 免费看不卡的av| 久久99精品国语久久久| av国产精品久久久久影院| 热re99久久精品国产66热6| 免费高清在线观看日韩| 日韩不卡一区二区三区视频在线| 欧美 亚洲 国产 日韩一| 赤兔流量卡办理| 成人黄色视频免费在线看| 亚洲精品日本国产第一区| kizo精华| 男人操女人黄网站| 欧美日韩国产mv在线观看视频| 男人添女人高潮全过程视频| 国产极品天堂在线| 午夜免费鲁丝| 国产乱人偷精品视频| 亚洲一区二区三区欧美精品| 日韩免费高清中文字幕av| 欧美人与善性xxx| av又黄又爽大尺度在线免费看| 最近2019中文字幕mv第一页| 亚洲成av片中文字幕在线观看 | 狠狠婷婷综合久久久久久88av| 一边亲一边摸免费视频| 亚洲四区av| 欧美人与善性xxx| 亚洲欧洲国产日韩| www日本在线高清视频| 亚洲,欧美,日韩| 亚洲欧美成人综合另类久久久| 人妻人人澡人人爽人人| 亚洲精品久久成人aⅴ小说| av福利片在线| 久久毛片免费看一区二区三区| 免费少妇av软件| 亚洲在久久综合| 国产色婷婷99| 精品一区二区三区四区五区乱码 | 最近手机中文字幕大全| 在线精品无人区一区二区三| 久久久国产一区二区| 一级黄片播放器| 中文字幕人妻丝袜制服| xxxhd国产人妻xxx| 午夜影院在线不卡| 国产一区二区在线观看日韩| 丝袜脚勾引网站| 国产av精品麻豆| 日本爱情动作片www.在线观看| 亚洲国产看品久久| 菩萨蛮人人尽说江南好唐韦庄| 国产无遮挡羞羞视频在线观看| 国产高清不卡午夜福利| 9色porny在线观看| h视频一区二区三区| 欧美bdsm另类| 水蜜桃什么品种好| 国产永久视频网站| 久久久久久久久久人人人人人人| 日韩大片免费观看网站| 韩国av在线不卡| 亚洲中文av在线| 国产国拍精品亚洲av在线观看| 性色av一级| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦中文免费视频观看日本| 老女人水多毛片| av播播在线观看一区| 日韩欧美精品免费久久| 国产av一区二区精品久久| av福利片在线| 97人妻天天添夜夜摸| 欧美激情极品国产一区二区三区 | 高清欧美精品videossex| av视频免费观看在线观看| 国产成人免费观看mmmm| 国产精品久久久久久久久免| 母亲3免费完整高清在线观看 | 日韩中文字幕视频在线看片| 激情视频va一区二区三区| 九九在线视频观看精品| 在线观看美女被高潮喷水网站| 少妇精品久久久久久久| 久久99蜜桃精品久久| 色婷婷久久久亚洲欧美| 色5月婷婷丁香| 久久这里只有精品19| 熟女人妻精品中文字幕| 欧美xxⅹ黑人| 欧美精品av麻豆av| 免费大片18禁| 老熟女久久久| 亚洲精品,欧美精品| 亚洲高清免费不卡视频| 伊人久久国产一区二区| av在线观看视频网站免费| 国产成人免费无遮挡视频| 巨乳人妻的诱惑在线观看| 欧美日韩精品成人综合77777| 人妻少妇偷人精品九色| 狠狠婷婷综合久久久久久88av| 国产亚洲最大av| 在线观看人妻少妇| 国产 一区精品| 最近的中文字幕免费完整| 精品熟女少妇av免费看| 久久精品久久精品一区二区三区| 蜜桃在线观看..| 搡女人真爽免费视频火全软件| 水蜜桃什么品种好| 全区人妻精品视频| 天美传媒精品一区二区| 国产淫语在线视频| 777米奇影视久久| 女人被躁到高潮嗷嗷叫费观| 一级,二级,三级黄色视频| 大片免费播放器 马上看| 熟女av电影| 免费av中文字幕在线| 成年女人在线观看亚洲视频| 国产一区二区在线观看日韩| 久久精品夜色国产| 国产免费福利视频在线观看| 黄片播放在线免费| 久久精品熟女亚洲av麻豆精品| freevideosex欧美| 免费日韩欧美在线观看| 久久久久久久久久人人人人人人| 天美传媒精品一区二区| 久久久欧美国产精品| 欧美精品国产亚洲| videos熟女内射| 欧美日韩精品成人综合77777| 国产精品久久久久久精品电影小说| 国产精品一区二区在线观看99| 亚洲 欧美一区二区三区| 久久亚洲国产成人精品v| 国产精品三级大全| 亚洲国产看品久久| 七月丁香在线播放| 国产精品蜜桃在线观看| 两性夫妻黄色片 | 夜夜爽夜夜爽视频| av免费在线看不卡| 久久久久久久国产电影| 国产黄色免费在线视频| 又黄又粗又硬又大视频| 亚洲成人av在线免费| 久久韩国三级中文字幕| 亚洲国产精品999| 高清av免费在线| 极品少妇高潮喷水抽搐| 国产国语露脸激情在线看| 国产精品蜜桃在线观看| 国产熟女午夜一区二区三区| 国产乱来视频区| 伊人亚洲综合成人网| 美女国产高潮福利片在线看| 人体艺术视频欧美日本| 91精品国产国语对白视频| 欧美人与善性xxx| av一本久久久久| 免费黄频网站在线观看国产| 色婷婷av一区二区三区视频| 婷婷色av中文字幕| a级毛片在线看网站| 欧美 亚洲 国产 日韩一| 美女脱内裤让男人舔精品视频| 国产成人免费无遮挡视频| 在线精品无人区一区二区三| 青春草亚洲视频在线观看| 80岁老熟妇乱子伦牲交| 国产成人精品婷婷| 久热这里只有精品99| 国产精品嫩草影院av在线观看| 国产午夜精品一二区理论片| 国产成人精品婷婷| 日本欧美视频一区| 国产一区亚洲一区在线观看| 在线观看人妻少妇| 国产精品三级大全| 国产在线免费精品| 国产片特级美女逼逼视频| 欧美精品一区二区免费开放| 91成人精品电影| 97在线视频观看| www.av在线官网国产| 亚洲四区av| 国产黄频视频在线观看| 男女午夜视频在线观看 | 在线天堂中文资源库| 一本大道久久a久久精品| 丝袜人妻中文字幕| 国产深夜福利视频在线观看| √禁漫天堂资源中文www| 亚洲精品美女久久av网站| 我要看黄色一级片免费的| 亚洲精品久久久久久婷婷小说| 九九爱精品视频在线观看| 一边亲一边摸免费视频| 免费大片黄手机在线观看| 久久国内精品自在自线图片| 黄片无遮挡物在线观看| 乱码一卡2卡4卡精品| av天堂久久9| 欧美最新免费一区二区三区| 免费黄色在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻在线不人妻| 日日摸夜夜添夜夜爱| 永久网站在线| 美女视频免费永久观看网站| 国产男人的电影天堂91| 久久久久精品人妻al黑| 五月天丁香电影| 免费观看a级毛片全部| 天天影视国产精品| 免费观看性生交大片5| 国产精品久久久久久久电影| 另类精品久久| 精品午夜福利在线看| 91精品三级在线观看| 99久国产av精品国产电影| 香蕉国产在线看| 在线观看一区二区三区激情| 看免费av毛片| 欧美精品一区二区大全| 成年av动漫网址| 免费av中文字幕在线| 久久久国产精品麻豆| 一二三四在线观看免费中文在 | 欧美亚洲日本最大视频资源| 国产男人的电影天堂91| 国产欧美亚洲国产| 9热在线视频观看99| 久久久久精品久久久久真实原创| 80岁老熟妇乱子伦牲交| 精品国产国语对白av| 午夜日本视频在线| 乱码一卡2卡4卡精品| 国产精品女同一区二区软件| 欧美xxⅹ黑人| 中文字幕亚洲精品专区| 国国产精品蜜臀av免费| 久久国产精品大桥未久av| 少妇熟女欧美另类| 亚洲精品乱久久久久久| 国产一区二区三区综合在线观看 | 欧美日韩成人在线一区二区| 一个人免费看片子| 看免费av毛片| 母亲3免费完整高清在线观看 | 人妻 亚洲 视频| 午夜福利视频精品| 最新的欧美精品一区二区| 新久久久久国产一级毛片| 亚洲情色 制服丝袜| 亚洲欧美精品自产自拍| 亚洲国产成人一精品久久久| 免费av中文字幕在线| 免费高清在线观看日韩| 性色av一级| 啦啦啦在线观看免费高清www| 亚洲精品美女久久av网站| 人妻少妇偷人精品九色| 亚洲伊人久久精品综合| 高清视频免费观看一区二区| 国产老妇伦熟女老妇高清| 国产欧美日韩综合在线一区二区| 秋霞伦理黄片| 丝袜美足系列| av在线观看视频网站免费| 天天影视国产精品| 卡戴珊不雅视频在线播放| 婷婷色麻豆天堂久久| 免费黄频网站在线观看国产| 久久精品久久精品一区二区三区| 精品人妻偷拍中文字幕| 国产在线一区二区三区精| 两个人看的免费小视频| 高清黄色对白视频在线免费看| 成年av动漫网址| 青春草亚洲视频在线观看| 亚洲av日韩在线播放| 免费黄网站久久成人精品| 久久影院123| 男女免费视频国产| 亚洲图色成人| 黑人欧美特级aaaaaa片| 天天操日日干夜夜撸| 高清在线视频一区二区三区| 视频中文字幕在线观看| 精品人妻偷拍中文字幕| 国产亚洲av片在线观看秒播厂| 久久精品久久精品一区二区三区| 91午夜精品亚洲一区二区三区| 久久久国产欧美日韩av| 亚洲,一卡二卡三卡| 2022亚洲国产成人精品| 精品一区二区三卡| 王馨瑶露胸无遮挡在线观看| 国产精品一国产av| 考比视频在线观看| av卡一久久| 如日韩欧美国产精品一区二区三区| 日韩一区二区视频免费看| 亚洲精品国产av成人精品| 又大又黄又爽视频免费|