Marmousi模型二维数据 240道接收 750采样

推导了声波方程空间二阶导数的隐式求解公式及差分系数的求解方法,讨论了该方法的数值频散特征。利用该方法分别对均匀介质及Marmousi模型进行了数值模拟,将其结果与传统的显式差分格式的模拟结果进行了对比分析。结果表明:该方法较传统的显式求解方法具有更低的数值频散、更高的计算精度。

研究了二维任意广角单程声波方程叠前逆时深度偏移的吸收边界条件。底边界处波场值的计算采用内层网格的波场延拓算法实现,顶边界处采用虚网格加刚性边界条件来消除边界反射,对于左右边界,文中给出了波场计算方法,并推导了该边界条件逆时延拓的有限差分格式。均匀介质中的脉冲响应试验表明,只需要边界处镶边3层即可有效吸收入射到截断边界的外行波,消除边界反射。Marmousi模型的偏移处理取得了满意效果。

全波形反演是一种全新的地震成像方法,主要利用全波形信息反演地下介质参数,通过非线性优化波场理论值和观测值的残差实现波形反演。基于时间域声波方程,建立了时间域波场残差目标函数,以分层模型为例,分别从波场精度、目标函数收敛性和运行时间3个方面,比较了共轭梯度(Conjugate Gradient,CG)算法和拟牛顿算法(Broyden-Fletcher-Goldfarb-Shanno,BFGS)反演的效果。同时,应用共轭梯度法对正、逆断层模型和Marmousi模型进行了速度结构反演。反演结果表明:共