[发明专利]基于初至波信息的近地表建模方法及系统

权利要求书

1.一种基于初至波信息的近地表建模方法，包括：

1)基于叠前数据，拾取初至波旅行时；

2)建立反映地下介质速度分布的初始速度模型；

3)基于所述初始速度模型，进行最短路径射线追踪，获得相应的模型旅行时，并与拾取的所述初至波旅行时做差，获得旅行时残差；

4)基于所述初始速度模型，获得初至射线追踪计算核函数；

5)基于所述旅行时残差和所述核函数，建立层析反演方程组，求解所述层析反演方程组获得参数更新量；

6)基于所述参数更新量，对所述初始速度模型进行修改，获得迭代的结果；

7)当所述迭代的结果满足精度要求时，获得近地表模型；

在各向同性条件下，旅行时可以表示为沿着射线路径的Radon变换，将x，y域变为s，r激发接收点域，

t(s,r)＝∫_L(s,r)s(x,y)dl (2)

其中，t(s,r)表示初至波旅行时，s和r分别表示激发点和接收点，L(s,r)表示从激发点s到接收点r的射线路径，s(x,y)表示慢度，dl表示沿射线路径的距离增量，在离散情况下上式变为：

其中，Δl_j表示当前射线在j网格内的长度，m_j表示j网格内的慢度；

其中，所述层析反演方程组表示为：

Δt＝L·Δm

其中，Δt表示旅行时残差，L表示计算射线路径长度，Δm表示慢度残差。

2.根据权利要求1所述的基于初至波信息的近地表建模方法，其中，还包括：

8)当所述迭代的结果不满足精度要求时，重复步骤2)-6)，直到所述迭代的结果满足精度要求。

3.根据权利要求1所述的基于初至波信息的近地表建模方法，其中，所述叠前数据包括观测系统数据，激发点数据以及接收点高程。

4.根据权利要求1所述的基于初至波信息的近地表建模方法，其中，在进行最短路径射线追踪时分别以炮点、检波点为激发点计算正、反向旅行时。

5.一种基于初至波信息的近地表建模系统，包括：

A)用于基于叠前数据，拾取初至波旅行时的单元；

B)用于建立反映地下介质速度分布的初始速度模型的单元；

C)用于基于所述初始速度模型，进行最短路径射线追踪，获得相应的模型旅行时，并与拾取的所述初至波旅行时做差，获得旅行时残差的单元；

D)用于基于所述初始速度模型，获得初至射线追踪计算核函数的单元；

E)用于基于所述旅行时残差和所述核函数，建立层析反演方程组，求解所述层析反演方程组获得参数更新量的单元；

F)用于基于所述参数更新量，对所述初始速度模型进行修改，获得迭代的结果的单元；

G)用于当所述迭代的结果满足精度要求时，获得近地表模型的单元；

在各向同性条件下，旅行时可以表示为沿着射线路径的Radon变换，将x，y域变为s，r激发接收点域，

t(s,r)＝∫_L(s,r)s(x,y)dl (2)

其中，t(s,r)表示初至波旅行时，s和r分别表示激发点和接收点，L(s,r)表示从激发点s到接收点r的射线路径，s(x,y)表示慢度，dl表示沿射线路径的距离增量，在离散情况下上式变为：

其中，Δl_j表示当前射线在j网格内的长度，m_j表示j网格内的慢度；

其中，所述层析反演方程组表示为：

Δt＝L·Δm

其中，Δt表示旅行时残差，L表示计算射线路径长度，Δm表示慢度残差。

6.根据权利要求5所述的基于初至波信息的近地表建模系统，其中，还包括：

H)用于当所述迭代的结果不满足精度要求时，重复单元B)-F)，直到所述迭代的结果满足精度要求的单元。

7.根据权利要求5所述的基于初至波信息的近地表建模系统，其中，所述叠前数据包括观测系统数据，激发点数据以及接收点高程。

8.根据权利要求5所述的基于初至波信息的近地表建模系统，其中，在进行最短路径射线追踪时分别以炮点、检波点为激发点计算正、反向旅行时。