利用瞬时谱属性进行薄互层分析.pdfVIP

利用瞬时谱属性进行薄互层分析 高静怀1 陈凤1 陈树民2 1． 西安交通大学电子与信息工程学院波动与信息研究所． 西安．710049 2． 大庆石油有限责任公司勘探开发研究院 黑龙江省大庆市．16371 1引言 在陆相油气勘探中，薄互层油气田研究占有十分重要的地位。例如我国大庆油田的中浅 层油气勘探、长庆油田的油气勘探等。常用的分析薄互层油气田地震资料的方法有：加窗 Fourier变换方法、小波变换方法等¨’zJ，这些方法是把薄互层地震记录变换到时一频域，然 后利用其能量在时一频域的分布规律，研究薄互层的结构，如：判断薄互层组是属于由厚逐 渐变薄型，还是属于其他类型llt2】，而无法研究薄互层内部的结构及岩性变化。 瞬时谱的属性(简称瞬时属性)是研究地下岩性分布及地层结构的有效工具【3，引，不同 的瞬时属性，反映出信号中所含的不同信息(见本文第二部分)。Barnes等人提出用瞬时带 宽、瞬时主频、瞬时品质因数等参数，与常用的瞬时频率、瞬时相位及瞬时振幅等参数相互 配合，分析反射地震资料的思想和方法瞄J。 然而，信号的瞬时属性对噪声是很敏感的【5l， 这使其在实际应用中受到限制。近年来发展起来的相空间瞬时属性计算方法，计算精度高， 抗噪能力强，在实际中见到了很好的效果¨．6t7．引。本文拟探讨瞬时属性用于薄互层分析的可 能性。第二部分给出了本文用到的瞬时属性的定义及算法；第三部分通过几个典型模型，研 究了不同瞬时属性对薄互层的响应；第四部分是实际资料算例；第五部分是结论。 2瞬时属性的定义及算法 设输入信号为SrO)，经HJlbcrt变换后结果为墨(r)，S，(f)对应的解析信号为 s(t)=Sr(f)+弘f(f)， 这里，J=4=-1。在实际应用中，为了克服分母趋于零时引起的不稳定性，常常采用阻尼瞬 时属性。阻尼瞬时频率、瞬时带宽、瞬时主频以及瞬时品质因数定义为【5】： 删2瓦Sl(丽t)Sr(≮t)-S丽r(t)S丽i(t) (2) k(f)s，(f)+s；(f)sf(f)I (3) m—ax上) 仃f(f)=』2x—(—Ai—(t—)2—+—F-A—i—(t)—2 fd(t)-厄丽葡万 (4) “垆研彘 (5) 上式中，￡o称为阻尼因子，下标max表示求极大值，如，g(f)一表示求g(f)的极大值， ，’(f)表示，(f)对时间的导数，J，(f)J表示，(f)的绝对值。下文算例中用到的瞬时频率、瞬时 带宽、瞬时主频以及瞬时品质因数均采用阻尼瞬时参数。上述瞬时属性与时一频谱的联系参 见文献‘21，其物理意义参见文献‘91．在实际资料分析中，尤其对低信噪比资料，可采用高静怀 等提出的相空间算法‘5。钔。本文不考虑噪声的影响，故采用Hilbert变换方法。 薄互层模型的瞬时属性分析 图1．(a)是反射介质模型：在500ms处反射 ¨ 系数为O．5；2-㈣是合成记录；1-(c)是该 记录的瞬时振幅；1-(d)是它的瞬时频率； 宝三圜望函 1-(e)是它的瞬时带宽；1-(f)是它的瞬时主 频；1-(g)是它的瞬时品质因数。 靼嘲望困：四 图1 图1是分析其它较复杂模型的基础。比较图l、图2：图2的瞬时频率、瞬时主频均有 两个峰，而图1的瞬时频率、瞬时主频均只有—个峰；图1的瞬时品质因数仅有—个峰值而 图2的瞬时品质因数有三个峰值，且中间的峰值较两边的峰值大很多。其他瞬时属性无明显 差异。 图2．(a)是反射介质模型：500ms和510ms 处，反射系数均为0．5；图2．(b)是合成记 录，图2-(c)是合成记录的瞬时振幅。图 2-(d)是它的瞬