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

    運用小波變換方法研究華北地區(qū)的地震活動性*

    2011-01-09 02:58:38張小濤馮向東張躍剛張雙鳳楊雅瓊韓麗萍
    地震研究 2011年1期
    關鍵詞:小波基華北地區(qū)時間尺度

    張小濤,馮向東,張躍剛,張雙鳳,楊雅瓊,韓麗萍

    (1.河北省地震局邯鄲中心臺,河北邯鄲 056001;2.吉林大學地球探測科技與技術學院,長春 130026;3.河北省地震局,石家莊 050021)

    運用小波變換方法研究華北地區(qū)的地震活動性*

    張小濤1,2,馮向東3,張躍剛3,張雙鳳1,楊雅瓊3,韓麗萍1

    (1.河北省地震局邯鄲中心臺,河北邯鄲 056001;2.吉林大學地球探測科技與技術學院,長春 130026;3.河北省地震局,石家莊 050021)

    運用小波變換方法,以不同的小波時間尺度分析了華北地區(qū) 1480~2009年M≥5.0和 1970~2009年M≥4.0地震的貝尼奧夫應變數(shù)據(jù),得到了不同時間尺度的地震活躍期和平靜期。分析表明:目前華北地區(qū)從 1480年以來的M≥5.0地震活動剛過活躍期,處于應力積累階段,而 1970年以來M≥4.0地震存在 7年左右的活動周期,目前正在進入 4級以上地震的活躍期。

    貝尼奧夫應變;小波變換;地震活躍期;地震平靜期;華北地區(qū)

    0 引言

    表征地震活動性的各種參數(shù)從不同角度對地震活動的特征進行了完整、系統(tǒng)的描述,成為地震預測預報的重要參考依據(jù),但是地震活動在時間軸上的分布是否有規(guī)律可循,多年來一直是個有爭議的問題。馬宗晉等(1982)、張國民等(1987,2000)、張培震等(2003)對我國的地震活動作了分析,發(fā)現(xiàn)在不同的時期,地震活動顯示了平靜和活躍變化;同一地震區(qū)內(nèi)的地震活動具有某種程度的自相似性,在時間分布上具有某種“似周期性”分布特征;在不同時期地震發(fā)生的頻度和強度有較明顯的不同(時振梁等,1997)。Benioff(1951)曾把 1900~1950年全球的地震活動分為 5個時段,指出地震活動存在地震活躍期和相對平靜期;馬宗晉等(1982)指出地震活動強弱隨時間分段現(xiàn)象在我國大陸地區(qū)是相當普遍的,而且強弱變化的時間范圍長短不一;張國民和張培震(2000)、傅征祥等(2000)討論了華北地區(qū)M≥7.0強震的時間分布,發(fā)現(xiàn)在每一個地震活躍期內(nèi)大地震的積累頻度曲線隨時間的變化大致呈指數(shù)函數(shù)形式;萬永革等(2003)也分析了中國大陸地區(qū)M≥5.0地震活動的周期性,得到了不同時間尺度的地震平靜期和地震活躍期;劉虹和蘇有錦(2007)研究了 1970年以來云南地區(qū)MS≥7.0地震的各種活動參數(shù),發(fā)現(xiàn)震源區(qū)地震活動水平在大地震前后有很大的差異;蘇有錦等(2001)研究了 20世紀以來云南地區(qū)MS≥5.0地震的活躍期與平靜期的交替過程??傊?活躍期和平靜期的劃分對地震預測預報,特別是強震的中短期分析預報甚至短臨預報工作具有十分重要的意義,是我國地震預測研究中一個不可缺少的環(huán)節(jié)(馬宏生等,2002,2005)。

    華北地區(qū)作為中國內(nèi)陸的地震多發(fā)區(qū),一直備受地震學家的關注。如何從不同尺度上劃分活躍期和平靜期,并在大尺度活躍期(或平靜期)中提取小尺度的活躍期和平靜期,進一步為地震預測預報提供依據(jù)?在傳統(tǒng)的M-t圖上難以找到地震活動演化過程突變點。小波分析方法具有能夠根據(jù)分析對象不同而自動調(diào)整有關參數(shù)的 “自適應性”和能夠根據(jù)觀測對象不同而自動“調(diào)焦”的特性,因而被廣泛應用于地震研究的各個領域(劉希強等,1998;吳忠良,1999;邵輝成等,2000;徐義賢,王家映,2000)。筆者采用小波變換方法,從不同時間尺度和震級尺度研究了華北地區(qū)的地震在時間上的分布規(guī)律,以期為地震預測預報提供參考。

    1 小波變換的原理

    一個地震時間序列 f(t)是由不同強度的地震構(gòu)成的,在 t時刻地震的強度可以表示為

    其中 a為奇異性指數(shù)或標度指數(shù)。奇異性是指 t→0時 a的數(shù)值(特別是 a<0)決定了 f(t)→0的速度。標度性是指若變量 t變成λt,則 f(t)和 f(λt)是自相似的,即

    式(2)常稱為標度律。自相似的標度律是多種尺度(即無特征尺度)現(xiàn)象的特征。地震資料的時間序列含有多種不同尺度。直觀地看 f(t)是雜亂無章的,但小波變換好比顯微鏡,調(diào)節(jié)其中的放大倍數(shù)就可以清楚地看出 f(t)在各個層次上的變化趨勢。小波變換 Tg是將 f(t)分解成具有局部特性的小波g(t):

    設θ(t)為具有低通性質(zhì)的平滑函數(shù),以它的一階、二階導數(shù)作為小波對 f(t)作小波變換,可以證明(Logan,1977)

    式中,*表示卷積。也就是說,用φ1(t)、φ2(t)對 f(t)作小波變換,分別相當于 f(t)被θ(t)平滑后對 t求一階或二階導數(shù)。因此,對某一固定 a值,f(t)*θa(t)的拐點既是Wφ1f(a,τ)的極值點,又是Wφ2f(a,τ)的過零點,由此可檢測地震活動信號的急劇變化之處(Mallat,1991;Mallat,Hwang,1992)。

    本研究中采用具有低通濾波性質(zhì)的平滑高斯函數(shù),根據(jù)它的一階導數(shù)、二階導數(shù)作為小波基函數(shù)進行突變點分析。小波函數(shù)的表達式為

    選擇該函數(shù)是因為它對稱、可微、可積,時域、頻域都是高斯型且呈平方型指數(shù)衰減,在時、頻兩域均具有很好的局域性。運用小波變換可以準確地識別出突變點,運用高斯函數(shù)作為小波基的小波變換相當于對原來的信號做平滑處理,而運用高斯函數(shù)的一階導數(shù)作為小波基的小波變換可精確地找到信號所對應的突變點:若脈沖向上則表示信號增大,反之表示信號減弱(萬永革,2007)。如果使用傅立葉變換,就只能從頻率軸來分析該信號是由哪些頻率的波組成的,它只是將這些信息鋪開到整個頻率軸上。但是地震資料 f(t)組成的序列并不是平穩(wěn)的,而是隨時間不斷變化,傅立葉變換并不能提取 f(t)中的奇異性和突變點的信息。

    2 地震資料選取

    應變與構(gòu)造應力和地震活動是聯(lián)系緊密的物理量,在地震分期研究中較為常用(夏浩明,1987)。為了表征某個地區(qū)在不同時期的地震活動特征,根據(jù)該地區(qū)的應變積累和釋放變化,可將強震活動劃分為若干周期,每一周期包括一個應變能釋放高潮和低潮。一個地區(qū)的地震能量釋放是局部狀態(tài)的表征,由于地震輻射能量的平方根與地震的應變釋放成正比,在研究中常用應變釋放來表征地震能量釋放。根據(jù)歷史地震目錄可以推算出地震應變釋放的大小。

    本文采用貝尼奧夫曲線(Benioff,1951)來表征華北地區(qū)地震應變能的釋放。計算貝尼奧夫應變時采用時振梁等(1997)的公式:

    式中,εb(t)為以 t時刻為中心的 0.1 a之內(nèi)的貝尼奧夫應變積累,即某一地區(qū) 0.1 a內(nèi)所有大于M級地震對應變釋放的貢獻;M為震級。地震目錄資料的選取對研究結(jié)果有很大的影響。華北地區(qū)的地震資料在 1480年以后是比較完整的(黃瑋瓊等,1994),因此筆者結(jié)合《中國歷史強震目錄(公元前 23世紀—公元 1911年)》 (中國地震局震害防御司,1995)和中國地震地震臺網(wǎng)中心地震目錄選取了 1480~2009年華北地區(qū)(110°~120°E,34°~42°N)1 449個 M ≥4.0的地震進行研究(圖 1)。

    圖1 研究區(qū)內(nèi)M≥4.0地震震中分布圖Fig.1 Epicenter distribution ofM≥4.0 earthquakes in North China

    3 地震資料的小波分析

    貝尼奧夫應變曲線表述了一個地區(qū)地震能量積累的特性,通過研究貝尼奧夫應變曲線可以反映該區(qū)內(nèi)地震能量的釋放變化過程,但是用貝尼奧夫應變不能明顯地劃分出不同時間尺度和層次上地震活動的平靜期和活躍期。筆者分別對華北地區(qū) 1480~2009年M≥5.0的地震和 1970~2009年M≥4.0的地震進行計算,得到了兩組貝尼奧夫應變數(shù)據(jù),然后分別以高斯函數(shù)及其一階導數(shù)為小波基進行小波變換,分析不同尺度的平靜期和活躍期。

    圖2為研究區(qū)內(nèi) 1480~2009年M≥5.0地震的貝尼奧夫應變曲線及其相應的小波變換。圖2a最上端一條是貝尼奧夫應變?nèi)ゾ档脑假Y料曲線,下邊為不同時間尺度以高斯函數(shù)為小波基的小波變換曲線;圖 2b是不同時間尺度以高斯函數(shù)一階導數(shù)為小波基的小波變換,圖中數(shù)字是對應的小波時間尺度。從小波變換曲線能夠看出,總體上華北地區(qū)自 1480年以來M≥5.0地震大約呈 300 a的周期性變化,無論采用哪種小波尺度,這種周期性都較為明顯。在過去的 530多年時間內(nèi)華北地區(qū)大致經(jīng)歷了兩個活動期,前一活躍期高峰在1668年郯城 8.5級地震前后,后一活躍期高峰在1976年唐山7.8級地震前后。從小波變換的曲線形態(tài)來看,目前華北地區(qū)處于應變積累階段,且處于應變較低狀態(tài)。由于歷史地震資料主要都是來源于《中國歷史地震目錄(公元前 23世紀—公元 1911年)》,所以筆者對這些地震資料的分析結(jié)果跟張國民和張培震(2000)、馬宏生等(2005)結(jié)果一致,說明小波變換的分析結(jié)果是可信的。華北地區(qū)剛剛經(jīng)歷強震活躍期,不會有巨大地震或群集性中強地震發(fā)生。從以高斯一階導數(shù)為小波基的分析還能看出,自 1800年以后,當小波時間尺度小于 5 a時,有 4個明顯的突變點對應 4個活躍期,第一個活躍期對應 1830年6月 12日河北磁縣 7.5級地震,第二個對應 1888年6月 13日渤海灣 7.5級地震,第三個對應1937年8月 1日山東菏澤 7.0級地震,第四個對應 1966年3月 8日河北 7.2級地震和 1978年7月 28日唐山 7.8級地震,當小波時間尺度大于5 a時,這兩個地震對應的曲線突變點融合為一個突變點。

    圖2 研究區(qū)內(nèi) 1480~2009年M≥5.0地震的貝尼奧夫應變曲線的小波變換(a)以高斯函數(shù)為小波基的小波變換曲線;(b)以高斯函數(shù)一階導數(shù)為小波基的小波變換Fig.2 Results of seismic data ofM≥5.0 earthquakes in North China from 1480 to 2009 bywavelet analysis(a)The wavelet transform in various scales of Gauss function as base wavelet.The first upper trace is the original data of deaveraged Benioff strain in the period of 0.1 year in North China.The scales ofwavelet are labeled;(b)Wavelet transfo rm in various scales using derivative with Gauss function as base wavelet.The scales ofwavelet are labeled

    從小波變換曲線還可以看出:①隨著時間尺度的增加,極值點逐漸變少,且相鄰的極值點逐漸融合成較大的極值點;②各極值點的幅值反映了以同尺度高斯函數(shù)為小波基變化對應點局部梯度的大小,極值點的正負反映了濾波后貝尼奧夫曲線的邊沿方向特性——正極大值點對應上升沿,負極大值點對應下降沿;③極值點隨尺度變化的幅值大小反映了曲線邊沿變化特性——幅值較大反映濾波后貝尼奧夫曲線較劇烈的變化趨勢(較陡的邊沿),較小的幅值對應濾波后貝尼奧夫曲線較平緩的變化趨勢(較平的邊沿),由此可得到地震活躍期的開始年代(地震活躍突變點);④當尺度增大時極值點消失或成對結(jié)合成新的極值點,這反映了信號中隱含的某種 “嵌套結(jié)構(gòu)”。

    圖3 研究區(qū)內(nèi) 1970~2009年M≥4.0地震的貝尼奧夫應變曲線及其相應的小波變換(a)以高斯函數(shù)為小波基的小波變換曲線;(b)以高斯函數(shù)一階導數(shù)為小波基的小波變換Fig.3 Results of seismic data ofM≥4.0 earthquakes in North China from 1970 to 2009 bywavelet analysis(a)Wavelet transform in various scales of Gauss function as base wavelet. The first upper trace is the original data of deaveraged Benioff strain in the period of 0.1 year. The scales ofwavelet are labeled.(b)Wavelet transform in various scales using derivative with Guass function as base wavelet.The scales ofwavelet are labeled.

    運用小波變換能夠較好地分離和刻畫不同時間尺度的地震活動性,探測序列突變。地震活動性的強弱,是相對于不同時間尺度而言,因此,地震活動的活躍和平靜分期的問題只有建立在具體的時間尺度上才有意義。不同時間尺度上的地震活動,可能與地殼塊體的整體應力應變、不均勻性,以及與鄰區(qū)相互作用有密切的關系。筆者選用華北地區(qū) 1970~2009年的 M≥4.0地震,是因為這個時間段內(nèi)M≥4.0的地震已被非常準確、系統(tǒng)、完整地記錄下來,且M≥4.0地震對中強震的發(fā)生有很大的影響,選取這個震級更符合實際的地震能量的釋放演化過程。

    從圖 3可以看出,自 1970年以來,在0.5~5 a的時間尺度上華北地區(qū)M≥4.0地震一共有 5個活躍期,大致呈現(xiàn)出 7 a左右的周期,而且其強度隨時間在變化,隨著小波尺度的增加,地震活躍期突跳更加平滑。第一個活躍期對應的是1976年河北唐山 7.8級地震前后,第二個對應1983年山東菏澤 5.9級地震前后,第三個對應1989和 1991年山西省大同 6.1、5.8級地震前后,第四個對應 1998年河北張北 6.2級地震前后,第五個對應 2004年河北邢臺 5.0級和 2006年河北文安 5.1級地震前后,但是這個活躍期比較弱,當小波時間尺度增加到 4 a時已經(jīng)很不明顯。根據(jù)曲線分析,1998~2009年研究區(qū)內(nèi)缺少 4級以上地震,在 5 a以下的時間尺度上,研究區(qū)內(nèi)目前應該正在進入 4級以上地震的活躍期。

    從圖 3還可以看出,每個活躍期的出現(xiàn)總是具有雙峰特點的,當小波尺度小于 3 a時,這個特點非常明顯,1970年以來一直很有規(guī)律性,且以高斯一階導數(shù)函數(shù)為小波基時比以高斯函數(shù)為小波基時更加明顯,這也說明在分析突變時,小波基的選擇對結(jié)果有一定的影響,選擇合適的小波基函數(shù)可以更好地突出地震活動的規(guī)律。

    4 結(jié)論與討論

    (1)小波分析方法能得出不同時間尺度的地震活動特征,反映地震演化過程。運用該方法不僅可以將隱含在序列中的各種隨時間變化的周期顯現(xiàn)出來,還能確定地震活動突變點的位置,該位置對應于地震活躍期和平靜期的轉(zhuǎn)折處。

    (2)華北地震活動在不同的時間尺度上表現(xiàn)出不同的特征,華北地區(qū)M≥5.0地震活動出現(xiàn)兩次高峰期,周期約為 300 a,且剛剛經(jīng)歷了 1個活躍期,處于應力積累的平靜階段,該地區(qū)目前不會有巨大地震或群集性中強地震發(fā)生。1970年以來,研究區(qū)M≥4.0地震活動的周期大約為 7 a左右,而且其強度隨時間在變化,目前正在進入 4級以上地震的活躍期。因此,對于地震活躍或平靜期的討論,只有建立在具體的時間尺度上才有意義。

    對中國地震局地球物理研究所萬永革研究員和陳化然研究員給予的指導表示衷心的感謝!

    傅征祥,劉杰,劉桂萍,等 .2000.中國大陸成組強震活動水平的三分法及其應用[J].地震,20(增刊):34-38.

    黃瑋瓊,李文香,曹學鋒 .1994.中圍地震資料完整性研究之二[J].地震學報,16(4):423-431.

    劉虹,蘇有錦 .2007.云南地區(qū) 7級大震震源區(qū)地震活動特征研究[J].地震研究,30(1):18-23.

    劉希強,周蕙蘭,鄭治真,等 .1998.基于小波包變換的弱震相識別方法[J].地震學報,20(4):373-380.

    馬宗晉,傅征祥,張郢珍,等 .1982.1966~1976年中國九大地震[M].北京:地震出版社 .

    馬宏生,劉杰,張國民,等 .2002.用應變積累釋放模型研究中國大陸地塊分區(qū)地震活動[J].地震學報,24(6):569-578.

    馬宏生,劉杰,張曉東,等 .2005.基于 G—R關系的應變積累釋放模型研究中國大陸強震的分區(qū)活動[J].地震學報 27(4):355-366.

    蘇有錦,李忠華,劉祖萌,等 .2001.20世紀云南地區(qū)MS≥5.0級地震活動的基本特征[J].地震研究,24(1):1-9.

    邵輝成,杜興信,金學申,等 .2000.小波分析在地震趨勢預測中的應用[J].中國地震,16(1):48-52.

    時振梁,汪良謀,傅征祥,等 .1997.中國大陸中長期強震危險性預測方法研究[M].北京:海洋出版社 .

    吳忠良 .1999.鑒別地震活動周期性和隨機性的廣義 Rydelek Sacks檢驗[J].地震學報,21(5):487-494.

    萬永革,齊福榮,孟曉春,等 .2003.中國大陸及華北地區(qū)地震資料的小波分析[J].大地測量與地球動力學,23(4):28-33.

    萬永革 .2007.數(shù)字信號處理的 matlab實現(xiàn) [M].北京:科學出版社.

    夏浩明 .1987.大華北M≥4.0地震階級性分區(qū)特征及其意義[J].地球物理學報,30(3):281-291.

    徐義賢,王家映 .2000.基于連續(xù)小波變換的大地電磁信號譜估計方法[J].地球物理學報,43(5):677~683.

    張國民 .1987.我國大陸強震活動的韻律特征[J].地震地質(zhì),9(2):27-37.

    張國民,張培震 .2000.大陸強震機理與預測中期學術進展[J].中國基礎科學,(10):4-10.

    張培震,鄧起東,張國民,等 .2003.中國大陸的強震活動與活動地塊[J].中國科學 (D輯),33(增刊):12-20.

    中國地震局震害防御司 .1995.中國歷史強震目錄 (公元前 23世紀—公元 1911年)[M].北京:地震出版社 .

    Benioff H.1951.Crustal strain characteristics derived from earthquake sequences[J].TransAmer GeophyUnion,32(4):203-208.

    Logan B.1977.Information in the zero-crossirig of bandpass signals[J].Bell Syst Tech J,56,510-532.

    Mallat S.1991.Zero-crossing of a wavelet transform[J].IEEE Trans,IT-37(1):1019-1033.

    Mallat S,Hwang W L.1992.Singularity detection and processing with wavelets[J].IEEE TransAMS,IT-38:617-643.

    Research on the Characteristics of Seism ic Activity in North China UsingWavelet Transform

    ZHANG X iao-tao1,2,FENG X iang-dong3,ZHANG Yue-gang3,ZHANG Shuang-feng1,YANG Ya-qiong3,HAN L i-ping1
    (1.Central Seism ic station of Handan,Earthquake Adm inistration of Hebei Province,Handan056001,Hebei,China)
    (2.J ilin U niversity,Changchun130026,J ilin,China)
    (3.Earthquake Adm inistration of Hebei Province,Shijiazhuang050021,Hebei,China)

    Referring to the M≥5.0earthquakes from1500to2000and the M≥4.0earthquakes from1970to2009in N orth China,w e analyze the Benioff strain using w avelet transform in various scales.A nalysis show s that on a century-scale the M≥5.0earthquake-activity period has just finished and now the stress is accum ulating.A nd since1970there has existed activity cycle of about7years for theM≥4.0earthquakes; now theM≥4.0significant earthquake-activity period is beginning in N orth China.

    B enioff strain;w avelet transform; seism ic active period; seism ic quiet period;N orth China

    P315.5

    A

    1000-0666(2011)01-0024-06

    2010-03-09.

    河北省科技廳項目“地震預測與方法研究”(09276904D)資助 .

    猜你喜歡
    小波基華北地區(qū)時間尺度
    時間尺度上非完整系統(tǒng)的Noether準對稱性與守恒量
    時間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    力學學報(2021年10期)2021-12-02 02:32:04
    交直流混合微電網(wǎng)多時間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    華北地區(qū)SY1井鉆井技術難點及對策
    利用小波變換分析電能質(zhì)量擾動問題中的電壓驟升影響
    小波閾值圖像去噪中小波基選擇
    大連市暴雨多時間尺度研究分析
    華北地區(qū)不同林分類型枯落物層持水性能研究
    小波非參數(shù)回歸分析方法的實現(xiàn)及比較研究*
    2014年度華北地區(qū)經(jīng)營工作交流會在河北召開
    午夜日本视频在线| 国产 精品1| 国产熟女欧美一区二区| 韩国高清视频一区二区三区| 色吧在线观看| 男女啪啪激烈高潮av片| 亚洲国产精品999| 99热这里只有精品一区| 男女高潮啪啪啪动态图| 满18在线观看网站| 免费黄频网站在线观看国产| 亚洲成人手机| 亚洲欧美成人综合另类久久久| 青春草国产在线视频| 日本wwww免费看| 日韩熟女老妇一区二区性免费视频| 久久久久久久国产电影| 丝袜在线中文字幕| 国产精品一二三区在线看| 亚洲欧美日韩另类电影网站| 国产精品嫩草影院av在线观看| av国产久精品久网站免费入址| 18禁裸乳无遮挡动漫免费视频| 永久网站在线| 国产精品99久久久久久久久| 国产成人freesex在线| 人人妻人人爽人人添夜夜欢视频| 精品亚洲乱码少妇综合久久| 热re99久久国产66热| 久久精品久久久久久噜噜老黄| 在线观看免费视频网站a站| 91午夜精品亚洲一区二区三区| 亚洲国产精品999| 九九爱精品视频在线观看| 九色亚洲精品在线播放| 人人妻人人添人人爽欧美一区卜| 满18在线观看网站| 成人午夜精彩视频在线观看| 欧美日韩国产mv在线观看视频| 男女边吃奶边做爰视频| 男女边摸边吃奶| 免费高清在线观看日韩| 日日摸夜夜添夜夜添av毛片| 精品视频人人做人人爽| 亚洲精品中文字幕在线视频| 曰老女人黄片| 亚洲婷婷狠狠爱综合网| 国产精品免费大片| 看非洲黑人一级黄片| 久久国产精品大桥未久av| 97超碰精品成人国产| 啦啦啦啦在线视频资源| 免费日韩欧美在线观看| xxxhd国产人妻xxx| 99久久中文字幕三级久久日本| 边亲边吃奶的免费视频| 国产免费现黄频在线看| av专区在线播放| 毛片一级片免费看久久久久| 久久综合国产亚洲精品| 亚洲综合色惰| 欧美精品人与动牲交sv欧美| 人妻 亚洲 视频| 亚洲美女搞黄在线观看| 少妇熟女欧美另类| 美女内射精品一级片tv| 午夜免费鲁丝| 97超视频在线观看视频| 久久久久精品久久久久真实原创| 日韩人妻高清精品专区| 免费高清在线观看日韩| 一区二区日韩欧美中文字幕 | 久久久国产精品麻豆| 80岁老熟妇乱子伦牲交| 狠狠精品人妻久久久久久综合| 少妇人妻精品综合一区二区| 少妇猛男粗大的猛烈进出视频| 国产亚洲午夜精品一区二区久久| av福利片在线| 日本与韩国留学比较| 全区人妻精品视频| 日韩欧美精品免费久久| 岛国毛片在线播放| 亚洲av在线观看美女高潮| 少妇 在线观看| 大香蕉97超碰在线| 国语对白做爰xxxⅹ性视频网站| 好男人视频免费观看在线| 中文字幕精品免费在线观看视频 | 97超碰精品成人国产| 老熟女久久久| 国产高清三级在线| 国产av一区二区精品久久| 建设人人有责人人尽责人人享有的| 亚洲综合精品二区| 黑人欧美特级aaaaaa片| 欧美亚洲 丝袜 人妻 在线| 色网站视频免费| 女人久久www免费人成看片| 国产探花极品一区二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 青春草国产在线视频| 久久久久精品性色| 一本大道久久a久久精品| 一区二区三区免费毛片| av国产精品久久久久影院| 精品久久国产蜜桃| 一级片'在线观看视频| 色吧在线观看| 国产精品秋霞免费鲁丝片| 成人手机av| 91久久精品电影网| 美女大奶头黄色视频| 91久久精品电影网| 久久久国产精品麻豆| 国产精品一区二区在线观看99| 赤兔流量卡办理| 新久久久久国产一级毛片| 亚洲精品456在线播放app| 久久人人爽人人片av| 精品国产国语对白av| 亚洲精品国产色婷婷电影| 午夜福利在线观看免费完整高清在| 国产av一区二区精品久久| 18禁动态无遮挡网站| 精品人妻熟女av久视频| 一级毛片我不卡| 日本猛色少妇xxxxx猛交久久| 91aial.com中文字幕在线观看| 蜜臀久久99精品久久宅男| 91久久精品电影网| 欧美精品一区二区大全| 精品熟女少妇av免费看| 人人妻人人澡人人爽人人夜夜| 看免费成人av毛片| 人人妻人人爽人人添夜夜欢视频| 两个人免费观看高清视频| 国产淫语在线视频| 亚洲第一av免费看| av在线老鸭窝| 丁香六月天网| 午夜福利视频在线观看免费| 亚洲精品乱码久久久v下载方式| 狂野欧美激情性bbbbbb| 欧美日韩av久久| 精品人妻熟女av久视频| 老司机亚洲免费影院| 麻豆乱淫一区二区| 内地一区二区视频在线| 日韩精品免费视频一区二区三区 | 欧美bdsm另类| 自拍欧美九色日韩亚洲蝌蚪91| 看非洲黑人一级黄片| 精品国产露脸久久av麻豆| 国产成人freesex在线| 美女脱内裤让男人舔精品视频| 日本av免费视频播放| 日韩制服骚丝袜av| 国产精品三级大全| 最新中文字幕久久久久| 欧美日韩精品成人综合77777| 色婷婷av一区二区三区视频| 欧美亚洲 丝袜 人妻 在线| 99国产综合亚洲精品| 国产又色又爽无遮挡免| 夜夜爽夜夜爽视频| 精品一品国产午夜福利视频| 久久久国产欧美日韩av| 成年女人在线观看亚洲视频| 国产av一区二区精品久久| 欧美老熟妇乱子伦牲交| 精品久久久久久电影网| 99精国产麻豆久久婷婷| 人人妻人人澡人人爽人人夜夜| 日本免费在线观看一区| 两个人的视频大全免费| 黄片无遮挡物在线观看| 成人国产麻豆网| 一本—道久久a久久精品蜜桃钙片| 一级a做视频免费观看| 国产黄片视频在线免费观看| 久久99热这里只频精品6学生| 亚洲婷婷狠狠爱综合网| 国产视频内射| 丁香六月天网| 国产av一区二区精品久久| 十八禁网站网址无遮挡| 久久久久久久久久成人| 女人久久www免费人成看片| 色94色欧美一区二区| 五月开心婷婷网| 成人毛片60女人毛片免费| 亚洲精品视频女| 哪个播放器可以免费观看大片| 国产熟女午夜一区二区三区 | 国产一区亚洲一区在线观看| 七月丁香在线播放| 国国产精品蜜臀av免费| 伦理电影免费视频| 国精品久久久久久国模美| 夜夜骑夜夜射夜夜干| 久久久久精品性色| 99热6这里只有精品| 九九爱精品视频在线观看| 国产色爽女视频免费观看| 久久婷婷青草| 欧美激情 高清一区二区三区| 看免费成人av毛片| .国产精品久久| 成年美女黄网站色视频大全免费 | 久久 成人 亚洲| 国模一区二区三区四区视频| 成年人午夜在线观看视频| 黄色视频在线播放观看不卡| 超碰97精品在线观看| 丝瓜视频免费看黄片| 午夜影院在线不卡| 日韩欧美一区视频在线观看| 人妻 亚洲 视频| 伦理电影大哥的女人| 波野结衣二区三区在线| 国产精品 国内视频| 欧美精品人与动牲交sv欧美| 久久久久国产网址| 男女边摸边吃奶| 免费高清在线观看视频在线观看| 精品久久久久久电影网| 2018国产大陆天天弄谢| 一级二级三级毛片免费看| 99久国产av精品国产电影| 国产精品久久久久久精品古装| 建设人人有责人人尽责人人享有的| 亚洲不卡免费看| 色5月婷婷丁香| 国产精品国产三级专区第一集| 黄片无遮挡物在线观看| 99热网站在线观看| 精品久久国产蜜桃| 国产成人a∨麻豆精品| 一本一本综合久久| 91精品国产九色| 男人添女人高潮全过程视频| 国产精品.久久久| 亚洲激情五月婷婷啪啪| 日本vs欧美在线观看视频| 久久午夜综合久久蜜桃| 春色校园在线视频观看| 超碰97精品在线观看| 亚洲三级黄色毛片| 亚洲欧美精品自产自拍| 啦啦啦在线观看免费高清www| 女性被躁到高潮视频| 久久婷婷青草| 亚洲精品亚洲一区二区| 亚洲情色 制服丝袜| 亚洲av综合色区一区| 亚洲久久久国产精品| 亚洲精品日韩av片在线观看| 天堂8中文在线网| 一级黄片播放器| 另类亚洲欧美激情| 91国产中文字幕| 黑人欧美特级aaaaaa片| 国产精品久久久久久精品古装| 国产免费现黄频在线看| www.色视频.com| 日韩免费高清中文字幕av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男女无遮挡免费网站观看| 日韩不卡一区二区三区视频在线| 欧美xxxx性猛交bbbb| 黑丝袜美女国产一区| 欧美日韩视频高清一区二区三区二| 免费少妇av软件| 少妇丰满av| 三级国产精品片| 亚洲第一区二区三区不卡| 免费观看无遮挡的男女| 蜜桃在线观看..| 久久久久视频综合| 国产乱人偷精品视频| 国产精品成人在线| 久久人人爽人人片av| 满18在线观看网站| 国内精品宾馆在线| 最近中文字幕2019免费版| 国产综合精华液| 国产免费一级a男人的天堂| 丰满迷人的少妇在线观看| 日韩视频在线欧美| 国产成人91sexporn| 香蕉精品网在线| 精品人妻偷拍中文字幕| 五月天丁香电影| 国产高清不卡午夜福利| av在线播放精品| 免费看不卡的av| 久久人人爽人人爽人人片va| 丁香六月天网| av免费在线看不卡| 成年av动漫网址| 国产精品无大码| 少妇人妻精品综合一区二区| 两个人的视频大全免费| 久久久久网色| 免费观看在线日韩| 人妻人人澡人人爽人人| 精品国产露脸久久av麻豆| 免费观看无遮挡的男女| 精品人妻熟女毛片av久久网站| 97在线视频观看| 一级爰片在线观看| 国产精品国产三级国产av玫瑰| 蜜桃久久精品国产亚洲av| 男女啪啪激烈高潮av片| 成年人午夜在线观看视频| 九九在线视频观看精品| 能在线免费看毛片的网站| 岛国毛片在线播放| 国产成人免费无遮挡视频| 高清视频免费观看一区二区| 欧美少妇被猛烈插入视频| 久久久久久久亚洲中文字幕| 日韩制服骚丝袜av| 亚洲少妇的诱惑av| 日韩制服骚丝袜av| 高清毛片免费看| 精品久久久久久久久亚洲| 啦啦啦啦在线视频资源| 国产午夜精品久久久久久一区二区三区| 亚洲欧美一区二区三区国产| 高清午夜精品一区二区三区| 丝袜美足系列| 亚洲av二区三区四区| 18禁在线播放成人免费| 午夜福利,免费看| 亚洲欧洲精品一区二区精品久久久 | 国产片特级美女逼逼视频| 在线天堂最新版资源| 久久韩国三级中文字幕| 成人免费观看视频高清| 毛片一级片免费看久久久久| 欧美日本中文国产一区发布| 欧美日韩视频高清一区二区三区二| 日韩在线高清观看一区二区三区| 丰满乱子伦码专区| av视频免费观看在线观看| 爱豆传媒免费全集在线观看| 亚洲成人一二三区av| 精品熟女少妇av免费看| 人成视频在线观看免费观看| 亚洲av在线观看美女高潮| 色网站视频免费| 黑人高潮一二区| 91久久精品国产一区二区成人| 亚洲人成网站在线观看播放| 99久久中文字幕三级久久日本| 欧美日韩av久久| 久久av网站| 亚洲综合色惰| 下体分泌物呈黄色| 国产午夜精品一二区理论片| 中文天堂在线官网| 97在线视频观看| 秋霞在线观看毛片| 青春草视频在线免费观看| 亚洲精品,欧美精品| 黑人猛操日本美女一级片| 韩国av在线不卡| 熟女av电影| 亚洲av欧美aⅴ国产| 亚洲欧美一区二区三区国产| videossex国产| 精品亚洲乱码少妇综合久久| 少妇高潮的动态图| 亚洲综合色惰| 下体分泌物呈黄色| 看免费成人av毛片| 又大又黄又爽视频免费| 热99国产精品久久久久久7| 久久久久国产网址| 日韩中字成人| 少妇 在线观看| 久久人妻熟女aⅴ| 美女福利国产在线| 久久鲁丝午夜福利片| 黄片播放在线免费| 午夜免费男女啪啪视频观看| 日韩中文字幕视频在线看片| 婷婷色综合www| 十八禁高潮呻吟视频| 少妇丰满av| 欧美xxxx性猛交bbbb| 在线精品无人区一区二区三| 亚洲伊人久久精品综合| 2021少妇久久久久久久久久久| 欧美日本中文国产一区发布| 久久久国产一区二区| 高清毛片免费看| 大话2 男鬼变身卡| 丝袜美足系列| kizo精华| 最近最新中文字幕免费大全7| 麻豆精品久久久久久蜜桃| 欧美精品一区二区大全| 丝袜美足系列| 高清毛片免费看| 免费高清在线观看视频在线观看| 国产免费一区二区三区四区乱码| 午夜视频国产福利| 26uuu在线亚洲综合色| 99久久精品国产国产毛片| 亚洲第一av免费看| 欧美bdsm另类| 插阴视频在线观看视频| 看免费成人av毛片| 日韩强制内射视频| 欧美3d第一页| 高清黄色对白视频在线免费看| 黑人欧美特级aaaaaa片| 97在线人人人人妻| av网站免费在线观看视频| 在线观看美女被高潮喷水网站| 亚洲美女搞黄在线观看| 一级毛片黄色毛片免费观看视频| 免费观看无遮挡的男女| 在线观看一区二区三区激情| 黄片无遮挡物在线观看| 少妇人妻精品综合一区二区| 成年美女黄网站色视频大全免费 | 欧美亚洲日本最大视频资源| www.av在线官网国产| www.色视频.com| 曰老女人黄片| 亚洲经典国产精华液单| av又黄又爽大尺度在线免费看| 久久精品国产a三级三级三级| av免费观看日本| 国产精品久久久久成人av| 色5月婷婷丁香| 欧美少妇被猛烈插入视频| 热re99久久国产66热| av专区在线播放| 久久久久久久久大av| 午夜激情av网站| 插逼视频在线观看| 亚洲欧美日韩卡通动漫| 男人添女人高潮全过程视频| 天天躁夜夜躁狠狠久久av| 国产高清有码在线观看视频| 国产成人一区二区在线| 精品国产一区二区久久| 中文乱码字字幕精品一区二区三区| 国产精品人妻久久久久久| 看非洲黑人一级黄片| 国产在视频线精品| 精品少妇黑人巨大在线播放| 一级片'在线观看视频| 亚洲图色成人| 久久久午夜欧美精品| 日韩亚洲欧美综合| 极品少妇高潮喷水抽搐| 久久久午夜欧美精品| 晚上一个人看的免费电影| 日本wwww免费看| 精品久久久噜噜| 精品一区在线观看国产| 亚洲国产日韩一区二区| 中文字幕精品免费在线观看视频 | 大片电影免费在线观看免费| 国产成人一区二区在线| 能在线免费看毛片的网站| 中文天堂在线官网| 晚上一个人看的免费电影| 女人精品久久久久毛片| 两个人免费观看高清视频| 国产精品女同一区二区软件| 欧美最新免费一区二区三区| 精品人妻熟女av久视频| 男人操女人黄网站| 国产亚洲精品第一综合不卡 | 99热全是精品| 亚洲美女视频黄频| 在线观看免费视频网站a站| 国产成人av激情在线播放 | 国产黄片视频在线免费观看| 美女脱内裤让男人舔精品视频| 精品人妻熟女av久视频| 3wmmmm亚洲av在线观看| 欧美一级a爱片免费观看看| 伦理电影免费视频| 伦理电影大哥的女人| 老司机影院成人| 亚洲美女黄色视频免费看| 亚洲精品乱码久久久v下载方式| 欧美精品国产亚洲| 99re6热这里在线精品视频| 王馨瑶露胸无遮挡在线观看| 91午夜精品亚洲一区二区三区| 亚洲av不卡在线观看| 爱豆传媒免费全集在线观看| 色5月婷婷丁香| 18禁在线无遮挡免费观看视频| 91久久精品电影网| 亚洲国产欧美日韩在线播放| 国产精品秋霞免费鲁丝片| 美女中出高潮动态图| 国产精品欧美亚洲77777| 岛国毛片在线播放| 久久久久精品性色| 精品国产乱码久久久久久小说| 国产成人精品一,二区| 搡老乐熟女国产| videosex国产| 人妻少妇偷人精品九色| 99热国产这里只有精品6| 久久狼人影院| 国产黄片视频在线免费观看| 寂寞人妻少妇视频99o| 超碰97精品在线观看| 成人国产av品久久久| 精品一区在线观看国产| 久久精品人人爽人人爽视色| 亚洲人与动物交配视频| 久久久久精品久久久久真实原创| 日韩一区二区视频免费看| 最近中文字幕高清免费大全6| 亚洲欧美精品自产自拍| 丝袜在线中文字幕| 午夜激情久久久久久久| 国产成人freesex在线| 在现免费观看毛片| 欧美人与性动交α欧美精品济南到 | 精品一区二区免费观看| 考比视频在线观看| 久久午夜福利片| 18在线观看网站| 国产日韩欧美亚洲二区| 黑人欧美特级aaaaaa片| 亚洲第一区二区三区不卡| 婷婷成人精品国产| 一级片'在线观看视频| 九草在线视频观看| 国产成人aa在线观看| 草草在线视频免费看| 精品人妻在线不人妻| 精品少妇内射三级| 亚洲av中文av极速乱| 久久精品国产亚洲av天美| 伦理电影免费视频| 精品国产露脸久久av麻豆| 韩国高清视频一区二区三区| 爱豆传媒免费全集在线观看| 精品人妻熟女毛片av久久网站| 五月天丁香电影| 欧美bdsm另类| 午夜视频国产福利| 精品久久久精品久久久| 欧美亚洲 丝袜 人妻 在线| 午夜福利视频在线观看免费| 人人妻人人添人人爽欧美一区卜| 国产免费一级a男人的天堂| 亚洲精品国产av蜜桃| 国产亚洲午夜精品一区二区久久| 亚洲经典国产精华液单| 国产精品国产三级国产专区5o| 美女国产视频在线观看| 日日摸夜夜添夜夜添av毛片| 国内精品宾馆在线| 欧美bdsm另类| 18+在线观看网站| 夜夜骑夜夜射夜夜干| 黄片无遮挡物在线观看| 久久狼人影院| 亚洲高清免费不卡视频| 亚洲精品久久久久久婷婷小说| av在线app专区| 色视频在线一区二区三区| 国产女主播在线喷水免费视频网站| 91精品国产九色| 日韩av免费高清视频| 国产免费现黄频在线看| 青春草国产在线视频| 精品午夜福利在线看| 久久鲁丝午夜福利片| 最近2019中文字幕mv第一页| 少妇人妻久久综合中文| 欧美激情国产日韩精品一区| 亚洲国产精品成人久久小说| 久久精品国产自在天天线| 久热这里只有精品99| 性色avwww在线观看| 91aial.com中文字幕在线观看| 男女无遮挡免费网站观看| 美女内射精品一级片tv| 中文欧美无线码| h视频一区二区三区| 大片免费播放器 马上看| 免费高清在线观看视频在线观看| 欧美一级a爱片免费观看看| 久久韩国三级中文字幕| 亚洲精品日韩在线中文字幕| 一级爰片在线观看| 毛片一级片免费看久久久久| 国产精品成人在线| 黄色欧美视频在线观看| 好男人视频免费观看在线| 欧美日韩在线观看h| 男人添女人高潮全过程视频| 王馨瑶露胸无遮挡在线观看| 欧美日韩在线观看h| 中文字幕久久专区| 亚洲婷婷狠狠爱综合网| 亚洲精品日韩av片在线观看| 亚洲人成网站在线观看播放|