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

    北斗單歷元基線解算與恒星日濾波算法

    2015-12-26 06:22:20韓曉飛,馬緒瀛,丁曉光
    測繪通報(bào) 2015年4期

    引文格式: 韓曉飛,馬緒瀛,丁曉光,等. 北斗單歷元基線解算與恒星日濾波算法[J].測繪通報(bào),2015(4):5-9.DOI:10.13474/j.cnki.11-2246.2015.0099

    北斗單歷元基線解算與恒星日濾波算法

    韓曉飛1,馬緒瀛2,丁曉光1,張永奇1,蘇利娜1

    (1. 陜西省地震局,陜西 西安 710068; 2. 同濟(jì)大學(xué),上海 200092)

    SingleEpochBaselineSolutionandSiderealFilteringofBeiDou

    HANXiaofei,MAXuying,DINGXiaoguang,ZHANGYongqi,SULina

    摘要:對單歷元相對定位解算方法進(jìn)行了研究,編制了相應(yīng)的北斗數(shù)據(jù)處理軟件,對兩條1km以內(nèi)的短基線數(shù)據(jù)進(jìn)行了測試。根據(jù)基線解算結(jié)果提取了測站的多路徑誤差時(shí)間序列,并通過最大相關(guān)性分析確定了多路徑誤差的周期,在此基礎(chǔ)上提出了適用于北斗系統(tǒng)的恒星日濾波算法。

    關(guān)鍵詞:北斗衛(wèi)星導(dǎo)航系統(tǒng);單歷元基線解算;多路徑誤差;恒星日濾波

    中圖分類號:P228.4

    收稿日期:2014-04-21

    作者簡介:韓曉飛(1982—),男,工程師,主要從事地震觀測、數(shù)據(jù)分析等工作。E-mail:voldemortpotter@sina.com

    一、引言

    北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDouNavigationSatelliteSystem,BDS)是中國自行研制的繼美國全球定位系統(tǒng)(GPS)和俄羅斯GLONASS之后第3個(gè)成熟的衛(wèi)星導(dǎo)航系統(tǒng)[1]。2012年底第13、14號衛(wèi)星的成功運(yùn)行標(biāo)志著覆蓋亞太地區(qū)的北斗區(qū)域衛(wèi)星導(dǎo)航系統(tǒng)已經(jīng)建設(shè)完成,偽距定位精度可達(dá)10m以內(nèi),并且初步具備高精度定位能力。GNSS系統(tǒng)具有高精度定位能力,逐漸在各類地形、地貌、地震預(yù)警和建筑物變形監(jiān)測等方面得到廣泛應(yīng)用。隨著我國北斗系統(tǒng)的完善與高頻(1Hz)、超高頻(20~50Hz)GNSS定位技術(shù)的進(jìn)一步發(fā)展,基于相位觀測值的北斗單歷元相對定位技術(shù)可以快速、準(zhǔn)確地獲得接收機(jī)的位置信息,它可以應(yīng)用于地震學(xué)、火山學(xué)等領(lǐng)域,越來越多的國內(nèi)外學(xué)者對其數(shù)據(jù)處理方法和應(yīng)用前景產(chǎn)生了濃厚的興趣[2]。相對定位技術(shù)可以大幅度削弱相關(guān)誤差,獲得高精度的定位結(jié)果。但是差分技術(shù)無法消除多路徑誤差,在局部監(jiān)測中,布設(shè)的基線較短,在消除其他強(qiáng)相關(guān)性的誤差后,多路徑誤差便成為影響監(jiān)測結(jié)果的主要誤差源[3-6]。GNSS系統(tǒng)星座多由軌道具有周期特性的衛(wèi)星構(gòu)成,靜態(tài)觀測站的多路徑等環(huán)境誤差將會呈現(xiàn)出周期特性,因此可利用前一天同時(shí)段的多路徑誤差改進(jìn)后一天定位結(jié)果的方法(稱為恒星日濾波)[5-6]。本文探索了北斗系統(tǒng)的多路徑誤差周期,在獲得其周期后,給出了基于北斗系統(tǒng)的恒星日濾波方法。試驗(yàn)結(jié)果顯示,該方法對北斗多路徑效應(yīng)的消除和衛(wèi)星地形監(jiān)測具有重要的意義。

    二、北斗單歷元基線解算

    單歷元基線解算屬于動(dòng)態(tài)相對定位的一種,其觀測方程與靜態(tài)相對定位是相同的[7],因?yàn)槭菃螝v元解算,考慮到單歷元解算若僅采用載波相位觀測值會導(dǎo)致法方程秩虧,本文同時(shí)加入偽距觀測值組建雙頻偽距及載波相位觀測方程。

    1. 數(shù)學(xué)模型

    進(jìn)行雙差后,可消除接收機(jī)和衛(wèi)星鐘差,其觀測方程為[7]

    (1)

    (2)

    式中,下標(biāo)A、B為測站;上標(biāo)i、j為觀測衛(wèi)星。其中,A為基準(zhǔn)站;i為參考衛(wèi)星;下標(biāo)組合表示對測站作差,上標(biāo)組合表示對衛(wèi)星作差;P為偽距觀測值;dX、dY和dZ分別為基線向量的3個(gè)坐標(biāo)改正數(shù);l、m和n分別為3個(gè)坐標(biāo)方向上的方向余弦;λf是f頻率載波所對應(yīng)的波長;N是載波信號的整周模糊度;I和T分別為電離層延遲改正和對流層延遲改正,當(dāng)基線長度較短時(shí)此兩項(xiàng)誤差可忽略不計(jì)。常數(shù)項(xiàng)為

    (3)

    式中,ρ為測站到衛(wèi)星之間的距離;MP、ML分別為偽距和相位的多路徑誤差;εP、εL分別為偽距和相位的觀測噪聲;雙差觀測值的隨機(jī)模型為[8]

    (4)

    式中,s表示觀測值類型;ns為該類觀測值的衛(wèi)星數(shù)。

    2. 模糊度固定方法

    進(jìn)行快速定位時(shí),雖然觀測時(shí)間較短,但只要能正確確定整周模糊度,依然能得到高精度的定位結(jié)果,因此,基于載波相位觀測值的單歷元基線解算的關(guān)鍵在于快速確定整周模糊度,文中采用LAMBDA法進(jìn)行模糊度搜索。

    本文采用了較為成熟的onthefly(OTF)技術(shù)進(jìn)行單歷元短基線動(dòng)態(tài)定位解算。該解算方法可以概括為兩部分:模糊度參數(shù)的初始化和固定模糊度后的單歷元定位解算[9]。

    3. 基線解算算例分析

    采用和芯星通(UNICORE)公司生產(chǎn)的UB240-CORS北斗/GPS雙系統(tǒng)四頻接收機(jī)(GPS:L1、L2;北斗:B1、B2)對兩條小于1km的短基線進(jìn)行了數(shù)據(jù)采集,一條位于北京的超短基線3.0669m,數(shù)據(jù)時(shí)間段為2012年11月8—10日;另一條位于上海同濟(jì)大學(xué)的基線,長度為470.3009m,數(shù)據(jù)時(shí)間段為2012年11月8—15日,數(shù)據(jù)的采樣間隔均為1s。

    利用自行編制的北斗數(shù)據(jù)處理軟件,按1s采樣率分別解算了同濟(jì)大學(xué)測站11月8—15日、北京測站11月8—10日的實(shí)測數(shù)據(jù)。同時(shí)利用Bernese軟件對兩個(gè)地區(qū)的GPS數(shù)據(jù)進(jìn)行了解算,將解算得到的觀測站坐標(biāo)作為真值開展后續(xù)分析。由于篇幅限制,本文只給出上海站11月8日一天觀測數(shù)據(jù)(采樣率1s,解算歷元86400個(gè))求得的北斗基線相對于標(biāo)準(zhǔn)值的誤差序列,如圖1所示。其他時(shí)間段觀測數(shù)據(jù)的解算結(jié)果如圖2所示。

    圖1 上海站2012年11月8日北斗基線解算誤差時(shí)間序列

    由圖1的誤差序列求得北斗基線在北(N)、東(E)、高程(U)各方向上的RMS值分別為0.41cm、0.30cm和0.71cm。

    為了評估北斗系統(tǒng)在上海、北京兩站的定位表現(xiàn),對多天解算結(jié)果進(jìn)行了統(tǒng)計(jì)分析,如圖2所示。經(jīng)統(tǒng)計(jì),上海站6天北斗解算結(jié)果N、E、U方向誤差序列的平均RMS值分別為0.42cm、0.31cm、0.80cm;北京站3天北斗解算結(jié)果N、E、U方向誤差序列的平均RMS值分別為0.21cm、0.18cm、0.55cm。可以看到,對于這兩條短基線數(shù)據(jù),北斗單歷元基線解算得到了較好的結(jié)果。

    圖2 上海(上)、北京(下)站北斗系統(tǒng)多天基線解算結(jié)果坐標(biāo)三分量RMS統(tǒng)計(jì)

    三、北斗多路徑誤差周期特性分析

    多路徑效應(yīng)是指北斗信號通過不止一條不同路徑到達(dá)接收機(jī)天線的現(xiàn)象,即對所期望信號經(jīng)反射或散射后的復(fù)制品的接收。多路徑誤差不僅與觀測環(huán)境有關(guān),同時(shí)還與衛(wèi)星的視運(yùn)動(dòng)有關(guān)。對于靜止的觀測站,其周邊環(huán)境長期保持不變,而導(dǎo)航衛(wèi)星繞地球軌道進(jìn)行周期性運(yùn)動(dòng),因此觀測站接收機(jī)的多路徑誤差具有周期特性,同時(shí)與衛(wèi)星軌道信息高度相關(guān)。

    1. 基于基線解算結(jié)果的北斗多路徑誤差提取

    定位解算后的時(shí)間序列主要包含載波相位多路徑誤差和隨機(jī)噪聲,因此該數(shù)據(jù)可以分析站點(diǎn)相位多路徑誤差的綜合影響。相對多路徑信號,隨機(jī)噪聲信號具有高頻特性,多路徑誤差具有低頻特性,利用隨機(jī)噪聲和多路徑信號的頻率特性進(jìn)行小波濾波提取多路徑誤差。本文采用db10小波進(jìn)行7層的小波包分解,剔除高頻噪聲信號后合成去噪后的新信號。

    2. 北斗多路徑誤差周期特性

    為了從理論上探討多路徑誤差的相關(guān)性,計(jì)算連續(xù)多天多路徑誤差時(shí)間序列的相關(guān)性。

    算例采用北京、上海兩站11月8—10日連續(xù)3天的基線數(shù)據(jù),先分別進(jìn)行北斗單歷元基線解算,依次獲得第1天、第2天、第3天N、E、U方向誤差的時(shí)間序列,然后對原始解算結(jié)果進(jìn)行小波去噪濾波提取相對純凈的多路徑誤差(見表1)。

    從表1結(jié)果可見,相鄰兩天在相關(guān)延遲約為240s附近出現(xiàn)最大相關(guān),略大于相關(guān)延遲的理論值236s,隨著間隔天數(shù)的增加,最大相關(guān)系數(shù)值會逐漸減小[10]。如第1天與第3天的相關(guān)值比第1天與第2天或第2天與第3天這類相鄰兩天的相關(guān)值都要低。說明同一地點(diǎn),隨著天數(shù)增加,衛(wèi)星幾何條件在發(fā)生變化。如上所述,對于靜止觀測站,北斗多路徑效應(yīng)的周期約為86160s,略少于24h。

    表1 北京、上海站北斗多路徑效應(yīng)相關(guān)延遲統(tǒng)計(jì)

    3. 多路徑誤差周期與衛(wèi)星軌道周期的關(guān)系

    本文介紹了兩種方法來確定各個(gè)衛(wèi)星的軌道周期:廣播星歷法與精密星歷法。

    (1) 通過廣播星歷計(jì)算軌道重復(fù)周期

    根據(jù)廣播星歷提供的長半軸的平方根和平均角速度攝動(dòng)參數(shù),由開普勒第三定律得

    (5)

    式中,a為衛(wèi)星軌道橢圓長半軸的平方根;GM為萬有引力常數(shù)G與地球總質(zhì)量M的乘積,其值為:GM=398600.4418×109m3/s2;Δn為平均運(yùn)動(dòng)角速度的改正量;T0為衛(wèi)星運(yùn)行一周的時(shí)間;TR為軌道重復(fù)周期。利用2012年6月13日至2012年8月1日這段時(shí)間的廣播星歷求得在這段時(shí)間內(nèi)北斗衛(wèi)星的軌道重復(fù)周期。

    圖3所示分別為GEO、IGSO和MEO三類衛(wèi)星的軌道重復(fù)周期,其中1—5號、6—10號和11、12號衛(wèi)星分別為GEO、IGSO和MEO衛(wèi)星。可以看出,每一類衛(wèi)星的軌道重復(fù)周期都各不相同,對于同一顆衛(wèi)星來說,每天的運(yùn)動(dòng)周期也不完全一致。但是對于每一類衛(wèi)星,其軌道重復(fù)周期都在某一常數(shù)附近波動(dòng)。根據(jù)統(tǒng)計(jì),GEO、IGSO和MEO衛(wèi)星的平均軌道重復(fù)周期分別為86163s、86162s、46391s,與以前求得的重復(fù)周期基本一致[10]。即衛(wèi)星每天提前時(shí)間分別為237s、238s和245s(北斗MEO衛(wèi)星與GPS系統(tǒng)的MEO衛(wèi)星有較大差異,與GPS衛(wèi)星11h 58min的運(yùn)行周期不同,北斗MEO衛(wèi)星7d繞地球運(yùn)行13圈,7d累計(jì)提前上一周1715s,平均每天提前245s,該提前時(shí)間與GPS的MEO衛(wèi)星結(jié)果基本一致)。

    圖3 通過2012年6月13日至8月1日廣播星歷獲得的全部北斗衛(wèi)星軌道重復(fù)周期

    (2) 通過精密星歷計(jì)算軌道重復(fù)周期

    利用廣播星歷可以方便、快速地得到衛(wèi)星的軌道重復(fù)周期,然而利用精密星歷求得的軌道重復(fù)周期更加準(zhǔn)確,同時(shí)還克服了廣播星歷計(jì)算時(shí)并未考慮測站位置信息影響的不足(廣播星歷計(jì)算的重復(fù)周期僅僅體現(xiàn)了衛(wèi)星運(yùn)行一周的時(shí)間)。對于長期靜止的觀測站而言,為了進(jìn)一步的數(shù)據(jù)處理(如恒星日濾波),我們更關(guān)心的是衛(wèi)星相鄰2天在測站上空同一位置出現(xiàn)的時(shí)間間隔(即地面跟蹤重復(fù)時(shí)間,Duncan Carr Agnew等人在《Finding the repeat times of the GPS constellation》中將它定義為Aspect repeat time,也就是視方位重復(fù)時(shí)間[10])。

    由精密星歷計(jì)算衛(wèi)星軌道重復(fù)周期的步驟為:

    1) 收集觀測站的測站坐標(biāo)(X0,Y0,Z0)及要處理的數(shù)據(jù)相鄰2d的精密星歷。

    2) 指定初始參考時(shí)間Tr,然后根據(jù)精密星歷計(jì)算Tr時(shí)刻各個(gè)北斗衛(wèi)星的空中位置坐標(biāo)(xi,yi,zi)及測站和衛(wèi)星構(gòu)成的各個(gè)方向向量ni。

    表2中通過精密星歷計(jì)算了2012年11月11日各顆北斗衛(wèi)星的軌道重復(fù)周期。其中,SatID為衛(wèi)星編號,Tp為精密星歷計(jì)算結(jié)果。有一點(diǎn)需要指出的是,由于北斗MEO衛(wèi)星7d繞地球運(yùn)行13圈(一周運(yùn)行時(shí)間大約12.92h),無法像GPS衛(wèi)星一樣以天進(jìn)行對比,一般以一周7d為周期進(jìn)行比較,因此在表2中未進(jìn)行對比。

    對比以上兩種方法的結(jié)果,基于廣播星歷與精密星歷計(jì)算得到的衛(wèi)星軌道重復(fù)周期與靜止觀測站北斗多路徑誤差周期大致相同,驗(yàn)證了多路徑誤差周期與衛(wèi)星軌道信息高度相關(guān)這一結(jié)論。

    表2 精密星歷計(jì)算的各顆衛(wèi)星軌道重復(fù)周期結(jié)果 s

    四、基于多路徑誤差周期的恒星日濾波

    GPS系統(tǒng)星座由多顆MEO衛(wèi)星構(gòu)成,衛(wèi)星軌道運(yùn)行周期為11h 58min。接收機(jī)位置不變時(shí),第二天提前236s時(shí)刻的衛(wèi)星與接收機(jī)天線相位中心的幾何構(gòu)形在理論上與前一天完全相同,靜態(tài)觀測站的多路徑等環(huán)境誤差將會呈現(xiàn)出周期特性,因此可利用前一天同時(shí)段的多路徑誤差改進(jìn)后一天定位結(jié)果的方法稱為恒星日濾波[5-6,11]。

    北斗系統(tǒng)衛(wèi)星星座的構(gòu)成與GPS系統(tǒng)差異較大,其星座由GEO、IGSO、MEO三種不同類型的衛(wèi)星構(gòu)成,導(dǎo)致其衛(wèi)星軌道的重復(fù)特性也與GPS相異,因此基于GPS星座特性的恒星日濾波方法不適用于北斗衛(wèi)星系統(tǒng)。根據(jù)上文北斗多路徑誤差周期計(jì)算的結(jié)果綜合考慮,北斗衛(wèi)星的平均提前時(shí)間約為240s。

    算例根據(jù)前文得到的平均重復(fù)周期來對基線解算結(jié)果進(jìn)行恒星日濾波。所用數(shù)據(jù)為上海、北京站連續(xù)2d的短基線數(shù)據(jù)。北斗系統(tǒng)濾波前后的坐標(biāo)三分量結(jié)果如圖4所示(由于篇幅限制,只繪制了上海站北斗系統(tǒng)的濾波結(jié)果圖),其他解算結(jié)果見表3、表4。

    圖4 2012年11月14日上海站濾波前后的坐標(biāo)三分量時(shí)間序列

    RMSNRMSERMSU原始數(shù)據(jù)解算結(jié)果/cm0.510.341.06恒星日濾波結(jié)果/cm0.430.270.82提升比例/(%)15.6920.5922.64

    表4 北京站濾波前后的結(jié)果統(tǒng)計(jì)

    可以看出,恒星日濾波消除了部分周期誤差,坐標(biāo)時(shí)間序列的精度有了明顯的提高,這說明此方法在北斗數(shù)據(jù)處理中是具有實(shí)際意義的。

    五、結(jié)論

    本文就北斗單歷元相對定位解算方法進(jìn)行了研究,對兩條1km以內(nèi)的短基線數(shù)據(jù)進(jìn)行了測試。根據(jù)基線解算結(jié)果提取了測站的多路徑誤差時(shí)間序列,并通過最大相關(guān)性分析確定了多路徑誤差的周期,在此基礎(chǔ)上提出了適用于北斗系統(tǒng)的恒星日濾波算法。得出以下結(jié)論:

    1) 對于1km以內(nèi)的短基線,北斗單歷元相對定位解算可以達(dá)到厘米級甚至毫米級的定位結(jié)果,北斗系統(tǒng)已經(jīng)具備高精度定位能力。

    2) 對于靜止觀測站,北斗多路徑效應(yīng)的周期約為86160s,與衛(wèi)星軌道信息高度相關(guān)。

    3) 北斗恒星日濾波算法消除了部分周期誤差,坐標(biāo)時(shí)間序列的精度有了10%~20%的改進(jìn),此方法在北斗數(shù)據(jù)處理中具有實(shí)際意義。

    參考文獻(xiàn):

    [1]中國衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室. 北斗衛(wèi)星導(dǎo)航系統(tǒng)空間信號接口控制文件(測試版)[R]. 北京: [s.n.], 2011.

    [2]殷海濤, 甘衛(wèi)軍, 肖根如. 恒星日濾波的修正以及對高頻GPS定位的影響研究[J]. 武漢大學(xué)學(xué)報(bào):信息科學(xué)版, 2011, 36(5): 609-611.

    [3]劉思慧, 樊婧琦, 歐鋼, 等. 多徑效應(yīng)的動(dòng)態(tài)特性對碼跟蹤環(huán)路的影響分析[J]. 測繪學(xué)報(bào), 2011, 40(S0): 89-94.

    [4]王夢麗, 王飛雪. 三頻電離層延遲改正中多路徑誤差和觀測噪聲的削弱算法[J]. 測繪學(xué)報(bào), 2008, 37(4): 418-422.

    [5]CHOI K, BILICH A, LARSON K M, et al. Modified Sidereal Filtering: Implications for High-Rate GPS Positioning[J]. Geophysical Research Letters,2004,31(22). DOI: 10.1029/2004GL021621.

    [6]RAGHEB E A, CLARKE P J, EDWARDS S J. GPS Sidereal Filtering: Coordinate and Carrier Phase Level Strategies[J]. Journal of Geodesy,2007,81(5): 325-335. DOI: 10.1007/s00190-006-0113-1.

    [7]李征航, 黃勁松. GPS測量與數(shù)據(jù)處理[M]. 武漢: 武漢大學(xué)出版社, 2005.

    [8]魏子卿, 葛茂榮. GPS相對定位的數(shù)學(xué)模型[M]. 北京: 測繪出版社, 1998.

    [9]胡叢瑋, 劉大杰. 單歷元確定GPS整周模糊度的分析[J]. 南京航空航天大學(xué)學(xué)報(bào), 2001, 33(3): 267-271.

    [10]MA X Y, SHEN Y Z. Multipath Error Analysis of COMPASS Triple Frequency Observations[J].Positioning,2014,5(1):12-21.DOI:10.42361/pos.2014.51002.

    [11]AGNEW D C, LARSON K M. Finding the Repeat Times of the GPS Constellation[J]. GPS Solutions, 2006, 11(3): 71-76.DOI: 10.1007/s10291-006-0038-4.

    2015北京國際城市管線展覽會將于6月在京舉辦

    [本刊訊]由中國城市規(guī)劃協(xié)會地下管線專業(yè)委員會和商務(wù)部經(jīng)濟(jì)技術(shù)交流中心主辦的2015北京國際城市管線展覽會(官方網(wǎng)站:www.dxguanxian.org),將于2015年6月16—18日在北京國際會議中心舉辦?!氨本﹪H城市管線展覽會—CIUP EXPO”立足中國管線產(chǎn)業(yè)發(fā)展,注重管線全產(chǎn)業(yè)鏈建設(shè),倡導(dǎo)管線市場健康繁榮,推動(dòng)行業(yè)自主技術(shù)創(chuàng)新和產(chǎn)業(yè)化。

    2014年11月14—16日曾在北京召開了首屆北京國際城市管線大會,并同期召開了首屆北京國際城市管線展覽會,會議主題為“智慧管線與城市安全”。在首屆大會上,一大批管線先進(jìn)企業(yè)集體亮相于展會,把新技術(shù)、新產(chǎn)品統(tǒng)一展示給中外賓客,全方位展示了我國管線行業(yè)國際化視野和戰(zhàn)略成果,打造管線行業(yè)全產(chǎn)業(yè)鏈展示與交流平臺。首屆展會近150家參展廠商,3天共接待專業(yè)觀眾1600多人,展會期間共舉辦了11場主題和專題報(bào)告會。

    2015城市管線大會將以“共建、共享=互聯(lián)網(wǎng)+時(shí)代下的智慧管網(wǎng)”為主題,全面展示我國管線行業(yè)的最新發(fā)展成果,涵蓋給水、排水、燃?xì)?、熱力、石油、石化、電力、通信管線以及地下工程等領(lǐng)域的城市管線綜合解決方案。而同期舉辦的城市管線展覽會,還將舉辦多場高水平、高層次、高規(guī)格的產(chǎn)業(yè)論壇,分別涵蓋給排水、燃?xì)?、石油石化、智慧生態(tài)、安全管理、工程施工等主題。

    (本刊編輯部)

    极品教师在线免费播放| 精品人妻1区二区| 国产三级中文精品| 久久中文看片网| 亚洲av第一区精品v没综合| 一本一本综合久久| 欧美一级a爱片免费观看看| 欧美三级亚洲精品| 他把我摸到了高潮在线观看| 亚洲精华国产精华精| 亚洲欧美精品综合久久99| 赤兔流量卡办理| 不卡一级毛片| 国产成人福利小说| 给我免费播放毛片高清在线观看| 久久国产精品人妻蜜桃| 人妻夜夜爽99麻豆av| 久久婷婷人人爽人人干人人爱| 欧美日韩乱码在线| 国产亚洲欧美98| 国产在线男女| 看十八女毛片水多多多| 特级一级黄色大片| 国产午夜精品论理片| 国产精品一区二区免费欧美| 国产老妇女一区| 中国美女看黄片| 亚洲av免费高清在线观看| 亚洲最大成人手机在线| 亚洲人成网站在线播| 国产69精品久久久久777片| 91麻豆av在线| 亚洲不卡免费看| 国产伦精品一区二区三区四那| 亚洲成av人片在线播放无| 成人av一区二区三区在线看| 琪琪午夜伦伦电影理论片6080| 亚洲最大成人手机在线| 国产激情偷乱视频一区二区| 欧美又色又爽又黄视频| 免费大片18禁| 性欧美人与动物交配| 欧美高清成人免费视频www| 国产高清视频在线播放一区| 色播亚洲综合网| 岛国在线免费视频观看| 中文字幕免费在线视频6| 天天一区二区日本电影三级| 亚洲精品成人久久久久久| 亚洲av一区综合| 美女高潮的动态| 国产亚洲精品av在线| 99国产极品粉嫩在线观看| 成年女人永久免费观看视频| 99热这里只有是精品在线观看| 亚洲午夜理论影院| 久久国产乱子免费精品| 精品久久久久久久末码| 乱码一卡2卡4卡精品| 免费av不卡在线播放| 男人舔女人下体高潮全视频| 精品乱码久久久久久99久播| 亚洲成人精品中文字幕电影| 亚洲不卡免费看| 亚洲精品456在线播放app | 精品人妻偷拍中文字幕| 天堂√8在线中文| 国产精品女同一区二区软件 | 久久精品国产亚洲av香蕉五月| av视频在线观看入口| 久久久国产成人精品二区| 久久久久久久久久成人| 亚洲乱码一区二区免费版| 乱人视频在线观看| or卡值多少钱| 哪里可以看免费的av片| 亚洲欧美激情综合另类| 欧美绝顶高潮抽搐喷水| 日韩精品有码人妻一区| 日韩高清综合在线| 欧美一级a爱片免费观看看| 麻豆久久精品国产亚洲av| 99久久精品国产国产毛片| xxxwww97欧美| 网址你懂的国产日韩在线| 久久亚洲精品不卡| 亚洲四区av| 亚洲在线观看片| av在线蜜桃| 少妇的逼好多水| 狂野欧美白嫩少妇大欣赏| 波多野结衣巨乳人妻| 久久精品人妻少妇| 一级毛片久久久久久久久女| 精品99又大又爽又粗少妇毛片 | 嫁个100分男人电影在线观看| 欧美一区二区精品小视频在线| 国产精品人妻久久久影院| 日本欧美国产在线视频| 欧美高清成人免费视频www| 亚洲自拍偷在线| 69av精品久久久久久| 欧美高清性xxxxhd video| 精品福利观看| 久久人人爽人人爽人人片va| 成年免费大片在线观看| 婷婷亚洲欧美| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 九色国产91popny在线| 久久久久免费精品人妻一区二区| 欧美潮喷喷水| 91av网一区二区| 精品久久久久久成人av| 老司机深夜福利视频在线观看| 少妇人妻一区二区三区视频| 18禁裸乳无遮挡免费网站照片| 日本一本二区三区精品| 欧美日本亚洲视频在线播放| 啪啪无遮挡十八禁网站| 高清在线国产一区| 12—13女人毛片做爰片一| 国产一区二区三区在线臀色熟女| 亚洲欧美精品综合久久99| 欧美最新免费一区二区三区| 亚洲内射少妇av| 日日摸夜夜添夜夜添av毛片 | 欧美日韩亚洲国产一区二区在线观看| 亚洲国产高清在线一区二区三| 日韩欧美一区二区三区在线观看| 在线观看美女被高潮喷水网站| 成人国产麻豆网| 日韩大尺度精品在线看网址| 国产av麻豆久久久久久久| 国产免费男女视频| 成年女人看的毛片在线观看| 无遮挡黄片免费观看| 欧美成人性av电影在线观看| 亚洲成a人片在线一区二区| 国内久久婷婷六月综合欲色啪| 免费搜索国产男女视频| 欧美一区二区精品小视频在线| av在线蜜桃| 午夜福利高清视频| 日韩欧美精品v在线| 日日撸夜夜添| 观看美女的网站| 一本一本综合久久| 直男gayav资源| 亚洲av.av天堂| av福利片在线观看| 亚洲精品亚洲一区二区| 国产午夜精品久久久久久一区二区三区 | 日本黄色视频三级网站网址| 国产精品无大码| 有码 亚洲区| www.www免费av| 亚洲精品乱码久久久v下载方式| 中文字幕av在线有码专区| 亚洲中文字幕日韩| 黄片wwwwww| 国国产精品蜜臀av免费| 久久午夜福利片| 国产精品伦人一区二区| 亚洲人成网站高清观看| 日本-黄色视频高清免费观看| 看十八女毛片水多多多| .国产精品久久| 内地一区二区视频在线| av专区在线播放| 变态另类成人亚洲欧美熟女| 国产乱人视频| 高清日韩中文字幕在线| 91在线精品国自产拍蜜月| 少妇裸体淫交视频免费看高清| 搡老妇女老女人老熟妇| 18禁黄网站禁片午夜丰满| 国产综合懂色| 亚洲va日本ⅴa欧美va伊人久久| 婷婷精品国产亚洲av在线| 哪里可以看免费的av片| av视频在线观看入口| 国产aⅴ精品一区二区三区波| 亚洲成人久久性| 国产精品嫩草影院av在线观看 | 人人妻人人澡欧美一区二区| 亚洲欧美清纯卡通| 小蜜桃在线观看免费完整版高清| 国产精品一及| 久久精品国产99精品国产亚洲性色| 国产三级中文精品| videossex国产| 大又大粗又爽又黄少妇毛片口| 人人妻人人澡欧美一区二区| 国产淫片久久久久久久久| 国产伦在线观看视频一区| 久久精品国产亚洲av香蕉五月| 国产精品美女特级片免费视频播放器| 免费无遮挡裸体视频| 亚洲国产欧洲综合997久久,| 亚洲国产精品久久男人天堂| 久久亚洲精品不卡| 精品99又大又爽又粗少妇毛片 | 男插女下体视频免费在线播放| netflix在线观看网站| 精品国产三级普通话版| 最新在线观看一区二区三区| 99视频精品全部免费 在线| 91久久精品国产一区二区三区| 国产久久久一区二区三区| 99久久久亚洲精品蜜臀av| 亚洲aⅴ乱码一区二区在线播放| 欧美一级a爱片免费观看看| 国产 一区精品| 毛片女人毛片| 日韩精品中文字幕看吧| 国产精品三级大全| 国产高清三级在线| 成人国产麻豆网| 国产美女午夜福利| 又紧又爽又黄一区二区| 淫秽高清视频在线观看| 成人av一区二区三区在线看| 久久久久久久久久成人| 免费在线观看影片大全网站| 午夜福利18| 亚洲国产欧美人成| 日韩欧美免费精品| 高清毛片免费观看视频网站| 美女高潮的动态| 精品久久久久久久人妻蜜臀av| 日本成人三级电影网站| 国产不卡一卡二| 国产伦人伦偷精品视频| 午夜激情欧美在线| 精品国内亚洲2022精品成人| 黄色欧美视频在线观看| 国产一区二区在线观看日韩| 亚洲精品久久国产高清桃花| 国产探花极品一区二区| 床上黄色一级片| 亚洲avbb在线观看| 久久久午夜欧美精品| 一个人观看的视频www高清免费观看| 99热精品在线国产| 亚洲av美国av| 91麻豆av在线| 日韩欧美国产一区二区入口| 一本一本综合久久| 国产亚洲精品久久久久久毛片| 欧美性感艳星| 日日摸夜夜添夜夜添小说| 午夜a级毛片| 国产精品美女特级片免费视频播放器| 永久网站在线| 欧美zozozo另类| 国内精品久久久久精免费| 日韩一本色道免费dvd| 偷拍熟女少妇极品色| 国内精品久久久久精免费| 久久精品夜夜夜夜夜久久蜜豆| 麻豆成人午夜福利视频| 亚洲人成网站在线播放欧美日韩| 欧美成人性av电影在线观看| 国产精品99久久久久久久久| 久久久久久大精品| 亚洲三级黄色毛片| avwww免费| 国产中年淑女户外野战色| 久久久久久伊人网av| 亚洲欧美日韩高清专用| 国产精品亚洲一级av第二区| 欧美zozozo另类| 日韩在线高清观看一区二区三区 | 精品久久久久久久末码| 少妇高潮的动态图| 日本免费a在线| 国产探花在线观看一区二区| 日韩欧美一区二区三区在线观看| 亚洲人成网站在线播| 国产精品福利在线免费观看| 国产成人av教育| 中文亚洲av片在线观看爽| 日日干狠狠操夜夜爽| 精品久久久久久久人妻蜜臀av| 亚洲av美国av| 精品无人区乱码1区二区| 九九爱精品视频在线观看| 人人妻人人澡欧美一区二区| 非洲黑人性xxxx精品又粗又长| 久久精品国产自在天天线| 国产麻豆成人av免费视频| 亚洲国产精品合色在线| 国产乱人视频| 日韩一区二区视频免费看| 亚洲国产高清在线一区二区三| АⅤ资源中文在线天堂| 国产亚洲精品久久久com| 人妻丰满熟妇av一区二区三区| 亚洲成人久久性| 日韩欧美免费精品| bbb黄色大片| 久久国产精品人妻蜜桃| 欧美一区二区精品小视频在线| 俄罗斯特黄特色一大片| 国产久久久一区二区三区| 亚洲一区二区三区色噜噜| a级毛片免费高清观看在线播放| 91精品国产九色| 亚洲中文字幕日韩| 九九热线精品视视频播放| 91精品国产九色| 欧美日韩黄片免| 一个人免费在线观看电影| 亚洲精品456在线播放app | 久久久久国内视频| 又黄又爽又免费观看的视频| 俄罗斯特黄特色一大片| av专区在线播放| 免费一级毛片在线播放高清视频| 床上黄色一级片| 禁无遮挡网站| 淫妇啪啪啪对白视频| 成人毛片a级毛片在线播放| 变态另类丝袜制服| 免费观看精品视频网站| 变态另类丝袜制服| 窝窝影院91人妻| 亚州av有码| 国内久久婷婷六月综合欲色啪| 五月伊人婷婷丁香| 国产精品av视频在线免费观看| 日韩强制内射视频| 日韩欧美国产在线观看| 色吧在线观看| 国产大屁股一区二区在线视频| 一a级毛片在线观看| 夜夜夜夜夜久久久久| 丰满乱子伦码专区| 日韩大尺度精品在线看网址| 国产亚洲91精品色在线| 可以在线观看毛片的网站| 亚州av有码| 国产综合懂色| 天堂动漫精品| 十八禁国产超污无遮挡网站| 欧美黑人欧美精品刺激| 国产aⅴ精品一区二区三区波| 国产精品久久久久久精品电影| 亚洲avbb在线观看| 黄色女人牲交| 亚洲第一电影网av| 一进一出好大好爽视频| 国产高清三级在线| 国产国拍精品亚洲av在线观看| 亚洲在线自拍视频| 亚洲精品粉嫩美女一区| 亚洲中文日韩欧美视频| 极品教师在线免费播放| 伊人久久精品亚洲午夜| 日本a在线网址| 欧美+日韩+精品| 中文字幕精品亚洲无线码一区| 亚洲va日本ⅴa欧美va伊人久久| 日本免费a在线| 国产伦在线观看视频一区| 偷拍熟女少妇极品色| 韩国av一区二区三区四区| 丰满乱子伦码专区| 国产午夜精品论理片| 免费大片18禁| 又黄又爽又免费观看的视频| 色精品久久人妻99蜜桃| 又紧又爽又黄一区二区| 午夜精品在线福利| 亚洲中文字幕日韩| 欧美激情在线99| av在线天堂中文字幕| 黄色日韩在线| 婷婷色综合大香蕉| 久久午夜福利片| 午夜精品一区二区三区免费看| av在线观看视频网站免费| 国产精品久久视频播放| 欧美高清成人免费视频www| 高清日韩中文字幕在线| 免费看av在线观看网站| 亚洲电影在线观看av| 最近中文字幕高清免费大全6 | 人人妻人人看人人澡| 日韩欧美一区二区三区在线观看| 韩国av在线不卡| 九九在线视频观看精品| 少妇裸体淫交视频免费看高清| 精品国内亚洲2022精品成人| 国产蜜桃级精品一区二区三区| 亚洲最大成人中文| 久久精品91蜜桃| 男人和女人高潮做爰伦理| 国产一区二区三区视频了| 好男人在线观看高清免费视频| 国产免费一级a男人的天堂| 久久午夜亚洲精品久久| 欧美一区二区精品小视频在线| 国产成人av教育| 亚洲真实伦在线观看| 成人二区视频| 日韩精品有码人妻一区| 又黄又爽又刺激的免费视频.| 亚洲国产欧洲综合997久久,| 午夜免费男女啪啪视频观看 | 少妇的逼水好多| 亚洲狠狠婷婷综合久久图片| 中国美白少妇内射xxxbb| 高清日韩中文字幕在线| 日韩av在线大香蕉| 国产高清视频在线观看网站| 亚洲欧美激情综合另类| 日韩欧美在线乱码| 久久久久久久久久成人| 久久精品国产自在天天线| www.www免费av| 欧美精品啪啪一区二区三区| 乱人视频在线观看| 欧美最黄视频在线播放免费| 午夜福利欧美成人| 国产精品一区www在线观看 | 久久久久久久午夜电影| 给我免费播放毛片高清在线观看| 51国产日韩欧美| 国产精品久久久久久久电影| 亚洲av第一区精品v没综合| 小蜜桃在线观看免费完整版高清| 久久久久久久久久黄片| xxxwww97欧美| 成人高潮视频无遮挡免费网站| 综合色av麻豆| 97超级碰碰碰精品色视频在线观看| av天堂中文字幕网| 国产高清激情床上av| 国产精品女同一区二区软件 | 91午夜精品亚洲一区二区三区 | 亚洲国产高清在线一区二区三| 欧美zozozo另类| 精品不卡国产一区二区三区| 综合色av麻豆| 欧美潮喷喷水| 制服丝袜大香蕉在线| 看十八女毛片水多多多| 好男人在线观看高清免费视频| 男人舔奶头视频| 亚洲一区高清亚洲精品| 精品国产三级普通话版| 中文字幕久久专区| 舔av片在线| 天天一区二区日本电影三级| 国产成人影院久久av| ponron亚洲| 亚洲成a人片在线一区二区| 色噜噜av男人的天堂激情| 又爽又黄a免费视频| 少妇的逼好多水| 久久久久免费精品人妻一区二区| 日韩av在线大香蕉| 精品一区二区三区人妻视频| 免费在线观看成人毛片| 久久精品夜夜夜夜夜久久蜜豆| av在线老鸭窝| 听说在线观看完整版免费高清| 欧美人与善性xxx| 91久久精品国产一区二区成人| 精品不卡国产一区二区三区| 51国产日韩欧美| 日韩欧美在线二视频| 国产午夜精品论理片| 国产精品久久久久久av不卡| 午夜福利在线观看免费完整高清在 | 久久精品国产亚洲av香蕉五月| 淫妇啪啪啪对白视频| 亚洲最大成人av| 欧美不卡视频在线免费观看| 在线免费观看不下载黄p国产 | 国产aⅴ精品一区二区三区波| 五月玫瑰六月丁香| 18+在线观看网站| 女人十人毛片免费观看3o分钟| 欧美区成人在线视频| 欧美色视频一区免费| 国产三级在线视频| 亚洲美女搞黄在线观看 | 能在线免费观看的黄片| 成人国产麻豆网| 欧美一级a爱片免费观看看| 国产精品精品国产色婷婷| 看免费成人av毛片| 亚洲欧美日韩高清专用| 精品人妻视频免费看| 中文在线观看免费www的网站| 久久精品综合一区二区三区| 日本熟妇午夜| 日韩亚洲欧美综合| 久久国产乱子免费精品| 精品一区二区三区av网在线观看| 一级毛片久久久久久久久女| 中文亚洲av片在线观看爽| x7x7x7水蜜桃| 丰满的人妻完整版| 午夜日韩欧美国产| 高清在线国产一区| 亚洲av日韩精品久久久久久密| 亚洲欧美日韩东京热| 毛片女人毛片| 亚洲欧美日韩无卡精品| 99热这里只有是精品50| 国产精品女同一区二区软件 | 国产成人aa在线观看| 黄片wwwwww| 啦啦啦观看免费观看视频高清| 日本免费a在线| 亚洲电影在线观看av| 色视频www国产| 国产精品,欧美在线| 在线国产一区二区在线| 亚洲精品影视一区二区三区av| 一进一出抽搐gif免费好疼| 久久久国产成人精品二区| 亚洲av美国av| 国产探花极品一区二区| 日韩亚洲欧美综合| 国产一区二区在线观看日韩| 成年女人毛片免费观看观看9| 亚洲成人免费电影在线观看| 日本撒尿小便嘘嘘汇集6| 久久婷婷人人爽人人干人人爱| 日韩人妻高清精品专区| 久久精品国产亚洲网站| 国产精品一及| 成人永久免费在线观看视频| a在线观看视频网站| 国产三级在线视频| 国产精品一区二区免费欧美| 嫩草影院新地址| 免费观看在线日韩| 亚洲精品456在线播放app | netflix在线观看网站| 国产高清有码在线观看视频| 成人高潮视频无遮挡免费网站| 亚洲国产精品sss在线观看| 91在线精品国自产拍蜜月| 日韩国内少妇激情av| 女人十人毛片免费观看3o分钟| 国产精品98久久久久久宅男小说| 免费电影在线观看免费观看| 国产成人av教育| 国产精品一区二区性色av| 亚洲人与动物交配视频| 亚洲专区中文字幕在线| 久久精品国产亚洲网站| 美女大奶头视频| 亚洲avbb在线观看| 久久久久久伊人网av| 亚洲人成网站高清观看| 琪琪午夜伦伦电影理论片6080| 国产精品不卡视频一区二区| 日韩欧美免费精品| 欧美不卡视频在线免费观看| 日韩欧美精品免费久久| 久久精品91蜜桃| 乱系列少妇在线播放| 免费看a级黄色片| 91精品国产九色| 久久久成人免费电影| 国产高清三级在线| 女生性感内裤真人,穿戴方法视频| 日日夜夜操网爽| 深爱激情五月婷婷| 欧美精品啪啪一区二区三区| 久久精品国产亚洲网站| 免费在线观看影片大全网站| 欧洲精品卡2卡3卡4卡5卡区| 成人特级av手机在线观看| av天堂中文字幕网| 美女高潮的动态| 国产国拍精品亚洲av在线观看| 精品福利观看| 久久99热6这里只有精品| 国产男人的电影天堂91| 精品福利观看| 日本黄色片子视频| 国产亚洲91精品色在线| 亚洲av二区三区四区| 两人在一起打扑克的视频| 97超视频在线观看视频| 国产精品乱码一区二三区的特点| 乱系列少妇在线播放| 日日撸夜夜添| 国产高潮美女av| 在线播放国产精品三级| 少妇裸体淫交视频免费看高清| 国产高清视频在线播放一区| 黄色丝袜av网址大全| 嫩草影院精品99| 九九在线视频观看精品| 女的被弄到高潮叫床怎么办 | 淫妇啪啪啪对白视频| 亚洲国产精品成人综合色| 亚洲av中文av极速乱 | 国产一区二区三区在线臀色熟女| 亚洲第一区二区三区不卡| 亚洲久久久久久中文字幕| 欧美一区二区亚洲| 久久久久久九九精品二区国产| 国产av一区在线观看免费| 亚洲av第一区精品v没综合| 97碰自拍视频|