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

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

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

    王長(zhǎng)在,吳建平,房立華,王未來(lái)

    中國(guó)地震局地球物理研究所,北京 100081

    1 引 言

    2010年4月14日07時(shí)49分37秒(北京時(shí)間),我國(guó)青海省玉樹(shù)縣發(fā)生了MS7.1級(jí)地震.根據(jù)中國(guó)地震臺(tái)網(wǎng)測(cè)定,震中位于33.2°N,96.6°E,震源深度14km.截至2010年5月9日15時(shí)00分玉樹(shù)地震共記錄到余震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等利用近場(chǎng)和遠(yuǎn)震數(shù)據(jù),估算了2010年玉樹(shù)MS7.1級(jí)地震的破裂速度,并反演得到了斷層面上的靜態(tài)滑動(dòng)量分布[7].Tobita等利用PALSAR干涉測(cè)量數(shù)據(jù),得到了2010年玉樹(shù)地震斷層模型和同震變形特征[8].

    為加強(qiáng)余震監(jiān)測(cè)和震源區(qū)的深部結(jié)構(gòu)研究,主震發(fā)生之后,青海省地震局在余震區(qū)布設(shè)了7個(gè)流動(dòng)地震臺(tái)站,中國(guó)地震局地球物理研究所在附近布設(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),并探討了玉樹(shù)—甘孜斷裂的展布形態(tài),速度結(jié)構(gòu)與余震分布的關(guān)系,為更深入地了解發(fā)震構(gòu)造提供了基礎(chǔ)信息.

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

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

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

    3 方法原理

    雙差層析成像方法是在雙差定位法的基礎(chǔ)上發(fā)展起來(lái)的,其中雙差定位法(hypoDD)已經(jīng)被國(guó)內(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)行簡(jiǎn)單闡述:理論到時(shí)與觀測(cè)到時(shí)的差(即殘差)rik與震源參數(shù)的擾動(dòng)量和波速之間的關(guān)系表示為:τi是地震i的發(fā)震時(shí)刻,u是慢度矢量,ds是路徑積分元.其中震源參數(shù)(x1,x2,x3)、發(fā)震時(shí)刻、慢度場(chǎng)、射線路徑是未知量.

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

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

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

    雙差層析成像方法是運(yùn)用絕對(duì)走時(shí)和相對(duì)走時(shí)資料來(lái)實(shí)現(xiàn)三維波速結(jié)構(gòu)和震源參數(shù)的聯(lián)合反演.該方法首先采用網(wǎng)格節(jié)點(diǎn)法進(jìn)行模型參數(shù)化,通過(guò)劃分空間三維網(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)定的解.

    在反演的過(guò)程中,先賦予絕對(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)的分辨率.在迭代過(guò)程中,通過(guò)調(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日玉樹(shù)M7.1級(jí)地震發(fā)生后,由于流動(dòng)地震臺(tái)的布設(shè),在一定程度上彌補(bǔ)了固定地震臺(tái)站分布的缺陷,從而為我們對(duì)震源區(qū)的研究提供了寶貴的數(shù)據(jù)及資料.利用中國(guó)地震局地球物理研究所數(shù)據(jù)備份中心[20]和中國(guó)地震局地球物理研究所震后架設(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)站分布見(jiàn)圖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 初始速度模型建立

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

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

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

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

    5.1 余震定位結(jié)果

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

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

    劉超等[25]利用全球地震臺(tái)網(wǎng)的寬頻帶地震波形資料,采用P波波形快速反演方法得到的震源機(jī)制解,反演得出,這次地震的矩震級(jí)為Mw7.0,是一次典型的左旋走滑型地震,并判定這次地震事件的發(fā)震斷層為玉樹(shù)斷層,走向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為L(zhǎng)i等[26]利用SAR、光學(xué)成像、地震體波等反演得到的斷層幾何特征和滑動(dòng)量分布.該結(jié)果也顯示出主震的破裂過(guò)程位移最大的部分(圖中為深顏色部分),一部分位于主震破裂點(diǎn)附近,另一部分位于距離破裂點(diǎn)30km左右的位置,并且滑動(dòng)量總體上從西北(0km)到東南方向(45km)呈現(xiàn)從深到淺破裂過(guò)程,即應(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)周圍射線的空間分布來(lái)進(jìn)行評(píng)估.Thurber和Phillips定義DWS(Derivative Weight Sum)反映一個(gè)模型參量周圍平均的相對(duì)射線密度,考慮了射線節(jié)點(diǎn)到節(jié)點(diǎn)的距離加權(quán)以及此節(jié)點(diǎn)周圍射線長(zhǎng)度,因此在反映射線密度方面,DWS比不加權(quán)的射線總和更有優(yōu)勢(shì)[27].雙差層析成像方法評(píng)估解的可靠性,就是通過(guò)DWS值來(lái)進(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é)果顯示玉樹(shù)縣城附近均位于高速異常體內(nèi),在玉樹(shù)縣城西南側(cè)存在一明顯的低速體異常,該區(qū)域地質(zhì)圖表明玉樹(shù)縣城主要位于碳酸巖鹽的基巖區(qū)域,故表現(xiàn)為高速異常體區(qū)域,在玉樹(shù)縣城西南側(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]

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

    圖7 AA′、BB′、CC′、DD′、EE′、FF′深度剖面的速度結(jié)構(gòu)圖像黑點(diǎn)表示余震分布,紅色五星為主震.白色虛線為推測(cè)斷層位置,白色實(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é)果顯示玉樹(shù)—甘孜斷裂的東北側(cè)主要呈現(xiàn)為高速體異常,西南側(cè)表現(xiàn)為低速體異常,余震主要分布在高速低速過(guò)渡帶區(qū)域.主震(即地震的初始破裂點(diǎn))發(fā)生在NWW—SEE向的高速體與低速體過(guò)渡帶區(qū)域并偏向高速體一側(cè),這一位置往往是應(yīng)力集中的部位,可能導(dǎo)致大地震的發(fā)生.

    AA′剖面結(jié)合BB′、CC′剖面的速度結(jié)構(gòu)結(jié)果顯示在玉樹(shù)—甘孜斷裂西北端(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ā)生了大量余震,且余震深度分布較深,推測(cè)為共軛方向次級(jí)斷層的破裂.

    國(guó)內(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)幾乎沒(méi)有地震分布.推測(cè)中上地殼的高速體通常具有較高的強(qiáng)度,可以積累較強(qiáng)的孕震能量,主震發(fā)生后,高速體內(nèi)積累的彈性能量向周邊釋放,可能是導(dǎo)致高速體周邊余震發(fā)生的主要原因.

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

    結(jié)合震源區(qū)速度結(jié)構(gòu)與余震分布推測(cè)了斷層面(圖7中白色虛線所示)的產(chǎn)狀,BB′剖面顯示斷層傾向南西,傾角接近80°,穿過(guò)主震的CC′剖面,斷層傾向南西,傾角接近85°,DD′剖面斷層傾向南西,傾角接近90°、EE′剖面斷層傾向北東,傾角接近85°和FF′剖面斷層傾向北東,傾角接近80°.綜上所述在震源區(qū)的玉樹(shù)—甘孜斷裂,從西南到東北方向,傾向由南西向逐漸反轉(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í)際觀測(cè)數(shù)據(jù)相同的各震相的理論走時(shí),將計(jì)算獲得的理論走時(shí)作為觀測(cè)值,基于初始模型,反演新的三維速度模型,通過(guò)對(duì)比新舊三維速度模型的差異,可以了解成像結(jié)果的可靠性.恢復(fù)性試驗(yàn)中采用的阻尼系數(shù)和光滑約束系數(shù)與利用實(shí)際資料進(jìn)行反演時(shí)完全相同,因此可以較好評(píng)估模型空間分辨率.這種恢復(fù)性試驗(yàn)來(lái)檢測(cè)模型空間分辨率的方法被廣泛應(yīng)用[35-38].

    圖8給出了玉樹(shù)地震余震區(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)站得到的觀測(cè)數(shù)據(jù),可以較好地揭示該震源區(qū)中上地殼的速度結(jié)構(gòu).

    6 結(jié) 論

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

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

    穿過(guò)主震的水平切片和縱剖面顯示,高低速體的存在對(duì)余震的分布有一定的控制作用:主震發(fā)生在高速體與低速體過(guò)渡帶區(qū)域并偏向高速體一側(cè),這一位置往往是應(yīng)力集中的部位,可能導(dǎo)致大地震的發(fā)生.速度結(jié)構(gòu)其他各剖面顯示表明大部分余震分布在高速體的外圍,在高速體內(nèi)幾乎沒(mé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é)果推測(cè)出從東北到西南方向斷層傾角發(fā)生了從緩到陡再緩的變化過(guò)程,其中在EE′剖面附近斷層的傾向從南西向反轉(zhuǎn)為北東向.

    致 謝 感謝張海江博士提供tomoDD程序,感謝中國(guó)地震臺(tái)網(wǎng)中心和中國(guó)地震局地球物理研究所數(shù)據(jù)備份中心提供震相數(shù)據(jù)和波形數(shù)據(jù).感謝鄭秀芬研究員在數(shù)據(jù)分析上給予的幫助和張?zhí)熘醒芯繂T、蔣長(zhǎng)勝博士、彭漢書博士在方法、討論方面給予的指導(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] 劉云華,單新建,屈春燕等.青海玉樹(shù)Ms7.1級(jí)地震地表形變場(chǎng)特征研究.中國(guó)科學(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] 陳正位,楊攀新,李智敏等.玉樹(shù)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日青海玉樹(shù)地震破裂過(guò)程快速反演.地震學(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] 李金臣,潘華,張志中.青海玉樹(shù)7.1級(jí)地震構(gòu)造背景.國(guó)際地震動(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] 彭華,馬秀敏,白嘉啟等.甘孜玉樹(shù)斷裂帶第四紀(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] 王未來(lái),吳建平,房立華等.2010年玉樹(shù)Ms7.1地震及其余震的雙差定位研究.中國(guó)科學(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] 刁桂苓,王曉山,高國(guó)英等.以震源機(jī)制類型劃分汶川、玉樹(shù)地震構(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)泰,鄭月軍等.雙差地震定位法在我國(guó)中西部地區(qū)地震精確定位中的應(yīng)用.中國(guó)科學(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] 王長(zhǎng)在,吳建平,房立華等.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ū)的寬頻地震探測(cè)與地殼上地幔結(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),張成科,段永紅等.玉樹(shù)強(qiáng)震區(qū)地殼結(jié)構(gòu)與深部孕震環(huán)境探測(cè)研究.中國(guó)地球物理學(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] 周榮軍,馬聲浩,蔡長(zhǎng)星.甘孜—玉樹(shù)斷裂帶的晚第四紀(jì)活動(dòng)特征.中國(guó)地震,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日青海玉樹(shù)地震快速矩張量解.地震學(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] 吳建平,楊婷,王未來(lái)等.小江斷裂帶周邊地區(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ī)律研究
    “超長(zhǎng)待機(jī)”的余震
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    論結(jié)構(gòu)
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    生死之間的靈魂救贖——《余震》和《云中記》的倫理問(wèn)題
    論《日出》的結(jié)構(gòu)
    本土化改編與再創(chuàng)——從小說(shuō)《余震》到電影《唐山大地震》
    三次8級(jí)以上大地震的余震活動(dòng)特征分析*
    地震研究(2015年4期)2015-12-25 05:33:44
    創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長(zhǎng)
    亚洲精品国产成人久久av| 亚洲国产精品一区三区| 国产在线视频一区二区| 伊人久久国产一区二区| 精品卡一卡二卡四卡免费| 亚洲av男天堂| 久久热精品热| 久久久久久久久久成人| 中文字幕免费在线视频6| 亚洲四区av| 成人18禁高潮啪啪吃奶动态图 | 免费观看的影片在线观看| 亚洲av.av天堂| 天天躁夜夜躁狠狠久久av| 欧美激情国产日韩精品一区| 久久精品国产鲁丝片午夜精品| 在线 av 中文字幕| 久久久久精品久久久久真实原创| 男女免费视频国产| 亚洲精品一二三| 精品视频人人做人人爽| 如何舔出高潮| 国产黄片视频在线免费观看| 国产精品人妻久久久影院| 在线观看免费日韩欧美大片 | 欧美3d第一页| 夜夜骑夜夜射夜夜干| 最近手机中文字幕大全| 一二三四中文在线观看免费高清| 国产亚洲午夜精品一区二区久久| 日韩大片免费观看网站| 女人精品久久久久毛片| 亚洲精品日韩在线中文字幕| h日本视频在线播放| 两个人免费观看高清视频 | 九九在线视频观看精品| 午夜视频国产福利| 成年av动漫网址| 色婷婷久久久亚洲欧美| 人人澡人人妻人| 日韩大片免费观看网站| 日本av免费视频播放| 男的添女的下面高潮视频| 成人亚洲精品一区在线观看| 精品人妻熟女毛片av久久网站| 两个人免费观看高清视频 | 精华霜和精华液先用哪个| 久久久午夜欧美精品| 国产成人精品福利久久| 两个人免费观看高清视频 | 一区二区av电影网| 国产亚洲5aaaaa淫片| 精品熟女少妇av免费看| 人人妻人人爽人人添夜夜欢视频 | 国产午夜精品一二区理论片| 国产91av在线免费观看| 熟女人妻精品中文字幕| av国产精品久久久久影院| 少妇人妻久久综合中文| 日本wwww免费看| 亚洲伊人久久精品综合| 色吧在线观看| 爱豆传媒免费全集在线观看| 国产成人91sexporn| 一级毛片 在线播放| 欧美日韩视频精品一区| 精品国产一区二区三区久久久樱花| 老女人水多毛片| 日本黄色片子视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品国产色婷婷电影| 日韩伦理黄色片| 我的老师免费观看完整版| 久久精品国产亚洲av涩爱| 免费人成在线观看视频色| 最近的中文字幕免费完整| 观看美女的网站| 日韩中文字幕视频在线看片| 免费观看a级毛片全部| 插逼视频在线观看| 岛国毛片在线播放| 少妇的逼水好多| 黄色欧美视频在线观看| 国产一级毛片在线| 成年人午夜在线观看视频| 久久国产亚洲av麻豆专区| 日韩强制内射视频| 亚洲av.av天堂| av福利片在线| 亚洲av国产av综合av卡| 国产成人精品久久久久久| 18+在线观看网站| 国产日韩欧美在线精品| 高清在线视频一区二区三区| 狂野欧美白嫩少妇大欣赏| 热re99久久国产66热| 成人18禁高潮啪啪吃奶动态图 | 黄色欧美视频在线观看| 蜜桃久久精品国产亚洲av| a级毛色黄片| 国产 精品1| 国内精品宾馆在线| 黄色毛片三级朝国网站 | 十八禁高潮呻吟视频 | 亚洲一区二区三区欧美精品| 你懂的网址亚洲精品在线观看| 看十八女毛片水多多多| 日日摸夜夜添夜夜添av毛片| 亚洲丝袜综合中文字幕| 国产乱人偷精品视频| 美女视频免费永久观看网站| 久久久久国产精品人妻一区二区| 人妻人人澡人人爽人人| 在线观看美女被高潮喷水网站| 精品一区在线观看国产| 免费久久久久久久精品成人欧美视频 | 亚洲,欧美,日韩| 午夜免费男女啪啪视频观看| 日本欧美视频一区| 丝袜喷水一区| 乱码一卡2卡4卡精品| 亚洲不卡免费看| 熟妇人妻不卡中文字幕| 日本黄色片子视频| 青春草国产在线视频| 水蜜桃什么品种好| av不卡在线播放| 插阴视频在线观看视频| 成人无遮挡网站| 久久午夜福利片| 黄色毛片三级朝国网站 | 熟女av电影| 两个人免费观看高清视频 | 亚洲欧洲国产日韩| 在线观看免费日韩欧美大片 | 国产黄色免费在线视频| 久久久久国产网址| 久久精品国产亚洲av涩爱| 少妇高潮的动态图| www.av在线官网国产| 日本黄色日本黄色录像| 99re6热这里在线精品视频| 国产探花极品一区二区| 久久亚洲国产成人精品v| 欧美人与善性xxx| 欧美激情国产日韩精品一区| 黄色一级大片看看| 亚洲中文av在线| 九草在线视频观看| 五月开心婷婷网| 成人毛片a级毛片在线播放| 国产又色又爽无遮挡免| 国产av一区二区精品久久| 免费人成在线观看视频色| 精品视频人人做人人爽| 人妻 亚洲 视频| 久热这里只有精品99| 成年人免费黄色播放视频 | 亚洲天堂av无毛| 免费人妻精品一区二区三区视频| 日韩av在线免费看完整版不卡| 久久99蜜桃精品久久| 久久久久精品久久久久真实原创| 少妇丰满av| 美女主播在线视频| 3wmmmm亚洲av在线观看| 亚洲伊人久久精品综合| 日本猛色少妇xxxxx猛交久久| 一本—道久久a久久精品蜜桃钙片| 99久久中文字幕三级久久日本| 十分钟在线观看高清视频www | 嫩草影院入口| 最近中文字幕高清免费大全6| 3wmmmm亚洲av在线观看| 精品视频人人做人人爽| 王馨瑶露胸无遮挡在线观看| 一级二级三级毛片免费看| 国产男女内射视频| 97在线视频观看| 成人影院久久| 亚洲精品国产色婷婷电影| 欧美性感艳星| 伦理电影免费视频| av天堂久久9| 最近中文字幕高清免费大全6| 日本av手机在线免费观看| 国产综合精华液| 亚洲精品日本国产第一区| 热99国产精品久久久久久7| 亚洲人与动物交配视频| 欧美丝袜亚洲另类| 丝瓜视频免费看黄片| 国产亚洲av片在线观看秒播厂| 久久亚洲国产成人精品v| 亚洲久久久国产精品| 欧美+日韩+精品| 精品酒店卫生间| 免费av不卡在线播放| 国产有黄有色有爽视频| 一级毛片黄色毛片免费观看视频| 久久久久久久久久人人人人人人| 肉色欧美久久久久久久蜜桃| 亚洲精品aⅴ在线观看| 日韩伦理黄色片| 国产成人免费观看mmmm| 一级片'在线观看视频| 97在线视频观看| 亚洲av电影在线观看一区二区三区| 91久久精品国产一区二区成人| 亚洲高清免费不卡视频| 插阴视频在线观看视频| 男人爽女人下面视频在线观看| 51国产日韩欧美| 国产91av在线免费观看| 99久久人妻综合| 少妇精品久久久久久久| 赤兔流量卡办理| 亚洲在久久综合| 大陆偷拍与自拍| 久热这里只有精品99| 欧美老熟妇乱子伦牲交| 伦精品一区二区三区| 成人亚洲精品一区在线观看| 精品人妻偷拍中文字幕| 亚洲一区二区三区欧美精品| 国产一区二区三区av在线| 观看av在线不卡| 成人国产av品久久久| 久久久亚洲精品成人影院| 亚洲欧洲国产日韩| 大片电影免费在线观看免费| 免费观看性生交大片5| 欧美性感艳星| av国产久精品久网站免费入址| 下体分泌物呈黄色| 欧美日韩亚洲高清精品| 精品一品国产午夜福利视频| 亚洲欧美日韩东京热| 永久网站在线| 国产真实伦视频高清在线观看| 亚州av有码| 黄色毛片三级朝国网站 | 欧美精品一区二区大全| 久久精品久久久久久噜噜老黄| 国产亚洲欧美精品永久| 国产在视频线精品| 成人特级av手机在线观看| 夫妻午夜视频| 99九九在线精品视频 | 插逼视频在线观看| 一区二区三区精品91| 日韩不卡一区二区三区视频在线| 国产亚洲午夜精品一区二区久久| 黄色怎么调成土黄色| 性色avwww在线观看| 国产成人aa在线观看| 三上悠亚av全集在线观看 | 少妇被粗大猛烈的视频| 爱豆传媒免费全集在线观看| 久久精品久久精品一区二区三区| 午夜福利影视在线免费观看| 丰满饥渴人妻一区二区三| 99久久人妻综合| 蜜桃在线观看..| 夜夜看夜夜爽夜夜摸| 成人黄色视频免费在线看| 爱豆传媒免费全集在线观看| 国产综合精华液| 色婷婷久久久亚洲欧美| 久久久久久人妻| 色5月婷婷丁香| 欧美日韩综合久久久久久| 国产黄色免费在线视频| 国产精品国产三级国产av玫瑰| 亚洲情色 制服丝袜| 国产老妇伦熟女老妇高清| 亚洲av成人精品一区久久| av有码第一页| 国产黄频视频在线观看| 九色成人免费人妻av| 亚洲婷婷狠狠爱综合网| 国产老妇伦熟女老妇高清| 中文字幕亚洲精品专区| 人人妻人人爽人人添夜夜欢视频 | 观看av在线不卡| 国产成人精品婷婷| 王馨瑶露胸无遮挡在线观看| 国产亚洲一区二区精品| 一级毛片aaaaaa免费看小| 精品国产一区二区三区久久久樱花| 欧美日韩在线观看h| 成人特级av手机在线观看| 色婷婷久久久亚洲欧美| 中文在线观看免费www的网站| 国产在线视频一区二区| 亚洲一区二区三区欧美精品| 大片电影免费在线观看免费| 午夜91福利影院| 亚洲精品中文字幕在线视频 | 人妻夜夜爽99麻豆av| 国产成人免费观看mmmm| 黄色怎么调成土黄色| 国产美女午夜福利| 一个人免费看片子| 夫妻性生交免费视频一级片| 日本av免费视频播放| 国产日韩一区二区三区精品不卡 | 少妇人妻久久综合中文| 在线观看美女被高潮喷水网站| 亚洲精品日本国产第一区| 免费看日本二区| 成人综合一区亚洲| 亚洲成色77777| 国产永久视频网站| 亚洲伊人久久精品综合| 国产高清有码在线观看视频| 国产老妇伦熟女老妇高清| 久久精品国产a三级三级三级| 亚洲无线观看免费| 亚洲av日韩在线播放| 菩萨蛮人人尽说江南好唐韦庄| 人人澡人人妻人| 亚洲欧美一区二区三区国产| 国产日韩欧美在线精品| .国产精品久久| 亚洲美女黄色视频免费看| 国国产精品蜜臀av免费| 欧美日韩av久久| 亚洲国产毛片av蜜桃av| 人人妻人人澡人人爽人人夜夜| 深夜a级毛片| 搡老乐熟女国产| 夜夜骑夜夜射夜夜干| 日本91视频免费播放| 女性生殖器流出的白浆| 亚洲熟女精品中文字幕| 日韩强制内射视频| 2021少妇久久久久久久久久久| 久久精品国产亚洲av天美| 国产精品一区二区三区四区免费观看| .国产精品久久| 在线观看av片永久免费下载| 国产视频首页在线观看| 女性被躁到高潮视频| 中文字幕免费在线视频6| 国产毛片在线视频| 亚洲真实伦在线观看| 国产精品人妻久久久影院| 国模一区二区三区四区视频| 亚洲精品国产av蜜桃| 国产成人freesex在线| 亚洲美女黄色视频免费看| 丰满迷人的少妇在线观看| 韩国高清视频一区二区三区| 免费少妇av软件| 亚洲欧美一区二区三区黑人 | 日韩伦理黄色片| av国产精品久久久久影院| 国产精品伦人一区二区| 99re6热这里在线精品视频| 久久久久人妻精品一区果冻| 日韩免费高清中文字幕av| 精品99又大又爽又粗少妇毛片| 亚洲在久久综合| 国产精品无大码| 亚洲国产精品一区三区| 啦啦啦啦在线视频资源| 久久这里有精品视频免费| 91久久精品电影网| 国国产精品蜜臀av免费| 亚洲欧美一区二区三区国产| 久热这里只有精品99| 国产精品成人在线| 一级毛片久久久久久久久女| 色视频www国产| 欧美老熟妇乱子伦牲交| 日本与韩国留学比较| 欧美日韩视频高清一区二区三区二| 久久99热6这里只有精品| 国产av一区二区精品久久| 男女边摸边吃奶| av一本久久久久| 免费黄色在线免费观看| 蜜桃在线观看..| 国产欧美日韩综合在线一区二区 | 亚洲国产精品成人久久小说| 国产乱人偷精品视频| 久久久国产精品麻豆| 99热网站在线观看| 春色校园在线视频观看| 日韩视频在线欧美| 这个男人来自地球电影免费观看 | 日韩在线高清观看一区二区三区| av网站免费在线观看视频| 久久久久精品久久久久真实原创| 黑人猛操日本美女一级片| 99热全是精品| 久久精品国产亚洲av涩爱| 一级毛片黄色毛片免费观看视频| 精品国产乱码久久久久久小说| 大片电影免费在线观看免费| 国产成人免费无遮挡视频| 高清毛片免费看| 你懂的网址亚洲精品在线观看| xxx大片免费视频| 国产片特级美女逼逼视频| 欧美变态另类bdsm刘玥| 国产黄频视频在线观看| 啦啦啦中文免费视频观看日本| 国产爽快片一区二区三区| 女性生殖器流出的白浆| 亚洲av.av天堂| 国产精品嫩草影院av在线观看| 欧美+日韩+精品| 精品国产乱码久久久久久小说| av国产久精品久网站免费入址| 国产成人精品一,二区| 交换朋友夫妻互换小说| 男女国产视频网站| 免费看日本二区| 黄色配什么色好看| 国产成人免费观看mmmm| 一本久久精品| av视频免费观看在线观看| 亚洲欧美精品专区久久| 在线观看免费高清a一片| 在线观看人妻少妇| 纵有疾风起免费观看全集完整版| 我要看日韩黄色一级片| 亚洲国产成人一精品久久久| 国产成人freesex在线| 国产午夜精品久久久久久一区二区三区| 嫩草影院新地址| av黄色大香蕉| 成人毛片60女人毛片免费| 欧美少妇被猛烈插入视频| 嘟嘟电影网在线观看| 国产av精品麻豆| 夫妻午夜视频| 中文字幕亚洲精品专区| 亚洲电影在线观看av| 91久久精品国产一区二区三区| 91精品国产九色| 精品国产露脸久久av麻豆| 国产亚洲午夜精品一区二区久久| 亚洲av综合色区一区| 大香蕉久久网| 中文精品一卡2卡3卡4更新| 亚洲成人一二三区av| 国产男女超爽视频在线观看| 国产亚洲欧美精品永久| 免费看日本二区| 中国三级夫妇交换| av播播在线观看一区| 在线观看国产h片| 精品亚洲乱码少妇综合久久| 肉色欧美久久久久久久蜜桃| 亚洲精品国产成人久久av| 久久久久久久大尺度免费视频| 国产在线视频一区二区| 一本色道久久久久久精品综合| 日韩伦理黄色片| 久久久久精品性色| 亚洲欧美成人精品一区二区| 99热国产这里只有精品6| 色婷婷久久久亚洲欧美| 欧美高清成人免费视频www| 熟妇人妻不卡中文字幕| 性高湖久久久久久久久免费观看| xxx大片免费视频| 精品亚洲乱码少妇综合久久| 色吧在线观看| 日日啪夜夜爽| 精品久久久久久久久亚洲| 一级黄片播放器| 国产精品一区二区性色av| 国产成人免费无遮挡视频| 亚洲无线观看免费| 色哟哟·www| 久久精品国产亚洲av天美| 久久免费观看电影| 日韩免费高清中文字幕av| 99国产精品免费福利视频| 一级片'在线观看视频| 亚洲国产精品国产精品| 这个男人来自地球电影免费观看 | 成人午夜精彩视频在线观看| 最近最新中文字幕免费大全7| 如何舔出高潮| 成人毛片a级毛片在线播放| 一区二区三区四区激情视频| 另类亚洲欧美激情| 丰满人妻一区二区三区视频av| 狠狠精品人妻久久久久久综合| 欧美精品一区二区大全| 91精品伊人久久大香线蕉| 尾随美女入室| 日韩中文字幕视频在线看片| 国产精品久久久久久精品古装| 亚洲成色77777| 成人漫画全彩无遮挡| 国产一区有黄有色的免费视频| 又大又黄又爽视频免费| 建设人人有责人人尽责人人享有的| 国产成人精品无人区| 最后的刺客免费高清国语| 日韩,欧美,国产一区二区三区| 少妇丰满av| 欧美丝袜亚洲另类| 天天躁夜夜躁狠狠久久av| 国产精品久久久久成人av| 三级经典国产精品| a级毛片免费高清观看在线播放| 黑人巨大精品欧美一区二区蜜桃 | 亚洲四区av| 女人久久www免费人成看片| 青春草亚洲视频在线观看| 一级爰片在线观看| 亚洲综合色惰| 特大巨黑吊av在线直播| 久久久亚洲精品成人影院| 一区二区三区乱码不卡18| 国产伦理片在线播放av一区| 2018国产大陆天天弄谢| 丁香六月天网| 国产精品秋霞免费鲁丝片| 91精品国产九色| 国产精品成人在线| 国产美女午夜福利| 国产成人精品久久久久久| 一区二区三区乱码不卡18| 国产黄片美女视频| 嫩草影院新地址| 日韩成人伦理影院| 国产亚洲最大av| 人妻 亚洲 视频| 啦啦啦在线观看免费高清www| 久久精品国产a三级三级三级| 久久精品国产自在天天线| 热re99久久精品国产66热6| 国产成人精品久久久久久| 卡戴珊不雅视频在线播放| 久久av网站| 成人特级av手机在线观看| 亚洲欧洲日产国产| 国产黄色视频一区二区在线观看| 精品国产露脸久久av麻豆| 啦啦啦视频在线资源免费观看| 日日摸夜夜添夜夜添av毛片| 搡老乐熟女国产| 免费看不卡的av| 亚洲精品国产成人久久av| 一级毛片aaaaaa免费看小| 亚洲精品色激情综合| 涩涩av久久男人的天堂| 99热网站在线观看| 高清在线视频一区二区三区| 精品国产一区二区久久| 欧美成人午夜免费资源| 最新中文字幕久久久久| 中文字幕制服av| 欧美bdsm另类| 高清毛片免费看| 七月丁香在线播放| 日本91视频免费播放| 99热国产这里只有精品6| 日日啪夜夜爽| 中文精品一卡2卡3卡4更新| 在线观看免费视频网站a站| 午夜久久久在线观看| 中文字幕久久专区| 亚洲精品色激情综合| 中文字幕人妻熟人妻熟丝袜美| 一区二区三区精品91| 三级国产精品片| 天天躁夜夜躁狠狠久久av| 曰老女人黄片| 在线看a的网站| 午夜91福利影院| 国产精品一区二区在线观看99| 少妇猛男粗大的猛烈进出视频| 黄色一级大片看看| 91久久精品国产一区二区成人| 国产成人精品福利久久| 国产免费一区二区三区四区乱码| 日本-黄色视频高清免费观看| 日韩免费高清中文字幕av| 最近中文字幕2019免费版| 丰满少妇做爰视频| 你懂的网址亚洲精品在线观看| 99九九在线精品视频 | 午夜老司机福利剧场| 午夜福利网站1000一区二区三区| 色视频www国产| 日韩制服骚丝袜av| 久久国产精品男人的天堂亚洲 | 欧美精品一区二区大全| 在线免费观看不下载黄p国产| 26uuu在线亚洲综合色| 国产日韩欧美亚洲二区| 国产精品欧美亚洲77777| 国产探花极品一区二区| 国产精品久久久久成人av| 日日啪夜夜撸| 蜜桃久久精品国产亚洲av| 自拍欧美九色日韩亚洲蝌蚪91 | 一级毛片电影观看| 一个人看视频在线观看www免费| 国产免费视频播放在线视频| 老司机影院毛片| 大香蕉97超碰在线| 男女边摸边吃奶| 免费看日本二区| 精华霜和精华液先用哪个| 久久久久久久久久久丰满|