本发明属于地球物理勘探数据处理技术领域，公开了一种有限差分数值模拟方法、系统、介质、计算机设备及应用，通过时间任意偶数阶精度有限差分逼近时间偏导数、空间隐式交错网格有限差分逼近空间偏导数、纵横波波场分离实现弹性P‑SV波高精度数值模拟。本发明的弹性P‑SV波高精度有限差分数值模拟方法，采用“十字+菱形”交错网格有限差分模板逼近弹性波方程中的时间偏导数；采样隐式交错网格有限差分法逼近弹性波方程中的空间偏导数；推导了相应的频散关系；基于Taylor级数展开求取差分系数；基于纵横波分离实现了P‑SV波高精度模拟。本发明可以大幅提高P‑SV波的模拟精度，提高弹性波逆时偏移和全波形反演的精度。