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

    基于動態(tài)時間規(guī)整ICA算法地震隨機噪聲壓制

    2018-11-07 01:17:50逯宇佳曹俊興田仁飛呂雪松
    石油物探 2018年5期
    關鍵詞:拉平同相軸規(guī)整

    逯宇佳,曹俊興,田仁飛,呂雪松,何 沂

    (成都理工大學地球物理學院,四川成都610059)

    獨立分量分析(ICA)是20世紀末出現(xiàn)的一種新型信號處理方法,該方法可在信源未知的條件下,僅依賴混合信號將混疊在一起的不同信號分離開來[1]。其基本思想是假設原始信號為相互獨立的非高斯信號,利用信號的高階統(tǒng)計量信息,提取觀測信號獨立分量成分,最終得到相互獨立的源信號的近似估計[2-7]。1994年,COMON[8]首次提出了基于特征值分解的獨立分量分析方法。1998年,CICHOCKI等[9]利用獨立分量分析對含有加性噪聲的混合信號進行解混,取得較好的效果。在地球物理領域中,ICA主要應用于地震數(shù)據(jù)去噪處理。2003年,劉喜武等[10]假設地震記錄中有效信號和隨機干擾在統(tǒng)計上獨立且服從非高斯分布,對相鄰兩道地震記錄進行ICA處理實現(xiàn)地震資料的去噪;LU等[11]提出基于獨立分量分析的最大峰度自適應衰減算法實現(xiàn)了模型數(shù)據(jù)的多次波壓制;2006年,MIRKO[12]將獨立分量分析應用于P波和S波的波場分離;2011年,DONNO[13]將獨立分量分析與最小二乘法相結合提出了改進的多次波消除算法;呂文彪等[14-16]先后利用改進的ICA算法對疊后數(shù)據(jù)、地震屬性以及疊前數(shù)據(jù)進行了去噪處理;2012年,王維強等[17]將經(jīng)驗模態(tài)分解(EMD)技術與ICA技術結合應用于地震信號隨機噪聲壓制;2016年,孫成禹等[18]通過構造度量數(shù)據(jù)非高斯性的目標函數(shù),將地震數(shù)據(jù)轉換到ICA域,結合閾值法實現(xiàn)了復雜地震資料的去噪處理。

    傳統(tǒng)的ICA算法如基于負熵的快速ICA算法只適用于地震同相軸較平緩的疊后數(shù)據(jù),且算法不夠穩(wěn)定,容易出現(xiàn)解混失敗現(xiàn)象,計算效率較低。針對這個問題,本文運用兩步奇異值分解法對傳統(tǒng)ICA算法的白化過程進行改進,將語音識別中的動態(tài)時間規(guī)整(DTW)與ICA算法相結合,利用動態(tài)時間規(guī)整進行地震道間的模式匹配,度量時間序列的相似程度,將地震數(shù)據(jù)的同相軸一一匹配,再利用ICA算法進行去噪處理。

    1 方法原理

    1.1 ICA算法

    獨立分量分析處理需要對觀測信號矩陣進行去均值、預白化、ICA迭代處理。其中傳統(tǒng)的白化過程需要利用類似于主成分分析(PCA)方法,如零時滯的協(xié)方差矩陣的特征值分解來完成。利用兩步奇異值分解改進白化過程。

    當觀察信號個數(shù)(m)大于源的個數(shù)(n)時,ICA算法中白化過程可利用兩步奇異值分解法(圖像處理中又稱AMUSE算法,即多個未知信號提取算法)[19]完成。這一過程能夠實現(xiàn)源信號與噪聲信號的分離。主要實現(xiàn)步驟如下。

    1) 觀測信號X減去它的均值,得到零均值觀測信號,實現(xiàn)X的中心化。

    2) 估計觀測信號X的相關矩陣,即為協(xié)方差矩陣RX(0):

    (1)

    3) 利用奇異值分解(SVD)算法對RX(0)進行分解:

    (2)

    式中:VS為前n個特征值ΛS(從大到小排序)相對應的特征向量;VN為ΛN中后面m-n個不重要特征值ΛN=diag{λn+1,λn+2,…,λm}相對應的特征向量,λ為特征值。

    5) 進行穩(wěn)健的預白化變換:

    (3)

    (4)

    (5)

    (6)

    經(jīng)過以上兩步奇異值分解,即可完成解混過程,避免了傳統(tǒng)方法中復雜的迭代過程,在提高計算效率的同時也增加了算法的穩(wěn)定性。

    1.2 動態(tài)時間規(guī)整(DTW)

    設時間序列X=(x1,x2,…,xn),Y=(y1,y2,…,yn),則X,Y的DTW距離[20-22]DDTW(xi,yj)定義為[23]:

    (7)

    式中:Dbase(xi,yj)表示向量點xi和yj之間的基距離,可以根據(jù)情況選擇不同的距離度量。不失一般性,本文使用歐氏距離作為度量。根據(jù)公式(7),得到如圖1所示的一個鄰接矩陣,然后采用回溯法在鄰接矩陣上對DDTW(xi,yj)值進行遞歸回溯搜索DTW彎曲路徑。

    圖1 DTW彎曲路徑示意

    DTW距離允許序列點自我復制后再進行對齊匹配,能夠很好地支持時間軸彎曲,并且它可以對非等長時間序列進行度量,也支持時間軸伸縮。DTW距離實際上就是確定序列X與Y上每個點之間的對齊匹配關系(點對匹配)。

    彎曲路徑必須滿足以下3個基本條件[24]。

    1) 邊界條件:路徑起始于點(x1,y1)、終止于點(xn,yn),它表示兩個序列的起始點和結束點對應匹配。

    2) 連續(xù)性:路徑上的任意兩個相鄰點(xi1,yj1)和(xi2,yj2)滿足條件0≤|i1-i2|≤1,0≤|j1-j2|≤1。

    3) 單調性:若(xi1,yj1)和(xi2,yj2)為路徑上前后兩個點,則須滿足i2-i1≥0,j2-j1≥0。

    滿足上述條件的彎曲路徑有很多,每條彎曲路徑都代表一種點對匹配關系。設彎曲路徑為F=(f1,f2,…,fk,…,fK),fk=(i,j)k是彎曲路徑上第k個元素,它表示xi和yj建立的匹配關系,路徑長度滿足max(n,m)≤K≤n+m-1。

    點對匹配關系中,點對基距離之和的最小值即為DTW距離,對應的彎曲路徑為最佳路徑。DTW距離表示為:

    (8)

    求解最佳路徑需要構造一個m行n列的累積距離矩陣Mm×n,矩陣中的每個元素γi,j定義為:

    γi,j=Dbase(xi,yj)+min{γi-1,j,γi-1,j-1,γi,j-1}

    (9)

    γi,j為序列X[1∶i]與序列Y[1∶j]的DTW距離,因此,DDTW(X,Y)=γm,n,γm,n可以用動態(tài)時間規(guī)整法求解[25]。

    1.3 基于動態(tài)時間規(guī)整的ICA算法實現(xiàn)步驟

    基于動態(tài)時間規(guī)整的ICA算法去噪過程主要分為以下3個步驟:

    1) 選取第1道地震記錄作為模型道(模型道可以任選),利用動態(tài)時間規(guī)整算法對正演記錄的同相軸進行特征匹配(具體實現(xiàn)過程如圖1所示),在保證上覆水平層不變的情況下對下層傾斜地層同相軸進行校正拉平;

    2) 利用兩步奇異值分解法進行去噪處理;

    3) 利用動態(tài)時間規(guī)整提取的道間時差將拉平且去噪后的地震數(shù)據(jù)還原為原始形態(tài),完成整個去噪過程。

    2 應用實例

    2.1 模型測試

    為了驗證基于動態(tài)時間規(guī)整的ICA算法在地震隨機噪聲壓制中的可行性,建立如圖2a所示的速度模型,模型包含水平界面、傾斜界面、斷層以及彎曲界面。通過正演模擬得到自激自收地震記錄如圖2b所示。在生成的正演記錄中加入信噪比為0.5的高斯隨機噪聲見圖3a。利用ICA算法直接(未使用DTW)對加噪數(shù)據(jù)進行去噪處理,圖3b為去噪后的地震剖面,圖3c為去除的噪聲。由圖3b可以看出,ICA算法在同相軸水平的位置去噪效果較好,而當同相軸傾斜或彎曲時,能夠明顯看到一部分有效波的信息被當成噪聲去除了(圖3c),有效信號損失嚴重,斷層位置處尤為明顯。

    利用動態(tài)時間規(guī)整ICA算法對加噪數(shù)據(jù)進行去噪處理。圖4a為動態(tài)時間規(guī)整特征匹配后的拉平剖面,圖4b為兩步奇異值分解法對拉平剖面的去噪結果,圖5a為基于動態(tài)時間規(guī)整ICA算法去噪的結果剖面,與圖3b對比可以看出,經(jīng)過動態(tài)時間規(guī)整處理的去噪結果中隨機噪聲得到壓制,且在復雜界面如斷層、傾斜和彎曲界面的有效波信號保留較好,差剖面中沒有殘留有效信號(圖5b)。

    2.2 疊前去噪應用

    高信噪比的地震數(shù)據(jù)是進行地震解釋的重要基礎。解釋工作開始前,往往需要對疊前CRP道集數(shù)據(jù)進行層拉平或去噪處理,但這一工作非常費時,因此在保證去噪效果的同時,算法的計算效率顯得尤為重要。為了驗證兩步奇異值分解法的實際效果,本文選用疊前CRP道集進行驗證。疊前CRP道集是來自地下同一反射點的一系列地震響應特征,其每一道數(shù)據(jù)之間有極大的相似性,因此可利用ICA算法對其進行去噪處理。

    圖2 速度模型(a)和正演記錄(b)

    圖3 利用ICA算法處理的去噪效果(未使用DTW)a 加噪剖面;b 去噪后剖面;c 去除的噪聲

    去噪效果如圖6所示。其中,圖6a為原始CRP道集,圖6b為基于負熵的快速ICA算法去噪結果,圖6c為兩步奇異值分解法的去噪剖面,圖6d為兩步奇異值分解法去除的噪聲。相比于原始記錄,兩種去噪方法都能夠有效去除原始剖面中的隨機噪聲,增強同相軸連續(xù)性。但是基于負熵的快速ICA算法處理結果中出現(xiàn)了由于迭代不收斂產(chǎn)生的壞道(圖6b紅色箭頭處),破壞了原有的有效波信息。而兩步奇異值分解法速度快,穩(wěn)定性高,不需要復雜的迭代過程,可以直接對矩陣進行分解,最終得到的剖面同相軸連續(xù),去噪效果有明顯的改善。圖7為圖6 中疊前CRP數(shù)據(jù)對應的疊加剖面。對比圖7a,圖7b和圖7c可以看出,利用兩步奇異值分解法對疊前資料進行去噪處理后所得疊后剖面同相軸連續(xù)性增強(紅色矩形框內(nèi)較為明顯),信噪比有所提高。

    圖4 基于動態(tài)時間規(guī)整ICA算法去噪效果(Ⅰ)a 動態(tài)時間規(guī)整特征匹配拉平剖面;b 兩步奇異值分解法對拉平剖面去噪結果

    圖5 基于動態(tài)時間規(guī)整ICA算法的去噪效果(Ⅱ)a 基于動態(tài)時間規(guī)整還原后的去噪結果剖面;b 去除的噪聲

    2.3 疊后應用實例

    將基于動態(tài)時間規(guī)整ICA算法應用于疊后地震資料隨機噪聲壓制,結果如圖8和圖9所示。其中,圖8a 是處理前的原始地震剖面;圖8b是經(jīng)過DTW處理后的層拉平剖面,由于原始地震資料的品質限制,中心位置拉平效果一般;圖8c是拉平后兩步奇異值算法去噪結果。圖9b是基于負熵的快速ICA算法去噪結果,圖9c是基于動態(tài)時間規(guī)整的ICA算法去噪結果。對比圖9a,圖9b和圖9c可以看出,針對地層傾角較大且有斷層存在的復雜剖面,基于負熵的快速ICA算法雖然能去除一部分的隨機噪聲,但是處理后的有效波同相軸能量變?nèi)?連續(xù)性反而變差,這是因為在處理過中相鄰地震道存在的時差使得地震道之間的相似性降低,從而導致算法不收斂產(chǎn)生壞道。采用基于動態(tài)時間規(guī)整的ICA算法去噪后,同相軸連續(xù)性明顯增強,弱振幅信號得到放大,整體能量更為均一,斷層信息更為突出,斷點清晰可見(圖9c紅色橢圓標注)。圖10對比了基于動態(tài)時間規(guī)整ICA算法去噪前、后的頻譜,可以看出,本文方法去噪前、后頻譜基本一致,未改變原始資料的頻譜特征。相比于傳統(tǒng)的方法,基于動態(tài)時間規(guī)整的ICA算法更適合復雜地震資料的隨機噪聲壓制,有較好的實用價值。

    圖6 疊前CRP道集去噪效果a 原始記錄;b 基于負熵的快速ICA算法去噪結果;c 兩步奇異值分解法去噪結果;d 兩步奇異值分解法去除的噪聲

    圖9 疊后剖面去噪效果a 原始剖面;b 基于負熵的快速ICA算法去噪結果;c 基于動態(tài)時間規(guī)整ICA算法去噪結果;d 基于動態(tài)時間規(guī)整ICA算法去除的噪聲

    圖10 基于動態(tài)時間規(guī)整ICA算法去噪前、后的頻譜對比

    3 結論

    本文將ICA算法與動態(tài)時間規(guī)整相結合,提出了基于動態(tài)時間規(guī)整ICA算法,實現(xiàn)了疊前疊后數(shù)據(jù)的隨機噪聲壓制處理,獲得以下結論和認識:

    1) 利用兩步奇異值分解法對矩陣進行穩(wěn)健白化,可將有效波和噪聲分離。提高了ICA算法的穩(wěn)定性和計算效率,隨機噪聲壓制效果明顯;

    2) 將動態(tài)時間規(guī)整與ICA算法相結合,解決了僅利用ICA算法對非平緩同相軸的地震資料去噪不適應性問題;

    3) 基于動態(tài)時間規(guī)整的ICA算法在實際資料應用中取得較好效果。既可以應用于疊前資料的隨機噪聲壓制,也可以應用于疊后地震資料,具有普適性。

    猜你喜歡
    拉平同相軸規(guī)整
    賣狗人
    虛同相軸方法及其在陸上地震層間多次波壓制中的應用
    300kt/a硫酸系統(tǒng)規(guī)整填料使用情況簡介
    一種改進的相關法自動拾取同相軸
    提高日用玻璃陶瓷規(guī)整度和表面光滑度的處理方法
    佛山陶瓷(2016年11期)2016-12-23 08:50:27
    電梯的建筑化藝術探索
    大觀(2016年9期)2016-11-16 10:31:30
    可重復使用飛行器進場著陸拉平縱向控制
    一種反射同相軸自動拾取算法
    基于發(fā)音機制的貪婪自適應語音時長規(guī)整算法
    計算機工程(2015年8期)2015-07-03 12:20:34
    基于同相軸追蹤的多次波衰減*
    亚洲图色成人| 国产老妇女一区| 亚洲精品久久久久久婷婷小说| 2022亚洲国产成人精品| 最近手机中文字幕大全| 亚洲精品亚洲一区二区| 啦啦啦韩国在线观看视频| 高清视频免费观看一区二区 | 久久久精品免费免费高清| 国产爱豆传媒在线观看| 国产高清有码在线观看视频| 久久99蜜桃精品久久| 欧美+日韩+精品| 亚洲人成网站在线播| 日本黄色片子视频| 国产高潮美女av| 国产伦在线观看视频一区| 国产色婷婷99| 国产欧美日韩精品一区二区| 国产精品一二三区在线看| 又粗又硬又长又爽又黄的视频| 纵有疾风起免费观看全集完整版 | 69av精品久久久久久| 精品久久久久久久久av| 国产 亚洲一区二区三区 | 国产极品天堂在线| 国产黄片美女视频| 成人漫画全彩无遮挡| 一个人看的www免费观看视频| 久久精品国产亚洲网站| 国产伦理片在线播放av一区| 亚洲精品乱久久久久久| 国产熟女欧美一区二区| 亚洲精品成人av观看孕妇| 九九久久精品国产亚洲av麻豆| av女优亚洲男人天堂| 3wmmmm亚洲av在线观看| 久久久久久久久中文| 少妇被粗大猛烈的视频| 国产片特级美女逼逼视频| 你懂的网址亚洲精品在线观看| 亚洲天堂国产精品一区在线| 国产精品美女特级片免费视频播放器| 一夜夜www| 精品久久久噜噜| 床上黄色一级片| 中文字幕人妻熟人妻熟丝袜美| 国产激情偷乱视频一区二区| 中文精品一卡2卡3卡4更新| 亚洲精品一二三| 亚洲欧美一区二区三区国产| 日日摸夜夜添夜夜爱| 日产精品乱码卡一卡2卡三| 欧美性感艳星| 成人综合一区亚洲| 亚洲国产色片| 午夜日本视频在线| 国产精品蜜桃在线观看| 一区二区三区高清视频在线| 十八禁国产超污无遮挡网站| 大香蕉97超碰在线| 国产黄a三级三级三级人| 亚洲性久久影院| 日韩 亚洲 欧美在线| 丰满人妻一区二区三区视频av| 国产综合懂色| 嫩草影院新地址| 欧美日韩在线观看h| 99久久九九国产精品国产免费| 中文在线观看免费www的网站| 青青草视频在线视频观看| 国产精品一区二区三区四区免费观看| 国产片特级美女逼逼视频| 国产免费又黄又爽又色| 国国产精品蜜臀av免费| 国产极品天堂在线| 欧美最新免费一区二区三区| 亚洲国产成人一精品久久久| 亚洲自偷自拍三级| 中文字幕免费在线视频6| 国产爱豆传媒在线观看| 男女那种视频在线观看| 嘟嘟电影网在线观看| 久久精品国产亚洲网站| 免费看av在线观看网站| 日韩欧美精品免费久久| 国产欧美另类精品又又久久亚洲欧美| 只有这里有精品99| 欧美日韩一区二区视频在线观看视频在线 | 亚洲人成网站高清观看| 国产一级毛片在线| 51国产日韩欧美| 啦啦啦中文免费视频观看日本| 日韩av免费高清视频| 七月丁香在线播放| 日韩人妻高清精品专区| 免费看不卡的av| 久久精品国产鲁丝片午夜精品| 国产乱人偷精品视频| 亚洲四区av| 国产av国产精品国产| 国产精品一区二区三区四区久久| 亚洲成人精品中文字幕电影| 纵有疾风起免费观看全集完整版 | 一二三四中文在线观看免费高清| 久久久国产一区二区| 少妇被粗大猛烈的视频| 亚洲乱码一区二区免费版| 一个人看的www免费观看视频| 国产淫语在线视频| 久久久久久久亚洲中文字幕| a级一级毛片免费在线观看| 欧美97在线视频| 国产 一区 欧美 日韩| 亚洲天堂国产精品一区在线| 久久精品国产鲁丝片午夜精品| 丝袜喷水一区| 伊人久久国产一区二区| 国产高清三级在线| 尤物成人国产欧美一区二区三区| 国产精品久久久久久久久免| 三级国产精品片| 久久久色成人| 国产成人a∨麻豆精品| 大片免费播放器 马上看| 最后的刺客免费高清国语| 午夜精品在线福利| 日本一二三区视频观看| 激情 狠狠 欧美| 亚洲欧洲日产国产| 一个人免费在线观看电影| 九色成人免费人妻av| 搡老乐熟女国产| 日日摸夜夜添夜夜爱| 有码 亚洲区| 岛国毛片在线播放| 欧美xxⅹ黑人| 国内精品美女久久久久久| 免费观看在线日韩| 亚洲成人久久爱视频| 日本-黄色视频高清免费观看| 三级经典国产精品| 99视频精品全部免费 在线| 亚洲经典国产精华液单| 成人亚洲精品一区在线观看 | 欧美日本视频| 国产成人精品福利久久| 伦理电影大哥的女人| 欧美激情国产日韩精品一区| 国产精品一区二区在线观看99 | 男女边摸边吃奶| 久久久久久久久大av| 一个人免费在线观看电影| 非洲黑人性xxxx精品又粗又长| 亚洲在线观看片| 国产一区二区三区综合在线观看 | 免费看美女性在线毛片视频| 亚洲一区高清亚洲精品| 久久亚洲国产成人精品v| 亚洲色图av天堂| 国产精品久久久久久久电影| 日韩欧美一区视频在线观看 | 亚洲国产高清在线一区二区三| av.在线天堂| 免费看光身美女| av黄色大香蕉| 午夜激情福利司机影院| 春色校园在线视频观看| 青春草亚洲视频在线观看| 在线免费观看的www视频| 亚洲精品自拍成人| 免费观看av网站的网址| 男的添女的下面高潮视频| 天堂中文最新版在线下载 | 国产精品一区二区三区四区免费观看| 中文字幕人妻熟人妻熟丝袜美| 蜜桃亚洲精品一区二区三区| 熟妇人妻不卡中文字幕| 国产v大片淫在线免费观看| 青青草视频在线视频观看| 国产精品一及| 成人毛片60女人毛片免费| 99re6热这里在线精品视频| 成人美女网站在线观看视频| 日韩av在线大香蕉| 成人国产麻豆网| 成人亚洲欧美一区二区av| 国产毛片a区久久久久| 国产毛片a区久久久久| 国产免费福利视频在线观看| 在线观看av片永久免费下载| 亚洲欧美日韩无卡精品| 男女下面进入的视频免费午夜| av又黄又爽大尺度在线免费看| 偷拍熟女少妇极品色| 午夜激情福利司机影院| 国产成人午夜福利电影在线观看| 22中文网久久字幕| 亚洲一级一片aⅴ在线观看| 欧美日韩国产mv在线观看视频 | 亚洲国产av新网站| 男女国产视频网站| 亚洲一区高清亚洲精品| 高清视频免费观看一区二区 | 男人和女人高潮做爰伦理| 又爽又黄a免费视频| 国产黄片美女视频| 免费av毛片视频| 色哟哟·www| 久久久久精品性色| 精品久久久久久久人妻蜜臀av| 亚洲激情五月婷婷啪啪| 亚洲18禁久久av| 午夜福利视频1000在线观看| 亚洲欧美清纯卡通| 亚洲精品一区蜜桃| 大话2 男鬼变身卡| 国产成人精品婷婷| 日本欧美国产在线视频| 真实男女啪啪啪动态图| 女人久久www免费人成看片| 亚洲av成人av| 精品99又大又爽又粗少妇毛片| 不卡视频在线观看欧美| av国产免费在线观看| 成人特级av手机在线观看| 成年人午夜在线观看视频 | 十八禁国产超污无遮挡网站| a级一级毛片免费在线观看| 亚洲欧美一区二区三区黑人 | 国产片特级美女逼逼视频| 看免费成人av毛片| 激情 狠狠 欧美| 在线观看人妻少妇| 日韩伦理黄色片| 成人av在线播放网站| 看十八女毛片水多多多| 午夜激情久久久久久久| 久久久色成人| 国产大屁股一区二区在线视频| 性插视频无遮挡在线免费观看| 人妻系列 视频| 18+在线观看网站| 哪个播放器可以免费观看大片| 亚洲精品久久午夜乱码| 国产人妻一区二区三区在| 午夜视频国产福利| av国产免费在线观看| 国产av国产精品国产| 国产精品蜜桃在线观看| 三级经典国产精品| 国产精品人妻久久久影院| 亚洲自偷自拍三级| 22中文网久久字幕| 男女国产视频网站| 久久人人爽人人片av| 波多野结衣巨乳人妻| 欧美日韩在线观看h| av在线播放精品| 色尼玛亚洲综合影院| 成年女人看的毛片在线观看| 成人特级av手机在线观看| 2021少妇久久久久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲av日韩在线播放| 七月丁香在线播放| 久久久久久久国产电影| 美女脱内裤让男人舔精品视频| 男女国产视频网站| 久久久久久久大尺度免费视频| 伊人久久精品亚洲午夜| 一级毛片 在线播放| 春色校园在线视频观看| 亚洲在线自拍视频| 亚洲精品aⅴ在线观看| av一本久久久久| 免费看日本二区| 人人妻人人澡人人爽人人夜夜 | 91久久精品国产一区二区成人| 欧美高清性xxxxhd video| 亚洲色图av天堂| 国产一区亚洲一区在线观看| 成人性生交大片免费视频hd| 午夜免费男女啪啪视频观看| 国产精品综合久久久久久久免费| 在线观看美女被高潮喷水网站| 只有这里有精品99| 黄片wwwwww| 欧美日韩视频高清一区二区三区二| 男人和女人高潮做爰伦理| 中文在线观看免费www的网站| 大陆偷拍与自拍| 亚洲av成人av| 美女高潮的动态| 精品酒店卫生间| 免费不卡的大黄色大毛片视频在线观看 | 韩国高清视频一区二区三区| 久久久亚洲精品成人影院| 日韩欧美精品免费久久| 网址你懂的国产日韩在线| 最近视频中文字幕2019在线8| 免费少妇av软件| 热99在线观看视频| 免费看av在线观看网站| 街头女战士在线观看网站| 亚洲欧美中文字幕日韩二区| 国产免费福利视频在线观看| 在现免费观看毛片| 三级男女做爰猛烈吃奶摸视频| 综合色丁香网| 国产午夜精品论理片| 九色成人免费人妻av| 最近最新中文字幕大全电影3| 婷婷色综合大香蕉| 99热这里只有是精品在线观看| 国内少妇人妻偷人精品xxx网站| 国产亚洲av嫩草精品影院| 国产精品一区二区在线观看99 | 免费av观看视频| 色哟哟·www| 成人综合一区亚洲| 午夜激情久久久久久久| 国产探花在线观看一区二区| 国产精品人妻久久久影院| av黄色大香蕉| 99久久中文字幕三级久久日本| 午夜精品国产一区二区电影 | 亚洲精品色激情综合| 伊人久久国产一区二区| 99热网站在线观看| 精品亚洲乱码少妇综合久久| 亚洲综合色惰| 尤物成人国产欧美一区二区三区| 亚洲av.av天堂| 狠狠精品人妻久久久久久综合| 2021天堂中文幕一二区在线观| 亚洲欧美一区二区三区国产| 亚洲av中文av极速乱| 在线 av 中文字幕| 三级国产精品欧美在线观看| 亚洲国产色片| 男女边摸边吃奶| 高清欧美精品videossex| 亚洲av二区三区四区| 日日干狠狠操夜夜爽| 熟妇人妻久久中文字幕3abv| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av涩爱| 舔av片在线| 人人妻人人看人人澡| 国国产精品蜜臀av免费| 亚洲精华国产精华液的使用体验| 亚洲一级一片aⅴ在线观看| 国产精品熟女久久久久浪| 美女被艹到高潮喷水动态| 久久精品久久久久久噜噜老黄| 国产人妻一区二区三区在| 一级二级三级毛片免费看| 小蜜桃在线观看免费完整版高清| 丰满乱子伦码专区| 国产欧美日韩精品一区二区| 亚洲最大成人av| 大又大粗又爽又黄少妇毛片口| 国产精品国产三级国产av玫瑰| 极品少妇高潮喷水抽搐| 亚洲精品乱久久久久久| 中文字幕久久专区| 色综合站精品国产| 日韩在线高清观看一区二区三区| 少妇裸体淫交视频免费看高清| 高清视频免费观看一区二区 | 熟女电影av网| 日本欧美国产在线视频| 午夜福利高清视频| 97精品久久久久久久久久精品| 亚洲精品亚洲一区二区| 美女高潮的动态| 亚洲国产精品sss在线观看| 少妇熟女欧美另类| 深夜a级毛片| 极品教师在线视频| 日韩精品有码人妻一区| 亚洲天堂国产精品一区在线| xxx大片免费视频| 大又大粗又爽又黄少妇毛片口| 2018国产大陆天天弄谢| 久久久欧美国产精品| 嫩草影院精品99| 日韩在线高清观看一区二区三区| 国产黄色小视频在线观看| 国产一区二区三区综合在线观看 | 国产av不卡久久| 免费观看精品视频网站| 国产伦精品一区二区三区视频9| 久久精品久久久久久久性| 色网站视频免费| 国产在视频线精品| 国产高潮美女av| 51国产日韩欧美| 汤姆久久久久久久影院中文字幕 | 两个人的视频大全免费| 亚洲成人久久爱视频| 亚洲av中文av极速乱| 久久99热6这里只有精品| 国产 一区精品| 午夜爱爱视频在线播放| 亚洲成人av在线免费| 夜夜看夜夜爽夜夜摸| 国产av在哪里看| 欧美日韩国产mv在线观看视频 | 97在线视频观看| 欧美成人午夜免费资源| 亚洲精品亚洲一区二区| 毛片女人毛片| 看十八女毛片水多多多| 嫩草影院精品99| 成人欧美大片| 嫩草影院精品99| 熟妇人妻久久中文字幕3abv| 六月丁香七月| 人妻夜夜爽99麻豆av| 欧美xxxx性猛交bbbb| 婷婷色麻豆天堂久久| 国产视频内射| 亚洲自拍偷在线| 99热6这里只有精品| 国产亚洲av片在线观看秒播厂 | 女人久久www免费人成看片| 偷拍熟女少妇极品色| 国产伦精品一区二区三区四那| 91av网一区二区| 六月丁香七月| 日本黄大片高清| av女优亚洲男人天堂| 男女啪啪激烈高潮av片| a级毛色黄片| 精品人妻视频免费看| 亚洲国产成人一精品久久久| 免费电影在线观看免费观看| 成年女人在线观看亚洲视频 | 免费黄色在线免费观看| 中文欧美无线码| 精品久久久久久成人av| 精品午夜福利在线看| 午夜视频国产福利| 91在线精品国自产拍蜜月| 白带黄色成豆腐渣| 亚洲欧美清纯卡通| 亚洲av中文av极速乱| 国产精品久久久久久久久免| 一级毛片 在线播放| 久久精品熟女亚洲av麻豆精品 | 成人美女网站在线观看视频| 搡老乐熟女国产| 久久鲁丝午夜福利片| 嫩草影院入口| 欧美精品一区二区大全| 婷婷色av中文字幕| 日韩制服骚丝袜av| 成人综合一区亚洲| 亚洲av一区综合| 人人妻人人澡人人爽人人夜夜 | 老司机影院毛片| 中文乱码字字幕精品一区二区三区 | 91精品伊人久久大香线蕉| 最近中文字幕2019免费版| 午夜精品在线福利| 日本与韩国留学比较| 精品99又大又爽又粗少妇毛片| 美女黄网站色视频| 小蜜桃在线观看免费完整版高清| 亚洲欧美日韩卡通动漫| 寂寞人妻少妇视频99o| 51国产日韩欧美| 亚洲成人中文字幕在线播放| 真实男女啪啪啪动态图| 我要看日韩黄色一级片| 亚洲经典国产精华液单| 亚洲欧美中文字幕日韩二区| 日本色播在线视频| 99热这里只有是精品50| 欧美一区二区亚洲| 波野结衣二区三区在线| 最近2019中文字幕mv第一页| 哪个播放器可以免费观看大片| 中文在线观看免费www的网站| 国产一区二区亚洲精品在线观看| 亚洲性久久影院| 97热精品久久久久久| 22中文网久久字幕| 成人亚洲精品av一区二区| 美女黄网站色视频| 欧美日本视频| 91精品伊人久久大香线蕉| 麻豆久久精品国产亚洲av| 99热全是精品| 国产精品久久久久久精品电影| a级一级毛片免费在线观看| av在线播放精品| 精品熟女少妇av免费看| 日韩一区二区三区影片| 草草在线视频免费看| 成人欧美大片| 91久久精品国产一区二区成人| 身体一侧抽搐| 日本黄色片子视频| 国产精品国产三级国产av玫瑰| 国产精品人妻久久久久久| 成人午夜高清在线视频| 亚洲欧美精品自产自拍| 国产精品伦人一区二区| 亚洲国产欧美在线一区| 国产成人freesex在线| 能在线免费看毛片的网站| 少妇猛男粗大的猛烈进出视频 | 一级片'在线观看视频| 欧美人与善性xxx| 午夜免费激情av| 久久久久久九九精品二区国产| 欧美xxxx性猛交bbbb| 汤姆久久久久久久影院中文字幕 | 亚洲,欧美,日韩| 欧美日韩视频高清一区二区三区二| a级一级毛片免费在线观看| 精品午夜福利在线看| 久久久久久久亚洲中文字幕| 爱豆传媒免费全集在线观看| 精品久久久噜噜| 九色成人免费人妻av| 黄色欧美视频在线观看| 高清日韩中文字幕在线| 又黄又爽又刺激的免费视频.| 我的女老师完整版在线观看| 精品久久久久久久久久久久久| 九九久久精品国产亚洲av麻豆| 亚洲久久久久久中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 午夜久久久久精精品| 亚洲精品色激情综合| 国产探花在线观看一区二区| 夫妻午夜视频| 国产成人a区在线观看| av又黄又爽大尺度在线免费看| 国产在视频线在精品| 国产一区亚洲一区在线观看| 免费观看精品视频网站| 亚洲欧洲日产国产| 免费看不卡的av| 超碰av人人做人人爽久久| 日韩国内少妇激情av| 一二三四中文在线观看免费高清| 岛国毛片在线播放| 春色校园在线视频观看| 精品人妻视频免费看| 亚洲精品国产av蜜桃| 国产精品国产三级国产av玫瑰| 国产精品99久久久久久久久| 99久久精品国产国产毛片| 久久国产乱子免费精品| 国产一区有黄有色的免费视频 | 国产成人精品一,二区| 国产黄片视频在线免费观看| 国产精品久久久久久久电影| 又爽又黄无遮挡网站| 久久这里只有精品中国| a级一级毛片免费在线观看| 精品久久久久久久末码| 亚洲精品,欧美精品| 亚洲精品一区蜜桃| 国产高清国产精品国产三级 | 又大又黄又爽视频免费| 久久久久网色| 两个人的视频大全免费| 一级毛片我不卡| 国产老妇女一区| 免费观看的影片在线观看| 寂寞人妻少妇视频99o| 亚洲人与动物交配视频| 人人妻人人澡欧美一区二区| 久久99热6这里只有精品| 日本色播在线视频| 直男gayav资源| 亚洲怡红院男人天堂| 亚洲人成网站在线观看播放| 97热精品久久久久久| 黄片无遮挡物在线观看| 观看免费一级毛片| 国产av国产精品国产| 男女啪啪激烈高潮av片| 最近最新中文字幕大全电影3| 久久99热这里只频精品6学生| 久久久久久九九精品二区国产| 国产乱人偷精品视频| 欧美日韩视频高清一区二区三区二| 成人毛片60女人毛片免费| 精品久久久久久成人av| 乱人视频在线观看| 1000部很黄的大片| 日日啪夜夜爽| 久久99热这里只频精品6学生| 久久热精品热| 真实男女啪啪啪动态图| 日日摸夜夜添夜夜爱| 蜜臀久久99精品久久宅男| 国产色婷婷99| 久久韩国三级中文字幕| av播播在线观看一区| 中文字幕久久专区| 成人无遮挡网站| 日本午夜av视频| 久久久欧美国产精品| 老司机影院成人| 免费av毛片视频| 看免费成人av毛片| 三级国产精品片|