【摘 要】
：
地震波走时计算在观测系统设计、偏移成像、速度模型走时反演和地震定位等方面起到重要作用.各向异性广泛存在于地球介质中,影响地震波传播的振幅和走时,忽略各向异性的影响将对成像、反演以及地震定位等造成一定的误差.因此对于高分辨率成像和反演,走时计算中考虑各向异性十分重要.快速扫描法不需要存储和追踪波前面信息,在各向异性初至波走时计算方面应用广泛.传统的方法通过将慢度四次方程转换为走时四次方程并结合快速扫描法求解走时.该方法没有对程函方程做近似,适用于强各向异性介质,但存在计算效率低的问题.对于求解qSV波走时,
【机 构】
：
南方科技大学地球与空间科学系,广东深圳 518055;深圳市深远海油气勘探技术重点实验室(南方科技大学),广东深圳 518055;南方科技大学地球与空间科学系,广东深圳 518055;南方海洋科学与工
论文部分内容阅读
地震波走时计算在观测系统设计、偏移成像、速度模型走时反演和地震定位等方面起到重要作用.各向异性广泛存在于地球介质中,影响地震波传播的振幅和走时,忽略各向异性的影响将对成像、反演以及地震定位等造成一定的误差.因此对于高分辨率成像和反演,走时计算中考虑各向异性十分重要.快速扫描法不需要存储和追踪波前面信息,在各向异性初至波走时计算方面应用广泛.传统的方法通过将慢度四次方程转换为走时四次方程并结合快速扫描法求解走时.该方法没有对程函方程做近似,适用于强各向异性介质,但存在计算效率低的问题.对于求解qSV波走时,本文发展了一种在局部解中将慢度四次方程简化为二次方程解析地快速求解走时的方法,极大地提高了计算效率.对于qSH波,慢度方程是二次的,可以直接解析求解.最后,本文用各向异性均匀模型和BP复杂模型进行测试,计算结果表明走时计算准确,验证了该方法的有效性.
其他文献
地震动旋转分量可以提供传统平动分量所欠缺的空间梯度信息,理论上可以提升对震源参数的约束.本文对旋转分量用于矩张量反演进行了较为全面的数值试验.首先分析了均匀无限空间中位错点源产生的旋转波场与平动波场在辐射花样上的互补性,由此推测旋转分量用于矩张量反演能提供额外的约束.随后基于层状介质模型合成的平动和旋转记录,开展了地震矩张量反演测试研究.我们通过对比反演得到的矩张量、矩张量差异角△φ和地震矩,定量评估在矩张量反演中引入旋转分量的作用,发现在合适的条件下,在原平动记录基础上加入旋转记录来联合反演会显著提升反
本文选取自2004年8月至2008年5月12日汶川地震发生前,紫坪铺库区发生的2.1＜ML＜4.0的65次小震,分别计算了历次地震在汶川地震震源处产生的应力扰动量,随后根据不同机构给出的汶川地震震源机制解求取震源区的静态库仑应力改变量,讨论了紫坪铺库区小震与汶川地震的关系.在不同计算参数下,小震震群在汶川震源处产生的静态库仑应力扰动量之和为6±3 Pa,尽管分析认为小震群可能在汶川震源处产生较大的动态应力扰动,约为百帕量级,但均远未达到触发阈值,认为小震群并未直接对汶川地震有触发作用.推测小震群之间相互作
木里—盐源地区地处青藏高原东南缘,属于古特提斯洋构造域,是松潘—甘孜地块及扬子地块的交接地带,是研究青藏高原东南缘构造演化过程的重要区域.本文介绍的是横穿木里—盐源地区的大地电磁剖面,自北西向南东依次跨越锦屏山断裂、木里弧形构造区、丽江—小金河断裂、盐源盆地、金河—箐河断裂等构造.维性分析表明木里弧形构造区和金河—箐河断裂都表现为较强的三维性,因此本文采用大地电磁三维反演技术,获得了木里—盐源地区的精细电性结构.电性模型显示,沿剖面可以划分为4个主要的电性构造单元.锦屏山断裂以北的川西北次级地块下方10～
川东褶皱带作为华南板块中部的“侏罗山式”褶皱,开始形成于晚古生代.自白垩纪晚期开始,受太平洋板块和印度板块对欧亚大陆挤压的影响,这一构造带乃至华南板块中部又叠加了新的构造变形.但是,目前对于川东褶皱带白垩纪以来的构造演化缺乏足够的认识.位于川东褶皱带东侧、雪峰造山带西麓的沅麻盆地形成于早白垩世.晚白垩世以来,沅麻盆地与川东褶皱带处于同一构造应力场中,因此对盆地内早白垩世红层的古地磁研究对于解释川东褶皱带中生代晚期以来的构造演化有着重要的意义.该研究在沅麻盆地早白垩世红层中开展的古地磁学研究获得了可靠的原生
微地震矩张量能够描述储层岩石破坏的细节,有助于了解水力压裂储层的地质力学特征,震源机制矩张量反演对于非常规油气开发具有重要的作用.水力压裂通常在页岩油气储层中进行,震源区介质表现出速度和衰减各向异性,需要研究衰减各向异性介质中微地震信号的传播规律.基于观测到的微地震记录,通过波形反演能够获得矩张量,忽视震源区介质的衰减和速度各向异性会导致震源机制反演结果的误差.本文提出了一种考虑地层各向异性吸收衰减作用的震源机制全波形反演方法,对于黏弹性VTI介质,采用GPU并行的伪谱法计算合成地震记录,通过全波形反演方
地震P波、S波到时是精确分析地震水平位置、深度与速度结构等的重要参数,如何准确拾取P波和S波到时是地震学的一项重要的基础工作.大数据量与强噪声环境给地震到时的自动拾取带来了很大挑战.在频率域中可将信号与噪声分离,但会造成震相的偏移.针对上述问题,本文在STA/LTA、AIC方法的基础上,引入了标准时频变换(Normal Time-Frequency Transform,NTFT),结合信号时间域与频率域特征,提出了基于NTFT的STA/LTA方法,以及基于NTFT的AIC方法来拾取P波和S波的到时.基于N
相似性是一种常用的衡量不同图像之间差异程度的属性,广泛应用于地震数据处理环节.由于地震数据本质上是非平稳的,局部相似性比全局相似性更适用于刻画地震数据的时-空变化特征.现有的局部相似属性可以通过正则化最小二乘问题进行计算,但是其计算过程需要大量的计算时间和数据存储空间,难以适应当前的海量数据处理任务.本文提出了一种基于快速流式算法的局部余弦相似性计算方法,其采取局部数据递推的模式,避免迭代算法带来的计算负担,在保证计算精度的前提下能够快速地表征不同数据之间的差异.流式局部余弦相似性可以用于解决不同的地震数
远震接收函数中包含了大量台站下方速度间断面所产生的Ps转换波及其多次反射波的信息,被广泛应用于反演精细的地壳及上地幔速度结构.但当地壳结构复杂如存在沉积层或高速层时,转换波和多次波震相的到时和振幅发生变化,从而导致接收函数反演存在强非唯一性.深度学习作为一种高效的特征提取方法,能够建立接收函数与横波速度在空间上的映射关系.因此,利用全球模型数据以及高质量观测接收函数建立样本库,设计了利用接收函数预测横波速度的卷积神经网络.测试集结果表明,合成数据预测的横波速度与模型拟合程度较高,在实际数据检测中预测横波速
Radon变换是一种稀疏变换,被广泛应用于地震数据处理,其中线性Radon和抛物Radon最为常用.在实际地震数据中,直达波和面波的同相轴形态为线性,反射波为双曲型,单独使用线性Radon或抛物Radon变换时,不能确保所有同相轴在变换域的系数都是稀疏的,影响地震数据处理效果.本文提出的多路径Radon变换联合了线性Radon变换和抛物Radon变换,每个同相轴都有两种不同的变换参数来与积分路径相适应,能够兼顾不同形态的同相轴;然后利用最小二乘稀疏反演方法对不同形态同相轴匹配最佳积分路径,使其自动分离到两
本文介绍了我国最新研制的首套吊舱式高温超导全张量磁梯度测量系统,系统主要由高精度超导量子干涉磁强计构成的全张量磁梯度探头、数据采集系统、惯性导航系统、无磁吊舱等组成.在江苏省丹阳市试验区开展了飞行试验,成功地获取试验区全张量磁梯度分布图,测试结果表明吊舱式高温超导全张量磁梯度测量仪各项指标满足航磁测量技术规范的要求,测量精度优于±30 pT·m-1,具备了实际生产的能力,推动了我国航磁勘探的进步与发展,能够为航磁多参量数据采集-处理-解释系统提供宝贵的实测数据,填补我国航空高温超导全张量磁梯度测量技术研究