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

    一種用于重力梯度動(dòng)態(tài)測(cè)量的載體環(huán)境引力梯度補(bǔ)償方法

    2020-06-13 05:42:50孟兆海王懷君
    關(guān)鍵詞:重力梯度張量測(cè)線

    嚴(yán) 飛,李 達(dá),李 中,3,孟兆海,3,4,5,王懷君

    (1.海裝北京局駐天津地區(qū)某軍事代表室,天津 300131;2.天津航海儀器研究所,天津 300131;3.中船集團(tuán)航海保障技術(shù)實(shí)驗(yàn)室,天津 300131;4 中國(guó)科學(xué)院地質(zhì)與地球物理研究所,北京 100029;5 中國(guó)科學(xué)院油氣資源研究重點(diǎn)實(shí)驗(yàn)室,北京 100029)

    重力梯度是重力矢量的空間變化率,重力梯度儀是用于測(cè)量重力梯度的精密設(shè)備,基于Bell Aerospace公司提出的旋轉(zhuǎn)加速度計(jì)測(cè)量原理的重力梯度儀是迄今唯一實(shí)用的近地表動(dòng)態(tài)重力梯度儀[1-3]。重力梯度儀的核心部件是重力梯度敏感器,動(dòng)態(tài)測(cè)量時(shí),重力梯度敏感器安裝在慣性穩(wěn)定平臺(tái)上,穩(wěn)定在地理坐標(biāo)系下,測(cè)量地理坐標(biāo)系下的重力梯度。在測(cè)量過(guò)程中,飛機(jī)、艦船等載體的姿態(tài)變化,將引起載體質(zhì)量分布與重力梯度敏感器的相對(duì)位置發(fā)生變化。載體質(zhì)量對(duì)重力梯度敏感器產(chǎn)生的引力梯度將隨著載體姿態(tài)變化而產(chǎn)生變化,形成載體環(huán)境引力梯度變化[4,5]。隨著重力梯度儀精度水平不斷提高,在動(dòng)態(tài)測(cè)量中載體環(huán)境引力梯度變化已成為高精度重力梯度儀的一項(xiàng)重要誤差來(lái)源[6]。

    本文基于重力梯度張量在不同坐標(biāo)系下的轉(zhuǎn)換關(guān)系建立了載體環(huán)境引力梯度變化的解析模型,以此為基礎(chǔ)建立了回歸方程。通過(guò)Tikhonov正則化方法抑制了該方程解的不適定問(wèn)題,得到載體環(huán)境引力梯度在載體坐標(biāo)系下的數(shù)值解,以此進(jìn)行載體環(huán)境引力梯度變化改正,通過(guò)實(shí)測(cè)海試數(shù)據(jù)驗(yàn)證該方法的有效性,為動(dòng)態(tài)重力梯度儀研制提供參考。

    1 載體環(huán)境引力梯度變化模型

    1.1 重力梯度儀測(cè)量原理

    重力梯度是重力位的空間二階導(dǎo)數(shù),表征重力矢量的空間變化率。在地理坐標(biāo)系中,重力矢量可以分解為x,y,z三個(gè)方向上的三個(gè)分量中,每一分量沿平行于坐標(biāo)軸方向均有一個(gè)梯度。因此,重力梯度張量共有9個(gè)分量,可將其表示為:

    式中,Γij(i,j=x,y,z)為梯度張量的分量,表示重力分量gi在j方向上的變化率。

    旋轉(zhuǎn)加速度計(jì)式重力梯度儀主要由重力梯度敏感器和慣性穩(wěn)定平臺(tái)兩個(gè)關(guān)鍵組件構(gòu)成。重力梯度敏感器主要用于完成重力梯度張量水平分量的測(cè)量;慣性穩(wěn)定平臺(tái)則用于承載重力梯度敏感器,為重力梯度動(dòng)態(tài)測(cè)量提供穩(wěn)定的動(dòng)力學(xué)環(huán)境。如圖1所示,重力梯度敏感器基于加速度計(jì)位置差分的測(cè)量原理,通過(guò)機(jī)械旋轉(zhuǎn)的方式將旋轉(zhuǎn)中心處的重力梯度張量水平分量調(diào)制到系統(tǒng)旋轉(zhuǎn)頻率的二倍頻處,加速度計(jì)四路和信號(hào)與重力梯度張量水平分量之間的關(guān)系可表示為[7]:

    式中,a1、a2、a3、a4是四只加速度計(jì)敏感軸方向的加速度,l是加速度計(jì)檢測(cè)質(zhì)心到旋轉(zhuǎn)中心的距離,Γxx、Γyy、Γxy是測(cè)量平面內(nèi)對(duì)應(yīng)方向的重力梯度張量分量,單位是 E(1 E=10-9s-2),ω是旋轉(zhuǎn)圓盤的旋轉(zhuǎn)角速率。進(jìn)行動(dòng)態(tài)測(cè)量時(shí),慣性穩(wěn)定平臺(tái)采用三環(huán)固定指北半解析式控制方案,在隔離載體角運(yùn)動(dòng)的同時(shí),將重力梯度敏感器穩(wěn)定在地理坐標(biāo)系,為測(cè)量提供坐標(biāo)基準(zhǔn)。

    圖1 重力梯度敏感器測(cè)量原理示意圖Fig.1 Schematic diagram of measuring principle of gravity gradient sensor

    1.2 重力梯度坐標(biāo)變換

    依據(jù)重力勢(shì)函數(shù)梯度和方向?qū)?shù)關(guān)系,得到重力梯度張量在不同坐標(biāo)系下的變換公式為[8]:

    式中:Γn是重力梯度張量在地理坐標(biāo)系下的投影;Γb是重力梯度張量在載體坐標(biāo)系下的投影;Cbn是從載體坐標(biāo)系到地理坐標(biāo)系的坐標(biāo)變換矩陣。定義地理坐標(biāo)系為“東-北-天”地理坐標(biāo)系時(shí),依據(jù)歐拉轉(zhuǎn)動(dòng)原理,Cbn可表示為:

    式中,ψ、θ和γ分別為載體的航向角、俯仰角和橫滾角。

    重力場(chǎng)是一個(gè)保守場(chǎng),表明重力場(chǎng)是一個(gè)無(wú)源場(chǎng),且場(chǎng)線在空間內(nèi)不閉合,有[9]:

    式(5)表示重力場(chǎng)的旋度為零,說(shuō)明重力梯度張量矩陣具有對(duì)稱性,即:

    式(6)表示重力場(chǎng)的散度為零,意味著重力梯度張量矩陣的跡為零,即:

    因此,得到重力梯度的坐標(biāo)變換公式為:

    1.3 載體環(huán)境引力梯度變化模型

    式(10)可簡(jiǎn)寫為:

    式中,Rbn為從載體坐標(biāo)系到地理坐標(biāo)系的重力梯度張量列向量的坐標(biāo)變換矩陣。

    在動(dòng)態(tài)測(cè)量過(guò)程中重力梯度敏感器始終穩(wěn)定在地理坐標(biāo)系,載體質(zhì)量的分布與敏感器的相對(duì)位置發(fā)生變化,形成載體環(huán)境引力梯度變化。結(jié)合式(11),得到重力梯度儀載體環(huán)境引力梯度變化的解析方程為:

    式中,δTen是由載體環(huán)境引力梯度變化引起的水平分量重力梯度測(cè)量誤差;Teb是載體質(zhì)量在載體坐標(biāo)系下對(duì)重力梯度敏感器旋轉(zhuǎn)中心位置產(chǎn)生的引力梯度。

    利用實(shí)測(cè)艦船姿態(tài)信息結(jié)合艦船結(jié)構(gòu)引力梯度正演計(jì)算,滿載排水量為5000 t的艦船在測(cè)線上載體環(huán)境引力梯度變化量約為30 E量級(jí)。

    2 基于Tikhonov正則化的載體環(huán)境梯度誤差補(bǔ)償方法

    2.1 回歸方程

    依據(jù)式(12),可以寫出考慮載體環(huán)境引力梯度變化的重力梯度儀測(cè)量方程為:

    利用載體環(huán)境引力梯度變化和載體姿態(tài)的解析關(guān)系建立回歸方程,可用多元回歸的方式求取載體質(zhì)量在載體坐標(biāo)系下對(duì)重力梯度敏感器旋轉(zhuǎn)中心位置引起的引力梯度變化,從而結(jié)合載體姿態(tài)信息對(duì)測(cè)量數(shù)據(jù)進(jìn)行載體環(huán)境引力梯度變化改正,提高測(cè)量精度。

    2.2 解的不適定分析

    利用最小二乘法求解式(14),可得[11]:

    式中是x的估計(jì)解,即載體質(zhì)量在載體坐標(biāo)系下對(duì)重力梯度敏感器旋轉(zhuǎn)中心位置引起引力梯度的估計(jì)值。

    在不考慮回歸算子K的擾動(dòng)時(shí),有估計(jì)解的擾動(dòng)式為[12]:

    式中cond(K)是算子K在2范數(shù)意義下的條件數(shù),其計(jì)算式為:

    式(16)表明了回歸方程式(14)中的非線性量δ以cond(K)的倍數(shù)進(jìn)行放大,因此為抑制擾動(dòng)量δ對(duì)估計(jì)解的擾動(dòng),要求算子K的條件數(shù)越小越好。

    在進(jìn)行重力梯度測(cè)量作業(yè)時(shí),通常由飛機(jī)或艦船搭載重力梯度儀,在固定區(qū)域內(nèi)沿網(wǎng)格化測(cè)線測(cè)量。在測(cè)線中保持勻速直線運(yùn)動(dòng),最大程度降低載體加速度對(duì)于重力梯度測(cè)量的干擾。因此在測(cè)量過(guò)程中載體俯仰角和橫滾角的波動(dòng)較小,只有航向角能產(chǎn)生較大變化。而在飛機(jī)爬升和艦船轉(zhuǎn)彎等俯仰角、橫滾角波動(dòng)較大的航行狀態(tài)下,載體運(yùn)動(dòng)加速度干擾同樣很大,此時(shí)重力梯度儀無(wú)法實(shí)現(xiàn)高精度測(cè)量。

    某次重力梯度儀船載試驗(yàn)中兩條反向的測(cè)線船體俯仰角、橫滾角和航向角如圖2,3所示(定義為測(cè)線1#和測(cè)線2#),兩條測(cè)線的航跡線如圖4所示。從中可以看出測(cè)量過(guò)程中船體俯仰角波動(dòng)峰峰值為2°,橫滾角波動(dòng)峰峰值為2°,航向角波動(dòng)峰峰值為5°。計(jì)算cond(K)的結(jié)果如表1所示。

    圖2 重力梯度儀船載試驗(yàn)中測(cè)線1#船體航姿曲線Fig.2 Line 1# ship attitude curve in the shipboard test of gravity gradiometer

    圖3 重力梯度儀船載試驗(yàn)中測(cè)線2#船體航姿曲線Fig.3 Line 2# ship attitude curve in the shipboard test of gravity gradiometer

    圖4 重力梯度儀船載試驗(yàn)中測(cè)線1#和測(cè)線2#航跡線Fig.4 The track lines of line 1# and line 2# in the shipborne test of gravity gradiometer

    表1 測(cè)線1#和測(cè)線2#的cond(K)計(jì)算結(jié)果Tab.1 Calculation results of cond(K)of line 1# and line 2# respectively

    從表1中可以看出,利用兩條反向的測(cè)線1#和2#上的船體姿態(tài)數(shù)據(jù),cond(K)為3×103~5×103量級(jí),仍是病態(tài)的,即具有嚴(yán)重的復(fù)共線性特性。因此考慮利用兩條正交的測(cè)線進(jìn)行回歸分析,降低cond(K)的數(shù)值。

    引入一條與測(cè)線1#正交的測(cè)線,定義為測(cè)線3#,其船體俯仰角、橫滾角和航向角如圖5所示,測(cè)線1#和測(cè)線3#的航跡線如圖6所示,利用1#、2#和3#測(cè)線的數(shù)據(jù)計(jì)算cond(K)的結(jié)果如表2所示。

    圖5 重力梯度儀船載試驗(yàn)中測(cè)線3#船體航姿曲線Fig.5 Line 3# ship attitude curve in the shipboard test of gravity gradiometer

    圖6 重力梯度儀船載試驗(yàn)中測(cè)線1#和測(cè)線3#航跡線Fig.6 The track lines of line 1# and line 3# in the shipborne test of gravity gradiometer

    表2 測(cè)線1#和測(cè)線3#的cond(K)計(jì)算結(jié)果Tab.2 Calculation results of cond(K)of line 1# and line 3# respectively

    從表2中可以看出,即使利用3條正交的測(cè)線(1#、2#和3#)上的船體姿態(tài)數(shù)據(jù)計(jì)算cond(K),數(shù)值仍在1×103量級(jí),仍是病態(tài)的,在式(14)中會(huì)放大測(cè)量信號(hào)yδ中的擾動(dòng)量δ。同時(shí),回歸方程的擾動(dòng)量δ不僅包含測(cè)量噪聲,還包含真實(shí)的地表重力梯度信號(hào),回歸方程的擾動(dòng)量是擬合量的5倍以上,為防止“過(guò)擬合”的情況,必須選擇合理的正則化方法,以求取回歸方程的精確估計(jì)解。

    2.3 基于Tikhonov正則化的回歸方法

    目前解決最小二乘回歸病態(tài)問(wèn)題應(yīng)用較廣的方法是Tikhonov正則化方法[13],它通過(guò)求解如式(18)所示的最小化問(wèn)題獲得正則數(shù)值解。

    可通過(guò)式(19)得到。

    式中,I是單位矩陣,α是控制參數(shù)(α> 0)。其中(KTK+αI )-1KT是正則化算子。Tikhonov正則化方法的基本思想是犧牲估計(jì)解的無(wú)偏性將系數(shù)矩陣K的病態(tài)問(wèn)題轉(zhuǎn)化為良性問(wèn)題,使得正則數(shù)值解更加逼近于精確解x。

    2.4 正則化參數(shù)確定

    目前已有大量文獻(xiàn)研究了Tikhonov正則化方法中正則化算子的構(gòu)造和相關(guān)正則化參數(shù)的選取[14,15],本文采用半物理仿真的方式選擇控制參數(shù)α的最佳數(shù)值。通過(guò)引力梯度正演的方法仿真得到該方程的理論解x,利用實(shí)測(cè)重力梯度數(shù)據(jù)和實(shí)測(cè)姿態(tài)信息結(jié)合式(14)分別計(jì)算δ、K和yδ。定義誤差比為??刂茀?shù)α在[0.0001,100]范圍內(nèi)分別求取回歸方程的估計(jì)解,計(jì)算步長(zhǎng)為0.0001,并計(jì)算誤差比,結(jié)果如圖7所示。作為對(duì)比,再使用常規(guī)最小二乘的方法求取估計(jì)解,并計(jì)算誤差比。典型計(jì)算結(jié)果見(jiàn)表3。

    圖7 不同控制參數(shù)α?xí)r的誤差比Fig.7 Error ratio of different control parameter α

    表3 不同條件下回歸方程估計(jì)解的計(jì)算結(jié)果Tab.3 Calculation results of estimated solution of regression equation under different conditions

    表3 不同條件下回歸方程估計(jì)解的計(jì)算結(jié)果Tab.3 Calculation results of estimated solution of regression equation under different conditions

    估計(jì)值/E b Γ b xx Γ b xy Γ b xz Γ b yy Γ 誤差比/%yz理論值 26.2225.874.01 -5.332.88 / α估=計(jì)解 860.7225.563.04829.86 -50.7823.220.0001估α計(jì)=解 15.9225.443.50 -13.909.681.989 α估計(jì)=解 15.3725.280.94 -14.594.331.9926 α估=計(jì)解 11.9422.27-0.71 -11.900.6515.141000最估小計(jì)二解乘 878.1725.552.71847.26 -52.6924.69

    結(jié)合圖7和表3中可知,與較常規(guī)最小二乘方法比,Tikhonov正則化的方法只要將控制參數(shù)α選取在9~26之間,能將補(bǔ)償量的誤差比從24.69%降低至2%以內(nèi),且相差不大。因此在本研究中選擇α= 15。

    3 試驗(yàn)驗(yàn)證

    3.1 載體環(huán)境引力梯度估計(jì)

    重力梯度儀船載試驗(yàn)中測(cè)線1# 和2# 兩條重復(fù)線的測(cè)量數(shù)據(jù)未進(jìn)行載體環(huán)境引力梯度補(bǔ)償時(shí),重力梯度測(cè)量信號(hào)結(jié)果如圖8和9所示。結(jié)合載體姿態(tài)數(shù)據(jù)使用式(19)進(jìn)行載體環(huán)境梯度估計(jì),分別使用常規(guī)最小二乘法和Tikhonov正則化方法進(jìn)行回歸估計(jì),計(jì)算結(jié)果如表4所示。

    圖8 重力梯度 Γuv分量原始輸出信號(hào)Fig.8 Original output signal of gravity gradient component Γuv

    圖9 重力梯度 Γxy分量原始輸出信號(hào)Fig.9 Original output signal of gravity gradient component Γxy

    表4 載體環(huán)境引力梯度估計(jì)結(jié)果Tab.4 Calculation results of carrier environmental gradient

    3.1 重力梯度測(cè)量結(jié)果及精度評(píng)價(jià)

    將Tikhonov正則化法和常規(guī)最小二乘法得到的載體環(huán)境梯度估計(jì)值代入式(14),進(jìn)行重力梯度信號(hào)載體環(huán)境引力梯度變化補(bǔ)償,使用Tikhonov正則化法補(bǔ)償后重力梯度測(cè)量信號(hào)結(jié)果如圖10和11所示,使用最小二乘法補(bǔ)償后重力梯度測(cè)量信號(hào)結(jié)果如圖12和13所示,分別計(jì)算補(bǔ)償前后重力梯度信號(hào)內(nèi)符合精度,如表5所示。

    圖10 Tikhonov正則化法載體環(huán)境引力梯度變化補(bǔ)償后重力梯度分量信號(hào)Fig.10. component signal of gravity gradient after compensation of environmental gravity gradient change by Tikhonov regularization method

    圖11 Tikhonov正則化法載體環(huán)境引力梯度變化補(bǔ)償后重力梯度 Γxy分量信號(hào)Fig.11 Γxy component signal of gravity gradient after compensation of environmental gravity gradient change by Tikhonov regularization method

    圖12 最小二乘法載體環(huán)境引力梯度變化補(bǔ)償后 重力梯度 Γxy分量信號(hào)Fig.12 Γuv component signal of gravity gradient after compensation of gravity gradient change of carrier environment by least square method

    圖13 最小二乘法法載體環(huán)境引力梯度變化補(bǔ)償后 重力梯度 Γxy分量信號(hào)Fig.13.Γxy component signal of gravity gradient after compensation of gravity gradient change of carrier environment by least square method

    表5 載體環(huán)境引力梯度變化補(bǔ)償前后內(nèi)符合中誤差對(duì)比Tab.5 Comparison of the mean square error of the internal coincidence before and after compensation for the change of the gravity gradient of the carrier environment

    4 結(jié) 論

    本文針對(duì)重力梯度儀動(dòng)態(tài)測(cè)量過(guò)程中載體環(huán)境引力梯度干擾的問(wèn)題,建立了動(dòng)態(tài)測(cè)量過(guò)程中載體環(huán)境引力梯度變化的傳遞模型,以此為基礎(chǔ)建立了回歸方程。通過(guò)Tikhonov正則化方法抑制了該方程解的不適定問(wèn)題,利用半實(shí)物仿真的方法得到最優(yōu)控制參數(shù)α的取值范圍,使補(bǔ)償量的誤差控制在2%。經(jīng)實(shí)測(cè)船載數(shù)據(jù)驗(yàn)證,該方法切實(shí)可行,補(bǔ)償后內(nèi)符合中誤差較原始重力梯度測(cè)量信號(hào)中誤差分別降低19 E和 21 E,較常規(guī)最小二乘法進(jìn)一步降低7 E和3 E,補(bǔ)償后重力梯度儀兩路測(cè)量信號(hào)內(nèi)符合中誤差均達(dá)到優(yōu)于10 E的精度水平。

    猜你喜歡
    重力梯度張量測(cè)線
    極地海洋多波束測(cè)量測(cè)線布設(shè)系統(tǒng)設(shè)計(jì)及實(shí)現(xiàn)
    基于動(dòng)態(tài)規(guī)劃的多波束測(cè)線布設(shè)模型
    偶數(shù)階張量core逆的性質(zhì)和應(yīng)用
    四元數(shù)張量方程A*NX=B 的通解
    擴(kuò)散張量成像MRI 在CO中毒后遲發(fā)腦病中的應(yīng)用
    旋轉(zhuǎn)加速度計(jì)重力梯度儀標(biāo)定方法
    利用地形數(shù)據(jù)計(jì)算重力梯度張量的直接積分法
    星載重力梯度儀的研究發(fā)展
    工程中張量概念的思考
    河南科技(2014年19期)2014-02-27 14:15:33
    隧洞中雷達(dá)探測(cè)地質(zhì)構(gòu)造的測(cè)線布置與三維地質(zhì)解譯
    不卡视频在线观看欧美| 人妻 亚洲 视频| 在线观看人妻少妇| 欧美97在线视频| 晚上一个人看的免费电影| 国产爽快片一区二区三区| 高清欧美精品videossex| 性色av一级| 国产一区有黄有色的免费视频| 99热这里只有是精品50| 成年免费大片在线观看| 在线观看人妻少妇| 国产精品一区二区性色av| 汤姆久久久久久久影院中文字幕| 边亲边吃奶的免费视频| 天天躁夜夜躁狠狠久久av| 七月丁香在线播放| 综合色丁香网| 国产在视频线精品| 人妻少妇偷人精品九色| 成年av动漫网址| 午夜视频国产福利| 日韩av不卡免费在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 久久国产亚洲av麻豆专区| 成年免费大片在线观看| 日韩成人伦理影院| 日产精品乱码卡一卡2卡三| 亚洲欧洲日产国产| 秋霞在线观看毛片| 中文资源天堂在线| 亚洲一区二区三区欧美精品| 欧美另类一区| 免费黄色在线免费观看| 亚洲无线观看免费| 观看av在线不卡| 天堂8中文在线网| 亚洲人与动物交配视频| 97超碰精品成人国产| 成人毛片60女人毛片免费| 一本—道久久a久久精品蜜桃钙片| 天天躁日日操中文字幕| 亚洲精品一区蜜桃| 免费观看在线日韩| 夫妻午夜视频| 韩国高清视频一区二区三区| 尾随美女入室| 久久99蜜桃精品久久| 搡老乐熟女国产| 国产黄色视频一区二区在线观看| 亚洲欧美成人精品一区二区| 中文欧美无线码| 青春草亚洲视频在线观看| 少妇人妻精品综合一区二区| 熟女av电影| 国产精品蜜桃在线观看| 三级经典国产精品| 大码成人一级视频| 久久精品夜色国产| 亚洲在久久综合| 精品人妻视频免费看| 超碰av人人做人人爽久久| 少妇 在线观看| 亚洲av中文字字幕乱码综合| 国产精品福利在线免费观看| 亚洲欧美日韩卡通动漫| 中文字幕人妻熟人妻熟丝袜美| 久久鲁丝午夜福利片| 欧美日韩视频精品一区| 亚洲高清免费不卡视频| 成人美女网站在线观看视频| 男的添女的下面高潮视频| 国产亚洲欧美精品永久| 精品人妻熟女av久视频| 人人妻人人爽人人添夜夜欢视频 | 美女国产视频在线观看| 久久国内精品自在自线图片| 久久这里有精品视频免费| av网站免费在线观看视频| 有码 亚洲区| 国产精品一区www在线观看| 男人添女人高潮全过程视频| 久久久久久九九精品二区国产| 日韩av免费高清视频| 人妻少妇偷人精品九色| 欧美日本视频| 国产一区二区在线观看日韩| 久久 成人 亚洲| 国产爽快片一区二区三区| 国产精品一区二区在线观看99| 一区二区三区精品91| 久久久久久久久大av| 免费久久久久久久精品成人欧美视频 | 午夜福利高清视频| 视频区图区小说| 在线观看一区二区三区| 欧美日韩综合久久久久久| 日本免费在线观看一区| 日韩一本色道免费dvd| 美女主播在线视频| 亚洲国产最新在线播放| 久久久a久久爽久久v久久| 最近中文字幕2019免费版| 国产爱豆传媒在线观看| 免费人妻精品一区二区三区视频| 亚洲欧美成人综合另类久久久| 久久国产亚洲av麻豆专区| 又爽又黄a免费视频| 精品国产露脸久久av麻豆| 免费不卡的大黄色大毛片视频在线观看| 国产av国产精品国产| 麻豆成人午夜福利视频| 中文字幕免费在线视频6| 99国产精品免费福利视频| 久久久久久九九精品二区国产| 亚洲天堂av无毛| 国产亚洲91精品色在线| 99热这里只有是精品在线观看| 日本黄色日本黄色录像| 各种免费的搞黄视频| 中文乱码字字幕精品一区二区三区| av专区在线播放| 国产高潮美女av| 免费播放大片免费观看视频在线观看| 亚洲国产欧美在线一区| 少妇 在线观看| 亚洲在久久综合| 搡女人真爽免费视频火全软件| 欧美日韩国产mv在线观看视频 | 多毛熟女@视频| 99热这里只有是精品在线观看| 久久久久久久久久成人| 久久精品国产自在天天线| 一级毛片我不卡| 欧美97在线视频| 欧美亚洲 丝袜 人妻 在线| 美女xxoo啪啪120秒动态图| 晚上一个人看的免费电影| 大陆偷拍与自拍| 在线观看三级黄色| 一本久久精品| 免费看日本二区| av在线老鸭窝| 尾随美女入室| 午夜福利网站1000一区二区三区| 色视频在线一区二区三区| 国产久久久一区二区三区| 亚洲电影在线观看av| 黑人猛操日本美女一级片| 久久精品久久久久久久性| 大香蕉久久网| 黑人猛操日本美女一级片| 国产一级毛片在线| 一级毛片黄色毛片免费观看视频| 久久久精品94久久精品| 亚洲精品一区蜜桃| 成人毛片a级毛片在线播放| 亚洲自偷自拍三级| 国产综合精华液| 国产精品99久久99久久久不卡 | 日韩欧美 国产精品| 激情 狠狠 欧美| 久久久久久久大尺度免费视频| 色网站视频免费| 男女边吃奶边做爰视频| 久久久久久人妻| 欧美xxⅹ黑人| 国产精品女同一区二区软件| 国产精品一二三区在线看| 老师上课跳d突然被开到最大视频| 国产精品麻豆人妻色哟哟久久| 免费看光身美女| 日本vs欧美在线观看视频 | 国产av一区二区精品久久 | 欧美精品亚洲一区二区| 网址你懂的国产日韩在线| 丰满人妻一区二区三区视频av| 国产精品99久久99久久久不卡 | 亚洲欧美日韩卡通动漫| 日韩一区二区视频免费看| 在线观看免费日韩欧美大片 | 国产精品女同一区二区软件| 日本欧美视频一区| 少妇被粗大猛烈的视频| 大香蕉97超碰在线| 成人高潮视频无遮挡免费网站| 黑人高潮一二区| 五月天丁香电影| videos熟女内射| 免费看不卡的av| 亚洲精品aⅴ在线观看| 欧美一区二区亚洲| 尤物成人国产欧美一区二区三区| 国产黄频视频在线观看| 麻豆成人午夜福利视频| 91精品一卡2卡3卡4卡| 男女啪啪激烈高潮av片| 99久久精品热视频| 欧美国产精品一级二级三级 | 九色成人免费人妻av| 18禁在线无遮挡免费观看视频| 纯流量卡能插随身wifi吗| 精品少妇黑人巨大在线播放| 免费高清在线观看视频在线观看| 久久ye,这里只有精品| 熟妇人妻不卡中文字幕| 成人免费观看视频高清| 日韩中字成人| 久久久久久久精品精品| 久久久a久久爽久久v久久| 亚洲精品自拍成人| 欧美日韩视频精品一区| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲欧美中文字幕日韩二区| 伦精品一区二区三区| 自拍偷自拍亚洲精品老妇| 国产精品人妻久久久影院| 国产亚洲精品久久久com| 国产午夜精品一二区理论片| 日韩大片免费观看网站| 极品教师在线视频| 久久久亚洲精品成人影院| 女人十人毛片免费观看3o分钟| 夜夜爽夜夜爽视频| 日韩一本色道免费dvd| 亚洲av电影在线观看一区二区三区| 国产乱来视频区| 91狼人影院| 国产 一区精品| 激情五月婷婷亚洲| 韩国av在线不卡| 美女脱内裤让男人舔精品视频| 丰满迷人的少妇在线观看| av黄色大香蕉| 哪个播放器可以免费观看大片| 成年人午夜在线观看视频| 不卡视频在线观看欧美| 久久影院123| 国产在线一区二区三区精| 午夜福利在线在线| 中文乱码字字幕精品一区二区三区| 大香蕉久久网| 色视频在线一区二区三区| 大又大粗又爽又黄少妇毛片口| 五月伊人婷婷丁香| 韩国av在线不卡| www.av在线官网国产| 男女下面进入的视频免费午夜| 久久6这里有精品| 久久久久久人妻| 这个男人来自地球电影免费观看 | 欧美性感艳星| 国产又色又爽无遮挡免| 精品久久国产蜜桃| 久久久久久久亚洲中文字幕| av国产免费在线观看| av女优亚洲男人天堂| 色综合色国产| av国产精品久久久久影院| 美女中出高潮动态图| 黄片wwwwww| 国产成人aa在线观看| 国产熟女欧美一区二区| 18禁在线无遮挡免费观看视频| 色5月婷婷丁香| 亚洲精品乱久久久久久| 亚洲精品中文字幕在线视频 | 欧美国产精品一级二级三级 | 下体分泌物呈黄色| 国产深夜福利视频在线观看| 九九在线视频观看精品| 亚洲精品国产色婷婷电影| 欧美日韩亚洲高清精品| 国产成人精品一,二区| 极品少妇高潮喷水抽搐| 在线播放无遮挡| 天堂中文最新版在线下载| 亚洲精品乱久久久久久| 亚洲综合精品二区| 成人综合一区亚洲| 久久久久久久久久人人人人人人| 欧美精品一区二区免费开放| 观看美女的网站| 久久久久国产精品人妻一区二区| av专区在线播放| 久久久亚洲精品成人影院| 欧美日韩视频高清一区二区三区二| 成人无遮挡网站| 亚洲美女搞黄在线观看| 久久久a久久爽久久v久久| 少妇的逼水好多| 永久免费av网站大全| 日韩精品有码人妻一区| 青青草视频在线视频观看| 日韩 亚洲 欧美在线| 国产亚洲一区二区精品| 中国美白少妇内射xxxbb| 高清黄色对白视频在线免费看 | 日韩人妻高清精品专区| 欧美精品国产亚洲| 男女边吃奶边做爰视频| 久久久色成人| 夜夜爽夜夜爽视频| 日韩中字成人| 国产美女午夜福利| 纯流量卡能插随身wifi吗| 成人综合一区亚洲| 男人舔奶头视频| 国产精品久久久久成人av| 夫妻性生交免费视频一级片| 18禁在线无遮挡免费观看视频| 亚洲最大成人中文| 国产综合精华液| 亚洲精华国产精华液的使用体验| 国产精品免费大片| 午夜免费观看性视频| 亚洲精品乱久久久久久| 在线观看一区二区三区激情| 777米奇影视久久| 免费黄频网站在线观看国产| 亚洲av男天堂| 丝袜喷水一区| 亚洲一级一片aⅴ在线观看| av线在线观看网站| 99热网站在线观看| av黄色大香蕉| 九九久久精品国产亚洲av麻豆| 中文字幕制服av| 欧美日韩视频精品一区| 18禁动态无遮挡网站| 亚洲熟女精品中文字幕| 超碰97精品在线观看| 成人毛片a级毛片在线播放| 亚洲av免费高清在线观看| 精品国产露脸久久av麻豆| 国产免费一区二区三区四区乱码| 99精国产麻豆久久婷婷| 丰满少妇做爰视频| 免费观看av网站的网址| 日韩亚洲欧美综合| 成人毛片a级毛片在线播放| 在现免费观看毛片| 日本黄大片高清| 最近手机中文字幕大全| 国产精品一区二区在线观看99| 久久 成人 亚洲| 免费大片黄手机在线观看| 18禁在线播放成人免费| 中国美白少妇内射xxxbb| 搡老乐熟女国产| 国产精品人妻久久久久久| 激情五月婷婷亚洲| 日日摸夜夜添夜夜添av毛片| 我要看日韩黄色一级片| 欧美激情国产日韩精品一区| 午夜免费鲁丝| 国产探花极品一区二区| 一二三四中文在线观看免费高清| 成人高潮视频无遮挡免费网站| 久久99热这里只频精品6学生| 精品亚洲成国产av| 只有这里有精品99| 亚洲精品视频女| 亚洲国产av新网站| 亚洲欧美中文字幕日韩二区| 成人高潮视频无遮挡免费网站| 亚洲av.av天堂| 下体分泌物呈黄色| 国产欧美亚洲国产| 美女国产视频在线观看| 精品人妻熟女av久视频| 一个人看的www免费观看视频| 国产精品三级大全| 国产在线男女| 成人影院久久| 亚洲无线观看免费| av在线播放精品| 国产在线免费精品| av卡一久久| 久久女婷五月综合色啪小说| 免费观看的影片在线观看| 精品国产露脸久久av麻豆| 观看美女的网站| 久久这里有精品视频免费| 免费观看a级毛片全部| 97热精品久久久久久| 国产精品伦人一区二区| 久久久久久久久久人人人人人人| 精品一区在线观看国产| 国产精品人妻久久久久久| 一级片'在线观看视频| 麻豆精品久久久久久蜜桃| 久久av网站| 各种免费的搞黄视频| 中文字幕亚洲精品专区| 永久免费av网站大全| 一本一本综合久久| 22中文网久久字幕| 韩国av在线不卡| 91aial.com中文字幕在线观看| 欧美成人精品欧美一级黄| 五月伊人婷婷丁香| 极品少妇高潮喷水抽搐| 大片电影免费在线观看免费| 97超碰精品成人国产| 久久久久人妻精品一区果冻| 亚洲第一av免费看| 美女主播在线视频| 日韩亚洲欧美综合| 亚洲久久久国产精品| 亚洲内射少妇av| av黄色大香蕉| 精华霜和精华液先用哪个| 国语对白做爰xxxⅹ性视频网站| 亚洲人成网站在线观看播放| 成人特级av手机在线观看| 久久av网站| 日本av手机在线免费观看| 女人久久www免费人成看片| 国模一区二区三区四区视频| 两个人的视频大全免费| 99热6这里只有精品| 欧美日韩在线观看h| 亚洲精品国产av蜜桃| 亚洲国产精品专区欧美| 成人美女网站在线观看视频| 国产精品欧美亚洲77777| 国产伦理片在线播放av一区| 午夜福利视频精品| 亚洲国产高清在线一区二区三| av女优亚洲男人天堂| 婷婷色麻豆天堂久久| 性色av一级| 精品国产三级普通话版| 亚洲国产精品国产精品| 九色成人免费人妻av| 97超碰精品成人国产| 国产黄片视频在线免费观看| 久久女婷五月综合色啪小说| 成年人午夜在线观看视频| freevideosex欧美| 女的被弄到高潮叫床怎么办| 国产成人91sexporn| 亚洲国产毛片av蜜桃av| 国产亚洲最大av| 欧美+日韩+精品| 欧美成人a在线观看| 国产av码专区亚洲av| 简卡轻食公司| 一区在线观看完整版| 尤物成人国产欧美一区二区三区| 精品久久久久久久久av| 在线精品无人区一区二区三 | 精品少妇黑人巨大在线播放| 性色avwww在线观看| 男的添女的下面高潮视频| 色视频www国产| 看非洲黑人一级黄片| 老女人水多毛片| 97在线视频观看| 日韩国内少妇激情av| 成人无遮挡网站| 中文天堂在线官网| 久热久热在线精品观看| 美女主播在线视频| 色视频www国产| 国产乱人视频| 女性被躁到高潮视频| 国产av码专区亚洲av| 久久久久精品性色| 国产精品无大码| 国产在线一区二区三区精| 国产国拍精品亚洲av在线观看| 国产又色又爽无遮挡免| 午夜福利视频精品| 亚洲精品第二区| 午夜日本视频在线| 免费黄网站久久成人精品| 五月玫瑰六月丁香| 尾随美女入室| 丰满少妇做爰视频| 亚洲精品自拍成人| 一级黄片播放器| 男的添女的下面高潮视频| 国产精品嫩草影院av在线观看| av.在线天堂| 少妇 在线观看| 日韩伦理黄色片| 国产高清三级在线| 一级二级三级毛片免费看| 性高湖久久久久久久久免费观看| 日韩中文字幕视频在线看片 | 久久久久久伊人网av| 狠狠精品人妻久久久久久综合| 亚洲欧美精品专区久久| 嘟嘟电影网在线观看| 十分钟在线观看高清视频www | 国产成人精品一,二区| 亚洲av成人精品一二三区| 午夜精品国产一区二区电影| 亚洲欧洲国产日韩| 久久婷婷青草| 久久6这里有精品| 大话2 男鬼变身卡| 色视频www国产| 久久久久久久久大av| 欧美变态另类bdsm刘玥| 啦啦啦中文免费视频观看日本| 99视频精品全部免费 在线| 日韩 亚洲 欧美在线| 99热国产这里只有精品6| 亚洲av成人精品一区久久| 搡老乐熟女国产| 久久久久久久久久人人人人人人| 亚洲精品第二区| 日本wwww免费看| 亚洲av中文字字幕乱码综合| 毛片一级片免费看久久久久| 男女啪啪激烈高潮av片| 99久久人妻综合| 一级毛片电影观看| 夫妻性生交免费视频一级片| 国产中年淑女户外野战色| 亚洲国产av新网站| 欧美高清成人免费视频www| 国产免费一区二区三区四区乱码| 日韩av在线免费看完整版不卡| 内地一区二区视频在线| 精品久久国产蜜桃| 国产黄片美女视频| 亚洲欧美日韩另类电影网站 | 国产v大片淫在线免费观看| av又黄又爽大尺度在线免费看| 国产毛片在线视频| 99久久精品国产国产毛片| 国产精品无大码| 国产日韩欧美在线精品| av在线播放精品| 永久免费av网站大全| 嘟嘟电影网在线观看| 日韩伦理黄色片| 高清在线视频一区二区三区| 久久久久视频综合| 精品久久久噜噜| 久久久国产一区二区| av专区在线播放| 舔av片在线| 高清午夜精品一区二区三区| 亚洲成人一二三区av| 各种免费的搞黄视频| 十八禁网站网址无遮挡 | 精品一区二区三卡| 一二三四中文在线观看免费高清| 日韩一本色道免费dvd| 直男gayav资源| 国产视频首页在线观看| 高清午夜精品一区二区三区| 一级av片app| 九色成人免费人妻av| 三级经典国产精品| 人妻系列 视频| 黑人高潮一二区| 看免费成人av毛片| 国产淫片久久久久久久久| 成人国产麻豆网| 欧美变态另类bdsm刘玥| 18禁裸乳无遮挡免费网站照片| 欧美精品一区二区免费开放| 青春草视频在线免费观看| 久久精品久久精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 直男gayav资源| 欧美性感艳星| 永久免费av网站大全| 丝瓜视频免费看黄片| 青春草亚洲视频在线观看| 欧美xxⅹ黑人| 成年av动漫网址| 久久久久人妻精品一区果冻| 成年美女黄网站色视频大全免费 | 亚洲欧美日韩另类电影网站 | 高清毛片免费看| 中国国产av一级| 成年美女黄网站色视频大全免费 | 欧美3d第一页| 国产男女内射视频| 欧美日韩综合久久久久久| 成人二区视频| 亚洲国产欧美人成| 国产高清不卡午夜福利| 一级av片app| 亚洲国产欧美人成| 国产有黄有色有爽视频| 一级av片app| 国国产精品蜜臀av免费| 大陆偷拍与自拍| 超碰97精品在线观看| 亚洲欧洲日产国产| 国产av码专区亚洲av| 女的被弄到高潮叫床怎么办| 男女国产视频网站| 午夜福利在线观看免费完整高清在| 亚洲人成网站在线播| 成人漫画全彩无遮挡| 亚洲精品aⅴ在线观看| 97在线人人人人妻| 欧美日韩视频高清一区二区三区二| 日韩欧美 国产精品| 高清毛片免费看| 亚洲欧洲日产国产| 九草在线视频观看| 97超视频在线观看视频| 久久久久国产精品人妻一区二区| 国产精品偷伦视频观看了| 精品99又大又爽又粗少妇毛片|