CN108363103B - 规则干扰压制方法及装置 - Google Patents
规则干扰压制方法及装置 Download PDFInfo
- Publication number
- CN108363103B CN108363103B CN201810101694.6A CN201810101694A CN108363103B CN 108363103 B CN108363103 B CN 108363103B CN 201810101694 A CN201810101694 A CN 201810101694A CN 108363103 B CN108363103 B CN 108363103B
- Authority
- CN
- China
- Prior art keywords
- apparent velocity
- time difference
- rule
- big gun
- regular
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000001629 suppression Effects 0.000 title claims abstract description 50
- 230000001788 irregular Effects 0.000 claims abstract description 24
- 238000005070 sampling Methods 0.000 claims abstract description 20
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 238000004590 computer program Methods 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 5
- 239000002131 composite material Substances 0.000 claims description 3
- 244000089409 Erythrina poeppigiana Species 0.000 claims 4
- 235000009776 Rathbunia alamosensis Nutrition 0.000 claims 4
- 239000000284 extract Substances 0.000 claims 2
- 238000012545 processing Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 238000003860 storage Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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
- G01V1/364—Seismic filtering
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)
Abstract
本申请实施例提供了一种规则干扰压制方法及装置，该方法包括：获取非规则采样地震数据的炮集记录；确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；将所述规则干扰区域数据变换为频率空间域地震数据；确定频率空间域规则干扰滤波算子；将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。本申请实施例不但能压制规则干扰，还能更好地保持振幅特征，从而能够达到在压制规则干扰的同时，提高地震数据信噪比和分辨率的目的。
Description
技术领域
本申请涉及物探地震数据处理技术领域，尤其是涉及一种规则干扰压制方法及装置。
背景技术
在物探地震数据处理技术中，规则干扰是影响地震数据叠加成像质量的一个重要因数，这类叠加剖面上观察到的规则噪音通常情况下在地震数据叠加前的炮记录或者共检波点记录上都能观察到。
地震资料中对这类噪音的压制通常是在地震数据叠前的炮记录上进行，常用的衰减办法是在FK根据用户选择的一个规则干扰的视速度范围进行切除压制，炮记录上处于这个视速度范围内的相干噪音就被衰减了，成像剖面上看不到规则干扰，这样地震数据信噪比得到了有效提高。但是这种方法不适应空间非规则采样的情况，噪声衰减后会产生不同程度的空间假频，影响叠加成像，甚至带来假象，给地震数据的后期处理、属性反演、油藏预测带来无法估量的损失。
发明内容
本申请实施例的目的在于提供一种规则干扰压制方法及装置，以实现在压制规则干扰的同时，还可以提高有用信号的保真度。
为达到上述目的，一方面，本申请实施例提供了一种规则干扰压制方法，包括：
获取非规则采样地震数据的炮集记录；
确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；
根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；
将所述规则干扰区域数据变换为频率空间域地震数据；
确定频率空间域规则干扰滤波算子；
将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；
根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
较佳的，所述确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差，包括：
确定所述炮集记录中横跨特定数量地震道的规则干扰波的最大视速度和最小视速度；
将所述规则干扰波的最大视速度和最小视速度分别除以所述规则干扰波所横跨地震道的特定数量，对应获得所述炮集记录中相邻地震道的视速度最大平均时差和视速度最大平均时差。
较佳的，所述根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据，包括：
在所述炮集记录上识别出规则干扰和有效信号的界线；
判断所述界线是否大于所述视速度最大平均时差或小于所述视速度最小平均时差；
如果述界线大于所述视速度最大平均时差或小于所述视速度最小平均时差，则根据所述非规则采样地震数据在所述炮集记录上的时间分布和空间范围以及所述界线，将所述炮集记录划分为规则干扰区域数据和非规则干扰区域数据，并从所述炮集记录中提取出所述规则干扰区域数据。
较佳的，所述将所述规则干扰区域数据变换为频率空间域地震数据，包括：
通过傅里叶变换将所述规则干扰区域数据变换为频率空间域地震数据。
较佳的，根据以下公式确定频率空间域规则干扰滤波算子：
其中，XN(ω)为频率空间域规则干扰滤波算子，RN(ω)为由给定的视速度时差参数在频率空间域计算的方程组的系数，2N+1为频率空间域规则干扰滤波算子的道数。
较佳的，所述系数RN(ω)包括实部系数Rrm(ω)和虚部系数Rim(ω)，所述实部系数Rrm(ω)和虚部系数Rim(ω)分别通过以下公式确定：
Rrm(ω)＝amp*sin A*cos(B)/A；
Rim(ω)＝amp*sin A*sin(B)/A；
其中，A＝(Dt2-Dt1)/2，B＝(Dt2+Dt1)/2，amp为地震数据规则干扰波的振幅，Dt1为视速度最大平均时差，Dt2为视速度最小平均时差。
较佳的，所述将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子，包括：
通过傅里叶反变换将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子。
较佳的，所述根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据，包括：
在时间域内，将所述时间域规则干扰滤波算子与所述规则干扰区域数据进行褶积运算，获得压制后的地震数据。
较佳的，在获得压制后的地震数据之后，还包括：
将压制后的地震数据替换所述炮集记录中的规则干扰区域数据，从而组合成新的炮集记录。
另一方面，本申请实施例提供了一种规则干扰压制装置，包括：
炮集记录获取模块，用于获取非规则采样地震数据的炮集记录；
极限时差确定模块，用于确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；
干扰区域确定模块，用于根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；
第一数据转换模块，用于将所述规则干扰区域数据变换为频率空间域地震数据；
滤波算子确定模块，用于确定频率空间域规则干扰滤波算子；
第二数据转换模块，用于将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；
规则干扰压制模块，用于根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
较佳的，所述确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差，包括：
确定所述炮集记录中横跨特定数量地震道的规则干扰波的最大视速度和最小视速度；
将所述规则干扰波的最大视速度和最小视速度分别除以所述规则干扰波所横跨地震道的特定数量，对应获得所述炮集记录中相邻地震道的视速度最大平均时差和视速度最大平均时差。
较佳的，所述根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据，包括：
在所述炮集记录上识别出规则干扰和有效信号的界线；
判断所述界线是否大于所述视速度最大平均时差或小于所述视速度最小平均时差；
如果述界线大于所述视速度最大平均时差或小于所述视速度最小平均时差，则根据所述非规则采样地震数据在所述炮集记录上的时间分布和空间范围以及所述界线，将所述炮集记录划分为规则干扰区域数据和非规则干扰区域数据，并从所述炮集记录中提取出所述规则干扰区域数据。
较佳的，根据以下公式确定频率空间域规则干扰滤波算子：
其中，XN(ω)为频率空间域规则干扰滤波算子，RN(ω)为由给定的视速度时差参数在频率空间域计算的方程组的系数，2N+1为频率空间域规则干扰滤波算子的道数。
较佳的，所述系数RN(ω)包括实部系数Rrm(ω)和虚部系数Rim(ω)，所述实部系数Rrm(ω)和虚部系数Rim(ω)分别通过以下公式确定：
Rrm(ω)＝amp*sin A*cos(B)/A；
Rim(ω)＝amp*sin A*sin(B)/A；
其中，A＝(Dt2-Dt1)/2，B＝(Dt2+Dt1)/2，amp为地震数据规则干扰波的振幅，Dt1为视速度最大平均时差，Dt2为视速度最小平均时差。
较佳的，所述装置还包括：
炮集记录组合模块，用于将压制后的地震数据替换所述炮集记录中的规则干扰区域数据，从而组合成新的炮集记录。
另一方面，本申请实施例还提供了另一种规则干扰压制装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：
获取非规则采样地震数据的炮集记录；
确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；
根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；
将所述规则干扰区域数据变换为频率空间域地震数据；
确定频率空间域规则干扰滤波算子；
将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；
根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
由此可见，本申请实施例的规则干扰压制方案利用了有效信号与规则干扰在传播时间上和视速度的差异，在频率内完成了规则干扰滤波算子，然后再将频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子，从而避免了因空间上采样不均匀给频率空间域去噪带来的不利影响，因而提高了规则干扰噪声压制精度，且去噪后有效信号的保真度比较高，从而可为后期地震数据高保真处理奠定了坚实的基础。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本申请一实施方式中规则干扰压制方法的流程图；
图2为本申请一实施方式中规则干扰压制前的单炮地震记录；
图3为图2所示单炮地震记录采用本申请实施方式的规则干扰压制方法压制后得到的单炮地震记录；
图4为图2所示单炮地震记录采用现有技术中的FK域规则干扰压制方法压制后得到的单炮地震记录；
图5为图2所示单炮地震记录采用本申请实施方式的规则干扰压制方法压制后输出的噪声记录；
图6为图2所示单炮地震记录采用现有技术中的FK域规则干扰压制方法压制后输出的噪声记录；
图7为本申请另一实施方式中规则干扰压制装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
参考图1所示，本申请实施例的规则干扰压制方法可以包括以下步骤：
S101、获取非规则采样地震数据的炮集记录。
在本申请一实施方式中，可按照地震道地理位置关系排列地震道，以获取非规则采样地震数据中的炮集记录，这样可以尽量保持规则干扰的规则特性。
S102、确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差。
在本申请一实施方式中，所述确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差可以包括：
1)、确定所述炮集记录中横跨特定数量地震道的规则干扰波的最大视速度和最小视速度。在本申请一实施方式中，所述特定数量可根据需要或依据经验设定，例如可以为10。
2)、将所述规则干扰波的最大视速度和最小视速度分别除以所述规则干扰波所横跨地震道的特定数量，对应获得所述炮集记录中相邻地震道的视速度最大平均时差和视速度最大平均时差。
S103、根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据。
在本申请一实施方式中，在提取所述炮集记录中的规则干扰区域数据后，可仅对所述炮集记录中的规则干扰区域数据进行规则干扰压制，而不对所述炮集记录中的非规则干扰区域数据进行处理，从而可以更好的保护无规则干扰存在部分原有地震数据的波形特征。其中，所述规则干扰区域数据即为规则干扰的炮集记录部分，所述非规则干扰区域数据即为不含规则干扰的炮集记录部分。
在本申请一实施方式中，所述根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据可以包括：
1)、在所述炮集记录上识别出规则干扰和有效信号的界线。
2)、判断所述界线是否大于所述视速度最大平均时差或小于所述视速度最小平均时差。上述界线是一个大致界线，为了准确确定界线是否满足要求，需要判断所述界线是否大于所述视速度最大平均时差或小于所述视速度最小平均时差，即判断所述界线所在直线的视速度平均时差是否大于所述视速度最大平均时差或小于所述视速度最小平均时差。具体的，可以从所述界线上选择两点，然后在这个两个点对应地震数据道头中找到对应的炮间距大小，并在剖面上找到对应的地震反射时间，然后确定这两个点两点所在直线的视速度平均时差。
3)、如果述界线大于所述视速度最大平均时差或小于所述视速度最小平均时差，则根据所述非规则采样地震数据在所述炮集记录上的时间分布和空间范围以及所述界线，将所述炮集记录划分为规则干扰区域数据和非规则干扰区域数据，并从所述炮集记录中提取出所述规则干扰区域数据。
在本申请另一实施方式中，如果述界线小于所述视速度最大平均时差或大于所述视速度最小平均时差，则可以从所述界线附近沿所述界线延伸方向重新确定新的界线，并判断其是否满足要求。
S104、将所述规则干扰区域数据变换为频率空间域地震数据。
在本申请一实施方式中，可通过例如傅里叶变换等将所述规则干扰区域数据变换为频率空间域地震数据。在本申请另一实施方式中，在将所述规则干扰区域数据变换为频率空间域地震数据之前还可以对述规则干扰区域数据进行时间空间窗划分。例如如果时间上选择512个样点，空间上选择9个地震道，则在进行规则干扰压制时，512 个样点和9个地震道是后续一次所处理数据量的默认值，具体可以根据实验情况进行调整。相应的，可通过傅里叶变换将所述规则干扰区域数据中的9个地震道在时间方向上分别做一次傅里叶变换，从而将所述规则干扰区域数据中的9个地震道在时间方向上变换为频率空间域地震数据。
S105、确定频率空间域规则干扰滤波算子。
在本申请一实施方式中，可根据以下公式确定频率空间域规则干扰滤波算子：
其中，XN(ω)为频率空间域规则干扰滤波算子，RN(ω)为由给定的视速度时差参数在频率空间域计算的方程组的系数，2N+1为频率空间域规则干扰滤波算子的道数。
在本申请一实施方式中，上述系数RN(ω)还可以包括实部系数Rrm(ω)和虚部系数Rim(ω)。其中，所述实部系数Rrm(ω)和虚部系数Rim(ω)分别通过以下公式确定：
Rrm(ω)＝amp*sin A*cos(B)/A；
Rim(ω)＝amp*sin A*sin(B)/A；
其中，A＝(Dt2-Dt1)/2，B＝(Dt2+Dt1)/2，amp为地震数据规则干扰波的振幅，Dt1为视速度最大平均时差，Dt2为视速度最小平均时差。
S106、将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子。
在本申请一实施方式中，可通过例如傅里叶反变换等将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子。
S107、根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
在本申请一实施方式中，所述根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制是指，在时间域内，将所述时间域规则干扰滤波算子与所述规则干扰区域数据进行褶积运算，获得压制后的地震数据。
S108、将压制后的地震数据替换所述炮集记录中的规则干扰区域数据，从而组合成新的炮集记录。
在本申请一实施方式中，可按照上述步骤完成每个单炮记录的规则干扰压制，进而可以完成对全工区内所述地震数据的规则干扰压制。而在完成每个单炮记录的规则干扰压制后，还可以将压制后的地震数据替换所述炮集记录中的规则干扰区域数据，从而组合成一个没有规则干扰的新的单炮记录。
下面介绍本申请一示例性实施方方式。参考图2示出了非均匀采样规则干扰的单炮记录，其对应的地震道是按照地理位置进行了排序。采用本申请的规则干扰压制方法可得到如图3所示的单炮地震记录，并还可输出如图5所示的噪声记录。为了便于比较，还可以对于图2所示的单炮记录采用现有技术中的FK域规则干扰压制方法进行压制，从而可获得如图4所示的单炮地震记录，并还可输出如图6所示的噪声记录。
对比图2、图3和图4可以发现，这两种方法对规则干扰压制都比较彻底，但是对比图2、图3和图5所示，采用本申请的规则干扰压制方法得到的单炮记录上的振幅特征相对于图2变化不大，说明这种去噪方法可以比较好的保护地震数据的振幅特征。对比图2、图4和图6所示，而采用现有技术中的FK域规则干扰压制方法得到的单炮记录上的振幅特征相对于图2变化较大，噪声衰减也不够彻底。由此可见，本申请实施方式的规则干扰压制方法明显优于现有技术中的FK域规则干扰压制方法。因此，本申请实施方式的规则干扰压制方法不但能压制规则干扰，还能更好地保持振幅特征，从而能够达到在压制规则干扰的同时，提高地震数据信噪比和分辨率的目的。
虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
本申请一实施方式的规则干扰压制装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：
获取非规则采样地震数据的炮集记录；
确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；
根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；
将所述规则干扰区域数据变换为频率空间域地震数据；
确定频率空间域规则干扰滤波算子；
将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；
根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
参考图7所示，本申请另一实施方式的规则干扰压制装置可以包括：
炮集记录获取模块71，可以用于获取非规则采样地震数据的炮集记录；
极限时差确定模块72，可以用于确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；
干扰区域确定模块73，可以用于根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；
第一数据转换模块74，可以用于将所述规则干扰区域数据变换为频率空间域地震数据；
滤波算子确定模块75，可以用于确定频率空间域规则干扰滤波算子；
第二数据转换模块76，可以用于将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；
规则干扰压制模块77，可以用于根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
本申请另一实施方式，所述装置还可以包括：
炮集记录组合模块78，可以用于将压制后的地震数据替换所述炮集记录中的规则干扰区域数据，从而组合成新的炮集记录。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/ 或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (16)
1.一种规则干扰压制方法，其特征在于，包括：
获取非规则采样地震数据的炮集记录；
确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；
根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；
将所述规则干扰区域数据变换为频率空间域地震数据；
确定频率空间域规则干扰滤波算子；
将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；
根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
2.如权利要求1所述的规则干扰压制方法，其特征在于，所述确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差，包括：
确定所述炮集记录中横跨特定数量地震道的规则干扰波的最大视速度和最小视速度；
将所述规则干扰波的最大视速度和最小视速度分别除以所述规则干扰波所横跨地震道的特定数量，对应获得所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差。
3.如权利要求1所述的规则干扰压制方法，其特征在于，所述根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据，包括：
在所述炮集记录上识别出规则干扰和有效信号的界线；
判断所述界线是否大于所述视速度最大平均时差或小于所述视速度最小平均时差；
如果述界线大于所述视速度最大平均时差或小于所述视速度最小平均时差，则根据所述非规则采样地震数据在所述炮集记录上的时间分布和空间范围以及所述界线，将所述炮集记录划分为规则干扰区域数据和非规则干扰区域数据，并从所述炮集记录中提取出所述规则干扰区域数据。
4.如权利要求1所述的规则干扰压制方法，其特征在于，所述将所述规则干扰区域数据变换为频率空间域地震数据，包括：
通过傅里叶变换将所述规则干扰区域数据变换为频率空间域地震数据。
5.如权利要求1所述的规则干扰压制方法，其特征在于，根据以下公式确定频率空间域规则干扰滤波算子：
其中，XN(ω)为频率空间域规则干扰滤波算子，RN(ω)为由给定的视速度时差参数在频率空间域计算的方程组的系数，2N+1为频率空间域规则干扰滤波算子的道数。
6.如权利要求5所述的规则干扰压制方法，其特征在于，所述系数RN(ω)包括实部系数Rrm(ω)和虚部系数Rim(ω)，所述实部系数Rrm(ω)和虚部系数Rim(ω)分别通过以下公式确定：
Rrm(ω)＝amp*sinA*cos(B)/A；
Rim(ω)＝amp*sinA*sin(B)/A；
其中，A＝(Dt2-Dt1)/2，B＝(Dt2+Dt1)/2，amp为地震数据规则干扰波的振幅，Dt1为视速度最大平均时差，Dt2为视速度最小平均时差。
7.如权利要求1所述的规则干扰压制方法，其特征在于，所述将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子，包括：
通过傅里叶反变换将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子。
8.如权利要求1所述的规则干扰压制方法，其特征在于，所述根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据，包括：
在时间域内，将所述时间域规则干扰滤波算子与所述规则干扰区域数据进行褶积运算，获得压制后的地震数据。
9.如权利要求1所述的规则干扰压制方法，其特征在于，在获得压制后的地震数据之后，还包括：
将压制后的地震数据替换所述炮集记录中的规则干扰区域数据，从而组合成新的炮集记录。
10.一种规则干扰压制装置，其特征在于，包括：
炮集记录获取模块，用于获取非规则采样地震数据的炮集记录；
极限时差确定模块，用于确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；
干扰区域确定模块，用于根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；
第一数据转换模块，用于将所述规则干扰区域数据变换为频率空间域地震数据；
滤波算子确定模块，用于确定频率空间域规则干扰滤波算子；
第二数据转换模块，用于将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；
规则干扰压制模块，用于根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
11.如权利要求10所述的规则干扰压制装置，其特征在于，所述确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差，包括：
确定所述炮集记录中横跨特定数量地震道的规则干扰波的最大视速度和最小视速度；
将所述规则干扰波的最大视速度和最小视速度分别除以所述规则干扰波所横跨地震道的特定数量，对应获得所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差。
12.如权利要求10所述的规则干扰压制装置，其特征在于，所述根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据，包括：
在所述炮集记录上识别出规则干扰和有效信号的界线；
判断所述界线是否大于所述视速度最大平均时差或小于所述视速度最小平均时差；
如果述界线大于所述视速度最大平均时差或小于所述视速度最小平均时差，则根据所述非规则采样地震数据在所述炮集记录上的时间分布和空间范围以及所述界线，将所述炮集记录划分为规则干扰区域数据和非规则干扰区域数据，并从所述炮集记录中提取出所述规则干扰区域数据。
13.如权利要求10所述的规则干扰压制装置，其特征在于，根据以下公式确定频率空间域规则干扰滤波算子：
其中，XN(ω)为频率空间域规则干扰滤波算子，RN(ω)为由给定的视速度时差参数在频率空间域计算的方程组的系数，2N+1为频率空间域规则干扰滤波算子的道数。
14.如权利要求13所述的规则干扰压制装置，其特征在于，所述系数RN(ω)包括实部系数Rrm(ω)和虚部系数Rim(ω)，所述实部系数Rrm(ω)和虚部系数Rim(ω)分别通过以下公式确定：
Rrm(ω)＝amp*sinA*cos(B)/A；
Rim(ω)＝amp*sinA*sin(B)/A；
其中，A＝(Dt2-Dt1)/2，B＝(Dt2+Dt1)/2，amp为地震数据规则干扰波的振幅，Dt1为视速度最大平均时差，Dt2为视速度最小平均时差。
15.如权利要求10所述的规则干扰压制装置，其特征在于，所述装置还包括：
炮集记录组合模块，用于将压制后的地震数据替换所述炮集记录中的规则干扰区域数据，从而组合成新的炮集记录。
16.一种规则干扰压制装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，其特征在于，所述计算机程序被所述处理器运行时执行如下步骤：
获取非规则采样地震数据的炮集记录；
确定所述炮集记录中相邻地震道的视速度最大平均时差和视速度最小平均时差；
根据所述视速度最大平均时差和视速度最小平均时差提取所述炮集记录中的规则干扰区域数据；
将所述规则干扰区域数据变换为频率空间域地震数据；
确定频率空间域规则干扰滤波算子；
将所述频率空间域规则干扰滤波算子变换为时间域规则干扰滤波算子；
根据所述时间域规则干扰滤波算子对所述规则干扰区域数据进行滤波压制，获得压制后的地震数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810101694.6A CN108363103B (zh) | 2018-02-01 | 2018-02-01 | 规则干扰压制方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810101694.6A CN108363103B (zh) | 2018-02-01 | 2018-02-01 | 规则干扰压制方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108363103A CN108363103A (zh) | 2018-08-03 |
| CN108363103B true CN108363103B (zh) | 2019-10-11 |
Family
ID=63004377
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810101694.6A Active CN108363103B (zh) | 2018-02-01 | 2018-02-01 | 规则干扰压制方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108363103B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110954952B (zh) * | 2019-12-05 | 2021-06-18 | 山东省科学院激光研究所 | 一种微震信号初至波类型判别和波速校正方法 |
| CN113126164B (zh) * | 2020-01-10 | 2024-04-30 | 中国石油天然气集团有限公司 | 地震数据去噪方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102681014A (zh) * | 2012-05-23 | 2012-09-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于多项式拟合的规则线性干扰压制方法 |
| CN104849758A (zh) * | 2015-05-05 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 针对地震数据中的规则干扰的压制方法 |
| CN105242309A (zh) * | 2015-10-19 | 2016-01-13 | 中国石油天然气集团公司 | 一种地震数据中规则干扰的压制方法及装置 |
-
2018
- 2018-02-01 CN CN201810101694.6A patent/CN108363103B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102681014A (zh) * | 2012-05-23 | 2012-09-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于多项式拟合的规则线性干扰压制方法 |
| CN104849758A (zh) * | 2015-05-05 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 针对地震数据中的规则干扰的压制方法 |
| CN105242309A (zh) * | 2015-10-19 | 2016-01-13 | 中国石油天然气集团公司 | 一种地震数据中规则干扰的压制方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 《三维转换波自适应矢量去噪技术》;韩世春等;《地球物理学进展》;20141215;第29卷(第6期);2642-2647 * |
| 《单道地震勘探应用中的干扰因素分析及压制》;卢胜周等;《四川地质学报》;20160626;第36卷(第2期);310-314 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108363103A (zh) | 2018-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107153216B (zh) | 确定地震波场的坡印廷矢量的方法、装置以及计算机存储介质 | |
| Liu et al. | Seismic data interpolation beyond aliasing using regularized nonstationary autoregression | |
| Bergamo et al. | Retrieving 2D structures from surface-wave data by means of space-varying spatial windowing | |
| CN108549106B (zh) | 混叠噪声压制方法及装置 | |
| CN107272058B (zh) | 成像方法、成像装置以及计算机存储介质 | |
| CN105242309B (zh) | 一种地震数据中规则干扰的压制方法及装置 | |
| CN109471171B (zh) | 一种混叠地震数据分离的方法、装置及系统 | |
| Haney et al. | Causal instrument corrections for short‐period and broadband seismometers | |
| CN105425289B (zh) | 确定低频波阻抗的方法和装置 | |
| CN108363103B (zh) | 规则干扰压制方法及装置 | |
| CN105911585B (zh) | 一种地震记录规则干扰波的提取方法及装置 | |
| Curry | Interpolation with Fourier-radial adaptive thresholding | |
| CN108897036B (zh) | 一种地震数据处理方法及装置 | |
| CN108445532B (zh) | 一种深度域反偏移方法及装置 | |
| CN104181600B (zh) | 地震数据线性噪声衰减方法及装置 | |
| CN107561588B (zh) | 一种地震数据噪声压制方法及装置 | |
| CN107656313B (zh) | 一种拉东域多次波压制方法及装置 | |
| CN105068119B (zh) | 低频地震数据中的面波的衰减方法及其装置 | |
| CN109597964A (zh) | 理论油藏渗流参数值确定方法、系统、设备、可读介质 | |
| CN107831541B (zh) | 基于高密度vsp数据的薄地层识别方法及装置 | |
| CN109856672B (zh) | 基于深度波数谱的瞬变波包提取方法、存储介质与终端 | |
| CN108957553A (zh) | 动校正量递推修正的无拉伸畸变动校正方法及装置 | |
| CN102338884B (zh) | 物探中的椭圆窗方向带通保幅滤波数据处理方法 | |
| Liu et al. | Seismic random noise attenuation based on adaptive nonlocal median filter | |
| CN109917459A (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 |