[发明专利]一种方位地震子波的融合方法

权利要求书

1.一种方位地震子波的融合方法，其特征在于：所述方法包括以下步骤：

(1)方位子波预处理：从同一个地震数据或CMP叠前数据中提取至少两个不同的方位子波；

(2)提取各方位子波的主频和主频相位；

(3)利用步骤(2)得到的各方位子波的主频和主频相位生成新子波。

2.根据权利要求1所述的方位地震子波的融合方法，其特征在于：所述步骤(2)具体如下：

每个方位子波的主频是这样得到的：将方位子波转换到频率域后，其频谱最大值点对应的频率即为主频；

所述主频相位是描述讯号波形变化的度量，以度或角度作为单位。

3.根据权利要求2所述的方位地震子波的融合方法，其特征在于：所述步骤(3)具体如下：

首先，利用公式(1)和公式(2)将所有方位子波的主频与主频相位做平均：

average_freq=1nΣk=1nfreqk]]>公式(1)

average_ph=1nΣk=1nphk]]>公式(2)

公式(1)中，average_freq为最终得到的平均主频，freq为每个方位子波的主频，n表示共有n个方位子波，k表示方位子波的序号；

公式(2)中average_ph为最终得到的平均主频相位，ph为每个方位子波的主频相位；

然后，利用公式(3)生成新子波：

b(t)＝e^-at2Sin2πft

                           公式(3)

公式(3)中，t是采样时间，b是新子波的样点数，a是子波衰减系数，f是主频，将由公式(1)得到的averagefreq作为f代入公式(3)，求得子波b(t)；

再利用公式(4)将子波b(t)用傅里叶变换求其频谱：

公式(4)

公式(4)中，|B(ω)|是子波振幅谱，是子波的相位谱，具体是这样得到的：

|B(ω)|是子波振幅谱，其公式为并将B(ω)求模即得到|B(ω)|；

是子波的相位谱，其中imag(B(w))为频谱的虚部，real(B(w))为频谱的实部；

最后，将average_ph作为代入公式(4)，之后再利用傅里叶逆变换将B(ω)生成新子波b1(t)，此子波即融合后的新子波。

4.根据权利要求1所述的方位地震子波的融合方法，其特征在于：所述公式(3)中，采样时间t、新子波的样点数b和是子波衰减系数a都是由用户输入的。