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

    聚類在風(fēng)廓線雷達(dá)數(shù)據(jù)特征分析中的應(yīng)用①

    2022-09-20 04:12:14趙玉娟李宗飛李曉波李祥海閆文月
    關(guān)鍵詞:風(fēng)廓均值降水

    趙玉娟, 王 彥, 李宗飛, 李曉波, 郭 陽, 李祥海, 閆文月

    1(天津市氣象信息中心, 天津 300074)

    2(天津市人工影響天氣辦公室, 天津 300074)

    3(天津市氣象探測中心, 天津 300061)

    風(fēng)廓線雷達(dá)是利用大氣湍流對電磁波的散射作用對大氣風(fēng)場進(jìn)行探測的設(shè)備, 可連續(xù)獲得測站上空較高時(shí)間、空間分辨率的垂直風(fēng)廓線資料, 研究其在不同天氣條件下的變化規(guī)律對于改善氣象預(yù)報(bào)服務(wù)質(zhì)量具有重要意義, 因此許多專家學(xué)者開展了風(fēng)廓線雷達(dá)數(shù)據(jù)特征的分析研究.

    楊引明等[1]討論了風(fēng)廓線雷達(dá)的垂直速度和溫度資料在強(qiáng)對流天氣中的應(yīng)用. 鄭石等[2]運(yùn)用風(fēng)廓線雷達(dá)數(shù)據(jù)對南京浦口區(qū)短時(shí)暴雨過程進(jìn)行分析, 發(fā)現(xiàn)了風(fēng)廓線雷達(dá)產(chǎn)品的水平風(fēng)廓線、垂直氣流和折射率結(jié)構(gòu)常數(shù)等在降水過程不同階段的精細(xì)特征. 董保舉等[3]利用云南大理的一次暴雨過程, 研究了風(fēng)廓線雷達(dá)資料在降水前和降水期間的不同特征. 陳楠等[4]利用2005-2008年兩部風(fēng)廓線雷達(dá)風(fēng)場數(shù)據(jù)對南京地區(qū)低空急流的活動(dòng)特征進(jìn)行了統(tǒng)計(jì)分析. 史珺等[5]基于天津地區(qū)一次短時(shí)強(qiáng)降雨過程, 利用3部風(fēng)廓線雷達(dá)資料和降水實(shí)況資料, 對比分析了降水發(fā)生、維持和消亡期間風(fēng)廓線雷達(dá)資料的變化特征, 探討了風(fēng)廓線雷達(dá)對降水天氣的監(jiān)測能力. 王彥等[6]聯(lián)合應(yīng)用5部風(fēng)廓線雷達(dá)、多普勒天氣雷達(dá)組網(wǎng)觀測資料、海河流域自動(dòng)氣象站資料詳細(xì)分析了2012年7月21日海河流域強(qiáng)降水發(fā)生、發(fā)展過程中溫度場、風(fēng)場演變特征.周芯玉等[7]利用風(fēng)廓線雷達(dá)對廣東前汛期短時(shí)強(qiáng)降水類暴雨過程低空風(fēng)場特征進(jìn)行了研究. 王佳津等[8]分析了成都兩次暴雨期間風(fēng)廓線雷達(dá)觀測的低空風(fēng)場特征.

    總體來看, 目前針對不同天氣條件下風(fēng)廓線雷達(dá)數(shù)據(jù)特征的研究多采用的是少量典型天氣個(gè)例和常規(guī)的統(tǒng)計(jì)分析, 利用挖掘技術(shù)的研究較少. 數(shù)據(jù)挖掘是從大量數(shù)據(jù)中發(fā)現(xiàn)有價(jià)值信息的有效手段, 通信、電力、氣象等領(lǐng)域的諸多專家學(xué)者將數(shù)據(jù)挖掘技術(shù)引入實(shí)際應(yīng)用, 取得了不錯(cuò)的效果. 林勤等[9]將雙聚類算法引入電信高價(jià)值客戶細(xì)分, 實(shí)驗(yàn)結(jié)果表明, 可挖掘出更多高價(jià)值客戶群體. 周笑天等[10]、史靜等[11]基于關(guān)聯(lián)規(guī)則、聚類, 提出了氣象數(shù)據(jù)質(zhì)量控制的新方法. 劉偉東等[12]利用K均值聚類方法, 將北京地區(qū)2007-2010年123個(gè)自動(dòng)氣象站夏半年小時(shí)降水分為4個(gè)區(qū), 進(jìn)而分析得到了北京地區(qū)降水的精細(xì)化時(shí)空分布特征.熊亞軍等[13]開展了KNN數(shù)據(jù)挖掘算法在北京地區(qū)霾等級預(yù)報(bào)中的應(yīng)用研究. 杜良敏等[14]提出了基于聚類分區(qū)的中國夏季降水預(yù)測模型. 雷景生等[15]利用模糊C均值算法得到了不同區(qū)域用戶的用電特征. 茍浩峰等[16]基于聚類方法得到了蘭州地區(qū)自動(dòng)站降水的精細(xì)化特征. 李宗飛等[17]基于聚類開展了X波段雙偏振雷達(dá)相態(tài)識(shí)別的研究, 取得了良好的效果. 蘭榮亨等[18]以斗魚直播平臺(tái)為例, 聚焦于直播平臺(tái)上的高消費(fèi)群體,通過構(gòu)建觀眾特征, 采用聚類方法分析高消費(fèi)群體的行為. 任禹丞等[19]將聚類分析技術(shù)應(yīng)用在電力客戶用電問題分析領(lǐng)域, 挖掘出了隱藏的客戶用電問題關(guān)鍵信息, 為改進(jìn)電力客服質(zhì)量與潛在服務(wù)風(fēng)險(xiǎn)預(yù)測提供了數(shù)據(jù)支撐.

    天津靜海邊界層風(fēng)廓線雷達(dá)站已積累了大量的觀測資料, 探索基于數(shù)據(jù)挖掘技術(shù)的風(fēng)廓線雷達(dá)數(shù)據(jù)特征分析方法, 揭示風(fēng)廓線雷達(dá)數(shù)據(jù)更具普遍性的規(guī)律,為氣象預(yù)報(bào)預(yù)測、防災(zāi)減災(zāi)提供更多參考, 具有非常重要的意義.

    本文主要研究了將聚類分析技術(shù)應(yīng)用在風(fēng)廓線雷達(dá)數(shù)據(jù)特征分析領(lǐng)域, 通過數(shù)據(jù)預(yù)處理及改進(jìn)的聚類方法, 對風(fēng)廓線雷達(dá)觀測數(shù)據(jù)進(jìn)行分析. 文中基于聚類技術(shù), 以天津靜海風(fēng)廓線雷達(dá)觀測數(shù)據(jù)為基礎(chǔ), 提出了針對風(fēng)廓線雷達(dá)數(shù)據(jù)分析的改進(jìn)聚類算法, 并通過實(shí)驗(yàn)驗(yàn)證了該方法可快速準(zhǔn)確地實(shí)現(xiàn)風(fēng)廓線雷達(dá)數(shù)據(jù)的自動(dòng)聚類, 挖掘出風(fēng)廓線雷達(dá)數(shù)據(jù)中更具普遍性的規(guī)律, 從而為風(fēng)廓線雷達(dá)數(shù)據(jù)特征分析提供了新思路, 為氣象預(yù)報(bào)服務(wù)提供了新參考依據(jù).

    1 聚類模型設(shè)計(jì)

    1.1 K-means算法原理及適用性分析

    聚類包括基于劃分、基于層次、基于密度、基于網(wǎng)格等多種方法, 本文采用基于劃分的聚類方法中較為經(jīng)典的K-means算法. K-means算法把類的形心定義為類的均值, 對包含N個(gè)對象的數(shù)據(jù)集D的聚類處理流程如下:

    Step 1. 從數(shù)據(jù)集中隨機(jī)選取K個(gè)樣本作為初始聚類中心C={c1, c2, …, ck}.

    Step 2. 針對數(shù)據(jù)集中每個(gè)樣本xi, 計(jì)算它到K個(gè)聚類中心的距離并將其分到距離最小的聚類中心所對應(yīng)的類中.

    Step 3. 針對每個(gè)類別ci, 重新計(jì)算它的聚類中心.

    Step 4. 重復(fù)第2步和第3步直到聚類中心不再變化.

    K-means方法適合處理數(shù)值型數(shù)據(jù)集, 可對大量數(shù)據(jù)根據(jù)屬性進(jìn)行高效分類. 風(fēng)廓線雷達(dá)數(shù)據(jù)的探測高度、垂直速度屬性均為數(shù)值型. 邊界層風(fēng)廓線雷達(dá)工作在L波段, 在無降水天氣條件下, 接收的信號(hào)主要是大氣湍流運(yùn)動(dòng)造成折射率分布不均勻而產(chǎn)生的散射;在降水天氣條件下, 接收的信號(hào)主要是雨滴在入射電磁波極化下作強(qiáng)迫的多極震蕩而產(chǎn)生的散射[20,21]. 不同天氣條件下, 風(fēng)廓線雷達(dá)探測能力差別很大. 風(fēng)廓線雷達(dá)的探測原理決定了其在不同天氣條件下的探測數(shù)據(jù)會(huì)存在較大差異, 已有個(gè)例研究也表明, 風(fēng)廓線雷達(dá)數(shù)據(jù)的探測高度、垂直速度在降水、晴天等天氣條件存在不同規(guī)律, 因此, 可將K-means方法應(yīng)用于風(fēng)廓線雷達(dá)數(shù)據(jù)探測高度、垂直速度屬性的聚類劃分.

    1.2 基于K-means改進(jìn)的風(fēng)廓線雷達(dá)數(shù)據(jù)聚類算法

    本文選用K-means算法對風(fēng)廓線雷達(dá)數(shù)據(jù)進(jìn)行聚類分析, 考慮到傳統(tǒng)K-means算法的初始聚類中心隨機(jī)生成, 容易陷入局部最優(yōu), 本文針對初始聚類中心的選擇進(jìn)行了改進(jìn), 采用K-means++算法初始化聚類中心, 具體步驟如下:

    (1) 從輸入樣本中隨機(jī)選擇一個(gè)樣本作為初始聚類中心Ur.

    (2) 首先計(jì)算每個(gè)樣本與當(dāng)前已有聚類中心間的最短距離(即與最近一個(gè)聚類中心的距離), 用D(x)表示; 接著計(jì)算每個(gè)樣本被選為下一個(gè)聚類中心的概率P, P的計(jì)算如式(1)所示. 最后, 按照輪盤法選擇出下一個(gè)聚類中心.

    (3) 重復(fù)第2步直到選擇出K個(gè)聚類中心.

    1.3 數(shù)據(jù)規(guī)范化方法

    對于基于距離度量的聚類方法, 數(shù)據(jù)標(biāo)準(zhǔn)化可以賦予所有屬性相等的權(quán)重, 避免值域較大的屬性與較小值域的屬性相比權(quán)重過大, 保障聚類效果, 包括最小-最大規(guī)范化、Z-score規(guī)范化、按小數(shù)定標(biāo)規(guī)范化等方法, 本文選用Z-score規(guī)范化方法對聚類指標(biāo)進(jìn)行預(yù)處理, Z-score規(guī)范化方法公式如下.

    首先計(jì)算均值絕對偏差和平均值:

    其中, Mf為第f個(gè)屬性的平均值; Sf為第f個(gè)屬性的均值絕對偏差; Xif表示第i條記錄的第f個(gè)屬性.

    然后計(jì)算標(biāo)準(zhǔn)化度量值, Zf為標(biāo)準(zhǔn)化后的第f個(gè)屬性值, 計(jì)算公式為:

    1.4 聚類最優(yōu)K值選取方法

    對于K-means算法, 聚類數(shù)K的選取非常重要. 為避免K值選取上的主觀性, 本研究運(yùn)用肘部法則和輪廓系數(shù)兩種方法確定最優(yōu)K值, 這兩項(xiàng)指標(biāo)綜合考慮了聚類結(jié)構(gòu)中不同類別樣本的離散性和同類別樣本的凝聚性.

    肘部法則(elbow method)基于如下觀察: 增加類個(gè)數(shù)有助于降低每個(gè)類的類內(nèi)誤差平方和(sum of the squared errors, SSE), 這是因?yàn)橛懈嗟念惪梢圆东@更細(xì)的數(shù)據(jù)對象類, 類內(nèi)對象間更為相似. 然而, 如果形成太多的類, 則降低SSE的邊際效應(yīng)可能下降, 因?yàn)榘岩粋€(gè)凝聚的類分裂成兩個(gè)引起SSE的降幅較小. 繪制SSE關(guān)于聚類數(shù)K的曲線, 曲線的第一個(gè)(或最顯著的)拐點(diǎn)即手肘處對應(yīng)最佳聚類數(shù).

    對于包含N個(gè)對象的數(shù)據(jù)集D, 假設(shè)D被劃分為K個(gè)類C1, C2, …, Ck. 對于每個(gè)對象o∈D, 其輪廓系數(shù)(sihouette coefficient) s(o)的定義如下:

    其中, a(o)代表對象o與o所屬類的其他對象間的平均距離, b(o)代表o與除自身所在類外對象的最小平均距離. a(o)反映了o所屬類的緊湊性, b(o)反映了o與其他簇的分離程度.

    由定義可知, 輪廓系數(shù)值在-1和1之間. 我們計(jì)算數(shù)據(jù)集中所有對象輪廓系數(shù)的平均值, 以評估聚類質(zhì)量. 輪廓系數(shù)值越接近1, 代表聚類效果越好.

    2 風(fēng)廓線雷達(dá)數(shù)據(jù)聚類及結(jié)果分析

    2.1 總體流程

    本研究針對風(fēng)廓線雷達(dá)數(shù)據(jù)聚類分析的步驟主要包括: (1)構(gòu)建風(fēng)廓線雷達(dá)觀測信息數(shù)據(jù)庫: 基于風(fēng)廓線雷達(dá)觀測的小時(shí)數(shù)據(jù)文件, 解析提取其中不同高度的觀測數(shù)據(jù), 建立風(fēng)廓線雷達(dá)觀測信息數(shù)據(jù)庫. (2)數(shù)據(jù)預(yù)處理: 對存儲(chǔ)在風(fēng)廓線雷達(dá)觀測信息數(shù)據(jù)庫中的數(shù)據(jù), 首先剔除不合理數(shù)據(jù), 并按時(shí)次統(tǒng)計(jì)得到逐時(shí)最大探測高度、最大垂直速度; 其次將最大探測高度、最大垂直速度進(jìn)行規(guī)范化處理, 得到聚類指標(biāo). (3)聚類分析及特征提取: 利用改進(jìn)后的K-means算法, 對預(yù)處理后的風(fēng)廓線雷達(dá)數(shù)據(jù)進(jìn)行聚類, 根據(jù)聚類結(jié)果、自動(dòng)站觀測的降水及人工觀測的云量、天氣現(xiàn)象等數(shù)據(jù)歸納總結(jié)出不同天氣條件下風(fēng)廓線雷達(dá)數(shù)據(jù)特征.

    2.2 數(shù)據(jù)來源及預(yù)處理

    研究選取天津靜海的風(fēng)廓線雷達(dá)觀測數(shù)據(jù), 同時(shí)結(jié)合自動(dòng)站觀測的降水及人工觀測的天氣現(xiàn)象、云量數(shù)據(jù)標(biāo)記各時(shí)次天氣條件. 天津靜海風(fēng)廓線雷達(dá)為CFL-03B型邊界層風(fēng)廓線雷達(dá), 雷達(dá)基本技術(shù)性能如表1所示. 靜海自動(dòng)站與風(fēng)廓線雷達(dá)站距離在1 km以內(nèi), 自動(dòng)站降水量數(shù)據(jù)可代表風(fēng)廓線雷達(dá)上空降水情況, 風(fēng)廓線雷達(dá)及自動(dòng)站皆選用1 h頻次數(shù)據(jù). 天津地區(qū)降水主要發(fā)生在5-9月, 故研究時(shí)段選取2017-2019年5-9月. 由于靜海站相關(guān)時(shí)段無云量人工觀測, 因此云量數(shù)據(jù)選用距靜海站最近的西青站人工觀測數(shù)據(jù).

    表1 風(fēng)廓線雷達(dá)基本性能參數(shù)

    靜海風(fēng)廓線雷達(dá)小時(shí)觀測數(shù)據(jù)包括66個(gè)高度層,2017-2019年5-9月小時(shí)產(chǎn)品數(shù)據(jù)文件解析入庫后得到逐時(shí)、分高度層觀測數(shù)據(jù)共計(jì)697 008條. 通過剔除無效數(shù)據(jù), 對風(fēng)廓線雷達(dá)小時(shí)觀測數(shù)據(jù)按時(shí)次統(tǒng)計(jì)得到10 553條記錄. 本文所使用的垂直速度是風(fēng)廓線雷達(dá)探測到的垂直速度(定義垂直速度向下為正, 向上為負(fù)), 未經(jīng)訂正, 所以降水時(shí)代表了空氣的垂直運(yùn)動(dòng)和降水粒子下沉運(yùn)動(dòng)的總和, 探測高度指距離地面的高度.

    按照GB/T 35663-2017中定義的天氣預(yù)報(bào)基本術(shù)語, 根據(jù)自動(dòng)站觀測的過去1 h降水量及人工觀測的天氣現(xiàn)象、云量, 將各時(shí)次劃分為晴、少云、多云、陰、降水形成(降水發(fā)生前1-6 h)、降水期間、降水結(jié)束(降水結(jié)束后1-3 h) 7種狀態(tài), 詳細(xì)劃分規(guī)則如下:

    (1)若當(dāng)前時(shí)次自動(dòng)站觀測的過去1 h降水量≥0.1 mm或人工觀測的天氣現(xiàn)象記錄有降水, 標(biāo)記為降水期間; 否則, 轉(zhuǎn)(2).

    (2)當(dāng)前時(shí)次未來1至6 h為降水階段, 標(biāo)記為降水形成階段. 否則, 轉(zhuǎn)(3).

    (3)當(dāng)前時(shí)次前推1至3 h為降水階段, 標(biāo)記為降水結(jié)束階段. 否則, 轉(zhuǎn)(4).

    (4)根據(jù)人工觀測的日均總云量標(biāo)記天氣條件. 日均總云量0-2成, 標(biāo)記為晴; 3-5成, 標(biāo)記為少云;6-8成, 標(biāo)記為多云; 9-10成標(biāo)記為陰.

    2.3 最佳聚類數(shù)K值確定

    輪廓系數(shù)值越大聚類效果越好, 肘部法則中理論上最佳 K 值在肘處取得. 聚類數(shù)K的范圍設(shè)定為[2, 9],輪廓系數(shù)及聚類誤差平方和(sum of the squared errors,SSE)隨聚類個(gè)數(shù)K的變化如圖1所示. 由圖1可知,聚類數(shù)K與誤差平方和的拐點(diǎn)出現(xiàn)在K=3處, 此時(shí)輪廓系數(shù)也較大, 綜合考慮輪廓系數(shù)、肘部法則, 確定聚類數(shù)K為3.

    圖1 不同K值的SSE及輪廓系數(shù)變化

    2.4 傳統(tǒng)K-means算法與改進(jìn)算法對比分析

    分別利用原K-means算法和改進(jìn)后的K-means算法對風(fēng)廓線雷達(dá)數(shù)據(jù)進(jìn)行聚類, K值選取最佳聚類數(shù)3, 每種算法運(yùn)行10次, 選取平均值進(jìn)行對比, 對比指標(biāo)包括輪廓系數(shù)、迭代次數(shù)、運(yùn)行時(shí)間3方面, 對比結(jié)果如表2所示. 由對比結(jié)果可見, 改進(jìn)后的算法在迭代次數(shù)和運(yùn)行時(shí)間均明顯降低的情況下, 可獲得與原K-means算法同等聚類效果.

    表2 K-means和改進(jìn)后的K-means算法聚類對比

    2.5 基于聚類的風(fēng)廓線雷達(dá)數(shù)據(jù)特征分析

    應(yīng)用改進(jìn)后的K-means算法對10 553個(gè)時(shí)次風(fēng)廓線雷達(dá)數(shù)據(jù)進(jìn)行聚類, K值選取最佳聚類數(shù)3, 聚類結(jié)果如圖2所示, 計(jì)算每類各特征統(tǒng)計(jì)描述, 結(jié)果如表3所示. 根據(jù)表3可知, class0包含7 669個(gè)時(shí)次數(shù)據(jù), 占比72.7%. class1包含615個(gè)時(shí)次數(shù)據(jù), 占比僅5.8%.class2包含2 269個(gè)時(shí)次數(shù)據(jù), class0的最大探測高度、最大垂直速度均最低, 最大探測高度均值低于3 500 m, 最大值低于5 000 m, 最大垂直速度均值在0.5 m/s. class1的最大探測高度和最大垂直速度皆最高, 最大探測高度均值為6 718 m, 該類近70%的時(shí)次最大探測高度可達(dá)靜海風(fēng)廓線雷達(dá)的最大探測高度7 080 m, 最大垂直速度均值在5.5 m/s. class2的最大探測高度和最大垂直速度介于class0和class1之間, 最大探測高度遠(yuǎn)高于class0, 均值亦在6 000 m以上, 但達(dá)到7 080 m的僅占38%, 最大垂直速度均值略高于class0, 為1 m/s, 明顯低于class1.

    表3 各類最大探測高度、最大垂直速度特征

    圖2 聚為3類的風(fēng)廓線雷達(dá)數(shù)據(jù)分布

    經(jīng)分析, class0主要為晴、少云天氣觀測數(shù)據(jù),class1主要為降水期間的觀測數(shù)據(jù). class2主要為多云、陰天及降水前1-6 h、降水后1-3 h的觀測數(shù)據(jù).

    探測高度方面, 晴天、少云天氣的最大探測高度均值在3 562 m, 多云、陰天的最大探測高度均值為4 414 m, 降水前1-6 h、降水后1-3 h的最大探測高度均值為5 480 m, 降水期間的最大探測高度均值為6 548 m. 多云、陰天探測高度高于晴天及少云天氣, 降水前1-6 h、降水后1-3 h風(fēng)廓線雷達(dá)的探測高度較多云、陰天更高, 降水期間探測高度最高且多數(shù)可達(dá)靜海風(fēng)廓線雷達(dá)的最大探測高度. 風(fēng)廓線雷達(dá)主要以晴空大氣作為探測對象, 利用大氣湍流對電磁波的散射作用來探測大氣風(fēng)場等要素. 其探測高度除了和雷達(dá)技術(shù)指標(biāo)有關(guān)之外, 受天氣狀況的影響很大[22], 當(dāng)有天氣系統(tǒng)的時(shí)候, 大氣湍流活動(dòng)加強(qiáng), 高層的水汽含量增加, 雷達(dá)探測高度會(huì)增加[23], 聚類結(jié)果表明風(fēng)廓線雷達(dá)的最大探測高度對于大氣高層水汽含量增大有較好的指示作用.

    垂直速度方面, 晴天、少云天氣的最大垂直速度均值為0.5 m/s, 多云、陰天的最大垂直速度均值為0.7 m/s,降水前1-6 h、降水后1-3 h的最大垂直速度均值為1.4 m/s, 降水期間數(shù)據(jù)最大垂直速度均值為5.5 m/s. 降水期間數(shù)據(jù)最大垂直速度大于等于4 m/s的占74%, 而其他天氣條件97%的時(shí)次最大垂直速度低于4 m/s. 楊引明等[1]分析發(fā)現(xiàn), 4 m/s的垂直速度反映了降水開始和結(jié)束(定義垂直速度向下為正, 向上為負(fù)), 本文研究結(jié)果與其研究結(jié)論基本一致.

    總體來說, 風(fēng)廓線雷達(dá)的最大探測高度在天氣狀況由晴天、少云轉(zhuǎn)變?yōu)槎嘣啤㈥幪?、降水天氣時(shí), 會(huì)有明顯的上升. 但最大垂直速度僅在降水開始時(shí)急劇增大, 在降水開始前增長趨勢不明顯.

    兩個(gè)典型天氣個(gè)例也印證了本文的挖掘結(jié)果. 圖3展示了2018年7月23日00時(shí)至2018年7月26日05時(shí)靜海站最大探測高度、最大垂直速度與小時(shí)雨量的時(shí)序變化對比. 23日屬于少云天氣, 24日00時(shí)至12時(shí), 靜海本站出現(xiàn)大暴雨, 累計(jì)降水量153.6 mm, 最大小時(shí)降水量27.1 mm.

    圖3 2018年7月23日至24日靜海風(fēng)廓線雷達(dá)最大探測高度、最大垂直速度與小時(shí)雨量對比圖

    由圖3可見, 降水前24 h到降水開始, 探測高度不斷上升. 降水前24-12 h最大探測高度均值在3.5 km左右, 降水前12-6 h最大探測高度均值升至 4.5 km, 之后探測高度有所下降, 23日22時(shí)后又開始升高, 到24日01時(shí)升至靜海風(fēng)廓線雷達(dá)最大探測高度7 km,3 h增幅達(dá)4.3 km, 強(qiáng)降水期間一直維持在最大探測高度. 說明隨著降水臨近, 大氣的湍流活動(dòng)明顯加強(qiáng)、大氣高層的水汽含量急劇增加. 降水結(jié)束后, 由于屬于多云天氣, 最大探測高度在4.8 km以上維持了一段時(shí)間,直到25日05時(shí)后, 才降低到3.5 km以下. 降水期間,最大垂直速度均在4 m/s以上, 降水結(jié)束后最大垂直速度迅速降低.

    圖4展示了2018年8月18日00時(shí)至2018年8月20日12時(shí)靜海風(fēng)廓線雷達(dá)在陰天及地面降水前后最大探測高度、最大垂直速度的時(shí)序變化.

    圖4 2018年8月18日至20日靜海風(fēng)廓線雷達(dá)最大探測高度、最大垂直速度與小時(shí)雨量對比圖

    8月18日靜海為陰天并在15時(shí)至16時(shí)發(fā)生1.4 mm小雨級別降水, 最大探測高度基本處于5 km以上,18日12時(shí)降至4.5 km后又在降水前1 h的14時(shí)再次升至6.5 km, 2 h躍升2 km, 18日20時(shí)后逐漸下降. 探測高度從19日02時(shí)的3.7 km又開始升高, 到19日04時(shí)降水前1 h升至最大探測高度7 km, 19日05時(shí)至22時(shí)靜海再次發(fā)生17 mm中雨級別降水, 降水期間探測高度基本維持在7 km, 20日03時(shí)之后開始逐漸下降. 從探測高度曲線圖可以看出, 降水前6 h最大探測高度基本在5 km及以上, 即使有降低至4 km以下的時(shí)次, 也會(huì)在降水前1 h回升, 這可能是地面降水前高空先有水汽到達(dá)的緣故, 降水結(jié)束后, 探測高度沒有立即降低, 在降水結(jié)束3 h或更長時(shí)間后才逐漸降至4 km以下, 可能是高空還有一定潮濕空氣的緣故, 降水期間最大探測高度多在7 080 m. 對比最大垂直速度和小時(shí)降水量曲線圖可知, 降水期間最大垂直速度多在4 m/s以上, 降水強(qiáng)度大時(shí)垂直速度較大, 但最大垂直速度值不一定出現(xiàn)在降雨最強(qiáng)時(shí)次. 垂直速度在降水開始會(huì)迅速增大, 降水結(jié)束會(huì)快速下降.

    3 總結(jié)與展望

    針對風(fēng)廓線雷達(dá)觀測數(shù)據(jù)特點(diǎn), 本文利用改進(jìn)的K-means聚類算法建立了風(fēng)廓線雷達(dá)數(shù)據(jù)特征聚類分析模型, 對天津靜海邊界層風(fēng)廓線雷達(dá)2017-2019年(5-9月)的觀測資料進(jìn)行聚類, 挖掘出了風(fēng)廓線雷達(dá)的最大探測高度、最大垂直速度在晴、陰天、降水前1-6 h、降水期間、降水后1-3 h等天氣條件下的一些變化規(guī)律, 為氣象預(yù)報(bào)服務(wù)提供了更具普遍性的參考依據(jù), 為風(fēng)廓線雷達(dá)數(shù)據(jù)特征分析提供了新思路. 但本文僅對探測高度、垂直速度特征進(jìn)行了分析, 風(fēng)廓線雷達(dá)還有水平風(fēng)向、水平風(fēng)速、折射率結(jié)構(gòu)常數(shù)等觀測要素, 下一步有必要對風(fēng)廓線雷達(dá)其他觀測信息繼續(xù)開展挖掘分析, 以獲取更多的數(shù)據(jù)特征為氣象預(yù)報(bào)服務(wù)提供支撐.

    猜你喜歡
    風(fēng)廓均值降水
    高郵邊界層風(fēng)廓線雷達(dá)數(shù)據(jù)獲取率分析
    一種綜合的風(fēng)廓線雷達(dá)數(shù)據(jù)質(zhì)量控制方法
    用L波段探空測風(fēng)雷達(dá)評估風(fēng)廓線雷達(dá)測風(fēng)準(zhǔn)確性
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    四川盆地風(fēng)廓線雷達(dá)大氣折射率結(jié)構(gòu)常數(shù)特征分析
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    降水現(xiàn)象儀模擬軟件設(shè)計(jì)與實(shí)現(xiàn)
    均值不等式失效時(shí)的解決方法
    均值與方差在生活中的應(yīng)用
    關(guān)于均值有界變差函數(shù)的重要不等式
    国内精品久久久久精免费| 身体一侧抽搐| 国产精品人妻久久久久久| 欧美日韩亚洲国产一区二区在线观看| 97热精品久久久久久| 精品一区二区三区视频在线| 国产v大片淫在线免费观看| 亚洲第一电影网av| av专区在线播放| 一本综合久久免费| 国产老妇女一区| 国产乱人视频| 91字幕亚洲| 国产91精品成人一区二区三区| 我要看日韩黄色一级片| 国产精品久久久久久人妻精品电影| 一级作爱视频免费观看| 简卡轻食公司| 一区二区三区免费毛片| 欧美成人一区二区免费高清观看| 最后的刺客免费高清国语| 欧美色视频一区免费| 国产精华一区二区三区| 一本久久中文字幕| 在线观看午夜福利视频| 一本一本综合久久| 很黄的视频免费| 女人被狂操c到高潮| 国内精品久久久久久久电影| 97超级碰碰碰精品色视频在线观看| 在线免费观看不下载黄p国产 | 91午夜精品亚洲一区二区三区 | 亚洲av五月六月丁香网| 18+在线观看网站| 欧美绝顶高潮抽搐喷水| 757午夜福利合集在线观看| 一本一本综合久久| 最近中文字幕高清免费大全6 | 久久亚洲精品不卡| 婷婷六月久久综合丁香| 久久久成人免费电影| 国内揄拍国产精品人妻在线| 91麻豆精品激情在线观看国产| 亚洲av美国av| 亚洲最大成人av| 国内毛片毛片毛片毛片毛片| 久久亚洲真实| 夜夜爽天天搞| 97热精品久久久久久| 国产成人影院久久av| 亚洲国产欧美人成| a级毛片a级免费在线| 乱人视频在线观看| 成人午夜高清在线视频| 亚洲精品久久国产高清桃花| 国产亚洲av嫩草精品影院| 国产成人av教育| 国产成人欧美在线观看| 国产高清三级在线| 国产精品野战在线观看| 日日摸夜夜添夜夜添av毛片 | 一区二区三区激情视频| 12—13女人毛片做爰片一| 国产男靠女视频免费网站| 最好的美女福利视频网| 热99在线观看视频| 观看美女的网站| 我要看日韩黄色一级片| 色播亚洲综合网| 亚洲精品久久国产高清桃花| av女优亚洲男人天堂| 很黄的视频免费| 在线看三级毛片| 久久精品国产99精品国产亚洲性色| 极品教师在线视频| 精品久久久久久成人av| 脱女人内裤的视频| 久久久久久久久中文| 熟妇人妻久久中文字幕3abv| 又黄又爽又免费观看的视频| 少妇高潮的动态图| 成人av一区二区三区在线看| 日本五十路高清| 一a级毛片在线观看| 欧美黄色片欧美黄色片| 色综合站精品国产| 搡老妇女老女人老熟妇| 欧美日韩乱码在线| 久久久久久久精品吃奶| 丝袜美腿在线中文| 国产亚洲精品久久久com| 午夜激情福利司机影院| 亚洲av电影不卡..在线观看| 高清在线国产一区| 可以在线观看毛片的网站| 日韩欧美国产在线观看| 好男人在线观看高清免费视频| 18禁黄网站禁片午夜丰满| 如何舔出高潮| av在线老鸭窝| 国产激情偷乱视频一区二区| 亚洲中文日韩欧美视频| 搞女人的毛片| 黄色日韩在线| 国产精品1区2区在线观看.| 麻豆久久精品国产亚洲av| 少妇人妻一区二区三区视频| 黄色女人牲交| 女生性感内裤真人,穿戴方法视频| a级毛片免费高清观看在线播放| 成人午夜高清在线视频| 亚洲人成电影免费在线| 欧美成人免费av一区二区三区| 久久久久亚洲av毛片大全| 男人舔女人下体高潮全视频| 久久这里只有精品中国| 天堂√8在线中文| 99国产综合亚洲精品| 夜夜爽天天搞| 最近在线观看免费完整版| 久久6这里有精品| a级毛片免费高清观看在线播放| 少妇裸体淫交视频免费看高清| 女人十人毛片免费观看3o分钟| 久久午夜亚洲精品久久| 成人性生交大片免费视频hd| 日本黄色片子视频| 国产精品一区二区免费欧美| 久久久久久久久中文| 亚洲经典国产精华液单 | 亚洲国产精品成人综合色| 亚洲五月婷婷丁香| 免费高清视频大片| 免费电影在线观看免费观看| 欧美日韩中文字幕国产精品一区二区三区| 欧美一区二区亚洲| 欧美黑人巨大hd| 一本一本综合久久| 色综合欧美亚洲国产小说| 精品一区二区三区av网在线观看| 色5月婷婷丁香| 欧美三级亚洲精品| 91av网一区二区| 床上黄色一级片| 老鸭窝网址在线观看| 一级作爱视频免费观看| 久久精品国产亚洲av天美| 成人欧美大片| 国语自产精品视频在线第100页| 国产高潮美女av| eeuss影院久久| 变态另类丝袜制服| 男女视频在线观看网站免费| 国产欧美日韩一区二区精品| 久久久国产成人精品二区| 直男gayav资源| 免费一级毛片在线播放高清视频| 伊人久久精品亚洲午夜| 亚洲电影在线观看av| 精品久久久久久,| www.999成人在线观看| 久久6这里有精品| 99热6这里只有精品| 色尼玛亚洲综合影院| 亚洲成av人片在线播放无| 99视频精品全部免费 在线| 日韩欧美精品v在线| 久久精品影院6| 性插视频无遮挡在线免费观看| 国产欧美日韩一区二区精品| 国产精品久久电影中文字幕| 99久久精品一区二区三区| 国产精品一区二区三区四区免费观看 | 两性午夜刺激爽爽歪歪视频在线观看| 国产免费一级a男人的天堂| avwww免费| 国产精品日韩av在线免费观看| 国产精品久久视频播放| 成人性生交大片免费视频hd| 国产综合懂色| 国产精品亚洲美女久久久| 免费在线观看影片大全网站| 男人舔女人下体高潮全视频| 久久久久久大精品| 亚洲最大成人手机在线| 欧美日韩亚洲国产一区二区在线观看| 窝窝影院91人妻| 午夜两性在线视频| 婷婷色综合大香蕉| 精品久久久久久久久av| 久久久久国内视频| 精品无人区乱码1区二区| 午夜久久久久精精品| 亚洲最大成人中文| 又粗又爽又猛毛片免费看| 色吧在线观看| 欧美日韩瑟瑟在线播放| 在线播放国产精品三级| 国产色婷婷99| 国产 一区 欧美 日韩| 少妇熟女aⅴ在线视频| 欧美+亚洲+日韩+国产| 我的女老师完整版在线观看| 岛国在线免费视频观看| 非洲黑人性xxxx精品又粗又长| 九色国产91popny在线| 精品福利观看| 一级a爱片免费观看的视频| 欧美一区二区国产精品久久精品| 最近在线观看免费完整版| 亚洲美女搞黄在线观看 | 精品久久国产蜜桃| a在线观看视频网站| 又爽又黄a免费视频| 亚洲专区国产一区二区| 亚洲美女黄片视频| 亚洲自偷自拍三级| 夜夜爽天天搞| 久久精品国产99精品国产亚洲性色| 免费一级毛片在线播放高清视频| 亚洲一区二区三区色噜噜| 中文字幕人成人乱码亚洲影| 亚洲av电影在线进入| av女优亚洲男人天堂| 人妻丰满熟妇av一区二区三区| 国产精品1区2区在线观看.| 久久精品国产自在天天线| 超碰av人人做人人爽久久| 在线免费观看不下载黄p国产 | 亚洲精品乱码久久久v下载方式| 特级一级黄色大片| 日韩免费av在线播放| 在线免费观看的www视频| 最后的刺客免费高清国语| 成人亚洲精品av一区二区| 特级一级黄色大片| 亚洲精品日韩av片在线观看| 午夜亚洲福利在线播放| 欧美色欧美亚洲另类二区| 淫妇啪啪啪对白视频| 在线观看舔阴道视频| 国产欧美日韩精品一区二区| 亚洲最大成人手机在线| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品综合久久久久久久免费| 精品免费久久久久久久清纯| 久久精品国产亚洲av香蕉五月| 99热只有精品国产| 99riav亚洲国产免费| 国产一区二区激情短视频| 搡老妇女老女人老熟妇| 久久久久性生活片| 国产精品人妻久久久久久| 精品人妻偷拍中文字幕| 亚洲成人精品中文字幕电影| 我的女老师完整版在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲av电影不卡..在线观看| 蜜桃久久精品国产亚洲av| 欧美日本视频| 国产精品98久久久久久宅男小说| 亚洲精品亚洲一区二区| 国产精品三级大全| 国产真实乱freesex| 国产成人a区在线观看| 小说图片视频综合网站| 99久久精品热视频| 两个人视频免费观看高清| 男人狂女人下面高潮的视频| 国产探花极品一区二区| 无遮挡黄片免费观看| 欧美日韩乱码在线| 首页视频小说图片口味搜索| 脱女人内裤的视频| 日本黄大片高清| 宅男免费午夜| a在线观看视频网站| 九色国产91popny在线| 999久久久精品免费观看国产| 97超视频在线观看视频| 久久九九热精品免费| 国产在线精品亚洲第一网站| 免费av不卡在线播放| 久久亚洲精品不卡| 亚洲黑人精品在线| 国产久久久一区二区三区| 亚洲av五月六月丁香网| 日韩成人在线观看一区二区三区| 我要搜黄色片| 欧美3d第一页| 欧美日韩国产亚洲二区| 亚洲人与动物交配视频| 又黄又爽又刺激的免费视频.| 国产精品久久久久久久电影| 国产精品一及| 在现免费观看毛片| 国产aⅴ精品一区二区三区波| 熟女电影av网| 精品一区二区免费观看| 一个人免费在线观看电影| 美女高潮的动态| 老司机深夜福利视频在线观看| 免费看a级黄色片| 观看美女的网站| 特大巨黑吊av在线直播| 欧美不卡视频在线免费观看| 热99在线观看视频| 欧美一区二区国产精品久久精品| 亚洲av熟女| 久久伊人香网站| 国产精品久久久久久久电影| 天堂动漫精品| 一级作爱视频免费观看| 精品国内亚洲2022精品成人| 日本黄大片高清| 久久精品影院6| 午夜福利视频1000在线观看| 中文字幕久久专区| 老熟妇仑乱视频hdxx| 国产欧美日韩精品亚洲av| 99在线视频只有这里精品首页| 中文在线观看免费www的网站| 国产老妇女一区| 日本与韩国留学比较| 欧美日韩乱码在线| 极品教师在线视频| 婷婷精品国产亚洲av| 一级a爱片免费观看的视频| 天天躁日日操中文字幕| 亚洲在线观看片| 国产亚洲精品久久久久久毛片| 亚洲综合色惰| 久久久国产成人免费| 高潮久久久久久久久久久不卡| 免费在线观看影片大全网站| 天堂动漫精品| 三级毛片av免费| 日日摸夜夜添夜夜添av毛片 | 成年版毛片免费区| 久久精品久久久久久噜噜老黄 | 黄色配什么色好看| 成年版毛片免费区| 久久久久久久久大av| 国产精华一区二区三区| 丁香六月欧美| 亚洲国产精品999在线| 亚洲av第一区精品v没综合| 蜜桃久久精品国产亚洲av| 欧美区成人在线视频| 舔av片在线| 欧美极品一区二区三区四区| av福利片在线观看| 亚州av有码| 少妇人妻精品综合一区二区 | 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 亚洲自偷自拍三级| 久99久视频精品免费| 国产精品电影一区二区三区| 亚洲自拍偷在线| 啪啪无遮挡十八禁网站| 级片在线观看| 老熟妇乱子伦视频在线观看| 婷婷色综合大香蕉| 亚洲成av人片免费观看| 麻豆一二三区av精品| 久久国产乱子伦精品免费另类| 亚洲五月天丁香| www.熟女人妻精品国产| 国内精品美女久久久久久| a级毛片a级免费在线| 首页视频小说图片口味搜索| 国产精品98久久久久久宅男小说| 亚洲精品色激情综合| 无人区码免费观看不卡| 日韩国内少妇激情av| 国产精品一及| 久久久久性生活片| 网址你懂的国产日韩在线| 精品久久国产蜜桃| xxxwww97欧美| 免费无遮挡裸体视频| 国产男靠女视频免费网站| 国产一区二区在线av高清观看| 国产视频内射| 欧美最新免费一区二区三区 | 乱人视频在线观看| 99热这里只有是精品在线观看 | 欧美zozozo另类| 国产色爽女视频免费观看| 久久欧美精品欧美久久欧美| 亚洲av第一区精品v没综合| 久久人妻av系列| 午夜福利在线在线| 91久久精品国产一区二区成人| 观看免费一级毛片| 日韩欧美国产一区二区入口| 女人十人毛片免费观看3o分钟| 夜夜躁狠狠躁天天躁| 亚洲自拍偷在线| 国产精华一区二区三区| 国产午夜精品论理片| 日韩欧美一区二区三区在线观看| 日本一二三区视频观看| 在线播放无遮挡| 在线免费观看的www视频| 国产中年淑女户外野战色| 国产精品亚洲一级av第二区| 亚洲内射少妇av| 麻豆av噜噜一区二区三区| 国产亚洲精品久久久久久毛片| av在线天堂中文字幕| 中文在线观看免费www的网站| 自拍偷自拍亚洲精品老妇| 国产成人aa在线观看| 国产精品1区2区在线观看.| 欧美性猛交╳xxx乱大交人| 美女被艹到高潮喷水动态| 亚洲av美国av| 亚洲第一电影网av| 97碰自拍视频| 欧美xxxx黑人xx丫x性爽| 午夜精品一区二区三区免费看| 国产三级在线视频| 夜夜爽天天搞| 国产精品野战在线观看| 丁香欧美五月| 久久人人爽人人爽人人片va | 中文字幕免费在线视频6| 女同久久另类99精品国产91| 国产精品免费一区二区三区在线| 国产麻豆成人av免费视频| 国产探花极品一区二区| 国产成年人精品一区二区| 成熟少妇高潮喷水视频| 亚洲人成网站在线播| 久久精品人妻少妇| 中亚洲国语对白在线视频| 好看av亚洲va欧美ⅴa在| 国产成年人精品一区二区| 欧美另类亚洲清纯唯美| 久久精品人妻少妇| 国产精品久久视频播放| 国产成+人综合+亚洲专区| 国产精品日韩av在线免费观看| 久久精品国产亚洲av天美| 亚洲精品一区av在线观看| 18美女黄网站色大片免费观看| 国产伦精品一区二区三区四那| 日本三级黄在线观看| 国产一区二区三区视频了| 日本a在线网址| 窝窝影院91人妻| 国产伦精品一区二区三区视频9| 美女 人体艺术 gogo| 日本成人三级电影网站| 日韩欧美免费精品| 亚洲av成人精品一区久久| 午夜a级毛片| 免费在线观看亚洲国产| 午夜两性在线视频| 成人午夜高清在线视频| 日日摸夜夜添夜夜添小说| 国产主播在线观看一区二区| 日本黄色片子视频| 亚洲专区中文字幕在线| 日韩高清综合在线| 亚洲成av人片在线播放无| 精品久久久久久久久av| 直男gayav资源| 国产一级毛片七仙女欲春2| 亚洲精华国产精华精| 一个人看视频在线观看www免费| 99在线人妻在线中文字幕| 午夜久久久久精精品| av在线观看视频网站免费| 麻豆国产av国片精品| 激情在线观看视频在线高清| 18禁在线播放成人免费| 日韩欧美国产在线观看| 亚洲成人中文字幕在线播放| 热99re8久久精品国产| 成人特级黄色片久久久久久久| 观看美女的网站| 可以在线观看毛片的网站| 国产精品,欧美在线| 国产成人欧美在线观看| 一区二区三区高清视频在线| 精品久久久久久久末码| 少妇熟女aⅴ在线视频| 俄罗斯特黄特色一大片| 99久久精品一区二区三区| 91麻豆精品激情在线观看国产| 国产v大片淫在线免费观看| 波多野结衣高清无吗| 日韩亚洲欧美综合| 亚洲成人免费电影在线观看| 69人妻影院| 欧美+日韩+精品| 哪里可以看免费的av片| 久久人妻av系列| 男女做爰动态图高潮gif福利片| 国产精品久久久久久久久免 | 男女下面进入的视频免费午夜| 免费在线观看亚洲国产| 12—13女人毛片做爰片一| 免费观看的影片在线观看| 亚洲经典国产精华液单 | 国产精品一区二区免费欧美| 日韩欧美一区二区三区在线观看| 老司机福利观看| 亚洲精品日韩av片在线观看| 国产视频内射| 人人妻人人澡欧美一区二区| 热99re8久久精品国产| 国产高清视频在线播放一区| 久久性视频一级片| 99视频精品全部免费 在线| 免费人成视频x8x8入口观看| 国产精品99久久久久久久久| 日韩欧美免费精品| 亚洲无线观看免费| 极品教师在线免费播放| 赤兔流量卡办理| 国产黄a三级三级三级人| 久久久精品大字幕| 亚洲,欧美精品.| 精品久久久久久久久久久久久| 亚洲自拍偷在线| 色视频www国产| 村上凉子中文字幕在线| 国产真实乱freesex| 一进一出抽搐gif免费好疼| 国产人妻一区二区三区在| 嫩草影视91久久| 国产麻豆成人av免费视频| 亚洲专区中文字幕在线| 精品熟女少妇八av免费久了| 午夜福利高清视频| 精品一区二区三区视频在线| 不卡一级毛片| 亚洲五月婷婷丁香| 午夜两性在线视频| 最后的刺客免费高清国语| 在线观看66精品国产| 日本 欧美在线| 99久久九九国产精品国产免费| 狠狠狠狠99中文字幕| 好看av亚洲va欧美ⅴa在| 18禁在线播放成人免费| 91狼人影院| 日本a在线网址| 亚洲av五月六月丁香网| 欧美黑人巨大hd| 国产精品亚洲一级av第二区| 欧美成狂野欧美在线观看| 色哟哟哟哟哟哟| 真人一进一出gif抽搐免费| 国内精品美女久久久久久| 午夜精品在线福利| 1000部很黄的大片| 精品人妻熟女av久视频| 亚洲精品在线观看二区| 国产单亲对白刺激| 男女下面进入的视频免费午夜| 亚洲不卡免费看| 免费电影在线观看免费观看| 久久久久久久亚洲中文字幕 | 久久伊人香网站| 一进一出好大好爽视频| 国产高清三级在线| 国产成年人精品一区二区| 成人欧美大片| 九九热线精品视视频播放| 韩国av一区二区三区四区| 超碰av人人做人人爽久久| 欧美激情久久久久久爽电影| 亚洲国产精品成人综合色| 亚洲经典国产精华液单 | 三级男女做爰猛烈吃奶摸视频| 99riav亚洲国产免费| 精品国产亚洲在线| 别揉我奶头~嗯~啊~动态视频| 97碰自拍视频| 欧美中文日本在线观看视频| 一a级毛片在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲18禁久久av| 赤兔流量卡办理| 精品免费久久久久久久清纯| 免费无遮挡裸体视频| 桃红色精品国产亚洲av| 国产69精品久久久久777片| 亚洲真实伦在线观看| 免费看日本二区| 久久亚洲精品不卡| 欧美国产日韩亚洲一区| 黄片小视频在线播放| 免费人成在线观看视频色| 亚洲狠狠婷婷综合久久图片| 久久精品国产亚洲av涩爱 | 国产真实乱freesex| 99国产极品粉嫩在线观看| 国产精品亚洲av一区麻豆| 久久精品国产亚洲av香蕉五月| 老女人水多毛片| 变态另类丝袜制服| 欧美日本视频| av视频在线观看入口| 又紧又爽又黄一区二区| 久久久色成人| 中出人妻视频一区二区| 中文字幕久久专区| 噜噜噜噜噜久久久久久91| 特级一级黄色大片|