CN115685320A - 一种地震炮集数据去噪方法和装置 - Google Patents
一种地震炮集数据去噪方法和装置 Download PDFInfo
- Publication number
- CN115685320A CN115685320A CN202211366969.1A CN202211366969A CN115685320A CN 115685320 A CN115685320 A CN 115685320A CN 202211366969 A CN202211366969 A CN 202211366969A CN 115685320 A CN115685320 A CN 115685320A
- Authority
- CN
- China
- Prior art keywords
- gather data
- seismic
- shot gather
- seismic shot
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000012937 correction Methods 0.000 claims abstract description 49
- 238000012545 processing Methods 0.000 claims abstract description 26
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 24
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 238000001914 filtration Methods 0.000 claims description 19
- 238000001514 detection method Methods 0.000 claims description 16
- 239000000523 sample Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 238000009795 derivation Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 6
- 238000004891 communication Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000036039 immunity Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005293 physical law Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种地震炮集数据去噪方法和装置，涉及地震勘探的技术领域，包括：获取待处理地理区域的地震炮集数据、动校正速度和叠后地震数据；基于动校正速度计算地震炮集数据对应的零偏移距信息；利用平面波分解算法对叠后地震数据进行处理，得到地表射线参数；基于零偏移距信息、动校正速度和地表射线参数构建地震炮集数据的局部倾角；基于地震炮集数据的局部倾角对地震炮集数据进行去噪处理，得到待处理地理区域去噪后的地震炮集数据。本发明方法是通过零偏移距信息、动校正速度和地表射线参数间接构建地震炮集数据的局部倾角场，所以受噪音影响较小，具有较强的抗噪性和较高的工业应用价值，可实现高质量的地震数据去噪效果。
Description
技术领域
本发明涉及地震勘探的技术领域，尤其是涉及一种地震炮集数据去噪方法和装置。
背景技术
地震数据经常被称为噪声的干扰能量破坏。这种不必要的噪声会干扰地震数据的连续性，并掩盖有效的信息。此外，许多后续处理任务如地震属性分析、振幅随偏移量变化(AVO)分析和自动解释都受到噪声的负面影响。因此，增强高信噪比地震事件对提高地震勘探效率至关重要。
现有技术中，局部倾角信息常用于地震数据去噪，然而常规方法是数据驱动的，受噪音影响较大，在信噪比下难以获得准确的局部倾角信息，进而导致地震数据去噪效果欠佳。
发明内容
本发明的目的在于提供一种地震炮集数据去噪方法和装置，以缓解了现有技术中地震炮集数据去噪方法存在的去噪效果不佳的技术问题。
第一方面，本发明提供一种地震炮集数据去噪方法，包括：获取待处理地理区域的地震炮集数据、动校正速度和叠后地震数据；基于所述动校正速度计算所述地震炮集数据对应的零偏移距信息；利用平面波分解算法对所述叠后地震数据进行处理，得到地表射线参数；基于所述零偏移距信息、所述动校正速度和所述地表射线参数构建所述地震炮集数据的局部倾角；基于所述地震炮集数据的局部倾角对所述地震炮集数据进行去噪处理，得到所述待处理地理区域去噪后的地震炮集数据。
在可选的实施方式中，计算所述地震炮集数据对应的零偏移距信息，包括：获取所述地震炮集数据中的炮点位置信息、检波点位置信息和非零偏移距走时；基于所述炮点位置信息和检波点位置信息计算所述地震炮集数据对应的共中心点位置信息和偏移距信息；基于所述非零偏移距走时、所述偏移距信息和所述动校正速度计算所述地震炮集数据对应的零偏移距走时；基于所述共中心点位置信息和所述零偏移距走时确定所述零偏移距信息。
在可选的实施方式中，利用平面波分解算法对所述叠后地震数据进行处理，得到地表射线参数，包括：对所述叠后地震数据进行平面波分解，得到所述叠后地震数据在空间方向的第一分量和在时间方向的第二分量；基于所述第一分量和所述第二分量计算所述地表射线参数。
在可选的实施方式中，基于所述零偏移距信息、所述动校正速度和所述地表射线参数构建所述地震炮集数据的局部倾角，包括：基于所述零偏移距信息、所述动校正速度和所述地表射线参数，构建检波点移动目标距离的情况下，地震炮集数据的共反射面元走时方程；基于所述共反射面元走时方程计算所述地震炮集数据的局部倾角。
在可选的实施方式中，基于所述共反射面元走时方程计算所述地震炮集数据的局部倾角，包括：对所述共反射面元走时方程进行求导数处理，得到求导方程；将求导方程中的所述目标距离设置0的结果作为所述地震炮集数据的局部倾角。
在可选的实施方式中，基于所述地震炮集数据的局部倾角对所述地震炮集数据进行去噪处理，包括：沿着所述地震炮集数据的局部倾角方向对所述地震炮集数据进行均值滤波处理，得到所述待处理地理区域去噪后的地震炮集数据。
在可选的实施方式中，沿着所述地震炮集数据的局部倾角方向对所述地震炮集数据进行均值滤波处理，包括：基于算式
对所述地震炮集数据进行均值滤波，得到所述待处理地理区域去噪后的地震炮集数据；其中，R(xs,xr,t)表示所述去噪后的地震炮集数据，U(xs,xr,t)表示所述地震炮集数据，xs表示所述地震炮集数据中的炮点位置信息，xr表示所述地震炮集数据中的检波点位置信息，t表示所述地震炮集数据中与(xs,xr)对应的非零偏移距走时，n表示均值滤波窗口的地震道数量，Δxi表示当前检波点位置xr相对第i道地震数据的位移，Δti表示当前检波点位置xr相对第i道地震数据的时移，σ表示所述地震炮集数据的局部倾角。
第二方面，本发明提供一种地震炮集数据去噪装置，包括：获取模块，用于获取待处理地理区域的地震炮集数据、动校正速度和叠后地震数据；计算模块，用于基于所述动校正速度计算所述地震炮集数据对应的零偏移距信息；分解模块，用于利用平面波分解算法对所述叠后地震数据进行处理，得到地表射线参数；构建模块，用于基于所述零偏移距信息、所述动校正速度和所述地表射线参数构建所述地震炮集数据的局部倾角；去噪模块，用于基于所述地震炮集数据的局部倾角对所述地震炮集数据进行去噪处理，得到所述待处理地理区域去噪后的地震炮集数据。
第三方面，本发明提供一种电子设备，包括存储器、处理器，所述存储器上存储有可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现前述实施方式中任一项所述的地震炮集数据去噪方法的步骤。
第四方面，本发明提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令被处理器执行时实现前述实施方式中任一项所述的地震炮集数据去噪方法。
本发明提供的地震炮集数据去噪方法，包括：获取待处理地理区域的地震炮集数据、动校正速度和叠后地震数据；基于动校正速度计算地震炮集数据对应的零偏移距信息；利用平面波分解算法对叠后地震数据进行处理，得到地表射线参数；基于零偏移距信息、动校正速度和地表射线参数构建地震炮集数据的局部倾角；基于地震炮集数据的局部倾角对地震炮集数据进行去噪处理，得到待处理地理区域去噪后的地震炮集数据。
本发明方法不直接依赖地震炮集数据进行倾角计算，而是通过零偏移距信息、动校正速度和地表射线参数间接构建地震炮集数据的局部倾角场，由于地表射线参数信息是基于信噪比较高的叠后数据中确定的，所以受噪音影响较小。因此，本发明提供的地震炮集数据去噪方法具有较强的抗噪性和较高的工业应用价值，可实现高质量的地震数据去噪效果。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种地震炮集数据去噪方法的流程图；
图2为本发明实施例提供的一种计算地震炮集数据对应的零偏移距信息的流程图；
图3为本发明实施例提供的一种地震炮集数据去噪装置的功能模块图；
图4为本发明实施例提供的一种电子设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合附图，对本发明的一些实施方式作详细说明。在不冲突的情况下，下述的实施例及实施例中的特征可以相互组合。
实施例一
图1为本发明实施例提供的一种地震炮集数据去噪方法的流程图，如图1所示，该方法具体包括如下步骤：
步骤S102，获取待处理地理区域的地震炮集数据、动校正速度和叠后地震数据。
具体的，首先获取待处理地理区域的地震炮集数据U(xs,xr,t)、动校正速度V(m,t0)和叠后地震数据M(m,t0)，其中，xs表示炮点位置信息(也即，坐标信息)，xr表示检波点位置信息，t表示非零偏移距走时，m表示共中心点位置信息信息，共中心点是指炮点和检波点的中点，t0表示零偏移距走时。
V(m,t0)和M(m,t0)中的共中心点位置信息m和零偏移距走时t0是上述两种数据体对应的两个维度，也即，空间方向和时间方向，类似于平面方程Z(x,y)中的x和y。而地震炮集数据中的t是非零偏移距走时，其虽然也是时间方向，但与t0的物理意义不同，且t和t0之间存在一定的联系(符合动校正公式)。
步骤S104，基于动校正速度计算地震炮集数据对应的零偏移距信息。
基于上文中的描述可知，零偏移距走时对于V(m,t0)和M(m,t0)是已知的，但是对于地震炮集数据U(xs,xr,t)却是待求的。在本发明实施例中，零偏移距信息包括：共中心点位置信息m和零偏移距走时t0，因此在已知地震炮集数据(包括：炮点位置信息、检波点位置信息和非零偏移距走时)、共中心点的定义，以及t和t0之间符合动校正公式的情况下，利用动校正速度即可计算出地震炮集数据对应的零偏移距信息。
步骤S106，利用平面波分解算法对叠后地震数据进行处理，得到地表射线参数。
本发明实施例不对平面波分解法进行具体限定，用户可以根据实际需求进行选择，例如，可使用希尔伯特变换算法。
步骤S108，基于零偏移距信息、动校正速度和地表射线参数构建地震炮集数据的局部倾角。
步骤S110，基于地震炮集数据的局部倾角对地震炮集数据进行去噪处理，得到待处理地理区域去噪后的地震炮集数据。
本发明实施例提供的方法使用了共反射面元走时理论，根据零偏移距信息、动校正速度信息和地表射线参数信息，可准确获取炮集数据中地震同相轴局部倾角信息，进而沿着地震炮集数据的局部倾角方向压制地震炮集数据中的噪音干扰。
本发明实施例所提供的地震炮集数据去噪方法不直接依赖地震炮集数据进行倾角计算，而是通过零偏移距信息、动校正速度和地表射线参数间接构建地震炮集数据的局部倾角场，动校正速度是地震数据处理中常用的速度，容易获取，且由于地表射线参数信息是基于信噪比较高的叠后数据中确定的，所以受噪音影响较小。因此，本发明提供的地震炮集数据去噪方法具有较强的抗噪性和较高的工业应用价值，可实现高质量的地震数据去噪效果。
在一个可选的实施方式中，如图2所示，上述步骤S104，计算地震炮集数据对应的零偏移距信息，具体包括如下步骤：
步骤S1041，获取地震炮集数据中的炮点位置信息、检波点位置信息和非零偏移距走时。
在获取到地震炮集数据U(xs,xr,t)之后，即可从中提取出炮点位置信息xs、检波点位置信息xr和非零偏移距走时t。
步骤S1042，基于炮点位置信息和检波点位置信息计算地震炮集数据对应的共中心点位置信息和偏移距信息。
已知共中心点表示炮点和检波点的中点，因此，根据炮点位置信息和检波点位置信息，即可计算出地震炮集数据对应的共中心点位置信息m和偏移距信息h，其中，
步骤S1043，基于非零偏移距走时、偏移距信息和动校正速度计算地震炮集数据对应的零偏移距走时。
具体的，动校正速度的公式为：
其中，V表示动校正速度V(m,t0)。因此，在已知动校正速度V(m,t0)，(xs,xr)对应的偏移距信息h和非零偏移距走时t的情况下，根据动校正公式即可求得地震炮集数据对应的零偏移距走时t0。
步骤S1044，基于共中心点位置信息和零偏移距走时确定零偏移距信息。
在一个可选的实施方式中，上述步骤S106，利用平面波分解算法对叠后地震数据进行处理，得到地表射线参数，具体包括如下步骤：
步骤S1061，对叠后地震数据进行平面波分解，得到叠后地震数据在空间方向的第一分量和在时间方向的第二分量。
步骤S1062，基于第一分量和第二分量计算地表射线参数。
已知地震平面波微分方程可表示为：
其中，M表示叠后地震数据M(m,t0)，p表示地表射线参数p(m,t0)，
其中，Hm(M)表示对叠后地震数据进行平面波分解之后，得到的叠后地震数据在空间方向的第一分量，Ht0(M)表示对叠后地震数据进行平面波分解之后，得到的叠后地震数据在时间方向的第二分量。因此，为了求出地表射线参数，若平面波分解算法采用希尔伯特变换，那么Hm(M)和Ht0(M)即对应的表示为对叠后地震数据进行希尔伯特变换在m和t0方向的分量。在得到第一分量和第二分量之后，根据上述算式，即可求得地表射线参数p(m,t0)。
在一个可选的实施方式中，上述步骤S108，基于零偏移距信息、动校正速度和地表射线参数构建地震炮集数据的局部倾角，具体包括如下步骤：
步骤S1081，基于零偏移距信息、动校正速度和地表射线参数，构建检波点移动目标距离的情况下，地震炮集数据的共反射面元走时方程。
步骤S1082，基于共反射面元走时方程计算地震炮集数据的局部倾角。
具体的，传统CRS走时方程为：
其中，TCRS表示共反射面元走时，Δm表示偏心距，偏心距是指共中心点m与参照共中心点m0之间的距离，也即，Δm＝m-m0，
其中。xr，xs分别是检波点和炮点信息。xr0，xs0分别是参照检波点和参照炮点信息。A，B，C为CRS走时计算的三个参数，A，B，C分别与地表射线参数，波前曲率和动校正速度有关。
接下来，对CRS走时方程进行数据域的变换，也即进行坐标转换，可获得CRS走时在炮集中的表达式：
该表达式的本质是用某个中心点周围的其他中心点对应的地震信息来提高该中心点对应信息的有效性。此时，该中心点即是参照中心点。
已知
m0表示(xr0，xs0)对应的共中心点，h0表示(xr0，xs0)对应的偏移距，当固定炮点，移动检波点(移动目标距离为Δx)时，可获得以下关系式：
也即，
因此，移动检波点后的走时，也即，地震炮集数据的共反射面元走时方程可表示为：
已知局部倾角与走时之间的关系，因此，在确定走时表达式的情况下，可推导出局部倾角算式。
在一个可选的实施方式中，上述步骤S1082，基于共反射面元走时方程计算地震炮集数据的局部倾角，具体包括如下步骤：
步骤S10821，对共反射面元走时方程进行求导数处理，得到求导方程。
步骤S10822，将求导方程中的目标距离设置0的结果作为地震炮集数据的局部倾角。
要计算地震炮集数据的局部倾角σ(xs,xr,t)，首先对地震炮集数据的共反射面元走时方程进行求导处理，得到：
根据以上求导方程可知，在xr0(也即，目标距离Δx＝0)处的导数即为当前位置(参照中心点m0)的局部倾角：
该算式中TCRS的取值即为参照中心点m0对应的非零偏移距走时t(xr0，xs0)，A＝2p，
其中，p表示地表射线参数p(m,t0)，V表示动校正速度V(m,t0)。以此类推，即可计算出地震炮集数据对应的所有共中心点的局部倾角，也即，地震炮集数据的局部倾角。
基于地下地质信息和地震波物理规律的去噪方法可以通过预测有效信号或干扰噪声来增强有效信号，其中，共反射面(CRS)叠加技术可以对地震多覆盖反射数据进行零偏移模拟。基于几何地震学，CRS叠加利用第一菲涅耳带的所有反射来提高信噪比。然而，传统的CRS参数估计繁琐、耗时且成本高昂。与现有技术相比，本发明实施例使用CRS参数来估计局部倾角，将所需的CRS参数数量从3个(A，B，C)减少到2个(A，C)，有效提高了计算效率。
在一个可选的实施方式中，上述步骤S110，基于地震炮集数据的局部倾角对地震炮集数据进行去噪处理，具体包括如下内容：
沿着地震炮集数据的局部倾角方向对地震炮集数据进行均值滤波处理，得到待处理地理区域去噪后的地震炮集数据。
可选地，沿着地震炮集数据的局部倾角方向对地震炮集数据进行均值滤波处理，具体包括：基于算式
对地震炮集数据进行均值滤波，得到待处理地理区域去噪后的地震炮集数据；其中，R(xs,xr,t)表示去噪后的地震炮集数据，U(xs,xr,t)表示地震炮集数据，xs表示地震炮集数据中的炮点位置信息，xr表示地震炮集数据中的检波点位置信息，t表示地震炮集数据中与(xs,xr)对应的非零偏移距走时，n表示均值滤波窗口的地震道数量，也即均值滤波的滤波窗口尺度，用户可以根据实际需求对n值进行设置，本发明实施例不对其进行具体限定。Δxi表示当前检波点位置xr相对第i道地震数据的位移，根据地震数据的道头信息即可获取；Δti表示当前检波点位置xr相对第i道地震数据的时移，且
也即，上述位移与时移均是沿着倾角σ(xs,xr,t)方向的时移量和位移量，σ表示地震炮集数据的局部倾角σ(xs,xr,t)。
针对地震事件在炮域振幅连续性大的特点，本发明实施例提供的地震炮集数据去噪方法实质是一种基于炮集局部倾角的去噪方法，与常规倾角计算方法不同，本发明实施例提出了炮集局部倾角构建的一种新思路，利用动校正速度信息和地表射线参数信息，可准确获取炮集数据中地震同相轴局部倾角信息，为地震去噪提供准确的方向信息，既能增强有效事件，又能抑制非相干和相干噪声等干扰噪声，实现地震数据的噪音压制，提高数据信噪比。基于共反射面(CRS)走时，推导和估计地震事件在炮域的局部斜率，为平面波预测提供结构信息。与传统的全CRS行程时间相比，本发明实施例提出的基于CRS求解局部倾角的方法依赖更少的参数，使其计算效率更高。
实施例二
本发明实施例还提供了一种地震炮集数据去噪装置，该地震炮集数据去噪装置主要用于执行上述实施例一所提供的地震炮集数据去噪方法，以下对本发明实施例提供的地震炮集数据去噪装置做具体介绍。
图3是本发明实施例提供的一种地震炮集数据去噪装置的功能模块图，如图3所示，该装置主要包括：获取模块10，计算模块20，分解模块30，构建模块40，去噪模块50，其中：
获取模块10，用于获取待处理地理区域的地震炮集数据、动校正速度和叠后地震数据。
计算模块20，用于基于动校正速度计算地震炮集数据对应的零偏移距信息。
分解模块30，用于利用平面波分解算法对叠后地震数据进行处理，得到地表射线参数。
构建模块40，用于基于零偏移距信息、动校正速度和地表射线参数构建地震炮集数据的局部倾角。
去噪模块50，用于基于地震炮集数据的局部倾角对地震炮集数据进行去噪处理，得到待处理地理区域去噪后的地震炮集数据。
本发明实施例所提供的地震炮集数据去噪装置，不直接依赖地震炮集数据进行倾角计算，而是通过零偏移距信息、动校正速度和地表射线参数间接构建地震炮集数据的局部倾角场，由于地表射线参数信息是基于信噪比较高的叠后数据中确定的，所以受噪音影响较小。因此，本发明提供的地震炮集数据去噪装置具有较强的抗噪性和较高的工业应用价值，可实现高质量的地震数据去噪效果。
可选地，计算模块20具体用于：
获取地震炮集数据中的炮点位置信息、检波点位置信息和非零偏移距走时。
基于炮点位置信息和检波点位置信息计算地震炮集数据对应的共中心点位置信息和偏移距信息。
基于非零偏移距走时、偏移距信息和动校正速度计算地震炮集数据对应的零偏移距走时。
基于共中心点位置信息和零偏移距走时确定零偏移距信息。
可选地，分解模块30具体用于：
对叠后地震数据进行平面波分解，得到叠后地震数据在空间方向的第一分量和在时间方向的第二分量。
基于第一分量和第二分量计算地表射线参数。
可选地，构建模块40包括：
构建单元，用于基于零偏移距信息、动校正速度和地表射线参数，构建检波点移动目标距离的情况下，地震炮集数据的共反射面元走时方程。
计算单元，用于基于共反射面元走时方程计算地震炮集数据的局部倾角。
可选地，计算单元具体用于：
对共反射面元走时方程进行求导数处理，得到求导方程。
将求导方程中的目标距离设置0的结果作为地震炮集数据的局部倾角。
可选地，去噪模块50包括：
滤波单元，用于沿着地震炮集数据的局部倾角方向对地震炮集数据进行均值滤波处理，得到待处理地理区域去噪后的地震炮集数据。
可选地，滤波单元具体用于：
基于算式
对地震炮集数据进行均值滤波，得到待处理地理区域去噪后的地震炮集数据；其中，R(xs,xr,t)表示去噪后的地震炮集数据，U(xs,xr,t)表示地震炮集数据，xs表示地震炮集数据中的炮点位置信息，xr表示地震炮集数据中的检波点位置信息，t表示地震炮集数据中与(xs,xr)对应的非零偏移距走时，n表示均值滤波窗口的地震道数量，Δxi表示当前检波点位置xr相对第i道地震数据的位移，Δti表示当前检波点位置xr相对第i道地震数据的时移，σ表示地震炮集数据的局部倾角。
实施例三
参见图4，本发明实施例提供了一种电子设备，该电子设备包括：处理器60，存储器61，总线62和通信接口63，所述处理器60、通信接口63和存储器61通过总线62连接；处理器60用于执行存储器61中存储的可执行模块，例如计算机程序。
其中，存储器61可能包含高速随机存取存储器(RAM，Random Access Memory)，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。通过至少一个通信接口63(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。
总线62可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图4中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
其中，存储器61用于存储程序，所述处理器60在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的过程定义的装置所执行的方法可以应用于处理器60中，或者由处理器60实现。
处理器60可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器60中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器60可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(Digital SignalProcessing，简称DSP)、专用集成电路(Application Specific Integrated Circuit，简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器61，处理器60读取存储器61中的信息，结合其硬件完成上述方法的步骤。
本发明实施例所提供的一种地震炮集数据去噪方法和装置的计算机程序产品，包括存储了处理器可执行的非易失的程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，或者是该发明产品使用时惯常摆放的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
此外，术语“水平”、“竖直”、“悬垂”等术语并不表示要求部件绝对水平或悬垂，而是可以稍微倾斜。如“水平”仅仅是指其方向相对“竖直”而言更加水平，并不是表示该结构一定要完全水平，而是可以稍微倾斜。
在本发明的描述中，还需要说明的是，除非另有明确的规定和限定，术语“设置”、“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种地震炮集数据去噪方法，其特征在于，包括：
获取待处理地理区域的地震炮集数据、动校正速度和叠后地震数据；
基于所述动校正速度计算所述地震炮集数据对应的零偏移距信息；
利用平面波分解算法对所述叠后地震数据进行处理，得到地表射线参数；
基于所述零偏移距信息、所述动校正速度和所述地表射线参数构建所述地震炮集数据的局部倾角；
基于所述地震炮集数据的局部倾角对所述地震炮集数据进行去噪处理，得到所述待处理地理区域去噪后的地震炮集数据。
2.根据权利要求1所述的地震炮集数据去噪方法，其特征在于，计算所述地震炮集数据对应的零偏移距信息，包括：
获取所述地震炮集数据中的炮点位置信息、检波点位置信息和非零偏移距走时；
基于所述炮点位置信息和检波点位置信息计算所述地震炮集数据对应的共中心点位置信息和偏移距信息；
基于所述非零偏移距走时、所述偏移距信息和所述动校正速度计算所述地震炮集数据对应的零偏移距走时；
基于所述共中心点位置信息和所述零偏移距走时确定所述零偏移距信息。
3.根据权利要求1所述的地震炮集数据去噪方法，其特征在于，利用平面波分解算法对所述叠后地震数据进行处理，得到地表射线参数，包括：
对所述叠后地震数据进行平面波分解，得到所述叠后地震数据在空间方向的第一分量和在时间方向的第二分量；
基于所述第一分量和所述第二分量计算所述地表射线参数。
4.根据权利要求1所述的地震炮集数据去噪方法，其特征在于，基于所述零偏移距信息、所述动校正速度和所述地表射线参数构建所述地震炮集数据的局部倾角，包括：
基于所述零偏移距信息、所述动校正速度和所述地表射线参数，构建检波点移动目标距离的情况下，地震炮集数据的共反射面元走时方程；
基于所述共反射面元走时方程计算所述地震炮集数据的局部倾角。
5.根据权利要求4所述的地震炮集数据去噪方法，其特征在于，基于所述共反射面元走时方程计算所述地震炮集数据的局部倾角，包括：
对所述共反射面元走时方程进行求导数处理，得到求导方程；
将求导方程中的所述目标距离设置0的结果作为所述地震炮集数据的局部倾角。
6.根据权利要求1所述的地震炮集数据去噪方法，其特征在于，基于所述地震炮集数据的局部倾角对所述地震炮集数据进行去噪处理，包括：
沿着所述地震炮集数据的局部倾角方向对所述地震炮集数据进行均值滤波处理，得到所述待处理地理区域去噪后的地震炮集数据。
7.根据权利要求6所述的地震炮集数据去噪方法，其特征在于，沿着所述地震炮集数据的局部倾角方向对所述地震炮集数据进行均值滤波处理，包括：
基于算式
对所述地震炮集数据进行均值滤波，得到所述待处理地理区域去噪后的地震炮集数据；其中，R(xs,xr,t)表示所述去噪后的地震炮集数据，U(xs,xr,t)表示所述地震炮集数据，xs表示所述地震炮集数据中的炮点位置信息，xr表示所述地震炮集数据中的检波点位置信息，t表示所述地震炮集数据中与(xs,xr)对应的非零偏移距走时，n表示均值滤波窗口的地震道数量，Δxi表示当前检波点位置xr相对第i道地震数据的位移，Δti表示当前检波点位置xr相对第i道地震数据的时移，σ表示所述地震炮集数据的局部倾角。
8.一种地震炮集数据去噪装置，其特征在于，包括：
获取模块，用于获取待处理地理区域的地震炮集数据、动校正速度和叠后地震数据；
计算模块，用于基于所述动校正速度计算所述地震炮集数据对应的零偏移距信息；
分解模块，用于利用平面波分解算法对所述叠后地震数据进行处理，得到地表射线参数；
构建模块，用于基于所述零偏移距信息、所述动校正速度和所述地表射线参数构建所述地震炮集数据的局部倾角；
去噪模块，用于基于所述地震炮集数据的局部倾角对所述地震炮集数据进行去噪处理，得到所述待处理地理区域去噪后的地震炮集数据。
9.一种电子设备，包括存储器、处理器，所述存储器上存储有可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至7中任一项所述的地震炮集数据去噪方法的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令被处理器执行时实现上述权利要求1至7中任一项所述的地震炮集数据去噪方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211366969.1A CN115685320B (zh) | 2022-11-02 | 2022-11-02 | 一种地震炮集数据去噪方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211366969.1A CN115685320B (zh) | 2022-11-02 | 2022-11-02 | 一种地震炮集数据去噪方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115685320A true CN115685320A (zh) | 2023-02-03 |
| CN115685320B CN115685320B (zh) | 2024-01-26 |
Family
ID=85048054
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211366969.1A Active CN115685320B (zh) | 2022-11-02 | 2022-11-02 | 一种地震炮集数据去噪方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115685320B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105242309A (zh) * | 2015-10-19 | 2016-01-13 | 中国石油天然气集团公司 | 一种地震数据中规则干扰的压制方法及装置 |
| CN106772583A (zh) * | 2017-01-10 | 2017-05-31 | 中国科学院地质与地球物理研究所 | 一种地震绕射波分离方法和装置 |
| CN107817523A (zh) * | 2017-10-30 | 2018-03-20 | 中国矿业大学(北京) | 绕射波偏移速度的分析方法及装置 |
| CN113640872A (zh) * | 2021-08-12 | 2021-11-12 | 中国矿业大学(北京) | 绕射波分离方法、装置和电子设备 |
-
2022
- 2022-11-02 CN CN202211366969.1A patent/CN115685320B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105242309A (zh) * | 2015-10-19 | 2016-01-13 | 中国石油天然气集团公司 | 一种地震数据中规则干扰的压制方法及装置 |
| CN106772583A (zh) * | 2017-01-10 | 2017-05-31 | 中国科学院地质与地球物理研究所 | 一种地震绕射波分离方法和装置 |
| WO2018129844A1 (zh) * | 2017-01-10 | 2018-07-19 | 中国科学院地质与地球物理研究所 | 一种地震绕射波分离方法和装置 |
| CN107817523A (zh) * | 2017-10-30 | 2018-03-20 | 中国矿业大学(北京) | 绕射波偏移速度的分析方法及装置 |
| CN113640872A (zh) * | 2021-08-12 | 2021-11-12 | 中国矿业大学(北京) | 绕射波分离方法、装置和电子设备 |
Non-Patent Citations (1)
| Title |
|---|
| 李正伟;张剑锋;刘伟;: "基于倾角-偏移距域道集的绕射波成像", 地球物理学报, vol. 61, no. 04, pages 1447 - 1459 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115685320B (zh) | 2024-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7295930B2 (en) | Method, device and software package for smoothing a subsurface property | |
| RU2300123C2 (ru) | Высокоразрешающее преобразование радона для обработки сейсмических данных | |
| CN113640872B (zh) | 绕射波分离方法、装置和电子设备 | |
| CN111929729B (zh) | 绕射波成像方法、装置和电子设备 | |
| EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
| CN113945969B (zh) | 不连续地质体检测方法、装置及电子设备 | |
| CN109116413B (zh) | 成像域立体层析速度反演方法 | |
| CN110297273B (zh) | 储层的叠前反演方法及装置 | |
| CN112099082B (zh) | 一种共面元共方位角道集的地震回折波走时反演方法 | |
| CN115685320A (zh) | 一种地震炮集数据去噪方法和装置 | |
| CN107451992B (zh) | 一种sar图像变化检测的方法与装置 | |
| US9442207B1 (en) | System and method for computing residual moveout from seismic images | |
| CN109387872B (zh) | 表面多次波预测方法 | |
| CN113945987B (zh) | 病害地质体检测的方法、装置及电子设备 | |
| CN111929730B (zh) | 小尺度地质异常体检测方法及装置 | |
| CN112986990B (zh) | 一种大气相位改正方法及系统 | |
| WO2014176116A1 (en) | Surface multiple prediction | |
| CN113608261B (zh) | 绕射波成像方法、装置和电子设备 | |
| CN111399038B (zh) | 斜率参数提取方法、装置及计算机可读存储介质 | |
| CN110873898B (zh) | 地震数据的滤波方法及装置 | |
| Waldeland et al. | Fast estimation of prestack common reflection surface parameters | |
| US9014985B2 (en) | System and method for compensating time and offset varying near-surface effects in seismic data background | |
| CN115903043A (zh) | 一种绕射波分离方法和装置 | |
| CN116381786B (zh) | 绕射波成像方法和装置 | |
| CN112379420B (zh) | 高精度弯曲测线叠前时间域成像方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |