CN112099091B - Static correction method and device for seismic exploration point data - Google Patents
Static correction method and device for seismic exploration point data Download PDFInfo
- Publication number
- CN112099091B CN112099091B CN201910520558.5A CN201910520558A CN112099091B CN 112099091 B CN112099091 B CN 112099091B CN 201910520558 A CN201910520558 A CN 201910520558A CN 112099091 B CN112099091 B CN 112099091B
- Authority
- CN
- China
- Prior art keywords
- point data
- seismic exploration
- grids
- seismic
- exploration point
- 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.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 132
- 230000003068 static effect Effects 0.000 title claims abstract description 128
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000004590 computer program Methods 0.000 claims description 13
- 230000005284 excitation Effects 0.000 claims description 13
- 238000012216 screening Methods 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 14
- 238000012545 processing Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000000638 solvent extraction Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
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
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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 provides a static correction method and device for seismic exploration point data, wherein the method comprises the following steps: dividing seismic survey point data into a plurality of grids; determining grids which correspond to the target seismic exploration point data and accord with a search range for each target seismic exploration point data, and obtaining the seismic exploration point data corresponding to the target seismic exploration point data; and carrying out static correction on all the obtained seismic exploration point data. The invention can carry out static correction on the seismic exploration point data, and has high efficiency.
Description
Technical Field
The invention relates to the field of seismic exploration, in particular to a method and a device for static correction of seismic exploration point data.
Background
The seismic exploration needs to arrange a plurality of seismic exploration points, including excitation points and receiving points, and static correction needs to be carried out on the data of the seismic exploration points, and before static correction is carried out, the seismic exploration points corresponding to the target seismic exploration points in a specified search range are firstly determined, namely the seismic exploration points are screened, so that the static correction is carried out on the seismic exploration points in the range. At present, when screening the seismic exploration points, traversing all the seismic exploration points once according to the coordinates of each seismic exploration point and a designated search radius, and judging whether the seismic exploration points are within the search radius or not by calculating the distance between each seismic exploration point and the seismic exploration point, so that a pile of the seismic exploration points within the search radius range is found out to carry out static correction. But screening and judging are carried out on a point-by-point basis, the efficiency is very low, the current seismic exploration steps into a big data era, the billions of scattered points even reach the billions level, the point-by-point screening efficiency is low, the seismic exploration points of the whole seismic working area are searched, the cycle is needed to be repeated for N times, the searching efficiency is very low (N is the total number of excitation points and receiving points), and the static correction efficiency of the seismic exploration point data is very low.
Disclosure of Invention
The embodiment of the invention provides a static correction method for seismic exploration point data, which is used for carrying out static correction on the seismic exploration point data and has high efficiency, and the method comprises the following steps:
dividing seismic survey point data into a plurality of grids;
determining grids which correspond to the target seismic exploration point data and accord with the search range for each target seismic exploration point data, and obtaining the seismic exploration point data which correspond to the target seismic exploration point data and accord with the search range;
and carrying out static correction on all the obtained seismic exploration point data.
The embodiment of the invention provides a seismic exploration point data static correction device which is used for carrying out static correction on seismic exploration point data and has high efficiency, and the device comprises:
the network dividing module is used for dividing the seismic exploration point data into a plurality of grids;
the screening module is used for determining grids which correspond to the target seismic exploration point data and accord with the search range for each target seismic exploration point data, and obtaining the seismic exploration point data which correspond to the target seismic exploration point data and accord with the search range;
and the static correction module is used for carrying out static correction on all the obtained seismic exploration point data.
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 computer program is executed by the processor to realize the seismic exploration point data static correction method.
The embodiment of the invention also provides a computer readable storage medium, which stores a computer program for executing the seismic exploration point data static correction method.
In the embodiment of the invention, firstly, seismic exploration point data are divided into a plurality of grids; then, determining grids which correspond to the target seismic exploration point data and accord with the search range for each target seismic exploration point data, and obtaining the seismic exploration point data which correspond to the target seismic exploration point data and accord with the search range; and finally, carrying out static correction on all the obtained seismic exploration point data. In the screening process, firstly, the grid which corresponds to the target seismic exploration point data and accords with the searching range can be obtained quickly, then all the seismic exploration point data in the grid can be obtained, and whether each seismic exploration point data is in the searching range or not is not needed to be judged, so that the searching times can be greatly reduced, the searching efficiency is improved, and the overall efficiency of static correction of the seismic exploration point data is improved.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art. In the drawings:
FIG. 1 is a flow chart of a method for statics correction of seismic survey point data in an embodiment of the invention;
FIG. 2 is a diagram of a meshing scheme in accordance with an embodiment of the present invention;
FIG. 3 is a diagram of a grid numbering scheme in an embodiment of the present invention;
FIG. 4 is a schematic diagram of a seismic survey point data search in an embodiment of the invention;
FIG. 5 is a detailed flow chart of a method for statics correction of seismic survey point data according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a seismic survey point data statics correction apparatus in accordance with an embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention will be described in further detail with reference to the accompanying drawings. The exemplary embodiments of the present invention and their descriptions herein are for the purpose of explaining the present invention, but are not to be construed as limiting the invention.
During processing of land seismic data, the seismic data is typically corrected to a uniform reference plane, which is typically a horizontal plane. The theory of seismic interpretation assumes that the excitation and reception points are in a horizontal plane and that the formation velocity is uniform. However, in practice, the ground is often uneven, and the depths of the excitation points may be different, so that the wave velocity in the low-speed zone and the wave velocity in the stratum are greatly different, and therefore, the actual measured time-distance curve shape must be affected. To eliminate these effects, terrain corrections, depth of excitation corrections, low-speed zone corrections, etc. are performed on the raw seismic data, which are constant for different seismic interfaces at the same observation point, and are therefore collectively referred to as statics corrections. Generalized static corrections also include phase corrections and corrections for instrument factor effects. Static correction can be used in specific application fields such as speed analysis, superposition imaging, signal-to-noise ratio and resolution analysis of seismic records, geometric forms of various geologic bodies, and the like. The purpose of static correction is to make the excitation point and the receiving point on the same plane or curved surface (reference plane) so that the time-distance curve of the reflected wave has hyperbolic form, and the seismic data after static correction is equivalent to the seismic data acquisition on the reference elevation. A seismic trace corresponds to an excitation point and a reception point, and the amount of statics is the sum of the amounts of statics of the emission point and the reception point.
Before static correction is performed, the corresponding seismic exploration points of the target seismic exploration points in the designated search range need to be determined, namely the seismic exploration points are screened, so that the static correction is performed on the seismic exploration points in the range. At present, when screening the seismic exploration points, traversing all the seismic exploration points once according to each seismic exploration point coordinate and a designated search radius, and judging whether the seismic exploration points are within the search radius by calculating the distance between each seismic exploration point and the seismic exploration point, so as to find a pile of the seismic exploration points within the search radius range to carry out static correction, but screening and judging point by point, and having low efficiency, therefore, the embodiment of the invention provides a seismic exploration point data static correction method which is used for carrying out static correction on the seismic exploration point data, and before carrying out static correction, the seismic exploration point data subjected to static correction is quickly obtained, so that the static correction efficiency is improved.
FIG. 1 is a flowchart of a method for statics correction of seismic survey point data according to an embodiment of the invention, as shown in FIG. 1, comprising:
step 101, dividing seismic exploration point data into a plurality of grids;
102, determining grids which correspond to the target seismic exploration point data and accord with the search range for each target seismic exploration point data, and obtaining the seismic exploration point data which correspond to the target seismic exploration point data and accord with the search range;
and 103, performing static correction on all the obtained seismic exploration point data.
In the embodiment of the invention, firstly, seismic exploration point data are divided into a plurality of grids; then, determining a grid which corresponds to the target seismic exploration point data and accords with the search range for each target seismic exploration point data, obtaining the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the search range, and finally, carrying out static correction on all the obtained seismic exploration point data. In the screening process, firstly, the grid which corresponds to the target seismic exploration point data and accords with the searching range can be obtained quickly, then all the seismic exploration point data in the grid can be obtained, and whether each seismic exploration point data is in the searching range or not is not needed to be judged, so that the searching times can be greatly reduced, the searching efficiency is improved, and the overall efficiency of static correction of the seismic exploration point data is improved.
Before step 101, a large amount of seismic exploration point data in a work area needs to be obtained, wherein the seismic exploration point comprises an excitation point and a receiving point, the excitation point can be a shot point, the receiving point can be a detection point, and the method screens out the seismic exploration point data taking at least one target seismic exploration point data as a center from the large amount of seismic exploration point data so as to meet the search range. In step 101, seismic survey point data is first divided into a plurality of grids; step 102 is then entered, for each target seismic survey point data, determining a grid corresponding to the target seismic survey point data and conforming to the search range, and obtaining seismic survey point data corresponding to the target seismic survey point data and conforming to the search range; finally, step 103 is entered to perform statics correction on all the obtained seismic survey point data, and the efficiency of the process of obtaining the seismic survey point data requiring statics correction is greatly improved before the statics correction, so that the overall statics correction efficiency is effectively improved.
In particular embodiments, there are a number of ways to divide the seismic survey point data into a plurality of grids, one of which is given below.
In one embodiment, partitioning seismic survey point data into a plurality of grids includes:
according to the coordinate range of the seismic exploration point data, the grid number in the horizontal direction and the grid number in the vertical direction are obtained;
drawing a plurality of grids according to the number of grids in the horizontal direction and the number of grids in the vertical direction;
numbering the grids, and determining the number of the grid where the seismic exploration point data are located.
In the above embodiment, the seismic survey point data is divided into a plurality of grids, that is, a large number of seismic survey point data are grouped, the distribution of the large number of seismic survey point data is within a specific work area, first, the coordinate ranges of all the seismic survey point data are obtained, assuming that the coordinates of N seismic survey point data are divided into (X1, Y1), (X2, Y2), (X3, Y3), … (Xn, yn), the maximum value Xmax and the minimum value Xmin of the X axis (horizontal direction) and the maximum value Ymax and the minimum value Ymin of the Y axis (vertical direction) in these coordinates are found, the coordinate ranges in the horizontal direction are Xmax-Xmin, the coordinate ranges in the vertical direction are Ymax-Ymin, fig. 2 is a schematic diagram of one grid division in the embodiment of the present invention, the maximum value Xmax and the minimum value Xmin of the X axis (horizontal direction) in fig. 2 are 4312 and 4437, and the maximum value Ymax and the minimum value Ymin of the Y axis (vertical direction) are 2945 and 3024, respectively; then, the number of grids in the horizontal direction and the number of grids in the vertical direction are determined according to the coordinate range in the horizontal direction and the coordinate range in the vertical direction, for example, in fig. 2, the number of grids in the horizontal direction is determined to be 25, and the number of grids in the vertical direction is determined to be 8;
then, drawing a plurality of grids according to the number of grids in the horizontal direction and the number of grids in the vertical direction, for example, when the number of grids in the horizontal direction is 25 and the number of grids in the vertical direction is 8, the drawn grids are 25×8=200;
finally, numbering a plurality of grids to determine the number of the grid where the seismic exploration point data is located, fig. 3 is a schematic diagram of the grid number in the embodiment of the invention, the grid number in fig. 3 adopts a two-dimensional representation method, the numbers start from the lower left corner of the grid, increment to the right first and then increment upwards, the number of each grid is (column serial number and row serial number), and the number of each grid where the seismic exploration point data is located is determined in the following manner:
wherein, (mx, my) is the number of the grid;
(Xi, yi) is the coordinates of the seismic survey point data;
xmin and Ymin are respectively a minimum value in the horizontal direction and a minimum value in the vertical direction;
dx and dy are respectively the grid pitch in the horizontal direction and the grid pitch in the vertical direction.
In specific implementation, according to the coordinate range of the seismic exploration point data, a plurality of methods for obtaining the number of grids in the horizontal direction and the number of grids in the vertical direction are provided, and one embodiment is given below.
In one embodiment, obtaining the number of grids in the horizontal direction and the number of grids in the vertical direction according to the coordinate range of the seismic exploration point data comprises:
according to the coordinate range of the seismic exploration point data, grid intervals in the horizontal direction and grid intervals in the vertical direction are obtained;
and obtaining the number of grids in the horizontal direction and the number of grids in the vertical direction according to the grid spacing in the horizontal direction and the grid spacing in the vertical direction.
In the above embodiment, after the grid intervals in the horizontal direction and the grid intervals in the vertical direction are obtained, the grid numbers in the horizontal direction and the grid numbers in the vertical direction can be obtained using the following formulas:
if Nx0 is an integer, then
Nx＝Nx0 (3)
If Nx0 is a fraction, then
Nx＝int(Nx0)+1 (4)
If Ny0 is an integer, then
Ny＝Ny0 (5)
If Ny0 is a decimal fraction then
Ny＝int(Ny0)+1 (6)
Wherein Nx and Ny are the number of grids in the horizontal direction and the number of grids in the vertical direction respectively;
x max and Y max are the maximum value in the horizontal direction and the maximum value in the vertical direction respectively;
dx and dy are respectively the grid pitch in the horizontal direction and the grid pitch in the vertical direction.
The formula can ensure that all the seismic exploration point data fall in the grid, and the condition that the existing seismic exploration point data fall outside the grid can not occur.
In specific implementation, there are various methods for obtaining the grid spacing in the horizontal direction and the grid spacing in the vertical direction according to the coordinate range of the seismic exploration point data, and one of the embodiments is given below.
In an embodiment, if the seismic survey point data are equally spaced apart in the horizontal direction and in the vertical direction, determining that the grid spacing in the horizontal direction is an integer multiple of the distance between the seismic survey point data in the horizontal direction, and determining that the grid spacing in the vertical direction is an integer multiple of the distance between the seismic survey point data in the vertical direction;
if the seismic exploration point data are distributed in a non-equidistant mode in the horizontal direction or the vertical direction, according to the number of the seismic exploration point data, the grid spacing in the horizontal direction and the grid spacing in the vertical direction are obtained according to the coordinate range of the seismic exploration point data.
In the above embodiment, if the seismic survey point data are equally spaced in the horizontal direction and the vertical direction, the seismic survey point data may be directly divided according to the characteristic of equidistant distribution, for example, the seismic survey point data are distributed every 2 in the horizontal direction, and the seismic survey point data are distributed every 3 in the vertical direction, and then the grid spacing in the horizontal direction is an integer multiple of 2, for example, 2,4,6, etc.; the grid spacing in the vertical direction is an integer multiple of 3, for example 3,6,9, etc.
If the seismic exploration point data are distributed in a non-equidistant mode in the horizontal direction or the vertical direction, according to the number of the seismic exploration point data, the grid spacing in the horizontal direction and the grid spacing in the vertical direction are obtained according to the coordinate range of the seismic exploration point data.
In an example, if the seismic prospecting point data are distributed in a non-equidistant manner in the horizontal direction or in the vertical direction, the following formula is adopted, and according to the number of the seismic prospecting point data, the grid spacing in the horizontal direction and the grid spacing in the vertical direction are obtained according to the coordinate range of the seismic prospecting point data:
wherein dx and dy are respectively the grid spacing in the horizontal direction and the grid spacing in the vertical direction;
xLen and yLen are the coordinate range of the seismic exploration point data in the horizontal direction and the coordinate range of the seismic exploration point data in the vertical direction, xlen=x max-X min, ylen=y max-Y min, respectively;
n is the number of seismic exploration point data.
Since the grid is a rectangular grid, the search range is also a rectangular range.
In one embodiment, the search range includes a lateral search radius and a longitudinal search radius.
In an embodiment, for each target seismic survey point data, determining a grid corresponding to the target seismic survey point data and conforming to a search range, and obtaining the seismic survey point data corresponding to the target seismic survey point data and conforming to the search range, includes:
determining the number of the transverse search grids according to the transverse search radius;
determining the number of the longitudinal search grids according to the longitudinal search radius;
determining the number of a grid where each target seismic exploration point data is located;
according to the number of grids where the target seismic exploration point data are located, transversely searching the number of grids, longitudinally searching the number of grids, and determining the number of grids which correspond to the target seismic exploration point data and accord with the searching range;
and obtaining the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the search range according to the serial numbers of the grids which correspond to the target seismic exploration point data and accord with the search range.
In the above embodiment, the number of the lateral search grids is determined according to the lateral search radius by adopting the following formula:
NSearchX＝dSearchX/dx (8)
wherein NSearchX is the number of the transverse search grids;
dSearchX is the lateral search radius.
According to the longitudinal search radius, the number of the longitudinal search grids is determined by adopting the following formula:
NSearchY＝dSearchY/dy (9)
wherein NSearchY is the number of the longitudinal search grids;
dSearchY is the longitudinal search radius.
Then, for each target seismic survey point data, the number of the grid in which the target seismic survey point data is located is determined according to formula (1). According to the number of grids where the target seismic exploration point data is located, the number of grids is searched transversely, the number of grids is searched longitudinally, the number of grids which correspond to the target seismic exploration point data and accord with the searching range is determined, fig. 4 is a schematic diagram of seismic exploration point data searching in the embodiment of the invention, as shown in fig. 4, the M points are target seismic exploration points, and the coarse line frames are screened grids, wherein 5 grids are transversely arranged, and 3 grids are longitudinally arranged.
And finally, according to the number of the grid which corresponds to the target seismic exploration point data and accords with the searching range, the coordinates of the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the searching range can be easily determined, so that the problem of inefficiency caused by global searching is avoided.
In practice, there are various methods for static correction of all the obtained seismic survey points, one of which is given below.
In one embodiment, static correction is performed on all seismic survey point data obtained, including:
determining a static correction amount of the seismic prospecting point data for each obtained seismic prospecting point data;
obtaining a high-frequency component and a low-frequency component of the seismic exploration point data according to the static correction value;
and carrying out static correction on the seismic exploration point data according to the high-frequency component and the low-frequency component of the seismic exploration point data.
In the above embodiment, the static correction amount, which is a function of the spatial positions of the transmitting point (for example, a shot point) and the receiving point (for example, a detector point), is a curve (plane) that varies in space, and can be decomposed into a low-frequency component and a high-frequency component.
The low frequency component, i.e., long wavelength (wavelength longer than alignment length), has less obvious effect on the superposition effect, but affects the exploration of low amplitude structures. The high frequency component, i.e., short wavelength (wavelength less than the alignment length), is static in amount, responsive to co-directional superposition.
Based on the static correction, after obtaining the high frequency component and the low frequency component of the seismic survey point data,
residual static correction of the earth surface consistency mainly solves the problem of short-wavelength static correction; the field static correction and refraction static correction can be used for carrying out static correction on the seismic exploration point data according to the high-frequency component and the low-frequency component of the seismic exploration point data.
With the development of digital processing technology, there are various methods of automatic static correction, such as reference plane static correction, residual static correction, and relative refraction static correction. The reference plane static correction method comprises a folding field Gao Chengjing correction method, a refraction static correction method and a first-arrival static correction method. Fold field Gao Chengjing correction is only applicable to areas where low speed bands are absent or where low speed bands are not laterally altered. The actual assumption condition is difficult to meet, so the method is only applicable to field collection and field processing and indoor processing quality monitoring. The refraction static correction is strictly a model static correction method, and the model is built by the method which is different from the traditional method, mainly comprises the steps of obtaining the speed and delay time of a refraction surface through solving equation inversion, building a speed model by means of surface speed, and completing calculation of a static correction value on the basis. As the first-arrival static correction is the same as the refraction static correction, the model static correction is strictly the model static correction, and the model static correction is not the accurate description of the geophysical attribute of the near-surface medium, a lamellar refraction model cannot be established, and the velocity field distribution of the near-surface layer is obtained through tomography inversion by using the first-arrival time, so that the most excellent model is obtained. Residual static correction generally refers to reflection static correction, and theoretical assumptions and implementation methods of reflection residual static correction determine its unavoidable limitations, so its application must be based on good reference plane static correction. The relative refraction static correction is a method between the static correction of a reference plane and the residual static correction, avoids two conditions which are required to be met by the refraction static correction, only focuses on partial offset refraction waves with better quality, and can obtain more accurate high-frequency components and partial medium-frequency components although an accurate model cannot be built. Different static correction methods can be selected according to actual situations, and are not described herein.
In one embodiment, the seismic survey point data includes excitation points and reception points of the seismic survey.
In the above embodiment, the excitation point may be a shot point, and the receiving point may be a detector point.
Based on the above embodiments, the present invention proposes the following embodiment to explain the detailed flow of the seismic prospecting point data static correction method, and fig. 5 is a detailed flow chart of the seismic prospecting point data static correction method according to the embodiment of the present invention, as shown in fig. 5, in one embodiment, the detailed flow of the seismic prospecting point data static correction method includes:
step 501, if the seismic prospecting point data are equally spaced in the horizontal direction and the vertical direction, determining that the grid spacing in the horizontal direction is an integer multiple of the distance between the seismic prospecting point data in the horizontal direction, and determining that the grid spacing in the vertical direction is an integer multiple of the distance between the seismic prospecting point data in the vertical direction;
step 502, if the seismic exploration point data are distributed in a non-equidistant manner in the horizontal direction or the vertical direction, according to the number of the seismic exploration point data, the coordinate range of the seismic exploration point data is used for obtaining the grid spacing in the horizontal direction and the grid spacing in the vertical direction;
step 503, obtaining the number of grids in the horizontal direction and the number of grids in the vertical direction according to the grid spacing in the horizontal direction and the grid spacing in the vertical direction;
step 504, drawing a plurality of grids according to the number of grids in the horizontal direction and the number of grids in the vertical direction;
step 505, numbering a plurality of grids, and determining the number of the grid where the seismic exploration point data are located;
step 506, determining the number of the transverse search grids according to the transverse search radius;
step 507, determining the number of the longitudinal search grids according to the longitudinal search radius;
step 508, for each target seismic exploration point data, determining the number of the grid where the target seismic exploration point data is located;
step 509, according to the number of the grids where the target seismic exploration point data are located, searching the number of grids horizontally and searching the number of grids longitudinally, and determining the number of the grids which correspond to the target seismic exploration point data and accord with the searching range;
step 510, obtaining the seismic exploration point data corresponding to the target seismic exploration point data and conforming to the search range according to the number of the grid conforming to the search range corresponding to the target seismic exploration point data;
step 511, for each obtained seismic prospecting point data, determining a static correction amount of the seismic prospecting point data;
step 512, obtaining a high-frequency component and a low-frequency component of the seismic prospecting point data according to the static correction value;
step 513, performing static correction on the seismic survey point data based on the high frequency component and the low frequency component of the seismic survey point data.
Of course, it can be understood that other variations of the detailed flow of the seismic prospecting point data static correction method can be provided, and all the related variations should fall within the protection scope of the present invention.
An embodiment is given below to illustrate a specific application of the seismic prospecting point data static correction method according to the embodiment of the present invention.
In the embodiment of the invention, the seismic exploration data points consisting of the shots and the geophones are located in the three-dimensional exploration data of a region in the west of China, the data of the seismic exploration points which need to be subjected to static correction are firstly obtained, the data comprise the shots and the geophones, and then the shots and the geophones are subjected to static correction.
The maximum value of the mass seismic survey point data in the horizontal direction is 4437, the minimum value is 4312, the maximum value in the vertical direction is 3024, and the minimum value is 2945, as shown in fig. 2.
As can be seen from fig. 2, the seismic survey point data are equally spaced in the horizontal direction and the vertical direction, and can be directly divided according to the characteristic of equidistant distribution, wherein the grid pitch in the horizontal direction is 5, and the grid pitch in the vertical direction is 10.
And (3) obtaining the number of grids in the horizontal direction and the number of grids in the vertical direction by adopting a formula (2) -a formula (6) according to the grid spacing in the horizontal direction and the grid spacing in the vertical direction, wherein the grid numbers are 25 and 8 respectively.
Then, drawing a plurality of grids according to the number of grids 25 in the horizontal direction and the number of grids 8 in the vertical direction to form 200 grids;
numbering 200 grids, and determining the number of the grid where the seismic exploration point data are located by adopting a formula (1);
the number of lateral search grids is determined according to a lateral search radius, wherein the lateral search radius is 25, and the number of lateral search grids is determined to be 25/5=5 according to formula (8).
The number of the vertical search grids is determined according to the vertical search radius, wherein the vertical search radius is 30, and the number of the horizontal search grids is determined to be 30/10=3 according to the formula (9).
For each target seismic survey point data, determining the number of the grid in which the target seismic survey point data is located, e.g., the coordinates of the target seismic survey point data are (3413.2, 2945.1), and determining the grid number of the target seismic survey point data as (0, 0) according to formula (1); for another example, the grid number where the M point in fig. 4 is located is (14, 4).
According to the number of the grids where the target seismic exploration point data are located, for example, (14, 4), the number of the grids is searched transversely, the number of the grids is searched longitudinally, the number of the grids which correspond to the target seismic exploration point data and accord with the searching range is determined, and the number of the grids corresponds to the M points.
And obtaining the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the search range according to the number of the grid which corresponds to the target seismic exploration point data and accords with the search range.
Determining a static correction amount of the seismic prospecting point data for each obtained seismic prospecting point data;
obtaining a high-frequency component and a low-frequency component of the seismic exploration point data according to the static correction value;
and carrying out static correction on the seismic exploration point data according to the high-frequency component and the low-frequency component of the seismic exploration point data.
In the seismic prospecting point data static correction method provided by the embodiment of the invention, firstly, the seismic prospecting point data is divided into a plurality of grids; then, determining grids which correspond to the target seismic exploration point data and accord with the search range for each target seismic exploration point data, and obtaining the seismic exploration point data which correspond to the target seismic exploration point data and accord with the search range; and finally, carrying out static correction on all the obtained seismic exploration point data. In the screening process, firstly, the grid which corresponds to the target seismic exploration point data and accords with the searching range can be obtained quickly, then all the seismic exploration point data in the grid can be obtained, and whether each seismic exploration point data is in the searching range or not is not needed to be judged, so that the searching times can be greatly reduced, the searching efficiency is improved, and the overall efficiency of static correction of the seismic exploration point data is improved.
Based on the same inventive concept, the embodiment of the invention also provides a seismic exploration point data static correction device, as described in the following embodiment. Since the principles of solving the problems are similar to those of the seismic prospecting point data static correction method, the implementation of the device can be referred to the implementation of the method, and the repetition is omitted.
FIG. 6 is a schematic diagram of a seismic survey point data statics correction apparatus according to an embodiment of the invention, as shown in FIG. 6, comprising:
a network partitioning module 601 for partitioning seismic survey point data into a plurality of grids;
the screening module 602 is configured to determine, for each target seismic exploration point data, a grid corresponding to the target seismic exploration point data and conforming to a search range, and obtain seismic exploration point data corresponding to the target seismic exploration point data and conforming to the search range;
the static correction module 603 is configured to perform static correction on all the obtained seismic prospecting point data.
In one embodiment, the network dividing module 601 is specifically configured to:
according to the coordinate range of the seismic exploration point data, the grid number in the horizontal direction and the grid number in the vertical direction are obtained;
drawing a plurality of grids according to the number of grids in the horizontal direction and the number of grids in the vertical direction;
numbering the grids, and determining the number of the grid where the seismic exploration point data are located.
In one embodiment, the network dividing module 601 is specifically configured to:
according to the coordinate range of the seismic exploration point data, grid intervals in the horizontal direction and grid intervals in the vertical direction are obtained;
and obtaining the number of grids in the horizontal direction and the number of grids in the vertical direction according to the grid spacing in the horizontal direction and the grid spacing in the vertical direction.
In one embodiment, the network dividing module 601 is specifically configured to:
if the seismic exploration point data are distributed at equal intervals in the horizontal direction and the vertical direction, determining that the grid intervals in the horizontal direction are integral multiples of the distance between the seismic exploration point data in the horizontal direction, and determining that the grid intervals in the vertical direction are integral multiples of the distance between the seismic exploration point data in the vertical direction;
if the seismic exploration point data are distributed in a non-equidistant mode in the horizontal direction or the vertical direction, according to the number of the seismic exploration point data, the grid spacing in the horizontal direction and the grid spacing in the vertical direction are obtained according to the coordinate range of the seismic exploration point data.
In one embodiment, the search range includes a lateral search radius and a longitudinal search radius.
In one embodiment, the screening module 602 is specifically configured to:
determining the number of the transverse search grids according to the transverse search radius;
determining the number of the longitudinal search grids according to the longitudinal search radius;
determining the number of a grid where each target seismic exploration point data is located;
according to the number of grids where the target seismic exploration point data are located, transversely searching the number of grids, longitudinally searching the number of grids, and determining the number of grids which correspond to the target seismic exploration point data and accord with the searching range;
and obtaining the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the search range according to the number of the grid which corresponds to the target seismic exploration point data and accords with the search range.
In one embodiment, the static correction module 603 is specifically configured to:
determining a static correction amount of the seismic prospecting point data for each obtained seismic prospecting point data;
obtaining a high-frequency component and a low-frequency component of the seismic exploration point data according to the static correction value;
and carrying out static correction on the seismic exploration point data according to the high-frequency component and the low-frequency component of the seismic exploration point data.
In one embodiment, the seismic survey point data includes excitation points and reception points of the seismic survey.
In one embodiment, the network partitioning module 601 is specifically configured to:
if the seismic exploration point data are distributed in a non-equidistant mode in the horizontal direction or the vertical direction, the following formula is adopted, and according to the number of the seismic exploration point data and the coordinate range of the seismic exploration point data, the grid spacing in the horizontal direction and the grid spacing in the vertical direction are obtained:
wherein dx and dy are respectively the grid spacing in the horizontal direction and the grid spacing in the vertical direction;
xLen and yLen are the coordinate range of the seismic exploration point data in the horizontal direction and the coordinate range in the vertical direction, respectively;
n is the number of seismic exploration point data.
In the seismic prospecting point data static correction device provided by the embodiment of the invention, firstly, the seismic prospecting point data is divided into a plurality of grids; then, determining grids which correspond to the target seismic exploration point data and accord with the search range for each target seismic exploration point data, and obtaining the seismic exploration point data corresponding to the target seismic exploration point data; and finally, carrying out static correction on all the obtained seismic exploration point data. In the screening process, firstly, the grid which corresponds to the target seismic exploration point data and accords with the searching range can be obtained quickly, then all the seismic exploration point data in the grid can be obtained, and whether each seismic exploration point data is in the searching range or not is not needed to be judged, so that the searching times can be greatly reduced, the searching efficiency is improved, and the overall efficiency of static correction of the seismic exploration point data is improved.
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, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the 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.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the invention, and is not meant to limit the scope of the invention, but to limit the invention to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (10)
1. A method for statics correction of seismic survey point data, comprising:
dividing seismic survey point data into a plurality of grids;
determining grids which correspond to the target seismic exploration point data and accord with the search range for each target seismic exploration point data, and obtaining the seismic exploration point data which correspond to the target seismic exploration point data and accord with the search range;
performing static correction on all obtained seismic exploration point data;
the search range comprises a transverse search radius and a longitudinal search radius;
for each target seismic exploration point data, determining a grid which corresponds to the target seismic exploration point data and accords with a search range, and obtaining the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the search range, wherein the method comprises the following steps: determining the number of the transverse search grids according to the transverse search radius; determining the number of the longitudinal search grids according to the longitudinal search radius; determining the number of a grid where each target seismic exploration point data is located; according to the number of grids where the target seismic exploration point data are located, transversely searching the number of grids, longitudinally searching the number of grids, and determining the number of grids which correspond to the target seismic exploration point data and accord with the searching range; and obtaining the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the search range according to the number of the grid which corresponds to the target seismic exploration point data and accords with the search range.
2. The method of statics correction of seismic survey point data of claim 1, wherein dividing the seismic survey point data into a plurality of grids comprises:
according to the coordinate range of the seismic exploration point data, the grid number in the horizontal direction and the grid number in the vertical direction are obtained;
drawing a plurality of grids according to the number of grids in the horizontal direction and the number of grids in the vertical direction;
numbering the grids, and determining the number of the grid where the seismic exploration point data are located.
3. The method of statics correction of seismic survey point data according to claim 2, wherein obtaining the number of grids in the horizontal direction and the number of grids in the vertical direction from the coordinate range of the seismic survey point data comprises:
according to the coordinate range of the seismic exploration point data, grid intervals in the horizontal direction and grid intervals in the vertical direction are obtained;
and obtaining the number of grids in the horizontal direction and the number of grids in the vertical direction according to the grid spacing in the horizontal direction and the grid spacing in the vertical direction.
4. A method of statics correction of seismic survey point data as claimed in claim 3, wherein obtaining the grid spacing in the horizontal direction and the grid spacing in the vertical direction from the coordinate range of the seismic survey point data comprises:
if the seismic exploration point data are distributed at equal intervals in the horizontal direction and the vertical direction, determining that the grid intervals in the horizontal direction are integral multiples of the distance between the seismic exploration point data in the horizontal direction, and determining that the grid intervals in the vertical direction are integral multiples of the distance between the seismic exploration point data in the vertical direction;
if the seismic exploration point data are distributed in a non-equidistant mode in the horizontal direction or the vertical direction, according to the number of the seismic exploration point data, the grid spacing in the horizontal direction and the grid spacing in the vertical direction are obtained according to the coordinate range of the seismic exploration point data.
5. The method for statics correction of seismic survey point data according to claim 1, wherein statics correction is performed on all obtained seismic survey point data, comprising:
determining a static correction amount of the seismic prospecting point data for each obtained seismic prospecting point data;
obtaining a high-frequency component and a low-frequency component of the seismic exploration point data according to the static correction value;
and carrying out static correction on the seismic exploration point data according to the high-frequency component and the low-frequency component of the seismic exploration point data.
6. The method of seismic survey point data statics correction of claim 1, wherein the seismic survey point data comprises excitation points and reception points of a seismic survey.
7. The method for statics correction of seismic survey point data according to claim 1, wherein if the seismic survey point data are distributed in a non-equidistant manner in a horizontal direction or in a vertical direction, the grid spacing in the horizontal direction and the grid spacing in the vertical direction are obtained from the coordinate range of the seismic survey point data based on the number of the seismic survey point data by using the following formula:
wherein dx and dy are respectively the grid spacing in the horizontal direction and the grid spacing in the vertical direction;
xLen and yLen are the coordinate range of the seismic exploration point data in the horizontal direction and the coordinate range in the vertical direction, respectively;
n is the number of seismic exploration point data.
8. A seismic survey point data statics correction apparatus comprising:
the network dividing module is used for dividing the seismic exploration point data into a plurality of grids;
the screening module is used for determining grids which correspond to the target seismic exploration point data and accord with the search range for each target seismic exploration point data, and obtaining the seismic exploration point data which correspond to the target seismic exploration point data and accord with the search range;
the static correction module is used for carrying out static correction on all the obtained seismic exploration point data;
the search range comprises a transverse search radius and a longitudinal search radius;
for each target seismic exploration point data, determining a grid which corresponds to the target seismic exploration point data and accords with a search range, and obtaining the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the search range, wherein the method comprises the following steps: determining the number of the transverse search grids according to the transverse search radius; determining the number of the longitudinal search grids according to the longitudinal search radius; determining the number of a grid where each target seismic exploration point data is located; according to the number of grids where the target seismic exploration point data are located, transversely searching the number of grids, longitudinally searching the number of grids, and determining the number of grids which correspond to the target seismic exploration point data and accord with the searching range; and obtaining the seismic exploration point data which corresponds to the target seismic exploration point data and accords with the search range according to the number of the grid which corresponds to the target seismic exploration point data and accords with the search range.
9. 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 the method of any of claims 1 to 7 when executing the computer program.
10. A computer readable storage medium, characterized in that the computer readable storage medium stores a computer program for executing the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910520558.5A CN112099091B (en) | 2019-06-17 | 2019-06-17 | Static correction method and device for seismic exploration point data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910520558.5A CN112099091B (en) | 2019-06-17 | 2019-06-17 | Static correction method and device for seismic exploration point data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112099091A CN112099091A (en) | 2020-12-18 |
| CN112099091B true CN112099091B (en) | 2024-03-01 |
Family
ID=73748407
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910520558.5A Active CN112099091B (en) | 2019-06-17 | 2019-06-17 | Static correction method and device for seismic exploration point data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112099091B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6075752A (en) * | 1995-01-23 | 2000-06-13 | Elf Aquitaine Production | Method for exploring an environment with a complex tectonics |
| CN101980055A (en) * | 2010-11-30 | 2011-02-23 | 中国石油集团川庆钻探工程有限公司 | Digital modeling method for seismic exploration |
| CN102819042A (en) * | 2011-06-10 | 2012-12-12 | 中国石油天然气集团公司 | Near-surface chromatography inversion static rectification method of rarefied seismic data |
| CN103605163A (en) * | 2013-10-28 | 2014-02-26 | 中国石油天然气集团公司 | Static-correction-value quantitative evaluation method and device |
| CN104122588A (en) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Spectral decomposition based post-stack seismic data resolution ratio increasing method |
| CN105759309A (en) * | 2016-02-29 | 2016-07-13 | 中国石油集团东方地球物理勘探有限责任公司 | Method and apparatus for making geological structure map |
| CN107561574A (en) * | 2017-07-20 | 2018-01-09 | 中国石油天然气集团公司 | Shot position, the detection method of receiving point position and device |
| CN107843928A (en) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | A kind of shot point method for correcting error based on grid-search algorithms |
| CN109581503A (en) * | 2017-09-28 | 2019-04-05 | 中国石油化工股份有限公司 | A kind of method and system solving desert exploratory area tomographic statics border issue |
| RU2690977C1 (en) * | 2018-07-20 | 2019-06-07 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Кубанский государственный аграрный университет имени И.Т. Трубилина" | Method for prediction of zones of development of secondary reservoirs of a fracture type in a sedimentary cover for search of hydrocarbon deposits |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190094397A1 (en) * | 2017-09-25 | 2019-03-28 | Real Time Geomechanics, Llc | Surface detection and location of microseismic events and earthquakes without the use of a velocity model |
-
2019
- 2019-06-17 CN CN201910520558.5A patent/CN112099091B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6075752A (en) * | 1995-01-23 | 2000-06-13 | Elf Aquitaine Production | Method for exploring an environment with a complex tectonics |
| CN101980055A (en) * | 2010-11-30 | 2011-02-23 | 中国石油集团川庆钻探工程有限公司 | Digital modeling method for seismic exploration |
| CN102819042A (en) * | 2011-06-10 | 2012-12-12 | 中国石油天然气集团公司 | Near-surface chromatography inversion static rectification method of rarefied seismic data |
| CN103605163A (en) * | 2013-10-28 | 2014-02-26 | 中国石油天然气集团公司 | Static-correction-value quantitative evaluation method and device |
| CN104122588A (en) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Spectral decomposition based post-stack seismic data resolution ratio increasing method |
| CN105759309A (en) * | 2016-02-29 | 2016-07-13 | 中国石油集团东方地球物理勘探有限责任公司 | Method and apparatus for making geological structure map |
| CN107843928A (en) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | A kind of shot point method for correcting error based on grid-search algorithms |
| CN107561574A (en) * | 2017-07-20 | 2018-01-09 | 中国石油天然气集团公司 | Shot position, the detection method of receiving point position and device |
| CN109581503A (en) * | 2017-09-28 | 2019-04-05 | 中国石油化工股份有限公司 | A kind of method and system solving desert exploratory area tomographic statics border issue |
| RU2690977C1 (en) * | 2018-07-20 | 2019-06-07 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Кубанский государственный аграрный университет имени И.Т. Трубилина" | Method for prediction of zones of development of secondary reservoirs of a fracture type in a sedimentary cover for search of hydrocarbon deposits |
Non-Patent Citations (7)
| Title |
|---|
| 《3D Multi-scale Tomostatics using First Arrivals》;Huawei Zhou et al.;IPTC;第1-17页 * |
| 《COMMON-REFLECTION-SURFACE STACK:ACCOUNTING FOR CONFLICTING DIP SITUATIONS BY CONSIDERING ALL POSSIBLES DIPS》;MEHRDAD SOLEIMANI et al.;JOURNAL OF SEISMIC EXPLORATION;第271-288页 * |
| 《一种高精度初至波二次定位新方法――搜索法》;倪成洲等;石油地球物理勘探;第43卷(第2期);第131-133、157页 * |
| 《地震反演与非线性随机优化方法》;尹成等;物探化探计算技术;第23卷(第1期);第6-10页 * |
| 《复杂区三维折射静校正技术与应用效果》;杨海申等;石油地球物理勘探;第40卷(第2期);第219-225页 * |
| 《粒子群算法炮点偏移校正技术》;朱海伟等;石油地球物理勘探;第51卷(第4期);第661-664页 * |
| 《非规则网带断层地震数据的网格化》;张美根等;地球物理学进展;第14卷(第1期);第69-76页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112099091A (en) | 2020-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110907995B (en) | Reverse time migration method and device for VSP seismic data in well | |
| CN105093301B (en) | The generation method and device of common imaging point angle of reflection angle gathers | |
| CN107462924B (en) | A kind of absolute wave impedance inversion method independent of well-log information | |
| EA029700B1 (en) | Method of processing seismic image of the subsurface | |
| EP2877880B1 (en) | System and method for migration velocity modeling | |
| WO2014164970A2 (en) | Methods and systems for locating seismic events | |
| CN104483705A (en) | Three-dimensional residual static correction method | |
| CN110031898B (en) | Data optimization method and integral method prestack depth migration method | |
| CN116774292B (en) | Seismic wave travel time determining method, system, electronic equipment and storage medium | |
| CN106959467B (en) | Seimic wave velocity inversion method and device | |
| CN112305591B (en) | Tunnel advanced geological prediction method and computer readable storage medium | |
| CN109633749A (en) | Non-linear Fresnel zone seismic traveltime tomography method based on scattering integral method | |
| CN113552621B (en) | Shale gas ground stress determination method and device | |
| CN106291698B (en) | Seismic facies sedimentary facies determines method and apparatus | |
| CN112099091B (en) | Static correction method and device for seismic exploration point data | |
| CN107831536A (en) | Desert surface seismic data processing floating datum method for building up | |
| CN107340539A (en) | A kind of earthquake data inversion method and system | |
| CN108375794B (en) | VSP (vertical seismic profiling) slit-hole diffraction imaging technical method based on symmetrical observation | |
| CN103513279A (en) | Lighting analysis calculation method and device based on wave equation of earthquake waves | |
| CN104570062B (en) | A kind of VSP observation system design methods centered on exciting | |
| CN115826038A (en) | VTI medium adjoint state method travel time multi-parameter tomography method and system | |
| CN109387872A (en) | Surface-related multiple prediction technique | |
| RU2221262C1 (en) | Process of seismic prospecting for examination of sedimentary mantle in presence of heavily broken-up acoustically stiff boundaries ( variants ) | |
| CN113917533A (en) | Systematic implementation method of double-linkage omnibearing imaging of TI medium | |
| CN104537246B (en) | A kind of observation system imaging capability evaluation method |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |