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

    青藏高原東南緣及鄰區(qū)地殼介質(zhì)散射強(qiáng)度變化特征

    2016-07-28 09:30:55范小平王俊菲楊從杰李清河何奕成畢雪梅
    地球物理學(xué)報(bào) 2016年7期
    關(guān)鍵詞:波場(chǎng)塊體斷裂帶

    范小平, 王俊菲, 楊從杰, 李清河, 何奕成, 畢雪梅

    江蘇省地震局, 南京 210014

    ?

    青藏高原東南緣及鄰區(qū)地殼介質(zhì)散射強(qiáng)度變化特征

    范小平, 王俊菲, 楊從杰, 李清河, 何奕成, 畢雪梅

    江蘇省地震局, 南京210014

    摘要基于遠(yuǎn)震擾動(dòng)場(chǎng)方法,利用47個(gè)數(shù)字地震觀測(cè)臺(tái)記錄的10次遠(yuǎn)震記錄對(duì)青藏高原東南緣及鄰區(qū)地殼介質(zhì)散射強(qiáng)度進(jìn)行了研究.結(jié)果表明研究區(qū)內(nèi)地殼介質(zhì)橫向非均勻性強(qiáng)烈,揚(yáng)子塊體地殼介質(zhì)散射強(qiáng)度與青藏塊體存在較為明顯的差異;地殼介質(zhì)散射強(qiáng)度梯度帶與主要斷裂帶分布趨勢(shì)一致,龍門山斷裂帶、安寧河斷裂帶與地殼介質(zhì)散射強(qiáng)度梯度帶具有較強(qiáng)的相關(guān)性,散射強(qiáng)度高值區(qū)偏向青藏塊體一側(cè);強(qiáng)震震中沿散射強(qiáng)度梯度帶分布,其震中偏向散射強(qiáng)度高值區(qū)一側(cè);散射強(qiáng)度呈現(xiàn)縱向非均勻性特征,下地殼介質(zhì)散射強(qiáng)度高于上地殼.研究區(qū)強(qiáng)、弱散射強(qiáng)度的空間分布可能與區(qū)內(nèi)強(qiáng)烈的地質(zhì)構(gòu)造運(yùn)動(dòng)、頻繁的地震及火山活動(dòng)所引起深部物質(zhì)破碎、熔融及其沿?cái)嗔褞У纳嫌颗c運(yùn)移有關(guān).

    關(guān)鍵詞青藏高原東南緣; 散射強(qiáng)度; 地震波散射; 遠(yuǎn)震

    1引言

    觀測(cè)表明,地球介質(zhì)中非均勻體可引起地震波的散射,這些散射攜帶了豐富的介質(zhì)結(jié)構(gòu)信息,因此可以用散射波來(lái)研究介質(zhì)非均勻性(Aki,1973;Wu and Aki,1988; Flatté and Wu,1988; Korn,1990, 1993, 1997).根據(jù)地震波波長(zhǎng)與非均勻體大小關(guān)系,介質(zhì)非均勻性分為大尺度非均勻性和小尺度非均勻性(Sato,1984, 1989; Sato and Fehler, 1998;Saito et al.,2002).大尺度非均勻性可以用傳統(tǒng)的介質(zhì)結(jié)構(gòu)參數(shù)來(lái)描述,如速度、電導(dǎo)率、密度等參數(shù)(Sato,1998),而小尺度非均勻性則采用速度擾動(dòng)率(σ)及非均勻體相關(guān)長(zhǎng)度(a)來(lái)描述(Sato,1984).雖σ和a并不解耦(Sato,1984;Sato and Fehler,1989;Saito et al.,2002),但σ2a卻與地震波散射強(qiáng)度呈線性相關(guān)(Ritter et al.,1998; Ritter and Rothert, 2000),其值反映了單位尺度內(nèi)介質(zhì)速度擾動(dòng)率變化,描述了介質(zhì)非均勻性程度(Sato,1984).遠(yuǎn)震擾動(dòng)場(chǎng)法則可以實(shí)現(xiàn)對(duì)參數(shù)σ2a進(jìn)行定量化估計(jì)(Ritter et al.,1998; Ritter and Rothert, 2000).

    對(duì)于同一遠(yuǎn)震事件,區(qū)域觀測(cè)臺(tái)記錄的P波波場(chǎng)包含了地震波傳播路徑上介質(zhì)結(jié)構(gòu)的信息,除各觀測(cè)臺(tái)站的下方外,地震波傳播路徑基本一致.相同的激發(fā)震源、基本相同的地震波傳播路徑,各臺(tái)站觀測(cè)到地震波場(chǎng)的差異主要反映了臺(tái)站下方介質(zhì)結(jié)構(gòu)的差異(Ritter et al.,1998; Ritter and Rothert, 2000).將地震波場(chǎng)分解為穩(wěn)定場(chǎng)和散射場(chǎng),穩(wěn)定場(chǎng)反映了地震傳播路徑上介質(zhì)結(jié)構(gòu)的“共性”信息,散射場(chǎng)則攜帶了觀測(cè)臺(tái)站下方介質(zhì)結(jié)構(gòu)的“個(gè)性”信息,通過(guò)分析散射場(chǎng)及穩(wěn)定場(chǎng)之間的關(guān)系,便可獲取觀測(cè)臺(tái)站下方地殼介質(zhì)結(jié)構(gòu)信息,這種方法稱為遠(yuǎn)震擾動(dòng)場(chǎng)法(Ritter et al.,1998; Ritter and Rothert, 2000).

    遠(yuǎn)震擾動(dòng)場(chǎng)法最初廣泛應(yīng)用于臺(tái)陣研究中,現(xiàn)逐漸被拓展應(yīng)用介質(zhì)非均勻性研究中.Mach利用LASA數(shù)據(jù)首次對(duì)遠(yuǎn)震P波波場(chǎng)進(jìn)行了系統(tǒng)性研究(Mach,1969),并根據(jù)LASA臺(tái)陣記錄與其它子臺(tái)陣記錄的遠(yuǎn)震P波波場(chǎng)振幅信息非相關(guān)性,發(fā)現(xiàn)了兩個(gè)臺(tái)陣之間地殼介質(zhì)結(jié)構(gòu)存在間斷面.Aki發(fā)現(xiàn)在LASA臺(tái)陣下方非均勻介質(zhì)層厚約60km,且其速度擾動(dòng)率約4%左右(Aki,1973).Langston(Langston,1989)、Hock(Hock et al.,2004)及Wager(Wagner and Langston,1992)等采用數(shù)值模擬方法驗(yàn)證了遠(yuǎn)震擾動(dòng)場(chǎng)方法研究介質(zhì)非均勻性的有效性.Ritter等(1998), Ritter和Rothert(2000)運(yùn)用遠(yuǎn)震擾動(dòng)場(chǎng)方法對(duì)法國(guó)Massif Central地區(qū)地殼上地幔介質(zhì)結(jié)構(gòu)參數(shù)進(jìn)行了研究,獲取該區(qū)地殼上地幔介質(zhì)速度擾動(dòng)率、非均勻體尺度以及地震波散射強(qiáng)度參數(shù).Korn(1990,1993)對(duì)中歐和北歐地區(qū)巖石圈介質(zhì)非均勻程度進(jìn)行了研究.

    青藏高原東南緣處于印度板塊與歐亞板塊碰撞接觸的前沿地帶,是青藏高原與揚(yáng)子地塊擠壓拼接的交匯部分之一,因其特殊的地理位置及其復(fù)雜的、漫長(zhǎng)的劇烈構(gòu)造運(yùn)動(dòng),致使青藏高原東南緣在研究印度板塊與歐亞板塊碰撞過(guò)程、認(rèn)識(shí)青藏高原與揚(yáng)子克拉通碰撞接觸關(guān)系及深部動(dòng)力過(guò)程等地球動(dòng)力學(xué)方面具有獨(dú)特的優(yōu)勢(shì)(鄧起東等,1994).國(guó)內(nèi)外許多學(xué)者對(duì)青藏高原東南緣及周邊地區(qū)深部介質(zhì)結(jié)構(gòu)特征進(jìn)行了研究,如深部地球物理勘探(嘉世旭和張先康,2008,盧占武等,2009)、三維速度結(jié)構(gòu)(劉建華等,1989;王椿鏞等,2002,2003a,2003b;劉啟元等,2009; 吳建平等,2009;雷建設(shè)等,2009;郭飚等,2009;楊海燕等,2009;李志偉等,2011;李永華等,2009)、Q值成像(周連慶等,2008)、橫波分裂與衰減(馬宏生等,2007;蘇有錦,2006;石玉濤等,2009)、深部電性結(jié)構(gòu)(王緒本等,2009)等,結(jié)果均表明該區(qū)域地殼上地幔存在強(qiáng)烈的介質(zhì)非均勻性.但區(qū)域地殼介質(zhì)非均勻性強(qiáng)弱程度及強(qiáng)弱非均勻體空間分布特征如何,到目前為止,還沒(méi)有學(xué)者對(duì)此問(wèn)題進(jìn)行量化研究.基于此,本文運(yùn)用遠(yuǎn)震擾動(dòng)場(chǎng)方法對(duì)青藏高原東南緣及鄰區(qū)地殼介質(zhì)散射強(qiáng)度變化特征進(jìn)行研究,分析其與深部斷裂構(gòu)造、中強(qiáng)震孕育之間的關(guān)系,并探討其可能蘊(yùn)含的地球動(dòng)力學(xué)意義.

    2原理與方法

    受非均勻體的作用,地震波穿過(guò)非均勻介質(zhì)層時(shí)會(huì)引起地震波的散射.散射地震波場(chǎng)可分為穩(wěn)定波場(chǎng)和散射波場(chǎng),而非均勻介質(zhì)可認(rèn)為由均勻介質(zhì)和非均勻介質(zhì)疊加構(gòu)成,因此可以認(rèn)為穩(wěn)定波場(chǎng)反映了均勻介質(zhì)的結(jié)構(gòu)信息,散射波場(chǎng)則包含了非均勻介質(zhì)的結(jié)構(gòu)信息.非均勻介質(zhì)中任意一點(diǎn)r處,地震波場(chǎng)(Ritter et al.,1998; Ritter and Rothert, 2000)可以表示為

    (1)

    其中Ut(r,t)、〈U(r,t)〉和Uf(r,t)分別代表在空間r(x,y,z)處、t時(shí)刻的總波場(chǎng)、穩(wěn)定波場(chǎng)和散射波場(chǎng),三角括號(hào)代表對(duì)波場(chǎng)進(jìn)行空間(統(tǒng)計(jì))平均.那么,散射波場(chǎng)的擾動(dòng)強(qiáng)度可用無(wú)量綱參數(shù)ε來(lái)描述,

    (2)

    如果不考慮介質(zhì)粘彈性對(duì)地震波的衰減,地震波的能量損失主要由非均勻介質(zhì)中非均勻體的散射所引起,那么根據(jù)能量守恒原理,地震波場(chǎng)強(qiáng)度可表示為

    (3)

    (4)

    由(2)—(4)式可以看出,若〈ε2〉?1,也就意味著穩(wěn)定波場(chǎng)強(qiáng)度遠(yuǎn)遠(yuǎn)大于散射波場(chǎng)強(qiáng)度,則表明介質(zhì)非均勻性程度相對(duì)較弱;若〈ε2〉?1,意味著散射場(chǎng)強(qiáng)度遠(yuǎn)大于穩(wěn)定波場(chǎng)強(qiáng)度,則表明介質(zhì)非均勻性程度相對(duì)較強(qiáng).因此,方程(4)可用來(lái)研究地殼上地幔介質(zhì)非均勻性強(qiáng)弱程度.

    對(duì)于高斯型或指數(shù)型隨機(jī)介質(zhì),〈ε2〉與非均勻介質(zhì)結(jié)構(gòu)參數(shù)(c、σ、a和L)之間的關(guān)系(Ritteretal.,1998;RitterandRothert, 2000)可表示為

    (5)

    其中c為背景介質(zhì)速度,a為非均勻介質(zhì)相關(guān)長(zhǎng)度,L為非均勻介質(zhì)中地震波傳播路徑,σ為介質(zhì)速度擾動(dòng)率.由(5)式可以看出,ln(〈ε2〉+1)與頻率f呈正相關(guān),通過(guò)最小二乘法在頻率域?qū)n(〈ε2〉+1)進(jìn)行擬合,擬合系數(shù)γ則為地震波的散射強(qiáng)度.對(duì)于某一特定的研究區(qū)域,散射層的厚度L及介質(zhì)P波傳播速度可根據(jù)寬角反射、折射資料或其它資料獲取,在計(jì)算出散射層結(jié)構(gòu)參數(shù)L/c2后,通過(guò)(5)式可求取σ2a的取值及變化范圍.σ2a反映了單位尺度內(nèi)介質(zhì)速度擾動(dòng)率的變化.在某一特定的非均勻體尺度a下,σ越大則表明介質(zhì)非均勻性程度越強(qiáng),σ越小則反映了介質(zhì)非均勻性程度相對(duì)較弱.

    3研究區(qū)構(gòu)造背景

    青藏高原東南緣及鄰近地區(qū),其主要構(gòu)造單元為松潘─甘孜褶皺帶和揚(yáng)子準(zhǔn)地臺(tái).圖1為青藏高原東緣及鄰區(qū)構(gòu)造示意圖.由圖1可以看出,北東向的龍門山斷裂帶(F1)是一條重要的構(gòu)造分界線,其西側(cè)為松潘─甘孜褶皺帶的一套淺變質(zhì)的古生界及淺變質(zhì)的中生界三疊系地層,經(jīng)逆沖或推覆到東側(cè)屬于揚(yáng)子準(zhǔn)地臺(tái)的古生界及中生界穩(wěn)定淺海相和陸相沉積地層上(鄧起東等,1994; 李志偉等,2011).揚(yáng)子準(zhǔn)地臺(tái)晚古生代以來(lái)沉積環(huán)境比較穩(wěn)定,具有較厚未變質(zhì)的沉積蓋層,在燕山運(yùn)動(dòng)和喜馬拉雅運(yùn)動(dòng)期間才出現(xiàn)褶皺運(yùn)動(dòng).鮮水河斷裂帶(F2)是第四紀(jì)強(qiáng)烈活動(dòng)的左旋走滑斷裂帶,近100年來(lái),沿該斷裂發(fā)生了一系列的強(qiáng)震,安寧河斷裂(F3)是形成于前震旦紀(jì)一條區(qū)域性深大斷裂,該斷裂特點(diǎn)是深度大,傾角大,有多期巖漿的侵入和噴發(fā),向南延伸到云南境內(nèi),并與具有強(qiáng)烈南北向新構(gòu)造活動(dòng)的小江斷裂相連(郭飚等,2009;楊海燕等,2009;李志偉等,2011;李永華等,2009).研究區(qū)內(nèi)發(fā)育了一系列次級(jí)活動(dòng)斷裂,如龍泉山斷裂(f1)、馬邊嚴(yán)津斷裂帶(f2)、蓮峰斷裂(f3)、玉龍溪斷裂(f4)及木拉斷裂(f5).

    圖1 青藏高原東南緣及鄰區(qū)主要地質(zhì)構(gòu)造及臺(tái)站位置圖 F1:龍門山斷裂帶, F2:鮮水河斷裂帶, F3:安寧河斷裂帶, f1:龍泉山斷裂, f2:馬邊嚴(yán)津斷裂帶, f3:蓮峰斷裂, f4:玉龍溪斷裂, f5:木拉斷裂.Fig.1 Geological structures and seismic stations in the southeast margin of Tibetan plateau and adjacent regions F1: Longmengshan fault zone, F2: Xianshuihe fault zone, F3: Anninghe fault zone, f1: Longquanshan fault, f2: Mabian-yanjin fault, f3: Lianfeng fault, f4: Yulongxi fault; f5: Mula fault.

    4數(shù)據(jù)收集

    本文收集了2007年1月─2010年12月之間四川省地震臺(tái)網(wǎng)47個(gè)地震臺(tái)站所記錄(鄭秀芬等,2009)的地震波形,這些記錄必須滿足以下幾個(gè)條件:

    (1) 震中距大于70°以便獲取垂直入射的P波波場(chǎng);

    (2) 震源深度大于70km,以便避開(kāi)PP震相;

    (3)P波初至清楚、簡(jiǎn)單,47個(gè)觀測(cè)臺(tái)站均記錄到該事件.

    同時(shí)為了避免其它震相遠(yuǎn)震P波波場(chǎng)的干擾,我們對(duì)每個(gè)地震事件進(jìn)行速度譜分析,將“速度振幅─時(shí)間域”記錄轉(zhuǎn)換到“慢度─時(shí)間域”,其慢度變化范圍為3.5s/ °~7.5s/°,步長(zhǎng)為0.2s/°;并根據(jù)“慢度─時(shí)間域”震相分布范圍,剔除P波尾波20秒內(nèi)有其它震相的地震事件.挑選后遠(yuǎn)震地震事件參數(shù)見(jiàn)表1.

    5數(shù)據(jù)處理及分析

    5.1數(shù)據(jù)處理

    為獲取研究區(qū)地殼介質(zhì)地震散射強(qiáng)度,需要對(duì)遠(yuǎn)震記錄進(jìn)行下面幾個(gè)步驟處理:

    (1) 首先對(duì)原始波形進(jìn)行0.1~6.0 Hz的帶通濾波,去除低頻和相對(duì)高頻干擾;

    表1 遠(yuǎn)震事件參數(shù)表

    (2) 對(duì)遠(yuǎn)震P波震相進(jìn)行P波初至對(duì)齊,將對(duì)齊后的P波波場(chǎng)進(jìn)行疊加,獲取穩(wěn)定場(chǎng);

    (3) 用各觀測(cè)臺(tái)原始記錄減去穩(wěn)定波場(chǎng),獲取散射波場(chǎng);

    (4) 對(duì)擾動(dòng)場(chǎng)和穩(wěn)定場(chǎng)進(jìn)行傅氏變換;

    (5) 用擾動(dòng)場(chǎng)譜除以穩(wěn)定場(chǎng)譜,獲取散射波場(chǎng)強(qiáng)度;

    (6) 在對(duì)數(shù)空間對(duì)散射波場(chǎng)能量上升段進(jìn)行最小二乘擬合,獲取反映地震波散射強(qiáng)度的參數(shù)γ;

    (7) 根據(jù)散射層厚度(L)及相應(yīng)層P波速度(C),計(jì)算觀測(cè)臺(tái)L/c2;

    (8) 根據(jù)公式(5)可獲取σ2a值.

    分別對(duì)47個(gè)觀測(cè)臺(tái)站重復(fù)執(zhí)行上述8個(gè)步驟,便可獲取一個(gè)地震事件的計(jì)算結(jié)果.為了增加結(jié)果的可靠性,減少隨機(jī)誤差,對(duì)表1中10個(gè)地震事件分別進(jìn)行計(jì)算,并將每個(gè)臺(tái)站10個(gè)地震事件計(jì)算結(jié)果平均,作為該臺(tái)站的最終結(jié)果,這樣便可獲取σ2a值的空間分布特征.

    為了獲取研究區(qū)47個(gè)觀測(cè)臺(tái)站下方地殼內(nèi)介質(zhì)地震波散射強(qiáng)度,建立雙層(上地殼和下地殼)散射模型,綜合區(qū)內(nèi)寬角反射/折射(王椿鏞等,2002,2003a,2003b)、遠(yuǎn)震接收函數(shù)(李永華等,2009)、天然地震資料統(tǒng)計(jì)分析(趙珠和張潤(rùn)生,1987)研究成果確定各觀測(cè)臺(tái)站的散射層的厚度(L)及彈性P波速度(C),各參數(shù)值見(jiàn)表2.

    表2 各觀測(cè)臺(tái)站地殼速度結(jié)構(gòu)及介質(zhì)地震散射強(qiáng)度

    5.2數(shù)據(jù)分析

    圖2a為研究區(qū)典型穩(wěn)定波場(chǎng)(a)和散射波場(chǎng)(b)記錄圖.圖2a中第一行為穩(wěn)定波場(chǎng)記錄,其余各行分別為各觀測(cè)臺(tái)站記錄的原始地震波波形.圖2b則為相應(yīng)各觀測(cè)臺(tái)記錄的散射波波場(chǎng).由圖2a可以看出,穩(wěn)定波場(chǎng)(stack)信噪比明顯高于各觀測(cè)臺(tái)記錄波形的信噪比,散射引起的高頻信息得到了有效壓制,主要震相得以突出顯示;由圖2b可以看出,無(wú)論是同一區(qū)域,還是不同區(qū)域,幾乎沒(méi)有相同的散射波場(chǎng),這表明各觀測(cè)臺(tái)站下方介質(zhì)結(jié)構(gòu)存在較為明顯的差異,介質(zhì)結(jié)構(gòu)不同,其對(duì)地震波作用強(qiáng)度則不同.如AXI臺(tái)與ZJG臺(tái)雖相距不足40 km(圖1),但兩臺(tái)記錄的散射波場(chǎng)無(wú)論是頻率還是波形形態(tài)均有明顯差異(圖2a、圖2b).

    圖3為研究區(qū)穩(wěn)定波場(chǎng)(紫色線)與散射波場(chǎng)(藍(lán)色線)頻譜圖.由圖3可以看出,穩(wěn)定波場(chǎng)與散射波場(chǎng)能量強(qiáng)度不同頻段存在明顯差異.在0.1~1.1 Hz區(qū)間,穩(wěn)定波場(chǎng)強(qiáng)度明顯大于散射波場(chǎng)強(qiáng)度;而在1.1~6.0 Hz頻段內(nèi),散射波場(chǎng)強(qiáng)度又明顯強(qiáng)于穩(wěn)定波場(chǎng)強(qiáng)度.

    圖3 研究區(qū)穩(wěn)定波場(chǎng)(紫色線)與 散射波場(chǎng)(藍(lán)色線)頻譜圖Fig.3 Frequency spectra of stable filed (purple line) and scattering wave filed (blue line) in study area

    圖4為研究區(qū)穩(wěn)定波場(chǎng)與典型散射波場(chǎng)頻譜圖.由圖4可以看出,榮縣臺(tái)(HMS)和汶川臺(tái)(WCH)散射波場(chǎng)能量與穩(wěn)定波場(chǎng)能量存在明顯差異.在0.1~0.75 Hz區(qū)間,穩(wěn)定波場(chǎng)能量均強(qiáng)于HMS臺(tái)和WCH臺(tái)散射波場(chǎng)能量,但WCH臺(tái)從0.75 Hz開(kāi)始散射波場(chǎng)能量強(qiáng)于穩(wěn)定波場(chǎng)能量,而HMS臺(tái)則從1.0 Hz開(kāi)始其能量開(kāi)始優(yōu)于穩(wěn)定波場(chǎng)能量強(qiáng)度;就散射波場(chǎng)能量強(qiáng)度而言,從1.2 Hz開(kāi)始,WCH臺(tái)散射波場(chǎng)能量強(qiáng)度強(qiáng)于HMS臺(tái)散射波場(chǎng)強(qiáng)度.由不同臺(tái)站記錄散射波場(chǎng)的差異性可以看出,雖具有相同的激發(fā)震源、幾乎相同的地震波傳播路徑,但由于臺(tái)站下方物質(zhì)結(jié)構(gòu)的差異,導(dǎo)致了不同的散射波場(chǎng)特征.WCH臺(tái)位于龍門山斷裂帶上,而HMS臺(tái)則處四川盆地邊緣,地質(zhì)構(gòu)造運(yùn)動(dòng)史及深部物質(zhì)的差異性,直接導(dǎo)致了二者散射波場(chǎng)及其波場(chǎng)強(qiáng)度轉(zhuǎn)換頻率的差異性.

    圖4 研究區(qū)穩(wěn)定波場(chǎng)(紫色線)與典型散射波場(chǎng) (HWS臺(tái):綠色線;WCH臺(tái):藍(lán)色線)頻譜圖Fig.4 Frequency spectra of stable field (purple line) and typical scattering wave fields (HWS :green line; WCH: blue line) in study area

    圖5為研究區(qū)典型散射波場(chǎng)強(qiáng)度圖.圖5a為小金臺(tái)(XJI)、峨眉山臺(tái)(EMS)及瀘沽湖(LGH)典型散射波場(chǎng)強(qiáng)度,圖5b分別為XJI、EMS及LGH臺(tái)散射系數(shù)擬合圖.由圖5可以看出,因臺(tái)站分布區(qū)域不同,散射波場(chǎng)強(qiáng)度不同,散射波場(chǎng)與穩(wěn)定波場(chǎng)能量也有差異.XJI臺(tái)、EMS臺(tái)及LGH臺(tái)穩(wěn)定波場(chǎng)與散射波場(chǎng)強(qiáng)度轉(zhuǎn)換頻率點(diǎn)分別為1.57、1.92、2.89 Hz;XJI臺(tái)、EMS臺(tái)及LGH臺(tái)散射系數(shù)分別為0.179、0.162及0.061.各觀測(cè)臺(tái)地震波散射系數(shù)的不同,體現(xiàn)了臺(tái)站下方介質(zhì)結(jié)構(gòu)對(duì)地震波的作用程度不同,進(jìn)一步說(shuō)明了各臺(tái)站下方介質(zhì)結(jié)構(gòu)的差異性,也說(shuō)明了研究區(qū)介質(zhì)非均勻性程度的差異性.

    圖5 研究區(qū)典型散射波場(chǎng)強(qiáng)度圖(a:散射波場(chǎng)強(qiáng)度;b:散射系數(shù)擬合)Fig.5 Typical scattering intensity of fluctuation wave field in study area (a: scattering intensity of scattering wave filed, b: fitting of scattering strength).

    6結(jié)果及分析

    表2為研究區(qū)內(nèi)各觀測(cè)臺(tái)站地殼速度結(jié)構(gòu)及介質(zhì)地震散射強(qiáng)度值.圖6為研究區(qū)地殼介質(zhì)地震波散射強(qiáng)度空間分布圖.圖6中不同顏色代表地殼介質(zhì)散射強(qiáng)度不同,藍(lán)色區(qū)域代表強(qiáng)散射區(qū),紅色區(qū)域代表相對(duì)較弱散射區(qū);介質(zhì)散射性越強(qiáng),代表介質(zhì)非均勻性程度較高,介質(zhì)較為破碎,而散射性越低,則表明介質(zhì)非均勻性程度較低,介質(zhì)較為均勻、完整.

    圖6 上地殼介質(zhì)地震波散射強(qiáng)度分布圖Fig.6 Map of seismic scattering strength of the upper crust medium in study area

    研究區(qū)內(nèi)地殼介質(zhì)散射強(qiáng)度呈明顯的橫向不均勻分布特征.上地殼(圖6),馬爾康(MEK)、小金(XJI)、汶川(WCH)及平武(PWU)地區(qū)同屬于松潘─甘孜褶皺帶,但馬爾康地區(qū)介質(zhì)散射強(qiáng)度明顯弱于小金、汶川及平武地區(qū);九龍(JLO)、瀘沽湖(LGH)、鹽源(YYU)及雅江(YJI)地區(qū)同屬于川滇菱形塊體,但鹽源及雅江地區(qū)介質(zhì)散射強(qiáng)度明顯高于九龍和瀘沽湖地區(qū);雖同屬于川中塊體(揚(yáng)子塊體西緣),但仁壽(YGD)、容縣(HMS)地區(qū)地殼介質(zhì)散射強(qiáng)度則強(qiáng)于塊體內(nèi)其它地區(qū).相比較而言,川中塊體上地殼介質(zhì)地震散射強(qiáng)度整體偏小.下地殼(圖7),沿龍門山斷裂帶仍為散射強(qiáng)度高值分布區(qū),分別在WCH和PWU地區(qū)存在兩個(gè)散射強(qiáng)度高值中心;沿汶川、都江堰、成都、自貢一帶存在一條北西向散射強(qiáng)度高值分布帶;沿安寧河斷裂帶散射強(qiáng)度呈相對(duì)高值分布,而其兩側(cè)散射強(qiáng)度則呈相對(duì)低值分布;散射強(qiáng)度低值區(qū)分布在黑水(HSH)、小金(XJI)及九龍(JLO)以西地區(qū).

    圖7 下地殼介質(zhì)地震波散射強(qiáng)度分布圖Fig.7 Map of seismic scattering strength of the lower crust medium in study area

    地殼介質(zhì)散射強(qiáng)度梯度帶基本上沿主要斷裂帶分布.上地殼(圖6),龍門山斷裂帶(F1)兩側(cè)介質(zhì)散射強(qiáng)度呈明顯差異,龍門山斷裂帶東側(cè)區(qū)域?yàn)樯⑸鋸?qiáng)度低值區(qū),以地震弱散射為主,而斷裂帶西側(cè)區(qū)域?yàn)樯⑸鋸?qiáng)度高值區(qū),以地震強(qiáng)散射為主.安寧河斷裂帶(F3)內(nèi)與其兩側(cè)地殼介質(zhì)散射強(qiáng)度也呈現(xiàn)出較強(qiáng)的差異性,安寧河斷裂帶東側(cè)及西側(cè)區(qū)域?yàn)樯⑸鋸?qiáng)度低值區(qū),以地震波弱散射為主,但安寧河斷裂帶(F3)內(nèi)部為散射強(qiáng)度高值區(qū),以地震波強(qiáng)散射為主.鮮水河斷裂帶(F2)兩側(cè)地殼介質(zhì)散射強(qiáng)度也有差異性,鮮水河斷裂帶兩側(cè)散射強(qiáng)度梯度明顯小于龍門山斷裂帶和安寧河斷裂帶散射強(qiáng)度梯度.鮮水河斷裂帶的東北側(cè)地震散射強(qiáng)度相對(duì)較強(qiáng),而西南側(cè)地震散射強(qiáng)度相對(duì)較弱.沿鮮水河斷裂帶地震散射強(qiáng)度似乎有分段性特征,道浮(DFU)東南地區(qū),地震散射系數(shù)相對(duì)較強(qiáng),而其北西一側(cè),地震散射系數(shù)相對(duì)較弱.下地殼(圖7),地殼介質(zhì)散射強(qiáng)度在不同區(qū)域其散射強(qiáng)度有明顯差異.沿龍門山斷裂帶仍表現(xiàn)出較強(qiáng)的散射強(qiáng)度梯度帶,龍門山斷裂帶西側(cè)散射強(qiáng)度衰減速度相對(duì)比東側(cè)高;沿安寧河斷裂帶仍具有較為明顯的散射強(qiáng)度梯度帶,散射強(qiáng)度高值區(qū)沿安寧河斷裂帶分布,其兩側(cè)則呈相對(duì)低值分布;沿鮮水河斷裂帶仍為散射強(qiáng)度梯度帶,但梯度值相對(duì)偏小.

    地殼介質(zhì)散射強(qiáng)度梯度帶基本上為強(qiáng)震孕育區(qū).由圖6和圖7可以看出,研究區(qū)內(nèi)6級(jí)以上地震(公元624年—2014年)震中基本上分布于散射強(qiáng)度梯度帶上,尤其是7級(jí)以上地震震中與散射強(qiáng)度梯度帶的相關(guān)性更強(qiáng),且地震震中偏向散射強(qiáng)度高值一側(cè).在上地殼,汶川8.0級(jí)地震震中位于散射強(qiáng)度梯度帶上(圖6),但梯度值并不明顯,但在下地殼(圖7)震中處于散射強(qiáng)度強(qiáng)梯度帶上,沿北東向偏向散射強(qiáng)度高值一側(cè).位于安寧河斷裂帶南段以及鮮水河斷裂帶西北段的地震震中,在上地殼(圖6)與介質(zhì)散射強(qiáng)度梯度帶具有較強(qiáng)相關(guān)性,而在下地殼,除安寧河斷裂帶南段部分地震震中外,其余地震震中與散射強(qiáng)度梯度帶相關(guān)性相對(duì)較差.從不同深度處散射強(qiáng)度梯度帶與地震震中之間的相關(guān)關(guān)系,大致可以判斷出地震震源孕育深度,汶川地震震中與下地殼介質(zhì)散射強(qiáng)度梯度帶有較強(qiáng)的相關(guān)性,精確定位結(jié)果顯示,汶川8.0級(jí)地震震源深度約20 km左右(劉啟元等,2009;雷建設(shè)等,2009),該深度處于下地殼的頂部區(qū)域.

    下地殼介質(zhì)散射強(qiáng)度均強(qiáng)于上地殼.圖8為上地殼相對(duì)于下地殼散射強(qiáng)度的擾動(dòng).由圖8可以看出,上地殼相對(duì)于下地殼散射強(qiáng)度的擾動(dòng)均為負(fù)擾動(dòng),說(shuō)明上地殼介質(zhì)散射強(qiáng)度均弱于下地殼.相比較而言,川中塊體下地殼散射強(qiáng)度幅值變化較為明顯,最大增幅可達(dá)到82%;松潘─甘孜褶皺帶、川滇菱形塊體下地殼散射強(qiáng)度幅值變化較弱,變化幅度最大為8%左右.從上、下地殼介質(zhì)散射強(qiáng)度變化可以看出,川中塊體與青藏塊體地殼介質(zhì)存在較為明顯的介質(zhì)結(jié)構(gòu)差異,這種差異可能與各構(gòu)造塊體深部介質(zhì)結(jié)構(gòu)、物質(zhì)成份及復(fù)雜的地質(zhì)構(gòu)造歷史有關(guān).川中塊體地殼厚度相對(duì)較薄,平均厚度約40 km左右(趙珠和張潤(rùn)生,1987),而青藏塊體地殼厚度較厚,平均厚度在60 km左右(趙珠和張潤(rùn)生,1987;王椿鏞等,2002,2003a,2003b;),同時(shí)川中塊體、青藏塊體的地殼速度結(jié)構(gòu)也有較為明顯的差異(盧占武等,2009;王椿鏞等,2002,2003a, 2003b).沿龍門山斷裂帶、安寧河斷裂帶呈現(xiàn)出較明顯的散射強(qiáng)度擾動(dòng)梯度帶(圖8),作為青藏塊體與揚(yáng)子塊體的邊界接觸帶,這條梯度帶進(jìn)一步說(shuō)明了青藏塊體與揚(yáng)子塊體深部物質(zhì)結(jié)構(gòu)的差異性,同時(shí)也說(shuō)明二者之間的接觸關(guān)系.

    圖8 上地殼相對(duì)于下地殼散射強(qiáng)度的擾動(dòng)Fig.8 The scattering strength perturbation of the upper crust relative to the lower crust in study area

    7討論與結(jié)論

    7.1本文運(yùn)用遠(yuǎn)震擾動(dòng)場(chǎng)法對(duì)研究區(qū)內(nèi)地殼介質(zhì)地震散射強(qiáng)度進(jìn)行了研究.研究結(jié)果表明地殼介質(zhì)散射強(qiáng)度呈現(xiàn)強(qiáng)烈的橫向非均勻性,揚(yáng)子塊體地殼介質(zhì)散射強(qiáng)度與青藏塊體存在較為明顯的差異;散射強(qiáng)度梯度與主要斷裂構(gòu)造帶分布趨勢(shì)一致;強(qiáng)震震中沿散射強(qiáng)度梯度帶分布,且震中偏向散射強(qiáng)度高值區(qū)一側(cè);散射強(qiáng)度表現(xiàn)出縱向非均勻性特征,下地殼介質(zhì)散射強(qiáng)度強(qiáng)于上地殼,揚(yáng)子塊體介質(zhì)散射強(qiáng)度與青藏塊體存在較為明顯的差異.

    7.2研究區(qū)散射強(qiáng)度高值區(qū)(地震波強(qiáng)散射區(qū))分布范圍與地殼介質(zhì)Q值異常(周連慶等,2008;蘇有錦等,2006)、P波速度異常(郭飚等,2009;楊海燕等,2009;李志偉等,2011)分布范圍具有較強(qiáng)的一致性.沿川滇塊體的東南邊界、龍門山斷裂帶以西區(qū)域均為介質(zhì)高衰減區(qū)域(周連慶等,2008;蘇有錦等,2006)、P波波速異常(郭飚等,2009;楊海燕等,2009;李志偉等,2011)及散射強(qiáng)度高值區(qū).無(wú)論是Q值異常、P波波速異常,還是散射強(qiáng)度異常,三者表現(xiàn)形式不同,但具有較為一致的物理本質(zhì).介質(zhì)非均勻性可引起地震波散射,介質(zhì)非均勻性程度越強(qiáng)其對(duì)地震波散射作用就越強(qiáng),因而散射強(qiáng)度就呈現(xiàn)高值,介質(zhì)非均性程度又與介質(zhì)完整性有關(guān),介質(zhì)完整性越差,則意味著介質(zhì)速度相對(duì)較低;介質(zhì)完整性較差,其對(duì)地震波散射、衰減作用就越強(qiáng),相應(yīng)其Q值就偏低.

    7.3研究區(qū)地殼介質(zhì)散射強(qiáng)度分布與活動(dòng)斷裂、強(qiáng)震活動(dòng)區(qū)域具較強(qiáng)的相關(guān)性,這種相關(guān)性可能與區(qū)內(nèi)地質(zhì)構(gòu)造運(yùn)動(dòng)、地震活動(dòng)以及火山活動(dòng)有關(guān).印度板塊對(duì)青藏高原的俯沖碰撞動(dòng)力下,青藏地塊區(qū)物質(zhì)橫向(東流)溢出,東流物質(zhì)在其東部受川中塊體的抵阻,進(jìn)而向東南方向移動(dòng),形成了川滇菱形塊體.龍門山斷裂帶、安寧河斷裂帶作為青藏塊體、川滇菱形塊體與川中塊體的接觸帶,強(qiáng)烈的地質(zhì)構(gòu)造運(yùn)動(dòng)、頻繁的地震活動(dòng)及火山活動(dòng)引起深部物質(zhì)沿?cái)嗔褞嫌颗c運(yùn)移,造成斷裂帶發(fā)育區(qū)深部介質(zhì)破碎與熔融.復(fù)雜的深部介質(zhì)結(jié)構(gòu)及強(qiáng)烈的介質(zhì)非均勻性,致使沿龍門山斷裂帶、安寧河斷裂帶成為地震散射強(qiáng)度梯度帶,也成為地震孕育帶.

    7.4對(duì)比圖6和圖7可以看出,在上地殼,沿安寧河斷裂帶、并在安寧河斷裂帶、龍門山斷裂帶以及鮮水河斷裂帶交匯區(qū)域,地殼介質(zhì)散射強(qiáng)度呈現(xiàn)高值區(qū),而在下地殼,沿龍門山斷裂帶地殼介質(zhì)散射強(qiáng)度呈現(xiàn)高值區(qū).地殼介質(zhì)散射強(qiáng)度的強(qiáng)弱差異反映了介質(zhì)非均勻性程度,散射強(qiáng)度高值區(qū)意味著強(qiáng)非均勻體的存在,而散射強(qiáng)度低值區(qū),則表明介質(zhì)非均勻性程度較弱.因此,上、下地殼介質(zhì)強(qiáng)、弱散射強(qiáng)度分布的差異性,可能反映了青藏塊體與揚(yáng)子塊體邊界的接觸位置及作用深度.

    研究區(qū)上、下地殼散射強(qiáng)度的差異性可能與區(qū)內(nèi)復(fù)雜的斷裂系統(tǒng)有關(guān).研究表明區(qū)內(nèi)巖石圈存在兩套斷裂系統(tǒng),一是以地殼表層脆性剪切帶為主的淺層斷裂系統(tǒng);另一是以切割莫霍界面或殼幔過(guò)渡帶的韌性剪切帶為主的深部斷裂系統(tǒng),稱為殼幔韌性剪切帶(蔡學(xué)林等,2008).區(qū)內(nèi)殼幔韌性剪切帶屬于擠壓型殼幔韌性剪切帶,在深、淺不同的斷裂系統(tǒng)作用下,造成了復(fù)雜的上、下地殼介質(zhì)結(jié)構(gòu).下地殼沿汶川、都江堰、成都、自貢、瀘州一帶分布的介質(zhì)散射強(qiáng)度高值區(qū),與深部人工地震剖面揭示的都江堰殼幔韌性剪切帶位置較為一致(蔡學(xué)林等,2008).

    致謝中國(guó)地震局地球物理研究所“國(guó)家數(shù)字測(cè)震臺(tái)網(wǎng)數(shù)據(jù)備份中心”為本研究提供地震波形數(shù)據(jù),在此表示感謝!感謝評(píng)審專家提出的寶貴意見(jiàn),使得本文論述更趨完善.

    References

    Aki K. 1973. Scattering of P waves under the Montana Lasa. J. Geophys. Res., 78(8): 1334-1346.

    Flatté S M, Wu R S. 1988. Small-scale structure in the lithosphere and asthenosphere deduced from arrival time and amplitude fluctuations at NORSAR. J. Geophys. Res., 93(B6): 6601-6614.

    Guo B, Liu Q Y, Chen J H, et al. 2009. Teleseismic P-wave tomography of the crust and upper mantle in Longmenshan area, west Sichuan. Chinese Journal of Geophysics (in Chinese), 52(2): 346-355. Hock S, Korn M, Ritter J R R, et al. 2004. Mapping random lithospheric heterogeneities in northern and central Europe. Geophys. J. Int., 157(1): 251-264.

    Jia S X, Zhang X K. 2008. Study on the crust phases of deep seismic sounding experiments and fine crust structures in the northeast margin of Tibetan plateau. Chinese Journal of Geophysics (in Chinese), 51(5): 1431-1443.

    Korn M. 1990. A modified energy flux model for lithospheric scattering of teleseismic body waves. Geophys. J. Int., 102(1): 165-175.Korn M. 1993. Determination of site-dependent scattering Q from P-wave coda analysis with an energy-flux model. Geophys. J. Int., 113(1): 54-72.

    Korn M. 1997. Modelling the teleseismic P coda envelope: depth dependent scattering and deterministic structure. Phys. Earth Planet. Int., 104(1-3): 23-36.

    Langston C A. 1989. Scattering of teleseismic body waves under Pasadena, California. J. Geophys. Res., 94(B2): 1935-1951.

    Lei J S, Zhao D P, Su J R, et al. 2009. Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuan earthquake. Chinese Journal of Geophysics (in Chinese), 52(2): 339-345.Li Y H, Wu Q J, Tian X B, et al. 2009. Crustal structure in the Yunnan region determined by modeling receiver functions. Chinese Journal of Geophysics (in Chinese), 52(1): 67-80.

    Liu J H, Loi F T, Wu H, et al. 1989. Three dimensional velocity images of the crust and upper mantle beneath north-south zone in China. Acta Geophysica Sinica (in Chinese), 32(2): 143-151.Liu Q Y, Li Y, Chen J H, et al. 2009. Wenchuan MS8.0 earthquake: preliminary study of the S-wave velocity structure of the crust and upper mantle. Chinese Journal of Geophysics (in Chinese), 52(2): 309-319.

    Lu Z W, Gao R, Li Q S, et al. 2009. Testing deep seismic reflection profiles across the central uplift of the Qiangtang terrane in the Tibetan Plateau. Chinese Journal of Geophysics (in Chinese), 52(8): 2008-2014.

    Ma H S, Wang S Y, Pei S P, et al. 2007. Q0 tomography of S wave attenuation in Sichuan-Yunnan and adjacent regions. Chinese Journal of Geophysics (in Chinese), 50(2): 465-471.

    Mach H. 1969. Nature of short-period P-wave signal variations at LASA. J. Geophys. Res, 74(12): 3161-3170. Ritter J R R, Shapiro S A, Schechinger B. 1998. Scattering parameters of the lithosphere below the Massif Central, France, from teleseismic wavefield records. Geophys. J. Int., 134(1): 187-198.

    Ritter J R R, Rothert E. 2000. Variations of the lithospheric seismic scattering strength below the Massif Central, France and the Frankonian Jura, SE Germany. Tectonophysics, 328(3-4): 297-305.

    Saito T, Sato H, Ohtake M. 2002. Envelope broadening of spherically outgoing waves in three-dimensional random media having power law spectra. J. Geophys. Res., 107(B5): ESE 3-1-ESE 3-15.

    Sato H. 1984. Attenuation and envelope formation of three-component seismograms of small local earthquakes in randomly inhomogeneous lithosphere. J. Geophys. Res., 89(B2): 1221-1241.

    Sato H. 1989. Broadening of seismogram envelopes in the randomly inhomogeneous lithosphere based on the parabolic approximation: Southeastern Honshu, Japan. J. Geophys. Res., 94(B12): 17735-17747. Sato H, Fehler M C. 1998. Seismic wave propagation and Scattering in the Heterogeneous Earth. New York: Springer, 1-3.

    Shi Y T, Tao Y, Zhao C P, et al. 2009. A study of seismic anisotropy of Wenchuan earthquake sequence. Chinese Journal of Geophysics (in Chinese), 52(2): 398-407.

    Wagner G S, Langston C A. 1992. A numerical investigation of scattering effects for teleseismic plane wave propagation in a heterogeneous layer over a homogeneous half-space. Geophys. J. Int., 110(3): 486-500.

    Wang X B, Zhu Y T, Zhao X K, et al. 2009. Deep conductivity characteristics of the Longmen Shan, Eastern Qinghai-Tibet Plateau. Chinese Journal of Geophysics (in Chinese), 52(2): 564-571.

    Wu J P, Huang Y, Zhang T Z, et al. 2009. Aftershock distribution of the MS8.0 Wenchuan earthquake and three dimensional P-wave velocity structure in and around source region. Chinese Journal of Geophysics (in Chinese), 52(2): 320-328.

    Wu R S, Aki K. 1985. Elastic wave scattering by a random medium and the small-scale inhomogeneities in the lithosphere. J. Geophys. Res., 90(B12): 10261-10273.

    Wu R S, Aki K. 1988. Introduction: seismic wave scattering in three-dimensionally heterogeneous earth. Pure Appl. Geophys., 128(1-2): 1-6. Yang H Y, Hu J F, Zhao H, et al. 2009. Crust-mantle structure and seismogenic background of Wenchuan MS8.0 earthquake in Western Sichuan area. Chinese Journal of Geophysics (in Chinese), 52(2): 356-364.

    Zheng X F, Ouyang B, Zhang D N, et al. 2009. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake. Chinese Journal of Geophysics (in Chinese), 52(5): 1412-1417, doi: 10.3969/j.issn.0001-5733.2009.05.031.

    Zhou L Q, Zhao C P, Xiu J G, et al. 2008. Tomography of QLgin Sichuan-Yunnan Zone. Chinese Journal of Geophysics (in Chinese), 51(6): 1745-1752.

    附中文參考文獻(xiàn)

    蔡學(xué)林, 曹家敏, 朱介壽等. 2008. 龍門山巖石圈地殼三維結(jié)構(gòu)及汶川大地震成因淺析. 城都理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 35(4): 357-365.

    鄧起東, 陳社發(fā), 趙小麟. 1994. 龍門山及其鄰區(qū)的構(gòu)造和地震活動(dòng)及動(dòng)力學(xué). 地震地質(zhì), 16(4): 389-403.

    郭飚, 劉啟元, 陳九輝等. 2009. 川西龍門山及鄰區(qū)地殼上地幔遠(yuǎn)震P波層析成像. 地球物理學(xué)報(bào), 52(2): 346-355.

    嘉世旭, 張先康. 2008. 青藏高原東北緣深地震測(cè)深震相研究與地殼細(xì)結(jié)構(gòu). 地球物理學(xué)報(bào), 51(5): 1431-1443.

    雷建設(shè), 趙大鵬, 蘇金蓉等. 2009. 龍門山斷裂帶地殼精細(xì)結(jié)構(gòu)與汶川地震發(fā)震機(jī)理. 地球物理學(xué)報(bào), 52(2): 339-345.

    李志偉, 胥頤, 黃潤(rùn)秋等. 2011. 龍門山地區(qū)的P波速度結(jié)構(gòu)與汶川地震的深部構(gòu)造特征. 中國(guó)科學(xué): 地球科學(xué), 41(3): 283-290.

    李永華, 吳慶舉, 田小波等. 2009. 用接收函數(shù)方法研究云南及其鄰區(qū)地殼上地幔結(jié)構(gòu). 地球物理學(xué)報(bào), 52(1): 67-80.

    劉建華, 劉福田, 吳華等. 1989. 中國(guó)南北帶地殼和上地幔的三維速度圖象. 地球物理學(xué)報(bào), 32(2): 143-151.

    劉啟元, 李昱, 陳九輝等. 2009. 汶川8.0地震: 地殼上地幔S波速度結(jié)構(gòu)的初步研究. 地球物理學(xué)報(bào), 52(2): 309-319.

    盧占武, 高 銳, 李秋生等. 2009. 橫過(guò)青藏高原羌塘地體中央隆起區(qū)的深反射地震試驗(yàn)剖面. 地球物理學(xué)報(bào), 52(8): 2008-2014.

    馬宏生, 汪素云, 裴順平等. 2007. 川滇及周邊地區(qū)地殼橫波衰減的成像研究. 地球物理學(xué)報(bào), 50(2): 465-471.

    石玉濤, 高原, 趙翠萍等. 2009. 汶川地震余震序列的地震各向異性. 地球物理學(xué)報(bào), 52(2): 398-407.

    蘇有錦, 劉杰, 鄭斯華等. 2006. 云南地區(qū)S波非彈性衰減Q值研究. 地震學(xué)報(bào), 28(2): 206-212.

    王椿鏞, Mooney W D, 王溪莉等. 2002. 川滇地區(qū)地殼上地幔三維速度結(jié)構(gòu)研究. 地震學(xué)報(bào), 24(1): 1-16.

    王椿鏞, 韓渭賓, 吳建平等. 2003a. 松潘—甘孜造山帶地殼速度結(jié)構(gòu). 地震學(xué)報(bào), 25(3): 229-311.

    王椿鏞, 吳建平, 樓海等. 2003b. 川西藏東地區(qū)的地殼P波速度結(jié)構(gòu). 中國(guó)科學(xué)(D輯), 33(S1): 181-189.

    王緒本, 朱迎堂, 趙錫奎等. 2009. 青藏高原東緣龍門山逆沖構(gòu)造深部電性結(jié)構(gòu)特征. 地球物理學(xué)報(bào), 52(2): 564-571.

    吳建平, 黃嬡, 張?zhí)熘械? 2009. 汶川MS8.0級(jí)地震余震分布及周邊區(qū)域P波三維速度結(jié)構(gòu)研究. 地球物理學(xué)報(bào), 52(2): 320-328.

    楊海燕, 胡家富, 趙宏等. 2009. 川西地區(qū)殼幔結(jié)構(gòu)與汶川MS8.0級(jí)地震的孕震背景. 地球物理學(xué)報(bào), 52(2): 356-364.

    趙珠, 張潤(rùn)生. 1987. 四川地區(qū)地殼上地幔速度結(jié)構(gòu)的初步研究. 地震學(xué)報(bào), 19(2): 154-166.

    鄭秀芬, 歐陽(yáng)飚, 張東寧等. 2009. “國(guó)家數(shù)字測(cè)震臺(tái)網(wǎng)數(shù)據(jù)備份中心”技術(shù)系統(tǒng)建設(shè)及其對(duì)汶川大地震研究的數(shù)據(jù)支撐. 地球物理學(xué)報(bào), 52(5): 1412-1417, doi: 10.3969/j.issn.0001-5733.2009.05.031.

    周連慶, 趙翠萍, 修濟(jì)剛等. 2008. 川滇地區(qū)Lg波Q值層析成像. 地球物理學(xué)報(bào), 51(6): 1745-1752.

    (本文編輯劉少華)

    基金項(xiàng)目國(guó)家自然科學(xué)基金(41074036和40974031)資助.

    作者簡(jiǎn)介范小平,男,1974年生,博士,正研級(jí)高級(jí)工程師,主要從事地震波散射和衰減研究. E-mail:nj_fxp@163.com

    doi:10.6038/cjg20160714 中圖分類號(hào)P315

    收稿日期2015-08-08,2016-02-02收修定稿

    Variations of the crust seismic scattering strength below the southeastern margin of the Tibetan plateau and adjacent regions

    FAN Xiao-Ping, WANG Jun-Fei, YANG Cong-Jie, LI Qing-He, HE Yi-Cheng, BI Xue-Mei

    EarthquakeAdministrationofJiangsuProvince,Nanjing210014,China

    AbstractThe ten broadband teleseismic waveform data recorded by forty-seven digital seismic networks were selected to study the seismic scattering strength of crust in the southeastern margin of the Tibetan plateau and adjacent regions by the teleseismic fluctuation wavefield method (TFWM). The results show that there exists strong transverse inhomogeneity beneath the southeastern Tibetan plateau and adjacent regions. The scattering strength of the Yangtsz block and the Tibetan block has obvious difference. The gradient zone of seismic scattering strength is along the main fault zone. The position of the Longmenshan fault and the Anninghe fault have strong correlation with the gradient zone of seismic scattering strength. In addition, the high value region of strong seismic scattering strength is almost turning to the Tibet side. The epicenters of strong earthquakes are along the gradient zone of seismic scattering strength, and lie the side of strong scattering strength. The scattering strength shows the characteristic of vertical inhomogeneity and the scattering of lower crust is stronger than the upper crust. The spacial distribution of the strong and the weak scattering strengths maybe related to the strong movement of geological tectonics, the broken medium caused by frequent tectonic and volcanic activity, and also with the heat flow upwelling along active faults.

    KeywordsSoutheastern of margin of Tibetan plateau; Scattering strength; Seismic scattering; Teleseismic

    范小平, 王俊菲, 楊從杰等. 2016. 青藏高原東南緣及鄰區(qū)地殼介質(zhì)散射強(qiáng)度變化特征. 地球物理學(xué)報(bào),59(7):2486-2497,doi:10.6038/cjg20160714.

    Fan X P, Wang J F, Yang C J, et al. 2016. Variations of the crust seismic scattering strength below the southeastern margin of the Tibetan plateau and adjacent regions. Chinese J. Geophys. (in Chinese),59(7):2486-2497,doi:10.6038/cjg20160714.

    猜你喜歡
    波場(chǎng)塊體斷裂帶
    冷凍斷裂帶儲(chǔ)層預(yù)測(cè)研究
    依蘭—伊通斷裂帶黑龍江段構(gòu)造運(yùn)動(dòng)特征
    地震研究(2021年1期)2021-04-13 01:05:12
    一種新型單層人工塊體Crablock 的工程應(yīng)用
    彈性波波場(chǎng)分離方法對(duì)比及其在逆時(shí)偏移成像中的應(yīng)用
    交錯(cuò)網(wǎng)格與旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離的影響分析
    基于Hilbert變換的全波場(chǎng)分離逆時(shí)偏移成像
    準(zhǔn)噶爾盆地西北緣克-夏斷裂帶構(gòu)造特征新認(rèn)識(shí)
    一種Zr 基塊體金屬玻璃的納米壓入蠕變行為研究
    上海金屬(2015年3期)2015-11-29 01:09:58
    郯廬斷裂帶及兩側(cè)地區(qū)強(qiáng)震異常特征分析
    塊體非晶合金及其應(yīng)用
    狂野欧美激情性bbbbbb| 久久久久久久大尺度免费视频| 午夜福利网站1000一区二区三区| 久久99蜜桃精品久久| 国产男女内射视频| 久久午夜福利片| 国产精品久久久久久精品古装| tube8黄色片| 欧美激情极品国产一区二区三区 | 有码 亚洲区| 春色校园在线视频观看| 老熟女久久久| 国产在线免费精品| 哪个播放器可以免费观看大片| 久久久久人妻精品一区果冻| 麻豆精品久久久久久蜜桃| 精品一品国产午夜福利视频| a级毛片在线看网站| 亚洲国产欧美在线一区| 亚洲国产色片| 成年女人在线观看亚洲视频| 高清黄色对白视频在线免费看| 搡女人真爽免费视频火全软件| 国产成人av激情在线播放 | 午夜福利在线观看免费完整高清在| 一区二区av电影网| av国产久精品久网站免费入址| 亚洲av不卡在线观看| 欧美激情 高清一区二区三区| 精品熟女少妇av免费看| 亚洲内射少妇av| 大香蕉97超碰在线| 国产精品麻豆人妻色哟哟久久| 免费高清在线观看日韩| 日韩大片免费观看网站| 亚洲四区av| 99久久人妻综合| 亚洲经典国产精华液单| 久久久久久久久大av| 亚洲成人手机| 老女人水多毛片| 国产熟女欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 性色avwww在线观看| 香蕉精品网在线| 亚洲国产精品999| 视频在线观看一区二区三区| 欧美日韩视频高清一区二区三区二| 国产精品国产三级国产专区5o| 精品人妻在线不人妻| 午夜福利网站1000一区二区三区| 国产片内射在线| 午夜激情福利司机影院| 人妻一区二区av| freevideosex欧美| 色视频在线一区二区三区| 日本黄色片子视频| 少妇高潮的动态图| 亚洲一级一片aⅴ在线观看| 久久午夜福利片| 亚洲国产精品国产精品| 欧美xxⅹ黑人| 黑人猛操日本美女一级片| 日韩 亚洲 欧美在线| 一级毛片aaaaaa免费看小| 我的老师免费观看完整版| 精品久久蜜臀av无| av女优亚洲男人天堂| 男女无遮挡免费网站观看| 男女免费视频国产| 国产男女内射视频| 九色亚洲精品在线播放| 久久久久久久久久成人| 最新中文字幕久久久久| 国产亚洲精品久久久com| 在线观看人妻少妇| 满18在线观看网站| 91久久精品国产一区二区三区| 九色亚洲精品在线播放| 精品人妻熟女毛片av久久网站| 制服诱惑二区| 卡戴珊不雅视频在线播放| 久久久久久久国产电影| 国产有黄有色有爽视频| 97在线视频观看| 少妇熟女欧美另类| 国产成人精品福利久久| 黑人高潮一二区| 美女国产高潮福利片在线看| 色5月婷婷丁香| 日本与韩国留学比较| 成人漫画全彩无遮挡| 久久久久国产精品人妻一区二区| 熟女av电影| 国产黄频视频在线观看| 毛片一级片免费看久久久久| 亚洲av成人精品一区久久| 在线精品无人区一区二区三| 又大又黄又爽视频免费| av不卡在线播放| 青春草视频在线免费观看| 纵有疾风起免费观看全集完整版| 在线观看美女被高潮喷水网站| 女的被弄到高潮叫床怎么办| 国产av码专区亚洲av| 97超视频在线观看视频| 女性被躁到高潮视频| h视频一区二区三区| 天天躁夜夜躁狠狠久久av| 国产亚洲午夜精品一区二区久久| 美女视频免费永久观看网站| 日本-黄色视频高清免费观看| 成人国产av品久久久| 精品国产一区二区三区久久久樱花| 国产国语露脸激情在线看| 男男h啪啪无遮挡| 十八禁高潮呻吟视频| 在线天堂最新版资源| 又粗又硬又长又爽又黄的视频| 水蜜桃什么品种好| 国产成人freesex在线| 日韩精品有码人妻一区| 成人无遮挡网站| 亚洲欧美色中文字幕在线| 久久热精品热| 亚洲,一卡二卡三卡| 黑人巨大精品欧美一区二区蜜桃 | 各种免费的搞黄视频| 亚洲精品日本国产第一区| 日韩欧美一区视频在线观看| 日产精品乱码卡一卡2卡三| 久久久久久久精品精品| 国产精品免费大片| 国产乱来视频区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久久久久亚洲中文字幕| 一区二区三区精品91| 亚洲精品久久午夜乱码| 精品人妻一区二区三区麻豆| 边亲边吃奶的免费视频| 久久人人爽av亚洲精品天堂| 91成人精品电影| 我的老师免费观看完整版| 中文字幕最新亚洲高清| 乱人伦中国视频| 国产免费福利视频在线观看| 人人澡人人妻人| 狠狠精品人妻久久久久久综合| 在线天堂最新版资源| 国产亚洲午夜精品一区二区久久| 国产爽快片一区二区三区| 欧美三级亚洲精品| 夫妻午夜视频| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av涩爱| 日韩制服骚丝袜av| 亚洲国产精品成人久久小说| 人人妻人人爽人人添夜夜欢视频| 99热网站在线观看| 国产精品国产三级国产专区5o| 免费大片黄手机在线观看| av免费在线看不卡| 在线观看www视频免费| 成人免费观看视频高清| 午夜久久久在线观看| 天堂俺去俺来也www色官网| av有码第一页| 国产亚洲最大av| 飞空精品影院首页| 伊人亚洲综合成人网| 国产成人av激情在线播放 | 成人黄色视频免费在线看| 99九九在线精品视频| 国产淫语在线视频| 99热国产这里只有精品6| 亚洲av免费高清在线观看| 18+在线观看网站| 精品久久久久久久久亚洲| 午夜福利影视在线免费观看| 亚洲久久久国产精品| 亚洲国产av新网站| 男人爽女人下面视频在线观看| 尾随美女入室| 免费不卡的大黄色大毛片视频在线观看| 国产精品99久久99久久久不卡 | 国内精品宾馆在线| 人妻制服诱惑在线中文字幕| 免费av中文字幕在线| 一边摸一边做爽爽视频免费| 在线天堂最新版资源| 国产视频首页在线观看| 国产日韩欧美亚洲二区| av播播在线观看一区| 精品少妇内射三级| 赤兔流量卡办理| 日韩精品有码人妻一区| 日韩一区二区三区影片| 制服丝袜香蕉在线| 日日爽夜夜爽网站| 男人添女人高潮全过程视频| 亚洲欧美成人综合另类久久久| 国产精品久久久久久精品电影小说| 99国产精品免费福利视频| 99热网站在线观看| 晚上一个人看的免费电影| 久久99一区二区三区| 最近中文字幕高清免费大全6| 国产高清三级在线| 在线观看免费视频网站a站| 免费黄网站久久成人精品| 大陆偷拍与自拍| 亚洲精品乱码久久久v下载方式| 日本免费在线观看一区| 高清不卡的av网站| 久久久欧美国产精品| 国产伦精品一区二区三区视频9| 简卡轻食公司| 亚洲国产日韩一区二区| 91国产中文字幕| 久久久国产精品麻豆| 精品一区二区三卡| 国产精品一区二区在线观看99| 成人毛片a级毛片在线播放| 成人毛片a级毛片在线播放| 亚洲av福利一区| 久久久久网色| 日韩电影二区| 亚洲精品aⅴ在线观看| 亚洲国产欧美在线一区| 亚洲国产色片| 久久毛片免费看一区二区三区| 久久99精品国语久久久| 成人18禁高潮啪啪吃奶动态图 | 在线观看人妻少妇| 永久网站在线| 九色成人免费人妻av| 欧美xxxx性猛交bbbb| 国产精品国产av在线观看| 超色免费av| 男女边吃奶边做爰视频| 午夜福利视频在线观看免费| 国产极品粉嫩免费观看在线 | 美女脱内裤让男人舔精品视频| 国产高清国产精品国产三级| 日本欧美视频一区| 亚洲伊人久久精品综合| 久久久精品区二区三区| 大话2 男鬼变身卡| 一本久久精品| 亚洲av二区三区四区| 久久影院123| 久久久久久久亚洲中文字幕| 少妇丰满av| 亚洲欧美日韩另类电影网站| 欧美精品一区二区大全| 久久国产亚洲av麻豆专区| 欧美成人午夜免费资源| 丝袜在线中文字幕| 国产精品熟女久久久久浪| 热99国产精品久久久久久7| 亚洲国产成人一精品久久久| 欧美日韩综合久久久久久| 免费久久久久久久精品成人欧美视频 | 欧美激情极品国产一区二区三区 | 中文字幕最新亚洲高清| 亚洲成人一二三区av| 美女内射精品一级片tv| 美女cb高潮喷水在线观看| 久久精品熟女亚洲av麻豆精品| 视频在线观看一区二区三区| 色视频在线一区二区三区| av线在线观看网站| 高清午夜精品一区二区三区| 午夜久久久在线观看| 天堂俺去俺来也www色官网| 黄色配什么色好看| 中国三级夫妇交换| 亚洲欧美一区二区三区国产| 亚洲精品一二三| 一区在线观看完整版| 中文字幕精品免费在线观看视频 | 日本欧美视频一区| 99国产精品免费福利视频| 最黄视频免费看| 日韩强制内射视频| 欧美成人精品欧美一级黄| 黑丝袜美女国产一区| 欧美精品人与动牲交sv欧美| 免费久久久久久久精品成人欧美视频 | 18+在线观看网站| 日产精品乱码卡一卡2卡三| av国产精品久久久久影院| 99热这里只有精品一区| 欧美xxⅹ黑人| 中国国产av一级| 九九爱精品视频在线观看| 我的老师免费观看完整版| 欧美日韩视频高清一区二区三区二| av线在线观看网站| 成人毛片60女人毛片免费| 97在线视频观看| 18禁动态无遮挡网站| 精品国产一区二区三区久久久樱花| 精品久久久久久久久av| 国产无遮挡羞羞视频在线观看| 久久久久久久久久人人人人人人| 亚洲综合色网址| 国产视频内射| 99国产综合亚洲精品| 欧美少妇被猛烈插入视频| 国产精品蜜桃在线观看| 插逼视频在线观看| 最近最新中文字幕免费大全7| 免费黄网站久久成人精品| 制服诱惑二区| 999精品在线视频| 午夜日本视频在线| 成年人午夜在线观看视频| 国产av码专区亚洲av| 亚洲高清免费不卡视频| 欧美最新免费一区二区三区| 在线观看免费高清a一片| 欧美+日韩+精品| 男女高潮啪啪啪动态图| 国产又色又爽无遮挡免| 亚洲欧洲精品一区二区精品久久久 | 人成视频在线观看免费观看| 五月天丁香电影| 免费观看av网站的网址| 免费日韩欧美在线观看| 国产成人精品久久久久久| 好男人视频免费观看在线| 九草在线视频观看| 少妇猛男粗大的猛烈进出视频| 欧美日本中文国产一区发布| 91在线精品国自产拍蜜月| 欧美激情国产日韩精品一区| 久久99热这里只频精品6学生| 久久久精品免费免费高清| 人人妻人人澡人人看| 久久99热这里只频精品6学生| 国产毛片在线视频| 午夜福利在线观看免费完整高清在| 大码成人一级视频| 麻豆乱淫一区二区| 成人亚洲精品一区在线观看| 国产亚洲午夜精品一区二区久久| 久久久久久久国产电影| 亚洲欧美精品自产自拍| 精品少妇久久久久久888优播| 欧美精品高潮呻吟av久久| 亚洲美女黄色视频免费看| 热99国产精品久久久久久7| 欧美人与善性xxx| 一本—道久久a久久精品蜜桃钙片| 欧美丝袜亚洲另类| 99视频精品全部免费 在线| 国产精品女同一区二区软件| 又黄又爽又刺激的免费视频.| 91国产中文字幕| 老司机亚洲免费影院| 国产成人freesex在线| 中文精品一卡2卡3卡4更新| 国产成人精品一,二区| 亚洲精品色激情综合| 国产一区有黄有色的免费视频| 一本—道久久a久久精品蜜桃钙片| 91成人精品电影| 一级毛片 在线播放| 涩涩av久久男人的天堂| 在线观看三级黄色| 成人午夜精彩视频在线观看| 亚洲国产精品一区三区| 久久精品国产亚洲av涩爱| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 人人妻人人添人人爽欧美一区卜| 国产精品一区二区三区四区免费观看| 美女福利国产在线| 99视频精品全部免费 在线| 性色avwww在线观看| av.在线天堂| 国产精品嫩草影院av在线观看| 天堂中文最新版在线下载| 一本一本综合久久| 午夜福利网站1000一区二区三区| 亚洲精品第二区| av在线观看视频网站免费| 国产精品国产三级国产专区5o| 日本与韩国留学比较| 亚洲av综合色区一区| 观看av在线不卡| 免费黄网站久久成人精品| 国产一区二区在线观看av| 涩涩av久久男人的天堂| 大片电影免费在线观看免费| kizo精华| 国精品久久久久久国模美| 欧美人与善性xxx| 亚洲色图综合在线观看| 色94色欧美一区二区| 日韩欧美精品免费久久| 国产成人精品无人区| 一本一本综合久久| 亚洲精品视频女| 国产亚洲一区二区精品| 久久久久精品久久久久真实原创| 乱码一卡2卡4卡精品| 精品久久蜜臀av无| 亚洲人成网站在线观看播放| 日韩成人av中文字幕在线观看| 成年人免费黄色播放视频| 午夜福利视频精品| 狠狠婷婷综合久久久久久88av| 免费日韩欧美在线观看| 午夜日本视频在线| 晚上一个人看的免费电影| 亚洲欧美成人综合另类久久久| 精品人妻一区二区三区麻豆| 2018国产大陆天天弄谢| 亚洲精华国产精华液的使用体验| 免费观看性生交大片5| 国产成人精品福利久久| 夜夜骑夜夜射夜夜干| 两个人的视频大全免费| 99热6这里只有精品| 欧美日韩在线观看h| 中文字幕制服av| www.色视频.com| 国产色爽女视频免费观看| 少妇人妻 视频| 满18在线观看网站| av一本久久久久| 日韩视频在线欧美| 又大又黄又爽视频免费| 亚洲四区av| 国产白丝娇喘喷水9色精品| 亚洲精品色激情综合| 999精品在线视频| 满18在线观看网站| 精品久久久精品久久久| 一级毛片黄色毛片免费观看视频| 水蜜桃什么品种好| 成人黄色视频免费在线看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲情色 制服丝袜| 色5月婷婷丁香| 91久久精品国产一区二区成人| 日本黄大片高清| 丝袜脚勾引网站| 亚洲精品自拍成人| 飞空精品影院首页| 日韩大片免费观看网站| 少妇被粗大的猛进出69影院 | 最近手机中文字幕大全| 亚洲精品aⅴ在线观看| 国产乱来视频区| 26uuu在线亚洲综合色| 高清在线视频一区二区三区| 大码成人一级视频| 亚洲国产精品999| 久久韩国三级中文字幕| 亚洲欧美一区二区三区黑人 | 五月天丁香电影| 亚洲精品久久久久久婷婷小说| 亚洲精品中文字幕在线视频| 国产成人免费观看mmmm| av专区在线播放| 9色porny在线观看| 亚洲av二区三区四区| 五月天丁香电影| 只有这里有精品99| 久久精品国产鲁丝片午夜精品| 国产亚洲欧美精品永久| 久久av网站| 高清不卡的av网站| 男女免费视频国产| 国产无遮挡羞羞视频在线观看| 九九久久精品国产亚洲av麻豆| 亚洲三级黄色毛片| 日本av手机在线免费观看| 九九在线视频观看精品| 亚洲欧美一区二区三区国产| 少妇丰满av| 日本欧美视频一区| 综合色丁香网| 国产av码专区亚洲av| 51国产日韩欧美| 一区在线观看完整版| 乱码一卡2卡4卡精品| a级毛片免费高清观看在线播放| 黄色配什么色好看| 老司机影院成人| 中文欧美无线码| 97在线人人人人妻| 伦理电影免费视频| 国产白丝娇喘喷水9色精品| 在线 av 中文字幕| 男人操女人黄网站| 一级爰片在线观看| 色视频在线一区二区三区| 插逼视频在线观看| a级毛色黄片| 亚洲怡红院男人天堂| 国产精品麻豆人妻色哟哟久久| 97超视频在线观看视频| 一本一本综合久久| 极品少妇高潮喷水抽搐| 精品国产乱码久久久久久小说| 欧美日韩av久久| 天天躁夜夜躁狠狠久久av| 亚洲精品第二区| 91aial.com中文字幕在线观看| 亚洲精品456在线播放app| 欧美日韩精品成人综合77777| 亚洲精品乱码久久久v下载方式| 亚洲国产av新网站| 新久久久久国产一级毛片| 一本色道久久久久久精品综合| 精品一区二区三卡| 一级毛片电影观看| 亚洲,欧美,日韩| 国产高清国产精品国产三级| 国产精品人妻久久久久久| 一级二级三级毛片免费看| 91精品三级在线观看| 午夜免费鲁丝| 国产在线免费精品| 亚洲婷婷狠狠爱综合网| 国产有黄有色有爽视频| 久久久久精品性色| 亚洲av国产av综合av卡| 韩国高清视频一区二区三区| 国产熟女午夜一区二区三区 | 我的老师免费观看完整版| 97精品久久久久久久久久精品| 欧美精品一区二区免费开放| 国产欧美日韩一区二区三区在线 | 亚洲av成人精品一区久久| 欧美国产精品一级二级三级| 免费看不卡的av| 最近中文字幕高清免费大全6| 91精品三级在线观看| 国产精品一区二区在线不卡| 久久久久久久久久成人| av国产精品久久久久影院| 美女脱内裤让男人舔精品视频| 久久97久久精品| 久久久久国产网址| 91精品国产国语对白视频| 久久久久久久久久人人人人人人| 在现免费观看毛片| 日韩av不卡免费在线播放| 国产乱来视频区| 亚洲精品视频女| 妹子高潮喷水视频| 久久精品国产自在天天线| 在线亚洲精品国产二区图片欧美 | 少妇熟女欧美另类| 欧美人与善性xxx| 国产高清国产精品国产三级| 亚洲性久久影院| 啦啦啦视频在线资源免费观看| 黄色一级大片看看| 久久97久久精品| 久久久久久久久久久免费av| 18禁在线无遮挡免费观看视频| 人成视频在线观看免费观看| 亚洲精品乱码久久久v下载方式| 精品久久国产蜜桃| 亚洲成色77777| 国产精品国产三级国产专区5o| 各种免费的搞黄视频| 亚洲久久久国产精品| 精品一区二区三区视频在线| 少妇人妻 视频| 在线播放无遮挡| 欧美 亚洲 国产 日韩一| 日韩中字成人| 在线观看免费视频网站a站| 黄色毛片三级朝国网站| 亚洲av欧美aⅴ国产| 亚洲精品一区蜜桃| 在线观看美女被高潮喷水网站| 午夜福利网站1000一区二区三区| 久久精品国产亚洲av天美| 久久午夜综合久久蜜桃| 欧美xxⅹ黑人| 久久久久久人妻| 国产高清国产精品国产三级| 永久网站在线| 九九爱精品视频在线观看| 亚洲国产精品一区三区| 高清黄色对白视频在线免费看| 亚洲不卡免费看| 久久综合国产亚洲精品| 天天操日日干夜夜撸| 亚洲精品乱码久久久v下载方式| 亚洲av电影在线观看一区二区三区| 亚洲欧美一区二区三区黑人 | 超碰97精品在线观看| 午夜久久久在线观看| 九九在线视频观看精品| 欧美日韩亚洲高清精品| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产av新网站| 人人妻人人爽人人添夜夜欢视频| 亚洲国产精品成人久久小说| 亚洲国产精品专区欧美| 免费看光身美女| 最近2019中文字幕mv第一页| 97超视频在线观看视频| 女性生殖器流出的白浆| 国产69精品久久久久777片| 高清av免费在线|