[发明专利]基于目标层可视性分析的激发位置确定方法

权利要求书

1.一种基于目标层可视性分析的激发位置确定方法，其特征在于，包括以下步骤：

（a）从地质模型中选取目标反射层；

（b）在地质模型的整个表面上均匀地布置多个勘探炮和多个检波器；

（c）使所有勘探炮分别产生地震波，并且当多个勘探炮中的一个勘探炮产生地震波时，多个检波器分别接收该地震波被目标反射层反射回来的反射波，同时记录每个检波器接收到的能量值，从而得到每个检波器分别针对当前产生地震波的勘探炮的炮检对的能量值；

（d）在所有勘探炮分别产生相同次数的地震波之后，基于多个检波器分别针对每个勘探炮的炮检对的能量值，绘制基于地质模型的炮检对可视性分析图；

（e）基于炮检对可视性分析图，确定在实际的地震勘探中应该增加勘探炮数量的范围，

其中，炮检对的能量值是指一个检波器接收到的一个地震波的能量值。

2.根据权利要求1所述的基于目标层可视性分析的激发位置确定方法，其特征在于，在步骤（c）中，利用有限差分算子，对当前产生的地震波进行正向延拓到目标反射层，然后再从目标反射层反向延拓到检波器，来模拟地震波被目标反射层反射回来的反射波，从而记录每个检波器接收到的反射波的能量值。

3.根据权利要求1所述的基于目标层可视性分析的激发位置确定方法，其特征在于，在步骤（c）中，当产生地震波时，地震波在频率域中以下面的公式（1）的形式传播，

（k²+Δ）U(r_s,r,ω)=-f(ω)δ(r-r_s)   公式（1）

其中，k=ω/c₀(r)为波数，ω表示地震波的频率，c₀(r)为位于目标反射层的点r=(x,y,z)处的背景场速度值，Δ为拉普拉斯算子，f(ω)为勘探炮函数，δ(r-r_s)是单位脉冲函数，r_s表示勘探炮的三维坐标。

4.根据权利要求3所述的基于目标层可视性分析的激发位置确定方法，其特征在于，在步骤（c）中，利用有限差分算子，对当前产生的地震波进行正向延拓到目标反射层，然后再从目标反射层反向延拓到检波器，

其中，检波器接收到的反向延拓的反射波的地震波场通过公式（2）来表示，

P(r_s,r_g,ω)=

2k²∫_V(r)G(r,r_s;ω)m(r)G(r,r_g;ω)d³r   公式（2）

其中，G(r,r_s;ω)和G(r,r_g;ω)分别是频率域中勘探炮和检波器位置的格林函数，r_g表示检波器位置的三维坐标，m(r)是位于目标反射层的点r=(x,y,z)处速度模型扰动，V(r)为目标反射层的所有点的集合。

5.根据权利要求4所述的基于目标层可视性分析的激发位置确定方法，其特征在于，多个检波器中的一个检波器针对当前产生地震波的勘探炮的炮检对的能量值通过公式（3）来表示，

P_sr(r_s,r_g,ω)=|P(r_s,r_g,ω)|   公式（3）。