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

    基于不同時間尺度風(fēng)切變指數(shù)推算風(fēng)資源的對比分析

    2022-11-05 13:44:14徐紅龔強(qiáng)朱玲沈歷都顧正強(qiáng)晁華王乙舒陳軍慶
    氣象與環(huán)境學(xué)報 2022年5期
    關(guān)鍵詞:風(fēng)塔測風(fēng)塔探空

    徐紅 龔強(qiáng) 朱玲 沈歷都 顧正強(qiáng) 晁華 王乙舒 陳軍慶

    (1.沈陽區(qū)域氣候中心,遼寧沈陽 110166;2.中國科學(xué)院沈陽應(yīng)用生態(tài)研究所,遼寧沈陽 110016;3.空軍 93161部隊,遼寧沈陽 110015)

    引言

    在常規(guī)能源短缺和環(huán)境污染的雙重壓力下,風(fēng)能資源因其自身獨(dú)有的優(yōu)點(diǎn),在很短時間內(nèi)便獲得巨大發(fā)展。近年來,隨著風(fēng)電事業(yè)逐步發(fā)展,中國風(fēng)電總裝機(jī)容量逐步增大,風(fēng)能可利用區(qū)域逐步縮減。隨著材料、工藝、運(yùn)輸?shù)仁侄蔚牟粩噙M(jìn)步,企業(yè)常用增加葉輪直徑、提升輪轂高度等技術(shù)提升效益。近年來,低風(fēng)速區(qū)域測風(fēng)塔高度率先向120—150 m發(fā)展。因風(fēng)切變特征不同,同種風(fēng)電機(jī)組提升相同高度所增加的發(fā)電量也不相同,一般情況下,風(fēng)切變指數(shù)越大,發(fā)電量提升越大。但是隨著高度的提升,高度對氣流速度的影響逐漸變小,導(dǎo)致風(fēng)速因高度提升造成的增幅逐漸變小,最終必然導(dǎo)致風(fēng)電機(jī)組發(fā)電量增幅隨著高度的提升而減少。同時,塔筒高度的提升造成制造、施工的難度也迅速上升。發(fā)電量收益增幅的降低和工程投入增幅的變大,必然使風(fēng)電機(jī)組塔筒高度存在一個可計算的平衡值,在這個平衡值上風(fēng)電場將實(shí)現(xiàn)發(fā)電量收益的最大化。

    風(fēng)切變指數(shù)受大氣穩(wěn)定度、下墊面屬性、天氣活動等因素影響,在不同時間尺度上呈一定程度的波動變化[1-5]。在建設(shè)風(fēng)電場前往往需要建立測風(fēng)塔,進(jìn)行為期一年的完整觀測,對風(fēng)電場風(fēng)資源狀況進(jìn)行前期評估??紤]測風(fēng)塔的建設(shè)成本以及主流風(fēng)機(jī)輪轂高度,一些風(fēng)電場早期設(shè)立測風(fēng)塔達(dá)不到現(xiàn)在輪轂高度要求,亟需在現(xiàn)有觀測基礎(chǔ)上對更高輪轂高度風(fēng)資源進(jìn)行合理的推算。風(fēng)資源推算方法有多種[6-9],常用低層風(fēng)速和風(fēng)切變指數(shù)推算高層風(fēng)資源狀況,推算結(jié)果直接決定了風(fēng)機(jī)選型、發(fā)電量估算等。為評估推算精確度,本文對遼寧省風(fēng)切變特征進(jìn)行統(tǒng)計分析,并采用多種風(fēng)切變指數(shù)計算方案,對無觀測的高層風(fēng)資源進(jìn)行推算,并與實(shí)測風(fēng)對比分析,旨在確定一種較好的方案,提高風(fēng)資源評估工作中風(fēng)資源推算精度。

    1 資料與方法

    1.1 資料來源

    中國目前共有119個長期觀測的探空站,其中遼寧省在沈陽和大連各有1個。沈陽探空站位于遼寧中部,為遼中平原和遼東丘陵交界處。大連探空站位于遼東半島南端沿海地區(qū)。探空站每日08時和20時定時觀測,觀測高度從地面至高空約30 km。取2014—2018年秒級探空數(shù)據(jù),包括溫度、氣壓、風(fēng)速、風(fēng)向、高度等。考慮對風(fēng)資源貢獻(xiàn),剔除300 m高度以下各觀測高度(探空儀器隨時間和地點(diǎn)不同,每秒上升高度不同)平均風(fēng)速≤1.5 m·s-1、最大風(fēng)速≤3.0 m·s-1或最大風(fēng)速超過40 m·s-1的觀測數(shù)值,剩余數(shù)據(jù)約占總數(shù)據(jù)量的80%。

    參照《全國風(fēng)能資源詳查和評價報告》(全國400座測風(fēng)塔,其中遼寧省26座測風(fēng)塔分布比較均勻)可知遼寧省風(fēng)能資源比較豐富[10]。龔強(qiáng)等[11]對遼寧省26座測風(fēng)塔風(fēng)切變以及企業(yè)測風(fēng)塔觀測結(jié)果研究表明,遼寧省風(fēng)切變有較為一致的日變化特征,故綜合考慮地理位置、地形及測風(fēng)高度等,選用2座代表性測風(fēng)塔,分別為1#和2#。測風(fēng)塔1#位于遼東半島,屬于東部沿海風(fēng)能資源豐富帶,適合風(fēng)能資源開發(fā)。測風(fēng)塔2#位于遼中北部平原到丘陵過渡地帶,為我國“三北”風(fēng)帶,是我國風(fēng)能資源豐富地區(qū)之一,測風(fēng)塔2#代表風(fēng)資源較好的丘陵地區(qū)。雖然遼寧省地表有山地、平原、海岸線等不同地形地貌,地面粗糙度不同,但是在測風(fēng)塔數(shù)據(jù)有限的情況下,本文選用的兩測風(fēng)塔對遼寧省還是相對具有一定的代表性。

    2座測風(fēng)塔均有5層風(fēng)速觀測。測風(fēng)塔1#(觀測時段為2010年6月至2011年5月)觀測高度分別為10 m(第1層)、30 m(第2層)、50 m(第3層)、60 m(第4層)、70 m(第5層),測風(fēng)塔2#(觀測時段為2016年2月至2017年1月)觀測高度分別為30 m(第1層)、70 m(第2層)、80 m(第3層)、90 m(第4層)、100 m(第5層)。觀測數(shù)據(jù)完整率高、真實(shí)可靠。原始數(shù)據(jù)為10 min觀測,經(jīng)處理變?yōu)橥暾晷r數(shù)據(jù),并參照就近氣象觀測站歷史風(fēng)序列[12]進(jìn)行測風(fēng)年度代表性分析與訂正,能較好地反映常年平均風(fēng)資源狀況。探空站和測風(fēng)塔位置見圖1。

    圖1 探空站和測風(fēng)塔位置Fig.1 Location of the tower station and sounding balloon

    1.2 研究方法

    雖然探空站每天僅有08時和20時觀測,但也可以作為垂直背景風(fēng)參考。

    在中性大氣層結(jié)下,對數(shù)和冪指數(shù)方程都可以較好地描述風(fēng)速垂直廓線,為綜合反映近地層風(fēng)切變狀況,選用冪指數(shù)擬合方法確定風(fēng)切變指數(shù)[13-18],風(fēng)隨高度變化公式如下

    風(fēng)切變指數(shù)α計算公式如下

    式(1)~式(2)中,α為風(fēng)切變指數(shù);Z為觀測高度;VZ為對應(yīng)高度的風(fēng)速;Zr為計算綜合切變指數(shù)的參考高度;Vr為參考高度的風(fēng)速。

    2 結(jié)果分析

    2.1 背景風(fēng)分析

    風(fēng)受多種因素影響,比如海陸屬性差異、地形、天氣過程等。選取2014—2018年共5 a的探空風(fēng)數(shù)據(jù),從個例分析到年、月尺度進(jìn)行統(tǒng)計分析??傮w來看,2014—2018年平均風(fēng)廓線比較平滑。從全國119個探空站平均垂直風(fēng)廓線統(tǒng)計結(jié)果來看(圖略),各探空站風(fēng)速隨高度增加的快慢不同,風(fēng)速平穩(wěn)不再增加的高度也不相同,有些站還會出現(xiàn)超低空急流,即風(fēng)速先增大再減小。圖2a為沈陽和大連探空站平均風(fēng)廓線??紤]地面50 m高度以下受下墊面影響較大,不再繪制對應(yīng)高度風(fēng)廓線,探空站風(fēng)廓線分析高度上至300 m。通過年平均風(fēng)廓線可以看出,沈陽站風(fēng)速隨高度增加較快,大約至180 m高度風(fēng)速基本保持平穩(wěn)不變。大連站風(fēng)廓線比較平直,風(fēng)速隨高度增加較慢,高低層風(fēng)速差異小,風(fēng)速無明顯的轉(zhuǎn)折高度。150 m高度以下大連站風(fēng)速大于沈陽站,150~300 m高度風(fēng)速小于沈陽地區(qū)。即遼寧省180 m高度以下基本可以采用同一種風(fēng)廓線擬合方式。各季節(jié)代表月的風(fēng)廓線顯示(圖2b和圖2c),沈陽站冬季和夏季風(fēng)廓線差異較小,春秋季風(fēng)速明顯增大。與沈陽站相比,大連站緯度偏低且受海洋性氣候影響,四季風(fēng)廓線差異偏小,除夏季外,其他季節(jié)差異相對更小。

    圖2 2014—2018年沈陽站和大連站年平均風(fēng)廓線(a)及沈陽站(b)和大連站(c)代表月平均風(fēng)廓線Fig.2 Annual wind profiles(a)and typically month wind profiles at Shenyang(b)and Dalian(c)during 2014-2018

    2.2 不同高度組合月風(fēng)切變指數(shù)擬合分析

    表1為測風(fēng)塔逐月平均風(fēng)速,各測風(fēng)塔不同高度風(fēng)速逐月變化趨勢較為一致,2座測風(fēng)塔相同高度風(fēng)速大小有一定的差異??傮w來看,2座測風(fēng)塔風(fēng)速逐月變化特征符合遼寧省風(fēng)分布特征。圖3為采用表1中測風(fēng)塔不同高度組合的月平均風(fēng)速擬合的風(fēng)切變指數(shù)圖。從風(fēng)切變指數(shù)月變化趨勢可以看出,月風(fēng)切變指數(shù)與下墊面屬性、風(fēng)速等變化有關(guān),雖然兩測風(fēng)塔風(fēng)切變大小有較大的差異,但存在較為一致的月變化特征,春季風(fēng)切變最小,夏季風(fēng)切變最大。測風(fēng)塔1#去掉最低層10 m觀測擬合風(fēng)切變指數(shù)范圍為0.00~0.12,其余方案風(fēng)切變指數(shù)范圍為0.02~0.08。測風(fēng)塔2#去掉最低層30 m高度的風(fēng)速擬合風(fēng)切變指數(shù)范圍為0.18~0.33,其余方案風(fēng)切變指數(shù)范圍為0.18~0.35,相比于測風(fēng)塔1#各方案差異偏小。

    表1 測風(fēng)塔1#、2#月平均風(fēng)速Table 1 Monthly wind speed at Tower 1 and Tower 2 m·s-1

    圖3 測風(fēng)塔1#(a)和2#(b)不同高度組合月平均風(fēng)切變指數(shù)Fig.3 Monthly wind shear indices of different height combinations at Tower 1#(a)and Tower 2#(b)

    2.3 不同高度組合日風(fēng)切變指數(shù)擬合分析

    由于風(fēng)有晝夜、月變化特征,對比風(fēng)切變指數(shù)日、月變化特征差異,一是為了分析不同時間尺度的風(fēng)切變變化特征;二是為采用不同風(fēng)切變指數(shù)擬合高層風(fēng)速和風(fēng)功率密度做鋪墊和解釋。表2為測風(fēng)塔逐時平均風(fēng)速,各測風(fēng)塔不同高度風(fēng)速逐時變化趨勢較為一致。從測風(fēng)塔各層風(fēng)速對比來看,底層風(fēng)速日變化最為顯著,隨著高度的增加,風(fēng)速日變化趨于減小。圖4為采用不同高度組合的小時平均風(fēng)速擬合風(fēng)切變指數(shù)圖。由圖4可以看出,風(fēng)切變指數(shù)日變化特征比月變化顯著。測風(fēng)塔1#去掉10 m高度的風(fēng)速擬合風(fēng)切變指數(shù)與其他方案有明顯差異,去掉最低層10 m觀測擬合風(fēng)切變指數(shù)范圍為-0.04~0.17,其余方案風(fēng)切變指數(shù)范圍為-0.01~0.11。測風(fēng)塔2#去掉最低層30 m高度的風(fēng)速擬合風(fēng)切變指數(shù)與其他方案有一定的差異,去掉最低層30 m觀測擬合風(fēng)切變指數(shù)范圍為0.04~0.40,其余方案風(fēng)切變指數(shù)范圍為0.08~0.39,相比于測風(fēng)塔1#差異偏小。

    表2 測風(fēng)塔1#、2#逐時平均風(fēng)速Table 2 Mean hourly wind speed at Tower 1#(a)and Tower 2#(b) m·s-1

    圖4 測風(fēng)塔1#(a)、2#(b)不同高度組合小時平均風(fēng)切變指數(shù)Fig.4 Hourly wind shear indices of different height combinations at Tower 1#(a)and Tower 2#(b)

    從2座測風(fēng)塔月和日變化對比可以看出,測風(fēng)塔1#去掉10 m高度擬合方案與其他方案有明顯差異,因?yàn)?0 m高度風(fēng)受下墊面影響較大,對擬合風(fēng)切變指數(shù)有顯著的影響。測風(fēng)塔2#最低層30 m高度風(fēng)速受下墊面影響明顯減小,對擬合風(fēng)切變指數(shù)影響也減小。

    結(jié)合表2和圖4可以看出,雖然2座測風(fēng)塔風(fēng)切變指數(shù)差異較大,但其變化趨勢表現(xiàn)出高度的一致性。大約10時開始,太陽輻射逐漸增強(qiáng),地表溫度升高,容易形成較強(qiáng)的不穩(wěn)定層結(jié),有利于高層大氣的動量下傳,上下層空氣動量交換造成了較小高低層風(fēng)速差異,風(fēng)切變較小。16時開始,太陽輻射減弱,大氣逐漸形成穩(wěn)定層結(jié),加上下墊面粗糙拖曳作用,上下層風(fēng)速差異逐漸增大,風(fēng)切變指數(shù)也逐漸增大,直到晚上20時,大氣穩(wěn)定層結(jié)趨于穩(wěn)定之后,風(fēng)切變一直保持比較平穩(wěn),直至次日06時。通常06—10時和16—20時為大氣層結(jié)穩(wěn)定與不穩(wěn)定的過渡期,風(fēng)切變指數(shù)變化比較快速。從兩測風(fēng)塔位置來看,測風(fēng)塔1#緯度低且位于沿海地區(qū),受海陸熱力性質(zhì)差異影響較大,風(fēng)切變指數(shù)日變化不如測風(fēng)塔2#顯著。

    整體來看,去掉最低層風(fēng)速的擬合方案與其他幾個方案結(jié)果差異較大,因而推斷認(rèn)為近地層梯度測風(fēng)并非層數(shù)越多越好,距地面10 m高度以上一定間隔的4層左右的風(fēng)觀測比較適宜。

    2.4 推算方案和結(jié)果分析

    根據(jù)上述風(fēng)切變指數(shù)變化特征設(shè)置3種推算方案。方案1:采用測風(fēng)塔低層月平均風(fēng)速擬合月風(fēng)切變指數(shù),根據(jù)次高層小時風(fēng)速和對應(yīng)月風(fēng)切變指數(shù)推算最高層小時風(fēng)速。方案2:采用測風(fēng)塔低層平均小時風(fēng)速擬合小時風(fēng)切變指數(shù),根據(jù)次高層小時風(fēng)速和對應(yīng)時次風(fēng)切變指數(shù)推算最高層小時風(fēng)速。方案3:采用測風(fēng)塔低層年平均風(fēng)速擬合年風(fēng)切變指數(shù),根據(jù)次高層小時風(fēng)速和年風(fēng)切變指數(shù)推算最高層小時風(fēng)速。

    表3為測風(fēng)塔低層年平均風(fēng)速及年風(fēng)切變指數(shù)。由表3可以看出,年平均風(fēng)切變指數(shù)接近于方案2計算的小時風(fēng)切變指數(shù)平均值。對比2座測風(fēng)塔不同時間尺度的風(fēng)切變指數(shù)來看,2座測風(fēng)塔因地理位置及地形地貌等因素的影響,使得風(fēng)切變指數(shù)有較大的差異。本文選用風(fēng)切變指數(shù)差異較大的測風(fēng)塔進(jìn)行相關(guān)分析研究,可排除因風(fēng)切變指數(shù)大小造成結(jié)果的不確定性,進(jìn)一步佐證方法的適用性。

    表3 方案3測風(fēng)塔年平均風(fēng)速和風(fēng)切變指數(shù)Table 3 Annual wind speed and wind shear index in Scheme 3

    采用上述3種方案推算最高層月平均風(fēng)速和風(fēng)功率密度并與實(shí)測結(jié)果對比(圖5)。由圖5可以看出,測風(fēng)塔3種推算月值風(fēng)速接近,與實(shí)測對比,測風(fēng)塔1#各月推算與實(shí)測持平或者偏小0.0~0.2 m·s-1(0.0~3.4%),與實(shí)測風(fēng)功率密度差值為-15~1 W·m-2(-5.9%~0.7%)。測風(fēng)塔2#各月推算與實(shí)測持平或者偏大0.0~0.2 m·s-1(0.0~4.3%),與實(shí)測風(fēng)功率密度差值為-8~30.5 W·m-2(-3.7%~6.5%)。綜合風(fēng)速和風(fēng)功率密度對比結(jié)果,方案2推算接近實(shí)測,優(yōu)于其他方案。

    圖5 測風(fēng)塔1#、2#最高層推算和實(shí)測月風(fēng)速(a、b)和月風(fēng)功率密度(c、d)分布Fig.5 Calculation at the highest layer of Tower 1#and Tower 2#,and observations of monthly wind speed(a,b)and wind power(c,d)

    采用上述3種方案推算最高層小時風(fēng)速和風(fēng)功率密度(圖6)并與實(shí)測結(jié)果對比可以看出,3種推算方案推算出風(fēng)速小時值存在一定差異,測風(fēng)塔1#小時推算值與實(shí)測差值約為-0.2~0.2 m·s-1(-3.0%~3.5%),與實(shí)測風(fēng)功率密度差值約為-17.9~10.1 W·m-2(-7.3%~3.4%)。測風(fēng)塔2#小時推算值與實(shí)測差值約-0.2~0.3 m·s-1(-3.3%~4.4%),與實(shí)測風(fēng)功率密度差值約為-1.9~38.3 W·m-2(-0.5%~7.5%)。2座測風(fēng)塔均為采用方案2推算結(jié)果與實(shí)測差異最小,明顯優(yōu)于其他方案推算結(jié)果。在風(fēng)切變較小的午間,推算風(fēng)功率密度均高于實(shí)測。測風(fēng)塔1#推算結(jié)果低于實(shí)測。整體來看,測風(fēng)塔2#推算結(jié)果高于實(shí)測。對比2座測風(fēng)塔還可以看出,由于風(fēng)切變指數(shù)有日、月變化特征,且日變化相較于月變化更為顯著,所以方案2即平均小時風(fēng)速擬合風(fēng)切變推算高層風(fēng)速和風(fēng)功率密度更接近實(shí)測。

    圖6 測風(fēng)塔1#(a)、2#(b)最高層推算和實(shí)測小時風(fēng)速1#(c)和2#(d)和風(fēng)功率密度分布Fig.6 Calculation at the highest layer of Tower 1#and Tower 2#,and observations of hourly wind speed(a,b)and wind power(c,d)

    表4為實(shí)測與推算年平均風(fēng)速和風(fēng)功率密度對比。各推算方案的年平均結(jié)果差異較小,其中方案2推算結(jié)果與實(shí)測最接近,風(fēng)速差異在1.6%以內(nèi),風(fēng)功率密度差異在3.7%以內(nèi)。

    表4 測風(fēng)塔1#、2#實(shí)測和推算年平均風(fēng)速、風(fēng)功率對比Table 4 Comparison of annual wind speed and wind power at Tower 1#and Tower 2#

    3 結(jié)論與討論

    (1)遼寧省月平均風(fēng)切變指數(shù)呈波動變化,春季風(fēng)切變指數(shù)最小,夏季風(fēng)切變指數(shù)最大,風(fēng)切變指數(shù)與風(fēng)速呈顯著的負(fù)相關(guān)。遼寧省風(fēng)切變指數(shù)存在明顯的日變化特征,小時風(fēng)切變指數(shù)比月平均風(fēng)切變指數(shù)變化更顯著。一天之內(nèi),10—16時風(fēng)切變指數(shù)最小,夜間20時至翌日06時風(fēng)切變指數(shù)穩(wěn)定且偏大。06—10時和16—20時為風(fēng)切變指數(shù)急劇變化的過渡時段。

    (2)從探空背景風(fēng)可知,遼寧省冪指數(shù)擬合規(guī)律可以延伸至180 m高度。受觀測條件限制,在非復(fù)雜地形情況下,近地層10 m高度以上間隔一定高度設(shè)立4層左右的風(fēng)觀測,基本可以滿足近地層風(fēng)資源評估需求。

    (3)不同方案風(fēng)切變指數(shù)推算測風(fēng)塔風(fēng)速和風(fēng)功率密度并與實(shí)測對比,月值差異在6.5%以內(nèi),小時差異在7.5%以內(nèi),年值差異在3.7%以內(nèi)。整體來講,各推算方案均比較合理可靠。因?yàn)轱L(fēng)切變指數(shù)日變化比月變化顯著,在利用不同時間尺度風(fēng)切變指數(shù)推算高層風(fēng)況時,利用小時風(fēng)切變指數(shù)推算的方案優(yōu)于采用月和年風(fēng)切變指數(shù)方案。同時,風(fēng)切變?nèi)兆兓斤@著,采用小時風(fēng)切變指數(shù)推算方案越優(yōu)于其他風(fēng)切變推算方案。

    猜你喜歡
    風(fēng)塔測風(fēng)塔探空
    探空氣球?yàn)楦呖諝庀蟆鞍衙}”
    一種自安裝海上測風(fēng)塔的運(yùn)輸和安裝穩(wěn)性分析
    WRF氣象模式在高海拔復(fù)雜山地風(fēng)電場中的應(yīng)用研究
    測風(fēng)塔法在棄風(fēng)電量評估中的應(yīng)用
    臺風(fēng)“威馬遜”影響廣東期間近地層風(fēng)特性
    TK-2GPS人影火箭探空數(shù)據(jù)與L波段探空數(shù)據(jù)對比分析
    廣西風(fēng)能資源分布高分辨率數(shù)值模擬應(yīng)用研究
    廣東電力(2016年6期)2016-07-16 04:28:40
    談風(fēng)塔的施工技術(shù)要點(diǎn)
    河南省風(fēng)能數(shù)據(jù)庫系統(tǒng)建設(shè)
    淺談凈舉力對探空氣球升速及施放高度的影響
    河南科技(2014年14期)2014-02-27 14:12:08
    亚洲av成人不卡在线观看播放网| 亚洲av成人精品一区久久| 91国产中文字幕| 小说图片视频综合网站| 老司机在亚洲福利影院| 日韩中文字幕欧美一区二区| 免费观看人在逋| 亚洲av第一区精品v没综合| 狠狠狠狠99中文字幕| av免费在线观看网站| 亚洲精品一卡2卡三卡4卡5卡| 在线观看66精品国产| 老司机靠b影院| 成人av一区二区三区在线看| 99久久精品热视频| 99在线人妻在线中文字幕| 舔av片在线| 免费搜索国产男女视频| 国产真实乱freesex| 国产成人aa在线观看| 精品熟女少妇八av免费久了| 美女 人体艺术 gogo| 中文字幕精品亚洲无线码一区| 麻豆成人午夜福利视频| 中文字幕精品亚洲无线码一区| 色av中文字幕| 亚洲欧洲精品一区二区精品久久久| 国产99久久九九免费精品| 国产伦一二天堂av在线观看| 黄色成人免费大全| 精品电影一区二区在线| av国产免费在线观看| 十八禁网站免费在线| 免费看美女性在线毛片视频| 一本久久中文字幕| 性欧美人与动物交配| 国产精品一区二区免费欧美| 欧美不卡视频在线免费观看 | 精品久久久久久久末码| 精品久久久久久久久久免费视频| 一区二区三区高清视频在线| 淫秽高清视频在线观看| 成年女人毛片免费观看观看9| 国产午夜精品论理片| 欧美日本亚洲视频在线播放| 久久久久久久精品吃奶| 国产aⅴ精品一区二区三区波| 国产精品99久久99久久久不卡| 国产又色又爽无遮挡免费看| 这个男人来自地球电影免费观看| 人妻丰满熟妇av一区二区三区| 不卡一级毛片| 精品久久久久久久毛片微露脸| 午夜激情福利司机影院| 久久欧美精品欧美久久欧美| 国产亚洲精品久久久久5区| 18禁黄网站禁片午夜丰满| 中文字幕久久专区| 国内少妇人妻偷人精品xxx网站 | 麻豆成人午夜福利视频| 欧美乱码精品一区二区三区| 色噜噜av男人的天堂激情| 男女之事视频高清在线观看| 欧美在线黄色| 无限看片的www在线观看| 亚洲av中文字字幕乱码综合| 黑人巨大精品欧美一区二区mp4| 日本五十路高清| 少妇被粗大的猛进出69影院| videosex国产| 色精品久久人妻99蜜桃| 国产成人一区二区三区免费视频网站| 成年免费大片在线观看| 精品乱码久久久久久99久播| 中文字幕熟女人妻在线| 母亲3免费完整高清在线观看| 欧美日韩福利视频一区二区| 又爽又黄无遮挡网站| 成人18禁在线播放| 国产成人精品久久二区二区免费| 亚洲专区国产一区二区| avwww免费| 久久人妻福利社区极品人妻图片| 久久性视频一级片| 亚洲人成电影免费在线| 精品久久久久久成人av| 怎么达到女性高潮| 成人精品一区二区免费| 日韩欧美在线乱码| 中文亚洲av片在线观看爽| 每晚都被弄得嗷嗷叫到高潮| 巨乳人妻的诱惑在线观看| 99在线人妻在线中文字幕| 亚洲黑人精品在线| 免费在线观看亚洲国产| 精品久久久久久久久久久久久| 9191精品国产免费久久| 亚洲电影在线观看av| 国产一区在线观看成人免费| 欧美性猛交黑人性爽| 久久性视频一级片| 99热这里只有精品一区 | www.www免费av| 99久久精品国产亚洲精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲午夜精品一区,二区,三区| 两个人的视频大全免费| 亚洲成人中文字幕在线播放| 亚洲成人国产一区在线观看| 51午夜福利影视在线观看| 欧美日韩福利视频一区二区| 午夜成年电影在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美 国产精品| 他把我摸到了高潮在线观看| 俄罗斯特黄特色一大片| 后天国语完整版免费观看| 亚洲国产精品久久男人天堂| 国产免费男女视频| 99久久久亚洲精品蜜臀av| 国产精品一及| 成人午夜高清在线视频| 最新在线观看一区二区三区| 窝窝影院91人妻| 日本黄大片高清| 国产精品久久久人人做人人爽| 国产精品自产拍在线观看55亚洲| 99riav亚洲国产免费| 欧美一级毛片孕妇| 91在线观看av| 麻豆久久精品国产亚洲av| 国产精品一及| 国产精品久久视频播放| 丁香六月欧美| 亚洲美女视频黄频| 操出白浆在线播放| 国产精品99久久99久久久不卡| 男女床上黄色一级片免费看| 免费在线观看成人毛片| 一级黄色大片毛片| 舔av片在线| 国产精品久久电影中文字幕| 亚洲 国产 在线| 舔av片在线| 午夜精品在线福利| 欧美一级毛片孕妇| 啦啦啦韩国在线观看视频| 日韩欧美在线乱码| 无限看片的www在线观看| 日韩欧美 国产精品| 国产精品综合久久久久久久免费| av欧美777| 久久久精品欧美日韩精品| tocl精华| svipshipincom国产片| 无限看片的www在线观看| 男人舔奶头视频| 国产三级黄色录像| 欧美色视频一区免费| 久久热在线av| 国产黄片美女视频| 嫩草影视91久久| 99在线人妻在线中文字幕| 久久午夜亚洲精品久久| 久久久久精品国产欧美久久久| 国产精品99久久99久久久不卡| 少妇人妻一区二区三区视频| 久久久久九九精品影院| 欧美乱色亚洲激情| 亚洲av成人av| 成人三级做爰电影| 日韩精品青青久久久久久| 国产成人精品久久二区二区91| 色av中文字幕| 成人国产综合亚洲| 99热只有精品国产| 色av中文字幕| 亚洲精品在线观看二区| 成年人黄色毛片网站| 亚洲国产日韩欧美精品在线观看 | 身体一侧抽搐| 日韩av在线大香蕉| 亚洲av成人一区二区三| 国产私拍福利视频在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲一区中文字幕在线| 90打野战视频偷拍视频| 两个人看的免费小视频| 一区二区三区高清视频在线| 亚洲国产欧美人成| tocl精华| 欧美日韩福利视频一区二区| 欧美一级毛片孕妇| 欧美绝顶高潮抽搐喷水| 女人高潮潮喷娇喘18禁视频| 波多野结衣高清作品| 毛片女人毛片| 夜夜躁狠狠躁天天躁| 我的老师免费观看完整版| 视频区欧美日本亚洲| 亚洲人成伊人成综合网2020| 国产精品亚洲美女久久久| 搡老岳熟女国产| 亚洲精品久久国产高清桃花| 可以在线观看毛片的网站| 可以免费在线观看a视频的电影网站| 精品高清国产在线一区| 亚洲自拍偷在线| av视频在线观看入口| 99久久99久久久精品蜜桃| 久久精品人妻少妇| 中文字幕人妻丝袜一区二区| 久99久视频精品免费| 国产精品野战在线观看| 一进一出抽搐gif免费好疼| 免费高清视频大片| 成人特级黄色片久久久久久久| 在线免费观看的www视频| 精品欧美国产一区二区三| 成人一区二区视频在线观看| 91字幕亚洲| 亚洲av日韩精品久久久久久密| 真人一进一出gif抽搐免费| 久久久精品大字幕| 亚洲电影在线观看av| 99热只有精品国产| 两性夫妻黄色片| 国产探花在线观看一区二区| 亚洲 国产 在线| 黄色视频,在线免费观看| 观看免费一级毛片| 国产精品av视频在线免费观看| 国产麻豆成人av免费视频| 久热爱精品视频在线9| or卡值多少钱| 正在播放国产对白刺激| 身体一侧抽搐| 欧美日韩中文字幕国产精品一区二区三区| 国产成人一区二区三区免费视频网站| 窝窝影院91人妻| 嫩草影视91久久| 最近视频中文字幕2019在线8| 精品无人区乱码1区二区| 超碰成人久久| 最近视频中文字幕2019在线8| 久久人妻av系列| 亚洲国产精品sss在线观看| 99精品久久久久人妻精品| 婷婷精品国产亚洲av| 欧美性猛交╳xxx乱大交人| 日本一本二区三区精品| 巨乳人妻的诱惑在线观看| 女警被强在线播放| 麻豆久久精品国产亚洲av| 51午夜福利影视在线观看| 法律面前人人平等表现在哪些方面| 熟女电影av网| 日本精品一区二区三区蜜桃| 可以在线观看的亚洲视频| 亚洲午夜理论影院| av天堂在线播放| 亚洲av第一区精品v没综合| 欧美另类亚洲清纯唯美| 变态另类成人亚洲欧美熟女| 国产精品免费一区二区三区在线| 一本一本综合久久| 日日摸夜夜添夜夜添小说| 亚洲真实伦在线观看| 青草久久国产| 黑人操中国人逼视频| 色噜噜av男人的天堂激情| 男女视频在线观看网站免费 | 亚洲乱码一区二区免费版| 亚洲av第一区精品v没综合| 国产久久久一区二区三区| 夜夜看夜夜爽夜夜摸| 国产黄a三级三级三级人| 免费在线观看亚洲国产| 国内精品一区二区在线观看| 国产精品久久久久久精品电影| 18禁裸乳无遮挡免费网站照片| 91老司机精品| 精品少妇一区二区三区视频日本电影| 国产精品,欧美在线| 久久久久国内视频| 可以在线观看毛片的网站| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 欧美久久黑人一区二区| 久久精品人妻少妇| 久久久久久久精品吃奶| 丁香六月欧美| 美女午夜性视频免费| 免费在线观看成人毛片| 亚洲一区二区三区色噜噜| 香蕉丝袜av| 亚洲国产精品久久男人天堂| 小说图片视频综合网站| 国产成年人精品一区二区| 狂野欧美白嫩少妇大欣赏| 欧美性长视频在线观看| 精品国产美女av久久久久小说| 在线观看免费午夜福利视频| 日本熟妇午夜| 亚洲熟妇中文字幕五十中出| 老司机午夜十八禁免费视频| 国产伦一二天堂av在线观看| 在线观看一区二区三区| 国产区一区二久久| 777久久人妻少妇嫩草av网站| 亚洲片人在线观看| 两性夫妻黄色片| 99久久99久久久精品蜜桃| 亚洲国产看品久久| 久久久久九九精品影院| 久久久久久久久中文| 午夜福利欧美成人| 真人一进一出gif抽搐免费| 日韩欧美在线乱码| 婷婷六月久久综合丁香| 久久久久九九精品影院| 一级毛片高清免费大全| 美女免费视频网站| 18禁裸乳无遮挡免费网站照片| 男男h啪啪无遮挡| avwww免费| 别揉我奶头~嗯~啊~动态视频| 日韩国内少妇激情av| 亚洲男人天堂网一区| 欧美黄色淫秽网站| 亚洲av成人不卡在线观看播放网| 精品无人区乱码1区二区| 日本免费一区二区三区高清不卡| 操出白浆在线播放| 成人国产一区最新在线观看| 麻豆成人午夜福利视频| 老熟妇仑乱视频hdxx| 欧美成狂野欧美在线观看| 久久香蕉精品热| 中国美女看黄片| 婷婷精品国产亚洲av| 免费人成视频x8x8入口观看| 国产精品免费视频内射| 中亚洲国语对白在线视频| 免费无遮挡裸体视频| 女生性感内裤真人,穿戴方法视频| 亚洲第一欧美日韩一区二区三区| 亚洲av五月六月丁香网| 午夜免费观看网址| 黄片大片在线免费观看| 搡老熟女国产l中国老女人| 国产高清激情床上av| 国产成人av教育| 婷婷精品国产亚洲av| 琪琪午夜伦伦电影理论片6080| 亚洲精品国产一区二区精华液| 成人av一区二区三区在线看| 一个人观看的视频www高清免费观看 | 最近视频中文字幕2019在线8| 精品久久久久久久久久久久久| av中文乱码字幕在线| 欧美成狂野欧美在线观看| 国内精品一区二区在线观看| 丰满人妻一区二区三区视频av | 一二三四在线观看免费中文在| 国产三级黄色录像| 亚洲真实伦在线观看| 99热6这里只有精品| 男人舔奶头视频| 非洲黑人性xxxx精品又粗又长| 又粗又爽又猛毛片免费看| 香蕉国产在线看| 国产精品国产高清国产av| 国产精品爽爽va在线观看网站| 亚洲精品国产一区二区精华液| 999久久久国产精品视频| 夜夜爽天天搞| 国产亚洲精品av在线| 国产成年人精品一区二区| 精品久久蜜臀av无| ponron亚洲| 久久久久国内视频| 女人爽到高潮嗷嗷叫在线视频| 天堂√8在线中文| 中文在线观看免费www的网站 | 级片在线观看| 天堂动漫精品| 一区二区三区国产精品乱码| 欧美性猛交╳xxx乱大交人| 成人18禁高潮啪啪吃奶动态图| 日韩成人在线观看一区二区三区| 老司机福利观看| 人人妻人人澡欧美一区二区| 久9热在线精品视频| 国产亚洲精品久久久久久毛片| 亚洲五月婷婷丁香| ponron亚洲| 精品无人区乱码1区二区| 国产男靠女视频免费网站| 欧美3d第一页| 国产精品av视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 此物有八面人人有两片| 十八禁人妻一区二区| 日韩精品中文字幕看吧| www.999成人在线观看| 免费看日本二区| 午夜福利高清视频| 又大又爽又粗| 成人高潮视频无遮挡免费网站| 麻豆国产97在线/欧美 | 国产爱豆传媒在线观看 | 国产一区在线观看成人免费| 高清在线国产一区| 日本 av在线| 久久久久国产一级毛片高清牌| 婷婷精品国产亚洲av在线| 两性夫妻黄色片| 老熟妇乱子伦视频在线观看| 久久久国产精品麻豆| 成人18禁在线播放| 日韩欧美精品v在线| 亚洲人成电影免费在线| 黑人巨大精品欧美一区二区mp4| 51午夜福利影视在线观看| 亚洲av片天天在线观看| 国产高清视频在线观看网站| www.精华液| 一级a爱片免费观看的视频| aaaaa片日本免费| 国产野战对白在线观看| 国产成人啪精品午夜网站| 免费av毛片视频| 国内精品久久久久精免费| 亚洲av日韩精品久久久久久密| 国产午夜精品论理片| 99国产精品一区二区蜜桃av| 亚洲电影在线观看av| 操出白浆在线播放| 在线视频色国产色| 久久久久精品国产欧美久久久| 亚洲片人在线观看| 性欧美人与动物交配| 亚洲av成人不卡在线观看播放网| 91国产中文字幕| 亚洲一区中文字幕在线| 国产精品av久久久久免费| 欧美一级a爱片免费观看看 | 国产97色在线日韩免费| 在线观看66精品国产| 欧美成人性av电影在线观看| 麻豆一二三区av精品| 午夜福利18| 精品欧美国产一区二区三| 国产亚洲av高清不卡| 制服诱惑二区| 久久久久亚洲av毛片大全| 欧美一级毛片孕妇| 国产精品,欧美在线| 黄色毛片三级朝国网站| 亚洲av第一区精品v没综合| 一级黄色大片毛片| 一级作爱视频免费观看| 国产高清videossex| 深夜精品福利| 欧美日本视频| 精品久久久久久久末码| 久久久久久人人人人人| 久久草成人影院| 国产精品久久久久久精品电影| 精品国内亚洲2022精品成人| 午夜亚洲福利在线播放| 国产亚洲精品久久久久久毛片| 久99久视频精品免费| 日韩成人在线观看一区二区三区| 中文字幕最新亚洲高清| 欧美日韩中文字幕国产精品一区二区三区| 少妇熟女aⅴ在线视频| 国产精品久久视频播放| 日韩免费av在线播放| 成年女人毛片免费观看观看9| 精品国产美女av久久久久小说| 国产激情欧美一区二区| 国产高清有码在线观看视频 | 欧美成人性av电影在线观看| 特级一级黄色大片| 国产亚洲精品久久久久5区| 国产伦在线观看视频一区| 亚洲狠狠婷婷综合久久图片| 免费一级毛片在线播放高清视频| 亚洲18禁久久av| 国产精品一区二区免费欧美| 国产精品一及| 久久久久精品国产欧美久久久| www日本黄色视频网| 久久久国产成人精品二区| 免费看a级黄色片| 99久久精品国产亚洲精品| 国产精品久久久久久亚洲av鲁大| 在线看三级毛片| 久久这里只有精品中国| 国产激情偷乱视频一区二区| cao死你这个sao货| 少妇粗大呻吟视频| 国产v大片淫在线免费观看| 久久婷婷人人爽人人干人人爱| 国产人伦9x9x在线观看| 99热只有精品国产| 淫妇啪啪啪对白视频| www.999成人在线观看| 中文字幕最新亚洲高清| 深夜精品福利| 又黄又粗又硬又大视频| 久久性视频一级片| 精品久久久久久,| 久久国产精品人妻蜜桃| 久久中文字幕一级| 亚洲 欧美 日韩 在线 免费| 夜夜夜夜夜久久久久| 精品国产美女av久久久久小说| 精品久久久久久久人妻蜜臀av| 亚洲狠狠婷婷综合久久图片| 一本综合久久免费| 搡老妇女老女人老熟妇| 又黄又爽又免费观看的视频| 亚洲一码二码三码区别大吗| 亚洲va日本ⅴa欧美va伊人久久| 岛国在线免费视频观看| 白带黄色成豆腐渣| 国产精品久久视频播放| 国产乱人伦免费视频| 久久精品91无色码中文字幕| 三级男女做爰猛烈吃奶摸视频| 丝袜人妻中文字幕| 国产av在哪里看| 午夜久久久久精精品| 国产精品自产拍在线观看55亚洲| 国产视频一区二区在线看| 日本撒尿小便嘘嘘汇集6| 午夜免费观看网址| 亚洲一区二区三区不卡视频| 久久九九热精品免费| 美女大奶头视频| 人妻丰满熟妇av一区二区三区| 18禁美女被吸乳视频| 我的老师免费观看完整版| 欧美成人免费av一区二区三区| 哪里可以看免费的av片| 亚洲色图av天堂| 日本一区二区免费在线视频| 97人妻精品一区二区三区麻豆| www.精华液| 国产私拍福利视频在线观看| xxx96com| 国产午夜精品久久久久久| 国产av麻豆久久久久久久| 欧美性猛交黑人性爽| 香蕉久久夜色| 久久久久国内视频| 亚洲美女黄片视频| 高潮久久久久久久久久久不卡| av视频在线观看入口| ponron亚洲| 丰满人妻熟妇乱又伦精品不卡| 亚洲人成77777在线视频| 国产视频一区二区在线看| 欧美丝袜亚洲另类 | 十八禁人妻一区二区| 亚洲乱码一区二区免费版| 50天的宝宝边吃奶边哭怎么回事| 成人午夜高清在线视频| 又粗又爽又猛毛片免费看| 级片在线观看| 国产黄色小视频在线观看| 少妇粗大呻吟视频| 成人18禁在线播放| 国内毛片毛片毛片毛片毛片| 色哟哟哟哟哟哟| 精品久久久久久成人av| or卡值多少钱| 久久久久久免费高清国产稀缺| 特大巨黑吊av在线直播| 亚洲av日韩精品久久久久久密| 国产亚洲av高清不卡| 狂野欧美激情性xxxx| 国产一区二区在线av高清观看| 久热爱精品视频在线9| 亚洲人成电影免费在线| 色在线成人网| www.熟女人妻精品国产| 99在线人妻在线中文字幕| 91九色精品人成在线观看| av视频在线观看入口| 亚洲中文字幕日韩| 日本一区二区免费在线视频| av中文乱码字幕在线| 欧美另类亚洲清纯唯美| 久久99热这里只有精品18| 国内精品久久久久精免费| 亚洲人成77777在线视频| 91九色精品人成在线观看| 国产高清有码在线观看视频 | 一区福利在线观看| 日韩欧美在线二视频| 日本精品一区二区三区蜜桃| 看黄色毛片网站| 亚洲av第一区精品v没综合| 在线观看免费视频日本深夜| 给我免费播放毛片高清在线观看| 亚洲全国av大片| 色综合欧美亚洲国产小说| 精品乱码久久久久久99久播| 国产av麻豆久久久久久久| 极品教师在线免费播放| 国产精品1区2区在线观看.|