【摘 要】
:
The wavefront traveltime is widely used for seismic ray migration and tomography,while the seismic wavefields computed with wave equations are used for reverse-time migration (RTM) and full waveform i
【机 构】
:
School of Geoscience &Technology, Southwest Petroleum University, Chengdu 610500, China
【出 处】
:
2018勘探地球物理学研究进展学术研讨会
论文部分内容阅读
The wavefront traveltime is widely used for seismic ray migration and tomography,while the seismic wavefields computed with wave equations are used for reverse-time migration (RTM) and full waveform inversion (FWI).In this paper we tried to consider the two kinds of information in seismic numerical simulation at the same time.In order to improve the efficiency of calculation,first-arrival travel time discriminance (FATTD) is first established to estimate wavefront traveltime through wave equation forward modeling.Then,the time is used to detect the interesting spatial area of seismic modeling,with the goal of reducing computational complexity.To address the mass storage requirement of RTM/FWI,the first arrival time was used both to determine the storage region of source wavefield and to compress them,thus preventing the I/O delay of disk.As exemplified by a Hess model,this paper applies the FATTD method to seismic modeling and RTM.Numerical results show that the computational efficiency is improved by 30% and the memory demand is reduced to 1/461 of the original.
其他文献
以炮检点坐标和时间作为维度,利用此维度可定义任意地震道之间的广义距离,通过广义距离可索引待重建地震道的周边数据并利用压缩感知原理对其进行重建.由于此方法是单独对需要重建的每一道地震数据进行重建,该方法可用于常规压缩感知原理不适用的非均匀非规则网格场景.另外,该方法也可用于全新的规则观测系统的数据重建.
提出了反褶积同步挤压广义S变换,通过数学推导,对广义S变换结果进行“挤压”和重排,再将同步挤压广义S变换谱与窗函数Wigner-Ville分布进行二维反褶积可得到更高分辨率的信号时频分布.通过数值模拟结果可知,反褶积同步挤压广义S变换兼顾了广义S变换和同步挤压变换的优势,不仅能够极大地提高信号时频变换的分辨率,而且保留了多分辨率的特征.
将地震数据的高斯束重建与高斯束偏移结合起来,利用高斯束具有位置、振幅、斜率、曲率和一定宽度等特征化的地震同相轴特点,基于平面波分解理论实现地震数据的高斯束分解与重建.首先,通过平面波分解滤波器把信号分解为局部平面波;然后,构建高斯束来表征不同类型的局部平面波,迭代更新使残差最小,完成记录的重建;在此基础上,分析在含噪声和缺失道的情况下,地震数据的高斯束重建效果,并对重建数据进行高斯束深度偏移.测试
根据地震波的反射走时信息采用层析技术反演构建深度域速度模型,是实际地震资料生产中大规模叠前深度偏移速度建模的主流技术。通过研究基于叠前深度偏移的偏移距道集的速度建模方法,形成了剩余速度自动分析和拾取、反射控制点自动生成、构造倾角和构造连续性提取、构造倾角约束、偏移距双差走时层析等技术,解决了高陡复杂区域大规模深度域速度建模难题,实现了高效高精度的深度域速度建模方法和软件研发.以HLC三维区块为例,
提出了基于三维弥散黏滞波动方程的地震成像新方法—DVWEM(diffusive-viscous wave equation migration).DVWE的弥散—黏滞系数直接与含流体地层中传播的地震波的吸收衰减特性有关,因此,来自流体饱和储层的地震反射波的偏移成像处理对流体预测十分重要.DVWEM的波场延拓分两步进行;第一步延拓在频率—波数域实现,使用空间平均参数值,参数包括波的传播速度与介质的弥
地震波场数值模拟是研究复杂地区地震资料采集处理和解释的有效辅助手段,其中波动方程法能得到更为全面的波场信息,很好地保留了地震波传播的动力学特征.文章采用一阶速度一应力弹性波方程,使用更先进的交错网格,对非均匀各向同性介质P-SV波进行高阶有限差分法的声波与弹性波数值模拟.在边界条件的处理上,采用PML边界,并获得相关的波场快照与地震记录,并对弹性波场进行两种方式的波场分离,且取得良好的实验结果.有
波动方程数值模拟是勘探地震学的重要基础研究内容,也是逆时偏移和全波形反演的基础和核心。创造性地提出利用常规直角坐标系和旋转直角坐标系中的网格点一起差分近似波动方程中的Laplace算子,将Laplace算子表示为常规直角坐标系中M个Laplace算子和旋转直角坐标系中N个Laplace算子的加权平均,构建了一种新的混合2M+N型有限差分格式(简称为M2M+N-FD).频散分析表明:与T2M-FD和
推导出了TTI介质下qP波控制方程,同时基于拟微分算子分解思路,导出更符合物理意义的TTI介质纯qP波控制方程.新方程能够在物理上保持波传播能量,因此相对于正则化有限横波方程具有更好的稳定性.在数值算法方面,采用交高阶有限差分进行数值求解,给出差分的稳定性格式.给出均匀各向异性介质和复杂Foothill TTI介质模型中纯qP波波场模拟的结果,证明了新方程的稳定性.采用这个方程进行复杂TTI模型正
龙门山山区地势陡峭、地表起伏剧烈,表层破碎,岩性复杂,表层结构纵、横向变化大,且地下构造复杂,地层倾角大,属典型的“双复杂”勘探区域,从而造成次生干扰严重,激发接收条件差,资料信噪比低,成像困难等.因此,迫切要求研究复杂地表复杂地下构造条件下地质建模及地震波场特征和传播规律.首先,采用多尺度复杂地表复杂构造地质一地震建模方法,建立了龙门山山前带复杂地表复杂构造地质—地震模型.其次,采用任意起伏地表
储层含油气以后与围岩之间的波阻抗增大,所以含油气区域储层能量比较强,但是含油气储层的下方往往会出现能量衰减以及主频降低的现象,为此通过含流体模型的数值模拟研究了这种现象产生的原因.低频伴影现象的影响机理比较复杂,在实际研究中发现油层厚度与低频伴影之间存在着一定的非线性关系,利用数值模拟半定量分析了油层厚度与低频伴影之间的关系.