• <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
    分量
    亚洲情色 制服丝袜| 狂野欧美激情性xxxx| 亚洲人成伊人成综合网2020| 久久久精品国产亚洲av高清涩受| 岛国视频午夜一区免费看| 亚洲国产看品久久| 国产精品免费一区二区三区在线| 亚洲av日韩精品久久久久久密| 国产黄a三级三级三级人| АⅤ资源中文在线天堂| 国产一区二区激情短视频| 久久精品国产清高在天天线| 成人免费观看视频高清| 少妇裸体淫交视频免费看高清 | 嫩草影院精品99| 亚洲成人国产一区在线观看| 亚洲精品美女久久久久99蜜臀| 国产一区在线观看成人免费| 久久精品人人爽人人爽视色| 男人的好看免费观看在线视频 | 在线十欧美十亚洲十日本专区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产成人精品二区| 校园春色视频在线观看| 真人做人爱边吃奶动态| 在线十欧美十亚洲十日本专区| 一进一出抽搐gif免费好疼| 国产三级在线视频| 亚洲激情在线av| 19禁男女啪啪无遮挡网站| 日韩欧美免费精品| 这个男人来自地球电影免费观看| 99久久精品国产亚洲精品| 制服诱惑二区| 中出人妻视频一区二区| 人妻久久中文字幕网| 国产亚洲精品av在线| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美日韩另类电影网站| 久久久久久人人人人人| 亚洲片人在线观看| 国产99白浆流出| 一区二区三区激情视频| 日日夜夜操网爽| 色尼玛亚洲综合影院| 亚洲成a人片在线一区二区| 欧美乱妇无乱码| 大型av网站在线播放| 亚洲狠狠婷婷综合久久图片| 国产精品,欧美在线| 美女免费视频网站| 日韩中文字幕欧美一区二区| 激情在线观看视频在线高清| 亚洲成人久久性| 亚洲av美国av| 欧美激情高清一区二区三区| 亚洲久久久国产精品| 黄片小视频在线播放| 欧美精品亚洲一区二区| 91精品国产国语对白视频| 欧美国产日韩亚洲一区| 亚洲成a人片在线一区二区| 99久久国产精品久久久| 国产精华一区二区三区| 午夜亚洲福利在线播放| 黑丝袜美女国产一区| 18禁国产床啪视频网站| 亚洲第一av免费看| 黑人巨大精品欧美一区二区蜜桃| 成人av一区二区三区在线看| 国产一级毛片七仙女欲春2 | 亚洲av美国av| 欧美黑人欧美精品刺激| 少妇粗大呻吟视频| 亚洲精品一卡2卡三卡4卡5卡| 韩国av一区二区三区四区| 亚洲avbb在线观看| 夜夜夜夜夜久久久久| 99精品久久久久人妻精品| 可以在线观看的亚洲视频| 亚洲人成电影观看| 黄频高清免费视频| 日本精品一区二区三区蜜桃| 日本精品一区二区三区蜜桃| 91麻豆精品激情在线观看国产| 精品一品国产午夜福利视频| 国产三级在线视频| 久久久久亚洲av毛片大全| 日本免费一区二区三区高清不卡 | 妹子高潮喷水视频| 欧美日韩瑟瑟在线播放| 久久精品亚洲精品国产色婷小说| 一二三四在线观看免费中文在| 久热这里只有精品99| 宅男免费午夜| 欧美人与性动交α欧美精品济南到| 怎么达到女性高潮| 国产精品影院久久| 美女高潮到喷水免费观看| 欧美人与性动交α欧美精品济南到| 看免费av毛片| 亚洲国产精品999在线| 天天添夜夜摸| 夜夜爽天天搞| 亚洲第一电影网av| 波多野结衣巨乳人妻| www.自偷自拍.com| 亚洲熟妇熟女久久| 国产精品 国内视频| 一进一出抽搐动态| 黄色丝袜av网址大全| www.熟女人妻精品国产| 午夜日韩欧美国产| 无遮挡黄片免费观看| 国产精品亚洲一级av第二区| 欧美+亚洲+日韩+国产| 免费看十八禁软件| 夜夜夜夜夜久久久久| 国产1区2区3区精品| 757午夜福利合集在线观看| 午夜视频精品福利| av有码第一页| 午夜精品久久久久久毛片777| 久久国产精品人妻蜜桃| 中文亚洲av片在线观看爽| 久久久久亚洲av毛片大全| 亚洲第一av免费看| 满18在线观看网站| 久久久久久亚洲精品国产蜜桃av| 色综合欧美亚洲国产小说| 亚洲色图av天堂| 此物有八面人人有两片| 欧美激情极品国产一区二区三区| 日韩高清综合在线| 久久久久精品国产欧美久久久| 天堂影院成人在线观看| 久久精品亚洲精品国产色婷小说| 亚洲国产毛片av蜜桃av| 日韩欧美免费精品| 人妻久久中文字幕网| 久久影院123| 午夜两性在线视频| av欧美777| 欧美日韩亚洲综合一区二区三区_| 欧美老熟妇乱子伦牲交| 亚洲人成77777在线视频| 欧美成狂野欧美在线观看| 中文字幕色久视频| 黄色丝袜av网址大全| 久久人人爽av亚洲精品天堂| 亚洲天堂国产精品一区在线| 色综合欧美亚洲国产小说| 在线播放国产精品三级| 久久久久久久久免费视频了| 怎么达到女性高潮| 精品高清国产在线一区| 欧美久久黑人一区二区| 欧美成人免费av一区二区三区| 成年女人毛片免费观看观看9| 成人三级做爰电影| 伦理电影免费视频| 黑人欧美特级aaaaaa片| 看免费av毛片| 9191精品国产免费久久| 不卡av一区二区三区| 妹子高潮喷水视频| 亚洲av美国av| 国产色视频综合| 自拍欧美九色日韩亚洲蝌蚪91| 日日爽夜夜爽网站| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久精品电影 | 日韩欧美免费精品| 黄片大片在线免费观看| 国产精品 欧美亚洲| 热re99久久国产66热| 女同久久另类99精品国产91| 亚洲五月色婷婷综合| 1024香蕉在线观看| 19禁男女啪啪无遮挡网站| 亚洲人成77777在线视频| 久久人妻福利社区极品人妻图片| 狠狠狠狠99中文字幕| 18禁裸乳无遮挡免费网站照片 | 国产一区在线观看成人免费| 国产精品av久久久久免费| 啪啪无遮挡十八禁网站| 久久久久国内视频| 真人一进一出gif抽搐免费| 极品人妻少妇av视频| 欧美丝袜亚洲另类 | 他把我摸到了高潮在线观看| 久久精品国产亚洲av香蕉五月| 宅男免费午夜| 亚洲激情在线av| 欧美一区二区精品小视频在线| 麻豆久久精品国产亚洲av| 最好的美女福利视频网| 亚洲欧美日韩无卡精品| 精品人妻在线不人妻| 日本免费一区二区三区高清不卡 | 岛国视频午夜一区免费看| 村上凉子中文字幕在线| 波多野结衣高清无吗| 97人妻天天添夜夜摸| 亚洲精品粉嫩美女一区| 19禁男女啪啪无遮挡网站| 十八禁网站免费在线| 国产av一区二区精品久久| 亚洲片人在线观看| 99国产精品一区二区三区| 禁无遮挡网站| 黄色女人牲交| 久久久久九九精品影院| 美国免费a级毛片| 成人国语在线视频| 黄色视频,在线免费观看| 国产成人精品无人区| 久久久国产成人免费| 校园春色视频在线观看| 欧美中文综合在线视频| 校园春色视频在线观看| 免费观看精品视频网站| 夜夜躁狠狠躁天天躁| 国产精品亚洲av一区麻豆| 性少妇av在线| 亚洲男人的天堂狠狠| 9191精品国产免费久久| 亚洲成人免费电影在线观看| 黄色成人免费大全| 免费在线观看影片大全网站| 999久久久国产精品视频| 一a级毛片在线观看| 在线观看一区二区三区| 亚洲片人在线观看| 岛国视频午夜一区免费看| 动漫黄色视频在线观看| 日韩欧美在线二视频| 最新美女视频免费是黄的| 波多野结衣巨乳人妻| 老司机福利观看| 欧美久久黑人一区二区| 亚洲va日本ⅴa欧美va伊人久久| 免费无遮挡裸体视频| 18禁观看日本| 一区二区日韩欧美中文字幕| 久久久久久久午夜电影| 老鸭窝网址在线观看| 精品熟女少妇八av免费久了| 成人三级做爰电影| 黄色毛片三级朝国网站| 亚洲色图av天堂| 91字幕亚洲| 亚洲激情在线av| 亚洲性夜色夜夜综合| 久久久久久久久久久久大奶| 国产精华一区二区三区| 两人在一起打扑克的视频| 日韩一卡2卡3卡4卡2021年| av中文乱码字幕在线| 婷婷精品国产亚洲av在线| 国产91精品成人一区二区三区| 亚洲人成77777在线视频| 大型av网站在线播放| 一本久久中文字幕| 久久人妻熟女aⅴ| 久久性视频一级片| 自拍欧美九色日韩亚洲蝌蚪91| 制服人妻中文乱码| 国产精品野战在线观看| 午夜福利影视在线免费观看| 亚洲国产欧美日韩在线播放| 午夜福利18| 欧美激情 高清一区二区三区| 国产精品精品国产色婷婷| 久久天躁狠狠躁夜夜2o2o| 黄色片一级片一级黄色片| 国产片内射在线| 色播在线永久视频| 一区二区三区高清视频在线| 啦啦啦 在线观看视频| 18禁黄网站禁片午夜丰满| 亚洲av片天天在线观看| 女人高潮潮喷娇喘18禁视频| 久久婷婷人人爽人人干人人爱 | 亚洲国产中文字幕在线视频| 日韩有码中文字幕| 99国产精品99久久久久| 亚洲第一青青草原| 三级毛片av免费| 99久久99久久久精品蜜桃| 国产成人系列免费观看| 19禁男女啪啪无遮挡网站| 操出白浆在线播放| 在线十欧美十亚洲十日本专区| 999久久久国产精品视频| 国产午夜福利久久久久久| 一级a爱视频在线免费观看| 高清黄色对白视频在线免费看| 丝袜人妻中文字幕| a在线观看视频网站| 嫁个100分男人电影在线观看| 免费在线观看亚洲国产| 成人免费观看视频高清| 亚洲色图 男人天堂 中文字幕| 日韩欧美在线二视频| 不卡av一区二区三区| 午夜免费成人在线视频| 精品国产亚洲在线| 国产精品久久久av美女十八| 男女床上黄色一级片免费看| 人成视频在线观看免费观看| 别揉我奶头~嗯~啊~动态视频| 老司机午夜福利在线观看视频| 国产伦一二天堂av在线观看| 午夜两性在线视频| 麻豆久久精品国产亚洲av| 国产私拍福利视频在线观看| 中文字幕精品免费在线观看视频| 亚洲精品国产区一区二| 午夜福利,免费看| 51午夜福利影视在线观看| 香蕉丝袜av| 久久久国产精品麻豆| 中文字幕高清在线视频| 好看av亚洲va欧美ⅴa在| 国产精品综合久久久久久久免费 | 在线观看一区二区三区| 中国美女看黄片| 9色porny在线观看| 久久精品91蜜桃| 丝袜美腿诱惑在线| 亚洲天堂国产精品一区在线| 国产av在哪里看| 久久久久国内视频| 国产亚洲精品综合一区在线观看 | 美女国产高潮福利片在线看| 国产欧美日韩精品亚洲av| 可以在线观看的亚洲视频| 黄色a级毛片大全视频| 国产xxxxx性猛交| 18禁裸乳无遮挡免费网站照片 | 黄片小视频在线播放| www.www免费av| 国产麻豆成人av免费视频| 精品福利观看| 黄色视频不卡| 中文字幕色久视频| 亚洲欧美精品综合一区二区三区| 女性生殖器流出的白浆| 91大片在线观看| 丝袜在线中文字幕| 搡老妇女老女人老熟妇| 国产激情欧美一区二区| 欧美av亚洲av综合av国产av| av超薄肉色丝袜交足视频| 国产精品一区二区在线不卡| 亚洲欧美激情在线| 久久久久久国产a免费观看| 女警被强在线播放| 首页视频小说图片口味搜索| 亚洲精品一卡2卡三卡4卡5卡| 啦啦啦免费观看视频1| av网站免费在线观看视频| 性少妇av在线| 午夜免费成人在线视频| 国产精品久久视频播放| 国产熟女xx| 久久久久亚洲av毛片大全| 狠狠狠狠99中文字幕| 黄色丝袜av网址大全| 正在播放国产对白刺激| 一区二区三区激情视频| 国产三级在线视频| 精品午夜福利视频在线观看一区| 高清黄色对白视频在线免费看| www.精华液| 国产在线精品亚洲第一网站| 国产精品二区激情视频| 国产男靠女视频免费网站| or卡值多少钱| 国产成人av教育| 欧美另类亚洲清纯唯美| 国产成人免费无遮挡视频| 老司机靠b影院| 99国产精品免费福利视频| 亚洲男人天堂网一区| 制服丝袜大香蕉在线| 亚洲国产精品成人综合色| 一区二区三区激情视频| 亚洲欧美日韩无卡精品| 色在线成人网| 黄色视频不卡| 91av网站免费观看| 欧美老熟妇乱子伦牲交| 国产一区二区三区在线臀色熟女| 国产精品 国内视频| 国内久久婷婷六月综合欲色啪| 日日夜夜操网爽| 色婷婷久久久亚洲欧美| 91大片在线观看| 级片在线观看| 色哟哟哟哟哟哟| 亚洲成人国产一区在线观看| 十八禁人妻一区二区| 最近最新免费中文字幕在线| 午夜成年电影在线免费观看| 99国产精品99久久久久| 国产av一区在线观看免费| 狂野欧美激情性xxxx| 成人三级做爰电影| 欧美日韩福利视频一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 夜夜看夜夜爽夜夜摸| 曰老女人黄片| 欧美日韩亚洲国产一区二区在线观看| 色av中文字幕| 亚洲中文日韩欧美视频| 91成人精品电影| 久久久水蜜桃国产精品网| 这个男人来自地球电影免费观看| 国产片内射在线| 国产精品98久久久久久宅男小说| 婷婷六月久久综合丁香| 美女高潮喷水抽搐中文字幕| 大陆偷拍与自拍| 日韩欧美免费精品| 欧美丝袜亚洲另类 | 国产成人av教育| 精品欧美一区二区三区在线| 欧美日韩瑟瑟在线播放| 男人舔女人的私密视频| 又紧又爽又黄一区二区| 亚洲第一青青草原| 一进一出抽搐gif免费好疼| 老汉色∧v一级毛片| 真人做人爱边吃奶动态| 成人国产综合亚洲| 免费看a级黄色片| 男女床上黄色一级片免费看| 18禁美女被吸乳视频| 亚洲伊人色综图| 日韩欧美一区二区三区在线观看| 国产成人精品久久二区二区91| 亚洲性夜色夜夜综合| 巨乳人妻的诱惑在线观看| 国产亚洲精品av在线| 午夜老司机福利片| 夜夜夜夜夜久久久久| 国产精品亚洲美女久久久| 成人亚洲精品一区在线观看| 欧美日本亚洲视频在线播放| 国产真人三级小视频在线观看| 国产真人三级小视频在线观看| 成人特级黄色片久久久久久久| 12—13女人毛片做爰片一| 国产精品爽爽va在线观看网站 | www.www免费av| 欧美午夜高清在线| 亚洲精品一卡2卡三卡4卡5卡| 日本 av在线| 18禁黄网站禁片午夜丰满| 成人av一区二区三区在线看| 黄片大片在线免费观看| 久久精品91蜜桃| 嫁个100分男人电影在线观看| 岛国在线观看网站| 国产黄a三级三级三级人| 色av中文字幕| 999久久久精品免费观看国产| 久热这里只有精品99| 日韩视频一区二区在线观看| 久久久国产成人精品二区| 国产精品免费视频内射| 男女床上黄色一级片免费看| aaaaa片日本免费| 色综合婷婷激情| 日韩欧美三级三区| 最新在线观看一区二区三区| 不卡一级毛片| 日本a在线网址| 校园春色视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品,欧美在线| 国产麻豆成人av免费视频| 精品久久久精品久久久| 日韩欧美一区视频在线观看| 美女 人体艺术 gogo| 欧美日韩精品网址| 国产黄a三级三级三级人| 此物有八面人人有两片| АⅤ资源中文在线天堂| 老司机福利观看| 黑人巨大精品欧美一区二区蜜桃| 露出奶头的视频| 老司机午夜福利在线观看视频| 久久久久国产精品人妻aⅴ院| 免费高清视频大片| av片东京热男人的天堂| 午夜两性在线视频| 极品教师在线免费播放| 久久午夜综合久久蜜桃| 波多野结衣高清无吗| 欧美黄色片欧美黄色片| 亚洲色图 男人天堂 中文字幕| 女人爽到高潮嗷嗷叫在线视频| 国产一卡二卡三卡精品| 99在线视频只有这里精品首页| 国产乱人伦免费视频| 日本免费a在线| 精品第一国产精品| 国产亚洲精品第一综合不卡| 在线观看午夜福利视频| 波多野结衣巨乳人妻| 亚洲成人国产一区在线观看| 久久香蕉激情| 757午夜福利合集在线观看| 国产精品二区激情视频| 老司机午夜福利在线观看视频| 久久久久久人人人人人| 免费搜索国产男女视频| 怎么达到女性高潮| 神马国产精品三级电影在线观看 | 精品无人区乱码1区二区| 亚洲专区中文字幕在线| 制服诱惑二区| 亚洲熟女毛片儿| 精品免费久久久久久久清纯| 欧美日本视频| 亚洲精品国产区一区二| videosex国产| 久久这里只有精品19| 亚洲国产高清在线一区二区三 | 久久精品亚洲熟妇少妇任你| 一区二区三区高清视频在线| av网站免费在线观看视频| 国产片内射在线| 一本综合久久免费| 精品国产美女av久久久久小说| 欧美一级毛片孕妇| 午夜精品国产一区二区电影| 很黄的视频免费| 一区二区日韩欧美中文字幕| 午夜福利成人在线免费观看| 伊人久久大香线蕉亚洲五| 中亚洲国语对白在线视频| 欧美 亚洲 国产 日韩一| 琪琪午夜伦伦电影理论片6080| 成熟少妇高潮喷水视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久久久99蜜臀| 日本在线视频免费播放| 成人av一区二区三区在线看| 欧美在线一区亚洲| 久久久久国产精品人妻aⅴ院| 精品久久久久久久毛片微露脸| 看免费av毛片| 国产精品一区二区三区四区久久 | 99国产极品粉嫩在线观看| 午夜精品久久久久久毛片777| 精品不卡国产一区二区三区| 欧美在线一区亚洲| 大型黄色视频在线免费观看| 国产精品久久视频播放| 国产视频一区二区在线看| 母亲3免费完整高清在线观看| 国产精品乱码一区二三区的特点 | 欧美久久黑人一区二区| 在线观看一区二区三区| 亚洲色图 男人天堂 中文字幕| 亚洲国产高清在线一区二区三 | 欧美绝顶高潮抽搐喷水| 免费高清视频大片| 国产成人av激情在线播放| 亚洲熟妇中文字幕五十中出| 亚洲欧美精品综合久久99| 亚洲国产欧美日韩在线播放| 最近最新免费中文字幕在线| 国产精品一区二区精品视频观看| 黄频高清免费视频| 九色国产91popny在线| 91成人精品电影| 婷婷丁香在线五月| 男人舔女人下体高潮全视频| 淫妇啪啪啪对白视频| 一夜夜www| 身体一侧抽搐| 成熟少妇高潮喷水视频| 日韩欧美三级三区| 国产精品电影一区二区三区| 欧美激情极品国产一区二区三区| 91大片在线观看| 18禁国产床啪视频网站| 国产麻豆成人av免费视频| 老司机靠b影院| 亚洲五月婷婷丁香| 老司机深夜福利视频在线观看| 日本a在线网址| 男男h啪啪无遮挡| 1024香蕉在线观看| 久久伊人香网站| 久久亚洲精品不卡| 亚洲欧美精品综合久久99| avwww免费| 亚洲天堂国产精品一区在线| 国产精品,欧美在线| 国产精品一区二区免费欧美| 欧美老熟妇乱子伦牲交| 两个人看的免费小视频| 国产91精品成人一区二区三区| 久久热在线av| 99国产精品免费福利视频| 女性被躁到高潮视频| 国产成人av教育| 在线观看一区二区三区|