[发明专利]基于N次根叠加的定向地震波数据处理方法

权利要求书

1.一种基于N次根叠加的定向地震波数据处理方法，其特征在于，包括以下步骤：

A、输入n炮炮点连续的地震记录，记为U＝{U₁(t,x)、U₂(t,x)、…、U_m(t,x)、…、U_n-1(t,x)、U_n(t,x)}，t代表地震记录时间，x代表接收阵列，n取值为大于2的奇数，其中m＝(n+1)/2，m代表中间地震记录编号，x＝{x₁,x₂,…x_i,…,x_p}，i为接收阵列中的各道采集序号，i＝1,2,…,p，p为单炮记录的采集道数；

B、对U中n炮地震记录按照公式(1)

${U_j}^{\′}(t,x)=\frac{U_j(t,x)}{\left|U_j(t,x)\right|}{\[\left|U_j(t,x)\right|\]}^{\frac{1}{N}}\mathrm{...}\left(1\right)$

进行开N次根处理，其中j＝1,2,…,n，U_j'(t,x)为开N次根处理后的第j炮地震记录，经过开N次根处理后新地震记录记为U'＝{U₁'(t,x)、U₂'(t,x)、…、U_m'(t,x)、…、U_n-1'(t,x)、U_n'(t,x)}；

C、根据地震勘探原理及勘探区域目标层位信息，确定定向地震波主波束方向轴向角θ和目标层上覆介质平均速度v，按公式(2)

${\mathrm{\τ}}_0=\frac{d_{1}cos\left(\mathrm{\θ}\right)}{v}\mathrm{...}\left(2\right)$

计算τ₀，其中τ₀为延时参数，d₁为炮间距；

D、以第m炮地震记录为零延时基准，对U'中各地震记录依次分别按照延时参数τ₀进行等间隔延时，得到地震记录集

E、令λ＝d₁/d₂，其中d₂为道间距，根据地震勘探观测系统设计要求d₁为d₂的整数倍，因此λ为整数，在R的各炮地震记录中分别选取对应R_m(t,x)中的第x_i道接收位置的地震信号，并将选取的n炮地震信号按照公式(3)

$\begin{array}{c}r(t,x_i)=R_1(t-(\frac{n-1}{2}){\mathrm{\τ}}_0,x_i+(\frac{n-1}{2}\left)\mathrm{\λ}\right)+R_2(t-(\frac{n-1}{2}-1){\mathrm{\τ}}_0,x_i+(\frac{n-1}{2}-1\left)\mathrm{\λ}\right)+...+\\ R_m(t,x_i)+\mathrm{...}\end{array}$

$+R_{n-1}(t+(\frac{n-1}{2}-1){\mathrm{\τ}}_0,x_i-(\frac{n-1}{2}-1\left)\mathrm{\λ}\right)+R_n(t+(\frac{n-1}{2}){\mathrm{\τ}}_0,x_i-(\frac{n-1}{2}\left)\mathrm{\λ}\right)\mathrm{...}\left(3\right)$

进行叠加，其中r(t,x_i)为叠加后地震记录中第i道地震信号，对应i＝1,2,…,p，按公式(3)求出R_m(t,x)地震记录中1,2,...,p道地震信号，并定义新的地震记录M(t,x)＝[r(t,x₁),r(t,x₂),…,r(t,x_p)]；

F、令

W(t,x)＝M(t,x)|M(t,x)|^N-1………(4)

W(t,x)即为生成的能够抑制主波束外畸变信号的定向地震记录。