• <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ì)解譯
    免费不卡黄色视频| 亚洲av电影在线进入| 韩国精品一区二区三区| 亚洲色图综合在线观看| 麻豆av在线久日| 日本91视频免费播放| 另类精品久久| 亚洲欧美成人综合另类久久久| 久久人人爽人人片av| 国产精品一区二区在线观看99| 一区二区三区四区激情视频| 高清视频免费观看一区二区| 久久久久久久久免费视频了| 欧美老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 丰满人妻熟妇乱又伦精品不卡| 欧美亚洲 丝袜 人妻 在线| 赤兔流量卡办理| 99热全是精品| 久久久久久久精品精品| 亚洲国产av新网站| 国产又爽黄色视频| 欧美国产精品一级二级三级| 亚洲人成网站在线观看播放| 久久久精品94久久精品| 日韩欧美一区视频在线观看| 久9热在线精品视频| 另类精品久久| 国产欧美日韩综合在线一区二区| 欧美日韩亚洲高清精品| 欧美日韩一级在线毛片| 好男人电影高清在线观看| 日本午夜av视频| 美女视频免费永久观看网站| 国产欧美日韩综合在线一区二区| 91字幕亚洲| 亚洲av男天堂| 久久久国产欧美日韩av| 日韩大片免费观看网站| 欧美精品高潮呻吟av久久| 一级片'在线观看视频| 青春草视频在线免费观看| 自拍欧美九色日韩亚洲蝌蚪91| av天堂在线播放| 老司机午夜十八禁免费视频| 精品人妻一区二区三区麻豆| 午夜老司机福利片| 亚洲av在线观看美女高潮| 中文字幕另类日韩欧美亚洲嫩草| 国语对白做爰xxxⅹ性视频网站| 一级毛片黄色毛片免费观看视频| 亚洲精品日本国产第一区| 中国美女看黄片| 女人被躁到高潮嗷嗷叫费观| 少妇人妻久久综合中文| 天天影视国产精品| 欧美精品av麻豆av| 久久久久精品人妻al黑| 美女福利国产在线| 免费观看av网站的网址| 亚洲精品一卡2卡三卡4卡5卡 | av网站在线播放免费| 考比视频在线观看| 校园人妻丝袜中文字幕| 亚洲中文av在线| 女性被躁到高潮视频| 国产成人av激情在线播放| 90打野战视频偷拍视频| 欧美 亚洲 国产 日韩一| 日韩人妻精品一区2区三区| 免费看av在线观看网站| 免费不卡黄色视频| 成人午夜精彩视频在线观看| 制服人妻中文乱码| 国产一区二区三区综合在线观看| 91九色精品人成在线观看| 国产成人精品无人区| 久久久久久免费高清国产稀缺| 国产成人一区二区三区免费视频网站 | 精品久久久久久电影网| 美女午夜性视频免费| 人人澡人人妻人| av在线老鸭窝| 精品久久蜜臀av无| 久久人人97超碰香蕉20202| 天天操日日干夜夜撸| 亚洲精品第二区| 18禁观看日本| 国产在线一区二区三区精| 国产淫语在线视频| 久久鲁丝午夜福利片| 欧美精品啪啪一区二区三区 | 欧美成人午夜精品| 午夜老司机福利片| 成人午夜精彩视频在线观看| 国产av国产精品国产| 性色av乱码一区二区三区2| 久久久精品区二区三区| 欧美黄色淫秽网站| 尾随美女入室| 久久天堂一区二区三区四区| 国产精品av久久久久免费| 久久亚洲精品不卡| 精品人妻1区二区| 韩国高清视频一区二区三区| 五月天丁香电影| 首页视频小说图片口味搜索 | 色综合欧美亚洲国产小说| 欧美激情极品国产一区二区三区| 国产老妇伦熟女老妇高清| 国产国语露脸激情在线看| 国产午夜精品一二区理论片| 久久精品国产综合久久久| 91国产中文字幕| 巨乳人妻的诱惑在线观看| 国产成人免费无遮挡视频| 久久国产亚洲av麻豆专区| 国产男女内射视频| 桃花免费在线播放| 欧美日韩视频精品一区| 久久久久久久久久久久大奶| 免费观看人在逋| 国产精品久久久人人做人人爽| 王馨瑶露胸无遮挡在线观看| 纵有疾风起免费观看全集完整版| 国产亚洲精品第一综合不卡| 免费一级毛片在线播放高清视频 | 99久久人妻综合| 精品一区二区三区av网在线观看 | cao死你这个sao货| 新久久久久国产一级毛片| 脱女人内裤的视频| 超色免费av| 一区福利在线观看| 无限看片的www在线观看| 美女国产高潮福利片在线看| 国产免费又黄又爽又色| 成人手机av| 国语对白做爰xxxⅹ性视频网站| 人妻 亚洲 视频| av福利片在线| 国产精品免费视频内射| 久久精品人人爽人人爽视色| 日本猛色少妇xxxxx猛交久久| 国产无遮挡羞羞视频在线观看| av片东京热男人的天堂| www.av在线官网国产| 黄网站色视频无遮挡免费观看| 亚洲精品国产一区二区精华液| 免费观看av网站的网址| 丝袜喷水一区| 9191精品国产免费久久| 久久这里只有精品19| 欧美激情极品国产一区二区三区| 国产成人免费无遮挡视频| 美女大奶头黄色视频| 亚洲第一av免费看| 亚洲欧美一区二区三区久久| 少妇人妻 视频| 看免费成人av毛片| 亚洲精品在线美女| 两个人看的免费小视频| 又大又爽又粗| 免费看不卡的av| 一级黄色大片毛片| 国产视频首页在线观看| 18在线观看网站| 国产伦理片在线播放av一区| 午夜两性在线视频| 成人影院久久| 国产精品国产三级专区第一集| 宅男免费午夜| 蜜桃国产av成人99| 青青草视频在线视频观看| 免费观看a级毛片全部| 欧美人与善性xxx| 国产精品久久久久久精品古装| 日本av手机在线免费观看| 中文欧美无线码| 高清视频免费观看一区二区| 亚洲国产欧美日韩在线播放| av国产精品久久久久影院| 99热国产这里只有精品6| 成人亚洲精品一区在线观看| 十八禁高潮呻吟视频| 一区二区三区乱码不卡18| 成人国语在线视频| 男女之事视频高清在线观看 | 人人妻人人爽人人添夜夜欢视频| 天堂俺去俺来也www色官网| 久久精品久久精品一区二区三区| 99香蕉大伊视频| 亚洲成色77777| 国产野战对白在线观看| 欧美国产精品va在线观看不卡| 中文精品一卡2卡3卡4更新| 高清黄色对白视频在线免费看| 日韩欧美一区视频在线观看| 亚洲三区欧美一区| 午夜免费鲁丝| 一区二区av电影网| 两性夫妻黄色片| 人体艺术视频欧美日本| 人体艺术视频欧美日本| 免费看不卡的av| 欧美+亚洲+日韩+国产| 午夜福利影视在线免费观看| av又黄又爽大尺度在线免费看| 母亲3免费完整高清在线观看| 免费看不卡的av| 999久久久国产精品视频| 亚洲欧美清纯卡通| 多毛熟女@视频| 如日韩欧美国产精品一区二区三区| 国产亚洲精品第一综合不卡| 成人国产一区最新在线观看 | 90打野战视频偷拍视频| 肉色欧美久久久久久久蜜桃| 99精国产麻豆久久婷婷| 大型av网站在线播放| 最新的欧美精品一区二区| 亚洲欧美一区二区三区黑人| 久久精品国产亚洲av涩爱| 国产免费又黄又爽又色| 老司机靠b影院| 国产老妇伦熟女老妇高清| 看免费av毛片| 国产男人的电影天堂91| 波多野结衣一区麻豆| 国产伦人伦偷精品视频| 国产成人一区二区在线| 97在线人人人人妻| 日韩一区二区三区影片| 亚洲国产最新在线播放| www.熟女人妻精品国产| 亚洲av成人精品一二三区| 国产精品久久久久久人妻精品电影 | 久久久欧美国产精品| 韩国精品一区二区三区| 久久久久国产精品人妻一区二区| 日日爽夜夜爽网站| 国产精品三级大全| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色麻豆天堂久久| 久久国产精品大桥未久av| 久久天堂一区二区三区四区| 国产野战对白在线观看| 国产无遮挡羞羞视频在线观看| 国产精品熟女久久久久浪| 久久青草综合色| 亚洲国产中文字幕在线视频| 18禁黄网站禁片午夜丰满| 亚洲精品自拍成人| 色网站视频免费| 尾随美女入室| 国产xxxxx性猛交| 自线自在国产av| 亚洲欧美激情在线| 午夜福利影视在线免费观看| 黄色a级毛片大全视频| 久久久国产一区二区| 又大又黄又爽视频免费| 麻豆av在线久日| 久久人妻熟女aⅴ| 啦啦啦 在线观看视频| 另类亚洲欧美激情| 搡老岳熟女国产| 国产深夜福利视频在线观看| 久久人人爽人人片av| 老司机深夜福利视频在线观看 | 日韩中文字幕视频在线看片| 中文字幕人妻熟女乱码| 18在线观看网站| 国产女主播在线喷水免费视频网站| 国产一区二区激情短视频 | 男人舔女人的私密视频| 国产视频一区二区在线看| 热99国产精品久久久久久7| 精品少妇内射三级| 国产精品一区二区免费欧美 | 91国产中文字幕| 久久免费观看电影| 少妇裸体淫交视频免费看高清 | 大片电影免费在线观看免费| 亚洲国产av新网站| 美女主播在线视频| 久久性视频一级片| 亚洲一区二区三区欧美精品| 久久精品国产亚洲av高清一级| 熟女av电影| 成人国产一区最新在线观看 | 国产亚洲av片在线观看秒播厂| 美女国产高潮福利片在线看| 一本综合久久免费| 在线观看免费高清a一片| 亚洲国产精品一区三区| 亚洲欧洲日产国产| 国产一区亚洲一区在线观看| 日韩一卡2卡3卡4卡2021年| 18禁黄网站禁片午夜丰满| 午夜av观看不卡| 成人免费观看视频高清| 人人妻人人爽人人添夜夜欢视频| 国产午夜精品一二区理论片| 一区福利在线观看| 久久狼人影院| 久久精品熟女亚洲av麻豆精品| 国产日韩欧美在线精品| 波多野结衣av一区二区av| 91成人精品电影| 亚洲国产精品国产精品| 亚洲专区中文字幕在线| 日韩精品免费视频一区二区三区| 成人亚洲欧美一区二区av| 制服诱惑二区| 午夜久久久在线观看| 亚洲欧美一区二区三区久久| 老司机靠b影院| 国产一卡二卡三卡精品| 亚洲av美国av| 在线天堂中文资源库| 国产在线视频一区二区| www.av在线官网国产| 色94色欧美一区二区| 国产av精品麻豆| 免费观看av网站的网址| 大话2 男鬼变身卡| 国产精品久久久久久人妻精品电影 | 日本a在线网址| 亚洲黑人精品在线| 日韩av在线免费看完整版不卡| 成人18禁高潮啪啪吃奶动态图| 在线 av 中文字幕| 涩涩av久久男人的天堂| 国产亚洲欧美在线一区二区| 国产黄频视频在线观看| 国产精品三级大全| 欧美性长视频在线观看| 搡老岳熟女国产| 两个人免费观看高清视频| 深夜精品福利| 日韩熟女老妇一区二区性免费视频| 人人妻人人澡人人爽人人夜夜| 久久女婷五月综合色啪小说| 免费女性裸体啪啪无遮挡网站| 天天操日日干夜夜撸| 欧美日韩成人在线一区二区| 婷婷色综合大香蕉| 免费不卡黄色视频| 日韩av免费高清视频| 丰满人妻熟妇乱又伦精品不卡| 久久综合国产亚洲精品| 国产av国产精品国产| 人妻人人澡人人爽人人| 久久久久网色| 秋霞在线观看毛片| 免费在线观看影片大全网站 | 国产av一区二区精品久久| 另类亚洲欧美激情| 成年人免费黄色播放视频| 国产熟女欧美一区二区| 亚洲国产欧美一区二区综合| 两性夫妻黄色片| 极品人妻少妇av视频| 亚洲精品国产一区二区精华液| 制服人妻中文乱码| 男男h啪啪无遮挡| 老司机影院成人| 日韩中文字幕欧美一区二区 | 亚洲黑人精品在线| 后天国语完整版免费观看| 人人妻人人澡人人看| 国产一区亚洲一区在线观看| 亚洲欧洲精品一区二区精品久久久| 久久久久久久精品精品| 美国免费a级毛片| 亚洲欧美日韩另类电影网站| 18禁观看日本| 青春草视频在线免费观看| 在线精品无人区一区二区三| 色网站视频免费| 色精品久久人妻99蜜桃| 免费在线观看黄色视频的| 亚洲成国产人片在线观看| 成人影院久久| 午夜91福利影院| 少妇 在线观看| 乱人伦中国视频| 赤兔流量卡办理| 国产精品一区二区在线观看99| 欧美在线一区亚洲| 精品国产乱码久久久久久男人| 欧美日韩亚洲综合一区二区三区_| 色精品久久人妻99蜜桃| 精品国产一区二区久久| 国产精品免费大片| 亚洲,一卡二卡三卡| 亚洲国产精品成人久久小说| 久久久久久久大尺度免费视频| 国产xxxxx性猛交| 看免费成人av毛片| 精品福利永久在线观看| 美女主播在线视频| 2018国产大陆天天弄谢| 一区二区三区四区激情视频| 美女国产高潮福利片在线看| 亚洲专区中文字幕在线| 少妇被粗大的猛进出69影院| 天天添夜夜摸| 色94色欧美一区二区| 久久久久久久久久久久大奶| 亚洲欧洲国产日韩| 一级毛片电影观看| 欧美乱码精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 伊人亚洲综合成人网| 美女福利国产在线| 国产又色又爽无遮挡免| 久久99一区二区三区| 黄色片一级片一级黄色片| 国产人伦9x9x在线观看| 老司机深夜福利视频在线观看 | 9191精品国产免费久久| 妹子高潮喷水视频| 在线观看一区二区三区激情| 美女高潮到喷水免费观看| 免费看十八禁软件| 777米奇影视久久| 精品高清国产在线一区| 久久久精品区二区三区| av欧美777| 黄频高清免费视频| bbb黄色大片| 黑人欧美特级aaaaaa片| www日本在线高清视频| 日本五十路高清| 一二三四社区在线视频社区8| 天天躁夜夜躁狠狠久久av| 女人精品久久久久毛片| 日日爽夜夜爽网站| 一区二区三区激情视频| 另类亚洲欧美激情| 男女下面插进去视频免费观看| 99热网站在线观看| 啦啦啦在线免费观看视频4| 一级黄片播放器| 久久亚洲精品不卡| 在线看a的网站| 久久久久久亚洲精品国产蜜桃av| 国产成人免费无遮挡视频| 侵犯人妻中文字幕一二三四区| 国产一区有黄有色的免费视频| 日韩制服骚丝袜av| 国产亚洲一区二区精品| 99精国产麻豆久久婷婷| 欧美另类一区| 亚洲精品国产av蜜桃| 肉色欧美久久久久久久蜜桃| 国产高清videossex| 晚上一个人看的免费电影| 久久99一区二区三区| 免费不卡黄色视频| 天天躁日日躁夜夜躁夜夜| 国产黄色视频一区二区在线观看| 在线亚洲精品国产二区图片欧美| 亚洲精品美女久久久久99蜜臀 | 亚洲欧美一区二区三区久久| 亚洲伊人色综图| 国产亚洲欧美在线一区二区| 免费在线观看黄色视频的| 久久 成人 亚洲| 91字幕亚洲| 国产成人精品在线电影| 亚洲国产av新网站| netflix在线观看网站| 精品人妻1区二区| 亚洲九九香蕉| 日韩大片免费观看网站| 高清欧美精品videossex| 国产精品麻豆人妻色哟哟久久| 亚洲精品国产av蜜桃| 老汉色av国产亚洲站长工具| 69精品国产乱码久久久| 国精品久久久久久国模美| 中文字幕亚洲精品专区| 成人免费观看视频高清| 精品一区二区三区四区五区乱码 | 国产精品久久久久久精品电影小说| 免费少妇av软件| 男的添女的下面高潮视频| 黑人猛操日本美女一级片| 久久人人爽人人片av| 黄片小视频在线播放| 国产色视频综合| 啦啦啦视频在线资源免费观看| 男人舔女人的私密视频| 色网站视频免费| 国产男人的电影天堂91| 欧美精品亚洲一区二区| 黄片小视频在线播放| 久9热在线精品视频| 国产成人免费无遮挡视频| 精品人妻1区二区| 最近最新中文字幕大全免费视频 | 久久人人97超碰香蕉20202| 99精国产麻豆久久婷婷| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品99久久99久久久不卡| 1024香蕉在线观看| 观看av在线不卡| 天天躁日日躁夜夜躁夜夜| 国产xxxxx性猛交| 久久精品成人免费网站| 久久青草综合色| 新久久久久国产一级毛片| 男女下面插进去视频免费观看| 成人影院久久| 亚洲专区中文字幕在线| 欧美激情高清一区二区三区| 国产精品三级大全| 啦啦啦啦在线视频资源| 人人妻人人澡人人爽人人夜夜| 国产一区亚洲一区在线观看| 一二三四在线观看免费中文在| 亚洲精品美女久久久久99蜜臀 | 一级黄色大片毛片| 欧美日本中文国产一区发布| 精品久久久久久电影网| 久久久久网色| 在线精品无人区一区二区三| 老司机午夜十八禁免费视频| 欧美日韩av久久| 国产一区二区三区综合在线观看| 日韩一卡2卡3卡4卡2021年| 美女国产高潮福利片在线看| 亚洲,欧美,日韩| 日韩一卡2卡3卡4卡2021年| 精品福利观看| 国产一区二区激情短视频 | 超色免费av| 久久精品国产综合久久久| 成人午夜精彩视频在线观看| 香蕉国产在线看| 亚洲国产最新在线播放| 国产欧美日韩一区二区三 | 国产亚洲av片在线观看秒播厂| 亚洲激情五月婷婷啪啪| 99国产精品99久久久久| 在线观看免费视频网站a站| 欧美日韩一级在线毛片| 男人爽女人下面视频在线观看| 亚洲久久久国产精品| 精品久久蜜臀av无| 超碰97精品在线观看| 男女之事视频高清在线观看 | 亚洲欧美清纯卡通| 午夜福利,免费看| 爱豆传媒免费全集在线观看| 黄网站色视频无遮挡免费观看| 亚洲九九香蕉| 50天的宝宝边吃奶边哭怎么回事| 国产精品麻豆人妻色哟哟久久| 久久精品亚洲av国产电影网| 超碰成人久久| www.av在线官网国产| 欧美另类一区| 国产亚洲av片在线观看秒播厂| 成人手机av| 久久精品亚洲熟妇少妇任你| 国产在线视频一区二区| 美女扒开内裤让男人捅视频| 伊人久久大香线蕉亚洲五| 又大又爽又粗| 日本欧美国产在线视频| av有码第一页| 久久影院123| 欧美av亚洲av综合av国产av| 69精品国产乱码久久久| 国产黄色视频一区二区在线观看| 久久国产精品大桥未久av| 人人澡人人妻人| 欧美黑人精品巨大| 久久久久久亚洲精品国产蜜桃av| 老汉色∧v一级毛片| 交换朋友夫妻互换小说| 中文字幕高清在线视频| 宅男免费午夜| 一边亲一边摸免费视频| 亚洲av成人不卡在线观看播放网 | 女人久久www免费人成看片| 国产片内射在线| 日韩 欧美 亚洲 中文字幕| 国产极品粉嫩免费观看在线| 国产精品一二三区在线看| 人成视频在线观看免费观看| 精品久久久精品久久久| 19禁男女啪啪无遮挡网站| 国产精品一区二区在线不卡| 美女中出高潮动态图| av片东京热男人的天堂| 国产欧美日韩一区二区三区在线| 亚洲精品久久久久久婷婷小说| www日本在线高清视频| 人体艺术视频欧美日本| 久久国产亚洲av麻豆专区| 人人妻人人爽人人添夜夜欢视频| 男女边吃奶边做爰视频| 国产精品久久久av美女十八| 少妇人妻久久综合中文| 免费在线观看影片大全网站 | 一区二区av电影网| 免费在线观看影片大全网站 | 操美女的视频在线观看| 嫩草影视91久久| 18禁黄网站禁片午夜丰满|