地震波勘探方法是目前进行能源资源勘探开发中最为广泛应用的一种方法。在现实中,主要能源资源大多数都埋藏在具有复杂几何构造和复杂物理属性的地球岩石圈中,在工程勘探中,往往把含资源地层中的复杂地质构造简化成较为简单便于数值模拟的模型,简化的过程中难免产生误差使数值模拟得到的结果不准确,与实际的地下情况不符。所以,模拟出更符合实际地质结构的数值模拟图是快速解决当今地震勘探领域中问题的关键。在地震勘探中,勘探层往往表现为双相介质。所以,对双相介质波动理论的研究可以解决勘探中的难题,本文基于改进BISQ模型利用交错网格有限差分方法进行数值模拟。改进BISQ模型是结合Biot理论机制与喷射流理论机制,对理论计算公式进行简化的双相介质模型,能够更全面、更具体地反应了地下介质的真实情况的同时,计算量相对较少。交错网格有限差分方法被广泛的应用在双相介质地震波数值模拟中。论文基于改进BISQ模型,对二维( x , z )情况下一阶弹性波动方程利用交错网格有限差分方法进行差分离散,采用一网格点处的变量取相邻两网格点变量之和求平均的差分离散格式。利用PML完全匹配层进行吸收边界反射,对加入阻尼因子之后的波动方程同样利用上述差分格式进行离散,通过对固相波场进行数值模拟,观测到PML完全匹配层有效的吸收了人工边界反射。对用交错网格有限差分方法离散波动方程过程中产生的数值频散现象产生的原因及影响进行了分析,并采用FCT(通量校正传输)方法进行压制数值频散,讨论FCT方法中漫射系数的选定标准,通过固相波场数值模拟观测压制数值频散的效果。最后,对单层、双层地质模型进行数值模拟。通过观测数值模拟得到的波场快照及合成地震记录图,分析、讨论波场特征,对固相、流相中的快纵波、慢纵波、S波的特征进行分析讨论。同时观测FCT方法压制数值频散与PML完全匹配层吸收边界反射的效果。