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

    長沙機場陣列天氣雷達地物識別算法

    2020-05-13 08:39:04魏萬益馬舒慶甄小瓊4呂寺煒
    應用氣象學報 2020年3期
    關鍵詞:子陣徑向速度雜波

    魏萬益 馬舒慶 楊 玲* 甄小瓊4) 呂寺煒

    1)(成都信息工程大學電子工程學院, 成都 610225)2)(中國氣象局大氣探測重點開放實驗室, 成都 610225)3)(中國氣象局氣象探測中心, 北京 100081)4)(中國科學院大氣物理研究所, 北京 100029)5)(雷象科技(北京)有限公司, 北京 100089)

    引 言

    地物雜波識別和去除是天氣雷達數(shù)據(jù)質(zhì)量控制研究中的一項基本工作。地物會影響雷達回波,掩蓋天氣回波,甚至直接影響雷達產(chǎn)品的可信度。

    為了消除地物雜波,國內(nèi)天氣雷達專家們提出了多種去除地物雜波的方法,這些方法可分為以下3類:①通過改變雷達的設計和布局位置減少地物雜波,②在信號處理端對雷達信號進行處理,③針對雷達基數(shù)據(jù)進行后期處理。Smith[1]和Mann等[2]通過選擇雷達的安放位置和采用更短的波長減少雜波及其相對強度。第2類方法通常是對收集到的所有數(shù)據(jù)采用時域[2-4]或頻域[5-7]濾波器去除速度接近0的回波數(shù)據(jù)。Li等[8]在頻譜上提取特征參數(shù)作為樸素貝葉斯算法的輸入,以判斷是否存在地物雜波并加以去除。針對雷達基數(shù)據(jù)的處理方法主要是從基數(shù)據(jù)中提取二維或三維地物識別特征參數(shù)。Steiner等[9]使用回波強度三維結(jié)構作為決策樹算法輸入?yún)?shù)構建自動化雜波抑制算法,發(fā)現(xiàn)最有效的3個參數(shù)是雷達回波的垂直范圍、回波強度的水平變化及回波強度的垂直梯度。Zhang等[10]使用類似的回波強度垂直梯度去除地物雜波,并在計算該參數(shù)時考慮仰角抬高對水平距離的影響。Lakshmanan等[11-12]采用神經(jīng)網(wǎng)絡方法,將回波強度的三維結(jié)構、速度和譜寬以及其他特征作為神經(jīng)網(wǎng)絡的輸入數(shù)據(jù)去除地物雜波。此外,Kessinger等[13]、Ferrer等[14]、Cho等[15]、劉黎平等[16]和江源等[17]在雷達基數(shù)據(jù)中提取各種特征,使用模糊邏輯算法識別并去除地物雜波。李豐等[18-19]采用類似的模糊邏輯算法識別非降水回波,并對比分析S波段和C波段多普勒天氣雷達各種特征的分布。

    Ruiz等[20]使用相控陣天氣雷達探測數(shù)據(jù),提取回波強度與時間的相關性紋理特征等參數(shù),并將其用于樸素貝葉斯算法,有效識別和去除地物雜波。但該算法運算效率較低,不能滿足快速掃描相控陣天氣雷達的產(chǎn)品生成需求。文中同時指出,與傳統(tǒng)拋物線天線天氣雷達相比,日本大阪大學的相控陣天氣雷達具有較強的旁瓣和較寬的波束,會導致雷達數(shù)據(jù)中存在更多的地物雜波,對地物雜波的識別和去除算法要求更高。

    為了進一步研究小尺度強對流天氣系統(tǒng),在國內(nèi)外網(wǎng)絡化和相控陣技術發(fā)展的背景下[21-23],2015年中國氣象局氣象探測中心設計并聯(lián)合有關廠家開始研制陣列天氣雷達,它結(jié)合網(wǎng)絡化雷達和分布式相控陣技術,具有高時空分辨率的特點。本文采用基于模糊邏輯的地物識別算法,參考日本大阪大學的相控陣天氣雷達地物雜波去除算法,提出了一個新的參數(shù)——回波強度時間變化量(time variability of reflectivity factor,TVR),并在不同天氣情況下驗證了地物雜波去除算法的有效性。

    1 陣列天氣雷達簡介

    陣列天氣雷達是一種分布式相控陣天氣雷達,由控制處理中心和至少3個相控陣接收發(fā)射子陣(簡稱收發(fā)子陣)構成,通過增加收發(fā)子陣數(shù)量擴大探測區(qū)域[24]。每3個相鄰的收發(fā)子陣為一組進行協(xié)同掃描,保證3個相鄰收發(fā)子陣在空間同一點的數(shù)據(jù)時差小于2 s。

    收發(fā)子陣是陣列天氣雷達的前端,由于缺少獨立的控制和數(shù)據(jù)處理部分,不是完整的雷達,故被稱為收發(fā)子陣而非相控陣雷達。收發(fā)子陣由相控陣天線陣列、收發(fā)模塊陣列、信號處理陣列、方位旋轉(zhuǎn)伺服、同步和通信模塊等部分組成。水平方向采用機械掃描方式覆蓋0°~360°,垂直方向發(fā)射4個寬度為22.5°的波束覆蓋0°~90°,接收到回波信號后,通過數(shù)字波束得到16個平均波束寬度為1.6°的波束。2018年3月在長沙機場布設的陣列天氣雷達,天線旋轉(zhuǎn)速度為30 °·s-1,完成1個體掃僅需要12 s,約為現(xiàn)在業(yè)務中使用的多普勒天氣雷達常用體掃模式1個體掃時間的1/30。每個收發(fā)子陣的最大探測距離為20.28 km,徑向分辨率為30 m。

    2 基于模糊邏輯的地物識別算法

    2.1 數(shù)據(jù)預處理

    由于雷達基數(shù)據(jù)中含有部分孤立的點狀或線狀回波(圖1a),本文采用Zhang等[10]提出的孤立點濾波算法對雷達基數(shù)據(jù)進行預處理。具體算法如下:以空間中點X為中心位置,在方位、徑向和俯仰上5×5×1個距離庫中有效回波點的百分比為PX,如果PX小于預設閾值(默認為75%),則認為該點為孤立點并將其剔除。圖1為孤立點濾波前后的回波對比。由圖1b可知,通過孤立點濾波可以消除大部分孤立回波點。

    圖1 孤立點濾波前(a)、濾波后(b)回波強度Fig.1 Intensity echo of isolated point before filtering(a) and after filtering(b)

    2.2 識別算法介紹

    陣列天氣雷達的時間分辨率空間分辨率均高于常規(guī)多普勒天氣雷達。本文在已有地物特征參數(shù)基礎上加入表征回波強度隨時間變化的參數(shù),回波強度時間變化量(TVR),作為模糊邏輯算法的輸入?yún)?shù)來識別地物雜波[25-27]。其他特征參數(shù)包括回波強度紋理(TDBZ)、回波強度垂直梯度(GDBZ)、徑向庫間回波強度變化程度(SPIN)、徑向速度平均值(MDVE)和速度譜寬平均值(MDSW)。這些參數(shù)的定義見式(1)~(7)[16]:

    (1)

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    其中,NA,NR表示在方位和徑向定義的計算范圍;i,j為范圍內(nèi)的索引;Z,V和SW分別為回波強度、徑向速度和速度譜寬;Zlow和Zup為對應的本層和上層PPI的回波強度;Zthresh為徑向上相鄰距離庫強度變化的閾值;ZN和ZL分別表示當前時刻和上一時刻的強度。在實際特征值計算過程中,設定NA=5,NR=5,Zthresh=2 dBZ,ZN和ZL的時間差為1 min。

    除少數(shù)特殊地物外,大部分地物位置固定不動且形狀大小不隨時間改變,所以地物雜波的徑向速度主要分布在0附近,譜寬相對較小,回波強度時間變化量(TVR)接近0。相比降水回波,地物雜波大部分出現(xiàn)在低仰角數(shù)據(jù)中,強度在空間上變化大,且隨著仰角的抬高,強度迅速減小。因此地物雜波的強度紋理(TDBZ)和徑向庫間強度變化程度(SPIN)較大,強度的垂直梯度(GDBZ)為負值且絕對值比較大。

    本文收集了2019年3—7月不同天氣情況下長沙機場陣列天氣雷達基數(shù)據(jù),根據(jù)地物雜波和天氣回波的特征預先確定回波類型。判斷依據(jù):地物主要出現(xiàn)在較低的仰角,且隨著仰角的抬高強度迅速減小,地物的徑向速度接近零,位置不隨時間改變;對于降水回波,可根據(jù)回波形狀、垂直結(jié)構、演變和移動等方面確定[16]。使用部分回波數(shù)據(jù),統(tǒng)計得到地物和降水各個參數(shù)的概率分布情況(圖2),其余回波數(shù)據(jù)驗證,檢驗算法的識別效果。

    由圖2可知,地物的徑向速度平均值(MDVE)集中分布在0附近,與降水的區(qū)別較大(圖2a)。地物和降水回波的速度譜寬平均值(MDSW)重合面積較大,區(qū)別較小,未達到預期效果,因此,在識別算法中剔除該參數(shù)(圖2b)。與降水回波相比,地物的回波強度紋理(TDBZ)分布更加離散,符合地物雜波空間分布差異較大的特征(圖2c)。地物雜波的強度垂直梯度(GDBZ)存在少數(shù)大于0的部分,與預期的效果存在差異,但與降水回波明顯不同,不影響作為識別地物的輸入?yún)?shù)(圖2d)。地物和降水回波的強度徑向庫間變化程度(SPIN)分布的波峰位置區(qū)別明顯,地物雜波的波峰位置更靠右,即在計算區(qū)域內(nèi),地物雜波強度波動超過2 dB的數(shù)據(jù)點更多(圖2e)?;夭◤姸葧r間變化量(TVR)在地物和降水回波中有明顯差異,地物雜波主要分布在0~2 dB范圍內(nèi),而降水回波的分布則相對離散(圖2f)。綜上所述,地物和降水回波的大部分參數(shù)存在較明顯差異,但也有部分重合區(qū)域,如在回波強度紋理(TDBZ)靠近0的區(qū)域也有部分地物雜波,徑向速度平均值(MDVE)較小的區(qū)域也有降水回波。因此,使用單一特征參數(shù)難以識別地物雜波,需綜合使用多種特征參數(shù)才能較好地區(qū)分地物和降水回波。

    為了簡化計算,本文根據(jù)地物和降水特征分布的范圍、交點和波峰位置確定采用梯形隸屬函數(shù)(圖3)。根據(jù)隸屬函數(shù),計算每個點各參數(shù)的判據(jù)M。然后,對每個點的所有判據(jù)(N個)進行累加平均(式(8)),可得到每個點地物的判斷值EM(0~1)。該值越大,代表該點是地物的可能性越大,反之,代表降水回波的可能性越大。最后將得到的判斷值EM與預先設定的閾值進行比較,大于閾值判斷為地物,反之判斷為降水。

    (8)

    圖2 地物和降水回波各特征參數(shù)的概率分布(a)MDVE,(b)MDSW,(c)TDBZ,(d)GDBZ,(e)SPIN,(f)TVRFig.2 Probability distributions of ground clutter and precipitation echo(a)MDVE,(b)MDSW,(c)TDBZ,(d)GDBZ,(e)SPIN,(f)TVR

    圖3 地物識別的隸屬函數(shù)(a)MDVE,(b)TDBZ,(c)GDBZ,(d)SPIN,(e)TVRFig.3 Membership functions of ground clutter detection(a)MDVE,(b)TDBZ,(c)GDBZ,(d)SPIN,(e)TVR

    3 地物識別效果

    不同閾值的回波強度時間變化量(TVR)參數(shù)對識別效果影響見表1。隨著閾值的提高,地物識別準確率和降水誤判率降低,在閾值不大于0.55時,TVR對改善地物識別準確率和降水識別誤判率有貢獻。當閾值大于0.55后,TVR不再對識別算法產(chǎn)生影響。

    雖然閾值為0.4時地物識別的準確率最高,但降水識別的誤判率也最大。綜合考慮地物識別的準確率和降水識別誤判率,最終選取閾值為0.45,并用于分析3個子陣在無降水、混合性降水和對流性降水天氣情況的地物識別算法效果。

    表1 地物識別準確率和降水識別誤判率Table 1 Accuracy of ground clutter detection and error rate of precipitation detection

    3.1 無降水情況地物識別效果

    為了檢驗算法在無降水情況下的地物識別效果,采用2019年3—7月無降水的雷達數(shù)據(jù),分析了3個子陣在無降水情況下算法識別效果(表2)。結(jié)果表明:地物識別算法在無降水情況下表現(xiàn)良好,地物識別的準確率高達96%。圖4為子陣1在2019年7月31日10:21地物識別前后的對比,經(jīng)過地物識別算法,大部分地物雜波可被有效識別并去除。

    表2 無降水情況下地物識別算法效果Table 2 Accuracy of ground clutter detection algorithm under no precipitation condition

    圖4 2019年7月31日10:21子陣1在1.4°仰角的回波強度(相鄰距離圈間距為5 km)(a)地物識別前,(b)地物識別后Fig.4 The echo intensity of subarray 1 at 1.4° elevation angle at 1027 BT 31 Jul 2019(the distance between adjacent rang rings is 5 km)(a)before ground clutter detection,(b)after ground clutter detection

    3.2 混合性降水情況地物識別效果

    分析混合性降水情況下地物識別算法效果(表3)。在此類降水情況下,各個子陣地物雜波識別的準確率均達到91%以上,而降水回波識別的誤判率最高僅10%。圖5為子陣1在2019年6月21日15:17地物識別前后的對比。由圖5可知,降水識別誤判率較高主要由降水回波邊緣部分的誤判引起。相對于無降水情況,地物識別準確率略有下降,從速度對比可以看到,與降水回波重合的地物雜波部分未被準確識別。

    表3 混合性降水算法識別效果Table 3 Accuracy of ground clutter detection and error rate of precipitation detection under mixed prcipitation condition

    圖5 2019年6月21日15:17子陣1在1.4°仰角的回波強度和徑向速度(相鄰距離圈間距為5 km)(a)識別前的回波強度,(b)識別后的回波強度,(c)識別前的徑向速度,(d)識別后的徑向速度Fig.5 The echo intensity and radial velocity of subarray 1 at 1.4°elevation angle at 1517 BT 21 Jun 2019(the distance between adjacent range rings is 5 km)(a)echo intensity before ground clutter detection,(b)echo intensity after ground clutter detection,(c)radial velocity before ground clutter detection,(d)radial velocity after ground clutter detection

    3.3 對流性降水情況地物識別效果

    分析對流性降水情況下地物識別算法效果見表4。由于子陣1和子陣3的個例中地物雜波和降水回波無重合部分,地物識別的準確率與無降水情況下差別不大,分別達到92%和94%。子陣2的大部分個例中降水和地物雜波重合區(qū)域中存在不能被準確識別部分,準確率僅91%。圖6和圖7分別是2019年7月21日13:57子陣2的1.4°和2.8°仰角地物識別前后的對比,2.8°仰角地物明顯減少,降水區(qū)域被更完整保留。因此,該類型降水情況下,算法識別處理后,大部分降水回波被較好保留,僅降水回波邊緣部分和西南方向一小塊孤立的降水回波被錯誤識別。

    表4 對流性降水算法識別效果Table 4 Accuracy of ground clutter detection and error rate of precipitation detection under convective precipitation condition

    圖6 2019年7月21日13:57子陣2在1.4°仰角的回波強度和徑向速度(相鄰距離圈間距為5 km)(a)識別前的回波強度,(b)識別后的回波強度,(c)識別前的徑向速度,(d)識別后的徑向速度Fig.6 The echo intensity and radial velocity of subarray 2 at 1.4°elevation angle at 1357 BT 21 Jul 2019(the distance between adjacent range rings is 5 km)(a)echo intensity before ground clutter detection,(b)echo intensity after ground clutter detection,(c)radial velocity before ground clutter detection,(d)radial velocity after ground clutter detection

    圖7 2019年7月21日13:57子陣2在2.8°仰角的回波強度和徑向速度(相鄰距離圈間距為5 km)(a)識別前的回波強度,(b)識別后的回波強度,(c)識別前的徑向速度,(d)識別后的徑向速度Fig.7 The echo intensity and radial velocity of subarray 2 at 2.8° elevation angle at 1357 BT 21 Jul 2019(the distance between adjacent range rings is 5 km)(a)echo intensity before ground clutter detection,(b)echo intensity after ground clutter detection,(c)radial velocity before ground clutter detection,(d)radial velocity after ground clutter detection

    續(xù)圖7

    4 小 結(jié)

    本文利用陣列天氣雷達高時空分辨率的特點,提出回波強度時間變化量(TVR)這一新參數(shù)。采用2019年3—7月長沙黃花機場X波段陣列天氣雷達數(shù)據(jù),分析地物識別算法中6個特征參數(shù)的統(tǒng)計特征,改進基于模糊邏輯的地物識別算法,并利用雷達基數(shù)據(jù)分析TVR對算法的貢獻和不同天氣情況下識別算法的效果,得到以下結(jié)論:

    1) 降水回波和地物雜波回波強度紋理(TDBZ)、回波強度垂直梯度(GDBZ)、徑向庫間回波強度變化程度(SPIN)、回波強度時間變化量(TVR)和徑向速度平均值(MDVE)5個參數(shù)特征分布差異較大,但速度譜寬平均值(MDSW)差異較小,因此,只采用前5個參數(shù)作為模糊邏輯算法的輸入特征。

    2) 對回波強度時間變化量(TVR)的貢獻程度進行定量評價,通過選擇最優(yōu)的閾值,采用該參數(shù)可將地物識別準確率提高4%,降水識別誤判率降低2%。

    3) 在不同天氣情況下分別對地物識別算法的表現(xiàn)進行驗證。無降水時,地物識別的準確率為94%;有降水時,地物識別準確率為92%,降水識別誤判率穩(wěn)定在10%左右。

    本算法在陣列天氣雷達地物識別中取得了較好的效果,為生成更加準確的氣象產(chǎn)品提供了保障。但還存在一些問題:降水回波的邊緣部分可能被誤判,且當降水回波與地物回波混合在一起時,去除地物回波也損失了降水回波信息。在后續(xù)的研究中,將考慮采用多子陣聯(lián)合的地物去除方法,進一步提高識別準確率,以及補償與地物一起去除的降水回波。

    猜你喜歡
    子陣徑向速度雜波
    低副瓣AiP 混合子陣稀布陣設計
    STAR2000型空管一次雷達雜波抑制淺析
    子陣劃分對相控陣設備性能影響
    非圓形光纖研究進展
    航空兵器(2017年6期)2018-01-24 15:00:10
    臺風威馬遜造成云南文山州強降水天氣雷達回波分析
    密集雜波環(huán)境下確定性退火DA-HPMHT跟蹤算法
    相關廣義復合分布雷達海雜波仿真
    遙測遙控(2015年2期)2015-04-23 08:15:22
    距離頻率ML方法無模糊估計動目標徑向速度
    微波雷達海雜波混沌特性分析
    一種平面陣的非均勻子陣劃分方法
    国产视频一区二区在线看| 国产成人欧美在线观看| 欧美又色又爽又黄视频| 国产高清视频在线播放一区| 黄片播放在线免费| 国产亚洲欧美98| 亚洲精品粉嫩美女一区| 国产精品九九99| 国产高清激情床上av| 国产一区二区三区视频了| 亚洲国产日韩欧美精品在线观看 | 老熟妇乱子伦视频在线观看| 男女做爰动态图高潮gif福利片| 国产亚洲精品久久久久5区| 午夜老司机福利片| 香蕉久久夜色| 国产不卡一卡二| 亚洲九九香蕉| 亚洲免费av在线视频| 窝窝影院91人妻| 亚洲精品一卡2卡三卡4卡5卡| 精品国产一区二区三区四区第35| 亚洲精品中文字幕一二三四区| 亚洲欧美激情综合另类| 一本大道久久a久久精品| 婷婷精品国产亚洲av在线| 国产午夜精品久久久久久| 亚洲成av人片免费观看| 日韩大尺度精品在线看网址| 精品电影一区二区在线| 国产黄a三级三级三级人| 香蕉国产在线看| 欧美成人免费av一区二区三区| 成人18禁在线播放| 一级毛片女人18水好多| 国产av一区在线观看免费| 午夜老司机福利片| 真人一进一出gif抽搐免费| 成人手机av| 久久精品国产亚洲av高清一级| 在线永久观看黄色视频| 波多野结衣av一区二区av| 亚洲真实伦在线观看| 亚洲第一av免费看| 午夜a级毛片| 搡老岳熟女国产| 12—13女人毛片做爰片一| 国产精品一区二区精品视频观看| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久久久久久久久 | 久久久精品国产亚洲av高清涩受| 一本一本综合久久| 三级毛片av免费| 亚洲成a人片在线一区二区| 这个男人来自地球电影免费观看| 亚洲精华国产精华精| 亚洲国产欧美一区二区综合| 免费人成视频x8x8入口观看| 久久亚洲真实| 久久国产精品人妻蜜桃| 午夜视频精品福利| 叶爱在线成人免费视频播放| 精品久久久久久成人av| 成人特级黄色片久久久久久久| 国产真人三级小视频在线观看| 亚洲一区二区三区色噜噜| 999精品在线视频| 亚洲第一av免费看| 日韩高清综合在线| 久99久视频精品免费| 欧美激情 高清一区二区三区| 色精品久久人妻99蜜桃| 欧美大码av| 亚洲成人免费电影在线观看| 久久久久久久午夜电影| 中出人妻视频一区二区| 黑人巨大精品欧美一区二区mp4| 他把我摸到了高潮在线观看| 国产精华一区二区三区| 久久久久久久久中文| 久久欧美精品欧美久久欧美| 成人午夜高清在线视频 | 日韩成人在线观看一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 777久久人妻少妇嫩草av网站| 亚洲国产看品久久| 两性夫妻黄色片| 高潮久久久久久久久久久不卡| 亚洲精品在线美女| 国产一区二区三区在线臀色熟女| 欧美中文日本在线观看视频| 国产精品一区二区精品视频观看| 在线天堂中文资源库| 一进一出抽搐动态| 757午夜福利合集在线观看| 亚洲片人在线观看| 1024手机看黄色片| 女人爽到高潮嗷嗷叫在线视频| 日韩成人在线观看一区二区三区| av视频在线观看入口| 欧美一级毛片孕妇| 久久精品91蜜桃| 波多野结衣高清作品| 午夜精品久久久久久毛片777| 亚洲精品av麻豆狂野| 亚洲精品在线观看二区| 久久天躁狠狠躁夜夜2o2o| 一本精品99久久精品77| 亚洲av成人不卡在线观看播放网| 欧美成人性av电影在线观看| 国产精品亚洲av一区麻豆| 亚洲va日本ⅴa欧美va伊人久久| 美女午夜性视频免费| 一级片免费观看大全| 最新美女视频免费是黄的| 国产一级毛片七仙女欲春2 | 亚洲专区国产一区二区| 女生性感内裤真人,穿戴方法视频| 国产成人一区二区三区免费视频网站| 亚洲av五月六月丁香网| 亚洲精品美女久久av网站| 国产极品粉嫩免费观看在线| 国产精品亚洲美女久久久| 国产一区在线观看成人免费| tocl精华| 国产成人av教育| 精品国产国语对白av| 久久人人精品亚洲av| 免费人成视频x8x8入口观看| 国产精品野战在线观看| 国产三级在线视频| 国产色视频综合| 国产1区2区3区精品| 亚洲自拍偷在线| 日韩av在线大香蕉| 99精品在免费线老司机午夜| 波多野结衣av一区二区av| 母亲3免费完整高清在线观看| 久久午夜综合久久蜜桃| 男女午夜视频在线观看| 亚洲精华国产精华精| 怎么达到女性高潮| 精品人妻1区二区| av视频在线观看入口| 91大片在线观看| 亚洲全国av大片| 一夜夜www| 中文字幕最新亚洲高清| 一区二区三区高清视频在线| 黄色视频,在线免费观看| 老司机福利观看| 国产精品亚洲一级av第二区| 美女高潮到喷水免费观看| 夜夜夜夜夜久久久久| 欧美日韩精品网址| 亚洲国产欧洲综合997久久, | 中文字幕人成人乱码亚洲影| 麻豆国产av国片精品| 在线观看舔阴道视频| 亚洲精品久久成人aⅴ小说| 两性夫妻黄色片| 黑人操中国人逼视频| 午夜福利欧美成人| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品乱码一区二三区的特点| 亚洲精品美女久久久久99蜜臀| 国产视频内射| 欧美成人午夜精品| 日韩欧美免费精品| 啦啦啦观看免费观看视频高清| 男女午夜视频在线观看| 后天国语完整版免费观看| av电影中文网址| 日韩精品免费视频一区二区三区| av福利片在线| 黑人巨大精品欧美一区二区mp4| av有码第一页| 亚洲国产精品久久男人天堂| 天天躁夜夜躁狠狠躁躁| 大型av网站在线播放| 淫妇啪啪啪对白视频| 天天躁夜夜躁狠狠躁躁| 欧美zozozo另类| 亚洲人成网站高清观看| 午夜福利在线在线| 91国产中文字幕| 欧美日本视频| 99热只有精品国产| 亚洲免费av在线视频| 一级a爱视频在线免费观看| 哪里可以看免费的av片| 十分钟在线观看高清视频www| 美女国产高潮福利片在线看| 国产午夜精品久久久久久| 欧美成人免费av一区二区三区| 欧美日韩黄片免| 日本免费a在线| 中文在线观看免费www的网站 | 国产三级黄色录像| 成人18禁在线播放| 欧美又色又爽又黄视频| 久久九九热精品免费| 91在线观看av| 他把我摸到了高潮在线观看| 国产亚洲欧美98| 中文亚洲av片在线观看爽| 听说在线观看完整版免费高清| 日本在线视频免费播放| 久久国产乱子伦精品免费另类| www.www免费av| 久久久久久久久中文| 少妇被粗大的猛进出69影院| 夜夜夜夜夜久久久久| 熟女电影av网| 满18在线观看网站| 国产熟女xx| 亚洲激情在线av| av片东京热男人的天堂| 每晚都被弄得嗷嗷叫到高潮| 免费看十八禁软件| 欧美乱色亚洲激情| 日本精品一区二区三区蜜桃| 亚洲第一电影网av| 2021天堂中文幕一二区在线观 | ponron亚洲| 69av精品久久久久久| 久久 成人 亚洲| 亚洲片人在线观看| 成年女人毛片免费观看观看9| 久久久久久免费高清国产稀缺| 国产黄色小视频在线观看| 91老司机精品| 巨乳人妻的诱惑在线观看| 两性夫妻黄色片| 亚洲一区中文字幕在线| 久久中文字幕一级| 热re99久久国产66热| 国产一卡二卡三卡精品| 99热6这里只有精品| 亚洲欧美精品综合一区二区三区| 国产黄色小视频在线观看| 免费在线观看日本一区| 久热这里只有精品99| 99久久国产精品久久久| 老司机靠b影院| 免费搜索国产男女视频| 国产成+人综合+亚洲专区| www.999成人在线观看| 亚洲男人的天堂狠狠| 91字幕亚洲| 午夜福利在线观看吧| 国产伦在线观看视频一区| 波多野结衣高清作品| 美女扒开内裤让男人捅视频| 亚洲国产毛片av蜜桃av| 亚洲成a人片在线一区二区| 不卡一级毛片| 日韩大码丰满熟妇| 亚洲欧美精品综合一区二区三区| av视频在线观看入口| 成人精品一区二区免费| 亚洲男人天堂网一区| 精品午夜福利视频在线观看一区| 村上凉子中文字幕在线| 日本撒尿小便嘘嘘汇集6| 亚洲精品中文字幕在线视频| 久久久久国内视频| 久久久久国产精品人妻aⅴ院| 久久久久久久午夜电影| 九色国产91popny在线| 精品国产一区二区三区四区第35| 国产单亲对白刺激| 亚洲美女黄片视频| 97超级碰碰碰精品色视频在线观看| 身体一侧抽搐| 桃红色精品国产亚洲av| 手机成人av网站| 国产精品av久久久久免费| 久久人妻av系列| 国产精品久久久av美女十八| 成人国语在线视频| 国产野战对白在线观看| 老司机午夜十八禁免费视频| 男人舔女人的私密视频| 国产欧美日韩一区二区精品| 亚洲中文字幕一区二区三区有码在线看 | 最好的美女福利视频网| 国产成人av教育| 999精品在线视频| 久久天躁狠狠躁夜夜2o2o| 老司机深夜福利视频在线观看| 亚洲av日韩精品久久久久久密| 国产又黄又爽又无遮挡在线| 久久天堂一区二区三区四区| 高清在线国产一区| 国产主播在线观看一区二区| 狂野欧美激情性xxxx| 97碰自拍视频| 日本一区二区免费在线视频| 午夜福利高清视频| 在线观看免费日韩欧美大片| 日韩中文字幕欧美一区二区| www.www免费av| 一本大道久久a久久精品| 午夜福利免费观看在线| 成人精品一区二区免费| videosex国产| 国产激情欧美一区二区| 亚洲男人天堂网一区| www.熟女人妻精品国产| 欧美黑人精品巨大| 法律面前人人平等表现在哪些方面| 在线视频色国产色| 久久草成人影院| 精品人妻1区二区| 一进一出抽搐动态| 韩国精品一区二区三区| 国产国语露脸激情在线看| 国产一区二区三区视频了| 日韩高清综合在线| 亚洲熟女毛片儿| av超薄肉色丝袜交足视频| 亚洲va日本ⅴa欧美va伊人久久| 听说在线观看完整版免费高清| 中文字幕av电影在线播放| 男女床上黄色一级片免费看| 精品国产乱子伦一区二区三区| 精品国产美女av久久久久小说| 亚洲avbb在线观看| 男人舔女人下体高潮全视频| 熟女少妇亚洲综合色aaa.| 国产成人系列免费观看| x7x7x7水蜜桃| 欧美日韩精品网址| 精品乱码久久久久久99久播| 日韩有码中文字幕| 亚洲中文日韩欧美视频| 波多野结衣av一区二区av| 色精品久久人妻99蜜桃| 黄色丝袜av网址大全| 不卡一级毛片| 嫩草影院精品99| 色综合欧美亚洲国产小说| 女同久久另类99精品国产91| 琪琪午夜伦伦电影理论片6080| 可以免费在线观看a视频的电影网站| 成人永久免费在线观看视频| 国内揄拍国产精品人妻在线 | 亚洲 欧美一区二区三区| 12—13女人毛片做爰片一| 两个人视频免费观看高清| 欧美久久黑人一区二区| 91字幕亚洲| 一夜夜www| 久久久久精品国产欧美久久久| 亚洲精品在线观看二区| 丝袜人妻中文字幕| 9191精品国产免费久久| 在线天堂中文资源库| 免费在线观看日本一区| 日韩三级视频一区二区三区| 麻豆一二三区av精品| 亚洲熟妇熟女久久| 黄色毛片三级朝国网站| 日韩视频一区二区在线观看| 国产成+人综合+亚洲专区| 久久中文字幕一级| 757午夜福利合集在线观看| 亚洲欧美一区二区三区黑人| 亚洲色图av天堂| 啪啪无遮挡十八禁网站| 国产精品亚洲一级av第二区| 又大又爽又粗| 欧美日韩黄片免| 男人操女人黄网站| 老熟妇仑乱视频hdxx| 久久久国产成人免费| 亚洲成人久久性| 麻豆成人午夜福利视频| 亚洲全国av大片| 久久这里只有精品19| 人人妻,人人澡人人爽秒播| 99久久精品国产亚洲精品| 精品人妻1区二区| e午夜精品久久久久久久| 免费看a级黄色片| 在线av久久热| 男女下面进入的视频免费午夜 | 午夜福利18| e午夜精品久久久久久久| 村上凉子中文字幕在线| 欧美日韩亚洲国产一区二区在线观看| 少妇被粗大的猛进出69影院| 淫妇啪啪啪对白视频| 国产av不卡久久| 国产激情偷乱视频一区二区| av在线天堂中文字幕| 欧美日本视频| 成人午夜高清在线视频 | 午夜精品久久久久久毛片777| 老熟妇仑乱视频hdxx| 国产国语露脸激情在线看| 人人妻,人人澡人人爽秒播| 欧美丝袜亚洲另类 | 黄色 视频免费看| 久99久视频精品免费| 日韩欧美国产一区二区入口| 成人国语在线视频| 国产主播在线观看一区二区| 国产在线精品亚洲第一网站| 在线天堂中文资源库| 精品久久久久久久人妻蜜臀av| 精品久久久久久久毛片微露脸| 欧美乱码精品一区二区三区| 满18在线观看网站| 久久香蕉精品热| 禁无遮挡网站| 亚洲专区中文字幕在线| 欧美中文日本在线观看视频| 久热这里只有精品99| 天堂影院成人在线观看| 天堂√8在线中文| 日日摸夜夜添夜夜添小说| or卡值多少钱| 在线十欧美十亚洲十日本专区| 69av精品久久久久久| 变态另类成人亚洲欧美熟女| 国产av又大| 嫩草影视91久久| e午夜精品久久久久久久| 欧美在线黄色| 国产人伦9x9x在线观看| 亚洲精品粉嫩美女一区| 亚洲专区字幕在线| 麻豆av在线久日| 精品一区二区三区视频在线观看免费| 女警被强在线播放| 国产成人系列免费观看| 正在播放国产对白刺激| 欧美激情极品国产一区二区三区| av电影中文网址| cao死你这个sao货| 国产99白浆流出| 女警被强在线播放| www日本黄色视频网| 亚洲国产精品久久男人天堂| 十八禁网站免费在线| 中文字幕av电影在线播放| 在线观看免费午夜福利视频| 亚洲全国av大片| 亚洲成国产人片在线观看| 日本 av在线| 男女那种视频在线观看| or卡值多少钱| 老司机在亚洲福利影院| 男人操女人黄网站| 久久久国产欧美日韩av| 国产精品 欧美亚洲| 天堂动漫精品| 欧美乱码精品一区二区三区| 国产精品99久久99久久久不卡| 久久九九热精品免费| 嫩草影视91久久| 精品欧美一区二区三区在线| 亚洲成人国产一区在线观看| 成人精品一区二区免费| 久久久久久久精品吃奶| 亚洲精品国产区一区二| aaaaa片日本免费| 国产精品香港三级国产av潘金莲| 国产精品一区二区免费欧美| 90打野战视频偷拍视频| 成年版毛片免费区| 午夜精品久久久久久毛片777| 国产三级在线视频| 国产高清videossex| 国产黄a三级三级三级人| 久久精品亚洲精品国产色婷小说| 免费看十八禁软件| 男人舔女人的私密视频| 校园春色视频在线观看| 国产成年人精品一区二区| 亚洲午夜理论影院| 国产精品 国内视频| 精品国产超薄肉色丝袜足j| 久久这里只有精品19| 久久狼人影院| 一本精品99久久精品77| 亚洲av电影在线进入| 激情在线观看视频在线高清| 国产熟女午夜一区二区三区| 国产成人系列免费观看| 日本一本二区三区精品| av在线天堂中文字幕| 国产精品98久久久久久宅男小说| 精品熟女少妇八av免费久了| 日韩 欧美 亚洲 中文字幕| 精品久久久久久久久久免费视频| 99国产综合亚洲精品| 欧美日韩精品网址| 可以在线观看的亚洲视频| 美女大奶头视频| 两个人看的免费小视频| 欧美日韩亚洲综合一区二区三区_| 久热爱精品视频在线9| 18禁黄网站禁片午夜丰满| 久久国产精品男人的天堂亚洲| 国产精品一区二区免费欧美| 久久人妻福利社区极品人妻图片| 在线观看免费视频日本深夜| 亚洲人成伊人成综合网2020| 丰满的人妻完整版| 精品久久久久久久人妻蜜臀av| 国产成人精品久久二区二区91| 999久久久精品免费观看国产| 久久久久久亚洲精品国产蜜桃av| 黄片小视频在线播放| 亚洲午夜理论影院| 又大又爽又粗| 免费看a级黄色片| 级片在线观看| 国产精品精品国产色婷婷| 国产精品一区二区精品视频观看| bbb黄色大片| 日韩精品免费视频一区二区三区| 亚洲精品av麻豆狂野| 国产成人精品无人区| 欧美亚洲日本最大视频资源| 欧美一级毛片孕妇| 亚洲精品美女久久av网站| 欧美精品啪啪一区二区三区| 一级作爱视频免费观看| 男女下面进入的视频免费午夜 | 黄片小视频在线播放| 在线观看舔阴道视频| 又大又爽又粗| 午夜影院日韩av| 精品第一国产精品| 亚洲九九香蕉| 国产一级毛片七仙女欲春2 | 深夜精品福利| 亚洲va日本ⅴa欧美va伊人久久| 欧美乱妇无乱码| 精品久久久久久,| 国产成人av激情在线播放| 制服丝袜大香蕉在线| 亚洲电影在线观看av| 久久久国产成人精品二区| 深夜精品福利| 色婷婷久久久亚洲欧美| 最近在线观看免费完整版| 巨乳人妻的诱惑在线观看| 亚洲 欧美一区二区三区| cao死你这个sao货| 嫩草影视91久久| 级片在线观看| 91av网站免费观看| 色哟哟哟哟哟哟| 成人亚洲精品av一区二区| 国产三级黄色录像| 侵犯人妻中文字幕一二三四区| 久久久精品国产亚洲av高清涩受| 久久精品国产清高在天天线| 国产熟女xx| 亚洲成av人片免费观看| 成人三级黄色视频| 黄片播放在线免费| 午夜福利免费观看在线| 国产精品乱码一区二三区的特点| 青草久久国产| 精品国产一区二区三区四区第35| 亚洲在线自拍视频| 亚洲av电影在线进入| 国产精品乱码一区二三区的特点| tocl精华| 一区福利在线观看| 99久久99久久久精品蜜桃| 一本久久中文字幕| 成人三级黄色视频| 中文字幕av电影在线播放| 最近最新免费中文字幕在线| 国产av在哪里看| 黄片播放在线免费| 日韩av在线大香蕉| 99国产极品粉嫩在线观看| 观看免费一级毛片| 满18在线观看网站| 此物有八面人人有两片| 制服人妻中文乱码| 国产亚洲精品久久久久5区| 少妇粗大呻吟视频| 精品久久久久久久久久久久久 | 一区二区三区国产精品乱码| 色精品久久人妻99蜜桃| 色在线成人网| 黄色毛片三级朝国网站| www.www免费av| 18禁黄网站禁片免费观看直播| 免费一级毛片在线播放高清视频| 中国美女看黄片| 亚洲国产精品999在线| 中文字幕av电影在线播放| www.www免费av| 啦啦啦免费观看视频1| 久久久久久久久中文| 成人手机av| 久久亚洲精品不卡| 美国免费a级毛片| 欧美人与性动交α欧美精品济南到| 黄色女人牲交| 美国免费a级毛片| 亚洲成人久久爱视频| 国产精品,欧美在线| 国产精品乱码一区二三区的特点| 99久久国产精品久久久|