姜迪迪,江為為,胥頤,郝天珧,胡衛(wèi)劍,余景鋒
1中國科學(xué)院地質(zhì)與地球物理研究所,中國科學(xué)院油氣資源研究重點(diǎn)實(shí)驗(yàn)室,北京 100029
2中國科學(xué)院大學(xué),北京 100049
中國大陸西部地區(qū),主要包括新疆、西藏和青海以及甘肅、四川西部地區(qū).地形上南高北低,南部主要是青藏高原,北部多為盆地和山脈.根據(jù)板塊構(gòu)造的觀點(diǎn),西部地區(qū)現(xiàn)在的構(gòu)造格局主要是印度板塊與歐亞板塊碰撞的結(jié)果(萬天豐,2004).板塊碰撞不僅造就了“世界屋脊”之稱的青藏高原,而且影響了青藏高原鄰近地區(qū)的構(gòu)造運(yùn)動(dòng),導(dǎo)致天山等古老造山帶的復(fù)活.目前,西部地區(qū)已經(jīng)成為研究巖石層結(jié)構(gòu)、殼幔物質(zhì)流動(dòng)以及板內(nèi)強(qiáng)震機(jī)制等前沿學(xué)科的熱點(diǎn)地區(qū)(劉光鼎,2007;徐果明等,2007;羅文行等,2008;滕吉文等,2012;張健和陳石,2013).本文利用重力、地震層析成像反演的結(jié)果、結(jié)合其他地質(zhì)與地球物理場資料,重點(diǎn)探討70°E—100°E,30°N—50°N范圍內(nèi)中國西部地區(qū)的地殼結(jié)構(gòu)與強(qiáng)震活動(dòng)的相關(guān)性,期望為開展相關(guān)研究提供參考.
中國西部具有復(fù)雜而獨(dú)特的地質(zhì)構(gòu)造特征(胥頤等,2000b;尹安,2001;萬天豐,2004)(見圖1).北部地區(qū)具有典型的三山夾兩盆構(gòu)造,塔里木盆地和準(zhǔn)噶爾盆地位于阿爾泰山、天山、昆侖山造山帶之間;天山山脈橫貫新疆全境、并向西延伸到哈薩克斯坦和吉爾吉斯境內(nèi),是一個(gè)由許多地塊和一系列造山帶組成的復(fù)式造山體系.西昆侖、阿爾金、祁連山等造山帶構(gòu)成青藏高原的西北邊緣,其中西昆侖是羌塘地體與塔里木盆地的構(gòu)造邊界、阿爾金是塔里木盆地和柴達(dá)木地體的構(gòu)造邊界、東昆侖是藏北羌塘地體和柴達(dá)木地體的構(gòu)造邊界.在南部地區(qū),喜馬拉雅、拉薩地體、羌塘地體、可可西里—甘孜—松潘地體、東昆侖柴達(dá)木地體、祁連山地體、西昆侖和喀喇昆侖山脈等構(gòu)造拼合成青藏高原的主體(尹安,2001).
圖1 研究區(qū)大地構(gòu)造概況(據(jù)尹安,2001,有修改)Fig.1 Tectonics of study area(modified from Yin,2001)
本文通過美國加利福尼亞大學(xué)圣地亞哥分校(University of California,San Diego,UCSD)斯克里普斯海洋所提供的衛(wèi)星重力數(shù)據(jù)(Sandwell and Smith,2009),得到了中國西部空間重力異常(圖2),網(wǎng)格精度為1′×1′.空間重力異常反映了實(shí)際的地球形狀和物質(zhì)分布與大地橢球體的偏差,從圖2可以看出,中國西部空間重力異常基本呈EW走向,低值區(qū)與地勢低的盆地相對應(yīng),高值區(qū)與隆起的地塊和山脈相對應(yīng),如塔里木盆地、準(zhǔn)噶爾盆地和柴達(dá)木盆地的空間重力異常值為-200~100mGal,而青藏中部和青藏南部地區(qū)的空間重力異常值為0~200mGal.較明顯的梯度帶主要分布在喜馬拉雅、昆侖、阿爾金、祁連和天山等山脈地區(qū),最高可達(dá)600mGal,空間重力異常的起伏與構(gòu)造特征吻合較好.
布格重力異常包含了殼內(nèi)各種偏離正常密度分布異常體和構(gòu)造的影響,也包括了地殼下界面起伏在橫向上相對上地幔質(zhì)量的虧損和盈余的影響(曾華霖,2005).中國西部地區(qū)布格重力異常整體特征為北高南低(圖3),變化范圍在-571~-50mGal之間.
青藏高原布格重力異常值較低,其中,在高原中部地區(qū)的唐古拉山脈可達(dá)-550mGal以下,是中國布格重力異常值最低的地區(qū).在高原的北部,柴達(dá)木盆地具有平緩的布格重力異常值,變化范圍為-425~-375mGal,在布格重力異常圖中大體可以看到盆地的形態(tài).
新疆地區(qū)的三山和兩盆是大致平行相間排列的,在形態(tài)上與布格重力異常圖有很好的對應(yīng).準(zhǔn)噶爾盆地和塔里木盆地具有相對較高的異常,重力異常在-225~-125mGal之間變化,較為平緩,大部分地區(qū)形成高值圈閉區(qū);天山和阿爾泰山等造山帶則有較低的布格異常,異常值在-275~-200mGal之間.
在重力異常圖中等值線平行排列的密集帶為重力異常梯級帶,幾條明顯的重力異常梯級帶分別圍繞著準(zhǔn)噶爾盆地、塔里木盆地、柴達(dá)木盆地和青藏高原周邊分布.主要有:(1)阿爾泰—天山帶,規(guī)模較?。唬?)昆侖—秦嶺帶,規(guī)模巨大,呈弧形展布,該帶西起昆侖山脈,向東到柴達(dá)木盆地分為兩支,北支沿阿爾金山,祁連山呈向北突出的弧形展布,又像一個(gè)頂點(diǎn)向上的三角形,它又向東南拐彎與六盤—龍門—烏蒙山梯度帶匯合,進(jìn)入緬甸;南支沿祁曼塔格山,經(jīng)阿尼瑪卿山向東南拐彎與岷山—大雪山梯度帶匯合,松潘地區(qū)受周圍梯級帶的包圍,呈倒三角形;(3)南部喜馬拉雅—橫斷山梯度帶,與中部梯度帶在中緬邊界匯合,將青藏高原與其他地區(qū)分開.這些梯級帶與中國西部地區(qū)大的構(gòu)造斷裂帶有很好的對應(yīng).
為了深入研究中國西部地區(qū)不同深度的構(gòu)造特征,利用位場小波變換法對研究區(qū)布格重力異常進(jìn)行了位場分離(圖4和圖5).1階小波變換細(xì)節(jié)中的串珠狀高頻異常反映了淺層的構(gòu)造特征,以及淺層物質(zhì)的橫向不均勻性.在天山、西昆侖、東昆侖和祁連山,以及塔里木盆地的西部巴楚隆起等,淺層地殼的重力異常表現(xiàn)為正異常,喜馬拉雅山,阿爾泰山及塔里木盆地北部凹陷則表現(xiàn)為負(fù)異常.
4階小波變換逼近主要反映區(qū)域場的趨勢,可反映深部物質(zhì)的特征.從圖5中可以看出青藏高原邊界存在梯度值較大的梯級帶,深部物質(zhì)的橫向密度在青藏高原邊界處開始發(fā)生較大的變化,產(chǎn)生了相對幅度300mGal左右的異常.
根據(jù)Parker-Oldenburg反演方法(Parker,1973;Oldenburg,1974),利用布格重力異常數(shù)據(jù)結(jié)合地震剖面(崔作舟等,1995;胥頤等,1996;陳俊湘和陳景亮,2003;王有學(xué)等,2004;王有學(xué)等,2005)控制反演計(jì)算了中國西部莫霍面深度圖(圖6),在反演計(jì)算時(shí)參考了郝天珧等(2014)的反演結(jié)果.從圖6可以看出,中國西部地區(qū)莫霍面深度特征與其地形起伏大致呈鏡像關(guān)系,且與盆地、造山帶分布有很大聯(lián)系.總體來講,盆地的莫霍面深度較小,山脈和高原的莫霍面深度較大.從莫霍面形態(tài)上看,西部地區(qū)的造山帶如天山造山帶、西昆侖—東昆侖造山帶、祁連山造山帶等,都屬于有山根的造山帶(熊小松,2010).另一方面,莫霍面深度的變化也反映了地殼運(yùn)動(dòng)的結(jié)果,青藏高原是中國大陸莫霍面深度最大的地區(qū),自50Ma年以來,印度大陸一直向北運(yùn)動(dòng)并向青藏高原的地殼擠入,造成高原地殼的縮短增厚,這一現(xiàn)象已經(jīng)被深地震測深的結(jié)果證實(shí)(Hauck et al.,1998;Chen et al.,2000).
圖2 研究區(qū)空間重力異常圖Fig.2 Free-air gravity anomalies of study area
圖3 研究區(qū)布格重力異常圖Fig.3 Bouguer gravity anomalies of study area
圖4 研究區(qū)布格重力異常小波變換1階細(xì)節(jié)Fig.4 First-degree detail wavelet of Bouguer gravity anomalies in study area
圖5 研究區(qū)布格重力異常小波變換4階逼近Fig.5 Fourth-degree approach wavelet of Bougure gravity anomalies in study area
中國大陸的莫霍面呈現(xiàn)出東淺西深的特征,而西部地區(qū)的莫霍面深度總體趨勢為南深北淺.研究區(qū)南部青藏高原地區(qū)莫霍面較深,平均深度大于60km,高原中部地區(qū)可達(dá)70km,呈向下陷的凹形.青藏高原復(fù)雜的多期演化造成了其莫霍面深度分布的特殊性,但對于青藏高原的具體形成過程至今仍存有很大的爭議.青藏高原的邊界出現(xiàn)莫霍面深度梯級帶,與地貌上的高度突變和深大斷裂構(gòu)造相對應(yīng),西南邊界喜馬拉雅山,北部昆侖—阿爾金—祁連山帶,都有較大的梯級帶,尤其以喜馬拉雅山處變化最大.喜馬拉雅碰撞帶是個(gè)非常突出的擠壓應(yīng)變帶,GPS觀測運(yùn)動(dòng)方向以向北為主,略偏向東,運(yùn)動(dòng)速率在35~42mm·a-1之間(王琪等,2001),地殼收縮量高,馮銳(1985)認(rèn)為這也是上地幔密度急劇變化較大的地方.
與南部青藏高原地區(qū)相比,新疆地區(qū)地殼厚度相對要小,塔里木盆地、準(zhǔn)噶爾盆地和吐—哈盆地等盆地和穩(wěn)定地塊都呈現(xiàn)莫霍面隆起,莫霍面深度為41~48km.塔里木盆地西部巴楚隆起莫霍面深度小于盆地中的凹陷地區(qū).準(zhǔn)噶爾盆地的莫霍面深度最淺,約41km.從莫霍面深度圖(圖6)中可以看到天山造山帶的地殼厚度從西天山到東天山的變化,東天山地殼厚度最大50km,西天山地殼厚度50~52km,而中天山莫霍面略淺于東西天山,反映了天山地區(qū)構(gòu)造運(yùn)動(dòng)的復(fù)雜性和分段性.熊小松(2010)分析莫霍面記錄并結(jié)合其他地質(zhì)地球物理資料認(rèn)為,天山西段和天山東段可能存在不同的隆升機(jī)制或處于造山作用的不同階段.東天山由于處于準(zhǔn)噶爾盆地東緣,碰撞的相向擠壓作用沒有西天山強(qiáng)烈,同時(shí)可能伴隨有上地幔的上涌,才導(dǎo)致東天山的伸展.
圖7為中國西部地區(qū)的地震層析成像結(jié)果(胥頤等,2000b),在不同深度的速度圖像中,正擾動(dòng)表示速度高于初始模型的參考值,負(fù)擾動(dòng)表示速度低于初始模型的參考值,計(jì)算時(shí)界面深度和各層的速度模型選取如表1所示.在本文中僅選取了相關(guān)的深度進(jìn)行探討.
圖6 研究區(qū)莫霍面深度圖地震剖面及編號:①阿爾金—龍門山(王有學(xué)等,2005);②阿爾泰—阿爾金(王有學(xué)等,2004);③格爾木—額濟(jì)納旗(崔作舟等,1995);④若羌—阿勒泰(陳俊湘和陳景亮,2003);⑤天山地區(qū)4條剖面(胥頤等,1996).Fig.6 Moho depth map of study area Seismic profiles:① Altun—Longmen Mountains(Wang et al.,2005);② Altay—Altun(Wang et al.,2004);③ Golmud—Ejin Qi(Cui et al.,1995);④ Ruoqiang—Altay(Chen and Chen,2003);⑤ Four profiles in Tian Shan region(Xu et al.,1996).
表1 界面深度和初始速度模型Table 1 Interface depths and initial velocity values
5+km深度的速度圖像(圖7a)反映了淺層地質(zhì)構(gòu)造的基本特征.阿爾泰、天山和準(zhǔn)噶爾西部界山均為高速區(qū);盆地和主要凹陷區(qū)如準(zhǔn)噶爾盆地、南天山的庫車坳陷、塔里木盆地西部及西南部、柴達(dá)木盆地均為低速區(qū).這一深度上的速度差異主要與巖性構(gòu)成有關(guān).造山帶主要是古生代甚至太古代的地層,通常具有較高的速度;而盆地和坳陷中相對松散、固結(jié)程度不高的中新生代沉積層速度普遍較低.
在27+km深度上(圖7b),阿爾泰山、天山和準(zhǔn)噶爾西部界山、塔里木盆地中部等地仍為高速區(qū),準(zhǔn)噶爾盆地、庫車坳陷、柴達(dá)木盆地、以及西昆侖、帕米爾等山前地帶仍為低速區(qū).27+km深度的速度分布與5+km相似,說明中地殼的結(jié)構(gòu)仍然與淺層地質(zhì)構(gòu)造有關(guān).東、西天山之間存在一條明顯的南北走向的低速帶,它與準(zhǔn)噶爾盆地南緣的低速帶相連,表明東、西天山在深部構(gòu)造上可能是不連續(xù)的.
35+km的速度圖像(圖7c)反映了下地殼的速度分布特征.阿爾泰山、準(zhǔn)噶爾盆地的北部、東部和西部一帶均為高速區(qū),準(zhǔn)噶爾盆地南部的低速區(qū)向南延伸進(jìn)入天山的中段,分隔了東、西天山高速區(qū),繼承了東、西天山構(gòu)造邊界的性質(zhì).塔里木盆地西南的巴楚隆起為高速區(qū),帕米爾東北側(cè)至塔里木盆地西南、塔里木盆地南緣至西昆侖山前兩個(gè)地區(qū)仍然為低速區(qū);甘肅的北山為高速區(qū),青藏高原的柴達(dá)木盆地為低速區(qū),上述這種速度分布從上地殼一直延續(xù)到下地殼,變化十分穩(wěn)定.
50+km的速度圖像(圖7d)反映出地殼厚度的變化以及莫霍面的大致形態(tài).通常負(fù)擾動(dòng)對應(yīng)于莫霍面坳陷區(qū),正擾動(dòng)對應(yīng)于莫霍面隆起區(qū),可以定性估計(jì)莫霍面的深淺程度.準(zhǔn)噶爾盆地中部和東部的北塔山、三塘湖一帶、哈薩克斯坦東南的阿拉湖盆地、吐魯番—哈密盆地直到甘肅的北山—敦煌、整個(gè)塔里木盆地北部和西南邊緣為高速區(qū),這些地區(qū)上地幔的頂部向上隆起,莫霍面的深度應(yīng)該小于50km.根據(jù)重力資料計(jì)算,以上大部分地區(qū)的地殼厚度為42~48km.在一條穿過準(zhǔn)噶爾盆地的人工地震測深剖面阿爾泰—阿爾金剖面上(王有學(xué)等,2004),準(zhǔn)噶爾盆地的莫霍面向上隆起,深度大約為46km;而在準(zhǔn)噶爾盆地南側(cè),莫霍面增加到50km左右,這些現(xiàn)象與本研究的結(jié)果基本一致.
圖7 研究區(qū)地震層析成像圖(a)5+km深度的速度圖像;(b)27+km深度的速度圖像;(c)35+km深度的速度圖像;(d)50+km深度的速度圖像;(e)71km深度的速度圖像.Fig.7 Seismic tomography of study area(a)Velocity image of 5+km;(b)Velocity image of 27+km;(c)Velocity image of 35+km;(d)Velocity image of 50+km;(e)Velocity image of 71km.
根據(jù)地震測深結(jié)果(邵學(xué)鐘和張家茹,1994;丁道桂等,1996;姜春發(fā),1997;張家茹等,1998;李秋生等,2002),塔里木盆地大部分地區(qū)的莫霍面深度在40~50km之間.由于上地幔頂部隆起,因此地震層析成像顯示出高速特征,但是塔里木盆地的南部表現(xiàn)出較大的差異.廣角反射地震資料(李秋生等,2002)表明塔里木地殼向西昆侖山下發(fā)生低角度俯沖,莫霍面從塔中開始由北向南逐漸加深,塔中隆起為42±2km,塔西南坳陷地殼厚度介于47~57km之間,西昆侖山下地殼最大厚度為62km.本文研究結(jié)果為解釋上述現(xiàn)象提供了有利的依據(jù).其余的莫霍面隆起區(qū)由于缺少資料很難進(jìn)行對比,但是這些隆起區(qū)大部分位于盆地和穩(wěn)定的地塊區(qū).參考前面的分析,它們的莫霍面深度應(yīng)該在45km左右.阿爾泰山、天山(包括吉爾吉斯境內(nèi)部分)、帕米爾、青藏高原的柴達(dá)木盆地均為低速區(qū),說明上述造山帶和青藏高原的平均地殼厚度超過50km,這個(gè)結(jié)果與已有的認(rèn)識(shí)基本上是一致的.沙雅—布爾津人工地震剖面(趙俊猛等,2001)穿過了天山的西段,表明天山的莫霍面深度為55km;地震轉(zhuǎn)換波測深(邵學(xué)鐘等,1996)也表明天山的莫霍面從兩側(cè)盆地邊緣的42km向山體中心逐漸加深到56km.因此,天山地區(qū)的莫霍面深度為55km左右已經(jīng)不存異議,根據(jù)地震層析成像反映的速度分布可以看出整個(gè)天山(包括吉爾吉斯境內(nèi))莫霍面坳陷的范圍.
71km深度的速度圖像(圖7e)反映了新疆地區(qū)上地幔頂部和青藏高原莫霍面附近的巖石層結(jié)構(gòu).西昆侖西段、塔里木盆地、天山和準(zhǔn)噶爾盆地的大部、阿爾金山、北山、阿爾泰東段為高速區(qū),這些地區(qū)具有相對穩(wěn)定的巖石層結(jié)構(gòu);阿爾泰山、吉爾吉斯天山、南天山至帕米爾一帶為低速區(qū),表示了上地幔頂部的巖石層異常狀態(tài).以阿爾金斷裂和西昆侖為界,整個(gè)青藏高原北部速度偏低,表明西昆侖、可可西里—甘孜—松潘地體、羌塘地體等青藏大部分地區(qū)的地殼厚度較大,可能略大于71km.
對于青藏高原的低速性質(zhì)存在兩種可能的解釋:一是反映了青藏高原年輕造山帶現(xiàn)今的構(gòu)造活動(dòng)狀況,由于深部熱物質(zhì)不斷侵入,在上地幔頂部形成低速區(qū),例如Zhao和Xie(1993)也發(fā)現(xiàn)青藏高原部分地區(qū)的上地幔頂部存在異常區(qū),宋仲和等(1991)則認(rèn)為青藏高原部分地區(qū)上地幔頂部低速層的深度大約為80~100km.二是青藏高原較大的地殼厚度造成的結(jié)果,71km深度已經(jīng)達(dá)到了塔里木盆地的上地幔頂部,但是僅相當(dāng)于青藏高原地殼的底部附近,這一差別也會(huì)引起青藏高原的速度偏低.比較有意義的是阿爾泰山上地幔頂部的低速性質(zhì),眾所周知,阿爾泰山是中國著名的貴重金屬礦產(chǎn)地,巖石層頂部的低速性質(zhì)在一定程度表明了上地幔物質(zhì)的熱活動(dòng)狀態(tài),它可能成為巖漿熱液上涌的源區(qū).在這個(gè)深度上,周圍的天山、準(zhǔn)噶爾西部界山等造山帶都沒有出現(xiàn)低速區(qū),因此這也許是阿爾泰山富含金屬礦床的深部構(gòu)造背景之一.
綜上所述,可以得出以下結(jié)論:中國西部的造山帶一般是莫霍面坳陷區(qū),盆地則是莫霍面隆起區(qū);盡管塔里木盆地北部的莫霍面向上隆起,但是塔里木盆地南部至西昆侖山前莫霍面的深度卻較大;沿著祁連山的南北兩側(cè),敦煌地塊和柴達(dá)木盆地的地殼速度截然不同,北部莫霍面隆起,南部莫霍面坳陷.
比較地震層析成像結(jié)果與重力反演計(jì)算結(jié)果可以看出:準(zhǔn)噶爾盆地中部區(qū)域?yàn)楦咚賲^(qū),地殼厚度小于50km,重力反演42~45km,這與人工地震測深結(jié)果45km相符;整個(gè)塔里木盆地北部和西南邊緣為高速區(qū),這些地區(qū)上地幔的頂部向上隆起,莫霍面的深度小于50km,重力反演結(jié)果表明地殼厚度46~49km,地震測深資料表明塔里木盆地大部分地區(qū)的莫霍面深度在40~50km之間;柴達(dá)木盆地主要為低速區(qū),地殼厚度應(yīng)大于50km,重力反演結(jié)果在59~64之間;在這些地區(qū)重力反演結(jié)果與地震層析成像結(jié)果十分吻合.而在一些地區(qū)也存在兩種方法結(jié)果不同的現(xiàn)象,天山地區(qū)為低速區(qū),地震層析成像結(jié)果表明地殼厚度應(yīng)大于50km,重力反演結(jié)果為49~52km,而人工地震探測剖面和地震轉(zhuǎn)換測深結(jié)果表明天山地區(qū)莫霍面深度為55km左右;阿爾泰山為低速區(qū),表明地殼厚度大于50km,而重力反演結(jié)果是46~50km;在71km深度的速度圖像中,青藏高原中部地區(qū)為低速區(qū),表明地殼厚度大于71km,而重力反演在60~70km左右,推測兩種反演結(jié)果差異較大主要由地勢較高地區(qū)測量平面的差異引起的.
中國西部是地震活動(dòng)十分強(qiáng)烈的地區(qū),地震密集分布在天山、阿爾泰山、帕米爾和西昆侖等構(gòu)造帶,這些地區(qū)不僅在歷史上發(fā)生過多次強(qiáng)烈的地震,至今地震活動(dòng)仍然十分頻繁.根據(jù)中國地震臺(tái)網(wǎng)(CSN)地震目錄,從1970年1月1日到2014年2月28日,研究區(qū)共發(fā)生5級以上地震1059個(gè),地震的震中和震源深度分布如圖8所示.
對比重力及地震層析成像的反演結(jié)果,從研究區(qū)地震震中和震源深度分布(圖8)中可以看出,中國西部地區(qū)5級以上地震的震源深度主要分布在40km以內(nèi),少量地震震源深度大于60km,多位于帕米爾高原與喜馬拉雅構(gòu)造帶等印度板塊與歐亞板塊的碰撞俯沖帶.在應(yīng)力場作用下,造山帶上地殼剛性較強(qiáng),以脆性破裂為主,相比之下,下地殼在壓力、溫度的影響下易于發(fā)生塑性變形,不利于強(qiáng)震的發(fā)生.近年來的研究表明,發(fā)生強(qiáng)震的構(gòu)造環(huán)境與介質(zhì)的橫向不均勻性密切相關(guān),中國西部的強(qiáng)震多發(fā)生在剛性地塊和褶皺帶的交接邊界,即垂直差異運(yùn)動(dòng)強(qiáng)烈的壓性逆斷裂帶和壓性剪切斷裂帶上,如天山、昆侖山、以及龍門山和鮮水河地震帶等(胥頤等,2000a;胥頤等,2006;田有等,2007).
在北部新疆地區(qū),地殼厚度一般在50km以內(nèi),天山及帕米爾地區(qū)地震幾乎分布在整個(gè)地殼厚度,而在西昆侖地區(qū),地震的分布也達(dá)到了莫霍面的深度,這與Keith Priestley等(2008)觀點(diǎn)相似.Keith Priestley等在分析印度、喜馬拉雅和藏南地區(qū)巖石圈結(jié)構(gòu)和地震的關(guān)系時(shí)認(rèn)為,在有些地區(qū),整個(gè)地殼厚度,甚至到莫霍面都是孕震的,大部分位于前寒武紀(jì)地盾的內(nèi)部或邊緣.他還通過分析地震分布與巖石圈的溫度結(jié)構(gòu)和速度結(jié)構(gòu)的關(guān)系,認(rèn)為地震一般發(fā)生在三種情況,溫度在350℃左右的“濕”的上地殼物質(zhì),具有較高溫度的干的麻粒巖相的下地殼,或者溫度低于600℃的地幔,具有非常重要的意義.本文重點(diǎn)從中國西部的地殼速度結(jié)構(gòu)及運(yùn)動(dòng)速率等其他地質(zhì)地球物理特征,分析地殼結(jié)構(gòu)與強(qiáng)震活動(dòng)的相關(guān)性.
地殼的速度結(jié)構(gòu)是反映地殼介質(zhì)屬性的重要基礎(chǔ),可以提供地震震源和介質(zhì)結(jié)構(gòu)的重要信息,與地震活動(dòng)的空間分布相結(jié)合,為認(rèn)識(shí)地震成因及發(fā)震構(gòu)造機(jī)理提供重要的依據(jù).對于天山地殼發(fā)震環(huán)境,本文通過穿過天山及帕米爾地區(qū)的3條地震層析成像剖面進(jìn)行分析,從地殼速度結(jié)構(gòu)方面分析該地區(qū)的發(fā)震機(jī)制.
對比地震層析成像結(jié)果和天山地區(qū)三條地震層析成像剖面(圖9)可以看出,大部分地震都發(fā)生在低速帶向相鄰高速塊體過渡的區(qū)域.天山北部的強(qiáng)震主要分布在天山和準(zhǔn)噶爾盆地之間的邊緣地帶,從剖面AA′圖和BB′圖中可以看到北天山和準(zhǔn)噶爾盆地地殼速度變化劇烈,準(zhǔn)噶爾盆地南部具有雙層地殼結(jié)構(gòu),上地殼和下地殼的速度差別很大:上地殼為高速層,中地殼沿水平方向出現(xiàn)的低速層是上地殼和下地殼之間的韌性滑脫帶,上地殼受南北方向的擠壓作用沿滑脫層向兩側(cè)盆地推覆逆沖,導(dǎo)致結(jié)合帶內(nèi)部地震頻發(fā).
從剖面圖BB′中可以看出,庫車坳陷基底深度可以達(dá)到15km左右,在20至30km之間出現(xiàn)低速層,該低速層向北延伸到天山的中地殼.庫車坳陷上地殼的低速反映了沉積層的特點(diǎn),中地殼和下地殼至莫霍面的低速估計(jì)與深斷裂的性質(zhì)有關(guān).天山南緣的深斷裂多為巖石圈斷裂,地表出露華力西期花崗巖,推測邊界附近上地幔熱物質(zhì)可能沿著斷裂的通道上侵,在地殼中部形成低速層,地殼熱狀態(tài)的改變影響了地殼的流變性質(zhì),進(jìn)而產(chǎn)生地殼局部的變形與地震的發(fā)生.
從剖面CC′中可以看到天山—帕米爾結(jié)合帶的地震活動(dòng)與地殼結(jié)構(gòu)的組成也有著內(nèi)在聯(lián)系,地震大都發(fā)生在塔里木西部邊緣的高速塊體周圍.沿東經(jīng)76°附近橫跨帕米爾—天山的最新GPS測量結(jié)果表明其現(xiàn)今地殼縮短速率為20~24mm·a-1,幾乎占印度板塊向北推擠速率的一半(李濤,2012).帕米爾和塔里木兩個(gè)剛性塊體的相對運(yùn)動(dòng)在構(gòu)造邊界附近產(chǎn)生較大的應(yīng)力差,構(gòu)造變形主要通過兩者之間的邊界斷裂喀什—葉城轉(zhuǎn)換帶的走滑運(yùn)動(dòng)進(jìn)行調(diào)節(jié),伽師、喀什一帶的低速帶正是起到了韌性剪切帶的作用,因此,塔里木和帕米爾兩個(gè)剛性塊體在喀什韌性剪切帶東、西兩側(cè)的運(yùn)動(dòng)差異是伽師地震的深部成因之一.
根據(jù)鄭洪偉(2006)對青藏高原沿88°E剖面的層析成像結(jié)果(剖面DD′,圖10),MCT(主中央逆沖斷層)以南的地殼表現(xiàn)為高速帶,該高速帶從地殼一直向地幔延伸.MCT和IYS(印度河—雅魯藏布江縫合帶)之間存在一個(gè)低速體,即喜馬拉雅低速帶,可能是高溫和局部熔融的反映,而班公—怒江縫合帶以北的羌塘地體表現(xiàn)為低速,被解釋為是印度板塊的向北俯沖下插引起的地幔物質(zhì)的上涌,在地殼內(nèi)部形成低速層.從圖10中也可以看到青藏高原地震的分布與地殼速度結(jié)構(gòu)有一定的聯(lián)系.南部高速的印度板塊的和低速的喜馬拉雅碰撞帶之間,以及喜馬拉雅低速帶與北部縫合帶之間強(qiáng)震頻發(fā).推測韌性剪切產(chǎn)生的局部熔融和板塊俯沖產(chǎn)生的構(gòu)造應(yīng)力,使喜馬拉雅的北部和南部邊緣都孕育了大量地震,尤其是該低速帶與南部高速帶相接的地方更容易發(fā)生強(qiáng)震.另外,地幔上涌使得羌塘地體與邊緣地帶的地殼結(jié)構(gòu)產(chǎn)生橫向不均勻,導(dǎo)致羌塘地體邊緣,即班公—怒江縫合帶和金沙江縫合帶附近成為強(qiáng)震易發(fā)地帶.
結(jié)合以往的研究及4條地殼結(jié)構(gòu)剖面,可以看出強(qiáng)震主要分布在地殼速度結(jié)構(gòu)差異較大的地區(qū).地質(zhì)演化,如斷裂形成,地殼組成和熱狀態(tài)、板塊運(yùn)動(dòng)和變形等都會(huì)導(dǎo)致地殼速度的不均勻分布.不管是熱流物質(zhì)上涌,或是地殼斷裂導(dǎo)致的地殼介質(zhì)產(chǎn)生強(qiáng)烈非均勻性,在構(gòu)造應(yīng)力作用下,這些地殼介質(zhì)非均勻地區(qū)易發(fā)生強(qiáng)震,這是中國西部造山帶和盆—山邊界附近頻發(fā)強(qiáng)震的構(gòu)造原因之一.
圖8 研究區(qū)地震震中和震源深度分布(圖中AA′、BB′、CC′為圖9截取3條剖面的位置,DD′為圖10中剖面)Fig.8 Epicenter and focal depth distribution of study area(AA′、BB′、and CC′are locations of three profiles in Fig.9,DD′in Fig.10)
圖9 地殼速度剖面圖(單位:km·s-1)(a)AA′北天山地殼速度剖面;(b)BB′庫車—烏蘇地殼速度剖面;(c)CC′新疆烏恰—伽師地震區(qū)地殼速度剖面.Fig.9 Crustalvelocity profiles(unit:km·s-1)(a)AA′North Tian Shan crustal velocity profile;(b)BB′Kuqa—Usu crustal velocity profile;(c)CC′Wuqia—Jiashi earthquake zone crustal velocity profile.
圖10 青藏高原沿88°E的剖面結(jié)果(據(jù)鄭洪偉,2006,有修改)Fig.10 Profile of major earthquakes and seismic tomography along 88°E in Tibetan Plateau(modified from Zheng,2006)
通過分析中國大陸西部地區(qū)的地殼結(jié)構(gòu)以及與地震活動(dòng)的相關(guān)性,取得了以下的認(rèn)識(shí):
(1)中國西部地區(qū)的莫霍面深度總體趨勢為南深北淺:南部青藏高原地區(qū)莫霍面深度較大,平均深度大于60km,高原中部地區(qū)可達(dá)70km以上,呈向下陷的凹形.與南部青藏高原地區(qū)相比,北部的新疆地區(qū)地殼厚度相對要小,特別是盆地地區(qū)更淺.總體上莫霍面起伏與重力異常特征的形態(tài)相對應(yīng).
(2)地殼結(jié)構(gòu)不均勻特征明顯,主要造山帶上地殼為高速區(qū),盆地和凹陷區(qū)多為低速區(qū).中地殼的結(jié)構(gòu)仍然與地表的地質(zhì)構(gòu)造有關(guān),其中東天山與西天山之間存在南北走向的低速帶,并且與準(zhǔn)噶爾盆地南緣的低速帶相連,表明東、西天山在深部構(gòu)造上是不連續(xù)的.青藏高原的柴達(dá)木盆地為低速區(qū),從上地殼一直延續(xù)到下地殼,變化十分穩(wěn)定.中國西部地殼厚度的變化很大,莫霍面形態(tài)復(fù)雜,造山帶一般是莫霍面坳陷區(qū),盆地則是莫霍面隆起區(qū).大部分地區(qū)莫霍面重力反演結(jié)果與地震層析成像結(jié)果十分吻合,但在部分地區(qū)存在不同的現(xiàn)象,如阿爾泰和青藏高原中部地區(qū).
(3)中國西部的地震活動(dòng)與地殼結(jié)構(gòu)的橫向不均勻密切相關(guān),強(qiáng)震多發(fā)生在地殼速度結(jié)構(gòu)差異較大的地區(qū),如天山、帕米爾、西昆侖等造山帶與毗鄰盆地的結(jié)合部位.地質(zhì)演化,如斷裂形成,地殼組成和熱狀態(tài)、板塊運(yùn)動(dòng)和變形等都會(huì)導(dǎo)致地殼速度的不均勻分布.不管是熱流物質(zhì)上涌,或是地殼斷裂導(dǎo)致的地殼介質(zhì)產(chǎn)生強(qiáng)烈非均勻性,在構(gòu)造應(yīng)力作用下,這些地殼介質(zhì)非均勻地區(qū)易發(fā)生強(qiáng)震,這是中國西部造山帶和盆-山邊界附近頻發(fā)強(qiáng)震的構(gòu)造原因之一.
致謝 本文是在劉光鼎院士中國大陸構(gòu)造格架思想指導(dǎo)下完成的研究,在此表示衷心地感謝.同時(shí)感謝審稿專家提供的寶貴修改意見.
Chen J X,Chen J L.2003.The crustal structure and velocity feature from Ruoqiang-Altai seismic profile in Xinjiang region.XinjiangPetroleumGeology(in Chinese),24(6):498-501.
Chen Z,Bruchfiel B C,Liu Y,et al.2000.Global positioning system measurements from eastern Tibet and their implications for India/Eurasia intercontinental deformation.J.Geophys.,Res.:SolidEarth(1978—2012),105(B7):16215-16227.
Cui Z Z,Li Q S,Wu C D,et al.1995.The crustal and deep structures in GOLMUD-EJIN QI GGT.ChineseJ.Geophys.(in Chinese),38(S2):15-28.
Ding D G,Wang G X,Liu W X,et al.1996.West Kunlunorogenic Belt and Basin (in Chinese).Beijing:Geological Publishing House.
Feng R.1985.The crustal thickness and distribution of upper mantle density of China.ActaSeismologicaSinica(in Chinese),31(7):529-536.
Hao T Y,Hu W J,Xing J,et al.2014.The Moho depth map(1∶5000000)in the land and seas of China and adjacent areas and its geological implications.ChineseJ.Geophys.(in Chinese),57(12):3869-3883,doi:10.6038/cjg20141202.
Hauck M L,Nelson K D,Brown L D,et al.1998.Crustal structure of the Himalayan orogen at~90°east longitude from Project INDEPTH deep reflection profiles.Tectonics,17(4):481-500.
Jiang C F.1997.Opening-closing tectonics of Tarim Platform.XinjiangGeology(in Chinese),15(3):193-202.
Li Q S,Gao R,Lu D Y,et al.2002.Tarim crustal subduction in the West Kunlun Mountains—the evidence of wide-angle seismic survey.ActaGeoscienticaSinica(in Chinese),23(5):442.
Li T.2012.Active thrusting and folding along the Pamir Frontal Thrust System(in Chinese)[Ph.D.thesis].Beijing:Institute of Geology,China Earthquake Administration.
Liu G D.2007.Geodynamical evolution and tectonic framework of China.EarthScienceFrontiers(in Chinese),14(3):39-46.
Luo W X,Li D W,Wang X F.Focal depth and mechanism of intraplate earthquakes in the Qinghai-Tibet Plateau.Earth Science(JournalofChinaUniversityofGeosciences) (in Chinese),33(5):618-626.
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.Geophys.J.Int.,31(4):447-455.
Priestley K,Jackson J,McKenzie D.2008.Lithospheric structure and deep earthquakes beneath India,the Himalaya and southern Tibet.Geophys.J.Int.,172(1):345-362.
Sandwell D T,Smith W H F.2009.Global marine gravity from retrackedGeosat and ERS-1altimetry:Ridge segmentation versus spreading rate.J.Geophys.Res.:SolidEarth(1978—2012),114(B1):B01411,doi:10.1029/2008JB006008.
Shao X Z,Zhang J R.1994.Preliminary results of study of deep structures in Tarim Basin by the method of converted waves of earthquakes.ChineseJ.Geophys.(in Chinese),37(6):836.
Shao X Z,Zhang J R,F(xiàn)an H J,et al.1996.The crust structures of Tianshanorogenic belt:A deep sounding work by converted waves of earthquakes alongürümqi-Korla profile.ActaGeophysica Sinica(in Chinese),39(3):336-346.
Song Z H,An C Q,Chen G Y,et al.1991.Study on 3Dvelocity structure and anisotropy beneath the west China from the Love wave dispersion.ActaGeophysicaSinica(in Chinese),34(6):694-707.
Teng J W,Yuan X M,Zhang Y Q,et al.2012.The stratificational velocity structure of crust and covering strata of upper mantle and the orbit of deep interaquifer substance locus of movement for Tibetan Plateau.ActaPetrologicaSinica(in Chinese),28(12):4077-4100
Tian Y,Zhao D P,Sun R M,et al.2007.The 1992Landers earthquake:Effect of crustal heterogeneity on earthquake generation.ChineseJ.Geophys.(in Chinese),50(5):1488-1496.
Wan T F.2004.Tectonics of China(in Chinese).Beijing:Geological Publishing House.
Wang Q,Zhang P Z,Niu Z J,et al.2001.Chinese contemporary crustal movement and tectonic deformation.ScienceinChina(in Chinese),31(7):529-536.
Wang Y X,Han G H,Jiang M,et al.2004.Crustal structure along the geosciences transect from Altay to AltunTagh.ChineseJ.Geophys.(in Chinese),47(2):240-249.
Wang Y X,Mooney W D,Han G H,et al.2005.Crustal P-wave velocity structure from AltynTagh to Longmen mountains along the Taiwan-Altay geoscience transect.ChineseJ.Geophys.(in Chinese),48(1):98-106.
Xiong X S.2010.Moho depth and variation of the continent in China and its geodynamic implications(in Chinese)[Ph.D.thesis].Beijing:China Academy of Geological Sciences.
Xu G M,Yao H J,Zhu L B,et al.2007.Shear wave velocity structure of the crust and upper mantle in western China and its adjacent area.ChineseJ.Geophys.(in Chinese),50(1):193-208.
Xu Y,F(xiàn)eng X Y,Shstsilov V I.1996.Deep velocity structure in the Tianshan region.SeismologyandGeology(in Chinese),18(4):375-381.
Xu Y,Liu F T,Liu J H,et al.2000a.Crustal structure and tectonic environment of strong earthquakes in the Tianshan Earthquake Belt.ChineseJ.Geophys.(in Chinese),43(2):184-193.
Xu Y,Liu F T,Liu J H,et al.2000b.Seismic tomography on northwest orogenic belt and adjacent basins of Chinese Mianland.ScienceinChina(SeriesD)(in Chinese),30(2):113-122.
Xu Y,Liu J H,Liu F T,et al.2006.Crustal velocity structure and seismic activity in the Tianshan-Pamir conjunctive zone.Chinese J.Geophys.(in Chinese),49(2):469-476.
Yin A.2001.Geologic evolution of the Himalayan-Tibetan orogen—Phanerozoic growth of Aisacontinetal.ActaGeoscientiaSinica(in Chinese),22(3):193-230.
Zhang J R,Shao X Z,F(xiàn)an H J.1998.Deep sounding survey by converted waves of earthquakes in central part of the Tarim Basin and its interpretation.SeismologyandGeology(in Chinese),20(1):34-42.
Zhang J,Chen S.2013.Gravity field analysis for crust and mantle structure of the Tibetan Plateau.Earthquake(in Chinese),33(4):11-18.
Zhao J M,Liu G D,Lu Z X,et al.2001.The crust-mantle transition zone beneath the Tianshanorogenic belt and the JunggarBasin and their geodynamic implications.Sciencein China(SeriesD)(in Chinese),31(4):272-282.
Zhao L S,Xie J k.1993.Lateral variations in compressional velocities beneath the Tibetan Plateau from Pntraveltime tomography.Geophys.J.Int.,115(3):1070-1084.
Zeng H L.2006.Gravity Field and Gravity Exploration (in Chinese).Beijing:Geological Publishing House.
Zheng H W.2006.3-D velocity structure of the crust and upper mantle in Tibet and its geodynamic effect(in Chinese)[Ph.D.thesis].Beijing:China Academy of Geological Sciences.
附中文參考文獻(xiàn)
陳俊湘,陳景亮.2003.若羌—阿勒泰區(qū)域剖面的速度特征和地殼結(jié)構(gòu).新疆石油地質(zhì),24(6):498-501.
崔作舟,李秋生,吳朝東等.1995.格爾木—額濟(jì)納旗地學(xué)斷面的地殼結(jié)構(gòu)與深部構(gòu)造.地球物理學(xué)報(bào),38(S2):15-28.
丁道桂,王桂軒,劉偉新等.1996.西昆侖造山帶與盆地.北京:地質(zhì)出版社.
馮銳.1985.中國地殼厚度及上地幔密度分布(三維重力反演結(jié)果).地震學(xué)報(bào),7(2):143-157.
郝天珧,胡衛(wèi)劍,邢健等.2014.中國海陸1∶500萬莫霍面深度圖及其所反映的地質(zhì)內(nèi)涵.地球物理學(xué)報(bào),57(12):3869-3883,doi:10.6038/cjg20141202.
姜春發(fā).1997.塔里木地臺(tái)開合構(gòu)造簡述.新疆地質(zhì),15(3):193-202.
李秋生,高銳,盧徳源等.2002.塔里木地殼俯沖在西昆侖山之下——寬角地震測量的證據(jù).地球?qū)W報(bào),23(5):442.
李濤.2012.帕米爾前緣逆沖推覆系活動(dòng)斷層和活動(dòng)褶皺作用[博士論文].北京:中國地震局地質(zhì)研究所.
劉光鼎.2007.中國大陸構(gòu)造格架的動(dòng)力學(xué)演化.地學(xué)前緣,14(3):39-46.
羅文行,李德威,汪校鋒.2008.青藏高原板內(nèi)地震震源深度分布規(guī)律及其成因.地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報(bào)),33(5):618-626.
邵學(xué)鐘,張家茹.1994.塔里木盆地深部構(gòu)造的地震轉(zhuǎn)換波探測結(jié)果.地球物理學(xué)報(bào),37(6):836.
邵學(xué)鐘,張家茹,范會(huì)吉等.1996.天山造山帶地殼結(jié)構(gòu)與構(gòu)造──烏魯木齊—庫爾勒地震轉(zhuǎn)換波測深剖面.地球物理學(xué)報(bào),39(3):336-346.
宋仲和,安昌強(qiáng),陳國英等.1991.中國西部三維速度結(jié)構(gòu)及其各向異性.地球物理學(xué)報(bào),34(6):694-707.
滕吉文,阮小敏,張永謙等.2012.青藏高原地殼與上地幔成層速度結(jié)構(gòu)與深部層間物質(zhì)的運(yùn)移軌跡.巖石學(xué)報(bào),28(12):4077-4100.
田有,趙大鵬,孫若昧等.2007.1992年美國加州蘭德斯地震——地殼結(jié)構(gòu)不均勻性對地震發(fā)生的影響.地球物理學(xué)報(bào),50(5):1488-1496.
萬天豐.2004.中國大地構(gòu)造學(xué)綱要.北京:地質(zhì)出版社.
王琪,張培震,牛之俊等.2001.中國大陸現(xiàn)今地殼運(yùn)動(dòng)和構(gòu)造變形.中國科學(xué),31(7):529-536.
王有學(xué),韓果花,姜枚等.2004.阿爾泰—阿爾金地學(xué)斷面地殼結(jié)構(gòu).地球物理學(xué)報(bào),47(2):240-249.
王有學(xué),Mooney W D,韓果花等.2005.臺(tái)灣—阿爾泰地學(xué)斷面阿爾金—龍門山剖面的地殼縱波速度結(jié)構(gòu).地球物理學(xué)報(bào),48(1):98-106.
熊小松.2010.中國大陸莫霍面深度與變化特征及其地球動(dòng)力學(xué)意義[博士論文].北京:中國地質(zhì)科學(xué)院.
徐果明,姚華建,朱良保等.2007.中國西部及其鄰域地殼上地幔橫波速度結(jié)構(gòu).地球物理學(xué)報(bào),50(1):193-208.
胥頤,馮先岳,Shstsilov V I.1996.天山地區(qū)的深部速度結(jié)構(gòu).地震地質(zhì),18(4):375-381.
胥頤,劉福田,劉建華等.2000a.天山地震帶的地殼結(jié)構(gòu)與強(qiáng)震構(gòu)造環(huán)境.地球物理學(xué)報(bào),43(2):184-193.
胥頤,劉福田,劉建華等.2000b.中國大陸西北造山帶及其毗鄰盆地的地震層析成像.中國科學(xué)(D輯:地球科學(xué)),30(2):113-122.
胥頤,劉建華,劉福田等.2006.天山—帕米爾結(jié)合帶的地殼速度結(jié)構(gòu)及地震活動(dòng)研究.地球物理學(xué)報(bào),49(2):469-476.
尹安.2001.喜馬拉雅—青藏高原造山帶地質(zhì)演化——顯生宙亞洲大陸生長.地球?qū)W報(bào),22(3):193-230.
張家茹,邵學(xué)鐘,范會(huì)吉.1998.塔里木盆地中部地震轉(zhuǎn)換波測深及其解釋.地震地質(zhì),20(1):34-42.
張健,陳石.2013.青藏高原重力場與殼幔結(jié)構(gòu)分析.地震,33(4):11-18.
趙俊猛,劉國棟,盧造勛等.2001.天山造山帶與準(zhǔn)噶爾盆地殼幔過渡帶及其動(dòng)力學(xué)含義.中國科學(xué)(D輯),31(4):272-282.
曾華霖.2005.重力場與重力勘探.北京:地質(zhì)出版社.
鄭洪偉.2006.青藏高原地殼上地幔三維速度結(jié)構(gòu)及其地球動(dòng)力學(xué)意義[博士論文].北京:中國地質(zhì)科學(xué)院.