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

    旋轉(zhuǎn)加速度計(jì)重力梯度儀標(biāo)定方法

    2016-04-19 09:08:38錢(qián)學(xué)武蔡體菁
    關(guān)鍵詞:重力梯度零位標(biāo)度

    錢(qián)學(xué)武,蔡體菁

    (東南大學(xué) 儀器科學(xué)與工程學(xué)院,南京 210096)

    旋轉(zhuǎn)加速度計(jì)重力梯度儀標(biāo)定方法

    錢(qián)學(xué)武,蔡體菁

    (東南大學(xué) 儀器科學(xué)與工程學(xué)院,南京 210096)

    提出一種旋轉(zhuǎn)加速度計(jì)重力梯度儀標(biāo)定方法,在重力梯度儀外側(cè)空間確定的四個(gè)位置上依次放置一定質(zhì)量的檢測(cè)質(zhì)量體,根據(jù)檢測(cè)質(zhì)量體在正交位置上引起的重力梯度大小相等,符號(hào)相反的關(guān)系確定出重力梯度標(biāo)度系數(shù),采用旋轉(zhuǎn)重力梯度儀本體方式確定出重力梯度零位。給出了重力梯度零位和標(biāo)度系數(shù)計(jì)算表達(dá)式,在重力梯度半物理仿真系統(tǒng)上進(jìn)行了仿真試驗(yàn)驗(yàn)證,仿真結(jié)果表明,標(biāo)定后重力梯度測(cè)量誤差小于1E。

    重力梯度儀;旋轉(zhuǎn)加速度計(jì);檢測(cè)質(zhì)量體;梯度標(biāo)定

    高精度重力梯度測(cè)量在地球科學(xué)、空間科學(xué)、地質(zhì)科學(xué)、能源勘探以及慣性導(dǎo)航等領(lǐng)域具有非常重要的應(yīng)用價(jià)值和廣泛的應(yīng)用前景。重力梯度測(cè)量相對(duì)于重力測(cè)量具有明顯的優(yōu)勢(shì),對(duì)于地質(zhì)結(jié)構(gòu)的形狀、尺寸和埋藏深度等信息,重力梯度測(cè)量能夠提供更高的解析度,特別適合于探測(cè)局部地質(zhì)結(jié)構(gòu)及其細(xì)節(jié)[1-4]。目前,國(guó)外正在進(jìn)行研究具有應(yīng)用前景的重力梯度儀主要有旋轉(zhuǎn)加速度計(jì)重力梯度儀、超導(dǎo)重力梯度儀、冷原子重力梯度儀等,其中旋轉(zhuǎn)加速度計(jì)重力梯度儀技術(shù)最成熟,也是目前唯一成功用于機(jī)載/船載動(dòng)機(jī)座并進(jìn)行商用的重力梯度測(cè)量系統(tǒng),已經(jīng)進(jìn)行了大量的勘探工作,并取得了良好的效果[5-8]。20世紀(jì)90年代國(guó)內(nèi)開(kāi)始研究重力梯度儀[9],目前有的科研單位已經(jīng)研制出了重力梯度儀實(shí)驗(yàn)室樣機(jī)[10-12]。

    旋轉(zhuǎn)加速度計(jì)重力梯度儀是一種高精密測(cè)量?jī)x器,在出廠(chǎng)前或使用前需要對(duì)其進(jìn)行標(biāo)定。目前,關(guān)于旋轉(zhuǎn)加速度計(jì)重力梯度儀的標(biāo)定方法在公開(kāi)刊物上沒(méi)有報(bào)道。

    重力梯度儀標(biāo)定一般是在實(shí)驗(yàn)室環(huán)境條件下進(jìn)行,在重力梯度儀外側(cè)放置形狀規(guī)則的均質(zhì)檢測(cè)質(zhì)量體,通過(guò)改變檢測(cè)質(zhì)量體與GGI圓盤(pán)中心的距離來(lái)改變GGI中心處的引力場(chǎng),變化的引力場(chǎng)被重力梯度敏感器敏感到,然后對(duì)重力梯度儀的輸出信號(hào)進(jìn)行重力梯度解調(diào),就可以獲得距離與GGI測(cè)量梯度之間的對(duì)應(yīng)關(guān)系,根據(jù)這種對(duì)應(yīng)關(guān)系與重力梯度理論值進(jìn)行對(duì)比,可以確定出重力梯度零位和標(biāo)度系數(shù),但由于引力梯度對(duì)距離敏感,檢測(cè)質(zhì)量體距GGI較近時(shí)會(huì)產(chǎn)生較大畸變,從而導(dǎo)致較大標(biāo)定誤差。為此,本文提出一種旋轉(zhuǎn)加速度計(jì)重力梯度儀部分張量的高精度標(biāo)定方法。

    1 重力梯度標(biāo)定方法

    重力梯度儀主要有全張量和部分張量?jī)煞N類(lèi)型,針對(duì)全張量重力梯度儀需要對(duì)五個(gè)獨(dú)立重力梯度元素進(jìn)行標(biāo)定,全張量重力梯度儀標(biāo)定方法與部分張量重力梯度儀標(biāo)定方法類(lèi)似,本文只對(duì)部分張量重力梯度儀標(biāo)定方法進(jìn)行分析。對(duì)于部分張量重力梯度儀,需要確定重力梯度分量(Гyy-Гxx)和 Гxy的零位和標(biāo)度系數(shù),重力梯度儀誤差模型方程為

    式中:Γc為重力梯度儀實(shí)際測(cè)量值;K為重力梯度標(biāo)度系數(shù);Γt為檢測(cè)質(zhì)量體引起的理論重力梯度,即重力梯度精確值;G0為周?chē)h(huán)境物體引起的重力梯度,Γ0為重力梯度儀儀器本身性能引起的重力梯度值,即為重力梯度零位,重力梯度零位主要由GGI動(dòng)力學(xué)擾動(dòng)、線(xiàn)路噪聲竄擾等因素引起,不考慮時(shí)變影響,它是一個(gè)定值。從式(1)可以得到,G0和 Γ0相互耦合,不易提取出重力梯度零位,而G0主要與環(huán)境物體相對(duì)于GGI的位置有關(guān),只要環(huán)境物體相對(duì)位置不改變,其值也不會(huì)改變。旋轉(zhuǎn)重力梯度儀本體相當(dāng)于環(huán)境物體相對(duì)于GGI位置的改變,若重力梯度儀本體旋轉(zhuǎn)整數(shù)倍周期,那么G0將以周期性呈現(xiàn),利用G0的周期性和Γ0的定值性,可以確定出重力梯度零位Γ0。

    以GGI圓盤(pán)中心為原點(diǎn)O,以其所在平面為xy平面建立圓盤(pán)坐標(biāo)系。假設(shè)檢測(cè)質(zhì)量體為均質(zhì)長(zhǎng)方體,其密度為ρ,檢測(cè)質(zhì)量體質(zhì)心坐標(biāo)為Q(W, D, H),檢測(cè)質(zhì)量體內(nèi)任意一點(diǎn)坐標(biāo)為P(x, y, z),檢測(cè)質(zhì)量體的長(zhǎng)、寬、高分別設(shè)為w、d、h,如圖1所示。

    圖1 檢測(cè)質(zhì)量體作用于重力梯度儀示意圖Fig.1 Schematic for proof mass’s acting on GGI

    根據(jù)萬(wàn)有引力定律,可得檢測(cè)質(zhì)量體對(duì)GGI圓盤(pán)中心的重力梯度分量計(jì)算表達(dá)式為:

    式中:G為萬(wàn)有引力常數(shù)。從式(2)可以得到重力梯度與檢測(cè)質(zhì)量體的關(guān)系。為了便于說(shuō)明檢測(cè)質(zhì)量體位置對(duì)重力梯度的影響,對(duì)式(2)變換為:

    式中:θ為GGI圓盤(pán)中心到檢測(cè)質(zhì)量體質(zhì)心矢量在圓盤(pán)平面上的投影與圓盤(pán)坐標(biāo)系x軸正方向間的夾角,其中,

    從式(3)可以看出,檢測(cè)質(zhì)量體相對(duì)于重力梯度儀旋轉(zhuǎn)θ角,相當(dāng)于重力梯度儀相對(duì)于檢測(cè)質(zhì)量體以反方向旋轉(zhuǎn)θ角。設(shè)重力梯度儀本體旋轉(zhuǎn)一周,旋轉(zhuǎn)角度間隔為 10°,在每次旋轉(zhuǎn)位置上測(cè)量一次重力梯度值,仿真結(jié)果如圖2所示,從圖2可以看出,通過(guò)旋轉(zhuǎn)重力梯度儀本體,環(huán)境物體引起的重力梯度G0以周期性呈現(xiàn),而重力梯度零位Γ0保持常值,采用這種方法可以分離出環(huán)境梯度 G0和重力梯度零位 Γ0,從而可得重力梯度零位Γ0的計(jì)算表達(dá)式為:

    從式(3)可以進(jìn)一步得到:當(dāng)θ=0°, 90°, 180°, 360°時(shí),即當(dāng)檢測(cè)質(zhì)量體質(zhì)心在x或y坐標(biāo)軸上時(shí),檢測(cè)質(zhì)量體引起的重力梯度分量 Гxy=0E;當(dāng)θ=45°, 135°, 225°, 315°時(shí),即當(dāng)檢測(cè)質(zhì)量體質(zhì)心在y=±x直線(xiàn)方程上時(shí),檢測(cè)質(zhì)量體引起的重力梯度分量(Гyy-Гxx)=0E;當(dāng)檢測(cè)質(zhì)量體質(zhì)心在θ=0°, 90°位置上時(shí),重力梯度分量(Гyy-Гxx)達(dá)到最大值,且兩位置處的重力梯度分量(Гyy-Гxx)互為相反數(shù),因此只要對(duì)兩位置處的重力梯度分量(Гyy-Гxx)相減然后取均值,可以消除環(huán)境梯度以及由儀器本身引起的重力梯度零位。

    圖2 環(huán)境物體引起的重力梯度Fig.2 Gravity gradient caused by ambient object

    同理,檢測(cè)質(zhì)量體質(zhì)心在θ=45°, 135°位置上時(shí),重力梯度分量 Гxy達(dá)到最大值,且兩位置處的重力梯度分量Гxy互為相反數(shù)。

    基于以上分析,給出重力梯度分量(Гyy-Гxx)和Гxy的標(biāo)定方式:當(dāng)對(duì)重力梯度分量(Гyy-Гxx)進(jìn)行標(biāo)定時(shí),設(shè)定檢測(cè)質(zhì)量體質(zhì)心在x和y坐標(biāo)軸上;當(dāng)對(duì)重力梯度分量Гxy進(jìn)行標(biāo)定時(shí),設(shè)定質(zhì)量體質(zhì)心在y=±x直線(xiàn)關(guān)系上,且4個(gè)位置上的檢測(cè)質(zhì)量體質(zhì)心到GGI圓盤(pán)中心的距離相等,檢測(cè)質(zhì)量體放置位置如圖3所示。

    圖3 檢測(cè)質(zhì)量體配置示意圖Fig.3 Setting of proof mass bodies

    檢測(cè)質(zhì)量體以圖3所示位置配置方式進(jìn)行放置,根據(jù)重力梯度分量計(jì)算式(3)可知,檢測(cè)質(zhì)量體在位置1與位置2處引起的重力梯度分量(Гyy-Гxx)大小相等,符號(hào)相反;同理,檢測(cè)質(zhì)量體在位置3與位置4處引起的重力梯度分量 Гxy大小相等,符號(hào)相反,根據(jù)上述關(guān)系給出4個(gè)位置檢測(cè)質(zhì)量體引起的重力梯度分量測(cè)量表達(dá)式為

    式中:Kyy-xx和Kxy分別表示重力梯度分量(Гyy-Гxx)和Гxy的標(biāo)度系數(shù),Gyy-xx和Gxy分別表示由環(huán)境物體引起的重力梯度分量分別表示檢測(cè)質(zhì)量體在位置1和位置2實(shí)際引起的理論重力梯度分量分別表示檢測(cè)質(zhì)量體在位置3和位置4實(shí)際引起的理論重力梯度分量Гxy,分別表示重力梯度分量零位。根據(jù)重力梯度零位計(jì)算表達(dá)式(5)可以得到重力梯度分量(Гyy-Гxx)和Гxy的零位計(jì)算表達(dá)式為

    位置1與位置2處的重力梯度分量(Гyy-Гxx)測(cè)量值相減,消除環(huán)境物體和零位引起的重力梯度分量(Гyy-Гxx),得到如下關(guān)系式:

    同理,位置3與位置4處的重力梯度分量Гxy測(cè)量值相減,消除環(huán)境物體和零位引起的重力梯度分量Гxy,得到重力梯度分量Гxy的標(biāo)度系數(shù)計(jì)算表達(dá)式為

    在進(jìn)行重力梯度測(cè)量時(shí),重力梯度儀真正測(cè)量的重力梯度信息是GGI周?chē)h(huán)境物體引起的重力梯度異常,因此重力梯度分量計(jì)算表達(dá)式為

    1)在沒(méi)有檢測(cè)質(zhì)量體下,旋轉(zhuǎn)重力梯度儀本體整數(shù)倍圈數(shù),在相同旋轉(zhuǎn)間隔位置處進(jìn)行重力梯度解調(diào),利用式(7)計(jì)算得到重力梯度零位Γ0;

    2)在正交位置處依次放置檢測(cè)質(zhì)量體,然后對(duì)重力梯度儀輸出信號(hào)進(jìn)行重力梯度解調(diào),利用式(9)和式(10)計(jì)算重力梯度標(biāo)度系數(shù)K;

    3)采用式(11)計(jì)算出重力梯度儀真正測(cè)量到的重力梯度分量(Гyy-Гxx)和Гxy。

    2 半物理仿真實(shí)驗(yàn)與結(jié)果分析

    采用上述給出的重力梯度分量零位和標(biāo)度系數(shù)計(jì)算方法,在重力梯度儀半物理仿真分析系統(tǒng)上進(jìn)行重力梯度標(biāo)定仿真實(shí)驗(yàn),仿真參數(shù)如下:加速度計(jì)標(biāo)度因數(shù)為10 mA/g,加速度計(jì)對(duì)基線(xiàn)距離為0.2 m,檢測(cè)質(zhì)量體邊長(zhǎng)為0.3 m的均質(zhì)正方體,密度為18 000 kg/m3,檢測(cè)質(zhì)量體質(zhì)心到GGI圓盤(pán)中心的距離為1 m。根據(jù)重力梯度分量計(jì)算式(2)可以計(jì)算出檢測(cè)質(zhì)量體在位置 1處引起的重力梯度分量(Гyy-Гxx)理論值為-97.53E,在位置3處引起的重力梯度分量Гxy理論值為48.76E。為了模擬環(huán)境物體對(duì)重力梯度儀的影響,設(shè)環(huán)境物體邊長(zhǎng)為 1 m的均質(zhì)正方體,密度為 2000 kg/m3,環(huán)境物體質(zhì)心位置坐標(biāo)為(2, 1, 0) m,此環(huán)境物體對(duì)GGI引起的重力梯度分量(Гyy-Гxx)和Гxy理論值分別為-21.63E和14.26E。加入2倍頻干擾信號(hào)用于模擬重力梯度半物理仿真系統(tǒng)性能引起的重力梯度零位,重力梯度分量(Гyy-Гxx)和 Гxy的零位分別為 200E和10E;檢測(cè)質(zhì)量體依次被放置在圖3所示的位置1、2、3、4上,然后分別對(duì)重力梯度儀仿真系統(tǒng)輸出信號(hào)進(jìn)行數(shù)據(jù)采集和重力梯度解調(diào)。設(shè)重力梯度儀本體旋轉(zhuǎn)一周,旋轉(zhuǎn)角度間隔為 20°,由環(huán)境物體引起的重力梯度分量輸出信號(hào)如圖4所示,從圖4可以看出,重力梯度以近似正弦波形式呈現(xiàn),這與理論分析結(jié)果一致。采用式(7)計(jì)算出半物理仿真系統(tǒng)重力梯度分量(Гyy-Гxx)和 Гxy的零位分別為190.93E和9.98E。檢測(cè)質(zhì)量體在不同位置處產(chǎn)生的重力梯度如表1所示,表1中上標(biāo)t為重力梯度理論值,c為實(shí)際測(cè)量重力梯度值。根據(jù)表1中的重力梯度分量數(shù)據(jù),采用式(9)和式(10)計(jì)算出重力梯度分量(Гyy-Гxx)和 Гxy的標(biāo)度系數(shù)分別為0.9578和0.9565。重力梯度零位和標(biāo)度系數(shù)是重力梯度測(cè)量系統(tǒng)固有特性,它只與重力梯度儀本身性能有關(guān),與外界因素?zé)o關(guān)。

    為了檢驗(yàn)重力梯度零位和標(biāo)度系數(shù)的正確性,假設(shè)檢驗(yàn)質(zhì)量體參數(shù)為:質(zhì)量體為均質(zhì)長(zhǎng)方體,長(zhǎng)、寬和高分別為0.2 m、0.5 m、0.4 m,密度為10 000 kg/m3,其質(zhì)心坐標(biāo)為(0.8, 0.4, 0) m,此檢驗(yàn)質(zhì)量體對(duì)GGI引起的重力梯度分量(Гyy-Гxx)和 Гxy理論值分別為-68.63E和 39.83E,標(biāo)定前后重力梯度測(cè)量值如表 2所示,從表2的結(jié)果可以得到,重力梯度標(biāo)定前,由于儀器性能等原因引起的重力梯度被引入到重力梯度信息中,從而造成重力梯度測(cè)量誤差,通過(guò)對(duì)重力梯度標(biāo)定,標(biāo)定后重力梯度測(cè)量值與重量梯度理論值之間的誤差小于1E,因此提高了重力梯度測(cè)量精度。

    圖4 半物理仿真系統(tǒng)環(huán)境梯度仿真結(jié)果Fig.4 Gravity gradient caused by ambient object on the Hardware-in-the-loop simulation platform

    表1 不同位置處的檢測(cè)質(zhì)量體產(chǎn)生的重力梯度Tab.1 Gravity gradient caused by proof mass in different positions

    表2 標(biāo)定前后重力梯度測(cè)量值Tab.2 Gravity gradiometry of GGI before and after calibration

    3 結(jié) 論

    本文從系統(tǒng)級(jí)角度提出了一種部分張量旋轉(zhuǎn)加速度計(jì)重力梯度儀的標(biāo)定方法,給出了重力梯度零位和標(biāo)度系數(shù)的計(jì)算表達(dá)式以及標(biāo)定步驟。通過(guò)重力梯度半物理仿真系統(tǒng)試驗(yàn),仿真結(jié)果表明該方法能夠?qū)崿F(xiàn)重力梯度儀的標(biāo)定。

    (References):

    [1] Dransfield M H, Christensen A N. Performance of airborne gravity gradiometers[J]. The Leading Edge, 2013, 32(8): 908-922.

    [2] DiFrancesco D, Grierson A, Kaputa D, et al. Gravity gradiometer systems - advances and challenges[J]. Geophysical Prospecting, 2009, 57(4): 615-623.

    [3] Welker T C, Pachter M, Huffman R E. Gravity gradiometer integrated inertial navigation[C]//2013 European Control Conference. Zurich, Switzerland, 2013: 846-851.

    [4] Hodges G, Dransfield M H, Shei T C. The falcon airborne gravity gradiometer for engineering applications[C]//23rdSymposium on the Application of Geophysics to Engineering and Environmental Problems 2010. Keystone, USA, 2010: 443-447.

    [5] Uzun S. Estimating parameters of subsurface structures from airborne gravity gradiometry data using a Monte-Carlo optimization method[R]. Geodetic Science, The Ohio State University, Columbus, Ohio, USA, 2013.

    [6] Rummel R, Yi W, Stummer C. GOCE gravitational gradiometry[J]. Journal of Geodesy, 2011, 85(11): 777-790.

    [7] Roberts D, Chowdhury P R, Lowe S J, et al. Airborne gravity gradiometer surveying of petroleum systems under Lake Tanganyika, Tanzania[C]//ASEG-PESA 2015. Perth, Australia, 2015: 1-5.

    [8] Christensen A N, Galder C V, Dransfield M. Improved resolution of fixed-wing airborne gravity gradiometer surveys[C]//2014 SEG Annual Meeting. Denver, Colorado, USA, 2014: 1319-1323.

    [9] 蔡體菁, 周百令. 重力梯度儀的現(xiàn)狀和前景[J]. 中國(guó)慣性技術(shù)學(xué)報(bào), 1999, 7(1): 39-42. Cai Ti-jing, Zhou Bai-ling. Status and prospects of gravity gradiometers[J]. Journal of Chinese Inertial Technology, 1999, 7(1): 39-42.

    [10] 劉昱, 張宇, 潘軍, 等. 重力梯度測(cè)量中光柵角編碼器安裝誤差的影響及分析[J]. 中國(guó)慣性技術(shù)學(xué)報(bào), 2015, 23(1): 6-8. Liu Yu, Zhang Yu, Pan Jun, et al. Influence and analysis of optical angular encoder’s installation error in gravity gradiometry[J]. Journal of Chinese Inertial Technology, 2015, 23(1): 6-8.

    [11] 李海兵, 楊慧, 馬存尊, 等. 用于重力梯度儀系統(tǒng)中的高精度微弱信號(hào)檢測(cè)電路[J]. 中國(guó)慣性技術(shù)學(xué)報(bào), 2013, 21(5): 581-584. Li Hai-bing, Yang Hui, Ma Cun-zun, et al. High precision and week signal detection circuits of gravity gradiometer system[J]. Journal of Chinese Inertial Technology, 2013, 21(5): 581-584.

    [12] Tu L C, Wang Z W, Liu J Q, et al. Implementation of the scale factor balance on two pairs of quartz-flexure capacitive accelerometers by trimming bias voltage[J]. Review of Scientific Instruments, 2014, 85: 095108-1-7.

    Calibration method for gravity gradiometer of rotating accelerometer

    QIAN Xue-wu, CAI Ti-jing
    (School of Instrument Science and Engineering, Southeast University, Nanjing 210096, China)

    A calibration method for rotating accelerometer gravity gradiometer instrument (GGI) is presented. The proof mass bodies are sequentially placed at four positions outside the GGI. The gravity gradient scale coefficient is determined based on the relationship at perpendicular position of gravity gradient component, i.e. the gradients caused by proof masses at orthogonal positions are the same in magnitude and opposite in sign. The gravity gradient bias is determined by the method of rotating the GGI body. The calculation expression of gravity gradient bias and scale coefficient are deduced in details. A validation test of calibration method is performed on the Hardware-in-the-loop simulation platform. Simulation results show that the measurement error of the calibrated gravity gradient is less than 1E.

    gravity gradiometer instrument; rotating accelerometer; proof mass; gradient calibration

    U666.1

    :A

    2016-04-27;

    :2016-07-08

    國(guó)家863計(jì)劃項(xiàng)目(2011AA060501)

    錢(qián)學(xué)武(1981—),男,博士研究生,從事精密儀器研究。E-mail: njqxwu@163.com

    聯(lián) 系 人:蔡體菁(1961—),男,教授,博士生導(dǎo)師。E-mail: caitij@seu.edu.cn

    1005-6734(2016)04-0427-04

    10.13695/j.cnki.12-1222/o3.2016.04.002

    猜你喜歡
    重力梯度零位標(biāo)度
    層次分析法中兩種標(biāo)度的對(duì)比分析
    航空重力梯度儀實(shí)時(shí)重力梯度解調(diào)方法
    物探與化探(2022年6期)2023-01-03 12:14:06
    無(wú)刷力矩電機(jī)零位設(shè)計(jì)技術(shù)研究
    電機(jī)零位對(duì)電動(dòng)汽車(chē)能耗影響研究
    利用地形數(shù)據(jù)計(jì)算重力梯度張量的直接積分法
    加權(quán)無(wú)標(biāo)度網(wǎng)絡(luò)上SIRS 類(lèi)傳播模型研究
    星載重力梯度儀的研究發(fā)展
    霍爾傳感器零位特性及補(bǔ)償方法研究
    分析起重機(jī)零位保護(hù)及檢驗(yàn)
    創(chuàng)新孵化網(wǎng)絡(luò)演化無(wú)標(biāo)度特征仿真分析
    久久久久久国产a免费观看| 狠狠狠狠99中文字幕| 国产精品,欧美在线| 中文字幕av在线有码专区| 精品熟女少妇av免费看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲第一电影网av| 狂野欧美激情性xxxx在线观看| 精品久久久久久久久久免费视频| 色播亚洲综合网| 一区二区三区四区激情视频 | 欧美人与善性xxx| 一本精品99久久精品77| 一卡2卡三卡四卡精品乱码亚洲| 麻豆av噜噜一区二区三区| 亚洲精品久久国产高清桃花| 一区二区三区四区激情视频 | 久久久精品94久久精品| 久久人人爽人人爽人人片va| 美女大奶头视频| 日产精品乱码卡一卡2卡三| 亚洲国产日韩欧美精品在线观看| 美女cb高潮喷水在线观看| 免费大片18禁| 精品乱码久久久久久99久播| 村上凉子中文字幕在线| 两个人的视频大全免费| 日日摸夜夜添夜夜添小说| 亚洲久久久久久中文字幕| 插阴视频在线观看视频| 少妇裸体淫交视频免费看高清| 国产私拍福利视频在线观看| 麻豆久久精品国产亚洲av| 成人永久免费在线观看视频| 直男gayav资源| 国产精品一二三区在线看| 毛片女人毛片| 激情 狠狠 欧美| 国产私拍福利视频在线观看| av免费在线看不卡| 亚洲精品在线观看二区| 日本黄大片高清| 色综合色国产| 看片在线看免费视频| av黄色大香蕉| 亚洲一区高清亚洲精品| 午夜精品在线福利| 噜噜噜噜噜久久久久久91| 成人三级黄色视频| 又粗又爽又猛毛片免费看| 亚洲成人av在线免费| 熟女电影av网| 久久久久国产精品人妻aⅴ院| 亚洲人成网站在线观看播放| 欧美高清性xxxxhd video| 久久精品国产99精品国产亚洲性色| 亚洲欧美日韩卡通动漫| 色哟哟·www| 亚洲自偷自拍三级| 91麻豆精品激情在线观看国产| 亚洲成人中文字幕在线播放| 自拍偷自拍亚洲精品老妇| 国内精品美女久久久久久| 两个人的视频大全免费| 看十八女毛片水多多多| 日韩在线高清观看一区二区三区| 国产在线男女| 久久人妻av系列| 国产淫片久久久久久久久| 久久综合国产亚洲精品| 国产一区亚洲一区在线观看| 嫩草影视91久久| 一级黄色大片毛片| 亚洲三级黄色毛片| 九九爱精品视频在线观看| 久99久视频精品免费| 日日摸夜夜添夜夜爱| 最近视频中文字幕2019在线8| 午夜a级毛片| 欧美日韩一区二区视频在线观看视频在线 | 天美传媒精品一区二区| 免费高清视频大片| 久久精品国产亚洲av香蕉五月| 插逼视频在线观看| 欧美激情国产日韩精品一区| 成人精品一区二区免费| 国产精品野战在线观看| 欧美高清性xxxxhd video| 别揉我奶头 嗯啊视频| 九九爱精品视频在线观看| 国内精品久久久久精免费| 精品一区二区三区av网在线观看| 天堂网av新在线| 精品人妻偷拍中文字幕| 亚洲中文日韩欧美视频| 久久精品国产鲁丝片午夜精品| 禁无遮挡网站| 一本久久中文字幕| 免费人成视频x8x8入口观看| 国产精品女同一区二区软件| 亚洲国产欧美人成| 变态另类丝袜制服| 国产精品人妻久久久影院| 最新中文字幕久久久久| 3wmmmm亚洲av在线观看| 91av网一区二区| 亚洲精品成人久久久久久| 日本熟妇午夜| 在线播放无遮挡| 99精品在免费线老司机午夜| 99久久精品一区二区三区| 直男gayav资源| 国产v大片淫在线免费观看| 久久精品国产亚洲av天美| 亚洲一区高清亚洲精品| 久久久久九九精品影院| 日韩亚洲欧美综合| 亚洲av免费高清在线观看| 亚洲欧美日韩高清专用| 麻豆成人午夜福利视频| 亚洲四区av| aaaaa片日本免费| 亚洲一区高清亚洲精品| 久久精品国产99精品国产亚洲性色| 国产中年淑女户外野战色| 成熟少妇高潮喷水视频| 国产精品久久电影中文字幕| 免费在线观看成人毛片| 成人二区视频| 欧美极品一区二区三区四区| 别揉我奶头 嗯啊视频| 精品一区二区三区av网在线观看| 久久精品国产亚洲av涩爱 | 亚洲av五月六月丁香网| 免费观看人在逋| 国产免费男女视频| 波多野结衣高清作品| 欧美在线一区亚洲| 国内揄拍国产精品人妻在线| 国产高清激情床上av| www.色视频.com| 亚洲精品亚洲一区二区| av在线观看视频网站免费| 亚洲乱码一区二区免费版| 亚洲精品日韩在线中文字幕 | 99视频精品全部免费 在线| h日本视频在线播放| 网址你懂的国产日韩在线| 国产午夜精品论理片| 精品久久久久久久人妻蜜臀av| 亚洲av中文字字幕乱码综合| 成人永久免费在线观看视频| 国产精品亚洲美女久久久| 久久精品国产亚洲av涩爱 | 亚洲天堂国产精品一区在线| 在现免费观看毛片| 最近在线观看免费完整版| av女优亚洲男人天堂| 乱系列少妇在线播放| 秋霞在线观看毛片| 干丝袜人妻中文字幕| 亚洲精品456在线播放app| 亚洲欧美精品综合久久99| 中出人妻视频一区二区| 内地一区二区视频在线| 内地一区二区视频在线| 内地一区二区视频在线| 久久精品国产鲁丝片午夜精品| 精品日产1卡2卡| 亚洲av不卡在线观看| 亚洲欧美精品综合久久99| 日日摸夜夜添夜夜添小说| 啦啦啦韩国在线观看视频| 美女被艹到高潮喷水动态| av专区在线播放| 91在线观看av| 中国美白少妇内射xxxbb| 天天一区二区日本电影三级| 看免费成人av毛片| 你懂的网址亚洲精品在线观看 | 69人妻影院| 男女下面进入的视频免费午夜| 99热网站在线观看| 久久精品国产亚洲av香蕉五月| 看十八女毛片水多多多| 熟女电影av网| 国产高潮美女av| 国产亚洲精品久久久久久毛片| 午夜福利在线观看免费完整高清在 | 中文字幕熟女人妻在线| 一区二区三区四区激情视频 | eeuss影院久久| 搞女人的毛片| 国产老妇女一区| 22中文网久久字幕| 成人亚洲欧美一区二区av| 久久久久精品国产欧美久久久| 欧美潮喷喷水| 国产伦在线观看视频一区| 国产精品一区二区三区四区免费观看 | 国产精品国产高清国产av| 国产视频一区二区在线看| 村上凉子中文字幕在线| 久久久成人免费电影| 色综合站精品国产| 美女内射精品一级片tv| 国内精品宾馆在线| 男人狂女人下面高潮的视频| 久久久a久久爽久久v久久| a级一级毛片免费在线观看| 亚洲,欧美,日韩| 日本-黄色视频高清免费观看| 神马国产精品三级电影在线观看| 国产精品久久久久久精品电影| 亚洲无线在线观看| 国产精品久久久久久av不卡| 三级毛片av免费| 午夜久久久久精精品| 免费av观看视频| 成人特级黄色片久久久久久久| 99国产极品粉嫩在线观看| 国产精品久久久久久精品电影| 国产精品精品国产色婷婷| 男插女下体视频免费在线播放| 欧美3d第一页| 国产成人a区在线观看| 亚洲国产精品合色在线| 国内精品美女久久久久久| 黄片wwwwww| 国产精品日韩av在线免费观看| 免费搜索国产男女视频| 少妇熟女欧美另类| 亚洲久久久久久中文字幕| 成人精品一区二区免费| 国产视频一区二区在线看| 精华霜和精华液先用哪个| 精品久久国产蜜桃| 久久午夜福利片| 非洲黑人性xxxx精品又粗又长| 别揉我奶头 嗯啊视频| 久久亚洲国产成人精品v| 三级经典国产精品| 婷婷精品国产亚洲av| av中文乱码字幕在线| 国产精品,欧美在线| 在线观看美女被高潮喷水网站| 国模一区二区三区四区视频| 久久精品国产亚洲网站| 中文字幕人妻熟人妻熟丝袜美| 欧美激情国产日韩精品一区| 老熟妇仑乱视频hdxx| 久久99热这里只有精品18| 欧美性感艳星| 可以在线观看的亚洲视频| 欧洲精品卡2卡3卡4卡5卡区| 99久久精品一区二区三区| 久久精品91蜜桃| 深夜a级毛片| 99热只有精品国产| av在线天堂中文字幕| 嫩草影院入口| 在线免费观看不下载黄p国产| 精品久久久久久久久久久久久| 中文字幕av在线有码专区| 亚洲色图av天堂| 久久久久久九九精品二区国产| 成人国产麻豆网| 少妇的逼好多水| 午夜免费激情av| 亚洲av电影不卡..在线观看| 欧美高清成人免费视频www| 寂寞人妻少妇视频99o| 日产精品乱码卡一卡2卡三| 国产亚洲91精品色在线| 狠狠狠狠99中文字幕| 嫩草影院入口| 国产成人aa在线观看| 久久精品国产鲁丝片午夜精品| 天天一区二区日本电影三级| av天堂在线播放| 亚洲自偷自拍三级| 最近中文字幕高清免费大全6| 国产精品美女特级片免费视频播放器| 国产精品永久免费网站| 又黄又爽又免费观看的视频| 神马国产精品三级电影在线观看| 天堂√8在线中文| 人人妻人人澡欧美一区二区| 精品一区二区免费观看| 日日摸夜夜添夜夜爱| 日本撒尿小便嘘嘘汇集6| 天堂av国产一区二区熟女人妻| 成年女人看的毛片在线观看| 国内揄拍国产精品人妻在线| 在线观看美女被高潮喷水网站| 欧美色视频一区免费| 色av中文字幕| 日韩,欧美,国产一区二区三区 | 最近视频中文字幕2019在线8| 能在线免费观看的黄片| 欧美一级a爱片免费观看看| 中文亚洲av片在线观看爽| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av免费高清在线观看| 真实男女啪啪啪动态图| 熟女电影av网| 精品国产三级普通话版| 日本爱情动作片www.在线观看 | 免费av不卡在线播放| 亚洲精品一区av在线观看| 久久久久久久久大av| 一区二区三区四区激情视频 | 日日摸夜夜添夜夜爱| 亚洲一区高清亚洲精品| a级一级毛片免费在线观看| 人妻夜夜爽99麻豆av| 日日撸夜夜添| 亚洲五月天丁香| 国产精品久久电影中文字幕| 特级一级黄色大片| 久久欧美精品欧美久久欧美| 国产一区二区激情短视频| 麻豆乱淫一区二区| 在线国产一区二区在线| 日韩欧美 国产精品| 少妇的逼好多水| 最近手机中文字幕大全| 乱码一卡2卡4卡精品| 亚洲专区国产一区二区| 嫩草影视91久久| 国产一区二区亚洲精品在线观看| 成人av在线播放网站| 人妻丰满熟妇av一区二区三区| 三级男女做爰猛烈吃奶摸视频| 草草在线视频免费看| 日韩av在线大香蕉| 老司机午夜福利在线观看视频| 久久精品国产亚洲av涩爱 | 麻豆久久精品国产亚洲av| 一进一出好大好爽视频| 18禁在线无遮挡免费观看视频 | 久久热精品热| 亚洲国产欧美人成| 国产人妻一区二区三区在| 中文字幕人妻熟人妻熟丝袜美| 女同久久另类99精品国产91| 国产伦精品一区二区三区视频9| 少妇人妻精品综合一区二区 | 插逼视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 欧美性感艳星| 在线免费观看不下载黄p国产| 黄色欧美视频在线观看| 日本黄大片高清| 亚洲av成人av| 久久久色成人| av在线老鸭窝| 在线观看免费视频日本深夜| 日本 av在线| 最后的刺客免费高清国语| 国产大屁股一区二区在线视频| 老熟妇乱子伦视频在线观看| 亚洲熟妇中文字幕五十中出| 99在线视频只有这里精品首页| 欧美潮喷喷水| 又粗又爽又猛毛片免费看| 亚洲av第一区精品v没综合| 免费看日本二区| 大香蕉久久网| 日韩制服骚丝袜av| 一本久久中文字幕| 黄色欧美视频在线观看| 久久热精品热| 免费看美女性在线毛片视频| 亚洲欧美精品综合久久99| 少妇高潮的动态图| avwww免费| 亚洲国产精品成人综合色| 欧美高清性xxxxhd video| 久久精品91蜜桃| 亚洲精品国产成人久久av| 狠狠狠狠99中文字幕| 最近的中文字幕免费完整| 国产精品不卡视频一区二区| eeuss影院久久| 天堂动漫精品| 一个人观看的视频www高清免费观看| 欧美色欧美亚洲另类二区| 久久人妻av系列| 免费在线观看影片大全网站| 国产一区二区激情短视频| 99热6这里只有精品| 联通29元200g的流量卡| 国产精品综合久久久久久久免费| 国产精品爽爽va在线观看网站| 日本黄色片子视频| 成年版毛片免费区| 国产蜜桃级精品一区二区三区| 国内精品宾馆在线| 联通29元200g的流量卡| 日本精品一区二区三区蜜桃| 少妇被粗大猛烈的视频| 精品免费久久久久久久清纯| 99热这里只有精品一区| 国内精品一区二区在线观看| 亚洲熟妇中文字幕五十中出| 99热全是精品| 免费看光身美女| 草草在线视频免费看| 亚洲最大成人手机在线| 免费一级毛片在线播放高清视频| 免费看a级黄色片| 日日摸夜夜添夜夜爱| 亚洲国产精品国产精品| 免费观看在线日韩| a级毛片免费高清观看在线播放| 日本黄色片子视频| 大香蕉久久网| 长腿黑丝高跟| 国产亚洲欧美98| 免费无遮挡裸体视频| 亚洲av.av天堂| 午夜激情福利司机影院| 国产成人freesex在线 | 国产精品嫩草影院av在线观看| 国产亚洲精品av在线| 婷婷色综合大香蕉| 午夜老司机福利剧场| 99久久精品一区二区三区| 久久国产乱子免费精品| 免费人成视频x8x8入口观看| 久久人妻av系列| 此物有八面人人有两片| 69av精品久久久久久| 观看免费一级毛片| 久久热精品热| 少妇裸体淫交视频免费看高清| 人妻久久中文字幕网| 夜夜看夜夜爽夜夜摸| 舔av片在线| 免费看av在线观看网站| 亚洲经典国产精华液单| 国产av不卡久久| 久久精品国产亚洲网站| 伦精品一区二区三区| 亚洲综合色惰| 国产麻豆成人av免费视频| 国产亚洲精品av在线| 久久精品国产亚洲av香蕉五月| 亚洲国产色片| 国产成人一区二区在线| 亚洲性久久影院| 免费观看精品视频网站| 精品久久久久久成人av| 长腿黑丝高跟| 午夜福利在线在线| 亚洲精品一区av在线观看| 国产不卡一卡二| 一进一出抽搐动态| 亚洲欧美成人精品一区二区| 97人妻精品一区二区三区麻豆| 俺也久久电影网| 国产精品亚洲美女久久久| 午夜a级毛片| aaaaa片日本免费| 成人亚洲欧美一区二区av| 成人漫画全彩无遮挡| 久久久久精品国产欧美久久久| 日本黄大片高清| 嫩草影院新地址| 欧美最黄视频在线播放免费| 一a级毛片在线观看| 精品人妻熟女av久视频| 成人高潮视频无遮挡免费网站| 亚洲熟妇中文字幕五十中出| 麻豆av噜噜一区二区三区| 国产成人精品久久久久久| 欧美性感艳星| 国产成人a区在线观看| www日本黄色视频网| 精品人妻熟女av久视频| 亚洲经典国产精华液单| 国产精品爽爽va在线观看网站| 成年av动漫网址| 12—13女人毛片做爰片一| 最近中文字幕高清免费大全6| 日韩大尺度精品在线看网址| 一a级毛片在线观看| 精品一区二区三区视频在线| 在线播放国产精品三级| 波野结衣二区三区在线| 高清日韩中文字幕在线| 午夜老司机福利剧场| 99久久成人亚洲精品观看| 亚洲色图av天堂| 亚洲精品成人久久久久久| av视频在线观看入口| 一卡2卡三卡四卡精品乱码亚洲| 一个人看视频在线观看www免费| 久久国产乱子免费精品| 亚洲av第一区精品v没综合| 亚洲欧美精品综合久久99| 黄色配什么色好看| 亚洲精品成人久久久久久| 久久6这里有精品| 亚洲国产精品国产精品| 亚洲av五月六月丁香网| 国产精品久久久久久精品电影| 久久午夜福利片| 欧美潮喷喷水| 自拍偷自拍亚洲精品老妇| 国产一区二区在线av高清观看| 久久国内精品自在自线图片| 亚洲欧美日韩高清在线视频| 麻豆久久精品国产亚洲av| 九色成人免费人妻av| 久久久国产成人免费| 可以在线观看毛片的网站| 两个人的视频大全免费| 久久久久久久午夜电影| 1024手机看黄色片| 大又大粗又爽又黄少妇毛片口| 国产女主播在线喷水免费视频网站 | 免费av不卡在线播放| or卡值多少钱| 成人美女网站在线观看视频| 久久国产乱子免费精品| 国产欧美日韩一区二区精品| 亚洲国产精品合色在线| 亚洲国产欧美人成| 内地一区二区视频在线| 久久国产乱子免费精品| 日韩 亚洲 欧美在线| 免费观看在线日韩| 校园人妻丝袜中文字幕| 成人欧美大片| 日日撸夜夜添| 欧美另类亚洲清纯唯美| 神马国产精品三级电影在线观看| 最近手机中文字幕大全| 看十八女毛片水多多多| 男女之事视频高清在线观看| 国产在线男女| 日韩欧美三级三区| 直男gayav资源| 少妇人妻一区二区三区视频| 少妇人妻精品综合一区二区 | 色av中文字幕| 欧美日本亚洲视频在线播放| 麻豆av噜噜一区二区三区| 亚洲激情五月婷婷啪啪| 中文字幕精品亚洲无线码一区| 国产精品久久久久久久电影| 亚洲av熟女| 久久精品国产亚洲av涩爱 | 国产蜜桃级精品一区二区三区| 老司机午夜福利在线观看视频| av天堂在线播放| 精品人妻熟女av久视频| 国产亚洲精品久久久久久毛片| 日韩av在线大香蕉| 熟妇人妻久久中文字幕3abv| 观看美女的网站| 伦精品一区二区三区| 黄色视频,在线免费观看| 日韩一区二区视频免费看| 亚洲精品色激情综合| 日本在线视频免费播放| 国产一区二区在线观看日韩| 男人狂女人下面高潮的视频| 波多野结衣高清无吗| 1024手机看黄色片| 日本黄大片高清| 给我免费播放毛片高清在线观看| 成人特级黄色片久久久久久久| 成人无遮挡网站| 免费无遮挡裸体视频| 嫩草影院新地址| 国产成人freesex在线 | 精品久久久噜噜| 免费不卡的大黄色大毛片视频在线观看 | 国产精品三级大全| 国产综合懂色| 亚洲av.av天堂| 久久久a久久爽久久v久久| 一进一出好大好爽视频| 国产久久久一区二区三区| 搡女人真爽免费视频火全软件 | 午夜日韩欧美国产| 久久久久久国产a免费观看| 亚洲第一电影网av| 婷婷六月久久综合丁香| 亚洲人成网站在线播| 国产不卡一卡二| 一卡2卡三卡四卡精品乱码亚洲| 无遮挡黄片免费观看| 免费在线观看成人毛片| 日日啪夜夜撸| 国产中年淑女户外野战色| 男人舔奶头视频| 99热这里只有是精品在线观看| 最好的美女福利视频网| aaaaa片日本免费| 国产av在哪里看| 97超级碰碰碰精品色视频在线观看| 欧美日本亚洲视频在线播放| 少妇熟女欧美另类| 日本与韩国留学比较| 久久99热6这里只有精品| 亚洲乱码一区二区免费版| 日韩 亚洲 欧美在线| 此物有八面人人有两片| 欧美色欧美亚洲另类二区| 欧美日本亚洲视频在线播放| 成人鲁丝片一二三区免费|