CN112748466A - 一种基于菲涅尔体的旅行时场数据处理方法及装置 - Google Patents

一种基于菲涅尔体的旅行时场数据处理方法及装置 Download PDF

Info

Publication number
CN112748466A
CN112748466A CN201911043097.3A CN201911043097A CN112748466A CN 112748466 A CN112748466 A CN 112748466A CN 201911043097 A CN201911043097 A CN 201911043097A CN 112748466 A CN112748466 A CN 112748466A
Authority
CN
China
Prior art keywords
field data
time field
travel
shot point
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.)
Granted
Application number
CN201911043097.3A
Other languages
English (en)
Other versions
CN112748466B (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 CN201911043097.3A priority Critical patent/CN112748466B/zh
Publication of CN112748466A publication Critical patent/CN112748466A/zh
Application granted granted Critical
Publication of CN112748466B publication Critical patent/CN112748466B/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
    • 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/282—Application of seismic models, synthetic seismograms
    • 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/303—Analysis for determining velocity profiles or travel times

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)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种基于菲涅尔体的旅行时场数据处理方法及装置，所述方法包括：根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据依次采用第一数据限定规则、第二数据限定规则和第三数据限定规则进行数据限定，最终获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。所述装置用于执行上述方法。本发明实施例提供的基于菲涅尔体的旅行时场数据处理方法及装置，提高了旅行时场数据的处理效率。

Description

一种基于菲涅尔体的旅行时场数据处理方法及装置
技术领域
本发明涉及地质勘探技术领域，具体涉及一种基于菲涅尔体的旅行时场数据处理方法及装置。
背景技术
目前，菲涅尔体层析方法在近地表速度反演和二氧化碳储藏的时移地震监测中的应用展现出了其巨大的潜力，是近地表建模非常重要的技术之一。
菲涅尔体射线追踪技术，不仅考虑了最优的高频射线路径上速度对地震波走势的影响，还结合频率信息考虑了在第一菲涅尔带范围内位置上的影响，更符合地震波实际传播的规律。并且，由于能弥补高频射线数量稀少、覆盖不均匀导致的层析反演系数矩阵稀疏，从而解决求解过程不稳定的问题，得到更加精确的反演结果，提高了反演稳定性。但是由于菲涅尔体射线追踪结果信息丰富，导致菲涅尔体的计算复杂程度高，耗时巨大，特别是在三维工区的应用过程中，由于三维初至信息陡增，导致菲涅尔体射线追踪计算量成几何数量级增长，计算效率十分低下。
因此，如何提出一种基于菲涅尔体的旅行时场数据处理方法，能够减少旅行时场数据处理量，提高旅行时场数据处理效率成为本领域需要解决的重要课题。
发明内容
针对现有技术中的问题，本发明实施例提供一种基于菲涅尔体的旅行时场数据处理方法及装置。
一方面，本发明提出一种基于菲涅尔体的旅行时场数据处理方法，包括：
根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；
根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；
根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；
根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
另一方面，本发明提供一种基于菲涅尔体的旅行时场数据处理装置，包括：
第一获得单元，用于根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；
第二获得单元，用于根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；
第三获得单元，用于根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；
第四获得单元，用于根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
再一方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一实施例所述基于菲涅尔体的旅行时场数据处理方法的步骤。
又一方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一实施例所述基于菲涅尔体的旅行时场数据处理方法的步骤。
本发明实施例提供的基于菲涅尔体的旅行时场数据处理方法及装置，能够根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据，接着根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据，再根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据，然后根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据，缩小了旅行时场数据的范围，减少了参与菲涅尔体计算的数量，提高了旅行时场数据的处理效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1是本发明一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图。
图2a是本发明另一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图。
图2b是本发明还一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图。
图3是本发明又一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图。
图4是本发明再一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图。
图5是本发明一实施例提供的基于菲涅尔体的旅行时场数据处理装置的结构示意图。
图6是本发明另一实施例提供的基于菲涅尔体的旅行时场数据处理装置的结构示意图。
图7是本发明又一实施例提供的基于菲涅尔体的旅行时场数据处理装置的结构示意图。
图8是本发明再一实施例提供的基于菲涅尔体的旅行时场数据处理装置的结构示意图。
图9是本发明一实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
图1是本发明一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图，如图1所示，本发明实施例提供的基于菲涅尔体的旅行时场数据处理方法，包括：
S101、根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；
具体地，服务器在获得目标区域的旅行时初至数据以及初始速度模型之后，可以根据所述旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据，例如可以采用Fast Marching Method(FMA)算法根据所述旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据。其中，所述旅行时初至数据包括每个炮点的坐标和高程，以及每个炮点对应的检波点的坐标和高程。所述初始速度模型是利用所述目标区域的旅行时初至数据中所有炮点和检波点的坐标与高程插值形成的，对于二维的初始速度模型在水平方向x和垂直方向z上具有一定网格数量，对于三维的初始速度模型在x和z以及垂直于xz平面的方向y上具有一定网格数量。可理解的是，所述目标区域的旅行时初至数据是所述目标区域的原始地震数据进行处理后获得的。其中，本发明实施例提供的超参数的自动确定方法的执行主体包括但不限于服务器。
例如，对于所述目标区域，为了获得二维的初始速度模型，可以在所述目标区域的地表沿一条直线设置第一预设数量的炮点，每个炮点对应第二预设数量的检波点，所述第二预设数量的炮点均匀分布在对应的炮点的两侧，利用设置在每个检波点的检波器采集原始地震数据，经过相关软件，例如Geoeast软件、tomodel软件、omega等行业内软件的处理，可以获得目标区域的二维旅行时初至数据，再对上述二维旅行时初至数据进行处理，可以二维的初始速度模型。所述服务器再根据上述二维旅行时初至数据和二维的初始速度模型，可以获得所述目标区域的每个炮点和每个检波点的原始旅行时场数据。其中，所述第一预设数量和所述第二预设数量根据实际需要进行设置，本发明实施例不做限定。
例如，对于所述目标区域，为了获得三维的初始速度模型，可以在所述目标区域的地表设置第三预设数量的炮点，每个炮点有多条接收线并且每个炮点位于多条接收线的交点，每条接收线上设置第四预设数量的检波点。利用设置在每个检波点的检波器采集原始地震数据，原始地震数据经过相关软件的处理，可以获得目标区域的三维旅行时初至数据，再对上述三维旅行时初至数据进行处理，可以三维的初始速度模型。所述服务器再根据上述三维旅行时初至数据和三维的初始速度模型，可以获得所述目标区域的每个炮点和每个检波点的原始旅行时场数据。其中，所述第三预设数量和所述第四预设数量根据实际需要进行设置，本发明实施例不做限定。
S102、根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；
具体地，所述服务器在获得每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据之后，根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，可以获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据。其中，所述第一数据限定规则可以是对每个炮点和每个检波点在偏移方向上、高度方向上或者深度方向上的原始旅行时场数据的限定。如果所述原始旅行时场数据为二维数据，那么所述偏移方向是指水平方向，所述高度方向是指垂直于水平方向向上的方向，所述深度方向是指垂直于水平方向向下的方向。如果所述原始旅行时场数据为三维数据，那么所述偏移方向是指水平面上方向，所述高度方向是指垂直于水平面向上的方向，所述深度方向是指垂直于水平面向下的方向。
例如，所述服务器可以根据每个炮点以及对应的各个所述检波点的网格坐标，获得每个炮点相对于检波点的最大偏移距离，然后根据每个炮点相对于检波点的最大偏移距离对每个炮点的原始旅行时场数据进行限定，从而获得每个炮点的第一限定旅行时场数据。所述服务器可以根据每个炮点以及对应的各个所述检波点，获得每个检波点对应的各个所述炮点，然后根据每个检波点以及对应的各个所述炮点的网格坐标，获得每个检波点相对于炮点的最大偏移距离，再根据每个检波点相对于炮点的最大偏移距离对每个检波点的原始旅行时场数据进行限定，获得每个检波点的第一限定旅行时场数据。
S103、根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；
具体地，所述服务器在获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据之后，可以根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据。其中，所述第二数据限定规则是对每个炮点和每个检波点在偏移方向上、高度方向上或者深度方向上的第一限定旅行时场数据的限定。可理解的是，所述第二数据限定规则与所述第一数据限定规则不同，即如果所述第一数据限定规则是在偏移方向上的数据限定，那么所述第二数据限定规则是在高度方向上或者深度方向上的数据限定，如果所述第一数据限定规则是在高度方向上的数据限定，那么所述第二数据限定规则是在偏移方向上或者深度方向上的数据限定；如果所述第一数据限定规则是在深度方向上的数据限定，那么所述第二数据限定规则是在高度方向上或者深度方向上的数据限定。
例如，所述服务器根据各个所述炮点以及各个所述检波点的在垂直方向上的网格坐标，获得高程限定边界，然后根据所述高程限定边界对每个炮点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据，并根据所述高程限定条件对每个检波点的第一限定旅行时场数据进行限定，获得每个检波点的第二限定旅行时场数据。
S104、根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
具体地，所述服务器在获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据之后，可以根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。其中，所述第三数据限定规则是对每个炮点和每个检波点在偏移方向上、高度方向上或者深度方向上的第一限定旅行时场数据的限定。可理解的是，所述第三数据限定规则与所述第一数据限定规则和所述第二数据限定规则都不同，即如果所述第三数据限定规则是在深度方向上的数据限定，那么所述第一数据限定规则和所述第二数据限定规则是在偏移方向上和高度方向上的数据限定；如果所述第三数据限定规则是在偏移方向上的数据限定，那么所述第一数据限定规则和所述第二数据限定规则是在深度方向上和高度方向上的数据限定；如果所述第三数据限定规则是在高度方向上的数据限定，那么所述第一数据限定规则和所述第二数据限定规则是在偏移方向上和深度方向上的限定。
例如，所述服务器可以根据每个炮点的第二限定旅行时场数据以及每个炮点的最大偏移距离对应的检波点的第二限定旅行时场数据，获得每个炮点的最大偏移距离炮检对的旅行时场数据，接着根据每个炮点的最大偏移距离炮检对的旅行时场数据，获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息，然后根据每个炮点的最大偏移距离炮检对的菲涅尔体路径信息，获得每个炮点的最大菲涅尔体深度，再根据每个炮点的最大菲涅尔体深度对每个炮点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据，并根据每个炮点的最大菲涅尔体深度对每个炮点对应的各个所述检波点的第二限定旅行时场数据进行限定，获得每个炮点对应的各个所述检波点的最终限定旅行时场数据。
本发明实施例提供的基于菲涅尔体的旅行时场数据处理方法，能够根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据，接着根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据，再根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据，然后根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据，缩小了旅行时场数据的范围，减少了参与菲涅尔体计算的数量，提高了旅行时场数据的处理效率。
图2a是本发明另一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图，图2b是本发明还一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图，如图2a和图2b所示，所述根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据包括：
S1021、根据每个炮点以及对应的各个所述检波点的网格坐标，获得每个炮点相对于检波点的最大偏移距离；
具体地，对于每个炮点的原始旅行时场数据，每个炮点以及对应的各个所述检波点都是以网格坐标表示的。所述服务器根据每个炮点以及对应的各个所述检波点的网格坐标，可以获得每个炮点与对应的每个检波点之间的偏移距离，上述偏移距离可以用网格数表示。所述服务器对比每个炮点与对应的各个所述检波点之间的偏移距离，从中获得最大的偏移距离，作为每个炮点相对于检波点的最大偏移距离，对于每个炮点都会有一个所述最大偏移距离。可理解的是，如果所述原始旅行时场数据为二维数据，那么炮点与对应的检波点的偏移距离是指炮点和对应的检波点在水平方向上的距离；如果所述原始旅行时场数据为三维数据，那么炮点与对应的检波点的偏移距离是指炮点和对应的检波点之间的直线距离在水平面上的投影。
S1022、根据每个炮点相对于检波点的最大偏移距离对每个炮点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据；
具体地，所述服务器在获得每个炮点相对于检波点的最大偏移距离之后，根据每个炮点相对于检波点的最大偏移距离对每个炮点的原始旅行时场数据进行限定，即保留每个炮点的原始旅行时场数据中与每个炮点的偏移距离在所述最大偏移距离范围内的数据，舍弃每个炮点的原始旅行时场数据中与每个炮点的偏移距离大于所述最大偏移距离的数据，获得每个炮点的第一限定旅行时场数据。
例如，如果所述炮点的原始旅行时场数据为二维数据，所述炮点的网格坐标为(x1，z1)，x1为水平方向坐标，z1为垂直方向坐标，所述炮点相对于检波点的最大偏移距离为d1，那么保留所述炮点的原始旅行时场数据中水平方向的坐标大于等于x1-d1且小于等于x1+d1的数据，作为所述炮点的第一限定旅行时场数据。
例如，如果所述炮点的原始旅行时场数据为三维数据，所述炮点的网格坐标为(x2，y2，z2)，x2为水平面上的第一方向上的坐标，y2为水平面上的第一方向上的坐标，第一方向与第二方向垂直，z2为垂直与水平面方向上的坐标，所述炮点相对于检波点的最大偏移距离为d2，那么保留所述炮点的原始旅行时场数据中第一方向的坐标大于等于x2-d2且小于等于x2+d2，第二方向的坐标大于等于y2-d2且小于等于y2+d2的数据，作为所述炮点的第一限定旅行时场数据。
S1023、根据每个炮点以及对应的各个所述检波点，获得每个检波点对应的各个所述炮点；
具体地，由于已知每个炮点以及对应的各个所述检波点，所述服务器可以根据每个检波点，在每个炮点以及对应的各个所述检波点中查找到每个检波点对应的各个炮点。
例如，每个检波点都有对应的编号，所述服务器根据检波点a的编号查找炮点b对应的各个所述检波点的编号，如果所述炮点对应的各个所述检波点的编号中存在与检波点a的编号相同的编号，那么炮点b为检波点a对应的炮点。对每个炮点对应的各个所述检波点的编号，重复上述过程，即可获得检波点a对应的各个所述炮点。
S1024、根据每个检波点以及对应的各个所述炮点的网格坐标，获得每个检波点相对于炮点的最大偏移距离；
具体地，本步骤的具体实现过程，与步骤S1021类似，此处不进行赘述。
S1025、根据每个检波点相对于炮点的最大偏移距离对每个检波点的原始旅行时场数据进行限定，获得每个检波点的第一限定旅行时场数据。
具体地，本步骤的具体实现过程，与步骤S1022类似，此处不进行赘述。
图3是本发明又一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图，如图3所示，在上述各实施例的基础上，进一步地，所述根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据包括：
S1031、根据各个所述炮点以及各个所述检波点的在垂直方向上的网格坐标，获得高程限定边界；
具体地，所述旅行时初至数据中每个炮点的高程转换成了每个炮点在垂直方向上的网格坐标，所述旅行时初至数据中每个检波点的高程转换成了每个检波点在垂直方向上的网格坐标。所述服务器根据各个所述炮点以及各个所述检波点的在垂直方向上的网格坐标，可以拟合出高程限定边界。其中，如果所述旅行时初至数据为二维数据，那么拟合出高程限定边界为一条曲线。如果所述旅行时初至数据为三维数据，那么拟合出高程限定边界为一个曲面。
S1032、根据所述高程限定边界对每个炮点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据，并根据所述高程限定条件对每个检波点的第一限定旅行时场数据进行限定，获得每个检波点的第二限定旅行时场数据。
具体地，所述服务器在获得所述高程限定边界之后，可以根据所述高程限定边界对每个炮点的第一限定旅行时场数据进行限定，即将每个炮点的第一限定旅行时场数据中所述高程限定边界以上的数据舍去，获得每个炮点的第二限定旅行时场数据。所述服务器还可以根据所述高程限定边界对每个检波点的第一限定旅行时场数据进行限定，即将每个检波点的第一限定旅行时场数据中所述高程限定边界以上的数据舍去，获得每个检波点的第二限定旅行时场数据。
图4是本发明再一实施例提供的基于菲涅尔体的旅行时场数据处理方法的流程示意图，如图4所示，在上述各实施例的基础上，进一步地，所述根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据包括：
S1041、根据每个炮点的第二限定旅行时场数据以及每个炮点的最大偏移距离对应的检波点的第二限定旅行时场数据，获得每个炮点的最大偏移距离炮检对的旅行时场数据；
具体地，每个炮点与对应的每个检波点之间存在偏移距离，所述服务器可以从每个炮点与对应的各个所述检波点之间的偏移距离中获得每个炮点相对于检波点的最大偏移距离，即每个炮点的最大偏移距离。所述服务器将每个炮点的第二限定旅行时场数据以及每个炮点的最大偏移距离对应的检波点的第二限定旅行时场数据进行叠加，可以获得每个炮点的最大偏移距离炮检对的旅行时场数据。其中，每个炮点与每个炮点最大偏移距离对应的检波点构成每个炮点的最大偏移距离炮检对。
例如，炮点c与对应的各个所述检波点之间的偏移距离中，炮点c与检波点d之间的偏移距离最大，那么炮点c的最大偏移距离为炮点c与检波点d之间的偏移距离，检波点d为炮点c的最大偏移距离对应的检波点，炮点c与检波点d构成了炮点c的最大偏移距离炮检对。炮点c的第二限定旅行时场数据与检波点d的第二限定旅行时场数据叠加之后，形成了炮点c的最大偏移距离炮检对的旅行时场数据。
S1042、根据每个炮点的最大偏移距离炮检对的旅行时场数据，获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息；
具体地，所述服务器在获得每个炮点的最大偏移距离炮检对的旅行时场数据之后，根据每个炮点的最大偏移距离炮检对的旅行时场数据以及经验公式，可以计算获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息。
例如，所述经验公式为：tL(g,r)+tL(r,s)-tLmin(g,s)≤1/(2f)，其中，g代表地面检波点位置，s代表地面炮点位置，tL(g,r)是指在检波点旅行时场中任意点r到检波点g位置的时间，tL(r,s)是指在炮点旅行时场中r点到炮点s的时间，tLmin(g,s)是炮点旅行时场与检波点旅行时场之和的总时场中最小值得位置，f代表频率，此处是预设的，比如25Hz。上述经验公式的意义就是如果某位置点到炮点的时间与到检波点时间之和与炮检旅行总时场的最小值之差小于等于给定的频率倒数的一半，那么认为该位置点在有效菲涅尔体范围内。
S1043、根据每个炮点的最大偏移距离炮检对的菲涅尔体路径信息，获得每个炮点的最大菲涅尔体深度；
具体地，所述服务器在获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息之后，可以根据每个炮点的最大偏移距离炮检对的菲涅尔体路径信息，计算获得每个炮点的最大菲涅尔体深度。
例如，在获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息之后，就可以获得每个炮点的最大偏移距离炮检对的菲涅尔体范围，上述菲涅尔体范围内各个位置点在深度方向上的坐标的绝对值的最大值，即为每个炮点的最大菲涅尔体深度。
S1044、根据每个炮点的最大菲涅尔体深度对每个炮点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据，并根据每个炮点的最大菲涅尔体深度对每个炮点对应的各个所述检波点的第二限定旅行时场数据进行限定，获得每个炮点对应的各个所述检波点的最终限定旅行时场数据。
具体地，所述服务器在获得每个炮点的最大菲涅尔体深度之后，可以根据每个炮点的最大菲涅尔体深度对每个炮点的第二限定旅行时场数据进行限定，即保留每个炮点的第二限定旅行时场数据中在每个炮点的最大菲涅尔体深度范围内的数据，舍弃每个炮点的第二限定旅行时场数据中超出每个炮点的最大菲涅尔体深度范围的数据，获得每个所述炮点的最终限定旅行时场数据。所述服务器还可以根据每个炮点的最大菲涅尔体深度对每个炮点对应的各个检波点的第二限定旅行时场数据进行限定，即保留每个炮点对应的每个检波点的第二限定旅行时场数据中在每个炮点的最大菲涅尔体深度范围内的数据，舍弃每个炮点对应的每个检波点的第二限定旅行时场数据中超出每个炮点的最大菲涅尔体深度范围的数据，获得每个炮点对应的各个所述检波点的最终限定旅行时场数据。
例如，所述炮点的原始旅行时场数据为二维数据，所述炮点的网格坐标为(x1，z1)，x1为水平方向坐标，z1为垂直方向坐标，所述炮点的最大菲涅尔体深度为h1，那么保留所述炮点的第二限定旅行时场数据中垂直方向的坐标大于等于z1-h1的数据，作为所述炮点的最终限定旅行时场数据。
例如，如果所述炮点的原始旅行时场数据为三维数据，所述炮点的网格坐标为(x2，y2，z2)，x2为水平面上的第一方向上的坐标，y2为水平面上的第一方向上的坐标，第一方向与第二方向垂直，z2为垂直于水平面方向上的坐标，所述炮点的最大菲涅尔体深度为h2，那么保留所述炮点的第二限定旅行时场数据中垂直于水平方向的坐标大于等于z2-h2的数据，作为所述炮点的最终限定旅行时场数据。
图5是本发明一实施例提供的基于菲涅尔体的旅行时场数据处理装置的结构示意图，如图5所示，本发明实施例提供的基于菲涅尔体的旅行时场数据处理装置包括第一获得单元501、第二获得单元502、第三获得单元503和第四获得单元504，其中：
第一获得单元501用于根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；第二获得单元502用于根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；第三获得单元503用于根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；第四获得单元504用于根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
具体地，在获得目标区域的旅行时初至数据以及初始速度模型之后，第一获得单元501可以根据所述旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据，例如可以采用FMA算法根据所述旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据。其中，所述旅行时初至数据包括每个炮点的坐标和高程，以及每个炮点对应的检波点的坐标和高程。所述初始速度模型是利用所述目标区域的旅行时初至数据中所有炮点和检波点的坐标与高程插值形成的，对于二维的初始速度模型在水平方向x和垂直方向z上具有一定网格数量，对于三维的初始速度模型在x和z以及垂直于xz平面的方向y上具有一定网格数量。可理解的是，所述目标区域的旅行时初至数据是所述目标区域的原始地震数据进行处理后获得的。
在获得每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据之后，第二获得单元502根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，可以获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据。其中，所述第一数据限定规则可以是对每个炮点和每个检波点在偏移方向上、高度方向上或者深度方向上的原始旅行时场数据的限定。如果所述原始旅行时场数据为二维数据，那么所述偏移方向是指水平方向，所述高度方向是指垂直于水平方向向上的方向，所述深度方向是指垂直于水平方向向下的方向。如果所述原始旅行时场数据为三维数据，那么所述偏移方向是指水平面上方向，所述高度方向是指垂直于水平面向上的方向，所述深度方向是指垂直于水平面向下的方向。
在获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据之后，第三获得单元503可以根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据。其中，所述第二数据限定规则是对每个炮点和每个检波点在偏移方向上、高度方向上或者深度方向上的第一限定旅行时场数据的限定。可理解的是，所述第二数据限定规则与所述第一数据限定规则不同，即如果所述第一数据限定规则是在偏移方向上的数据限定，那么所述第二数据限定规则是在高度方向上或者深度方向上的数据限定，如果所述第一数据限定规则是在高度方向上的数据限定，那么所述第二数据限定规则是在偏移方向上或者深度方向上的数据限定；如果所述第一数据限定规则是在深度方向上的数据限定，那么所述第二数据限定规则是在高度方向上或者深度方向上的数据限定。
在获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据之后，第四获得单元504可以根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。其中，所述第三数据限定规则是对每个炮点和每个检波点在偏移方向上、高度方向上或者深度方向上的第一限定旅行时场数据的限定。可理解的是，所述第三数据限定规则与所述第一数据限定规则和所述第二数据限定规则都不同，即如果所述第三数据限定规则是在深度方向上的数据限定，那么所述第一数据限定规则和所述第二数据限定规则是在偏移方向上和高度方向上的数据限定；如果所述第三数据限定规则是在偏移方向上的数据限定，那么所述第一数据限定规则和所述第二数据限定规则是在深度方向上和高度方向上的数据限定；如果所述第三数据限定规则是在高度方向上的数据限定，那么所述第一数据限定规则和所述第二数据限定规则是在偏移方向上和深度方向上的限定。
本发明实施例提供的基于菲涅尔体的旅行时场数据处理装置，能够根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据，接着根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据，再根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据，然后根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据，缩小了旅行时场数据的范围，减少了参与菲涅尔体计算的数量，提高了旅行时场数据的处理效率。
图6是本发明另一实施例提供的基于菲涅尔体的旅行时场数据处理装置的结构示意图，如图6所示，在上述各实施例的基础上，进一步地，第二获得单元502包括第一获得子单元5021、第二获得子单元5022、第三获得子单元5023、第四获得子单元5024和第五获得子单元5025，其中：
第一获得子单元5021用于根据每个炮点以及对应的各个所述检波点的网格坐标，获得每个炮点相对于检波点的最大偏移距离；第二获得子单元5022用于根据每个炮点相对于检波点的最大偏移距离对每个炮点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据；第三获得子单元5023用于根据每个炮点以及对应的各个所述检波点，获得每个检波点对应的各个所述炮点；第四获得子单元5024用于根据每个检波点以及对应的各个所述炮点的网格坐标，获得每个检波点相对于炮点的最大偏移距离；第五获得子单元5025用于根据每个检波点相对于炮点的最大偏移距离对每个检波点的原始旅行时场数据进行限定，获得每个检波点的第一限定旅行时场数据。
具体地，对于每个炮点的原始旅行时场数据，每个炮点以及对应的各个所述检波点都是以网格坐标表示的。第一获得子单元5021根据每个炮点以及对应的各个所述检波点的网格坐标，可以获得每个炮点与对应的每个检波点之间的偏移距离，上述偏移距离可以用网格数表示。第一获得子单元5021对比每个炮点与对应的各个所述检波点之间的偏移距离，从中获得最大的偏移距离，作为每个炮点相对于检波点的最大偏移距离，对于每个炮点都会有一个所述最大偏移距离。可理解的是，如果所述原始旅行时场数据为二维数据，那么炮点与对应的检波点的偏移距离是指炮点和对应的检波点在水平方向上的距离；如果所述原始旅行时场数据为三维数据，那么炮点与对应的检波点的偏移距离是指炮点和对应的检波点之间的直线距离在水平面上的投影。
在获得每个炮点相对于检波点的最大偏移距离之后，第二获得子单元5022根据每个炮点相对于检波点的最大偏移距离对每个炮点的原始旅行时场数据进行限定，即保留每个炮点的原始旅行时场数据中与每个炮点的偏移距离在所述最大偏移距离范围内的数据，舍弃每个炮点的原始旅行时场数据中与每个炮点的偏移距离大于所述最大偏移距离的数据，获得每个炮点的第一限定旅行时场数据。
由于已知每个炮点以及对应的各个所述检波点，第三获得子单元5023可以根据每个检波点，在每个炮点以及对应的各个所述检波点中查找到每个检波点对应的各个炮点。
第四获得子单元5024用于根据每个检波点以及对应的各个所述炮点的网格坐标，获得每个检波点相对于炮点的最大偏移距离的具体实现过程与第一获得子单元5021用于根据每个炮点以及对应的各个所述检波点的网格坐标，获得每个炮点相对于检波点的最大偏移距离的具体过程类似，此处不进行赘述。
第五获得子单元5025用于根据每个检波点相对于炮点的最大偏移距离对每个检波点的原始旅行时场数据进行限定，获得每个检波点的第一限定旅行时场数据具体实现过程，与第二获得子单元5022用于根据每个炮点相对于检波点的最大偏移距离对每个炮点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据的具体过程类似，此处不进行赘述。
图7是本发明又一实施例提供的基于菲涅尔体的旅行时场数据处理装置的结构示意图，如图7所示，在上述各实施例的基础上，进一步地，第三获得单元503包括第六获得子单元5031和第七获得子单元5032，其中：
第六获得子单元5031用于根据各个所述炮点以及各个所述检波点的在垂直方向上的网格坐标，获得高程限定边界；第七获得子单元5032用于根据所述高程限定边界对每个炮点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据，并根据所述高程限定条件对每个检波点的第一限定旅行时场数据进行限定，获得每个检波点的第二限定旅行时场数据。
具体地，所述旅行时初至数据中每个炮点的高程转换成了每个炮点在垂直方向上的网格坐标，所述旅行时初至数据中每个检波点的高程转换成了每个检波点在垂直方向上的网格坐标。第六获得子单元5031根据各个所述炮点以及各个所述检波点的在垂直方向上的网格坐标，可以拟合出高程限定边界。其中，如果所述旅行时初至数据为二维数据，那么拟合出高程限定边界为一条曲线。如果所述旅行时初至数据为三维数据，那么拟合出高程限定边界为一个曲面。
具体地，在获得所述高程限定边界之后，第七获得子单元5032可以根据所述高程限定边界对每个炮点的第一限定旅行时场数据进行限定，即将每个炮点的第一限定旅行时场数据中所述高程限定边界以上的数据舍去，获得每个炮点的第二限定旅行时场数据。所述服务器还可以根据所述高程限定边界对每个检波点的第一限定旅行时场数据进行限定，即将每个检波点的第一限定旅行时场数据中所述高程限定边界以上的数据舍去，获得每个检波点的第二限定旅行时场数据。
图8是本发明再一实施例提供的基于菲涅尔体的旅行时场数据处理装置的结构示意图，如图8所示，在上述各实施例的基础上，进一步地，第四获得单元504包括第八获得子单元5041、第九获得子单元5042、第十获得子单元5043和第十一获得子单元5044，其中：
第八获得子单元5041用于根据每个炮点的第二限定旅行时场数据以及每个炮点的最大偏移距离对应的检波点的第二限定旅行时场数据，获得每个炮点的最大偏移距离炮检对的旅行时场数据；第九获得子单元5042用于根据每个炮点的最大偏移距离炮检对的旅行时场数据，获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息；第十获得子单元5043用于根据每个炮点的最大偏移距离炮检对的菲涅尔体路径信息，获得每个炮点的最大菲涅尔体深度；第十一获得子单元5044用于根据每个炮点的最大菲涅尔体深度对每个炮点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据，并根据每个炮点的最大菲涅尔体深度对每个炮点对应的各个所述检波点的第二限定旅行时场数据进行限定，获得每个炮点对应的各个所述检波点的最终限定旅行时场数据。
具体地，每个炮点与对应的每个检波点之间存在偏移距离，第八获得子单元5041可以从每个炮点与对应的各个所述检波点之间的偏移距离中获得每个炮点相对于检波点的最大偏移距离，即每个炮点的最大偏移距离。第八获得子单元5041将每个炮点的第二限定旅行时场数据以及每个炮点的最大偏移距离对应的检波点的第二限定旅行时场数据进行叠加，可以获得每个炮点的最大偏移距离炮检对的旅行时场数据。其中，每个炮点与每个炮点最大偏移距离对应的检波点构成每个炮点的最大偏移距离炮检对。
在获得每个炮点的最大偏移距离炮检对的旅行时场数据之后，第九获得子单元5042根据每个炮点的最大偏移距离炮检对的旅行时场数据以及经验公式，可以计算获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息。
在获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息之后，第十获得子单元5043可以根据每个炮点的最大偏移距离炮检对的菲涅尔体路径信息，计算获得每个炮点的最大菲涅尔体深度。
在获得每个炮点的最大菲涅尔体深度之后，第十一获得子单元5044可以根据每个炮点的最大菲涅尔体深度对每个炮点的第二限定旅行时场数据进行限定，即保留每个炮点的第二限定旅行时场数据中在每个炮点的最大菲涅尔体深度范围内的数据，舍弃每个炮点的第二限定旅行时场数据中超出每个炮点的最大菲涅尔体深度范围的数据，获得每个所述炮点的最终限定旅行时场数据。第十一获得子单元5044还可以根据每个炮点的最大菲涅尔体深度对每个炮点对应的各个检波点的第二限定旅行时场数据进行限定，即保留每个炮点对应的每个检波点的第二限定旅行时场数据中在每个炮点的最大菲涅尔体深度范围内的数据，舍弃每个炮点对应的每个检波点的第二限定旅行时场数据中超出每个炮点的最大菲涅尔体深度范围的数据，获得每个炮点对应的各个所述检波点的最终限定旅行时场数据。
本发明实施例提供的装置的实施例具体可以用于执行上述各方法实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
图9是本发明一实施例提供的电子设备的实体结构示意图，如图9所示，该电子设备可以包括：处理器(processor)901、通信接口(Communications Interface)902、存储器(memory)903和通信总线904，其中，处理器901，通信接口902，存储器903通过通信总线904完成相互间的通信。处理器901可以调用存储器903中的逻辑指令，以执行如下方法：根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
此外，上述的存储器903中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本实施例公开一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例所提供的方法，例如包括：根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
本实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机程序，所述计算机程序使所述计算机执行上述各方法实施例所提供的方法，例如包括：根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在本说明书的描述中，参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于菲涅尔体的旅行时场数据处理方法，其特征在于，包括：
根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；
根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；
根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；
根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
2.根据权利要求1所述的方法，其特征在于，所述根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据包括：
根据每个炮点以及对应的各个所述检波点的网格坐标，获得每个炮点相对于检波点的最大偏移距离；
根据每个炮点相对于检波点的最大偏移距离对每个炮点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据；和
根据每个炮点以及对应的各个所述检波点，获得每个检波点对应的各个所述炮点；
根据每个检波点以及对应的各个所述炮点的网格坐标，获得每个检波点相对于炮点的最大偏移距离；
根据每个检波点相对于炮点的最大偏移距离对每个检波点的原始旅行时场数据进行限定，获得每个检波点的第一限定旅行时场数据。
3.根据权利要求1所述的方法，其特征在于，所述根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据包括：
根据各个所述炮点以及各个所述检波点的在垂直方向上的网格坐标，获得高程限定边界；
根据所述高程限定边界对每个炮点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据，并根据所述高程限定条件对每个检波点的第一限定旅行时场数据进行限定，获得每个检波点的第二限定旅行时场数据。
4.根据权利要求1至3任一项所述的方法，其特征在于，所述根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据包括：
根据每个炮点的第二限定旅行时场数据以及每个炮点的最大偏移距离对应的检波点的第二限定旅行时场数据，获得每个炮点的最大偏移距离炮检对的旅行时场数据；
根据每个炮点的最大偏移距离炮检对的旅行时场数据，获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息；
根据每个炮点的最大偏移距离炮检对的菲涅尔体路径信息，获得每个炮点的最大菲涅尔体深度；
根据每个炮点的最大菲涅尔体深度对每个炮点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据，并根据每个炮点的最大菲涅尔体深度对每个炮点对应的各个所述检波点的第二限定旅行时场数据进行限定，获得每个炮点对应的各个所述检波点的最终限定旅行时场数据。
5.一种基于菲涅尔体的旅行时场数据处理装置，其特征在于，包括：
第一获得单元，用于根据旅行时初至数据以及初始速度模型，获得每个炮点和每个检波点的原始旅行时场数据；
第二获得单元，用于根据第一数据限定规则对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据；
第三获得单元，用于根据第二数据限定规则对每个炮点的第一限定旅行时场数据和每个检波点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据；
第四获得单元，用于根据第三数据限定规则对每个炮点的第二限定旅行时场数据和每个检波点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。
6.根据权利要求5所述的装置，其特征在于，所述第二获得单元包括：
第一获得子单元，用于根据每个炮点以及对应的各个所述检波点的网格坐标，获得每个炮点相对于检波点的最大偏移距离；
第二获得子单元，用于根据每个炮点相对于检波点的最大偏移距离对每个炮点的原始旅行时场数据进行限定，获得每个炮点的第一限定旅行时场数据；
第三获得子单元，用于根据每个炮点以及对应的各个所述检波点，获得每个检波点对应的各个所述炮点；
第四获得子单元，用于根据每个检波点以及对应的各个所述炮点的网格坐标，获得每个检波点相对于炮点的最大偏移距离；
第五获得子单元，用于根据每个检波点相对于炮点的最大偏移距离对每个检波点的原始旅行时场数据进行限定，获得每个检波点的第一限定旅行时场数据。
7.根据权利要求5所述的装置，其特征在于，所述第三获得单元包括：
第六获得子单元，用于根据各个所述炮点以及各个所述检波点的在垂直方向上的网格坐标，获得高程限定边界；
第七获得子单元，用于根据所述高程限定边界对每个炮点的第一限定旅行时场数据进行限定，获得每个炮点的第二限定旅行时场数据，并根据所述高程限定条件对每个检波点的第一限定旅行时场数据进行限定，获得每个检波点的第二限定旅行时场数据。
8.根据权利要求5至7任一项所述的装置，其特征在于，所述第四获得单元包括：
第八获得子单元，用于根据每个炮点的第二限定旅行时场数据以及每个炮点的最大偏移距离对应的检波点的第二限定旅行时场数据，获得每个炮点的最大偏移距离炮检对的旅行时场数据；
第九获得子单元，用于根据每个炮点的最大偏移距离炮检对的旅行时场数据，获得每个炮点的最大偏移距离炮检对的菲涅尔体路径信息；
第十获得子单元，用于根据每个炮点的最大偏移距离炮检对的菲涅尔体路径信息，获得每个炮点的最大菲涅尔体深度；
第十一获得子单元，用于根据每个炮点的最大菲涅尔体深度对每个炮点的第二限定旅行时场数据进行限定，获得每个炮点的最终限定旅行时场数据，并根据每个炮点的最大菲涅尔体深度对每个炮点对应的各个所述检波点的第二限定旅行时场数据进行限定，获得每个炮点对应的各个所述检波点的最终限定旅行时场数据。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至4任一项所述方法的步骤。
CN201911043097.3A 2019-10-30 2019-10-30 一种基于菲涅尔体的旅行时场数据处理方法及装置 Active CN112748466B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911043097.3A CN112748466B (zh) 2019-10-30 2019-10-30 一种基于菲涅尔体的旅行时场数据处理方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911043097.3A CN112748466B (zh) 2019-10-30 2019-10-30 一种基于菲涅尔体的旅行时场数据处理方法及装置

Publications (2)

Publication Number Publication Date
CN112748466A true CN112748466A (zh) 2021-05-04
CN112748466B CN112748466B (zh) 2024-03-26

Family

ID=75640366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911043097.3A Active CN112748466B (zh) 2019-10-30 2019-10-30 一种基于菲涅尔体的旅行时场数据处理方法及装置

Country Status (1)

Country Link
CN (1) CN112748466B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102692644A (zh) * 2012-05-30 2012-09-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 生成深度域成像道集的方法
CN102937721A (zh) * 2012-11-07 2013-02-20 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用初至波走时的有限频层析成像方法
CN103675894A (zh) * 2013-12-24 2014-03-26 中国海洋石油总公司 一种基于三维高斯束射线追踪和频率域合成地震记录方法
CN104849751A (zh) * 2015-05-15 2015-08-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 叠前地震资料成像的方法
CN106814391A (zh) * 2015-11-27 2017-06-09 中国石油化工股份有限公司 基于菲涅尔体层析反演的地面微地震事件定位方法
US20170276818A1 (en) * 2016-03-24 2017-09-28 Saudi Arabian Oil Company Simultaneous wavefield reconstruction and receiver deghosting of seismic streamer data using an l1 inversion
CN107329169A (zh) * 2017-07-28 2017-11-07 中国石油天然气股份有限公司 一种角度道集的提取方法及装置
CN109144405A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 一种旅行时数据缓存方法与装置
CN110095773A (zh) * 2019-06-03 2019-08-06 中南大学 探地雷达多尺度全波形双参数反演方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102692644A (zh) * 2012-05-30 2012-09-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 生成深度域成像道集的方法
CN102937721A (zh) * 2012-11-07 2013-02-20 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用初至波走时的有限频层析成像方法
CN103675894A (zh) * 2013-12-24 2014-03-26 中国海洋石油总公司 一种基于三维高斯束射线追踪和频率域合成地震记录方法
CN104849751A (zh) * 2015-05-15 2015-08-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 叠前地震资料成像的方法
CN106814391A (zh) * 2015-11-27 2017-06-09 中国石油化工股份有限公司 基于菲涅尔体层析反演的地面微地震事件定位方法
US20170276818A1 (en) * 2016-03-24 2017-09-28 Saudi Arabian Oil Company Simultaneous wavefield reconstruction and receiver deghosting of seismic streamer data using an l1 inversion
CN109144405A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 一种旅行时数据缓存方法与装置
CN107329169A (zh) * 2017-07-28 2017-11-07 中国石油天然气股份有限公司 一种角度道集的提取方法及装置
CN110095773A (zh) * 2019-06-03 2019-08-06 中南大学 探地雷达多尺度全波形双参数反演方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
ERIKA BARISON ET AL: "《Wave equation datuming applied to marine OBS data and to land high resolution seismic profiling》", JOURNAL OF APPLIED GEOPHYSICS, pages 267 - 277 *
刘小民等: "《潜水波胖射线走时层析速度反演及其在深度偏移速度建模中的应用》", 石油物探, vol. 56, no. 5, pages 718 - 726 *
张东等: "《三维旅行时场B样条插值射线追踪方法》", 石油地球物理勘探, vol. 48, no. 4, pages 559 - 566 *
张亨等: "《快速菲涅尔体初至波层析的实现与应用》", 中国石油学会2021年物探技术研讨会论文集, pages 125 - 128 *
张建中等: "《Fresnel层析成像并行算法研究》", 计算机研究与发展, vol. 44, no. 10, pages 1661 - 1666 *
杨旭明: "《菲涅耳体初至层析静校正技术的研究与应用》", 复杂油气藏, vol. 4, no. 1, pages 28 - 32 *
王宏伟: "《菲涅耳体地震层析》", 石油地球物理勘探, vol. 31, no. 3, pages 309 - 326 *
王祥春等: "《起伏地表二维声波方程地震波场延拓》", 中国矿业大学学报, vol. 38, no. 5, 30 September 2009 (2009-09-30), pages 702 - 706 *

Also Published As

Publication number Publication date
CN112748466B (zh) 2024-03-26

Similar Documents

Publication Publication Date Title
US11281229B2 (en) Method and apparatus for outputting obstacle information
CN105277978B (zh) 一种确定近地表速度模型的方法及装置
US20150109290A1 (en) Device and method for removing noise points in point clouds
CN106896403B (zh) 弹性高斯束偏移成像方法和系统
CA2838857C (en) Seismic data processing including true-azimuth three-dimensional internal multiple attentuation without subsurface information
CN104850682A (zh) 基于位置的多点地质统计学建模方法
CN111596347A (zh) 一种快速求取表层纵横波速度比的方法及装置
CN104181593A (zh) 一种三维无射线追踪回折波层析成像方法及装置
CN112099082B (zh) 一种共面元共方位角道集的地震回折波走时反演方法
CA2953100A1 (en) System and method for rock property estimation of subsurface geologic volumes
CN106908836B (zh) 采集脚印压制方法及系统
CN112748466A (zh) 一种基于菲涅尔体的旅行时场数据处理方法及装置
CN112904424B (zh) 一种异常体位置的确定方法、装置及可读存储介质
CN109188527B (zh) 滩浅海地区快速建立三维近海底速度模型的方法
CN113866827B (zh) 一种解释性速度建模地震成像方法、系统、介质和设备
CN104570091B (zh) 一种获得初至波射线的方法和装置
CN110390702B (zh) 基于地质走向分析的变面元平面插值方法
CN110376642B (zh) 一种基于锥面波的三维地震速度反演方法
CN111538077A (zh) 基于倾角约束的叠前深度偏移方法及装置
CN112083486A (zh) 一种低速层速度获取方法及装置
CN112379420B (zh) 高精度弯曲测线叠前时间域成像方法及装置
CN112305600B (zh) 一种污染道横波初至获取方法及装置
CN112394415B (zh) 海底节点水平方位角的确定方法及装置
CN104133238B (zh) 一种基于预计算插值的射线走时计算方法及系统
CN115761533B (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