简介： 本发明涉及一种任意复杂介质地震波传播矩阵模拟方法，该方法包括以下步骤：1)将水平地表的非均匀地质模型剖分成厚度为Δz的若干水平薄层；2)在每个水平薄层内对广义Lippmann‑Schwinger积分方程进行傅里叶变换，推导反射/透射矩阵；3)利用水平薄层之间的边界连续条件，通过反射/透射矩阵推导广义反射/透射矩阵；4)使用层间迭代关系，使用广义反射/透射矩阵在频率域下进行地震波场的求解；5)使用傅里叶反变换将地震波场转换到时间域，达成所求目标。本发明将传统广义反射/透射矩阵法(推广到界面起伏的分层非均匀地质结构模型，推导的反射/透射矩阵中含有层内非均匀介质的影响，可以在保证精度的同时大大节省计算资源。