• <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
    日韩欧美一区二区三区在线观看| 国产精品久久电影中文字幕| 国产精品乱码一区二三区的特点| 亚洲真实伦在线观看| 日韩欧美三级三区| 黄片大片在线免费观看| 亚洲天堂国产精品一区在线| 变态另类成人亚洲欧美熟女| 国产蜜桃级精品一区二区三区| 精品免费久久久久久久清纯| 免费看日本二区| 一级毛片女人18水好多| 国产极品精品免费视频能看的| 他把我摸到了高潮在线观看| 国产精华一区二区三区| 亚洲在线观看片| 国产精品国产高清国产av| 99国产极品粉嫩在线观看| 亚洲中文av在线| 黄片大片在线免费观看| 淫妇啪啪啪对白视频| 亚洲精品一卡2卡三卡4卡5卡| 97人妻精品一区二区三区麻豆| 国产精品爽爽va在线观看网站| 99在线视频只有这里精品首页| 成人亚洲精品av一区二区| 成人18禁在线播放| 亚洲专区国产一区二区| 精品国产三级普通话版| 91久久精品国产一区二区成人 | 青草久久国产| 亚洲在线观看片| 亚洲精品在线美女| 哪里可以看免费的av片| 成人性生交大片免费视频hd| netflix在线观看网站| 亚洲 国产 在线| 成年女人永久免费观看视频| 97超级碰碰碰精品色视频在线观看| 日日夜夜操网爽| 久久精品91蜜桃| 国产精品,欧美在线| 免费看十八禁软件| 啦啦啦观看免费观看视频高清| 午夜福利在线观看免费完整高清在 | 国产精品野战在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲九九香蕉| 午夜福利成人在线免费观看| 变态另类成人亚洲欧美熟女| 长腿黑丝高跟| x7x7x7水蜜桃| 国产成人精品久久二区二区免费| 婷婷精品国产亚洲av在线| 欧美一区二区国产精品久久精品| 成人鲁丝片一二三区免费| 香蕉久久夜色| 国产精品香港三级国产av潘金莲| 又大又爽又粗| 可以在线观看毛片的网站| 两个人的视频大全免费| a级毛片a级免费在线| 一进一出抽搐动态| 不卡av一区二区三区| 日韩欧美精品v在线| 免费看美女性在线毛片视频| 欧美中文综合在线视频| 国产亚洲av高清不卡| 色综合亚洲欧美另类图片| 国产精品久久视频播放| 中文在线观看免费www的网站| 在线观看66精品国产| 热99在线观看视频| 国产亚洲精品久久久com| 黑人欧美特级aaaaaa片| 久久久久久人人人人人| 偷拍熟女少妇极品色| 国产97色在线日韩免费| 精品福利观看| 久久天堂一区二区三区四区| 禁无遮挡网站| 日日夜夜操网爽| www.自偷自拍.com| 国产黄色小视频在线观看| 免费看日本二区| 成人国产综合亚洲| 美女黄网站色视频| 免费在线观看成人毛片| 亚洲无线观看免费| 午夜福利欧美成人| 十八禁人妻一区二区| 九色国产91popny在线| 一区二区三区激情视频| 欧美一区二区国产精品久久精品| 成人亚洲精品av一区二区| 国产精品98久久久久久宅男小说| 国产伦精品一区二区三区视频9 | av片东京热男人的天堂| 日本熟妇午夜| 色吧在线观看| 欧美色视频一区免费| 99久久精品一区二区三区| 男人舔女人的私密视频| av视频在线观看入口| 亚洲男人的天堂狠狠| 日本成人三级电影网站| 国产精品九九99| 美女 人体艺术 gogo| 韩国av一区二区三区四区| 亚洲欧美日韩无卡精品| 日韩高清综合在线| 人人妻人人看人人澡| 亚洲成人精品中文字幕电影| 黄片大片在线免费观看| 中文在线观看免费www的网站| 又黄又粗又硬又大视频| 亚洲av五月六月丁香网| 精品久久久久久久久久久久久| 国产熟女xx| 日韩欧美精品v在线| av在线天堂中文字幕| 99久久精品国产亚洲精品| 国产亚洲精品综合一区在线观看| 在线a可以看的网站| 亚洲成人中文字幕在线播放| 露出奶头的视频| 99在线视频只有这里精品首页| 国产精品电影一区二区三区| 午夜视频精品福利| 日本精品一区二区三区蜜桃| 18禁裸乳无遮挡免费网站照片| 精品无人区乱码1区二区| 可以在线观看毛片的网站| 三级国产精品欧美在线观看 | 成年女人毛片免费观看观看9| 在线国产一区二区在线| av在线蜜桃| 国产野战对白在线观看| 欧美黄色片欧美黄色片| 美女高潮喷水抽搐中文字幕| 性色av乱码一区二区三区2| 老司机午夜福利在线观看视频| 韩国av一区二区三区四区| 久久久久久九九精品二区国产| 亚洲国产精品sss在线观看| 波多野结衣高清作品| 叶爱在线成人免费视频播放| 国产精品 国内视频| 午夜激情福利司机影院| 一进一出抽搐动态| 亚洲国产精品久久男人天堂| 黑人操中国人逼视频| 一二三四在线观看免费中文在| 我的老师免费观看完整版| 亚洲avbb在线观看| 久久伊人香网站| 欧美精品啪啪一区二区三区| 午夜免费成人在线视频| 在线免费观看的www视频| 亚洲欧美一区二区三区黑人| 国产成人欧美在线观看| 手机成人av网站| 免费av毛片视频| 一边摸一边抽搐一进一小说| 日本撒尿小便嘘嘘汇集6| 精品一区二区三区av网在线观看| 欧美大码av| 白带黄色成豆腐渣| 亚洲欧美日韩无卡精品| 亚洲熟妇熟女久久| av视频在线观看入口| 欧美另类亚洲清纯唯美| 精品电影一区二区在线| 日日干狠狠操夜夜爽| 老司机午夜十八禁免费视频| 岛国在线观看网站| 久久香蕉国产精品| 法律面前人人平等表现在哪些方面| 精品国产超薄肉色丝袜足j| 亚洲狠狠婷婷综合久久图片| 欧美大码av| 男人舔奶头视频| 亚洲aⅴ乱码一区二区在线播放| 观看免费一级毛片| 日韩成人在线观看一区二区三区| 欧美成狂野欧美在线观看| 在线国产一区二区在线| 精品乱码久久久久久99久播| 婷婷丁香在线五月| 此物有八面人人有两片| 天天躁日日操中文字幕| 久久久久免费精品人妻一区二区| 欧美黄色片欧美黄色片| 欧美日韩亚洲国产一区二区在线观看| 亚洲专区中文字幕在线| 丝袜人妻中文字幕| 搡老熟女国产l中国老女人| 婷婷六月久久综合丁香| 久久性视频一级片| 国产熟女xx| 香蕉av资源在线| 毛片女人毛片| 欧美中文综合在线视频| 亚洲真实伦在线观看| 欧美午夜高清在线| 夜夜夜夜夜久久久久| 一个人看视频在线观看www免费 | 无遮挡黄片免费观看| 午夜免费观看网址| 97超级碰碰碰精品色视频在线观看| 又紧又爽又黄一区二区| 亚洲av电影不卡..在线观看| 午夜日韩欧美国产| 久久久国产欧美日韩av| 成年人黄色毛片网站| 在线国产一区二区在线| 波多野结衣高清作品| 久久久久国产一级毛片高清牌| 老熟妇乱子伦视频在线观看| 在线观看美女被高潮喷水网站 | 国模一区二区三区四区视频 | 欧美xxxx黑人xx丫x性爽| 成人精品一区二区免费| 亚洲18禁久久av| 欧美3d第一页| 男人舔奶头视频| 99riav亚洲国产免费| 欧美乱码精品一区二区三区| 99在线人妻在线中文字幕| 在线看三级毛片| 久久精品91蜜桃| 亚洲欧美一区二区三区黑人| 黑人操中国人逼视频| 美女黄网站色视频| 亚洲精品456在线播放app | 性欧美人与动物交配| 国内精品美女久久久久久| 国产成人精品久久二区二区91| 亚洲在线观看片| 一a级毛片在线观看| 国产精品美女特级片免费视频播放器 | 岛国视频午夜一区免费看| 99热6这里只有精品| 国产乱人伦免费视频| 露出奶头的视频| 最近最新中文字幕大全电影3| 波多野结衣高清无吗| www.www免费av| 熟妇人妻久久中文字幕3abv| 国产91精品成人一区二区三区| 首页视频小说图片口味搜索| 999精品在线视频| 制服丝袜大香蕉在线| 欧美最黄视频在线播放免费| 欧美午夜高清在线| 久久精品国产99精品国产亚洲性色| 色哟哟哟哟哟哟| 亚洲精品中文字幕一二三四区| 国产精品电影一区二区三区| 欧美黑人欧美精品刺激| 老司机午夜福利在线观看视频| 99久久久亚洲精品蜜臀av| 久久久久久久精品吃奶| 午夜日韩欧美国产| 国产亚洲欧美在线一区二区| 国内精品久久久久精免费| 久久这里只有精品中国| 久久草成人影院| 黄色 视频免费看| 91久久精品国产一区二区成人 | 国产黄片美女视频| 俄罗斯特黄特色一大片| 国产精品电影一区二区三区| 亚洲精品美女久久久久99蜜臀| 亚洲av电影不卡..在线观看| 美女大奶头视频| 18禁黄网站禁片午夜丰满| 全区人妻精品视频| 中国美女看黄片| 天天一区二区日本电影三级| 在线观看美女被高潮喷水网站 | 国产日本99.免费观看| 久久精品91蜜桃| 一个人观看的视频www高清免费观看 | 国产野战对白在线观看| 免费在线观看亚洲国产| 国产成人影院久久av| 日韩精品中文字幕看吧| 精品国产乱码久久久久久男人| 亚洲 欧美一区二区三区| 1024香蕉在线观看| 国产成人啪精品午夜网站| 亚洲欧美日韩无卡精品| 欧美3d第一页| 午夜福利18| 别揉我奶头~嗯~啊~动态视频| 久久香蕉精品热| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产日韩欧美精品在线观看 | 日韩有码中文字幕| 一二三四在线观看免费中文在| 国产黄色小视频在线观看| 在线国产一区二区在线| 90打野战视频偷拍视频| 国产av麻豆久久久久久久| 少妇丰满av| 亚洲一区高清亚洲精品| 亚洲精品一区av在线观看| 成人性生交大片免费视频hd| 国产精品综合久久久久久久免费| 午夜福利视频1000在线观看| h日本视频在线播放| 51午夜福利影视在线观看| 91九色精品人成在线观看| 亚洲精品美女久久av网站| 国产人伦9x9x在线观看| 国产伦精品一区二区三区视频9 | 丰满人妻熟妇乱又伦精品不卡| 美女免费视频网站| 久久午夜亚洲精品久久| 身体一侧抽搐| 18禁裸乳无遮挡免费网站照片| 免费无遮挡裸体视频| а√天堂www在线а√下载| 国产亚洲精品av在线| 国产精品九九99| 午夜精品一区二区三区免费看| 国产欧美日韩一区二区三| 国产成人啪精品午夜网站| 精华霜和精华液先用哪个| 久久久久性生活片| 国产三级在线视频| 首页视频小说图片口味搜索| 亚洲国产精品999在线| 精品一区二区三区四区五区乱码| 好男人电影高清在线观看| 91久久精品国产一区二区成人 | 国产99白浆流出| 免费在线观看影片大全网站| 日韩精品青青久久久久久| 久久久久久国产a免费观看| 在线观看免费午夜福利视频| 国产成人精品久久二区二区91| 国产高清三级在线| 亚洲 欧美 日韩 在线 免费| 国产精品免费一区二区三区在线| 在线观看免费午夜福利视频| 超碰成人久久| 国产精品久久久久久精品电影| 欧美色视频一区免费| 看黄色毛片网站| 少妇丰满av| 成年人黄色毛片网站| 韩国av一区二区三区四区| 午夜福利欧美成人| 给我免费播放毛片高清在线观看| 黑人操中国人逼视频| 热99在线观看视频| 日本一二三区视频观看| 人妻丰满熟妇av一区二区三区| 亚洲成人精品中文字幕电影| 亚洲国产欧美一区二区综合| 成人欧美大片| 亚洲真实伦在线观看| 亚洲五月天丁香| 国产av麻豆久久久久久久| 久久亚洲精品不卡| 精品无人区乱码1区二区| 国产亚洲精品久久久久久毛片| 最近视频中文字幕2019在线8| 国产av麻豆久久久久久久| 亚洲精品粉嫩美女一区| 亚洲中文日韩欧美视频| 精品国产乱码久久久久久男人| 成人永久免费在线观看视频| 亚洲av熟女| 欧美不卡视频在线免费观看| 这个男人来自地球电影免费观看| 国产99白浆流出| 亚洲欧美日韩高清专用| 18禁观看日本| 午夜激情欧美在线| 午夜免费成人在线视频| 99久国产av精品| 97人妻精品一区二区三区麻豆| 成年女人毛片免费观看观看9| 美女免费视频网站| 久久伊人香网站| 草草在线视频免费看| 熟女少妇亚洲综合色aaa.| 草草在线视频免费看| 国产av一区在线观看免费| 男女下面进入的视频免费午夜| 婷婷丁香在线五月| 国产伦精品一区二区三区视频9 | 午夜精品久久久久久毛片777| 一夜夜www| 99久国产av精品| 麻豆久久精品国产亚洲av| 国产一区二区在线av高清观看| 在线观看免费午夜福利视频| 变态另类成人亚洲欧美熟女| 天天添夜夜摸| 真实男女啪啪啪动态图| 国产真人三级小视频在线观看| 国产美女午夜福利| 小说图片视频综合网站| 最近视频中文字幕2019在线8| 国产精品一区二区免费欧美| 一a级毛片在线观看| 国产视频内射| 12—13女人毛片做爰片一| 桃红色精品国产亚洲av| 欧美激情久久久久久爽电影| 国产一区二区激情短视频| 制服人妻中文乱码| 欧美av亚洲av综合av国产av| 午夜激情福利司机影院| 国产一区二区三区视频了| 亚洲五月天丁香| 亚洲欧美日韩东京热| 亚洲九九香蕉| 啦啦啦韩国在线观看视频| 免费电影在线观看免费观看| 丁香六月欧美| 欧美性猛交╳xxx乱大交人| 搡老熟女国产l中国老女人| 日日干狠狠操夜夜爽| 一二三四在线观看免费中文在| 国产精华一区二区三区| 又粗又爽又猛毛片免费看| 最新在线观看一区二区三区| av欧美777| 最新美女视频免费是黄的| 哪里可以看免费的av片| netflix在线观看网站| 精品免费久久久久久久清纯| 国产精品亚洲一级av第二区| 日韩有码中文字幕| 少妇人妻一区二区三区视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩综合久久久久久 | 一卡2卡三卡四卡精品乱码亚洲| 日本 av在线| 熟女少妇亚洲综合色aaa.| 免费av毛片视频| 看免费av毛片| 俄罗斯特黄特色一大片| 黄色成人免费大全| 久久午夜亚洲精品久久| 最好的美女福利视频网| 色尼玛亚洲综合影院| www.999成人在线观看| 看片在线看免费视频| 少妇裸体淫交视频免费看高清| 国内精品一区二区在线观看| 日本一二三区视频观看| 波多野结衣高清无吗| 色综合欧美亚洲国产小说| 黄色女人牲交| 人妻丰满熟妇av一区二区三区| 久久精品影院6| 中文资源天堂在线| 老鸭窝网址在线观看| 亚洲 欧美 日韩 在线 免费| 免费在线观看成人毛片| 国产不卡一卡二| 久久人妻av系列| 琪琪午夜伦伦电影理论片6080| 欧美绝顶高潮抽搐喷水| 欧美成人一区二区免费高清观看 | 成人鲁丝片一二三区免费| 精品国产三级普通话版| 免费在线观看日本一区| 视频区欧美日本亚洲| 少妇熟女aⅴ在线视频| 高清毛片免费观看视频网站| 999久久久国产精品视频| 亚洲第一欧美日韩一区二区三区| 免费av毛片视频| 99久久精品一区二区三区| 一夜夜www| 女同久久另类99精品国产91| 熟女人妻精品中文字幕| 午夜福利欧美成人| 精品国产乱子伦一区二区三区| 叶爱在线成人免费视频播放| 在线观看免费视频日本深夜| 中文字幕精品亚洲无线码一区| 91av网一区二区| 久久人妻av系列| 噜噜噜噜噜久久久久久91| 白带黄色成豆腐渣| 99精品久久久久人妻精品| 成人三级黄色视频| 九色国产91popny在线| or卡值多少钱| 日韩精品中文字幕看吧| 精品一区二区三区视频在线 | 性色av乱码一区二区三区2| 久久香蕉精品热| 精品一区二区三区四区五区乱码| 久久精品国产亚洲av香蕉五月| 日日夜夜操网爽| а√天堂www在线а√下载| 成年女人看的毛片在线观看| 国产69精品久久久久777片 | 中文资源天堂在线| 免费看十八禁软件| 午夜福利在线在线| 在线播放国产精品三级| 成人特级av手机在线观看| 老熟妇乱子伦视频在线观看| 日韩有码中文字幕| 国内揄拍国产精品人妻在线| 一级毛片精品| 久久亚洲真实| ponron亚洲| 熟女少妇亚洲综合色aaa.| 亚洲熟妇熟女久久| 欧美不卡视频在线免费观看| 欧美日韩福利视频一区二区| 一a级毛片在线观看| 99精品在免费线老司机午夜| 免费看美女性在线毛片视频| 亚洲 欧美 日韩 在线 免费| 欧美丝袜亚洲另类 | 精品久久久久久久人妻蜜臀av| 亚洲aⅴ乱码一区二区在线播放| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品456在线播放app | 老司机深夜福利视频在线观看| 这个男人来自地球电影免费观看| av中文乱码字幕在线| 中文字幕最新亚洲高清| 叶爱在线成人免费视频播放| 久久精品91蜜桃| 国产成人欧美在线观看| 一二三四社区在线视频社区8| 亚洲最大成人中文| 成人av在线播放网站| 国产精品香港三级国产av潘金莲| 真人做人爱边吃奶动态| 国产精品久久视频播放| 国产av在哪里看| 黄色丝袜av网址大全| 天天添夜夜摸| 欧美三级亚洲精品| 国产v大片淫在线免费观看| 老汉色av国产亚洲站长工具| 无限看片的www在线观看| 成年人黄色毛片网站| 一级a爱片免费观看的视频| 人妻丰满熟妇av一区二区三区| 两个人视频免费观看高清| 一进一出抽搐gif免费好疼| 精品日产1卡2卡| 国产三级黄色录像| 婷婷精品国产亚洲av| 国产麻豆成人av免费视频| 色综合亚洲欧美另类图片| 91九色精品人成在线观看| 国产亚洲av嫩草精品影院| 亚洲成人中文字幕在线播放| 一级毛片精品| 免费看十八禁软件| 亚洲欧美日韩无卡精品| 亚洲 欧美一区二区三区| 欧美黑人巨大hd| 婷婷精品国产亚洲av在线| 男女那种视频在线观看| 制服丝袜大香蕉在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产精品久久男人天堂| 极品教师在线免费播放| 身体一侧抽搐| 国产成人啪精品午夜网站| 舔av片在线| 国产一级毛片七仙女欲春2| 精品午夜福利视频在线观看一区| 亚洲精品美女久久av网站| 黄色日韩在线| 国产69精品久久久久777片 | 少妇熟女aⅴ在线视频| 国产成人aa在线观看| 久久精品国产清高在天天线| 国产午夜精品论理片| 精品一区二区三区视频在线 | 国模一区二区三区四区视频 | 久久久久久久久免费视频了| 亚洲欧美日韩卡通动漫| 波多野结衣高清无吗| 亚洲人成网站高清观看| 69av精品久久久久久| 国产av在哪里看| 啦啦啦观看免费观看视频高清| 亚洲国产欧美一区二区综合| 久久精品人妻少妇| 性欧美人与动物交配| 久久久成人免费电影| 国产探花在线观看一区二区| 国产精品精品国产色婷婷| 欧美日韩瑟瑟在线播放| 18美女黄网站色大片免费观看| 夜夜看夜夜爽夜夜摸| 狂野欧美白嫩少妇大欣赏| 国产成人系列免费观看| www.熟女人妻精品国产| 亚洲精品美女久久久久99蜜臀| netflix在线观看网站| 99国产精品一区二区蜜桃av| 动漫黄色视频在线观看| 亚洲欧美激情综合另类| 男人舔女人下体高潮全视频| 啦啦啦观看免费观看视频高清|