【技术实现步骤摘要】
一种基于最大熵魏格纳-威尔分布的含气性检测方法
[0001]本专利技术涉及石油、天然气勘探开发领域,特别是一种基于最大熵魏格纳-威尔分布的含气性检测方法。
技术介绍
[0002]地震信号是非平稳信号,其瞬时频率是随时间变化的,传统的瞬时频率计算是利用希尔伯特变换的复数道来实现。近年来,随着时频分析方法技术的出现,瞬时频率通过精度更高的时频分析方法来求取,对非平稳地震动信号的常用时频方法有短时傅里叶变换、小波变换和魏格纳-威尔分布、S变换,希尔伯特黄变换及一些改进的算法。这些方法都已广泛应用,但这些方法仍然存在一定的不足,短时傅里叶变换由于其时间和频率分辨率很低;小波变换存在小波变换的冗余度大、小波基的问题;魏格纳-威尔分布时频分布集中但在时频平面上则存在交叉项干扰问题,S变换分辨率较高但存在频率向高频漂移问题;希尔伯特黄变换出现,频率出现畸变、负值振荡等问题。
[0003]通常瞬时频率通过相位求导来获得,当信号存在噪声时,噪声对计算结果影响大。瞬时频率也可通过时频谱的一阶矩来求取,并且一阶矩方法是通过积分来实现的,对
【技术保护点】
【技术特征摘要】
1.一种基于最大熵魏格纳-威尔分布的含气性检测方法,其特征在于,包括:S100对地震信号,通过最大熵魏格纳-威尔分布的方法,得到时频谱;最大熵魏格纳-威尔分布具体包括:采用魏格纳-威尔时频分布,并利用MEM方法对魏格纳-威尔时频分布中核函数进行外推填补观察区间外的数据,以消除魏格纳-威尔时频分布的交叉项,得到时频谱;S200对时频谱,采用一阶矩的方法通过积分获取瞬时频率;S300在过井剖面上,将瞬时频率与设定的阈值进行判断:当储层的瞬时频率大于或等于阈值,则认为储层为产气井;当储层的瞬时频率小于阈值,则认为储层为干井。2.根据权利要求1所述的一种基于最大熵魏格纳-威尔分布的含气性检测方法,其特征在于,所述步骤S100中,所述魏格纳-威尔时频分布为:其中z(t)是地震信号x(t)的分析信号,t是时间,ω是频率,τ是延迟,其中zz*为核函数;由于地震信号是多分量数据,考虑魏格纳-威尔时频分布的二次性,会受到干扰项的限制;用单频率分量z(t)和g(t)来描述,魏格纳-威尔时频分布表示如下:W
z+s
=W
z
(t,f)+W
s
(t,f)+2Re[W
z,s
(t,f)]其中W
z
(t,f)和W
s
(t,f)分别是分量z(t)和g(t)的魏格纳-威尔时频分布,Re[W
z,s
【专利技术属性】
技术研发人员:许多,董霞,喻勤,王浩,钱玉贵,杨国鹏,刘红爱,王鹏,范红娟,高顺华,
申请(专利权)人:中国石油化工股份有限公司西南油气分公司,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。