一种利用地震数据瞬时频率属性进行油气检测的方法技术

本发明专利技术目的是利用地震数据瞬时频率属性进行油气检测的方法，包括采集记录地震原始数据；对数据在频率域进行广义Ｓ变换得到其每个样点的时频谱，确定时频谱中总能量的１５％和３５％分别对应的频率为计算低频段吸收衰减梯度的起始频率和终止频率；对起始频率和终止频率间的时频谱ｐ（Ｔｋ，ｆ）取自然对数后利用最小二乘法进行拟合，得到该样点的低频段吸收衰减梯度；低频段吸收衰减梯度强的位置对应含油气有利区域，确定油气检测结果。本发明专利技术对信噪比的要求不高，时频分辨能力强，算法简单，计算效率高。

【技术实现步骤摘要】

本专利技术涉及地震勘探数据处理技术，是一种利用地震数据瞬时频率属性进行油气 检测的方法。
技术介绍
含油气层的内部及其下部，地震波的能量将发生明显的高频衰减，利用地震能量 衰减进行油气预测成为可能。衰减是地震波在地下介质传播中总能量的损失，是介质内在的属性。引起地震波 衰减的因素有内部和外部两种因素引起地震波衰减的内部因素是介质中固体与固体、固 体与液体、液体与液体界面之间的能量耗损；外部因素主要来自大尺度的不均勻性介质引 起的散射。当前对于吸收衰减梯度的计算主要利用短时傅里叶变换和小波变换对地震波的 高频段进行分析计算。都是先对地震波进行时频分解，然后利用总能量的65%和85%之 间对应的时频谱进行曲线拟合，其得到的是高频段吸收衰减梯度。由于地震波传播过程中 要经受地震的强烈吸收，高频能量损失很快，在一定的信噪比情况下，当高频噪音干扰较强 时，高频段吸收衰减梯度更多的是反映噪音的影响，不能完全指示油气。同时，基于短时傅里叶变换的吸收衰减梯度计算方法存在两个严重不足(1)短 时傅里叶变换的时窗固定，时频网格等宽，不能随频率的需要而改变，对高低频的分辨都不 利。(2)当时窗选取过大，尤其是时窗内包含强反射层时，气层的衰减异常特征会减弱以至 消失；因为时窗选取过大，气层的衰减异特征受强反射层的影响而产生“屏蔽”效应。实际资料处理中的小波变换都是离散小波变换，其尺度因子与频率之间的关系不 确定，也会造成时频分辨能力不足的问题。广义S变换采用时窗宽度随频率f呈反比变化高斯窗函数，在低频段时窗较宽，获 得较高的频率分辨率；高频段时窗较窄，获得很高的时间分辨率。可以根据实际资料的分析 需要，调整参数，完全满足实际资料的时频分析。
技术实现思路
本专利技术目的是提供一种利用广义S变换对地震资料进行时频分析，并对低频段的 时频谱进行拟合，得到低频段吸收衰减梯度，以更好的指示油气的利用地震数据瞬时频率 属性进行油气检测的方法。本专利技术具体步骤包括1)采集记录地震原始数据；2)对数据在频率域进行广义S变换得到其每个样点的时频谱，确定时频谱中总能 量的15%和35%分别对应的频率为计算低频段吸收衰减梯度的起始频率和终止频率；3)对起始频率和终止频率间的时频谱ρ (Tk，f)取自然对数后利用最小二乘法进行 拟合，得到该样点的低频段吸收衰减梯度；3步骤3)所述的吸收衰减梯度的计算公式为 其中，P(Tk，f)为Tk时间，对应频率f的地震信号的时频谱； p(o, f)为地震子波的波峰处，对应频率f的地震信号的时频谱； Tk^ τ τ ’ Qeq (Tk)是Tk处的等效Q值；QeqD =I"Q(r)Q(T)为时间τ处对应的层Q s和r均为比例因子，取正数； fmain为地震子波的主频；4)重复步骤2)_3)，直到得到所有道中每个样点的低频段吸收衰减梯度；5)低频段吸收衰减梯度强的位置对应含油气有利区域，确定油气检测结果。 本专利技术对信噪比的要求不高，时频分辨能力强，对每个样点的时频谱进行曲线拟合时可消除部分较弱的噪音，同时避开高频干扰的影响；本专利技术算法简单，计算效率高。附图说明图1计算低频段吸收衰减梯度示意图；图2动校后的VSP上行波；图3低频段吸收衰减梯度图；图4瞬时主频；图5层速度；具体实施例方式本专利技术是在地震勘探资料中利用广义S变换，将地震资料变换到时频域；再结合 地震波吸收衰减在瞬时频率上的差异，在时频域对低频段的时频谱进行拟合，得到低频段 吸收衰减梯度，以更好的指示油气。以下结合附图详细说明本专利技术。本专利技术具体实现步骤如下1)采集记录地震原始数据；图2是某井的VSP资料经过处理后的上行波拉平剖 面，据钻井显示，含气层出现在1.78-1.79S处;2)对数据在频率域进行广义S变换得到其每个样点的时频谱，确定时频谱中总能 量的15%和35%分别对应的频率为计算低频段吸收衰减梯度的起始频率和终止频率；3)对起始频率和终止频率间的时频谱ρ (Tk，f)取自然对数后利用最小二乘法进行 拟合，计算该样点的低频段吸收衰减梯度；如利用图1中总能量的15% -35%之间对应的时 频谱进行拟合；步骤3)所述的吸收衰减梯度的计算公式为4 其中，p(Tk，f)为Tk时间，对应频率f的地震信号的时频谱；p(0，f)为地震子波的波峰处，对应频率f的地震信号的时频谱； ，Qeq (Tk)是 Tk 处的等效 Q 值； Α Q(T)Q(T)为时间τ处对应的层Qs和r均为比例因子，取正数；fmain为地震子波的主频；4)重复步骤2)_3)，得到所有道中每个样点的低频段吸收衰减梯度；5)确定低频段吸收衰减梯度强的位置，对应含油气有利区域，进行油气检测。如图3中在1. 78-1. 79S处，低频段吸收衰减梯度明显变大，反应出地层很强的吸 收，图4瞬时主频是时频谱中能量最大值对应的频率，在1. 78-1. 79S处，地震波的瞬 时主频约为20hz，相对围岩而言，明显降低。图5层速度在1. 78-1. 79S处，明显降低到4800m/s，变化趋势吻合很好，指示了油 气的存在。权利要求，其特征是具体步骤包括1)采集记录地震原始数据；2)对数据在频率域进行广义S变换得到其每个样点的时频谱，确定时频谱中总能量的15％和35％分别对应的频率为计算低频段吸收衰减梯度的起始频率和终止频率；3)对起始频率和终止频率间的时频谱p(Tk，f)取自然对数后利用最小二乘法进行拟合，得到该样点的低频段吸收衰减梯度；4)重复步骤2) 3)，直到得到所有道中每个样点的低频段吸收衰减梯度；5)低频段吸收衰减梯度强的位置对应含油气有利区域，确定油气检测结果。2.根据权利要求1的方法，特征是步骤3)所述的吸收衰减梯度的计算公式为 )其中，p(tk，f)为tk时间，对应频率f的地震信号的时频谱； p(o, f)为地震子波的波峰处，对应频率f的地震信号的时频谱； q(t)为时间τ处对应的层q s和r均为比例因子，取正数； fmain为地震子波的主频。全文摘要本专利技术目的是利用地震数据瞬时频率属性进行油气检测的方法，包括采集记录地震原始数据；对数据在频率域进行广义S变换得到其每个样点的时频谱，确定时频谱中总能量的15％和35％分别对应的频率为计算低频段吸收衰减梯度的起始频率和终止频率；对起始频率和终止频率间的时频谱p(Tk，f)取自然对数后利用最小二乘法进行拟合，得到该样点的低频段吸收衰减梯度；低频段吸收衰减梯度强的位置对应含油气有利区域，确定油气检测结果。本专利技术对信噪比的要求不高，时频分辨能力强，算法简单，计算效率高。文档编号G01V1/28GK101923176SQ20101022706公开日2010年12月22日 申请日期2010年7月6日 优先权日2010年7月6日专利技术者容娇君, 张固澜, 张彦斌, 彭继新, 李彦鹏 申请人:中国石油天然气集团公司;中国石油集团东方地球物理勘探有限责任公司本文档来自技高网...

【技术保护点】
一种利用地震数据瞬时频率属性进行油气检测的方法，其特征是具体步骤包括：１）采集记录地震原始数据；２）对数据在频率域进行广义Ｓ变换得到其每个样点的时频谱，确定时频谱中总能量的１５％和３５％分别对应的频率为计算低频段吸收衰减梯度的起始频率和终止频率；３）对起始频率和终止频率间的时频谱ｐ（Ｔ↓［ｋ］，ｆ）取自然对数后利用最小二乘法进行拟合，得到该样点的低频段吸收衰减梯度；４）重复步骤２）－３），直到得到所有道中每个样点的低频段吸收衰减梯度；５）低频段吸收衰减梯度强的位置对应含油气有利区域，确定油气检测结果。

【技术特征摘要】

【专利技术属性】
技术研发人员：张固澜，李彦鹏，彭继新，容娇君，张彦斌，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：11[中国|北京]