CN109254323A - 基于能量包络的地震波初至拾取方法及计算机可读存储介质 - Google Patents
基于能量包络的地震波初至拾取方法及计算机可读存储介质 Download PDFInfo
- Publication number
- CN109254323A CN109254323A CN201710576613.3A CN201710576613A CN109254323A CN 109254323 A CN109254323 A CN 109254323A CN 201710576613 A CN201710576613 A CN 201710576613A CN 109254323 A CN109254323 A CN 109254323A
- Authority
- CN
- China
- Prior art keywords
- seismic
- breaks
- data
- channel data
- energy
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000007704 transition Effects 0.000 claims description 32
- 238000004590 computer program Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 description 4
- 241001269238 Data Species 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003908 quality control method Methods 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/362—Effecting static or dynamic corrections; Stacking
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
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
技术领域
本发明涉及地球物理勘探技术领域，更具体地，涉及一种基于能量包络的地震波初至拾取方法，其可应用于油气地震勘探及非常规地震勘探的资料处理中。
背景技术
地震勘探资料中存在大量复杂地表的勘探资料，在对地震勘探资料进行处理时，其静校正问题十分复杂，层析静校正等方法的应用效果都依赖于初至拾取的质量。当前初至拾取的方法一般是采用软件自动拾取，然后通过人工予以修正。但由于地表条件复杂，初至波受多种因素影响，自动拾取的质量往往不是很高，很大程度上需要人工拾取修正，工作量巨大，处理周期延长。为此，期待一种能够提高地震波初至拾取的质量和效率的地震波初至拾取方法。
发明内容
本发明的目的是提出一种基于能量包络的地震波初至拾取方法，其可克服现有的初至拾取方法工作量巨大、处理周期长的缺点。
本发明提出了一种基于能量包络的地震波初至拾取方法，包括：
对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；
去除所述首次拾取的地震波初至中的异常值，缩小所述首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；
基于所述原始地震资料计算振幅能量包络的地震道数据；
在所述二次初至拾取时窗内，对所述振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。
优选地，所述基于能量包络的地震波初至拾取方法还包括：
输入所述原始地震资料和空间属性定义，对所述原始地震资料进行滤波和增益处理。
优选地，所述基于能量包络的地震波初至拾取方法还包括：
对所述最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。
优选地，所述基于原始地震资料计算振幅能量包络的地震道数据包括：
基于所述原始地震资料获得地震道数据，对所述地震道数据的数值做平方运算，得到当前地震道的能量；
将所述地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；
将所述当前地震道的能量与所述相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。
优选地，通过炮相关方法对所述原始地震资料进行首次初至拾取。
本发明还提出一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：
对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；
去除所述首次拾取的地震波初至中的异常值，缩小所述首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；
基于所述原始地震资料计算振幅能量包络的地震道数据；
在所述二次初至拾取时窗内，对所述振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。
优选地，所述程序被处理器执行时还实现以下步骤：
输入所述原始地震资料和空间属性定义，对所述原始地震资料进行滤波和增益处理。
优选地，所述程序被处理器执行时还实现以下步骤：
对所述最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。
优选地，所述基于原始地震资料计算振幅能量包络的地震道数据包括：
基于所述原始地震资料获得地震道数据，对所述地震道数据的数值做平方运算，得到当前地震道的能量；
将所述地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；
将所述当前地震道的能量与所述相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。
优选地，通过炮相关方法对所述原始地震资料进行首次初至拾取。
本发明的有益效果在于：振幅能量包络的地震道数据的抗干扰能力增强，易于进行初至识别，因此根据本发明的基于能量包络的地震波初至拾取方法能够提高地震波初至拾取的效率和准确率，为提高后续的静校正精度奠定了良好的基础。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的附图标记通常代表相同部件。
图1示出了根据本发明实施例的基于能量包络的地震波初至拾取方法的流程图。
图2示出了根据本发明实施例的地震波初至拾取方法中的地震道数据。
图3出了根据本发明实施例的地震波初至拾取方法中的地震道的能量。
图4示出了根据本发明实施例的地震波初至拾取方法中将地震道数据的相位转换90°之后的地震道数据。
图5出了根据本发明实施例的地震波初至拾取方法中的相位转换后的地震道能量。
图6示出了根据本发明实施例的地震波初至拾取方法中的振幅能量包络的地震道数据。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明实施例的地震波初至拾取方法的流程图，如图1所示，根据本发明实施例的基于能量包络的地震波初至拾取方法包括以下步骤：
步骤1：对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至。
可使用炮相关等方法对原始地震资料进行快速的初至拾取，得到首次拾取的地震波初至，首次拾取的地震波初至可为后续的二次初至拾取提供拾取时窗依据；
步骤2：去除首次拾取的地震波初至中的异常值，并缩小首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗。
针对首次拾取的地震波初至进行质量控制，去除异常值，去除异常值后，地震波初至的范围有一定程度的缩小，缩小后的范围即可作为二次初至拾取时窗，在该范围内，不包含异常值。特别地，针对相邻两个地震道的初至值发生异常变化(例如相邻两个地震道的初至值的差异超过预定值)的不合理情况，需要去除异常值。
步骤3：基于原始地震资料计算振幅能量包络的地震道数据。
可以按照以下步骤基于原始地震资料计算振幅能量包络的地震道数据：
子步骤301：基于原始地震资料获得地震道数据，对地震道数据的数值做平方运算，得到当前地震道的能量；
子步骤302：将地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；
子步骤303：将当前地震道的能量与相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。
振幅能量包络的地震道数据的抗干扰能力增强，易于进行初至识别。
步骤4：在二次初至拾取时窗内，对振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。
在步骤2得到的二次初至拾取时窗内，对振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。二次初至拾取可通过现有软件来实现。
在一个示例中，根据本发明实施例的基于能量包络的地震波初至拾取方法还包括以下步骤：
输入原始地震资料和空间属性定义，对原始地震资料进行滤波和增益处理。对原始地震资料的预先处理时在不影响初至形态的前提下进行的。空间属性定义是指在地震资料处理中地震资料与空间位置的对应关系，将地震资料与空间位置相对应，确定二者之间的对应关系，是地震资料处理中必不可少的环节，是本领域技术人员根据现有技术能够实现的。
在一个示例中，根据本发明实施例的基于能量包络的地震波初至拾取方法还包括以下步骤：
对最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。
本发明实施例还提出一种计算机可读存储介质，其上存储有计算机程序，其中，程序被处理器执行时实现以下步骤：
对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；
去除首次拾取的地震波初至中的异常值，缩小首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；
基于原始地震资料计算振幅能量包络的地震道数据；
在二次初至拾取时窗内，对振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。
在一个示例中，程序被处理器执行时还实现以下步骤：
输入原始地震资料和空间属性定义，对原始地震资料进行滤波和增益处理。
在一个示例中，程序被处理器执行时还实现以下步骤：
对最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。
在一个示例中，基于原始地震资料计算振幅能量包络的地震道数据包括：
基于原始地震资料获得地震道数据，对地震道数据的数值做平方运算，得到当前地震道的能量；
将地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；
将当前地震道的能量与相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。
在一个示例中，通过炮相关方法对原始地震资料进行首次初至拾取。
实施例
在一个实施例中，基于能量包络的地震波初至拾取方法包括以下步骤：
步骤1：对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；
步骤2：去除首次拾取的地震波初至中的异常值，缩小首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；
步骤3：基于原始地震资料计算振幅能量包络的地震道数据，该步骤包括以下三个子步骤：
子步骤301：基于原始地震资料获得地震道数据，对地震道数据的数值做平方运算，得到当前地震道的能量，图2显示了地震道数据，图3显示了地震道的能量；
子步骤302：将地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量，图4显示了将地震道数据的相位转换90°之后的地震道数据，图5显示了相位转换后的地震道能量；
子步骤303：将当前地震道的能量与相位转换后的地震道能量相加，得到振幅能量包络的地震道数据，图6显示了振幅能量包络的地震道数据。
由于振幅能量包络的地震道数据的抗干扰能力增强，易于进行初至识别，因此根据本发明的基于能量包络的地震波初至拾取方法能够提高地震波初至拾取的效率和准确率，为提高后续的静校正精度奠定了良好的基础。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种基于能量包络的地震波初至拾取方法，包括：
对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；
去除所述首次拾取的地震波初至中的异常值，缩小所述首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；
基于所述原始地震资料计算振幅能量包络的地震道数据；
在所述二次初至拾取时窗内，对所述振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。
2.根据权利要求1所述的基于能量包络的地震波初至拾取方法，还包括：
输入所述原始地震资料和空间属性定义，对所述原始地震资料进行滤波和增益处理。
3.根据权利要求1所述的基于能量包络的地震波初至拾取方法，还包括：
对所述最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。
4.根据权利要求1所述的基于能量包络的地震波初至拾取方法，其中，所述基于原始地震资料计算振幅能量包络的地震道数据包括：
基于所述原始地震资料获得地震道数据，对所述地震道数据的数值做平方运算，得到当前地震道的能量；
将所述地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；
将所述当前地震道的能量与所述相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。
5.根据权利要求1所述的基于能量包络的地震波初至拾取方法，其中，通过炮相关方法对所述原始地震资料进行首次初至拾取。
6.一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：
对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；
去除所述首次拾取的地震波初至中的异常值，缩小所述首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；
基于所述原始地震资料计算振幅能量包络的地震道数据；
在所述二次初至拾取时窗内，对所述振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。
7.根据权利要求6所述的计算机可读存储介质，其中，所述程序被处理器执行时还实现以下步骤：
输入所述原始地震资料和空间属性定义，对所述原始地震资料进行滤波和增益处理。
8.根据权利要求6所述的计算机可读存储介质，其中，所述程序被处理器执行时还实现以下步骤：
对所述最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。
9.根据权利要求6所述的计算机可读存储介质，其中，所述基于原始地震资料计算振幅能量包络的地震道数据包括：
基于所述原始地震资料获得地震道数据，对所述地震道数据的数值做平方运算，得到当前地震道的能量；
将所述地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；
将所述当前地震道的能量与所述相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。
10.根据权利要求6所述的计算机可读存储介质，其中，通过炮相关方法对所述原始地震资料进行首次初至拾取。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710576613.3A CN109254323A (zh) | 2017-07-14 | 2017-07-14 | 基于能量包络的地震波初至拾取方法及计算机可读存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710576613.3A CN109254323A (zh) | 2017-07-14 | 2017-07-14 | 基于能量包络的地震波初至拾取方法及计算机可读存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109254323A true CN109254323A (zh) | 2019-01-22 |
Family
ID=65051885
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710576613.3A Pending CN109254323A (zh) | 2017-07-14 | 2017-07-14 | 基于能量包络的地震波初至拾取方法及计算机可读存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109254323A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113806673A (zh) * | 2020-06-16 | 2021-12-17 | 中国石油化工股份有限公司 | 初至质控方法、装置、计算机设备和存储介质 |
| US11467302B1 (en) | 2021-05-26 | 2022-10-11 | Saudi Arabian Oil Company | Seismic first breaks onset times determination by interface tracking using level-sets method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140219054A1 (en) * | 2013-02-04 | 2014-08-07 | King Fahd University Of Petroleum And Minerals | Method of first arrival picking of seismic refraction data |
| CN105204065A (zh) * | 2015-10-09 | 2015-12-30 | 中国石油天然气集团公司 | 一种拾取初至波的方法和装置 |
| CN105487114A (zh) * | 2015-12-08 | 2016-04-13 | 中南大学 | 一种微震信号p波初至点综合拾取方法 |
-
2017
- 2017-07-14 CN CN201710576613.3A patent/CN109254323A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140219054A1 (en) * | 2013-02-04 | 2014-08-07 | King Fahd University Of Petroleum And Minerals | Method of first arrival picking of seismic refraction data |
| CN105204065A (zh) * | 2015-10-09 | 2015-12-30 | 中国石油天然气集团公司 | 一种拾取初至波的方法和装置 |
| CN105487114A (zh) * | 2015-12-08 | 2016-04-13 | 中南大学 | 一种微震信号p波初至点综合拾取方法 |
Non-Patent Citations (5)
| Title |
|---|
| 付小宁等: "《数据处理及模型化方法》", 30 June 2016, 西北农林科技大学出版社 * |
| 大港油田科技丛书编委会: "《地震勘探资料处理和解释技术》", 30 September 1999, 北京：石油工业出版社 * |
| 张伟等: "地震道瞬时强度比法拾取初至波", 《地球物理学进展》 * |
| 张玉芬: "《反射波地震勘探原理和资料解释》", 31 July 2007, 北京：地质出版社 * |
| 徐钰等: "浅层折射波勘探中初至自动拾取新算法", 《石油地球物理勘探》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113806673A (zh) * | 2020-06-16 | 2021-12-17 | 中国石油化工股份有限公司 | 初至质控方法、装置、计算机设备和存储介质 |
| US11467302B1 (en) | 2021-05-26 | 2022-10-11 | Saudi Arabian Oil Company | Seismic first breaks onset times determination by interface tracking using level-sets method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105204065B (zh) | 一种拾取初至波的方法和装置 | |
| CN109738946B (zh) | 两步法剔除异常地震初至的方法及系统 | |
| CN109254323A (zh) | 基于能量包络的地震波初至拾取方法及计算机可读存储介质 | |
| CN105093319A (zh) | 基于三维地震数据的地面微地震静校正方法 | |
| CN102338885B (zh) | 三分量vsp资料初至时间自动拾取方法 | |
| CN117031539A (zh) | 一种自监督深度学习地震数据低频重建方法及系统 | |
| CN104597486B (zh) | 利用初至时间快速检测炮集数据的方法 | |
| CN106324671A (zh) | 一种检查炮点偏移的方法及装置 | |
| CN102478666B (zh) | 一种提高可控震源地震记录初至拾取精度的方法 | |
| CN100349005C (zh) | 地震数据多炮组合叠前深度偏移方法 | |
| CN115456040A (zh) | 一种基于卷积神经网络的p波捡拾算法 | |
| CN113627051B (zh) | 一种重力异常场分离方法、系统、存储介质和电子设备 | |
| CN104849684A (zh) | 一种示波器、校正装置及其自动校正水平中心的方法 | |
| CN116009077A (zh) | 一种基于谱比法的近地表q值建模方法、设备及介质 | |
| CN110568499B (zh) | 一种vsp地震资料的初至波时差校正方法及装置 | |
| CN104280766B (zh) | 一种利用局部数据同相轴斜率的直接偏移方法 | |
| CN109188526B (zh) | 一种改进型地表一致性剩余静校正方法 | |
| CN113721296A (zh) | 远震数据处理方法及装置 | |
| CN107895102B (zh) | 一种生长指令检查方法及装置 | |
| CN112379429A (zh) | 地震数据的振幅补偿方法及装置 | |
| US20230041249A1 (en) | Method and apparatus for extracting downgoing wavelet and attenuation parameters by using vertical seismic data | |
| CN105510974B (zh) | 折射波静校正方法和装置 | |
| CN109765614A (zh) | 一种地震前兆观测数据异常识别方法 | |
| CN108196306A (zh) | 全自动层析静校正 | |
| CN109387871A (zh) | P-sv转换波组合静校正方法及系统 |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190122 |
|
| RJ01 | Rejection of invention patent application after publication |