• <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ì)算
    午夜精品在线福利| 国产精品久久久久久久电影 | 他把我摸到了高潮在线观看| 欧美成人一区二区免费高清观看 | 国产免费av片在线观看野外av| 国产精品自产拍在线观看55亚洲| 亚洲色图 男人天堂 中文字幕| 精品午夜福利视频在线观看一区| 成人高潮视频无遮挡免费网站| 精品电影一区二区在线| 亚洲无线在线观看| 一进一出抽搐动态| 99久久精品一区二区三区| 亚洲人成网站在线播放欧美日韩| 久久精品aⅴ一区二区三区四区| 日本免费a在线| 国产淫片久久久久久久久 | 午夜福利在线观看吧| 久久久久国产一级毛片高清牌| 国产精品电影一区二区三区| 97碰自拍视频| 国产精品av久久久久免费| 国产麻豆成人av免费视频| 看免费av毛片| 久久精品人妻少妇| av在线天堂中文字幕| 三级毛片av免费| 两个人看的免费小视频| 色哟哟哟哟哟哟| 欧美黄色片欧美黄色片| 久久久国产成人免费| 精品久久久久久久久久免费视频| 中文亚洲av片在线观看爽| 草草在线视频免费看| 欧美成人性av电影在线观看| 日本在线视频免费播放| 小说图片视频综合网站| 亚洲av美国av| 日本一二三区视频观看| 久久热在线av| 国产精品99久久久久久久久| 中文字幕人成人乱码亚洲影| 夜夜爽天天搞| 美女扒开内裤让男人捅视频| 变态另类成人亚洲欧美熟女| 久久99热这里只有精品18| 老鸭窝网址在线观看| 亚洲avbb在线观看| 精品不卡国产一区二区三区| 亚洲熟妇中文字幕五十中出| 国产成人aa在线观看| 日日摸夜夜添夜夜添小说| 成人av一区二区三区在线看| 成年女人永久免费观看视频| 日本黄色片子视频| www日本在线高清视频| 国产成人系列免费观看| 夜夜爽天天搞| 国产精品久久电影中文字幕| 国产亚洲精品久久久久久毛片| 99re在线观看精品视频| 十八禁人妻一区二区| 久久午夜综合久久蜜桃| 一进一出抽搐gif免费好疼| 长腿黑丝高跟| 久9热在线精品视频| 高潮久久久久久久久久久不卡| 欧美日韩综合久久久久久 | 免费看日本二区| 欧美绝顶高潮抽搐喷水| www.999成人在线观看| 亚洲av第一区精品v没综合| 亚洲五月天丁香| 香蕉国产在线看| 两人在一起打扑克的视频| 国产精品 欧美亚洲| 脱女人内裤的视频| 国产伦在线观看视频一区| 老司机在亚洲福利影院| 亚洲 欧美 日韩 在线 免费| 嫩草影视91久久| 午夜免费观看网址| a级毛片a级免费在线| 亚洲黑人精品在线| 桃红色精品国产亚洲av| 高清毛片免费观看视频网站| 国产视频内射| 欧美国产日韩亚洲一区| 窝窝影院91人妻| 国产精品亚洲美女久久久| 成人特级黄色片久久久久久久| 午夜免费成人在线视频| 久久精品影院6| а√天堂www在线а√下载| 99久久无色码亚洲精品果冻| 亚洲成av人片免费观看| 成人三级做爰电影| 怎么达到女性高潮| 国产乱人伦免费视频| 国产视频一区二区在线看| 国产野战对白在线观看| 国产一区二区在线av高清观看| 国产欧美日韩一区二区精品| 免费在线观看日本一区| 亚洲电影在线观看av| 色吧在线观看| 日本 欧美在线| 亚洲第一欧美日韩一区二区三区| 小蜜桃在线观看免费完整版高清| 亚洲电影在线观看av| 97人妻精品一区二区三区麻豆| 亚洲av成人av| 亚洲熟妇熟女久久| 草草在线视频免费看| 五月玫瑰六月丁香| 亚洲欧美日韩卡通动漫| 毛片女人毛片| x7x7x7水蜜桃| 欧美另类亚洲清纯唯美| 免费一级毛片在线播放高清视频| 亚洲色图av天堂| 欧美日韩综合久久久久久 | 每晚都被弄得嗷嗷叫到高潮| 久久精品夜夜夜夜夜久久蜜豆| 听说在线观看完整版免费高清| 精品一区二区三区视频在线 | 欧美绝顶高潮抽搐喷水| 男女床上黄色一级片免费看| 性欧美人与动物交配| 国产精品一区二区三区四区免费观看 | 国产男靠女视频免费网站| 嫩草影院精品99| 久久午夜亚洲精品久久| 亚洲男人的天堂狠狠| 亚洲国产精品久久男人天堂| 亚洲一区二区三区不卡视频| 国产99白浆流出| 男人舔女人下体高潮全视频| 日本一二三区视频观看| 久久久久久久精品吃奶| 国产视频内射| 亚洲欧洲精品一区二区精品久久久| 亚洲av中文字字幕乱码综合| 老司机在亚洲福利影院| 午夜视频精品福利| 男人舔女人下体高潮全视频| 国产成人福利小说| 免费看十八禁软件| 国产精品亚洲av一区麻豆| 亚洲欧美日韩高清在线视频| 精品福利观看| 亚洲欧美日韩无卡精品| 久久久水蜜桃国产精品网| 91久久精品国产一区二区成人 | 亚洲aⅴ乱码一区二区在线播放| 国产aⅴ精品一区二区三区波| 国产一级毛片七仙女欲春2| 久久久久久大精品| 天堂影院成人在线观看| 女人高潮潮喷娇喘18禁视频| 美女cb高潮喷水在线观看 | 69av精品久久久久久| 久久国产精品人妻蜜桃| 99国产极品粉嫩在线观看| 成人午夜高清在线视频| 亚洲在线自拍视频| 美女免费视频网站| 欧美日韩一级在线毛片| 韩国av一区二区三区四区| 国产精品av视频在线免费观看| 国产亚洲av嫩草精品影院| 色综合婷婷激情| 少妇丰满av| 午夜精品在线福利| 欧美不卡视频在线免费观看| 免费高清视频大片| 欧美另类亚洲清纯唯美| 久久久久亚洲av毛片大全| 欧美在线一区亚洲| 99久久99久久久精品蜜桃| 午夜a级毛片| 国产淫片久久久久久久久 | 老司机福利观看| 亚洲中文av在线| 国产麻豆成人av免费视频| 观看免费一级毛片| 欧美乱码精品一区二区三区| 久久久久久九九精品二区国产| 中文字幕人妻丝袜一区二区| 免费搜索国产男女视频| 午夜精品久久久久久毛片777| 免费在线观看成人毛片| 亚洲人成网站高清观看| 婷婷精品国产亚洲av| 巨乳人妻的诱惑在线观看| 亚洲av日韩精品久久久久久密| 国产成人影院久久av| 亚洲avbb在线观看| 少妇熟女aⅴ在线视频| 丁香欧美五月| 国产黄a三级三级三级人| 国产亚洲av嫩草精品影院| 国产综合懂色| 欧美日韩精品网址| 神马国产精品三级电影在线观看| 亚洲av美国av| 黄片大片在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 无限看片的www在线观看| 成年女人看的毛片在线观看| 啪啪无遮挡十八禁网站| 黑人操中国人逼视频| 欧美午夜高清在线| 特大巨黑吊av在线直播| 18禁观看日本| 免费观看的影片在线观看| 欧美日韩国产亚洲二区| 国产成人福利小说| 国产成人啪精品午夜网站| 18禁观看日本| 久久性视频一级片| 网址你懂的国产日韩在线| 久久久久久国产a免费观看| 欧美成人免费av一区二区三区| 亚洲国产精品合色在线| 国产精品一及| 亚洲欧美一区二区三区黑人| 视频区欧美日本亚洲| 日本撒尿小便嘘嘘汇集6| 日韩国内少妇激情av| 美女免费视频网站| 99国产综合亚洲精品| av女优亚洲男人天堂 | 日韩欧美精品v在线| 露出奶头的视频| 麻豆国产av国片精品| 色尼玛亚洲综合影院| 欧美日韩精品网址| 欧美午夜高清在线| 中亚洲国语对白在线视频| 久久伊人香网站| 久久国产精品影院| 国产免费av片在线观看野外av| 免费看十八禁软件| xxx96com| 一区二区三区激情视频| 亚洲18禁久久av| 精品无人区乱码1区二区| 亚洲中文字幕日韩| 母亲3免费完整高清在线观看| 99国产精品99久久久久| 免费高清视频大片| 国产高清三级在线| 少妇人妻一区二区三区视频| 99热6这里只有精品| 久久精品91蜜桃| 欧美在线一区亚洲| 国产亚洲精品一区二区www| 十八禁网站免费在线| 午夜成年电影在线免费观看| www日本黄色视频网| 国产麻豆成人av免费视频| 亚洲人成网站在线播放欧美日韩| 亚洲五月天丁香| 国产伦在线观看视频一区| 长腿黑丝高跟| 国产成人av激情在线播放| 国内揄拍国产精品人妻在线| 久久草成人影院| 999久久久国产精品视频| 在线观看午夜福利视频| 亚洲欧美日韩东京热| 视频区欧美日本亚洲| 国产精品一区二区免费欧美| 婷婷六月久久综合丁香| 精品日产1卡2卡| 久久精品人妻少妇| 亚洲成人久久爱视频| 亚洲熟女毛片儿| а√天堂www在线а√下载| 色综合婷婷激情| 亚洲国产欧洲综合997久久,| 无人区码免费观看不卡| 国产 一区 欧美 日韩| 成年女人毛片免费观看观看9| 成人国产综合亚洲| 99在线人妻在线中文字幕| 特大巨黑吊av在线直播| 久久精品91蜜桃| 观看美女的网站| 中文字幕久久专区| 2021天堂中文幕一二区在线观| 久久久久久久精品吃奶| 国产高清videossex| 亚洲熟女毛片儿| 丰满的人妻完整版| 一边摸一边抽搐一进一小说| 欧洲精品卡2卡3卡4卡5卡区| 日日夜夜操网爽| 亚洲专区中文字幕在线| 久久精品国产清高在天天线| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 国产蜜桃级精品一区二区三区| 午夜福利视频1000在线观看| av黄色大香蕉| 国产精品综合久久久久久久免费| 国产三级黄色录像| 性欧美人与动物交配| 黄色成人免费大全| 变态另类丝袜制服| 久久久久性生活片| 天天一区二区日本电影三级| 久久婷婷人人爽人人干人人爱| 国产欧美日韩一区二区精品| 日日夜夜操网爽| 少妇丰满av| 欧美色欧美亚洲另类二区| 亚洲成人免费电影在线观看| 俄罗斯特黄特色一大片| 久久国产精品人妻蜜桃| 91麻豆精品激情在线观看国产| 一卡2卡三卡四卡精品乱码亚洲| 日韩成人在线观看一区二区三区| 国产成人欧美在线观看| 日韩成人在线观看一区二区三区| 亚洲精品456在线播放app | 久久这里只有精品19| 国产精品亚洲av一区麻豆| 黄色日韩在线| 变态另类丝袜制服| 老司机午夜十八禁免费视频| 在线国产一区二区在线| 村上凉子中文字幕在线| 色播亚洲综合网| 一本精品99久久精品77| 久久久久久九九精品二区国产| а√天堂www在线а√下载| 九色国产91popny在线| 久久精品人妻少妇| 国产精品日韩av在线免费观看| 成人三级做爰电影| 久久香蕉精品热| 综合色av麻豆| 美女 人体艺术 gogo| 天天一区二区日本电影三级| 成年版毛片免费区| 五月伊人婷婷丁香| 在线观看午夜福利视频| 小蜜桃在线观看免费完整版高清| 精品国产三级普通话版| 日本在线视频免费播放| 久久久成人免费电影| 亚洲国产中文字幕在线视频| 久久精品国产亚洲av香蕉五月| 亚洲男人的天堂狠狠| 色噜噜av男人的天堂激情| 精品午夜福利视频在线观看一区| 国产一区在线观看成人免费| 蜜桃久久精品国产亚洲av| 午夜a级毛片| 亚洲av电影不卡..在线观看| 国产主播在线观看一区二区| 亚洲天堂国产精品一区在线| 国产成年人精品一区二区| 男女那种视频在线观看| 欧美午夜高清在线| 亚洲在线自拍视频| 女生性感内裤真人,穿戴方法视频| 成人一区二区视频在线观看| 99在线人妻在线中文字幕| 精品国内亚洲2022精品成人| 亚洲美女视频黄频| 午夜激情福利司机影院| 成年版毛片免费区| 18禁国产床啪视频网站| 一进一出抽搐动态| 757午夜福利合集在线观看| 天堂动漫精品| 成人三级做爰电影| 老熟妇仑乱视频hdxx| 成人无遮挡网站| 成人精品一区二区免费| 久久精品aⅴ一区二区三区四区| 久久这里只有精品19| 99视频精品全部免费 在线 | 亚洲av中文字字幕乱码综合| 国产又色又爽无遮挡免费看| 国产日本99.免费观看| 激情在线观看视频在线高清| 舔av片在线| 两个人看的免费小视频| 色在线成人网| 亚洲国产色片| 超碰成人久久| 99在线人妻在线中文字幕| 亚洲男人的天堂狠狠| 男人舔女人的私密视频| 欧美成人一区二区免费高清观看 | 91av网站免费观看| 免费观看精品视频网站| 婷婷六月久久综合丁香| e午夜精品久久久久久久| 国产高潮美女av| 俺也久久电影网| 黄色成人免费大全| а√天堂www在线а√下载| 99国产精品99久久久久| 国产精品久久久久久精品电影| xxxwww97欧美| 熟女电影av网| 男女之事视频高清在线观看| 日日摸夜夜添夜夜添小说| 久久亚洲精品不卡| 精品欧美国产一区二区三| 又黄又爽又免费观看的视频| 欧美性猛交黑人性爽| 黑人欧美特级aaaaaa片| 国产高清激情床上av| 国内精品久久久久久久电影| 亚洲中文字幕日韩| 岛国在线免费视频观看| 中文字幕av在线有码专区| 久久天堂一区二区三区四区| 久久久国产成人精品二区| 久久久久久人人人人人| 在线观看午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久| 色综合站精品国产| xxx96com| 最新中文字幕久久久久 | 搞女人的毛片| a级毛片a级免费在线| 最新美女视频免费是黄的| 国产高清激情床上av| 亚洲第一欧美日韩一区二区三区| 国产免费av片在线观看野外av| 久久久久久久久免费视频了| 欧美一区二区国产精品久久精品| 日韩国内少妇激情av| 男女视频在线观看网站免费| 国内精品美女久久久久久| 亚洲色图 男人天堂 中文字幕| 大型黄色视频在线免费观看| 中出人妻视频一区二区| 欧美午夜高清在线| 天堂av国产一区二区熟女人妻| 国产精品av久久久久免费| 成人国产综合亚洲| 欧美乱妇无乱码| 成人国产一区最新在线观看| 婷婷丁香在线五月| 99热只有精品国产| 99久久99久久久精品蜜桃| 51午夜福利影视在线观看| 变态另类成人亚洲欧美熟女| 精品熟女少妇八av免费久了| 国产成人系列免费观看| 国产综合懂色| 露出奶头的视频| 神马国产精品三级电影在线观看| 亚洲人成伊人成综合网2020| 美女午夜性视频免费| 精品久久久久久久久久久久久| АⅤ资源中文在线天堂| 国产伦在线观看视频一区| 啪啪无遮挡十八禁网站| 两个人的视频大全免费| 国产97色在线日韩免费| 久久久成人免费电影| 啦啦啦韩国在线观看视频| 欧美性猛交╳xxx乱大交人| 亚洲av五月六月丁香网| 男人和女人高潮做爰伦理| 国产三级中文精品| 成年免费大片在线观看| 亚洲精品一区av在线观看| 热99re8久久精品国产| 国内精品美女久久久久久| 欧美黑人巨大hd| 欧美绝顶高潮抽搐喷水| 国产精品av久久久久免费| 亚洲av中文字字幕乱码综合| 国产精品av视频在线免费观看| 成年免费大片在线观看| 国产男靠女视频免费网站| 成人性生交大片免费视频hd| 日韩欧美 国产精品| 国产淫片久久久久久久久 | 亚洲精品中文字幕一二三四区| 真人做人爱边吃奶动态| 老汉色av国产亚洲站长工具| 特大巨黑吊av在线直播| 日本五十路高清| 亚洲欧美激情综合另类| 夜夜看夜夜爽夜夜摸| 在线观看午夜福利视频| 岛国在线观看网站| 欧美性猛交╳xxx乱大交人| 国产激情久久老熟女| ponron亚洲| 国产精品美女特级片免费视频播放器 | 精品电影一区二区在线| 99在线视频只有这里精品首页| 久久香蕉精品热| 亚洲欧美日韩卡通动漫| 天天一区二区日本电影三级| 欧美性猛交黑人性爽| 99久久精品一区二区三区| 亚洲,欧美精品.| 法律面前人人平等表现在哪些方面| 搡老岳熟女国产| 亚洲第一欧美日韩一区二区三区| 久久国产乱子伦精品免费另类| 99热6这里只有精品| 国产三级中文精品| 国产成人av教育| 久久精品91蜜桃| 美女cb高潮喷水在线观看 | 十八禁人妻一区二区| 日韩精品中文字幕看吧| 91久久精品国产一区二区成人 | 99精品在免费线老司机午夜| 国产av不卡久久| 少妇的逼水好多| 啪啪无遮挡十八禁网站| 男人舔女人的私密视频| 欧洲精品卡2卡3卡4卡5卡区| 色综合站精品国产| 免费av不卡在线播放| 午夜激情欧美在线| 免费av毛片视频| 午夜日韩欧美国产| 中国美女看黄片| 人人妻,人人澡人人爽秒播| 久久久久精品国产欧美久久久| 天堂√8在线中文| 小蜜桃在线观看免费完整版高清| a级毛片a级免费在线| 搞女人的毛片| 少妇裸体淫交视频免费看高清| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩乱码在线| 亚洲国产精品合色在线| 免费人成视频x8x8入口观看| 婷婷亚洲欧美| 国内毛片毛片毛片毛片毛片| 日韩欧美在线乱码| 偷拍熟女少妇极品色| 两性午夜刺激爽爽歪歪视频在线观看| 中文字幕熟女人妻在线| 搡老熟女国产l中国老女人| 国产精品女同一区二区软件 | 久久久久久九九精品二区国产| 熟女人妻精品中文字幕| 男人舔女人的私密视频| 非洲黑人性xxxx精品又粗又长| 男人舔女人的私密视频| 久久久久性生活片| 天天躁狠狠躁夜夜躁狠狠躁| 欧美中文日本在线观看视频| 搡老岳熟女国产| 日韩欧美在线二视频| 淫秽高清视频在线观看| 亚洲成人免费电影在线观看| 成年女人毛片免费观看观看9| 国产欧美日韩精品亚洲av| 手机成人av网站| 国产麻豆成人av免费视频| 他把我摸到了高潮在线观看| 亚洲一区二区三区不卡视频| 后天国语完整版免费观看| 国产不卡一卡二| 真人做人爱边吃奶动态| www.www免费av| 日韩欧美精品v在线| 免费观看精品视频网站| 嫩草影院入口| 精品人妻1区二区| 精品久久久久久久久久免费视频| 最近最新中文字幕大全免费视频| ponron亚洲| 久久久久精品国产欧美久久久| 久久中文字幕一级| 每晚都被弄得嗷嗷叫到高潮| 日本在线视频免费播放| 国内毛片毛片毛片毛片毛片| 两个人的视频大全免费| 成人18禁在线播放| 在线观看美女被高潮喷水网站 | 真人做人爱边吃奶动态| 又大又爽又粗| 亚洲成人中文字幕在线播放| 一级毛片精品| 一边摸一边抽搐一进一小说| 亚洲国产精品999在线| 国产精品野战在线观看| 丝袜人妻中文字幕| 在线视频色国产色| 国产激情偷乱视频一区二区| 一本久久中文字幕| 午夜福利欧美成人| 9191精品国产免费久久| 精品欧美国产一区二区三| 男女做爰动态图高潮gif福利片| 欧美日韩国产亚洲二区| 在线观看日韩欧美| 国产欧美日韩一区二区三| 免费看光身美女| 国产成人一区二区三区免费视频网站| 在线免费观看的www视频| 国产毛片a区久久久久| aaaaa片日本免费| 免费高清视频大片| 日本成人三级电影网站|