地震波层析成像.ppt

《地震波层析成像.ppt》由会员分享，可在线阅读，更多相关《地震波层析成像.ppt（56页珍藏版）》请在三一办公上搜索。

1、第六章 地震层析成像,Seismic wave is currently the only effective tool that can penetrate the entire earth Structural information of the Earth,From IRIS,Seismic waves,Jeffreys-Bullen1-D Earth Model,1939:Jeffreys&Bullen First travel-time tables：Jeffreys-Bullen Seismological Tables 1D Earth model,地球内部结构,PREM一维全

2、球速度模型,1.PREM参考地球模型：Preliminary Reference Earth Model(Dziewonski&Anderson,1981),2.IASP91速度模型(Kennett and Engdahl,1991),3.AK135 速度模型Kennett et al.(1995)；Montagner and Kennett(1996).,另外 还有：MC35，STW105，TNA/SNA模型等,莫霍诺维齐（前南斯拉夫）莫霍面（壳幔边界）地壳与地幔分界面的发现者古登堡（美国）核幔边界 地幔与地核分界面的发现者莱曼（丹麦，女）（内外核边界）外核与内核分界面的发现者,Topogr

3、aphy,地球一维结构是远远不够的！,mantle convection,Plate tectonics,Shearer,2009,Travel time table from ak135 model,Travel time picks,3-D variations of Earths Structure from Seismic Tomography,Seismic waves in the Earth,Inverse problem,Researchers at MIT and Harvard,led by Keiti Aki and Adam Dziewonski in late 197

4、0s and 1980s,pioneered the technique of seismic tomography.,PREM模型给出了地球的一维结构，而地球内部三维结构需要更精细的刻画。地震层析成像方法是给出地球内部三维结构的最重要的方法。某种意义上说，地震是照亮地球内部的明灯。地震层析成像方法可以给出：全球地球结构的横向不均匀性；典型地球动力学过程的三维结构：俯冲带、地幔柱、大洋中脊等；小尺度的构造（断层等）；地震分布特征。,全球地震层析成像,地幔柱,区域地震层析成像,区域地震层析成像,大洋中脊,区域地震层析成像,俯冲带,局部地震层析成像,Where it all began:Radon

5、 transform:(Johan Radon,1917):integral of function over a straight line segment.,where p is the radon transform of f(x,y),and is a Dirac Delta Function(an infinite spike at 0 with an integral area of 1)p is also called sinogram,and it is a sine wave when f(x,y)is a point value.,Radon transform,Tomo

6、Greek for“tomos”(body),graphy-study or subject,地震层析成像的基础Radon 变换,Back projection of the function is a way to solve f()from p()(“Inversion”):,Shepp-Logan Phantom(human cerebral),Input,Radon Projected,Recovered(output),Different“generations”of X-Ray Computed Tomography(angled beams are used to increas

7、e resolution).Moral:good coverage&cross-crossing,Cunningham&Jurdy,2000,A few of the early medical tomo setups,Parallel beam,Fan beam,Multi-receiver,Moves in big steps,Broader fan beam,Coupled,moving source receivers,fast moving,Broader fan beam,Moving source,fixed receivers,fast moving(1976),18,Brai

8、n Scanning Cool Fact:According to an earlier report,the best valentines gift to your love ones is a freshly taken brainogram.The spots of red shows your love,not your words!,Present Generation of models:Dense receiver sets,all rotating,great coverage and crossing rays.,19,20 世纪60 年代初期，美国科学家Cormack 从

9、数学和实验结果证实了根据X射线的投影可以唯一地确定人体内部结构,从而奠定了医学诊断上图像重建的理论基础,即X 射线CT(X Ray Computer Tomography)。,地震波层析成像的发展历史,在数学方法上出现了本质上与奥地利数学家1917年提出的Rndon逆变换方法相同的褶积投影方法，Chapman首先从理论上证明了地震学中的-P 变换即是Radon 变换(Chapman,1981)。,地震波层析成像首先由Aki等提出，并给出了小尺度（Aki and Lee,1976）和区域尺度（Aki et al.,1977）远震体波层析成像（Teleseismic body-wave tomog

10、raphy）。Dziewonski等在1977年给出了全球尺度的体波层析成像成果（Dziewonski et al.,1977）。,体波层析成像：主流，对于远震距离，其可沿着多种穿过地球核%幔等路径传播面波层析成像：仅能对上地幔及其以上部分提供约束，适合于台站稀疏和缺乏地震的区域，,按成像的范围来分主要有三种：全球层析成像，区域层析成像和局部层析成像。,全球层析成像,全球间断面，全球波速场分布等；,区域层析成像,研究地幔柱，洋中脊，板块俯冲带等,局部层析成像,石油、天然气、煤矿、矿产等勘探，地下水流向等,1.体波层析成像2.面波层析成像3.噪音层析成像4.衰减层析成像5.有限频层析成像,按成像

11、的方法可以分为：,体波速度：P，SV，SH,面波速度：Rayleigh，Love,类似于面波成像,品质因子，衰减系数,非射线理论，考虑地震波频率,1.体波走时层析成像,2.有限频带走时层析成像,Fat-ray tomography(Husen and Kissling,2001)Banana-doughnut tomography(Dahlen et al.,2000),这两种方法主要是利用了地震波传播时，影响地震波传播的不是简单的射线，而是形状与香蕉类似（即具有一定直径的弯曲传播路径）的区域，而震源和记录台站是该香蕉的两个端点。这两种方法优点是在依据地震波传播理论的基础上，考虑到了地震波频率

12、与模型介质尺度间的影响，同时也用到震相走时信息。,3.地震波衰减层析成像方法(Attenuation tomography,Rietbrock,2001),该方法是利用了地震波在传播过程中，地震波的固有频率会随着周期变化成幂指数衰减的关系（即传播介质的品质因子）及其与波速的关系，对周围介质进行成像，其优点是能够给出近地表介质的更为精细结构特征。,体波具有相对小的振幅和尖脉冲特征（Shearer，1999）,finite frequency kernels for travel time perturbations,对地震波的传统认识认为，震发生后所记录的有效地震事件波形为有效信号。除此外，其余

13、部分被称为噪音。这样，形成噪音的声源广泛地分布在地球表面的每个角落，如海浪、风暴和公路上的车流等，统被称为随机分布的波场(Yang and Ritzwoller,2008)在一个特定环境中，只要是记录这些随机分布的波场足够长，这些噪音也是有规律的，就可以很好地认识区域噪音规律，对区域进行结构成像，跟面波成像方法相似。,4.噪音层析成像（Ambient seismic noise tomography,Shapiro et al.,2005;Sabra et al.,2005;Pollitz and Flecher,2005）,Ambient noise is enriched at short

14、 periods.Better constraints on crustal and uppermost mantle structure than information from earthquakes.Particularly useful in aseismic areas;e.g.,continental interiors.For temporary deployments-do not have to wait for earthquakes to occur.Measurements are repeatable:rigorous uncertainty estimates.,

15、R.Weaver,Science,2005,Processing Steps:Remove instrument response,de-mean,de-trend,bandpass filter,time-domain normalization,spectral whitening Cross-correlation:1 day at a time.Stack over many days.Waveform selection(SNR)for tomography,Mike Ritzwoller et al.,2008,体波走时层析成像,1.块体模型方法（Aki and Lee,1976）

16、,该方法是将整个反演模型用多个均匀的六面体表示，而每个六面体中心的速度来表示该六面体的整体速度。,块体=Block,Roecher（1982）提出了模型空间由多个尺度大小不同的六面体来描述，即可变块体方法。该方法在我国首先由刘福田等（1989）运用。,后来者在块体模型中引入了不连续界面（Wintlinger et al.,1998,2004;胥颐等，2000）,2.格点模型方法,Thurber（1983）提出了六面体内的速度变化用其八个顶点的速度变化来表示的技术。,之后有人提出利用四个顶点组成的四面体代替六面体（Lin and Roecker,1997）,在六面体格点模型的基础上，通过引入Sn

17、ell定律和射线伪弯曲法，Zhao(1992,1994)利用相关区域的先验信息资料（如莫霍面、康拉德面以及板块边界等）构建较为接近真实的模型结构。,Travel time(or slowness)inversions:,慢度的扰动，作为反演的基本变量,30,Menke,1989,Least-Squares Solutions,Suppose we have a simple set of linear equations,A X=d,We can define a simple scalar quantity E,Mean square error(or total error),Error

18、function,We want to minimize the total error,to do so,find first derivative of function E and set to 0.,So,do,we should have,This is known as the system of normal equations.,So this involves the inversion of the term ATA,This matrix is often called the inner-product matrix,or Toeplitz matrix.The sol

19、ution is called the least-squares solution,while X=A-1 d is not a least squares solution.,最小二乘法适合于观测数据个数多于未知数个数。,Left multiply by,Pre-conditioning for ill-conditioned inverse problem(damping,smoothing,regularization.Purpose:Stabilize,enhance smoothness/simplicity),Lets use the same definition,Define

20、:An objective function J where,where m is the damping or regularization parameter.,I is identity matrix,Minimize the above by,(Damped Least Squares solution),阻尼最小二乘法，适合于观测数据个数少于未知数个数。,Solve for the inverse problem,(1)Standard Least Squares Solution,GTG may be singular or ill-conditioned singular val

21、ue decomposition(SVD),(2)Damped Least Squares Solution,minimize,Ray-based traveltime tomography,minimize,L:Laplacian operator,Solution:m=(GTG+2LTL)-1GTd,3.Smooth model,Combined norm and Laplacian regularization,misfit function,Crust correction,Data misfit,Model roughness,Model norm,在Zhao 等(1992,1994

22、)方法基础上,双差走时定位(double difference location;Waldhauser and Ellsw orth,2000;Wolfe,2002)方法被引入到地震波层析成像反演中,即地震波双差走时层析成像方法(Zhang and Thurber,2003,2005),基本概念震相对：如果两个时间到同一个站台有相同震相的走时数据，那么对于这两个事件来说这两个震相就是一个震相对。2.事件对：两个事件有一定数目的震相对，并且震中之间的距离在一定的范围内。,Double difference method,汪锐等，2013,参考模型（Reference model）：计算理论走时。

23、地震波走时自动拾取：获得观测走时.,Crust correction:using 3-D Crust 2.0 as the reference crust model,Crust 2.0,Input model,1-D crust reference model,3-D crust reference model,检测版试验是为了检查模型的分辨率。,Checkerboard resolution tests,Checkerboard and synthetic resolution tests,地震波走时自动拾取问题,相邻道互相关方法(Gelchinsky，1983)、能量比较方法(Coppe

24、ns，1985)、改变褶积算子宽度的方法(Ramananantoandro，1987)，分形理论的Divider 方法和Hurst 方法(Boschetti，1996)。中国科学家提出的基于hausdorff 分维算法的地震波走时全自动拾取方法(Chang，1999)。,走时是地震层析成像的最基础数据，这个数据量一般都是很大的，需要研究走时的自动拾取技术。,A dense path coverage minimizes the amount of a priori information needed,射线覆盖密度,射线覆盖越密，反演结果越好！,P波隐区104-143度,一次较大的地震激发的地

25、震波，可以在很长的距离上观测到，结果发现在1040突然消失了，到了1430重又出现!也就是说，1040 到1430“销声匿迹”。这是因为地球内存在一个速度很低的“核”，“影区”就是地核的影子。,S波隐区104-180度,from:,P波隐区：104-143度S波隐区：104-180度两侧皆是,Mantle tomography,E.g.,Bijwaard,Spakman,Engdahl,1999.,Subducting Farallonslab,Montelli et al.,Science,303,2004,Plume clusters or superplumes,Montelli et al.,Science,303,2004,Montelli et al.,Science,303,2004,Dapeng Zhao,2012,