• <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
    水面艦船風尾流效應減弱的模擬研究
    欧美少妇被猛烈插入视频| 久久鲁丝午夜福利片| 久久精品久久久久久久性| 久久国内精品自在自线图片| 久久97久久精品| 国产综合精华液| 欧美日韩视频精品一区| 亚洲伊人久久精品综合| 男女边吃奶边做爰视频| 免费观看无遮挡的男女| 综合色丁香网| 一本大道久久a久久精品| 水蜜桃什么品种好| 欧美3d第一页| 日韩精品免费视频一区二区三区 | 水蜜桃什么品种好| 飞空精品影院首页| 九九爱精品视频在线观看| 男女无遮挡免费网站观看| 欧美精品高潮呻吟av久久| 亚洲人成网站在线观看播放| 伦精品一区二区三区| 亚洲中文av在线| 国产免费一级a男人的天堂| 永久免费av网站大全| 狂野欧美白嫩少妇大欣赏| 9色porny在线观看| 丰满迷人的少妇在线观看| 精品久久蜜臀av无| 亚洲精品乱久久久久久| 一级二级三级毛片免费看| 久久午夜福利片| 麻豆精品久久久久久蜜桃| 国产在线一区二区三区精| 精品亚洲乱码少妇综合久久| 人妻 亚洲 视频| 九草在线视频观看| 欧美精品人与动牲交sv欧美| 精品人妻偷拍中文字幕| 在线观看免费日韩欧美大片 | 91精品国产九色| 午夜免费观看性视频| 在线天堂最新版资源| 大香蕉久久成人网| 日韩制服骚丝袜av| 成人漫画全彩无遮挡| 春色校园在线视频观看| 女的被弄到高潮叫床怎么办| 国产色爽女视频免费观看| 午夜免费男女啪啪视频观看| 日韩熟女老妇一区二区性免费视频| 中文天堂在线官网| 日本av免费视频播放| 一本色道久久久久久精品综合| 亚洲情色 制服丝袜| 青青草视频在线视频观看| 国产在视频线精品| 一级黄片播放器| 在现免费观看毛片| 婷婷色综合www| 女性被躁到高潮视频| 一本一本综合久久| 简卡轻食公司| 三级国产精品片| 热re99久久精品国产66热6| 蜜桃久久精品国产亚洲av| 人妻制服诱惑在线中文字幕| 交换朋友夫妻互换小说| 日本与韩国留学比较| 亚洲综合色惰| 91精品伊人久久大香线蕉| 一本大道久久a久久精品| 精品国产国语对白av| 国产精品嫩草影院av在线观看| 九色亚洲精品在线播放| 嘟嘟电影网在线观看| 97在线人人人人妻| 老熟女久久久| 午夜视频国产福利| 国产在线视频一区二区| 国产免费现黄频在线看| 国产永久视频网站| 国产永久视频网站| 成人黄色视频免费在线看| 亚洲av国产av综合av卡| 国产精品成人在线| 久久99蜜桃精品久久| 69精品国产乱码久久久| 亚洲一区二区三区欧美精品| 哪个播放器可以免费观看大片| 制服诱惑二区| 免费观看a级毛片全部| 七月丁香在线播放| 婷婷色综合大香蕉| 91国产中文字幕| 欧美精品一区二区免费开放| 97超视频在线观看视频| 免费观看无遮挡的男女| 中文乱码字字幕精品一区二区三区| 最近的中文字幕免费完整| 欧美 日韩 精品 国产| 亚洲图色成人| 久久精品国产鲁丝片午夜精品| 丝袜喷水一区| 亚洲欧美成人精品一区二区| 永久网站在线| 日产精品乱码卡一卡2卡三| 国产欧美另类精品又又久久亚洲欧美| 国产综合精华液| 日本欧美国产在线视频| 久久精品夜色国产| 久久精品国产亚洲av天美| 乱人伦中国视频| 亚洲丝袜综合中文字幕| 免费av不卡在线播放| 亚洲精品视频女| xxx大片免费视频| 嫩草影院入口| 如何舔出高潮| kizo精华| 日本黄大片高清| 亚洲av男天堂| 精品少妇久久久久久888优播| 亚洲人与动物交配视频| 老女人水多毛片| 欧美日本中文国产一区发布| 91午夜精品亚洲一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 久久午夜福利片| 纯流量卡能插随身wifi吗| 亚洲欧美成人精品一区二区| 两个人的视频大全免费| 日韩大片免费观看网站| 热re99久久国产66热| 久热久热在线精品观看| 国产成人a∨麻豆精品| 老司机影院成人| 成人二区视频| 欧美日本中文国产一区发布| 国产熟女午夜一区二区三区 | 狂野欧美激情性bbbbbb| 18禁观看日本| 91精品三级在线观看| freevideosex欧美| 久久久精品94久久精品| 亚洲国产色片| h视频一区二区三区| 亚洲精品aⅴ在线观看| 中国三级夫妇交换| 亚洲综合精品二区| 另类精品久久| 青春草视频在线免费观看| 精品国产一区二区三区久久久樱花| 国产免费福利视频在线观看| 日韩三级伦理在线观看| 亚洲精品日本国产第一区| 又黄又爽又刺激的免费视频.| 免费av不卡在线播放| 街头女战士在线观看网站| av免费在线看不卡| 亚洲情色 制服丝袜| 男男h啪啪无遮挡| 美女cb高潮喷水在线观看| 亚洲精品成人av观看孕妇| 飞空精品影院首页| 老女人水多毛片| 欧美激情极品国产一区二区三区 | 欧美+日韩+精品| 校园人妻丝袜中文字幕| 亚洲精品成人av观看孕妇| kizo精华| 18在线观看网站| 成人18禁高潮啪啪吃奶动态图 | 久久精品国产a三级三级三级| 国产欧美亚洲国产| 亚洲欧美成人综合另类久久久| 久久综合国产亚洲精品| 熟女人妻精品中文字幕| 色婷婷久久久亚洲欧美| 777米奇影视久久| 中文精品一卡2卡3卡4更新| 国产成人av激情在线播放 | 久久人人爽人人爽人人片va| 成人漫画全彩无遮挡| 国产成人aa在线观看| 久久久国产欧美日韩av| 国产午夜精品一二区理论片| 岛国毛片在线播放| 久久免费观看电影| 五月伊人婷婷丁香| 久久久精品94久久精品| 久久97久久精品| 九色亚洲精品在线播放| 亚洲人成网站在线播| a级片在线免费高清观看视频| 国产精品一区二区三区四区免费观看| 亚洲精品亚洲一区二区| 哪个播放器可以免费观看大片| 自拍欧美九色日韩亚洲蝌蚪91| 国产男女内射视频| 国产精品久久久久久av不卡| 美女cb高潮喷水在线观看| 少妇人妻精品综合一区二区| 国产69精品久久久久777片| 自线自在国产av| av视频免费观看在线观看| 看免费成人av毛片| 欧美最新免费一区二区三区| 亚洲av中文av极速乱| 69精品国产乱码久久久| 中文字幕最新亚洲高清| 国产免费视频播放在线视频| 五月天丁香电影| 麻豆精品久久久久久蜜桃| 国产色爽女视频免费观看| 五月开心婷婷网| 哪个播放器可以免费观看大片| 欧美另类一区| 欧美精品亚洲一区二区| 欧美成人午夜免费资源| 亚洲国产精品999| 亚洲精品,欧美精品| 成人影院久久| 啦啦啦啦在线视频资源| 91精品一卡2卡3卡4卡| freevideosex欧美| 久久影院123| 建设人人有责人人尽责人人享有的| 久久亚洲国产成人精品v| 成人黄色视频免费在线看| 久久 成人 亚洲| 黑人巨大精品欧美一区二区蜜桃 | 日日摸夜夜添夜夜添av毛片| 精品国产露脸久久av麻豆| 18禁裸乳无遮挡动漫免费视频| 免费观看的影片在线观看| 在线观看一区二区三区激情| 高清黄色对白视频在线免费看| 精品国产乱码久久久久久小说| av卡一久久| 国产精品麻豆人妻色哟哟久久| 三级国产精品欧美在线观看| 国产精品免费大片| 狂野欧美白嫩少妇大欣赏| 少妇的逼水好多| 久久久亚洲精品成人影院| 一边亲一边摸免费视频| 亚洲av免费高清在线观看| 乱码一卡2卡4卡精品| 全区人妻精品视频| 黄色视频在线播放观看不卡| 欧美97在线视频| 男女免费视频国产| 亚洲精品乱久久久久久| 高清欧美精品videossex| 精品人妻一区二区三区麻豆| 亚洲经典国产精华液单| 久久精品国产自在天天线| av在线播放精品| 乱人伦中国视频| 大陆偷拍与自拍| 我的女老师完整版在线观看| 国模一区二区三区四区视频| 十八禁高潮呻吟视频| 国产成人精品在线电影| 国产极品粉嫩免费观看在线 | 亚洲国产成人一精品久久久| 在线精品无人区一区二区三| 日韩人妻高清精品专区| 色哟哟·www| 亚洲四区av| 久久精品熟女亚洲av麻豆精品| 国产女主播在线喷水免费视频网站| 高清av免费在线| 日韩 亚洲 欧美在线| 成人午夜精彩视频在线观看| 久久久久久久久久久久大奶| 午夜久久久在线观看| 国产黄色视频一区二区在线观看| 免费人成在线观看视频色| 成人毛片60女人毛片免费| 又粗又硬又长又爽又黄的视频| 亚洲国产色片| 亚洲人成网站在线观看播放| 欧美日韩成人在线一区二区| 青春草亚洲视频在线观看| 久久人人爽人人片av| 亚洲欧美成人精品一区二区| 成人手机av| 精品久久久精品久久久| 欧美一级a爱片免费观看看| 人妻夜夜爽99麻豆av| 一级片'在线观看视频| 在线免费观看不下载黄p国产| 国产高清不卡午夜福利| 久久午夜福利片| 欧美日本中文国产一区发布| 午夜免费鲁丝| 极品人妻少妇av视频| 亚洲人与动物交配视频| 日韩人妻高清精品专区| 国产精品一区www在线观看| 国产精品久久久久久久电影| 亚洲综合精品二区| 国产av码专区亚洲av| 美女福利国产在线| 看十八女毛片水多多多| 久久国产亚洲av麻豆专区| 伦精品一区二区三区| 国产精品欧美亚洲77777| 精品亚洲成国产av| 色吧在线观看| 少妇精品久久久久久久| 亚洲精品色激情综合| 免费黄网站久久成人精品| 欧美日韩在线观看h| 韩国av在线不卡| 纵有疾风起免费观看全集完整版| 亚洲av免费高清在线观看| 99视频精品全部免费 在线| 中文字幕av电影在线播放| 国产成人一区二区在线| 日韩成人av中文字幕在线观看| 国产又色又爽无遮挡免| 久久久精品94久久精品| 午夜免费男女啪啪视频观看| 亚洲精品国产色婷婷电影| 久久婷婷青草| 有码 亚洲区| 国产免费又黄又爽又色| 大片电影免费在线观看免费| 一级毛片aaaaaa免费看小| 中国三级夫妇交换| 亚洲少妇的诱惑av| 国产精品国产三级国产av玫瑰| 亚洲av在线观看美女高潮| 日韩伦理黄色片| 蜜桃国产av成人99| 只有这里有精品99| 国产亚洲一区二区精品| 国产午夜精品久久久久久一区二区三区| 国产精品99久久99久久久不卡 | 亚洲精品中文字幕在线视频| 2022亚洲国产成人精品| 99热这里只有是精品在线观看| 婷婷色麻豆天堂久久| 国产黄色视频一区二区在线观看| 亚洲精品亚洲一区二区| 一级黄片播放器| 亚洲精品视频女| 国产黄色视频一区二区在线观看| 午夜日本视频在线| 麻豆精品久久久久久蜜桃| 18在线观看网站| 欧美xxxx性猛交bbbb| 久久热精品热| 久久精品国产亚洲网站| 91aial.com中文字幕在线观看| 国产亚洲一区二区精品| 免费观看a级毛片全部| 国产精品久久久久久精品古装| 汤姆久久久久久久影院中文字幕| 在线观看www视频免费| 18在线观看网站| 韩国av在线不卡| 最新的欧美精品一区二区| 亚洲欧美日韩另类电影网站| 亚洲精品久久久久久婷婷小说| 少妇 在线观看| 国产精品 国内视频| 亚洲国产精品成人久久小说| 男女啪啪激烈高潮av片| 国产视频内射| 美女主播在线视频| 乱人伦中国视频| 亚洲伊人久久精品综合| 一级爰片在线观看| 国产国拍精品亚洲av在线观看| 男人添女人高潮全过程视频| 黑人猛操日本美女一级片| 一级毛片aaaaaa免费看小| 精品久久国产蜜桃| 欧美激情 高清一区二区三区| 中国国产av一级| 色5月婷婷丁香| 视频中文字幕在线观看| 国产精品一区二区三区四区免费观看| 一区二区三区免费毛片| 精品熟女少妇av免费看| 水蜜桃什么品种好| 汤姆久久久久久久影院中文字幕| 亚洲精品456在线播放app| 美女主播在线视频| 高清在线视频一区二区三区| 国产精品人妻久久久影院| 精品一区二区三卡| 91久久精品国产一区二区三区| 亚洲人成网站在线观看播放| 综合色丁香网| 国产免费现黄频在线看| 我要看黄色一级片免费的| 日本-黄色视频高清免费观看| 亚洲四区av| 天天操日日干夜夜撸| 狂野欧美白嫩少妇大欣赏| av在线观看视频网站免费| 久久久精品区二区三区| 国产极品粉嫩免费观看在线 | 日本色播在线视频| 亚洲av在线观看美女高潮| 色吧在线观看| 一区二区日韩欧美中文字幕 | videossex国产| 日本免费在线观看一区| 亚洲欧美一区二区三区国产| 汤姆久久久久久久影院中文字幕| 日韩中字成人| 久久久久久久国产电影| 十分钟在线观看高清视频www| 在线观看一区二区三区激情| 男女国产视频网站| 久久鲁丝午夜福利片| 卡戴珊不雅视频在线播放| 99精国产麻豆久久婷婷| 亚洲国产精品成人久久小说| 久久99精品国语久久久| 亚洲av不卡在线观看| 又大又黄又爽视频免费| 亚洲av成人精品一区久久| 又粗又硬又长又爽又黄的视频| 考比视频在线观看| 人人澡人人妻人| 亚洲中文av在线| 午夜福利影视在线免费观看| 久久热精品热| 欧美亚洲日本最大视频资源| 美女国产视频在线观看| 99久久精品一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 丝袜喷水一区| 丰满乱子伦码专区| 99国产精品免费福利视频| 久久韩国三级中文字幕| 国产成人精品久久久久久| 亚洲欧美一区二区三区国产| 一区在线观看完整版| 久久精品国产a三级三级三级| 国产欧美另类精品又又久久亚洲欧美| 久久毛片免费看一区二区三区| 青青草视频在线视频观看| 成人免费观看视频高清| 日韩亚洲欧美综合| 18在线观看网站| 欧美日韩亚洲高清精品| 爱豆传媒免费全集在线观看| 国产又色又爽无遮挡免| 啦啦啦视频在线资源免费观看| 国产熟女欧美一区二区| 午夜激情福利司机影院| 国产极品粉嫩免费观看在线 | 搡女人真爽免费视频火全软件| tube8黄色片| 色哟哟·www| 在线精品无人区一区二区三| 一级毛片我不卡| 亚洲av欧美aⅴ国产| 插阴视频在线观看视频| 在线天堂最新版资源| 日韩 亚洲 欧美在线| 成人黄色视频免费在线看| 精品亚洲乱码少妇综合久久| 亚洲欧洲日产国产| 看免费成人av毛片| 免费观看的影片在线观看| 亚洲精品亚洲一区二区| 亚洲人成网站在线播| 婷婷色综合www| 你懂的网址亚洲精品在线观看| 国产黄色视频一区二区在线观看| 精品久久久久久久久亚洲| 丝袜喷水一区| 99九九线精品视频在线观看视频| 女人久久www免费人成看片| 欧美 亚洲 国产 日韩一| 久久久午夜欧美精品| 久久久久久久国产电影| 又粗又硬又长又爽又黄的视频| 欧美激情 高清一区二区三区| 91精品伊人久久大香线蕉| 在线看a的网站| 人妻系列 视频| 免费播放大片免费观看视频在线观看| 男女无遮挡免费网站观看| 欧美日韩视频高清一区二区三区二| 国产高清不卡午夜福利| 狂野欧美激情性bbbbbb| 成人国产麻豆网| 日韩精品免费视频一区二区三区 | 午夜福利影视在线免费观看| 亚洲精品日韩在线中文字幕| 亚洲美女搞黄在线观看| 亚洲内射少妇av| 欧美日本中文国产一区发布| 色婷婷久久久亚洲欧美| 哪个播放器可以免费观看大片| tube8黄色片| 丰满迷人的少妇在线观看| 久久国产精品大桥未久av| 2021少妇久久久久久久久久久| 精品熟女少妇av免费看| 亚洲国产色片| 少妇丰满av| 国产免费又黄又爽又色| 一级a做视频免费观看| 国产av一区二区精品久久| 久久综合国产亚洲精品| 中文字幕人妻丝袜制服| 热re99久久精品国产66热6| 色吧在线观看| 人人妻人人爽人人添夜夜欢视频| 一本大道久久a久久精品| 免费人妻精品一区二区三区视频| 91成人精品电影| 伦理电影免费视频| 亚洲成人av在线免费| 秋霞在线观看毛片| 菩萨蛮人人尽说江南好唐韦庄| 欧美变态另类bdsm刘玥| 王馨瑶露胸无遮挡在线观看| 久久久午夜欧美精品| 十八禁高潮呻吟视频| 中文字幕最新亚洲高清| 插逼视频在线观看| 天堂俺去俺来也www色官网| 女人精品久久久久毛片| 中国美白少妇内射xxxbb| 精品久久久精品久久久| 久久99蜜桃精品久久| tube8黄色片| 人人妻人人添人人爽欧美一区卜| 日本黄大片高清| 国产精品蜜桃在线观看| 欧美精品国产亚洲| 国产在线免费精品| 最近2019中文字幕mv第一页| 精品久久久久久电影网| 国产精品不卡视频一区二区| 美女cb高潮喷水在线观看| 国产一区二区三区av在线| 亚洲第一av免费看| 一区二区三区四区激情视频| 22中文网久久字幕| 亚洲国产精品国产精品| 成人漫画全彩无遮挡| 肉色欧美久久久久久久蜜桃| 美女视频免费永久观看网站| 亚洲内射少妇av| 国产免费现黄频在线看| 国产日韩欧美视频二区| 爱豆传媒免费全集在线观看| 日韩,欧美,国产一区二区三区| av女优亚洲男人天堂| 国产极品粉嫩免费观看在线 | 亚洲av福利一区| 国产成人freesex在线| 简卡轻食公司| 九九久久精品国产亚洲av麻豆| 亚洲精华国产精华液的使用体验| 亚洲欧美清纯卡通| 伊人亚洲综合成人网| 国产视频内射| av播播在线观看一区| 美女视频免费永久观看网站| 国产熟女午夜一区二区三区 | 亚洲av.av天堂| 男人添女人高潮全过程视频| 搡老乐熟女国产| 国产精品.久久久| 国语对白做爰xxxⅹ性视频网站| 国产有黄有色有爽视频| 91久久精品国产一区二区成人| 免费av不卡在线播放| 在现免费观看毛片| 超色免费av| 国产爽快片一区二区三区| 成人国产av品久久久| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 免费看av在线观看网站| 亚洲欧美成人综合另类久久久| 亚洲美女搞黄在线观看| 日韩制服骚丝袜av| 国产色爽女视频免费观看| 免费久久久久久久精品成人欧美视频 | 精品国产乱码久久久久久小说| 中文字幕免费在线视频6| 亚洲美女视频黄频| 中文字幕人妻熟人妻熟丝袜美| 一本色道久久久久久精品综合| 亚洲三级黄色毛片| 丰满乱子伦码专区| 国产一区二区三区综合在线观看 | 国产精品一区www在线观看| 国产精品偷伦视频观看了| 亚洲图色成人| 一边亲一边摸免费视频| 国产成人一区二区在线| 王馨瑶露胸无遮挡在线观看| 午夜福利网站1000一区二区三区| 秋霞在线观看毛片| 日韩精品有码人妻一区| 国产色婷婷99| 午夜免费鲁丝| 人体艺术视频欧美日本| 中文乱码字字幕精品一区二区三区|