提高复杂构造偏移数据信噪比的方法
【技术领域】
[0001] 本发明设及一种提高复杂构造偏移数据信噪比的方法，用于石油地震勘探，属于 地震勘探资料处理与解释领域，是一种对低信噪比资料中的倾角相干噪声进行压制，提高 地震资料的信噪比，能够为后续解释提供质量更高的地震记录。
【背景技术】
[0002] 通常，复杂山地地震资料噪声来源十分复杂，噪声对有效波有较强的干扰，尤其是 各种类型的规则干扰和不规则干扰发育强烈，导致地震资料的信噪比低，严重影响了地震 数据的预处理和偏移处理效果，也给后续的解释带来了困扰。
[0003] 例如，通常随着油气勘探由平原转向丘陵、山地等复杂构造地区，运些地区地表高 程起伏大，低降速层的速度和厚度变化大，地层产状在横向和纵向上变化显著，而在运些地 区采集的地震资料与平原地震资料相比，具有各种类型的规则干扰和不规则干扰发育强 烈，资料的信噪比很低，即使经过前期噪声压制的预处理，但在偏移处理中还会残留大量的 线性干扰，影响了偏移的正确成像，低信噪比的偏移剖面也给后续的解释带来了困扰。
[0004] 由于复杂采集条件的影响，偏移剖面上往往存在预处理无法完全压制的睹倾角相 干噪声，而运些睹倾角相干噪声的倾角有时与有效信号的倾角非常接近，因此用常规的F-X 域滤波、F-K域滤波或K-L滤波等方法无法达到满意的噪声压制效果，因为运些方法在处理 过程中用同样的参数对各个角度的倾角相干噪声进行压制，不能进行空变，在压制噪声的 同时，也会牺牲部分倾角构造的有效信号，运显然违背了噪声压制的初衷。
【发明内容】
[0005] 本发明提出的一种提高复杂构造偏移数据信噪比的方法，能够实现倾角相干噪声 压制效果与倾角有效信号特征保持的最佳折中。
[0006] 为了实现上述目的，本发明提供了一种提高复杂构造偏移数据信噪比的方法。所 述方法包括步骤:A、采集地震数据，经过预处理、偏移处理偏移后得到偏移剖面;B、将偏移 剖面的初至W前时间样点设置为零值，并增大数据时间样点数;C、对偏移数据进行构造层 位文件拾取，并进行全Ξ维层位线性插值处理;D、给定校正拉平处理时间点值，计算层位每 道时间与其给定时间点的差值，作为每道拉平校正的校正量，对整个数据进行拉平处理;E、 对步骤D的数据进行分段倾角滤波压制相干噪声处理;F、对步骤E中分段处理数据进行分段 时变混波处理;G、对步骤F中分段数据进行合并，并按道和偏移距重新抽取数据，并进行校 正反拉平处理，得到压制相干噪声后偏移数据。
[0007] 与现有技术相比，本发明的有益效果包括:通过对偏移剖面中拾取构造层位，通过 层位插值对数据进行校正拉平处理，使得有效信号变成近似水平同相轴，而相干噪声还保 留一定的倾角，通过分段选取数据，分段进行倾角滤波压制相干噪声，就能很好保护水平状 有效信号，并通过时变混波处理增强其保真性，最后通过数据的合并抽取恢复为去噪后的 偏移数据，实现倾角相干噪声压制效果与倾角有效信号特征保持的最佳折中；提高了复杂 构造偏移数据信噪比的处理效果。
【附图说明】
[0008] 通过下面结合附图进行的描述，本发明的上述和其他目的和特点将会变得更加清 楚，其中：
[0009] 图1示出了实际地震数据偏移剖面；
[0010] 图2示出了采用本发明的一个示例性实施例对图1中的实际地震数据偏移剖面进 行去噪处理后的偏移剖面；
[0011] 图3示出了采用常规方法对图1的实际地震数据偏移剖面进行去噪处理后的偏移 剖面。
【具体实施方式】
[0012] 在下文中，将结合示例性实施例来详细说明本发明的提高复杂构造偏移数据信噪 比的方法。
[0013] 本发明旨在提出一种能够针对复杂构造偏移数据相干噪声进行压制，进而提高其 信噪比的方法。
[0014] 本发明的核屯、内容是，对偏移剖面中拾取构造层位，通过层位插值对数据进行校 正拉平处理，使得有效信号变成近似水平同相轴，而相干噪声还保留一定的倾角，通过分段 选取数据，分段进行倾角滤波压制相干噪声，就能很好保护水平状有效信号，并通过时变混 波处理增强其保真性，最后通过数据的合并抽取恢复为去噪后的偏移数据，实现倾角相干 噪声压制效果与倾角有效信号特征保持的最佳折中。
[0015] 在本发明的一个示例性实施例中，提高复杂构造偏移数据信噪比的方法可通过W 下步骤来实现：
[0016] (1)采集地震数据，经过预处理、偏移处理偏移后得到偏移剖面。运里，偏移剖面可 W为时间偏移剖面或者深度偏移剖面。
[0017] 本步骤可按下列方式实施：
[0018] 对野外采集到的地震数据，经过静校正、叠前去噪、振幅补偿、速度分析等预处理 后，利用偏移算法(例如，时间偏移或速度偏移)对地震数据进行偏移归位，得到偏移剖面。 运里，静校正、叠前去噪、振幅补偿、速度分析等偏移之前的预处理、W及偏移算法属于常规 操作。
[0019] (2)将偏移剖面的初至时间及初至时间之前的时间样点数的振幅值均设置为零 值，并增大数据时间样点数，且增加的时间样点数的振幅值均为零。例如，增大的数据时间 样点数可W为未增加前数据样点总数的1/7至1/4,其跟随在未增加前数据之后。增大样点 数能够有效避免后面的数据拉平处理时可能造成的数据丢失。
[0020] 本步骤可按下列方式实施：
[0021] (a)由于初至时间采样点如果有非零值，则会对后续的噪声压制产生异值影响，所 W必须对初至W前时间样点设置为零。
[0022] 对偏移数据按道计算数据，当数据振幅值在-0.000001与0.000001之间，则可判断 为初至前样点值，则把该振幅值赋值为零，完成所有道的计算，所有道初至前样点值赋值为 零。
[0023] (b)由于后续校正拉平处理过程中，可能存在拉平数据超出数据本身样点数，而损 失样点值的风险，所W必须在校正前，至少增加1000个样点数，而增加的1000样点值的振幅 值赋值为零。例如，数据本身为5s数据，增加 Is数据变为6s数据，5s-6s之间数据的振幅值均 赋值为零。
[0024] (3)对偏移数据进行构造层位(即，层位时间值)文件拾取，W拾取拟保护的复杂构 造同相轴层位（即，同相轴层位时间值），并进行全Ξ维层位线性插值处理。对于Ξ维(3D)偏 移数据，需要拾取两条或Ξ条W上主测线方向上的层位值，并对已拾取的多条主测线(两条 或Ξ条W上的主测线）中所有的彼此相邻的两条主测线进行线间线性插值处理，随后对每 一条主测线进行本主测线线内层位线性插值处理。对于二维(2D)偏移数据，则仅需对该2D 数据的测线进行测线内部的层位线性插值处理即可。
[0025] 本步骤可按下列方式实施：
[0026] (a)对偏移数据进行构造层位文件拾取，拾取的原则为拾取需要保护的复杂构造 同相轴层位，W免后续的噪声压制伤害该同相轴。如果数据为3D数据，则需要多拾取几条主 测线方向上的层位值，例如，拾取所有主测线方向上的层位值。
[0027] (b)如果偏移数据为2D数据，则只需对拾取的层位值进行线性插值处理，插值出该 层位中每一道的层位值，用于后续计算其校正量。
[0028] 如果偏移数据为3D数据，则首先利用已拾取层位值的两条主测线进行线间线性插 值处理，直到所有的主测线均插值有同拾取层位个数相同的层位点值，然后对每条主测线 进行