复杂介质地震波场数值模拟的研究.pdfVIP

中国地球物理2011 (13)地球介质各向异性 复杂介质地震波场数值模拟研究 NumericalSimulationofSeismicWave in TectonicArea PropagationComplex 石玉涛1 高 原1+刘澜波2 ShiYuta01GaoYuanlLiuLanb02 and ofCivilEnvironmental 1．中国地震局地震预测研究所北京100036；2．Department Engineering， of Connecticut，Storrs，Connecticut06269，USA University 地震波对传播路径上的各向异性介质十分敏感，剪切波分裂是方位各向异性介质的最主要特征之一。 剪切波分裂得到地震各向异性参数，能有效地反映出地壳介质的裂隙的走向以及裂隙间的结构，以及区域 1991)。目前，地震波场的数值模拟研究是对于各向异性介质结构的研究是最有效的方法之一。 地壳介质中广泛存在各向异性特征，地震各向异性研究关系到地球结构、介质特性、应力应变等科学 Transverse 问题。由于成因机理不同，各向异性主要包括：周期性互薄层各向异性(VTI：Vertical Isotropy) 是具有竖向对称轴的横向各向同性：及EDA裂隙(Extensive 水平向对称轴的裂隙(HTI：HorizontalTransverse Isotropy)，称为方位各向异性或水平横向各向同性。在地 壳内，尤其在活跃的构造变形区，充满EDA裂隙的介质可以认为是一种很常见的介质。本文通过给出网 格节点的各向异性参数，建立各向异性模型，研究地震波在各向异性介质中的传播问题。 地震波数值模拟是对地下介质结构和相对应的物理参数已知的情况下，模拟研究地震波在地下各种介 质中的传播规律。本研究在地震波在的均匀、各向异性介质的传播理论的基础上，根据剪切波分裂结果建 TimeDomain 立各向异性(HTI)介质的模型，采用伪谱法(PseudospectralMethod)数值模拟地震波在各 向异性介质中传播过程(Liuetal。2006)。该方法精度高、节省计算时间并且对于计算区域内的介质性质变 化不用特别处理。在建立各向异性的模型上，采用五个Thomsen各向异性参数和非均匀介质的速度变化系 化描述了实际介质的各向异性性质。研究认为，利用非均匀各向同性介质表示各向异性介质是一个有效及 灵活的地震波场数值模拟的方法。 本研究根据汶川地震源区的台站实际观测资料得到剪切波分裂参数结果确立各向异性介质的模型，数 值模拟进行汶川地震源区地壳各向异性介质反演。结合实际观测的地震波形与模拟波形进行拟合，调整地 壳的介质结构，获得一个相对可信的地壳介质的速度模型。 在EDA理论模型的基础上，研究得到汶川地震源区各个台站的各向异性特征(石玉涛等，2009)，剪切波 分裂研究表明该地区存在强烈的各向异性特征，在震源区的北东段的快剪切波的偏振方向为北东向，而西南段 的快剪切波的偏振方向为北西向，各向异性程度设置在1．5％--4．5％之间。根据赵珠等人得到的龙门山地区地壳 介质的速度结构，确立汶川地震震源区地壳介质的初步模型。研究采用伪谱法进行地震波场的数值模拟，得到 各个台站的理论地震图。各个台站的理论地震图与实际观测资料拟合，得到观测体波和理论体波的到时差，以 及理论快、慢剪切波的时间延迟和观测快慢剪切