山前带地震去噪方法技术

公开了一种山前带地震去噪方法，包括：步骤1：对原始地震数据分类，获得多个地震数据组与属性分布；步骤2：计算噪音特征，获得噪音特征参数；步骤3：基于噪音特征参数，获得去噪因子与去噪参数组；步骤4：建立数据矩阵并求解信噪比因子，获得优化地震数据；步骤5：将优化地震数据作为原始地震数据，重复执行步骤1‑4，获得优化信噪比因子与再优化地震数据；步骤6：比较信噪比因子与优化信噪比因子，如果差大于预定值，返回步骤5将再优化地震数据作为原始地震数据，如果差小于预定值，输出再优化地震数据为去噪地震数据。本发明专利技术通过递进迭代消除噪音，提高数据的信号能量，具有较强的稳定性和适应性，提高了信噪比。

【技术实现步骤摘要】
山前带地震去噪方法
本专利技术涉及石油地球物理勘探
，更具体地，涉及一种山前带地震去噪方法。
技术介绍
在山前带油气勘探工程中，提高地震数据信噪比的问题是整个地震处理过程的关键问题。常规信噪分离方法主要是依据信号和噪音的比值(S/N)。为此提高信噪比的方案包括在噪音不变的基础上，放大增强信号；再一种就是在信号不变的前提下，缩小噪音值。或者，同时放大信号与缩小噪音值，也能达到提高数据信噪比的目的。在常规地区，这种去噪方法相对成熟，对干扰特征明显的干扰和噪音，去除效果较好。但是，山前带干扰因素复杂多变，地震信号的激发和接收千差万别，各种干扰源以及次生干扰影响因素多而杂，有效信号的能量非常弱，频带完全被噪音淹没，分离非常困难。因此，有必要开发一种适用于山前带的地震去噪方法。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种山前带地震去噪方法，其根据山前带地震原始数据的特点，与高精度静校正方法结合，递进迭代消除噪音，提高数据的信号能量，具有较强的稳定性和适应性，提高了山前带地震数据的信噪比。根据本专利技术的一方面，提出了一种山前带地震去噪方法。所述方法可以包括：步骤1：对原始地震数据进行分类，获得多个地震数据组，进而获得所述多个地震数据组的属性分布；步骤2：基于所述多个地震数据组的属性分布，计算所述多个地震数据组的噪音特征，获得每个地震数据组的噪音特征参数；步骤3：针对每个地震数据组，基于所述地震数据组的噪音特征参数，获得所述地震数据组对应的去噪因子，并将所述去噪因子与所述原始地震数据进行互相关运算，获得所述地震数据组对应的去噪参数组；步骤4：基于所述原始地震数据、所述多个地震数据组的属性分布、噪音特征参数与去噪参数组，建立数据矩阵并求解，获得信噪比因子，将所述信噪比因子应用于所述原始地震数据，获得优化地震数据；步骤5：将所述优化地震数据作为原始地震数据，重复执行步骤1-4，获得优化信噪比因子，进而获得再优化地震数据；步骤6：比较所述信噪比因子与所述优化信噪比因子，如果所述信噪比因子与所述优化信噪比因子之差大于预定值，则返回步骤5并将所述再优化地震数据作为原始地震数据，如果二者之差小于所述预定值，则输出所述再优化地震数据为去噪地震数据。优选地，还包括：判断所述信噪比因子是否在预定范围内，如果不在所述预定范围内，则排查所述多个地震数据组的噪音特征的异常点并删除。优选地，所述噪音特征包括：振幅、能量、速度与频率。优选地，对原始地震数据进行分类包括，获得多个地震数据组：根据所述原始地震数据的噪音属性，将所述原始地震数据分为随机干扰与规则干扰，根据振幅、频率、速度与影响范围，对所述规则干扰进一步划分，获得所述多个地震数据组。优选地，所述多个地震数据组包括：线性干扰组、面波干扰组、浅层多次干扰组、随机干扰组。优选地，所述信噪比因子为：其中，Sn表示信噪比因子，F(x)表示原始地震数据的地震数据组，A(x)表示属性分布，B(x)表示噪音特征参数，N(x)表示去噪参数组，n表示地震数据组的个数，x表示原始地震数据。本专利技术的有益效果在于：根据山前带地震原始数据的特点，与高精度静校正方法结合，递进迭代消除噪音，方法可靠，计算思路明确，能够较好的控制山前带单炮信噪比的变化，避免了单一去噪方法的局限性，计算易于收敛，结果可靠。本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的山前带地震去噪方法的步骤的流程图。图2示出了根据本专利技术的一个实施例的山前带工区近地表高程变化的示意图。图3示出了根据本专利技术的一个实施例的近地表地质模型的示意图。图4a、图4b、图4c和图4d分别示出了根据本专利技术的一个实施例的山前带平原区地形的面波、过渡带地形的线性波、过渡带地形的浅层多次波和山上地形的随机噪音的原始单炮数据和对应高度的示意图。图5a、图5c和图5e分别示出了根据本专利技术的一个实施例的山前带平原区地形、过渡带地形和山上地形的原始单炮数据的示意图，图5b、图5d和图5f分别示出了对应的信噪比的示意图。图6示出了根据本专利技术的一个实施例的山前带信噪比平面分布的示意图。图7a和图7b分别示出了根据本专利技术的一个实施例的原始单炮数据和信号与噪音分离频谱的示意图。图8a、图8b和图8c分别示出了平原区地形的原始地震数据与根据本专利技术的一个实施例的去噪地震数据及去除的噪音的示意图。图9a、图9b和图9c分别示出了过渡带地形的原始地震数据与根据本专利技术的一个实施例的去噪地震数据及去除的噪音的示意图。图10a、图10b和图10c分别示出了山上地形的原始地震数据与根据本专利技术的一个实施例的去噪地震数据及去除的噪音的示意图。图11a和图11b分别示出了根据本专利技术的一个实施例的共中心点道集与去噪后的共中心点道集的示意图。图12a和图12b分别示出了根据本专利技术的一个实施例的共中心点叠加剖面与去噪后的共中心点叠加剖面的示意图。图13a和图13b分别示出了现有技术的去噪地震数据与根据本专利技术的一个实施例的去噪地震数据的示意图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。图1示出了根据本专利技术的山前带地震去噪方法的步骤的流程图。根据本专利技术的山前带地震去噪方法可以包括：步骤1：对原始地震数据进行分类，获得多个地震数据组，进而获得多个地震数据组的属性分布，包括地震数据组的信噪比水平、信号能量与噪音类型。在一个示例中，对原始地震数据进行分类包括，获得多个地震数据组：根据原始地震数据的噪音属性，将原始地震数据分为随机干扰与规则干扰，根据振幅、频率、速度与影响范围，对规则干扰进一步划分，进而根据不同的干扰类型划分多个地震数据组。在一个示例中，多个地震数据组包括：线性干扰组、面波干扰组、浅层多次干扰组、随机干扰组。具体地，在地震勘探中，接收到的地震数据中叠合了各种干扰，只有包含地下地层的反射信息的波信号，被称为有效信号，其它所有妨碍有效信号辨识和追踪的地震波，都被称为噪音。由于工程采集条件以及地下地质的差别，使得接收到的噪音具有多样性。由于噪音和有效信号叠合在一起，不能一次性把噪音从地震波中去除，只有根据不同噪音的不同特征，采用有针对性的手段，分别衰减或者消除，并不损失有效信号。根据噪音的特点在时间上随机性和相干性的特点，将具有随机性的特征噪音称为随机干扰，具有规则相干性的噪音称为规则干扰，根据噪音在振幅、频率、速度等特性的差异以及影响的范围的不同，对规则干扰进一步划分，可以包括线性波、面波、浅层多本文档来自技高网...

【技术保护点】
1.一种山前带地震去噪方法，包括：步骤1：对原始地震数据进行分类，获得多个地震数据组，进而获得所述多个地震数据组的属性分布；步骤2：基于所述多个地震数据组的属性分布，计算所述多个地震数据组的噪音特征，获得每个地震数据组的噪音特征参数；步骤3：针对每个地震数据组，基于所述地震数据组的噪音特征参数，获得所述地震数据组对应的去噪因子，并将所述去噪因子与所述原始地震数据进行互相关运算，获得所述地震数据组对应的去噪参数组；步骤4：基于所述原始地震数据、所述多个地震数据组的属性分布、噪音特征参数与去噪参数组，建立数据矩阵并求解，获得信噪比因子，将所述信噪比因子应用于所述原始地震数据，获得优化地震数据；步骤5：将所述优化地震数据作为原始地震数据，重复执行步骤1‑4，获得优化信噪比因子，进而获得再优化地震数据；步骤6：比较所述信噪比因子与所述优化信噪比因子，如果所述信噪比因子与所述优化信噪比因子之差大于预定值，则返回步骤5并将所述再优化地震数据作为原始地震数据，如果二者之差小于所述预定值，则输出所述再优化地震数据为去噪地震数据。

【技术特征摘要】
1.一种山前带地震去噪方法，包括：步骤1：对原始地震数据进行分类，获得多个地震数据组，进而获得所述多个地震数据组的属性分布；步骤2：基于所述多个地震数据组的属性分布，计算所述多个地震数据组的噪音特征，获得每个地震数据组的噪音特征参数；步骤3：针对每个地震数据组，基于所述地震数据组的噪音特征参数，获得所述地震数据组对应的去噪因子，并将所述去噪因子与所述原始地震数据进行互相关运算，获得所述地震数据组对应的去噪参数组；步骤4：基于所述原始地震数据、所述多个地震数据组的属性分布、噪音特征参数与去噪参数组，建立数据矩阵并求解，获得信噪比因子，将所述信噪比因子应用于所述原始地震数据，获得优化地震数据；步骤5：将所述优化地震数据作为原始地震数据，重复执行步骤1-4，获得优化信噪比因子，进而获得再优化地震数据；步骤6：比较所述信噪比因子与所述优化信噪比因子，如果所述信噪比因子与所述优化信噪比因子之差大于预定值，则返回步骤5并将所述再优化地震数据作为原始地震数据，如果二者之差小于所述预定值，则输出所述再优化地震...

【专利技术属性】
技术研发人员：孙志锋，刘思思，蒋波，高艳霞，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11