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

    玉樹地震震源區(qū)速度結(jié)構(gòu)與余震分布的關(guān)系

    2013-10-08 01:01:44王長在吳建平房立華王未來
    地球物理學(xué)報(bào) 2013年12期
    關(guān)鍵詞:余震結(jié)構(gòu)

    王長在,吳建平,房立華,王未來

    中國地震局地球物理研究所,北京 100081

    1 引 言

    2010年4月14日07時(shí)49分37秒(北京時(shí)間),我國青海省玉樹縣發(fā)生了MS7.1級(jí)地震.根據(jù)中國地震臺(tái)網(wǎng)測定,震中位于33.2°N,96.6°E,震源深度14km.截至2010年5月9日15時(shí)00分玉樹地震共記錄到余震1712個(gè),其中6.0~6.9級(jí)地震1個(gè).此次地震導(dǎo)致了2698人遇難,經(jīng)濟(jì)損失達(dá)到8000億.是繼汶川地震之后又一次有極大社會(huì)影響的地震事件.

    地震發(fā)生之后,許多學(xué)者對(duì)此次地震的地質(zhì)背景、發(fā)震構(gòu)造、發(fā)震機(jī)理等方面進(jìn)行了研究,取得了許多重要的研究成果[1-6].Wang等利用近場和遠(yuǎn)震數(shù)據(jù),估算了2010年玉樹MS7.1級(jí)地震的破裂速度,并反演得到了斷層面上的靜態(tài)滑動(dòng)量分布[7].Tobita等利用PALSAR干涉測量數(shù)據(jù),得到了2010年玉樹地震斷層模型和同震變形特征[8].

    為加強(qiáng)余震監(jiān)測和震源區(qū)的深部結(jié)構(gòu)研究,主震發(fā)生之后,青海省地震局在余震區(qū)布設(shè)了7個(gè)流動(dòng)地震臺(tái)站,中國地震局地球物理研究所在附近布設(shè)了14個(gè)流動(dòng)地震臺(tái)站.這21個(gè)流動(dòng)地震臺(tái)站與周圍的固定臺(tái)站組成了一個(gè)比較好的小孔徑臺(tái)網(wǎng),提高了地震定位的精度,為進(jìn)行精定位及震源區(qū)速度結(jié)構(gòu)研究提供了寶貴的數(shù)據(jù).本文利用這批數(shù)據(jù),應(yīng)用雙差層析成像方法,得到了精定位結(jié)果和震源區(qū)速度結(jié)構(gòu),并探討了玉樹—甘孜斷裂的展布形態(tài),速度結(jié)構(gòu)與余震分布的關(guān)系,為更深入地了解發(fā)震構(gòu)造提供了基礎(chǔ)信息.

    2 區(qū)域地質(zhì)概況

    玉樹震區(qū)處于青藏高原東緣,在巴顏喀拉地塊與羌塘地塊交界的金沙江縫合帶附近(圖1).青藏高原東部由6個(gè)地塊以及中間以5條縫合線相隔,自北而南分別為:祁連—阿爾金縫合帶,祁連山地體,青海南山及北淮陽斷裂帶,東昆侖—柴達(dá)木地體,東昆侖—阿尼瑪卿縫合帶,松潘—甘孜—可可西里地體,金沙江印支期縫合帶,羌塘地體,班公湖—怒江縫合帶,拉薩地體,雅魯藏布江喜馬拉雅期縫合帶和喜馬拉雅地體[9].青藏高原5條狹長的縫合帶是各時(shí)期洋盆消減和地體碰撞的結(jié)果.縫合帶的時(shí)代由北向南依次變新,顯示了亞洲大陸逐漸向南增生.印度次大陸亞洲大陸碰撞,使雅魯藏布江一帶的特提斯殘余海最終封閉,并開始了緩慢的分階段的構(gòu)造運(yùn)動(dòng)及地面隆升過程[10].自中新世中晚期開始,青藏高原內(nèi)部發(fā)生了較大范圍區(qū)域性的走滑運(yùn)動(dòng),地殼內(nèi)發(fā)生了東向水平運(yùn)動(dòng)和垂直隆升運(yùn)動(dòng),從而青藏高原內(nèi)部的構(gòu)造塊體逐漸被擠出,沿東昆侖斷裂帶、玉樹—甘孜—鮮水河斷裂帶、嘉黎斷裂帶等邊界斷裂帶向青藏高原東緣滑移.伴隨這些塊體運(yùn)動(dòng),青藏高原7級(jí)以上強(qiáng)震(例如2001年昆侖山口8.2級(jí)地震、2010年M7.1級(jí)玉樹地震)主要沿大型塊體邊界斷裂帶分布.

    甘孜—玉樹斷裂西起青海治多縣那王草曲塘,經(jīng)當(dāng)江、玉樹、鄧柯、玉隆,至四川甘孜縣城南,全長約500km.斷裂整體呈北西向展布,僅在當(dāng)江附近走向北西西,斷層傾向以北東為主(僅擋拖一帶傾向南西),傾角70°~85°.甘孜—玉樹斷裂帶在第四紀(jì)的強(qiáng)烈活動(dòng),造成了新沉積物的強(qiáng)烈變形.根據(jù)探槽揭露,這種新沉積物的變形大多顯示了脆性破裂的特點(diǎn),是斷裂帶上古地震形變的證據(jù)[11].玉樹地震震區(qū)位于巴顏喀拉地塊南邊界,周圍斷裂構(gòu)造發(fā)育,其中以北西向左旋走滑的斷裂為主.區(qū)域地震活動(dòng)與活動(dòng)斷裂構(gòu)造關(guān)系密切,6級(jí)以上地震主要分布在晚第四紀(jì)以來活動(dòng)斷裂附近.其中這次玉樹M7.1級(jí)地震就是甘孜—玉樹斷裂突然錯(cuò)動(dòng)的結(jié)果.

    3 方法原理

    雙差層析成像方法是在雙差定位法的基礎(chǔ)上發(fā)展起來的,其中雙差定位法(hypoDD)已經(jīng)被國內(nèi)、外地震學(xué)家廣泛地應(yīng)用到地震定位中[14-17],它在確定地震之間相對(duì)位置方面具有很高的精度,是研究特定地區(qū)地震活動(dòng)特征、活動(dòng)斷層空間展布等的重要手段.但雙差層析成像方法由于考慮了介質(zhì)速度結(jié)構(gòu)的空間變化,克服了雙差定位對(duì)臺(tái)站到事件對(duì)之間路徑為恒定速度的假設(shè),因此得到的定位結(jié)果更加精確.

    將雙差層析成像的基本原理進(jìn)行簡單闡述:理論到時(shí)與觀測到時(shí)的差(即殘差)rik與震源參數(shù)的擾動(dòng)量和波速之間的關(guān)系表示為:τi是地震i的發(fā)震時(shí)刻,u是慢度矢量,ds是路徑積分元.其中震源參數(shù)(x1,x2,x3)、發(fā)震時(shí)刻、慢度場、射線路徑是未知量.

    若地震j也被臺(tái)站k所記錄,則有:

    則這兩個(gè)事件與計(jì)算理論走時(shí)差的殘差即雙差:

    通過聯(lián)合反演可以得到三維速度結(jié)構(gòu)、震源的相對(duì)位置和絕對(duì)位置[18].

    雙差層析成像方法是運(yùn)用絕對(duì)走時(shí)和相對(duì)走時(shí)資料來實(shí)現(xiàn)三維波速結(jié)構(gòu)和震源參數(shù)的聯(lián)合反演.該方法首先采用網(wǎng)格節(jié)點(diǎn)法進(jìn)行模型參數(shù)化,通過劃分空間三維網(wǎng)格節(jié)點(diǎn),采用偽彎曲射線追蹤法找到地震波的最小走時(shí)路徑,并計(jì)算理論走時(shí)及走時(shí)對(duì)震源位置和慢度的偏導(dǎo)數(shù),聯(lián)合使用絕對(duì)走時(shí)、雙差走時(shí)數(shù)據(jù)進(jìn)行反演,雙差數(shù)據(jù)主要用于確定震源區(qū)的精細(xì)結(jié)構(gòu),絕對(duì)走時(shí)數(shù)據(jù)主要確定震源區(qū)以外區(qū)域的速度結(jié)構(gòu).采用阻尼最小二乘分解算法求解,在三個(gè)方向采取相同的光滑權(quán)重對(duì)模型進(jìn)行光滑約束,多次迭代直至得到穩(wěn)定的解.

    在反演的過程中,先賦予絕對(duì)走時(shí)較高的權(quán)重,給P波絕對(duì)走時(shí)的權(quán)重為1.0,給予差分?jǐn)?shù)據(jù)P波的權(quán)重為0.1.可以在一維速度模型的基礎(chǔ)上,建立一個(gè)比較大區(qū)域的三維速度結(jié)構(gòu)的結(jié)果,在幾次迭代之后,將絕對(duì)走時(shí)的權(quán)重降低為0.1,差分?jǐn)?shù)據(jù)的權(quán)重增加為1.0,以提高震源區(qū)速度結(jié)構(gòu)的分辨率.在迭代過程中,通過調(diào)整阻尼值,使方程求得的解穩(wěn)定,獲得震源區(qū)重定位和速度結(jié)構(gòu)的結(jié)果.

    該方法聯(lián)合使用絕對(duì)走時(shí)和相對(duì)走時(shí),因此在震源區(qū)外可以得到與傳統(tǒng)層析成像相同的結(jié)果,由于增加了雙差方程,從而能夠反演震源區(qū)精細(xì)的速度結(jié)構(gòu)及地震重新定位結(jié)果,因此可以揭示比傳統(tǒng)方法更多的細(xì)結(jié)構(gòu)信息[19].

    4 余震資料選取及反演模型的建立

    4.1 資料選取

    2010年4月14日玉樹M7.1級(jí)地震發(fā)生后,由于流動(dòng)地震臺(tái)的布設(shè),在一定程度上彌補(bǔ)了固定地震臺(tái)站分布的缺陷,從而為我們對(duì)震源區(qū)的研究提供了寶貴的數(shù)據(jù)及資料.利用中國地震局地球物理研究所數(shù)據(jù)備份中心[20]和中國地震局地球物理研究所震后架設(shè)的流動(dòng)地震臺(tái)站提供的波形數(shù)據(jù),選取2010年4月14日到2010年6月15日期間M≥1.0級(jí)地震的事件波形,進(jìn)行了震相拾取和地震初定位.選擇的每個(gè)地震至少被4個(gè)臺(tái)站記錄到,最終有1108個(gè)地震事件參與反演,其中參加反演的P波絕對(duì)到時(shí)6437個(gè),反演的S波絕對(duì)到時(shí)5035個(gè),相對(duì)到時(shí)資料P波33550個(gè),S波28539個(gè),參與的臺(tái)站有30個(gè).臺(tái)站分布見圖2.

    圖3 反演網(wǎng)格分布圖圖為平面網(wǎng)格劃分,“+”表示劃分的網(wǎng)格節(jié)點(diǎn),三角形代表臺(tái)站,X,Y代表建立的坐標(biāo)方向.Fig.3 Grid map of tomographyTriangle denotes stations.The coordinate system (X,Y)shown in this figure is used throughout this study.

    在研究區(qū)域,建立坐標(biāo)系:X軸垂直于余震方向,Y軸平行于余震方向(圖3),坐標(biāo)原點(diǎn)為106.06°E、25.56°N.在X方向的水平網(wǎng)格節(jié)點(diǎn)劃分為-488.0,-60.0,-40.0,-20.0,0.0,20.0,40.0,60.0,499.0km;在Y方向劃分為 -488.0,-80.0,-60.0,-40.0,-20.0,0.0,20.0,40.0,60.0,499.0km;在Z方向(垂直向下)劃分為-150.0,0.0,5.0,10.0,15.0,20.0,25.0,30.0,449.0km.

    4.2 初始速度模型建立

    為確定玉樹周圍地區(qū)的地殼厚度和波速比,選取了青海省地震局玉樹地震臺(tái)站在2007年1月1日至2010年1月1日觀測期間,在30°~90°震中距范圍內(nèi),震級(jí)在6級(jí)以上的波形數(shù)據(jù),采用最大熵反褶積方法[21]提取遠(yuǎn)震體波接收函數(shù).提取過程中采用的高斯濾波系數(shù)為2.5,水準(zhǔn)量系數(shù)為0.01.經(jīng)過嚴(yán)格的挑選,最終獲得了233條高質(zhì)量的接收函數(shù).將地震事件的震中距和深度轉(zhuǎn)換到水平層慢度,以0.001的慢度步長對(duì)接收函數(shù)進(jìn)行疊加,得到各慢度范圍內(nèi)的平均接收函數(shù).對(duì)玉樹臺(tái)站的多道平均接收函數(shù)進(jìn)行疊加搜索,獲得此臺(tái)站下方的平均地殼厚度及泊松比值.結(jié)果表明,玉樹臺(tái)站附近地殼厚度約為70km,波速比1.72.

    本文反演采用的初始模型參考了該地區(qū)人工地震觀測的結(jié)果[22-23].表1給出了 P波初始速度模型.表中的上界面深度表示在分層速度模型中各層的頂面深度,其層厚度為下一個(gè)上界面深度與該上界面深度之差,左側(cè)為該層對(duì)應(yīng)的速度值.海平面以上部分與第一層的速度相同.地殼厚度和波速比根據(jù)玉樹臺(tái)的接收函數(shù)H-K疊加結(jié)果給出P波速度與S波速度的比值為1.72,地殼厚度為70km.

    表1 玉樹地區(qū)P波初始速度結(jié)構(gòu)Table 1 Initial P wave velocity structure of Yushu

    5 雙差層析成像的結(jié)果與討論

    5.1 余震定位結(jié)果

    共1108次地震參加重新定位,最終獲得995次地震的精定位結(jié)果,重新定位后均方根殘差平均值由原來的0.5s降為了0.1s,通對(duì)觀測到時(shí)和理論到時(shí)殘差的理論估計(jì),震源位置的測定誤差在E—W方向平均為0.3km,在N—S方向平均為0.3km,在垂直方向平均為0.5km.

    圖4a、4b分別為玉樹地震序列地震定位前后的結(jié)果.初定位結(jié)果顯示余震分布相對(duì)分散,重定位后的結(jié)果顯示余震沿NW向成窄條狀分布在斷層的兩側(cè),清晰地勾勒出地震活動(dòng)圖像,揭示了脆性破裂應(yīng)力釋放主要集中于一個(gè)狹窄的區(qū)域內(nèi).余震分布表明在西北端96.2°左右,余震不完全沿甘孜—玉樹斷裂分布,而在甘孜—玉樹斷裂垂直方向也有分布.分析認(rèn)為在該位置可能存在一條與甘孜—玉樹斷裂相垂直的次級(jí)斷裂.周榮軍等[24]研究表明,玉樹—甘孜斷裂在西北端96.2°左右,斷裂發(fā)生了錯(cuò)斷,斷裂并不是在原來的位置繼續(xù)向北西向,而是斷裂向南西向偏移了10多公里.可能正是這一與甘孜—玉樹斷裂相垂直的次級(jí)斷裂的某個(gè)歷史時(shí)期的錯(cuò)動(dòng),導(dǎo)致了玉樹—甘孜斷裂在這一位置發(fā)生偏移.

    劉超等[25]利用全球地震臺(tái)網(wǎng)的寬頻帶地震波形資料,采用P波波形快速反演方法得到的震源機(jī)制解,反演得出,這次地震的矩震級(jí)為Mw7.0,是一次典型的左旋走滑型地震,并判定這次地震事件的發(fā)震斷層為玉樹斷層,走向119°、傾角83°、滑動(dòng)角-2°.余震的展布方向?yàn)?19°與主震的發(fā)震斷層吻合.

    余震深度分布剖面圖(圖5)表明,初始定位震源深度主要分布在3~15km之間,在5~10km之間出現(xiàn)了明顯的橫向排列.重定位之后,地震主要分布在2~30km之間,從西北到東南(-60~60km)震源深度總體上由深變淺,在西北端(-60~-30km)絕大部分余震震源深度分布在10~20km之間,并與(-30~60km)那一段在深度變化上有一個(gè)明顯的區(qū)別.圖5b為Li等[26]利用SAR、光學(xué)成像、地震體波等反演得到的斷層幾何特征和滑動(dòng)量分布.該結(jié)果也顯示出主震的破裂過程位移最大的部分(圖中為深顏色部分),一部分位于主震破裂點(diǎn)附近,另一部分位于距離破裂點(diǎn)30km左右的位置,并且滑動(dòng)量總體上從西北(0km)到東南方向(45km)呈現(xiàn)從深到淺破裂過程,即應(yīng)力在西北方向主要是在深部的調(diào)整,而到東南方向主要是淺部的應(yīng)力釋放.而在(-30km到-60km),震源深度較深,可能為其他次級(jí)斷層破裂導(dǎo)致的.

    5.2 反演得到的速度結(jié)構(gòu)

    5.2.1 速度結(jié)構(gòu)的解釋

    三維速度結(jié)構(gòu)反演結(jié)果的精度通常用節(jié)點(diǎn)周圍射線的空間分布來進(jìn)行評(píng)估.Thurber和Phillips定義DWS(Derivative Weight Sum)反映一個(gè)模型參量周圍平均的相對(duì)射線密度,考慮了射線節(jié)點(diǎn)到節(jié)點(diǎn)的距離加權(quán)以及此節(jié)點(diǎn)周圍射線長度,因此在反映射線密度方面,DWS比不加權(quán)的射線總和更有優(yōu)勢[27].雙差層析成像方法評(píng)估解的可靠性,就是通過DWS值來進(jìn)行評(píng)估的,許多研究者驗(yàn)證當(dāng)DWS>100時(shí),得到的結(jié)果具有較高的可靠性[28-29].圖6、圖7為應(yīng)用雙差層析成像方法反演得到的速度結(jié)構(gòu).其中被白線包圍的部分為DWS>100的區(qū)域,認(rèn)為得到的速度結(jié)構(gòu)結(jié)果可靠,垂直剖面圖中色標(biāo)均表示絕對(duì)速度,為了方便在同一個(gè)水平面對(duì)比速度變化,水平層析成像的色標(biāo)表示速度的擾動(dòng)量.

    深度為z=1km和z=3km淺層速度結(jié)構(gòu)結(jié)果顯示玉樹縣城附近均位于高速異常體內(nèi),在玉樹縣城西南側(cè)存在一明顯的低速體異常,該區(qū)域地質(zhì)圖表明玉樹縣城主要位于碳酸巖鹽的基巖區(qū)域,故表現(xiàn)為高速異常體區(qū)域,在玉樹縣城西南側(cè)有一片較大區(qū)域的第四紀(jì)沉積層,與其對(duì)應(yīng)的位置表現(xiàn)為明顯的低速異常體.上述地表地質(zhì)特征與淺層速度結(jié)構(gòu)的一致性,在一定程度上佐證了雙差層析成像較高的分辨率和可靠性.

    圖4 (a)重定位前地震分布;(b)重定位后地震分布綠點(diǎn)表示地震,黑線為斷裂.Fig.4 (a)Distribution of earthquakes before relocation;(b)Distribution of earthquakes after relocation Dots denote the aftershocks,black lines show the location of the large fault.

    圖5 (a)重定位前余震深度分布;(b)重定位后余震深度分布及Li等[26]反演得到的斷層面上靜態(tài)滑動(dòng)量分布Fig.5 (a)Distributions of aftershock sequences before earthquake relocation;(b)Distributions of aftershock sequences after earthquake relocation.The slip distribution is from an inversion by Li et al.[26]

    玉樹震區(qū)處于巴顏喀拉地塊與羌塘地塊交界的金沙江縫合帶附近,玉樹—甘孜斷裂為巴顏喀拉地塊南邊界.因此玉樹—甘孜斷裂的北部具有巴顏喀拉地塊性質(zhì),南側(cè)代表了羌塘地塊特征.深度為z=6km、z=9km和z=12km速度結(jié)構(gòu)結(jié)果顯示,在中上地殼巴顏喀拉地塊表現(xiàn)為高速體異常,羌塘地塊表現(xiàn)為低速異常.玉樹—甘孜斷裂處于高低速異常過渡帶區(qū)域.

    圖7 AA′、BB′、CC′、DD′、EE′、FF′深度剖面的速度結(jié)構(gòu)圖像黑點(diǎn)表示余震分布,紅色五星為主震.白色虛線為推測斷層位置,白色實(shí)線包圍的區(qū)域?yàn)镈WS值大于100區(qū)域.震源位置剖面給出的是離剖面10km距離內(nèi)的地震.Fig.7 Vertical cross-sections of VPalong lines AA′,BB′,CC′,DD′,EE′,and FF′Dots denote the aftershocks,stars denote epicenters of the main shock.The white dotted line indicates earthquake fault,the region surrounded by the solid white line shows the region of DWS value more than 100.The location profile shows the earthquakes within 10km from the profile.

    5.2.2 余震分布與速度結(jié)構(gòu)的關(guān)系

    深度為9km和12km水平層析成像結(jié)果顯示玉樹—甘孜斷裂的東北側(cè)主要呈現(xiàn)為高速體異常,西南側(cè)表現(xiàn)為低速體異常,余震主要分布在高速低速過渡帶區(qū)域.主震(即地震的初始破裂點(diǎn))發(fā)生在NWW—SEE向的高速體與低速體過渡帶區(qū)域并偏向高速體一側(cè),這一位置往往是應(yīng)力集中的部位,可能導(dǎo)致大地震的發(fā)生.

    AA′剖面結(jié)合BB′、CC′剖面的速度結(jié)構(gòu)結(jié)果顯示在玉樹—甘孜斷裂西北端(AA′剖面即坐標(biāo)-60~-30km位置),在深度為3~8km附近存在一個(gè)較大范圍的高速體異常,余震主要發(fā)生在這一高速異常體的下方.一般高速體代表較堅(jiān)硬的巖體,震后應(yīng)力調(diào)整時(shí),由于受到前方及上方該堅(jiān)硬巖體的阻擋,從而應(yīng)力在這一位置向共軛方向和向下調(diào)整,導(dǎo)致了共軛方向和下方的巖石破裂,隨后在該位置發(fā)生了一次M5.7級(jí)中強(qiáng)地震,之后又發(fā)生了大量余震,且余震深度分布較深,推測為共軛方向次級(jí)斷層的破裂.

    國內(nèi)外研究發(fā)現(xiàn),層析成像得到的高速異常體代表地震發(fā)生的凹凸體[30-32].一般認(rèn)為高速體與地殼的較脆、較強(qiáng)部分巖體有關(guān),這些部分能夠積累更大的孕震能量.Kato等人研究日本新潟縣中部Mw6.6級(jí)地震時(shí),發(fā)現(xiàn)發(fā)生在斷層附近的大部分余震分布在高速體的外圍[32].圖7垂直于斷層方向的BB′、CC′、DD′、EE′和FF′剖面顯示余震均分布在高速體的外圍,在高速體內(nèi)幾乎沒有地震分布.推測中上地殼的高速體通常具有較高的強(qiáng)度,可以積累較強(qiáng)的孕震能量,主震發(fā)生后,高速體內(nèi)積累的彈性能量向周邊釋放,可能是導(dǎo)致高速體周邊余震發(fā)生的主要原因.

    速度結(jié)構(gòu)剖面圖均顯示了在地殼中18~25km深度范圍內(nèi)存在一低速層,與姜枚等的人工探測剖面結(jié)果基本一致[22].Huang和Zhao研究表明,大地震發(fā)生在高低速交界部位偏于高速區(qū)的一側(cè),震源下方存在的低速異常體可能促進(jìn)了地震成核[33-34].產(chǎn)生這一現(xiàn)象的原因可能是低速體難以積累能量,卻易于傳遞能量,足夠尺度的高速塊體有可能積累足以引發(fā)大地震的應(yīng)變能.高低速異常體的過渡地帶,既是應(yīng)力集中的地方,又是介質(zhì)相對(duì)比較脆弱的地方,這樣的環(huán)境具備了積累大量應(yīng)變能的介質(zhì)條件,容易發(fā)生破裂,易于釋放應(yīng)力,因而容易引發(fā)大的地震.

    結(jié)合震源區(qū)速度結(jié)構(gòu)與余震分布推測了斷層面(圖7中白色虛線所示)的產(chǎn)狀,BB′剖面顯示斷層傾向南西,傾角接近80°,穿過主震的CC′剖面,斷層傾向南西,傾角接近85°,DD′剖面斷層傾向南西,傾角接近90°、EE′剖面斷層傾向北東,傾角接近85°和FF′剖面斷層傾向北東,傾角接近80°.綜上所述在震源區(qū)的玉樹—甘孜斷裂,從西南到東北方向,傾向由南西向逐漸反轉(zhuǎn)為北東向,其中反轉(zhuǎn)的位置在EE′剖面附近.

    5.2.3 結(jié)果的分辨率及可靠性分析

    為了檢驗(yàn)速度結(jié)構(gòu)反演的可靠性及空間分辨率,我們進(jìn)行了恢復(fù)性試驗(yàn),恢復(fù)性試驗(yàn)方法的基本原理是:根據(jù)層析成像反演獲得的三維速度模型,計(jì)算與實(shí)際觀測數(shù)據(jù)相同的各震相的理論走時(shí),將計(jì)算獲得的理論走時(shí)作為觀測值,基于初始模型,反演新的三維速度模型,通過對(duì)比新舊三維速度模型的差異,可以了解成像結(jié)果的可靠性.恢復(fù)性試驗(yàn)中采用的阻尼系數(shù)和光滑約束系數(shù)與利用實(shí)際資料進(jìn)行反演時(shí)完全相同,因此可以較好評(píng)估模型空間分辨率.這種恢復(fù)性試驗(yàn)來檢測模型空間分辨率的方法被廣泛應(yīng)用[35-38].

    圖8給出了玉樹地震余震區(qū)速度結(jié)構(gòu)的恢復(fù)性試驗(yàn)結(jié)果,恢復(fù)性試驗(yàn)結(jié)果與用實(shí)際資料反演得到的三維速度結(jié)構(gòu)相比,總體上可以得到良好的恢復(fù),速度異常的形態(tài)基本相同,只是部分異常體的幅度有較小的差異.恢復(fù)性試驗(yàn)結(jié)果表明,利用現(xiàn)有的地震、臺(tái)站得到的觀測數(shù)據(jù),可以較好地揭示該震源區(qū)中上地殼的速度結(jié)構(gòu).

    6 結(jié) 論

    應(yīng)用雙差層析成像方法反演得到了玉樹地震余震的重新定位結(jié)果.重定位結(jié)果顯示余震主要沿NW向成窄帶狀分布在斷層的兩側(cè),但是在斷層西北端(經(jīng)度為96.2°附近),余震偏離玉樹—甘孜斷裂分布,在SW向也有分布,推測可能與南西向次級(jí)斷裂有關(guān).地震定位結(jié)果縱剖面顯示從主震發(fā)生的位置到玉樹縣城,呈現(xiàn)出震源深度從深到淺的變化過程,在玉樹縣城附近發(fā)生的余震較淺,表明玉樹—甘孜斷裂在玉樹縣城附近主要是淺層破裂.

    反演得到的淺層速度結(jié)構(gòu)顯示玉樹縣城位于高速異常體內(nèi),玉樹縣城西南側(cè)存在一明顯的低速體異常,該區(qū)域地表地質(zhì)資料表明玉樹縣城主要位于碳酸巖鹽的基巖區(qū)域,在玉樹縣城西南側(cè)有一片較大區(qū)域的第四紀(jì)沉積層.說明了該地表地質(zhì)特性與淺層速度結(jié)構(gòu)的一致性.6km以下的速度結(jié)構(gòu)表明巴顏喀拉地塊表現(xiàn)為高速體異常,羌塘地塊表現(xiàn)為低速體異常,玉樹—甘孜斷裂處于高低速過渡帶附近.速度結(jié)構(gòu)的剖面圖均顯示了在地殼18km到25km深度范圍內(nèi)存在一低速層.

    穿過主震的水平切片和縱剖面顯示,高低速體的存在對(duì)余震的分布有一定的控制作用:主震發(fā)生在高速體與低速體過渡帶區(qū)域并偏向高速體一側(cè),這一位置往往是應(yīng)力集中的部位,可能導(dǎo)致大地震的發(fā)生.速度結(jié)構(gòu)其他各剖面顯示表明大部分余震分布在高速體的外圍,在高速體內(nèi)幾乎沒有地震分布.這意味著高速體內(nèi)積累的彈性能量并未在高速體內(nèi)釋放,而是在周圍環(huán)境釋放.

    圖8 BB′EE′剖面和深度為3、6、9、12km速度結(jié)構(gòu)水平切片恢復(fù)性試驗(yàn)的結(jié)果Fig.8 Results of RRT test along profile BB′and EE′;Plan view of dVp%image at depths of 3,6,9,and 12km from RRT test

    震源區(qū)的速度結(jié)構(gòu)表明由于受到前方及上方堅(jiān)硬巖體的阻擋,導(dǎo)致了共軛方向和下方的巖石破裂,發(fā)生了一次M5.7級(jí)中強(qiáng)地震,隨后產(chǎn)生了大量余震.這就解釋了余震不沿NW向分布,而在垂直方向分布,余震深度分布在西側(cè)比東南側(cè)要深一些的原因.

    根據(jù)速度結(jié)構(gòu)和重定位結(jié)果推測出從東北到西南方向斷層傾角發(fā)生了從緩到陡再緩的變化過程,其中在EE′剖面附近斷層的傾向從南西向反轉(zhuǎn)為北東向.

    致 謝 感謝張海江博士提供tomoDD程序,感謝中國地震臺(tái)網(wǎng)中心和中國地震局地球物理研究所數(shù)據(jù)備份中心提供震相數(shù)據(jù)和波形數(shù)據(jù).感謝鄭秀芬研究員在數(shù)據(jù)分析上給予的幫助和張?zhí)熘醒芯繂T、蔣長勝博士、彭漢書博士在方法、討論方面給予的指導(dǎo).

    [1] Lin A,Jia D,Rao G,et al.Recurrent morphogenic earthquakes in the past millennium along the strike-slip Yushu fault,Central Tibetan Plateau.Bull.Seism.Soc.Am.,2011,101(6):2755-2764.

    [2] Pei S P,Chen Y S J.Link between seismic velocity structure and the 2010Ms7.1Yushu earthquake,Qinghai,China:evidence from aftershock tomography.Bull.Seism.Soc.Am.,2012,102(1):445-450.

    [3] Zha X J,Dai Z Y.Fault geometry and slip distribution of the 2010Yushu earthquakes inferred from inSAR measurement.Bull.Seism.Soc.Am.,2011,101(4):1951-1958.

    [4] 劉云華,單新建,屈春燕等.青海玉樹Ms7.1級(jí)地震地表形變場特征研究.中國科學(xué)D輯:地球科學(xué),2010,40(10):1310-1320.Liu Y H,Shan X J,Qu C Y,et al.Earthquake deformation field characteristics associated with the 2010YushuMs7.1 earthquake.ScienceinChinaSeriesD:EarthSciences(in Chinese),2010,40(10):1310-1320.

    [5] 陳正位,楊攀新,李智敏等.玉樹7.1級(jí)地震斷裂特征與地震地表破裂帶.第四紀(jì)研究,2010,30(3):628-631.Chen Z W,Yang P X,Li Z M,et al.Characteristics ofMs7.1Yushu earthquake fault and the surface rupture.QuaternarySciences(in Chinese),2010,30(3):628-631.

    [6] 張勇,許力生,陳運(yùn)泰.2010年4月14日青海玉樹地震破裂過程快速反演.地震學(xué)報(bào),2010,32(3):361-365.Zhang Y,Xu L S,Chen Y T.Fast inversion of rupture process for 14April 2010Yushu,Qinghai,earthquake.Acta SeismologicaSinica(in Chinese),2010,32(3):361-365.

    [7] Wang D,Mori J.The 2010Qinghai,China,Earthquake:A moderate earthquake with supershear rupture.Bull.Seism.Soc.Am.,2012,102(1):301-308.

    [8] Tobita M,Nishimura T,Kobayashia T,et al.Estimation of coseismic deformation and a fault model of the 2010Yushu earthquake using PALSAR interferometry data.Earthand PlanetaryScienceLetters,2011,307(3-4):430-438.

    [9] 李建彪,甘衛(wèi)軍,冉勇康等.青藏高原東部構(gòu)造塊體的運(yùn)動(dòng)學(xué)及形變特征分析.西北地震學(xué)報(bào),2006,28(2):97-103.Li J B,Gan W J,Ran Y K,et al.Analysis on kinematics and deformation characteristic of tectonic blocks in Eastern Tibet Plateau.NorthwesternSeismologicalJournal(in Chinese),2006,28(2):97-103.

    [10] 李金臣,潘華,張志中.青海玉樹7.1級(jí)地震構(gòu)造背景.國際地震動(dòng)態(tài),2010,(5):1-5.Li J C,Pan H,Zhang Z Z.The seismotectonic setting of YushuM7.1earthquake in Qinghai Province.RecentDevelopments inWorldSeismology(in Chinese),2010,(5):1-5.

    [11] 彭華,馬秀敏,白嘉啟等.甘孜玉樹斷裂帶第四紀(jì)活動(dòng)特征.地質(zhì)力學(xué)學(xué)報(bào),2006,12(3):295-304.Peng H,Ma X M,Bai J Q,et al.Characteristics of Quaternary activities of the Garzê-Yushu fault zone.JournalofGeomechanics(in Chinese),2006,12(3):295-304.

    [12] 王未來,吳建平,房立華等.2010年玉樹Ms7.1地震及其余震的雙差定位研究.中國科學(xué)D輯:地球科學(xué),2013,42(7):1037-1046.Wang W L,Wu J P,F(xiàn)ang L H,et al.Relocation of the YushuMs7.1earthquake and its aftershocks in 2010from HypoDD.ScienceChina:EarthSciences,2013,56(2):182-191.

    [13] 刁桂苓,王曉山,高國英等.以震源機(jī)制類型劃分汶川、玉樹地震構(gòu)造塊體歸屬.地球物理學(xué)報(bào),2010,53(8):1778-1783.Diao G L,Wang X S,Gao G Y,et al.Tectonic block attribution of Wenchuan and Yushu earthquakes distinguished by focal mechanism type.ChineseJ.Geophys.(in Chinese),2010,53(8):1778-1783.

    [14] Waldhauser F,Ellsworth W.A double-difference earthquake location algorithm:Method and application to the northern Hayward fault,California.Bull.Seism.Soc.Am.,2000,90(6):1353-1368.

    [15] 楊智嫻,陳運(yùn)泰,鄭月軍等.雙差地震定位法在我國中西部地區(qū)地震精確定位中的應(yīng)用.中國科學(xué)(D輯),2003,33(S1):129-134.Yang Z X,Chen Y T,Zheng Y J,et al.Accurate relocation of earthquakes in central-western China using the double difference earthquake location algorithm.ScienceinChina(SeriesD),2003,46(S2):181-188.

    [16] 黃嬡,楊建思,張?zhí)熘?2003年新疆巴楚伽—師地震序列的雙差法重新定位研究.地球物理學(xué)報(bào),2006,49(1):162-169.Huang Y,Yang J S,Zhang T Z.Relocation of the Bachu-Jiashi,Xinjiang earthquake sequence in 2003using the double-difference location algorithm.ChineseJ.Geophys.(in Chinese),2006,49(1):162-169.

    [17] 張?zhí)熘?,武巴特爾,黃媛等.近臺(tái)資料對(duì)近震相對(duì)定位算法的影響.地球物理學(xué)報(bào),2007,50(4):1123-1130.Zhang T Z,Wu B T R,Huang Y,et al.Effect of the data recorded at nearby stations on earthquake relative location.ChineseJ.Geophys.(in Chinese),2007,50(4):1123-1130.

    [18] Zhang H J,Thurber C.Double-difference tomography:the method and its application to the Hayward fault,California.Bull.Seism.Soc.Am.,2003,93(5):1875-1889.

    [19] 王長在,吳建平,房立華等.2009年姚安地震序列定位及震源區(qū)三維P波速度結(jié)構(gòu)研究.地震學(xué)報(bào),2011,33(2):123-133.Wang C Z,Wu J P,F(xiàn)ang L H,et al.Relocation of aftershocks of the 2009YaoanMs6.0earthquake and 3-D P wave velocity structure around its source region.Acta SeismologicaSinica(in Chinese),2011,33(2):123-133.

    [20] Zheng X F,Yao Z X,Liang J H,et al.The role played and opportunities provided by IGP DMC of China National Seismic Network in Wenchuan earthquake disaster relief and researches.Bull.Seism.Soc.Am.,2010,100(5B):2866-2872,doi 10.1785/0120090257.

    [21] 吳慶舉,田小波,張乃鈴等.計(jì)算臺(tái)站接收函數(shù)的最大熵譜反褶積方法.地震學(xué)報(bào),2003,25(4):382-389.Wu Q J,Tian X B,Zhang N L,et al.Receiver function estimated by maximum entropy deconvolution.ActaSeismologica Sinica(in Chinese),2003,25(4):382-389.

    [22] 姜枚,王有學(xué),錢輝等.造山的高原——青藏高原及其鄰區(qū)的寬頻地震探測與地殼上地幔結(jié)構(gòu).北京:地質(zhì)出版社,2009:58-85.Jiang M,Wang Y X,Qian H,et al.Orogenic Plateau-Broadband Seismic Investigations of Crust and Upper Mantle Structure in the Qinghai-Tibet Plateau and Adjacent Areas(in Chinese).Beijing:Geological Publishing House,2009:58-85.

    [23] 王夫運(yùn),張成科,段永紅等.玉樹強(qiáng)震區(qū)地殼結(jié)構(gòu)與深部孕震環(huán)境探測研究.中國地球物理學(xué)會(huì)第二十七屆年會(huì)論文集,2011:97-98.Wang F Y,Zhang C K,Duan Y H,et al.Study on crustal structure and seismogenic condition of Yushu strong earthquake zone(in Chinese).Twenty-Seventh Annual Meeting Symposium of the Chinese Geophysical Society,2011:97-98.

    [24] 周榮軍,馬聲浩,蔡長星.甘孜—玉樹斷裂帶的晚第四紀(jì)活動(dòng)特征.中國地震,1996,12(3):250-260.Zhou R J,Ma S H,Cai C X.Late Quaternary active features of the Ganzi-Yushu fault zone.EarthquakeResearchinChina(in Chinese),1996,12(3):250-260.

    [25] 劉超,許力生,陳運(yùn)泰.2010年4月14日青海玉樹地震快速矩張量解.地震學(xué)報(bào),2010,32(3):366-368.Liu C,Xu L S,Chen Y T.Quick moment tensor solution for 14April 2010Yushu,Qinghai,earthquake.ActaSeismologica Sinica(in Chinese),2010,32(3):366-368.

    [26] Li Z H,Elliott J R,F(xiàn)eng W P,et al.The 2010Mw6.8 Yushu(Qinghai,China)earthquake:Constraints provided by InSAR and body wave seismology.J.Geophys.Res.,2011,116(B10),doi:10.1029/2011JB008358.

    [27] Thurber C,Eberhart-Phillips D.Local earthquake tomography with flexible gridding.Comput.Geosci.,1999,25(7):809-818.

    [28] Thurber C H,Brocher T M,Zhang H J,et al.Threedimensional P wave velocity model for the San Francisco Bay region,California.J.Geophys.Res.,2007,112(B07),doi:10.1029/2006JB004682.

    [29] Scarfi L,Giampiccolo E,Musumeci C,et al.New insights on 3Dcrustal structure in southeastern Sicily (Italy)and tectonic implications from an adaptive mesh seismic tomography.Phys.EarthPlanetInter.,2007,161(1-2):74-85.

    [30] Michael A J,Eberhart-Phillips D.Relations among fault behavior,subsurface geology,and three-dimensional velocity models.Science,1991,253(5020):651-654.

    [31] Chiarabba C,Amato A.VpandVp/Vsimages in theMw6.0 Colfiorito fault region(central Italy):A contribution to the understanding of seismotectonic and seismogenic processes.J.Geophys.Res.,2003,108(B5),doi 10.1029/2001JB001665.

    [32] Kato A,Miyatake T,Hirata N.Asperity and barriers of the 2004Mid-Niigata prefecture earthquake revealed by highly dense seismic observations.Bull.Seism.Soc.Am.,2010,100(1):298-306.

    [33] Huang J,Zhao D.Crustal heterogeneity and seismotectonics of the region around Beijing,China.Tectonophysics,2004,385(1-4):159-180.

    [34] 黃金莉,趙大鵬,鄭斯華.川滇活動(dòng)構(gòu)造區(qū)地震層析成像.地球物理學(xué)報(bào),2001,44(S1):127-136.Huang J L,Zhao D P,Zheng S H.Seismic tomography of the Sichuan-Yunnan active tectonic region.ChineseJ.Geophys.(in Chinese),2001,44(S1):127-136.

    [35] 雷建設(shè),趙大鵬,蘇金蓉等.龍門山斷裂帶地殼精細(xì)結(jié)構(gòu)與汶川地震發(fā)震機(jī)理.地球物理學(xué)報(bào),2009,52(2):339-345.Lei J S,Zhao D P,Su J R,et al.Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuan earthquake.ChineseJ.Geophys.(in Chinese),2009,52(2):339-345.

    [36] Okada T,Hasegawa A,Suganomataa J,at al.Imaging the heterogeneous source area of the 2003M6.4northern Miyagi earthquake,NE Japan,by double-difference tomography.Tectonophysics,2007,430(1-4):67-81.

    [37] Zhao D,Hasegawa A,Horiuchi S.Tomographic imaging of P and S wave velocity structure beneath northeastern Japan.J.Geophys.Res.,2005,97(B13):19909-19928.

    [38] 吳建平,楊婷,王未來等.小江斷裂帶周邊地區(qū)三維P波速度結(jié)構(gòu)及其構(gòu)造意義.地球物理學(xué)報(bào),2013,56(7):2257-2267.Wu J P,Yang T,Wang W L,et al.Three dimensional P-wave velocity structure around Xiaojiang fault system and its tectonic implications.ChineseJ.Geophys.(in Chinese),2013,56(7):2257-2267.

    猜你喜歡
    余震結(jié)構(gòu)
    基于指數(shù)函數(shù)的川滇地區(qū)余震序列衰減規(guī)律研究
    “超長待機(jī)”的余震
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    生死之間的靈魂救贖——《余震》和《云中記》的倫理問題
    阿來研究(2019年2期)2019-03-03 13:35:00
    論《日出》的結(jié)構(gòu)
    本土化改編與再創(chuàng)——從小說《余震》到電影《唐山大地震》
    三次8級(jí)以上大地震的余震活動(dòng)特征分析*
    地震研究(2015年4期)2015-12-25 05:33:44
    創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長
    亚洲国产欧美在线一区| 亚洲国产日韩一区二区| 麻豆av在线久日| 成人国产av品久久久| 国产高清videossex| 国产成+人综合+亚洲专区| 午夜精品国产一区二区电影| 欧美日韩中文字幕国产精品一区二区三区 | 欧美黑人欧美精品刺激| 麻豆av在线久日| 亚洲精品国产一区二区精华液| 色播在线永久视频| 每晚都被弄得嗷嗷叫到高潮| 精品人妻熟女毛片av久久网站| 国产黄色免费在线视频| 制服人妻中文乱码| 妹子高潮喷水视频| 十八禁网站网址无遮挡| 国产熟女午夜一区二区三区| 精品一区二区三区av网在线观看 | 另类精品久久| 99re6热这里在线精品视频| 亚洲精华国产精华精| 大码成人一级视频| 在线十欧美十亚洲十日本专区| 91成人精品电影| 久久女婷五月综合色啪小说| 日本撒尿小便嘘嘘汇集6| 国产成人免费无遮挡视频| 超碰97精品在线观看| 欧美激情 高清一区二区三区| 美女大奶头黄色视频| 久久久久久亚洲精品国产蜜桃av| 高清欧美精品videossex| 新久久久久国产一级毛片| 国产三级黄色录像| 免费在线观看黄色视频的| 精品少妇久久久久久888优播| 国产精品国产三级国产专区5o| 午夜免费观看性视频| 高清黄色对白视频在线免费看| 最新在线观看一区二区三区| 久久久久网色| 久久精品国产亚洲av香蕉五月 | 久久人妻福利社区极品人妻图片| 少妇的丰满在线观看| 天天躁日日躁夜夜躁夜夜| 久久精品国产a三级三级三级| av有码第一页| 亚洲免费av在线视频| 亚洲欧美一区二区三区久久| 免费少妇av软件| 永久免费av网站大全| 美女主播在线视频| 久久精品国产综合久久久| 乱人伦中国视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精华国产精华精| 夫妻午夜视频| 久久久久国产精品人妻一区二区| 在线精品无人区一区二区三| 亚洲国产av影院在线观看| 国产1区2区3区精品| 国产高清视频在线播放一区 | 夜夜骑夜夜射夜夜干| 12—13女人毛片做爰片一| 91精品三级在线观看| 男女边摸边吃奶| 亚洲综合色网址| 老司机影院成人| 久久国产精品男人的天堂亚洲| 亚洲欧美一区二区三区黑人| 欧美精品av麻豆av| 一级黄色大片毛片| 亚洲成av片中文字幕在线观看| 亚洲天堂av无毛| 久久青草综合色| 老汉色∧v一级毛片| 日韩,欧美,国产一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产99久久九九免费精品| 12—13女人毛片做爰片一| 大香蕉久久网| 成人国语在线视频| 国产老妇伦熟女老妇高清| 欧美97在线视频| 日韩人妻精品一区2区三区| 国产精品 欧美亚洲| 一级黄色大片毛片| 日本精品一区二区三区蜜桃| 亚洲国产av新网站| 老司机影院成人| 久久精品国产亚洲av香蕉五月 | 天天操日日干夜夜撸| 欧美中文综合在线视频| 国产精品国产三级国产专区5o| 男女床上黄色一级片免费看| 在线观看免费午夜福利视频| 久久精品国产综合久久久| 波多野结衣av一区二区av| 男女之事视频高清在线观看| 国产精品欧美亚洲77777| 久久人人爽人人片av| a级片在线免费高清观看视频| 久久九九热精品免费| 午夜免费鲁丝| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品自拍成人| 热re99久久国产66热| av在线app专区| 久久人妻熟女aⅴ| 大片电影免费在线观看免费| 另类精品久久| 狠狠婷婷综合久久久久久88av| 亚洲国产欧美在线一区| 999久久久精品免费观看国产| 日韩,欧美,国产一区二区三区| 精品国产一区二区久久| 在线十欧美十亚洲十日本专区| 人妻 亚洲 视频| 汤姆久久久久久久影院中文字幕| 亚洲成人手机| 男女无遮挡免费网站观看| 男女边摸边吃奶| 亚洲欧美日韩高清在线视频 | 老司机影院毛片| 各种免费的搞黄视频| a 毛片基地| 国产一卡二卡三卡精品| av片东京热男人的天堂| 少妇 在线观看| 深夜精品福利| 菩萨蛮人人尽说江南好唐韦庄| 一本综合久久免费| 另类亚洲欧美激情| 成人黄色视频免费在线看| 中文字幕高清在线视频| 考比视频在线观看| 后天国语完整版免费观看| 精品高清国产在线一区| 国产主播在线观看一区二区| 国产极品粉嫩免费观看在线| 欧美日韩av久久| 欧美激情 高清一区二区三区| 在线观看免费午夜福利视频| 777久久人妻少妇嫩草av网站| 亚洲男人天堂网一区| 黄色片一级片一级黄色片| 国产在视频线精品| 人妻一区二区av| 人人妻人人澡人人看| 国产在视频线精品| www.999成人在线观看| 精品久久久久久电影网| 亚洲精品第二区| 欧美日本中文国产一区发布| 欧美日韩亚洲国产一区二区在线观看 | av福利片在线| 一边摸一边抽搐一进一出视频| 国产伦理片在线播放av一区| bbb黄色大片| 两性午夜刺激爽爽歪歪视频在线观看 | 狂野欧美激情性xxxx| a 毛片基地| 亚洲精品美女久久久久99蜜臀| 亚洲av男天堂| 亚洲精品在线美女| 99国产精品一区二区蜜桃av | 黑人猛操日本美女一级片| 超碰97精品在线观看| 最近中文字幕2019免费版| 亚洲精品国产一区二区精华液| 精品国产超薄肉色丝袜足j| 精品久久久久久久毛片微露脸 | 狠狠精品人妻久久久久久综合| 狠狠婷婷综合久久久久久88av| 男女国产视频网站| av天堂在线播放| 母亲3免费完整高清在线观看| 亚洲精品久久午夜乱码| 国产精品久久久av美女十八| 国产黄色免费在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 一级片免费观看大全| 日本av手机在线免费观看| 极品少妇高潮喷水抽搐| 亚洲一区二区三区欧美精品| 99久久99久久久精品蜜桃| 男女无遮挡免费网站观看| 亚洲熟女精品中文字幕| 日韩免费高清中文字幕av| 国产精品久久久av美女十八| 黄频高清免费视频| 国产男女内射视频| 亚洲中文字幕日韩| 国产精品自产拍在线观看55亚洲 | 人妻一区二区av| 后天国语完整版免费观看| 少妇猛男粗大的猛烈进出视频| 国产亚洲精品第一综合不卡| 久久久久国内视频| 91九色精品人成在线观看| 一区福利在线观看| 在线观看免费午夜福利视频| 超碰成人久久| 91精品三级在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品久久久久久精品电影小说| 亚洲精品一二三| 精品国产乱码久久久久久小说| 国产免费一区二区三区四区乱码| 一区二区日韩欧美中文字幕| 视频区图区小说| 欧美精品啪啪一区二区三区 | 亚洲人成电影免费在线| 不卡av一区二区三区| 国产欧美日韩一区二区三 | 99久久99久久久精品蜜桃| 精品一区二区三卡| 久久香蕉激情| 777久久人妻少妇嫩草av网站| 国产精品二区激情视频| 老司机午夜福利在线观看视频 | 中文字幕制服av| 国产成人精品无人区| 看免费av毛片| 久久天堂一区二区三区四区| 亚洲国产av影院在线观看| 97精品久久久久久久久久精品| 美女高潮到喷水免费观看| 日本猛色少妇xxxxx猛交久久| 欧美日韩亚洲高清精品| 中亚洲国语对白在线视频| 少妇猛男粗大的猛烈进出视频| 欧美变态另类bdsm刘玥| 亚洲欧美成人综合另类久久久| 免费日韩欧美在线观看| 美女脱内裤让男人舔精品视频| 制服人妻中文乱码| 亚洲第一av免费看| 国产极品粉嫩免费观看在线| 国产亚洲一区二区精品| 亚洲七黄色美女视频| 国产免费视频播放在线视频| 国产精品国产av在线观看| 99久久精品国产亚洲精品| 91老司机精品| 亚洲精品av麻豆狂野| 伊人亚洲综合成人网| 欧美日韩成人在线一区二区| e午夜精品久久久久久久| 一级毛片女人18水好多| 亚洲专区国产一区二区| av视频免费观看在线观看| 精品国产乱码久久久久久男人| 美女主播在线视频| 淫妇啪啪啪对白视频 | 午夜免费成人在线视频| videosex国产| 99精品久久久久人妻精品| 色视频在线一区二区三区| 后天国语完整版免费观看| 窝窝影院91人妻| 狂野欧美激情性bbbbbb| 丰满人妻熟妇乱又伦精品不卡| 亚洲一区中文字幕在线| 在线观看人妻少妇| 天堂中文最新版在线下载| 日本av免费视频播放| 久久久久精品国产欧美久久久 | 狂野欧美激情性bbbbbb| av视频免费观看在线观看| 性高湖久久久久久久久免费观看| 麻豆av在线久日| 黑人猛操日本美女一级片| 国产视频一区二区在线看| 国产欧美日韩一区二区三 | 欧美久久黑人一区二区| 动漫黄色视频在线观看| 91精品三级在线观看| 国产黄色免费在线视频| 欧美日韩亚洲国产一区二区在线观看 | 在线亚洲精品国产二区图片欧美| 香蕉国产在线看| 亚洲午夜精品一区,二区,三区| 老汉色av国产亚洲站长工具| 色综合欧美亚洲国产小说| 亚洲精品国产一区二区精华液| 一级黄色大片毛片| 久久久久精品人妻al黑| 久久免费观看电影| 宅男免费午夜| 91国产中文字幕| 制服人妻中文乱码| 老汉色∧v一级毛片| 一二三四社区在线视频社区8| 国产av又大| 久久久精品94久久精品| 国内毛片毛片毛片毛片毛片| 日韩一卡2卡3卡4卡2021年| 免费日韩欧美在线观看| 亚洲久久久国产精品| av片东京热男人的天堂| 国产精品久久久久久精品电影小说| 青春草视频在线免费观看| 久久免费观看电影| svipshipincom国产片| 黑人巨大精品欧美一区二区mp4| 午夜福利乱码中文字幕| 亚洲av欧美aⅴ国产| 桃花免费在线播放| 午夜福利视频精品| 国产在线观看jvid| 欧美在线一区亚洲| 久久久久精品国产欧美久久久 | 成人av一区二区三区在线看 | 亚洲国产欧美网| 高清视频免费观看一区二区| 欧美大码av| 日韩精品免费视频一区二区三区| 一级片'在线观看视频| 国产亚洲av片在线观看秒播厂| 日本一区二区免费在线视频| 亚洲 国产 在线| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久人人人人人| 久久精品久久久久久噜噜老黄| 性色av一级| 这个男人来自地球电影免费观看| 午夜两性在线视频| 中文欧美无线码| 亚洲精品一卡2卡三卡4卡5卡 | 91精品国产国语对白视频| 日韩精品免费视频一区二区三区| 狠狠狠狠99中文字幕| 黄色视频不卡| 天天操日日干夜夜撸| 亚洲国产欧美在线一区| 99热全是精品| 亚洲精品粉嫩美女一区| 18禁裸乳无遮挡动漫免费视频| 99久久综合免费| 国产真人三级小视频在线观看| 99国产精品一区二区蜜桃av | 国产一区二区三区综合在线观看| 黄色片一级片一级黄色片| 我要看黄色一级片免费的| 1024香蕉在线观看| 精品福利永久在线观看| 中文精品一卡2卡3卡4更新| 一区二区三区乱码不卡18| 亚洲男人天堂网一区| 爱豆传媒免费全集在线观看| 国产一区有黄有色的免费视频| 成人国产av品久久久| 午夜福利一区二区在线看| 亚洲国产中文字幕在线视频| 欧美97在线视频| 狂野欧美激情性xxxx| 国产成人系列免费观看| 久久久久久亚洲精品国产蜜桃av| 丝袜喷水一区| 2018国产大陆天天弄谢| 久久毛片免费看一区二区三区| 少妇的丰满在线观看| 无遮挡黄片免费观看| 桃花免费在线播放| 久久久久久久国产电影| 色婷婷av一区二区三区视频| 国产在线观看jvid| 91麻豆av在线| 人妻 亚洲 视频| 国产亚洲午夜精品一区二区久久| 国产精品国产av在线观看| 精品一区在线观看国产| 高清视频免费观看一区二区| 亚洲精品av麻豆狂野| 国产三级黄色录像| 女人爽到高潮嗷嗷叫在线视频| 久久精品人人爽人人爽视色| 99久久综合免费| 国产欧美日韩精品亚洲av| 高清av免费在线| 在线观看www视频免费| 国产成人精品在线电影| 狠狠精品人妻久久久久久综合| 久久精品成人免费网站| 下体分泌物呈黄色| 深夜精品福利| 欧美另类一区| 午夜免费成人在线视频| 久久人妻福利社区极品人妻图片| 成人18禁高潮啪啪吃奶动态图| 国产亚洲av片在线观看秒播厂| 免费观看a级毛片全部| 亚洲自偷自拍图片 自拍| 十八禁网站网址无遮挡| 亚洲情色 制服丝袜| 国产不卡av网站在线观看| 伊人亚洲综合成人网| 丰满少妇做爰视频| 亚洲中文av在线| 99香蕉大伊视频| 亚洲专区国产一区二区| 肉色欧美久久久久久久蜜桃| 亚洲欧美日韩高清在线视频 | 国产成人欧美| 男女边摸边吃奶| 我要看黄色一级片免费的| 国产欧美日韩精品亚洲av| www.熟女人妻精品国产| 亚洲成人免费av在线播放| 最黄视频免费看| 欧美精品亚洲一区二区| 大陆偷拍与自拍| 我的亚洲天堂| 另类亚洲欧美激情| 日韩有码中文字幕| 99热全是精品| 69av精品久久久久久 | 久久国产亚洲av麻豆专区| 亚洲第一青青草原| 婷婷丁香在线五月| 欧美午夜高清在线| 啦啦啦在线免费观看视频4| 视频区欧美日本亚洲| 国产亚洲精品一区二区www | 欧美xxⅹ黑人| 日韩 亚洲 欧美在线| 大型av网站在线播放| 精品久久久精品久久久| 久久狼人影院| 亚洲精品美女久久av网站| 日韩有码中文字幕| 97精品久久久久久久久久精品| 国产精品一区二区在线观看99| 大码成人一级视频| 久久久久久人人人人人| 最近最新免费中文字幕在线| 日韩免费高清中文字幕av| 国产一区二区三区综合在线观看| 国产成人啪精品午夜网站| 人人妻,人人澡人人爽秒播| 亚洲美女黄色视频免费看| 丝袜在线中文字幕| 性少妇av在线| av欧美777| av天堂久久9| 19禁男女啪啪无遮挡网站| 王馨瑶露胸无遮挡在线观看| 亚洲精品中文字幕一二三四区 | 在线 av 中文字幕| 欧美午夜高清在线| 十八禁网站免费在线| 99久久国产精品久久久| 亚洲欧美清纯卡通| 精品少妇黑人巨大在线播放| 久久久精品94久久精品| 日本欧美视频一区| 热99re8久久精品国产| 在线观看www视频免费| 欧美亚洲 丝袜 人妻 在线| 亚洲专区字幕在线| 亚洲性夜色夜夜综合| 精品高清国产在线一区| 亚洲中文av在线| 中文字幕av电影在线播放| 久久久久精品人妻al黑| 久久久久久久精品精品| 亚洲男人天堂网一区| 久久综合国产亚洲精品| 亚洲男人天堂网一区| 欧美久久黑人一区二区| 午夜福利,免费看| 91老司机精品| 91精品三级在线观看| 高潮久久久久久久久久久不卡| 国产精品二区激情视频| 黄色怎么调成土黄色| 十八禁网站网址无遮挡| 亚洲国产精品一区三区| 日韩精品免费视频一区二区三区| 老司机午夜福利在线观看视频 | 天堂8中文在线网| 亚洲国产看品久久| 国产一区二区在线观看av| 国产视频一区二区在线看| 男女之事视频高清在线观看| av在线播放精品| 久久中文看片网| 亚洲免费av在线视频| 丝袜在线中文字幕| 亚洲一区二区三区欧美精品| 男人爽女人下面视频在线观看| 女人久久www免费人成看片| 国产极品粉嫩免费观看在线| 免费在线观看影片大全网站| 欧美日韩成人在线一区二区| 曰老女人黄片| 天堂中文最新版在线下载| 国产亚洲精品一区二区www | 91九色精品人成在线观看| 国产欧美亚洲国产| 青春草亚洲视频在线观看| 国产精品国产av在线观看| 一级片'在线观看视频| 在线观看www视频免费| 亚洲精品日韩在线中文字幕| 在线 av 中文字幕| 久久中文字幕一级| 午夜福利影视在线免费观看| 亚洲九九香蕉| 黑人操中国人逼视频| 亚洲性夜色夜夜综合| 99久久99久久久精品蜜桃| 大陆偷拍与自拍| 亚洲欧美一区二区三区久久| 国产伦人伦偷精品视频| 国产老妇伦熟女老妇高清| 最近最新免费中文字幕在线| 国产亚洲一区二区精品| 如日韩欧美国产精品一区二区三区| 丰满饥渴人妻一区二区三| 麻豆国产av国片精品| 亚洲欧美激情在线| 人人妻人人爽人人添夜夜欢视频| 国产人伦9x9x在线观看| 国产精品香港三级国产av潘金莲| 各种免费的搞黄视频| 丝瓜视频免费看黄片| bbb黄色大片| 亚洲全国av大片| 日韩有码中文字幕| 亚洲 国产 在线| 久久久精品国产亚洲av高清涩受| 亚洲精品国产av蜜桃| 欧美精品人与动牲交sv欧美| 曰老女人黄片| 久久久久久久久久久久大奶| a 毛片基地| 亚洲精华国产精华精| 亚洲欧美日韩另类电影网站| 这个男人来自地球电影免费观看| 黄色视频在线播放观看不卡| 亚洲国产精品一区二区三区在线| 热99国产精品久久久久久7| 成年美女黄网站色视频大全免费| 高清欧美精品videossex| 亚洲欧美一区二区三区黑人| 老司机亚洲免费影院| 麻豆av在线久日| 另类精品久久| 三级毛片av免费| 另类亚洲欧美激情| 美国免费a级毛片| 国产黄频视频在线观看| 免费少妇av软件| 中文字幕av电影在线播放| 国产精品 国内视频| 亚洲精品一区蜜桃| 国产精品一区二区免费欧美 | 我要看黄色一级片免费的| 波多野结衣一区麻豆| 国产淫语在线视频| 国产人伦9x9x在线观看| 中文字幕人妻丝袜制服| 亚洲精品国产区一区二| 久久久国产一区二区| 欧美激情极品国产一区二区三区| 曰老女人黄片| 菩萨蛮人人尽说江南好唐韦庄| 黄片小视频在线播放| 国产精品熟女久久久久浪| 日韩有码中文字幕| 久久天堂一区二区三区四区| 亚洲自偷自拍图片 自拍| 99精品久久久久人妻精品| 久久久精品免费免费高清| 欧美乱码精品一区二区三区| 国产精品成人在线| 精品视频人人做人人爽| 精品卡一卡二卡四卡免费| 国产成人影院久久av| 午夜影院在线不卡| 99国产精品99久久久久| 免费高清在线观看视频在线观看| 天天影视国产精品| 日日爽夜夜爽网站| 亚洲精品久久久久久婷婷小说| 天天躁狠狠躁夜夜躁狠狠躁| 一级毛片精品| 1024视频免费在线观看| 国产亚洲av片在线观看秒播厂| 脱女人内裤的视频| 免费看十八禁软件| 国产精品麻豆人妻色哟哟久久| av在线播放精品| 亚洲av国产av综合av卡| 丝袜喷水一区| 国产成人免费无遮挡视频| 欧美亚洲 丝袜 人妻 在线| 亚洲专区中文字幕在线| 少妇的丰满在线观看| 少妇被粗大的猛进出69影院| 欧美大码av| 在线天堂中文资源库| 91九色精品人成在线观看| 午夜福利在线观看吧| 亚洲三区欧美一区| 日韩制服骚丝袜av| 男女无遮挡免费网站观看| 亚洲欧美日韩高清在线视频 | 亚洲七黄色美女视频| 国产亚洲一区二区精品| 国产成人免费观看mmmm|