FDMOD–声波方程有限差分正演模拟二维

FDMOD –声波方程有限差分正演模拟（二维）

格式：

fdmod wfile nx= nz= tmax= xs= zs= [optional parameters] 必需的参数：

>wfile 波场输出文件（包含每个时间步的波场值wave[nx][nz]）nx= x采样点个数(第二维)

nz= z采样点个数(第一维)

xs= 炮点x坐标

dxs= 炮点x坐标间隔

zs= 炮点z坐标

dzs= 炮点z坐标间隔

ns= 炮点个数

tmax= 最大记录时间

可选参数：

nt=1+tmax/dt 时间采样点数(dt决定结果的稳定度)

mt=1 波场输出时间切片的时间步长间隔

dx=1.0 x采样间隔

fx=0.0 x起始值

dz=1.0 z采样间隔

fz=0.0 z起始值

fmax = vmin/(10.0*h) 震源子波的最高频率

fpeak=0.5*fmax 雷克子波的峰值频率

dfile= 密度输入文件（包含密度值d[nx][nz]）

vsx= 垂直测线的x坐标

hsz= 水平测线的z坐标

rsx= 水平测线的起始检波器x坐标

rlen= 水平测线长度

rivl= 水平测线检波器采样间隔

vsfile= 垂直测线的输出文件data[nz][nt]

hsfile= 水平测线的输出文件data[nx][nt]

ssfile= 震源点检波器的输出文件data[nt]

verbose=0 =1 显示输出信息=2 更多输出信息

abs=1,1,1,1 模型的顶,底,左,右使用吸收边界条件

=0,1,1,1 顶部使用自由边界条件

PML 参数:

pml_max=1000.0 完全匹配边界条件参数

pml_thick=0 完全匹配半层厚度(0 = 不使用PML边界条件)

注意:

程序使用普通显式二维差分方法