专利名称:一种基于Morlet谱勘测油气的方法
技术领域:
本发明涉 及能源勘探开发领域,特别是关于一种基于Morlet谱勘测油气的方法。
背景技术:
在油气地震勘探过程中,地震波在粘弹性介质中传播时,地震波的弹性能量逐渐被介质所吸收,转化为热能,弹性能转换成热能的过程称之为吸收。而含气地层能引起地震波高频能量的强烈吸收,通过计算它的衰减系数可以用来检测油气的存在与否。传统吸收系数计算都是利用窗口傅立叶变换求振幅谱,用高低频振幅谱比计算地震信号的高频衰减,即吸收系数。用傅立叶分析这个工具,虽能观察到高频能量的衰减现象,但难以描绘能量的细节变化。如要计算某一点的吸收系数,一般先要计算以该点为中心的某个时窗的傅立叶变换,但基于时窗的傅立叶变换计算衰减存在着如下问题若时窗太大,计算得到的频谱在时窗内有平均效应,频谱的时间局部性很差,从而影响吸收系数的精度,难以准确刻画出高频衰减的砂层位置,即所求得的结果不能满足地震记录的时间局部性要求;反之,若时窗太小,或者说期望有更高的时间分辨率,则由傅立叶变换的性质,会产生严重的边界效应,即吉布斯效应,所求得的频谱由于受到边界的污染,会有很多假象,影响了吸收系数的精度,也不能很好地反映地震信号的真实情况。
发明内容
针对上述问题,本发明的目的是提供一种能够描绘能量的细节变化,勘测结果准确、稳定,不存在人为异常的基于Morlet谱勘测油气的方法。为实现上述目的,本发明采取以下技术方案一种基于Morlet谱勘测油气的方法,其包括以下步骤1)对地震记录道目的层中的地震信号x(t)进行傅里叶分析,得到地震信号x(t)的频谱范围[f\,fh];根据频谱范围,确定进行分析的尺度个数;2)利用连续小波变换对地震信号x(t)进行分析,得到小波变换系数Wx(b,s)
权利要求
1. 一种基于Morlet谱勘测油气的方法,其包括以下步骤1)对地震记录道目的层中的地震信号x(t)进行傅里叶分析,得到地震信号x(t)的频谱范围[f\,fh];根据频谱范围,确定进行分析的尺度个数;2)利用连续小波变换对地震信号x(t)进行分析,得到小波变换系数Wx(b,s)
2.如权利要求1所述的一种基于Morlet谱勘测油气的方法,其特征在于所述步骤4) 中,Fl的具体计算方法为由所述步骤幻得到所有时间点、的点谱极大值所对应的尺度参数%,令
3.如权利要求1所述的一种基于Morlet谱勘测油气的方法,其特征在于所述步骤4) 中,&的具体计算方法为&为点谱在区间[Clpd2]上的振幅积分,其中,
4.如权利要求1所述的一种基于Morlet谱勘测油气的方法,其特征在于所述步骤4) 中,&的具体计算方法为&为点谱在区间[%,ae]上的振幅积分,其中,尺度参数 根据频谱形状确定, >%, ae在附近取值。
5.如权利要求1或2或3或4所述的一种基于Morlet谱勘测油气的方法, 其特征在于所述步骤幻中,小波函数Φ (t)满足
全文摘要
本发明涉及一种基于Morlet谱勘测油气的方法,其包括以下步骤1)对地震记录道目的层中的地震信号进行傅里叶分析,得到地震信号的频谱范围;2)利用连续小波变换对地震信号进行分析,得到小波变换系数Wx(b,s);3)定义不同尺度因子s下的瞬时振幅为,对于每一时间点b0,得到不同尺度因子ai所对应的瞬时振幅值{Ai(b0,1/ai)i=0,1,...,N},利用Lagrange样条进行插值,b点处任一频率的振幅值都可由对应的a值代入得到,这些离散点的连线,就构成了点b的频谱,即Morlet谱;4)根据步骤3)中得到的任一时间点处的瞬时振幅,得到每一点的吸收系数;如果该点的吸收系数ξI相比于围岩的吸收系数低百分之十以上,则表示该位置含有气层。
文档编号G01V1/30GK102169188SQ20101060615
公开日2011年8月31日 申请日期2010年12月15日 优先权日2010年12月15日
发明者刘兰峰, 曹思远, 李绪宣, 汪小将, 袁殿, 陈宝书, 高云峰 申请人:中国海洋石油总公司, 中国石油大学(北京), 中海石油研究中心