[发明专利]一种基于最小二乘反演的高频静校正量提取方法

说明书

本发明涉及一种基于最小二乘反演的高频静校正量提取方法，包括以下步骤：1)对原始地震数据进行预处理，并对预处理后的地震数据分选共炮道集和共检波点道集；2)拾取所有共炮道集和共检波点道集中每一道的初至波到达时；3)统计整个观测系统中不重复的炮点总数m、不重复的检波点总数n以及每个不重复炮点i和不重复的检波点j的坐标；4)分别构建检波点端和炮点端高频静校正量反演基本方程组；5)利用阻尼最小二乘法求解方程组分别获得检波点端和炮点端高频静校正量；6)对预处理过的地震数据施加地表一致性高频校正量，实现剩余静校正。包括以下步骤：与现有技术相比，本发明具有校正效果好、原理简单、实施简单等优点。

技术领域

本发明涉及地震勘探的数据处理领域，尤其是涉及一种基于最小二乘反演的高频静校正量提取方法。

背景技术

静校正是地震数据处理的关键步骤，静校正量提取的准确与否直接影响后续地震数据处理结果的优劣，静校正量可以分为低频分量与高频分量，低频分量通常与模型及地表起伏形态的低波数成份有关，其数值一般较大，低频分量相对来说比较容易获得，通常根据表层速度层析成像结果将观测数据校正至一平滑起伏面上来实现低频静校正分量的计算与施加；高频分量通常与模型及地表起伏形态的高频分量有关，其数值一般较小。

受地下速度模型反演精度的影响，高频分量往往难以准确得到，通常在静校正之后，或叠加之后，通过平滑反射同相轴实现高频分量的施加，但这样做存在一定的随意性与经验性，且无法实现地表一致性高频静校正。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于最小二乘反演的高频静校正量提取方法。

本发明的目的可以通过以下技术方案来实现：

一种基于最小二乘反演的高频静校正量提取方法，包括以下步骤：

1)对原始地震数据进行预处理，并对预处理后的地震数据分选共炮道集和共检波点道集；

2)拾取所有共炮道集和共检波点道集中每一道的初至波到达时；

3)统计整个观测系统中不重复的炮点总数m、不重复的检波点总数n以及每个不重复炮点i和不重复的检波点j的坐标；

4)分别构建检波点端和炮点端高频静校正量反演基本方程组；

5)利用阻尼最小二乘法求解方程组分别获得检波点端和炮点端高频静校正量；

6)对预处理过的地震数据施加地表一致性高频校正量，实现剩余静校正。

所述的步骤1)中，预处理包括滤波、去噪和反褶积。

所述的步骤4)中，根据共炮道集中拾取的初至波到达时，构建检波点端高频静校正量反演基本方程组，则有：

简化为：

AR＝BT_s

其中，r_j为不重复的检波点j对应的地表一致性高频静校正量，为不重复的炮点i的炮道集中对应不重复的检波点j的初至波到达时，A为检波点高频静校正量二阶差分系数矩阵，R为检波点高频静校正量列向量，B为炮道集初至波到达时二阶差分系数矩阵，T_s为炮道集初至波到达时列向量。

所述的步骤4)中，根据共检波点道集中拾取的初至波到达时，构建炮点端高频静校正量反演基本方程组，则有：

简化为：

CS＝DT_r