[发明专利]一种变偏移距VSP资料时差校正方法

权利要求书

1.一种变偏移距VSP资料时差校正方法，其特征在于：所述方法包括：

(1)解编采集到的原始变偏移距VSP资料，得到三分量波场数据；

(2)输入VSP观测井所在区域地面地震资料处理所用的浅层替换速度V_替换和地震资料处理基准面高程Z_基准面；

(3)在各偏移距垂直分量(Z)上分别拾取直达波初至时间得到时间-深度数据T₀(n,L_k)；

(4)计算单一偏移距内时差；

(5)计算变偏移距间时差；

(6)总时差校正：计算变偏移距各观测接收深度点的累计时差，并将得到的累计时差值校正到相应偏移距的接收深度点的初至时间上，得到经过时差校正后的变偏移距VSP三分量资料；

(7)校正误差分析：

在所述经过时差校正后的变偏移距VSP三分量资料上，抽取任一共深度点的垂直分量道集，判断直达波的初至时间曲线随偏移距的变化是否符合时差变化规律，如果是，则进入步骤(8)，如果否，则检查并重新输入误差点所在偏移距的激发炮井的高程、井深、坐标数据这三个数据，然后返回步骤(4)；

(8)结束；

其中，步骤(1)中的所述三分量波场数据包括VSP井坐标(x_w,y_w,z_w)、井台高h_W、各偏移距激发炮点坐标和激发井深为这四个数据，

其中，步骤(4)是这样实现的：

首先，选定激发炮群的中心点，以该中心点的井口高程为偏移距的统一高程面；然后，将炮群内的激发炮点校正到激发中心点的井口高程面；

分别计算出激发炮点到达相应接收深度点直达波的旅程D₁及激发中心点的井口高程面到达相应接收深度点直达波的旅程D₂，再计算出旅程差D₁-D₂，进一步计算出波到达同一接收点的旅程差产生的时差

其中，

$D_1(n_i,L_{k_i})=\sqrt{{(x_k-x_{n_i})}^2+{(y_k-y_{n_i})}^2+{\left(\right(L_{k_i}-h_w)-(z_w-z_{n_i}+h_{n_i}\left)\right)}^2}$

$D_2(n_i,L_{k_i})=\sqrt{{(x_k-x_{n_0})}^2+{(y_k-y_{n_0})}^2+{\left(\right(L_{k_i}-h_w)-(z_w-z_{n_0}\left)\right)}^2}$

x_w,y_w,z_w为观测井井口地面坐标；为第n个偏移距，第i个激发炮点的坐标；为激发炮群中心点的井口坐标；

以井口为0开始起算的电缆长度，计为L_k(k＝1,……N),N为VSP观测设计的最大观测点，井中各观测深度点的坐标可计为(x_k,y_k,L_k-h_w)；

n为偏移距序列，为第n个偏移距，第i个激发炮点到达相应接收深度点直达波的旅程，为第n个偏移距激发中心点井口平面到达相应接收深度点直达波的旅程；

当观测井为直井时，各观测深度点x_k＝x_w，y_k＝y_w；当观测井为斜井时，各观测深度点的x_k，y_k通过井斜数据计算得到。