摘要
叠前逆时偏移需要计算成像点的初至时间。射线追踪法因为计算效率高,常被用作旅行时的计算,但难以处理速度剧变和非规则网格情形。常用的有限元法在求解波动方程时,也会像有限差分法一样存在数据频散问题。本文采用有限元解波动方程获取地震波初至时间,与一般有限元法的主要区别在于利用集中质量矩阵压制时间频散使波前形态保真,并给出了根据震源子波第一极值出现时间求初至时间的方法。数值实验结果表明该方法能适用于非规则网格和速度剧变界面情形。
Prestack reverse time migration requires computing the first breaks time of imaging points. The ray-tracing approach is usually used for the computation of travel time because of higher computational efficiency, but difficult to handle the rapid velocity variation and irregular grid. Common-used finite-element algorithm also has same data dispersion problem as the finite-difference algorithm when solving the wave equation. Different from the common finite-element algorithm, the finite-element algorithm presented in the paper for solving the wave equation and gaining seismic first breaks time is using lumped mass matrix to suppress time dispersion,which makes the feature of wavefront hi-fi. The paper gives the method for computation of first breaks time according to first extreme appearing time of source wavelet. The numeric simulation results showed that the approach is suitable for the irregular grid and rapid-changed velocity interface.
出处
《石油地球物理勘探》
EI
CSCD
北大核心
2007年第3期259-262,共4页
Oil Geophysical Prospecting
基金
国家重点研究规划"973"项目(2001CB209105)资助
关键词
波动方程
有限元法
初至时间
集中质量矩阵
数据处理
逆时偏移
wave equation, finite element algorithm, first breaks time, lumped mass matrix, data processing, reverse time migration
