【seismic unix地震波时间域-有限差分正演】

弹性波正演简介

弹性波正演模拟是地震勘探和工程物探中的基础技术，用于模拟地下介质中弹性波的传播规律。Seismic Unix（SU）是一个开源地震数据处理软件包，其中sufdmod2函数可用于二维弹性波正演模拟。

安装与准备

确保系统已安装Seismic Unix软件包。可从官方仓库或GitHub克隆最新版本。编译时需要Fortran和C编译器支持。安装完成后，检查环境变量CWPROOT是否配置正确。

参数文件配置

sufdmod2需要输入参数文件（通常为.par文件）控制模拟参数。以下为关键参数示例：

n1=1000       # 时间采样点数
d1=0.001      # 时间采样间隔（秒）
n2=500        # x方向网格数
d2=5          # x方向网格间距（米）
n3=300        # z方向网格数
d3=5          # z方向网格间距（米）
vp=3000       # 纵波速度（米/秒）
vs=1500       # 横波速度（米/秒）
rho=2000      # 密度（kg/m³）

速度模型构建

创建包含介质参数的三维数组文件（.H格式）。使用unif2生成均匀模型，或通过math工具自定义复杂模型。示例命令：

unif2 < vp.H n1=500 n2=300 d1=5 d2=5 v=3000

震源设置

定义震源时间函数（如雷克子波）和位置参数。使用makewave生成震源信号：

makewave dt=0.001 nt=1000 fpeak=25 | suaddhead ns=1000 > wavelet.su

运行正演模拟

调用sufdmod2执行模拟，指定参数文件和输出文件：

sufdmod2 par=elastic.par vp=vp.H vs=vs.H rho=rho.H source=wavelet.su > output.su

结果可视化

使用SU的绘图工具查看结果。显示位移分量或波场快照：

suximage < output.su perc=99 title="Elastic Wavefield"

高级参数调整

对于复杂模型，可调整以下参数：

• abc=1 启用吸收边界条件
• nb=20 设置边界层厚度
• freesurf=1 模拟自由表面效应
• srctype=1 选择震源类型（0=力源，1=爆炸源）

并行计算优化

大规模模拟可使用MPI并行版本sufdmod2_mpi。运行前需配置MPI环境：

mpirun -np 4 sufdmod2_mpi par=elastic.par vp=vp.H vs=vs.H rho=rho.H > output.su

常见问题处理

若出现数值不稳定，尝试：

• 减小时间步长dt
• 检查CFL条件：dt ≤ min(dx,dz)/sqrt(vp² + vs²)
• 增加阻尼系数damp=0.01

注意：实际使用时需根据具体模型调整参数，建议通过小规模测试验证参数合理性后再进行完整模拟。