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

    基于單-雙高斯模型擬合法的測風(fēng)激光雷達海上風(fēng)電機組尾流特征分析

    2021-02-26 01:17:58劉清媛吳松華張凱臨李榮忠翟曉春
    關(guān)鍵詞:尾流風(fēng)場損失率

    劉清媛,吳松華,2,3?,張凱臨,3,李榮忠,3,翟曉春

    (1 中國海洋大學(xué)信息科學(xué)與工程學(xué)院,海洋遙感研究所, 山東 青島 266100;2 青島海洋科學(xué)與技術(shù)試點國家實驗室,區(qū)域海洋動力學(xué)與數(shù)值模擬功能實驗室, 山東 青島 266237;3 青島鐳測創(chuàng)芯科技有限公司, 山東 青島 266101)

    0 引 言

    近年來,隨著綠色低碳轉(zhuǎn)型政策的深入推進,新能源裝機和發(fā)電比重進一步上升,風(fēng)能作為一種清潔的可再生能源,受到了人們的廣泛關(guān)注。國家能源局報告指出,截至2019 年9 月底,我國風(fēng)電裝機容量達到1.98 億千瓦,風(fēng)電發(fā)展進入了快車道。我國領(lǐng)土廣闊,風(fēng)能資源豐富,有利于大型風(fēng)電場的建設(shè)和并網(wǎng)發(fā)電。在風(fēng)電場中,風(fēng)電機組的發(fā)電效率受到大氣風(fēng)場環(huán)境和鄰近風(fēng)電機組尾流效應(yīng)的影響,對風(fēng)電機組尾流的觀測和研究有利于優(yōu)化風(fēng)電機組布局,實現(xiàn)風(fēng)電場經(jīng)濟效益最大化。

    目前風(fēng)電行業(yè)對風(fēng)電機組尾流的研究大致分為三類: 尾流模型研究、數(shù)值模擬研究和基于實測數(shù)據(jù)的實驗研究。其中尾流模型和數(shù)值模擬研究較為成熟, 并且普遍應(yīng)用于風(fēng)電場的數(shù)據(jù)分析工作中。Jensen[1]提出的一維線性尾流模型(Park 模型)能較好地模擬平坦地形的尾流情況,是目前風(fēng)能資源評估軟件WAsP 中使用的尾流模型。Lissaman[2]在實測數(shù)據(jù)的基礎(chǔ)上提出針對風(fēng)電機組高低位置不同的尾流模型。Ainslie[3]提出渦黏性尾流模型,采用軸對稱坐標下N-S 方程求解計算尾流速度虧損場。Larsen 模型是半解析性模型,由Prandtl 旋轉(zhuǎn)對稱湍流邊界層公式推導(dǎo)得出[4]。張曉東、張鎮(zhèn)[5]結(jié)合適用于近場尾流分析的無黏尾流模型和修正后更適用于遠場尾流分析的Jensen 模型建立全場尾流模型,提高尾流模型分析的可靠性和真實性。楊祥生等[6]在傳統(tǒng)Park 模型的基礎(chǔ)上提出的Park-Gauss 模型可以很好模擬尾流區(qū)的風(fēng)速,不僅在精度上與試驗結(jié)果接近,還能分析尾流橫向的徑向風(fēng)速分布特征。Nedjari 等[7]利用CFD 對復(fù)雜地形和平坦地形下的風(fēng)機尾跡演化進行研究,CFD 在一定程度上可以模擬風(fēng)電機組周圍復(fù)雜流場。大氣風(fēng)場復(fù)雜多變,風(fēng)電機組的偏航、變槳對下風(fēng)處尾流產(chǎn)生一定影響,運用聲雷達、多普勒測風(fēng)激光雷達等新型測風(fēng)設(shè)備可以對真實風(fēng)電場中的風(fēng)電機組尾流進行觀測和研究。H¨ogstr¨om 等[8]利用高分辨率的聲雷達分析了尾流區(qū)2 D~4 D 范圍內(nèi)的風(fēng)速廓線和湍流結(jié)構(gòu)。Krishnamurthy 等[9]通過掃描型多普勒激光雷達對海上風(fēng)電機組尾流進行探測,獲取了尾流中心風(fēng)速和自由流場風(fēng)速數(shù)據(jù)計算風(fēng)速損失率。中國海洋大學(xué)Feng 等[10]運用兩臺脈沖相干多普勒激光雷達對潮間帶風(fēng)電場風(fēng)電機組尾流進行同步觀測,分析了尾流區(qū)小尺度的大氣運動和尾流擴散運動。中國海洋大學(xué)于曉慶、吳松華[11]運用脈沖相干多普勒激光雷達在新疆風(fēng)電場開展風(fēng)場探測實驗,利用單高斯擬合算法處理尾流風(fēng)速數(shù)據(jù),并分析了尾流特征參數(shù)。

    在尾流模型和數(shù)值模擬研究過程中,通過改進尾流模型和模型疊加等方式可提高尾流特征研究的準確性。在基于實測數(shù)據(jù)的實驗研究過程中,在不同風(fēng)電場開展風(fēng)場探測實驗,從實測數(shù)據(jù)中提取尾流特征參數(shù),部分學(xué)者采用單一高斯擬合模型處理尾流橫向的速度,沒有將高斯擬合算法細化。Magnusson[12]通過理論計算和測量結(jié)果提出在尾流形成初期,由于葉輪對氣流的切割作用,中間機艙對氣流的影響相對較弱,尾流橫向風(fēng)速呈雙高斯分布,隨著尾流擴散距離的增加,尾流橫向風(fēng)速逐漸呈單高斯分布。因此,僅運用傳統(tǒng)單高斯擬合算法處理風(fēng)電機組的尾流風(fēng)速數(shù)據(jù),對緊鄰風(fēng)電機組的尾流風(fēng)速分析存在誤差較大的情況。本文提出一種單–雙高斯擬合算法,結(jié)合單高斯擬合與雙高斯擬合對風(fēng)電機組尾流數(shù)據(jù)分析的優(yōu)勢,對傳統(tǒng)單高斯擬合算法進行改進,修正擬合誤差較大的結(jié)果。在江蘇某海上風(fēng)電場開展風(fēng)場探測實驗,通過對尾流寬度、風(fēng)速損失率、尾流長度和尾流中心位置等尾流特征分析,驗證了單–雙高斯擬合算法的可行性和準確性。

    1 算法原理

    運用相干多普勒測風(fēng)激光雷達獲取海上風(fēng)場的徑向風(fēng)速信息,通過VAD(Velocity azimuth display)風(fēng)場反演算法可以得到真實風(fēng)場的主風(fēng)向和背景風(fēng)速,可以進一步獲取目標風(fēng)電機組的尾流數(shù)據(jù);運用高斯擬合算法可以實現(xiàn)尾流垂直截面上橫向風(fēng)速數(shù)據(jù)的擬合分析,從而獲取尾流特征參數(shù)。

    1.1 VAD 風(fēng)場反演算法

    在相干多普勒測風(fēng)激光雷達實測數(shù)據(jù)處理過程中,假設(shè)不同高度大氣風(fēng)場是均勻的,選取同一高度上的徑向風(fēng)速數(shù)據(jù),通過正弦擬合得出激光雷達徑向風(fēng)速隨方位角的變化信息,該方法為VAD 風(fēng)場反演方法[13],其表達式為

    式中: y 為水平風(fēng)速, x 為激光雷達的掃描方位角, A、φ 為計算真實風(fēng)場的風(fēng)速和風(fēng)向的擬合參數(shù)。以圖1擬合數(shù)據(jù)為例,徑向風(fēng)速最小值ymin對應(yīng)橫坐標的激光雷達掃描方位角約為297?,此時風(fēng)場主風(fēng)向為東偏南風(fēng),徑向風(fēng)速最大值ymax計算風(fēng)場的背景風(fēng)速為6.2 m·s?1。

    圖1 VAD 風(fēng)場反演算法擬合曲線Fig.1 Fitting curve of VAD wind field inversion algorithm

    1.2 高斯模型擬合算法

    定義主風(fēng)向為風(fēng)電機組尾流的縱向,尾流沿縱向擴散,尾流橫向上的風(fēng)速為水平并且垂直于主風(fēng)向的風(fēng)速。對尾流橫向風(fēng)速進行高斯擬合之前,需將尾流區(qū)分成若干個垂直截面,獲取每個垂直截面的橫向風(fēng)速,對每一組尾流風(fēng)速進行高斯擬合處理。

    1.2.1 單高斯擬合式(2)為單高斯擬合函數(shù),式中: uLOS為尾流垂直截面的橫向風(fēng)速,u 為高斯函數(shù)的擬合參數(shù),表示尾流垂直截面外自由流場風(fēng)速,a0為擬合函數(shù)最大值與最小值之差,b0為高斯擬合峰值所在位置,即對應(yīng)尾流中心的位置,c0為高斯擬合函數(shù)的標準差。根據(jù)定義,單高斯擬合計算的尾流中心風(fēng)速為uwake= u ?a0,風(fēng)速損失率的計算公式為

    單高斯函數(shù)的標準差c0是計算尾流寬度的重要參數(shù),定義95%置信區(qū)間為尾流寬度,尾流寬度[14]約為

    1.2.2 雙高斯擬合

    在尾流形成初期,由于葉片對兩邊氣流的切割作用和中間機艙的存在,造成尾流橫向風(fēng)速值呈雙高斯分布,傳統(tǒng)尾流模型算法認為尾流橫向風(fēng)速分布為雙峰對稱的雙高斯線型,即雙高斯擬合函數(shù)的左右標準差相同。對于工作狀態(tài)的風(fēng)電機組而言,來流風(fēng)向與風(fēng)電機組掃風(fēng)平面的夾角和風(fēng)電機組自身變槳都將影響葉片后方尾流風(fēng)速分布和尾流大小,導(dǎo)致尾流中心兩邊的風(fēng)速損失不均,緊鄰風(fēng)電機組尾流垂直截面上風(fēng)速分布是不對稱、不規(guī)則的雙高斯線型,因此對傳統(tǒng)對稱雙高斯模型做出改進,具體的表達式為

    式中: a1、a2分別為雙高斯函數(shù)中兩個峰值與自由流場風(fēng)速的差值,b1、b2分別雙高斯擬合兩個峰值所在位置,c1、c2分別為左右高斯線型標準差。雙高斯擬合定義尾流風(fēng)速為uwake= u ?max(a1,a2),風(fēng)速損失率和尾流寬度計算公式分別為

    圖2 為對單個尾流垂直截面上尾流橫向風(fēng)速進行單高斯擬合和雙高斯擬合的對比結(jié)果,橫坐標為尾流橫向距離,用幾倍葉輪直徑n D(D 為風(fēng)機機組葉輪直徑,下同)表示,縱坐標為尾流橫向風(fēng)速。

    圖2 單高斯擬合與雙高斯擬合結(jié)果對比Fig.2 Comparison of single Gaussian fitting and double Gaussian fitting results

    1.3 單–雙高斯擬合算法實現(xiàn)

    單–雙高斯擬合算法在傳統(tǒng)單高斯擬合的基礎(chǔ)上,對單高斯擬合進行改進,同時加入雙高斯擬合算法,對單高斯擬合精度不足的單個尾流垂直截面數(shù)據(jù)進行修正,提高擬合精度,減小擬合誤差,單–雙高斯擬合算法流程如圖3 所示。

    1)首先將風(fēng)速數(shù)據(jù)進行歸一化處理,根據(jù)式(2)對參數(shù)進行初始化: u 定為尾流橫向上的最大風(fēng)速;a0為尾流橫向風(fēng)速最大值與最小值之差;b0為初始化尾流中心位置,設(shè)為0;令初始尾流寬度為D,c0為尾流寬度的1/4,即0.25 D;令p 為均方根誤差,當(dāng)p >0.25 時,清空初始值,對該組數(shù)據(jù)再次進行單高斯擬合。

    2)通過實測數(shù)據(jù)預(yù)處理與分析,定義單高斯擬合誤差q >0.2 或尾流寬度大于3 D 時,采用雙高斯擬合對該組尾流風(fēng)速重新處理。雙高斯擬合均方根誤差小于0.1,尾流寬度小于3 D 時擬合結(jié)束,輸出擬合參數(shù)。對于擬合均方根誤差無法小于0.1 的情況,通過循環(huán)擬合,選取最小擬合誤差的結(jié)果進行后續(xù)尾流參數(shù)計算與分析。

    圖3 算法流程圖Fig.3 Flow chart of algorithm

    2 實驗設(shè)置與數(shù)據(jù)分析

    2.1 實驗設(shè)置

    中國海洋大學(xué)激光雷達團隊于2017 年10-11 月在江蘇某海上風(fēng)電場開展風(fēng)電機組尾流觀測實驗,實驗期間平均風(fēng)速為6.29 m·s?1,激光雷達與目標風(fēng)電機組的輪轂中心相對高度為73.6 m,水平距離為663 m。

    相干多普勒測風(fēng)激光雷達利用光的多普勒效應(yīng)實現(xiàn)大氣風(fēng)場探測,激光器發(fā)出的激光束與氣溶膠粒子相互作用產(chǎn)生回波信號,回波信號光與本振光在探測器表面混頻得到多普勒頻移,由多普勒頻移即可得到探測范圍內(nèi)的徑向風(fēng)速[15]。實驗采用相干多普勒測風(fēng)激光雷達的PPI (Plane position indicator)掃描模式進行風(fēng)電機組尾流探測,PPI 模式是固定俯仰角,改變方位角的探測方式,根據(jù)實驗的需求設(shè)置方位角和俯仰角的范圍[16],數(shù)據(jù)反演結(jié)果是以激光雷達為中心的極坐標扇形圖(如圖4),色標顏色深淺表示徑向風(fēng)速大小,風(fēng)向與激光出射方向相同時,風(fēng)速為正值,反之風(fēng)速為負值。

    本次實驗設(shè)置PPI 掃描方位角為263?~360?(如圖5 所示),俯仰角為1.0?~3.6?,本次實驗數(shù)據(jù)處理過程中主要分析探測俯仰角為3.4?的尾流數(shù)據(jù),相干多普勒測風(fēng)激光雷達的技術(shù)指標如表1 所示。

    圖4 PPI 掃描模式(a)與風(fēng)場反演結(jié)果(b)Fig.4 PPI scanning mode(a)and the inversion result of wind farm(b)

    圖5 多普勒測風(fēng)激光雷達探測俯視圖Fig.5 Top view of Doppler wind lidar detection

    表1 多普勒激光雷達技術(shù)指標Table 1 Technical specifications of Doppler lidar

    2.2 單–雙高斯擬合算法分析

    選取2017 年11 月6 日第一組尾流實測數(shù)據(jù)進行風(fēng)電機組尾流的可視化分析,背景風(fēng)場主風(fēng)向為東偏南風(fēng), 結(jié)合風(fēng)電機組的偏航信息,確定風(fēng)電機組尾流的擴散方向。數(shù)據(jù)處理分析過程中, 主要研究尾流區(qū)0~5 D 范圍內(nèi)的速度場變化,把尾流區(qū)分成15 個垂直截面,每個截面距離為0.34 D。不同距離處的尾流垂直截面風(fēng)速分布如圖6 所示,橫坐標表示尾流橫向距離,坐標點(0,0)為初始尾流中心位置。從圖中可以看出尾流形成初期,尾流中心風(fēng)速值驟減,在0~2 D 距離內(nèi)尾流風(fēng)速損失較大,隨著距離的增加尾流風(fēng)速逐漸增大,尾流中心位置發(fā)生左右偏移。

    圖6 尾流垂直截面風(fēng)速變化規(guī)律Fig.6 Variation of wind speed in the vertical section of wake

    圖7 風(fēng)電機組尾流PPI 圖。(a)傳統(tǒng)單高斯擬合算法;(b)單–雙高斯擬合算法Fig.7 PPI detection of wind turbine. (a)Traditional single Gaussian fitting,(b)single-double Gaussian fitting

    如圖7 所示,風(fēng)電機組下風(fēng)處5 D 范圍內(nèi)的尾流蹤跡較為清晰,PPI 探測模式俯仰角為3.4?的掃描平面接近尾流中心所在平面,本次實驗選取風(fēng)電機組下風(fēng)處0~5 D 距離內(nèi)的尾流特征參數(shù)進行分析。圖7(a)為傳統(tǒng)單高斯擬合算法的尾流可視化示意圖,由于緊鄰風(fēng)電機組處尾流橫向風(fēng)速呈雙高斯分布,用傳統(tǒng)單高斯擬合算法提取尾流特征參數(shù)在個別尾流垂直截面存在擬合偏差較大的情況,如對第一、第二、第四個尾流垂直截面所示,由于擬合誤差較大或者產(chǎn)生無效擬合數(shù)據(jù),導(dǎo)致單個尾流截面可視化空缺和偏差,尾流寬度計算偏大。圖7(b)為運用單–雙高斯擬合算法進行尾流特征提取的可視化圖,緊鄰風(fēng)電機組的尾流垂直截面中,單–雙高斯擬合對擬合誤差較大的結(jié)果進行訂正,提高了尾流特征參數(shù)提取的有效性,從而保證了尾流特征數(shù)據(jù)的完整性。

    圖8 風(fēng)速損失率對比分析。(a)傳統(tǒng)單高斯擬合算法;(b)單–雙高斯擬合算法Fig.8 Comparative analysis of velocity deficit. (a)Traditional single Gaussian fitting,(b)single-double Gaussian fitting

    圖8 為每個尾流垂直截面對應(yīng)的風(fēng)速損失率對比分析,橫坐標為風(fēng)電機組下風(fēng)處距離,通常用幾倍的葉輪直徑表示(n D),風(fēng)電機組尾流在3 D 距離內(nèi)風(fēng)速損失率較為平穩(wěn),平均風(fēng)速損失率為83.27%,3 D 距離處開始風(fēng)速損失率下降,在5 D 距離達到22.35%。如圖8(a)所示,由于傳統(tǒng)單高斯擬合算法將擬合誤差較大的結(jié)果作為無效數(shù)據(jù),影響后期風(fēng)速損失率的統(tǒng)計分析。

    圖9 緊鄰風(fēng)電機組的尾流風(fēng)速高斯擬合結(jié)果。(a)第一個尾流垂直截面(0.34 D);(b)第二個尾流垂直截面(0.68 D);(c)第三個尾流垂直截面(1.02 D)Fig.9 Gaussian fitting of the near wake of wind turbine. (a)The first wake vertical section(0.34 D),(b)the second wake vertical section(0.68 D),(c)the third wake vertical section(1.02 D)

    選取07:11 探測的尾流數(shù)據(jù)作個例分析,前三個尾流垂直截面單高斯擬合和雙高斯擬合結(jié)果對比如圖9所示。由圖9 可以看出,尾流速度損失具有明顯的不對稱性,隨著距離的增大,尾流逐漸融合,因為風(fēng)電機組葉輪對氣流的切割作用,葉片后方尾流速度虧損較大。表2 為單高斯擬合和雙高斯擬合對緊鄰風(fēng)電機組尾流數(shù)據(jù)的尾流特征擬合結(jié)果與誤差分析,在緊鄰風(fēng)電機組處,雙高斯擬合對尾流橫向風(fēng)速數(shù)據(jù)的擬合效果優(yōu)于單高斯擬合。

    表2 尾流特征參數(shù)與誤差指標Table 2 Wake characteristic and error index

    圖10 傳統(tǒng)單高斯擬合和單–雙高斯擬合提取的參考風(fēng)速與尾流風(fēng)速對比分析Fig.10 Comparative analysis of reference wind speed and wake speed by traditional single Gaussian fitting and single-double Gaussian fitting

    參考風(fēng)速是尾流垂直截面橫向尾流區(qū)外的自由流場風(fēng)速;尾流風(fēng)速是尾流中心處的風(fēng)速值,通常情況下定義實測數(shù)據(jù)中,尾流橫向最小風(fēng)速值為尾流風(fēng)速。根據(jù)式(2)和式(5)得到參考風(fēng)速和尾流中心風(fēng)速的計算參數(shù),其中,參考風(fēng)速為u,單高斯擬合提取的尾流風(fēng)速為uwake= u ?a0,雙高斯擬合提取的尾流風(fēng)速為uwake= u ?max(a1,a2)。圖10 為運用傳統(tǒng)單高斯擬合算法和單–雙高斯擬合算法提取的參考風(fēng)速和尾流風(fēng)速對比分析。其中,Vrefer1、Vmean1、Uwake1、Umean1分別為傳統(tǒng)單高斯擬合算法提取的參考風(fēng)速、平均參考風(fēng)速、尾流風(fēng)速、平均尾流中心風(fēng)速;Vrefer2、Vmean2、Uwake2、Umean2分別為單–雙高斯擬合算法提取的參考風(fēng)速、平均參考風(fēng)速、尾流風(fēng)速、平均尾流中心風(fēng)速。尾流區(qū)外的自由流場風(fēng)速較為穩(wěn)定,尾流形成初期,由于風(fēng)電機組葉片對氣流的切割作用,葉片后方湍流強度增大,氣壓驟減,尾流風(fēng)速損失較大,隨著尾流擴散距離增加,風(fēng)電機組對下游尾流場的速度影響減小,尾流中心風(fēng)速逐漸增大。圖中可以看出,目標風(fēng)電機組尾流區(qū)0~3 D 范圍內(nèi),尾流風(fēng)速小于2 m·s?1,從3 D 距離處開始尾流風(fēng)速逐漸增大,在5 D 距離處尾流中心風(fēng)速約為3 m·s?1。

    通過平均絕對誤差(Mean absolute error,MAE)可以評估傳統(tǒng)單高斯擬合算法與單–雙高斯擬合算法對尾流中心風(fēng)速的計算偏差,其計算公式為

    圖11 尾流風(fēng)速計算誤差Fig.11 Calculation error of wake speed

    圖12 風(fēng)速損失率統(tǒng)計分析Fig.12 Statistical analysis of velocity deficit

    2.3 海上風(fēng)電機組尾流特征參數(shù)分析

    2.3.1 風(fēng)速損失率分析

    風(fēng)速損失率隨著尾流運動距離的增加而呈遞減趨勢,其計算方式參考式(6)。對00:11-16:47 時段的風(fēng)速損失率進行統(tǒng)計擬合,結(jié)果如圖12 所示,在尾流區(qū)0~2 D 范圍內(nèi),風(fēng)速損失率沒有明顯下降,反而有一定程度的攀升,根據(jù)上述尾流橫向風(fēng)速分布規(guī)律可知,此時尾流橫向風(fēng)速呈雙高斯分布,在2 D 距離之后,風(fēng)速損失率隨著距離的增加呈緩慢減小的趨勢。

    實際探測過程中,尾流長度約為6 D~8.4 D。在海上風(fēng)電場中,風(fēng)電機組下風(fēng)處3 D 范圍內(nèi)風(fēng)速損失率沒有明顯下降,平均風(fēng)速損失率為84.77%;在3 D 距離之后隨著尾流中心風(fēng)速不斷恢復(fù),風(fēng)速損失率逐漸減小,6 D 距離處風(fēng)速損失率為36.85%。定義風(fēng)速損失率減小到10%時對應(yīng)的距離為尾流長度,單–雙高斯擬合算法估算的平均尾流長度為7.92 D,由于末端尾流風(fēng)速逐漸與自由風(fēng)場風(fēng)速融合,越接近尾流末端,其風(fēng)速損失率下降越緩慢。結(jié)合尾流實測數(shù)據(jù)可視化圖與數(shù)據(jù)分析可知,在風(fēng)電機組下風(fēng)處距離大于8 D 時,尾流風(fēng)速恢復(fù)自由流場風(fēng)速,尾流效應(yīng)減弱,尾流特性不顯著,測風(fēng)激光雷達無法觀測這部分尾流特性。

    2.3.2 尾流寬度分析

    理想尾流模型假設(shè)尾流寬度隨著距離的增加呈線性增大,尾流區(qū)呈圓臺狀。實際尾流觀測實驗顯示,尾流寬度受到背景風(fēng)速和大氣湍流的影響,尾流寬度呈先增大后減小的變化趨勢,表現(xiàn)為尾流邊緣氣流與背景風(fēng)場氣流混合。圖13 為單–雙高斯擬合算法計算的每個尾流垂直截面上的尾流寬度。尾流寬度的計算參考式(7),初始尾流寬度約為1 D,尾流區(qū)0~5 D 范圍內(nèi)的平均尾流寬度為1.42 D。在風(fēng)電機組下風(fēng)處0~2 D 距離內(nèi),尾流寬度逐漸增大,2 D 距離處平均最大尾流寬度為1.57 D。尾流運動距離大于2 D 時,尾流寬度隨著距離的增大呈減小或波動趨勢,尾流運動末端寬度約為尾流形成初期寬度的1.51 倍。

    圖13 風(fēng)電機組尾流寬度統(tǒng)計分析Fig.13 Statistical analysis of wake width

    通過單–雙高斯擬合算法提取不同尾流垂直截面外的背景風(fēng)速,尾流擴散過程受背景風(fēng)場的影響,不同風(fēng)速條件下,尾流風(fēng)速損失率和尾流寬度存在差異性。以0~3 D 距離內(nèi)的尾流特征參數(shù)為例,11 月6 日實驗期間平均風(fēng)速為6.7 m·s?1,平均風(fēng)速約為5.5、6.5、7.4、8.8 m·s?1時,風(fēng)速損失率分別為82.14%、85.53%、87.59%、86.41%,平均尾流寬度分別為1.35 D、1.49 D、1.52 D、1.48 D,如圖14 所示,背景風(fēng)速小于6 m·s?1時,風(fēng)電機組尾流風(fēng)速損失較小,平均尾流寬度較小。

    2.3.3 尾流中心位置校正

    在尾流模型和數(shù)值模擬研究中,風(fēng)電機組位置固定,當(dāng)氣流經(jīng)過風(fēng)輪時,在其下游拖出的尾渦在空間上呈螺旋形渦線。在均勻大氣條件下,若將尾流看成均勻的圓臺狀,則尾流中心為圓臺中心,傳統(tǒng)定義中尾流橫向風(fēng)速最小值所在位置為尾流中心位置。實際風(fēng)場分析中發(fā)現(xiàn),緊鄰風(fēng)電機組處,尾流橫向風(fēng)速呈雙高斯分布,尾流中心為左后兩個速度峰值位置中心;隨著尾流的擴散,尾流橫向風(fēng)速呈單高斯分布,此時尾流中心為風(fēng)速最小值所在位置。

    圖14 風(fēng)速損失率與背景風(fēng)速的關(guān)系(a)和尾流寬度與背景風(fēng)速的關(guān)系(b)Fig.14 Relationship between velocity deficit and background wind speed(a)and relationship between wake width and background wind speed(b)

    圖15 為2017 年11 月6 日07:11 的風(fēng)電機組尾流數(shù)據(jù)個例分析。如圖15(a)所示,在第二、第四、第六個尾流垂直截面的尾流中心位置存在一定偏差;如圖15(b)所示,單–雙高斯擬合算法對傳統(tǒng)定義的尾流中心位置進行修正,對尾流中心位置計算較為準確。海上風(fēng)場較為穩(wěn)定,在風(fēng)電機組尾流擴散過程中,0~5 D距離內(nèi)的尾流中心不會出現(xiàn)較大的左右偏移,而在尾流末端(距離大于5 D)由于尾流中心風(fēng)速逐漸恢復(fù),尾流受周圍自由流場的氣流影響出現(xiàn)明顯的彎曲變化。

    圖15 修正尾流中心位置。(a)傳統(tǒng)定義的尾流中心位置;(b)單-雙高斯擬合計算的尾流中心位置Fig.15 Correction of wake center position. (a)Traditional single Gaussian fitting,(b)single-double Gaussian fitting

    3 結(jié) 論

    運用相干多普勒測風(fēng)激光雷達在江蘇某海上風(fēng)電場開展了風(fēng)電機組尾流觀測實驗,通過PPI 掃描模式對風(fēng)電機組尾流速度場進行了可視化探測與分析。針對傳統(tǒng)單高斯擬合算法擬合精度不足的問題,提出單–雙高斯擬合算法訂正單次擬合誤差較大的數(shù)據(jù),修正尾流中心位置,提高尾流特征參數(shù)分析的準確性,同時對海上風(fēng)電機組尾流特征進行統(tǒng)計分析,實驗觀測結(jié)果顯示:

    1)由于風(fēng)電機組葉輪對氣流的切割作用明顯,中間機艙對氣流的影響較弱,氣流經(jīng)過風(fēng)電機組掃風(fēng)平面時,下風(fēng)處氣壓驟降、湍流增強,緊鄰風(fēng)電機組葉輪的尾流橫向風(fēng)速呈雙峰不對稱的雙高斯分布;隨著距離增加,氣壓逐漸回升,尾流橫向風(fēng)速呈單高斯分布,尾流中心風(fēng)速逐漸增大。傳統(tǒng)單高斯擬合算法處理緊鄰風(fēng)電機組尾流橫向風(fēng)速數(shù)據(jù)出現(xiàn)擬合誤差較大的情況,雙高斯擬合算法處理緊鄰風(fēng)電機組尾流數(shù)據(jù)具有一定優(yōu)勢。本文提出單-雙高斯擬合算法優(yōu)化傳統(tǒng)高斯擬合算法,對傳統(tǒng)單高斯擬合算法進行參數(shù)優(yōu)化,同時對單個尾流垂直截面擬合誤差較大的結(jié)果進行修正,擬合曲線與實測尾流橫向風(fēng)速分布的相關(guān)系數(shù)可達到0.99,通過均方根誤差和尾流寬度評估一組尾流橫向風(fēng)速數(shù)據(jù)的擬合結(jié)果,可以實現(xiàn)近尾流場的尾流特征參數(shù)的提取,提高尾流速度場分析的準確性。

    2)在海上風(fēng)電機組尾流特征分析中可知,在0~3 D 范圍內(nèi)風(fēng)電機組尾流風(fēng)速損失率較大,此時尾流寬度逐漸增加,背景風(fēng)場風(fēng)速為6.5 m·s?1時,平均風(fēng)速損失率為85.53%,平均尾流寬度為1.49 D。3 D~5 D 范圍內(nèi)受周圍大氣風(fēng)場影響,尾流中心風(fēng)速逐漸增大,風(fēng)速損失率逐漸下降,在5 D 距離處的平均風(fēng)速損失率約為58.21%。在風(fēng)電機組尾流擴散過程中,0~5 D 距離內(nèi)的尾流中心無明顯偏移。通過風(fēng)速損失率的統(tǒng)計擬合分析可知,單-雙高斯擬合算法估計尾流長度約為8 D,8 D 距離后尾流逐漸恢復(fù)成自由流場風(fēng)速。

    猜你喜歡
    尾流風(fēng)場損失率
    基于FLUENT的下?lián)舯┝魅S風(fēng)場建模
    農(nóng)業(yè)農(nóng)村部印發(fā)《意見》提出到2025年農(nóng)產(chǎn)品加工環(huán)節(jié)損失率降到5%以下
    帶有治療函數(shù)及免疫損失率的SIRS流行病模型的動力學(xué)分析
    “最美風(fēng)場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    飛機尾流的散射特性與探測技術(shù)綜述
    12部使用一年后最廉價轉(zhuǎn)售車
    海外星云(2016年19期)2016-10-24 11:53:42
    2014~2015年冬季美國蜂群損失調(diào)查
    側(cè)向風(fēng)場中無人機的飛行研究
    錐形流量計尾流流場分析
    水面艦船風(fēng)尾流效應(yīng)減弱的模擬研究
    亚洲美女视频黄频| 2018国产大陆天天弄谢| 男女免费视频国产| 欧美人与善性xxx| 乱码一卡2卡4卡精品| 三上悠亚av全集在线观看| 麻豆乱淫一区二区| 大香蕉97超碰在线| 热re99久久国产66热| 亚洲精品国产色婷婷电影| 五月玫瑰六月丁香| 国产欧美另类精品又又久久亚洲欧美| 亚洲美女视频黄频| 国产熟女午夜一区二区三区| 男女啪啪激烈高潮av片| 最近最新中文字幕大全免费视频 | 在线亚洲精品国产二区图片欧美| 国产av一区二区精品久久| 新久久久久国产一级毛片| 久热久热在线精品观看| 在线免费观看不下载黄p国产| 国产免费视频播放在线视频| 午夜免费鲁丝| 制服人妻中文乱码| av天堂久久9| 爱豆传媒免费全集在线观看| 七月丁香在线播放| 日本vs欧美在线观看视频| 看非洲黑人一级黄片| 22中文网久久字幕| 国产福利在线免费观看视频| 亚洲成人av在线免费| 你懂的网址亚洲精品在线观看| av有码第一页| 免费高清在线观看日韩| 欧美激情极品国产一区二区三区 | 成人国产麻豆网| av播播在线观看一区| 日产精品乱码卡一卡2卡三| 纵有疾风起免费观看全集完整版| 久久人妻熟女aⅴ| 99国产综合亚洲精品| www.av在线官网国产| 亚洲精品美女久久av网站| 另类精品久久| 欧美人与性动交α欧美软件 | 不卡视频在线观看欧美| h视频一区二区三区| 国产又爽黄色视频| 99久久人妻综合| 亚洲国产成人一精品久久久| 国产欧美日韩一区二区三区在线| 黄网站色视频无遮挡免费观看| 99热6这里只有精品| 51国产日韩欧美| 国产黄色视频一区二区在线观看| 黄色视频在线播放观看不卡| 女性被躁到高潮视频| 国产精品嫩草影院av在线观看| 丰满少妇做爰视频| 久久午夜福利片| 在线观看免费高清a一片| 精品卡一卡二卡四卡免费| 国产精品女同一区二区软件| 狂野欧美激情性xxxx在线观看| 青青草视频在线视频观看| 中国国产av一级| 久久精品久久精品一区二区三区| 免费大片18禁| 亚洲国产日韩一区二区| 国产精品久久久av美女十八| 亚洲精品国产av蜜桃| 99久久人妻综合| 新久久久久国产一级毛片| 男人添女人高潮全过程视频| 赤兔流量卡办理| 男人舔女人的私密视频| 黄网站色视频无遮挡免费观看| 少妇熟女欧美另类| 人妻少妇偷人精品九色| 97精品久久久久久久久久精品| 亚洲五月色婷婷综合| 97人妻天天添夜夜摸| 啦啦啦啦在线视频资源| 欧美成人午夜精品| 丝袜在线中文字幕| 18禁动态无遮挡网站| 久久久久久人妻| 午夜福利视频在线观看免费| 97超碰精品成人国产| 一区二区三区乱码不卡18| 亚洲婷婷狠狠爱综合网| 中文字幕av电影在线播放| 卡戴珊不雅视频在线播放| 青青草视频在线视频观看| 91在线精品国自产拍蜜月| 亚洲,一卡二卡三卡| 午夜免费鲁丝| 国产成人免费无遮挡视频| 人人妻人人澡人人看| 国产精品 国内视频| 香蕉丝袜av| 国产av码专区亚洲av| 欧美日韩综合久久久久久| av有码第一页| 精品国产一区二区三区久久久樱花| 免费高清在线观看日韩| 国产av码专区亚洲av| 最近最新中文字幕大全免费视频 | 国产精品不卡视频一区二区| 99re6热这里在线精品视频| 在线天堂中文资源库| 少妇被粗大的猛进出69影院 | 五月开心婷婷网| 国产高清不卡午夜福利| 国产福利在线免费观看视频| 日韩av不卡免费在线播放| 成人手机av| 91aial.com中文字幕在线观看| 最新中文字幕久久久久| 一本—道久久a久久精品蜜桃钙片| 日韩中文字幕视频在线看片| 少妇的丰满在线观看| 成人国产av品久久久| 美女脱内裤让男人舔精品视频| 男女免费视频国产| 亚洲精品久久久久久婷婷小说| 国产熟女午夜一区二区三区| 日韩大片免费观看网站| 国产男人的电影天堂91| 亚洲精品国产色婷婷电影| 久久国产亚洲av麻豆专区| 中文字幕精品免费在线观看视频 | 成人亚洲欧美一区二区av| 亚洲av日韩在线播放| 两性夫妻黄色片 | 国产高清国产精品国产三级| 天美传媒精品一区二区| 精品久久久精品久久久| 久久精品夜色国产| 日日啪夜夜爽| 国产av一区二区精品久久| 亚洲,欧美,日韩| 最近中文字幕2019免费版| 中文天堂在线官网| 观看av在线不卡| xxx大片免费视频| 国产精品女同一区二区软件| 精品卡一卡二卡四卡免费| 人妻少妇偷人精品九色| 亚洲国产成人一精品久久久| 国产成人精品一,二区| 久久精品国产亚洲av涩爱| 国产av精品麻豆| 亚洲精品视频女| 亚洲欧美成人综合另类久久久| 免费av中文字幕在线| 美女内射精品一级片tv| 男人爽女人下面视频在线观看| 青春草国产在线视频| 日韩制服丝袜自拍偷拍| 一区二区日韩欧美中文字幕 | 久久99精品国语久久久| 亚洲av成人精品一二三区| 九色成人免费人妻av| 一二三四中文在线观看免费高清| a级毛片黄视频| 午夜福利,免费看| 欧美日韩av久久| 女的被弄到高潮叫床怎么办| 18禁观看日本| 秋霞伦理黄片| 国产一区二区三区综合在线观看 | 日日摸夜夜添夜夜爱| 久久久久国产精品人妻一区二区| 国产爽快片一区二区三区| 精品一品国产午夜福利视频| 少妇高潮的动态图| 国产精品国产三级专区第一集| 少妇人妻久久综合中文| 久久久久人妻精品一区果冻| 日韩欧美精品免费久久| 久久人人爽人人片av| 亚洲少妇的诱惑av| 免费在线观看完整版高清| 欧美日韩综合久久久久久| 日韩av在线免费看完整版不卡| 日韩大片免费观看网站| 最近中文字幕高清免费大全6| 最近2019中文字幕mv第一页| 晚上一个人看的免费电影| 精品一品国产午夜福利视频| 精品人妻熟女毛片av久久网站| 最近最新中文字幕大全免费视频 | 国产极品天堂在线| 少妇人妻久久综合中文| 成人漫画全彩无遮挡| 黄网站色视频无遮挡免费观看| 日日啪夜夜爽| 深夜精品福利| 欧美精品一区二区大全| 国产精品秋霞免费鲁丝片| 亚洲欧洲精品一区二区精品久久久 | av片东京热男人的天堂| 国产日韩一区二区三区精品不卡| 欧美 亚洲 国产 日韩一| 久热久热在线精品观看| 国产精品国产av在线观看| 国产精品久久久久久久电影| 成人黄色视频免费在线看| 2018国产大陆天天弄谢| 美女福利国产在线| 亚洲五月色婷婷综合| 十八禁高潮呻吟视频| 久久久久久久久久成人| 美女大奶头黄色视频| 日韩制服骚丝袜av| 久久午夜综合久久蜜桃| 中文乱码字字幕精品一区二区三区| 中文字幕制服av| 亚洲精品一二三| 一级,二级,三级黄色视频| 亚洲一级一片aⅴ在线观看| 少妇猛男粗大的猛烈进出视频| 久久精品久久久久久久性| 美女中出高潮动态图| 国产熟女欧美一区二区| 国产 一区精品| 亚洲精品国产色婷婷电影| 国产亚洲精品久久久com| 纯流量卡能插随身wifi吗| 久久国产精品男人的天堂亚洲 | 免费观看性生交大片5| 亚洲五月色婷婷综合| 巨乳人妻的诱惑在线观看| 你懂的网址亚洲精品在线观看| 男人操女人黄网站| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美一区二区三区国产| 乱码一卡2卡4卡精品| 黄色配什么色好看| 国产黄频视频在线观看| 亚洲av福利一区| 日韩免费高清中文字幕av| 丝瓜视频免费看黄片| 亚洲精品美女久久av网站| 亚洲精品久久久久久婷婷小说| 国产国拍精品亚洲av在线观看| 久久99热这里只频精品6学生| 丝袜脚勾引网站| 免费少妇av软件| 大陆偷拍与自拍| 成人黄色视频免费在线看| 久久热在线av| 久久99热6这里只有精品| 男人添女人高潮全过程视频| 欧美日韩一区二区视频在线观看视频在线| 街头女战士在线观看网站| 亚洲精品国产色婷婷电影| 最近2019中文字幕mv第一页| 欧美97在线视频| 欧美激情极品国产一区二区三区 | 精品一区二区三区视频在线| 亚洲av在线观看美女高潮| 看免费成人av毛片| 男人操女人黄网站| 黑人欧美特级aaaaaa片| 熟女人妻精品中文字幕| 成年动漫av网址| 久久人妻熟女aⅴ| 日日撸夜夜添| 一二三四中文在线观看免费高清| 超碰97精品在线观看| av视频免费观看在线观看| 满18在线观看网站| 亚洲成国产人片在线观看| 日韩欧美精品免费久久| av黄色大香蕉| 久久久国产精品麻豆| 一级片免费观看大全| 少妇的丰满在线观看| 久久 成人 亚洲| 欧美国产精品一级二级三级| 大香蕉97超碰在线| 乱码一卡2卡4卡精品| 国产成人精品在线电影| 欧美精品人与动牲交sv欧美| 久久久久精品久久久久真实原创| 久久97久久精品| 午夜91福利影院| 天堂俺去俺来也www色官网| kizo精华| 国产精品一国产av| 青春草亚洲视频在线观看| 亚洲av综合色区一区| 国产精品不卡视频一区二区| 午夜日本视频在线| 欧美精品av麻豆av| 99热全是精品| 国产精品不卡视频一区二区| 丝袜人妻中文字幕| 在线观看国产h片| 秋霞在线观看毛片| 考比视频在线观看| 国产黄色免费在线视频| 色吧在线观看| 一级毛片我不卡| 97在线人人人人妻| 欧美精品一区二区大全| 亚洲精品av麻豆狂野| 国产色婷婷99| 看非洲黑人一级黄片| a级毛片黄视频| 十八禁网站网址无遮挡| 亚洲综合色网址| 久久久久久久国产电影| 天天躁夜夜躁狠狠久久av| 亚洲高清免费不卡视频| 午夜免费男女啪啪视频观看| 极品人妻少妇av视频| 一边摸一边做爽爽视频免费| 日本免费在线观看一区| 天堂8中文在线网| 免费少妇av软件| 久久久久精品久久久久真实原创| 亚洲人与动物交配视频| 欧美最新免费一区二区三区| 91久久精品国产一区二区三区| 久久久久久久国产电影| 欧美人与善性xxx| 国产欧美日韩综合在线一区二区| 嫩草影院入口| av线在线观看网站| 人体艺术视频欧美日本| av又黄又爽大尺度在线免费看| 亚洲国产精品成人久久小说| 午夜免费观看性视频| 婷婷色麻豆天堂久久| 国产综合精华液| 女人被躁到高潮嗷嗷叫费观| 成人亚洲精品一区在线观看| 欧美日韩精品成人综合77777| 精品一区在线观看国产| 在线观看www视频免费| 成年人免费黄色播放视频| 久久精品国产自在天天线| 午夜福利视频精品| 男女午夜视频在线观看 | 啦啦啦视频在线资源免费观看| 亚洲图色成人| 熟女电影av网| 国产成人精品婷婷| 黄色视频在线播放观看不卡| 日韩成人伦理影院| 亚洲精品久久久久久婷婷小说| 日本av手机在线免费观看| 高清毛片免费看| 你懂的网址亚洲精品在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 三级国产精品片| 免费大片黄手机在线观看| 国产麻豆69| 日本欧美国产在线视频| 日韩一区二区视频免费看| 亚洲精品久久久久久婷婷小说| 最近中文字幕高清免费大全6| 九九爱精品视频在线观看| 熟女电影av网| 日韩av免费高清视频| 韩国高清视频一区二区三区| 久久精品人人爽人人爽视色| 一本色道久久久久久精品综合| 亚洲国产色片| 蜜桃在线观看..| 国产极品粉嫩免费观看在线| 99久国产av精品国产电影| 人妻 亚洲 视频| 欧美精品一区二区免费开放| 成人无遮挡网站| 成年女人在线观看亚洲视频| 日本av免费视频播放| 嫩草影院入口| 精品国产一区二区三区四区第35| 亚洲av中文av极速乱| 免费观看无遮挡的男女| 内地一区二区视频在线| 国产成人精品无人区| 国产片内射在线| 亚洲国产成人一精品久久久| 欧美日本中文国产一区发布| 日本wwww免费看| 国产片特级美女逼逼视频| 18禁国产床啪视频网站| 亚洲激情五月婷婷啪啪| 午夜免费鲁丝| 97在线人人人人妻| 久久免费观看电影| 91在线精品国自产拍蜜月| 在线观看人妻少妇| 欧美精品人与动牲交sv欧美| 啦啦啦中文免费视频观看日本| 丁香六月天网| 久久久亚洲精品成人影院| 观看美女的网站| av黄色大香蕉| 男女啪啪激烈高潮av片| 亚洲 欧美一区二区三区| 一级爰片在线观看| 2022亚洲国产成人精品| 美女福利国产在线| 精品一区二区三卡| 亚洲天堂av无毛| 亚洲av免费高清在线观看| 国产在线视频一区二区| 亚洲一级一片aⅴ在线观看| 国产成人aa在线观看| 久久久国产精品麻豆| 国产精品三级大全| 亚洲美女搞黄在线观看| 美女大奶头黄色视频| 精品少妇久久久久久888优播| 久久久a久久爽久久v久久| 久久99精品国语久久久| av卡一久久| 久久亚洲国产成人精品v| 极品人妻少妇av视频| 在线免费观看不下载黄p国产| 亚洲av电影在线进入| 男女国产视频网站| 中文字幕亚洲精品专区| 天堂中文最新版在线下载| 观看美女的网站| 91成人精品电影| 日本色播在线视频| 夜夜骑夜夜射夜夜干| 亚洲久久久国产精品| 精品久久蜜臀av无| 少妇人妻精品综合一区二区| 波野结衣二区三区在线| 国产视频首页在线观看| 天天躁夜夜躁狠狠躁躁| 天堂中文最新版在线下载| 久久人人爽av亚洲精品天堂| 亚洲国产毛片av蜜桃av| 免费人妻精品一区二区三区视频| 9191精品国产免费久久| 97超碰精品成人国产| 日韩制服骚丝袜av| 免费黄色在线免费观看| 亚洲中文av在线| 一本—道久久a久久精品蜜桃钙片| 建设人人有责人人尽责人人享有的| 99国产综合亚洲精品| 精品第一国产精品| 亚洲国产毛片av蜜桃av| av播播在线观看一区| 亚洲欧美清纯卡通| 久久久欧美国产精品| 亚洲综合色惰| 国产乱人偷精品视频| 日产精品乱码卡一卡2卡三| 久久午夜福利片| 久久亚洲国产成人精品v| 一本色道久久久久久精品综合| 久久国产精品男人的天堂亚洲 | a 毛片基地| 亚洲成色77777| 久久热在线av| 一级爰片在线观看| 视频区图区小说| 丰满饥渴人妻一区二区三| 国产av精品麻豆| 搡老乐熟女国产| 国产淫语在线视频| 69精品国产乱码久久久| 亚洲国产精品成人久久小说| 成年人免费黄色播放视频| 国产精品欧美亚洲77777| 纯流量卡能插随身wifi吗| 日产精品乱码卡一卡2卡三| 卡戴珊不雅视频在线播放| 美女中出高潮动态图| 色婷婷av一区二区三区视频| 蜜桃国产av成人99| 母亲3免费完整高清在线观看 | 又大又黄又爽视频免费| 97人妻天天添夜夜摸| 国产精品一区二区在线观看99| 成人18禁高潮啪啪吃奶动态图| 久久久a久久爽久久v久久| 国产黄频视频在线观看| 亚洲国产精品成人久久小说| 国产免费现黄频在线看| 天美传媒精品一区二区| 国产av一区二区精品久久| av片东京热男人的天堂| 国产成人欧美| 日韩人妻精品一区2区三区| 婷婷色综合www| 成人亚洲精品一区在线观看| 人人妻人人添人人爽欧美一区卜| 国产av一区二区精品久久| 国产一区亚洲一区在线观看| 免费黄网站久久成人精品| 日本与韩国留学比较| 一区二区三区乱码不卡18| 亚洲中文av在线| 青春草亚洲视频在线观看| 成年女人在线观看亚洲视频| 少妇人妻 视频| 晚上一个人看的免费电影| 欧美日韩成人在线一区二区| 熟女人妻精品中文字幕| 午夜福利网站1000一区二区三区| 丝袜美足系列| 久久久久国产网址| 久久综合国产亚洲精品| 热99久久久久精品小说推荐| 亚洲av福利一区| 国产一区有黄有色的免费视频| 精品一区二区免费观看| 如何舔出高潮| 日本av免费视频播放| 99热全是精品| 久久免费观看电影| 中文字幕精品免费在线观看视频 | av一本久久久久| 久久精品熟女亚洲av麻豆精品| 日韩三级伦理在线观看| 中文字幕制服av| 日本午夜av视频| 青春草国产在线视频| 午夜日本视频在线| 国产成人aa在线观看| 国产成人精品无人区| 春色校园在线视频观看| 五月开心婷婷网| 高清不卡的av网站| 精品少妇黑人巨大在线播放| 国产精品一区二区在线不卡| 夜夜爽夜夜爽视频| 18禁裸乳无遮挡动漫免费视频| 狠狠精品人妻久久久久久综合| 69精品国产乱码久久久| 国产亚洲一区二区精品| 母亲3免费完整高清在线观看 | 观看av在线不卡| 婷婷成人精品国产| 亚洲精品,欧美精品| 麻豆乱淫一区二区| 精品福利永久在线观看| 在线 av 中文字幕| 免费观看在线日韩| 一级毛片我不卡| 亚洲经典国产精华液单| 久久韩国三级中文字幕| 亚洲精品美女久久av网站| www.熟女人妻精品国产 | 国产精品国产三级国产av玫瑰| 久久久久久人妻| 国产国语露脸激情在线看| 飞空精品影院首页| 大香蕉久久成人网| 在线免费观看不下载黄p国产| 91午夜精品亚洲一区二区三区| 日本欧美国产在线视频| 香蕉丝袜av| av电影中文网址| 国产精品久久久久久久电影| 国产成人a∨麻豆精品| 看十八女毛片水多多多| 9色porny在线观看| 欧美+日韩+精品| 欧美日韩视频高清一区二区三区二| 欧美老熟妇乱子伦牲交| 国产综合精华液| 精品少妇黑人巨大在线播放| 成人影院久久| 免费在线观看完整版高清| 黄色怎么调成土黄色| 99国产精品免费福利视频| 久久青草综合色| 亚洲中文av在线| 免费观看无遮挡的男女| 国产成人精品福利久久| 丝袜美足系列| 久久久久久久久久人人人人人人| 国产在视频线精品| 亚洲欧美精品自产自拍| 老司机影院成人| 精品人妻偷拍中文字幕| 少妇熟女欧美另类| 1024视频免费在线观看| 99热全是精品| 久久99精品国语久久久| 国产无遮挡羞羞视频在线观看| freevideosex欧美| 亚洲精品成人av观看孕妇| 久久免费观看电影| a 毛片基地| 天天操日日干夜夜撸| 99久久精品国产国产毛片| 满18在线观看网站| 美女国产高潮福利片在线看| 男的添女的下面高潮视频| 欧美成人午夜精品| 成人黄色视频免费在线看| 中文字幕免费在线视频6| 日韩制服骚丝袜av| 乱码一卡2卡4卡精品| 国产精品人妻久久久影院| 成年女人在线观看亚洲视频| 精品一区在线观看国产| 午夜福利,免费看| 久久精品久久精品一区二区三区|