基于非稳态反褶积的变Q值求算方法技术

本发明专利技术的一种基于非稳态反褶积的变Q值求算方法是，在对数域中充分利用并统计分析双曲条带内的能量关系，通过利用线性拟合的斜率信息建立了Q求算的基本技术思路，并且借助迭代机制不断压制地震子波的不利影响从而增强Q值估算的纵横向稳定性及实用性。本发明专利技术首次实现了利用能量统计关系求算Q值，与传统方法相比，能更充分地利用整个时频域的数据，信息量充足，估算结果更可靠，与地质背景吻合情况更好。强调利用迭代线性拟合的思路不断压制地震子波影响，进而实现对地下Q值的精度逼近。

【技术实现步骤摘要】
基于非稳态反褶积的变Q值求算方法
本专利技术涉及属于地球物理勘探领域，尤其涉及一种基于非稳态反褶积的变Q值求算方法。
技术介绍
非稳态特征是地震资料处理领域中的一个重要研究方面。现行的大多数反褶积均建立在经典的Robinson褶积模型的基础之上，然而此褶积模型近似地认为地震子波不随时间变化，对地震信号固有的非稳态性特征认识不足，因此这类反褶积又称作静态反褶积，其实用价值在很大程度上受到限制。不同于静态反褶积技术，Gabor反褶积以非稳态褶积模型为基础，重点考虑了吸收衰减影响，通过同时移除吸收衰减与震源子波(即传播子波)的方式移除地震信号的非稳态特征，以此Gabor反褶积又称作非稳态反褶积。在过去十年中，大量学者与研究人员对Gabor反褶积的诸多技术方面进行了改善，增强其稳定性和可靠性，内容包括时窗函数，时频谱估计，相位校正以及更合理的平滑方法。这些技术层面也巩固了 Gabor反裙积在地震资料处理领域中的地位和领域。 与传统静态反褶积不同，Gabor反褶积能够解决由吸收衰减所造成的非稳态特征。根据使用最为广泛的Q理论，吸收衰减表现为随时间和频率增加，振幅呈指数衰减 a (t, f)=e^tf/(MHUtf/Q) (I) 其中， t，f分别代表时间与频率； Q代表地层品质因子； H 为 Hi Ibert 变换。 一般认为，吸收衰减具有最小相位性质，所以上式中的虚部提供最小相位信息，H则代表Hilbert变换。实质上，非平稳地震道s (t)是在吸收衰减a (t, f)与反射系数r(t) 动态褶积的基础上，再与震源子波>.?’(/)进行静态褶积的结果，即本文档来自技高网...

【技术保护点】
一种基于非稳态反褶积的变Q值求算方法，其特征在于：在对数域中分析双曲条带内的能量统计关系，建立线性拟合方程，从获得的斜率信息实现对Q值的求算，并且利用去野值、插值、平滑及直方图校正等手段增强Q值估算的纵横向稳定性及实用性。

【技术特征摘要】
1.一种基于非稳态反褶积的变Q值求算方法，其特征在于:在对数域中分析双曲条带内的能量统计关系，建立线性拟合方程，从获得的斜率信息实现对Q值的求算，并且利用去野值、插值、平滑及直方图校正等手段增强Q值估算的纵横向稳定性及实用性。2.根据权利要求1所述的基于非稳态反褶积的变Q值求算方法，其特征在于；该方法步骤如下: 步骤1:在对未进行衰减补偿的地震道进行Gabor变换的基础上取对数，进入对数时频域； 步骤2:在对数时频域沿时间方向按实际情况确定时窗，在时窗划分双曲条带并统计其能量关系，形成衰减曲线； 步骤3:对衰减曲线进行的主要部分进行线性拟合，根据拟合斜率确定该分析时窗中...

【专利技术属性】
技术研发人员：孙学凯，孙赞东，韩剑发，
申请(专利权)人：孙赞东，
类型：发明
国别省市：北京;11