CN112444879B - 地震勘探数据的处理方法及存储介质 - Google Patents

地震勘探数据的处理方法及存储介质 Download PDF

Info

Publication number
CN112444879B
CN112444879B CN201910817078.5A CN201910817078A CN112444879B CN 112444879 B CN112444879 B CN 112444879B CN 201910817078 A CN201910817078 A CN 201910817078A CN 112444879 B CN112444879 B CN 112444879B
Authority
CN
China
Prior art keywords
shot
gather data
point
correction value
wavelength
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
Application number
CN201910817078.5A
Other languages
English (en)
Other versions
CN112444879A (zh
Inventor
余青露
邹少峰
陈林谦
马江林
刘晗
李进
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201910817078.5A priority Critical patent/CN112444879B/zh
Publication of CN112444879A publication Critical patent/CN112444879A/zh
Application granted granted Critical
Publication of CN112444879B publication Critical patent/CN112444879B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • 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

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

地震勘探数据的处理方法及存储介质
技术领域
本公开涉及地震勘探数据处理技术领域，特别地涉及一种地震勘探数据的处理方法及存储介质。
背景技术
由于黄土塬地区受表层结构特征(高程、风化层的速度、厚度和潜水面埋深)在区域上的不均匀性等诸多因素的影响，除了存在严重的短波长静校正问题外，还存在一定的中、长波长静校正问题，静校正问题严重且极其复杂是黄土塬地区资料的主要特点，它严重影响到均方根速度的准确求取和构造的落实。
发明内容
本公开提供一种地震勘探数据的处理方法及存储介质，以解决相关技术中地震资料中的中波长静校正量难以精确提取的技术问题。
为实现上述目的，本公开实施例的第一方面，提供一种地震勘探数据的处理方法，所述方法包括：
步骤S11：在获得长波长静校正处理后的地震勘探数据后，从所述地震勘探数据中抽取共检波点道集数据；
步骤S12：处理所述共检波点道集数据，以获得检波点的中波长静校正量并应用于所述地震勘探数据；
步骤S13：从经过所述检波点的中波长静校正量应用后的所述地震勘探数据中抽取炮点道集数据；
步骤S14：处理所述炮点道集数据，以获得炮点的中波长静校正量并应用于所述地震勘探数据；
步骤S15：从经过所述炮点的中波长静校正量应用后的所述地震勘探数据中抽取共检波点道集数据，以通过处理所述共检波点道集数据获得检波点的中波长静校正量；
步骤S16：判断所述炮点的中波长静校正量和步骤S15中的所述检波点的中波长静校正量是否满足预设条件；当所述炮点的中波长静校正量和所述检波点的中波长静校正量不满足预设条件时，迭代步骤S13至步骤S15，直至所述炮点的中波长静校正量和所述检波点的中波长静校正量满足预设条件；
步骤S17：当所述炮点的中波长静校正量和所述检波点的中波长静校正量满足预设条件时，输出满足预设条件的炮点的中波长静校正量和检波点的中波长静校正量并应用。
可选地，处理所述共检波点道集数据，以获得检波点的中波长静校正量，包括：
对所述共检波点道集数据进行线性动校正；
对线性动校正后的所述共检波点道集数据进行叠加，以获得叠加后的共检波点道集数据；
从叠加后的所述共检波点道集数据中提取检波点的中波长静校正量。
可选地，所述对所述共检波点道集数据进行线性动校正，包括：
利用速度扫描方式，根据初至波校平的情况，选取线性动校正的速度以对所述共检波点道集数据进行线性动校正。
可选地，对线性动校正后的所述共检波点道集数据进行叠加，以获得叠加后的共检波点道集数据，包括：
避开近炮检距的直达波部分，保留稳定的折射波部分，以对线性动校正后的所述共检波点道集数据进行叠加。
可选地，从叠加后的所述共检波点道集数据中提取检波点的中波长静校正量，包括：
基于所述共检波点道集数据叠加后的初至波趋势，拾取用于提取中波长剩余静校正的零时间线和延迟时间线；
将拾取的延迟时间线减去零时间线，以获得检波点的中波长静校正量。
可选地，处理所述炮点道集数据，以获得炮点的中波长静校正量，包括：
对所述炮点道集数据进行线性动校正；
对线性动校正后的所述炮点道集数据进行叠加，以获得叠加后的炮点道集数据；
从叠加后的所述炮点道集数据中提取炮点的中波长静校正量。
可选地，对所述炮点道集数据进行线性动校正，包括：
利用速度扫描方式，根据初至波校平的情况，选取线性动校正的速度以对所述炮点道集数据进行线性动校正。
可选地，对线性动校正后的所述炮点道集数据进行叠加，以获得叠加后的炮点道集数据，包括：
避开近炮检距的直达波部分，保留稳定的折射波部分，以对线性动校正后的所述炮点道集数据进行叠加。
可选地，从叠加后的所述炮点道集数据中提取炮点的中波长静校正量，包括：
基于所述炮点道集数据叠加后的初至波趋势，拾取用于提取中波长剩余静校正的零时间线和延迟时间线；
将拾取的延迟时间线减去零时间线，以获得炮点的中波长静校正量。
本公开实施例的第二方面，提供一种地震勘探数据的处理装置，所述装置包括：
存储器，其上存储有计算机程序；以及
处理器，用于执行所述存储器中的所述计算机程序，以实现上述第一方面中任一项所述方法的步骤。
本公开实施例的第三方面，提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述第一方面中任一项所述方法的步骤。
采用上述技术方案，至少能够达到如下技术效果：
针对复杂的地震勘探数据的处理，本公开通过采取交替迭代提取中波长静校正量的方法，可提高拾取的精度，进而提高处理成果构造的可靠性和资料的信噪比，为解释提供可靠的基础资料，使得本公开在黄土塬等静校正复杂探区地震数据处理中具有很好的应用前景。
本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本公开的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本公开，但并不构成对本公开的限制。在附图中：
图1是本公开一示例性实施例示出的一种地震勘探数据的处理方法流程图。
图2是本公开一示例性实施例示出的一种地震勘探数据的处理方法包括的步骤处理共检波点道集数据的流程图。
图3是本公开一示例性实施例示出的一种地震勘探数据的处理方法包括的步骤处理炮点道集数据的流程图。
图4是本公开一示例性实施例示出的共检波点道集数据的示意图。
图5是本公开一示例性实施例示出的线性动校正后的共检波点道集数据的示意图。
图6是本公开一示例性实施例示出的叠加后的共检波点道集数据的示意图。
图7是本公开一示例性实施例示出的检波点中波长静校正量应用后叠加的共检波点道集数据示意图。
图8是本公开一示例性实施例示出的线性动校正后的炮点道集数据示意图。
图9是本公开一示例性实施例示出的叠加后的炮点道集数据示意图。
图10是本公开一示例性实施例示出的炮点中波长静校正量应用后叠加的共检波点道集数据示意图。
图11是本公开一示例性实施例示出的另一共检波点道集数据示意图。
图12是本公开一示例性实施例示出的只有长波长静校正量处理后的单炮记录示意图。
图13是本公开一示例性实施例示出的炮点和检波点的中波长静校正量处理后的单炮记录示意图。
图14是本公开一示例性实施例示出的没有经过炮点和检波点的中波长静校正量应用后的叠加剖面示意图。
图15是本公开一示例性实施例示出的经过炮点和检波点的中波长静校正量应用后的叠加剖面示意图。
图16是本公开一示例性实施例示出的一种地震勘探数据的处理装置框图。
具体实施方式
以下将结合附图及实施例来详细说明本公开的实施方式，借此对本公开如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本公开的保护范围之内。
发明人经研究发现，如何确保静校正的真实可靠应是黄土塬地区资料处理中的最大难点。当前黄土塬地震资料处理中，静校正处理技术还不是很成熟，在长波长静校正问题后，还存在较大的中、短波长的静校正问题，尤其是中波长静校正量的精确提取问题仍然有待研究解决，本公开针对性的提出了一种地震勘探数据的处理方法，该方法提取的静校正量精度高，可较好解决中波长静校正量准确提取难题，为进一步短波长的剩余静校正处理奠定良好的基础。
实施例一
图1是本公开一示例性实施例示出的一种地震勘探数据的处理方法流程图，以解决相关技术中地震资料中的中波长静校正量难以精确提取的技术问题。如图1所示，该地震资料中的中波长静校正量难以精确提取方法可以包括以下步骤：
步骤S11，在获得长波长静校正处理后的地震勘探数据后，从所述地震勘探数据中抽取共检波点道集数据。
步骤S12，处理所述共检波点道集数据，以获得检波点的中波长静校正量并应用于所述地震勘探数据。
步骤S13，从经过所述检波点的中波长静校正量应用后的所述地震勘探数据中抽取炮点道集数据。
步骤S14，处理所述炮点道集数据，以获得炮点的中波长静校正量并应用于所述地震勘探数据。
步骤S15，从经过所述炮点的中波长静校正量应用后的所述地震勘探数据中抽取共检波点道集数据，以通过处理所述共检波点道集数据获得检波点的中波长静校正量。
步骤S16，判断步骤S14中的炮点的中波长静校正量和步骤S15中的检波点的中波长静校正量是否满足预设条件。当所述炮点的中波长静校正量和所述检波点的中波长静校正量不满足预设条件时，迭代步骤S13至步骤S15，直至所述炮点的中波长静校正量和所述检波点的中波长静校正量满足预设条件。
步骤S17，当所述炮点的中波长静校正量和所述检波点的中波长静校正量满足预设条件时，输出满足预设条件的炮点的中波长静校正量和检波点的中波长静校正量并应用。
其中，在步骤S11中，所述地震勘探数据是经过长波长静校正处理后的数据，将经过长波长静校正处理后的数据抽取成共检波点道集数据。在抽取到共检波点道集数据后执行步骤S12，处理所述共检波点道集数据，以获得检波点的中波长静校正量并应用于所述地震勘探数据。接下来，请参考图2，所述处理所述共检波点道集数据可以包括以下步骤：
步骤S121，对所述共检波点道集数据进行线性动校正。
行线性动校正可以是利用速度扫描方式，根据初至波校平的情况，选取线性动校正的速度以对所述共检波点道集数据进行线性动校正。
步骤S122，对线性动校正后的所述共检波点道集数据进行叠加，以获得叠加后的共检波点道集数据。
共检波点道集数据进行线性动校正后，避开近炮检距的直达波部分，保留稳定的折射波部分，以便同相叠加，获得叠加后的共检波点道集数据。
步骤S123，从叠加后的所述共检波点道集数据中提取检波点的中波长静校正量。
基于所述共检波点道集数据叠加后的初至波趋势，拾取用于提取中波长剩余静校正的零时间线和延迟时间线，将拾取的延迟时间线减去零时间线，即可获得检波点的中波长静校正量。
提取检波点的中波长静校正量后，将检波点的中波长静校正量应用于步骤S11中的所述地震勘探数据中，接着，执行步骤S13，从经过所述检波点的中波长静校正量应用后的所述地震勘探数据中抽取炮点道集数据。
为了获得炮点的中波长静校正量，在经过步骤S12中的检波点中波长静校正应用于所述地震勘探数据的基础上，将该地震勘探数据抽取成共炮点道集数据。
在抽取到共检波点道集数据后执行步骤S14，处理所述炮点道集数据，以获得炮点的中波长静校正量并应用于所述地震勘探数据。接下来，请参考图3，所述处理所述炮点道集数据可以包括以下步骤：
步骤S141，对所述炮点道集数据进行线性动校正。
行线性动校正可以是利用速度扫描方式，根据初至波校平的情况，选取线性动校正的速度以对所述炮点道集数据进行线性动校正。
步骤S142，对线性动校正后的所述炮点道集数据进行叠加，以获得叠加后的炮点道集数据。
炮点道集数据进行线性动校正后，避开近炮检距的直达波部分，保留稳定的折射波部分，以便同相叠加，获得叠加后的炮点道集数据。
步骤S143，从叠加后的所述炮点道集数据中提取炮点的中波长静校正量。
基于所述炮点道集数据叠加后的初至波趋势，拾取用于提取中波长剩余静校正的零时间线和延迟时间线，将拾取的延迟时间线减去零时间线，即可获得炮点的中波长静校正量。
提取炮点的中波长静校正量后，将炮点的中波长静校正量应用于步骤S11中的所述地震勘探数据中，接着，执行步骤S15，从经过所述炮点的中波长静校正量应用后的所述地震勘探数据中抽取共检波点道集数据，以通过处理所述共检波点道集数据获得检波点的中波长静校正量。
为了获得更精确的中波长的静校正量，再次在经过步骤S14中的炮点的中波长静校正应用于所述地震勘探数据的基础上，将该地震勘探数据抽取成共检波点道集数据。
接着，执行步骤S16，判断所述炮点的中波长静校正量和步骤S15中的所述检波点的中波长静校正量是否满足预设条件。当所述炮点的中波长静校正量和所述检波点的中波长静校正量不满足预设条件时，迭代步骤S13至步骤S15，直至所述炮点的中波长静校正量和所述检波点的中波长静校正量满足预设条件。
步骤S14中获得的炮点的中波长静校正量和步骤S15中获得的检波点的中波长静校正量是否满足预设条件，可以是通过查看经过炮点和检波点的中波长静校正量处理后的单炮记录判定。
当然，在其它的实施例中，步骤S14中获得的炮点的中波长静校正量和步骤S15中获得的检波点的中波长静校正量是否满足预设条件，也可以是通过查看经过炮点和检波点的中波长静校正量应用于地震勘探数据后的叠加剖面判定。
在判定炮点的中波长静校正量和检波点的中波长静校正量满足预设条件后，执行步骤S17，输出满足预设条件的炮点的中波长静校正量和检波点的中波长静校正量并应用。应用可包括以下几个方面，在可控震源激发无井深或炸药激发存在井深的应用情境下，将炮点的中波长静校正量应用于炮点，将检波点的中波长量应用于检波点。
针对复杂的地震勘探数据的处理，本公开通过采取交替迭代提取中波长静校正量的方法，可提高拾取的精度，进而提高处理成果构造的可靠性和资料的信噪比，为解释提供可靠的基础资料，使得本公开在黄土塬等静校正复杂探区地震数据处理中具有很好的应用前景。
值得说明的是，对于图1所示的方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本公开并不受所描述的动作顺序的限制。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于优选实施例，所涉及的动作并不一定是本公开所必须的。
实施例二
接下来，以黄土塬地震资料的数据为例，对本公开的地震勘探数据的处理方法进行阐述。
首先，在黄土塬地震资料的数据进行长波长静校正处理后，将处理后的黄土塬地震资料的数据抽取成如图4所示的共检波点道集数据。
接着，利用速度扫描方式，根据初至波校平的情况，选取线性动校正的速度以对图4所示的共检波点道集数据进行线性动校正，线性动校正后的所述共检波点道集数据如图5所示。
然后，在图5所示的所述共检波点道集数据基础上，避开近炮检距的直达波部分，保留稳定的折射波部分，进行同相叠加，获得叠加后的共检波点道集数据，叠加后的共检波点道集数据如图6所示。
接着，图6给出了基于共检波点初至波叠加剖面，如图6所示，基于所述共检波点道集数据叠加后的初至波趋势，拾取用于提取中波长剩余静校正的零时间线和延迟时间线，将拾取的延迟时间线减去零时间线，即可获得检波点的中波长静校正量，参考图7，图7为检波点中波长静校正量应用后叠加的共检波点道集数据示意图。
然后，将检波点的中波长静校正量应用于经过长波长静校正处理后的黄土塬地震资料的数据，并将应用后的黄土塬地震资料的数据抽取成共炮点道集数据，以便利用速度扫描方式，根据初至波校平的情况，选取线性动校正的速度以对所述炮点道集数据进行线性动校正。所述炮点道集数据进行线性动校正后如图8所示。
接着，在图8所示的所述炮点道集数据的基础上，避开近炮检距的直达波部分，保留稳定的折射波部分，进行同相叠加，获得叠加后的炮点道集数据。叠加后的炮点道集数据如图9所示。
然后，图9给出了基于炮点初至波叠加剖面，如图9所示，基于所述炮点道集数据叠加后的初至波趋势，拾取用于提取中波长剩余静校正的零时间线和延迟时间线，将拾取的延迟时间线减去零时间线，即可获得炮点的中波长静校正量，参考图10，图10为炮点中波长静校正量应用后叠加的炮点道集数据示意图。
接着，将获得的炮点的中波长静校正量应用于经过长波长静校正处理后的黄土塬地震资料的数据，并将应用后的黄土塬地震资料的数据再次抽取成共检波点道集数据，再次抽取后的共检波点道集数据如图11所示。
然后，判断再次获得的检波点的中波长静校正量和炮点的中波长静校正量是否满足预设条件。以下以两个判定方法为例：
其一、通过查看经过炮点和检波点的中波长静校正量处理后的单炮记录判定为例，请参考图12和图13，图12是只有长波长静校正量处理后的单炮记录，图13是长波长静校正量、炮点和检波点的中波长静校正量处理后的单炮记录，可以看出反射波同相轴的双曲线规律变的明显，即可以判定再次获得的检波点的中波长静校正量和炮点的中波长静校正量满足预设条件，进而输出满足预设条件的炮点的中波长静校正量和检波点的中波长静校正量并应用。
其二、通过查看经过炮点和检波点的中波长静校正量应用于地震勘探数据后的叠加剖面判定为例，请参考图14和图15，图14为没有经过炮点和检波点的中波长静校正量应用后的叠加剖面，该叠加剖面的中出现了串层现象，波组特征不理想。图15为经过炮点和检波点的中波长静校正量应用后的叠加剖面，可以看出该叠加剖面品质得到了明显提升，串层现象明显减少，波组特征理想，即可以判定再次获得的检波点的中波长静校正量和炮点的中波长静校正量满足预设条件，进而输出满足预设条件的炮点的中波长静校正量和检波点的中波长静校正量并应用。
当所述炮点的中波长静校正量和所述检波点的中波长静校正量不满足预设条件时，通过迭代上述抽取炮点和检波点的中波长静校正量的步骤，直至所述炮点和所述检波点的中波长静校正量满足预设条件。通过图13和图15可以看出，本公开提供的地震勘探数据的处理方法可在一定程度上解决了巨厚黄土塬地区的中波长静校正量难以精确提取的问题，获得了构造形态可靠且信噪比较高的处理成果。
实施例三
本公开还提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述任一项可选实施例所述的地震勘探数据的处理方法步骤。其中，在所述处理器上运行的煤矿监控系统的测试方法的计算机程序被执行时所实现的方法可参照本公开地震勘探数据的处理方法的具体实施例，此处不再赘述。
所述处理器可以是一种集成电路芯片，具有信息处理能力。所述处理器可以是通用处理器，包括中央处理器(Central Processing Unit，CPU)、网络处理器(NetworkProcessor，NP)等。
实施例四
本公开还提供一种地震勘探数据的处理装置，包括：
存储器，其上存储有计算机程序；以及
处理器，用于执行所述存储器中的所述计算机程序，以实现上述任一项可选实施例所述的地震勘探数据的处理方法步骤。
图16是根据一示例性实施例示出的一种的装置400的框图。如图16所示，该装置400可以包括：处理器401，存储器402，多媒体组件403，输入/输出(I/O)接口404，以及通信组件405。
其中，处理器401用于控制该装置400的整体操作，以完成上述的地震勘探数据的处理方法中的全部或部分步骤。存储器402用于存储各种类型的数据以支持在该装置400的操作，这些数据例如可以包括用于在该装置400上操作的任何应用程序或方法的指令，以及应用程序相关的数据。该存储器402可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(Static Random Access Memory，简称SRAM)，电可擦除可编程只读存储器(Electrically Erasable Programmable Read-Only Memory，简称EEPROM)，可擦除可编程只读存储器(Erasable Programmable Read-Only Memory，简称EPROM)，可编程只读存储器(Programmable Read-Only Memory，简称PROM)，只读存储器(Read-Only Memory，简称ROM)，磁存储器，快闪存储器，磁盘或光盘。多媒体组件403可以包括屏幕和音频组件。其中屏幕例如可以是触摸屏，音频组件用于输出和/或输入音频信号。例如，音频组件可以包括一个麦克风，麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器402或通过通信组件405发送。音频组件还包括至少一个扬声器，用于输出音频信号。I/O接口404为处理器401和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。通信组件405用于该装置400与其他设备之间进行有线或无线通信。无线通信，例如Wi-Fi，蓝牙，近场通信(Near Field Communication，简称NFC)，2G、3G或4G，或它们中的一种或几种的组合，因此相应的该通信组件405可以包括：Wi-Fi模块，蓝牙模块，NFC模块。
在一示例性实施例中，装置400可以被一个或多个应用专用集成电路(Application Specific Integrated Circuit，简称ASIC)、数字信号处理器(DigitalSignal Processor，简称DSP)、数字信号处理设备(Digital Signal Processing Device，简称DSPD)、可编程逻辑器件(Programmable Logic Device，简称PLD)、现场可编程门阵列(Field Programmable Gate Array，简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述的地震勘探数据的处理方法。
在另一示例性实施例中，还提供了一种包括程序指令的计算机可读存储介质，例如包括程序指令的存储器402，上述程序指令可由装置400的处理器401执行以完成上述的地震勘探数据的处理方法。
以上结合附图详细描述了本公开的优选实施方式，但是，本公开并不限于上述实施方式中的具体细节，在本公开的技术构思范围内，可以对本公开的技术方案进行多种简单变型，这些简单变型均属于本公开的保护范围。
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本公开对各种可能的组合方式不再另行说明。
此外，本公开的各种不同的实施方式之间也可以进行任意组合，只要其不违背本公开的思想，其同样应当视为本公开所公开的内容。

Claims (10)

1.一种地震勘探数据的处理方法，其特征在于，所述方法包括：
步骤S11：在获得长波长静校正处理后的地震勘探数据后，从所述地震勘探数据中抽取共检波点道集数据；
步骤S12：处理所述共检波点道集数据，以获得检波点的中波长静校正量并应用于所述地震勘探数据；
步骤S13：从经过所述检波点的中波长静校正量应用后的所述地震勘探数据中抽取炮点道集数据；
步骤S14：处理所述炮点道集数据，以获得炮点的中波长静校正量并应用于所述地震勘探数据；
步骤S15：从经过所述炮点的中波长静校正量应用后的所述地震勘探数据中抽取共检波点道集数据，以通过处理所述共检波点道集数据获得检波点的中波长静校正量；
步骤S16：判断所述炮点的中波长静校正量和步骤S15中的所述检波点的中波长静校正量是否满足预设条件；当所述炮点的中波长静校正量和所述检波点的中波长静校正量不满足预设条件时，迭代步骤S13至步骤S15，直至所述炮点的中波长静校正量和所述检波点的中波长静校正量满足预设条件；
步骤S17：当所述炮点的中波长静校正量和所述检波点的中波长静校正量满足预设条件时，输出满足预设条件的炮点的中波长静校正量和检波点的中波长静校正量并应用。
2.根据权利要求1所述的方法，其特征在于，处理所述共检波点道集数据，以获得检波点的中波长静校正量，包括：
对所述共检波点道集数据进行线性动校正；
对线性动校正后的所述共检波点道集数据进行叠加，以获得叠加后的共检波点道集数据；
从叠加后的所述共检波点道集数据中提取检波点的中波长静校正量。
3.根据权利要求2所述的方法，其特征在于，对所述共检波点道集数据进行线性动校正，包括：
利用速度扫描方式，根据初至波校平的情况，选取线性动校正的速度以对所述共检波点道集数据进行线性动校正。
4.根据权利要求2所述的方法，其特征在于，对线性动校正后的所述共检波点道集数据进行叠加，以获得叠加后的共检波点道集数据，包括：
避开近炮检距的直达波部分，保留稳定的折射波部分，以对线性动校正后的所述共检波点道集数据进行叠加。
5.根据权利要求2所述的方法，其特征在于，从叠加后的所述共检波点道集数据中提取检波点的中波长静校正量，包括：
基于所述共检波点道集数据叠加后的初至波趋势，拾取用于提取中波长剩余静校正的零时间线和延迟时间线；
将拾取的延迟时间线减去零时间线，以获得检波点的中波长静校正量。
6.根据权利要求1所述的方法，其特征在于，处理所述炮点道集数据，以获得炮点的中波长静校正量，包括：
对所述炮点道集数据进行线性动校正；
对线性动校正后的所述炮点道集数据进行叠加，以获得叠加后的炮点道集数据；
从叠加后的所述炮点道集数据中提取炮点的中波长静校正量。
7.根据权利要求6所述的方法，其特征在于，对所述炮点道集数据进行线性动校正，包括：
利用速度扫描方式，根据初至波校平的情况，选取线性动校正的速度以对所述炮点道集数据进行线性动校正。
8.根据权利要求6所述的方法，其特征在于，对线性动校正后的所述炮点道集数据进行叠加，以获得叠加后的炮点道集数据，包括：
避开近炮检距的直达波部分，保留稳定的折射波部分，以对线性动校正后的所述炮点道集数据进行叠加。
9.根据权利要求6所述的方法，其特征在于，从叠加后的所述炮点道集数据中提取炮点的中波长静校正量，包括：
基于所述炮点道集数据叠加后的初至波趋势，拾取用于提取中波长剩余静校正的零时间线和延迟时间线；
将拾取的延迟时间线减去零时间线，以获得炮点的中波长静校正量。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1至9中任一项所述方法的步骤。
CN201910817078.5A 2019-08-30 2019-08-30 地震勘探数据的处理方法及存储介质 Active CN112444879B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910817078.5A CN112444879B (zh) 2019-08-30 2019-08-30 地震勘探数据的处理方法及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910817078.5A CN112444879B (zh) 2019-08-30 2019-08-30 地震勘探数据的处理方法及存储介质

Publications (2)

Publication Number Publication Date
CN112444879A CN112444879A (zh) 2021-03-05
CN112444879B true CN112444879B (zh) 2022-10-14

Family

ID=74734836

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910817078.5A Active CN112444879B (zh) 2019-08-30 2019-08-30 地震勘探数据的处理方法及存储介质

Country Status (1)

Country Link
CN (1) CN112444879B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133245A (zh) * 2014-07-22 2014-11-05 中国石油天然气集团公司 一种地震资料的静校正方法及系统
CN105093326A (zh) * 2014-05-23 2015-11-25 中国石油化工股份有限公司 一种地震勘探资料中波长剩余静校正处理方法
CN107526101A (zh) * 2017-07-07 2017-12-29 中国地质调查局油气资源调查中心 一种获取地震反射波的采集和处理方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1128373C (zh) * 2001-01-21 2003-11-19 周熙襄 用初至波对地震勘探资料进行短波长静校正处理的方法
US7257492B2 (en) * 2005-08-26 2007-08-14 Westerngeco L.L. Handling of static corrections in multiple prediction
CN101201409B (zh) * 2006-12-14 2011-05-25 中国石油天然气集团公司 一种地震数据变相位校正方法
US20150212223A1 (en) * 2014-01-28 2015-07-30 Cgg Services Sa Methods and apparatus for determining residual static corrections using individual ranges

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105093326A (zh) * 2014-05-23 2015-11-25 中国石油化工股份有限公司 一种地震勘探资料中波长剩余静校正处理方法
CN104133245A (zh) * 2014-07-22 2014-11-05 中国石油天然气集团公司 一种地震资料的静校正方法及系统
CN107526101A (zh) * 2017-07-07 2017-12-29 中国地质调查局油气资源调查中心 一种获取地震反射波的采集和处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
三维炮检域迭代直接静校正技术及应用;黄明忠等;《勘探地球物理进展》;20100228;第33卷(第1期);第48-49页 *
复杂地表低信噪比地震数据处理研究;马义忠等;《地球物理学进展》;20081231;第23卷(第6期);第1849-1858页 *

Also Published As

Publication number Publication date
CN112444879A (zh) 2021-03-05

Similar Documents

Publication Publication Date Title
CN104122587B (zh) 一种基于炮检距域的异常初至识别方法以及系统
CN107450108B (zh) 甜点区的确定方法和装置
CN104199103A (zh) 转换波静校正方法及装置
CN104330827A (zh) 一种地表模型静校正的处理方法
CN112444879B (zh) 地震勘探数据的处理方法及存储介质
CN103487834A (zh) 转换波共检波点叠加静校正方法
CN111399037B (zh) 高速顶界面提取的方法和装置
CN109633744B (zh) 地震子波的提取方法、装置、设备及存储介质
EP4390466A1 (en) Vsp-based level calibration method and apparatus for depth-domain seismic profile
CN111208558A (zh) 超深低幅度三维地质构造的建立方法及装置
CN112255679B (zh) 地震资料绕射深度偏移处理方法及装置
CN116027431A (zh) 储层含气性预测方法、装置、存储介质及电子设备
CN111596363B (zh) 初至时间飘移校正的方法和装置
CN106772596A (zh) 一种确定叠前时间偏移速度场的方法及装置
CN103837890B (zh) 获取地震数据的方法及设备
CN105789081A (zh) 一种加速wat测试的系统和方法
CN112462421A (zh) 储层信息预测方法、装置、电子设备及存储介质
CN104749627A (zh) 基于相似性的微地震信号凸显方法
Changliang et al. VTI Media Kirchhoff Pre-Stack Depth Migration Technology and Its Application
CN108260359A (zh) 一种地震信号处理方法、装置和系统
CN104199107A (zh) 基于垂直地震多波数据的钻前深度预测方法以及系统
CN105242308A (zh) 一种宽方位地震数据的叠加方法和装置
CN114114417B (zh) 地震成像方法及装置
CN112824939B (zh) 地震数据的剩余静校正方法、装置、设备及可读存储介质
CN112505748B (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