[发明专利]一种确定地震数据反射界面位置的方法

摘要

本发明提供了一种确定地震数据反射界面位置的方法，属于石油地球物理勘探领域。本方法包括：(1)设地震波为f(t)，对其进行Hilbert变换，求取实地震信号的解析信号，得到地震波f(t)的瞬时振幅、瞬时相位剖面θ(t)和瞬时频率；(2)基于地震波瞬时相位剖面θ(t)进行一阶差分运算，得到瞬时相位一阶差分G(t)；(3)对地震数据f(t)进行极性反转处理得到F(t)＝-f(t)，对F(t)进行Hilbert变换，得到瞬时相位剖面θ₁(t)，再对瞬时相位剖面θ₁(t)进行一阶差分运算，得到极性反转后的地震波瞬时相位一阶差分j(t)，对其求相反数得到J(t)＝-j(t)。

主权项

一种确定地震数据反射界面位置的方法，其特征在于：所述方法包括：(1)设地震波为f(t)，对其进行Hilbert变换，求取实地震信号的解析信号，得到地震波f(t)的瞬时振幅、瞬时相位剖面θ(t)和瞬时频率；(2)基于地震波瞬时相位剖面θ(t)进行一阶差分运算，得到瞬时相位一阶差分G(t)，则G(t)＝θ(t+1)‑θ(t)；(3)对地震数据f(t)进行极性反转处理得到F(t)＝‑f(t)，对F(t)进行Hilbert变换，得到瞬时相位剖面θ₁(t)，再对瞬时相位剖面θ₁(t)进行一阶差分运算，得到极性反转后的地震波瞬时相位一阶差分j(t)，对其求相反数得到J(t)＝‑j(t)；(4)将步骤(2)和步骤(3)的结果作为输入，利用下式对两个结果进行加权相加运算：M(t)＝a₁*G(t)+a₂*J(t)其中a₁+a₂＝1，a₁代表地震数据中波谷占的权重，a₂代表地震数据中波峰占的权重；(5)将步骤(4)的结果M(t)与原始地震波f(t)相乘得到新的属性Q(t)：Q(t)＝M(t)*f(t)，并选取阈值滤去新的属性Q(t)中的较低值，获得准确的反射界面位置。