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

    旋轉(zhuǎn)加速度計重力梯度儀數(shù)據(jù)處理方法

    2016-09-21 03:37:40錢學(xué)武蔡體菁
    關(guān)鍵詞:重力梯度加速度計小波

    錢學(xué)武  蔡體菁

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

    ?

    旋轉(zhuǎn)加速度計重力梯度儀數(shù)據(jù)處理方法

    錢學(xué)武 蔡體菁

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

    為了有效去除旋轉(zhuǎn)加速度計重力梯度儀輸出信號中的各種干擾噪聲,提出了一種提取重力梯度信號的有效方法.首先對重力梯度儀輸出信號進行故障診斷,然后采用基于Chebyshev最佳一致逼近法設(shè)計的帶通濾波器對故障診斷后的信號進行濾波,并對濾波后的重力梯度信號進行梯度解調(diào),最后采用dmey小波基函數(shù)強制閾值方法對解調(diào)后的重力梯度信息做進一步去噪處理.在重力梯度半物理仿真平臺上進行了仿真試驗測試,結(jié)果表明,所提方法可以有效降低有用信號均方差,不會造成數(shù)據(jù)丟失和信號偏移,提高了重力梯度測量精度.

    重力梯度儀;旋轉(zhuǎn)加速度計;梯度解調(diào);帶通濾波;小波去噪

    高精度重力梯度測量對于空間科學(xué)、地球科學(xué)、地質(zhì)科學(xué)、能源勘探、慣性導(dǎo)航等領(lǐng)域的研究具有重要意義[1-4].目前,重力梯度儀主要有旋轉(zhuǎn)加速度計重力梯度儀、超導(dǎo)重力梯度儀、冷原子重力梯度儀、靜電重力梯度儀,以及基于微機械加工(MEMS)技術(shù)的重力梯度儀等,其中旋轉(zhuǎn)加速度計重力梯度儀是目前唯一成功用于機載/船載動基座上的重力梯度測量儀器[5-8].

    旋轉(zhuǎn)加速度計重力梯度儀由多種機械裝置組成,加速度計性能、加速度計安裝精度、旋轉(zhuǎn)機構(gòu)穩(wěn)定性、溫控性能、平臺穩(wěn)定性、電路系統(tǒng)穩(wěn)定性等因素都會對重力梯度信號造成嚴重影響[9-10],這些因素會使重力梯度儀輸出信號包含各種噪聲,降低了測量精度,因此研究重力梯度信號處理方法愈顯重要.國外學(xué)者已對航空重力梯度儀數(shù)據(jù)處理方法進行了相關(guān)研究[11-13],但都是采用基于計算機模擬仿真的處理方法,可信度不高.由于傳統(tǒng)數(shù)字濾波器具有平滑效應(yīng),會對重力梯度信息起到平滑作用,并且無法濾除重力梯度同頻帶內(nèi)的噪聲,相比而言,小波分析提供了一種自適應(yīng)時域和頻域同時局部化的多分辨率分析方法,可以很好地體現(xiàn)出信號的非平穩(wěn)特性,為信號處理提供了一種全新途徑[14-16].本文針對重力梯度信號特點,提出了一種基于小波分析的去除重力梯度干擾噪聲方法,仿真實驗結(jié)果表明該處理方法可以有效去除重力梯度信號中的干擾噪聲,提高了重力梯度測量精度.

    1 重力梯度信號處理方法

    在實際重力梯度測量過程中,由于儀器性能和外部環(huán)境的干擾,重力梯度儀輸出信號中往往含有多種較大幅度干擾噪聲,并且存在數(shù)據(jù)丟失或數(shù)據(jù)錯誤等情況.針對重力梯度信號特點,本文提出一種重力梯度信號有效處理方法:首先對重力梯度儀輸出信號進行故障診斷處理,消除故障數(shù)據(jù)點;然后采用數(shù)字帶通濾波器對故障診斷處理后的重力梯度信號進行濾波處理并進行重力梯度解調(diào);最后采用小波分解方法對含噪重力梯度信息進行去噪處理,得到高精度重力梯度信息.重力梯度信號處理過程如圖1所示.

    圖1 重力梯度信號處理框圖

    1.1FIR數(shù)字帶通濾波器

    重力梯度信號信噪比非常低,為了去除噪聲獲得有用信號,需要設(shè)計合適的帶通濾波器.現(xiàn)有的數(shù)字濾波器設(shè)計方法中,在保證濾波器所有性能指標滿足的前提下,采用Chebyshev最佳一致逼近方法設(shè)計的數(shù)字濾波器階數(shù)最?。碚撋?FIR濾波器的通帶帶寬越窄,階數(shù)越高,濾波效果越好,但濾波器階數(shù)太高會出現(xiàn)濾波效果不明顯,甚至出現(xiàn)濾波不穩(wěn)定現(xiàn)象.濾波器階數(shù)可根據(jù)Kaiser[17]提出的經(jīng)驗公式進行近似估計:

    (1)

    式中,N為濾波器階數(shù);過渡帶帶寬Δf=(ωs-ωp)/(2π),ωp和ωs分別為濾波器通帶和阻帶截止角頻率;δp和δs分別為濾波器通帶和阻帶的紋波峰值.本文采用Remez算法設(shè)計濾波器,根據(jù)重力梯度信號特征,并反復(fù)進行實驗測試驗證,最終得到一組濾波效果較好的帶通濾波器參數(shù).本文選取的濾波器參數(shù)為:中心頻率為0.5Hz,通帶帶寬為0.02Hz,過渡帶帶寬為0.01Hz,濾波器階數(shù)為800.

    1.2重力梯度解調(diào)

    旋轉(zhuǎn)加速度計重力梯度儀理想輸出信號表達式為

    Eout=2RKKI[(Γyy-Γxx)sin2ωt+2Γxycos2ωt]

    (2)

    式中,Eout為重力梯度儀輸出信號;R為加速度計質(zhì)量中心到GGI圓盤中心的距離;K為重力梯度信號放大增益;KI為加速度計標度因數(shù);ω為GGI圓盤旋轉(zhuǎn)角頻率;Γyy,Γxx和Γxy分別為相應(yīng)方向上的重力梯度張量元素.只要對式(2)在2ω頻率上進行幅值解調(diào)就可以得到重力梯度分量Γyy-Γxx和Γxy,解調(diào)公式為

    (3)

    式中, T為重力梯度儀圓盤旋轉(zhuǎn)周期.

    1.3小波去噪

    小波分析屬于時頻分析的一種,它能夠同時給出信號的時域和頻域特征,有效區(qū)分出信號中的突變成分和干擾噪聲,從而實現(xiàn)信號分離和降噪.設(shè)重力梯度信號為Γ(t),被噪聲污染的重力梯度信號為G(t),重力梯度信號模型可表示為

    G(t)=Γ(t)+e(t)

    (4)

    (5)

    式中,j,k分別為離散化的縮放因子和平移因子;ψ(t)為小波基函數(shù);ψ*(t)為ψ(t)的共軛.利用式(5)可得到信號G(t)在第j尺度上的小波系數(shù).采用Mallat算法對信號G(t)進行小波分解和信號重構(gòu).

    重力梯度信號本質(zhì)上是一種非平穩(wěn)低頻隨機信號,經(jīng)過小波分解后,在某個分解尺度上的小波系數(shù)具有較大幅值,而噪聲頻帶比重力梯度信號頻帶大得多.針對這種情況,在重力梯度數(shù)據(jù)去噪階段采用強制閾值去噪方法.該方法僅對某個尺度下的低頻系數(shù)進行小波變換來重建信號,直接舍棄其他尺度上的系數(shù)分量,從而達到信號降噪的目的.在強制閾值濾波過程中,采用dmey小波基函數(shù),這是由于該小波基具有很好的正則性、緊支撐性和平滑性,并且適合重力梯度信號的特點.對重力梯度信號進行10層分解,經(jīng)過比較相鄰層信號與真實信號的變化趨勢發(fā)現(xiàn),第6層低頻信號與真實信號比較接近,其他層與真實信號相差較大,故選擇第6層低頻信號作為重力梯度信息的最終處理結(jié)果.

    2 半物理仿真實驗與結(jié)果分析

    根據(jù)旋轉(zhuǎn)加速度計重力梯度儀工作原理,采用高性能計算機、可編程高精度電流源、低噪聲電流放大器、低噪聲電壓放大器、多路切換開關(guān)和高精度數(shù)字電壓表等高精密儀器組建了重力梯度半物理仿真系統(tǒng)[18].利用該仿真系統(tǒng)可以完成信號產(chǎn)生、信號放大、信號濾波、重力梯度解調(diào)等仿真實驗.仿真系統(tǒng)采用LabVIEW軟件開發(fā).重力梯度半物理仿真系統(tǒng)實物圖如圖2所示.

    圖2 重力梯度半物理仿真系統(tǒng)實物圖

    采用高密度的鎢合金均質(zhì)圓球作為重力梯度引力裝置,利用所開發(fā)的重力梯度半物理仿真系統(tǒng)對重力梯度信號進行仿真測試. 實驗參數(shù)如下:圓球密度為18 000kg/m3,質(zhì)量為1 000kg,圓球質(zhì)心坐標為(0.8,0,0)m,圓盤旋轉(zhuǎn)周期為4s,加速度計對的基線距離為0.2m,加速度計標度因數(shù)為10mA/g,信號放大增益為5×106,信號采樣頻率為10Hz.根據(jù)萬有引力定律,可以計算出在上述參數(shù)下圓球引起的重力梯度分量理論值分別為Γyy-Γxx=-393E,Γxy=0E.將從旋轉(zhuǎn)加速度計重力梯度儀樣機上采集到的原始加速度計數(shù)據(jù)作為重力梯度半物理仿真系統(tǒng)信號源,對重力梯度半物理仿真系統(tǒng)輸出信號進行數(shù)據(jù)采集和數(shù)據(jù)故障診斷處理,然后采用帶通濾波器對故障診斷后的重力梯度信號進行濾波.濾波前后的重力梯度信號時域波形如圖3所示,對應(yīng)的頻譜如圖4所示.

    圖3 重力梯度信號濾波前后信號波形

    圖4 重力梯度信號濾波前后信號頻譜

    從圖3可看出,重力梯度信號全部被噪聲淹沒,完全看不到重力梯度信號,濾波后的信號幅值大幅降低.從圖4可看出,0.5Hz的重力梯度信號較為明顯,除0.5Hz外的低高頻干擾諧波幅度較大,而帶通濾波后的重力梯度信號頻譜只含有0.5Hz窄帶信號和通帶帶寬內(nèi)殘留的少量噪聲信號,0.5Hz外的高幅值低高頻噪聲基本被濾除.利用式(3)對帶通濾波后的重力梯度信號進行梯度解調(diào),然后對解調(diào)得到的重力梯度信息進行100s平滑濾波、200s平滑濾波和小波強制閾值去噪處理,不同處理方法的結(jié)果如表1所示,波形對比如圖5所示.

    從表1和圖5可看出,采用平滑濾波處理方法

    表1 重力梯度不同處理方法結(jié)果對比 E

    (b) 重力梯度分量Γxy

    雖然可以有效去除噪聲,但隨著平滑時間的延長,重力梯度信息丟失越多,且平滑濾波后的均值發(fā)生了明顯偏移.而采用小波去噪方法不會出現(xiàn)數(shù)據(jù)丟失和均值偏移現(xiàn)象,并且可以明顯降低重力梯度信息均方差.同時還發(fā)現(xiàn),重力梯度分量Γyy- Γxx和Γxy均值分別約為-430和-22E,而理論值分別為-393和0E,重力梯度測量值與理論值相差較大,這主要是由于在0.5Hz上存在信號串擾.造成信號串擾的原因很多,可通過采用重力梯度標定方法進行重力梯度補償解決.

    3 結(jié)語

    本文針對重力梯度信號特點,提出了一種去除重力梯度干擾噪聲的有效處理方法.首先對重力梯度信號進行故障診斷處理,然后采用帶通濾波器對故障診斷后的信號進行濾波,此濾波方法可以把重力梯度信號頻率外的低高頻干擾噪聲濾除,提高了重力梯度信號信噪比;接著對濾波后的重力梯度信號進行梯度解調(diào);最后采用小波去噪方法對解調(diào)結(jié)果做進一步處理.為了驗證所提處理方法的有效性,在重力梯度半物理仿真系統(tǒng)上進行了仿真實驗測試.仿真結(jié)果表明,小波去噪方法可以顯著提高重力梯度測量精度,并且不會造成數(shù)據(jù)丟失和信號偏移,為實現(xiàn)高精度重力梯度測量提供了一種方法.

    References)

    [1]Difrancesco D. Advances and challenges in the development and deployment of gravity gradiometer systems[C/OL]//EGM2007InternationalWorkshop. Capri, Italy, 2007. http://www.earthdoc.org/publication/publicationdetails/?publication=41199.

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

    [3]Christensen A N, Galder C V, Dransfield M. Improved resolution of fixed-wing airborne gravity gradiometer surveys[C]//2014SEGAnnualMeeting. Denver, Colorado, USA, 2014: 1319-1323. DOI: 10.1190/segam2014-0713.1.

    [4]DiFrancesco D, Meyer T, Christensen A, et al. Gravity gradiometry—Today and tomorrow[C]//11thSAGABiennialTechnicalMeetingandExhibition. Swaziland, 2009: 80-83.

    [5]Hodges G, Dransfield M H, Shei T C. The Falcon airborne gravity gradiometer for engineering applications[C]//SymposiumontheApplicationofGeophysicstoEngineeringandEnvironmentalProblems2010. Keystone, Colorado, USA, 2010: 443-447. DOI:10.4133/1.3445467.

    [6]Christensen A N, Hodges G. HeliFALCON?airborne gravity gradiometer data acquisition in rugged terrain[C]//Proceedingsofthe11thSEGJInternationalSymposium. Yokohama, Japan, 2013:140-145. DOI:10.1190/segj112013-036.

    [7]Difrancesco D, Grierson A, Kaputa D, et al. Gravity gradiometer systems—Advances and challenges[J].GeophysicalProspecting, 2009, 57(4): 615-623.

    [8]Dransfield M H, Christensen A N. Performance of airborne gravity gradiometers[J].TheLeadingEdge, 2013, 32(8):908-922. DOI:10.1190/tle32080908.1.

    [9]Jekeli C. The gravity gradiometer survey system (GGSS)[J].EosTransactionsAmericanGeophysicalUnion, 1988, 69(8):105, 116-117.

    [10]Jekeli C. Statistical analysis of moving-base gravimetry and gravity gradiometry[R]. Columbus,Ohio,USA: Geodetic and GeoInfomation Science, The Ohio State University, 2003.

    [11]Cesar J, Lyrio S O. Wavelet denoising of gravity gradiometry data[C]//SEGInternationalExpositionand71stAnnualMeeting. San Antonio, USA, 2001:1474-1477. DOI:10.1190/1.1816384.

    [12]Carlos D U, Braga M A, Galbiatti H F, et al. Airborne gravity gradiometry-data processing and interpretation [J].RevistaBrasileiraDeGeofǐsica, 2013, 31(3):427-453.

    [13]Christensen A N, Dransfield M H, Galder C V. Noise and repeatability of airborne gravity gradiometry [J].FirstBreak, 2015, 33:55-63.

    [14]Barnes G, Lumley J. Processing gravity gradient data [J].Geophysics, 2011, 76(2):133-147.

    [15]Jekeli C, Abt T L. The statistical performance of the matched filter for anomaly detection using gravity gradients[R]. Columbus, Ohio, USA: Division of Geodetic Science, Ohio State University, 2010.

    [16]de Oliveira Lyrio J C S, Tenorio L, Li Y, et al. Efficient automatic denoising of gravity gradiometry data[J].Geophysics, 2004, 69(3):772-782. DOI:10.1190/1.1759463.

    [17]Kaiser J F. Nonrecursive digital filter design using the Io-sinh window function[C]//IEEEInternationalSymposiumonCircuits&Systems. San Francisco, CA, USA, 1974:20-23.

    [18]蔡體菁,錢學(xué)武,丁昊.旋轉(zhuǎn)加速度計重力梯度儀重力梯度信號仿真[J].物探與化探,2015,39(S1):76-79.

    Cai Tijing, Qian Xuewu, Ding Hao. Signal simulation of gravity gradiometer of rotating accelerometer[J].GeophysicalandGeochemicalExploration, 2015, 39(S1):76-79. (in Chinese)

    DOI:10.3969/j.issn.1001-0505.2016.04.007

    Data processing method for rotating accelerometer gravity gradiometer

    Qian Xuewu Cai Tijing

    (School of Instrument Science and Engineering, Southeast University, Nanjing 210096, China)

    In order to remove the noise in the output signals of the rotating accelerometer gravity gradiometer instrument (GGI), an effective method for extracting the gravity gradient signal is presented. First, the fault diagnosis is carried out to preprocess the output signals of the GGI. Then a band-passing filter based on the Chebyshev optimal uniform approximation method is used to filtrate the diagnosed signals, and the output signals of the band-pass filter are demodulated in the next step. Finally, the gravity gradient is further denoised by the compulsory threshold method based on the demy wavelet function. The proposed method is tested on the gravity gradient hard-in-the-loop simulation platform. The simulation results show that the method can effectively reduce the mean square deviation of the useful signals and avoid data loss or signal offset. Therefore, the method can improve the precision of the gravity gradiometer.

    gravity gradiometer; rotating accelerometer; gradient demodulation; band-passing filtering; wavelet de-noising

    10.3969/j.issn.1001-0505.2016.04.006

    2015-12-28.作者簡介: 錢學(xué)武(1981—),男,博士生;蔡體菁(聯(lián)系人),男,博士,教授,博士生導(dǎo)師,caitij@seu.edu.cn.

    國家高技術(shù)研究發(fā)展計劃(863計劃)資助項目(2011AA060501).

    10.3969/j.issn.1001-0505.2016.04.006.

    U666.1

    A

    1001-0505(2016)04-0708-05

    引用本文: 錢學(xué)武,蔡體菁.旋轉(zhuǎn)加速度計重力梯度儀數(shù)據(jù)處理方法[J].東南大學(xué)學(xué)報(自然科學(xué)版),2016,46(4):708-712.

    猜你喜歡
    重力梯度加速度計小波
    基于加速度計的起重機制動下滑量測量
    構(gòu)造Daubechies小波的一些注記
    科技風(2021年19期)2021-09-07 14:04:29
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    基于改進的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    基于遺傳算法的加速度計免轉(zhuǎn)臺標定方法
    常溫下硅微諧振加速度計零偏穩(wěn)定性的提高
    旋轉(zhuǎn)加速度計重力梯度儀標定方法
    利用地形數(shù)據(jù)計算重力梯度張量的直接積分法
    星載重力梯度儀的研究發(fā)展
    四加速度計組合的測試方法研究
    51国产日韩欧美| 一级毛片aaaaaa免费看小| 国产精品不卡视频一区二区| 麻豆久久精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆| 久久综合国产亚洲精品| 三级国产精品片| 亚洲在线自拍视频| 伦理电影大哥的女人| 欧美精品一区二区大全| 免费不卡的大黄色大毛片视频在线观看 | 国语自产精品视频在线第100页| 亚洲成人精品中文字幕电影| 偷拍熟女少妇极品色| 国产老妇伦熟女老妇高清| 高清在线视频一区二区三区 | 欧美性猛交黑人性爽| 久久久色成人| 国产淫语在线视频| 久久精品人妻少妇| 中文字幕熟女人妻在线| 亚洲一区高清亚洲精品| a级一级毛片免费在线观看| 日本-黄色视频高清免费观看| 成人国产麻豆网| 久久精品久久精品一区二区三区| 中文字幕制服av| 国产伦一二天堂av在线观看| 色噜噜av男人的天堂激情| 国产精品女同一区二区软件| 久久草成人影院| 伦精品一区二区三区| 久久久午夜欧美精品| 国产精品三级大全| 一级黄片播放器| 成人欧美大片| 2021天堂中文幕一二区在线观| 淫秽高清视频在线观看| 免费看日本二区| 国产午夜精品论理片| 国产精品国产三级专区第一集| 高清日韩中文字幕在线| 久久99热这里只有精品18| 中文字幕av成人在线电影| 美女脱内裤让男人舔精品视频| 亚洲欧美精品自产自拍| 亚洲精品乱码久久久久久按摩| 黄片无遮挡物在线观看| 日本wwww免费看| 九色成人免费人妻av| 深爱激情五月婷婷| 久久午夜福利片| 99在线视频只有这里精品首页| 欧美3d第一页| 午夜久久久久精精品| 欧美日本视频| 在线a可以看的网站| 国产三级在线视频| 日本爱情动作片www.在线观看| 干丝袜人妻中文字幕| 国产亚洲精品av在线| 欧美又色又爽又黄视频| 欧美一区二区亚洲| 亚洲av不卡在线观看| 国产探花极品一区二区| 婷婷六月久久综合丁香| 99久久精品国产国产毛片| av在线老鸭窝| 在线观看av片永久免费下载| 免费看美女性在线毛片视频| 久久久久久久国产电影| 午夜福利在线观看免费完整高清在| 久久久成人免费电影| 精品免费久久久久久久清纯| 精品久久久久久久久久久久久| 亚洲成人av在线免费| 99久久精品热视频| 国产精品久久视频播放| 男人的好看免费观看在线视频| 少妇熟女欧美另类| 久热久热在线精品观看| 免费av不卡在线播放| 精品一区二区免费观看| 1024手机看黄色片| 国产淫语在线视频| 2021天堂中文幕一二区在线观| 一级毛片久久久久久久久女| 内地一区二区视频在线| 日本免费在线观看一区| 91久久精品国产一区二区三区| 三级国产精品片| 国产精品99久久久久久久久| 亚洲av成人精品一二三区| 18禁在线无遮挡免费观看视频| 亚洲精品一区蜜桃| 老司机影院毛片| 美女黄网站色视频| 国产乱来视频区| 国产伦在线观看视频一区| 国产 一区精品| 69人妻影院| 中文字幕免费在线视频6| 三级经典国产精品| 高清av免费在线| 亚洲成人久久爱视频| 18+在线观看网站| 免费看a级黄色片| 最近中文字幕高清免费大全6| 如何舔出高潮| 成年女人看的毛片在线观看| 淫秽高清视频在线观看| 亚洲美女搞黄在线观看| 三级国产精品欧美在线观看| 综合色丁香网| 亚洲精品一区蜜桃| 亚洲欧美日韩东京热| 午夜福利在线在线| 我要看日韩黄色一级片| 在线免费十八禁| 精品久久久久久久人妻蜜臀av| 亚洲精品久久久久久婷婷小说 | 国产日韩欧美在线精品| 国产一区二区三区av在线| 99国产精品一区二区蜜桃av| 乱码一卡2卡4卡精品| 两个人的视频大全免费| 校园人妻丝袜中文字幕| 国产大屁股一区二区在线视频| 身体一侧抽搐| 观看美女的网站| 99热这里只有精品一区| 22中文网久久字幕| 一区二区三区四区激情视频| 国产在视频线在精品| 亚洲中文字幕一区二区三区有码在线看| 国产国拍精品亚洲av在线观看| 51国产日韩欧美| 日韩国内少妇激情av| 搡女人真爽免费视频火全软件| 一区二区三区乱码不卡18| 欧美一级a爱片免费观看看| 国内少妇人妻偷人精品xxx网站| 亚洲最大成人av| 18禁动态无遮挡网站| 白带黄色成豆腐渣| 99在线人妻在线中文字幕| 久久精品久久久久久久性| 黄色一级大片看看| 99在线人妻在线中文字幕| 日本黄色片子视频| 在线天堂最新版资源| 成年女人永久免费观看视频| 如何舔出高潮| 男的添女的下面高潮视频| 日本猛色少妇xxxxx猛交久久| 韩国av在线不卡| 白带黄色成豆腐渣| 亚洲欧洲日产国产| 国产乱人偷精品视频| 1000部很黄的大片| 日本免费a在线| 亚洲美女搞黄在线观看| 蜜臀久久99精品久久宅男| 岛国在线免费视频观看| 一级毛片aaaaaa免费看小| 亚洲成人精品中文字幕电影| 国产精品嫩草影院av在线观看| 国产在线男女| 亚洲怡红院男人天堂| 午夜精品在线福利| 日韩成人伦理影院| 精品午夜福利在线看| 三级国产精品片| 亚洲经典国产精华液单| 成人漫画全彩无遮挡| 亚洲一区高清亚洲精品| 国产一级毛片七仙女欲春2| 美女黄网站色视频| 国产精华一区二区三区| 老师上课跳d突然被开到最大视频| 亚洲丝袜综合中文字幕| 熟女电影av网| 国产一区二区在线观看日韩| 超碰av人人做人人爽久久| 人妻少妇偷人精品九色| 最近最新中文字幕大全电影3| 亚洲经典国产精华液单| 国产伦精品一区二区三区视频9| 青春草亚洲视频在线观看| 亚洲国产精品成人久久小说| 观看免费一级毛片| h日本视频在线播放| 成人三级黄色视频| 嫩草影院精品99| 99热6这里只有精品| 黄色欧美视频在线观看| 精品国产露脸久久av麻豆 | www日本黄色视频网| 少妇裸体淫交视频免费看高清| 国产精品久久久久久久久免| 91午夜精品亚洲一区二区三区| 色播亚洲综合网| 欧美一级a爱片免费观看看| 欧美日韩在线观看h| 亚洲国产最新在线播放| 亚洲,欧美,日韩| 国产精品国产三级国产av玫瑰| 国产视频内射| 精品一区二区三区人妻视频| 亚洲五月天丁香| 国产探花在线观看一区二区| 伊人久久精品亚洲午夜| 欧美精品国产亚洲| 亚洲久久久久久中文字幕| 亚洲av日韩在线播放| 观看美女的网站| 成人美女网站在线观看视频| 最新中文字幕久久久久| 国产亚洲一区二区精品| 午夜福利高清视频| 丰满少妇做爰视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区在线av高清观看| 成年av动漫网址| 一区二区三区高清视频在线| 午夜福利在线观看吧| 国产精品人妻久久久久久| 99热网站在线观看| 99久国产av精品| 午夜福利高清视频| 日日摸夜夜添夜夜爱| 免费人成在线观看视频色| 欧美bdsm另类| 99热这里只有是精品在线观看| 日本免费在线观看一区| 嫩草影院入口| 看黄色毛片网站| 久久久成人免费电影| 青春草国产在线视频| 小说图片视频综合网站| 国产黄片视频在线免费观看| 欧美变态另类bdsm刘玥| 最近最新中文字幕免费大全7| 桃色一区二区三区在线观看| 久久久久久九九精品二区国产| 嫩草影院入口| 亚洲国产精品合色在线| 成人漫画全彩无遮挡| 国产日韩欧美在线精品| 99热网站在线观看| 免费黄网站久久成人精品| 久久久久久九九精品二区国产| 亚洲婷婷狠狠爱综合网| 天天一区二区日本电影三级| 亚洲最大成人av| 男女那种视频在线观看| 亚洲国产精品专区欧美| 91精品国产九色| 国产女主播在线喷水免费视频网站 | 两个人视频免费观看高清| 可以在线观看毛片的网站| 久久婷婷人人爽人人干人人爱| 99热6这里只有精品| 精品熟女少妇av免费看| 日韩av在线大香蕉| 国产精品99久久久久久久久| 狂野欧美激情性xxxx在线观看| 久久久a久久爽久久v久久| 欧美日韩精品成人综合77777| 简卡轻食公司| 亚洲中文字幕一区二区三区有码在线看| 一二三四中文在线观看免费高清| 熟妇人妻久久中文字幕3abv| 九色成人免费人妻av| 精品人妻一区二区三区麻豆| 男人舔女人下体高潮全视频| 如何舔出高潮| 2021天堂中文幕一二区在线观| 好男人在线观看高清免费视频| 网址你懂的国产日韩在线| or卡值多少钱| 小蜜桃在线观看免费完整版高清| 久久精品影院6| 精品99又大又爽又粗少妇毛片| 国产av在哪里看| 午夜精品在线福利| 国产高清不卡午夜福利| 免费观看性生交大片5| 日韩一区二区视频免费看| 日日撸夜夜添| 中文字幕人妻熟人妻熟丝袜美| 五月伊人婷婷丁香| 美女高潮的动态| 中文字幕av成人在线电影| 国产亚洲av嫩草精品影院| 看免费成人av毛片| 综合色丁香网| 久久精品国产亚洲av涩爱| 欧美性猛交黑人性爽| 天堂网av新在线| 亚洲怡红院男人天堂| 嘟嘟电影网在线观看| 日韩av在线免费看完整版不卡| 男人和女人高潮做爰伦理| 欧美日韩综合久久久久久| 久久99热这里只频精品6学生 | 亚洲18禁久久av| 99久久成人亚洲精品观看| 色噜噜av男人的天堂激情| 干丝袜人妻中文字幕| 久久久亚洲精品成人影院| 国产亚洲精品久久久com| av福利片在线观看| videossex国产| 国产精品乱码一区二三区的特点| 99九九线精品视频在线观看视频| 人人妻人人看人人澡| 九色成人免费人妻av| 一级爰片在线观看| 国产精品一区二区三区四区久久| 国产日韩欧美在线精品| 国产男人的电影天堂91| 99热这里只有是精品50| av又黄又爽大尺度在线免费看 | eeuss影院久久| 成人鲁丝片一二三区免费| 日韩av不卡免费在线播放| 国产淫片久久久久久久久| 精品国产露脸久久av麻豆 | 亚洲欧美一区二区三区国产| 男的添女的下面高潮视频| 男人和女人高潮做爰伦理| 插逼视频在线观看| 美女cb高潮喷水在线观看| 精品不卡国产一区二区三区| 99国产精品一区二区蜜桃av| 国产精品.久久久| 日本免费a在线| 色5月婷婷丁香| 亚洲av福利一区| 极品教师在线视频| av福利片在线观看| 亚洲欧洲国产日韩| av卡一久久| 免费人成在线观看视频色| 国产极品天堂在线| 一区二区三区四区激情视频| 亚洲在线观看片| 欧美日本亚洲视频在线播放| 一本久久精品| 国产精品人妻久久久久久| 日韩制服骚丝袜av| 久99久视频精品免费| 亚洲av二区三区四区| 国产亚洲午夜精品一区二区久久 | 五月玫瑰六月丁香| 中文字幕制服av| 日韩人妻高清精品专区| 成人午夜高清在线视频| 久久精品国产亚洲av涩爱| 我要搜黄色片| 精品少妇黑人巨大在线播放 | 国产精品一区二区在线观看99 | 欧美xxxx性猛交bbbb| 22中文网久久字幕| 非洲黑人性xxxx精品又粗又长| 亚洲精品亚洲一区二区| 国产中年淑女户外野战色| 亚洲av中文av极速乱| 精品人妻熟女av久视频| 国产女主播在线喷水免费视频网站 | 最近2019中文字幕mv第一页| 汤姆久久久久久久影院中文字幕 | 国产精品永久免费网站| 色综合亚洲欧美另类图片| 欧美高清性xxxxhd video| 国产精品久久久久久精品电影小说 | 欧美性猛交黑人性爽| 国产在线男女| 免费看美女性在线毛片视频| 看免费成人av毛片| 欧美xxxx黑人xx丫x性爽| 人人妻人人澡人人爽人人夜夜 | 亚洲人与动物交配视频| 国产精品国产三级专区第一集| 人妻系列 视频| 欧美日本亚洲视频在线播放| 亚洲欧美精品综合久久99| 欧美日本亚洲视频在线播放| 女的被弄到高潮叫床怎么办| 久久精品夜夜夜夜夜久久蜜豆| av线在线观看网站| 国产精品,欧美在线| 日本一本二区三区精品| 亚洲国产欧美人成| 成人午夜精彩视频在线观看| 国产一区二区在线av高清观看| 精品久久久久久久人妻蜜臀av| or卡值多少钱| 日韩,欧美,国产一区二区三区 | 人妻夜夜爽99麻豆av| 国产人妻一区二区三区在| 男女啪啪激烈高潮av片| 免费看av在线观看网站| 精品国产三级普通话版| 国产高清三级在线| 国产成年人精品一区二区| 国产免费又黄又爽又色| 日本猛色少妇xxxxx猛交久久| 国产伦精品一区二区三区四那| 一级毛片电影观看 | 长腿黑丝高跟| 99久久精品一区二区三区| 淫秽高清视频在线观看| 国产免费男女视频| 毛片一级片免费看久久久久| 色网站视频免费| 午夜福利视频1000在线观看| 午夜福利高清视频| 人妻系列 视频| 国产真实伦视频高清在线观看| 极品教师在线视频| 日韩制服骚丝袜av| 久久久精品大字幕| 色哟哟·www| 精品久久久久久久人妻蜜臀av| 中文字幕熟女人妻在线| 日韩,欧美,国产一区二区三区 | ponron亚洲| 国产黄色视频一区二区在线观看 | 国产精品永久免费网站| 国产成年人精品一区二区| 久久精品熟女亚洲av麻豆精品 | 舔av片在线| 搞女人的毛片| 午夜福利在线观看吧| 精品熟女少妇av免费看| 狂野欧美白嫩少妇大欣赏| 2021天堂中文幕一二区在线观| 嘟嘟电影网在线观看| 久久亚洲精品不卡| 国产精品麻豆人妻色哟哟久久 | 国产精品一区www在线观看| 天天躁夜夜躁狠狠久久av| 麻豆精品久久久久久蜜桃| 国产精品无大码| 如何舔出高潮| 国产私拍福利视频在线观看| 成人美女网站在线观看视频| 嫩草影院精品99| 中文字幕人妻熟人妻熟丝袜美| 少妇被粗大猛烈的视频| 蜜桃久久精品国产亚洲av| 我的女老师完整版在线观看| 久久久久久久亚洲中文字幕| www.色视频.com| 又黄又爽又刺激的免费视频.| 国产成人精品婷婷| 看片在线看免费视频| 国产精品精品国产色婷婷| 亚洲国产精品合色在线| videossex国产| 精品熟女少妇av免费看| 国产黄色视频一区二区在线观看 | 中文在线观看免费www的网站| 精品午夜福利在线看| 午夜免费男女啪啪视频观看| 偷拍熟女少妇极品色| 国产精品.久久久| 久久久久九九精品影院| 日韩一本色道免费dvd| 少妇熟女欧美另类| 男人舔女人下体高潮全视频| 久久久久久久亚洲中文字幕| 综合色av麻豆| 国内精品一区二区在线观看| 午夜福利视频1000在线观看| 波多野结衣高清无吗| 中文字幕制服av| 日韩欧美三级三区| 国产精品人妻久久久久久| 日日摸夜夜添夜夜爱| 亚洲第一区二区三区不卡| 一级黄片播放器| 男女下面进入的视频免费午夜| 一级爰片在线观看| 最近视频中文字幕2019在线8| 亚洲欧美日韩东京热| 午夜福利视频1000在线观看| 久久久午夜欧美精品| 中文字幕免费在线视频6| 国产精品人妻久久久影院| 国产精品人妻久久久久久| 亚洲国产精品成人综合色| 夜夜看夜夜爽夜夜摸| 精品无人区乱码1区二区| 久久精品夜夜夜夜夜久久蜜豆| 国产免费又黄又爽又色| 超碰av人人做人人爽久久| 成人亚洲欧美一区二区av| 久久久国产成人免费| 精品人妻视频免费看| 老女人水多毛片| 97在线视频观看| 搞女人的毛片| 午夜福利视频1000在线观看| 免费看a级黄色片| 国产av一区在线观看免费| 国产精品日韩av在线免费观看| 国产精品99久久久久久久久| 亚洲综合精品二区| 神马国产精品三级电影在线观看| 欧美成人午夜免费资源| av又黄又爽大尺度在线免费看 | 久久精品综合一区二区三区| 久久鲁丝午夜福利片| 午夜免费激情av| 久久久久久九九精品二区国产| 亚洲欧美精品综合久久99| 精品午夜福利在线看| 久久久久久久久久成人| 99久久精品国产国产毛片| 欧美日本亚洲视频在线播放| 午夜久久久久精精品| 狂野欧美激情性xxxx在线观看| h日本视频在线播放| 国产精品美女特级片免费视频播放器| 色噜噜av男人的天堂激情| 啦啦啦韩国在线观看视频| 国产黄片视频在线免费观看| 亚洲av日韩在线播放| 91精品国产九色| 欧美精品国产亚洲| 久久鲁丝午夜福利片| 久久久国产成人免费| 精华霜和精华液先用哪个| 啦啦啦观看免费观看视频高清| 免费观看性生交大片5| 三级国产精品欧美在线观看| 99久国产av精品| 亚洲国产欧洲综合997久久,| 久久鲁丝午夜福利片| 国产淫片久久久久久久久| eeuss影院久久| 国产欧美日韩精品一区二区| 国产成人a∨麻豆精品| 国产一级毛片七仙女欲春2| 免费看日本二区| 欧美精品国产亚洲| 岛国毛片在线播放| 亚洲欧美一区二区三区国产| 日韩制服骚丝袜av| 在线播放无遮挡| 最近手机中文字幕大全| 精品人妻视频免费看| 国产色婷婷99| 狠狠狠狠99中文字幕| 三级国产精品欧美在线观看| 欧美极品一区二区三区四区| 97在线视频观看| 在线a可以看的网站| 桃色一区二区三区在线观看| av又黄又爽大尺度在线免费看 | 欧美性猛交╳xxx乱大交人| 国产精品久久久久久久电影| 嘟嘟电影网在线观看| 国产黄a三级三级三级人| 国语自产精品视频在线第100页| 精品无人区乱码1区二区| 最近最新中文字幕免费大全7| 搞女人的毛片| 成人性生交大片免费视频hd| 久久精品夜夜夜夜夜久久蜜豆| 看黄色毛片网站| 久久久久九九精品影院| 亚洲成人中文字幕在线播放| 啦啦啦韩国在线观看视频| 性插视频无遮挡在线免费观看| 国产极品天堂在线| 亚洲av成人精品一区久久| 亚洲欧美日韩卡通动漫| 男人舔奶头视频| 91精品伊人久久大香线蕉| 波多野结衣高清无吗| av在线老鸭窝| 99国产精品一区二区蜜桃av| 午夜日本视频在线| 又爽又黄a免费视频| 青青草视频在线视频观看| 国产女主播在线喷水免费视频网站 | 国产日韩欧美在线精品| 亚洲aⅴ乱码一区二区在线播放| 国产三级中文精品| 日韩欧美国产在线观看| 精品免费久久久久久久清纯| 日韩人妻高清精品专区| 老师上课跳d突然被开到最大视频| 欧美一区二区亚洲| 亚洲国产高清在线一区二区三| 18禁动态无遮挡网站| 国产亚洲最大av| 欧美三级亚洲精品| 久久草成人影院| 欧美一区二区国产精品久久精品| 全区人妻精品视频| 久久草成人影院| 成人特级av手机在线观看| 日本免费在线观看一区| 午夜精品国产一区二区电影 | 黄片无遮挡物在线观看| 成人一区二区视频在线观看| 亚洲国产最新在线播放| 亚洲精品成人久久久久久| 天堂av国产一区二区熟女人妻| 亚洲精品国产av成人精品|