中国东北地区上地幔P波方位各向异性研究

地震体波走时层析成像是以地震体波到时作为观测数据,反演地下介质的地震波传播速度分布的一种反演技术方法.基于弱各向异性介质的假定,通过在走时反演方程中引入各向异性参数,可以同时进行介质的速度扰动和各向异性属性的联合反演.对三维介质的各向异性分析,有利于对区域结构构造演化、深部物质变形、动力学模式等问题的进一步解释.本文基...     

中国大陆及邻区海域地壳上地幔各向异性研究

利用分别由Love波和Rayleigh波得到的S波速度结构的差值（VSH-VSV）对中国大陆及邻区海域（70°E～145°E,15°N～55°N）地壳上地幔中的偏振各向异性进行研究.初步研究结果表明,各向异性在空间分布上存在不均匀性：（1）在小于150 km的深度范围内,VSH>VSV的各向异性体占主导地位,反映出在地球的浅部岩石圈内的水平应力作用及软流圈顶部物质的水平向流动对各向异性的形成起主导作用.在大陆地区,各向异性的强度随深度有显著变化.上地壳和上地幔盖层中的各向异性普遍较弱,而在流变性较强的下地壳和软流圈存在较大范围的各向异性.这一现象说明下地壳在岩石圈变形中可能有解耦作用.（2）在大于200 km深度的软流圈下部主要表现为VSHVSV的各向异性比较显著,而软流圈中VSHVSV各向异性.     

用地震面波研究上地幔各向异性

本文论述了上地幔各向异性的参数描述和主要证据,以及与地球动力学的关系。并就近几年来国外利用地震面波研究上地幔各向异性的重要成果和主要方法作了综合评述。     

中国东部及其邻区上地幔顶部Pn波速度结构及各向异性

中国东部地处欧亚板块东南部,紧邻西太平洋板块俯冲带,有着复杂的地质构造和深部结构,是国内外学者研究的重点区域.本文采用Pn波层析成像方法反演得到了中国东部及其邻区上地幔顶部Pn波速度结构及各向异性.研究中应用的Pn到时数据主要来自多种地震观测报告,并特别补充了东北流动台阵、华北流动台阵以及区域地震台网记录的地震事件,拾取了大量高精度的Pn到时数据,最终挑选出2049个台站记录的24072次地震的240814条Pn波到时数据.结果显示,中国东部地区上地幔顶部平均速度为8.06 km·s^-1,速度变化范围为7.81~8.32 km·s^-1.东北地区东部表现为显著的低速异常,华北克拉通中、东部Pn呈现低速,而西部地区则表现为高速异常,华南地块主体表现为高速.反演结果还揭示江汉盆地、下辽河盆地、二连盆地及海拉尔盆地都显示出高速,而四川盆地和松辽盆地内部则呈现出不均匀的结构特征.四川盆地的高速异常显示出明显的分块现象,这可能是该盆地在沉积前具有不同的基底物质;松辽盆地的北部呈现为高速,而南部却表现为低速异常,这一特征与松辽盆地南、北部分别为高、低热流相对应,暗示盆地南部的岩石圈已经历了改造.研究进一步揭示Pn波低速区和高、低速过渡带的各向异性也较为强烈,而大部分强震都发生在这些区域之上的地壳内,说明这些部位容易发生变形而应力集中或产生应力差.

     

中国大陆及邻区上地幔各向异性研究

为了深化对中国大陆及邻区上地幔介质形变强度和方向的空间变化图象、各向异性和剪切波分裂、各向异性以及应变、应变和构造过程关系的认识，我们以剪切波在分层各向异性介质中的传播理论为指导，以分裂剪切波观测为基础，借助数字化地震资料和高分辨率的信号处理技术，给出了对136次地震的分析结果，得到了20个台站下面的上地幔各向异性介质的变形强度和方向的空间变化图象；并结合已有研究结果，对中国大陆及邻区上地幔各向异性特性及其起源问题进行了讨论．      

华南上地幔P波速度结构

通过拟合15°-30°内的长周期P波走时及波形资料,得到了华南地区上地幔P波速度结构模型SC.研究结果表明,在华南地区的上地幔内无低速层存在,在405km和660km深度处仃一级间断面存在,速度跳跃分别为5.7％和4.6％.通过与欧洲西部上地幔模型K8,欧洲西北部上地幔模型S8和西藏地区上地幔模型QX8比较,我们发现各个构造区P波速度的过渡区是一致,但间断面的绝对深度不一样.     

华南上地幔P波速度结构

通过拟合15°—30°内的长周期P波走时及波形资料,得到了华南地区上地幔P波速度结构模型SC.研究结果表明,在华南地区的上地幔内无低速层存在,在405km和660km深度处仃一级间断面存在,速度跳跃分别为5.7％和4.6％.通过与欧洲西部上地幔模型K8,欧洲西北部上地幔模型S8和西藏地区上地幔模型QX8比较,我们发现各个构造区P波速度的过渡区是一致,但间断面的绝对深度不一样.     

四川及邻区上地幔各向异性研究

对布设在四川及邻区的固定和流动宽频带地震台网共94个台站记录作远震SKS波形资料偏振分析,采用最小切向能量的网格搜索和叠加分析求得每一个台站的SKS快波偏振方向和快、慢波的延迟时间,获得四川及邻区的上地幔各向异性图像．从得到结果整体看,各向异性快波方向在研究区东北部为NWW-SEE方向,到中部的NW-SE方向,再到西部的近NS方向,有顺时针旋转的趋势,主体以NW—SE方向为主,快、慢波延迟时间是0．47～1．68s．各向异性的快波偏振方向与GPS测量的地壳运动速度场方向变化相一致．通过分析各向异性图像,认为四川及邻区上地幔物质在区域构造应力场的作用下发生了变形,它使橄榄岩中晶格排列方向平行于物质变形方向,研究区壳幔变形可能存在垂直连贯变形特征．汶川地震可能是高原东缘地壳或上地壳物质受上地幔物质SE向的拖曳力作用向东挤压,这种挤压受到四川盆地下的刚性地块的顽强阻挡,经过长期的构造应力能量的积累,最终在龙门山映秀地区突然释放,导致汶川8．0级地震发生．     

中国云南地区上地幔地震各向异性

本文分析云南地区中周期ＳＫＳ偏振，得出观测台附近上地幔Ｓ波快速方向和分裂时间为：昆明，８５°，１．０秒；洱源，８９°，１．０秒；腾冲，１６５°，０．７秒；昭通１６７°，０．９秒。按上地幔地震各向异性是应力作用下矿物优势取向所致的观点，腾冲和昭通的各向异性观测与现代应力场作用相符。但昆明和洱源的观测与现代应力场作用矛盾，它可能是古生代至中生代间攀西裂谷张裂时期，地幔物质流动引起矿物有序排列的结果。     

川西龙门山及邻区地壳上地幔各向异性地震P波走时层析成像研究

经典的地震层析成像假定介质为各向同性,通过走时反演确定地球介质的非均匀结构,得到的仅是近似地球模型。如果各向异性较强时,应用基于各向同性假设的层析成像算法,可能导致错误的结果。非均匀与各向异性效应是耦合的,如果仅考虑非均匀效应则各向异性结构可能映射到非均匀结构中造成非均匀速度结构的误差,反之亦然。因此,高分辨率的走时层析成像必须同时考虑非均匀和各向异性两种效应。同时反演非均匀和各向异性结构(即各向异性层析成像),不仅可以获得可靠的各向同性速度结构,同时可以确定各向异性结构。相对于S波偏振分析,体波走时方法基于不同的数据集,能提供独立的地震各向异性观测证据。利用地震走时层析成像方法同时研究介质的非均匀性和各向异性,对于认识地球的结构及动力学过程都具有非常重要的意义。     

大别—苏鲁及邻区上地幔的各向异性

大别—苏鲁是扬子与华北的碰撞造山带,对该地区上地幔各向异性的研究有助于了解该区的地幔动力学机制.本文选用了中国数字化地震台网和区域数字地震台网(山东、安徽、江苏、河南、湖北)三分量宽频带的远震地震波形数据,分别采用最小能量法和旋转相关法,对大别—苏鲁及邻区进行剪切波偏振分析,计算了研究区台站下方介质的各向异性分裂参数:快波偏振方向(Φ)和快慢波延迟时间(δt).本文研究结果发现,研究区内快、慢波延迟时间0.5~1.63 s,推测各向异性层深度为57.5~187.6 km,由软流圈和岩石圈地幔的各向异性共同作用引起.快波偏振方向在4个不同构造区表现出不同的特点:华北板块快波偏振方向为近E-W向,根据地质资料,我们分析认为华北板块的各向异性受地幔软流圈流动的影响明显;大别造山带各向异性平行于大别主构造,反映造山过程中岩石圈物质沿大别造山轴部NW-SE向迁移的特点; 在大别南侧和东侧的扬子板块快波偏振方向分别表现为近垂直于造山带走向和NEE-SWW,苏鲁造山带各向异性结果为NEE-SWW,与地表构造有一定的夹角,同时与板块运动方向相差较大,分析认为扬子板块和苏鲁造山带各向异性是由地幔软流圈流动和印支—燕山期构造运动残留在岩石圈地幔的"化石各向异性"共同作用的结果.     

大别—苏鲁及邻区上地幔的各向异性

大别—苏鲁是扬子与华北的碰撞造山带,对该地区上地幔各向异性的研究有助于了解该区的地幔动力学机制.本文选用了中国数字化地震台网和区域数字地震台网(山东、安徽、江苏、河南、湖北)三分量宽频带的远震地震波形数据,分别采用最小能量法和旋转相关法,对大别—苏鲁及邻区进行剪切波偏振分析,计算了研究区台站下方介质的各向异性分裂参数:快波偏振方向(Φ)和快慢波延迟时间(δt).本文研究结果发现,研究区内快、慢波延迟时间0.5~1.63 s,推测各向异性层深度为57.5~187.6 km,由软流圈和岩石圈地幔的各向异性共同作用引起.快波偏振方向在4个不同构造区表现出不同的特点:华北板块快波偏振方向为近E-W向,根据地质资料,我们分析认为华北板块的各向异性受地幔软流圈流动的影响明显;大别造山带各向异性平行于大别主构造,反映造山过程中岩石圈物质沿大别造山轴部NW-SE向迁移的特点; 在大别南侧和东侧的扬子板块快波偏振方向分别表现为近垂直于造山带走向和NEE-SWW,苏鲁造山带各向异性结果为NEE-SWW,与地表构造有一定的夹角,同时与板块运动方向相差较大,分析认为扬子板块和苏鲁造山带各向异性是由地幔软流圈流动和印支—燕山期构造运动残留在岩石圈地幔的"化石各向异性"共同作用的结果.     

Azimuthal anisotropy of Rayleigh waves beneath the Tibetan Plateau and adjacent areas

The crustal and upper mantle azimuthal anisotropy of the Tibetan Plateau and adjacent areas was studied by Rayleigh wave tomography. We collected sufficient broadband digital seismograms trav-ersing the Tibetan Plateau and adjacent areas from available stations, including especially some data from the temporary stations newly deployed in Yunnan, eastern Tibet, and western Sichuan. They made an adequate path coverage in most regions to achieve a reasonable resolution for the inversion. The model resolution tests show that the anisotropic features of scope greater than 400 km and strength greater than 2% are reliable. The azimuthal anisotropy pattern inside the Tibetan Plateau was similar to the characteristic of tectonic partition. The crustal anisotropy strength is greater than 2% in most re-gions of East Tibet, and the anisotropy shows clockwise rotation surrounding the eastern Himalayan syntaxis. Vertically, the anisotropy direction indicates a coherent pattern within the upper crust, lower crust, and lithosphere mantle of the Tibetan Plateau, which also is consistent with GPS velocity field and SKS fast polarization directions. The result supports that the crust-mantle deformation beneath the Tibetan Plateau is vertically coherent. The anisotropy strength of crust and lithospheric upper mantle in Yunnan outside the Tibetan Plateau is lower than 2%, so SKS splitting from core-mantle boundary to station should largely be attributed to the anisotropy of asthenosphere.     

华北上地幔各向异性研究

对华北地震科学台阵的200个宽频带和甚宽带地震台站所记录的远震SKS(SKKS)波形资料作偏振分析,采用最小切向能量的网格搜索法和叠加分析方法求得每一个台站的SKS(SKKS)快波偏振方向和快、慢波的时间延迟,并结合已发表的固定台站的结果,获得了华北上地幔各向异性图像.从得到结果看,华北东部各向异性快波方向基本为NWW-SEE方向,而西部的快波方向转到NW-SE或NNW-SSE.快、慢波时间延迟范围是0.50~1.47 s,华北西部的平均快、慢波时间延迟小于华北东部.在华北东部,快波方向与绝对板块运动(APM)方向基本一致,预示了NWW向的软流圈地幔流是引起该区域上地幔各向异性的主要原因,它使得上地幔橄榄岩等晶体的晶格优势取向沿地幔物质流动方向,从而导致了NWW趋向的快波方向.然而,在稳定的西部,快波方向既不与绝对板块运动方向一致,也不与构造走向一致,这种弱各向异性很可能是遗留在古老克拉通的厚的岩石圈内的"化石"各向异性.     

中国大陆地壳与上地幔地震各向异性研究

地壳与上地幔各向异性研究在地球动力学领域有广泛的应用和重要的科学意义,其研究有助于许多地质和地球物理基本问题的解释.中国大陆的各向异性研究,在地壳与上地幔两个尺度上都有许多重要进展,并在基本理论和岩石实验研究中取得进步.文章认为,各向异性研究揭示的科学问题将会是广泛的和深远的.     

中国东北地区北部上地幔各向异性及其动力学意义

中国东北地区广泛发育新生代板内火山,晚中生代以来岩石圈遭受过多期拉张作用.作为中国唯一的深震孕育区,中国东北地区受到太平洋板块的西向俯冲,使得其成为研究岩石圈变形、板块俯冲和板内火山成因及其相互作用关系的天然实验室.通过分析架设在中国东北地区北部的147个流动和固定台站的SKS波形数据,共计得到了377对各向异性参数和251个无效分裂结果.结果表明,中国东北地区东西两侧具有不同的各向异性分布:西部地区各向异性方向变化范围为N143-199°E,平均N169°E,与晚中生代岩石圈伸展方向一致;其各向异性延迟时间平均值约为0.8s,说明来自地幔的各向异性比较微弱,主要由残留在岩石圈中的古老变形所引起.同时,在松辽盆地和佳木斯地块部分区域,观测到延迟时间较小的各向异性(~0.4s),可能是由于岩石圈的拆沉和热地幔物质的上涌侵蚀了保留在岩石圈的古老形变所致.在研究区东部,NNW-SSE朝向的各向异性被观测到,并伴随较大的延迟时间(大于1.0s),可能与太平洋板块撕裂回撤而产生的地幔流动有关.此外,近W-E方向的各向异性只在佳木斯地块被观测到,而太平洋板块在地幔过渡带中的俯冲可能是其产生的主要成因.     

中国东北地区SKS分裂的上地幔各向异性结构与动力学

本研究收集了中国东北地区2008—2016年九年时间内207个固定地震台站和127个NECESSArray流动地震台站的波形资料,利用SKS波分裂的最小切向能量网格搜索方法获得了243个台站的有效分裂结果.研究结果显示,尽管研究区各向异性快波方向基本以NW-SE向为主,但无论是在快波方向上还是快慢波时间延迟上不同构造单元内部与不同构造单元之间均存在着较大差别.大兴安岭造山带北部的各向异性快波方向自北向南由NNE-SSW向转变为NNW-SSE向,在中部以NW-SE向为主,而南部自北而南由NE-SW向逐渐转变为近E-W向;松辽盆地的各向异性快波方向在北部自西向东主要表现为由NNW-SSE向逐渐转变为NW-SE向,在中部自西向东由NE-SW向转变为近E-W向,而在南部既有NE-SW向又有NW-SE向;佳木斯地块各向异性方向由西部的NW-SE转变为东部的NNW-SSE,同时快慢波时间延迟逐渐变大;长白山造山带北部自北向南由NW-SE向逐渐转变为近E-W向,中部各向异性快波方向为NNW-SSE向,且快慢波时间延迟较大,而南部以NW-SE向为主;燕山造山带的各向异性快波方向主要沿E-W向分布,基本平行于燕山造山带的走向.这些结果说明,尽管复杂的各向异性快波方向与局部岩石圈拆沉和热物质上涌有关,但更重要是与"大地幔楔"中物质水平流等动力过程密切相关,也有待将来结合更多地震资料如面波不同深度的特征各向异性进行分析.在阿巴嘎火山群、哈拉哈火山群、长白山火山、龙岗火山和镜泊湖火山区及五大连池火山区等特殊构造区的周边地区,各向异性快波方向围绕这些构造区随方位均发生明显变化,暗示了火山区下方热物质上涌可能影响了"大地幔楔"中的软流圈物质水平流方向.     

Upper mantle anisotropy and crust-mantle deformation pattern beneath the Chinese mainland

Over the past 10 years,the number of broadband seismic stations in China has increased significantly.The broadband seismic records contain information about shear-wave splitting which plays an important role in revealing the upper mantle anisotropy in the Chinese mainland.Based on teleseismic SKS and SKKS phases recorded in the seismic stations,we used the analytical method of minimum transverse energy to determine the fast wave polarization direction and delay time of shear-wave splitting.We also collected results of shear-wave splitting in China and the surrounding regions from previously published papers.From the combined dataset we formed a shear-wave splitting dataset containing 1020 parameter pairs.These splitting parameters reveal the complexity of the upper mantle anisotropy image.Our statistical analysis indicates stronger upper mantle anisotropy in the Chinese mainland,with an average shear-wave time delay of 0.95 s;the anisotropy in the western region is slightly larger(1.01 s)than in the eastern region(0.92 s).On a larger scale,the SKS splitting and surface deformation data in the Tibetan Plateau and the Tianshan region jointly support the lithospheric deformation mode,i.e.the crust-lithospheric mantle coherent deformation.In eastern China,the average fast-wave direction is approximately parallel to the direction of the absolute plate motion;thus,the upper mantle anisotropy can be attributed to the asthenospheric flow.The area from the Ordos block to the Sichuan Basin in central China is the transition zone of deformation modes between the east and the west regions,where the anisotropy images are more complicated,exhibiting"fossil"anisotropy and/or two-layer anisotropy.The collision between the Indian Plate and the Eurasian Plate is the main factor of upper mantle anisotropy in the western region of the Chinese mainland,while the upper mantle anisotropy in the eastern region is related to the subduction of the Pacific Plate and the Philippine Sea Plate beneath the Eurasian Plate.     

Any P-wave kinematic algorithm for vertical transversely isotropic media can be ADAPTed to moderately anisotropic media of arbitrary symmetry type

Most of P-wave anisotropic kinematic algorithms (modeling, processing, and inversion) have been developed for the case of Transverse Isotropy (TI). Does it mean that when dealing with more complex symmetry types (Arbitrarily tilted TI, orthorhombic, monoclinic or even triclinic), all these algorithms are irrelevant? In fact, not at all. It has recently been demonstrated that in 2D geometry any qP-wave TI kinematic algorithm can be simply generalized to the case of monoclinic symmetry using the so-called Azimuthally Dependent Anisotropy Parameter Transformation (ADAPT), assuming moderate anisotropy. The extension of the technique to the case of arbitrary anisotropy type (triclinic) is achieved in this paper. The method is successfully checked for seismic modeling in a full 2D model with layers of contrasted anisotropy types and with arbitrary vertical and horizontal velocity variations (non-constant gradient). Typically, the approximate travel times using ADAPT differ from the exact travel times by a few milliseconds for total travel times of the order of a few seconds. Applications to seismic processing are also described. The simplicity of the procedure and the generality of the applicability of the ADAPT recipe are striking and very convenient for practical applications. They certainly deserve further analysis.