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

权利要求书

1.一种时间域地震层速度反演方法，其特征在于：所述方法首先通过地震测井获取初始瞬时速度，然后建立瞬时速度与均方根速度之间的关系，将所述初始瞬时速度转换成均方根速度，再用实际拾取的均方根速度进行约束构建目标函数，然后求解目标函数得到更新速度模型所需的变量，通过多次迭代后，当目标函数的值小于阀值ε时，得到最终的层速度。

2.根据权利要求1所述的时间域地震层速度反演方法，其特征在于：所述方法包括以下步骤：

第一步：获取初始瞬时速度：

将测井速度或平均瞬时速度作为初始瞬时速度；

第二步：构建速度更新方程如下：

V_n(t)＝V_n-1+ΔV*e^-kt             (1)

(1)式中，V_n(t)是待求解的层速度，ΔV是速度差，n代表第n次迭代，k为速度变化梯度，k的初始值设为零，*表示相乘，e为自然指数；

第三步：求取均方根速度，计算公式为：

Vrmscal(ΔV,k,t)=1t·∫0tVn2(τ)dτ=1t·∫0t(Vn2+2VnΔV·e-2kτ+ΔV2·e-2kt)dτ---(2)]]>

(2)式中，n的初始值为1，当n＝1时，设ΔV＝0，k＝0，则V₁＝V₀，V₀为第一步得到的初始瞬时速度；

第四步：构建目标函数如下：

Φ(ΔV,k)=Σi=1I[Vrmscal(ΔV,k,i)-Vidata]2→Min---(3)]]>

(3)式中，I为速度的时间采样点数，为初始拾取的均方根速度；通过求解所述目标函数求得变量ΔV和k；

第五步，将第四步求得的ΔV和k代入公式(1)，得到新的层速度，同时判断目标函数值是否小于阀值ε，如果是，则所述新的层速度即为最终的层速度数据体，然后转入第六步，否则，使n＝n+1，然后返回第二步；

第六步，输出所述层速度数据体。