[发明专利]三维初至波剩余静校正方法及装置

权利要求书

1.一种三维初至波剩余静校正方法，其特征在于，包括：

利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，所述线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；

将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；

对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。

2.如权利要求1所述的方法，其特征在于，利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组，包括：

假设折射面水平，折射层为均匀介质，构建表征低速层厚度和三维折射初至波的旅行时之间对应关系的三维折射初至波的时距方程；

应用二层水平介质折射模型，建立炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系；

根据所述时距方程、炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系，建立三维折射初至波的旅行时与炮点静校正量、检波点静校正量和折射层慢度之间的关联关系；

对所述关联关系应用最小二乘法，建立以炮检距和三维折射初至时间为已知量，以炮点静校正量、检波点静校正量和综合折射层慢度为未知量的线性方程组。

3.如权利要求1所述的方法，其特征在于，所述线性方程组为：

其中，代表在给定的炮检距范围内，第n个炮点所包括的全部检波点之和，个数为N_n；代表在给定的炮检距范围内，第m个检波点所包括的全部炮点之和，个数为M_m；ST_n代表第n个炮点的炮点静校正量；ST_j表示第j个炮点的炮点静校正量；RT_j代表第j个检波点的检波点静校正量；RT_m代表第m个检波点的检波点静校正量；X_nj代表第n个炮点与第j个检波点之间的炮检距；X_im代表第i个炮点与第m个检波点之间的炮检距；P_n代表第n个炮点附近的综合折射层慢度；P_i代表第i个炮点附近的综合折射层慢度；t_nj代表第n炮、第j个检波点拾取的折射初至时间；t_im代表第i炮、第m个检波点拾取的折射初至时间。

4.如权利要求1所述的方法，其特征在于，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量，包括：

将迭代求解方程组进行变换，形成用于求解综合折射层慢度的慢度求解迭代公式；

对综合折射层慢度进行去高频分量处理，结合慢度求解迭代公式，确定三维观测平面上任意炮点处的综合折射层慢度值；

根据三维观测平面上任意炮点处的综合折射层慢度值和迭代求解方程组，进行迭代求解，得到炮点静校正量和检波点静校正量。

5.一种三维初至波剩余静校正装置，其特征在于，包括：

旅行时拟合模块，用于利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，所述线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；

迭代求解模块，用于将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；

去低频模块，用于对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。