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

    HY-2B衛(wèi)星散射計神經(jīng)網(wǎng)絡(luò)多區(qū)間風(fēng)速反演

    2021-05-21 07:58:52上海海洋大學(xué)海洋科學(xué)學(xué)院上海201306上海河口海洋測繪工程技術(shù)研究中心上海201306
    海洋科學(xué)進展 2021年2期
    關(guān)鍵詞:散射系數(shù)風(fēng)場反演

    (1.上海海洋大學(xué) 海洋科學(xué)學(xué)院,上海201306;2.上海河口海洋測繪工程技術(shù)研究中心,上海201306)

    郭 鑫1,韓 震1,2*,張雪薇1,周瑋辰1

    海面風(fēng)場是研究海洋運動和氣候變化重要的物理參數(shù),獲取高質(zhì)量的海洋風(fēng)場資料具有重要的意義。星載散射計是一種搭載在衛(wèi)星上的主動微波傳感器,是獲得大范圍、高時空分辨率海表風(fēng)矢量最有效的方法之一。散射計發(fā)射出的電磁波束到達海面后,電磁波會發(fā)生后向散射現(xiàn)象。風(fēng)通過調(diào)制海面的粗糙度來影響散射計的后向散射系數(shù),而后向散射系數(shù)是風(fēng)速、風(fēng)向、入射角、方位角和極化方式的函數(shù),可以建立后向散射系數(shù)與海面風(fēng)矢量之間的關(guān)系,實現(xiàn)對海面風(fēng)場的反演[1]。目前運行的散射計主要工作在C和K波段,很多地球物理模型函數(shù)(Geophysical Model Function,GMF)也是基于這2個波段建立的,比如CMOD系列和NSCAT系列[2-5]。人們對風(fēng)矢量和海面幾何狀況之間的關(guān)系,以及電磁波與海表面之間的相互作用機制還未研究透徹,因此這些GMF都是半經(jīng)驗?zāi)P秃瘮?shù)。在臺風(fēng)等極端情況,受海表面物理特性以及環(huán)境的復(fù)雜影響,后向散射系數(shù)與風(fēng)場之間呈高度非線性關(guān)系,致使所建立的地球物理模型函數(shù)的普適性并不理想。

    由于利用神經(jīng)網(wǎng)絡(luò)方法反演風(fēng)速不依賴復(fù)雜的GMF,因此具有簡潔高效和適應(yīng)不同海況的能力[6]。Thiria等[7]于1993年將神經(jīng)網(wǎng)絡(luò)方法運用于散射計風(fēng)場反演,證明了復(fù)雜傳遞函數(shù)可以通過使用神經(jīng)網(wǎng)絡(luò)有效地建模。宋新改等[8]于2006年基于ERS-1,ERS-2數(shù)據(jù),采用神經(jīng)網(wǎng)絡(luò)算法對風(fēng)場進行反演,推動了國內(nèi)對神經(jīng)網(wǎng)絡(luò)反演海面風(fēng)場的研究。Xie等[9]于2008年在考慮了海面溫度對雷達后向散射系數(shù)的影響下,建立了一種基于神經(jīng)網(wǎng)絡(luò)的地球物理模型函數(shù)。Stiles等[10]于2014年開發(fā)了一種神經(jīng)網(wǎng)絡(luò)技術(shù),用于從Ku波段散射計的測量數(shù)據(jù)中提取精確的12.5 km空間分辨率的風(fēng)速。王婧等[11]在2017年使用BP神經(jīng)網(wǎng)絡(luò)對高風(fēng)速風(fēng)場GMF建模,并通過神經(jīng)網(wǎng)絡(luò)的訂正來修正降雨對風(fēng)速反演的影響。Xu等[12]2018年基于神經(jīng)網(wǎng)絡(luò)方法,聯(lián)合HY-2A衛(wèi)星主被動傳感器(散射計和輻射計),實現(xiàn)了臺風(fēng)條件下的風(fēng)場反演。

    近年來,神經(jīng)網(wǎng)絡(luò)相關(guān)算法不斷被優(yōu)化,反向傳播算法從梯度下降法、Momentum、RMSprop到Adam,網(wǎng)絡(luò)訓(xùn)練的效果不斷提升。隨著同化技術(shù)的發(fā)展,再分析資料的質(zhì)量也在提升。研究發(fā)現(xiàn),歐洲中期天氣預(yù)報中心(European Centre for Medium-Range Weather Forecasts,ECWMF)推出的ERA5再分析資料的適用性總體上要優(yōu)于ERA-Interim[13]。本文針對當(dāng)西太平洋地區(qū)存在臺風(fēng)的情況,基于HY-2B衛(wèi)星散射計的L2A數(shù)據(jù)和ECWMF的ERA5數(shù)據(jù),利用反向傳播神經(jīng)網(wǎng)絡(luò)對其進行訓(xùn)練,建立多風(fēng)速段的風(fēng)速反演模型,并與基于NSCAT-4 GMF得的L2B風(fēng)速(以下簡稱L2B風(fēng)速,vL2B)進行了比較分析,以期通過建立合適的反演模型提高臺風(fēng)存在情況下的衛(wèi)星散射計反演風(fēng)速質(zhì)量。

    1 數(shù)據(jù)處理

    1.1 HY-2B衛(wèi)星數(shù)據(jù)

    HY-2B衛(wèi)星于2018-10-25成功發(fā)射升空,是我國第二顆海洋動力環(huán)境衛(wèi)星。所搭載的散射計為筆形波束圓錐掃描體制,具有內(nèi)、外兩個波束,內(nèi)波束為HH極化,入射角為41.4°,外波束為VV極化,入射角為48.5°。Wang等[14]發(fā)現(xiàn)HY-2B散射計風(fēng)場精度優(yōu)于OSCAT2和HY-2A散射計,與ASCAT相近。HY-2B散射計的L2數(shù)據(jù)由L2A和L2B數(shù)據(jù)構(gòu)成,兩者都以軌道為單位進行存儲。L2A數(shù)據(jù)包含了風(fēng)矢量單元(Wind Vector Cells,WVC)中的后向散射系數(shù),以及相對應(yīng)的入射角、方位角等信息。L2B數(shù)據(jù)包含了L2A數(shù)據(jù)基于最大似然估計方法反演得到的海面風(fēng)場信息。本文中,L2A數(shù)據(jù)經(jīng)過預(yù)處理后,用于神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練;L2B數(shù)據(jù)用于模型效果的檢驗。對于L2A數(shù)據(jù),本文做兩步預(yù)處理:1)統(tǒng)一WVC參數(shù)數(shù)量。由于每個WVC所觀測的后向散射系數(shù)數(shù)目不同,為了滿足神經(jīng)網(wǎng)絡(luò)輸入特征一致,根據(jù)觀測情況(內(nèi)波束前視、內(nèi)波束后視、外波束前視和外波束后視)對每個WVC的后向散射系數(shù)、入射角、以及方位角進行求和平均,使得每個單元格中具有相等數(shù)量的觀測值以及儀器狀態(tài)信息。2)剔除不滿足實驗要求數(shù)據(jù)。對于后向散射系數(shù)數(shù)目小于4以及沒有同時擁有4種觀測情況的單元格選擇剔除;在外部條帶(風(fēng)矢量單元格的列索引值為[1,10]和[65,76]),只有VV極化測量,與條帶的其他部分相比,質(zhì)量下降,選擇剔除。經(jīng)過處理后,將L2A數(shù)據(jù)與L2B數(shù)據(jù)進行匹配,最終得到以風(fēng)速、風(fēng)向、后向散射系數(shù)、電磁波方位角、入射角為一組的同一時刻和經(jīng)緯度的單元格數(shù)據(jù)。

    1.2 ERA5風(fēng)速數(shù)據(jù)

    衛(wèi)星海面風(fēng)場數(shù)據(jù)和現(xiàn)場觀測海面風(fēng)場數(shù)據(jù)各有其利弊。衛(wèi)星數(shù)據(jù)具有獲得大范圍、高時空分辨率觀測值的特點,但精度不如現(xiàn)場觀測數(shù)據(jù);而現(xiàn)場觀測數(shù)據(jù)雖然精度高,但數(shù)據(jù)量小(如浮標(biāo)數(shù)據(jù)),在時間和空間上連續(xù)性差,很難滿足神經(jīng)網(wǎng)絡(luò)訓(xùn)練時的數(shù)據(jù)量要求。

    ERA5再分析數(shù)據(jù)是歐洲中期天氣預(yù)報中心(European Centre for Medium Range Weather Forecasts,ECMWF)發(fā)布的第四代全球大氣再分析產(chǎn)品,同化了衛(wèi)星數(shù)據(jù)、現(xiàn)場觀測數(shù)據(jù)等多種數(shù)據(jù)產(chǎn)品,是目前最為先進的再分析數(shù)據(jù)產(chǎn)品之一。陳克海等[15]對HY-2B風(fēng)場數(shù)據(jù)與ECMWF再分析風(fēng)場數(shù)據(jù)、浮標(biāo)風(fēng)場數(shù)據(jù)進行了對比,發(fā)現(xiàn)HY-2B風(fēng)速與他們之間的誤差特性基本一致。使用ERA5再分析數(shù)據(jù)的海平面10 m風(fēng)速作為模型訓(xùn)練時的真值(以下簡稱ERA5風(fēng)速,vERA5),具有較高精度和較高時空分辨率。由于ERA5再分析數(shù)據(jù)并沒有同化HY-2B散射計數(shù)據(jù),因此利用ERA5作為訓(xùn)練真值能保證數(shù)據(jù)的相對獨立。

    本文使用ERA5再分析數(shù)據(jù)的10 m經(jīng)向風(fēng)速和10 m緯向風(fēng)速,其空間分辨率與時間分辨率分別為1 h和0.25°,經(jīng)矢量合成后得到海平面10 m風(fēng)速,再與L2A數(shù)據(jù)匹配后形成訓(xùn)練集,用于本文的神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練。

    1.3 L2B數(shù)據(jù)與ERA5數(shù)據(jù)對比分析

    本文選取2019年發(fā)生在西太平洋的臺風(fēng)“海貝思”為研究對象,首先對ERA5風(fēng)速和L2B風(fēng)速進行匹配,匹配時,在空間上選用的是雙線性插值方法,在時間上選用的是三節(jié)點插值方法?!昂X愃肌弊顝姇r期為10月7日至10月10日,通過獲取這期間每天上午經(jīng)過西太平洋地區(qū)的L2B風(fēng)速,然后與ERA5風(fēng)速進行時空插值,一共得到36 482組數(shù)據(jù)。

    對L2B風(fēng)速與ERA5風(fēng)速的相關(guān)性進行分風(fēng)速段分析(圖1a~圖1d)可知:當(dāng)vERA5<4 m/s時,L2B風(fēng)速與ERA5風(fēng)速的均方根誤差(Root Mean Square Error,RMSE)為1.98 m/s,相關(guān)系數(shù)(R)為0.27,雖然RMSE滿足HY-2B散射計的技術(shù)指標(biāo)(2 m/s),但相關(guān)性很差。當(dāng)vERA5>24 m/s時,L2B風(fēng)速與ERA5風(fēng)速的RMSE為4.15 m/s,R為0.17,在相關(guān)性很差的同時,誤差也超過了HY-2B散射計的技術(shù)指標(biāo)。造成這種情況的主要原因是在低風(fēng)速條件下(2~4 m/s),海面較為平靜,對電磁波的反射主要為鏡面反射,其后向散射能量較小,從而導(dǎo)致儀器測量信噪比較低,后向散射系數(shù)測量誤差較大;在高風(fēng)速段,海面后向散射系數(shù)隨風(fēng)速的變化越來越小,并逐漸趨于飽和[15]。因此,散射計風(fēng)速反演被限制在了一定的范圍內(nèi),HY-2B散射計的風(fēng)速有效觀測范圍為4~24 m/s。當(dāng)風(fēng)速為中低風(fēng)速時(4 m/s≤vERA5<15 m/s),L2B風(fēng)速與ERA5風(fēng)速的RMSE為1.44 m/s,R為0.90,當(dāng)風(fēng)速為中高風(fēng)速時(15 m/s≤vERA5≤24 m/s),L2B風(fēng)速與ERA5風(fēng)速的RMSE為1.84 m/s,R為0.74。中低風(fēng)速和中高風(fēng)速都滿足HY-2B散射計的技術(shù)指標(biāo),其中低風(fēng)速的RMSE比中高風(fēng)速小,且R大于中高風(fēng)速。

    圖1 L2B風(fēng)速與ERA5風(fēng)速相關(guān)分析Fig.1 Correlation analysis between wind speed of L2B and ERA5 wind speed

    對L2B風(fēng)速與ERA5風(fēng)速的相關(guān)性進行總體分析(圖1e)可知:在中低風(fēng)速階段,L2B風(fēng)速比ERA5風(fēng)速偏高,在中高風(fēng)速階段,L2B風(fēng)速總體比ERA5風(fēng)速偏低。此現(xiàn)象主要由降水所致。一般情況下,降水會使后向散射系數(shù)增大,從而使反演的風(fēng)速偏大;而在大風(fēng)情況下(如臺風(fēng)),衰減影響又起主導(dǎo)作用,導(dǎo)致反演的風(fēng)速偏小[12]。

    HY-2B衛(wèi)星散射計的L2B風(fēng)場數(shù)據(jù)由最大似然估計方法反演所得,該算法需要以某個已知的GMF為基礎(chǔ),L2B風(fēng)場數(shù)據(jù)采用的是NSCAT-4 GMF。NSCAT-4來源于NSCAT-2,它提高了浮標(biāo)與反演Ku波段散射風(fēng)的一致性,特別是對于風(fēng)速大于15 m/s的風(fēng)。但由于當(dāng)時可用的高風(fēng)速數(shù)據(jù)有限,浮標(biāo)觀測風(fēng)速和ECWMF風(fēng)速也會存在一定的誤差,導(dǎo)致NSCAT-4 GMF反演得到的風(fēng)速誤差在中高風(fēng)速階段大于中低風(fēng)速階段,這與圖1b和圖1c的分析結(jié)果一致。

    通過對L2B風(fēng)速與ERA5風(fēng)速的對比分析,可以發(fā)現(xiàn),在不同風(fēng)速段,L2B風(fēng)速的精度是不同的,并且影響風(fēng)速精度的因素也有所不同。為了有效地避免不同因素之間的相互影響,減少非線性擬合的復(fù)雜程度,可以分風(fēng)速段進行研究。在以往的研究中,可用的高風(fēng)速數(shù)據(jù)有限,導(dǎo)致NSCAT-4 GMF對中高風(fēng)速的反演精度依舊有待提高,為此,可以考慮使用其他方法進行風(fēng)速反演。

    2 研究方法

    2.1 風(fēng)速反演模型方法

    人工神經(jīng)網(wǎng)絡(luò)簡稱神經(jīng)網(wǎng)絡(luò),是機器學(xué)習(xí)諸多算法中的一種,能夠處理高度非線性化問題。本文選取反向傳播神經(jīng)網(wǎng)絡(luò)對訓(xùn)練集進行訓(xùn)練,進而建立風(fēng)速反演模型。

    反向傳播神經(jīng)網(wǎng)絡(luò)在從輸入層到輸出層的過程中,通過線性求和函數(shù)和激活函數(shù)建立層與層之間的聯(lián)系,不同層之間激活函數(shù)可以不同。多特征多樣本量的線性求和函數(shù)為

    式中,Z l為l層的線性輸出;l為層數(shù);W l為(l-1)層到l層之間各節(jié)點的權(quán)重;A l-1為(l-1)層到l層的輸入值;b l為(l-1)層到l層的偏置。

    激活函數(shù)是對線性輸出后進行非線性運算,將非線性特性引入網(wǎng)絡(luò)中。常用的激活函數(shù)有Sigmoid函數(shù)、tanh函數(shù)和Re LU函數(shù),3個函數(shù)的公式為

    經(jīng)過激活函數(shù),Z l轉(zhuǎn)換為A l后再作為下一層的輸入值,進行下一層的線性函數(shù)和激活函數(shù)處理,直到輸出層輸出結(jié)果。得到向前傳播結(jié)果后,在進行反向傳播前,利用成本函數(shù)計算網(wǎng)絡(luò)輸出值與真實值之間總體誤差。在回歸問題中,常采用均方誤差(Mean Squared Error,MSE)作為神經(jīng)網(wǎng)絡(luò)的成本函數(shù),成本函數(shù)的計算公式為

    式中,A l-last為神經(jīng)網(wǎng)絡(luò)的輸出值,y∧為真實值,N為樣本個數(shù)。

    確定成本函數(shù)后,利用誤差反向傳播算法對網(wǎng)絡(luò)的權(quán)重和偏置(統(tǒng)稱為參數(shù)θ)進行更新,誤差反向傳播算法公式為

    式中,θ為更新前參數(shù),θ′為更新后參數(shù),η為學(xué)習(xí)率,?θ為梯度算子。

    通過前向傳播和誤差反向傳播的交替進行,不斷更新各個節(jié)點的權(quán)重和偏置,使得成本函數(shù)最小化:

    反向傳播神經(jīng)網(wǎng)絡(luò)算法通過找到使成本函數(shù)最小化的參數(shù)θ,從而建立使輸出結(jié)果與真實結(jié)果誤差最小的網(wǎng)絡(luò)模型,以達到非線性擬合的目的。

    2.2 模型檢驗指標(biāo)

    采用均方誤差(MSE)、平均絕對值誤差(MAE)和相關(guān)性系數(shù)(R)等統(tǒng)計指標(biāo)對模型的性能進行評估。MSE,MAE和R的具體計算公式如下:

    式中,N為樣本量為第i個模型輸出值,y i為第i個真實參考值為模型輸出值的平均值為真實參考值的平均值。

    3 網(wǎng)絡(luò)訓(xùn)練

    3.1 建立網(wǎng)絡(luò)框架

    反向傳播神經(jīng)網(wǎng)絡(luò)一般由輸入層、隱藏層和輸出層組成。其中輸入層一共有20個節(jié)點,每個風(fēng)矢量單元格包括4組數(shù)據(jù),每組數(shù)據(jù)包括后向散射系數(shù)(σ)、電磁波入射角的正余弦值(sinλ,cosλ)、電磁波方位角與風(fēng)向的夾角的正余弦值(sinχ,cosχ)這5個特征。輸出層為1個節(jié)點,即風(fēng)速值,如圖2所示。

    圖2 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)Fig.2 Structure of neural network

    神經(jīng)網(wǎng)絡(luò)中的參數(shù)設(shè)置依據(jù)隱藏層數(shù)、隱藏層節(jié)點數(shù)、學(xué)習(xí)率、mini-batch大小和激活函數(shù)。

    1)隱藏層數(shù)。當(dāng)隱藏層數(shù)為2層,通過激活函數(shù)可以表示各種精度的決策邊界,并且可以擬合各種精度的平滑映射。

    2)隱藏層節(jié)點數(shù)。隱藏層節(jié)點數(shù)的確定采用了經(jīng)驗公式[16],該經(jīng)驗公式為

    式中:Ni是輸入層神經(jīng)元個數(shù);No是輸出層神經(jīng)元個數(shù);Ns訓(xùn)練集的樣本數(shù);α為可調(diào)參數(shù),通常取2~10。經(jīng)過調(diào)試,發(fā)現(xiàn)當(dāng)α=8時,即當(dāng)節(jié)點數(shù)為173時,訓(xùn)練效果最好。

    3)學(xué)習(xí)率。通過對不同學(xué)習(xí)率的網(wǎng)絡(luò)進行訓(xùn)練,發(fā)現(xiàn)學(xué)習(xí)率為0.01時,模型能夠收斂到最小值附近,且收斂的速度較快。

    4)mini-batch大小。當(dāng)mini-batch大小是2n,代碼運行較快。

    5)激活函數(shù)。隱藏層采用Sigmoid函數(shù),輸出層采用ReLU函數(shù)。Re LU函數(shù)使得模型的輸出值大于0(風(fēng)速值為正),避免了梯度爆炸和梯度消失問題,防止了參數(shù)無法更新,并加快了網(wǎng)絡(luò)訓(xùn)練速度。

    3.2 數(shù)據(jù)集的劃分

    當(dāng)5 m/s<vERA5<15 m/s時,海面對電磁波的后向散射和海面風(fēng)場之間有很好的對應(yīng)關(guān)系[17]。當(dāng)風(fēng)速超過15 m/s時,可能產(chǎn)生波浪破裂并與海面氣流分離等現(xiàn)象,對海面粗糙度造成影響,從而使影響海面后向散射系數(shù)的因素多樣化[18]。因此,根據(jù)HY-2B散射計的有效風(fēng)速范圍(4~24 m/s),本文將進行3次不同的神經(jīng)網(wǎng)絡(luò)訓(xùn)練,分別為全風(fēng)速神經(jīng)網(wǎng)絡(luò)訓(xùn)練(4 m/s≤vERA5≤24 m/s)、中低風(fēng)速神經(jīng)網(wǎng)絡(luò)訓(xùn)練(4 m/s≤vERA5<15 m/s)以及中高風(fēng)速神經(jīng)網(wǎng)絡(luò)訓(xùn)練(15 m/s≤vERA5≤24 m/s)。全風(fēng)速共有29 199組數(shù)據(jù)、中低風(fēng)速共有25 739組數(shù)據(jù),中高風(fēng)速共有3 460組數(shù)據(jù)。訓(xùn)練集和測試集以8∶2的比例進行隨機抽取。在訓(xùn)練過程中,訓(xùn)練數(shù)據(jù)每遍歷網(wǎng)絡(luò)一次就從訓(xùn)練集中隨機抽取20%數(shù)據(jù)作為驗證集進行驗證。

    本文使用核密度估計[19]方法來查看數(shù)據(jù)的分布情況。核密度估計是非參數(shù)估計中的一種,相較于參數(shù)估計,它不加入任何先驗知識,而是根據(jù)數(shù)據(jù)本身的特點、性質(zhì)來擬合分布。由圖3(縱坐標(biāo)為觀測數(shù)據(jù)落在某一風(fēng)速區(qū)間的概率值,橫坐標(biāo)為風(fēng)速區(qū)間)可知,3種神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練集和測試集的核密度分布相似,符合訓(xùn)練要求。

    為了保證網(wǎng)絡(luò)的收斂效果,在進行訓(xùn)練前,對數(shù)據(jù)進行標(biāo)準(zhǔn)化。標(biāo)準(zhǔn)化可以在不清楚各個維度的相對重要程度之前,讓輸入的各個維度分布相近。本文使用Z-score標(biāo)準(zhǔn)化[20]對訓(xùn)練集進行處理。

    圖3 不同風(fēng)速范圍下訓(xùn)練集和測試集的核密度估計Fig.3 Nuclear density estimation of training set and test set under different range of wind speeds

    3.3 訓(xùn)練結(jié)果

    進行網(wǎng)絡(luò)框架的構(gòu)建、參數(shù)設(shè)置,以及數(shù)據(jù)集的預(yù)處理后,開始對網(wǎng)絡(luò)進行訓(xùn)練。隨著訓(xùn)練數(shù)據(jù)遍歷網(wǎng)絡(luò)次數(shù)(Epoch)的增加,模型輸出值與真實值的總體誤差將逐漸減少。在模型訓(xùn)練過程中,如果連續(xù)30次驗證集的均方誤差(MSE)不再下降,則停止訓(xùn)練。由各風(fēng)速模型訓(xùn)練時絕對值誤差與均方誤差的走勢(圖4)可知,隨著訓(xùn)練次數(shù)的增加,訓(xùn)練集和驗證集的均方誤差和絕對值誤差逐漸下降,并最終趨于穩(wěn)定。

    圖4 各風(fēng)速模型訓(xùn)練時絕對值誤差與均方誤差的走勢Fig.4 The trend of absolute value error and mean square error during training of different wind speed model

    訓(xùn)練結(jié)束后,采用平均絕對誤差和均方誤差評估模型性能。在中高風(fēng)速中,L2B風(fēng)速與ERA5風(fēng)速的平均絕對值誤差(MAE)為1.29 m/s,均方誤差(MSE)為3.25 m/s;模型反演風(fēng)速與ERA5風(fēng)速的MAE為0.27 m/s,MSE為0.18 m/s。在中低風(fēng)速中,L2B風(fēng)速與ERA5風(fēng)速的MAE為0.95 m/s,MSE為2.10 m/s;模型反演風(fēng)速與ERA5風(fēng)速的MAE為0.24 m/s,MSE為0.14 m/s。在全風(fēng)速中,L2B風(fēng) 速 與ERA5風(fēng) 速MAE為0.99 m/s,MSE為2.22 m/s;模型反演風(fēng)速與ERA5風(fēng)速的MAE為0.34 m/s,MSE為0.32 m/s,如表1所示。

    表1 各模型訓(xùn)練集反演風(fēng)速誤差情況Table 1 Wind speed inversion errors from training set different each models

    由L2B風(fēng)速和各模型訓(xùn)練集反演風(fēng)速與ERA5風(fēng)速的相關(guān)分析(圖5)可知,在中高風(fēng)速中,L2B風(fēng)速與ERA5風(fēng)速相關(guān)系數(shù)(R)為0.75;模型反演風(fēng)速與ERA5風(fēng)速R為0.99。在中低風(fēng)速中,L2B風(fēng)速與ERA5風(fēng)速R為0.89;模型反演風(fēng)速與ERA5風(fēng)速R為0.99。在全風(fēng)速中,L2B風(fēng)速與ERA5風(fēng)速R為0.94,模型反演風(fēng)速與ERA5風(fēng)速R為0.99。

    從模型的訓(xùn)練效果來看,不同模型反演得到的風(fēng)速與ERA5風(fēng)速相關(guān)系數(shù)(R)都達到0.99。全風(fēng)速模型的MSE最高(0.32 m/s),其次是中高風(fēng)速模型(0.18 m/s),中低風(fēng)速模型最小(0.14 m/s)。由此可見,對于在不同風(fēng)速段下建立的神經(jīng)網(wǎng)絡(luò)模型,網(wǎng)絡(luò)性能的表現(xiàn)有所不同。造成這種現(xiàn)象的原因可能是:在不同的風(fēng)速段下,影響后向散射系數(shù)與風(fēng)速之間的映射關(guān)系的因素有所不同,將不同風(fēng)速段放在一起訓(xùn)練時,后向散射系數(shù)與風(fēng)速之間的映射關(guān)系變得更加復(fù)雜。中低風(fēng)速模型的MSE和MAE比中高風(fēng)模型風(fēng)速小,這可能是由中高風(fēng)速模型包含了更多的臺風(fēng)信息、擬合情況更加復(fù)雜所致。

    圖5 在訓(xùn)練集中v L2B與v ERA5的相關(guān)分析以及各模型反演風(fēng)速與v ERA5的相關(guān)分析Fig.5 The correlation analysis between v L2B and v ERA5,as well as the correlation analysis between wind speed inversion resutts of each model and v ERA5 in the training set

    3.4 測試結(jié)果

    為進一步驗證模型的有效性,利用測試集對模型的泛化能力進行測試。結(jié)果(表2)表明,在中高風(fēng)速段,L2B風(fēng)速與ERA5風(fēng)速的MAE為1.38 m/s,MSE為3.99m/s;模型反演風(fēng)速與ERA5風(fēng)速的MAE為0.48 m/s,MSE為0.54 m/s。在中低風(fēng)速段,L2B風(fēng)速與ERA5風(fēng)速MAE為0.95 m/s,MSE為2.00 m/s;模型反演風(fēng)速與ERA5風(fēng)速的MAE為0.35 m/s,MSE為0.27 m/s。在全風(fēng)速段,L2B風(fēng) 速 與ERA5風(fēng) 速MAE為1.02 m/s,MSE為2.29 m/s;模型反演風(fēng)速與ERA5風(fēng)速的MAE為0.42 m/s,MSE為0.46 m/s。

    表2 各模型測試集反演風(fēng)速誤差情況Table 2 Wind speed inversion errors from test set of different models

    由L2B風(fēng)速和各模型測試集反演風(fēng)速與ERA5風(fēng)速的相關(guān)分析(圖6)可知:在中高風(fēng)速中,L2B風(fēng)速與ERA5風(fēng)速相關(guān)系數(shù)(R)為0.71;模型反演風(fēng)速與ERA5風(fēng)速R為0.95。在中低風(fēng)速中,L2B風(fēng)速與ERA5風(fēng)速R為0.90;模型反演風(fēng)速與ERA5風(fēng)速R為0.98。在全風(fēng)速中,L2B風(fēng)速與ERA5風(fēng)速R為0.94,模型反演風(fēng)速與ERA5風(fēng)速R為0.99。全風(fēng)速和中低風(fēng)速模型反演的風(fēng)速與ERA5風(fēng)速的R都超過了0.98,而中高風(fēng)速模型的R從訓(xùn)練集時的0.99下降到了0.95。在MSE方面,中高風(fēng)速模型MSE最高(0.54 m/s),其次是全風(fēng)速模型(0.46 m/s),中低風(fēng)速模型最低(0.27 m/s)。在測試集中,全風(fēng)速模型的MSE低于中高風(fēng)速模型,與訓(xùn)練集的情況相反;中低風(fēng)速模型的MSE依舊最低;綜合相關(guān)系數(shù)和均方誤差來看,高風(fēng)速模型的泛化能力弱于其他2種模型,這可能由高風(fēng)速下的訓(xùn)練和測試樣本較少所致。整體來看,3種網(wǎng)絡(luò)模型都有很好的泛化能力。

    圖6 在測試集中v L2B與v ERA5的相關(guān)分析以及各模型反演風(fēng)速與v ERA5的相關(guān)分析Fig.6 The correlation analysis between v L2B and v ERA5,as well as the correlation analysis between wind speed inverison results of each model and v ERA5 in the test set

    4 結(jié) 論

    針對西太平洋地區(qū)有臺風(fēng)存在的情況,以ERA5風(fēng)速為真實風(fēng)速參考值,基于HY-2B衛(wèi)星散射計數(shù)據(jù),利用神經(jīng)網(wǎng)絡(luò)方法進行風(fēng)速反演。受物理機制的制約,在vERA5<4 m/s和vERA5>24 m/s情況下,利用后向散射進行風(fēng)速反演不再有效。作為神經(jīng)網(wǎng)絡(luò)訓(xùn)練時的輸入特征,過低或過高風(fēng)速下的后向散射值會對模型訓(xùn)練造成不利影響,因此選擇在有效風(fēng)速內(nèi)的后向散射值進行模型訓(xùn)練。由于不同風(fēng)速段內(nèi)影響風(fēng)速反演的環(huán)境因子不同,因此分別建立全風(fēng)速(4 m/s≤vEAR5≤24 m/s)、中低風(fēng)速(4 m/s≤vEAR5<15 m/s)、以及中高風(fēng)速(15 m/s≤vEAR5≤24 m/s)三種不同的風(fēng)速反演模型。通過對模型反演結(jié)果的對比分析,得到結(jié)論:

    1)3種風(fēng)速反演模型的MSE和MAE優(yōu)于基于NSCAT-4 GMF得到的L2B風(fēng)速的MSE和AME,訓(xùn)練結(jié)果和測試結(jié)果都有良好的表現(xiàn),其中中低風(fēng)速模型的MSE和MAE最低,中高風(fēng)速模型的MSE和MAE降低幅度最大。

    2)對比訓(xùn)練集和測試集,雖然3種反演模型在測試集中的MSE和MAE都有所升高,但依舊優(yōu)于L2B風(fēng)速。

    3)相比于全風(fēng)速,對風(fēng)速進行分區(qū)間所建立的反演模型效果更好,雖然中高風(fēng)速模型在測試集中的誤差大于全風(fēng)速模型,但這可能是由中高速模型測試集數(shù)據(jù)較少所致。

    相比于L2B風(fēng)速,利用本文模型反演的風(fēng)速MSE和MAE都有所下降,同時相關(guān)系數(shù)(R)有所增加,主要得益于3種因素:神經(jīng)網(wǎng)絡(luò)強大的非線性擬合以及泛化能力;通過分區(qū)間進行反演,避免了影響不同區(qū)間后向散射因素的相互影響;將ERA5作為真實風(fēng)速參考值,具有較高精度和較高時空分辨率的特點,能夠更好地模擬真實風(fēng)速。然而,由于中高風(fēng)速模型的訓(xùn)練數(shù)據(jù)較少,因此泛化能力相比于其他2種模型較弱,在今后的工作中,我們將增加中高風(fēng)速的樣本數(shù)量以優(yōu)化中高風(fēng)速模型。

    猜你喜歡
    散射系數(shù)風(fēng)場反演
    等離子體層嘶聲波對輻射帶電子投擲角散射系數(shù)的多維建模*
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    基于FLUENT的下?lián)舯┝魅S風(fēng)場建模
    北部灣后向散射系數(shù)的時空分布與變化分析
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    “最美風(fēng)場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    側(cè)向風(fēng)場中無人機的飛行研究
    疊前同步反演在港中油田的應(yīng)用
    一維帶限Weierstrass分形粗糙面電磁散射的微擾法研究
    欧美色欧美亚洲另类二区| 午夜福利18| 亚洲,欧美精品.| 精品国产乱码久久久久久男人| 日本一区二区免费在线视频| 午夜精品一区二区三区免费看| www日本在线高清视频| 欧美+亚洲+日韩+国产| 搡老岳熟女国产| 99国产综合亚洲精品| 一级作爱视频免费观看| 久久人妻av系列| 亚洲一区中文字幕在线| 日韩 欧美 亚洲 中文字幕| 欧美性猛交黑人性爽| 亚洲一卡2卡3卡4卡5卡精品中文| 国产人伦9x9x在线观看| 日韩中文字幕欧美一区二区| 亚洲自拍偷在线| 久久久精品欧美日韩精品| 久久午夜亚洲精品久久| 日韩欧美国产在线观看| 一二三四在线观看免费中文在| 巨乳人妻的诱惑在线观看| 淫秽高清视频在线观看| 日韩av在线大香蕉| 一级黄色大片毛片| 国产高清videossex| 欧美日韩一级在线毛片| 日韩大码丰满熟妇| av福利片在线| 狂野欧美白嫩少妇大欣赏| 久久国产精品人妻蜜桃| 欧美 亚洲 国产 日韩一| 日本五十路高清| 一二三四社区在线视频社区8| 欧美在线一区亚洲| 亚洲精品国产一区二区精华液| 亚洲最大成人中文| 国产真人三级小视频在线观看| a级毛片a级免费在线| bbb黄色大片| 日日爽夜夜爽网站| www日本在线高清视频| 变态另类成人亚洲欧美熟女| 国产精品1区2区在线观看.| 亚洲电影在线观看av| 日韩欧美 国产精品| 日本黄色视频三级网站网址| 老熟妇仑乱视频hdxx| 国产爱豆传媒在线观看 | 99热只有精品国产| 每晚都被弄得嗷嗷叫到高潮| 在线观看免费视频日本深夜| 国语自产精品视频在线第100页| 天堂av国产一区二区熟女人妻 | 999久久久国产精品视频| 午夜视频精品福利| 久久香蕉精品热| 动漫黄色视频在线观看| 国产精品免费视频内射| 人人妻人人澡欧美一区二区| 18禁美女被吸乳视频| 免费在线观看亚洲国产| 成人永久免费在线观看视频| 亚洲欧美日韩无卡精品| 亚洲国产欧美一区二区综合| 午夜精品一区二区三区免费看| 妹子高潮喷水视频| 19禁男女啪啪无遮挡网站| 成人av在线播放网站| 国产99久久九九免费精品| 两个人的视频大全免费| or卡值多少钱| 精品日产1卡2卡| 精品福利观看| 一区福利在线观看| 又大又爽又粗| 国产激情欧美一区二区| 国产高清视频在线播放一区| 欧美性猛交黑人性爽| 老司机午夜十八禁免费视频| 国产精品野战在线观看| 男人舔奶头视频| 久久久精品欧美日韩精品| 97超级碰碰碰精品色视频在线观看| 中文资源天堂在线| 国产在线观看jvid| 一个人观看的视频www高清免费观看 | 中国美女看黄片| 欧洲精品卡2卡3卡4卡5卡区| 小说图片视频综合网站| 床上黄色一级片| 女警被强在线播放| 窝窝影院91人妻| 身体一侧抽搐| 天堂√8在线中文| 麻豆国产av国片精品| 亚洲在线自拍视频| 欧美+亚洲+日韩+国产| 亚洲av日韩精品久久久久久密| 国产欧美日韩一区二区精品| 亚洲av中文字字幕乱码综合| 美女大奶头视频| 69av精品久久久久久| 婷婷丁香在线五月| 麻豆成人av在线观看| 国产蜜桃级精品一区二区三区| 国产成人影院久久av| 欧美大码av| 熟女电影av网| 99精品欧美一区二区三区四区| 大型av网站在线播放| 99精品久久久久人妻精品| 国产99久久九九免费精品| 两个人的视频大全免费| a级毛片a级免费在线| e午夜精品久久久久久久| 欧美性长视频在线观看| 草草在线视频免费看| www.www免费av| 精品一区二区三区av网在线观看| 色哟哟哟哟哟哟| 久久久久久久午夜电影| 99热这里只有精品一区 | 不卡一级毛片| 国产亚洲欧美在线一区二区| av片东京热男人的天堂| 精品不卡国产一区二区三区| 亚洲无线在线观看| svipshipincom国产片| 欧美三级亚洲精品| 在线观看美女被高潮喷水网站 | 亚洲真实伦在线观看| 男女之事视频高清在线观看| 熟女少妇亚洲综合色aaa.| 变态另类成人亚洲欧美熟女| 成人国产一区最新在线观看| 九色成人免费人妻av| 操出白浆在线播放| 国产精品99久久99久久久不卡| 国产精品香港三级国产av潘金莲| 久久午夜亚洲精品久久| 中文字幕久久专区| 欧美av亚洲av综合av国产av| 这个男人来自地球电影免费观看| 怎么达到女性高潮| 12—13女人毛片做爰片一| 国产精品久久久久久精品电影| 久久久久九九精品影院| 免费一级毛片在线播放高清视频| 亚洲国产精品sss在线观看| 亚洲黑人精品在线| 一卡2卡三卡四卡精品乱码亚洲| 久久久久国产精品人妻aⅴ院| 1024香蕉在线观看| 欧美丝袜亚洲另类 | www日本黄色视频网| 亚洲五月天丁香| 国产精品 欧美亚洲| 身体一侧抽搐| 精品国产乱码久久久久久男人| 国内久久婷婷六月综合欲色啪| 夜夜看夜夜爽夜夜摸| 午夜福利免费观看在线| 国产人伦9x9x在线观看| 精品久久久久久,| 免费搜索国产男女视频| 久久久久久久久免费视频了| 巨乳人妻的诱惑在线观看| 一本久久中文字幕| 亚洲欧美日韩东京热| 香蕉av资源在线| 热99re8久久精品国产| 国产亚洲精品一区二区www| 国产av一区二区精品久久| 好男人在线观看高清免费视频| 欧美乱码精品一区二区三区| www.999成人在线观看| 欧美成狂野欧美在线观看| 亚洲欧美日韩高清在线视频| 亚洲国产精品999在线| 日本一二三区视频观看| 日韩国内少妇激情av| 国产精品一及| 99久久综合精品五月天人人| 又黄又爽又免费观看的视频| 一进一出好大好爽视频| 午夜福利在线观看吧| 国产又黄又爽又无遮挡在线| 三级国产精品欧美在线观看 | 日本撒尿小便嘘嘘汇集6| 超碰成人久久| 国产黄a三级三级三级人| 成年免费大片在线观看| 12—13女人毛片做爰片一| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久| 亚洲一区二区三区色噜噜| 丝袜美腿诱惑在线| 可以在线观看毛片的网站| 亚洲欧美日韩无卡精品| 午夜精品在线福利| 国产成+人综合+亚洲专区| 久久这里只有精品中国| 可以在线观看的亚洲视频| 手机成人av网站| 亚洲,欧美精品.| 国产成人一区二区三区免费视频网站| 国语自产精品视频在线第100页| 久久香蕉国产精品| 欧美成人免费av一区二区三区| 动漫黄色视频在线观看| 岛国视频午夜一区免费看| 国产单亲对白刺激| 亚洲一区二区三区色噜噜| 51午夜福利影视在线观看| 日韩大尺度精品在线看网址| 国产精品久久久人人做人人爽| x7x7x7水蜜桃| 亚洲 欧美 日韩 在线 免费| 免费搜索国产男女视频| 国产伦在线观看视频一区| 美女黄网站色视频| 国产精品久久电影中文字幕| 国产野战对白在线观看| 成人av在线播放网站| 免费无遮挡裸体视频| 欧美一区二区国产精品久久精品 | 国产精品久久久av美女十八| tocl精华| 激情在线观看视频在线高清| 精品久久久久久久毛片微露脸| 狂野欧美激情性xxxx| 亚洲精品中文字幕在线视频| 欧美zozozo另类| 亚洲 欧美一区二区三区| 男女下面进入的视频免费午夜| 亚洲人成电影免费在线| 国产精品免费视频内射| 久久中文字幕一级| 久久久精品大字幕| 91麻豆精品激情在线观看国产| 欧美丝袜亚洲另类 | 99在线人妻在线中文字幕| 欧美黄色淫秽网站| 又黄又粗又硬又大视频| 日本黄色视频三级网站网址| 丁香六月欧美| 国语自产精品视频在线第100页| 2021天堂中文幕一二区在线观| 亚洲国产精品成人综合色| 观看免费一级毛片| 亚洲一区二区三区色噜噜| 国产精品爽爽va在线观看网站| 午夜福利成人在线免费观看| 婷婷精品国产亚洲av| 在线永久观看黄色视频| 91成年电影在线观看| 日本一二三区视频观看| 一级毛片精品| 亚洲国产高清在线一区二区三| 在线永久观看黄色视频| 日本 欧美在线| 国内精品久久久久精免费| 亚洲av成人不卡在线观看播放网| 中文字幕人妻丝袜一区二区| 久久中文字幕一级| 成人国语在线视频| 亚洲午夜精品一区,二区,三区| 国产爱豆传媒在线观看 | 午夜激情福利司机影院| 啦啦啦韩国在线观看视频| 亚洲中文字幕日韩| 亚洲国产看品久久| 亚洲欧美一区二区三区黑人| 夜夜夜夜夜久久久久| 欧美乱色亚洲激情| 国产男靠女视频免费网站| 少妇人妻一区二区三区视频| 丝袜美腿诱惑在线| 一本精品99久久精品77| 手机成人av网站| 精品久久久久久,| 日本免费a在线| 香蕉国产在线看| 欧美绝顶高潮抽搐喷水| 神马国产精品三级电影在线观看 | xxx96com| 成人国产一区最新在线观看| 天堂动漫精品| 久久久国产成人免费| 免费在线观看日本一区| 99久久久亚洲精品蜜臀av| 一个人免费在线观看电影 | 日韩欧美国产在线观看| 一级黄色大片毛片| 欧美黑人精品巨大| a在线观看视频网站| 夜夜躁狠狠躁天天躁| 国产99白浆流出| 免费高清视频大片| 欧美成人一区二区免费高清观看 | 国产激情久久老熟女| 99国产精品一区二区三区| 欧美午夜高清在线| 久久久久国内视频| 日韩欧美在线二视频| 国产午夜福利久久久久久| 国产精品九九99| 欧美人与性动交α欧美精品济南到| 国产伦一二天堂av在线观看| 欧美黑人欧美精品刺激| 国产麻豆成人av免费视频| 亚洲色图av天堂| 免费av毛片视频| 操出白浆在线播放| 欧美性长视频在线观看| 男女下面进入的视频免费午夜| 亚洲国产欧美一区二区综合| www.精华液| 国产视频一区二区在线看| 亚洲激情在线av| 久久精品91无色码中文字幕| 伦理电影免费视频| 老汉色∧v一级毛片| 国产精品精品国产色婷婷| 最新美女视频免费是黄的| 亚洲电影在线观看av| 亚洲狠狠婷婷综合久久图片| 精品欧美一区二区三区在线| 亚洲最大成人中文| 亚洲18禁久久av| 免费看日本二区| 国产亚洲av高清不卡| 国产三级黄色录像| 久久久久久亚洲精品国产蜜桃av| 国产乱人伦免费视频| 国产精品久久久久久精品电影| 精品欧美国产一区二区三| 男男h啪啪无遮挡| av欧美777| cao死你这个sao货| 高清在线国产一区| 日本成人三级电影网站| 亚洲av日韩精品久久久久久密| 99riav亚洲国产免费| 欧美最黄视频在线播放免费| 国产真人三级小视频在线观看| 老鸭窝网址在线观看| 免费观看精品视频网站| 在线观看一区二区三区| 黄片大片在线免费观看| 久久午夜综合久久蜜桃| 国产精品综合久久久久久久免费| 成在线人永久免费视频| 国产激情久久老熟女| 国产视频一区二区在线看| a级毛片a级免费在线| 国产高清videossex| 国产爱豆传媒在线观看 | 好看av亚洲va欧美ⅴa在| 免费看a级黄色片| 亚洲国产精品成人综合色| 精品久久久久久久人妻蜜臀av| av在线播放免费不卡| 麻豆av在线久日| 欧美日韩国产亚洲二区| 日本精品一区二区三区蜜桃| 亚洲人成电影免费在线| 国产亚洲精品综合一区在线观看 | 黄色片一级片一级黄色片| 欧美乱码精品一区二区三区| 久久婷婷人人爽人人干人人爱| 亚洲人与动物交配视频| 欧美精品亚洲一区二区| 久久人妻福利社区极品人妻图片| 成人18禁在线播放| 国产精品综合久久久久久久免费| 亚洲性夜色夜夜综合| 精品日产1卡2卡| 欧美成人午夜精品| 天天添夜夜摸| 天堂影院成人在线观看| 亚洲乱码一区二区免费版| 欧美色视频一区免费| 曰老女人黄片| 国产亚洲欧美在线一区二区| 在线观看免费午夜福利视频| 国产视频一区二区在线看| 国产精品九九99| www.精华液| 亚洲欧美日韩无卡精品| 亚洲欧美日韩高清专用| 国产黄片美女视频| 精品国产乱码久久久久久男人| 国产精品久久电影中文字幕| 久久人人精品亚洲av| 在线观看免费日韩欧美大片| 在线十欧美十亚洲十日本专区| 村上凉子中文字幕在线| 一级毛片精品| 搡老岳熟女国产| 亚洲精华国产精华精| 99国产精品99久久久久| 国产又色又爽无遮挡免费看| 操出白浆在线播放| 亚洲欧美精品综合一区二区三区| av在线天堂中文字幕| 国产高清videossex| 日韩精品青青久久久久久| 午夜福利在线观看吧| 午夜日韩欧美国产| 十八禁网站免费在线| 欧美成人免费av一区二区三区| 精品人妻1区二区| 一区二区三区国产精品乱码| 精华霜和精华液先用哪个| 后天国语完整版免费观看| 成年版毛片免费区| 亚洲中文av在线| 久久国产精品人妻蜜桃| 色在线成人网| 毛片女人毛片| av免费在线观看网站| 嫩草影院精品99| 免费观看精品视频网站| 亚洲国产看品久久| 国产亚洲欧美98| 悠悠久久av| 一本久久中文字幕| 岛国在线观看网站| 色综合站精品国产| 久久久精品大字幕| 欧美黑人精品巨大| 首页视频小说图片口味搜索| 我的老师免费观看完整版| 一级片免费观看大全| 久久精品亚洲精品国产色婷小说| 免费高清视频大片| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一区av在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 一二三四在线观看免费中文在| 国内毛片毛片毛片毛片毛片| 日本一本二区三区精品| 亚洲精品久久成人aⅴ小说| 国产69精品久久久久777片 | 国产黄片美女视频| 国产伦人伦偷精品视频| 日韩欧美在线乱码| 女生性感内裤真人,穿戴方法视频| 欧美大码av| 97碰自拍视频| 香蕉国产在线看| 国产精品免费视频内射| 日本熟妇午夜| 五月玫瑰六月丁香| svipshipincom国产片| 少妇人妻一区二区三区视频| 精品一区二区三区四区五区乱码| 亚洲一区中文字幕在线| 制服丝袜大香蕉在线| 日韩欧美在线乱码| 亚洲熟妇中文字幕五十中出| 日韩欧美国产在线观看| 午夜a级毛片| 久久久久久亚洲精品国产蜜桃av| 亚洲人成网站在线播放欧美日韩| 日日干狠狠操夜夜爽| 国产午夜精品论理片| 天天添夜夜摸| 黄色毛片三级朝国网站| 88av欧美| 男女那种视频在线观看| 99久久精品热视频| 欧美日韩亚洲国产一区二区在线观看| 99国产极品粉嫩在线观看| 成人手机av| 久久香蕉精品热| 黄频高清免费视频| 中文资源天堂在线| 露出奶头的视频| 最近最新中文字幕大全电影3| 久久亚洲精品不卡| 欧美丝袜亚洲另类 | 免费一级毛片在线播放高清视频| 88av欧美| 久久人妻av系列| 啦啦啦免费观看视频1| 中文字幕精品亚洲无线码一区| ponron亚洲| 两性夫妻黄色片| 久久久久久免费高清国产稀缺| 青草久久国产| 国产精品久久久久久久电影 | 好男人电影高清在线观看| 久久九九热精品免费| 日韩精品中文字幕看吧| 性色av乱码一区二区三区2| 色综合婷婷激情| 老汉色av国产亚洲站长工具| 精品久久久久久久久久免费视频| 在线永久观看黄色视频| 老司机福利观看| 大型黄色视频在线免费观看| 一本一本综合久久| 精品国产亚洲在线| 国产精品久久视频播放| 日韩精品中文字幕看吧| av视频在线观看入口| 99国产精品一区二区三区| 亚洲性夜色夜夜综合| 精品高清国产在线一区| 免费高清视频大片| 一级毛片女人18水好多| 亚洲,欧美精品.| 91大片在线观看| 日韩欧美一区二区三区在线观看| 午夜福利成人在线免费观看| 色噜噜av男人的天堂激情| 婷婷丁香在线五月| 久久久国产欧美日韩av| 精品欧美一区二区三区在线| 国内少妇人妻偷人精品xxx网站 | 午夜精品久久久久久毛片777| 在线a可以看的网站| 久久精品成人免费网站| 亚洲一区中文字幕在线| 人成视频在线观看免费观看| 老熟妇乱子伦视频在线观看| 亚洲av第一区精品v没综合| 国产熟女xx| 我要搜黄色片| 国产精品乱码一区二三区的特点| 国产精品自产拍在线观看55亚洲| 香蕉久久夜色| 一二三四社区在线视频社区8| 国产v大片淫在线免费观看| 一a级毛片在线观看| 人妻久久中文字幕网| 伊人久久大香线蕉亚洲五| 精品熟女少妇八av免费久了| 成人特级黄色片久久久久久久| 日韩av在线大香蕉| 一二三四在线观看免费中文在| 日韩av在线大香蕉| 久热爱精品视频在线9| www.www免费av| 精品一区二区三区视频在线观看免费| 亚洲成av人片免费观看| 一进一出抽搐gif免费好疼| 18禁黄网站禁片午夜丰满| 国产精品久久久久久亚洲av鲁大| 中出人妻视频一区二区| 色噜噜av男人的天堂激情| 色尼玛亚洲综合影院| 国产主播在线观看一区二区| 日韩欧美在线二视频| 国产精品一及| 中文字幕久久专区| 色综合站精品国产| 熟妇人妻久久中文字幕3abv| 18禁美女被吸乳视频| 在线十欧美十亚洲十日本专区| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 亚洲一码二码三码区别大吗| 国产精品国产高清国产av| www.www免费av| 每晚都被弄得嗷嗷叫到高潮| 久9热在线精品视频| 欧美一区二区精品小视频在线| 狂野欧美白嫩少妇大欣赏| 欧美日本亚洲视频在线播放| 婷婷丁香在线五月| 在线观看一区二区三区| 日韩欧美国产在线观看| 亚洲真实伦在线观看| 成人亚洲精品av一区二区| 欧美午夜高清在线| 99久久综合精品五月天人人| 啦啦啦免费观看视频1| 久久久国产成人精品二区| 国产一区二区在线观看日韩 | 此物有八面人人有两片| 国产亚洲欧美98| 久久久国产欧美日韩av| 熟女少妇亚洲综合色aaa.| 99久久精品热视频| 怎么达到女性高潮| 亚洲精品粉嫩美女一区| 精品高清国产在线一区| av在线天堂中文字幕| 欧美不卡视频在线免费观看 | 一边摸一边做爽爽视频免费| 国产亚洲av嫩草精品影院| 亚洲专区中文字幕在线| 中文字幕精品亚洲无线码一区| 美女大奶头视频| 亚洲欧洲精品一区二区精品久久久| 亚洲,欧美精品.| 免费一级毛片在线播放高清视频| 一个人免费在线观看电影 | 色哟哟哟哟哟哟| 亚洲欧美激情综合另类| 亚洲欧美日韩无卡精品| 午夜成年电影在线免费观看| 在线永久观看黄色视频| 黑人欧美特级aaaaaa片| 精品久久久久久久毛片微露脸| 9191精品国产免费久久| 99久久无色码亚洲精品果冻| 国产午夜福利久久久久久| 99久久精品国产亚洲精品|