一种地震记录规则干扰波的提取方法及装置的制造方法
【专利摘要】本发明涉及地震资料处理领域，尤其涉及一种地震记录规则干扰波的提取方法及装置。该方法包括：获取预设滑动窗内的叠前地震记录；对所述叠前地震记录进行干扰波时移排齐处理，获取时移排齐处理后的地震记录；计算所述时移排齐处理后的地震记录中地震道的每个采样点对应的相似权系数和能量权系数；根据所述地震道的每个采样点对应的相似权系数和能量权系数，对所述时移排齐处理后的地震记录中地震道进行加权叠加处理，获得所述时移排齐处理后的地震记录中心道对应的叠加地震道；根据所述叠前地震记录和所述叠加地震道，获取剩余地震记录。本申请实施例的方法能有效的实现对干扰波的提取，为地震资料的解释奠定了基础。
【专利说明】
一种地震记录规则干扰波的提取方法及装置
技术领域
[0001] 本发明涉及地震资料处理领域，尤其涉及一种地震记录规则干扰波的提取方法及 装置。
【背景技术】
[0002] 在采集野外实际地震数据包含着丰富的信息，除了能够反映地下构造的有效反射 波，通常地震记录上还包含着各种各样的干扰波。干扰波可分为两种，具有明显传播规律的 规则干扰波，如声波、面波，多次波等，和没有明显传播规律性的随机干扰波。干扰波的存在 会降低有效波的质量，影响我们对地下地质目标的认识，所以，如何在对地震记录上的有效 反射波没有明显伤害的同时，将地震记录中的干扰波去除，从而提高有效信号的分辨率，始 终是地震资料处理中的一个关键问题。
[0003] 目前，根据不同的地震波属性，有很多种规则干扰波提取或压制方法，如中值滤波 法、F-K滤波法、ip变换法、矢量合成法、极化分析法、带通滤波法等。其中，中值滤波法、F-K 滤波法、T-p变换法主要是利用了规则干扰波和有效波的视速度差异;矢量合成法和极化分 析法利用了规则干扰波和有效波的偏振特性的不同；带通滤波法则是根据规则干扰波和有 效波分布的频率不同进行分离。
[0004] 上述规则干扰波提取方法都是基于规则干扰波和有效波在视速度、偏振特性和频 率等属性方面的明显差异，当规则干扰波和有效波的上述属性差异很小或者没有差异，上 述方法就不能对干扰波进行有效提取或压制。因此，如何有效的提取地震记录中的规则干 扰波仍是目前地震资料处理中亟需解决的问题。
【发明内容】
[0005] 本申请实施例提供了一种地震记录规则干扰波的提取方法及装置，以有效的提取 地震记录中的规则干扰波。
[0006] 为达到上述目的，一方面，本申请实施例提供了一种地震记录规则干扰波的提取 方法，所述方法包括：
[0007] 获取预设滑动窗内的叠前地震记录；
[0008] 对所述叠前地震记录进行干扰波时移排齐处理，获取时移排齐处理后的地震记 录；
[0009] 计算所述时移排齐处理后的地震记录中地震道每个采样点对应的相似权系数和 能量权系数；
[0010] 根据所述地震道的每个采样点对应的相似权系数和能量权系数，对所述时移排齐 处理后的地震记录中地震道进行加权叠加处理，获得所述时移排齐处理后的地震记录中心 道对应的叠加地震道;其中，所述叠加地震道为提取的规则干扰波；
[0011] 根据所述叠前地震记录和所述叠加地震道，获取剩余地震记录;其中，所述剩余地 震记录为去除规则干扰波的地震记录。
[0012] 进一步地，在根据所述叠前地震记录和所述叠加地震道，获取剩余地震记录之后， 所述方法还包括：
[0013] 逐步减小所述滑动窗的大小，对所述剩余地震记录进行规则干扰波去除处理。
[0014] 进一步地，采用以下计算公式获得所述时移排齐处理后地震记录中心道对应的叠 加地震道：
[0015] ^(j) = ^ii,(t)W(t)PV(t) k=l.
[0016] 其中，s(t)表示叠加地震道；表示地震道每个采样点对应的相似权系数； 表示地震道每个采样点对应的能量权系数;uk(t)表示时移排齐处理后的地震记录上 第k道地震道在时间样点序号为t时的振幅值;k表示地震道的序号;Mi表示滑动窗内地震记 录的道数;t表示时间样点序号，N为每道地震记录的时间采样点数。
[0017] 进一步地，根据所述叠前地震记录和所述叠加地震道，获取剩余地震记录，包括：
[0018] 对所述叠加地震道进行反时移校正处理，获得反时移校正处理后的叠加地震道；
[0019] 根据所述叠前地震记录和所述反时移校正处理后的叠加地震道，获取振幅匹配系 数；
[0020] 根据所述振幅匹配系数对所示叠加地震道进行振幅匹配处理，利用所述叠前地震 记录减去振幅匹配处理后的叠加地震道，获得剩余地震记录。
[0021] 进一步地，采用以下公式计算振幅匹配系数：
[0022] 乞 V-、. I i=k-
[0023] 其中，a」为第j道的振幅匹配系数，Ujl为第j道叠前地震记录第i个采样点的振幅 值，Sl为第j道地震记录对应的叠加地震道第i个采样点的振幅值，i b和h为振幅计算时窗的 起止样点序号。
[0024]进一步地，计算所述时移排齐处理后的地震记录中地震道每个采样点对应的相似 权系数，包括：
[0025]计算所述时移排齐处理后的地震记录中地震道的相关系数；
[0026]根据所述地震道的相关系数计算所述地震道中心采样点的相似权系数；
[0027]对所述地震道中心采样点的相似权系数进行归一化处理，并根据所述归一化处理 后的相似权系数计算所述地震道每个采样点对应的相似权系数。
[0028] 进一步地，采用以下公式计算所述地震道的相关系数：
[0029] , = /^,^!< (0' Uj(/), k, j = 1,? ? ?,
[0030] 其中，k和j是滑动窗内地震道的序号，Mi是滑动窗内地震记录道数，i是滑动窗内 的时间样点序号，ib和分别是当前滑动窗第一个和最后一个时间样点序号，u k(i)为第k道 第i个时间样点的振幅值，^⑴为第j道第i个时间样点的振幅值，Cu为第k道和第j道地震 记录的相关系数。
[0031] 进一步地，采用以下公式计算所述地震道中心采样点的相似权系数：
[0033]其中，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数，Ckj为第k道和 第j道的互相关系数，Cw和Ckk分别为第j道和第k道的自相关系数，i。为中心采样点的时间序 号，wib为第k道在i。时间样点处的相似权系数。
[0034]进一步地，采用以下公式对所述地震道中心样点的相似权系数进行归一化处理：
[0036]其中，为第k道在i。时间样点处的相似权系数，为第j道在i。时间样点 处的相似权系数，为第k道在i。时间样点处的归一化后的相似权系数，i。为当前滑动 窗中心采样点的时间序号，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数。 [0037]进一步地，计算所述时移排齐处理后的地震记录中地震道每个采样点对应的能量 权系数，包括：
[0038] 根据所述时移排齐处理后的地震记录获取能量滤波参数；
[0039] 根据所述能量滤波参数获得能量滤波算子；
[0040] 根据所述能量滤波算子计算所述时移排齐处理后的地震记录中地震道中心采样 点的能量权系数；
[0041] 对所述地震道中心采样点的能量权系数进行归一化处理，并根据所述归一化处理 后的能量权系数计算所地震道每个采样点对应的能量权系数。
[0042] 进一步地，根据所述时移排齐处理后的地震记录获取能量滤波参数，包括：
[0043]计算所述时移排齐处理后的地震记录中每一地震道的能量；
[0044] 根据预设的能量比例阈值，将所述时移排齐处理后的地震记录中能量小于阈值的 地震道去除；
[0045] 计算所述时移排齐处理后的地震记录中剩余地震道的平均能量值，并根据所述每 一地震道的能量与所述平均能量值计算所述时移排齐处理后的地震记录的能量标准偏差；
[0046] 根据所述能量标准偏差获取能量滤波参数。
[0047] 进一步地，采用以下公式计算所述能量滤波参数：
[0048] \ " ^ e e$^ef"+l.5rae e4^em + 3-<tc
[0049] 其中，ei、e2、e3和e4表示能量滤波参数，em表示平均能量值，〇e表示能量标准偏差。
[0050] 进一步地，采用以下公式计算所述能量滤波算子： 1 < ec < e3 、 0,5-().5*cos(/T(e. -e,)/(e. -e,)) e： < e(. < e.
[0051 ] f{e,-cn..a；)^< ^ " ! ' 1 ~ 0.5-0,^-cos(/T(ec ~e4)/(e. -e_J) e. < ec < e4 , 〇 ec <ei,ec >e4
[0052] 其中，ei、e2、e3和e4表示能量滤波参数，em表示平均能量值，e c表示地震道能量值。
[0053]进一步地，采用以下公式对所述地震道中心样点的能量权系数进行归一化处理：
[0055] 其中，为第k道在i。时间样点处的能量权系数，Wje)为第j道在i。时间样点 处的能量权系数，〇}ie)为第k道在i。时间样点处的归一化后的能量权系数，i。为当前滑动 窗中心采样点的时间序号，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数。
[0056] 另一方面，本申请实施例还提供了一种地震记录规则干扰波的提取装置，所述装 置包括：
[0057] 获取单元，用于获取预设滑动窗内的叠前地震记录；
[0058]时移排齐处理单元，用于对所述叠前地震记录进行干扰波时移排齐处理，获取时 移排齐处理后的地震记录；
[0059] 计算单元，用于计算所述时移排齐处理后的地震记录中地震道的每个采样点对应 的相似权系数和能量权系数；
[0060] 加权叠加处理单元，用于根据所述地震道的每个采样点对应的相似权系数和能量 权系数，对所述时移排齐处理后的地震记录中地震道进行加权叠加处理，获得所述时移排 齐处理后的地震记录中心道对应的叠加地震道;其中，所述叠加地震道为提取的规则干扰 波；
[0061] 剩余地震记录获取单元，用于根据所述叠前地震记录和所述叠加地震道，获取剩 余地震记录;其中，所述剩余地震记录为去除规则干扰波的地震记录。
[0062] 本申请实施例提供了一种地震记录规则干扰波的提取方法及装置，该方法采用多 道相关的方法计算获得滑动窗内地震道每个样点对应的相似权系数和能量权系数，然后再 利用多相关叠加技术提取规则干扰波，从而可以获得去除规则干扰波的有效的地震记录。 现有技术是利用干扰波和有效波在属性方面的明显差异实现干扰波的提取，本申请实施例 的方法利用了干扰波自身的波形和能量的相关性，即使在规则干扰波和有效波的属性差异 性很小或者没有差异时，也能实现对规则干扰波的提取。本申请实施例的方法能有效的实 现了规则干扰波的提取，为地震资料的解释奠定了基础。
【附图说明】
[0063] 为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提 下，还可以根据这些附图获得其他的附图。
[0064] 图1是本申请实施例的地震记录规则干扰波的提取方法的示意图；
[0065] 图2是本申请实施例的地震记录规则干扰波的提取装置的示意图；
[0066] 图3是本申请一实施例的合成叠前记录的共炮点道集图；
[0067] 图4是本申请一实施例的干扰波时移排齐后的地震记录图；
[0068]图5是本申请一实施例的时移排齐后的地震记录做镶边处理后的图；
[0069] 图6是本申请一实施例的提取的干扰波的图；
[0070] 图7是本申请一实施例的提取的干扰波做反时移校正后的图；
[0071] 图8是本申请一实施例的去除规则干扰波后的剩余地震记录的图；
[0072] 图9是本申请一实施例的野外实际地震记录的图；
[0073] 图10是本申请一实施例的实际地震记录提取的线性干扰波场的图；
[0074] 图11是本申请一实施例的实际地震记录去除线性干扰波场后的剩余地震记录的 图。
【具体实施方式】
[0075] 为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实 施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护 的范围。
[0076]下面结合附图，对本申请实施例的【具体实施方式】作进一步的详细说明。
[0077] 参考图1，本申请实施例提供了一种地震记录规则干扰波的提取方法，该方法包 括：
[0078] S1、获取预设滑动窗内的叠前地震记录。
[0079] 本申请实施例中，首先预设滑动窗的大小，即确定滑动窗内地震道的道数，然后通 过射线追踪正演模拟的方法得到合成的叠前地震记录u r(t)，其中r = l，…，当前滑 动窗内地震记录道数，t为地震记录时间序号，t=l, ...，N，N为地震记录的时间采样点数。 如图3所示为正演模拟的合成叠前地震记录的共炮点道集，其中地震记录道数为240道，时 间采样点数为2000个。在本申请实施例中，叠前地震记录可以为野外的实际地震记录，也可 以为模拟合成的地震记录，本申请实施例中并不限定。
[0080] S2、对所述叠前地震记录进行干扰波时移排齐处理，获取时移排齐处理后的地震 记录。
[0081] 本申请实施例中，对所述叠前地震记录进行干扰波时移排齐处理，获取时移排齐 处理后的地震记录，包括：
[0082]确定所述叠前地震记录的干扰波的时距曲线；
[0083]根据所述干扰波的时距曲线对所述叠前地震记录进行处理，获得时移排齐处理后 的地震记录。
[0084]本申请实施例中，所述干扰波的时距曲线可以根据不同类型的干扰波的特点，采 用如下几种方式确定：
[0085] (1)如果干扰波出现的时间与空间位置无关，则干扰波的时距曲线T (r)可由时间 常量T。来确定，即T(r)=Tc，r为地震道的序号，T。为干扰波在地震记录上出现的时间常量； (2)如果干扰波为线性波场（如直达波），则干扰波的时距曲线可由线性函数
来确定，其中X(r)为序号为r的地震道对应的炮检距，f为干扰波的视速 度，To为零炮检距对应的干扰波的出现时间；（3)如果干扰波为双曲型波场(如反射波），则 干扰波的时距曲线可由双曲函数
来确定，其中X(r)为序号为r的地震 道对应的炮检距，V为干扰波的传播速度，To为零炮检距对应的干扰波在地震记录上出现的 时间；（4)如果干扰波没有上述特征，则可以通过拾取每道地震记录上的干扰波同相轴来确 定干扰波的时距曲线T(r)，其中r为地震道的序号。
[0086]根据确定的干扰波的时距曲线，对每道地震记录ur(t)做时移，每道地震记录的时 移量AT(r) = T(r)-Td，Td为时间常量，时移排齐后的干扰波同相轴被校正到该时间上。在本 申请一实施例汇总，将Td设为150,则时移排齐后的干扰波同相轴被校正到该时间上，如图4 为干扰波时移排齐后的地震记录。
[0087]在本申请另一实施例中，为减小边界效应，在步骤S2后，还包括对时移排齐处理后 的地震记录进行镶边处理。如果采用的滑动窗道数为施，则镶边道数为(Mi-D/2道，也就是 在地震的第一道前面和最后一道后面，分别镶嵌(MrD/2道。地震记录的前(MrD/2道关于 第1道镜像放置形成前边的镶边地震道;地震记录的后(Mi-D/2道关于最后一道镜像放置 形成后边的镶边地震道，图5为对时移排齐后的地震记录做镶边处理的结果。
[0088] S3、计算所述时移排齐处理后的地震记录中地震道的每个采样点对应的相似权系 数和能量权系数。
[0089] 计算所述时移排齐处理后的地震记录中地震道每个采样点对应的相似权系数，包 括：
[0090] 计算所述时移排齐处理后的地震记录计算中地震道的相关系数；
[0091] 根据所述地震道的相关系数计算所述地震道中心采样点的相似权系数；
[0092] 对所述地震道中心采样点的相似权系数进行归一化处理，并根据所述归一化处理 后的相似权系数计算所述地震道每个采样点对应的相似权系数。
[0093] 本申请实施例中，可以采用以下公式计算所述地震道的相关系数：
[0094] CLi： = ^uk(/) ? UjO'), k, j - 1,- ? A/,
[0095]其中，k和j是滑动窗内地震道的序号，Mi是滑动窗内地震记录道数，i是滑动窗内 的时间样点序号，ib和分别是当前滑动窗第一个和最后一个时间样点序号，uk(i)为第k道 第i个时间样点的振幅值，^⑴为第j道第i个时间样点的振幅值，Cu为第k道和第j道地震 记录的相关系数。
[0096]本申请实施例中，在获得地震道的相关系数后，可以采用以下公式计算所述地震 道中心采样点的相似权系数：
[0098] 其中，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数，Ckj为第k道和 第j道的互相关系数，Cw和Ckk分别为第j道和第k道的自相关系数，i。为当前滑动窗中心采样 点的时间序号，wjl ；)为第k道在i。时间样点处的相似权系数。
[0099] 本申请实施例中，可以采用以下公式对所述地震道中心采样点的相似权系数进行 归一化处理：
[0101] 其中，wjb为第k道在i。时间样点处的相似权系数，1为第j道在i。时间样点 处的相似权系数，为第k道在i。时间样点处的归一化后的相似权系数，i。为当前滑动 窗中心采样点的时间序号，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数。
[0102] 本申请实施例中，在时间方向上，滑动窗按照给定的时间步长滑动，逐个计算每一 个滑动窗内每一道时间中心样点对应的相似权系数，在两个时窗的时间中心样点之间，利 用线性插值计算每个样点对应的相似权系数，其中，t为时间样点序号，t = l,…，N，N 为每道地震记录的时间采样点数，k是当前滑动窗内地震道的序号，k=l，…％屬是当前滑 动窗内地震道记录道数。
[0103] 本申请实施例中，计算所述时移排齐处理后的地震记录计算中地震道每个采样点 对应的能量权系数，包括：
[0104] 根据所述时移排齐处理后的地震记录获取能量滤波参数；
[0105] 根据所述能量滤波参数获得能量滤波算子；
[0106] 根据所述能量滤波算子计算所述时移排齐处理后的地震记录中地震道中心样点 的能量权系数；
[0107] 对所述地震道中心采样点的能量权系数进行归一化处理，并根据所述归一化处理 后的能量权系数计算所述地震道每个采样点对应的能量权系数。
[0108] 本申请实施例中，根据所述时移排齐处理后的地震记录获取能量滤波参数，包括：
[0109] 计算所述时移排齐处理后的地震记录中每一地震道的能量；
[0110] 根据预设的能量比例阈值，将所述时移排齐处理后的地震记录中能量小于阈值的 地震道去除；
[0111] 计算所述时移排齐处理后的地震记录中剩余地震道的平均能量值，并根据所述每 一地震道的能量与所述平均能量值计算所述时移排齐处理后的地震记录的能量标准偏差；
[0112] 根据所述能量标准偏差获取能量滤波参数。
[0113] 本申请实施例中，可以采用以下公式计算所述滑动窗的能量标准偏差％，
[0115] 其中，< 表示地震道的能量，1 = 1，~沁，111是当前滑动窗内地震道记录道数，如表 示平均能量值。
[0116] 本申请实施例中，可以采用以下公式计算所述能量滤波参数： ^1 - _ 3 CT 「 n -1-5 <j
[0117] - T c " ei^em+\.5l-(je e4=em + 3-cre
[0118] 其中，ei、e2、e3和e4表示能量滤波参数，em表示平均能量值，〇e表示能量标准偏差。
[0119] 本申请实施例中，可以采用以下公式计算所述能量滤波算子： 1 e2<ec <e3
[0120] ^ ' 0.5- ()_5- cos(;r(e.」e4)/(，,- e4)) e-; < ec < e_t 、 0 ec <el5ec >e4
[0121] 其中，f (ec^em^e)表示能量滤波算子，ei、e2、e3和e4表示能量滤波参数，e。表示当前 地震道的能量，~表示能量标准偏差。
[0122] 本申请实施例中，可以采用以下公式计算所述滑动窗内地震道中心样点的能量权 系数：
[0123] V^:.)=/(d〇；.) k = ~-(k.)
[0124] 其中，为第k道在ic时间样点处的能量权系数，f(ec;em，〇e)表示能量滤波算 子，i。为当前滑动窗中心样点的时间序号，k为滑动窗内地震道的序号，Mi是当前滑动窗内地 震记录道数。
[0125] 本申请实施例中，可以采用以下公式对所述地震道中心采样点的能量权系数进行 归一化处理：
[0127] 其中，Wefe)为第k道在i。时间样点处的能量权系数，W^}为第j道在i。时间样点 处的能量权系数，^为第k道在i。时间样点处的归一化后的能量权系数，i。为当前滑动 窗中心采样点的时间序号，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数。
[0128] 本申请实施例中，在时间方向上，滑动窗按照给定的时间步长滑动，逐个计算每一 个滑动窗内每一道时间中心样点对应的能量权系数，在两个时窗的时间中心样点之间，利 用线性插值方法计算每个采样点对应的能量权系数其中，t为时间样点序号，t = 1，…，N，N为每道地震记录的时间采样点数，k是当前滑动窗内地震道的序号， 是当前滑动窗内地震道总道数。
[0129] S4、根据所述地震道的每个采样点对应的相似权系数和能量权系数，对所述时移 排齐处理后的地震记录中地震道进行加权叠加处理，获得所述时移排齐处理后的地震记录 中心道对应的叠加地震道;其中，所述叠加地震道为提取的规则干扰波。
[0130] 本申请实施例中，可以采用以下计算公式获得所述时移排齐处理后地震记录中心 道对应的叠加地震道： ik) ik)
[0131] s(r) = (，)价:(r)%,(，）t=-\，…N k=l
[0132] 其中，s(t)表示叠加地震道；表示地震道每个采样点对应的相似权系数； 表示地震道每个采样点对应的能量权系数;Uk(t)表示时移排齐处理后的地震记录上 第k道地震道在时间样点序号为t时的振幅值;k表示地震道的序号;Mi表示滑动窗内地震记 录的道数;t表示时间样点序号，N为每道地震记录的时间采样点数。
[0133] 根据地震波理论可知，地震记录来自同一激发源的规则干扰波，除了在视速度、偏 振特性和频率等属性上具有相关性外，不同道上的干扰波的波形和能量也具有相关性。本 申请实施例根据规则干扰波的这一特点，采用多道相关的方法计算出每道的波形相似权系 数和能量权系数，然后应用多相关加权叠加技术从而可以有效的提取规则干扰波。如图6所 示为本申请一实施例得到的叠加地震记录，即需要提取的干扰波。
[0134] S5、根据所述叠前地震记录和所述叠加地震道，获取剩余地震记录;其中，所述剩 余地震记录为去除规则干扰波的地震记录。
[0135] 本申请实施例中，根据所述叠前地震记录和所述叠加地震道，获取剩余地震记录， 包括：
[0136] 对所述叠加地震道进行反时移校正处理，获得反时移校正处理后的叠加地震道；
[0137] 根据所述叠前地震记录和所述反时移校正处理后的叠加地震道，获取振幅匹配系 数；
[0138] 根据所述振幅匹配系数对所述叠加地震道进行振幅匹配处理，利用所述叠前地震 记录减去振幅匹配处理后的叠加地震道，获得剩余地震记录。
[0139] 本申请实施例中，可以根据步骤S2中确定的干扰波的时距曲线，对所示叠加地震 道进行反时移校正，获得反时移校正处理后的叠加地震道。其中，时移量为A T(r)=Td_T (r)，本申请实施例中，Td为150,如图7所示为提取的干扰波做反时移校正后的结果。
[0140] 本申请实施例中，可以采用以下公式计算振幅匹配系数：
[0141] ?；/ /乞'.， i =k / i-=ib
[0142] 其中，a」为第j道的振幅匹配系数，Ujl为第j道叠前地震记录第i个采样点的振幅 值，Sl为第j道地震记录对应的叠加地震道第i个采样点的振幅值，i b和h为振幅计算时窗的 起止样点序号。
[0143] 本申请实施例中，根据所述振幅匹配系数对所述叠加地震道进行振幅匹配处理， 然后用叠前地震记录减去振幅匹配处理后的叠加地震道，得到剩余地震记录，可以采用以 下公式计算剩余地震记录：
[0144] rj(t) =uj(t)-a ? s(t),tb^t^te
[0145] 其中，a表示振幅匹配系数，s(t)表示叠加地震道，Uj(t)表示叠前地震记录，〇(t) 表示剩余地震记录，j是滑动窗内地震道的序号，t b和te为地震记录相减应用时窗的起止时 间。这样，就实现了将规则干扰波从叠前地震记录中分离出来的目的。如图8所述为将规则 干扰波去除后的剩余地震记录。
[0146] 在本申请一实施例中，在步骤S5之后，还包括：
[0147] 逐步减小所述滑动窗的大小，对所述剩余地震记录进行规则干扰波去除处理。
[0148] 本申请实施例中干扰波的提取是在预设的滑动窗内进行的，在开始进行干扰波提 取时可以采用较大的滑动窗，从而将叠前地震记录上的干扰波的主要能量提取出来，然后 再逐步减小滑动窗的大小，将剩余地震记录中的干扰波的残余能量逐步地提取出来，从而 得到质量更好的有效波。如图9所示为野外实际地震记录，图10所示为本申请实施例的方法 提取的线性干扰波场，如图11所示为本申请实施例的方法得到的去除线性干扰波后的剩余 地震记录，从图中可以看出，本申请实施例可以有效的提取干扰波，从而为地震资料解释奠 定基础。
[0149] 另一方面，参考图2,本申请实施例还提供了一种地震记录规则干扰波的提取装 置，所述装置包括：
[0150] 获取单元21，用于获取预设滑动窗内的叠前地震记录；
[0151]时移排齐处理单元22,用于对所述叠前地震记录进行干扰波时移排齐处理，获取 时移排齐处理后的地震记录；
[0152]计算单元23,用于计算所述时移排齐处理后的地震记录中地震道的每个采样点对 应的相似权系数和能量权系数；
[0153]加权叠加处理单元24,用于根据所述地震道的每个采样点对应的相似权系数和能 量权系数，对所述时移排齐处理后的地震记录中地震道进行加权叠加处理，获得所述时移 排齐处理后的地震记录中心道对应的叠加地震道;其中，所述叠加地震道为提取的规则干 扰波；
[0154] 剩余地震记录获取单元25,用于根据所述叠前地震记录和所述叠加地震道，获取 剩余地震记录;其中，所述剩余地震记录为去除规则干扰波的地震记录。
[0155] 本实施例的装置的各组成部分分别用于实现前述实施例的方法的各步骤，由于在 方法实施例中，已经对各步骤进行了详细说明，在此不再赘述。
[0156] 本申请实施例提供了一种地震记录规则干扰波的提取方法及装置，该方法采用多 道相关的方法计算获得滑动窗内地震道每个样点对应的相似权系数和能量权系数，然后再 利用多相关叠加技术及振幅平滑相减技术提取规则干扰波，从而获得去除规则干扰波的有 效的地震记录。现有技术是利用规则干扰波和有效波在属性方面的差异性实现对规则干扰 波的提取，当规则干扰波和有效波的属性差异很小或没有差异时，现有技术的方法就不能 对干扰波进行有效提取或压制;而本申请实施例的方法是利用了干扰波自身的波形和能量 的相关性，因此即使规则干扰波和有效波的属性差异很小或没有差异，也可以利用干扰波 自身波形和能量的相关性提取出规则干扰波。本申请实施例的方法能有效的实现对干扰波 的提取，从而为地震资料的解释奠定了基础。
[0157] 在一个或多个示例性的设计中，本申请实施例所描述的上述功能可以在硬件、软 件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的 媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电 脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以 是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但 不限于1^11、1?(观、￡￡?1?(通、^-1?(通或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任 何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理 器读取形式的程序代码的媒介。
[0158] 以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本申请实施例的具体实施例而已，并不用于限定本申 请的保护范围，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包 含在本申请的保护范围之内。
【主权项】
1. 一种地震记录规则干扰波的提取方法，其特征在于，所述方法包括： 获取预设滑动窗内的叠前地震记录； 对所述叠前地震记录进行干扰波时移排齐处理，获取时移排齐处理后的地震记录； 计算所述时移排齐处理后的地震记录中地震道的每个采样点对应的相似权系数和能 量权系数； 根据所述地震道的每个采样点对应的相似权系数和能量权系数，对所述时移排齐处理 后的地震记录中地震道进行加权叠加处理，获得所述时移排齐处理后的地震记录中屯、道对 应的叠加地震道;其中，所述叠加地震道为提取的规则干扰波； 根据所述叠前地震记录和所述叠加地震道，获取剩余地震记录;其中，所述剩余地震记 录为去除规则干扰波的地震记录。2. 如权利要求1所述的方法，其特征在于，在根据所述叠前地震记录和所述叠加地震 道，获取剩余地震记录之后，所述方法还包括： 逐步减小所述滑动窗的大小，对所述剩余地震记录进行规则干扰波提取处理。3. 如权利要求1所述的方法，其特征在于，采用W下计算公式获得所述时移排齐处理后 地震记录中屯、道对应的叠加地震道：其中，s(t)表示叠加地震道;表示地震道每个采样点对应的相似权系数;表 示地震道每个采样点对应的能量权系数;Uk(t)表示时移排齐处理后的地震记录上第k道地 震道在时间样点序号为t时的振幅值;k表示地震道的序号;Mi表示滑动窗内地震记录的道 数;t表示时间样点序号，N为每道地震记录的时间采样点数。4. 如权利要求1所述的方法，其特征在于，根据所述叠前地震记录和所述叠加地震道， 获取剩余地震记录，包括： 对所述叠加地震道进行反时移校正处理，获得反时移校正处理后的叠加地震道； 根据所述叠前地震记录和所述反时移校正处理后的叠加地震道，获取振幅匹配系数； 根据所述振幅匹配系数对所述叠加地震道进行振幅匹配处理，利用所述叠前地震记录 减去振幅匹配处理后的叠加地震道，获得剩余地震记录。5. 如权利要求4所述的方法，S賠化#干.要田W下/A式计算振幅匹配系数：其中，W为第j道的振幅匹配系数，UW为第j道叠前地震记录第i个采样点的振幅值，Si为 第j道地震记录对应的叠加地震道第i个采样点的振幅值，ib和ie为振幅计算时窗的起止样 点序号。6. 如权利要求1所述的方法，其特征在于，计算所述时移排齐处理后的地震记录中地震 道每个采样点对应的相似权系数，包括： 计算所述时移排齐处理后的地震记录中地震道的相关系数； 根据所述地震道的相关系数计算所述地震道中屯、采样点的相似权系数； 对所述地震道中屯、采样点的相似权系数进行归一化处理，并根据所述归一化处理后的 相似权系数计算所述地震道每个采样点对应的相似权系数。7. 如权利要求6所述的方法，其特征在于，采用W下公式计算所述地震道的相关系数：其中，k和j是滑动窗内地震道的序号，Mi是滑动窗内地震记录道数，i是滑动窗内的时间 样点序号，ib和ie分别是当前滑动窗第一个和最后一个时间样点序号，Uk(i)为第k道第i个 时间样点的振幅值，11如）为第j道第i个时间样点的振幅值，Cm为第k道和第j道地震记录 的相关系数。8. 如权利要求6所述的方法，其特征在于，采用W下公式计算所述地震道中屯、采样点的 相似权系数：其中，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数，CkJ为第k道和第j 道的互相关系数，Cw和Ckk分别为第j道和第k道的自相关系数，ic为当前滑动窗中屯、采样点 的时间序号，^\^]^为第4道在1。时间样点处的相似权系数。9. 如权利要求6所述的方法，其特征在于，采用W下公式对所述地震道中屯、采样点的相 似权系数进行归一化处理：其中，W:策：)为第k道在1。时间样点处的相似权系数，\滿。；> 为第j道在1。时间样点处的 相似权系数，为第k道在i。时间样点处的归一化后的相似权系数，i。为当前滑动窗中 屯、采样点的时间序号，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数。10. 如权利要求1所述的方法，其特征在于，计算所述时移排齐处理后的地震记录中地 震道每个采样点对应的能量权系数，包括： 根据所述时移排齐处理后的地震记录获取能量滤波参数； 根据所述能量滤波参数获得能量滤波算子； 根据所述能量滤波算子计算所述时移排齐处理后的地震记录中地震道中屯、采样点的 能量权系数； 对所述地震道中屯、采样点的能量权系数进行归一化处理，并根据所述归一化处理后的 能量权系数计算所述地震道每个采样点对应的能量权系数。11. 如权利要求10所述的方法，其特征在于，根据所述时移排齐处理后的地震记录获取 能量滤波参数，包括： 计算所述时移排齐处理后的地震记录中每一地震道的能量； 根据预设的能量比例阔值，将所述时移排齐处理后的地震记录中能量小于阔值的地震 道去除； 计算所述时移排齐处理后的地震记录中剩余地震道的平均能量值，并根据所述每一地 震道的能量与所述平均能量值计算所述时移排齐处理后的地震记录的能量标准偏差； 根据所述能量标准偏差获取能量滤波参数。12.如权利要求11所述的方法，其据佈在干，采巧Pi下公式计算所述能量滤波参数： 其中，ei、62、63和64表示能量滤波参数，em表示平均能量值，Oe表示能量标准偏差。 13 .如权利要求10所述的方法，其特征在于，采用W下公式计算所述能量滤波算子：其中，ei、62、63和64表示能量滤波参数，Gm表示平均能量值，Gc表示地震道能量值。14. 如权利要求10所述的方法，其特征在于，采用W下公式对所述地震道中屯、采样点的 能量权系数进行归一化处理：其中，为第k道在1。时间样点处的能量权系数，为第j道在1。时间样点处的 能量权系数，\\5；^；)为第4道在1。时间样点处的归一化后的能量权系数，1。为当前滑动窗中 屯、采样点的时间序号，k为滑动窗内地震道的序号，Mi是当前滑动窗内地震记录道数。15. -种地震记录规则干扰波的提取装置，其特征在于，所述装置包括： 获取单元，用于获取预设滑动窗内的叠前地震记录； 时移排齐处理单元，用于对所述叠前地震记录进行干扰波时移排齐处理，获取时移排 齐处理后的地震记录； 计算单元，用于计算所述时移排齐处理后的地震记录中地震道的每个采样点对应的相 似权系数和能量权系数； 加权叠加处理单元，用于根据所述地震道的每个采样点对应的相似权系数和能量权系 数，对所述时移排齐处理后的地震记录中地震道进行加权叠加处理，获得所述时移排齐处 理后的地震记录中屯、道对应的叠加地震道;其中，所述叠加地震道为提取的规则干扰波； 剩余地震记录获取单元，用于根据所述叠前地震记录和所述叠加地震道，获取剩余地 震记录;其中，所述剩余地震记录为去除规则干扰波的地震记录。
【文档编号】G01V1/28GK105911585SQ201610522410
【公开日】2016年8月31日
【申请日】2016年7月5日
【发明人】张文波, 李合群, 张文栋
【申请人】中国石油集团东方地球物理勘探有限责任公司, 中油油气勘探软件国家工程研究中心有限公司