用于地震数据中噪声衰减的系统和方法
【技术领域】
[0001]本发明一般涉及用于处理地震数据的方法和系统，并且具体涉及用于衰减地震数据中的相干和非相干噪声的方法和系统。
【背景技术】
[0002]在地震数据的帮助下可以有效地实现油气藏(hydrocarbon reservoir)的勘探与开发，所述地震数据必须被恰当处理以允许解释地下特征(subsurface feature)。通过激活生成地震能量的地震源可以获得所述地震数据，所述地震能量传播穿过所述地下并且由地震接收器阵列记录。在实践中，地震数据经常被可以在获取时由运动(kinematics)引起的噪声污染，诸如记录的多个反射能量，或者被所使用的处理技术中的缺陷污染，诸如成像伪影。
[0003]需要用于衰减地震数据中的噪声的高效且有效的方法，以改善最终地震图像并且允许对地下特征的恰当解释。
【发明内容】
[0004]在此描述的是用于来自地震数据的地下特性描述的计算机实施的方法的各种方案的实施方式。
[0005]在一种实施例中，公开了一种用于处理地震数据的方法。所述方法包括:接收叠前地震数据集，所述叠前地震数据集用代表时间或深度的第一轴、代表诸如偏移或角度之类的平滑变化空间的至少一个冗余轴、以及至少一个其他轴来安排；转置所述叠前地震数据集以便所述第一轴变为所述冗余轴或其他轴；将所述转置的地震数据集切片(slice)以生成多个时间或深度切片，其中所述切片具有所述至少一个冗余轴和所述至少一个其他轴；处理所述多个切片以生成多个处理的切片；以及转置所述多个处理的切片以生成处理的叠前地震数据集，所述处理的叠前地震数据集用变为代表时间或深度的第一轴来安排。
[0006]在另一个实施例中，所述冗余轴可以代表两个或更多个时间延迟地震勘测之间的时间或方位角。
[0007]在还有的实施例中，在所述切片上执行的所述处理可以是2D、2.5D、3D或N_D滤波。所述滤波可以被设计成衰减噪声。
[0008]还有的实施例可以包括一种执行所述方法的系统。所述系统可以包括数据存储设备、配置成执行被设计成执行所述方法的步骤的计算机模块的处理器、以及用户接口。
[0009]还有的另一个实施例可以包括一种制品，所述制品包括在其上具有计算机可读代码的计算机可读介质，所述计算机可读代码被配置成实施所述方法。
[0010]提供上述概要部分以简化的形式介绍概念选集，所述概念在以下的详细描述部分被进一步描述。所述概要并不旨在确定所要求保护的主题的关键特征或本质特征，也不旨在用于限制所要求保护的主题的范围。而且，所要求保护的主题不限于解决本公开的任何部分中提及的任何或所有缺点的实施方案。
【附图说明】
[0011]参考下述描述、权利要求和附图，本发明的这些和其他特征会变得更好理解:
[0012]图1是2-D地震勘测的图；
[0013]图2是3-D地震勘测的图；
[0014]图3是说明根据本发明的一种方法的流程图；
[0015]图4A和4B是说明根据本发明的一种实施例的正被转变的示例地震数据集的图；
[0016]图5说明了在使用本发明的一种实施例之前和之后的地震数据部分；
[0017]图6说明了在使用本发明的一种实施例之前和之后的地震深度切片和所述切片之间的差别；以及
[0018]图7概要说明了一种用于执行根据本发明的一种实施例的方法的系统。
【具体实施方式】
[0019]本发明可以在系统和要由计算机执行的计算机方法的一般语境中描述和实施。此类计算机可执行指令可以包括可被用于执行特定任务和处理抽象数据类型的程序、例程、对象、组件、数据结构、和计算机软件技术。本发明的软件实施可以针对不同计算平台和环境中的应用以不同的语言编码。将要理解，本发明的范围和基础原理不限于任何特定的计算机软件技术。
[0020]而且，本领域技术人员会理解本发明可以使用硬件和软件配置的任意一个或组合来实践，包括但不限于具有单个和/或多个处理器计算机的系统、手持设备、可编程消费电子产品、迷你计算机、大型主机等。本发明也可以在分布式计算环境中实践，其中任务由通过一个或多个数据通信网络链接的服务器或其他处理设备执行。在分布式计算环境中，程序模块可以位于包括存储器存储设备的本地和远程计算机存储介质中。本发明还可以实践为井下传感器(down-hole sensor)或测量设备的一部分或实践为实验室测量设备的一部分。
[0021 ] 此外，与计算机处理器一同使用的诸如⑶、预刻录盘或其他等价设备之类的制品，可以包括其上记录有用于引导计算机处理器便于本发明的实施和实践的计算机程序存储介质和程序装置。此类设备和制品也落入到本发明的精神和范围内。
[0022]现在参考所述附图，将描述本发明的实施例。本发明可以以多种方式实施，包括例如作为系统(包括计算机处理系统)、方法(包括计算机实施的方法)、装置、计算机可读介质、计算机程序产品、图形用户界面、门户网站、或可触知地固定在计算机可读存储器中的数据结构。以下讨论本发明的若干实施例。所述附图仅说明了本发明的典型实施例并且因此不被认为是对其范围和广度的限制。
[0023]本发明涉及衰减由地震数据中的成像伪影和/或多个反射能量引起的噪声。在此描述的地震数据是通过记录从活动地震源(例如，振动检波器(Vitooseis)、气枪、炸药等)穿过感兴趣的地下区域传播到地震接收器的地震能量来获得。地震勘测可以包括同时和/或顺序启动的多个地震源，和多个地震接收器。在被记录之后，所述地震数据可以经受任意数目的处理步骤，例如包括正常时差校正(normal moveout correct1n)或叠前偏移(prestack migrat1n)。这些处理步骤不旨在作为限制；本领域技术人员会认识到在应用本发明的实施例之前存在可以使用的多个地震数据处理选项。
[0024]图1示出了使用单个发射(shot)12和四个接收器14的简单2D地震勘测。发射12生成由在此的从表面10反射并且返回到接收器14的射线16代表的地震波。图2示出了也使用发射12、接收器14和反射表面10的简单3D地震勘测。在此情况中，仅示出了射线16中的四条以简化所述图，但是本领域技术人员会认识到所有接收器会记录来自所述发射的能量。在上述两种勘测中，记录的数据可以由指示所述接收器到所述发射的关系的地理位置(例如，共中点(common midpoint, CMP)、主测线位置(inline posit1n)、联络线位置(crossline posit1n)等)和相对位置(例如，偏移、方位角等)来索引。所述地理位置也可以称为地理轴。还可以存在多个地理轴(例如，CMP_X, CMP_Y) ο所述相对位置可以称为相对轴并且还可以具有多个维度(例如，offset_x, offset_y)。这些例子不意味着限制；地理位置和相对位置都可以是N维的并且因此具有N个轴。
[0025]当所述地震数据由校正接收器之间的传播时间差的算法处理时，来自反射表面的地震能量会在沿着所述相对轴的每个轨迹的相同时间或深度处发生。实际上，地震事件沿着相对轴变平，以便沿着特定事件的幅度变化是平滑的。所述事件的平滑变化本质使得沿着那个轴的数据冗余。针对本文的目的，所述相对轴称为冗余轴。
[0026]除了上述的相对轴，所述冗余轴还可能可以代表来自两个或更多个时间延迟数据集的数据。在此实例中，冗余轴会指示勘测之间的时间差别。
[0027]图3示出了本发明的一种实施例的流程图。图4A和4B的地震数据的图也说明了方法30。在操作32接收到的叠前地震数据集可以看起来类似于地震数据体40。所述叠前地震数据集具有至少一个冗余轴。在图4A中，所述地震数据体的轴被定向为由位于所述数据体左侧的方向箭头指示。这些方向箭头是示出地理轴的40X、示出冗余轴的40H、以及指示深度或时间轴的40Z。对于地震数据体40，所述时间或深度轴40Z也通常称为竖直或第一轴。地震数据体40具有三个事件(反射表面)40A、40B、和40C。沿着地理轴40X，事件40A是倾斜反射层(dipping reflector)并且事件40B和40C是非对称背斜层(asymmetricalanticline)。沿着冗余轴40H，所有三个事件40A、40B、和40C都是平的。地震数据体40也含有噪声41，示出为40Z-40H平面中的俯冲线(swooping line)和40X-40H平面中的点。此噪声是简单的例子并且不意味着限制。本领域技术人员会理解存在着可以污染地震数据的多种类型的相干和非相干噪声。
[0028]在方法30的操作33处，叠前地震数据集被转置。所述转置将地震数据体的布置转换以便时间或深度轴不再是第一轴。图4A示出了转置的地震数据集40T，其具有用于40X、40H和40Z的新的方向箭头集合。在此例子中，第一轴现在是地理轴40X。它也可以是冗余轴40H。事件40A、40B和40C被标识为噪声41。
[0029]在叠前地震数据集已经被转置以生成转置的地震数据集之后，所述转置的地震数据集在操作34处被切片为时间或深度切片。每个此类切片具有至少一个冗余轴。在图4A中，转置的地震数据集40T已经被切片成切片42A-F。每个此类切片具有地理轴40X(竖直轴)和冗余轴40H(水平轴)。所述切片不必具有地理轴；所述地震数据体可能可以具有多个冗余轴以便所述切片可以具有冗余轴。尽管术语“切片”暗