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

    CALMET模式不同參數(shù)化方案對江西省山地 風場模擬的對比分析

    2020-04-18 02:44:30姚琳徐衛(wèi)民彭王敏子溫新龍沈競
    關鍵詞:高度層風塔測風塔

    姚琳 徐衛(wèi)民 彭王敏子 溫新龍 沈競

    (江西省氣象科學研究所,南昌 330096)

    0 引言

    風能作為一種可再生的清潔能源,已成為許多國家推進能源轉型的核心發(fā)展內容和應對氣候變化的重要途徑。江西省風能資源主要集中在環(huán)鄱陽湖地區(qū)與高海拔山地,相對規(guī)?;虚_發(fā)大型風電場而言,江西省更適合開發(fā)分散式風電,因此江西省山地風場的精細化準確模擬對提高風電選址和風資源開發(fā)有著重要意義。

    目前對風能資源風速預測主要有統(tǒng)計法和數(shù)值模擬法,統(tǒng)計法由于需要的數(shù)據(jù)量大等缺點沒有得到廣泛應用[1-3],計算機技術的發(fā)展使數(shù)值模擬法廣泛應用于風能資源評估方面[4-8],隨著數(shù)值模式水平的不斷提高,分辨率更高更為精細的小尺度模式在特定的區(qū)域風場模擬方面得到了很好的應用,其中CALMET模式是美國EPA推薦的由美國地球技術公司(Earth Tech Inc)開發(fā)的可用于復雜地形的邊界層微尺度氣象診斷模式,基于CALMET模式的動力降尺度方法是國內外風資源評估的手段之一[9-13]。CALMET模式的風場選項參數(shù)包含了地形動力學效應、下坡氣流效應、地形阻塞效應和O’Brien垂直速度調整等,由于數(shù)值模式在不同氣候背景、地形地貌區(qū)域對風場模擬的效果與參數(shù)化方案的適應性緊密相關,所以選取適合區(qū)域特點的參數(shù)化方案是數(shù)值模式模擬的關鍵[14-17]。

    本文利用CALMET模式對江西省境內高海拔山地風電場測風塔進行1月和7月的數(shù)值模擬,對比分析不同參數(shù)化方案對山地風場逐時風速的模擬效果,從而選出最優(yōu)參數(shù)化方案,并且對最優(yōu)參數(shù)化方案進一步進行整1年的模擬效果研究。

    1 模式及方案設計介紹

    1.1 模式介紹及模擬方案設置

    CALMET模式是中尺度與小尺度結合的模式系統(tǒng),利用中尺度數(shù)值模式WRF模式最內層格點輸出結果(3 km×3 km)作為CALMET的背景驅動場(初始猜測場),通過進一步的動力降尺度得到分辨率為0.5 km×0.5 km的診斷風場。

    文中WRF模式是由美國NOAA、NCEP、Air Force等多個研究部門聯(lián)合開發(fā)的中尺度數(shù)值預報模式,是一個全可壓非靜力模式。本文選用WRF3.6版本,模擬區(qū)域的中心點位于(25.6°N,115.2°E),27 km×9 km×3 km的3層嵌套網(wǎng)格,網(wǎng)格數(shù)分別為:52×55,82×91,91×121,三重嵌套網(wǎng)格垂直方向均分為不等距30層,模式頂約為19 km,至少15層在1 km以下。地形數(shù)據(jù)分別為美國地質勘探局(USGS)的全球2′、1′和30″地形數(shù)據(jù),模式初始場為全球1°×1°、6 h一次的FNL/NCEP再分析數(shù)據(jù)。通過多種參數(shù)化方案模擬試驗分析,選取模擬結果與實測分析接近的方案,最后確定為:WSM3微物理過程方案+MRF邊界層方案+Noah陸面過程方案[18]。WRF模式對研究區(qū)域進行為期1年的逐日模擬,設置為每天運行36 h,起止時間為每日12時(世界時)至第三日00時,模擬結果取每天模擬的后24 h逐時輸出結果。

    CALMET模式是美國EPA推薦的由美國地球技術公司(Earth Tech Inc)開發(fā)的邊界層微尺度氣象診斷模式,CALMET模式包括診斷風場模塊和微氣象模塊,它通過質量守恒連續(xù)方程對風場進行診斷,在輸入模式所需的常規(guī)氣象觀測資料和中尺度氣象模式輸出場后,進行地形動力學、下坡氣流、地形阻塞效應調整,并通過插值、平滑處理、垂直速度計算、輻散最小化等產(chǎn)生最終風場。本文計算采用的是CALMET v6.334的版本,CALMET模式水平網(wǎng)格數(shù)為200×200,網(wǎng)格大小為0.5 km×0.5 km,垂直分為11層,地形資料采用水平分辨率3″的SRTM3資料,下墊面類型資料為30″水平分辨率的USGS資料。

    為了研究適用于江西山地風電場的CALMET模式最優(yōu)參數(shù)化方案,選取了位于江西省境內高海拔山地風電場的2座測風塔資料對多種參數(shù)化方案的模擬結果進行效果評估。2座測風塔為NRG測風系統(tǒng),均位于興國縣北部丘陵山區(qū)大水山一帶,周圍均無遮障并暴露于強烈高空風中,海拔高度分別為879 m和831 m,均為5層測風塔,每層各有一套風速儀,在80 m、10 m高度層安裝了風向儀;塔層高度分別為10 m、30 m、50 m、70 m和80 m,記錄為10 min風速平均,觀測時間為2016年6月—2017年5月,測風塔實測數(shù)據(jù)使用時進行了嚴格的質量控制,按照國家標準《風電場風能資源評估方法》(GB/T 18710—2002)對測風塔原始數(shù)據(jù)進行完整性、合理性及有效性分析檢驗,有效完整均達80%以上。模式區(qū)域及測風塔位置如圖1。

    圖1 模式區(qū)域(紅色實線為WRF,黑色實線為CALMET)及測風塔位置 Fig. 1 The simulation experiment domain and the location of the wind tower

    1.2 參數(shù)化方案組設計

    通過CALMET模式準確地診斷風場,精細化的地形資料和下墊面類型、高質量和高密度的觀測資料均至關重要,但并不是任何區(qū)域能具備這樣的條件,特別是對于山地風場,氣象站分布往往距測風塔較遠、稀疏且地形差距較大,且測風高度也僅有10 m,因此,此次CALMET模式地面及高空氣象輸入資料均由中尺度氣象模式WRF的模擬結果提供,并且主要針對山地風場50~80 m高度層進行模擬研究。

    CALMET模式的風場選項參數(shù)有地形動力學效應、下坡氣流效應、地形阻塞效應和O’Brien垂直速度調整等。地形動力學效應是指地形對氣團速度的強迫作用,主要通過計算整個區(qū)域的風來獲得受地形影響的垂直風速,并滿足大氣穩(wěn)定度遞減指數(shù)函數(shù),對初始猜測風場重復執(zhí)行輻散最小化方法,直到三維輻散小于閾值,以獲得水平方向風分量所受到的地形動力學影響。下坡氣流效應是利用地形坡面、坡高、時間等參數(shù)計算,其風分量調入風場調整空氣動力學影響。地形阻塞效應是指地形對風場的熱力學阻塞效應通過局地弗勞德(Froude)數(shù)計算,如果網(wǎng)格點計算值小于臨界Froude數(shù),且風有上坡風量,則風向調整為與地形相切的方向,風速不變;如果超過臨界Froude數(shù),則不需要調整。O’Brien垂直速度調整是根據(jù)Godden和Lurmann的理論假定模擬區(qū)域的層頂垂直風速為0 m·s-1[19]。此次對CALMET模式的以上4個風場選項設計了5組參數(shù)化方案,模擬了2016年7月和2017年1月江西某山地風電場2座測風塔逐時風速風向,對50 m、70 m和80 m高度層模擬結果與測風塔數(shù)據(jù)進行比對分析,方案設計見表1。

    表1 CALMET模式參數(shù)化方案 Table 1 Comparison of parameterization schemes in the CALMET model

    2 CALMET模式參數(shù)化方案模擬效果對比

    2.1 全風速模擬效果對比

    平均絕對誤差能反映出模擬值精度的實際情況,均方根誤差反映了模擬結果與觀測結果的離散程度,相關系數(shù)則可說明模擬值與觀測值的變化趨勢逼近情況。表2和表3分別給出了CALMET模式各組參數(shù)化方案對2座測風塔逐小時風速模擬結果與實測風速相關系數(shù)(R)、平均絕對誤差(MAE)和均方根誤差(RMSE)[20]。由表2結果可知,7月相關性比1月更好,數(shù)值上各組結果差距不大。圖2為測風塔逐小時風速模擬檢驗結果,從圖中可以看出,1月模擬效果優(yōu)于天氣過程較多的7月,特別是表現(xiàn)在2#塔模擬效果中,1月各高度層的均方根誤差和平均絕對誤差均小于7月。各組方案檢驗結果上看,第3組模擬效果最差,說明CALMET模式采用地形動力學時,地形強迫產(chǎn)生的垂直速度對水平風場調整不夠充分,會導致模擬結果誤差較大,特別是天氣過程較多的7月份。為了防止模式格點的高層產(chǎn)生異常的大垂直風速,通過采用O’Brien對垂直風速進行調整強迫模式模擬區(qū)域頂層的垂直風速為0 m·s-1,進而讓垂直風速更合理的對水平風速進行調整,通過比對方案檢驗結果第5組在均方根誤差和平均絕對誤差上均為最小,為最優(yōu)方案,即動力降尺度CALMET模式對山地風場風速模擬時,不宜采用地形動力效應參數(shù),采用下坡氣流效應調整對50 m以上高度層風場影響較小,關閉Froude數(shù)調整以及采用O’Brien垂直風速調整更利于提高對山地風場的模擬效果。

    表2 測風塔逐小時風速模擬相關性結果 Table 2 Correlation coefficients of simulated hourly wind speed from the wind tower

    2.2 風向模擬效果對比

    利用CALMET模式各組參數(shù)化方案模擬2017年1月和2016年7月山地風場1#測風塔逐時風向,通過風向統(tǒng)計結果得到各月風向玫瑰對比圖(圖3~圖4),由圖可知,1月各組模擬結果主導風向均與實測主導風向一致,第1、2、4和5組主導風向頻率較實測值為3%~13%,第3組方案組主導風向頻率較實測值偏小約5%,7月各組模擬結果主導風向均與實測主導風向有一個方位的偏差,主導風頻偏大16%~24%,總體上各組方案均能較好地模擬出各測風塔的主導風向。

    3 最優(yōu)參數(shù)化方案效果檢驗

    利用上述所選出的CALMET模式最優(yōu)方案第5組對山地風場兩座測風塔進行2016年6月至2017年5月整一年的逐小時風資源模擬,對模擬效果進行進一步分析。

    3.1 指標檢驗

    表3給出了CALMET模式最優(yōu)方案對山地2座測風塔連續(xù)一年的逐小時風速模擬結果,1#、2#測風塔全年模擬結果相關系數(shù)分別介于0.62~0.79和0.53~0.81,均方根誤差分別介于2.0~2.7 m·s-1和2.1~2.7 m·s-1,絕對平均偏差均介于1.6~2.0 m·s-1。CALMET模式對山地風場整一年的各層風速模擬年均相關系數(shù)為0.69,年均均方根誤差和絕對平均偏差分別為2.3 m·s-1和1.8 m·s-1。

    表3 山地風場測風塔逐小時風速模擬效果檢驗結果 Table 3 Simulated results of hourly wind speed of wind tower in the mountain wind field

    圖2 測風塔逐小時風速模擬效果檢驗結果 (a)1#塔均方根誤差;(b)2#塔均方根誤差;(c)1#塔平均絕對誤差;(d)2#塔平均絕對誤差 Fig. 2 Test results of hourly wind speed simulation of wind tower (a) Root mean square error of 1# wind tower; (b) Root mean square error of 2# wind tower; (c) Mean absolute error of 1# wind tower; (d) Mean absolute error of 2# wind tower

    從總體平均結果來看,CALMET模式能較好的模擬出山地風場測風塔逐時風速,四季模擬效果差別較小,夏季相關系數(shù)優(yōu)于其他季節(jié),秋季誤差較其他季節(jié)最小。

    3.2 風速模擬效果檢驗

    1#塔實測80 m、70 m和50 m年均風速分別為6.4 m·s-1、6.2 m·s-1和6.0 m·s-1,模擬值年均風速分別為6.3 m·s-1、6.2 m·s-1和5.9 m·s-1,2#塔實測80 m、70 m和50 m年均風速分別為5.6 m·s-1、5.6 m·s-1和5.5 m·s-1,模擬值年均風速分別為5.8 m·s-1、5.9 m·s-1和5.5 m·s-1。

    圖5給出了2座測風塔80 m高度層實測與模擬風速月均變化,由圖可知,CALMET模式對測風塔的模擬結果與實測值風速月均變化趨勢基本一致,特別是1#測風塔風速較大月份4月、6月、7月和較小月份5月、8月,以及2#測風塔風速較大月份6月、7月和較小月份5月、8月均有很好的模擬效果。

    圖3 1#測風塔80 m高度層2017年1月風向玫瑰圖 Fig. 3 Wind rose diagram of 1# wind tower, 80 m high, January 2017

    圖4 1#測風塔80 m高度層2016年7月風向玫瑰圖 Fig. 4 Wind rose diagram of 1# wind tower, 80 m high, July 2016

    圖5 1#(a)、2#(b)測風塔80 m高度層風速月均變化 Fig. 5 Monthly variation of wind speed at 80 m height layer of 1#(a) and 2#(b) wind towers

    圖6 1#(a)、2#(b)測風塔全年80 m高度層風速分布 Fig. 6 Distribution of wind speed section of 80 m in the 1#(a) and 2#(b) wind towers

    圖6給出了2座測風塔80 m高度層實測值和模擬值風速段分布,由圖可知,CALMET模式對風速段分布的模擬效果與實測基本一致,1#測風塔實測值和模擬值風速分布均主要集中在3~9 m·s-1,但在風速峰值區(qū)模擬值概率偏大且與實測值差別較大,有2%~4%的偏差。2#測風塔實測值和模擬值風速分布均主要集中在2~8 m·s-1,在風速峰值區(qū)實測值偏大且與模擬值概率差別較大,有2%~5%的偏差。

    3.3 風向檢驗

    圖7給出了各測風塔模式模擬及實測風向玫瑰圖,由圖可知,模擬與實測結果主導風向及風能分布有一定的差異,約有一個方位的偏差,1#塔實測值主導風向為N方向,模擬結果主導風向為NNE方向,頻率上約有7%的偏差;2#塔實測值主導風向為SSW方向,模擬結果主導風向為S方向,頻率上約有8%的偏差。

    圖7 1#(a)、2#(b)測風塔80 m高度層全年風向玫瑰圖 Fig.7 Annual wind rose of 1# (a) and 2# (b) wind towers at 80 m

    4 結論

    基于不同參數(shù)化方案的CALMET模式對江西省境內山地風場的模擬,優(yōu)選出最佳參數(shù)化方案并進一步進行全年山地風場模擬,得出以下結論。

    1)采用動力降尺度CALMET模式對江西省山地風場2016年7月和2017年1月風速進行模擬時,通過5組參數(shù)化方案模擬結果的比對分析得到:第5組參數(shù)化方案最優(yōu),即CALMET模式以不采用地形動力效應參數(shù)調整和Froude數(shù)調整,采用下坡氣流效應調整和O’Brien垂直風速調整時,對江西省境內山地風場50 m以上高度層風場模擬效果最佳。

    2)采用CALMET模式最優(yōu)參數(shù)化方案對江西省山地風場進行整1年的風速進行模擬,模擬結果表明:CALMET模式能較好的模擬出山地風場測風塔逐時風速,全年四個季節(jié)模擬指標效果差別較小,夏季相關系數(shù)優(yōu)于其他季節(jié),秋季誤差較其他季節(jié)最小。其不僅能較好地模擬出風速及風能密度月均變化,對全年風速段分布模擬也與實測值較為一致,在峰值區(qū)分布概率有2%~5%的偏差。風向分布方面,CALMET模式能較好地模擬出實際測風塔全年主導風向,但模擬與實測結果主導風向分布約有一個方位的偏差,主導風向頻率有7%~8%的偏差。

    猜你喜歡
    高度層風塔測風塔
    一種自安裝海上測風塔的運輸和安裝穩(wěn)性分析
    WRF氣象模式在高海拔復雜山地風電場中的應用研究
    基于高度層的航路短時利用率模型研究
    測風塔法在棄風電量評估中的應用
    臺風“威馬遜”影響廣東期間近地層風特性
    廣西風能資源分布高分辨率數(shù)值模擬應用研究
    廣東電力(2016年6期)2016-07-16 04:28:40
    談風塔的施工技術要點
    河南省風能數(shù)據(jù)庫系統(tǒng)建設
    飛機最佳航路爬升時機研究
    國務院、中央軍委關于修改《中華人民共和國飛行基本規(guī)則》的決定
    亚洲精品在线美女| 午夜亚洲福利在线播放| 久久九九热精品免费| 色av中文字幕| 老司机深夜福利视频在线观看| 嫩草影视91久久| 精品国产美女av久久久久小说| 偷拍熟女少妇极品色| 给我免费播放毛片高清在线观看| 最近最新中文字幕大全电影3| 激情在线观看视频在线高清| 精品国产三级普通话版| 国产黄色小视频在线观看| 亚洲欧美激情综合另类| 9191精品国产免费久久| 12—13女人毛片做爰片一| 小说图片视频综合网站| 精品一区二区三区四区五区乱码| 俺也久久电影网| 免费观看人在逋| 亚洲aⅴ乱码一区二区在线播放| 欧美成人一区二区免费高清观看 | 色视频www国产| 成人三级做爰电影| www.熟女人妻精品国产| 久久国产精品人妻蜜桃| 亚洲国产欧美网| 欧美日韩瑟瑟在线播放| 日本成人三级电影网站| 1000部很黄的大片| 免费看a级黄色片| 亚洲片人在线观看| 日韩欧美精品v在线| 无遮挡黄片免费观看| 亚洲国产欧洲综合997久久,| 在线观看日韩欧美| 精品久久久久久久人妻蜜臀av| 日韩欧美 国产精品| 亚洲av免费在线观看| 午夜免费激情av| 国产av一区在线观看免费| 国产乱人伦免费视频| 午夜福利欧美成人| 色精品久久人妻99蜜桃| 精品久久久久久久久久免费视频| 12—13女人毛片做爰片一| 长腿黑丝高跟| 精品无人区乱码1区二区| 小蜜桃在线观看免费完整版高清| 女人高潮潮喷娇喘18禁视频| 黄频高清免费视频| 熟女少妇亚洲综合色aaa.| 岛国在线观看网站| 黄色女人牲交| 日本一二三区视频观看| 午夜精品一区二区三区免费看| 成年女人看的毛片在线观看| 国产午夜精品久久久久久| 国产真实乱freesex| 啦啦啦韩国在线观看视频| 亚洲成人久久性| 亚洲精品美女久久久久99蜜臀| 禁无遮挡网站| 午夜精品一区二区三区免费看| 国产黄片美女视频| 婷婷精品国产亚洲av在线| 五月玫瑰六月丁香| 亚洲七黄色美女视频| 1024手机看黄色片| 18禁黄网站禁片午夜丰满| 午夜免费激情av| 国产97色在线日韩免费| 黑人欧美特级aaaaaa片| 国产成+人综合+亚洲专区| 全区人妻精品视频| 国产精品1区2区在线观看.| 久久香蕉精品热| 成人高潮视频无遮挡免费网站| 三级男女做爰猛烈吃奶摸视频| 日韩精品青青久久久久久| 久久99热这里只有精品18| 90打野战视频偷拍视频| 在线观看午夜福利视频| 成年女人毛片免费观看观看9| 成人欧美大片| 999精品在线视频| 免费电影在线观看免费观看| 黄色视频,在线免费观看| 青草久久国产| 成人18禁在线播放| 国产极品精品免费视频能看的| 亚洲18禁久久av| 精品欧美国产一区二区三| 国产av麻豆久久久久久久| 日韩人妻高清精品专区| 免费观看人在逋| 色综合欧美亚洲国产小说| 亚洲av五月六月丁香网| 亚洲熟妇中文字幕五十中出| 色视频www国产| 国内精品一区二区在线观看| 亚洲av美国av| 亚洲欧美激情综合另类| 变态另类成人亚洲欧美熟女| 国产精品一区二区三区四区久久| e午夜精品久久久久久久| 老司机深夜福利视频在线观看| 午夜影院日韩av| 五月玫瑰六月丁香| 亚洲自偷自拍图片 自拍| 国产又黄又爽又无遮挡在线| 亚洲男人的天堂狠狠| 国产人伦9x9x在线观看| 国产亚洲精品av在线| 在线a可以看的网站| 成年女人看的毛片在线观看| 丁香欧美五月| 亚洲熟妇熟女久久| 国产av不卡久久| 亚洲av电影在线进入| 麻豆国产97在线/欧美| 在线国产一区二区在线| 国产亚洲精品综合一区在线观看| 老熟妇乱子伦视频在线观看| 神马国产精品三级电影在线观看| 国产精品免费一区二区三区在线| 亚洲熟妇中文字幕五十中出| av在线天堂中文字幕| 亚洲人成伊人成综合网2020| 看片在线看免费视频| 熟女电影av网| 免费看a级黄色片| 五月玫瑰六月丁香| 俺也久久电影网| 真实男女啪啪啪动态图| 一区福利在线观看| 久久这里只有精品中国| 一个人免费在线观看电影 | 久久伊人香网站| 又黄又粗又硬又大视频| 国产精品1区2区在线观看.| 18禁观看日本| 男女床上黄色一级片免费看| 成年人黄色毛片网站| 国产成人啪精品午夜网站| 免费在线观看视频国产中文字幕亚洲| 国产精品影院久久| 久久久久久大精品| 欧美色欧美亚洲另类二区| 欧美一级a爱片免费观看看| 免费看光身美女| 18禁观看日本| 天堂动漫精品| 这个男人来自地球电影免费观看| 国产亚洲欧美98| 小蜜桃在线观看免费完整版高清| av天堂在线播放| 一个人免费在线观看的高清视频| 国产美女午夜福利| 亚洲在线自拍视频| 两个人的视频大全免费| 成年女人永久免费观看视频| aaaaa片日本免费| 在线国产一区二区在线| www.自偷自拍.com| 亚洲精品粉嫩美女一区| 给我免费播放毛片高清在线观看| 亚洲国产中文字幕在线视频| 亚洲专区字幕在线| 亚洲 欧美 日韩 在线 免费| 国产欧美日韩一区二区三| 国产成人影院久久av| 老司机午夜福利在线观看视频| 成人欧美大片| 色哟哟哟哟哟哟| 国产午夜精品久久久久久| 欧美中文日本在线观看视频| 国产乱人伦免费视频| 国产精品免费一区二区三区在线| 老鸭窝网址在线观看| 国产69精品久久久久777片 | 曰老女人黄片| 亚洲成人中文字幕在线播放| 午夜福利欧美成人| 狠狠狠狠99中文字幕| 三级毛片av免费| 日日摸夜夜添夜夜添小说| 国产精品亚洲一级av第二区| 精品电影一区二区在线| 精品久久久久久久毛片微露脸| 午夜免费激情av| 精品一区二区三区四区五区乱码| 九色国产91popny在线| 精品一区二区三区视频在线观看免费| 精品国产亚洲在线| av中文乱码字幕在线| 亚洲精品中文字幕一二三四区| 欧美日韩福利视频一区二区| 亚洲在线观看片| 国产精品国产高清国产av| 国产91精品成人一区二区三区| 免费看十八禁软件| 熟女少妇亚洲综合色aaa.| 看黄色毛片网站| 最好的美女福利视频网| 久久中文字幕一级| 成人性生交大片免费视频hd| 欧美一级a爱片免费观看看| 国产精品亚洲一级av第二区| 久久中文看片网| cao死你这个sao货| 日本撒尿小便嘘嘘汇集6| 成人一区二区视频在线观看| 亚洲欧美日韩无卡精品| 淫秽高清视频在线观看| 老司机午夜十八禁免费视频| 久久午夜综合久久蜜桃| 午夜视频精品福利| av欧美777| а√天堂www在线а√下载| 黄色片一级片一级黄色片| 日韩欧美在线乱码| 国产一级毛片七仙女欲春2| 亚洲国产欧洲综合997久久,| 日韩欧美免费精品| 午夜福利在线在线| 免费一级毛片在线播放高清视频| 欧美三级亚洲精品| 久久精品影院6| 成人性生交大片免费视频hd| 亚洲人成伊人成综合网2020| 啪啪无遮挡十八禁网站| 丰满人妻一区二区三区视频av | 亚洲成av人片在线播放无| 大型黄色视频在线免费观看| 哪里可以看免费的av片| 毛片女人毛片| 日本三级黄在线观看| 一区二区三区高清视频在线| 男女下面进入的视频免费午夜| 日韩高清综合在线| 亚洲 欧美一区二区三区| 一区二区三区高清视频在线| 亚洲国产欧洲综合997久久,| 中文字幕精品亚洲无线码一区| 91在线观看av| 久久国产乱子伦精品免费另类| 国模一区二区三区四区视频 | 日本黄色视频三级网站网址| 99国产极品粉嫩在线观看| 国产精品av久久久久免费| 欧美中文日本在线观看视频| 免费一级毛片在线播放高清视频| 日本五十路高清| 一夜夜www| 最近最新中文字幕大全免费视频| 日本撒尿小便嘘嘘汇集6| 免费在线观看视频国产中文字幕亚洲| 亚洲精品在线美女| 午夜久久久久精精品| 国产视频一区二区在线看| 久久久水蜜桃国产精品网| 久久久久性生活片| 日韩国内少妇激情av| 国产不卡一卡二| 欧美成狂野欧美在线观看| 国产精品 欧美亚洲| 婷婷亚洲欧美| 久久精品国产清高在天天线| 18禁裸乳无遮挡免费网站照片| 久久国产精品人妻蜜桃| 老熟妇仑乱视频hdxx| or卡值多少钱| 18禁黄网站禁片午夜丰满| 久久久久免费精品人妻一区二区| 亚洲美女黄片视频| 狠狠狠狠99中文字幕| 国产精品久久电影中文字幕| 国产视频一区二区在线看| 国产av不卡久久| 午夜激情欧美在线| 精品久久久久久久人妻蜜臀av| 亚洲乱码一区二区免费版| 日日干狠狠操夜夜爽| 欧美国产日韩亚洲一区| 亚洲成人免费电影在线观看| 美女黄网站色视频| 成人av在线播放网站| 中文资源天堂在线| 在线观看一区二区三区| 久久亚洲真实| 国产 一区 欧美 日韩| 免费在线观看影片大全网站| 欧美日韩中文字幕国产精品一区二区三区| 国产成人一区二区三区免费视频网站| 国产亚洲精品av在线| 一a级毛片在线观看| 毛片女人毛片| 日韩欧美在线乱码| 欧美在线黄色| 免费观看的影片在线观看| 久久久久久久久中文| 国产成人系列免费观看| 亚洲精品久久国产高清桃花| 一级毛片精品| 美女大奶头视频| 18禁裸乳无遮挡免费网站照片| 看片在线看免费视频| 午夜福利欧美成人| 久久久国产成人免费| 美女被艹到高潮喷水动态| 2021天堂中文幕一二区在线观| 日韩av在线大香蕉| 性欧美人与动物交配| 黄色成人免费大全| 波多野结衣高清作品| 国产美女午夜福利| 免费在线观看视频国产中文字幕亚洲| 国产高清激情床上av| 男插女下体视频免费在线播放| 国产探花在线观看一区二区| 无遮挡黄片免费观看| 麻豆av在线久日| 国产激情偷乱视频一区二区| 搡老岳熟女国产| 中亚洲国语对白在线视频| 中文字幕久久专区| 久久久国产欧美日韩av| 中国美女看黄片| 丁香六月欧美| 亚洲片人在线观看| 日本三级黄在线观看| 国产不卡一卡二| 欧美成人免费av一区二区三区| 无遮挡黄片免费观看| 国产视频内射| 村上凉子中文字幕在线| 亚洲精品456在线播放app | 一本一本综合久久| 老汉色∧v一级毛片| 久久九九热精品免费| 久久香蕉精品热| 在线观看午夜福利视频| 全区人妻精品视频| 成年人黄色毛片网站| 日本黄色视频三级网站网址| 好男人在线观看高清免费视频| 久久精品91蜜桃| 嫁个100分男人电影在线观看| 国产精品98久久久久久宅男小说| 在线观看午夜福利视频| 夜夜看夜夜爽夜夜摸| 99riav亚洲国产免费| 一二三四在线观看免费中文在| 99久久无色码亚洲精品果冻| 国产精品1区2区在线观看.| 男女视频在线观看网站免费| 麻豆成人av在线观看| 欧美成人一区二区免费高清观看 | 日韩有码中文字幕| 97碰自拍视频| 国产亚洲欧美98| 午夜影院日韩av| 欧美丝袜亚洲另类 | 日本三级黄在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av电影不卡..在线观看| 草草在线视频免费看| svipshipincom国产片| 久久99热这里只有精品18| 亚洲精品乱码久久久v下载方式 | 久久性视频一级片| 我的老师免费观看完整版| 毛片女人毛片| 久久久久久九九精品二区国产| 午夜精品久久久久久毛片777| 麻豆久久精品国产亚洲av| 欧美乱妇无乱码| 999久久久精品免费观看国产| 成年人黄色毛片网站| av在线蜜桃| 精品熟女少妇八av免费久了| 国产精品久久久久久亚洲av鲁大| 久久久国产成人精品二区| 亚洲欧美日韩无卡精品| 亚洲国产中文字幕在线视频| 18禁国产床啪视频网站| 中文字幕久久专区| a级毛片a级免费在线| 色综合婷婷激情| 亚洲性夜色夜夜综合| 欧美日本视频| 日日干狠狠操夜夜爽| 白带黄色成豆腐渣| 十八禁人妻一区二区| 亚洲 欧美 日韩 在线 免费| 在线观看免费午夜福利视频| 免费看a级黄色片| 狂野欧美激情性xxxx| 午夜福利欧美成人| 欧美大码av| 97超视频在线观看视频| 啦啦啦韩国在线观看视频| 久久久久久国产a免费观看| 亚洲国产欧美网| 91av网站免费观看| 国产黄片美女视频| 成年人黄色毛片网站| 午夜视频精品福利| 国产亚洲欧美98| 一级a爱片免费观看的视频| 午夜激情欧美在线| 亚洲中文字幕日韩| 日韩av在线大香蕉| 国产精品乱码一区二三区的特点| 最近最新免费中文字幕在线| 99视频精品全部免费 在线 | 亚洲国产精品sss在线观看| 亚洲av电影不卡..在线观看| 国产亚洲欧美98| 亚洲成人中文字幕在线播放| 黄色丝袜av网址大全| 黄片小视频在线播放| 高潮久久久久久久久久久不卡| 19禁男女啪啪无遮挡网站| 亚洲av成人一区二区三| 国产高清有码在线观看视频| 九九在线视频观看精品| 亚洲国产色片| 亚洲乱码一区二区免费版| 午夜福利成人在线免费观看| 国产91精品成人一区二区三区| 夜夜躁狠狠躁天天躁| 成年女人永久免费观看视频| 又黄又粗又硬又大视频| 国模一区二区三区四区视频 | 美女高潮的动态| 男人舔女人下体高潮全视频| 久久国产精品人妻蜜桃| 90打野战视频偷拍视频| 国产精品国产高清国产av| 又大又爽又粗| 国产爱豆传媒在线观看| 狂野欧美激情性xxxx| 99久久精品热视频| 国产精品电影一区二区三区| 一个人免费在线观看电影 | 在线观看免费午夜福利视频| 日韩成人在线观看一区二区三区| 国产精品永久免费网站| 国产一区二区在线av高清观看| 9191精品国产免费久久| 男女视频在线观看网站免费| 我的老师免费观看完整版| 免费在线观看成人毛片| 亚洲片人在线观看| 精品国内亚洲2022精品成人| 国产99白浆流出| 久久久久九九精品影院| 此物有八面人人有两片| 精品久久久久久成人av| 一二三四社区在线视频社区8| 波多野结衣巨乳人妻| 少妇裸体淫交视频免费看高清| 成人高潮视频无遮挡免费网站| 久久久久国内视频| 伊人久久大香线蕉亚洲五| 国产精品亚洲av一区麻豆| 黄色女人牲交| 看片在线看免费视频| 精品免费久久久久久久清纯| 亚洲熟女毛片儿| 99久久国产精品久久久| 一级a爱片免费观看的视频| 天堂影院成人在线观看| 午夜福利免费观看在线| 老熟妇乱子伦视频在线观看| 欧美色视频一区免费| 午夜影院日韩av| 久久伊人香网站| 久久久水蜜桃国产精品网| 国产伦一二天堂av在线观看| 18禁裸乳无遮挡免费网站照片| 男人舔女人的私密视频| 精品99又大又爽又粗少妇毛片 | 久久精品亚洲精品国产色婷小说| 色精品久久人妻99蜜桃| 亚洲 国产 在线| 成人三级做爰电影| 长腿黑丝高跟| 国产精品乱码一区二三区的特点| 久久性视频一级片| 久久久国产欧美日韩av| 真实男女啪啪啪动态图| 最近视频中文字幕2019在线8| 色老头精品视频在线观看| 欧美成狂野欧美在线观看| 2021天堂中文幕一二区在线观| 欧美国产日韩亚洲一区| 亚洲av日韩精品久久久久久密| 欧美日韩瑟瑟在线播放| 99在线视频只有这里精品首页| 亚洲欧美日韩无卡精品| 日本一本二区三区精品| 我要搜黄色片| 国产午夜福利久久久久久| av片东京热男人的天堂| 国产欧美日韩一区二区三| 看黄色毛片网站| 久久亚洲真实| 伊人久久大香线蕉亚洲五| 国产激情久久老熟女| 国产黄a三级三级三级人| 不卡av一区二区三区| 给我免费播放毛片高清在线观看| 久久国产乱子伦精品免费另类| 最新美女视频免费是黄的| 天天添夜夜摸| 熟女电影av网| 欧美丝袜亚洲另类 | 欧美高清成人免费视频www| 99久久无色码亚洲精品果冻| 男人舔女人的私密视频| 日本三级黄在线观看| 成人特级av手机在线观看| 亚洲色图av天堂| bbb黄色大片| 一a级毛片在线观看| 可以在线观看毛片的网站| 69av精品久久久久久| 黄色 视频免费看| 少妇的丰满在线观看| 日韩欧美三级三区| 两性午夜刺激爽爽歪歪视频在线观看| 久久性视频一级片| 日本免费一区二区三区高清不卡| 熟女人妻精品中文字幕| 波多野结衣高清无吗| 精品久久蜜臀av无| 一区福利在线观看| 观看免费一级毛片| АⅤ资源中文在线天堂| 波多野结衣巨乳人妻| 久久久精品欧美日韩精品| 国产亚洲精品综合一区在线观看| 国产精品精品国产色婷婷| 老司机深夜福利视频在线观看| 久久久久久久精品吃奶| 麻豆成人av在线观看| 两性夫妻黄色片| 午夜视频精品福利| 97超级碰碰碰精品色视频在线观看| 成人三级黄色视频| 欧美中文综合在线视频| 色吧在线观看| 成人欧美大片| 国产精品av久久久久免费| 成人高潮视频无遮挡免费网站| 亚洲国产欧美一区二区综合| 日韩中文字幕欧美一区二区| 别揉我奶头~嗯~啊~动态视频| 国产精品久久久久久亚洲av鲁大| 欧美中文日本在线观看视频| 久久亚洲精品不卡| 国产成人av激情在线播放| 成人特级av手机在线观看| 国产伦在线观看视频一区| 啪啪无遮挡十八禁网站| 亚洲成人免费电影在线观看| 成人性生交大片免费视频hd| 人人妻人人澡欧美一区二区| 色av中文字幕| 黄片小视频在线播放| 中文字幕av在线有码专区| 琪琪午夜伦伦电影理论片6080| 三级男女做爰猛烈吃奶摸视频| 成年免费大片在线观看| 天堂动漫精品| 精品福利观看| 日本黄色片子视频| 天堂√8在线中文| 婷婷精品国产亚洲av| 亚洲成a人片在线一区二区| 亚洲18禁久久av| 亚洲最大成人中文| 午夜激情欧美在线| 岛国视频午夜一区免费看| 欧美激情久久久久久爽电影| 在线免费观看不下载黄p国产 | 国产熟女xx| 精品欧美国产一区二区三| 热99在线观看视频| 国产亚洲精品一区二区www| 一本一本综合久久| 国产精品av视频在线免费观看| 91九色精品人成在线观看| 好男人在线观看高清免费视频| 国产高清激情床上av| 亚洲一区二区三区不卡视频| 精品一区二区三区视频在线 | 午夜福利视频1000在线观看| 国产精品久久久av美女十八| 真人做人爱边吃奶动态| 又大又爽又粗| 国产一区二区在线观看日韩 | 在线观看66精品国产| 国产成人欧美在线观看| 久久久久久久午夜电影| 欧美三级亚洲精品| 法律面前人人平等表现在哪些方面| 久久中文字幕人妻熟女| 老熟妇仑乱视频hdxx| 最新美女视频免费是黄的| 日韩av在线大香蕉| 日本与韩国留学比较|