CN116413774A - Common-center-point first-arrival gather generation method, calculation residual static correction method and device - Google Patents
Common-center-point first-arrival gather generation method, calculation residual static correction method and device Download PDFInfo
- Publication number
- CN116413774A CN116413774A CN202111651249.5A CN202111651249A CN116413774A CN 116413774 A CN116413774 A CN 116413774A CN 202111651249 A CN202111651249 A CN 202111651249A CN 116413774 A CN116413774 A CN 116413774A
- Authority
- CN
- China
- Prior art keywords
- arrival
- point
- common
- center point
- bin
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000012937 correction Methods 0.000 title claims abstract description 30
- 230000003068 static effect Effects 0.000 title claims abstract description 30
- 238000004364 calculation method Methods 0.000 title abstract description 15
- 238000003860 storage Methods 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 abstract description 14
- 230000032683 aging Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 12
- 230000001174 ascending effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a common-center point first arrival gather generation method, a residual static correction calculation method and a residual static correction calculation device. The method and the system mesh the positions of the common-center points, determine the sequence of each first-arrival seismic channel in the common-center point, count the number of the first-arrival seismic channels of the surface element where each common-center point is positioned and the storage position of the first-arrival seismic channel of each surface element according to the sequence, and thus, when the common-center-point first-arrival gather file is generated, the storage position of the first-arrival seismic channel of the surface element and the sequence of the first-arrival seismic channel can be utilized, the offset of the addresses of all the first-arrival seismic channels of each surface element can be rapidly determined, and therefore, the writing of the data of all the first-arrival seismic channels of the common-center points corresponding to each surface element is rapidly completed, and the common-center-point first-arrival gather file is generated. The method shortens the processing time of the first-arrival seismic channel data, improves the generation efficiency of the first-arrival gather of the common center point, and better meets the aging requirement of seismic data processing.
Description
Technical Field
The invention relates to the technical field of geophysical exploration data processing, in particular to a common-center point first-arrival gather generation method, a calculation residual static correction method and a calculation residual static correction device.
Background
In seismic data, a plurality of shot-geophone pairs share the center points of shot and geophone, such center points being referred to as co-centerpoints.
The common-center first-arrival gather is a data body commonly used for seismic data processing, and the field data is the common-shot gather, so that the common-center first-arrival gather needs to be generated autonomously according to SPS files (including shot files, detector files and relation files defining shot-detector relations).
The simplest generation method is to search from shots one by one, extract related information, and then generate a first arrival gather of common center points, but the calculated amount is the largest because of the large number of shots, especially three-dimensional, the calculated amount is quite striking, and the calculated time is long.
Currently, in the published literature, no published method of generating a common-center first-arrival gather is retrieved. Therefore, according to the simple analysis, the calculation amount for generating the three-dimensional common-center first arrival gathers is huge, the time required on a common microcomputer is quite long, and the aging requirement of seismic data processing is difficult to meet in actual production.
Disclosure of Invention
The present invention has been made in view of the above-mentioned problems, and has as its object to provide a common-center first-arrival gather generation method, a calculation residual static correction method and an apparatus which overcome or at least partially solve the above-mentioned problems.
In a first aspect, an embodiment of the present invention provides a method for generating a first arrival gather at a common center point, including:
according to the azimuth angle of the measuring line, rotating the first coordinates of each detecting point and each shot point to obtain the second coordinates in a new coordinate system taking the measuring line as the positive east or the positive north;
grid division is carried out on the new coordinate system according to the preset bin size, and bins where common center points of all detection points-shot point pairs are located are determined according to second coordinates of all detection points and all shot points;
determining a common center point of a first arrival seismic channel of each shot and the position sequence of the first arrival seismic channel in all the first arrival seismic channels sharing the common center point;
counting the total channel number of the first arrival seismic channels contained in the common center point of each grid surface element according to the identification of the surface element where the common center point of each first arrival seismic channel is located;
according to the total channel number of the first-arrival seismic channels contained in the common-center point of each surface element, distributing the recording position of the first-arrival seismic channel in the common-center point first-arrival channel set file as the initial recording position of the common-center point according to the sequence of each surface element in the grid;
For each first-arrival seismic channel, determining the actual recording position of each first-arrival seismic channel in a first-arrival gather file of the common center point according to the initial recording position of the common center point where the first-arrival seismic channel is located and the position sequence of the first-arrival seismic channel in the common center point, and writing the related information of the first-arrival seismic channel into the actual recording position to generate the first-arrival gather file of the common center point.
Further, according to the azimuth of the measuring line, the first coordinates of each detecting point and each shot point are rotated to obtain the second coordinates in the new coordinate system taking the measuring line as the positive east or the positive north, and the method specifically comprises the following steps:
rotating the first coordinate in a clockwise direction;
the second coordinates are obtained as follows:
X R ＝(X-X 0 )*cos(90-θ)+(Y-Y 0 )*sin(90-θ)+X 0
Y R ＝(Y-Y 0 )*cos(90-θ)-(X-X 0 )*sin(90-θ)+Y 0
wherein:
X R the east coordinate in the second coordinate obtained after rotation;
Y R the north coordinate in the second coordinate obtained after rotation;
x is the east coordinate in the first coordinate before rotation;
y is the north coordinate of the first coordinate before rotation;
X 0 an east coordinate for the specified origin;
Y 0 north coordinates for a specified origin;
θ is the azimuth of the line.
Further, according to the azimuth of the measuring line, the first coordinates of each detecting point and each shot point are rotated to obtain the second coordinates in the new coordinate system taking the measuring line as the positive east or the positive north, and the method specifically comprises the following steps:
Rotating the first coordinate in a counterclockwise direction;
the second coordinates obtained are:
X R ＝(X-X 0 )*cos(θ)-(Y-Y 0 )*sin(θ)+X 0
Y R ＝(Y-Y 0 )*cos(θ)+(X-X 0 )*sin(θ)+Y 0
wherein:
X R the east coordinate in the second coordinate obtained after rotation;
Y R the north coordinate in the second coordinate obtained after rotation;
x is the east coordinate in the first coordinate before rotation;
y is the north coordinate of the first coordinate before rotation;
X 0 an east coordinate for the specified origin;
Y 0 north coordinates for a specified origin;
θ is the azimuth of the line.
Further, according to a preset bin size, the new coordinate system is divided into grid shapes, which specifically includes:
according to the size of a bin preset by an observation system, performing grid division on a new coordinate system where the second coordinates of each detection point and each shot point are located, wherein the minimum unit of the grid obtained after division is a bin;
the method comprises the steps of marking each bin in a line number-point number mode sequentially from the bin with the smallest abscissa to the bin with the smallest ordinate according to the descending order of the abscissa.
Further, the bin in which the common center point is located is determined by:
if the direction of the rotated measuring line is the north direction, the line number and the point number of the bin where the common center point of the first arrival seismic channel is located are respectively:
i＝Round((RX R +SX R )/2-cmpXzbmin)/ΔIz)+1；
j＝Round((RY R +SY R )/2-cmpYzbmin)/ΔCz)+1；
Wherein:
i is the line number of the bin where the calculated common center point is located, and the value range is from 1 to M;
j is the point number of the calculated common-center point face element, and the value range is from 1 to N;
round is a rounding function;
RX R the east coordinate in the second coordinate is obtained after the wave detection point rotates;
SX R the east coordinates in the second coordinates are obtained after the shot point rotates;
RY R the north coordinate in the second coordinate obtained after the rotation of the wave detector;
SY R the north coordinates in the second coordinates obtained after the shot point rotates;
cmpXzbmin is the line number-the point number is the east coordinate of the (1, 1) bin center;
cmpYzbmin is the line number-the point number is the north coordinate of the (1, 1) bin center;
ΔIz is the length of one bin of the line in the line direction in the new coordinate system after rotation;
Δcz is the length of one bin in the Crossline direction perpendicular to the Inline direction in the new coordinate system after rotation.
Further, the identification of the bin where the common center point of the first arrival seismic trace of each shot is located is determined by the following method:
if the direction of the rotated survey line is the forward direction, the line number and the point number of the bin where the common center point of the first arrival seismic channel is located are respectively:
i＝Round((RY R +SY R )/2-cmpYzbmin)/ΔIz)+1；
j＝Round((RX R +SX R )/2-cmpXzbmin)/ΔCz)+1；
wherein:
i is the line number of the bin where the calculated common center point is located, and the value range is from 1 to M;
j is the point number of the calculated common-center point face element, and the value range is from 1 to N;
round is a rounding function;
RX R the east coordinate in the second coordinate is obtained after the wave detection point rotates;
SX R the east coordinates in the second coordinates are obtained after the shot point rotates;
RY R the north coordinate in the second coordinate obtained after the rotation of the wave detector;
SY R the north coordinates in the second coordinates obtained after the shot point rotates;
cmpXzbmin is the line number-the point number is the east coordinate of the (1, 1) bin center;
cmpYzbmin is the line number-the point number is the north coordinate of the (1, 1) bin center;
ΔIz is the length of one bin of the line in the line direction in the new coordinate system after rotation;
Δcz is the length of one bin in the Crossline direction perpendicular to the Inline direction in the new coordinate system after rotation.
Further, for a first arrival seismic trace at a common center point of each bin, assigning its recording location in the common center point first arrival gather file, comprising:
the recording position of the first arrival seismic channel of each common center point in the common center point first arrival gather file is calculated according to the following mode:
T 1,1 ＝1；
wherein:
T 1,1 representing the recording position of a first arrival seismic trace of a common center point of a line number-point number (1, 1) bin in the common center point first arrival gather file;
T i,j Representing the recording position of a first arrival seismic trace of a common center point of a line number-point number (i, j) bin in the common center point first arrival gather file;
C l,k representing the total number of first-arrival seismic traces with line numbers-points numbered as the common center point of the (l, k) bin;
i and l represent line numbers of the cells where the common center point is located, i takes values from 1 to M, l takes values from 1 to i-1, and M represents the maximum line number of the cells in the grid.
j and k represent the point numbers of the cells where the common center point is located, j takes values from 1 to N, k takes values from 1 to j-1, and N represents the maximum point number of the cells in the grid.
Further, determining an actual recording position of each first-arrival seismic trace in the common-center-point first-arrival gather file according to the initial recording position of the common-center point where the first-arrival seismic trace is located and the position sequence of the first-arrival seismic trace in the common-center point, specifically including:
the actual recording position is calculated by the following formula:
R i,j ＝T i,j +N s,i,j -1；
wherein: r is R i,j Representing the actual recording position of the related information of the first-arrival seismic channel in the first-arrival gather file of the common center point;
T i,j a start recording position representing a line number-point number of a common center point of the (i, j) bin;
N s,i,j representing the order of locations in the common center point (i, j) of the first-arrival seismic trace on-line number-point number (i, j) bin;
S is the shot number.
Further, before rotating the first coordinates of each detector point and each shot point according to the azimuth of the survey line, the method further includes:
according to the sequence of the shots and the wave detectors, respectively establishing indexes for each shot and each wave detector;
and carrying out pickup operation of the first arrival time according to the seismic data collected in the field to obtain data of the first arrival seismic channel.
In a second aspect, an embodiment of the present invention provides a method for calculating a residual static correction of a first arrival wave, where a common-center first arrival gather used in the method for calculating a residual static correction of a first arrival wave is obtained by using a method for generating a common-center first arrival gather as described above.
In a third aspect, an embodiment of the present invention provides a device for generating a first arrival gather at a common center point, including:
the rotation module is used for rotating the first coordinates of each detection point and each shot point according to the azimuth angle of the measuring line to obtain the second coordinates in a new coordinate system taking the measuring line as the positive east or the positive north;
the grid division module is used for carrying out grid division on the new coordinate system according to the preset bin size, and determining the bin where the common center point of each detector-shot point pair is located according to the second coordinates of each detector and each shot point;
The seismic channel position sequence determining module is used for determining the common center point of the first arrival seismic channel of each shot and the position sequence of the first arrival seismic channel in all the first arrival seismic channels sharing the common center point;
the total channel number determining module is used for counting the total channel number of the first arrival seismic channels contained in the common center point of each grid surface element according to the identification of the surface element where the common center point of the first arrival seismic channel of each gun is located;
the initial position determining module is used for distributing the recording position of the first arrival seismic trace of the common center point of each surface element in the first arrival seismic trace of the common center point of each surface element as the initial recording position of the common center point according to the total trace number of the first arrival seismic trace contained in the common center point of each surface element and the sequence of each surface element in the grid;
the system comprises a common-center first-arrival gather file generation module, a first-arrival gather file generation module and a first-arrival gather file generation module, wherein the common-center first-arrival gather file generation module is used for determining the actual recording position of each first-arrival gather in the common-center first-arrival gather file according to the initial recording position of the common-center where the first-arrival gather is located and the position sequence of the first-arrival gather in the common-center, and writing the related information of the first-arrival gather into the actual recording position to generate the common-center first-arrival gather file.
In a fourth aspect, an embodiment of the present invention provides a device for calculating a residual static correction of a first arrival wave, where a common-center first arrival gather used by the device for calculating the residual static correction of the first arrival wave is obtained by using a method for generating the common-center first arrival gather as described above.
In a fifth aspect, an embodiment of the present invention provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, where the processor implements a method for generating a common-center first-arrival gather as described above or a method for calculating a first-arrival-wave residual static correction as described above when executing the computer program.
In a sixth aspect, embodiments of the present invention provide a computer-readable storage medium storing a computer program for executing the method of generating a common-center first-arrival gather as described above, or a computer program for implementing the method of calculating a first-arrival-wave residual static correction as described above.
The technical scheme provided by the invention has the beneficial effects that at least:
in the method and the device for generating the first arrival channel set file of the common center point, the positions of all the first arrival channels of each surface element are meshed, the sequence of each first arrival channel in the common center point is determined, the number of the first arrival channels of the surface element where each common center point is located and the storage position of the first arrival channel of each surface element are counted according to the sequence, so that the storage position of the first arrival channel of the surface element and the sequence of the first arrival channel can be utilized when the first arrival channel set file of the common center point is generated, and the offset of the addresses of all the first arrival channels of each surface element can be rapidly determined, and the writing of the data of all the first arrival channels of each common center point corresponding to each surface element is rapidly completed, so that the first arrival channel set file of the common center point is generated. The method shortens the processing time of the first-arrival seismic channel data, improves the generation efficiency of the first-arrival gather of the common center point, and better meets the aging requirement of seismic data processing.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims thereof as well as the appended drawings.
The technical scheme of the invention is further described in detail through the drawings and the embodiments.
Drawings
The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate the invention and together with the embodiments of the invention, serve to explain the invention. In the drawings:
FIG. 1 is a flowchart of a method for generating a common-point first-arrival gather according to an embodiment of the present invention;
FIGS. 2A-2C are schematic diagrams of an original coordinate system and a new coordinate system after rotation in different rotation modes according to an embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating meshing of a new rotated coordinate system in an embodiment of the present invention;
FIG. 4 is a schematic illustration of common center point locations of different offset pairs within bins in an embodiment of the invention;
FIGS. 5A and 5B are schematic diagrams of bin line numbers and point numbers representing common center points in an embodiment of the invention;
FIG. 6 is a diagram of the total trace number of first-arrival seismic traces contained in each common center point in an embodiment of the invention;
FIG. 7 is a schematic diagram of storage locations of first-arrival seismic traces contained in respective centerpoints in accordance with an embodiment of the present invention;
FIG. 8 is a bar graph comparing prior art and example embodiments of the present invention when used;
fig. 9 is a block diagram of a device for generating a first-arrival gather at a common center point according to an embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below 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 embodied 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.
In order to shorten the generation time of the first-arrival gather file at the common center and improve the processing efficiency of the data of the first-arrival gather file at the common center, the embodiment of the invention provides a method for generating the first-arrival gather at the common center, which is shown by referring to fig. 1, and comprises the following steps:
S11, according to azimuth angles of the measuring lines, rotating first coordinates of each detection point and each shot point to obtain second coordinates in a new coordinate system taking the measuring lines as the positive east or the positive north;
s12, performing grid division on the new coordinate system according to the preset bin size, and determining bins of common center points of all detection point-shot point pairs according to second coordinates of all detection points and all shot points;
s13, determining a common center point of a first-arrival seismic channel of each gun and the position sequence of the first-arrival seismic channel in all the first-arrival seismic channels sharing the common center point;
s14, counting the total channel number of the first arrival seismic channels contained in the common center point of each grid surface element according to the identification of the surface element where the common center point of the first arrival seismic channel of each gun is located;
s15, according to the total channel number of the first-arrival seismic channels contained in the common center point of each surface element, distributing the recording position of the first-arrival seismic channel in the common center point first-arrival gather file as the initial recording position of the common center point according to the sequence of each surface element in the grid;
s16, determining the actual recording position of each first-arrival seismic channel in the common-center-point first-arrival gather file according to the initial recording position of the common center point of the first-arrival seismic channel and the position sequence of the first-arrival seismic channel in the common center point, and writing the related information of the first-arrival seismic channel into the actual recording position to generate the common-center-point first-arrival gather file.
According to the method for generating the first arrival channel set file of the common center point, provided by the embodiment of the invention, the positions of the common center points are meshed, the sequence of each first arrival channel in the common center point is determined, the number of first arrival channels of the surface element where each common center point is located and the storage position of the first arrival channel of each surface element are counted according to the sequence, so that when the first arrival channel set file of the common center point is generated, the storage position of the first arrival channel of the surface element and the sequence of the first arrival channel can be utilized, the offset of the addresses of all the first arrival channels of each surface element can be rapidly determined, and therefore the writing of the data of all the first arrival channels of the common center points corresponding to each surface element can be rapidly completed, and the first arrival channel set file of the common center point is generated. The method shortens the processing time of the first-arrival seismic channel data, improves the generation efficiency of the first-arrival gather of the common center point, and better meets the aging requirement of seismic data processing.
Before executing the steps S11 to S16, the embodiment of the present invention further arranges the SPS files in order, establishes indexes for the detectors and shots, picks up the first arrival time, and generates data of the first arrival seismic trace.
Indexing the spots and shots, i.e. arranging the spot files, shot files in ascending order, and indexing each spot and shot according to its sequential position in all spots and shots, e.g. numbering the shots: shot 1, shot 2, shot 3 … …, and the detector number is: pick-up point 1, pick-up point 2, pick-up point 3 … ….
The first arrival time is picked up, data of first arrival seismic channels are generated, the first arrival time is picked up in a manual interaction mode or an automatic pickup mode according to the seismic data collected in the field, the quality and the quantity of first arrival pickup are guaranteed, and the data of the first arrival seismic channels are generated.
Further, for convenience of explanation, the original coordinates of each detector point and each shot point in a field actual coordinate system are referred to as "first coordinates", and the coordinates obtained after rotation are referred to as "second coordinates".
In the step S11, the first coordinates of each detector and each shot are rotated according to the azimuth of the measuring line, so as to obtain the second coordinates in the new coordinate system with the measuring line being the positive east or the positive north. The following are respectively described:
mode one:
rotating the first coordinate in a clockwise direction;
the second coordinates are obtained as follows:
X R ＝(X-X 0 )*cos(90-θ)+(Y-Y 0 )*sin(90-θ)+X 0
Y R ＝(Y-Y 0 )*cos(90-θ)-(X-X 0 )*sin(90-θ)+Y 0
wherein:
X R the east coordinate in the second coordinate obtained after rotation;
Y R for north in the second coordinates obtained after rotationMarking;
x is the east coordinate in the first coordinate before rotation;
y is the north coordinate of the first coordinate before rotation;
X 0 An east coordinate for the specified origin;
Y 0 north coordinates for a specified origin;
θ is the azimuth of the line.
Mode two:
rotating the first coordinate in a counterclockwise direction;
the second coordinates obtained are:
X R ＝(X-X 0 )*cos(θ)-(Y-Y 0 )*sin(θ)+X 0
Y R ＝(Y-Y 0 )*cos(θ)+(X-X 0 )*sin(θ)+Y 0
wherein:
X R the east coordinate in the second coordinate obtained after rotation;
Y R the north coordinate in the second coordinate obtained after rotation;
x is the east coordinate in the first coordinate before rotation;
y is the north coordinate of the first coordinate before rotation;
X 0 an east coordinate for the specified origin;
Y 0 north coordinates for a specified origin;
θ is the azimuth of the line.
Referring to fig. 2A to 2C, fig. 2A is a schematic view of an original coordinate system before rotation, fig. 2B is a schematic view of a coordinate system after clockwise rotation, and fig. 2C is a schematic view of a coordinate system after counterclockwise rotation. Fig. 2A to 2C each include a large number of shots and detectors (different areas are shown to be different).
Further, in the step S12, the new coordinate system is gridded according to the preset bin size, specifically, the new coordinate system where the second coordinates of each detection point and each shot point are located is gridded according to the bin size preset by the observation system, and the minimum unit of the grid obtained after the division is a bin; the method comprises the steps of marking each bin in a line number-point number mode sequentially from the bin with the smallest abscissa to the bin with the smallest ordinate according to the descending order of the abscissa.
The purpose of the aforementioned step S11 is to enable the rotated new coordinate system to more quickly locate the position of the bin where the common center point is located, that is, to quickly locate one common center point (each common center point will only fall into one bin, that is, only correspond to one bin) by calculating the bin position and the corresponding number thereof only through the east coordinate and the north coordinate. Thereby improving the calculation speed and calculation efficiency.
In order to quickly locate the common center point of the bin, the size of the bin needs to be determined by fully considering the distribution density condition of the common center point formed by each detector-shot point pair during grid division.
For example, the new coordinate system after rotation may be meshed according to a specification of a bin size (10 meters in a line direction and 20 meters in a Crossline direction perpendicular to the line direction), and the meshed mesh is shown in fig. 3. The surface element is a square grid (rectangle or square), and the specific size is not limited in the embodiment of the invention.
Referring to FIG. 4, the shared center point (common center point) of offset pairs from different directions and different offset magnitudes is located within a bin.
After meshing, each bin is numbered, each bin representing a common center point, and thus the bins, i.e., the common center points in the bins, are numbered. The number of the bin in which each common center point is located can be determined as follows:
If the direction of the measuring line in the new coordinate system after rotation is the north direction, the line number and the point number of the bin where the common center point of the first arrival seismic channel is located are respectively as follows:
i＝Round((RX R +SX R )/2-cmpXzbmin)/ΔIz)+1；
j＝Round((RY R +SY R )/2-cmpYzbmin)/ΔCz)+1；
wherein:
i is the line number of the bin where the calculated common center point is located, and the value range is from 1 to M;
j is the point number of the calculated common-center point face element, and the value range is from 1 to N;
round is a rounding function;
RX R the east coordinate in the second coordinate is obtained after the wave detection point rotates;
SX R the east coordinates in the second coordinates are obtained after the shot point rotates;
RY R the north coordinate in the second coordinate obtained after the rotation of the wave detector;
SY R the north coordinates in the second coordinates obtained after the shot point rotates;
cmpXzbmin is the line number-the point number is the east coordinate of the (1, 1) bin center;
cmpYzbmin is the line number-the point number is the north coordinate of the (1, 1) bin center;
ΔIz is the length of one bin of the line direction (i.e., the line direction) of the line in the new coordinate system after rotation;
Δcz is the length of one bin in the Crossline direction perpendicular to the Inline direction in the new coordinate system after rotation.
(II) if the direction of the measuring line in the new coordinate system after rotation is the forward direction, the line number and the point number of the bin where the common center point of the first arrival seismic channel is located are respectively as follows:
i＝Round((RY R +SY R )/2-cmpYzbmin)/ΔIz)+1；
j＝Round((RX R +SX R )/2-cmpXzbmin)/ΔCz)+1；
wherein:
i is the line number of the bin where the calculated common center point is located, and the value range is from 1 to M;
j is the point number of the calculated common-center point face element, and the value range is from 1 to N;
round is a rounding function;
RX R the east coordinate in the second coordinate is obtained after the wave detection point rotates;
SX R the east coordinates in the second coordinates are obtained after the shot point rotates;
RY R the north coordinate in the second coordinate obtained after the rotation of the wave detector;
SY R the north coordinates in the second coordinates obtained after the shot point rotates;
cmpXzbmin is the line number-the point number is the east coordinate of the (1, 1) bin center;
cmpYzbmin is the line number-the point number is the north coordinate of the (1, 1) bin center;
ΔIz is the length of one bin of the line in the line direction in the new coordinate system after rotation;
Δcz is the length of one bin in the Crossline direction perpendicular to the Inline direction in the new coordinate system after rotation.
Referring to fig. 5A, the number of each bin (representing one common center point) is composed of a line number and a point number, and the dotted line in fig. 5A indicates the order of the line numbers and the solid line indicates the order of the point numbers. The same line number has a plurality of different points to represent different surface elements, and the line number of the different surface elements is formed at the intersection of the solid line of the same point number and the dotted line represented by the different line numbers. The range of line numbers is: 1 to M, the range of point numbers is: 1 to N.
The bin where the common center point is located is numbered (line number, point number) from the bin with the smallest abscissa in the grid. The line number starts from 1, the increment is 1, M is the maximum line number of the upper element of the grid; the point numbers on each line are numbered starting from 1, the increment is 1, and N is the maximum point number on each line on the grid.
Of course, in the implementation of the embodiment of the present invention, the line number and the point number are not necessarily all numbers from 1, referring to fig. 5B, each bin number is composed of 7 digits, and numbers are sequentially numbered from 1001001, where the number of the leftmost column of bins is 1001001 ～ 1001010, the first four digits are the dot line number, and the last three digits are the point number.
After numbering is completed, statistics can be performed on the first arrival seismic traces of the common center point of each bin, and the total trace number of the first arrival seismic traces contained in each common center point (namely in the bin) can be determined. An example of the total number of seismic traces contained in each common center point is shown with reference to fig. 6.
Aiming at the seismic channels with first arrivals in each gun, calculating the line numbers and the point numbers of grids of the common center points of the first arrival seismic channels in the grids of the common center points channel by channel, and simultaneously accumulating and recording the total channel number of the first arrival seismic channels of the common center points.
Referring to the example shown in FIG. 4, the common center point contains a number of first-arrival seismic traces of 5 (assuming that there is one first-arrival seismic trace data per shot).
The common-center bin is calculated first in bin order, in which the order of each offset pair in the bin is also recorded, as numbered in fig. 4, with bin 1 and bin 1 representing the first bin information for that offset pair, bin 2 and bin 2 representing the second bin information for that offset pair, and so on. After the bin is assigned to the starting position in the stored file, the stored position of the offset pair in the file can be quickly calculated according to the starting position, which is also an important factor for improving the operation speed in the embodiment of the invention.
The initial position is the recording position of the first arrival seismic trace of the common center point of the bin (the first arrival seismic trace of the smallest offset pair when the shot numbers are ordered in ascending order) in the common center point first arrival gather file.
Further, in the above step S15, the recording position of the first arrival seismic trace of each common center point in the common center point first arrival gather file is calculated and allocated as follows:
T 1,1 ＝1；
wherein:
T 1,1 representing the recording position of a first arrival seismic trace of a common center point of a line number-point number (1, 1) bin in the common center point first arrival gather file;
T i,j representing the recording position of a first arrival seismic trace of a common center point of a line number-point number (i, j) bin in the common center point first arrival gather file;
C l,k representing the total number of first-arrival seismic traces with line numbers-points numbered as the common center point of the (l, k) bin;
i and l represent line numbers of the cells where the common center point is located, i takes values from 1 to M, l takes values from 1 to i-1, and M represents the maximum line number of the cells in the grid.
j and k represent the point numbers of the cells where the common center point is located, j takes values from 1 to N, k takes values from 1 to j-1, and N represents the maximum point number of the cells in the grid.
Based on the initial position, all the first-arrival seismic trace data of the common center point can be quickly written into the corresponding position of the first-arrival gather file according to the offset.
Specifically, the actual recording position in S16 can be calculated by the following formula:
R i,j ＝T i,j +N s,i,j -1；
wherein: r is R i,j Representing the actual recording position of the related information of the first-arrival seismic channel in the first-arrival gather file of the common center point;
T i,j A start recording position representing a line number-point number of a common center point of the (i, j) bin;
N s,i,j representing the order of locations in the common center point (i, j) of the first-arrival seismic trace on-line number-point number (i, j) bin;
s is the shot number.
Referring to fig. 7, according to statistics of the number of first-arrival seismic traces of each bin, storage locations (addresses) in a first-arrival gather file of a bin (common-center point) are then sequentially allocated in order of size of the bin line number-point number. For example, 1-5 in FIG. 7 indicates that the storage locations of the first-arrival seismic trace data for the 1001001 bin in FIG. 6 are 1-5 in the common-center first-arrival gather file, and so on, 137-148 indicates that the storage locations of the first-arrival seismic trace data for the 1002005 bin are 137-148 in the common-center first-arrival gather file.
The first arrival seismic trace data includes: index of the detection point of the first-arrival seismic channel, index of the shot point of the first-arrival seismic channel, offset of the first-arrival seismic channel, first-arrival time of the first-arrival seismic channel and the like.
The method for generating the first arrival gather of the common center point is briefly described below by using an example, and the flow of the example includes the following steps:
1) Sequentially arranging SPS and establishing an index;
the indexes of the first wave detection point and the gun point are respectively 1, and the indexes of the following wave detection point and gun point are respectively accumulated according to increment 1.
2) Picking first arrival time
And carrying out pickup operation of the first arrival time according to the seismic data collected in the field to obtain data of the first arrival seismic channel.
3) Rotation of coordinates
The coordinate rotation in the step 3) is to rotate all the detection points and the shot coordinates according to the appointed origin of coordinates and the clockwise rotation direction, wherein the azimuth angle of the measuring line is 37.50, so as to obtain new coordinates of a coordinate system which takes the positive east direction as the horizontal coordinate and the positive north direction as the vertical coordinate after all the detection points and the shot rotate, and the direction of the measuring line after rotation is the positive north direction;
the coordinate clockwise rotation formula is:
X R ＝(X-X 0 )*cos(90-37.5)+(Y-Y 0 )*sin(90-37.5)+X 0
Y R ＝(Y-Y 0 )*cos(90-37.5)-(X-X 0 )*sin(90-37.5)+Y 0
wherein: x is X R Is the rotated east coordinate;
Y R is the rotated north coordinates;
x is the east coordinate before rotation;
y is the north coordinate before rotation;
X 0 an east coordinate for the specified origin;
Y 0 north coordinates for a specified origin;
4) Common-center point bin layout
The arrangement of the common center point in the step 4) is to divide the common center point into grids according to the new coordinates obtained in the step 3) according to the preset size of the surface element of the observation system (10 meters in the line direction and 20 meters in the cross line direction), and start with the surface element of the minimum horizontal and vertical coordinates in the grids, and number the surface element corresponding to the common center point according to the line number and the point number. The line number starts from 1, the increment is 1, 650 is the maximum line number of the concentric point grid; the point numbers on each line are numbered starting from 1 with an increment of 1, 1200 being the maximum point number on each line of the concentric point grid.
5) Co-center first-arrival data initialization
The initial value of the number of the first arrival channels of the earthquake contained in the surface element where each common center point is located is 0.
6) Co-center point bin information calculation
When the direction of the rotated measuring line is in the north direction, the calculating formulas of the line number and the point number of the common center point of the first arrival seismic channel are respectively as follows:
i＝Round((RX R +SX R )/2-cmpXzbmin)/10)+1
j＝Round((RY R +SY R )/2-cmpYzbmin)/20)+1
wherein: i is the line number of the grid of the calculation common center points, and the value range is from 1 to 650;
j is the point number of the calculated common-center point grid, and the value range is from 1 to 1200;
round is a rounding function;
RX R the east coordinate of the wave detector after rotation;
SX R the east coordinates of the shot after rotation;
RY R the north coordinate of the wave detector after rotation;
SY R the north coordinates of the shot after rotation;
cmpXzbmin is the line number-the point number is the east coordinate of the (1, 1) bin center;
cmpYzbmin is the line number-the point number is the north coordinate of the (1, 1) bin center;
7) Allocation of co-centered point bin addresses
In the face information calculation result obtained in the step 6), the initial record position, which is to be distributed in the first arrival file of the common center point, of the first arrival seismic channel data of each common center point is calculated in a statistics mode. Starting from the surface element (1, 1), calculating the initial position of the first arrival data of the grid in the first arrival file of the common center point to be generated according to the sequence of line numbers and point numbers for all grids, wherein the calculation mode is as follows:
T 1,1 ＝1
wherein: t (T) 1,1 Representing the recording position of a first arrival seismic trace of a common center point of a line number-point number (1, 1) bin in the common center point first arrival gather file;
T i,j representing the recording position of a first arrival seismic trace of a common center point of a line number-point number (i, j) bin in the common center point first arrival gather file;
C l,k representing the total number of first-arrival seismic traces with line numbers-points numbered as the common center point of the (l, k) bin;
i and l represent line numbers of the cells where the common center point is located, i takes on values from 1 to M, and l takes on values from 1 to i-1. Where m=650.
j and k represent the point numbers of the cells where the common center point is located, j takes the value from 1 to N, and k takes the value from 1 to j-1. Where n=1200.
8) Generating a common-center first-arrival file
And starting from the first shot, writing the related information of the seismic channels one by one into a common-center first arrival file for the seismic channels with the first arrival. The related information includes: the indexes obtained by the detection point of the seismic channel in the step 1), the indexes obtained by the shot point of the seismic channel in the step 1), the offset of the seismic channel and the first arrival time of the seismic channel.
9) Co-centrality first-arrival file application
And the generated common-center point first-arrival data file is used for calculating the residual static correction of the first-arrival wave or for monitoring the static correction application effect.
FIG. 8 is a graph comparing the time duration of a method according to the prior art and the method according to the embodiment of the present invention, which are run on the same microcomputer, for generating a three-dimensional common-center first-arrival seismic gather file in a target area. As can be seen from comparison of FIG. 8, the method in the prior art only takes 3 hours, while the method in the embodiment of the invention only takes 90 hours, so that the operation speed is improved by 30 times after the method provided by the embodiment of the invention is adopted, and the time requirement of three-dimensional treatment can be well met.
The embodiment of the invention also provides a method for calculating the residual static correction of the first arrival wave, which is obtained by adopting the method for generating the common-center first arrival gather, wherein the common-center first arrival gather is used in the calculation of the residual static correction of the first arrival wave.
The specific implementation process of calculating the residual static correction of the first arrival wave can refer to the prior art, and the embodiment of the invention is not limited to this.
Based on the same inventive concept, the embodiment of the invention also provides a device for generating the first arrival gather of the common center point and a device for calculating the residual static correction of the first arrival wave, and because the principle of solving the problems by the devices is similar to the method for generating the first arrival gather of the common center point and the method for calculating the residual static correction of the first arrival wave, the implementation of the device can refer to the implementation of the method, and the repetition is omitted.
The device for generating the first arrival gather at the common center point according to the embodiment of the present invention, as shown in fig. 9, includes:
the rotation module 91 is configured to rotate the first coordinates of each detection point and each shot point according to the azimuth of the measurement line, so as to obtain a second coordinate in a new coordinate system with the measurement line being the positive east or the positive north;
the meshing module 92 is configured to mesh the new coordinate system according to a preset bin size, and determine a bin where a common center point of each detector-shot pair is located according to second coordinates of each detector and each shot;
the seismic channel position sequence determining module 93 is used for determining a common center point of the first arrival seismic channel of each shot and the position sequence of the first arrival seismic channel in all the first arrival seismic channels sharing the common center point;
the total channel number determining module 94 is configured to count the total channel number of the first arrival seismic channels included in the common center point of each grid bin according to the identifier of the bin where the common center point of the first arrival seismic channel of each gun is located;
the initial position determining module 95 is configured to allocate, according to the total number of first-arrival seismic traces included in the common-center point of each bin, a recording position of the first-arrival seismic trace in the common-center point first-arrival gather file of each bin as a starting recording position of the common-center point according to the sequence of each bin in the grid;
The common-center first-arrival gather file generation module 96 is configured to determine, for each first-arrival seismic trace, an actual recording position of each first-arrival seismic trace in the common-center first-arrival gather file according to the initial recording position of the common-center where the first-arrival seismic trace is located and the position sequence of the first-arrival seismic trace in the common-center, and write related information of the first-arrival seismic trace into the actual recording position to generate the common-center first-arrival gather file.
The embodiment of the invention also provides a device for calculating the residual static correction of the first arrival wave, wherein the common-center first arrival gather used by the device for calculating the residual static correction of the first arrival wave is obtained by adopting the generation method of the common-center first arrival gather.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the generation method of the first arrival gather of the common center point or the method for calculating the residual static correction of the first arrival wave when executing the computer program.
The embodiment of the invention also provides a computer readable storage medium storing a computer program for executing the method for generating the common-center first-arrival gather or a computer program for realizing the method for calculating the residual static correction of the first-arrival wave.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, magnetic disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (14)
1. A method for generating a common-center first-arrival gather, comprising:
according to the azimuth angle of the measuring line, rotating the first coordinates of each detecting point and each shot point to obtain the second coordinates in a new coordinate system taking the measuring line as the positive east or the positive north;
grid division is carried out on the new coordinate system according to the preset bin size, and bins where common center points of all detection points-shot point pairs are located are determined according to second coordinates of all detection points and all shot points;
determining a common center point of a first arrival seismic channel of each shot and the position sequence of the first arrival seismic channel in all the first arrival seismic channels sharing the common center point;
counting the total channel number of the first arrival seismic channels contained in the common center point of each grid surface element according to the identification of the surface element where the common center point of each first arrival seismic channel is located;
according to the total channel number of the first-arrival seismic channels contained in the common-center point of each surface element, distributing the recording position of the first-arrival seismic channel in the common-center point first-arrival channel set file as the initial recording position of the common-center point according to the sequence of each surface element in the grid;
for each first-arrival seismic channel, determining the actual recording position of each first-arrival seismic channel in a first-arrival gather file of the common center point according to the initial recording position of the common center point where the first-arrival seismic channel is located and the position sequence of the first-arrival seismic channel in the common center point, and writing the related information of the first-arrival seismic channel into the actual recording position to generate the first-arrival gather file of the common center point.
2. The method of claim 1, wherein the rotating the first coordinates of each detector point and each shot point according to the azimuth of the line to obtain the second coordinates in the new coordinate system with the line being the positive east or the positive north, specifically comprises:
rotating the first coordinate in a clockwise direction;
the second coordinates are obtained as follows:
X R ＝(X-X 0 )*cos(90-θ)+(Y-Y 0 )*sin(90-θ)+X 0
Y R ＝(Y-Y 0 )*cos(90-θ)-(X-X 0 )*sin(90-θ)+Y 0
wherein:
X R the east coordinate in the second coordinate obtained after rotation;
Y R the north coordinate in the second coordinate obtained after rotation;
x is the east coordinate in the first coordinate before rotation;
y is the north coordinate of the first coordinate before rotation;
X 0 an east coordinate for the specified origin;
Y 0 north coordinates for a specified origin;
θ is the azimuth of the line.
3. The method of claim 1, wherein the rotating the first coordinates of each detector point and each shot point according to the azimuth of the line to obtain the second coordinates in the new coordinate system with the line being the positive east or the positive north, specifically comprises:
rotating the first coordinate in a counterclockwise direction;
the second coordinates obtained are:
X R ＝(X-X 0 )*cos(θ)-(Y-Y 0 )*sin(θ)+X 0
Y R ＝(Y-Y 0 )*cos(θ)+(X-X 0 )*sin(θ)+Y 0
wherein:
X R the east coordinate in the second coordinate obtained after rotation;
Y R the north coordinate in the second coordinate obtained after rotation;
X is the east coordinate in the first coordinate before rotation;
y is the north coordinate of the first coordinate before rotation;
X 0 an east coordinate for the specified origin;
Y 0 north coordinates for a specified origin;
θ is the azimuth of the line.
4. The method of claim 1, wherein meshing the new coordinate system according to a preset bin size, specifically comprises:
according to the size of a bin preset by an observation system, performing grid division on a new coordinate system where the second coordinates of each detection point and each shot point are located, wherein the minimum unit of the grid obtained after division is a bin;
the method comprises the steps of marking each bin in a line number-point number mode sequentially from the bin with the smallest abscissa to the bin with the smallest ordinate according to the descending order of the abscissa.
5. The method of claim 4, wherein the bin in which the common center point is located is determined by:
if the direction of the rotated measuring line is the north direction, the line number and the point number of the bin where the common center point of the first arrival seismic channel is located are respectively:
i＝Round((RX R +SX R )/2-cmpXzbmin)/ΔIz)+1；
j＝Round((RY R +SY R )/2-cmpYzbmin)/ΔCz)+1；
wherein:
i is the line number of the bin where the calculated common center point is located, and the value range is from 1 to M;
j is the point number of the calculated common-center point face element, and the value range is from 1 to N;
Round is a rounding function;
RX R the east coordinate in the second coordinate is obtained after the wave detection point rotates;
SX R the east coordinates in the second coordinates are obtained after the shot point rotates;
RY R the north coordinate in the second coordinate obtained after the rotation of the wave detector;
SY R the north coordinates in the second coordinates obtained after the shot point rotates;
cmpXzbmin is the line number-the point number is the east coordinate of the (1, 1) bin center;
cmpYzbmin is the line number-the point number is the north coordinate of the (1, 1) bin center;
ΔIz is the length of one bin of the line in the line direction in the new coordinate system after rotation;
Δcz is the length of one bin in the Crossline direction perpendicular to the Inline direction in the new coordinate system after rotation.
6. The method of claim 4, wherein the identity of the bin at which the common center point of the first-arrival seismic trace for each shot is located is determined by:
if the direction of the rotated survey line is the forward direction, the line number and the point number of the bin where the common center point of the first arrival seismic channel is located are respectively:
i＝Round((RY R +SY R )/2-cmpYzbmin)/ΔIz)+1；
j＝Round((RX R +SX R )/2-cmpXzbmin)/ΔCz)+1；
wherein:
i is the line number of the bin where the calculated common center point is located, and the value range is from 1 to M;
j is the point number of the calculated common-center point face element, and the value range is from 1 to N;
round is a rounding function;
RX R The east coordinate in the second coordinate is obtained after the wave detection point rotates;
SX R the east coordinates in the second coordinates are obtained after the shot point rotates;
RY R the north coordinate in the second coordinate obtained after the rotation of the wave detector;
SY R the north coordinates in the second coordinates obtained after the shot point rotates;
cmpXzbmin is the line number-the point number is the east coordinate of the (1, 1) bin center;
cmpYzbmin is the line number-the point number is the north coordinate of the (1, 1) bin center;
ΔIz is the length of one bin of the line in the line direction in the new coordinate system after rotation;
Δcz is the length of one bin in the Crossline direction perpendicular to the Inline direction in the new coordinate system after rotation.
7. The method of claim 4, wherein assigning the first arrival seismic trace at the common center point for each bin its recorded location in the common center point first arrival gather file comprises:
the recording position of the first arrival seismic channel of each common center point in the common center point first arrival gather file is calculated according to the following mode:
T 1,1 ＝1；
wherein:
T 1,1 representing the recording position of a first arrival seismic trace of a common center point of a line number-point number (1, 1) bin in the common center point first arrival gather file;
T i,j Representing the recording position of a first arrival seismic trace of a common center point of a line number-point number (i, j) bin in the common center point first arrival gather file;
C l,k representing the total number of first-arrival seismic traces with line numbers-points numbered as the common center point of the (l, k) bin;
i and l represent line numbers of the cells where the common center point is located, i takes values from 1 to M, l takes values from 1 to i-1, and M represents the maximum line number of the cells in the grid.
j and k represent the point numbers of the cells where the common center point is located, j takes values from 1 to N, k takes values from 1 to j-1, and N represents the maximum point number of the cells in the grid.
8. The method of claim 4, wherein determining the actual recording position of each first-arrival seismic trace in the common-point first-arrival gather file according to the starting recording position of the common-point where the first-arrival seismic trace is located and the position order of the first-arrival seismic trace in the common-point, specifically comprises:
the actual recording position is calculated by the following formula:
R i,j ＝T i,j +N s,i,j -1；
wherein: r is R i,j Representing the actual recording position of the related information of the first-arrival seismic channel in the first-arrival gather file of the common center point;
T i,j a start recording position representing a line number-point number of a common center point of the (i, j) bin;
N s,i,j Representing the order of locations in the common center point (i, j) of the first-arrival seismic trace on-line number-point number (i, j) bin;
s is the shot number.
9. The method of any of claims 1-8, wherein prior to rotating the first coordinates of each detector point and each shot point based on the azimuth of the survey line, the method further comprises:
according to the sequence of the shots and the wave detectors, respectively establishing indexes for each shot and each wave detector;
and carrying out pickup operation of the first arrival time according to the seismic data collected in the field to obtain data of the first arrival seismic channel.
10. A method for calculating the residual static correction of a first arrival wave, wherein the common-point first arrival 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-point first arrival gather according to any one of claims 1 to 9.
11. A device for generating a set of first arrival tracks at a common center point, comprising:
the rotation module is used for rotating the first coordinates of each detection point and each shot point according to the azimuth angle of the measuring line to obtain the second coordinates in a new coordinate system taking the measuring line as the positive east or the positive north;
the grid division module is used for carrying out grid division on the new coordinate system according to the preset bin size, and determining the bin where the common center point of each detector-shot point pair is located according to the second coordinates of each detector and each shot point;
The seismic channel position sequence determining module is used for determining the common center point of the first arrival seismic channel of each shot and the position sequence of the first arrival seismic channel in all the first arrival seismic channels sharing the common center point;
the total channel number determining module is used for counting the total channel number of the first arrival seismic channels contained in the common center point of each grid surface element according to the identification of the surface element where the common center point of the first arrival seismic channel of each gun is located;
the initial position determining module is used for distributing the recording position of the first arrival seismic trace of the common center point of each surface element in the first arrival seismic trace of the common center point of each surface element as the initial recording position of the common center point according to the total trace number of the first arrival seismic trace contained in the common center point of each surface element and the sequence of each surface element in the grid;
the system comprises a common-center first-arrival gather file generation module, a first-arrival gather file generation module and a first-arrival gather file generation module, wherein the common-center first-arrival gather file generation module is used for determining the actual recording position of each first-arrival gather in the common-center first-arrival gather file according to the initial recording position of the common-center where the first-arrival gather is located and the position sequence of the first-arrival gather in the common-center, and writing the related information of the first-arrival gather into the actual recording position to generate the common-center first-arrival gather file.
12. A device for calculating residual static correction of first-arrival waves, wherein the common-point first-arrival gather used by the device for calculating residual static correction of first-arrival waves is obtained by adopting the method for generating the common-point first-arrival 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 implements a method of generating a set of first-arrival-point first-arrival-gathers as claimed in any one of claims 1 to 9 or a method of calculating a first-arrival-wave residual static correction as claimed in claim 10 when executing the computer program.
14. A computer-readable storage medium storing a computer program for executing the method of generating a common-center first-arrival gather according to any one of claims 1 to 9, or a computer program for implementing the method of calculating a first-arrival-wave residual static correction according to claim 10.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111651249.5A CN116413774A (en) | 2021-12-30 | 2021-12-30 | Common-center-point first-arrival gather generation method, calculation residual static correction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111651249.5A CN116413774A (en) | 2021-12-30 | 2021-12-30 | Common-center-point first-arrival gather generation method, calculation residual static correction method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116413774A true CN116413774A (en) | 2023-07-11 |
Family
ID=87054950
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111651249.5A Pending CN116413774A (en) | 2021-12-30 | 2021-12-30 | Common-center-point first-arrival gather generation method, calculation residual static correction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116413774A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117519610A (en) * | 2024-01-05 | 2024-02-06 | 中国石油集团东方地球物理勘探有限责任公司 | Method, device and medium for creating CMP data body by post-stack seismic data |
-
2021
- 2021-12-30 CN CN202111651249.5A patent/CN116413774A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117519610A (en) * | 2024-01-05 | 2024-02-06 | 中国石油集团东方地球物理勘探有限责任公司 | Method, device and medium for creating CMP data body by post-stack seismic data |
| CN117519610B (en) * | 2024-01-05 | 2024-03-22 | 中国石油集团东方地球物理勘探有限责任公司 | Method, device and medium for creating CMP data body by post-stack seismic data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103984024B (en) | A kind of three-component geophone horizontal component data auto-correction method | |
| CN116413774A (en) | Common-center-point first-arrival gather generation method, calculation residual static correction method and device | |
| CN105093319B (en) | Ground micro-seismic static correcting method based on 3D seismic data | |
| WO2017107862A1 (en) | Method and device for determining overall measurement of seismic observation system repeatability | |
| CN108828668A (en) | A kind of pre-stack time migration data processing method and device | |
| CN111638551A (en) | Seismic first-motion wave travel time chromatography method and device | |
| CN104199103A (en) | Static correction method and static correction device for converted wave | |
| CN109298450A (en) | Method and device for analyzing attributes of seismic acquisition observation system | |
| CN105319576A (en) | Concave polygon region observation system automatic full coverage quick layout method | |
| CN110411360A (en) | A kind of smart phone surveys tree device and its tree survey method | |
| CN113534259B (en) | Real-time prestack time migration imaging method for efficient collection of controllable seismic source | |
| CN104570098B (en) | A kind of acquisition methods and its device of azimuth geological data altogether | |
| CN104502967A (en) | Method and device for quickly obtaining surface element information of seismic exploration and observation system | |
| CN112764102B (en) | Seismic data acquisition method and device | |
| CN108089229A (en) | A kind of detector arrangement location determining method, device and computer storage media | |
| CN110824559B (en) | Method and device for generating co-coordinate point gather and computer-readable storage medium | |
| CN114002742B (en) | Euler Gaussian beam offset imaging method and device | |
| CN107422375B (en) | The determination method and apparatus of the CRP gather of subsea node | |
| CN110187386A (en) | It is a kind of quickly to identify architectonic DTW seismic volume attributes analysis method automatically | |
| CN106814395B (en) | Method and device for determining background cannon | |
| CN104199104B (en) | 3D seismic data surface wave attenuation method and device | |
| CN104880732B (en) | A kind of construction method and device of cross subset | |
| CN110873899B (en) | Seismic data super gather generation method based on offset distance grouping | |
| CN105259568B (en) | A kind of method and apparatus of determination exploratory area shot point maximum offset | |
| CN112444876A (en) | Seismic data processing method and device |
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 |