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

    地震資料的Hilbert Huang變換及其應(yīng)用

    2012-03-02 15:37:22劉浩杰中石化勝利油田分公司物探研究院山東東營257022
    石油天然氣學(xué)報 2012年12期
    關(guān)鍵詞:模態(tài)信號

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

    地震資料的Hilbert Huang變換及其應(yīng)用

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

    地震信號能量和頻率隨著時間的推移而改變,僅僅從時間域或頻率域來分析信號已不能充分地描述地震信號的特征。基于Hilbert-Huang變換 (HHT),把非穩(wěn)態(tài)的地震信號分解為不同分辨尺度的基本模態(tài)分量,利用不同模態(tài)分量的瞬時時頻特征精細(xì)刻畫地震信號的時間-頻率-能量之間的關(guān)系,探討了HHT在面波壓制方面的應(yīng)用。應(yīng)用結(jié)果表明,HHT在時頻域可以較好地達(dá)到衰減面波的目的。

    地震信號;Hilbert-Huang變換 (HHT);時頻分析;時間分辨率;頻率分辨率

    地震信號具有短時、突變等特點(diǎn),其能量和頻率隨著時間的推移而變化,僅僅從時間域或頻率域來分析信號已不能充分地描述信號的特征。時頻分析技術(shù)將地震數(shù)據(jù)變換到二維時間-頻率域,得到信號的頻譜、能量特征隨時間的變化,在時間-頻率域中研究地震資料的各種特征。時頻分析方法有很多,主要有短時傅里葉變換、小波變換、S變換以及廣義S變換等方法[1~3]。傳統(tǒng)的傅里葉變換建立在穩(wěn)態(tài)信號處理基礎(chǔ)上,分析信號總體所包含的各種頻率成分,但不能給出不同頻率成分的時間局部化信息;短時傅里葉變換以傅里葉變換為基礎(chǔ),通過窗函數(shù)作用,分析固定分辨率條件下信號的時頻特征;小波變換在時域和頻域都具有很好的局部化性質(zhì),提高了頻譜分辨率和時間定位精度,但本質(zhì)上仍是一種窗口可調(diào)的傅里葉變換;S變換和廣義S變換改進(jìn)了小波的形式,分辨率與頻率 (尺度)有關(guān),但本質(zhì)上仍是窗口可調(diào)的傅里葉變換,分辨率受窗口長度的限制。常規(guī)時頻分析方法在分辨率和信號平穩(wěn)性方面的限制影響了對地震信號時頻特征刻畫的精確性。

    Hilbert-Huang變換 (HHT)是一種新興的時頻分析方法,被認(rèn)為是近年來對以傅里葉變換為基礎(chǔ)進(jìn)行線性和穩(wěn)態(tài)譜分析的一個重大突破。它最初是由美國國家宇航局的Huang[4]于1998年提出的。該方法不受傅里葉分析的全局性限制,將任意信號依據(jù)信號本身的時間尺度特征,進(jìn)行多尺度時頻特征的刻畫,得到的時頻譜有較高的時間分辨率與頻率分辨率,也使得信號瞬時頻率具有實(shí)際的物理意義。HHT已應(yīng)用到地震力學(xué)、地球物理學(xué)等領(lǐng)域,并取得了較好的效果[5~7]。該次研究首先介紹了HHT的原理,研究了HHT頻譜和時頻譜,探討了HHT變換在面波壓制等方面的應(yīng)用。

    1 Hilbert-Huang變換 (HHT)

    HHT包括兩部分[4]:第1部分是經(jīng)驗?zāi)B(tài)分解 (empirical mode decomposition,EMD);第2部分是Hilbert譜分析 (Hilbert spectrum analysis,HSA)。

    為了利用HHT來分析非線性和非平穩(wěn)信號以獲得有明確物理意義的瞬時頻率,必須保證信號處理方法的完備性、正交性、局部性和自適應(yīng)性。基于此,Huang提出了本征模態(tài)函數(shù) (intrinsic mode function,IMF)的概念。本征模態(tài)函數(shù)必須滿足2個條件:一是在信號的有限時間內(nèi),其零交點(diǎn)的數(shù)目與極值點(diǎn)的數(shù)目必須相等或者至多相差一個;二是在任意時刻,由極大值得到的上包絡(luò)線和由極小值得到的下包絡(luò)線的平均值為零。非平穩(wěn)信號并不滿足IMF分量的兩個條件,因此必須先將信號分解成符合上述條件的IMF分量,然后再對各分量進(jìn)行HHT?;贗MF的定義,Huang提出了一種將任意信號分解成IMF分量的方法,即經(jīng)驗?zāi)B(tài)分解方法。其具體步驟如下:

    1)首先找出原始信號S(t)的所有的極大值點(diǎn)與極小值點(diǎn),然后分別用三次樣條函數(shù)插值擬合得到原始數(shù)據(jù)序列的上、下包絡(luò)線。計算上、下包絡(luò)線的均值得到平均包絡(luò)線m10(t),用S(t)減去m10(t)后,就得到一個去掉低頻的新序列h10(t),公式表述如下:

    一般來說,h10(t)不一定是符合IMF分量的條件,因此需要將其作為原始信號再重復(fù)上述過程。如果h10(t)的平均包絡(luò)線為m11(t),則去除該包絡(luò)線所包含的低頻成分后得到的序列為h11(t):

    重復(fù)上述過程,經(jīng)過k次循環(huán)后,得到平均包絡(luò)線m1k(t)。當(dāng)m1k(t)趨向于零時,此時的h1k(t)滿足IMF分量的條件,h1k(t)即為第1個IMF序列,將其定義為x1(t),它代表的是原始信號S(t)中頻率最高的成分。

    2)用S(t)減去x1(t),得到一個去掉高頻成分的新序列r1(t);對r1(t)進(jìn)行式(1)、(2)中的分解與計算,得到第2個IMF分量x2(t)。重復(fù)上述過程直到最后一個序列rn(t)成為一個單調(diào)函數(shù)不能再分解出IMF分量為止。上述過程可以表述為:

    圖1(a)為一個原始理論信號和經(jīng)過EMD分解得到的2個IMF分量信號。原始信號為一個頻率線性增加的正弦信號與單一頻率的余弦信號的疊加,其中正弦信號的頻率為線性的,而余弦信號的頻率為一固定值,公式為:S(t)=sin[2πt(6+6t)]+cos(2πt);經(jīng)過EMD分解后的信號分別相當(dāng)于原始信號的2個部分。從圖1(a)中可以看出,利用HHT可以準(zhǔn)確地分離變頻信號,表明了HHT對非線性、非穩(wěn)態(tài)信號處理的有效性。

    圖1 理論信號的IMF 分解及時間頻率曲線

    通過EMD分解得到的線性穩(wěn)態(tài)IMF分量非常適合做HHT,得到瞬時頻譜特征,即HHT譜。HHT是一種線性變換,如果輸入信號是平穩(wěn)的,那么輸出信號也應(yīng)該是平穩(wěn)的,且HHT強(qiáng)調(diào)局部屬性,這避免了傅里葉變換為擬合原序列而產(chǎn)生的許多多余的、事實(shí)上并不存在的高、低頻成分,也賦予了信號瞬時頻率分量的具體物理意義。圖1(b)為對圖1(a)中信號的IMF分量作HHT后得到的時間-頻率曲線。從圖1(b)中可以看出,通過HHT能夠得到準(zhǔn)確的瞬時頻率,其中頻率隨時間推移線性增長的是第1個IMF分量 (IMF1)的瞬時頻率,另外一個是第2個IMF分量 (IMF2)的瞬時頻率值。

    2 HHT邊際譜和時頻特征分析

    得到信號本征模態(tài)函數(shù)之后,就可以利用HHT得到每個本征模態(tài)信號的瞬時頻率。如果忽略HHT的殘余量,信號就可以表示成如下形式:

    式中,X(t)為原始信號;a(t)為分解的IMF分量;下標(biāo)j為IMF分量的序號,j=1,2,…,n;ω(t)為特定頻率,Hz;t為時間,s。

    式(4)把信號的幅度和頻率統(tǒng)一表示為時間的函數(shù)。信號的傅里葉變換表達(dá)式為:

    式(4)與式(5)比較表明,信號的HHT實(shí)質(zhì)上是傅里葉變換的擴(kuò)展。而且,通過HHT能夠很容易地實(shí)現(xiàn)信號幅度和頻率的調(diào)整,這種變化的幅度和瞬時頻率不僅極大地改進(jìn)了變換的效率,也使得信號的擴(kuò)展適應(yīng)于非穩(wěn)態(tài)信號。因此,這種擴(kuò)展突破了傅里葉變換固定幅度和固定頻率的限制。

    2.1 HHT邊際譜分析

    HHT能夠以三維圖的形式表示幅度、瞬時頻率和時間的關(guān)系。這種幅度的頻率-時間分布稱之為Hilbert譜,幅度的頻率分布定義為HHT邊際譜:

    式中,h(ω)為HHT邊際譜;H(ω,t)為Hilbert譜;T為信號持續(xù)時間。

    這種HHT邊際譜所表示的意義與傅里葉譜 (FFT振幅譜)不完全一樣。對于FFT振幅譜,特定頻率ω處的能量是表示該頻率的正弦或余弦波在整個信號持續(xù)時間范圍內(nèi)的貢獻(xiàn);而對于HHT邊際譜,特定頻率ω處的能量是表示該頻率的振動在整個信號持續(xù)時間范圍內(nèi)的貢獻(xiàn)。也就是說,F(xiàn)FT振幅譜是把信號分解為不同分量的正弦或余弦信號后進(jìn)行能量分布的分析,而HHT邊際譜是直接分析信號瞬時頻率的能量分布。因此,通過瞬時頻率的引入,使得HHT邊際譜有明確的物理意義,而且也去除了常規(guī)FFT振幅譜為擬合原始信號而引入的不存在的諧波成分。

    圖2為圖1中信號HHT邊際譜與FFT振幅譜的比較。從圖2中可以看出,由于原始信號是正、余弦信號的疊加,所以HHT邊際譜和FFT振幅譜基本一致,只存在一些細(xì)微的差別。

    圖3(a)為實(shí)際地震道,采樣間隔為2ms;圖3(b)為其對應(yīng)的HHT邊際譜和FFT振幅譜。從對比中可以看出,實(shí)際地震信號的HHT邊際譜與FFT振幅譜存在較大差別。從物理意義上,HHT邊際譜更準(zhǔn)確地表明了不同頻率分量在信號持續(xù)時間內(nèi)的作用。

    圖2 HHT邊際譜與FFT振幅譜比較

    圖3 實(shí)際地震道及HHT邊際譜和FFT振幅譜的比較

    2.2 HHT時頻特征分析

    如式 (4)所述,HHT可得到能量關(guān)于時間、頻率的函數(shù),即進(jìn)行時頻特征分析。HHT作為一種新興的時頻分析技術(shù),已經(jīng)在許多領(lǐng)域得到了應(yīng)用。為了說明HHT時頻分析技術(shù)的特征,將其與幾種常規(guī)時頻分析方法進(jìn)行了對比。

    選取圖3所示的實(shí)際地震道作為試驗數(shù)據(jù),對該道數(shù)據(jù)分別做短時傅里葉變換 (STFT)、連續(xù)小波變換 (CWT)、廣義S變換 (GST)和HHT,得到其時頻譜圖 (圖4)。比較而言,3種常規(guī)時頻分析方法在地震頻帶范圍內(nèi)的效果相似,盡管能表示出不同頻率分量能量隨時間的變化趨勢,但其時頻分辨率仍然比較低。與常規(guī)時頻分析相比,HHT時頻分析的時間分辨率還是頻率分辨率都有較大的提高,能夠很清楚地刻畫不同頻率分量能量的時間-頻率分布。

    圖4 實(shí)際地震道的常規(guī)時頻分析和HHT時頻分析比較

    3 HHT應(yīng)用探討

    在地震勘探中,作為一種規(guī)則干擾波,面波視速度較低且存在頻散現(xiàn)象,具有能量強(qiáng)、頻率低、振動延續(xù)時間長等特點(diǎn),在炮集上呈線性分布,降低了地震資料的信噪比。如何有效去除面波一直是地震勘探領(lǐng)域的難點(diǎn)。

    對疊前炮集記錄的HHT分析表明,信號不同模態(tài)的分量能夠精細(xì)刻畫某一時段和頻段的不同信號分量,即能夠清晰地反映有效信號和干擾信號的能量以及頻率隨時間的變化。因此,可以在HHT基礎(chǔ)上設(shè)計時頻濾波器達(dá)到壓制面波的目的。首先,對原始炮集地震記錄進(jìn)行EMD分解,得到不同模態(tài)信號的IMF分量;然后對所有分量分別作HHT,得到時頻譜圖,確定面波所在頻段;結(jié)合常規(guī)面波消除方法,設(shè)計時頻濾波器對包含面波頻率成分的分量作處理或直接去除,得到新的IMF分量;最后將所有新得到的IMF分量相加得到去除面波的地震記錄。

    圖5為面波壓制前、后的效果圖,經(jīng)過對比可以看到,壓制面波后的剖面基本不存在面波成分,而且有效波也得到了保留;有效反射波同相軸的連續(xù)性增強(qiáng),信噪比提高。因此,HHT在時頻域壓制面波方面可以較好地達(dá)到衰減面波的目的,并有很好的效果。

    圖5 疊前炮集數(shù)據(jù)壓制面波前 (a)、后 (b)的對比圖

    4 結(jié) 論

    在HHT理論分析的基礎(chǔ)上,利用理論信號和實(shí)際地震資料,研究了HHT邊際譜和時頻譜的特點(diǎn),證實(shí)了HHT理論的完備性以及物理意義的明確性,探討了HHT在地震勘探領(lǐng)域的應(yīng)用,得到了以下認(rèn)識:

    1)地震信號的HHT通過對非平穩(wěn)地震信號的線性化,建立了以瞬時頻率為表征信號變化的基本量,以固有模態(tài)信號為基本時域信號的新時頻分析方法體系,實(shí)現(xiàn)了具有物理意義的地震信號的線性瞬時時頻分析,更加精細(xì)刻畫了地震信號的時間-頻率-能量之間的關(guān)系。這對于從地震資料中提取反映地層特征信息的高精度和高分辨屬性參數(shù)具有重要的意義。

    2)HHT邊際譜與FFT振幅譜所表示的意義不完全一樣。對于傅里葉譜,特定頻率處的能量是表示該頻率的正弦或余弦波在整個信號持續(xù)時間范圍內(nèi)的貢獻(xiàn);而對于HHT邊際譜,特定頻率處的能量是表示該頻率的振動在整個信號持續(xù)時間范圍內(nèi)的貢獻(xiàn)。也就是說,通過瞬時頻率的引入,使得HHT邊際譜有明確的物理意義,也去除了常規(guī)傅里葉譜為擬合原始信號而引入的不存在的諧波成分。

    3)HHT的多尺度和高分辨特性使得該變換在地震勘探領(lǐng)域有潛在的重要應(yīng)用價值。

    [1]王培茂 .地震信號的時頻特征表示方法及應(yīng)用 [D].長春:吉林大學(xué),2008.

    [2]鄒文,陳愛萍,顧漢明 .地震信號的時頻分析方法 [J].世界地質(zhì),2004,23(1):95~99.

    [3]李振春,刁瑞,韓文功,等 .線性時頻分析方法綜述 [J].勘探地球物理進(jìn)展,2010,33(4):239~246.

    [4]Huang N E,Shen Zheng,Long S R,etal.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J].Proc R Soc Lond A,1998,454(1971):903~906.

    [5]張郁山 .希爾伯特-黃變換 (HHT)與地震動時程的希爾伯特譜 [D].北京:中國地震局地球物理研究所,2003.

    [6]胡勁松,楊世錫,任達(dá)千 .一種基于EMD的振動信號時頻分析新方法研究 [J].振動與沖擊,2008,27(8):71~73.

    [7]公茂盛,謝禮立 .HHT方法在地震工程中的應(yīng)用之初步探討 [J].世界地震工程,2003,19(3):39~43.Dongying257022,Shan-

    [編輯] 龍 舟

    Hilbert-Huang Transform for Seismic Data and Its Application

    LIU Hao-jie (Author' s Address:Geophysical Research Institute,Shengli Oilfield Company,SINOPEC,dong,China)

    The energy and the frequency of seismic signal were varied with time,so it was difficult to characterize seismic signal accurately and subtly only in the field of time or frequency.Based on Hilbert-Huang transform,the non-stationary seismic signal could be decomposed into a finite and small number of intrinsic mode functions and the relation between time-frequency-energy of seismic signal could then be identified accurately by using the instantaneous time-frequency analysis of different model signal.The application of Hilbert-Huang transform in compressing wave is studied.The result shows that Hilbert-Huang transform can attenuate surface wave accurately in time-frequency field.

    seismic signal;Hilbert-Huang transform;time-frequency analysis;time resolution;frequency resolution

    P631.44

    A

    1000-9752(2012)12-0076-05

    2012-08-11

    國家 “863”計劃項目 (2011AA060302)。

    劉浩杰 (1976-),男,1998年江漢石油學(xué)院畢業(yè),博士,研究員,現(xiàn)主要從事油藏綜合地球物理研究。

    猜你喜歡
    模態(tài)信號
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    孩子停止長個的信號
    車輛CAE分析中自由模態(tài)和約束模態(tài)的應(yīng)用與對比
    基于LabVIEW的力加載信號采集與PID控制
    國內(nèi)多模態(tài)教學(xué)研究回顧與展望
    一種基于極大似然估計的信號盲抽取算法
    高速顫振模型設(shè)計中顫振主要模態(tài)的判斷
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    由單個模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
    計算物理(2014年2期)2014-03-11 17:01:39
    国产日韩欧美视频二区| 丝袜脚勾引网站| 男男h啪啪无遮挡| 色婷婷久久久亚洲欧美| 另类精品久久| 国产在线一区二区三区精| 亚洲国产成人一精品久久久| 老汉色∧v一级毛片| 秋霞在线观看毛片| 国产精品久久久久久精品古装| 亚洲 国产 在线| 亚洲图色成人| 美女午夜性视频免费| 亚洲视频免费观看视频| a级片在线免费高清观看视频| 十八禁网站网址无遮挡| 久久av网站| 夫妻性生交免费视频一级片| 精品国产一区二区三区久久久樱花| 久久免费观看电影| 久久九九热精品免费| 国产精品三级大全| 一级黄片播放器| 午夜福利乱码中文字幕| 亚洲av日韩在线播放| 一级毛片 在线播放| 在线观看免费高清a一片| 999精品在线视频| 国产一卡二卡三卡精品| 午夜av观看不卡| 久久影院123| 在线 av 中文字幕| 一区二区日韩欧美中文字幕| 国产黄色免费在线视频| 一级毛片 在线播放| 69精品国产乱码久久久| 国产又色又爽无遮挡免| 亚洲少妇的诱惑av| 久久天堂一区二区三区四区| 91精品三级在线观看| 色网站视频免费| 亚洲黑人精品在线| 亚洲欧美精品自产自拍| 香蕉丝袜av| 欧美激情 高清一区二区三区| 青春草视频在线免费观看| 精品人妻在线不人妻| 欧美变态另类bdsm刘玥| 女性生殖器流出的白浆| 色精品久久人妻99蜜桃| 国产精品久久久久久人妻精品电影 | 午夜福利免费观看在线| 久久久久精品人妻al黑| 19禁男女啪啪无遮挡网站| 在线天堂中文资源库| 成年人免费黄色播放视频| 成在线人永久免费视频| 大型av网站在线播放| 日韩av在线免费看完整版不卡| 国产女主播在线喷水免费视频网站| 成年美女黄网站色视频大全免费| 久久99一区二区三区| 欧美精品人与动牲交sv欧美| 男女边摸边吃奶| 两性夫妻黄色片| 观看av在线不卡| 一级,二级,三级黄色视频| 交换朋友夫妻互换小说| 免费观看av网站的网址| 成人亚洲欧美一区二区av| 亚洲欧美精品自产自拍| 国产精品 国内视频| 国产精品国产三级专区第一集| 成年美女黄网站色视频大全免费| cao死你这个sao货| 一区二区三区精品91| 欧美人与性动交α欧美精品济南到| 国产成人av教育| 亚洲精品乱久久久久久| 国产精品欧美亚洲77777| 丰满迷人的少妇在线观看| 亚洲精品国产色婷婷电影| 天堂8中文在线网| 国产成人免费无遮挡视频| 在现免费观看毛片| 亚洲欧美激情在线| 飞空精品影院首页| 两个人免费观看高清视频| 人人妻人人爽人人添夜夜欢视频| 国产伦人伦偷精品视频| 久久精品久久久久久噜噜老黄| 亚洲欧洲日产国产| 国产日韩一区二区三区精品不卡| 亚洲av片天天在线观看| 国产精品久久久久久精品电影小说| 丰满少妇做爰视频| 久久精品久久久久久久性| 日韩 亚洲 欧美在线| 欧美 日韩 精品 国产| 97精品久久久久久久久久精品| 午夜两性在线视频| 亚洲 国产 在线| 久久国产精品男人的天堂亚洲| 欧美av亚洲av综合av国产av| 日日摸夜夜添夜夜爱| 美女福利国产在线| 五月天丁香电影| 午夜福利影视在线免费观看| 久久久久国产精品人妻一区二区| 欧美乱码精品一区二区三区| 我要看黄色一级片免费的| 久久热在线av| 91九色精品人成在线观看| 久久久精品区二区三区| 秋霞在线观看毛片| 97精品久久久久久久久久精品| 九色亚洲精品在线播放| 老司机影院成人| 久久精品久久精品一区二区三区| 一区二区三区乱码不卡18| 麻豆乱淫一区二区| 亚洲精品国产av蜜桃| 美女视频免费永久观看网站| 97精品久久久久久久久久精品| 国产爽快片一区二区三区| 日韩电影二区| avwww免费| 欧美97在线视频| 久久久欧美国产精品| 夫妻性生交免费视频一级片| 欧美精品一区二区大全| 亚洲精品美女久久av网站| 国产在线视频一区二区| 熟女av电影| 国产精品.久久久| 免费少妇av软件| 精品一品国产午夜福利视频| 国产日韩欧美在线精品| 亚洲精品一二三| 极品少妇高潮喷水抽搐| 亚洲欧美色中文字幕在线| 在线观看免费视频网站a站| 亚洲人成网站在线观看播放| 欧美日韩综合久久久久久| 中文字幕人妻丝袜一区二区| 成年人免费黄色播放视频| 99精国产麻豆久久婷婷| 制服诱惑二区| 成人三级做爰电影| av又黄又爽大尺度在线免费看| 狠狠婷婷综合久久久久久88av| 美女脱内裤让男人舔精品视频| 日韩av不卡免费在线播放| 亚洲午夜精品一区,二区,三区| 另类精品久久| 久久久国产欧美日韩av| 一边摸一边做爽爽视频免费| 成人国产一区最新在线观看 | 色精品久久人妻99蜜桃| 亚洲欧美成人综合另类久久久| 男女下面插进去视频免费观看| 免费高清在线观看视频在线观看| 日韩伦理黄色片| 国产男女内射视频| 国产高清国产精品国产三级| 久久综合国产亚洲精品| 午夜91福利影院| 亚洲欧美激情在线| 精品人妻1区二区| 狠狠精品人妻久久久久久综合| h视频一区二区三区| 欧美人与善性xxx| 日韩一本色道免费dvd| 最新的欧美精品一区二区| 肉色欧美久久久久久久蜜桃| 日韩av在线免费看完整版不卡| 成人国产av品久久久| 亚洲精品国产av蜜桃| 久久天堂一区二区三区四区| 肉色欧美久久久久久久蜜桃| 男女边摸边吃奶| 久久免费观看电影| 99香蕉大伊视频| 黑人猛操日本美女一级片| 真人做人爱边吃奶动态| 国产亚洲一区二区精品| 国产精品久久久久久精品古装| 精品人妻一区二区三区麻豆| 日本黄色日本黄色录像| 嫁个100分男人电影在线观看 | 久久精品亚洲熟妇少妇任你| 亚洲国产精品成人久久小说| 欧美成人午夜精品| av天堂在线播放| 考比视频在线观看| 国产真人三级小视频在线观看| 欧美老熟妇乱子伦牲交| 成在线人永久免费视频| 一区在线观看完整版| 一级,二级,三级黄色视频| 免费人妻精品一区二区三区视频| 汤姆久久久久久久影院中文字幕| 免费在线观看黄色视频的| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜av观看不卡| 欧美日韩亚洲国产一区二区在线观看 | 嫩草影视91久久| a级片在线免费高清观看视频| 久久女婷五月综合色啪小说| 中文字幕人妻熟女乱码| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品一区三区| 热99国产精品久久久久久7| 日韩制服丝袜自拍偷拍| 两个人免费观看高清视频| 国产精品一区二区在线观看99| 嫁个100分男人电影在线观看 | 国产精品久久久久成人av| 亚洲精品在线美女| 日韩欧美一区视频在线观看| 波野结衣二区三区在线| 黑丝袜美女国产一区| 国产精品九九99| 国产精品国产三级专区第一集| 少妇的丰满在线观看| netflix在线观看网站| 国产深夜福利视频在线观看| 老汉色av国产亚洲站长工具| 少妇裸体淫交视频免费看高清 | 国产精品九九99| 在线观看免费日韩欧美大片| 中文字幕最新亚洲高清| 亚洲精品国产色婷婷电影| 国产精品久久久久久精品电影小说| 亚洲伊人色综图| 国产成人精品在线电影| 亚洲色图 男人天堂 中文字幕| 这个男人来自地球电影免费观看| 精品国产乱码久久久久久小说| 国产成人av激情在线播放| 免费看av在线观看网站| 国产一卡二卡三卡精品| videosex国产| 国产精品一国产av| 久久久久久久久免费视频了| 一级黄色大片毛片| 在线观看一区二区三区激情| 国产亚洲精品久久久久5区| 成人手机av| 美女主播在线视频| 国产精品人妻久久久影院| 99九九在线精品视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲人成77777在线视频| 女警被强在线播放| 欧美日韩亚洲国产一区二区在线观看 | 黑人猛操日本美女一级片| 亚洲欧洲国产日韩| 免费一级毛片在线播放高清视频 | 亚洲一区中文字幕在线| 国产亚洲精品第一综合不卡| 夫妻性生交免费视频一级片| 欧美日韩成人在线一区二区| 视频区图区小说| 国产成人免费无遮挡视频| 国产爽快片一区二区三区| videosex国产| 国产97色在线日韩免费| 高潮久久久久久久久久久不卡| 97精品久久久久久久久久精品| 赤兔流量卡办理| 亚洲国产欧美一区二区综合| 国产精品欧美亚洲77777| 国产精品熟女久久久久浪| 午夜精品国产一区二区电影| 91麻豆精品激情在线观看国产 | 亚洲欧美激情在线| 国产成人一区二区在线| 欧美日韩视频精品一区| 99国产精品99久久久久| 亚洲精品一区蜜桃| 一二三四社区在线视频社区8| 国产女主播在线喷水免费视频网站| 亚洲伊人久久精品综合| 国语对白做爰xxxⅹ性视频网站| av在线app专区| 宅男免费午夜| 91麻豆精品激情在线观看国产 | 久久精品亚洲av国产电影网| 久久久久视频综合| 亚洲av片天天在线观看| 老司机在亚洲福利影院| 久久精品aⅴ一区二区三区四区| 丝袜人妻中文字幕| 成人国产av品久久久| 精品一区二区三卡| 黄频高清免费视频| 亚洲国产av新网站| 一级毛片女人18水好多 | 精品人妻一区二区三区麻豆| 精品亚洲乱码少妇综合久久| 久久亚洲精品不卡| 1024视频免费在线观看| av一本久久久久| 日韩 亚洲 欧美在线| 亚洲av美国av| 后天国语完整版免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 免费日韩欧美在线观看| 日韩 欧美 亚洲 中文字幕| 国产免费福利视频在线观看| 午夜激情久久久久久久| 在线观看免费高清a一片| 777久久人妻少妇嫩草av网站| 久久毛片免费看一区二区三区| 后天国语完整版免费观看| 国产成人精品久久二区二区免费| 亚洲专区中文字幕在线| 国产又爽黄色视频| 男女免费视频国产| www.精华液| 一区二区三区激情视频| 日韩熟女老妇一区二区性免费视频| 亚洲国产av影院在线观看| 欧美精品av麻豆av| 亚洲一区中文字幕在线| 精品人妻1区二区| 交换朋友夫妻互换小说| 制服人妻中文乱码| 女人被躁到高潮嗷嗷叫费观| 久久久精品免费免费高清| 狠狠婷婷综合久久久久久88av| 免费不卡黄色视频| 亚洲一区中文字幕在线| 欧美黄色片欧美黄色片| 国产精品偷伦视频观看了| 脱女人内裤的视频| 亚洲欧美中文字幕日韩二区| h视频一区二区三区| 国产一区二区三区av在线| 日日爽夜夜爽网站| 狠狠精品人妻久久久久久综合| 国产欧美亚洲国产| 欧美精品一区二区免费开放| 丁香六月欧美| av在线老鸭窝| 日韩大码丰满熟妇| 国产1区2区3区精品| 水蜜桃什么品种好| 波多野结衣av一区二区av| 一级毛片黄色毛片免费观看视频| 日韩制服骚丝袜av| 一级毛片女人18水好多 | 99精国产麻豆久久婷婷| 亚洲伊人久久精品综合| 国产片特级美女逼逼视频| 1024香蕉在线观看| 欧美日韩一级在线毛片| 亚洲国产最新在线播放| 女人被躁到高潮嗷嗷叫费观| 亚洲av电影在线进入| 一级毛片电影观看| 久久久久久久大尺度免费视频| 精品国产超薄肉色丝袜足j| 久久久久久人人人人人| 亚洲中文av在线| 国产麻豆69| a级毛片黄视频| 一二三四社区在线视频社区8| 永久免费av网站大全| 国产精品免费大片| 国产亚洲一区二区精品| 90打野战视频偷拍视频| 2018国产大陆天天弄谢| 水蜜桃什么品种好| 成人午夜精彩视频在线观看| 国产精品偷伦视频观看了| 久久精品久久精品一区二区三区| 美女中出高潮动态图| 国产熟女午夜一区二区三区| 午夜福利影视在线免费观看| 国产精品麻豆人妻色哟哟久久| 免费高清在线观看日韩| 国产成人啪精品午夜网站| 婷婷成人精品国产| 亚洲欧美一区二区三区黑人| 黄色a级毛片大全视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品av久久久久免费| 国产男女内射视频| 纯流量卡能插随身wifi吗| 国产亚洲精品久久久久5区| 久久影院123| 97在线人人人人妻| 国产成人影院久久av| 99九九在线精品视频| 考比视频在线观看| 久久精品久久久久久噜噜老黄| 国产在视频线精品| 热99久久久久精品小说推荐| 一区二区av电影网| 老司机在亚洲福利影院| 99国产精品免费福利视频| 日本vs欧美在线观看视频| 侵犯人妻中文字幕一二三四区| 欧美日韩综合久久久久久| 日日爽夜夜爽网站| 国产精品国产av在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 伊人久久大香线蕉亚洲五| 免费高清在线观看视频在线观看| 欧美 亚洲 国产 日韩一| 久久鲁丝午夜福利片| 国产成人啪精品午夜网站| 精品亚洲成a人片在线观看| 波多野结衣一区麻豆| 久热爱精品视频在线9| 免费女性裸体啪啪无遮挡网站| 日日爽夜夜爽网站| 亚洲第一青青草原| 国产亚洲欧美精品永久| 韩国精品一区二区三区| 中文欧美无线码| 日韩 亚洲 欧美在线| 国产伦人伦偷精品视频| 亚洲中文日韩欧美视频| 国产精品二区激情视频| 免费久久久久久久精品成人欧美视频| 黑人巨大精品欧美一区二区蜜桃| 欧美xxⅹ黑人| 精品国产国语对白av| 欧美日本中文国产一区发布| 久久ye,这里只有精品| 精品国产一区二区三区四区第35| 国产1区2区3区精品| 一级毛片女人18水好多 | 首页视频小说图片口味搜索 | 亚洲av成人精品一二三区| 久久精品久久久久久久性| 精品福利永久在线观看| 午夜91福利影院| 国产在线免费精品| 在线观看国产h片| avwww免费| e午夜精品久久久久久久| 在线亚洲精品国产二区图片欧美| 国产黄频视频在线观看| 少妇 在线观看| 大话2 男鬼变身卡| 久久精品国产综合久久久| 岛国毛片在线播放| 9色porny在线观看| 午夜福利乱码中文字幕| 成年女人毛片免费观看观看9 | 纯流量卡能插随身wifi吗| 久久99精品国语久久久| 亚洲精品第二区| 男的添女的下面高潮视频| 亚洲成人免费电影在线观看 | 国产日韩欧美在线精品| 欧美精品av麻豆av| 看十八女毛片水多多多| 亚洲一区中文字幕在线| 制服诱惑二区| 极品少妇高潮喷水抽搐| 黑人猛操日本美女一级片| 在线精品无人区一区二区三| 一级毛片我不卡| 纵有疾风起免费观看全集完整版| 亚洲av综合色区一区| 在线观看免费午夜福利视频| 少妇精品久久久久久久| 青春草亚洲视频在线观看| 五月开心婷婷网| 99国产精品一区二区三区| 亚洲av国产av综合av卡| 亚洲 国产 在线| 国产男人的电影天堂91| 午夜福利,免费看| 精品卡一卡二卡四卡免费| 亚洲av日韩在线播放| 亚洲中文字幕日韩| 性少妇av在线| 国产午夜精品一二区理论片| 国产1区2区3区精品| 日本a在线网址| 国产视频一区二区在线看| 欧美日韩福利视频一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av电影在线观看一区二区三区| 五月开心婷婷网| 后天国语完整版免费观看| 97在线人人人人妻| 丰满人妻熟妇乱又伦精品不卡| 91精品三级在线观看| 人妻人人澡人人爽人人| 18在线观看网站| 成人国产一区最新在线观看 | 亚洲男人天堂网一区| 亚洲国产欧美网| 电影成人av| 女人被躁到高潮嗷嗷叫费观| 中文精品一卡2卡3卡4更新| 两个人免费观看高清视频| 亚洲精品av麻豆狂野| 黄频高清免费视频| 热re99久久精品国产66热6| 日韩一卡2卡3卡4卡2021年| 久久精品aⅴ一区二区三区四区| 国产99久久九九免费精品| 性色av乱码一区二区三区2| 亚洲成人免费电影在线观看 | 欧美变态另类bdsm刘玥| 美女福利国产在线| 国产精品 国内视频| 久久久久久久国产电影| 国产精品免费大片| 国产一区二区激情短视频 | 日韩精品免费视频一区二区三区| 校园人妻丝袜中文字幕| 国产成人精品久久久久久| 亚洲中文字幕日韩| 欧美黄色淫秽网站| 午夜福利,免费看| 国语对白做爰xxxⅹ性视频网站| 亚洲av国产av综合av卡| 日本五十路高清| 欧美精品av麻豆av| 99国产综合亚洲精品| 欧美成人午夜精品| 亚洲免费av在线视频| 欧美激情 高清一区二区三区| 欧美日韩国产mv在线观看视频| 久久精品成人免费网站| 亚洲av成人不卡在线观看播放网 | 欧美精品亚洲一区二区| 国产欧美日韩综合在线一区二区| 久久精品熟女亚洲av麻豆精品| 日日夜夜操网爽| 大型av网站在线播放| 国产高清视频在线播放一区 | 国产精品99久久99久久久不卡| 国产欧美日韩综合在线一区二区| 久久精品aⅴ一区二区三区四区| 久久国产亚洲av麻豆专区| 一级片免费观看大全| 国产欧美日韩一区二区三区在线| 熟女少妇亚洲综合色aaa.| 欧美日本中文国产一区发布| 男人爽女人下面视频在线观看| 99精品久久久久人妻精品| 男女无遮挡免费网站观看| 啦啦啦视频在线资源免费观看| 男的添女的下面高潮视频| 好男人电影高清在线观看| 男女边吃奶边做爰视频| 免费高清在线观看日韩| 黄网站色视频无遮挡免费观看| 国产精品亚洲av一区麻豆| 亚洲久久久国产精品| 国产精品熟女久久久久浪| 美女主播在线视频| 日韩欧美一区视频在线观看| 好男人视频免费观看在线| 久久99精品国语久久久| 在线观看一区二区三区激情| 亚洲,一卡二卡三卡| 一级片免费观看大全| av线在线观看网站| 久久国产精品人妻蜜桃| 欧美xxⅹ黑人| 久久久久久久国产电影| 亚洲精品国产色婷婷电影| 国产有黄有色有爽视频| 激情视频va一区二区三区| 精品欧美一区二区三区在线| 19禁男女啪啪无遮挡网站| 亚洲欧美精品综合一区二区三区| 极品人妻少妇av视频| 国产91精品成人一区二区三区 | 少妇猛男粗大的猛烈进出视频| 日韩av在线免费看完整版不卡| videos熟女内射| 欧美日韩视频高清一区二区三区二| 精品免费久久久久久久清纯 | 各种免费的搞黄视频| 亚洲国产日韩一区二区| a 毛片基地| 操美女的视频在线观看| 91精品三级在线观看| 老司机在亚洲福利影院| 99久久精品国产亚洲精品| 中文字幕人妻丝袜一区二区| 亚洲欧美中文字幕日韩二区| 亚洲国产av影院在线观看| 国产成人一区二区在线| 一边摸一边抽搐一进一出视频| 国产麻豆69| 国产91精品成人一区二区三区 | 50天的宝宝边吃奶边哭怎么回事| 人人妻人人澡人人爽人人夜夜| 国产老妇伦熟女老妇高清| 欧美日韩福利视频一区二区| 国产不卡av网站在线观看| 亚洲av片天天在线观看| 欧美精品一区二区大全| 亚洲精品中文字幕在线视频| 精品欧美一区二区三区在线| videosex国产| 妹子高潮喷水视频| 国产亚洲午夜精品一区二区久久| 美女脱内裤让男人舔精品视频| 每晚都被弄得嗷嗷叫到高潮| 男女高潮啪啪啪动态图| 精品免费久久久久久久清纯 |