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

    中法海洋衛(wèi)星微波散射計(jì)風(fēng)場(chǎng)反演殘差特性研究

    2021-12-13 03:46:06王冰花董曉龍林文明郎姝燕朱迪云日升
    海洋學(xué)報(bào) 2021年11期
    關(guān)鍵詞:概率模型風(fēng)場(chǎng)殘差

    王冰花,董曉龍*,林文明,郎姝燕,朱迪,云日升

    ( 1. 中國(guó)科學(xué)院大學(xué),北京 100049;2. 中國(guó)科學(xué)院國(guó)家空間科學(xué)中心 微波遙感技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京 100190;3. 南京信息工程大學(xué) 海洋科學(xué)學(xué)院,江蘇 南京 210044;4. 國(guó)家衛(wèi)星海洋應(yīng)用中心,北京 100081)

    1 引言

    海面風(fēng)場(chǎng)是海洋和大氣科學(xué)研究與應(yīng)用中的重要物理參數(shù)。星載微波散射計(jì)作為當(dāng)前獲取全球海面風(fēng)場(chǎng)信息最主要的遙感儀器,其觀測(cè)結(jié)果在數(shù)值天氣預(yù)報(bào)、海洋災(zāi)害監(jiān)測(cè)、海洋環(huán)境數(shù)值預(yù)報(bào)、氣象預(yù)報(bào)和氣候研究等方面發(fā)揮了重要的作用[1-2]。

    散射計(jì)觀測(cè)風(fēng)場(chǎng)質(zhì)量對(duì)數(shù)值天氣預(yù)報(bào)的精度有著重要影響。因此,在海面風(fēng)場(chǎng)的諸多應(yīng)用中,對(duì)風(fēng)場(chǎng)質(zhì)量的定量化描述尤為重要,通常用風(fēng)場(chǎng)質(zhì)量敏感因子來(lái)實(shí)現(xiàn),其中最常用的一種質(zhì)量敏感因子是風(fēng)場(chǎng)反演代價(jià)函數(shù)的殘差[3]。衛(wèi)星散射計(jì)海面風(fēng)場(chǎng)反演最常用的算法是最大似然估計(jì)法(Maximum Likelihood Estimation, MLE),風(fēng)場(chǎng)反演的殘差也稱作MLE值,表示一組測(cè)量的后向散射系數(shù)( σ0)與地球物理模型函數(shù)(Geophysical Model Function, GMF)構(gòu)成多維參考面之間的距離[4]。通常情況下,散射計(jì)測(cè)量的海面 σ0與GMF模擬的 σ0最為接近,反演的風(fēng)場(chǎng)質(zhì)量高、MLE值小。然而,當(dāng)風(fēng)單元(Wind Vector Cell, WVC)中的雷達(dá)后向散射信號(hào)主要受到除風(fēng)以外的其他地球物理?xiàng)l件的影響時(shí),實(shí)際測(cè)量的 σ0和GMF仿真得到的σ0值會(huì)有較大的差異,導(dǎo)致MLE值偏大[5],由此可見(jiàn),MLE是一個(gè)很好的風(fēng)場(chǎng)質(zhì)量指示標(biāo)志。另一方面,在風(fēng)場(chǎng)去模糊處理過(guò)程中,MLE值越大,相應(yīng)的模糊解被選中為“真實(shí)解”的可能性的就越小,反之亦然。因此,MLE值還包含模糊解為真實(shí)解的概率,廣泛用于數(shù)據(jù)同化和二維變分分析模糊去除(2D_Var)[6]。

    圖1為中法海洋衛(wèi)星散射計(jì)(CSCAT)風(fēng)場(chǎng)反演流程圖,圖1中輸入的L1B級(jí)數(shù)據(jù)包括條帶后向散射系數(shù)及其輔助數(shù)據(jù),其中條帶后向散射系數(shù)要經(jīng)過(guò)面元匹配、外定標(biāo)等過(guò)程變換成網(wǎng)格后向散射系數(shù)才能進(jìn)行風(fēng)場(chǎng)反演。面元匹配是將按時(shí)序存儲(chǔ)的后向散射數(shù)據(jù)以及后向散射系數(shù)對(duì)應(yīng)的觀測(cè)幾何信息等參數(shù),重采樣到風(fēng)矢量單元。外定標(biāo)的目的是校正內(nèi)定標(biāo)環(huán)路之外的儀器不確定因素引起的后向散射系數(shù)偏差,對(duì)每個(gè)獨(dú)立的測(cè)量都要進(jìn)行外定標(biāo),外定標(biāo)的精度會(huì)影響殘差的數(shù)值結(jié)果,外定標(biāo)精度越低,風(fēng)單元多次觀測(cè)之間的不一致性越大,殘差值就越大,風(fēng)場(chǎng)反演的質(zhì)量越低。通過(guò)面元匹配和外定標(biāo)處理生成的L2A級(jí)網(wǎng)格后向散射系數(shù),為后續(xù)風(fēng)矢量反演提供輸入?yún)?shù)。

    圖1 CSCAT風(fēng)場(chǎng)反演流程圖Fig. 1 Wind retrieval flow chart of CSCAT

    L2A級(jí)后向散射系數(shù)通過(guò)最大似然估計(jì)得到風(fēng)場(chǎng)模糊解。風(fēng)場(chǎng)模糊解在模糊去除之前需通過(guò)模糊解似然概率模型函數(shù)將反演的MLE值轉(zhuǎn)換為真實(shí)解的概率。概率模型函數(shù)的參數(shù)與MLE值特性有關(guān),目前CSCAT標(biāo)準(zhǔn)的風(fēng)場(chǎng)反演流程采用的是筆形波束旋轉(zhuǎn)掃描散射計(jì)的經(jīng)驗(yàn)概率模型函數(shù)。由于CSCAT散射計(jì)采用了新的觀測(cè)幾何(圖2),這種觀測(cè)幾何大大增加了對(duì)風(fēng)單元的觀測(cè)次數(shù)。以往的散射計(jì)對(duì)同一海面風(fēng)矢量單元的觀測(cè)次數(shù)一般為2~4次,CSCAT對(duì)刈幅內(nèi)每一個(gè)風(fēng)單元的觀測(cè)次數(shù)高達(dá)4~16次,每個(gè)觀測(cè)視數(shù)包含25~50個(gè)觀測(cè)條帶后向散射信息,提供了遠(yuǎn)比其他微波散射計(jì)更多的后向散射信息。由于CSCAT的觀測(cè)幾何與以往散射計(jì)的都不相同,其反演殘差特性也會(huì)有所不同,因此本文對(duì)CSCAT反演殘差特性進(jìn)行了詳細(xì)的研究并建立了與風(fēng)單元位置相關(guān)的似然概率模型函數(shù)。

    圖2 CSCAT觀測(cè)幾何示意圖Fig. 2 Schematic diagram of observation geometry of CSCAT

    CSCAT在數(shù)據(jù)處理中采用二維變分分析法進(jìn)行模糊解去除[7],二維變分分析是將數(shù)據(jù)同化到數(shù)值天氣預(yù)報(bào)(Numerical Weather Prediction, NWP)模型中常用的方法[6]。由于降雨、海冰、未知海況等會(huì)影響后向散射系數(shù)從而降低風(fēng)場(chǎng)反演質(zhì)量[8-9],所以要對(duì)模糊去除之后的風(fēng)場(chǎng)進(jìn)行質(zhì)量控制以檢測(cè)和剔除不合格的風(fēng)矢量單元,最終獲得L2B級(jí)風(fēng)場(chǎng)產(chǎn)品。

    本文基于最大似然估計(jì)風(fēng)場(chǎng)反演算法,研究CSCAT觀測(cè)幾何下風(fēng)場(chǎng)反演代價(jià)函數(shù)的殘差特性。

    2 數(shù)據(jù)與模型

    2.1 數(shù)據(jù)

    CSCAT是國(guó)際上第一個(gè)在軌運(yùn)行的扇形波束旋轉(zhuǎn)掃描微波散射計(jì),與固定扇形波束和旋轉(zhuǎn)筆形波束散射計(jì)相比,CSCAT可以進(jìn)行大入射角范圍、多方位向、覆蓋范圍廣的后向散射系數(shù)觀測(cè)[10]。這種新的觀測(cè)機(jī)制為海面風(fēng)場(chǎng)反演提供了比以往散射計(jì)更多的觀測(cè)信息,在提高風(fēng)場(chǎng)反演質(zhì)量的同時(shí),也給風(fēng)場(chǎng)反演代價(jià)函數(shù)的構(gòu)建和求解帶來(lái)新的挑戰(zhàn)。Liu等[7]研究發(fā)現(xiàn),CSCAT在2019年期間的數(shù)據(jù)都比較穩(wěn)定,所以我們采用CSCAT 2019年3月和5月3.0版本L2B級(jí)風(fēng)場(chǎng)數(shù)據(jù)研究CSCAT風(fēng)場(chǎng)反演代價(jià)函數(shù)的殘差特性。

    CSCAT的刈幅在沿軌方向劃分成一系列行,每一行包括42個(gè)分辨率為25 km的風(fēng)單元,風(fēng)單元編號(hào)從最左端開始向最右端(以衛(wèi)星前進(jìn)方向?yàn)閰⒖迹┮来芜f增,如圖3所示。其中,1~5列和38~42列為刈幅遠(yuǎn)端,6~17列和26~37列為刈幅中間部分,18~25列為星下點(diǎn)區(qū)域。

    圖3 CSCAT地面風(fēng)單元?jiǎng)澐諪ig. 3 CSCAT ground wind vector cell meshing

    風(fēng)場(chǎng)數(shù)據(jù)的輔助數(shù)據(jù)主要有海冰、海陸及歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium Range Weather Forecasts,ECMWF)數(shù)據(jù)。海陸、海冰數(shù)據(jù)用于地面類型識(shí)別,以海陸分布圖和海冰分布圖作為基本判斷標(biāo)準(zhǔn),結(jié)合每個(gè)風(fēng)矢量單位(Wind Vector Cell,WVC)中心的地理位置,對(duì)WVC所對(duì)應(yīng)的地面類型進(jìn)行判斷。其中,海陸分布圖使用的是ECMWF提供的高分辨率海陸掩膜。海冰掩膜是由歐洲氣象衛(wèi)星組織海洋海冰應(yīng)用中心提供的每日南北半球的海冰邊緣線融合產(chǎn)品。參考風(fēng)場(chǎng)采用了相應(yīng)時(shí)間范圍的ECMWF預(yù)報(bào)風(fēng)場(chǎng),ECMWF每天在0點(diǎn)和12點(diǎn)發(fā)布預(yù)報(bào)風(fēng)場(chǎng),預(yù)報(bào)風(fēng)場(chǎng)時(shí)間間隔為3 h,空間網(wǎng)格分辨率為0.125°,通過(guò)空間雙線性插值和時(shí)間樣條插值獲取散射計(jì)觀測(cè)位置和時(shí)間的匹配風(fēng)場(chǎng),進(jìn)行時(shí)間插值時(shí)應(yīng)有兩個(gè)背景風(fēng)場(chǎng)數(shù)據(jù)在散射計(jì)測(cè)量時(shí)間點(diǎn)之前,一個(gè)背景風(fēng)場(chǎng)數(shù)據(jù)在散射計(jì)測(cè)量點(diǎn)之后。考慮到ECMWF預(yù)報(bào)數(shù)據(jù)獲取的時(shí)效性,本研究使用的預(yù)報(bào)時(shí)間在3~18 h之間。

    2.2 CSCAT風(fēng)場(chǎng)反演模型

    隨著散射計(jì)技術(shù)的不斷發(fā)展和完善,形成了多種風(fēng)場(chǎng)反演方法,其中最常用的是最大似然估計(jì)[11],CSCAT最大似然估計(jì)的代價(jià)函數(shù)定義為[12]

    式中,N是獨(dú)立觀測(cè)的次數(shù);i是獨(dú)立觀測(cè)的序號(hào);和分 別表示第i個(gè)視圖觀測(cè)的后向散射系數(shù)和仿真的后向散射系數(shù);Kpi是后向散射測(cè)量的歸一化標(biāo)準(zhǔn)差,可以根據(jù)下式進(jìn)行估計(jì)[13]:

    式中,M(通常為20~30)是分到第i個(gè)視圖的條帶總數(shù);是第j個(gè)條帶的后向散射值。

    每個(gè)風(fēng)矢量單元在反演時(shí)都可以得到一個(gè)代價(jià)函數(shù),MLE代價(jià)函數(shù)曲線的形狀很大程度上決定了反演風(fēng)矢量的準(zhǔn)確性,圖4為CSCAT第39列某一風(fēng)單元的MLE代價(jià)函數(shù)與風(fēng)向風(fēng)速曲線關(guān)系的示例。MLE與風(fēng)向曲線關(guān)系中存在4個(gè)極小值點(diǎn),由此可以確定,在MLE反演過(guò)程檢索到對(duì)應(yīng)的4個(gè)模糊解,代價(jià)函數(shù)曲線取得極小值時(shí)附近的曲線形狀和模糊解的真實(shí)性有關(guān)。第1個(gè)與第3個(gè)極小值附近的曲線比較寬,這時(shí)反演得到的風(fēng)矢量偏差較大。因?yàn)榇藭r(shí)極小值附近的風(fēng)矢量是真實(shí)風(fēng)矢量的概率和極小值對(duì)應(yīng)的概率相當(dāng),但是這些概率相當(dāng)?shù)腗LE中只有一個(gè)被選中。

    圖4 MLE風(fēng)場(chǎng)反演曲線Fig. 4 Curve of MLE wind inversion

    為了避免代價(jià)函數(shù)曲線極小值點(diǎn)較寬導(dǎo)致的風(fēng)場(chǎng)反演偏差較大的情況,CSCAT結(jié)合多解方案(Multiple Solution Scheme, MSS)[14]進(jìn)行風(fēng)場(chǎng)反演,MSS反演中通常將風(fēng)向從0°到360°每隔2.5°進(jìn)行劃分,對(duì)于每個(gè)風(fēng)向,通過(guò)對(duì)代價(jià)函數(shù)求極小值計(jì)算得到與之相對(duì)應(yīng)的風(fēng)矢量。每一個(gè)風(fēng)單元可以得到144個(gè)風(fēng)矢量解,其中只有一個(gè)是真實(shí)解,其余的稱為模糊解。為了得到真實(shí)解還要進(jìn)行風(fēng)向的模糊去除。

    3 殘差特性分析

    3.1 CSCAT的殘差分布特性

    為了分析CSCAT的MLE殘差分布特性,本文研究了MLE平均值及標(biāo)準(zhǔn)差隨風(fēng)速和WVC列數(shù)的變化,并根據(jù)衛(wèi)星數(shù)據(jù)產(chǎn)品的質(zhì)量標(biāo)識(shí)符,剔除了陸地、海冰區(qū)域以及其他質(zhì)量不合格的風(fēng)矢量單元數(shù)據(jù)。然后計(jì)算相同節(jié)點(diǎn)相同風(fēng)速下所有被選中MLE值的平均值和標(biāo)準(zhǔn)差(圖5),通過(guò)風(fēng)場(chǎng)反演會(huì)得到4個(gè)模糊風(fēng)矢量及其對(duì)應(yīng)的MLE值,通過(guò)模糊去除選擇真實(shí)風(fēng)矢量的模糊解對(duì)應(yīng)的MLE值,即被選中的MLE值。從圖中可以明顯看出MLE的分布與風(fēng)速有很強(qiáng)的相關(guān)性。隨著風(fēng)速的減小,MLE均值和標(biāo)準(zhǔn)差均呈上升趨勢(shì);由于CSCAT的觀測(cè)幾何隨著節(jié)點(diǎn)數(shù)變化較快,不同列的風(fēng)單元之間也存在明顯差異,在刈幅遠(yuǎn)端的風(fēng)速偏差和標(biāo)準(zhǔn)差都較大,因?yàn)樨追h(yuǎn)端為大入射角觀測(cè),σ0的精度降低[7],導(dǎo)致觀測(cè)的 σ0和GMF得到的仿真值差別較大。星下點(diǎn)區(qū)域在中等風(fēng)速時(shí)偏差和標(biāo)準(zhǔn)差明顯增大是因?yàn)樵搮^(qū)域觀測(cè)入射角信息豐富,觀測(cè)方位角集中在180°左右,由于非風(fēng)場(chǎng)因素的地球物理噪聲的影響及質(zhì)量控制的效果在中等風(fēng)速出現(xiàn)的概率最大,質(zhì)量控制雖然能剔除一部分海冰、降雨等因素引起的質(zhì)量較差的觀測(cè),但無(wú)法完全識(shí)別和剔除,所以在中等風(fēng)速時(shí)偏差和標(biāo)準(zhǔn)差明顯增大。

    圖5 平均MLE隨風(fēng)速和WVC列數(shù)的變化(a)及MLE標(biāo)準(zhǔn)差隨風(fēng)速和WVC列數(shù)的變化(b)Fig. 5 Mean MLE versus wind speed and WVC number (a) and MLE standard deriation versus wind speed and WVC number (b)

    圖6、圖7是同一時(shí)間內(nèi)不同節(jié)點(diǎn)下CSCAT風(fēng)速、風(fēng)向相對(duì)于ECMWF風(fēng)速、風(fēng)向的偏差和標(biāo)準(zhǔn)差隨平均MLE變化的曲線。對(duì)比偏差和標(biāo)準(zhǔn)差可以明顯看出不同節(jié)點(diǎn)的風(fēng)速、風(fēng)向偏差存在明顯差異,MLE值較大時(shí),不同節(jié)點(diǎn)的風(fēng)速差異可達(dá)到0.2 m/s。風(fēng)速、風(fēng)向的標(biāo)準(zhǔn)差隨MLE的增大呈上升趨勢(shì),不同節(jié)點(diǎn)之間略有不同,所以要區(qū)分不同的節(jié)點(diǎn)研究反演的殘差特性。

    圖6 CSCAT風(fēng)速、風(fēng)向相對(duì)于ECMWF風(fēng)速、風(fēng)向的偏差隨平均MLE變化的曲線Fig. 6 The variation curve of CSCAT wind speed and direction bias relative to ECMWF wind speed and direction with mean MLEa. 風(fēng)速偏差;b. 風(fēng)向偏差;不同線型的曲線表示不同的WVC列數(shù)a. Biases of wind speed; b. biases of wind direction; lines in different styles are for different WVC number

    圖7 CSCAT風(fēng)速、風(fēng)向相對(duì)于ECMWF風(fēng)速、風(fēng)向的標(biāo)準(zhǔn)差隨平均MLE變化的曲線Fig. 7 The variation curve of CSCAT wind speed and direction standard deviation relative to ECMWF wind speed and direction with mean MLEa. 風(fēng)速標(biāo)準(zhǔn)差;b. 風(fēng)向標(biāo)準(zhǔn)差;不同線型的曲線表示不同的WVC列數(shù)a. Standard deviation of wind speed differences; b. standard deviation of wind direction differences; lines in different styles are for different WVC number

    3.2 似然概率模型函數(shù)

    反演得到的模糊解要通過(guò)模糊去除選擇一個(gè)模糊解作為真實(shí)解,模糊解被選為真實(shí)解的概率即模糊解似然概率,根據(jù)貝葉斯理論[15]和公式(1),散射計(jì)反演的模糊風(fēng)矢量是真實(shí)風(fēng)矢量的概率可以表示為與MLE有關(guān)的函數(shù):

    式中,v表示“真實(shí)”風(fēng)速; σ0表示一系列散射計(jì)后向散射系數(shù)測(cè)量值;k是標(biāo)準(zhǔn)化比例因子。因此,該概率函數(shù)理論上是一個(gè)和MLE相關(guān)的指數(shù)函數(shù)關(guān)系。從上式可以看出MLE值越小,模糊解成為真實(shí)解的概率越大。模糊風(fēng)矢量越接近真實(shí)風(fēng)矢量,MLE值越小,由公式(3)得到的概率值就越大。概率值較大的模糊風(fēng)矢量也越接近真實(shí)風(fēng)矢量,因此概率函數(shù)模型與MLE值的含義是一致的。在實(shí)際應(yīng)用中,由于存在一些觀測(cè)誤差,指數(shù)函數(shù)的形狀可能與理論值不同。

    Portabella和Stoffelen[16]首先提出在QuickSCAT質(zhì)量控制中采用歸一化的MLE值或歸一化的殘差替代MLE,目的是解決MLE反演方法過(guò)程中存在的問(wèn)題和噪聲估計(jì),Portabella和Stoffelen[16]通過(guò)統(tǒng)計(jì)分析得出QuikSCAT的概率模型函數(shù)表達(dá)式為

    式中,x代表歸一化殘差。

    本研究中,設(shè)CSCAT的概率模型函數(shù)表達(dá)式為

    式中,k為歸一化因子;l為要根據(jù)經(jīng)驗(yàn)推導(dǎo)的參數(shù);Rn是歸一化的殘差。為了獲得經(jīng)驗(yàn)公式(5)的具體形式,可以忽略先驗(yàn)知識(shí)對(duì)指數(shù)關(guān)系的影響,并假設(shè)存在一個(gè)函數(shù)關(guān)系ps(x),對(duì)于某一風(fēng)矢量單元,如果有N個(gè)模糊風(fēng)矢量vi和相應(yīng)的歸一化最大似然值Rnj,則其中第j個(gè)風(fēng)矢量最接近真實(shí)風(fēng)矢量的概率可以表示為

    為了得到ps(x)函數(shù)的具體形式,本文對(duì)CSCAT L2B級(jí)數(shù)據(jù)進(jìn)行了分析,研究只有兩個(gè)模糊解的風(fēng)矢量單元。選擇與ECMWF風(fēng)場(chǎng)最接近的模糊解作為選擇解。統(tǒng)計(jì)每一階的模糊解作為選擇解的個(gè)數(shù),可以得到每階模糊解為真實(shí)風(fēng)解的概率。在公式中令N為2,可以的到第一階模糊解成為真實(shí)解的概率是:

    通過(guò)公式(7)的運(yùn)算可以得到每階統(tǒng)計(jì)概率,進(jìn)而得到ps(Rn2)/ps(Rn1) 的估計(jì)值。圖8是Rn1為0.1時(shí)概率函數(shù)曲線隨節(jié)點(diǎn)位置的變化情況,可以看出星下點(diǎn)下降速率最慢,刈幅遠(yuǎn)端下降速率最快,表明隨著MLE的增大,模糊解成為真實(shí)解的概率急劇減小,說(shuō)明較小的MLE值對(duì)應(yīng)的解成為真實(shí)解的概率更大,有更多的風(fēng)單元選擇一階模糊解作為真實(shí)風(fēng)矢量解。

    圖8 不同節(jié)點(diǎn)下( Rn2?Rn1 ) 值與 p(Rn2)/p(Rn1) 函數(shù)關(guān)系的指數(shù)擬合( R n1=0.1)Fig. 8 The exponential fit to ratio R n2 and R n1 as a founction of (R n2?Rn1) value for different WVC number( Rn1=0.1)

    圖9為通過(guò)實(shí)驗(yàn)數(shù)據(jù)確定的第8列所有風(fēng)單元的Rn1和Rn2組合的比率函數(shù)關(guān)系。盡管曲線較為離散,但是總體明確表明當(dāng)Rn1小 于2時(shí),P(Rn2)/P(Rn1)比值是關(guān)于Rn1和Rn2不變的指數(shù)關(guān)系函數(shù),符合公式(5)假設(shè)的曲線關(guān)系。當(dāng)Rn1是 一個(gè)常數(shù)時(shí),ps(Rn1)也是常數(shù),圖9表現(xiàn)了ps(x)的曲線形狀。

    從圖9中可以看出,指數(shù)函數(shù)的下降速率與Rn1的取值有關(guān),當(dāng)Rn1較大時(shí)不再符合指數(shù)關(guān)系,為了找出原因?qū)n1的 分布特性進(jìn)行了統(tǒng)計(jì)(圖10),從Rn1分布直方圖可以看出Rn1的絕大多數(shù)值分布集中在0~1之間,Rn1較 大時(shí),由于在此范圍的風(fēng)單元數(shù)量很少,ps(x)不 符合指數(shù)函數(shù)關(guān)系。

    圖9 不同 R n1取 值下 p(Rn2)/p(Rn1)與 ( R n2?Rn1)值函數(shù)關(guān)系的指數(shù)擬合(第8列WVC)Fig. 9 The exponential fit to ratio of R n2 and R n1 as a founction of (R n2?Rn1 ) value for different R n1 (WVC number 8)

    圖10 R n1分布直方圖(節(jié)點(diǎn)數(shù)為8)Fig. 10 The distribution histogram of R n1 (the number of node is 8)

    3.3 結(jié)果與檢驗(yàn)

    通過(guò)以上分析可知,概率函數(shù)的指數(shù)是與風(fēng)單元列數(shù)有關(guān)的,由于不同位置的風(fēng)單元觀測(cè)數(shù)量不一樣,MLE代價(jià)函數(shù)的維度也不一樣,所以不同位置的風(fēng)單元概率函數(shù)的指數(shù)也不一樣。根據(jù)刈幅不同區(qū)域位置和Rn值的不同分布情況,可以得到每一列指數(shù)擬合的系數(shù)和指數(shù),結(jié)果如圖11所示。刈幅不同區(qū)域的概率函數(shù)差別較大,整體呈現(xiàn)一定的對(duì)稱性,指數(shù)l分布在?0.4~?1.8,與筆形波束散射計(jì)的指數(shù)不同,概率函數(shù)的系數(shù)分布在1左右。

    圖11 概率模型函數(shù)的指數(shù)及系數(shù)隨節(jié)點(diǎn)的變化Fig. 11 The exponents and coefficients of the probabilistic model versus node number

    表2 預(yù)測(cè)概率/觀測(cè)概率的分布對(duì)比(刈幅中間)Table 2 Distribution comparision of predicted probability/observed probability (sweet swath)

    為了驗(yàn)證之前的假設(shè)是否正確,能否推廣到有多個(gè)解的情況,本文對(duì)2019年5月份的150軌CSCAT散射計(jì)L2B級(jí)數(shù)據(jù)進(jìn)行了分析,并推廣到有多個(gè)解的情況,統(tǒng)計(jì)每一級(jí)模糊解成為真解的預(yù)測(cè)概率。表1至表3分別對(duì)不同區(qū)域每一級(jí)模糊解的預(yù)測(cè)概率和觀測(cè)概率進(jìn)行的對(duì)比。表中第一行表示在所有研究數(shù)據(jù)中風(fēng)矢量單元模糊解個(gè)數(shù)為2、3、4的風(fēng)矢量單元的總數(shù)。第2行至第5行中左側(cè)數(shù)據(jù)表示每一階的所有模糊解概率值的平均值,即預(yù)測(cè)概率,右側(cè)數(shù)據(jù)表示統(tǒng)計(jì)每一階模糊解是真解(最接近ECMWF風(fēng)場(chǎng))的個(gè)數(shù)占所有模糊解個(gè)數(shù)的百分比,即實(shí)際觀測(cè)概率。通過(guò)對(duì)比每一列的預(yù)測(cè)概率和觀測(cè)概率可以看出明顯的相關(guān)性,刈幅遠(yuǎn)端和中間區(qū)域存在微小差別,這是因?yàn)榍拔臄M合曲線時(shí),曲線本身是有波動(dòng)的,但是總體對(duì)比表明一致性很顯著。星下點(diǎn)差別較大,這是因?yàn)樾窍曼c(diǎn)進(jìn)行擬合時(shí)不是標(biāo)準(zhǔn)的指數(shù)關(guān)系。因此,可以證明之前的假設(shè)是正確的,擬合的與風(fēng)單元位置相關(guān)的概率模型函數(shù)可以用來(lái)計(jì)算某一模糊解為真實(shí)風(fēng)矢量的概率。

    表1 預(yù)測(cè)概率/觀測(cè)概率的分布對(duì)比(刈幅遠(yuǎn)端)Table 1 Distribution comparision of predicted probability/observed probability (far swath)

    表3 預(yù)測(cè)概率/觀測(cè)概率的分布對(duì)比 (星下點(diǎn)區(qū)域)Table 3 Distribution comparision of predicted probability/observed probability (nadir swath)

    為了檢驗(yàn)上文擬合得到的似然概率模型函數(shù)反演風(fēng)場(chǎng)的效果,分別采用與風(fēng)單元位置相關(guān)的概率模型函數(shù)與標(biāo)準(zhǔn)反演過(guò)程中的概率模型函數(shù)進(jìn)行風(fēng)場(chǎng)反演,得到相應(yīng)的風(fēng)產(chǎn)品,并收集美國(guó)浮標(biāo)數(shù)據(jù)中心的浮標(biāo)數(shù)據(jù)與反演風(fēng)場(chǎng)進(jìn)行對(duì)比。圖12和圖13是兩種方法反演的風(fēng)場(chǎng)與浮標(biāo)風(fēng)場(chǎng)的散點(diǎn)圖對(duì)比結(jié)果??梢钥闯霾捎门c列數(shù)相關(guān)的概率模型函數(shù)反演的風(fēng)場(chǎng)偏差減小了0.03 m/s,而標(biāo)準(zhǔn)差未發(fā)生變化,表明風(fēng)速的離散度未發(fā)生變化,而風(fēng)速的準(zhǔn)確性有提高。對(duì)比風(fēng)向變化可以發(fā)現(xiàn)風(fēng)向偏差未發(fā)生變化但風(fēng)向偏差減小了0.2°,風(fēng)向標(biāo)準(zhǔn)差減小了0.5°,表明改進(jìn)的概率模型函數(shù)反演的風(fēng)場(chǎng)精度有明顯的提高。

    圖12 改進(jìn)前概率模型函數(shù)反演的CSCAT風(fēng)場(chǎng)與浮標(biāo)測(cè)量風(fēng)場(chǎng)Fig. 12 CSCAT wind filed and buoy wind filed retrieved by the likelihood probability model before improvement

    圖13 改進(jìn)的概率模型函數(shù)反演的CSCAT風(fēng)場(chǎng)與浮標(biāo)測(cè)量風(fēng)場(chǎng)Fig. 13 CSCAT wind filed and buoy wind filed retrieved by the improved likelihood probability model

    4 總結(jié)與展望

    本文基于最大似然估計(jì)風(fēng)場(chǎng)反演算法,探討了CSCAT新的觀測(cè)體制下反演殘差分布特性。在理論推導(dǎo)的基礎(chǔ)上,利用CSCAT觀測(cè)數(shù)據(jù),定量擬合了似然概率模型函數(shù)與WVC列數(shù)的關(guān)系。

    CSCAT反演殘差與風(fēng)速和WVC列數(shù)有關(guān),MLE平均值、標(biāo)準(zhǔn)差隨風(fēng)速和列數(shù)有較為明顯的變化,反映了CSCAT數(shù)據(jù)用于風(fēng)場(chǎng)反演時(shí)自身的一些特點(diǎn)。不同WVC列數(shù)的CSCAT風(fēng)速風(fēng)向與ECMWF風(fēng)速風(fēng)向差別較大,并且風(fēng)速風(fēng)向標(biāo)準(zhǔn)差變化趨勢(shì)比較一致,偏差變化明顯不同。

    在風(fēng)場(chǎng)反演時(shí),MLE代價(jià)函數(shù)被轉(zhuǎn)化成似然概率代價(jià)函數(shù)用于風(fēng)場(chǎng)的模糊去除,在第3節(jié)本文研究了模糊解是真實(shí)解的概率問(wèn)題,并根據(jù)不同區(qū)域的殘差對(duì)所有節(jié)點(diǎn)進(jìn)行指數(shù)擬合,發(fā)現(xiàn)不同區(qū)域的概率模型函數(shù)參數(shù)有明顯的差別。為了對(duì)本文擬合結(jié)果進(jìn)行驗(yàn)證,本文利用2019年5月份的數(shù)據(jù)對(duì)實(shí)際觀測(cè)概率與通過(guò)本文的概率模型函數(shù)得到的預(yù)測(cè)概率進(jìn)行對(duì)比。結(jié)果表明,本文提出的概率模型函數(shù)能夠有效預(yù)測(cè)模糊解成為真實(shí)解的概率。最后對(duì)改進(jìn)的概率模型函數(shù)反演的風(fēng)場(chǎng)進(jìn)行檢驗(yàn),可以發(fā)現(xiàn)風(fēng)場(chǎng)反演的質(zhì)量有明顯的改進(jìn)。

    分析結(jié)果為后續(xù)CSCAT風(fēng)場(chǎng)質(zhì)量控制和2D_Var算法的精細(xì)化調(diào)整提供了重要的參考。CFOSAT從2018年發(fā)射以來(lái),已經(jīng)獲取了大量寶貴的海面風(fēng)場(chǎng)觀測(cè)數(shù)據(jù),隨著數(shù)據(jù)的積累和算法的改進(jìn),CFOSAT微波散射計(jì)必將為海洋觀測(cè)和應(yīng)用做出更大的貢獻(xiàn)。

    猜你喜歡
    概率模型風(fēng)場(chǎng)殘差
    基于雙向GRU與殘差擬合的車輛跟馳建模
    基于FLUENT的下?lián)舯┝魅S風(fēng)場(chǎng)建模
    在精彩交匯中,理解兩個(gè)概率模型
    基于殘差學(xué)習(xí)的自適應(yīng)無(wú)人機(jī)目標(biāo)跟蹤算法
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    基于停車服務(wù)效率的選擇概率模型及停車量仿真研究
    “最美風(fēng)場(chǎng)”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    側(cè)向風(fēng)場(chǎng)中無(wú)人機(jī)的飛行研究
    一類概率模型的探究與應(yīng)用
    平穩(wěn)自相關(guān)過(guò)程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    欧美日韩av久久| 国产免费一级a男人的天堂| 黄色怎么调成土黄色| 午夜免费观看性视频| 久久 成人 亚洲| 久久久久久久精品精品| 亚洲av.av天堂| 内地一区二区视频在线| 一级毛片久久久久久久久女| 亚洲中文av在线| 少妇被粗大的猛进出69影院 | 毛片一级片免费看久久久久| 国产一区亚洲一区在线观看| 亚洲综合色惰| 国产精品秋霞免费鲁丝片| 亚洲人成网站在线播| 蜜桃久久精品国产亚洲av| 久久人人爽av亚洲精品天堂| 一级二级三级毛片免费看| 久久99一区二区三区| 在线精品无人区一区二区三| 十八禁高潮呻吟视频 | 久久久久国产精品人妻一区二区| 老司机影院成人| 亚洲,一卡二卡三卡| 熟女电影av网| 一级毛片aaaaaa免费看小| 国产91av在线免费观看| av在线播放精品| 啦啦啦在线观看免费高清www| 亚洲真实伦在线观看| 国产熟女欧美一区二区| 日韩在线高清观看一区二区三区| 国产日韩欧美视频二区| 成人无遮挡网站| 国产精品三级大全| 欧美bdsm另类| 男人添女人高潮全过程视频| 久久久久精品性色| 久久青草综合色| 亚州av有码| 男人舔奶头视频| 亚洲成色77777| 蜜桃在线观看..| 精品人妻一区二区三区麻豆| 国产成人aa在线观看| 久久青草综合色| 久久久国产一区二区| 婷婷色综合大香蕉| 人妻夜夜爽99麻豆av| 日韩,欧美,国产一区二区三区| 乱码一卡2卡4卡精品| 97精品久久久久久久久久精品| 国产成人a∨麻豆精品| 777米奇影视久久| 一区在线观看完整版| 亚洲精品aⅴ在线观看| 色94色欧美一区二区| 18+在线观看网站| 欧美 日韩 精品 国产| 亚洲国产色片| 啦啦啦啦在线视频资源| 热re99久久精品国产66热6| 亚洲婷婷狠狠爱综合网| 久久久精品94久久精品| 观看av在线不卡| 免费播放大片免费观看视频在线观看| 日韩中字成人| 亚洲精品自拍成人| 国内少妇人妻偷人精品xxx网站| av在线观看视频网站免费| 日本色播在线视频| 国产极品粉嫩免费观看在线 | 婷婷色综合www| av天堂久久9| 日本猛色少妇xxxxx猛交久久| 亚洲无线观看免费| 国产毛片在线视频| 欧美97在线视频| 国产成人精品一,二区| 亚洲国产毛片av蜜桃av| 中文字幕免费在线视频6| 麻豆乱淫一区二区| 亚洲成人一二三区av| 99热这里只有是精品在线观看| 一区二区三区乱码不卡18| 亚洲av.av天堂| 一区二区av电影网| 日韩不卡一区二区三区视频在线| 国产精品秋霞免费鲁丝片| 哪个播放器可以免费观看大片| 久久国产精品大桥未久av | 亚洲精品乱码久久久久久按摩| 永久网站在线| 高清午夜精品一区二区三区| 岛国毛片在线播放| 不卡视频在线观看欧美| 97超碰精品成人国产| 久久久国产欧美日韩av| 极品教师在线视频| 中文字幕亚洲精品专区| 99热6这里只有精品| 高清欧美精品videossex| 亚洲精品久久午夜乱码| 国产精品欧美亚洲77777| 成人国产av品久久久| 午夜免费观看性视频| 亚洲成人一二三区av| 曰老女人黄片| 日韩 亚洲 欧美在线| 国内精品宾馆在线| a 毛片基地| 99久久精品国产国产毛片| 一区二区三区乱码不卡18| 亚洲国产最新在线播放| 我要看日韩黄色一级片| 久久99精品国语久久久| 日韩一本色道免费dvd| av播播在线观看一区| 精品国产乱码久久久久久小说| 国产女主播在线喷水免费视频网站| 国产成人精品福利久久| 美女脱内裤让男人舔精品视频| 欧美激情国产日韩精品一区| 18+在线观看网站| 国产精品99久久久久久久久| 日韩熟女老妇一区二区性免费视频| 国产成人一区二区在线| 男女无遮挡免费网站观看| 亚洲精品成人av观看孕妇| 国产成人aa在线观看| 在线观看美女被高潮喷水网站| 亚洲精品亚洲一区二区| 国产精品99久久久久久久久| 欧美老熟妇乱子伦牲交| a级片在线免费高清观看视频| 观看免费一级毛片| 又粗又硬又长又爽又黄的视频| 看非洲黑人一级黄片| 成人黄色视频免费在线看| av在线观看视频网站免费| 国产精品国产三级专区第一集| 午夜日本视频在线| 精品人妻偷拍中文字幕| 日韩人妻高清精品专区| 天天躁夜夜躁狠狠久久av| 日本黄大片高清| 国产免费福利视频在线观看| 精品熟女少妇av免费看| 欧美日韩一区二区视频在线观看视频在线| 国产男女超爽视频在线观看| 久久久久久久久久久免费av| 性高湖久久久久久久久免费观看| 亚洲欧美日韩东京热| 久久久久国产网址| 大香蕉久久网| 精品酒店卫生间| 久久久国产精品麻豆| 亚洲电影在线观看av| 赤兔流量卡办理| 成人午夜精彩视频在线观看| av专区在线播放| 大香蕉久久网| 日韩一本色道免费dvd| 汤姆久久久久久久影院中文字幕| 亚洲精品自拍成人| av在线观看视频网站免费| 如日韩欧美国产精品一区二区三区 | 丰满人妻一区二区三区视频av| 中文精品一卡2卡3卡4更新| 综合色丁香网| 国语对白做爰xxxⅹ性视频网站| 国精品久久久久久国模美| 久久久精品94久久精品| 亚洲精品国产色婷婷电影| 我要看日韩黄色一级片| 噜噜噜噜噜久久久久久91| 在线观看免费高清a一片| 大片免费播放器 马上看| 夜夜骑夜夜射夜夜干| 欧美3d第一页| 亚洲精品乱码久久久v下载方式| 日本欧美视频一区| 岛国毛片在线播放| 精品熟女少妇av免费看| 两个人的视频大全免费| 如日韩欧美国产精品一区二区三区 | 日韩一区二区三区影片| 亚洲av综合色区一区| 毛片一级片免费看久久久久| 久久精品久久精品一区二区三区| 大香蕉97超碰在线| 国产69精品久久久久777片| 国产一级毛片在线| 国产乱人偷精品视频| 精品久久久噜噜| 午夜免费观看性视频| 国产毛片在线视频| 久久国产精品大桥未久av | 亚洲色图综合在线观看| 男人和女人高潮做爰伦理| a 毛片基地| 麻豆精品久久久久久蜜桃| 内地一区二区视频在线| 国产高清三级在线| 亚洲国产毛片av蜜桃av| 人妻 亚洲 视频| 亚洲无线观看免费| 免费看光身美女| av福利片在线观看| 国产成人精品久久久久久| 久久热精品热| 午夜av观看不卡| 菩萨蛮人人尽说江南好唐韦庄| 99久久中文字幕三级久久日本| 又粗又硬又长又爽又黄的视频| a级毛片免费高清观看在线播放| 成人漫画全彩无遮挡| 最近中文字幕高清免费大全6| 国内精品宾馆在线| 青春草国产在线视频| 欧美 日韩 精品 国产| 精品亚洲成国产av| 亚洲无线观看免费| 久久精品国产亚洲av涩爱| 成人亚洲欧美一区二区av| 男女免费视频国产| 亚洲美女黄色视频免费看| 九色成人免费人妻av| 热re99久久国产66热| 熟女人妻精品中文字幕| 精品99又大又爽又粗少妇毛片| 欧美高清成人免费视频www| 国产成人精品久久久久久| 久久精品国产自在天天线| 亚洲伊人久久精品综合| 国产成人精品婷婷| 久久ye,这里只有精品| 国产精品久久久久久av不卡| 在线观看www视频免费| 少妇的逼水好多| 午夜福利在线观看免费完整高清在| 只有这里有精品99| 国产欧美日韩综合在线一区二区 | 夫妻午夜视频| 免费观看性生交大片5| 亚洲人与动物交配视频| 国产深夜福利视频在线观看| 看免费成人av毛片| 欧美日韩综合久久久久久| 免费av不卡在线播放| 中文欧美无线码| 一级片'在线观看视频| 成人无遮挡网站| 国产69精品久久久久777片| 成年av动漫网址| 成人午夜精彩视频在线观看| 日韩熟女老妇一区二区性免费视频| a级毛片在线看网站| 99久久精品热视频| 亚洲情色 制服丝袜| 亚洲精品日本国产第一区| 一本—道久久a久久精品蜜桃钙片| 大片电影免费在线观看免费| 国产日韩一区二区三区精品不卡 | 免费人妻精品一区二区三区视频| 少妇人妻 视频| 精品一区在线观看国产| av国产久精品久网站免费入址| 最近的中文字幕免费完整| 欧美成人午夜免费资源| 80岁老熟妇乱子伦牲交| 久久久久久久久久人人人人人人| 亚洲一区二区三区欧美精品| 精品一区二区免费观看| 免费大片18禁| 黄色视频在线播放观看不卡| 日韩在线高清观看一区二区三区| 国产女主播在线喷水免费视频网站| 成人黄色视频免费在线看| 美女xxoo啪啪120秒动态图| 最近的中文字幕免费完整| 中文字幕亚洲精品专区| 欧美老熟妇乱子伦牲交| 免费观看在线日韩| 最黄视频免费看| 久久精品久久精品一区二区三区| 亚洲精品第二区| 国产在线男女| 婷婷色麻豆天堂久久| 多毛熟女@视频| 五月天丁香电影| 日韩视频在线欧美| 精品熟女少妇av免费看| 国产精品国产三级国产av玫瑰| av国产久精品久网站免费入址| 特大巨黑吊av在线直播| 亚洲人成网站在线观看播放| 日韩,欧美,国产一区二区三区| 亚洲精品亚洲一区二区| 伦理电影免费视频| 一级黄片播放器| 国产又色又爽无遮挡免| 91午夜精品亚洲一区二区三区| 一级片'在线观看视频| 最近手机中文字幕大全| 国产日韩欧美亚洲二区| 不卡视频在线观看欧美| 如日韩欧美国产精品一区二区三区 | 嘟嘟电影网在线观看| 精品人妻熟女毛片av久久网站| 国产黄色免费在线视频| 在线观看三级黄色| 91精品国产国语对白视频| 亚洲国产欧美日韩在线播放 | 国产亚洲欧美精品永久| 精品卡一卡二卡四卡免费| 人妻系列 视频| 国产av一区二区精品久久| 国产女主播在线喷水免费视频网站| 少妇人妻久久综合中文| 热re99久久国产66热| 亚洲精品国产av蜜桃| 国产探花极品一区二区| 中文欧美无线码| 中国国产av一级| 国产成人aa在线观看| 久久久精品免费免费高清| 日本黄大片高清| 一边亲一边摸免费视频| 水蜜桃什么品种好| 日韩电影二区| 亚洲欧美精品自产自拍| 亚洲国产av新网站| 中文字幕av电影在线播放| 建设人人有责人人尽责人人享有的| 秋霞伦理黄片| 亚洲国产精品一区二区三区在线| 亚洲美女视频黄频| 国产高清有码在线观看视频| 99热全是精品| 两个人免费观看高清视频 | 天美传媒精品一区二区| 97精品久久久久久久久久精品| 青春草视频在线免费观看| 亚洲第一区二区三区不卡| 桃花免费在线播放| 精品国产国语对白av| 国产精品国产三级国产专区5o| h视频一区二区三区| 欧美三级亚洲精品| 99久久中文字幕三级久久日本| 午夜福利视频精品| 国产一级毛片在线| 麻豆成人av视频| 一区二区三区精品91| 噜噜噜噜噜久久久久久91| 99热这里只有是精品在线观看| 国产在线免费精品| 一区在线观看完整版| 卡戴珊不雅视频在线播放| 国产精品蜜桃在线观看| 中文资源天堂在线| 久久免费观看电影| 亚洲av男天堂| 精品一区二区三区视频在线| 亚洲人成网站在线观看播放| 国产精品福利在线免费观看| 国产黄片视频在线免费观看| 国产91av在线免费观看| 黄色日韩在线| 18禁裸乳无遮挡动漫免费视频| 日本欧美国产在线视频| 久久久欧美国产精品| 成人黄色视频免费在线看| 久久97久久精品| 久久久久久久久久久丰满| 波野结衣二区三区在线| 亚洲怡红院男人天堂| 亚洲四区av| 熟女av电影| 免费播放大片免费观看视频在线观看| 丰满少妇做爰视频| 日本欧美视频一区| 十分钟在线观看高清视频www | 亚洲精品国产av成人精品| 在线观看人妻少妇| 午夜福利在线观看免费完整高清在| 大片电影免费在线观看免费| 亚洲欧美日韩另类电影网站| 国产精品不卡视频一区二区| 精华霜和精华液先用哪个| 久久6这里有精品| 亚洲成人一二三区av| 岛国毛片在线播放| 一本—道久久a久久精品蜜桃钙片| 丝袜脚勾引网站| 在线精品无人区一区二区三| 哪个播放器可以免费观看大片| 国精品久久久久久国模美| 涩涩av久久男人的天堂| 免费观看在线日韩| 青春草亚洲视频在线观看| 欧美3d第一页| 亚洲国产精品一区二区三区在线| 国产国拍精品亚洲av在线观看| 国产乱来视频区| 五月开心婷婷网| 日韩中文字幕视频在线看片| 熟女人妻精品中文字幕| 亚洲精华国产精华液的使用体验| 欧美日韩视频精品一区| 欧美日韩国产mv在线观看视频| 精品一品国产午夜福利视频| 亚洲精品视频女| av国产久精品久网站免费入址| 日本欧美视频一区| 日韩av不卡免费在线播放| 9色porny在线观看| 午夜精品国产一区二区电影| 99精国产麻豆久久婷婷| 又粗又硬又长又爽又黄的视频| 国产中年淑女户外野战色| 免费看光身美女| 亚洲综合色惰| 国产精品99久久99久久久不卡 | 国产精品人妻久久久影院| 狂野欧美激情性bbbbbb| 黄色配什么色好看| 一本—道久久a久久精品蜜桃钙片| 丝袜脚勾引网站| 高清视频免费观看一区二区| 亚洲欧洲国产日韩| av不卡在线播放| 久久久欧美国产精品| 少妇被粗大猛烈的视频| 亚洲精品国产av成人精品| 久久久精品免费免费高清| 少妇人妻一区二区三区视频| 18禁在线无遮挡免费观看视频| 成人美女网站在线观看视频| 亚洲av国产av综合av卡| 免费观看无遮挡的男女| 又黄又爽又刺激的免费视频.| 日韩av在线免费看完整版不卡| 免费观看性生交大片5| 久久久国产一区二区| 亚洲欧美清纯卡通| 另类亚洲欧美激情| 国产伦精品一区二区三区视频9| 国产精品99久久99久久久不卡 | 精华霜和精华液先用哪个| 欧美 亚洲 国产 日韩一| 精品人妻一区二区三区麻豆| 久久女婷五月综合色啪小说| 有码 亚洲区| 亚洲精品国产av蜜桃| 人妻夜夜爽99麻豆av| 亚洲欧美日韩卡通动漫| 两个人免费观看高清视频 | 亚洲精品456在线播放app| av网站免费在线观看视频| 久久韩国三级中文字幕| 国产在线男女| 国产成人91sexporn| 99热全是精品| 美女脱内裤让男人舔精品视频| 97在线视频观看| 久久精品国产鲁丝片午夜精品| 国产极品天堂在线| 精品久久久精品久久久| 欧美+日韩+精品| 啦啦啦在线观看免费高清www| 欧美 日韩 精品 国产| 夜夜看夜夜爽夜夜摸| 18禁在线播放成人免费| 观看免费一级毛片| 午夜av观看不卡| 国产精品久久久久久精品电影小说| 亚洲经典国产精华液单| 久久人人爽av亚洲精品天堂| 午夜免费观看性视频| 又大又黄又爽视频免费| 简卡轻食公司| 欧美最新免费一区二区三区| 亚洲国产精品一区三区| 欧美精品高潮呻吟av久久| 日本午夜av视频| 国产精品一区www在线观看| 老司机影院毛片| 视频中文字幕在线观看| 成年女人在线观看亚洲视频| 亚洲av成人精品一二三区| 久久99热6这里只有精品| 男女免费视频国产| 精品少妇内射三级| 亚洲精华国产精华液的使用体验| 国产男人的电影天堂91| 亚洲国产精品国产精品| 深夜a级毛片| 国模一区二区三区四区视频| 亚洲,欧美,日韩| 一级片'在线观看视频| 中文字幕免费在线视频6| 一级,二级,三级黄色视频| 国模一区二区三区四区视频| 午夜福利,免费看| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美日韩综合久久久久久| 美女xxoo啪啪120秒动态图| 午夜影院在线不卡| 国产欧美日韩综合在线一区二区 | 久久久a久久爽久久v久久| 美女主播在线视频| av福利片在线观看| 人人澡人人妻人| 黄色欧美视频在线观看| 亚洲欧洲国产日韩| 色视频www国产| 永久免费av网站大全| 一区二区三区免费毛片| 国产成人一区二区在线| 日韩av免费高清视频| 欧美成人午夜免费资源| 街头女战士在线观看网站| 黄色视频在线播放观看不卡| 欧美日韩综合久久久久久| 人人妻人人添人人爽欧美一区卜| 午夜免费观看性视频| 国产精品人妻久久久影院| 精品亚洲乱码少妇综合久久| 日本午夜av视频| 青春草国产在线视频| 国产成人精品无人区| 大香蕉97超碰在线| 欧美日韩一区二区视频在线观看视频在线| 免费观看无遮挡的男女| 国产一区二区三区综合在线观看 | 久久久精品免费免费高清| 免费久久久久久久精品成人欧美视频 | 水蜜桃什么品种好| 亚洲国产精品一区二区三区在线| 亚洲综合精品二区| 国产有黄有色有爽视频| 国产免费一区二区三区四区乱码| 伦理电影大哥的女人| 中国国产av一级| 两个人免费观看高清视频 | 国产免费又黄又爽又色| 国产精品三级大全| 午夜视频国产福利| av天堂久久9| a级片在线免费高清观看视频| 久久ye,这里只有精品| 亚洲人成网站在线观看播放| 亚洲av电影在线观看一区二区三区| 秋霞伦理黄片| 夫妻午夜视频| 人妻 亚洲 视频| 日韩av免费高清视频| 国产一区有黄有色的免费视频| 国产一区二区在线观看av| 熟妇人妻不卡中文字幕| 国产熟女午夜一区二区三区 | 亚洲av不卡在线观看| 日韩中字成人| 久久久久久久国产电影| 又大又黄又爽视频免费| 免费黄网站久久成人精品| 26uuu在线亚洲综合色| 国产精品伦人一区二区| 国产精品福利在线免费观看| 久久精品久久久久久久性| 国产真实伦视频高清在线观看| 三级国产精品片| 亚洲av中文av极速乱| 免费黄色在线免费观看| 亚洲精品视频女| 久久人人爽人人片av| 日日啪夜夜撸| 麻豆精品久久久久久蜜桃| 国产精品国产三级专区第一集| 久久99热6这里只有精品| 只有这里有精品99| 性色av一级| 国产真实伦视频高清在线观看| 麻豆精品久久久久久蜜桃| 男人舔奶头视频| 99热全是精品| 国产精品偷伦视频观看了| 色视频www国产| 国产熟女午夜一区二区三区 | 在线免费观看不下载黄p国产| 午夜日本视频在线| 日韩一区二区三区影片| 国产精品福利在线免费观看| 欧美3d第一页| 美女主播在线视频| 亚洲不卡免费看| 一级毛片 在线播放| 久久久久久久亚洲中文字幕| 大码成人一级视频| 精品一品国产午夜福利视频| av一本久久久久| 国产伦在线观看视频一区| 国产高清有码在线观看视频| 久久久国产精品麻豆| 一级毛片我不卡| 制服丝袜香蕉在线| 人妻少妇偷人精品九色| 在线观看av片永久免费下载| 少妇熟女欧美另类| 国产 一区精品| 内地一区二区视频在线| 七月丁香在线播放|