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

臧春娟  王明梁  刘永刚  唐红峰 《地球物理学报》2014,57(7):2099-2112

应用超声波反射-透射法,在最高压力为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均符合中地壳的速度,可能为中地壳中的一种岩石.  相似文献   

一阶弹性波方程数值模拟中的混合吸收边界条件   总被引：2，自引：1，他引：1       下载免费PDF全文

任志明  刘洋 《地球物理学报》2014,57(2):595-606

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

间断有限元方法的数值频散分析及其波场模拟          下载免费PDF全文

贺茜君  杨顶辉  吴昊 《地球物理学报》2014,57(3):906-917

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

黑龙江及邻区地壳应力场初步研究   总被引：1，自引：0，他引：1  

吴宝峰  徐建权  张永刚  郝永梅 《地震地磁观测与研究》2014,35(5):72-75

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

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

荣代潞 《地震地磁观测与研究》2014,(Z2)

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

影区地震的识别与定位     

周建东  ;舒优良  ;高自明  ;马世虎  ;闫俊义 《地震地磁观测与研究》2014,(1):50-54

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

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

董菲蕃  陈国兴  金丹丹 《地震学刊》2014,(2):154-160

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

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

王涛  张延腾  白超英 《地球物理学进展》2014,(2)

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

多次波压制技术在南海北部陆缘深水区的应用     

贾连凯  吴时国  赵昌垒 《地球物理学进展》2014,(2)

南海北部洋陆过渡带水深变化剧烈,海底地形崎岖,构造复杂多变,地震勘探中多次波十分发育,常规处理难以取得很好的成像效果.由于方法的局限性,任何一种多次波压制方法都很难彻底地压制全部多次波.本文通过对国家基金委共享航次在南海北部陆缘采集的二维反射地震资料的分析,总结了不同多次波压制技术的优缺点和适用范围,结合滤波类方法和波动理论类方法,综合应用了表面相关多次波衰减法(SRME,Surface Related Multiple Elimination)、τ-p域预测反褶积、高分辨率抛物线Radon变换以及叠后压制残余多次波等方法去除多次波,有效突出了一次反射波,证明了方法的有效性和实用性,同时也形成了针对南海深水资料多次波压制的一种有效的处理流程和方法组合.  相似文献   

地震波初至走时的计算方法综述   总被引：3，自引：0，他引：3  

赵烽帆  马婷  徐涛 《地球物理学进展》2014,(3)

在地震波场中,初至波到时信息由于初至震相可追踪、易识别性,在地震学领域占有重要的位置,广泛地应用于叠前偏移、叠前速度分析、地震走时层析成像及地震定位等.本文主要介绍了四类具有代表性的计算初至波走时的方法:(1)基于高频近似射线理论方法,如最短路径方法(SPM),及修正后的最短路径方法(MSPM);(2)基于程函方程的数值解方法,如有限差分方法(FD)、快速推进法(FMM)和快速扫描法(FSM);(3)基于惠更斯原理的波前构建法(WFC);(4)基于频率域波动方程数值解法(FWQ).最短路径方法计算精度较高,稳定性较好,但其需要采用更多的网格节点,因此计算效率低;程函方程数值解法无需计算射线路径,具有计算效率高、稳定性较好、易于实现等优势,但其计算精度较低,可以通过引入高阶差分格式得到提高;波前构建法计算精度高,稳定性好,但其需要在射线网格和规则网格之间做网格转换,因此计算效率较低;频率域波动方程方法能适应任意复杂介质,但其计算精度和计算效率较低.  相似文献