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

    2012年唐山ML5.2地震前后的速度結(jié)構(gòu)對(duì)比分析

    2019-05-04 11:23:56郭蕾宮猛王寧
    中國(guó)地震 2019年4期
    關(guān)鍵詞:雙差層析成像走時(shí)

    郭蕾 宮猛 王寧

    河北省地震局,石家莊 050021

    0 引言

    唐山地區(qū)位于陰山-燕山隆起帶東段的燕山隆起和華北平原拗陷帶北部的冀渤凹陷的結(jié)合部位,在燕山隆起一側(cè),區(qū)內(nèi)淺層活動(dòng)性斷裂發(fā)育,地質(zhì)構(gòu)造復(fù)雜,地震活動(dòng)頻繁(徐杰等,1996;羅艷等,2008)。1976年唐山7.8級(jí)大地震發(fā)生在唐山斷裂附近,震后由于老震區(qū)地殼介質(zhì)相對(duì)破碎,不易于積累更大能量,唐山地區(qū)地震序列強(qiáng)度和頻次總體隨時(shí)間呈現(xiàn)階段性衰減的趨勢(shì),近10年來僅發(fā)生一次ML≥5地震,即2012年5月28日唐山ML5.2地震。此次地震是在老震區(qū)5級(jí)地震平靜17年背景下發(fā)生的,且與地震序列強(qiáng)度衰減趨勢(shì)相比震級(jí)明顯偏高。研究表明地震孕育與發(fā)生是應(yīng)力積累與釋放的過程,在這個(gè)過程中通常會(huì)引起地殼介質(zhì)的物性改變,這種改變可通過地震波速的變化體現(xiàn)(Birch,1961)。因此,獲得該地震前后不同時(shí)間段的唐山老震區(qū)的三維速度結(jié)構(gòu),總結(jié)速度結(jié)構(gòu)的變化特征,對(duì)分析和探討該地震的形成、產(chǎn)生以及誘發(fā)因素具有意義。

    諸多學(xué)者對(duì)唐山地區(qū)的速度結(jié)構(gòu)進(jìn)行了研究,并取得大量成果。陳立華等(1990)根據(jù)構(gòu)造特征將華北地區(qū)劃分為4個(gè)區(qū)域,研究該區(qū)P波速度結(jié)構(gòu),結(jié)果表明橫向不均勻性和低速的存在可能是構(gòu)造活動(dòng)區(qū)的最主要特征;朱露培等(1990)利用空間頻率域離散參數(shù)化反演方法,計(jì)算京津唐張地區(qū)的地殼和上地幔深部結(jié)構(gòu),發(fā)現(xiàn)研究區(qū)內(nèi)大多數(shù)強(qiáng)震發(fā)生在低速區(qū)向高速區(qū)過渡的梯度帶;孫若昧等(1999)采用分類表格數(shù)據(jù)的AIC分析方法,對(duì)華北地區(qū)強(qiáng)震發(fā)生位置與速度結(jié)構(gòu)之間的關(guān)系進(jìn)行定量化評(píng)估;盧造勛等(2002)采用正交投影法重建了研究區(qū)(110°~126°E,36°~44°N)的地殼和上地幔的三維速度結(jié)構(gòu),結(jié)果顯示了“易震層”的速度分布特征,強(qiáng)震多發(fā)生在低速層上方的高速層內(nèi);王志鑠等(2008)利用LSQR方法反演了華北及鄰區(qū)地殼上地幔三維速度結(jié)構(gòu),發(fā)現(xiàn)唐山地區(qū)在25km、42km和60km深處存在連續(xù)的低速異常;賴曉玲等(2013)利用唐山臨時(shí)地震臺(tái)網(wǎng)數(shù)據(jù)5000多個(gè)P波到時(shí)進(jìn)行速度成像,結(jié)果表明唐山斷裂帶為速度陡變帶,其兩側(cè)的速度存在顯著差異,西北側(cè)為高速區(qū),東南側(cè)為低速區(qū);楊峰(2019)采用雙差層析成像的方法反演得到華北北部地區(qū)地殼三維P波速度結(jié)構(gòu),發(fā)現(xiàn)唐山震源區(qū)下方存在明顯的低速異常區(qū)。

    雖然前人在唐山地區(qū)速度結(jié)構(gòu)方面開展了大量研究,但是上述研究基本是靜態(tài)成像,未分析地震前后波速的動(dòng)態(tài)演化,然而正是這種動(dòng)態(tài)演化過程,可能有助于深入探索地震的發(fā)震機(jī)理(周龍泉等,2007;雷建設(shè)等,2010)。本研究采用結(jié)合波形互相關(guān)的雙差層析成像方法,以2年為時(shí)間間隔對(duì)唐山地區(qū)2010年5月~2016年5月的地震數(shù)據(jù)進(jìn)行反演,獲得2012年唐山ML5.2地震前后介質(zhì)的三維地震P波、S波速度的空間分布,結(jié)合此次地震的精定位震源參數(shù),探討唐山地區(qū)中強(qiáng)震發(fā)生前該區(qū)介質(zhì)的波速特征。

    1 理論方法與數(shù)據(jù)

    1.1 雙差層析成像方法

    雙差層析成像方法(Zhang et al,2003、2007)結(jié)合絕對(duì)走時(shí)和差分走時(shí)資料對(duì)三維速度結(jié)構(gòu)和震源參數(shù)的聯(lián)合反演。假設(shè)有2個(gè)相鄰地震的距離遠(yuǎn)小于其到同一臺(tái)站的距離,則可忽略速度差異,認(rèn)為它們到同1個(gè)臺(tái)站的路徑相似。對(duì)地震對(duì)的走時(shí)殘差做差即可得到雙差。

    根據(jù)射線理論,地震事件i到地震臺(tái)站k的體波到時(shí)可以表示為積分形式

    (1)

    式中,τi為地震事件i的發(fā)震時(shí)刻;u代表慢度;其中震源坐標(biāo)(x1,x2,x3)、發(fā)震時(shí)刻、射線路徑及慢度未知。將式(1)在震源處展開為一階Taylor級(jí)數(shù)得到線性方程

    (2)

    (3)

    式(3)也可表示為

    (4)

    雙差層析成像利用絕對(duì)走時(shí)數(shù)據(jù)確定震源區(qū)外的速度結(jié)構(gòu),利用到時(shí)差數(shù)據(jù)確定震源區(qū)的精細(xì)速度結(jié)構(gòu)。其通過速度結(jié)構(gòu)的適當(dāng)調(diào)整提高地震定位的精度,而地震定位精度的提高又可以降低由于射線路徑差異造成的速度結(jié)構(gòu)的誤差,進(jìn)而提高速度結(jié)構(gòu)反演的精度(呂子強(qiáng)等,2016)。因此,雙差層析成像方法可以得到比標(biāo)準(zhǔn)層析成像更精確的定位和速度結(jié)構(gòu)結(jié)果。

    1.2 地震數(shù)據(jù)與模型

    唐山老震區(qū)小震序列豐富,且臺(tái)站分布較為合理,能有效地監(jiān)測(cè)該區(qū)的地震活動(dòng),為雙差層析成像提供了良好的數(shù)據(jù)基礎(chǔ)。本文研究區(qū)域?yàn)?9°00′~40°30′N、117°30′~119°30′E,收集河北地震臺(tái)網(wǎng)2010年5月~2016年5月唐山地區(qū)的地震事件波形及震相報(bào)告,根據(jù)觀測(cè)數(shù)據(jù),2010年5月~2012年5月共1398個(gè)地震事件,P波走時(shí)16822條,S波走時(shí)17616條;2012年5月~2014年5月共1628個(gè)地震事件,P波走時(shí)18958條,S波走時(shí)20546條;2014年5月~2016年5月共3219個(gè)地震事件,P波走時(shí)22937條,S波走時(shí)23697條。2014年5月~2016年5月地震事件的數(shù)量與前2個(gè)時(shí)段有明顯差別,這可能是因?yàn)榻陙淼卣鹉夸浘忍岣?,存在一定?shù)量的單臺(tái)地震,為保證對(duì)比效果,需對(duì)數(shù)據(jù)進(jìn)行篩選。

    圖 1 唐山地區(qū)2010年5月~2016年5月的地震震相走時(shí)曲線紅色為P波走時(shí)曲線;綠色為S波走時(shí)曲線

    結(jié)合走時(shí)曲線(圖1)剔除誤差較大的震相數(shù)據(jù),且挑選出滿足每個(gè)地震需要的最少記錄震相數(shù)(8個(gè)),去掉單臺(tái)地震,篩選后,2010年5月~2012年5月可用的地震事件1271個(gè),P波走時(shí)16280條,S波走時(shí)16598條;2012年5月~2014年5月可用的地震事件1423個(gè),P波走時(shí)18274條,S波走時(shí)19396條;2014年5月~2016年5月可用的地震事件1520個(gè),P波走時(shí)19246條,S波走時(shí)20001條。3個(gè)時(shí)段臺(tái)站分布、地震數(shù)量及到時(shí)資料數(shù)基本相當(dāng)。

    在雙差層析成像中,事件對(duì)之間的距離對(duì)反演結(jié)果至關(guān)重要。該閾值越小則反演結(jié)果精度越高,但取值越小則地震對(duì)的數(shù)量越少,不利于反演運(yùn)算,因此確定合適的事件對(duì)距離是雙差層析成像計(jì)算過程中重要的一環(huán)。通過統(tǒng)計(jì)研究區(qū)所有臺(tái)站記錄到的事件對(duì)的震源距和相應(yīng)的波形互相關(guān)系數(shù)(圖2),確定事件對(duì)距離為30km。經(jīng)過篩選后的地震震中主要集中在唐山斷裂附近,震源深度多數(shù)小于30km,為確保每個(gè)網(wǎng)格有足夠的射線通過,將網(wǎng)格間距設(shè)置為0.2°(圖3)。根據(jù)于湘?zhèn)サ?2010)的一維速度模型研究成果(表1),垂直向下的深度分別取0、5km、10km、15km、20km、25km和30km,與速度模型的分層一致。

    圖 2 所有臺(tái)站地震對(duì)和互相關(guān)系數(shù)統(tǒng)計(jì)(a)、(b)為P波數(shù)據(jù);(c)、(d)為S波數(shù)據(jù)

    圖 3 研究區(qū)雙差層析成像地震事件震中分布、深度剖面和統(tǒng)計(jì)圖

    2 結(jié)果與討論

    結(jié)合初始一維速度模型,采用帶阻尼的最小二乘正交分解法,通過8組23次迭代獲得唐山地區(qū)的三維速度結(jié)構(gòu),采用DWS(Derivative Weighted Sum)對(duì)反演結(jié)果的可靠性進(jìn)行檢驗(yàn)。

    表 1 一維分層速度模型

    2.1 結(jié)果可靠性分析

    DWS能很好地衡量節(jié)點(diǎn)的射線密度(Thurber et al,1999),前人的研究認(rèn)為在DWS大于100的區(qū)域具有很高的模型分辨率(Scarfì et al,2007)。本文采用DWS檢測(cè)3個(gè)時(shí)段反演速度模型的分辨率。根據(jù)反演結(jié)果得到了P波DWS分布圖(圖4~6),結(jié)果表明,反演網(wǎng)格大部分區(qū)域的DWS大于100,具有較好的模型分辨率,3個(gè)時(shí)段的震源深度均在5~20km,尤其在唐山斷裂附近,具有很高的射線密度分布,能保證得到較好的分辨率。

    圖 4 2010年5月~2012年5月研究區(qū)P波DWS分布

    圖 5 2012年5月~2014年5月研究區(qū)P波DWS分布

    圖 6 2014年5月~2016年5月研究區(qū)P波DWS分布

    表 2 精定位地震基本參數(shù)

    2.2 2012年唐山ML5.2地震前后三維速度結(jié)構(gòu)對(duì)比分析

    精定位后2012年唐山ML5.2地震的基本參數(shù)如表2 所示。精定位后該地震的震源深度為11km,且震中位于唐山斷裂附近。DWS結(jié)果表明在5~20km的震源深度,尤其在唐山斷裂附近具有較高的模型分辨率,且前后3個(gè)時(shí)間段的地震分布基本一致,獲得的射線密度分布差異不大,反演過程從數(shù)據(jù)處理到反演閾值的選取均保持一致,因此認(rèn)為2012年唐山ML5.2地震前后三維速度結(jié)構(gòu)的對(duì)比結(jié)果具有一定的參考意義。選取前后3個(gè)時(shí)期10km的水平速度結(jié)構(gòu)切片進(jìn)行觀察比較。

    由圖7 可以看出,在10km 深度,研究區(qū)P波平均速度結(jié)構(gòu)橫向不均性顯著,存在2個(gè)NNW向的高低速交界區(qū),分別與灤縣-樂亭斷裂和薊運(yùn)河斷裂位置對(duì)應(yīng)。沿NE向的唐山斷裂帶具有串珠狀高低速的線性特征,反映了唐山斷裂帶是研究區(qū)內(nèi)的主要構(gòu)造帶,這與賴曉玲等(2013)的研究成果基本一致。2012唐山ML5.2地震前,震中位于在P波高低速過渡帶偏向高速的位置(圖7),出現(xiàn)這種現(xiàn)象的一種解釋是,低速體難以積累能量,而高速體有可能積累足以引發(fā)大地震的應(yīng)變能(孫若昧等,1995)。震后震中附近的波速有一定程度的降低,西北方向高速塊體范圍減小;而震中東南方向的低速區(qū)變?yōu)楦咚賲^(qū),形成了1個(gè)高速塊體(圖8);之后,震中附近高低速過渡帶逐漸消失。

    圖 7 2010年5月~2012年5月(a)、2012年5月~2014年5月(b)及2014年5月~2016年5月(c)的P波速度結(jié)構(gòu)白色實(shí)線為DWS=100

    圖 8 2012年5月~2014年5月與2010年5月~2012年5月(a)、2014年5月~2016年5月與2012年5月~2014年5月(b)的P波速度變化白色實(shí)線為DWS=100

    圖 9 剖面位置示意圖

    2012唐山ML5.2地震發(fā)生在唐山斷裂附近,根據(jù)震源機(jī)制計(jì)算結(jié)果發(fā)現(xiàn),有1個(gè)節(jié)面與唐山斷裂走向一致,因此沿唐山斷裂方向(AA′)和垂直唐山斷裂方向(BB′)做剖面(圖9),觀察該地震前后在深度上的速度結(jié)構(gòu)變化。由圖10、11 可以看出,5~15km深度為P波低速區(qū)與高速區(qū)對(duì)比強(qiáng)烈的交界地帶,唐山地區(qū)中下地殼存在低速異常,這與王椿鏞等(2017)的研究成果一致。

    王椿鏞等(2017)的地震層析成像結(jié)果表明,華北大震大多發(fā)生在高速和低速區(qū)的過渡帶上。水平速度結(jié)構(gòu)切片與垂直剖面速度結(jié)構(gòu)切片均顯示,此次地震發(fā)生前,震中位置處于高低速過渡帶,震后震中附近的波速出現(xiàn)一定程度的降低,與王椿鏞等(2017)的研究成果保持一致。2012年唐山ML5.2地震前,震中位置處在高低速過渡帶附近,震后震中附近的波速出現(xiàn)的小范圍的降低,震中下方的低速區(qū)范圍逐漸增大,震中附近高低速過渡帶逐漸消失。

    圖 10 沿AA′方向3個(gè)時(shí)段研究區(qū)P波和S波速度剖面

    圖 11 沿BB′方向3個(gè)時(shí)段研究區(qū)P波和S波速度剖面

    圖 12 2010年5月~2012年5月(a)、2012年5月~2014年5月(b)及2014年5月~2016年5月(c)的波速比白色實(shí)線為DWS=100

    波速比能夠綜合反映P波和S波變化的相對(duì)大小,單純地反映介質(zhì)的泊松比,所受影響比P波、S波速度小很多。圖12 給出了10km 深度的研究區(qū)波速比分布,唐山ML5.2地震前震源區(qū)附近存在明顯的波速比低值,震后波速比值回升。王亞茹等(2016)研究結(jié)果表明,2012年唐山ML5.2地震前,震源區(qū)附近波速比出現(xiàn)了下降—低值的過程。由于本研究在震前僅計(jì)算1個(gè)時(shí)段的數(shù)據(jù),無法判斷震前是否存在下降—低值的過程,但是低值—恢復(fù)的過程與已有研究結(jié)果(王亞茹等,2016)一致。

    3 結(jié)論

    采用波形互相關(guān)的雙差層析成像方法對(duì)唐山地區(qū)2010年5月~2016年5月的地震數(shù)據(jù)以2年為時(shí)間間隔進(jìn)行反演,獲得2012年唐山ML5.2地震前后介質(zhì)的三維地震P波、S波速度的空間分布。結(jié)果表明:

    (1)在10km 深度,研究區(qū)P波平均速度結(jié)構(gòu)橫向不均性明顯,沿NE向的唐山斷裂帶有串珠狀高低速的線性特征,反映了唐山斷裂帶為研究區(qū)內(nèi)的主要構(gòu)造帶。

    (2)P波和S波水平速度結(jié)構(gòu)切片與垂直剖面速度結(jié)構(gòu)切片均表明,2012唐山ML5.2地震前,震中位置處于高低速過渡帶,震后震中附近的波速出現(xiàn)一定程度的降低,原有的高低速過渡帶逐漸消失。

    (3)唐山ML5.2地震前,震源區(qū)波速比存在低值異常,震后波速比值回升。

    地震的形成是一個(gè)復(fù)雜的過程,用層析成像的方法探討2012唐山ML5.2地震前后波速演化過程時(shí),雖然盡量保持不同時(shí)段地震數(shù)量、臺(tái)站分布及反演參數(shù)的一致性,但是地震射線覆蓋的差別,仍會(huì)對(duì)成像結(jié)果造成一定影響。因而,本研究中的速度變化是否僅僅受此次地震的影響而產(chǎn)生,仍需要進(jìn)一步研究。

    致謝:中國(guó)科學(xué)技術(shù)大學(xué)張海江教授提供了雙差層析成像程序,審稿專家對(duì)本文提出了寶貴的修改意見,在此一并感謝。

    猜你喜歡
    雙差層析成像走時(shí)
    虛擬地震臺(tái)陣雙差測(cè)深法及應(yīng)用
    BDS中長(zhǎng)基線三頻RTK算法研究
    基于大數(shù)據(jù)量的初至層析成像算法優(yōu)化
    BDS參考站間低高度角衛(wèi)星三頻整周模糊度解算方法
    基于快速行進(jìn)法地震層析成像研究
    來了晃一圈,走時(shí)已鍍金 有些掛職干部“假裝在基層”
    基于雙差的波動(dòng)方程反射波旅行時(shí)反演方法
    基于分布式無線網(wǎng)絡(luò)的無線電層析成像方法與實(shí)驗(yàn)研究
    基于多級(jí)小波域變換的時(shí)域擴(kuò)散熒光層析成像方法
    成人三级黄色视频| 偷拍熟女少妇极品色| 精品欧美国产一区二区三| 国内久久婷婷六月综合欲色啪| 在线观看舔阴道视频| 精品一区二区三区视频在线观看免费| 亚洲av美国av| 日韩欧美精品免费久久| 国产在视频线在精品| 亚洲,欧美,日韩| 国内少妇人妻偷人精品xxx网站| 黄色女人牲交| 国产精品久久久久久av不卡| 一本精品99久久精品77| 人妻制服诱惑在线中文字幕| 少妇被粗大猛烈的视频| 美女高潮喷水抽搐中文字幕| 亚洲最大成人av| 一区二区三区四区激情视频 | 免费一级毛片在线播放高清视频| 淫妇啪啪啪对白视频| 人妻丰满熟妇av一区二区三区| 精品久久久久久久久久免费视频| 观看免费一级毛片| 别揉我奶头 嗯啊视频| 欧美色欧美亚洲另类二区| 在线看三级毛片| 一级av片app| 99热这里只有精品一区| 亚洲av成人精品一区久久| 久久久精品大字幕| 免费av观看视频| 成人永久免费在线观看视频| 成人国产综合亚洲| 1000部很黄的大片| 亚洲电影在线观看av| 久久精品国产亚洲av天美| 精品久久久久久成人av| 中文字幕av在线有码专区| 成人特级黄色片久久久久久久| 久久精品国产自在天天线| 国产精品爽爽va在线观看网站| 少妇猛男粗大的猛烈进出视频 | 国产熟女欧美一区二区| 久久精品人妻少妇| 亚洲av免费在线观看| 最近视频中文字幕2019在线8| av视频在线观看入口| 久久久久九九精品影院| 亚洲一区高清亚洲精品| 国产中年淑女户外野战色| 成人亚洲精品av一区二区| 啦啦啦啦在线视频资源| 国产精品一区二区免费欧美| 亚洲午夜理论影院| 久久精品国产亚洲av涩爱 | 一本一本综合久久| 性色avwww在线观看| 99在线视频只有这里精品首页| 午夜福利在线观看吧| 一级黄片播放器| 蜜桃久久精品国产亚洲av| 免费一级毛片在线播放高清视频| 欧美激情久久久久久爽电影| 欧美日韩综合久久久久久 | 亚洲狠狠婷婷综合久久图片| 99九九线精品视频在线观看视频| 一个人观看的视频www高清免费观看| 俺也久久电影网| 成人特级av手机在线观看| 欧美高清性xxxxhd video| 国产又黄又爽又无遮挡在线| 亚洲精华国产精华液的使用体验 | 亚洲av中文av极速乱 | 日本在线视频免费播放| 午夜视频国产福利| 日本精品一区二区三区蜜桃| 亚洲熟妇熟女久久| 男女之事视频高清在线观看| 观看美女的网站| 精品国产三级普通话版| or卡值多少钱| 亚洲最大成人手机在线| 亚洲不卡免费看| 91在线观看av| 99久久久亚洲精品蜜臀av| 男人舔女人下体高潮全视频| 波多野结衣高清作品| 久久久久九九精品影院| 亚洲天堂国产精品一区在线| 欧美激情久久久久久爽电影| 午夜福利在线观看免费完整高清在 | 网址你懂的国产日韩在线| 中文资源天堂在线| 国产极品精品免费视频能看的| 97人妻精品一区二区三区麻豆| 久久久久久大精品| 国产欧美日韩一区二区精品| 国产色爽女视频免费观看| 网址你懂的国产日韩在线| 无遮挡黄片免费观看| 久久精品影院6| a在线观看视频网站| 亚洲美女搞黄在线观看 | 久久香蕉精品热| 亚洲精品成人久久久久久| 国产一区二区在线观看日韩| 欧美精品国产亚洲| 欧美zozozo另类| 国产精品美女特级片免费视频播放器| 欧美3d第一页| 日日摸夜夜添夜夜添小说| 国产乱人伦免费视频| 久久欧美精品欧美久久欧美| 久久久久久久久大av| 国产精品日韩av在线免费观看| 亚洲av熟女| 嫁个100分男人电影在线观看| 久久久久性生活片| 精品日产1卡2卡| 中文字幕免费在线视频6| 中文在线观看免费www的网站| 色综合婷婷激情| 久久久久久久久久成人| 欧美日本视频| 看免费成人av毛片| 在线观看av片永久免费下载| 国产主播在线观看一区二区| 国产一区二区亚洲精品在线观看| 天天躁日日操中文字幕| x7x7x7水蜜桃| 午夜久久久久精精品| 免费看光身美女| 国产精品综合久久久久久久免费| 精品一区二区三区视频在线观看免费| 熟女电影av网| 少妇人妻一区二区三区视频| 国产一区二区在线av高清观看| 久久久久国内视频| 一卡2卡三卡四卡精品乱码亚洲| 日本在线视频免费播放| 免费观看精品视频网站| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美国产在线观看| 搞女人的毛片| 亚洲四区av| 小说图片视频综合网站| 日本三级黄在线观看| 国产久久久一区二区三区| 中国美女看黄片| 国产91精品成人一区二区三区| 日韩人妻高清精品专区| 中文字幕高清在线视频| 国产精品久久久久久久久免| 午夜免费激情av| 国产精品日韩av在线免费观看| 午夜福利在线在线| eeuss影院久久| 九九在线视频观看精品| 久久久精品大字幕| 99久国产av精品| 久久亚洲真实| 亚洲av美国av| 在线观看一区二区三区| 亚洲专区国产一区二区| 窝窝影院91人妻| 十八禁国产超污无遮挡网站| 成人av在线播放网站| 欧美性感艳星| 哪里可以看免费的av片| 久久人人精品亚洲av| 日韩精品中文字幕看吧| 午夜视频国产福利| 女的被弄到高潮叫床怎么办 | 午夜激情福利司机影院| 嫩草影院新地址| 搡老妇女老女人老熟妇| 欧美成人一区二区免费高清观看| bbb黄色大片| 亚洲欧美精品综合久久99| 免费观看精品视频网站| 国产精品1区2区在线观看.| 亚洲 国产 在线| 国产不卡一卡二| 日本免费一区二区三区高清不卡| 18+在线观看网站| 国产成人a区在线观看| 国产视频内射| 亚洲va在线va天堂va国产| a级毛片a级免费在线| 亚洲精品成人久久久久久| 黄色日韩在线| 亚洲无线观看免费| 久久久久久久久久久丰满 | 色尼玛亚洲综合影院| 日韩亚洲欧美综合| 制服丝袜大香蕉在线| 国产私拍福利视频在线观看| av天堂在线播放| av天堂中文字幕网| 91狼人影院| 亚洲欧美日韩高清在线视频| 很黄的视频免费| 97热精品久久久久久| 制服丝袜大香蕉在线| 97超视频在线观看视频| 亚洲精品在线观看二区| 久久久久久九九精品二区国产| 日韩欧美精品v在线| 舔av片在线| 国产视频内射| 精品午夜福利视频在线观看一区| 精品人妻熟女av久视频| 男人的好看免费观看在线视频| 1024手机看黄色片| 精品久久久久久久久久免费视频| 久久香蕉精品热| 黄色一级大片看看| 九九在线视频观看精品| 国产毛片a区久久久久| 中国美白少妇内射xxxbb| 看十八女毛片水多多多| 国产精品99久久久久久久久| 午夜免费男女啪啪视频观看 | 精品人妻一区二区三区麻豆 | 99国产精品一区二区蜜桃av| 久久久久久久午夜电影| 久久久成人免费电影| 中文字幕熟女人妻在线| 国产精品三级大全| 九色成人免费人妻av| 哪里可以看免费的av片| 亚洲精品色激情综合| 淫秽高清视频在线观看| 精品午夜福利在线看| 亚洲自偷自拍三级| 欧美一区二区国产精品久久精品| 国产av不卡久久| 国产国拍精品亚洲av在线观看| 国国产精品蜜臀av免费| 波多野结衣高清作品| 丝袜美腿在线中文| 亚洲欧美日韩卡通动漫| 久久久久久久精品吃奶| 男女视频在线观看网站免费| 91久久精品国产一区二区成人| 亚洲欧美清纯卡通| 午夜福利在线观看吧| 性色avwww在线观看| 久久精品国产亚洲网站| 国产中年淑女户外野战色| 日本一本二区三区精品| 国产又黄又爽又无遮挡在线| 亚洲中文字幕日韩| 国产免费av片在线观看野外av| 亚洲av一区综合| 亚洲成a人片在线一区二区| 啦啦啦韩国在线观看视频| av福利片在线观看| 午夜影院日韩av| 亚洲av美国av| 毛片一级片免费看久久久久 | 午夜精品一区二区三区免费看| 国产成人一区二区在线| 国内少妇人妻偷人精品xxx网站| 免费av毛片视频| 超碰av人人做人人爽久久| 国产亚洲av嫩草精品影院| 真实男女啪啪啪动态图| 国产精品久久久久久久电影| 日本精品一区二区三区蜜桃| 日韩av在线大香蕉| 午夜精品在线福利| 国产v大片淫在线免费观看| 18+在线观看网站| 日本黄大片高清| 精品人妻视频免费看| 神马国产精品三级电影在线观看| 精品人妻1区二区| 在线免费观看的www视频| 岛国在线免费视频观看| 免费在线观看成人毛片| 窝窝影院91人妻| 小说图片视频综合网站| 精品国产三级普通话版| 九色成人免费人妻av| 校园人妻丝袜中文字幕| 嫁个100分男人电影在线观看| 两人在一起打扑克的视频| 黄色丝袜av网址大全| 国产亚洲精品av在线| 欧美区成人在线视频| 深夜a级毛片| 床上黄色一级片| 日韩av在线大香蕉| 91精品国产九色| 赤兔流量卡办理| 一区二区三区激情视频| 午夜a级毛片| 亚洲av一区综合| 桃色一区二区三区在线观看| 赤兔流量卡办理| 欧美一区二区精品小视频在线| 日日撸夜夜添| 一区二区三区免费毛片| 国产日本99.免费观看| 精品久久久久久久久久免费视频| 亚洲经典国产精华液单| 欧美bdsm另类| 日本爱情动作片www.在线观看 | 亚洲狠狠婷婷综合久久图片| 久久久久国内视频| 内地一区二区视频在线| 日韩,欧美,国产一区二区三区 | av天堂中文字幕网| 国产不卡一卡二| 欧美zozozo另类| 极品教师在线视频| 日韩,欧美,国产一区二区三区 | 日日啪夜夜撸| 九九热线精品视视频播放| 别揉我奶头~嗯~啊~动态视频| 中文字幕免费在线视频6| 国产精品99久久久久久久久| 成人三级黄色视频| 日韩欧美免费精品| 欧美3d第一页| 国产探花极品一区二区| 日韩欧美 国产精品| 日本 av在线| 久久久久久国产a免费观看| 亚洲人成伊人成综合网2020| 成熟少妇高潮喷水视频| 不卡视频在线观看欧美| 国产三级在线视频| 丝袜美腿在线中文| 校园春色视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 色在线成人网| 欧美精品国产亚洲| 最近最新免费中文字幕在线| 欧美高清性xxxxhd video| 国产黄片美女视频| 国产精品av视频在线免费观看| 色播亚洲综合网| 欧美成人免费av一区二区三区| 国产成人一区二区在线| 国产精品无大码| 欧美激情在线99| 九九爱精品视频在线观看| 一级黄色大片毛片| 97超视频在线观看视频| 久久久久国内视频| 久久久精品大字幕| 欧美一区二区精品小视频在线| 欧美一区二区国产精品久久精品| 狂野欧美白嫩少妇大欣赏| 中文资源天堂在线| 日韩中字成人| 欧美激情在线99| 真实男女啪啪啪动态图| 亚洲国产精品合色在线| 精品一区二区三区视频在线观看免费| 精品国产三级普通话版| 九九热线精品视视频播放| 亚洲人成伊人成综合网2020| 久久久久久久久久久丰满 | 国产精品国产高清国产av| 人人妻人人澡欧美一区二区| 亚洲av.av天堂| 久久久久久久精品吃奶| 少妇人妻精品综合一区二区 | 伦理电影大哥的女人| 国产高清视频在线播放一区| 国产私拍福利视频在线观看| 男女之事视频高清在线观看| 久久99热6这里只有精品| 精品一区二区三区av网在线观看| 国产精品乱码一区二三区的特点| 日韩欧美 国产精品| 丰满人妻一区二区三区视频av| 搡老妇女老女人老熟妇| 国产69精品久久久久777片| 成人美女网站在线观看视频| 久久香蕉精品热| 国产不卡一卡二| 成人午夜高清在线视频| 欧美xxxx性猛交bbbb| 欧美绝顶高潮抽搐喷水| 日韩欧美免费精品| 午夜精品久久久久久毛片777| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区免费观看 | 一级av片app| 麻豆一二三区av精品| 老熟妇仑乱视频hdxx| www.色视频.com| 一区二区三区免费毛片| 国产精品电影一区二区三区| 男人舔女人下体高潮全视频| 可以在线观看的亚洲视频| 黄色日韩在线| 美女免费视频网站| 男女之事视频高清在线观看| 日本撒尿小便嘘嘘汇集6| 成人精品一区二区免费| 精品一区二区三区视频在线观看免费| 久久99热这里只有精品18| 日本免费a在线| 女人被狂操c到高潮| 欧美高清成人免费视频www| 日韩欧美 国产精品| 美女高潮喷水抽搐中文字幕| www.www免费av| 日本在线视频免费播放| 欧美不卡视频在线免费观看| 97人妻精品一区二区三区麻豆| 亚洲最大成人av| 简卡轻食公司| 可以在线观看毛片的网站| 身体一侧抽搐| 欧美黑人欧美精品刺激| 俺也久久电影网| 国产综合懂色| 久久精品91蜜桃| 国产精品免费一区二区三区在线| 一本久久中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 在线观看66精品国产| 狂野欧美激情性xxxx在线观看| 欧美最黄视频在线播放免费| 日本免费一区二区三区高清不卡| 亚洲七黄色美女视频| aaaaa片日本免费| 亚洲性久久影院| 日韩人妻高清精品专区| 亚洲最大成人av| 韩国av在线不卡| 在线观看舔阴道视频| 2021天堂中文幕一二区在线观| 亚洲狠狠婷婷综合久久图片| 日韩大尺度精品在线看网址| 国产一级毛片七仙女欲春2| 91在线精品国自产拍蜜月| 欧美最黄视频在线播放免费| 1000部很黄的大片| 国产成年人精品一区二区| 午夜免费激情av| 十八禁国产超污无遮挡网站| 日本黄大片高清| 久久人人爽人人爽人人片va| 91精品国产九色| 少妇熟女aⅴ在线视频| www.色视频.com| 免费人成在线观看视频色| 成人av一区二区三区在线看| 美女高潮喷水抽搐中文字幕| 亚洲国产色片| 日日摸夜夜添夜夜添小说| 大型黄色视频在线免费观看| 久久久久久伊人网av| 欧美成人a在线观看| 久久久久久久久久久丰满 | 日本五十路高清| 一级毛片久久久久久久久女| 女人被狂操c到高潮| 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 日韩 亚洲 欧美在线| 亚洲欧美日韩东京热| 国产高清激情床上av| 成人无遮挡网站| 日本一本二区三区精品| 22中文网久久字幕| av视频在线观看入口| 国产在视频线在精品| 村上凉子中文字幕在线| 可以在线观看的亚洲视频| 波野结衣二区三区在线| АⅤ资源中文在线天堂| 日本在线视频免费播放| 午夜影院日韩av| 日韩亚洲欧美综合| 看十八女毛片水多多多| 国产精品精品国产色婷婷| 又黄又爽又刺激的免费视频.| 五月伊人婷婷丁香| 国产中年淑女户外野战色| 亚洲电影在线观看av| 一区二区三区四区激情视频 | 色哟哟·www| 国产又黄又爽又无遮挡在线| 少妇丰满av| 热99re8久久精品国产| av.在线天堂| 嫩草影院入口| 天天一区二区日本电影三级| 免费看美女性在线毛片视频| 亚洲午夜理论影院| 麻豆国产97在线/欧美| 淫妇啪啪啪对白视频| 久久亚洲精品不卡| 中出人妻视频一区二区| 亚洲真实伦在线观看| bbb黄色大片| 12—13女人毛片做爰片一| videossex国产| 国产精品99久久久久久久久| 女生性感内裤真人,穿戴方法视频| 国产在线精品亚洲第一网站| 成人国产麻豆网| 干丝袜人妻中文字幕| 男人舔奶头视频| 久9热在线精品视频| 国产亚洲精品久久久久久毛片| 亚洲欧美日韩卡通动漫| 在线播放无遮挡| 国产精品日韩av在线免费观看| 国产成人a区在线观看| 日韩欧美三级三区| 亚洲第一电影网av| 色噜噜av男人的天堂激情| 在线看三级毛片| 成人国产一区最新在线观看| 久久久色成人| 全区人妻精品视频| 69人妻影院| 18+在线观看网站| 亚洲国产欧洲综合997久久,| 亚洲国产精品久久男人天堂| 国产av在哪里看| 久久久国产成人免费| 麻豆成人av在线观看| 亚洲av第一区精品v没综合| 国产精品野战在线观看| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| av在线观看视频网站免费| 简卡轻食公司| 精品一区二区三区视频在线观看免费| 毛片一级片免费看久久久久 | 91久久精品电影网| 国内精品宾馆在线| 成人特级黄色片久久久久久久| 亚洲五月天丁香| 99久久九九国产精品国产免费| 欧美日韩乱码在线| 麻豆一二三区av精品| 免费高清视频大片| 淫妇啪啪啪对白视频| 高清日韩中文字幕在线| 国产精品永久免费网站| 搡女人真爽免费视频火全软件 | 99热网站在线观看| 午夜福利视频1000在线观看| 天天一区二区日本电影三级| 成人特级av手机在线观看| 99热这里只有精品一区| 欧美+亚洲+日韩+国产| 日韩 亚洲 欧美在线| 国产爱豆传媒在线观看| 国产精品一区www在线观看 | 99riav亚洲国产免费| 网址你懂的国产日韩在线| 亚洲黑人精品在线| 成人高潮视频无遮挡免费网站| 国产高清有码在线观看视频| 特级一级黄色大片| 美女xxoo啪啪120秒动态图| 亚洲国产日韩欧美精品在线观看| 床上黄色一级片| 大又大粗又爽又黄少妇毛片口| 色综合亚洲欧美另类图片| 真实男女啪啪啪动态图| 少妇的逼好多水| 最近最新中文字幕大全电影3| 亚州av有码| 联通29元200g的流量卡| 日韩欧美三级三区| 黄色欧美视频在线观看| 少妇的逼好多水| 99久久精品一区二区三区| 国产成人福利小说| 看片在线看免费视频| 床上黄色一级片| 国产欧美日韩精品一区二区| 女生性感内裤真人,穿戴方法视频| 午夜视频国产福利| 我要看日韩黄色一级片| 欧美+日韩+精品| 日本爱情动作片www.在线观看 | 国产激情偷乱视频一区二区| 亚洲无线在线观看| 午夜免费男女啪啪视频观看 | 国产精品日韩av在线免费观看| 午夜日韩欧美国产| 欧美一区二区亚洲| 日本熟妇午夜| 一个人观看的视频www高清免费观看| 久久香蕉精品热| 国产麻豆成人av免费视频| 精品久久久久久,| 国产 一区 欧美 日韩| 亚洲av成人av| 国产精品人妻久久久久久| 国产一区二区三区在线臀色熟女| 搡老熟女国产l中国老女人| 少妇裸体淫交视频免费看高清| 99久久九九国产精品国产免费| www.www免费av| 91午夜精品亚洲一区二区三区 | 午夜免费成人在线视频| 最新在线观看一区二区三区| 亚洲乱码一区二区免费版| 亚洲欧美精品综合久久99|