• <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
    一区在线观看完整版| 久久午夜福利片| 在线观看免费日韩欧美大片 | 午夜福利视频在线观看免费| 欧美xxⅹ黑人| 热re99久久精品国产66热6| 午夜av观看不卡| 少妇 在线观看| 热re99久久国产66热| 中国美白少妇内射xxxbb| tube8黄色片| 99久久综合免费| 69精品国产乱码久久久| 色5月婷婷丁香| 蜜桃国产av成人99| 天堂8中文在线网| 国产精品一国产av| 亚洲人成77777在线视频| av线在线观看网站| 亚洲精品久久午夜乱码| 夜夜看夜夜爽夜夜摸| 乱码一卡2卡4卡精品| 国产成人免费无遮挡视频| 少妇人妻久久综合中文| 超碰97精品在线观看| 亚洲欧美日韩卡通动漫| 欧美97在线视频| 五月玫瑰六月丁香| 成人手机av| 在线精品无人区一区二区三| 欧美日韩亚洲高清精品| 日本爱情动作片www.在线观看| 亚洲精品日韩av片在线观看| 国产亚洲精品第一综合不卡 | 精品久久久久久久久av| 欧美一级a爱片免费观看看| 亚洲色图综合在线观看| 国产精品国产三级专区第一集| 大码成人一级视频| 日韩熟女老妇一区二区性免费视频| 国产欧美日韩一区二区三区在线 | 97精品久久久久久久久久精品| 亚洲精品,欧美精品| 91精品三级在线观看| 久久人人爽人人片av| 欧美少妇被猛烈插入视频| 成人18禁高潮啪啪吃奶动态图 | 99视频精品全部免费 在线| av黄色大香蕉| 五月天丁香电影| 如日韩欧美国产精品一区二区三区 | 多毛熟女@视频| 这个男人来自地球电影免费观看 | 成人手机av| 色吧在线观看| 亚洲人成网站在线观看播放| 校园人妻丝袜中文字幕| 永久网站在线| 欧美丝袜亚洲另类| 国产亚洲最大av| 午夜久久久在线观看| 精品国产国语对白av| 国产极品粉嫩免费观看在线 | 国产伦理片在线播放av一区| 国产亚洲精品第一综合不卡 | 777米奇影视久久| 亚洲av日韩在线播放| 男女无遮挡免费网站观看| 91精品伊人久久大香线蕉| 婷婷色麻豆天堂久久| 夫妻性生交免费视频一级片| 乱码一卡2卡4卡精品| 亚洲成人手机| 亚洲第一区二区三区不卡| 精品视频人人做人人爽| 国产精品国产三级国产专区5o| 精品国产乱码久久久久久小说| kizo精华| 久久精品熟女亚洲av麻豆精品| 91久久精品国产一区二区成人| 成年人免费黄色播放视频| 国产免费一区二区三区四区乱码| 女人久久www免费人成看片| 午夜日本视频在线| 国产精品一区www在线观看| 国产成人精品无人区| 18禁动态无遮挡网站| 亚洲人与动物交配视频| 日韩强制内射视频| 永久免费av网站大全| 看十八女毛片水多多多| 黄色欧美视频在线观看| 国产乱来视频区| 一区二区三区免费毛片| 黄片无遮挡物在线观看| 中文天堂在线官网| 国产亚洲精品第一综合不卡 | 色吧在线观看| a级毛片在线看网站| 免费观看av网站的网址| av专区在线播放| 国产视频内射| 少妇猛男粗大的猛烈进出视频| 制服诱惑二区| 午夜福利网站1000一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一级,二级,三级黄色视频| 亚洲国产精品一区三区| 日韩强制内射视频| 亚洲精品aⅴ在线观看| 国产乱人偷精品视频| 亚洲av免费高清在线观看| 黄色毛片三级朝国网站| 国产精品一区二区在线不卡| 少妇被粗大猛烈的视频| 久久99一区二区三区| 日韩,欧美,国产一区二区三区| 久久久a久久爽久久v久久| 女的被弄到高潮叫床怎么办| 2022亚洲国产成人精品| 熟女av电影| av在线app专区| videossex国产| xxx大片免费视频| 69精品国产乱码久久久| 欧美精品一区二区免费开放| 亚洲av日韩在线播放| 老司机影院毛片| 国产一区二区三区综合在线观看 | 十分钟在线观看高清视频www| 美女脱内裤让男人舔精品视频| 18禁动态无遮挡网站| 一本久久精品| 2018国产大陆天天弄谢| 男女边摸边吃奶| 久久久久久久久久成人| 日韩电影二区| 内地一区二区视频在线| 亚洲国产精品成人久久小说| 国产黄片视频在线免费观看| 亚洲av中文av极速乱| 久久久久久久久久久免费av| 久久精品国产a三级三级三级| 国产av国产精品国产| 国产亚洲午夜精品一区二区久久| 色哟哟·www| 亚洲欧洲国产日韩| 婷婷色麻豆天堂久久| 中文字幕人妻丝袜制服| 一个人看视频在线观看www免费| 涩涩av久久男人的天堂| 国产亚洲精品第一综合不卡 | 美女内射精品一级片tv| 午夜免费鲁丝| 观看美女的网站| 一级二级三级毛片免费看| 亚洲色图 男人天堂 中文字幕 | 蜜桃国产av成人99| 国产在视频线精品| 日本色播在线视频| 一区二区av电影网| 久久精品国产亚洲av涩爱| 日韩成人伦理影院| 精品人妻一区二区三区麻豆| 国产精品久久久久久精品电影小说| 久久99热这里只频精品6学生| √禁漫天堂资源中文www| 国产成人精品婷婷| 大片电影免费在线观看免费| 永久网站在线| 狠狠精品人妻久久久久久综合| 久久精品久久精品一区二区三区| 亚洲,一卡二卡三卡| av卡一久久| 大香蕉97超碰在线| 大片电影免费在线观看免费| 亚洲国产欧美在线一区| 少妇的逼好多水| 97在线视频观看| xxxhd国产人妻xxx| 国产免费视频播放在线视频| 91精品伊人久久大香线蕉| 久久国内精品自在自线图片| 亚洲欧美清纯卡通| 校园人妻丝袜中文字幕| 91精品国产国语对白视频| av天堂久久9| 人人妻人人添人人爽欧美一区卜| 国产高清国产精品国产三级| 中文字幕最新亚洲高清| 精品卡一卡二卡四卡免费| 国产黄色视频一区二区在线观看| 亚洲精品,欧美精品| 午夜激情久久久久久久| 久久99蜜桃精品久久| 国产午夜精品久久久久久一区二区三区| videosex国产| 欧美激情极品国产一区二区三区 | 精品一区二区三区视频在线| 国产精品久久久久久av不卡| 热re99久久国产66热| 9色porny在线观看| 久久影院123| 99九九线精品视频在线观看视频| 国产av国产精品国产| 免费观看无遮挡的男女| 久久午夜福利片| 男人爽女人下面视频在线观看| 久久久久久久国产电影| 性色av一级| 一区二区日韩欧美中文字幕 | av电影中文网址| 亚洲美女视频黄频| 久久女婷五月综合色啪小说| 一边摸一边做爽爽视频免费| 波野结衣二区三区在线| 老司机影院毛片| 亚洲成人一二三区av| 欧美精品一区二区大全| 亚洲国产日韩一区二区| 91精品国产九色| 99热全是精品| 婷婷色综合www| 国模一区二区三区四区视频| 久久精品久久久久久噜噜老黄| 国产国语露脸激情在线看| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品一区三区| 久久国内精品自在自线图片| 少妇人妻精品综合一区二区| 97在线人人人人妻| 中文字幕人妻熟人妻熟丝袜美| 999精品在线视频| 欧美丝袜亚洲另类| 日韩免费高清中文字幕av| av电影中文网址| 51国产日韩欧美| 九九爱精品视频在线观看| 欧美97在线视频| 在线播放无遮挡| 人人妻人人澡人人爽人人夜夜| 人妻夜夜爽99麻豆av| 五月伊人婷婷丁香| 国产一区二区在线观看av| 各种免费的搞黄视频| 伦理电影免费视频| av有码第一页| 在线看a的网站| 69精品国产乱码久久久| 中文精品一卡2卡3卡4更新| 男女边吃奶边做爰视频| 高清不卡的av网站| 2021少妇久久久久久久久久久| 国产片特级美女逼逼视频| 特大巨黑吊av在线直播| 中文字幕久久专区| 国产男女超爽视频在线观看| 久久热精品热| 美女国产视频在线观看| 久久久久人妻精品一区果冻| 国产亚洲精品久久久com| 美女视频免费永久观看网站| 国产高清三级在线| 欧美激情国产日韩精品一区| 久久久国产精品麻豆| 男女边摸边吃奶| 各种免费的搞黄视频| 男男h啪啪无遮挡| 国产亚洲一区二区精品| 国产日韩欧美亚洲二区| 人妻 亚洲 视频| 欧美最新免费一区二区三区| 视频在线观看一区二区三区| 成人国语在线视频| 久久综合国产亚洲精品| 久久青草综合色| 国产精品女同一区二区软件| 亚洲高清免费不卡视频| 麻豆精品久久久久久蜜桃| 国产日韩欧美亚洲二区| 午夜视频国产福利| 日韩在线高清观看一区二区三区| 插逼视频在线观看| 视频在线观看一区二区三区| 亚洲国产av影院在线观看| 天堂8中文在线网| 永久免费av网站大全| 视频在线观看一区二区三区| 爱豆传媒免费全集在线观看| 王馨瑶露胸无遮挡在线观看| 在线免费观看不下载黄p国产| 免费黄频网站在线观看国产| av有码第一页| 国产精品无大码| 国产亚洲欧美精品永久| 免费不卡的大黄色大毛片视频在线观看| 高清av免费在线| 亚洲成人手机| 国产精品一区www在线观看| 成人二区视频| 国产av一区二区精品久久| 免费av不卡在线播放| 国产免费一级a男人的天堂| www.色视频.com| 在线观看美女被高潮喷水网站| 亚洲欧美一区二区三区国产| 丁香六月天网| 日韩在线高清观看一区二区三区| 一级毛片黄色毛片免费观看视频| 亚洲天堂av无毛| 91久久精品国产一区二区三区| 嫩草影院入口| 我要看黄色一级片免费的| 国产精品久久久久久精品电影小说| 亚洲国产精品国产精品| 色5月婷婷丁香| 永久免费av网站大全| 国产精品国产三级国产av玫瑰| 日韩成人av中文字幕在线观看| 最近的中文字幕免费完整| 日日啪夜夜爽| 国产精品99久久久久久久久| 国产女主播在线喷水免费视频网站| 最新中文字幕久久久久| 色视频在线一区二区三区| 免费黄色在线免费观看| 大又大粗又爽又黄少妇毛片口| 高清午夜精品一区二区三区| 亚洲欧美成人综合另类久久久| 制服人妻中文乱码| 免费观看的影片在线观看| 如何舔出高潮| 久久亚洲国产成人精品v| 满18在线观看网站| 极品人妻少妇av视频| 成人漫画全彩无遮挡| 美女国产高潮福利片在线看| 高清毛片免费看| 国产一区有黄有色的免费视频| 热99久久久久精品小说推荐| 免费大片黄手机在线观看| 18禁观看日本| 亚洲国产最新在线播放| 99久久综合免费| 亚洲综合精品二区| 欧美日韩亚洲高清精品| 日日啪夜夜爽| 国产又色又爽无遮挡免| 欧美日韩在线观看h| 国产精品人妻久久久影院| 久久人人爽人人爽人人片va| 免费大片18禁| 最近2019中文字幕mv第一页| 秋霞伦理黄片| 成年av动漫网址| 国产亚洲av片在线观看秒播厂| 妹子高潮喷水视频| 亚洲婷婷狠狠爱综合网| a 毛片基地| 九九久久精品国产亚洲av麻豆| 亚洲国产av新网站| 爱豆传媒免费全集在线观看| 国产男女内射视频| a级毛片免费高清观看在线播放| av不卡在线播放| 成年人午夜在线观看视频| 狂野欧美白嫩少妇大欣赏| av专区在线播放| 不卡视频在线观看欧美| 18+在线观看网站| 国产有黄有色有爽视频| 免费观看av网站的网址| 日韩不卡一区二区三区视频在线| 黄色怎么调成土黄色| av国产精品久久久久影院| 久久久久人妻精品一区果冻| 99久国产av精品国产电影| 在线观看免费高清a一片| 国产精品久久久久久精品电影小说| 久久精品久久久久久久性| 国产精品嫩草影院av在线观看| a级毛片免费高清观看在线播放| 少妇丰满av| 亚洲少妇的诱惑av| 岛国毛片在线播放| 在线观看免费视频网站a站| 老熟女久久久| 丝袜在线中文字幕| 国产深夜福利视频在线观看| 免费少妇av软件| 久久人妻熟女aⅴ| 人妻夜夜爽99麻豆av| 欧美人与性动交α欧美精品济南到 | 成人午夜精彩视频在线观看| av黄色大香蕉| 国产精品久久久久久久久免| www.av在线官网国产| 91午夜精品亚洲一区二区三区| av免费在线看不卡| 日韩强制内射视频| 欧美日韩国产mv在线观看视频| 99热这里只有精品一区| 久久毛片免费看一区二区三区| 91aial.com中文字幕在线观看| 亚洲av男天堂| 亚洲精品日韩在线中文字幕| 婷婷色综合大香蕉| 国产伦精品一区二区三区视频9| 中文精品一卡2卡3卡4更新| 男男h啪啪无遮挡| 简卡轻食公司| 少妇精品久久久久久久| 国产无遮挡羞羞视频在线观看| 性色av一级| 啦啦啦啦在线视频资源| 男人添女人高潮全过程视频| 伦理电影大哥的女人| 亚洲av二区三区四区| av网站免费在线观看视频| 久久精品人人爽人人爽视色| 欧美激情国产日韩精品一区| 国产女主播在线喷水免费视频网站| 一边亲一边摸免费视频| 亚洲av电影在线观看一区二区三区| 秋霞在线观看毛片| 精品一区二区三卡| 国产精品偷伦视频观看了| 国产男人的电影天堂91| a 毛片基地| 男人添女人高潮全过程视频| 精品人妻一区二区三区麻豆| 亚洲av日韩在线播放| 极品人妻少妇av视频| 亚洲av综合色区一区| 一本一本综合久久| tube8黄色片| 尾随美女入室| 高清午夜精品一区二区三区| 精品国产一区二区久久| 男人操女人黄网站| 搡老乐熟女国产| 人妻少妇偷人精品九色| 欧美人与善性xxx| 亚洲熟女精品中文字幕| 高清av免费在线| 精品久久久精品久久久| 男人添女人高潮全过程视频| 亚洲av日韩在线播放| 青春草亚洲视频在线观看| 男人爽女人下面视频在线观看| 18+在线观看网站| 中文天堂在线官网| 精品久久久久久久久亚洲| 18禁在线播放成人免费| 午夜激情福利司机影院| 国产毛片在线视频| 久久免费观看电影| 桃花免费在线播放| 99热全是精品| 久久毛片免费看一区二区三区| 两个人免费观看高清视频| 亚洲国产毛片av蜜桃av| 亚洲,欧美,日韩| 国产精品一国产av| 精品少妇久久久久久888优播| 国产日韩欧美亚洲二区| 老司机影院毛片| 激情五月婷婷亚洲| a级毛片黄视频| 久久久国产欧美日韩av| 国产黄色免费在线视频| 国产视频首页在线观看| 国产色爽女视频免费观看| 亚洲人成77777在线视频| 少妇的逼好多水| 男人爽女人下面视频在线观看| 国产色婷婷99| 日韩伦理黄色片| 91精品国产国语对白视频| 亚洲人成网站在线播| 人妻 亚洲 视频| 久久久久国产网址| 亚洲色图 男人天堂 中文字幕 | 狂野欧美白嫩少妇大欣赏| 内地一区二区视频在线| 99国产精品免费福利视频| 丝袜脚勾引网站| 乱码一卡2卡4卡精品| 亚洲人成77777在线视频| 18禁裸乳无遮挡动漫免费视频| 久久久久久人妻| 免费观看在线日韩| 人人妻人人添人人爽欧美一区卜| 如何舔出高潮| 国产男女超爽视频在线观看| 久久精品久久精品一区二区三区| 亚洲人成网站在线播| 国产午夜精品一二区理论片| 久久女婷五月综合色啪小说| 菩萨蛮人人尽说江南好唐韦庄| 插阴视频在线观看视频| 纯流量卡能插随身wifi吗| 插逼视频在线观看| 亚洲国产av影院在线观看| 建设人人有责人人尽责人人享有的| 国产精品 国内视频| 黄色怎么调成土黄色| 在线观看免费视频网站a站| 国产免费一级a男人的天堂| 国产成人a∨麻豆精品| 蜜桃久久精品国产亚洲av| 国产成人精品在线电影| 亚洲精品乱久久久久久| 最近手机中文字幕大全| 九九爱精品视频在线观看| 国产高清不卡午夜福利| 成人漫画全彩无遮挡| 国产精品蜜桃在线观看| 成人毛片a级毛片在线播放| 国产精品女同一区二区软件| 免费少妇av软件| 亚洲av欧美aⅴ国产| 丰满乱子伦码专区| 国产乱人偷精品视频| 在线天堂最新版资源| 久久青草综合色| 久久精品国产亚洲av天美| 男的添女的下面高潮视频| 亚洲国产精品一区三区| 免费少妇av软件| 这个男人来自地球电影免费观看 | 欧美日本中文国产一区发布| 精品一区二区免费观看| 能在线免费看毛片的网站| 欧美激情国产日韩精品一区| 亚洲美女视频黄频| 免费看av在线观看网站| 男人添女人高潮全过程视频| 最黄视频免费看| 日本av手机在线免费观看| 免费人成在线观看视频色| 午夜日本视频在线| 国产精品熟女久久久久浪| 成人午夜精彩视频在线观看| 少妇人妻久久综合中文| 亚洲欧美中文字幕日韩二区| 九九在线视频观看精品| av又黄又爽大尺度在线免费看| 美女国产高潮福利片在线看| 免费av不卡在线播放| 2018国产大陆天天弄谢| 中文字幕人妻熟人妻熟丝袜美| 一级毛片黄色毛片免费观看视频| av在线播放精品| 亚洲av成人精品一区久久| 成人国产麻豆网| 亚洲国产精品成人久久小说| 成人亚洲精品一区在线观看| 亚州av有码| 国产欧美日韩一区二区三区在线 | 一级a做视频免费观看| 国产免费视频播放在线视频| 国产深夜福利视频在线观看| 91久久精品电影网| 2021少妇久久久久久久久久久| 一个人免费看片子| 日韩不卡一区二区三区视频在线| 亚洲精品日本国产第一区| 99国产综合亚洲精品| 国产片内射在线| 成年女人在线观看亚洲视频| 久久99蜜桃精品久久| 国产精品一区www在线观看| 免费看光身美女| 久久精品久久精品一区二区三区| 午夜视频国产福利| 少妇人妻 视频| 中文字幕免费在线视频6| 亚洲一区二区三区欧美精品| 热re99久久国产66热| 青春草国产在线视频| 在线天堂最新版资源| av电影中文网址| 赤兔流量卡办理| 精品一区在线观看国产| 亚洲精品久久成人aⅴ小说 | 九草在线视频观看| 99久久综合免费| 伦精品一区二区三区| 国产黄片视频在线免费观看| 晚上一个人看的免费电影| 国产精品偷伦视频观看了| 久久午夜综合久久蜜桃| 考比视频在线观看| 26uuu在线亚洲综合色| 亚洲精品中文字幕在线视频| 久久久亚洲精品成人影院| 亚洲精品国产色婷婷电影| 久久久欧美国产精品| 国产亚洲午夜精品一区二区久久| 久久久久久久久久久丰满| 亚洲国产av影院在线观看| 国产精品久久久久久精品电影小说| 精品熟女少妇av免费看| 在线观看三级黄色| 日韩电影二区| 亚洲内射少妇av| 午夜91福利影院| 久久精品国产鲁丝片午夜精品| 成人国产麻豆网| 18禁在线无遮挡免费观看视频| 欧美激情极品国产一区二区三区 | 大又大粗又爽又黄少妇毛片口| xxx大片免费视频| 久久久久久久大尺度免费视频|