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

    機載GNSS反射信號海面測高模型的研究

    2020-05-22 09:50:00張云張楊陽孟婉婷楊樹瑚韓彥嶺
    海洋學(xué)報 2020年3期
    關(guān)鍵詞:仰角海面反演

    張云,張楊陽,孟婉婷,楊樹瑚,韓彥嶺

    (1.上海海洋大學(xué) 信息學(xué)院,上海 201306;2.上海航天電子技術(shù)研究所,上海 201109)

    1 引言

    海面高度是物理海洋、地球物理學(xué)、大地測量學(xué)等領(lǐng)域的基礎(chǔ)數(shù)據(jù)。全球海面高度數(shù)據(jù)可以用于大洋環(huán)流的確定、大洋潮汐模型的建立以及中尺度氣候模型的研究等,因此獲取高精度海面高度對氣象學(xué)和海洋學(xué)都非常重要[1]。GNSS-R(Global Navigation Satellite Systems Reflectometry)技術(shù)是目前海面測高中的新興遙感技術(shù),屬于被動接收技術(shù),通過接收經(jīng)由反射面反射的反射信號與直射信號之間的差異進(jìn)行物理量的反演。該技術(shù)在1993年,由Martin-Neira教授首次提出[2],并迅速在國內(nèi)外成為研究熱點。GNSS-R技術(shù)有眾多優(yōu)點,例如擁有大量信號源、不受天氣(云、霧)影響、低成本實時監(jiān)測等[3-4]。

    機載GNSS-R較岸基GNSS-R相比,其優(yōu)勢在于其空間分辨率高、監(jiān)測范圍廣,可對大范圍的區(qū)域進(jìn)行監(jiān)測,同時機載測高已經(jīng)成為了星載測高技術(shù)的預(yù)研技術(shù),因此GNSS-R機載測高正在被得到廣泛的研究。目前GNSS-R機載測高技術(shù)主要包括碼測高[5-7]、載波相位測高[8-13]、干涉測高[14]、時延多普勒圖(Delay-Doppler-Map,DDM)測高[15-16]及信噪比測高[17-18]等。國外團(tuán)隊對機載GNSS-R測高進(jìn)行了相關(guān)研究,取得了一定的研究成果[19-20],而國內(nèi)大部分測高研究大部分集中基于岸基開展的[17,21],關(guān)于機載GNSS-R海面測高方面的研究較少,缺少必要的實驗成果及分析。

    碼測高技術(shù)是目前最常用的GNSS-R機載測高技術(shù),相對于其他技術(shù),具有實現(xiàn)方便,抗干擾能力強的特點,因此在GNSS-R移動平臺的海面監(jiān)測領(lǐng)域有著廣泛的應(yīng)用。本文在GNSS-R岸基碼測高技術(shù)研究[22]的基礎(chǔ)上,分析了2011年11月11日,CSICIEEC(Consejo Superior de Investigaciones Científicas(CSIC)—Institut d'Estudis Espacials de Catalunya(IEEC))在芬蘭波羅的海進(jìn)行的機載實驗數(shù)據(jù),在岸基反演的基礎(chǔ)上,針對大氣延遲、天線距離、海洋潮汐模型等進(jìn)行修正,優(yōu)化了GNSS-R機載海面測高模型,成功地反演了海面高度,證明了機載GNSS-R技術(shù)海面測高的可行性,在分析過程中,關(guān)于反射仰角的測高精度的影響做了定性的分析,得出了低反射仰角對測高結(jié)果精度產(chǎn)生較大影響的結(jié)論。

    2 機載GNSS-R海面高度反演原理

    2.1 GNSS-R碼測高技術(shù)

    碼測高技術(shù)對直射信號的本地碼進(jìn)行適當(dāng)延遲和多普勒頻移補償后,再生成反射信號[23],它通過反射信號與直射信號之間的路徑延遲進(jìn)行高度反演。GNSS-R碼測高技術(shù)基本概念如圖1所示,其中延遲距離為Dp。

    圖1 GNSS-R碼測高技術(shù)基本概念Fig.1 GNSS-R group delay altimetry basic concept

    2.2 延遲路徑模型

    圖1中的反射信號相對于直射信號的延遲路徑ρ由下式給出,

    式中,Dwin為直射信號開始窗口至反射信號開始窗口距離;Ddir為直射信號峰值位置;Dref為反射信號導(dǎo)數(shù)最大位置;Datm為大氣延遲;Dins為反射天線和直射天線之間的距離。

    通過對機載實驗中使用反射信號峰值位置及導(dǎo)數(shù)最大位置反演出的結(jié)果進(jìn)行對比,結(jié)果表明,在機載實驗中,采用反射信號導(dǎo)數(shù)最大值位置時的誤差比采用直射信號峰值位置時的誤差小[19]。因此,本文中選用反射信號導(dǎo)數(shù)最大位置作為反射信號的延遲位置。

    在機載測高模型中,將考慮到大氣延遲的干擾。由于本次實驗中飛機高度處于大約3 km的高度,所以不對電離層進(jìn)行考慮,只考慮對流層延遲進(jìn)行誤差修正。本文中,大氣層的延遲Datm使用下式計算[24],

    式中,e為鏡面反射點處仰角;Ha為飛機上方直射接收天線所在的高度。

    與岸基GNSS-R反演測高不同的是,由于直射接收天線與反射接收天線分別位于飛機的上方與下方,直射天線與反射天線位置不在同一個水平面上,所以在進(jìn)行偽距測量時,對兩個天線之間的距離(Dins)進(jìn)行誤差修正。由于缺少機載平臺的姿態(tài)角信息,本文只選取飛機直線飛行時的數(shù)據(jù)進(jìn)行分析,刪除了飛機轉(zhuǎn)向時的實驗數(shù)據(jù),在計算過程中,忽略了姿態(tài)角對天線距離的影響。

    2.3 機載測高模型

    結(jié)合延遲路徑模型,機載測高模型如圖2所示。

    圖2 機載測高模型Fig.2 Airborne altimetry model

    本文在岸基GNSS-R測高模型[22]的基礎(chǔ)上,優(yōu)化了機載GNSS-R測高模型,根據(jù)機載特點,增加了對大氣延遲以及天線之間距離的修正,如圖2中的Datm及Dins所示。同時,增加了對低仰角下反演結(jié)果的閾值判斷。

    根據(jù)公式(1)計算出延遲路徑之后,通過下式得出GNSS-R反演的直射天線到海面的高度H-R,

    通過接收機接收到的位置信息,可以得到飛機上方直射接收天線的高度(Ha)。通過計算飛機上方直射接收天線的高度(Ha)與GNSS-R反演的直射天線到海面的高度(H-R)之間的距離差,以及通過潮汐模型HOTM的修正,可以得到海面高度(Hsea),

    3 實驗介紹

    3.1 實驗環(huán)境介紹

    本次實驗數(shù)據(jù)為2011年11月11日,CSIC-IEEC在芬蘭波羅的海進(jìn)行的機載實驗數(shù)據(jù)。本次機載實驗采用接收機GOLD-RTR(GPS Open Loop Differential Real-Time Receiver)進(jìn)行數(shù)據(jù)收集,它的采樣率為1 ms。其中,直射信號接收天線位于飛機上方,反射信號接收天線位于飛機下方,天線增益均為15 -dBi。其中在飛機上方,安裝了Novatel接收機用于提供時間、位置等信息。

    3.2 實驗數(shù)據(jù)介紹

    本文采用了兩個時間段的數(shù)據(jù)進(jìn)行了分析,同時選取鏡面反射點仰角大于75°的時間段作為高仰角,小于60°的時間段作為低仰角進(jìn)行分析。

    第一個實驗數(shù)據(jù)時間段全球定位系統(tǒng)(Global Positioning System,GPS)時間為 452 008~453 350 s,飛機的軌跡圖如圖3a所示。高仰角數(shù)據(jù)選取偽隨機噪聲碼(Pseudo Random Noise code,PRN)12 號衛(wèi)星,鏡面反射點仰角為80.51°~84.29°;低仰角數(shù)據(jù)選取PRN25號衛(wèi)星,鏡面反射點仰角為46.60°~56.78°。為了避免飛機轉(zhuǎn)向所造成的影響,去除了GPS時間452 608~452 750 s的轉(zhuǎn)向時間,只選取飛機沿直線飛行時的數(shù)據(jù)作為本次實驗分析數(shù)據(jù),此時間段仰角變化如圖4a所示。

    第二個實驗數(shù)據(jù)時間段GPS時間為456 700~458 099 s,飛機的軌跡圖如圖3b所示。高仰角數(shù)據(jù)選取PRN25號衛(wèi)星,鏡面反射點仰角為76.82°~79.93°;低仰角數(shù)據(jù)選取PRN12號衛(wèi)星,鏡面反射點仰角為47.30°~58.30°。為了避免飛機轉(zhuǎn)向所造成的影響,去除了GPS時間457 300~457 499 s的轉(zhuǎn)向時間,選取飛機沿直線飛行時的數(shù)據(jù)作為本次實驗分析數(shù)據(jù),此時間段仰角變化如圖4b所示。

    兩個實驗數(shù)據(jù)總結(jié)如表1所示。

    3.3 模型驗證數(shù)據(jù)介紹

    在進(jìn)行機載測高分析時,GNSS-R獲取的是瞬時海面高度,因此需要消除波浪影響,進(jìn)行潮位改正,獲得海面高度(Hsea)。本文選取丹麥技術(shù)大學(xué)(Technical University of Denmark,DTU10)全球海洋潮汐模型(Global Ocean Tide Model)作為本次實驗的潮汐模型HOTM[25]。

    同時,本文選取DTU10海洋平均海面高度[26-28]作為本次實驗的參考值HDTU。公式(4)求得的結(jié)果(Hsea)與DTU10海洋平均海面高度HDTU進(jìn)行對比,可以驗證GNSS-R反演的海面高度結(jié)果(Hsea)的精度。

    4 反演結(jié)果分析

    4.1 鏡面反射點

    圖3 飛行軌跡圖Fig.3 Flight trajectory

    圖4 PRN12號衛(wèi)星(紅)及PRN25號衛(wèi)星(藍(lán))仰角變化Fig.4 Variation of PRN12 satellite (red) and PRN25 satellite (blue) elevation angle

    表1 實驗數(shù)據(jù)結(jié)果Table 1 Result of experimental data

    第一段實驗時間期間,PRN12和PRN25的鏡面反射點軌跡如圖5a所示,第二段實驗時間期間,PRN12和PRN25鏡面反射點軌跡如圖5b所示。從圖5可知,所有的反射信號數(shù)據(jù)都落在海面范圍內(nèi)。

    4.2 海面高度反演

    以GPS時間452 008 s,PRN12的數(shù)據(jù)分析為例。首先,獲取直射信號與反射信號原始I/Q波形之后,對1 s的數(shù)據(jù)長度進(jìn)行篩選。本次實驗中數(shù)據(jù)采樣頻率為1 kHz,基于大量數(shù)據(jù)分析所得到的結(jié)論,為了考慮數(shù)據(jù)的完整性,選取1 s數(shù)據(jù)包含大于800個原始數(shù)據(jù)進(jìn)行分析,這個范圍可以剔除掉由于波形異常從而導(dǎo)致結(jié)果不準(zhǔn)確的情況。在選取數(shù)據(jù)之后,對它們進(jìn)行1 ms相干累加及1 s非相干累加操作。至此,得出1 s的直射/反射信號功率波形,接著對直射/反射信號波形分別計算直射信號峰值位置Ddir及反射信號導(dǎo)數(shù)最大位置Dref。此時通過1 ms相干累加及1 s非相干累加操作以后的直射信號/反射信號圖如圖1中的直射信號及反射信號所示。

    通過了誤差模型,以及已知的設(shè)置,利用公式(4),得到了該時刻通過反射信號反演的海面高度是18.04 m,該時刻的海面高度值是16.79 m。

    4.3 第一段實驗結(jié)果

    第一段實驗時間期間,在高仰角情況下,即PRN12號衛(wèi)星作為觀測衛(wèi)星,通過反射信號反演的海面高度如圖6a所示。在低仰角情況下,即PRN25號衛(wèi)星作為觀測衛(wèi)星,通過反射信號反演的海面高度如圖6b所示。得出1 s結(jié)果之后,對結(jié)果進(jìn)行20 s平均處理。進(jìn)行平均操作之后的結(jié)果如圖7所示,圖7a是20 s平均后(PRN12號衛(wèi)星)反演的高度圖,圖7b是20 s平均后(PRN25號衛(wèi)星)反演的高度圖。

    4.4 第二段實驗結(jié)果

    第二段實驗時間期間,發(fā)現(xiàn)在低仰角的情況下,GPS時間段456 700~457 299 s(PRN12號衛(wèi)星)期間,反演的海面高度出現(xiàn)幾百米誤差的異常值,異常反演海面高度如圖8a所示。經(jīng)過分析,發(fā)現(xiàn)在這個時間段PRN12號衛(wèi)星接收到的直射信號異常,導(dǎo)致所得的波形無法計算出準(zhǔn)確的結(jié)果。

    圖9a是直射信號正常功率波形,圖9b是GPS時刻456 793 s的直射信號功率波形,通過比較圖9a與圖9b的波形,發(fā)現(xiàn)在圖9b時刻的直射功率波形異常,推測在低仰角時,直射天線可能受到了機身的干擾,導(dǎo)致直射信號異常。實驗采用的接收機的直射信號采樣總共有64個相關(guān)器,采樣步長為15 m,通常情況下直射信號的相關(guān)峰值位置在第34個相關(guān)器正負(fù)1個步長距離間波動,直射信號的峰值位置范圍(Ddir)為480~510 m。所以,本文在低仰角數(shù)據(jù)反演前,對直射信號的峰值位置范圍進(jìn)行閾值判斷(480<Ddir<510),不滿足條件的數(shù)據(jù)將不參與高度反演計算,經(jīng)過篩選之后的海面高度如圖8b所示,與圖8a比較后可以發(fā)現(xiàn),經(jīng)過篩選后,異常值得到了有效的刪除,反演的可靠性得到了提高。

    在高仰角情況下,即PRN25號衛(wèi)星作為觀測衛(wèi)星,通過反射信號反演的海面高度如圖10a所示。在低仰角情況下,即PRN12號衛(wèi)星作為觀測衛(wèi)星,經(jīng)過對異常數(shù)據(jù)刪除之后的反演高度如圖10b所示。在得出1 s結(jié)果之后,對數(shù)據(jù)進(jìn)行20 s平均處理。進(jìn)行平均操作之后的結(jié)果如圖11所示。圖11a是20 s平均后(PRN25號衛(wèi)星)反演的高度圖,圖11b是20 s平均后(PRN12號衛(wèi)星)反演的高度圖。

    圖5 PRN12(紅色)和PRN25(藍(lán)色)衛(wèi)星鏡面反射點軌跡Fig.5 PRN12 (red) and PRN25 (blue) satallite specular reflection point trace

    圖6 第一個數(shù)據(jù)時間段1 s高度對比Fig.6 The first data period 1 s height comparison

    圖7 第一個數(shù)據(jù)時間段20 s平均后高度對比Fig.7 The first data period after 20 s average height comparison

    4.5 數(shù)據(jù)分析總結(jié)

    本文將以偏差(Bias)、平均絕對誤差(Mean Absolute Error, MAE)和平均絕對誤差的標(biāo)準(zhǔn)偏差(Standard Deviations, STD)作為實驗的評估標(biāo)準(zhǔn)。三者的計算公式如下式所示,

    圖8 反演高度異常值(a)及篩選過后高度值(b)Fig.8 Inversion of height anomaly values(a)and filtered height values(b)

    表2為1 s平均,10 s平均,15 s平均和20 s平均的實驗結(jié)果。通過表2中的結(jié)果顯示,在低仰角情況下,由于接收到的信號會受到干擾,所以反演的海面高度結(jié)果比較不穩(wěn)定,在經(jīng)過異常值處理之后,仍然與參考值相比普遍偏大,且離散度高;高仰角下受到的干擾則較小,通過反射信號反演的海面高度與參考值相比誤差較小且離散較低,通過一定時間的平均之后可以消除波浪的影響,提高精度,其中20 s平均的結(jié)果精度最高,可以達(dá)到亞米級的精度,特別是第一階段的數(shù)據(jù),由于高仰角在80°以上,甚至達(dá)到了厘米級的精度,證明了GNSS-R機載海面測高的可行性。

    5 結(jié)論

    本文在岸基GNSS-R海面測高基礎(chǔ)上,建立了機載測高模型,依據(jù)碼延遲測高原理,通過分析2011年11月11日的機載數(shù)據(jù),分別對兩個時間段的試驗數(shù)據(jù)段進(jìn)行分析,對高仰角(75°以上)和低仰角(60°以下)的結(jié)果進(jìn)行了對比分析。

    圖9 直射信號正常功率波形(a)及異常功率波形(b)Fig.9 Direct signal normal(a)and abnormal power waveform(b)

    圖10 第二個數(shù)據(jù)時間段1 s高度對比Fig.10 The second data period 1 s height comparison

    圖11 第二個數(shù)據(jù)時間段20 s平均后高度對比Fig.11 The second data period after 20 s average height comparison

    表2 實驗數(shù)據(jù)結(jié)果Table 2 Result of experimental data

    經(jīng)過數(shù)據(jù)分析,證明了GNSS-R機載技術(shù)可以實現(xiàn)高時間分辨率以及較大范圍的海面高度反演,并且通過采用20 s的數(shù)據(jù)平均方法以后,可以有效地增加反演精度,達(dá)到亞米級精度。

    通過分析同一時刻下,不同反射仰角的反射信號反演的海面高度的結(jié)果分析,高仰角的反演結(jié)果更加穩(wěn)定、準(zhǔn)確?;趯嶒灲Y(jié)果,進(jìn)行機載測高反演時,選取仰角大于75°的數(shù)據(jù)可以減少噪聲干擾以及機身對直射信號的干擾,得到更高精度的反演結(jié)果。同時在低仰角下對直射信號異?,F(xiàn)象進(jìn)行了判斷,對直射信號原始數(shù)據(jù)進(jìn)行篩選后,提高了低仰角下反演結(jié)果的可靠性。本論文中所提出的GNSS-R海面測高模型也可以作為星載測高的理論基礎(chǔ)。

    致謝:感謝CSIC-IEEC提供本次論文的分析數(shù)據(jù),感謝北京航空航天大學(xué)楊東凱教授、CSIC-IEEC的李偉強博士、上海航天電子技術(shù)研究所周勃高級工程師、秦瑾工程師在GNSS-R機載數(shù)據(jù)分析方法上提供的建議。

    猜你喜歡
    仰角海面反演
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    海面床,輕輕搖
    第六章 邂逅“胖胖號”
    用銳角三角函數(shù)解決仰角、俯角問題
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    海面上的“一千座埃菲爾鐵塔”
    分段三次Hermite插值計算GNSS系統(tǒng)衛(wèi)星仰角
    科技視界(2016年15期)2016-06-30 12:48:10
    疊前同步反演在港中油田的應(yīng)用
    欧美不卡视频在线免费观看| 香蕉久久夜色| 国产精品99久久久久久久久| 在线观看免费视频日本深夜| www日本黄色视频网| 精品久久久久久久毛片微露脸| 十八禁人妻一区二区| 欧美丝袜亚洲另类 | 最近最新免费中文字幕在线| 美女 人体艺术 gogo| 黄色丝袜av网址大全| 色老头精品视频在线观看| 成人国产一区最新在线观看| 亚洲av成人av| 黄片大片在线免费观看| 色尼玛亚洲综合影院| 岛国在线观看网站| 亚洲成人精品中文字幕电影| 亚洲狠狠婷婷综合久久图片| 欧美日韩一级在线毛片| 国内精品久久久久精免费| 女同久久另类99精品国产91| 午夜a级毛片| 69人妻影院| 日本黄色片子视频| 日韩欧美在线乱码| 国产三级中文精品| 成年女人永久免费观看视频| 国产欧美日韩一区二区精品| 国产高清三级在线| 波野结衣二区三区在线 | 国产日本99.免费观看| 窝窝影院91人妻| 国产99白浆流出| 国产精品,欧美在线| 中亚洲国语对白在线视频| 欧美成人免费av一区二区三区| 久久久久免费精品人妻一区二区| 精品电影一区二区在线| 哪里可以看免费的av片| 一本一本综合久久| 真实男女啪啪啪动态图| 亚洲成人免费电影在线观看| 黄色女人牲交| 国产精品 国内视频| 欧美极品一区二区三区四区| 免费看光身美女| 精品久久久久久,| 乱人视频在线观看| 非洲黑人性xxxx精品又粗又长| 可以在线观看毛片的网站| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久久大av| 一个人观看的视频www高清免费观看| 国产亚洲欧美98| 在线国产一区二区在线| 亚洲av成人不卡在线观看播放网| 亚洲av日韩精品久久久久久密| 亚洲成人久久爱视频| 啦啦啦韩国在线观看视频| 久久精品91蜜桃| 亚洲精品色激情综合| 久久性视频一级片| 久久久久久久亚洲中文字幕 | 国产精品免费一区二区三区在线| 99riav亚洲国产免费| 高清在线国产一区| 99视频精品全部免费 在线| 国产蜜桃级精品一区二区三区| 亚洲精品影视一区二区三区av| 精品99又大又爽又粗少妇毛片 | 18+在线观看网站| 99国产极品粉嫩在线观看| 舔av片在线| 免费在线观看成人毛片| 亚洲一区二区三区不卡视频| 日韩欧美精品v在线| 欧美日韩乱码在线| 99久久综合精品五月天人人| 岛国在线观看网站| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久人妻精品电影| 亚洲av不卡在线观看| 老熟妇乱子伦视频在线观看| 搡女人真爽免费视频火全软件 | 久久精品亚洲精品国产色婷小说| 搡老岳熟女国产| 成人午夜高清在线视频| 亚洲av电影不卡..在线观看| 观看免费一级毛片| 久久精品91蜜桃| www.999成人在线观看| 亚洲五月婷婷丁香| 天堂网av新在线| 1000部很黄的大片| 婷婷六月久久综合丁香| a级毛片a级免费在线| 国内精品久久久久久久电影| 欧美乱色亚洲激情| 国产一区二区在线观看日韩 | 两性午夜刺激爽爽歪歪视频在线观看| 色噜噜av男人的天堂激情| 日本 欧美在线| 尤物成人国产欧美一区二区三区| 成年女人看的毛片在线观看| 女生性感内裤真人,穿戴方法视频| 热99re8久久精品国产| 亚洲国产日韩欧美精品在线观看 | 久久人妻av系列| 1000部很黄的大片| 久久久精品欧美日韩精品| 国产视频内射| 亚洲国产精品合色在线| 黄色女人牲交| 精品一区二区三区视频在线 | 一个人观看的视频www高清免费观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人免费电影在线观看| 亚洲精品粉嫩美女一区| 亚洲人与动物交配视频| 免费看光身美女| 九九久久精品国产亚洲av麻豆| 禁无遮挡网站| 国产激情偷乱视频一区二区| 九九在线视频观看精品| 悠悠久久av| 一级黄片播放器| 真实男女啪啪啪动态图| 国产亚洲精品一区二区www| 精品人妻偷拍中文字幕| 十八禁网站免费在线| a级毛片a级免费在线| 国产淫片久久久久久久久 | 亚洲五月婷婷丁香| 级片在线观看| 九九在线视频观看精品| 黄片小视频在线播放| 色综合站精品国产| 国产精品 国内视频| 亚洲国产精品成人综合色| 午夜福利免费观看在线| 十八禁网站免费在线| 色吧在线观看| 久久亚洲精品不卡| 国产一区在线观看成人免费| av中文乱码字幕在线| 一区二区三区高清视频在线| 老汉色∧v一级毛片| 日韩中文字幕欧美一区二区| 亚洲色图av天堂| 亚洲av二区三区四区| 老熟妇乱子伦视频在线观看| 久久久久国产精品人妻aⅴ院| 国产亚洲精品av在线| 亚洲国产欧美网| 午夜日韩欧美国产| 精品国内亚洲2022精品成人| 免费观看精品视频网站| 亚洲国产日韩欧美精品在线观看 | 99在线人妻在线中文字幕| 国产成人影院久久av| 青草久久国产| 特大巨黑吊av在线直播| 波多野结衣巨乳人妻| av在线蜜桃| 欧美在线黄色| 亚洲熟妇中文字幕五十中出| 国产精品一及| 久久亚洲精品不卡| 男女做爰动态图高潮gif福利片| 亚洲精品成人久久久久久| eeuss影院久久| 亚洲自拍偷在线| 91av网一区二区| 一a级毛片在线观看| 搡女人真爽免费视频火全软件 | 岛国视频午夜一区免费看| 国产蜜桃级精品一区二区三区| 国产主播在线观看一区二区| 欧美bdsm另类| av欧美777| 亚洲片人在线观看| 国产精品av视频在线免费观看| 久久精品国产亚洲av香蕉五月| 国产91精品成人一区二区三区| 性色avwww在线观看| 在线播放国产精品三级| 欧美日韩精品网址| 精品国产三级普通话版| 一a级毛片在线观看| 久久久久久大精品| 久久久色成人| 亚洲欧美日韩高清在线视频| 亚洲精品在线美女| 99精品欧美一区二区三区四区| 欧美在线一区亚洲| 美女高潮的动态| 国产精品乱码一区二三区的特点| 五月伊人婷婷丁香| 国产亚洲欧美在线一区二区| 18禁裸乳无遮挡免费网站照片| 精品免费久久久久久久清纯| 欧美3d第一页| 国产伦在线观看视频一区| 欧美乱色亚洲激情| 在线视频色国产色| 亚洲av一区综合| 国产欧美日韩一区二区精品| 午夜视频国产福利| 嫩草影视91久久| 床上黄色一级片| 亚洲国产色片| 中文字幕久久专区| 91字幕亚洲| 精品久久久久久久毛片微露脸| x7x7x7水蜜桃| 国产亚洲精品综合一区在线观看| 日韩人妻高清精品专区| 国产精品三级大全| 99在线视频只有这里精品首页| 午夜精品一区二区三区免费看| 麻豆久久精品国产亚洲av| 国产国拍精品亚洲av在线观看 | 中文字幕av在线有码专区| 人人妻,人人澡人人爽秒播| 黄色视频,在线免费观看| 久久久国产成人精品二区| 天天一区二区日本电影三级| 一区福利在线观看| 最后的刺客免费高清国语| 久久久精品大字幕| 久久久久久九九精品二区国产| 99在线人妻在线中文字幕| 国内精品久久久久久久电影| 国产探花极品一区二区| 精品国产三级普通话版| 亚洲五月婷婷丁香| 欧美最黄视频在线播放免费| 日本五十路高清| 制服丝袜大香蕉在线| 亚洲男人的天堂狠狠| 12—13女人毛片做爰片一| 久久久久久国产a免费观看| 日本撒尿小便嘘嘘汇集6| 成人18禁在线播放| 欧美日韩黄片免| 日韩欧美免费精品| 最近最新中文字幕大全电影3| 熟女电影av网| 国产欧美日韩一区二区三| 51国产日韩欧美| 99国产精品一区二区三区| 草草在线视频免费看| 久久久色成人| 真人做人爱边吃奶动态| 日韩欧美国产一区二区入口| 91久久精品国产一区二区成人 | av天堂中文字幕网| 波多野结衣高清作品| 亚洲一区高清亚洲精品| 高潮久久久久久久久久久不卡| 免费无遮挡裸体视频| 欧美日本亚洲视频在线播放| 国产黄色小视频在线观看| svipshipincom国产片| 精品人妻偷拍中文字幕| 最近最新中文字幕大全免费视频| 国产精品永久免费网站| 动漫黄色视频在线观看| 乱人视频在线观看| 校园春色视频在线观看| 精品一区二区三区视频在线观看免费| 久久久国产精品麻豆| 18美女黄网站色大片免费观看| 国产亚洲av嫩草精品影院| 日日夜夜操网爽| 亚洲专区中文字幕在线| 久久久久久久精品吃奶| 欧美成人免费av一区二区三区| 亚洲在线自拍视频| 免费在线观看影片大全网站| 久久久国产精品麻豆| 色综合欧美亚洲国产小说| 亚洲欧美日韩东京热| 国产精品av视频在线免费观看| 午夜福利视频1000在线观看| 人妻丰满熟妇av一区二区三区| 日韩欧美国产一区二区入口| 亚洲七黄色美女视频| 99久国产av精品| 国内精品美女久久久久久| 欧美日韩乱码在线| 精品久久久久久久久久免费视频| 最好的美女福利视频网| 国产精品久久久久久人妻精品电影| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本黄色片子视频| 欧美最新免费一区二区三区 | 亚洲成人久久性| 19禁男女啪啪无遮挡网站| 91麻豆av在线| 成人特级黄色片久久久久久久| 九色成人免费人妻av| 99久久综合精品五月天人人| 成年版毛片免费区| 国产国拍精品亚洲av在线观看 | av福利片在线观看| 老汉色av国产亚洲站长工具| 又爽又黄无遮挡网站| 亚洲 欧美 日韩 在线 免费| 我要搜黄色片| xxx96com| 久久欧美精品欧美久久欧美| 国产精品女同一区二区软件 | 两个人的视频大全免费| 99在线人妻在线中文字幕| 午夜激情欧美在线| 欧美精品啪啪一区二区三区| 国产成年人精品一区二区| 18禁黄网站禁片午夜丰满| 亚洲欧美日韩东京热| 国产精品,欧美在线| 18禁黄网站禁片午夜丰满| www日本在线高清视频| 国产伦精品一区二区三区视频9 | 可以在线观看毛片的网站| 亚洲乱码一区二区免费版| 中亚洲国语对白在线视频| 亚洲av中文字字幕乱码综合| 色噜噜av男人的天堂激情| 真人一进一出gif抽搐免费| 长腿黑丝高跟| 国产日本99.免费观看| 99视频精品全部免费 在线| 69av精品久久久久久| 国产aⅴ精品一区二区三区波| 国产精品影院久久| x7x7x7水蜜桃| 麻豆成人av在线观看| 久久久久久久午夜电影| 美女被艹到高潮喷水动态| 长腿黑丝高跟| 成人一区二区视频在线观看| 国产精品一区二区三区四区久久| 精品久久久久久久久久免费视频| 免费电影在线观看免费观看| 一本精品99久久精品77| 免费电影在线观看免费观看| 最新美女视频免费是黄的| 婷婷六月久久综合丁香| 亚洲国产精品久久男人天堂| 99久久精品热视频| av黄色大香蕉| 国产欧美日韩一区二区三| www.熟女人妻精品国产| 国产欧美日韩一区二区三| 国产av一区在线观看免费| 国产高清激情床上av| 蜜桃久久精品国产亚洲av| 午夜福利欧美成人| 精品久久久久久成人av| 搡女人真爽免费视频火全软件 | 日韩欧美精品v在线| 97超级碰碰碰精品色视频在线观看| 亚洲精品在线观看二区| 国产视频内射| 日韩欧美精品v在线| 女人被狂操c到高潮| 久久人人精品亚洲av| 欧美色欧美亚洲另类二区| 波多野结衣高清无吗| 真人一进一出gif抽搐免费| 色哟哟哟哟哟哟| 国产欧美日韩一区二区精品| 久久久久久九九精品二区国产| 精品熟女少妇八av免费久了| 又粗又爽又猛毛片免费看| 一进一出抽搐动态| 成人国产一区最新在线观看| 欧美色欧美亚洲另类二区| 最好的美女福利视频网| 老司机午夜福利在线观看视频| 亚洲精品在线观看二区| 亚洲av免费高清在线观看| 中文字幕人妻熟人妻熟丝袜美 | 成人一区二区视频在线观看| 国产精品影院久久| 十八禁人妻一区二区| www日本在线高清视频| 亚洲精品粉嫩美女一区| 久久伊人香网站| 午夜福利18| 国产成人av教育| 国产精品亚洲一级av第二区| 免费看美女性在线毛片视频| 日韩大尺度精品在线看网址| 亚洲精品成人久久久久久| 男人和女人高潮做爰伦理| 99精品在免费线老司机午夜| 97超视频在线观看视频| 老司机午夜福利在线观看视频| 波多野结衣高清无吗| 国产亚洲精品久久久com| 欧美日韩瑟瑟在线播放| 舔av片在线| 国产成人a区在线观看| 五月伊人婷婷丁香| 18禁黄网站禁片午夜丰满| 18+在线观看网站| 国产真实伦视频高清在线观看 | 国产毛片a区久久久久| 国产伦精品一区二区三区视频9 | 国产熟女xx| 少妇的逼好多水| 高清在线国产一区| 国产真实伦视频高清在线观看 | 九九久久精品国产亚洲av麻豆| 热99re8久久精品国产| 国产99白浆流出| 一夜夜www| 免费av观看视频| 国产成年人精品一区二区| 亚洲avbb在线观看| 亚洲精品在线观看二区| 国产淫片久久久久久久久 | 日韩大尺度精品在线看网址| 一区福利在线观看| 99久久精品国产亚洲精品| 成年女人看的毛片在线观看| 啪啪无遮挡十八禁网站| 国产亚洲精品综合一区在线观看| 美女黄网站色视频| 成人鲁丝片一二三区免费| 身体一侧抽搐| 狠狠狠狠99中文字幕| 亚洲av中文字字幕乱码综合| 国产精品免费一区二区三区在线| 亚洲avbb在线观看| 国产久久久一区二区三区| 国产午夜精品论理片| 国产一区二区在线av高清观看| 精品久久久久久久久久久久久| 亚洲精品日韩av片在线观看 | 91九色精品人成在线观看| 女生性感内裤真人,穿戴方法视频| 男女那种视频在线观看| 国产高清激情床上av| 老司机午夜福利在线观看视频| 国产精品影院久久| 亚洲成人久久爱视频| 999久久久精品免费观看国产| 最近最新免费中文字幕在线| 亚洲国产精品合色在线| 午夜免费激情av| 精品电影一区二区在线| 亚洲成a人片在线一区二区| 国产一区在线观看成人免费| 两个人看的免费小视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久久久人妻精品电影| 一区二区三区激情视频| 美女 人体艺术 gogo| 国产真实伦视频高清在线观看 | 99久国产av精品| www.色视频.com| 夜夜躁狠狠躁天天躁| or卡值多少钱| 国产黄片美女视频| 一级作爱视频免费观看| 丁香欧美五月| 99热这里只有是精品50| 国产一区二区三区在线臀色熟女| 又黄又爽又免费观看的视频| 毛片女人毛片| 日本与韩国留学比较| 国产久久久一区二区三区| 国产精品一区二区三区四区久久| 国产伦在线观看视频一区| 亚洲av一区综合| 国产精品亚洲美女久久久| 超碰av人人做人人爽久久 | 床上黄色一级片| 丰满的人妻完整版| x7x7x7水蜜桃| 身体一侧抽搐| 午夜精品久久久久久毛片777| 国产午夜福利久久久久久| 午夜老司机福利剧场| 国产精品98久久久久久宅男小说| 俄罗斯特黄特色一大片| 午夜日韩欧美国产| 中文字幕熟女人妻在线| 国产伦精品一区二区三区四那| 久久精品影院6| 免费在线观看日本一区| 五月玫瑰六月丁香| 亚洲av成人精品一区久久| 国产私拍福利视频在线观看| 99久久成人亚洲精品观看| 丰满的人妻完整版| 精品久久久久久久久久久久久| 欧美午夜高清在线| 中文在线观看免费www的网站| 久久久久久久久大av| 国产又黄又爽又无遮挡在线| 在线观看一区二区三区| 内地一区二区视频在线| 国产高清激情床上av| www.熟女人妻精品国产| 久久精品国产清高在天天线| 级片在线观看| 一区福利在线观看| av黄色大香蕉| 亚洲性夜色夜夜综合| 日本成人三级电影网站| 在线看三级毛片| 一进一出好大好爽视频| 在线播放无遮挡| 免费看十八禁软件| av中文乱码字幕在线| 精品人妻一区二区三区麻豆 | 婷婷丁香在线五月| 久久精品影院6| 欧美精品啪啪一区二区三区| 99在线视频只有这里精品首页| 97人妻精品一区二区三区麻豆| 国产精品一区二区免费欧美| 欧美日韩乱码在线| 精品一区二区三区人妻视频| 久久久国产成人精品二区| 亚洲av美国av| 国产探花极品一区二区| 麻豆国产av国片精品| 丁香欧美五月| 欧美最黄视频在线播放免费| 91久久精品国产一区二区成人 | 十八禁人妻一区二区| 国产男靠女视频免费网站| 高清在线国产一区| 国产精品精品国产色婷婷| 国产亚洲精品久久久久久毛片| 一级黄片播放器| 久久99热这里只有精品18| 国产伦精品一区二区三区四那| 色播亚洲综合网| 成人国产一区最新在线观看| 尤物成人国产欧美一区二区三区| 嫩草影院精品99| 小说图片视频综合网站| 亚洲欧美激情综合另类| 亚洲激情在线av| 国产熟女xx| 亚洲最大成人手机在线| 18禁黄网站禁片免费观看直播| 国产探花在线观看一区二区| 免费人成视频x8x8入口观看| ponron亚洲| 国产成人啪精品午夜网站| 淫秽高清视频在线观看| 夜夜看夜夜爽夜夜摸| 老司机午夜福利在线观看视频| 亚洲av成人av| 成人国产综合亚洲| 一级黄片播放器| 国产精品,欧美在线| 欧美绝顶高潮抽搐喷水| 99精品久久久久人妻精品| 婷婷六月久久综合丁香| 欧美乱码精品一区二区三区| 久久午夜亚洲精品久久| 最新中文字幕久久久久| 波多野结衣巨乳人妻| 99riav亚洲国产免费| 亚洲国产精品999在线| 亚洲欧美日韩高清在线视频| 免费人成视频x8x8入口观看| 久久久久亚洲av毛片大全| 国产 一区 欧美 日韩| 我要搜黄色片| 成人av一区二区三区在线看| 国产精品永久免费网站| 日韩欧美免费精品| 日韩欧美 国产精品| 成人18禁在线播放| 日韩欧美免费精品| 一区二区三区高清视频在线| 亚洲第一电影网av| 超碰av人人做人人爽久久 | 午夜老司机福利剧场| 亚洲欧美日韩无卡精品| av专区在线播放| 最后的刺客免费高清国语| 国产成人系列免费观看| 天堂av国产一区二区熟女人妻| 日韩欧美精品免费久久 | 亚洲aⅴ乱码一区二区在线播放| 亚洲精品亚洲一区二区| 黑人欧美特级aaaaaa片| 亚洲欧美一区二区三区黑人| 大型黄色视频在线免费观看| 国产 一区 欧美 日韩| 女人十人毛片免费观看3o分钟| 男人舔奶头视频| 精品一区二区三区av网在线观看| 久久久久久久精品吃奶| 人人妻,人人澡人人爽秒播| 亚洲欧美日韩高清在线视频| 欧美成人a在线观看| 色在线成人网| 国产伦精品一区二区三区四那| 俺也久久电影网| 色噜噜av男人的天堂激情| 中文资源天堂在线| 最近在线观看免费完整版|