• <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+在线观看网站| 熟妇人妻不卡中文字幕| 丝袜脚勾引网站| 有码 亚洲区| 97超视频在线观看视频| 一区二区三区精品91| 大话2 男鬼变身卡| 日韩在线高清观看一区二区三区| 十八禁高潮呻吟视频| 狂野欧美激情性bbbbbb| 亚洲高清免费不卡视频| 搡老乐熟女国产| 久久精品国产自在天天线| 亚洲国产精品一区二区三区在线| 精品人妻在线不人妻| 精品亚洲成a人片在线观看| 免费黄网站久久成人精品| 亚洲精品第二区| 国产老妇伦熟女老妇高清| 日韩免费高清中文字幕av| 亚洲高清免费不卡视频| 国产黄片视频在线免费观看| 夫妻午夜视频| xxx大片免费视频| 国产成人免费观看mmmm| av福利片在线| 中文欧美无线码| 国产极品天堂在线| 国产免费福利视频在线观看| 高清午夜精品一区二区三区| 国产免费现黄频在线看| 国产又色又爽无遮挡免| 人妻 亚洲 视频| 少妇的逼水好多| 美女中出高潮动态图| 成人午夜精彩视频在线观看| 精品人妻在线不人妻| 日韩一区二区视频免费看| 日韩一区二区视频免费看| 国产精品国产av在线观看| 成年女人在线观看亚洲视频| 一级毛片aaaaaa免费看小| 亚洲国产毛片av蜜桃av| av在线app专区| 久久久国产欧美日韩av| 精品亚洲成a人片在线观看| 老女人水多毛片| 免费黄频网站在线观看国产| 新久久久久国产一级毛片| 国产亚洲av片在线观看秒播厂| 日本91视频免费播放| 97在线人人人人妻| 免费少妇av软件| 亚洲第一av免费看| 亚洲三级黄色毛片| 精品久久久噜噜| 中国美白少妇内射xxxbb| 在线播放无遮挡| 人人妻人人爽人人添夜夜欢视频| 精品少妇黑人巨大在线播放| 免费看光身美女| 大香蕉久久网| 色网站视频免费| 天堂中文最新版在线下载| 你懂的网址亚洲精品在线观看| 永久网站在线| 91久久精品电影网| 日韩一本色道免费dvd| 久久久午夜欧美精品| 日韩av在线免费看完整版不卡| av.在线天堂| 大陆偷拍与自拍| 99视频精品全部免费 在线| 在线天堂最新版资源| 91午夜精品亚洲一区二区三区| 久久精品成人免费网站| 中亚洲国语对白在线视频| 中文欧美无线码| 99精品久久久久人妻精品| 高清毛片免费观看视频网站 | 99精品久久久久人妻精品| tocl精华| 国产熟女午夜一区二区三区| 国产熟女午夜一区二区三区| 国产黄频视频在线观看| 男人操女人黄网站| 日本av免费视频播放| 亚洲色图 男人天堂 中文字幕| 久久午夜亚洲精品久久| 欧美人与性动交α欧美精品济南到| 18禁美女被吸乳视频| 少妇裸体淫交视频免费看高清 | 伊人久久大香线蕉亚洲五| 91字幕亚洲| 十八禁高潮呻吟视频| tube8黄色片| 一个人免费看片子| 国产又爽黄色视频| 两个人免费观看高清视频| 91麻豆精品激情在线观看国产 | 无人区码免费观看不卡 | 欧美成人午夜精品| 一级a爱视频在线免费观看| av免费在线观看网站| 国产一区二区三区视频了| 黑人巨大精品欧美一区二区mp4| 亚洲欧美日韩另类电影网站| 亚洲欧美激情在线| 999久久久精品免费观看国产| 色老头精品视频在线观看| 国产精品久久久久成人av| 老司机在亚洲福利影院| 97人妻天天添夜夜摸| 97人妻天天添夜夜摸| 黄网站色视频无遮挡免费观看| 999久久久国产精品视频| 国产真人三级小视频在线观看| 国产一区有黄有色的免费视频| 伊人久久大香线蕉亚洲五| 三上悠亚av全集在线观看| 亚洲欧美日韩另类电影网站| 欧美亚洲日本最大视频资源| 男人舔女人的私密视频| 亚洲成人免费av在线播放| 日韩视频在线欧美| 免费在线观看完整版高清| 捣出白浆h1v1| 久久中文字幕一级| 免费女性裸体啪啪无遮挡网站| 最近最新免费中文字幕在线| tocl精华| 老司机午夜十八禁免费视频| 免费久久久久久久精品成人欧美视频| 一级毛片精品| 久久国产精品人妻蜜桃| 最新美女视频免费是黄的| 我的亚洲天堂| 丝袜美足系列| 中文字幕精品免费在线观看视频| 免费日韩欧美在线观看| 欧美精品啪啪一区二区三区| 一级黄色大片毛片| www.自偷自拍.com| 啦啦啦免费观看视频1| 99re6热这里在线精品视频| 亚洲熟妇熟女久久| 在线看a的网站| 桃红色精品国产亚洲av| 99精品久久久久人妻精品| 在线av久久热| 午夜久久久在线观看| 搡老乐熟女国产| 9191精品国产免费久久| a级毛片黄视频| 久久久久久久国产电影| 国产成人系列免费观看| 蜜桃在线观看..| 亚洲九九香蕉| 亚洲熟女精品中文字幕| 老司机午夜十八禁免费视频| www.熟女人妻精品国产| 亚洲色图综合在线观看| 国产老妇伦熟女老妇高清| 黄色 视频免费看| 免费少妇av软件| 欧美变态另类bdsm刘玥| 亚洲熟妇熟女久久| 在线看a的网站| 18禁黄网站禁片午夜丰满| 大片电影免费在线观看免费| 脱女人内裤的视频| av电影中文网址| 亚洲三区欧美一区| 99久久人妻综合| 精品国产一区二区三区四区第35| 亚洲精品一二三| 99国产综合亚洲精品| av超薄肉色丝袜交足视频| 午夜老司机福利片| 黑人猛操日本美女一级片| 亚洲av片天天在线观看| 亚洲精品一卡2卡三卡4卡5卡| 精品亚洲乱码少妇综合久久| 国产91精品成人一区二区三区 | 国产日韩一区二区三区精品不卡| 精品国内亚洲2022精品成人 | 免费在线观看黄色视频的| 亚洲第一欧美日韩一区二区三区 | 国产高清视频在线播放一区| 国产在线视频一区二区| 欧美另类亚洲清纯唯美| 国产一区二区三区视频了| 精品国产乱码久久久久久男人| 正在播放国产对白刺激| 国产男靠女视频免费网站| 亚洲av国产av综合av卡| 中文字幕av电影在线播放| 色在线成人网| 99热国产这里只有精品6| 免费日韩欧美在线观看| 日本黄色日本黄色录像| 丰满人妻熟妇乱又伦精品不卡| 久久久久久人人人人人| 午夜福利在线免费观看网站| 黄色丝袜av网址大全| 日韩制服丝袜自拍偷拍| 亚洲国产欧美一区二区综合| 日日爽夜夜爽网站| 久热这里只有精品99| 在线观看免费日韩欧美大片| 亚洲色图av天堂| 19禁男女啪啪无遮挡网站| 久久久国产一区二区| 亚洲国产欧美日韩在线播放| 男女边摸边吃奶| 热re99久久国产66热| 国产野战对白在线观看| 黄片小视频在线播放| 老熟妇仑乱视频hdxx| 久久国产精品大桥未久av| 最近最新免费中文字幕在线| 午夜免费鲁丝| 黑人巨大精品欧美一区二区蜜桃| 91av网站免费观看| 久久午夜综合久久蜜桃| 国产精品99久久99久久久不卡| 久久亚洲精品不卡| 亚洲精品中文字幕一二三四区 | 97人妻天天添夜夜摸| 两个人免费观看高清视频| 男人操女人黄网站| 国产成人精品久久二区二区免费| 啦啦啦免费观看视频1| 欧美日韩成人在线一区二区| 午夜福利视频精品| 99久久人妻综合| 国产福利在线免费观看视频| 两个人看的免费小视频| 黑人巨大精品欧美一区二区蜜桃| 精品免费久久久久久久清纯 | 丰满人妻熟妇乱又伦精品不卡| 欧美成人午夜精品| 国产91精品成人一区二区三区 | 9热在线视频观看99| 精品国产一区二区三区四区第35| 国产免费福利视频在线观看| 亚洲精品国产一区二区精华液| 九色亚洲精品在线播放| 国产欧美日韩一区二区三| 亚洲人成电影观看| 夜夜骑夜夜射夜夜干| www.999成人在线观看| 免费久久久久久久精品成人欧美视频| 久久久国产欧美日韩av| 国产黄频视频在线观看| 精品久久久精品久久久| a级片在线免费高清观看视频| 亚洲成国产人片在线观看| 少妇粗大呻吟视频| 91国产中文字幕| 曰老女人黄片| 性色av乱码一区二区三区2| 日日爽夜夜爽网站| 亚洲久久久国产精品| 国产激情久久老熟女| 日本五十路高清| 香蕉国产在线看| 午夜91福利影院| 国产不卡av网站在线观看| 老司机午夜福利在线观看视频 | 天堂中文最新版在线下载| 人人妻人人添人人爽欧美一区卜| 两个人免费观看高清视频| 欧美黄色淫秽网站| 91字幕亚洲| 99riav亚洲国产免费| 色婷婷久久久亚洲欧美| 国产男女超爽视频在线观看| 精品一区二区三区av网在线观看 | 桃花免费在线播放| 国产在线一区二区三区精| 超色免费av| 成年人午夜在线观看视频| 久久国产精品大桥未久av| 久久青草综合色| kizo精华| 成人18禁在线播放| 国产精品香港三级国产av潘金莲| 久久久久国产一级毛片高清牌| 午夜精品国产一区二区电影| 中文字幕人妻熟女乱码| 国产精品久久久久久精品古装| 男女边摸边吃奶| 99国产精品一区二区三区| videosex国产| 啦啦啦视频在线资源免费观看| 18禁国产床啪视频网站| 制服人妻中文乱码| 宅男免费午夜| 国产精品一区二区在线不卡| 在线永久观看黄色视频| 在线播放国产精品三级| 国产精品国产av在线观看| 日韩中文字幕视频在线看片| 国产野战对白在线观看| 国产97色在线日韩免费| 色在线成人网| 亚洲成国产人片在线观看| 乱人伦中国视频| 久9热在线精品视频| 日韩中文字幕欧美一区二区| 嫁个100分男人电影在线观看| 女性被躁到高潮视频| 免费高清在线观看日韩| 午夜91福利影院| 久久精品国产亚洲av香蕉五月 | 精品国内亚洲2022精品成人 | 国产亚洲精品第一综合不卡| 老司机午夜十八禁免费视频| 精品人妻1区二区| 免费在线观看黄色视频的| 欧美日韩黄片免| 高清毛片免费观看视频网站 | 91大片在线观看| 精品视频人人做人人爽| 熟女少妇亚洲综合色aaa.| 国产在线观看jvid| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲少妇的诱惑av| 十分钟在线观看高清视频www| 嫩草影视91久久| 99精品久久久久人妻精品| 超碰97精品在线观看| 亚洲一区二区三区欧美精品| 亚洲五月色婷婷综合| 日韩人妻精品一区2区三区| 色老头精品视频在线观看| 性色av乱码一区二区三区2| 香蕉久久夜色| 日韩欧美免费精品| 在线十欧美十亚洲十日本专区| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品一区二区三区| 亚洲av成人一区二区三| 亚洲av日韩精品久久久久久密| 啦啦啦在线免费观看视频4| 老汉色av国产亚洲站长工具| 丁香六月欧美| 人人妻人人澡人人爽人人夜夜| 一级毛片精品| 老汉色av国产亚洲站长工具| 黄色毛片三级朝国网站| 国产亚洲欧美在线一区二区| 丝袜美足系列| 久久人妻福利社区极品人妻图片| 一本大道久久a久久精品| 香蕉丝袜av| 91九色精品人成在线观看| 日韩视频一区二区在线观看| 久久精品亚洲精品国产色婷小说| 狠狠婷婷综合久久久久久88av| 黑人操中国人逼视频| 国产一区二区三区在线臀色熟女 | 精品一区二区三区av网在线观看 | svipshipincom国产片| 午夜福利欧美成人| 午夜久久久在线观看| 国产主播在线观看一区二区| 成人影院久久| 免费在线观看完整版高清| 国产一区二区 视频在线| 国产伦理片在线播放av一区| 国产aⅴ精品一区二区三区波| 欧美老熟妇乱子伦牲交| 精品少妇久久久久久888优播| 欧美黑人欧美精品刺激| 免费在线观看日本一区| 岛国毛片在线播放| 欧美日韩国产mv在线观看视频| 极品教师在线免费播放| 老司机亚洲免费影院| 99国产精品一区二区蜜桃av | 久久久欧美国产精品| 国产片内射在线| 18禁国产床啪视频网站| 看免费av毛片| 久久久水蜜桃国产精品网| 中文字幕av电影在线播放| 亚洲av日韩精品久久久久久密| 男男h啪啪无遮挡| 久久久国产精品麻豆| 精品卡一卡二卡四卡免费| 国产精品一区二区在线观看99| 丝袜美足系列| 国产有黄有色有爽视频| 曰老女人黄片| 黄片大片在线免费观看| 国产成人欧美| 老司机午夜十八禁免费视频| 搡老岳熟女国产| 国产亚洲欧美精品永久| 久久久欧美国产精品| 久久热在线av| 大香蕉久久网| 亚洲专区字幕在线| 母亲3免费完整高清在线观看| 亚洲欧美日韩另类电影网站| 波多野结衣一区麻豆| 午夜精品国产一区二区电影| 亚洲免费av在线视频| 日本欧美视频一区| 2018国产大陆天天弄谢| 黄片小视频在线播放| 精品国产乱子伦一区二区三区| 岛国毛片在线播放| 午夜福利影视在线免费观看| 人成视频在线观看免费观看| 国产亚洲一区二区精品| 99久久国产精品久久久| 国产成人系列免费观看| 在线播放国产精品三级| 老汉色av国产亚洲站长工具| 成人精品一区二区免费| 51午夜福利影视在线观看| 精品福利观看| 大香蕉久久成人网| 亚洲专区中文字幕在线| 可以免费在线观看a视频的电影网站| 国产精品香港三级国产av潘金莲| 成人18禁高潮啪啪吃奶动态图| 精品人妻熟女毛片av久久网站| 亚洲精品av麻豆狂野| 色94色欧美一区二区| 久久精品亚洲精品国产色婷小说| av超薄肉色丝袜交足视频| 国产亚洲一区二区精品| 午夜老司机福利片| 日韩视频在线欧美| 久久热在线av| 国产精品一区二区在线观看99| 无限看片的www在线观看| 满18在线观看网站| 中文亚洲av片在线观看爽 | www.999成人在线观看| 淫妇啪啪啪对白视频| 日韩熟女老妇一区二区性免费视频| 一本久久精品| 交换朋友夫妻互换小说| 国产日韩欧美视频二区| 日韩有码中文字幕| 丝袜喷水一区| www.熟女人妻精品国产| 久久久久精品人妻al黑| 咕卡用的链子| 母亲3免费完整高清在线观看| 天堂俺去俺来也www色官网| 精品国产乱码久久久久久小说| 91成人精品电影| 国产亚洲午夜精品一区二区久久| 国产精品二区激情视频| 王馨瑶露胸无遮挡在线观看| 国产福利在线免费观看视频| 黄网站色视频无遮挡免费观看| av福利片在线| 欧美av亚洲av综合av国产av| av一本久久久久| 国产在线视频一区二区| 香蕉久久夜色| 麻豆av在线久日| 少妇粗大呻吟视频| 首页视频小说图片口味搜索| 欧美变态另类bdsm刘玥| 啦啦啦在线免费观看视频4| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品一卡2卡三卡4卡5卡| 成人黄色视频免费在线看| 精品人妻1区二区| 国产精品久久电影中文字幕 | 精品国产乱码久久久久久男人| 777久久人妻少妇嫩草av网站| 国产av又大| 男女午夜视频在线观看| 一级毛片电影观看| 亚洲精品美女久久久久99蜜臀| 天天躁日日躁夜夜躁夜夜| 人成视频在线观看免费观看| 亚洲精品中文字幕在线视频| 免费不卡黄色视频| 考比视频在线观看| 免费一级毛片在线播放高清视频 | 精品国产亚洲在线| 日本vs欧美在线观看视频| 国产精品香港三级国产av潘金莲| 午夜免费鲁丝| 国产亚洲精品久久久久5区| 精品少妇一区二区三区视频日本电影| 午夜精品久久久久久毛片777| 欧美av亚洲av综合av国产av| 国产一区二区 视频在线| 精品久久久久久久毛片微露脸| 老司机在亚洲福利影院| 中文字幕人妻熟女乱码| 美女午夜性视频免费| 我的亚洲天堂| 少妇裸体淫交视频免费看高清 | 啦啦啦中文免费视频观看日本| 叶爱在线成人免费视频播放| 亚洲五月色婷婷综合| 精品国产一区二区久久| 久久久精品免费免费高清| 丁香六月欧美| 天天操日日干夜夜撸| 12—13女人毛片做爰片一| 日韩大片免费观看网站| 蜜桃在线观看..| 亚洲av片天天在线观看| 午夜福利在线观看吧| 伊人久久大香线蕉亚洲五| 国产成人啪精品午夜网站| 国产欧美日韩精品亚洲av| 无限看片的www在线观看| 在线观看免费日韩欧美大片| videos熟女内射| 伦理电影免费视频| 国产男女超爽视频在线观看| 亚洲精品成人av观看孕妇| 国产成人av教育| 久热这里只有精品99| 美女福利国产在线| av网站免费在线观看视频| 三上悠亚av全集在线观看| 熟女少妇亚洲综合色aaa.| 国产欧美亚洲国产| 老熟妇乱子伦视频在线观看| 午夜精品久久久久久毛片777| 老熟妇乱子伦视频在线观看| 国产欧美日韩综合在线一区二区| 国产淫语在线视频| 免费一级毛片在线播放高清视频 | 久久久久久人人人人人| 国产精品久久电影中文字幕 | 精品福利永久在线观看| 大型黄色视频在线免费观看| 一级毛片电影观看| 国产精品自产拍在线观看55亚洲 | 丁香欧美五月| 免费在线观看日本一区| cao死你这个sao货| 91字幕亚洲| 丰满人妻熟妇乱又伦精品不卡| 亚洲av电影在线进入| 国产亚洲av高清不卡| 操美女的视频在线观看| 一本一本久久a久久精品综合妖精| 国产精品久久久久成人av| 亚洲熟女精品中文字幕| 久久精品熟女亚洲av麻豆精品| 丁香六月欧美| 久久天躁狠狠躁夜夜2o2o| 十八禁高潮呻吟视频| 久久影院123| 叶爱在线成人免费视频播放| 999久久久精品免费观看国产| 亚洲国产欧美在线一区| 亚洲精品美女久久av网站| 又黄又粗又硬又大视频| 欧美人与性动交α欧美精品济南到| av在线播放免费不卡| 99国产精品99久久久久| 99riav亚洲国产免费| tube8黄色片| 免费在线观看日本一区| 国产精品 国内视频| 老熟女久久久| 亚洲情色 制服丝袜| 亚洲av日韩在线播放| 日韩免费av在线播放| svipshipincom国产片| 国产av精品麻豆| 国产一区二区三区视频了| 欧美另类亚洲清纯唯美| 国产一区二区三区在线臀色熟女 | 午夜激情av网站| 日韩一卡2卡3卡4卡2021年| 美女高潮喷水抽搐中文字幕| www日本在线高清视频| 精品熟女少妇八av免费久了| 精品人妻在线不人妻| 女人爽到高潮嗷嗷叫在线视频| 国产av又大| 一区二区日韩欧美中文字幕| 99久久99久久久精品蜜桃| 嫁个100分男人电影在线观看| 国产日韩欧美在线精品| 成人精品一区二区免费| 国产精品偷伦视频观看了| 午夜福利欧美成人| 王馨瑶露胸无遮挡在线观看| 欧美成狂野欧美在线观看| 成年动漫av网址| 一级毛片精品| 十八禁人妻一区二区| 日本欧美视频一区| 色综合婷婷激情| 国产色视频综合| 精品一区二区三卡| 99九九在线精品视频| 国产欧美日韩一区二区精品| 日韩欧美一区二区三区在线观看 | 国产精品亚洲一级av第二区| 久久精品国产99精品国产亚洲性色 | 91精品三级在线观看| 激情在线观看视频在线高清 |