[发明专利]一种利用数据自适应性时窗估算几何属性的方法

说明书

本发明公开了一种利用数据自适应性时窗估算几何属性的方法，通过平均振幅能量数据F(t，f)，计算自适应时窗h(t，f)；沿构造走向提取时窗内地震数据；计算时窗内各地震道的平均地震道u_avg(t)；计算时窗内各地震道与平均地震道的相关系数r(t，x，y)，并使用平均地震道与各地震道相关系数的乘积，得到各地震道的分析地震道u′(t，x，y)；计算各分析地震道的能量之和，得到相干能量值；计算时窗内各原始地震道的能量之和，得到总体能量之；相干能量值与总体能量值的比值即为相干值，最后计算整个三维地震数据体的相干体并输出。本发明的有益效果是使用拟合地震道来代替原始地震道可以减少相干体计算过程中的随机数据。

技术领域

本发明属于地质学技术领域，涉及一种利用数据自适应性时窗估算几何属性的方法。

背景技术

地震几何属性如相干属性通常用来描述地震解释和油藏描述中的断层、河道和其他地质特征。传统上，我们使用单一自定义固定分析时窗计算整个数据体的属性。一般，较小的时窗产生更清晰的地质边界，但对噪音却非常敏感。相反，较大的时窗可以减少随机噪音的影响，但在横向上很难区分断层与河道的边界，垂向上容易将断层边界与地层界面相混淆。这主要是因为地震数据具有较低的信噪比，在地震数据的主周期范围内，在相干体中地层边界一般随着时窗的增大而增强。当时窗大于主周期时，时窗稍微提高就会导致浅部和深部的地层成像结果发生混淆。由于地震信号的衰减损失和地层速度的增加，三维地震数据的垂向和横向分辨率随着深度的增加而发生变化。因此，一个在浅部地层可以得到最佳成像结果的时窗在深部地层中往往就会太小。常用的解决这个问题的方案是在一个合适的固定时窗内计算地震属性，然后拼接计算结果，但常常要冒着数据体垂向连续性减小的危险。这些算法在实现过程中大多数使用固定数量的道和采样数去分析整个地震数据体。然而，由于表层频率的损失，随着深部的增加地震速度逐渐增加，而地震波入射角逐渐减小，地震数据的时间分辨率和横向分辨率逐渐减小。因此，浅部地层最优的固定分析时窗在深部地层可能不是最优的。

发明内容

本发明的目的在于提供一种利用数据自适应性时窗估算几何属性的方法，解决了地震几何属性固定分析时窗在浅部地层可以得到最佳成像结果的时窗应用到深部地层中往往就会太小的问题。

本发明所采用的技术方案是按照以下步骤进行：

步骤1：通过平均振幅能量数据F(t，f)，计算自适应时窗h(t，f)；

其中，t为时间，f为频率；

步骤2：沿构造走向提取时窗内地震数据；

通过自适应时窗、构造倾角与倾向、还有地震振幅，沿构造走向按照计算的自适应时窗大小提取时窗内地震数据，即地震数据道u(t，x，y)，其中x,y分别为三维地震工区x和y方向时窗大小；

步骤3：计算时窗内各地震道的平均地震道u_avg(t)；

其中，J为参与属性计算的地震道数；

步骤4：计算时窗内各地震道与平均地震道的相关系数r(t，x，y)，并使用平均地震道与各地震道相关系数的乘积，得到各地震道的分析地震道u′(t，x，y)；

u′(t，x，y)＝u_avg(t)·r(t，x，y)；

步骤5：计算各分析地震道的能量之和，得到相干能量值；计算时窗内各原始地震道的能量之和，得到总体能量值；相干能量值与总体能量值的比值即为相干值；

步骤6：循环步骤1到步骤5，计算整个三维地震数据体的相干体并输出。