CN111352157B - Shear wave static correction method and system - Google Patents
Shear wave static correction method and system Download PDFInfo
- Publication number
- CN111352157B CN111352157B CN201811567259.9A CN201811567259A CN111352157B CN 111352157 B CN111352157 B CN 111352157B CN 201811567259 A CN201811567259 A CN 201811567259A CN 111352157 B CN111352157 B CN 111352157B
- Authority
- CN
- China
- Prior art keywords
- data
- set data
- shot gather
- time shift
- shift amount
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000001514 detection method Methods 0.000 claims description 32
- 230000009466 transformation Effects 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000002994 Monte Carlo simulated annealing Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a shear wave static correction method and a shear wave static correction system. The shear wave static correction method comprises the following steps: obtaining p set data according to the original shot set data; obtaining model channel data according to the p sets of data; selecting a plurality of first time windows according to the transverse wave data in the p set data; acquiring p set data and model channel data in each first time window, performing cross correlation to obtain a first time shift amount of each first time window, and further obtaining a second time shift amount; correcting the p-set data according to the second time shift amount, and obtaining new shot gather data according to the corrected p-set data; selecting a plurality of second time windows according to converted wave data in the new shot gather data; acquiring original shot gather data and new shot gather data in each second time window, performing cross correlation to obtain a third time shift amount of each second time window, and further obtaining a fourth time shift amount; and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data. The invention has more accurate calculation result and can effectively maintain the fidelity of the data.
Description
Technical Field
The invention relates to the technical field of oil exploration, in particular to a shear wave static correction method and a shear wave static correction system.
Background
In seismic exploration, time shift is caused between recorded in-phase axes due to the influence of changes in the elevation of the earth surface, the thickness of a weathered layer, the velocity of the weathered layer and a reference surface, so that subsequent velocity analysis and structural imaging are influenced. The correction can be used to eliminate the time-shifting effect of the near-surface on the recording of the same-phase axis.
In multi-wave, multi-component seismic data acquisition, the longitudinal and transverse waves that reach the surface are recorded, requiring correction for each wave. If a compressional source is used in seismic data acquisition, most shear energy recorded at the surface is the result of compressional-shear conversion, which occurs where rock impedance changes. Therefore, in order to eliminate all the influence of the near-surface on the longitudinal and transverse waves, it is necessary to calculate the longitudinal wave static correction amount of the shot point and the transverse wave static correction amount of the demodulator probe.
The velocity of the cross wave is low due to the low near surface formation pressure. Therefore, the influence of near-surface shear waves on the travel of longitudinal and shear waves is large. In general, when longitudinal statics correction is performed, it is assumed that the vertical ray path and the near-surface are fast in speed change, but these simplifications are not appropriate when calculating the shear statics correction amount. The shear wave velocity variation characteristics near the surface are different from the longitudinal wave velocity variation characteristics because the shear wave velocity is not significantly affected by the fluid in the pores. The near-surface shear wave velocity gradient is lower than the longitudinal wave velocity gradient. Theory and practice has shown that a large range of incident angles is required to observe strong cross-wave energy, which is inconsistent with the assumption of a near-surface vertical ray path. All of these features make shear wave statics a complex problem.
In practical applications, the most common and simple static correction method for transverse waves is a proportionality coefficient method, that is, the static correction amount of longitudinal waves is multiplied by a proportionality coefficient, and the proportionality coefficient is the velocity ratio of the longitudinal waves and the transverse waves at the near-surface position. The velocity ratio may be derived from surface wave data analysis, shallow borehole data, or reflected wave data. The literature mentions that the range of variation of the ratio of the longitudinal to the transverse wave velocities is between 2 and 6. The results with this method are inaccurate due to the particularity of the propagation of each wave mode, which results in the correction of the longitudinal and transverse waves being independent of each other.
Under the assumption of surface consistency, there are several methods for solving shear wave static correction. These methods include: and a manual picking correction method, a time difference refraction analysis method, generalized interchange inversion and a Monte Carlo simulated annealing method are superposed on the common detection wave points. All of these methods rely on time picking for either the shear wave refractive or the shear wave reflective in-phase axes. At the beginning of data processing, it is difficult to identify these in-phase axes, and it takes a lot of time.
Under the assumption of non-surface consistency, transformation-based methods such as radial track transformation shear wave static correction method and tau-p domain shear wave static correction method are proposed. When the change of the surface speed is large, the adaptability of the radial track transformation shear wave static correction method is poor. For the existing tau-p domain shear wave static correction method, after positive and negative tau-p conversion and smoothing and convolution processing, the fidelity of data cannot be effectively maintained.
In recent years, much attention has been paid to a method for performing shear static correction by inverting the frequency dispersion of a surface wave to obtain near-surface features. Performing multi-spectral analysis (MASW) on the surface wave to establish a near-surface transverse wave velocity model, and further calculating a correction value. However, surface wave multispectral analysis lacks sufficient resolution when the transverse wave velocity varies strongly laterally or when the near-surface is non-horizontally stratified. Also, the surface wave inversion depth is limited to the lowest frequency that is reliable in the data.
In the above methods, the different methods are limited differently, and the calculation result has the disadvantages of low precision, time-consuming manual picking, low resolution, low fidelity rate, etc., which seriously affects the subsequent data processing.
Disclosure of Invention
The embodiment of the invention mainly aims to provide a transverse wave static correction method and a transverse wave static correction system, which are more accurate in calculation result, high in calculation speed and high in transverse resolution and can effectively keep the fidelity of data.
In order to achieve the above object, an embodiment of the present invention provides a shear wave static correction method, including:
obtaining p set data according to the original shot set data;
obtaining model channel data corresponding to the p set data according to the p set data;
selecting a plurality of first time windows according to the transverse wave data in the p set data;
acquiring p set data and model channel data in each first time window;
performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift quantity of each first time window;
obtaining a second time shift amount according to the first time shift amount of each first time window;
correcting the p set data according to the second time shift amount to obtain corrected p set data;
obtaining new shot gather data according to the corrected p set data;
selecting a plurality of second time windows according to converted wave data in the new shot gather data;
acquiring original shot gather data and new shot gather data in each second time window;
performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
obtaining a fourth time shift amount according to the third time shift amount of each second time window;
and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
The embodiment of the present invention further provides a shear wave static correction system, including:
the p set data module is used for obtaining p set data according to the original shot set data;
the model track data module is used for obtaining model track data corresponding to the p set data according to the p set data;
the first time window module is used for selecting a plurality of first time windows according to the transverse wave data in the p set data;
the first acquisition module is used for acquiring the p set data and the model channel data in each first time window;
the first time shift module is used for performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift of each first time window;
the second time shift module is used for obtaining a second time shift according to the first time shift of each first time window;
the first correction module is used for correcting the p set data according to the second time shift amount to obtain corrected p set data;
the new shot gather data module is used for obtaining new shot gather data according to the corrected p set data;
the second time window module is used for selecting a plurality of second time windows according to the converted wave data in the new shot gather data;
the second acquisition module is used for acquiring the original shot gather data and the new shot gather data in each second time window;
the third time shift module is used for performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift of each second time window;
the fourth time shift module is used for obtaining a fourth time shift according to the third time shift of each second time window;
and the second correction module is used for correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
An embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, where the processor implements the following steps when executing the computer program:
obtaining p set data according to the original shot set data;
obtaining model channel data corresponding to the p set data according to the p set data;
selecting a plurality of first time windows according to shear wave data in the p sets of data;
acquiring p set data and model channel data in each first time window;
performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift quantity of each first time window;
obtaining a second time shift amount according to the first time shift amount of each first time window;
correcting the p set data according to the second time shift amount to obtain corrected p set data;
obtaining new shot gather data according to the corrected p set data;
selecting a plurality of second time windows according to converted wave data in the new shot gather data;
acquiring original shot gather data and new shot gather data in each second time window;
performing cross correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
obtaining a fourth time shift amount according to the third time shift amount of each second time window;
and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the following steps:
obtaining p set data according to the original shot set data;
obtaining model channel data corresponding to the p set data according to the p set data;
selecting a plurality of first time windows according to the transverse wave data in the p set data;
acquiring p set data and model channel data in each first time window;
performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift amount of each first time window;
obtaining a second time shift amount according to the first time shift amount of each first time window;
correcting the p set data according to the second time shift amount to obtain corrected p set data;
obtaining new shot gather data according to the corrected p set data;
selecting a plurality of second time windows according to converted wave data in the new shot gather data;
acquiring original shot gather data and new shot gather data in each second time window;
performing cross correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
obtaining a fourth time shift amount according to the third time shift amount of each second time window;
and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
According to the shear wave static correction method and system, p set data are obtained according to original shot set data, a second time shift amount is obtained according to the p set data, the p set data are corrected according to the second time shift amount, and new shot set data are obtained according to the corrected p set data; and obtaining a fourth time shift amount according to the new shot gather data and the original shot gather data, and finally correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without creative efforts.
FIG. 1 is a flow chart of a shear wave static correction method according to an embodiment of the present invention;
fig. 2 is a flowchart of S101 in the embodiment of the present invention;
FIG. 3 is a flowchart of S108 in an embodiment of the present invention;
FIG. 4 is a schematic illustration of raw shot gather data undergoing a preprocessing operation in an embodiment of the present invention;
FIG. 5 is a diagram of p-set data in an embodiment of the present invention;
FIG. 6 is a schematic illustration of corrected p-set data in an embodiment of the invention;
FIG. 7 is a schematic illustration of corrected shot gather data in an embodiment of the present invention;
fig. 8 is a block diagram of the structure of the shear wave static correction system in the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As will be appreciated by one skilled in the art, embodiments of the present invention may be embodied as a system, apparatus, device, method, or computer program product. Accordingly, the present disclosure may be embodied in the form of: entirely hardware, entirely software (including firmware, resident software, micro-code, etc.), or a combination of hardware and software.
In view of the defects of low precision, time-consuming manual picking, low resolution, low fidelity and the like of the calculation result obtained in the prior art, the embodiment of the invention provides a transverse wave static correction method, which has the advantages of more accurate calculation result, high calculation speed, high transverse resolution and capability of effectively maintaining the fidelity of data. The present invention will be described in detail below with reference to the accompanying drawings.
FIG. 1 is a flow chart of a shear wave static correction method according to an embodiment of the present invention. As shown in fig. 1, the shear wave static correction method includes:
s101: and obtaining p set data according to the original shot set data.
S102: and obtaining model trace data corresponding to the p set data according to the p set data.
S103: and selecting a plurality of first time windows according to the shear wave data in the p sets of data.
S104: and acquiring p set data and model trace data in each first time window.
S105: and performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift of each first time window.
S106: and obtaining a second time shift amount according to the first time shift amount of each first time window. The second time shift amount may be an average value of the plurality of first time shift amounts.
S107: and correcting the p set data according to the second time shift amount to obtain corrected p set data.
S108: and obtaining new shot gather data according to the corrected p set data.
S109: and selecting a plurality of second time windows according to the converted wave data in the new shot gather data.
S110: and acquiring the original shot gather data and the new shot gather data in each second time window.
S111: and performing cross correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window.
S112: and obtaining a fourth time shift amount according to the third time shift amount of each second time window. The fourth time shift amount may be an average value of a plurality of third time shift amounts.
S113: and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
The main body for executing the shear wave static correction method shown in fig. 1 may be a computer. As can be seen from the process shown in fig. 1, in the shear wave static correction method according to the embodiment of the present invention, p-set data is obtained according to the original shot gather data, a second time shift amount is obtained according to the p-set data, the p-set data is corrected according to the second time shift amount, and then new shot gather data is obtained according to the corrected p-set data; and obtaining a fourth time shift amount according to the new shot gather data and the original shot gather data, and finally correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
Fig. 2 is a flowchart of S101 in the embodiment of the present invention. As shown in fig. 2, S101 includes:
s201: and sorting the original shot gather data to obtain first detection point set data.
S202: the first set of detection points data is subjected to a tau-p transform.
S203: and sorting the first detection point set data after tau-p conversion to obtain p set data.
S201 specifically comprises the steps of carrying out preprocessing operations such as defining an observation system and denoising on the original shot gather data, and then sorting the preprocessed original shot gather data to obtain first detection point set data.
After executing S101, further comprising: and denoising the p set data according to the signal-to-noise ratio of each channel of data in the p set data.
In an embodiment, S102 specifically includes:
performing superposition processing on the p set data to obtain model channel data corresponding to the p set data; or the like, or, alternatively,
and smoothing the p set data to obtain model trace data corresponding to the p set data.
Fig. 3 is a flowchart of S108 in the embodiment of the present invention. As shown in fig. 3, S108 includes:
s301: sorting the corrected p set data to obtain second detection point set data.
S302: and carrying out tau-p inverse transformation on the second detection point set data.
S303: and sorting the second detection point set data subjected to tau-p inverse transformation to obtain new shot set data.
One of the specific embodiments of the present invention is as follows:
1. and carrying out pretreatment operations such as defining an observation system and denoising on the original shot gather data, and then sorting the pretreated original shot gather data to obtain first detection point set data. FIG. 4 is a schematic illustration of raw shot gather data undergoing a preprocessing operation in an embodiment of the present invention. As shown in FIG. 4, the horizontal axis in FIG. 4 is the track number, and is dimensionless; the vertical axis is time in milliseconds.
2. And carrying out tau-p conversion on the first detection point set data, and sorting the tau-p converted first detection point set data to obtain p set data. FIG. 5 is a diagram of p-set data in an embodiment of the invention. As shown in FIG. 5, the horizontal axis in FIG. 5 is the track number, and is dimensionless; the vertical axis is time in milliseconds.
3. And denoising the p set data according to the signal-to-noise ratio of each channel of data in the p set data.
4. And performing superposition processing or smoothing processing on the p set data subjected to denoising processing to obtain model channel data corresponding to the p set data, and selecting a plurality of first time windows according to transverse wave data in the p set data.
5. And acquiring p set data and model trace data in each first time window, and performing cross correlation on the p set data and the model trace data in each first time window to obtain a first time shift of each first time window.
6. And obtaining a second time shift amount according to the first time shift amount of each first time window, and correcting the p set data according to the second time shift amount to obtain corrected p set data. FIG. 6 is a schematic of corrected p-set data in an embodiment of the invention. As shown in FIG. 6, the horizontal axis in FIG. 6 is the track number, and is dimensionless; the vertical axis is time in milliseconds.
The second time shift amount may be an average value of the plurality of first time shift amounts.
7. Sorting the corrected p set data to obtain second detection point set data.
8. And carrying out tau-p inverse transformation on the second detection point set data, and sorting the second detection point set data subjected to tau-p inverse transformation to obtain new shot gather data.
9. And selecting a plurality of second time windows according to the converted wave data in the new shot gather data, and acquiring the original shot gather data and the new shot gather data in each second time window.
10. And performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window, and obtaining a fourth time shift amount according to the third time shift amount of each second time window. The fourth time shift amount may be an average value of a plurality of third time shift amounts.
11. And correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data. FIG. 7 is a schematic illustration of corrected shot gather data in an embodiment of the invention. As shown in FIG. 7, the horizontal axis in FIG. 7 is the track number, dimensionless; the vertical axis is time in milliseconds.
To sum up, the shear wave static correction method of the embodiment of the invention obtains p-set data according to the original shot gather data, obtains a second time shift amount according to the p-set data, corrects the p-set data according to the second time shift amount, and obtains new shot gather data according to the corrected p-set data; and obtaining a fourth time shift amount according to the new shot gather data and the original shot gather data, and finally correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data. Therefore, the invention has wide application prospect in processing and explaining the multi-wave and multi-component earthquake.
Based on the same inventive concept, the embodiment of the invention also provides a shear wave static correction system, and as the principle of solving the problems of the system is similar to that of the shear wave static correction method, the implementation of the system can refer to the implementation of the method, and repeated parts are not described again.
Fig. 8 is a block diagram of the structure of the shear wave static correction system in the embodiment of the present invention. As shown in fig. 8, the shear wave static correction system includes:
the p set data module is used for obtaining p set data according to the original shot set data;
the model track data module is used for obtaining model track data corresponding to the p set data according to the p set data;
the first time window module is used for selecting a plurality of first time windows according to the transverse wave data in the p set data;
the first acquisition module is used for acquiring the p set data and the model channel data in each first time window;
the first time shift module is used for performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift of each first time window;
the second time shift module is used for obtaining a second time shift according to the first time shift of each first time window;
the first correction module is used for correcting the p set data according to the second time shift amount to obtain corrected p set data;
the new shot gather data module is used for obtaining new shot gather data according to the corrected p set data;
the second time window module is used for selecting a plurality of second time windows according to the converted wave data in the new shot gather data;
the second acquisition module is used for acquiring the original shot gather data and the new shot gather data in each second time window;
the third time shift amount module is used for performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
the fourth time shift module is used for obtaining a fourth time shift according to the third time shift of each second time window;
and the second correction module is used for correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
In one embodiment, the p-set data module is specifically configured to:
sorting original shot gather data to obtain first detection point set data;
carrying out tau-p transformation on the first detection point set data;
and sorting the first detection point set data after tau-p conversion to obtain p set data.
In one embodiment, the new shot gather data module is specifically configured to:
sorting the corrected p set data to obtain second detection point set data;
carrying out tau-p inverse transformation on the second detection point set data;
and sorting the second detection point set data subjected to tau-p inverse transformation to obtain new shot set data.
In one embodiment, the model trace data module is specifically configured to:
performing superposition processing on the p set data to obtain model channel data corresponding to the p set data; or the like, or a combination thereof,
and smoothing the p set data to obtain model channel data corresponding to the p set data.
To sum up, the shear wave static correction system of the embodiment of the invention obtains p-set data according to the original shot gather data, obtains a second time shift amount according to the p-set data, corrects the p-set data according to the second time shift amount, and obtains new shot gather data according to the corrected p-set data; and obtaining a fourth time shift amount according to the new shot gather data and the original shot gather data, and finally correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data. Therefore, the invention has wide application prospect in processing and explaining the multi-wave and multi-component earthquake.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and when the processor executes the computer program, the following steps are implemented:
obtaining p set data according to the original shot set data;
obtaining model channel data corresponding to the p set data according to the p set data;
selecting a plurality of first time windows according to shear wave data in the p sets of data;
acquiring p set data and model channel data in each first time window;
performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift amount of each first time window;
obtaining a second time shift amount according to the first time shift amount of each first time window;
correcting the p set data according to the second time shift amount to obtain corrected p set data;
obtaining new shot gather data according to the corrected p set data;
selecting a plurality of second time windows according to converted wave data in the new shot gather data;
acquiring original shot gather data and new shot gather data in each second time window;
performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
obtaining a fourth time shift amount according to the third time shift amount of each second time window;
and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
To sum up, the computer device of the embodiment of the present invention obtains p-set data according to the original shot gather data, obtains a second time shift amount according to the p-set data, corrects the p-set data according to the second time shift amount, and obtains new shot gather data according to the corrected p-set data; and obtaining a fourth time shift amount according to the new shot gather data and the original shot gather data, and finally correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data. Therefore, the invention has wide application prospect in processing and explaining the multi-wave and multi-component earthquake.
An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the following steps:
obtaining p set data according to the original shot set data;
obtaining model channel data corresponding to the p set data according to the p set data;
selecting a plurality of first time windows according to the transverse wave data in the p set data;
acquiring p set data and model channel data in each first time window;
performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift quantity of each first time window;
obtaining a second time shift amount according to the first time shift amount of each first time window;
correcting the p set data according to the second time shift amount to obtain corrected p set data;
obtaining new shot gather data according to the corrected p set data;
selecting a plurality of second time windows according to converted wave data in the new shot gather data;
acquiring original shot gather data and new shot gather data in each second time window;
performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
obtaining a fourth time shift amount according to the third time shift amount of each second time window;
and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
To sum up, the computer-readable storage medium according to the embodiment of the present invention first obtains p-set data according to the original shot gather data, then obtains a second time shift amount according to the p-set data, corrects the p-set data according to the second time shift amount, and then obtains new shot gather data according to the corrected p-set data; and obtaining a fourth time shift amount according to the new shot gather data and the original shot gather data, and finally correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data. Therefore, the invention has wide application prospect in processing and explaining the multi-wave and multi-component earthquake.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A shear wave static correction method is characterized by comprising the following steps:
obtaining p set data according to the original shot set data;
obtaining model channel data corresponding to the p set data according to the p set data;
selecting a plurality of first time windows according to the shear wave data in the p sets of data;
acquiring p set data and model channel data in each first time window;
performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift amount of each first time window;
obtaining a second time shift amount according to the first time shift amount of each first time window;
correcting the p set data according to the second time shift amount to obtain corrected p set data;
obtaining new shot gather data according to the corrected p set data;
selecting a plurality of second time windows according to converted wave data in the new shot gather data;
acquiring original shot gather data and new shot gather data in each second time window;
performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
obtaining a fourth time shift amount according to the third time shift amount of each second time window;
and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
2. The shear wave statics correction method of claim 1, wherein deriving p-set data from raw shot-set data comprises:
sorting the original shot gather data to obtain first detection point set data;
carrying out tau-p transformation on the first detection point set data;
and sorting the first detection point set data after tau-p conversion to obtain p set data.
3. The shear wave statics method of claim 1, wherein deriving new shot gather data from said corrected p-set data comprises:
sorting the corrected p set data to obtain second detection point set data;
carrying out tau-p inverse transformation on the second detection point set data;
and sorting the second detection point set data subjected to tau-p inverse transformation to obtain the new shot set data.
4. The shear wave static correction method of claim 1, wherein obtaining model trace data corresponding to the p-set data from the p-set data comprises:
performing superposition processing on the p set data to obtain model channel data corresponding to the p set data; or the like, or, alternatively,
and smoothing the p set data to obtain model channel data corresponding to the p set data.
5. A shear wave static correction system, comprising:
the p set data module is used for obtaining p set data according to the original shot set data;
the model track data module is used for obtaining model track data corresponding to the p set data according to the p set data;
the first time window module is used for selecting a plurality of first time windows according to the transverse wave data in the p set data;
the first acquisition module is used for acquiring the p set data and the model channel data in each first time window;
the first time shift module is used for performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift of each first time window;
the second time shift module is used for obtaining a second time shift according to the first time shift of each first time window;
the first correction module is used for correcting the p set data according to the second time shift amount to obtain corrected p set data;
the new shot gather data module is used for obtaining new shot gather data according to the corrected p set data;
the second time window module is used for selecting a plurality of second time windows according to converted wave data in the new shot gather data;
the second acquisition module is used for acquiring the original shot gather data and the new shot gather data in each second time window;
the third time shift module is used for performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift of each second time window;
the fourth time shift module is used for obtaining a fourth time shift according to the third time shift of each second time window;
and the second correction module is used for correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
6. The shear wave static correction system of claim 5, wherein the p-set data module is specifically configured to:
sorting the original shot gather data to obtain first detection point set data;
carrying out tau-p transformation on the first detection point set data;
and sorting the first detection point set data after tau-p conversion to obtain p set data.
7. The shear wave statics correction system of claim 5, wherein the new shot gather data module is specifically configured to:
sorting the corrected p set data to obtain second detection point set data;
carrying out tau-p inverse transformation on the second detection point set data;
and sorting the second detection point set data subjected to tau-p inverse transformation to obtain the new shot set data.
8. The shear wave static correction system of claim 5, wherein the model track data module is specifically configured to:
performing superposition processing on the p set data to obtain model channel data corresponding to the p set data; or the like, or, alternatively,
and smoothing the p set data to obtain model channel data corresponding to the p set data.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the following steps when executing the computer program:
obtaining p set data according to the original shot set data;
obtaining model channel data corresponding to the p set data according to the p set data;
selecting a plurality of first time windows according to the transverse wave data in the p set data;
acquiring p set data and model channel data in each first time window;
performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift amount of each first time window;
obtaining a second time shift amount according to the first time shift amount of each first time window;
correcting the p set data according to the second time shift amount to obtain corrected p set data;
obtaining new shot gather data according to the corrected p set data;
selecting a plurality of second time windows according to converted wave data in the new shot gather data;
acquiring original shot gather data and new shot gather data in each second time window;
performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
obtaining a fourth time shift amount according to the third time shift amount of each second time window;
and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
10. A computer-readable storage medium, on which a computer program is stored, which computer program, when being executed by a processor, carries out the steps of:
obtaining p set data according to the original shot set data;
obtaining model channel data corresponding to the p set data according to the p set data;
selecting a plurality of first time windows according to the transverse wave data in the p set data;
acquiring p set data and model channel data in each first time window;
performing cross correlation on the p set data and the model channel data in each first time window to obtain a first time shift amount of each first time window;
obtaining a second time shift amount according to the first time shift amount of each first time window;
correcting the p set data according to the second time shift amount to obtain corrected p set data;
obtaining new shot gather data according to the corrected p set data;
selecting a plurality of second time windows according to converted wave data in the new shot gather data;
acquiring original shot gather data and new shot gather data in each second time window;
performing cross-correlation on the original shot gather data and the new shot gather data in each second time window to obtain a third time shift amount of each second time window;
obtaining a fourth time shift amount according to the third time shift amount of each second time window;
and correcting the original shot gather data according to the fourth time shift amount to obtain corrected shot gather data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811567259.9A CN111352157B (en) | 2018-12-20 | 2018-12-20 | Shear wave static correction method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811567259.9A CN111352157B (en) | 2018-12-20 | 2018-12-20 | Shear wave static correction method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111352157A CN111352157A (en) | 2020-06-30 |
| CN111352157B true CN111352157B (en) | 2022-08-30 |
Family
ID=71196765
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811567259.9A Active CN111352157B (en) | 2018-12-20 | 2018-12-20 | Shear wave static correction method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111352157B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112558157B (en) * | 2020-12-08 | 2021-10-15 | 中国地质大学（北京） | Static correction method and device based on multi-wave combination |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1101797A (en) * | 1965-07-06 | 1968-01-31 | Mobil Oil Corp | Seismic corrections |
| GB1380886A (en) * | 1973-04-13 | 1975-01-15 | Chevron Res | Static corrections for seismic traces by cross-correlation method |
| GB8410296D0 (en) * | 1983-04-20 | 1984-05-31 | Chevron Res | Seismic processing |
| WO2000073820A1 (en) * | 1999-05-27 | 2000-12-07 | Phillips Petroleum Company | Seismic velocity analysis for class ii sands |
| CN101072084A (en) * | 2006-05-11 | 2007-11-14 | 中国石油化工股份有限公司 | Wideband digital signal orthogonal beam frequency division method and its use |
| CN103728665A (en) * | 2012-10-12 | 2014-04-16 | 中国石油化工股份有限公司 | Fine residual static correction method |
| CN107193047A (en) * | 2016-03-15 | 2017-09-22 | 中国石油化工股份有限公司 | Reflection static corr ections time-varying dynamic residual static corrections and device |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5128899A (en) * | 1991-07-16 | 1992-07-07 | Conoco Inc. | Method for sharpening prestack depth migration images by removing residual moveout from common image point gathers before stacking |
| CN101201409B (en) * | 2006-12-14 | 2011-05-25 | 中国石油天然气集团公司 | Method for revising earthquake data phase |
| CN101907727B (en) * | 2010-08-17 | 2012-05-30 | 中国科学院地质与地球物理研究所 | Multi-component converted wave static correction method by using surface waves |
| CN102565857B (en) * | 2011-12-16 | 2014-09-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Automatic remaining dynamic correction method |
| US9291736B2 (en) * | 2012-01-05 | 2016-03-22 | Cggveritas Services Sa | Surface-consistent amplitude and deconvolution simultaneous joined inversion |
| CN104656139B (en) * | 2013-11-21 | 2018-04-06 | 中国石油天然气集团公司 | The static correcting method and equipment of a kind of converted wave |
| CN103630937B (en) * | 2013-12-09 | 2016-08-03 | 中国石油天然气集团公司 | The method and apparatus that earthquake record is carried out shear wave static correction process |
| CN104914469B (en) * | 2014-03-13 | 2017-08-04 | 中国石油化工股份有限公司 | Static correcting method and device in a kind of converted shear wave |
| US11307320B2 (en) * | 2015-01-15 | 2022-04-19 | Jyoti Behura | Expedient processing and waveform inversion of seismic data |
| CN108919345B (en) * | 2018-05-16 | 2020-07-07 | 中国海洋石油集团有限公司 | Submarine cable land detection noise attenuation method |
-
2018
- 2018-12-20 CN CN201811567259.9A patent/CN111352157B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1101797A (en) * | 1965-07-06 | 1968-01-31 | Mobil Oil Corp | Seismic corrections |
| GB1380886A (en) * | 1973-04-13 | 1975-01-15 | Chevron Res | Static corrections for seismic traces by cross-correlation method |
| GB8410296D0 (en) * | 1983-04-20 | 1984-05-31 | Chevron Res | Seismic processing |
| WO2000073820A1 (en) * | 1999-05-27 | 2000-12-07 | Phillips Petroleum Company | Seismic velocity analysis for class ii sands |
| CN101072084A (en) * | 2006-05-11 | 2007-11-14 | 中国石油化工股份有限公司 | Wideband digital signal orthogonal beam frequency division method and its use |
| CN103728665A (en) * | 2012-10-12 | 2014-04-16 | 中国石油化工股份有限公司 | Fine residual static correction method |
| CN107193047A (en) * | 2016-03-15 | 2017-09-22 | 中国石油化工股份有限公司 | Reflection static corr ections time-varying dynamic residual static corrections and device |
Non-Patent Citations (3)
| Title |
|---|
| Computing near-surface velocity models for S-wave static corrections using raypath interferometry;Cova.R等;《GEOPHYSICS》;20180913;第U23-24页 * |
| 几种转换波静校正方法讨论;张旭东;《岩性油气藏》;20101215;第22卷(第4期);第95-99页 * |
| 基于低频模型道的剩余静校正方法在长垣油田地震资料处理中的应用;赵忠华;《油气藏评价与开发》;20140226;第4卷(第1期);第16-20页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111352157A (en) | 2020-06-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (en) | Dolomite oil and gas reservoir distribution prediction method and device | |
| CN102937721B (en) | Limited frequency tomography method for utilizing preliminary wave travel time | |
| CN106405651B (en) | Full waveform inversion initial velocity model construction method based on logging matching | |
| CN109669212B (en) | Seismic data processing method, stratum quality factor estimation method and device | |
| CN110018517B (en) | Multi-scale ground micro-seismic reverse-time interference positioning method | |
| CN107065013B (en) | A kind of interval velocity under earthquake scale determines method and device | |
| CN113625337B (en) | Ultra-shallow water high-precision seismic data rapid imaging method | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| CN103713324B (en) | Self-adaption curvelet threshold value earthquake denoising method based on local variance analysis | |
| CN111352157B (en) | Shear wave static correction method and system | |
| CN108562936B (en) | Crack prediction method, system, storage medium and terminal | |
| KR101693310B1 (en) | Full-waveform inversion using plane wave method in anisotropic media | |
| CN109143345B (en) | Quality factor Q nonlinear inversion method and system based on simulated annealing | |
| CN114428343A (en) | Markenko imaging method and system based on normalized cross-correlation | |
| CN111123361B (en) | Method and device for regularly reconstructing seismic data of vertical seismic profile and storage medium | |
| CN107918152B (en) | A kind of seismic coherence chromatography imaging method | |
| Mostafanejad et al. | Velocity structure of the northern Mississippi Embayment sediments, Part I: Teleseismic P‐wave spectral ratios analysis | |
| US20140249755A1 (en) | Method and device for calculating time-shifts and time-strains in seismic data | |
| CN115079257A (en) | Q value estimation and seismic attenuation compensation method based on fusion network | |
| CN110824555B (en) | Seismic energy equalization method and device and computer readable storage medium | |
| Chen et al. | Least-squares decomposition with time–space constraint for denoising microseismic data | |
| CN112415601A (en) | Method and device for determining surface quality factor Q value | |
| Huang et al. | Shannon entropy-based seismic local correlation measure and enhancement | |
| CN113433588B (en) | Near-surface velocity analysis method based on offset scanning superposition | |
| CN110837119A (en) | Seismic data processing method and system for enhancing inverse Q value compensation stability |
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 |