[发明专利]近地表初至层析静校正方法及计算机可读存储介质

权利要求书

1.一种近地表初至层析静校正方法，包括：

获得地震叠前数据，建立初始近地表速度模型；

通过射线追踪计算初至旅行时和灵敏度矩阵L；

对计算初至旅行时进行时差校正，并获取旅行时残差Δt；

基于所述灵敏度矩阵L和所述旅行时残差Δt，建立层析方程组；

基于所述层析方程组，进行初至层析反演，迭代更新所述近地表速度模型。

2.根据权利要求1所述的近地表初至层析静校正方法，其中，进一步包括：

基于射线传播深度与偏移距的关系，实现对所述近地表速度模型的校正。

3.根据权利要求1所述的近地表初至层析静校正方法，其中，通过射线追踪计算初至旅行时为：

其中，s和r分别是炮点和检波点；

Δl_j为当前射线在j网格内的长度；

m_j为j网格内的慢度。

4.根据权利要求3所述的近地表初至层析静校正方法，其中，通过射线追踪计算初至旅行时和所述灵敏度矩阵L：

t＝Lm

其中，t为计算初至旅行时；

m为慢度，是速度的倒数。

5.根据权利要求3所述的近地表初至层析静校正方法，其中，通过校正后的初至旅行时与记录或观测的初至旅行时做残差获取所述旅行时残差Δt为：

Δt＝t_obs-t_cal

其中，t_obs为观测或记录的旅行时；

t_cal为校正后的计算初至旅行时。

6.根据权利要求5所述的近地表初至层析静校正方法，其中，时差校正后的计算初至旅行时为：

t_cal＝t(s,r)+Δt_s+Δt_r

其中，Δt_s为炮点时差校正；

Δt_r为检波点时差校正。

7.根据权利要求5所述的近地表初至层析静校正方法，其中，所述基于射线理论，获取炮点时差校正：

其中，d＝offset/2，为偏移距；

h为射线穿透深度，h＝radio*d，radio为射线追踪穿透深度与偏移距之间的比例因子；

Δh_s为炮点处真地表与平滑地表高程差；

v为模型速度。

8.根据权利要求5所述的近地表初至层析静校正方法，其中，所述基于射线理论，获取检波点时差校正：

其中，Δh_r为检波点处真地表与平滑地表高程差。

9.根据权利要求6所述的近地表初至层析静校正方法，其中，通过所述旅行时残差Δt和所述灵敏度矩阵L，建立所述层析方程组为：

Δt＝L_truem_true-Lm＝L·Δm

其中，t_obs为观测或记录的旅行时；

t_cal为校正后的计算初至旅行时；

m_true为观测或记录的慢度；

L_true＝L为灵敏度矩阵；

Δm为慢度的扰动。

10.一种计算机可读存储介质，其上存储计算机程序，其中，所述程序被处理器执行时实现以下步骤：

获得地震叠前数据，建立初始近地表速度模型；

通过射线追踪计算初至旅行时和灵敏度矩阵L；

对计算得到的初至旅行时进行时差校正，并获取旅行时残差Δt；

基于所述灵敏度矩阵L和所述旅行时残差Δt，建立层析方程组；

基于所述层析方程组，进行初至层析反演，迭代更新所述近地表速度模型。