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

    海氣邊界層內(nèi)風電場尾流及功率特性的中尺度數(shù)值模擬研究

    2024-01-05 09:24:14樂可定郁冶王異成吳春雷王強羅坤樊建人
    海洋預報 2023年6期
    關鍵詞:尾流風力機海面

    樂可定,郁冶,王異成,吳春雷,王強,羅坤*,樊建人

    (1.浙江浙能國電投嵊泗海上風力發(fā)電有限公司,浙江舟山 202450;2.國家電投集團浙江新能源有限公司,浙江杭州 310016;3.杭州意能電力技術有限公司,浙江杭州 310027;4.浙江大學能源清潔利用國家重點實驗室,浙江杭州 310027)

    0 引言

    在能源領域,風能具有巨大的開發(fā)潛力,它可以提供大量的清潔與可再生能源,是中國實現(xiàn)碳達峰和碳中和的重要手段[1]。截止2019 年底,全球風力發(fā)電量達4 057 億千瓦時,占總發(fā)電量的5.5%。與陸上風資源相比,海上風資源的風能密度普遍較高,平均風速更大且風能穩(wěn)定性更好,是未來風能發(fā)展的重要方向。雖然海上相較于陸上受地形影響小,但是海氣運動的復雜性(例如海氣邊界的自然通量交換、動態(tài)海面粗糙度等因素)使得科學研究具有更大的挑戰(zhàn)性。為了更好地發(fā)展和規(guī)劃海上風電產(chǎn)業(yè),亟需開展針對海上風電場的運行特性研究。

    海上風電場位于海洋大氣邊界層(Marine Atmospheric Boundary Layer,MABL)內(nèi)[2],大氣邊界層存在大氣和海洋表面之間的多物理、多尺度耦合作用過程,其可以通過改變海洋近表面多種物理通量顯著影響大氣邊界層的結構[3],因此,在近海大氣運動的研究中,海洋與大氣的相互作用被認為是必不可少的考量因素。此外,海浪的存在會影響邊界層湍流分布以及近海面的垂直風速剖面,因此,研究海上大氣運動還需對海面粗糙度(z0)進行評估。CHARNOCK[4]考慮了海面風速對z0的影響。另外,海面粗糙度還取決于表面風速和波峰傳播速度的關系[5],一般來說,波浪越高,海洋表面粗糙度越長,z0與波浪高度和陡度的依賴關系可以用無量綱法表示[6]。總的來說,海浪的影響可通過物理參數(shù)化的方式近似表征。例如,GOLBAZI 等[7]應用CHARNOCK 參數(shù)化方法估算了美國東北海岸海上風能評估的z0,并揭示了其對輪轂高度速度預測的重要性。BAHAMONDE 等[8]進一步結合大氣穩(wěn)定性來計算發(fā)電量,通過與現(xiàn)場數(shù)據(jù)進行比較,得到了較好的評估結果。研究證明,獨立的大氣模型不足以模擬海洋大氣邊界層的大氣運動,需要進一步耦合海洋模型來考慮海面和大氣之間的相互作用。

    除了海氣環(huán)境本身的復雜性以外,風電場中風機的分布也會改變海氣邊界層中的大氣動量和湍流通量。觀測研究表明,在較穩(wěn)定的大氣邊界層內(nèi),風電場的尾跡可以延伸幾十千米[9],且下風向近下墊面的溫度也會發(fā)生相應變化[10-11]。與海浪參數(shù)化方式類似,風電場的影響也可以表示為空氣動力學粗糙度的增加。通過采用定性增加z0的方法,研究模擬得到了風電場區(qū)域存在1~2 K 的溫度變化[12-13]。WANG 等[14]采用大氣-海洋耦合模型,通過修改表面粗糙度和位勢高度來表示風電場,發(fā)現(xiàn)陸上風電場的溫度升高1 ℃,海上風電場的溫度降低1 ℃。這種方法可以有效地估計具有數(shù)千臺風力機的大型風電場甚至超大型風場群的影響,但是,z0的增量是根據(jù)經(jīng)驗確定的,這種方法無法捕捉垂直風廓線的變化以及風力機群之間的相互干擾。

    近幾年,將風力機表示為平均來流上的動量匯的參數(shù)化方法被廣泛用于模擬和量化風電場對大氣的作用。風電場參數(shù)化(Wind Farm Parameterization,WFP)由BLAHAK 等[15]和FITCH 等[16]開發(fā)和發(fā) 展,并 由JIMéNEZ 等[17]、ARCHER 等[18]和WU等[19]等先后進行改進。WFP 是研究風電場大氣影響的有效工具,利用適量的計算資源便可以得到垂直方向上風電場的氣象影響以及風電場群水平方向的相互干涉。WANG 等[20]將WFP 與天氣研究與預測(Weather Research and Forecast,WRF)模型相結合,研究了中國新疆兩個風電場之間的尾流干擾以及功率虧損。PLATIS 等[21]采用相同方法衡量了風電場后方的尾跡長度,發(fā)現(xiàn)此模型高估了尾跡的長度。然而,目前風電場中尺度數(shù)值研究通常將海面和風電場的影響解耦,使得結果通常具有一定的誤差與局限性。風電場參數(shù)化的模擬通常忽略了大氣和海洋的相互作用,而針對海上風資源的研究則忽略了風電場的影響。例如,WU 等[22]發(fā)現(xiàn),海氣邊界條件可能會改變海上風能密度(約5%),但相鄰風電場的影響卻未評估,因此會高估海上風能密度[23]。

    因此,為了研究海氣邊界條件下的風電場運行特性,我們在真實沿海環(huán)境下模擬了兩個風電場,分別為設計的虛擬風電場和真實的風電場?;谥谐叨葦?shù)值模型并結合海上多物理建模方法,評估了海上風電場的運行性能。通過開展多組對比模擬實驗及分析,分別研究了兩個風電場下的尾流特性以及功率輸出特性。通過比較有無海面的模擬結果和兩個風電場的不同表現(xiàn),進而得到海氣作用對風電場的影響以及不同風力機排布下海氣作用的變化。

    1 研究方法

    1.1 研究對象

    為了展現(xiàn)較符合實際環(huán)境的研究結果,選取已規(guī)劃并將持續(xù)發(fā)展海上風電的杭州灣地區(qū)作為研究區(qū)域,對兩個風電場進行模擬與評估,風電場位置見圖1。1號風電場為假設的理想風電場,風力機排布方式均勻,即風力機具有相同的水平間距。風力機放置在WRF 水平網(wǎng)格點內(nèi),由于WRF 的水平網(wǎng)格同經(jīng)緯度平行,因此風電場整體與經(jīng)緯度平行,目的是為了更加清晰地觀察到不同下墊面下風力機的尾流特征和功率分布特征。具體地,1 號風電場的中心地理坐標設定為30.46°N,121.61°E,面積為10 km×10 km,100 臺風力機均勻分布其中,即徑向與展向距離均為1 km;風力機根據(jù)當前典型的海上風力發(fā)電機的推力和功率系數(shù)進行參數(shù)化,額定功率為4.0 MW,輪轂高度為99 m,葉片直徑為127 m;切入和切出風速分別為3 m/s和30 m/s,低于或高于這兩個數(shù)值時,風力機不運行,靜態(tài)推力系數(shù)為0.13。2 號風電場為某實際運行的風電場,中心地理坐標為30.68°N,121.64°E,其具體位置以及風力機排布位置可見圖1、圖2,運行風力機額定功率為6.25 MW,輪轂高度為110 m,葉片直徑為172 m,切入和切出風速分別為3 m/s 和25 m/s。

    圖1 研究區(qū)域以及1號風電場和2號風電場位置Fig.1 Study area and the locations of wind farm#1 and#2

    圖2 2號風電場的尾流特性Fig.2 The wake flow of wind farm#2

    本研究所使用的模型為上層海洋-大氣-風電場耦合模型。通過耦合海面-大氣模型實現(xiàn)海面動量、熱量、水汽量等物理通量的相互交換;通過耦合海面粗糙度模型表征海面波浪的動態(tài)變化;通過耦合風電場參數(shù)化模型表征風力運行對大氣的作用。各模型將在以下章節(jié)詳細介紹。

    1.2 風場參數(shù)化模型

    為了表征風電場對MABL 的作用,采用了FITCH 等[24]提出的WFP 模型。該模型將風電場內(nèi)每臺風力機視為致動盤,其對MABL 的作用表征為大氣動量的匯以及湍流的源,其中,動量變化項和湍動能變化項作為交互變量。方程如下:

    式中:Nt為網(wǎng)格內(nèi)風力機的數(shù)量;CT為推力系數(shù);CTKE=CT-CP,其中CP為功率系數(shù);v、u代表水平速度分量;|V|為速度標量;A為風力機風輪旋轉(zhuǎn)面積;Δz為垂直層的距離。CT、CP和A由風力機設計參數(shù)決定。

    1.3 海氣邊界模型

    本研究采用大氣模式WRF 和3DPWP(3D Price-Weller-Pinkel)上層海洋環(huán)流模式[3]構建海洋大氣環(huán)境。已有研究表明,該模型顯著提高了近海氣象過程模擬的準確性,且與WRF耦合的可靠性得到驗證[25]。此模型基于靜水力學模型,通過求解動量、熱量和鹽度平衡方程,評估上層海洋對地表附近流動的響應。對應的方程如下:

    式中:T和S分別為溫度和鹽度;P是靜水壓力;H、E和τ分別是熱量、鹽度和動量通量;V為水平流速;W是速度的垂直分量;f是科里奧利參數(shù)。

    海洋模型設置:水平方向上,海洋模型和大氣模型(見1.5 節(jié))具有相同的水平方向網(wǎng)格;垂直方向上,從海面延伸至海下400 m 共應用了30 層垂直層,其中,頂部劃分為20 層,垂直間距為10 m,下層劃分為10 層,分辨率為20 m。根據(jù)LEE 等[3]的建議,同時采用了較高分辨率的先進微波掃描輻射計(Advanced Microwave Scanning Radiometer - 2,AMSR-2)海洋數(shù)據(jù)集[26]為海洋模型提供初始邊界條件。

    1.4 海面粗糙度模型

    海面粗糙度因海浪的存在而動態(tài)變化,z0間接決定了海面對近海面的垂直風廓線分布。研究發(fā)現(xiàn)3DPWP 海洋模型對海面粗糙度長度參數(shù)化方案較為敏感[27],因此,海面粗糙度模型的選擇對模型的準確性至關重要。GREESHMA等[27]研究了3DPWP對粗糙度參數(shù)化方案的敏感性,通過與實測的對比,推薦了最新的Donelan 方案,它比WRF 中默認的Garratt 方案能提供更好的模擬結果,因此,本研究也采用Donelan方案。診斷方程為:

    式中:u*為摩擦速度;g為重力加速度。

    基于海氣邊界層模型3DPWP 與海面粗糙度模型Donelan,可以構建出適用于中尺度WRF 模式下的大氣海洋環(huán)境,從而為MABL 條件下風電場尾流及功率輸出特性模擬提供良好的海氣環(huán)境。

    1.5 大氣模型

    大氣模型采用WRF(Version 4.3)模式。針對兩個風電場研究區(qū)域,模型采用如下設置:采用雙向嵌套的方式,共設置3層計算域,最外兩層設計150×150個水平網(wǎng)格單元,最內(nèi)層為250×250,其水平分辨率為500 m,父域為子域網(wǎng)格大小的3 倍,即由外到內(nèi)網(wǎng)格的分辨率分別為4.5 km,1.5 km 和0.5 km,最內(nèi)層網(wǎng)格中心設置在風電場中心。垂直方向網(wǎng)格劃分為71 層,400 m 以下有25 個垂直層,其中15個垂直層用于劃分風力機風輪旋轉(zhuǎn)面。

    WRF 模式采用美國國家環(huán)境預報中心(National Centers for Environmental Prediction,NCEP)提供的FNL 數(shù)據(jù)集(Final Reanalysis Data)為模型提供初始邊界條件,分辨率為6 h,模擬日期為10 月23 日00 時(世界時,下同)—10 月31 日00時。最內(nèi)層時間步長設置為3 s。大氣邊界層參數(shù)化方案選擇MMYN2(Mellor-Yamada-Nakanishi-Niino Level-2.5),微物理方案采用WSM-6,長波和短波輻射過程分別采用RRTM 模型(Rapid Radiative Transfer Model)和Goddard 模型,其他物理模型多采用默認模型設置。

    至此,上層海洋-大氣-風電場耦合的多物理模型框架搭建完成。為了分析和發(fā)現(xiàn)海洋環(huán)境下的風資源和風電場運行特性,本研究設計了對照實驗,對照組(CTRL)為平坦表面,實驗組(SS)為具有波浪的海面。平坦表面不采用海洋模型,設置表面粗糙度為z0=0.000 1(即WRF 模型中海面粗糙度的初始默認z0值)。為控制變量,CTRL與SS僅存在海面粗糙度不同和是否采用海洋模型的差別。模型的設置見表1。

    表1 模型設置Tab.1 Model configuration

    2 結果與討論

    不同風向以及熱力條件均會對大氣的湍流特性產(chǎn)生影響。目前,已有的研究表明風電場僅在穩(wěn)定大氣邊界層下具有較為明顯的尾流效應,不穩(wěn)定邊界層下風電場尾流與周圍大氣混合劇烈,風電場造成的大氣影響較為混亂,多種非線性影響較難分析,因此多數(shù)研究通常在穩(wěn)定或中性大氣邊界層下分析風電場的中尺度大氣作用。同樣,本研究選取了模擬時間段內(nèi)主風向下連續(xù)4 h 穩(wěn)定大氣邊界層下的結果,并通過取平均的方法研究了海上風電場的尾流和功率輸出。

    2.1 風電場尾流特性

    對于均勻布局的1號風電場(見圖3),隨著上游越來越多的風力機消耗來流空氣動能轉(zhuǎn)化為電能,風速虧損(Wind Speed Deficit,WSD)沿風向逐漸增大,最終在風場后方形成數(shù)十公里的尾流,隨后,通過從周圍空氣中補充動能,尾流不斷減弱至最終消失。通過平坦表面和海面的對比可以發(fā)現(xiàn),在平坦下墊面上,風電場區(qū)域內(nèi)風力機的尾跡相對獨立,每排風力機尾跡較為清晰,即使在風電場數(shù)十公里后依舊保留部分尾跡分支;然而,在MABL 中,上游風力機后方的尾跡具有明顯的展向擴散,導致尾流之間易于混合,最大WSD 可達1.5 m/s,同時,尾流區(qū)與周圍空氣的動量交換更劇烈,加快了尾流的恢復,最終導致尾跡延伸距離變短,但寬度有所增加。PLATIS 等[9]也發(fā)現(xiàn)了類似結果,F(xiàn)ITCH 模型低估了尾跡的恢復率,高估了尾跡的長度。因此,忽視海洋對大氣的反饋作用會高估風電場的尾流長度,低估尾流的展向擴展。

    進一步地,通過關閉海洋模型,并且只設置一個粗糙度可變的表面,尾跡與對照組幾乎沒有差異(尾流長度僅略減?。Q言之,粗糙度的變化在一定程度上有助于尾跡區(qū)風速的恢復,但不是引起尾流混合和加快耗散的主要原因,而海氣相互作用是導致此現(xiàn)象的主要因素。

    對于實際的海上風電場,風力機布局相對復雜,造成的尾流效應也相對復雜。位于左右兩側(cè)的風力機群下游形成兩條尾流(見圖2),對比平坦表面以及海氣環(huán)境下的尾流特性,差距較為明顯。在平坦表面上,尾流效應呈現(xiàn)出較強的不穩(wěn)定性以及非局地性,在下游30 km 外依舊能產(chǎn)生大于0.5 m/s的速度虧損,WSD連續(xù)性較差,尾流長度約為45 km。相反,在海氣環(huán)境下的風電場有著連續(xù)的尾流擴散與衰減,下風電場區(qū)域形成較強尾流區(qū),隨著距離的增加,尾流逐漸耗散,尾流長度約為35 km。對比1號和2號海上風電場,真實風電場下的尾流發(fā)展平穩(wěn),而理想排布下的風電場尾流存在絮狀發(fā)展的情況。

    通過兩種不同邊界層的對比來看,實際風電場與均勻布局風電場都有著一致的結果,即海洋邊界層下的尾流具有更強的展向擴散,在近風電場區(qū)域能形成寬度更大、強度更強的速度虧損區(qū),最終形成的尾流具有較短的耗散長度。除此之外,從圖2中還可以發(fā)現(xiàn),平坦表面上的實際風電場由于產(chǎn)生更多湍流導致了尾流的不連續(xù)性以及非局地性,然而在真實的海洋環(huán)境下,由于海面吸收了一部分湍流,會大大增強尾流的穩(wěn)定性與延展性[28]。由此我們得出,實際海上風電場的模擬對海氣環(huán)境有著更強的依賴性,海氣模型的選用尤為關鍵。

    2.2 風電場功率輸出特性

    由于尾流效應,下風向的風力機的性能將受到上游風力機的影響,其發(fā)電功率會相應降低,圖4和圖5 中展示了與圖3 和圖2 相對應的風力機功率輸出分布。在不同下墊面下,1 號均勻分布電場的風力機具有不同的功率輸出空間分布(見圖4)。平坦表面上,下游風力機的平均功率輸出為整場功率的96%左右,即約4%的功率受到上游風力機的抑制,隨著尾流干涉沿著風速方向不斷累加,最下游風力機的功率降低約10%。從圖4b 可以看出,海面條件下的尾流對功率輸出的影響更為明顯,下游風力機的功率輸出較正常平均功率下降約10%,最下游風力機的功率降低約15%,即比在平坦表面上進一步下降了5%。

    圖4 1號風電場的功率輸出Fig.4 The power output of wind farm#1

    圖5 2號風電場功率輸出Fig.5 The power output of wind farm#2

    在實際風電場模擬中,尾流對功率輸出的影響也較為明顯,下游風力機風功率低于2.2 MW,上游風力機在2.6 MW 以上(見圖5),差距約為15%,結果與均勻排布的風力機平均功率輸出結果相類似。此風向下,海面吸收了部分大氣動能,降低了風力機的功率輸出。MABL 下的風電場依舊顯現(xiàn)出較強的尾流干涉,但由于真實風電場的風力機展向間距布置較大,因此相較于均勻排布,耦合多物理模型與CTRL模型在評估尾流效應對功率影響時的差距不大。

    總的來說,海上風電場尾流效應對整場功率輸出具有較大的影響,在海面的作用下海上風電場的尾流更易展向擴張,導致下游風力機受到的尾流干擾增強,因此,在進行海上風力機布局時應重視風力機展向間距。另一方面,在海上風電場的數(shù)值模擬中,評估海氣反饋對風電場功率影響尤為重要,使用耦合海洋-大氣-風電場的多物理模型是一種有效的方法。

    3 結論

    本文基于中尺度大氣模型WRF,結合上層海洋反饋和風電場參數(shù)化的數(shù)值方法,研究了海上風電場在海面上的風資源與運行特性。結果表明,使用耦合海洋-大氣-風電場模型尤為關鍵。海氣反饋主要會造成以下影響:

    ①尾流展向擴散增強,風力機尾流具有易混合特性。

    ②減弱風力機的湍流效應,使得尾流延伸的穩(wěn)定性增強,風電場非局地性影響減弱。

    ③下游風力機功率輸出受上游尾流的疊加影響,功率衰減可達15%,相較于平坦表面,海上風力機的尾流干涉更強。

    以上現(xiàn)象的產(chǎn)生主要來源于耦合海洋模型導致的海氣通量變化,這種海氣相互作用需在海上風電中尺度模擬中得到足夠重視,同時海浪造成的風切變變化對于功率評估也是不可忽略的因素。本研究充分發(fā)掘了海洋-大氣-風電場之間的相互作用對風電場運行特性的影響。此外,風力機的布局在本研究中也為影響因素之一,但并未深入研究。后續(xù)研究需針對風力機布局和裝機容量等風力機的物理特性對海上風電場的運行影響進一步深入挖掘,為海上風電場的設計提供更多的理論指導。

    猜你喜歡
    尾流風力機海面
    海面床,輕輕搖
    第六章 邂逅“胖胖號”
    小學科學(2020年6期)2020-06-22 13:13:41
    基于UIOs的風力機傳動系統(tǒng)多故障診斷
    測控技術(2018年5期)2018-12-09 09:04:38
    海面上的“一千座埃菲爾鐵塔”
    飛機尾流的散射特性與探測技術綜述
    雷達學報(2017年6期)2017-03-26 07:53:06
    錐形流量計尾流流場分析
    大型風力機整機氣動彈性響應計算
    小型風力機葉片快速建模方法
    太陽能(2015年6期)2015-02-28 17:09:35
    水面艦船風尾流效應減弱的模擬研究
    九九热线精品视视频播放| 日韩一本色道免费dvd| 亚洲欧美精品综合久久99| 国产精品嫩草影院av在线观看 | 一夜夜www| 中文字幕熟女人妻在线| 偷拍熟女少妇极品色| 久久久久久久午夜电影| 免费高清视频大片| 久久精品久久久久久噜噜老黄 | 热99在线观看视频| 国产精品福利在线免费观看| 国产三级中文精品| 亚洲专区中文字幕在线| 午夜福利高清视频| 一级黄色大片毛片| 国内精品一区二区在线观看| 欧美日本视频| 观看免费一级毛片| 99久国产av精品| 亚洲精品456在线播放app | 在线观看av片永久免费下载| 老熟妇仑乱视频hdxx| 国产午夜精品论理片| 蜜桃久久精品国产亚洲av| 一区二区三区免费毛片| 亚洲熟妇熟女久久| 精品人妻偷拍中文字幕| 国产亚洲精品av在线| 亚洲成a人片在线一区二区| 嫩草影院新地址| 婷婷色综合大香蕉| 欧美不卡视频在线免费观看| 亚洲图色成人| 亚洲欧美日韩高清专用| 国产精品三级大全| 日日啪夜夜撸| 淫妇啪啪啪对白视频| 国产淫片久久久久久久久| 淫秽高清视频在线观看| 亚洲自拍偷在线| 婷婷亚洲欧美| 亚洲天堂国产精品一区在线| 国产精品嫩草影院av在线观看 | 露出奶头的视频| 国产成人aa在线观看| 久久精品影院6| 91午夜精品亚洲一区二区三区 | 日韩精品有码人妻一区| 又黄又爽又免费观看的视频| 欧洲精品卡2卡3卡4卡5卡区| 在线观看午夜福利视频| 精品午夜福利在线看| 神马国产精品三级电影在线观看| 熟妇人妻久久中文字幕3abv| 久久久久精品国产欧美久久久| 国产人妻一区二区三区在| 狠狠狠狠99中文字幕| 国产精品一区二区三区四区久久| 国产三级中文精品| 两性午夜刺激爽爽歪歪视频在线观看| 高清在线国产一区| 午夜爱爱视频在线播放| 欧美黑人巨大hd| 国产精品一区二区三区四区久久| 搡老岳熟女国产| 日本黄色视频三级网站网址| 内射极品少妇av片p| 国产成人a区在线观看| 最新在线观看一区二区三区| 高清毛片免费观看视频网站| 九九热线精品视视频播放| 日本 av在线| 国产女主播在线喷水免费视频网站 | 给我免费播放毛片高清在线观看| 美女xxoo啪啪120秒动态图| 亚洲国产欧美人成| 搡女人真爽免费视频火全软件 | 少妇人妻精品综合一区二区 | 麻豆成人午夜福利视频| 偷拍熟女少妇极品色| 精品午夜福利在线看| 亚洲精品日韩av片在线观看| 亚洲国产高清在线一区二区三| 伦精品一区二区三区| 美女被艹到高潮喷水动态| 国产精品永久免费网站| 成年女人毛片免费观看观看9| 久久国内精品自在自线图片| 岛国在线免费视频观看| 亚洲精品粉嫩美女一区| 搡女人真爽免费视频火全软件 | 日韩大尺度精品在线看网址| 欧美中文日本在线观看视频| 最后的刺客免费高清国语| 国产一区二区激情短视频| 亚洲国产精品sss在线观看| 草草在线视频免费看| 香蕉av资源在线| 国产高清三级在线| 蜜桃久久精品国产亚洲av| 精品99又大又爽又粗少妇毛片 | 国产成人a区在线观看| 国产精品一区二区三区四区免费观看 | 99久久精品国产国产毛片| 女的被弄到高潮叫床怎么办 | 国内少妇人妻偷人精品xxx网站| 国产高潮美女av| 国产午夜精品久久久久久一区二区三区 | 老司机午夜福利在线观看视频| 性欧美人与动物交配| 淫妇啪啪啪对白视频| 欧美绝顶高潮抽搐喷水| 日本黄色视频三级网站网址| 搡女人真爽免费视频火全软件 | 深夜a级毛片| 美女高潮喷水抽搐中文字幕| 亚洲av熟女| 午夜爱爱视频在线播放| 国产亚洲欧美98| 精品久久久噜噜| 少妇熟女aⅴ在线视频| 精品不卡国产一区二区三区| 亚洲国产精品合色在线| 97碰自拍视频| 深夜a级毛片| 哪里可以看免费的av片| 深爱激情五月婷婷| 欧美成人一区二区免费高清观看| x7x7x7水蜜桃| 国产aⅴ精品一区二区三区波| 欧美激情国产日韩精品一区| 国产免费一级a男人的天堂| 1000部很黄的大片| 国产三级中文精品| 搡老岳熟女国产| 久久这里只有精品中国| av在线天堂中文字幕| 舔av片在线| av视频在线观看入口| 久久久久九九精品影院| 成年人黄色毛片网站| 美女 人体艺术 gogo| 欧美日韩亚洲国产一区二区在线观看| 成熟少妇高潮喷水视频| 精品久久久噜噜| 黄色配什么色好看| 国产三级在线视频| 美女黄网站色视频| 亚洲国产精品成人综合色| 久久精品国产99精品国产亚洲性色| 人妻丰满熟妇av一区二区三区| 夜夜看夜夜爽夜夜摸| 欧美日韩精品成人综合77777| 免费高清视频大片| 九色成人免费人妻av| 欧美在线一区亚洲| xxxwww97欧美| 赤兔流量卡办理| 不卡一级毛片| 在线免费观看的www视频| 校园人妻丝袜中文字幕| 成人性生交大片免费视频hd| 又黄又爽又刺激的免费视频.| 少妇的逼水好多| 在线观看舔阴道视频| 日韩,欧美,国产一区二区三区 | 成人欧美大片| 啦啦啦韩国在线观看视频| av专区在线播放| 免费一级毛片在线播放高清视频| 欧美日本视频| 国产精品国产高清国产av| 精品一区二区三区视频在线观看免费| 成人国产一区最新在线观看| 精品久久久久久久久亚洲 | 乱码一卡2卡4卡精品| 免费无遮挡裸体视频| 久久精品夜夜夜夜夜久久蜜豆| 欧美一区二区亚洲| av女优亚洲男人天堂| 欧美高清性xxxxhd video| av在线亚洲专区| 国产免费一级a男人的天堂| 色吧在线观看| 99久久无色码亚洲精品果冻| 成年免费大片在线观看| 啦啦啦观看免费观看视频高清| 亚洲专区国产一区二区| 男人舔奶头视频| 国产成年人精品一区二区| 男女视频在线观看网站免费| 99热这里只有是精品在线观看| 少妇的逼水好多| 日本黄大片高清| 男女下面进入的视频免费午夜| 欧美日韩中文字幕国产精品一区二区三区| 日本一二三区视频观看| 成人亚洲精品av一区二区| 国产色婷婷99| 悠悠久久av| 3wmmmm亚洲av在线观看| 欧美日韩乱码在线| 国产日本99.免费观看| 人妻丰满熟妇av一区二区三区| 国产视频一区二区在线看| 国产不卡一卡二| 欧美一区二区国产精品久久精品| 岛国在线免费视频观看| 欧美丝袜亚洲另类 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲无线在线观看| 亚洲av免费高清在线观看| 中文字幕人妻熟人妻熟丝袜美| 两性午夜刺激爽爽歪歪视频在线观看| av天堂在线播放| 亚洲无线在线观看| 狂野欧美白嫩少妇大欣赏| a级一级毛片免费在线观看| 午夜福利在线观看免费完整高清在 | 麻豆成人午夜福利视频| 久久久午夜欧美精品| 亚洲中文字幕一区二区三区有码在线看| 成人毛片a级毛片在线播放| 亚洲一级一片aⅴ在线观看| 深爱激情五月婷婷| 动漫黄色视频在线观看| 国产主播在线观看一区二区| 免费高清视频大片| 欧美黑人欧美精品刺激| 日韩欧美免费精品| 成人性生交大片免费视频hd| 欧美另类亚洲清纯唯美| 国产色婷婷99| 国产精品久久久久久久电影| 成人国产一区最新在线观看| 97人妻精品一区二区三区麻豆| 91麻豆精品激情在线观看国产| 黄色视频,在线免费观看| 国产成人av教育| 如何舔出高潮| 久久久久精品国产欧美久久久| 级片在线观看| 赤兔流量卡办理| 又紧又爽又黄一区二区| 91在线观看av| 波多野结衣高清无吗| 3wmmmm亚洲av在线观看| 免费电影在线观看免费观看| 精品午夜福利在线看| 欧美性猛交黑人性爽| 18禁在线播放成人免费| 一个人观看的视频www高清免费观看| 欧美丝袜亚洲另类 | 九色国产91popny在线| 免费人成在线观看视频色| 伊人久久精品亚洲午夜| 人妻制服诱惑在线中文字幕| 国产精品1区2区在线观看.| 老熟妇仑乱视频hdxx| 少妇被粗大猛烈的视频| 久久精品国产清高在天天线| 成人国产麻豆网| 精品久久国产蜜桃| 毛片一级片免费看久久久久 | 久久久成人免费电影| 国产三级在线视频| 国产精品福利在线免费观看| 亚洲熟妇中文字幕五十中出| 高清毛片免费观看视频网站| 欧美色欧美亚洲另类二区| 一级av片app| 亚洲在线观看片| av天堂中文字幕网| 精品人妻视频免费看| 亚洲av.av天堂| 久久精品综合一区二区三区| 国产人妻一区二区三区在| 中文字幕免费在线视频6| 亚洲欧美日韩高清在线视频| 亚洲最大成人中文| 色噜噜av男人的天堂激情| 长腿黑丝高跟| 极品教师在线免费播放| 亚洲va在线va天堂va国产| 亚洲男人的天堂狠狠| 91狼人影院| 日本黄大片高清| 亚洲av成人av| 天天躁日日操中文字幕| 在线免费观看不下载黄p国产 | 人人妻人人澡欧美一区二区| 精品乱码久久久久久99久播| 色综合站精品国产| 一区二区三区免费毛片| 色视频www国产| 男女那种视频在线观看| 亚洲中文字幕日韩| 精华霜和精华液先用哪个| АⅤ资源中文在线天堂| 欧美一区二区国产精品久久精品| 99国产精品一区二区蜜桃av| 人妻少妇偷人精品九色| 又爽又黄a免费视频| 国产探花极品一区二区| 精品人妻偷拍中文字幕| 嫩草影院精品99| 久久午夜亚洲精品久久| 亚洲精品乱码久久久v下载方式| 亚洲专区国产一区二区| 国产 一区精品| 久久久久久久午夜电影| 久久天躁狠狠躁夜夜2o2o| av在线老鸭窝| 国产精品98久久久久久宅男小说| 国产爱豆传媒在线观看| 久久久久久久久久黄片| 美女 人体艺术 gogo| 真实男女啪啪啪动态图| 国产视频一区二区在线看| 国产精品国产三级国产av玫瑰| 久久6这里有精品| 成人二区视频| 色综合色国产| 国产高清视频在线播放一区| 九九久久精品国产亚洲av麻豆| 亚洲五月天丁香| 成人二区视频| 真人做人爱边吃奶动态| 在线天堂最新版资源| 欧美bdsm另类| 亚洲成人中文字幕在线播放| 欧美精品国产亚洲| 男人的好看免费观看在线视频| 日本一本二区三区精品| 亚洲天堂国产精品一区在线| 免费av观看视频| 国产 一区精品| 深爱激情五月婷婷| 国产精品一区www在线观看 | 色综合婷婷激情| 身体一侧抽搐| avwww免费| 日本一二三区视频观看| 欧洲精品卡2卡3卡4卡5卡区| 成人精品一区二区免费| 国产 一区精品| 国产伦人伦偷精品视频| 精品一区二区三区视频在线| 在线免费观看不下载黄p国产 | 亚洲成人中文字幕在线播放| 亚洲欧美日韩卡通动漫| xxxwww97欧美| 麻豆一二三区av精品| 在线a可以看的网站| 99热网站在线观看| 色精品久久人妻99蜜桃| 亚洲av不卡在线观看| 一级黄片播放器| 亚洲经典国产精华液单| 91av网一区二区| 春色校园在线视频观看| 国内久久婷婷六月综合欲色啪| 熟妇人妻久久中文字幕3abv| 国产乱人视频| 日本a在线网址| 成年人黄色毛片网站| 十八禁国产超污无遮挡网站| 国产高清不卡午夜福利| 亚洲成人精品中文字幕电影| 欧美性感艳星| 麻豆精品久久久久久蜜桃| 人人妻人人看人人澡| 日本 av在线| 少妇裸体淫交视频免费看高清| 又爽又黄无遮挡网站| 国内精品久久久久精免费| 国产精品亚洲美女久久久| 人人妻人人澡欧美一区二区| 亚洲18禁久久av| 精品人妻熟女av久视频| 99精品在免费线老司机午夜| 国产一级毛片七仙女欲春2| 一a级毛片在线观看| 免费观看人在逋| 欧美xxxx性猛交bbbb| av在线老鸭窝| 国产爱豆传媒在线观看| 搡老岳熟女国产| 一本一本综合久久| 久久久久久久久久久丰满 | 中文资源天堂在线| 国内精品一区二区在线观看| 成人国产一区最新在线观看| 色综合站精品国产| 无人区码免费观看不卡| 日韩精品中文字幕看吧| 无遮挡黄片免费观看| 狠狠狠狠99中文字幕| av专区在线播放| 两个人的视频大全免费| 亚洲国产欧美人成| 免费看av在线观看网站| 在线播放无遮挡| 国产高清视频在线播放一区| 国产精品三级大全| 中出人妻视频一区二区| 少妇熟女aⅴ在线视频| 国产精品无大码| 亚洲精品粉嫩美女一区| 亚洲成a人片在线一区二区| 色综合婷婷激情| 午夜福利在线观看免费完整高清在 | 亚洲av不卡在线观看| 欧美日韩综合久久久久久 | 大型黄色视频在线免费观看| 日日啪夜夜撸| 97热精品久久久久久| 国产aⅴ精品一区二区三区波| 国产精品1区2区在线观看.| 国国产精品蜜臀av免费| 精品久久久久久久末码| 久久久久性生活片| 成人精品一区二区免费| 日本黄大片高清| 少妇的逼好多水| 人妻少妇偷人精品九色| 久久久久久久久中文| 亚洲精品一区av在线观看| 国产一区二区亚洲精品在线观看| 久久久国产成人免费| 亚洲人成网站高清观看| 亚洲国产精品合色在线| 久久亚洲精品不卡| 高清毛片免费观看视频网站| av黄色大香蕉| 99视频精品全部免费 在线| 免费电影在线观看免费观看| 一本久久中文字幕| 日韩一区二区视频免费看| 精品一区二区三区人妻视频| 国产精品国产高清国产av| 神马国产精品三级电影在线观看| 午夜福利在线观看吧| 国产成人a区在线观看| 日韩欧美国产一区二区入口| 成人无遮挡网站| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美在线二视频| 韩国av一区二区三区四区| 精品99又大又爽又粗少妇毛片 | 999久久久精品免费观看国产| 两个人视频免费观看高清| 欧洲精品卡2卡3卡4卡5卡区| 日本熟妇午夜| 日韩精品青青久久久久久| 丰满乱子伦码专区| 人妻久久中文字幕网| 成人特级黄色片久久久久久久| 女生性感内裤真人,穿戴方法视频| 中出人妻视频一区二区| 一级黄色大片毛片| 天堂√8在线中文| 日本欧美国产在线视频| 99久久精品一区二区三区| 午夜福利高清视频| 亚洲av成人av| 成人性生交大片免费视频hd| 国产黄片美女视频| 99久久中文字幕三级久久日本| 一卡2卡三卡四卡精品乱码亚洲| 一区二区三区激情视频| 日本爱情动作片www.在线观看 | 亚洲性夜色夜夜综合| 亚洲天堂国产精品一区在线| 12—13女人毛片做爰片一| АⅤ资源中文在线天堂| 尤物成人国产欧美一区二区三区| av.在线天堂| 丰满乱子伦码专区| 乱系列少妇在线播放| 欧美日本视频| 国产伦精品一区二区三区视频9| 18禁黄网站禁片午夜丰满| 午夜免费男女啪啪视频观看 | 色5月婷婷丁香| 男女边吃奶边做爰视频| 国产69精品久久久久777片| 九九在线视频观看精品| 一夜夜www| 天堂网av新在线| 大型黄色视频在线免费观看| 禁无遮挡网站| 色5月婷婷丁香| 少妇丰满av| 国产69精品久久久久777片| h日本视频在线播放| 国产真实伦视频高清在线观看 | 国内精品久久久久精免费| 99久久精品热视频| 我要搜黄色片| 免费在线观看成人毛片| 日本熟妇午夜| 亚洲乱码一区二区免费版| 噜噜噜噜噜久久久久久91| 亚洲av不卡在线观看| 国产高清视频在线播放一区| 久9热在线精品视频| 欧美日韩国产亚洲二区| 亚洲精品国产成人久久av| 欧美极品一区二区三区四区| 熟女电影av网| 97超视频在线观看视频| 尤物成人国产欧美一区二区三区| 自拍偷自拍亚洲精品老妇| 韩国av在线不卡| 亚洲经典国产精华液单| 日韩av在线大香蕉| 免费在线观看日本一区| 最好的美女福利视频网| 少妇人妻精品综合一区二区 | 免费看光身美女| 久久婷婷人人爽人人干人人爱| 在线免费十八禁| 伦理电影大哥的女人| 一个人观看的视频www高清免费观看| 成人精品一区二区免费| 精品久久国产蜜桃| 51国产日韩欧美| 亚洲中文字幕日韩| 欧美成人一区二区免费高清观看| 看黄色毛片网站| 在线观看av片永久免费下载| 听说在线观看完整版免费高清| 午夜福利成人在线免费观看| 国产一区二区三区视频了| 免费观看人在逋| 99久久精品国产国产毛片| 国产精品电影一区二区三区| 黄片wwwwww| 午夜爱爱视频在线播放| 亚洲成人久久性| 国产伦人伦偷精品视频| 尾随美女入室| 长腿黑丝高跟| 成人一区二区视频在线观看| 国产 一区精品| 能在线免费观看的黄片| 男女啪啪激烈高潮av片| 午夜老司机福利剧场| 久久久精品大字幕| 亚洲成人精品中文字幕电影| 极品教师在线免费播放| 久久人妻av系列| 成人国产一区最新在线观看| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播放欧美日韩| a级毛片免费高清观看在线播放| 国产69精品久久久久777片| 美女大奶头视频| 日韩高清综合在线| 久久久久久久久中文| 亚洲av日韩精品久久久久久密| 成人综合一区亚洲| 久久人人精品亚洲av| av在线天堂中文字幕| 精品一区二区三区av网在线观看| 亚洲精品日韩av片在线观看| 欧美日韩亚洲国产一区二区在线观看| 99视频精品全部免费 在线| 成年免费大片在线观看| 国产精品一区二区性色av| 久久6这里有精品| 国产毛片a区久久久久| 国国产精品蜜臀av免费| a级一级毛片免费在线观看| 黄片wwwwww| 一本一本综合久久| 免费观看在线日韩| 色哟哟·www| 色综合婷婷激情| 精品乱码久久久久久99久播| 黄色欧美视频在线观看| 亚洲成人中文字幕在线播放| 婷婷丁香在线五月| 男女啪啪激烈高潮av片| 有码 亚洲区| 欧美最黄视频在线播放免费| 91麻豆精品激情在线观看国产| 九色国产91popny在线| 亚洲精品日韩av片在线观看| 亚洲最大成人手机在线| 一区二区三区高清视频在线| 欧美在线一区亚洲| 免费一级毛片在线播放高清视频| 亚洲狠狠婷婷综合久久图片| 日韩欧美国产在线观看| 中文字幕av成人在线电影| 俺也久久电影网| 熟妇人妻久久中文字幕3abv| 日本爱情动作片www.在线观看 | 村上凉子中文字幕在线| 特大巨黑吊av在线直播| 欧美又色又爽又黄视频| 男人狂女人下面高潮的视频| 久久精品91蜜桃| 一个人观看的视频www高清免费观看| 一个人看的www免费观看视频| 欧美日韩乱码在线| 伦理电影大哥的女人| 我要看日韩黄色一级片| 真实男女啪啪啪动态图| 国产一区二区激情短视频| 尤物成人国产欧美一区二区三区| 88av欧美|