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

权利要求书

1.一种利用波动方程获取地震波路径的方法，包括以下步骤：

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)继续进行计算；若频率循环完成，即已累加上所求频段的最后一个频率值，则存储地震波路径。

2.如权利要求1所述的一种利用波动方程获取地震波路径的方法，其特征在于：所述步骤3)和4)中，所使用的波动方程为以下二者之一：

▿2u(x)+k2u(x,ω)=0---(1)]]>

(∂∂z±ikz)u(x,ω)=0---(2)]]>

其中，x＝(x,y,z)为待勘探的复杂非均匀地下介质模型某一点处的空间位置，ω为地震波场的频率，u(x,ω)为x点处的频率域地震波场，为拉普拉斯算子，k为波数，k_z为垂直波数，i为虚数单位。

3.如权利要求1或2所述的一种利用波动方程获取地震波路径的方法，其特征在于：所述步骤6)中，地震波路径为