[go: up one dir, main page]

CN111929730B - 小尺度地质异常体检测方法及装置 - Google Patents

小尺度地质异常体检测方法及装置 Download PDF

Info

Publication number
CN111929730B
CN111929730B CN202010846124.7A CN202010846124A CN111929730B CN 111929730 B CN111929730 B CN 111929730B CN 202010846124 A CN202010846124 A CN 202010846124A CN 111929730 B CN111929730 B CN 111929730B
Authority
CN
China
Prior art keywords
point
diffraction
diffracted wave
gather
shot
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
CN202010846124.7A
Other languages
English (en)
Other versions
CN111929730A (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 University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN202010846124.7A priority Critical patent/CN111929730B/zh
Publication of CN111929730A publication Critical patent/CN111929730A/zh
Application granted granted Critical
Publication of CN111929730B publication Critical patent/CN111929730B/zh
Priority to US17/506,577 priority patent/US12242007B2/en
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/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/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
    • 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

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

小尺度地质异常体检测方法及装置
技术领域
本发明涉及地震勘探高分辨率成像领域，尤其是涉及一种小尺度地质异常体检测方法及装置。
背景技术
小尺度地质异常体，如断层、陷落柱、裂缝等，往往和矿产资源分布具有紧密的关系，对非均匀不连续地质异常体的精确定位，可有效提高勘探成功率，降低成本，避免可能发生的地质灾害，规避风险。绕射波是小尺度地质异常体的地震响应，包含了小尺度地质异常体的构造信息，可以用来精确定位非均匀不连续地质异常体，提供更强的地下空间照明度。但相对反射波而言，绕射波在传播过程中衰减较快，能量较弱，容易被强能量的反射波掩盖。因此，需要将绕射波从反射波中分离，以增强绕射响应，从而进行绕射波高精度成像，对小尺度地质体精确定位。
现有方法中，一般在共偏移距域或叠后域进行绕射波分离，破坏了波场的一致性，且需要利用偏移算法单独进行绕射波偏移过程，处理流程较为繁琐。该方法基于绕射波传播的几何光学特点，由单炮数据构建共绕射点道集数据，有利于保持波场的一致性和完整性特点，利用同一绕射点不同炮集绕射波时差为常数的特点进行绕射波速度分析，同时完成了绕射波成像过程，不需要再进行波场偏移，处理流程方便快捷。
发明内容
有鉴于此，本发明的目的在于提供一种小尺度地质异常体检测方法及装置，以改善上述技术问题。
第一方面，本发明实施例提供了一种小尺度地质异常体检测方法，所述方法包括以下步骤：
获取在待处理区域进行采集得到的绕射波炮集数据，并在所述绕射波炮集数据中确定与中心点位置之间的距离为预定距离的目标单炮数据；
计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离；
基于所述中心点对应的第一水平距离和第二水平距离，构建所述中心点的共绕射点道集；
利用绕射波同向轴走时校正算法，对所述中心点的共绕射点道集进行处理，得到绕射波成像剖面。
结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，利用绕射波同向轴走时校正算法，对所述中心点的共绕射点道集进行处理，得到绕射波成像剖面包括：
利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集；
对所属于相同中心点的拉平校正之后的共绕射点道集进行叠加，得到所述绕射波成像剖面。
结合第一方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所述中心点的数量为多个，且每个所述中心点对应一个或多个所述炮点和检波点；
计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离包括：
计算每个所述中心点所对应的目标单炮数据中每个炮点和每个所述中心点之间的第一水平距离；
计算每个所述中心点和其所对应的一个或多个所述检波点之间的第二水平距离。
结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第三种可能的实施方式，其中，利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集包括：
获取真实速度和校正速度下的绕射波同向轴走时公式，并根据所述真实速度和所述校正速度下的绕射波同向轴走时公式，获取绕射波同向轴走时校正公式；
根据所述绕射波同向轴走时校正公式，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集。
结合第一方面的第三种可能的实施方式，本发明实施例提供了第一方面的第四种可能的实施方式，其中，根据所述真实速度和所述校正速度下的绕射波同向轴走时公式，获取绕射波同向轴走时校正公式包括：
将所述真实速度下的绕射波同向轴走时公式减去所述校正速度下的绕射波同向轴走时公式，得到所述绕射波同向轴走时校正公式；
所述绕射波同向轴走时校正公式表示为：
Figure BDA0002641576750000031
其中，t0为在每个网格的中心点位置处的垂直双程走时；Sd为第一水平距离，且Rd为第二水平距离；V0为真实速度，且V1为校正速度。
第二方面，本发明实施例还提供一种小尺度地质异常体检测装置，所述装置包括：
数据获取模块，用于获取在待处理区域进行采集得到的绕射波炮集数据，并在所述绕射波炮集数据中确定与中心点位置之间的距离为预定距离的目标单炮数据；
距离计算模块，用于计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离；
道集构建模块，用于基于所述中心点对应的第一水平距离和第二水平距离，构建所述中心点的共绕射点道集；
道集处理模块，用于利用绕射波同向轴走时校正算法，对所述中心点的共绕射点道集进行处理，得到绕射波成像剖面。
第三方面，本发明实施例还提供一种服务器，所述服务器包括：处理器和存储器，所述存储器存储有能够被所述处理器执行的计算机可执行指令，所述处理器执行所述计算机可执行指令以实现上文所述的方法。
结合第二方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所述道集处理模块用于：
利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集。
对所属于相同中心点的拉平校正之后的共绕射点道集进行叠加，得到所述绕射波成像剖面。
结合第二方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所述中心点的数量为多个，且每个所述中心点对应一个或多个所述炮点和检波点；
所述距离计算模块用于：
计算每个所述中心点所对应的目标单炮数据中每个炮点和每个所述中心点之间的第一水平距离；
计算每个所述中心点和其所对应的一个或多个所述检波点之间的第二水平距离。
第四方面，本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机可执行指令，所述计算机可执行指令在被处理器调用和执行时，所述计算机可执行指令促使处理器实现上文所述的方法。
本发明实施例带来了以下有益效果：本发明实施例提供的一种小尺度地质异常体检测方法及装置，获取绕射波炮集数据并确定与中心点位置之间的距离为预定距离的目标单炮数据；计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离；利用绕射波同向轴走时校正算法对中心点的共绕射点道集进行处理得到绕射波成像剖面。本发明实施例的小尺度地质异常体检测方法及装置基于绕射波传播的几何光学特点，由绕射波炮集数据构建共绕射点道集数据，并利用绕射波对于速度的敏感性对绕射波进行直接速度分析并达到成像目的，定位地下空间异常构造以降低事故风险。
本发明的其他特征和优点将在随后的说明书中阐述，并且部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构中实现和获得。
为使本发明的上述目的、特征和优点能够更加明显易懂，下文特举优选实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明的具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图进行简单的介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种小尺度地质异常体检测方法的流程图；
图2为本发明实施例提供的另一种小尺度地质异常体检测方法的流程图；
图3为本发明实施例提供的第一水平距离和第二水平距离的示意图；
图4本发明实施例提供的一种小尺度地质异常体检测装置的结构框图；
图5为本发明实施例提供的一种服务器的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前，在对绕射波进行成像的过程中，大多数的现存方法是在共偏移距域或叠后域对绕射波进行分离，这不仅破坏了绕射波场的一致性，而且需要利用偏移算法单独进行绕射波偏移过程，这使得成像的处理流程较为繁琐。基于此，本发明实施例提供了一种小尺度地质异常体检测方法及装置，以缓解上述问题。
为了便于对本实施例进行理解，首先对本发明实施例所公开的一种小尺度地质异常体检测方法进行详细介绍。
在一种可能的实施方式中，本发明提供了一种小尺度地质异常体检测方法。如图1所示为本发明实施例提供的一种小尺度地质异常体检测方法的流程图，该方法包括以下步骤：
步骤S102：获取在待处理区域进行采集得到的绕射波炮集数据，并在所述绕射波炮集数据中确定与中心点位置之间的距离为预定距离的目标单炮数据。
步骤S104：计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离。
步骤S106：基于所述中心点对应的第一水平距离和第二水平距离，构建所述中心点的共绕射点道集。
步骤S108：利用绕射波同向轴走时校正算法，对所述中心点的共绕射点道集进行处理，得到绕射波成像剖面。
本发明实施例带来了以下有益效果：本发明实施例通过一种小尺度地质异常体检测方法，获取绕射波炮集数据并确定与中心点位置之间的距离为预定距离的目标单炮数据；计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离；利用绕射波同向轴走时校正算法对中心点的共绕射点道集进行处理得到绕射波成像剖面。本发明实施例的小尺度地质异常体检测方法及装置基于绕射波传播的几何光学特点，由绕射波炮集数据构建共绕射点道集数据，并利用绕射波对于速度的敏感性对绕射波进行直接速度分析并达到成像目的，定位地下空间异常构造以降低事故风险果。
在实际使用时，为了对图1中的方法的过程进行更加详细的描述，本发明实施例在图2中示出了本发明实施例提供的另一种小尺度地质异常体检测方法的流程图，该方法包括以下步骤：
步骤S202：获取在待处理区域进行采集得到的绕射波炮集数据，并在所述绕射波炮集数据中确定与中心点位置之间的距离为预定距离的目标单炮数据。
需要进一步说明的是，按照目标层位深度，最大可选取炮点和该中心点距离与目标层位深度相同内的目标炮点数据。
步骤S204：计算每个所述中心点所对应的目标单炮数据中每个炮点和每个所述中心点之间的第一水平距离。
步骤S206：计算每个所述中心点和其所对应的一个或多个所述检波点之间的第二水平距离。
其中，为了对第一水平距离和第二水平距离进行更为直观的呈现，本发明实施例在图3中对两者进行了描述，其中，图3为本发明实施例提供的第一水平距离和第二水平距离的示意图。
其中，O点为中心点，炮点S为该中心点的目标单炮数据，R1和R2为该炮点的检波点，D为地下绕射点，其垂直投影为O点，计算炮点与中心点的第一水平距离Sd、以及中心点和两个检波点第二水平距离Rd。
步骤S208：基于所述中心点对应的第一水平距离和第二水平距离，构建所述中心点的共绕射点道集。
步骤S210：利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集。
具体地，利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集包括：
获取真实速度和校正速度下的绕射波同向轴走时公式，并根据所述真实速度和所述校正速度下的绕射波同向轴走时公式，获取绕射波同向轴走时校正公式；
根据所述绕射波同向轴走时校正公式，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集。
更进一步地，根据所述真实速度和所述校正速度下的绕射波同向轴走时公式，获取绕射波同向轴走时校正公式包括：
将所述真实速度下的绕射波同向轴走时公式减去所述校正速度下的绕射波同向轴走时公式，得到所述绕射波同向轴走时校正公式；
所述绕射波同向轴走时校正公式表示为：
Figure BDA0002641576750000091
其中，t0为在每个网格的中心点位置处的垂直双程走时；Sd为第一水平距离，且Rd为第二水平距离；V 0为真实速度，且V1为校正速度。
步骤S212：对所属于相同中心点的拉平校正之后的共绕射点道集进行叠加，得到所述绕射波成像剖面。
综上所述，本发明的小尺度地质异常体检测方法及装置，获取绕射波炮集数据并确定与中心点位置之间的距离为预定距离的目标单炮数据；计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离；利用绕射波同向轴走时校正算法对中心点的共绕射点道集进行处理得到绕射波成像剖面。本发明实施例的小尺度地质异常体检测方法及装置基于绕射波传播的几何光学特点，由绕射波炮集数据构建共绕射点道集数据，并利用绕射波对于速度的敏感性对绕射波进行直接速度分析并达到成像目的，定位地下空间异常构造以降低事故风险。
在另一种可能的实施方式中，对应于上述实施方式提供的小尺度地质异常体检测方法，本发明实施例还提供了一种小尺度地质异常体检测装置，图4本发明实施例提供的一种小尺度地质异常体检测装置的结构框图。如图4所示，该装置包括：
数据获取模块401，用于获取在待处理区域进行采集得到的绕射波炮集数据，并在所述绕射波炮集数据中确定与中心点位置之间的距离为预定距离的目标单炮数据；
距离计算模块402，用于计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离；
道集构建模块403，用于基于所述中心点对应的第一水平距离和第二水平距离，构建所述中心点的共绕射点道集；
道集处理模块404，用于利用绕射波同向轴走时校正算法，对所述中心点的共绕射点道集进行处理，得到绕射波成像剖面。
在实际使用时，所述道集处理模块404用于：
利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集；
对所属于相同中心点的拉平校正之后的共绕射点道集进行叠加，得到所述绕射波成像剖面。
在实际使用时，所述中心点的数量为多个，且每个所述中心点对应一个或多个所述炮点和检波点；
所述距离计算模块402用于：
计算每个所述中心点所对应的目标单炮数据中每个炮点和每个所述中心点之间的第一水平距离；
计算每个所述中心点和其所对应的一个或多个所述检波点之间的第二水平距离。
在又一种可能的实施方式中，本发明实施例还提供了一种服务器，图5示出了本发明实施例提供的一种服务器的结构示意图，参见图5，该服务器包括：处理器500、存储器501、总线502和通信接口503，该处理器500、存储器501、通信接口503和通过总线502连接；处理器500用于执行存储器501中存储的可执行模块，例如计算机程序。
其中，存储器501存储有能够被处理器500执行的计算机可执行指令，处理器500执行计算机可执行指令以实现上文所述的方法。
进一步地，存储器501可能包含高速随机存取存储器(RAM，Random AccessMemory)，也可能还包括非不稳定的存储器(non-volatile memory)，例如至少一个磁盘存储器。通过至少一个通信接口503(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。
总线502可以是ISA总线、PCI总线或EISA总线等。该总线可以分为地址总线、数据总线、控制总线等。为便于表示，图4中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
其中，存储器501用于存储程序，处理器500在接收到程序执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的小尺度地质异常体检测方法可以应用于处理器500中，或者由处理器500实现。
此外，处理器500可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器500中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器500可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(DigitalSignal Processing，简称DSP)、专用集成电路(Application Specific IntegratedCircuit，简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器501，处理器500读取存储器501中的信息，结合其硬件完成上述方法的步骤。
在又一种可能的实施方式中，本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机可执行指令，所述计算机可执行指令在被处理器调用和执行时，计算机可执行指令促使处理器实现上文所述的方法。
本发明实施例提供的小尺度地质异常体检测装置，与上述实施例提供的小尺度地质异常体检测方法具有相同的技术特征，所以也能解决相同的技术问题，达到相同的技术效果。
本发明实施例所提供的小尺度地质异常体检测方法及装置的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的装置的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
另外，在本发明实施例的描述中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，ReaD-Only Memory)、随机存取存储器(RAM，RanDom Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
最后应说明的是：以上实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

Claims (7)

1.一种小尺度地质异常体检测方法，其特征在于，所述方法包括以下步骤：
获取在待处理区域进行采集得到的绕射波炮集数据，并在所述绕射波炮集数据中确定与中心点位置之间的距离为预定距离的目标单炮数据；
计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离；
基于所述中心点对应的第一水平距离和第二水平距离，构建所述中心点的共绕射点道集；
利用绕射波同向轴走时校正算法，对所述中心点的共绕射点道集进行处理，得到绕射波成像剖面；
其中，利用绕射波同向轴走时校正算法，对所述中心点的共绕射点道集进行处理，得到绕射波成像剖面包括：
利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集；
对所属于相同中心点的拉平校正之后的共绕射点道集进行叠加，得到所述绕射波成像剖面；
其中，利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集包括：
获取真实速度和校正速度下的绕射波同向轴走时公式，并根据所述真实速度和所述校正速度下的绕射波同向轴走时公式，获取绕射波同向轴走时校正公式；
根据所述绕射波同向轴走时校正公式，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集。
2.根据权利要求1所述的方法，其特征在于，所述中心点的数量为多个，且每个所述中心点对应一个或多个所述炮点和检波点；
计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离包括：
计算每个所述中心点所对应的目标单炮数据中每个炮点和每个所述中心点之间的第一水平距离；
计算每个所述中心点和其所对应的一个或多个所述检波点之间的第二水平距离。
3.根据权利要求1所述的方法，其特征在于，根据所述真实速度和所述校正速度下的绕射波同向轴走时公式，获取绕射波同向轴走时校正公式包括：
将所述真实速度下的绕射波同向轴走时公式减去所述校正速度下的绕射波同向轴走时公式，得到所述绕射波同向轴走时校正公式；
所述绕射波同向轴走时校正公式表示为：
Figure FDA0002966454780000021
其中，t0为在每个网格的中心点位置处的垂直双程走时；Sd为第一水平距离，且Rd为第二水平距离；V0为真实速度，且V1为校正速度。
4.一种小尺度地质异常体检测装置，其特征在于，所述装置包括：
数据获取模块，用于获取在待处理区域进行采集得到的绕射波炮集数据，并在所述绕射波炮集数据中确定与中心点位置之间的距离为预定距离的目标单炮数据；
距离计算模块，用于计算所述目标单炮数据中的每个炮点与所述中心点的第一水平距离，并计算每个所述炮点对应的检波点与所述中心点的第二水平距离；
道集构建模块，用于基于所述中心点对应的第一水平距离和第二水平距离，构建所述中心点的共绕射点道集；
道集处理模块，用于利用绕射波同向轴走时校正算法，对所述中心点的共绕射点道集进行处理，得到绕射波成像剖面；
其中，所述道集处理模块用于：
利用绕射波同向轴走时校正算法，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集；
对所属于相同中心点的拉平校正之后的共绕射点道集进行叠加，得到所述绕射波成像剖面；
其中，所述道集处理模块还用于：
获取真实速度和校正速度下的绕射波同向轴走时公式，并根据所述真实速度和所述校正速度下的绕射波同向轴走时公式，获取绕射波同向轴走时校正公式；
根据所述绕射波同向轴走时校正公式，在所述中心点的共绕射点道集中对绕射波在同向轴方向进行拉平校正处理，得到拉平校正之后的共绕射点道集。
5.根据权利要求4所述的装置，其特征在于，所述中心点的数量为多个，且每个所述中心点对应一个或多个所述炮点和检波点；
所述距离计算模块用于：
计算每个所述中心点所对应的目标单炮数据中每个炮点和每个所述中心点之间的第一水平距离；
计算每个所述中心点和其所对应的一个或多个所述检波点之间的第二水平距离。
6.一种服务器，其特征在于，包括处理器和存储器，所述存储器存储有能够被所述处理器执行的计算机可执行指令，所述处理器执行所述计算机可执行指令以实现权利要求1至3任一项所述的方法。
7.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机可执行指令，所述计算机可执行指令在被处理器调用和执行时，所述计算机可执行指令促使处理器实现权利要求1至3任一项所述的方法。
CN202010846124.7A 2020-08-20 2020-08-20 小尺度地质异常体检测方法及装置 Active CN111929730B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010846124.7A CN111929730B (zh) 2020-08-20 2020-08-20 小尺度地质异常体检测方法及装置
US17/506,577 US12242007B2 (en) 2020-08-20 2021-10-20 Small-scale geological anomalous body detection method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010846124.7A CN111929730B (zh) 2020-08-20 2020-08-20 小尺度地质异常体检测方法及装置

Publications (2)

Publication Number Publication Date
CN111929730A CN111929730A (zh) 2020-11-13
CN111929730B true CN111929730B (zh) 2021-04-16

Family

ID=73304873

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010846124.7A Active CN111929730B (zh) 2020-08-20 2020-08-20 小尺度地质异常体检测方法及装置

Country Status (2)

Country Link
US (1) US12242007B2 (zh)
CN (1) CN111929730B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112904424B (zh) * 2021-01-21 2021-10-15 中国科学院地理科学与资源研究所 一种异常体位置的确定方法、装置及可读存储介质

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5293352A (en) * 1993-01-07 1994-03-08 Western Atlas International, Inc. Method for removing noise due to near surface scatterers
US6151556A (en) * 1999-06-18 2000-11-21 Mobil Oil Corporation Method and apparatus for doppler smear correction in marine seismology measurements
US8209125B2 (en) * 2007-03-12 2012-06-26 Geomage (2003) Ltd. Method for identifying and analyzing faults/fractures using reflected and diffracted waves
WO2013012610A2 (en) * 2011-07-19 2013-01-24 Halliburton Energy Services, Inc. System and method for moment tensor migration imaging
US11391857B2 (en) * 2013-12-30 2022-07-19 Pgs Geophysical As Methods and systems for attenuating residual acoustic energy in seismic data
CN105301649B (zh) * 2015-06-03 2017-07-28 中国地质大学（北京） 转换波偏移速度建模装置
CN105607121B (zh) * 2016-02-02 2016-12-21 中国矿业大学(北京) 一种煤炭陷落柱识别方法和装置
CN106772583B (zh) * 2017-01-10 2018-09-04 中国科学院地质与地球物理研究所 一种地震绕射波分离方法和装置
CN107918147B (zh) * 2017-11-20 2018-12-14 中国矿业大学(北京) 绕射波成像方法和装置
CN108051858B (zh) * 2017-12-13 2018-11-13 中国矿业大学(北京) 优化的地震勘探方法及装置
CN110031900B (zh) * 2019-05-27 2020-05-01 中国矿业大学(北京) 绕射波成像方法、装置、电子设备和存储介质
US11681043B2 (en) * 2019-09-03 2023-06-20 Saudi Arabian Oil Company Diffraction imaging using pseudo dip-angle gather
US11313988B2 (en) * 2019-12-13 2022-04-26 Saudi Arabian Oil Company Identifying geologic features in a subterranean formation using seismic diffraction imaging
US11402529B2 (en) * 2020-01-09 2022-08-02 Saudi Arabian Oil Company Identifying geologic features in a subterranean formation using seismic diffraction and refraction imaging

Also Published As

Publication number Publication date
US12242007B2 (en) 2025-03-04
CN111929730A (zh) 2020-11-13
US20220196865A1 (en) 2022-06-23

Similar Documents

Publication Publication Date Title
CN110084116B (zh) 路面检测方法、装置、计算机设备和存储介质
CN111238352B (zh) 墙面虚拟靠尺检测方法、系统、设备和存储介质
WO2019126972A1 (zh) 基于InSAR的形变信息提取方法、终端及存储介质
US12352914B2 (en) Stereoscopic identification method and apparatus for disturbance stress evolution process of underground cave surrounding rock
CN118823095B (zh) 基于视觉模型和局部单元法的结构构件截面尺寸检测方法
CN111929730B (zh) 小尺度地质异常体检测方法及装置
CN109490951B (zh) 绕射波成像方法、装置和电子设备
CN110260918A (zh) 建筑监测方法、装置、电子设备及计算机可读存储介质
Zhao et al. Automatic identification of concealed dangerous rock blocks on high-steep slopes considering finite-sized discontinuity intersections
CN109405740A (zh) 一种基于收割机升运器的谷物体积计算方法及装置
CN116434072A (zh) 基于多源数据的地质灾害早期识别方法、装置
CN115903043A (zh) 一种绕射波分离方法和装置
CN119669706A (zh) 多尺度黄土地基稳定性预测方法、系统、存储介质及设备
CN111210120B (zh) 极震区和影响场方向的确定方法、装置、设备及存储介质
CN113608261B (zh) 绕射波成像方法、装置和电子设备
WO2019056536A1 (zh) 井眼截面形状测量方法、装置和终端
CN112379420A (zh) 高精度弯曲测线叠前时间域成像方法及装置
CN112465983A (zh) 一种高程地表模型点线状障碍物修正方法和装置
CN113267815B (zh) 重复断棱数据的滤除方法及装置
CN117741744A (zh) 地震预测方法、装置、电子设备及存储介质
CN115951336A (zh) 确定激光雷达误差的方法、装置、设备及存储介质
CN118068410A (zh) 一种散射波的偏移成像方法、装置、设备及存储介质
CN120802208B (zh) 激光点云平面检测与距离测量方法、装置、设备及介质
CN113031073B (zh) 叠前深度偏移速度确定方法、装置及设备
CN119535588B (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