用于衰减地震数据中的噪声的系统和方法
【技术领域】
[0001]本发明一般涉及用于处理地震数据的方法和系统，并且具体涉及用于衰减地震数据中的噪声的方法和系统。
【背景技术】
[0002]借助地震数据，可以有效率地进行油气储层的勘探和开采，所述地震数据必须恰当处理以便允许对地下特征(surface feature)的解释。通常，通过使用有源地震源(active seismic source)将地震能量注入到地下，所述地震能量接着被地下特征折射或反射并且在地震接收器处被记录来获得地震数据。在实践中，地震数据经常被可以在本质上是相干的或非相干(例如，随机)的噪声污染。
[0003]为改善最终地震图像和允许对地下特征的恰当解释，需要用于衰减地震数据中的噪声的有效率的且有效果的方法。
【发明内容】
[0004]本文描述了用于地震数据中噪声衰减的计算机实施的方法的各种方法的实施方式。
[0005]公开了一种用于衰减地震数据中的噪声的计算机实施的方法，所述地震数据代表感兴趣的地下区域。所述方法包括:接收代表地震信号和地震噪声的第一地震数据集和代表地震信号或噪声的第二地震数据集，将所述地震数据集变换到其中它们具有稀疏或可压缩表示形式的域，比较变换后的系数集合以在变换后的信号和噪声数据集中识别期望系数，选择变换后的信号和噪声数据集的期望系数以得到改善的系数集合，以及将所改善的系数集合逆变换以得到修改的地震数据集。依赖于所选择的系数，所修改的地震数据集可以代表信号或噪声。如果修改的地震数据集代表噪声，则可以从原始信号和噪声数据集中减去它以产生代表所述信号的数据集。
[0006]在另一种实施例中，公开了一种用以实施用于衰减地震数据中的噪声的方法的计算机系统，所述计算机系统包括数据源或存储设备、至少一个计算机处理器和用户界面。
[0007]在又一种实施例中，公开了一种包括计算机可读介质的制成品，所述计算机可读介质上包括计算机可读代码，所述计算机可读代码被配置成实施用于衰减地震数据中的噪声的方法，所述地震数据代表感兴趣的地下区域。
[0008]提供以上概要部分以便以简单形式介绍对概念的选择，所述概念在以下详细描述部分被进一步描述。所述概要不旨在识别所要求主题的关键特征或必要特征，也不旨在用以限制所要求主题的范围。而且，所要求主题不限于解决本公开的任何部分中指出的任何或全部不足的实施方式。
【附图说明】
[0009]参考以下描述、权利要求和附图，本发明的这些和其他特征会变得更好理解:
[0010]图1是说明根据本发明实施例的方法的流程图；
[0011 ]图2说明本发明实施例中的步骤；
[0012]图3示出了本发明衰减随机噪声的实施例的应用；
[0013]图4A示出了本发明衰减多个反射的一种实施例的应用；
[0014]图4B示出了本发明一种实施例的应用；
[0015]图5不出了本发明衰减随机噪声而不损害信号的另一种实施例的应用；
[0016]图6示意性说明了用于执行根据本发明实施例的方法的系统。
【具体实施方式】
[0017]本发明可以在要由计算机执行的计算机方法和系统的一般背景中描述和实施。这种计算机可执行指令可以包括可被用以执行特定任务和处理抽象数据类型的程序、例程、对象、组件、数据结构、和计算机软件技术。本发明的软件实施方式可以以用于各种计算平台和环境中的应用的不同语言编码。将要理解，本发明的范围和基本原理不限于任何具体的计算机软件技术。
[0018]而且，本领域技术人员会理解本发明可以使用硬件和软件配置的任何一种或组合来实践，包括但不限于，具有单个和/或多个处理器计算机的系统、手持设备、平板设备、可编程消费电子产品、迷你计算机、大型计算机等。本发明还可以在分布式计算环境中实践，其中任务由通过一个或多个数据通信网络链接的服务器或其他处理设备执行。在分布式计算环境中，程序模块可以位于包括存储器存储设备的本地和远程计算机存储介质这两者中。本发明还可被实践为井下传感器或测量设备的一部分或实践为实验室测量设备的一部分。
[0019]另外，诸如CD、预录盘或其他等价设备的与计算机处理器一同使用的制成品，可以包括其上记录有用于引导所述计算机处理器有助于本发明的实践和实施的有形计算机程序存储介质和程序装置。这种设备和制成品也落入到本发明的精神和范围内。
[0020]现在参考附图，将描述本发明的实施例。本发明可以多种方式被实施为，包括例如系统(包括计算机处理系统)、方法(包括计算机实施的方法)、装置、计算机可读介质、计算机程序产品、图形用户界面、门户网站、或有形地固定在计算机可读存储器中的数据结构。以下讨论本发明的若干实施例。所述附图仅说明了本发明的典型实施例并且因此不被认为限制本发明的范围和广度。
[0021]本发明涉及衰减地震数据中的噪声。本发明的一种实施例被示出为图1中的方法100。在操作12处，接收两个地震数据集。这些数据集中的一个数据集代表地震信号和噪声，例如在海底传感器的地震检波器处记录的垂直运动分量，并且可被称为第一地震数据集。另一个数据集仅代表地震信号，例如从海底传感器取得的水中检波器记录，并且可被称为第二地震数据集。这些例子不意在限制；有至少两个地震数据集可用的任何实例可被用以作为该方法的输入，其中至少一个数据集与另一数据集相比被认为是基本杜绝了噪声或基本杜绝了信号的。而且，输入数据集可以以各种方式被安排和/或预处理，作为例子，包括但不限于多个时移数据集、地震数据的堆栈或部分堆栈、基于诸如多个反射等的地震波的预期行为或基于地震处理算法的已知行为而生成的噪声或信号模型。所述数据集可以是炮道集(shot gather)、共检波点道集(common receiver gather)、共偏移距道集(commonoffset gather)、偏移向量片(offset vector tile)、共成像道集(common image gather)(角度或偏移)，并且可以以诸如纵向、横向、或深度切片/时间切片等的不同方向被安排；也可以使用这些的组合。本领域技术人员会理解，所述数据集的其他安排和预处理是可能的并且也可用作操作12的输入。此外，所述地震数据集可以是使用诸如气枪的主动源(active source)或被动源的记录。所述记录可以通过例如拖缆、海底电缆、海底节点、或诸如包括例如2-D线路勘测、3-D勘测、宽方位角和全方位角勘测的任意数目的接收器阵列配置中的地震检波器或加速计等的陆基传感器来得到。主动源可以同时地或顺序地以线性源几何或以诸如圆形采集技术(coil shooting)等的替代几何来点火。可以使用不同的源或接收器类型的组合。数据集可以是诸如基线和监控勘测等的时移数据。此外，一个或多个地震数据集可以是合成数据。本领域技术人员会理解，有许多生成适于所述第一和/或第二地震数据集的合成地震数据的方式。
[0022]在实施例中，可以有多于两个的输入数据集。一个输入数据集会代表信号和噪声并且与前述的第一地震数据集相同。另一个数据集可以代表刚才所述的信号或刚才所述的噪声的不同模型。在此实施例中，贯穿所述方法始终，会以如前述的对待第二地震数据集一样的方式来对待额外的信号或噪声模型。
[0023]在方法100的操作13处，第一和第二地震数据集被变换到其中它们具有稀疏或可压缩表示形式的域。所述变换可以使用多尺度多方向变换来完成。所述变换可以在数据的3-D体上或诸如纵向或横向部分或者时间或深度切片等的2-D部分上执行。所述数据集可被变换到曲波域或小波域。这些例子不意在限制；其中所变换的数据具有稀疏或可压缩表示形式的任何域都可在此方法中使用。此外，本领域技术人员会明白，也可能将1-D轨迹变换到其中所变换的数据具有稀疏或可压缩表示形式的域。
[0024]在操作14处，所变换的第一和第二数据集的代表系数互相比较。第一信号和噪声地震数据集中与第二信号数据集的代表系数接近的代表系数可被认为代表第二地震数据集中的信号。第一地震数据集中与第二地震数据集的代表系数接近的代表系数被认为是期望的。
[0025]图2示出了用于执行操作14的过程。在此，第一地震数据集(信号+噪声)已经被变换到曲波域并且在图形14A中绘出了它的系数。代表信号的两个其他输入数据集也已经被变换并且被表示成图形14B处的信号模型I和图形14C处的信号模型2。图形14B和14C示出了以粗虚线绘出的信号模型系数覆盖在以细实线示出的来自信号+噪声数据的对应系数上。图形14B和14C以水平虚线指示每个信号模型系数的限定阈值。所述限定阈值可以是某些默认阈值(例如±10%)，可以基于信号模型之间的系数大小的分布，或者可以是用户指定的。信号+噪声系数在图形14D中被指示为细线。在系数落入到基于信号模型的范围内时，这些系数被认为是期望的而其他系数被认定为是不期望的，在图形14D中用粗X指示。每个信号+噪声系数都与每个信号模型的对应系数比较。用户可以选择只接受充分接近所有信号模型的系数，或可以选择接受充分接近至少一个信号模型的系数。本领域技术人员会理解，使用一个或多个噪声模型而不是信号模型来执行此步骤也是可能的；信号模型或噪声模型都在本发明的范围内。
[0026]在操作15处，选择第一地震数据集的期望系数。这可以通过将不期望系数设置为零来完成，将不期望系数设置为零具有从第一地震数据集中移除与噪声有关的系数的效果。用于选择期望系数的其他方法是可能的，作为例子，包括但不限于，以这种方式修改不期望系数使得它们不同于期望系数或者修改期望系数。可以完成对期望系数的修改以从不期望系数中区分出它们或者强调所述期望系数的特定属性。在实施例中，第一地震数据集的期望代表系数是那些与所述信号有关的系数。将第一地震数据集的系数分成期望系数和不期望系数这两个系数集合也是可能的，并且将两个集合都传递给下一操作以便可观测到第一地震数据集的不期望部分。
[0027]在操作16处将第一地震数据集的期望系数逆变换以生成噪声衰减的第一地震数据集。如果不期望系数被分成单独的集合而不是被置为零，则操作16还可独立地变换不期望系数。
[0028]本领域技术人员还将理解在操作14处，将对不期望系数的指派改变成接近第二地震数据集的代表系数的那些系数也是可能的。这具有将第一地震数据集中的信号称为是不期望信号的效果，因此所述信号在操作15处通过置为零移除并且操作16的逆变换会生成噪声模型。
[