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

摘要

本发明提供了一种方位地震子波的融合方法，属于油气及煤层气地震勘探与开发领域。本发明方法包括以下步骤(1)方位子波预处理从同一个地震数据或CMP叠前数据中提取至少两个不同的方位子波；(2)提取各方位子波的主频和主频相位；(3)利用步骤(2)得到的各方位子波的主频和主频相位生成新子波。利用本发明方法能够使一个地震子波可以包含多个地震子波的特性。

主权项

一种方位地震子波的融合方法，其特征在于：所述方法包括以下步骤：(1)方位地震子波预处理：从同一个地震数据或CMP叠前数据中提取至少两个不同的方位地震子波；(2)提取各方位地震子波的主频和主频相位；(3)利用步骤(2)得到的各方位地震子波的主频和主频相位生成新子波，所述步骤(3)具体如下：首先，利用公式(1)和公式(2)将所有方位地震子波的主频与主频相位做平均：公式(1)中，average_freq为最终得到的平均主频，freq为每个方位地震子波的主频，n表示共有n个方位地震子波，k表示方位地震子波的序号；公式(2)中average_ph为最终得到的平均主频相位，ph为每个方位地震子波的主频相位；然后，利用公式(3)生成新子波：公式(3)中，t是采样时间，a是子波衰减系数，f是主频，将由公式(1)得到的average_freq作为f代入公式(3)，输入新子波的样点数，求得新子波b(t)；再利用公式(4)将新子波b(t)用傅里叶变换求其频谱：公式(4)中，|B(ω)|是子波振幅谱，是子波的相位谱，具体是这样得到的：|B(ω)|是子波振幅谱，其公式为并将B(ω)求模即得到|B(ω)|；是子波的相位谱，其中imag(B(w))为频谱的虚部，real(B(w))为频谱的实部；最后，将average_ph作为代入公式(4)，之后再利用傅里叶逆变换将B(ω)生成新子波b1(t)，此子波即融合后的新子波。