CN111766631B - 一种面波压制方法及装置 - Google Patents
一种面波压制方法及装置 Download PDFInfo
- Publication number
- CN111766631B CN111766631B CN201910255825.0A CN201910255825A CN111766631B CN 111766631 B CN111766631 B CN 111766631B CN 201910255825 A CN201910255825 A CN 201910255825A CN 111766631 B CN111766631 B CN 111766631B
- Authority
- CN
- China
- Prior art keywords
- surface wave
- shot
- super
- determining
- frequency domain
- 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
- 238000003825 pressing Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 53
- 230000001629 suppression Effects 0.000 claims abstract description 21
- 230000001131 transforming effect Effects 0.000 claims abstract description 9
- 238000001514 detection method Methods 0.000 claims description 18
- 238000009826 distribution Methods 0.000 claims description 15
- 238000001228 spectrum Methods 0.000 claims description 9
- 238000005520 cutting process Methods 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 description 7
- 239000000523 sample Substances 0.000 description 7
- 238000005070 sampling Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 239000010410 layer Substances 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005055 memory storage 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
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 239000002344 surface layer Substances 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/30—Noise handling
- G01V2210/32—Noise reduction
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)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本说明书实施例提供一种面波压制方法及装置。所述方法包括:对采集到的地震数据中的共炮点道集进行分组得到超炮集;确定所述超炮集中的面波特征信息;根据所述面波特征信息,在变换至频域的超炮集中确定面波频域能量;将所述面波频域能量变换为面波时域能量;利用所述面波特征信息,在所述面波时域能量中确定面波;从所述超炮集中去除面波。利用本申请实施例,能够实现对于地震波中面波的压制,提高地震资料的信噪比。
Description
技术领域
本说明书实施例涉及地质勘探技术领域,特别涉及一种面波压制方法及装置。
背景技术
地震学的主要内容之一是研究地震波带来的信息。而地震波中又存在着很多与有效信号重叠的干扰信号。其中一种沿自由表面或岩层分界面传播的波,其能量主要限制在一个表面薄层内,被称为面波。面波是一种主要的相干噪声,具有频率低、能量强、影响范围广的特点。它产生于近地表,是由瑞雷波的垂直分量组成。地震波在地下传播,由于地层对地震体波的吸收衰减作用使得地震信号的高频成分随传播距离的增加迅速衰减,导致地震有效信号的频率降低而能量减小。在中深层,有效波频带因高频吸收作用与面波的频率重叠严重。而且面波能量衰减较为缓慢,与中深层的有效信号相比,通常能量是有效波能量的几倍甚至几十倍。因此,面波干扰是地震记录中深层信噪比低的重要原因。
根据面波的特征,目前已经提出了许多种压制面波的方法,如低截频滤波、区域滤波、内切除滤波、F-K滤波,Radon变换等。利用其中的某些方法可以实现根据面波的局部线性特征进行面波衰减。
但是根据面波的线性特征进行面波衰减对于面波同向轴的线性规律要求较高,目前普遍应用的三维地震勘探中面波同向轴为曲线型。曲线型的面波在曲率较大的位置,利用线性倾角扫描的方法识别比较困难,面波的衰减效果较差。以上所述的Radon变换、F-K滤波在处理这些问题时都存在不足之处。同时,由于空间采样不足,上述方法在分离面波的时候可能会出现假频影响,从而干扰最终的面波压制结果。因此,目前亟需一种能够有效压制面波的方法。
发明内容
本说明书实施例的目的是提供一种面波压制方法及装置,以解决现有技术中压制面波时准确率较差,面波衰减效果较差的问题。
为解决上述问题,本申请实施例提供一种面波压制的方法及装置是这样实现的:
一种面波压制方法,包括:
对采集到的地震数据中的共炮点道集进行分组得到超炮集;
确定所述超炮集中的面波特征信息;
根据所述面波特征信息,在变换至频域的超炮集中确定面波频域能量;
将所述面波频域能量变换为面波时域能量;
利用所述面波特征信息,在所述面波时域能量中确定面波;
去除所述超炮集中的面波。
一种面波压制装置,包括:
超炮集获取模块,用于对采集到的地震数据中的共炮点道集进行分组得到超炮集;
面波特征信息确定模块,用于确定所述超炮集中的面波特征信息;
面波频域能量确定模块,用于利用所述面波特征信息,在变换至频域的超炮集中确定面波频域能量;
面波时域能量获取模块,用于将面波频域能量变换为面波时域能量;
面波确定模块,用于利用所述面波特征信息,在所述面波时域能量中确定面波;
面波去除模块,用于去除所述超炮集中的面波。
由以上本说明书实施例提供的技术方案可见,本说明书实施例将共炮点道集进行分组得到超炮集后,通过确定超炮集中的面波特征信息,再根据所述面波特征信息,在频域和时域中针对面波能量进行确定,从而得到准确的面波分布区域。根据上述实施例,将炮集进行整合,提高了空间采样率,降低了假频影响,从而实现效果更好的面波压制。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本说明书实施例一种面波压制方法的流程图;
图2为本说明书实施例一种面波压制装置的模块图;
图3为本说明书实施例利用面波压制方法得到的超炮集的示意图;
图4为利用十字排列法压制面波的地震记录示意图;
图5为本说明书实施例利用面波压制方法压制面波的地震记录示意图;
图6为利用十字排列法压制面波的频谱示意图;
图7为本说明书实施例利用面波压制方法压制面波的频谱示意图。
具体实施方式
下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本说明书一部分实施例,而不是全部的实施例。基于本说明书中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本说明书保护的范围。
为了对地震数据中的面波进行压制,去除面波对所观测到的地震数据所造成的干扰,本说明书实施例提出了一种面波压制方法及装置。以下结合附图1说明本申请一种面波压制方法的实施例,所述方法的执行主体为服务器,所述方法的具体步骤如下:
S100:对采集到的地震数据中的共炮点道集进行分组得到超炮集。
地震数据可以包括炮点、检波点、地表高程、成像点的坐标、成像点处的速度等数据。
共炮点道集是同一个炮点激发,不同检波点接受到的地震道的集合,即同一道集中所有检波点共有的激发炮点。
超炮集即为将同一排列片内的共炮点道集进行组合所得到的道集的集合。如图1所示的排列片内,存在有五个共炮点道集,可以将这五个共炮点道集进行组合得到超炮集。通过将共炮点道集分组得到超炮集,既可以通过提高覆盖次数提高信噪比,又能缩小空间采样间隔,防止空间假频。
地震数据中还包含有炮点位置和检波点位置。一般情况下,检波点分别在不同的水平线上进行排列,从而形成检波线。根据所述炮点位置和检波点位置,能够获取炮点距和检波线距,分别表示炮点和检波线之间的距离。利用所述炮点距和检波线距,计算得到超炮集炮点数,并获取数量为超炮集炮点数的共炮点道集,作为超炮集。具体的,可以利用公式n=gl_d/sp_d计算超炮集炮点数,其中,n为超炮集炮点数,gl_d为检波线距,sp_d为炮点距。
S200:确定所述超炮集中的面波特征信息。
面波特征信息用于表示超炮集中的面波的分布范围、速度、频率等信息。在面波特征信息中,至少包括面波频率范围、面波速度范围、面波时域范围中的一种。
在确定所述超炮集中的面波特征信息之前,还可以基于超炮集建立相对坐标系。相对坐标系包括相对坐标系原点和坐标轴。相对坐标系原点可以是针对超炮集中炮点的坐标求取均值,作为相对坐标系原点。在一个具体的示例中,如图3所示,针对排列片内的五个炮点,对其坐标求取均值之后,确定第三个炮点所在的位置为相对坐标系原点。为了便于统计数据,可以设置炮线和检波线方向为坐标轴方向。炮线和检波线分别为炮点和检波点所排列形成的线。在一个具体的示例中,可以设置检波线方向为相对坐标系的横轴,炮线方向为纵轴,结合相对坐标系原点,形成相对坐标系。
在建立相对坐标系之后,还可以对炮点进行排序。针对炮点进行排序可以是针对炮点与检波线之间的距离,对炮号进行编号,从而得到各个炮点的炮号。具体的,如图3所示,在五炮八线的排列片中,依次选取出与检波线距离较近的炮点,并进行编号,最终将上述五个炮点依次分配1-5的炮号。通过对炮点分配相应的炮号。通过对炮点进行编号,能够便于在后续步骤中对炮点进行排序,从而提高空间采样率。
在得到各个炮点的炮号之后,再结合所设定的相对坐标系和炮号的先后顺序对炮点进行排序。炮号的排序决定了信号在空间中的分布,确定了面波的速度、波数域的特征及空间域的特征。本申请中利用相邻炮点的相似性进行排序,可以在保证速度、波数及空间域特征不变的前提下,提高空间采样率,从而降低假频现象。
面波时域范围表示面波在时域中的分布范围。在相对坐标系中针对超炮集截取切面,获取区域面积最大的切面,根据所述切面在相对坐标系中的位置和数据参数,通过观测确定面波的最大炮检距Off_max、最小炮检距Off_min以及时间分布范围T_max、T_min,从而确定了面波时域范围。面波速度范围表示超炮集中面波的视速度的分布范围,可以包括面波的最大速度Vmax和最小速度Vmin。在截取得到上述切面之后,结合面波时域范围,可以求取出面波速度范围。由于面波在时域分布在Off_min/Vmin、Off_min/Vmax、Off_max/Vmin、Off_min/Vmax所确定的范围内,因此结合上述面波时域范围的参数,可以得到面波的速度范围。
面波频率范围表示面波在频域的分布范围,包括面波的最高频率和最低频率,可以根据超炮集的频谱获取。具体的,可以选取面波发育区域,开时窗进行频谱分析,从而确定面波的最低频率和最高频率。
S300:根据所述面波特征信息,在变换至频域的超炮集中确定面波频域能量。
为了更好地针对面波的相关特性,确定面波分布区域,可以将超炮集变换至频域进行后续操作。在一个具体的示例中,可以利用三维傅里叶变化将超炮集变换至频域,得到频域超炮集。
在得到频域超炮集后,可以利用线性倾角扫描的方式,确定面波频率范围中斜率在面波速度范围内的区域,作为面波频域能量。具体的,可以在面波频率范围内首选选取出多个频率。在选定出的频率处,再分别选取不同的群速度,利用公式进行线性扫描,其中,f为扫描频率,f0为散射频率,v为面波速度范围内的群速度,为波数向量。散射频率和波数向量可以预先进行设置,从而确保更好的实验效果,通过不同群速度的选取及计算,得到一系列扫描频率,并选取出其中最大的扫描频率。扫描频率的大小,反映了面波在频域的能量大小。通过确定最大的扫描频率,即可确定为面波的主要频域能量分布区域,方便进行后续的筛选。
依照上述方法,在不同的频率处依次扫描出的扫描频率f最大值时的面波能量,整合所有得到的面波能量,作为面波频域能量,即可确定面波在频域的大致分布区域。
S400:将所述面波频域能量变换为面波时域能量。
若仅仅在时域或频域中对面波进行扫描确定其分布区域,并不能完全将面波与有效信号进行分离,即所得到的面波频域能量中仍然包含有有效信号的成分。为了更加准确地确定出超炮集中的面波,需要将面波频域能量变换为面波时域能量。在一个具体的示例中,可以利用三维傅里叶反变换将所述面波频域能量变换为面波时域能量。
S500:利用所述面波特征信息,在所述面波时域能量中确定面波。
在将所述面波频域能量转换为面波时域能量之后,可以结合步骤S300中在频域中扫描出面波频域能量的方法,在所述面波时域能量中确定面波。在面波时域范围内的面波时域能量中,利用线性倾角扫描确定斜率在面波速度范围内的区域,作为面波。具体的,可以是在面波时域范围内,在不同的时间点处,选取面波速度范围内的速度,利用线性倾角扫描,确定出在不同的时间点处能够选取出的最大能量。再整合所有时间点处的最大能量,即为确定的面波。
根据面波的的时域分布、频率分布以及速度等特征信息,在频域进行又一次的线性倾角扫描之后,使得面波时域能量中包含的有效信号能够被分离出去。通过将时域与频域进行串联并执行线性倾角扫描,能够确保所得到的面波包含尽可能少的有效信号,提高了压制面波的准确度。
S600:去除所述超炮集中的面波。
当确定了超炮集中的面波后,可以从所述超炮集中去除面波,从而消除面波对地震数据的干扰。
在去除所述超炮集中的面波之前,还可以将面波转化为与超炮集维度相同的面波数据体,再从超炮集中去除所述面波数据体。通过生成与超炮集维度相同的面波数据体,可以直接从超炮集中减去所述面波数据体,方便了计算,能够更加快捷地得到相应的结果。
下面结合一个具体的场景示例对本申请实施例进行进一步说明。利用本申请面波压制方法对地震数据中的面波进行压制,得到压制后的地震记录示意图和频谱图如图5和图7所示。图4和图6是利用现有技术中的十字排列法得到的地震记录示意图和频谱图。通过将图4和图5进行对比可以看出,本申请的面波压制方法相较于现有技术中去除面波的方法,能够更加有效地去除噪声的干扰,得到展示效果更好的地震记录数据。再将图6和图7得到的频谱图进行对比可以看出,本申请面波压制方法的实施例中,所得到频谱相较于现有技术所得到的频谱,结果更为清晰,突出了相应的细节。可以看出,本申请面波压制方法的实施例取得了更好的面波压制效果。
通过上述本申请面波压制方法实施例方法步骤的描述及场景示例的对比,可以看出,本申请的面波压制方法提高了采样率,降低了假频影响,能够取得更好地面波压制效果,提高了测量的准确度,更有利于后续过程中需要利用压制面波后的地震数据的步骤的开展。
以下介绍本申请一种面波压制装置的实施例,如图2所示,该面波压制装置包括:
超炮集获取模块210,用于对采集到的地震数据中的共炮点道集进行分组得到超炮集;
面波特征信息确定模块220,用于确定所述超炮集中的面波特征信息;
面波频域能量确定模块230,用于利用所述面波特征信息,在变换至频域的超炮集中确定面波频域能量;
面波时域能量获取模块240,用于将面波频域能量变换为面波时域能量;
面波确定模块250,用于利用所述面波特征信息,在所述面波时域能量中确定面波;
面波去除模块260,用于去除所述超炮集中的面波。
所述地震数据,包括炮点位置和检波点位置;所述超炮集获取模块210,包括:
数据获取子单元211,用于根据所述炮点位置和检波点位置,获取炮点距和检波线距;
超炮集炮点数计算子单元212,用于利用所述炮点距和检波线距,计算得到超炮集炮点数;
超炮集获取子单元213,用于获取数量为超炮集炮点数的共炮点道集,作为超炮集。
所述超炮集炮点数计算子单元212,包括利用公式n=gl_d/sp_d计算超炮集炮点数,其中,n为超炮集炮点数,gl_d为检波线距,sp_d为炮点距。
所述面波特征信息,包括以下至少一种:面波频率范围、面波速度范围、面波时域范围。
所述面波特征信息确定模块220,包括:
相对坐标系建立子单元221,用于基于超炮集建立相对坐标系;
速度及时域范围确定子单元222,用于根据所述相对坐标系获取面波速度范围和面波时域范围;
面波频率范围确定子单元223,用于根据超炮集的频谱确定面波频率范围。
所述坐标系建立子单元221,包括针对超炮集中炮点的坐标求取均值,作为相对坐标系原点;设置炮线方向和检波线方向为坐标轴方向。
所述面波特征信息确定模块220,还包括:
炮号确立子单元224,用于根据炮点与排列片之间的距离,对炮点进行编号,得到各个炮点的炮号;
排序子单元225,用于根据所述相对坐标系和所述炮号对炮点进行排序。
所述速度及时域范围确定子单元222,包括结合所述相对坐标系在超炮集中截取切面;根据所述切面的分布范围确定面波速度范围和面波时域范围。
所述面波频域能量确定模块230,包括:
频域超炮集获取子单元231,用于利用三维傅里叶变换将超炮集变换至频域,得到频域超炮集;
面波频域能量获取子单元232,用于利用线性倾角扫描确定面波频率范围中斜率在面波速度范围内的区域,作为面波频域能量。
所述面波频域能量获取子单元232,包括:
确定不同频率处扫描出的扫描频率f的最大值为面波能量,整合所述面波能量得到面波频域能量。
所述面波时域能量获取模块240,包括:
变换子单元241,用于利用三维傅里叶反变换将所述面波频域能量变换为面波时域能量。
所述面波确定模块250,包括:
面波获取子单元251,用于在面波时域范围内的面波时域能量中,利用线性倾角扫描确定斜率在面波速度范围内的区域,作为面波。
所述装置,还包括:
面波数据体获取模块270,用于将面波转化为与超炮集维度相同的面波数据体;
相应的,所述面波去除模块260,包括:
面波数据体去除子单元261,用于从所述超炮集中去除面波数据体。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本说明书可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本说明书的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本说明书各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本说明书可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本说明书,本领域普通技术人员知道,本说明书有许多变形和变化而不脱离本说明书的精神,希望所附的权利要求包括这些变形和变化而不脱离本说明书的精神。
Claims (11)
1.一种面波压制方法,其特征在于,所述方法包括:
对采集到的地震数据中的共炮点道集进行分组得到超炮集;
确定所述超炮集中的面波特征信息;所述面波特征信息,包括以下至少一种:面波频率范围、面波速度范围、面波时域范围;
根据所述面波特征信息,在变换至频域的超炮集中确定面波频域能量;利用所述面波特征信息,在变换至频域的超炮集中确定面波频域能量,包括:利用三维傅里叶变换将超炮集变换至频域,得到频域超炮集;利用线性倾角扫描确定面波频率范围中斜率在面波速度范围内的区域,作为面波频域能量;所述利用线性倾角扫描确定面波频率范围中斜率在面波速度范围内的区域,作为面波频域能量,包括:在面波频率范围内各个频率处选取不同的群速度,根据公式进行线性扫描,其中,f为扫描频率,f0为散射频率,v为面波速度范围内的群速度,为波数向量;确定不同频率处扫描出的扫描频率f的最大值为面波能量,整合所述面波能量得到面波频域能量;
将所述面波频域能量变换为面波时域能量;
利用所述面波特征信息,在所述面波时域能量中确定面波;
去除所述超炮集中的面波。
2.如权利要求1所述的方法,其特征在于,所述地震数据,包括炮点位置和检波点位置;所述对采集到的地震数据中的共炮点道集进行分组得到超炮集,包括:
根据所述炮点位置和检波点位置,获取炮点距和检波线距;
利用所述炮点距和检波线距,计算得到超炮集炮点数;
获取数量为超炮集炮点数的共炮点道集,作为超炮集。
3.如权利要求2所述的方法,其特征在于,所述利用所述炮点距和检波线距,计算得到超炮集炮点数,包括:
利用公式n=gl_d/sp_d计算超炮集炮点数,其中,n为超炮集炮点数,gl_d为检波线距,sp_d为炮点距。
4.如权利要求1所述的方法,其特征在于,所述确定所述超炮集中的面波特征信息,包括:
基于超炮集建立相对坐标系;
根据所述相对坐标系获取面波速度范围和面波时域范围;
根据超炮集的频谱确定面波频率范围。
5.如权利要求4所述的方法,其特征在于,所述基于超炮集建立相对坐标系,包括:
针对超炮集中炮点的坐标求取均值,作为相对坐标系原点;
设置炮线方向和检波线方向为坐标轴方向。
6.如权利要求4所述的方法,其特征在于,所述基于超炮集建立相对坐标系之后,还包括:
根据炮点与检波线之间的距离,对炮点进行编号,得到各个炮点的炮号;
根据所述相对坐标系和所述炮号对炮点进行排序。
7.如权利要求4所述的方法,其特征在于,所述根据所述相对坐标系获取面波速度范围和面波时域范围,包括:
结合所述相对坐标系在超炮集中截取切面;
根据所述切面的分布范围确定面波速度范围和面波时域范围。
8.如权利要求1所述的方法,其特征在于,所述将所述面波频域能量变换为面波时域能量,包括:
利用三维傅里叶反变换将所述面波频域能量变换为面波时域能量。
9.如权利要求1所述的方法,其特征在于,所述利用所述面波特征信息,在所述面波时域能量中确定面波,包括:
在面波时域范围内的面波时域能量中,利用线性倾角扫描确定斜率在面波速度范围内的区域,作为面波。
10.如权利要求1所述的方法,其特征在于,所述去除所述超炮集中的面波之前,还包括:
将面波转化为与超炮集维度相同的面波数据体;
相应地,从所述超炮集中去除面波,包括:
从所述超炮集中去除面波数据体。
11.一种面波压制装置,其特征在于,包括:
超炮集获取模块,用于对采集到的地震数据中的共炮点道集进行分组得到超炮集;
面波特征信息确定模块,用于确定所述超炮集中的面波特征信息;所述面波特征信息,包括以下至少一种:面波频率范围、面波速度范围、面波时域范围;
面波频域能量确定模块,用于利用所述面波特征信息,在变换至频域的超炮集中确定面波频域能量;利用所述面波特征信息,在变换至频域的超炮集中确定面波频域能量,包括:利用三维傅里叶变换将超炮集变换至频域,得到频域超炮集;利用线性倾角扫描确定面波频率范围中斜率在面波速度范围内的区域,作为面波频域能量;所述利用线性倾角扫描确定面波频率范围中斜率在面波速度范围内的区域,作为面波频域能量,包括:在面波频率范围内各个频率处选取不同的群速度,根据公式进行线性扫描,其中,f为扫描频率,f0为散射频率,v为面波速度范围内的群速度,为波数向量;确定不同频率处扫描出的扫描频率f的最大值为面波能量,整合所述面波能量得到面波频域能量;
面波时域能量获取模块,用于将面波频域能量变换为面波时域能量;
面波确定模块,用于利用所述面波特征信息,在所述面波时域能量中确定面波;
面波去除模块,用于去除所述超炮集中的面波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910255825.0A CN111766631B (zh) | 2019-04-01 | 2019-04-01 | 一种面波压制方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910255825.0A CN111766631B (zh) | 2019-04-01 | 2019-04-01 | 一种面波压制方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111766631A CN111766631A (zh) | 2020-10-13 |
CN111766631B true CN111766631B (zh) | 2022-12-02 |
Family
ID=72718832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910255825.0A Active CN111766631B (zh) | 2019-04-01 | 2019-04-01 | 一种面波压制方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111766631B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112255683A (zh) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | 地震资料的噪声压制方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007044858A1 (en) * | 2005-10-11 | 2007-04-19 | Steady State Imaging | Frequency swept excitation for magnetic resonance |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2326971A4 (en) * | 2008-08-11 | 2017-06-14 | Exxonmobil Upstream Research Company | Removal of surface-wave noise in seismic data |
CN101930081B (zh) * | 2009-06-26 | 2012-05-30 | 中国石油集团东方地球物理勘探有限责任公司 | 地震记录单频干扰的自动识别与压制方法 |
CN102176056B (zh) * | 2011-02-18 | 2012-11-28 | 中国石油化工股份有限公司 | 基于时频分析的地震有效波多域能量补偿方法 |
CN102305945B (zh) * | 2011-06-20 | 2013-09-25 | 电子科技大学 | 一种线性噪声消除方法 |
US9442209B2 (en) * | 2012-07-10 | 2016-09-13 | Pgs Geophysical As | Methods and systems for reconstruction of low frequency particle velocity wavefields and deghosting of seismic streamer data |
US10126453B2 (en) * | 2015-01-21 | 2018-11-13 | Baker Hughes, A Ge Company, Llc | Processing seismic data to remote noise |
CN106707342B (zh) * | 2016-12-12 | 2018-09-04 | 中国石油天然气集团公司 | 共炮点道集多级面波压制方法和装置 |
-
2019
- 2019-04-01 CN CN201910255825.0A patent/CN111766631B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007044858A1 (en) * | 2005-10-11 | 2007-04-19 | Steady State Imaging | Frequency swept excitation for magnetic resonance |
Also Published As
Publication number | Publication date |
---|---|
CN111766631A (zh) | 2020-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110658557B (zh) | 基于生成对抗网络的地震数据面波压制方法 | |
CN108983287B (zh) | 一种基于凸集投影算法的曲波变换反假频地震数据重建方法 | |
CN107193045B (zh) | 一种地震数据处理方法及装置 | |
CN102681014A (zh) | 基于多项式拟合的规则线性干扰压制方法 | |
CN106707342B (zh) | 共炮点道集多级面波压制方法和装置 | |
CN105242309A (zh) | 一种地震数据中规则干扰的压制方法及装置 | |
CN111766631B (zh) | 一种面波压制方法及装置 | |
CN111708087A (zh) | 一种基于DnCNN神经网络对地震数据噪声压制的方法 | |
CN110907989A (zh) | 重建拟地面地震反射波成像方法及系统 | |
CN107561588B (zh) | 一种地震数据噪声压制方法及装置 | |
CN109581494B (zh) | 叠前偏移方法及装置 | |
CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
CN111929732B (zh) | 一种地震数据去噪方法、装置及设备 | |
CN109917459A (zh) | 一种压制地震噪声的方法、装置及系统 | |
CN112630837A (zh) | 一种地震品质因子估算方法、装置、设备及存储介质 | |
Ma et al. | A Global and Multi-Scale Denoising Method Based on Generative Adversarial Network for DAS VSP Data | |
Liu et al. | Research on a matching detection method for magnetic anomaly of underwater targets | |
CN113391351B (zh) | 一种基于被动源地震波场分析提取矿集区结构的方法 | |
CN113721297B (zh) | 利用速度特征压制沙丘鸣震的方法 | |
CN110231651A (zh) | 一种基于振动信号矢量分解原理的叠前地震数据去噪方法 | |
CN112444878B (zh) | 地震数据分频初至层析静校正方法、存储介质和计算设备 | |
CN113640882B (zh) | 采集脚印噪声去除方法、电子设备及计算机可读存储介质 | |
CN113126163B (zh) | 五维地震数据噪声衰减方法及装置 | |
CN117434569B (zh) | 基于小波去噪的固定基线长度动态解算方法 | |
CN112379427A (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 |