非均匀各向异性介质的地震P波走时层析成像研究

国际地震动态　２０１０经　非均匀各向异性介质的地震　Ｐ波走时层析成像研究　郭飚　１０００２９）　（中国地震局地质研究所，北京中图分类号：Ｐ３１５；　文献标识码：Ａ；ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２３５－４９７５．２０１０．０３．０１４　大陆动力学已经成为当今固体地球物理各领域研究的主导方向。大陆动力学涉及问题　非常广泛，但核心问题是大陆形变及其动力学。地震各向异性是地球动力学过程的指示器，　根据地震各向异性的研究结果，可以推断上地幔物质的流动或变形，有助于了解地球内部的　应力状态和地球的演化过程。　地震层析成像是研究地球地壳上地幔结构的主要地球物理方法。地震观测技术的进步、　数据的积累以及计算技术的发展使得利用体波走时数据反演各向异性参数成为可能。利用地　震走时层析成像方法同时研究介质的非均匀性和各向异性，对于认识地球的结构及动力学过　程都具有非常重要的意义。　本文系统地讨论了弱各向异性条件下Ｐ波速度扰动的近似表达式，并在此基础上构建了　地震各向异性Ｐ波走时层析成像算法。利用多个不同理论模型对本文提出的方法进行了数值　检验，系统地讨论了该方法的应用条件，并将本文提出的理论和方法成功地应用在中国境内　天山地区和川西龙门山地区，获得了上述地区地壳上地幔的Ｐ波速度结构和Ｐ波各向异性参数。　本论文的研究成果主要包括以下４个部分：　１　理论构建　在Ｂａｃｋｕｓ¨　的弹性张量调和展开理论基础上，导出了弱各向异性介质中的Ｐ波相速度扰　动调和展开的的各阶系数，校正了Ｓｍｉｔｈ＆Ｄａｈｌｅｎ－２　和Ｂｏｋｅｌｍａｎｎｌ３　文献中的错误。根据Ｐ波　速度扰动的调和展开公式，导出了在不同观测系统和不同各向异性介质情况下的简化公式。　根据Ｂｏｎｄ旋转矩阵证明了绕坐标轴旋转时弹性矩阵元素变化规律，证明了在正交对称系垂　直对称轴情况下，仅需４个参数即可模拟方位各向异性，而在任意各向异性情况下，仅需６　个参数即可描述三维各向异性。　２技术实现　基于弱各向异性近似条件下弹性张量的调和展开理论和地震走时层析成像方法，构建了　地震各向异性Ｐ波走时层析成像算法。该算法利用均匀的网格节点和三次Ｂ样条插值函数描　述各向同性速度场；利用各向异性块描述各向异性速度扰动场。这种模型参数化方法的优点　在于：（１）模型光滑，在反演时不需要加入光滑项；（２）各向同性和各向异性的分辨率可以　分别考虑，以减少介质的非均匀性与各向异性间的耦合。另外，总结对比了地震波走时计算　的各种方法，构建了各向同性介质中快速行进波前扩展程序，并在此基础上构建了适合弱各　●　第３期　中国地震局地质研究所２００９届博士学位论文摘要　４７　向异性介质的各向异性快速行进算法，利用互易性定理减少了相应的走时计算的计算量。在　该算法中，首先计算以台站位置为假想震源的走时场，再计算台站与近震震源或远震射线与　模型边界的交点间的射线路径。这种走时计算方式特别适用于多源多台站的观测方式，可有　效地减少正演计算量。反演方法采用Ｔａｒａｎｔｏｌａ＆Ｖａｌｅｔｔａ　和Ｔａｒａｎｔｏｌａ　的非线性反演方法，　对于非均匀和各向异性参数可以分别设置不同的阻尼因子。　３数值检验　本文给出的多个理论模型测试结果表明：（１）本文给出的各向同性和弱各向异性条件下　的快速行进算法的计算精度可以满足地震走时层析成像的要求（误差小于震相的检测误差）；　（２）各向异性的强度与非均匀结构是耦合的，而各向异性方向与非均匀结构是可以解耦的；　（３）在数据分布较好的情况下，Ｐ波走时方法可以分辨多层各向异性；（４）对于远震，只要　数据的方位覆盖大于１８０。即可较好地反演各向异性参数；（５）三维各向异性的确定需要大于　５０。的震中距覆盖。　４实际数据解释　（１）利用横跨中国境内天山的库车一奎屯宽频带流动地震台阵和区域地震台网记录的近　震和远震Ｐ波走时数据和Ｚｈａｏ等　’　地震层析成像方法，重建了沿该地震台阵剖面４００　ｋｍ　深度范围内地壳上地幔的各向同性Ｐ波速度结构。在此基础上，利用弱各向异性介质条件下　的Ｐ波速度扰动公式和线性反演方法，进一步解释了剩余残差，获得了该地震剖面上的Ｐ波　快波方向。我们的结果表明：　（ａ）沿新疆库车一奎屯剖面，天山地壳具有明显的横向分块结构，且南、北天山地壳显　示了较为强烈的横向变形特征，表明塔里木地块对天山地壳具有强烈的侧向挤压作用。在塔　里木和准噶尔地块上地幔顶部有厚度约６０～９０　ｋｍ的高速异常体，塔里木一南天山下方的高　速异常体产生了较为明显的弯曲变形，而准噶尔一北天山下方的高速异常体向南一直俯冲到　中天山南侧边界下方３００　ｋｍ的深度。两者形成了不对称对冲构造。在塔里木和准噶尔地块下　方１５０～４００　ｋｍ深度存在上地幔低速体。其中，塔里木地块一侧的上地幔低速体上涌到南天　山地块的下方。在塔里木一南天山２００～３００　ｋｍ深度范围的上地幔存在高速异常体，它可能　是地幔热物质向上迁移过程融断的塔里木岩石圈拆离体。　（ｂ）塔里木地块的俯冲可能涉及整个岩石圈深度，但其前缘仅限于南天山的北缘；青藏　高原隆升的远程效应可能驱动了塔里木岩石圈向北俯冲，同时还造成天山造山带南侧上地幔　物质的涌入。天山造山带上地幔广泛存在的低速异常有助于上地幔的变形。上地幔物质的强　烈非均匀性应有助于推动天山造山带上地幔范围小尺度地幔对流的形成。根据研究区地壳上　地幔速度结构特征可以推断，新近纪以来，天山快速隆升的主要力源来自青藏高原快速隆升　的远程效应，相对软弱的上地幔为加速天山造山带的变形和隆升创造了必要条件。　（ｃ）塔里木盆地和南天山的Ｐ波快波方向为近南北向，这与塔里木地块对天山南北向挤　压变形作用是一致的。进入中天山褶皱带和北天山，Ｐ波快波方向逐渐偏转为北西向，北天　山山前褶皱带和准噶尔盆地前缘的Ｐ波快波方向为近东西向，总体上，该区Ｐ波快波轴方向　平行于天山造山带的走向，垂直于板块相对运动方向。塔里木盆地和南天山的壳幔为连续一　致变形，而中天山、北天山和准噶尔上地幔的各向异性方向可能受到地幔小尺度对流和构造　国际地震动态　２０１０年　方向的影响，壳幔的形变发生解耦。这表明，由于两侧盆地的挤压，天山山体处于挤压状态，　而上地幔物质的运移沿着山体平行的方向，两者之间构成了立交模型。根据研究区地壳上地　幔速度结构特征和各向异性方向推断，上地幔物质的流动变形在天山造山带的变形和隆升过　程中起着重要作用。　（２）利用在龙门山及邻区布设的川西地震台阵记录到的远震Ｐ波走时数据和本文提出的　各向异性层析成像方法，同时反演了该区地壳上地幔远震Ｐ波各向同性速度结构和Ｐ波快波　方向。我们的结果表明：　（ａ）研究区地壳上地幔Ｐ波速度结构具有较为明显的分区特征，松潘一甘孜地块和川滇　地块速度较低，龙门山断裂带及四川盆地西部速度较高。四川盆地西部地壳上地幔的高速异　常厚度从南侧２５０　ｋｍ向北逐渐减薄至１００　ｋｍ，推测这个高速异常体可能代表四川盆地的岩　石圈。相对于四川盆地，松潘一甘孜地块的地壳上地幔较为软弱，但不存在四川盆地向松　潘一甘孜地块的俯冲。松潘一甘孜地块的抬升可能与地幔上涌有关。四川盆地与松潘一甘孜　地块和川滇地块间的动力学机制完全不同，川滇地块和四川盆地仅是垂直接触关系，而在四　川盆地北部靠近龙门山一侧，发现四川盆地前缘自东向西减薄的现象，推测松潘一甘孜的上　地幔物质侵蚀了四川盆地下方的岩石圈。鲜水河断裂带和龙门山断裂带都为超壳的深大断　裂，鲜水河断裂带在深部结构上显示为较宽的低速特征。　（ｂ）龙门山断裂带与四川盆地的地壳均显示为高速特征，四川盆地与松潘一甘孜地块之　间的构造边界可能是汶川一茂县断裂。以汶川为界，龙门山断裂带被从松潘～甘孜地块侵入　的低速物质分为两段，南北两段均显示为高速特征。汶川大地震及余震序列均发生在龙门山　北段的高速介质区域内。这种深部结构特点对汶川大地震孕育和余震序列发育过程具有明显　的控制作用。根据本文给出的松潘一甘孜地壳上地幔速度结构，我们可以进一步推断，松　潘一甘孜地块的抬升应与地幔物质上涌有关，在坚硬的四川盆地的阻挡作用下，青藏高原向　东挤压和地幔上涌的双重作用造成松潘一甘孜地块隆升，并造成了龙门山断裂带的逆冲推　覆。龙门山南北两段间的地壳强度较低，在长期的缓慢变形过程中，易于在龙门山北段的坚　硬的上地壳内形成巨大的应力积累，而汶川处于龙门山北段的最南端，应力容易在此集中，　并成为汶川大地震破裂的起始点。　（ｃ）Ｐ波各向异性研究表明川滇地块的上地幔流动方向与ＧＰＳ的观测结果相吻合，反映　了川滇地块地壳和上地幔的形变一致，属于壳幔耦合型。松潘一甘孜地块的上地幔流动方向　为北东向，而ＧＰＳ的观测结果为近东西向，说明松潘甘孜地块壳幔的形变不一致，属于壳幔　解耦型。四川盆地北部的Ｐ波快波方向与ＳＫＳ的结果较为一致，而四川盆地南部的Ｐ波快波　方向平行于主要构造方向，与ＳＫＳ的结果相互垂直。四川盆地内部南北两部分具有不同的壳　幔变形特征，这与速度结构给出的结论相一致。　（ｄ）本文的研究结果并不支持四ＪＩｉ盆地向松潘一甘孜地块的俯冲的动力学模型。与其结　果相反，上地幔的物质流动在该区的演化变形过程中起主导作用，在上地幔深度上，松潘一　甘孜地块的上地幔物质可能侵入四川盆地下方。　关键词地震各向异性；地震层析成像；弹性张量的调和展开；地震台阵；天山造山带；　龙门山断裂带；汶川地震　（作者电子信箱，郭飚：Ｇｕｏｂｉａｏ７４＠ｈｏｔｍａｉｌ．ｃｏｍ）　（下转第２２页）　国际地震动态　２０１０年　络的地震台网手册（包括临时布局），使每一　于网络的帮助平台，提供有关地震参数数据　台网的资料都得以清晰显示，并对其从前和　的处理和上报的所有必要的指导性文档。主　现在的状况、分析流程、职员、目前的问题　持编写《地震观测实践手册》也是ＩＳＣ对地震　及发展规划等基本信息加以表述。　学方面的教育工作所作出的一份特殊贡献。　（２）重新设计ＩＳＣ台站手册，使其包含　（５）向信誉卓著的机构提供ＩＳＣ通报选　（６）向数据贡献者提供他们所关心的领　图片以及ＩＳＣ数据库中可利用的参数数据等　集，供其审阅。　简要信息。可对每一台站报告的参数数据进　行分析。　域的及时反馈。　（３）建立一个隶属ＣＯＳＯＩ／ＩＡＳＰＥＩ之下、　与美国国家地震信息中心（ＮＥＩＣ）共同管理　资料来源：ｔｉｐ：／／ｆｔｐ．ｉｓｃ．ａｃ．ｕｋ／ｐｕｂ／ｄｏｃｓ／ＩＳＣ—　的新的“国际地震网络台站档案”。该档案可　ｄｅｖｅｌｏｐ＿０２０９－ｌ１．ｐａｌ　使ＩＳＣ对网络所有人和运营者以及参数数据　原题：Ｓｔｒａｔｅｇｉｅｓ　ｆｏｒ　ｔｈｅ　ＩＳＣ　Ｄｅｖｅｌｏｐｍｅｎｔ　２００９－　编写者和发布者等所作出的贡献给予肯定。　２０１ｌ　有了这样的档案，ＩＳＣ还可以有效地应对常　（中国地震局地球物理研究所左玉玲译）　见的情况，即根据不同的数据提供者报告给　（译者电子信箱，左玉玲：ｙｕｌｉｎｇｚｕｏ＠ｙａｈｏｏ．　ｃｏｉｎ．ｃａ）　ＩＳＣ的相同免费波形复制参数数据等。　（４）与ＣＯＳＯＩ／ＩＡＳＰＥＩ联手设立一个基　（上接第４８页）　参考文献　［１］Ｂａｃｋｕｓ　Ｇ　Ｅ．Ａ　ｇｅｏｍｅｍｃａｌ　ｐｉｃｔｕｒｅ　ｏｆ　ａｎｉｓｏｔｒｏｐｉｃ　ｅｌａｓｔｉｃ　ｔｅｎｓｏｒｓ．Ｒｅｖ．Ｇｅｏｐｈｙｓ．Ｓｐａｃｅ　Ｐｈｙｓ．，１９７０，８：６３３—　６７１　［２］Ｓｍｉｔｈ　Ｍ　Ｌ，Ｄａｈｌｅｎ　Ｆ　Ａ．Ｔｈｅ￣ｉｍｕｔｈａｌ　ｄｅｐｅｎｄｅｎｃｅ　ｆｏ　Ｌｏｖｅ　ａｎｄ　Ｒａｙｌｅｉｇｈ　ｗａｖｅ　ｐｒｏｐａｇａｔｉｏｎ　ｉｎ　ａ　ｓｌｉ【ｓｈｆｌｙ　ａｎｉｓｏ－　ｔｒｏｐｉｃ　ｍｅｄｉｕｍ．Ｊ．Ｇｅｏｐｈｙｓ．Ｒｅｓ．，１９７３，７８（１７）：３　３２１－３　３３３　［３］Ｂｏｋｅｌｍｅｎｎ　Ｇ　Ｈ．Ｃｏｎｖｅｃｔｉｏｎ－ｄｒｉｖｅｎ　ｍｏｔｉｏｎ　ｏｆ　ｔｈｅ　Ｎｏｒｔｈ　Ａｍｅｉｒｃｎａ　ｃｒａｔｏｎ：Ｅｖｉｄｅｎｃｅ　ｆｒｏｍ　Ｐ－ｗａｖｅ　ａｎｉｓｏｔｒｏｐｙ．　Ｇｅｏｐｈｙｓ．Ｊ．Ｉｎｔ．，２００２，１４８：２７８－２８７　［４］Ｔａｒａｎｔｏｌａ　Ａ，Ｖａｌｅｔｔｅ　Ｂ．Ｇｅｎｅｒａｌｉｚｅｄ　ｎｏｎｌｉｎｅａｒ　ｉｎｖｅｒｓｅ　ｐｒｏｂｌｅｍｓ　ｓｏｌｖｅｄ　ｕｓｉｎｇ　ｔｈｅ　ｌｅａｓｔ　ｓｑｕａｒｅｓ　ｃｒｉｔｅｒｉｏｎ．Ｒｅｖ．　Ｇｅｏｐｈｙｓ．Ｓｐａｃｅ　Ｐｈｙｓ．，１９８２，２０：２１９－２３２　［５］Ｔａｒａｎｔｏｌａ　Ａ．Ｉｎｖｅｒｓｅ　Ｐｒｏｂｌｅｍ　ｈＴｅｏｒｙ：Ｍｅｔｈｏｄｓ　ｆｏｒ　Ｄａｔａ　Ｆｉｔｔｉｎｇ　ａｎｄ　Ｍｏｄｅｌ　Ｐａｒａｍｅｔｅｒ　Ｅｓｔｉｍ￣ｉｏｎ．Ａｍｓｔｅｒｄａｍ：　Ｅｌｓｅｖｉｅｒ，１９８７　［６］Ｚｈａｏ　Ｄ　Ｐ，Ｈａｓｅｇａｗａ　Ａ，Ｈｏｒｉｕｃｈｉ　Ｓ．Ｔｏｍｏｇｒａｐｈｉｃ　ｉｍ　ｎｇ　ｆｏ　Ｐ　ａｎｄ　Ｓ　ｗａｖｅ　ｖｅｌｏｃｉｔｙ　ｓｔｒｕｃｔｕｒｅ　ｂｅｎｅａｔｈ　Ｎｏｒｔｈｅｓａｔ—　ｅｖｉｌ　Ｊａｐａｎ．Ｊ．Ｇｅｏｐｈｙｓ．Ｒｅｓ．，１９９２，９７（Ｂ３）：１９９０９－１９９２８　［７］Ｚｈａｏ　Ｄ　Ｐ，Ｈａｓｅｇａｗａ　Ａ，Ｋａｎａｍｏｒｉ　Ｈ．Ｄｅｅｐ　ｓｔｒｕｃｔｕｒｅ　ｆｏ　Ｊａｐａｎ　ｓｕｂｄｕｅｔｉｏｎ　ｚｏｎｅ　ａｓ　ｄｅｒｉｖｅｄ　ｆｒｏｍ　ｌｏｃａｌ，ｒｅｇｉｏｎａｌ，　ｎａｄ　ｔｅｌｅｓｅｉｓｍｉｃ　ｅｖｅｎｔｓ．Ｊ．Ｇｅｏｐｈｙｓ．Ｒｅｓ．，１９９４，９９（Ｂ１１），２２３１３－２２３２９