[发明专利]由叠前时间偏移直接扫描确定地层叠加品质因子方法

说明书

技术领域

本发明属于油气田和矿产资源地震勘探中反射地震资料处理技术领域，是获取高分辨率叠前地震成像重要输入参数场的一种技术手段。

背景技术

叠加Q值是时域反Q滤波、粘性介质叠前时间偏移所需的一个重要参数，精确反演叠加Q值对现场的高分辨地震数据处理具有重要的应用价值。

Q值反演方法很多，包括时域算法和频域算法。时域算法Q值反演得精度较差，而频域算法中最具代表性的两种方法是：对数谱比法和频率移动法，二者都是建立在粘性吸收是地震反射波谱唯一的影响因素假设基础上的。频移法通常只使用峰值或平均频率等一种频率信息，因此其稳定性和有效性均受到质疑；对数谱比法是进行Q值估计的一种经典的方法，因为其与源的特征无关，可以消除一切与频率无关的影响因素；但是对数谱比法对于噪声和谱估计误差十分敏感，通常也会受到能量泄露的影响，使得算法并不稳定。

从数据来源看，Q值可由地震透射数据，如VSP(垂直地震剖面)资料、井间地震数据、声波测井数据求取。但这些资料都要受到观测井位置和数量的限制，无法给出整个工区的Q值场在时间和空间的变化情况。

从实际应用的角度来讲，由地面地震资料来反演Q值场更加具有意义。Dasgupta and Clark首次提出了基于对数谱比法由共中心点道集(CMP)地面数据反演Q值场；Zhang and Ulrych采用频移法反演CMP道集的Q值场，取得了一定的效果。

现有的由地面地震数据反演Q值场的方法存在诸多问题，首先，叠加后的共中心点道集是不同偏移距数据信号(其代表了不同的传播路径，传播距离、入射角和不同炮-检对)的一种折中效应，叠加后将导致叠加道出现不可逆转的高频损失；对于水平层状地层，其可以反映地下同一深度点的反射情况。但是当地下情况不满足这一假设时，则不同偏移距的数据是来自不同反射点的反射衰减数据，叠加后的反射波谱并没有真正反映地下同一反射点的真实衰减情况。

其次，对于大多数沉积地层，薄层的调谐作用和短周期的层间多次波对反射波谱具有重要影响。当地层厚度小于1/4波长时，会由于建设性干涉作用和破坏性干涉作用而使振幅谱的幅值增大或者减小。除了幅度变化以外，反射地震波的视频率和相位也要发生变化，在此种情况下，基于频谱变化的对数谱比法和频移法在Q值估计上都会存在很大的误差。

第三，相比于透射数据资料(如VSP资料)，反射数据具有更长的传播路径而使幅度进一步衰减，其数据更易受到噪声的影响。

第四，精确的信号谱估计是提高Q值反演精度的关键因素。以往地震波谱获取是通过短时傅立叶变换来实现的，但实际地震数据是非平稳信号，使得该方法所获取的频谱与真实的子波谱相比产生很大扭曲。

对于由地面地震数据反演叠加Q值场所存在的问题，需要一种由叠前时间偏移直接扫描确定反演叠加Q值场的新方法，该方法需要克服现有叠加Q值分析方法的缺点，可以给出具有较高精度的、具有时变和空变特性的、适用于粘性介质叠前时间地震偏移成像处理的叠加Q值场。

发明内容

本发明的目的是提供一种由叠前时间偏移直接扫描确定地层叠加品质因子方法。该方法通过采用地面叠前反射数据反演分析叠加Q值，可以克服由VSP资料、井间地震数据、声波测井等资料无法给出整个工区的时变和空变Q值场的缺陷；通过采用补偿成像后的共反射点道集，使不同偏移距数据能够反映地下同一反射点的黏性吸收补偿情况；采用广义S变换精确表征非平稳地震信号的瞬时频谱，保证算法可有效地处理局部拉伸校正、局部频谱滤波和薄层调谐作用的影响，克服传统傅立叶谱在表征上述问题中的缺陷；采用扫描分析方法建立起扫描Q值和深层频谱宽度恢复情况之间的对应关系；通过引入对数比值谱的平均导数，可有效克服噪声和薄层调谐的影响，使算法对于实际地震资料更具鲁棒性，适应性。

本发明采用的技术方案是：由叠前时间偏移直接扫描确定地层叠加品质因子方法，具体步骤包括：

1)用拖缆或测线记录由人工震源激发的经地下地层反射的地震信号，记录到磁带上；

2)从磁带上读取地震信号，根据地层性质和地震信号衰减特征初步确定Q值扫描步长和待扫描Q值；

3)对于每一个扫描Q值，应用常Q反Q滤波快速补偿叠前地震数据，时变地拓宽地震波频带；

4)对补偿后的叠前地震数据进行保幅叠前时间偏移，形成与该扫描Q值相对应的粘性补偿偏移后的数据体；

5)对于每一个粘性补偿数据体，沿着空间方向，抽取彼此相邻的3或5个共反射点道集，叠加后形成超道集；