• <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)度特征仿真分析
    丰满饥渴人妻一区二区三| 久热这里只有精品99| 免费在线观看视频国产中文字幕亚洲| 亚洲avbb在线观看| 97超级碰碰碰精品色视频在线观看| 国产av精品麻豆| 日本精品一区二区三区蜜桃| 777久久人妻少妇嫩草av网站| 亚洲一区高清亚洲精品| 中文字幕另类日韩欧美亚洲嫩草| av在线播放免费不卡| 亚洲三区欧美一区| 可以在线观看毛片的网站| 国产麻豆69| 满18在线观看网站| 黄色女人牲交| 国产一区二区三区视频了| 又黄又爽又免费观看的视频| 国产高清videossex| 女警被强在线播放| 大香蕉久久成人网| 免费人成视频x8x8入口观看| videosex国产| 国产在线观看jvid| 亚洲成人久久性| 丝袜美足系列| 成人永久免费在线观看视频| 国产精品野战在线观看 | 黄色毛片三级朝国网站| 亚洲人成电影免费在线| 中文字幕人妻丝袜制服| 午夜免费观看网址| 欧美亚洲日本最大视频资源| 国产亚洲精品久久久久久毛片| 午夜两性在线视频| 黄色a级毛片大全视频| 午夜成年电影在线免费观看| 女人精品久久久久毛片| 午夜精品在线福利| 亚洲精华国产精华精| 夜夜爽天天搞| 亚洲国产欧美一区二区综合| 在线天堂中文资源库| 成在线人永久免费视频| 亚洲av熟女| 国产av精品麻豆| 国产又色又爽无遮挡免费看| 亚洲一区二区三区色噜噜 | 一边摸一边做爽爽视频免费| 国产亚洲精品一区二区www| 亚洲 欧美一区二区三区| 久久伊人香网站| 中亚洲国语对白在线视频| 亚洲男人天堂网一区| 丰满饥渴人妻一区二区三| 亚洲色图综合在线观看| av片东京热男人的天堂| 国产精品 国内视频| 亚洲男人天堂网一区| 午夜免费成人在线视频| 在线十欧美十亚洲十日本专区| 超色免费av| 亚洲av成人av| 亚洲伊人色综图| 精品国产乱子伦一区二区三区| 老鸭窝网址在线观看| 午夜免费观看网址| 亚洲av成人不卡在线观看播放网| 精品一区二区三区视频在线观看免费 | 国产成人系列免费观看| 欧美日韩亚洲高清精品| 夜夜看夜夜爽夜夜摸 | 亚洲精品在线美女| 亚洲三区欧美一区| 色播在线永久视频| 日韩国内少妇激情av| 老熟妇乱子伦视频在线观看| 欧美 亚洲 国产 日韩一| 久久人人爽av亚洲精品天堂| 亚洲情色 制服丝袜| 国产精品电影一区二区三区| av网站在线播放免费| 亚洲国产中文字幕在线视频| 亚洲情色 制服丝袜| 午夜老司机福利片| 午夜精品在线福利| 免费在线观看黄色视频的| 久久九九热精品免费| 欧美日韩福利视频一区二区| 啪啪无遮挡十八禁网站| 欧美日韩亚洲高清精品| 午夜福利影视在线免费观看| 免费在线观看日本一区| 欧美中文综合在线视频| 色婷婷av一区二区三区视频| 亚洲人成电影观看| 高清黄色对白视频在线免费看| 一二三四社区在线视频社区8| 亚洲国产看品久久| 国产精品影院久久| 精品一区二区三区视频在线观看免费 | 中出人妻视频一区二区| 国产欧美日韩综合在线一区二区| 如日韩欧美国产精品一区二区三区| 亚洲av第一区精品v没综合| 国产一区二区在线av高清观看| 18禁国产床啪视频网站| 男人舔女人的私密视频| 亚洲专区中文字幕在线| 国产av在哪里看| 国产不卡一卡二| 亚洲久久久国产精品| 亚洲九九香蕉| 中文字幕人妻丝袜制服| 人人妻人人爽人人添夜夜欢视频| 18禁观看日本| 热re99久久精品国产66热6| 亚洲色图av天堂| www日本在线高清视频| av电影中文网址| 热re99久久国产66热| 一级毛片精品| 最新美女视频免费是黄的| 国产激情久久老熟女| 精品久久久久久电影网| 国产亚洲欧美98| 国产人伦9x9x在线观看| 男男h啪啪无遮挡| 激情在线观看视频在线高清| 嫩草影视91久久| 久久狼人影院| 不卡一级毛片| 精品久久久久久,| 18禁国产床啪视频网站| 国产精品日韩av在线免费观看 | 无遮挡黄片免费观看| 天堂俺去俺来也www色官网| 黄色成人免费大全| 自线自在国产av| 在线视频色国产色| 亚洲第一青青草原| 欧美av亚洲av综合av国产av| 久久久久久人人人人人| av天堂在线播放| 欧美在线黄色| 99热只有精品国产| av视频免费观看在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美一区二区三区在线观看| 91麻豆精品激情在线观看国产 | 国产成人精品无人区| 国产亚洲精品久久久久5区| 国产真人三级小视频在线观看| 色哟哟哟哟哟哟| 一级毛片高清免费大全| 亚洲一区二区三区不卡视频| 热re99久久精品国产66热6| 亚洲欧美激情在线| 少妇 在线观看| 亚洲色图av天堂| 欧美日韩av久久| 免费日韩欧美在线观看| 日韩国内少妇激情av| 欧美人与性动交α欧美软件| 在线十欧美十亚洲十日本专区| 搡老乐熟女国产| 国产1区2区3区精品| 水蜜桃什么品种好| 日韩免费av在线播放| 精品一区二区三卡| 亚洲精品国产精品久久久不卡| 成在线人永久免费视频| 成人免费观看视频高清| 大型av网站在线播放| 欧美乱妇无乱码| 99香蕉大伊视频| 熟女少妇亚洲综合色aaa.| 精品国产乱子伦一区二区三区| 欧美日韩瑟瑟在线播放| 国产成人欧美| 日韩 欧美 亚洲 中文字幕| av在线播放免费不卡| 黄片小视频在线播放| 久久久水蜜桃国产精品网| 久久久水蜜桃国产精品网| 国产精品爽爽va在线观看网站 | 最近最新免费中文字幕在线| 欧美人与性动交α欧美精品济南到| 午夜老司机福利片| 97超级碰碰碰精品色视频在线观看| svipshipincom国产片| 精品国产国语对白av| 夜夜看夜夜爽夜夜摸 | 免费女性裸体啪啪无遮挡网站| 日日干狠狠操夜夜爽| 国产精品av久久久久免费| 中文字幕另类日韩欧美亚洲嫩草| 精品免费久久久久久久清纯| 亚洲欧美一区二区三区黑人| 男女之事视频高清在线观看| 交换朋友夫妻互换小说| 亚洲一区二区三区欧美精品| 男人的好看免费观看在线视频 | 99精国产麻豆久久婷婷| 亚洲av第一区精品v没综合| 国产熟女午夜一区二区三区| 国产人伦9x9x在线观看| 亚洲自偷自拍图片 自拍| 亚洲av片天天在线观看| av电影中文网址| 精品国产亚洲在线| av在线播放免费不卡| 午夜福利免费观看在线| 超色免费av| 亚洲国产中文字幕在线视频| 久久人妻av系列| 国产成+人综合+亚洲专区| 999久久久精品免费观看国产| 亚洲少妇的诱惑av| 免费少妇av软件| 他把我摸到了高潮在线观看| 国产三级黄色录像| 国产成人精品无人区| 免费在线观看黄色视频的| 正在播放国产对白刺激| 男女高潮啪啪啪动态图| 在线观看舔阴道视频| 一二三四在线观看免费中文在| 丝袜在线中文字幕| 国产精品久久久人人做人人爽| 午夜福利欧美成人| 男女做爰动态图高潮gif福利片 | 在线国产一区二区在线| 免费在线观看日本一区| 久久国产亚洲av麻豆专区| 亚洲自偷自拍图片 自拍| 国产精品1区2区在线观看.| 亚洲三区欧美一区| 一进一出抽搐gif免费好疼 | 91国产中文字幕| 久9热在线精品视频| 亚洲欧美激情综合另类| av网站在线播放免费| 久久影院123| 中文字幕精品免费在线观看视频| 欧美激情久久久久久爽电影 | 黄色 视频免费看| 久久国产精品男人的天堂亚洲| 黄色视频不卡| 日韩视频一区二区在线观看| 性色av乱码一区二区三区2| 欧美在线一区亚洲| 两性夫妻黄色片| 日本一区二区免费在线视频| 两个人看的免费小视频| 窝窝影院91人妻| 成熟少妇高潮喷水视频| 国产亚洲欧美98| 丁香六月欧美| 亚洲国产欧美日韩在线播放| 久久精品国产亚洲av香蕉五月| 日韩欧美免费精品| 国产麻豆69| 亚洲国产欧美网| 伊人久久大香线蕉亚洲五| 天堂动漫精品| 99国产精品一区二区蜜桃av| 国产精华一区二区三区| 老司机福利观看| netflix在线观看网站| 成人特级黄色片久久久久久久| 大型黄色视频在线免费观看| 久热这里只有精品99| 天堂中文最新版在线下载| 无人区码免费观看不卡| 亚洲色图综合在线观看| 日本黄色视频三级网站网址| 国产xxxxx性猛交| 日韩欧美三级三区| 男女下面插进去视频免费观看| 国产极品粉嫩免费观看在线| 两个人看的免费小视频| 精品国产超薄肉色丝袜足j| 精品国产美女av久久久久小说| 最近最新免费中文字幕在线| www.精华液| 丁香六月欧美| 国产精品国产高清国产av| 国产99白浆流出| 最好的美女福利视频网| 国产精品久久久人人做人人爽| 一区福利在线观看| 99在线人妻在线中文字幕| 一区二区三区精品91| 国产精品一区二区免费欧美| 久久久久国产精品人妻aⅴ院| 女人高潮潮喷娇喘18禁视频| 天堂动漫精品| 99精国产麻豆久久婷婷| 成年人黄色毛片网站| 黑人猛操日本美女一级片| 国产亚洲精品一区二区www| 99国产极品粉嫩在线观看| 五月开心婷婷网| 精品国产一区二区久久| 免费一级毛片在线播放高清视频 | 丰满迷人的少妇在线观看| 又黄又爽又免费观看的视频| 午夜影院日韩av| 久久中文字幕人妻熟女| 啪啪无遮挡十八禁网站| 无人区码免费观看不卡| 真人做人爱边吃奶动态| 国产深夜福利视频在线观看| 国产午夜精品久久久久久| 夜夜爽天天搞| 久久人妻av系列| 亚洲中文字幕日韩| 19禁男女啪啪无遮挡网站| 久久天躁狠狠躁夜夜2o2o| 桃色一区二区三区在线观看| 搡老熟女国产l中国老女人| 亚洲精品在线观看二区| 欧美日韩乱码在线| 精品久久蜜臀av无| 欧美日韩瑟瑟在线播放| 亚洲第一青青草原| 亚洲五月婷婷丁香| 桃色一区二区三区在线观看| 黑人欧美特级aaaaaa片| 不卡av一区二区三区| 亚洲欧美精品综合一区二区三区| 男人舔女人的私密视频| 亚洲久久久国产精品| 国产精品av久久久久免费| 99riav亚洲国产免费| √禁漫天堂资源中文www| 97人妻天天添夜夜摸| 午夜福利一区二区在线看| av超薄肉色丝袜交足视频| 亚洲精品国产色婷婷电影| 欧美乱码精品一区二区三区| 欧美日韩av久久| 国产av一区二区精品久久| 久久久久国内视频| 午夜成年电影在线免费观看| 欧美国产精品va在线观看不卡| 日本免费一区二区三区高清不卡 | 亚洲午夜精品一区,二区,三区| 韩国精品一区二区三区| 欧美激情极品国产一区二区三区| 精品熟女少妇八av免费久了| 国产激情欧美一区二区| 中文亚洲av片在线观看爽| 免费在线观看完整版高清| 国产熟女午夜一区二区三区| 久久 成人 亚洲| 国产精品亚洲一级av第二区| 巨乳人妻的诱惑在线观看| 国产精品 国内视频| 成熟少妇高潮喷水视频| 久久精品国产综合久久久| 国产av在哪里看| 久久精品亚洲av国产电影网| 成人三级做爰电影| 精品久久久精品久久久| 欧美亚洲日本最大视频资源| 在线观看一区二区三区| 黄频高清免费视频| 亚洲欧美激情综合另类| 可以在线观看毛片的网站| 欧美日韩黄片免| 另类亚洲欧美激情| 免费人成视频x8x8入口观看| 国产精品亚洲av一区麻豆| 色播在线永久视频| 黑人巨大精品欧美一区二区蜜桃| 91精品三级在线观看| 久久精品成人免费网站| 国产精品 欧美亚洲| 亚洲一区中文字幕在线| 国产在线观看jvid| 51午夜福利影视在线观看| 中出人妻视频一区二区| 日韩国内少妇激情av| 欧美国产精品va在线观看不卡| 丝袜在线中文字幕| 久久久久久久午夜电影 | 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 91成年电影在线观看| 午夜免费激情av| 欧美大码av| 亚洲成人免费电影在线观看| 国产精品免费一区二区三区在线| 青草久久国产| 欧美久久黑人一区二区| 亚洲欧美精品综合一区二区三区| 午夜a级毛片| 国产精品久久视频播放| 国产高清国产精品国产三级| 日韩中文字幕欧美一区二区| 午夜福利,免费看| 日本撒尿小便嘘嘘汇集6| 国产成+人综合+亚洲专区| avwww免费| tocl精华| 国产亚洲精品综合一区在线观看 | 神马国产精品三级电影在线观看 | 国产精品久久久av美女十八| 一级片免费观看大全| 天堂影院成人在线观看| 久久人人97超碰香蕉20202| 中文字幕人妻丝袜制服| 午夜久久久在线观看| 欧美日韩亚洲高清精品| 欧美日韩亚洲综合一区二区三区_| 精品久久久久久电影网| 久久午夜综合久久蜜桃| 亚洲欧美日韩另类电影网站| 日韩欧美三级三区| 99在线视频只有这里精品首页| 男男h啪啪无遮挡| 精品电影一区二区在线| 高清黄色对白视频在线免费看| 天堂影院成人在线观看| 亚洲全国av大片| 午夜福利,免费看| 狠狠狠狠99中文字幕| 熟女少妇亚洲综合色aaa.| 丁香六月欧美| 免费看a级黄色片| 一夜夜www| 老司机午夜福利在线观看视频| 国产精品国产高清国产av| 在线观看免费视频网站a站| 满18在线观看网站| 国产97色在线日韩免费| 欧美激情久久久久久爽电影 | 精品国产美女av久久久久小说| 精品久久久久久电影网| 又黄又爽又免费观看的视频| 久久精品aⅴ一区二区三区四区| 精品福利永久在线观看| 国产精品秋霞免费鲁丝片| 好男人电影高清在线观看| 又紧又爽又黄一区二区| 99国产精品免费福利视频| 97超级碰碰碰精品色视频在线观看| 老司机在亚洲福利影院| 午夜免费观看网址| 99在线视频只有这里精品首页| 超碰成人久久| 黄片小视频在线播放| 69精品国产乱码久久久| 免费在线观看视频国产中文字幕亚洲| 国产欧美日韩一区二区三| 日韩一卡2卡3卡4卡2021年| 国产精品98久久久久久宅男小说| 长腿黑丝高跟| 欧美精品一区二区免费开放| 成人国语在线视频| 最新美女视频免费是黄的| 黄色成人免费大全| 欧美黑人欧美精品刺激| 一个人观看的视频www高清免费观看 | 亚洲少妇的诱惑av| 日本a在线网址| 国产日韩一区二区三区精品不卡| 大型黄色视频在线免费观看| 男女之事视频高清在线观看| 丰满迷人的少妇在线观看| 男女午夜视频在线观看| 大香蕉久久成人网| 亚洲专区字幕在线| 黄色视频不卡| 国产精品永久免费网站| 男人舔女人的私密视频| 最近最新中文字幕大全免费视频| av天堂久久9| 女人高潮潮喷娇喘18禁视频| 大香蕉久久成人网| 美女福利国产在线| 啦啦啦免费观看视频1| 成人黄色视频免费在线看| 一边摸一边抽搐一进一出视频| 亚洲精品粉嫩美女一区| 欧美丝袜亚洲另类 | 男女午夜视频在线观看| 999久久久精品免费观看国产| 91在线观看av| 99久久久亚洲精品蜜臀av| 丰满的人妻完整版| 满18在线观看网站| 十分钟在线观看高清视频www| 国产精品国产高清国产av| 99国产极品粉嫩在线观看| 天天影视国产精品| av免费在线观看网站| 一二三四在线观看免费中文在| 日韩精品青青久久久久久| 亚洲 国产 在线| 亚洲三区欧美一区| 亚洲国产精品999在线| 精品福利永久在线观看| 性少妇av在线| 美女国产高潮福利片在线看| 亚洲精品久久成人aⅴ小说| 亚洲专区中文字幕在线| 午夜福利,免费看| 后天国语完整版免费观看| av国产精品久久久久影院| 波多野结衣av一区二区av| 久久久久久久久久久久大奶| 久久青草综合色| 久久精品aⅴ一区二区三区四区| 97超级碰碰碰精品色视频在线观看| 超碰成人久久| 免费高清视频大片| 国产亚洲欧美精品永久| 黄网站色视频无遮挡免费观看| 韩国av一区二区三区四区| 天堂俺去俺来也www色官网| 国产成人精品在线电影| 中文字幕最新亚洲高清| 天天躁夜夜躁狠狠躁躁| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 村上凉子中文字幕在线| 丝袜美腿诱惑在线| 十分钟在线观看高清视频www| 在线观看舔阴道视频| 亚洲专区中文字幕在线| 亚洲精品美女久久av网站| 18禁裸乳无遮挡免费网站照片 | 欧美另类亚洲清纯唯美| 中文字幕精品免费在线观看视频| 99久久久亚洲精品蜜臀av| 亚洲中文av在线| 99久久综合精品五月天人人| 免费在线观看日本一区| 精品一区二区三区视频在线观看免费 | 亚洲aⅴ乱码一区二区在线播放 | 国产精品久久久人人做人人爽| 国产一区二区在线av高清观看| 这个男人来自地球电影免费观看| 欧美另类亚洲清纯唯美| 亚洲专区国产一区二区| 国产又爽黄色视频| 五月开心婷婷网| 日本一区二区免费在线视频| 亚洲精品国产色婷婷电影| 亚洲欧美激情在线| 久久久国产精品麻豆| 极品人妻少妇av视频| 一本大道久久a久久精品| 亚洲自拍偷在线| 久久人妻福利社区极品人妻图片| 国产麻豆69| 精品国产国语对白av| 国产精品久久久av美女十八| 欧美在线一区亚洲| 亚洲欧美精品综合久久99| 亚洲午夜精品一区,二区,三区| 精品国产美女av久久久久小说| 久久久久九九精品影院| 久久人人爽av亚洲精品天堂| 精品一区二区三区视频在线观看免费 | 欧美 亚洲 国产 日韩一| 精品一区二区三卡| 在线看a的网站| 法律面前人人平等表现在哪些方面| 黄色成人免费大全| 久久天堂一区二区三区四区| 丝袜人妻中文字幕| 脱女人内裤的视频| 精品国产亚洲在线| 啦啦啦 在线观看视频| 日本一区二区免费在线视频| 一本综合久久免费| 国产不卡一卡二| 国产精华一区二区三区| 啦啦啦在线免费观看视频4| 最好的美女福利视频网| 午夜亚洲福利在线播放| av网站在线播放免费| 在线观看66精品国产| 88av欧美| 国产精品亚洲av一区麻豆| 天堂俺去俺来也www色官网| 黑丝袜美女国产一区| 亚洲精品国产一区二区精华液| 18禁观看日本| 久久99一区二区三区| 高潮久久久久久久久久久不卡| 97碰自拍视频| 免费人成视频x8x8入口观看| av天堂久久9| 亚洲精品一区av在线观看| 69av精品久久久久久| 日韩大尺度精品在线看网址 | 黄色怎么调成土黄色| 在线观看一区二区三区激情| 国产精品偷伦视频观看了| av电影中文网址| 黄色成人免费大全| videosex国产| 日日夜夜操网爽| 神马国产精品三级电影在线观看 | 中文欧美无线码| 精品一区二区三区av网在线观看| 精品电影一区二区在线| 欧美色视频一区免费| 成人av一区二区三区在线看| av超薄肉色丝袜交足视频|