地震数据去噪方法及系统技术方案

公开了一种地震数据去噪方法及系统。该方法可包括以下步骤：对时域的地震数据进行变换以获得频域地震数据；对频域地震数据求取预测算子，并基于预测算子对频域地震数据进行去噪处理；以及对经去噪的频域地震数据进行变换以获得经去噪的时域地震数据，其中，通过如下公式进行去噪处理，i,j不能同时为0：其中，表示某道地震数据在频率点f的频域地震数据Do,p(f)的经去噪的频域地震数据，1≤o≤M，1≤p≤N，M、N分别为地震数据在x、y方向的总道数，Ai,j(f)表示预测算子的各个分量，lx＝(lfx-1)/2，ly＝(lfy-1)/2，lfx和lfy分别为预测算子Ai,j(f)在x、y方向的长度，Do+i,p+j(f)表示上述某道周围相邻道中，与频域地震数据Do,p(f)同频率的频域地震数据。

【技术实现步骤摘要】

本公开涉及地球物理勘探领域，更具体地，涉及一种地震数据去噪方法及系统。
技术介绍
在地球物理勘探领域，地震油气勘探的难度和复杂程度不断增加，对地震勘探精度的要求也越来越高。在对地震资料的处理中，如何提高地震资料的分辨率是一项重要工作，而提高分辨率的关键是实现高精度的去噪。目前，在去噪方面常用的一种技术是基于高覆盖次数的多次叠加。该技术要求高的覆盖次数，覆盖次数越高，对高频段信噪比的改善越明显。然而提高覆盖次数要求在设计观测系统时，增加接收点数、加大排列长度，提高空间采样密度，加大了野外采集的成本。在去噪方面常用的另外一种技术是预测算子外推技术。该技术对于信噪比较高的中频段数据，可以直接求取中频预测算子，实现中频段去噪；而对于信噪比较低的高频段数据，可以利用信噪较高的中频段数据外推出高频预测算子，再实现高频段去噪。专利技术人发现，现有的地震数据去噪方法在对三维地震数据进行去噪处理时，经去噪的地震数据的信噪比不高，尤其在保护弯曲同相轴，断点和弱信号等方面存在不足。因此，有必要开发一种针对三维地震数据的高精度的地震数据去噪方法。公开于本公开
技术介绍
部分的信息仅仅旨在加深对本公开的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本公开提供了一种地震数据去噪方法及系统，该方法可以充分利用地震数据的三维信息实现去噪，从而有效改善地震数据的信噪比，提高处理精度。根据本公开的一方面，提出了一种地震数据去噪方法，该方法可以包括以下步骤：对时域的地震数据进行变换以获得频域地震数据；对频域地震数据求取预测算子，并基于预测算子对频域地震数据进行去噪处理；以及对经去噪的频域地震数据进行变换以获得经去噪的时域地震数据，其中，通过如下公式进行去噪处理，i,j不能同时为0：D^o,p(f)=Σi=-lxlxΣj=-lylyAi,j(f)·Do+i,p+j(f)]]>其中，表示某道地震数据的在频率点f的频域地震数据Do,p(f)的经去噪的频域地震数据，1≤o≤M，1≤p≤N，M、N分别为地震数据在x、y方向的总道数，Ai,j(f)表示预测算子的各个分量，lx＝(lfx-1)/2，ly＝(lfy-1)/2，lfx和lfy分别为预测算子Ai,j(f)在x、y方向的长度，Do+i,p+j(f)表示上述某道周围相邻道中，与频域地震数据Do,p(f)同频率的频域地震数据。根据本公开的另一方面，提出了一种地震数据去噪系统，该系统可以包括以下单元：用于对时域的地震数据进行变换以获得频域地震数据的单元；用于对频域地震数据求取预测算子，并基于预测算子对频域地震数据进行去噪处理的单元；以及用于对经去噪的频域地震数据进行变换以获得经去噪的时域地震数据的单元，其中，通过如下公式进行去噪处理，，i,j不能同时为0：D^o,p(f)=Σi=-lxlxΣj=-lylyAi,j(f)·Do+i,p+j(f)]]>其中，表示某道地震数据的在频率点f的频域地震数据Do,p(f)的经去噪的频域地震数据，1≤o≤M，1≤p≤N，M、N分别为地震数据在x、y方向的总道数，Ai,j(f)表示预测算子的各个分量，lx＝(lfx-1)/2，ly＝(lfy-1)/2，lfx和lfy分别为预测算子Ai,j(f)在x、y方向的长度，Do+i,p+j(f)表示上述某道周围相邻道中，与频域地震数据Do,p(f)同频率的频域地震数据。本公开的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本公开的特定原理。附图说明通过结合附图对本公开的示例性实施例进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开的示例性实施例中，相同的参考标号通常代表相同部件。图1是示出了根据本公开的示例性实施例的地震数据去噪方法的步骤的框图。图2a-2f是示出了本公开的地震数据去噪方法与现有技术的比较的示意图，其中，图2a为有效信号，图2b为信噪比谱，图2c为含噪数据，图2d为正确的算子，图2e为f-x域预测滤波求得的算子，以及图2f为本公开的方法求得的算子。图3a-3d是示出了本公开的地震数据去噪方法与现有技术的50-80Hz频段比较的示意图，其中，图3a为有效信号，图3b为含噪数据，图3c为f-x域预测滤波结果，图3d为本公开的方法求得的算子。图4a-4d，图5a-5d和图6a-6d是示出了本公开的地震数据去噪方法与现有技术的比较的示意图，其中，图4a为理论合成无噪记录，图4b为含噪记录，图4c为二维滤波结果，图4d为本公开的方法的三维非因果滤波结果，图5a为理论合成无噪记录，图5b为含噪记录，图5c为二维滤波结果，图5d为本公开的方法的滤波结果，图6a为理论合成无噪记录，图6b为含噪记录，图6c为二维滤波结果，以及图6d为本公开的方法的滤波结果。图7a-图7e、图8a-8c、图9a-9b以及图10a-10c是示出了本公开的地震数据去噪方法与现有技术的比较的示意图，其中，图7a为某段单炮记录(经过炮域动校正和压制面波的处理)，图7b为RNA处理结果，图7c为本公开的处理结果，图7d为图7b虚线框内数据的放大显示，图7e为图7c虚线框内数据的放大显示，图8a为图7a的频谱，图8b为图7b的频谱，图8c为图7c的频谱，图9a为图7a减去图7b，图9b为图7a减去图7c，图10a为图7a经50-80Hz带通滤波，图10b为图7b经50-80Hz带通滤波，以及图10c为图7c经50-80Hz带通滤波。具体实施方式下面将参照附图更详细地描述本公开的优选实施例。虽然附图中显示了本公开的优选实施例，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。实施例1图1示出了根据本公开的示例性实施例的地震数据去噪方法的流程图。在该实施例中，地震数据去噪方法可以包括以下步骤：步骤101，对时域的地震数据进行变换以获得频域地震数据；步骤102，对频域地震数据求取预测算子，并基于预测算子对频域地震数据进行去噪处理；以及步骤103，对经去噪的频域地震数据进行变换以获得经去噪的时域地震数据，其中，通过如下的公式(1)进行去噪处理，i,j不能同时为0：D^o,p(f)=Σi=-lxlxΣj=-lylyAi,j(f)·Do+i,p+j(f)---(1)]]>其中，表示某道地震数据的在频率点f的频域地震数据Do,p(f)的经去噪的频域地震数据，1≤o≤M，1≤p≤N，M、N分别为地震数据在x、y方向的总道数，Ai,j(f)表示预测算子的各个分量，lx＝(lfx-1)/2，ly＝(lfy-1)/2，lfx和lfy分别为预测算子Ai,j(f)在x、y方向的长度，Do+i,p+j(f)表示与频域地震数据Do,p(f)周围相邻本文档来自技高网...

【技术保护点】
一种地震数据去噪方法，所述方法包括以下步骤：对时域的地震数据进行变换以获得频域地震数据；对所述频域地震数据求取预测算子，并基于所述预测算子对频域地震数据进行去噪处理；以及对经去噪的频域地震数据进行变换以获得经去噪的时域地震数据，其中，通过如下公式进行去噪处理，i,j不能同时为0：D^o,p(f)=Σi=-lxlxΣj=-lylyAi,j(f)·Do+i,p+j(f)]]>其中，表示某道地震数据在频率点f的频域地震数据Do,p(f)的经去噪的频域地震数据，1≤o≤M，1≤p≤N，M、N分别为地震数据在x、y方向的总道数，Ai,j(f)表示预测算子的各个分量，lx＝(lfx‑1)/2，ly＝(lfy‑1)/2，lfx和lfy分别为预测算子Ai,j(f)在x、y方向的长度，Do+i,p+j(f)表示上述某道周围相邻道中，与频域地震数据Do,p(f)同频率的频域地震数据。

【技术特征摘要】
1.一种地震数据去噪方法，所述方法包括以下步骤：对时域的地震数据进行变换以获得频域地震数据；对所述频域地震数据求取预测算子，并基于所述预测算子对频域地震数据进行去噪处理；以及对经去噪的频域地震数据进行变换以获得经去噪的时域地震数据，其中，通过如下公式进行去噪处理，i,j不能同时为0：D^o,p(f)=Σi=-lxlxΣj=-lylyAi,j(f)·Do+i,p+j(f)]]>其中，表示某道地震数据在频率点f的频域地震数据Do,p(f)的经去噪的频域地震数据，1≤o≤M，1≤p≤N，M、N分别为地震数据在x、y方向的总道数，Ai,j(f)表示预测算子的各个分量，lx＝(lfx-1)/2，ly＝(lfy-1)/2，lfx和lfy分别为预测算子Ai,j(f)在x、y方向的长度，Do+i,p+j(f)表示上述某道周围相邻道中，与频域地震数据Do,p(f)同频率的频域地震数据。2.根据权利要求1所述的地震数据去噪方法，进一步包括，确定地震数据的中频段和高频段，以针对所述中频段和所述高频段分别求取预测算子。3.根据权利要求2所述的地震数据去噪方法，其中，对所述中频段和所述高频段分别求取预测算子包括：对所述中频段的频域地震数据直接求取预测算子，以及对所述高频段的频域地震数据通过预测算子外推求取预测算子。4.根据权利要求1所述的地震数据去噪方法，进一步包括，设定滑动窗体，通过所述滑动窗体在地震数据中滑动而分段进行地震数据去噪。5.一种基于算子外推的三维地震数据...

【专利技术属性】
技术研发人员：魏巍，黄忠玉，季玉新，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京;11