CN112255679B - 地震资料绕射深度偏移处理方法及装置 - Google Patents

地震资料绕射深度偏移处理方法及装置 Download PDF

Info

Publication number
CN112255679B
CN112255679B CN202011154201.9A CN202011154201A CN112255679B CN 112255679 B CN112255679 B CN 112255679B CN 202011154201 A CN202011154201 A CN 202011154201A CN 112255679 B CN112255679 B CN 112255679B
Authority
CN
China
Prior art keywords
sampling point
position information
offset
determining
point
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
CN202011154201.9A
Other languages
English (en)
Other versions
CN112255679A (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 National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN202011154201.9A priority Critical patent/CN112255679B/zh
Publication of CN112255679A publication Critical patent/CN112255679A/zh
Application granted granted Critical
Publication of CN112255679B publication Critical patent/CN112255679B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/70—Other details related to processing
    • G01V2210/74—Visualisation of seismic data

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

地震资料绕射深度偏移处理方法及装置
技术领域
本申请涉及地质勘探领域，具体涉及一种地震资料绕射深度偏移处理方法及装置。
背景技术
绕射是地震波传播过程中产生的一种地震波，在断裂发育的地层，绕射波非常规发育，尤其是地层断点位置，绕射波异常发育。
发明人发现，在现有技术中，地震资料偏移成像处理采用的是基于反射波的偏移算法，偏移成像是反射波的成像，这就使断点产生的绕射波绝大部分能量不能被正确利用，往往能量最强的部分无法得到有效的成像，而是使用相对较弱的反射部分来偏移成像，由此导致偏移的断点成像不精确。
发明内容
针对现有技术中的问题，本申请提供一种地震资料绕射深度偏移处理方法及装置，能够有效利用断点产生的绕射波进行偏移成像，改善地震资料断层断点的成像精度。
为了解决上述问题中的至少一个，本申请提供以下技术方案：
第一方面，本申请提供一种地震资料绕射深度偏移处理方法，包括：
根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息；
确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息；
根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像。
进一步地，所述根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息包括：
根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度，确定各层出射角，其中，所述预设三维坐标系是由地面网格和纵向深度构建的；
根据所述各层出射角和目标采样点在预设三维坐标系内对应的纵波深度域各层厚度，确定各层横向距离；
根据所述各层横向距离，确定参照采样点与所述目标采样点的横纵坐标变化量，其中，所述参照采样点为所述目标采样点在所述预设三维坐标系中以法线为中心，围绕所述法线方向以至少一个方位角和至少一个出射角发射射线并到达地面时所述地面上的点；
根据所述横纵坐标变化量，确定所述参照采样点位置信息。
进一步地，所述确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息，包括：
确定地面上以所述参照采样点位置信息为圆心、以设定长度为半径的圆形范围内的各炮检点位置信息，并获取地震资料中各所述炮检点位置信息对应的旅行时信息。
进一步地，所述根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像，包括：
对与所述炮检点位置信息对应的地震道数据在所述旅行时信息中各旅行时总和的时间位置处的样点值进行取样叠加处理；
根据经过所述取样叠加处理后的样点值，确定所述目标采样点的绕射深度偏移成像。
第二方面，本申请提供一种地震资料绕射深度偏移处理装置，包括：
参照采样点位置信息确定模块，用于根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息；
炮检点位置信息和旅行时信息确定模块，用于确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息；
绕射深度偏移成像模块，用于根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像。
进一步地，所述参照采样点位置信息确定模块包括：
各层出射角确定单元，用于根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度，确定各层出射角，其中，所述预设三维坐标系是由地面网格和纵向深度构建的；
各层横向距离确定单元，用于根据所述各层出射角和目标采样点在预设三维坐标系内对应的纵波深度域各层厚度，确定各层横向距离；
横纵坐标变化量确定单元，用于根据所述各层横向距离，确定参照采样点与所述目标采样点的横纵坐标变化量，其中，所述参照采样点为所述目标采样点在所述预设三维坐标系中以法线为中心，围绕所述法线方向以至少一个方位角和至少一个出射角发射射线并到达地面时所述地面上的点；
参照采样点位置信息确定单元，用于根据所述横纵坐标变化量，确定所述参照采样点位置信息。
进一步地，所述炮检点位置信息和旅行时信息确定模块包括：
炮检点筛选单元，用于确定地面上以所述参照采样点位置信息为圆心、以设定长度为半径的圆形范围内的各炮检点位置信息，并获取地震资料中各所述炮检点位置信息对应的旅行时信息。
进一步地，所述绕射深度偏移成像模块包括：
取样叠加处理单元，用于对与所述炮检点位置信息对应的地震道数据在所述旅行时信息中各旅行时总和的时间位置处的样点值进行取样叠加处理；
绕射深度偏移成像单元，用于根据经过所述取样叠加处理后的样点值，确定所述目标采样点的绕射深度偏移成像。
第三方面，本申请提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的地震资料绕射深度偏移处理方法的步骤。
第四方面，本申请提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的地震资料绕射深度偏移处理方法的步骤。
由上述技术方案可知，本申请提供一种地震资料绕射深度偏移处理方法及装置，通过目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息，并根据与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息，确定目标采样点的绕射深度偏移成像，以此通过同相叠加实现了绕射波偏移成像，改善了断层、断点成像精度，打破了常规基于反射波偏移射线数量的限制，极大的提高了参与叠加的地震道数，提高了信噪比。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例中的地震资料绕射深度偏移处理方法的流程示意图之一；
图2为本申请实施例中的地震资料绕射深度偏移处理方法的流程示意图之二；
图3为本申请实施例中的地震资料绕射深度偏移处理方法的流程示意图之三；
图4为本申请实施例中的地震资料绕射深度偏移处理装置的结构图之一；
图5为本申请实施例中的地震资料绕射深度偏移处理装置的结构图之二；
图6为本申请实施例中的地震资料绕射深度偏移处理装置的结构图之三；
图7为本申请实施例中的地震资料绕射深度偏移处理装置的结构图之四；
图8为本申请一具体实施例中的绕射偏移原理示意图；
图9为本申请实施例中的电子设备的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
考虑到在现有技术中地震资料偏移成像处理采用的是基于反射波的偏移算法，偏移成像是反射波的成像，这就使断点产生的绕射波绝大部分能量不能被正确利用，往往能量最强的部分无法得到有效的成像，而是使用相对较弱的反射部分来偏移成像，由此导致偏移的断点成像不精确的问题，本申请提供一种地震资料绕射深度偏移处理方法及装置，通过目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息，并根据与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息，确定目标采样点的绕射深度偏移成像，以此通过同相叠加实现了绕射波偏移成像，改善了断层、断点成像精度，打破了常规基于反射波偏移射线数量的限制，极大的提高了参与叠加的地震道数，提高了信噪比。
参见图8绕射偏移原理图可知，通常断层、断点绕射异常发育，从原理图上可以看出常规偏移是基于反射波的偏移，偏移算法只利用了反射射线一条射线，能量较低，而绕射偏移则利用了选取范围内所有地震道的绕射能量进行叠加，能够极大的提高了断层、断点偏移成像精度。
为了能够有效利用断点产生的绕射波进行偏移成像，改善地震资料断层断点的成像精度，本申请提供一种地震资料绕射深度偏移处理方法的实施例，参见图1，所述地震资料绕射深度偏移处理方法具体包含有如下内容：
步骤S101：根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息。
可选的，首先可以采用现有层析反演处理方法获得精确的纵波深度域各层速度、各层厚度和常规纵波深度偏移剖面，在深度偏移叠加剖面上还可获得地层的倾角信息。准备好用于偏移的道集。
可选的，在地面网格和纵向深度构成的三维空间(也即三维坐标系)内，取网格中的一点(即目标采样点)作为激发点，围绕地层的法线方向从该点向上发射曲线，直到追踪到地面(地面上对应的点即所述参照采样点)。
步骤S102：确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息。
可选的，在地面在一定误差范围内遇到炮点或检波点则记录下该射线对应炮点或检波点的旅行时，给定最初射线和法线的夹角范围来确定需要追踪射线的最大范围。
步骤S103：根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像。
可选的，在地面网格寻找与法线出射点距离最近的检波点，以这个检波点为中心划定一个范围圆，找出圆范围内所有的地震道，取地震道炮检点旅行时之和处样点值进行叠加。叠加样点值置于地下激发点(即所述目标采样点)处，由此实现快速绕射波深度偏移成像。
从上述描述可知，本申请实施例提供的地震资料绕射深度偏移处理方法，能够通过目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息，并根据与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息，确定目标采样点的绕射深度偏移成像，以此通过同相叠加实现了绕射波偏移成像，改善了断层、断点成像精度，打破了常规基于反射波偏移射线数量的限制，极大的提高了参与叠加的地震道数，提高了信噪比。
为了能够准确确定目标采样点的绕射在相关地震道存在的位置，在本申请的地震资料绕射深度偏移处理方法的一实施例中，参见图2，上述步骤S101还可以具体包含如下内容：
步骤S201：根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度，确定各层出射角，其中，所述预设三维坐标系是由地面网格和纵向深度构建的。
步骤S202：根据所述各层出射角和目标采样点在预设三维坐标系内对应的纵波深度域各层厚度，确定各层横向距离。
步骤S203：根据所述各层横向距离，确定参照采样点与所述目标采样点的横纵坐标变化量，其中，所述参照采样点为所述目标采样点在所述预设三维坐标系中以法线为中心，围绕所述法线方向以至少一个方位角和至少一个出射角发射射线并到达地面时所述地面上的点。
步骤S204：根据所述横纵坐标变化量，确定所述参照采样点位置信息。
可选的，首先，通过纵波均方根速度转化成层速度，可以确定每一层的厚度和速度。
然后，在地下cmp某一样点M(X1,Y1)围绕法线方向，以法线为中心，沿不同方位角α(0到359度)，和不同出射角β(0到指定角度，小于90度)向上发射射线，在每一个地层，使用公式依次计算各层的出射角，直到地面参照采样点N(X2,Y2)。
再然后，通过各层的厚度和出射角度，使用公式x＝h*tanβ来计算每一层的横向距离，进而，通过公式来计算M点到N点的距离；通过公式Δx＝X*sinα来计算M点到N点的横坐标变化量；通过公式X2＝X1+ΔX计算出地面N点的横坐标；通过公式Δy＝X*cosα来计算M点到N点的纵坐标变化量；通过公式Y2＝Y1+ΔY计算出地面N点的纵坐标。
此外，还可以通过公式li＝hi/cosβi来计算每一层的旅行距离；通过公式ti＝li/vi来分别计算每一层的旅行时。
为了能够准确确定地震道绕射能量的选取范围，在本申请的地震资料绕射深度偏移处理方法的一实施例中，上述步骤S102还可以具体包含如下内容：
确定地面上以所述参照采样点位置信息为圆心、以设定长度为半径的圆形范围内的各炮检点位置信息，并获取地震资料中各所述炮检点位置信息对应的旅行时信息。
可选的，本申请通过给定一定误差值，在误差值内寻找和N点坐标吻合的炮点和检波点。通过改变方位角和出射角扫描出所有范围内的地面炮点和检波点，记录相应的旅行时。
具体的，在地面网格寻找与法线出射点距离最近的检波点，以这个检波点为中心划定一个范围圆，找出圆范围内所有的地震道及对应的旅行时。
为了能够准确通过同相叠加实现了绕射波偏移成像，在本申请的地震资料绕射深度偏移处理方法的一实施例中，参见图3，上述步骤S103还可以具体包含如下内容：
步骤S301：对与所述炮检点位置信息对应的地震道数据在所述旅行时信息中各旅行时总和的时间位置处的样点值进行取样叠加处理。
步骤S302：根据经过所述取样叠加处理后的样点值，确定所述目标采样点的绕射深度偏移成像。
可选的，把寻找到的和M点有关的所有地震道对应的炮点和检波点的旅行时之和的时间位置处的样点值叠加到一起放置到M点位置处，就完成了M点的绕射深度偏移处理。把地下所有cmp样点值都按照以上步骤进行取样叠加，就完成了整个地震数据的绕射深度偏移处理。
为了能够有效利用断点产生的绕射波进行偏移成像，改善地震资料断层断点的成像精度，本申请提供一种用于实现所述地震资料绕射深度偏移处理方法的全部或部分内容的地震资料绕射深度偏移处理装置的实施例，参见图4，所述地震资料绕射深度偏移处理装置具体包含有如下内容：
参照采样点位置信息确定模块10，用于根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息。
炮检点位置信息和旅行时信息确定模块20，用于确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息。
绕射深度偏移成像模块30，用于根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像。
从上述描述可知，本申请实施例提供的地震资料绕射深度偏移处理装置，能够通过目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息，并根据与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息，确定目标采样点的绕射深度偏移成像，以此通过同相叠加实现了绕射波偏移成像，改善了断层、断点成像精度，打破了常规基于反射波偏移射线数量的限制，极大的提高了参与叠加的地震道数，提高了信噪比。
为了能够准确确定目标采样点的绕射在相关地震道存在的位置，在本申请的地震资料绕射深度偏移处理装置的一实施例中，参见图5，所述参照采样点位置信息确定模块包括：
各层出射角确定单元11，用于根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度，确定各层出射角，其中，所述预设三维坐标系是由地面网格和纵向深度构建的。
各层横向距离确定单元12，用于根据所述各层出射角和目标采样点在预设三维坐标系内对应的纵波深度域各层厚度，确定各层横向距离。
横纵坐标变化量确定单元13，用于根据所述各层横向距离，确定参照采样点与所述目标采样点的横纵坐标变化量，其中，所述参照采样点为所述目标采样点在所述预设三维坐标系中以法线为中心，围绕所述法线方向以至少一个方位角和至少一个出射角发射射线并到达地面时所述地面上的点。
参照采样点位置信息确定单元14，用于根据所述横纵坐标变化量，确定所述参照采样点位置信息。
为了能够准确确定地震道绕射能量的选取范围，在本申请的地震资料绕射深度偏移处理装置的一实施例中，参见图6，所述炮检点位置信息和旅行时信息确定模块20包括：
炮检点筛选单元21，用于确定地面上以所述参照采样点位置信息为圆心、以设定长度为半径的圆形范围内的各炮检点位置信息，并获取地震资料中各所述炮检点位置信息对应的旅行时信息。
为了能够准确通过同相叠加实现了绕射波偏移成像，在本申请的地震资料绕射深度偏移处理装置的一实施例中，参见图7，所述绕射深度偏移成像模块30包括：
取样叠加处理单元31，用于对与所述炮检点位置信息对应的地震道数据在所述旅行时信息中各旅行时总和的时间位置处的样点值进行取样叠加处理。
绕射深度偏移成像单元32，用于根据经过所述取样叠加处理后的样点值，确定所述目标采样点的绕射深度偏移成像。
为了更进一步说明本方案，本申请还提供一种应用上述地震资料绕射深度偏移处理装置实现地震资料绕射深度偏移处理方法的具体应用实例，具体包含有如下内容：
步骤1)纵波均方根速度转化成层速度。确定每一层的厚度和速度。
步骤2)在地下cmp某一样点M(X1,Y1)围绕法线方向，以法线为中心，沿不同方位角α(0到359度)，和不同出射角β(0到指定角度，小于90度)向上发射射线。
步骤3)在每一个地层，使用公式依次计算各层的出射角，直到地面某点N(X2,Y2)。
步骤4)通过各层的厚度和出射角度，使用公式x＝h*tanβ来计算每一层的横向距离。
步骤5)通过公式来计算M点到N点的距离。
步骤6)通过公式Δx＝X*sinα来计算M点到N点的横坐标变化量。
步骤7)通过公式X2＝X1+ΔX计算出地面N点的横坐标。
步骤8)通过公式Δy＝X*cosα来计算M点到N点的纵坐标变化量。
步骤9)通过公式Y2＝Y1+ΔY计算出地面N点的纵坐标。
步骤10)通过公式li＝hi/cosβi来计算每一层的旅行距离。
步骤11)通过公式ti＝li/vi来分别计算每一层的旅行时。
步骤12)给定一定误差值，在误差值内寻找和N点坐标吻合的炮点和检波点。
步骤13)通过改变方位角和出射角扫描出所有范围内的地面炮点和检波点，记录相应的旅行时。
步骤14)把寻找到的和M点有关的所有地震道对应的炮点和检波点的旅行时之和的时间位置处的样点值叠加到一起放置到M点位置处，就完成了M点的绕射深度偏移处理。
步骤15)把地下所有cmp样点值都按照以上步骤进行取样叠加，就完成了整个地震数据的绕射深度偏移处理。
由上述描述可知，本申请至少还可以实现如下技术效果：
在纵波精确成像的基础上通过扫描来确定地下某点的绕射在相关地震道存在的位置，通过同相叠加实现了绕射波偏移成像，改善了断层、断点成像精度，该算法打破了常规基于反射波偏移射线数量的限制，极大的提高了参与叠加的地震道数，提高了信噪比。
从硬件层面来说，为了能够有效利用断点产生的绕射波进行偏移成像，改善地震资料断层断点的成像精度，本申请提供一种用于实现所述地震资料绕射深度偏移处理方法中的全部或部分内容的电子设备的实施例，所述电子设备具体包含有如下内容：
处理器(processor)、存储器(memory)、通信接口(Communications Interface)和总线；其中，所述处理器、存储器、通信接口通过所述总线完成相互间的通信；所述通信接口用于实现地震资料绕射深度偏移处理装置与核心业务系统、用户终端以及相关数据库等相关设备之间的信息传输；该逻辑控制器可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该逻辑控制器可以参照实施例中的地震资料绕射深度偏移处理方法的实施例，以及地震资料绕射深度偏移处理装置的实施例进行实施，其内容被合并于此，重复之处不再赘述。
可以理解的是，所述用户终端可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中，所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。
在实际应用中，地震资料绕射深度偏移处理方法的部分可以在如上述内容所述的电子设备侧执行，也可以所有的操作都在所述客户端设备中完成。具体可以根据所述客户端设备的处理能力，以及用户使用场景的限制等进行选择。本申请对此不作限定。若所有的操作都在所述客户端设备中完成，所述客户端设备还可以包括处理器。
上述的客户端设备可以具有通信模块(即通信单元)，可以与远程的服务器进行通信连接，实现与所述服务器的数据传输。所述服务器可以包括任务调度中心一侧的服务器，其他的实施场景中也可以包括中间平台的服务器，例如与任务调度中心服务器有通信链接的第三方服务器平台的服务器。所述的服务器可以包括单台计算机设备，也可以包括多个服务器组成的服务器集群，或者分布式装置的服务器结构。
图9为本申请实施例的电子设备9600的系统构成的示意框图。如图9所示，该电子设备9600可以包括中央处理器9100和存储器9140；存储器9140耦合到中央处理器9100。值得注意的是，该图9是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。
一实施例中，地震资料绕射深度偏移处理方法功能可以被集成到中央处理器9100中。其中，中央处理器9100可以被配置为进行如下控制：
步骤S101：根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息。
步骤S102：确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息。
步骤S103：根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像。
从上述描述可知，本申请实施例提供的电子设备，通过目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息，并根据与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息，确定目标采样点的绕射深度偏移成像，以此通过同相叠加实现了绕射波偏移成像，改善了断层、断点成像精度，打破了常规基于反射波偏移射线数量的限制，极大的提高了参与叠加的地震道数，提高了信噪比。
在另一个实施方式中，地震资料绕射深度偏移处理装置可以与中央处理器9100分开配置，例如可以将地震资料绕射深度偏移处理装置配置为与中央处理器9100连接的芯片，通过中央处理器的控制来实现地震资料绕射深度偏移处理方法功能。
如图9所示，该电子设备9600还可以包括：通信模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是，电子设备9600也并不是必须要包括图9中所示的所有部件；此外，电子设备9600还可以包括图9中没有示出的部件，可以参考现有技术。
如图9所示，中央处理器9100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。
其中，存储器9140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序，以实现信息存储或处理等。
输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器9140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142，该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。
存储器9140还可以包括数据存储部9143，该数据存储部9143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通信模块(发送机/接收机)9110耦合到中央处理器9100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块9110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132，以经由扬声器9131提供音频输出，并接收来自麦克风9132的音频输入，从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器9130还耦合到中央处理器9100，从而使得可以通过麦克风9132能够在本机上录音，且使得可以通过扬声器9131来播放本机上存储的声音。
本申请的实施例还提供能够实现上述实施例中的执行主体为服务器或客户端的地震资料绕射深度偏移处理方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的执行主体为服务器或客户端的地震资料绕射深度偏移处理方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
步骤S101：根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息。
步骤S102：确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息。
步骤S103：根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像。
从上述描述可知，本申请实施例提供的计算机可读存储介质，通过目标采样点在预设三维坐标系内对应的纵波深度域各层速度和各层厚度，确定参照采样点位置信息，并根据与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息，确定目标采样点的绕射深度偏移成像，以此通过同相叠加实现了绕射波偏移成像，改善了断层、断点成像精度，打破了常规基于反射波偏移射线数量的限制，极大的提高了参与叠加的地震道数，提高了信噪比。
本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(装置)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种地震资料绕射深度偏移处理方法，其特征在于，所述方法包括：
根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度，确定各层出射角，其中，所述预设三维坐标系是由地面网格和纵向深度构建的；
根据所述各层出射角和目标采样点在预设三维坐标系内对应的纵波深度域各层厚度，确定各层横向距离；
根据所述各层横向距离，确定参照采样点与所述目标采样点的横纵坐标变化量，其中，所述参照采样点为所述目标采样点在所述预设三维坐标系中以法线为中心，围绕所述法线方向以至少一个方位角和至少一个出射角发射射线并到达地面时所述地面上的点；
根据所述横纵坐标变化量，确定所述参照采样点位置信息；
确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息；
根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像。
2.根据权利要求1所述的地震资料绕射深度偏移处理方法，其特征在于，所述确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息，包括：
确定地面上以所述参照采样点位置信息为圆心、以设定长度为半径的圆形范围内的各炮检点位置信息，并获取地震资料中各所述炮检点位置信息对应的旅行时信息。
3.根据权利要求1所述的地震资料绕射深度偏移处理方法，其特征在于，所述根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像，包括：
对与所述炮检点位置信息对应的地震道数据在所述旅行时信息中各旅行时总和的时间位置处的样点值进行取样叠加处理；
根据经过所述取样叠加处理后的样点值，确定所述目标采样点的绕射深度偏移成像。
4.一种地震资料绕射深度偏移处理装置，其特征在于，包括：
参照采样点位置信息确定模块，用于：
根据目标采样点在预设三维坐标系内对应的纵波深度域各层速度，确定各层出射角，其中，所述预设三维坐标系是由地面网格和纵向深度构建的；
根据所述各层出射角和目标采样点在预设三维坐标系内对应的纵波深度域各层厚度，确定各层横向距离；
根据所述各层横向距离，确定参照采样点与所述目标采样点的横纵坐标变化量，其中，所述参照采样点为所述目标采样点在所述预设三维坐标系中以法线为中心，围绕所述法线方向以至少一个方位角和至少一个出射角发射射线并到达地面时所述地面上的点；
根据所述横纵坐标变化量，确定所述参照采样点位置信息；
炮检点位置信息和旅行时信息确定模块，用于确定与所述参照采样点位置信息对应的炮检点位置信息以及与所述炮检点位置信息对应的旅行时信息；
绕射深度偏移成像模块，用于根据所述炮检点位置信息和所述旅行时信息，确定所述目标采样点的绕射深度偏移成像。
5.根据权利要求4所述的地震资料绕射深度偏移处理装置，其特征在于，所述炮检点位置信息和旅行时信息确定模块包括：
炮检点筛选单元，用于确定地面上以所述参照采样点位置信息为圆心、以设定长度为半径的圆形范围内的各炮检点位置信息，并获取地震资料中各所述炮检点位置信息对应的旅行时信息。
6.根据权利要求4所述的地震资料绕射深度偏移处理装置，其特征在于，所述绕射深度偏移成像模块包括：
取样叠加处理单元，用于对与所述炮检点位置信息对应的地震道数据在所述旅行时信息中各旅行时总和的时间位置处的样点值进行取样叠加处理；
绕射深度偏移成像单元，用于根据经过所述取样叠加处理后的样点值，确定所述目标采样点的绕射深度偏移成像。
7.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至3任一项所述的地震资料绕射深度偏移处理方法的步骤。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至3任一项所述的地震资料绕射深度偏移处理方法的步骤。
CN202011154201.9A 2020-10-26 2020-10-26 地震资料绕射深度偏移处理方法及装置 Active CN112255679B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011154201.9A CN112255679B (zh) 2020-10-26 2020-10-26 地震资料绕射深度偏移处理方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011154201.9A CN112255679B (zh) 2020-10-26 2020-10-26 地震资料绕射深度偏移处理方法及装置

Publications (2)

Publication Number Publication Date
CN112255679A CN112255679A (zh) 2021-01-22
CN112255679B true CN112255679B (zh) 2023-09-26

Family

ID=74262964

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011154201.9A Active CN112255679B (zh) 2020-10-26 2020-10-26 地震资料绕射深度偏移处理方法及装置

Country Status (1)

Country Link
CN (1) CN112255679B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115616660B (zh) * 2022-12-20 2023-05-19 广州海洋地质调查局三亚南海地质研究所 利用绕射波监测海域碳封存项目二氧化碳泄露情况的方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237940A (zh) * 2014-09-29 2014-12-24 中国石油天然气股份有限公司 一种基于动力学特征的绕射波成像方法及装置
CN106772583A (zh) * 2017-01-10 2017-05-31 中国科学院地质与地球物理研究所 一种地震绕射波分离方法和装置
CN106772592A (zh) * 2016-11-10 2017-05-31 中国矿业大学(北京) 绕射波聚焦能量的分析方法及装置
CN107505653A (zh) * 2017-08-30 2017-12-22 中国石油集团东方地球物理勘探有限责任公司 一种确定叠前偏移时间结果的方法和装置
CN107861154A (zh) * 2017-11-08 2018-03-30 中国矿业大学(北京) 断层走向识别方法和装置
CN108845358A (zh) * 2018-09-26 2018-11-20 中国矿业大学(北京) 断层及构造异常体识别方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8209125B2 (en) * 2007-03-12 2012-06-26 Geomage (2003) Ltd. Method for identifying and analyzing faults/fractures using reflected and diffracted waves
CN105607121B (zh) * 2016-02-02 2016-12-21 中国矿业大学(北京) 一种煤炭陷落柱识别方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237940A (zh) * 2014-09-29 2014-12-24 中国石油天然气股份有限公司 一种基于动力学特征的绕射波成像方法及装置
CN106772592A (zh) * 2016-11-10 2017-05-31 中国矿业大学(北京) 绕射波聚焦能量的分析方法及装置
CN106772583A (zh) * 2017-01-10 2017-05-31 中国科学院地质与地球物理研究所 一种地震绕射波分离方法和装置
CN107505653A (zh) * 2017-08-30 2017-12-22 中国石油集团东方地球物理勘探有限责任公司 一种确定叠前偏移时间结果的方法和装置
CN107861154A (zh) * 2017-11-08 2018-03-30 中国矿业大学(北京) 断层走向识别方法和装置
CN108845358A (zh) * 2018-09-26 2018-11-20 中国矿业大学(北京) 断层及构造异常体识别方法及装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Diffraction imaging and time-migration velocity analysis using oriented velocity continuation;Luke Decker等;GEOPHYSICS;全文 *
倾角域逆时偏移绕射波成像方法;汪天池等;石油地球物理勘探;全文 *
倾角成像道集中反射波和绕射波特征分析及成像质量改善方法研究;徐辉;;石油物探(第02期);全文 *
徐辉.倾角成像道集中反射波和绕射波特征分析及成像质量改善方法研究.石油物探.2015,全文. *

Also Published As

Publication number Publication date
CN112255679A (zh) 2021-01-22

Similar Documents

Publication Publication Date Title
CN105116448B (zh) 一种转换波方位各向异性校正方法及装置
CN107576986B (zh) 一种确定海底地震反射波旅行时的方法和装置
CN104656139A (zh) 一种转换波的静校正方法及设备
CN112255679B (zh) 地震资料绕射深度偏移处理方法及装置
US20160320502A1 (en) Electroformed nickel-chromium alloy
CN112147688A (zh) 一种储层油气识别方法和装置
US20150153466A1 (en) Source Start Time Determination
US11899147B2 (en) Method and system for seismic denoising using omnifocal reformation
CN112462421B (zh) 储层信息预测方法、装置、电子设备及存储介质
CN112394408A (zh) 各向异性介质裂缝预测方法及装置
CN114578428A (zh) 一种速度场高效构建方法及装置
CN112379433A (zh) 基于双井微测井数据的近地表品质因子获取方法及装置
CN111679324A (zh) 地震数据零相位化处理方法、装置、设备和存储介质
CN112505748B (zh) 炮域反射波拾取方法及装置
EP3987317A1 (en) Marine survey source firing control
CN114755720B (zh) 一种页岩储层地震预测方法、装置、设备及介质
Petrov et al. Three-dimensional modelling of underwater noise produced by a bulk carrier vessel and estimation of its environmental impact
CN110632650B (zh) 全波形反演数据的分时窗提取方法及装置
CN114442169B (zh) 地震资料中近源信号压制方法及装置
CN112764107B (zh) 地表高速覆盖地区的地震资料静校正方法及装置
US11391858B2 (en) Marine survey source firing control
CN112748458B (zh) 用于现场精细变观的在线实时面元获取方法及装置
CN112285782B (zh) 近地表地震波吸收衰减调查方法及装置
US20240302512A1 (en) Systems and methods for an heads up display for use with an underwater handheld computer
CN112114364B (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