大地电磁测深与地震初至波走时交叉梯度反演

在地球物理勘探过程中,观测数据有限且存在误差,导致反演结果存在非唯一解。为了解决非唯一性问题,笔者研究了交叉梯度数学理论,实现了用交叉梯度约束2D 大地电磁和地震初至波走时的联合反演,并建立了2个模型。通过各自模型实验和不同模型间的对比,证明了交叉梯度联合反演比单独反演修复几何形态和物性参数的效果好,且交叉梯度与岩石物性无关;利用得到的交叉梯度图进一步证明了联合反演模型相似度高,也从交叉绘图中得到联合反演的物性关系更好的结论,反向证明了联合反演的正确性。     

PS转换波界面二次源法射线追踪

界面二次源法是最近提出的一种最小走时射线追踪方法,尤其适合层状介质中走时和射线路径的计算.该方法相对于传统的最小走时树方法（如Moser法）,仅在物性界面上设置二次源,射线路径的方向只在层界面处发生改变,该方法最大程度地消除了射线路径的锯齿状现象,同时也避免了低变速区的射线路径多值现象,因此,它具有更高的追踪精度和效率.本文采用界面二次源法在各向同性介质中实现了PS转换波射线追踪,理论模型的计算证实了界面二次源法追踪PS转换波的准确性和高效性,同时该方法在各向异性介质中也很好地追踪出分离的PSV波和PSH波,因此该方法有利于横波分裂在地震勘探中的研究和应用     

垂向非均匀介质中转换波转换点计算

转换点位置的计算是转换波资料处理中的一个关键问题. 本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法. 研究了通过速度拟合、走时近似和相似系数谱三种方式选择合适的等效速度方法. 结合理论模型对非均匀介质转换点计算方法、渐进转换点计算方法、Thomsen近似公式和均匀介质解析计算方法的误差进行了分析,结果表明非均匀介质转换点计算方法能更准确地计算转换点位置.     

地震波初至走时的计算方法综述

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

中国地震速报台网走时残差分析与走时修正

使用了３０个速报台站的观测数据,计算了１００多个中国境内外的地震事件初至波走时残差。在此基础上分析各台所记录的所有地震事件的残差随方位角和震中距的分布,及产生这些残差的地球深部非流体静力平衡状态和地球表面构造等方面的原因,并给出了不同区域地震走时的平均校正值,从而提高了地震定位精度。     

川滇地区地震动预测模型及其对2021年漾濞6.4级地震的适用性

随着川滇地区强震记录的不断增加,为了建立更符合该区域地震动特征的预测模型,文中基于该区域现有的地震动数据,通过随机效应回归模型建立适用于川滇地区的地震动预测模型;2021年5月21日,云南省大理州漾濞县发生6.4级地震,为了分析文中预测模型对漾濞地震的适用性,首先根据预测模型的适用范围选取合适的漾濞地震数据,计算真实记...     

多震相菲涅耳体射线走时同时反演成像

高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线附近一定范围(称其为第一菲涅耳体)内的速度异常分布有关.鉴于此,本文提出了计算多震相地震波菲涅耳体有限频射线的方法,并定义了走时敏感核函数,同时给出了利用多震相菲涅耳体有限频射线进行速度模型和反射界面同时反演成像的公式.利用多震相走时资料,使用传统射线层析成像方法与有限频射线层析成像方法进行了速度和界面的同时反演成像.结果表明,当射线密度较小时,无论是对速度模型的重建还是对反射界面几何形状的更新,有限频射线层析成像方法均优于传统射线层析成像方法, 而变频有限频射线层析成像则是实际地震层析成像的首选反演算法.      

地震叠前时间偏移的一种图形处理器提速实现方法

新近发展的图形处理器（GPU,Graphic Processing Unit）通用计算技术,现已日趋实用成型,并获得诸多应用领域的广泛关注.对油气勘探专项资料处理技术的运用而言,概因GPU与中央处理器（CPU）的计算性能的甚大差异,致使GPU这一通用计算技术在石油工业中的应用研究正在有效开展.本文仅借助于油气勘探中广泛使用的叠前时间偏移,旨在于扼要阐明其基于GPU应用的有效性;文中还提出一种利用GPU实现地震叠前时间偏移的软件构件方法,并针对非对称走时叠前时间偏移所拓展的应用软件提供一种具体实现架构.与以往用个人计算机（PC,Personal Computer）或者PC集群所用的叠前时间偏移相比,本文方法可甚大地提高计算效率,从而在石油物探资料处理中可显著地节约计算成本和维护费用.文中实际例证也表明,基于GPU进行高性能并行计算,当是适应目前石油工业中大规模计算需求的一个重要发展途径.     

地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.     

华北北部不同构造单元地壳速度模型——均匀分层介质反射波走时反演的结果

从广义反演理论出发,对多层反射波的走时联合反演,可以同时获得所有的各层速度及厚度参数。其优点是在不需要把问题线性化过程的同时,统一考虑了各层反射波走时的离散对整个地壳模型的影响,又可兼顾观测数据的不同误差分布。本文将该方法进行了数值模拟并应用于华北北部的几个不同构造单元,反演得出了其反射P波速度结构,并与正演结果进行了对比。