ntraces);% 初始化地震数据矩阵%zr = 200; %检波器的位置在第200行zr=2;%检波器的位置在第2行% 模拟波传播forit=1:nt% 在源点注入 Ricker 波u(x0,z0)=u(x0,z0)+rick(it);% 将源函数注入波场% 计算波场的下一时刻(有限差分法,2D波动方程)u_new(:,:)=0;% 每次计算时先清空新波场
一、三维波动方程有限差分正演方法原理 三维波动方程的一般形式可以表示为: \[ \frac{\partial^2 p}{\partial t^2} = \nabla^2 p + f(x,y,z,t) \] 其中$p$表示波场,$f(x,y,z,t)$表示源项函数,$\nabla^2$表示拉普拉斯算子。对于三维波动方程的有限差分正演方法,其基本的数值离散形式如下: \[ ...
三维波动方程可以表示为:  其中,u(x,y,z,t)为波动场,c(x,y,z)为介质速度,ρ(x,y,z)为介质密度,t为时间。 有限差分正演方法通过将空间和时间离散化,将三维波动方程转化为差分方程,进而求解波动场在不同时刻的数值解。具体而言,有限差分正演方法将空间和时间分别离散化...
内容提示: 2015年 4月 油气地球物理 PETROLEUM GEOPHYSICS 第 13卷 第2期 波动方程法地震波正演数值模拟研究综述 张睿璇 ,廉西猛 (胜利油田分公司物探研究院,山东 东营 257022) 摘要 :波动方程法地震波正演数值模拟主要包括有限差分法、连续有限元法,谱方法以及新兴的间断有限元法。 本文总结了这些方法近年来的...
1、高阶弹性波波动方程正演模拟及逆时偏移成像 本文主要研究二维、三维弹性波高阶正演数值模拟和高精度弹性波叠前、叠后逆时深度偏移成像两个方面内容。弹性波正演数值模拟是基于弹性波动理论,能够保留地震波的运动学、动力学特征,是认识地震波传播规律的重要手段,可为地震波场的识别提供理论依据,同时也是弹性波叠前、...
( E 0 P H Y S IC A L P R 0 S P E CT l N G F O R P E T R OL E U M S ep. J9 90 一维波动方程两种正演数值方法的 关系 及反 演唐立 民刘 迎 曦姚 曙光 ( 大连理 工大学 ) _ 擂i 要 本文证 明 一维波动方程 的特征线法正演 ‘ 与合成地震记录 0 的一致性;给出了 一...
声波波动方程正演模拟程序总结声波波动方程正演模拟程序 程序介绍: 第一部分:加载震源,此处选用雷克子波当作震源。 编写震源程序后,我将输出的数据复制,然后我用excel做成了图片,以检验程序编写是否正确。以下为雷克子波公式部分的程序: for(it=0;it<Nt;it++) { t1=it*dt; t2=t1-t0; source[it]=(1.0-2.0...
实验题目:叠前地震记录的相移波动方程正演模拟实验 姓名 学号 班级 专业 勘测技术与工程(石油物探) 院(系) 地球物理学院 地球探测与信息技术系 单项成绩 内容理解 写作结构 程序设计 模型设计 计算结果 结果分析 总成绩 实验报告 一、实验题目 叠前地震记录的相移波动方程正演模拟实验 二、实验目的 掌握各向同性介质...
343 波动方程偏移方法二 波动方程偏移的成像原理 热度: 二维弹性波动方程F—K正演与偏移 奇恒定痒.正演可 生步长廷拓I儡移是特震源脉冲和地面记录波蝎同时下延,燕后在每 十深度作时f时事延迟互 相关来完直.由于建拓和成簟在囊率嵌盘置(F—K)中进行,目而算法 ...
三维波动方程延拓方法实现了三维地质模型的快速叠后正演。该方法可以采用相位移加插值方法处理一定的 横向变速情况,可以更加灵活方便地模拟地下复杂的三维地质体。首先进行了三维French模型数值模拟,得到 了和实际物理模型实验结果相一致的正演记录,并对比分析了三维偏移剖面和二维偏移剖面的偏移效果;然后 ...