CN105093329B - 地震数据的低频信息的自动补偿方法及装置 - Google Patents
地震数据的低频信息的自动补偿方法及装置 Download PDFInfo
- Publication number
- CN105093329B CN105093329B CN201510524870.3A CN201510524870A CN105093329B CN 105093329 B CN105093329 B CN 105093329B CN 201510524870 A CN201510524870 A CN 201510524870A CN 105093329 B CN105093329 B CN 105093329B
- Authority
- CN
- China
- Prior art keywords
- frequency
- low
- geological data
- initial
- cut
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例提供了一种地震数据的低频信息的自动补偿方法及装置,该方法包括:利用地震数据估算地震子波的归一化振幅谱;在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,所述起始频率小于所述截止频率,0<A<1;将幂运算得到的结果除以所述归一化振幅谱,得到所述地震数据在所述起始频率至所述截止频率区间的补偿算子的频率响应曲线;根据补偿算子的频率响应曲线对地震数据进行低频补偿,本申请实施例实现了不通过采用特殊采集设备也能对地震数据进行低频补偿,而且本申请的补偿算子是从地震数据中直接提取得到的,更加客观。其相对于现有技术而言,本申请是完全数据驱动的、基于子波估计的自动补偿方法。
Description
技术领域
本申请涉及低频地震勘探技术领域,尤其是涉及一种地震数据的低频信息的自动补偿方法及装置。
背景技术
低频地震勘探是近年来地球物理勘探技术领域的热点,这是由于:1)低频信息在油气藏检测与指示中具有重要意义,例如油气藏低频伴影现象;2)低频信号影响地震勘探的垂向分辨率,地震剖面的分辨率是由地震子波的带宽决定,同时需要高频和低频信息;3)低频信号在深层油藏探测(如火成岩地区勘探)中具有重要作用,随着勘探目的层深度的加大,高频信号衰减、散射、频散严重,而低频信号保留相对完整;4)低频信息在地震数据处理的反演(如全波形反演FWI、以及常规的波阻抗反演等)过程中具有决定性作用。
低频地震勘探的关键在于两个方面:一是野外地震勘探采集时的低频震源激发和对应的检波器接收;另一方面是室内能够进行保护低频信号的处理,如低频面波噪声压制过程中保留有效的低频信息。在实际的操作过程中,从野外采集到室内处理,地震数据中的低频信号都可能、或者说不可避免的受到损伤并被压制。例如,陆上低频可控震源的设计和制造本身就是难题,这是因为检波器组合、频响也会对接收到低频信号带来损失。此外,既要保持常规震源的高频能量,又要提高震源低频信号的下传能力,在大多数采集环境下,噪声强度将随频率降低呈指数增加。在室内处理中,低频的面波是主要的干扰源,它与低频信号在频带上存在重叠,压制干扰的过程中低频信号极易受到损伤。
目前,由于硬件成本等各方面条件的限制,采用特殊设备的低频地震采集并未普及。因而,对于采用常规震源和检波器采集的地震数据,如何低成本地利用现有信息增强、补偿数据中的低频信号以满足低频地震勘探的需求,具有重要的现实意义。
发明内容
本申请实施例的目的在于提供一种地震数据的低频信息的自动补偿方法及装置,以低成本实现地震数据的低频信息的自动补偿。
为达到上述目的,一方面,本申请实施例提供了一种地震数据的低频信息的自动补偿方法,包括以下步骤:
利用地震数据估算地震子波的归一化振幅谱;
在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,所述起始频率小于所述截止频率,0<A<1;
将幂运算得到的结果除以所述归一化振幅谱,得到所述地震数据在所述起始频率至所述截止频率区间的补偿算子的频率响应曲线;
根据所述补偿算子的频率响应曲线对所述地震数据进行低频补偿。
本申请实施例的地震数据的低频信息的自动补偿方法,在所述根据补偿算子的频率响应曲线对所述地震数据进行低频补偿之前,还包括:
对所述补偿算子的频率响应曲线进行平滑处理。
本申请实施例的地震数据的低频信息的自动补偿方法,所述对补偿算子的频率响应曲线进行平滑处理,具体包括:
将所述起始频率减去预设的平滑斜坡长度得到第一频率阈值F1;
判断F1是否小于零,如果小于零,取F1为零;
从所述起始频率至F1区间,将各频率对应的幅值定为由所述起始频率对应的幅值渐变至1;
将所述截止频率加上所述平滑斜坡长度得到第二频率阈值F2;
从所述截止频率至F2区间,将各频率对应的幅值定为由所述截止频率对应的幅值渐变至1;
将F1至F2区间以外的频率幅值设定为1。
本申请实施例的地震数据的低频信息的自动补偿方法,所述预设的截止频率通过以下方式得到:
定义初始截止频率;
从所述起始频率起,搜索所述归一化振幅谱的第一个不小于预设阈值的峰值,并将所述第一个不小于所述预设阈值的峰值所对应的频率定义为第一频率,其中所述预设阈值小于1;
判断所述初始截止频率是否大于所述第一频率;
如果大于,则将所述第一频率确定为预设的截止频率。
本申请实施例的地震数据的低频信息的自动补偿方法,还包括:
如果所述初始截止频率不大于所述第一频率,则将所述初始截止频率确定为预设的截止频率。
本申请实施例的地震数据的低频信息的自动补偿方法,所述预设阈值的取值范围为0.6~1。
本申请实施例的地震数据的低频信息的自动补偿方法,所述截止频率取值范围为0~30Hz。
另一方面,本申请实施例还提供了一种地震数据的低频信息的自动补偿装置,包括:
振幅谱获取模块,用于利用地震数据估算地震子波的归一化振幅谱;
幂运算模块,用于在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,所述起始频率小于所述截止频率,0<A<1;
补偿算子获取模块,用于将幂运算得到的结果除以所述归一化振幅谱,得到所述地震数据在所述起始频率至所述截止频率区间的补偿算子的频率响应曲线;
低频补偿模块,用于根据所述补偿算子的频率响应曲线对所述地震数据进行低频补偿。
本申请实施例的地震数据的低频信息的自动补偿装置,其特征在于,还包括:
平滑处理模块,用于在所述根据补偿算子的频率响应曲线对所述地震数据进行低频补偿之前,对所述补偿算子的频率响应曲线进行平滑处理。
本申请实施例的地震数据的低频信息的自动补偿装置,所述平滑处理模块具体包括:
第一频率阈值获取子模块,用于将所述起始频率减去预设的平滑斜坡长度得到第一频率阈值F1;
第一判断子模块,用于判断F1是否小于零,如果小于零,取F1为零;
第一渐变子模块,用于并从所述起始频率至F1区间,将各频率对应的幅值定为由所述起始频率对应的幅值渐变至1;
第二频率阈值获取子模块,用于将所述截止频率加上所述平滑斜坡长度得到第二频率阈值F2;
第二渐变子模块,用于从所述截止频率至F2区间,将各频率对应的幅值定为由所述截止频率对应的幅值渐变至1;
第一设定子模块,用于将F1至F2区间以外的频率幅值设定为1。
本申请实施例的地震数据的低频信息的自动补偿装置,所述预设的截止频率通过以下方式得到:
定义初始截止频率;
从所述起始频率起,搜索所述归一化振幅谱的第一个不小于预设阈值的峰值,并将所述第一个不小于所述预设阈值的峰值所对应的频率定义为第一频率,其中所述预设阈值小于1;
判断所述初始截止频率是否大于所述第一频率;
如果大于,则将所述第一频率确定为预设的截止频率。
本申请实施例的地震数据的低频信息的自动补偿装置,还包括:
如果所述初始截止频率不大于所述第一频率,则将所述初始截止频率确定为预设的截止频率。
本申请实施例的地震数据的低频信息的自动补偿装置,所述预设阈值的取值范围为0.6~1。
本申请实施例的地震数据的低频信息的自动补偿装置,所述截止频率取值范围为0~30Hz。
本申请实施例利用地震数据估算地震子波的归一化振幅谱;其次在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,0<A<1;然后将得到的结果除以归一化振幅谱,得到地震数据在起始频率至截止频率区间的补偿算子的频率响应曲线;最后根据补偿算子的频率响应曲线对地震数据进行低频补偿,从而实现了不通过采用特殊设备也能对地震数据进行低频补偿,而且本申请的补偿算子是从地震数据中直接提取得到的,更加客观,其相对于现有技术而言,能够更好地防止过补偿而造成的低频地震数据损伤。
附图说明
此处所说明的附图用来提供对本申请实施例的进一步理解,构成本申请实施例的一部分,并不构成对本申请实施例的限定。在附图中:
图1为本申请实施例中地震数据的低频信息的自动补偿方法的流程图;
图2为本申请实施例中地震数据振幅谱、地震子波振幅谱及开放运算结果示意图;
图3为本申请实施例中未经平滑处理及已经平滑处理的补偿算子的频率响应曲线的振幅谱示意图;
图4为本申请实施例中地震数据及低频补偿后地震数据的振幅谱示意图;
图5为本申请实施例中未经低频补偿处理的地震数据的示意图;
图6为本申请实施例中低频补偿处理后的地震数据的示意图;
图7为本申请实施例中地震数据的低频信息的自动补偿装置的逻辑结构框图;
图8为本申请实施例中平滑处理模块的逻辑结构框图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本申请实施例做进一步详细说明。在此,本申请实施例的示意性实施例及其说明用于解释本申请实施例,但并不作为对本申请实施例的限定。
下面结合附图,对本申请实施例的具体实施方式作进一步的详细说明。
参考图1所示,本申请实施例的地震数据的低频信息的自动补偿方法包括以下步骤:
S1、利用地震数据估算地震子波的归一化振幅谱。本申请实施例中,获取地震子波振幅谱的方法可以为通用的子波振幅谱估算方法:例如用谱模拟方法估算地震子波振幅谱,或用地震数据自相关来估算地震子波振幅谱。估算地震子波振幅谱,可利用单道地震数据,也可利用多道地震数据,并且可以选择时窗(如图2所示,地震数据的振幅谱为多道数据的叠加振幅谱,且是局部时窗数据,用以提高数据的信噪比)。获得地震子波的振幅谱要进行归一化处理的目的是为了简化运算。
S2、在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,起始频率小于截止频率,0<A<1。而A的大小决定了低频补偿的强弱,一般的,A越小,补偿越强,对应的,要耗费的运算越多。
S3、将幂运算得到的结果除以所述归一化振幅谱,得到所述地震数据在所述起始频率至所述截止频率区间的补偿算子的频率响应曲线。
S4、根据所述补偿算子的频率响应曲线对所述地震数据进行低频补偿。
本申请实施例中,所述截止频率取值范围为0~30Hz;而所述预设的截止频率通过以下方式得到:
首先,定义初始截止频率;
其次,从所述起始频率起,搜索所述归一化振幅谱的第一个不小于预设阈值的峰值,并将所述第一个不小于所述预设阈值的峰值所对应的频率定义为第一频率,其中所述预设阈值的取值范围为0.6~1;
然后,判断所述初始截止频率是否大于所述第一频率;
最后,如果大于,则将所述第一频率确定为预设的截止频率,以防止低频补过;而如果所述初始截止频率不大于所述第一频率,则将所述初始截止频率确定为预设的截止频率。
本申请另一实施例中,在所述根据补偿算子的频率响应曲线对所述地震数据进行低频补偿之前,还包括对所述补偿算子的频率响应曲线进行平滑处理的步骤。该平滑处理具体包括如下步骤:
(1)、将所述起始频率减去预设的平滑斜坡长度得到第一频率阈值F1。一般的,平滑斜坡长度的取值范围为1~10Hz。
(2)、判断F1是否小于零。
(3)、如果小于,则从所述起始频率至F1区间,将各频率对应的幅值定为由所述起始频率对应的幅值渐变至1将各频率对应的幅值定为由所述起始频率对应的幅值渐变至1。当然,如果不小于零,则将0频率对应的幅值设定为1。
(4)、将所述截止频率加上所述平滑斜坡长度得到第二频率阈值F2。
(5)、从所述截止频率至F2区间,将各频率对应的幅值定为由所述截止频率对应的幅值渐变至1。
(6)、将F1至F2区间以外的频率幅值设定为1。
上述平滑处理可以降低或减弱吉布斯效应。此外,上述这种平滑处理可使用简单的线性变化;当然也可使用较复杂算法以获取更好的平滑效果,例如使补偿算子的频率响应曲线在所述第一频率阈值、所述起始频率、所述截止频率、所述第二频率阈值等端点所对应的幅值一阶可导。
本申请实施例中,图2中显示的地震数据的振幅谱为多道数据的叠加振幅谱,且是局部时窗数据,用以提高数据的信噪比;地震子波振幅谱是利用谱模拟方法直接从地震数据振幅谱自动估算而得。为简化运算,图中的振幅谱都进行了归一化处理。
本申请实施例中,图3中显示的是利用子波振幅谱开2次方(即A=0.5)再除以子波振幅谱得到的低频补偿算子的振幅谱,实线为最终使用的低频补偿算子的振幅谱。注意:在定义的补偿范围BF~EF之间(其中,BF为起始频率,EF为终止频率),实线与虚线重合;在BF之前及EF之后,算子振幅谱渐变至1;在BF及EF位置附近,算子振幅谱进行了平滑处理。
本申请实施例中,图4中显示的是低频补偿前后地震数据的振幅谱。可以看出,在定义的低频补偿范围中,地震数据的振幅谱能量得到明显提升,在范围之外的中高频,振幅谱能量则保持不变。注意:这里的地震数据振幅谱是多道数据的叠加振幅谱,且使用的是整道数据,未分时窗,因此原始输入数据的振幅谱与图1(估算子波振幅谱使用的是局部时窗数据)略有差异。
本申请实施例中,图5、图6分别为输入的原始地震数据以及低频补偿后的地震数据。仔细对比可以发现,地震数据中的低频信息得到增强,特别是中深层;数据的信噪比没有明显降低;中深层的、低频的同相轴更为清晰、连续。
本申请实施例中,图4~图6说明了本申请实施例具有良好的、稳健的结果,与本申请实施例的预期目标一致。
本申请实施例利用地震数据估算地震子波的归一化振幅谱;其次在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,0<A<1;然后将得到的结果除以归一化振幅谱,得到地震数据在起始频率至截止频率区间的补偿算子的频率响应曲线;最后根据补偿算子的频率响应曲线对地震数据进行低频补偿,从而实现了不通过采用特殊设备也能对地震数据进行低频补偿,而且本申请的补偿算子是从地震数据中提取得到的,其相对于现有技术而言,能够更好地防止过补偿而造成的低频地震数据损伤。
参考图7所示,与上述低频地震勘探中低频地震数据的自动补偿对应,本申请实施例的地震数据的低频信息的自动补偿装置包括:
振幅谱获取模块71,用于利用地震数据估算地震子波的归一化振幅谱。本申请实施例中,获取地震子波振幅谱的方式可以为通用的子波振幅谱估算方式:例如用谱模拟方法估算地震子波振幅谱,或用地震数据自相关来估算地震子波振幅谱。估算地震子波振幅谱,可利用单道地震数据,也可利用多道地震数据,并且可以选择时窗。获得地震子波的振幅谱要进行归一化处理的目的是为了简化运算。
幂运算模块72,用于在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,所述起始频率小于所述截止频率,0<A<1。而A的大小决定了低频补偿的强弱,一般的,A越小,补偿越强,对应的,要耗费的运算越多。
补偿算子获取模块73,用于将幂运算得到的结果除以所述归一化振幅谱,得到所述地震数据在所述起始频率至所述截止频率区间的补偿算子频率响应曲。
低频补偿模块74,用于根据所述补偿算子的频率响应曲线对所述地震数据进行低频补偿。本申请实施例中,所述截止频率取值范围为0~30Hz;而所述预设的截止频率通过以下方式得到:
首先,定义初始截止频率;
其次,从所述起始频率起,搜索所述归一化振幅谱的第一个不小于预设阈值的峰值,并将所述第一个不小于所述预设阈值的峰值所对应的频率定义为第一频率,其中所述预设阈值的取值范围为0.6~1;
然后,判断所述初始截止频率是否大于所述第一频率;
最后,如果大于,则将所述第一频率确定为预设的截止频率,以防止低频补过;而如果所述初始截止频率不大于所述第一频率,则将所述初始截止频率确定为预设的截止频率。
在本申请的另一实施例中,地震数据的低频信息的自动补偿装置还可以包括:
平滑处理模块75,用于在所述根据补偿算子的频率响应曲线对所述地震数据进行低频补偿之前,对所述补偿算子的频率响应曲线进行平滑处理。
结合图8所示,上述平滑处理模块75具体包括:
第一频率阈值获取子模块751,用于将所述起始频率减去预设的平滑斜坡长度得到第一频率阈值F1;
第一判断子模块752,用于判断F1是否小于零,如果小于零,取F1为零;
第一渐变子模块753,用于并从所述起始频率至F1区间,将各频率对应的幅值定为由所述起始频率对应的幅值渐变至1;
第二频率阈值获取子模块754,用于将所述截止频率加上所述平滑斜坡长度得到第二频率阈值F2;
第二渐变子模块755,用于从所述截止频率至F2区间,将各频率对应的幅值定为由所述截止频率对应的幅值渐变至1;
第一设定子模块756,用于将F1至F2区间以外的频率幅值设定为1。
上述平滑处理可以降低或减弱吉布斯效应。此外,上述这种平滑处理可使用简单的线性变化;当然也可使用较复杂算法以获取更好的平滑效果,例如使补偿算子的频率响应曲线在所述第一频率阈值、所述起始频率、所述截止频率、所述第二频率阈值等端点所对应的幅值一阶可导。
本申请实施例利用地震数据估算地震子波的归一化振幅谱;其次在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,0<A<1;然后将得到的结果除以归一化振幅谱,得到地震数据在起始频率至截止频率区间的补偿算子的频率响应曲线;最后根据补偿算子的频率响应曲线对地震数据进行低频补偿,从而实现了不通过采用特殊设备也能对地震数据进行低频补偿,而且本申请的补偿算子是从地震数据中直接提取得到的,更加客观,其相对于现有技术而言,能够更好地防止过补偿而造成的低频地震数据损伤。
本领域技术人员还可以了解到本申请实施例列出的各种说明性逻辑块、单元和步骤可以通过硬件、软件或两者的结合来实现。至于是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本申请实施例保护的范围。
本申请实施例中所描述的各种说明性的逻辑块,或单元都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本申请实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本申请实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
以上所述的具体实施例,对本申请的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本申请实施例的具体实施例而已,并不用于限定本申请的保护范围,凡在本申请的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (14)
1.一种地震数据的低频信息的自动补偿方法,其特征在于,包括以下步骤:
利用地震数据估算地震子波的归一化振幅谱;
在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,所述起始频率小于所述截止频率,0<A<1;
将幂运算得到的结果除以所述归一化振幅谱,得到所述地震数据在所述起始频率至所述截止频率区间的补偿算子的频率响应曲线;
根据所述补偿算子的频率响应曲线对所述地震数据进行低频补偿。
2.根据权利要求1所述的地震数据的低频信息的自动补偿方法,其特征在于,在所述根据补偿算子的频率响应曲线对所述地震数据进行低频补偿之前,还包括:
对所述补偿算子的频率响应曲线进行平滑处理。
3.根据权利要求2所述的地震数据的低频信息的自动补偿方法,其特征在于,所述对补偿算子的频率响应曲线进行平滑处理,具体包括:
将所述起始频率减去预设的平滑斜坡长度得到第一频率阈值F1;
判断F1是否小于零,如果小于零,取F1为零;
从所述起始频率至F1区间,将各频率对应的幅值定为由所述起始频率对应的幅值渐变至1;
将所述截止频率加上所述平滑斜坡长度得到第二频率阈值F2;
从所述截止频率至F2区间,将各频率对应的幅值定为由所述截止频率对应的幅值渐变至1;
将F1至F2区间以外的频率幅值设定为1。
4.根据权利要求1所述的地震数据的低频信息的自动补偿方法,其特征在于,所述预设的截止频率通过以下方式得到:
定义初始截止频率;
从所述起始频率起,搜索所述归一化振幅谱的第一个不小于预设阈值的峰值,并将所述第一个不小于所述预设阈值的峰值所对应的频率定义为第一频率,其中所述预设阈值小于1;
判断所述初始截止频率是否大于所述第一频率;
如果大于,则将所述第一频率确定为预设的截止频率。
5.根据权利要求4所述的地震数据的低频信息的自动补偿方法,其特征在于,还包括:
如果所述初始截止频率不大于所述第一频率,则将所述初始截止频率确定为预设的截止频率。
6.根据权利要求4或5所述的地震数据的低频信息的自动补偿方法,其特征在于,所述预设阈值的取值范围为0.6~1。
7.根据权利要求1所述的地震数据的低频信息的自动补偿方法,其特征在于,所述截止频率取值范围为0~30Hz。
8.一种地震数据的低频信息的自动补偿装置,其特征在于,包括:
振幅谱获取模块,用于利用地震数据估算地震子波的归一化振幅谱;
幂运算模块,用于在预设的起始频率至预设的截止频率区间,计算所述归一化振幅谱的幅值的A次幂,其中,0<A<1;
补偿算子获取模块,用于将幂运算得到的结果除以所述归一化振幅谱,得到所述地震数据在所述起始频率至所述截止频率区间的补偿算子的频率响应曲线;
低频补偿模块,用于根据所述补偿算子的频率响应曲线对所述地震数据进行低频补偿。
9.根据权利要求8所述的地震数据的低频信息的自动补偿装置,其特征在于,还包括:
平滑处理模块,用于在所述根据补偿算子的频率响应曲线对所述地震数据进行低频补偿之前,对所述补偿算子的频率响应曲线进行平滑处理。
10.根据权利要求9所述的地震数据的低频信息的自动补偿装置,其特征在于,所述平滑处理模块,具体包括:
第一频率阈值获取子模块,用于将所述起始频率减去预设的平滑斜坡长度得到第一频率阈值F1;
第一判断子模块,用于判断F1是否小于零,如果小于零,取F1为零;
第一渐变子模块,用于并从所述起始频率至F1区间,将各频率对应的幅值定为由所述起始频率对应的幅值渐变至1;
第二频率阈值获取子模块,用于将所述截止频率加上所述平滑斜坡长度得到第二频率阈值F2;
第二渐变子模块,用于从所述截止频率至F2区间,将各频率对应的幅值定为由所述截止频率对应的幅值渐变至1;
第一设定子模块,用于将F1至F2区间以外的频率幅值设定为1。
11.根据权利要求8所述的地震数据的低频信息的自动补偿装置,其特征在于,所述预设的截止频率通过以下方式得到:
定义初始截止频率;
从所述起始频率起,搜索所述归一化振幅谱的第一个不小于预设阈值的峰值,并将所述第一个不小于所述预设阈值的峰值所对应的频率定义为第一频率,其中所述预设阈值小于1;
判断所述初始截止频率是否大于所述第一频率;
如果大于,则将所述第一频率确定为预设的截止频率。
12.根据权利要求11所述的地震数据的低频信息的自动补偿装置,其特征在于,还包括:
如果所述初始截止频率不大于所述第一频率,则将所述初始截止频率确定为预设的截止频率。
13.根据权利要求11或12所述的地震数据的低频信息的自动补偿装置,其特征在于,所述预设阈值的取值范围为0.6~1。
14.根据权利要求8所述的地震数据的低频信息的自动补偿装置,其特征在于,所述截止频率取值范围为0~30Hz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510524870.3A CN105093329B (zh) | 2015-08-25 | 2015-08-25 | 地震数据的低频信息的自动补偿方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510524870.3A CN105093329B (zh) | 2015-08-25 | 2015-08-25 | 地震数据的低频信息的自动补偿方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093329A CN105093329A (zh) | 2015-11-25 |
CN105093329B true CN105093329B (zh) | 2017-07-07 |
Family
ID=54574183
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510524870.3A Active CN105093329B (zh) | 2015-08-25 | 2015-08-25 | 地震数据的低频信息的自动补偿方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105093329B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866838B (zh) * | 2016-05-17 | 2018-10-16 | 中国石油天然气股份有限公司 | 一种地震数据低频信息补偿方法及装置 |
CN106772574B (zh) * | 2016-11-28 | 2019-02-15 | 中国石油天然气集团有限公司 | 一种保持低频信息的谱延拓方法 |
CN110050205B (zh) * | 2016-12-02 | 2022-03-11 | Bp北美公司 | 使用偏移道集的潜波照明 |
CN112034510B (zh) * | 2019-06-04 | 2023-04-07 | 中国石油天然气股份有限公司 | 低频补偿方法及低频补偿装置 |
CN111427088A (zh) * | 2020-03-11 | 2020-07-17 | 王仰华 | 用于识别薄互储层的地震数据低频补偿方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4875166A (en) * | 1987-10-09 | 1989-10-17 | Input/Output, Inc. | Bandwidth enhancing seismic acquisition system and method |
MY129095A (en) * | 2001-02-13 | 2007-03-30 | Exxonmobil Upstream Res Co | Method for spectral balancing of near-and far-offset seismic data. |
CN103323875B (zh) * | 2012-03-20 | 2015-08-19 | 中国石油化工股份有限公司 | 一种振幅归一化处理方法 |
CN103645502B (zh) * | 2013-12-11 | 2016-08-17 | 中国海洋石油总公司 | 一种曲波域中地震波衰减补偿方法 |
-
2015
- 2015-08-25 CN CN201510524870.3A patent/CN105093329B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105093329A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093329B (zh) | 地震数据的低频信息的自动补偿方法及装置 | |
US11598203B2 (en) | Reducing resonant noise in seismic data acquired using a distributed acoustic sensing system | |
US11435493B2 (en) | Enhanced waveform analysis for target modes of borehole waves | |
CN104395779B (zh) | 用于估计和减弱地震数据中的噪声的系统和方法 | |
CN106199702B (zh) | 地震数据存储磁带掉磁粉产生的异常振幅压制方法 | |
US11194070B2 (en) | Wavelet transform-based coherent noise reduction in distributed acoustic sensing | |
CN106646616B (zh) | 一种补偿地层吸收衰减的方法及装置 | |
CN104808245A (zh) | 道集优化处理方法及其装置 | |
Ozgun Konca et al. | Kinematic inversion of physically plausible earthquake source models obtained from dynamic rupture simulations | |
CN111060961B (zh) | 基于多信息约束反演的品质因子确定方法、装置及系统 | |
CN105044772A (zh) | 一种实现时变谱模拟反褶积的方法和装置 | |
CN103675904A (zh) | 一种井震匹配解释性目标处理方法及装置 | |
CN110989020A (zh) | 一种音频大地电磁数据噪声干扰的滤波方法及系统 | |
CN105301656B (zh) | 一种压制虚反射信号的方法及装置 | |
CN103954994B (zh) | 基于连续小波变换的地震信号增强方法及装置 | |
US20190025451A1 (en) | Identifying and visually presenting formation slowness based on low-frequency dispersion asymptotes | |
WO2021072072A1 (en) | Look-ahead vsp workflow that uses a time and depth variant q to reduce uncertainties in depth estimation ahead of a drilling bit | |
CN105425297B (zh) | 一种压制虚反射信号的方法及装置 | |
US20190219718A1 (en) | Dipole Shear Velocity Estimation | |
CN104063619B (zh) | 一种累计声发射信号kaiser点的识别方法 | |
US11366241B2 (en) | Spatially adaptive vibrator sweep parameter selection during seismic data acquisition | |
US12000972B2 (en) | Attenuation of interface waves using single component seismic data | |
CN110873900A (zh) | 频率域叠前地震道q补偿方法及系统 | |
KR102505997B1 (ko) | 지하수 수위 변동에 기반한 지진 전조 분석 방법 및 이를 수행하기 위한 컴퓨팅 장치 | |
CN112305616B (zh) | 光纤井中地震数据剖面的获取方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |