主权项 |
1.一种时频分解地震流体识别方法,其特征在于包括,步骤1,根据m(t)=exp[-β·f<sup>2</sup>(t-τ)<sup>2</sup>]exp[i(2πf(t-τ)+φ)]的morlet小波函数建立时频原子库D,其中τ为中心时间;f为主频;φ相位;β为能量衰减因子,所述β用于调节小波的时间延续度和衰减速度;所述时频原子库D表示为<img file="FDA00003025482400013.GIF" wi="470" he="92" /><img file="FDA00003025482400014.GIF" wi="71" he="63" />为小波原子,γ<sub>n</sub>={τ<sub>n</sub>,f<sub>n</sub>,φ<sub>n</sub>,β<sub>n</sub>},τ<sub>n</sub>、f<sub>n</sub>、φ<sub>n</sub>和β<sub>n</sub>分别为第n个小波原子的中心时间、主频、相位和衰减因子;步骤2,将输入的地震道的数据进行希尔伯特黄变换(HHT)得到复地震道,根据所述地震道和复地震道计算初始中心时间τ、主频f和相位φ,得到上述morlet小波函数的初始匹配原子,其中β在计算初始匹配原子时取值β=4ln2;步骤3,对所述初始匹配原子以所述时频原子库D为约束进行迭代匹配,得到最优匹配原子;步骤4,针对输入的地震道数据分时窗进行扫描进行上述步骤2-3,得到多个最优匹配原子,多个最优匹配原子的线性组合得到重构的地震道数据;步骤5,根据所述重构的地震道数据,与对应的地震道数据做差,得到残差;步骤6,当不满足预定的停止条件时,将该残差作为地震道数据重复步骤2-5,直到满足预定停止条件时停止,此时的最优匹配原子为最终匹配原子;步骤7,将所述最终匹配原子变换到时频域,获取地震道时频谱分布;步骤8,针对所有地震道数据进行上述步骤2-7;步骤9,从地震资料的时频谱上直接提取目的层层段的地震流体活动性属性;步骤10,根据所述流体活动性属性预测气藏的分布范围和空间展布。 |