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

    川西高原重磁異常特征與構(gòu)造背景分析

    2015-03-01 01:37:06高玲舉張健董淼
    地球物理學(xué)報 2015年8期
    關(guān)鍵詞:鮮水河川西高原理塘

    高玲舉, 張健*, 董淼

    1 中國科學(xué)院計算地球動力學(xué)重點實驗室, 北京 100049 2 中國科學(xué)院大學(xué), 北京 100049

    ?

    川西高原重磁異常特征與構(gòu)造背景分析

    高玲舉1,2, 張健1,2*, 董淼1,2

    1 中國科學(xué)院計算地球動力學(xué)重點實驗室, 北京 100049 2 中國科學(xué)院大學(xué), 北京 100049

    川西高原位于青藏高原東緣,是我國大陸地殼構(gòu)造變形及地震活動最強烈的區(qū)域.利用最新重力、航磁資料,通過異常分析和反演計算,研究了該區(qū)鮮水河斷裂、理塘斷裂、金沙江斷裂的重磁異常特征、莫霍面特征、居里面特征,分析得出了這些斷裂的深部地質(zhì)結(jié)構(gòu)與構(gòu)造背景.計算表明:川西高原莫霍面東南淺、西北深,地殼厚度在43~63 km之間.居里面特征表現(xiàn)為條帶狀,深度在17~23 km之間.其中,鮮水河斷裂帶對應(yīng)莫霍面深度梯度帶,居里面為高低起伏圈閉.理塘斷裂帶北段莫霍面局部隆坳相間,南段莫霍面逐漸抬升,居里面呈現(xiàn)由西向東加深的梯度帶.金沙江斷裂帶,居里面形成局部抬升,深部可能存在高溫地?zé)岙惓T?綜合分析認為,川西高原地殼結(jié)構(gòu)主要特點為:增厚的下地殼,熱-塑性變形的中地殼,脆性變形的上地殼.

    川西高原; 川滇菱形塊體; 地殼結(jié)構(gòu); 重磁反演

    1 引言

    川西高原地處川滇活動塊體北部,是我國新構(gòu)造運動最活躍、現(xiàn)代地震最密集的地區(qū).川西高原的主體是川西北活動塊體,是川滇菱形塊體的重要組成部分(徐錫偉等,2003;張燕等,2013;喬學(xué)軍等,2004;吳建平等,2006;徐錫偉等,2005).川西北活動塊體內(nèi)部,發(fā)育有一系列近南北向弧形活動斷裂帶.這些活動斷裂帶規(guī)模不等、力學(xué)性質(zhì)不同,對歷史、現(xiàn)今和未來破壞性地震具有明顯的空間限制作用,比如,2014年11月22日的康定地震,就發(fā)生在鮮水河斷裂帶上.

    川西高原的地震活動與斷裂帶密切相關(guān),而斷裂帶的性質(zhì)、規(guī)模、活動方式受深部地質(zhì)結(jié)構(gòu)與構(gòu)造條件控制.前人依據(jù)川滇塊體及鄰近地區(qū)大量觀測數(shù)據(jù)和理論研究(張培震,2008;王椿墉等,2003;王閻昭等,2008;王慶良等,2008;黃金莉等,2001;李飛等,2011;崔仲雄和裴順平,2009;韓渭賓和蔣國芳,2004;程佳等,2009;聞學(xué)澤,2000;羅鈞等,2014;徐朝繁等,2008;丁志峰等,1999;李本亮等,2011;樓海和王椿鏞,2005;周榮軍等,2005;劉援朝等,2010;劉蓓莉,1993;吳建平等,2006;唐文清等,2005),對該區(qū)構(gòu)造樣式、活動構(gòu)造基本特征、主要斷裂現(xiàn)今活動速率、構(gòu)造變形與應(yīng)變、地殼垂直運動、強震活動分布特征與構(gòu)造背景、震群發(fā)震構(gòu)造及動力來源、震源機制和應(yīng)力場特征、地震破裂特征等開展了大量研究,特別在川滇活動構(gòu)造區(qū)地震層析成像、地殼上地幔三維速度結(jié)構(gòu)、Pn波速度結(jié)構(gòu)和各向異性、大地電磁觀測地殼物質(zhì)流、重力異常小波分解研究地殼深部構(gòu)造等方面開展了大量卓有成效研究.

    上述研究,依據(jù)地震活動特點和規(guī)律,在塊體邊界、塊體相互關(guān)系、運動學(xué)特征方面等獲得了豐富的認識.但由于塊體內(nèi)部變形較小,地震活性相對較弱,資料缺乏,塊體內(nèi)部深部結(jié)構(gòu)的研究往往被忽視.本文利用最新重磁資料(本文的重力資料來自于加州大學(xué)圣迭戈分校斯克里普斯海洋研究所(http:∥topex.ucsd.edu/),數(shù)據(jù)精度為1′×1′,重力模型基于V23.1;航磁資料來源于中國國土資源部航空物探遙感中心,數(shù)據(jù)的網(wǎng)度為2 km×2 km),開展聯(lián)合反演分析,不但研究塊體邊界斷層,而且針對活動塊體內(nèi)部地殼結(jié)構(gòu),開展計算分析.在重磁反演基礎(chǔ)上,分析川滇菱形塊體內(nèi)部的川西北次級塊體的綜合地質(zhì)、地球物理特征,重點探討地殼結(jié)構(gòu)背景與活動斷裂帶構(gòu)造關(guān)系,以及塊體內(nèi)部地質(zhì)結(jié)構(gòu).

    2 新生代構(gòu)造活動特點及地球物理特征

    川西高原位于青藏高原東緣,具有獨特的構(gòu)造變形特征.青藏高原東緣自第四紀以來的大規(guī)模水平剪切變形和強烈隆升,導(dǎo)致川西高原地區(qū)被一系列南北向、北西向活動斷裂分割(圖1a).

    川西高原地處龍門山、茶坪山、夾金山、貢嘎山、錦屏山以西,平均海撥3000~4000 m(圖1b),屬青藏高原東延部分.研究區(qū)內(nèi),各種體系的構(gòu)造和斷裂,互相重合、交接、穿插,并有多期巖漿巖侵入,是我國著名的地震活動區(qū).川西高原內(nèi)的川西北塊體,四周被北西向甘孜—玉樹斷裂帶、鮮水河斷裂帶、北東向麗江—小金河斷裂帶、金沙江斷裂帶和德欽—中甸—大具斷裂等圍限.甘孜—玉樹斷裂帶自青海玉樹向東南至四川甘孜,長約270 km,沿線斷錯地貌顯示強烈左旋走滑運動.鮮水河斷裂帶也是左旋走滑型斷裂,走向N40°W,傾角近直立,略呈向北東凸出的弧形彎曲.起自四川甘孜,與甘孜—玉樹斷裂帶左階羽列,在階列區(qū)形成寬22~35 km的甘孜拉分區(qū),向南經(jīng)爐霍、道孚,錯斷一系列山脊水系和地質(zhì)體,過康定后沿大渡河谷延伸,全長約400 km,300年來沿該斷裂發(fā)生過4次7級以上強震.近期地震活動十分活躍,如2014年11月22日康定地震(圖1a中紅色五角星標注).麗江—小金河斷裂帶是在中新生代龍門山—錦屏山—玉龍雪山推覆構(gòu)造帶南西段基礎(chǔ)上形成的一條北東向活動斷裂帶,自南西劍川盆地,向北東經(jīng)麗江、木里等地,止于安寧河的冕寧以北,由多條斜列次級斷層組成.北東段左旋拉張、中段左旋走滑、南西段擠壓逆沖.整體走向N45—50°E,傾向北西,傾角60°~80°.以左旋走滑為主,兼有垂(逆)向抬升,北西(上)盤相對東南(下)盤抬升,是逆傾左滑型斷裂.金沙江斷裂帶和德欽—中甸—大具斷裂現(xiàn)今活動性很弱,是一組斷續(xù)、分散、滑動速率較低的右旋走滑斷裂(徐錫偉等,2003;張燕等,2013;喬學(xué)軍等,2004;吳建平等,2006;徐錫偉等,2005).

    川西北塊體的內(nèi)部還發(fā)育理塘斷裂、巴塘斷裂、玉農(nóng)希斷裂等次級活動斷裂.其中,理塘斷裂是一條與鮮水河斷裂近于平行的左旋走滑活動斷裂,北西起自金沙江畔,向東南延伸,經(jīng)毛埡壩、理塘、甲洼、德巫,至木里以北消失,全長約385 km.斷裂走向N40—50°W,總體傾向NE,傾角較陡,控制毛埡壩、理塘、甲洼、德巫等盆地的發(fā)展,并使古近紀—第四紀地層普遍遭受褶皺或斷錯作用.理塘盆地以東主要是由3條50~65 km長、具明顯活動性的次級斷層(即理塘斷裂、擦忠斷裂和木拉斷裂)斜列而成.巴塘斷裂是一條右旋走滑活動斷裂,北東起自莫西附近,向南西延伸,經(jīng)巴塘、過金沙江、穿莽嶺、至瀾滄江邊消失,全長約200 km.斷裂走向約N30°E,總體傾向NW,傾角較陡,右旋走滑斜切金沙江構(gòu)造帶,系陸內(nèi)變形作用的產(chǎn)物,生成時間稍晚于金沙江構(gòu)造帶.玉農(nóng)希斷裂是一條與麗江—小金河斷裂帶北東段近于平行的逆傾右滑型斷裂,北東起自鮮水河斷裂上的康定,向東南延伸,至理塘斷裂附近消失.玉農(nóng)希斷裂是貢嘎山第四紀強烈斷塊隆起的西邊界,具有較強的活動性,貢嘎山及其周邊地區(qū)存在整體抬升的趨勢,來自北西側(cè)地殼物質(zhì)向南東方向的運動在貢嘎山地區(qū)受到阻擋,造成玉農(nóng)希斷裂的擠壓及其附近地殼的縮短,引起貢嘎山的隆升(周榮軍等,2005;劉援朝等,2010;劉蓓莉,1993;吳建平等,2006;唐文清等,2005).

    上述活動斷裂按運動學(xué)性質(zhì)各不相同,例如,塊體東界的鮮水河斷裂帶為正走滑型斷裂,塊體內(nèi)部的理塘斷裂帶為正傾滑型斷裂,塊體西界的金沙江斷裂帶為逆傾滑型斷裂.不僅運動學(xué)性質(zhì)不同,地球物理場特征也有明顯區(qū)別.

    重磁資料表明,布格異?;咎卣魇俏鞅钡?、東南高,多為串珠狀的異常(圖1c).總體上,川西次級塊體造山帶密度相對較低.川西北次級塊體出露的地層幾乎全是三疊系,沉積厚度巨大.大部分地區(qū)缺失侏羅系和白堊系地層,第三系直接覆蓋于三疊系之上,構(gòu)成密度差值達0.31~0.8 g·cm-3的明顯界面.三疊系上統(tǒng)內(nèi)的界面密度差值在0.11 g·cm-3左右.二疊系上統(tǒng)與下統(tǒng)之間界面密度差值為0.13 g·cm-3(劉蓓莉,1993).鮮水河斷裂帶處于密度平穩(wěn)地區(qū),斷裂帶兩側(cè)密度差異不明顯,沿斷裂帶也沒有強度較高的小局部異常存在,表明沿鮮水河斷裂帶沒有發(fā)生過大規(guī)模的深部巖漿上涌活動.沿金沙江斷裂分布著一系列局部高低密度異常,表明沿斷裂有不同密度的物質(zhì)分布,有深部高密度物質(zhì)侵入中上地殼中.從航磁異常圖上可以看出(圖1d),本區(qū)的磁場特征從西向東具有明顯分帶的特征,大致以甘孜—理塘和魯通—灣壩一線為界可劃分白玉—德欽劇烈變化磁場區(qū)、雅江—九龍平緩磁場區(qū)和瀘定—冕寧強磁場區(qū).整體上看,該區(qū)具一定規(guī)模的北西向構(gòu)造,在一定深度上被超過磁性基底的東西向斷裂構(gòu)造切割.

    圖1 川西高原區(qū)域地質(zhì)—地球物理圖(a) 區(qū)域地質(zhì)背景圖; (b) 地形圖; (c) 布格重力異常圖; (d) 航磁ΔT化極圖.Fig.1 Western Sichuan Plateau regional geological and geophysical maps(a) Regional geological background map; (b) Topographic map; (c) Bouguer gravity anomaly map; (d) Aeromagnetic polarization map.

    二維地殼速度結(jié)構(gòu)表明(王椿墉等,2003),川西地殼15~30 km深度之間,存在厚度8~10 km的低速層,P波速度為5.8 km·s-1.該低速層對應(yīng)電阻率在幾十到上百歐姆米的高導(dǎo)層.殼內(nèi)低速層通常認為與巖石的孔隙壓力和含水性質(zhì)有關(guān),也可能是由于介質(zhì)熱狀態(tài)變化或特殊礦物成分所引起的.殼內(nèi)高導(dǎo)層則被認為是由部分熔融物質(zhì)所組成.因此,川西上地殼內(nèi)的低速、高導(dǎo)層可能與介質(zhì)熱狀態(tài)有關(guān).川西地震震源深度集中在上地殼上部0~15 km范圍.中地殼15~30 km深度范圍內(nèi),普遍存在厚約5 km的缺震層.下地殼30~55 km深度范圍,地震活動重新出現(xiàn).此外,理塘斷裂帶表現(xiàn)為向西傾斜的鏟形高、低速分界帶,向深部歸并到上地殼低速—高導(dǎo)層內(nèi),即低速—高導(dǎo)層為上述2條逆斷裂帶的推覆滑脫體.

    上述地球物理場特征表明,川西北塊體5~15 km的上地殼是地殼脆性孕震層位,5~30 km的中地殼內(nèi)存在低速高導(dǎo)層,可能存在溫度相對較高的局部熔融體或糜棱巖化的韌性變形體(層),下地殼增厚并與上地幔形成明顯的熱聯(lián)系.這種深部熱過程驅(qū)動、中部塑性流變、淺部脆性地殼內(nèi)應(yīng)力積累的構(gòu)造格架,是川西斷裂構(gòu)造活動、地震多發(fā)的原因.

    3 計算方法與結(jié)果討論

    3.1 計算方法

    3.1.1 計算居里面的功率譜法

    功率譜法是目前計算居里面深度的常用方法,由于功率譜法估算磁性體深度不需要物性資料就能

    完成,減少了對先驗資料的要求,所以該方法具有方便、快捷的優(yōu)點,在國內(nèi)外得到了廣泛的發(fā)展和應(yīng)用(Bhattacharyya,1966;Spector and Grant,1970;Bhattacharyya and Leu,1975;Blakely,1988;Tanaka et al.,1999;郝書儉等,1982;申寧華,1985;侯重初和李保國,1985b;趙百民和郝天珧,2006).

    計算居里面深度的功率譜方法:將計算窗口內(nèi)地下深處的磁性體視為矩形棱柱體,則其磁異常頻譜F(u,v)為

    (1)

    式中,Ip為磁化強度,l、m、n為地磁場強度的方向余弦,L、M、N為總磁化強度的方向余弦,2a、2b分別為長和寬,x0,y0為磁性體中心水平坐標,u與v分別為沿x軸與y軸方向的圓頻率.

    令h1=h2-t=h2-αΔh,其中t為棱柱體的延深,α為大于1的正數(shù),則有

    (2)

    在a、b很小的時候,且在低頻區(qū)可以得到關(guān)于徑向頻率s的方程:

    LnH(s)-4Lns=A-2sh2.

    (3)

    本文采用重磁數(shù)據(jù)處理解釋軟件(RGIS)中的功率譜法反演居里面的深度.反演參數(shù):窗口寬度為8,滑動距離為4,低頻數(shù)為0.1,高頻數(shù)為0.5.

    3.1.2 計算莫霍面的Parker方法

    R.L.Parker提出了一種快速、精確的重磁界面正反演算法.因為它能計算物性橫向變化的連續(xù)界面,速度快,在國內(nèi)得到了廣泛的發(fā)展和應(yīng)用(Parker,1973;Oldenburg,1974;申寧華,1990;王萬銀和潘作樞,1993).

    設(shè)有地下物質(zhì)界面s,其上的剩余密度ρ為零,其下的剩余密度ρ不為零,平均深度為D,其在地面產(chǎn)生的重力異常為

    (4)

    (5)

    所以有

    (6)

    本文采用重磁數(shù)據(jù)處理解釋軟件(RGIS),反演密度界面.反演參數(shù):平均深度為55km,剩余密度為0.3g·m-3(考慮到下地殼與上地幔頂部的密度差).3.2 結(jié)果討論

    3.2.1 鮮水河斷裂帶

    鮮水河斷裂帶計算結(jié)果如圖2所示.其中,圖2a是沿斷裂分布的布格重力異常圖,圖2b是沿斷裂分布的航磁化極異常圖,圖2c是反演計算的莫霍面深度圖,圖2d是反演計算的居里面深度圖.

    在布格重力圖上(圖2a),鮮水河斷裂帶對應(yīng)一系列串珠狀重力異常.其西北端,道孚縣境內(nèi)斷層沿鮮水河延伸.中段,道孚至康定,斷裂帶走向向弧形改變,串珠狀重力異常的軸線也明顯隨斷層走向改變而變化.東南端,康定至石棉,斷裂帶沿大渡河延伸,重力異常受東部盆地影響,較為雜亂,近似梯級變化,斷裂帶以東為梯度變化較大的高值異常區(qū),斷裂帶以西為梯度變化較平緩的低值異常區(qū).航磁ΔT化極圖上(圖2b),鮮水河斷裂帶構(gòu)成磁場分界線.道孚以南,康定以北,斷裂帶為正負磁場分界線,

    圖2 鮮水河斷裂帶計算結(jié)果(a) 布格重力異常圖; (b) 航磁ΔT化極異常圖; (c) 莫霍面深度圖; (d) 居里面深度圖.Fig.2 The results in Xianshui River maps(a) Bouguer gravity anomaly map; (b) Aeromagnetic polarization map; (c) The depth of Moho surface map; (d) The depth of Curie isothermal surface map.

    圖3 理塘斷裂帶計算結(jié)果(a) 布格重力異常圖; (b) 航磁ΔT化極異常圖; (c) 莫霍面深度圖; (d) 居里面深度圖.Fig.3 The results in Litang maps(a) Bouguer gravity anomaly map; (b) Aeromagnetic polarization map; (c) The depth of Moho surface map; (d) The depth of Curie isothermal surface map.

    西側(cè)為平靜磁場,東側(cè)為雜亂串珠狀磁異常.康定以南,磁場明顯增強,在瀘定形成強磁區(qū),并被鮮水河斷裂帶分為東北高值強磁區(qū),西南低值強磁區(qū).

    反演得到的莫霍面深度圖中(2c),鮮水河斷裂帶對應(yīng)莫霍面深度梯度帶.莫霍面深度梯度帶的梯度方向與斷裂帶走向一致,深度由西北向東南逐漸降低,在道孚約53 km,在康定約48 km,在石棉約45 km.反演得到的居里面深度圖中(2d),鮮水河斷裂帶對應(yīng)一系列居里面高低起伏圈閉,表明地下19 km深度地溫沿鮮水河斷裂帶起伏變化明顯.總體上,沿鮮水河斷裂帶,居里面深度由西北向東南在19 km上下起伏,在道孚略微隆起,在康定以北略微隆起、在瀘定明顯隆起、在石棉明顯凹陷.

    區(qū)域地質(zhì)資料表明(四川省地質(zhì)礦產(chǎn)局,1991;裴錫瑜等,1985),鮮水河具有長期活動的歷史.鮮水河斷裂帶在早二疊世開始活動,形成海槽.晚二疊世至晚三疊世早-中期,形成兩套玄武巖建造并夾有復(fù)理石沉積.晚三疊中-晚期,由張裂轉(zhuǎn)為擠壓,兩側(cè)形成北西-北北西向褶皺帶.晚第三紀,鮮水河斷裂帶性質(zhì)發(fā)生轉(zhuǎn)變,形成北西-北北西向左行走滑斷層,水平錯距達15~76 km,控制了晚第三紀晚期至第四紀地震活動.現(xiàn)代階段,沿斷裂帶有頻繁的破壞性地震發(fā)生,1973年在爐霍發(fā)生地震的震級達7.9級,以及2014年11月22日發(fā)生康定地震,說明斷裂至今仍在活動.

    3.2.2 理塘斷裂帶

    理塘斷裂帶屬川西北次級塊體的內(nèi)部斷裂,呈近南北向延伸.理塘斷裂帶計算結(jié)果如圖3所示.其中,圖3a是沿斷裂分布的布格重力異常圖,圖3b是沿斷裂分布的航磁化極異常圖,圖3c是反演計算的莫霍面深度圖,圖3d是反演計算的居里面深度圖.

    布格重力異常圖上(圖3a),理塘斷裂帶以理塘縣為界,南北差異明顯.北段沿斷裂走向重力異常分布凌亂,規(guī)律性不強,只是在新龍縣,出現(xiàn)小的串珠狀重力異常沿斷裂走向分布.南段沿雅礱江延伸,在孟底溝鄉(xiāng)一帶,斷裂帶兩側(cè)重力高、中間重力低,形成沿斷裂帶走向的明顯重力“凹槽”.航磁ΔT化極圖上(圖3b),理塘斷裂帶也是一條明顯的磁場分界線.理塘縣以北,斷裂帶西側(cè)為強磁異常條帶區(qū),斷裂帶東側(cè)為低緩磁場區(qū),局地甚至無磁異常顯示.理塘縣以南,磁異常帶逐漸變寬,顯示出斷裂分叉特征,且斷裂帶西側(cè)的強磁異常條帶順斷層走向明顯加強,斷裂帶東側(cè)低緩磁場仍然較為平靜.

    反演得到的莫霍面深度圖中(圖3c),理塘斷裂帶南段孟底溝鄉(xiāng)一帶,斷裂帶與莫霍面深度等值線接近垂直相交,由木拉區(qū)南到孟底溝鄉(xiāng)南,莫霍面深度逐漸由52 km上升到45 km.理塘斷裂帶北段由理塘縣到新龍縣一帶,莫霍面形成局部微小的隆坳,在理塘縣深度約為55 km,新龍縣深度約為56 km.反演得到的居里面深度圖中(圖3d),總體上,沿理塘斷裂帶,居里面深度由西北向東南逐漸變深,西北端約18.5 km,東南端約19.5 km.但具體變化細節(jié),理塘斷裂帶南、北兩端特點仍然差別明顯.理塘縣以北,斷裂帶對應(yīng)居里面東深西淺的梯度帶,理塘縣以南,在木拉區(qū),斷裂帶穿過上下起伏的居里面,在孟底溝鄉(xiāng),斷裂帶則對應(yīng)居里面明顯下凹的“凹槽”.反演結(jié)果表明,理塘縣南北附近,居里面高低起伏較明顯,表明理塘斷裂帶在理塘縣附近深部地溫有較大起伏變化.

    區(qū)域地質(zhì)資料表明(四川省地質(zhì)礦產(chǎn)局,1991;裴錫瑜等,1985),理塘斷裂帶對應(yīng)甘孜—理塘縫合帶.沿理塘斷裂帶,蛇綠巖、蛇綠混雜巖、以及混雜堆積等十分發(fā)育.板塊構(gòu)造研究認為,理塘斷裂帶西側(cè)發(fā)育弧-盆體系,甘孜—理塘洋向西俯沖消減形成甘孜—理塘縫合帶.沿該斷裂帶分布的重磁異常,清楚地反應(yīng)了蛇綠巖、蛇綠混雜巖(基性—超基性巖)的巖漿弧特征.

    3.2.3 金沙江斷裂帶

    金沙江斷裂帶是川西高原重要的西邊界,呈近NS向分布.金沙江斷裂帶計算結(jié)果如圖4所示.其中,圖4a是沿斷裂分布的布格重力異常圖,圖4b是沿斷裂分布的航磁化極異常圖,圖4c是反演計算的莫霍面深度圖,圖4d是反演計算的居里面深度圖.

    布格重力圖上(圖4a),金沙江斷裂帶對應(yīng)重力高、低異常區(qū)轉(zhuǎn)換帶.其中,巴塘縣以北,斷裂帶兩側(cè)為高值區(qū),斷裂帶對應(yīng)重力低值條帶,這一特征在霍熱拉喀鄉(xiāng)一帶尤其明顯.巴塘縣以南,斷裂帶東側(cè)為一系列沿斷裂帶分布的串珠狀高值異常圈閉,斷裂帶西側(cè)為一系列沿斷裂帶分布的串珠狀低值異常圈閉,斷裂帶上則是高、低異常的過渡區(qū),這一特征在竹巴龍區(qū)以南尤其明顯.航磁ΔT化極圖上(圖4b),金沙江斷裂帶是一條不太明顯的磁場分界線.斷裂帶西側(cè)以大面積團塊狀正磁異常為主,東側(cè)以平緩的負磁場為主,零星分布一些橢圓狀、不規(guī)則狀局部磁異常.

    反演得到的莫霍面深度圖中(4c),金沙江斷裂帶北段,巴塘縣至白玉縣莫霍面局部隆坳相間,在55~57 km之間上下起伏.金沙江斷裂帶南段,巴塘縣至得榮縣,莫霍面逐漸抬升,由55 km上升至51 km,最后在得榮縣附近形成高臺圈閉.反演得到的居里面深度圖中(4d),總體上,金沙江斷裂帶穿過居里面由西向東加深的梯度帶上,斷裂帶東側(cè)居里面深度大多為18.5 km,西側(cè)居里面深度大多為19 km.沿金沙江斷裂帶,在巴塘縣南側(cè)的竹巴龍區(qū)、北側(cè)的霍熱拉喀鄉(xiāng),居里面形成局部的抬升高地,表明此處地溫出現(xiàn)高異常區(qū).

    圖4 金沙江斷裂帶計算結(jié)果(a) 布格重力異常圖; (b) 航磁ΔT化極異常圖; (c) 莫霍面深度圖; (d) 居里面深度圖.Fig.4 The results in Jinsha River maps(a) Bouguer gravity anomaly map; (b) Aeromagnetic polarization map; (c) The depth of Moho surface map; (d) The depth of Curie isothermal surface map.

    金沙江是長江上游干流,在川西流經(jīng)石渠、德格、白玉、巴塘、得榮五縣,為四川與西藏的界河.金沙江斷裂帶由一系列南北向斷裂組成,在重力、磁力異常圖和反演結(jié)果圖上,著名的金沙江蛇綠混雜巖群基本沒有顯示.說明該蛇綠混雜巖可能僅是一些小型推覆體,現(xiàn)有航磁和重力資料精度還難以反映這種小規(guī)模的構(gòu)造.而在巴塘縣附近,斜切金沙江構(gòu)造帶主體的巴塘斷裂在航磁和重力資料上反映卻十分清楚,說明該斷裂切割較深,影響較大.巴塘斷裂是一條走向約N30°E,傾向NW,傾角較陡的北東向斷裂;它起于莫多附近,向SW延伸經(jīng)巴塘、斜穿金沙江后繼續(xù)延伸至瀾滄江邊消失;巴塘斷裂生成時間晚于金沙江構(gòu)造帶,晚第四紀以來具有明顯的活動性,曾發(fā)生過1870年巴塘7級地震,該地震的地表破裂在一些地段上現(xiàn)今仍依稀可辨.

    3.2.4 巴塘—康定剖面

    以川西藏東地區(qū)的地殼P波速度結(jié)構(gòu)剖面(王椿墉等,2003)為參照依據(jù),我們沿北緯30°線,與人工地震探測二維剖面重合,計算了一條重磁聯(lián)合反演剖面,結(jié)果如圖5所示.

    該剖面始于金沙江附近的巴塘,向東經(jīng)理塘、雅江、康定至瀘定,全長320 km.反演結(jié)果如圖5所示.其中,圖5a是重力布格異常與反演擬合結(jié)果,圖5b是航磁異常與反演結(jié)果.圖5c是反演計算的剖面地質(zhì)結(jié)構(gòu),包括淺層巖性結(jié)構(gòu)、居里面和莫霍面.圖5d和圖5e是二維人工地震結(jié)果(王椿墉等,2003),是重磁反演的初始模型.

    巴塘—康定剖面跨越了金沙江斷裂,理塘斷裂和鮮水河斷裂.反演的密度和磁化強度按照實測物性資料給出(劉蓓莉,1993;四川省地質(zhì)礦產(chǎn)局,1991;裴錫瑜等,1985),片麻巖的密度為2.69 g·cm-3,磁化率為45×10-5SI,花崗巖的密度也為2.69 g·cm-3,磁化率為42×10-5SI,粉砂巖的密度為2.2 g·cm-3,磁化率為9×10-5SI,砂巖的密度為2.4 g·cm-3,磁化率為6×10-5SI,二疊系的平均密度為2.5 g·cm-3,平均磁化率28×10-5SI,三疊系的平均密度為2.45 g·cm-3,平均磁化率12×10-5SI.居里面與莫霍面之間的平均密度為2.90 g·cm-3,上地幔的平均密度為3.2 g·cm-3.反演的初始模型根據(jù)地震資料確定的上地殼結(jié)構(gòu)特征進行設(shè)計(王椿墉等,2003).

    圖5 典型剖面的反演解釋圖(a) 布格重力異常剖面圖; (b) 航磁異常剖面圖; (c) 地層結(jié)構(gòu)反演圖; (d) 地層結(jié)構(gòu); (e) 地殼人工地震剖面速度結(jié)構(gòu)(王椿墉等,2003),圖中數(shù)字單位為km·s-1.Fig.5 Inversion and interpretation of a typical profile maps(a) Bouguer gravity anomaly profile; (b) Aeromagnetic anomaly profile; (c) Strata structure inversion map; (d) Strata structure map; (e) Crustal artificial seismic velocity structure profile map(Wang Ch Y, 2003).

    反演結(jié)果表明(圖5c),剖面上,地殼厚度大致在60 km至49 km之間,莫霍面在巴塘以東和雅江以東呈上隆狀,兩地最淺埋深分別為50 km、49 km.而在理塘附近莫霍面達到最深,為60 km.這些與地震資料解釋結(jié)果(王椿墉等,2003)相近.反演計算的居里面起伏變化不大,大致由西向東逐漸抬升,由巴塘之下的22 km,到瀘定之下的17 km.地震資料解釋結(jié)果(王椿墉等,2003)為川西高原居里面在13~30 km深度之間,存在厚度8~10 km的低速層,速度5.8 km·s-1.低速層與居里面反演深度相符.地表淺層的巖性分布,以淺變質(zhì)的古生界及中生界三疊系地層為主,其下為二疊系片麻巖、砂巖.地層中有花崗巖體,分布在巴塘、理塘之間,以及康定附近.與地震資料(王椿墉等,2003)給出的初始模型(圖5d)形態(tài)一致,但深度和厚度略大.3.2.5 區(qū)域地質(zhì)背景

    利用Parker算法反演莫霍面深度、功率譜法反演居里面深度,是研究區(qū)域地質(zhì)背景的重要方法.我們的反演結(jié)果如圖6所示.其中,圖6a是莫霍面深度圖,圖6b是居里面深度圖.

    從莫霍面深度分布中(圖6a)可以看出,川西高原具有造山帶型地殼結(jié)構(gòu),地殼厚度介于43~63 km之間,莫霍面平均深度為55 km.從莫霍面等值線圖(圖6a)中可以看出,莫霍面基本特征表現(xiàn)為東南淺、西北深,反應(yīng)了該區(qū)地殼厚度逐漸加厚的趨勢.根據(jù)莫霍面的特征,可以將該區(qū)分為三個臺階,在丹巴—九龍以東,莫霍面深度較淺,為地殼厚度正常區(qū).道孚—雅江—稻城一帶,由東南至西北,莫霍面深度由淺變深,地殼厚度處于逐漸加厚的過渡區(qū).新龍—理塘—芒康以西,莫霍面深度進一步變深,為地殼厚度增厚區(qū).

    從居里面深度分布(圖6b)中可以看出,川西居里面深度介于17~23 km之間,平均深度為20 km,與前人(袁學(xué)誠,1996;張昌達,2003)對該地區(qū)的反演結(jié)果相近.居里面分布具有明顯的條帶狀特征,沿著丹巴—瀘定—冕寧以東,居里面深度較淺.雅礱江東西兩側(cè),居里面深度特征截然不同,雅礱江以東至貢嘎山之間,居里面深度較大.雅礱江以西,居里面的深度變淺.金沙江南北兩段,居里面深度特征也截然不同.金沙江北段,巴塘至白玉之間,東側(cè)居里面深度較小,西側(cè)居里面深度較大.金沙江南段,巴塘至得榮之間,東側(cè)居里面深度較大,西側(cè)居里面深度較小.

    利用本文反演結(jié)果,結(jié)合前人研究成果(裴錫瑜等,1985;熊熊和滕吉文,2002;盧占武等,2006),本文給出川西高原區(qū)域地質(zhì)背景(圖7).該圖以莫霍面、居里面反演結(jié)果為依據(jù),利用莫霍面的物質(zhì)成分屬性、居里面的物理溫度屬性,構(gòu)造深部地質(zhì)背景.莫霍面是因為物質(zhì)成分組成的差別形成的速度不連續(xù)面,居里面是因為介質(zhì)磁性隨溫度變化而與磁效應(yīng)相關(guān)的等溫面.這兩個互不相干的面,卻都與深部地質(zhì)背景密切相關(guān).圖7中,莫霍面構(gòu)成了川西高原的地幔與地殼的深部形態(tài),居里面則給出了地殼一定深度的可能溫度分布.居里面的確定主要依據(jù)兩個數(shù)據(jù),即磁鐵礦的居里溫度580 ℃、地殼地溫梯度25~30 ℃·km-1.由于川西高原地殼內(nèi)的這兩個數(shù)值并不容易確定,所以我們以400~600 ℃確定居里面的分布范圍,如圖7中所示.

    從圖7可知,川西高原莫霍面深度在43~63 km之間,居里面深度在17~23 km之間,與前人研究的地殼結(jié)構(gòu)特征相符.前人資料表明(裴錫瑜等,1985;熊熊和滕吉文,2002;盧占武等,2006),川西高原地殼在20±5 km深度范圍內(nèi)普遍存在低速高導(dǎo)層,速度一般為5.6~5.8 km·s-1,電阻率約為1~10 Ωm,厚度一般為5~10 km.雖然低速層與高導(dǎo)層位置不完全重合,但深度和厚度變化趨勢一致.居里面給出的結(jié)果表明殼內(nèi)低速高導(dǎo)層可能成因是地殼內(nèi)部分熔融.圖7中,川西高原西界的金沙江一線,與居里面相關(guān)的深部溫度出現(xiàn)兩個異常增大區(qū).這既可能導(dǎo)致沿金沙江斷裂帶的地震活動,也會導(dǎo)致巖漿活動和地?zé)釡厝c的廣泛分布.此外,地殼內(nèi)的弱物質(zhì)將沿邊界斷層發(fā)生塑性擠出,從而形成川西高原現(xiàn)今的構(gòu)造格局.

    圖6 川西高原地區(qū)重磁數(shù)據(jù)反演結(jié)果圖(a) 莫霍面深度分布; (b) 居里面深度分布.Fig.6 Western Sichuan Plateau gravity and magnetic datum inversion images(a) The distribution of the Moho depth; (b) The distribution of the Curie depth.

    圖7 地質(zhì)模型(圖中數(shù)字為溫度,單位℃)Fig.7 Geological model

    整體上,川西高原發(fā)生的絕大多數(shù)地震震源深度在5~15 km之間,居里面深度層與地震活動特征相符.居里面之上,為脆性地殼,屬于脆性介質(zhì)的上地殼范圍內(nèi).斷裂帶上的地震應(yīng)力在脆性地殼內(nèi)集中,這是強地震發(fā)生的深部介質(zhì)條件.上地殼低速、低阻層,與居里面推斷的高溫深層韌性剪切帶一致.逆沖斷裂帶使堅硬揚子塊體與軟弱的川西塊體相撞,產(chǎn)生強烈變形,使得下地殼相對增厚,并在地表形成不同性質(zhì)的斷裂帶.圖7的結(jié)果表明,川西高原的地殼結(jié)構(gòu)特征主要表現(xiàn)在下地殼的增厚、中地殼高溫弱物質(zhì)流的塑性變形、上地殼的塊體走滑變形.

    4 結(jié)論

    川西高原是中國大陸的主要地震活動區(qū)之一.巴塘地震區(qū)、理塘地震區(qū)和康定地震區(qū)分別與金沙江斷裂帶、理塘斷裂帶和鮮水河斷裂帶相聯(lián)系.川西高原發(fā)生的絕大多數(shù)地震震源深度在上地殼內(nèi),中深部地殼結(jié)構(gòu)明顯與居里面起伏反映的熱異常過程相關(guān).通過對川西重磁異常分析研究,主要結(jié)論如下:

    (1)川西高原布格重力異常特征是西北低、東南高,多為串珠狀異常.磁場特征表現(xiàn)為明顯的分帶性,從西向東以甘孜—理塘為界,可劃分為白玉—德欽劇烈變化磁場區(qū)、雅江—九龍平緩磁場區(qū)、瀘定—冕寧強磁場區(qū).

    (2)重力反演結(jié)果表明,該區(qū)莫霍面東南淺、西北深,反映川西高原地殼厚度由東南向西北逐漸加厚,地殼厚度在43~63 km之間,平均厚度55 km.莫霍面變化具有階梯式變化的特點,從西南向東北依次加深.

    (3)居里面結(jié)果表明,該區(qū)居里面特征表現(xiàn)為條帶狀,深度在17~23 km之間,平均深度為20 km.丹巴—瀘定—冕寧以東,居里面深度較淺.雅礱江東側(cè)居里面加深、西側(cè)居里面變淺.金沙江南段,東側(cè)居里面深度較大,西側(cè)居里面深度較小.金沙江北段,東側(cè)居里面深度較小,西側(cè)居里面深度較大.

    (4)巴塘—瀘定剖面的二維反演結(jié)果表明,不同構(gòu)造單元地殼結(jié)構(gòu)縱向分布特征不同.巴塘以東50 km處,莫霍面上隆,地殼最淺處只有50 km.理塘之下,莫霍面下凹,地殼最深處為60 km.居里面起伏變化不大,大致由西向東逐漸抬升,由巴塘之下的22 km,到瀘定之下的17 km.

    (5)總體上,川西高原居里面呈NS向條帶狀分布特征,與該地區(qū)斷裂的走向特征一致,斷裂兩側(cè)的居里面起伏相間.居里面位于莫霍面之上,莫霍面的深度從東向西逐漸加深,但是居里面深度變化不大.

    (6)川西高原深部呈現(xiàn)下地殼以增厚為主,中地殼以熱塑性變形為主,上地殼以塊體走滑變形為主的地殼結(jié)構(gòu)特點.

    Bhattacharyya B K, Leu L K. 1975. Spectral analysis of gravity and magnetic anomalies due to two-dimensional structures.Geophysics, 40(6): 933-1031.Bhattacharyya B K. 1966. Continuous spectrum of the total-magnetic-field anomaly due to a rectangular prismatic body.Geophysics, 31(1): 97-121.

    Blakely R J. 1988. Curie temperature isotherm analysis and tectonic implications of aeromagnetic data from Nevada.JournalofGeophysicalResearch, 93(B10): 11817-11832.

    Cheng J, Liu J, Gan W J, et al. 2009. Influence of coseismic deformation of the Wenchuan earthquake on the occurrence of earthquakes on active faults in Sichuan-Yunnan region.ActaSeismologicaSinica, 31(5): 477-490.

    Cui Z X, Pei S P. 2009. Study on Pn velocity and anisotropy in the uppermost mantle of the eastern Himalayan syntaxes and surrounding regions.ChineseJ.Geophys.(in Chinese), 52(9): 2245-2254, doi: 10.3969/j.issn.0001-5733.2009.09.008.

    Ding Z F, He Z Q, Sun W G, et al. 1999. 3-D crust and upper mantle velocity structure in eastern Tibetan Plateau and its surrounding areas.ChineseJ.Geophys.(in Chinese), 42(2):197-205.

    Han W B, Jiang G F. 2004. Study on distribution characteristics of strong earthquakes in Sichuan-Yunnan area and their geological tectonic background.ActaSeismologicaSinica, 26(2): 211-222.

    Hao S J, Gao H G, Wang C H. 1982. Isotheral surface of Curie temperature around the Beijing area and its relation to earthquakes.ChineseJ.Geophys.(in Chinese), 25(3): 264-269.

    Hou Z C, Li B G. 1985b. Power spectrum methods for direct estimating the depth to the bottom of magnetic source.ComputingTechniquesforGeophysicalandGeochemicalExploration, 7(3): 179-187.

    Huang J L, Zhao D P, Zheng S H. 2001. Seismic tomography of the Sichuan-Yunan active tectonic region.ChineseJ.Geophys.(in Chinese), 44(s1): 127-135. doi: 10.3321/j.issn:0001-5733.2001.z1.016.

    Li B L, Lei Y L, Chen Y X, et al. 2011. The neo-tectonic deformation features in the eastern segment of Circum-Tibetan Plateau Basin and Range System: Take the western Sichuan basin as an example.ActaPetrologicaSinica, 27(3): 636-644.

    Li F, Zhou S Y, Su Y G, et al. 2011. Study on Pn-wave velocity structure and anisotropy in the Sichuan-Yunnan and its adjacent areas.ChineseJ.Geophys.(in Chinese), 2011, 54(1): 44-54, doi: 10.3969/j.issn.0001-5733.2011.01.006.

    Liu B L. 1993. Analysis and application of density data of rocks from Sichuan Province.GeophysialandGeochemaicalExploration, 18(3):232-237.Liu Y C, Wu D C, Wang D Y, et al. 2010. The structural deformation characteristics of western Sichuan Luding fault.JournalofChengduUniversityofTechnology(Science&TechnologyEdition), 37(3): 262-267.

    Lou H, Wang C Y. 2005. Wavelet analysis and interpretation of gravity data in Sichuan-Yunnan region, China.ActaSeismologicalSinica, 27(5):515-523.Lu Z W, Gao R, Li Q S, et al. 2006. Deep geophysical probe and geodynamic study on the Qinghai-Tibet P1ateau (1958—2004).ChineseJ.Geophys.(in Chinese), 49(3):753-770.

    Luo G, Zhao C P, Zhou L Q. 2014. Characteristics of focal mechanisms and stress field of the Chuan-Dian rhombic block and its adjacent regions.SeismologyandGeology, 26(2): 405-421, doi:10.3969 /j.issn.0253-4967.2014.02.011.

    Oldenburg D W. 1974. The inversion and interpretation of gravity anomalies.Geophysics, 39(4): 526-536.

    Parker R L. 1973. The rapid calculation of potential anomalies.GeophysicsJournalInternationa, 31(4): 447-455.

    Pei X Y, Chen W D, Zhang J T, et al. 1985. A preliminary study on the use of gravity and magnetic data deep crustal structure and its relationship with the earthquake in Sichuan Province in the western region.SichuanEarthquake, (2): 16-22.

    Qiao X J, Wang Q, Du R L. 2004. Characteristics of current crustal deformation of active blocks in the Sichuan-Yunnan region.ChineseJ.Geophys.(in Chinese), 47(5): 805-811.

    Shen N H. 1985. Introduction to principleand methods of Curie point calculation by aeromagnetic data.ComputingTechniquesforGeophysicalandGeochemicalExploration, 7(2):89-98.

    Shen N H. 1990. The rapid caculation of potential anomalies and its application.ComputingTechniquesforGeophysicalandGeochemicalExploration, 12(1):4-12.Sichuan Provincial Bureau of Geology and Mineral Resources. 1991. Regional Geology of Sichuan. Beijing:Geological Publishing House.

    Spector A, Grant F S. 1970. Statistical models for interpreting aeromagnetic data.Geophysics, 35(2): 293-302.

    Tanaka A, Okubo Y, Matsubayashi O. 1999. Curie point depth based on spectrum analysis of the magnetic anomaly data in East and Southeast Asia.Tectonophysics, 306(3-4): 461-470.

    Tang W Q, Chen Z L, Liu Y P, et al. 2005. Present-day tectonics activity in the intersection area of the Xianshuihe fault and Longmenshan fault on the eastern margin of the Qinghai-Tibet Plateau.GeologicalBulletinofChina, 24(12): 1169-1172.

    Wang C Y, Wu J P, Lou H, et al. 2003. P-wave crustal velocity structure of in western Sichuan and eastern Tibetan region.ScienceinChinaSeriesD-EarthScience, 46(Suppl.): 254-265.

    Wang W Y, Pan Z S. 1993. Fast solution of forward and inverse problems for gravity field in a dual interface model.GeophysicalProspectingforPetroleum, 32(2):81-87.

    Wang Y Z, Wang E N, Chen Z K, et al. 2008. GPS-constrained inversion of present-day slip rates along major faults of the Sichuan-Yunnan region, China.ScienceinChinaSeriesD:EarthSciences, 51(9): 1267-1283.

    Wen X Z. 2000. Character of rupture segmentation of the Xianshuihe-Anninghe-Zemuhe fault zone, western Sichuan.SeismologyandGeology, 22(3): 239-249.

    Wu J P, Ming Y H, Wang C Y. 2006. Regional waveform inversion for crustal and upper mantle velocity structure below Chuandian region.ChineseJ.Geophys.(in Chinese), 49(5): 1369-1376.

    Xiong X, Teng J W. 2002. Study on crustal movement and deep process in eastern Qinghai-Xizang plateau.ChineseJ.Geophys.(in Chinese), 45(4):507-515.

    Xue X W, Wen X Z, Cheng R Z, et al. 2003b. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region.ScienceinChinaSeriesD:EarthSciences, 51(9): 210-220.

    Xü C F, Pan J S, Wang F Y, et al. 2008. Deep geophysical exploration of Longmenshan and its adjacent area.JournalofGeodesyandGeodynamics, 28(6):31-37.

    Xv X W, Zhang P Z, Wen X Z, et al. 2005. Features of active tectonics and recurrence behaviors of strong earthquakes in the western Sichuan Province and its adjacent regions.SeismologyandGeology, 27(3): 446-461.

    Yuan X C. 1996. Geophysics of Atlas of China. Beijing:Geological Publishing House.

    Zhang C D. 2003. Deduction of magnetic characteristics of lithosphere in China from results on satellite and aeromagnetic measurements.ProgressinGeophysics, 18(1): 103-110.

    Zhang Y, Cheng Sh Y, Zhao B K, et al. 2013. The feature of tectonics in the Tibet Plateau from new egional gravity signals.ChineseJ.Geophys.(in Chinese),56(4): 1369-1380, doi: 10.6038/cjg20130431.

    Zhao B M, Hao T Y. 2006. The geological significance and the inversion method of the magnetic geologic interface.ProgressinGeophysics, 21(2): 353-359.

    Zhou R J, Chen G X, Li Y, et al. 2005. Research on active faults in Litang-Batang region, western Sichuan Province, and the seism-genic structures of the 1989 BatangM6.7 earthquake swarm.SeismologyandGeology, 27(1):31-43.

    附中文參考文獻

    程佳, 劉杰, 甘衛(wèi)軍等. 2009. 汶川地震同震形變場對川滇地區(qū)主要活動斷裂地震發(fā)生趨勢的影響. 地震學(xué)報, 31(5): 477-490.

    崔仲雄, 裴順平. 2009. 青藏高原東構(gòu)造結(jié)及周邊地區(qū)上地幔頂部Pn速度結(jié)構(gòu)和各向異性研究. 地球物理學(xué)報, 52(9): 2245-2254, doi: 10.3969/j.issn.0001-5733.2009.09.008.

    丁志峰, 何正勤, 孫為國等. 1999. 青藏高原東部及其邊緣地區(qū)的地殼上地幔三維速度結(jié)構(gòu). 地球物理學(xué)報, 42(2): 197-205.

    韓渭賓, 蔣國芳. 2004. 川滇地區(qū)強震活動分布特征及其與地殼塊體構(gòu)造背景關(guān)系的研究. 地震學(xué)報, 26(2): 211-222.

    郝書儉, 高華根, 王春華. 1982. 京津唐地區(qū)居里等溫面及其與地震的關(guān)系. 地球物理學(xué)報, 25(3): 264-269.

    侯重初, 李保國. 1985b. 直接計算磁性下界面深度的功率譜法. 物化探計算技術(shù), 7(3): 179-187.

    黃金莉, 趙大鵬, 鄭斯華. 2001. 川滇活動構(gòu)造區(qū)地震層析成像. 地球物理學(xué)報, 44(s1), 127-135, doi: 10.3321/j.issn:0001-5733.2001.z1.016.

    李本亮, 雷永良, 陳竹新等. 2011. 環(huán)青藏高原盆山體系東段新構(gòu)造變形特征——以川西為例. 巖石學(xué)報, 27(3): 636-644.

    李飛, 周仕勇, 蘇有錦等. 2011. 川滇及鄰區(qū)Pn波速度結(jié)構(gòu)和各向異性研究. 地球物理學(xué)報, 54(1): 44-54, doi: 10.3969/j.issn.0001-5733.2011.01.006.

    劉蓓莉. 1993. 四川省巖石密度數(shù)據(jù)的分析及應(yīng)用. 物探與化探, 18(3): 232-237.

    劉援朝, 吳德超, 王道永等. 2010. 四川西部瀘定斷裂帶構(gòu)造變形特征. 成都理工大學(xué)學(xué)報(自然科學(xué)版), 37(3): 262-267.

    樓海, 王椿鏞. 2005. 川滇地區(qū)重力異常的小波分解與解釋. 地震學(xué)報, 27(5): 515-523.

    盧占武, 高銳, 李秋生等. 2006. 中國青藏高原深部地球物理探測與地球動力學(xué)研究(1958-2004). 地球物理學(xué)報, 49(3): 753-770.

    羅鈞, 趙翠萍, 周連慶. 2014. 川滇塊體及周邊區(qū)域現(xiàn)今震源機制和應(yīng)力場特征. 地震地質(zhì), 6(2): 405-421, doi:10.3969 /j.issn.0253-4967.2014.02.011.

    裴錫瑜, 陳文德, 張家濤等. 1985. 利用重磁資料對四川西部地區(qū)地殼深部構(gòu)造及其與地震關(guān)系的初步研究. 四川地震, (2): 16-22.

    喬學(xué)軍, 王琪, 杜瑞林. 2004. 川滇地區(qū)活動地塊現(xiàn)今地殼形變特征. 地球物理學(xué)報, 47(5): 805-811.

    申寧華. 1985. 用航磁數(shù)據(jù)計算居里點深度的原理及方法. 物化探計算技術(shù), 7(2): 89-98.

    申寧華. 1990. 界面位場異常的快速正反演計算及其應(yīng)用. 物化探計算技術(shù), 12(1): 4-12.

    四川省地質(zhì)礦產(chǎn)局. 1991. 四川省區(qū)域地質(zhì)志. 北京: 地質(zhì)出版社.

    唐文清, 陳智梁, 劉宇平等. 2005. 青藏高原東緣鮮水河斷裂與龍門山斷裂交會區(qū)現(xiàn)今的構(gòu)造活動. 地質(zhì)通報, 24(12): 1169-1172.

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

    王慶良, 崔篤信, 王文萍等. 2008. 川西地區(qū)現(xiàn)今垂直地殼運動研究. 中國科學(xué)(D輯): 地球科學(xué), 38(5): 598-610.

    王萬銀, 潘作樞. 1993. 雙界面模型重力場快速正反演問題. 石油物探, 32(2): 81-87.

    王閻昭, 王恩寧, 沈正康等. 2008. 基于GPS資料約束反演川滇地區(qū)主要斷裂現(xiàn)今活動速率. 中國科學(xué) D 輯: 地球科學(xué), 38(5): 582-597.

    聞學(xué)澤. 2000. 四川西部鮮水河—安寧河—則木河斷裂帶的地震破裂分段特征. 地震地質(zhì), 22(3): 239-249.

    吳建平, 明躍紅, 王椿鏞. 2006. 川滇地區(qū)速度結(jié)構(gòu)的區(qū)域地震波形反演研究. 地球物理學(xué)報, 49(5): 1369-1376.

    熊熊, 滕吉文. 2002. 青藏高原東緣地殼運動與深部過程的研究. 地球物理學(xué)報, 45(4): 507-515.

    徐朝繁, 潘紀順, 王夫運等. 2008. 龍門山及其鄰近地區(qū)深部地球物理探測. 大地測量與地球動力學(xué), 28(6): 31-37.

    徐錫偉, 聞學(xué)澤, 鄭榮章等. 2003b. 川滇地區(qū)活動塊體最新構(gòu)造變動樣式及其動力來源. 中國科學(xué)(D輯), 33(s1): 151-162.

    徐錫偉, 張培震, 聞學(xué)澤等. 2005. 川西及其鄰近地區(qū)活動構(gòu)造基本特征與強震復(fù)發(fā)模型. 地震地質(zhì), 27(3): 446-461.

    袁學(xué)誠. 1996. 中國地球物理圖集. 北京: 地質(zhì)出版社.

    張昌達. 2003. 由衛(wèi)星和航空磁測成果推斷中國巖石圈的磁性特征. 地球物理學(xué)進展, 18(1): 103-110.

    張培震. 2008. 青藏高原東緣川西地區(qū)的現(xiàn)今構(gòu)造變形、應(yīng)變分配與深部動力過程. 中國科學(xué)(D輯): 地球科學(xué), 38(9): 1041-1056.

    張燕, 程順有, 趙炳坤等. 2013. 青藏高原構(gòu)造結(jié)構(gòu)特點: 新重力異常成果的啟示. 地球物理學(xué)報, 56(4): 1369-1380, doi: 10.6038/cjg20130431.

    趙百民, 郝天珧. 2006. 反演磁性地質(zhì)界面的意義與方法. 地球物理學(xué)進展, 21(2): 353-359.

    周榮軍, 陳國星, 李勇等. 2005. 四川西部理塘—巴塘地區(qū)的活動斷裂與1989年巴塘6.7級震群發(fā)震構(gòu)造研究. 地震地質(zhì), 27(1): 31-43.

    (本文編輯 胡素芳)

    The study of gravity-magnetic anomaly and tectonic background in Sichuan west region

    GAO Ling-Ju1,2, ZHANG Jian1,2*, DONG Miao1,2

    1KeyLaboratoryofComputationalGeodynamics,ChineseAcademyofScience,Beijing100049,China2UniversityofChineseAcademyofScience,Beijing100049,China

    The Western Sichuan Plateau, the eastern margin of the Tibetan Plateau, is the important part of the Sichuan-Yunnan rhombus block and the most intense tectonic deformation of the continental crust and high-seismicity region. The Batang earthquake zone, the Litang earthquake zone and the Kangding earthquake zone are related to the Jinsha River fault, the Litang fault and the Xianshui River fault respectively. Most of focal depth distributed in the upper crust of the Western Sichuan Plateau. The middle and deep crustal structure obviously related to the thermal abnormal process which revealed from the Curie isothermal surface undulation.The characteristics of gravity and magnetic anomalies, the Moho depth and the Curie isothermal depth in the Xianshui River fault, the Litang fault and the Jinsha River fault were studied. The analytical and inversion methods were used to the latest gravity and aeromagnetic data in the paper. The characteristic of Bouguer gravity anomalies of the Western Sichuan Plateau are higher in the southeast than in the northwest, showing a beads shaped. The band shaped magnetic field can be divided into three districts from west to east with Ganzi-Litang boundary: the strong changing district in Baiyu-Deqin, the quiet district in Ya River-Jiulong, the strong district in Luding-Mianning. The Moho depth is shallower in the southeast and characterized with chair-step shape from the gravity inversion results. The crust thickness in the Western Sichuan Plateau varies from 43 to 63 km and the average thickness is 55 km. The Curie isothermal surface depth ranges from 17 km to 23 km and with an average of 20 km. Xianshui River fault corresponding to the Moho surface gradient zone and the amplitude of the Curie isothermal surface trap. The Moho surface characterized with uplift and depression in the northern of Litang River fault, and uplift gradually in the southern. The Curie isothermal surface deepens from east to west. Along Jinsha River fault, the Curie isothermal surface shows local uplift, where the heat geothermal anomalies source maybe existing under this area. The two-dimensional inversion results of the Batang-Luding profile indicate the vertical crustal structure variant among different tectonic units. The Moho surface shows uplift and the minimum crust thickness is only 50 km in the 50 km to the east of Batang. The Moho depth below Litang is 60 km with a concave trend. The Curie isothermal surface shows weak variation and uplift gradually from west to east between 22 km under Batang and 17 km under Luding. In general, the Curie isothermal surface in the Western Sichuan Plateau is in zoning distribution from north to south, which in accordance with the fault strike. The Curie isothermal surface shows ups and downs on either side of the fracture. The Moho depth gradually deepens from west to east, but the Curie isothermal surface above the Moho surface shows weak variation. We concluded that the Western Sichuan Plateau can be characterized with lower crust thickening, mid-crust thermal deforming and upper crust brittle deforming from comprehensive analytical approach.

    The Western Sichuan Plateau; The Sichuan-Yunnan rhombus block;Crustal structure; Gravity and magnetic inversion

    國家自然科學(xué)基金項目(41574074,41174085,41430319),中國科學(xué)院創(chuàng)新團隊項目(KZZD-EW-TZ-19),中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(XDA1103010102)聯(lián)合資助;中國地調(diào)局項目(12120113101400)資助.

    高玲舉,1988年生,女,博士生,地球物理專業(yè).E-mail:gaolingju88@126.com

    *通訊作者張健,1963年生,男,教授,博士生導(dǎo)師,地球物理專業(yè).E-mail:zhangjian@ucas.ac.cn

    10.6038/cjg20150831.Gao L J, Zhang J, Dong M. 2015. The study of gravity-magnetic anomaly and tectonic background in Sichuan west region.ChineseJ.Geophys. (in Chinese),58(8):2996-3008,doi:10.6038/cjg20150831.

    10.6038/cjg20150831

    P631

    2015-01-19,2015-06-17收修定稿

    高玲舉, 張健, 董淼.2015.川西高原重磁異常特征與構(gòu)造背景分析.地球物理學(xué)報,58(8):2996-3008,

    猜你喜歡
    鮮水河川西高原理塘
    理塘札記
    基于PRAC法則的理塘抖音短視頻營銷策略研究
    中國市場(2022年14期)2022-06-07 06:01:05
    從“丁真走紅”談脫貧攻堅新路徑
    汶川MS8.0大震前后鮮水河斷裂帶滑動及其西部地殼形變的動態(tài)變化
    鮮水河斷裂帶中南段同位素年代學(xué)及其地質(zhì)意義
    川西甘孜-理塘結(jié)合帶碳酸鹽巖碳氧同位素特征
    川西高原東部兩次連續(xù)強降雪過程對比分析
    川西高原一次持續(xù)性暖區(qū)強降水分析
    約旦野生二棱大麥在川西高原的物候期和農(nóng)藝性狀分析
    鮮水河斷裂的幾何形態(tài)對地震發(fā)生的影響1
    av线在线观看网站| 国产成人午夜福利电影在线观看| 国产白丝娇喘喷水9色精品| 日韩av在线免费看完整版不卡| 国内揄拍国产精品人妻在线| 日韩成人av中文字幕在线观看| 国产精品一区www在线观看| 中文在线观看免费www的网站| 国产成人freesex在线| 精品人妻偷拍中文字幕| 丰满乱子伦码专区| 热99re8久久精品国产| 免费观看精品视频网站| 一区二区三区高清视频在线| 精品久久久久久久久av| 国产亚洲av嫩草精品影院| 亚洲美女搞黄在线观看| 三级经典国产精品| 国产成人91sexporn| 久久6这里有精品| 免费看av在线观看网站| av在线亚洲专区| 在线观看av片永久免费下载| 综合色av麻豆| 久久久久性生活片| 国产美女午夜福利| 成人亚洲精品av一区二区| 亚洲性久久影院| 日韩av不卡免费在线播放| 国产高清视频在线观看网站| av女优亚洲男人天堂| av在线老鸭窝| 午夜福利网站1000一区二区三区| 日韩欧美三级三区| 久久午夜福利片| 午夜福利成人在线免费观看| .国产精品久久| 欧美bdsm另类| 18禁在线无遮挡免费观看视频| 国产v大片淫在线免费观看| 尾随美女入室| 亚洲乱码一区二区免费版| 久久久精品欧美日韩精品| 一级黄片播放器| 日韩一区二区视频免费看| 人妻夜夜爽99麻豆av| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品成人久久小说| kizo精华| 亚洲成人中文字幕在线播放| 嫩草影院入口| 国产高清三级在线| 国内少妇人妻偷人精品xxx网站| 欧美精品国产亚洲| 午夜激情欧美在线| 一个人看视频在线观看www免费| 村上凉子中文字幕在线| 视频中文字幕在线观看| 九九爱精品视频在线观看| 黄色一级大片看看| 国产精品一及| 一本久久精品| 精品国内亚洲2022精品成人| 2022亚洲国产成人精品| 国产黄片视频在线免费观看| 亚洲国产成人一精品久久久| 欧美日韩综合久久久久久| 亚洲国产精品专区欧美| 中文乱码字字幕精品一区二区三区 | 国产中年淑女户外野战色| 99久久精品一区二区三区| 中文天堂在线官网| 深爱激情五月婷婷| 内地一区二区视频在线| av.在线天堂| 亚洲精品乱久久久久久| 国产黄色视频一区二区在线观看 | 男女国产视频网站| 乱人视频在线观看| 亚洲色图av天堂| 国产精品野战在线观看| 久久久欧美国产精品| av专区在线播放| 久久久久久久久久成人| 日韩国内少妇激情av| 岛国毛片在线播放| 亚洲中文字幕一区二区三区有码在线看| 精品人妻偷拍中文字幕| 欧美bdsm另类| 日韩三级伦理在线观看| 国产亚洲最大av| 99国产精品一区二区蜜桃av| 男人的好看免费观看在线视频| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久久久久久久| 小说图片视频综合网站| 日韩三级伦理在线观看| 三级国产精品片| av在线亚洲专区| 免费观看a级毛片全部| 一区二区三区四区激情视频| 丝袜喷水一区| 神马国产精品三级电影在线观看| 亚洲国产精品成人综合色| 直男gayav资源| 高清视频免费观看一区二区 | 久久久精品94久久精品| av在线播放精品| 69av精品久久久久久| 18禁裸乳无遮挡免费网站照片| 水蜜桃什么品种好| 婷婷色av中文字幕| 亚洲欧美日韩高清专用| 午夜爱爱视频在线播放| 国产乱来视频区| 亚洲人成网站在线播| 天堂网av新在线| 91狼人影院| 国产精品国产高清国产av| 日韩欧美 国产精品| 欧美bdsm另类| 国产精品三级大全| 99国产精品一区二区蜜桃av| 国产精品嫩草影院av在线观看| 青春草国产在线视频| 特大巨黑吊av在线直播| 看黄色毛片网站| 久久久亚洲精品成人影院| 我要搜黄色片| 国产 一区 欧美 日韩| 人妻夜夜爽99麻豆av| 国产视频首页在线观看| 国产乱人偷精品视频| 欧美三级亚洲精品| 美女xxoo啪啪120秒动态图| 午夜精品国产一区二区电影 | 直男gayav资源| 亚洲中文字幕一区二区三区有码在线看| 久久人人爽人人爽人人片va| 成人二区视频| 国产极品精品免费视频能看的| 色哟哟·www| a级毛色黄片| 天天一区二区日本电影三级| 国产成年人精品一区二区| 午夜激情欧美在线| 在线a可以看的网站| 别揉我奶头 嗯啊视频| 精品国产一区二区三区久久久樱花 | av国产久精品久网站免费入址| 亚洲在线观看片| 波多野结衣高清无吗| 69人妻影院| 亚洲图色成人| 亚洲av中文av极速乱| 搡老妇女老女人老熟妇| 国产av一区在线观看免费| 国内精品美女久久久久久| 久久久久久久午夜电影| 99热这里只有精品一区| 午夜福利高清视频| 男女视频在线观看网站免费| 人人妻人人看人人澡| 91av网一区二区| av卡一久久| 99热精品在线国产| 久久久久精品久久久久真实原创| 日本一二三区视频观看| 国产伦一二天堂av在线观看| 亚洲国产欧美在线一区| 国产v大片淫在线免费观看| 国产av码专区亚洲av| 97在线视频观看| 中文字幕制服av| 亚洲电影在线观看av| 一级毛片电影观看 | 日日啪夜夜撸| 国产在视频线精品| 国产成人freesex在线| 丰满人妻一区二区三区视频av| 在线免费观看不下载黄p国产| 少妇熟女aⅴ在线视频| 色吧在线观看| 亚洲精品456在线播放app| 美女脱内裤让男人舔精品视频| 久久久久久国产a免费观看| 蜜桃亚洲精品一区二区三区| 一区二区三区高清视频在线| 午夜福利在线观看吧| 麻豆一二三区av精品| www日本黄色视频网| 免费观看的影片在线观看| 国产精品人妻久久久久久| 久久这里只有精品中国| 国产精品国产三级国产专区5o | 国产精品美女特级片免费视频播放器| 国产探花在线观看一区二区| 蜜桃亚洲精品一区二区三区| 性插视频无遮挡在线免费观看| 中文乱码字字幕精品一区二区三区 | 亚洲自偷自拍三级| 午夜福利在线观看吧| 老司机福利观看| 国产精品不卡视频一区二区| 国产精品综合久久久久久久免费| 中文欧美无线码| 床上黄色一级片| 毛片女人毛片| 久久久久久久亚洲中文字幕| 亚洲国产高清在线一区二区三| 亚洲av成人精品一二三区| 少妇猛男粗大的猛烈进出视频 | 乱系列少妇在线播放| 国产黄色视频一区二区在线观看 | 久久久久精品久久久久真实原创| 观看美女的网站| 韩国av在线不卡| 午夜视频国产福利| 一级黄色大片毛片| 亚洲欧洲国产日韩| 91在线精品国自产拍蜜月| 水蜜桃什么品种好| 国产精品国产三级国产专区5o | 国产亚洲午夜精品一区二区久久 | 国产高清国产精品国产三级 | 91久久精品国产一区二区三区| 亚洲国产成人一精品久久久| 日韩av在线免费看完整版不卡| videossex国产| 97超碰精品成人国产| 日韩大片免费观看网站 | 卡戴珊不雅视频在线播放| 日本免费a在线| 国产精品久久久久久久久免| 欧美日韩国产亚洲二区| 亚洲精品成人久久久久久| 毛片一级片免费看久久久久| 热99re8久久精品国产| 中文在线观看免费www的网站| 大又大粗又爽又黄少妇毛片口| 美女脱内裤让男人舔精品视频| 日本色播在线视频| 在线观看美女被高潮喷水网站| 狂野欧美白嫩少妇大欣赏| 看黄色毛片网站| 亚洲精品自拍成人| 久久久精品94久久精品| 国产精品美女特级片免费视频播放器| 亚洲三级黄色毛片| 日本免费一区二区三区高清不卡| 国产av不卡久久| 99在线视频只有这里精品首页| 亚洲av成人精品一区久久| 永久网站在线| 麻豆国产97在线/欧美| 蜜桃亚洲精品一区二区三区| 久久99精品国语久久久| 久久国内精品自在自线图片| 精品国产一区二区三区久久久樱花 | 国模一区二区三区四区视频| 女人久久www免费人成看片 | 天堂中文最新版在线下载 | 国产精品av视频在线免费观看| 精品无人区乱码1区二区| 最近中文字幕2019免费版| 免费看美女性在线毛片视频| 国产亚洲精品av在线| 亚州av有码| 晚上一个人看的免费电影| 国产精品精品国产色婷婷| 国产 一区 欧美 日韩| eeuss影院久久| 非洲黑人性xxxx精品又粗又长| 精品熟女少妇av免费看| av卡一久久| 日韩在线高清观看一区二区三区| 欧美一区二区国产精品久久精品| 国产又黄又爽又无遮挡在线| 26uuu在线亚洲综合色| 热99re8久久精品国产| 精品熟女少妇av免费看| 亚洲综合精品二区| 网址你懂的国产日韩在线| 久久久精品欧美日韩精品| 成人综合一区亚洲| 亚洲精品色激情综合| 亚洲国产精品久久男人天堂| 自拍偷自拍亚洲精品老妇| 在线免费十八禁| 熟女电影av网| 一级毛片我不卡| 久久鲁丝午夜福利片| 最近视频中文字幕2019在线8| 联通29元200g的流量卡| 国产成人精品婷婷| 久久久精品大字幕| 亚洲精品,欧美精品| videossex国产| 精品久久久久久久久久久久久| 欧美成人a在线观看| 日韩欧美国产在线观看| 亚洲av中文av极速乱| 久久精品夜夜夜夜夜久久蜜豆| 最近手机中文字幕大全| 国产伦一二天堂av在线观看| 免费播放大片免费观看视频在线观看 | 亚洲伊人久久精品综合 | 亚洲天堂国产精品一区在线| h日本视频在线播放| 爱豆传媒免费全集在线观看| 亚洲欧美日韩高清专用| a级一级毛片免费在线观看| 成年女人看的毛片在线观看| АⅤ资源中文在线天堂| 午夜激情福利司机影院| 欧美高清成人免费视频www| 日韩欧美精品v在线| 乱人视频在线观看| 久久久国产成人精品二区| av天堂中文字幕网| 欧美精品国产亚洲| 桃色一区二区三区在线观看| 成人亚洲欧美一区二区av| 69av精品久久久久久| 国产探花极品一区二区| 观看美女的网站| 日韩成人伦理影院| av在线天堂中文字幕| 久久久久久久久中文| 黑人高潮一二区| 久久精品久久久久久久性| 久久久欧美国产精品| 久久精品久久久久久久性| 99久久精品一区二区三区| 日韩av不卡免费在线播放| 观看美女的网站| 人人妻人人看人人澡| 亚洲国产欧美人成| 午夜视频国产福利| 舔av片在线| 天堂网av新在线| 国产精品国产三级国产av玫瑰| 精品久久久久久久久久久久久| 久久这里只有精品中国| av黄色大香蕉| 亚州av有码| 观看免费一级毛片| 97超视频在线观看视频| 夜夜爽夜夜爽视频| 床上黄色一级片| 联通29元200g的流量卡| 久久久久久久午夜电影| 亚洲国产日韩欧美精品在线观看| 一级二级三级毛片免费看| 免费av观看视频| 一级毛片电影观看 | 欧美日本亚洲视频在线播放| 麻豆精品久久久久久蜜桃| 91久久精品电影网| 成人性生交大片免费视频hd| 国产精品一区二区在线观看99 | av在线播放精品| 免费看光身美女| 国产精品女同一区二区软件| 高清在线视频一区二区三区 | 人妻少妇偷人精品九色| 老司机影院毛片| 亚洲精品国产成人久久av| 一级毛片aaaaaa免费看小| 变态另类丝袜制服| 精品久久久久久久人妻蜜臀av| 床上黄色一级片| 国产欧美日韩精品一区二区| 噜噜噜噜噜久久久久久91| 国产午夜精品一二区理论片| 欧美一区二区亚洲| 久久久久国产网址| 丝袜美腿在线中文| 伦理电影大哥的女人| 丝袜美腿在线中文| 精品国产一区二区三区久久久樱花 | 免费黄网站久久成人精品| 3wmmmm亚洲av在线观看| 久久久久久久国产电影| 免费av毛片视频| videos熟女内射| 色吧在线观看| 成人一区二区视频在线观看| 久久久久久国产a免费观看| 赤兔流量卡办理| 激情 狠狠 欧美| 国产又色又爽无遮挡免| 国产成人91sexporn| 少妇猛男粗大的猛烈进出视频 | 美女脱内裤让男人舔精品视频| 国产精品av视频在线免费观看| 成人性生交大片免费视频hd| 日韩在线高清观看一区二区三区| 99久久精品热视频| 亚洲国产最新在线播放| 麻豆一二三区av精品| 精品一区二区免费观看| 内地一区二区视频在线| 日日撸夜夜添| 联通29元200g的流量卡| 嫩草影院精品99| 国产午夜精品论理片| 午夜亚洲福利在线播放| 看十八女毛片水多多多| 日本一二三区视频观看| 国产老妇女一区| 在线观看美女被高潮喷水网站| 在线免费观看不下载黄p国产| 能在线免费看毛片的网站| 亚洲av成人精品一区久久| 看非洲黑人一级黄片| 最近手机中文字幕大全| 九草在线视频观看| 一边亲一边摸免费视频| 九九爱精品视频在线观看| 亚洲精品亚洲一区二区| 日韩av在线大香蕉| 免费人成在线观看视频色| 亚洲国产精品合色在线| 精品久久久久久久久亚洲| 久久久色成人| 久久久久久久久久久免费av| 丰满人妻一区二区三区视频av| 一级av片app| 寂寞人妻少妇视频99o| 有码 亚洲区| 久久久久久久久久久丰满| 91午夜精品亚洲一区二区三区| 国产精品久久视频播放| 国产伦一二天堂av在线观看| 久久久精品94久久精品| 免费在线观看成人毛片| 男人狂女人下面高潮的视频| 老女人水多毛片| 精华霜和精华液先用哪个| 三级毛片av免费| 日本黄色视频三级网站网址| 自拍偷自拍亚洲精品老妇| 亚洲综合精品二区| 精品酒店卫生间| 国产精品女同一区二区软件| 国产激情偷乱视频一区二区| 爱豆传媒免费全集在线观看| 直男gayav资源| 国产一区二区在线观看日韩| 欧美丝袜亚洲另类| 国内少妇人妻偷人精品xxx网站| 日本一二三区视频观看| 三级经典国产精品| 午夜福利视频1000在线观看| 亚洲av中文av极速乱| 久热久热在线精品观看| 七月丁香在线播放| 小蜜桃在线观看免费完整版高清| 国产高清有码在线观看视频| 久久精品91蜜桃| 日韩 亚洲 欧美在线| 日本与韩国留学比较| 欧美高清性xxxxhd video| 亚洲自偷自拍三级| 久久久久免费精品人妻一区二区| 免费搜索国产男女视频| 麻豆久久精品国产亚洲av| 人妻系列 视频| 午夜a级毛片| 国产成人精品一,二区| 国产精品爽爽va在线观看网站| 成人高潮视频无遮挡免费网站| 一区二区三区乱码不卡18| 免费观看在线日韩| 亚洲综合精品二区| 97超碰精品成人国产| 亚洲av二区三区四区| 欧美zozozo另类| 国产白丝娇喘喷水9色精品| 最近中文字幕高清免费大全6| 少妇的逼水好多| 国产成人福利小说| 久久久久久久亚洲中文字幕| 丰满人妻一区二区三区视频av| 精品人妻一区二区三区麻豆| 黄片wwwwww| 久久99热6这里只有精品| 边亲边吃奶的免费视频| 99久久精品国产国产毛片| 成年版毛片免费区| 精品久久国产蜜桃| 精品国内亚洲2022精品成人| 在线免费观看的www视频| 国产精品久久电影中文字幕| 精品久久久久久久久亚洲| 婷婷色av中文字幕| 亚洲成人精品中文字幕电影| 麻豆精品久久久久久蜜桃| 午夜福利在线在线| 久久人人爽人人片av| 国产三级中文精品| 欧美xxxx性猛交bbbb| 91久久精品电影网| 国产人妻一区二区三区在| 一区二区三区高清视频在线| 老司机影院成人| 欧美精品一区二区大全| 午夜爱爱视频在线播放| 国产爱豆传媒在线观看| 人妻少妇偷人精品九色| 22中文网久久字幕| 午夜爱爱视频在线播放| 22中文网久久字幕| 狂野欧美白嫩少妇大欣赏| 久久精品人妻少妇| 非洲黑人性xxxx精品又粗又长| 男的添女的下面高潮视频| 国产精品久久久久久久久免| 国产高清视频在线观看网站| 日本熟妇午夜| 国产精品熟女久久久久浪| 国产午夜精品久久久久久一区二区三区| 99久久成人亚洲精品观看| 九九热线精品视视频播放| 网址你懂的国产日韩在线| 水蜜桃什么品种好| 婷婷六月久久综合丁香| 看黄色毛片网站| 18禁在线无遮挡免费观看视频| 又粗又爽又猛毛片免费看| 欧美丝袜亚洲另类| 少妇人妻一区二区三区视频| 精品久久久久久成人av| 高清av免费在线| 欧美激情久久久久久爽电影| 国产人妻一区二区三区在| 久久精品久久久久久噜噜老黄 | 亚洲熟妇中文字幕五十中出| 久久鲁丝午夜福利片| 久久久久网色| 精品久久久久久久久亚洲| 啦啦啦韩国在线观看视频| 99国产精品一区二区蜜桃av| 精品人妻偷拍中文字幕| 99久久人妻综合| 国产视频首页在线观看| 黄色配什么色好看| 日韩高清综合在线| 大话2 男鬼变身卡| 在现免费观看毛片| 日本黄大片高清| 成人性生交大片免费视频hd| 午夜福利在线在线| 午夜福利成人在线免费观看| 亚洲一级一片aⅴ在线观看| 日韩一区二区视频免费看| 日本免费a在线| 亚洲国产精品合色在线| 亚洲最大成人av| 尾随美女入室| 校园人妻丝袜中文字幕| www.av在线官网国产| 欧美成人一区二区免费高清观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品456在线播放app| 看非洲黑人一级黄片| 又爽又黄a免费视频| 国产亚洲91精品色在线| 久久久久久久久中文| 国产黄色小视频在线观看| 黄色日韩在线| 国内精品宾馆在线| 日韩av在线免费看完整版不卡| 国产亚洲一区二区精品| 婷婷色综合大香蕉| 国产大屁股一区二区在线视频| 99在线人妻在线中文字幕| 免费av观看视频| 岛国在线免费视频观看| 美女脱内裤让男人舔精品视频| 久久久久网色| 亚洲欧美成人综合另类久久久 | 午夜福利在线观看吧| 国产精品福利在线免费观看| 大又大粗又爽又黄少妇毛片口| 99国产精品一区二区蜜桃av| 午夜日本视频在线| 又黄又爽又刺激的免费视频.| av视频在线观看入口| 男人舔奶头视频| 日本一本二区三区精品| 欧美激情在线99| 久久久午夜欧美精品| 国产大屁股一区二区在线视频| 精品免费久久久久久久清纯| 亚洲av日韩在线播放| 99热这里只有是精品在线观看| 婷婷色综合大香蕉| 亚洲最大成人av| 欧美一级a爱片免费观看看| 精品久久国产蜜桃| 国产又色又爽无遮挡免| 欧美激情在线99| 国产熟女欧美一区二区| 成人av在线播放网站| 国产精品福利在线免费观看| 直男gayav资源| 亚洲欧洲日产国产| 国产成人freesex在线| 亚洲综合精品二区| 男人舔奶头视频| 中文字幕制服av| 久久人人爽人人爽人人片va| 欧美人与善性xxx|