Method and device for determining seismic wave travel time, computer equipment and storage medium

Abstract

The application provides a method and device for determining seismic travel time, computer equipment and storage medium, and belongs to the technical field of seismic exploration. The method comprises the following steps: obtaining exploration data in a work area to be determined, wherein the exploration data comprises a plurality of seismic trace data; determining the seismic travel time corresponding to a plurality of simulation shot points on a preset grid in the work area according to the plurality of seismic trace data; determining the zero-offset travel time corresponding to a plurality of imaging points in the work area according to the seismic travel time corresponding to the plurality of simulation shot points; and determining the seismic travel time corresponding to a plurality of shot points and a plurality of geophones according to the plurality of seismic trace data, the seismic travel time corresponding to the plurality of simulation shot points and the zero-offset travel time corresponding to the plurality of imaging points. Since the zero-offset seismic travel time corresponds to the travel time when the planar coordinates of the imaging point and the simulation shot point are the same, the seismic travel time of the exploration point is prevented from exceeding the accuracy range, and the accuracy of the determined seismic travel time is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of seismic exploration, and particularly relates to a method and device for determining seismic wave travel time, computer equipment and a storage medium. BACKGROUND

[0002] At present, seismic wave migration imaging technology is an important technology for determining the geological structure of a reservoir. The seismic wave migration imaging technology determines the stacking profile of the reservoir according to the travel time of the seismic wave at the positions of the exploration points such as the shot points and the receiver points, and the geological structure of the reservoir is represented by the stacking profile, and therefore, the travel time of the seismic wave at the positions of the exploration points needs to be determined.

[0003] In the related art, the travel time of the seismic wave at the positions of the exploration points is determined by a travel time interpolation method. The method includes: determining the travel time of a plurality of simulation shot points on a preset grid; determining the grid area in which the exploration point is located according to the position of the exploration point; and determining the travel time of the exploration point according to the positions and the travel times of the four simulation shot points that constitute the grid area.

[0004] However, in the above related art, when the distances between the exploration points and the imaging points are all less than the distances from the four simulation shot points to the imaging points, the accuracy range of the travel time interpolation method is exceeded, and the accuracy of the determined travel time of the exploration point is low. SUMMARY

[0005] Embodiments of the present application provide a method and device for determining seismic wave travel time, computer equipment and a storage medium, which can improve the accuracy of drilling compliance information in the exploration area. The technical solution is as follows:

[0006] In one aspect, the present application provides a method for determining seismic wave travel time, which includes:

[0007] Obtaining exploration data in a work area to be determined, the exploration data including a plurality of seismic trace data, one seismic trace data corresponding to one shot point and one receiver point among a plurality of shot points and a plurality of receiver points arranged in the work area;

[0008] Determining the seismic wave travel time corresponding to a plurality of simulation shot points on a preset grid in the work area according to the plurality of seismic trace data, the preset grid including a plurality of grid areas, and the seismic wave travel time being used to represent the propagation time of the seismic wave signal from the simulation shot point to a plurality of imaging points in the work area;

[0009] Determining the zero offset travel time corresponding to a plurality of imaging points in the work area according to the seismic wave travel time corresponding to the plurality of simulation shot points, the zero offset travel time being used to represent the propagation time of the seismic wave signal from the simulation shot point to the imaging point when the planar coordinates of the imaging point and the simulation shot point are the same;

[0010] According to the plurality of seismic trace data, the seismic wave travel times corresponding to the plurality of simulated shot points, and the zero-offset travel times corresponding to the plurality of imaging points, the seismic wave travel times corresponding to the plurality of shot points and the plurality of geophones are determined.

[0011] In a possible implementation, the determining, according to the plurality of seismic trace data, the seismic wave travel times corresponding to the plurality of simulated shot points, and the zero-offset travel times corresponding to the plurality of imaging points, the seismic wave travel times corresponding to the plurality of shot points and the plurality of geophones comprises:

[0012] For each seismic trace data, a first position of a shot point corresponding to the seismic trace data and a second position of a geophone corresponding to the seismic trace data are determined.

[0013] Based on the first position, the seismic wave travel times corresponding to the plurality of simulated shot points, and the zero-offset travel times corresponding to the plurality of imaging points, a travel time of the shot point corresponding to the seismic trace data to the plurality of imaging points is determined, to obtain the seismic wave travel time corresponding to the shot point, and based on the second position, the travel times corresponding to the plurality of simulated shot points, and the zero-offset travel times corresponding to the plurality of imaging points, a travel time of the geophone corresponding to the seismic trace data to the plurality of imaging points is determined, to obtain the seismic wave travel time corresponding to the geophone.

[0014] In another possible implementation, the determining, based on the first position, the seismic wave travel times corresponding to the plurality of simulated shot points, and the zero-offset travel times corresponding to the plurality of imaging points, the travel time of the shot point corresponding to the seismic trace data to the plurality of imaging points comprises:

[0015] Based on the first position, a first target grid region in which the shot point corresponding to the seismic trace data is located is determined, four first simulated shot point corresponding seismic wave travel times constituting a first target grid region of the shot point corresponding to the seismic trace data are determined from the seismic wave travel times corresponding to the plurality of simulated shot points, and the first simulated shot point corresponding seismic wave travel time is used to represent a travel time of a seismic wave signal from the first simulated shot point to the plurality of imaging points.

[0016] For each imaging point, a target zero-offset travel time corresponding to the imaging point is determined from the zero-offset travel times corresponding to the plurality of imaging points, and for each first simulated shot point, a target travel time of the first simulated shot point to the imaging point is determined from the first simulated shot point corresponding seismic wave travel time, to obtain target travel times of the four first simulated shot points.

[0017] Based on the first position, the target zero-offset travel time corresponding to the imaging point, and the target travel times of the four first simulated shot points, the travel time of the shot point corresponding to the seismic trace data to the imaging point is determined, to obtain the travel time of the shot point to the plurality of imaging points.

[0018] In another possible implementation, the determining of the travel time of the shot point corresponding to the seismic trace data to the imaging point based on the first position, the target zero offset travel time corresponding to the imaging point and the target travel times of the four first simulated shot points comprises:

[0019] determining positions of the imaging point and the four first simulated shot points, determining four triangular regions based on the positions of the imaging point and the four first simulated shot points, and one triangular region including one imaging point and two adjacent first simulated shot points;

[0020] determining a target triangular region in which the first position is located from the four triangular regions based on the first position, the position of the imaging point and the positions of the four first simulated shot points;

[0021] determining the travel time of the shot point corresponding to the seismic trace data to the imaging point based on the target travel times of the two adjacent first simulated shot points corresponding to the target triangular region and the target zero offset travel time corresponding to the imaging point.

[0022] In another possible implementation, the determining of the zero offset travel times corresponding to the multiple imaging points in the work area according to the seismic wave travel times corresponding to the multiple simulated shot points comprises:

[0023] determining the zero offset travel times corresponding to the multiple simulated shot points from the seismic wave travel times of the multiple simulated shot points;

[0024] for each imaging point, determining a second target grid region in which the imaging point is located, and determining the zero offset travel times corresponding to four second simulated shot points constituting the second target grid region from the zero offset travel times corresponding to the multiple simulated shot points;

[0025] determining the zero offset travel time of the imaging point based on the zero offset travel times corresponding to the four second simulated shot points, to obtain the zero offset travel times corresponding to the multiple imaging points in the work area.

[0026] In another possible implementation, the determining of the seismic wave travel times corresponding to the multiple simulated shot points on the preset grid in the work area according to the multiple seismic trace data comprises:

[0027] determining a velocity parameter of seismic wave signal propagation in a reservoir in the work area and an anisotropy parameter of the reservoir according to the multiple seismic trace data;

[0028] determining the seismic wave travel times corresponding to the multiple simulated shot points on the preset grid in the work area according to the velocity parameter and the anisotropy parameter.

[0029] In another possible implementation, the method further includes:

[0030] Based on the seismic wave travel times corresponding to the multiple shot points and the multiple receiver points, the migration imaging information of the work area is determined.

[0031] On the other hand, this application provides a device for determining the travel time of seismic waves, the device comprising:

[0032] The acquisition module is used to acquire exploration data within the work area to be determined. The exploration data includes multiple seismic traces, and one seismic trace corresponds to one of the multiple shot points and one of the multiple receiver points deployed within the work area.

[0033] The first determining module is used to determine the seismic wave travel time corresponding to multiple simulated shot points on a preset grid within the work area based on the multiple seismic trace data. The preset grid includes multiple grid regions, and the seismic wave travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to multiple imaging points within the work area.

[0034] The second determining module is used to determine the zero-shot distance travel time corresponding to multiple imaging points in the work area based on the seismic wave travel time corresponding to the multiple simulated shot points. The zero-shot distance travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to the imaging point when the plane coordinates of the imaging point and the simulated shot point are the same.

[0035] The third determining module is used to determine the seismic wave travel times corresponding to the multiple shot points and the multiple receiver points based on the multiple seismic trace data, the seismic wave travel times corresponding to the multiple simulated shot points, and the zero-shot distance travel times corresponding to the multiple imaging points.

[0036] In one possible implementation, the third determining module is configured to, for each seismic trace data, determine a first position of the shot point and a second position of the receiver point corresponding to the seismic trace data; based on the first position, the seismic wave travel time corresponding to the plurality of simulated shot points and the zero-shot distance travel time corresponding to the plurality of imaging points, determine the travel time from the shot point corresponding to the seismic trace data to the plurality of imaging points, thereby obtaining the seismic wave travel time corresponding to the shot point; and based on the second position, the travel time corresponding to the plurality of simulated shot points and the zero-shot distance travel time corresponding to the plurality of imaging points, determine the travel time from the receiver point corresponding to the seismic trace data to the plurality of imaging points, thereby obtaining the seismic wave travel time corresponding to the receiver point.

[0037] In another possible implementation, the third determining module is configured to, based on the first location, determine the first target grid region where the shot point corresponding to the seismic trace data is located; determine the seismic wave travel times corresponding to the four first simulated shot points constituting the first target grid region from the seismic wave travel times corresponding to the plurality of simulated shot points, wherein the seismic wave travel times corresponding to the first simulated shot points are used to represent the travel time of the seismic wave signal from the first simulated shot points to the plurality of imaging points; for each imaging point, determine the target zero-shot distance travel time corresponding to the imaging point from the zero-shot distance travel times corresponding to the plurality of imaging points; and for each first simulated shot point, determine the target travel time from the first simulated shot point to the imaging point from the seismic wave travel times corresponding to the first simulated shot point, thereby obtaining the target travel time of the four first simulated shot points; and, based on the first location, the target zero-shot distance travel time corresponding to the imaging point, and the target travel time of the four first simulated shot points, determine the travel time from the shot point corresponding to the seismic trace data to the imaging point, thereby obtaining the travel time from the shot point to the plurality of imaging points.

[0038] In another possible implementation, the third determining module is used to determine the position of the imaging point and the positions of the four first simulated shot points, and based on the positions of the imaging point and the four first simulated shot points, determine four triangular regions, each triangular region including one imaging point and two adjacent first simulated shot points; based on the first position, the position of the imaging point, and the positions of the four first simulated shot points, determine the target triangular region where the first position is located from the four triangular regions; based on the target travel time of the two adjacent first simulated shot points corresponding to the target triangular region and the target zero-shot distance travel time corresponding to the imaging point, determine the travel time from the shot point corresponding to the seismic trace data to the imaging point.

[0039] In another possible implementation, the second determining module is configured to determine the zero-range travel time corresponding to the plurality of simulated shot points from the seismic wave travel times of the plurality of simulated shot points; for each imaging point, determine the second target grid region where the imaging point is located, and determine the zero-range travel time corresponding to the four second simulated shot points constituting the second target grid region from the zero-range travel times corresponding to the plurality of simulated shot points; based on the zero-range travel times corresponding to the four second simulated shot points, determine the zero-range travel time of the imaging point, thereby obtaining the zero-range travel time corresponding to the plurality of imaging points within the work area.

[0040] In another possible implementation, the first determining module is used to determine the velocity parameters of the seismic wave signal propagating in the reservoir of the work area and the anisotropy parameters of the reservoir based on the multiple seismic trace data; and to determine the seismic wave travel time corresponding to multiple simulated shot points on a preset grid in the work area based on the velocity parameters and the anisotropy parameters.

[0041] In another possible implementation, the device further includes:

[0042] The fourth determining module is used to determine the migration imaging information of the work area based on the seismic wave travel time corresponding to the plurality of shot points and the plurality of receiver points.

[0043] On the other hand, embodiments of this application provide a computer device, the computer device including: a processor and a memory, the memory storing at least one piece of program code, the at least one piece of program code being loaded and executed by the processor to implement the operations performed in the method for determining seismic wave travel time as described in any of the above possible implementations.

[0044] On the other hand, embodiments of this application provide a computer-readable storage medium storing at least one piece of program code, which is loaded and executed by a processor to perform the operations performed in the method for determining seismic wave travel time described in any of the above possible implementations.

[0045] On the other hand, embodiments of this application provide a computer program product comprising at least one piece of program code, which is loaded and executed by a processor to implement the operations performed in the method for determining seismic wave travel time described in any of the possible implementations above.

[0046] The beneficial effects of the technical solutions provided in this application include at least the following:

[0047] This application provides a method for determining seismic wave travel time. By using the seismic wave travel times corresponding to multiple simulated shot points and the zero-range travel times corresponding to multiple imaging points, the seismic wave travel times corresponding to multiple shot points and multiple receiver points are determined. The zero-range travel time is the travel time corresponding to the same planar coordinates between the imaging point and the simulated shot point. At this time, the distance between the imaging point and the simulated shot point is less than the distance from the shot point and receiver point to the imaging point. This avoids the seismic wave travel time of the exploration point from exceeding the accuracy range of the travel time interpolation method, thus improving the accuracy of the travel time of the determined exploration point location. Attached Figure Description

[0048] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0049] Figure 1 This is a flowchart illustrating a method for determining the travel time of seismic waves according to an exemplary embodiment;

[0050] Figure 2 This is a seismic profile determined by the seismic wave travel time determined by this application, according to an exemplary embodiment.

[0051] Figure 3 This is a seismic profile determined by seismic wave travel time determined by the prior art, according to an exemplary embodiment.

[0052] Figure 4 This is a block diagram illustrating a device for determining the travel time of seismic waves according to an exemplary embodiment;

[0053] Figure 5 This is a structural block diagram of a computer device according to an exemplary embodiment. Detailed Implementation

[0054] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0055] Figure 1 This is a flowchart illustrating a method for determining the travel time of seismic waves according to an exemplary embodiment, executed by a computer device. See also Figure 1 The method includes:

[0056] 101. Computer equipment acquires exploration data within the work area to be determined. The exploration data includes multiple seismic traces. One seismic trace corresponds to one of the multiple shot points and one of the multiple receiver points deployed within the work area.

[0057] In one possible implementation, multiple shot points and multiple geophones are deployed within the work area to be determined; shot points are used to transmit seismic wave signals, and geophones are used to receive seismic wave signals. Accordingly, this step involves: multiple geophones receiving seismic wave signals transmitted from multiple shot points, obtaining multiple seismic trace data, where each seismic trace data includes a seismic wave signal received by one geophone from one shot point; the multiple seismic trace data are uploaded to a computer device, which then acquires the multiple seismic trace data for the work area, obtaining exploration data corresponding to the multiple shot points and multiple geophones. The seismic trace data includes changes in seismic wave signals within a preset time period.

[0058] In this embodiment, the preset duration is not specifically limited and can be set as needed. Optionally, the preset duration is any value between 1s and 20s, for example, a preset duration of 6s, 8s, 10s, etc.

[0059] 102. The computer equipment determines the seismic wave travel time corresponding to multiple simulated shot points on a preset grid within the work area based on multiple seismic trace data. The preset grid includes multiple grid regions, and each grid region is a quadrilateral region composed of four adjacent simulated shot points. The seismic wave travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to multiple imaging points within the work area.

[0060] In this step, the distance between any two adjacent simulated shot points can be the same or different. In one possible implementation, the grid area is rectangular, the quadrilateral area formed by four adjacent simulated shot points is rectangular, the distance between any two adjacent simulated shot points is the same, and this distance is a first preset distance. In this embodiment, the value of the first preset distance is not specifically limited and can be set and modified as needed. Optionally, the first preset distance is any value between 10m and 300m, for example, the first preset distance is 10m, 50m, 100m, etc.

[0061] In one possible implementation, the computer equipment determines the seismic wave travel time corresponding to multiple simulated shot points based on exploration data received from actual shot points. Correspondingly, the step of the computer equipment determining the seismic wave travel time corresponding to multiple simulated shot points on a preset grid within the work area based on multiple seismic trace data is as follows: the computer equipment determines the velocity parameters of the seismic wave signal propagating within the reservoir of the work area and the anisotropy parameters of the reservoir based on the multiple seismic trace data; based on the velocity parameters and anisotropy parameters, it determines the seismic wave travel time corresponding to multiple simulated shot points on a preset grid within the work area.

[0062] It should be noted that the imaging point can be any location within the reservoir. In one possible implementation, the distance between any two adjacent imaging points is the same, and this distance is a second preset distance. In the embodiments of this application, the value of the second preset distance is not specifically limited and can be set and modified as needed. Optionally, the second preset distance is any value between 10m and 100m, for example, the second preset distance is 10m, 50m, 80m, etc.

[0063] The density of the grid formed by multiple imaging points can be greater than or less than the preset grid density. Correspondingly, when the distance between any two adjacent imaging points is less than the distance between any two adjacent simulated shot points, the density of the grid formed by the imaging points is greater than the preset grid density; when the distance between any two adjacent imaging points is greater than the distance between any two adjacent simulated shot points, the density of the grid formed by multiple imaging points is less than the preset grid density.

[0064] 103. The computer equipment determines the zero-shot distance travel time corresponding to multiple imaging points in the work area based on the seismic wave travel time of multiple simulated shot points. The zero-shot distance travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to the imaging point when the plane coordinates of the imaging point and the simulated shot point are the same.

[0065] In one possible implementation, the computer device determines the zero-range travel time corresponding to multiple imaging points based on the zero-range travel time of multiple simulated shot points. Accordingly, the computer device determining the zero-range travel time corresponding to multiple imaging points within the work area may include the following steps (11) to (13):

[0066] (11) The computer equipment determines the zero-distance travel time corresponding to multiple simulated shot points from the seismic wave travel time of multiple simulated shot points.

[0067] In one possible implementation, this step is as follows: For each simulated shot point, the computer device determines a target imaging point with the same location as the simulated shot point from among the multiple imaging points, based on the location of the simulated shot point and the locations of multiple imaging points; the travel time from the simulated shot point to the target imaging point is selected from the seismic wave travel time of the simulated shot point to obtain the zero-range travel time corresponding to the simulated shot point.

[0068] (12) For each imaging point, the computer device determines the second target grid region where the imaging point is located, and determines the zero-distance travel time of the four second simulated gun points that make up the second target grid region from the zero-distance travel time of the multiple simulated gun points.

[0069] In one possible implementation, the computer device stores the coordinate ranges of multiple grid regions included in a preset grid. Accordingly, the step of the computer device determining the second target grid region where each imaging point is located is as follows: for each imaging point, the computer device acquires the coordinates of the imaging point; based on the coordinates of the imaging point, determines the coordinate range in which the coordinates lie; and determines the second target grid region corresponding to the coordinate range from the coordinate ranges of the multiple grid regions. Optionally, the coordinates of the imaging point are planar coordinates; for example, the coordinates of the imaging point can be represented by (x3, y3), where x3 represents the planar abscissa of the imaging point, and y3 represents the planar ordinate of the imaging point.

[0070] In one possible implementation, a grid region is a quadrilateral region composed of four adjacent simulated shot points, and the computer device stores the correspondence between the grid region and the four adjacent simulated shot points. Accordingly, the step of the computer device determining the zero-distance travel time corresponding to the four second simulated shot points constituting the second target grid region is as follows: the computer device determines the four second simulated shot points corresponding to the second target grid region, and determines the zero-distance travel time corresponding to the four second simulated shot points from the zero-distance travel times corresponding to multiple simulated shot points.

[0071] It should be noted that the zero-distance travel time is the travel time from the second simulated shot point to the imaging point when the planar coordinates of the second simulated shot point and the imaging point are the same. The imaging points with the same planar coordinates as the second simulated shot point include multiple imaging points with different depths but the same planar coordinates; correspondingly, the zero-distance travel time includes the travel time from the second simulated shot point to multiple imaging points with different depths but the same planar coordinates.

[0072] (13) The computer equipment determines the zero-range travel time of the imaging point based on the zero-range travel time corresponding to the four second simulated shot points, and obtains the zero-range travel time corresponding to multiple imaging points in the work area.

[0073] In one possible implementation, the computer device determines the zero-range travel time of the imaging point using bilinear interpolation. Accordingly, this step involves: the computer device determining the zero-range travel time of the imaging point along the horizontal axis based on the zero-range travel times corresponding to the four second simulated shot points, the x-coordinates of the four second simulated shot points, and the x-coordinate of the imaging point; then, based on the zero-range travel time of the imaging point along the horizontal axis, the y-coordinates of the four second simulated shot points, and the x and y coordinates of the imaging point, the computer device determines the zero-range travel time of the imaging point.

[0074] 104. The computer equipment determines the seismic wave travel times corresponding to multiple shot points and multiple receiver points based on multiple seismic trace data, seismic wave travel times of multiple simulated shot points, and zero-shot distance travel times of multiple imaging points.

[0075] In one possible implementation, this step may include the following steps (21) to (23):

[0076] (21) For each seismic trace data, the computer equipment determines the first position of the shot point and the second position of the receiver point corresponding to the seismic trace data.

[0077] In one possible implementation, a seismic trace corresponds to a seismic wave signal emitted from a shot point and received by a geophone. The computer device stores the correspondence between the seismic trace data and the shot point and geophone locations. Accordingly, this step is as follows: for each seismic trace, the computer device determines the first position of the shot point and the second position of the geophone corresponding to that seismic trace data from the stored correspondence between the seismic trace data and the shot point and geophone locations.

[0078] (22) The computer device determines the travel time from the shot point to the multiple imaging points corresponding to the seismic trace data based on the first position, the seismic wave travel time corresponding to the multiple simulated shot points and the zero shot distance travel time corresponding to the multiple imaging points, and obtains the seismic wave travel time corresponding to the shot point.

[0079] In one possible implementation, this step can be achieved through the following steps (a) to (c):

[0080] (a) The computer device determines the first target grid region where the shot point corresponding to the seismic trace data is located based on the first location, and determines the seismic wave travel time corresponding to the four first simulated shot points that make up the first target grid region from the seismic wave travel time corresponding to multiple simulated shot points. The seismic wave travel time corresponding to the first simulated shot points is used to represent the travel time of the seismic wave signal from the first simulated shot point to multiple imaging points.

[0081] In one possible implementation, the computer device stores the coordinate ranges of multiple grid regions included in a preset grid. Accordingly, the step of the computer device determining the first target grid region where the shot point corresponding to the seismic trace data is located based on the first position is as follows: the computer device determines the coordinates of the shot point corresponding to the seismic trace data based on the first position; determines the first target grid region corresponding to the coordinate range of the shot point corresponding to the seismic trace data; the computer device determines four first simulated shot points corresponding to the first target grid region, and determines the seismic wave travel times corresponding to the four first simulated shot points from the seismic wave travel times corresponding to the multiple simulated shot points.

[0082] It should be noted that the coordinates of the shot points corresponding to the seismic trace data are planar coordinates. For example, the coordinates of the shot points corresponding to the seismic trace data can be represented by (x, y), where x represents the planar abscissa of the shot point corresponding to the seismic trace data, and y represents the planar ordinate of the shot point corresponding to the seismic trace data.

[0083] (b) For each imaging point, the computer device determines the target zero-range travel time corresponding to the imaging point from the zero-range travel time corresponding to multiple imaging points, and for each first simulated shot point, determines the target travel time from the first simulated shot point to the imaging point from the seismic wave travel time corresponding to the first simulated shot point, thus obtaining the target travel time of four first simulated shot points.

[0084] In one possible implementation, the target exploration location comprises multiple surface elements, with each surface element corresponding to an imaging point location. The computer stores the correspondence between the locations of multiple imaging points and their zero-range travel times. Accordingly, this step involves the computer determining the target zero-range travel time corresponding to each imaging point based on the target location of that imaging point and the correspondence between the location and the zero-range travel time.

[0085] It should be noted that the seismic wave travel time corresponding to the first simulated shot point includes the travel time from the first simulated shot point to multiple imaging points. For each imaging point, the target travel time is used to represent the travel time from the first simulated shot point to that imaging point. In one possible implementation, the computer stores the correspondence between the locations of the simulated shot points, the locations of the imaging points, and their travel times. Accordingly, the step of the computer device determining the target travel time corresponding to the four first simulated shot points from the travel times corresponding to the four first simulated shot points is as follows: For each first simulated shot point, based on the location of the first simulated shot point and the location of the imaging point, the computer device determines the target travel time corresponding to the first simulated shot point from the correspondence between the locations of the simulated shot points, the locations of the imaging points, and their travel times, thus obtaining the target travel times corresponding to the four first simulated shot points.

[0086] (c) The computer device determines the travel time from the shot point to the imaging point corresponding to the first position, the target zero-shot distance travel time corresponding to the imaging point, and the target travel time of the four first simulated shot points, and obtains the travel time from the shot point to multiple imaging points.

[0087] In one possible implementation, this step is as follows: the computer device determines the location of the imaging point and the locations of four first simulated shot points; based on the location of the imaging point and the locations of the four first simulated shot points, four triangular regions are determined, each triangular region including one imaging point and two adjacent first simulated shot points; based on the first location, the location of the imaging point, and the locations of the four first simulated shot points, the target triangular region where the first location is located is determined from the four triangular regions; based on the target travel time of the two adjacent first simulated shot points corresponding to the target triangular region and the target zero-shot distance travel time corresponding to the imaging point, the travel time from the shot point to the imaging point corresponding to the seismic trace data is determined; the travel time from the shot point to multiple imaging points is obtained.

[0088] It should be noted that the imaging point can be located within the first grid area formed by the four first simulated shot points, or it can be located outside the first grid area. The first grid area includes four sides, and the imaging point and one side form a triangular area, resulting in four triangular areas.

[0089] Optionally, the step of the computer device determining the target triangular region where the first position is located from the four triangular regions is as follows: For each triangular region, the computer device determines the target parameters according to the first position, the position of the imaging point, and the positions of the two adjacent first simulated shot points corresponding to the triangular region, using the following formula 1; when the target parameters are positive, the shot point corresponding to the seismic trace data is determined to be within the triangular region, and the triangular region is determined to be the target triangular region where the first position is located.

[0090] Formula 1: (a1,b1,c1)=(x-x1,y-y1,0)×(x-x2,y-y2,0)

[0091] (a2,b2,c2)=(x-x2,y-y2,0)×(x-x3,y-y3,0)

[0092] (a3,b3,c3)=(x-x3,y-y3,0)×(x-x1,y-y1,0)

[0093] A = c1 * c2 * c3

[0094] Where A represents the target parameters, (x,y) represents the first position, (x1,y1) and (x2,y2) represent the positions of two adjacent first simulated shot points, and (x3,y3) represents the position of the imaging point.

[0095] It should be noted that since the shot points corresponding to the seismic trace data are located within a grid area composed of four first simulated shot points, and the four triangular areas are obtained by dividing the grid area, the shot points corresponding to the seismic trace data are within a triangular area. That is to say, only one target triangular area has a positive value for A, while the other three triangular areas have negative values ​​for A.

[0096] In this embodiment, since the target triangular region where the shot point corresponding to the seismic trace data is located is determined by the location of the imaging point and the location of the first simulated shot point, the distance from the first simulated shot point to the imaging point is greater than the distance from the first location to the imaging point. This ensures that the accuracy range corresponding to the target triangular region is greater than the seismic wave travel time of the shot point corresponding to the first location, thus improving the accuracy of the determined seismic wave travel time.

[0097] Optionally, the computer equipment determines the seismic wave travel time of the shot point corresponding to the seismic trace data by the following steps: Based on the target travel time of two adjacent first simulated shot points corresponding to the target triangular region and the target zero-shot distance travel time corresponding to the imaging point, the computer equipment determines the seismic wave travel time of the shot point corresponding to the seismic trace data by the following formula 2.

[0098] Formula 2: t = t1*c1 + t2*c2 + t3*c3

[0099] c1=(ys3*x 23 -ys3*x 23 ) / (y 13 *x 23 -x 13 *y 23 )

[0100] c2=(xs3*y 13 -ys3*x 13 ) / (y 13 *x 23 -x 13 *y 23 )

[0101] c3 = 1 - c1 - c2

[0102] Where t represents the seismic wave travel time of the shot point corresponding to the seismic trace data, t1 and t2 represent the target travel time of the two adjacent first simulated shot points corresponding to the target triangular region, and t3 represents the target zero-shot distance travel time corresponding to the imaging point.

[0103] In this embodiment of the application, since the seismic wave travel time of the shot point corresponding to the seismic trace data is determined based on the target travel time of the first simulated shot point and the target zero-range travel time of the imaging point, and the distance from the first simulated shot point to the imaging point is greater than the distance from the shot point to the imaging point corresponding to the seismic trace data, it is ensured that the target travel time of the first simulated shot point covers the seismic wave travel time of the shot point, thereby improving the accuracy of the determined seismic wave travel time.

[0104] (23) The computer equipment determines the travel time from the receiver point to the multiple imaging points based on the second position, the travel time corresponding to multiple simulated shot points and the zero-shot distance travel time corresponding to multiple imaging points, and obtains the seismic wave travel time corresponding to the receiver point.

[0105] It should be noted that the geophones corresponding to the seismic trace data and the shot points corresponding to the seismic trace data are only different in location. The method for computer equipment to determine the travel time corresponding to multiple geophones is the same as the method for determining the travel time corresponding to multiple shot points, and will not be elaborated here.

[0106] Another point to note is that after determining the travel times corresponding to multiple shot points and multiple receiver points, the computer equipment determines the migration imaging information of the work area based on the seismic wave travel times of the multiple shot points and multiple receiver points. In one possible implementation, the step of the computer equipment determining the migration imaging information of the work area based on the seismic wave travel times of multiple shot points and multiple receiver points is as follows: the computer equipment performs superimposed imaging on the seismic wave travel times of multiple shot points and multiple receiver points using the seismic wave travel times of multiple shot points and multiple receiver points, respectively, to obtain a superimposed profile of the work area, and determines this superimposed profile as the migration imaging information of the work area. Optionally, the computer equipment performs superimposed imaging on the seismic wave travel times of multiple shot points and multiple receiver points to obtain a superimposed profile of the work area, as shown below. Figure 2 As shown.

[0107] In one possible implementation, for the same work area, the superimposed profile obtained during seismic wave travel is determined using methods in the prior art, such as... Figure 3 As shown. By Figure 2 and Figure 3 The comparison shows that, in shallow strata, the accuracy of the overlay profile obtained by this application is significantly higher than that of the overlay profile obtained by the prior art.

[0108] This application provides a method for determining seismic wave travel time. By using the seismic wave travel times corresponding to multiple simulated shot points and the zero-range travel times corresponding to multiple imaging points, the seismic wave travel times corresponding to multiple shot points and multiple receiver points are determined. The zero-range travel time is the travel time corresponding to the same planar coordinates between the imaging point and the simulated shot point. At this time, the distance between the imaging point and the simulated shot point is less than the distance from the shot point and receiver point to the imaging point. This avoids the seismic wave travel time of the exploration point from exceeding the accuracy range of the travel time interpolation method, thus improving the accuracy of the travel time of the determined exploration point location.

[0109] Figure 4 This is a block diagram illustrating a device for determining the travel time of seismic waves according to an exemplary embodiment. See also... Figure 4 The device includes:

[0110] The acquisition module 401 is used to acquire exploration data within the work area to be determined. The exploration data includes multiple seismic traces. One seismic trace corresponds to one of the multiple shot points and one of the multiple receiver points deployed within the work area.

[0111] The first determining module 402 is used to determine the seismic wave travel time corresponding to multiple simulated shot points on a preset grid within the work area based on multiple seismic trace data. The preset grid includes multiple grid regions, and the seismic wave travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to multiple imaging points within the work area.

[0112] The second determining module 403 is used to determine the zero-shot distance travel time of multiple imaging points in the work area based on the seismic wave travel time corresponding to multiple simulated shot points. The zero-shot distance travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to the imaging point when the plane coordinates of the imaging point and the simulated shot point are the same.

[0113] The third determining module 404 is used to determine the seismic wave travel times corresponding to multiple shot points and multiple receiver points based on multiple seismic trace data, seismic wave travel times corresponding to multiple simulated shot points, and zero-shot distance travel times corresponding to multiple imaging points.

[0114] In one possible implementation, the third determining module 404 is used to determine, for each seismic trace data, a first position of the shot point and a second position of the receiver point corresponding to the seismic trace data; based on the first position, the seismic wave travel time corresponding to multiple simulated shot points and the zero-shot distance travel time corresponding to multiple imaging points, determine the travel time from the shot point to multiple imaging points corresponding to the seismic trace data, thereby obtaining the seismic wave travel time corresponding to the shot point; and based on the second position, the travel time corresponding to multiple simulated shot points and the zero-shot distance travel time corresponding to multiple imaging points, determine the travel time from the receiver point to multiple imaging points corresponding to the seismic trace data, thereby obtaining the seismic wave travel time corresponding to the receiver point.

[0115] In another possible implementation, the third determining module 404 is used to determine, based on the first location, the first target grid region where the shot point corresponding to the seismic trace data is located; determine the seismic wave travel times corresponding to the four first simulated shot points that make up the first target grid region from the seismic wave travel times corresponding to multiple simulated shot points, wherein the seismic wave travel times corresponding to the first simulated shot points are used to represent the travel time of the seismic wave signal from the first simulated shot points to multiple imaging points; for each imaging point, determine the target zero-shot distance travel time corresponding to the imaging point from the zero-shot distance travel times corresponding to multiple imaging points; and for each first simulated shot point, determine the target travel time from the first simulated shot point to the imaging point from the seismic wave travel times corresponding to the first simulated shot point, thereby obtaining the target travel time of the four first simulated shot points; and determine the travel time from the shot point corresponding to the seismic trace data to the imaging point based on the first location, the target zero-shot distance travel time corresponding to the imaging point, and the target travel time of the four first simulated shot points, thereby obtaining the travel time from the shot point to multiple imaging points.

[0116] In another possible implementation, the third determining module 404 is used to determine the location of the imaging point and the locations of the four first simulated shot points. Based on the location of the imaging point and the locations of the four first simulated shot points, four triangular regions are determined, each triangular region including one imaging point and two adjacent first simulated shot points. Based on the first location, the location of the imaging point, and the locations of the four first simulated shot points, the target triangular region where the first location is located is determined from the four triangular regions. Based on the target travel time of the two adjacent first simulated shot points corresponding to the target triangular region and the target zero-shot distance travel time corresponding to the imaging point, the travel time from the shot point to the imaging point corresponding to the seismic trace data is determined.

[0117] In another possible implementation, the second determining module 403 is used to determine the zero-range travel time corresponding to multiple simulated shot points from the seismic wave travel time of multiple simulated shot points; for each imaging point, the second target grid region where the imaging point is located is determined, and the zero-range travel time corresponding to the four second simulated shot points constituting the second target grid region is determined from the zero-range travel time corresponding to the multiple simulated shot points; based on the zero-range travel time corresponding to the four second simulated shot points, the zero-range travel time of the imaging point is determined, thus obtaining the zero-range travel time corresponding to multiple imaging points in the work area.

[0118] In another possible implementation, the first determining module 402 is used to determine the velocity parameters of the seismic wave signal propagating in the reservoir of the work area and the anisotropy parameters of the reservoir based on multiple seismic trace data; and to determine the seismic wave travel time corresponding to multiple simulated shot points on a preset grid in the work area based on the velocity parameters and anisotropy parameters.

[0119] In another possible implementation, the device further includes:

[0120] The fourth determination module is used to determine the migration imaging information of the work area based on the seismic wave travel times corresponding to multiple shot points and multiple receiver points.

[0121] This application provides a device for determining seismic wave travel time. By using the seismic wave travel time corresponding to multiple simulated shot points and the zero-range travel time corresponding to multiple imaging points, the seismic wave travel time corresponding to multiple shot points and multiple receiver points is determined. The zero-range travel time is the travel time corresponding to the same planar coordinates between the imaging point and the simulated shot point. At this time, the distance between the imaging point and the simulated shot point is less than the distance from the shot point and receiver point to the imaging point, thus avoiding the seismic wave travel time of the exploration point from exceeding the accuracy range of the travel time interpolation method. Therefore, the accuracy of the travel time of the determined exploration point location is improved.

[0122] Figure 5This diagram illustrates a structural block diagram of a computer device 500 provided in an exemplary embodiment of the present invention. The computer device 500 may be a smartphone, tablet computer, MP3 player (Moving Picture Experts Group Audio Layer III), MP4 player (Moving Picture Experts Group Audio Layer IV), laptop computer, or desktop computer. The computer device 500 may also be referred to as a user device, portable computer device, laptop computer device, desktop computer device, or other names.

[0123] Typically, computer device 500 includes a processor 501 and a memory 502.

[0124] Processor 501 may include one or more processing cores, such as a quad-core processor, an octa-core processor, etc. Processor 501 may be implemented using at least one hardware form selected from DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). Processor 501 may also include a main processor and a coprocessor. The main processor, also known as a CPU (Central Processing Unit), is used to process data in the wake-up state; the coprocessor is a low-power processor used to process data in the standby state. In some embodiments, processor 501 may integrate a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the screen. In some embodiments, processor 501 may also include an AI (Artificial Intelligence) processor, which is used to handle computational operations related to machine learning.

[0125] Memory 502 may include one or more computer-readable storage media, which may be non-transitory. Memory 502 may also include high-speed random access memory and non-volatile memory, such as one or more disk storage devices or flash memory devices. In some embodiments, the non-transitory computer-readable storage media in memory 502 is used to store at least one instruction, which is executed by processor 501 to implement the method for determining seismic wave travel time provided in the method embodiments of this application.

[0126] In some embodiments, the computer device 500 may also optionally include a peripheral device interface 503 and at least one peripheral device. The processor 501, memory 502, and peripheral device interface 503 can be connected via a bus or signal line. Each peripheral device can be connected to the peripheral device interface 503 via a bus, signal line, or circuit board. Specifically, the peripheral device includes at least one of the following: a radio frequency circuit 504, a display screen 505, a camera 506, an audio circuit 507, a positioning component 508, and a power supply 509.

[0127] Peripheral device interface 503 can be used to connect at least one I / O (Input / Output) related peripheral device to processor 501 and memory 502. In some embodiments, processor 501, memory 502 and peripheral device interface 503 are integrated on the same chip or circuit board; in some other embodiments, any one or two of processor 501, memory 502 and peripheral device interface 503 can be implemented on separate chips or circuit boards, which is not limited in this embodiment.

[0128] The radio frequency (RF) circuit 504 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The RF circuit 504 communicates with communication networks and other communication devices via electromagnetic signals. The RF circuit 504 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals back into electrical signals. Optionally, the RF circuit 504 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a user identity module card, etc. The RF circuit 504 can communicate with other computer devices via at least one wireless communication protocol. This wireless communication protocol includes, but is not limited to: metropolitan area networks (MANs), various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks (WLANs), and / or WiFi (Wireless Fidelity) networks. In some embodiments, the RF circuit 504 may also include circuitry related to NFC (Near Field Communication), which is not limited in this application.

[0129] Display screen 505 is used to display a UI (User Interface). This UI may include graphics, text, icons, videos, and any combination thereof. When display screen 505 is a touch display screen, it also has the ability to collect touch signals on or above its surface. These touch signals can be input as control signals to processor 501 for processing. In this case, display screen 505 can also be used to provide virtual buttons and / or a virtual keyboard, also known as soft buttons and / or a soft keyboard. In some embodiments, there may be one display screen 505, which serves as the front panel of the computer device 500; in other embodiments, there may be at least two display screens, respectively disposed on different surfaces of the computer device 500 or in a folded design; in still other embodiments, display screen 505 may be a flexible display screen, disposed on a curved or folded surface of the computer device 500. Furthermore, display screen 505 may be configured as a non-rectangular, irregular shape, i.e., a non-rectangular screen. Display screen 505 may be made of materials such as LCD (Liquid Crystal Display) or OLED (Organic Light-Emitting Diode).

[0130] The camera assembly 506 is used to acquire images or videos. Optionally, the camera assembly 506 includes a front-facing camera and a rear-facing camera. Typically, the front-facing camera is located on the front panel of the computer device, and the rear-facing camera is located on the back of the computer device. In some embodiments, there are at least two rear-facing cameras, which are any one of a main camera, a depth-sensing camera, a wide-angle camera, and a telephoto camera, to achieve background blurring by fusion of the main camera and the depth-sensing camera, panoramic shooting by fusion of the main camera and the wide-angle camera, VR (Virtual Reality) shooting, or other fusion shooting functions. In some embodiments, the camera assembly 506 may also include a flash. The flash can be a single-color temperature flash or a dual-color temperature flash. A dual-color temperature flash is a combination of a warm-light flash and a cool-light flash, which can be used for light compensation at different color temperatures.

[0131] The audio circuit 507 may include a microphone and a speaker. The microphone is used to collect sound waves from the user and the environment, converting the sound waves into electrical signals that are input to the processor 501 for processing, or input to the radio frequency circuit 504 for voice communication. For stereo sound acquisition or noise reduction purposes, multiple microphones may be used, each located in a different part of the computer device 500. The microphone may also be an array microphone or an omnidirectional microphone. The speaker is used to convert the electrical signals from the processor 501 or the radio frequency circuit 504 into sound waves. The speaker may be a conventional diaphragm speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can convert electrical signals not only into audible sound waves but also into inaudible sound waves for purposes such as distance measurement. In some embodiments, the audio circuit 507 may also include a headphone jack.

[0132] The positioning component 508 is used to locate the current geographical location of the computer device 500 in order to enable navigation or LBS (Location Based Service). The positioning component 508 can be a positioning component based on the US GPS (Global Positioning System), China's BeiDou system, Russia's Granas system, or the European Union's Galileo system.

[0133] Power supply 509 is used to supply power to the various components in computer device 500. Power supply 509 can be AC ​​power, DC power, a disposable battery, or a rechargeable battery. When power supply 509 includes a rechargeable battery, the rechargeable battery can support wired or wireless charging. The rechargeable battery can also be used to support fast charging technology.

[0134] In some embodiments, the computer device 500 further includes one or more sensors 510. The one or more sensors 510 include, but are not limited to: an accelerometer 511, a gyroscope 512, a pressure sensor 513, a fingerprint sensor 514, an optical sensor 515, and a proximity sensor 516.

[0135] Accelerometer 511 can detect the magnitude of acceleration along the three coordinate axes of a coordinate system established by computer device 500. For example, accelerometer 511 can be used to detect the components of gravitational acceleration along the three coordinate axes. Processor 501 can control display screen 505 to display the user interface in either a landscape or portrait view based on the gravitational acceleration signal acquired by accelerometer 511. Accelerometer 511 can also be used for games or for acquiring user motion data.

[0136] The gyroscope sensor 512 can detect the orientation and rotation angle of the computer device 500. The gyroscope sensor 512, in conjunction with the accelerometer sensor 511, can collect 3D motion data from the user on the computer device 500. Based on the data collected by the gyroscope sensor 512, the processor 501 can perform the following functions: motion sensing (e.g., changing the UI based on the user's tilt), image stabilization during shooting, game control, and inertial navigation.

[0137] The pressure sensor 513 can be disposed on the side bezel of the computer device 500 and / or on the lower layer of the display screen 505. When the pressure sensor 513 is disposed on the side bezel of the computer device 500, it can detect the user's grip signal on the computer device 500, and the processor 501 can perform left / right hand recognition or quick operation based on the grip signal collected by the pressure sensor 513. When the pressure sensor 513 is disposed on the lower layer of the display screen 505, the processor 501 can control the operable controls on the UI interface based on the user's pressure operation on the display screen 505. The operable controls include at least one of button controls, scroll bar controls, icon controls, and menu controls.

[0138] The fingerprint sensor 514 is used to collect a user's fingerprint. The processor 501 identifies the user based on the fingerprint collected by the fingerprint sensor 514, or vice versa. When the user's identity is verified as trusted, the processor 501 authorizes the user to perform relevant sensitive operations, including unlocking the screen, viewing encrypted information, downloading software, making payments, and changing settings. The fingerprint sensor 514 can be located on the front, back, or side of the computer device 500. When the computer device 500 has physical buttons or a manufacturer's logo, the fingerprint sensor 514 can be integrated with the physical buttons or the manufacturer's logo.

[0139] An optical sensor 515 is used to collect ambient light intensity. In one embodiment, the processor 501 can control the display brightness of the display screen 505 based on the ambient light intensity collected by the optical sensor 515. Specifically, when the ambient light intensity is high, the display brightness of the display screen 505 is increased; when the ambient light intensity is low, the display brightness of the display screen 505 is decreased. In another embodiment, the processor 501 can also dynamically adjust the shooting parameters of the camera assembly 506 based on the ambient light intensity collected by the optical sensor 515.

[0140] The proximity sensor 516, also known as a distance sensor, is typically located on the front panel of the computer device 500. The proximity sensor 516 is used to detect the distance between the user and the front of the computer device 500. In one embodiment, when the proximity sensor 516 detects that the distance between the user and the front of the computer device 500 is gradually decreasing, the processor 501 controls the display screen 505 to switch from a screen-on state to a screen-off state; when the proximity sensor 516 detects that the distance between the user and the front of the computer device 500 is gradually increasing, the processor 501 controls the display screen 505 to switch from a screen-off state to a screen-on state.

[0141] Those skilled in the art will understand that Figure 5 The structure shown does not constitute a limitation on the computer device 500, and may include more or fewer components than shown, or combine certain components, or use different component arrangements.

[0142] This application also provides a computer-readable storage medium storing at least one piece of program code, which is loaded and executed by a processor to implement the operations performed in the method for determining the travel time of seismic waves in this application.

[0143] This application also provides a computer program product, which includes at least one piece of program code. The at least one piece of program code is loaded and executed by a processor to implement the operation performed by the method for determining the travel time of seismic waves in this application.

[0144] In some embodiments, the computer program involved in the present application embodiments may be deployed and executed on a computer device, or executed on multiple computer devices located in one location, or executed on multiple computer devices distributed in multiple locations and interconnected through a communication network. Multiple computer devices distributed in multiple locations and interconnected through a communication network may constitute a blockchain system.

[0145] The above description is only for the purpose of enabling those skilled in the art to understand the technical solution of this application, and is not intended to limit this application. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.

Claims

1. A method for determining the travel time of seismic waves, characterized in that, The method includes: Obtain exploration data within the work area to be determined. The exploration data includes multiple seismic traces, with one seismic trace corresponding to one of the multiple shot points and multiple receiver points deployed within the work area. Based on the multiple seismic trace data, the seismic wave travel time corresponding to multiple simulated shot points on a preset grid within the work area is determined. The preset grid includes multiple grid regions, and the seismic wave travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to multiple imaging points within the work area. Based on the seismic wave travel time corresponding to the multiple simulated shot points, the zero shot distance travel time corresponding to the multiple imaging points in the work area is determined. The zero shot distance travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to the imaging point when the plane coordinates of the imaging point and the simulated shot point are the same. For each seismic trace, determine the first position of the shot point and the second position of the receiver point corresponding to the seismic trace; Based on the first location, the seismic wave travel time corresponding to the plurality of simulated shot points, and the zero-shot distance travel time corresponding to the plurality of imaging points, the travel time from the shot point corresponding to the seismic trace data to the plurality of imaging points is determined, thereby obtaining the seismic wave travel time corresponding to the shot point. Furthermore, based on the second location, the travel time corresponding to the plurality of simulated shot points, and the zero-shot distance travel time corresponding to the plurality of imaging points, the travel time from the receiver point corresponding to the seismic trace data to the plurality of imaging points is determined, thereby obtaining the seismic wave travel time corresponding to the receiver point.

2. The method according to claim 1, characterized in that, The determination of the travel time from the shot point corresponding to the seismic trace data to the multiple imaging points based on the first location, the seismic wave travel time corresponding to the multiple simulated shot points, and the zero-shot distance travel time corresponding to the multiple imaging points includes: Based on the first location, the first target grid region where the shot point corresponding to the seismic trace data is located is determined. The seismic wave travel time corresponding to the four first simulated shot points that make up the first target grid region is determined from the seismic wave travel time corresponding to the multiple simulated shot points. The seismic wave travel time corresponding to the first simulated shot point is used to represent the travel time of the seismic wave signal from the first simulated shot point to the multiple imaging points. For each imaging point, the target zero-range travel time corresponding to the imaging point is determined from the zero-range travel time corresponding to the plurality of imaging points, and for each first simulated shot point, the target travel time from the first simulated shot point to the imaging point is determined from the seismic wave travel time corresponding to the first simulated shot point, thus obtaining the target travel time of the four first simulated shot points. Based on the first location, the target zero-range travel time corresponding to the imaging point, and the target travel time of the four first simulated shot points, the travel time from the shot point corresponding to the seismic trace data to the imaging point is determined, and the travel time from the shot point to the plurality of imaging points is obtained.

3. The method according to claim 2, characterized in that, The determination of the travel time from the shot point corresponding to the imaging point based on the first location, the target zero-range travel time corresponding to the imaging point, and the target travel time of the four first simulated shot points includes: The positions of the imaging point and the four first simulated shot points are determined. Based on the positions of the imaging point and the four first simulated shot points, four triangular regions are determined. Each triangular region includes one imaging point and two adjacent first simulated shot points. Based on the first position, the position of the imaging point, and the positions of the four first simulated gun points, the target triangular region where the first position is located is determined from the four triangular regions; Based on the target travel time of two adjacent first simulated shot points corresponding to the target triangular region and the target zero-shot distance travel time corresponding to the imaging point, the travel time from the shot point corresponding to the seismic trace data to the imaging point is determined.

4. The method according to claim 1, characterized in that, The step of determining the zero-range travel time corresponding to multiple imaging points within the work area based on the seismic wave travel time corresponding to the multiple simulated shot points includes: Determine the zero-distance travel time corresponding to the multiple simulated shot points from the seismic wave travel times of the multiple simulated shot points; For each imaging point, determine the second target grid region where the imaging point is located, and determine the zero-range travel time corresponding to the four second simulated gun points that make up the second target grid region from the zero-range travel time corresponding to the plurality of simulated gun points; Based on the zero-range travel time corresponding to the four second simulated shot points, the zero-range travel time of the imaging point is determined, and the zero-range travel time corresponding to multiple imaging points in the work area is obtained.

5. The method according to claim 1, characterized in that, The step of determining the seismic wave travel time corresponding to multiple simulated shot points on a preset grid within the work area based on the multiple seismic trace data includes: Based on the multiple seismic traces, the velocity parameters of the seismic wave signal propagation within the reservoir in the work area and the anisotropy parameters of the reservoir are determined. Based on the velocity parameters and the anisotropy parameters, the seismic wave travel times corresponding to multiple simulated shot points on a preset grid within the work area are determined.

6. The method according to claim 1, characterized in that, The method further includes: Based on the seismic wave travel times corresponding to the multiple shot points and the multiple receiver points, the migration imaging information of the work area is determined.

7. A device for determining the travel time of seismic waves, characterized in that, The device includes: The acquisition module is used to acquire exploration data within the work area to be determined. The exploration data includes multiple seismic traces, and one seismic trace corresponds to one of the multiple shot points and one of the multiple receiver points deployed within the work area. The first determining module is used to determine the seismic wave travel time corresponding to multiple simulated shot points on a preset grid within the work area based on the multiple seismic trace data. The preset grid includes multiple grid regions, and the seismic wave travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to multiple imaging points within the work area. The second determining module is used to determine the zero-shot distance travel time corresponding to multiple imaging points in the work area based on the seismic wave travel time corresponding to the multiple simulated shot points. The zero-shot distance travel time is used to represent the propagation time of the seismic wave signal from the simulated shot point to the imaging point when the plane coordinates of the imaging point and the simulated shot point are the same. The third determining module is used to determine the first position of the shot point and the second position of the receiver point corresponding to each seismic trace data. Based on the first location, the seismic wave travel time corresponding to the plurality of simulated shot points, and the zero-shot distance travel time corresponding to the plurality of imaging points, the travel time from the shot point corresponding to the seismic trace data to the plurality of imaging points is determined, thereby obtaining the seismic wave travel time corresponding to the shot point. Furthermore, based on the second location, the travel time corresponding to the plurality of simulated shot points, and the zero-shot distance travel time corresponding to the plurality of imaging points, the travel time from the receiver point corresponding to the seismic trace data to the plurality of imaging points is determined, thereby obtaining the seismic wave travel time corresponding to the receiver point.

8. A computer device, characterized in that, The computer device includes: A processor and a memory, wherein the memory stores at least one line of program code, which is loaded and executed by the processor to perform the operations performed in the method for determining the travel time of seismic waves according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one piece of program code, which is loaded and executed by a processor to perform the operations performed in the method for determining the travel time of seismic waves as described in any one of claims 1 to 6.

Images