[发明专利]一种时频分解地震流体识别方法

摘要

本发明涉及石油勘探技术领域，具体而言涉及一种时频分解地震流体识别方法，根据m(t)＝exp[-β·f2(t-τ)2]exp[i(2πf(t-τ)+φ)]的morlet小波函数建立时频原子库D，由地震道和复地震道方法计算得到morlet小波函数的初始匹配原子；针对地震道匹配分解，在初始匹配原子的邻域内以时频原子库D为约束进行迭代优选最佳的匹配原子，当达到预定的停止条件时停止匹配分解，这样就可以把原始地震道表示为一系列的morlet小波原子的线性组合；将这些最优匹配的morlet原子变换到时频域，从而可获取原始地震道的时频谱分布；在地震资料的时频谱上直接提取目的层层段的地震流体活动性属性；根据流体活动性属性预测气藏的分布范围和空间展布。本发明实施例可以准确的预测气藏的分布范围和空间展布，为天然气勘探的有利目标优选提供技术保障。

主权项

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