CN112782764B - 线性干扰衰减方法及系统 - Google Patents
线性干扰衰减方法及系统 Download PDFInfo
- Publication number
- CN112782764B CN112782764B CN201911085616.2A CN201911085616A CN112782764B CN 112782764 B CN112782764 B CN 112782764B CN 201911085616 A CN201911085616 A CN 201911085616A CN 112782764 B CN112782764 B CN 112782764B
- Authority
- CN
- China
- Prior art keywords
- ray
- linear
- ray parameter
- gather
- data volume
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 229910052704 radon Inorganic materials 0.000 claims abstract description 56
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 56
- 238000012545 processing Methods 0.000 claims abstract description 26
- 238000004590 computer program Methods 0.000 claims description 15
- 238000009499 grossing Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 8
- 239000000523 sample Substances 0.000 description 14
- 238000002474 experimental method Methods 0.000 description 7
- 230000002238 attenuated effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000002547 anomalous effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种线性干扰衰减方法及系统。该线性干扰衰减方法包括:对原始地震数据中的每个道集进行分频处理;对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅;确定每个时间样点的多个振幅中的最大值对应的射线参数;根据每个道集中的每个频带的射线参数形成第一射线参数数据体;获取其中一个道集的线性干扰速度和该道集对应的射线参数;根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体;根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。本发明可以提高线性干扰去噪的效率和精度。
Description
技术领域
本发明涉及地震数据处理技术领域,具体地,涉及一种线性干扰衰减方法及系统。
背景技术
线性干扰是地震原始道集中的常见干扰,对线性干扰进行噪声衰减是地震数据处理过程中的常规步骤,线性干扰衰减需要提前获取线性干扰的频带和速度范围参数。
为了有效地衰减线性干扰,处理人员需要在工区范围内在不同的位置人工统计线性干扰的特征,统计过程效率低、误差大,受制于时间和人力限制也无法对全工区每个道集进行统计。
发明内容
本发明实施例的主要目的在于提供一种线性干扰衰减方法及系统,以提高线性干扰去噪的效率和精度。
为了实现上述目的,本发明实施例提供一种线性干扰衰减方法,包括:
对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;原始地震数据包括多个道集;
对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅;
确定每个时间样点的多个振幅中的最大值对应的射线参数;
根据每个道集中的每个频带的射线参数形成第一射线参数数据体;
获取其中一个道集的线性干扰速度和该道集对应的射线参数;
根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体;
根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。
本发明实施例还提供一种线性干扰衰减系统,包括:
分频单元,用于对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;原始地震数据包括多个道集;
线性拉东变换单元,用于对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅;
射线参数确定单元,用于确定每个时间样点的多个振幅中的最大值对应的射线参数;
第一射线参数数据体单元,用于根据每个道集中的每个频带的射线参数形成第一射线参数数据体;
射线参数获取单元,用于获取其中一个道集的线性干扰速度和该道集对应的射线参数;
第二射线参数数据体单元,用于根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体;
线性干扰衰减单元,用于根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现所述的线性干扰衰减方法的步骤。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现所述的线性干扰衰减方法的步骤。
本发明实施例的线性干扰衰减方法及系统先对原始地震数据中的每个道集进行分频处理,再对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果,接着确定每个时间样点的多个振幅中的最大值对应的射线参数,根据每个道集中的每个频带的射线参数形成第一射线参数数据体;然后获取其中一个道集的线性干扰速度和该道集对应的射线参数,并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体,最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据,无需参数实验且易于实施,可以提高线性干扰去噪的效率和精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例中线性干扰衰减方法的流程图;
图2是本发明实施例中异常能量干扰衰减前的不同时窗及属性平面图;
图3是本发明实施例中异常能量干扰衰减后的不同时窗及属性平面图;
图4是本发明实施例中线性干扰衰减系统的结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本领域技术人员知道,本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此,本公开可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
鉴于现有技术的统计过程效率低、误差大,受制于时间和人力限制也无法对全工区每个道集进行统计,本发明实施例提供一种线性干扰衰减方法,无需参数实验且易于实施,可以提高线性干扰去噪的效率和精度。以下结合附图对本发明进行详细说明。
图1是本发明实施例中线性干扰衰减方法的流程图。如图1所示,线性干扰衰减方法包括:
S101:对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;原始地震数据包括多个道集。
其中,道集可以为炮集、检波点道集、共中心点道集或十字排列道集,道集可以分为左右排列分别处理。频带宽度可以为10HZ-20HZ。
S102:对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅。
其中,线性拉东变换结果中的射线参数大于原始地震数据中的最大射线参数,或线性拉东变换结果中的射线参数小于原始地震数据中的最小射线参数。线性拉东变换结果的纵坐标为时间,横坐标为射线参数。
S103:确定每个时间样点的多个振幅中的最大值对应的射线参数。
S104:根据每个道集中的每个频带的射线参数形成第一射线参数数据体。
S105:获取其中一个道集的线性干扰速度和该道集对应的射线参数。
S106:根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体。
其中,在执行S106之前,还可以包括:对第一射线参数数据体沿空间方向进行平滑处理以消除部分异常结果。平滑半径可以为道集位置横向间隔的3-5倍。
此时S106包括:根据线性干扰速度、该道集对应的射线参数和经过平滑处理的第一射线参数数据体,获得第二射线参数数据体。
具体实施时,第二射线参数数据体为经过平滑处理的第一射线参数数据体与比例系数的积;比例系数为线性干扰速度与该道集对应的射线参数的商。即,线性干扰速度除以该道集对应的射线参数,乘以经过平滑处理的第一射线参数数据体,即可得到第二射线参数数据体。
S107:根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。
图1所示的线性干扰衰减方法的执行主体可以为计算机。由图1所示的流程可知,本发明实施例的线性干扰衰减方法先对原始地震数据中的每个道集进行分频处理,再对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果,接着确定每个时间样点的多个振幅中的最大值对应的射线参数,根据每个道集中的每个频带的射线参数形成第一射线参数数据体;然后获取其中一个道集的线性干扰速度和该道集对应的射线参数,并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体,最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据,无需参数实验且易于实施,可以提高线性干扰去噪的效率和精度。
本发明实施例的具体流程如下:
1、对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;原始地震数据包括多个道集。
2、对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅。线性拉东变换结果中的射线参数大于原始地震数据中的最大射线参数,或线性拉东变换结果中的射线参数小于原始地震数据中的最小射线参数。
3、确定每个时间样点的多个振幅中的最大值对应的射线参数,根据每个道集中的每个频带的射线参数形成第一射线参数数据体。
4、对第一射线参数数据体沿空间方向进行平滑处理以消除部分异常结果。平滑半径可以为道集位置横向间隔的3-5倍。
5、获取其中一个道集的线性干扰速度和该道集对应的射线参数。
6、令线性干扰速度除以该道集对应的射线参数,并乘以经过平滑处理的第一射线参数数据体,得到第二射线参数数据体。
7、根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。
图2是本发明实施例中线性干扰衰减前的不同时窗及射线参数属性平面图;图3是本发明实施例中线性干扰衰减后的不同时窗及射线参数属性平面图。由图2可以看出在线性干扰衰减前,线性干扰的倍数分布不均匀且横向差异较大,如果通过人工调查来获得线性干扰视速度数据,不但工作量大且精度低。由图3可以看出经过线性干扰衰减后,地震数据射线参数属性横向趋于一致,说明地震数据中在线性干扰视速度范围内已经没有明显的线性信号,达到了线性干扰衰减的目的。
综上,本发明实施例的线性干扰衰减方法先对原始地震数据中的每个道集进行分频处理,再对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果,接着确定每个时间样点的多个振幅中的最大值对应的射线参数,根据每个道集中的每个频带的射线参数形成第一射线参数数据体;然后获取其中一个道集的线性干扰速度和该道集对应的射线参数,并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体,最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据,无需参数实验且易于实施,可以提高线性干扰去噪的效率和精度。
基于同一发明构思,本发明实施例还提供了一种线性干扰衰减系统,由于该系统解决问题的原理与线性干扰衰减方法相似,因此该系统的实施可以参见方法的实施,重复之处不再赘述。
图4是本发明实施例中线性干扰衰减系统的结构框图。如图4所示,线性干扰衰减系统包括:
分频单元,用于对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;原始地震数据包括多个道集;
线性拉东变换单元,用于对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅;
射线参数确定单元,用于确定每个时间样点的多个振幅中的最大值对应的射线参数;
第一射线参数数据体单元,用于根据每个道集中的每个频带的射线参数形成第一射线参数数据体;
射线参数获取单元,用于获取其中一个道集的线性干扰速度和该道集对应的射线参数;
第二射线参数数据体单元,用于根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体;
线性干扰衰减单元,用于根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。
在其中一种实施例中,线性拉东变换结果中的射线参数大于原始地震数据中的最大射线参数,或线性拉东变换结果中的射线参数小于原始地震数据中的最小射线参数。
在其中一种实施例中,还包括:
平滑单元,用于对第一射线参数数据体进行平滑处理;
第二射线参数数据体单元具体用于:
根据线性干扰速度、该道集对应的射线参数和经过平滑处理的第一射线参数数据体,获得第二射线参数数据体。
在其中一种实施例中,道集为炮集、检波点道集、共中心点道集或十字排列道集。
综上,本发明实施例的线性干扰衰减系统先对原始地震数据中的每个道集进行分频处理,再对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果,接着确定每个时间样点的多个振幅中的最大值对应的射线参数,根据每个道集中的每个频带的射线参数形成第一射线参数数据体;然后获取其中一个道集的线性干扰速度和该道集对应的射线参数,并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体,最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据,无需参数实验且易于实施,可以提高线性干扰去噪的效率和精度。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时可以实现基于线性干扰衰减方法的全部或部分内容,例如,处理器执行计算机程序时可以实现如下内容:
对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;原始地震数据包括多个道集;
对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅;
确定每个时间样点的多个振幅中的最大值对应的射线参数;
根据每个道集中的每个频带的射线参数形成第一射线参数数据体;
获取其中一个道集的线性干扰速度和该道集对应的射线参数;
根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体;
根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。
综上,本发明实施例的计算机设备先对原始地震数据中的每个道集进行分频处理,再对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果,接着确定每个时间样点的多个振幅中的最大值对应的射线参数,根据每个道集中的每个频带的射线参数形成第一射线参数数据体;然后获取其中一个道集的线性干扰速度和该道集对应的射线参数,并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体,最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据,无需参数实验且易于实施,可以提高线性干扰去噪的效率和精度。
本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时可以实现线性干扰衰减方法的全部或部分内容,例如,处理器执行计算机程序时可以实现如下内容:
对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;原始地震数据包括多个道集;
对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅;
确定每个时间样点的多个振幅中的最大值对应的射线参数;
根据每个道集中的每个频带的射线参数形成第一射线参数数据体;
获取其中一个道集的线性干扰速度和该道集对应的射线参数;
根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体;
根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据。
综上,本发明实施例的计算机可读存储介质先对原始地震数据中的每个道集进行分频处理,再对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果,接着确定每个时间样点的多个振幅中的最大值对应的射线参数,根据每个道集中的每个频带的射线参数形成第一射线参数数据体;然后获取其中一个道集的线性干扰速度和该道集对应的射线参数,并根据线性干扰速度、该道集对应的射线参数和第一射线参数数据体,获得第二射线参数数据体,最后根据第二射线参数数据体和原始地震数据获得经过线性干扰衰减的地震数据,无需参数实验且易于实施,可以提高线性干扰去噪的效率和精度。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrative logical block),单元,和步骤可以通过电子硬件、电脑软件,或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability),上述的各种说明性部件(illustrative components),单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块,或单元,或装置都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
Claims (8)
1.一种线性干扰衰减方法,其特征在于,包括:
对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;所述原始地震数据包括多个道集;
对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;所述线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅;
确定每个时间样点的多个振幅中的最大值对应的射线参数;
根据每个道集中的每个频带的射线参数形成第一射线参数数据体;
获取其中一个道集的线性干扰速度和该道集对应的射线参数;
根据所述线性干扰速度、该道集对应的射线参数和所述第一射线参数数据体,获得第二射线参数数据体;
根据所述第二射线参数数据体和所述原始地震数据获得经过线性干扰衰减的地震数据;
所述线性干扰衰减方法还包括:
对所述第一射线参数数据体进行平滑处理;
获得第二射线参数数据体包括:
所述线性干扰速度除以该道集对应的射线参数,乘以经过平滑处理的第一射线参数数据体,获得第二射线参数数据体。
2.根据权利要求1所述的线性干扰衰减方法,其特征在于,所述线性拉东变换结果中的射线参数大于所述原始地震数据中的最大射线参数,或所述线性拉东变换结果中的射线参数小于所述原始地震数据中的最小射线参数。
3.根据权利要求1所述的线性干扰衰减方法,其特征在于,所述道集为炮集、检波点道集、共中心点道集或十字排列道集。
4.一种线性干扰衰减系统,其特征在于,包括:
分频单元,用于对原始地震数据中的每个道集进行分频处理,每个道集均包括多个频带;所述原始地震数据包括多个道集;
线性拉东变换单元,用于对每个频带的原始地震数据进行线性拉东变换,得到每个频带的线性拉东变换结果;所述线性拉东变换结果包括多个时间样点,每个时间样点对应多个振幅;
射线参数确定单元,用于确定每个时间样点的多个振幅中的最大值对应的射线参数;
第一射线参数数据体单元,用于根据每个道集中的每个频带的射线参数形成第一射线参数数据体;
射线参数获取单元,用于获取其中一个道集的线性干扰速度和该道集对应的射线参数;
第二射线参数数据体单元,用于根据所述线性干扰速度、该道集对应的射线参数和所述第一射线参数数据体,获得第二射线参数数据体;
线性干扰衰减单元,用于根据所述第二射线参数数据体和所述原始地震数据获得经过线性干扰衰减的地震数据;
所述线性干扰衰减系统还包括:
平滑单元,用于对所述第一射线参数数据体进行平滑处理;
所述第二射线参数数据体单元具体用于:
所述线性干扰速度除以该道集对应的射线参数,乘以经过平滑处理的第一射线参数数据体,获得第二射线参数数据体。
5.根据权利要求4所述的线性干扰衰减系统,其特征在于,所述线性拉东变换结果中的射线参数大于所述原始地震数据中的最大射线参数,或所述线性拉东变换结果中的射线参数小于所述原始地震数据中的最小射线参数。
6.根据权利要求4所述的线性干扰衰减系统,其特征在于,所述道集为炮集、检波点道集、共中心点道集或十字排列道集。
7.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至3任一项所述的线性干扰衰减方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至3任一项所述的线性干扰衰减方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911085616.2A CN112782764B (zh) | 2019-11-08 | 2019-11-08 | 线性干扰衰减方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911085616.2A CN112782764B (zh) | 2019-11-08 | 2019-11-08 | 线性干扰衰减方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112782764A CN112782764A (zh) | 2021-05-11 |
CN112782764B true CN112782764B (zh) | 2022-10-04 |
Family
ID=75748135
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911085616.2A Active CN112782764B (zh) | 2019-11-08 | 2019-11-08 | 线性干扰衰减方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112782764B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5138583A (en) * | 1991-05-22 | 1992-08-11 | Halliburton Geophysical Services, Inc. | Noise attenuation method |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7065004B2 (en) * | 2004-04-22 | 2006-06-20 | Pgs Americas, Inc. | Horizon-based residual depth migration velocity analysis |
US7561491B2 (en) * | 2005-03-04 | 2009-07-14 | Robinson John M | Radon transformations for removal of noise from seismic data |
CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
US20110141849A1 (en) * | 2009-12-16 | 2011-06-16 | John Brittan | Method for attenuating interference noise in dual-sensor seismic data |
CN102141634A (zh) * | 2010-12-17 | 2011-08-03 | 西南交通大学 | 基于曲波变换的叠前地震信号中线性干扰的压制方法 |
CN102681014B (zh) * | 2012-05-23 | 2014-08-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于多项式拟合的规则线性干扰压制方法 |
CN104808245A (zh) * | 2015-05-19 | 2015-07-29 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 道集优化处理方法及其装置 |
-
2019
- 2019-11-08 CN CN201911085616.2A patent/CN112782764B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5138583A (en) * | 1991-05-22 | 1992-08-11 | Halliburton Geophysical Services, Inc. | Noise attenuation method |
Also Published As
Publication number | Publication date |
---|---|
CN112782764A (zh) | 2021-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3171363B1 (en) | Voice activity detection methods and apparatuses | |
CN107748734B (zh) | 一种解析-经验模态分解方法 | |
CN112415078B (zh) | 一种质谱数据谱图信号校准方法及装置 | |
CN107070940B (zh) | 一种从流式登录日志中判断恶意登录ip地址的方法及装置 | |
CN111583132B (zh) | 一种去除遥感影像异常条带式噪声的方法、装置、设备及介质 | |
CN112782764B (zh) | 线性干扰衰减方法及系统 | |
CN111460934B (zh) | 一种超导高场磁共振的中心频率计算方法及系统 | |
CN104133248A (zh) | 一种高保真声波干扰压制方法 | |
CN112183407B (zh) | 一种基于时频域谱减法的隧道地震波数据去噪方法及系统 | |
CN108040315A (zh) | 一种计算机可读存储介质和应用该介质的测试机 | |
CN116660176A (zh) | 傅里叶光谱自动基线校正方法、装置和存储介质 | |
CN113375065B (zh) | 管道泄漏监测中趋势信号的消除方法及装置 | |
CN112782765B (zh) | 异常能量干扰衰减方法及系统 | |
CN112674776B (zh) | 一种基于分段多项式拟合ecg信号的平滑滤波方法 | |
CN112014884B (zh) | 压制近炮点强能量噪声的方法及装置 | |
CN109726667B (zh) | 质谱数据处理方法和装置、计算机设备、计算机存储介质 | |
CN109446572B (zh) | 一种量子测控数据的处理方法及装置 | |
CN111353860A (zh) | 产品信息推送方法及系统 | |
CN113900147B (zh) | 沙丘鸣震压制方法及系统 | |
CN113900146A (zh) | 面波压制方法及系统 | |
CN109857973B (zh) | 红外数据滤波方法及装置 | |
CN113050167B (zh) | 地震数据的噪声压制方法及装置 | |
CN116629028B (zh) | 一种花瓣形发射台导流槽参数的确定方法及装置 | |
CN111291862A (zh) | 用于模型压缩的方法和装置 | |
CN110470396B (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 |