阅读说明：本技术 一种基于菲涅尔体的旅行时场数据处理方法及装置 (Travel time field data processing method and device based on Fresnel body ) 是由 张亨 彭文 张华� 金德刚 刘鸿 李乐 于 2019-10-30 设计创作，主要内容包括：本发明提供一种基于菲涅尔体的旅行时场数据处理方法及装置,所述方法包括：根据旅行时初至数据以及初始速度模型,获得每个炮点和每个检波点的原始旅行时场数据；对每个炮点的原始旅行时场数据和每个检波点的原始旅行时场数据依次采用第一数据限定规则、第二数据限定规则和第三数据限定规则进行数据限定,最终获得每个炮点的最终限定旅行时场数据和每个检波点的最终限定旅行时场数据。所述装置用于执行上述方法。本发明实施例提供的基于菲涅尔体的旅行时场数据处理方法及装置,提高了旅行时场数据的处理效率。(The invention provides a traveling time field data processing method and a traveling time field data processing device based on a Fresnel body, wherein the method comprises the following steps: acquiring original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model; and sequentially adopting a first data definition rule, a second data definition rule and a third data definition rule to carry out data definition on the original travel time field data of each shot point and the original travel time field data of each demodulator probe, and finally obtaining the final definition travel time field data of each shot point and the final definition travel time field data of each demodulator probe. The device is used for executing the method. The method and the device for processing the travel time field data based on the Fresnel body improve the processing efficiency of the travel time field data.)

1. A travel time field data processing method based on a Fresnel body is characterized by comprising the following steps:

acquiring original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model;

defining the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe according to a first data definition rule to obtain first defined travel-time field data of each shot point and first defined travel-time field data of each demodulator probe;

defining the first defined travel-time field data of each shot point and the first defined travel-time field data of each demodulator probe according to a second data definition rule to obtain second defined travel-time field data of each shot point and second defined travel-time field data of each demodulator probe;

and defining the second defined travel-time field data of each shot point and the second defined travel-time field data of each demodulator probe according to a third data definition rule, and obtaining the final defined travel-time field data of each shot point and the final defined travel-time field data of each demodulator probe.

2. The method of claim 1, wherein the defining the original travel-time field data for each shot and the original travel-time field data for each demodulator probe according to a first data definition rule, and wherein obtaining the first defined travel-time field data for each shot and the first defined travel-time field data for each demodulator probe comprises:

obtaining the maximum offset distance of each shot point relative to the demodulator probe according to the grid coordinates of each shot point and each corresponding demodulator probe;

limiting the original travel time field data of each shot point according to the maximum offset distance of each shot point relative to the demodulator probe to obtain first limited travel time field data of each shot point; and

obtaining each shot point corresponding to each wave detection point according to each shot point and each corresponding wave detection point;

obtaining the maximum offset distance of each demodulator probe relative to the shot point according to each demodulator probe and the corresponding grid coordinate of each shot point;

and defining the original travel-time field data of each demodulator probe according to the maximum offset distance of each demodulator probe relative to the shot point to obtain first defined travel-time field data of each demodulator probe.

3. The method of claim 1, wherein the defining the first defined traveltime field data for each shot and the first defined traveltime field data for each demodulator probe according to a second data definition rule, the obtaining the second defined traveltime field data for each shot and the second defined traveltime field data for each demodulator probe comprising:

obtaining an elevation limit boundary according to the grid coordinates of each shot point and each demodulator probe in the vertical direction;

and defining the first limited travel-time field data of each shot point according to the elevation defining boundary to obtain second limited travel-time field data of each shot point, and defining the first limited travel-time field data of each wave detection point according to the elevation defining condition to obtain the second limited travel-time field data of each wave detection point.

4. A method according to any one of claims 1 to 3, wherein the defining of the second defined traveltime field data for each shot and the second defined traveltime field data for each demodulator probe according to a third data definition rule, the obtaining of the final defined traveltime field data for each shot and the final defined traveltime field data for each demodulator probe comprising:

acquiring travel time field data of the shot-survey pair with the maximum offset distance of each shot point according to the second limited travel time field data of each shot point and the second limited travel time field data of the corresponding demodulator probe with the maximum offset distance of each shot point;

acquiring Fresnel volume path information of the maximum offset distance shot-inspection pair of each shot point according to the travel time field data of the maximum offset distance shot-inspection pair of each shot point;

acquiring the maximum Fresnel body depth of each shot point according to the Fresnel body path information of the maximum offset distance shot-check pair of each shot point;

and limiting the second limited travel time field data of each shot point according to the maximum Fresnel body depth of each shot point to obtain the final limited travel time field data of each shot point, and limiting the second limited travel time field data of each demodulator probe corresponding to each shot point according to the maximum Fresnel body depth of each shot point to obtain the final limited travel time field data of each demodulator probe corresponding to each shot point.

5. A travel time field data processing device based on a Fresnel body is characterized by comprising:

the first obtaining unit is used for obtaining original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial velocity model;

a second obtaining unit, configured to define, according to a first data definition rule, the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe, and obtain first defined travel-time field data of each shot point and first defined travel-time field data of each demodulator probe;

a third obtaining unit, configured to define the first defined travel-time field data of each shot point and the first defined travel-time field data of each geophone point according to a second data definition rule, and obtain the second defined travel-time field data of each shot point and the second defined travel-time field data of each geophone point;

and a fourth obtaining unit, configured to define the second defined travel-time field data of each shot point and the second defined travel-time field data of each geophone point according to a third data definition rule, and obtain the final defined travel-time field data of each shot point and the final defined travel-time field data of each geophone point.

6. The apparatus of claim 5, wherein the second obtaining unit comprises:

the first obtaining subunit is configured to obtain a maximum offset distance of each shot point relative to a corresponding demodulator probe according to each shot point and the corresponding grid coordinate of each demodulator probe;

the second obtaining subunit is used for limiting the original travel time field data of each shot point according to the maximum offset distance of each shot point relative to the demodulator probe, and obtaining the first limited travel time field data of each shot point;

a third obtaining subunit, configured to obtain, according to each shot point and each corresponding demodulator probe, each shot point corresponding to each demodulator probe;

a fourth obtaining subunit, configured to obtain, according to each demodulator probe and the grid coordinate of each corresponding shot point, a maximum offset distance of each demodulator probe with respect to a shot point;

and the fifth obtaining subunit is used for limiting the original travel time field data of each demodulator probe according to the maximum offset distance of each demodulator probe relative to the shot point, and obtaining the first limited travel time field data of each demodulator probe.

7. The apparatus of claim 5, wherein the third obtaining unit comprises:

a sixth obtaining subunit, configured to obtain an elevation-limited boundary according to the grid coordinates of each shot point and each geophone point in the vertical direction;

and the seventh obtaining subunit is configured to define the first defined travel-time field data of each shot point according to the elevation defining boundary, obtain second defined travel-time field data of each shot point, and define the first defined travel-time field data of each geophone point according to the elevation defining condition, so as to obtain the second defined travel-time field data of each geophone point.

8. The apparatus according to any one of claims 5 to 7, wherein the fourth obtaining unit comprises:

the eighth obtaining subunit is configured to obtain travel time field data of a shot-survey pair with the maximum offset distance of each shot point according to the second limited travel time field data of each shot point and the second limited travel time field data of a corresponding demodulator probe with the maximum offset distance of each shot point;

a ninth obtaining subunit, configured to obtain fresnel volume path information of the maximum offset distance shot-inspection pair of each shot point according to the travel time field data of the maximum offset distance shot-inspection pair of each shot point;

a tenth obtaining subunit, configured to obtain, according to the fresnel volume path information of the maximum offset distance shot-check pair of each shot point, a maximum fresnel volume depth of each shot point;

and the eleventh obtaining subunit is configured to define the second limited travel time field data of each shot point according to the maximum fresnel volume depth of each shot point, obtain the final limited travel time field data of each shot point, and define the second limited travel time field data of each demodulator probe corresponding to each shot point according to the maximum fresnel volume depth of each shot point, so as to obtain the final limited travel time field data of each demodulator probe corresponding to each shot point.

9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method of any of claims 1 to 4 are implemented when the computer program is executed by the processor.

10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 4.

Technical Field

The invention relates to the technical field of geological exploration, in particular to a traveling time field data processing method and device based on a Fresnel body.

Background

At present, the application of the Fresnel body chromatography method in near-surface velocity inversion and time-lapse seismic monitoring of carbon dioxide reservoir shows great potential, and is one of the very important technologies for near-surface modeling.

The Fresnel volume ray tracing technology considers the influence of the speed on the optimal high-frequency ray path on the seismic wave trend, also considers the influence on the position in the range of the first Fresnel zone in combination with frequency information, and better accords with the law of the actual propagation of the seismic waves. Moreover, the sparse chromatographic inversion coefficient matrix caused by the rare number of high-frequency rays and the uneven coverage can be compensated, so that the problem of unstable solving process is solved, a more accurate inversion result is obtained, and the inversion stability is improved. However, the Fresnel body ray tracing result information is rich, so that the computation complexity of the Fresnel body is high, the time consumption is huge, especially in the application process of a three-dimensional work area, the computation amount of the Fresnel body ray tracing is increased in geometric orders of magnitude due to the steep increase of the three-dimensional first arrival information, and the computation efficiency is very low.

Therefore, how to provide a method for processing data in a travel time field based on a fresnel body, which can reduce the data processing amount in the travel time field and improve the data processing efficiency in the travel time field becomes an important issue to be solved in the field.

Disclosure of Invention

Aiming at the problems in the prior art, the embodiment of the invention provides a method and a device for processing travel time field data based on a Fresnel body.

On one hand, the invention provides a traveling time field data processing method based on a Fresnel body, which comprises the following steps:

acquiring original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model;

defining the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe according to a first data definition rule to obtain first defined travel-time field data of each shot point and first defined travel-time field data of each demodulator probe;

defining the first defined travel-time field data of each shot point and the first defined travel-time field data of each demodulator probe according to a second data definition rule to obtain second defined travel-time field data of each shot point and second defined travel-time field data of each demodulator probe;

and defining the second defined travel-time field data of each shot point and the second defined travel-time field data of each demodulator probe according to a third data definition rule, and obtaining the final defined travel-time field data of each shot point and the final defined travel-time field data of each demodulator probe.

In another aspect, the present invention provides a data processing apparatus for a travel time field based on a fresnel body, including:

the first obtaining unit is used for obtaining original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial velocity model;

a second obtaining unit, configured to define, according to a first data definition rule, the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe, and obtain first defined travel-time field data of each shot point and first defined travel-time field data of each demodulator probe;

a third obtaining unit, configured to define the first defined travel-time field data of each shot point and the first defined travel-time field data of each geophone point according to a second data definition rule, and obtain the second defined travel-time field data of each shot point and the second defined travel-time field data of each geophone point;

and a fourth obtaining unit, configured to define the second defined travel-time field data of each shot point and the second defined travel-time field data of each geophone point according to a third data definition rule, and obtain the final defined travel-time field data of each shot point and the final defined travel-time field data of each geophone point.

In another aspect, the present invention provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the program, the steps of the fresnel volume-based travel time field data processing method according to any one of the embodiments described above are implemented.

In still another aspect, the present invention provides a computer-readable storage medium, on which a computer program is stored, wherein the computer program, when executed by a processor, implements the steps of the fresnel body-based travel time field data processing method according to any one of the above embodiments.

The method and the device for processing the travel time field data based on the Fresnel body, provided by the embodiment of the invention, can obtain the original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model, then limit the original travel time field data of each shot point and the original travel time field data of each demodulator probe according to the first data limiting rule, obtain the first limited travel time field data of each shot point and the first limited travel time field data of each demodulator probe, then limit the first limited travel time field data of each shot point and the first limited travel time field data of each demodulator probe according to the second data limiting rule, obtain the second limited travel time field data of each shot point and the second limited travel time field data of each demodulator probe, and then limit the second limited travel time field data of each shot point and the second limited travel time field data of each demodulator probe according to the third data limiting rule, the final travel time field limiting data of each shot point and the final travel time field limiting data of each demodulator probe are obtained, the range of the travel time field data is narrowed, the number of Fresnel body calculation is reduced, and the processing efficiency of the travel time field data is improved.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:

fig. 1 is a schematic flow chart of a travel time field data processing method based on a fresnel body according to an embodiment of the present invention.

Fig. 2a is a schematic flow chart of a fresnel volume-based travel time field data processing method according to another embodiment of the present invention.

Fig. 2b is a schematic flow chart of a fresnel-volume-based travel time field data processing method according to still another embodiment of the present invention.

Fig. 3 is a schematic flow chart of a fresnel volume-based travel time field data processing method according to yet another embodiment of the present invention.

Fig. 4 is a schematic flow chart of a fresnel volume-based travel time field data processing method according to still another embodiment of the present invention.

Fig. 5 is a schematic structural diagram of a fresnel volume-based travel time field data processing apparatus according to an embodiment of the present invention.

Fig. 6 is a schematic structural diagram of a fresnel volume-based travel time field data processing apparatus according to another embodiment of the present invention.

Fig. 7 is a schematic structural diagram of a fresnel volume-based travel time field data processing apparatus according to yet another embodiment of the present invention.

Fig. 8 is a schematic structural diagram of a fresnel volume-based travel time field data processing apparatus according to still another embodiment of the present invention.

Fig. 9 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention. It should be noted that the embodiments and features of the embodiments in the present application may be arbitrarily combined with each other without conflict.

Fig. 1 is a schematic flow chart of a travel time field data processing method based on a fresnel body according to an embodiment of the present invention, and as shown in fig. 1, the travel time field data processing method based on a fresnel body according to an embodiment of the present invention includes:

s101, obtaining original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model;

specifically, after obtaining the travel-time first arrival data and the initial velocity model of the target area, the server may obtain the original travel-time field data of each shot point and each geophone according to the travel-time first arrival data and the initial velocity model, for example, the original travel-time field data of each shot point and each geophone may be obtained according to the travel-time first arrival data and the initial velocity model by using Fast Marching Method (FMA) algorithm. The travel time first arrival data comprises the coordinates and the elevation of each shot point and the coordinates and the elevation of a corresponding demodulator probe of each shot point. The initial velocity model is formed by utilizing coordinates and elevation interpolation of all shot points and demodulator probes in travel-time first-arrival data of the target area, and has a certain grid number in the horizontal direction x and the vertical direction z for a two-dimensional initial velocity model and a certain grid number in the directions x and z and the direction y perpendicular to the xz plane for a three-dimensional initial velocity model. It will be appreciated that the travel-time first-arrival data for the target area is obtained after processing the raw seismic data for the target area. The execution subject of the automatic super-parameter determining method provided by the embodiment of the invention includes, but is not limited to, a server.

For example, for the target area, in order to obtain a two-dimensional initial velocity model, a first preset number of shot points may be set on the ground surface of the target area along a straight line, each shot point corresponds to a second preset number of wave detection points, the second preset number of shot points are uniformly distributed on both sides of the corresponding shot point, a geophone arranged at each wave detection point is used to collect original seismic data, and two-dimensional travel-time first arrival data of the target area may be obtained through processing by relevant software, such as Geoeast software, tomodel software, omega and other software in the industry, and then the two-dimensional travel-time first arrival data may be processed to obtain a two-dimensional initial velocity model. And the server can obtain the original travel time field data of each shot point and each demodulator probe of the target area according to the two-dimensional travel time first arrival data and the two-dimensional initial velocity model. The first preset quantity and the second preset quantity are set according to actual needs, and the embodiment of the invention is not limited.

For example, for the target area, in order to obtain a three-dimensional initial velocity model, a third preset number of shots may be set on the surface of the target area, each shot having a plurality of receiving lines and each shot being located at an intersection of the plurality of receiving lines, and a fourth preset number of detector points being set on each receiving line. The method comprises the steps of collecting original seismic data by using a detector arranged at each wave detection point, processing the original seismic data by related software to obtain three-dimensional travel-time first-arrival data of a target area, and processing the three-dimensional travel-time first-arrival data to obtain a three-dimensional initial velocity model. And the server can obtain the original travel time field data of each shot point and each demodulator probe of the target area according to the three-dimensional travel time first arrival data and the three-dimensional initial velocity model. The third preset number and the fourth preset number are set according to actual needs, and the embodiment of the invention is not limited.

S102, limiting the original travel time field data of each shot point and the original travel time field data of each demodulator probe according to a first data limiting rule to obtain first limited travel time field data of each shot point and first limited travel time field data of each demodulator probe;

specifically, after obtaining the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe, the server limits the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe according to a first data limiting rule, and can obtain the first limited travel-time field data of each shot point and the first limited travel-time field data of each demodulator probe. Wherein the first data definition rule may be a definition of original travel-time field data for each shot and each geophone point in an offset direction, an elevation direction, or a depth direction. If the original travel time field data is two-dimensional data, the offset direction refers to a horizontal direction, the height direction refers to an upward direction perpendicular to the horizontal direction, and the depth direction refers to a downward direction perpendicular to the horizontal direction. If the original travel time field data is three-dimensional data, the offset direction refers to an up-horizontal direction, the height direction refers to a direction perpendicular to a horizontal plane and the depth direction refers to a direction perpendicular to a horizontal plane and downward.

For example, the server may obtain a maximum offset distance of each shot from the geophone point according to the grid coordinates of each shot and the corresponding respective geophone point, and then define the original travel-time field data of each shot according to the maximum offset distance of each shot from the geophone point, thereby obtaining first defined travel-time field data of each shot. The server may obtain each shot point corresponding to each pickup point according to each shot point and each corresponding pickup point, then obtain the maximum offset distance of each pickup point relative to the shot point according to the grid coordinates of each pickup point and each corresponding shot point, and then limit the original travel time field data of each pickup point according to the maximum offset distance of each pickup point relative to the shot point to obtain the first limited travel time field data of each pickup point.

S103, limiting the first limited travel time field data of each shot point and the first limited travel time field data of each demodulator probe according to a second data limiting rule to obtain second limited travel time field data of each shot point and second limited travel time field data of each demodulator probe;

specifically, after obtaining the first defined travel-time field data of each shot and the first defined travel-time field data of each geophone, the server may define the first defined travel-time field data of each shot and the first defined travel-time field data of each geophone according to a second data definition rule, and obtain the second defined travel-time field data of each shot and the second defined travel-time field data of each geophone. Wherein the second data definition rule is a definition of first defined travel-time-field data for each shot and each geophone point in the offset direction, the elevation direction, or the depth direction. It is to be understood that the second data definition rule is different from the first data definition rule, i.e. if the first data definition rule is a data definition in the offset direction, the second data definition rule is a data definition in the height direction or in the depth direction, if the first data definition rule is a data definition in the height direction, the second data definition rule is a data definition in the offset direction or in the depth direction; if the first data definition rule is a data definition in the depth direction, the second data definition rule is a data definition in the height direction or in the depth direction.

For example, the server obtains an elevation limit boundary according to grid coordinates of each shot point and each geophone point in the vertical direction, then limits first limit travel time field data of each shot point according to the elevation limit boundary, obtains second limit travel time field data of each shot point, and limits the first limit travel time field data of each geophone point according to the elevation limit condition, so as to obtain the second limit travel time field data of each geophone point.

And S104, limiting the second limited travel time field data of each shot point and the second limited travel time field data of each demodulator probe according to a third data limiting rule to obtain the final limited travel time field data of each shot point and the final limited travel time field data of each demodulator probe.

Specifically, after obtaining the second defined travel-time field data of each shot and the second defined travel-time field data of each geophone, the server may define the second defined travel-time field data of each shot and the second defined travel-time field data of each geophone according to a third data definition rule, and obtain the final defined travel-time field data of each shot and the final defined travel-time field data of each geophone. Wherein the third data definition rule is a definition of first defined travel-time-field data for each shot and each geophone point in the offset direction, the elevation direction, or the depth direction. It is to be understood that the third data definition rule is different from both the first data definition rule and the second data definition rule, i.e. if the third data definition rule is a data definition in the depth direction, the first data definition rule and the second data definition rule are data definitions in the offset direction and in the height direction; if the third data definition rule is a data definition in an offset direction, the first data definition rule and the second data definition rule are data definitions in a depth direction and a height direction; if the third data defining rule is a data definition in the height direction, the first data defining rule and the second data defining rule are definitions in the offset direction and in the depth direction.

For example, the server may obtain the travel time field data of the maximum offset distance shot-check pair of each shot point according to the second limited travel time field data of each shot point and the second limited travel time field data of the corresponding demodulator probe of the maximum offset distance of each shot point, then obtain the fresnel body path information of the maximum offset distance shot-check pair of each shot point according to the travel time field data of the maximum offset distance shot-check pair of each shot point, then obtain the maximum fresnel body depth of each shot point according to the fresnel body path information of the maximum offset distance shot-check pair of each shot point, then limit the second limited travel time field data of each shot point according to the maximum fresnel body depth of each shot point, obtain the final limited travel time field data of each shot point, and limit the second limited travel time field data of each corresponding demodulator probe of each shot point according to the maximum fresnel body depth of each shot point, and acquiring final limited travel time field data of the demodulator probes corresponding to each shot point.

The method for processing the travel time field data based on the Fresnel body, provided by the embodiment of the invention, can obtain the original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model, then limit the original travel time field data of each shot point and the original travel time field data of each demodulator probe according to the first data limiting rule, obtain the first limited travel time field data of each shot point and the first limited travel time field data of each demodulator probe, then limit the first limited travel time field data of each shot point and the first limited travel time field data of each demodulator probe according to the second data limiting rule, obtain the second limited travel time field data of each shot point and the second limited travel time field data of each demodulator probe, and then limit the second limited travel time field data of each shot point and the second limited travel time field data of each demodulator probe according to the third data limiting rule, the final travel time field limiting data of each shot point and the final travel time field limiting data of each demodulator probe are obtained, the range of the travel time field data is narrowed, the number of Fresnel body calculation is reduced, and the processing efficiency of the travel time field data is improved.

Fig. 2a is a schematic flow chart of a fresnel-body-based travel time field data processing method according to another embodiment of the present invention, and fig. 2b is a schematic flow chart of a fresnel-body-based travel time field data processing method according to yet another embodiment of the present invention, where as shown in fig. 2a and fig. 2b, the defining the original travel time field data of each shot point and the original travel time field data of each demodulator probe according to a first data defining rule, and obtaining the first defined travel time field data of each shot point and the first defined travel time field data of each demodulator probe includes:

s1021, obtaining the maximum offset distance of each shot point relative to the demodulator probe according to each shot point and the grid coordinates of the corresponding demodulator probes;

specifically, for the raw travel-time-field data for each shot, each shot and the corresponding respective demodulator probe are represented in grid coordinates. The server can obtain the offset distance between each shot point and each corresponding demodulator probe according to the grid coordinates of each shot point and each corresponding demodulator probe, and the offset distance can be represented by the grid number. And the server compares the offset distance between each shot point and each corresponding demodulator probe to obtain the maximum offset distance from the offset distances, and the maximum offset distance is used as the maximum offset distance of each shot point relative to the demodulator probes, and one maximum offset distance is provided for each shot point. Understandably, if the original travel-time field data is two-dimensional data, the offset distance between a shot point and a corresponding demodulator probe refers to the distance between the shot point and the corresponding demodulator probe in the horizontal direction; if the original travel time field data is three-dimensional data, the offset distance between the shot point and the corresponding demodulator probe refers to the projection of the straight-line distance between the shot point and the corresponding demodulator probe on the horizontal plane.

S1022, limiting the original travel time field data of each shot point according to the maximum offset distance of each shot point relative to the demodulator probe to obtain first limited travel time field data of each shot point;

specifically, after obtaining the maximum offset distance of each shot point relative to a demodulator probe, the server limits the original travel time field data of each shot point according to the maximum offset distance of each shot point relative to the demodulator probe, namely, data of the original travel time field data of each shot point, of which the offset distance from each shot point is within the maximum offset distance range, is reserved, data of the original travel time field data of each shot point, of which the offset distance from each shot point is larger than the maximum offset distance, is discarded, and first limited travel time field data of each shot point is obtained.

For example, if the original travel-time field data for the shot is two-dimensional data, the grid coordinates of the shot are (x)₁，z₁)，x₁As a horizontal coordinate, z₁The maximum offset distance of the shot point relative to the demodulator probe is d as a vertical coordinate₁Then the horizontal coordinate in the original travel time field data of the shot is retained to be greater than or equal to x₁-d₁And is less than or equal to x₁+d₁As first defined travel-time-horizon data for the shot.

For example, if the original travel-time field data for the shot is three-dimensional data, the grid coordinates of the shot are (x)₂，y₂，z₂)，x₂Is a coordinate in a first direction on a horizontal plane, y₂Is a coordinate in a first direction on a horizontal plane, the first direction being perpendicular to a second direction, z₂The maximum offset distance of the shot point relative to the demodulator probe is d₂Then the coordinates of the first direction in the original travel time field data for that shot are retained to be greater than or equal to x₂-d₂And is less than or equal to x₂+d₂The coordinate in the second direction is greater than or equal to y₂-d₂And is less than or equal to y₂+d₂As first defined travel-time-horizon data for the shot.

S1023, obtaining each shot point corresponding to each wave detection point according to each shot point and each corresponding wave detection point;

specifically, since each shot point and each corresponding demodulator probe are known, the server may find each shot point corresponding to each demodulator probe in each shot point and each corresponding demodulator probe according to each demodulator probe.

For example, each demodulator probe has a corresponding number, the server searches the number of each demodulator probe corresponding to the shot point b according to the number of the demodulator probe a, and if the number of each demodulator probe corresponding to the shot point is the same as the number of the demodulator probe a, the shot point b is the shot point corresponding to the demodulator probe a. And repeating the process for the serial numbers of the wave detection points corresponding to each shot point to obtain each shot point corresponding to the wave detection point a.

S1024, obtaining the maximum offset distance of each demodulator probe relative to the shot point according to each demodulator probe and the corresponding grid coordinate of each shot point;

specifically, the specific implementation process of this step is similar to step S1021, and is not described herein again.

S1025, limiting the original travel time field data of each demodulator probe according to the maximum offset distance of each demodulator probe relative to the shot point, and obtaining first limited travel time field data of each demodulator probe.

Specifically, the specific implementation process of this step is similar to step S1022, and is not described herein again.

Fig. 3 is a schematic flow chart of a fresnel-volume-based travel time field data processing method according to yet another embodiment of the present invention, and as shown in fig. 3, based on the foregoing embodiments, further, the defining the first defined travel time field data of each shot point and the first defined travel time field data of each demodulator probe according to a second data definition rule, and obtaining the second defined travel time field data of each shot point and the second defined travel time field data of each demodulator probe includes:

s1031, obtaining elevation limit boundaries according to grid coordinates of the shot points and the demodulator probes in the vertical direction;

specifically, the elevation of each shot in the travel-time first arrival data is converted into grid coordinates of each shot in the vertical direction, and the elevation of each geophone in the travel-time first arrival data is converted into grid coordinates of each geophone in the vertical direction. And the server can fit an elevation limit boundary according to the grid coordinates of each shot point and each demodulator probe in the vertical direction. And if the travel time first arrival data are two-dimensional data, fitting an elevation limit boundary into a curve. And if the travel first arrival data are three-dimensional data, fitting the elevation limit boundary into a curved surface.

S1032, limiting the first limited travel time field data of each shot point according to the elevation limiting boundary to obtain second limited travel time field data of each shot point, and limiting the first limited travel time field data of each demodulator probe according to the elevation limiting condition to obtain the second limited travel time field data of each demodulator probe.

Specifically, after obtaining the elevation-defined boundary, the server may define the first defined travel-time field data of each shot according to the elevation-defined boundary, that is, truncating data above the elevation-defined boundary in the first defined travel-time field data of each shot, and obtaining the second defined travel-time field data of each shot. The server can also limit the first limited travel-time field data of each wave detection point according to the elevation limit boundary, namely, the data above the elevation limit boundary in the first limited travel-time field data of each wave detection point is omitted, and the second limited travel-time field data of each wave detection point is obtained.

Fig. 4 is a schematic flow chart of a fresnel-volume-based travel time field data processing method according to still another embodiment of the present invention, and as shown in fig. 4, based on the foregoing embodiments, further, the defining the second limited travel time field data of each shot point and the second limited travel time field data of each demodulator probe according to a third data definition rule, and obtaining the final limited travel time field data of each shot point and the final limited travel time field data of each demodulator probe includes:

s1041, obtaining travel time field data of the shot-survey pair with the maximum offset distance of each shot point according to the second limited travel time field data of each shot point and the second limited travel time field data of the corresponding wave detection point with the maximum offset distance of each shot point;

specifically, there is an offset distance between each shot point and each corresponding receiver point, and the server may obtain a maximum offset distance of each shot point relative to a receiver point, that is, a maximum offset distance of each shot point, from the offset distances between each shot point and the corresponding respective receiver points. The server superposes the second limited travel-time field data of each shot point and the second limited travel-time field data of the corresponding demodulator probe of the maximum offset distance of each shot point, so that the travel-time field data of the shot-receiver pair of the maximum offset distance of each shot point can be obtained. And each shot point and the corresponding demodulator probe of the maximum offset distance of each shot point form a shot detection pair of the maximum offset distance of each shot point.

For example, if the offset distance between the shot point c and the corresponding respective geophone points is the largest, the maximum offset distance of the shot point c is the offset distance between the shot point c and the geophone point d, and the geophone point d is the corresponding geophone point of the maximum offset distance of the shot point c, and the shot point c and the geophone point d form the maximum offset distance shot pair of the shot point c. The second defined traveltime field data for shot c is superimposed with the second defined traveltime field data for demodulator probe d to form the traveltime field data for the maximum offset distance geophone pair for shot c.

S1042, acquiring Fresnel body path information of the maximum offset distance shot-check pair of each shot point according to the travel time field data of the maximum offset distance shot-check pair of each shot point;

specifically, after obtaining the travel time field data of the maximum offset distance shot-check pair of each shot point, the server may calculate and obtain the fresnel volume path information of the maximum offset distance shot-check pair of each shot point according to the travel time field data of the maximum offset distance shot-check pair of each shot point and an empirical formula.

For example, the empirical formula is: t is t_L(g,r)+t_L(r,s)-t_Lmin(g, s) is less than or equal to 1/(2f), wherein g represents the ground wave detection point position, s represents the ground shot point position, t_L(g, r) is the time from any point r to the position of the detection point g in the field when the detection point travels, t_L(r, s) refers to the time from point r to point s in the point travel time field, t_Lmin(g, s) is the minimum worth position in the total time field of the sum of the shot and geophone travel time fields, and f represents the frequency, here preset, e.g. 25 Hz. The significance of the above empirical formula is that if the difference between the sum of the time from a location point to the shot point and the time to the demodulator probe and the minimum value of the total time field of the shot-survey journey is less than or equal to half of the given reciprocal frequency, the location point is considered to be within the range of the effective Fresnel volume.

S1043, obtaining the maximum Fresnel body depth of each shot point according to the Fresnel body path information of the maximum offset distance shot-check pair of each shot point;

specifically, after the server obtains the fresnel body path information of the maximum offset distance shot-test pair of each shot point, the server may calculate and obtain the maximum fresnel body depth of each shot point according to the fresnel body path information of the maximum offset distance shot-test pair of each shot point.

For example, after obtaining the fresnel volume path information of the maximum offset distance shot pair for each shot point, the fresnel volume range of the maximum offset distance shot pair for each shot point can be obtained, and the maximum value of the absolute value of the coordinates of each position point in the depth direction in the fresnel volume range is the maximum fresnel volume depth of each shot point.

S1044, limiting the second limited travel time field data of each shot point according to the maximum Fresnel body depth of each shot point to obtain the final limited travel time field data of each shot point, and limiting the second limited travel time field data of each demodulator probe corresponding to each shot point according to the maximum Fresnel body depth of each shot point to obtain the final limited travel time field data of each demodulator probe corresponding to each shot point.

Specifically, after obtaining the maximum fresnel depth of each shot point, the server may limit the second limited travel time field data of each shot point according to the maximum fresnel depth of each shot point, that is, data in the second limited travel time field data of each shot point within the maximum fresnel depth of each shot point is retained, data in the second limited travel time field data of each shot point exceeding the maximum fresnel depth of each shot point is discarded, and the final limited travel time field data of each shot point is obtained. The server can also limit the second limited travel time field data of each demodulator probe corresponding to each shot point according to the maximum fresnel body depth of each shot point, namely, data in the maximum fresnel body depth range of each shot point in the second limited travel time field data of each demodulator probe corresponding to each shot point is reserved, data exceeding the maximum fresnel body depth range of each shot point in the second limited travel time field data of each demodulator probe corresponding to each shot point is discarded, and finally limited travel time field data of each demodulator probe corresponding to each shot point is obtained.

For example, the original travel-time field data of the shot is two-dimensional data, and the grid coordinate of the shot is (x)₁，z₁)，x₁As a horizontal coordinate, z₁The maximum Fresnel body depth of the shot point is h as a vertical coordinate₁Then the vertical coordinate in the second defined travel time field data holding the shot is greater than or equal to z₁-h₁As final defined travel time field data for the shot.

For example, if the original travel-time field data for the shot is three-dimensional data, the grid coordinates of the shot are (x)₂，y₂，z₂)，x₂Is a coordinate in a first direction on a horizontal plane, y₂Is a coordinate in a first direction on a horizontal plane, the first direction being perpendicular to a second direction, z₂The maximum Fresnel body depth of the shot point is h which is a coordinate in the direction vertical to the horizontal plane₂Then the coordinate perpendicular to the horizontal in the second defined travel time field data for that shot is retained is greater than or equal to z₂-h₂As final defined travel time field data for the shot.

Fig. 5 is a schematic structural diagram of a fresnel-volume-based travel time field data processing apparatus according to an embodiment of the present invention, and as shown in fig. 5, the fresnel-volume-based travel time field data processing apparatus according to an embodiment of the present invention includes a first obtaining unit 501, a second obtaining unit 502, a third obtaining unit 503, and a fourth obtaining unit 504, where:

the first obtaining unit 501 is configured to obtain original travel-time field data of each shot point and each demodulator probe according to the travel-time first arrival data and the initial velocity model; the second obtaining unit 502 is configured to define the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe according to a first data definition rule, and obtain first defined travel-time field data of each shot point and first defined travel-time field data of each demodulator probe; the third obtaining unit 503 is configured to define the first defined travel-time field data of each shot point and the first defined travel-time field data of each geophone point according to a second data definition rule, and obtain the second defined travel-time field data of each shot point and the second defined travel-time field data of each geophone point; the fourth obtaining unit 504 is configured to define the second defined travel-time field data of each shot point and the second defined travel-time field data of each geophone point according to a third data definition rule, and obtain the final defined travel-time field data of each shot point and the final defined travel-time field data of each geophone point.

Specifically, after obtaining the travel-time first arrival data and the initial velocity model of the target area, the first obtaining unit 501 may obtain the original travel-time field data of each shot point and each geophone according to the travel-time first arrival data and the initial velocity model, and may obtain the original travel-time field data of each shot point and each geophone according to the travel-time first arrival data and the initial velocity model by using, for example, an FMA algorithm. The travel time first arrival data comprises the coordinates and the elevation of each shot point and the coordinates and the elevation of a corresponding demodulator probe of each shot point. The initial velocity model is formed by utilizing coordinates and elevation interpolation of all shot points and demodulator probes in travel-time first-arrival data of the target area, and has a certain grid number in the horizontal direction x and the vertical direction z for a two-dimensional initial velocity model and a certain grid number in the directions x and z and the direction y perpendicular to the xz plane for a three-dimensional initial velocity model. It will be appreciated that the travel-time first-arrival data for the target area is obtained after processing the raw seismic data for the target area.

After obtaining the original travel-time field data of each shot and the original travel-time field data of each geophone, the second obtaining unit 502 defines the original travel-time field data of each shot and the original travel-time field data of each geophone according to the first data defining rule, and can obtain the first defined travel-time field data of each shot and the first defined travel-time field data of each geophone. Wherein the first data definition rule may be a definition of original travel-time field data for each shot and each geophone point in an offset direction, an elevation direction, or a depth direction. If the original travel time field data is two-dimensional data, the offset direction refers to a horizontal direction, the height direction refers to an upward direction perpendicular to the horizontal direction, and the depth direction refers to a downward direction perpendicular to the horizontal direction. If the original travel time field data is three-dimensional data, the offset direction refers to an up-horizontal direction, the height direction refers to a direction perpendicular to a horizontal plane and the depth direction refers to a direction perpendicular to a horizontal plane and downward.

After obtaining the first defined travel-time field data of each shot and the first defined travel-time field data of each geophone, the third obtaining unit 503 may define the first defined travel-time field data of each shot and the first defined travel-time field data of each geophone according to a second data definition rule, and obtain the second defined travel-time field data of each shot and the second defined travel-time field data of each geophone. Wherein the second data definition rule is a definition of first defined travel-time-field data for each shot and each geophone point in the offset direction, the elevation direction, or the depth direction. It is to be understood that the second data definition rule is different from the first data definition rule, i.e. if the first data definition rule is a data definition in the offset direction, the second data definition rule is a data definition in the height direction or in the depth direction, if the first data definition rule is a data definition in the height direction, the second data definition rule is a data definition in the offset direction or in the depth direction; if the first data definition rule is a data definition in the depth direction, the second data definition rule is a data definition in the height direction or in the depth direction.

After obtaining the second defined travel-time field data for each shot and the second defined travel-time field data for each geophone point, the fourth obtaining unit 504 may define the second defined travel-time field data for each shot and the second defined travel-time field data for each geophone point according to a third data definition rule, and obtain final defined travel-time field data for each shot and final defined travel-time field data for each geophone point. Wherein the third data definition rule is a definition of first defined travel-time-field data for each shot and each geophone point in the offset direction, the elevation direction, or the depth direction. It is to be understood that the third data definition rule is different from both the first data definition rule and the second data definition rule, i.e. if the third data definition rule is a data definition in the depth direction, the first data definition rule and the second data definition rule are data definitions in the offset direction and in the height direction; if the third data definition rule is a data definition in an offset direction, the first data definition rule and the second data definition rule are data definitions in a depth direction and a height direction; if the third data defining rule is a data definition in the height direction, the first data defining rule and the second data defining rule are definitions in the offset direction and in the depth direction.

The fresnel-volume-based travel time field data processing device provided by the embodiment of the invention can obtain original travel time field data of each shot point and each demodulator probe according to travel time first arrival data and an initial speed model, then define the original travel time field data of each shot point and the original travel time field data of each demodulator probe according to a first data definition rule, obtain first defined travel time field data of each shot point and first defined travel time field data of each demodulator probe, further define the first defined travel time field data of each shot point and the first defined travel time field data of each demodulator probe according to a second data definition rule, obtain second defined travel time field data of each shot point and second defined travel time field data of each demodulator probe, and then define the second defined travel time field data of each shot point and the second defined travel time field data of each demodulator probe according to a third data definition rule, the final travel time field limiting data of each shot point and the final travel time field limiting data of each demodulator probe are obtained, the range of the travel time field data is narrowed, the number of Fresnel body calculation is reduced, and the processing efficiency of the travel time field data is improved.

Fig. 6 is a schematic structural diagram of a traveling time field data processing apparatus based on a fresnel body according to another embodiment of the present invention, and as shown in fig. 6, on the basis of the foregoing embodiments, further, the second obtaining unit 502 includes a first obtaining subunit 5021, a second obtaining subunit 5022, a third obtaining subunit 5023, a fourth obtaining subunit 5024, and a fifth obtaining subunit 5025, where:

the first obtaining subunit 5021 is configured to obtain a maximum offset distance of each shot point relative to a demodulator probe according to each shot point and the grid coordinates of each corresponding demodulator probe; the second obtaining subunit 5022 is configured to define the original travel-time field data of each shot according to the maximum offset distance of each shot relative to the demodulator probe, and obtain first defined travel-time field data of each shot; the third obtaining subunit 5023 is configured to obtain each shot point corresponding to each wave detection point according to each shot point and each corresponding wave detection point; the fourth obtaining subunit 5024 is configured to obtain a maximum offset distance of each demodulator probe relative to a shot point according to each demodulator probe and the corresponding grid coordinate of each shot point; a fifth obtaining subunit 5025 is configured to define the original travel-time field data for each of the demodulator probes according to the maximum offset distance of each of the demodulator probes relative to the shot, and obtain first defined travel-time field data for each of the demodulator probes.

Specifically, for the raw travel-time-field data for each shot, each shot and the corresponding respective demodulator probe are represented in grid coordinates. The first obtaining subunit 5021 may obtain an offset distance between each shot point and each corresponding demodulator probe according to the grid coordinates of each shot point and each corresponding demodulator probe, where the offset distance may be represented by a grid number. The first obtaining subunit 5021 compares the offset distances between each shot point and the corresponding respective demodulator probes to obtain a maximum offset distance therefrom, which is the maximum offset distance of each shot point relative to the demodulator probes, and there is one maximum offset distance for each shot point. Understandably, if the original travel-time field data is two-dimensional data, the offset distance between a shot point and a corresponding demodulator probe refers to the distance between the shot point and the corresponding demodulator probe in the horizontal direction; if the original travel time field data is three-dimensional data, the offset distance between the shot point and the corresponding demodulator probe refers to the projection of the straight-line distance between the shot point and the corresponding demodulator probe on the horizontal plane.

After obtaining the maximum offset distance of each shot from the demodulator probe, the second obtaining subunit 5022 defines the original travel time field data of each shot according to the maximum offset distance of each shot from the demodulator probe, that is, retains data of the original travel time field data of each shot, of which the offset distance from each shot is within the maximum offset distance range, discards data of each shot, of which the offset distance from each shot is greater than the maximum offset distance, and obtains the first defined travel time field data of each shot.

Since each shot and each corresponding demodulator probe are known, the third obtaining subunit 5023 may find each shot corresponding to each demodulator probe in each shot and each corresponding demodulator probe according to each demodulator probe.

The specific implementation process of the fourth obtaining subunit 5024, which is configured to obtain the maximum offset distance of each detector point relative to a shot point according to each detector point and the corresponding grid coordinates of each shot point, is similar to the specific process of the first obtaining subunit 5021, which is configured to obtain the maximum offset distance of each shot point relative to a detector point according to each shot point and the corresponding grid coordinates of each detector point, and is not described here in detail.

The fifth obtaining subunit 5025 is configured to define the original travel-time field data of each probe point according to the maximum offset distance of each probe point relative to a shot point, to obtain a specific implementation process of the first defined travel-time field data of each probe point, and is similar to the specific process of the second obtaining subunit 5022, which is configured to define the original travel-time field data of each shot point according to the maximum offset distance of each shot point relative to a probe point, to obtain the first defined travel-time field data of each shot point, and details of the process are omitted here.

Fig. 7 is a schematic structural diagram of a traveling time field data processing device based on a fresnel body according to yet another embodiment of the present invention, and as shown in fig. 7, on the basis of the foregoing embodiments, further, the third obtaining unit 503 includes a sixth obtaining sub-unit 5031 and a seventh obtaining sub-unit 5032, where:

the sixth obtaining subunit 5031 is configured to obtain an elevation-defined boundary according to the grid coordinates of each shot point and each geophone point in the vertical direction; the seventh obtaining subunit 5032 is configured to define the first defined travel-time field data of each shot according to the elevation defining boundary, obtain the second defined travel-time field data of each shot, and define the first defined travel-time field data of each geophone according to the elevation defining condition, and obtain the second defined travel-time field data of each geophone.

Specifically, the elevation of each shot in the travel-time first arrival data is converted into grid coordinates of each shot in the vertical direction, and the elevation of each geophone in the travel-time first arrival data is converted into grid coordinates of each geophone in the vertical direction. The sixth obtaining subunit 5031 may fit an elevation-defining boundary according to the grid coordinates in the vertical direction of each shot point and each geophone point. And if the travel time first arrival data are two-dimensional data, fitting an elevation limit boundary into a curve. And if the travel first arrival data are three-dimensional data, fitting the elevation limit boundary into a curved surface.

Specifically, after obtaining the elevation-defined boundary, the seventh obtaining subunit 5032 may define the first defined travel-time field data of each shot according to the elevation-defined boundary, that is, truncate the data above the elevation-defined boundary in the first defined travel-time field data of each shot, and obtain the second defined travel-time field data of each shot. The server can also limit the first limited travel-time field data of each wave detection point according to the elevation limit boundary, namely, the data above the elevation limit boundary in the first limited travel-time field data of each wave detection point is omitted, and the second limited travel-time field data of each wave detection point is obtained.

Fig. 8 is a schematic structural diagram of a fresnel-volume-based travel time field data processing apparatus according to still another embodiment of the present invention, and as shown in fig. 8, on the basis of the above embodiments, further, the fourth obtaining unit 504 includes an eighth obtaining subunit 5041, a ninth obtaining subunit 5042, a tenth obtaining subunit 5043, and an eleventh obtaining subunit 5044, where:

the eighth obtaining subunit 5041 is configured to obtain travel time field data of the shot-survey pair with the maximum offset distance for each shot according to the second limited travel time field data of each shot and the second limited travel time field data of the corresponding demodulator probe with the maximum offset distance for each shot; the ninth obtaining subunit 5042 is configured to obtain fresnel volume path information of the maximum offset distance geophone pair for each shot point according to the travel time field data of the maximum offset distance geophone pair for each shot point; the tenth obtaining subunit 5043 is configured to obtain, according to the fresnel volume path information of the maximum offset distance shot-check pair of each shot point, a maximum fresnel volume depth of each shot point; the eleventh obtaining subunit 5044 is configured to define the second limited travel time field data of each shot according to the maximum fresnel volume depth of each shot, obtain the final limited travel time field data of each shot, and define the second limited travel time field data of each geophone corresponding to each shot according to the maximum fresnel volume depth of each shot, so as to obtain the final limited travel time field data of each geophone corresponding to each shot.

Specifically, there is an offset distance between each shot point and each corresponding geophone point, and the eighth obtaining subunit 5041 may obtain, from the offset distances between each shot point and the corresponding respective geophone points, a maximum offset distance of each shot point with respect to the geophone point, that is, a maximum offset distance of each shot point. The eighth obtaining subunit 5041 may obtain the travel-time field data of the shot-survey pair with the maximum offset distance for each shot by superimposing the second limited travel-time field data of each shot and the second limited travel-time field data of the corresponding demodulator probe with the maximum offset distance for each shot. And each shot point and the corresponding demodulator probe of the maximum offset distance of each shot point form a shot detection pair of the maximum offset distance of each shot point.

After obtaining the travel-time field data of the maximum offset shot pair for each shot, the ninth obtaining subunit 5042 may calculate fresnel volume path information of the maximum offset shot pair for each shot according to the travel-time field data of the maximum offset shot pair for each shot and an empirical formula.

After obtaining the fresnel volume path information of the maximum offset distance shot pair for each shot point, the tenth obtaining subunit 5043 may calculate and obtain the maximum fresnel volume depth for each shot point according to the fresnel volume path information of the maximum offset distance shot pair for each shot point.

After obtaining the maximum fresnel body depth of each shot point, the eleventh obtaining subunit 5044 may define the second limited travel-time field data of each shot point according to the maximum fresnel body depth of each shot point, that is, retain data in the range of the maximum fresnel body depth of each shot point in the second limited travel-time field data of each shot point, discard data exceeding the range of the maximum fresnel body depth of each shot point in the second limited travel-time field data of each shot point, and obtain the final limited travel-time field data of each shot point. The eleventh obtaining subunit 5044 may further limit the second limited travel time field data of each demodulator probe corresponding to each shot point according to the maximum fresnel body depth of each shot point, that is, data in the maximum fresnel body depth range of each shot point in the second limited travel time field data of each demodulator probe corresponding to each shot point is retained, data exceeding the maximum fresnel body depth range of each shot point in the second limited travel time field data of each demodulator probe corresponding to each shot point is discarded, and the final limited travel time field data of each demodulator probe corresponding to each shot point is obtained.

The embodiment of the apparatus provided in the embodiment of the present invention may be specifically configured to execute the processing flows of the above method embodiments, and the functions of the apparatus are not described herein again, and refer to the detailed description of the above method embodiments.

Fig. 9 is a schematic physical structure diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 9, the electronic device may include: a processor (processor)901, a communication Interface (Communications Interface)902, a memory (memory)903 and a communication bus 904, wherein the processor 901, the communication Interface 902 and the memory 903 are communicated with each other through the communication bus 904. The processor 901 may call logic instructions in the memory 903 to perform the following method: acquiring original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model; defining the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe according to a first data definition rule to obtain first defined travel-time field data of each shot point and first defined travel-time field data of each demodulator probe; defining the first defined travel-time field data of each shot point and the first defined travel-time field data of each demodulator probe according to a second data definition rule to obtain second defined travel-time field data of each shot point and second defined travel-time field data of each demodulator probe; and defining the second defined travel-time field data of each shot point and the second defined travel-time field data of each demodulator probe according to a third data definition rule, and obtaining the final defined travel-time field data of each shot point and the final defined travel-time field data of each demodulator probe.

In addition, the logic instructions in the memory 903 may be implemented in a software functional unit and stored in a computer readable storage medium when the logic instructions are sold or used as a separate product. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.

The present embodiment discloses a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions which, when executed by a computer, enable the computer to perform the method provided by the above-mentioned method embodiments, for example, comprising: acquiring original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model; defining the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe according to a first data definition rule to obtain first defined travel-time field data of each shot point and first defined travel-time field data of each demodulator probe; defining the first defined travel-time field data of each shot point and the first defined travel-time field data of each demodulator probe according to a second data definition rule to obtain second defined travel-time field data of each shot point and second defined travel-time field data of each demodulator probe; and defining the second defined travel-time field data of each shot point and the second defined travel-time field data of each demodulator probe according to a third data definition rule, and obtaining the final defined travel-time field data of each shot point and the final defined travel-time field data of each demodulator probe.

The present embodiment provides a computer-readable storage medium, which stores a computer program, where the computer program causes the computer to execute the method provided by the above method embodiments, for example, the method includes: acquiring original travel time field data of each shot point and each demodulator probe according to the travel time first arrival data and the initial speed model; defining the original travel-time field data of each shot point and the original travel-time field data of each demodulator probe according to a first data definition rule to obtain first defined travel-time field data of each shot point and first defined travel-time field data of each demodulator probe; defining the first defined travel-time field data of each shot point and the first defined travel-time field data of each demodulator probe according to a second data definition rule to obtain second defined travel-time field data of each shot point and second defined travel-time field data of each demodulator probe; and defining the second defined travel-time field data of each shot point and the second defined travel-time field data of each demodulator probe according to a third data definition rule, and obtaining the final defined travel-time field data of each shot point and the final defined travel-time field data of each demodulator probe.

As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

In the description herein, reference to the description of the terms "one embodiment," "a particular embodiment," "some embodiments," "for example," "an example," "a particular example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.