[发明专利]一种全波形反演方法和系统

权利要求书

1.一种全波形反演方法，其特征在于，包括以下步骤：

基于陆上地震资料进行反演获得地震剖面的背景速度场；

根据已知井的测井资料获得地震剖面中的井速度矢量；

利用井速度矢量对背景速度场进行插值处理，得到在低频范围内与测井资料耦合的初始速度模型；

基于该初始速度模型进行正演计算得到速度扰动模型，根据速度扰动模型更新初始速度模型，确定在背景速度场中低频信息缺失情况下的反演模型；

其中，所述利用井速度矢量对背景速度场进行插值处理，得到在低频范围内与测井资料耦合的初始速度模型包括：

确定地震剖面中位置u处的井速度矢量在位置i处的高斯权重值G(i)；

选取用于调节井速度矢量在背景速度场中权重的权重系数λ；

根据该高斯权重值G(i)和权重系数λ分别得到井速度矢量和背景速度场的组合分量；

对井速度矢量和背景速度场的组合分量求和得到初始速度模型m(i,j)，j为地震剖面中的纵坐标。

2.根据权利要求1所述的方法，其特征在于，所述地震剖面中位置u处的井速度矢量在位置i处的高斯权重值G(i)为：

G(i)=1σ2πexp[-(i-u)22σ2],i=1...nx,u=1...nx;]]>

其中，u为已知井在地震剖面中的横坐标，i为地震剖面中的横坐标，σ²为方差，n_x为地震剖面中的横坐标最大值。

3.根据权利要求2所述的方法，其特征在于，所述初始速度模型为：

m(i,j)＝λG(i)m_log(j)+[1-λG(i)m_bg(i,j)]，i＝1…n_x,j＝1…n_z；

其中，λG(i)m_log(j)表示井速度矢量的组合分量，[1-λG(i)m_bg(i,j)]表示背景速度场的组合分量，m_bg(i,j)为背景速度场，m_log(j)为已知井速度矢量，λ为权重系数，i为地震剖面中的横坐标，j为地震剖面中的纵坐标，n_x为地震剖面中的横坐标最大值，n_z为地震剖面中的纵坐标最大值。

4.根据权利要求3所述的方法，其特征在于，所述基于该初始速度模型进行正演计算得到速度扰动模型，根据速度扰动模型更新初始速度模型，确定在该背景速度场中低频信息缺失情况下的反演模型包括：迭代执行以下步骤；

判断步骤，判断初始速度模型是否满足精度要求，若不能满足精度要求，则执行扰动模型计算步骤和模型更新步骤，若能够满足精度要求，则将更新后的初始速度模型确定为反演模型；

扰动模型计算步骤，基于初始速度模型进行正演计算，根据计算数据与陆上地震资料的数据残差得到速度扰动模型；

模型更新步骤，根据速度扰动模型更新初始速度模型。