讲地震波层析成像ppt课件
复杂构造地震成像PPT课件
53
基于绕射理论的主能量偏移成像 （反射+绕射）
• 绕射偏移原理
• 走时：T=TSD+TDG
54
走时计算
• 走时计算：以射线参数p、介质速度v为参数
dz
dt
t(p)
,
v2 (1v2 p2) v (1v2 p2)
x(p) pdz
vp dz v2 p dt
20
Tm
21
• 常用的叠前时间偏移走时计算是双平方根公式， 其实就是分段单平方根走时计算之和：
t tm 2 vxr22 tm 2 vxs22
22
叠前深度偏移
• 叠前时间偏移的成像质量较叠后偏移大幅 提高，但没办法解决速度横向变化情况下 其下地层变形问题和不成像的问题。而且， 在横向变速快的情况下，在波动方程中， Tm不能代替Z。
• 绕射波成像原理 原理等同于基于绕射理论的叠前偏移，所 不同的是输入数据中不能包含反射波。
• 绕射波成像的关键 反射波与绕射波的分离
• 实现方式 主能量部分主要包含反射波，非主能量不 能主要是绕射波。反射波能量一致性强， 变化慢，绕射波能量各道变化快。
84
过夏72井inline剖面
85
过夏72井绕射波成像剖面-石炭-二叠系
11
时间偏移与深度偏移
• 波动方程：记录的波由空间域到时间域
2P2P2P12P 2x 2y 2z v2 2t
Tm
dz v
22P x22P yv222P Tmv12 22P t
12
速度横无变化
x
x
垂向单位变换
y y
z
Tm
时间与深度变换其实就是坐标轴变化，都是指向地 下（深度）方向，单位不同
地震层析成像技术 ppt课件
二、地震层析成像方法面临的主要问题
2.1 地震波走时自动拾取问题 在地震层析成像的研究中 ,可获得的观测数据是地震 记录 .从地震记录中可以获得地震波的走时、振幅和 频率 ,其中最关键的是地震波走时 .随着数字地震技 术的发展 ,观测数据的数量迅速增加 ,准确地进行地 震波走时的拾取越来越成为一项重要且繁重的工作 . 为此 ,走时的自动拾取成为人们研究与关注的对象 .
二、地震层析成像方法面临的主要问题
2.2 三维波动方程有限差分算法模拟地震波场的问题 开展非弹性介质和完全弹性介质有限差分法三维
一、地震层析成像研究发展概况
地震层析成像的研究在70年代首先以井间速度结构 调查为研究对象（Bois et al.1972）。1979年， Dines和Lytle首先对地震层析成像坐了大量数值模 拟，并公布了利用弯曲的地震射线进行地下地震波 速度成像的结果，并首先将层析成像 （ Computerized Geophysical Tomography ） 这 一 名 词 用 于 论 文 的 标 题 。 1984 年 ， 美 国 的 Anderson利用天然地震数据着手全 化、密度结构、地幔物质流动有了新的认识。
二、地震层析成像方法面临的主要问题
2.2 三维波动方程有限差分算法模拟地震波场的问题 不论是天然地震还是人工地震 (即使是二维观测方 式 )的观测数据都是在三维空间介质中形成 .由于地 下地质结构的千变万化 ,理论数据的正演计算只有在 三维空间中实现才更具有实际意义 .而目前大多采用 二维计算 ,使得理论数据与观测数据之间的误差不仅 由地质模型形成而且还由计算方法的数学模型形成 . 三维波动方程的有限差分解是获取地震波三维波场 的有效方法 .
一、地震层析成像研究发展概况
20世纪60年代初期，美国科学家Cormack从数学和 实验结果证实了根据X射线的投影可以唯一地确定人 体内部结构，从而奠定了医学诊断上图像重建的理 论 基 础 ， 即 X 射 线 CT(X Ray Computer Tomography). 60年代中期和70年代中期,随着数 学图像重建方法在射电天文学和电子显微学方面的 应用和发展，在数学方法上出现了本质上与奥地利 数学家1917年提出的Rndon逆变换方法相同的褶积 投影方法，Chapman，1981）。此后，地学界借 助医学CT思想，利用地震波的传播对地壳乃至上地 幔结构开始进行半定量研究。从此，低着层析成像 成为地球物理学研究的一个新领域。
地震层析成像技术
地震层析成像技术《地震层析成像技术：探索地球内部的神奇“X光”》地震层析成像技术，听起来就像是一种很神秘又高大上的东西呢。
其实呀，咱们可以把地球想象成一个超级大的苹果，这个苹果内部有着各种各样不同的结构，就像苹果有果核、果肉、果皮一样。
而地震层析成像技术呢，就像是一种特殊的手段，能让我们看到这个大苹果内部到底是啥样的。
咱们知道地震会产生地震波，这些地震波就像一个个小信使，在地球内部跑来跑去。
地震层析成像技术就是抓住这些小信使，从它们身上挖掘地球内部的秘密。
不同的物质对地震波的传播速度、方向等会有不同的影响。
比如说，坚硬的岩石可能就会让地震波跑得比较快，而那些软软的、像岩浆一样的物质可能就会让地震波的速度变慢。
这就好比在不同的路况下开车，在平坦的高速公路上汽车就能开得快，在坑坑洼洼的土路上汽车就只能慢慢开。
科学家们在地球的不同地方放置很多的监测仪器，就像在苹果周围放了好多小耳朵。
当地震发生的时候，这些小耳朵就能听到地震波的动静。
然后通过收集到的这些地震波的信息，利用复杂的数学计算和算法，就能构建出地球内部的图像。
这图像就像是地球内部的一张大地图，能告诉我们哪里是高山的根基在地下的延伸，哪里可能藏着神秘的地下洞穴，又或者哪里有大片的岩浆在涌动。
这技术在很多方面都超级有用呢。
比如说寻找石油。
石油就像地球这个大仓库里的宝藏，藏得很深。
以前找石油就像在黑夜里摸瞎，现在有了地震层析成像技术，就像是给找石油的人戴上了夜视镜。
通过这个技术看到地下的结构，就能大概知道哪里可能有石油的踪迹，这样就能更有针对性地去钻探，节省很多的人力、物力和时间。
还有在研究板块运动方面。
地球的板块就像漂浮在大锅里的几块大饼，一直在动来动去。
地震层析成像技术能让我们看到板块在地下是怎么相互作用的。
是紧紧地挤压在一起，还是有一部分已经插到另一个板块下面去了。
这就像看大饼之间是怎么碰撞、叠加的一样。
通过这个技术，我们就能更好地预测地震的发生。
地震偏移成像基本原理ppt课件(共114张PPT)
三大处理技术:
反褶积、叠加、和偏移成像
反褶积和叠加引自其它相关学科 偏移成像基于古典技术
偏移成像: 1.具有地震勘探本身的特征。
过渡到地震波动力学特 征
§1.1 偏移成像的基本原理
一．偏移成像的概念
偏移
反偏移
反射地震方法:
1.激发弹性波，2.记录反射波, 3.研究地质岩层结构和物 性特征。是一种反散射问题。
（1.1.24）
此时反向外推遇到倏逝波，正向外推发生耗损波。分别表示为：
考虑到我们的边界条件是线性的，可以求出反射系数：
40a〕式可改写为：
38〕公式组可以看出，后两种展开是等价的。
9〕式得出F-k域的向下外推公式
20〕式完全相同，因此〔1.
z方向上差分网格向下外推时不重叠，速度变化可稍大些。
根据这个公式可以进行地震记录的向下半空间延拓，求出地下任何一点的波场，实现地震波偏移的目的。
这个方程可用来模拟下行波的地震记录。
（2〕下行波反向外推公式 下行波的反向外推是指沿负z方向的外 推。其外推式为：
（1.1.20）
上式可用来从下行波场进行反向求源的计算工作。
下面分析波场本身的条件对外推结果的影响
（1.1.21）
当
时， 为正或负的实数，这时所有外推公式中存在虚指
数。说明在外推过程中波场发生相位变化。一般都能得出正确的结
2). Kirchhoff积分法波动方程法偏移：70年代中期，French和 Schneider等在绕射偏移法的基础上使用了波动方程解的Kirchhoff积 分公式，发展为地震偏移的波动方程积分法。使绕射偏移建立在可 靠的波的基本原理上。因而改善了偏移剖面，取得了良好的效果。
地震波层析成像和电磁波层析成像
地震波层析成像和电磁波层析成像地震波层析成像和电磁波层析成像1.地震波CT地震层析成像的主要目标是确定地球内部的精细结构和局部不均匀性。
这不仅可以促进地球科学的发展，而且还可以解决许多地质勘探和矿产资源开发中的难题。
第一个原因是岩石地震波与岩性性质有比较稳定的相关性，易于对地球内部成像，反之，对找水活确定流体性质时，电磁波层析成像较好。
第二个原因是对于主要频段的电磁波，其衰减比地震波大。
对于地址勘探、采矿工程、勘察工程等来说目标提一般为几米到几百米，对应波长为几十米，频率为数十赫兹。
这种的地震波在不松散的岩石中传播为几公里后耍贱一般不超过120dB，接收起来不费力。
反而相应波长的电磁波在岩石中传播几十米后就可能衰减100dB，难以穿透几百米的岩层。
第三个原因是电磁波速度太快，反映波速的到时参数难以测量。
地震波波速为每秒几千米，振幅、到时都易于测量，而且在地震记录上可以区分不同的震相，从而得到丰富地质信息。
1.井间地震波数据的采集方法一般地层观测排列均匀布置在风化层一下，以使提高成像分辨率。
一般采集方法及对应的观测方式有：1.共激发点道集数据采集方法单点激发，多点接收的观测方式采集地震数据。
这种方法比较适用于在震源连续性能较差且接收为多道检波系统的情况下使用。
这种方法有采集快，效率高的特点。
但要求至少有一口井的井深超过目的层且满足目的层覆盖要求。
2.共接收点道集数据采集方法这种方法以移动式多点源激发，单点接收的观测方式采集地震数据。
适合在震源连续激发性能较好且接收器为单级检波器系统情况下使用。
但施工效率不高，也有井深要求。
3.YO-YO道集数据采集这种方法采用激发点和接收点反向移动的观测方式采集地震数据。
要求震源系统具有良好的连续激发性能，获得道集多用于反射波成像。
适合井深不符合透射层析成像要求的目的层成像问题。
4.井间地震连续测井方法这种方法采用激发点和接收点等间距同向移动的观测方式采集地震数据。
《地震波》PPT课件教案模板
解：4/Sin30=5/Sin(b1) b1=38.7 临界角a1(P波折射角为90） 4/Sin(a1)=5/Sin90 a1=53.1 临界角a2(P波反射角为90） 4/Sin(a2)=6/Sin90 a2=41.8
作业：求P波反射角c1和Sv波折射角b
page19
Institute of Disaster Prevention
地震波
Add the author and the accompanying title
2.2 地震波
➢ 波的特点 ➢ 地震波的特点 ➢ 波形转换 ➢ 斯奈尔定律 ➢ 面波的特点 ➢ 频散现象 ➢ 地震波序列：震相
page2
Institute of Disaster Prevention
2021/8/9
page27
Institute of Disaster Prevention
2021/8/9
P震相和S震相
在震中距为 105°的范围以内，P 震相是地震图上的初 至震相。其后是S震 相，其振幅、周期都 比P震相大，质点运 动垂直于传播方向。
page28
Institute of Disaster Prevention
2021/8/9
界面波
波动入射至界面，还会发生更复杂的
转换现象。例如，当折射波或反射波的波 速大于入射波波速时，折射角或反射角将 大于入射角，90°的折射角或反射角对应 的入射角称为临界入射角。当入射角大于 临界入射角时，将生成沿界面传播的能量 集中于界面附近的非均匀平面波，称为界 面波，地震学和地震工程学中称其为面波。 地震面波有瑞利波、拉夫波和斯通利波三 种。
page20
Institute of Disaster Prevention
地震波偏移成像课件
地震波偏移成像的展望
技术创新 随着科技的不断进步，未来将有更多的新技术应用于地震 波偏移成像，如人工智能、大数据等，有望提高偏移成像 的精度和效率。
多学科融合 未来地震波偏移成像将更加注重多学科的融合，如物理学、 数学、计算机科学等，通过多学科的交叉融合，有望突破 现有技术的限制。
应用拓展 随着地震波偏移成像技术的不断成熟，其应用领域也将不 断拓展，如矿产资源勘探、环境监测、城市工程等领域。
地震波偏移成像
03
例分析
实例一：某地区地震波偏移成像结果
总结词
该实例展示了地震波偏移成像技术在某地区的应用，通过偏移成像结果，可以清晰地反映出地下结构的形态和特 征。
详细描述
该实例选取了一个具有复杂地质结构的研究区域，利用地震波探测技术采集地震数据，并采用偏移成像方法对数 据进行处理和分析。最终得到的偏移成像结果能够清晰地反映出地下岩层的分布、界面起伏以及断层结构等信息， 为地质勘探和资源开发提供了重要的依据。
05
总结地震波偏移成像的主要内容
基本原理
应用领域
地震波偏移成像是一种通过分析地震 波在地下的传播规律，反演地下结构 的地球物理方法。它利用地震波在地 下介质中的传播速度差异，通过记录 和分析地震波的传播路径和时间，推 断地下介质的分布和性质。
地震波偏移成像广泛应用于地质勘探、 矿产资源评估、油气田开发、环境工 程、城市地下空间开发等领域。通过 地震波偏移成像，可以揭示地下结构 Байду номын сангаас形态、分布和性质，为相关领域的 科学研究和实践提供重要的数据支持。
地震波数据处理技术
数据预处理
数据预处理包括数据格式转换、 噪声去除、滤波和振幅补偿等， 以提高数据质量和信噪比。
地震波CT成像系统PPT课件
地震波CT成像系统可以用于评估地下矿产 资源的分布和储量，为矿产资源的开发提 供科学依据。
地下水研究
工程地质勘察
地震波CT成像系统可以用于探测地下水的 水位、流向、储量等信息，为水资源管理 和开发提供支持。
地震波CT成像系统可以用于工程地质勘察 ，如隧道、桥梁、建筑等工程的场地勘察 ，为工程设计和施工提供地质依据。
的图像。
面临的挑战与解决方案
数据采集难度大
由于地震波在传播过程中会受到多种因素的影响，导致数 据采集难度较大。解决方案是采用高精度的地震波采集设 备，并优化采集参数。
图像分辨率有限
由于地震波的传播特性，其成像分辨率受到一定限制。解 决方案是采用更高频率的采集设备，以提高图像分辨率。
数据处理时间长
地震波数据处理需要耗费大量的时间和计算资源。解决方 案是采用高性能计算机和优化算法，以提高数据处理速度 。
地震波CT成像系统PPT课件
contents
目录
• 引言 • 系统工作原理 • 系统硬件结构 • 系统性能评估 • 实际应用案例 • 系统优势与挑战 • 未来发展趋势
01 引言
地震波CT成像系统的定义
地震波CT成像系统是一种利用地震波进行地下结构成像的探 测技术。它通过向地下发射地震波，并记录地震波在地下介 质中的传播和反射信息，利用计算机技术对记录的地震波数 据进行处理和分析，以重建地下结构的图像。
04 系统性能评估
分辨率与对比度
分辨率
地震波CT成像系统的分辨率决定了图像中物体细节的识别能力。高分辨率的图 像能够更好地展示物体的细微结构，有助于准确诊断。
对比度
对比度决定了图像中不同组织或结构的区分能力。良好的对比度能够使医生更 容易地区分不同的组织或结构，从而提高诊断的准确性。
地震映像方法PPT课件
3、岩溶塌陷
桂林市中心广场
3、岩溶塌陷
黎塘水泥厂试验剖面
3、岩溶塌陷
溶沟
4、基岩面起伏：利用反射波
贵港某办公楼
4、基岩面起伏：利用反射波
贵港某办公楼
4、基岩面起伏：利用折射波
云南楚雄垃圾场试验剖面
4、基岩面起伏：利用折射波
云南楚雄垃圾场地震映像波形图
5、管线探测：非金属管道
6、大型混凝土构件质量检测
2、反射波
2、反射波
资料解释：
主要也是根据反射波同相轴的变化 反射波的传播时间主要是与界面深度有关
T 4z2 L2 V1
2、反射波
适用条件：
适用于存在波阻抗差异的各种地质条件 当界面深度发生变化时，反射波的传播时间会发
生变化，如在断层两侧表现为突变；如果是倾斜 界面，反射点的位置会偏离记录点向界面的上倾 方向移动。同样可以可根据反射波同相轴的变化 情况定性推断界面的起伏情况。
结束语
当你尽了自己的最大努力时，失败也是伟大的， 所以不要放弃，坚持就是正确的。
When You Do Your Best, Failure Is Great, So Don'T Give Up, Stick To The End
感谢聆听
不足之处请大家批评指导
Please Criticize And Guide The Shortcomings
地震映像法
地震映像方法及其原理 地震映像法的野外工作方法 地震映像方法的应用
地震映像方法及其原理
地震映像法的概念：
地震映像（又称高密度地震勘探和地震多波勘 探），是基于反射波法中的最佳偏移距技术发展 起来的一种常用浅地层勘探方法。
可利用的信息包括：
地震层析成像概论
文档从互联网中收集，已重新修正排版，word格式支持编辑，如有帮助欢迎下载支持。
《地震层析成像概论》大作业张义蜜，01272016-01-04目录1简述用于地震走时成像方法中的射线追踪算法及原理。
............................... 错误!未定义书签。
1.1打靶法........................................................................................................ 错误!未定义书签。
.................................................................................................................. 错误!未定义书签。
.................................................................................................................. 错误!未定义书签。
1.2弯曲(调整)法 ............................................................................................. 错误!未定义书签。
.................................................................................................................. 错误!未定义书签。
.................................................................................................................. 错误!未定义书签。
地震层析成像科普
地震层析成像科普嘿，朋友们！今天咱来聊聊地震层析成像，这可真是个神奇又有趣的玩意儿！你说地震就像地球这个大家伙偶尔发发脾气，那地震层析成像呢，就好比是给地球做了一次超级详细的“体检”。
想象一下，医生给我们做检查，用各种仪器来看看我们身体里面的情况，地震层析成像也是这样，只不过它检查的对象是地球！它是咋工作的呢？就好像我们拿着超级厉害的“透视眼”，能透过地球的层层“皮肉”，看到它里面的结构。
通过接收地震波，就像听地球内部传来的“声音”，然后分析这些“声音”的特点，从而了解地球内部的情况，比如哪里有不同的物质啦，哪里的结构比较特别啦。
这可太重要了呀！为啥这么说呢？你想啊，如果我们能清楚地知道地球内部的结构，那不就像我们知道了家里的布局一样，心里有底呀！对于地震的研究和预测来说，这可是大宝贝呢！它能让我们更好地理解地震是怎么发生的，说不定以后还能更准确地告诉我们啥时候可能会有地震呢，那可就太棒啦！而且啊，地震层析成像的应用可不止于此呢。
它就像一把万能钥匙，能打开好多知识的大门。
比如说，它能帮助我们了解地球的演化历史，就像我们看自己小时候的照片能知道自己是怎么长大的一样。
它还能让我们对地球内部的各种过程有更深入的认识，像地球内部的热传导啦、物质循环啦等等。
你说这是不是很神奇？这就好比我们突然有了一双能看穿地球的眼睛，能看到那些我们以前想都想不到的东西。
它让我们对地球这个神秘的大家伙有了更多的了解，也让我们对大自然的力量充满了敬畏。
所以啊，地震层析成像真的是个了不起的东西呢！它就像一个默默工作的科学家，不断地给我们带来新的发现和惊喜。
咱可得好好感谢那些研究地震层析成像的科学家们，是他们让我们对地球有了更深的认识。
怎么样，现在是不是对地震层析成像有了更清楚的认识啦？是不是觉得它特别厉害？反正我是这么觉得的！希望以后它能给我们带来更多的好消息，让我们能更好地和地球这个大家伙相处呀！。
(第九讲)地震层析成象
Si
其中 S i 表示小区域 Si 之面积， 向量 f f1 , f 2 ,..., f I T被称为图像向量。
三、地震走时层析成像算法
3．1离散图像重建
设射线 Li 与小区域 S i 相交部分之长度为 a ji ，根据 Radon变换，函数 f 沿射线 Li 的投影函数为
j f
二、地震层析成像方法面临的主要问题
2.4地震反演解的可靠性问题 由于震源和检波器位臵分布及连续问题的离散化 , 地震层析成像反演将遇到方程的不适定问题 . 若方 程组是欠定的 , 解可能不存在 , 或者没有唯一解 . 当条件数很大时 , 反演问题将是不稳定的 , 所用算 法也可能不稳定 . 若方程是超定的 , 说明方程组中 的一个或几个方程是其它方程的线性组合 , 或者所 有的方程中某些变量是其它变量的同一线性组合 , 这两种情况都得不出唯一解 ( 杨文采 , 1 997; 刘福 田 , 1989) .
Af τ
三、地震走时层析成像算法
3．2重建算法 方程组（ 3 ）中的系数矩阵 A 是极其稀疏的，因为它 的每一行有 J 各元素，而每条地震波只通过所有 I 个 像元中的一小部分，因此矩阵A中的大部分元素为零。 根据系数矩阵稀疏的特点，对方程组（3）多采用迭 代方法求解。
一、地震层析成像研究发展概况
地震层析成像的研究在 70 年代首先以井间速度结构 调查为研究对象 （ Bois et al.1972 ）。 1979 年， Dines 和 Lytle 首先对地震层析成像坐了大量数值模 拟，并公布了利用弯曲的地震射线进行地下地震波 速 度 成 像 的 结 果 ， 并 首 先 将 层 析 成 像 （ Computerized Geophysical Tomography ）这一名 词用于论文的标题。 1984 年，美国的 Anderson 利用 天然地震数据着手全球构造研究，并公布了全球三 维速度结构。从而使人们对重力场变化、密度结构、 地幔物质流动有了新的认识。
第十四讲地震波层析成像
mantle convection
Travel time table from ak135 model
Travel time picks
Shearer, 2009
3-D variations of Earth’s Structure from Seismic Tomography
Seismic waves in the Earth
Parallel beam
Fan beam, Multi-receiver, Moves in big steps
Cunningham & Jurdy, 2000
Broader fan beam, Coupled, moving source receivers, fast moving
Broader fan beam, Moving source, fixed receivers, fast moving (1976)
Traveltime/waveform
3-D wave speeds Inverse problem
Researchers at MIT and Harvard, led by Keiti Aki and Adam Dziewonski in late 1970’s and 1980’s, pioneered the technique of seismic tomography.
R. Weaver,Science, 2005
Processing Steps:
Remove instrument response, de-mean, detrend, bandpass filter, time-domain normalization, spectral whitening Cross-correlation: 1 day at a time. Stack over many days. Waveform selection (SNR) for tomography
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
全球地震层析成像
10
地幔柱
区域地震层析成像
11
大洋中脊
区域地震层析成像
12
俯冲带
区域地震层析成像
13
局部地震层析成像
14
地震层析成像的基础——Radon 变换 Tomo— Greek for “tomos” (body), graphy --- study or subject Where it all began: Radon transform: (Johan Radon, 1917): integral of function over a straight line segment. Radon transform
Cunningham & Jurdy, 2000
Broader fan beam, Moving source, fixed receivers, fast moving (1976)
Different “generations” of X-Ray Computed Tomography (angled beams are used to increase resolution). Moral: good coverage & cross-crossing 17
第六章 地震层析成像
1
Seismic waves
From IRIS
• Seismic wave is currently the only effective
tool that can penetrate the entire earth
Structural information of the Earth
f (x, y) p(x cos y sin),)d 0
16
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
2
1939:
Jeffreys & Bullen
First travel-time tables：JeffreysBullen Seismological Tables
→ 1D Earth model
Jeffreys-Bullen 1-D Earth Model
3
地球内部结构
1. PREM参考地球模型： Preliminary Reference Earth Model (Dziewonski & Anderson, 1981)
Present Generation of models: Dense receiver sets, all rotating, great coverage and crossing rays.
Brain Scanning Cool Fact: According to an earlier report, the best valentine’s gift to your love ones is a freshly taken brainogram. The spots of red shows your love, not your words!
mantle convection
6
Travel time table from ak135 model
Travel time picks
7
Shearer, 2009
3-D variations of Earth’s Structure from Seismic Tomography
Seismic waves in the Earth
15
Shepp-Logan Phantom (human cerebral)
Input
Radon Projected
Recovered (output)
Back projection of the function is a way to solve f() from p() (“Inversion”):
PREM一维全球速度模型
2. IASP91速度模型 (Kennett and Engdahl, 1991)
3. AK135 速度模型 Kennett et al. (1995)；
Montagner and Kennett (1996).
另外 还有： MC35， STW105， TNA/SNA模型等
4
8
PREM模型给出了地球的一维结构，而地球内部三维结构需要更精细的刻画。 地震层析成像方法是给出地球内部三维结构的最重要的方法。 某种意义上说，地震是照亮地球内部的明灯。 地震层析成像方法可以给出：
全球地球结构的横向不均匀性； 典型地球动力学过程的三维结构：俯冲带、地幔柱、大洋中脊等； 小尺度的构造（断层等）； 地震分布特征。
• • 莫霍诺维齐（前南斯拉夫）——莫霍面（壳幔边界） • －－地壳与地幔分界面的发现者
• 古登堡（美国）——核幔边界 • －－地幔与地核分界面的发现者
• 莱曼（丹麦，女）—— （内外核边界）
•
－－外核与内核分界的！ Plate tectonics
Topography
p(s,) f (x, y)(xcos ysin s)dxdy
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.
Traveltime/waveform
3-D wave speeds
Inverse problem
Researchers at MIT and Harvard, led by Keiti Aki and Adam Dziewonski in late 1970’s
and 1980’s, pioneered the technique of seismic tomography.