CN104865601B - 确定地震勘探中采集的地震记录的类型的方法 - Google Patents
确定地震勘探中采集的地震记录的类型的方法 Download PDFInfo
- Publication number
- CN104865601B CN104865601B CN201510274159.7A CN201510274159A CN104865601B CN 104865601 B CN104865601 B CN 104865601B CN 201510274159 A CN201510274159 A CN 201510274159A CN 104865601 B CN104865601 B CN 104865601B
- Authority
- CN
- China
- Prior art keywords
- mrow
- seismic
- unit
- wave
- msub
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种确定地震勘探中采集的地震记录的类型的方法，所述方法包括：(A)获取地震记录，其中，地震记录为一个炮点对应的所有检波器在预定时间内按预定采样间隔接收到的地震数据，地震数据包括地震振幅；(B)根据所述地震记录中的地震数据，确定所述所有检波器中设置在一个排列上的全部检波器接收到地震波的初至时间；(C)判断接收到的地震波的初至时间是否满足预定条件；(D)当接收到的地震波的初至时间满足预定条件时，所述地震记录为单炮记录。根据本发明的示例性实施例的方法能够准确地识别地震记录的类型。
Description
技术领域
本发明属于地震勘探领域，更具体地说，涉及一种确定地震勘探中采集的地震记录的类型的方法。
背景技术
在地震勘探生产过程中，每进行几次单炮记录后就要进行一次背景记录(即，背景噪声的记录)。这样，获取的地震记录中就夹杂着不同类型的地震记录，一种是用来确定地下地质构造的单炮记录，另一种是用来分析背景噪声的背景记录。然而，随着地震数据采集的高效化，迫切需要对野外现场采集的地震记录进行分类以快速提取所需的地震记录。目前，现有的确定地震记录的类型的方法无法准确地识别地震记录的类型。
发明内容
本发明的示例性实施例的目的在于提供一种确定地震勘探中采集的地震记录的类型的方法，以准确地识别地震记录的类型。
本发明的示例性实施例提供一种确定地震勘探中采集的地震记录的类型的方法，所述方法包括：(A)获取地震记录，其中，地震记录为一个炮点对应的所有检波器在预定时间内按预定采样间隔接收到的地震数据，地震数据包括地震振幅；(B)根据所述地震记录中的地震数据，确定所述所有检波器中设置在一个排列上的全部检波器接收到地震波的初至时间；(C)判断接收到的地震波的初至时间是否满足预定条件；(D)当接收到的地震波的初至时间满足预定条件时，所述地震记录为单炮记录。
可选地，所述方法还包括：(E)当接收到的地震波的初至时间不满足预定条件时，所述地震记录为背景记录。
可选地，所述预定条件为所述设置在一个排列上的全部检波器中有连续的预定数量的检波器接收到的地震波的初至时间递增或递减。
可选地，确定一个检波器接收到地震波的初至时间的步骤包括：选取所述一个检波器接收到的所有地震振幅中对应的时间位于预定时窗内的多个地震振幅；通过将预定长度的滑动时窗在所述多个地震振幅上进行滑动，将所述多个地震振幅划分为多个单元，其中，在所述滑动时窗的各个位置上的滑动时窗内的地震振幅构成一个单元，滑动步长为预定时间间隔，所述预定长度为所述预定采样间隔的奇数倍；将所述多个单元按单元的起始时间从小到大的顺序排列，分别计算排列后的多个单元中相邻的两个单元之间的能量比值，其中，能量比值为后一单元的能量值比上前一单元的能量值；确定计算的所述能量比值中最大的能量比值对应的单元，所述单元中位于中间位置的地震振幅对应的时间为所述一个检波器接收到地震波的初至时间。
可选地，所述预定时间间隔为所述预定采样间隔。
可选地，通过如下公式计算所述能量比值：
其中，Ai+1为与第i+1单元对应的能量比值，xi+1(t)为第i+1单元内的时间为t时的地震振幅，Ti+1,1为第i+1单元的起始时间，Ti+1,2为第i+1单元的终止时间，xi(t)为第i单元内的时间为t时的地震振幅，Ti,1为第i单元的起始时间，Ti,2为第i单元的终止时间，i为大于0的整数，其中，第i单元与第i+1单元为相邻的两个单元，第i单元为前一单元，第i+1单元为后一单元。
可选地，在地震勘探为三维地震勘探时，步骤(B)中所述一个排列为所述所有检波器设置到的多个排列中距离所述一个炮点最近的一个检波器所在的一个排列。
根据本发明的示例性实施例的一种确定地震勘探中采集的地震记录的类型的方法，能够准确地识别地震记录的类型。
将在接下来的描述中部分阐述本发明另外的方面和/或优点，还有一部分通过描述将是清楚的，或者可以经过本发明的实施而得知。
附图说明
通过下面结合附图进行的对实施例的描述，本发明的上述和/或其它目的和优点将会变得更加清楚，其中：
图1是示出根据本发明示例性实施例的确定地震勘探中采集的地震记录的类型的方法的流程图。
具体实施方式
现将详细描述本发明的示例性实施例，所述实施例的示例在附图中示出，其中，相同的标号指示相同的部分。以下将通过参照附图来说明所述实施例，以便解释本发明。
图1是示出根据本发明示例性实施例的确定地震勘探中采集的地震记录的类型的方法的流程图。
如图1所示，在步骤S100，获取地震记录。这里，地震记录为一个炮点对应的所有检波器在预定时间内按预定采样间隔接收到的地震数据，地震数据包括地震振幅。换言之，地震记录就是在预定时间内使用对应一个炮点的所有检波器分别对地震波进行采样而得到的地震数据，采样间隔为所述预定采样间隔。所述地震波可以为人工产生的地震波(即在所述所有检波器对应的炮点产生的地震波)，也可以为由干扰源产生的地震波(即背景噪声)。
在步骤S200，根据步骤S100获取的地震记录中的地震数据，确定所述所有检波器中设置在一个排列上的全部检波器接收到地震波的初至时间。在二维地震勘探中，对应一个炮点的所有检波器都设置在一个排列中，则所述全部检波器为所述所有检波器。在三维地震勘探中，对应一个炮点的所有检波器设置在多个排列中，此时，所述全部检波器可以为所述所有检波器中的设置在多个排列中任一个排列上的检波器。优选地，在三维地震勘探时，步骤S200中所述一个排列为所述所有检波器设置到的多个排列中距离所述所有检波器对应的炮点最近的一个检波器所在的一个排列。由此，可提高后述步骤S300的执行效率，进而提高识别地震记录类型的效率。
可以理解，这里，可使用各种现有的确定初至时间的方法来确定所述全部检波器接收到地震波的初至时间。作为示例，确定一个检波器接收到地震波的初至时间的步骤包括：
首先，选取所述一个检波器接收到的所有地震振幅中对应的时间位于预定时窗内的多个地震振幅。即，从所述所有地震振幅中选取出采样时间在预定时窗内的多个地震振幅。这里，预定时窗位于所述预定时间内，可以根据经验来设定所述预定时窗的位置和所述预定时窗的长度。
然后，通过将预定长度的滑动时窗在所述多个地震振幅上进行滑动，将所述多个地震振幅划分为多个单元。这里，在所述滑动时窗的各个位置上的所述滑动时窗内的地震振幅构成一个单元，滑动步长为预定时间间隔，所述预定长度为所述预定采样间隔的奇数倍。可以理解，所述预定长度小于所述预定时窗的长度。优选地，所述预定时间间隔为所述预定采样间隔，即，所述滑动时窗每次滑动一个采样点(即，地震振幅)。所述预定时间间隔还可设定为所述预定采样间隔的大于1的整数倍。另外，所述滑动可以从所述多个地震振幅中采样时间最小的地震振幅开始滑动到采样时间最大的地震振幅，也可以与前述相反从所述多个地震振幅中采样时间最大的地震振幅开始滑动到采样时间最小的地震振幅。
其次，将所述多个单元按单元的起始时间从小到大的顺序排列，分别计算排列后的多个单元中相邻的两个单元之间的能量比值。这里，能量比值为后一单元的能量值比上前一单元的能量值。这里，所述单元的起始时间为所述单元所对应的位置上的滑动时窗的起始时间。这里，可通过现有的各种方法计算所述能量比值。作为示例，可通过如下公式计算所述能量比值：
这里，Ai+1为与第i+1单元对应的能量比值，xi+1(t)为第i+1单元内的时间(即，对应的时间)为t时的地震振幅，Ti+1,1为第i+1单元的起始时间，Ti+1,2为第i+1单元的终止时间，xi(t)为第i单元内的时间(即，对应的时间)为t时的地震振幅，Ti,1为第i单元的起始时间，Ti,2为第i单元的终止时间，i为大于0的整数。这里，第i单元与第i+1单元为相邻的两个单元，第i单元为前一单元，第i+1单元为后一单元。
最后，确定计算的所述能量比值中最大的能量比值对应的单元，所述单元中位于中间位置的地震振幅对应的时间为所述一个检波器接收到地震波的初至时间。这里，能量比值对应的单元为计算所述能量比值时的后一单元。
在步骤S300，判断步骤S200接收到的地震波的初至时间是否满足预定条件。这里，所述预定条件可以为任何能够确定所述地震记录的类型的条件。优选地，所述预定条件为所述设置在一个排列上的全部检波器中有连续的预定数量的检波器接收到的地震波的初至时间递增或递减。这里，所述预定数量可以根据试验来确定。举例说来，在所述预定数量设定为5时，则所述预定条件为所述一个排列中有5个连续的检波器所接收到地震波的初至时间递增或递减。
具体地，在步骤S300进行所述判断时，要依次比较所述一个排列上的全部检波器接收到的地震波的初至时间，当判断接收到的地震波的初至时间满足预定条件时，无论是否已经对所述全部检波器接收到的地震波的初至时间都进行比较，都可停止比较，并进入下一步。
在步骤S400，当步骤S200接收到的地震波的初至时间满足预定条件时，所述地震记录为单炮记录。
在步骤S500，当步骤S200接收到的地震波的初至时间不满足预定条件时，所述地震记录为背景记录。
根据本发明的示例性实施例的一种确定地震勘探中采集的地震记录的类型的方法，能够准确地识别地震记录的类型。
根据本发明的上述方法可以被实现为计算机可读记录介质中的计算机代码。本领域技术人员可以根据对上述方法的描述来实现所述计算机代码。当所述计算机代码在计算机中被执行时实现本发明的上述方法。
本发明的以上实施例仅仅是示例性的，而本发明并不受限于此。本领域技术人员应该理解：在不脱离本发明的原理和精神的情况下，可对这些实施例进行改变，其中，本发明的范围在权利要求及其等同物中限定。
Claims (6)
1.一种确定地震勘探中采集的地震记录的类型的方法，所述方法包括：
(A)获取地震记录，其中，地震记录为一个炮点对应的所有检波器在预定时间内按预定采样间隔接收到的地震数据，地震数据包括地震振幅，所述地震记录分为背景记录和单炮记录两种类型；
(B)根据所述地震记录中的地震数据，确定所述所有检波器中设置在一个排列上的全部检波器接收到地震波的初至时间；
(C)判断接收到的地震波的初至时间是否满足预定条件，其中，预定条件为所述设置在一个排列上的全部检波器中有连续的预定数量的检波器接收到的地震波的初至时间递增或递减；
(D)当接收到的地震波的初至时间满足预定条件时，所述地震记录为单炮记录。
2.根据权利要求1所述的方法，其特征在于，所述方法还包括：
(E)当接收到的地震波的初至时间不满足预定条件时，所述地震记录为背景记录。
3.根据权利要求1所述的方法，其特征在于，确定一个检波器接收到地震波的初至时间的步骤包括：
选取所述一个检波器接收到的所有地震振幅中对应的时间位于预定时窗内的多个地震振幅；
通过将预定长度的滑动时窗在所述多个地震振幅上进行滑动，将所述多个地震振幅划分为多个单元，其中，在所述滑动时窗的各个位置上的滑动时窗内的地震振幅构成一个单元，滑动步长为预定时间间隔，所述预定长度为所述预定采样间隔的奇数倍；
将所述多个单元按单元的起始时间从小到大的顺序排列，分别计算排列后的多个单元中相邻的两个单元之间的能量比值，其中，能量比值为后一单元的能量值比上前一单元的能量值；
确定计算的所述能量比值中最大的能量比值对应的单元，所述单元中位于中间位置的地震振幅对应的时间为所述一个检波器接收到地震波的初至时间。
4.根据权利要求3所述的方法，其特征在于，所述预定时间间隔为所述预定采样间隔。
5.根据权利要求3所述的方法，其特征在于，通过如下公式计算所述能量比值：
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>{</mo>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</munderover>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</munderover>
<msup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
其中，Ai+1为与第i+1单元对应的能量比值，xi+1(t)为第i+1单元内的时间为t时的地震振幅，Ti+1,1为第i+1单元的起始时间，Ti+1,2为第i+1单元的终止时间，xi(t)为第i单元内的时间为t时的地震振幅，Ti,1为第i单元的起始时间，Ti,2为第i单元的终止时间，i为大于0的整数，
其中，第i单元与第i+1单元为相邻的两个单元，第i单元为前一单元，第i+1单元为后一单元。
6.根据权利要求1所述的方法，其特征在于，在地震勘探为三维地震勘探时，步骤(B)中所述一个排列为所述所有检波器设置到的多个排列中距离所述一个炮点最近的一个检波器所在的一个排列。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510274159.7A CN104865601B (zh) | 2015-05-18 | 2015-05-18 | 确定地震勘探中采集的地震记录的类型的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510274159.7A CN104865601B (zh) | 2015-05-18 | 2015-05-18 | 确定地震勘探中采集的地震记录的类型的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104865601A CN104865601A (zh) | 2015-08-26 |
| CN104865601B true CN104865601B (zh) | 2017-08-25 |
Family
ID=53911559
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510274159.7A Active CN104865601B (zh) | 2015-05-18 | 2015-05-18 | 确定地震勘探中采集的地震记录的类型的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104865601B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108897040B (zh) * | 2018-05-10 | 2020-05-08 | 中国石油天然气集团有限公司 | 一种确定陆检与海底之间的耦合性的方法及装置 |
| CN109101615B (zh) * | 2018-08-07 | 2021-03-30 | 中国石油天然气集团有限公司 | 地震勘探数据处理方法及装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101609160A (zh) * | 2008-06-19 | 2009-12-23 | 中国石油天然气集团公司 | 一种抗背景噪音强的地震数据初至波自动拾取方法 |
| EP2594964A2 (en) * | 2011-11-21 | 2013-05-22 | CGGVeritas Services S.A. | Device and method for computing depth velocity variations |
| CN103439738A (zh) * | 2013-08-27 | 2013-12-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地震勘探单炮记录异常道识别方法 |
| CN103616722A (zh) * | 2013-11-28 | 2014-03-05 | 中国石油天然气股份有限公司 | 一种初至波旅行时拾取方法及装置 |
-
2015
- 2015-05-18 CN CN201510274159.7A patent/CN104865601B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101609160A (zh) * | 2008-06-19 | 2009-12-23 | 中国石油天然气集团公司 | 一种抗背景噪音强的地震数据初至波自动拾取方法 |
| EP2594964A2 (en) * | 2011-11-21 | 2013-05-22 | CGGVeritas Services S.A. | Device and method for computing depth velocity variations |
| CN103439738A (zh) * | 2013-08-27 | 2013-12-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地震勘探单炮记录异常道识别方法 |
| CN103616722A (zh) * | 2013-11-28 | 2014-03-05 | 中国石油天然气股份有限公司 | 一种初至波旅行时拾取方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 利用能量比法拾取地震初至的一种改进方法;左国平 等;《石油物探》;20040731;第43卷(第4期);第345-347页 * |
| 地震资料品质自动评价系统;段云卿;《勘探地球物理进展》;20060630;第29卷(第3期);第220-224页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104865601A (zh) | 2015-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108805269B (zh) | 一种基于lstm循环神经网络拾取震相到时的方法 | |
| Verdon et al. | Microseismic monitoring using a fiber-optic distributed acoustic sensor array | |
| CN110133715B (zh) | 一种基于初至时差和波形叠加的微地震震源定位方法 | |
| Yang et al. | Simultaneous earthquake detection on multiple stations via a convolutional neural network | |
| CN103439738B (zh) | 地震勘探单炮记录异常道识别方法 | |
| CN104155691B (zh) | 转换波各向异性速度分析方法及装置 | |
| CN107728210B (zh) | 多仪器采集数据中缺失道的确定方法和装置 | |
| CN100429527C (zh) | 三维地震资料处理质量监控技术 | |
| CN105425292A (zh) | 一种油气预测方法及装置 | |
| CN105093299A (zh) | 一种基于炮检距向量片技术优化观测系统的方法及装置 | |
| CN107402405B (zh) | 静相位虚源道集构建方法 | |
| CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
| CN105467460A (zh) | 电磁勘探方法和装置 | |
| Maity et al. | Neuro-evolutionary event detection technique for downhole microseismic surveys | |
| CN105005075B (zh) | 基于地震频率信息的多波匹配方法 | |
| CN104865601B (zh) | 确定地震勘探中采集的地震记录的类型的方法 | |
| CN105986819B (zh) | 用于测井资料自动处理与综合解释的方法和装置 | |
| CN106324666A (zh) | 横向各向同性介质转换波角道集抽取方法和装置 | |
| CN105319587A (zh) | 一种基于初至子波提取的可控震源整形滤波方法 | |
| CN103513279B (zh) | 一种基于地震波波动方程的照明分析计算方法及计算装置 | |
| CN106291733B (zh) | 一种便携式石油检测器 | |
| CN102288995A (zh) | 地震勘探采集质量量化分析评价系统 | |
| CN103201649A (zh) | 用于脉冲中子俘获西格马反演的方法和系统 | |
| CN103777244A (zh) | 一种地震裂缝属性体的定量分析方法 | |
| CN105445788B (zh) | 一种基于模型和全局寻优的速度谱自动解释方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180207 Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp. Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology Patentee before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200921 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: BGP Inc., China National Petroleum Corp. Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee before: BGP Inc., China National Petroleum Corp. |
|
| TR01 | Transfer of patent right |