立即咨詢
您當前的位置:職稱驛站 > 論文 > 理工論文 > 環境科學論文職稱驛站 期刊論文發表 權威認證機構

基于近震P波走時的南加州地區速度結構層析成像

職稱驛站所屬分類:環境科學論文發布時間:2020-10-27 08:39:41瀏覽:1

采用2015—2018年南加州地區密集分布的寬頻帶地震臺站記錄的505個近震事件,利用AIC準則拾取了35 600個初至P波到時,應用FMTOMO軟件進行走時層析成像,獲得了該地區0~40 km深度的三維P波地殼速度結構

   摘要:采用2015—2018年南加州地區密集分布的寬頻帶地震臺站記錄的505個近震事件,利用AIC準則拾取了35 600個初至P波到時,應用FMTOMO軟件進行走時層析成像,獲得了該地區0~40 km深度的三維P波地殼速度結構,其分辨率達0.5°×0.5°×5 km。不同深度的速度結構表明:南加州地區P波速度不僅隨深度的變化而變化,且該地區地殼具有明顯的橫向不均勻性。反演結果顯示:中地殼以上的速度結構與地殼巖層和斷層系統密切相關,下地殼以下的速度結構存在大范圍低速區,呈現橫向不均勻性。

  關鍵詞:南加州;P波走時;三維地殼速度結構模型;地震層析成像

  中圖分類號:P315.2文獻標識碼:A文章編號:1000-0666(2020)04-0674-06

0a1e5b3f013e4db17507498510eb3aa6.jpg

  《華北地震科學》(季刊)1983年創刊,本刊是以報導、交流地震科研成果為主的綜合性學術季刊,主要刊登與地震科學有關的實驗、觀測、考察、理論和科研成果、問題討論、防震減災等方面的論文,以及國內外地震科學研究狀況的介紹和評論性文章。

  0引言

  美國南加州地區位于太平洋板塊和北美板塊的交匯地帶,是世界上地震活動最活躍的地區之一。該地區自中生代以來經歷了強烈的地質構造活動,經過多次強烈的擠壓和伸展作用,形成了特殊的地質構造格局,出現了板塊俯沖、島弧形成、大陸增生等現象(田有,2008)。歷史上南加州地區多次發生中強地震,如舊金山1906年8.6級大地震和1989年6.9級地震、1994年洛杉磯6.6級地震等。2019年南加州發生7.1級強震,刷新了近20年來最強烈地震的記錄,再次引起了社會各界廣泛的關注。據統計,南加州地區大部分中強地震發生于地殼層,對人類的生命和財產安全危害極大。因此研究南加州地區的地殼速度結構對于認識該地區的地質構造、輔助地震的分析預報等具有十分重要的作用。

  地震層析成像技術是近年來獲取地球內部非均勻結構的一種重要手段,可使用的地震資料豐富,理論方法多樣。如利用地震體波、面波資料或者地震接收臺站之間的噪聲記錄都可以重建研究區的速度結構模型(Yang,Forsyth,2006;Barak,2015)。相較于地震的其他信息,地震體波的走時信息具有更穩健、更可靠、更易提取且數據精度高的優點,因此地震體波走時層析成像技術的應用最為廣泛,在南加州地區的地下結構研究中取得了大量成果(Lin et al,2007;Tape et al,2010;Barak et al,2015)。如Hauksson(2000)利用體波走時資料獲得了南加州地區0~22 km深度的地殼速度結構,但該模型對于南加州地區地殼的成像深度不夠;Tian等(2007)聯合近震和遠震數據對南加州地區進行體波走時層析成像,獲得了該地區地殼上地幔的三維速度結構,但垂向分辨率較低。

  本文利用南加州地區近震P波的初至到時,采用FMTOMO軟件進行地震走時層析成像,得到南加州地區分辨率為0.5°×0.5°×5 km的三維地殼P波速度結構。

  1數據資料與研究方法

  1.1數據資料

  本文研究區位于美國西部的南加州地區(121°~114°W,32°~38°N),使用2015—2018年研究區內239個寬頻帶地震臺站記錄的505個近震事件的波形數據資料進行研究,臺站和地震分布如圖1所示。事件挑選原則為:①首先挑選震中距為0~5°、震級大于3級、震源深度為0~40 km的近震事件;②對數據進行重采樣、去毛刺、去均值、去線性、滅尖、帶通濾波等預處理;③使用全球一維速度模型AK135,計算各地震事件理論到時,并將波形截取為理論到時的前10 s至后20 s;④利用AIC法(Maeda,1985)實現實際到時的自動拾取,并進行波形的質量篩選,刪除信噪比低的波形數據;⑤保證每個地震事件至少具備10個以上的臺站到時,剔除不符合此要求的地震事件。由于近震震中距較小,后續震相不便于自動拾取,本文僅使用P波的初至到時進行層析成像,最終挑選了35 600條初至P波和Pn波到時。

  1.2研究方法

  地震層析成像過程主要包括模型參數化、正演、反演以及解的評價(羅炬等,2011)。本文采用Rawlinson等(2006)發展的FMTOMO軟件進行近震P波走時層析成像,該軟件具有僅需地震走時信息、追蹤多種震相、可聯合多種數據集進行反演等優點(趙啟光等,2011),是近年來研究地球內部結構較流行的一種方法。

  地震研究43卷第4期周茜茜等:基于近震P波走時的南加州地區速度結構層析成像該軟件采用網格結點法進行模型參數化,采用的網格結點包括速度網格結點、界面網格結點和傳播網格結點。其中傳播網格結點代表正演計算得到的走時場,速度網格結點和界面網格結點構成反演初始模型。本文采用的初始速度模型為Hadley和Kanamori(1977)提出的一維速度模型,如圖2所示。對于界面網格結點,在深度方向上以29 km的平均莫霍面深度為界將研究區分為兩層,并在深度方向上以每2 km設置速度網格結點,在經度和緯度方向以15 km×15 km設置速度網格結點。本文反演過程中保持震源初始位置不變,僅對P波速度進行反演。

  FMTOMO軟件正演方法采用多步快速行進法 (Kool et al,2006),該方法根據給定的初始速度模型,利用迎風有限差分算法求解程函方程,從而得到各傳播網格結點的走時值。反演方法采用子空間反演法(Rawlinson et al,2006),該反演方法將非線性問題線性化進行迭代求解,通過求解當前模型的走時數據與觀測走時數據構建的目標函數的極小值獲得速度模型的擾動量,進行速度場的解算。

  2解的評價

  2.1檢測板測試

  在分析結果之前,為了驗證層析成像結果的可靠性,對層析成像結果進行評價,可了解層析成像結果的可靠性、分辨率及誤差等信息(雷棟,胡祥云,2006)。本文采用檢測板分辨率測試:對初始模型添加±0.4 km/s的速度擾動值,保持地震接收臺站和震源的數量以及位置分布不變,計算出一組走時代替實際觀測走時,進行層析成像,反演得到的速度結構可直觀地分析層析成像結果的可靠性。

  本文采用水平格網結點間距約為0.5°,垂向格網結點間距約為5 km的檢測板進行測試,測試結果如圖3所示:1 km水平切片的檢測板結果顯示,邊緣地區由于臺站和地震稀少,檢測板還原效果較差,而臺站和地震密集分布的區域,還原效果較好;11 km水平切片的檢測板還原效果較好,這是由于地震事件震源深度大部分都位于20 km以上,該深度上射線密集;22 km水平切片的圖3檢測板測試結果

  Fig.3Checkerboard test results檢測板結果顯示,由于所選地震事件大多發生于上地殼,且本文僅采用了初至的直達P波和Pn波,所以該深度射線密度不足,此處檢測板還原效果較差,但依然可以看出檢測板還原出正負相間的速度擾動;35 km水平切片的檢測板還原效果較好,因為此處Pn波射線能較好地交叉穿過且射線密集。對于垂向剖面的檢測板,在0~15 km以及30~40 km深度還原效果較好,而15~30 km深度還原效果稍差。

  2.2反演殘差分析

  本文經過5次迭代計算,速度結構模型的解逐漸收斂。觀測走時數據與反演走時數據走時差方差從2.015 s2下降到1.361 s2,改善了近32%;卡方值從222.259下降到119.801,改善了46%。

  3層析成像結果與討論

  3.1層析成像結果

  由反演結果(圖4a)可以看出,在1 km深度,研究區內一些低地形區域如圣華金山谷、文圖拉盆地和莫哈韋沙漠南部呈現低速分布,這是由于圣華金山谷堆積著新近紀的沉積物,文圖拉盆地下方堆積著較厚的上新世的沉積物(Shaw et al,2015),莫哈韋沙漠南部在1~3 km的深度存在沉積盆地(Tape et al,2010)。一些高山山脈呈現明顯高速,如內華達山脈、橫斷山脈、半島山脈、海岸山脈、圣蓋博山脈、圣貝納迪諾山脈。中生代的花崗巖在半島山脈、橫斷山脈和內華達山脈南部的基巖中占據主要地位(Phil,2013),圣貝納迪諾山脈巖層分布以侏羅紀花崗巖侵入性火成巖為主(宋國政,2011),從沉積物到火成巖、花崗巖的過渡,定義了南加州大部分地區橫向速度的不均勻性。圣安德烈斯斷層兩端速度分布存在明顯差異,其WS側速度較高,EN側速度偏低。由圖1可以看出,大多數地震分布在圣安德烈斯斷層附近,尤其是其WS側地震活動較多,該現象與圣安德烈斯斷層附近的速度分布情況相符。由此可見,該地區淺部的P波速度分布與地表的地質構造以及斷層系統的分布具有較好的一致性。

  由圖4b可見,在11 km深度,大部分地區呈現高速分布,但異常幅度較1 km深度小。圣安德烈斯斷層兩端速度分布依舊存在差異,半島山脈持續高速,大陸邊緣帶、圣華金山谷、思爾頓凹槽、圣瑪麗亞盆地呈現高速分布,而洛杉磯、文圖拉等盆地的低速特征消失,呈現高速分布。在南加州地區的沉積盆地如洛杉磯盆地、文圖拉盆地、圣瑪麗亞盆地等在地殼中形成了重要的速度結構。這些盆地的基底巖通常被沉積物覆蓋,其中洛杉磯盆地的基底巖由火成巖構成(Wright,1987),文圖拉盆地和圣瑪麗亞盆地的基底巖由早期的第三紀和中生代的變質巖構成(Yerkes et al,1987)。這種由沉積物到基底巖的分布,反映了南加州盆地地區垂向上速度的不連續性。Fuis等(1984)的研究成果表明,在該深度上,半島山脈和思爾頓凹槽下方分布著鎂鐵質巖層,地震波較同深度其他巖層傳播速度快。在莫哈韋沙漠東南部的蘭德斯地區,有一小范圍區域呈現高速異常,并被低速塊體截斷,這樣的速度分布與該地區地震頻發情況相符。前人的研究也表明,地震多發生于高速和低速塊體的過渡部位或偏向高速塊體(田有等,2007)。該地區中地殼的速度結構不僅與地殼巖層構成有關,與以往的地質構造活動也相關。

  在22 km深度,反演結果(圖4c)顯示了下地殼速度的橫向不均勻性,半島山脈、思爾頓凹槽、圣瑪麗亞盆地、文圖拉盆地及蘭德斯地區持續高速。圣瑪麗亞盆地和文圖拉盆地的高速分布受新近紀橫斷山脈順時針旋轉的影響,為俯沖捕獲的法拉隆洋殼所導致(Tape et al,2010)。內華達山脈呈現低速分布,多學科綜合研究揭示內華達山脈巖石圈地幔以及下地殼在8~3.5 Ma范圍內發生拆沉,而拆沉將導致該深度上P波低速(曾令森等,2006;Cin-Ty等,2001)。圣蓋博山脈和海岸山脈呈現低速分布,該地區下地殼存在早期第三紀的沉積巖和變質沉積巖(Shaw et al,2015)。

  在35 km深度,圣華金山谷、內華達南部、半島山脈和莫哈韋沙漠呈現明顯低速,橫斷山脈、海岸山脈、思爾頓凹槽、文圖拉盆地、圣瑪麗亞盆地和洛杉磯盆地呈現明顯高速(圖4d)。該結果與Zhu和Kanamori(2000)利用遠震接收函數得到的南加州莫霍面深度結果相對應,其研究結果顯示在圣華金山谷、內華達山脈、半島山脈和莫哈韋沙漠地區地殼深度較深,在該深度上這些區域顯示的是下地殼的速度特性,因而速度較低;而思爾頓凹槽、橫斷山脈、文圖拉盆地、圣瑪麗亞盆地和洛杉磯盆地地殼厚度較薄。

  3.2對比分析

  為了進一步分析本文層析成像結果的可靠性,將本文結果與Tian等(2007)聯合近震與遠震的P波走時資料獲得的相似區域的層析成像結果進行對比分析。由于所使用的數據集、模型參數化以及從屬格網的劃分不同,再加上反演問題的多解性,本文層析成像結果雖與其存在一些差異,但仍然具有較好的一致性,具體為:地表層P波速度分布與地表的地質構造、地形起伏具有較好的一致性,如半島山脈、圣蓋博山脈以及海岸山脈等高山山脈的高速異常,圣華金山谷和文圖拉盆地的低速異常。在中地殼大部分區域都顯示高速異常,如半島山脈和圣華金山谷等。在下地殼,如半島山脈、文圖拉盆地等呈現高速異常,內華達山脈呈現低速異常。在上地幔,大部分區域呈現低速異常,如圣華金山谷、內華達山脈、半島山脈等,而在研究區的西部,如海岸山脈、沉積盆地等表現為高速異常。與Tian等(2007)結果相比,本文還揭示了蘭德斯地區較清晰的速度異常。

  4結論

  本文利用2015—2018年南加州臺網239個地震臺站記錄的共505個近震波形數據,采用AIC法提取了35 600條P波初至到時信息,應用FMTOMO軟件進行地震走時層析成像,確定了南加州地區分辨率達0.5°×0.5°×5 km的0~40 km深度的三維P波速度結構,得到以下結論:

  (1)本文僅需引入近震波形數據的絕對到時信息參與層析成像,所需數據相較于其它方法具有易獲取、精度高、成像方法簡單快速的優點。得到的層析成像結果與前人研究成果總體上具有較好的一致性,且垂向分辨率得到了較大改善,有助于準確預測強地面震動,更有助于進一步認識該地區的地殼構造和地質構造活動。因此該方法可用來獲取區域高分辨率的地殼速度結構。

  (2)層析成像結果反映了南加州地區P波速度隨深度發生變化,且地殼存在明顯橫向不均勻性。淺部的速度結構較好地恢復了該地區地表地形的大致位置,如洛杉磯盆地、文圖拉盆地等沉積盆地和半島山脈、海岸山脈等高山山脈的位置,與地表的地質構造和巖性密切相關,深部的速度結構存在大面積的低速區,呈現橫向不均勻性。

  由于數據的限制,某些深度上層析成像結果的分辨率稍差,有待進一步改善,在后續的研究中需要搜集更多的數據,如通過提取后續震相的信息來增加地震射線的密度,從而獲得更高分辨率的三維速度結構。

  感謝審稿專家和編輯部的幫助,感謝Rawlinson提供的FMTOMO層析成像軟件。

  參考文獻:

  雷棟,胡祥云.2006.地震層析成像方法綜述[J].地震研究,29(4):104-112.

  羅炬,李志海,王海濤.2011.地震層析成像研究方法綜述[J].內陸地震,25(4):27-36.

  宋國政.2011.美國加州圣本那迪諾郡奧德山地區成礦地質條件[J].科技資訊,(29):122-122.

  田有,趙大鵬,孫若昧,等.2007.1992年美國加州蘭德斯地震——地殼結構不均勻性對地震發生的影響[J].地球物理學報,50(5):1488-1496.

  田有.2008.中國華北和美國加州地區地震層析成像研究[D].北京:中國科學院地質與地球物理研究所.

  曾令森,劉靜,Jason B S.2006.大型花崗巖基形成和演化的深部動力學過程:滴水構造、鉀質火山作用與地表地質過程[J].地質通報,25(11):21-37.

  趙啟光,樓海,王椿鏞.2011.青藏高原東緣深部速度結構遠震層析成像[J].CT理論與應用研究,20(3):311-320.

  Cin-Ty L,任向文,吳福元.2001.加利福尼亞州內華達山脈巖石圈地幔中生代移離的鋨同位素證據[J].世界地質,20(4):327-331.

  Barak S,Lawrence J F,Klemperer S L.2015.San andreas fault dip,peninsular ranges mafic lower crust and partial melt in the Salton Trough,Southern California,from ambient-noise tomography[J].Geochemistry,Geophysics,Geosystems,16(11):3946-3972.

  Kool M D,Rawlinson N,Sambridge M.2006.A practical grid-based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous media[J].Geophysical Journal of the Royal Astronomical Society,167(1):253-270.

  Fuis G S,Mooney W D,Healy J H,et al.1984.A seismic refraction survey of the Imperial Valley Region,California[J].Journal of Geophysical Research:Solid Earth,89(B2):1165-1189.

  Hadley D M,Kanamori H.1977.Seismic structure of the transverse ranges,California[J].Geological Society of America Bulletin,88(10):1469-1478.

  Hauksson E.2000.Crustal structure and seismicity distribution adjacent to the Pacific and North America Plate Boundary in Southern California[J].Journal of Geophysical Research,105(B6):13875-13903.

  Lin G,Shearer P M,Hauksson E,et al.2007.A three-dimensional crustal seismic velocity model for Southern California from a composite event method[J].Journal of Geophysical Research:Solid Earth,112:B11306.

  Maeda N.1985.A method for reading and checking phase times in autoprocessing system of seismic wave data[J].Zisin,38(3):365-379.

  Phil S.2013.Geologic history of california[EB/OL].(2013-09-17) [2020-01-01].http://www.geologycafe.com/geologic_history/index.html.

  Rawlinson N,Reading A M,Kennett B L N.2006.Lithospheric structure of Tasmania from a novel form of teleseismic tomography[J].Journal of Geophysical Resrarch,111:B02301.

  Shaw J H,Plesch A,Tape C,et al.2015.Unified structural representation of the southern california crust and upper mantle[J].Earth and Planetary Science Letters,415:1-15.

  Tape C,Liu Q,Maggi A,et al.2010.Seismic tomography of the Southern California crust based on spectral-element and adjoint methods[J].Geophysical Journal of the Royal Astronomical Society,180(1):433-462.

  Tian Y,Zhao D,Teng J.2007.Deep structure of southern California[J].Physics of the Earth and Planetary Interiors,165(1-2):93-113.

  Wright T L.1987.Structural geology and tectonic evolution of the Los Angeles Basin[M].United States:AAPG Bull.

  Yang Y,Forsyth D W.2006.Regional tomographic inversion of the amplitude and phase of rayleigh waves with 2-D sensitivity kernels[J].Geophysical Journal International,166(3):1148-1160.

  Yerkes R F,Sarna-Wojcicki A M,Lajoie K R.1987.Geology and quaternary deformation of the ventura area[J].U S Geological Survey Professional Paper,1339:169-178.

  Zhu L,Kanamori H.2000.Moho Depth variation in Southern California from teleseismic receiver functions[J].Journal of Geophysical Research,105(B2):2969-2980.

  Velocity Structure Tomography in Southern California Area Based

  on Local Earthquakes P-wave Travel-time

  ZHOU Xixi,CHEN Qiang,ZHANG Yijun,HUANG Xiaomei

  (Faculty of Geosciences and Environmental Engineering,Southwest Jiaotong University,Chengdu 611756,Sichuan,China)

  Abstract

  In this study,basing on 505 local-seismic events recorded by densely distributed broadband seismic stations in Southern California between Jan.,2015 and Dec.,2018,we pick up 35 600 first-arrival travel-time of P-waves by adopted the AIC method.Then using the FMTOMO software,we obtained a three-dimensional P-wave crust velocity structure with a depth range of 0~40 km in this area through travel-time tomography,and the resolution of the velocity structure can reach 0.5°×0.5°×5 km.The velocity structure at different depth revealed that the P-wave velocity in Southern California changed with the depth,and the crust in this region had obvious lateral heterogeneity.The inversion results show that the velocity structure above the mid-crust is closely related to the rock formations and fault systems,and that below the lower crust has a large area of low-speed,showing lateral heterogeneity.

  Keywords:Southern California;P-wave travel-times;three-dimensional crustal velocity structure;seismic tomography

《基于近震P波走時的南加州地區速度結構層析成像》

本文由職稱驛站首發,一個權威專業的職稱論文發表網

文章名稱:基于近震P波走時的南加州地區速度結構層析成像

文章地址:http://www.guantaorf.com/lunwen/ligong/hjkx/43325.html

'); })(); 可以免费观看的av毛片