CN116136605A - Method and device for determining travel time of seismic waves, computer equipment and storage medium - Google Patents

Method and device for determining travel time of seismic waves, computer equipment and storage medium Download PDF

Info

Publication number
CN116136605A
CN116136605A CN202111356613.5A CN202111356613A CN116136605A CN 116136605 A CN116136605 A CN 116136605A CN 202111356613 A CN202111356613 A CN 202111356613A CN 116136605 A CN116136605 A CN 116136605A
Authority
CN
China
Prior art keywords
travel time
seismic
determining
shots
shot
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111356613.5A
Other languages
Chinese (zh)
Inventor
张巍毅
王成祥
吴建鲁
岳玉波
李文静
张建磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN202111356613.5A priority Critical patent/CN116136605A/en
Publication of CN116136605A publication Critical patent/CN116136605A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures

Abstract

The application provides a method and a device for determining travel time of seismic waves, computer equipment and a storage medium, and belongs to the technical field of seismic exploration. The method comprises the following steps: acquiring exploration data in a work area to be determined, wherein the exploration data comprises a plurality of seismic channel data; determining seismic wave travel time corresponding to a plurality of simulated shots on a preset grid in a work area according to the plurality of seismic channel data; determining zero-offset travel time corresponding to a plurality of imaging points in a work area according to the seismic wave travel time corresponding to a plurality of simulated shots; and determining the seismic wave travel time corresponding to the plurality of shot points and the plurality of detection points according to the plurality of seismic channel data, the seismic wave travel time corresponding to the plurality of simulated shot points and the zero-shot travel time corresponding to the plurality of imaging points. As the zero-offset seismic wave travel time is the travel time corresponding to the time when the imaging point is the same as the plane coordinate of the simulated shot point, the seismic wave travel time of the exploration point is prevented from exceeding the precision range, and the accuracy of the determined seismic wave travel time is improved.

Description

Method and device for determining travel time of seismic waves, computer equipment and storage medium
Technical Field
The present disclosure relates to the field of seismic exploration, and in particular, to a method and apparatus for determining travel time of seismic waves, a computer device, and a storage medium.
Background
Currently, seismic migration imaging techniques are an important technique for determining the geologic structure of a reservoir. The seismic wave migration imaging technology is to determine the superposition profile of a reservoir according to the seismic wave travel time of exploration point positions such as a shot point, a wave detection point and the like, and characterize the geological structure of the reservoir through the superposition profile, so that the seismic wave travel time of the exploration point positions needs to be determined.
In the related art, determining the seismic wave travel time of the exploration point position by a travel time interpolation method; the method comprises the following steps: determining travel time of a plurality of simulated shots on a preset grid; determining a grid area where the exploration point is located according to the position of the exploration point; determining travel time of the exploration point position by linear interpolation according to the positions and travel time of four simulated shots composing the grid area.
However, in the above related art, when the distances between the exploration point and the imaging point are smaller than the distances from the four simulated shots to the imaging point, the accuracy of the interpolation method at the time of travel is exceeded, resulting in low accuracy at the time of travel of the determined exploration point position.
Disclosure of Invention
The embodiment of the application provides a method, a device, computer equipment and a storage medium for determining seismic wave travel, which can improve the accuracy of drilling coincidence information of a detection zone. The technical scheme is as follows:
in one aspect, the present application provides a method of determining when a seismic wave travels, the method comprising:
acquiring exploration data in a work area to be determined, wherein the exploration data comprises a plurality of seismic channel data, and one seismic channel data corresponds to one shot point and one wave detection point of a plurality of shot points and a plurality of wave detection points which are distributed in the work area;
determining seismic wave travel time corresponding to a plurality of simulated shots on a preset grid in the work area according to the plurality of seismic channel data, wherein the preset grid comprises a plurality of grid areas, and the seismic wave travel time is used for representing the propagation time of seismic wave signals from the simulated shots to a plurality of imaging points in the work area;
determining zero-offset travel corresponding to a plurality of imaging points in the work area according to the seismic wave travel corresponding to the plurality of simulated shots, wherein the zero-offset travel is used for representing the propagation time of a seismic wave signal from the simulated shots to the imaging points when the plane coordinates of the imaging points and the simulated shots are the same;
And determining the seismic wave travel time corresponding to the plurality of shot points and the plurality of detection points according to the plurality of seismic channel data, the seismic wave travel time corresponding to the plurality of simulated shot points and the zero-offset travel time corresponding to the plurality of imaging points.
In a possible implementation manner, the determining the seismic wave travel time corresponding to the plurality of shots and the plurality of detectors according to the seismic trace data, the seismic wave travel time corresponding to the plurality of simulated shots and the zero-offset travel time corresponding to the plurality of imaging points includes:
for each seismic trace data, determining a first position of a shot point and a second position of a wave detection point corresponding to the seismic trace data;
determining travel time of the shot corresponding to the seismic channel data to travel time of the imaging points based on the first position, the travel time of the seismic waves corresponding to the simulated shots and the travel time of zero distances corresponding to the imaging points, obtaining the travel time of the seismic waves corresponding to the shots, and determining travel time of the detector corresponding to the seismic channel data to travel time of the imaging points based on the second position, the travel time of the simulated shots and the travel time of zero distances corresponding to the imaging points.
In another possible implementation manner, the determining the travel time of the shot corresponding to the seismic trace data to the imaging points based on the first position, the travel time of the seismic waves corresponding to the simulated shots and the travel time of zero distances corresponding to the imaging points includes:
determining a first target grid area where the shot corresponding to the seismic channel data is located based on the first position, and determining the travel time of the seismic waves corresponding to four first simulation shots forming the first target grid area from the travel time of the seismic waves corresponding to the simulation shots, wherein the travel time of the seismic waves corresponding to the first simulation shots is used for representing the travel time of the seismic wave signals from the first simulation shots to the imaging points;
for each imaging point, determining a target zero-offset travel time corresponding to the imaging point from zero-offset travel times corresponding to the imaging points, and for each first simulation shot, determining a target travel time from the first simulation shot to the imaging point from a seismic wave travel time corresponding to the first simulation shot, thereby obtaining target travel times of the four first simulation shots;
And determining travel time from the shot corresponding to the seismic channel data to the imaging point based on the first position, the target zero-shot travel time corresponding to the imaging point and the target travel time of the four first simulation shots, and obtaining travel time from the shot to the imaging points.
In another possible implementation manner, the determining the travel of the shot corresponding to the seismic trace data to the imaging point based on the first position, the target zero-shot travel time corresponding to the imaging point, and the target travel time of the four first simulated shots includes:
determining the positions of the imaging point and the positions of the four first simulation shots, and determining four triangular areas based on the positions of the imaging point and the positions of the four first simulation shots, wherein one triangular area comprises one imaging point and two adjacent first simulation shots;
determining a target triangle area where the first position is located from the four triangle areas based on the first position, the position of the imaging point and the positions of the four first simulation shots;
and determining travel time from the shot corresponding to the seismic channel data to the imaging point based on the travel time of the target of the two adjacent first simulation shots corresponding to the target triangle area and the travel time of the target zero-offset corresponding to the imaging point.
In another possible implementation manner, the determining, according to the seismic wave travel time corresponding to the plurality of simulated shots, zero-offset travel time corresponding to a plurality of imaging points in the work area includes:
determining zero-offset travel time corresponding to the plurality of simulated shots from the seismic wave travel time of the plurality of simulated shots;
for each imaging point, determining a second target grid area where the imaging point is located, and determining zero-offset travel times corresponding to four second simulation shots composing the second target grid area from zero-offset travel times corresponding to the simulation shots;
and determining zero-offset travel time of the imaging point based on the zero-offset travel time corresponding to the four second simulation shot points, and obtaining the zero-offset travel time corresponding to a plurality of imaging points in the work area.
In another possible implementation manner, the determining, according to the plurality of seismic trace data, when the seismic waves corresponding to the plurality of simulated shots on the preset grid in the work area travel, includes:
determining a speed parameter of seismic wave signal propagation in a reservoir of the work area and an anisotropy parameter of the reservoir according to the plurality of seismic channel data;
And determining seismic wave travel time corresponding to a plurality of simulated shots on a preset grid in the work area according to the speed parameter and the anisotropy parameter.
In another possible implementation, the method further includes:
and determining offset imaging information of the work area based on the plurality of shot points and the seismic wave travel time corresponding to the plurality of detection points.
In another aspect, the present application provides a device for determining when a seismic wave travels, the device comprising:
the acquisition module is used for acquiring exploration data in a work area to be determined, wherein the exploration data comprises a plurality of seismic channel data, and one seismic channel data corresponds to one shot point and one detection point of a plurality of shot points and a plurality of detection points which are distributed in the work area;
the first determining module is used for determining seismic wave travel time corresponding to a plurality of simulated shots on a preset grid in the work area according to the plurality of seismic channel data, wherein the preset grid comprises a plurality of grid areas, and the seismic wave travel time is used for representing the propagation time of a seismic wave signal from the simulated shots to a plurality of imaging points in the work area;
the second determining module is used for determining zero-offset travel corresponding to a plurality of imaging points in the work area according to the seismic wave travel time corresponding to the plurality of simulated shots, and the zero-offset travel time is used for representing the propagation time length of a seismic wave signal from the simulated shots to the imaging points when the plane coordinates of the imaging points are the same as those of the simulated shots;
And the third determining module is used for determining the seismic wave travel time corresponding to the plurality of shot points and the detection points according to the plurality of seismic channel data, the seismic wave travel time corresponding to the plurality of simulated shot points and the zero-offset travel time corresponding to the plurality of imaging points.
In a possible implementation manner, the third determining module is configured to determine, for each seismic trace data, a first position of a shot point and a second position of a detector point corresponding to the seismic trace data; determining travel time of the shot corresponding to the seismic channel data to travel time of the imaging points based on the first position, the travel time of the seismic waves corresponding to the simulated shots and the travel time of zero distances corresponding to the imaging points, obtaining the travel time of the seismic waves corresponding to the shots, and determining travel time of the detector corresponding to the seismic channel data to travel time of the imaging points based on the second position, the travel time of the simulated shots and the travel time of zero distances corresponding to the imaging points.
In another possible implementation manner, the third determining module is configured to determine, based on the first location, a first target grid area where a shot corresponding to the seismic trace data is located, and determine, from among travel times of seismic waves corresponding to the plurality of simulated shots, travel times of seismic waves corresponding to four first simulated shots that constitute the first target grid area, where travel times of seismic waves corresponding to the first simulated shots are used to represent travel times of seismic wave signals from the first simulated shot to the plurality of imaging points; for each imaging point, determining a target zero-offset travel time corresponding to the imaging point from zero-offset travel times corresponding to the imaging points, and for each first simulation shot, determining a target travel time from the first simulation shot to the imaging point from a seismic wave travel time corresponding to the first simulation shot, thereby obtaining target travel times of the four first simulation shots; and determining travel time from the shot corresponding to the seismic channel data to the imaging point based on the first position, the target zero-shot travel time corresponding to the imaging point and the target travel time of the four first simulation shots, and obtaining travel time from the shot to the imaging points.
In another possible implementation manner, the third determining module is configured to determine a position of the imaging point and positions of the four first simulated shots, and determine four triangular areas based on the positions of the imaging point and the positions of the four first simulated shots, where one triangular area includes one imaging point and two adjacent first simulated shots; determining a target triangle area where the first position is located from the four triangle areas based on the first position, the position of the imaging point and the positions of the four first simulation shots; and determining travel time from the shot corresponding to the seismic channel data to the imaging point based on the travel time of the target of the two adjacent first simulation shots corresponding to the target triangle area and the travel time of the target zero-offset corresponding to the imaging point.
In another possible implementation manner, the second determining module is configured to determine zero-offset travel times corresponding to the multiple simulated shots from the seismic wave travel times of the multiple simulated shots; for each imaging point, determining a second target grid area where the imaging point is located, and determining zero-offset travel times corresponding to four second simulation shots composing the second target grid area from zero-offset travel times corresponding to the simulation shots; and determining zero-offset travel time of the imaging point based on the zero-offset travel time corresponding to the four second simulation shot points, and obtaining the zero-offset travel time corresponding to a plurality of imaging points in the work area.
In another possible implementation, the first determining module is configured to determine, from the plurality of seismic trace data, a velocity parameter of propagation of a seismic wave signal within a reservoir of the work area and an anisotropy parameter of the reservoir; and determining seismic wave travel time corresponding to a plurality of simulated shots on a preset grid in the work area according to the speed parameter and the anisotropy parameter.
In another possible implementation, the apparatus further includes:
and the fourth determining module is used for determining offset imaging information of the work area based on the plurality of shot points and the seismic waves corresponding to the plurality of wave detection points when traveling.
In another aspect, embodiments of the present application provide a computer device, the computer device comprising: a processor and a memory having stored therein at least one program code loaded and executed by the processor to perform the operations performed in the method of determining travel of a seismic wave as described in any of the possible implementations above.
In another aspect, embodiments of the present application provide a computer-readable storage medium having at least one program code stored therein, the at least one program code loaded and executed by a processor to implement operations performed in a method of determining when a seismic wave travels as described in any one of the possible implementations above.
In another aspect, embodiments of the present application provide a computer program product comprising at least one program code loaded and executed by a processor to implement operations performed in a method of determining when a seismic wave is traveling as described in any one of the possible implementations.
The beneficial effects of the technical scheme provided by the embodiment of the application at least comprise:
the embodiment of the application provides a method for determining when a seismic wave travels, because when the seismic wave travels corresponding to a plurality of simulated shots and the seismic wave travels corresponding to a plurality of imaging points are zero-offset, the seismic wave travels corresponding to a plurality of shots and a plurality of detection points are determined, and when the zero-offset travels are the travel corresponding to the same plane coordinates of the imaging points and the simulated shots, the distance between the imaging points and the simulated shots is smaller than the distance from the exploration points such as the shots and the detection points to the imaging points, the seismic wave travel of the exploration points is avoided from exceeding the accuracy range of interpolation methods when traveling, and therefore the accuracy of the determined travel of the exploration point position is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart illustrating a method of determining when a seismic wave is traveling, according to an exemplary embodiment;
FIG. 2 is a seismic profile view illustrating a determination of travel time of a seismic wave determined by the present application, according to an exemplary embodiment;
FIG. 3 is a seismic profile diagram illustrating a seismic travel time determination determined by the prior art according to an exemplary embodiment;
FIG. 4 is a block diagram illustrating a seismic travel determination apparatus according to an exemplary embodiment;
fig. 5 is a block diagram of a computer device, according to an example embodiment.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present application more apparent, the embodiments of the present application will be described in further detail below with reference to the accompanying drawings.
FIG. 1 is a flowchart of a method of determining when a seismic wave travels, performed by a computer device, according to an exemplary embodiment. Referring to fig. 1, the method includes:
101. the computer device obtains exploration data in a work area to be determined, wherein the exploration data comprises a plurality of seismic trace data, and one seismic trace data corresponds to one shot point and one wave detection point of a plurality of shot points and a plurality of wave detection points which are distributed in the work area.
In one possible implementation, a plurality of shots and a plurality of detectors are arranged in the work area to be determined; the shot point is used for transmitting the seismic wave signals, and the wave detection point is used for receiving the seismic wave signals. Correspondingly, the method comprises the following steps: receiving seismic wave signals emitted by a plurality of shot points by a plurality of wave detectors to obtain a plurality of seismic channel data, wherein one seismic channel data comprises the seismic wave signals emitted by one shot point received by one wave detector; uploading the plurality of seismic channel data to computer equipment, and obtaining the plurality of seismic channel data of the work area by the computer equipment to obtain exploration data corresponding to a plurality of shots and a plurality of detection points. The seismic channel data comprise seismic wave signal change data within a preset time length.
In the embodiment of the present application, the value of the preset duration is not specifically limited, and may be set as required. Optionally, the preset duration is any value between 1s and 20s, for example, the preset duration is 6s, 8s, 10s, etc.
102. The computer equipment determines the travel time of the earthquake waves corresponding to a plurality of simulated shots on a preset grid in the work area according to the plurality of earthquake channel data, wherein the preset grid comprises a plurality of grid areas, one grid area is a quadrilateral area formed by four adjacent simulated shots, and the travel time of the earthquake waves is used for representing the propagation time of the earthquake wave signals from the simulated shots to a plurality of imaging points in the work area.
In this step, the distance between any two adjacent simulated shots may be the same or different. In one possible implementation, the grid area is rectangular, the quadrilateral area formed by four adjacent simulated shots is rectangular, the distance between any two adjacent simulated shots is the same, and the distance is a first preset distance. In the embodiment of the present application, the value of the first preset distance is not specifically limited, and may be set and modified as required. Optionally, the first preset distance is any value between 10m and 300m, for example, the first preset distance is 10m, 50m, 100m, etc.
In one possible implementation, the computer device determines seismic travel times for a plurality of simulated shots based on survey data received from actual shots. Correspondingly, the computer equipment determines the earthquake wave travel time corresponding to a plurality of simulated shots on a preset grid in the work area according to the plurality of earthquake channel data, and comprises the following steps: the computer equipment determines the speed parameter of the propagation of the seismic wave signals in the reservoir of the work area and the anisotropy parameter of the reservoir according to the plurality of seismic channel data; and determining the travel time of the seismic waves corresponding to the plurality of simulated shots on the preset grid in the work area according to the speed parameter and the anisotropy parameter.
It should be noted that the imaging point may be any point in the reservoir. In one possible implementation, the distance between any two adjacent imaging points is the same, and the distance is a second preset distance. In the embodiment of the present application, the value of the second preset distance is not specifically limited, and may be set and modified as required. Optionally, the second preset distance is any value between 10m and 100m, for example, the second preset distance is 10m, 50m, 80m, etc.
The density of the grid formed by the plurality of imaging points can be greater than that of the preset grid or less than that of the preset grid. Correspondingly, when the distance between any two adjacent imaging points is smaller than the distance between any two adjacent simulated shots, the density of the grid formed by the imaging points is greater than that of a preset grid; when the distance between any two adjacent imaging points is larger than the distance between any two adjacent simulated shots, the density of the grid formed by the plurality of imaging points is smaller than that of the preset grid.
103. The computer equipment determines the propagation time length of the seismic wave signals from the simulated shot point to the imaging point when the zero-offset travel corresponding to the imaging points in the work area is performed according to the seismic wave travel time of the simulated shot points and the plane coordinates of the simulated shot points are the same when the zero-offset travel time is used for representing the imaging points to be the same as the plane coordinates of the simulated shot points.
In one possible implementation, the computer device determines zero-offset travel times for the plurality of imaging points based on the zero-offset travel times for the plurality of simulated shots. Accordingly, the computer apparatus may include the following steps (11) to (13) when determining zero offset travel corresponding to a plurality of imaging points within the work area:
(11) The computer device determines zero-offset travel times corresponding to the plurality of simulated shots from among the seismic travel times of the plurality of simulated shots.
In one possible implementation, this step is: for each simulated shot, determining a target imaging point which is the same as the position of the simulated shot from a plurality of imaging points according to the position of the simulated shot and the positions of the plurality of imaging points by the computer equipment; and selecting the travel time from the simulated shot point to the target imaging point from the seismic wave travel time of the simulated shot point, and obtaining the zero-shot travel time corresponding to the simulated shot point.
(12) For each imaging point, the computer equipment determines a second target grid area where the imaging point is located, and determines zero-offset travel time corresponding to four second simulation shots which form the second target grid area from zero-offset travel time corresponding to a plurality of simulation shots.
In one possible implementation, the computer device stores a coordinate range of a plurality of grid areas included in the preset grid. Correspondingly, for each imaging point, the computer device determines the second target grid area where the imaging point is located by: the computer equipment acquires the coordinates of each imaging point for the imaging point; determining from the coordinates of the imaging point The coordinate range in which the coordinate is located; and determining a second target grid region corresponding to the coordinate range from the coordinate ranges of the grid regions. Alternatively, the coordinates of the imaging point are plane coordinates, for example, the coordinates of the imaging point may be represented by (x) 3 ，y 3 ) Representation, where x 3 Plane abscissa, y representing imaging point 3 Representing the planar ordinate of the imaging point.
In one possible implementation, one grid area is a quadrilateral area formed by four adjacent simulated shots, and the correspondence between the grid area and the four adjacent simulated shots is stored in the computer device. Correspondingly, the computer equipment determines the zero-offset travel time corresponding to the four second simulation shots composing the second target grid area by the steps of: the computer device determines four second simulated shots corresponding to the second target grid region, and determines zero-offset travel times corresponding to the four second simulated shots from the zero-offset travel times corresponding to the plurality of simulated shots.
It should be noted that, when the zero-offset travel time is the travel time from the second simulated shot point to the imaging point when the plane coordinates of the second simulated shot point and the imaging point are the same. Wherein the imaging points with the same plane coordinates as the second simulated shot point comprise a plurality of imaging points with different depths and the same plane coordinates; correspondingly, the zero-offset travel time comprises travel time from the second simulated shot to a plurality of imaging points with different depths and same plane coordinates.
(13) The computer equipment determines zero-offset travel time of the imaging point based on zero-offset travel time corresponding to four second simulation shots, and obtains zero-offset travel time corresponding to a plurality of imaging points in the work area.
In one possible implementation, the computer device determines a zero-offset travel time for the imaging point by the bilinear difference. Correspondingly, the method comprises the following steps: the computer equipment determines zero-offset travel time of the imaging point along the abscissa direction according to zero-offset travel time corresponding to the four second simulation shots, the abscissa of the four second simulation shots and the abscissa of the imaging point; and then determining the zero-offset travel time of the imaging point according to the zero-offset travel time of the imaging point along the abscissa direction, the ordinate of the four second simulation shots and the longitudinal and transverse coordinates of the imaging point.
104. The computer equipment determines the seismic wave travel time corresponding to the plurality of shot points and the plurality of detection points according to the plurality of seismic channel data, the seismic wave travel time of the plurality of simulated shot points and the zero-shot travel time of the plurality of imaging points.
In one possible implementation, this step may include the following steps (21) to (23):
(21) For each seismic trace data, the computer device determines a first location of a shot point and a second location of a pickup point corresponding to the seismic trace data.
In one possible implementation, one seismic trace data corresponds to a seismic wave signal emitted by one shot point received by one detector point, and the corresponding relationship between the seismic trace data and the shot point position and the detector point position is stored in the computer device. Correspondingly, the method comprises the following steps: for each seismic trace data, the computer device determines a first position of a shot and a second position of a geophone corresponding to the seismic trace data from stored correspondence of the seismic trace data to the shot and the geophone positions.
(22) The computer equipment determines travel time of the shot corresponding to the seismic channel data to a plurality of imaging points based on the first position, the seismic wave travel time corresponding to the plurality of simulated shots and the zero-shot travel time corresponding to the plurality of imaging points, and obtains the seismic wave travel time corresponding to the shot.
In one possible implementation, this step may be implemented by the following steps (a) to (c):
(a) The computer equipment determines a first target grid area where the shot corresponding to the seismic channel data is located based on the first position, and determines the travel time of the seismic waves corresponding to four first simulation shots forming the first target grid area from the travel time of the seismic waves corresponding to the simulation shots, wherein the travel time of the seismic waves corresponding to the simulation shots is used for representing the travel time of the seismic wave signals from the simulation shots to the imaging points.
In one possible implementation, the computer device stores a coordinate range of a plurality of grid areas included in the preset grid. Correspondingly, the computer equipment determines a first target grid area where the shot corresponding to the seismic channel data is based on the first position, and the steps are as follows: the computer equipment determines coordinates of a shot point corresponding to the seismic channel data based on the first position; determining a first target grid area corresponding to a coordinate range where the coordinates are located according to the coordinates of the shot points corresponding to the seismic channel data; the computer equipment determines four first simulation shots corresponding to the first target grid area, and determines the seismic wave travel time corresponding to the four first simulation shots from the seismic wave travel time corresponding to the simulation shots.
It should be noted that, the coordinates of the shot corresponding to the seismic trace data are plane coordinates, for example, the coordinates of the shot corresponding to the seismic trace data may be represented by (x, y), where x represents the plane abscissa of the shot corresponding to the seismic trace data and y represents the plane ordinate of the shot corresponding to the seismic trace data.
(b) For each imaging point, the computer equipment determines the target zero-offset travel time corresponding to the imaging point from the zero-offset travel time corresponding to the imaging points, and determines the target travel time from the first simulation shot to the imaging point from the seismic wave travel time corresponding to the first simulation shot for each first simulation shot, so as to obtain the target travel time of four first simulation shots.
In one possible implementation, the target survey location includes a plurality of bins, one bin corresponding to each imaging point location; the computer stores the corresponding relation between the positions of a plurality of imaging points and zero offset travel. Correspondingly, the method comprises the following steps: for each imaging point, the computer equipment determines the corresponding target zero-offset travel time of the imaging point from the corresponding relation between the position and the zero-offset travel time according to the target position of the imaging point.
It should be noted that, when the travel time of the seismic wave corresponding to the first simulated shot includes travel time from the first simulated shot to a plurality of imaging points, for each imaging point, the target travel time is used to represent travel time from the first simulated shot to the imaging point. In one possible implementation, the computer stores a correspondence between the location of the simulated shot, the location of the imaging point, and the travel time. Correspondingly, the computer equipment determines the target travel time corresponding to the four first simulation shots from travel time corresponding to the four first simulation shots by the steps of: and the computer equipment determines the target travel time corresponding to the first simulated shot points from the corresponding relation among the positions of the simulated shot points, the positions of the imaging points and the travel time according to the positions of the first simulated shot points and the positions of the imaging points for each first simulated shot point, and obtains the target travel time corresponding to the four first simulated shot points.
(c) The computer equipment determines travel time of the shot corresponding to the seismic channel data to the imaging point based on the first position, the target zero-shot travel time corresponding to the imaging point and the target travel time of the four first simulated shots, and obtains travel time of the shot to a plurality of imaging points.
In one possible implementation, this step is: the computer device determines the position of the imaging point and the positions of the four first simulated shots, determines four triangular areas based on the position of the imaging point and the positions of the four first simulated shots, and one triangular area comprises one imaging point and two adjacent first simulated shots; determining a target triangle area where the first position is located from the four triangle areas based on the first position, the position of the imaging point and the positions of the four first simulated shots; determining travel time from shot corresponding to the seismic channel data to an imaging point based on target travel time of two adjacent first simulation shot corresponding to the target triangle area and target zero-shot travel time corresponding to the imaging point; a travel time of the shot to a plurality of imaging points is obtained.
It should be noted that the imaging point may be located in the first grid area formed by the four first simulated shots, or may be located outside the first grid area. The first grid area comprises four sides, and the imaging point and one side form a triangle area to obtain four triangle areas.
Optionally, the step of determining the target triangle area where the first position is located from the four triangle areas by the computer device is: for each triangle area, the computer equipment determines target parameters according to a first position, the position of an imaging point and the positions of two adjacent first simulation shots corresponding to the triangle area through the following formula I; under the condition that the target parameter is a positive value, determining that a shot corresponding to the seismic channel data is in the triangular area, and determining that the triangular area is a target triangular area where the first position is located;
equation one: (a) 1 ,b 1 ,c 1 )＝(x-x 1 ,y-y 1 ,0)×(x-x 2 ,y-y 2 ,0)
(a 2 ,b 2 ,c 2 )＝(x-x 2 ,y-y 2 ,0)×(x-x 3 ,y-y 3 ,0)
(a 3 ,b 3 ,c 3 )＝(x-x 3 ,y-y 3 ,0)×(x-x 1 ,y-y 1 ,0)
A＝c 1 *c 2 *c 3
Wherein A represents a target parameter, (x, y) represents a first position, (x) 1 ,y 1 )、(x 2 ,y 2 ) Representing the positions of two adjacent first simulated shots, (x) 3 ,y 3 ) Representing the position of the imaging point.
It should be noted that, since the shots corresponding to the seismic trace data are located in the grid area formed by the four first simulated shots, and the four triangle areas are obtained by dividing the grid area, the shots corresponding to the seismic trace data are located in one three shape areas, that is, only one target triangle area, the corresponding value of a is positive, and the corresponding values of a are negative.
In the embodiment of the application, the target triangle area where the shot corresponding to the seismic trace data is located is determined by the position of the imaging point and the position of the first simulated shot, so that the distance from the first simulated shot to the imaging point is larger than that from the first position to the imaging point, the accuracy range corresponding to the target triangle area is ensured to be larger than that of the seismic wave travel time of the shot corresponding to the first position, and the accuracy of the determined seismic wave travel time is improved.
Optionally, the computer device determines the travel time of the seismic waves of the shot corresponding to the seismic channel data as follows: the computer equipment determines the seismic wave travel time of the shot corresponding to the seismic channel data based on the target travel time of two adjacent first simulation shots corresponding to the target triangle area and the target zero-shot travel time corresponding to the imaging point through the following formula II;
formula II: t=t 1 *c 1 +t 2 *c 2 +t 3 *c 3
c 1 ＝(ys 3 *x 23 -ys 3 *x 23 )/(y 13 *x 23 -x 13 *y 23 )
c 2 ＝(xs 3 *y 13 -ys 3 *x 13 )/(y 13 *x 23 -x 13 *y 23 )
c 3 ＝1-c 1 -c 2
Wherein t represents the travel time of the seismic waves of the shot corresponding to the seismic channel data, t 1 、t 2 Representing target travel time of two adjacent first simulation shots corresponding to the target triangle area; t is t 3 Representing the corresponding target zero-gun travel time of the imaging point.
In the embodiment of the application, the seismic wave travel time of the shot corresponding to the seismic channel data is determined according to the target travel time of the first simulated shot and the target zero-distance travel time of the imaging point, and the distance from the first simulated shot to the imaging point is larger than the distance from the shot corresponding to the seismic channel data to the imaging point, so that the target travel time of the first simulated shot is ensured to cover the seismic wave travel time of the shot, and the accuracy of the determined seismic wave travel time is further improved.
(23) The computer equipment determines travel time of the detection point corresponding to the seismic channel data to travel time of the imaging points based on the second position, travel time corresponding to the simulated shot points and zero-shot travel time corresponding to the imaging points, and obtains seismic wave travel time corresponding to the detection point.
It should be noted that, the positions of the detection points corresponding to the seismic trace data and the shots corresponding to the seismic trace data are only different, and the method of determining the travel corresponding to the plurality of detection points by the computer device is the same as the method of determining the travel corresponding to the plurality of shots, and will not be described in detail herein.
Another point to be noted is that the computer device determines offset imaging information of the work area based on the seismic wave travel time of the plurality of shots and the plurality of detectors after determining the travel time of the plurality of shots and the plurality of detectors. In one possible implementation, the step of determining offset imaging information for the work area based on the seismic travel time of the plurality of shots and the plurality of detectors by the computer device is: the computer equipment performs superposition imaging on the seismic wave travel time of the plurality of shot points and the seismic wave travel time of the plurality of detection points to obtain a superposition profile of the work area, and determines the superposition profile as offset imaging information of the work area. Optionally, the computer device performs superposition imaging on the seismic travel time of the plurality of shots, and performs superposition imaging on the seismic travel time of the plurality of detection points, so as to obtain a superposition cross-sectional view of the work area as shown in fig. 2.
In one possible implementation, the superimposed cross-section obtained when the seismic travel is determined by the prior art method is shown in fig. 3 for the same work area. As can be seen from a comparison of fig. 2 and 3, in shallow formation regions, the accuracy of the superimposed profile obtained by the present application is significantly higher than that obtained in the prior art.
The embodiment of the application provides a method for determining when a seismic wave travels, because when the seismic wave travels corresponding to a plurality of simulated shots and the seismic wave travels corresponding to a plurality of imaging points are zero-offset, the seismic wave travels corresponding to a plurality of shots and a plurality of detection points are determined, and when the zero-offset travels are the travel corresponding to the same plane coordinates of the imaging points and the simulated shots, the distance between the imaging points and the simulated shots is smaller than the distance from the exploration points such as the shots and the detection points to the imaging points, the seismic wave travel of the exploration points is avoided from exceeding the accuracy range of interpolation methods when traveling, and therefore the accuracy of the determined travel of the exploration point position is improved.
FIG. 4 is a block diagram illustrating a seismic travel determination apparatus according to an exemplary embodiment. Referring to fig. 4, the apparatus includes:
The acquisition module 401 is configured to acquire exploration data in a work area to be determined, where the exploration data includes a plurality of seismic trace data, and one seismic trace data corresponds to one shot point and one detector point of a plurality of shots and a plurality of detectors arranged in the work area;
a first determining module 402, configured to determine, according to the plurality of seismic trace data, a travel time of a seismic wave corresponding to a plurality of simulated shots on a preset grid in the work area, where the preset grid includes a plurality of grid areas, where the travel time of the seismic wave is used to represent a propagation time of the seismic wave signal from the simulated shots to a plurality of imaging points in the work area;
a second determining module 403, configured to determine, according to the seismic wave travel times corresponding to the plurality of simulated shots, a propagation duration of the seismic wave signal from the simulated shots to the imaging points when the zero-offset travel time corresponding to the plurality of imaging points in the work area is used to represent that the planar coordinates of the imaging points are the same as the planar coordinates of the simulated shots;
the third determining module 404 is configured to determine, according to the plurality of seismic trace data, the seismic wave travel time corresponding to the plurality of simulated shots, and the zero-offset travel time corresponding to the plurality of imaging points, the seismic wave travel time corresponding to the plurality of shots and the plurality of detectors.
In one possible implementation, the third determining module 404 is configured to determine, for each seismic trace data, a first position of a shot point and a second position of a detector point corresponding to the seismic trace data; and determining travel time of the shot corresponding to the seismic channel data to the plurality of imaging points based on the first position, the travel time of the seismic wave corresponding to the plurality of simulated shots and the travel time of the zero distances corresponding to the plurality of imaging points, obtaining the travel time of the seismic wave corresponding to the shot, and determining travel time of the detector corresponding to the seismic channel data to the travel time of the zero distances corresponding to the plurality of imaging points based on the second position, the travel time corresponding to the plurality of simulated shots and the travel time of the zero distances corresponding to the plurality of imaging points, and obtaining the travel time of the seismic wave corresponding to the detector.
In another possible implementation manner, the third determining module 404 is configured to determine, based on the first location, a first target grid area where the shot corresponding to the seismic trace data is located, and determine, from among the travel times of the seismic waves corresponding to the plurality of simulated shots, travel times of the seismic waves corresponding to the four first simulated shots that constitute the first target grid area, where the travel times of the seismic waves corresponding to the first simulated shots are used to represent travel times of the seismic wave signals from the first simulated shots to the plurality of imaging points; for each imaging point, determining a target zero-offset travel time corresponding to the imaging point from zero-offset travel times corresponding to a plurality of imaging points, and for each first simulation shot, determining a target travel time from the first simulation shot to an imaging point from a seismic wave travel time corresponding to the first simulation shot, so as to obtain target travel times of four first simulation shots; determining travel time of the shot corresponding to the seismic channel data to the imaging point based on the first position, the travel time of the target zero-shot corresponding to the imaging point and the travel time of the target of the four first simulation shots, and obtaining travel time of the shot to a plurality of imaging points.
In another possible implementation, the third determining module 404 is configured to determine a position of the imaging point and positions of the four first simulated shots, determine four triangle areas based on the position of the imaging point and the positions of the four first simulated shots, and determine one triangle area including one imaging point and two adjacent first simulated shots; determining a target triangle area where the first position is located from the four triangle areas based on the first position, the position of the imaging point and the positions of the four first simulated shots; and determining travel time from the shot corresponding to the seismic channel data to the imaging point based on the target travel time of the two adjacent first simulation shots corresponding to the target triangle area and the target zero-shot travel time corresponding to the imaging point.
In another possible implementation manner, the second determining module 403 is configured to determine zero-offset travel times corresponding to the multiple simulated shots from the seismic travel times of the multiple simulated shots; for each imaging point, determining a second target grid area where the imaging point is located, and determining zero-offset travel time corresponding to four second simulation shots forming the second target grid area from zero-offset travel time corresponding to a plurality of simulation shots; and determining zero-offset travel time of the imaging points based on the zero-offset travel time corresponding to the four second simulation shots, and obtaining the zero-offset travel time corresponding to the plurality of imaging points in the work area.
In another possible implementation, a first determining module 402 is configured to determine, from the plurality of seismic trace data, a velocity parameter of propagation of the seismic wave signal within a reservoir of the work area and an anisotropy parameter of the reservoir; and determining the travel time of the seismic waves corresponding to the plurality of simulated shots on the preset grid in the work area according to the speed parameter and the anisotropy parameter.
In another possible implementation, the apparatus further includes:
and the fourth determining module is used for determining offset imaging information of the work area based on the seismic waves corresponding to the plurality of shots and the plurality of detection points when traveling.
The embodiment of the application provides a determination device for seismic wave travel, because the seismic wave travel time corresponding to a plurality of simulated shots and the zero-offset travel time corresponding to a plurality of imaging points are used for determining the seismic wave travel time corresponding to a plurality of shots and a plurality of detection points, and the zero-offset travel time is the travel time corresponding to the same plane coordinates of the imaging points and the simulated shots, at the moment, the distance between the imaging points and the simulated shots is smaller than the distance from the exploration points such as the shots and the detection points to the imaging points, the seismic wave travel time of the exploration points is prevented from exceeding the accuracy range of interpolation methods during travel, and therefore the accuracy of the determined travel time of the exploration point position is improved.
Fig. 5 shows a block diagram of a computer device 500 according to an exemplary embodiment of the invention. The computer device 500 may be: a smart phone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III, motion picture expert compression standard audio plane 3), an MP4 (Moving Picture Experts Group Audio Layer IV, motion picture expert compression standard audio plane 4) player, a notebook computer, or a desktop computer. Computer device 500 may also be referred to by other names of user devices, portable computer devices, laptop computer devices, desktop computer devices, and the like.
In general, the computer device 500 includes: a processor 501 and a memory 502.
Processor 501 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and the like. The processor 501 may be implemented in at least one hardware form of DSP (Digital Signal Processing ), FPGA (Field-Programmable Gate Array, field programmable gate array), PLA (Programmable Logic Array ). The processor 501 may also include a main processor and a coprocessor, the main processor being a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit ); a coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 501 may integrate a GPU (Graphics Processing Unit, image processor) for rendering and drawing of content required to be displayed by the display screen. In some embodiments, the processor 501 may also include an AI (Artificial Intelligence ) processor for processing computing operations related to machine learning.
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, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 502 is used to store at least one instruction for execution by processor 501 to implement the method of determining when a seismic wave travels provided by a method embodiment in the present application.
In some embodiments, the computer device 500 may further optionally include: a peripheral interface 503 and at least one peripheral. The processor 501, memory 502, and peripheral interface 503 may be connected by buses or signal lines. The individual peripheral devices may be connected to the peripheral device interface 503 by buses, signal lines or circuit boards. Specifically, the peripheral device includes: at least one of radio frequency circuitry 504, a display 505, a camera 506, audio circuitry 507, a positioning component 508, and a power supply 509.
Peripheral interface 503 may be used to connect at least one Input/Output (I/O) related peripheral to processor 501 and memory 502. In some embodiments, processor 501, memory 502, and peripheral interface 503 are integrated on the same chip or circuit board; in some other embodiments, either or both of the processor 501, memory 502, and peripheral interface 503 may be implemented on separate chips or circuit boards, which is not limited in this embodiment.
The Radio Frequency circuit 504 is configured to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuitry 504 communicates with a communication network and other communication devices via electromagnetic signals. The radio frequency circuit 504 converts an electrical signal into an electromagnetic signal for transmission, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 504 includes: antenna systems, RF transceivers, one or more amplifiers, tuners, oscillators, digital signal processors, codec chipsets, subscriber identity module cards, and so forth. The radio frequency circuitry 504 may communicate with other computer devices via at least one wireless communication protocol. The wireless communication protocol includes, but is not limited to: metropolitan area networks, various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and/or WiFi (Wireless Fidelity ) networks. In some embodiments, the radio frequency circuitry 504 may also include NFC (Near Field Communication ) related circuitry, which is not limited in this application.
The display screen 505 is used to display a UI (user interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display 505 is a touch display, the display 505 also has the ability to collect touch signals at or above the surface of the display 505. The touch signal may be input as a control signal to the processor 501 for processing. At this time, the display 505 may also be used to provide virtual buttons and/or virtual keyboards, also referred to as soft buttons and/or soft keyboards. In some embodiments, the display 505 may be one, providing a front panel of the computer device 500; in other embodiments, the display 505 may be at least two, respectively disposed on different surfaces of the computer device 500 or in a folded design; in still other embodiments, the display 505 may be a flexible display disposed on a curved surface or a folded surface of the computer device 500. Even more, the display 505 may be arranged in a non-rectangular irregular pattern, i.e., a shaped screen. The display 505 may be made of LCD (Liquid Crystal Display ), OLED (Organic Light-Emitting Diode) or other materials.
The camera assembly 506 is used to capture images or video. Optionally, the camera assembly 506 includes a front camera and a rear camera. Typically, the front camera is disposed on a front panel of the computer device and the rear camera is disposed on a rear surface of the computer device. In some embodiments, the at least two rear cameras are any one of a main camera, a depth camera, a wide-angle camera and a tele camera, so as to realize that the main camera and the depth camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize a panoramic shooting and Virtual Reality (VR) shooting function or other fusion shooting functions. In some embodiments, camera assembly 506 may also include a flash. The flash lamp can be a single-color temperature flash lamp or a double-color temperature flash lamp. The dual-color temperature flash lamp refers to a combination of a warm light flash lamp and a cold light flash lamp, and can be used for light compensation under different color temperatures.
The audio circuitry 507 may include a microphone and a speaker. The microphone is used for collecting sound waves of users and environments, converting the sound waves into electric signals, and inputting the electric signals to the processor 501 for processing, or inputting the electric signals to the radio frequency circuit 504 for voice communication. The microphone may be provided in a plurality of different locations of the computer device 500 for stereo acquisition or noise reduction purposes. The microphone may also be an array microphone or an omni-directional pickup microphone. The speaker is used to convert electrical signals from the processor 501 or the radio frequency circuit 504 into sound waves. The speaker may be a conventional thin film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only the electric signal can be converted into a sound wave audible to humans, but also the electric signal can be converted into a sound wave inaudible to humans for ranging and other purposes. In some embodiments, audio circuitry 507 may also include a headphone jack.
The location component 508 is used to locate the current geographic location of the computer device 500 to enable navigation or LBS (Location Based Service, location-based services). The positioning component 508 may be a positioning component based on the United states GPS (Global Positioning System ), the Beidou system of China, the Granati system of Russia, or the Galileo system of the European Union.
The power supply 509 is used to power the various components in the computer device 500. The power supply 509 may be an alternating current, a direct current, a disposable battery, or a rechargeable battery. When the power supply 509 comprises a rechargeable battery, the rechargeable battery may support wired or wireless charging. The rechargeable battery may also be used to support fast charge technology.
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 acceleration sensor 511, a gyro sensor 512, a pressure sensor 513, a fingerprint sensor 514, an optical sensor 515, and a proximity sensor 516.
The acceleration sensor 511 can detect the magnitudes of accelerations on three coordinate axes of the coordinate system established with the computer device 500. For example, the acceleration sensor 511 may be used to detect components of gravitational acceleration on three coordinate axes. The processor 501 may control the display 505 to display a user interface in a landscape view or a portrait view according to a gravitational acceleration signal acquired by the acceleration sensor 511. The acceleration sensor 511 may also be used for acquisition of motion data of a game or a user.
The gyro sensor 512 may detect a body direction and a rotation angle of the computer device 500, and the gyro sensor 512 may collect a 3D motion of the user on the computer device 500 in cooperation with the acceleration sensor 511. The processor 501 may implement the following functions based on the data collected by the gyro sensor 512: motion sensing (e.g., changing UI according to a tilting operation by a user), image stabilization at shooting, game control, and inertial navigation.
The pressure sensor 513 may be disposed on a side frame of the computer device 500 and/or on an underlying layer of the display 505. When the pressure sensor 513 is disposed on the side frame of the computer device 500, a grip signal of the computer device 500 by a user may be detected, and the processor 501 performs left-right hand recognition or quick operation according to the grip signal collected by the pressure sensor 513. When the pressure sensor 513 is disposed at the lower layer of the display screen 505, the processor 501 controls the operability control on the UI interface according to the pressure operation of the user on the display screen 505. The operability controls include at least one of a button control, a scroll bar control, an icon control, and a menu control.
The fingerprint sensor 514 is used for collecting the fingerprint of the user, and the processor 501 identifies the identity of the user according to the fingerprint collected by the fingerprint sensor 514, or the fingerprint sensor 514 identifies the identity of the user according to the collected fingerprint. Upon recognizing that the user's identity is a trusted identity, the user is authorized by the processor 501 to perform relevant sensitive operations including unlocking the screen, viewing encrypted information, downloading software, paying for and changing settings, etc. The fingerprint sensor 514 may be provided on the front, back or side of the computer device 500. When a physical key or vendor Logo is provided on the computer device 500, the fingerprint sensor 514 may be integrated with the physical key or vendor Logo.
The optical sensor 515 is used to collect the ambient light intensity. In one embodiment, the processor 501 may control the display brightness of the display screen 505 based on the intensity of ambient light collected by the optical sensor 515. Specifically, when the intensity of the ambient light is high, the display brightness of the display screen 505 is turned up; when the ambient light intensity is low, the display brightness of the display screen 505 is turned down. In another embodiment, the processor 501 may also dynamically adjust the shooting parameters of the camera assembly 506 based on the ambient light intensity collected by the optical sensor 515.
A proximity sensor 516, also referred to as a distance sensor, is typically provided on the front panel of the computer device 500. The proximity sensor 516 is used to collect the distance between the user and the front of the computer device 500. In one embodiment, when the proximity sensor 516 detects a gradual decrease in the distance between the user and the front of the computer device 500, the processor 501 controls the display 505 to switch from the bright screen state to the off screen state; when the proximity sensor 516 detects that the distance between the user and the front of the computer device 500 gradually increases, the display 505 is controlled by the processor 501 to switch from the off-screen state to the on-screen state.
Those skilled in the art will appreciate that the architecture shown in fig. 5 is not limiting as to the computer device 500, and may include more or fewer components than shown, or may combine certain components, or employ a different arrangement of components.
Embodiments of the present application also provide a computer readable storage medium having at least one program code stored therein, the at least one program code being loaded and executed by a processor to implement operations performed in the method for determining travel of seismic waves in embodiments of the present application.
Embodiments of the present application also provide a computer program product comprising at least one piece of program code loaded and executed by a processor to perform the operations performed by the method for determining when a seismic wave travels in embodiments of the present application.
In some embodiments, the computer program related to the embodiments of the present application may be deployed to be executed on one computer device or on multiple computer devices located at one site, or on multiple computer devices distributed across multiple sites and interconnected by a communication network, where the multiple computer devices distributed across multiple sites and interconnected by a communication network may constitute a blockchain system.
The foregoing is merely for facilitating understanding of the technical solutions of the present application by those skilled in the art, and is not intended to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.

Claims (10)

1. A method of determining when a seismic wave travels, the method comprising:
acquiring exploration data in a work area to be determined, wherein the exploration data comprises a plurality of seismic channel data, and one seismic channel data corresponds to one shot point and one wave detection point of a plurality of shot points and a plurality of wave detection points which are distributed in the work area;
determining seismic wave travel time corresponding to a plurality of simulated shots on a preset grid in the work area according to the plurality of seismic channel data, wherein the preset grid comprises a plurality of grid areas, and the seismic wave travel time is used for representing the propagation time of seismic wave signals from the simulated shots to a plurality of imaging points in the work area;
determining zero-offset travel corresponding to a plurality of imaging points in the work area according to the seismic wave travel corresponding to the plurality of simulated shots, wherein the zero-offset travel is used for representing the propagation time of a seismic wave signal from the simulated shots to the imaging points when the plane coordinates of the imaging points and the simulated shots are the same;
and determining the seismic wave travel time corresponding to the plurality of shot points and the plurality of detection points according to the plurality of seismic channel data, the seismic wave travel time corresponding to the plurality of simulated shot points and the zero-offset travel time corresponding to the plurality of imaging points.
2. The method of claim 1, wherein determining the seismic travel time for the plurality of shots and the plurality of detectors based on the plurality of seismic trace data, the seismic travel time for the plurality of simulated shots, and the zero-offset travel time for the plurality of imaging points comprises:
for each seismic trace data, determining a first position of a shot point and a second position of a wave detection point corresponding to the seismic trace data;
determining travel time of the shot corresponding to the seismic channel data to travel time of the imaging points based on the first position, the travel time of the seismic waves corresponding to the simulated shots and the travel time of zero distances corresponding to the imaging points, obtaining the travel time of the seismic waves corresponding to the shots, and determining travel time of the detector corresponding to the seismic channel data to travel time of the imaging points based on the second position, the travel time of the simulated shots and the travel time of zero distances corresponding to the imaging points.
3. The method of claim 2, wherein the determining travel time of the shot corresponding to the seismic trace data to the plurality of imaging points based on the first location, the travel time of the seismic waves corresponding to the plurality of simulated shots, and the travel time of the zero shot corresponding to the plurality of imaging points comprises:
Determining a first target grid area where the shot corresponding to the seismic channel data is located based on the first position, and determining the travel time of the seismic waves corresponding to four first simulation shots forming the first target grid area from the travel time of the seismic waves corresponding to the simulation shots, wherein the travel time of the seismic waves corresponding to the first simulation shots is used for representing the travel time of the seismic wave signals from the first simulation shots to the imaging points;
for each imaging point, determining a target zero-offset travel time corresponding to the imaging point from zero-offset travel times corresponding to the imaging points, and for each first simulation shot, determining a target travel time from the first simulation shot to the imaging point from a seismic wave travel time corresponding to the first simulation shot, thereby obtaining target travel times of the four first simulation shots;
and determining travel time from the shot corresponding to the seismic channel data to the imaging point based on the first position, the target zero-shot travel time corresponding to the imaging point and the target travel time of the four first simulation shots, and obtaining travel time from the shot to the imaging points.
4. The method of claim 3, wherein determining travel of the shot corresponding to the seismic trace data to the imaging point based on the first location, the target zero-shot travel time corresponding to the imaging point, and the target travel time of the four first simulated shots comprises:
Determining the positions of the imaging point and the positions of the four first simulation shots, and determining four triangular areas based on the positions of the imaging point and the positions of the four first simulation shots, wherein one triangular area comprises one imaging point and two adjacent first simulation shots;
determining a target triangle area where the first position is located from the four triangle areas based on the first position, the position of the imaging point and the positions of the four first simulation shots;
and determining travel time from the shot corresponding to the seismic channel data to the imaging point based on the travel time of the target of the two adjacent first simulation shots corresponding to the target triangle area and the travel time of the target zero-offset corresponding to the imaging point.
5. The method of claim 1, wherein determining zero-offset travel time for a plurality of imaging points within the work area from the seismic travel time for the plurality of simulated shots comprises:
determining zero-offset travel time corresponding to the plurality of simulated shots from the seismic wave travel time of the plurality of simulated shots;
for each imaging point, determining a second target grid area where the imaging point is located, and determining zero-offset travel times corresponding to four second simulation shots composing the second target grid area from zero-offset travel times corresponding to the simulation shots;
And determining zero-offset travel time of the imaging point based on the zero-offset travel time corresponding to the four second simulation shot points, and obtaining the zero-offset travel time corresponding to a plurality of imaging points in the work area.
6. The method of claim 1, wherein determining, from the plurality of seismic trace data, travel of seismic waves corresponding to a plurality of simulated shots on a predetermined grid within the work area comprises:
determining a speed parameter of seismic wave signal propagation in a reservoir of the work area and an anisotropy parameter of the reservoir according to the plurality of seismic channel data;
and determining seismic wave travel time corresponding to a plurality of simulated shots on a preset grid in the work area according to the speed parameter and the anisotropy parameter.
7. The method according to claim 1, wherein the method further comprises:
and determining offset imaging information of the work area based on the plurality of shot points and the seismic wave travel time corresponding to the plurality of detection points.
8. A seismic travel time determination apparatus, the apparatus comprising:
the acquisition module is used for acquiring exploration data of a work area to be determined, wherein the exploration data comprises a plurality of seismic channel data, and one seismic channel data corresponds to one shot point and one detection point of a plurality of shot points and a plurality of detection points which are distributed in the work area;
The first determining module is used for determining seismic wave travel time corresponding to a plurality of simulated shots on a preset grid in the work area according to the plurality of seismic channel data, wherein the preset grid comprises a plurality of grid areas, and the seismic wave travel time is used for representing the propagation time of a seismic wave signal from the simulated shots to a plurality of imaging points in the work area;
the second determining module is used for determining zero-offset travel corresponding to a plurality of imaging points in the work area according to the seismic wave travel time corresponding to the plurality of simulated shots, and the zero-offset travel time is used for representing the propagation time length of a seismic wave signal from the simulated shots to the imaging points when the plane coordinates of the imaging points are the same as those of the simulated shots;
and the third determining module is used for determining the seismic wave travel time corresponding to the plurality of shot points and the detection points according to the plurality of seismic channel data, the seismic wave travel time corresponding to the plurality of simulated shot points and the zero-offset travel time corresponding to the plurality of imaging points.
9. A computer device, the computer device comprising:
a processor and a memory having stored therein at least one program code loaded and executed by the processor to perform the operations performed in the method of determining travel of seismic waves of any of claims 1 to 7.
10. A computer readable storage medium having stored therein at least one program code loaded and executed by a processor to perform the operations performed in the method of determining when a seismic wave is traveling as claimed in any one of claims 1 to 7.
CN202111356613.5A 2021-11-16 2021-11-16 Method and device for determining travel time of seismic waves, computer equipment and storage medium Pending CN116136605A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111356613.5A CN116136605A (en) 2021-11-16 2021-11-16 Method and device for determining travel time of seismic waves, computer equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111356613.5A CN116136605A (en) 2021-11-16 2021-11-16 Method and device for determining travel time of seismic waves, computer equipment and storage medium

Publications (1)

Publication Number Publication Date
CN116136605A true CN116136605A (en) 2023-05-19

Family

ID=86332987

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111356613.5A Pending CN116136605A (en) 2021-11-16 2021-11-16 Method and device for determining travel time of seismic waves, computer equipment and storage medium

Country Status (1)

Country Link
CN (1) CN116136605A (en)

Similar Documents

Publication Publication Date Title
CN111754386B (en) Image area shielding method, device, equipment and storage medium
CN109886208B (en) Object detection method and device, computer equipment and storage medium
CN112922589B (en) Pinch-out line determining method, pinch-out line determining device, terminal and storage medium
CN110174699B (en) Method, device and storage medium for determining fault closure
CN109696706B (en) Salt structure detection method, device and storage medium
CN111754564B (en) Video display method, device, equipment and storage medium
CN112882094B (en) First-arrival wave acquisition method and device, computer equipment and storage medium
CN110673214B (en) Method and device for predicting depths of inlet target point and end point of horizontal well
CN109343113B (en) Method and device for predicting oil and gas reservoir position
CN116136605A (en) Method and device for determining travel time of seismic waves, computer equipment and storage medium
CN113032712B (en) Determination method and device for drilling coincidence information, computer equipment and storage medium
CN112698391B (en) Sand body describing method, device, equipment and storage medium
CN110058306B (en) Method, apparatus and computer storage medium for determining three-dimensional velocity volume
CN112329355B (en) Method and device for determining single-well control area, computer equipment and storage medium
CN113009573B (en) Seismic data acquisition method, device, computer equipment and storage medium
CN111650637B (en) Seismic horizon interpretation method and device
CN115263244B (en) Method and device for controlling drainage of air-water layer and computer storage medium
CN112096373B (en) Method and device for selecting shale gas area and computer storage medium
CN112347604B (en) Method and device for determining vehicle path set
CN112214645B (en) Method and device for storing track data
CN109116415B (en) Seismic wave data separation method, device and storage medium
CN114254687A (en) Method, device and equipment for determining matching degree of drilling track and storage medium
CN115496828A (en) Method and device for determining stratum inclination angle parameters, computer equipment and storage medium
CN116136609A (en) Method and device for determining cloudization parameters, computer equipment and storage medium
CN115390157A (en) Method, device, equipment and storage medium for determining uniform parameters of exploration system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination