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

    基于精密單點(diǎn)定位技術(shù)的遠(yuǎn)程GPS驗(yàn)潮

    2016-11-07 05:43:17李圳章傳銀單瑞路媛琦
    全球定位系統(tǒng) 2016年4期
    關(guān)鍵詞:歷元潮位單點(diǎn)

    李圳,章傳銀,單瑞,路媛琦

    (1.山東科技大學(xué) 測(cè)繪科學(xué)與工程學(xué)院,青島 266510;2.中國(guó)測(cè)繪科學(xué)研究院,北京 100830;3.國(guó)土資源部海洋油氣資源與環(huán)境地質(zhì)重點(diǎn)實(shí)驗(yàn)室,青島 266071;4.青島海洋地質(zhì)研究所,青島 266071)

    ?

    基于精密單點(diǎn)定位技術(shù)的遠(yuǎn)程GPS驗(yàn)潮

    李圳1,2,章傳銀2,單瑞3,4,路媛琦1,2

    (1.山東科技大學(xué) 測(cè)繪科學(xué)與工程學(xué)院,青島 266510;2.中國(guó)測(cè)繪科學(xué)研究院,北京 100830;3.國(guó)土資源部海洋油氣資源與環(huán)境地質(zhì)重點(diǎn)實(shí)驗(yàn)室,青島 266071;4.青島海洋地質(zhì)研究所,青島 266071)

    精密單點(diǎn)定位(PPP)技術(shù)應(yīng)用于遠(yuǎn)程海水驗(yàn)潮,可有效彌補(bǔ)傳統(tǒng)驗(yàn)潮方法受距離限制等缺點(diǎn)。本文將單天的IGS觀測(cè)數(shù)據(jù)基于Bernese軟件進(jìn)行了單歷元PPP解算,探討了PPP應(yīng)用于海水驗(yàn)潮的有效性。并采用門限濾波方法提取潮汐信息,通過(guò)分析模擬潮位數(shù)據(jù)及西沙實(shí)測(cè)數(shù)據(jù)的內(nèi)符合精度證明,PPP驗(yàn)潮能夠達(dá)到厘米級(jí)的精度。

    精密單點(diǎn)定位;精度;潮汐觀測(cè);門限濾波

    0 引 言

    沿海地區(qū)是人類進(jìn)行生產(chǎn)生活活動(dòng)最頻繁的區(qū)域,然而在這一地帶潮汐現(xiàn)象顯著地影響著人們的生活,例如,沿海地區(qū)的農(nóng)田排灌,圍海造田,水產(chǎn)養(yǎng)殖,潮能發(fā)電,航行筑港等,與潮汐的漲落有著密切的聯(lián)系,因而,潮汐十分具有研究意義[1]。

    傳統(tǒng)潮位獲取方法有修建驗(yàn)潮站、拋投壓力式驗(yàn)潮儀或根據(jù)潮位模型推算得到,但它們都受成本過(guò)高或精度不夠等缺點(diǎn)的制約[2]。而PPP驗(yàn)潮是近幾年興起的潮位測(cè)量方法,已經(jīng)有相當(dāng)數(shù)量的學(xué)者做了相關(guān)研究。馮義楷等分析了影響GPS驗(yàn)潮精度的各種因素,提出了基于PPP技術(shù)的遠(yuǎn)程潮汐測(cè)量[3];范士杰等進(jìn)行了海上船載PPP驗(yàn)潮試驗(yàn),并借助自主研制的PPP軟件,取得了良好的驗(yàn)潮數(shù)據(jù)[4]。此外,國(guó)內(nèi)外學(xué)者對(duì)動(dòng)態(tài)PPP定位精度及其在海水驗(yàn)潮中的應(yīng)用均進(jìn)行了相關(guān)研究[5-7]。

    本文研究了基于Bernese軟件的單歷元PPP解算精度,并將門限濾波應(yīng)用于PPP驗(yàn)潮,從內(nèi)符合精度上驗(yàn)證了其可行性。

    1 PPP驗(yàn)潮

    1.1單歷元精密單點(diǎn)定位觀測(cè)方程

    由于潮位是實(shí)時(shí)變化著的,要想獲得某一瞬間的潮位數(shù)據(jù),就需要逐歷元解算GPS觀測(cè)數(shù)據(jù)。在單歷元情況下,由于受到更多的誤差影響,使得定位精度有所下降[8],精密單點(diǎn)定位通常使用載波相位觀測(cè)值,其線性化后的觀測(cè)方程為

    λφi=-liVx-miVy-niVz-cVtR-λNi+

    (1)

    1.2基于PPP的驗(yàn)潮模型

    GPS定位技術(shù)的發(fā)展,給PPP應(yīng)用于海水驗(yàn)潮提供了可能性,PPP驗(yàn)潮可以有效彌補(bǔ)傳統(tǒng)驗(yàn)潮方法受距離限制等缺陷。PPP驗(yàn)潮的基本流程及原理如圖1、圖2所示。

    圖1 PPP驗(yàn)潮的流程圖

    圖2 PPP驗(yàn)潮方法原理圖

    載體上安放接收機(jī),其正常高和大地高分別為hkm、hkg,h為接收機(jī)天線高。潮位通常以海圖高來(lái)描述,因而需要對(duì)高程基準(zhǔn)進(jìn)行轉(zhuǎn)換,則流動(dòng)站的瞬時(shí)潮位為

    T=hkm-h-h0=hkg-ξk-h-h0,

    (2)

    式中:h0為水尺零點(diǎn)高程,通常借助潮汐表獲取;ξk為高程異常。

    由于受波浪影響,載體姿態(tài)也會(huì)隨之發(fā)生相應(yīng)改變,從而對(duì)PPP驗(yàn)潮獲取的數(shù)據(jù)產(chǎn)生一定影響[11]。GPS測(cè)量大地高,它所測(cè)得的瞬時(shí)海面并不代表實(shí)際潮位面T(t),為此需要進(jìn)行門限濾波,用以消除波浪ω的影響如圖3所示。

    不考慮其他因素影響,瞬時(shí)海面Ts可看作由潮汐T和波浪ω組成,表達(dá)式為

    Ts=T(t)+ω(t) ,

    (3)

    (4)

    G(V0i+ui)-gi) ,

    (5)

    式中:ai,bi為波浪i在正交方向上的分振幅; ?i為波浪i的周期;MSL0為長(zhǎng)期平均水位高度;fi為分潮i的交點(diǎn)因子;Hi為分潮i的平均振幅;qi為分潮i的角速率;V0i為分潮i的格林尼治零時(shí)天文初相角;ui為分潮i點(diǎn)訂正角;gi為分潮的區(qū)時(shí)專用遲角。

    由于波浪的涌動(dòng),載體發(fā)生變化,其影響較大,根據(jù)式(5),潮位周期一般大于1 h;涌浪周期一般為1~30 s.兩項(xiàng)的周期變化如圖4所示。對(duì)式(3)兩邊在[0,Δt]區(qū)間取積分平均:

    (6)

    圖4說(shuō)明,只要選擇適當(dāng)?shù)摩,便可以實(shí)現(xiàn)涌浪的濾除。

    圖3 涌浪影響示意圖

    圖4 水位及波浪周期示意圖

    2 實(shí)驗(yàn)與精度分析

    2.1單歷元精密單點(diǎn)定位精度分析

    由于本文的研究重點(diǎn)是PPP在海水驗(yàn)潮中的應(yīng)用,此時(shí)載體是不斷運(yùn)動(dòng)變化著的,探討的是接收機(jī)在某一時(shí)間的瞬時(shí)位置,故采用了單歷元解算模式,分析了動(dòng)態(tài)PPP的定位精度。

    為了驗(yàn)證單歷元精密單點(diǎn)定位精度,本文選取了位于中國(guó)境內(nèi)的TNML(新竹)IGS跟蹤站2015年3月5日的觀測(cè)數(shù)據(jù),并結(jié)合IGS提供的精密星歷、鐘差等,基于Bernese5.0軟件對(duì)數(shù)據(jù)進(jìn)行了單歷元解算,其中測(cè)站單歷元解與全天解坐標(biāo)之差的統(tǒng)計(jì)結(jié)果如表1所示,N、E、U三個(gè)坐標(biāo)分量相對(duì)于全天所有歷元解算結(jié)果平均值的偏差如圖5所示。

    表1 測(cè)站單歷元解與全天解坐標(biāo)之差的統(tǒng)計(jì)結(jié)果

    圖5 N、E、U方向上測(cè)站單歷元解結(jié)果偏差圖

    根據(jù)表1和圖5中可以看出,動(dòng)態(tài)PPP在N、E、U三個(gè)坐標(biāo)分量上的精度基本均勻,TNML測(cè)站的單歷元PPP定位精度為5.7 cm,其中,高程方向上精度優(yōu)于5 cm,據(jù)此可知單歷元PPP的解算結(jié)果略低于GPS RTK技術(shù)的精度。但單歷元PPP只需要一臺(tái)接收機(jī),無(wú)需進(jìn)行數(shù)據(jù)差分,且其精度均勻,所以將單歷元PPP應(yīng)用于海水驗(yàn)潮,特別是遠(yuǎn)離大陸的遠(yuǎn)程海水驗(yàn)潮是切實(shí)可行的。

    2.2模擬潮位數(shù)據(jù)的構(gòu)造設(shè)計(jì)與門限濾波處理

    由于潮位數(shù)據(jù)主要由潮汐和波浪及其他一些噪聲組成,所以本文設(shè)潮位數(shù)據(jù)由兩個(gè)周期分別為12 h和6 h的潮位,周期為15 s的波浪和其他各種環(huán)境影響造成的服從(0,0.1)分布的隨機(jī)噪聲組成,則潮位可用下式表示:

    f1(t)=10×sin[2π/(3600×12)×t+π/3],

    (7)

    f2(t)=5×sin[2π/(3600×6)×t+π/6],

    (8)

    fw=0.5×sin(2π/15×t),

    (9)

    fh(t)=0.1×N(0,1),

    (10)

    f(t)=f1(t)+f2(t)+fw(t)+fh(t),

    (11)

    其中: f1、f2為潮汐函數(shù); fw為波浪函數(shù); fh為噪聲。

    本文基于MATLAB,根據(jù)式(7)~式(11)構(gòu)建了模擬潮位數(shù)據(jù)f(t),對(duì)f(t)共采樣了1 024個(gè)數(shù)據(jù),采樣頻率為1 Hz,結(jié)果如圖6所示。為驗(yàn)證門限濾波效果,通過(guò)編程實(shí)現(xiàn)了模擬數(shù)據(jù)的門限濾波,繪制了濾波后的潮位變化曲線和理論潮位變化曲線,其結(jié)果分別如圖7、圖8所示,表2示出了兩者之間差值的統(tǒng)計(jì)結(jié)果。

    圖6 模擬潮位數(shù)據(jù)f(t)

    圖7 門限濾波處理后的模擬潮位數(shù)據(jù)

    圖8 理論潮位數(shù)據(jù)

    參數(shù)Max/mMin/mMean/mRms/m差值大小0.1536-0.17850.01730.0510

    通過(guò)對(duì)圖7和圖8的對(duì)比可以發(fā)現(xiàn),門限濾波很好的過(guò)濾了干擾。由表2可知,濾波結(jié)果與理論潮位值的標(biāo)準(zhǔn)差為5.1 cm,平均差異為1.7 cm,這反映了在濾波過(guò)程中對(duì)波浪的干擾濾除作用較好。

    2.3PPP驗(yàn)潮實(shí)驗(yàn)

    本文利用3臺(tái)GPS浮標(biāo)觀測(cè)裝置,于2014年5月,在我國(guó)西沙群島某島嶼附近同一海域使用三臺(tái)接收機(jī)進(jìn)行了同步PPP驗(yàn)潮實(shí)驗(yàn),采樣間隔為60 s,浮標(biāo)相對(duì)位置關(guān)系如圖9所示。

    圖9 浮標(biāo)相對(duì)位置關(guān)系示意圖

    考慮浮標(biāo)較容易受海面異常狀態(tài)的影響,故先對(duì)原始驗(yàn)潮數(shù)據(jù)進(jìn)行了粗差剔除,粗差的剔除主要依靠設(shè)置閾值實(shí)現(xiàn),閾值如果過(guò)大則粗差探測(cè)不靈敏,剔除不干凈,如果過(guò)小則可能會(huì)剔除一些比較干凈的數(shù)據(jù),因此應(yīng)當(dāng)反復(fù)實(shí)驗(yàn)獲得一個(gè)恰當(dāng)?shù)拈撝?。在剔除粗差?為了減小波浪影響,對(duì)潮位數(shù)據(jù)進(jìn)行了門限濾波。在數(shù)據(jù)預(yù)處理工作完成之后,由于點(diǎn)云數(shù)據(jù)不易觀察分析,故采用了八階高斯擬合處理。

    此外,三號(hào)接收機(jī)在150~250歷元的高程有異?,F(xiàn)象,出現(xiàn)了15 cm左右的系統(tǒng)誤差,此處應(yīng)當(dāng)為浮標(biāo)錨定不穩(wěn)造成的。為了分析該誤差,對(duì)100~300歷元間的驗(yàn)潮數(shù)據(jù)進(jìn)行了相鄰歷元作差,其結(jié)果如圖10所示。

    圖10 相鄰歷元作差結(jié)果

    從圖中可以看出,該浮標(biāo)在155歷元和250歷元附近發(fā)生了15 cm左右的較大變化,但在155至250歷元之間浮標(biāo)的變化卻較為平穩(wěn)。這說(shuō)明了浮標(biāo)在155歷元時(shí)因外力影響使錨定異常,在250歷元時(shí)及時(shí)發(fā)現(xiàn)并修復(fù),從而使得這一期間的測(cè)量數(shù)據(jù)出現(xiàn)了誤差,故在分析時(shí)將這一時(shí)段的潮位數(shù)據(jù)減小了15 cm.原始潮位數(shù)據(jù)及預(yù)處理后驗(yàn)潮數(shù)據(jù)的擬合曲線如圖11所示。

    圖11 三臺(tái)儀器所觀測(cè)的原始潮位數(shù)據(jù)及預(yù)處理后擬合曲線(a)一號(hào)接收機(jī); (b)二號(hào)接收機(jī); (c)三號(hào)接收機(jī)

    驗(yàn)潮數(shù)據(jù)經(jīng)過(guò)預(yù)處理后,得到比較干凈的觀測(cè)值,之后將觀測(cè)時(shí)間同步,三臺(tái)接收機(jī)所得潮位變化值如圖12所示。可以從圖中看出,三臺(tái)儀器所測(cè)量的漲潮時(shí)間、落潮時(shí)間與潮差基本相同,潮位變化相當(dāng)一致,將三臺(tái)儀器測(cè)量的潮高均值作為真值,研究每一歷元的潮高,分析其內(nèi)符合精度,則各儀器與真值的歷元之差間統(tǒng)計(jì)結(jié)果如表3所示。

    圖12 三臺(tái)接收機(jī)所測(cè)潮位變化值

    表3 各接收機(jī)與真值的歷元之差間統(tǒng)計(jì)結(jié)果表

    可以看出,三臺(tái)儀器的潮高數(shù)據(jù)內(nèi)符合精度很好,其平均偏差僅有1 cm左右,均方根在6 cm之內(nèi),但第一臺(tái)接收機(jī)測(cè)得數(shù)據(jù)的RMS值較大,這在圖13中也可以比較清晰的看出其潮高與其他兩組數(shù)據(jù)相比有明顯的差異,在高潮和低潮階段更加明顯,這應(yīng)該與浮標(biāo)分布的不同海域潮汐存在微小差異有關(guān)。

    3 結(jié)束語(yǔ)

    本文研究了單歷元PPP應(yīng)用于海水驗(yàn)潮的可能性,并將門限濾波應(yīng)用于驗(yàn)潮數(shù)據(jù)分析,于西沙群島某海域進(jìn)行了相關(guān)實(shí)驗(yàn),得出了如下結(jié)論:

    1) Bernese軟件解算測(cè)站單歷元解的精度在平面上優(yōu)于6 cm,高程方向上優(yōu)于5 cm,略低于GPS RTK的定位精度,但是其點(diǎn)位誤差不受距離變化的影響,也無(wú)需陸地基站提供差分信號(hào)。

    2) 門限濾波對(duì)數(shù)學(xué)函數(shù)構(gòu)造的模擬潮位數(shù)據(jù)起到了較好去噪作用,研究西沙實(shí)驗(yàn)數(shù)據(jù)進(jìn)行粗差剔除與門限濾波之后的結(jié)果發(fā)現(xiàn),其內(nèi)符合精度很高,平均偏差僅有1 cm左右,均方根在6 cm之內(nèi)。

    綜上所述,利用單歷元PPP定位技術(shù),并通過(guò)相關(guān)的數(shù)據(jù)處理之后進(jìn)行海水驗(yàn)潮可以達(dá)到厘米級(jí)的精度,與現(xiàn)在廣泛采用的GPS RTK驗(yàn)潮技術(shù)所達(dá)到的精度基本相同。但PPP的優(yōu)勢(shì)在于其不依賴陸地基準(zhǔn)站,精度均勻,因而PPP應(yīng)用于遠(yuǎn)程海水驗(yàn)潮具有一定的實(shí)用價(jià)值。

    [1] 汪連賀.基于GPS PPK技術(shù)的遠(yuǎn)距離高精度驗(yàn)潮方法研究[J].海洋測(cè)繪,2014,34(4):24-27.

    [2] 趙建虎.現(xiàn)代海洋測(cè)繪[M]. 武漢:武漢大學(xué)出版社, 2007.

    [3] 馮義楷,李杰,楊龍,等.遠(yuǎn)程GPS驗(yàn)潮方法研究[J].海洋測(cè)繪,2010,30(1):4-6.

    [4] 范士杰,秦學(xué)彬,吳紹玉,等.不同星歷和鐘差產(chǎn)品的PPP驗(yàn)潮試驗(yàn)及結(jié)果分析[J].海洋測(cè)繪,2014,34(4):43-46.

    [5] MUELLERSCHOEN R J,BERTIGER W I,LOUGH M,etal. An internet-based global differential GPS system,initial result [C]//Proceedings of ION GPS 2000,Salt Lake City: UT, September,2000.

    [6] 張小紅,李星星,郭斐,等.GPS單頻精密單點(diǎn)定位軟件實(shí)現(xiàn)與精度分析[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2008,33(8):783-787.

    [7] 程世來(lái),張小紅.基于PPP技術(shù)的GPS浮標(biāo)海嘯預(yù)警模擬研究[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2007,32(9):764-766.

    [8] 胡艷,胡劍.GPS單歷元精密單點(diǎn)定位精度分析[J].山東科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,32(1):74-78.

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

    [10]李施佳,蔡昌盛,李佳霖,等.對(duì)流層投影函數(shù)對(duì)精密單點(diǎn)定位精度的影響分析[J].全球定位系統(tǒng),2013,38(5):22-26.

    [11]陽(yáng)凡林,趙建虎.GPS驗(yàn)潮中波浪的誤差分析和消除[J].海洋測(cè)繪,2003,23(3):1-4.

    Remote GPS Tide Observing Method Based on the Precise Point Positioning

    LI Zhen1,2,ZHANG Chuanyin2,SHAN Rui3,4,LU Yuanqi1,2

    (1.CollegeofGeomatics,ShandongUniversityofScienceandTechnology,Qingdao266510,China; 2.ChineseAcademyofSurveyingandMapping,Beijing100830,China; 3.KeyLaboratoryofMarineHydrocarbonResourcesandEnvironmentalGeology,Qingdao266071,China; 4.MinistryofLandandResources,QingdaoInstituteofMarineGeology,Qingdao266071,China)

    The precise point position (PPP) technology is applied to the remote tide observing, which can effectively make up the disadvantage of traditional tidal methods limited by the distance. The paper makes the single epoch PPP calculation for the observation data of single IGS based on Bernese software to discuss the feasibility of PPP applied on the tide observing. PPP tide observing can achieve centimeter level accuracy by the tide information extracted from threshold filtering method and analyzing the analysis of simulated tidal data and the inner precision demonstration of Xisha measured data.

    Precise Point Position; accuracy; tide observing; threshold filtering

    10.13442/j.gnss.1008-9268.2016.04.012

    2016-03-07

    國(guó)家自然科學(xué)基金(批準(zhǔn)號(hào):41406115)

    P228.4

    A

    1008-9268(2016)04-0053-06

    李圳(1993-),男,碩士生,研究方向?yàn)榭臻g大地測(cè)量學(xué)。

    章傳銀(1968-),男,研究員,主要從事地球重力場(chǎng)與大地測(cè)量學(xué)研究。

    單瑞(1985-),男,助理研究員,主要從事海洋測(cè)繪研究。

    路媛琦(1992-),女,山東青島人,碩士生,研究方向?yàn)樾l(wèi)星定位與導(dǎo)航。

    聯(lián)系人: 李圳 E-mail: sdkjlizhen@163.com

    猜你喜歡
    歷元潮位單點(diǎn)
    基于距離倒數(shù)加權(quán)的多站潮位改正方法可行性分析
    唐山市警戒潮位標(biāo)志物維護(hù)研究
    歷元間載波相位差分的GPS/BDS精密單點(diǎn)測(cè)速算法
    超薄異型坯連鑄機(jī)非平衡單點(diǎn)澆鑄實(shí)踐與分析
    山東冶金(2019年5期)2019-11-16 09:09:10
    多潮位站海道地形測(cè)量潮位控制方法研究
    基于改進(jìn)的OLS-RBF模型的感潮河段潮位預(yù)測(cè)研究
    Recent advances of TCM treatment of childhood atopic dermatitis
    數(shù)字電視地面?zhèn)鬏斢脝晤l網(wǎng)與單點(diǎn)發(fā)射的效果比較
    Clinical observation of Huatan Huoxue Formula in treating coronary heart disease with hyperlipidemia
    Mechanism of sex hormone level in biological clock disorder induced acne and analysis of TCM Pathogenesis
    亚洲在线观看片| 亚洲自偷自拍三级| 97超视频在线观看视频| 精品久久国产蜜桃| 亚洲最大成人手机在线| 国产亚洲精品久久久com| 亚洲精品乱码久久久v下载方式| 乱码一卡2卡4卡精品| 亚洲丝袜综合中文字幕| 精品人妻偷拍中文字幕| 久久这里只有精品中国| 日韩高清综合在线| 国产日韩欧美在线精品| 国产免费一级a男人的天堂| av免费在线看不卡| 久久国内精品自在自线图片| 国产高清三级在线| 欧美最新免费一区二区三区| 午夜a级毛片| 免费电影在线观看免费观看| 亚洲av免费高清在线观看| 99热全是精品| 搞女人的毛片| 国产一区二区在线观看日韩| 五月玫瑰六月丁香| 亚洲国产精品成人久久小说| 天堂av国产一区二区熟女人妻| 精品一区二区三区人妻视频| 亚洲av成人av| 亚洲va在线va天堂va国产| 国产精品伦人一区二区| 日韩国内少妇激情av| 黄片无遮挡物在线观看| 欧美zozozo另类| 嘟嘟电影网在线观看| 亚洲无线观看免费| 婷婷色综合大香蕉| 国模一区二区三区四区视频| 国产av一区在线观看免费| 色综合色国产| 久久久久久久久大av| 乱人视频在线观看| 如何舔出高潮| 白带黄色成豆腐渣| 美女黄网站色视频| 能在线免费观看的黄片| 久久久久免费精品人妻一区二区| av免费观看日本| 少妇裸体淫交视频免费看高清| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 久久国产乱子免费精品| 六月丁香七月| 久久久久久久久久久丰满| 网址你懂的国产日韩在线| 少妇丰满av| 日韩av在线大香蕉| 久久精品国产亚洲av天美| 白带黄色成豆腐渣| 国产精品久久久久久精品电影| 18禁在线无遮挡免费观看视频| 亚洲久久久久久中文字幕| 黄色一级大片看看| 久久久久久久久久黄片| 麻豆成人午夜福利视频| 少妇的逼水好多| 亚洲精品亚洲一区二区| 天美传媒精品一区二区| 国产精品永久免费网站| 国产av码专区亚洲av| 亚洲精品乱久久久久久| 51国产日韩欧美| 波多野结衣巨乳人妻| 国产三级中文精品| 免费观看性生交大片5| 成人鲁丝片一二三区免费| 特级一级黄色大片| 九九久久精品国产亚洲av麻豆| 18禁在线无遮挡免费观看视频| 在线播放国产精品三级| 欧美一区二区精品小视频在线| 97在线视频观看| 热99re8久久精品国产| 国产午夜精品一二区理论片| 国产真实伦视频高清在线观看| 一夜夜www| 成年av动漫网址| 国产精品一区二区在线观看99 | 91久久精品电影网| 深爱激情五月婷婷| 久热久热在线精品观看| 99视频精品全部免费 在线| 丰满人妻一区二区三区视频av| 国产亚洲av嫩草精品影院| 国产成人精品一,二区| 大话2 男鬼变身卡| 精品久久久久久久久av| 91精品国产九色| 国语自产精品视频在线第100页| 久久婷婷人人爽人人干人人爱| 午夜老司机福利剧场| 国产精品麻豆人妻色哟哟久久 | 变态另类丝袜制服| 成人二区视频| 99久久人妻综合| 九草在线视频观看| 国产高清视频在线观看网站| 欧美丝袜亚洲另类| 少妇熟女aⅴ在线视频| 久久久久精品久久久久真实原创| 最近手机中文字幕大全| 免费av观看视频| 亚洲自拍偷在线| 日日摸夜夜添夜夜爱| 如何舔出高潮| 人妻夜夜爽99麻豆av| 综合色av麻豆| 91精品国产九色| 黄色日韩在线| 国产一区二区亚洲精品在线观看| 人妻系列 视频| 国产成人freesex在线| av免费观看日本| 久久午夜福利片| 国产高清三级在线| 日韩精品青青久久久久久| 18禁在线播放成人免费| 99久久成人亚洲精品观看| 成年女人看的毛片在线观看| 99九九线精品视频在线观看视频| 直男gayav资源| 亚洲第一区二区三区不卡| 国产av不卡久久| 亚洲四区av| 精品久久久久久久人妻蜜臀av| 国产一区二区三区av在线| 69av精品久久久久久| 天堂影院成人在线观看| 日产精品乱码卡一卡2卡三| 两个人视频免费观看高清| 天堂影院成人在线观看| 亚洲精品久久久久久婷婷小说 | 日韩欧美在线乱码| 国产伦精品一区二区三区视频9| 春色校园在线视频观看| 欧美+日韩+精品| 久久久久精品久久久久真实原创| 久99久视频精品免费| 高清视频免费观看一区二区 | 亚洲高清免费不卡视频| 成人欧美大片| 最新中文字幕久久久久| 麻豆成人午夜福利视频| 国产免费男女视频| 国产一区亚洲一区在线观看| 99久国产av精品| 国产白丝娇喘喷水9色精品| 亚洲av免费在线观看| 久久精品久久精品一区二区三区| 亚洲精品乱码久久久v下载方式| 国产精品一区二区在线观看99 | 精品一区二区免费观看| 免费黄色在线免费观看| 欧美区成人在线视频| 一级毛片电影观看 | av在线观看视频网站免费| 少妇猛男粗大的猛烈进出视频 | 亚洲精品乱码久久久v下载方式| 国产成人精品久久久久久| or卡值多少钱| 久久久久久久久久久丰满| 国产精品综合久久久久久久免费| 观看美女的网站| 又爽又黄a免费视频| 欧美一区二区精品小视频在线| 精品免费久久久久久久清纯| 亚洲三级黄色毛片| 免费观看a级毛片全部| 一区二区三区乱码不卡18| 赤兔流量卡办理| 3wmmmm亚洲av在线观看| 午夜激情欧美在线| 精品国产一区二区三区久久久樱花 | 五月玫瑰六月丁香| 中文字幕制服av| 三级国产精品欧美在线观看| 欧美zozozo另类| 欧美97在线视频| 亚洲欧美日韩无卡精品| 麻豆成人午夜福利视频| 欧美最新免费一区二区三区| 欧美性猛交╳xxx乱大交人| 国产久久久一区二区三区| 永久免费av网站大全| 久久99热这里只频精品6学生 | 草草在线视频免费看| 国产淫片久久久久久久久| 高清视频免费观看一区二区 | 国产精品熟女久久久久浪| 久久久久免费精品人妻一区二区| 精品午夜福利在线看| 久久久国产成人精品二区| 你懂的网址亚洲精品在线观看 | 女人久久www免费人成看片 | 亚洲成人av在线免费| 久久人人爽人人爽人人片va| 高清在线视频一区二区三区 | 99久久精品国产国产毛片| 亚洲欧美成人精品一区二区| av又黄又爽大尺度在线免费看 | 91午夜精品亚洲一区二区三区| 免费观看性生交大片5| 免费观看精品视频网站| 精品国产三级普通话版| 三级毛片av免费| www.av在线官网国产| 亚洲国产欧美在线一区| 中文字幕免费在线视频6| 日本免费一区二区三区高清不卡| 一本久久精品| 国产午夜精品论理片| 日日撸夜夜添| 如何舔出高潮| 免费看a级黄色片| 禁无遮挡网站| av天堂中文字幕网| 天美传媒精品一区二区| 中文字幕久久专区| 国产精品,欧美在线| 国产乱来视频区| 九九在线视频观看精品| 色播亚洲综合网| 欧美97在线视频| 婷婷色综合大香蕉| 色视频www国产| 亚洲电影在线观看av| 18禁在线播放成人免费| 午夜视频国产福利| 五月玫瑰六月丁香| 老师上课跳d突然被开到最大视频| 国产精品国产三级专区第一集| 能在线免费看毛片的网站| 亚洲无线观看免费| 成人特级av手机在线观看| 国产av一区在线观看免费| 永久免费av网站大全| 久久综合国产亚洲精品| 床上黄色一级片| 极品教师在线视频| 久久99蜜桃精品久久| 老司机影院成人| 色综合色国产| 91久久精品国产一区二区三区| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩无卡精品| 免费黄色在线免费观看| 一区二区三区高清视频在线| 亚洲精品影视一区二区三区av| 在线播放无遮挡| 97超视频在线观看视频| 级片在线观看| 日韩精品有码人妻一区| 男人狂女人下面高潮的视频| 久久久午夜欧美精品| 99久久精品国产国产毛片| 国产探花极品一区二区| 国产精品日韩av在线免费观看| 午夜久久久久精精品| 2022亚洲国产成人精品| 国产精品美女特级片免费视频播放器| 国产午夜福利久久久久久| 一个人观看的视频www高清免费观看| 91精品一卡2卡3卡4卡| 亚洲精品乱久久久久久| 高清在线视频一区二区三区 | 最近2019中文字幕mv第一页| 国产高潮美女av| 日韩欧美精品免费久久| 亚洲成色77777| 一级毛片久久久久久久久女| 亚洲欧美成人精品一区二区| 久久久久久大精品| 可以在线观看毛片的网站| 色综合亚洲欧美另类图片| 国产一级毛片在线| 自拍偷自拍亚洲精品老妇| 午夜福利在线观看吧| 欧美人与善性xxx| 色综合站精品国产| 国产探花在线观看一区二区| av免费在线看不卡| 中文字幕人妻熟人妻熟丝袜美| 日本色播在线视频| 免费看av在线观看网站| av在线观看视频网站免费| 国产精品久久久久久精品电影小说 | 五月玫瑰六月丁香| 日本一本二区三区精品| 热99在线观看视频| 亚洲精品乱码久久久v下载方式| 亚洲自拍偷在线| 熟女人妻精品中文字幕| 成人毛片a级毛片在线播放| 国模一区二区三区四区视频| 一级黄色大片毛片| 卡戴珊不雅视频在线播放| 国产高清有码在线观看视频| 人人妻人人看人人澡| 久久久久久久久大av| 久久久久久久国产电影| 久久精品久久久久久噜噜老黄 | 亚洲综合色惰| 欧美高清性xxxxhd video| 国产极品天堂在线| 午夜激情福利司机影院| 麻豆久久精品国产亚洲av| 97热精品久久久久久| 熟妇人妻久久中文字幕3abv| 日本av手机在线免费观看| 国产成人午夜福利电影在线观看| 色哟哟·www| 久久精品国产亚洲av天美| 婷婷色综合大香蕉| 特大巨黑吊av在线直播| 黄片wwwwww| 麻豆av噜噜一区二区三区| ponron亚洲| 国产精品嫩草影院av在线观看| 最近中文字幕2019免费版| 色播亚洲综合网| 亚洲欧美日韩东京热| 久久久精品欧美日韩精品| 亚洲三级黄色毛片| 亚洲欧洲日产国产| 国产伦精品一区二区三区四那| 国产成年人精品一区二区| 国产又黄又爽又无遮挡在线| 97超视频在线观看视频| 亚洲在线观看片| 97超视频在线观看视频| 男女那种视频在线观看| av又黄又爽大尺度在线免费看 | 国产一区二区在线观看日韩| 国产91av在线免费观看| 日本wwww免费看| 欧美极品一区二区三区四区| 亚洲美女视频黄频| 99久久人妻综合| 精品国产三级普通话版| 亚洲精品一区蜜桃| 美女黄网站色视频| 黄色欧美视频在线观看| 亚洲国产欧美人成| 美女被艹到高潮喷水动态| 高清av免费在线| 久久久精品大字幕| 久久精品影院6| 免费观看精品视频网站| 97在线视频观看| 午夜视频国产福利| a级毛片免费高清观看在线播放| 国产精品野战在线观看| 亚洲国产欧美人成| 亚洲成av人片在线播放无| 亚洲最大成人av| 日韩在线高清观看一区二区三区| 青春草亚洲视频在线观看| 欧美高清成人免费视频www| 淫秽高清视频在线观看| 日韩成人伦理影院| 日日干狠狠操夜夜爽| 搞女人的毛片| 一二三四中文在线观看免费高清| 国产高清不卡午夜福利| 最近手机中文字幕大全| 亚洲精品aⅴ在线观看| 国产精品熟女久久久久浪| 大话2 男鬼变身卡| 狂野欧美白嫩少妇大欣赏| 午夜日本视频在线| 超碰97精品在线观看| 国产真实伦视频高清在线观看| 国产单亲对白刺激| 国内揄拍国产精品人妻在线| 国产在线一区二区三区精 | 如何舔出高潮| 久久婷婷人人爽人人干人人爱| 久久精品国产鲁丝片午夜精品| 色综合亚洲欧美另类图片| 天堂影院成人在线观看| 欧美精品国产亚洲| 欧美变态另类bdsm刘玥| 偷拍熟女少妇极品色| 国产伦在线观看视频一区| 一级二级三级毛片免费看| 欧美日韩国产亚洲二区| 最近中文字幕高清免费大全6| 成人高潮视频无遮挡免费网站| av视频在线观看入口| 看十八女毛片水多多多| 男插女下体视频免费在线播放| 国产免费福利视频在线观看| 热99re8久久精品国产| 蜜臀久久99精品久久宅男| 国产69精品久久久久777片| 久久国产乱子免费精品| 国产精品永久免费网站| 国产高清视频在线观看网站| 小蜜桃在线观看免费完整版高清| 18禁在线播放成人免费| 女人久久www免费人成看片 | 日本与韩国留学比较| 色视频www国产| 欧美日韩精品成人综合77777| 免费黄网站久久成人精品| 亚洲av男天堂| 成人综合一区亚洲| 最近视频中文字幕2019在线8| 国产伦一二天堂av在线观看| av视频在线观看入口| 欧美最新免费一区二区三区| 国产精品嫩草影院av在线观看| 国产爱豆传媒在线观看| 亚洲一区高清亚洲精品| 精品国产三级普通话版| 亚洲在线自拍视频| 国产精品日韩av在线免费观看| 中文资源天堂在线| 国产成人a∨麻豆精品| 精品国产露脸久久av麻豆 | 久热久热在线精品观看| 美女国产视频在线观看| 联通29元200g的流量卡| 成人漫画全彩无遮挡| 国产伦理片在线播放av一区| 国产午夜福利久久久久久| 日日摸夜夜添夜夜爱| 欧美日本视频| 国产精品日韩av在线免费观看| 秋霞在线观看毛片| 免费大片18禁| 99久国产av精品| 日本免费在线观看一区| or卡值多少钱| 久久精品熟女亚洲av麻豆精品 | 免费av不卡在线播放| 国产精品女同一区二区软件| 亚洲丝袜综合中文字幕| 久久精品影院6| 国产淫语在线视频| av线在线观看网站| 小蜜桃在线观看免费完整版高清| 免费看a级黄色片| 水蜜桃什么品种好| 国产精品一区二区在线观看99 | 中文天堂在线官网| 久久午夜福利片| 黄色一级大片看看| 亚洲av电影在线观看一区二区三区 | 亚洲成人精品中文字幕电影| 在线免费观看不下载黄p国产| 亚洲第一区二区三区不卡| 亚洲婷婷狠狠爱综合网| 两个人的视频大全免费| 国产亚洲一区二区精品| 日本wwww免费看| 午夜精品一区二区三区免费看| 欧美一区二区国产精品久久精品| 熟妇人妻久久中文字幕3abv| 国语对白做爰xxxⅹ性视频网站| 国产高潮美女av| 狂野欧美白嫩少妇大欣赏| 亚洲欧美日韩无卡精品| 日韩欧美国产在线观看| 床上黄色一级片| 夜夜爽夜夜爽视频| 久久久久久久久久黄片| 麻豆一二三区av精品| 超碰av人人做人人爽久久| 赤兔流量卡办理| 五月玫瑰六月丁香| 国产亚洲精品av在线| av卡一久久| 三级国产精品欧美在线观看| 超碰97精品在线观看| 男女下面进入的视频免费午夜| 久久久a久久爽久久v久久| 男女国产视频网站| 国语自产精品视频在线第100页| 男人和女人高潮做爰伦理| 日日啪夜夜撸| 精品99又大又爽又粗少妇毛片| 日韩三级伦理在线观看| 久久久欧美国产精品| 最近中文字幕高清免费大全6| 插逼视频在线观看| 一级黄片播放器| 久久这里只有精品中国| 可以在线观看毛片的网站| 一夜夜www| 国产熟女欧美一区二区| 桃色一区二区三区在线观看| 亚洲伊人久久精品综合 | 搡女人真爽免费视频火全软件| 日本一本二区三区精品| 一级毛片久久久久久久久女| 一个人看的www免费观看视频| 亚洲国产精品成人久久小说| av国产免费在线观看| 国产高清三级在线| 美女xxoo啪啪120秒动态图| 免费av不卡在线播放| 国产成人aa在线观看| 在线观看一区二区三区| 黄色欧美视频在线观看| 嫩草影院新地址| 亚洲电影在线观看av| 日韩欧美精品免费久久| 中文亚洲av片在线观看爽| 国产在视频线精品| 十八禁国产超污无遮挡网站| 在线天堂最新版资源| 欧美xxxx黑人xx丫x性爽| 看十八女毛片水多多多| av在线播放精品| 色综合色国产| 18禁在线无遮挡免费观看视频| 成年女人看的毛片在线观看| 成人av在线播放网站| 婷婷色麻豆天堂久久 | 久久鲁丝午夜福利片| 毛片一级片免费看久久久久| 亚洲色图av天堂| 色播亚洲综合网| 色网站视频免费| 春色校园在线视频观看| 亚洲中文字幕一区二区三区有码在线看| 1024手机看黄色片| 日本免费一区二区三区高清不卡| 国产黄色视频一区二区在线观看 | 国产片特级美女逼逼视频| 男的添女的下面高潮视频| 乱系列少妇在线播放| 男人舔奶头视频| 久久精品久久久久久久性| 亚洲欧美精品专区久久| 中文字幕免费在线视频6| 久久婷婷人人爽人人干人人爱| 99久久成人亚洲精品观看| 日韩欧美三级三区| 亚洲三级黄色毛片| 淫秽高清视频在线观看| 免费一级毛片在线播放高清视频| 国产精品精品国产色婷婷| 内射极品少妇av片p| 国产精品久久久久久av不卡| 免费av不卡在线播放| 春色校园在线视频观看| 午夜福利网站1000一区二区三区| 免费黄色在线免费观看| 中文欧美无线码| 成人亚洲精品av一区二区| 99久久精品热视频| 亚洲在线观看片| 久久精品综合一区二区三区| 国产伦理片在线播放av一区| 69av精品久久久久久| 欧美区成人在线视频| 亚洲国产精品成人久久小说| 亚洲精品日韩在线中文字幕| 成人二区视频| 久久99热这里只有精品18| 99热这里只有是精品50| 男女边吃奶边做爰视频| 亚洲性久久影院| 国产精品久久久久久av不卡| 国产免费一级a男人的天堂| 亚洲国产精品成人综合色| 黄色日韩在线| 男人和女人高潮做爰伦理| 亚洲激情五月婷婷啪啪| 美女高潮的动态| 少妇丰满av| 国产乱人偷精品视频| 一边亲一边摸免费视频| 亚洲人成网站在线观看播放| 狠狠狠狠99中文字幕| 最近中文字幕2019免费版| 久99久视频精品免费| 只有这里有精品99| 国产91av在线免费观看| 在线免费观看的www视频| 麻豆一二三区av精品| 超碰av人人做人人爽久久| 亚洲成色77777| 国产亚洲精品久久久com| 久久欧美精品欧美久久欧美| 成人亚洲精品av一区二区| 日韩三级伦理在线观看| 亚洲人成网站在线观看播放| 成人午夜高清在线视频| 欧美三级亚洲精品| 午夜激情福利司机影院| 午夜a级毛片| av线在线观看网站| 精品久久久久久久久av| 高清视频免费观看一区二区 | 国产精华一区二区三区| 久久久久久九九精品二区国产| 国产亚洲一区二区精品| 一边摸一边抽搐一进一小说| 人妻夜夜爽99麻豆av| 国产精品一二三区在线看| 久久婷婷人人爽人人干人人爱| 夜夜爽夜夜爽视频|