[发明专利]水平沙层地表的地震剖面模拟方法、装置、设备及介质

权利要求书

1.一种水平沙层地表的地震剖面模拟方法，其特征在于，包括：

生成单道含沙丘鸣震多次波信号；

基于深度域层速度模型计算时间域反射系数序列；

将时间域反射系数模型逐道与单道沙丘鸣震多次波信号做褶积运算，合成水平沙层地表条件下的地震剖面。

2.根据权利要求1所述的水平沙层地表的地震剖面模拟方法，其特征在于，其中由沙层厚度、沙层速度、潜水面反射系数、地震波振幅衰减因子和一次波信号，生成单道含沙丘鸣震多次波信号。

3.根据权利要求2所述的水平沙层地表的地震剖面模拟方法，其特征在于，所述单道含沙丘鸣震多次波信号为：

X(t)＝X_s(t)+(-1)rαX_s(t-τ)+(-1)²r²α²X_s(t-2τ)+(-1)³r³α³X_s(t-3τ)+… (2)

其中，自由地表反射系数近似为-1，潜水面反射系数为r，一个多次波周期内地震波振幅衰减因子为α，一次波信号为f(t)，X_s(t)表示震源端一次波和多次波信号，且有：

X_s(t)＝f(t)+(-1)(1-r)²αf(t-τ)+(-1)²(1-r)²rα²f(t-2τ)+(-1)³(1-r)²r²α³f(t-3τ)+… (3)

其中，τ为多次波周期。

4.根据权利要求3所述的水平沙层地表的地震剖面模拟方法，其特征在于，所述多次波周期为：

τ＝2h/c (1)

其中，h是沙层厚度，c是沙层速度。

5.根据权利要求1所述的水平沙层地表的地震剖面模拟方法，其特征在于，所述计算时间域反射系数序列包括：

对深度域层速度模型进行深时转换，将深度域层速度模型转换为时间域层速度模型；

由时间域层速度模型计算时间域反射系数序列。

6.根据权利要求5所述的水平沙层地表的地震剖面模拟方法，其特征在于，所述时间域反射系数序列中某点(i,j)的反射系数为：

其中，v_i,j表示点(i,j)的速度，v_i,j+1表示点(i,j)时间方向下一个点(i,j+1)的速度。

7.根据权利要求6所述的水平沙层地表的地震剖面模拟方法，其特征在于，所述时间域反射系数模型逐道与单道沙丘鸣震信号做褶积运算表示为：

S(x,t)＝R(x,t)*X(t) (5)

其中，R(x,t)表示空间x处的地下反射系数序列。

8.一种水平沙层地表的地震剖面模拟装置，其特征在于，包括：

信号单元，生成单道含沙丘鸣震多次波信号；

计算单元，基于深度域层速度模型计算时间域反射系数序列；

合成单元，将时间域反射系数模型逐道与单道沙丘鸣震多次波信号做褶积运算，合成水平沙层地表条件下的地震剖面。

9.一种电子设备，其特征在于，所述电子设备包括：

存储器，存储有可执行指令；

处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-7中任一项所述的水平沙层地表的地震剖面模拟方法。

10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的水平沙层地表的地震剖面模拟方法。