摘要：由于重力引起的岩石压实效应,一般来说,地震波传播速度由浅入深整体逐渐增大.梯形坐标系设计可耦合速度由浅入深逐渐增大的变化,该坐标系中均匀网格采样所对应的物理直角坐标系网格由浅入深逐渐增大,也即浅部低速区对应细网格,深部高速区对应粗网格.在梯形坐标系表征波动方程后利用有限差分求解,本文实现一种深度均匀采样、横向采样间隔随深度增加逐渐线性增大的有限差分地震波模拟方法.梯形坐标系波动方程离散后,仍采用常规均匀网格有限差分算法对其求解.由于横向网格大小由浅入深线性增加,本方法可避免不同大小网格区域过渡所产生的虚假反射.梯形坐标系波场模拟浅层精度高,深层横向响应范围广,可有效减少有限差分网格数量.本文提出的方法是在更广义的坐标系下利用有限差分求解波动方程,正交坐标系仅为该梯形坐标系之特例.本文旨在为大速度动态范围深地高效高精度地震波场模拟提供一种思路.

关键词：有限差分 变网格 梯形坐标变换 波动方程 

单位：西安交通大学数学与统计学院; 西安710049; 西安交通大学电子与信息工程学院; 西安710049; 海洋石油勘探国家工程实验室; 西安710049; 中国石化地球物理重点实验室; 南京211103