• <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ù)的重要不等式
    狠狠婷婷综合久久久久久88av| 亚洲成人免费电影在线观看| 精品第一国产精品| 婷婷丁香在线五月| 亚洲国产av新网站| 视频区图区小说| 亚洲国产av新网站| 老司机午夜十八禁免费视频| 91国产中文字幕| 国产黄频视频在线观看| 一区二区三区精品91| 欧美日韩成人在线一区二区| 中亚洲国语对白在线视频| 欧美日韩亚洲综合一区二区三区_| 欧美激情高清一区二区三区| 老司机影院毛片| 欧美日本中文国产一区发布| 精品人妻熟女毛片av久久网站| 久久亚洲精品不卡| 久久国产精品男人的天堂亚洲| 人人妻,人人澡人人爽秒播| 中文字幕制服av| 999精品在线视频| 亚洲中文字幕日韩| 手机成人av网站| 欧美激情久久久久久爽电影 | 又大又爽又粗| 久久人妻福利社区极品人妻图片| 中文字幕另类日韩欧美亚洲嫩草| 国产成人精品在线电影| 日韩视频一区二区在线观看| 国产av精品麻豆| 国产欧美亚洲国产| 免费在线观看完整版高清| 19禁男女啪啪无遮挡网站| 一二三四社区在线视频社区8| 黄色视频不卡| 视频区欧美日本亚洲| 99精国产麻豆久久婷婷| 性高湖久久久久久久久免费观看| 国产精品.久久久| 肉色欧美久久久久久久蜜桃| 永久免费av网站大全| 久久女婷五月综合色啪小说| 国产成人精品在线电影| 久久精品人人爽人人爽视色| 国产精品一区二区在线观看99| 91成年电影在线观看| 日韩欧美免费精品| 久久精品国产亚洲av高清一级| 亚洲精华国产精华精| a在线观看视频网站| 久久久国产成人免费| 国产熟女午夜一区二区三区| 久9热在线精品视频| 成人免费观看视频高清| 美女高潮喷水抽搐中文字幕| 久久久久网色| 老司机深夜福利视频在线观看 | 俄罗斯特黄特色一大片| 免费在线观看完整版高清| 国产在线免费精品| av有码第一页| 精品国产一区二区三区久久久樱花| 99国产综合亚洲精品| 日韩欧美免费精品| 麻豆av在线久日| 亚洲精品久久久久久婷婷小说| 在线观看免费日韩欧美大片| 男人操女人黄网站| 视频区图区小说| 亚洲成国产人片在线观看| 777米奇影视久久| 亚洲国产中文字幕在线视频| 中国国产av一级| 一区二区三区四区激情视频| 王馨瑶露胸无遮挡在线观看| 在线观看舔阴道视频| 免费日韩欧美在线观看| 狂野欧美激情性bbbbbb| 一边摸一边做爽爽视频免费| 飞空精品影院首页| 无遮挡黄片免费观看| 国产免费av片在线观看野外av| 亚洲久久久国产精品| 成在线人永久免费视频| 亚洲一码二码三码区别大吗| 国产成人影院久久av| 国产欧美日韩综合在线一区二区| 极品人妻少妇av视频| 狂野欧美激情性xxxx| 中国国产av一级| 一区二区三区四区激情视频| 日韩免费高清中文字幕av| 亚洲国产精品成人久久小说| 亚洲精品在线美女| 91字幕亚洲| 久久精品国产a三级三级三级| 麻豆乱淫一区二区| 99热全是精品| 好男人电影高清在线观看| 欧美精品高潮呻吟av久久| 午夜免费观看性视频| 国产又爽黄色视频| 日韩欧美一区视频在线观看| 精品乱码久久久久久99久播| 成人亚洲精品一区在线观看| av网站免费在线观看视频| 一区二区三区精品91| 亚洲国产av影院在线观看| 色综合欧美亚洲国产小说| 精品一区在线观看国产| 人人澡人人妻人| 亚洲av成人一区二区三| 下体分泌物呈黄色| 亚洲专区中文字幕在线| 巨乳人妻的诱惑在线观看| 欧美在线黄色| 久久av网站| 亚洲性夜色夜夜综合| 满18在线观看网站| 国产精品偷伦视频观看了| 精品久久久久久电影网| 亚洲欧洲日产国产| 青青草视频在线视频观看| 精品国产国语对白av| 亚洲午夜精品一区,二区,三区| 午夜两性在线视频| 韩国精品一区二区三区| 国产精品麻豆人妻色哟哟久久| av视频免费观看在线观看| 岛国毛片在线播放| 国产成人精品无人区| 日本av免费视频播放| av片东京热男人的天堂| 欧美中文综合在线视频| 久久久久久久精品精品| 一本久久精品| 婷婷成人精品国产| av视频免费观看在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 精品视频人人做人人爽| 91麻豆精品激情在线观看国产 | 永久免费av网站大全| 别揉我奶头~嗯~啊~动态视频 | 丁香六月天网| 国产精品久久久人人做人人爽| 女性生殖器流出的白浆| 亚洲欧美一区二区三区黑人| 国产精品成人在线| 国产成人欧美在线观看 | 亚洲国产日韩一区二区| 久久99一区二区三区| 欧美精品一区二区免费开放| 久久久久国产一级毛片高清牌| 国产高清国产精品国产三级| 亚洲色图综合在线观看| 色婷婷久久久亚洲欧美| 婷婷色av中文字幕| 免费黄频网站在线观看国产| 99国产精品免费福利视频| 男女之事视频高清在线观看| 在线观看免费日韩欧美大片| 亚洲国产欧美在线一区| 看免费av毛片| 欧美精品亚洲一区二区| 91字幕亚洲| 免费在线观看完整版高清| 国产成人精品无人区| 日本一区二区免费在线视频| 久久久久视频综合| 久久青草综合色| 老熟妇乱子伦视频在线观看 | 欧美日本中文国产一区发布| 交换朋友夫妻互换小说| 国产一区二区三区av在线| 亚洲精品乱久久久久久| 精品国内亚洲2022精品成人 | 免费在线观看完整版高清| 久久精品国产亚洲av高清一级| 69精品国产乱码久久久| 亚洲精品国产一区二区精华液| 国产免费视频播放在线视频| 午夜激情av网站| 国产深夜福利视频在线观看| 久久精品aⅴ一区二区三区四区| www.av在线官网国产| 中亚洲国语对白在线视频| 国产麻豆69| 久久久久久久久免费视频了| 色老头精品视频在线观看| 亚洲全国av大片| 国产一卡二卡三卡精品| 欧美日韩亚洲综合一区二区三区_| 一区在线观看完整版| 亚洲欧洲日产国产| 免费日韩欧美在线观看| 久久人人97超碰香蕉20202| 免费在线观看完整版高清| 男女下面插进去视频免费观看| 国产成人欧美在线观看 | 青春草视频在线免费观看| 日韩三级视频一区二区三区| 男女边摸边吃奶| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利视频精品| 国产视频一区二区在线看| 欧美精品高潮呻吟av久久| 捣出白浆h1v1| 永久免费av网站大全| 人人妻人人澡人人看| 人人妻人人添人人爽欧美一区卜| 老熟妇乱子伦视频在线观看 | 在线天堂中文资源库| av网站免费在线观看视频| 精品一区二区三区四区五区乱码| 91大片在线观看| 亚洲欧美一区二区三区久久| 999精品在线视频| 国产精品国产av在线观看| 国产伦人伦偷精品视频| 欧美人与性动交α欧美精品济南到| 夫妻午夜视频| 又黄又粗又硬又大视频| 亚洲av成人不卡在线观看播放网 | 日本猛色少妇xxxxx猛交久久| 亚洲av成人一区二区三| 国产精品久久久久久精品古装| 午夜福利一区二区在线看| 日韩欧美免费精品| 亚洲国产av新网站| 国精品久久久久久国模美| 国产精品免费视频内射| 午夜免费成人在线视频| 十八禁网站免费在线| 黄片大片在线免费观看| www.熟女人妻精品国产| 老司机深夜福利视频在线观看 | 国产精品成人在线| 美女扒开内裤让男人捅视频| 超碰成人久久| 涩涩av久久男人的天堂| 国产三级黄色录像| 91国产中文字幕| 久久影院123| 欧美日韩av久久| 在线观看一区二区三区激情| 深夜精品福利| 中国美女看黄片| 黑人猛操日本美女一级片| 不卡av一区二区三区| 久久久久久久国产电影| 91精品三级在线观看| 亚洲美女黄色视频免费看| 国产精品 欧美亚洲| www日本在线高清视频| 国产成人系列免费观看| 亚洲性夜色夜夜综合| 国产麻豆69| 9191精品国产免费久久| av线在线观看网站| 日韩,欧美,国产一区二区三区| 超色免费av| 亚洲欧洲精品一区二区精品久久久| 超碰成人久久| av在线app专区| 美女高潮到喷水免费观看| 欧美成狂野欧美在线观看| 久久人人爽av亚洲精品天堂| 啦啦啦啦在线视频资源| 久久狼人影院| 亚洲一码二码三码区别大吗| 2018国产大陆天天弄谢| 老司机亚洲免费影院| 涩涩av久久男人的天堂| 国产极品粉嫩免费观看在线| 人妻久久中文字幕网| 一本久久精品| 精品亚洲成国产av| 99精国产麻豆久久婷婷| 成人国产av品久久久| 最新在线观看一区二区三区| 男男h啪啪无遮挡| 美女视频免费永久观看网站| 午夜福利一区二区在线看| 国产精品久久久久久人妻精品电影 | 久久久国产精品麻豆| 国产成人一区二区三区免费视频网站| 久久国产亚洲av麻豆专区| 亚洲欧美精品综合一区二区三区| 日本欧美视频一区| 香蕉丝袜av| 日韩大片免费观看网站| 人人妻人人添人人爽欧美一区卜| 97精品久久久久久久久久精品| 男女无遮挡免费网站观看| 在线看a的网站| 爱豆传媒免费全集在线观看| 丁香六月欧美| 热re99久久精品国产66热6| 亚洲 国产 在线| 美女午夜性视频免费| 男女床上黄色一级片免费看| 午夜福利在线免费观看网站| 黑丝袜美女国产一区| 日本a在线网址| 黄片大片在线免费观看| 国产一区二区三区av在线| 国产黄频视频在线观看| 热99re8久久精品国产| www.熟女人妻精品国产| 啦啦啦在线免费观看视频4| 99久久99久久久精品蜜桃| 国产老妇伦熟女老妇高清| 九色亚洲精品在线播放| 欧美日韩亚洲高清精品| 欧美在线黄色| 欧美黑人欧美精品刺激| 亚洲黑人精品在线| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看舔阴道视频| 少妇精品久久久久久久| 亚洲九九香蕉| 亚洲第一欧美日韩一区二区三区 | 99热国产这里只有精品6| 免费一级毛片在线播放高清视频 | 中文字幕色久视频| 男人爽女人下面视频在线观看| 国产成人系列免费观看| 日韩欧美免费精品| 法律面前人人平等表现在哪些方面 | 丝瓜视频免费看黄片| 免费高清在线观看日韩| e午夜精品久久久久久久| av电影中文网址| 久久久国产欧美日韩av| 高清黄色对白视频在线免费看| 国产亚洲午夜精品一区二区久久| 欧美精品亚洲一区二区| 99re6热这里在线精品视频| 久9热在线精品视频| 美女福利国产在线| 国产在线观看jvid| 中亚洲国语对白在线视频| 国产高清videossex| 又紧又爽又黄一区二区| 99精品欧美一区二区三区四区| 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| 新久久久久国产一级毛片| 亚洲欧美一区二区三区黑人| 午夜激情久久久久久久| 一级黄色大片毛片| 热99久久久久精品小说推荐| 男女边摸边吃奶| 国产精品一区二区在线不卡| 精品国产超薄肉色丝袜足j| 一二三四社区在线视频社区8| 在线观看免费午夜福利视频| 亚洲自偷自拍图片 自拍| 国产精品久久久av美女十八| 99国产精品一区二区蜜桃av | 少妇被粗大的猛进出69影院| 青草久久国产| 亚洲精品av麻豆狂野| 国产深夜福利视频在线观看| 91av网站免费观看| 久久久精品94久久精品| 丰满饥渴人妻一区二区三| 国产一区二区三区在线臀色熟女 | 一区二区av电影网| 我要看黄色一级片免费的| 国产av一区二区精品久久| 汤姆久久久久久久影院中文字幕| 欧美国产精品一级二级三级| 999久久久精品免费观看国产| 国产一区二区三区综合在线观看| 一区在线观看完整版| 99国产精品99久久久久| 亚洲国产毛片av蜜桃av| 欧美日韩亚洲综合一区二区三区_| 少妇人妻久久综合中文| 国产免费视频播放在线视频| 亚洲国产毛片av蜜桃av| 丝袜人妻中文字幕| 淫妇啪啪啪对白视频 | 蜜桃在线观看..| 99久久国产精品久久久| 亚洲av成人一区二区三| 巨乳人妻的诱惑在线观看| 久久久久国内视频| 黄色a级毛片大全视频| av福利片在线| 国产亚洲精品一区二区www | 亚洲情色 制服丝袜| 嫩草影视91久久| 老鸭窝网址在线观看| 色视频在线一区二区三区| 女人被躁到高潮嗷嗷叫费观| 亚洲精品美女久久av网站| 精品久久蜜臀av无| 亚洲五月色婷婷综合| 欧美黑人欧美精品刺激| 夜夜夜夜夜久久久久| a 毛片基地| 欧美av亚洲av综合av国产av| 丝袜人妻中文字幕| 久久人人爽av亚洲精品天堂| 精品第一国产精品| 亚洲精品一区蜜桃| 伊人亚洲综合成人网| 一级毛片精品| 99久久人妻综合| 亚洲一区二区三区欧美精品| 2018国产大陆天天弄谢| 在线十欧美十亚洲十日本专区| 亚洲国产精品999| 一区二区av电影网| 日韩一区二区三区影片| 成人亚洲精品一区在线观看| 日本猛色少妇xxxxx猛交久久| 激情视频va一区二区三区| 精品国产一区二区三区久久久樱花| 亚洲国产中文字幕在线视频| 亚洲三区欧美一区| 一二三四社区在线视频社区8| 国产激情久久老熟女| 亚洲精品在线美女| 精品国产乱码久久久久久男人| 国产精品 欧美亚洲| 乱人伦中国视频| 欧美精品人与动牲交sv欧美| 国产成人免费无遮挡视频| 美女午夜性视频免费| 成人手机av| 欧美人与性动交α欧美软件| kizo精华| 99久久精品国产亚洲精品| 亚洲精品成人av观看孕妇| 色94色欧美一区二区| 久久久精品94久久精品| 自线自在国产av| 精品福利永久在线观看| 99国产精品99久久久久| 精品视频人人做人人爽| 久久久精品国产亚洲av高清涩受| 男人爽女人下面视频在线观看| 国产精品国产av在线观看| 超碰成人久久| 精品人妻一区二区三区麻豆| 99国产精品一区二区蜜桃av | 国产精品一区二区在线观看99| 夫妻午夜视频| 新久久久久国产一级毛片| 欧美日韩视频精品一区| 色老头精品视频在线观看| 亚洲性夜色夜夜综合| 欧美久久黑人一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲人成电影免费在线| 久久国产精品影院| 精品国产乱码久久久久久男人| 精品国内亚洲2022精品成人 | 成年人黄色毛片网站| av电影中文网址| 亚洲成人国产一区在线观看| av片东京热男人的天堂| 黄网站色视频无遮挡免费观看| 欧美xxⅹ黑人| 高潮久久久久久久久久久不卡| 在线观看一区二区三区激情| 亚洲天堂av无毛| 97人妻天天添夜夜摸| 黄色视频,在线免费观看| 欧美少妇被猛烈插入视频| 涩涩av久久男人的天堂| 这个男人来自地球电影免费观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品美女久久久久99蜜臀| 中文字幕色久视频| 欧美亚洲日本最大视频资源| 国产精品 国内视频| 大片免费播放器 马上看| 一二三四社区在线视频社区8| 婷婷丁香在线五月| 久久久久久人人人人人| 纵有疾风起免费观看全集完整版| 精品少妇久久久久久888优播| 80岁老熟妇乱子伦牲交| 日韩大片免费观看网站| 欧美 日韩 精品 国产| 午夜福利免费观看在线| 国产精品一区二区免费欧美 | 亚洲国产av影院在线观看| 热re99久久精品国产66热6| 欧美在线一区亚洲| av不卡在线播放| 99精品欧美一区二区三区四区| 亚洲国产欧美在线一区| 女人被躁到高潮嗷嗷叫费观| 亚洲激情五月婷婷啪啪| 99九九在线精品视频| 国产精品九九99| 制服人妻中文乱码| av天堂在线播放| 亚洲黑人精品在线| 天天躁日日躁夜夜躁夜夜| 亚洲国产日韩一区二区| 欧美少妇被猛烈插入视频| 欧美成人午夜精品| 国产精品二区激情视频| 中文字幕最新亚洲高清| 欧美av亚洲av综合av国产av| 中文字幕人妻熟女乱码| 热99国产精品久久久久久7| 大香蕉久久网| 精品一区在线观看国产| 国产精品久久久人人做人人爽| 精品久久蜜臀av无| 在线观看人妻少妇| 国产国语露脸激情在线看| 一区二区三区激情视频| 老鸭窝网址在线观看| 亚洲欧美精品综合一区二区三区| 免费观看人在逋| 美国免费a级毛片| 国产一区二区三区在线臀色熟女 | 精品国产国语对白av| 美女高潮到喷水免费观看| 91老司机精品| h视频一区二区三区| 2018国产大陆天天弄谢| 新久久久久国产一级毛片| 黄色片一级片一级黄色片| 中文字幕制服av| 波多野结衣一区麻豆| 天堂俺去俺来也www色官网| a级毛片在线看网站| 国产在视频线精品| 国产成人精品在线电影| 91精品伊人久久大香线蕉| 99精品欧美一区二区三区四区| 中文字幕最新亚洲高清| 精品人妻在线不人妻| 美女扒开内裤让男人捅视频| tube8黄色片| 日本vs欧美在线观看视频| 久久久久久久大尺度免费视频| 国产一区二区三区av在线| 一区二区三区激情视频| 这个男人来自地球电影免费观看| 我要看黄色一级片免费的| 国产精品九九99| 国产激情久久老熟女| 一边摸一边做爽爽视频免费| 久久亚洲国产成人精品v| 日韩人妻精品一区2区三区| 色综合欧美亚洲国产小说| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久久国产电影| 亚洲人成电影观看| 国产精品99久久99久久久不卡| 国产有黄有色有爽视频| 欧美日韩国产mv在线观看视频| 19禁男女啪啪无遮挡网站| 国产精品一区二区免费欧美 | 国产精品秋霞免费鲁丝片| www.自偷自拍.com| 女警被强在线播放| 亚洲少妇的诱惑av| 两人在一起打扑克的视频| 午夜福利视频在线观看免费| 欧美黄色淫秽网站| 丝袜喷水一区| 91麻豆av在线| 午夜影院在线不卡| 91麻豆av在线| 精品视频人人做人人爽| 久久久久国产一级毛片高清牌| 国产日韩欧美亚洲二区| 日韩一卡2卡3卡4卡2021年| 少妇粗大呻吟视频| 女警被强在线播放| 在线亚洲精品国产二区图片欧美| tocl精华| 国产欧美日韩一区二区三区在线| tocl精华| 男人舔女人的私密视频| 亚洲美女黄色视频免费看| 热re99久久精品国产66热6| 大香蕉久久成人网| 不卡一级毛片| 成年人黄色毛片网站| 国产精品久久久久久精品古装| 男人操女人黄网站| 久久精品国产a三级三级三级| 日本wwww免费看| 日韩欧美国产一区二区入口| 欧美黄色淫秽网站| av视频免费观看在线观看| 91大片在线观看| 99re6热这里在线精品视频| 午夜免费鲁丝| 久久女婷五月综合色啪小说| 一本综合久久免费| 亚洲一码二码三码区别大吗| 777久久人妻少妇嫩草av网站| 欧美精品一区二区免费开放| 精品少妇一区二区三区视频日本电影| 1024视频免费在线观看| 动漫黄色视频在线观看| 亚洲国产av新网站| 亚洲成人免费av在线播放| 中文字幕色久视频| 日本欧美视频一区|