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

    基于SRTM地形數(shù)據(jù)的天氣雷達(dá)電磁波 非常規(guī)遮擋回波補(bǔ)償技術(shù)研究

    2021-11-06 10:10:54胡啟元戚友存王楠朱自偉
    關(guān)鍵詞:仰角方位角電磁波

    胡啟元 戚友存 王楠 朱自偉

    (1 陜西省氣象臺,西安 710014;2 秦嶺和黃土高原生態(tài)環(huán)境氣象重點(diǎn)實(shí)驗(yàn)室,陜西省氣象局,西安 710014;3 中國科學(xué)院地理科學(xué)與資源研究所,陸地水循環(huán)實(shí)驗(yàn)室,北京 100101;4 中國科學(xué)院大學(xué),北京 100864)

    0 引言

    基于多普勒天氣雷達(dá)觀測數(shù)據(jù)的定量降水估測(QPE)產(chǎn)品是0~2小時定量短時臨近(以下簡稱“短臨”)降水預(yù)報(QPF)和強(qiáng)降水短臨預(yù)警的重要基礎(chǔ),但在日常業(yè)務(wù)中雷達(dá)觀測數(shù)據(jù)往往存在缺失或不完整現(xiàn)象,其原因一方面是由于雷達(dá)站點(diǎn)附近存在海拔高于雷達(dá)高度的復(fù)雜地形,其低仰角電磁波束可能會部分或全部被某處地形所遮擋,另一方面由于城市建設(shè)程度日益加快,雷達(dá)站周邊增加諸多如高樓大廈、通信塔或綠化樹木等遮擋雷達(dá)電磁波的人造建筑物,使得雷達(dá)電磁波束發(fā)生遮擋,造成雷達(dá)反射率測量出現(xiàn)偏低或完全遮擋,進(jìn)而導(dǎo)致定量降水估測產(chǎn)品降水率的低估或缺測。因此,需要通過研究雷達(dá)無縫混合仰角計(jì)算方法及雷達(dá)非常規(guī)遮擋訂正技術(shù),來正確計(jì)算雷達(dá)QPE所需的雷達(dá)回波反射率場。

    最低混合有效仰角掃描反射率是指由雷達(dá)各方位角上未被嚴(yán)重遮擋的最低有效仰角反射率因子所組成的反射率平面分布,它能夠比較準(zhǔn)確地反映出最接近地面的真實(shí)降水強(qiáng)度和空間分布。美國國家強(qiáng)風(fēng)暴實(shí)驗(yàn)室(NSSL)Fulton等在1997年研究WSR-88D單偏振天氣雷達(dá)降水估計(jì)算法流程時利用數(shù)字地面模型(DTM)數(shù)據(jù)設(shè)計(jì)了美國單偏振雷達(dá)混合仰角反射率計(jì)算方法,通過對輸入反射率進(jìn)行基于地形的混合仰角掃描預(yù)處理,盡可能降低地物雜波出現(xiàn)頻次以及由于地形阻擋造成的降水反射率數(shù)據(jù)丟失情況;Bech等利用西班牙巴塞羅那地區(qū)探空觀測數(shù)據(jù)模擬雷達(dá)在低層大氣中的傳播特征,發(fā)現(xiàn)在不同雷達(dá)波束傳播條件下采用對應(yīng)的觀測或預(yù)報的垂直折射率梯度(VRG)信息,能夠有效校正單偏振雷達(dá)電磁波束被地形部分遮擋情況;Kucera等利用地理信息系統(tǒng)(GIS)和高分辨率數(shù)字高程模型(DEM)進(jìn)行部分或完全受遮擋雷達(dá)波束的識別和能量損失校正;Langston等介紹了美國國家強(qiáng)風(fēng)暴實(shí)驗(yàn)室開發(fā)的能夠根據(jù)地形、雷達(dá)功率密度函數(shù)和波束傳播路徑來計(jì)算雷達(dá)混合掃描仰角的HybScan自動算法,并利用該算法及1弧秒DEM數(shù)據(jù)為全美所有WSR-88D雷達(dá)制作了波束遮擋區(qū)域及混合仰角參數(shù)集;Zhang等在對美國Multi-Radar Multi-Sensor (MRMS)系統(tǒng)定量降水估測應(yīng)用研究中介紹了無縫混合掃描反射率高度(SHSRH)概念及其在MRMS系統(tǒng)中對垂直反射率廓線(VPR)的校正作用,而后者是近十年國外用來改進(jìn)定量降水估測產(chǎn)品的效果較好的方法之一;Tang等針對雷達(dá)電磁波受到的非常規(guī)遮擋情況設(shè)計(jì)制作了三種遮擋修訂方法和用于業(yè)務(wù)使用的非常規(guī)遮擋修訂參數(shù)方案,并且已經(jīng)在美國本土130多部雷達(dá)進(jìn)行應(yīng)用,評估結(jié)果表明該方案能夠有效減少和消除美國降水產(chǎn)品中各種受遮擋而不連續(xù)情況。國內(nèi)也有部分學(xué)者利用混合仰角反射率進(jìn)行雷達(dá)觀測受遮擋修訂的改進(jìn)工作。肖艷姣等在2008年利用1:25萬的DEM數(shù)據(jù)計(jì)算湖北SB雷達(dá)波束阻擋率,以60%阻擋率為阻擋閾值得出湖北雷達(dá)混合掃描仰角參數(shù),并通過等射束高度拼圖數(shù)據(jù)計(jì)算混合反射率拼圖產(chǎn)品,有效解決了反射率受地物阻擋問題,但在受非地形阻擋區(qū)域回波仍明顯偏弱;戚友存等2011年在研究不同雷達(dá)垂直反射率廓線(VPR)訂正方案對雷達(dá)定量降水估測準(zhǔn)確率的改進(jìn)作用時設(shè)計(jì)了基于VPR訂正方法的QPE算法,該算法提出了一個將兩個相鄰未遮擋仰角上的反射率數(shù)據(jù)進(jìn)行融合來消除遮擋不連續(xù)的無縫混合反射率(SHSR)計(jì)算方法;楊瀧等利用杭州雷達(dá)受遮擋區(qū)相鄰的無遮擋區(qū)的平均VPR遮擋模型,將高仰角無遮擋反射率因子線性插值到低仰角受遮擋區(qū)域,填補(bǔ)前、后結(jié)果評估結(jié)果表明填補(bǔ)后QPE降水量與觀測量有很好一致性,填補(bǔ)后降水估算效果優(yōu)于填補(bǔ)前;勾亞彬等通過分析浙江6部雷達(dá)觀測數(shù)據(jù)概率特征,提出了一個不依賴地形數(shù)據(jù)以及大氣狀態(tài)的基于概率統(tǒng)計(jì)的雷達(dá)部分遮擋識別方法,在雷達(dá)組網(wǎng)拼圖計(jì)算時應(yīng)用該方法,能夠有效選擇在同一區(qū)域中受遮擋較輕的雷達(dá)回波反射率因子進(jìn)行拼圖,增強(qiáng)雷達(dá)重疊區(qū)域內(nèi)的組網(wǎng)拼圖數(shù)據(jù)的質(zhì)量。

    目前國內(nèi)外雷達(dá)業(yè)務(wù)系統(tǒng)中的雷達(dá)估測降水產(chǎn)品例如美國早期StageⅡ及StageⅣ產(chǎn)品和目前MRMS系統(tǒng)Q3產(chǎn)品、中國全國組網(wǎng)定量降水(MQPE)產(chǎn)品和SWAN系統(tǒng)QPE產(chǎn)品等都已運(yùn)用混合掃描反射率生成技術(shù)進(jìn)行QPE業(yè)務(wù)改進(jìn)工作,但與美國等發(fā)達(dá)國家相比我國仍存在一定差距。美國Q3產(chǎn)品針對遮擋不連續(xù)情況不僅采用基于地形和無縫處理結(jié)合方法,并且通過計(jì)算融化層高度和訂正雷達(dá)垂直廓線來找到存在有效反射率數(shù)據(jù)的最低有效仰角,進(jìn)一步提高混合反射率精準(zhǔn)度。國內(nèi)2005年研發(fā)的MQPE產(chǎn)品雖然應(yīng)用了混合反射率訂正技術(shù)但沒有利用地形信息,只是固定采集0.5°~3.4°四個仰角反射率,若遇到低仰角地形阻擋時采用阻擋物后方距離庫上的反射率數(shù)值,因此結(jié)果存在較大誤差。SWAN系統(tǒng)QPE產(chǎn)品較MQPE產(chǎn)品增加了不同地形數(shù)據(jù)對混合反射率生成的應(yīng)用,但未進(jìn)行反射率無縫技術(shù)處理,而且QPE業(yè)務(wù)產(chǎn)品仍存在由于非地形因素導(dǎo)致部分地區(qū)反射率不連續(xù)情況。因此本文研究針對因人造建筑物、樹木等非地形因素造成的雷達(dá)電磁波受遮擋不連續(xù)現(xiàn)象,形成基于最優(yōu)無縫混合掃描仰角的遮擋修訂參考方案,通過對典型個例或汛期降水過程的定量、定性統(tǒng)計(jì)特征來分析參考方案的應(yīng)用效果,解決陜西雷達(dá)因觀測被遮擋而造成定量降水估測產(chǎn)品缺失或錯誤問題,并在陜西短臨預(yù)報業(yè)務(wù)系統(tǒng)中進(jìn)行業(yè)務(wù)應(yīng)用。

    1 數(shù)據(jù)與方法

    1.1 選取數(shù)據(jù)

    本文使用數(shù)據(jù)主要包括DEM地形數(shù)據(jù)、雷達(dá)體掃觀測和地面雨量計(jì)降水觀測數(shù)據(jù)。DEM地形資料由美國航空航天局(NASA)和國防部國家測繪局(NIMA)聯(lián)合測量的中國境內(nèi)航天飛機(jī)雷達(dá)地形測繪使命(Shuttle Radar Topography Mission,SRTM)數(shù)據(jù),分辨率為90 m×90 m;雷達(dá)體掃觀測數(shù)據(jù)選取2017—2019年5—10月陜西7部C波段業(yè)務(wù)雷達(dá)逐6分鐘實(shí)時基數(shù)據(jù)(圖1a),雷達(dá)最遠(yuǎn)探測距離400 km,反射率和徑向速度分辨率分別為0.5和0.125 km/h,采用VCP21體掃模式,9個有效仰角(0.5°,1.5°,2.4°,3.4°,4.3°,6.0°,9.9°,14.6°,19.5°);地面雨量計(jì)降水觀測數(shù)據(jù)選取2017—2019年5—10月陜西1872個加密氣象觀測站逐小時降水量數(shù)據(jù)(圖1b)。

    圖1 陜西雷達(dá)站點(diǎn)分布與地形圖(a)和地面雨量計(jì)觀測站分布圖(b) Fig. 1 Distribution of radar stations and topographic (a) and distribution of surface rainfall observation stations (b) in Shaanxi Province

    1.2 雷達(dá)電磁波非常規(guī)遮擋回波補(bǔ)償技術(shù)

    多普勒雷達(dá)低層仰角電磁波束被遮擋原因主要包括因復(fù)雜地形遮擋和地物(樹木、高樓、高速路、高壓線電塔等人造建筑物)遮擋,本文將地形遮擋稱為雷達(dá)電磁波常規(guī)遮擋因素(STANDARD BLOCKAGE,SB),地物(非地形遮擋)遮擋稱為雷達(dá)電磁波非常規(guī)遮擋(NON-STANDARD BLOCKAGE,NSB)。為解決這兩類因素對雷達(dá)電磁波造成的遮擋問題,本文基于陜西省7部多普勒天氣雷達(dá),研發(fā)了雷達(dá)電磁波遮擋訂正算法(圖2),首先利用2017—2019年5—10月逐6分鐘雷達(dá)原始基數(shù)據(jù)通過

    Z

    -

    R

    關(guān)系計(jì)算得出雷達(dá)各仰角的累積QPE空間分布,結(jié)合由陜西DEM地形數(shù)據(jù)及雷達(dá)不同仰角波束高度信息形成的最低有效仰角分布,計(jì)算基于地形的混合掃描仰角的累積QPE空間分布;其次在單雷達(dá)QPE累積分布上查找識別出NSB區(qū)域,根據(jù)NSB區(qū)域特征和大小研究制定三種修訂方案進(jìn)行分類修訂處理,并形成各雷達(dá)NSBM參考方案,將參考方案再應(yīng)用到QPE累積分布修訂流程中,得出基于地形遮擋和NSB共同修訂的最優(yōu)無縫混合仰角反射率及QPE產(chǎn)品;最后通過在典型降水個例和所有降水過程中,對各雷達(dá)的最優(yōu)無縫混合仰角反射率QPE產(chǎn)品與觀測范圍內(nèi)的地面雨量計(jì)實(shí)況資料進(jìn)行定性、定量統(tǒng)計(jì)分析,評估修訂算法改進(jìn)優(yōu)劣和算法性能。

    圖2 陜西多普勒天氣雷達(dá)電磁波非常規(guī)遮擋區(qū)域回波補(bǔ)償方案技術(shù)流程示意圖 Fig. 2 Flow chart of compensated scheme for electromagnetic beam non-standard blockage area in Shaanxi Doppler weather radar.

    (1)基于陜西DEM地形數(shù)據(jù)計(jì)算雷達(dá)最低有效仰角

    計(jì)算雷達(dá)混合掃描仰角反射率的前提是計(jì)算出雷達(dá)不同仰角上電磁波受遮擋的百分比;再根據(jù)計(jì)算得到的雷達(dá)電磁波遮擋百分比,確定雷達(dá)在各方位角上的最低有效仰角。本研究中,把單仰角電磁波束被地形遮擋的電磁波能量占總能量的百分比稱為遮擋率,并且在雷達(dá)不同方位角上選取的雷達(dá)波束被遮擋率小于50%的最低仰角作為該方位角上的最低有效仰角,由此得出該雷達(dá)掃描平面內(nèi)最低有效仰角分布。

    (2)NSB區(qū)域的識別定位

    基于最低有效仰角分布,計(jì)算得到的雷達(dá)最低混合仰角反射率分布可以很好地改善由于復(fù)雜地形造成的SB情況。但在實(shí)際業(yè)務(wù)應(yīng)用中仍存在由于地物導(dǎo)致的NSB情況,因此需要對NSB區(qū)域進(jìn)行識別定位并進(jìn)行訂正。通常NSB區(qū)域可以從雷達(dá)各仰角單時次混合反射率PPI分布上識別出來,但由于有些情況下雷達(dá)電磁波遮擋不嚴(yán)重,通過單時次雷達(dá)反射率很難識別出來。因此本研究中,采用長時段的雷達(dá)混合反射率對應(yīng)的QPE累積分布(采用

    Z

    =300×

    R

    固定

    Z

    -

    R

    關(guān)系)來識別NSB區(qū)域。

    (3)NSB區(qū)域修訂方法

    識別出雷達(dá)NSB區(qū)域位置信息后,根據(jù)NSB區(qū)域不同特征采用三種修訂方法進(jìn)行遮擋不連續(xù)修訂。

    ① 針對方位角寬度小于5°的NSB區(qū)域,采用交叉方位角插值法進(jìn)行處理,算法見式(1)。造成這種方位角寬度<5°NSB區(qū)域的原因通常是由于孤立的高層建筑或高大樹木阻擋了低層仰角反射率,并且相鄰NSB區(qū)域兩側(cè)的電磁波束能量沒有受到任何遮擋影響。通過該方法進(jìn)行遮擋修正后反射率或降水分布圖像更加平滑,降水系統(tǒng)表現(xiàn)更加完整,但前提是假設(shè)降水回波強(qiáng)度在NSB區(qū)域的水平變化較小。

    其中:

    A

    A

    為NSB區(qū)域的上下邊界方位角,

    A

    A

    A

    之間的某個方位角, 為方位角

    A

    上的反射率值。

    ② 針對方位角寬度大于5°的NSB區(qū)域,采用更高一層仰角反射率數(shù)據(jù)替代NSB區(qū)域數(shù)據(jù)方法,算法見式(2)。該類NSB在地物遮擋情況中占大多數(shù),當(dāng)出現(xiàn)該類NSB情況時說明該區(qū)域雷達(dá)低層電磁波束被高大建筑群或茂密植被帶嚴(yán)重遮擋,并且由于無法忽視降水系統(tǒng)的大范圍水平空間變化而不能采用方法①進(jìn)行修訂,因此需要采用更高一層接近地面的反射率場進(jìn)行填補(bǔ)替代。通過該方法可以有效地修正雷達(dá)大范圍反射率或降水分布存在的遮擋不連續(xù),但前提是在NSB區(qū)域中降水系統(tǒng)內(nèi)部垂直變化較小,例如在取得的更高一層反射率場上不存在明顯亮帶等。

    ③ 混合仰角反射率場進(jìn)行SB處理和方法②修訂后,可能還會存在由于相鄰方位角采用不同仰角反射率而造成的回波強(qiáng)度不連續(xù)梯度,需要采用線性加權(quán)滑動平均法對其進(jìn)行平滑處理,算法見式(3)。該方法并不直接對NSB區(qū)域進(jìn)行修正,而是對基于DEM數(shù)據(jù)的標(biāo)準(zhǔn)混合仰角掃描反射率或者方法②生成出來的混合仰角反射率進(jìn)行非自然邊界平滑處理。

    其中:

    A

    為NSB區(qū)域邊界的方位角,

    i

    A

    相鄰方位角, -

    N

    i

    N

    ,選取

    N

    =5(方位角

    A

    前后5個相鄰點(diǎn)),

    W

    為線性權(quán)重,

    R

    為方位角

    A

    上對應(yīng)仰角的反射率值。

    (4)定量統(tǒng)計(jì)分析方法

    為評估遮擋修訂前后算法改進(jìn)優(yōu)劣,選取典型個例和汛期降水過程進(jìn)行均方根誤差(

    RMSE

    )、平均絕對誤差(

    MAE

    )和相關(guān)系數(shù)(

    CC

    )(公式4~6)統(tǒng)計(jì)特征分析,對基于最優(yōu)無縫混合掃描仰角反射率的定量降水估測產(chǎn)品與地面雨量計(jì)進(jìn)行修訂算法評估。

    其中:

    r

    為雷達(dá)定量估測降水量,

    g

    為地面雨量計(jì)觀測降水量。

    2 基于最優(yōu)無縫混合掃描仰角的雷達(dá)遮擋區(qū)域修訂算法

    2.1 陜西雷達(dá)受常規(guī)遮擋(SB)情況及改進(jìn)結(jié)果評估

    陜西7部CB波段多普勒天氣雷達(dá)(榆林、延安、西安、寶雞、漢中、安康、商洛)站點(diǎn)布設(shè)位置附近均受到不同程度的地形遮擋情況,其中商洛雷達(dá)低層電磁波束受遮擋最嚴(yán)重,榆林雷達(dá)受遮擋程度最小。以商洛雷達(dá)為例,圖3a為商洛雷達(dá)2019年5月—10月0.5°~3.4°基于QPE累積降水的各方位角上最大徑向探測距離(

    L

    )分布,可以看出雷達(dá)0.5°、1.5°仰角上

    L

    >300 km區(qū)域占總區(qū)域9%和31%, 在2.4°仰角上

    L

    >200 km和3.4°仰角上

    L

    >150 km區(qū)域占比分別提高至61%和78%,可見雷達(dá)有效可探測區(qū)域隨著仰角抬高迅速增大,大部分地形遮擋情況存在1.5°仰角以下,但在商洛雷達(dá)東北部(約25°~80°)和西北部(約330°~350°)在3.4°仰角上仍存在明顯遮擋不連續(xù)區(qū)域,直到抬高至6.0°仰角才能明顯改善(圖略),這種復(fù)雜地形環(huán)境造成無法直接使用0.5°或1.5°仰角反射率進(jìn)行定量降水估測,需要進(jìn)行混合仰角處理后才能計(jì)算出雷達(dá)最接近地面的有效降水回波。

    圖3 商洛雷達(dá)2019年5—10月QPE累積降水各方位角上最大徑向距離分布 (a)0.5°—3.4°仰角;(b)混合仰角 Fig. 3 Maximum radial distance distribution of accumulated QPE in Shangluo radar from May to October 2019 (a) 0.5°to 3.4°elevation angle, (b) hybrid elevation angle

    通過最低有效仰角參數(shù)選取方法計(jì)算得出的雷達(dá)基于地形混合仰角參數(shù)分布(圖4),可以看出除榆林雷達(dá)大部地區(qū)可直接采用0.5°仰角反射率進(jìn)行QPE計(jì)算以外,其他雷達(dá)都需要取更高一層或以上仰角:延安雷達(dá)北部、西南部、東南部和西北部局地需采用1.5°仰角;西安雷達(dá)東南部、西南部、西北部及寶雞雷達(dá)西北部、西部需采用1.5°仰角,南部需采用2.4°仰角才能觀測到有效反射率;陜南三市由于地處秦巴山區(qū)中,低層仰角反射率受遮擋更加嚴(yán)重,其中安康雷達(dá)大部地區(qū)需采用2.4°仰角,南部受山峰遮擋需采用4.3°仰角,漢中雷達(dá)除東部及西部局地以外均需采用1.5°以上仰角,商洛雷達(dá)各方位角上采用了從1.5°~3.4°差別明顯的混合仰角,并且在漢中北部和商洛西北部、東北部有效最低仰角高度達(dá)到了4.3°或6.0°仰角。通過各雷達(dá)基于地形的混合仰角參數(shù)分布,能夠有效訂正受地形遮擋情況,對比圖3a可以看出,商洛雷達(dá)進(jìn)行混合仰角處理后其5—10月QPE累積降水(圖3b)在大部分方位角上都能觀測到有效回波,尤其在西北部受遮擋區(qū)域改進(jìn)效果明顯。

    圖4 陜西7部雷達(dá)基于地形資料的混合仰角參數(shù)分布 Fig. 4 Distribution of hybird elevation angle parameters of seven radars based on terrain data in Shaanxi Province

    雖然理論上通過基于地形的混合仰角參數(shù)能夠有效解決雷達(dá)低層降水回波受遮擋情況,但實(shí)際業(yè)務(wù)中計(jì)算出的混合仰角反射率分布仍存在部分區(qū)域缺失和局地不連續(xù)現(xiàn)象。圖5給出了安康雷達(dá)2019年5—10月在進(jìn)行SB情況修訂前的2.4°仰角QPE累積分布(0.5°和1.5°仰角可觀測數(shù)據(jù)過少故不采用)和進(jìn)行SB情況修訂后的混合仰角QPE累積分布,可以看出采用混合仰角處理比直接采用2.4°仰角能夠觀測到更多降水,南部、東北部原本缺失的降水區(qū)域得到有效填補(bǔ),但修訂后的混合仰角QPE分布依然存在大量不連續(xù)區(qū)域:1) 西南部、偏西部存在方位角寬度約5°~10°的空白區(qū)域,其原因是由于過時DEM地形數(shù)據(jù)或者密集建筑群造成;2) 多處方位角寬度<5°的分散的楔形不連續(xù)區(qū)域,多是由于孤立建筑物或高大樹木造成;3) 2.4°仰角QPE分布上部分區(qū)域如偏南部(如170°或190°)雖然有降水?dāng)?shù)據(jù),但明顯低于附近方位角降水量,可能是由于電磁波束被地形部分遮擋導(dǎo)致觀測數(shù)據(jù)偏低造成,這一部分本文不做討論,按照遮擋不連續(xù)情況進(jìn)行修訂。

    圖5 2019年5—10月安康雷達(dá)QPE累積分布 (a)2.4°仰角QPE分布;(b)經(jīng)過SB修訂后的混合仰角QPE分布 Fig. 5 QPE cumulative distribution of Ankang radar from May to October 2019 (a) QPE distribution at 2.4°elevation, (b) Hybird elevation QPE distribution revised by SB scheme

    2.2 陜西雷達(dá)受非常規(guī)遮擋(NSB)區(qū)域修訂方法及改進(jìn)結(jié)果評估

    對雷達(dá)NSB區(qū)域進(jìn)行修訂之前需要按照NSB修訂方法進(jìn)行分類,然后分區(qū)域進(jìn)行處理。由于雷達(dá)各體掃模式的起始體掃高度存在差異,對不同體掃模式下需采用對應(yīng)的修訂處理,陜西CB雷達(dá)普遍采用VCP21模式,因此本文針對7部雷達(dá)分別制作適用于VCP21模式的NSB修訂參考方案,實(shí)現(xiàn)在業(yè)務(wù)應(yīng)用中能夠?qū)崟r自動修訂而不需要再進(jìn)行人工干預(yù)。利用經(jīng)過SB修訂的陜西7部雷達(dá)2019年5—10月逐6分鐘QPE累積分布,計(jì)算確定出NSB區(qū)域位置信息,按順時針方向定義

    A

    為NSB區(qū)域起始方位角,

    A

    為NSB區(qū)域終止方位角,

    R

    為起始徑向距離,

    R

    為終止徑向距離,根據(jù)2.2節(jié)NSB修訂方法,分類采用不同執(zhí)行方案,其中針對方法2在NSB修訂參考方案上定義EleNum為所取更高仰角的數(shù)量(0代表采用其他方法,1代表需在原仰角上再抬高一個仰角,以此類推),由此可以得出每部雷達(dá)的NSB修訂參考方案(表1,以安康為例)。

    表1 安康雷達(dá)NSB修訂參考方案 Table 1 Reference scheme for NSB revision of Ankang radar

    按照安康雷達(dá)NSB修訂參考方案,繼續(xù)對2.1節(jié)經(jīng)過SB修訂的QPE累積分布進(jìn)行NSB修訂,得出最終基于最優(yōu)無縫混合掃描仰角遮擋修訂算法的定量降水估測產(chǎn)品(圖6),結(jié)果表明修訂效果明顯,修訂后水平方向不連續(xù)區(qū)域明顯減少,可觀測到回波降水范圍基本覆蓋所有地區(qū)。通過安康雷達(dá)2019年6月5日8時—11時逐小時混合反射率在NSB修訂前、后對比結(jié)果來分析改進(jìn)效果,圖7a、7c、7e、7g為未經(jīng)過NSB修訂混合反射率分布,隨著回波從西南向東北發(fā)展,雷達(dá)西北部和西南部狹長的楔形遮擋不連續(xù)現(xiàn)象非常明顯,而對應(yīng)時次的圖7b、7d、7f、7h中上述不連續(xù)區(qū)域經(jīng)過NSB修訂得到了明顯改進(jìn),使整個過程中混合仰角反射率表現(xiàn)出了較好的連續(xù)和平滑特征。值得注意的是圖7c和7d上09時次在雷達(dá)東南部120°~135°附近回波相比緊鄰兩側(cè)回波強(qiáng)度偏弱,如果只從單時次的回波特征來看可能是由于電磁波束部分被遮擋造成的回波強(qiáng)度不連續(xù),但從下幾個時次(10時和11時)可以看出隨著回波移動該區(qū)域回波強(qiáng)度與兩側(cè)基本一致,實(shí)際上并不存在不連續(xù)情況。上述現(xiàn)象說明從單時次反射率或QPE降水分布上無法準(zhǔn)確識別出NSB區(qū)域,因此需要采用長時段反射率分布或累積QPE分布來進(jìn)行識別遮擋不連續(xù)位置。

    圖6 經(jīng)過NSB修訂后的安康雷達(dá)2019年5—10月QPE累積分布 Fig. 6 QPE cumulative distribution of Ankang radar from May to October 2019 revised by NSB scheme

    圖7 安康雷達(dá)2019年6月5日08—11時逐小時修訂前混合仰角反射率(a,c,e,g)和修訂后混合仰角反射率(b,d,f,h)(單位:dBz) Fig. 7 The comparison of Ankang radar hourly hybird reflectivity before (a,c,e,g) and after (b,d,f,h) revision from 08:00 to 11:00 on 5 June 2019 (unit: dBz)

    3 修訂效果評估

    3.1 降水個例改進(jìn)效果

    從2.2節(jié)安康一次降水過程雷達(dá)反射率演變的定性評估可以看出,經(jīng)過最優(yōu)無縫混合掃描仰角處理遮擋不連續(xù)后的降水回波覆蓋面和平滑程度較改進(jìn)前有明顯提高,為了進(jìn)一步定量評估改進(jìn)前、后效果,選取了延安、寶雞、西安、漢中、安康、商洛6部雷達(dá)在2017年—2019年6次對流降水或穩(wěn)定降水過程進(jìn)行效果評估。選取的對流降水過程包括延安2019年7月21—22日、漢中2019年8月9日,穩(wěn)定降水過程包括西安2019年9月14日、寶雞2017年6月3—4日、商洛2017年6月4日、安康2018年6月18日。

    對上述過程以均方根誤差(RMSE)、平均絕對誤差(MAE)和相關(guān)系數(shù)(CC)為檢驗(yàn)指標(biāo),以使用固定

    Z

    =300×

    R

    關(guān)系得出的雷達(dá)估測降水量與對應(yīng)的地面降水量作為檢驗(yàn)要素,以各雷達(dá)存在的非常規(guī)遮擋區(qū)域作為評估空間范圍,進(jìn)行逐小時混合反射率修訂前、后評估對比分析。通過圖8可以看出,各雷達(dá)非常規(guī)遮擋區(qū)域內(nèi)大部分時間RMSE、MAE、CC在修訂后比修訂前表現(xiàn)更好,其中漢中、延安、安康三市改進(jìn)效果最為明顯,平均RMSE從修訂前6.214減小到4.236,平均MAE由0.861降低至0.627,相關(guān)系數(shù)CC由0.238提高至0.514,非常規(guī)遮擋區(qū)域的QPE雨量與地面雨量偏差離散度得到有效抑制,兩者相關(guān)性提高53%。評估結(jié)果表明在固定

    Z

    -

    R

    關(guān)系下,經(jīng)過NSB修訂后的雷達(dá)估測降水量與地面降水量之間的誤差更小,估測降水分布比修訂前更接近于地面降水分布,同時在延安和漢中對流降水過程中相同

    Z

    -

    R

    關(guān)系下修訂后的RMSE、MAE和CC比修訂前有更好表現(xiàn),這說明NSB修訂處理能有效改善強(qiáng)對流降水過程中由于遮擋導(dǎo)致強(qiáng)降水回波缺失現(xiàn)象,使得更多與地面降水量線性相關(guān)度高的降水回波被填補(bǔ)進(jìn)混合反射率中。

    圖8 各雷達(dá)在典型個例過程中均方根誤差(a)、平均絕對誤差(b)和相關(guān)系數(shù)(c)在修訂前、后對比情況 (實(shí)線為修訂前,虛線為修訂后) Fig. 8 The comparison before and after revision of the root mean square error (a), mean absolute error (b) and correlation coefficient (c) of each radar (The solid line is before revision, the dotted line is after revision)

    為更客觀評估雷達(dá)估測降水的改進(jìn)效果,將延安、漢中、安康個例過程中回波覆蓋面積最大時刻(延安7月21日23時、漢中8月9日4時、安康6月18日0時)的修訂前、后估測降水量與實(shí)況降水量進(jìn)行定量對比。圖9a、9b、9c顯示了修訂前、后發(fā)生變化的估測降水站點(diǎn),修訂后站點(diǎn)QPE降水更趨近分布在等比例線附近,與地面雨量線性相關(guān)程度更好,對應(yīng)的趨勢線(黑色虛線為修訂后,灰色虛線為修訂前)在修訂后擬合程度更高;QPE偏小而地面雨量偏大的點(diǎn)顯著減小,說明修訂后的降水回波在不受遮擋影響情況下更接近理想情況下代表地面降水的最低降水回波。各雷達(dá)在修訂后可觀測到QPE站點(diǎn)數(shù)超過修訂前站點(diǎn)數(shù),進(jìn)一步說明原本被遮擋的降水回波在進(jìn)行NSB修訂后可以被觀測到。

    圖9 延安、漢中、安康雷達(dá)1小時雷達(dá)估計(jì)降水量與地面降水量比較 (a)延安7月21日23時;(b)漢中8月9日4時;(c)安康6月18日0時 (灰方塊為修訂前結(jié)果,黑圓點(diǎn)為修訂后結(jié)果;灰虛線為修訂前R2趨勢線,黑虛線為修訂后R2趨勢線) Fig. 9 The comparison of 1-hour radar estimated precipitation with surface precipitation in Yan’an, Hanzhong and Ankang (a) 23:00 BT, July 21, Yan’an, (b) 04:00 BT on August 9, Hanzhong, (c) 00:00 BT, June 18, Ankang (the gray box is the result before revision, the black dot is the result after revision, the gray dotted line is the R2 trend line before revision, the black dotted line is the R2 trend line after revision)

    3.2 所有降水過程定量評估

    對2019年5月—10月6部雷達(dá)觀測范圍內(nèi)所有降水過程進(jìn)行統(tǒng)計(jì)分析,同樣以各雷達(dá)存在的非常規(guī)遮擋區(qū)域作為評估空間范圍,各雷達(dá)長時段序列統(tǒng)計(jì)結(jié)果(圖略)顯示在穩(wěn)定性降水過程中修訂后RMSE明顯優(yōu)于修訂前,在對流性降水過程修訂后RMSE有所減小但不明顯;MAE統(tǒng)計(jì)結(jié)果顯示經(jīng)過遮擋修訂后大部分降水過程中部分時段的高誤差現(xiàn)象得到了有效抑制,并且約有75%降水過程MAE較修訂前降低;相關(guān)系數(shù)CC由修訂前的69.7%降水過程為正相關(guān)提升至修訂后的93.5%降水過程為正相關(guān)。表2為全部降水過程計(jì)算得出的修訂前、后平均統(tǒng)計(jì)結(jié)果,漢中平均RMSE和平均MAE減小最多,改進(jìn)提高效果最明顯,延安平均CC提高幅度最大,該統(tǒng)計(jì)結(jié)果可以應(yīng)用于業(yè)務(wù)逐6分鐘混合仰角反射率產(chǎn)品檢驗(yàn)中,作為參考樣本實(shí)時評估產(chǎn)品改進(jìn)效果。

    表2 2019年5—10月陜西各雷達(dá)平均RMSE、MAE和CC Table 2 Average RMSE, MAE and CC of radar in Shaanxi from May to October 2019

    4 結(jié)論

    本文通過分析陜西多普勒天氣雷達(dá)因人造建筑物、森林等非地形因素造成的降水觀測空間不連續(xù)問題,基于精細(xì)的地形數(shù)據(jù)信息,提出了一種雷達(dá)電磁波非常規(guī)遮擋回波補(bǔ)償技術(shù),結(jié)合長時段雷達(dá)估測降水累積分布和地面實(shí)況雨量分布制定了各雷達(dá)遮擋修訂方案,并利用降水典型個例和汛期所有降水過程檢驗(yàn)評估了遮擋修訂方案對用于雷達(dá)定量降水估測的提升效果,最終得到以下結(jié)論:

    1)陜西復(fù)雜地形環(huán)境導(dǎo)致雷達(dá)觀測資料需要進(jìn)行最優(yōu)無縫混合仰角處理遮擋現(xiàn)象后才能正確計(jì)算用于降水估測的降水回波。首先利用SRTM地形數(shù)據(jù)和雷達(dá)體掃信息計(jì)算得出各雷達(dá)最低有效混合仰角參數(shù)信息;其次利用長時間段雷達(dá)QPE累積分布識別非地形遮擋(NSB)導(dǎo)致的降水空間不連續(xù)區(qū)域,經(jīng)過交叉方位角插值、更高一層仰角數(shù)據(jù)填補(bǔ)和線性加權(quán)滑動平均等方法,最終形成最優(yōu)無縫非常規(guī)電磁波遮擋補(bǔ)償技術(shù)。

    2)修訂前、后的雷達(dá)混合仰角反射率及QPE分布對比結(jié)果表明,經(jīng)過雷達(dá)電磁波非常規(guī)遮擋補(bǔ)償技術(shù)修訂后可以有效填補(bǔ)反射率分布和QPE分布中大范圍或局地楔形的遮擋不連續(xù)區(qū)域,較為準(zhǔn)確地還原雷達(dá)實(shí)際降水回波分布。經(jīng)過最優(yōu)無縫混合掃描仰角處理后的降水回波覆蓋面和平滑程度較改進(jìn)前有明顯提高,采用的遮擋不連續(xù)修訂參考方案可直接應(yīng)用到業(yè)務(wù)觀測中,實(shí)時自動修訂得出混合仰角反射率產(chǎn)品。

    3)典型個例和汛期降水過程檢驗(yàn)評估顯示在固定

    Z

    -

    R

    關(guān)系下,經(jīng)過遮擋不連續(xù)修訂后的雷達(dá)估測降水量與地面降水量之間均方根誤差、平均絕對誤差明顯減小,兩者具有更好的線性相關(guān)性,QPE分布比修訂前更接近于地面降水分布,并且在強(qiáng)對流降水過程中能有效改善由于遮擋導(dǎo)致強(qiáng)降水回波缺失現(xiàn)象。

    由于部分降水過程中存在低層仰角電磁波束部分被遮擋現(xiàn)象,造成低層仰角反射率強(qiáng)度偏小但不能忽視,或者存在降水系統(tǒng)相對淺薄、云頂高度較低導(dǎo)致抬高仰角后電磁波束超過云頂?shù)惹闆r,并且NSB修訂參考表格隨城市發(fā)展需要及時更新,因此后續(xù)仍需繼續(xù)開展更多研究以便進(jìn)一步提高改進(jìn)陜西雷達(dá)遮擋修訂工作及QPE精度。

    猜你喜歡
    仰角方位角電磁波
    聚焦電磁波和相對論簡介
    電磁波和相對論簡介考點(diǎn)解讀
    探究無線電方位在無線電領(lǐng)航教學(xué)中的作用和意義
    卷宗(2021年2期)2021-03-09 07:57:24
    近地磁尾方位角流期間的場向電流增強(qiáng)
    用銳角三角函數(shù)解決仰角、俯角問題
    用有源音箱驗(yàn)證電磁波的發(fā)射和接收
    分段三次Hermite插值計(jì)算GNSS系統(tǒng)衛(wèi)星仰角
    科技視界(2016年15期)2016-06-30 12:48:10
    向量內(nèi)外積在直線坐標(biāo)方位角反算中的應(yīng)用研究
    河南科技(2015年18期)2015-11-25 08:50:14
    平行透刺聯(lián)合電磁波治療肩周炎32例
    基于對話框編寫實(shí)現(xiàn)坐標(biāo)方位角正算與反算的程序
    城市勘測(2013年5期)2013-03-06 05:44:13
    国产亚洲欧美在线一区二区| 亚洲男人的天堂狠狠| 免费一级毛片在线播放高清视频| 黄色成人免费大全| 99re在线观看精品视频| 精品国产乱码久久久久久男人| 精品熟女少妇八av免费久了| 亚洲av中文字字幕乱码综合 | 丝袜人妻中文字幕| 校园春色视频在线观看| 亚洲自偷自拍图片 自拍| 大型黄色视频在线免费观看| 国产黄a三级三级三级人| 桃红色精品国产亚洲av| 国产精品永久免费网站| www日本在线高清视频| 精品人妻1区二区| 亚洲av美国av| 一本一本综合久久| 人人澡人人妻人| 亚洲熟女毛片儿| 久久精品国产亚洲av香蕉五月| 午夜亚洲福利在线播放| 成熟少妇高潮喷水视频| 欧美乱码精品一区二区三区| 久久久久亚洲av毛片大全| 哪里可以看免费的av片| 一进一出抽搐动态| 亚洲精品av麻豆狂野| 在线观看www视频免费| 丰满的人妻完整版| 91国产中文字幕| 国产av一区在线观看免费| 亚洲国产精品合色在线| 久久99热这里只有精品18| 满18在线观看网站| 久久伊人香网站| 曰老女人黄片| 久久精品国产清高在天天线| 视频在线观看一区二区三区| 精品久久久久久久久久免费视频| 老司机靠b影院| 亚洲精品色激情综合| 男人的好看免费观看在线视频 | 国产精品99久久99久久久不卡| 丁香欧美五月| 日本 欧美在线| 亚洲国产中文字幕在线视频| 精品少妇一区二区三区视频日本电影| 亚洲第一青青草原| 精品熟女少妇八av免费久了| 亚洲一区二区三区不卡视频| 午夜福利欧美成人| 日韩高清综合在线| 在线播放国产精品三级| 欧美三级亚洲精品| 99热只有精品国产| 日韩有码中文字幕| av电影中文网址| 男女午夜视频在线观看| 91av网站免费观看| 午夜福利成人在线免费观看| 国产黄片美女视频| 久久久久久久久久黄片| 在线观看日韩欧美| 美国免费a级毛片| 亚洲成人精品中文字幕电影| 69av精品久久久久久| 精品午夜福利视频在线观看一区| 国产精品1区2区在线观看.| 亚洲精品久久成人aⅴ小说| 99国产精品一区二区蜜桃av| 好看av亚洲va欧美ⅴa在| 欧美又色又爽又黄视频| 免费看美女性在线毛片视频| 国产v大片淫在线免费观看| 搞女人的毛片| 黄色女人牲交| 一区二区三区精品91| 制服人妻中文乱码| 国产精品综合久久久久久久免费| 久久中文字幕人妻熟女| 人妻久久中文字幕网| 色老头精品视频在线观看| 精品国产一区二区三区四区第35| 男人的好看免费观看在线视频 | 亚洲精品国产区一区二| 露出奶头的视频| 非洲黑人性xxxx精品又粗又长| 免费一级毛片在线播放高清视频| 此物有八面人人有两片| 大型av网站在线播放| 99久久精品国产亚洲精品| 中文资源天堂在线| 婷婷六月久久综合丁香| 麻豆国产av国片精品| 99精品在免费线老司机午夜| 97碰自拍视频| 亚洲,欧美精品.| 18禁美女被吸乳视频| 又大又爽又粗| 国产高清激情床上av| 91麻豆精品激情在线观看国产| 久久久国产精品麻豆| 亚洲av日韩精品久久久久久密| 成年人黄色毛片网站| 国产激情久久老熟女| 欧美绝顶高潮抽搐喷水| 欧美+亚洲+日韩+国产| 在线永久观看黄色视频| 国产精品爽爽va在线观看网站 | 美女 人体艺术 gogo| 亚洲avbb在线观看| 男男h啪啪无遮挡| 国产一级毛片七仙女欲春2 | 一卡2卡三卡四卡精品乱码亚洲| 在线国产一区二区在线| 熟妇人妻久久中文字幕3abv| 亚洲成人久久性| 久久精品夜夜夜夜夜久久蜜豆 | 美女午夜性视频免费| 亚洲精品色激情综合| 久久久久久大精品| 亚洲精品中文字幕一二三四区| 一级作爱视频免费观看| www.精华液| 国产精品亚洲美女久久久| 一级黄色大片毛片| 老熟妇乱子伦视频在线观看| 精品少妇一区二区三区视频日本电影| 男人舔奶头视频| 国产主播在线观看一区二区| 久久99热这里只有精品18| 亚洲成av片中文字幕在线观看| 国产精品av久久久久免费| 免费搜索国产男女视频| 男人的好看免费观看在线视频 | 免费观看人在逋| cao死你这个sao货| 久久这里只有精品19| 欧美另类亚洲清纯唯美| 精品久久久久久久人妻蜜臀av| 欧美绝顶高潮抽搐喷水| 身体一侧抽搐| 亚洲专区中文字幕在线| 亚洲精品久久成人aⅴ小说| 欧美激情极品国产一区二区三区| 亚洲国产欧美一区二区综合| bbb黄色大片| 正在播放国产对白刺激| 亚洲中文日韩欧美视频| 中文亚洲av片在线观看爽| 亚洲精品一区av在线观看| 男女床上黄色一级片免费看| 叶爱在线成人免费视频播放| 国产色视频综合| 日韩一卡2卡3卡4卡2021年| 日本三级黄在线观看| 怎么达到女性高潮| 婷婷亚洲欧美| 亚洲真实伦在线观看| 久久精品国产亚洲av香蕉五月| 欧美日韩乱码在线| a在线观看视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美精品综合久久99| 久久国产精品人妻蜜桃| 久久久久久久久中文| 又紧又爽又黄一区二区| 黄网站色视频无遮挡免费观看| 国产精品久久久久久亚洲av鲁大| 国产高清激情床上av| 欧美黑人欧美精品刺激| 视频在线观看一区二区三区| 夜夜夜夜夜久久久久| 男男h啪啪无遮挡| 黑人巨大精品欧美一区二区mp4| 国产亚洲精品第一综合不卡| 亚洲熟妇熟女久久| 免费观看人在逋| 成人手机av| 色综合欧美亚洲国产小说| 日日摸夜夜添夜夜添小说| 妹子高潮喷水视频| 搡老岳熟女国产| 美女高潮喷水抽搐中文字幕| 非洲黑人性xxxx精品又粗又长| 久久久国产成人精品二区| 国产在线精品亚洲第一网站| 亚洲欧美精品综合一区二区三区| 国产主播在线观看一区二区| 狠狠狠狠99中文字幕| 免费在线观看黄色视频的| 中文字幕av电影在线播放| 久久精品91蜜桃| 国产精品久久久久久精品电影 | 欧美绝顶高潮抽搐喷水| 黄色 视频免费看| 欧美日本亚洲视频在线播放| 一二三四社区在线视频社区8| 老鸭窝网址在线观看| 夜夜夜夜夜久久久久| 午夜免费鲁丝| 久久国产精品影院| 男女午夜视频在线观看| 免费高清在线观看日韩| 国产真实乱freesex| 精品国内亚洲2022精品成人| 手机成人av网站| 国产三级在线视频| 哪里可以看免费的av片| 亚洲一码二码三码区别大吗| 精品欧美一区二区三区在线| 国产伦在线观看视频一区| 动漫黄色视频在线观看| 99国产极品粉嫩在线观看| 哪里可以看免费的av片| av福利片在线| 男女那种视频在线观看| 少妇 在线观看| 黄色女人牲交| tocl精华| 精品久久久久久久久久免费视频| 一级黄色大片毛片| 欧美激情 高清一区二区三区| 麻豆av在线久日| 亚洲国产精品久久男人天堂| 熟妇人妻久久中文字幕3abv| 日韩欧美 国产精品| 亚洲成av片中文字幕在线观看| 午夜老司机福利片| 50天的宝宝边吃奶边哭怎么回事| 欧美成人午夜精品| 国产欧美日韩精品亚洲av| 国产精品美女特级片免费视频播放器 | 在线观看舔阴道视频| 亚洲成人久久性| 国产一卡二卡三卡精品| 两个人视频免费观看高清| 国产欧美日韩一区二区精品| 看片在线看免费视频| 神马国产精品三级电影在线观看 | 亚洲av成人av| 亚洲国产精品成人综合色| 精品无人区乱码1区二区| 国产高清videossex| 久久国产亚洲av麻豆专区| 国内少妇人妻偷人精品xxx网站 | 中文字幕久久专区| 12—13女人毛片做爰片一| 久久久久久久久久黄片| 亚洲一码二码三码区别大吗| 国产成人一区二区三区免费视频网站| 91老司机精品| 亚洲成人久久性| 天堂动漫精品| 美女 人体艺术 gogo| 热re99久久国产66热| 少妇粗大呻吟视频| 美女午夜性视频免费| 欧美日韩乱码在线| 亚洲欧洲精品一区二区精品久久久| 老熟妇乱子伦视频在线观看| 久久狼人影院| 久久精品国产亚洲av高清一级| 免费看日本二区| 亚洲专区字幕在线| 一级a爱片免费观看的视频| 亚洲国产精品成人综合色| www.精华液| 欧美日韩瑟瑟在线播放| 黄色女人牲交| 亚洲五月天丁香| 久久伊人香网站| 亚洲成人免费电影在线观看| 久久精品影院6| 国产极品粉嫩免费观看在线| 成人精品一区二区免费| 神马国产精品三级电影在线观看 | 长腿黑丝高跟| 色精品久久人妻99蜜桃| 国产成人啪精品午夜网站| 日韩精品中文字幕看吧| 久久婷婷人人爽人人干人人爱| 国产亚洲av嫩草精品影院| 精品久久久久久久久久久久久 | 国产真人三级小视频在线观看| 18美女黄网站色大片免费观看| 久久午夜综合久久蜜桃| 在线观看www视频免费| www.999成人在线观看| 欧美日韩福利视频一区二区| 色播在线永久视频| 亚洲国产精品sss在线观看| 日本 欧美在线| 精品国内亚洲2022精品成人| 日韩欧美一区二区三区在线观看| www国产在线视频色| 国产爱豆传媒在线观看 | 制服人妻中文乱码| 精品国产美女av久久久久小说| 精品久久蜜臀av无| 搡老岳熟女国产| 99久久精品国产亚洲精品| 一进一出抽搐gif免费好疼| 亚洲五月婷婷丁香| 国产视频一区二区在线看| 欧美成人性av电影在线观看| 少妇粗大呻吟视频| 中亚洲国语对白在线视频| 搡老熟女国产l中国老女人| 亚洲精品国产精品久久久不卡| 校园春色视频在线观看| 男女下面进入的视频免费午夜 | 嫩草影视91久久| 一本一本综合久久| 欧美亚洲日本最大视频资源| www.www免费av| 欧美日韩精品网址| 国产成人精品久久二区二区91| 国产熟女xx| 成人亚洲精品一区在线观看| 热re99久久国产66热| 日韩三级视频一区二区三区| 免费在线观看日本一区| 国产私拍福利视频在线观看| 校园春色视频在线观看| 久久亚洲精品不卡| 午夜久久久在线观看| 久久国产精品男人的天堂亚洲| 十分钟在线观看高清视频www| 伊人久久大香线蕉亚洲五| 国产熟女午夜一区二区三区| 999精品在线视频| 大香蕉久久成人网| 在线观看免费视频日本深夜| 精品久久久久久久末码| 国产亚洲av高清不卡| 久久天躁狠狠躁夜夜2o2o| 激情在线观看视频在线高清| 亚洲国产精品成人综合色| 女性被躁到高潮视频| 搡老妇女老女人老熟妇| 亚洲va日本ⅴa欧美va伊人久久| 成在线人永久免费视频| 亚洲国产日韩欧美精品在线观看 | 色老头精品视频在线观看| 男女之事视频高清在线观看| 欧美日韩瑟瑟在线播放| 国产精品久久久久久人妻精品电影| 不卡av一区二区三区| 叶爱在线成人免费视频播放| 国产免费男女视频| 18禁美女被吸乳视频| 熟妇人妻久久中文字幕3abv| 满18在线观看网站| 国产精品99久久99久久久不卡| 免费搜索国产男女视频| 一级黄色大片毛片| 午夜福利视频1000在线观看| 老汉色∧v一级毛片| 一本一本综合久久| 日韩免费av在线播放| 亚洲成人久久爱视频| 老司机靠b影院| 国产欧美日韩一区二区精品| 亚洲专区中文字幕在线| 日韩一卡2卡3卡4卡2021年| 午夜福利在线在线| 国产97色在线日韩免费| 亚洲第一av免费看| √禁漫天堂资源中文www| 男人舔女人下体高潮全视频| 亚洲国产精品久久男人天堂| 久久中文看片网| 九色国产91popny在线| 欧美激情 高清一区二区三区| 制服人妻中文乱码| 亚洲av成人av| 亚洲中文日韩欧美视频| 中文字幕另类日韩欧美亚洲嫩草| 无限看片的www在线观看| 黄频高清免费视频| 国产乱人伦免费视频| 色播亚洲综合网| av免费在线观看网站| 在线观看午夜福利视频| 国产一区二区三区视频了| 1024手机看黄色片| 男人舔女人的私密视频| 国内精品久久久久久久电影| 91成人精品电影| 狂野欧美激情性xxxx| 欧美在线黄色| 丰满人妻熟妇乱又伦精品不卡| 他把我摸到了高潮在线观看| 正在播放国产对白刺激| 亚洲av成人av| 亚洲av成人一区二区三| 18禁美女被吸乳视频| 美女国产高潮福利片在线看| 中国美女看黄片| 女同久久另类99精品国产91| 欧美乱妇无乱码| 国产高清视频在线播放一区| xxxwww97欧美| 中文字幕人成人乱码亚洲影| 正在播放国产对白刺激| 亚洲精品美女久久av网站| 成人精品一区二区免费| 免费女性裸体啪啪无遮挡网站| 久久久久国产精品人妻aⅴ院| 日韩av在线大香蕉| 国产精品99久久99久久久不卡| 在线观看午夜福利视频| 黄色 视频免费看| 黄色片一级片一级黄色片| 亚洲av成人不卡在线观看播放网| 欧美中文日本在线观看视频| 婷婷六月久久综合丁香| 午夜免费鲁丝| 免费高清在线观看日韩| 亚洲一区中文字幕在线| 男人的好看免费观看在线视频 | 欧美激情极品国产一区二区三区| 国产熟女xx| 国产高清有码在线观看视频 | 日韩成人在线观看一区二区三区| 久久草成人影院| 天天躁狠狠躁夜夜躁狠狠躁| 丰满人妻熟妇乱又伦精品不卡| 黄片小视频在线播放| 一级毛片女人18水好多| 国产精品精品国产色婷婷| 色综合婷婷激情| 黄频高清免费视频| 精品少妇一区二区三区视频日本电影| 久9热在线精品视频| www.自偷自拍.com| 国产v大片淫在线免费观看| 男女做爰动态图高潮gif福利片| 观看免费一级毛片| 久久久久久久午夜电影| 天天一区二区日本电影三级| 1024视频免费在线观看| 久久久久久人人人人人| 狂野欧美激情性xxxx| 草草在线视频免费看| 级片在线观看| 久久久久久久午夜电影| www日本在线高清视频| 免费看十八禁软件| 日本免费一区二区三区高清不卡| www国产在线视频色| 午夜福利一区二区在线看| xxx96com| or卡值多少钱| 女性被躁到高潮视频| 日韩大码丰满熟妇| 国产爱豆传媒在线观看 | 熟妇人妻久久中文字幕3abv| 久久精品成人免费网站| 国产爱豆传媒在线观看 | 欧美性长视频在线观看| 高清毛片免费观看视频网站| 日韩精品青青久久久久久| 欧美中文日本在线观看视频| 免费av毛片视频| 国产精品电影一区二区三区| 国产黄色小视频在线观看| 一二三四在线观看免费中文在| 国产亚洲av嫩草精品影院| 亚洲av成人不卡在线观看播放网| 国产亚洲欧美精品永久| 欧美黑人欧美精品刺激| 好看av亚洲va欧美ⅴa在| av超薄肉色丝袜交足视频| 999久久久精品免费观看国产| 国产麻豆成人av免费视频| 亚洲第一av免费看| 亚洲人成网站高清观看| 嫩草影院精品99| 白带黄色成豆腐渣| 后天国语完整版免费观看| 亚洲真实伦在线观看| 日韩欧美在线二视频| 天堂√8在线中文| 成人永久免费在线观看视频| 中国美女看黄片| 18禁观看日本| 欧美另类亚洲清纯唯美| 亚洲精品美女久久av网站| 久久精品亚洲精品国产色婷小说| 99riav亚洲国产免费| 97碰自拍视频| √禁漫天堂资源中文www| 免费高清视频大片| 免费在线观看视频国产中文字幕亚洲| 黄片大片在线免费观看| 最好的美女福利视频网| 午夜久久久久精精品| 波多野结衣高清无吗| 亚洲国产中文字幕在线视频| 国产一区二区在线av高清观看| 欧美性长视频在线观看| 久久午夜亚洲精品久久| 亚洲国产精品sss在线观看| 最新美女视频免费是黄的| 亚洲精品美女久久久久99蜜臀| 啦啦啦观看免费观看视频高清| 精品国产乱码久久久久久男人| 精品不卡国产一区二区三区| 亚洲人成伊人成综合网2020| 国产免费av片在线观看野外av| 国产激情偷乱视频一区二区| 欧美成人性av电影在线观看| 免费在线观看影片大全网站| 色综合欧美亚洲国产小说| 亚洲熟妇中文字幕五十中出| 亚洲一区二区三区不卡视频| 桃红色精品国产亚洲av| 亚洲美女黄片视频| 久久久久久久久免费视频了| 国产伦人伦偷精品视频| 久久久久久免费高清国产稀缺| 亚洲免费av在线视频| 亚洲片人在线观看| 久久国产亚洲av麻豆专区| 色播亚洲综合网| 欧美一区二区精品小视频在线| 久久久国产成人精品二区| 欧美黄色片欧美黄色片| 搡老熟女国产l中国老女人| 18禁裸乳无遮挡免费网站照片 | 欧美乱妇无乱码| 欧美av亚洲av综合av国产av| 国产在线精品亚洲第一网站| 性色av乱码一区二区三区2| 精品一区二区三区视频在线观看免费| 性欧美人与动物交配| 国产真人三级小视频在线观看| 精品一区二区三区av网在线观看| 国产伦人伦偷精品视频| 久久婷婷成人综合色麻豆| www国产在线视频色| 波多野结衣高清作品| 在线十欧美十亚洲十日本专区| 欧美色视频一区免费| 巨乳人妻的诱惑在线观看| 高潮久久久久久久久久久不卡| 免费女性裸体啪啪无遮挡网站| 999久久久精品免费观看国产| 老司机午夜福利在线观看视频| 久久精品影院6| 亚洲午夜精品一区,二区,三区| 日韩有码中文字幕| 国产蜜桃级精品一区二区三区| 亚洲av成人一区二区三| 桃红色精品国产亚洲av| 精品久久久久久,| 国内精品久久久久精免费| 成人午夜高清在线视频 | www国产在线视频色| 亚洲av成人av| 老司机在亚洲福利影院| 久9热在线精品视频| 久久精品夜夜夜夜夜久久蜜豆 | 白带黄色成豆腐渣| 99热6这里只有精品| 国产精品,欧美在线| 久久精品亚洲精品国产色婷小说| 真人一进一出gif抽搐免费| 中文字幕人成人乱码亚洲影| 亚洲av中文字字幕乱码综合 | 久久精品人妻少妇| 午夜久久久久精精品| 免费观看人在逋| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美在线二视频| 国产精品自产拍在线观看55亚洲| 一边摸一边做爽爽视频免费| 狂野欧美激情性xxxx| 19禁男女啪啪无遮挡网站| 一级a爱视频在线免费观看| 免费在线观看亚洲国产| 深夜精品福利| 色婷婷久久久亚洲欧美| 国产成人一区二区三区免费视频网站| 亚洲激情在线av| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美色视频一区免费| 久久亚洲精品不卡| 欧美乱码精品一区二区三区| 精品欧美国产一区二区三| 亚洲精品久久成人aⅴ小说| 黑人操中国人逼视频| 久久九九热精品免费| 国内揄拍国产精品人妻在线 | 在线天堂中文资源库| 一二三四在线观看免费中文在| 久久久久久亚洲精品国产蜜桃av| 亚洲精品美女久久久久99蜜臀| 国产精品一区二区精品视频观看| 成年免费大片在线观看| АⅤ资源中文在线天堂| 一级黄色大片毛片| 狂野欧美激情性xxxx| 亚洲av日韩精品久久久久久密| 亚洲国产高清在线一区二区三 | 午夜免费成人在线视频| 国产在线精品亚洲第一网站| 色播亚洲综合网| 国产亚洲精品第一综合不卡| 亚洲 欧美一区二区三区| www.www免费av|