一种地震波吸收衰减补偿方法技术

本发明专利技术涉及地球物理勘探数据处理的地震波吸收衰减补偿方法，在VSP原始数据选择一道作为参考道，进行广义S变换，选择采样点作为参考点，并记录时间及其时频谱，对每个样点数据进行广义S变换，将样点时频谱除以参考点时频谱对应频率的时频谱得到吸收衰减曲线，利用吸收衰减信号在时频域的频率差异公式，取自然对数算出相位补偿算子，进行时频域吸收衰减补偿和相位补偿，得到精细的VSP波场剖面。本发明专利技术对吸收衰减曲线拟合时可消除噪音，实现无能量损失的广义S反变换，补偿后可提高资料的信噪比及分辨率，实现记录分辨率控制，参考点选择灵活，算法简单，效率高。

【技术实现步骤摘要】

本专利技术涉及地震勘探数据处理技术，是一种利用地震数据的瞬时频率属性进行地 震波吸收衰减补偿方法。
技术介绍
在地震勘探中，地层对地震波的吸收造成地震子波时变、资料分辨率降低。特别是 在深层，地震资料不能反映地层的真实情况，使处理解释碰到很多难题。因此，有关吸收衰 减是需要解决的重点问题。当前吸收衰减补偿主要有以下几种方法1)常规反Q滤波。根据Futterman模型，振幅衰减满足权利要求1.，其特征在于具体步骤包括1)采用人工激发并由检波器记录地震波，取得原始地震数据，输入VSP原始数据，输入 和初至并按初至拉平；2)选择一道振幅较大、信噪比、主频较高资料作为参考道，对其进行广义S变换，在其 时频谱图上，选择瞬时频率较高，频带较宽的采样点作为参考点，并记录下该采样点的时间 及其时频谱；3)对每个样点数据进行广义S变换，并将该样点的时频谱除以参考点时频谱中对应频 率的时频谱得到吸收衰减曲线；4)利用吸收衰减信号在时频域的频率差异公式，对吸收衰减曲线取自然对数后利用最 小二乘法进行直线拟合，求得等效反射系及等效Q值，并计算出相位补偿算子；5)利用求取的等效反射系数和相位补偿算子，对该道数据中的每个采样点信号进行时 频域吸收衰减补偿和相位补偿；6)对该道数据补偿后的时频谱进行无能量损失的广义S反变换，变换到时间域，初至 反拉平后按道输出结果；7)重复步骤幻-6)，直到所有道都补偿完毕，经过图形处理得到准确描述精细的构造 特征VSP波场剖面。2.根据权利要求1所述的，其特征在于步骤幻所述的参 考道和参考点选择是地震记录中的某一时刻的信号，或某分辨率更高的子波，甚至是脉冲 信号中的脉冲点。3.根据权利要求1所述的，其特征在于步骤4)所述的 等效反射系数是直线拟合中的常数项，为无衰减的地震子波与反射系数之乘积；等效反射 系数对于一个地层反射系数，子波有多少个采样点，等效反射系数就有多少个采样点；等效 Q值是直线拟合中的一次项系数的倒数。4.根据权利要求1所述的，其特征在于步骤6)所述的变 换到时间域，若原始信号为零，则补偿后的数据依然为零；初至反拉平是初至前数据补零。5.根据权利要求1所述的，其特征在于步骤5)所述的吸 收衰减补偿采用如下公式计算PiTk，/) = ^p(T0J)oxPcxPr0Qeq(Tk)其中，P(Tk，f)为Tk时间，对应频率f的地震信号的时频谱；ρ(τ0, f)为参考点Ttl时间，对应频率f的地震信号的时频谱；Qeq{Tk) = l· ^ dT，Qeq (Tk)是 Tk 处的等效 Q 值； ^Q(T)rk为在没有吸收衰减情况下，Tk时间地层对应的反射系数；r0为在没有吸收衰减情况下，T0时间地层对应的反射系数；k P(ToJ) r0Qeq(Tk) }相位补偿算子为全文摘要本专利技术涉及地球物理勘探数据处理的地震波吸收衰减补偿方法，在VSP原始数据选择一道作为参考道，进行广义S变换，选择采样点作为参考点，并记录时间及其时频谱，对每个样点数据进行广义S变换，将样点时频谱除以参考点时频谱对应频率的时频谱得到吸收衰减曲线，利用吸收衰减信号在时频域的频率差异公式，取自然对数算出相位补偿算子，进行时频域吸收衰减补偿和相位补偿，得到精细的VSP波场剖面。本专利技术对吸收衰减曲线拟合时可消除噪音，实现无能量损失的广义S反变换，补偿后可提高资料的信噪比及分辨率，实现记录分辨率控制，参考点选择灵活，算法简单，效率高。文档编号G01V1/30GK102109612SQ20091024375公开日2011年6月29日 申请日期2009年12月23日 优先权日2009年12月23日专利技术者张固澜, 李彦鹏 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司本文档来自技高网...

【技术保护点】
１．一种地震波吸收衰减补偿方法，其特征在于具体步骤包括：１）采用人工激发并由检波器记录地震波，取得原始地震数据，输入ＶＳＰ原始数据，输入和初至并按初至拉平；２）选择一道振幅较大、信噪比、主频较高资料作为参考道，对其进行广义Ｓ变换，在其时频谱图上，选择瞬时频率较高，频带较宽的采样点作为参考点，并记录下该采样点的时间及其时频谱；３）对每个样点数据进行广义Ｓ变换，并将该样点的时频谱除以参考点时频谱中对应频率的时频谱得到吸收衰减曲线；４）利用吸收衰减信号在时频域的频率差异公式，对吸收衰减曲线取自然对数后利用最小二乘法进行直线拟合，求得等效反射系及等效Ｑ值，并计算出相位补偿算子；５）利用求取的等效反射系数和相位补偿算子，对该道数据中的每个采样点信号进行时频域吸收衰减补偿和相位补偿；６）对该道数据补偿后的时频谱进行无能量损失的广义Ｓ反变换，变换到时间域，初至反拉平后按道输出结果；７）重复步骤２）－６），直到所有道都补偿完毕，经过图形处理得到准确描述精细的构造特征ＶＳＰ波场剖面。

【技术特征摘要】

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