[发明专利]一种利用波动方程获取地震波路径的方法

摘要

本发明涉及一种利用波动方程获取地震波路径的方法，其包括以下步骤：1)给定初始参数；2)将待勘探的复杂非均匀地下介质模型划分为网状格点；3)将震源置于模型的震源点S处，逐点计算整个空间内所求频段内某个频率下的地震波场；4)将震源置于模型的接收点R处，逐点计算整个空间内所求频段内某个频率下的地震波场；5)在频率域内，将同一频率下在震源点S和接收点R处计算的地震波场相乘并与前一频率下得到的相乘结果进行累加；6)在所求频段内进行频率循环，若频率循环完成，则存储地震波路径。本发明避免了传统方法的高频近似和多路径问题，可以广泛应用于石油勘探复杂构造的地震模拟、地震照明和地震偏移技术中。

主权项

一种利用波动方程获取地震波路径的方法，包括以下步骤：1)给定初始参数，包括待勘探的复杂非均匀地下介质模型、地震波场的频率以及震源点S和接收点R的空间位置；2)将待勘探的复杂非均匀地下介质模型，以地表为起始点向着待勘探的方向按照不同的深度间隔和水平间隔划分为网状格点；3)将震源置于地下介质模型的震源点S处，利用波动方程逐点计算整个空间内所求频段内某个频率下震源点S处的地震波场u_S(x,ω)；4)将震源置于地下介质模型的接收点R处，利用波动方程逐点计算整个空间内所求频段内某个频率下接收点R处的地震波场u_R(x,ω)；5)在频率域内，将同一频率下震源点S处的地震波场u_S(x,ω)和接收点R处的地震波场u_R(x,ω)相乘并与前一频率下得到的相乘结果进行累加；6)在所求频段内进行频率循环，即重复步骤3)‑5)：若频率循环未完成，即未累加上所求频段的最后一个频率值，则返回步骤3)继续进行计算；若频率循环完成，即已累加上所求频段的最后一个频率值，则存储地震波路径。