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

    基于ICEEMD-RAKSVD的地震弱信號去噪方法

    2022-05-09 02:41:32朱劍兵
    物探化探計(jì)算技術(shù) 2022年2期
    關(guān)鍵詞:正則字典噪音

    朱劍兵

    (中國石化 勝利油田分公司物探研究院,東營 257022)

    0 引言

    眾所周知,地震勘探由于噪音來源多,一些弱信號往往會(huì)淹沒在噪音里,常規(guī)方法很難去識別弱信號中的構(gòu)造、儲層信息。為了有效識別弱信號,前人引入了很多方法,其中包括:基于非線性動(dòng)力系統(tǒng)的混沌理論[1]、基于獨(dú)立分量分析的盲源分離技術(shù)[2]、基于奇異值分解(SVD)的方法[3]、基于經(jīng)驗(yàn)?zāi)J椒纸?EMD)[4]和基于高階統(tǒng)計(jì)量[5]等識別技術(shù)。

    Aharon[6]提出了KSVD算法;韓文功等[7]提出的SVD分解方法可以較好進(jìn)行弱信號的識別和去噪。KSVD算法由于需要進(jìn)行多次SVD分解的操作,因此,耗時(shí)過長, Rubinstein[8]在KSVD的基礎(chǔ)上提出了AKSVD算法(Approximate KSVD),但是其效果比較差;Dumitrescud等[9]在KSVD基礎(chǔ)上引入了正則化,得到RKSVD與RAKSVD(正則化近似KSVD)方法,并進(jìn)行了對比,RAKSVD在效果幾乎與KSVD一致的情況下,大大提高了效率。在模態(tài)分解去噪方面,Huang等[10]在1998年首次提出了經(jīng)驗(yàn)?zāi)B(tài)分解理論(EMD)。然而EMD也并非完美,仍存有一定的缺陷性,諸如部分模態(tài)之間的混疊現(xiàn)象、端點(diǎn)效應(yīng)以及篩分準(zhǔn)則等;Wu等[11]引入了總體經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD),對EMD算法進(jìn)行了的改進(jìn),但又出現(xiàn)了新的問題,恢復(fù)信號中仍含有一定程度的殘余噪音成分,無法將其消除干凈;Torres[12]對EEMD算法進(jìn)行了改進(jìn),引入的完備總體經(jīng)驗(yàn)?zāi)B(tài)分解(CEEMD)擁有高精度去噪的優(yōu)勢,但也不能盲目去除高頻模態(tài),其中還保存著很多有效信息;楊濤[13]結(jié)合KSVD和CEEMD的特點(diǎn),利用KSVD對弱信號的識別優(yōu)勢及CEEMD的高效去噪,設(shè)計(jì)了一個(gè)基于KSVD的CEEMD去噪識別方法,去噪效果較好,但計(jì)算效率相對較低。筆者將改進(jìn)的CEEMD(ICEEMD)與正則化近似KSVD(RAKSVD)有機(jī)結(jié)合起來,形成了基于ICEEMD-RAKVSD的去噪方法,該方法能夠?qū)Φ卣鹑跣盘栍行ёR別,取得了較好的弱信號去噪效果,同時(shí)計(jì)算效率也得到了提升。

    1 RAKSVD方法原理

    由于KSVD中每次原子更新都要進(jìn)行一次SVD分解,從而導(dǎo)致了計(jì)算效率較低。Rubinstein等[8]提出了一種近似KSVD算法(AKSVD),該算法能夠明顯提高計(jì)算效率,但由于引入了近似,SVD分解結(jié)果存在一定的誤差,因而影響了分解效果。針對這一問題,Irofti等[14]提出了正則化近似KSVD(RAKSVD)算法,建立嶺估計(jì)模型,解決稀疏分解中的病態(tài)問題,從而使SVD分解結(jié)果的誤差得到很大程度上的控制,提高了稀疏分解效果。

    1.1 KSVD字典

    KSVD字典是學(xué)習(xí)型字典之一,具有更適合反映數(shù)據(jù)結(jié)構(gòu)特點(diǎn)的優(yōu)勢,能獲得較好的去噪效果。KSVD算法通過對初始化字典進(jìn)行逐列更新,用來尋找一個(gè)最佳字典的過程,假定訓(xùn)練的信號為Y∈Rm×n,其中m為信號的行數(shù),n為信號的列數(shù),每列對應(yīng)于一個(gè)訓(xùn)練樣本。學(xué)習(xí)型字典模型可以描述為式(1)。

    (1)

    式中:D∈Rm×p為待學(xué)習(xí)字典,p為字典的列數(shù);X∈Rp×n為稀疏表示系數(shù);‖·‖0表示l0范數(shù);‖·‖2表示l2范數(shù);T為每列個(gè)數(shù)的最大值。

    首先固定字典D,利用正交匹配追蹤算法(OMP)對已知樣本數(shù)據(jù)進(jìn)行分解,獲取初始稀疏分解系數(shù)矩陣。

    字典更新實(shí)施步驟如下:

    1)記當(dāng)前更新的第k列字典原子為dk,X中使用到dk原子的索引集合為Ik。

    2)計(jì)算殘差:

    (2)

    式中:xj為x的第j列;Ek為去除dk外的殘差。

    (3)

    式中:U表示左奇異矩陣;V表示右奇異矩陣;Λ是奇異值合集;r表示非零奇異值的數(shù)量;ui表示左奇異矩陣的第i行;σi表示第i個(gè)奇異值;νi表示右奇異矩陣的第i列。

    4)原子和稀疏系數(shù)更新為:

    xk=σ1ν1,dk=u1

    (4)

    1.2 近似KSVD

    Rubinstein等[8]在KSVD的基礎(chǔ)上提出了AKSVD算法,主要是在上述步驟3)中采用近似計(jì)算更新原子,AKSVD算法更新后的原子和相關(guān)系數(shù)分別為式(5)和式(6)。

    (5)

    (6)

    1.3 正則化AKSVD

    稀疏分解是一種不適定問題,而正則化方法能有效提高分解效果。正則化字典習(xí)模型可描述為式(7)。

    (7)

    式中,μ>0為正則化參數(shù)。同樣的,記當(dāng)前更新的字典原子為dk,X中使用到原子的索引集合為Ik,則公式(2)可以變?yōu)槭?8)。

    μ‖Xk,Ik‖2

    (8)

    由于相關(guān)系數(shù)里除了第k行外,都為固定項(xiàng)。

    由公式(2)則可得到:

    (9)

    (10)

    根據(jù)式(10),從KSVD算法的正則化[9]可以得出,正則化約束后的更新原子方法和原算法一致,而正則化參數(shù)是應(yīng)用到了稀疏系數(shù)的更新,使系數(shù)變?yōu)?/(1+μ)倍,通過給稀疏系數(shù)一個(gè)權(quán)值,每次通過這個(gè)權(quán)值影響后續(xù)的原子更新,最終完成字典的更新。所以RAKSVD的原子和稀疏系數(shù)更新[14]為式(11)和式(12)。

    (11)

    (12)

    1.4 效果對比

    圖1為KSVD、AKSVD、RAKSVD三種訓(xùn)練字典的誤差對比圖,其中訓(xùn)練集Y的稀疏度為5,字典維數(shù)為64*256,正則化參數(shù)μ為0.1,整個(gè)分解過程重復(fù)迭代120次,取平均值得到均方根誤差RMSE(單位無量綱)隨迭代次數(shù)的關(guān)系。

    圖1 誤差對比圖Fig.1 Error comparison

    由圖1可以看出,這三種KSVD訓(xùn)練字典在迭代20次后RMSE基本趨于穩(wěn)定,其中AKSVD訓(xùn)練效果里相對最差,RAKSVD與KSVD訓(xùn)練效果接近,經(jīng)前人研究AKSVD相比于KSVD效率提高,故而RAKSVD在提高效率的同時(shí),也能達(dá)到較好的分解效果。

    2 ICEEMD方法原理

    首先,EEMD算法在解決EMD存在的模態(tài)混疊、端點(diǎn)效應(yīng)等缺陷的同時(shí),也導(dǎo)致了分解的IMF分量存在噪音殘余,而CEEMD通過在不同階段加入不同高斯白噪音進(jìn)行調(diào)整,比例大小根據(jù)待分解信號自適應(yīng)計(jì)算得到,同時(shí)每個(gè)符合條件的模態(tài)分量對應(yīng)的計(jì)算殘差值是固定唯一的。不僅解決了模態(tài)混疊問題,同時(shí)使原始信號復(fù)原的精度變得更高。

    按照Torres[12]所設(shè)計(jì)的CEEMD方法步驟對信號做處理,原始信號中各個(gè)頻率成分的信號,基本上都能夠得到較好的分解,但該算法仍存在可能出現(xiàn)虛假模態(tài)的問題,同時(shí)由于高斯白噪聲的加入,分解的模態(tài)中會(huì)混入一定程度的噪聲。為了解決上述問題,改進(jìn)的CEEMD(Improved CEEMD,ICEEMD)算法重新定義了各個(gè)模態(tài)局部均值的求取。

    假設(shè):Ej(·)為CEEEMD分解得到的第j階模態(tài),s(t)為原始地震信號,ωi為加入的第i個(gè)高斯白噪聲,白噪聲加入的比例為εk。則改進(jìn)后的ICEEMD實(shí)現(xiàn)步驟如下:

    1)將不同正負(fù)成對的白噪聲分別加入原始信號中,合成得到T個(gè)新信號(T為偶數(shù)),計(jì)算集合平均得到IMF1(t)。

    (13)

    2)令r0(t)=s(t),對k=1、…、K,計(jì)算k階殘差rk(t)為式(14)。

    rk(t)=rk-1(t)-IMFk(t)

    (14)

    3)對k階殘差添加不同比例系數(shù)的白噪聲:rk(t)+εkEk+1(ωi(t)),經(jīng)EMD分解得到IMF1(t),計(jì)算總體平均并將其作為IMF(k+1)(t)。

    (15)

    4)重復(fù)上述步驟,直到達(dá)到迭代條件,得到最終的殘差E(t)為式(16)。

    (16)

    CEEMD去噪往往是直接去除IMF1或IMF1和IMF2分量,但其中往往包含有效信息,所以,雖能有效去除高頻噪音,但對于弱信號去噪的問題,容易丟掉真正的有效信號,使同相軸不連續(xù)。

    相比于CEEMD算法,改進(jìn)后的CEEMD(ICEEMD)算法可以加入更少的白噪聲個(gè)數(shù),就能使信號得到很好的模態(tài)分離,同時(shí),在相同參數(shù)的條件下,ICEEMD分解的效果更好,計(jì)算效率更高。

    3 去噪方法

    KSVD字典稀疏去噪方法是屬于對壓縮感知原理的改進(jìn)算法,通過訓(xùn)練樣本求取冗余字典,較好地刻畫出信號特征,從而更加精確地重構(gòu)原信號,實(shí)現(xiàn)噪聲去除和保留有效信號的目標(biāo)。KSVD字典雖然能將噪音基本壓制和消除,但在信號的波峰波谷處仍會(huì)存有殘留一些毛刺。ICEEMD算法不需要提前設(shè)置基函數(shù),它是以信號本身的特點(diǎn)作為基礎(chǔ)來自動(dòng)分解和重構(gòu)信號,打破了測不準(zhǔn)原理的限制,達(dá)到了高精度識別弱信號的能力。

    具體實(shí)現(xiàn)步驟為:

    1)利用ICEEMD算法對地震信號進(jìn)行分解得到固有模態(tài)函數(shù)集,根據(jù)分解效果,將其劃分為噪聲主導(dǎo)模態(tài)、信號主導(dǎo)模態(tài)和過渡模態(tài)。

    2)對于噪聲主導(dǎo)模態(tài),則直接舍棄。

    3)對于信號主導(dǎo)模態(tài),將其與原始信號疊加,形成主導(dǎo)信號集,通過RAKSVD方法進(jìn)行主導(dǎo)信號集的分解和去噪。

    4)對于過渡模態(tài),先將其疊加,再進(jìn)行ICEEMD分解;舍棄噪音模態(tài)后,與原始信號疊加,重構(gòu)過渡信號集,通過RAKSVD方法進(jìn)行過渡信號集的分解和去噪。

    5)將步驟3)和步驟4)中的兩個(gè)分解結(jié)果重構(gòu)并平均,得到最終的去噪結(jié)果。

    4 模型測試

    含有弱信號模型如圖2(a)所示,模型共有25道,每道300個(gè)采樣點(diǎn),Ricker子波主頻25 Hz。在200 ms處有一弱信號。圖2(b)為含噪模型,加入噪聲后的峰值信噪比(PSNR)為1.2,圖2中強(qiáng)軸能夠明顯看到,但弱軸在噪聲的埋沒中已經(jīng)不明顯了。對圖2(b)分別采用五種算法進(jìn)行去噪處理,得到的結(jié)果如圖3所示,圖4為模型的15道的單道方法去噪結(jié)果,表1是各方法的去噪效果與時(shí)間對比。

    表1 五種方法針對模型數(shù)據(jù)的去噪效果對比Tab.1 Comparison of the effects of five methods for denoising the model data

    圖3 五種方法針對模型數(shù)據(jù)的去噪結(jié)果Fig.3 The result of five methods for denoising the model data(a)KSVD去噪;(b)RAKSVD去噪;(c)CEEMD去噪;(d)CEEMD-KSVD去噪;(e)ICEEMD-RAKSVD去噪

    圖4 單道去噪方法結(jié)果對比Fig.4 Comparison of single channel denoising results

    從圖3、圖4和表1可以看到,如果CEEMD單獨(dú)使用,則對弱信號的識別和去噪效果最差;KSVD和RAKSVD對模型中弱信號的識別和去噪效果較好,說明該類方法相比于CEEMD,對弱信號識別有著獨(dú)特的優(yōu)勢,但仍有一些毛刺的殘留,在實(shí)際應(yīng)用中,當(dāng)?shù)叵陆橘|(zhì)不均勻或由薄互層形成復(fù)波時(shí),其缺點(diǎn)會(huì)被放大;CEEMD-KSVD和ICEEMD-RAKSVD是將模態(tài)分解和字典訓(xùn)練聯(lián)合在一起進(jìn)行去噪,將模態(tài)分解的信號分解為主導(dǎo)模態(tài)和過渡模態(tài),利用KSVD和RAKSVD對弱信號識別的優(yōu)勢,最終使含噪弱信號中的噪聲得到了較好的壓制,且保留了有效弱信號,ICEEMD-RAKSVD方法能夠達(dá)到最好的去噪效果。由表1可以看到,由于CEEMD算法需要多次迭代,每次迭代需要加入大量成對的高斯白噪聲,因此,其計(jì)算效率低。KSVD和RAKSVD計(jì)算時(shí)間少,RAKSVD較KSVD效率提高了60%。由于CEEMD的引入,CEEMD-KSVD和ICEEMD-RAKSVD算法明顯增加了計(jì)算時(shí)間,ICEEMD-RAKSVD 相對于CEEMD-KSVD效率提高6%,對于簡單的稀疏模型,CEEMD的計(jì)算時(shí)間占主導(dǎo),因此,計(jì)算效率提升幅度不是很大。

    5 實(shí)際資料處理

    針對實(shí)際資料,對本方法的去噪效果進(jìn)行了進(jìn)一步測試。圖5(a)為原始局部疊前地震道集,圖5中隨機(jī)噪聲十分發(fā)育,信噪比非常低,其中弱同相軸不清晰、不連續(xù)。上述五種方法的去噪結(jié)果如圖5(b)~圖5(f)所示,圖6為單道去噪結(jié)果對比展示,圖7為去噪后殘差剖面效果對比,去噪效果與時(shí)間對比見表2。

    表2 五種方法針對實(shí)際資料的去噪效果對比Tab.2 Comparison of the effects of five methods for denoising the actual data

    圖5 五種方法針對實(shí)際資料的去噪結(jié)果Fig.5 The result of five methods for denoising the actual data(a)原始疊前地震道集;(b)KSVD去噪;(c)RAKSVD去噪;(d)CEEMD去噪;(e)CEEMD-KSVD去噪;(f)ICEEMD-RAKSVD去噪

    圖6 單道去噪方法結(jié)果對比Fig.6 Comparison of single channel denoising results

    從圖5和表2可以看出,對于弱地震信號,實(shí)際資料的處理結(jié)果與模型測試得到的結(jié)論相一致:CEEMD去噪效果較差,去噪后PSNR最低,將部分弱同相軸的有效信號充當(dāng)噪音去除,弱同相軸損失嚴(yán)重;RAKSVD與KSVD稀疏去噪較明顯,但仍存留了許多噪音。CEEMD-KSVD和ICEEMD-RAKSVD算法都能夠很好地壓制噪聲,紅方框標(biāo)記可以看出ICEEMD-RAKSVD有效弱信號得以保留。從圖5(a)可以更為清晰看出五種方法對比,而圖6可以得出 CEEMD-KSVD和ICEEMD-RAKSVD算法去噪效果幾乎一致,能有效去除噪音,圖7藍(lán)圓圈中可以看出,CEEMD-KSVD去除的殘余還仍然保留些有效信號,而本文方法去除的噪音沒有保留有效信號,去噪效果更好,同時(shí)由于實(shí)際地震數(shù)據(jù)非常復(fù)雜,訓(xùn)練字典的時(shí)間大大增加,從表2可以看出,RAKSVD較KSVD計(jì)算時(shí)間縮短了65%,下降幅度較為明顯,ICEEMD-RAKSVD相對于CEEMD-KSVD時(shí)間縮短了40%。對實(shí)際地震資料的去噪處理,更能突顯本文方法在提高計(jì)算效率方面的優(yōu)勢。

    圖7 去噪后殘差剖面效果對比Fig.7 Comparison of residual profile effects after denoising(a)CEEMD_KSVD;(b)ICEEMD_RAKSVD

    6 總結(jié)

    由于弱地震信息往往受噪音影響更大,通過去噪能更好地識別弱信號。CEEMD去噪方法直接舍棄高頻模態(tài),容易造成高頻中有效信息的丟失,而完整地保留了其他含有噪音模態(tài),KSVD去噪方法和RAKSVD去噪方法,對于弱信號的識別效果較好,但去噪后殘留毛刺,其中RAKSVD計(jì)算效率得到明顯提高,CEEMD-KSVD去噪方法能達(dá)到較好的去噪效果,但計(jì)算效率相對較低,這里指出的基于ICEEMD-RAKSVD去噪方法能在保證較好的去噪效果情況下,還能有效提升了計(jì)算效率,在實(shí)際地震處理中提升效率尤為重要,該方法具有一定的應(yīng)用價(jià)值。

    猜你喜歡
    正則字典噪音
    開心字典
    家教世界(2023年28期)2023-11-14 10:13:50
    開心字典
    家教世界(2023年25期)2023-10-09 02:11:56
    噪音,總是有噪音!
    無法逃避的噪音
    剩余有限Minimax可解群的4階正則自同構(gòu)
    類似于VNL環(huán)的環(huán)
    噪音的小把戲
    白噪音的三種用法
    Coco薇(2017年9期)2017-09-07 22:09:28
    我是小字典
    正版字典
    讀者(2016年14期)2016-06-29 17:25:50
    国产一卡二卡三卡精品 | 国产欧美亚洲国产| 亚洲国产精品成人久久小说| xxxhd国产人妻xxx| 最新在线观看一区二区三区 | 久久人人97超碰香蕉20202| 2021少妇久久久久久久久久久| 黑人欧美特级aaaaaa片| 亚洲精品在线美女| 久久 成人 亚洲| 国产亚洲欧美精品永久| 久久性视频一级片| 在线天堂最新版资源| 亚洲精品美女久久久久99蜜臀 | 亚洲精品视频女| 女人高潮潮喷娇喘18禁视频| 99热全是精品| 一级毛片我不卡| av在线播放精品| 国产精品一区二区在线不卡| 久久久精品国产亚洲av高清涩受| 久久久久久久久久久久大奶| 国产av码专区亚洲av| 欧美激情极品国产一区二区三区| 80岁老熟妇乱子伦牲交| 欧美97在线视频| 熟妇人妻不卡中文字幕| 麻豆精品久久久久久蜜桃| 国产日韩欧美视频二区| 国产老妇伦熟女老妇高清| 99久国产av精品国产电影| 欧美乱码精品一区二区三区| 国产1区2区3区精品| 日本av手机在线免费观看| 中国国产av一级| 国产精品一区二区在线不卡| 又粗又硬又长又爽又黄的视频| 人人妻人人澡人人看| 一级,二级,三级黄色视频| 国产精品国产三级专区第一集| 国产女主播在线喷水免费视频网站| 成人亚洲欧美一区二区av| 亚洲久久久国产精品| 赤兔流量卡办理| 91成人精品电影| 亚洲图色成人| 亚洲精品中文字幕在线视频| 国产在视频线精品| 免费在线观看视频国产中文字幕亚洲 | 久久免费观看电影| 欧美激情高清一区二区三区 | 美女福利国产在线| 日日摸夜夜添夜夜爱| 天美传媒精品一区二区| 黑人欧美特级aaaaaa片| 在线观看国产h片| 国产片内射在线| 美女扒开内裤让男人捅视频| 国产精品亚洲av一区麻豆 | 国产免费又黄又爽又色| 天天影视国产精品| 国产日韩欧美在线精品| 国产精品麻豆人妻色哟哟久久| 高清在线视频一区二区三区| 久久久久久人人人人人| 国产精品蜜桃在线观看| 免费黄网站久久成人精品| 国产精品蜜桃在线观看| 久久久久久久国产电影| 亚洲美女搞黄在线观看| 老司机靠b影院| 18禁观看日本| 亚洲美女搞黄在线观看| 涩涩av久久男人的天堂| 十八禁高潮呻吟视频| 久久久精品免费免费高清| 熟女av电影| 国产在线视频一区二区| 久久久久精品人妻al黑| 女人久久www免费人成看片| 亚洲国产日韩一区二区| 19禁男女啪啪无遮挡网站| 亚洲成人av在线免费| 亚洲成人av在线免费| 中文字幕另类日韩欧美亚洲嫩草| 天堂中文最新版在线下载| 国产女主播在线喷水免费视频网站| 久久免费观看电影| 久久精品亚洲av国产电影网| 观看av在线不卡| 日本vs欧美在线观看视频| 日韩不卡一区二区三区视频在线| 亚洲美女视频黄频| h视频一区二区三区| 九九爱精品视频在线观看| 亚洲少妇的诱惑av| 午夜免费观看性视频| 免费在线观看完整版高清| 亚洲精品美女久久av网站| 美女脱内裤让男人舔精品视频| 1024视频免费在线观看| 欧美日韩一级在线毛片| 色婷婷av一区二区三区视频| 在线观看三级黄色| 国产精品 欧美亚洲| 9色porny在线观看| 啦啦啦视频在线资源免费观看| 在现免费观看毛片| 日韩中文字幕欧美一区二区 | 国产精品一二三区在线看| 麻豆精品久久久久久蜜桃| 女人高潮潮喷娇喘18禁视频| avwww免费| 欧美黄色片欧美黄色片| 女人久久www免费人成看片| 香蕉丝袜av| 亚洲精品久久久久久婷婷小说| 国产 一区精品| 少妇人妻 视频| 免费在线观看黄色视频的| 免费不卡黄色视频| 国产成人91sexporn| 色综合欧美亚洲国产小说| 制服人妻中文乱码| 大片免费播放器 马上看| 女性被躁到高潮视频| 欧美成人午夜精品| 女的被弄到高潮叫床怎么办| 亚洲精品在线美女| 亚洲av在线观看美女高潮| 国产亚洲精品第一综合不卡| 午夜福利一区二区在线看| 久久精品国产亚洲av涩爱| 亚洲七黄色美女视频| 久久99一区二区三区| 色精品久久人妻99蜜桃| 久久99精品国语久久久| 丝袜人妻中文字幕| 国产麻豆69| 18禁裸乳无遮挡动漫免费视频| 日本爱情动作片www.在线观看| 2021少妇久久久久久久久久久| 99热全是精品| 国产亚洲av片在线观看秒播厂| 日本黄色日本黄色录像| 在线看a的网站| 两性夫妻黄色片| av网站在线播放免费| 久久97久久精品| 老司机亚洲免费影院| 女人被躁到高潮嗷嗷叫费观| 少妇精品久久久久久久| 深夜精品福利| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av中文av极速乱| 又大又黄又爽视频免费| 免费久久久久久久精品成人欧美视频| 青青草视频在线视频观看| 亚洲婷婷狠狠爱综合网| 国产亚洲av高清不卡| 日韩制服丝袜自拍偷拍| 久久久欧美国产精品| 色网站视频免费| 汤姆久久久久久久影院中文字幕| xxx大片免费视频| 精品一区二区三区av网在线观看 | 午夜久久久在线观看| 国产免费视频播放在线视频| 亚洲精品日本国产第一区| 亚洲国产成人一精品久久久| 国产精品香港三级国产av潘金莲 | 国产欧美亚洲国产| 欧美激情高清一区二区三区 | a 毛片基地| 悠悠久久av| 午夜福利视频在线观看免费| 精品一品国产午夜福利视频| 久久av网站| 波野结衣二区三区在线| 最近最新中文字幕大全免费视频 | 久久天躁狠狠躁夜夜2o2o | 欧美日韩成人在线一区二区| 精品一区二区免费观看| 久久久国产一区二区| 午夜免费鲁丝| 又黄又粗又硬又大视频| 韩国高清视频一区二区三区| 在线观看人妻少妇| 啦啦啦 在线观看视频| 午夜福利在线免费观看网站| 亚洲欧美一区二区三区黑人| av免费观看日本| svipshipincom国产片| av卡一久久| 亚洲精品第二区| 极品少妇高潮喷水抽搐| 热99久久久久精品小说推荐| 亚洲国产日韩一区二区| 你懂的网址亚洲精品在线观看| 婷婷色麻豆天堂久久| 国产精品秋霞免费鲁丝片| 日本猛色少妇xxxxx猛交久久| 成人漫画全彩无遮挡| 国产在视频线精品| 国产爽快片一区二区三区| 精品视频人人做人人爽| 亚洲精品美女久久av网站| 亚洲精品自拍成人| 99久久精品国产亚洲精品| 人人澡人人妻人| 亚洲,一卡二卡三卡| 欧美人与性动交α欧美软件| 亚洲三区欧美一区| 黄色 视频免费看| 国产xxxxx性猛交| 久久av网站| 亚洲欧美一区二区三区黑人| 99热国产这里只有精品6| 久久久久人妻精品一区果冻| 欧美日韩视频高清一区二区三区二| 99九九在线精品视频| 9热在线视频观看99| 久久 成人 亚洲| 777久久人妻少妇嫩草av网站| 亚洲精品一区蜜桃| 欧美黄色片欧美黄色片| 国产精品成人在线| 人人妻人人添人人爽欧美一区卜| 麻豆精品久久久久久蜜桃| 日韩电影二区| 久久久国产精品麻豆| 一区二区三区精品91| 日韩av免费高清视频| 校园人妻丝袜中文字幕| 男人爽女人下面视频在线观看| 人人妻,人人澡人人爽秒播 | 美女国产高潮福利片在线看| 婷婷成人精品国产| 男女午夜视频在线观看| 国产精品国产av在线观看| 精品国产露脸久久av麻豆| 精品一区二区免费观看| 亚洲色图综合在线观看| 亚洲欧美激情在线| 日韩精品免费视频一区二区三区| 日韩人妻精品一区2区三区| 国产激情久久老熟女| 日韩 亚洲 欧美在线| 丁香六月天网| 捣出白浆h1v1| 别揉我奶头~嗯~啊~动态视频 | 久久av网站| 1024视频免费在线观看| www.熟女人妻精品国产| 狂野欧美激情性bbbbbb| 久久狼人影院| 51午夜福利影视在线观看| 可以免费在线观看a视频的电影网站 | 看免费av毛片| 熟妇人妻不卡中文字幕| 欧美日本中文国产一区发布| 最黄视频免费看| 日本爱情动作片www.在线观看| 国产成人欧美| 国产精品.久久久| 哪个播放器可以免费观看大片| 久久热在线av| 97人妻天天添夜夜摸| 美女中出高潮动态图| 精品国产露脸久久av麻豆| 亚洲免费av在线视频| 久久精品久久精品一区二区三区| 麻豆乱淫一区二区| 国产精品二区激情视频| 亚洲精品日韩在线中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久成人av| 亚洲国产成人一精品久久久| 亚洲久久久国产精品| 亚洲一区二区三区欧美精品| av网站在线播放免费| 老司机深夜福利视频在线观看 | 99国产综合亚洲精品| 中文字幕色久视频| 2021少妇久久久久久久久久久| 一级爰片在线观看| 国产高清不卡午夜福利| 操出白浆在线播放| 热re99久久国产66热| 婷婷成人精品国产| 99re6热这里在线精品视频| 中文欧美无线码| 久久久久久久久免费视频了| 纯流量卡能插随身wifi吗| 久久毛片免费看一区二区三区| 国产精品久久久久久人妻精品电影 | 日韩av不卡免费在线播放| 最新在线观看一区二区三区 | 精品久久久久久电影网| 人人妻人人添人人爽欧美一区卜| 亚洲国产欧美日韩在线播放| 亚洲,欧美,日韩| 欧美日韩视频精品一区| 人人妻人人添人人爽欧美一区卜| 国产精品久久久久久精品古装| 国产探花极品一区二区| 一区在线观看完整版| 欧美久久黑人一区二区| 亚洲欧美成人精品一区二区| 国产 一区精品| 久久 成人 亚洲| 久久久久久免费高清国产稀缺| kizo精华| 免费高清在线观看视频在线观看| 国产一区二区 视频在线| 超碰97精品在线观看| www.av在线官网国产| 又大又黄又爽视频免费| av卡一久久| 涩涩av久久男人的天堂| 亚洲一码二码三码区别大吗| 亚洲欧洲日产国产| 亚洲欧美精品综合一区二区三区| 熟女av电影| 成人手机av| 国产成人免费无遮挡视频| 成年av动漫网址| 国产男女超爽视频在线观看| 最近最新中文字幕免费大全7| 丝袜美足系列| 搡老岳熟女国产| 婷婷色麻豆天堂久久| 久久综合国产亚洲精品| 亚洲国产欧美日韩在线播放| 欧美日本中文国产一区发布| 视频区图区小说| 美女国产高潮福利片在线看| 久久毛片免费看一区二区三区| 亚洲天堂av无毛| 一级毛片我不卡| 男女午夜视频在线观看| 黄频高清免费视频| 免费久久久久久久精品成人欧美视频| 亚洲精品成人av观看孕妇| 国产精品二区激情视频| 亚洲欧洲国产日韩| 亚洲国产欧美网| 男人添女人高潮全过程视频| 最近最新中文字幕大全免费视频 | 国产高清国产精品国产三级| 久久精品久久久久久久性| 日韩,欧美,国产一区二区三区| 成人黄色视频免费在线看| 久久韩国三级中文字幕| 极品少妇高潮喷水抽搐| 国产精品久久久久久人妻精品电影 | 亚洲国产精品一区二区三区在线| 国产精品 欧美亚洲| 又黄又粗又硬又大视频| 波多野结衣av一区二区av| 肉色欧美久久久久久久蜜桃| 亚洲综合色网址| 日韩一区二区三区影片| 久久鲁丝午夜福利片| 一本—道久久a久久精品蜜桃钙片| 欧美国产精品va在线观看不卡| a级片在线免费高清观看视频| 操出白浆在线播放| 女人爽到高潮嗷嗷叫在线视频| 欧美精品高潮呻吟av久久| 欧美在线一区亚洲| 国产精品一区二区在线不卡| 午夜精品国产一区二区电影| 免费观看人在逋| 人人妻人人添人人爽欧美一区卜| 国产爽快片一区二区三区| 亚洲成人av在线免费| 黑人猛操日本美女一级片| 婷婷成人精品国产| 欧美最新免费一区二区三区| 免费久久久久久久精品成人欧美视频| 久久久久国产精品人妻一区二区| 亚洲,欧美精品.| 亚洲欧洲国产日韩| 亚洲精品日本国产第一区| 亚洲在久久综合| 国产成人欧美在线观看 | 各种免费的搞黄视频| 国产乱来视频区| 免费观看av网站的网址| 国产视频首页在线观看| 亚洲av综合色区一区| 欧美黄色片欧美黄色片| 亚洲成人国产一区在线观看 | 成人国语在线视频| 天天操日日干夜夜撸| 亚洲成av片中文字幕在线观看| 91成人精品电影| 国产精品二区激情视频| 欧美日韩视频精品一区| 久久97久久精品| 欧美av亚洲av综合av国产av | 久久久久精品性色| 亚洲国产精品999| 一级毛片 在线播放| 久久久久网色| 99国产精品免费福利视频| 国产国语露脸激情在线看| 男人爽女人下面视频在线观看| 免费在线观看完整版高清| 色婷婷av一区二区三区视频| 亚洲精品日本国产第一区| 深夜精品福利| 99热网站在线观看| 亚洲欧美日韩另类电影网站| 啦啦啦 在线观看视频| 大码成人一级视频| 国产日韩欧美视频二区| 夫妻性生交免费视频一级片| 国产精品欧美亚洲77777| 欧美日韩精品网址| 1024香蕉在线观看| 亚洲av男天堂| 日韩中文字幕欧美一区二区 | 黑丝袜美女国产一区| 亚洲国产av新网站| 97人妻天天添夜夜摸| 啦啦啦 在线观看视频| 国产亚洲欧美精品永久| 国产精品一二三区在线看| 99国产综合亚洲精品| 国产一区二区 视频在线| 免费黄网站久久成人精品| 一级片免费观看大全| 成人漫画全彩无遮挡| 国产在线一区二区三区精| 日韩av不卡免费在线播放| 国产av一区二区精品久久| 久久久精品国产亚洲av高清涩受| 午夜久久久在线观看| 五月天丁香电影| 色94色欧美一区二区| 黄片播放在线免费| 另类亚洲欧美激情| 色婷婷av一区二区三区视频| 老汉色∧v一级毛片| 1024视频免费在线观看| 咕卡用的链子| 精品亚洲乱码少妇综合久久| 欧美日韩成人在线一区二区| 女性被躁到高潮视频| 久久天堂一区二区三区四区| 狂野欧美激情性xxxx| 国产精品蜜桃在线观看| 亚洲成av片中文字幕在线观看| 欧美中文综合在线视频| videos熟女内射| 日本一区二区免费在线视频| 妹子高潮喷水视频| 亚洲一区二区三区欧美精品| 老司机影院成人| 香蕉丝袜av| netflix在线观看网站| 亚洲中文av在线| 久久综合国产亚洲精品| 午夜福利一区二区在线看| 亚洲国产欧美在线一区| 丝袜喷水一区| 性少妇av在线| 99久久人妻综合| 久久精品亚洲av国产电影网| 久久久久国产精品人妻一区二区| 亚洲一码二码三码区别大吗| 中文字幕精品免费在线观看视频| 午夜福利在线免费观看网站| 日日爽夜夜爽网站| 在线天堂中文资源库| 久久影院123| 国产精品国产av在线观看| 夫妻午夜视频| 日韩 欧美 亚洲 中文字幕| 日韩大片免费观看网站| 老熟女久久久| 亚洲av成人精品一二三区| 熟妇人妻不卡中文字幕| 日本av免费视频播放| 亚洲精品国产区一区二| 久久性视频一级片| 在线免费观看不下载黄p国产| 电影成人av| 9热在线视频观看99| 男的添女的下面高潮视频| 久久精品亚洲熟妇少妇任你| 成人亚洲欧美一区二区av| 亚洲国产欧美网| 久久精品人人爽人人爽视色| 久久 成人 亚洲| 亚洲中文av在线| 久久国产亚洲av麻豆专区| 大香蕉久久网| 免费观看性生交大片5| 十八禁网站网址无遮挡| 久久人人爽av亚洲精品天堂| 80岁老熟妇乱子伦牲交| 亚洲欧洲日产国产| 日本欧美视频一区| 国产一区亚洲一区在线观看| 欧美精品高潮呻吟av久久| 可以免费在线观看a视频的电影网站 | 自线自在国产av| 99热国产这里只有精品6| 女的被弄到高潮叫床怎么办| 一区二区三区精品91| 韩国av在线不卡| 久久性视频一级片| 超碰97精品在线观看| 亚洲天堂av无毛| 飞空精品影院首页| 日本一区二区免费在线视频| 熟妇人妻不卡中文字幕| 久久这里只有精品19| 午夜福利视频在线观看免费| av福利片在线| 免费黄网站久久成人精品| 日韩制服骚丝袜av| 丁香六月天网| 中文字幕色久视频| 精品国产一区二区三区四区第35| 爱豆传媒免费全集在线观看| 免费黄频网站在线观看国产| 成年女人毛片免费观看观看9 | 在线观看国产h片| 久热这里只有精品99| 久久ye,这里只有精品| 亚洲成人国产一区在线观看 | 久久国产亚洲av麻豆专区| 国产精品.久久久| netflix在线观看网站| 国产熟女欧美一区二区| 一本大道久久a久久精品| 久久精品国产亚洲av高清一级| 国产老妇伦熟女老妇高清| 午夜福利网站1000一区二区三区| 男人爽女人下面视频在线观看| 国产成人精品久久二区二区91 | 免费看不卡的av| 国产人伦9x9x在线观看| av视频免费观看在线观看| 在线观看国产h片| 人人妻人人爽人人添夜夜欢视频| 久久ye,这里只有精品| 国产一卡二卡三卡精品 | 免费看av在线观看网站| 日本色播在线视频| 亚洲国产欧美日韩在线播放| 香蕉国产在线看| 亚洲av综合色区一区| 哪个播放器可以免费观看大片| 晚上一个人看的免费电影| 少妇被粗大猛烈的视频| 老熟女久久久| 亚洲精品乱久久久久久| 久久国产精品男人的天堂亚洲| 中国三级夫妇交换| 国产免费又黄又爽又色| 2018国产大陆天天弄谢| 国产精品国产av在线观看| 久久久久久久久免费视频了| 啦啦啦 在线观看视频| 午夜91福利影院| 中文字幕另类日韩欧美亚洲嫩草| 超碰成人久久| 午夜福利网站1000一区二区三区| 亚洲五月色婷婷综合| 街头女战士在线观看网站| 亚洲精品国产av蜜桃| 99热网站在线观看| 美女午夜性视频免费| 久久久精品94久久精品| kizo精华| 黄色毛片三级朝国网站| 性高湖久久久久久久久免费观看| 欧美乱码精品一区二区三区| 亚洲七黄色美女视频| a 毛片基地| 久久久久久久久久久久大奶| 亚洲五月色婷婷综合| 久久精品亚洲av国产电影网| 黄色 视频免费看| 国产精品秋霞免费鲁丝片| 国产野战对白在线观看| 亚洲综合色网址| 宅男免费午夜| 啦啦啦在线免费观看视频4| 99国产综合亚洲精品| 咕卡用的链子| 午夜福利网站1000一区二区三区| 一级黄片播放器| 哪个播放器可以免费观看大片| 亚洲精品一区蜜桃| av网站在线播放免费| 国产不卡av网站在线观看| 9色porny在线观看| 啦啦啦 在线观看视频| 国产精品久久久久成人av| 色精品久久人妻99蜜桃| 悠悠久久av| e午夜精品久久久久久久| 国产精品久久久av美女十八| 国产亚洲午夜精品一区二区久久| avwww免费| 国产一区二区 视频在线| 女性生殖器流出的白浆| 国产精品久久久久成人av| 1024视频免费在线观看|