郯庐断裂带中南段及邻区Pn波速度结构与各向异性

书 关于书的成语关于读书的排比句社区图书漂流公约怎么写关于读书的小报汉书pdf 第５９卷第２期２０１６年２月地　球　物　理　学　报ＣＨＩＮＥＳＥ　ＪＯＵＲＮＡＬ　ＯＦ　ＧＥＯＰＨＹＳＩＣＳＶｏｌ．５９，Ｎｏ．２Ｆｅｂ．，２０１６顾勤平，丁志峰，康清清等．２０１６．郯庐断裂带中南段及邻区Ｐｎ波速度结构与各向异性．地球物理学报，５９（２）：５０４５１５，ｄｏｉ：１０．６０３８／ｃｊｇ２０１６０２１０．ＧｕＱＰ，ＤｉｎｇＺＦ，ＫａｎｇＱＱ，ｅｔａｌ．２０１６．ＰｎｗａｖｅｖｅｌｏｃｉｔｙａｎｄａｎｉｓｏｔｒｏｐｙｉｎｔｈｅｍｉｄｄｌｅｓｏｕｔｈｅｒｎｓｅｇｍｅｎｔｏｆｔｈｅＴａｎＬｕｆａｕｌｔｚｏｎｅａｎｄａｄｊａｃｅｎｔｒｅｇｉｏｎ．犆犺犻狀犲狊犲犑．犌犲狅狆犺狔狊．（ｉｎＣｈｉｎｅｓｅ），５９（２）：５０４５１５，ｄｏｉ：１０．６０３８／ｃｊｇ２０１６０２１０．郯庐断裂带中南段及邻区犘狀波速度结构与各向异性顾勤平１，２，丁志峰１，康清清２，赵启光２１中国地震局地球物理研究所，北京　１０００８１２江苏省地震局，南京　２１００１４摘要　郯庐断裂带是一条纵贯我国大陆东部ＮＮＥ走向的巨型深断裂，其中南段及邻区（１１５°Ｅ—１２２°Ｅ，２９°Ｎ—３８°Ｎ）跨越了华北断块区、扬子断块区和华南褶皱系三大一级构造单元，由于其重要性和复杂性，长期以来一直是地学家们研究的热点．本文从国际地震中心（ＩＳＣ）、中国地震台网及区域地震台网的地震观测报告中精心挑选出６３８１个Ｐｎ震相数据，用Ｐｎ波时间项层析成像法反演得到了郯庐断裂带中南段及邻区上地幔顶部Ｐｎ波速度结构和各向异性．结果显示，研究区上地幔顶部具有显著的横向非均匀性，相对于７．９５ｋｍ·ｓ－１的平均速度而言，Ｐｎ波速度值在７．６８～８．２４ｋｍ·ｓ－１范围内变化．Ｐｎ波速度分布在郯庐断裂带中段和南段具有分段性：沿中段及周边存在一ＮＥ向低速异常带，低速可能是由于岩石圈的减薄和软流圈的高温物质沿郯庐带上涌导致；沿南段 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 现为一ＮＮＥ向弱高波速异常带，作为高低速的边界带清晰地勾勒出了华北与扬子这两个不同块体，该边界在江苏域向华北地块ＮＷ方向凹进．Ｐｎ波速度各向异性的强弱与速度分布存在一定的相关性．总体上，如鲁西隆起及以南等低速区、茅山断裂附近的高低速过渡带，其速度各向异性较为强烈；而在具有高速异常的苏北盆地、合肥盆地等稳定区域下方其各向异性较弱．本文通过Ｐｎ波震相基本未能探测到郯庐断裂带中段的方位各向异性，推测是上地幔顶部被“冻结”下来的各向异性痕迹被软流圈热物质上涌这一强烈构造运动削弱所导致．南段具有与断裂伸展方向近乎平行的快波速方向．Ｐｎ波速度横向变化和强震活动存在一定关联．强震主要发生在Ｐｎ波低速异常区或高低速过渡带上．郯城８．５级地震震中位于中段和南段高低速过渡带，该区域也是速度横向变化最大的地方，最容易集中应力和产生应力差．关键词　郯庐断裂带；中南段；Ｐｎ波速度；各向异性；郯城地震ｄｏｉ：１０．６０３８／ｃｊｇ２０１６０２１０中图分类号　Ｐ３１５收稿日期２０１５０８０１，２０１５１２１４收修定稿基金项目　江苏省地震局青年科学基金项目“江苏及邻区上地幔顶部Ｐｎ波速度及各向异性研究”（１４４０４）和中国地震局震害防御司“縀茶河断裂活动性鉴定”项目（１５３０６０１）联合资助．作者简介　顾勤平，男，中国地震局地球物理研究所在读博士研究生，主要从事天然地震层析成像、浅层人工地震勘探及地震安全性评价研究工作．Ｅｍａｉｌ：ｇｑｐ１２２１＠１６３．ｃｏｍ通讯作者　丁志峰，男，研究员，主要从事地震学方面的研究．Ｅｍａｉｌ：ｄｉｎｇｚｆ＠ｃｅａｉｇｐ．ａｃ．ｃｎ犘狀狑犪狏犲狏犲犾狅犮犻狋狔犪狀犱犪狀犻狊狅狋狉狅狆狔犻狀狋犺犲犿犻犱犱犾犲狊狅狌狋犺犲狉狀狊犲犵犿犲狀狋狅犳狋犺犲犜犪狀犔狌犳犪狌犾狋狕狅狀犲犪狀犱犪犱犼犪犮犲狀狋狉犲犵犻狅狀ＧＵＱｉｎＰｉｎｇ１，２，ＤＩＮＧＺｈｉＦｅｎｇ１，ＫＡＮＧＱｉｎｇＱｉｎｇ２，ＺＨＡＯＱｉＧｕａｎｇ２１犐狀狊狋犻狋狌狋犲狅犳犌犲狅狆犺狔狊犻犮狊，犆犺犻狀犪犈犪狉狋犺狇狌犪犽犲犃犱犿犻狀犻狊狋狉犪狋犻狅狀，犅犲犻犼犻狀犵１０００８１，犆犺犻狀犪２犈犪狉狋犺狇狌犪犽犲犃犱犿犻狀犻狊狋狉犪狋犻狅狀狅犳犑犻犪狀犵狊狌犘狉狅狏犻狀犮犲，犖犪狀犼犻狀犵２１００１４，犆犺犻狀犪犃犫狊狋狉犪犮狋　ＴｈｅＴａｎｃｈｅｎｇＬｕｊｉａｎｇｆａｕｌｔｚｏｎｅｉｓｔｈｅｌａｒｇｅｓｔｄｅｅｐＮＮＥｔｒｅｎｄｉｎｇｆａｕｌｔｚｏｎｅｉｎｔｈｅＥａｓｔｅｒｎＣｈｉｎａ，ｗｈｉｃｈｈａｓａｌｗａｙｓｂｅｅｎｆｏｃｕｓｅｄｂｙｍａｎｙｇｅｏｓｃｉｅｎｔｉｓｔｓｂｅｃａｕｓｅｏｆｉｔｓｉｍｐｏｒｔａｎｃｅａｎｄｃｏｍｐｌｅｘｉｔｙ．ＴｈｅｍｉｄｄｌｅｓｏｕｔｈｅｒｎｓｅｇｍｅｎｔｏｆｔｈｅＴａｎＬｕｆａｕｌｔａｎｄｓｕｒｒｏｕｎｄｉｎｇａｒｅａｓ（１１５°Ｅ—１２２°Ｅ，　２期顾勤平等：郯庐断裂带中南段及邻区Ｐｎ波速度结构与各向异性２９°Ｎ—３８°Ｎ）ｉｓｌｏｃａｔｅｄｉｎｔｈｅｃｒｏｓｓｉｎｇｐａｒｔｏｆｔｈｒｅｅＣｈｉｎｅｓｅｇｅｏｌｏｇｉｃａｌｂｌｏｃｋｓ（ＮｏｒｔｈＣｈｉｎａｂｌｏｃｋ，ＹａｎｇｔｚｅＣｒａｔｏｎ，ＳｏｕｔｈＣｈｉｎａＦｏｌｄｓ）．Ｗｅｃａｒｅｆｕｌｙｓｅｌｅｃｔｅｄ６３８１ＰｎｔｒａｖｅｌｔｉｍｅｐｉｃｋｓｆｒｏｍｔｈｅｅａｒｔｈｑｕａｋｅｂｕｌｅｔｉｎｓｏｆｔｈｅＩｎｔｅｒｎａｔｉｏｎａｌＳｅｉｓｍｏｌｏｇｉｃａｌＣｅｎｔｅｒ（ＩＳＣ），ＣｈｉｎａｎａｔｉｏｎａｌｓｅｉｓｍｉｃｎｅｔｗｏｒｋａｎｄｒｅｇｉｏｎａｌｓｅｉｓｍｉｃｎｅｔｗｏｒｋｔｏｉｎｖｅｒｔｆｏｒＰｎｖｅｌｏｃｉｔｙｖａｒｉａｔｉｏｎａｎｄａｎｉｓｏｔｒｏｐｙｏｆｔｈｅｕｐｐｅｒｍｏｓｔｍａｎｔｌｅｂｅｎｅａｔｈｔｈｅｍｉｄｄｌｅｓｏｕｔｈｅｒｎｓｅｇｍｅｎｔｏｆｔｈｅＴａｎＬｕｆａｕｌｔａｎｄｉｔｓａｄｊａｃｅｎｔａｒｅａｓｂｙＰｎｗａｖｅｔｉｍｅｔｅｒｍｔｏｍｏｇｒａｐｈｙｍｅｔｈｏｄ．ＴｈｅｒｅｓｕｌｔｉｎｄｉｃａｔｅｓｔｈａｔｏｂｖｉｏｕｓｌａｔｅｒａｌｈｅｔｅｒｏｇｅｎｅｉｔｉｅｓｅｘｉｓｔｉｎｔｈｅｕｐｐｅｒｍｏｓｔｍａｎｔｌｅｂｅｎｅａｔｈｔｈｅｓｔｕｄｙａｒｅａａｎｄｔｈｅＰｎｖｅｌｏｃｉｔｙｖａｒｉｅｓｆｒｏｍ７．６８～８．２４ｋｍ·ｓ－１ｃｏｍｐａｒｅｄｔｏｔｈｅａｖｅｒａｇｅｖｅｌｏｃｉｔｙｏｆ７．９５ｋｍ·ｓ－１．ＴｈｅＰｎｖｅｌｏｃｉｔｙｓｔｒｕｃｔｕｒｅａｌｏｎｇＴａｎＬｕｆａｕｌｔｈａｓｓｅｇｍｅｎｔａｔｉｏｎ．ＴｈｅｍｉｄｄｌｅｓｅｇｍｅｎｔｏｆｔｈｅＴａｎＬｕｆａｕｌｔｚｏｎｅｓｈｏｗｓｏｂｖｉｏｕｓＮＥｔｒｅｎｄｉｎｇｌｏｗｖｅｌｏｃｉｔｙｚｏｎｅ．Ｌｏｗｖｅｌｏｃｉｔｙｍａｙｒｅｓｕｌｔｆｒｏｍｔｈｅｓｔｒｏｎｇｌｙｔｈｉｎｎｉｎｇｏｆｔｈｅｌｉｔｈｏｓｐｈｅｒｅａｎｄｔｈｅｕｐｗｅｌｉｎｇｏｆｔｈｅａｓｔｈｅｎｏｓｐｈｅｒｅｈｏｔｍａｔｅｒｉａｌ．ＴｈｅｓｏｕｔｈｅｒｎｐａｒｔｏｆｔｈｅＴａｎＬｕｆａｕｌｔｚｏｎｅｈａｓＮＮＥｔｒｅｎｄｉｎｇｗｅａｋｈｉｇｈｖｅｌｏｃｉｔｙａｎｏｍａｌｙｚｏｎｅ，ｗｈｉｃｈｉｓｃｏｎｓｉｄｅｒｅｄａｓａｈｉｇｈｌｏｗｖｅｌｏｃｉｔｙｂｏｕｎｄａｒｙｏｆｔｗｏｄｉｆｆｅｒｅｎｔｔｅｃｔｏｎｉｃｂｌｏｃｋｓ，ＮｏｒｔｈＣｈｉｎａａｎｄＹａｎｇｔｚｅ．ＷｅｈｏｌｄｔｈａｔＴａｎＬｕｈａｓｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｔｅｃｔｏｎｉｃｂｌｏｃｋｂｏｕｎｄａｒｙ．ＴｈｅｂｏｕｎｄａｒｙｃｕｒｖｅｔｏｔｈｅＮｏｒｔｈＣｈｉｎａｆｒｏｍＳＥｔｏＮＷｉｎＪｉａｎｇｓｕｄｏｍａｉｎ．ＴｈｅｉｎｔｅｎｓｉｔｙｄｉｓｔｒｉｂｕｔｉｏｎｏｆＰｎｖｅｌｏｃｉｔｙａｎｉｓｏｔｒｏｐｙｉｓｃｏｒｒｅｌａｔｅｄｗｉｔｈｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆｖｅｌｏｃｉｔｙａｎｏｍａｌｙ．Ｏｎｔｈｅｗｈｏｌｅ，ｖｅｌｏｃｉｔｙａｎｉｓｏｔｒｏｐｙｉｓｓｔｒｏｎｇｉｎｌｏｗｖｅｌｏｃｉｔｙａｎｏｍａｌｉｅｓｌｉｋｅＬｕｘｉｕｐｌｉｆｔａｎｄｓｏｕｔｈｅｒｎｒｅｇｉｏｎｏｒｔｈｅｔｒａｎｓｉｔｉｏｎｚｏｎｅｆｒｏｍｈｉｇｈｔｏｌｏｗｖｅｌｏｃｉｔｉｅｓｎｅａｒｔｈｅＭａｏｓｈａｎｆａｕｌｔｚｏｎｅ．ＡｎｉｓｏｔｒｏｐｙｉｓｒｅｌａｔｉｖｅｌｙｗｅａｋｕｎｄｅｒｓｔａｂｌｅｂｌｏｃｋａｎｄｂａｓｉｎｌｉｋｅＳｕｂｅｉｂａｓｉｎａｎｄＨｅｆｅｉｂａｓｉｎ，ｗｈｉｃｈｈａｓｈｉｇｈｖｅｌｏｃｉｔｙａｎｏｍａｌｉｅｓ．ＰｎｖｅｌｏｃｉｔｙａｎｉｓｏｔｒｏｐｙｉｓｗｅａｋａｌｏｎｇＴａｎＬｕｆａｕｌｔｚｏｎｅａｎｄｃａｎｎｏｔｂｅｄｅｔｅｃｔｅｄｍａｉｎｌｙｂｙＰｎｗａｖｅｅｓｐｅｃｉａｌｙｉｎｔｈｅｍｉｄｄｌｅｓｅｇｍｅｎｔ．Ｓｔｒｏｎｇｔｅｃｔｏｎｉｃｍｏｖｅｍｅｎｔｒｅｓｕｌｔｅｄｆｒｏｍｔｈｅｕｐｗｅｌｉｎｇｏｆａｓｔｈｅｎｏｓｐｈｅｒｅｈｏｔｍａｔｅｒｉａｌａｌｏｎｇｔｈｅｗｅａｋｔｅｃｔｏｎｉｃｂｅｌｔ（ＴａｎＬｕ）ｆａｕｌｔｗｅａｋｅｎｅｄｔｈｅａｎｉｓｏｔｒｏｐｙｔｒａｃｅ，ｗｈｉｃｈｌｅｆｔｆｒｏｍｅａｒｌｙｔｅｃｔｏｎｉｃｍｏｖｅｍｅｎｔｉｎｔｈｅｌｏｗｅｒｐａｒｔｏｆｔｈｅｌｉｔｈｏｓｐｈｅｒｅ．ＴｈｅｆａｓｔｗａｖｅｄｉｒｅｃｔｉｏｎｎｅａｒｓｏｕｔｈｅｒｎｓｅｇｍｅｎｔｏｆＴａｎｌｕｉｓｃｏｎｓｉｓｔｅｎｔｗｉｔｈｔｈｅｆｒａｃｔｕｒｅｅｘｔｅｎｓｉｏｎｄｉｒｅｃｔｉｏｎ．ＴｈｅＰｎｖｅｌｏｃｉｔｙｌａｔｅｒａｌｖａｒｉａｔｉｏｎｈａｓａｃｅｒｔａｉｎｃｏｒｒｅｌａｔｉｏｎｗｉｔｈｓｔｒｏｎｇｅａｒｔｈｑｕａｋｅａｃｔｉｖｉｔｙ．Ｍｏｓｔｏｆｍａｊｏｒｅａｒｔｈｑｕａｋｅｓｏｃｃｕｒｅｄｏｎｔｈｅｅｄｇｅｓｏｆｌｏｗｖｅｌｏｃｉｔｙａｎｏｍａｌｉｅｓｏｒｔｈｅｔｒａｎｓｉｔｉｏｎｚｏｎｅｆｒｏｍｈｉｇｈｔｏｌｏｗｖｅｌｏｃｉｔｉｅｓ．ＴｈｅｅｐｉｃｅｎｔｅｒｏｆｔｈｅＴａｎｃｈｅｎｇ８．５ｅａｒｔｈｑｕａｋｅｌｏｃａｔｅｏｎｔｈｅｔｒａｎｓｉｔｉｏｎｚｏｎｅｆｒｏｍｈｉｇｈｔｏｌｏｗｖｅｌｏｃｉｔｉｅｓｏｆｔｈｅｍｉｄｄｌｅｓｅｇｍｅｎｔａｎｄｓｏｕｔｈｅｒｎｓｅｇｍｅｎｔ，ｗｈｅｒｅｔｈｅｌａｔｅｒａｌｖａｒｉａｔｉｏｎｏｆｖｅｌｏｃｉｔｙｉｓｂｉｇｇｅｓｔａｎｄｓｔｒｅｓｓｍｏｓｔｃｏｎｃｅｎｔｒａｔｅｅａｓｉｌｙａｎｄｐｒｏｄｕｃｅｓｔｒｅｓｓｄｉｆｆｅｒｅｎｃｅ．犓犲狔狑狅狉犱狊　ＴａｎＬｕｆａｕｌｔ；Ｍｉｄｄｌｅｓｏｕｔｈｅｒｎｓｅｇｍｅｎｔ；Ｐｎｗａｖｅｖｅｌｏｃｉｔｙ；Ａｎｉｓｏｔｒｏｐｙ；Ｔａｎｃｈｅｎｇｅａｒｔｈｑｕａｋｅ１　引言郯庐断裂是一条纵贯我国大陆东部、ＮＮＥ走向的巨型深断裂带，对东部现代地震活动、区域构造和矿产资源的形成和分布都起着非常重要的控制作用．根据地震活动性、构造特点及演化历史等特征可将其划分为北、中及南三段，分段的结点为渤海湾和江苏新沂（王小凤等，２０００）．郯庐断裂带中南段及邻区（１１５°Ｅ—１２２°Ｅ，２９°Ｎ—３８°Ｎ）按照前人有关地质构造单元的划分（郑剑东等，１９８９），自南往北跨越了华南褶皱系、扬子断块区及华北断块区这３个一级地质构造单元．此外，该区域还包含了被地学家公认为全世界规模最大的苏鲁大别超高压变质带（Ｘｕｅｔａｌ．，１９９８），整个研究区域的地质构造条件和演化非常复杂．由于该地区的复杂性和重要性，长期以来吸引了许多地球物理学家围绕郯庐断裂带的形成、演化和活动机制等问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 展开研究，且获得了卓有成效的研究成果．人工地震测深结果显示（马杏垣等，１９９１；陈沪生等，１９９３；赵志新和徐纪人，２００９；刘保金等，２０１５），郯庐断裂带为一个切割莫霍面、深达岩石圈５０５地球物理学报（ＣｈｉｎｅｓｅＪ．Ｇｅｏｐｈｙｓ．）５９卷　并控制地幔高温热物质上涌及岩浆侵入的深大断裂带．大地电磁测深结果表明（肖骑彬等，２００８；张继红等，２０１０），郯庐断裂带深部为穿透Ｍｏｈｏ面，向下延伸至上地幔的走滑构造．以上研究表明郯庐断裂带中南段是切割整个地壳并深入上地幔的断裂破碎带和构造薄弱带，是上地幔高温高压热物质上涌的通道．天然地震的区域三维速度反演结果表明（徐佩芬等，２０００；Ｃｈｅｎｅｔａｌ．，２００６；黄忠贤等，２００９；黄耘等，２０１１；范小平等，２０１５），郯庐断裂带中南段速度结构呈现分段性，地壳内地震波速度跨郯庐两侧出现明显的变化，南段两侧Ｍｏｈｏ面分界线清晰．重力研究表明（唐新功等，２００６），郯庐断裂带中南段东侧处于布格重力高值异常区，西侧处于低值异常区；磁法结果揭示（李春峰等，２００９），郯庐带构成东侧负磁异常区和西侧正磁异常区的边界．由此可见，郯庐带作为高低值过渡带清晰地勾勒出了不同块体的边界，该边界已切穿地壳直抵上地幔．这些研究成果虽然使用了不同的资料和 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 ，但从各个角度揭示了郯庐断裂带及周边地区不同埋藏深度的特征，为研究郯庐断裂带及邻区的演化、深部结构特征提供了丰富的地震学依据，但是没有专门针对郯庐断裂带上地幔顶部速度结构或者各向异性进行研究．Ｍｏｈｏ面下方暨上地幔顶部是地壳和岩石圈地幔物质交换和能量传递的关键区域．Ｐｎ震相是地震波沿着莫霍界面滑行产生的折波，在上地幔顶部比其他震相（李大虎等，２０１５）能够得到更加密集的射线覆盖，从而分辨率和精度优于其他震相．因此，它成为研究上地幔顶部Ｐ波速度结构和波速方位各向异性的首选震相．近年来，一些学者使用Ｐｎ波震相研究了我国上地幔顶部的速度结构（丁志峰等，１９９２；汪素云等，２００１，２００３；李娟，２００３；裴顺平等，２００４；Ｌｉａｎｇｅｔａｌ．，２００４；Ｐｅｉｅｔａｌ．，２００７；胥颐等，２００８；李志伟等，２０１１；李飞等，２０１１），结果表明波速存在明显的横向变化，并存在随空间变化的波速各向异性．但是，这些成果或者因为研究尺度较大导致分辨尺度不够；或者因为数据分布范围局限，不足以对郯庐断裂带中南段及郯城地震的深部构造细节进行研究．王良书等（１９９０）利用的Ｐｎ震相对本文研究区有所涉及，但其采用的参考台法应用范围比较小，且不能提供二维平面结果图，结果局限于参与计算的台站下方．截至目前，有关郯庐断裂带中、南段的分段结点仍存在争议，仍然未能见到专门针对郯庐断裂带中南段及邻区上地幔顶部体波速度结构和波速方位各向异性的结果，本文采用Ｐｎ波时间项层析成像法试获得上述问题的一些地震学依据．基于上述原因，本文利用国际地震中心（ＩＳＣ）、中国地震台网和区域地震台网观测积累的Ｐｎ波资料，尤其注重补充２００８年以来记录的高质量的Ｐｎ资料，采用Ｈｅａｒｎ和Ｎｉ（１９９４）提出的Ｐｎ波各向异性的时间项层析成像方法，重建了郯庐断裂带中南段及邻区上地幔顶部１°×１°分辨尺度的体波速度结构和２°×２°分辨尺度的方位各向异性．研究结果为郯庐断裂带中南段的分段特征、作为不同块体的分界及郯城８．５级地震的深部动力成因提供了地震学依据；郯庐断裂带中南段及邻区上地幔顶部结构的研究对于研究扬子地块与华北地块的陆陆碰撞过程及其两者的边界、超高压变质作用及大陆动力学等方面具有重要意义．２　资料和方法为使得本文研究的目标区能被更多的地震射线有效覆盖，本文扩大资料搜集范围为１１４°Ｅ—１２４°Ｅ、２８°Ｎ—３９°Ｎ，该区自北向南跨越了华北断块、胶辽断块、苏鲁断块、下扬子断块、大别褶皱带及华南褶皱系．研究区内的主要地质构造和断裂见图１，图件采用ＧＭＴ软件完成（ＷｅｓｓｅｌａｎｄＳｍｉｔｈ，１９９８），底图为研究区的地形地貌图．本文搜集了ＩＳＣ（１９６４—２０１４年）的地震观测报告、中国地震台网（１９８６—２０１４年）及江苏（１９８４—２０１４年）、山东（２００８—２０１４年）、安徽（１９７６—２０１４年）、浙江（１９８８—２０１４年）、江西（１９８４—２０１４年）、湖北（１９８０—２０１４年）、河南（１９８１—２０１４年）、河北（２００８—２０１４年）、上海（１９８３—２０１４年）九个省级地震台网的地震观测报告，从中严格挑选出符合条件的Ｐｎ波震相记录，为获得高质量的成像结果奠定了坚实的基础．挑选Ｐｎ波震相的原则是：（１）震中距大于１．２°，以保证最先到达的是Ｐｎ波；（２）震源深度小于３３ｋｍ；（３）每个地震能被３个以上台站同时观测获得Ｐｎ震相；（４）每个台站能记录到３个以上的Ｐｎ波震相；（５）删掉丛集地震（Ｌｉａｎｇｅｔａｌ．，２００４），分组距离定为４ｋｍ；（６）相对于７．９５ｋｍ·ｓ－１的Ｐｎ波平均速度，折合走时残差在３～９ｓ之间（图２ａ）．为了检查数据的准确性，做出Ｐｎ震相的震中距走时时距关系图２ｂ．为了减少将Ｐｇ震相错判为Ｐｎ震相的概率，将走时残差控制在±２．５ｓ以内（图２ｃ）．该走时残差通过以下模型获得：地壳厚度３３ｋｍ，地壳６０５　２期顾勤平等：郯庐断裂带中南段及邻区Ｐｎ波速度结构与各向异性图１　研究区的地形地貌特征及主要构造①郯城—庐江断裂，②江山—绍兴断裂，③淮阴—响水口断裂，④烟台—五莲断裂，⑤茅山断裂，⑥广济—襄樊断裂，⑦肥中断裂，⑧寿县—定远断裂．Ｆｉｇ．１　Ｔｈｅｔｏｐｏｇｒａｐｈｉｃａｎｄｇｅｏｍｏｒｐｈｉｃｆｅａｔｕｒｅｓａｎｄｍａｉｎｔｅｃｔｏｎｉｃｂｌｏｃｋｓｉｎｔｈｅｓｔｕｄｙａｒｅａｓ①ＴａｎｃｈｅｎｇＬｕｊｉａｎｇｆａｕｌｔ，②ＪｉａｎｇｓｈａｎＳｈａｏｘｉｎｇｆａｕｌｔ，③ＨｕａｉｙｉｎＸｉａｎｇｓｈｕｉｋｏｕｆａｕｌｔ，④ＹａｎｔａｉＷｕｌｉａｎｆａｕｌｔ，⑤Ｍａｏｓｈａｎｆａｕｌｔ，⑥ＧｕａｎｇｊｉＸｉａｎｇｆａｎｆａｕｌｔ，⑦Ｆｅｉｚｈｏｎｇｆａｕｌｔ，⑧ＳｈｏｕｘｉａｎＤｉｎｇｙｕａｎｆａｕｌｔ．平均速度为６．２ｋｍ·ｓ－１，用迭代最小二乘法对Ｐｎ震相走时随震中距的变化进行线性拟合，拟合直线斜率的倒数为初始速度模型，最终得到本文目标区上地幔顶部的平均速度为７．９５ｋｍ·ｓ－１，该值略低于全国水平，与前人获得该区域的结果一致（汪素云等，２００３；胥颐等，２００８）．按照上述挑选Ｐｎ震相的原则，共精心挑选出符合条件的４４３个地震事件和２０８个台站，合计６３８１条地震射线满足本文研究的条件．图３给出了符合条件的地震震中、台站以及射线分布图，由图可见，总体上本文目标区即郯庐断裂带中南段及邻区（１１５°Ｅ—１２２°Ｅ，２９°Ｎ—３８°Ｎ）有较好的交叉覆盖，而在研究区边缘和海域，受地震或台站的分布影响，射线覆盖则不太理想．将研究目标区上地幔顶部沿经纬度水平方向划分成狀个１５′×１５′同等大小的网格，采用Ｈｅａｒｎ反演Ｐｎ波速度和各向异性的方法，反演计算Ｐｎ速度结构的同时将波速方位各向异性考虑在内，即考虑不同地震射线通过某个网格时，射线穿过的方位角不同速度不同，将Ｐｎ走时残差表示如下：狋犻犼＝犪犻＋犫犼＋∑犱犻犼犽（狊犽＋犃犽ｃｏｓ２＋犅犽ｓｉｎ２），式中狋犻犼为台站观测走时与理论计算走时之差；犪犻为第犻个台站的走时延迟扰动；犫犼为第犼个地震事件的走时延迟扰动；犱犻犼犽为第犻个台站接收第犼个地震时在通过第犽个网格内的距离；狊犽为第犽个网格的慢度扰动；犃犽、犅犽为第犽个速度网格内的各向异性系数；为地震震中相对接收台站的反方位角．Ｐｎ波速各向异性的大小用（犃２犽＋犅２犽）１／２表示，快波方向为１／２ａｒｃｔａｎ（犅犽／犃犽）＋９０°．震相数据的随机误差和地震射线的不均匀分布导致了所求方程具有病态特征和奇异值产生跳跃，为使得反演结果稳定我们加入了平滑约束，采用带阻尼的ＬＳＱＲ法求解以上方程中的犃犽、犅犽、犪犻、犫犼和狊犽这五个未知参数．我们对迭代次数、速度扰动及各向异性系数的阻尼系数做了多次对比试算，最终确定迭代次数为１２０次，速度扰动及各向异性系数的阻尼都取为３００．３　反演结果利用Ｈｅａｒｎ和Ｎｉ（１９９４）提出的Ｐｎ波各向异性的时间项层析成像方法，反演得到了郯庐断裂带中南段及邻区上地幔顶部相对于平均速度７．９５ｋｍ·ｓ－１的Ｐｎ波速度横向变化及波速方位各向异性分布图像．图２ｃ和２ｄ分别为反演前、后走时残差分布的统计结果图，走时残差的标准差由反演前的１．２１ｓ降低为０．７５ｓ，经过反演走时残差得到了一定程度的收敛．３．１　分辨率为验证反演结果的稳定性和可靠性，本文采用检测板（ｃｈｅｃｋｂｏａｒｄ）方法对数据体的分辨能力进行了测验．检测板测试就是给定初始速度模型，加上一定的速度扰动作为观测走时，再采用同样的控制参数和反演方法进行反演计算，看反演的结果能否恢复事先给定的初始速度模型（裴顺平，２００２）．限于篇幅，这里仅给出网格大小为１．５°×１．５°、１．０°×１．０°的速度异常恢复图像（图４）和２°×２°的各向异性恢复图像（图５）．反演结果的分辨率和可靠性依赖于数据的数量及质量，本文特别注重补充研究范围内２００８年以来记录的大量高质量的Ｐｎ波震相资料，由于数据的充实，本次获得的速度结构所达到的１°×１°分辨率比前人（汪素云等，２００３；Ｐｅｉｅｔａｌ．，２００７；胥颐等，２００８）在该地区采用Ｐｎ波速度反演所获得的分辨率有所提高．各向异性分辨尺度相对速度结构的分辨尺度而言，由于反演计算时增加了７０５地球物理学报（ＣｈｉｎｅｓｅＪ．Ｇｅｏｐｈｙｓ．）５９卷　图２　资料走时曲线及残差分布情况（ａ）相对７．９５ｋｍ·ｓ－１的平均速度得到的初始走时残差分布；（ｂ）经挑选后的数据体绘制的走时曲线；（ｃ）反演前基于初始模型的走时残差分布；（ｄ）反演后走时残差分布．Ｆｉｇ．２　Ｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｒａｖｅｌｔｉｍｅｃｕｒｖｅｓａｎｄｒｅｓｉｄｕａｌｓ（ａ）Ｉｎｉｔｉａｌｔｒａｖｅｌｔｉｍｅｒｅｓｉｄｕａｌｓｒｅｌａｔｉｖｅｔｏａｖｅｒａｇｅｖｅｌｏｃｉｔｙｏｆ７．９５ｋｍ·ｓ－１；（ｂ）Ｔｒａｖｅｌｔｉｍｅｃｕｒｖｅｓｏｆｃｈｏｓｅｎｄａｔａ；（ｃ）Ｔｒａｖｅｌｔｉｍｅｒｅｓｉｄｕａｌｓｄｉｓｔｒｉｂｕｔｉｏｎｂａｓｅｄｏｎｔｈｅｉｎｉｔｉａｌｖｅｌｏｃｉｔｙｍｏｄｅｌｂｅｆｏｒｅｉｎｖｅｒｓｉｏｎ；（ｄ）Ｔｉｍｅｒｅｓｉｄｕａｌｓａｆｔｅｒｉｎｖｅｒｓｉｏｎ．图３　Ｐｎ波射线路径及台站（白色三角）、地震（白色小圆圈）分布Ｆｉｇ．３　Ｐｎｗａｖｅｒａｙｐａｔｈｓａｎｄｄｉｓｔｒｉｂｕｔｉｏｎｏｆｓｔａｔｉｏｎｓ（ｗｈｉｔｅｔｒｉａｎｇｌｅｓ）ａｎｄｅａｒｔｈｑｕａｋｅｓ（ｗｈｉｔｅｃｉｒｃｌｅｓ）２个各向异性未知量即犃犽、犅犽的求解，在同等射线数的条件下分辨尺度不如速度结构，但由图５可见，总体上２°×２°分辨较好．由图４、５可以看出，研究范围内大部分地区的速度和波速各向异性能够较好恢复，南黄海由于台站布设限制、射线角度单一，导致分辨率有所下降．这也说明了可靠、稳定的反演结果须有较密集的射线覆盖和合理的射线分布．整体上，本文获得的反演结果在主要地区的速度结构及各向异性分布形态与已有的Ｐｎ波层析成像结果或者其他层析方法获得的结果总体一致（Ｗａｎｇｅｔａｌ．，２００３；胥颐等，２００８；赵志新和徐纪人，２００９），但本文给出了郯庐断裂带中南段及邻区上地幔顶部体波速度结构和波速各向异性的更多细节，可以为 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 研究目标区的深部结构特点提供更多新的信息．３．２　上地幔顶部犘狀波速度结构物质的组成成分、含水量以及压强等的变化都８０５　２期顾勤平等：郯庐断裂带中南段及邻区Ｐｎ波速度结构与各向异性图４　Ｐｎ波速度的检测板测试结果：（ａ）１．５°×１．５°和（ｂ）１．０°×１．０°Ｆｉｇ．４　ＣｈｅｃｋｂｏａｒｄｔｅｓｔｓｆｏｒＰｎｖｅｌｏｃｉｔｙｗｉｔｈｒｅｓｏｌｕｔｉｏｎｓｏｆ（ａ）１．５°×１．５°ａｎｄ（ｂ）１．０°×１．０°图５　Ｐｎ波各向异性２°×２°的检测板测试结果Ｆｉｇ．５　ＣｈｅｃｋｂｏａｒｄｔｅｓｔｓｆｏｒＰｎａｎｉｓｏｔｒｏｐｙｗｉｔｈｒｅｓｏｌｕｔｉｏｎｓｏｆ２°×２°能够引起Ｐｎ波速度的变化，但温度对其影响较大（Ｈｅａｒｎ，１９８４；ＢｅｇｈｏｕｌａｎｄＢａｒａｚａｎｇｉ，１９９５；Ｈｅａｒｎｅｔａｌ．，２００４）．图６给出了利用Ｐｎ波时间项层析成像法反演得到的郯庐断裂带中南段及邻区上地幔顶部Ｐｎ波速度结构分布图像，相对于平均速度７．９５ｋｍ·ｓ－１而言，其值在７．６８～８．２４ｋｍ·ｓ－１范围内变化．图中红色区域表示体波速度小于平均速度，蓝色区域表图６　Ｐｎ波速度横向变化黑色圆圈标识的是１９００年以来大于６级的地震事件，黑色圆点标识地名，红色五角星代表１６６８年郯城８．５级地震．Ｆｉｇ．６　ＰｎｖｅｌｏｃｉｔｙｌａｔｅｒａｌｖａｒｉａｔｉｏｎｓＴｈｅｂｌａｃｋｃｉｒｃｌｅｓｄｅｎｏｔｅｌａｒｇｅｅａｒｔｈｑｕａｋｅｗｉｔｈｍａｇｎｉｔｕｄｅｓｅｑｕａｌｔｏｏｒｇｒｅａｔｅｒｔｈａｎ犕６．０ｔｈａｔｏｃｃｕｒｒｅｄｓｉｎｃｅ１９００，ｂｌａｃｋｃｉｒｃｌｅｓａｎｄｄｏｔｓｄｅｎｏｔｅｈｏｍｏｎｙｍ，ｒｅｄｓｔａｒｓｄｅｎｏｔｅｔｈｅ１６６８Ｔａｎｃｈｅｎｇ（８．５）ｅａｒｔｈｑｕａｋｅ．示大于平均速度．由图可以看出，研究范围内速度结构表现出明显的横向非均匀性．总体而言，波速分布９０５地球物理学报（ＣｈｉｎｅｓｅＪ．Ｇｅｏｐｈｙｓ．）５９卷　形态不仅与前人采用相同方法（汪素云等，２００３；Ｐｅｉｅｔａｌ．，２００７；胥颐等，２００８）所得结果一致性较好，而且得到了由广角反射法探测获得的中国东部地壳三维Ｐ波在同等深度上结果的支持（赵志新和徐纪人，２００９）．广角反射分辨率和精度都优于天然地震，其所获得的结果具有较高的可信度，两者较好的契合关系表明了笔者所用的资料及方法反演获得的结果具有稳定性和可靠性．郯庐断裂带两侧的速度异常和分布形态明显不同，断裂西侧的华北地块总体速度偏低，速度异常为北西方向展布，而断裂东侧的扬子地块速度较高，异常为北东或近东西向．Ｐｎ波速度沿郯庐断裂带表现出分段特征．沿中段及周边存在一平行于郯庐断裂带的明显ＮＥ向低速异常带，其值约为７．７～７．９ｋｍ·ｓ－１，极低速异常位于１１９．０°Ｅ，３５．５°Ｎ．跨郯庐带中段布设的响水—满都拉地震剖面揭示郯庐带是曾伴随有深部物质上侵的巨大构造边界，中段切割Ｍｏｈｏ面、Ｍｏｈｏ面明显上隆但断裂两侧莫霍面起伏却不大（马杏垣等，１９９１）．接收函数和面波结果显示（Ｃｈｅｎｅｔａｌ．，２００６；黄忠贤等，２００９），中段软流圈顶面及上地幔顶部均显著隆起，而断裂两侧莫霍面基本处于同一水平面．本文沿中段较低的Ｐｎ波速度进一步证实，该段剪切作用已深达上地幔顶部．郯庐断裂带中段Ｍｏｈｏ面附近电性结构也明显，存在高温、高导体（肖骑彬等，２００８；张继红等，２０１０），成为高温高压热物质上涌通道．郯庐断裂带热流值在５３．９～８１．６ｍＷ·ｍ－２之间变动（马杏垣等，１９９１；胡圣标等，２００１），表现为高热流特征，与其长期活动性有关．郯庐断裂带在磁异常图上（李春峰等，２００９）显示的线性串珠状暗示着地幔上穹与岩浆灌入．牛漫兰等（２００５）通过对郯庐带中南段新生代的幔源玄武岩喷发的研究，认为这些时期断裂带的活动影响到了岩石圈地幔．综上所述，导致中段具有低速分布特征最有可能是岩石圈减薄和软流圈高温物质沿郯庐带中段上涌至地壳，而断裂两侧莫霍面基本连续、横向差异不明显是造成断裂两侧Ｐｎ波速没有明显变化的原因．南段表现为ＮＮＥ向弱高波速异常带，清晰地勾勒出了华北地块与扬子地块两个不同块体的边界，该边界在重磁异常分布图上（张鹏等，２００７）及面波层析反演获得的Ｍｏｈｏ面埋深分布图上（黄忠贤等，２００９）十分明显．边界以西的华北地块Ｐｎ波速度明显偏低，而边界以东的下扬子地区Ｐｎ波速度相对偏高，它们之间的分界大体上对应于郯庐断裂带南段，分界线被上地幔顶部的Ｐｎ震相探测获得，表明华北块体和扬子块体之间的边界带已深切至上地幔顶部．郯庐断裂带南段高、低速分界线在江苏域向华北地块ＮＷ方向凹近，这一特征在Ｗａｎｇ等（２００３）、赵志新和徐纪人（２００９）采用广角反射法所获得的结果中也曾有所体现．这一凹进高速区（江苏新沂至泗洪段）也正是郯庐带中南段地震活动水平最低的区域．凹进现象从一定程度上揭示出了扬子块体向北西俯冲于华北地块之下．凹进弯曲的高速异常区与中段的低速区边界成为郯庐带Ｐｎ波速度横向变化最大的区域，在此容易产生应力集中和应力差，１６６８年郯城８．５级大地震就发生于该处．迄今为止，以上２个不同块体的分界特征也已由精度较高的深地震反射得以验证：中国东部灵璧—奉贤ＨＱ１３地学断面（陈沪生等，１９９３）揭示，郯庐断裂带南段两侧莫霍面起伏明显，自东往西有下倾加深的趋势，断面近于直立．刘保金等（２０１５）采用深地震反射探测方法跨南段布设了一条长为１００ｋｍ的地震剖面，结果表明Ｍｏｈｏ面和岩石圈底界面在郯庐带南段两侧分别表现出存在约０．５ｓ和１．７ｓ的整体双程走时差，断裂两侧垂直落差明显．精度较高的地震深反射剖面证实了断裂两侧Ｍｏｈｏ面性质不连续、横向差异明显，加之断裂两侧块体的不同历史演化及性质，我们认为是导致Ｐｎ波速度沿郯庐带南段表现出明显的高低速分界带的原因．李飞等（２０１１）曾经得出这样的结论：在以走滑性质为主的鲜水河断裂两侧Ｐｎ波速度没有显著变化；而在具有倾滑性质的龙门山逆冲断层两侧Ｐｎ波速度表现出强烈差异．综上所述，我们推测郯庐带南段新生代以来具有倾滑性质的逆冲特征，这一特征与郯庐带晚第三纪以来遭受强烈挤压的过程相符（王小凤等，２０００；刘保金等，２０１５）；而中段两侧莫霍面起伏不明显，表现出了显著的走滑构造特征．Ｐｎ波速度分布形态与地壳浅表地质构造单元具有一定的相关性．整个研究范围内最为明显的高速异常为苏北盆地、南黄海一带，达到８．１ｋｍ·ｓ－１以上．苏北盆地未发现大规模热流活动的迹象和地幔扰动（胡圣标，２００１；Ｗａｎｇｅｔａｌ．，２００３），结构基本完整且具有较强的刚性特征．位于其周边的合肥盆地、固镇盆地等也为明显的弱高波速特征，这说明盆地地区上地幔顶部温度低于周边造山带区域，内部构造活动较弱，表现出较大的应变强度．鲁西隆起及以南地区、大别—皖南地区及苏南隆起表现出明显的低速异常，在７．７５～７．９ｋｍ·ｓ－１左右，表明隆０１５　２期顾勤平等：郯庐断裂带中南段及邻区Ｐｎ波速度结构与各向异性起造山带岩石层地幔强度较弱、具有韧性，易于发生形变和流动．前人研究结果显示（汪素云等，２００３；Ｐｅｉｅｔａｌ．，２００７），郯庐带中段以西的鲁西地体上地幔顶部表现为低速异常，本研究在补充最新的高质量数据后，发现中段以西的鲁西地体并不表现为一个整体低速异常区，实际情况是一近南北向的弱高波速异常带夹杂在郯庐带中段及鲁西地体中间；前人结果显示郯庐带中段东侧高速，而本文显示为低速异常区，该低速异常在广角反射法（赵志新和徐纪人，２００９）获得的同等深度上的结果中也较为清楚．华北地块总体上表现为低速异常，低速可能与太平洋板块俯冲导致我国东部地壳的拉伸减薄和伴生的岩浆活动有关．华北ＮＷ向的低速异常区与弱高波速异常区相间排列，这一现象与北西或北北西向隆起、盆地相间排列的地形地貌基本一致．南黄海北部与胶东半岛都具有低Ｐｎ波速度异常且相连，两者岩石层地幔性质可能一致．南黄海其余地方Ｐｎ波速度较高，大于８．０ｋｍ·ｓ－１，不存在明显的低速异常分布，推测上地幔顶部不存在明显的热异常和大规模的地幔扰动．沿苏鲁造山带存在一个狭长的ＮＥ向低速带区域，该低速特征在广角反射下地壳３０ｋｍ深度处也存在（赵志新和徐纪人，２００９），下地壳有可能存在着地壳山根．３．３　犘狀波速度各向异性图７为研究区Ｐｎ波速各向异性方向分布图，将各向异性结果投影在速度结果图像上以便作对比讨论．由图可见，研究范围内Ｐｎ波速各向异性具有随空间变化的差异性，表明岩石圈地幔的构造形变较为复杂．Ｐｎ波各向异性通常被认为是由于上地幔顶部橄榄石晶格因变形导致其优势定向排列（ＬＰＯ）造成的（ＮｉｃｏｌａｓａｎｄＣｈｒｉｓｔｅｎｓｅｎ，１９８７）．Ｐｎ波速方位各向异性探测得到的是上地幔顶部最近一次构造运动被“冻结”的形变痕迹．获取上地幔各向异性的方法较多，如ＳＫＳ波分裂方法（吴萍萍等，２０１２）、面波的方位各向异性和偏振各向异性（顾勤平等，２０１０；彭艳菊等，２００７）、Ｐｎ波层析技术（Ｐｅｉｅｔａｌ．，２００７）等，本文使用Ｐｎ波层析技术获得了研究区上地幔顶部Ｐ波波速方位各向异性结果．在Ｐｎ波方位各向异性较为明显的区域，快波方向与地壳浅表运动方向（杨国华等，２００２）、最大主压应力方向一致性较好（许忠淮，２００１；徐纪人和赵志新，２００６），表明物质运移在地壳和上地幔顶部之间存在一定的继承性．Ｐｎ波速方位各向异性的强弱图７　Ｐｎ波速度各向异性图中黑色短线长度表示各向异性大小，短线方向表示各向异性快波方向．Ｆｉｇ．７　ＴｏｍｏｇｒａｐｈｉｃｉｍａｇｅｏｆＰｎｖｅｌｏｃｉｔｙａｎｉｓｏｔｒｏｐｙＢｌａｃｋｌｉｎｅｓｅｇｍｅｎｔｓａｒｅｄｒａｗｎｐａｒａｌｅｌｔｏｔｈｅｄｉｒｅｃｔｉｏｎｏｆｆａｓｔＰｎｖｅｌｏｃｉｔｙ，ｗｉｔｈｔｈｅｉｒｌｅｎｇｔｈｐｒｏｐｏｒｔｉｏｎａｌｔｏｔｈｅｍａｇｎｉｔｕｄｅｏｆａｎｉｓｏｔｒｏｐｙ．分布与速度结构分布相关联．总体上，Ｐｎ波低速异常区或者高低速异常过渡带，其方位各向异性较为强烈，如具有Ｐｎ波低速分布的鲁西地体、大别—皖南及苏南隆起一带，这些造山带地区韧性大，容易发生变形；在鲁西块体北侧、茅山断裂西侧等高低速过渡带地区Ｐｎ波速各向异性较强，这些地方容易应力集中和产生应力差．而在稳定地块和盆地下方各向异性较弱，如研究区内具有最为明显的高波速异常的苏北盆地，其各向异性强度为整个研究区最弱，几乎不能被Ｐｎ波探测获得，表明苏北盆地岩石层地幔的构造变形较弱，具有构造稳定地区的特点，与该地区长期以来较低的地震活动水平相吻合．Ｈｅａｒｎ（１９９６）指出，在简单剪切变形情况下，如板块边界和穿透地壳的剪切带以及岩石圈地幔的剪切拖曳带，各向异性的快波方向平行于剪切带；在纯剪切变形情况下，如岩石圈中没有断裂活动的地区，快波方向与最大拉伸方向一致．郯庐断裂带南段作为扬子断块与华北地块的边界带，在该处具有与断裂伸展方向近平行的ＮＥ向Ｐｎ波快波速方向，表明１１５地球物理学报（ＣｈｉｎｅｓｅＪ．Ｇｅｏｐｈｙｓ．）５９卷　南段以剪切变形为主，该地区的各向异性主要受南段对Ｍｏｈｏ面切割变形的影响．已有研究表明（姜刚仁等，１９８９），茅山断裂带的埋深也已深达岩石圈地幔，该处也存在着平行于断裂的Ｐｎ波快波方向．上地幔顶部快波方向与断裂伸展方向的一致说明这些地方地壳和上地幔存在一定的耦合关系．ＳＫＳ波分裂是上地幔各向异性的综合效应，以上两处Ｐｎ波快波方向也与ＳＫＳ波各向异性结果吻合（吴萍萍等，２０１２），表明这些地方地幔上涌引起的物质流动已运移至上地幔顶部．沿郯庐断裂带方位各向异性总体较弱，尤其是中段，基本无法被Ｐｎ震相探测得到，这可能是上地幔顶部被“冻结”下来的各向异性痕迹被软流圈热物质沿郯庐带上涌这一强烈构造运动给削弱了．扬子块体受郯庐带的剪切平移作用向北下插中朝块体，在南黄海海域形成多个ＮＥ—ＮＮＥ向伸展的凹陷、隆起及断裂分布，在该地区Ｐｎ波快波方向以ＮＥ向为主、规律性较强，两者的契合表明南黄海扬子块体向北侧中朝块体的下插作用已深达上地幔顶部．华北地块Ｐｎ波快波方向主要是ＮＷ—ＳＥ和Ｎ—Ｓ向，该结果除与Ｐｎ波层析法获得的快波方向一致外（Ｗａｎｇｅｔａｌ．，２００３），也得到地幔ＳＫＳ波分裂结果的支持（吴萍萍等，２０１２）．鉴于各向异性及研究区构造演化的复杂性，加之本次数据体和方法的片面性，在此不再对各向异性做进一步的讨论．３．４　犘狀波速度结构与强震之间的关系图６中黑色空心圆表示１９００年以来仪器测定的６级以上的强震震中分布，郯城８．５级地震因其震级大、影响范围广，在讨论有仪器记录以来的强震与Ｐｎ波速度结构的关系时也将其考虑在内．由图可见，Ｐｎ波速度的分布与强震的发生具有一定的相关性．强震大部分发生在上地幔顶部Ｐｎ波高、低速的过渡地带或者低速区域内，高低速的边界附近应力容易集中和产生应力差，低速区可以认为容易发生韧性形变而诱发地震，这些区域可能更利于地下能量和物质的传递．上述结果表明，强震的孕育及发生与上地幔顶部介质结构的横向非均匀性及介质的运动紧密相关．郯庐断裂带地震活动水平具有分段特征，地震活动总体呈现出中段强、南段弱的特征（张鹏等，２００７；刘保金等，２０１５），本文获得的Ｐｎ波速度在中段及南段具有截然不同的分布特征为这一认识提供了有利的地震学依据．沿郯庐断裂带中段表现出的低速异常带推测是由岩石圈地幔的高温物质沿构造薄弱带郯庐断裂上涌所形成，热物质的上涌导致地壳中应力的不均匀分布、降低了速度和刚度，它是该区域强震发生的深部动力来源，与该段地震活动水平高、强度大及历史强震主要发生在该段有着较好的契合关系．郯庐带南段成为高、低Ｐｎ波速度的分界，应力容易在此集中并产生应力差，我们认为自有仪器记录以来该段仍然处于应力集中阶段、未达到地震破裂程度，因此未曾记录到强震．郯庐断裂带江苏新沂至泗洪段现今地震活动微弱，仪器记录的小震也很少，李家灵等（１９９４）认为断层深部处于闭锁状态．本文获得的该段上地幔顶部Ｐｎ波高速特征为这一现象提供了地震学证据．郯城８．５级地震的震中（图６中红色五角星）位于中段和南段之间速度结构横向差异最大的区域，该区域的强横向非均匀结构导致构造应力最易于集中和最易于产生应力差，震中深部具有强烈的非均匀性和十分明显的深部构造特征．此外，２０１２年江苏高邮市、宝应县交界（１１９．６°Ｅ，３３．０°Ｎ）发生的犕Ｓ４．９破坏性地震，其震中位置位于两个高速区之间的相对低速部分，两个高速区之间相对低速部分同样也容易产生应力集中，本文结果说明此次破坏性地震同样存在着明显的深部结构特征．４　结论与讨论本文搜集了国际地震中心（ＩＳＣ）、中国地震台网和区域地震台网的Ｐｎ波到时资料，采用Ｈｅａｒｎ和Ｎｉ（１９９４）提出的Ｐｎ波时间项层析成像法，反演获得了郯庐断裂带中南段及邻区上地幔顶部Ｐｎ波速度和各向异性．反演结果揭示了研究区上地幔顶部速度结构、方位各向异性的横向变化及速度分布特征与强震活动之间的关系．（１）利用６３８１条Ｐｎ射线的走时资料，反演获得了郯庐断裂带中南段及邻区上地幔顶部Ｐｎ波速度结构及方位各向异性，结果表明研究区内Ｐｎ波平均速度为７．９５ｋｍ·ｓ－１，横向变化量从－０．２７～＋０．２９ｋｍ·ｓ－１，具有明显的空间差异性．Ｐｎ波速度分布形态沿郯庐断裂带具有分段性，中段及周边存在一明显的ＮＥ向低速异常带，推测构造薄弱带郯庐断裂带中段成为深部地幔高温物质上涌的通道，上涌形成的岩浆底侵或热侵蚀作用造成了岩石圈的减薄、降低了上地幔顶部的速度和刚度．郯庐断裂带南段表现为一ＮＮＥ向Ｐｎ波弱高波速异常带，清晰地勾勒出了华北地块与扬子地块２个不同块体２１５　２期顾勤平等：郯庐断裂带中南段及邻区Ｐｎ波速度结构与各向异性之间的分界，表明南段已切割至莫霍面、断裂两侧横向差异明显，具有明显的构造块体边界特征．南段作为华北与扬子２个不同块体的分界向ＮＷ方向凹进，这一凹进高速区（江苏新沂至泗洪段）也正是郯庐带中南段地震活动水平最低的区域．（２）速度结构分布形态与地壳浅表地质构造的地形面貌有着较好的契合关系．在地质构造稳定的盆地和地块，如苏北盆地、合肥盆地及固镇盆地等，盆地内部构造变形及岩浆活动较弱，Ｐｎ波速度偏高表明盆地上地幔顶部温度较低，是构造上稳定的地区；在盆地周边构造活动强烈的造山带则速度偏低，如鲁西隆起及以南、大别—皖南及苏南隆起一带．这说明了稳定块体与活动造山带的岩石层地幔特性不同．（３）在Ｐｎ波速方位各向异性明显的地方，Ｐｎ波快波方向与地壳浅表运动方向、最大主压应力方向在宏观上的基本一致，表明物质运移在地壳和上地幔之间是耦合的．Ｐｎ波速度各向异性强弱分布与速度异常分布相关联．整体上，如鲁西隆起及以南、大别—皖南等低速区或者茅山断裂带附近、鲁北等高低速过渡带，其速度各向异性较为强烈；而在具有高速或弱高波速的苏北盆地、合肥盆地等稳定地块和盆地下方各向异性相对弱些．郯庐断裂带附近方位各向异性较弱，尤其是中段，基本无法被Ｐｎ波探测获得，推测是上地幔顶部“冻结”下来的各向异性痕迹被软流圈热物质的上涌这一强烈构造运动削弱所致．南段具有与断裂伸展方向一致的ＮＥ向快波速方向，表明南段以剪切变形为主且已切割至莫霍面．（４）强震的空间分布与Ｐｎ波速度分布存在一定相关性，强震主要发生在上地幔顶部Ｐｎ波速度的低速异常区或高低波速异常过渡带，表明强震的孕育及发生可能与上地幔顶部介质的横向不均匀性及地幔物质的运移密切相关．南段高低速边界在江苏域向华北地块ＮＷ一侧凹进弯曲，郯城８．５级地震的震中位于郯庐带中段和南段这一速度横向变化最大的区域，也是构造应力最易于集中和应力差最易于产生的地区，表明郯城地震具有十分明显的深部结构特点．郯庐带在纬度３３°—３４．５°范围仪器记录地震少、现今地震活动微弱，本文获得的该区间上地幔顶部Ｐ波高波速异常特征为这一现象提供了地震学依据．致谢　本文图件采用ＧＭＴ绘制，地震数据来源于国际地震中心（ＩＳＣ）、中国地震局台网中心和区域地震台网，特此致谢！感谢汪素云老师在研究过程中的帮助，感谢梁春涛博士、李大虎博士、范小平博士、张鹏博士对本文有益的讨论．犚犲犳犲狉犲狀犮犲狊ＢｅｇｈｏｕｌＮ，ＢａｒａｚａｎｇｉＭ．１９９５．Ｃｏｍｍｅｎｔｏｎ“ｌａｔｅｒａｌｖａｒｉａｔｉｏｎｓａｎｄａｚｉｍｕｔｈａｌｉｓｏｔｒｏｐｙｏｆＰｎｖｅｌｏｃｉｔｉｅｓｂｅｎｅａｔｈＢａｓｉｎａｎｄＲａｎｇｅＰｒｏｖｉｎｃｅ”ｂｙＬｉａｎｓｈｅＺｈａｏ．犑．犌犲狅狆犺狔狊．犚犲狊．，１００（Ｂ７）：１２４５９１２４６１．ＣｈｅｎＨＳ，ＺｈｏｕＸＱ，ＬｉＤＧ，ｅｔａｌ．１９９３．ＧｅｏｓｃｉｅｎｃｅｓＴｒａｎｓｅｃｔｏｆＬｉｎｇｂｉＦｅｎｇｘｉａｎ（ＨＱ１３）ｉｎＥａｓｔｅｒｎＣｈｉｎａ（ｉｎＣｈｉｎｅｓｅ）．Ｂｅｉｊｉｎｇ：ＧｅｏｌｏｇｉｃａｌＰｕｂｌｉｓｈｉｎｇＨｏｕｓｅ．ＣｈｅｎＬ，ＺｈｅｎｇＴＹ，ＸｕＷＷ．２００６．ＡｔｈｉｎｎｅｄｌｉｔｈｏｓｐｈｅｒｉｃｉｍａｇｅｏｆｔｈｅＴａｎｌｕＦａｕｌｔｚｏｎｅ，ｅａｓｔｅｒｎＣｈｉｎａ：Ｃｏｎｓｔｒｕｃｔｅｄｆｒｏｍｗａｖｅｅｑｕａｔｉｏｎｂａｓｅｄｒｅｃｅｉｖｅｒｆｕｎｃｔｉｏｎｍｉｇｒａｔｉｏｎ．犑．犌犲狅狆犺狔狊．犚犲狊．，１１１（Ｂ９）：Ｂ０９３１２．ＤｉｎｇＺＦ，ＺｅｎｇＲＳ，ＷｕＤＭ．１９９２．ＴｈｅＰｎｗａｖｅｖｅｌｏｃｉｔｉｅｓａｎｄｔｈｅｒｅｌｉｅｆｏｆＭｏｈｏｉｎｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕ．犃犮狋犪犛犲犻狊犿狅犾狅犵犻犮犪犛犻狀犻犮犪（ｉｎＣｈｉｎｅｓｅ），１１（Ｓ１）：５９２５９９．ＦａｎＸＰ，ＬｉＸＢ，ＺｈａｏＱＧ，ｅｔａｌ．２０１５．Ｄｅｔｅｃｔｉｎｇｂｕｒｉｅｄｆａｕｌｔｓｔｒｕｃｔｕｒｅｓｂｙｍｉｃｒｏｔｒｅｍｏｒｓｕｒｖｅｙｍｅｔｈｏｄ．犃犮狋犪犛犲犻狊犿狅犾狅犵犻犮犪犛犻狀犻犮犪（ｉｎＣｈｉｎｅｓｅ），３７（１）：１３４１４３．ＧｕＱＰ，ＺｈｕＪＳ，ＫａｎｇＱＱ，ｅｔａｌ．２０１０．ＲａｙｌｅｉｇｈｗａｖｅｇｒｏｕｐｖｅｌｏｃｉｔｙａｎｄｉｔｓａｚｉｍｕｔｈａｌａｎｉｓｏｔｒｏｐｙｉｎｔｈｅＥｕｒａｓｉａａｎｄｗｅｓｔｅｒｎＰａｃｉｆｉｃｍａｒｇｉｎａｌｓｅａｒｅｇｉｏｎ．犃犮狋犪犛犲犻狊犿狅犾狅犵犻犮犪犛犻狀犻犮犪（ｉｎＣｈｉｎｅｓｅ），３２（１）：１２２２．ＨｅａｒｎＴＭ．１９８４．ＰｎｔｒａｖｅｌｔｉｍｅｓｉｎＳｏｕｔｈｅｒｎＣａｌｉｆｏｒｎｉａ．犑．犌犲狅狆犺狔狊．犚犲狊．，８９（Ｂ３）：１８４３１８５５．ＨｅａｒｎＴＭ，ＮｉＪＦ．１９９４．Ｐｎｖｅｌｏｃｉｔｉｅｓｂｅｎｅａｔｈｃｏｎｔｉｎｅｎｔａｌｃｏｌｉｓｉｏｎｚｏｎｅｓ：ｔｈｅＴｕｒｋｉｓｈＩｒａｎｉａｎＰｌａｔｅａｕ．犌犲狅狆犺狔狊．犑．犐狀狋．，１１７（２）：２７３２８３．ＨｅａｒｎＴＭ．１９９６．ＡｎｉｓｏｔｒｏｐｙＰｎｔｏｍｏｇｒａｐｈｙｉｎｔｈｅｗｅｓｔｅｒｎＵｎｉｔｅｄＳｔａｔｅｓ．犑．犌犲狅狆犺狔狊．犚犲狊．，１０１（Ｂ４）：８４０３８４１４．ＨｅａｒｎＴＭ，ＷａｎｇＳＹ，ＮｉＪＦ，ｅｔａｌ．２００４．ＵｐｐｅｒｍｏｓｔｍａｎｔｌｅｖｅｌｏｃｉｔｉｅｓｂｅｎｅａｔｈＣｈｉｎａａｎｄｓｕｒｒｏｕｎｄｉｎｇｒｅｇｉｏｎｓ．犑．犌犲狅狆犺狔狊．犚犲狊．，１０９（Ｂ１１）：Ｂ１１３０１，ｄｏｉ：１０．１０２９／２００３ＪＢ００２８７４．ＨｕＳＢ，ＨｅＬＪ，ＷａｎｇＪＹ．２００１．ＣｏｍｐｉｌａｔｉｏｎｏｆｈｅａｔｆｌｏｗｄａｔａｉｎｔｈｅＣｈｉｎａｃｏｎｔｉｎｅｎｔａｌａｒｅａ（３ｒｄｅｄｉｔｉｏｎ）．犆犺犻狀犲狊犲犑．犌犲狅狆犺狔狊．（ｉｎＣｈｉｎｅｓｅ），４４（５）：６１１６２６．ＨｕａｎｇＹ，ＬｉＱＨ，ＺｈａｎｇＹＳ，ｅｔａｌ．２０１１．ＣｒｕｓｔａｌｖｅｌｏｃｉｔｙｓｔｒｕｃｔｕｒｅｂｅｎｅａｔｈｔｈｅＳｈａｎｄｏｎｇＪｉａｎｇｓｕＡｎｈｕｉｓｅｇｍｅｎｔｏｆｔｈｅＴａｎｃｈｅｎｇＬｕｊｉａｎｇＦａｕｌｔＺｏｎｅａｎｄａｄｊａｃｅｎｔａｒｅａｓ．犆犺犻狀犲狊犲犑．犌犲狅狆犺狔狊．（ｉｎＣｈｉｎｅｓｅ），５４（１０）：２５４９２５５９，ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０００１５７３３．２０１１．１０．０１２．ＨｕａｎｇＺＸ，ＸｕＹ，ＨａｏＴＹ，ｅｔａｌ．２００９．ＳｕｒｆａｃｅｗａｖｅｔｏｍｏｇｒａｐｈｙｏｆｌｉｔｈｏｓｐｈｅｒｉｃｓｔｒｕｃｔｕｒｅｉｎｔｈｅｓｅａｓｏｆｅａｓｔＣｈｉｎａ．犆犺犻狀犲狊犲犑．犌犲狅狆犺狔狊．（ｉｎＣｈｉｎｅｓｅ），５２（３）：６５３６６２．ＪｉａｎｇＧＲ，ＺｈａｎｇＤＱ，ＺｈｕＱＬ，ｅｔａｌ．１９８９．ＳｅｉｓｍｏｔｅｃｔｏｎｉｃｓｏｆＪｉａｎｇｓｕＰｒｏｖｉｎｃｅ，ｉｎｃｌｕｄｉｎｇＳｈａｎｇｈａｉ．∥ＭａＸＹｅｄ．ＬｉｔｈｏｓｐｈｅｒｉｃＤｙｎａｍｉｃｓＡｔｌａｓｏｆＣｈｉｎａ（ｉｎＣｈｉｎｅｓｅ）．Ｂｅｉｊｉｎｇ：ＣｈｉｎａＣａｒｔｏｇｒａｐｈｉｃＰｕｂｌｉｓｈｉｎｇＨｏｕｓｅ．３１５地球物理学报（ＣｈｉｎｅｓｅＪ．Ｇｅｏｐｈｙｓ．）５９卷　ＬｉＣＦ，ＣｈｅｎＢ，ＺｈｏｕＺＹ．２００９．ＤｅｅｐｃｒｕｓｔａｌｓｔｒｕｃｔｕｒｅｓｏｆｅａｓｔｅｒｎＣｈｉｎａａｎｄａｄｊａｃｅｎｔｓｅａｓｒｅｖｅａｌｅｄｂｙｍａｇｎｅｔｉｃｄａｔａ．犛犮犻犲狀犮犲犻狀犆犺犻狀犪犛犲狉犻犲狊犇：犈犪狉狋犺犛犮犻犲狀犮犲狊，５２（７）：９８４