• <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ī)則》的決定
    91字幕亚洲| 首页视频小说图片口味搜索| 窝窝影院91人妻| 九色国产91popny在线| 色综合欧美亚洲国产小说| 亚洲第一电影网av| 美女被艹到高潮喷水动态| 手机成人av网站| 白带黄色成豆腐渣| 三级男女做爰猛烈吃奶摸视频| 色老头精品视频在线观看| 国产激情偷乱视频一区二区| 少妇的逼水好多| 色综合婷婷激情| 99视频精品全部免费 在线| 国产野战对白在线观看| 午夜a级毛片| 国产成人福利小说| 天美传媒精品一区二区| 国产精品 国内视频| av女优亚洲男人天堂| 日韩有码中文字幕| 91麻豆av在线| 久久久久九九精品影院| 丰满乱子伦码专区| 国产一区二区在线av高清观看| 欧美性猛交╳xxx乱大交人| 深爱激情五月婷婷| 一本久久中文字幕| 中文字幕人妻熟人妻熟丝袜美 | 十八禁网站免费在线| www.www免费av| 精品99又大又爽又粗少妇毛片 | 精品国内亚洲2022精品成人| 一进一出抽搐gif免费好疼| 亚洲,欧美精品.| 国产黄色小视频在线观看| 日本 欧美在线| 偷拍熟女少妇极品色| 日韩av在线大香蕉| 欧美性猛交黑人性爽| 他把我摸到了高潮在线观看| 国产不卡一卡二| 亚洲中文字幕一区二区三区有码在线看| 欧美日韩综合久久久久久 | 日本 av在线| 亚洲中文日韩欧美视频| 久久久精品大字幕| 国产成人欧美在线观看| 舔av片在线| 精华霜和精华液先用哪个| 精品国产三级普通话版| 精品国产美女av久久久久小说| 十八禁网站免费在线| 免费高清视频大片| 亚洲精品成人久久久久久| 亚洲性夜色夜夜综合| 国产麻豆成人av免费视频| 免费观看精品视频网站| 久久欧美精品欧美久久欧美| 亚洲av第一区精品v没综合| 又爽又黄无遮挡网站| 欧美日本视频| 五月伊人婷婷丁香| 日本 欧美在线| 亚洲国产中文字幕在线视频| 一区二区三区激情视频| 九九热线精品视视频播放| 成人性生交大片免费视频hd| 免费人成在线观看视频色| 国产极品精品免费视频能看的| 美女黄网站色视频| 国产淫片久久久久久久久 | 久久久久久久久中文| 大型黄色视频在线免费观看| 免费大片18禁| 成人18禁在线播放| 精品不卡国产一区二区三区| 天堂√8在线中文| 麻豆国产av国片精品| 国内精品久久久久精免费| 色哟哟哟哟哟哟| 深夜精品福利| 成人av在线播放网站| 99热精品在线国产| 久久欧美精品欧美久久欧美| 日韩av在线大香蕉| 国产精品一区二区免费欧美| 日韩人妻高清精品专区| 日本免费a在线| 性欧美人与动物交配| 成人永久免费在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 69av精品久久久久久| 精品久久久久久成人av| av福利片在线观看| 国产一区二区亚洲精品在线观看| 两人在一起打扑克的视频| 国产亚洲精品综合一区在线观看| 国内揄拍国产精品人妻在线| 日本 欧美在线| 国产精品永久免费网站| 在线观看免费视频日本深夜| 90打野战视频偷拍视频| 国产v大片淫在线免费观看| 国内精品一区二区在线观看| avwww免费| 欧美黑人巨大hd| 黄色丝袜av网址大全| 一区二区三区激情视频| 午夜福利免费观看在线| 最新在线观看一区二区三区| 欧美日韩精品网址| 天堂动漫精品| aaaaa片日本免费| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久久大av| 美女cb高潮喷水在线观看| 国产在线精品亚洲第一网站| 中文亚洲av片在线观看爽| 日本 av在线| 在线视频色国产色| 国产黄片美女视频| 午夜日韩欧美国产| 成年版毛片免费区| 国产三级中文精品| 亚洲18禁久久av| 中亚洲国语对白在线视频| 国产精品嫩草影院av在线观看 | 久久精品国产综合久久久| 精品国产亚洲在线| 午夜免费成人在线视频| 久久久久久人人人人人| 99久久综合精品五月天人人| 国产美女午夜福利| 麻豆成人午夜福利视频| 可以在线观看的亚洲视频| 性欧美人与动物交配| 熟女少妇亚洲综合色aaa.| 99在线视频只有这里精品首页| 精品熟女少妇八av免费久了| 色尼玛亚洲综合影院| 亚洲欧美日韩高清在线视频| 国产精品久久久人人做人人爽| 亚洲狠狠婷婷综合久久图片| 午夜老司机福利剧场| 美女高潮的动态| 欧美另类亚洲清纯唯美| 欧美不卡视频在线免费观看| 中文字幕人成人乱码亚洲影| 久久6这里有精品| 成年女人毛片免费观看观看9| 脱女人内裤的视频| 一个人看的www免费观看视频| 亚洲中文字幕日韩| 无遮挡黄片免费观看| 桃红色精品国产亚洲av| 午夜老司机福利剧场| 色av中文字幕| 亚洲第一电影网av| 少妇熟女aⅴ在线视频| 亚洲真实伦在线观看| 亚洲av美国av| 欧美黑人巨大hd| 亚洲av不卡在线观看| 国产伦一二天堂av在线观看| 日本黄大片高清| 国产91精品成人一区二区三区| 日本a在线网址| 久久久久性生活片| 看黄色毛片网站| 国产高清视频在线播放一区| 国产av不卡久久| 午夜精品在线福利| 中文在线观看免费www的网站| 天天一区二区日本电影三级| 成人av一区二区三区在线看| 97碰自拍视频| 亚洲国产精品sss在线观看| 天堂√8在线中文| 天天躁日日操中文字幕| 亚洲av第一区精品v没综合| 三级毛片av免费| 国产aⅴ精品一区二区三区波| 黄色日韩在线| 在线观看一区二区三区| 在线免费观看的www视频| 母亲3免费完整高清在线观看| 三级毛片av免费| 中文在线观看免费www的网站| 老司机在亚洲福利影院| 在线观看一区二区三区| 欧美成人性av电影在线观看| av女优亚洲男人天堂| 国内久久婷婷六月综合欲色啪| 成人欧美大片| 精品福利观看| 欧美又色又爽又黄视频| 国产三级黄色录像| 欧美色欧美亚洲另类二区| 美女高潮的动态| 亚洲第一电影网av| 精品国产亚洲在线| 男女床上黄色一级片免费看| 欧美在线一区亚洲| 美女高潮的动态| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 免费观看人在逋| 国产成人影院久久av| 亚洲av电影不卡..在线观看| 免费av不卡在线播放| 女人高潮潮喷娇喘18禁视频| 亚洲国产欧美人成| 女人被狂操c到高潮| 免费无遮挡裸体视频| 欧美黄色片欧美黄色片| bbb黄色大片| 中文字幕av在线有码专区| 国产真实伦视频高清在线观看 | 免费无遮挡裸体视频| 欧美黄色片欧美黄色片| 午夜免费激情av| 国内久久婷婷六月综合欲色啪| 国产精品1区2区在线观看.| 精品一区二区三区视频在线观看免费| 亚洲av成人不卡在线观看播放网| 一个人免费在线观看电影| 国内精品久久久久久久电影| 99久久综合精品五月天人人| 免费观看精品视频网站| av天堂中文字幕网| 国产免费男女视频| 性欧美人与动物交配| 亚洲国产精品久久男人天堂| 亚洲国产日韩欧美精品在线观看 | 午夜福利视频1000在线观看| 国产视频内射| 欧美高清成人免费视频www| 国产精品av视频在线免费观看| 午夜精品在线福利| 九九在线视频观看精品| 亚洲av电影在线进入| 欧洲精品卡2卡3卡4卡5卡区| 色在线成人网| 亚洲国产色片| 国产精品亚洲一级av第二区| 国产精品精品国产色婷婷| 午夜久久久久精精品| 久久久久亚洲av毛片大全| 在线观看午夜福利视频| 国产亚洲精品综合一区在线观看| 亚洲18禁久久av| 亚洲欧美精品综合久久99| 日本成人三级电影网站| 岛国在线免费视频观看| 18+在线观看网站| 18禁国产床啪视频网站| 亚洲国产中文字幕在线视频| 亚洲精品一区av在线观看| 成人无遮挡网站| 欧美乱码精品一区二区三区| 国产精品三级大全| 国产高清三级在线| 男人舔女人下体高潮全视频| 亚洲一区高清亚洲精品| 日本免费一区二区三区高清不卡| 国产亚洲av嫩草精品影院| 国产精品乱码一区二三区的特点| 黄片小视频在线播放| 精品欧美国产一区二区三| 99热这里只有精品一区| 国产欧美日韩精品一区二区| 91久久精品电影网| 国产伦精品一区二区三区视频9 | 国产精品女同一区二区软件 | 国产高清有码在线观看视频| 亚洲性夜色夜夜综合| 国产视频一区二区在线看| 手机成人av网站| 亚洲精品日韩av片在线观看 | 一本精品99久久精品77| 精品一区二区三区人妻视频| 天美传媒精品一区二区| 亚洲欧美精品综合久久99| 亚洲国产色片| 九色国产91popny在线| 小蜜桃在线观看免费完整版高清| 亚洲成人免费电影在线观看| 亚洲美女黄片视频| 搡老妇女老女人老熟妇| 狂野欧美激情性xxxx| 国产熟女xx| 91av网一区二区| 亚洲av电影在线进入| 国产又黄又爽又无遮挡在线| 他把我摸到了高潮在线观看| 色尼玛亚洲综合影院| 啦啦啦免费观看视频1| 非洲黑人性xxxx精品又粗又长| 老汉色∧v一级毛片| 嫩草影院入口| 国产乱人伦免费视频| 天美传媒精品一区二区| 又黄又粗又硬又大视频| 一个人免费在线观看电影| 日韩欧美在线乱码| 嫩草影院精品99| 丁香六月欧美| 丰满乱子伦码专区| 久久久久久久久中文| 九色成人免费人妻av| 国产色爽女视频免费观看| 国产主播在线观看一区二区| 亚洲精品乱码久久久v下载方式 | 99热精品在线国产| 97超视频在线观看视频| 国产成人a区在线观看| 欧美午夜高清在线| 国产又黄又爽又无遮挡在线| 黄片小视频在线播放| 欧洲精品卡2卡3卡4卡5卡区| 国内揄拍国产精品人妻在线| 欧美另类亚洲清纯唯美| 国产精品香港三级国产av潘金莲| 欧美极品一区二区三区四区| 啦啦啦免费观看视频1| 国产亚洲欧美98| 国产爱豆传媒在线观看| 好男人在线观看高清免费视频| 小说图片视频综合网站| 午夜福利在线观看吧| 又黄又爽又免费观看的视频| 岛国在线免费视频观看| 久久精品国产亚洲av涩爱 | 蜜桃亚洲精品一区二区三区| 亚洲乱码一区二区免费版| 18禁黄网站禁片午夜丰满| 麻豆久久精品国产亚洲av| 精品久久久久久久末码| 变态另类成人亚洲欧美熟女| 免费看十八禁软件| 黄色成人免费大全| 亚洲精品久久国产高清桃花| 又紧又爽又黄一区二区| www.www免费av| 黄色成人免费大全| 日韩欧美精品免费久久 | 黄色女人牲交| 欧美极品一区二区三区四区| 亚洲片人在线观看| 丰满人妻熟妇乱又伦精品不卡| av欧美777| 国产伦精品一区二区三区四那| 五月玫瑰六月丁香| 国产精品电影一区二区三区| 九色成人免费人妻av| 免费在线观看日本一区| 亚洲精品国产精品久久久不卡| 精品久久久久久久久久免费视频| 亚洲av免费在线观看| 中文在线观看免费www的网站| 他把我摸到了高潮在线观看| 国产精品一区二区免费欧美| 天堂√8在线中文| 精品人妻一区二区三区麻豆 | 久久人妻av系列| 午夜日韩欧美国产| 麻豆国产97在线/欧美| 99热这里只有精品一区| 亚洲精品一区av在线观看| 国产精品久久久人人做人人爽| 成人性生交大片免费视频hd| av在线蜜桃| 成年人黄色毛片网站| 成年女人看的毛片在线观看| 一本久久中文字幕| 搞女人的毛片| 一个人观看的视频www高清免费观看| 国产精品国产高清国产av| 丰满人妻一区二区三区视频av | 午夜福利成人在线免费观看| 在线观看av片永久免费下载| 久久久久性生活片| 国产精品亚洲一级av第二区| 久久久久国内视频| 久久精品91蜜桃| 男人和女人高潮做爰伦理| 又爽又黄无遮挡网站| 一级a爱片免费观看的视频| 亚洲第一欧美日韩一区二区三区| 美女大奶头视频| 毛片女人毛片| 每晚都被弄得嗷嗷叫到高潮| 非洲黑人性xxxx精品又粗又长| 精品人妻1区二区| 嫁个100分男人电影在线观看| 精品国产超薄肉色丝袜足j| 久久九九热精品免费| 欧美日本视频| 亚洲欧美日韩卡通动漫| 欧美日韩精品网址| www.色视频.com| 精品久久久久久久久久久久久| 日韩欧美一区二区三区在线观看| 男女那种视频在线观看| 啦啦啦免费观看视频1| 校园春色视频在线观看| 国产私拍福利视频在线观看| 又爽又黄无遮挡网站| 日本一本二区三区精品| 亚洲国产精品久久男人天堂| 欧美日韩一级在线毛片| 欧美中文日本在线观看视频| 免费观看精品视频网站| 国产免费av片在线观看野外av| 色精品久久人妻99蜜桃| 国产亚洲精品综合一区在线观看| 午夜免费男女啪啪视频观看 | 18禁美女被吸乳视频| 色尼玛亚洲综合影院| 一区二区三区高清视频在线| 久久久久久久久大av| 国产成人啪精品午夜网站| 久久精品国产亚洲av涩爱 | 一级作爱视频免费观看| 日日摸夜夜添夜夜添小说| 精品不卡国产一区二区三区| 18禁美女被吸乳视频| 精品国产美女av久久久久小说| 成人性生交大片免费视频hd| 操出白浆在线播放| 国产精品亚洲一级av第二区| 国产一区二区激情短视频| 最新美女视频免费是黄的| 国产亚洲精品久久久久久毛片| 搡老妇女老女人老熟妇| 亚洲一区高清亚洲精品| 日韩欧美 国产精品| 久久6这里有精品| 国产97色在线日韩免费| 亚洲成人久久性| 亚洲,欧美精品.| 亚洲久久久久久中文字幕| av专区在线播放| 天堂网av新在线| 久久久成人免费电影| 成年人黄色毛片网站| 在线观看免费视频日本深夜| 国产三级黄色录像| 免费人成在线观看视频色| 91麻豆精品激情在线观看国产| bbb黄色大片| 国产一区二区三区视频了| 五月玫瑰六月丁香| 国产精品久久电影中文字幕| 夜夜夜夜夜久久久久| 蜜桃久久精品国产亚洲av| 又黄又粗又硬又大视频| 欧美黑人欧美精品刺激| 最新美女视频免费是黄的| 女同久久另类99精品国产91| 麻豆成人av在线观看| 少妇裸体淫交视频免费看高清| 99久久综合精品五月天人人| 好男人电影高清在线观看| 国模一区二区三区四区视频| 精品一区二区三区av网在线观看| 搡老岳熟女国产| 免费搜索国产男女视频| 中文字幕av在线有码专区| 色在线成人网| 51国产日韩欧美| 在线十欧美十亚洲十日本专区| 免费看十八禁软件| 757午夜福利合集在线观看| 观看美女的网站| 久久人妻av系列| 成人特级av手机在线观看| 可以在线观看毛片的网站| 久久久久久大精品| 激情在线观看视频在线高清| 99国产综合亚洲精品| 老司机在亚洲福利影院| 别揉我奶头~嗯~啊~动态视频| 欧美高清成人免费视频www| 91久久精品国产一区二区成人 | 日韩成人在线观看一区二区三区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 又粗又爽又猛毛片免费看| 亚洲五月天丁香| 成年版毛片免费区| 啪啪无遮挡十八禁网站| 日本一本二区三区精品| 午夜福利欧美成人| 精品一区二区三区av网在线观看| 免费大片18禁| 99视频精品全部免费 在线| 夜夜爽天天搞| 久久国产精品影院| 一级黄片播放器| 中文字幕av在线有码专区| 亚洲av一区综合| av天堂中文字幕网| 亚洲欧美一区二区三区黑人| 国产一级毛片七仙女欲春2| 色视频www国产| av天堂在线播放| 精品久久久久久成人av| 九色成人免费人妻av| 国产精品乱码一区二三区的特点| 黄片大片在线免费观看| 欧美日韩精品网址| 国产美女午夜福利| 热99在线观看视频| 欧美一级a爱片免费观看看| 国产高清激情床上av| 丁香六月欧美| 啪啪无遮挡十八禁网站| 日本三级黄在线观看| 热99在线观看视频| 成人精品一区二区免费| 中亚洲国语对白在线视频| 亚洲精品乱码久久久v下载方式 | 国产中年淑女户外野战色| 久久久久久九九精品二区国产| 午夜福利成人在线免费观看| 九九在线视频观看精品| 十八禁人妻一区二区| 一级毛片女人18水好多| 亚洲天堂国产精品一区在线| 两个人的视频大全免费| 久久中文看片网| 最新在线观看一区二区三区| 久久精品人妻少妇| 午夜亚洲福利在线播放| 少妇熟女aⅴ在线视频| 99国产精品一区二区蜜桃av| 午夜免费男女啪啪视频观看 | 欧美最黄视频在线播放免费| 久久久久久久午夜电影| 婷婷亚洲欧美| 国产精品精品国产色婷婷| 亚洲成人免费电影在线观看| 亚洲一区高清亚洲精品| 精华霜和精华液先用哪个| 免费在线观看日本一区| 宅男免费午夜| 啦啦啦韩国在线观看视频| 欧美在线一区亚洲| 免费无遮挡裸体视频| 久久久色成人| 男女下面进入的视频免费午夜| 亚洲成av人片免费观看| 3wmmmm亚洲av在线观看| 国产探花在线观看一区二区| 一个人观看的视频www高清免费观看| 亚洲无线观看免费| 午夜精品在线福利| 色噜噜av男人的天堂激情| 精品久久久久久久人妻蜜臀av| 国产欧美日韩精品亚洲av| 国产精品一区二区三区四区免费观看 | 久久久久精品国产欧美久久久| 国产亚洲欧美在线一区二区| 99国产精品一区二区三区| av女优亚洲男人天堂| 国产乱人视频| 亚洲男人的天堂狠狠| 亚洲精品色激情综合| 日日夜夜操网爽| 日韩av在线大香蕉| 久久精品人妻少妇| 少妇的丰满在线观看| 亚洲精品色激情综合| 久久久精品大字幕| 免费人成视频x8x8入口观看| 日韩欧美 国产精品| 久99久视频精品免费| 又粗又爽又猛毛片免费看| aaaaa片日本免费| 亚洲国产精品sss在线观看| 免费人成在线观看视频色| 国产精品 国内视频| 2021天堂中文幕一二区在线观| 国模一区二区三区四区视频| 精品久久久久久,| 中文字幕久久专区| 中亚洲国语对白在线视频| 色在线成人网| 国产欧美日韩精品一区二区| 听说在线观看完整版免费高清| 午夜两性在线视频| 特大巨黑吊av在线直播| 午夜福利免费观看在线| 精品国产超薄肉色丝袜足j| 可以在线观看的亚洲视频| 午夜a级毛片| 91九色精品人成在线观看| 1024手机看黄色片| 国产高潮美女av| 国产av一区在线观看免费| 日本免费a在线| 波野结衣二区三区在线 | 在线看三级毛片| 亚洲av免费在线观看| 白带黄色成豆腐渣| 日韩欧美精品免费久久 | 在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 淫秽高清视频在线观看| 色综合婷婷激情| tocl精华| 欧美+亚洲+日韩+国产| 久久久久久久亚洲中文字幕 | 欧美日韩中文字幕国产精品一区二区三区|