基于信号低频强异常自适应提取的地震资料流体识别方法

摘要

本发明提供了一种基于信号低频强异常自适应提取的地震资料流体识别方法，属于石油天然气地球物理勘探领域。本方法包括:步骤1：查看钻井岩心资料确定含流体储层的地下深度位置，并对井旁地震道的地震信号做谱分解，然后观察含流体储层相比围岩是否具有低频强振幅的地震反射特征，如果具有该特征，则转入步骤2;步骤2:对一道地震信号做基于低频强振幅异常的谱分解得到该道地震信号的时频谱;步骤3:获取步骤2得到的时频谱中每个时刻的瞬时最大能量;步骤4:对下一道地震信号重复步骤2和步骤3，直到所有地震道部完成为止;步骤5:利用所述低频强振幅的地震反射特征进行地震资料流体识别。利用本发明能快速自适应提取地震信号的低频异常。

主权项

基于信号低频强异常自适应提取的地震资料流体识别方法，其特征在于：所述方法包括：步骤1：查看钻井岩心资料确定含流体储层的地下深度位置，并对井旁地震道的地震信号做谱分解，然后观察含流体储层相比围岩是否具有低频强振幅的地震反射特征，如果具有该特征，则转入步骤2；步骤2：对一道地震信号做基于低频强振幅异常的谱分解得到该道地震信号的时频谱；步骤3：获取步骤2得到的时频谱中每个时刻的瞬时最大能量；步骤4：对下一道地震信号重复步骤2和步骤3，直到所有地震道都完成为止；步骤5：利用所述低频强振幅的地震反射特征进行地震资料流体识别，其中，所述步骤2中的做基于低频强振幅异常的谱分解得到该道地震信号的时频谱是利用下面的公式实现的：$HfWT(a,b)=\underset{-\infty }{\overset{+\infty }{\int }}f\left(t\right)w^{*}\left(\frac{t-b}{a}\right)adt---\left(2\right)$其中，f(t)为地震信号，WT为地震信号的小波变换，w^*(t)为子波的共轭复数，b为时间延迟，a为小波尺度因；设(2)式分解得到的单尺度信号为：式中j为虚数单位，为尺度a_i下分解的单频信号的初相位，对(3)式的结果取模得：|HfWT(a_i,b)|＝Aa_i      (4)尺度与频率的关系为$f_i=\frac{f_0}{{\mathrm{\Delta ta}}_i}$ f₀为小波变换中子波的主频，Δt为信号的时间采样间隔，公式(4)的结果为：$\left|HfWT(a_i,b)\right|={\mathrm{Aa}}_i=A\frac{f_0}{f_i\Delta t},---\left(5\right)$记为常数c则 |HfWT(a_i,b)|＝Ac/f_i        (6)。