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

    基于小波變換的高邊坡變形監(jiān)測數(shù)據(jù)誤差處理方法

    2021-01-05 10:57:34雷孟飛周俊華湯金毅
    導(dǎo)航定位學(xué)報 2020年6期
    關(guān)鍵詞:監(jiān)測數(shù)據(jù)小波監(jiān)測點(diǎn)

    雷孟飛,周俊華,湯金毅

    (湖南聯(lián)智科技股份有限公司,長沙 410073)

    0 引言

    北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou navigation satellite system, BDS)是我國自主研發(fā)的全天時、全天候提供高精度位置信息的導(dǎo)航定位系統(tǒng)。隨著BDS的逐步完善,BDS 定位精度不斷提高,BDS 高精度定位技術(shù)也逐漸被應(yīng)用到各種工程領(lǐng)域施工期和維護(hù)期的變形監(jiān)測中[1]。利用BDS 進(jìn)行變形監(jiān)測,可以克服人工測量所存在的受自然條件限制嚴(yán)重、不能實現(xiàn)自動化監(jiān)測、監(jiān)測成本高等弊端,因此在變形監(jiān)測中的應(yīng)用越來越廣泛。

    BDS 自動化監(jiān)測系統(tǒng)采用靜態(tài)相對定位的算法進(jìn)行高精度解算,由于基準(zhǔn)站和監(jiān)測站的距離相對較近,解算結(jié)果水平方向的精度通常能夠達(dá)到±3 mm,高程方向能達(dá)到±5 mm。監(jiān)測結(jié)果中的大多誤差已經(jīng)通過差分算法消除,但是對于一些周期性的誤差,比如多路徑誤差、周跳以及設(shè)備噪聲造成的高頻誤差和異常跳動則無法剔除干凈,特別對于一些高邊坡監(jiān)測項目,通常信號遮擋嚴(yán)重、衛(wèi)星星座分布不均,周期性誤差和高頻誤差明顯,并不能得到有效剔除[2]。這些周期性誤差和高頻誤差為邊坡監(jiān)測帶來很大困難,技術(shù)人員很難從各項誤差中提取有用信號來判斷邊坡是否發(fā)生形變以及是否需要形變預(yù)警,因此有必要對BDS變形監(jiān)測的數(shù)據(jù)進(jìn)行降噪處理,獲得監(jiān)測點(diǎn)的真實位移坐標(biāo)時間序列。

    小波變換是1 種新的變換分析方法,它克服了傅里葉變換窗口大小不能隨頻率變化等缺點(diǎn),提供隨頻率改變的“時間-頻率”窗口,可以探測信號中的瞬態(tài)成分及頻率成分[3]。在BDS 高邊坡變形監(jiān)測過程中,可以利用小波變換將監(jiān)測結(jié)果中的高頻噪聲信號分離出來并進(jìn)行去除,提取有效的形變信息,提高監(jiān)測精度和穩(wěn)定性。國內(nèi)外學(xué)者已經(jīng)將小波變換較多地運(yùn)用在全球定位系統(tǒng)(global positioning system, GPS)坐標(biāo)時間序列分析中。文獻(xiàn)[4]提出用小波變換方法提取 GPS 觀測數(shù)據(jù)中的多路徑信號,對觀測結(jié)果進(jìn)行修正;文獻(xiàn)[5]利用小波分析對全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system, GNSS)高頻低頻多路徑誤差進(jìn)行濾波降噪,效果良好;文獻(xiàn)[6]提出了1 種基于小波濾波與主成分分析相結(jié)合的GPS 去噪方法,構(gòu)建GPS多路徑誤差的誤差模型,并對GPS 多路徑誤差進(jìn)行修正;文獻(xiàn)[7]采用小波變換,提取了中國區(qū)域8 個國際GNSS 服務(wù)組織(International GNSS Service,IGS)站的坐標(biāo)時間序列中的非線性運(yùn)動信息,對季節(jié)性周期變化進(jìn)行了分析研究。上述研究中所采用的GPS 坐標(biāo)時序數(shù)據(jù)大多為IGS 站點(diǎn)數(shù)據(jù),數(shù)據(jù)質(zhì)量較高,沒有考慮實際應(yīng)用中衛(wèi)星被遮擋、其他信號干擾的情況,因此小波變換在BDS 高邊坡監(jiān)測中的應(yīng)用效果并不明確,并且上述研究中并沒有針對性地對坐標(biāo)時間序列中不同頻率的誤差項進(jìn)行單獨(dú)改正,監(jiān)測結(jié)果容易丟失有用的細(xì)節(jié)部分。

    針對上述問題,本文以廣西某高速公路高邊坡為實驗測區(qū),首先利用小波變換算法對監(jiān)測結(jié)果進(jìn)行分解,利用低頻層恢復(fù)測試數(shù)據(jù),剔除監(jiān)測數(shù)據(jù)中的粗差數(shù)據(jù),然后對監(jiān)測數(shù)據(jù)進(jìn)行2 次小波變換,并對變換后的各層系數(shù)分別進(jìn)行改正,最后恢復(fù)監(jiān)測數(shù)據(jù),分析誤差剔除效果。

    1 高邊坡變形監(jiān)測數(shù)據(jù)誤差處理方法

    1.1 小波變換

    小波變換主要分為連續(xù)小波變換和離散小波變換 2 類[8]。本文采用較常用的離散小波變換(discrete wavelet transform, DWT)進(jìn)行處理,離散小波變換基本公式[9]為

    式中:φ ( t)為小波母函數(shù),對其進(jìn)行平移和伸縮得到小波基函數(shù)集 φa,b(t );a 為尺度因子;b 為平移因子;R 為實數(shù)集。通常取

    式中:Z 為整數(shù)集; a0為尺度因子擴(kuò)展步長;0b 為平移因子擴(kuò)展步長;m 為尺度因子擴(kuò)展系數(shù);n 為平移因子擴(kuò)展系數(shù)。

    將式(2)代入式(1)可得[10]

    離散小波變換公式為

    式中: fDWT( a , b )為連續(xù)小波變換系數(shù);*表示共軛[11]。

    1.2 變形監(jiān)測數(shù)據(jù)誤差處理流程

    基于小波變換的高邊坡變形監(jiān)測數(shù)據(jù)誤差處理方法其基本思想是:首先利用小波變換對監(jiān)測結(jié)果的粗差進(jìn)行剔除;然后對處理后的時間序列進(jìn)行2 次小波變換;最后將變換后的系數(shù)分為高頻、中頻、低頻3 類,分別對每層細(xì)節(jié)系數(shù)采用不同的濾波算法進(jìn)行濾波,利用各層濾波后的細(xì)節(jié)系數(shù)恢復(fù)監(jiān)測時間序列。其方法流程如圖1 所示,其中各層系數(shù)改正的處理流程如圖2 所示。

    圖1 數(shù)據(jù)處理流程

    1)數(shù)據(jù)預(yù)處理。選擇監(jiān)測點(diǎn)一段時間的監(jiān)測數(shù)據(jù)進(jìn)行處理,假設(shè)監(jiān)測數(shù)據(jù)采樣頻率為s,監(jiān)測數(shù)據(jù)時間序列為Y,數(shù)據(jù)預(yù)處理采用3 次樣條插值算法對缺失數(shù)據(jù)進(jìn)行插值;

    2)小波變換。采用小波變換對監(jiān)測點(diǎn)數(shù)據(jù)進(jìn)行分解,分解后得到細(xì)節(jié)信號(監(jiān)測結(jié)果的細(xì)節(jié)信息)以及趨勢信號(包含監(jiān)測點(diǎn)的位移趨勢);

    3)分解系數(shù)分類。①分別計算第2 步的分層結(jié)果的頻帶范圍,計算公式[12]為

    圖2 各層系數(shù)改正流程

    式中:Y ' 為監(jiān)測結(jié)果的采樣頻率,通用的單位采樣頻率為1 s/次,而本文中變形監(jiān)測的采樣頻率較低,因此,本文中定義的單位采樣頻率為1 h/次;j 為分解后的第j 層;式(5)為細(xì)節(jié)系數(shù)的頻帶范圍,式(6)為最后一層分解后的近似系數(shù)的頻帶范圍。②根據(jù)不同頻段的閾值,將上一步驟計算得到的各層的頻帶范圍進(jìn)行區(qū)分,其中頻段的閾值為各分層系數(shù)高頻、中頻和低頻的區(qū)分值,高頻信號定義為頻率大于24 h/次,中頻信號約為24 h/次即周期約為24 h,低頻信號頻率小于24 h/次,其具體數(shù)值根據(jù)監(jiān)測數(shù)據(jù)的長度和采樣頻率確定。

    4)粗差識別修復(fù)。①利用第3 步得到的低頻層和中頻系數(shù)采用小波逆變換恢復(fù)監(jiān)測數(shù)據(jù),得到監(jiān)測數(shù)據(jù) Y ';②利用3 倍中誤差原則識別粗差數(shù)據(jù),具體為:|Yi- Yi'| > 3σ的數(shù)據(jù)為粗差數(shù)據(jù),其中 Yi為第i 個原始監(jiān)測數(shù)據(jù), Yi'為第i 個低頻系數(shù)恢復(fù)后的原始數(shù)據(jù);③利用臨近插值法對識別出的粗差數(shù)據(jù)進(jìn)行修復(fù),得到粗差修復(fù)后的監(jiān)測數(shù)據(jù) Y'' 。

    5)高頻系數(shù)濾波。將步第4 步2 次小波變換后得到的各層變換系數(shù)按照第3 步的方法進(jìn)行分類,同樣分為高頻和低頻2 部分,高頻系數(shù)部分采用自適應(yīng)卡爾曼濾波算法對高頻系數(shù)進(jìn)行濾波。

    6)低頻系數(shù)濾波。低頻系數(shù)部分主要包含了監(jiān)測點(diǎn)的趨勢信息,因此低頻系數(shù)部分采用滑動平均算法進(jìn)行簡單濾波,去除部分噪點(diǎn),最大限度保留監(jiān)測點(diǎn)的趨勢信息,算法公式為

    7)中頻系數(shù)改正。中頻系數(shù)改正采用中誤差原則進(jìn)行改正,首先計算中頻系數(shù)的中誤差Lσ ,然后按如下公式進(jìn)行改正,即

    式中:Mi為第i 個中頻系數(shù);Lσ 為中頻系數(shù)中誤差;為第i 個中頻系數(shù)改正值。

    8)監(jiān)測數(shù)據(jù)恢復(fù)。利用處理后各層的系數(shù)恢復(fù)監(jiān)測數(shù)據(jù)。

    2 實驗與結(jié)果分析

    2.1 數(shù)據(jù)來源

    本文采用的監(jiān)測數(shù)據(jù)來自廣西陽朔至鹿寨高速公路邊坡2 號監(jiān)測點(diǎn)(如圖3 所示),該邊坡高約60 m,五級邊坡,植被茂密,邊坡后方為陡坡向上。數(shù)據(jù)采集設(shè)備為本單位自己研發(fā)的BDS 高精度接收機(jī),解算軟件亦為本單位自己研發(fā)的高精度變形監(jiān)測后處理軟件,數(shù)據(jù)解算間隔為1 h,數(shù)據(jù)時間為2020-02-01—2020-05-19,另外在數(shù)據(jù)中加入3 個隨機(jī)粗差來檢驗粗差剔除效果。

    圖3 監(jiān)測點(diǎn)概況

    2.2 處理結(jié)果

    1)粗差修復(fù)。數(shù)據(jù)的小波變換函數(shù)選擇“db08”,分解層數(shù)為6 層,a6 為第6 層分解后的近似系數(shù),計算得到的各層頻帶范圍如表1 所示。

    表1 頻帶范圍表單位:Hz

    計算得到上文所述高頻、中頻和低頻的閾值為0.041 7 Hz,高頻層為d1、d2、d3,中頻層為d4,低頻層為d5、d6。利用低頻層進(jìn)行小波逆變換,恢復(fù)監(jiān)測數(shù)據(jù),并進(jìn)行粗差識別和修復(fù),修復(fù)效果如圖4所示。

    圖4 粗差修復(fù)效果

    從圖4 可以看出,3 處粗差數(shù)據(jù)都能夠得到較好修復(fù),修復(fù)效果良好。

    2)細(xì)節(jié)系數(shù)濾波。將上一步驟經(jīng)過粗差修復(fù)后的數(shù)據(jù)進(jìn)行2 次小波變換,然后對各項系數(shù)進(jìn)行改正,改正效果如圖5 所示。

    圖5 各層細(xì)節(jié)系數(shù)濾波效果

    從圖5 各層細(xì)節(jié)系數(shù)濾波效果圖中可以看出:d1、d2、d3、d4、d5及d6各層的細(xì)節(jié)系數(shù)都有比較大的改正,其中d1、d2、d3及d4的改正效果最好,系數(shù)的頻率和最大、最小值都有一定的減??;d4周期信號經(jīng)過改正以后,振幅明顯減小,周期性更加明顯;d5和d6層改正前后細(xì)節(jié)信息差別不是很大,主要原因是為了保留監(jiān)測點(diǎn)位移信息僅采用了滑動平均方法進(jìn)行濾波,能夠看到濾波前后d5和d62 層的整體趨勢都沒有改變,但是改正后的曲線更加平滑,細(xì)節(jié)抖動去除明顯。

    3)數(shù)據(jù)恢復(fù)。利用濾波后的各層系數(shù)進(jìn)行小波逆變換,變換后的效果如圖6 所示。

    圖6 誤差削弱效果

    從圖6 的誤差削弱效果中可以看出:處理后的觀測數(shù)據(jù)時間序列的平穩(wěn)性有很大提高,其中以天為周期的周期性位移信號有了較大改善,數(shù)據(jù)振幅更小;處理以后,噪聲信息明顯減少,可以清晰地看出點(diǎn)位的位移趨勢,有利于更精確地分析監(jiān)測站的位移趨勢;在3 月30 日以后監(jiān)測點(diǎn)發(fā)生20 mm 左右位移,經(jīng)過小波變換處理以后監(jiān)測點(diǎn)的位移曲線更加明顯,位移量與處理前保持一致,位移開始時間并沒有延遲;另外處理后的監(jiān)測點(diǎn)數(shù)據(jù)3 月20—30 日之間的小幅度位移也得到了有效保留;統(tǒng)計了3 月20 日以前原始數(shù)據(jù)數(shù)據(jù)的中誤差為3.58 mm,處理后數(shù)據(jù)的中誤差降為0.92 mm;從上述分析可以看出本文小波變換方法處理后的監(jiān)測數(shù)據(jù)時間序列的平穩(wěn)性和數(shù)據(jù)質(zhì)量都有明顯提高。

    3 結(jié)束語

    本文通過廣西陽朔至鹿寨高速公路高邊坡進(jìn)行小波變換數(shù)據(jù)處理,可以得到如下結(jié)論:

    1)利用本文的小波變換方法可以有效去除監(jiān)測結(jié)果中的粗差數(shù)據(jù),通過3 次樣條插值算法,能夠很好的對粗差數(shù)據(jù)進(jìn)行修復(fù);

    2)本文利用2 次小波變化對剔除粗差后的監(jiān)測數(shù)據(jù)進(jìn)行處理,能夠有效去除監(jiān)測結(jié)果時間序列中的高頻誤差信號和以天為周期的周期性信號,提高數(shù)據(jù)的平穩(wěn)性;

    3)經(jīng)過數(shù)據(jù)處理以后,處理結(jié)果中能夠有效保留監(jiān)測點(diǎn)的整體位移趨勢,位移量的大小并沒有損失,并且小位移量的細(xì)節(jié)部分也能夠有效保留,為技術(shù)人員判斷邊坡的形變狀態(tài)提供精確的數(shù)據(jù)支持,可以降低滑坡災(zāi)害,減少生命和財產(chǎn)損失。

    猜你喜歡
    監(jiān)測數(shù)據(jù)小波監(jiān)測點(diǎn)
    天津南港LNG接收站沉降監(jiān)測點(diǎn)位布設(shè)
    煤氣與熱力(2022年4期)2022-05-23 12:44:56
    撫河流域綜合治理監(jiān)測布局優(yōu)化
    構(gòu)造Daubechies小波的一些注記
    全站儀極坐標(biāo)法監(jiān)測點(diǎn)穩(wěn)定性分析方法研究
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    GSM-R接口監(jiān)測數(shù)據(jù)精確地理化方法及應(yīng)用
    基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    我省舉辦家畜血吸蟲病監(jiān)測點(diǎn)培訓(xùn)班
    GPS異常監(jiān)測數(shù)據(jù)的關(guān)聯(lián)負(fù)選擇分步識別算法
    基于小波函數(shù)對GNSS監(jiān)測數(shù)據(jù)降噪的應(yīng)用研究
    脱女人内裤的视频| 久久人妻福利社区极品人妻图片| 一本久久中文字幕| 身体一侧抽搐| 久久午夜综合久久蜜桃| 宅男免费午夜| 国产精品一区二区精品视频观看| 91麻豆av在线| 国产精品九九99| 亚洲天堂国产精品一区在线| 一区二区三区高清视频在线| 制服丝袜大香蕉在线| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av第一区精品v没综合| 欧美亚洲日本最大视频资源| x7x7x7水蜜桃| 99精品欧美一区二区三区四区| 搡老熟女国产l中国老女人| 可以免费在线观看a视频的电影网站| 一个人免费在线观看的高清视频| 国产熟女xx| 99国产综合亚洲精品| 老司机午夜福利在线观看视频| 99国产精品一区二区三区| 国内毛片毛片毛片毛片毛片| 日本一本二区三区精品| 黄片小视频在线播放| 国产精品亚洲一级av第二区| 亚洲成人久久爱视频| 日韩欧美三级三区| 国产一区二区三区视频了| 十八禁人妻一区二区| 精品卡一卡二卡四卡免费| 狠狠狠狠99中文字幕| 少妇被粗大的猛进出69影院| 国产黄片美女视频| 国产v大片淫在线免费观看| 嫩草影院精品99| 91九色精品人成在线观看| 精品欧美国产一区二区三| 亚洲成人久久性| 这个男人来自地球电影免费观看| 久久亚洲精品不卡| 国产精品综合久久久久久久免费| 欧美日本视频| 色尼玛亚洲综合影院| av福利片在线| 叶爱在线成人免费视频播放| 巨乳人妻的诱惑在线观看| 国产亚洲精品久久久久久毛片| 脱女人内裤的视频| 国产精品免费一区二区三区在线| 免费一级毛片在线播放高清视频| 一边摸一边抽搐一进一小说| 亚洲国产精品999在线| 黑人欧美特级aaaaaa片| 欧美三级亚洲精品| 婷婷六月久久综合丁香| 夜夜躁狠狠躁天天躁| 夜夜爽天天搞| 日本免费一区二区三区高清不卡| 我的亚洲天堂| 久久久国产成人精品二区| 亚洲五月天丁香| 免费在线观看影片大全网站| 国产av一区在线观看免费| 午夜精品在线福利| 亚洲精品国产区一区二| 在线观看免费日韩欧美大片| 嫩草影视91久久| 欧美一级毛片孕妇| 一二三四社区在线视频社区8| 精品国内亚洲2022精品成人| 一区二区三区精品91| 丰满的人妻完整版| 丝袜在线中文字幕| 欧美最黄视频在线播放免费| 两性夫妻黄色片| 伊人久久大香线蕉亚洲五| 久久久久国产一级毛片高清牌| 午夜福利一区二区在线看| 亚洲精品av麻豆狂野| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久久久免费视频了| 在线观看午夜福利视频| 欧美激情极品国产一区二区三区| 男女视频在线观看网站免费 | 久久99热这里只有精品18| 777久久人妻少妇嫩草av网站| 首页视频小说图片口味搜索| 老鸭窝网址在线观看| 无遮挡黄片免费观看| 黑人巨大精品欧美一区二区mp4| 欧美乱码精品一区二区三区| 三级毛片av免费| 中亚洲国语对白在线视频| 国产免费av片在线观看野外av| 国产男靠女视频免费网站| 超碰成人久久| 亚洲第一av免费看| 一级毛片女人18水好多| 欧美日韩中文字幕国产精品一区二区三区| 手机成人av网站| 99国产综合亚洲精品| 男人舔女人的私密视频| 九色国产91popny在线| 精品欧美国产一区二区三| 国产精品日韩av在线免费观看| 人人妻,人人澡人人爽秒播| 久久久久免费精品人妻一区二区 | 中国美女看黄片| 亚洲第一av免费看| 国产乱人伦免费视频| 午夜免费鲁丝| 亚洲精品国产一区二区精华液| 亚洲一卡2卡3卡4卡5卡精品中文| 88av欧美| 国语自产精品视频在线第100页| 亚洲国产中文字幕在线视频| 青草久久国产| 午夜影院日韩av| 国产精品乱码一区二三区的特点| 久久伊人香网站| 久久久久国产一级毛片高清牌| 黄色视频不卡| 久久久久久免费高清国产稀缺| 国产久久久一区二区三区| 欧美性猛交╳xxx乱大交人| 十八禁网站免费在线| 亚洲国产高清在线一区二区三 | 草草在线视频免费看| 久久人人精品亚洲av| 91字幕亚洲| 中文亚洲av片在线观看爽| 可以在线观看毛片的网站| 中文在线观看免费www的网站 | 在线天堂中文资源库| 免费无遮挡裸体视频| 亚洲第一欧美日韩一区二区三区| www.999成人在线观看| 国产成人av激情在线播放| 色精品久久人妻99蜜桃| 国产成年人精品一区二区| 久久久久久久久中文| 精品无人区乱码1区二区| 丝袜美腿诱惑在线| 午夜福利欧美成人| 男人舔奶头视频| 变态另类成人亚洲欧美熟女| 色播在线永久视频| 国产精品野战在线观看| 一级毛片女人18水好多| 久久国产精品人妻蜜桃| 级片在线观看| 亚洲成人久久爱视频| 中文字幕另类日韩欧美亚洲嫩草| 丝袜在线中文字幕| 成人18禁在线播放| 精品免费久久久久久久清纯| 久热这里只有精品99| 麻豆久久精品国产亚洲av| 757午夜福利合集在线观看| 怎么达到女性高潮| 欧美日韩亚洲国产一区二区在线观看| 欧美成狂野欧美在线观看| 欧美日韩亚洲国产一区二区在线观看| 色播亚洲综合网| 性色av乱码一区二区三区2| 免费av毛片视频| 国产一区二区三区视频了| svipshipincom国产片| 最近在线观看免费完整版| 一级片免费观看大全| 免费人成视频x8x8入口观看| 免费观看人在逋| 国产精品精品国产色婷婷| 国产精品九九99| 国产一区在线观看成人免费| 亚洲一区二区三区不卡视频| 国内久久婷婷六月综合欲色啪| 美女 人体艺术 gogo| 亚洲色图 男人天堂 中文字幕| 少妇被粗大的猛进出69影院| 亚洲自拍偷在线| 村上凉子中文字幕在线| 久久久久亚洲av毛片大全| 亚洲av成人不卡在线观看播放网| 色哟哟哟哟哟哟| 麻豆久久精品国产亚洲av| 免费在线观看影片大全网站| 丝袜在线中文字幕| 91九色精品人成在线观看| 脱女人内裤的视频| 国产国语露脸激情在线看| 国产不卡一卡二| 老汉色∧v一级毛片| 午夜成年电影在线免费观看| 欧美在线一区亚洲| 亚洲自拍偷在线| 欧美乱码精品一区二区三区| 老司机午夜十八禁免费视频| 视频区欧美日本亚洲| 美女午夜性视频免费| 巨乳人妻的诱惑在线观看| 12—13女人毛片做爰片一| 啦啦啦观看免费观看视频高清| 成年人黄色毛片网站| 女性生殖器流出的白浆| 琪琪午夜伦伦电影理论片6080| 亚洲熟妇熟女久久| 国产欧美日韩精品亚洲av| 日韩欧美一区二区三区在线观看| 国产不卡一卡二| 久久精品国产亚洲av香蕉五月| 18禁国产床啪视频网站| 日韩欧美三级三区| 亚洲av成人av| 欧美在线黄色| 两个人免费观看高清视频| 国产精品一区二区三区四区久久 | 精品久久久久久成人av| 欧美丝袜亚洲另类 | 国产人伦9x9x在线观看| 丰满人妻熟妇乱又伦精品不卡| 丰满人妻熟妇乱又伦精品不卡| 高清在线国产一区| 免费在线观看日本一区| 一级作爱视频免费观看| 黄色视频,在线免费观看| 亚洲精品美女久久久久99蜜臀| 天堂√8在线中文| 亚洲无线在线观看| 国产亚洲精品久久久久5区| 亚洲自拍偷在线| 免费看a级黄色片| 久热爱精品视频在线9| 国产午夜福利久久久久久| 国产精品电影一区二区三区| 可以在线观看的亚洲视频| 国产精华一区二区三区| 国产一区二区三区在线臀色熟女| 男人舔女人下体高潮全视频| 国产不卡一卡二| www日本黄色视频网| 美女大奶头视频| 日韩精品免费视频一区二区三区| 免费av毛片视频| 在线天堂中文资源库| a级毛片a级免费在线| 欧美日韩中文字幕国产精品一区二区三区| 日本精品一区二区三区蜜桃| 亚洲成av片中文字幕在线观看| 亚洲av美国av| 19禁男女啪啪无遮挡网站| 一二三四社区在线视频社区8| 天天躁狠狠躁夜夜躁狠狠躁| 国产单亲对白刺激| 99精品久久久久人妻精品| 久久精品国产亚洲av香蕉五月| 国产男靠女视频免费网站| 人人妻人人澡欧美一区二区| av天堂在线播放| 亚洲av成人不卡在线观看播放网| 国产黄片美女视频| 麻豆成人av在线观看| 国产精品1区2区在线观看.| 国语自产精品视频在线第100页| 亚洲国产精品合色在线| 亚洲午夜理论影院| 久久伊人香网站| 国产麻豆成人av免费视频| 无人区码免费观看不卡| videosex国产| 男女之事视频高清在线观看| 国产一区二区激情短视频| 日韩大码丰满熟妇| 欧美成人免费av一区二区三区| 99国产精品99久久久久| 中文亚洲av片在线观看爽| www日本黄色视频网| 国产一区二区三区在线臀色熟女| 午夜福利视频1000在线观看| 制服诱惑二区| 美女 人体艺术 gogo| 狂野欧美激情性xxxx| 人人澡人人妻人| 老熟妇乱子伦视频在线观看| 午夜免费成人在线视频| 男人舔女人下体高潮全视频| 欧美午夜高清在线| 91av网站免费观看| 一夜夜www| 亚洲成av人片免费观看| 十分钟在线观看高清视频www| 久久久久久亚洲精品国产蜜桃av| 91麻豆av在线| 最近最新中文字幕大全电影3 | 18禁黄网站禁片午夜丰满| 人人妻人人澡人人看| a级毛片在线看网站| 日本 欧美在线| 免费女性裸体啪啪无遮挡网站| 亚洲久久久国产精品| 国产欧美日韩精品亚洲av| 一级a爱片免费观看的视频| 男人舔女人的私密视频| 国内少妇人妻偷人精品xxx网站 | 一边摸一边做爽爽视频免费| 一区二区三区高清视频在线| 亚洲精品久久成人aⅴ小说| 十八禁网站免费在线| 欧美在线一区亚洲| 免费在线观看黄色视频的| svipshipincom国产片| 黑人操中国人逼视频| 搞女人的毛片| 国产精品98久久久久久宅男小说| 在线观看www视频免费| 丝袜在线中文字幕| 麻豆成人av在线观看| 亚洲精品在线美女| 久久欧美精品欧美久久欧美| 91字幕亚洲| 久久久国产欧美日韩av| 巨乳人妻的诱惑在线观看| 一级a爱片免费观看的视频| 国产精品免费一区二区三区在线| 啦啦啦 在线观看视频| 国产黄色小视频在线观看| 国产亚洲av嫩草精品影院| 99久久国产精品久久久| 亚洲人成伊人成综合网2020| 十分钟在线观看高清视频www| 亚洲人成电影免费在线| 嫁个100分男人电影在线观看| 99国产精品一区二区蜜桃av| 两个人视频免费观看高清| 丝袜人妻中文字幕| 最新在线观看一区二区三区| 成人18禁高潮啪啪吃奶动态图| 国产伦在线观看视频一区| 久久久久久免费高清国产稀缺| 精品国产一区二区三区四区第35| 欧美日韩中文字幕国产精品一区二区三区| 亚洲国产欧美一区二区综合| 国产精品久久视频播放| 国产精品美女特级片免费视频播放器 | 大型黄色视频在线免费观看| 久久国产精品男人的天堂亚洲| 男女做爰动态图高潮gif福利片| 亚洲精品美女久久av网站| 亚洲aⅴ乱码一区二区在线播放 | 国产激情偷乱视频一区二区| 国产精品免费视频内射| 999久久久精品免费观看国产| 欧美+亚洲+日韩+国产| 日韩欧美三级三区| 亚洲aⅴ乱码一区二区在线播放 | 国产成+人综合+亚洲专区| 亚洲va日本ⅴa欧美va伊人久久| 丝袜在线中文字幕| 亚洲精品久久成人aⅴ小说| 男女做爰动态图高潮gif福利片| 国产精品一区二区三区四区久久 | 欧美亚洲日本最大视频资源| 欧美三级亚洲精品| 欧美日韩黄片免| 亚洲av五月六月丁香网| 岛国在线观看网站| 色老头精品视频在线观看| 伊人久久大香线蕉亚洲五| 亚洲自拍偷在线| 国产一区在线观看成人免费| 天天添夜夜摸| 国产精品久久视频播放| 18禁黄网站禁片午夜丰满| 夜夜爽天天搞| 精品欧美国产一区二区三| 91麻豆av在线| 精品久久久久久久人妻蜜臀av| 女同久久另类99精品国产91| 成人三级黄色视频| 亚洲,欧美精品.| 亚洲九九香蕉| 欧美性猛交黑人性爽| 欧美日韩乱码在线| 色精品久久人妻99蜜桃| 久久久久九九精品影院| 久9热在线精品视频| 黄色片一级片一级黄色片| 欧美日韩一级在线毛片| √禁漫天堂资源中文www| 国产一区在线观看成人免费| 免费在线观看完整版高清| 香蕉丝袜av| a在线观看视频网站| 香蕉国产在线看| 亚洲人成伊人成综合网2020| 黄片播放在线免费| a级毛片a级免费在线| 成人18禁在线播放| 午夜激情av网站| 国产成年人精品一区二区| av电影中文网址| 伊人久久大香线蕉亚洲五| 国产亚洲精品久久久久5区| 成人三级做爰电影| 久久香蕉精品热| 中文字幕人妻丝袜一区二区| 国产麻豆成人av免费视频| 亚洲国产日韩欧美精品在线观看 | 精品久久蜜臀av无| 女性被躁到高潮视频| 18禁观看日本| 久久久久久久午夜电影| 又大又爽又粗| 亚洲第一电影网av| 无限看片的www在线观看| 男女午夜视频在线观看| 88av欧美| 制服诱惑二区| 日本三级黄在线观看| 美女国产高潮福利片在线看| 男女午夜视频在线观看| 精品久久久久久久毛片微露脸| 两人在一起打扑克的视频| 熟妇人妻久久中文字幕3abv| 在线永久观看黄色视频| av福利片在线| 国产亚洲精品一区二区www| 久久婷婷人人爽人人干人人爱| 欧美日韩黄片免| ponron亚洲| 欧美激情极品国产一区二区三区| 亚洲人成伊人成综合网2020| 手机成人av网站| 亚洲国产欧洲综合997久久, | 亚洲人成网站高清观看| 欧美精品亚洲一区二区| 97人妻精品一区二区三区麻豆 | 好男人在线观看高清免费视频 | 免费在线观看完整版高清| 最近在线观看免费完整版| 日日爽夜夜爽网站| 日本一区二区免费在线视频| 国产精品九九99| 欧美中文日本在线观看视频| 在线观看免费视频日本深夜| 亚洲精品色激情综合| 欧美日韩亚洲综合一区二区三区_| 亚洲人成77777在线视频| 欧美亚洲日本最大视频资源| 亚洲精品在线观看二区| 免费一级毛片在线播放高清视频| 黄色成人免费大全| 国产高清激情床上av| 久久久久国内视频| 欧美一级a爱片免费观看看 | 免费观看精品视频网站| 国产精品综合久久久久久久免费| 最近最新中文字幕大全免费视频| 精品国产亚洲在线| 久久国产精品男人的天堂亚洲| 日本精品一区二区三区蜜桃| e午夜精品久久久久久久| 亚洲av成人一区二区三| 色哟哟哟哟哟哟| 国内揄拍国产精品人妻在线 | 久久亚洲真实| 国产人伦9x9x在线观看| 嫩草影视91久久| 国产精品久久视频播放| 脱女人内裤的视频| 国产三级在线视频| 欧美zozozo另类| 久久热在线av| 亚洲精品一区av在线观看| 亚洲精华国产精华精| 成人欧美大片| 久久精品夜夜夜夜夜久久蜜豆 | 悠悠久久av| 欧美亚洲日本最大视频资源| 国产精品久久久久久人妻精品电影| 亚洲第一欧美日韩一区二区三区| 亚洲av五月六月丁香网| 国产一卡二卡三卡精品| 午夜精品久久久久久毛片777| 日韩国内少妇激情av| 日韩欧美一区二区三区在线观看| 最近最新中文字幕大全免费视频| 国产精品 欧美亚洲| 欧美一区二区精品小视频在线| 欧美日韩亚洲综合一区二区三区_| 亚洲avbb在线观看| 色哟哟哟哟哟哟| 免费女性裸体啪啪无遮挡网站| 一本精品99久久精品77| 国产av在哪里看| 免费看十八禁软件| 免费一级毛片在线播放高清视频| 欧美在线黄色| 麻豆国产av国片精品| 国产在线精品亚洲第一网站| 亚洲av五月六月丁香网| 国产欧美日韩精品亚洲av| 欧美中文日本在线观看视频| 夜夜看夜夜爽夜夜摸| 欧美成人性av电影在线观看| 50天的宝宝边吃奶边哭怎么回事| 天天躁夜夜躁狠狠躁躁| 51午夜福利影视在线观看| 无人区码免费观看不卡| 一区二区三区国产精品乱码| 国内少妇人妻偷人精品xxx网站 | 美女 人体艺术 gogo| 国产亚洲精品久久久久久毛片| 国内毛片毛片毛片毛片毛片| 亚洲天堂国产精品一区在线| 亚洲国产精品sss在线观看| 久9热在线精品视频| 国产精品98久久久久久宅男小说| 嫩草影院精品99| 中文字幕最新亚洲高清| 在线视频色国产色| 又大又爽又粗| 亚洲av片天天在线观看| 国产又色又爽无遮挡免费看| av超薄肉色丝袜交足视频| 啦啦啦韩国在线观看视频| avwww免费| e午夜精品久久久久久久| 亚洲一区二区三区色噜噜| 精品免费久久久久久久清纯| 亚洲欧美日韩高清在线视频| 久久狼人影院| 国产熟女午夜一区二区三区| 白带黄色成豆腐渣| 日韩视频一区二区在线观看| 黄色丝袜av网址大全| 日韩成人在线观看一区二区三区| 久久香蕉激情| 一边摸一边做爽爽视频免费| 亚洲av第一区精品v没综合| 久久久久久久久中文| 99精品欧美一区二区三区四区| 黄色丝袜av网址大全| 丁香欧美五月| 在线看三级毛片| 亚洲av成人不卡在线观看播放网| 久久九九热精品免费| 国产亚洲av高清不卡| 日韩av在线大香蕉| 免费在线观看视频国产中文字幕亚洲| 亚洲精品粉嫩美女一区| 一本精品99久久精品77| 日韩精品免费视频一区二区三区| 99久久精品国产亚洲精品| 亚洲精品在线美女| 精品久久久久久久末码| 激情在线观看视频在线高清| 亚洲成人国产一区在线观看| 日韩一卡2卡3卡4卡2021年| 国产精品1区2区在线观看.| 久久久久久国产a免费观看| 欧美中文综合在线视频| 天天躁夜夜躁狠狠躁躁| 久久久久国内视频| 欧美精品啪啪一区二区三区| 日本精品一区二区三区蜜桃| 老司机午夜福利在线观看视频| 欧美一区二区精品小视频在线| 人成视频在线观看免费观看| 中文字幕久久专区| 黄色女人牲交| 亚洲免费av在线视频| 97碰自拍视频| 久久午夜亚洲精品久久| 亚洲熟女毛片儿| 日日干狠狠操夜夜爽| 韩国av一区二区三区四区| 91成年电影在线观看| 亚洲精品av麻豆狂野| 国产av不卡久久| 欧美久久黑人一区二区| 国产精品久久久人人做人人爽| 亚洲欧美精品综合久久99| 国产精品久久久久久亚洲av鲁大| 麻豆国产av国片精品| 在线观看www视频免费| 桃红色精品国产亚洲av| 黑人巨大精品欧美一区二区mp4| 国产成人精品久久二区二区91| 好看av亚洲va欧美ⅴa在| 丝袜美腿诱惑在线| 亚洲av中文字字幕乱码综合 | 在线看三级毛片| 久久欧美精品欧美久久欧美| 国产熟女午夜一区二区三区| 欧美成人午夜精品| 久久精品成人免费网站| 亚洲午夜精品一区,二区,三区| 极品教师在线免费播放| www.精华液| 午夜免费成人在线视频| 日本一区二区免费在线视频| 精品福利观看| 国产精品 国内视频| 欧美不卡视频在线免费观看 | 国产亚洲av高清不卡| 亚洲国产欧美一区二区综合| 亚洲av中文字字幕乱码综合 | 日韩有码中文字幕| 国内少妇人妻偷人精品xxx网站 | 国产视频一区二区在线看| 人成视频在线观看免费观看|