三维波动方程时空域混合网格有限差分数值模拟方法

常规高阶和时空域高阶有限差分方法广泛应用于三维标量波动方程的数值模拟,这两种差分方法仅利用笛卡尔坐标系中的坐标轴网格点构建三维Laplace差分算子,相应的差分离散波动方程本质上仅具有2阶差分精度,模拟精度低.本文将三维笛卡尔坐标系中非坐标轴网格点分为两类:坐标平面内的非坐标轴网格点和坐标平面外的非坐标轴网格点,系统推...     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

弹性波方程正演的粗粒度并行算法

波动方程正演在地震资料采集、处理、解释与反演中均能发挥重要作用,但现有的基于求解地震波动方程的正演算法由于受庞大计算量的制约而难于大规模应用于工业生产.本文从二维弹性波动方程出发,研究了利用有限差分法并行求解该方程的基本思路与方法,给出了适于并行求解的计算空间划分方法与通信方案,分析了不同参数条件下并行程序的运行时间、加速比与效率.引入消息传递接口（MPI）实现了弹性波方程的并行求解,极大地提高了数值求解弹性波方程的计算效率.     

把握世界物探技术发展现状,促进地震资料处理、解释技术进步

第82届SEG年会成果表明:计算机技术和地震装备技术的发展,引领了国际地球物理勘探技术的快速发展;宽方位、宽频带、高密度等采集技术,各向异性逆时偏移、弹性波的全波形反演等处理技术,多波处理与解释、4D地震等技术,很快将成为物探的主导技术.为了更好地把握技术的发展,需要在地震资料处理、解释领域重点开展如下研究:发展宽频带、各向异性、最小二乘、弹性波的逆时偏移技术;层析、逆时偏移角度域道集速度分析、全波形反演的速度估计方法、是进一步提高成像精度重点发展方向.岩石物理是地震解释、储层反演、裂缝检测的基础,重点深化在实验室测量、弹性各向异性的测量、油藏描述和监测方向的研究.全波形反演技术是高精度成像建立速度模型和油藏描述重要手段,在理论上、在应用上仍是下一步发展的重要方向.     

二维弹性多波时空域高斯束偏移方法

随着我国勘探开发难度逐步增大,勘探目标开始向裂缝油气藏、岩性油气藏等复杂探区转移,研究高精度、适应性强的多波多分量深度偏移算法在后续的地震解释、属性分析及储层预测中具有重要意义.针对多波多分量地震数据,本文提出了一种二维弹性波时空域高斯束偏移方法.时空域高斯束沿中心射线传播时能够面向成像目标描述局部波场,且对振幅和频率可调制的Gabor基函数有天然的适应性,因而将基于Gabor分解的子波重构方法应用于震源波场构建,从而得到任意点源函数产生的时空域高斯束波场.该方法由于直接在时间域进行计算,可以避开频率域中出现的假频和边缘截断效应等问题.基于各向同性弹性波动方程的Kirchhoff-Helmholtz积分解,利用矢量时空域高斯束传播算子构建格林函数和格林位移张量,并结合上行射线追踪策略,实现了检波点波场的反向延拓.针对矢量波成像问题,本文借鉴弹性波逆时偏移方法从矢量延拓波场中分离出纯纵波分量和纯横波分量,进而采用修改后的内积成像条件产生具有明确物理意义的PP、PS成像结果,避免了转换波成像的极性反转问题.最后利用简单两层模型和不含盐体构造的部分Sigsbee2a模型的成像结果,并将其与应用近似纵横波成像条件、标量和矢量势成像条件的偏移剖面进行对比,验证了本文方法的正确性和有效性.

     

旋转网格和常规网格混合的时空域声波有限差分正演

压制数值频散,提高正演模拟精度,一直是有限差分正演模拟研究的重要内容.基于时空域频散关系的有限差分法,比基于空间域频散关系的传统有限差分法,模拟精度更高.时空域声波方程数值模拟,普遍采用常规十字交叉型高阶有限差分格式.而在频率-空间域,普遍采用旋转网格和常规网格混合的有限差分格式,有效提高了模拟精度和计算效率.本文将频率-空间域混合网格有限差分的思想引入到时空域,提出了时空域混合网格2M+N型声波方程有限差分方法.首先推导出基于时空域频散关系的混合网格差分系数计算方法,然后进行频散分析、稳定性分析,并和传统高阶、时空域高阶有限差分法对比,结果表明：计算量相同时,新方法能有效压制数值频散,显著提高模拟精度;新方法相比传统2M阶有限差分法,稳定性增强,与时空域2M阶有限差分法稳定性基本相当.最后利用新方法进行均匀介质、层状介质、盐丘模型的数值模拟和盐丘模型的逆时偏移,模拟效果和成像质量进一步证实了该方法的有效性和普遍适用性.     

基于机器学习的时间反转散射体检测方法

地下小尺度散射体的检测和识别对于提高地震勘探的分辨率具有重要意义,目前业界普遍采用绕射波分离及成像方法检测地下散射体,而绕射波成像的质量主要取决于绕射波和反射波波场分离的程度.本文将被动源震源定位问题中常用的时间反转原理引入到地下散射体检测中,首先通过分析被动源和主动源模型反传波场的聚焦状态,验证了时间反转原理应用于地...     

三维复杂黄土塬区变加密不等距网格初至走时计算与特征分析

三维复杂黄土塬区初至走时计算与特征分析,对于分析复杂黄土塬区地震波运动学特征、观测系统设计、初至拾取及层析反演等均具有重要意义.针对该问题：首先,采用三维分区局部变加密不等距网格剖分复杂黄土塬速度模型;然后,建立基于费马原理和迎风格式的双线性插值局部走时算法,采用修正后的群推进法作为波前扩展方式;最后,以理论模型考察了算法的精度与效率、分析了三维复杂黄土塬区地震初至波的走时特征.结果表明：变加密网格策略能在增加很少计算量的情况下,大幅降低因震源奇异性引起的计算误差问题,同时还能无条件稳定且灵活的适应三维复杂黄土塬区的复杂地震地质条件;和层状均匀模型相比,黄土塬区地震初至波的走时场和时距曲线较为复杂.