地球物理实验报告1用接收函数反演地壳结构.docxVIP

用接收函数反演地壳结构——P陆文玮一、实验目的本实验我们主要为了了解接收函数的概念，如何从地震数据中提取接收函数，如何利用得到的一维分层模型计算(理论)接收函数，以及利用已有程序进行简单的反演过程计算。实验程序为CharlesJ.Ammon,PennState,1997 写的，详细内容在下面这个网址里查询到。/~cammon/HTML/RftnDocs/rft n01.html 二、实验原理1.接收函数的原理：远震体波受到震源时间函数 、传播路径 、接收台站下方的介质结构以及仪器响应等多种复杂因素的共同作用 ,而从地震波中提取详细的地球内部结构一直是地震学研究的目标之一。接收函数是用远震P波波形的垂直分量对径向分量和切向分量作反褶积处理后得到的时间序列。简单来讲就是去除震源 、地震波传播路径以及仪器响应等因素后的时间序列。它主要包含地震台站下方地壳和上地幔速度间断面所产生的转换波及其多次反射波的信息。如上图所示，可以通过远震P波波形（接收函数）来反演右边的S 波波速结构。2.接收函数的意义：P 波接收函数可以非常有效地获取地壳内部S 波速度结构、地壳厚度及物质成分组成和地幔过渡带的厚度变化。通过接收函数的联合反演方法可以非常有效地降低了结果的非唯一性。同时，利用P波接收函数研究壳幔结构的横向非均匀性。三、实验内容1. 使用pwaveqn程序从数据中提取接收函数原始数据:uln_1995_113_b_057_d_045.[enz]，如图所示：图一、原始数据图输出: ULN_TEST.aftn ULN_TEST.eqr ULN_TEST.eqt，如图所示：图二、处理图2.模型计算并计算模型对应的接收函数模型计算我们采用的程序是icmod，对于最简单的情况选第5 项,使用默认泊松比,即每层只用输入vp 和层厚。然后我们再计算模型对应的接收函数，采用的程序是respknt，输出文件为smp.mod_sp.[trz]，我们同样需要用pwaveqn计算它们的接收函数。3.反演程序的使用反演程序有snglinv, smthinv, manginv, 我们采用的反演程序是snglinv。初始模型和最终模型如图所示：图三、初始和最终模型的vp值由初始模型到最终的模型我们一共设置了5 个迭代步骤，中间迭代过程中的模型如图所示：在实验的过程中使用了多次迭代反演图五、迭代波形图只选择初始与最终的波形图画出来相比较图六、初始与最终的比较图从以上的图可以看出：初始模型和最终反演结果的差异还是很大，特别是在浅部（20km），修正幅度特别大。四、结论我们采用了单个台站的数据，对其下方的速度结构进行了反演。通过这一过程，从数据出发得到了地球内部的信息。所得模型能较好地解释观测数据。思考题：如果射线图中所示的界面并非水平，而是倾斜的，会对结果有何种影响？利用接收函数反演地壳结构并不要求界面是水平的，因此界面倾斜应该无影响。