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

    利用RTK進(jìn)行多波束水深測(cè)量時(shí)GPS高程數(shù)據(jù)的處理

    2015-08-18 01:45:22劉福剛
    山東國(guó)土資源 2015年9期
    關(guān)鍵詞:潮位橢球換能器

    劉福剛

    (中石化石油工程設(shè)計(jì)有限公司,山東 東營(yíng) 257000)

    ?

    利用RTK進(jìn)行多波束水深測(cè)量時(shí)GPS高程數(shù)據(jù)的處理

    劉福剛

    (中石化石油工程設(shè)計(jì)有限公司,山東 東營(yíng)257000)

    利用RTK進(jìn)行多波束水深測(cè)量時(shí),GPS高程在提供瞬時(shí)垂直基準(zhǔn)方面起著重要的作用。海上實(shí)際作業(yè)過程中,GPS信號(hào)經(jīng)常會(huì)發(fā)生失鎖現(xiàn)象,造成部分歷元高程數(shù)據(jù)突變,導(dǎo)致海底地形失真。該文以茂名某項(xiàng)目為例,首先介紹了根據(jù)EGM2008大地水準(zhǔn)面模型,利用高程擬合法將多波束水深測(cè)量時(shí)GPS橢球高轉(zhuǎn)換為正常高的方法;然后針對(duì)GPS信號(hào)失鎖現(xiàn)象,提出了利用RTK高程中的固定解來求算海水潮位面的方法,同時(shí)闡述了該方法的優(yōu)勢(shì)。

    RTK;水深測(cè)量;多波束;橢球高;正常高;高程數(shù)據(jù)處理

    引文格式:劉福剛.利用RTK進(jìn)行多波束水深測(cè)量時(shí)GPS高程數(shù)據(jù)的處理[J].山東國(guó)土資源,2015,31(9):90-93.LIU Fugang.Elevation Data Processing of GPS in Multibeam Bathymetric Measuring by Using RTK Technology[J].Shandong Land and Resources,2015,31(9):90-93.

    0 引言

    多波束測(cè)深系統(tǒng)是水聲、計(jì)算機(jī)、導(dǎo)航定位和數(shù)字化傳感器等多種技術(shù)的高度集成。是利用安裝于船底或拖體上的聲基陣向與航向垂直的海底發(fā)射超寬聲波束,接收海底反向散射信號(hào),經(jīng)過模擬/數(shù)字信號(hào)處理,形成多個(gè)波束,同時(shí)獲得幾十個(gè)甚至上百個(gè)海底條帶上采樣點(diǎn)的水深數(shù)據(jù)。其測(cè)量條帶覆蓋范圍為水深的2~10倍,與現(xiàn)場(chǎng)采集的導(dǎo)航定位及姿態(tài)數(shù)據(jù)相結(jié)合,能夠繪制出高精度、高分辨率的數(shù)字成果圖。與單波束回聲測(cè)深儀相比,多波束測(cè)深系統(tǒng)具有測(cè)量范圍大、測(cè)量速度快、精度高的優(yōu)點(diǎn)。

    對(duì)發(fā)射聲波瞬間換能器的定位、定姿是多波束水深測(cè)量的重要一步。目前,可利用差分GPS(DGPS)系統(tǒng)或RTK的方式求出換能器平面位置,利用電羅經(jīng)、姿態(tài)傳感器求出換能器瞬時(shí)姿態(tài),定位定姿中的導(dǎo)航延時(shí)、橫搖、縱搖、艏偏通過檢校場(chǎng)檢校計(jì)算得出。傳統(tǒng)多波束測(cè)量中,換能器瞬時(shí)高程是由潮位、吃水參數(shù)和涌浪參數(shù)(Heave)聯(lián)合提供,理想條件下,該高程是準(zhǔn)確無誤的。但是實(shí)際測(cè)量過程中,潮位是由測(cè)區(qū)多個(gè)驗(yàn)潮站數(shù)據(jù)線性內(nèi)插得到,吃水參數(shù)易受船速等因素影響,Heave在船速突變時(shí)會(huì)有長(zhǎng)周期信號(hào)異常,使得計(jì)算出的換能器瞬時(shí)高程值精度大受影響[1]。

    GPS-RTK定位技術(shù)是基于載波相位觀測(cè)值的實(shí)時(shí)動(dòng)態(tài)定位技術(shù)。在RTK解固定時(shí),通過量測(cè)、計(jì)算GPS天線相位中心與換能器間的相對(duì)關(guān)系,能實(shí)時(shí)求出換能器瞬時(shí)高程。

    在已知GPS天線高程的情況下,只需量測(cè)GPS天線位置到換能器之間的高差值h1,便可求出換能器瞬時(shí)高程,從而消除了潮位模型選擇、動(dòng)態(tài)吃水測(cè)量不準(zhǔn)確及船速突變時(shí)Heave異常帶來的影響(圖1)。

    圖1 GPS接收機(jī)天線與換能器相對(duì)位置示意圖

    1 工程應(yīng)用

    茂名某海洋測(cè)量項(xiàng)目,由岸邊往海中東南方向延伸約15km,水深20m左右,需要測(cè)量管線兩側(cè)各250m內(nèi)水下地形。項(xiàng)目采用WGS-84 坐標(biāo)系統(tǒng)(UTM投影),1985國(guó)家高程基準(zhǔn)。該次采用設(shè)備有天寶GPS 2臺(tái)、Geo Swath Plus多波束測(cè)深系統(tǒng)(250Hz)和處理軟件GS+。

    GPS測(cè)得的是點(diǎn)位的橢球高,工程需要的是該點(diǎn)的正常高,兩者間的差值(高程異常)不是一個(gè)常數(shù)。同時(shí)海上作業(yè)過程中經(jīng)常出現(xiàn)失鎖現(xiàn)象,GPS信號(hào)失鎖部分的高程無法利用。以該項(xiàng)目為例介紹橢球高到正常高的轉(zhuǎn)換方法及利用GPS高程計(jì)算潮位數(shù)據(jù)的方法來消除GPS信號(hào)失鎖的影響。

    1.1GPS橢球高到正常高的轉(zhuǎn)換

    橢球高到正常高的轉(zhuǎn)換,一般有2種方法:

    1.1.1七參數(shù)法

    測(cè)區(qū)范圍內(nèi)取3個(gè)以上已知控制點(diǎn),通過輸入已知點(diǎn)的WGS84坐標(biāo)及當(dāng)?shù)刈鴺?biāo),求出3個(gè)平移參數(shù)、3個(gè)旋轉(zhuǎn)參數(shù)及1個(gè)尺度比參數(shù),從而在2個(gè)橢球間建立聯(lián)系,實(shí)現(xiàn)由橢球高到正常高的轉(zhuǎn)換。此方法要求測(cè)區(qū)在已知控制點(diǎn)范圍內(nèi)或邊緣,不要超出控制區(qū)域太遠(yuǎn),否則得出的高程值不具可信性。同時(shí),該方法無法處理測(cè)區(qū)內(nèi)高程異常值凹下或凸起的區(qū)域。

    1.1.2高程擬合法

    在傳統(tǒng)的大地測(cè)量中,高程異常是通過重力測(cè)量和天文測(cè)量的方法確定的,但對(duì)大多數(shù)測(cè)量單位來說,并不具備這2種作業(yè)的條件。目前國(guó)內(nèi)外常用的方法是采用數(shù)學(xué)擬合的方法[2]。

    通過EGM2008大地水準(zhǔn)面模型,結(jié)合模型內(nèi)不同點(diǎn)的高程異常,以線性內(nèi)插的方式求出待測(cè)點(diǎn)高程異常值,從而實(shí)現(xiàn)橢球高到正常高的轉(zhuǎn)變。此方法在測(cè)區(qū)周圍控制點(diǎn)不多或不能有效覆蓋測(cè)區(qū)時(shí),可以獲取高可靠性的RTK高程。同時(shí),該方法降低了測(cè)區(qū)范圍內(nèi)局部地區(qū)高程異常值凹下或凸起帶來的危害。 由于海上作業(yè)的一側(cè)是海洋,很難搜集到覆蓋測(cè)區(qū)的控制點(diǎn)資料,在采用RTK測(cè)量海深的作業(yè)模式上,不推薦用七參數(shù)法。該次文采用高程擬合法。圖2為該次作業(yè)GS+導(dǎo)出的某條測(cè)線的部分導(dǎo)航數(shù)據(jù)。

    圖2 某測(cè)線的部分導(dǎo)航數(shù)據(jù)

    圖2第一行中的“4”表示固定解,“17.777”表示正常高,“-10.323”表示大地水準(zhǔn)面差距,兩者之和為橢球高。此處的大地水準(zhǔn)面差距是由集成在GPS接收機(jī)中粗略的大地水準(zhǔn)面模型得出,誤差較大,同時(shí)GS+的算法較簡(jiǎn)單,在每條測(cè)線只取了1個(gè)大地水準(zhǔn)面差距值。因此需要對(duì)導(dǎo)航數(shù)據(jù)進(jìn)行后處理,以便求出正確的正常高。

    基本思路為:首先利用公式H=h+ζ(式中H代表該點(diǎn)橢球高;h代表該點(diǎn)正常高;ζ代表該點(diǎn)高程異常,下同),恢復(fù)橢球高H,然后根據(jù)EGM2008模型查找測(cè)量點(diǎn)位區(qū)域內(nèi)四周點(diǎn)的ζ值,利用線性內(nèi)插的方式求出待測(cè)點(diǎn)位的ζ值,最后利用公式h=H-ζ求出待測(cè)點(diǎn)位的正常高。

    理論上,要求得測(cè)量點(diǎn)位的正常高,有如下兩種方法可用:

    方法一:直接法。要想求得點(diǎn)2的正常高,利用該點(diǎn)的橢球高直接減去該點(diǎn)的高程異常值,公式為:

    (1)

    方法二:間接法。若要求得點(diǎn)2的正常高,利用已知點(diǎn)1的正常高加上兩者間橢球高的差再減去兩者高程異常的差,公式為:

    (2)

    圖3 基準(zhǔn)站點(diǎn)、流動(dòng)站點(diǎn)橢球高與正常高關(guān)系圖

    如圖3所示,公式(1)(2)中h2為RTK測(cè)量點(diǎn)位正常高;h1即為已知點(diǎn)(基準(zhǔn)站點(diǎn))正常高;H1為基站點(diǎn)橢球高;H2為導(dǎo)航數(shù)據(jù)中正常高與大地水準(zhǔn)面差距的和;ζ1,ζ2由選用的EGM2008模型求得(公式中的正常高、橢球高、在實(shí)際工作中要考慮儀器高、天線高等)。

    方法一得出的正常高為EGM2008高程基準(zhǔn),要想得出1985國(guó)家高程基準(zhǔn)的正常高,還需減去兩者間的差距。另外,該方法需準(zhǔn)確搜集到已知點(diǎn)的高精度WGS-84坐標(biāo)(B,L,H)。方法二消除了WGS84坐標(biāo)H絕對(duì)精度低及2個(gè)高程基準(zhǔn)面差距的影響,實(shí)際工作中主要應(yīng)用此方法。

    該項(xiàng)目15km測(cè)區(qū)內(nèi),Δζ(ζ2—ζ1)從岸邊(基準(zhǔn)站)往測(cè)區(qū)深處基本呈線性逐步增大,在5km處值約為16cm,在10km處約為34cm,在15km處約為46cm。Δζ通常情況下表現(xiàn)為使測(cè)區(qū)海底傾斜,離基站越遠(yuǎn),傾斜(誤差)程度越深。因?yàn)樗且粋€(gè)系統(tǒng)誤差,通過布設(shè)重合、交叉測(cè)線的方式無法將其檢測(cè)出來。所以在用RTK結(jié)合多波束進(jìn)行水深測(cè)量時(shí),一定要考慮Δζ的影響。

    1.2RTK高程計(jì)算潮汐數(shù)據(jù)

    實(shí)際海上測(cè)量中,受風(fēng)浪、船體操縱等影響,GPS信號(hào)會(huì)出現(xiàn)失鎖現(xiàn)象,此時(shí)的GPS高程不能直接利用。目前,國(guó)內(nèi)主要研究采用對(duì)GPS高程數(shù)據(jù)進(jìn)行修正的方式來確保GPS高程數(shù)據(jù)的質(zhì)量,修正后的GPS高程數(shù)據(jù)直接用于求算換能器瞬時(shí)高程。其中常見的有Kalman濾波、Heave修正等。

    項(xiàng)目中,首先利用GPS高程求算隨船移動(dòng)的潮位面數(shù)據(jù),然后根據(jù)潮位面數(shù)據(jù)去求算換能器瞬時(shí)高程。

    1.2.1GPS高程常用修正法

    Kalman濾波法:Kalman濾波是一種對(duì)動(dòng)態(tài)測(cè)量數(shù)據(jù)實(shí)施濾波的常用方法,隨著觀測(cè)數(shù)據(jù)的增加,狀態(tài)估計(jì)應(yīng)愈來愈精確。但在實(shí)際應(yīng)用中,當(dāng)GPS觀測(cè)存在連續(xù)異常時(shí),則采用傳統(tǒng)的Kalman濾波易出現(xiàn)發(fā)散問題[3]。

    Heave修正法:Heave修正法利用的是正常的GPS高程時(shí)序和Heave時(shí)序反應(yīng)的都是船體的垂直運(yùn)動(dòng),在短時(shí)間內(nèi)應(yīng)該具有一致性的原理通過姿態(tài)傳感器測(cè)得的Heave值來對(duì)GPS高程進(jìn)行修正。但是,Heave具有零均值特征,其變化反映的是船體的高頻運(yùn)動(dòng)特征,無法表現(xiàn)船體的全周期垂直運(yùn)動(dòng),即潮位、波浪和船體操縱等產(chǎn)生的綜合垂直運(yùn)動(dòng);GPS高程序列不但可以呈現(xiàn)船體的瞬時(shí)變化,還可以呈現(xiàn)其中、長(zhǎng)周期垂直運(yùn)動(dòng)波。因此,Heave對(duì)GPS RTK高程的修正只能在短時(shí)間段內(nèi)進(jìn)行。根據(jù)潮位、波浪和船體操縱的周期特征,這種修正的有效時(shí)間段最長(zhǎng)不能大于2min[4]。

    經(jīng)過修正后的GPS高程數(shù)據(jù)直接用于求算換能器瞬時(shí)高程。得出的高程起著提供瞬時(shí)垂直基準(zhǔn)的重要作用,每個(gè)歷元的觀測(cè)數(shù)據(jù)均須可靠穩(wěn)定,否則將會(huì)引起水下地形的異常[4]。因此,該方法要求GPS原始數(shù)據(jù)信號(hào)失鎖時(shí)間不能太長(zhǎng),同時(shí)處理過程中對(duì)GPS高程成果的質(zhì)量要嚴(yán)格把控。

    1.2.2利用RTK高程求算潮位數(shù)據(jù)文件

    鑒于該項(xiàng)目多條測(cè)線出現(xiàn)了持續(xù)2~3min的RTK浮點(diǎn)解或單點(diǎn)解數(shù)據(jù),利用Kalman濾波、Heave修正無法消除GPS失鎖對(duì)水深成果造成的影響,提出了另一個(gè)思路。利用GPS中的固定解數(shù)據(jù)反算海平面,得到一個(gè)“架設(shè)在船上的驗(yàn)潮站”潮位文件,并利用該文件結(jié)合涌浪參數(shù)及吃水求算換能器瞬時(shí)高度。

    求算潮位文件的過程為: 首先,劃分一個(gè)合理的時(shí)間段。該時(shí)間段要包含多個(gè)涌浪周期(Heave-Bandwidth)及足夠多的高程數(shù)據(jù),同時(shí)考慮2次失鎖情況,該周期不能過大。該項(xiàng)目GPS輸出頻率為20Hz,涌浪周期為7S左右,每個(gè)周期內(nèi)GPS輸出140次,1 min內(nèi)輸出1200次,該時(shí)間段(區(qū)間)設(shè)置為1min。

    其次,設(shè)置區(qū)間內(nèi)有效數(shù)據(jù)閾值。設(shè)置區(qū)間內(nèi)GPS固定解的數(shù)量的閾值為x,如果固定解數(shù)量小于x,則將該區(qū)間內(nèi)所有數(shù)據(jù)舍棄;否則保留該區(qū)間數(shù)據(jù)。該項(xiàng)目處理過程中在將區(qū)間段設(shè)置為1 min的情況下,該閾值設(shè)置為1000。失鎖嚴(yán)重的測(cè)線,也可將該值適當(dāng)縮小,并對(duì)該時(shí)間段內(nèi)潮位文件進(jìn)行人工檢核。

    最后,計(jì)算該區(qū)間時(shí)刻海水潮位面。經(jīng)過第二步的過濾后,將每個(gè)區(qū)間內(nèi)所有固定解高程取平均,結(jié)合天線高及換能器吃水值,從而得出該區(qū)間時(shí)刻海水潮位面數(shù)據(jù)。

    獲取海水潮位面的目的是為了求算換能器高程,所以在第三步中,也可結(jié)合GPS天線與換能器的高差(圖1中的h1),直接求算“換能器潮位面”數(shù)據(jù)。

    2 結(jié)論

    與傳統(tǒng)驗(yàn)潮方法相比,利用RTK進(jìn)行多波束水深測(cè)量在高程獲取方面是完全不同的一種方法。該文通過茂名某海洋測(cè)量項(xiàng)目為例,針對(duì)高程轉(zhuǎn)換及RTK信號(hào)失鎖情況,在GPS高程數(shù)據(jù)處理方面,提出了以下幾個(gè)建議:

    (1)GGA原始數(shù)據(jù)中的大地水準(zhǔn)面差距值是由內(nèi)置在GPS接收機(jī)中粗略的大地水準(zhǔn)面模型給出,誤差較大,要想得出高精度水深數(shù)據(jù),首先要恢復(fù)點(diǎn)位的橢球高。

    (2)在GPS橢球高向正常高轉(zhuǎn)變過程中,優(yōu)先使用EGM2008模型進(jìn)行高程擬合,測(cè)區(qū)內(nèi)不同點(diǎn)位的高程異常差值Δζ(ζ2-ζ1)在10km范圍內(nèi),可達(dá)30cm以上。忽略Δζ會(huì)對(duì)測(cè)深結(jié)果產(chǎn)生非常大的系統(tǒng)差值,且簡(jiǎn)單的布設(shè)重復(fù)測(cè)線無法檢測(cè)出此項(xiàng)誤差。

    (3)在測(cè)量過程中,GPS存在短時(shí)間失鎖情況時(shí),可以通過利用GPS高程數(shù)據(jù)求算潮位數(shù)據(jù)的方式,消弱失鎖對(duì)成果的影響。該方法利用的是GPS高程數(shù)據(jù)里面的有效高程值,且可對(duì)生成的潮位數(shù)據(jù)圖像化顯示,利用海水周期性漲落的特點(diǎn)進(jìn)行人工檢核;而對(duì)GPS高程數(shù)據(jù)進(jìn)行Kalman濾波、Heave修正等方法,需要修正后的每個(gè)歷元的觀測(cè)數(shù)據(jù)均必須可靠、穩(wěn)定。所以利用GPS高程數(shù)據(jù)求算潮位數(shù)據(jù)的方法,對(duì)GPS數(shù)據(jù)質(zhì)量要求更寬松,求算出的水深結(jié)果更可靠。

    綜上所述,利用提出的對(duì)RTK-GPS高程數(shù)據(jù)的處理,可以求算隨船移動(dòng)的潮位面數(shù)據(jù),降低了短時(shí)間內(nèi)GPS信號(hào)失鎖對(duì)多波束水深測(cè)量數(shù)據(jù)的影響。

    [1]陽凡林,趙建虎,張紅梅,等.RTK高程和Heave信號(hào)的融合及精度分析[J].武漢大學(xué)學(xué)報(bào)信息科學(xué)版,2007,32(3):225-228.

    [2]程海帆,冷曦晨,郭寶君. GPS高程兩種擬合方法在山區(qū)公路勘測(cè)中的應(yīng)用研究[J].公路交通科技,2008,25(9):336.

    [3]暴景陽,劉雁春.海道水位控制方法研究[J].測(cè)繪科學(xué),2006,31(6):49-51.

    [4]張紅梅,趙建虎. 精密多波束測(cè)量中GPS高程誤差的綜合修正法[J]. 測(cè)繪學(xué)報(bào),2009,38(1):24-26.

    Elevation Data Processing of GPS in Multibeam Bathymetric Measuring by Using RTK Technology

    LIU Fugang

    (Sinopec Petroleum Engineering Corporation, Shandong Dongying 257000, China)

    In the period of measuring multi-beam by using RTK technology, GPS elevation plays an important role in providing instantaneous vertical datum. During the actual operation process in the sea, GPS signals often happen the problem of out of lock. It will cause mutations of some epoch elevation data, and distortion of submarine topography. In this paper, setting Maoming project as an example, according to EGM2008 geoid model, by using elevation fitting method, in the process of measuring multibeam bathymetric, GPS ellipsoidal height will convert to normal height. Then pointing to the phenomenon of losing and locking GPS signals, the method of using elevation fixed solution to calculate sea tide surface by using RTK technology has been put forward. The advantages of this method have been explained as well.

    RTK;multibeam; ellipsoid height; normal height; elevation data processing

    2015-06-24;

    2015-07-31;編輯:陶衛(wèi)衛(wèi)

    劉福剛(1983—),男,山東昌樂人,工程師,現(xiàn)從事海洋測(cè)量方面的研究;E-mail:slecclfg@163.com

    P228.42

    B

    猜你喜歡
    潮位橢球換能器
    基于距離倒數(shù)加權(quán)的多站潮位改正方法可行性分析
    獨(dú)立坐標(biāo)系橢球變換與坐標(biāo)換算
    橢球槽宏程序編制及其Vericut仿真
    智能制造(2021年4期)2021-11-04 08:54:44
    唐山市警戒潮位標(biāo)志物維護(hù)研究
    多潮位站海道地形測(cè)量潮位控制方法研究
    橢球精加工軌跡及程序設(shè)計(jì)
    基于外定界橢球集員估計(jì)的純方位目標(biāo)跟蹤
    基于改進(jìn)的OLS-RBF模型的感潮河段潮位預(yù)測(cè)研究
    一種推挽式驅(qū)動(dòng)Ⅶ型彎張換能器的探討
    IV型彎張換能器輻射阻抗的等效球模型計(jì)算
    五月伊人婷婷丁香| 亚洲欧美精品自产自拍| 久久99蜜桃精品久久| 国产欧美另类精品又又久久亚洲欧美| 嫩草影院入口| 精品99又大又爽又粗少妇毛片| 亚洲国产精品999| 日韩欧美 国产精品| 亚洲怡红院男人天堂| 精品少妇久久久久久888优播| 午夜爱爱视频在线播放| 免费黄频网站在线观看国产| 久久久久久久久久久免费av| 草草在线视频免费看| 精品少妇久久久久久888优播| 麻豆乱淫一区二区| 在线播放无遮挡| 国产女主播在线喷水免费视频网站| 99热网站在线观看| 春色校园在线视频观看| 成人美女网站在线观看视频| 亚洲国产高清在线一区二区三| 精品国产三级普通话版| 国产免费一区二区三区四区乱码| 高清在线视频一区二区三区| 日韩国内少妇激情av| 蜜臀久久99精品久久宅男| 一边亲一边摸免费视频| 国产美女午夜福利| 亚洲精品一二三| 日韩一本色道免费dvd| 久久人人爽人人爽人人片va| 在线免费观看不下载黄p国产| 国产免费又黄又爽又色| 男的添女的下面高潮视频| 大码成人一级视频| 熟女人妻精品中文字幕| 欧美老熟妇乱子伦牲交| 一区二区三区精品91| 大片电影免费在线观看免费| 五月天丁香电影| 精品亚洲乱码少妇综合久久| 欧美最新免费一区二区三区| 亚洲怡红院男人天堂| 91在线精品国自产拍蜜月| 国产精品99久久99久久久不卡 | 欧美成人精品欧美一级黄| 国产伦精品一区二区三区四那| 美女脱内裤让男人舔精品视频| 精品久久久久久久久亚洲| 亚洲最大成人av| 久久久久精品性色| 成人综合一区亚洲| 亚洲精品aⅴ在线观看| 七月丁香在线播放| 国产免费一区二区三区四区乱码| 国产日韩欧美亚洲二区| 欧美三级亚洲精品| 精品午夜福利在线看| 一个人观看的视频www高清免费观看| 亚洲精品456在线播放app| 另类亚洲欧美激情| 精品人妻偷拍中文字幕| 91久久精品国产一区二区三区| 日本欧美国产在线视频| 在线观看一区二区三区激情| 亚洲成人av在线免费| 黄片wwwwww| 一级二级三级毛片免费看| 丝袜脚勾引网站| 男女那种视频在线观看| 三级国产精品欧美在线观看| 国产成人精品婷婷| 免费黄色在线免费观看| 青春草视频在线免费观看| 九九在线视频观看精品| 国产欧美另类精品又又久久亚洲欧美| 麻豆成人午夜福利视频| 亚洲三级黄色毛片| 国内少妇人妻偷人精品xxx网站| 亚洲精品色激情综合| 亚洲精品国产av成人精品| 亚洲不卡免费看| 麻豆成人av视频| 亚洲av中文av极速乱| 亚洲精品乱久久久久久| 久久国内精品自在自线图片| 看黄色毛片网站| 午夜福利在线在线| 久久精品久久精品一区二区三区| av在线app专区| 最近中文字幕2019免费版| 国产精品一及| 国产成人aa在线观看| 看黄色毛片网站| 街头女战士在线观看网站| 综合色丁香网| tube8黄色片| 永久网站在线| 亚洲最大成人av| 一本一本综合久久| 熟妇人妻不卡中文字幕| 成年av动漫网址| 一区二区三区免费毛片| 亚洲国产精品成人久久小说| 国产精品蜜桃在线观看| 在线观看一区二区三区激情| 亚洲,欧美,日韩| 亚洲国产高清在线一区二区三| 91aial.com中文字幕在线观看| 免费不卡的大黄色大毛片视频在线观看| 黄色怎么调成土黄色| 国产精品嫩草影院av在线观看| 亚洲精品456在线播放app| 黄色欧美视频在线观看| 五月开心婷婷网| 汤姆久久久久久久影院中文字幕| 最近最新中文字幕免费大全7| 久久97久久精品| 蜜桃久久精品国产亚洲av| 国产高清不卡午夜福利| 国产精品三级大全| 蜜桃亚洲精品一区二区三区| 国产精品久久久久久久电影| 少妇人妻 视频| 99久久人妻综合| 夜夜看夜夜爽夜夜摸| 国产男人的电影天堂91| 少妇 在线观看| 亚洲国产精品999| 日韩中字成人| 欧美日韩国产mv在线观看视频 | 久久久精品94久久精品| 又爽又黄无遮挡网站| 精品酒店卫生间| 国产精品一及| 韩国高清视频一区二区三区| 国产黄片美女视频| 国产成人福利小说| 蜜桃久久精品国产亚洲av| 日本-黄色视频高清免费观看| 精品久久久久久久久av| 中文字幕制服av| 一级毛片黄色毛片免费观看视频| 欧美另类一区| 欧美成人午夜免费资源| 22中文网久久字幕| 日日摸夜夜添夜夜爱| 精品人妻一区二区三区麻豆| 男人舔奶头视频| 午夜亚洲福利在线播放| 99热这里只有是精品50| 国产黄色免费在线视频| videossex国产| 精品视频人人做人人爽| 五月伊人婷婷丁香| 国产人妻一区二区三区在| 日日摸夜夜添夜夜添av毛片| 高清日韩中文字幕在线| 成年免费大片在线观看| 日本免费在线观看一区| 少妇裸体淫交视频免费看高清| 欧美成人午夜免费资源| 男的添女的下面高潮视频| 日本av手机在线免费观看| 极品少妇高潮喷水抽搐| 欧美亚洲 丝袜 人妻 在线| 男男h啪啪无遮挡| 国产成人a区在线观看| 国产亚洲5aaaaa淫片| 91久久精品国产一区二区成人| 又黄又爽又刺激的免费视频.| 国产乱人视频| 久久午夜福利片| 午夜免费观看性视频| 深夜a级毛片| 一区二区三区免费毛片| 国产精品不卡视频一区二区| 国产欧美另类精品又又久久亚洲欧美| 男人狂女人下面高潮的视频| av在线蜜桃| 18禁在线播放成人免费| 中文字幕久久专区| 国产伦精品一区二区三区视频9| 亚洲综合精品二区| 搞女人的毛片| 久久精品综合一区二区三区| 精品人妻一区二区三区麻豆| 亚洲成人中文字幕在线播放| 简卡轻食公司| 夫妻性生交免费视频一级片| 欧美3d第一页| 国产精品一区二区性色av| 一级二级三级毛片免费看| 日本三级黄在线观看| 国产爱豆传媒在线观看| 久久久久九九精品影院| 国产成人精品福利久久| 嫩草影院精品99| 亚洲不卡免费看| 日韩av不卡免费在线播放| 国产成人a∨麻豆精品| 亚洲丝袜综合中文字幕| www.色视频.com| 国产成人aa在线观看| 白带黄色成豆腐渣| 欧美日韩视频高清一区二区三区二| 久久久久久久国产电影| 亚洲精品自拍成人| 王馨瑶露胸无遮挡在线观看| 不卡视频在线观看欧美| 九色成人免费人妻av| 极品教师在线视频| 日韩av在线免费看完整版不卡| 欧美激情在线99| 直男gayav资源| 久久久午夜欧美精品| 欧美成人a在线观看| 看黄色毛片网站| 亚洲av.av天堂| 99热全是精品| 日本黄大片高清| 亚洲,一卡二卡三卡| 日韩视频在线欧美| 国产精品久久久久久精品电影| 97精品久久久久久久久久精品| h日本视频在线播放| 国产亚洲最大av| 国产在线一区二区三区精| 亚洲精品久久久久久婷婷小说| 精品久久久久久久末码| 欧美人与善性xxx| 亚州av有码| 又爽又黄无遮挡网站| 三级国产精品欧美在线观看| 性色avwww在线观看| 欧美激情国产日韩精品一区| 欧美精品国产亚洲| 亚洲国产精品国产精品| 久久99蜜桃精品久久| 国产精品熟女久久久久浪| 在线精品无人区一区二区三 | 亚洲婷婷狠狠爱综合网| 免费看av在线观看网站| 美女被艹到高潮喷水动态| 中文欧美无线码| av一本久久久久| 亚洲aⅴ乱码一区二区在线播放| 久热久热在线精品观看| 国产成人freesex在线| 国内揄拍国产精品人妻在线| 男男h啪啪无遮挡| 99久久人妻综合| 免费观看在线日韩| 九色成人免费人妻av| 51国产日韩欧美| 日韩在线高清观看一区二区三区| 国产亚洲5aaaaa淫片| 中国国产av一级| 欧美丝袜亚洲另类| 人妻一区二区av| 亚洲精品成人久久久久久| 一二三四中文在线观看免费高清| 中文资源天堂在线| 夫妻午夜视频| 久久精品国产鲁丝片午夜精品| 国产精品福利在线免费观看| 色婷婷久久久亚洲欧美| 国产精品一区二区三区四区免费观看| 欧美一区二区亚洲| 小蜜桃在线观看免费完整版高清| 中文精品一卡2卡3卡4更新| 亚洲精品国产av蜜桃| 亚洲精品亚洲一区二区| 在线 av 中文字幕| 国产成人精品久久久久久| 日日摸夜夜添夜夜爱| 欧美日本视频| 综合色av麻豆| 国产午夜精品一二区理论片| 一级片'在线观看视频| 亚洲美女搞黄在线观看| 国产爱豆传媒在线观看| 国精品久久久久久国模美| 99久久中文字幕三级久久日本| 中国国产av一级| 国产一级毛片在线| 在线a可以看的网站| 我的女老师完整版在线观看| 国产亚洲精品久久久com| 免费少妇av软件| 亚洲精品日韩在线中文字幕| 午夜亚洲福利在线播放| 丝袜喷水一区| 亚洲精品久久久久久婷婷小说| 亚洲va在线va天堂va国产| 少妇的逼水好多| 搡女人真爽免费视频火全软件| 丝袜脚勾引网站| 国产白丝娇喘喷水9色精品| 久久精品国产自在天天线| 午夜福利视频1000在线观看| 人妻 亚洲 视频| 中文欧美无线码| 午夜福利视频精品| 尤物成人国产欧美一区二区三区| 人妻夜夜爽99麻豆av| 麻豆国产97在线/欧美| 人妻系列 视频| 国产成人精品福利久久| 日韩三级伦理在线观看| 最近2019中文字幕mv第一页| 亚洲欧美日韩东京热| 日日撸夜夜添| 国产片特级美女逼逼视频| 人人妻人人看人人澡| 街头女战士在线观看网站| av卡一久久| 国产成人精品一,二区| 蜜臀久久99精品久久宅男| 大片免费播放器 马上看| 在线观看一区二区三区激情| 亚洲精品久久久久久婷婷小说| 亚洲国产精品成人综合色| 中文在线观看免费www的网站| 深夜a级毛片| 日韩一区二区三区影片| 午夜福利在线在线| 色视频www国产| 日韩av不卡免费在线播放| 亚洲国产色片| 中文字幕av成人在线电影| 国产片特级美女逼逼视频| 亚洲欧美精品自产自拍| 成年人午夜在线观看视频| 亚洲成色77777| 国产亚洲午夜精品一区二区久久 | 国语对白做爰xxxⅹ性视频网站| 国产高清三级在线| 国产免费一级a男人的天堂| 噜噜噜噜噜久久久久久91| 亚洲人成网站在线播| 精品熟女少妇av免费看| 久久ye,这里只有精品| 亚洲久久久久久中文字幕| 大话2 男鬼变身卡| 久久影院123| 伊人久久国产一区二区| 男女边摸边吃奶| 国产精品久久久久久av不卡| 国产探花在线观看一区二区| 久久99蜜桃精品久久| 一级毛片我不卡| 97超碰精品成人国产| 一级片'在线观看视频| 夜夜爽夜夜爽视频| 一级片'在线观看视频| 亚洲三级黄色毛片| 一级毛片电影观看| 亚洲精品一二三| 久久影院123| 日日摸夜夜添夜夜添av毛片| 日本爱情动作片www.在线观看| 内射极品少妇av片p| 成人漫画全彩无遮挡| 国产伦理片在线播放av一区| 偷拍熟女少妇极品色| 人人妻人人爽人人添夜夜欢视频 | 日韩欧美精品v在线| 啦啦啦啦在线视频资源| 3wmmmm亚洲av在线观看| 久久精品国产a三级三级三级| 国产伦精品一区二区三区视频9| 亚洲精品乱码久久久v下载方式| 视频区图区小说| 内地一区二区视频在线| 亚洲伊人久久精品综合| 老师上课跳d突然被开到最大视频| 色吧在线观看| 国产色婷婷99| 亚洲精品国产成人久久av| 国产在线男女| 亚洲av电影在线观看一区二区三区 | 国产精品一区二区性色av| 在线 av 中文字幕| 亚洲av国产av综合av卡| 美女高潮的动态| 日本色播在线视频| 直男gayav资源| 我的女老师完整版在线观看| 成人综合一区亚洲| 亚洲欧美成人综合另类久久久| 秋霞在线观看毛片| 欧美日韩一区二区视频在线观看视频在线 | 欧美性感艳星| av在线app专区| 亚洲精品成人久久久久久| 亚洲精品成人av观看孕妇| 国产熟女欧美一区二区| av在线天堂中文字幕| av在线老鸭窝| 草草在线视频免费看| 亚洲精品第二区| 插逼视频在线观看| 国内少妇人妻偷人精品xxx网站| 久久人人爽av亚洲精品天堂 | 久久这里有精品视频免费| 国产精品久久久久久精品电影| 久久精品久久久久久久性| 亚洲av日韩在线播放| 日日摸夜夜添夜夜爱| 91精品一卡2卡3卡4卡| 草草在线视频免费看| 热re99久久精品国产66热6| 国内精品美女久久久久久| 国精品久久久久久国模美| 国产又色又爽无遮挡免| 国产成年人精品一区二区| 成人亚洲精品av一区二区| 国产永久视频网站| 亚洲av免费高清在线观看| 日韩,欧美,国产一区二区三区| 亚洲国产欧美人成| 99久久人妻综合| 欧美日韩视频精品一区| 亚洲国产精品成人久久小说| 欧美日韩综合久久久久久| 91在线精品国自产拍蜜月| 熟女av电影| 久久久国产一区二区| 又爽又黄a免费视频| 精华霜和精华液先用哪个| 欧美xxxx性猛交bbbb| 九九在线视频观看精品| 久久6这里有精品| 国产又色又爽无遮挡免| 九色成人免费人妻av| 亚洲av在线观看美女高潮| 91aial.com中文字幕在线观看| 乱系列少妇在线播放| 亚洲国产高清在线一区二区三| 69av精品久久久久久| 亚洲,欧美,日韩| 久久精品综合一区二区三区| 国产日韩欧美在线精品| 国产精品秋霞免费鲁丝片| 欧美日韩亚洲高清精品| 亚洲av国产av综合av卡| 色婷婷久久久亚洲欧美| 插逼视频在线观看| 亚洲婷婷狠狠爱综合网| 高清视频免费观看一区二区| 日韩强制内射视频| 国产色爽女视频免费观看| 九九久久精品国产亚洲av麻豆| 国产中年淑女户外野战色| 男的添女的下面高潮视频| 国语对白做爰xxxⅹ性视频网站| 免费大片黄手机在线观看| 激情五月婷婷亚洲| 国产成人a区在线观看| 久久久久久伊人网av| 亚洲国产成人一精品久久久| 乱码一卡2卡4卡精品| 2021天堂中文幕一二区在线观| 成年人午夜在线观看视频| 全区人妻精品视频| 王馨瑶露胸无遮挡在线观看| 黄色一级大片看看| 中国国产av一级| 看黄色毛片网站| 色哟哟·www| 亚洲在线观看片| 一本色道久久久久久精品综合| 成年av动漫网址| eeuss影院久久| 免费看日本二区| 亚洲欧美精品专区久久| 韩国高清视频一区二区三区| 2022亚洲国产成人精品| 国产中年淑女户外野战色| 22中文网久久字幕| 亚洲国产精品成人久久小说| 寂寞人妻少妇视频99o| 日本一二三区视频观看| 欧美最新免费一区二区三区| 色播亚洲综合网| 免费看光身美女| 少妇裸体淫交视频免费看高清| 国产伦理片在线播放av一区| 久久ye,这里只有精品| 亚洲三级黄色毛片| 自拍欧美九色日韩亚洲蝌蚪91 | 久久综合国产亚洲精品| 色5月婷婷丁香| a级毛片免费高清观看在线播放| 午夜免费观看性视频| 国产精品三级大全| 精品久久久久久久人妻蜜臀av| 国产黄a三级三级三级人| 亚洲三级黄色毛片| 免费不卡的大黄色大毛片视频在线观看| 亚洲天堂国产精品一区在线| 国产伦精品一区二区三区视频9| 美女视频免费永久观看网站| 欧美日韩精品成人综合77777| 亚洲国产成人一精品久久久| 免费观看a级毛片全部| 欧美成人精品欧美一级黄| 国产真实伦视频高清在线观看| 久久精品久久久久久久性| 深夜a级毛片| 国产精品蜜桃在线观看| 免费观看a级毛片全部| 国产精品久久久久久久久免| 国产 一区精品| 男人狂女人下面高潮的视频| 深夜a级毛片| 亚洲精品国产av蜜桃| 国产成人91sexporn| 日韩 亚洲 欧美在线| 国产成人精品久久久久久| 国产精品嫩草影院av在线观看| 夜夜看夜夜爽夜夜摸| 一本久久精品| 欧美日韩一区二区视频在线观看视频在线 | 舔av片在线| 麻豆成人av视频| 午夜亚洲福利在线播放| 国产成人精品福利久久| 一区二区三区乱码不卡18| 国产午夜福利久久久久久| 国产精品久久久久久精品电影小说 | 免费av观看视频| 免费大片18禁| 国产一区二区三区综合在线观看 | 亚洲国产欧美在线一区| 精品视频人人做人人爽| av一本久久久久| 国产伦在线观看视频一区| 久久99热这里只频精品6学生| 成人免费观看视频高清| 丰满少妇做爰视频| 你懂的网址亚洲精品在线观看| 简卡轻食公司| 亚洲精品,欧美精品| 亚洲精品国产av蜜桃| 国产成人91sexporn| 久久久久久久亚洲中文字幕| 五月玫瑰六月丁香| 国产欧美日韩精品一区二区| 男女国产视频网站| 国产探花极品一区二区| 欧美zozozo另类| 亚洲经典国产精华液单| 国产在视频线精品| 亚洲精品aⅴ在线观看| 亚洲精品第二区| 亚洲熟女精品中文字幕| 精品99又大又爽又粗少妇毛片| 精品久久久久久电影网| 国产一区有黄有色的免费视频| 国产成人福利小说| 亚洲精华国产精华液的使用体验| 免费大片18禁| 精品酒店卫生间| 国产美女午夜福利| 久久99蜜桃精品久久| 18禁动态无遮挡网站| 亚洲伊人久久精品综合| 亚洲在线观看片| 久久久精品94久久精品| 国产成人一区二区在线| 亚洲国产精品成人综合色| 人体艺术视频欧美日本| 一区二区三区乱码不卡18| 人妻夜夜爽99麻豆av| 国产中年淑女户外野战色| 国产 一区精品| 亚洲精品一区蜜桃| 国产亚洲91精品色在线| 日韩在线高清观看一区二区三区| 久久久久久久精品精品| 老司机影院毛片| 男人爽女人下面视频在线观看| 婷婷色麻豆天堂久久| 欧美成人精品欧美一级黄| 亚洲综合色惰| 国产探花在线观看一区二区| 我要看日韩黄色一级片| 内射极品少妇av片p| 王馨瑶露胸无遮挡在线观看| 日本黄大片高清| 日韩一区二区三区影片| 少妇丰满av| 久久久久久久亚洲中文字幕| 国产亚洲最大av| 色综合色国产| 99久国产av精品国产电影| 亚洲精品日韩在线中文字幕| 国产v大片淫在线免费观看| 免费av毛片视频| 青春草视频在线免费观看| 99九九线精品视频在线观看视频| 大香蕉久久网| 日韩免费高清中文字幕av| 日韩av不卡免费在线播放| 国产v大片淫在线免费观看| 高清日韩中文字幕在线| 成人午夜精彩视频在线观看| 嫩草影院新地址| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品456在线播放app| 亚洲一级一片aⅴ在线观看| 小蜜桃在线观看免费完整版高清|