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

    應用iGMAS超快速星歷的實時精密單點定位研究

    2021-03-16 09:15:32徐偉證康國華彭攀周宏濤劉宗強趙騰
    中國空間科學技術(shù) 2021年1期
    關(guān)鍵詞:鐘差單點軌道

    徐偉證 ,康國華,*,彭攀 ,周宏濤,劉宗強,趙騰

    1. 南京航空航天大學 航天學院,南京 210016 2. 上海衛(wèi)星工程研究所,上海 201109

    為了進一步推動多模GNSS系統(tǒng)兼容和互操作,中國從2012年啟動國際GNSS監(jiān)測評估系統(tǒng)(International GNSS Monitoring & Assessment System, iGMAS)建設。iGMAS旨在建立我國自主的BDS、GPS、GLONASS和Galileo導航衛(wèi)星全弧段、多重覆蓋的全球近實時跟蹤網(wǎng),監(jiān)控GNSS導航衛(wèi)星的運行狀況、信號質(zhì)量和服務性能[1-2]。iGMAS能夠向全球用戶提供精密星歷、鐘差、地球定向參數(shù)等產(chǎn)品,為衛(wèi)星導航技術(shù)試驗提供支持,服務于科學研究以及各類應用,包括精密單點定位(Precious Point Positioning, PPP)技術(shù)[3]。

    相對于事后處理的精密單點定位,實時精準定位才能滿足5G、物聯(lián)網(wǎng)和無人駕駛等技術(shù)發(fā)展的需求,其中高質(zhì)量精密星歷、鐘差改正產(chǎn)品的實時獲取和應用是難點之一[4]。目前國際上開展實時PPP應用主要基于IGS實時、近實時產(chǎn)品,包括IGS實時定位服務(real-time service)、IGDG(Internet-Based Global Differential GPS)實時產(chǎn)品和BNC(BKG Ntrip Client)軟件等,以及西班牙GMV公司、Trimble公司等推出的商業(yè)實時PPP服務[5-7]。與國際上比較成熟的IGS系統(tǒng)相比,我國主導的iGMAS建設周期短,跟蹤站、分析中心和數(shù)據(jù)中心數(shù)量相對較少,因此iGMAS產(chǎn)品(特別是超快速產(chǎn)品)的實際應用仍需進一步發(fā)展。

    目前,國內(nèi)外應用iGMAS超快速產(chǎn)品開展的實時精密單點定位研究,主要依靠跟蹤站的觀測數(shù)據(jù)模擬實時PPP解算[6, 8, 9],未能解決iGMAS超快速數(shù)據(jù)的實時獲取和普遍實用問題。針對iGMAS產(chǎn)品的應用推廣和實時PPP對超快速精密數(shù)據(jù)的需求,本文對iGMAS超快速星歷產(chǎn)品精度和穩(wěn)定性進行分析,并通過高精度的IGS最終產(chǎn)品進行評估。本文以評估結(jié)果作為依據(jù)選擇精密改正數(shù),設計了iGMAS-Download實時下載程序和實時PPP算法,通過NovAtel雙頻接收機驗證了其在實時PPP中的性能,為iGMAS產(chǎn)品的應用提供借鑒思路。

    1 iGMAS超快速星歷質(zhì)量分析

    1.1 超快速星歷精度分析

    iGMAS和IGS衛(wèi)星精密軌道均包括三種不同精度和時延的數(shù)據(jù),即超快速、快速和最終產(chǎn)品。其中超快速產(chǎn)品包含24 h觀測數(shù)據(jù)和24 h預測數(shù)據(jù),更新頻率為6 h,分別于當日03、09、15、21時左右發(fā)布,雖然存在約3h滯后時間,但基本可以滿足實時定位的需求[10]。

    IGS 已經(jīng)在全球建立507個跟蹤站,建有充足的分析中心、全球關(guān)聯(lián)分析中心、區(qū)域關(guān)聯(lián)分析中心及分析中心協(xié)調(diào)組織來協(xié)同處理數(shù)據(jù),其最終產(chǎn)品軌道精度約2.5 cm,鐘差精度約0.075 ns[4,8]。本文以精度較高的IGS MGEX(The Multi-GNSS Experiment)最終產(chǎn)品為參考標準,結(jié)合實時精密單點定位需求,選取iGMAS和IGS系統(tǒng)在2019年7月7日至13日(北斗周0705)的精密星歷數(shù)據(jù),計算一周時間內(nèi)iGMAS超快速星歷相對IGS最終產(chǎn)品參考真值的均方根誤差(Root Mean Square, RMS),評估其數(shù)據(jù)的精度和穩(wěn)定性。其中IGS最終產(chǎn)品從CDDIS(The Crustal Dynamics Data Information System)獲取,iGMAS超快速產(chǎn)品從武漢數(shù)據(jù)中心獲取[11, 12]。

    iGMAS超快速軌道、鐘差產(chǎn)品和IGS最終產(chǎn)品采樣間隔均為15 min,但IGS產(chǎn)品采用GPS時間系統(tǒng)(GPST),iGMAS產(chǎn)品采用BDS時間系統(tǒng)(BDT),進行質(zhì)量評估需進行時間基準的協(xié)調(diào)。本研究在精密軌道和鐘差計算時以GPST作為時間基準,BDT與GPST的轉(zhuǎn)換關(guān)系為[13]:

    GPST=BDT+14+τUTC(USUN)-UTC(NTSC)

    (1)

    其中τUTC(USUN)-UTC(NTSC)表示GPS時和北斗時維持的協(xié)調(diào)世界時之間的微小量差值。從精度、階數(shù)和歷元間隔三方面衡量,9階擬合多項式求解GPS衛(wèi)星坐標可達毫米級精度,且效果最好[14]。因此本文采用9階切比雪夫多項式擬合來計算軌道數(shù)據(jù),選取北斗周0705期間內(nèi)iGMAS超快速星歷與IGS最終星歷。

    本文計算iGMAS超快速軌道與IGS最終軌道RMS差值,計算方式如下式所示[15]:

    分別對iGMAS超快速星歷和IGS最終星歷合對GPS、GLONASS、BDS和Galileo系統(tǒng)在北斗周0705期間的軌道數(shù)據(jù),各顆衛(wèi)星軌道誤差RMS和系統(tǒng)整體誤差結(jié)果如圖 1所示。

    1)GPS/GLONASS/Galileo/BDS(MEO/IGSO)/BDS(GEO)和系統(tǒng)在北斗周0705期間,系統(tǒng)軌道誤差均值分別為5.7 cm、9.8 cm、14.2 cm、18.0 cm、393.7 cm;BDS(GEO)衛(wèi)星相對地面站空間幾何構(gòu)型較差,參數(shù)相關(guān)性強而很難精確求解,故精度較低。

    2)由圖1(a)~(e)可見,各系統(tǒng)衛(wèi)星的軌道誤差RMS分布較為均勻,系統(tǒng)軌道精度依次是GPS、GLONASS、Galileo、BDS(MEO/IGSO)、BDS(GEO); GPS衛(wèi)星軌道精度穩(wěn)定在厘米級,且誤差RMS更為集中,主要得益于目前IGS地面跟蹤站較多,更為豐富的觀測數(shù)據(jù)有助于減小隨機誤差。

    3)由圖1(b)(e)可知,iGMAS數(shù)據(jù)目前提供的北斗衛(wèi)星精密數(shù)據(jù)遠少于實際在軌北斗衛(wèi)星數(shù),且系統(tǒng)間衛(wèi)星軌道精度較為離散。主要原因是北斗導航系統(tǒng)近兩年衛(wèi)星部署較多,在軌調(diào)試需要時間,而iGMAS跟蹤站數(shù)量少且分布集中,對北斗導航系統(tǒng)的監(jiān)測質(zhì)量和完善仍有待提升。

    圖1 iGMAS超快速星歷與IGS星歷各系統(tǒng)軌道差值Fig.1 Four navigation systems orbit error between iGMAS ultra-fast ephemeris and IGS final product

    通過對四大導航系統(tǒng)iGMAS超快速星歷的評估,可見GPS衛(wèi)星星歷精度在厘米級。衛(wèi)星鐘差同樣是影響實時精密單點定位的重要因素,因此本文針對GPS系統(tǒng)的衛(wèi)星軌道和鐘差做了進一步誤差分析,各顆衛(wèi)星在ECEF坐標系下X、Y、Z和3D方向誤差如圖 2所示,各顆衛(wèi)星和系統(tǒng)均值的鐘差如圖 3所示。在評估周期內(nèi),各顆衛(wèi)星的X、Y、Z方向的軌道誤差RMS均在2.0~5.5 cm范圍內(nèi),且Z方向軌道誤差略優(yōu)于X、Y方向,各顆衛(wèi)星3D誤差均值為5.78 cm。各顆GPS衛(wèi)星的鐘差集中在2.94~3.1 ns區(qū)間內(nèi),系統(tǒng)均值優(yōu)于3 ns,需要基于超快速軌道對鐘差進一步估計才能滿足厘米級定位需求。

    圖2 iGMAS超快速與IGS最終產(chǎn)品GPS衛(wèi)星差值RMSFig.2 GPS orbit error between iGMAS Ultra-fast and IGS final product

    圖3 iGMAS超快速與IGS最終產(chǎn)品的GPS衛(wèi)星鐘差Fig. 3 GPS satellite clock error between iGMAS Ultra-fast and IGS final product

    1.2 超快速星歷穩(wěn)定性分析

    實時PPP應用性能與精密星歷的穩(wěn)定性緊密相關(guān)。本文主要從系統(tǒng)軌道精度在一定周期內(nèi)的變化趨勢,以及超快速星歷觀測和預報部分的精度變化趨勢來評估其數(shù)據(jù)穩(wěn)定性。

    在北斗周0705期間四大導航系統(tǒng)每天衛(wèi)星及星座軌道誤差RMS,在擬合期間變化趨勢如圖 4所示,圖 5表示GPS衛(wèi)星軌道變化。

    1)由圖4可見在擬合周期內(nèi)GPS星座系統(tǒng)穩(wěn)定性最優(yōu),軌道整體誤差和仿真周期內(nèi)波動均相對最??;GLONASS穩(wěn)定性次之,存在小范圍波動,應與當日觀測數(shù)據(jù)質(zhì)量有關(guān);BDS(MEO/IGSO)以及Galileo導航系統(tǒng)相對波動較大。

    2)由圖4可見在擬合周期內(nèi),GPS、GLONASS和Galileo系統(tǒng)不同衛(wèi)星的軌道誤差RMS變化趨勢相同,軌道誤差同步增大或減??;BDS(MEO/IGSO)存在系統(tǒng)間定軌精度離散性相對較大的情況,更說明iGMAS需要建設全球分布監(jiān)測站或低軌星基監(jiān)測站,以提升北斗導航系統(tǒng)整體穩(wěn)定性。

    3)由圖5可見GPS衛(wèi)星在擬合周期內(nèi),所有衛(wèi)星的軌道誤差范圍均在3.4~8.2 cm,不同衛(wèi)星在同一天的軌道誤差差別較小,誤差范圍集中且系統(tǒng)整體定軌精度高;GPS系統(tǒng)有完善的全球分布的跟蹤站,可以實現(xiàn)衛(wèi)星的全弧段多重覆蓋。

    圖4 四系統(tǒng)軌道誤差在擬合周期內(nèi)的變化趨勢Fig.4 GPS、GLONASS、BDS(MEO/IGSO)、Galileo orbit error

    圖5 iGMAS超快速軌道相對IGS最終軌道差值變化Fig.5 Orbit difference changes between iGMAS ultrafast orbit and IGS final product(GPS)

    iGMAS超快速星歷在一周擬合周期內(nèi)GPS系統(tǒng)定軌精度在厘米級且最穩(wěn)定,目前依靠iGMAS超快速星歷開展實時精密單點定位試驗可優(yōu)先選取GPS衛(wèi)星觀測值參與定位解算,衛(wèi)星數(shù)量不足的情況下,可以用GLONASS系統(tǒng)聯(lián)立組合精密單點定位。

    iGMAS超快速星歷觀測和預測部分相對IGS最終產(chǎn)品的誤差分別如圖6(a)(b)所示。在擬合周期內(nèi),觀測部分星歷相對精度高且穩(wěn)定,除個別衛(wèi)星外,X、Y、Z方向誤差大致在±6 cm范圍內(nèi)。預測部分的軌道誤差相對較大,部分衛(wèi)星超過10 cm范圍,且隨預測時間越長軌道精度越低,因此實時PPP試驗需要實時獲取最新的iGMAS超快速星歷,以實現(xiàn)最優(yōu)精度的衛(wèi)星端誤差消除。

    圖6 iGMAS超快速星歷觀測、預測部分與IGS 最終精密星歷差值變化Fig.6 Error change of iGMAS ultra-fast ephemeris observation, prediction part and IGS final ephemeris

    2 基于iGMAS產(chǎn)品的實時 PPP定位

    2.1 實時PPP定位方案設計

    iGMAS超快速星歷的評估結(jié)果分析表明,目前GPS衛(wèi)星精度最高且穩(wěn)定性好,GLONASS次之,BDS(MEO/IGSO)和Galileo系統(tǒng)均存在一定范圍的波動,BDS(MEO/IGSO)各顆衛(wèi)星定軌精度較為離散。iGMAS超快速星歷觀測部分精度和穩(wěn)定性優(yōu)于預測部份,且隨預測時間越長精度越低,選取最近的超快速星歷有助于降低衛(wèi)星軌道和時鐘誤差。因此,本文設計iGAMS產(chǎn)品實時/事后下載軟件iGMAS-Download,實時獲取最新的超快速星歷,采用NovAtel雙頻接收機獲取GPS衛(wèi)星觀測數(shù)據(jù)開展實時精密單點定位應用研究,以驗證iGMAS超快速星歷在PPP試驗中的實用性能。

    本研究采用的實時PPP定位驗證方案如圖 7所示,通過自主設計的iGMAS-Download應用程序聯(lián)網(wǎng)實時獲取最新的iGMAS超快速星歷、鐘差。NovAtel OEM 617接收機雙頻偽距、載波相位觀測數(shù)據(jù)、廣播星歷和原始定位結(jié)果,實時經(jīng)串口轉(zhuǎn)USB通道輸入到計算機中,計算機進行觀測數(shù)據(jù)預處理、精密數(shù)據(jù)擬合以及實時PPP算法,實時輸出PPP定位結(jié)果。最后以北斗伴侶M2 RTK固定解為參考真值,結(jié)合接收機原始定位結(jié)果和實時PPP定位結(jié)果進行多通道定位結(jié)果的誤差對比分析,從而驗證iGMAS超快速星歷在實時PPP應用中的性能和可用性。

    圖7 實時PPP技術(shù)方案設計Fig.7 Real-time PPP technology solution

    2.2 iGMAS數(shù)據(jù)實時獲取和應用

    iGMAS已建成長沙、武漢、西安3個數(shù)據(jù)中心,用戶可通過任一數(shù)據(jù)中心或直接通過iGMAS官網(wǎng)(http:∥124.205.50.178)獲取相關(guān)產(chǎn)品。實時PPP技術(shù)應用需要獲取最新的超快速精密星歷產(chǎn)品,為了滿足本研究實時PPP算法的需求以及推廣iGMAS數(shù)據(jù)產(chǎn)品的實際應用,本文基于Windows平臺設計了實時/事后下載iGMAS精密數(shù)據(jù)的應用程序iGMAS-Download,程序下載文件的流程設計如圖 8所示。該程序包括實時和事后下載兩種模式,實時下載數(shù)據(jù)主要面向iGMAS超快速星歷,其文件命名格式為isuwwwwd_HH.sp3/clk.Z(wwww為北斗周;d為星期,0為星期日,1~6表示星期一至六;HH為小時,分為00,06,12,18;sp3/clk分別表示星歷和鐘差文件)。程序?qū)崟r刷新iGMAS數(shù)據(jù)中心產(chǎn)品列表,根據(jù)系統(tǒng)時間匹配最新的超快速產(chǎn)品文件,下載到指定文件夾并自動解壓,以備實時PPP算法讀取應用。事后模式可以根據(jù)用戶需求批量下載其他產(chǎn)品,包括衛(wèi)星軌道鐘差、跟蹤站地心坐標、地球自轉(zhuǎn)參數(shù)、大氣環(huán)境參數(shù)、頻間偏差信息、電離層閃爍指數(shù)、民用監(jiān)測評估結(jié)果、完好性產(chǎn)品等。

    圖8 iGMAS數(shù)據(jù)實時下載程序設計流程Fig.8 iGMAS data download program flowchart

    2.3 實時PPP算法設計

    本文采用PPP常規(guī)模型進行解算,采用雙頻偽距、載波相位觀測值的無電離層組合作為函數(shù)模型,其基礎模型表達式為[16-17]:

    (2)

    為保證高精度定位的質(zhì)量,本文首先剔除低高度角衛(wèi)星,并依靠高度角權(quán)重模型進行定權(quán),定權(quán)模型如式(3)所示[18]:

    δ2(e)=δ2/sin2(e)

    (3)

    式中:δ2是天頂方向的非差消電離層相位觀測值的方差;δ2(e)是在高度角e處的相位觀測值的方差。偽距觀測值也在得到天頂方向觀測值方差后,采用式(3)進行任意高度角映射,進而計算觀測值權(quán)陣[18]。對觀測數(shù)據(jù)進行實時預處理時,先剔除觀測值存在粗差的衛(wèi)星,再進行鐘跳探測與修復,避免將接收機鐘跳引起的觀測值跳變誤判為周跳。然后聯(lián)合使用寬巷Melbourne-Wubbena(MW)和Geometry-Free(GF)組合法進行周跳探測,對于出現(xiàn)周跳的衛(wèi)星進行模糊度初始化。通過數(shù)據(jù)的實時預處理,來消除異常觀測數(shù)據(jù)對PPP定位結(jié)果的影響[19-20]。

    本文采用擴展卡爾曼濾波進行參數(shù)估計,待估參數(shù)為X=(x,y,z,dtr,zpd,Ni(i=1~n))T,依次為接收機位置參數(shù)、接收機鐘差、天頂對流層濕延遲、可見衛(wèi)星模糊度參數(shù)。精密軌道和鐘差等產(chǎn)品采用iGMAS超快速產(chǎn)品,衛(wèi)星和接收機端天線相位中心偏差(PCO)和天線相位中心變化(PCV)分別使用IGS ANTEX文件進行改正。實時定位過程同步使用IGS超快速產(chǎn)品進行數(shù)據(jù)處理和結(jié)果對比。各項誤差采用模型進行處理,具體策略如表 1所示[9, 17]。

    表1 PPP各項誤差處理方法

    3 實時PPP試驗與分析

    3.1 實時PPP定位試驗環(huán)境搭建

    2019年7月28日14:00-17:30在南京航空航天大學明故宮校區(qū)操場開展了實測數(shù)據(jù)采集和靜態(tài)精密單點定位試驗。當日天氣狀況良好,實驗環(huán)境視野開闊無遮擋,數(shù)據(jù)采樣間隔為1秒。使用NovAtel OEM617雙頻接收機(單點定位精度:水平5 m,高程10 m)實時接收GPS觀測數(shù)據(jù),截止高度角設為15°,使用北斗伴侶M2 RTK固定解(定位精度:水平0.02 m,高程0.04 m)作為定位結(jié)果參考標準。在定位參考點的設計上,通過將NovAtel天線和北斗伴侶M2架設到兩臺三腳架并調(diào)整抵消其高度差,預先測量水平和高度誤差并在算法中予以修正。實驗過程同步使用兩臺筆記本電腦進行數(shù)據(jù)讀取和程序運算,以提高程序運算結(jié)果的可靠性。

    3.2 PPP試驗結(jié)果分析

    本文使用iGMAS超快速星歷和鐘差產(chǎn)品,對GPS實測數(shù)據(jù)進行實時精密單點定位實驗,其定位結(jié)果與NovAtel接收機自身定位誤差對比如圖 9所示,分別使用iGMAS和IGS超快速星歷的實時PPP定位結(jié)果相對于參考真值在ENU坐標系下收斂情況如圖 10、圖11所示。

    圖9 iGMAS超快速星歷實時PPP與接收機定位誤差對比Fig.9 Comparison of positioning error between iGMAS ultra-fast ephemeris real-time PPP and receiver

    圖10 iGMAS超快速星歷實時PPP定位誤差收斂情況Fig.10 iGMAS ultra-fast ephemeris real-time PPP convergence

    圖11 IGS超快速星歷實時PPP定位誤差收斂情況Fig.11 IGS ultra-fast ephemeris real-time PPP positioning error convergence

    1)由圖 9可知,NovAtel接收機自身定位誤差波動大,且方程方向誤差較大;與接收機原始定位相比,應用iGMAS超快速星歷的實時精密單點定位可以有效提升定位精度和穩(wěn)定性。

    2)由圖 10和圖 11可見,分別使用iGMAS和IGS超快速星歷進行精密單點定位試驗,其收斂時間和定位精度接近,連續(xù)觀測約20 min后,PPP定位誤差即可收斂到分米級,最終可收斂到20 cm以內(nèi)。

    3)由圖 10和圖 11,以及表 2可知,iGMAS超快速星歷在實時PPP試驗中可實現(xiàn)與IGS超快速星歷相近的定位結(jié)果,應用iGMAS超快速星歷的實時精密單點定位可以有效提升定位精度和穩(wěn)定性,驗證了iGMAS超快速星歷的實用性能。

    表2 接收機和實時PPP定位誤差均方根統(tǒng)計

    4 結(jié)束語

    衛(wèi)星精密星歷是實現(xiàn)PPP的基礎數(shù)據(jù),其精度和可用性直接影響導航定位的性能。本文以IGS最終產(chǎn)品為參考標準,評估了iGMAS超快速星歷的精度和穩(wěn)定性,設計了iGMAS-Download應用程序?qū)崟r獲取最新的iGMAS超快速星歷,以其中精度和穩(wěn)定性較高的GPS衛(wèi)星數(shù)據(jù),結(jié)合NovAtel雙頻接收機開展了實時PPP定位算法研究和試驗,得到如下結(jié)論:

    1)iGAMS超快速星歷中,GPS衛(wèi)星精度和穩(wěn)定性最高,具有較高的可用性,GLONASS次之,BDS(MEO/IGSO、GEO)和Galileo系統(tǒng)均存在一定范圍的波動。GPS衛(wèi)星星歷得益于數(shù)量眾多且全球分布的IGS跟蹤站支持,可以通過全弧段多重覆蓋的觀測數(shù)據(jù)來提高精度。

    2)北斗導航系統(tǒng)精密改正數(shù)據(jù),存在衛(wèi)星數(shù)量不足、衛(wèi)星精度相對較低(特別是GEO衛(wèi)星),以及系統(tǒng)間衛(wèi)星軌道誤差離散性較大的問題。未來有望通過建設更多全球分布的地面監(jiān)測站或低軌星基監(jiān)測站,以及星間鏈路技術(shù)來提升北斗產(chǎn)品的精度。

    3)本文設計的iGMAS-Download軟件可以實時獲取最新iGMAS超快速星歷,基于該數(shù)據(jù)的GPS實時PPP算法在結(jié)合NovAtel雙頻接收機試驗中,與應用IGS超快速星歷定位精度和收斂時間接近,驗證了iGMAS超快速產(chǎn)品的應用性能,有助于推動我國iGMAS產(chǎn)品的實際應用。

    猜你喜歡
    鐘差單點軌道
    歷元間載波相位差分的GPS/BDS精密單點測速算法
    基于單純形法的TLE軌道確定
    CryoSat提升軌道高度與ICESat-2同步運行
    朝美重回“相互羞辱軌道”?
    超薄異型坯連鑄機非平衡單點澆鑄實踐與分析
    山東冶金(2019年5期)2019-11-16 09:09:10
    IGS快速/超快速衛(wèi)星鐘差精度評定與分析
    實時干涉測量中對流層延遲與鐘差精修正建模
    載人航天(2016年4期)2016-12-01 06:56:24
    數(shù)字電視地面?zhèn)鬏斢脝晤l網(wǎng)與單點發(fā)射的效果比較
    基于拉格朗日的IGS精密星歷和鐘差插值分析
    16噸單點懸掛平衡軸的優(yōu)化設計
    人体艺术视频欧美日本| 国产一区亚洲一区在线观看| 91精品国产国语对白视频| 亚洲欧美日韩另类电影网站 | 久久久久久久久久久免费av| 国产真实伦视频高清在线观看| 黑人高潮一二区| 久久青草综合色| 国产乱来视频区| 国产淫语在线视频| 日韩人妻高清精品专区| 2018国产大陆天天弄谢| 成年人午夜在线观看视频| 中文字幕制服av| 色5月婷婷丁香| 日韩av在线免费看完整版不卡| 精品午夜福利在线看| 亚洲欧美日韩卡通动漫| 91精品国产九色| 亚洲国产精品999| 这个男人来自地球电影免费观看 | 亚洲精品自拍成人| 欧美精品国产亚洲| 春色校园在线视频观看| 欧美一级a爱片免费观看看| 老司机影院毛片| 男女下面进入的视频免费午夜| 人人妻人人添人人爽欧美一区卜 | 熟女电影av网| 各种免费的搞黄视频| 午夜老司机福利剧场| 九九爱精品视频在线观看| 内射极品少妇av片p| 国产极品天堂在线| 国产亚洲最大av| 国产美女午夜福利| 日日摸夜夜添夜夜添av毛片| 亚洲欧美日韩无卡精品| 国精品久久久久久国模美| 99热6这里只有精品| 青春草视频在线免费观看| 一个人看的www免费观看视频| 欧美日韩精品成人综合77777| 美女xxoo啪啪120秒动态图| 最近2019中文字幕mv第一页| 亚洲国产精品999| 国产极品天堂在线| av国产精品久久久久影院| 久久精品久久久久久久性| 人人妻人人看人人澡| 日韩av不卡免费在线播放| 一级a做视频免费观看| 国产片特级美女逼逼视频| 亚洲人成网站高清观看| av在线蜜桃| xxx大片免费视频| 99久久精品国产国产毛片| tube8黄色片| 亚洲电影在线观看av| 搡女人真爽免费视频火全软件| 国产成人免费无遮挡视频| 美女视频免费永久观看网站| 人人妻人人澡人人爽人人夜夜| 搡女人真爽免费视频火全软件| 91久久精品国产一区二区成人| 欧美日韩亚洲高清精品| 九草在线视频观看| 日韩精品有码人妻一区| 成人无遮挡网站| 久久精品国产亚洲av天美| 女性生殖器流出的白浆| 最近最新中文字幕免费大全7| 一本久久精品| 日韩视频在线欧美| 久久99精品国语久久久| 国产免费视频播放在线视频| 欧美日韩亚洲高清精品| 国产男女超爽视频在线观看| 国产一区亚洲一区在线观看| 韩国av在线不卡| 欧美高清成人免费视频www| 久久精品人妻少妇| 久久精品久久精品一区二区三区| 中文字幕免费在线视频6| 国产精品99久久久久久久久| 精品熟女少妇av免费看| 欧美bdsm另类| 亚洲欧美日韩卡通动漫| 又黄又爽又刺激的免费视频.| 全区人妻精品视频| 精品久久国产蜜桃| 国产无遮挡羞羞视频在线观看| 欧美激情国产日韩精品一区| 日日撸夜夜添| 少妇 在线观看| 久久亚洲国产成人精品v| 日本与韩国留学比较| 午夜免费鲁丝| 99久久人妻综合| .国产精品久久| 国产成人免费观看mmmm| 下体分泌物呈黄色| 国产成人午夜福利电影在线观看| 街头女战士在线观看网站| 深爱激情五月婷婷| 亚洲精品乱码久久久v下载方式| .国产精品久久| 免费不卡的大黄色大毛片视频在线观看| 一本一本综合久久| 在线观看一区二区三区激情| 亚洲真实伦在线观看| 日本av手机在线免费观看| 蜜桃久久精品国产亚洲av| 色网站视频免费| 亚洲国产欧美在线一区| 午夜激情福利司机影院| 国产成人精品福利久久| 亚洲av成人精品一二三区| 日韩免费高清中文字幕av| 日韩视频在线欧美| 亚洲精品视频女| 一级二级三级毛片免费看| 免费看日本二区| 免费av不卡在线播放| 日韩大片免费观看网站| 免费在线观看成人毛片| 精品酒店卫生间| 亚洲国产精品专区欧美| 国产亚洲精品久久久com| 久久婷婷青草| 亚洲美女黄色视频免费看| 国精品久久久久久国模美| 国产乱人偷精品视频| 欧美xxⅹ黑人| 观看美女的网站| 成年av动漫网址| 国产精品久久久久久久电影| 国产男女内射视频| 亚洲国产成人一精品久久久| 天堂中文最新版在线下载| 夫妻午夜视频| 亚洲精品国产av蜜桃| 精品国产一区二区三区久久久樱花 | 亚洲成人中文字幕在线播放| 一本一本综合久久| 看免费成人av毛片| 国产成人精品婷婷| 欧美成人精品欧美一级黄| 国产av精品麻豆| 欧美精品国产亚洲| 纯流量卡能插随身wifi吗| 老司机影院毛片| 一边亲一边摸免费视频| 欧美 日韩 精品 国产| 午夜福利影视在线免费观看| 蜜桃在线观看..| 亚洲熟女精品中文字幕| 毛片女人毛片| 免费看光身美女| 少妇人妻久久综合中文| 免费黄色在线免费观看| 成人午夜精彩视频在线观看| 最新中文字幕久久久久| 国产精品久久久久久精品电影小说 | 久久久成人免费电影| 成年人午夜在线观看视频| 色5月婷婷丁香| 亚洲av成人精品一二三区| 日韩亚洲欧美综合| 亚洲精华国产精华液的使用体验| 亚洲伊人久久精品综合| 国产精品久久久久久精品电影小说 | 久久久久视频综合| 亚洲国产色片| 国产黄片美女视频| 一级a做视频免费观看| 国产一级毛片在线| av网站免费在线观看视频| 日本av手机在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美 国产精品| 99九九线精品视频在线观看视频| 在线观看人妻少妇| 日本色播在线视频| 国产精品一区www在线观看| 十八禁网站网址无遮挡 | videossex国产| 久久这里有精品视频免费| 婷婷色综合大香蕉| 成年人午夜在线观看视频| 久久久久久久久久成人| 最近最新中文字幕免费大全7| 少妇 在线观看| 国产白丝娇喘喷水9色精品| 国产男人的电影天堂91| 久久久久久久国产电影| 亚洲精品aⅴ在线观看| 日本av手机在线免费观看| 免费黄频网站在线观看国产| 精品久久久久久久久亚洲| 美女cb高潮喷水在线观看| 国产伦理片在线播放av一区| xxx大片免费视频| 欧美3d第一页| 交换朋友夫妻互换小说| 丰满少妇做爰视频| 日韩,欧美,国产一区二区三区| 大码成人一级视频| 熟妇人妻不卡中文字幕| 国产精品福利在线免费观看| 美女cb高潮喷水在线观看| 日日摸夜夜添夜夜添av毛片| 99久久中文字幕三级久久日本| 人妻系列 视频| 精品国产一区二区三区久久久樱花 | 精品一区二区三区视频在线| av福利片在线观看| 亚洲精品国产色婷婷电影| 久久久久久久久久久丰满| 国产乱人视频| 性色avwww在线观看| 欧美高清成人免费视频www| 日本vs欧美在线观看视频 | 最近最新中文字幕大全电影3| 亚洲精品日韩在线中文字幕| 夜夜骑夜夜射夜夜干| 国产成人午夜福利电影在线观看| 国产一区亚洲一区在线观看| 99久久人妻综合| 婷婷色麻豆天堂久久| 欧美成人一区二区免费高清观看| 日韩在线高清观看一区二区三区| 国产精品人妻久久久影院| 性色av一级| 国产精品久久久久久久电影| 在线观看免费高清a一片| 亚洲丝袜综合中文字幕| 一本色道久久久久久精品综合| 性色av一级| 午夜福利在线观看免费完整高清在| 午夜免费鲁丝| 亚洲精华国产精华液的使用体验| 精品午夜福利在线看| 18禁动态无遮挡网站| 亚洲av日韩在线播放| av在线播放精品| a级一级毛片免费在线观看| 少妇人妻久久综合中文| 纯流量卡能插随身wifi吗| 男女边摸边吃奶| 欧美日韩国产mv在线观看视频 | 好男人视频免费观看在线| 日韩成人伦理影院| 黑人猛操日本美女一级片| 日日啪夜夜爽| 看非洲黑人一级黄片| 精品人妻熟女av久视频| 国产精品久久久久久久电影| 亚洲高清免费不卡视频| 国产一区二区三区综合在线观看 | 国产精品麻豆人妻色哟哟久久| 久久99蜜桃精品久久| 久久久久网色| 少妇高潮的动态图| 超碰av人人做人人爽久久| 免费少妇av软件| 亚洲天堂av无毛| av女优亚洲男人天堂| 国产精品.久久久| 国产男女内射视频| 在线观看一区二区三区激情| 日韩精品有码人妻一区| 欧美日韩亚洲高清精品| 国产精品无大码| 亚洲国产色片| 黑丝袜美女国产一区| 天堂8中文在线网| 国产伦精品一区二区三区四那| www.色视频.com| 国产精品国产三级国产专区5o| 午夜福利在线观看免费完整高清在| 91久久精品国产一区二区成人| 亚洲欧美成人精品一区二区| 伦理电影免费视频| 久久精品熟女亚洲av麻豆精品| 国产成人91sexporn| 亚洲av综合色区一区| 寂寞人妻少妇视频99o| 嫩草影院入口| 伦理电影大哥的女人| 18禁在线播放成人免费| 国产又色又爽无遮挡免| 欧美精品国产亚洲| 亚洲精品一区蜜桃| 又大又黄又爽视频免费| 久久久久久久国产电影| 国产一区二区在线观看日韩| 国产精品久久久久久精品古装| 精品人妻偷拍中文字幕| 在线免费十八禁| 最近中文字幕2019免费版| 日韩成人av中文字幕在线观看| 国产午夜精品一二区理论片| a级毛片免费高清观看在线播放| 亚洲国产毛片av蜜桃av| 亚洲国产精品专区欧美| 伊人久久国产一区二区| 搡女人真爽免费视频火全软件| 亚洲成人一二三区av| 91午夜精品亚洲一区二区三区| 夜夜骑夜夜射夜夜干| 亚洲成色77777| av在线播放精品| 亚洲成人手机| 久久久久久久久久久免费av| 国产成人a∨麻豆精品| 熟女av电影| 97在线视频观看| 精品久久久噜噜| 亚洲综合色惰| 丝瓜视频免费看黄片| av免费观看日本| 国产精品一区www在线观看| 免费观看在线日韩| 国产亚洲精品久久久com| 国产精品一区二区在线观看99| 国产高清有码在线观看视频| 久久精品国产亚洲av涩爱| 中国国产av一级| 免费人成在线观看视频色| 亚洲国产精品专区欧美| 新久久久久国产一级毛片| 久久久久精品久久久久真实原创| 中文字幕免费在线视频6| 精品人妻视频免费看| 亚洲精品成人av观看孕妇| 边亲边吃奶的免费视频| 人妻少妇偷人精品九色| 久久精品人妻少妇| 午夜激情福利司机影院| 国产极品天堂在线| 国产成人aa在线观看| 国产免费一级a男人的天堂| 我要看黄色一级片免费的| 大陆偷拍与自拍| 欧美一区二区亚洲| 十八禁网站网址无遮挡 | 欧美3d第一页| 哪个播放器可以免费观看大片| 秋霞在线观看毛片| 国产精品偷伦视频观看了| 97精品久久久久久久久久精品| 男的添女的下面高潮视频| av不卡在线播放| 肉色欧美久久久久久久蜜桃| 久久久色成人| av网站免费在线观看视频| 国产片特级美女逼逼视频| 国产精品一区二区性色av| 日本欧美国产在线视频| 日韩欧美 国产精品| 男人和女人高潮做爰伦理| 亚洲电影在线观看av| 夜夜骑夜夜射夜夜干| 一级毛片黄色毛片免费观看视频| 国内精品宾馆在线| 男女免费视频国产| 水蜜桃什么品种好| 亚洲av欧美aⅴ国产| 下体分泌物呈黄色| 亚洲精品日韩在线中文字幕| 色视频www国产| 街头女战士在线观看网站| 韩国高清视频一区二区三区| 一区在线观看完整版| 欧美 日韩 精品 国产| 精华霜和精华液先用哪个| 亚洲av综合色区一区| 亚洲美女黄色视频免费看| 啦啦啦啦在线视频资源| 最黄视频免费看| 中文字幕免费在线视频6| 国内揄拍国产精品人妻在线| 伊人久久国产一区二区| 亚洲自偷自拍三级| 一级毛片aaaaaa免费看小| 色哟哟·www| 午夜视频国产福利| 中国国产av一级| 国产伦理片在线播放av一区| av国产精品久久久久影院| 中文资源天堂在线| 国产av国产精品国产| 成人高潮视频无遮挡免费网站| 嫩草影院新地址| 欧美激情极品国产一区二区三区 | 成年美女黄网站色视频大全免费 | 国产一区亚洲一区在线观看| 久久婷婷青草| 中国国产av一级| 91狼人影院| 亚洲精品久久午夜乱码| 久久女婷五月综合色啪小说| 我要看黄色一级片免费的| 男女免费视频国产| 少妇 在线观看| 欧美亚洲 丝袜 人妻 在线| 国产亚洲一区二区精品| 精品国产乱码久久久久久小说| 熟妇人妻不卡中文字幕| 国产精品国产三级国产av玫瑰| 国产高清三级在线| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美精品人与动牲交sv欧美| 少妇人妻精品综合一区二区| 国产精品无大码| 中文字幕亚洲精品专区| 国产精品一及| 精品亚洲成a人片在线观看 | 久久久久久久国产电影| 日韩制服骚丝袜av| 国产一区有黄有色的免费视频| 国产一区二区在线观看日韩| 日本黄色片子视频| 又大又黄又爽视频免费| 又黄又爽又刺激的免费视频.| 国产亚洲欧美精品永久| 国产精品不卡视频一区二区| 亚洲精品久久午夜乱码| 成人亚洲欧美一区二区av| 免费观看的影片在线观看| 中文天堂在线官网| 热re99久久精品国产66热6| 国产一区二区在线观看日韩| 国产大屁股一区二区在线视频| 免费不卡的大黄色大毛片视频在线观看| 日韩视频在线欧美| 成人漫画全彩无遮挡| 欧美激情国产日韩精品一区| 嫩草影院入口| 插逼视频在线观看| 亚洲人与动物交配视频| 日本午夜av视频| 亚洲精品一区蜜桃| 午夜激情久久久久久久| 亚洲精品一二三| 91在线精品国自产拍蜜月| av福利片在线观看| 亚洲国产精品成人久久小说| av播播在线观看一区| 欧美日韩亚洲高清精品| 中国三级夫妇交换| 人人妻人人看人人澡| 日本av手机在线免费观看| 国产精品一区二区在线观看99| 久久精品人妻少妇| 国产精品一二三区在线看| 看非洲黑人一级黄片| 中文在线观看免费www的网站| 中文精品一卡2卡3卡4更新| 午夜免费男女啪啪视频观看| 欧美zozozo另类| 欧美日韩综合久久久久久| av女优亚洲男人天堂| 欧美亚洲 丝袜 人妻 在线| 午夜精品国产一区二区电影| 一级毛片 在线播放| av专区在线播放| 深爱激情五月婷婷| 成人亚洲欧美一区二区av| 精品一区二区三区视频在线| 欧美日韩在线观看h| 直男gayav资源| 777米奇影视久久| 97超碰精品成人国产| 久久青草综合色| 国产av精品麻豆| 丝袜脚勾引网站| 欧美精品一区二区免费开放| 亚洲不卡免费看| 久久人人爽人人爽人人片va| 欧美激情国产日韩精品一区| 亚洲va在线va天堂va国产| 国产爱豆传媒在线观看| 女的被弄到高潮叫床怎么办| 日韩中文字幕视频在线看片 | 午夜福利视频精品| 永久免费av网站大全| 国产69精品久久久久777片| 亚洲欧洲国产日韩| 男女无遮挡免费网站观看| 偷拍熟女少妇极品色| 国产精品一区二区在线观看99| 女性被躁到高潮视频| 自拍偷自拍亚洲精品老妇| 老熟女久久久| 黑人猛操日本美女一级片| 黄色配什么色好看| 女性被躁到高潮视频| 中文欧美无线码| 亚洲av在线观看美女高潮| 少妇的逼好多水| 一个人看视频在线观看www免费| 国产免费一区二区三区四区乱码| 中文字幕av成人在线电影| 极品教师在线视频| 亚洲自偷自拍三级| 国产精品99久久久久久久久| 夜夜骑夜夜射夜夜干| 汤姆久久久久久久影院中文字幕| 亚洲人成网站在线播| 啦啦啦在线观看免费高清www| 欧美老熟妇乱子伦牲交| 大片免费播放器 马上看| 欧美精品人与动牲交sv欧美| 韩国高清视频一区二区三区| 一级毛片 在线播放| 国产精品人妻久久久久久| 免费观看无遮挡的男女| 永久免费av网站大全| 欧美成人精品欧美一级黄| 亚州av有码| 国产伦理片在线播放av一区| 一区二区三区乱码不卡18| 午夜日本视频在线| 高清毛片免费看| 97在线视频观看| 看十八女毛片水多多多| 亚洲人与动物交配视频| 99热国产这里只有精品6| 免费人成在线观看视频色| 国产高清有码在线观看视频| 午夜福利网站1000一区二区三区| 国产精品蜜桃在线观看| 少妇的逼好多水| www.av在线官网国产| 麻豆成人午夜福利视频| 五月伊人婷婷丁香| 美女中出高潮动态图| 国产av国产精品国产| 两个人的视频大全免费| 99热这里只有是精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 另类亚洲欧美激情| av国产久精品久网站免费入址| 国产黄色免费在线视频| 中文字幕人妻熟人妻熟丝袜美| 舔av片在线| 人体艺术视频欧美日本| 性色av一级| 国产欧美另类精品又又久久亚洲欧美| 在线精品无人区一区二区三 | 国产亚洲5aaaaa淫片| 新久久久久国产一级毛片| 久久久久人妻精品一区果冻| 高清午夜精品一区二区三区| 国产精品国产三级国产专区5o| 亚洲伊人久久精品综合| 熟女av电影| 国产精品久久久久成人av| 国产午夜精品久久久久久一区二区三区| 国产色爽女视频免费观看| 视频中文字幕在线观看| 我的老师免费观看完整版| 午夜精品国产一区二区电影| 久久久精品94久久精品| 99视频精品全部免费 在线| 久久精品久久久久久久性| 大码成人一级视频| 亚洲欧美日韩东京热| 五月天丁香电影| 最新中文字幕久久久久| 小蜜桃在线观看免费完整版高清| 日韩欧美精品免费久久| 街头女战士在线观看网站| 精品久久久噜噜| 啦啦啦视频在线资源免费观看| 国内少妇人妻偷人精品xxx网站| av线在线观看网站| 国产精品一区二区在线不卡| av天堂中文字幕网| 成年免费大片在线观看| 亚洲av福利一区| 看十八女毛片水多多多| 日韩一区二区三区影片| 中国三级夫妇交换| 人妻制服诱惑在线中文字幕| 亚洲人与动物交配视频| 色吧在线观看| 精品少妇黑人巨大在线播放| 色婷婷久久久亚洲欧美| 日韩不卡一区二区三区视频在线| 久久久欧美国产精品| 99久久精品国产国产毛片| 少妇被粗大猛烈的视频| 欧美97在线视频| 六月丁香七月| 亚洲精品久久久久久婷婷小说| 国产精品麻豆人妻色哟哟久久| 国精品久久久久久国模美| 亚洲av中文av极速乱| 久久久久久久亚洲中文字幕| 亚洲,一卡二卡三卡| 国产黄色免费在线视频| 国产成人a区在线观看| 久久韩国三级中文字幕| 国产精品偷伦视频观看了| 男人添女人高潮全过程视频| 亚洲高清免费不卡视频| 91aial.com中文字幕在线观看| 精品视频人人做人人爽| 久久午夜福利片| 国产 一区精品| 少妇的逼水好多| 国内精品宾馆在线|