一种利用斜井直达地震波走时计算地层地震波速度的方法

摘要

本发明公开了一种利用斜井直达地震波走时计算地层地震波速度的方法，该方法假设激发地震波的位置为炮点S，其二维坐标为(x_s,z_s)，沿钻井井轨迹有N个地震波接收点，这些接收点从上到下按从小到大的顺序依次记为R₁,R₂,R₃,...,R_N，相应的接收点坐标记为(x₁,z₁),(x₂,z₂),(x₃,z₃),...,(x_N,z_N)，各接收点所接收到直达地震波的实际走时依次记为t₁,t₂,t₃,...,t_N。本发明适用于斜井、并可估算出光滑的地层地震波速度。

主权项

一种利用斜井直达地震波走时计算地层地震波速度的方法，其特征在于，包括以下步骤：1)设激发地震波的位置为炮点S，其二维坐标为(x_s,z_s)，沿钻井井轨迹有N个地震波接收点，这些接收点从上到下按从小到大的顺序依次记为R₁,R₂,R₃,...,R_N，相应的接收点坐标记为(x₁,z₁),(x₂,z₂),(x₃,z₃),...,(x_N,z_N)，各接收点所接收到直达地震波的实际走时依次记为t₁,t₂,t₃,...,t_N；设地下有N个水平地层，从上到下按从小到大的顺序从1开始依次编号，各地层底界面纵坐标从上到下依次对应接收点的纵坐标z₁,z₂,z₃,...,z_N，各地层速度从上到下依次为v₁,v₂,v₃,...,v_N；2)按地震波以直线传播的假设从上到下逐层计算地层速度，计算公式为：$v_i=(t_i-\underset{j=1}{\overset{i-1}{\Sigma }}{L}_j·v_j)/L_i---\left(1\right)$其中，L_j表示从炮点S到接收点R_j的地震波传播路径在第j个地层里的长度，其计算公式为：3)对以上计算得到的层速度v_i，采用(2M+1)点滑动窗口平均进行光滑处理，其中i＝1,2,3,...,N，计算公式如下：$v_i=\left\{\begin{array}{cc}\underset{j=1}{\overset{1+2M}{\Sigma }}{L}_j·v_j/(2M+1),& 1\le i\le M\\ \underset{j=1i-M}{\overset{i+M}{\Sigma }}{L}_j·v_j/(2M+1),& M<i<N-M\\ \underset{j=N-2M}{\overset{N}{\Sigma }}{L}_j·v_j/(2M+1),& N-M\le i\le N\end{array}\right.---\left(3\right)$其中，M的大小反映计算结果的光滑程度，M的值在M<(N‑1)/2条件下任意选取；令迭代次数iter的值为1；4)根据地震波Snell定律计算在速度v_i下的理论直达地震波走时T_i，并计算理论走时T_i和实际走时t_i的均方根误差rms_error，其中i＝1,2,3,...N；理论直达地震波走时T_i的计算方法采用射线追踪法；5)设定一个正数eps以及一个最大迭代次数Nmax，例如eps＝0.01及Nmax＝100，若计算的理论走时T_i(i＝1,2,3,...,N)和实际走时t_i(i＝1,2,3,...,N)的均方差rms_error>eps，并且迭代次数iter<Nmax，则按下式更新速度：${\mathrm{dv}}_i=\frac{(i-1)(T_i-t_i)}{N\underset{j=1}{\overset{i}{\Sigma }}{L}_j},i=1,2,3,\mathrm{...},N---\left(10\right)$然后令迭代次数iter的值增1并返回步骤4)；否则，输出v_i(i＝1,2,3,...,N)作为最终计算结果，计算结束。