[发明专利]局部地震道内插方法和装置

权利要求书

1.一种局部地震道内插方法，其特征在于，包括：

获取距离待插值的地震道最近的且满足预定条件的三个已有地震道；

计算所述三个已有地震道中各个地震道的一阶时间导数；

将所述三个已有地震道中位于中间的地震道作为中心地震道，对所述中心地震道每个采样点计算同一时刻的所述三个已有地震道中除所述中心地震道之外的两个地震道采样点的加权差值和所述三个已有地震道的一阶时间导数的加权和；

对所述中心地震道的每个采样点，利用给定时窗内的所有的所述加权差值和所述加权和，计算该采样点对应的信号的瞬时道间时差；

计算所述中心地震道与所述待插值的地震道之间的空间距离；

根据所述中心地震道每个采样点对应的信号的瞬时道间时差和所述空间距离，将所述中心地震道的数值外推到所述待插值的地震道；

其中，所述预定条件为：

空间坐标满足h_k＞x_k＞y_k＞z_k且x_k-y_k＝y_k-z_k＝Δx，或者，x_k＞y_k＞z_k＞h_k且x_k-y_k＝y_k-z_k＝Δx

其中，h_k表示待插值的地震道的空间坐标，y_k表示中心地震道的空间坐标，x_k和z_k分别表示所述三个已有地震道中除所述中心地震道之外的其它两个地震道的空间坐标，k表示待插值的地震道的序号，Δx表示所述三个已有地震道中相邻地震道的道间距。

2.如权利要求1所述的方法，其特征在于，按照以下公式计算所述三个已有地震道中各个地震道的一阶时间导数：

$U_k(x_k,t_n)=\frac{B_k(x_k,t_{n+1})-B_k(x_k,t_{n-1})}{2\mathrm{\Δ}t}$

$V_k(y_k,t_n)=\frac{C_k(y_k,t_{n+1})-C_k(y_k,t_{n-1})}{2\mathrm{\Δ}t}$

$W_k(z_k,t_n)=\frac{D_k(z_k,t_{n+1})-D_k(z_k,t_{n-1})}{2\mathrm{\Δ}t}$

其中，V_k(y_k,t_n)表示中心地震道的一阶时间导数的值，U_k(x_k,t_n)和W_k(z_k,t_n)分别表示所述三个已有地震道中除所述中心地震道之外的其它两个地震道的一阶时间导数的值，C_k表示中心地震道的数值，B_k和D_k分别表示所述三个已有地震道中除所述中心地震道之外的其它两个地震道的数值，y_k表示中心地震道的空间坐标，x_k和z_k分别表示所述三个已有地震道中除所述中心地震道之外的其它两个地震道的空间坐标，n表示地震道的时间采样点的序号，t_n表示第n个采样点的采样时间，Δt表示地震道的时间采样间隔，k表示待插值的地震道的序号。