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

    雙樹復(fù)小波變換在川藏鐵路拉林段某隧道超前地質(zhì)預(yù)報(bào)中的應(yīng)用

    2021-12-23 07:28:40楊丹李偉魏永梁宋斌
    物探與化探 2021年6期
    關(guān)鍵詞:雙樹探地掌子面

    楊丹,李偉,魏永梁,宋斌

    (1.西南交通大學(xué) 地球科學(xué)與環(huán)境工程學(xué)院,四川 成都 611756;2.中國中鐵科學(xué)研究院有限公司,四川 成都 610032)

    0 引言

    為探明掌子面前方地質(zhì)條件,查清該隧道施工掌子面前方斷層、節(jié)理裂隙分布等工程地質(zhì)情況,保證施工安全,采用探地雷達(dá)進(jìn)行超前地質(zhì)預(yù)報(bào)。但由于掌子面不平整、隧道所在地層地質(zhì)條件極其復(fù)雜,探測條件較差,隧道內(nèi)施工環(huán)境復(fù)雜等原因,導(dǎo)致采集的數(shù)據(jù)中存在強(qiáng)干擾覆蓋有效信號的現(xiàn)象,嚴(yán)重影響數(shù)據(jù)采集質(zhì)量,干擾工程人員的判斷,因此需要進(jìn)行信號處理技術(shù)研究,以改善數(shù)據(jù)質(zhì)量,增強(qiáng)有效回波信號,消除雜波及多次回波干擾的影響,提高超前地質(zhì)預(yù)報(bào)的準(zhǔn)確性。

    傳統(tǒng)的探地雷達(dá)(GPR)去噪方法以傅里葉變換理論為基礎(chǔ),通過數(shù)學(xué)計(jì)算實(shí)現(xiàn)濾波,能夠在干擾信號和有效信號的頻譜范圍存在明顯分界時(shí)有效濾除指定頻率(干擾波頻譜)的信號。但探地雷達(dá)信號屬于非平穩(wěn)信號,用傳統(tǒng)的濾波方法去噪的同時(shí)也會拓寬波形,損失重要的信號信息,且無法同時(shí)體現(xiàn)信號頻率與時(shí)域的信息,難以滿足復(fù)雜地質(zhì)條件和復(fù)雜工程環(huán)境下的信號處理。近年來,專家、學(xué)者將小波變換理論引入探地雷達(dá)信號處理中,此舉在繼承和發(fā)展了短時(shí)傅里葉變換的局部化思想的同時(shí),又克服了窗口大小不隨頻率變換的缺點(diǎn),可利用小波變換的調(diào)焦功能和頻率域—時(shí)間域雙重局部性來壓制噪聲[2-6]。多數(shù)學(xué)者選用離散小波變換處理探地雷達(dá)信號數(shù)據(jù)[7-10],但離散小波變換存在平移敏感性、頻率混疊和方向性選擇少的缺點(diǎn),會影響數(shù)據(jù)處理的準(zhǔn)確性,而雙樹復(fù)小波變換(dual-tree complex wavelet transform,DTCWT)承襲了傳統(tǒng)離散小波變換優(yōu)越的時(shí)頻分析和多分辨率分析能力,同時(shí)又克服了傳統(tǒng)離散小波變換的平移敏感性、頻率混疊和缺少方向選擇性等不足[11]。目前,國內(nèi)外學(xué)者利用雙樹復(fù)小波對探地雷達(dá)信號進(jìn)行處理的研究還相對較少。Behrooz Oskooi等人在研究利用雙樹復(fù)小波對探地雷達(dá)信號進(jìn)行去噪處理時(shí),在閾值的選擇上僅考慮了通用閾值[12],而通用閾值存在“過扼殺”小波系數(shù)的情況,會丟失有用的高頻信號[13]。

    本文以川藏鐵路拉林段某隧道為例,在對探測信號實(shí)施雙樹復(fù)小波變換的基礎(chǔ)上,采用通用閾值、無偏似然估計(jì)閾值、啟發(fā)式閾值、極大極小閾值4種閾值選擇方法,結(jié)合軟閾值、硬閾值、半軟閾值、非負(fù)消減4種閾值處理方法,對小波高頻系數(shù)進(jìn)行閾值量化處理后,進(jìn)行雙樹復(fù)小波逆變換重構(gòu)信號,進(jìn)而達(dá)到對探地雷達(dá)信號進(jìn)行去噪的目的,以提高超前地質(zhì)預(yù)報(bào)識別精度。

    1 掌子面測線布置

    根據(jù)前期地質(zhì)勘探資料,結(jié)合隧道施工需求,在隧道掌子面上布置2條測線(圖1),采用SIR-20型探地雷達(dá)和100 MHz天線,對掌子面前方30 m范圍的斷層、節(jié)理裂隙發(fā)育情況等進(jìn)行探測。為提高探地雷達(dá)探測信號質(zhì)量,利用雙樹復(fù)小波變換進(jìn)行4層分解,分別采用4種不同的閾值選擇方法并結(jié)合4種處理方法對小波系數(shù)進(jìn)行閾值處理,用雙樹復(fù)小波逆變換進(jìn)行重構(gòu),并對重構(gòu)后的信號進(jìn)行去噪效果評價(jià)。

    圖1 探地雷達(dá)測線布置示意

    2 基于雙樹復(fù)小波的信號處理方法

    2.1 雙樹復(fù)小波變換

    探地雷達(dá)向被探測的介質(zhì)發(fā)射寬帶電磁波脈沖信號,其中一部分電磁波在遇到不同介電常數(shù)的介質(zhì)界面時(shí)被反射回來,而這些反射回來的電磁波信號就攜帶了大量的信息,工程人員通過分析這些信息,解讀、判識被探測物體內(nèi)部信息;但在探測過程中存在大量環(huán)境干擾,這些回波信息中存在各種噪聲信息,在后期數(shù)據(jù)處理時(shí)需要壓制和去除干擾信號,解讀有效信號。文中采用雙樹復(fù)小波方法進(jìn)行信號數(shù)據(jù)變換,利用雙樹結(jié)構(gòu)、采用滿足希爾伯特變換的濾波器進(jìn)行探地雷達(dá)隧道超前地質(zhì)預(yù)報(bào)信號處理與降噪。

    探地雷達(dá)接收的雷達(dá)信號x(t)可以分解為:

    x(t)=s(t)+n(t),

    (1)

    雙樹復(fù)小波變換對信號x(t)進(jìn)行兩路離散小波變換,即采用2個(gè)獨(dú)立的小波來計(jì)算信號的實(shí)部與虛部。如圖2所示,樹a為實(shí)部,樹b為虛部,由實(shí)部和虛部2個(gè)樹狀的共軛正交濾波器組構(gòu)成一個(gè)復(fù)小波,h0(n)和h1(n)分別對應(yīng)實(shí)部的低通濾波器和高通濾波器,g0(n)和g1(n)分別對應(yīng)虛部的低通濾波器和高通濾波器,且滿足:

    圖2 雙樹復(fù)小波結(jié)構(gòu)示意

    G0(ejw)≈H0(ejw)e-j0.5w,

    (2)

    式中:G0(ejw)和H0(ejw)分別是g0(n)和h0(n)的傅里葉變換[12]。雙樹復(fù)小波的逆變換如圖3所示。

    圖3 雙樹復(fù)小波逆變換

    2.2 小波閾值去噪

    2.2.1 小波去噪信號處理過程

    小波閾值去噪過程分為3個(gè)階段:分解、閾值量化和重構(gòu),在雙樹復(fù)小波變換的基礎(chǔ)上,將含噪信號進(jìn)行分解,保留最底層的低頻系數(shù),僅對各個(gè)分解尺度下的高頻系數(shù)選擇一個(gè)閾值進(jìn)行量化處理,將處理過的小波系數(shù)通過雙樹復(fù)小波逆變換進(jìn)行重構(gòu),進(jìn)而獲得去除噪聲后的信號(圖4)。這個(gè)過程中最關(guān)鍵的是如何選擇閾值,并利用該閾值對小波系數(shù)進(jìn)行量化處理,它關(guān)系到信號消噪的質(zhì)量。

    圖4 雙樹復(fù)小波去噪流程示意

    2.2.2 閾值選擇及閾值處理

    閾值選擇直接關(guān)系著去噪結(jié)果,過大的閾值導(dǎo)致過多小波系數(shù)萎縮,去噪后信號過于平滑;而過小的閾值能最大限度保留信號的形狀,卻無法抑制噪聲。各閾值選擇方法各有優(yōu)缺點(diǎn),本文采用通用閾值(universal threshold)、無偏似然估計(jì)閾值(rigrsure)、啟發(fā)式閾值(heursure)、極大極小閾值(minimaxi)[3]4種方法進(jìn)行閾值選擇。

    1)通用閾值

    Donoho在高斯噪聲模型下,應(yīng)用多維獨(dú)立正態(tài)變量決策理論,得到了通用的閾值:

    (3)

    式中:σ為噪聲的標(biāo)準(zhǔn)差,N為信號長度。

    2)無偏似然估計(jì)閾值

    無偏似然估計(jì)閾值是一種基于 Stein 的無偏似然估計(jì)原理的自適應(yīng)閾值選擇,是一種軟件閾值估計(jì)器。

    3)啟發(fā)式閾值

    啟發(fā)式閾值是前兩種閾值的綜合,所選擇的是最優(yōu)預(yù)測變量閾值,視信號情況選定通用閾值或無偏似然估計(jì)閾值。

    4)極大極小閾值

    極值閾值也是一種固定閾值選擇形式,它所產(chǎn)生的是一個(gè)最小均方差的極值,不是無誤差。主要應(yīng)用極大極小原理選擇閾值[3]。

    (4)

    式中:wi為尺度i上各點(diǎn)的小波系數(shù)。當(dāng)噪聲為非白噪聲時(shí),可以在每個(gè)分解尺度上都估計(jì)噪聲的方差,從而計(jì)算不同分解尺度上的閾值。由于探地雷達(dá)噪聲受環(huán)境干擾較大,存在非白噪聲干擾,因此對不同分解層次的小波系數(shù)采用不同的閾值,以達(dá)到優(yōu)化去噪效果。

    在把信號進(jìn)行分解后,要將得到的高頻小波系數(shù)w(i)與閾值T比較。為了獲得準(zhǔn)確的數(shù)據(jù),在采用軟、硬閾值處理方法的同時(shí),又引入了半軟閾值處理(firm shrinkage)[14]、非負(fù)消減處理(nonnegative garrote shrinkage)[15]2種改進(jìn)的閾值處理方法。

    1)硬閾值處理

    (5)

    硬閾值處理將小波系數(shù)中絕對值小于閾值T的小波系數(shù)置零,因此會產(chǎn)生不連續(xù)的情況。

    2)軟閾值處理

    (6)

    軟閾值處理后的小波系數(shù)雖然整體連續(xù)性好,但是對所有的小波系數(shù)都進(jìn)行了抑制,造成高頻信息損失,因此給重構(gòu)帶了不小的誤差[13]。

    5) 將 gmchcfg.h 文件放置在C:Tornado221 argethugldrivergraphicsintel 路徑下(本文假設(shè). VxWorks 開發(fā)環(huán)境安裝在 C:Tornado221 路徑下);

    3)半軟閾值處理

    (7)

    Gao和Bruce于1997年提出半軟閾值處理[14],其中T2=2T。

    4)非負(fù)消減處理

    (8)

    Breiman于1995年提出了非負(fù)消減處理[15],這種方法在軟、硬閾值處理之間做了折中。

    2.3 去噪效果評價(jià)

    為評價(jià)雙樹復(fù)小波變換去噪重構(gòu)后的性能,采用以下3個(gè)參數(shù)進(jìn)行評價(jià)。

    1)歸一化均方根誤差NRMSE:

    (9)

    2)信噪比SNR:

    (10)

    3)能量比P:

    (11)

    式中:x(i)是探地雷達(dá)采集到的原始信號,x′(i)為雙樹復(fù)小波去噪重構(gòu)后的信號,ux為原始信號x(i)的平均值。歸一化均方根誤差描述了去噪信號x′(i)與原始信號x(i)之間的相似程度[5],信噪比指的是有效波與干擾波相對強(qiáng)弱的比較,能量比為降噪后的信號在原信號所占的能量比例;歸一化均方根誤差越小,信噪比越大,能量比越大,說明信號失真越小,信號去噪效果越好[4]。

    本文對探地雷達(dá)數(shù)據(jù)進(jìn)行雙樹復(fù)小波變換去噪,采用4種閾值選擇和4種閾值處理方法兩兩組合,通過對去噪效果的評價(jià),并結(jié)合去噪前后波形圖及堆積圖的對比,找到最佳的閾值選擇和閾值處理方法的組合。

    3 實(shí)測數(shù)據(jù)處理效果

    3.1 1#掌子面應(yīng)用情況

    在隧道1#掌子面上水平布置2條測線,其中測線2數(shù)據(jù)共計(jì)415道,單道數(shù)據(jù)采樣點(diǎn)數(shù)1 024,記錄長度600 ns。對測線2數(shù)據(jù)進(jìn)行雙樹復(fù)小波變換去噪并重構(gòu),以第250道數(shù)據(jù)為例進(jìn)行評價(jià)分析。

    從表1可以看出: 去噪后的歸一化均方差NRMSE值都比較小,在0.1左右,但采用無偏似然估計(jì)或啟發(fā)式閾值選擇方法,結(jié)合非負(fù)消減閾值處理方法時(shí),得到的歸一化均方差值最小,僅為0.04;去噪后的SNR信號信噪比較高,均在75~86.7 dB之間,但采用無偏似然估計(jì)或啟發(fā)式閾值選擇方法,結(jié)合非負(fù)消減閾值處理方法時(shí),得到的信噪比值最大;去噪后的信號P值信號能量比相差較大,采用無偏似然估計(jì)閾值選擇方法,結(jié)合非負(fù)消減閾值處理方法時(shí),得到的信號能量比最大,達(dá)到99.06%。圖5為第250道采用無偏似然估計(jì)閾值選擇方法和非負(fù)消減閾值處理方法重構(gòu)的信號波形,圖中顯示多次波及大部分噪聲被明顯去除。

    圖5 第250道單道原始信號與雙樹復(fù)小波變換去噪信號對比

    表1 第250道數(shù)據(jù)去噪效果評價(jià)

    通過歸一化均方根誤差、信噪比、能量比的結(jié)果分析,采用無偏似然估計(jì)閾值選擇方法和非負(fù)消減閾值處理方法對1#掌子面信號進(jìn)行雙樹復(fù)小波變換后的去噪效果最佳。

    從整個(gè)掌子面探測剖面信號的處理情況看(圖6),采用無偏似然估計(jì)閾值選擇方法,結(jié)合非負(fù)消減閾值處理方法進(jìn)行雙樹復(fù)小波變換去噪處理之后,噪聲干擾被有效消除,消減了多次反射波,異常體的邊界及位置清晰可見,比傳統(tǒng)濾波方法去噪效果更好。圖6中可以清晰看到250~350道之間的強(qiáng)反射區(qū)波形存在明顯畸變,同相軸連續(xù)性較差,結(jié)合地質(zhì)勘察資料與掌子面情況,可以推斷該段圍巖較破碎,節(jié)理裂隙較發(fā)育,穩(wěn)定性較差,施工時(shí)應(yīng)加強(qiáng)支護(hù)。

    圖6 1#掌子面?zhèn)鹘y(tǒng)濾波與雙樹復(fù)小波去噪結(jié)果對比

    3.2 2#掌子面應(yīng)用情況

    在隧道2#掌子面上水平布置2條測線,其中測線1數(shù)據(jù)共計(jì)800道,單道數(shù)據(jù)采樣點(diǎn)數(shù)1 024,記錄長度600 ns。對測線1數(shù)據(jù)進(jìn)行雙樹復(fù)小波變換去噪并重構(gòu),以第700道數(shù)據(jù)重構(gòu)信號為例進(jìn)行去噪效果評價(jià)。

    從表2可以看出:采用無偏似然估計(jì)或啟發(fā)式閾值選擇方法,結(jié)合非負(fù)消減閾值處理方法得到的歸一化均方差NRMSE值最小,還不到0.03,達(dá)到了較好的去噪效果,而采用其他方法時(shí),歸一化均方差NRMSE值相對較大,去噪效果不理想;采用無偏似然估計(jì)閾值或啟發(fā)式閾值選擇方法,結(jié)合非負(fù)消減閾值處理方法時(shí),去噪后的信噪比較大,可達(dá)到60 dB以上,信號能量遠(yuǎn)大于噪聲能量,去噪效果較好,而其他方法去噪后的信噪比只有27~30 dB,去噪效果一般;采用無偏似然估計(jì)閾值選擇方法,結(jié)合非負(fù)消減閾值處理方法時(shí),能量比值P最大,達(dá)到98.88%,去噪效果較好,而采用其他方法得到的能量比值P相對較小,去噪效果不理想。采用無偏似然估計(jì)閾值選擇方法和非負(fù)消減閾值處理方法重構(gòu)信號的波形圖如圖7所示,圖中顯示多次波及大部分噪聲被明顯去除。

    圖7 第700道原始信號與雙樹復(fù)小波變換去噪信號對比

    表2 第700道數(shù)據(jù)去噪效果評價(jià)

    通過歸一化均方根誤差、信噪比、能量比的結(jié)果分析,采用無偏似然估計(jì)閾值選擇方法和非負(fù)消減閾值處理方法對2#掌子面信號進(jìn)行雙樹復(fù)小波變換后的去噪效果最佳。從整個(gè)掌子面探測剖面信號的處理情況看(圖8),采用無偏似然估計(jì)閾值選擇方法并結(jié)合非負(fù)消減閾值處理方法進(jìn)行雙樹復(fù)小波變換去噪處理之后,噪聲干擾被有效消除,消減了多次反射波,異常體的邊界及位置清晰可見,比傳統(tǒng)濾波方法去噪效果更好。從圖8中可以清晰看到600~750道之間的強(qiáng)反射區(qū)波形存在明顯畸變,同相軸連續(xù)性較差;結(jié)合地質(zhì)勘察資料與掌子面情況,可以推斷該段圍巖較破碎,節(jié)理裂隙較發(fā)育,穩(wěn)定性較差,施工時(shí)應(yīng)加強(qiáng)支護(hù)。

    圖8 2#掌子面?zhèn)鹘y(tǒng)濾波與雙樹復(fù)小波去噪結(jié)果對比

    4 結(jié)論

    受川藏鐵路極其復(fù)雜的地質(zhì)條件和施工環(huán)境的影響,探地雷達(dá)地質(zhì)超前預(yù)報(bào)技術(shù)采集的信號噪聲干擾較大,給信號處理與數(shù)據(jù)分析帶來較大困難。在對探測信號實(shí)施雙樹復(fù)小波變換的基礎(chǔ)上,采用通用閾值、無偏似然估計(jì)閾值、啟發(fā)式閾值、極大極小閾值4種閾值選擇方法,結(jié)合軟閾值、硬閾值、半軟閾值、非負(fù)消減4種閾值處理方法,對小波高頻系數(shù)進(jìn)行閾值量化處理后進(jìn)行雙樹復(fù)小波逆變換重構(gòu)信號,達(dá)到了較好的去噪效果。對比分析歸一化均方差、信噪比、能量比3個(gè)評價(jià)因子后發(fā)現(xiàn),采用無偏似然估計(jì)閾值選擇方法,結(jié)合非負(fù)消減閾值處理方法,可以獲得最佳的去噪效果,信號得到很大的改善,噪聲殘留較少,多次反射波被減弱,提高了超前地質(zhì)預(yù)報(bào)中異常體的識別精度,能夠?yàn)榇ú罔F路隧道施工提供更好的安全保障。

    猜你喜歡
    雙樹探地掌子面
    基于超表面的探地雷達(dá)增強(qiáng)探測研究
    全極化探地雷達(dá)系統(tǒng)
    斜井掌子面形狀與傾角對隧道開挖面穩(wěn)定性影響
    構(gòu)造破碎帶隧道掌子面穩(wěn)定性及加固效果分析
    鐵道建筑(2019年7期)2019-08-08 07:36:58
    一個(gè)村莊的紅色記憶
    飽水砂性地層隧道掌子面穩(wěn)定性分析
    上軟下硬地層隧道掌子面穩(wěn)定性及塌方形態(tài)
    基于雙樹復(fù)小波的色譜重疊峰分解方法研究
    基于探地雷達(dá)法的地下管線探測頻譜分析
    婆羅雙樹樣基因2干擾對宮頸癌HeLa細(xì)胞增殖和凋亡的影響
    校园人妻丝袜中文字幕| 免费无遮挡裸体视频| 天堂影院成人在线观看| 能在线免费观看的黄片| 性插视频无遮挡在线免费观看| 久久久久久国产a免费观看| 亚洲综合色惰| 精品国产三级普通话版| 国产男靠女视频免费网站| 国产精品综合久久久久久久免费| 搡老熟女国产l中国老女人| 啦啦啦观看免费观看视频高清| 日韩欧美在线二视频| 一级av片app| 日韩欧美精品免费久久| 深爱激情五月婷婷| 国产精品自产拍在线观看55亚洲| 色吧在线观看| 99久久精品热视频| 97人妻精品一区二区三区麻豆| 国产乱人视频| 婷婷精品国产亚洲av| 久久久久国产精品人妻aⅴ院| 国产亚洲精品综合一区在线观看| 亚洲在线自拍视频| 午夜久久久久精精品| 女的被弄到高潮叫床怎么办 | 给我免费播放毛片高清在线观看| 露出奶头的视频| 国产一区二区三区在线臀色熟女| 久久久久久国产a免费观看| 欧美黑人欧美精品刺激| 草草在线视频免费看| 少妇丰满av| 免费看a级黄色片| 国内少妇人妻偷人精品xxx网站| 久久精品91蜜桃| 亚洲电影在线观看av| 国产精品久久久久久久电影| 变态另类成人亚洲欧美熟女| 看十八女毛片水多多多| 国产精品亚洲美女久久久| 18+在线观看网站| a级一级毛片免费在线观看| 欧美潮喷喷水| 天堂网av新在线| 免费人成在线观看视频色| 中文字幕久久专区| 国内精品久久久久精免费| 看片在线看免费视频| 中国美白少妇内射xxxbb| 国产探花在线观看一区二区| 亚洲av成人精品一区久久| 真实男女啪啪啪动态图| 级片在线观看| 十八禁国产超污无遮挡网站| 亚洲图色成人| 又爽又黄a免费视频| 中文字幕人妻熟人妻熟丝袜美| 蜜桃久久精品国产亚洲av| 免费看日本二区| 97超级碰碰碰精品色视频在线观看| 国产人妻一区二区三区在| 真人一进一出gif抽搐免费| 午夜免费激情av| 日韩高清综合在线| 精品久久久久久久末码| a在线观看视频网站| 亚洲欧美清纯卡通| 中国美女看黄片| 91久久精品国产一区二区成人| 日韩欧美免费精品| 我要搜黄色片| 岛国在线免费视频观看| 亚洲熟妇熟女久久| 亚洲国产欧洲综合997久久,| 国产伦精品一区二区三区视频9| 亚洲欧美日韩无卡精品| 亚洲成人精品中文字幕电影| 日本免费一区二区三区高清不卡| 免费观看在线日韩| 欧美日韩瑟瑟在线播放| 中国美白少妇内射xxxbb| 国产精品不卡视频一区二区| 成人美女网站在线观看视频| 国产又黄又爽又无遮挡在线| 亚洲天堂国产精品一区在线| 蜜桃久久精品国产亚洲av| 精品日产1卡2卡| 国产视频内射| 欧美丝袜亚洲另类 | 欧美+日韩+精品| 午夜免费男女啪啪视频观看 | 俺也久久电影网| 一区二区三区四区激情视频 | 非洲黑人性xxxx精品又粗又长| 99riav亚洲国产免费| 男女视频在线观看网站免费| 天美传媒精品一区二区| av在线天堂中文字幕| 男女之事视频高清在线观看| 国产精品自产拍在线观看55亚洲| 成人特级av手机在线观看| or卡值多少钱| 91久久精品国产一区二区三区| 亚洲av美国av| 亚洲综合色惰| 天堂动漫精品| 哪里可以看免费的av片| 国产精品,欧美在线| 久久亚洲精品不卡| 少妇人妻一区二区三区视频| 2021天堂中文幕一二区在线观| 国产精品无大码| 欧美日本亚洲视频在线播放| 12—13女人毛片做爰片一| 免费观看在线日韩| 一进一出抽搐gif免费好疼| 久久久久国产精品人妻aⅴ院| 看十八女毛片水多多多| 久久久久国产精品人妻aⅴ院| 婷婷精品国产亚洲av| 九九在线视频观看精品| 日本黄色视频三级网站网址| 亚洲第一区二区三区不卡| 永久网站在线| 老熟妇仑乱视频hdxx| 欧美又色又爽又黄视频| 亚洲真实伦在线观看| av国产免费在线观看| 国产一区二区在线观看日韩| 日韩欧美一区二区三区在线观看| 国产精品,欧美在线| 亚洲av成人av| 国产精品自产拍在线观看55亚洲| av视频在线观看入口| 日本一二三区视频观看| 国产精华一区二区三区| 日本欧美国产在线视频| 男插女下体视频免费在线播放| 日韩 亚洲 欧美在线| 非洲黑人性xxxx精品又粗又长| 免费看av在线观看网站| 特大巨黑吊av在线直播| 午夜福利成人在线免费观看| 精品人妻一区二区三区麻豆 | 制服丝袜大香蕉在线| 国产精品自产拍在线观看55亚洲| 男女视频在线观看网站免费| 桃红色精品国产亚洲av| 夜夜夜夜夜久久久久| 九九爱精品视频在线观看| 大又大粗又爽又黄少妇毛片口| 午夜视频国产福利| 老熟妇乱子伦视频在线观看| 成人特级黄色片久久久久久久| 亚洲专区中文字幕在线| 男女下面进入的视频免费午夜| 色综合亚洲欧美另类图片| 精品午夜福利在线看| 免费大片18禁| 亚洲欧美日韩卡通动漫| 最近在线观看免费完整版| 老司机午夜福利在线观看视频| av在线老鸭窝| 国产伦在线观看视频一区| 一个人看的www免费观看视频| 亚洲无线在线观看| 少妇被粗大猛烈的视频| 日韩大尺度精品在线看网址| 亚洲人成网站高清观看| 亚洲avbb在线观看| 又黄又爽又刺激的免费视频.| 免费av观看视频| 亚洲欧美日韩东京热| 日本黄色视频三级网站网址| 国产免费一级a男人的天堂| 婷婷色综合大香蕉| 深夜精品福利| 美女xxoo啪啪120秒动态图| 夜夜爽天天搞| 人妻少妇偷人精品九色| 欧美国产日韩亚洲一区| 成人综合一区亚洲| 国产成人一区二区在线| 日本黄大片高清| 91av网一区二区| 深夜精品福利| 观看美女的网站| 嫁个100分男人电影在线观看| 日韩 亚洲 欧美在线| 如何舔出高潮| 人妻久久中文字幕网| 午夜精品久久久久久毛片777| 亚洲精华国产精华液的使用体验 | 久久精品国产亚洲av天美| 国产成人影院久久av| 午夜福利视频1000在线观看| 成人鲁丝片一二三区免费| 亚洲中文日韩欧美视频| 日本成人三级电影网站| 男女下面进入的视频免费午夜| 91久久精品国产一区二区三区| 毛片女人毛片| 日本三级黄在线观看| 国产伦在线观看视频一区| 身体一侧抽搐| 搡女人真爽免费视频火全软件 | 中文亚洲av片在线观看爽| 亚洲av中文av极速乱 | 九九爱精品视频在线观看| netflix在线观看网站| 不卡一级毛片| 日韩精品有码人妻一区| 国内毛片毛片毛片毛片毛片| 1024手机看黄色片| 直男gayav资源| 亚洲,欧美,日韩| 22中文网久久字幕| 国产男人的电影天堂91| 午夜福利在线观看吧| 一级黄片播放器| 国产高清三级在线| h日本视频在线播放| 精品久久久久久,| 欧美日韩中文字幕国产精品一区二区三区| 国产精品久久久久久av不卡| 如何舔出高潮| 欧美日本视频| 日韩人妻高清精品专区| 亚洲精品色激情综合| 久久天躁狠狠躁夜夜2o2o| 成年女人毛片免费观看观看9| 99精品在免费线老司机午夜| 欧美三级亚洲精品| 国产伦精品一区二区三区视频9| 午夜亚洲福利在线播放| 亚洲人成伊人成综合网2020| 99热精品在线国产| 白带黄色成豆腐渣| 午夜福利18| 综合色av麻豆| 一边摸一边抽搐一进一小说| 别揉我奶头~嗯~啊~动态视频| 国产aⅴ精品一区二区三区波| 午夜免费男女啪啪视频观看 | 国产伦一二天堂av在线观看| 尾随美女入室| 国产三级中文精品| 日韩欧美免费精品| 久久久久性生活片| 免费观看精品视频网站| 三级男女做爰猛烈吃奶摸视频| 久久热精品热| 久久久久久大精品| 日本 欧美在线| 九九爱精品视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一区高清亚洲精品| 又爽又黄a免费视频| 尾随美女入室| 日韩欧美 国产精品| 亚洲av美国av| 久久国内精品自在自线图片| 免费黄网站久久成人精品| 特级一级黄色大片| 欧美又色又爽又黄视频| 国产成年人精品一区二区| 国模一区二区三区四区视频| 日本撒尿小便嘘嘘汇集6| 可以在线观看毛片的网站| 人妻夜夜爽99麻豆av| 十八禁网站免费在线| 毛片一级片免费看久久久久 | 成人欧美大片| 日韩中文字幕欧美一区二区| 久久久成人免费电影| 国产精品国产三级国产av玫瑰| 欧美在线一区亚洲| 女的被弄到高潮叫床怎么办 | 我要搜黄色片| 亚洲人与动物交配视频| 午夜视频国产福利| 成人av一区二区三区在线看| 亚洲国产欧美人成| 午夜爱爱视频在线播放| 啦啦啦韩国在线观看视频| 久久香蕉精品热| 丰满的人妻完整版| 岛国在线免费视频观看| 啦啦啦啦在线视频资源| 搡老熟女国产l中国老女人| 97碰自拍视频| 中文字幕熟女人妻在线| 非洲黑人性xxxx精品又粗又长| 国产伦精品一区二区三区四那| 老熟妇仑乱视频hdxx| 国产一区二区亚洲精品在线观看| 成年人黄色毛片网站| 国产 一区 欧美 日韩| 国产精品1区2区在线观看.| 男女下面进入的视频免费午夜| 啪啪无遮挡十八禁网站| 欧美日韩亚洲国产一区二区在线观看| 婷婷亚洲欧美| 亚洲av.av天堂| 禁无遮挡网站| 国产在线精品亚洲第一网站| 麻豆久久精品国产亚洲av| 色综合站精品国产| 国产精品美女特级片免费视频播放器| 日韩欧美在线乱码| 波多野结衣高清无吗| 18禁在线播放成人免费| 精品久久久久久成人av| 少妇高潮的动态图| 日本 av在线| 国产精品久久视频播放| 久久精品国产亚洲av天美| 欧美日韩综合久久久久久 | 一区二区三区四区激情视频 | 日韩欧美一区二区三区在线观看| 亚洲黑人精品在线| 亚洲第一电影网av| 白带黄色成豆腐渣| 亚洲av免费在线观看| 在线免费十八禁| netflix在线观看网站| 精品人妻1区二区| 极品教师在线免费播放| 麻豆av噜噜一区二区三区| 欧美绝顶高潮抽搐喷水| 精品一区二区免费观看| 久久国内精品自在自线图片| 国产在线男女| 婷婷精品国产亚洲av在线| 听说在线观看完整版免费高清| 日本撒尿小便嘘嘘汇集6| 色综合色国产| 国产精品,欧美在线| 精品人妻熟女av久视频| 少妇的逼水好多| 免费看a级黄色片| av视频在线观看入口| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久久电影| 十八禁网站免费在线| 亚洲成人久久爱视频| 国产精品不卡视频一区二区| 午夜福利18| 一个人观看的视频www高清免费观看| 色尼玛亚洲综合影院| 中文亚洲av片在线观看爽| 精品久久久久久久久av| 97超级碰碰碰精品色视频在线观看| 18+在线观看网站| 此物有八面人人有两片| 噜噜噜噜噜久久久久久91| 欧美中文日本在线观看视频| 亚洲最大成人中文| 午夜福利18| 亚洲av熟女| 国产 一区 欧美 日韩| 国产极品精品免费视频能看的| 国产91精品成人一区二区三区| 亚洲最大成人手机在线| 日韩大尺度精品在线看网址| 国产一级毛片七仙女欲春2| 极品教师在线免费播放| 欧美极品一区二区三区四区| 亚洲专区国产一区二区| 久久精品影院6| 亚洲欧美精品综合久久99| .国产精品久久| 亚洲精品国产成人久久av| 一卡2卡三卡四卡精品乱码亚洲| 网址你懂的国产日韩在线| 久久亚洲精品不卡| 国内精品宾馆在线| 亚洲成av人片在线播放无| 久久精品91蜜桃| 国产精品爽爽va在线观看网站| 午夜福利成人在线免费观看| 亚洲人成网站在线播放欧美日韩| 亚洲av免费在线观看| 1024手机看黄色片| 能在线免费观看的黄片| 国产探花在线观看一区二区| 日本五十路高清| 午夜老司机福利剧场| av专区在线播放| 亚洲乱码一区二区免费版| 久久亚洲真实| 国产综合懂色| 国产精品国产高清国产av| 久久久成人免费电影| 久久6这里有精品| 国产v大片淫在线免费观看| 夜夜看夜夜爽夜夜摸| 亚洲狠狠婷婷综合久久图片| 日本三级黄在线观看| 国产女主播在线喷水免费视频网站 | 看免费成人av毛片| 他把我摸到了高潮在线观看| 简卡轻食公司| 网址你懂的国产日韩在线| 久久久久久久午夜电影| 如何舔出高潮| 1000部很黄的大片| 欧美在线一区亚洲| 色综合婷婷激情| 丰满人妻一区二区三区视频av| 国产日本99.免费观看| 国产高清不卡午夜福利| 国产一区二区激情短视频| 亚洲国产精品久久男人天堂| av天堂在线播放| 一本久久中文字幕| 国产精品亚洲美女久久久| 欧美精品国产亚洲| 91麻豆av在线| 国产一区二区三区视频了| 亚洲av第一区精品v没综合| 国产成人a区在线观看| 国产乱人视频| 中亚洲国语对白在线视频| 麻豆一二三区av精品| 搡女人真爽免费视频火全软件 | 亚洲国产欧美人成| 亚洲精品乱码久久久v下载方式| 久久热精品热| 中国美女看黄片| 色5月婷婷丁香| 亚洲av免费在线观看| 最新在线观看一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 乱人视频在线观看| 又黄又爽又刺激的免费视频.| 午夜影院日韩av| 麻豆国产av国片精品| 伦理电影大哥的女人| 中文字幕熟女人妻在线| 精品人妻一区二区三区麻豆 | 欧美激情国产日韩精品一区| 免费不卡的大黄色大毛片视频在线观看 | 99久久九九国产精品国产免费| 亚洲国产欧美人成| 日韩欧美国产在线观看| 搡老岳熟女国产| 51国产日韩欧美| 日本黄大片高清| 久久天躁狠狠躁夜夜2o2o| 国产伦在线观看视频一区| 日本精品一区二区三区蜜桃| 人妻丰满熟妇av一区二区三区| 熟妇人妻久久中文字幕3abv| 亚洲国产欧洲综合997久久,| 亚洲av不卡在线观看| 亚洲中文字幕一区二区三区有码在线看| 亚洲av第一区精品v没综合| 欧美潮喷喷水| 色综合亚洲欧美另类图片| 久久久久性生活片| 国产午夜福利久久久久久| 国产黄色小视频在线观看| 日韩大尺度精品在线看网址| 一a级毛片在线观看| 亚洲人成伊人成综合网2020| 麻豆成人av在线观看| 国产高清有码在线观看视频| 99精品在免费线老司机午夜| 国产精品av视频在线免费观看| 97超级碰碰碰精品色视频在线观看| 国产一区二区三区视频了| 91久久精品国产一区二区成人| 欧美高清成人免费视频www| 美女黄网站色视频| 97超级碰碰碰精品色视频在线观看| 动漫黄色视频在线观看| av视频在线观看入口| 国产精品电影一区二区三区| 国内揄拍国产精品人妻在线| 亚洲人成伊人成综合网2020| 国产一级毛片七仙女欲春2| 天天躁日日操中文字幕| 国产精品亚洲美女久久久| 亚洲内射少妇av| 在线播放国产精品三级| 国产 一区精品| 国产成人aa在线观看| 夜夜看夜夜爽夜夜摸| 色播亚洲综合网| 不卡视频在线观看欧美| a级毛片a级免费在线| 久久久精品大字幕| 精华霜和精华液先用哪个| 三级毛片av免费| 免费在线观看影片大全网站| 午夜老司机福利剧场| 国产探花极品一区二区| 99在线人妻在线中文字幕| 中文字幕熟女人妻在线| 免费看av在线观看网站| 自拍偷自拍亚洲精品老妇| 久久精品国产亚洲网站| 日韩中文字幕欧美一区二区| 一区二区三区四区激情视频 | 国产高清有码在线观看视频| 国产精品不卡视频一区二区| 99热网站在线观看| 老司机福利观看| 国产真实伦视频高清在线观看 | 国内精品美女久久久久久| 午夜福利在线观看吧| 久久久久国产精品人妻aⅴ院| 国产精品女同一区二区软件 | 亚洲国产欧洲综合997久久,| 国内久久婷婷六月综合欲色啪| 亚洲精品日韩av片在线观看| 99热6这里只有精品| 一区二区三区高清视频在线| 亚洲欧美清纯卡通| 三级男女做爰猛烈吃奶摸视频| 国产高清有码在线观看视频| 亚洲性久久影院| or卡值多少钱| 18+在线观看网站| 亚洲avbb在线观看| 亚洲美女黄片视频| 亚洲欧美激情综合另类| 午夜福利高清视频| 日本欧美国产在线视频| 国产精品亚洲一级av第二区| 啪啪无遮挡十八禁网站| 亚洲va日本ⅴa欧美va伊人久久| 观看美女的网站| 91狼人影院| 欧美性感艳星| 欧美日韩亚洲国产一区二区在线观看| 欧美一区二区精品小视频在线| 好男人在线观看高清免费视频| 亚洲专区国产一区二区| av视频在线观看入口| 欧美绝顶高潮抽搐喷水| 国产一区二区激情短视频| 国产精品野战在线观看| 少妇的逼好多水| 99热精品在线国产| 哪里可以看免费的av片| 国产色婷婷99| 国产亚洲91精品色在线| 联通29元200g的流量卡| 五月玫瑰六月丁香| 黄色欧美视频在线观看| 国产精品国产三级国产av玫瑰| 九色国产91popny在线| 香蕉av资源在线| 黄色丝袜av网址大全| av在线观看视频网站免费| 国产中年淑女户外野战色| 91午夜精品亚洲一区二区三区 | 国产高清视频在线播放一区| 午夜激情福利司机影院| 国产三级在线视频| 亚洲久久久久久中文字幕| www.色视频.com| 亚洲天堂国产精品一区在线| 我的女老师完整版在线观看| 亚洲久久久久久中文字幕| 草草在线视频免费看| 久久精品国产鲁丝片午夜精品 | 国产精品久久电影中文字幕| 搡老岳熟女国产| 亚洲男人的天堂狠狠| 国产精品av视频在线免费观看| 久久精品国产亚洲av天美| 男人舔女人下体高潮全视频| 亚洲不卡免费看| 色综合亚洲欧美另类图片| 国产精品久久电影中文字幕| 黄色女人牲交| 国产精品永久免费网站| 久久久久久伊人网av| 又爽又黄a免费视频| 亚洲一级一片aⅴ在线观看| 久久精品国产鲁丝片午夜精品 | 国产精品久久电影中文字幕| 一区福利在线观看| 99久久成人亚洲精品观看| 亚洲一级一片aⅴ在线观看| 国产综合懂色| 国产男人的电影天堂91| 波多野结衣巨乳人妻| 国产精品自产拍在线观看55亚洲| 在线观看舔阴道视频| 岛国在线免费视频观看| 伦理电影大哥的女人| 在线观看午夜福利视频| 国产一级毛片七仙女欲春2| av女优亚洲男人天堂| 人人妻,人人澡人人爽秒播| 国产亚洲欧美98| 日本色播在线视频| 日本在线视频免费播放| 久久久久久久久大av| 亚洲三级黄色毛片| 成人特级av手机在线观看| 国产伦一二天堂av在线观看| 女人十人毛片免费观看3o分钟| 国产毛片a区久久久久| 久久久久九九精品影院| 不卡一级毛片| 韩国av在线不卡| 婷婷色综合大香蕉| 国产在视频线在精品|