地壳与上地幔探测的主要地震学方法

7、接收函数的筛选
接收函数的筛选过程实质是剔除坏数据的过程
筛选原则：接收函数的直达P波幅值要明显， Ps转换震相相对清晰，高频干扰较小。
8、接收函数动校正 把来自不同射线参数的地震事件的Ps转换震相时间延迟校正 到同一参考射线参数的Ps转换震相到时上。
接收函数的主要用途：
1、时深转换
2、接收函数偏移成像
地震学的主要研究方法介绍
一、接收函数方法
接收函数是根据三分量地震数据计算得到的一个时间 序列，它反映了近接收器下方地球介质的响应特性。 接收函数波形是地震台站下方介质扰动所导致的Ps转 换波的组合。
对接收函数的振幅和到时进行模拟，可以给出台站下 方介质非常重要的约束信息。一般情况下，我们将地 球介质近似为水平层状介质，这种情况下，远震体波 穿过地幔，以基本保持不变的水平相速度入射到接收 台站下方时，会在每个速度界面产生透射和强弱不等 的转换波震相以及地表与各速度界面之间的多次反射 震相。
P波之后的后续震相的相对振幅及频率成分取决于速 度转换带的性质，比如速度变化是平缓还是突变？
接收函数提取
Langston[1979]：震源等效化
DV t I t St EV t DR t I t St ER t DT t I t St ET t
远震事件：
EV (t)
DV t I t St
应用新的全球地震层析成像 方法(Zhao, 2001, PEPI),采用多种震相 (P、PP、pP、PcP、Pdiff)
初始速度模型
由标准的IASPEI91模型修改
中国大陆三维地幔速度结构
台站分布 地震分布
2% 0% -2%
中国大陆及邻区
P波速度扰动平面图(1°× 1°)
(Huang & Zhao, JGR, 2006)
36秒的瑞利面波群速度分布图与地壳厚度的对比
青藏高原不同深度S波速度结构
三、体波走时层析成像
Dapeng Zhao, Akira Hasegawa and Hiroo Kanamori. Deep structure of Japan subduction zone as derived from local, regional and teleseismic events. Journal of Geophysical Research, 1994, 99(B11):22313-22329.
3、低通滤波 对地震事件进行5Hz的低通滤波。
5Hz低通 滤波
0.25Hz低 通滤波
0.1Hz低 通滤波
4、截取P波波形 一般截取P波到达之前20秒、到达之后80秒进行计算。
5、坐标旋转 将N-S、E-W分量旋转到R、T分量。 Ps转换波在R分量上振幅最强。
6、接收函数计算 一般采用频率域反褶积方法计算。
四、震源机制
五、地震定位
Hypoinverse2000 HypoDD Joint Inversion
（Xiaofeng Liang，2009）
六、有限频层析成像
Finite Frequency Tomography
有限频率层析成像技术：最近几年沈旸及其科研团队利用有限 频率地震波的敏感核开发了一种新的层析成像技术。这是过去 20年来地震层析成像技术研究和开发领域上取得的一个重大突 破，从而使这一技术不再依赖于传统上的射线理论。射线理论 只对于无限频率假设才成立，而观测到的地震波都是有限频率 的。这种新方法可以考虑波阵面复原以及其它地震波的绕射效 应。结果显示，这一技术对于地震结构的成像有非常大的改善， 可以应用于从勘探到地震灾害研究等各种层析成像领域。
Байду номын сангаас
太平洋俯冲带地区上地幔结构
3、地壳厚度和泊松比变化
4、S波速度结构的反演
二、利用背景噪声研究地下结构
地震噪声是由多种不同的,空间分布互不相 关的连续源所产生的杂乱的记录,是一种平 稳随机过程,并且没有确定的相谱.
最近的科学研究表明，噪声中包含有大量有 用的信息，是完全可以利用的。
• 通过对长时间段内的噪声进行互相关运算, 可以提取接收点间的格林函数. 其实验证 明已经在太阳地震学和声学领域给出,在海 洋声学以及地震学中开始初步应用.
Crust and Lithosphere: Multi-resolution surface wave tomography
Seismic interferometry estimate data from background “noise”
背景噪声成像研究结果
16秒的瑞利面波群速度分布图于沉积层厚度的对比
如果三分量地震仪的脉冲响应一致:
ER ()
DR () I ()S()
DR () DV ()
ET
()
DT () I ()S()
DT DV
() ()
接收函数数据处理流程：
1、数据格式转换 将seed格式的数据或Reftek格式的数据转换为sac格 式，需要在LINUX下进行。
2、地震事件选取 选择震中距在30—95°、震级大于5.5的地震事件。
Yang Ting,2006, EPSL
七、各向异性研究
横波分裂 方位各向异性
径向各向异性
Silver and Chan,1991
Rayleigh波群速度横向不均匀性和 方位各向异性图像
• 本文主要介绍其在地震学领域中的发展. 通过对地面上两个地震台站记录的背景地 震噪声信号在长时间段内进行互相关运算, 提取格林函数,从而进行地震成像,获取对 地下结构的认识.
背景噪声成像方法流程图
叠加时间长度对相 关结果的影响。瑞 利面波的出现随着 时间序列叠加的增 长而趋于明显。
面波的频散特征： 不同周期的波有着 不同的传播速度， 长周期的波快于短 周期的波。
远震在接收介质中的传播射线路径示意图及接收响应
接收函数中各个震相的振幅取决于入射P波的入射角 及产生Ps转换波及多次转换波（PsPs, PpSs, PpPs） 的速度差异的强弱。
转换震相和多次震相的到时取决于如下的因素：速度 界面的深度，速度界面与地表之间P波与S波的速度， P波入射角或射线参数。