本申请公开了基于DFL粘弹性方程的瑞雷波数值模拟方法及系统，属于地震波模拟技术领域，方法包括以下步骤：根据地震波在衰减介质中的传播特性，建立DFL粘弹性方程；根据面波产生条件和消除人工边界反射方法，设置真空自由表面条件和CPML边界条件；以所述DFL粘弹性方程、真空形式和CPML作为模型，采用SGPS法数值实现瑞雷波正演模拟。本申请以DFL粘弹性方程描述体波的传播，真空形式作为自由表面边界条件，CPML为吸收边界条件，采用SGPS法对浅层各向同性衰减介质进行了包括瑞雷波和体波在内的全波场模拟，该方法既继承了真空形式的简单易实现特点，又避免了使用真空形式时困扰基于GSLS方程的数值频散问题。