Common midpoint first arrival gather generation method, and device for calculating residual static correction

Abstract

The application discloses a common midpoint first arrival gather generating method and a residual static correction method and device. The application grids the position of each common midpoint, determines the order of each first arrival seismic trace in the common midpoint, and counts the number of the first arrival seismic trace of each bin where the common midpoint is located and the storage position of the first first arrival seismic trace of each bin. Thus, when generating the common midpoint first arrival gather file, the storage position of the first first arrival seismic trace of each bin and the order thereof can be used to quickly determine the offset of the address of all the first arrival seismic traces of each bin, so that the data writing of all the first arrival seismic traces of each common midpoint corresponding to each bin is quickly completed, and the common midpoint first arrival gather file is generated. The application shortens the processing time of the first arrival seismic trace data, improves the generation efficiency of the common midpoint first arrival gather, and better meets the time efficiency requirement of seismic data processing.

Description

Technical Field

[0001] This invention relates to the field of geophysical exploration data processing technology, and in particular to a method for generating common center point initial arrival gathers, a method and apparatus for calculating residual static correction. Background Technology

[0002] In seismic data, multiple shot-receiver pairs share the center point of the shot and the receiver; such a center point is called a common center point.

[0003] Common center point first arrival gathers are a commonly used data volume in seismic data processing. Field data is common shot point gathers. Therefore, common center point first arrival gathers need to be generated automatically based on the SPS file (which includes shot point file, receiver point file, and relationship file defining the shot point-receiver point relationship).

[0004] The simplest method is to search for each shot point, extract relevant information, and then generate the common center point initial arrival gather. However, due to the large number of shot points, this method has the largest computational load, especially in three dimensions, where the computational load is astonishing and the computation time is long.

[0005] Currently, no publicly available methods for generating common centroid first-arrival gathers have been found in the literature. Therefore, based on the above simple analysis, the computational workload for generating three-dimensional common centroid first-arrival gathers is enormous, requiring a considerable amount of time on a regular microcomputer, making it difficult to meet the timeliness requirements of seismic data processing in actual production. Summary of the Invention

[0006] In view of the above problems, the present invention is proposed to provide a method for generating common center point initial arrival set, a method and apparatus for calculating residual static correction to overcome or at least partially solve the above problems.

[0007] In a first aspect, embodiments of the present invention provide a method for generating a common center point initial arrival set, comprising:

[0008] Based on the azimuth of the survey line, the first coordinates of each receiver point and each shot point are rotated to obtain the second coordinates in a new coordinate system with the survey line as due east or due north.

[0009] According to the preset area size, the new coordinate system is divided into grids, and the area where the common center point of each detector point-shot point pair is located is determined according to the second coordinates of each detector point and each shot point.

[0010] Determine the common center point of the first arrival seismic trace for each shot, and the positional order of the first arrival seismic trace among all first arrival seismic traces sharing the common center point;

[0011] Based on the identifier of the surface cell containing the common center point of the first arrival seismic traces of each shot, count the total number of first arrival seismic traces contained in the common center point of each grid surface cell.

[0012] Based on the total number of first arrival seismic traces contained in the common center point of each facet, and in the order of each facet in the grid, the first first arrival seismic trace of the common center point of each facet is assigned a record position in the common center point first arrival trace set file as the starting record position of the common center point.

[0013] For each first-arrival seismic trace, the actual recording position of each first-arrival seismic trace in the common center point first-arrival trace set file is determined according to the starting recording position of the common center point where the first-arrival seismic trace is located and the position order of the first-arrival seismic trace in the common center point. The relevant information of the first-arrival seismic trace is written into the actual recording position to generate the common center point first-arrival trace set file.

[0014] Furthermore, based on the azimuth of the survey line, the first coordinates of each receiver point and each shot point are rotated to obtain the second coordinates in a new coordinate system with the survey line as due east or due north, specifically including:

[0015] Rotate the first coordinate in a clockwise direction;

[0016] The second coordinate is obtained as follows:

[0017] X R =(X-X0)*cos(90-θ)+(Y-Y0)*sin(90-θ)+X0

[0018] Y R =(Y-Y0)*cos(90-θ)-(X-X0)*sin(90-θ)+Y0

[0019] in:

[0020] X R The coordinates are the east coordinates in the second coordinate system obtained after rotation.

[0021] Y R The north coordinate is the second coordinate obtained after rotation;

[0022] X is the east coordinate in the first coordinate system before rotation;

[0023] Y is the north coordinate in the first coordinate system before rotation;

[0024] X0 is the east coordinate of the specified origin;

[0025] Y0 is the north coordinate of the specified origin;

[0026] θ is the azimuth angle of the survey line.

[0027] Furthermore, based on the azimuth of the survey line, the first coordinates of each receiver point and each shot point are rotated to obtain the second coordinates in a new coordinate system with the survey line as due east or due north, specifically including:

[0028] Rotate the first coordinate in a counterclockwise direction;

[0029] The second coordinate obtained is:

[0030] X R =(X-X0)*cos(θ)-(Y-Y0)*sin(θ)+X0

[0031] Y R =(Y-Y0)*cos(θ)+(X-X0)*sin(θ)+Y0

[0032] in:

[0033] X R The coordinates are the east coordinates in the second coordinate system obtained after rotation.

[0034] Y R The north coordinate is the second coordinate obtained after rotation;

[0035] X is the east coordinate in the first coordinate system before rotation;

[0036] Y is the north coordinate in the first coordinate system before rotation;

[0037] X0 is the east coordinate of the specified origin;

[0038] Y0 is the north coordinate of the specified origin;

[0039] θ is the azimuth angle of the survey line.

[0040] Furthermore, based on a preset element size, the new coordinate system is divided into a grid, specifically including:

[0041] According to the preset cell size of the observation system, the new coordinate system where the second coordinates of each detector point and each shot point are located is divided into a grid. The smallest unit of the grid obtained after division is a cell.

[0042] Following the order of horizontal and vertical coordinates from smallest to largest, starting with the cell with the smallest horizontal and vertical coordinates, each cell is labeled sequentially using the line number followed by the point number.

[0043] Furthermore, the surface element containing the common center point is determined in the following way:

[0044] If the rotated survey line direction is due north, the line number and point number of the surface element containing the common center point of the first arrival seismic trace are as follows:

[0045] i = Round((RX R +SX R ) / 2-cmpXzbmin) / ΔIz)+1;

[0046] j = Round((RY R +SY R ) / 2-cmpYzbmin) / ΔCz)+1;

[0047] in:

[0048] i represents the line number of the surface element containing the common center point, with a value ranging from 1 to M;

[0049] j is the point number for calculating the common center point element, and its value ranges from 1 to N;

[0050] Round is the rounding function;

[0051] RX R The east coordinate in the second coordinate system is obtained by rotating the detector point;

[0052] SX R The east coordinate in the second coordinate system is obtained after rotating the gun point;

[0053] RY R The north coordinate is the second coordinate obtained after rotating the detector point.

[0054] SY R The north coordinate is the second coordinate obtained after rotating the shot point;

[0055] cmpXzbmin is the east coordinate of the center of the line number-point number (1,1) element;

[0056] cmpYzbmin is the north coordinate of the center of the (1,1) surface element, which is the line number - point number.

[0057] ΔIz is the length of one surface element in the Inline direction of the survey line in the new coordinate system after rotation;

[0058] ΔCz is the length of a surface element in the Crossline direction, which is perpendicular to the Inline direction of the survey line, in the new coordinate system after rotation.

[0059] Furthermore, the identifier of the surface element containing the common center point of the initial arrival seismic trace for each shot is determined in the following manner:

[0060] If the rotated survey line direction is due east, the line number and point number of the surface element containing the common center point of the initial arrival seismic trace are as follows:

[0061] i = Round((RY R +SYR ) / 2-cmpYzbmin) / ΔIz)+1;

[0062] j = Round((RX R +SX R ) / 2-cmpXzbmin) / ΔCz)+1;

[0063] in:

[0064] i represents the line number of the surface element containing the common center point, with a value ranging from 1 to M;

[0065] j is the point number for calculating the common center point element, and its value ranges from 1 to N;

[0066] Round is the rounding function;

[0067] RX R The east coordinate in the second coordinate system is obtained by rotating the detector point;

[0068] SX R The east coordinate in the second coordinate system is obtained after rotating the gun point;

[0069] RY R The north coordinate is the second coordinate obtained after rotating the detector point.

[0070] SY R The north coordinate is the second coordinate obtained after rotating the shot point;

[0071] cmpXzbmin is the east coordinate of the center of the line number-point number (1,1) element;

[0072] cmpYzbmin is the north coordinate of the center of the (1,1) surface element, which is the line number - point number.

[0073] ΔIz is the length of one surface element in the Inline direction of the survey line in the new coordinate system after rotation;

[0074] ΔCz is the length of a surface element in the Crossline direction, which is perpendicular to the Inline direction of the survey line, in the new coordinate system after rotation.

[0075] Furthermore, for the first first-arrival seismic trace at the common center point of each surface element, its recording position in the common center point first-arrival trace set file is assigned, including:

[0076] The record position of the first first arrival seismic trace for each common center point in the common center point first arrival trace set file is calculated as follows:

[0077] T 1,1 =1;

[0078] When i = 1, j > 1;

[0079] When i>1, j=1;

[0080] When i>1, j>1;

[0081] in:

[0082] T 1,1 This indicates the record position of the first first arrival seismic trace of the common center point of the (1,1) surface element in the common center point first arrival trace set file.

[0083] T i,j This indicates the record position of the first first arrival seismic trace of the common center point of the surface element with line number-point number (i,j) in the common center point first arrival trace set file;

[0084] C l,k This represents the total number of first-arrival seismic traces at the common center point of the (l,k) surface element, where the line number is minus the point number.

[0085] i and l represent the line number of the face cell containing the common center point. i takes values ​​from 1 to M, l takes values ​​from 1 to i-1, and M represents the maximum line number of the face cell in the grid.

[0086] j and k represent the point numbers of the face cells containing the common center point. j takes values ​​from 1 to N, k takes values ​​from 1 to j-1, and N represents the maximum point number of the face cells in the grid.

[0087] Further, according to the starting recording position of the common center point where the first-arrival seismic traces are located and the positional order of the first-arrival seismic traces within the common center point, the actual recording position of each first-arrival seismic trace in the common center point first-arrival trace set file is determined, specifically including:

[0088] The actual recorded position is calculated using the following formula:

[0089] R i,j =T i,j +N s,i,j -1;

[0090] Where: R i,j This indicates the actual record location of the relevant information of the first arrival seismic trace in the common center point first arrival gather file;

[0091] T i,j This indicates the starting record position of the common center point of the surface element with line number and point number (i,j);

[0092] N s,i,jThis indicates the positional order of the common center point (i,j) among the surface elements with the initial arrival seismic trace online number and point number (i,j);

[0093] S is the firing point number.

[0094] Furthermore, before rotating the first coordinates of each receiver point and each shot point according to the azimuth angle of the survey line, the method further includes:

[0095] Create an index for each shot point and each receiver point in the order of shot point and receiver point;

[0096] Based on the seismic data collected in the field, the first arrival time is picked up to obtain the first arrival seismic trace data.

[0097] Secondly, embodiments of the present invention provide a method for calculating the residual static correction of the first arrival wave. The common center point first arrival gather used in the method for calculating the residual static correction of the first arrival wave is obtained by the common center point first arrival gather generation method as described above.

[0098] Thirdly, embodiments of the present invention provide an apparatus for generating a common center point initial arrival set, comprising:

[0099] The rotation module is used to rotate the first coordinates of each receiver point and each shot point according to the azimuth angle of the survey line, so as to obtain the second coordinates in a new coordinate system with the survey line as due east or due north.

[0100] The mesh generation module is used to divide the new coordinate system into a mesh shape according to the preset cell size, and to determine the cell containing the common center point of each detector point-shot point pair according to the second coordinates of each detector point and each shot point.

[0101] The seismic trace location sequence determination module is used to determine the common center point of the first arrival seismic trace of each shot, and the position sequence of the first arrival seismic trace among all first arrival seismic traces sharing the common center point;

[0102] The total number of traces determination module is used to count the total number of first-arrival traces contained in the common center point of each grid cell based on the identifier of the cell containing the common center point of the first-arrival traces of each shot.

[0103] The initial position determination module is used to assign the first initial arrival seismic trace of the common center point of each facet element to the record position in the common center point initial arrival trace set file according to the total number of first arrival seismic traces contained in the common center point of each facet element in the order of each facet element in the grid, and to the common center point initial arrival trace as the starting record position of the common center point.

[0104] The common center point first arrival gather file generation module is used to determine the actual recording position of each first arrival seismic trace in the common center point first arrival gather file according to the starting recording position of the common center point where the first arrival seismic trace is located and the position order of the first arrival seismic trace in the common center point, and write the relevant information of the first arrival seismic trace into the actual recording position to generate the common center point first arrival gather file.

[0105] Fourthly, embodiments of the present invention provide an apparatus for calculating the residual static correction of the first arrival wave. The common center point first arrival gather used by the apparatus for calculating the residual static correction of the first arrival wave is obtained by the common center point first arrival gather generation method described above.

[0106] Fifthly, embodiments of the present invention provide a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements either the aforementioned method for generating common center point first arrival gathers or the aforementioned method for calculating the residual static correction of first arrival waves.

[0107] In a sixth aspect, embodiments of the present invention provide a computer-readable storage medium storing a computer program that performs the aforementioned method for generating common center point first arrival sets, or a computer program that implements the aforementioned method for calculating the residual static correction of first arrival waves.

[0108] The beneficial effects of the above-mentioned technical solution provided by the present invention include at least the following:

[0109] The method, residual static correction calculation method, and apparatus for generating common center point first-arrival gather files provided by this invention grid the positions of each common center point, determine the order of each first-arrival seismic trace within that common center point, and statistically count the number of first-arrival seismic traces in each surface cell containing the common center point, as well as the storage location of the first first-arrival seismic trace in each surface cell. Thus, when generating the common center point first-arrival gather file, the storage location of the first first-arrival seismic trace in each surface cell, along with its own order, can quickly determine the address offsets of all first-arrival seismic traces in each surface cell, thereby rapidly completing the writing of data for all first-arrival seismic traces corresponding to each common center point in each surface cell, thereby generating the common center point first-arrival gather file. This invention shortens the processing time for first-arrival seismic trace data, improves the generation efficiency of common center point first-arrival gathers, and better meets the timeliness requirements of seismic data processing.

[0110] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.

[0111] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0112] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings:

[0113] Figure 1 This is a flowchart of the method for generating the common center point initial arrival set in an embodiment of the present invention;

[0114] Figure 2A-2C This is a schematic diagram of the original coordinate system and the new coordinate system after rotation under different rotation methods in embodiments of the present invention;

[0115] Figure 3 This is a schematic diagram illustrating the mesh generation process for the new coordinate system after rotation, as described in an embodiment of the present invention.

[0116] Figure 4 This is a schematic diagram showing that the common center point of different shot-receiver pairs is located within a surface element in an embodiment of the present invention;

[0117] Figure 5A and Figure 5B This is a schematic diagram of the surface element line numbering and point numbering representing the common center point in an embodiment of the present invention;

[0118] Figure 6 This is a schematic diagram showing the total number of first-arrival seismic traces included at each common center point in an embodiment of the present invention;

[0119] Figure 7 This is a schematic diagram showing the storage location of the first arrival seismic traces contained in each common center point in an embodiment of the present invention;

[0120] Figure 8 This is a bar chart comparing the time taken by prior art and that of the present invention in the embodiments of the present invention;

[0121] Figure 9 This is a structural block diagram of the device for generating the common center point initial arrival set in an embodiment of the present invention. Detailed Implementation

[0122] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0123] To shorten the generation time of common center point initial arrival gather files and improve the processing efficiency of common center point initial arrival gather file data, this invention provides a method for generating common center point initial arrival gathers, referring to... Figure 1 As shown, the method includes the following steps:

[0124] S11. Based on the azimuth of the survey line, rotate the first coordinates of each receiver point and each shot point to obtain the second coordinates in a new coordinate system with the survey line as due east or due north.

[0125] S12. According to the preset area size, the new coordinate system is divided into grids, and the area where the common center point of each detector point and each shot point is located is determined according to the second coordinates of each detector point and each shot point.

[0126] S13. Determine the common center point of the first arrival seismic trace of each shot, and the positional order of the first arrival seismic trace among all first arrival seismic traces sharing the common center point;

[0127] S14. Based on the identifier of the surface cell containing the common center point of the first arrival seismic traces of each shot, count the total number of first arrival seismic traces contained in the common center point of each grid surface cell.

[0128] S15. Based on the total number of first arrival seismic traces contained in the common center point of each facet, and in the order of each facet in the grid, the first first arrival seismic trace of the common center point of each facet is assigned a record position in the common center point first arrival trace set file as the starting record position of the common center point.

[0129] S16. For each first-arrival seismic trace, determine the actual recording position of each first-arrival seismic trace in the common center point first-arrival trace set file according to the starting recording position of the common center point where the first-arrival seismic trace is located and the position order of the first-arrival seismic trace in the common center point, and write the relevant information of the first-arrival seismic trace into the actual recording position to generate the common center point first-arrival trace set file.

[0130] The method for generating common center point first-arrival gather files provided in this invention involves gridding the positions of each common center point, determining the order of each first-arrival seismic trace within that common center point, and statistically analyzing the number of first-arrival seismic traces in each facet cell containing that common center point, as well as the storage location of the first first-arrival seismic trace in each facet cell. Thus, when generating the common center point first-arrival gather file, the storage location of the first first-arrival seismic trace in each facet cell, along with its own order, allows for the rapid determination of the address offsets of all first-arrival seismic traces in each facet cell. This enables the rapid writing of data for all first-arrival seismic traces corresponding to each common center point in each facet cell, thereby generating the common center point first-arrival gather file. This invention shortens the processing time for first-arrival seismic trace data, improves the generation efficiency of common center point first-arrival gathers, and better meets the timeliness requirements of seismic data processing.

[0131] Before performing the above steps S11 to S16, this embodiment of the invention also includes arranging SPS files in sequence, establishing indexes for receiver points and shot points, and picking up first arrival times to generate first arrival seismic trace data.

[0132] Creating an index for the receiver and shot points involves arranging the receiver and shot point files in ascending order, and creating an index for each receiver and shot point based on its sequential position among all receivers and shot points. For example, the shot points are numbered as: Shot Point 1, Shot Point 2, Shot Point 3... and the receivers are numbered as: Receiver Point 1, Receiver Point 2, Receiver Point 3...

[0133] The process of picking first arrival times and generating first arrival seismic traces involves using field-collected seismic data and employing either manual or automated picking methods to capture the first arrival times, ensuring both the quality and quantity of the data captured, and then generating the first arrival seismic traces.

[0134] Furthermore, for ease of explanation, the original coordinates of each geophone point and each shot point in the actual field coordinate system are called the "first coordinates", and the coordinates obtained after rotation are called the "second coordinates".

[0135] In step S11 above, based on the azimuth of the survey line, the first coordinates of each receiver point and each shot point are rotated to obtain the second coordinates in a new coordinate system with the survey line as due east or due north. In practice, this can be divided into two cases: clockwise rotation and counterclockwise rotation. These are explained below:

[0136] Method 1:

[0137] Rotate the first coordinate in a clockwise direction;

[0138] The second coordinate is obtained as follows:

[0139] X R =(X-X0)*cos(90-θ)+(Y-Y0)*sin(90-θ)+X0

[0140] Y R =(Y-Y0)*cos(90-θ)-(X-X0)*sin(90-θ)+Y0

[0141] in:

[0142] X R The coordinates are the east coordinates in the second coordinate system obtained after rotation.

[0143] Y R The north coordinate is the second coordinate obtained after rotation;

[0144] X is the east coordinate in the first coordinate system before rotation;

[0145] Y is the north coordinate in the first coordinate system before rotation;

[0146] X0 is the east coordinate of the specified origin;

[0147] Y0 is the north coordinate of the specified origin;

[0148] θ is the azimuth angle of the survey line.

[0149] Method 2:

[0150] Rotate the first coordinate in a counterclockwise direction;

[0151] The second coordinate obtained is:

[0152] X R =(X-X0)*cos(θ)-(Y-Y0)*sin(θ)+X0

[0153] Y R =(Y-Y0)*cos(θ)+(X-X0)*sin(θ)+Y0

[0154] in:

[0155] X R The coordinates are the east coordinates in the second coordinate system obtained after rotation.

[0156] Y R The north coordinate is the second coordinate obtained after rotation;

[0157] X is the east coordinate in the first coordinate system before rotation;

[0158] Y is the north coordinate in the first coordinate system before rotation;

[0159] X0 is the east coordinate of the specified origin;

[0160] Y0 is the north coordinate of the specified origin;

[0161] θ is the azimuth angle of the survey line.

[0162] Reference Figures 2A-2C As shown, Figure 2A This is a schematic diagram of the original coordinate system before rotation. Figure 2B This is a schematic diagram of the coordinate system after clockwise rotation. Figure 2C The diagram shows the coordinate system after counterclockwise rotation. Figures 2A-2C The diagram contains many shot points and receiver points (distinguished by different depths in the diagram).

[0163] Furthermore, in step S12 above, the new coordinate system is divided into grids according to the preset cell size. Specifically, the new coordinate system in which the second coordinates of each receiver point and each shot point are located is divided into grids according to the preset cell size of the observation system. The smallest unit of the grid obtained after division is a cell. Starting from the cell with the smallest horizontal and vertical coordinates, each cell is identified in sequence according to the line numbering-point numbering method, in order of ascending horizontal and vertical coordinates.

[0164] The purpose of step S11 is to enable the rotated new coordinate system to more quickly locate the position of the facet containing the common center point. That is, by calculating the position of the facet and its corresponding number using the east and north coordinates, a common center point can be quickly located (each common center point will only fall into one facet, that is, it will only correspond to one facet). This improves the calculation speed and efficiency.

[0165] In order for the surface elements to quickly locate the common center point, the distribution density of the common center point formed by each receiver-shot point pair needs to be fully considered during mesh generation to determine the size of the surface elements.

[0166] For example, the new coordinate system after rotation can be meshed according to the specification of surface size (10 meters in the inline direction and 20 meters in the crossline direction perpendicular to the inline). The mesh diagram after meshing can be referenced. Figure 3 As shown. The face elements are square grids (rectangles or squares), and their specific size is not limited in this embodiment of the invention.

[0167] Reference Figure 4 As shown, the shared center point (common center point) of shot-receiver pairs from different directions and with different offset distances is located within a certain surface element.

[0168] After meshing, each facet needs to be numbered. Each facet represents a common midpoint; therefore, numbering the facets is equivalent to numbering the common midpoints within those facets. The numbering of the facet containing each common midpoint can be determined as follows:

[0169] (i) If the direction of the survey line in the new coordinate system after rotation is due north, the line number and point number of the surface element containing the common center point of the initial arrival seismic trace are as follows:

[0170] i = Round((RX R +SX R ) / 2-cmpXzbmin) / ΔIz)+1;

[0171] j = Round((RY R +SY R ) / 2-cmpYzbmin) / ΔCz)+1;

[0172] in:

[0173] i represents the line number of the surface element containing the common center point, with a value ranging from 1 to M;

[0174] j is the point number for calculating the common center point element, and its value ranges from 1 to N;

[0175] Round is the rounding function;

[0176] RX R The east coordinate in the second coordinate system is obtained by rotating the detector point;

[0177] SX R The east coordinate in the second coordinate system is obtained after rotating the gun point;

[0178] RY R The north coordinate is the second coordinate obtained after rotating the detector point.

[0179] SY R The north coordinate is the second coordinate obtained after rotating the shot point;

[0180] cmpXzbmin is the east coordinate of the center of the line number-point number (1,1) element;

[0181] cmpYzbmin is the north coordinate of the center of the (1,1) surface element, which is the line number - point number.

[0182] ΔIz is the length of one surface element in the inline direction (i.e., the direction of the survey line) of the new coordinate system after rotation;

[0183] ΔCz is the length of a surface element in the Crossline direction, which is perpendicular to the Inline direction of the survey line, in the new coordinate system after rotation.

[0184] (II) If the direction of the survey line in the new coordinate system after rotation is due east, the line number and point number of the surface element containing the common center point of the initial arrival seismic trace are as follows:

[0185] i = Round((RY R +SY R ) / 2-cmpYzbmin) / ΔIz)+1;

[0186] j = Round((RX R +SX R ) / 2-cmpXzbmin) / ΔCz)+1;

[0187] in:

[0188] i represents the line number of the surface element containing the common center point, with a value ranging from 1 to M;

[0189] j is the point number for calculating the common center point element, and its value ranges from 1 to N;

[0190] Round is the rounding function;

[0191] RX R The east coordinate in the second coordinate system is obtained by rotating the detector point;

[0192] SX R The east coordinate in the second coordinate system is obtained after rotating the gun point;

[0193] RY R The north coordinate is the second coordinate obtained after rotating the detector point.

[0194] SY R The north coordinate is the second coordinate obtained after rotating the shot point;

[0195] cmpXzbmin is the east coordinate of the center of the line number-point number (1,1) element;

[0196] cmpYzbmin is the north coordinate of the center of the (1,1) surface element, which is the line number - point number.

[0197] ΔIz is the length of one surface element in the Inline direction of the survey line in the new coordinate system after rotation;

[0198] ΔCz is the length of a surface element in the Crossline direction, which is perpendicular to the Inline direction of the survey line, in the new coordinate system after rotation.

[0199] Reference Figure 5A As shown, the numbering of each face element (representing a common center point) consists of a line number and a point number. Figure 5A Dashed lines indicate the order of line numbers, while solid lines indicate the order of point numbers. On the same line number, there are several different points to represent different face elements. Similarly, on the same solid line with a different line number, the intersections with dashed lines representing different face elements form line numbers for those face elements. The range of line numbers is 1 to M, and the range of point numbers is 1 to N.

[0200] Starting with the cell with the smallest x and y coordinates in the grid, number the cells containing the common center point according to (line number, point number). The line number starts from 1 and increments by 1, where M is the largest line number of the cell in the grid; the point number on each line starts from 1 and increments by 1, where N is the largest point number on each line in the grid.

[0201] Of course, in specific implementations of this invention, the line numbers and point numbers do not necessarily start from 1. (Refer to...) Figure 5BAs shown, each face element is numbered by 7 digits, starting from 1001001. The leftmost column of face elements is numbered from 1001001 to 1001010, with the first four digits representing the dot and line numbers and the last three digits representing the dot numbers.

[0202] After numbering is complete, the first-arrival seismic traces at the common center point of each area cell can be counted to determine the total number of first-arrival seismic traces contained within each common center point (i.e., within the area cell). For an example of the total number of seismic traces contained within each common center point, please refer to [link to example]. Figure 6 As shown.

[0203] For each shot with a first-arrival seismic trace, calculate the line number and point number of the common center point of that first-arrival seismic trace in the common center point grid, and simultaneously accumulate and record the total number of first-arrival seismic traces at that common center point.

[0204] Reference Figure 4 In the example shown, the common center point contains 5 first-arrival traces (assuming one first-arrival trace per shot).

[0205] The common center point element is calculated first according to the shot point order. During the calculation, the order of each shot-receiver pair in that element must also be recorded, such as... Figure 4 In the diagram, shot point 1 and receiver point 1 indicate that the shot-receiver pair is the first element of the array, shot point 2 and receiver point 2 indicate that the shot-receiver pair is the second element of the array, and so on. After the array is assigned a starting position in the storage file, the storage position of the shot-receiver pair in the file can be quickly calculated based on that starting position. This is an important factor in improving the calculation speed in this embodiment of the invention.

[0206] The aforementioned starting position, which is the record position of the first first arrival seismic trace of the common center point of the surface element (the first arrival seismic trace of the smallest shot-receiver pair when the shot point numbers are sorted in ascending order) in the common center point first arrival trace set file.

[0207] Further, in step S15 above, the recording position of the first first arrival seismic trace of each common center point in the common center point first arrival trace set file is calculated and assigned in the following manner:

[0208] T 1,1 =1;

[0209] (When i = 1, j > 1);

[0210] (When i>1, j=1);

[0211] (When i>1, j>1);

[0212] in:

[0213] T 1,1 This indicates the record position of the first first arrival seismic trace of the common center point of the (1,1) surface element in the common center point first arrival trace set file.

[0214] T i,j This indicates the record position of the first first arrival seismic trace of the common center point of the surface element with line number-point number (i,j) in the common center point first arrival trace set file;

[0215] C l,k This represents the total number of first-arrival seismic traces at the common center point of the (l,k) surface element, where the line number is minus the point number.

[0216] i and l represent the line number of the face cell containing the common center point. i takes values ​​from 1 to M, l takes values ​​from 1 to i-1, and M represents the maximum line number of the face cell in the grid.

[0217] j and k represent the point numbers of the face cells containing the common center point. j takes values ​​from 1 to N, k takes values ​​from 1 to j-1, and N represents the maximum point number of the face cells in the grid.

[0218] Based on the aforementioned initial position, all first-arrival seismic trace data of the common center point can be quickly written to the corresponding position in the first-arrival gather file according to the offset.

[0219] Specifically, the actual recording position in S16 above can be calculated, for example, using the following formula:

[0220] R i,j =T i,j +N s,i,j -1;

[0221] Where: R i,j This indicates the actual record location of the relevant information of the first arrival seismic trace in the common center point first arrival gather file;

[0222] T i,j This indicates the starting record position of the common center point of the surface element with line number and point number (i,j);

[0223] N s,i,j This indicates the positional order of the common center point (i,j) among the surface elements with the initial arrival seismic trace online number and point number (i,j);

[0224] S is the firing point number.

[0225] Reference Figure 7 As shown, based on the statistics of the number of first-arrival seismic traces for each facet element, the traces are then allocated to storage locations (addresses) in the common center point first-arrival trace set file in descending order of facet element (common center point) line number - point number. For example... Figure 7 1-5 represent Figure 6 The first arrival seismic trace data of the 1001001 surface element is stored in locations 1-5 in the common center point first arrival trace gathering file, and so on. 137-148 indicates that the first arrival seismic trace data of the 1002005 surface element is stored in locations 137-148 in the common center point first arrival trace gathering file.

[0226] First arrival seismic trace data includes: the index of the receiver point of the first arrival seismic trace, the index of the shot point of the first arrival seismic trace, the shot-receiver distance of the first arrival seismic trace, and the first arrival time of the first arrival seismic trace.

[0227] The following example briefly illustrates the method for generating the common center point initial arrival set. The process of this example includes the following steps:

[0228] 1) Arrange the SPS in sequence and create an index;

[0229] The indexes of the first geophone and shot points are 1, and the indices of the subsequent geophone and shot points are incremented by 1.

[0230] 2) Pick up the initial arrival time

[0231] Based on the seismic data collected in the field, the first arrival time is picked up to obtain the first arrival seismic trace data.

[0232] 3) Coordinate rotation

[0233] Step 3) The coordinate rotation is performed based on the specified coordinate origin and clockwise rotation direction, with the survey line azimuth angle being 37.50. The coordinates of all receiver points and shot points are rotated to obtain a new coordinate system with due east as the abscissa and due north as the ordinate. The direction of the survey line after rotation is due north.

[0234] The formula for clockwise rotation of coordinates is:

[0235] X R =(X-X0)*cos(90-37.5)+(Y-Y0)*sin(90-37.5)+X0

[0236] Y R =(Y-Y0)*cos(90-37.5)-(X-X0)*sin(90-37.5)+Y0

[0237] Where: X R The coordinates are the rotated east coordinates;

[0238] Y R The north coordinates are the rotated coordinates;

[0239] X is the east coordinate before rotation;

[0240] Y is the north coordinate before rotation;

[0241] X0 is the east coordinate of the specified origin;

[0242] Y0 is the north coordinate of the specified origin;

[0243] 4) Arrangement of elements with common center point

[0244] Step 4) describes the common center point arrangement. Based on the preset cell size of the observation system (10 meters in the inline direction and 20 meters in the crossline direction), the common center points are divided into grids according to the new coordinates obtained in Step 3). Starting with the cell with the smallest x and y coordinates in the grid, the cells corresponding to the common center point are numbered according to (line number, point number). Line numbers start from 1, incrementing by 1, with 650 being the maximum line number in the common center point grid. Point numbers on each line start from 1, incrementing by 1, with 1200 being the maximum point number on each line in the common center point grid.

[0245] 5) Initialization of data upon arrival at the common center point

[0246] The initial value of the number of first arrival traces of earthquakes contained in the surface element containing each common center point is assigned to 0.

[0247] 6) Calculation of common center point element information

[0248] When the rotated survey line is due north, the formulas for calculating the line number and point number of the common center point of the first arrival seismic trace are as follows:

[0249] i = Round((RX R +SX R ) / 2-cmpXzbmin) / 10)+1

[0250] j = Round((RY R +SY R ) / 2-cmpYzbmin) / 20)+1

[0251] Where: i is the line number of the grid for calculating the common center point, and its value ranges from 1 to 650;

[0252] j is the point number for calculating the common center point grid, and its value ranges from 1 to 1200;

[0253] Round is the rounding function;

[0254] RX R The coordinates of the receiver point after rotation are the east coordinates.

[0255] SX R The eastern coordinates of the firing point after rotation;

[0256] RY R The coordinates of the receiver point after rotation are the north coordinates.

[0257] SY R The north coordinates of the firing point after rotation;

[0258] cmpXzbmin is the east coordinate of the center of the line number-point number (1,1) element;

[0259] cmpYzbmin is the north coordinate of the center of the (1,1) surface element, which is the line number - point number.

[0260] 7) Allocation of common center point cell addresses

[0261] In the area information calculation results obtained in step 6), the starting record position of the first arrival trace data of each common center point in the common center point first arrival file is calculated. Starting from area (1,1), the starting position of the first arrival data of this grid in the soon-to-be-generated common center point first arrival file is calculated for all grids in line number-point number order. The calculation method is as follows:

[0262] T 1,1 =1

[0263] (When i = 1, j > 1)

[0264] (When i>1, j=1)

[0265] (When i>1, j>1);

[0266] Wherein: T 1,1 This indicates the record position of the first first arrival seismic trace of the common center point of the (1,1) surface element in the common center point first arrival trace set file.

[0267] T i,j This indicates the record position of the first first arrival seismic trace of the common center point of the surface element with line number-point number (i,j) in the common center point first arrival trace set file;

[0268] C l,k This represents the total number of first-arrival seismic traces at the common center point of the (l,k) surface element, where the line number is minus the point number.

[0269] i and l represent the line numbers of the face elements containing the common center point. i takes values ​​from 1 to M, and l takes values ​​from 1 to i-1. Here, M = 650.

[0270] j and k represent the point numbers of the surface elements containing the common center point, with j ranging from 1 to N and k ranging from 1 to j-1. Here, N = 1200.

[0271] 8) Generate common center point initial arrival file

[0272] Starting with the first shot, for each seismic trace with a first arrival, the relevant information for that trace is written into the common center point first arrival file. The relevant information includes: the index of the receiver point of this seismic trace obtained in step 1), the index of the shot point of this seismic trace obtained in step 1), the shot-receiver distance of this seismic trace, and the first arrival time of this seismic trace.

[0273] 9) Application of common center point initial arrival file

[0274] The generated common center point first arrival data file can be used to calculate the residual static correction of the first arrival wave or to monitor the effect of static correction application.

[0275] Figure 8 A comparison chart shows the time taken to generate a 3D common center point first-arrival seismic gather file for a target area using existing methods and the method provided in this embodiment of the invention, both running on the same microcomputer. The comparison... Figure 8 As can be seen above, the existing method requires 90 hours, while the method of this invention only requires 3 hours. Therefore, the running speed is increased by 30 times after adopting the method provided by this invention, which can well meet the time requirements of three-dimensional processing.

[0276] This invention also provides a method for calculating the residual static correction of the first arrival wave. The common center point first arrival gather used in this method for calculating the residual static correction of the first arrival wave is obtained by the method for generating the common center point first arrival gather as provided in the foregoing embodiments.

[0277] The specific implementation process for calculating the residual static correction of the first arrival wave can refer to the existing technology, and the embodiments of the present invention do not limit it.

[0278] Based on the same inventive concept, this invention also provides an apparatus for generating common center point first arrival gathers and an apparatus for calculating the residual static correction of first arrival waves. Since the principles by which these apparatuses solve problems are similar to the aforementioned methods for generating common center point first arrival gathers and calculating the residual static correction of first arrival waves, the implementation of these apparatuses can refer to the implementation of the aforementioned methods, and repeated details will not be repeated.

[0279] This invention provides an apparatus for generating common center point initial arrival gathers, referring to... Figure 9 As shown, it includes:

[0280] The rotation module 91 is used to rotate the first coordinates of each receiver point and each shot point according to the azimuth angle of the survey line to obtain the second coordinates in a new coordinate system with the survey line as due east or due north.

[0281] The meshing module 92 is used to divide the new coordinate system into a mesh shape according to the preset surface size, and to determine the surface element where the common center point of each detector point-shot point pair is located according to the second coordinates of each detector point and each shot point.

[0282] The seismic trace location sequence determination module 93 is used to determine the common center point where the first arrival seismic trace of each shot is located, and the position sequence of the first arrival seismic trace among all first arrival seismic traces sharing the common center point;

[0283] The total number of traces determination module 94 is used to count the total number of first-arrival traces contained in the common center point of each grid cell based on the identifier of the cell containing the common center point of the first-arrival traces of each shot.

[0284] The initial position determination module 95 is used to assign the first initial arrival seismic trace of the common center point of each facet to the record position in the common center point initial arrival trace set file according to the total number of first arrival seismic traces contained in the common center point of each facet, in the order of each facet in the grid, as the starting record position of the common center point.

[0285] The common center point first arrival gather file generation module 96 is used to determine the actual recording position of each first arrival seismic trace in the common center point first arrival gather file according to the starting recording position of the common center point where the first arrival seismic trace is located and the position order of the first arrival seismic trace in the common center point, and write the relevant information of the first arrival seismic trace into the actual recording position to generate the common center point first arrival gather file.

[0286] This invention also provides an apparatus for calculating the residual static correction of the first arrival wave. The common center point first arrival gather used by the apparatus for calculating the residual static correction of the first arrival wave is obtained by the common center point first arrival gather generation method described above.

[0287] This invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements either the method for generating the common center point first arrival gather as described above or the method for calculating the residual static correction of the first arrival wave as described above.

[0288] This invention also provides a computer-readable storage medium storing a computer program that performs the method for generating common center point first arrival gathers as described above, or a computer program that implements the method for calculating the residual static correction of first arrival waves as described above.

[0289] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage and optical storage) containing computer-usable program code.

[0290] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0291] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0292] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0293] Obviously, those skilled in the art can make various modifications and variations to this invention without departing from its spirit and scope. Therefore, if these modifications and variations fall within the scope of the claims of this invention and their equivalents, this invention also intends to include these modifications and variations.

Claims

1. A method for generating common midpoint first arrival gathers, characterized in that, The method comprises the following steps: According to the azimuth angle of the survey line, the first coordinates of each geophone and each shot point are rotated to obtain second coordinates in a new coordinate system with the survey line as the positive east or the positive north; According to the preset bin size, the new coordinate system is grid-divided, and according to the second coordinates of each geophone and each shot point, the bin in which the common midpoint of each geophone-shot point pair is located is determined; The common midpoint of the first arrival seismic trace of each shot is determined, and the position sequence of the first arrival seismic trace in all first arrival seismic traces sharing the common midpoint is determined; According to the identification of the bin in which the common midpoint of each first arrival seismic trace is located, the total number of first arrival seismic traces contained by the common midpoint of each grid bin is counted; According to the total number of first arrival seismic traces contained by the common midpoint in each bin, the first first arrival seismic trace of the common midpoint of each bin is sequentially assigned its recording position in the common midpoint first arrival trace gather file as the starting recording position of the common midpoint, according to the sequence of each bin in the grid; For each first arrival seismic trace, the actual recording position of each first arrival seismic trace in the common midpoint first arrival trace gather file is determined according to the starting recording position of the common midpoint in which the first arrival seismic trace is located and the position sequence of the first arrival seismic trace in the common midpoint, and the relevant information of the first arrival seismic trace is written into the actual recording position to generate the common midpoint first arrival trace gather file.

2. The method of claim 1, wherein, According to the azimuth angle of the survey line, the first coordinates of each geophone and each shot point are rotated to obtain second coordinates in a new coordinate system with the survey line as the positive east or the positive north, specifically comprising: The first coordinates are rotated in the clockwise direction; The second coordinates obtained are: X R = (X - X0) * cos(90 - Θ) + (Y - Y0) * sin(90 - Θ) + X0 Y R = (Y - Y0) * cos(90 - Θ) - (X - X0) * sin(90 - Θ) + Y0 Wherein: X R is the east coordinate in the second coordinate obtained after rotation; Y R is the north coordinate in the second coordinate obtained after rotation; X is the east coordinate in the first coordinates before rotation; Y is the north coordinate in the first coordinates before rotation; X0 is the east coordinate of the specified origin; Y0 is the north coordinate of the specified origin; θ is the azimuth angle of the survey line.

3. The method of claim 1, wherein, According to the azimuth angle of the survey line, the first coordinates of each geophone and each shot point are rotated to obtain second coordinates in a new coordinate system with the survey line as the positive east or the positive north, specifically comprising: The first coordinates are rotated in the counterclockwise direction; The second coordinates obtained are: X R = (X - X0) * cos(θ) - (Y - Y0) * sin(θ) + X0 Y R = (Y - Y0) * cos(θ) + (X - X0) * sin(θ) + Y0 Wherein: X R is the east coordinate in the second coordinate obtained after rotation; Y R is the north coordinate in the second coordinate obtained after rotation; X is the east coordinate in the first coordinates before rotation; Y is the north coordinate in the first coordinates before rotation; X0 is the east coordinate of the specified origin; Y0 is the north coordinate of the specified origin; θ is the azimuth angle of the survey line.

4. The method of claim 1, wherein, According to the preset bin size, the new coordinate system is grid-divided, specifically comprising: According to the preset bin size of the observation system, the new coordinate system in which the second coordinates of each geophone and each shot point are located is grid-divided, and the smallest unit of the grid obtained after the division is a bin; According to the order from small to large of the horizontal and vertical coordinates, the bins are sequentially identified in the order of line number-point number from the bin with the smallest horizontal and vertical coordinates.

5. The method of claim 4, wherein, The bin in which the common midpoint is located is determined in the following manner: If the direction of the rotated survey line is the positive north direction, the line number and the point number of the bin in which the common midpoint of the first arrival seismic trace is located are: i = Round((RX R + SX R ) / 2 - cmpXzbmin) / Δlz) + 1; j = Round((RY R + SY R / 2 - cmpYzbmin) / ΔCz) + 1; Wherein: i is the line number of the bin in which the common midpoint is calculated, and the value range is from 1 to M; j is the point number of the common midpoint bin, and the value range is from 1 to N; Round is a rounding function; RX R East coordinate in the second coordinate system after rotating the detection point; SX R is the east coordinate in the second coordinate system after the shot point is rotated; RY R N is the north coordinate in the second coordinate obtained after rotating the detection point; SY R N is the north coordinate in the second coordinate obtained after the shot point is rotated; cmpXzbmin is the east coordinate of the center of the bin with line number 1 and point number 1; cmpYzbmin is the north coordinate of the center of the bin with line number 1 and point number 1; ΔIz is the length of a bin in the direction of the survey line Inline in the new coordinate system after rotation; ΔCz is the length of a bin in the direction of the cross line Crossline perpendicular to the survey line Inline in the new coordinate system after rotation.

6. The method of claim 4, wherein, The identification of the bin in which the common midpoint of each first arrival seismic trace is located is determined in the following manner: If the survey line direction after rotation is the positive east direction, the line number and the point number of the bin in which the common midpoint of the first arrival seismic trace is located are respectively: i = Round((RY R + SY R / 2 - cmpYzbmin) / Δlz) + 1; j = Round((RX R + SX R ) / 2 - cmpXzbmin) / ACz) + 1; wherein: i is the line number of the bin in which the common midpoint is located, and the value range is from 1 to M; j is the point number of the common midpoint bin, and the value range is from 1 to N; Round is a rounding function; RX R is the east coordinate in the second coordinate system after the rotation of the detection point; SX R is the east coordinate in the second coordinate system after the shot point is rotated; RY R N is the north coordinate in the second coordinate obtained after rotating the detection point; SY R N is the north coordinate in the second coordinate obtained after the shot point is rotated; cmpXzbmin is the east coordinate of the center of the bin with line number 1 and point number 1; cmpYzbmin is the north coordinate of the center of the bin with line number 1 and point number 1; ΔIz is the length of a bin in the direction of the survey line Inline in the new coordinate system after rotation; ΔCz is the length of a bin in the direction of the cross line Crossline perpendicular to the survey line Inline in the new coordinate system after rotation.

7. The method of claim 4, wherein, For each first arrival seismic trace of the common midpoint, its recording position in the common midpoint first arrival trace gather file is assigned, including: The recording position of each first arrival seismic trace of the common midpoint in the common midpoint first arrival trace gather file is calculated in the following manner: T 1,1 ＝1； when i = 1, j > 1; when i > 1, j = 1 ; when i > 1, j > 1; wherein: T 1,1 denotes the record position of the first arrival seismic trace of the common center point of the (1,1) bin in the common center point first arrival trace gather file; T i,j denotes the record position of the first arrival seismic trace of the common center point of the (i,j) bin in the common center point first arrival trace gather file; C l,k N (l, k) represents the total number of the first seismic traces of the common center point of the (l, k) surface element; i and l represent the line number of the bin in which the common midpoint is located, i takes a value from 1 to M, and l takes a value from 1 to i-1, M representing the maximum line number of the bin in the grid; j and k represent the point number of the bin in which the common midpoint is located, j takes a value from 1 to N, and k takes a value from 1 to j-1, N representing the maximum point number of the bin in the grid.

8. The method of claim 4, wherein, The actual recording position of each first arrival seismic trace in the common midpoint first arrival trace gather file is determined according to the starting recording position of the common midpoint in which the first arrival seismic trace is located and the position sequence of the first arrival seismic trace in the common midpoint, and specifically includes: The actual recording position is calculated by the following formula: R i,j = T i,j + N s,i,j - 1 ; wherein: R i,j represents the actual recording position of the relevant information of the first arrival seismic trace in the common center point first arrival gather file; T i,j represents the starting record position of the co-center point of the (i, j) bin; N s,i,j denotes the position order in the common point (i,j) of the first arrival seismic trace online number-point number (i,j) bin; S is the shot point number.

9. The method according to any one of claims 1 to 8, wherein, Before rotating the first coordinates of each geophone and each shot point according to the azimuth of the survey line, the method further includes: indexes are respectively established for each shot point and each geophone according to the order of the shot points and the geophones; According to the seismic data collected in the field, a first arrival time picking operation is performed to obtain the data of the first arrival seismic trace.

10. A method of calculating first break residual statics, characterized by, The common midpoint first arrival trace gather used in the method for calculating the residual static correction of the first arrival wave is obtained by using the method for generating the common midpoint first arrival trace gather according to any one of claims 1-9.

11. A device for generating a common center point initial arrival set, characterized in that, It comprises: a rotating module, configured to rotate the first coordinates of each geophone and each shot point according to the azimuth of the survey line to obtain second coordinates in a new coordinate system in which the survey line is the positive east or the positive north; a grid dividing module, configured to divide the new coordinate system into grids according to a preset bin size, and determine a bin where a common point of each receiver-point pair is located according to the second coordinates of each receiver and each shot point; a seismic trace position sequence determining module, configured to determine a common point where a first arrival seismic trace of each shot is located, and a position sequence of the first arrival seismic trace in all first arrival seismic traces sharing the common point; a total trace number determining module, configured to count a total number of first arrival seismic traces contained by a common point in each grid bin according to an identifier of the bin where the common point of each shot is located; an initial position determining module, configured to assign a recording position of a first first arrival seismic trace of each common point in a grid bin in a common point first arrival trace gather file to the common point as a starting recording position of the common point according to the total number of first arrival seismic traces contained by the common point in each bin, and according to an order of the bins in the grid; a common point first arrival trace gather file generating module, configured to determine an actual recording position of each first arrival seismic trace in a common point first arrival trace gather file according to the starting recording position of the common point where the first arrival seismic trace is located and the position sequence of the first arrival seismic trace in the common point, and write related information of the first arrival seismic trace into the actual recording position to generate the common point first arrival trace gather file.

12. An apparatus for computing first break residual statics, characterized in that, The common point first arrival trace gather used by the device for calculating first arrival wave residual static correction is obtained by using the method for generating a common point first arrival trace gather according to any one of claims 1 to 9.

13. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor executes the computer program to implement the method for generating a common point first arrival trace gather according to any one of claims 1 to 9 or the method for calculating first arrival wave residual static correction according to claim 10.

14. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a computer program for executing the method for generating a common point first arrival trace gather according to any one of claims 1 to 9 or a computer program for implementing the method for calculating first arrival wave residual static correction according to claim 10.

Images