[发明专利]一种地震速度反演方法

权利要求书

1.一种地震速度反演方法，其特征在于，包括以下步骤：

S1，输入初始的观测地震数据、初始的速度模型与极性编码矩阵；

S2，对观测炮集利用极性编码矩阵合成编码超级炮集，其中，编码超级炮集包含多个超级炮记录；

S3，利用编码超级炮集与极性编码矩阵生成加权编码的参考道；

S4，利用参考道和目标子波对超级炮记录采用维纳滤波生成目标子波对应的超级炮；

S5，利用目标子波和极性编码矩阵与速度场生成正演超级炮；

S6，利用正演超级炮与步骤S4生成的超级炮记录做差反传，互相关得到反演梯度；

S7，利用抛物插值计算最佳步长，并对速度场进行更新；

S8，重复步骤S2到S7直到满足迭代次数或者速度场要求，输出最终的速度场。

2.根据权利要求1所述的方法，其特征在于，所述步骤S2合成编码超级炮集，具体公式如下：

其中，u_obs,n表示炮集第n炮的炮记录，e_n,k表示极性编码矩阵第k次迭代第n炮的极性值，其中编码矩阵是一个K×N_s的矩阵，其中值为随机生成的正负1,K为总迭代次数，N_s为编码的总炮数，U_su,k为第k次迭代合成的超级炮。

3.根据权利要求2所述的方法，其特征在于，所述步骤S3中的编码超级炮集与极性编码矩阵生成加权编码的参考道，具体公式如下：

U_su,k(x(n))为步骤S2中第k次迭代合成的超级炮，x(n)为第n个炮记录对应的炮点位置，U_ref,k为第k次迭代生成的加权编码的参考道。

4.根据权利要求3所述的一种地震速度反演方法，其特征在于，所述步骤S4中，利用参考道和目标子波对超级炮记录采用维纳滤波生成目标子波对应的超级炮，具体公式：

S_tar(w)为主频为w的目标子波，U_tar,k(w)为第k次迭代目标子波S_tar(w)所对应的超级炮，U_ref,k为应用步骤S3得到的第k次迭代的参考道，U_su,k为应用步骤S2中的公式生成的第k次迭代超级炮，ε为小于1的正数。