东亚及西太平洋地壳上地幔P波和S波三维速度结构.pdfVIP

东亚及西太平洋地壳上地幔P波和S波三维速度结构 瞿辰周蕙兰 (中国科学院研究生院计算地球动力学重点实验室北京100049) 本文研究的区域是东亚及西太平洋地区，为经度60。E至140。E、纬度10。s至60。N范围内的区域。该 区域内的地貌和地质构造以其复杂多样性而倍受研究者关注。同时，这种浅部的复杂性与深部复杂性的关联程度， 或浅部的复杂性延伸到了多大深度，更是引起很多学者的兴趣。而深部的复杂性的研究目前还主要靠地震波。地 震波给出的三维速度结构，在一定程度上可以揭示出介质弹性性质的不均匀性，进而被用来推断物性和温度的变 化，以及探讨引起这些变化的更深入的原因。该区域深部的结构复杂性及其它多方面的研究工作虽已有不少成果， 但关于其地壳上地幔结构的区域性地震层析成像仍需要进一步做工作。本文收集了最新的P波和s波的到时资料， 采用地震走时层析成像方法，反演该地区深度从0至800km范围内的三维V，，V。和V。。，，V。结构。 der 范围内的P波和S波到时数据。其中EHB数据是由由Engdahl，Vall 行地震重新定位后得到的，其震源基本参数更加准确。因此本文使用EHB提供的震源参数，ISC报告、中国基本测 震台网和NEIC台网提供的台站信息，并剔除ISC报告和中国台网报告中重叠的地震记录。这样，得到P波震相数据： 在同一个震源点，在不同时间可能发生多次地震，而这些地震都可能被相同的台站接收。所以，实际地震记 录中存在大量重复射线。在迭代解法中这些重复射线会造成对相同网格节点的波速多次修改，进而造成反演结果 中的局部奇异点；同时占用了大量内存和计算时间，影响计算速度。所以，必须进行射线的归并。本文将地震记 录中震源和接收台站位置相近的射线归类，选取其中走时残差最小的记录。因此最后归并的结果既保证了尽可能 多的地震射线被使用，又大大减少了路径几乎相同的射线。归并后得到由2946个台站记录到132302个事件的 1668133条P波数据。 S波的数据处理方法与P波完全相同。 本文主要采用Zhao等人(1994)的层析成像方法(嗍D)。该崩壁绷j网格节点侈鲤参数化，正演计算采用近似弯曲 快速射缀醵亲i去；反演采用DLSQR法。首先将研究区速度摸型参数‰从地袭至800虹深度，设置了3维网格节氨水平 方向上以2。X l血深度划 分。模型中只设置了两个间断爵一康氏界面和莫霍界面。初始速度模型取自iasp91速度模型，E述两个间断面深度取自 Bassin等人(2000)的2。×2。的全球地壳模型CRUST2．0。 因为在东巫及西太平洋地区莫霍晃面的起伏非常大，所以本文的初始模型中某些水平向节点层不可避免的与间 断面相交。在同一深度层上的节点(如30km、60km、70k击深度的节点)可能处在地壳地幔的不同位置，故它们的 初始速度值也不同。据此，初始速度模型放置成为简单的三维模型。 从初始三维模型起算，使用近f以弯曲快速射线追踪法CZhao，1992)计拿陆擞路径和理论走时，用DLSOR(Paige， Sauders。1982)方法求解反演方程，分别得到P波和S波的速度模型，并计算Vr／v。的分布。 结果的可靠性分析主要从检测板分辨率实验结果、节点周阑的射线分布和反演结果中的速度误差分布等兰个方面进行分 析。从得到的结构模型沿横向和纵向截取速度剖面，并对各们塞度剖面进行解释，综合地质学，地球物理，地球化学等多学 科多手段对层析成像结果进行解释，进—步研究东亚及西太平洋地区速度结构及其地质地球物理意义。 本项目得到国家自然科学基金资助(项目批准40374009) · · 387