利用接收函数方法研究蒙古中南部地区地壳结构

2011年8月至2013年7月中国地震局地球物理研究所与蒙古科学院天文与地球物理研究中心在蒙古中南部区域布设了宽频带流动地震台阵,这为开展远东地区深部结构的精细探测提供了有利的数据基础.利用台阵记录的远震地震事件,采用P波接收函数的H-κ叠加分析和共转换点（CCP）叠加方法获得了台站下方的地壳厚度及平均波速比.结果显示研究区的地壳厚度介于39 km至45 km之间.整体上Moho面埋深从西北往东南方向逐渐变浅.在蒙古主线性构造两侧地壳厚度呈现区域性变化特征,东南部地区地壳厚度较薄,约为39 km,而西北部地区地壳较厚,达45 km,为此推测蒙古主构造线可能是地壳的一个陡变带.此外,研究地区地壳的平均波速比值（V_P/V_S）在1.70到1.79之间,均值为1.75,低于全球大陆的平均值1.78,这可能暗示着该区其地壳是缺少铁镁质的.研究还发现测线的西北与东南地区其地壳波速比值较高,推测是古生代铁镁质地壳的残留或是新生代岩浆底侵的反映.     

新疆东准噶尔花岗岩类岩石高温高压弹性波速度及其对地壳结构的约束

应用超声波反射-透射法,在最高压力为1.0 GPa（室温）,最高温度为700℃（1.0 GPa）的条件下对新疆东准噶尔地区的卡拉麦里花岗岩带和野马泉岩体的典型花岗岩类岩石（碱长花岗岩、碱性花岗岩、花岗闪长岩、二长花岗岩和石英闪长岩）的纵波速度（V_P）和横波速度（V_S）进行了测量.结果显示,在常温、压力0.4~1.0 GPa条件下,东准噶尔地区花岗岩类岩石的V_P和V_S均随压力呈线性增加,说明在这个压力段岩石中的微裂隙已基本闭合.室温、1.0 GPa时花岗岩类岩石的V_P是5.79~6.84 km·s^-1,V_S是3.26~3.85 km·s^-1.依据压力与V_P及压力与V_S的线性关系,拟合得到常温常压下花岗岩类岩石的纵波和横波压力系数分别是0.1568~0.4078 km/（s·GPa）和0.0722~0.3271 km/（s·GPa）,V_P0和V_S0分别是5.62~6.47 km·s^-1和3.15~3.75 km·s^-1.恒压1.0 GPa、室温到700℃条件下,花岗岩类岩石的V_P和V_S均随温度的升高呈线性降低,温度系数分别为（-3.41~-4.96）×10^-4 km/（s·℃）和（-0.88~-3.22）×10^-4 km/（s·℃）.利用实验获得的花岗岩类岩石的V_P0、V_S0及温度系数和压力系数,结合东准噶尔地区的地热资料,建立了V_P和V_S随深度变化的剖面.将获得的V_P和V_S-深度剖面与该区地球物理探测结果对比,发现东准噶尔地区的碱长花岗岩、碱性花岗岩、二长花岗岩和部分花岗闪长岩的V_P和V_S与该区上地壳速度吻合很好,同时这几种岩石的平均泊松比也与上地壳泊松比一致,因此我们认为这几种类型的岩石是该区上地壳的重要组成部分.另外,石英闪长岩的V_P和V_S均符合中地壳的速度,可能为中地壳中的一种岩石.     

一阶弹性波方程数值模拟中的混合吸收边界条件

Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性.     

间断有限元方法的数值频散分析及其波场模拟

数值求解波动方程是大尺度正演波场模拟、基于波动方程的地震偏移和反演成像的关键.本文针对求解二维声波方程的Runge-Kutta 间断有限元（RKDG）方法的数值频散问题,从理论推导和数值分析的角度进行了深入研究,并将其与近似解析离散化方法（Optimal Nearly Analytic Discrete Method,简称ONAD 方法）、Lax-Wendroff 修正方法、交错网格（Staggered-Grid,简称SG）方法的数值频散进行了比较研究.结果表明：RKDG方法以及近似解析离散化方法在压制数值频散方面要好于上述其他方法,特别是空间精度为3阶的RKDG方法,即使当空间步长取波长的一半,即一个波长内取2个网格点时,最大的频散误差也不超过1.67%.同时,我们也通过波场模拟对比研究了不同数值方法的数值频散问题,进一步直观地验证了数值频散的理论分析结果.     

黑龙江及邻区地壳应力场初步研究

采用综合震源机制方法,利用5 002条P波极性数据反演黑龙江地区二维应力结构,求得该区0.5°×0.5°综合震源机制解.初步研究表明：21世纪以来,构造应力场总特征基本一致,具有明显的继承性,主压应力和主张应力的优势方向明显,与板块运动关系具有成因上的联系.     

浙江长兴深井地电台建设

随着人类的活动,地电阻率和地电场观测受到更多干扰。介绍浙江省长兴深井地电台建设过程,并讨论建设中存在的问题,初步分析地电观测资料,证明井下地电阻率和地电场观测在一定程度上可以屏蔽干扰,由此可见建设深井地电台的重要性。     

祁连山中东段地区中小地震震源机制和区域应力场

利用哥伦比亚大学GCMT目录给出的祁连山中东段地区中强地震震源机制资料,研究较大区域(34°-41°N,100°-106°E)的应力场;利用该地区布设的中法微震数字监测台网多年监测资料和甘肃数字监测台网资料,使用P波和S波初动及振幅比联合反演方法,反演中小地震震源机制解和发震应力场。结果表明,地区构造应力大致为北东40°-45°水平向压应力;景泰地区主压应力方向约北东45°,绝大多数地震为走滑型。天祝-古浪地区有相当部分的逆断层地震分布,主压应力方向约60°,P轴仰角在10°左右优势分布,大致为水平应力场。这与大区域构造应力场和断层实际分布基本一致。     

影区地震的识别与定位

通过分析西安地震台近几年各个方向记录到的影区地震,详细阐述 Lg1、Lg2波的确认方法,运用 Matlab对数据进行分析,给出S波与 Lg1、Lg2波的走时关系式,提高影区地震的分析速度和准确度,为单台大震速报和日常分析提供借鉴经验。     

泉州盆地地震效应特征的一维等效线性和二维非线性分析的比较

采用等效线性动粘弹性模型描述土的动力非线性特性,基于一维等效线性波传法,对泉州盆地地震效应进行了分析;同时,采用修正Martin-Seed-Davidenkov动粘弹塑性模型描述土的动力非线性特性,对泉州盆地非线性地震效应进行了大尺度二维精细化有限元分析,研究了地形地貌和土层横向不均匀性对地震效应的影响。将两种分析结果进行对比,结果表明：①随着基岩输入地震动强度增大,地表峰值加速度PGA放大效应总体呈现减小趋势,中震与小震、大震与小震的地表PGA放大系数之比依次为0.83~0.99、0.72~0.97;②该盆地Ⅲ类场地处,基岩、地表起伏不大,且土层横向分布较均匀,两种方法计算得到的地震效应特征类似;基岩或地表起伏剧烈、土层横向分布明显不均匀的Ⅱ类场地上,二维非线性分析给出的地表PGA放大系数明显大于一维等效线性结果,两种方法得到的地表加速度反应谱及PGA随土层深度的变化特征存在显著差异,二维非线性分析给出的地表加速度反应谱大多呈现双峰甚至多峰现象,且PGA在土层特定深度处存在聚集效应,使PGA随土层深度的变化呈现非单调性。     

2D各向异性TTI介质中地震波场的射线法和波动方程数值模拟对比

尽管射线理论是基于波动方程高频假设下的一种近似理论,两者本应相互佐证共同发展完善,然而实际上确是平行独立的.鉴于此,本文提出了一种可以相互验证波场数值模拟或射线追踪数值模拟方法正确性和有效性的佐证方法.为了实现上述目的,我们首先将用于各向同性介质中多次波预测的分区多步不规则最短路径算法推广至各向异性TTI介质中,使其具备追踪各向异性TTI介质中多次波的能力.然后,通过伪谱法波场数值模拟和分区多步不规则最短路径算法预测得到的波场快照图、单炮地震记录剖面、以及合成理论地震图的分析对比,不但可以相互验证两种算法的正确性和有效性,同时对地震波在各向异性TTI介质中的传播规律和特点有了深刻的理解和认识.结果表明:本文提出的各向异性介质中多次波射线追踪算法是正确有效的.