专利名称：用于地震道分析的系统和方法
技术领域：
本发明一般涉及地质数据的处理，尤其涉及用于减小采样数据中的混叠影响的系 统。
背景技术：
地震勘测往往产生稀疏和/或非均勻采样的数据。这可能源自非均勻的震源和接 收器位置、系统和导航不完善、地物障碍以及大型仪器阵列的过高费用。另一方面，许多数 据处理和可视化系统和方法需要密集和规则采样的数据作为输入。因此，常常需要变换数 据，使得数据具有适当的采样基础。
发明内容
本发明的实施例的一些方面提供了一种处理代表与地下区域有关的信息的采样 地震道的方法，包括使用离散傅立叶变换来估计地震道的傅立叶系数；通过将相干准则 应用于估计的傅立叶系数来选择估计的傅立叶系数，以减小所述傅立叶系数所代表的混叠 能量的量；在k_空间中重新正交化所述傅立叶系数；以及使用重新正交化的傅立叶系数来 产生地下区域的地震图像。本发明的实施例的一些方面包括一种成像地下特征的方法，包括接收包括与地 下特征有关的信息的采样地震道；使用离散傅立叶变换来估计地震道的傅立叶系数；通过 将相干准则应用于估计的傅立叶系数来选择估计的傅立叶系数，以减小所述傅立叶系数所 代表的混叠能量的量；在k_空间中重新正交化所述傅立叶系数；以及应用逆傅立叶变换来 产生规则化地震数据。本发明的实施例的一些方面可以包括编有用于执行上述方法或控制上述系统的 计算机可执行指令的计算机可读介质。本发明的实施例的一些方面可以包括一种合并有上述系统并被配置和安排成依 照上述方法来提供对所述系统的控制的系统。这样的系统可以合并有例如被编程成允许用 户依照所述方法或其它方法来控制设备的计算机。通过参照形成本说明书的一部分、相同标号在各个图形中表示相应部件的附图对 如下描述和所附权利要求书加以研究，本发明的这些和其它目的、特征、和特性，以及操作 方法、结构的相关元件的功能、部件的组合、和制造成本将变得更加显而易见。但是，不言而 喻，这些附图只用于例示和描述的目的，而无意作为限制本发明的定义。正如用在说明书和 权利要求书中的那样，除非上下文另有明确指明，单数形式“一个”、“一种”、和“该”也包括 复数指示物。
图1例示了采样道的集合的例子；图2例示了图1的集合的f_k谱；
4[oo10] 图3例示了在应用依照本发明一个实施例的方法之后图2的[—k谱；[oo11] 图4例示了图3的[—k谱的逆傅立叶变换；以及[oo12] 图5是执行依照本发明实施例的方法的系统的实施例的示意性例示图。
具体实施方式
[oo13] 地震数据可被编译成二维断面或集合。在图l中例示了一个例子。一般说来，可以将一组这样的数据表示成P(t，X)，其中，t和X分别代表时间和空间。在图l的例子中，X轴代表以米为单位的距离(空间)，而y轴代表以秒为单位的时间。对于这个例子，样本的间隔是不规则的，并且是稀疏采样的。为了清楚起见，这些道大部分是无噪声的，只示出了许多记录的事件lo，而在每个大振幅事件之外，道基本上是平的，具有零振幅。而且，已经将数据限制在窗口内，使得地震事件相当线性。[oo14] 对于被良好采样(即，足够的采样频率和在时间上足够规则)的数据，可以依照方程l，使用快速傅立叶变换([[T)来计算傅立叶变换[oo]5] P(‘o，t)一[T+[P(t，X)](1)[oo16] 其中，[T+是时间傅立叶变换算子。[oo17] 对于每个频率0，把数据看作X的函数，可以写出方程2[oo]8] h(X)一P(0，X)(2)[oo19] 如上所述，可以作出数据被均匀采样的假设，使得可以将这个函数的均匀采样序列表达成方程3
h。一h(n厶X)，n—o，l，．．．，N—l，(3)
可以将正向离散傅立叶变换(D[T)表达成
权利要求
1.一种处理代表与地下区域有关的信息的采样地震道的方法，包含 使用离散傅立叶变换来估计地震道的傅立叶系数；通过将相干准则应用于估计的傅立叶系数来选择估计的傅立叶系数，以减小所述傅立 叶系数所代表的混叠能量的量；在k_空间中重新正交化所述傅立叶系数；以及 使用重新正交化的傅立叶系数来产生地下区域的地震图像。
2.按照权利要求1所述的方法，进一步包含对于每个频率，选择混叠抑制函数值最大和离散傅立叶系数振幅最大的傅立叶系数。
3.按照权利要求1所述的方法，其中，所述相干准则包含抑制函数
4.按照权利要求1所述的方法，其中，所述相干准则包含计算沿着穿过k空间中的原点 的线段的相干性的函数。
5.按照权利要求1所述的方法，其中，所述相干准则包含振幅谱的叠加
6.按照权利要求1所述的方法，进一步包含在重新正交化傅立叶系数之后，应用逆傅 立叶变换以获得规则化地震数据。
7.一种成像地下特征的方法，包含接收包括与地下特征有关的信息的采样地震道； 使用离散傅立叶变换来估计地震道的傅立叶系数；通过将相干准则应用于估计的傅立叶系数来选择估计的傅立叶系数，以减小所述傅立 叶系数所代表的混叠能量的量；在k-空间中重新正交化所述傅立叶系数；以及 应用逆傅立叶变换来产生规则化地震数据。
8.按照权利要求7所述的方法，进一步包含根据所述规则化地震数据来显示地震道图像。
9.按照权利要求7所述的方法，进一步包含对于每个频率，选择混叠抑制函数值最大和快速傅立叶系数振幅最大的傅立叶系数。
10.按照权利要求7所述的方法，其中，所述相干准则包含抑制函数
11.按照权利要求7所述的方法，其中，所述相干准则包含计算沿着穿过k空间中的原 点的线段的相干性的函数。
12.按照权利要求7所述的方法，进一步包含在确定在所述地下特征附近存在烃资源 的概率的过程中使用所述规则化地震数据。
13. 一种成像与地下特征有关的地震数据的系统，包含具有计算机可读数据的数据存储设备，所述计算机可读数据包括与采样地震道有关的 数据，所述采样地震道包括与地下特征有关的信息；以及处理器，被配置和安排成执行存储在处理器可访问的存储器中的机器可执行指令，以 便执行包含如下步骤的方法使用离散傅立叶变换来估计地震道的傅立叶系数；通过将相干准则应用于估计的傅立叶系数来选择估计的傅立叶系数，以减小所述傅立 叶系数所代表的混叠能量的量；在k-空间中重新正交化所述傅立叶系数；以及 应用逆傅立叶变换来产生规则化地震数据。
全文摘要
处理地震数据以减小或消除由例如稀疏或不规则采样引起的混叠。迭代方法包括与估计傅立叶系数的振幅的函数结合使用的抑制函数，这两个函数一起作用以减小混叠能量的影响并优先选择真正能量。计算步骤主要在k空间中执行，而不返回到x空间，从而减小计算成本。
文档编号G01V1/32GK102105814SQ200980128648
公开日2011年6月22日 申请日期2009年6月15日 优先权日2008年6月19日
发明者L·L·张 申请人:雪佛龙美国公司