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

权利要求书

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

1)对原始地震数据进行预处理，并对预处理后的地震数据分选共炮道集和共检波点道集，预处理包括滤波、去噪和反褶积；

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

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

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

简化为：

AR＝BT_s

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

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

简化为：

CS＝DT_r

其中，s_i表示不重复的炮点i对应的地表一致性高频静校正量，为不重复的炮点i的炮道集中对应不重复的检波点j的初至波到达时，C为炮点高频静校正量二阶差分系数矩阵，S为炮点高频静校正量列向量，D为共检波点道集初至波到达时二阶差分系数矩阵，T_r为共检波点道集初至波到达时列向量；

5)利用阻尼最小二乘法求解方程组分别获得检波点端和炮点端高频静校正量，根据线性反演理论，检波点端高频静校正量反演基本方程组的阻尼最小二乘解为：

R＝(A^TA+λI)^-1A^TBT_s

根据线性反演理论，炮点端高频静校正量反演基本方程组的阻尼最小二乘解为：

S＝(C^TC+λI)^-1C^TDT_r

其中，λ为阻尼因子；

6)对预处理过的地震数据施加地表一致性高频校正量，实现剩余静校正，地表一致性高频静校正量由炮点端的地表一致性高频静校正量与检波点端的地表一致性高频静校正量求和得到，则有：

t_i,j＝s_i+r_j

其中，t_i,j为不重复的炮点i对于不重复的检波点j的地表一致性高频静校正量，s_i为不重复的炮点i对应的地表一致性高频静校正量，r_j为不重复的检波点j对应的地表一致性高频静校正量；

实现该高频静校正量提取方法的系统包括处理器(2)以及分别与处理器(2)连接的数据采集器(1)、输入设备(4)和显示器(3)，所述的数据采集器(1)用以采集原始地震数据，并对原始地震数据进行预处理，所述的输入设备(4)用以提供分选之后的共炮道集与共检波点道集，并输入拾取的初至波到达时，所述的处理器(2)根据初至波到达时建波点端和炮点端高频静校正量反演基本方程组，并采用阻尼最小二乘法进行求解获得检波点端和炮点端的地表一致性高频静校正量，在地震数据上施加地表一致性高频静校正量后，通过显示器(3)显示校正后结果。