[发明专利]地震纵波分量多次波总照明度指标分析方法及装置

权利要求书

1.一种地震纵波分量多次波总照明度指标分析方法，其特征在于，所述方法包括：

获取工区中已确定的地震采集信息和弹性地震地质模型，对所述弹性地震地质模型进行正演预处理，得到离散化模型，所述地震采集信息包括震源子波信息和地震观测系统信息；

遍历地震观测系统中所有炮，得到所有炮的单炮纵波分量多次波照明度指标；

将地震观测系统中所有炮的单炮纵波分量多次波照明度指标进行叠加，得到纵波分量多次波总照明度指标；

其中，每一炮按如下方式得到单炮纵波分量多次波照明度指标：

根据所述离散化模型和所述震源子波信息，对单炮进行双程弹性波地震激发数值模拟，确定单炮在指定时刻的第一弹性波正应力场；

根据所述第一弹性波正应力场，确定单炮的每个接收点的纵波分量合成记录；

将每个接收点的纵波分量合成记录作为虚震源同时激发，进行双程弹性波动方程地震激发数值模拟，确定指定时刻的第二弹性波正应力场；

根据第二弹性波正应力场，确定指定时刻中每个时刻虚震源激发的单炮纵波分量多次波照明度指标。

2.根据权利要求1所述的方法，其特征在于，所述根据所述离散化模型和所述震源子波信息，对单炮进行双程弹性波地震激发数值模拟，确定单炮在指定时刻的第一弹性波正应力场，包括：

设置边界条件、初始条件和正演时窗范围，其中，所述边界条件用于限制所述离散化模型的边界，在所述初始条件时，第一弹性波正应力场数值为0；

根据所述离散化模型和所述震源子波信息，对单炮进行双程弹性波动方程地震激发数值模拟，得到所述正演时窗范围内指定时刻的第一弹性波正应力场。

3.根据权利要求1所述的方法，其特征在于，在将每个接收点的纵波分量合成记录作为虚震源同时激发之前，所述方法还包括：

重置第一弹性波正应力场和第二弹性波正应力场的数值为0。

4.根据权利要求1所述的方法，其特征在于，所述对所述弹性地震地质模型进行正演预处理，得到离散化模型，包括：

从所述弹性地震地质模型中获取纵波速度、横波速度和密度；

利用公式λ(r)＝ρ(r)(v_p²(r)-2v_s²(r))以及μ(r)＝ρ(r)v_s²(r)计算拉梅参数λ(r)和μ(r)；

其中，r为弹性地震地质模型内网格上任一点的坐标向量，ρ(r)为r点对应的密度，v_p(r)为r点对应的纵波速度，v_s(r)为r点对应的横波速度；

其中，纵波速度、横波速度、密度和拉梅参数构成离散化模型。

5.根据权利要求4所述的方法，其特征在于，确定第一弹性波正应力场时使用的震源为各向同性压力源，满足f(r_s,r,t)＝δ(r-r_s)s(t)，其中，f(r_s,r,t)为各向同性压力源，δ(r-r_s)为仅在r_s点不为0的空间脉冲函数，r_s为目标炮的炮点坐标矢量，r为弹性地震地质模型内网格上任一点的坐标向量，t用于表示采样时刻，s(t)为震源信号。

6.根据权利要求5所述的方法，其特征在于，所述根据所述第一弹性波正应力场，确定单炮的每个接收点的纵波分量合成记录，包括：

利用公式计算每个接收点的纵波分量合成记录；

其中，r_g为目标炮对应的任意一个接收点，P(r_s,r_g,t)为r_g点的纵波分量合成记录，λ(r_g)和μ(r_g)为r_g点的拉梅参数，τ_xx(r_s,r_g,t)、τ_zz(r_s,r_g,t)分别为t时刻的第一弹性波正应力场在x方向和z方向的分量。