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

    基于小波和CEEMDAN的地震信號(hào)噪聲壓制

    2021-05-28 03:02:48方根顯
    物探化探計(jì)算技術(shù) 2021年3期
    關(guān)鍵詞:小波分量重構(gòu)

    殷 明, 方根顯

    (東華理工大學(xué) 地球物理與測(cè)控技術(shù)學(xué)院,南昌 330013)

    0 前言

    淺層地震勘探由于其精度高,在城市地下空間勘探中得到了廣泛地應(yīng)用,長(zhǎng)期以來(lái)為工程建設(shè)、城市規(guī)劃等提供了重要的支撐信息[1],因此在城市地質(zhì)調(diào)查中被作為常用的探測(cè)技術(shù)。由于城市特殊的地理環(huán)境的影響,地震信號(hào)不僅受各種各樣的人為因素干擾,而且還有天然的環(huán)境因素干擾,因此淺層地震勘探的探測(cè)效果在復(fù)雜的環(huán)境中受干擾影響大。對(duì)于非平穩(wěn)非線性的地震信號(hào), huang等[2]提出了EMD(經(jīng)驗(yàn)?zāi)B(tài)分解)克服了傅里葉變換的局限性,但使用EMD容易出現(xiàn)模態(tài)混淆[3]問(wèn)題,在使用后期對(duì)其進(jìn)行了改進(jìn)算法-集合經(jīng)驗(yàn)?zāi)B(tài)分解 (EEMD, ensemble empirical mode decomposition)[4],但由此引入的白噪聲并不能完全剔除。所以為了解決這一問(wèn)題,CEEMDAN(自適應(yīng)噪聲完備集合經(jīng)驗(yàn)?zāi)B(tài)分解)便被提出,這種方法的提出有效地克服了EEMD的不足。

    針對(duì)地震信號(hào)中的噪聲干擾不易壓制的問(wèn)題,使用小波-CEEMDAN相結(jié)合的方式,對(duì)其信號(hào)進(jìn)行濾波。首先對(duì)信號(hào)進(jìn)行小波分解,對(duì)其高頻系數(shù)中存在的噪聲信號(hào),使用CEEMDAN進(jìn)行分解,隨后對(duì)分解的IMF分量求其自相關(guān)函數(shù)并判斷噪聲主要存在的分量并進(jìn)行剔除,在對(duì)其與小波低頻系數(shù)進(jìn)行重構(gòu)并求出信噪比和均方差,判斷濾波后的效果,最后通過(guò)實(shí)際地震信號(hào)去噪對(duì)比效果,驗(yàn)證本文小波-CEEMDAN噪聲壓制方法的有效性。

    1 淺層地震信號(hào)干擾及噪聲壓制方法

    1.1 噪聲干擾

    在城市地震勘探中噪聲干擾是不確定的,有規(guī)則干擾波與不規(guī)則干擾波。規(guī)則干擾波可以用常規(guī)的方法去壓制噪聲,如帶通濾波、二維濾波、反褶積等[5],但這些方法對(duì)于不規(guī)則干擾卻不再適用,所以一般采用常規(guī)手段和小波變換相結(jié)合的方式進(jìn)行噪聲壓制。

    1)50 Hz工業(yè)電干擾。由于地震測(cè)線的上方存在高壓輸電線路時(shí),會(huì)導(dǎo)致檢波器感應(yīng)到50 Hz的電壓,出現(xiàn)電干擾[6-8],這種50 Hz的電干擾,可以用單頻干擾消除[8]。

    2)地下建筑產(chǎn)生干擾次波。地下建筑物和松散層之間存在明顯的物性差異,當(dāng)?shù)卣鸩▊鞑サ竭@些界面時(shí),會(huì)導(dǎo)致干擾波的形成(如折射,反射等),所以通過(guò)改變激發(fā)和觀察位置能明顯減少其影響。

    3)高、低頻背景干擾。當(dāng)?shù)卣鸩ㄔ谒缮⒔橘|(zhì)中激發(fā)時(shí),低頻干擾(10 Hz~30 Hz)會(huì)隨震源振動(dòng)而形成。震源在硬介質(zhì)中激發(fā)時(shí),會(huì)導(dǎo)致80 Hz到200 Hz的高頻干擾。減弱或消除高、低頻干擾波對(duì)有效波的干擾,可以采用帶通濾波截至高低頻帶的方法。

    4)虛反射。它是由于地震波向上傳播,當(dāng)遇到地面后或者低速帶底部時(shí)又轉(zhuǎn)為向下反射傳播,最后從下方反射界面再次反射到地面的現(xiàn)象[9]。由于該干擾波頻帶較寬, 可以利用反褶積來(lái)減弱或消除對(duì)有效波的干擾。

    5)聲波。它常??梢栽诘卣鹩涗浿姓业?,它的頻率較高(大于100 Hz),特點(diǎn)是波速穩(wěn)定(330 m/s ~340 m/s)。在淺層地震勘探中,特別是地面震源的聲波干擾較為嚴(yán)重[10],對(duì)此可以利用二維速度濾波方法消除。

    6)面波。主要是地表傳播的瑞雷面波,它在淺層地震勘探中是一種主要的干擾波。面波的主要特點(diǎn)是視速度小,頻率低(10 Hz~30 Hz)、衰減慢、能量強(qiáng),對(duì)此可通過(guò)提高低截止頻率來(lái)衰減干擾。

    7)微動(dòng)。是由人為因素所導(dǎo)致的一系列不規(guī)則的振動(dòng)或者大自然當(dāng)中的一些噪聲。這類(lèi)型的噪聲頻帶很寬(1 Hz~200 Hz),并且難以區(qū)分,對(duì)此可以通過(guò)小波變換與CEEMDAN進(jìn)行噪聲壓制。

    1.2 小波變換

    對(duì)于時(shí)頻信號(hào)分析,小波變換是一種較好的方法,它具有信噪分離、提取弱信號(hào)特征、多分辨率的特性[11],可以避免因傅氏變換出現(xiàn)的頻泄現(xiàn)象,并且在時(shí)域和頻域分析都有較好的效果,所以對(duì)于地震波這種非穩(wěn)定信號(hào)的分析處理中可以使用小波變換。

    對(duì)于任何的地震信號(hào)f(t),小波變換都可以將其表示為小波函數(shù)和地震信號(hào)f(t)內(nèi)積:

    (1)

    在信號(hào)處理中常采用Mallat算法實(shí)現(xiàn)對(duì)信號(hào)的小波變換和信號(hào)重構(gòu),即將a尺度空間的剩余系數(shù)dj,k,經(jīng)過(guò)濾波器系數(shù)h0(n)、h1(n)加權(quán)求和就得到j(luò)+1長(zhǎng)度空間的剩余系數(shù)dj+1,k和小波系數(shù)Cj+1,k。可以表示為:

    dj+1,k=∑mh0(m-2k)dj,m

    (2)

    Cj+1,k=∑mh1(m-2k)Cj,m

    (3)

    重構(gòu)公式為:

    Cj,n=∑mCj+1,kh1(m-2k)+

    ∑mdj+1,kh0(m-2k)

    (4)

    式中:h0、h1為濾波器系數(shù);dj+1,k、Cj+1,k為剩余系數(shù)和小波系數(shù)。

    從重構(gòu)信號(hào)可見(jiàn),選取小波系數(shù)和濾波手段,是壓制噪聲的關(guān)鍵。選用什么樣的小波基函數(shù)和處理等手段是決定小波降噪效果的關(guān)鍵。

    1.3 CEEMDAN原理

    EMD(經(jīng)驗(yàn)?zāi)B(tài)分解)是一種分析非平穩(wěn)非線性的自適應(yīng)方法[12],但是容易導(dǎo)致模態(tài)混淆問(wèn)題,所以Huang[13]提出了EEMD(集合經(jīng)驗(yàn)?zāi)B(tài)分解),通過(guò)添加白噪聲來(lái)讓EMD導(dǎo)致的模態(tài)混淆問(wèn)題減少,然而在EEMD重構(gòu)信號(hào)的同時(shí)會(huì)存在摻雜的噪聲。在EEMD的基礎(chǔ)上Torres[14]提出了CEEMDAN(自適應(yīng)噪聲完備集合經(jīng)驗(yàn)?zāi)B(tài)分解),CEEMDAN其實(shí)是在EEMD的基礎(chǔ)上進(jìn)行補(bǔ)充的一種信號(hào)處理方法。其主要在EEMD的算法上加入了有限次的自適應(yīng)白噪聲,在分解的每一個(gè)階段都添加自適應(yīng)白噪聲,再通過(guò)計(jì)算獲取每個(gè)IMF在比較少的次數(shù)下實(shí)現(xiàn)重構(gòu)地震信號(hào)并且達(dá)到幾乎沒(méi)有誤差,保證了分解的完整性同時(shí)也減少了計(jì)算的時(shí)間。

    CEEMDAN分解如下:

    1)通過(guò)對(duì)原始信號(hào)S(t)在添加不同的白噪聲ωi(t),原信號(hào)S(t)變?yōu)閤(t)+ε0ωi(t),ε為噪聲系數(shù)。用EMD獲取原始信號(hào)S(t)分解的第一個(gè)IMF分量公式表示為:

    E1[x(t)+ε0ωi(t) ]

    (5)

    2)第一個(gè)計(jì)算殘差信號(hào)r1為:

    r1(t)=s(t)-IMF1[t]

    (6)

    將r1(t)和經(jīng)過(guò)EMD分解的E1[ωi(t)]噪聲分量相加,然后在對(duì)其結(jié)果進(jìn)行EMD分解,即可以得到第二個(gè)IMF分量即:

    (7)

    式中:r1(t)為第一個(gè)殘差分量;E1[ωi(t)]為信號(hào)通過(guò)EMD分解得到的第一個(gè)模態(tài)函數(shù)。

    3)當(dāng)IMFk[t],K=2、3、…、k時(shí),計(jì)算K的殘差便重復(fù)以上過(guò)程。

    4)直到殘差不能在進(jìn)行分解便得到K個(gè)模態(tài)函數(shù)分量,其中K為完全自適應(yīng)無(wú)需人為設(shè)置,最后的結(jié)果R(t)為式(8)。

    (8)

    5)重構(gòu)信號(hào)x(t):

    (9)

    從上述步驟可以看出,CEEMDAN的分解能準(zhǔn)確地重構(gòu)目標(biāo)信號(hào),克服了EMD分解導(dǎo)致的模態(tài)混淆,有效地獲得了分離譜,計(jì)算效率也較大地提高。

    2 小波變換-CEEMDAN濾波

    2.1 高頻IMF分量的定位和選取

    當(dāng)使用傳統(tǒng)的EMD類(lèi)方法用于去噪時(shí),根據(jù)以往的知識(shí),即隨機(jī)噪聲信號(hào)通常分布在高頻分量中,則直接丟棄最前的1個(gè)到2個(gè)高頻IMF分量。 再用剩下的IMF分量重構(gòu)獲得去噪后的信號(hào)。但是依靠經(jīng)驗(yàn)選擇幾個(gè)高頻IMF分量會(huì)都受主觀因素的影響,高頻分量舍棄過(guò)多或不足都會(huì)導(dǎo)致去噪效果的不佳。

    故先通過(guò)CEEMDAN自適應(yīng)分解得到最優(yōu)的分解層數(shù),然后在計(jì)算高頻細(xì)節(jié)系數(shù)與分解的各個(gè)IMF分量之間的相關(guān)性系數(shù)的方法來(lái)確定高頻IMF分量的準(zhǔn)確位置,相關(guān)性系數(shù)大的IMF分量為地震信號(hào)中包含的噪聲干擾將其剔除,其他的分量則為有效信號(hào)將其保留,從而實(shí)現(xiàn)噪聲的準(zhǔn)確定位。

    相關(guān)系數(shù)計(jì)算方法:

    設(shè)Si(n)為第i個(gè)IMF分量Ui(n)與X(n)的相關(guān)系數(shù),則Si(n)可以表達(dá)為:

    (10)

    式中:Ui(n)、X(n)為IMF分量Ui(t)和高頻細(xì)節(jié)信號(hào)X(t)的離散化表示,n=1、2、3、…、m。

    在大多數(shù)情況下求得的相關(guān)系數(shù)的取值范圍應(yīng)在0~1之間,但通過(guò)實(shí)驗(yàn)發(fā)現(xiàn),個(gè)別的IMF分量在與高頻細(xì)節(jié)信號(hào)求得的相關(guān)系數(shù)會(huì)出現(xiàn)小于“0”的情況,由于欲通過(guò)各階IMF的相關(guān)系數(shù)來(lái)反應(yīng)噪聲信號(hào)在信號(hào)中的變化趨勢(shì),所以當(dāng)相關(guān)系數(shù)為負(fù)值時(shí),對(duì)其取絕對(duì)值。

    在剔除IMF分量的選取上,第1階IMF分量到第i階IMF分量Ui(n)和小波高頻細(xì)節(jié)系數(shù)X(n)之間的相關(guān)系數(shù)Si(n),總體來(lái)說(shuō)是呈現(xiàn)單調(diào)遞減的函數(shù)關(guān)系,且它們各點(diǎn)之間的斜率穩(wěn)定且大體不會(huì)發(fā)生較大的變化。直到第i+1階IMF分量時(shí)斜率出現(xiàn)折點(diǎn)且驟然變小之后趨于平滑。這說(shuō)明在第1階到i階,IMF分量中主要為噪聲信號(hào)。從第i+1階開(kāi)始有效信號(hào)為主要信號(hào)。此時(shí)剔除前i個(gè)IMF分量,從第i+1階開(kāi)始重構(gòu)信號(hào)。

    2.2 地震信號(hào)噪聲壓制處理流程

    小波變換降噪方式種類(lèi)太多,大都采用軟閥值、硬閥子。對(duì)此閥子的選取對(duì)于信號(hào)的降噪表現(xiàn)有很大影響,如果選擇不好,要么消噪太多,要么就是不夠,所以選取一個(gè)合適的閥子便顯得格外重要。為了擺脫這種問(wèn)題筆者將小波分解和CEEMDAN相結(jié)合對(duì)地震信號(hào)進(jìn)行降噪。整體降噪實(shí)現(xiàn)步驟如下(圖1):

    1)對(duì)地震原始信號(hào)S(t)進(jìn)行小波分解,將分解后的每一層高頻細(xì)節(jié)系數(shù)進(jìn)行CEEMDAN分解。

    2)使用CEEMDAN對(duì)高頻細(xì)節(jié)系數(shù)分解出的IMF分量,從頻率高低進(jìn)行排列。并計(jì)算出各個(gè)IMF分量與經(jīng)過(guò)小波分解后的高頻系數(shù)的相關(guān)性系數(shù)。

    3)對(duì)求出的相關(guān)性系數(shù),分析系數(shù)之間的線性關(guān)系。

    4)通過(guò)相關(guān)性分析,選擇斜率較大的數(shù)據(jù)進(jìn)行剔除,將保留下來(lái)的IMF分量進(jìn)行重構(gòu)得到處理后的小波高頻分量。

    5)對(duì)剩下的小波高頻分量重復(fù)步驟2)、步驟3)、步驟4)。

    6)重構(gòu)原始信號(hào)。將所有處理后的小波高頻細(xì)節(jié)系數(shù)與低頻近似系數(shù)進(jìn)行重構(gòu),得到去噪信號(hào)。

    圖1 小波—CEEMDAN降噪模型Fig.1 Wavelet-CEEMDAN noise reduction model

    圖2 地震原始信號(hào)Fig.2 The original seismic signal

    圖3 小波分解的高、低頻部分Fig.3 High and low frequency part of wavelet decomposition(a)第一層高頻細(xì)節(jié)系數(shù);(b)第二層高頻細(xì)節(jié)系數(shù);(c)第三層高頻細(xì)節(jié)系數(shù);(d)第三層低頻近似系數(shù)

    3 實(shí)際應(yīng)用

    為了體現(xiàn)本方法的可行性,筆者選用了地震實(shí)測(cè)數(shù)據(jù),數(shù)據(jù)來(lái)源于2018年南昌市多要素地質(zhì)調(diào)查項(xiàng)目。圖2是原始地震單道信號(hào),可以看出在單道信號(hào)中夾雜者較多噪聲干擾,不能明顯地分辨有效反射波的準(zhǔn)確位置。選擇bior2.4小波基[15]對(duì)其進(jìn)行小波分解。圖3是經(jīng)過(guò)小波分解后得到的地震低頻近似信號(hào)和各層的高頻細(xì)節(jié)信號(hào)。經(jīng)過(guò)小波分解后的高頻細(xì)節(jié)部分大都存在噪聲,所以在對(duì)前三層高頻信號(hào)進(jìn)行CEEMDAN分解,圖4是第一層高頻細(xì)節(jié)信號(hào)經(jīng)過(guò)CEEMDAN后自適應(yīng)分解得到的12個(gè)IMF分量。通過(guò)對(duì)分解得到的12個(gè)IMF分量求出各階自相關(guān)函數(shù)和相關(guān)性系數(shù),準(zhǔn)確定位高頻IMF分量。

    圖5為各階IMF自相關(guān)函數(shù)。通過(guò)圖5的觀察發(fā)現(xiàn)IMF1的自相關(guān)函數(shù)圖像有用信號(hào)完全被噪聲淹沒(méi)。大量的噪聲存在于IMF1分量中,而IMF2-MIF4分量中主要為有效信號(hào),干擾噪聲占比很小。表1為各階IMF與高頻細(xì)節(jié)系數(shù)的相關(guān)性系數(shù)。從表1可以看出IMF1的相關(guān)性系數(shù)高達(dá)0.986 4,此為極高相關(guān)分量。而之后分量的相關(guān)性系數(shù)大都比IMF1低2個(gè)數(shù)量級(jí)為極低相關(guān)分量。反應(yīng)出大量的噪聲信號(hào)存在于IMF1中。圖6為相關(guān)性系數(shù)之間的線性關(guān)系。通過(guò)圖6可以更加直觀地反應(yīng)出各相關(guān)性系數(shù)之間的關(guān)系,IMF1為要剔除的高頻分量。在對(duì)剩下的IMF余量進(jìn)行相加在與小波低頻信號(hào)進(jìn)行重構(gòu),得到處理后信號(hào)見(jiàn)圖7。從圖7上看,小波-CEEMDAN濾波取得了很好效果,對(duì)原始信號(hào)的干擾噪聲有了有效的壓制,突出了有效反射波,原始信號(hào)的主要變化特征也得到了很好的保留。

    圖4 CEEMDAN分解結(jié)果Fig.4 CEEMDAN decomposition results

    圖5 各階IMF自相關(guān)函數(shù)Fig.5 Autocorrelation function of each IMF

    表1 各階IMF與高頻細(xì)節(jié)系數(shù)的相關(guān)性系數(shù)

    圖6 相關(guān)性系數(shù)之間的線性關(guān)系Fig.6 Linear relationship between correlation coefficients

    圖7 小波-CEEMDAN濾波后信號(hào)Fig.7 Wavelet-CEEMDAN filtered signal

    為了體現(xiàn)小波-CEEMDAN濾波的效果,筆者對(duì)信號(hào)進(jìn)行信噪比(SNR)和均方差(RMSE)的量化比較評(píng)價(jià):信噪比數(shù)值越大則信號(hào)濾波效果越好,均方差數(shù)值越小濾波效果越好。

    SNR=10*lg{∑x(t)2/∑[x(t)-x(t)′]2}

    (10)

    (11)

    式中:x(t)、x(t)′分別為原始信號(hào)和去噪后信號(hào)。

    由式(10)、式(11)計(jì)算的信噪比和均方差結(jié)果見(jiàn)表2。

    表2 信噪比和均方差結(jié)果

    圖8 使用小波-CEEMDAN濾波地震剖面圖Fig.8 Seismic profile using wavelet- CEEMDAN filtering

    圖9 未使用小波-CEEMDAN濾波地震剖面圖Fig.9 Seismic profile without wavelet- CEEMDAN filtering

    從表2可以看出,使用小波-CEEMDAN相結(jié)合的濾波方法,比單獨(dú)使用CEEMDAN進(jìn)行濾波效果更好,大大提高了地震信號(hào)的信噪比,也降低了信號(hào)的均方差數(shù)值。因此可以說(shuō)明地震信號(hào)通過(guò)小波-CEEMDAN濾波,在有效地去除信號(hào)的噪聲信號(hào)的同時(shí),也極大地保留了原始信號(hào)的主要特征。使用小波-CEEMDAN壓制噪聲效果要優(yōu)于單獨(dú)使用CEEMDAN方法。所以我們對(duì)整個(gè)剖面的地震信號(hào)進(jìn)行小波-CEEMDAN濾波,濾波前后結(jié)果如圖8、圖9所示,可以看出,地震信號(hào)使用波-CEEMDAN濾波后,T0波組連貫性更好,更加平滑,橫向分辨率也大大增加。T1波組反射波信號(hào)得到加強(qiáng),反射界面更加突出。而未使用波-CEEMDAN濾波的地震信號(hào)T0波組地震信號(hào)橫向分辨率較弱,連續(xù)性較弱。T1波地震反射信號(hào)較弱,連續(xù)性不好。從圖8和圖9的對(duì)比中可以看出,地震信號(hào)經(jīng)過(guò)小波-CEEMDAN濾波后信號(hào)能量得到加強(qiáng),消除了同相軸周邊的噪聲干擾,信噪比也大大提高,橫向分辨率更好。

    4 結(jié)論

    1)對(duì)小波高頻細(xì)節(jié)系數(shù)使用CEEMDAN分解,再通過(guò)IMF分量算取各自信號(hào)的自相關(guān)函數(shù),能很好地對(duì)高頻系數(shù)進(jìn)行噪聲去除,使用這種方法重構(gòu)的小波信號(hào)能在提高數(shù)據(jù)擬合度的同時(shí)對(duì)信號(hào)的噪聲進(jìn)行壓制。

    2)基于小波-CEEMDAN處理的非平穩(wěn)非線性地震信號(hào)比單獨(dú)使用CEEMDAN處理的地震信號(hào)效果更好,使用小波-CEEMDAN濾波相結(jié)合的方式既可以避免小波去噪多種方式的選擇和去噪閾值函數(shù)的選擇,還可以最大限度的保留原始信號(hào)的特征特點(diǎn),避免信號(hào)的失真,使處理后的信號(hào)更加真實(shí),為地震信號(hào)噪聲壓制提供一種方法。

    3)對(duì)于在強(qiáng)干擾地區(qū)的實(shí)測(cè)城市淺層地震記錄經(jīng)過(guò)小波-CEEMDAN濾波處理后可以看出,其對(duì)噪聲壓制取得了比較滿意的效果,并且準(zhǔn)確地反映了有效反射波的位置,數(shù)據(jù)的信噪比也有明顯地提升,橫向分辨率提高為以后的解釋研究工作也提供了便利。

    猜你喜歡
    小波分量重構(gòu)
    長(zhǎng)城敘事的重構(gòu)
    攝影世界(2022年1期)2022-01-21 10:50:14
    構(gòu)造Daubechies小波的一些注記
    帽子的分量
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    一物千斤
    智族GQ(2019年9期)2019-10-28 08:16:21
    北方大陸 重構(gòu)未來(lái)
    論《哈姆雷特》中良心的分量
    基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動(dòng)軸承故障診斷
    北京的重構(gòu)與再造
    商周刊(2017年6期)2017-08-22 03:42:36
    分量
    看十八女毛片水多多多| 欧美中文综合在线视频| 午夜福利影视在线免费观看| 人人妻人人澡人人看| 亚洲一区二区三区欧美精品| 看十八女毛片水多多多| 18禁国产床啪视频网站| 国产亚洲av片在线观看秒播厂| 亚洲成人一二三区av| 咕卡用的链子| 欧美精品一区二区免费开放| 日韩在线高清观看一区二区三区| 啦啦啦啦在线视频资源| 久久久久网色| 大片免费播放器 马上看| 性少妇av在线| 这个男人来自地球电影免费观看 | 我要看黄色一级片免费的| 大片电影免费在线观看免费| 色网站视频免费| 91成人精品电影| 亚洲精品在线美女| 久久久久久久精品精品| 亚洲欧洲国产日韩| 国产免费福利视频在线观看| 亚洲经典国产精华液单| 亚洲成av片中文字幕在线观看 | av视频免费观看在线观看| 亚洲av在线观看美女高潮| 一级黄片播放器| av电影中文网址| 亚洲 欧美一区二区三区| 18禁裸乳无遮挡动漫免费视频| 高清黄色对白视频在线免费看| 尾随美女入室| 新久久久久国产一级毛片| 少妇人妻精品综合一区二区| 久久国产亚洲av麻豆专区| 国产亚洲欧美精品永久| 制服诱惑二区| 人成视频在线观看免费观看| 国产欧美亚洲国产| 男人添女人高潮全过程视频| 侵犯人妻中文字幕一二三四区| 免费观看无遮挡的男女| 校园人妻丝袜中文字幕| 亚洲内射少妇av| 亚洲一区二区三区欧美精品| 久久鲁丝午夜福利片| 国产又爽黄色视频| 新久久久久国产一级毛片| 亚洲av在线观看美女高潮| 侵犯人妻中文字幕一二三四区| 国语对白做爰xxxⅹ性视频网站| 卡戴珊不雅视频在线播放| 男女边吃奶边做爰视频| 国产成人精品久久二区二区91 | 青青草视频在线视频观看| 多毛熟女@视频| 色婷婷av一区二区三区视频| 成人毛片a级毛片在线播放| av有码第一页| 亚洲久久久国产精品| 亚洲第一青青草原| 国产精品不卡视频一区二区| 精品国产一区二区三区四区第35| 精品少妇久久久久久888优播| 这个男人来自地球电影免费观看 | 妹子高潮喷水视频| 91精品国产国语对白视频| 在线观看国产h片| 夫妻午夜视频| 成人国产av品久久久| 国产一区二区激情短视频 | 麻豆乱淫一区二区| 国产在视频线精品| 啦啦啦视频在线资源免费观看| 男女边吃奶边做爰视频| 女人精品久久久久毛片| 七月丁香在线播放| 欧美精品一区二区大全| 九草在线视频观看| 日韩欧美一区视频在线观看| 在线观看免费视频网站a站| 精品一区二区三卡| 久久久久久人人人人人| 三上悠亚av全集在线观看| av线在线观看网站| 日本91视频免费播放| 超碰97精品在线观看| 亚洲国产精品一区三区| 又粗又硬又长又爽又黄的视频| 男人操女人黄网站| 视频区图区小说| 精品午夜福利在线看| 欧美少妇被猛烈插入视频| 丰满乱子伦码专区| 卡戴珊不雅视频在线播放| 久久精品夜色国产| 新久久久久国产一级毛片| 免费在线观看完整版高清| 国产有黄有色有爽视频| 日韩中字成人| 日本猛色少妇xxxxx猛交久久| 99久久综合免费| 国产成人aa在线观看| 香蕉精品网在线| 日韩精品免费视频一区二区三区| 欧美日本中文国产一区发布| 亚洲精品aⅴ在线观看| 亚洲欧洲国产日韩| 国产亚洲最大av| 在线观看免费高清a一片| 丝袜美足系列| 侵犯人妻中文字幕一二三四区| 狠狠精品人妻久久久久久综合| 日韩欧美一区视频在线观看| 大陆偷拍与自拍| 亚洲天堂av无毛| 国产精品三级大全| 91精品伊人久久大香线蕉| 亚洲欧美精品综合一区二区三区 | 精品人妻在线不人妻| 女人久久www免费人成看片| 欧美亚洲 丝袜 人妻 在线| 色吧在线观看| 91在线精品国自产拍蜜月| 9热在线视频观看99| 青春草国产在线视频| 精品亚洲乱码少妇综合久久| 老女人水多毛片| 精品一品国产午夜福利视频| 人妻一区二区av| 18禁国产床啪视频网站| a级片在线免费高清观看视频| 80岁老熟妇乱子伦牲交| av视频免费观看在线观看| 捣出白浆h1v1| 国产男女超爽视频在线观看| 制服丝袜香蕉在线| 永久网站在线| 日韩精品免费视频一区二区三区| 久久国产精品男人的天堂亚洲| 亚洲第一av免费看| 久久这里有精品视频免费| 久久人人爽人人片av| 国语对白做爰xxxⅹ性视频网站| 一区二区三区乱码不卡18| 亚洲精品国产av成人精品| 婷婷色综合www| 国产精品嫩草影院av在线观看| 丝袜脚勾引网站| 精品国产国语对白av| 大片免费播放器 马上看| 啦啦啦中文免费视频观看日本| 久久久久人妻精品一区果冻| 男女边吃奶边做爰视频| 国产片内射在线| 国产视频首页在线观看| 国产亚洲午夜精品一区二区久久| 又粗又硬又长又爽又黄的视频| 熟女少妇亚洲综合色aaa.| 在线观看三级黄色| 久久久国产一区二区| 少妇的丰满在线观看| 国产精品人妻久久久影院| 大码成人一级视频| 大码成人一级视频| 日本猛色少妇xxxxx猛交久久| 欧美中文综合在线视频| 激情五月婷婷亚洲| 成人影院久久| 国产人伦9x9x在线观看 | 欧美av亚洲av综合av国产av | 久久精品国产鲁丝片午夜精品| 国产日韩欧美亚洲二区| 另类精品久久| 在线亚洲精品国产二区图片欧美| 三上悠亚av全集在线观看| 丝袜喷水一区| 在线观看免费视频网站a站| 精品一区二区三卡| 亚洲人成网站在线观看播放| 男人添女人高潮全过程视频| 亚洲精品美女久久av网站| 制服丝袜香蕉在线| 久久精品国产亚洲av高清一级| 免费观看av网站的网址| 国产极品天堂在线| 久久久久视频综合| 亚洲国产精品一区三区| av天堂久久9| 大片电影免费在线观看免费| 成人18禁高潮啪啪吃奶动态图| 热re99久久国产66热| 久久久久国产精品人妻一区二区| 久久久久久久国产电影| 一个人免费看片子| 极品少妇高潮喷水抽搐| 午夜福利视频在线观看免费| 欧美97在线视频| 伦精品一区二区三区| 在线天堂最新版资源| 免费人妻精品一区二区三区视频| 成人黄色视频免费在线看| 国产熟女欧美一区二区| 精品少妇黑人巨大在线播放| 天堂中文最新版在线下载| 欧美日韩精品成人综合77777| 国产成人精品久久二区二区91 | 老汉色∧v一级毛片| h视频一区二区三区| 国产av一区二区精品久久| 在线天堂中文资源库| 肉色欧美久久久久久久蜜桃| 日韩av在线免费看完整版不卡| 最新中文字幕久久久久| 国产精品亚洲av一区麻豆 | 成人手机av| 制服诱惑二区| 菩萨蛮人人尽说江南好唐韦庄| 国产精品亚洲av一区麻豆 | 中文字幕色久视频| 国产片内射在线| 啦啦啦在线免费观看视频4| 最近手机中文字幕大全| 一二三四在线观看免费中文在| 亚洲精品中文字幕在线视频| 极品少妇高潮喷水抽搐| 精品人妻在线不人妻| 亚洲第一青青草原| 国产免费现黄频在线看| 母亲3免费完整高清在线观看 | 大码成人一级视频| 2018国产大陆天天弄谢| av在线播放精品| 美国免费a级毛片| 中文字幕制服av| 国产成人午夜福利电影在线观看| 人人妻人人澡人人爽人人夜夜| av.在线天堂| 最新中文字幕久久久久| 亚洲av电影在线观看一区二区三区| 国产精品熟女久久久久浪| 如何舔出高潮| 极品少妇高潮喷水抽搐| 欧美97在线视频| 成人亚洲精品一区在线观看| 永久网站在线| 日韩不卡一区二区三区视频在线| 97在线视频观看| 人人妻人人添人人爽欧美一区卜| 国语对白做爰xxxⅹ性视频网站| 女的被弄到高潮叫床怎么办| 秋霞伦理黄片| 久久久国产精品麻豆| 亚洲精品日韩在线中文字幕| 亚洲综合色惰| 搡老乐熟女国产| www.熟女人妻精品国产| 免费观看无遮挡的男女| av不卡在线播放| 成年av动漫网址| 亚洲第一av免费看| 热99久久久久精品小说推荐| 国产av精品麻豆| 国产精品一国产av| √禁漫天堂资源中文www| 99久久综合免费| 日韩中文字幕欧美一区二区 | 制服丝袜香蕉在线| 亚洲美女视频黄频| 国产片内射在线| 91成人精品电影| 蜜桃在线观看..| 国产av精品麻豆| 街头女战士在线观看网站| 免费不卡的大黄色大毛片视频在线观看| 欧美中文综合在线视频| 亚洲成人手机| 一边摸一边做爽爽视频免费| 国产亚洲精品第一综合不卡| 国产在线免费精品| 国产成人精品久久二区二区91 | 伦精品一区二区三区| 男女啪啪激烈高潮av片| 大香蕉久久成人网| 国产精品亚洲av一区麻豆 | 有码 亚洲区| 日本欧美国产在线视频| 久久久久久久亚洲中文字幕| 国产白丝娇喘喷水9色精品| 国产一区二区三区av在线| 亚洲欧洲日产国产| 最近2019中文字幕mv第一页| 国产精品香港三级国产av潘金莲 | 2018国产大陆天天弄谢| 亚洲,一卡二卡三卡| 老司机亚洲免费影院| 亚洲国产精品国产精品| 久久精品国产亚洲av涩爱| 免费人妻精品一区二区三区视频| 亚洲内射少妇av| 亚洲欧美一区二区三区久久| 三上悠亚av全集在线观看| 80岁老熟妇乱子伦牲交| 黄频高清免费视频| 久久精品国产自在天天线| 日日摸夜夜添夜夜爱| 18禁动态无遮挡网站| 黄片小视频在线播放| 成人免费观看视频高清| 欧美精品一区二区免费开放| 欧美精品国产亚洲| 我要看黄色一级片免费的| 丝袜在线中文字幕| 国产av一区二区精品久久| 免费人妻精品一区二区三区视频| 亚洲色图综合在线观看| 国产日韩欧美亚洲二区| 久久久久久久久久久免费av| 亚洲天堂av无毛| 色视频在线一区二区三区| 97精品久久久久久久久久精品| 精品一区二区三卡| 久久精品国产亚洲av高清一级| 99久久中文字幕三级久久日本| 午夜日韩欧美国产| 在线亚洲精品国产二区图片欧美| 国产成人一区二区在线| freevideosex欧美| 国产亚洲欧美精品永久| 亚洲欧美清纯卡通| 日日啪夜夜爽| 国产免费现黄频在线看| 最黄视频免费看| 天美传媒精品一区二区| 人妻系列 视频| 欧美变态另类bdsm刘玥| 国产成人精品无人区| 18禁裸乳无遮挡动漫免费视频| 精品少妇内射三级| 男人舔女人的私密视频| 色94色欧美一区二区| 一级毛片 在线播放| 亚洲精品久久午夜乱码| 国产无遮挡羞羞视频在线观看| 秋霞伦理黄片| 老鸭窝网址在线观看| 久久久精品94久久精品| 人妻系列 视频| 国产欧美亚洲国产| 精品国产乱码久久久久久小说| 高清不卡的av网站| www.自偷自拍.com| 亚洲中文av在线| 亚洲欧美成人综合另类久久久| 午夜日韩欧美国产| 91成人精品电影| 久久久久久久久久久免费av| 国产极品粉嫩免费观看在线| 一区二区三区四区激情视频| 久久久久久久国产电影| 国产有黄有色有爽视频| 高清av免费在线| 亚洲精品第二区| 亚洲精华国产精华液的使用体验| 婷婷色麻豆天堂久久| 在线观看美女被高潮喷水网站| 99热全是精品| 看免费成人av毛片| 国产一区二区在线观看av| 国产在视频线精品| 天天影视国产精品| 国产精品 国内视频| 欧美另类一区| 亚洲精品日本国产第一区| 国产一区有黄有色的免费视频| 黄色一级大片看看| 亚洲国产av影院在线观看| 一本大道久久a久久精品| 99久久人妻综合| 日韩熟女老妇一区二区性免费视频| 少妇人妻久久综合中文| 99久国产av精品国产电影| 久久婷婷青草| 国产精品国产三级专区第一集| 日本av手机在线免费观看| 欧美成人午夜免费资源| 欧美日韩一区二区视频在线观看视频在线| 可以免费在线观看a视频的电影网站 | videosex国产| 高清av免费在线| 麻豆乱淫一区二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91国产中文字幕| 精品少妇久久久久久888优播| 国产精品蜜桃在线观看| 日韩熟女老妇一区二区性免费视频| 久久这里只有精品19| 国产在视频线精品| 久久国产精品男人的天堂亚洲| 永久免费av网站大全| 久久久国产一区二区| 天堂俺去俺来也www色官网| 久久久久久久久久久免费av| 在线观看三级黄色| 青春草国产在线视频| av有码第一页| 久久久久精品久久久久真实原创| 久久 成人 亚洲| 午夜福利视频在线观看免费| 国产成人精品久久久久久| 国产人伦9x9x在线观看 | 午夜日本视频在线| 最近中文字幕高清免费大全6| 日日啪夜夜爽| 午夜福利网站1000一区二区三区| 欧美日韩视频高清一区二区三区二| 97在线人人人人妻| 婷婷色麻豆天堂久久| 91精品伊人久久大香线蕉| 新久久久久国产一级毛片| 久久精品人人爽人人爽视色| 国产一区二区在线观看av| 老司机影院成人| 亚洲第一青青草原| 亚洲国产av影院在线观看| 97人妻天天添夜夜摸| 成年人免费黄色播放视频| 亚洲精品乱久久久久久| 少妇人妻精品综合一区二区| 亚洲中文av在线| 国产精品三级大全| 久久久久久久大尺度免费视频| 亚洲成色77777| 一本—道久久a久久精品蜜桃钙片| 人妻系列 视频| 丝袜美腿诱惑在线| 9色porny在线观看| 精品亚洲成国产av| 最新的欧美精品一区二区| 在线观看国产h片| 晚上一个人看的免费电影| 最近中文字幕2019免费版| 欧美少妇被猛烈插入视频| 国产深夜福利视频在线观看| 1024视频免费在线观看| 久久午夜综合久久蜜桃| 久热久热在线精品观看| 欧美xxⅹ黑人| 国精品久久久久久国模美| kizo精华| 国产97色在线日韩免费| 中文欧美无线码| 在线观看国产h片| 亚洲伊人色综图| 久久久亚洲精品成人影院| xxxhd国产人妻xxx| 亚洲国产毛片av蜜桃av| 日产精品乱码卡一卡2卡三| 美女福利国产在线| 久久婷婷青草| 人人妻人人澡人人看| 人妻少妇偷人精品九色| 久久久国产精品麻豆| 亚洲欧洲精品一区二区精品久久久 | 纵有疾风起免费观看全集完整版| 久久久久久久久免费视频了| 高清在线视频一区二区三区| 日韩一本色道免费dvd| 热re99久久国产66热| 尾随美女入室| 久久久国产一区二区| 亚洲av在线观看美女高潮| 91aial.com中文字幕在线观看| 亚洲国产色片| 在线天堂中文资源库| 熟女少妇亚洲综合色aaa.| 晚上一个人看的免费电影| a级毛片在线看网站| 亚洲精品国产av成人精品| 久久鲁丝午夜福利片| 国产亚洲午夜精品一区二区久久| 成年av动漫网址| 一级,二级,三级黄色视频| 黄色毛片三级朝国网站| 久久国产精品男人的天堂亚洲| 亚洲精品久久久久久婷婷小说| 国产免费一区二区三区四区乱码| 亚洲,欧美精品.| 日本猛色少妇xxxxx猛交久久| 男女午夜视频在线观看| av免费观看日本| 色婷婷久久久亚洲欧美| 大香蕉久久网| 国产高清国产精品国产三级| 国产精品麻豆人妻色哟哟久久| 亚洲美女黄色视频免费看| 9热在线视频观看99| 成人亚洲欧美一区二区av| 青青草视频在线视频观看| 亚洲精品久久成人aⅴ小说| av免费在线看不卡| 亚洲国产成人一精品久久久| 亚洲激情五月婷婷啪啪| 一本大道久久a久久精品| 国产成人一区二区在线| 日韩三级伦理在线观看| 999精品在线视频| 视频在线观看一区二区三区| 色播在线永久视频| 国产探花极品一区二区| 97在线视频观看| 91成人精品电影| 人人妻人人澡人人看| 热re99久久精品国产66热6| 久久免费观看电影| 99久国产av精品国产电影| 波野结衣二区三区在线| 精品少妇久久久久久888优播| 国产精品国产av在线观看| 欧美精品国产亚洲| videos熟女内射| 免费播放大片免费观看视频在线观看| 韩国高清视频一区二区三区| 成人毛片a级毛片在线播放| av电影中文网址| 久久久久国产网址| 久久影院123| 亚洲精品久久成人aⅴ小说| 一本色道久久久久久精品综合| 久热久热在线精品观看| 国产黄色免费在线视频| 五月天丁香电影| 国产精品av久久久久免费| 女性被躁到高潮视频| 看免费成人av毛片| 国产日韩一区二区三区精品不卡| 亚洲国产看品久久| 国产成人免费观看mmmm| 久久久久久人妻| 两性夫妻黄色片| 亚洲精品成人av观看孕妇| 最近最新中文字幕大全免费视频 | 亚洲精华国产精华液的使用体验| 国产男人的电影天堂91| 女人高潮潮喷娇喘18禁视频| 久久精品国产自在天天线| 国产无遮挡羞羞视频在线观看| 男女啪啪激烈高潮av片| 在现免费观看毛片| 丝袜美腿诱惑在线| 精品卡一卡二卡四卡免费| 不卡视频在线观看欧美| 亚洲av电影在线观看一区二区三区| 日韩精品免费视频一区二区三区| 亚洲一级一片aⅴ在线观看| 精品亚洲乱码少妇综合久久| 国产欧美日韩综合在线一区二区| 久久久久久免费高清国产稀缺| 亚洲精品日韩在线中文字幕| 欧美日韩av久久| 亚洲一码二码三码区别大吗| 一本大道久久a久久精品| 老司机影院毛片| 热re99久久国产66热| 精品酒店卫生间| 日日爽夜夜爽网站| 老汉色av国产亚洲站长工具| 美女xxoo啪啪120秒动态图| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频| 国产成人精品一,二区| 国产97色在线日韩免费| 天天躁狠狠躁夜夜躁狠狠躁| 七月丁香在线播放| av在线app专区| 少妇的丰满在线观看| 国产又爽黄色视频| 2018国产大陆天天弄谢| 久久这里有精品视频免费| 免费av中文字幕在线| 宅男免费午夜| 久久午夜综合久久蜜桃| 另类亚洲欧美激情| www.熟女人妻精品国产| 国产又色又爽无遮挡免| 青春草亚洲视频在线观看| 尾随美女入室| 日韩成人av中文字幕在线观看| 亚洲男人天堂网一区| 我的亚洲天堂| 精品人妻熟女毛片av久久网站| 最近中文字幕高清免费大全6| 精品人妻熟女毛片av久久网站| 黄色怎么调成土黄色| 欧美+日韩+精品| 啦啦啦中文免费视频观看日本| 伦理电影大哥的女人| 午夜福利视频精品| 亚洲人成77777在线视频| 亚洲欧美一区二区三区久久| 国产av码专区亚洲av| 国产精品嫩草影院av在线观看| 午夜福利在线观看免费完整高清在| 最近手机中文字幕大全| 人体艺术视频欧美日本| 欧美日韩国产mv在线观看视频| 欧美激情高清一区二区三区 | 中文字幕精品免费在线观看视频| 成人毛片a级毛片在线播放| 少妇 在线观看| 两性夫妻黄色片| 免费黄网站久久成人精品|