[发明专利]一种基于互相关道的散射波成像方法及系统

权利要求书

1.一种基于互相关道的散射波成像方法，其特征在于，该方法包括：

输入地震数据道集、速度模型以及初始化成像参数；

正演计算检波点到成像点的旅行时；

将地震数据道集转化为互相关道集；

利用互相关道集进行散射波成像；

输出散射波成像结果；

其中，将地震数据道集转化为互相关道集包括：

读入地震数据道集，道集中数据道数为N；

对于第i道数据道，求取第j道数据与第i道数据道间互相关函数，其中i＝1,2,...,N-1，j＝i+1,i+2,...,N；

改变i值，重复求取第j道数据与第i道数据道间互相关函数，直到任意两道间的互相关函数计算完毕，得到互相关道集；

其中，利用公式c_ij(t)＝∫s_i(τ)·s_j(τ+t)dτ求取第j道数据与第i道数据道间互相关函数，其中，c_ij(t)表示互相关道集，s_i(t)和s_j(t)分别表示第i和j道地震数据；

其中，利用互相关道集进行散射波成像包括：

读入任意一道互相关数据c_ij(t)，确定生成该互相关道数据的地震数据的两个检波器坐标r_i和r_j；

计算成像范围内任意成像点x到两个检波器间的旅行时差Δt_ij(x)；

对于任意成像点x，利用计算出的旅行时差Δt_ij(x)，在互相关数据道中提取对应时刻的采样点值；

将采样点值乘以振幅校正因子1/(r_i-x)(r_j-x)，输出和叠加到成像点x处；

重复上述步骤，直到相关道集中所有数据道计算完毕。

2.根据权利要求1所述的基于互相关道的散射波成像方法，其特征在于，正演计算检波点到成像点的旅行时包括：

根据地震数据道集的检波点坐标和成像范围，采用旅行时正演算法，计算检波点到所有成像点的旅行时。

3.一种基于互相关道的散射波成像系统，其特征在于，所述系统包括：

存储器，存储器存储地震数据道集、速度模型、初始化成像参数以及可执行指令；

处理器，所述处理器调用存储器存储的数据，执行所述可执行指令完成以下步骤：

正演计算检波点到成像点的旅行时；

将地震数据道集转化为互相关道集；

利用互相关道集进行散射波成像；

输出散射波成像结果；

其中，将地震数据道集转化为互相关道集包括：

读入地震数据道集，道集中数据道数为N；

对于第i道数据道，求取第j道数据与第i道数据道间互相关函数，其中i＝1,2,...,N-1，j＝i+1,i+2,...,N；

改变i值，重复求取第j道数据与第i道数据道间互相关函数，直到任意两道间的互相关函数计算完毕，得到互相关道集；

其中，利用公式c_ij(t)＝∫s_i(τ)·s_j(τ+t)dτ求取第j道数据与第i道数据道间互相关函数，其中，c_ij(t)表示互相关道集，s_i(t)和s_j(t)分别表示第i和j道地震数据；

其中，利用互相关道集进行散射波成像包括：

读入任意一道互相关数据c_ij(t)，确定生成该互相关道数据的地震数据的两个检波器坐标r_i和r_j；

计算成像范围内任意成像点x到两个检波器间的旅行时差Δt_ij(x)；

对于任意成像点x，利用计算出的旅行时差Δt_ij(x)，在互相关数据道中提取对应时刻的采样点值；

将采样点值乘以振幅校正因子1/(r_i-x)(r_j-x)，输出和叠加到成像点x处；

重复上述步骤，直到相关道集中所有数据道计算完毕。

4.根据权利要求3所述的基于互相关道的散射波成像系统，其特征在于，正演计算检波点到成像点的旅行时包括：

根据地震数据道集的检波点坐标和成像范围，采用旅行时正演算法，计算检波点到所有成像点的旅行时。